Kirchhoff prestack depth migration: Imaging conditions and amplitude recovery

Author(s):  
David E. Lumley
2018 ◽  
Author(s):  
Jianguo Li ◽  
Jianhua Huang ◽  
Yanpeng Li ◽  
Yuanzhong Chenand Junjun Wu

Geophysics ◽  
1992 ◽  
Vol 57 (12) ◽  
pp. 1608-1622 ◽  
Author(s):  
Scott MacKay ◽  
Ray Abma

Prestack depth migration uses two imaging conditions, zero time and zero offset, during downward continuation to form a migrated depth section. When the migration velocities are exact, the two imaging conditions act in a complementary fashion to yield a focused image. When the migration velocities are in error, reflected energy collapses to zero offset at depths that are inconsistent with the zero‐time imaging condition. The result is a deteriorated seismic image. However, by interpreting the nonzero times at which focusing actually occurs, the migration velocities can be updated iteratively in a process called depth‐focusing analysis. To produce a well‐focused seismic image, the goal of depth‐focusing analysis must be the elimination of focusing errors; however, practical considerations can prevent this goal from being achieved. Therefore, to relax the sensitivity of the migrated image to focusing errors, we introduce a nonzero‐time imaging condition by extracting the data along the interpreted surface of focusing from the depth‐focusing analysis volume. This method, called focal‐surface imaging, estimates the results of prestack depth migration using the updated velocities. Depth‐focusing analysis is shown to be a robust approach to velocity estimation and imaging. Limitations arising from constant‐velocity and low‐dip approximations are reduced in the presence of increasing velocities with depth. Lateral velocity errors, sources of exaggerated focusing errors and diverging velocity solutions, can also be addressed by applying a damping factor to the interpreted depth errors. Velocity estimation and focal‐surface imgaging, using iterative prestack depth migration, were applied to a southern North Sea data set. Starting with a regional velocity function, the first iteration provided an updated velocity field that more accurately conformed to the known lithologies. The focal‐surface image, formed from the same iteration, contained significantly more focused energy than the conventional section formed by prestack depth migration. However, structural differences between the two sections indicated the need for another iteration of migration using the updated velocities. The second iteration indicated smaller velocity errors and enough similarity between the migrated section and the new focal‐surface image to indicate that further iterations were unnecessary.


Geophysics ◽  
1997 ◽  
Vol 62 (4) ◽  
pp. 1226-1237 ◽  
Author(s):  
Irina Apostoiu‐Marin ◽  
Andreas Ehinger

Prestack depth migration can be used in the velocity model estimation process if one succeeds in interpreting depth events obtained with erroneous velocity models. The interpretational difficulty arises from the fact that migration with erroneous velocity does not yield the geologically correct reflector geometries and that individual migrated images suffer from poor signal‐to‐noise ratio. Moreover, migrated events may be of considerable complexity and thus hard to identify. In this paper, we examine the influence of wrong velocity models on the output of prestack depth migration in the case of straight reflector and point diffractor data in homogeneous media. To avoid obscuring migration results by artifacts (“smiles”), we use a geometrical technique for modeling and migration yielding a point‐to‐point map from time‐domain data to depth‐domain data. We discover that strong deformation of migrated events may occur even in situations of simple structures and small velocity errors. From a kinematical point of view, we compare the results of common‐shot and common‐offset migration. and we find that common‐offset migration with erroneous velocity models yields less severe image distortion than common‐shot migration. However, for any kind of migration, it is important to use the entire cube of migrated data to consistently interpret in the prestack depth‐migrated domain.


2021 ◽  
Author(s):  
Olaf Hellwig ◽  
Stefan Buske

<p>The polymetallic, hydrothermal deposit of the Freiberg mining district in the southeastern part of Germany is characterised by ore veins that are framed by Proterozoic orthogneiss. The ore veins consist mainly of quarz, sulfides, carbonates, barite and flourite, which are associated with silver, lead and tin. Today the Freiberg University of Mining and Technology is operating the shafts Reiche Zeche and Alte Elisabeth for research and teaching purposes with altogether 14 km of accessible underground galleries. The mine together with the most prominent geological structures of the central mining district are included in a 3D digital model, which is used in this study to study seismic acquisition geometries that can help to image the shallow as well as the deeper parts of the ore-bearing veins. These veins with dip angles between 40° and 85° are represented by triangulated surfaces in the digital geological model. In order to import these surfaces into our seismic finite-difference simulation code, they have to be converted into bodies with a certain thickness and specific elastic properties in a first step. In a second step, these bodies with their properties have to be discretized on a hexahedral finite-difference grid with dimensions of 1000 m by 1000 m in the horizontal direction and 500 m in the vertical direction. Sources and receiver lines are placed on the surface along roads near the mine. A Ricker wavelet with a central frequency of 50 Hz is used as the source signature at all excitation points. Beside the surface receivers, additional receivers are situated in accessible galleries of the mine at three different depth levels of 100 m, 150 m and 220 m below the surface. Since previous mining activities followed primarily the ore veins, there are only few pilot-headings that cut through longer gneiss sections. Only these positions surrounded by gneiss are suitable for imaging the ore veins. Based on this geometry, a synthetic seismic data set is generated with our explicit finite-difference time-stepping scheme, which solves the acoustic wave equation with second order accurate finite-difference operators in space and time. The scheme is parallelised using a decomposition of the spatial finite-difference grid into subdomains and Message Passing Interface for the exchange of the wavefields between neighbouring subdomains. The resulting synthetic seismic shot gathers are used as input for Kirchhoff prestack depth migration as well as Fresnel volume migration in order to image the ore veins. Only a top mute to remove the direct waves and a time-dependent gain to correct the amplitude decay due to the geometrical spreading are applied to the data before the migration. The combination of surface and in-mine acquisition helps to improve the image of the deeper parts of the dipping ore veins. Considering the limitations for placing receivers in the mine, Fresnel volume migration as a focusing version of Kirchhoff prestack depth migration helps to avoid migration artefacts caused by this sparse and limited acquisition geometry.</p>


2006 ◽  
Author(s):  
Scott MacKay ◽  
Héctor Ramírez Jiménez ◽  
Jorge San Martín Romero ◽  
Mark Morford

Geophysics ◽  
2009 ◽  
Vol 74 (4) ◽  
pp. S67-S74 ◽  
Author(s):  
Jun Cao ◽  
Ru-Shan Wu

Wave-equation-based acquisition aperture correction in the local angle domain can improve image amplitude significantly in prestack depth migration. However, its original implementation is inefficient because the wavefield decomposition uses the local slant stack (LSS), which is demanding computationally. We propose a faster method to obtain the image and amplitude correction factor in the local angle domain using beamlet decomposition in the local wavenumber domain. For a given frequency, the image matrix in the local wavenumber domain for all shots can be calculated efficiently. We then transform the shot-summed image matrix from the local wavenumber domain to the local angle domain (LAD). The LAD amplitude correction factor can be obtained with a similar strategy. Having a calculated image and correction factor, one can apply similar acquisition aperture corrections to the original LSS-based method. For the new implementation, we compare the accuracy and efficiency of two beamlet decompositions: Gabor-Daubechies frame (GDF) and local exponential frame (LEF). With both decompositions, our method produces results similar to the original LSS-based method. However, our method can be more than twice as fast as LSS and cost only twice the computation time of traditional one-way wave-equation-based migrations. The results from GDF decomposition are superior to those from LEF decomposition in terms of artifacts, although GDF requires a little more computing time.


Sign in / Sign up

Export Citation Format

Share Document