High-resolution seismic velocity analysis by sign-based weighted semblance

Geophysics ◽  
2021 ◽  
pp. 1-35
Author(s):  
M. Javad Khoshnavaz

Building an accurate velocity model plays a vital role in routine seismic imaging workflows. Normal-moveout-based seismic velocity analysis is a popular method to make the velocity models. However, traditional velocity analysis methodologies are not generally capable of handling amplitude variations across moveout curves, specifically polarity reversals caused by amplitude-versus-offset anomalies. I present a normal-moveout-based velocity analysis approach that circumvents this shortcoming by modifying the conventional semblance function to include polarity and amplitude correction terms computed using correlation coefficients of seismic traces in the velocity analysis scanning window with a reference trace. Thus, the proposed workflow is suitable for any class of amplitude-versus-offset effects. The approach is demonstrated to four synthetic data examples of different conditions and a field data consisting a common-midpoint gather. Lateral resolution enhancement using the proposed workflow is evaluated by comparison between the results from the workflow and the results obtained by the application of conventional semblance and three semblance-based velocity analysis algorithms developed to circumvent the challenges associated with amplitude variations across moveout curves, caused by seismic attenuation and class II amplitude-versus-offset anomalies. According to the obtained results, the proposed workflow is superior to all the presented workflows in handling such anomalies.

Geophysics ◽  
1993 ◽  
Vol 58 (1) ◽  
pp. 91-100 ◽  
Author(s):  
Claude F. Lafond ◽  
Alan R. Levander

Prestack depth migration still suffers from the problems associated with building appropriate velocity models. The two main after‐migration, before‐stack velocity analysis techniques currently used, depth focusing and residual moveout correction, have found good use in many applications but have also shown their limitations in the case of very complex structures. To address this issue, we have extended the residual moveout analysis technique to the general case of heterogeneous velocity fields and steep dips, while keeping the algorithm robust enough to be of practical use on real data. Our method is not based on analytic expressions for the moveouts and requires no a priori knowledge of the model, but instead uses geometrical ray tracing in heterogeneous media, layer‐stripping migration, and local wavefront analysis to compute residual velocity corrections. These corrections are back projected into the velocity model along raypaths in a way that is similar to tomographic reconstruction. While this approach is more general than existing migration velocity analysis implementations, it is also much more computer intensive and is best used locally around a particularly complex structure. We demonstrate the technique using synthetic data from a model with strong velocity gradients and then apply it to a marine data set to improve the positioning of a major fault.


Geophysics ◽  
2008 ◽  
Vol 73 (5) ◽  
pp. VE35-VE38 ◽  
Author(s):  
Jonathan Liu ◽  
Lorie Bear ◽  
Jerry Krebs ◽  
Raffaella Montelli ◽  
Gopal Palacharla

We have developed a new method to build seismic velocity models for complex structures. In our approach, we use a spatially nonuniform parameterization of the velocity model in tomography and a uniform grid representation of the same velocity model in ray tracing to generate the linear system of tomographic equations. Subsequently, a matrix transformation is applied to the system of equations to produce a new linear system of tomographic equations using nonuniform parameterization. In this way, we improved the stability of tomographic inversion without adding computing costs. We tested the effectiveness of our process on a 3D synthetic data example.


Geophysics ◽  
2018 ◽  
Vol 83 (5) ◽  
pp. B241-B252 ◽  
Author(s):  
Daniele Colombo ◽  
Diego Rovetta ◽  
Ersan Turkoglu

Seismic imaging in salt geology is complicated by highly contrasted velocity fields and irregular salt geometries, which cause complex seismic wavefield scattering. Although the imaging challenges can be addressed by advanced imaging algorithms, a fundamental problem remains in the determination of robust velocity fields in high-noise conditions. Conventional migration velocity analysis is often ineffective, and even the most advanced methods for depth-domain velocity analysis, such as full-waveform inversion, require starting from a good initial estimate of the velocity model to converge to a correct result. Nonseismic methods, such as electromagnetics, can help guide the generation of robust velocity models to be used for further processing. Using the multiphysics data acquired in the deepwater section of the Red Sea, we apply a controlled-source electromagnetic (CSEM) resistivity-regularized seismic velocity inversion for enhancing the velocity model in a complex area dominated by nappe-style salt tectonics. The integration is achieved by a rigorous approach of multiscaled inversions looping over model dimensions (1D first, followed by 3D), variable offsets and increasing frequencies, data-driven and interpretation-supported approaches, leading to a hierarchical inversion guided by a parameter sensitivity analysis. The final step of the integration consists of the inversion of seismic traveltimes subject to CSEM model constraints in which a common-structure coupling mechanism is used. Minimization is performed over the seismic data residuals and cross-gradient objective functions without inverting for the resistivity model, which is used as a reference for the seismic inversion (hierarchical approach). Results are demonstrated through depth imaging in which the velocity model derived through CSEM-regularized hierarchical inversion outperforms the results of a seismic-only derived velocity model.


