一种基于两参数的多波联合叠前属性提取方法及系统与流程

本发明涉及叠前ava属性提取领域，更具体地，涉及一种基于两参数的多波联合叠前属性提取方法及系统。

背景技术：

叠前地震反演是进行储层岩性识别和流体预测的主要技术之一，在研究储层性质时，从地震波资料中提取的储层岩性或含油气有关的属性可直接进行储层岩性与流体预测。由于复杂油气藏的分布日趋分散复杂，勘探难度越来越大多。多波多分量地震勘探技术可以为我们提供地下弹性波的矢量波场，全面地提供地下储层的信息，可以作为研究复杂油气藏的有效手段之一，多波地震资料由于拥有更多的地震波信息被广泛用于储层岩性识别和提取含油气有关的联合属性中，在地震资料解释中应用广泛。地震属性技术广泛应用于地震构造解释、地层分析、油藏特征描述以及油藏动态监测等各个领域，为油田开发后期的准确油气预测提供强有力的支持。

在叠前地震勘探技术中，人们研发了属性提取方法检测油气。流体指示因子对油气具有一定的敏感性而被广泛应用，杨氏模量能够凸显页岩气储层的页岩特征，表征岩石的脆性。因此，有必要开发一种基于两参数的多波联合叠前属性提取方法及系统。

公开于本发明背景技术部分的信息仅仅旨在加深对本发明的一般背景技术的理解，而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。

技术实现要素：

本发明提出了一种基于两参数的多波联合叠前属性提取方法及系统，其能够通过对敏感指示因子选取进而建立反射系数近似方程，利用方程结果最终完成储层识别和油气预测，对avo叠前反演储层识别及后续处理有重要意义。

根据本发明的一方面，提出了一种基于两参数的多波联合叠前属性提取方法。所述方法可以包括：

1)选取敏感的流体指示因子ρf和岩性指示因子eρ；

2)基于步骤1)中获得的所述流体指示因子和所述岩性指示因子获得纵波反射系数近似方程和转换横波反射系数近似方程；

3)基于zoeppritz方程计算得到的纵波和转换横波反射系数与雷克子波褶积进行正演得到纵波和转换横波角道集，与步骤2)中获得的纵波反射系数近似方程和转换横波反射系数近似方程进行联合反演，获得包含有流体指示因子和岩性指示因子的反演结果；

4)基于步骤3)中的反演结果，根据流体指示因子变化率的异常状态和岩性指示因子变化率的投影结果进行储层识别和油气预测。

优选地，在步骤2)中，基于所述流体指示因子和所述岩性指示因子获得的纵波反射系数方程、转换横波反射系数方程用如下：

其中，ρ＝(ρ1+ρ2)/2是平均密度，δρ＝ρ2-ρ1，ρ1、ρ2分别为界面上、下层的密度；θ＝(θ1+θ2)/2是p波入射角和透射角的平均，θ1、θ2分别为p波入射角和透射角；是psv波反射角和透射角的平均，分别为psv转换波反射角和透射角，γ表示纵横波速度比，γdy表示干岩石下纵横波速度比。

优选地，在步骤3)中，所述雷克子波的频率为35hz。

优选地，在步骤4)中，流体因子变化率呈负异常状态，判定所述异常状态处含有油气。

优选地，在步骤4)中，岩性因子变化率投影至属性剖面根据属性剖面的结果判定岩性。

根据本发明的另一方面，提出了一种基于两参数的多波联合叠前属性提取系统，其上存储有计算机程序，其中，所述程序被处理器执行时实现以下步骤：

步骤1：选取敏感的流体指示因子ρf和岩性指示因子eρ；

步骤2：基于步骤1中获得的所述流体指示因子和所述岩性指示因子获得纵波反射系数近似方程和转换横波反射系数近似方程；

步骤3：基于zoeppritz方程计算得到的纵波和转换横波反射系数与雷克子波褶积进行正演得到纵波和转换横波角道集，与步骤2中获得的纵波反射系数近似方程和转换横波反射系数近似方程进行联合反演，获得包含有流体指示因子和岩性指示因子的反演结果；

步骤4：基于步骤3中的反演结果，根据流体指示因子变化率的异常状态和岩性指示因子变化率的投影结果进行储层识别和油气预测。

优选地，在步骤2中，基于所述流体指示因子和所述岩性指示因子获得的纵波反射系数方程、转换横波反射系数方程用如下：

其中，ρ＝(ρ1+ρ2)/2是平均密度，δρ＝ρ2-ρ1，ρ1、ρ2分别为界面上、下层的密度；θ＝(θ1+θ2)/2是p波入射角和透射角的平均，θ1、θ2分别为p波入射角和透射角；是psv波反射角和透射角的平均，分别为psv转换波反射角和透射角，γ表示纵横波速度比，γdy表示干岩石下纵横波速度比。

优选地，在步骤3中，所述雷克子波的频率为35hz。

