初至波波形反演近地表速度模型方法
【专利摘要】本发明公开了初至波波形反演近地表速度模型方法，该方法包括声波方程波场正演和最速下降法波形反演技术，其方法步骤为：①提取时间域初至波波形记录和初始模型；②利用声波方程交错网格有限差分正演模拟计算模拟波场及波场残差；③将波场残差反向传播得到回传波场；④利用回传波场和正向传播波场计算目标函数的梯度，并计算更新步长；⑤对速度模型进行更新；⑥检验是否满足迭代终止条件，满足则输出速度模型，否则返回②步骤，继续迭代更新；本发明借鉴了基于波动方程理论的全波形反演技术，利用了能量较强、波形比较稳定的初至波进行反演，降低了全波形反演的多解性，提高了反演的稳定性和计算效率；提高了静校正和浅层深度成像的精度。
【专利说明】初至波波形反演近地表速度模型方法
【技术领域】
[0001]本发明属于石油地震勘探地震资料处理【技术领域】，旨在提高地震资料处理中的近地表速度模型精度，进而提高静校正精度和浅层深度成像精度，具体地说是一种基于初至波波形(有别于传统的初至走时)反演，建立高精度近地表速度模型的新技术、新方法。
【背景技术】
[0002]目前反演近地表速度模型，绝大部分是基于射线理论，利用初至波旅行时间反演近地表速度模型，如折射静校正、层析反演等等，由于射线理论是对地震波场的高频近似，只能反演出近地表速度模型的低波数成分(模型尺度远大于地震波长)，反演的速度模型分辨率受到限制。基于波动理论的全波形反演方法，利用全波场信息，反演地下速度模型，虽理论比较完善，但受限于数据的信噪比、反演方法本身的多解性和计算效率等问题，实际应用受到很大限制。
[0003]为了提高近地表速度模型的反演精度，借鉴全波形反演方法的技术路线，研究了利用初至波波形(旅行时间、振幅)信息反演近地表速度模型方法。
【发明内容】
[0004]鉴于已有技术的不足，本发明的目的是提供一种基于初至波波形反演，建立高精度近地表速度模型的新技术、新方法；该方法基于波动方程理论，和传统的基于射线理论的走时层析反演相比，能够反演出近地表速度模型的高波数成分(模型尺度接近或小于地震波长)，提高了近地表速度模型的分辨率；和全波形反演方法相比，该方法只利用能量较强、波形比较稳定的初至波，降低了全波形反演方法的多解性，提高了计算精度和计算效率；该方法具有提高静校正精度和浅层深度成像精度的特点。
[0005]实现解决技术问题的技术方案是:初至波波形反演近地表速度模型方法，该方法包括声波方程交错网格有限差分正演模拟技术和最速下降法波形反演技术，其实现过程的方法步骤如下:
[0006]( I )、提取时间域初至波波形记录Ptjbs和初始速度模型；
[0007](2)、利用声波方程交错网格有限差分正演模拟计算模拟波场PMl，并计算波场残差 5 P = Pobs-Pcal ;
[0008](3)、将波场残差反向传播得到回传波场P'；
[0009](4)、利用回传波场和正向传播波场计算目标函数的梯度，即速度的更新方向，并计算更新步长；
[0010](5)、对速度模型进行更新；
[0011](6)、检验是否满足迭代终止条件，满足则输出速度模型的反演结果，否则返回步骤(2)，继续迭代更新。
[0012]所述声波方程交错网格有限差分正演模拟技术，利用一阶声波方程有限差分方法，数值模拟炮点激发的正向传播波场；非均匀介质一阶声波方程:应力-质点速度方程；时间2阶、空间12阶的高阶有限差分，PML边界条件，交错网格。
[0013]所述最速下降法波形反演技术，⑴地面观测记录加时窗函数一初至波窗函数，作为一阶声波方程的初值条件，类似于炮点激发正向传播波场正演技术，数值模拟初至波反向传播的地震波场和数值模拟的正向传播波场求残差，利用L2模建立反演目标函数；(2)利用正向传播波场和波场残差函数，计算目标函数的最速下降方向即梯度，确定出速度扰动的最大变化方向利用梯度和Frechet导数计算出速度扰动的迭代步长利用速度更新方向和更新步长更新速度。
[0014]本发明相比现有技术的显著效果是:
[0015]I)本发明基于波动方程理论而没有高频近似假设，能够反演出近地表速度模型的高波数成分(模型尺度接近或小于地震波长)，提高了速度模型的分辨率，进而提高静校正的精度，并为近地表和地下深度偏移奠定了良好的基础；
[0016]2)本发明借鉴了基于波动方程理论的全波形反演技术，利用了能量较强、波形比较稳定的初至波进行反演，降低了全波形反演的多解性，提高了反演的稳定性和计算效率；
[0017]3)本发明采用时间2阶、空间12阶的高精度交错网格有限差分正演模拟方法，并利用PML (完全匹配层)边界条件消除边界效应，提高了正演模拟精度，有效地减小了人工边界反射对有效波场的影响；
[0018]4)本发明采用声波方程正演模拟技术，和传统的标量方程规则网格有限差分相t匕，该技术可以降低频散，减少边界效应，从而提高数值模拟的精度和稳定性；
[0019]5)本发明采用最速下降法速度反演技术，并通过正演和反演的多次迭代，直至目标函数完全收敛，就完成了初至波波形反演近地表速度模型；为了提高反演的效率和精度，初始模型可采用常规射线理论层析反演的近地表模型；
[0020]6)本发明采用波形反演的同时利用振幅和走时信息，比传统基于射线理论的旅行时层析反演方法具有更高的反演精度，能够获取较为精确的境地表速度模型，为静校正和浅层深度域成像精度创造了条件。
【专利附图】
【附图说明】
[0021]图1是本发明近地表速度异常体模型的示意图
[0022]图2是本发明声波方程有限差分正演模拟结果示意图
[0023]图2a为规则网格有限差分正演模拟得到的单炮波形记录图
[0024]图2c为规则网格有限差分正演模拟得到的200ms波场快照图
[0025]图2b为本方法正演模拟得到的单炮波形记录示意图
[0026]图2d为本方法正演模拟得到的200ms波场快照示意图
[0027]图3是本发明实际速度模型示意图
[0028]图4是本发明传统旅行时层析反演结果示意图
[0029]图5是本发明初至波波形反演结果示意图
[0030]图5a为第I次迭代反演结果示意图
[0031]图5b为第100次迭代反演结果示意图
[0032]图5c为第500次迭代反演结果示意图[0033]图6是本发明初至波波形反演目标函数梯度示意图
[0034]上述图1-图6均为计算机绘画软件suffer绘制的附图
【具体实施方式】
[0035]下面结合附图和实例对该技术方案作进一步描述
[0036]初至波波形反演近地表速度模型方法，该方法包括声波方程交错网格有限差分正演模拟和最速下降法波形反演两大技术，其实现过程的方法步骤如下:
[0037]⑴提取时间域初至波波形记录Pobs和初始速度模型；
[0038]⑵利用声波方程交错网格有限差分正演模拟计算模拟波场Pcal,并计算波场残差
5P = Pobs-Pcal ；
[0039]⑶将波场残差反向传播得到回传波场P'；
[0040]⑷利用回传波场和正向传播波场计算目标函数的梯度，即速度的更新方向，并计算更新步长；
[0041 ] (5)对速度模型进行更新；
[0042](6)检验是否满足迭代终止条件，满足则输出速度模型的反演结果，否则返回步骤2，继续迭代更新。
[0043]所述1、声波方程交错网格有限差分正演模拟技术
[0044]利用一阶声波方程有限差分方法，数值模拟炮点激发的正向传播波场；
[0045]技术要点:非均匀介质一阶声波方程(应力-质点速度方程，公式(I));高阶有限差分(时间2阶、空间12阶)，PML边界条件，交错网格；与传统的标量方程规则网格有限差分相比，该技术可以降低频散，减少边界效应，从而提高数值模拟的精度和稳定性。
[0046]
【权利要求】
1.一种初至波波形反演近地表速度模型方法，该方法包括声波方程交错网格有限差分正演模拟和最速下降法波形反演技术，其特征在于:实现过程的方法步骤如下: (I )、提取时间域初至波波形记录Ptjbs和初始速度模型； (2)、利用声波方程交错网格有限差分正演模拟计算模拟波场PMl，并计算波场残差Sp=P -P.1 obs 1 cal 9 (3)、将波场残差反向传播得到回传波场P'； (4)、利用回传波场和正向传播波场计算目标函数的梯度，即速度的更新方向，并计算更新步长； (5)、对速度模型进行更新； (6)、检验是否满足迭代终止条件，满足则输出速度模型的反演结果，否则返回步骤(2)，继续迭代更新。
2.根据权利要求1所述的初至波波形反演近地表速度模型方法，其特征在于，所述声波方程交错网格有限差分正演模拟技术，利用一阶声波方程有限差分方法，数值模拟炮点激发的正向传播波场；非均勻介质一阶声波方程:应力-质点速度方程；时间2阶、空间12阶的高阶有限差分，PML边界条件，交错网格。
3.根据权利要求1所述的初至波波形反演近地表速度模型方法，其特征在于，所述最速下降法波形反演技术，⑴地面观测记录加时窗函数一初至波窗函数，作为一阶声波方程的初值条件，类似于炮点激发正向传播波场正演技术，数值模拟初至波反向传播的地震波场和数值模拟的正向传播波场求残差，利用L2模建立反演目标函数利用正向传播波场和波场残差函数，计算目标函数的最速下降方向即梯度，确定出速度扰动的最大变化方向；⑶利用梯度和Frechet导数计算出速度扰动的迭代步长#)利用速度更新方向和更新步长更新速度。
【文档编号】G01V1/28GK103499835SQ201310476248
【公开日】2014年1月8日 申请日期:2013年10月13日 优先权日:2013年10月13日
【发明者】胡自多, 韩令贺, 王宇超, 王述江, 雍运动, 寇龙江 申请人:中国石油集团西北地质研究所