Migration velocity analysis using the common image cube (CIC)

Geophysics ◽  
2011 ◽  
Vol 76 (5) ◽  
pp. WB127-WB134 ◽  
Author(s):  
Saleh M. Al-Saleh ◽  
Jianwu Jiao ◽  
Adam J. Fox

Migration velocity analysis (MVA) is commonly performed in the image domain in conjunction with ray-based tomography to update the velocity model. This approach can be challenging in the presence of large velocity errors as it may require many MVA iterations before converging to a model that can focus the events in the image domain. We introduced a downward continuation-based domain for carrying out MVA that is more flexible than conventional domains. This approach consists of two steps: (1) forming the common image cube (CIC) and (2) modeling the Green’s functions. In the first step, the cross-correlation imaging condition is relaxed to produce more than the zero lag common image gather (CIG). Slicing these data at different lags forms a series of CIGs, whereas a conventional CIG can be obtained by slicing the cube at the zero lag. When the velocity model used for the migration differs from the true velocity model, properly flattened events may occur in CIGs other than the zero lag. In the second step, for each event on the CIG, we picked the cross-correlation lag and depth at which it flattens best. For each event, we modeled a Green’s function by seeding a source at the focusing depth using one-way wave-equation modeling. This process is then repeated for other events at different lateral positions. The result is a set of Green’s functions whose wavefield approximates the ones that would have been generated if the correct velocity model was used to simulate these gathers. The updated Green functions are easier to work with than the raw data as they have less noise. Wavefield tomography can then be applied on these data-driven, modeled Green’s functions to build the final velocity model. Tests on synthetic and real 2D data confirm the method’s effectiveness in building velocity models in complex structural areas with large lateral velocity variations.

Geophysics ◽  
2008 ◽  
Vol 73 (5) ◽  
pp. VE161-VE171 ◽  
Author(s):  
J. Schleicher ◽  
J. C. Costa ◽  
A. Novais

Image-wave propagation or velocity continuation describes the variation of the migrated position of a seismic event as a function of migration velocity. Image-wave propagation in the common-image gather (CIG) domain can be combined with residual-moveout analysis for iterative migration velocity analysis (MVA). Velocity continuation of CIGs leads to a detection of those velocities in which events flatten. Although image-wave continuation is based on the assumption of a constant migration velocity, the procedure can be applied in inhomogeneous media. For this purpose, CIGs obtained by migration with an inhomogeneous macrovelocity model are continued starting from a constant reference velocity. The interpretation of continued CIGs, as if they were obtained from residual migrations, leads to a correction formula that translates residual flattening velocities into absolute time-migration velocities. In this way, the migration velocity model can be improved iteratively until a satisfactory result is reached. With a numerical example, we found that MVA with iterative image continuation applied exclusively to selected CIGs can construct a reasonable migration velocity model from scratch, without the need to build an initial model from a previous conventional normal-moveout/dip-moveout velocity analysis.


Geophysics ◽  
2018 ◽  
Vol 83 (5) ◽  
pp. R475-R495
Author(s):  
Emmanuel Cocher ◽  
Hervé Chauris ◽  
René-Édouard Plessix

Migration velocity analysis is a family of methods aiming at automatically recovering large-scale trends of the velocity model from primary reflection data. We studied an image-domain version, in which the model is extended with the subsurface offset and we use the differential semblance optimization objective function. To incorporate first-order surface multiples in this method, the standard migration step is replaced by a least-squares iterative scheme aiming at determining an extended reflectivity model explaining primaries and multiples. Hence, this iterative migration velocity analysis strategy takes the form of a nested optimization problem, with gradient-based minimization techniques for the inner (migration part) and outer loops (macromodel estimation). The behavior of the outer loop gradient is unstable, depending on the number of iterations of the inner loop. This problem is addressed by slightly modifying the outer loop objective function: A “filter” operator attenuating unwanted energy in the extended reflectivity is applied before evaluating the focusing of reflectivity images. Simple synthetic numerical examples illustrate that this modification improves the stability of the gradient. In addition, a less expensive outer gradient computation is proposed, without harming the background velocity updates.


Geophysics ◽  
2002 ◽  
Vol 67 (4) ◽  
pp. 1202-1212 ◽  
Author(s):  
Hervé Chauris ◽  
Mark S. Noble ◽  
Gilles Lambaré ◽  
Pascal Podvin

