储层地球物理敏感参数优选方法及系统技术方案

公开了一种储层地球物理敏感参数优选方法及系统。该方法可以包括：基于地质数据、测井曲线与地震资料，获取弹性参数曲线与参考曲线；依次对所有弹性参数对进行坐标旋转，计算新参数曲线与参考曲线的相关系数；获取最佳弹性参数对，储层敏感参数曲线，最佳旋转角度；获取最佳弹性参数对中的两个弹性参数数据体，对弹性参数数据体进行地震反演，通过坐标旋转公式与最佳旋转角度，计算敏感参数数据体；设置门槛值，对敏感参数数据体进行判别，预测优质储层。本发明专利技术可以消除储层地球物理敏感参数优选过程中人为因素引起的误差，准确、快速地优选出对储层最为敏感的参数，进而提升储层预测精度，为勘探开发提供更有力的技术支持。

Optimization method and system for reservoir geophysical sensitive parameters

A method and system for optimizing the parameters of reservoir geophysical sensitive parameters is disclosed. The method can include geological data, well logging and seismic data acquisition based on elastic parameter curve and the reference curve; turn on all elastic parameters of coordinate rotation, calculated the correlation coefficient of the new parameter curve and the reference curve; get the best elastic parameters on reservoir parameters, sensitivity curve, the optimal rotation angle of two elastic; the parameters in the data body to obtain the optimal elastic parameters of seismic inversion of elastic parameters data, through the coordinate rotation formula and the optimal rotation angle calculation, sensitive parameter data; set the threshold value, to determine the sensitive parameters data, the prediction of high-quality reservoir. The invention can eliminate the errors caused by human factors in the optimization process of geophysical sensitive parameters of reservoirs, accurately and quickly select the most sensitive parameters for reservoirs, and further improve the prediction accuracy of reservoirs, and provide more powerful technical support for exploration and development.

