Structure constrained least-squares migration of total reflection and its applications

2020 ◽  
Author(s):  
Cheng Cheng ◽  
Yang He ◽  
Jian Mao ◽  
Bin Wang
Geophysics ◽  
2000 ◽  
Vol 65 (4) ◽  
pp. 1195-1209 ◽  
Author(s):  
Bertrand Duquet ◽  
Kurt J. Marfurt ◽  
Joe A. Dellinger

Because of its computational efficiency, prestack Kirchhoff depth migration is currently one of the most popular algorithms used in 2-D and 3-D subsurface depth imaging. Nevertheless, Kirchhoff algorithms in their typical implementation produce less than ideal results in complex terranes where multipathing from the surface to a given image point may occur, and beneath fast carbonates, salt, or volcanics through which ray‐theoretical energy cannot penetrate to illuminate underlying slower‐velocity sediments. To evaluate the likely effectiveness of a proposed seismic‐acquisition program, we could perform a forward‐modeling study, but this can be expensive. We show how Kirchhoff modeling can be defined as the mathematical transpose of Kirchhoff migration. The resulting Kirchhoff modeling algorithm has the same low computational cost as Kirchhoff migration and, unlike expensive full acoustic or elastic wave‐equation methods, only models the events that Kirchhoff migration can image. Kirchhoff modeling is also a necessary element of constrained least‐squares Kirchhoff migration. We show how including a simple a priori constraint during the inversion (that adjacent common‐offset images should be similar) can greatly improve the resulting image by partially compensating for irregularities in surface sampling (including missing data), as well as for irregularities in ray coverage due to strong lateral variations in velocity and our failure to account for multipathing. By allowing unstacked common‐offset gathers to become interpretable, the additional cost of constrained least‐squares migration may be justifiable for velocity analysis and amplitude‐variation‐with‐offset studies. One useful by‐product of least‐squares migration is an image of the subsurface illumination for each offset. If the data are sufficiently well sampled (so that including the constraint term is not necessary), the illumination can instead be calculated directly and used to balance the result of conventional migration, obtaining most of the advantages of least‐squares migration for only about twice the cost of conventional migration.


Geophysics ◽  
2003 ◽  
Vol 68 (1) ◽  
pp. 262-273 ◽  
Author(s):  
Henning Kühl ◽  
Mauricio D. Sacchi

We present an acoustic migration/inversion algorithm that uses extended double‐square‐root wave‐equation migration and modeling operators to minimize a constrained least‐squares data misfit function (least‐squares migration). We employ an imaging principle that allows for the extraction of ray‐parameter‐domain common image gathers (CIGs) from the propagated seismic wavefield. The CIGs exhibit amplitude variations as a function of half‐offset ray parameter (AVP) closely related to the amplitude variation with reflection angle (AVA). Our least‐squares wave‐equation migration/inversion is constrained by a smoothing regularization along the ray parameter. This approach is based on the idea that rapid amplitude changes or discontinuities along the ray parameter axis result from noise, incomplete wavefield sampling, and numerical operator artifacts. These discontinuities should therefore be penalized in the inversion. The performance of the proposed algorithm is examined with two synthetic examples. In the first case, we generated acoustic finite difference data for a horizontally layered model. The AVP functions based on the migrated/inverted ray parameter CIGs were converted to AVA plots. The AVA plots were then compared to the true acoustic AVA of the reflectors. The constrained least‐squares inversion compares favorably with the conventional migration, especially when incompleteness compromises the data. In the second example, we use the Marmousi data set to test the algorithm in complex media. The result shows that least‐squares migration can mitigate kinematic artifacts in the ray‐parameter domain CIGs effectively.


2019 ◽  
Author(s):  
Bruno Dias ◽  
Cláudio Guerra ◽  
André Bulcão ◽  
Roberto Dias

2020 ◽  
Author(s):  
Lian Duan ◽  
Alejandro Valenciano ◽  
Nizar Chemingui

1978 ◽  
Vol 3 (1) ◽  
pp. 43-47 ◽  
Author(s):  
S.A. Coons

Sign in / Sign up

Export Citation Format

Share Document