摘要
在传统的有限频层析理论中,采用Born或Rytov近似推导地震波旅行时对模型参数的一阶Fréchet微商,隐含了弱散射近似假设。对于强速度扰动,有限频层析正问题的线性化近似不再成立。近年来发展的广义Rytov近似理论可以较为精确地预测散射波的相位扰动,因此更适用于有限频旅行时层析。为突破弱散射近似的限制,将广义Rytov近似应用于传统的有限频层析理论,导出了基于广义Rytov近似理论的有限频旅行时敏感度核函数。由数值试验可知,对于前向小角度散射,无论速度扰动强度如何,新的旅行时核函数均可较为准确地预测地震波的旅行时扰动。针对旅行时层析反问题,提出了一种基于隐式矩阵向量乘的高斯—牛顿反演算法,具有准二阶收敛速度且无需显式计算和存储Hessian矩阵。复杂模型数据测试结果表明,提出的广义Rytov近似旅行时层析可以建立高精度的近地表速度模型,且收敛速度较快。
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 matrixvector 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 realized in a matrix-free fashion,reducing the computer memory and disk occupancy significantly.Numerical tests demonstrated that the proposed GRAbased traveltime tomography can estimate the nearsurface velocity model with high resolution and at a fast convergence rate.
作者
冯波
WU Ru-Shan
罗飞
许荣伟
王华忠
FENG Bo;WU Ru-Shan;LUO Fei;XU Rongwei;WANG Huazhong(Wave Phenomena Intelligent Inversion Imaging Group(WPI),School of Ocean and Earth Science,Tongji University,Shanghai 200092 China;Institute for Advanced Study,Tongji University,Shanghai 200092 China;Modeling and Imaging Laboratory,University of California,Santa Cruz,California 95060,USA)
出处
《石油地球物理勘探》
EI
CSCD
北大核心
2021年第1期92-99,I0010,I0011,共10页
Oil Geophysical Prospecting
基金
国家重点研发计划“深海关键技术与装备”重点专项“深海油气地震勘探高精度速度建模技术”(2019YFC0312004)
国家自然科学基金项目“强散射介质中的波动方程线性化及速度反演”(42074143)和“针对勘探地震介质及波场特征的高波数参数估计方法研究”(41774126)
国家重点研发计划“变革性技术关键科学问题”重点专项“多元信息联合驱动的地震成像方法研究”(2018YFA0702503)
上海市浦江人才计划(20PJ141300)
南方海洋科学与工程广东省实验室(湛江)项目“特征波域地震成像方法研究”(ZJW-2019-04)
中国石化地球物理重点实验室项目“基于波动方程的高精度近地表速度建模方法研究”(33550006-19-FW0399-0041)联合资助。
关键词
有限频旅行时敏感度核函数
广义Rytov近似
初至层析
高斯—牛顿反演算法
finite-frequency traveltime sensitivity kernel
generalized Rytov approximation
first-arrival traveltime tomography
Gauss-Newton algorithm