Ray-based seismic modeling of geologic models: Understanding and analyzing seismic images efficiently

2015 ◽  
Vol 3 (4) ◽  
pp. SAC71-SAC89 ◽  
Author(s):  
Isabelle Lecomte ◽  
Paul Lubrano Lavadera ◽  
Ingrid Anell ◽  
Simon J. Buckley ◽  
Daniel W. Schmid ◽  
...  
Keyword(s):  
Normal Incidence ◽  
Velocity Model ◽  
Convolution Operators ◽  
Seismic Modeling ◽  
Lateral Velocity ◽  
Velocity Models ◽  
Convolution Model ◽  
Seismic Sections ◽  
Efficient Option ◽  
Ray Methods

Often, interpreters only have access to seismic sections and, at times, well data, when making an interpretation of structures and depositional features in the subsurface. The validity of the final interpretation is based on how well the seismic data are able to reproduce the actual geology, and seismic modeling can help constrain that. Ideally, modeling should create complete seismograms, which is often best achieved by finite-difference modeling with postprocessing to produce synthetic seismic sections for comparison purposes. Such extensive modeling is, however, not routinely affordable. A far more efficient option, using the simpler 1D convolution model with reflectivity logs extracted along verticals in velocity models, generates poor modeling results when lateral velocity variations are expected. A third and intermediate option is to use the various ray-based approaches available, which are efficient and flexible. However, standard ray methods, such as the normal-incidence point for unmigrated poststack sections or image rays for simulating time-migrated poststack results, cannot deal with complex and detailed targets, and will not reproduce the realistic (3D) resolution effects of seismic imaging. Nevertheless, ray methods can also be used to estimate 3D spatial prestack convolution operators, so-called point-spread functions. These are functions of the survey, velocity model, and wavelet, among others, and therefore they include 3D angle-dependent illumination and resolution effects. Prestack depth migration images are thus rapidly simulated by spatial convolution with detailed 3D reflectivity models, which goes far beyond the limits of 1D convolution modeling. This 3D convolution modeling should allow geologists to better assess their interpretations and draw more definitive conclusions.

Seismic modeling and expression of common fold-thrust structures

2020 ◽  
Vol 8 (1) ◽  
pp. T55-T65
Author(s):  
Jianjun Li ◽  
Shankar Mitra
Keyword(s):  
Velocity Model ◽  
Fault Geometry ◽  
Seismic Modeling ◽  
Fault Propagation ◽  
Lateral Velocity ◽  
Slip Ratio ◽  
Velocity Models ◽  
Time Migration ◽  
The Common

We have conducted seismic modeling of common fold-thrust structures to understand the common geologic parameters influencing seismic data and to understand the common pitfalls associated with interpreting prestack time migration (PSTM) and prestack depth migration (PSDM) data. Mode 1 fault-bend folds are generally well-imaged in PSTM data, provided the correct migration velocities are used for the dipping back and front limbs. Seismic pull-ups of the footwall related to lateral velocity variations can result in problems in interpreting the fault geometry and the subthrust area underlying the crest. Fault-tip fault-propagation folds also show significant footwall pull-ups and show poor to no imaging of the steep front limbs. The geometry of trishear fault-propagation folds is dependent on the maximum slip on the fault (S) and the fault propagation to slip ratio (P/S ratio). We found that the slip has a strong influence on the dip of the front limb and therefore the quality of imaging whereas the P/S ratio, which controls the degree of folding versus thrust faulting, has only a secondary effect. For the front limb, only the area near the synclinal axial plan is well-imaged, so that the fault geometry and extent of propagation are typically difficult to interpret. The front limb dips are also sensitive to the accuracy of the rms velocity model used for migration. Lower velocities result in steeper dipping reflectors, whereas higher velocities result in shallower dips. In general, PSDM provides better imaging of the structures; however, the accuracy and quality of the image are dependent on the velocity models and interpretation derived from the PSTM data.

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.

scholarly journals Moving from 1-D to 3-D velocity model: automated waveform-based earthquake moment tensor inversion in the Los Angeles region

2019 ◽  
Vol 220 (1) ◽  
pp. 218-234 ◽  
Author(s):  
Xin Wang ◽  
Zhongwen Zhan
Keyword(s):  
Los Angeles ◽  
Moment Tensor ◽  
Velocity Model ◽  
Primary Control ◽  
Lateral Velocity ◽  
Real Time System ◽  
Moment Tensors ◽  
Velocity Models ◽  
Double Couple ◽  
The Moment

