储层孔隙度定量预测方法及装置制造方法及图纸

本发明专利技术公开了一种储层孔隙度定量预测方法及装置，该方法包括：根据目标储层段的岩石基质物理参数基于HM岩石物理理论模型计算纵波阻抗参数以及纵横波速度比参数；利用二重对数域的纵波阻抗参数和二重对数域的纵横波速度比参数构建近线性岩石物理解释模板；基于所述近线性岩石物理解释模板拟合孔隙度趋势线，确定孔隙度与二重对数域的纵波阻抗参数、二重对数域的纵横波速度比参数的关系；根据所述纵波阻抗地震数据体、所述纵横波速度比地震数据体以及所述孔隙度与二重对数域的纵波阻抗参数、二重对数域的纵横波速度比参数的关系确定目标储层段的孔隙度。本发明专利技术降低了储层孔隙度预测的多解性，提高了储层孔隙度预测的精度。

【技术实现步骤摘要】
储层孔隙度定量预测方法及装置
本专利技术涉及石油地球物理勘探领域，具体而言，涉及一种储层孔隙度定量预测方法及装置。
技术介绍
随着石油天然气开发程度的提高，石油天然气勘探的难度不断增大，勘探目标的复杂性越来越强，使得对储层预测的精度提出越来越高的要求，促使储层描述逐渐由定性检测发展到定量预测。常规的储层预测技术通常基于属性优选、AVO分析以及地震反演等方法实现，这些方法在一定程度上解决了储层定性检测的问题，对于相对简单的、单一的构造油藏具有较好的应用效果，比如利用属性优选结合井约束的方法，可以分析优选出与井参数吻合度最高的地震属性，根据该属性的平面分布特征可以实现对储层物性参数的平面展布预测，圈定有利区，提高储层预测的精度。同样的，基于AVO分析以及地震反演方法可以得到具有特定意义且与某种储层参数吻合的属性参数，并利用属性参数表征储层参数分布。在勘探初期，上述方法能够满足勘探精度的要求，然而当进入精细勘探阶段，剩余圈闭类型复杂化、油气水系统多样化、孔隙流体非均匀化等问题开始涌现，常规的定性检测技术已经无法达到准确评估油气藏潜力的目的，需要新的定量化的检测手段对复杂储层进行全面的精细的刻画。储层预测的核心问题是识别地下介质中的孔隙以及孔隙中的流体，其中孔隙在很大程度上决定了储层的规模以及潜力，准确预测储层孔隙分布规律是储层预测的重点。随着物探技术的发展，尤其是叠前反演技术的成熟，基于叠前反演数据体直接预测孔隙度的方法受到越来越多的重视，现有的方法通常根据测井数据拟合的孔隙度与地震物理参数之间的经验关系计算储层孔隙度分布，在拟合经验关系式的过程中人为因素影响严重，导致预测的结果因人而异并且误差难以控制，大大增加了储层预测的多解性，降低了预测精度。
技术实现思路
本专利技术为了解决上述技术问题中的至少一个，提出了一种储层孔隙度定量预测方法及装置。为了实现上述目的，根据本专利技术的一个方面，提供了一种储层孔隙度定量预测方法，该方法包括：根据目标储层段的岩石基质物理参数基于HM岩石物理理论模型计算纵波阻抗参数以及纵横波速度比参数，其中，孔隙度为所述纵波阻抗参数和所述纵横波速度比参数中的变量；利用二重对数域的纵波阻抗参数和二重对数域的纵横波速度比参数构建近线性岩石物理解释模板；基于所述近线性岩石物理解释模板拟合孔隙度趋势线，确定孔隙度与二重对数域的纵波阻抗参数、二重对数域的纵横波速度比参数的关系；获取反演得到的目标储层段的纵波阻抗地震数据体和纵横波速度比地震数据体；根据所述纵波阻抗地震数据体、所述纵横波速度比地震数据体以及所述孔隙度与二重对数域的纵波阻抗参数、二重对数域的纵横波速度比参数的关系确定目标储层段的孔隙度。可选的，所述岩石基质物理参数包括：岩石基质体积模量、岩石基质剪切模量、岩石基质密度、参考孔隙度、泊松比、孔隙压力以及孔隙流体的密度。可选的，所述利用二重对数域的纵波阻抗参数和二重对数域的纵横波速度比参数构建近线性岩石物理解释模板，具体包括：分别对所述纵波阻抗参数和所述纵横波速度比参数取二重对数，得到二重对数域的纵波阻抗参数和二重对数域的纵横波速度比参数；将所述二重对数域的纵波阻抗参数作为横坐标，将所述二重对数域的纵横波速度比参数作为纵坐标，构建近线性岩石物理解释模板。可选的，所述基于所述近线性岩石物理解释模板拟合孔隙度趋势线，包括：根据所述线性岩石物理解释模板中孔隙度趋势线的分布特征，拟合不同孔隙度对应的孔隙度趋势线的线性表达式，其中，不同孔隙度对应的线性表达式为：为第i个孔隙度，Xi为第i个孔隙度对应的二重对数域的纵波阻抗参数，Yi为第i个孔隙度对应的二重对数域的纵横波速度比参数，k1、k2、C为系数。可选的，所述确定孔隙度与二重对数域的纵波阻抗参数、二重对数域的纵横波速度比参数的关系，具体包括：基于最小二乘法拟合不同孔隙度对应的孔隙度趋势线的线性表达式，使拟合的结果满足关系式：其中，min(SST)表示误差平方和达到最小，Yi为不同孔隙度对应的二重对数域的纵横波速度比参数，为不同孔隙度对应的二重对数域的纵横波速度比参数的平均值；根据拟合的结果确定系数k1、k2、C，以确定孔隙度对二重对数域的纵波阻抗参数、二重对数域的纵横波速度比参数的表达式，其中表达式为：其中，为孔隙度、Y为二重对数域的纵横波速度比参数、X为二重对数域的纵波阻抗参数。为了实现上述目的，根据本专利技术的另一方面，提供了一种储层孔隙度定量预测装置，该装置包括：纵波阻抗及纵横波速度比计算单元，用于根据目标储层段的岩石基质物理参数基于HM岩石物理理论模型计算纵波阻抗参数以及纵横波速度比参数，其中，孔隙度为所述纵波阻抗参数和所述纵横波速度比参数中的变量；近线性岩石物理解释模板构建单元，用于利用二重对数域的纵波阻抗参数和二重对数域的纵横波速度比参数构建近线性岩石物理解释模板；孔隙度关系确定单元，用于基于所述近线性岩石物理解释模板拟合孔隙度趋势线，确定孔隙度与二重对数域的纵波阻抗参数、二重对数域的纵横波速度比参数的关系；纵波阻抗及纵横波速度比地震数据体获取单元，用于获取反演得到的目标储层段的纵波阻抗地震数据体和纵横波速度比地震数据体；储层孔隙度计算单元，用于根据所述纵波阻抗地震数据体、所述纵横波速度比地震数据体以及所述孔隙度与二重对数域的纵波阻抗参数、二重对数域的纵横波速度比参数的关系确定目标储层段的孔隙度。可选的，所述岩石基质物理参数包括：岩石基质体积模量、岩石基质剪切模量、岩石基质密度、参考孔隙度、泊松比、孔隙压力以及孔隙流体的密度。可选的，所述近线性岩石物理解释模板构建单元包括：二重对数计算模块，用于分别对所述纵波阻抗参数和所述纵横波速度比参数取二重对数，得到二重对数域的纵波阻抗参数和二重对数域的纵横波速度比参数；坐标系建立模块，用于将所述二重对数域的纵波阻抗参数作为横坐标，将所述二重对数域的纵横波速度比参数作为纵坐标，构建近线性岩石物理解释模板。可选的，所述孔隙度关系确定单元包括：线性表达式建立模块，用于根据所述线性岩石物理解释模板中孔隙度趋势线的分布特征，拟合不同孔隙度对应的孔隙度趋势线的线性表达式，其中，不同孔隙度对应的线性表达式为：为第i个孔隙度，Xi为第i个孔隙度对应的二重对数域的纵波阻抗参数，Yi为第i个孔隙度对应的二重对数域的纵横波速度比参数，k1、k2、C为系数。可选的，所述孔隙度关系确定单元还包括：拟合模块，用于基于最小二乘法拟合不同孔隙度对应的孔隙度趋势线的线性表达式，使拟合的结果满足关系式：其中，min(SST)表示误差平方和达到最小，Yi为不同孔隙度对应的二重对数域的纵横波速度比参数，为不同孔隙度对应的二重对数域的纵横波速度比参数的平均值；孔隙度表达式确定模块，用于根据拟合的结果确定系数k1、k2、C，以确定孔隙度对二重本文档来自技高网...

