Full-traveltime inversion

Geophysics ◽  
2016 ◽  
Vol 81 (5) ◽  
pp. R261-R274 ◽  
Author(s):  
Yi Luo ◽  
Yue Ma ◽  
Yan Wu ◽  
Hongwei Liu ◽  
Lei Cao
Keyword(s):  
Wave Equation ◽  
Seismic Data ◽  
Low Frequency ◽  
Velocity Model ◽  
New Wave ◽  
Severe Problem ◽  
Reflected Waves ◽  
Traveltime Inversion ◽  
Velocity Models ◽  
Kinematic Information

Many previously published wave-equation-based methods, which attempt to automatically invert traveltime or kinematic information in seismic data or migrated gathers for smooth velocities, suffer a common and severe problem — the inversions are involuntarily and unconsciously hijacked by amplitude information. To overcome this problem, we have developed a new wave-equation-based traveltime inversion methodology, referred to as full-traveltime (i.e., fully dependent on traveltime) inversion (FTI), to automatically estimate a kinematically accurate velocity model from seismic data. The key idea of FTI is to make the inversion fully dependent on traveltime information, and thus prevent amplitude interference during inversion. Under the assumption that velocity perturbations cause only traveltime changes, we have derived the FTI method in the data and image domains, which are applicable to transmitted arrivals and reflected waves, respectively. FTI does not require an accurate initial velocity model or low-frequency seismic data. Synthetic and field data tests demonstrate that FTI produces satisfactory inversion results, even when using constant velocity models as initials.

Traveltime information-based wave-equation inversion

Geophysics ◽  
2009 ◽  
Vol 74 (6) ◽  
pp. WCC27-WCC36 ◽  
Author(s):  
Yu Zhang ◽  
Daoliu Wang
Keyword(s):  
Wave Equation ◽  
Conventional Method ◽  
Seismic Data ◽  
Waveform Inversion ◽  
Back Propagation ◽  
Velocity Model ◽  
Data Fitting ◽  
Inversion Method ◽  
New Wave ◽  
Seismic Record

We propose a new wave-equation inversion method that mainly depends on the traveltime information of the recorded seismic data. Unlike the conventional method, we first apply a [Formula: see text] transform to the seismic data to form the delayed-shot seismic record, back propagate the transformed data, and then invert the velocity model by maximizing the wavefield energy around the shooting time at the source locations. Data fitting is not enforced during the inversion, so the optimized velocity model is obtained by best focusing the source energy after a back propagation. Therefore, inversion accuracy depends only on the traveltime information embedded in the seismic data. This method may overcome some practical issues of waveform inversion; in particular, it relaxes the dependency of the seismic data amplitudes and the source wavelet.

Inversion of seismic refraction and reflection data for building long-wavelength velocity models

Geophysics ◽  
2015 ◽  
Vol 80 (2) ◽  
pp. R81-R93 ◽  
Author(s):  
Haiyang Wang ◽  
Satish C. Singh ◽  
Francois Audebert ◽  
Henri Calandra
Keyword(s):  
Wave Equation ◽  
Short Wavelength ◽  
Velocity Model ◽  
Density Model ◽  
Data Set ◽  
Reflected Waves ◽  
Long Wavelength ◽  
Velocity Models ◽  
Computation Efficiency ◽  
Refracted Waves

Long-wavelength velocity model building is a nonlinear process. It has traditionally been achieved without appealing to wave-equation-based approaches for combined refracted and reflected waves. We developed a cascaded wave-equation tomography method in the data domain, taking advantage of the information contained in the reflected and refracted waves. The objective function was the traveltime residual that maximized the crosscorrelation function between real and synthetic data. To alleviate the nonlinearity of the inversion problem, refracted waves were initially used to provide vertical constraints on the velocity model, and reflected waves were then included to provide lateral constraints. The use of reflected waves required scale separation. We separated the long- and short-wavelength subsurface structures into velocity and density models, respectively. The velocity model update was restricted to long wavelengths during the wave-equation tomography, whereas the density model was used to absorb all the short-wavelength impedance contrasts. To improve the computation efficiency, the density model was converted into the zero-offset traveltime domain, where it was invariant to changes of the long-wavelength velocity model. After the wave-equation tomography has derived an optimized long-wavelength velocity model, full-waveform inversion was used to invert all the data to retrieve the short-wavelength velocity structures. We developed our method in two synthetic tests and then applied it to a marine field data set. We evaluated the results of the use of refracted and reflected waves, which was critical for accurately building the long-wavelength velocity model. We showed that our wave-equation tomography strategy was robust for the real data application.

