基于PG属性的火山岩烃类检测方法、存储介质及电子设备

技术领域

本发明涉及火山岩烃类检测技术领域，具体地涉及一种基于PG属性的火山岩烃类检测方法、一种存储介质和一种电子设备。

背景技术

AVO技术是烃类检测的关键技术手段。早在20世纪80年代，Ostrander利用AVO技术识别流体。随后，Smilth等提出了流体因子的概念。Goodway利用λρ－μρ(λ、μ为拉梅系数，ρ为密度)交会分析检测烃类。Russell利用纵波波阻抗与横波波阻抗构建了新的流体指示因子。Quakenbush等探讨了基于泊松阻抗检测流体的方法。宁忠华等通过引入纵波、横波波阻抗的高阶项提高流体因子的敏感性。王栋利用纵波、横波波阻抗多次幂的组合形式压制噪声。张玉洁等基于喷流效应对AVO计算P(截距)和G(梯度)、P*G、流体因子进行了推广和应用。姜仁等利用Russell流体因子检测致密储层的含气性。郑静静等对各类流体因子进行了分类与系统总结。谢玉洪等认为“暗点”型油气藏的流体指示因子可以放大Ⅱ类AVO异常的差异。李坤等利用混合概率模型驱动的叠前地震反演方法刻画复杂岩性油气藏、直接检测流体。

但在实际应用过程中存在如下问题：

1、PG属性往往适应于层状介质的砂泥岩储层的烃类检测，不适用于与横向非均质性强的火山岩储层；

2、P*G流体因子往往受诸多因素影响存在“流体假象”。

发明内容

本发明实施方式的目的是提供一种基于PG属性的火山岩烃类检测方法，该方法在AVO技术的基础上，设计含气因子数据体来对火山岩的烃类进行检测，以解决P*G属性和流体因子等不适用于火山岩烃类检测的问题。

为了实现上述目的，本发明第一方面提供一种基于PG属性的火山岩烃类检测方法，包括：

分析火山岩的AVO特征，基于火山岩的截距属性P、梯度属性G、以及临界入射角，设计含气因子计算公式，得到含气因子数据体；

根据计算结果分析识别火山岩烃类分布情况。

本方法针对火山岩岩性复杂、储层横向变化快等特点，由PG属性衍生设计含气因子对火山岩烃类的分布进行检测，解决了常规PG属性不适用于火山岩烃类检测的问题。

可选的，还包括：

采集火山岩CRP道集，基于所述CRP道集分析火山岩的临界入射角，计算CRP道集的截距属性P和梯度属性G。

基于采集的CRP道集地震数据，分析计算临界入射角、截距属性P和梯度属性G，保证数据的准确性。

可选的，基于所述CRP道集分析火山岩的临界入射角，具体包括：

分析CRP道集，获取目的层深度和最大偏移距，基于目的层深度和最大偏移距分析计算临界入射角。

临界入射角的计算公式如下：

θ＝arctan(offset/2/H)；

其中，θ表示临界入射角，H表示目的层深度，offset表示最大偏移距。

可选的，计算CRP道集的截距属性P的公式如下：

其中，所述P为截距，Vp为纵波速度，ρ为密度。

可选的，计算CRP道集的梯度属性G的公式如下：

μ＝ρVs

其中，G为梯度，Vp为纵波速度；Vs为横波速度；ρ为密度；μ为剪切模量。

可选的，所述设计含气因子计算公式，包括：

通过临界入射角对截距属性P和梯度属性G进行坐标旋转，建立含气因子计算公式。

所述含气因子的计算公式如下：

Gas＝P*cos(θ/180)+G*sin(θ/180)；

其中，Gas表示含气因子，G为梯度，P为截距，θ为临界入射角。

通过临界入射角对截距属性和梯度属性进行坐标旋转来放大截距值和减小梯度值，旋转后的P、G相加放大二者的差异，实现识别烃类的目的，计算的结果值越大，含气性越强，

本发明第二方面提供一种计算机可读存储介质，所述计算机可读存储介质上存储有计算机指令，所述计算机指令在计算机上执行时，使得所述计算机执行第一方面所述的方法。

本发明第三方面提供一种电子设备，包括：至少一个处理器、存储器；

所述存储器存储计算机执行指令；

所述至少一个处理器执行所述存储器存储的计算机执行指令，使得所述电子设备执行第一方面所述的方法。

本发明针对火山岩岩性复杂、储层横向变化快等特点，提出一种基于P、G和临界入射角的含气因子识别算法对火山岩储层进行烃类检测，该项技术具有地震属性确定、实现快速、识别准确的特点。

本发明实施方式的其它特征和优点将在随后的具体实施方式部分予以详细说明。

附图说明

附图是用来提供对本发明实施方式的进一步理解，并且构成说明书的一部分，与下面的具体实施方式一起用于解释本发明实施方式，但并不构成对本发明实施方式的限制。在附图中：

图1是本发明一种实施方式提供的基于PG属性的火山岩烃类检测方法流程图。

具体实施方式

