一种自适应波动方程波场延拓静校正方法

摘要

本发明公开了一种自适应波动方程波场延拓静校正方法。使用层析反演获得的速度场信息，利用原始单炮地震数据进行波动方程延拓，将地表采集的地震数据延拓到某一个基准面上。从共炮集记录开始使用有限差分法进行延拓，从地表延拓到高速层顶界面，然后采用相移法从高速层顶界面延拓到固定基准面上。将延拓后的数据进行道集抽取，抽取成共检波点道集数据，然后再采用有限差分法和相移法对共检波点道集记录进行波场延拓。完成共检波点数据延拓后将数据重新抽取为共炮道集数据。在进行波场延拓的过程中，加入随速度场复杂度自动调整的速度场平滑参数，使用经过平滑处理后的速度场进行延拓，可以解决有限差分法在大倾角下准确度低的问题。

主权项

一种自适应波动方程波场延拓静校正方法，其特征在于：包括以下步骤：步骤1：采用原始的单炮地震数据；步骤2：使用常规的初至波层析方法获得近地表的速度结构；步骤3：采用有限差分法将地表获得的波场延拓到高速层顶界面常规有限差分法延拓公式为：$\frac{\partial p}{\partial z}=\frac{\mathrm{i\omega }}{v(x,z)}·\frac{a_j\frac{v^2(x,z)}{{\omega }^2}\frac{{\partial }^2}{{\partial x}^2}}{1+b_j\frac{v^2(x,z)}{{\omega }^2}\frac{{\partial }^2}{\partial x^2}}p$   (式1)式中：p为波场值，v为速度，a_j、b_j为系数，j＝1，2，…，n；改进后的有限差分法延拓方程为：$\frac{\partial p}{\partial z}=\frac{\mathrm{i\omega }}{\bar{v}(x,z)}·\frac{a_j\frac{{\bar{v}}^2(x,z)}{{\omega }^2}\frac{{\partial }^2}{{\partial x}^2}}{1+b_j\frac{{\bar{v}}^2(x,z)}{{\omega }^2}\frac{{\partial }^2}{{\partial x}^2}}p$   (式2)式中：为模型自适应平滑速度，在进行计算前，需用对当前选取速度点三维空间一定范围内速度进行分析，获得其横向变化系数r＝|(v₀‑v₁)/v₀|，当r<0.2时，不修改原速度场，当0.2≤r时，按照w＝5*(r‑0.2)计算平滑范围；采用改进后的延拓方程，在延拓的过程中判断周围速度场的变化程度，自动进行平滑处理，使有限差分法在遇到速度突变时造成的空间假频现象得到压制；步骤4：采用相移法将波场从高速层顶界面延拓到地表，对数据进行延拓时采用变步长延拓法进行延拓；步骤5：对延拓后的数据进行分选，按照检波点选排为共检波点数据集；步骤6：再采用有限差分法和相移法对共检波点道集记录进行波场延拓；处理完成后，重新选排为共炮道集数据。