Joint wave-equation traveltime inversion of diving/direct and reflected waves for the P- and S-wave velocity macromodel building

Geophysics ◽  
2021 ◽  
pp. 1-145
Author(s):  
Zhiming Ren ◽  
Qianzong Bao ◽  
Bingluo Gu
Keyword(s):  
Wave Velocity ◽  
Velocity Model ◽  
S Wave ◽  
Direct Wave ◽  
Reflected Waves ◽  
Traveltime Inversion ◽  
Long Wavelength ◽  
Velocity Models ◽  
Reflected Wave ◽  
Background Models

Full waveform inversion (FWI) suffers from the local minima problem and requires a sufficiently accurate starting model to converge to the correct solution. Wave-equation traveltime inversion (WETI) is an effective tool to retrieve the long-wavelength components of the velocity model. We develop a joint diving/direct and reflected wave WETI (JDRWETI) method to build the P- and S-wave velocity macromodels. We estimate the traveltime shifts of seismic events (diving/direct waves, PP and PS reflections) through the dynamic warping scheme and construct a misfit function using both the time shifts of diving/direct and reflected waves. We derive the adjoint wave equations and the gradients with respect to the background models based on the joint misfit function. We apply the kernel decomposition scheme to extract the kernel of the diving/direct wave and the tomography kernels of PP and PS reflections. For an explosive source, the kernels of diving/direct wave and PP reflections and the kernel of PS reflections are used to compute the P- and S-wave gradients of the background models, respectively. We implement JDRWETI by a two-stage inversion workflow: first invert the P- and S-wave velocity models using the P-wave gradients and then improve the S-wave velocity model using the S-wave gradients. Numerical tests on synthetic and field datasets reveal that the JDRWETI method successfully recovers the long-wavelength components of P- and S-wave velocity models, which can be used for an initial model for the subsequent elastic FWI. Moreover, the proposed JDRWETI method prevails over the existing reflection WETI method and the cascaded diving/direct and reflected wave WETI method, especially when large velocity errors are present in the shallow part of the starting models. The JDRWETI method with the two-stage inversion workflow can give rise to reasonable inversion results even for the model with different P- and S-wave velocity structures.

Wave-equation reflection traveltime inversion with dynamic warping and full-waveform inversion

Geophysics ◽  
2013 ◽  
Vol 78 (6) ◽  
pp. R223-R233 ◽  
Author(s):  
Yong Ma ◽  
Dave Hale
Keyword(s):  
Wave Equation ◽  
Waveform Inversion ◽  
Velocity Model ◽  
Full Waveform Inversion ◽  
Local Minima ◽  
Low Frequencies ◽  
Traveltime Inversion ◽  
Velocity Models ◽  
Full Waveform ◽  
Dynamic Warping

In reflection seismology, full-waveform inversion (FWI) can generate high-wavenumber subsurface velocity models but often suffers from an objective function with local minima caused mainly by the absence of low frequencies in seismograms. These local minima cause cycle skipping when the low-wavenumber component in the initial velocity model for FWI is far from the true model. To avoid cycle skipping, we discovered a new wave-equation reflection traveltime inversion (WERTI) to update the low-wavenumber component of the velocity model, while using FWI to only update high-wavenumber details of the model. We implemented the low- and high-wavenumber inversions in an alternating way. In WERTI, we used dynamic image warping (DIW) to estimate the time shifts between recorded data and synthetic data. When compared with correlation-based techniques often used in traveltime estimation, DIW can avoid cycle skipping and estimate the time shifts accurately, even when shifts vary rapidly. Hence, by minimizing traveltime shifts estimated by dynamic warping, WERTI reduces errors in reflection traveltime inversion. Then, conventional FWI uses the low-wavenumber component estimated by WERTI as a new initial model and thereby refines the model with high-wavenumber details. The alternating combination of WERTI and FWI mitigates the velocity-depth ambiguity and can recover subsurface velocities using only high-frequency reflection data.

Frequency-domain wave-equation traveltime inversion with a monofrequency component

