scholarly journals Wave-equation migration velocity analysis using plane-wave common-image gathers

Geophysics ◽  
2017 ◽  
Vol 82 (5) ◽  
pp. S327-S340 ◽  
Author(s):  
Bowen Guo ◽  
Gerard T. Schuster
Keyword(s):  
Wave Equation ◽  
Plane Wave ◽  
Computational Cost ◽  
Migration Velocity ◽  
Memory Storage ◽  
Velocity Analysis ◽  
Data Set ◽  
Image Domain ◽  
Migration Velocity Analysis ◽  
Wave Migration

Wave-equation migration velocity analysis (WEMVA) based on subsurface-offset, angle domain, or time-lag common-image gathers (CIGs) requires significant computational and memory resources because it computes higher dimensional migration images in the extended image domain. To mitigate this problem, we have developed a WEMVA method using plane-wave CIGs. Plane-wave CIGs reduce computational cost and memory storage because they are directly calculated from prestack plane-wave migration and the number of plane waves is often much smaller than the number of shots. In the case of an inaccurate migration velocity, the moveout of plane-wave CIGs is automatically picked by a semblance analysis method, which is then linked to the migration velocity update by a connective function. Numerical tests on two synthetic data sets and a field data set validate the efficiency and effectiveness of this method.

Joint PP and PS plane-wave wave-equation migration-velocity analysis

Geophysics ◽  
2019 ◽  
Vol 84 (4) ◽  
pp. R507-R525 ◽  
Author(s):  
Zongcai Feng ◽  
Bowen Guo ◽  
Lianjie Huang
Keyword(s):  
Wave Equation ◽  
Plane Wave ◽  
Wave Velocity ◽  
Migration Velocity ◽  
Relative Depth ◽  
Velocity Analysis ◽  
Analysis Method ◽  
S Wave ◽  
Velocity Models ◽  
Migration Velocity Analysis

Conventional joint PP and PS velocity analysis is based on ray tomography. We develop a joint PP and PS wave-equation migration-velocity-analysis method using plane-wave common-image gathers (CIGs) to produce accurate P- and S-wave velocity models. The objective function of our new method consists of three terms: The first and second terms penalize the moveout residuals computed from PP and PS plane-wave CIGs, respectively, and the third term constrains the nonzero relative depth shifts between the PP and PS migration images. The moveout of plane-wave CIGs is automatically picked using a semblance analysis method, and the relative depth shifts between the PP and PS images are automatically computed using dynamic warping or manually picking the depths of certain primary reflectors. The moveout residuals and the relative depth shifts are transformed into weighted image perturbations, and they are then projected into the velocity models to update the P- and S-wave velocity models using the scalar-wave equations and their linearized forms. Numerical tests with synthetic and multicomponent field data demonstrate that our method can simultaneously invert for accurate P- and S-wave velocity models for elastic migration.

scholarly journals Automatic wave-equation migration velocity analysis by focusing subsurface virtual sources

Geophysics ◽  
2018 ◽  
Vol 83 (2) ◽  
pp. U1-U8 ◽  
Author(s):  
Bingbing Sun ◽  
Tariq Alkhalifah
Keyword(s):  
Wave Equation ◽  
Model Building ◽  
Waveform Inversion ◽  
Synthetic Data ◽  
Real Data ◽  
Migration Velocity ◽  
Velocity Analysis ◽  
Data Set ◽  
Migration Velocity Analysis ◽  
Band Limited

