[发明专利]储层孔隙度定量预测方法及装置

权利要求书

1.一种储层孔隙度定量预测方法，其特征在于，包括：

根据目标储层段的岩石基质物理参数基于HM岩石物理理论模型计算纵波阻抗参数以及纵横波速度比参数，其中，孔隙度为所述纵波阻抗参数和所述纵横波速度比参数中的变量；

利用二重对数域的纵波阻抗参数和二重对数域的纵横波速度比参数构建近线性岩石物理解释模板；

基于所述近线性岩石物理解释模板拟合孔隙度趋势线，确定孔隙度与二重对数域的纵波阻抗参数、二重对数域的纵横波速度比参数的关系；

获取反演得到的目标储层段的纵波阻抗地震数据体和纵横波速度比地震数据体；

根据所述纵波阻抗地震数据体、所述纵横波速度比地震数据体以及所述孔隙度与二重对数域的纵波阻抗参数、二重对数域的纵横波速度比参数的关系确定目标储层段的孔隙度。

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

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

分别对所述纵波阻抗参数和所述纵横波速度比参数取二重对数，得到二重对数域的纵波阻抗参数和二重对数域的纵横波速度比参数；

将所述二重对数域的纵波阻抗参数作为横坐标，将所述二重对数域的纵横波速度比参数作为纵坐标，构建近线性岩石物理解释模板。

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

根据所述线性岩石物理解释模板中孔隙度趋势线的分布特征，拟合不同孔隙度对应的孔隙度趋势线的线性表达式，其中，不同孔隙度对应的线性表达式为：

为第i个孔隙度，X_i为第i个孔隙度对应的二重对数域的纵波阻抗参数，Y_i为第i个孔隙度对应的二重对数域的纵横波速度比参数，k₁、k₂、C为系数。

5.根据权利要求4所述的储层孔隙度定量预测方法，其特征在于，所述确定孔隙度与二重对数域的纵波阻抗参数、二重对数域的纵横波速度比参数的关系，具体包括：

基于最小二乘法拟合不同孔隙度对应的孔隙度趋势线的线性表达式，使拟合的结果满足关系式：

其中，min(SST)表示误差平方和达到最小，Y_i为不同孔隙度对应的二重对数域的纵横波速度比参数，为不同孔隙度对应的二重对数域的纵横波速度比参数的平均值；

根据拟合的结果确定系数k₁、k₂、C，以确定孔隙度对二重对数域的纵波阻抗参数、二重对数域的纵横波速度比参数的表达式，其中表达式为：

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

6.一种储层孔隙度定量预测装置，其特征在于，包括：

纵波阻抗及纵横波速度比计算单元，用于根据目标储层段的岩石基质物理参数基于HM岩石物理理论模型计算纵波阻抗参数以及纵横波速度比参数，其中，孔隙度为所述纵波阻抗参数和所述纵横波速度比参数中的变量；

近线性岩石物理解释模板构建单元，用于利用二重对数域的纵波阻抗参数和二重对数域的纵横波速度比参数构建近线性岩石物理解释模板；

孔隙度关系确定单元，用于基于所述近线性岩石物理解释模板拟合孔隙度趋势线，确定孔隙度与二重对数域的纵波阻抗参数、二重对数域的纵横波速度比参数的关系；

纵波阻抗及纵横波速度比地震数据体获取单元，用于获取反演得到的目标储层段的纵波阻抗地震数据体和纵横波速度比地震数据体；

储层孔隙度计算单元，用于根据所述纵波阻抗地震数据体、所述纵横波速度比地震数据体以及所述孔隙度与二重对数域的纵波阻抗参数、二重对数域的纵横波速度比参数的关系确定目标储层段的孔隙度。