1. 同济大学海洋与地球科学学院, 上海 200092; 2. 同济大学海洋高等研究院, 上海 200092; 3. Modeling and Imaging Laboratory, University of California, Santa Cruz, CA 95060, USA
Wave-equation traveltime tomography using the ge-neralized Rytov approximation
FENG Bo1,2, WU Ru-Shan3, LUO Fei1, XU Rongwei1, WANG Huazhong1
1. Wave Phenomena Intelligent Inversion Imaging Group(WPI), School of Ocean and Earth Science, Tongji University, Shanghai 200092 China; 2. Institute for Advanced Study, Tongji University, Shanghai 200092 China; 3. Modeling and Imaging Laboratory, University of California, Santa Cruz, California 95060, USA
Abstract:The conventional finite-frequency tomography is usually derived from the Born or Rytov approximation which implies weak-scattering assumption. Therefore, the linearized forward problem in the finite-frequency theory is not satisfied for strong velocity perturbations. In the case of forward scattering and small-angle propagation, the generalized Rytov approximation (GRA) method recently developed can achieve improving phase accuracy of forward-scattered wavefield, making it more suitable for traveltime tomography. In this paper, we combine the conventional finite-frequency theory with GRA and propose a GRA-based traveltime sensitivity kernel, which works well regardless of the magnitude of velocity perturbations. Numerical examples show that the traveltime perturbation of forward-scattered waves can be correctly handled by the GRA-based traveltime sensitivity kernel. Then we propose an implicit matrix-vector product strategy which can calculate the Hessian matrix-vector product without explicitly forming the Hessian matrix, making it more attractive for 3-D problems. We solve the traveltime inverse problem with the Gauss-Newton method, where the Hessian matrix-vector product is obtained by the proposed implicit matrix-vector product method. Consequently, the Gauss-Newton method can be rea-lized in a matrix-free fashion, reducing the compu-ter memory and disk occupancy significantly. Numerical tests demonstrated that the proposed GRA-based traveltime tomography can estimate the near-surface velocity model with high resolution and at a fast convergence rate.
Slaney M,Kak A C,Larsen L. Limitations of imaging with first-order diffraction tomography[J]. IEEE Transactions On Microwave Theory and Techniques,1984,32(8):860-873.
[2]
Wu R S,Toksöz M N. Diffraction tomography and multisource holography applied to seismic imaging[J]. Geophysics,1987,52(1):11-25.
[3]
Luo Y,Schuster G T. Wave-equation travel time inversion[J]. Geophysics,1991,56(5):645-653.
[4]
Woodward M J. Wave-equation tomography[J]. Geo-physics,1992,57(1):15-26.
[5]
Lailly P. The seismic inverse problem as a sequence of before stack migrations[C]. Expanded Abstracts of Conference on Inverse Scattering,Theory and Application,1983,206-220.
[6]
Liu Q,Tromp J. Finite-frequency kernels based upon adjoint methods[J]. Bulletin of the Seismological Society of America,2006,96(6):2383-2397.
[7]
Liu Q,Tromp J. Finite-frequency sensitivity kernels for global seismic wave propagation based upon adjoint methods[J]. Geophysical Journal International,2008,174(1):265-286.
[8]
Tape C,Liu Q,Tromp J. Finite frequency tomography using adjoint methods-methodology and examples using membrane surface waves[J]. Geophysical Journal International,2007,168(3):1105-1129.
[9]
Tarantola A. Inversion of seismic reflection data in the acoustic approximation[J]. Geophysics,1984,49(8):1259-1266.
李海山,杨午阳,雍学善. 三维一阶速度-应力声波方程全波形反演[J]. 石油地球物理勘探,2018,53(4):730-736.LI Haishan,YANG Wuyang,YONG Xueshan. Three-dimensional full waveform inversion based on the first-order velocity-stress acoustic wave equation[J]. Oil Geophysical Prospecting,2018,53(4):730-736.
[12]
Marquering H,Dahlen F A,Nolet G. Three-dimensional sensitivity kernels for finite-frequency traveltimes:The banana-doughnut paradox[J]. Geophysical Journal International,1999, 137(3):805-815.
[13]
Dahlen F A,Hung S H,Nolet G. Fréchet kernels for finite-frequency traveltimes,I:Theory[J]. Geophy-sical Journal International,2000,141(1):157-174.
[14]
Hung S H,Dahlen F A,Nolet G. Fréchet kernels for finite-frequency traveltime,Ⅱ:Examples[J]. Geophysical Journal International,2000,141(1):175-203.
[15]
Zhao L,Jordan T H,Chapman C H. Three-dimensional Fréchet differential kernels for seismic delay times[J]. Geophysical Journal International,2000,141(3):558-576.
谢春,刘玉柱,董良国,等. 基于声波方程的有限频伴随状态法初至波旅行时层析[J]. 石油地球物理勘探,2015,50(2):274-282.XIE Chun,LIU Yuzhu,DONG Liangguo,et al. First arrival tomography with finite frequency adjoint-state method based on acoustic wave equation[J]. Oil Geophysical Prospecting,2015,50(2):274-282.
[19]
Wu R S. Wave propagation,scattering and imaging using dual-domain one-way and one-return propagators[J]. Pure and Applied Geophysics,2003,160(3-4):509-539.
[20]
Spetzler J,Snieder R. The effects of small-scale he-terogeneity on the arrival time of waves[J]. Geo-physical Journal International,2001,145(3):786-796.
[21]
Spetzler J,Snieder R. The Fresnel volume and transmitted waves[J]. Geophysics,2004,69(3):653-663.
[22]
Jocker J,Spetzler J,Smeulders D,et al. Validation of first-order diffraction theory for the traveltimes and amplitudes of propagating waves[J]. Geophysics,2006,71(6):T167-T177.
[23]
Xie X B,Yang H. The finite-frequency sensitivity kernel for migration residual moveout and its applications in migration velocity analysis[J]. Geophysics,2008,73(6):S241-S249.
[24]
刘玉柱,董良国,李培明,等. 初至波菲涅尔体地震层析成像[J]. 地球物理学报,2009,52(9):2310-2320.LIU Yuzhu,DONG Liangguo,LI Peiming,et al. Fresnel volume tomography based on the first arrival of the seismic wave[J]. Chinese Journal of Geophy-sics,2009,52(9):2310-2320.
[25]
Xu W,Xie X B,Geng J. Validity of the Rytov approximation in the form of finite-frequency sensitivity kernels[J]. Pure and Applied Geophysics,2015,172(6):1409-1427.
[26]
冯波,罗飞,王华忠. 一阶Rytov近似有限频旅行时层析[J]. 地球物理学报,2019,62(6):2217-2226.FENG Bo,LUO Fei,WANG Huazhong. Wave equation traveltime tomography using Rytov approximation[J]. Chinese Journal of Geophysics,2019,62(6):2217-2226.
[27]
刘定进,胡光辉,蔡杰雄,等. 高斯束层析与全波形反演联合速度建模[J]. 石油地球物理勘探,2019,54(5):1046-1056.LIU Dingjin,HU Guanghui,CAI Jiexiong,et al. A joint velocity model building method with Gaussian beam tomography and full waveform inversion for land seismic data[J]. Oil Geophysical Prospecting,2019,54(5):1046-1056.
[28]
Beydoun W B,Tarantola A. First Born and Rytov approximations:Modeling and inversion conditions in a canonical example[J]. The Journal of the Acoustical Society of America,1988,83(3),1045-1055.
[29]
罗飞,王华忠,冯波,等. 透射波旅行时Beam层析成像方法[J]. 石油物探,2019,58(3):356-370.LUO Fei,WANG Huazhong,FENG Bo,et al. Beam tomography based on transmission traveltime[J]. Geo-physical Prospecting for Petroleum,2019,58(3):356-370.
[30]
Feng B,Wu R S,Wang H. A generalized Rytov approximation for small scattering-angle wave propagation and strong perturbation media[J]. Geophysical Journal International,2019,219(2):968-974.
[31]
Snieder R,Lomax A. Wavefield smoothing and the effect of rough velocity perturbations on arrival times and amplitudes[J]. Geophysical Journal International,1996,125(3):796-812.
[32]
Luo F,Wang H,et al. Automatic traveltime picking using structure-enhanced high-dimensional time-frequency transform[C]. Extended Abstracts of 81st EAGE Conference and Exhibition,2019,1-5.
[33]
Yuan Y O,Simons F J,Tromp J. Double-difference adjoint seismic tomography[J]. Geophysical Journal International,2016,206(3):1599-1618.
[34]
Zhang H,Thurber C H. Double-difference tomography:The method and its application to the Hayward Fault,California[J]. Bulletin of the Seismological Society of America,2003,93(5):1875-1889.