Joint transmission and reflection traveltime tomography using the fast sweeping method and the adjoint-state technique

Huang, Junwei; Bellefleur, Gilles

doi:10.1111/j.1365-246x.2011.05273.x

Cited by 49 publications

(18 citation statements)

References 49 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…(4), one needs to calculate the gradient of this misfit function with respect to the velocity model m , which can be done using the adjoint-state technique (Huang, 2012;Taillandier, 2009;Plessix, 2006). Redistribution subject to SEG license or copyright; see Terms of Use at http://library.seg.org/ where obs t is the observed first-arrival traveltime, pred t is the predicted traveltime based on the current velocity model m .…”

Section: Gradient-based Eikonal Traveltime Tomographymentioning

confidence: 99%

A hybrid tomography method for crosswell seismic inversion

Cao

Liao

2014

Beijing 2014 International Geophysical Conference &Amp; Exposition, Beijing, China, 21-24 April 2014

View full text Add to dashboard Cite

The waveform tomography, or full waveform inversion (FWI), is in general better than the ray tomography, however a reliable initial model is usually required to ensure the success. In this work we designed a cascade-like hybrid tomography technology to solve the crosswell seismic inversion problem. The new method is a combination of several widely used tomography technologies. We start from the Radon-transform based back projection (BP) method, which produces a smooth initial velocity model. Next, the Linear Iterative Reconstruction (LIR) method is adopted to update this initial model, then it is further improved by the nonlinear Gradient-based Eikonal Traveltime Tomography (GETT) method. The velocity model reconstructed from the previous multi ray tomography methods is sufficiently reliable to serve as the initial model for the computational intensive waveform tomography, from which the accurate velocity model is obtained. The numerical example shows that this hybrid tomography method has great potential in reconstructing accurate acoustic velocity or other highresolution reservoir characterization between two wells. It is noticeable that this hybrid method is able to obtain an accurate velocity even when the recorded seismic data is in poor coverage at spatial and the signal-to-noise ratio is low.

show abstract

Section: Gradient-based Eikonal Traveltime Tomographymentioning

confidence: 99%

A hybrid tomography method for crosswell seismic inversion

Cao

Liao

2014

Beijing 2014 International Geophysical Conference &Amp; Exposition, Beijing, China, 21-24 April 2014

View full text Add to dashboard Cite

show abstract

“…; Rawlinson and Sambridge ; Taillandier et al . ; Huang and Bellefleur ). Up to now, two methods are mainly used to compute travel times, namely, ray tracing and finite difference (FD) approximations to the eikonal equation.…”

Section: Introductionmentioning

confidence: 99%

An upwind fast sweeping scheme for calculating seismic wave first‐arrival travel times for models with an irregular free surface

Lan

Chen

2017

Geophysical Prospecting

View full text Add to dashboard Cite

The topography‐dependent eikonal equation formulated in a curvilinear coordinate system has recently been established and revealed as being effective in calculating first‐arrival travel times of seismic waves in an Earth model with an irregular free surface. The Lax–Friedrichs sweeping scheme, widely used in previous studies as for approximating the topography‐dependent eikonal equation viscosity solutions, is more dissipative and needs a much higher number of iterations to converge. Furthermore, the required number of iterations grows with the grid refinement and results in heavy computation in dense grids, which hampers the application of the Lax–Friedrichs sweeping scheme to seismic wave travel‐time calculation and high‐resolution imaging. In this paper, we introduce a new upwind fast sweeping solver by discretising the Legendre transform of the numerical Hamiltonian of the topography‐dependent eikonal equation using an explicit formula. The minimisation related to the Legendre transform in the sweeping scheme is solved analytically, which proved to be much more efficient than the Lax–Friedrichs algorithm in solving the topography‐dependent eikonal equation. Several numerical experiments demonstrate that the new upwind fast sweeping method converges and achieves much better accuracy after a finite number of iterations, independently of the mesh size, which makes it an efficient and robust tool for calculating travel times in the presence of a non‐flat free surface.

show abstract

“…To overcome this high‐frequency limitation, several studies have developed methods, such as fat‐ray tomography, Fresnel‐volume tomography (Cervený and Soares ; Schuster and Quintus‐Bosz ; Vasco, Peterson and Majer ), and wave‐equation‐based tomography methods (Luo and Schuster ; Leung and Qian ; Taillandier et al . ; Huang and Belleeur ; Bretaudeau et al . ).…”

Section: Introductionmentioning

confidence: 99%

Applying joint traveltime and waveform inversion to image the Sichuan Basin, China

Zhang

2017

Near Surface Geophysics

View full text Add to dashboard Cite

We apply a joint first-arrival traveltime and early-arrival waveform inversion method to image complex near-surface structures in the Sichuan Basin, China. The area includes rugged topography and large near-surface velocity variations. Due to the near-surface effects, it is difficult to produce highquality reflection images of the deep subsurface. First-arrival traveltime tomography is often applied for near-surface imaging, but the results may not be sufficiently accurate because of ray assumptions and the limited traveltime information. Waveform inversion should allow complex structures to be resolved; however, it may fall into local minima because of cycle skipping issues and may also produce artefacts in very shallow areas that are associated with rugged topography. Therefore, we combine the advantages of the two methods and mitigate their problems by performing joint inversion of the two types of data. We demonstrate the effectiveness of the joint inversion method using synthetic and real data from Sichuan, China. In the real data example, we compare the velocity models resolved from waveform inversion alone with those resolved from the joint inversion. We calculate long-wavelength static corrections and apply them to the data processing. The common midpoint stacking results show that the joint inversion method produces a more effective statics solution. 1995). FWI can resolve the details of a velocity model because it is based on the waveform information; thus, FWI should be able to image geologically complex areas (Tarantola 1984;Sheng 2006). The FWI solution should also provide a higher resolution velocity image than traveltime tomography (Gauthier, Virieux and Tarantola 1986).However, the objective function of waveform inversion is highly nonlinear and may include many local minima. The inversion may fall into a local minima due to cycle skipping between the predicted and observed data (Virieux and Operto 2009). Several methods have been proposed recently to overcome or mitigate the cycle skipping issue. Ma and Hale (2013) proposed a hybrid waveform inversion method in which wave-equation reflection traveltime inversion was used to update the lowwavenumber component, and FWI was employed to update the high-wavenumber details of the velocity model. Wu, Luo and Wu (2013) developed a waveform envelope inversion method to utilise low-frequency information from the data. The envelope inversion method can avoid the cycle skipping problem and retrieve the long-wavelength background velocity model for FWI. Warner and Guasch (2014) proposed adaptive waveform inversion, which reformulates the objective functions in terms of matching filters. The difference between the observed and predicted data is projected onto the time lag of the filter. Other

show abstract

Joint transmission and reflection traveltime tomography using the fast sweeping method and the adjoint-state technique

Cited by 49 publications

References 49 publications

A hybrid tomography method for crosswell seismic inversion

A hybrid tomography method for crosswell seismic inversion

An upwind fast sweeping scheme for calculating seismic wave first‐arrival travel times for models with an irregular free surface

Applying joint traveltime and waveform inversion to image the Sichuan Basin, China

Contact Info

Product

Resources

About