优选地，在步骤4中，流体因子变化率呈负异常状态，判定所述异常状态处含有油气。

优选地，在步骤4中，岩性因子变化率投影至属性剖面根据属性剖面的结果判定岩性。

本发明具有其它的特性和优点，这些特性和优点从并入本文中的附图和随后的具体实施方式中将是显而易见的，或者将在并入本文中的附图和随后的具体实施方式中进行详细陈述，这些附图和具体实施方式共同用于解释本发明的特定原理。

附图说明

通过结合附图对本发明示例性实施例进行更详细的描述，本发明的上述以及其它目的、特征和优势将变得更加明显，其中，在本发明示例性实施例中，相同的参考标号通常代表相同部件。

图1示出了根据本发明的一种基于两参数的多波联合叠前属性提取方法的步骤的流程图；

图2示出了流体指示因子和岩性指示因子的交会图；

图3(a)示出了ostrander模型界面二的纵波反射系数曲线图；

图3(b)示出了ostrander模型界面二的转换横波反射系数曲线图；

图4(a)示出了纵波反演流体指示因子反射率反演结果、纵波和转换横波反演流体指示因子反演结果流体指示因子反射率反演结果与反演结果真实值的对比图；

图4(b)示出了纵波反演岩性指示因子反射率反演结果、纵波和转换横波反演岩性指示因子反演结果流体指示因子反射率反演结果与反演结果真实值的对比图；

图5(a)示出了1000cdp点处的纵波角道集；

图5(b)示出了1000cdp点处的转换横波角道集；

图6示出了目标区放大显示的纵波与转换横波资料联合属性提取结果；

图7示出了目标区放大显示的纵波与转换横波资料联合属性提取结果。

具体实施方式

下面将参照附图更详细地描述本发明。虽然附图中显示了本发明的优选实施例，然而应该理解，可以以各种形式实现本发明而不应被这里阐述的实施例所限制。相反，提供这些实施例是为了使本发明更加透彻和完整，并且能够将本发明的范围完整地传达给本领域的技术人员。

图1示出了根据本发明的一种基于两参数的多波联合叠前属性提取方法的步骤的流程图。

在该实施例中，根据本发明的一种基于两参数的多波联合叠前属性提取方法可以包括：

步骤101：选取敏感的流体指示因子ρf和岩性指示因子eρ；

其中，e为杨氏模量，ρ为密度，f为russell等2003年提出的流体指示因子。

在一个示例中，castagna和smith(1994)收集到的25组含水砂岩、页岩、含气砂岩的速度与密度数据代表了世界范围内的采样，包含了常见的三类含气砂岩异常。利用上述数据进行流体指示因子和岩性指示因子的交会图制作，如图2所示，可以很好地在图中将含气砂岩、含水砂岩与页岩进行区分，进而证明了所选取的敏感流体指示因子和岩性指示因子的正确性，能够达到将含气砂岩、含水砂岩和页岩进行区分。

步骤102：基于步骤101中获得的所述流体指示因子和所述岩性指示因子获得纵、横波反射系数近似方程；

在一个示例中，通过以下方程进行表示纵波反射系数及转换波反射系数：

其中，vp＝(vp1+vp2)/2是平均p波速度，δvp＝vp2-vp1，其中vp1、vp2分别为界面上、下层的纵波速度；vs＝(vs1+vs2)/2是平均sv波速度，δvs＝vs2-vs1，vs1、vs2分别为界面上、下层的sv波速度；ρ＝(ρ1+ρ2)/2是平均密度，δρ＝ρ2-ρ1，ρ1、ρ2分别为界面上、下层的密度；θ＝(θ1+θ2)/2是p波入射角和透射角的平均，θ1、θ2分别为p波入射角和透射角；是psv波反射角和透射角的平均，分别为psv转换波反射角和透射角。

在一个示例中，基于流体指示因子和岩性指示因子获得的纵波反射系数方程、转换横波反射系数方程用如下：

其中，rpp是纵波反射系数，rps是转换横波反射系数，ρ＝(ρ1+ρ2)/2是平均密度，δρ＝ρ2-ρ1，ρ1、ρ2分别为界面上、下层的密度；θ＝(θ1+θ2)/2是p波入射角和透射角的平均，θ1、θ2分别为p波入射角和透射角；是psv波反射角和透射角的平均，分别为psv转换波反射角和透射角，γ表示纵横波速度比，γdy表示干岩石下纵横波速度比。

在示例性实施例中，步骤102中获得的纵波反射系数近似方程和转换横波反射系数近似方程还需要与已有近似方程进行精度对比，本发明利用ostrander的三层模型计算纵、横波反射系数，进行精度的对比，进行精度对比的方程为：zoeppritz方程、aki-richards近似式、russell近似式、smith和gidlow近似式。如图3(a)所示，为对ostrander三层模型通过正演计算得到的纵波反射系数曲线，如图3(b)所示，为对ostrander三层模型通过正演计算得到的转换横波反射系数曲线，从两张图的整体而言，曲线与标准曲线误差在0.1左右，趋势临近，利用经验近似式进行简化对新推导的纵横波反射系数近似式的精度影响作用较小。

