3‐D prestack depth migration on real data

1990 ◽  
Author(s):  
Philippe Julien ◽  
Steve Cole ◽  
Stefan Kaculini
Keyword(s):  
Real Data ◽  
Prestack Depth Migration ◽  
Depth Migration

Slowness-driven Gaussian-beam prestack depth migration for low-fold seismic data

Geophysics ◽  
2009 ◽  
Vol 74 (6) ◽  
pp. WCA35-WCA45 ◽  
Author(s):  
Chaoshun Hu ◽  
Paul L. Stoffa
Keyword(s):  
Gaussian Beam ◽  
Seismic Data ◽  
Fresnel Zone ◽  
Seismic Energy ◽  
Real Data ◽  
Ocean Bottom ◽  
Prestack Depth Migration ◽  
Seismic Reflection Data ◽  
Depth Migration ◽  
Reflection Data

Subsurface images based on low-fold seismic reflection data or data with geometry acquisition limitations, such as obtained from ocean-bottom seismography (OBS), are often corrupted by migration swing artifacts. Incorporating prestack instantaneous slowness information into the imaging condition can significantly reduce these migration swing artifacts and improve image quality, especially for areas with poor illumination. We combine the horizontal surface slowness information of observed seismic data with Gaussian-beam depth migration to implement a new slowness-driven Gaussian-beam prestack depth migration whereby Fresnel weighting is combined naturally with beam summation. The prestack instantaneous slowness information is extracted from the original OBS or shot gathers using local slant stacks and is combined with a local semblance analysis. During migration, we propagate the seismic energy downward, knowing its instantaneous slowness information. At each image location, the beam summation is localized in a resolution-dependent Fresnel zone; the instantaneous slowness information controls the beam summation. Synthetic and real data examples confirm that slowness-driven Gaussian-beam migration can suppress most noise from inadequate stacking and give a clearer migration result.

VTI anisotropic corrections and effective parameter estimation after isotropic prestack depth migration

Geophysics ◽  
2006 ◽  
Vol 71 (3) ◽  
pp. D35-D43 ◽  
Author(s):  
Moshe Reshef ◽  
Murray Roth
Keyword(s):  
Parameter Estimation ◽  
Vertical Velocity ◽  
Transverse Isotropy ◽  
Real Data ◽  
P Wave ◽  
Prestack Depth Migration ◽  
Anisotropic Parameter ◽  
Depth Migration ◽  
Residual Correction ◽  
Dip Angle

In the method for applying anisotropic corrections after isotropic prestack depth migration (PSDM), the correction, which is calculated and implemented in the depth domain, is defined as a time difference between isotropic and anisotropic traveltimes, under the assumption that the vertical velocity is known. The definition of this correction uses a special postmigration common-image-gather (CIG) ordering, which collects the migrated data according to the input-trace’s source and receiver distance from the surface CIG location. In this postmigration domain, the dip of the events can be directly related to their horizontal position in the CIG, called the imaging offset, and the separation of flat and dipping reflectors becomes easy to perform. The dependency of the seismic anisotropic effect on the subsurface dip angle is well pronounced in these CIGs. After application of an isotropic PSDM, effective anisotropic-parameter estimation is performed at selected CIG locations by using a simple two-parameter scan procedure. The optimal anisotropic parameters can be used to perform a final anisotropic PSDM or to apply a residual correction to the isotropically migrated data. We demonstrate the method for P-wave data in 2D media with vertical transverse isotropy (VTI) symmetry by using both synthetic and real data. We also present a strategy for handling the ambiguity between the vertical velocity and the anisotropic parameters.

Curved-wave prestack depth migration

Geophysics ◽  
2011 ◽  
Vol 76 (3) ◽  
pp. S115-S129
Author(s):  
Xiumei Chen ◽  
Zaitian Ma ◽  
Huazhong Wang ◽  
Guohua Nie
Keyword(s):  
Synthesis Method ◽  
Real Data ◽  
Target Area ◽  
High Quality ◽  
Prestack Depth Migration ◽  
Kernel Operator ◽  
Depth Migration ◽  
High Quality Image ◽  
Quality Image ◽  
Source Excitation

