[发明专利]初至波波形反演近地表速度模型方法

说明书

技术领域

本发明属于石油地震勘探地震资料处理技术领域，旨在提高地震资料处理中的近地表速度模型精度，进而提高静校正精度和浅层深度成像精度，具体地说是一种基于初至波波形（有别于传统的初至走时）反演，建立高精度近地表速度模型的新技术、新方法。

背景技术

目前反演近地表速度模型，绝大部分是基于射线理论，利用初至波旅行时间反演近地表速度模型，如折射静校正、层析反演等等，由于射线理论是对地震波场的高频近似，只能反演出近地表速度模型的低波数成分（模型尺度远大于地震波长），反演的速度模型分辨率受到限制。基于波动理论的全波形反演方法，利用全波场信息，反演地下速度模型，虽理论比较完善，但受限于数据的信噪比、反演方法本身的多解性和计算效率等问题，实际应用受到很大限制。

为了提高近地表速度模型的反演精度，借鉴全波形反演方法的技术路线，研究了利用初至波波形（旅行时间、振幅）信息反演近地表速度模型方法。

发明内容

鉴于已有技术的不足，本发明的目的是提供一种基于初至波波形反演，建立高精度近地表速度模型的新技术、新方法；该方法基于波动方程理论，和传统的基于射线理论的走时层析反演相比，能够反演出近地表速度模型的高波数成分（模型尺度接近或小于地震波长），提高了近地表速度模型的分辨率；和全波形反演方法相比，该方法只利用能量较强、波形比较稳定的初至波，降低了全波形反演方法的多解性，提高了计算精度和计算效率；该方法具有提高静校正精度和浅层深度成像精度的特点。

实现解决技术问题的技术方案是：初至波波形反演近地表速度模型方法，该方法包括声波方程交错网格有限差分正演模拟技术和最速下降法波形反演技术，其实现过程的方法步骤如下：

（1）、提取时间域初至波波形记录P_obs和初始速度模型；

（2）、利用声波方程交错网格有限差分正演模拟计算模拟波场P_cal，并计算波场残差δp＝P_obs-P_cal；

（3）、将波场残差反向传播得到回传波场P′；

（4）、利用回传波场和正向传播波场计算目标函数的梯度，即速度的更新方向，并计算更新步长；

（5）、对速度模型进行更新；

（6）、检验是否满足迭代终止条件，满足则输出速度模型的反演结果，否则返回步骤（2），继续迭代更新。

所述声波方程交错网格有限差分正演模拟技术，利用一阶声波方程有限差分方法，数值模拟炮点激发的正向传播波场；非均匀介质一阶声波方程：应力-质点速度方程；时间2阶、空间12阶的高阶有限差分，PML边界条件，交错网格。

所述最速下降法波形反演技术，⑴地面观测记录加时窗函数—初至波窗函数,作为一阶声波方程的初值条件，类似于炮点激发正向传播波场正演技术，数值模拟初至波反向传播的地震波场和数值模拟的正向传播波场求残差，利用L2模建立反演目标函数；⑵利用正向传播波场和波场残差函数，计算目标函数的最速下降方向即梯度，确定出速度扰动的最大变化方向；⑶利用梯度和Frechet导数计算出速度扰动的迭代步长；⑷利用速度更新方向和更新步长更新速度。

本发明相比现有技术的显著效果是：

1）本发明基于波动方程理论而没有高频近似假设，能够反演出近地表速度模型的高波数成分（模型尺度接近或小于地震波长），提高了速度模型的分辨率，进而提高静校正的精度，并为近地表和地下深度偏移奠定了良好的基础；

2）本发明借鉴了基于波动方程理论的全波形反演技术，利用了能量较强、波形比较稳定的初至波进行反演，降低了全波形反演的多解性，提高了反演的稳定性和计算效率；

3）本发明采用时间2阶、空间12阶的高精度交错网格有限差分正演模拟方法，并利用PML（完全匹配层）边界条件消除边界效应，提高了正演模拟精度，有效地减小了人工边界反射对有效波场的影响；

4）本发明采用声波方程正演模拟技术，和传统的标量方程规则网格有限差分相比，该技术可以降低频散，减少边界效应，从而提高数值模拟的精度和稳定性；

5）本发明采用最速下降法速度反演技术，并通过正演和反演的多次迭代，直至目标函数完全收敛，就完成了初至波波形反演近地表速度模型；为了提高反演的效率和精度，初始模型可采用常规射线理论层析反演的近地表模型；

6）本发明采用波形反演的同时利用振幅和走时信息，比传统基于射线理论的旅行时层析反演方法具有更高的反演精度，能够获取较为精确的境地表速度模型，为静校正和浅层深度域成像精度创造了条件。