一种叠前地震储层岩性参数的反演方法技术

本发明专利技术涉及地震油气勘探、开发领域，提供了一种叠前地震储层岩性参数的反演方法，其包括地震梯度阻抗储层岩性参数反演，可以预测储层发育。

【技术实现步骤摘要】

本专利技术涉及地震油气勘探、开发领域，特别涉及一种叠前地震岩性参数的反演方法。
技术介绍
AVO(AmplitudeVersusOffset，振幅随偏移距的变化)技术用于研究地震反射振幅随炮点与检波器之间的距离即炮检距(或入射角)的变化特征来探讨反射系数响应随炮检距(或入射角)的变化，进而确定反射界面上覆、下伏介质的岩性特征及物性参数。油气勘探开发技术人员借助AVO技术可以更好地评估油气藏储层的岩石属性。与叠后地震含油气信息相比，由于叠前地震振幅包含了丰富的地层岩性与油气信息，经过许多研究者技术攻关，近几十年AVO技术得到了迅速的发展。在地震波Zoeppritz方程的基础上，发展了许多应用简便的AVO的近似表达式，如Aki和Rechards(1980)、Shuey(1985)、Bortfeid(1987)、Mallick(1993)近似表达式等，极大地推动了AVO技术的应用。但是由于影响岩石物性(如速度、泊淞比)变化的因素较多，再加上地震子波分辨率的影响，造成AVO储层预测技术存在多解性。为了提高储层的描述精度，虽然发展了井震约束的各种叠前储层参数反演技术，但应用条件苛刻，计算成本很高，应用受到限制。简便、有效的应用AVO叠前地震属性进行岩性与含油气检测技术有待进一步研究发展。虽然叠前AVO属性参数很多，它们对储层岩性和含油气性的敏感性各不相同，但在AVO属性分析中最常用的是Shuey的表达式，即AVO截距P和梯度G属性。截距P表示反射界面两侧的P波阻抗变化的响应，而斜率G表示振幅随入射角(或炮检距)的变化率，它是岩层弹性参数的综合反映，包含了反射界面两侧的横波速度、纵波速度和密度变化的信息，是岩性识别的敏感参数。然而，上述AVO技术存在地震识别岩性的精度低，油气勘探开发的风险高的问题。因此，有必要针对该技术进行改进，以提高地震识别岩性的精度，降低油气勘探开发的风险。
技术实现思路
目前在实际生产科研中被广泛应用地震地层纵波反射系数R(θ)是Shuey的AVO二项近似式：R(θ)＝P+Gsin2θ，其中，θ是地震纵波的入射角，P是地震AVO的截距属性，G是地震AVO的梯度属性。本专利技术针对AVO截距属性P对应储层岩性参数为纵波阻抗AI＝VPρ(其中，AI表示纵波的声波阻抗，Vp是储层的纵波速度，ρ为储层的密度)，而提出了对应梯度G属性的储层岩性参数梯度阻抗的反演技术，储层的地震梯度阻抗是AVO技术和叠前弹性参数反演技术的继承和发展，可以提高地震识别岩性的精度，降低油气勘探开发的风险。因此，本专利技术之一提供了一种叠前地震储层岩性参数的反演方法，其包括地震梯度阻抗储层岩性参数反演。在一个具体实施例中，所述地震梯度阻抗储层岩性参数反演为梯度阻抗的构造函数，且所述梯度阻抗函数中包括储层的S波(横波)速度、P波(纵波)速度和密度的信息。在一个具体实施例中，所述构造函数为其中，GI为梯度阻抗，Vp是储层的P波速度，Vs是储层的S波速度，ρ为储层的密度，并且在一个具体实施例中，依据储层的声波阻抗的值域对所述梯度阻抗GI进行标准化。其目的是使储层的梯度阻抗与P波阻抗具有相同的值域。在一个具体实施例中，对井储层的梯度阻抗的岩性识别进行统计分析，然后对所述井储层的梯度阻抗进行频谱分析，根据井储层的梯度阻抗的频谱和井附近储层的叠前地震梯度属性数据的频谱加权，求取算子O的频谱，使算子O的频谱具有井梯度阻抗的低频信息、高频信息及地震梯度的主要频谱信息；而后再将求得的算子O应用于叠前地震属性的梯度数据体，通过公式Z＝O*S可得到叠前地震属性的梯度阻抗；其中，Z表示储层岩性参数梯度阻抗，O表示褶积算子，*为褶积算子符号，S表示叠前地震属性梯度数据体的地震道。其中，低频信息是指总能量低于地震波总能量10％的低频成分的地震能量。地震低频也是一个相对概念，不同的地震数据体的低频成分不同。高频信息是指总能量低于地震波总能量10％的高频成分的地震能量。地震高频也是一个相对概念，不同的地震数据体的高频成分不同。地震梯度的主要频谱信息是指以地震主频为中心频带内的地震能量大于80％地震总能量的地震频率成分。其中的地震主频是指在地震频谱分析中频谱分析曲线最大值所对应的峰值频率。本专利技术之二提供了一种如上所述的方法在预测储层发育中的应用。本专利技术之三提供了一种使用如上所述的方法预测储层发育的方法，其包括如下步骤：(1)测定井储层的Vp、Vs和ρ，然后根据测定的所述的Vp、Vs和ρ构建井储层梯度阻抗GI曲线，其中，Vp是储层的P波速度，Vs是储层的S波速度，ρ为储层的密度，GI为梯度阻抗；(2)针对井储层进行梯度阻抗GI的岩性识别敏感性分析；(3)井储层的梯度阻抗数据在对数域的振幅频率交汇分析；(4)线性拟合梯度阻抗谱与频率的关系式L＝Fα，其中L表示储层的梯度阻抗，F表示频率，α为梯度阻抗关于变量频率F变化的衰减系数；(5)联合应用井储层梯度阻抗数据的频谱L和井附近储层地震梯度的频谱加权平均，作为储层梯度阻抗的平均频谱，计算反演算子O的频谱；(6)将(5)求得的算子O的频谱进行-90°相移，最终形成一个时间域的反演算子O，然后进行地震梯度体褶积运算得到梯度阻抗反演；(7)在井震标定和层位解释的基础上，求取储层层间的梯度阻抗均方根值的平面分布图，根据储层的岩性特征预测储层发育有利区。其中，步骤(5)中的L与上述O的不同之处在于：算子O的频谱具有L的低频和高频信息，同时保留了地震梯度的有效频宽和地震梯度包含的主要频谱信息。在预测储层发育的一个具体实施例中，且并且在预测储层发育的一个具体实施例中，步骤(6)中的井储层梯度阻抗数据的频谱L和井附近储层地震梯度的频谱均应用傅里叶变换得到。本专利技术提出的一种叠前地震储层岩性参数的反演方法，是井储层岩性信息约束的叠前地震梯度属性的拓广算法，通过井储层岩性信息与地震梯度属性有效融合，构建出一种具有储层岩性指示敏感的参数，提高了地震识别储层岩性的精度。该技术涉及参数少，计算成本低，应用简便，能够较好地识别砂体，适用于地震品质较好、AVO技术有效的油气田的勘探开发。该技术方法试验性地在西非海上某油田开发中，发现因其优化了油田开发方案，因而能够取得很好的应用效果，取得良好的经济效益，例如增加了可采原油550万桶，其产生的经济效益可达35亿多人民币，因而具有极大的应用前景。附图说明图1为某油田S井储层的梯度阻抗曲线、测井伽马曲线、P波阻抗曲线(灰色为砂岩)。其中，图1-A为梯度阻抗曲线(GI曲线)；图1-B为测井伽马曲线；图1-C为P波阻抗曲线(AI曲线)。图2为某油田S井的P波阻抗、梯度阻抗岩性统计直方图(粗黑线为砂岩，细黑线为泥岩)。图2-A为P波阻抗的岩性统计直方图；图2-B为梯度阻抗的岩性统计直方图。图3为某油田S井储层的梯度阻抗的对数域频谱分析。图4为求得的算子时间响应特征(左图)及算子的频谱特征(右图)。图5为过井(测井伽马曲线)的地震反射剖面(白色表示低值，黑色表示高值)。图6为求得与图5对应的过井(测井伽马曲线)梯度阻抗岩性剖面(白色表示砂岩，黑色表示泥岩)。图7为根据地震振幅求得目标的砂体平面分布图(白色表示砂体发育区)。图8为根据梯度阻抗体求得目标的砂体平面分布图(白色表示砂体发育区)。具体本文档来自技高网...