Conventional prestack depth migration (PSDM) based on point-source shot data involves a huge number of time-consuming wave extrapolations. Areal-source prestack data generated by a synthesis procedure can provide an alternative migration scheme that is efficient and accurate. We have developed a curved-wave synthesis method to implement curved-wave PSDM technology. A base-kernel synthesis operator is defined as a relative curved source-excitation datum. A general synthesis operator is constructed by disturbing the illumination directions of a base-kernel operator. The curved wavefields are synthesized by synthesis operators and then migrated by combination with a common-shot depth-migration algorithm. Curved wavefields can be constructed according to the configuration of a geologic structure to obtain a high-quality image of the target. Curved-wave records arising from different curved-source wavefields with different illumination properties are stacked after depth migration to produce an accurate image. A synthesis operator can be designed by simulating the configuration of target structures, or by velocity constraint in a target-oriented way, to efficiently achieve a high-quality image of a complex target area. Characteristic of curved excitation surface and space-variant illumination, the method is more adaptable to a complex medium than the plane-wave technique and more efficient than shot migration. Numerical demonstrations on synthetic and real data have proved effective and given good results.

Separation and imaging of seismic diffractions using migrated dip-angle gathers

Geophysics ◽  
2012 ◽  
Vol 77 (6) ◽  
pp. S131-S143 ◽  
Author(s):  
Alexander Klokov ◽  
Sergey Fomel
Keyword(s):  
Radon Transform ◽  
Real Data ◽  
Separation Procedure ◽  
Shape Difference ◽  
Depth Migration ◽  
Time Migration ◽  
Reflection Angle ◽  
Radon Space ◽  
Dip Angle ◽  
Analytical Expressions

Common-reflection angle migration can produce migrated gathers either in the scattering-angle domain or in the dip-angle domain. The latter reveals a clear distinction between reflection and diffraction events. We derived analytical expressions for events in the dip-angle domain and found that the shape difference can be used for reflection/diffraction separation. We defined reflection and diffraction models in the Radon space. The Radon transform allowed us to isolate diffractions from reflections and noise. The separation procedure can be performed after either time migration or depth migration. Synthetic and real data examples confirmed the validity of this technique.

Kinematic interpretation in the prestack depth‐migrated domain

Geophysics ◽  
1997 ◽  
Vol 62 (4) ◽  
pp. 1226-1237 ◽  
Author(s):  
Irina Apostoiu‐Marin ◽  
Andreas Ehinger
Keyword(s):  
Signal To Noise Ratio ◽  
Velocity Model ◽  
Point Of View ◽  
Prestack Depth Migration ◽  
Depth Migration ◽  
Velocity Models ◽  
Time Domain Data ◽  
And Migration ◽  
Point To Point ◽  
Homogeneous Media

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.

Simulation of a seismic survey with combined surface and in-mine data acquisition for imaging steeply dipping ore veins

2021 ◽  
Author(s):  
Olaf Hellwig ◽  
Stefan Buske
Keyword(s):  
Finite Difference ◽  
Message Passing Interface ◽  
Seismic Survey ◽  
Difference Operators ◽  
Mining District ◽  
Prestack Depth Migration ◽  
Data Set ◽  
Depth Migration ◽  
Difference Grid ◽  
Triangulated Surfaces

<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>

Calibrating prestack depth migration volumes with well control

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

Fast acquisition aperture correction in prestack depth migration using beamlet decomposition

Geophysics ◽  
2009 ◽  
Vol 74 (4) ◽  
pp. S67-S74 ◽  
Author(s):  
Jun Cao ◽  
Ru-Shan Wu
Keyword(s):  
Wave Equation ◽  
Correction Factor ◽  
Computing Time ◽  
Computation Time ◽  
Prestack Depth Migration ◽  
Depth Migration ◽  
Wavenumber Domain ◽  
Amplitude Correction ◽  
Local Wavenumber ◽  
Local Angle

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