Linearized wave-equation migration velocity analysis by image warping

Geophysics ◽  
2014 ◽  
Vol 79 (2) ◽  
pp. S35-S46 ◽  
Author(s):  
Francesco Perrone ◽  
Paul Sava ◽  
Clara Andreoletti ◽  
Nicola Bienati
Keyword(s):  
Wave Equation ◽  
Velocity Model ◽  
Migration Velocity ◽  
Background Model ◽  
Data Migration ◽  
Physical Parameters ◽  
Model Parameters ◽  
Single Shot ◽  
Velocity Analysis ◽  
Migration Velocity Analysis

Seismic imaging produces images of contrasts in physical parameters in the subsurface, e.g., velocity or impedance. To build such images, a background model describing the wave kinematics in the earth is necessary. In practice, the structural image and background velocity model are unknown and have to be estimated from the acquired data. Migration velocity analysis deals with estimation of the background model in the framework of seismic migration and relies on two main elements: data redundancy and invariance of the structures with respect to different seismic experiments. Because all the experiments probe the same model, the reflectors must be invariant in suitable domains (e.g., shots or reflection angle); the semblance principle is the tool used to measure the invariance of a set of multiple images. We measure the similarity of the structural features between pairs of single-shot migrated images obtained from adjacent experiments. By using the estimated warping vector field between two migrated images, we construct an image perturbation which describes the difference in reflectivity observed by two shots. We derive an expression for the image perturbation that drives a migration velocity analysis procedure based on a linearization of the wave-equation with respect to the model parameters. Synthetic 2D examples show promising results in retrieving errors in the velocity model. This methodology can be directly applied to 3D.

Joint inversion of the reflectivity and the velocity model

Geophysics ◽  
2021 ◽  
pp. 1-68
Author(s):  
Alejandro Cabrales-Vargas ◽  
Rahul Sarkar ◽  
Biondo L. Biondi ◽  
Robert G. Clapp
Keyword(s):  
Wave Equation ◽  
Field Data ◽  
Joint Inversion ◽  
Waveform Inversion ◽  
Velocity Model ◽  
Migration Velocity ◽  
Background Model ◽  
Seismic Events ◽  
Velocity Analysis ◽  
Migration Velocity Analysis

During linearized waveform inversion, the presence of small inaccuracies in the background subsurface model can lead to unfocused seismic events in the final image. The effect on the amplitude can mislead the interpretation. We present a joint inversion scheme in the model domain of the reflectivity and the background velocity model. The idea is to unify the inversion of the background and the reflectivity model into a single framework instead of treating them as decoupled problems. We show that with this method, we can obtain a better estimate of the reflectivity than that obtained with conventional linearized waveform inversion. Conversely, the background model is improved by the joint inversion with the reflectivity in comparison with wave-equation migration velocity analysis. We perform tests on 2D synthetics and 3D field data that demonstrate both benefits.

The finite-frequency sensitivity kernel for migration residual moveout and its applications in migration velocity analysis

Geophysics ◽  
2008 ◽  
Vol 73 (6) ◽  
pp. S241-S249 ◽  
Author(s):  
Xiao-Bi Xie ◽  
Hui Yang
Keyword(s):  
Wave Equation ◽  
Synthetic Data ◽  
Velocity Model ◽  
Migration Velocity ◽  
Data Sets ◽  
Velocity Analysis ◽  
Depth Migration ◽  
Sensitivity Kernel ◽  
Migration Velocity Analysis ◽  
Tomography Method