【技术实现步骤摘要】
储层地球物理敏感参数优选方法及系统
本专利技术涉及石油天然气勘探开发领域，更具体地，涉及一种储层地球物理敏感参数优选方法及系统。
技术介绍
地震相分析、地震属性分析以及地震反演等叠后储层预测技术是油藏描述中应用最为广泛的勘探技术，在实际应用中取得了较好的效果。然而随着石油天然气勘探和开发的逐步深入，有利目标隐蔽性越来越强，储层与非储层往往难以形成明显的速度、密度、阻抗或者弹性参数差异，因此使用单一的参数区判断储层与非储层往往存在较大的误差。如申请号为201310072499.2的专利基于多敏感参数的含气性储层预测方法及装置所提出的方法，解决的是相对比较简单的问题，即单一使用纵横波速度比或泊松比即可识别含气储层。但是实际问题往往更为复杂，尤其是针对较为致密的储层或者非常规储层，单一参数进行储层识别已经不能满足复杂岩性储层预测的要求。在利用多参数进行储层预测方面，通常包括两种方法：第一种为利用岩石物理定量解释模板来对岩性、物性、含气性进行定量解释，这在PerAvseth等人的专著《QuantitativeSeismicInterpretation》中有许多应用实例，然而这类实例往往解决的是高孔储层的定量预测问题，当储层较为致密时岩石物理建模技术却难以奏效；第二种方法如申请号为201310363844.8的专利提出的一种利用地震物相体预测储层的方法或如Quakenbush等人在文献《PoissonImpedance》提出的方法，其核心观点是利用坐标旋转的原理将现有坐标系统进行转换来构建新的能够进行岩性、物性、含气性识别的因子，进行坐标旋转的具体角度需人为确定，因此该方法存在人为因素造成的误差。因此，有必要开发一种储层地球物理敏感参数优选方法及系统。公开于本专利技术
技术介绍
部分的信息仅仅旨在加深对本专利技术的一般
技术介绍
的理解，而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。
技术实现思路
本专利技术提出了一种储层地球物理敏感参数优选方法及系统，其能够消除储层地球物理敏感参数优选过程中人为因素引起的误差，准确、快速地优选出对储层最为敏感的参数，进而提升储层预测精度，为勘探开发提供更有力的技术支持。根据本专利技术的一方面，提出了一种储层地球物理敏感参数优选方法。所述方法可以包括：基于地质数据、测井曲线与地震资料，获取弹性参数曲线与参考曲线；从弹性参数曲线中任意选择两条作为弹性参数对，依次对所有所述弹性参数对通过坐标旋转公式进行坐标旋转，分别计算每增加1度旋转后的新参数曲线与参考曲线的相关系数；获取相关系数最大对应的弹性参数对作为最佳弹性参数对，相关系数最大对应的新参数曲线作为储层敏感参数曲线，相关系数最大对应的旋转角度作为最佳旋转角度；获取所述最佳弹性参数对中的两个弹性参数数据体，对所述弹性参数数据体进行地震反演，通过所述坐标旋转公式与所述最佳旋转角度，计算敏感参数数据体；设置门槛值，对所述敏感参数数据体进行判别，预测优质储层。优选地，所述坐标旋转公式为：x＝x1×cosθ+x2×sinθ(1)其中，θ表示旋转角度，x1和x2分别为第一弹性参数曲线和第二弹性参数曲线，x为新参数曲线。优选地，所述相关系数为：其中，N为曲线采样点数目，xi为新参数曲线的采样点，yi为参考曲线的采样点。优选地，所述地震反演包括叠前地震反演与叠后地震反演。优选地，设置门槛值，对所述敏感参数数据体的孔隙度、含气饱和度、含油饱和度及非常规储层地化指标进行判别，预测优质储层。优选地，所述测井曲线包括纵波速度、横波速度、密度测井曲线以及矿物组分测井解释曲线、孔隙度测井解释曲线、含气饱和度曲线。优选地，还包括：基于所述纵波速度、横波速度、密度测井曲线计算阻抗曲线与弹性参数曲线。优选地，所述地震资料包括叠前地震道集与叠后成果数据。根据本专利技术的另一方面，提出了一种储层地球物理敏感参数优选系统，可以包括：存储器，存储有计算机可执行指令；处理器，所述处理器运行所述存储器中的计算机可执行指令，执行以下步骤：基于地质数据、测井曲线与地震资料，获取弹性参数曲线与参考曲线；从弹性参数曲线中任意选择两条作为弹性参数对，依次对所有所述弹性参数对通过坐标旋转公式进行坐标旋转，分别计算每增加1度旋转后的新参数曲线与参考曲线的相关系数；获取相关系数最大对应的弹性参数对作为最佳弹性参数对，相关系数最大对应的新参数曲线作为储层敏感参数曲线，相关系数最大对应的旋转角度作为最佳旋转角度；获取所述最佳弹性参数对中的两个弹性参数数据体，对所述弹性参数数据体进行地震反演，通过所述坐标旋转公式与所述最佳旋转角度，计算敏感参数数据体；设置门槛值，对所述敏感参数数据体进行判别，预测优质储层。优选地，所述坐标旋转公式为：x＝x1×cosθ+x2×sinθ(1)其中，θ表示旋转角度，x1和x2分别为第一弹性参数曲线和第二弹性参数曲线，x为新参数曲线。本专利技术的方法和装置具有其它的特性和优点，这些特性和优点从并入本文中的附图和随后的具体实施方式中将是显而易见的，或者将在并入本文中的附图和随后的具体实施方式中进行详细陈述，这些附图和具体实施方式共同用于解释本专利技术的特定原理。附图说明通过结合附图对本专利技术示例性实施例进行更详细的描述，本专利技术的上述以及其它目的、特征和优势将变得更加明显，其中，在本专利技术示例性实施例中，相同的参考标号通常代表相同部件。图1示出了根据本专利技术的储层地球物理敏感参数优选方法的步骤的流程图。图2示出了根据本专利技术的一个实施例的敏感参数优选程序的示意图。图3示出了根据本专利技术的一个实施例的弹性参数对参考系数随旋转角度的变化示意图。图4a和图4b分别示出了根据本专利技术的一个实施例的弹性参数交会图和新参数与弹性参数交会图的示意图。图5示出了根据本专利技术的一个实施例的新参数识别剖面的示意图。具体实施方式下面将参照附图更详细地描述本专利技术。虽然附图中显示了本专利技术的优选实施例，然而应该理解，可以以各种形式实现本专利技术而不应被这里阐述的实施例所限制。相反，提供这些实施例是为了使本专利技术更加透彻和完整，并且能够将本专利技术的范围完整地传达给本领域的技术人员。图1示出了根据本专利技术的储层地球物理敏感参数优选方法的步骤的流程图。在该实施例中，根据本专利技术的储层地球物理敏感参数优选方法可以包括：步骤101，基于地质数据、测井曲线与地震资料，获取弹性参数曲线与参考曲线；在一个示例中，测井曲线包括纵波速度、横波速度、密度测井曲线以及矿物组分测井解释曲线、孔隙度测井解释曲线、含气饱和度曲线。在一个示例中，还包括：基于纵波速度、横波速度、密度测井曲线计算阻抗曲线与弹性参数曲线。在一个示例中，地震资料包括叠前地震道集与叠后成果数据。步骤102，从弹性参数曲线中任意选择两条作为弹性参数对，依次对所有弹性参数对通过坐标旋转公式进行坐标旋转，分别计算每增加1度旋转后的新参数曲线与参考曲线的相关系数。在一个示例中，坐标旋转公式为：x＝x1×cosθ+x2×sinθ(1)其中，θ表示旋转角度，x1和x2分别为第一弹性参数曲线和第二弹性参数曲线，x为新参数曲线，该曲线用于与参考曲线计算相关性。步骤103，获取相关系数最大对应的弹性参数对作为最佳弹性参数对，相关系数最大对应的新参数曲线作为储层敏感参数曲线，相关系数最大对应的旋转本文档来自技高网...

