scholarly journals Machine Learning Enabled Traveltime Inversion Based on the Horizontal Source Location Perturbation

Geophysics ◽  
2021 ◽  
pp. 1-40
Author(s):  
Isa Eren Yildirim ◽  
Tariq Alkhalifah ◽  
Ertugrul Umut Yildirim

Gradient based traveltime tomography, which aims to minimize the difference between modeled and observed first arrival times, is a highly non-linear optimization problem. Stabilization of this inverse problem often requires employing regularization. While regularization helps avoid local minima solutions, it might cause low resolution tomograms because of its inherent smoothing property. On the other hand, although conventional ray-based tomography can be robust in terms of the uniqueness of the solution, it suffers from the limitations inherent in ray tracing, which limits its use in complex media. To mitigate the aforementioned drawbacks of gradient and ray-based tomography, we approach the problem in a completely novel way leveraging data-driven inversion techniques based on training deep convolutional neural networks (DCNN). Since DCNN often face challenges in detecting high level features from the relatively smooth traveltime data, we use this type of network to map horizontal changes in observed first arrival traveltimes caused by a source shift to lateral velocity variations. The relationship between them is explained by a linearized eikonal equation. Construction of the velocity models from this predicted lateral variation requires information from, for example, a vertical well-log in the area. This vertical profile is then used to build a tomogram from the output of the network. Both synthetic and field data results verify that the suggested approach estimates the velocity models reliably. Because of the limited depth penetration of first arrival traveltimes, the method is particularly favorable for near-surface applications.

2013 ◽  
Vol 18 (3) ◽  
pp. 183-194 ◽  
Author(s):  
C. A. Zelt ◽  
S. Haines ◽  
M. H. Powers ◽  
J. Sheehan ◽  
S. Rohdewald ◽  
...  

2021 ◽  
Author(s):  
Siegfried Rohdewald

<p>We demonstrate improved resolution in P-wave velocity tomograms obtained by inversion of the synthetic SAGEEP 2011 refraction traveltime data (Zelt 2010) using Wavepath-Eikonal Traveltime Inversion (WET; Schuster 1993) and Wavelength-Dependent Velocity Smoothing (WDVS; Zelt and Chen 2016). We use a multiscale inversion approach and a Conjugate-Gradient based search method. Our default starting model is a 1D-gradient model obtained directly from the traveltime first arrivals assuming diving waves (Sheehan, 2005). As a second approach, we map the first breaks to assumed refractors and obtain a layered starting model using the Plus-Minus refraction method (Hagedoorn, 1959). We compare tomograms obtained using WDVS to smooth the current velocity model grid before forward modeling traveltimes vs. tomograms obtained without WDVS. Results show that WET images velocity layer boundaries more sharply when engaging WDVS. We determine the optimum WDVS frequency iteratively by trial-and-error. We observe that the lower the used WDVS frequency, the stronger the imaged velocity contrast at the top-of-basement. Using a WDVS frequency that is too low makes WDVS based WET inversion unstable exhibiting increasing RMS error, too high modeled velocity contrast and too shallow imaged top-of-basement. To speed up WDVS, we regard each nth node only when scanning the velocity along straight scan lines radiating from the current velocity grid node. Scanned velocities are weighted with a Cosine-Squared function as described by (Zelt and Chen, 2016). We observe that activating WDVS allows decreasing WET regularization (smoothing and damping) to a higher degree than without WDVS.</p><p>References:</p><p><span>Hagedoorn, J.G., 1959, </span><span>The Plus-Minus method of interpreting seismic refraction sections, Geophysical Prospecting</span><span>, Volume 7, 158-182.</span></p><p><span>Rohdewald, S.R.C., 2021, SAGEEP11 data interpretation, https://rayfract.com/tutorials/sageep11_16.pdf.</span></p><p>Schuster, G.T., Quintus-Bosz, A., 1993, <span>Wavepath eikonal traveltime inversion: Theory</span>. Geophysics, Volume 58, 1314-1323.</p><p><span>Sheehan, J.R., Doll, W.E., Mandell, W., 2005, </span><span>An evaluation of methods and available software for seismic refraction tomography analysis</span><span>, JEEG, Volume 10(1), 21-34.</span></p><p>Shewchuk, J.R., 1994, An Introduction to the Conjugate Gradient Method Without the Agonizing Pain, <span>http://www.cs.cmu.edu/~quake-papers/painless-conjugate-gradient.pdf</span><span>. </span></p><p>Zelt, C.A., 2010, Seismic refraction shootout: blind test of methods for obtaining velocity models from first-arrival travel times, <span>http://terra.rice.edu/department/faculty/zelt/sageep2011</span>.</p><p><span>Zelt, C.A., Haines, S., Powers, M.H. et al. 2013, </span><span>Blind Test of Methods for Obtaining 2-D Near-Surface Seismic Velocity Models from First-Arrival Traveltimes</span><span>, JEEG, Volume 18(3), 183-194. </span></p><p><span>Zelt, C.A., Chen, J., 2016, </span><span>Frequency-dependent traveltime tomography for near-surface seismic refraction data</span><span>, Geophys. J. Int., Volume 207, 72-88. </span></p>


Geophysics ◽  
1992 ◽  
Vol 57 (2) ◽  
pp. 353-362 ◽  
Author(s):  
Livia J. Squires ◽  
Samuel N. Blakeslee ◽  
Paul L. Stoffa