We have derived a broadband sensitivity kernel that relates the residual moveout (RMO) in prestack depth migration (PSDM) to velocity perturbations in the migration-velocity model. We have compared the kernel with the RMO directly measured from the migration image. The consistency between the sensitivity kernel and the measured sensitivity map validates the theory and the numerical implementation. Based on this broadband sensitivity kernel, we propose a new tomography method for migration-velocity analysis and updating — specifically, for the shot-record PSDM and shot-index common-image gather. As a result, time-consuming angle-domain analysis is not required. We use a fast one-way propagator and multiple forward scattering and single backscattering approximations to calculate the sensitivity kernel. Using synthetic data sets, we can successfully invert velocity perturbations from the migration RMO. This wave-equation-based method naturally incorporates the wave phenomena and is best teamed with the wave-equation migration method for velocity analysis. In addition, the new method maintains the simplicity of the ray-based velocity analysis method, with the more accurate sensitivity kernels replacing the rays.

Subsalt velocity estimation by target-oriented wave-equation migration velocity analysis: A 3D field-data example

Geophysics ◽  
2013 ◽  
Vol 78 (1) ◽  
pp. U19-U29 ◽  
Author(s):  
Yaxun Tang ◽  
Biondo Biondi
Keyword(s):  
Wave Equation ◽  
Field Data ◽  
Velocity Model ◽  
Migration Velocity ◽  
Velocity Estimation ◽  
Velocity Analysis ◽  
Data Set ◽  
Generalized Born ◽  
Migration Velocity Analysis ◽  
Recorded Data

We apply target-oriented wave-equation migration velocity analysis to a 3D field data set acquired from the Gulf of Mexico. Instead of using the original surface-recorded data set, we use a new data set synthesized specifically for velocity analysis to update subsalt velocities. The new data set is generated based on an initial unfocused target image and by a novel application of 3D generalized Born wavefield modeling, which correctly preserves velocity kinematics by modeling zero and nonzero subsurface-offset-domain images. The target-oriented inversion strategy drastically reduces the data size and the computation domain for 3D wave-equation migration velocity analysis, greatly improving its efficiency and flexibility. We apply differential semblance optimization (DSO) using the synthesized new data set to optimize subsalt velocities. The updated velocity model significantly improves the continuity of subsalt reflectors and yields flattened angle-domain common-image gathers.

Numeric implementation of wave-equation migration velocity analysis operators

Geophysics ◽  
2008 ◽  
Vol 73 (5) ◽  
pp. VE145-VE159 ◽  
Author(s):  
Paul Sava ◽  
Ioan Vlad
Keyword(s):  
Wave Equation ◽  
Velocity Model ◽  
Migration Velocity ◽  
Velocity Estimation ◽  
Velocity Analysis ◽  
Velocity Models ◽  
Image Optimization ◽  
Migration Velocity Analysis ◽  
Zero Offset ◽  
Wavefield Extrapolation

Wave-equation migration velocity analysis (MVA) is a technique similar to wave-equation tomography because it is designed to update velocity models using information derived from full seismic wavefields. On the other hand, wave-equation MVA is similar to conventional, traveltime-based MVA because it derives the information used for model updates from properties of migrated images, e.g., focusing and moveout. The main motivation for using wave-equation MVA is derived from its consistency with the corresponding wave-equation migration, which makes this technique robust and capable of handling multipathing characterizing media with large and sharp velocity contrasts. The wave-equation MVA operators are constructed using linearizations of conventional wavefield extrapolation operators, assuming small perturbations relative to the background velocity model. Similar to typical wavefield extrapolation operators, the wave-equation MVA operators can be implemented in the mixed space-wavenumber domain using approximations of differentorders of accuracy. As for wave-equation migration, wave-equation MVA can be formulated in different imaging frameworks, depending on the type of data used and image optimization criteria. Examples of imaging frameworks correspond to zero-offset migration (designed for imaging based on focusing properties of the image), survey-sinking migration (designed for imaging based on moveout analysis using narrow-azimuth data), and shot-record migration (also designed for imaging based on moveout analysis, but using wide-azimuth data). The wave-equation MVA operators formulated for the various imaging frameworks are similar because they share elements derived from linearizations of the single square-root equation. Such operators represent the core of iterative velocity estimation based on diffraction focusing or semblance analysis, and their applicability in practice requires efficient and accurate implementation. This tutorial concentrates strictly on the numeric implementation of those operators and not on their use for iterative migration velocity analysis.