【技术保护点】
一种叠前地震储层岩性参数的反演方法，其包括地震梯度阻抗储层岩性参数反演。

【技术特征摘要】
1.一种叠前地震储层岩性参数的反演方法，其包括地震梯度阻抗储层岩性参数反演。2.根据权利要求1所述的方法，其特征在于，所述地震梯度阻抗储层岩性参数反演为梯度阻抗的构造函数，且所述梯度阻抗函数中包括储层的S波速度、P波速度和密度的信息。3.根据权利要求2所述的方法，其特征在于，所述构造函数为其中，GI为梯度阻抗，Vp是储层的P波速度，Vs是储层的S波速度，ρ为储层的密度，并且4.根据权利要求1-3任意一项所述的方法，其特征在于，依据储层的声波阻抗的值域对所述梯度阻抗GI进行标准化。5.根据权利要求1-4任意一项所述的方法，其特征在于，对井储层的梯度阻抗的岩性识别进行统计分析，然后对所述井储层的梯度阻抗进行频谱分析，根据井储层的梯度阻抗的频谱和井附近储层的叠前地震梯度属性数据的频谱加权，求取算子O的频谱，使算子O的频谱具有井梯度阻抗的低频信息、高频信息及地震梯度的主要频谱信息；而后再将求得的算子O应用于叠前地震属性的梯度数据体，通过公式Z＝O*S可得到叠前地震属性的梯度阻抗；其中，Z表示储层岩性参数梯度阻抗，O表示褶积算子，*为褶积算子符号，S表示叠前地震属性梯度数据体的地震道。6.一种如权利要求1-5任意一项所述的方法在...

【专利技术属性】
技术研发人员：苑书金，高君，陈立莹，
申请(专利权)人：中国石油化工股份有限公司，中国石油化工股份有限公司石油勘探开发研究院，
类型：发明
国别省市：北京;11