更具体的，还需要利用mahmoudianandmargrave的四层层状模型测试纵波反射系数近似方程和转换横波反射系数近似方程。

步骤103：基于zoeppritz方程计算得到的纵波和转换横波反射系数与雷克子波褶积进行正演得到纵波和转换横波角道集，与步骤2中获得的纵波反射系数近似方程和转换横波反射系数近似方程进行联合反演，获得包含有流体指示因子和岩性指示因子的反演结果；

在示例性实施例中，将经过zoeppritz精确方程计算的反射系数，与35hz雷克子波与得到的纵波和转换横波反射系数进行褶积，获得纵波和转换波角道集，或者也可以由实际资料预处理得到纵波和转换波角道集。

在示例性实施例中，可以将每个采样点处的模型数据，代入zoeppritz方程进而得到不超过临界角在深度域内的随角度变化的纵波反射系数和转换横波反射系数。

如图4(a)所示，为纵波反演流体指示因子反射率反演结果、纵波和转换横波反演流体指示因子反演结果流体指示因子反射率反演结果与真实反演结果的对比图，可以看出纵波和转换横波联合avo反演的横波速度和密度反射率精度高于仅用纵波反演的结果，如图(b)所示，纵波反演岩性指示因子反射率反演结果、纵波和转换横波反演岩性指示因子反演结果流体指示因子反射率反演结果与真实反演结果的对比图，可以看出纵波和转换横波联合avo反演的横波速度和密度反射率精度高于仅用纵波反演的结果。

步骤104：基于步骤103中的反演结果，根据流体指示因子变化率的异常状态和岩性指示因子变化率的投影结果进行储层识别和油气预测。

在一个示例中，流体因子变化率呈负异常状态，判定所述异常状态处含有油气；岩性因子变化率投影至属性剖面根据属性剖面的结果判定岩性。

本发明还提出了一种基于两参数的多波联合叠前属性提取系统，其上存储有计算机程序，其中，所述程序被处理器执行时实现以下步骤：

步骤1：选取敏感的流体指示因子ρf和岩性指示因子eρ；

步骤2：基于步骤1中获得的所述流体指示因子和所述岩性指示因子获得纵、横波反射系数近似方程；

步骤3：基于zoeppritz方程计算得到的纵波和转换横波反射系数与雷克子波褶积进行正演得到纵波和转换横波角道集，与步骤2中获得的纵波反射系数近似方程和转换横波反射系数近似方程进行联合反演，获得包含有流体指示因子和岩性指示因子的反演结果；

步骤4：基于步骤3中的反演结果，根据流体指示因子变化率的异常状态和岩性指示因子变化率的投影结果进行储层识别和油气预测。

应用示例

为便于理解本发明实施例的方案及其效果，以下给出一个具体应用示例。本领域技术人员应理解，该示例仅为了便于理解本发明，其任何具体细节并非意在以任何方式限制本发明。

为了检验多波联合ava属性在实际资料的应用效果，采用中国lj地区a油田m测线的纵波和转换横波数据进行测试。该区在地质上位于凹陷南缓斜坡，受燕山运动和喜马拉雅运动的影响，形成了向北倾没的大型鼻状构造，油气藏类型主要为受构造控制的断块油气藏，兼有岩性、泥岩裂缝等特殊油气藏。采用的地震资料pp波主频接近35hz，ps波主频接近15hz。在时间域内将纵波地震数据与转换横波地震数据匹配一致得到纵波角道集如图5(a)所示，转换横波角道集如图5(b)所示。

经过基于两参数的多波联合叠前属性提取的流体指示因子变化率结果剖面如图6所示，对有井目的层进行局部放大并将过测线的生产井投影到属性剖面上，含有油气时，纵波速度减小，而横波速度对流体相对不敏感，造成流体因子相对变化率负异常，体现在属性分析结果剖面中的蓝色异常值区域，证明了流体指示因子变化率属性结果剖面对油气的识别效果较好。

经过基于两参数的多波联合叠前属性提取的岩性指示因子变化率结果剖面如图7所示，有井目的层进行局部放大显示，生物灰岩位于深灰区域，泥页岩位于剖面的浅色区域，属性对泥页岩、生物灰岩具有较好的识别效果。

综上所述，本发明通过反应目标层的储层特征，较好的指示了油气并识别了岩性，理论严谨可靠且操作流程简单实用，为ava属性提取与油气预测提供了帮助。

本领域技术人员应理解，上面对本发明的实施例的描述的目的仅为了示例性地说明本发明的实施例的有益效果，并不意在将本发明的实施例限制于所给出的任何示例。

以上已经描述了本发明的各实施例，上述说明是示例性的，并非穷尽性的，并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。