一种计算砂砾岩体储层孔隙度方法技术

本发明专利技术公开了一种计算砂砾岩体储层孔隙度方法，其特征在于，包括以下步骤：S1、提出按照岩性类别建立砂砾岩体储层孔隙度模型，开展岩性类别划分；S2、提取对孔隙度敏感的测井响应曲线；S3、按砂砾岩体储层岩石类别建立岩心实测孔隙度与测井响应曲线数据表；S4、采用多元线性回归方法建立孔隙度计算模型，将其分为三个计算模型。与体积模型计算孔隙度比较，本发明专利技术技术方案优势在于充分利用多元测井信息优选对孔隙度敏感测井响应曲线，强化“岩性刻度测井技术”理念，克服体积模型对骨架参数的依赖性，计算精度明显提高，为复杂砂砾岩体储层提供了较为精准的孔隙度评价方法。层提供了较为精准的孔隙度评价方法。层提供了较为精准的孔隙度评价方法。

【技术实现步骤摘要】
一种计算砂砾岩体储层孔隙度方法


[0001]本专利技术涉及计算孔隙度
，具体为一种计算砂砾岩体储层孔隙度方法。

技术介绍

[0002]利津北带砂砾岩体储层岩性复杂，岩性包括砾岩、细砾岩、含砾砂岩、粗砂岩、中砂岩、细砂岩、粉砂岩、泥质粉砂岩、泥岩以及碳酸盐岩等，主要矿物成分为长石和石英，颗粒之间填隙物由杂基和胶结物组成，胶结物以白云石为主。目前砂砾岩体储层孔隙度主要采用体积模型来计算，计算的关键就是要准确确定岩石骨架参数值，但是由于砂砾岩体储层矿物组分复杂，导致岩石骨架参数值不易准确确定，从而影响体积模型计算孔隙度的精度。
[0003]因此需要一种有效准确确定孔隙度的方法，提高砂砾岩体储层孔隙度的计算精度。

技术实现思路