【技术保护点】
1.一种储层孔隙度定量预测方法，其特征在于，包括：/n根据目标储层段的岩石基质物理参数基于HM岩石物理理论模型计算纵波阻抗参数以及纵横波速度比参数，其中，孔隙度为所述纵波阻抗参数和所述纵横波速度比参数中的变量；/n利用二重对数域的纵波阻抗参数和二重对数域的纵横波速度比参数构建近线性岩石物理解释模板；/n基于所述近线性岩石物理解释模板拟合孔隙度趋势线，确定孔隙度与二重对数域的纵波阻抗参数、二重对数域的纵横波速度比参数的关系；/n获取反演得到的目标储层段的纵波阻抗地震数据体和纵横波速度比地震数据体；/n根据所述纵波阻抗地震数据体、所述纵横波速度比地震数据体以及所述孔隙度与二重对数域的纵波阻抗参数、二重对数域的纵横波速度比参数的关系确定目标储层段的孔隙度。/n

【技术特征摘要】
1.一种储层孔隙度定量预测方法，其特征在于，包括：
根据目标储层段的岩石基质物理参数基于HM岩石物理理论模型计算纵波阻抗参数以及纵横波速度比参数，其中，孔隙度为所述纵波阻抗参数和所述纵横波速度比参数中的变量；
利用二重对数域的纵波阻抗参数和二重对数域的纵横波速度比参数构建近线性岩石物理解释模板；
基于所述近线性岩石物理解释模板拟合孔隙度趋势线，确定孔隙度与二重对数域的纵波阻抗参数、二重对数域的纵横波速度比参数的关系；
获取反演得到的目标储层段的纵波阻抗地震数据体和纵横波速度比地震数据体；
根据所述纵波阻抗地震数据体、所述纵横波速度比地震数据体以及所述孔隙度与二重对数域的纵波阻抗参数、二重对数域的纵横波速度比参数的关系确定目标储层段的孔隙度。


