基于Rytov近似的波动方程走时层析方法研究1Wave-equation-based.PDF

1 基于Rytov 近似的波动方程走时层析方法研究 冯波, 罗飞, 王华忠 同济大学海洋与地球科学学院，波现象与反演成像研究组，上海，200092 【摘要】传统的波动方程走时核函数（或走时Fréchet 导数）多基于互相关时差测量方式及地 震波场的一阶Born 近似导出，其成立条件非常苛刻。然而，地震波走时与大尺度的速度结构 具有良好的线性关系，对于小角度的前向散射波场，Rytov 近似优于Born 近似。因此，本文基 于Rytov 近似和互相关时差测量方式，导出了基于Rytov 近似的有限频走时敏感度核函数的两 种等价形式：频率积分和时间积分表达式。在此基础之上，本文提出了一种隐式矩阵向量乘方 法，可以直接计算 Hessian 矩阵或者核函数与向量的乘积，而无需显式计算和存储核函数及 Hessian 矩阵。基于隐式矩阵向量乘方法，本文利用共轭梯度法求解法方程实现了一种高效的 Gauss-Newton 反演算法求解走时层析反问题。与传统的敏感度核函数反演方法相比，本文方法 在每次迭代过程中，无需显式计算和存储核函数，极大降低了存储需求。与基于梯度计算的伴 随状态方法相比，本文方法具有准二阶的收敛速度，收敛速度更快。数值实验证明了本文方法 的有效性。 【关键词】Rytov 近似，有限频走时敏感度核函数，波动方程走时层析，初至波，隐式矩阵向 量乘，Gauss-Newton 方法 Wave-equation-based traveltime tomography using Rytov approximation FENG Bo, LUO Fei, WANG Hua-Zhong Wave phenomena and inversion imaging group (WPI), School of Ocean and Earth Science, Tongji University, Shanghai, 200092, China Abstract: The conventional wave-equation traveltime sensitivity kernel (TSK) or traveltime Fréchet derivative is derived from the Born approximation and cross-correlation measurement, which has a very narrow valid condition. In fact, the seismic traveltime has a more linear relationship with the large-scale velocity structure. For small-angle forward scattered wavefield, Rytov approximation is proved to be superior to Born approximation. Based on the Rytov approximation and cross-correlation measurement, a new wave-equation traveltime sensitivity kernel is derived. Meanwhile, an implicit matrix-vector product method is proposed, which can directly calculate the product of a matrix (TSK) and a model-space vector as well as the product of a matrix transpose and a data-space vector, eliminating the need of calculatin