检索规则说明:AND代表“并且”;OR代表“或者”;NOT代表“不包含”;(注意必须大写,运算符两边需空一格)
检 索 范 例 :范例一: (K=图书馆学 OR K=情报学) AND A=范并思 范例二:J=计算机应用与软件 AND (U=C++ OR U=Basic) NOT M=Visual
作 者:冯波[1,2] WU Ru-Shan[3] 罗飞 许荣伟 王华忠 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)
机构地区:[1]同济大学海洋与地球科学学院,上海200092 [2]同济大学海洋高等研究院,上海200092 [3]Modeling and Imaging Laboratory,University of California,Santa Cruz,CA 95060,USA
出 处:《石油地球物理勘探》2021年第1期92-99,I0010,I0011,共10页Oil Geophysical Prospecting
基 金:国家重点研发计划“深海关键技术与装备”重点专项“深海油气地震勘探高精度速度建模技术”(2019YFC0312004);国家自然科学基金项目“强散射介质中的波动方程线性化及速度反演”(42074143)和“针对勘探地震介质及波场特征的高波数参数估计方法研究”(41774126);国家重点研发计划“变革性技术关键科学问题”重点专项“多元信息联合驱动的地震成像方法研究”(2018YFA0702503);上海市浦江人才计划(20PJ141300);南方海洋科学与工程广东省实验室(湛江)项目“特征波域地震成像方法研究”(ZJW-2019-04);中国石化地球物理重点实验室项目“基于波动方程的高精度近地表速度建模方法研究”(33550006-19-FW0399-0041)联合资助。
摘 要:在传统的有限频层析理论中,采用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.
关 键 词:有限频旅行时敏感度核函数 广义Rytov近似 初至层析 高斯—牛顿反演算法
分 类 号:P631[天文地球—地质矿产勘探]
正在载入数据...
正在载入数据...
正在载入数据...
正在载入数据...
正在载入数据...
正在载入数据...
正在载入数据...
正在链接到云南高校图书馆文献保障联盟下载...
云南高校图书馆联盟文献共享服务平台 版权所有©
您的IP:3.138.60.117