Migration velocity analysis using traveltime wavefield tomography

Geophysics ◽  
2012 ◽  
Vol 77 (5) ◽  
pp. U73-U85 ◽  
Author(s):  
Saleh M. Al-Saleh ◽  
Jianwu Jiao
Keyword(s):  
Wave Equation ◽  
Velocity Model ◽  
Migration Velocity ◽  
Equation Modeling ◽  
Velocity Analysis ◽  
Lateral Velocity ◽  
Traveltime Tomography ◽  
Velocity Models ◽  
Migration Velocity Analysis ◽  
Complex Structural

We introduce an integrated wave-equation technique for migration velocity analysis (MVA) that consists of three steps: (1) forming the extended data, (2) approximating the correct transmitted wavefield, and (3) using wavefield tomography to update the velocity model. In the first step, the crosscorrelation imaging condition is relaxed to produce other nonzero-lag common image gathers (CIG) that, combined, form a common image cube (CIC). Slicing the CIC at different crosscorrelation lags forms a series of CIGs. Flattened events will occur in the CIGs at a lag other than the zero-lag when an incorrect velocity model is used in the migration. In the second step, for each event on the CIG, we pick the focusing depth and crosscorrelation lag at which it is flattest. We then model a Green’s function by seeding a source at the focusing depth using one-way wave equation modeling, then shift the modeled wavefield with the focusing crosscorrelation lag. This process is repeated for the other primary events at different lateral and vertical positions. The result is a set of modeled data whose wavefield approximates the wavefield that would have been generated if the correct velocity model had been used to simulate these gathers. We then apply wavefield tomography on these data-driven modeled data to update the velocity model. Our inversion scheme is based on wave-equation traveltime tomography that can update the velocity model in the presence of large velocity errors and a complex environment. Tests on synthetic and real 2D seismic data confirm the method’s effectiveness in building velocity models in complex structural areas that have large lateral velocity variations.

Migration velocity analysis from locally coherent events in 2‐D laterally heterogeneous media, Part I: Theoretical aspects

Geophysics ◽  
2002 ◽  
Vol 67 (4) ◽  
pp. 1202-1212 ◽  
Author(s):  
Hervé Chauris ◽  
Mark S. Noble ◽  
Gilles Lambaré ◽  
Pascal Podvin
Keyword(s):  
Cost Function ◽  
Velocity Model ◽  
Migration Velocity ◽  
Velocity Estimation ◽  
Velocity Analysis ◽  
Seismic Reflection Data ◽  
Velocity Models ◽  
Migration Velocity Analysis ◽  
Paraxial Ray ◽  
The Cost

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.

Migration velocity analysis based on focusing diffractions generated using the CRS wavefront attributes: Application to real land data

Geophysics ◽  
2021 ◽  
pp. 1-50
Author(s):  
German Garabito ◽  
José Silas dos Santos Silva ◽  
Williams Lima
Keyword(s):  
Time Domain ◽  
Real Data ◽  
Normal Incidence ◽  
Velocity Model ◽  
Migration Velocity ◽  
Velocity Analysis ◽  
Velocity Models ◽  
Time Migration ◽  
Migration Velocity Analysis ◽  
The Time Domain

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.

Time-migration velocity analysis by image-wave propagation of common-image gathers

Geophysics ◽  
2008 ◽  
Vol 73 (5) ◽  
pp. VE161-VE171 ◽  
Author(s):  
J. Schleicher ◽  
J. C. Costa ◽  
A. Novais
Keyword(s):  
Wave Propagation ◽  
Seismic Event ◽  
Velocity Model ◽  
Migration Velocity ◽  
Velocity Analysis ◽  
Reference Velocity ◽  
Time Migration ◽  
Migration Velocity Analysis ◽  
Velocity Image ◽  
The Common

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.

Sign in / Sign up

Export Citation Format

Share Document