SUMMARY Earthquake focal mechanisms put primary control on the distribution of ground motion, and also bear on the stress state of the crust. Most routine focal mechanism catalogues still use 1-D velocity models in inversions, which may introduce large uncertainties in regions with strong lateral velocity heterogeneities. In this study, we develop an automated waveform-based inversion approach to determine the moment tensors of small-to-medium-sized earthquakes using 3-D velocity models. We apply our approach in the Los Angeles region to produce a new moment tensor catalogue with a completeness of ML ≥ 3.5. The inversions using the Southern California Earthquake Center Community Velocity Model (3D CVM-S4.26) significantly reduces the moment tensor uncertainties, mainly owing to the accuracy of the 3-D velocity model in predicting both the phases and the amplitudes of the observed seismograms. By comparing the full moment tensor solutions obtained using 1-D and 3-D velocity models, we show that the percentages of non-double-couple components decrease dramatically with the usage of 3-D velocity model, suggesting that large fractions of non-double-couple components from 1-D inversions are artifacts caused by unmodelled 3-D velocity structures. The new catalogue also features more accurate focal depths and moment magnitudes. Our highly accurate, efficient and automatic inversion approach can be expanded in other regions, and can be easily implemented in near real-time system.

Datum correction by wave‐equation extrapolation

Geophysics ◽  
1988 ◽  
Vol 53 (10) ◽  
pp. 1311-1322 ◽  
Author(s):  
V. Shtivelman ◽  
A. Canning
Keyword(s):  
Wave Equation ◽  
Synthetic Data ◽  
Real Data ◽  
Static Solution ◽  
Velocity Model ◽  
Lateral Velocity ◽  
Static Correction ◽  
Slow Velocity ◽  
Velocity Variations ◽  
Seismic Sections

Seismic sections are usually datum corrected by static shifting. For small differences in elevation and slow velocity variations between the input datum and the output datum, static shifting is a sufficiently accurate datum correction procedure. However, for significant differences in elevations and a more complicated velocity model, the accuracy of the static solution may prove to be insufficient; and a more exact method should be used. In this paper, we study the limitations of the static method of datum correction and develop simple and effective extrapolation schemes based on the wave equation, schemes which lead to more accurate datum correction. The distortions of seismic events caused by static correction are illustrated by a number of simple examples. To reduce the distortions, we propose a number of extrapolation schemes based on the asymptories of the Kirchhoff integral solution of the 2-D scalar wave equation. Application of the extrapolation algorithms to synthetic data shows that they provide accurate datum corrections even for a nonplanar input datum and vertical and lateral velocity variations. The algorithms have been successfully applied to real data.

Exploration beyond seismic: The role of electromagnetics and gravity gradiometry in deep water subsalt plays of the Red Sea

2014 ◽  
Vol 2 (3) ◽  
pp. SH33-SH53 ◽  
Author(s):  
Daniele Colombo ◽  
Gary McNeice ◽  
Nickolas Raterman ◽  
Mike Zinger ◽  
Diego Rovetta ◽  
...  
Keyword(s):  
Red Sea ◽  
Model Building ◽  
Three Dimensional ◽  
Velocity Model ◽  
Seismic Exploration ◽  
Gravity Gradiometry ◽  
3D Inversion ◽  
Lateral Velocity ◽  
Velocity Models ◽  
Basement Structures

The Red Sea is characterized by thick salt sequences representing a seal for potential hydrocarbon accumulations within Tertiary formations deposited over deep basement structures. The Red Sea “salt” is characterized by halite concentrations embedded in layered evaporite sequences composed of evaporite and clastic lithologies. Salt complicates seismic exploration efforts in the Red Sea by generating vertical and lateral velocity variations that are difficult to estimate by seismic methods alone. In these conditions, the exploration challenges of independently imaging the subsalt section and provide enhanced velocity model building capabilities were addressed by a multigeophysics strategy involving marine electromagnetics (magnetotellurics and controlled source electromagnetics [CSEM]) and gravity gradiometry surveys colocated with wide azimuth seismic. Three-dimensional inversion of MT and CSEM is performed first with minimal a priori constraints and then by including variable amounts of interpretation in the starting models. The internal variations in the evaporitic overburden, the subsalt, and the basement structures are independently imaged by combined electromagnetic methods and confirmed by new drilling results. CSEM, in particular, provides unprecedented detail of the internal structures within the salt overburden while magnetotellurics provides excellent reconstruction of the base of salt and basement. Gravity gradiometry shows primary sensitivity to the basement and the corresponding 3D inversion provides density distributions structurally consistent with the resistivity volumes. The common-structure, multiparameter models obtained from 3D inversion deliver additional aid to seismic interpreters to further derisk exploration in the Red Sea and provide additional detail to depth imaging velocity models. The reciprocal consistency of the obtained results show promises for extending the work to more analytical integration with seismic such as provided by joint geophysical inversion.

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.

First-arrival traveltime inversion of seismic diving waves observed on undulant surface

2021 ◽  
Vol 225 (2) ◽  
pp. 1020-1031
Author(s):  
Huachen Yang ◽  
Jianzhong Zhang ◽  
Kai Ren ◽  
Changbo Wang
Keyword(s):  
Turning Points ◽  
Analytical Formula ◽  
Velocity Model ◽  
Western China ◽  
Inversion Method ◽  
Mountain Area ◽  
Traveltime Inversion ◽  
Static Correction ◽  
Velocity Models ◽  
First Arrival

