[发明专利]基于散射积分法的非线性菲涅尔体地震走时层析成像方法

摘要

本发明涉及一种基于散射积分法的非线性菲涅尔体地震走时层析成像方法，包括以下步骤：1)对原始地震数据进行去噪和反褶积的预处理；2)在预处理过的地震数据炮记录上进行初至拾取，获得拾取走时；3)根据地下先验信息建立初始表层速度模型，并设定反演参数，包括反演频率范围和间隔；4)根据初至波到达时和初始表层速度模型进行散射积分法非线性菲涅尔体地震层析成像反演，获得最终表层速度模型并成像。与现有技术相比，本发明具有占用的内存小、计算效率高、易于并行、方便预条件、反演精度高、过程稳定等优点。

主权项

1.一种基于散射积分法的非线性菲涅尔体地震走时层析成像方法，其特征在于，包括以下步骤：/n1)对原始地震数据进行去噪和反褶积的预处理；/n2)在预处理过的地震数据炮记录上进行初至拾取，获得拾取走时；/n3)根据地下先验信息建立初始表层速度模型，并设定反演参数，包括反演频率范围和间隔；/n4)根据初至波到达时和初始表层速度模型进行散射积分法非线性菲涅尔体地震层析成像反演，获得最终表层速度模型并成像，具体包括以下步骤：/n41)在当前表层速度模型上进行射线追踪获取理论合成初至走时以及模型空间的旅行时场；/n42)判断理论合成走时与拾取走时是否匹配，若是，则将当前表层速度模型作为最终表层速度模型，若否，则执行步骤43)；/n43)采用LU分解算法获取所有频率的炮点端波场与检波点端格林函数；/n44)根据炮点端波场和格林函数计算每一个炮检对所对应的核函数，并根据旅行时场确定菲涅尔体的范围，获得带核函数特征的菲涅尔体；/n45)通过核函数-标量乘获得每一个炮检对所对应的梯度，累加形成整个观测系统所对应的梯度，同时计算核函数自身元素的模平方并沿着所有炮检对累加求和获得梯度的预条件算子；/n46)计算预条件最速下降方向、方向更新步长和模型更新量后更新当前表层速度模型；/n47)重复步骤41-46)，直到获得匹配的最终表层速度模型。/n