【技术保护点】
一种储层地球物理敏感参数优选方法，包括：基于地质数据、测井曲线与地震资料，获取弹性参数曲线与参考曲线；从弹性参数曲线中任意选择两条作为弹性参数对，依次对所有所述弹性参数对通过坐标旋转公式进行坐标旋转，分别计算每增加1度旋转后的新参数曲线与参考曲线的相关系数；获取相关系数最大对应的弹性参数对作为最佳弹性参数对，相关系数最大对应的新参数曲线作为储层敏感参数曲线，相关系数最大对应的旋转角度作为最佳旋转角度；获取所述最佳弹性参数对中的两个弹性参数数据体，对所述弹性参数数据体进行地震反演，通过所述坐标旋转公式与所述最佳旋转角度，计算敏感参数数据体；设置门槛值，对所述敏感参数数据体进行判别，预测优质储层。

【技术特征摘要】
1.一种储层地球物理敏感参数优选方法，包括：基于地质数据、测井曲线与地震资料，获取弹性参数曲线与参考曲线；从弹性参数曲线中任意选择两条作为弹性参数对，依次对所有所述弹性参数对通过坐标旋转公式进行坐标旋转，分别计算每增加1度旋转后的新参数曲线与参考曲线的相关系数；获取相关系数最大对应的弹性参数对作为最佳弹性参数对，相关系数最大对应的新参数曲线作为储层敏感参数曲线，相关系数最大对应的旋转角度作为最佳旋转角度；获取所述最佳弹性参数对中的两个弹性参数数据体，对所述弹性参数数据体进行地震反演，通过所述坐标旋转公式与所述最佳旋转角度，计算敏感参数数据体；设置门槛值，对所述敏感参数数据体进行判别，预测优质储层。2.根据权利要求1所述的储层地球物理敏感参数优选方法，其中，所述坐标旋转公式为：x＝x1×cosθ+x2×sinθ(1)其中，θ表示旋转角度，x1和x2分别为第一弹性参数曲线和第二弹性参数曲线，x为新参数曲线。3.根据权利要求1所述的储层地球物理敏感参数优选方法，其中，所述相关系数为：其中，N为曲线采样点数目，xi为新参数曲线的采样点，yi为参考曲线的采样点。4.根据权利要求1所述的储层地球物理敏感参数优选方法，其中，所述地震反演包括叠前地震反演与叠后地震反演。5.根据权利要求1所述的储层地球物理敏感参数优选方法，其中，设...

【专利技术属性】
技术研发人员：郭旭升，李文成，王良军，苏克露，蒲勇，朱岩松，劳牡丹，陈超，李雷涛，
申请(专利权)人：中国石油化工股份有限公司，中国石油化工股份有限公司勘探分公司，
类型：发明
国别省市：北京,11