SUMMARY A non-iterative first-arrival traveltime inversion method (NFTI) is proposed for building smooth velocity models using seismic diving waves observed on irregular surface. The new ray and traveltime equations of diving waves propagating in smooth media with undulant observation surface are deduced. According to the proposed ray and traveltime equations, an analytical formula for determining the location of the diving-wave turning points is then derived. Taking the influence of rough topography on first-arrival traveltimes into account, the new equations for calculating the velocities at turning points are established. Based on these equations, a method is proposed to construct subsurface velocity models from the observation surface downward to the bottom using the first-arrival traveltimes in common offset gathers. Tests on smooth velocity models with rugged topography verify the validity of the established equations, and the superiority of the proposed NFTI. The limitation of the proposed method is shown by an abruptly-varying velocity model example. Finally, the NFTI is applied to solve the static correction problem of the field seismic data acquired in a mountain area in the western China. The results confirm the effectivity of the proposed NFTI.

scholarly journals Earthquake location in tectonic structures of the Alpine Chain: the case of the Constance Lake (Central Europe) seismic sequence

2021 ◽  
Author(s):  
Francesca D’Ajello Caracciolo ◽  
Rodolfo Console
Keyword(s):  
Central Europe ◽  
Velocity Model ◽  
Moho Depth ◽  
Seismic Sequence ◽  
Tectonic Structures ◽  
Study Region ◽  
Waveform Data ◽  
Seismic Stations ◽  
Velocity Models ◽  
Master Event

AbstractA set of four magnitude Ml ≥ 3.0 earthquakes including the magnitude Ml = 3.7 mainshock of the seismic sequence hitting the Lake Constance, Southern Germany, area in July–August 2019 was studied by means of bulletin and waveform data collected from 86 seismic stations of the Central Europe-Alpine region. The first single-event locations obtained using a uniform 1-D velocity model, and both fixed and free depths, showed residuals of the order of up ± 2.0 s, systematically affecting stations located in different areas of the study region. Namely, German stations to the northeast of the epicenters and French stations to the west exhibit negative residuals, while Italian stations located to the southeast are characterized by similarly large positive residuals. As a consequence, the epicentral coordinates were affected by a significant bias of the order of 4–5 km to the NNE. The locations were repeated applying a method that uses different velocity models for three groups of stations situated in different geological environments, obtaining more accurate locations. Moreover, the application of two methods of relative locations and joint hypocentral determination, without improving the absolute location of the master event, has shown that the sources of the four considered events are separated by distances of the order of one km both in horizontal coordinates and in depths. A particular attention has been paid to the geographical positions of the seismic stations used in the locations and their relationship with the known crustal features, such as the Moho depth and velocity anomalies in the studied region. Significant correlations between the observed travel time residuals and the crustal structure were obtained.

Interferometric imaging condition for wave-equation migration

Geophysics ◽  
2008 ◽  
Vol 73 (2) ◽  
pp. S47-S61 ◽  
Author(s):  
Paul Sava ◽  
Oleg Poliannikov
Keyword(s):  
Large Scale ◽  
Wigner Distribution ◽  
Random Noise ◽  
Velocity Model ◽  
Distribution Functions ◽  
Small Scale ◽  
Interferometric Imaging ◽  
Imaging Data ◽  
Imaging Condition ◽  
Velocity Models

The fidelity of depth seismic imaging depends on the accuracy of the velocity models used for wavefield reconstruction. Models can be decomposed in two components, corresponding to large-scale and small-scale variations. In practice, the large-scale velocity model component can be estimated with high accuracy using repeated migration/tomography cycles, but the small-scale component cannot. When the earth has significant small-scale velocity components, wavefield reconstruction does not completely describe the recorded data, and migrated images are perturbed by artifacts. There are two possible ways to address this problem: (1) improve wavefield reconstruction by estimating more accurate velocity models and image using conventional techniques (e.g., wavefield crosscorrelation) or (2) reconstruct wavefields with conventional methods using the known background velocity model but improve the imaging condition to alleviate the artifacts caused by the imprecise reconstruction. Wedescribe the unknown component of the velocity model as a random function with local spatial correlations. Imaging data perturbed by such random variations is characterized by statistical instability, i.e., various wavefield components image at wrong locations that depend on the actual realization of the random model. Statistical stability can be achieved by preprocessing the reconstructed wavefields prior to the imaging condition. We use Wigner distribution functions to attenuate the random noise present in the reconstructed wavefields, parameterized as a function of image coordinates. Wavefield filtering using Wigner distribution functions and conventional imaging can be lumped together into a new form of imaging condition that we call an interferometric imaging condition because of its similarity to concepts from recent work on interferometry. The interferometric imaging condition can be formulated both for zero-offset and for multioffset data, leading to robust, efficient imaging procedures that effectively attenuate imaging artifacts caused by unknown velocity models.

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.

Sign in / Sign up

Export Citation Format

Share Document