Seismic first arrival times from crosshole, VSP, and reversed VSP (RVSP) experiments are collectively inverted by least‐squares for the velocity distribution between two boreholes. The tomographic reconstruction exhibits a large lateral velocity contrast that is not supported by the surface reflection data from the same location. After examining the traveltime residuals from the three tomographic datasets separately, we conclude that the velocity contrast is due primarily to static delays in the RVSP first arrival times. When a first‐order correction is made for the statics, tomographic inversion results in a velocity reconstruction that is more consistent with the surface reflection data. To isolate the velocity errors produced by the RVSP statics, we compute a residual tomogram by subtracting the statics adjusted tomogram from the original. The residual tomogram shows that the statics introduce errors not only in the region sampled by the RVSP rays, but they indirectly contaminate other regions of the tomogram as well. We reproduce this velocity error distribution as part of a model study designed to simulate the effects of statics on tomographic velocity reconstructions. Results indicate that traveltime errors on the order of 2 percent can result in tomographic velocity errors of up to 7 percent.


Geophysics ◽  
2019 ◽  
Vol 84 (6) ◽  
pp. U45-U57 ◽  
Author(s):  
Lianlian Hu ◽  
Xiaodong Zheng ◽  
Yanting Duan ◽  
Xinfei Yan ◽  
Ying Hu ◽  
...  

In exploration geophysics, the first arrivals on data acquired under complicated near-surface conditions are often characterized by significant static corrections, weak energy, low signal-to-noise ratio, and dramatic phase change, and they are difficult to pick accurately with traditional automatic procedures. We have approached this problem by using a U-shaped fully convolutional network (U-net) to first-arrival picking, which is formulated as a binary segmentation problem. U-net has the ability to recognize inherent patterns of the first arrivals by combining attributes of arrivals in space and time on data of varying quality. An effective workflow based on U-net is presented for fast and accurate picking. A set of seismic waveform data and their corresponding first-arrival times are used to train the network in a supervised learning approach, then the trained model is used to detect the first arrivals for other seismic data. Our method is applied on one synthetic data set and three field data sets of low quality to identify the first arrivals. Results indicate that U-net only needs a few annotated samples for learning and is able to efficiently detect first-arrival times with high precision on complicated seismic data from a large survey. With the increasing training data of various first arrivals, a trained U-net has the potential to directly identify the first arrivals on new seismic data.


Geophysics ◽  
2018 ◽  
Vol 83 (6) ◽  
pp. U63-U77
Author(s):  
Bernard K. Law ◽  
Daniel Trad

An accurate near-surface velocity model is critical for weathering statics correction and initial model building for depth migration and full-waveform inversion. However, near-surface models from refraction inversion often suffer from errors in refraction data, insufficient sampling, and over-simplified assumptions used in refraction algorithms. Errors in refraction data can be caused by picking errors resulting from surface noise, attenuation, and dispersion of the first-arrival energy with offset. These errors are partially compensated later in the data flow by reflection residual statics. Therefore, surface-consistent residual statics contain information that can be used to improve the near-surface velocity model. We have developed a new dataflow to automatically include median and long-wavelength components of surface-consistent reflection residual statics. This technique can work with any model-based refraction solution, including grid-based tomography methods and layer-based methods. We modify the cost function of the refraction inversion by adding model and data weights computed from the smoothed surface-consistent residual statics. By using an iterative inversion, these weights allow us to update the near-surface velocity model and to reject first-arrival picks that do not fit the updated model. In this nonlinear optimization workflow, the refraction model is derived from maximizing the coherence of the reflection energy and minimizing the misfit between model arrival times and the recorded first-arrival times. This approach can alleviate inherent limitations in shallow refraction data by using coherent reflection data.


2018 ◽  
Vol 16 (5) ◽  
pp. 507-526 ◽  
Author(s):  
Amin Khalaf ◽  
Christian Camerlynck ◽  
Nicolas Florsch ◽  
Ana Schneider

Geophysics ◽  
2010 ◽  
Vol 75 (6) ◽  
pp. U39-U47 ◽  
Author(s):  
Hui Liu ◽  
Hua-wei Zhou ◽  
Wenge Liu ◽  
Peiming Li ◽  
Zhihui Zou

First-arrival traveltime tomography is a popular approach to building the near-surface velocity models for oil and gas exploration, mining, geoengineering, and environmental studies. However, the presence of velocity-inversion interfaces (VIIs), across which the overlying velocity is higher than the underlying velocity, might corrupt the tomographic solutions. This is because most first-arrival raypaths will not traverse along any VII, such as the top of a low-velocity zone. We have examined the impact of VIIs on first-arrival tomographic velocity model building of the near surface using a synthetic near-surface velocity model. This examination confirms the severe impact of VIIs on first-arrival tomography. When the source-to-receiver offset is greater than the lateral extent of the VIIs, good near-surface velocity models can still be established using a multiscale deformable-layer tomography (DLT), which uses a layer-based model parameterization and a multiscale scheme as regularization. Compared with the results from a commercial grid-based tomography, the DLT delivers much better near-surface statics solutions and less error in the images of deep reflectors.


2021 ◽  
Vol 225 (2) ◽  
pp. 1020-1031
Author(s):  
Huachen Yang ◽  
Jianzhong Zhang ◽  
Kai Ren ◽  
Changbo Wang

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.


Sign in / Sign up

Export Citation Format

Share Document