Geophysics ◽  
1997 ◽  
Vol 62 (6) ◽  
pp. 1825-1838 ◽  
Author(s):  
Jun Ji

In areas with structurally complex geology, tomographic velocity analysis is often required to estimate velocities. In this paper I describe an algorithm for tomographic velocity estimation that uses plane‐wave synthesis imaging as a prestack migration. The classical iterative two‐step process (measures the traveltime errors with the current velocity model and then update the velocity model) is performed as follows. The events are picked in the image space after prestack migration with surface‐oriented plane‐wave synthesis imaging. the traveltime deviations are measured through residual‐moveout (RMO) velocity analysis in common‐surface‐location (CSL) gathers obtained by reflector‐oriented plane‐wave synthesis imaging, and the velocity update is calculated by inverting the traveltime deviations through a conjugate gradient. The results from synthetic data indicate that the tomographic method successfully estimates interval‐velocity models that lead to depth‐migrated images with no residual moveout.


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.


2018 ◽  
Vol 2018 ◽  
pp. 1-11
Author(s):  
Vladimir Sabinin

Some new computational techniques are suggested for estimating symmetry axis azimuth of fractures in the viscoelastic anisotropic target layer in the framework of QVOA analysis (Quality factor Versus Offset and Azimuth). The different QVOA techniques are compared using synthetic viscoelastic surface reflected data with and without noise. I calculated errors for these techniques which depend on different sets of azimuths and intervals of offsets. Superiority of the high-order “enhanced general” and “cubic” techniques is shown. The high-quality QVOA techniques are compared with one of the high-quality AVOA techniques (Amplitude Versus Offset and Azimuth) in the synthetic data with noise and attenuation. Results are comparable.


Geophysics ◽  
2003 ◽  
Vol 68 (3) ◽  
pp. 1008-1021 ◽  
Author(s):  
Frederic Billette ◽  
Soazig Le Bégat ◽  
Pascal Podvin ◽  
Gilles Lambaré

Stereotomography is a new velocity estimation method. This tomographic approach aims at retrieving subsurface velocities from prestack seismic data. In addition to traveltimes, the slope of locally coherent events are picked simultaneously in common offset, common source, common receiver, and common midpoint gathers. As the picking is realized on locally coherent events, they do not need to be interpreted in terms of reflection on given interfaces, but may represent diffractions or reflections from anywhere in the image. In the high‐frequency approximation, each one of these events corresponds to a ray trajectory in the subsurface. Stereotomography consists of picking and analyzing these events to update both the associated ray paths and velocity model. In this paper, we describe the implementation of two critical features needed to put stereotomography into practice: an automatic picking tool and a robust multiscale iterative inversion technique. Applications to 2D reflection seismic are presented on synthetic data and on a 2D line extracted from a 3D towed streamer survey shot in West Africa for TotalFinaElf. The examples demonstrate that the method requires only minor human intervention and rapidly converges to a geologically plausible velocity model in these two very different and complex velocity regimes. The quality of the velocity models is verified by prestack depth migration results.


Geophysics ◽  
2021 ◽  
pp. 1-73
Author(s):  
Hani Alzahrani ◽  
Jeffrey Shragge

Data-driven artificial neural networks (ANNs) offer a number of advantages over conventional deterministic methods in a wide range of geophysical problems. For seismic velocity model building, judiciously trained ANNs offer the possibility of estimating high-resolution subsurface velocity models. However, a significant challenge of ANNs is training generalization, which is the ability of an ANN to apply the learning from the training process to test data not previously encountered. In the context of velocity model building, this means learning the relationship between velocity models and the corresponding seismic data from a set of training data, and then using acquired seismic data to accurately estimate unknown velocity models. We ask the following question: what type of velocity model structures need be included in the training process so that the trained ANN can invert seismic data from a different (hypothetical) geological setting? To address this question, we create four sets of training models: geologically inspired and purely geometrical, both with and without background velocity gradients. We find that using geologically inspired training data produce models with well-delineated layer interfaces and fewer intra-layer velocity variations. The absence of a certain geological structure in training models, though, hinders the ANN's ability to recover it in the testing data. We use purely geometric training models consisting of square blocks of varying size to demonstrate the ability of ANNs to recover reasonable approximations of flat, dipping, and curved interfaces. However, the predicted models suffer from intra-layer velocity variations and non-physical artifacts. Overall, the results successfully demonstrate the use of ANNs in recovering accurate velocity model estimates, and highlight the possibility of using such an approach for the generalized seismic velocity inversion problem.


Sign in / Sign up

Export Citation Format

Share Document