[发明专利]一种基于锥面波的三维地震速度反演方法

摘要

本申请公开了基于锥面波的三维地震速度反演方法。本发明首先利用参考震源点与射线参数控制生成锥面波，进而对三维地震数据利用锥面波函数进行编码，得到合成锥面波超级炮记录，常规三维平面波编码需要两个射线参数，而本文方法只需要一个射线参数p，并且将平面编码变成圆锥曲面编码，对同一个射线参数可以选取不同的参考点，这样既可以达到单射线参数反演，也可以压制串扰噪音，从而达到反演不同尺度的速度场的目的，通过调整p的角度和平滑处理完成多尺度反演，从而达到在没有低频的地震数据下完成三维地震速度场的反演，从而提高了地震波形反演背景速度场的精度，从而反演的最终速度场更加精确。

主权项

1.一种基于锥面波的三维地震速度反演方法，其特征在于，包括以下步骤：S1，获取观测数据、初始速度场与震源子波；S2，获取多个射线参数和多个参考点，利用以下步骤进行迭代处理，直至每个射线参数迭代时选取参考点的次数满足预设要求时，输出目标速度场：S201，确定多个射线参数的迭代顺序；S203，基于射线参数的迭代顺序，对于任一确定的射线参数，迭代预设次数生成该确定的射线参数所对应的背景速度场，其中，该背景速度场用于下一射线参数进行迭代的初始速度场；S205，当最后一个射线参数迭代完毕时，输出目标速度场；其中，所述S203中，对于任一已经确定的射线参数，该确定的射线参数所对应的背景速度场由如下方式得到：S2031，基于本次迭代中的射线参数和参考点构造锥面波，并用锥面波对观测数据编码，得到合成锥面波超级炮；其中，锥面波的生成函数为：其中，△t表示时间延迟量，(x0,y0)表示在地震数据观测系统之内随机选取的参考点，x0与y0分别表示三维空间上的x轴(主侧线方向)位置与y轴(联络侧线方向)的位置，p为射线参数，其值越大表示锥面与地面夹角越大，第k炮(x_k,y_k)需要编码的炮的炮点位置，合成锥面波超级炮为：其中：K表示一个合成锥面波超级炮里面含有的总炮数，D_k(t)为第k炮t时刻的炮记录值，当小于0时，取值为0；S2032，利用锥面波对震源子波编码，基于当前的背景速度场进行三维有限差分正演得到正演锥面波超级炮；S2033，利用正演锥面波超级炮与合成锥面波超级炮做差反传，互相关确定反演梯度；S2034，根据反演梯度，利用抛物插值计算最佳步长，对当前的背景速度场迭代；S2035，随机选取另一参考点，重复S2031到S2034，直到同一射线参数下的重复次数达到预设次数，生成该确定的射线参数所对应的背景速度场。