We present a new method based on migration velocity analysis (MVA) to estimate 2‐D velocity models from seismic reflection data with no assumption on reflector geometry or the background velocity field. Classical approaches using picking on common image gathers (CIGs) must consider continuous events over the whole panel. This interpretive step may be difficult—particularly for applications on real data sets. We propose to overcome the limiting factor by considering locally coherent events. A locally coherent event can be defined whenever the imaged reflectivity locally shows lateral coherency at some location in the image cube. In the prestack depth‐migrated volume obtained for an a priori velocity model, locally coherent events are picked automatically, without interpretation, and are characterized by their positions and slopes (tangent to the event). Even a single locally coherent event has information on the unknown velocity model, carried by the value of the slope measured in the CIG. The velocity is estimated by minimizing these slopes. We first introduce the cost function and explain its physical meaning. The theoretical developments lead to two equivalent expressions of the cost function: one formulated in the depth‐migrated domain on locally coherent events in CIGs and the other in the time domain. We thus establish direct links between different methods devoted to velocity estimation: migration velocity analysis using locally coherent events and slope tomography. We finally explain how to compute the gradient of the cost function using paraxial ray tracing to update the velocity model. Our method provides smooth, inverted velocity models consistent with Kirchhoff‐type migration schemes and requires neither the introduction of interfaces nor the interpretation of continuous events. As for most automatic velocity analysis methods, careful preprocessing must be applied to remove coherent noise such as multiples.


Geophysics ◽  
2021 ◽  
pp. 1-50
Author(s):  
German Garabito ◽  
José Silas dos Santos Silva ◽  
Williams Lima

In land seismic data processing, the prestack time migration (PSTM) image remains the standard imaging output, but a reliable migrated image of the subsurface depends on the accuracy of the migration velocity model. We have adopted two new algorithms for time-domain migration velocity analysis based on wavefield attributes of the common-reflection-surface (CRS) stack method. These attributes, extracted from multicoverage data, were successfully applied to build the velocity model in the depth domain through tomographic inversion of the normal-incidence-point (NIP) wave. However, there is no practical and reliable method for determining an accurate and geologically consistent time-migration velocity model from these CRS attributes. We introduce an interactive method to determine the migration velocity model in the time domain based on the application of NIP wave attributes and the CRS stacking operator for diffractions, to generate synthetic diffractions on the reflection events of the zero-offset (ZO) CRS stacked section. In the ZO data with diffractions, the poststack time migration (post-STM) is applied with a set of constant velocities, and the migration velocities are then selected through a focusing analysis of the simulated diffractions. We also introduce an algorithm to automatically calculate the migration velocity model from the CRS attributes picked for the main reflection events in the ZO data. We determine the precision of our diffraction focusing velocity analysis and the automatic velocity calculation algorithms using two synthetic models. We also applied them to real 2D land data with low quality and low fold to estimate the time-domain migration velocity model. The velocity models obtained through our methods were validated by applying them in the Kirchhoff PSTM of real data, in which the velocity model from the diffraction focusing analysis provided significant improvements in the quality of the migrated image compared to the legacy image and to the migrated image obtained using the automatically calculated velocity model.


Geophysics ◽  
2019 ◽  
Vol 84 (5) ◽  
pp. S437-S447 ◽  
Author(s):  
Jean-Philippe Montel ◽  
Gilles Lambaré

Common-image gathers are a useful output of the migration process. Their kinematic behavior (i.e., the way they curve up or down) is an indicator of the quality of the velocity model used for migration. Traditionally, when used for migration velocity analysis, we pick structural dips in the common attribute panels (offset, angle, etc.) and residual moveout (RMO) in the gathers. The measured RMO will then tell us how much we need to update the velocity model to improve the gather’s flatness. Understanding the kinematics of the picked events is the key to an accurate model update. This point has been widely underestimated in many cases. For example, when dealing with angle gathers, there is a general assumption that the associated tomographic rays are fully defined by the picked structural dips and the gather opening and azimuth angle, and that if the velocity model is correctly updated down to a given horizon, it is not necessary to shoot the tomographic rays upward through this horizon. We find through an original theoretical analysis that both of these assumptions have to be modified when the gathers exhibit RMO. Using a kinematic analysis, we determine that knowledge of the RMO slopes is necessary to compute the tomographic rays.


Geophysics ◽  
1998 ◽  
Vol 63 (4) ◽  
pp. 1200-1209 ◽  
Author(s):  
Jinming Zhu ◽  
Larry Lines ◽  
Sam Gray

Reliable seismic depth migrations require an accurate input velocity model. Inaccurate velocity estimates will distort point diffractors into smiles or frowns on a depth section. For both poststack and prestack migrated sections, high velocities cause deep smiles while low velocities cause shallow frowns on migrated gathers. However, for prestack images in the offset domain, high velocities cause deep frowns while low velocities cause shallow smiles. If the velocity is correct, there will be no variation in the depth migration as a function of offset and no smiles or frowns in the offset domain. We explain migration responses both mathematically and graphically and thereby provide the basis for depth migration velocity analysis.


2017 ◽  
Vol 21 (4) ◽  
pp. 759-780 ◽  
Author(s):  
Emmanuel Cocher ◽  
Hervé Chauris ◽  
René-Édouard Plessix

Sign in / Sign up

Export Citation Format

Share Document