Geophysics ◽  
2021 ◽  
Vol 86 (6) ◽  
pp. R913-R926
Author(s):  
Jianhua Wang ◽  
Jizhong Yang ◽  
Liangguo Dong ◽  
Yuzhu Liu
Keyword(s):  
Wave Equation ◽  
Frequency Domain ◽  
Time Domain ◽  
Computational Cost ◽  
Velocity Model ◽  
Receiver Side ◽  
Data Set ◽  
Traveltime Inversion ◽  
Velocity Models ◽  
The Time Domain

Wave-equation traveltime inversion (WTI) is a useful tool for background velocity model building. It is generally formulated and implemented in the time domain, in which the gradient is calculated by temporally crosscorrelating the source- and receiver-side wavefields. The time-domain source-side snapshots are either stored in memory or are reconstructed through back propagation. The memory requirements and computational cost of WTI are thus prohibitively expensive, especially for 3D applications. To partially alleviate this problem, we provide an implementation of WTI in the frequency domain with a monofrequency component. Because only one frequency is used, it is affordable to directly store the source- and receiver-side wavefields in memory. There is no need for wavefield reconstruction during gradient calculation. In such a way, we have dramatically reduced the memory requirements and computational cost compared with the traditional time-domain WTI realization. For practical implementation, the frequency-domain wavefield is calculated by time-domain finite-difference forward modeling and is transformed to the frequency domain by an on-the-fly discrete Fourier transform. Numerical examples on a simple lateral periodic velocity model and the Marmousi model demonstrate that our method can obtain accurate background velocity models comparable with those from time-domain WTI and frequency-domain WTI with multiple frequencies. A field data set test indicates that our method obtains a background velocity model that well predicts the seismic wave traveltime.

Automatic traveltime inversion via sparse decomposition of seismic data

Geophysics ◽  
2018 ◽  
Vol 83 (6) ◽  
pp. R659-R668 ◽  
Author(s):  
Bo Feng ◽  
Huazhong Wang ◽  
Ru-Shan Wu
Keyword(s):  
Seismic Data ◽  
Low Frequency ◽  
Velocity Model ◽  
Inversion Method ◽  
Traveltime Inversion ◽  
Sparse Decomposition ◽  
Reflection Seismic ◽  
New Strategy ◽  
Gradient Based ◽  
Intersection Criterion

We have developed an automatic traveltime inversion (ATI) method to estimate the macrovelocity model from reflection seismic data. First, we extract the kinematic information (i.e., source/receiver ray parameters, traveltime, and source/receiver coordinates) of locally coherent events using a sparse-decomposition method. And then we evaluate a new strategy to calculate the reflection traveltime residual based on a ray-intersection criterion, eliminating the influence of seismic amplitude to the estimation of the traveltime residual. The velocity model can be updated iteratively by minimizing the traveltime residual functional with a gradient-based method. To obtain a smooth gradient free of artifacts, we first estimate the high-wavenumber components of the functional gradient with a total variation (TV) regularization method and then subtract it from the full gradient. Because the reflection traveltime residual calculation and velocity update are fully automated procedures, the proposed traveltime inversion method is referred to as ATI. We determine with 2D synthetic and field examples that ATI does not need a good starting model. Furthermore, it requires neither low-frequency seismic data nor long-offset acquisition. Nevertheless, the proposed traveltime residual calculation strategy is only valid for the 2D case, which limits its 3D applicability. We explore a possible solution for 3D extension.

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.

Laplace-domain full-waveform inversion of seismic data lacking low-frequency information

Geophysics ◽  
2012 ◽  
Vol 77 (5) ◽  
pp. R199-R206 ◽  
Author(s):  
Wansoo Ha ◽  
Changsoo Shin
Keyword(s):  
Seismic Data ◽  
Field Data ◽  
Waveform Inversion ◽  
Low Frequency ◽  
Single Frequency ◽  
Frequency Content ◽  
Laplace Domain ◽  
Frequency Information ◽  
Velocity Models ◽  
Gaussian Source

