[发明专利]高精度道积分计算相对波阻抗的方法

说明书

本发明涉及地震勘探的储层预测中的地震资料储层反演技术，尤其涉及高精度道积分计算相对波阻抗的方法，包括构建Robinson多级算子，计算地层反射系数、计算地层反射系数积分。本发明提供的高精度道积分计算相对波阻抗的方法，只需要常规地震数据，可以快速的得到高精度的相对波阻抗数据，从而实现快速的评价反映地下介质的岩性信息的波阻抗的相对变化情况。实现了高精度相对波阻抗计算方法，有力的促进了相对波阻抗反演技术在实际地震勘探中的应用，对于地震储层解释中的岩性参数在空间的分布规律有着重要的参考价值，对指导油气的勘探开发，提高油气勘探成功率有极大的帮助。

技术领域

本发明涉及地震勘探的储层预测中的地震资料储层反演技术，尤其涉及高精度道积分计算相对波阻抗的方法。

背景技术

随着地震勘探的逐渐深入,地震勘探逐步从构造油气藏向复杂岩性油气藏以及隐蔽油气藏转变,地震储层预测的难度越来越大，对地震资料储层反演的要求也越来越高。道积分反演是在地下地层波阻抗是深度变量连续且可微函数的假设条件下，通过地震道自相关方法求取地层反射系数，再对反射系数进行深度积分计算出地层波阻抗。由于道积分反演实际上只是对反射系数的积分，反演波阻抗数据体是相对值，道积分算子的精度严重影响了相对波阻抗反演的准确性和可靠性。

在油气地震勘探中，当地下地质体的岩性、物性发生变化时，经常会引起岩石密度和地震波传播速度的变化，导致波阻抗产生变化，波阻抗在岩性、物性界面处的变化会生成一个反射系数脉冲，当地震波传播到波阻抗界面时，就会在其界面处产生反射波。因此，若已知地下介质的波阻抗，通过求取反射系数与并地震子波褶积运算得到地震波通过地下介质后观测记录的过程就是地震波阻抗正演。当我们根据观测到的地震记录，通过地震波反射的褶积理论等求取反映地下介质的岩性信息的波阻抗等信息的过程就是地震波阻抗反演。地震资料是反射界面的表现，是一种界面型信息，只能间接表达地下地层的地质特征，不能直接进行储层的描述。而作为储层预测重要工具的地震反演技术可以将地震数据转换成波阻抗数据，波阻抗数据由于是地层型剖面，可以和测井地质信息直接对比，可以直观的进行储层的识别和描述。

目前常用的相对波阻抗计算方法基于均匀弹性地球模型，地震响应可以看成是震源子波与地层脉冲响应的褶积，包含在地震数据中的地震子波与地层的反射率函数是相互独立的。道积分反演根据地震波传播理论，用数学计算工具消除地震记录中的子波、弱化白化噪音，反射界面的反射系数序列，再通过积分方法求取地下波阻抗。通过积分方法把地震道记录直接转换为波阻抗数据，反演过程具有积分误差小、计算简单、对计算机要求较低的优点，但是道积分反演过程仅仅依赖于地震资料，地震资料的品质和带宽对反演结果影响很大。因为地震资料由于采集、处理等原因，缺失低频和高频信息，频带较窄，所以该方法求取的波阻抗纵向分辨率低、精度小。道积分反演方法无法使用测井等资料约束，因此反演结果是相对值、纵向分辨率不高，适用于地震勘探初期以及没有井数据的研究区。本发明提出了应用多级子波分解技术的高精度积分来计算相对波阻抗的道积分方法，将地震剖面“转换”成分辨率更高的相对波阻抗剖面，得到的相对波阻抗剖面不仅便于解释人员将地震资料与测井资料连接对比，而且能有效地对地层岩性参数的变化进行研究，从而得到岩性参数在空间的分布规律，指导油气的勘探开发。

发明的内容

本发明的目的是克服现有道积分计算相对波阻抗基于简单均匀弹性介质的缺点，提出了一种可以适用于复杂介质的多级子波高精度道积分计算相对波阻抗的技术，用来提高相对波阻抗反演的准确性和可靠性。

具体技术方案为：

高精度道积分计算相对波阻抗的方法，包括以下过程：

构建Robinson地震褶积三级模型：

x(t)＝f₃(f₂(f₁(r(t)*b₁(t))*b₂(t))*b₃(t)) (1)