2.根据权利要求1所述的储层孔隙度定量预测方法，其特征在于，所述岩石基质物理参数包括：岩石基质体积模量、岩石基质剪切模量、岩石基质密度、参考孔隙度、泊松比、孔隙压力以及孔隙流体的密度。


3.根据权利要求1所述的储层孔隙度定量预测方法，其特征在于，所述利用二重对数域的纵波阻抗参数和二重对数域的纵横波速度比参数构建近线性岩石物理解释模板，具体包括：
分别对所述纵波阻抗参数和所述纵横波速度比参数取二重对数，得到二重对数域的纵波阻抗参数和二重对数域的纵横波速度比参数；
将所述二重对数域的纵波阻抗参数作为横坐标，将所述二重对数域的纵横波速度比参数作为纵坐标，构建近线性岩石物理解释模板。


4.根据权利要求1所述的储层孔隙度定量预测方法，其特征在于，所述基于所述近线性岩石物理解释模板拟合孔隙度趋势线，包括：
根据所述线性岩石物理解释模板中孔隙度趋势线的分布特征，拟合不同孔隙度对应的孔隙度趋势线的线性表达式，其中，不同孔隙度对应的线性表达式为：




为第i个孔隙度，Xi为第i个孔隙度对应的二重对数域的纵波阻抗参数，Yi为第i个孔隙度对应的二重对数域的纵横波速度比参数，k1、k2、C为系数。


5.根据权利要求4所述的储层孔隙度定量预测方法，其特征在于，所述确定孔隙度与二重对数域的纵波阻抗参数、二重对数域的纵横波速度比参数的关系，具体包括：
基于最小二乘法拟合不同孔隙度对应的孔隙度趋势线的线性表达式，使拟合的结果满足关系式：



其中，min(SST)表示误差平方和达到最小，Yi为不同孔隙度对应的二重对数域的纵横波速度比参数，为不同孔隙度对应的二重对数域的纵横波速度比参数的平均值；
根据拟合的结果确定系数k1、k2、C，以确定孔隙度对二重对数域的纵波阻抗参数、二重对数域的纵横波速度比参数的表达式，其中表达式为：



其中，为孔隙度、Y为二重对数域的纵横波速度比参数、X为二重对数域的纵波阻抗参数。


6.一种储层孔隙度定量预测装置，其特征在于，包括：
纵波阻抗及纵横波速度比计算单元，用于根据目标储层段的岩石基质物理参数基于HM岩石物理理论模型计算纵波阻抗参数以及纵横波速度比参数，其中，孔隙度为所述纵波阻抗参数和所述纵横波速度比参...

【专利技术属性】
技术研发人员：王磊，石兰亭，方乐华，陈彬滔，徐中华，史忠生，薛罗，马轮，史江龙，
申请(专利权)人：中国石油天然气股份有限公司，
类型：发明
国别省市：北京;11