专利名称：基于波动方程外推算子的vti介质中叠前时间偏移方法
技术领域：
本发明涉及一种叠前时间偏移技术，尤其涉及一种一种VTI介质中叠前时间偏移方法。
背景技术：
在地质构造不甚复杂的区域，叠前时间偏移后的成像道集上的AVO(振幅随偏移距的变化)/AVA(振幅随入射角的变化)分析是油气储层预测的关键技术。叠前时间偏移技术在提高地震成像精度，为叠前反演的实现等方面发挥了重要作用。目前的VTI (具有垂直对称轴的横向各向异性)介质时间偏移基本上是使用Kirchhoff积分法VTI叠前时间偏移算法，采用的是用等效各项异性参数表示的四阶方程(Allchalifah，Τ, 1997, Velocity analysis using nonhyperbolic moveout in transversely isotropic media,. Geophysics，62，1839-1854)，其本质上还是基于层状介质假设下推导的时距关系。通常所采用的商用软件系统中基本上都是用Kirchhoff积分法叠前时间偏移给出Offset成像道集，并在此基础上进行AVO分析。实际上，KirchhofT偏移方法使用的是均方根速度场，适应速度横向变化的能力较弱，当速度场横向变化时，基于常速介质中的格林函数计算得到的振幅加权因子通常是不精确的。因此，Kirchhoff积分法给出的成像道集振幅不保真、成像道集中存在假象干扰，这样的成像道集不利于后续的AVO分析和储层预测。吸收衰减和各向异性的对成像道集和AVO分析的影响也要在成像过程中加以消除。还有，随着各油田探区油气勘探程度的不断加深，对偏移成像精度的要求越来越高，各向同性处理技术已经难以满足勘探开发的要求。迄今为止，还未见有基于波动方程外推算子的叠前时间偏移算法。与KirchhofT 偏移方法相比而言，波动方程外推算子对振幅的处理更加自然，并且对速度横向变化更强的适应能力。因此，有必要发展一套基于波动方程外推算子的VTI介质中叠前时间偏移方法。
发明内容
针对上述现有技术，从适应缓横向变速、提高成像分辨率(吸收衰减补偿)和考虑各项异性等方面综合考虑，本发明提供一种基于波动方程外推算子的VTI介质中叠前时间偏移方法。是将偏移距平面波波动方程叠前时间偏移方法至VTI介质中，从而得到VTI介质中基于波动方程的叠前时间偏移算法，本发明具有更好的保振幅能力，可以提供保真的成像道集，为后续的AV0/AVA分析提供更好的基础数据，有利于更精确的油气储层岩性预测。 本发明是一种更精确高效的波动方程叠前时间偏移方法，能补偿大地的吸收衰减和处理各向异性介质中的成像，可以提供很好的叠前时间偏移叠加剖面，输出射线参数或者角度域成像道集用于后续的剩余速度分析和AVO分析。为了解决上述技术问题，本发明一种基于波动方程外推算子的VTI介质中叠前时间偏移方法予以实现的技术方案包括以下步骤
步骤1 输入CMP道集以及时间域速度场；步骤2 对CMP道集进行局部倾斜叠加，将CMP道集映射到截距时间-射线参数域的平面波数据体；并按照射线参数对平面波数据体分选，形成共射线参数剖面；步骤3 对射线参数循环；步骤4 进程判断所在计算机节点是否为主控节点，由此进入主-从并行模式；即 若进程在主控节点上运行，则执行下述(5-A)过程，若进程未在主控节点上运行，而是在计算节点上运行，则执行下述(5-B)过程；(5-A)首先，读取当前射线参数对应的共射线参数剖面，然后沿时间方向作一维快速傅立叶变换，生成频率域的共射线参数平面波数据体；与此同时，设定总任务数为总频率数目与计算节点数目之中的最小值；主控节点控制着并行任务的分发与收集，并执行以下步骤(5-A1)对总任务数循环；(5-A2)将一个频率切片的平面波数据体发给计算节点，发送任务的计数器自增；(5-A3)判断任务发送是否结束，若否，转向(5-A1)；若是，转向(5_A4)；(5-A4)对总频率数循环；(5-A5)接收计算节点返回的任务标识符；(5-A6)检查已经发送给计算节点的任务序号是否小于总频率数目；若是，转向 (5-A7)；若否，转向(5-A8)；(5-A7)将当前频率切片的平面波数据体发给计算节点，发送任务的计数器自增， 转向(5-A4)；(5-A8)向所有计算节点发送任务结束标识符；(5-B)首先，判断该计算节点号是否大于总的频率数目；若是，转向步骤6 ；若否， 转向(5-B1)；(5-B1)进程进入死循环；(5-B2)进程等待接收从主控节点发送来的数据；(5-B3)进程判断主控节点是否发送结束标识符；若是，转向(5-B9)；若否，转向 (5-B4)；(5-B4)沿时间方向进行波场外推；(5-B5)计算偏移距平面波传播的相速度求解如下隐式方程
权利要求
1. 一种基于波动方程外推算子的VTI介质中叠前时间偏移方法，其特征在于包括以下步骤步骤1 输入CMP道集以及时间域速度场；步骤2 对CMP道集进行局部倾斜叠加，将CMP道集映射到截距时间-射线参数域的平面波数据体；并按照射线参数对平面波数据体分选，形成共射线参数剖面； 步骤3 对射线参数循环；步骤4 进程判断所在计算机节点是否为主控节点，由此进入主-从并行模式； 即若进程在主控节点上运行，则执行下述(5-A)过程，若进程未在主控节点上运行， 而是在计算节点上运行，则执行下述(5-B)过程；(5-A)首先，读取当前射线参数对应的共射线参数剖面，然后沿时间方向作一维快速傅立叶变换，生成频率域的共射线参数平面波数据体；与此同时，设定总任务数为总频率数目与计算节点数目之中的最小值；主控节点控制着并行任务的分发与收集，并执行以下步骤(5-A1)对总任务数循环；(5-A2)将一个频率切片的平面波数据体发给计算节点，发送任务的计数器自增； (5-A3)判断任务发送是否结束，若否，转向(5-A1)；若是，转向(5-A4)； (5-A4)对总频率数循环； (5-A5)接收计算节点返回的任务标识符；(5-A6)检查已经发送给计算节点的任务序号是否小于总频率数目；若是，转向 (5-A7)；若否，转向(5-A8)；(5-A7)将当前频率切片的平面波数据体发给计算节点，发送任务的计数器自增，转向 (5-A4)；(5-A8)向所有计算节点发送任务结束标识符；(5-B)首先，判断该计算节点号是否大于总的频率数目；若是，转向步骤6;若否，转向 (5-B1)；(5-B1)进程进入死循环；(5-B2)进程等待接收从主控节点发送来的数据；(5-B3)进程判断主控节点是否发送结束标识符；若是，转向(5-B9)；若否，转向 (5-B4)；(5-B4)沿时间方向进行波场外推；(5-B5)计算偏移距平面波传播的相速度求解如下隐式方程m^i sin^ -Ph^(η公式⑴中，平面波相速度νρ( θ ;νρ0, ε , δ)可以表示为VqP (f) =- + gsin26> + -J(l + 2gsin2θ)1 — 2(ε — sin22Θ (2) νρ0 2 2公式O)中，ε、δ为VTI介质各向异性参数，θ为相角，Vptl SVTI介质中沿对称轴 (Ζ轴)方向地震波的相速度；通过公式(1) 0)，可以得出由射线参数Ph表征的偏移距平面波传播的相速度vp (ph)；(5-B6)对VTI介质中的偏移距平面波方程采用隐式有限差分格式进行波场外推，实现偏移距平面波叠前时间偏移；所采用的偏移方程为
全文摘要
本发明公开了一种基于波动方程外推算子的VTI介质中叠前时间偏移方法，包括(1)输入CMP道集以及时间域速度场；(2)对CMP道集进行局部倾斜叠加，将CMP道集映射到截距时间-射线参数域的平面波数据体；并按照射线参数对平面波数据体分选，形成共射线参数剖面；(3)对射线参数循环；(4)进程判断所在计算机节点是否为主控节点，由此进入主-从并行模式；(6)所有进程同步；(7)主控节点将当前射线参数对应的成像数据体输出；(8)判断射线参数循环是否完成；若否，转向步骤3；(9)输出时间偏移剖面以及成像道集。本发明可以提供保真的成像道集，为后续的AVO/AVA分析提供更好的基础数据，有利于更精确的油气储层岩性预测。
文档编号G01V1/36GK102313903SQ20111018273
公开日2012年1月11日 申请日期2011年7月1日 优先权日2011年7月1日
发明者冯波, 史浩, 吕丁友, 吴奎, 孙书滨, 明君, 王华忠, 田立新, 耿建华, 贺电波 申请人:中国海洋石油总公司, 中海石油(中国)有限公司天津分公司, 同济大学