Macro-velocity model building is important for subsequent prestack depth migration and full-waveform inversion. Wave-equation migration velocity analysis uses the band-limited waveform to invert for velocity. Normally, inversion would be implemented by focusing the subsurface offset common-image gathers. We reexamine this concept with a different perspective: In the subsurface offset domain, using extended Born modeling, the recorded data can be considered as invariant with respect to the perturbation of the position of the virtual sources and velocity at the same time. A linear system connecting the perturbation of the position of those virtual sources and velocity is derived and solved subsequently by the conjugate gradient method. In theory, the perturbation of the position of the virtual sources is given by the Rytov approximation. Thus, compared with the Born approximation, it relaxes the dependency on amplitude and makes the proposed method more applicable for real data. We determined the effectiveness of the approach by applying the proposed method on isotropic and anisotropic vertical transverse isotropic synthetic data. A real data set example verifies the robustness of the proposed method.

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.

Wave-equation migration velocity analysis for VTI models

Geophysics ◽  
2014 ◽  
Vol 79 (3) ◽  
pp. WA59-WA68 ◽  
Author(s):  
Yunyue Li ◽  
Biondo Biondi ◽  
Robert Clapp ◽  
Dave Nichols
Keyword(s):  
Wave Equation ◽  
Synthetic Data ◽  
Rock Physics ◽  
Migration Velocity ◽  
Velocity Analysis ◽  
Sources Of Information ◽  
Data Set ◽  
Anisotropic Models ◽  
Migration Velocity Analysis ◽  
Power Maximization

Anisotropic models are needed for wave simulation and inversion where a complex geologic environment exists. We extended the theory of wave equation migration velocity analysis to build vertical transverse isotropic models. Because of the ambiguity between depth and [Formula: see text] in the acoustic regime, we assumed [Formula: see text] can be accurately obtained from other sources of information, and inverted for the NMO slowness and the anellipticity parameter [Formula: see text]. We combined the differential semblance optimization objective function with the stacking power maximization to evaluate the focusing of the prestack image in the subsurface-offset domain. To regularize the multiparameter inversion, we built a framework to adapt the geologic and the rock physics information to guide the updates in NMO slowness and [Formula: see text]. This regularization step was crucial to stabilize the inversion and to produce geologically meaningful results. We tested the proposed approach on a synthetic data set and a 2D Gulf of Mexico data set starting with a fairly good initial anisotropic model. The inversion results revealed shallow anomalies collocated in NMO velocity and [Formula: see text] and improved the continuity and the resolution of the final stacked images.

Angle-domain common-image gathers from plane-wave least-squares reverse time migration

Geophysics ◽  
2021 ◽  
pp. 1-60
Author(s):  
Chuang Li ◽  
Zhaoqi Gao ◽  
Jinghuai Gao ◽  
Feipeng Li ◽  
Tao Yang
Keyword(s):  
Inverse Problem ◽  
Plane Wave ◽  
Least Squares ◽  
Migration Velocity ◽  
Velocity Analysis ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Time Migration ◽  
Migration Velocity Analysis ◽  
Amplitude Versus

Angle-domain common-image gathers (ADCIGs) that can be used for migration velocity analysis and amplitude versus angle analysis are important for seismic exploration. However, because of limited acquisition geometry and seismic frequency band, the ADCIGs extracted by reverse time migration (RTM) suffer from illumination gaps, migration artifacts, and low resolution. We have developed a reflection angle-domain pseudo-extended plane-wave least-squares RTM method for obtaining high-quality ADCIGs. We build the mapping relations between the ADCIGs and the plane-wave sections using an angle-domain pseudo-extended Born modeling operator and an adjoint operator, based on which we formulate the extraction of ADCIGs as an inverse problem. The inverse problem is iteratively solved by a preconditioned stochastic conjugate gradient method, allowing for reduction in computational cost by migrating only a subset instead of the whole dataset and improving image quality thanks to preconditioners. Numerical tests on synthetic and field data verify that the proposed method can compensate for illumination gaps, suppress migration artifacts, and improve resolution of the ADCIGs and the stacked images. Therefore, compared with RTM, the proposed method provides a more reliable input for migration velocity analysis and amplitude versus angle analysis. Moreover, it also provides much better stacked images for seismic interpretation.

Sign in / Sign up

Export Citation Format

Share Document