The lack of the low-frequency information in field data prohibits the time- or frequency-domain waveform inversions from recovering large-scale background velocity models. On the other hand, Laplace-domain waveform inversion is less sensitive to the lack of the low frequencies than conventional inversions. In theory, frequency filtering of the seismic signal in the time domain is equivalent to a constant multiplication of the wavefield in the Laplace domain. Because the constant can be retrieved using the source estimation process, the frequency content of the seismic data does not affect the gradient direction of the Laplace-domain waveform inversion. We obtained inversion results of the frequency-filtered field data acquired in the Gulf of Mexico and two synthetic data sets obtained using a first-derivative Gaussian source wavelet and a single-frequency causal sine function. They demonstrated that Laplace-domain inversion yielded consistent results regardless of the frequency content within the seismic data.

Two robust imaging methodologies for challenging environments: Wave-equation dispersion inversion of surface waves and guided waves and supervirtual interferometry + tomography for far-offset refractions

2018 ◽  
Vol 6 (4) ◽  
pp. SM27-SM37 ◽  
Author(s):  
Jing Li ◽  
Kai Lu ◽  
Sherif Hanafy ◽  
Gerard Schuster
Keyword(s):  
Wave Equation ◽  
Velocity Distribution ◽  
Surface Waves ◽  
Guided Waves ◽  
Velocity Model ◽  
Dispersion Curves ◽  
P Wave ◽  
Head Wave ◽  
Near Surface ◽  
Velocity Models

Two robust imaging technologies are reviewed that provide subsurface geologic information in challenging environments. The first one is wave-equation dispersion (WD) inversion of surface waves and guided waves (GW) for the shear-velocity (S-wave) and compressional-velocity (P-wave) models, respectively. The other method is traveltime inversion for the velocity model, in which supervirtual refraction interferometry (SVI) is used to enhance the signal-to-noise ratio of far-offset refractions. We have determined the benefits and liabilities of both methods with synthetic seismograms and field data. The benefits of WD are that (1) there is no layered-medium assumption, as there is in conventional inversion of dispersion curves. This means that 2D or 3D velocity models can be accurately estimated from data recorded by seismic surveys over rugged topography, and (2) WD mostly avoids getting stuck in local minima. The liability is that WD for surface waves is almost as expensive as full-waveform inversion (FWI) and, for Rayleigh waves, only recovers the S-velocity distribution to a depth no deeper than approximately 1/2 to 1/3 wavelength of the lowest-frequency surface wave. The limitation for GW is that, for now, it can estimate the P-velocity model by inverting the dispersion curves from GW propagating in near-surface low-velocity zones. Also, WD often requires user intervention to pick reliable dispersion curves. For SVI, the offset of usable refractions can be more than doubled, so that traveltime tomography can be used to estimate a much deeper model of the P-velocity distribution. This can provide a more effective starting velocity model for FWI. The liability is that SVI assumes head-wave first arrivals, not those from strong diving waves.

scholarly journals Passive seismic event estimation using multiscattering waveform inversion

Geophysics ◽  
2019 ◽  
Vol 84 (3) ◽  
pp. KS59-KS69 ◽  
Author(s):  
Chao Song ◽  
Zedong Wu ◽  
Tariq Alkhalifah
Keyword(s):  
Seismic Data ◽  
Waveform Inversion ◽  
Velocity Model ◽  
Source Image ◽  
Time Inversion ◽  
Multiple Sources ◽  
Origin Time ◽  
Secondary Sources ◽  
Velocity Models ◽  
Passive Seismic

Passive seismic monitoring has become an effective method to understand underground processes. Time-reversal-based methods are often used to locate passive seismic events directly. However, these kinds of methods are strongly dependent on the accuracy of the velocity model. Full-waveform inversion (FWI) has been used on passive seismic data to invert the velocity model and source image, simultaneously. However, waveform inversion of passive seismic data uses mainly the transmission energy, which results in poor illumination and low resolution. We developed a waveform inversion using multiscattered energy for passive seismic to extract more information from the data than conventional FWI. Using transmission wavepath information from single- and double-scattering, computed from a predicted scatterer field acting as secondary sources, our method provides better illumination of the velocity model than conventional FWI. Using a new objective function, we optimized the source image and velocity model, including multiscattered energy, simultaneously. Because we conducted our method in the frequency domain with a complex source function including spatial and wavelet information, we mitigate the uncertainties of the source wavelet and source origin time. Inversion results from the Marmousi model indicate that by taking advantage of multiscattered energy and starting from a reasonably acceptable frequency (a single source at 3 Hz and multiple sources at 5 Hz), our method yields better inverted velocity models and source images compared with conventional FWI.

Sign in / Sign up

Export Citation Format

Share Document