[0004]本专利技术的目的在于提供一种计算砂砾岩体储层孔隙度方法，以解决上述
技术介绍
中提出的问题。
[0005]为了解决上述技术问题，本专利技术提供如下技术方案：一种计算砂砾岩体储层孔隙度方法，其特征在于，包括以下步骤：
[0006]S1、提出按照岩性类别建立砂砾岩体储层孔隙度模型，开展岩性类别划分；
[0007]S2、提取对孔隙度敏感的测井响应曲线；
[0008]S3、按砂砾岩体储层岩石类别建立岩心实测孔隙度与测井响应曲线数据表；
[0009]S4、采用多元线性回归方法建立孔隙度计算模型，将其分为三个计算模型：
[0010](1)砂砾岩体砾岩类储层孔隙度计算模型，砾岩类储层孔隙度计算模型为POR＝a*AC+b*CNL+c*DEN+d；
[0011](2)砂砾岩体砂泥岩类储层孔隙度计算模型，砂泥岩类储层孔隙度计算模型为POR＝a*AC+b*CNL+c*DEN+d；
[0012](3)砂砾岩体碳酸盐岩类储层孔隙度计算模型，碳酸盐岩类储层孔隙度计算模型为POR＝a*AC+b*CNL+c*DEN+d。
[0013]与现有技术相比，本专利技术所达到的有益效果是：
[0014]与体积模型计算孔隙度比较，本专利技术技术方案优势在于充分利用多元测井信息优选对孔隙度敏感测井响应曲线，强化“岩性刻度测井技术”理念，克服体积模型对骨架参数的依赖性，计算精度明显提高，为复杂砂砾岩体储层提供了较为精准的孔隙度评价方法。
附图说明
[0015]附图用来提供对本专利技术的进一步理解，并且构成说明书的一部分，与本专利技术的实施例一起用于解释本专利技术，并不构成对本专利技术的限制。在附图中：
[0016]图1是本专利技术整体的流程示意图。
[0017]图1.1是孔隙度与声波时差交会检验重叠图。
[0018]图1.2是孔隙度与中子交会检验重叠图。
[0019]图1.3是孔隙度与密度交会检验重叠图。
[0020]图2.1是孔隙度与声波时差交会检验重叠图。
[0021]图2.2是孔隙度与中子交会检验重叠图。
[0022]图2.3是孔隙度与密度交会检验重叠图。
[0023]图3.1是孔隙度与声波时差交会检验重叠图。
[0024]图3.2是孔隙度与中子交会检验重叠图。
[0025]图3.3是孔隙度与密度交会检验重叠图。
具体实施方式
[0026]下面将结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本专利技术一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本专利技术保护的范围。
[0027]一种计算砂砾岩体储层孔隙度方法，其特征在于，包括以下步骤：
[0028]S1、提出按照岩性类别建立砂砾岩体储层孔隙度模型，开展岩性类别划分，按照岩性类别将砂砾岩体储层划分为砾岩类储层、砂泥岩类储层和碳酸盐岩类储层。砾岩类储层主要岩性包括砾岩、砂砾岩、细砾岩、砾状砂岩、含砾砂岩；砂泥岩类储层主要岩性包括粗砂岩、中砂岩、细砂岩、粉砂岩、泥质(粉)砂岩；碳酸盐岩类储层主要岩性包括白云岩、泥质白云岩、灰质白云岩、白云质(粉、细)砂岩、灰质砂岩等；
[0029]S2、提取对孔隙度敏感的测井响应曲线，对孔隙度(POR)敏感的测井响应曲线主要包括声波时差(AC)、密度(DEN)和中子(CNL)。其中孔隙度计量单位采用“％”，声波时差(AC)计量单位采用“μs/ft”，中子(CNL) 计量单位采用“％”，密度(DEN)计量单位采用“g/cm3”；
[0030]S3、按砂砾岩体储层岩石类别建立岩心实测孔隙度与测井响应曲线数据表利用岩心物性数据，结合测井资料，通过岩心归位，建立实测孔隙度与测井曲线数据表，表头见下表：
[0031][0032]S4、采用多元线性回归方法建立孔隙度计算模型，将其分为三个计算模型：
[0033](1)砂砾岩体砾岩类储层孔隙度计算模型，砾岩类储层孔隙度计算模型为POR＝a*AC+b*CNL+c*DEN+d；
[0034](2)砂砾岩体砂泥岩类储层孔隙度计算模型，砂泥岩类储层孔隙度计算模型为POR＝a*AC+b*CNL+c*DEN+d；
[0035](3)砂砾岩体碳酸盐岩类储层孔隙度计算模型，碳酸盐岩类储层孔隙度计算模型为POR＝a*AC+b*CNL+c*DEN+d。
[0036]实施例一：
[0037]砂砾岩体砾岩类储层孔隙度计算模型：
[0038]利用5口井90个砂砾岩体砾岩类储层样本，将样本实测“孔隙度”与“声波时差”、“中子”、“密度”进行多元线性拟合，得到孔隙度多元线性回归方程，砾岩类储层孔隙度计算模型如下：
[0039]POR＝a*AC+b*CNL+c*DEN+d
[0040]其中a＝0.050，b＝0.146，c＝
‑
20.327，d＝56.390，Multiple R＝0.904，MultipleR为多元相关系数，孔隙度平均绝对误差为1.68％，与传统体积模型计算孔隙度方法相比精度提高1.32％。
[0041]制作分析了该模型计算孔隙度与岩心实测孔隙度检验重叠图。即：分别制作了“孔隙度与声波时差”、“孔隙度与中子”、“孔隙度与密度”交会检验重叠图(见图1.1、图1.2、图1.3)，图中显示用模型计算孔隙度与岩心实测孔隙度重叠性高度一致，模型精度高。
[0042]实施例二：
[0043]砂砾岩体砂泥岩类储层孔隙度计算模型：
[0044]利用5口井139个砂砾岩体砂泥岩类储层样本，将样本实测“孔隙度”与“声波时差”、“中子”、“密度”进行多元线性拟合，得到孔隙度多元线性回归方程，砂泥岩类储层孔隙度计算模型如下：
[0045]POR＝a*AC+b*CNL+c*DEN+d，
[0046]其中a＝0.085，b＝
‑
0.246，c＝
‑
36.452，d＝56.390101.287，MultipleR＝0.828，Multiple R为多元相关系数，孔隙度平均绝对误差为1.本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种计算砂砾岩体储层孔隙度方法，其特征在于，包括以下步骤：S1、提出按照岩性类别建立砂砾岩体储层孔隙度模型，开展岩性类别划分；S2、提取对孔隙度敏感的测井响应曲线；S3、按砂砾岩体储层岩石类别建立岩心实测孔隙度与测井响应曲线数据表；S4、采用多元线性回归方法建立孔隙度计算模型，将其分为三个计算模型：(1)砂砾岩体砾岩类储层孔隙...

【专利技术属性】
技术研发人员：陈志湘，朱孟高，霍爱民，杨兆刚，吕维刚，
申请(专利权)人：扬州华元油气技术有限公司，
类型：发明
国别省市：