3D first-arrival traveltime tomography with linearized eikonal equation and shaping regularization
Yang Guoquan1,2, Ding Pengcheng1,2, Li Zhenchun1,2, Zhang Kai1,2, Wang Chong1,2
1. School of Geosciences, China University of Petroleum(East China), Qingdao, Shandong 266580, China;
2. Function Laboratory of Marine Geo-resource Evaluation and Exploration Technology, National Marine Laboratory, Qingdao, Shandong 266071, China
Abstract:First-arrival traveltime tomography is an important tool for estimating near-surface long-wavenumber seismic velocities.However,conventional algorithms utilized for traveltime tomography usually encounter problems when handling large 3D seismic data.For instance the raytracing has low computational efficiency and instability when strong velocity contrast exists in velocity models,and conventional regularization algorithms lie in their low iterative convergence rate caused by the insufficient control on models estimated by Tikhonov's regularization,which could be detrimental for large 3D scale problems,when only very few iterations are affordable.Therefore,on the basis of linearized eikonal equation,we propose a new algorithm to calculate the Fréchet derivative matrix with the upwind differences scheme without using cumbersome raytracing.Additionally,we incorporate a shaping regularization into conjugate gradient algorithm to iteratively minimize the data misfit.Finally,numerical experiments on 3D models demonstrate that the proposed algorithm has higher accuracy and faster convergence compared with the conventional algorithms.
Sang Yunyun,Li Zhenchun,Zhang Kai. Shortest path ray tracing based on parabolic travel-time interpolation. OGP,2013,48(3):403-409.
[3]
Zhu X,Sixta D P,Angstman B G. Tomostatics:turning-ray tomography+static corrections. The Leading Edge,1992,11(12):15-23.
[4]
Zhang J,Toksöz M N. Nonlinear refraction traveltime tomography. Geophysics,1998,63(5):1726-1737.
[5]
Zelt C A,Barton P J. Three-dimensional seismic refraction tomography:A comparison of two methods applied to data from the Faeroe Basin.Journal of Geo-physical Research,1998,103(B4):7187-7210.
[6]
Berryman J. Stable iterative reconstruction algorithm for nonlinear traveltime tomography. Inverse Problems,1990,6(1):21-42.
Liu Yuzhu,Dong Liangguo,Li Peiming et al. Fresnel volume tomography based on the first arrival of the seismic wave. Chinese Journal of Geophysics,2009,52(9):2310-2320.
Fomel S. Traveltime computation with the linearized eikonal equation. Geophysical Prospecting,2002,50(4):373-382.
[10]
Franklin J B,Harris J M. A high-order fast marching scheme for the linearized eikonal equation. Journal of Computational Acoustics,2001,9(3):1095-1109.
[11]
Taillandier C,Noble M,Chauris H et al. First-arrival traveltime tomography based on the adjoint-state method. Geophysics,2009,74(6):WCB57-WCB66.
[12]
Leung S,Qian J. An adjoint state method for three-dimensional transmission traveltime tomography using first arrivals. Communications in Mathematical Sciences,2006,4(1):249-266.
Cui Yan,Wang Yanfei. Tikhonov regularization and gradient descent algorithms for tomography using first-arrival seismic traveltimes. Chinese Journal of Geophysics,2015,58(4):1367-1377.
[19]
Fomel S. Shaping regularization in geophysical estimation problems. Geophysics,2007,72(2):R29-R36.
[20]
Liu Y,Yang J,Chi B et al. An improved scattering-integral approach for frequency-domain full waveform inversion. Geophysical Journal International,2015,202(3):1827-1842.
Li Yongde,Dong Liangguo,Liu Yuzhu. First arrival traveltime tomography based on an improved scattering-integral algorithm. Chinese Journal of Geophy-sics,2016,59(10):3820-3828.
Wang Fei,Qu Xinxin,Liu Sixin et al. A new ray tracing approach based on both multi-stencils fast marching and the steepest descent. OGP,2014,49(6):1109-1114.