Redatuming and Quantifying Attenuation from Reflection Data Using the Marchenko Equation

Abstract

Marchenko Imaging is a new technology in geophysics which enables to retrieve Green's functions at any point in the subsurface having only reflection data. This method is based on the extension of the 1D Gelfand-Levitan-Marchenko equation to a 3D medium. One of the assumptions of the Marchenko method is that the medium is lossless. If the lossy reflection response is used in the Marchenko scheme, some artefacts in the Green's functions as well as in the seismic image are present. One way to circumvent this assumption is to find a compensation parameter for the lossy reflection series so that the lossless Marchenko scheme can be applied. The main tasks of this thesis are to: [1] use the Marchenko equation to estimate the attenuation in the subsurface, [2] find a compensation parameter for the lossy reflection series so that the lossless Marchenko scheme can be applied, and [3] to create an upscaling method for wave propagation. The Artefact Removal Method was created which makes it possible to calculate an effective temporal Q-factor of the medium between a virtual source in the subsurface and receivers at the surface. This method is based on the minimization of the artefacts produced by the lossless Marchenko scheme. The minimization was performed in three ways: [1] in the space-time domain, [2] in the frequency domain and [3] to the scales of the wavelet transform applied to the artefacts. This method can also be used to find the layers with high attenuation. The upscaling method which can be used to construct macro-scale homogenized viscoelastic properties of the medium from the micro-scale properties of a heterogeneous medium was developed. This is done through linking the macro- and micro- scale Lippmann-Schwinger equations which describe the wave field and the strain field scattering in an inhomogeneous medium, respectively. In this thesis, the macro-scale homogenized viscoelastic properties were calculated by using the T-matrix Approach and the Generalized Dvorkin-Mavko Attenuation Model. All theoretical results are supported by synthetic 1D modeling. The theoretical part of the thesis and the general work flow can be used for a very complex medium.