以下结合附图对本发明的具体实施方式进行详细说明。应当理解的是，此处所描述的具体实施方式仅用于说明和解释本发明，并不用于限制本发明。

图1是本发明一种实施方式提供的基于PG属性的火山岩烃类检测方法流程图。如图1所示，本发明实施方式提供一种基于PG属性的火山岩烃类检测方法，所述方法包括：

收集火山岩的CRP道集信息，分析火山岩的临界入射角，并计算截距属性和梯度属性；针对火山岩的AVO特征，设计一个由P(截距)、G(梯度)属性衍生的算法，该算法通过临界入射角对PG属性进行坐标旋转实现，该算法的计算结果值越大，表明火山岩含气性越强，将该算法的计算结果数据体命名为含气因子，含气因子本质是一种基于叠前地震数据AVO的特征响应计算的P、G属性组合参数；将火山岩的临界入射角、截距属性和梯度属性带入算法中进行计算，基于计算结果识别火山岩烃类分布。

数据体是一类多维矩阵，让用户从多个角度探索和分析数据集，通常是一次同时考虑三个因素(维度)。数据体是二维表格的多维扩展，如同几何学中立方体是正方形的三维扩展一样，便于基于计算结果进行分析。本实施方式记载的含气因子数据体是一个时间域三维数据体(平面上是XY表示距离，单位是m，纵向上为时间，单位ms)。

具体的，在本实施例中，火山岩的临界入射角的计算公式如下：

θ＝arctan(offset/2/H) (1)

其中，θ表示最佳临界入射角，H表示目的层深度，offset表示最大偏移距。

例如，本实施例中，收集松南断陷长岭地区火山岩CRP道集，收集的数据中，目的层深度H为定值4000米，最大偏移距offset为一个范围值：3860-5240米，在对最大偏移距设置的时候需要考虑在各个地震道上覆盖次数分布要均匀，当最大偏移距为4620米时，各个地震道上的覆盖次数为92次，且均匀分布，通过公式(1)计算得到最佳临界入射角θ为30。

具体的，在本实施例中，截距属性和梯度属性的计算公式如下：

μ＝ρVs

其中，P为截距；G为梯度；θ为临界入射角；Vp为纵波速度；Vs为横波速度；ρ为密度；μ为剪切模量。

上述公式(2)、(3)、(4)基于AVO理论得出，AVO的理论基础是Zoepritz方程，由于该方程较复杂，涉及参数多，使用起来比较困难，因此，许多学者对其进行了简化，目前使用较多的近似公式是Aki-Richards公式和Shuey公式,其中Shuey公式在AVO属性提取中使用较多，其表达式如下：

R(θ)≈P+Gsin

其中，R为地震的反射系数，由公式(5)即可推导出公式(2)、公式(3)、公式(4)。

火山岩储层横向变化快，纵向多期叠置，属于各向异性介质，对于火山岩的AVO特征为：当火山岩储集层含水时，P(截距)为负值，G(梯度)为正值，当火山岩储集层含气时，P为正值，G为负值。

针对火山岩的这个特点，具体的，在本实施例中，通过临界入射角θ对P、G进行坐标旋转，达到放大P值和减小G值的目的，坐标旋转后的P与G相加，得到含气因子的计算公式，如下：

Gas＝P*cos(θ/180)+G*sin(θ/180)(6)

其中，Gas表示含气因子，G为梯度，P为截距，θ为临界入射角。计算得到的结果值越大，表示火山岩的含气性就越强。基于得到的含气因子数据体，即可分析得出火山岩烃类的分布情况，达到烃类检测的目的。

本实施例还提供一种计算机可读存储介质，所述计算机可读存储介质上存储有计算机指令，所述计算机指令在计算机上执行时，使得所述计算机执行本实施例提供的一种基于PG属性的火山岩烃类检测方法。

本实施例还提供一种电子设备，包括：至少一个处理器、存储器；

所述存储器存储计算机执行指令；

所述至少一个处理器执行所述存储器存储的计算机执行指令，使得所述电子设备执行本实施例提供的一种基于PG属性的火山岩烃类检测方法。

本领域技术人员可以理解实现上述实施方式的方法中的全部或部分步骤是可以通过程序来指令相关的硬件来完成，该程序存储在一个存储介质中，包括若干指令用以使得单片机、芯片或处理器(processor)执行本发明各个实施方式所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM，Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。

以上结合附图详细描述了本发明的可选实施方式，但是，本发明实施方式并不限于上述实施方式中的具体细节，在本发明实施方式的技术构思范围内，可以对本发明实施方式的技术方案进行多种简单变型，这些简单变型均属于本发明实施方式的保护范围。另外需要说明的是，在上述具体实施方式中所描述的各个具体技术特征，在不矛盾的情况下，可以通过任何合适的方式进行组合。为了避免不必要的重复，本发明实施方式对各种可能的组合方式不再另行说明。

此外，本发明的各种不同的实施方式之间也可以进行任意组合，只要其不违背本发明实施方式的思想，其同样应当视为本发明实施方式所公开的内容。