随着油气勘探的不断深入，珠江口盆地东部已经发现越来越多的粉细砂岩和钙质砂岩油藏^[1]，这类油藏通常表现为岩性复杂、薄互层的特征，储层的孔隙结构和导电性质存在明显的差异性，在宏观上表现为孔渗关系复杂、物性变化范围广、测井响应特征复杂、非均质性特征较强。测井解释建模的传统做法是分区分层建立模型^[2]，其前提是假设同一区块、同一层段的储层是相对均质的，这就忽略了储层内部孔隙结构和渗流特性的差异。如果采用这种方法来评价珠江口盆地惠州凹陷复杂砂岩储层，通常会导致其孔隙度、渗透率和饱和度出现较大误差，而这些参数是落实油藏有效性的重要基础。目前通过针对特定区域开展复杂储层岩石物理相研究，基于储层岩石物理相分类建立测井解释模型，能够快速有效地获得准确的储层参数，并应用于开展油藏精细描述^[3-5]。本文以珠江口盆地惠州凹陷22口取心井珠江组和珠海组复杂砂岩储层的岩心资料、测井资料和分析化验资料为基础，分析不同类型岩石物理相对储层岩性、物性和电性的控制作用，确定基于岩石物理相的储层分类方法，采用Fisher判别分析方法确定岩石物理相的识别标准；在此基础上，利用岩心资料刻度测井信息，建立不同类型储层对应的孔隙度、渗透率和岩电参数模型。

珠江口盆地惠州凹陷位于珠一坳陷中部，面积约为10000km²，是珠江口盆地一个重要的富烃凹陷，其沉积地层自下而上依次为始新统文昌组和恩平组，渐新统珠海组，中新统珠江组、韩江组和粤海组，上新统万山组（图 1）及第四系。惠州凹陷产层主要分布在珠江组和珠海组，其中珠江组为大套海相三角洲沉积，珠海组为典型的海陆过渡相沉积，主要发育三角洲平原、三角洲前缘及滨岸相^[6]。根据前人研究^[7]，研究区目的层主要分布在1800~3500m范围内，相当于早成岩阶段B期到中成岩阶段A期，总体上经历的成岩改造程度较弱。储层岩性以中—细、极细长石石英砂岩和中—细含灰质砂岩为主，颗粒间主要为线接触，孔隙类型以原生粒间孔为主，次生溶孔包括粒间溶孔和粒内溶孔。储层孔隙度主要为8%~28%，渗透率为0.1~1000mD，但受颗粒粗细、泥质含量和灰质胶结物的影响，目的层孔隙结构多样、孔渗关系复杂，不同类型储层的物性与饱和度计算模型和参数差异较大，给解释评价带来很大挑战。

图 1

图 1 珠江口盆地惠州凹陷构造位置（a）及纵向地层格架（b） Fig. 1 Structural location(a)and stratigraphic framework(b)of Huizhou Sag in the Pearl River Mouth Basin

研究区储层岩性较为复杂，物性分布范围广泛，从特低孔、特低渗到高孔、高渗均有涵盖，如图 2所示。从图中可以看出，尽管孔渗之间拟合的相关系数较高，但同一孔隙度对应的渗透率相差可达两个数量级，这说明岩石的渗流特性虽然与孔隙度密切相关，但无法用单一的孔渗关系来表征。岩心压汞曲线形态多样，孔隙结构参数分布范围广泛：岩心样品的排驱压力为0.01~1.37MPa，最大连通孔喉半径为0.54~61.87μm，平均孔喉半径为0.12~18.97μm，孔喉半径中值为0.01~16.61μm，孔喉分选系数为0.11~5.35μm，反映了储层类型的多样性和孔隙结构的复杂性。为了精确描述储层特征，为后续勘探开发提供可靠的基础数据，需开展储层分类研究。

图 2

图 2 惠州凹陷珠江组、珠海组岩心分析孔隙度与渗透率关系图 Fig. 2 Relationship between porosity and permeability by core analysis of Zhujiang & Zhuhai Formations in Huizhou Sag

国外学者Spain D R于1992年提出在单井剖面上划分岩石物理类型^[8]，国内学者代金友等指出岩石物理相是具有一定岩石物理特性的储层成因单元，是沉积作用、成岩作用及后期构造作用的综合效应，最终表现为现今的孔隙几何学特征^[9]。从岩石物理相的角度对储层进行分类评价，可以取得良好的应用效果^[10-12]。前期关于岩石物理相的划分主要是基于流动单元^[13]和模糊聚类^[14]的方法，但这些方法多基于某个单一参数或划分的岩石物理相不具有明确的地质意义；后期在成岩相、沉积相和裂缝相研究的基础上，采用综合分析和命名的方法开展岩石物理相分类研究^[15-17]。本次研究依据研究区实际资料，采用岩石粒度和矿物成分表征沉积微相，采用孔隙度、渗透率和孔隙结构参数表征岩石物理相，珠江组和珠海组砂岩储层的裂缝基本不发育，所以不考虑裂缝相的影响。

根据研究区储层地质特征，选取能够表征岩石粒度、矿物成分和物性的6项参数开展岩石物理相综合分类，具体包括粒度均值、钾长石含量、方解石含量、孔隙度、渗透率和孔隙结构特征参数R35（进汞饱和度为35%时对应的孔喉半径值）。在对22口取心井的316块岩心样品实验结果进行综合分析的基础上，将储层划分为如下4类岩石物理相。

Ⅰ类：长石石英中—细砂岩相，岩性以中砂岩、细砂岩为主，中砂岩、细砂岩分别占该类岩石物理相的69.57%和23.91%。矿物成分以石英和钾长石为主，石英和钾长石的平均含量分别为76.9%和14.1%，方解石和白云石含量较少。该类储层颗粒粒度偏粗、分选好、抗压实作用强、胶结程度弱；储层主要发育粒间孔，连通性好，压汞曲线为宽平台型，R35平均为22.4μm；孔隙度一般大于12%，主值区间孔隙度平均为16.9%，渗透率一般大于10mD，主值区间渗透率平均为72.6mD，属于中孔中渗储层，是该区最优质的储层。该类储层对应测井响应特征为“三低一高”，即低自然伽马、低密度、低光电吸收截面指数和高中子孔隙度。

Ⅱ类：含黏土长石石英细砂岩相，岩性以细砂岩为主，细砂岩占该类岩石物理相的80.95%。矿物成分以石英和钾长石为主，平均含量分别为67.9%和11.3%；黏土矿物含量较高，平均为16.2%；方解石和白云石含量较低。相对于Ⅰ类岩石物理相而言，粒度较细。储层孔隙发育中等、连通性较好，压汞曲线为缓坡型，R35平均为2.59μm；孔隙度一般大于10%，主值区间孔隙度平均为16.3%，渗透率一般大于1mD，主值区间渗透率平均为21.6mD，储层以中孔低渗为主。该类储层对应测井响应特征为“三中一高”，即中自然伽马、中密度、中光电吸收截面指数、高中子孔隙度。

Ⅲ类：含黏土长石石英粉砂岩相，以粉砂岩、细砂岩为主，其中粉砂岩占75%。矿物成分以石英、斜长石和钾长石为主，平均含量分别为65.6%、7.1%和8.9%；黏土矿物含量为17.9%，方解石和白云石含量较低。相较于Ⅱ类岩相而言，粒度更细、泥质含量更高，储层更容易受到压实和胶结作用的改造，孔隙连通性较差，压汞曲线为斜坡型，R35平均为1.47μm，主值区间孔隙度平均为15.4%，主值区间渗透率平均为1.95mD，属于中孔特低渗储层。该类储层对应测井响应特征为“二高二中”，即高自然伽马、高中子孔隙度、中密度和中光电吸收截面指数。

Ⅳ类：钙质石英中—细砂岩相，以钙质中砂岩、钙质细砂岩为主。矿物成分以石英和方解石为主，平均含量分别为52.9%和23.2%；含有少量的长石和白云石。该类储层分选差，碳酸盐胶结程度较为严重，孔隙连通性很差；压汞曲线为斜坡型，R35平均为0.57μm；主值区间孔隙度平均为14.2%，主值区间渗透率平均为1.07mD，属于低孔特低渗储层，是研究区品质最差的储层。该类储层对应测井响应特征为“两低两高”，即低自然伽马、低中子孔隙度、高密度、高光电吸收截面指数。

从上述岩石物理相分类的描述可以看出，岩心资料揭示研究区储层品质主要受控于粒度、钾长石含量和方解石含量。在矿物成分分布相近的前提下，颗粒越粗，抗压实作用越强，孔喉连通性越好；钾长石在一定程度上能改善储层内部的渗流特性，其含量越高，储层渗透性越好；而作为胶结物形式出现的方解石，其含量越高，胶结程度越强，储层的渗透性和导电性也越差。

由于岩心获取难度大、分布不连续、分析化验成本高等因素，如何从测井资料实现对该区储层的岩石物理相准确分类就成为亟待解决的问题。考虑到目的层的测井系列，选取对粒度、钾长石含量和方解石含量较为敏感的自然伽马、密度、中子孔隙度与光电吸收截面指数4条测井曲线，以岩心分析数据作为标定，研究建立基于测井曲线的岩石物理相分类方法。

研究区测井系列主要采用ECLIPS-5700系统，大部分井主要为21世纪初期完成资料采集，不同井之间的测井曲线系统误差不大，因此不用考虑资料标准化。但是，不同测井曲线的量纲、数值分布区间及受孔隙流体的影响差异必须考虑，研究采用规则化资料预处理方法，提取主要反映储层骨架特征的参数^[18]，包括正则化自然伽马（NGR）、密度孔隙度与中子孔隙度差值（Δϕ）及骨架岩性指数（P）3个参数。利用岩心标定后的4类岩石物理相对应的3个测井参数响应图版见图 3，具体数值范围见表 1，其中：

图 3

图 3 惠州凹陷珠江组、珠海组储层4类岩石物理相测井响应特征参数图版 Fig. 3 Plate of logging response characteristic parameters of the four types of petrophysical facies of Zhujiang & Zhuhai Formations in Huizhou Sag

表 1

表 1 惠州凹陷珠江组、珠海组储层4类岩石物理相测井响应特征参数统计表 Table 1 Statistical table of logging response characteristic parameters of the four types of petrophysical facies of Zhujiang & Zhuhai Formations in Huizhou Sag 岩石物理相类型 分布 GR/API DEN/ (g·cm-3) ϕN Pe/ (bar·e-1) NGR Δϕ P/ (bar·g-1) Ⅰ类 范围 45.0~85.6 2.18~2.42 0.06~0.30 2.0~2.9 0.45~0.89 -0.33~0.57 0.70~1.18 平均值 66.2 2.28 0.21 2.4 0.62 0.26 0.87 Ⅱ类 范围 80.8~119.9 2.22~2.47 0.14~0.27 2.0~3.6 0.33~0.50 -0.58~0.28 0.69~1.21 平均值 102.2 2.36 0.22 2.8 0.40 -0.11 0.92 Ⅲ类 范围 107.0~165.9 2.28~2.51 0.17~0.34 2.5~3.5 0.24~0.37 -1.07~-0.03 0.87~1.19 平均值 132.3 2.39 0.26 3.1 0.31 -0.52 1.05 Ⅳ类 范围 47.0~94.6 2.36~2.58 0.15~0.24 2.9~4.7 0.42~0.85 -0.52~-0.04 1.01~1.56 平均值 64.9 2.46 0.17 3.8 0.63 -0.24 1.28

表 1 惠州凹陷珠江组、珠海组储层4类岩石物理相测井响应特征参数统计表 Table 1 Statistical table of logging response characteristic parameters of the four types of petrophysical facies of Zhujiang & Zhuhai Formations in Huizhou Sag

式中  GR——自然伽马，API；

GR_min——自然伽马最小值，API；

ϕ_D——石灰岩刻度的视密度孔隙度；

ϕ_N——中子孔隙度；

U——岩石体积光电吸收截面指数，bar/cm³；

U_f——流体体积光电吸收截面指数，取0.36bar/

cm³；

DEN——测井密度，g/cm³；

DEN_f——流体密度，取1.0g/cm³。

进一步利用Fisher判别分析方法构建不同岩石物理相类型识别函数，其中第1、第2和第3典型判别函数分别为F₁、F₂和F₃，表达式如下：

典型函数的特征值和贡献率如表 2，其中第1和第2典型函数的贡献率分别为63.6%和36.1%，累计达到99.7%，基本包含绝大部分的测井信息，故选择第1和第2典型函数作为岩石物理相类型判别的特征变量。将第1和第2典型函数做交会图（图 4），4类岩石物理相之间的界限明显，判别精度较高。

表 2

表 2 惠州凹陷珠江组、珠海组储层典型岩石物理相Fisher函数判别特征值和贡献率分布特征 Table 2 Fisher function eigenvalue and distribution characteristics of the contribution ratio of typical petrophysical facies of the Zhujiang & Zhuhai Formations in Huizhou Sag 典型函数 特征值 贡献率/ % 累计贡献率/ % 相关性 F1 2.886 63.6 63.6 0.862 F2 1.636 36.1 99.7 0.788 F3 0.016 0.3 100.0 0.125

表 2 惠州凹陷珠江组、珠海组储层典型岩石物理相Fisher函数判别特征值和贡献率分布特征 Table 2 Fisher function eigenvalue and distribution characteristics of the contribution ratio of typical petrophysical facies of the Zhujiang & Zhuhai Formations in Huizhou Sag

图 4

图 4 惠州凹陷珠江组、珠海组储层岩石物理相Fisher判别结果图 Fig. 4 Fisher discrimination graph of petrophysical facies of Zhujiang & Zhuhai Formations inHuizhou Sag

通过上述Fisher降维法得到的典型函数F（X），根据最小平方误差准则函数的梯度下降算法^[19]，经转化可得到4类储层的判别函数f：

Ⅰ类长石石英中—细砂岩相：

Ⅱ类含黏土长石石英细砂岩相：

Ⅲ类含黏土长石石英粉砂岩相：

Ⅳ类钙质石英中—细砂岩相：

对于未知储层，把其对应的测井参数代入到公式（7）~（10）分别计算其f值，判别函数最大值所属的类别即为该储层的岩石物理相类型。

利用上述判别函数对22口取心井的157个小层进行岩石物理相类型判识，符合率为92.4%（表 3）。统计结果表明不同岩石物理相时空分布大致具有如下规律：Ⅰ类和Ⅱ类岩石物理相主要分布在珠江组下段和珠海组上段，平均厚度分别为196m、348.6m。Ⅲ类岩石物理相主要分布在珠江组上段，平均厚度为161.2m。Ⅳ类相主要分布在珠江组下段，平均厚度为48.8m。

表 3

表 3 惠州凹陷珠江组、珠海组储层岩石物理相类型Fisher判别结果统计表 Table 3 Statistical table of Fisher discrimination results of petrophysical facies types of Zhujiang & Zhuhai Formations in Huizhou Sag 岩石物理相类型 Fisher判别分析识别结果 符合率/% 符合 不符合 Ⅰ类 73 3 96.1 Ⅱ类 29 5 85.3 Ⅲ类 22 4 84.6 Ⅳ类 21 0 100.0 合计 145 12 92.4

表 3 惠州凹陷珠江组、珠海组储层岩石物理相类型Fisher判别结果统计表 Table 3 Statistical table of Fisher discrimination results of petrophysical facies types of Zhujiang & Zhuhai Formations in Huizhou Sag

如前所述，研究区不同类型储层岩石物理特征差异较大，在准确分类的基础上，可以充分利用岩心实验数据针对不同岩石物理相建立对应的储层参数精细解释模型。

首先考虑孔隙度建模。由于岩心分析孔隙度与测井孔隙度曲线的纵向分辨率不同，直接利用岩心孔隙度数据与测井曲线值建立相关关系会带来误差。为提高测井解释模型的精度，需要对岩心分析数据和测井数据进行纵向匹配，本文采用划分样本层的方法来实现两者分辨率的匹配。

将研究区22口取心井的2204个岩心孔隙度样本点划分为453个样本小层，通过分析各样本小层岩心孔隙度与密度孔隙度、中子孔隙度和自然伽马等测井参数的相关性，采用多参数回归的方法，分别建立4类岩石物理相对应的孔隙度解释模型（表 4）。

表 4

表 4 惠州凹陷珠江组、珠海组储层不同岩石物理相测井解释参数模型及参数表 Table 4 Statistical table of log interpretation models and parameters for different petrophysical facies of Zhujiang & Zhuhai Formations in Huizhou Sag 岩石物理相类型 孔隙度计算模型 覆压孔渗关系 岩电参数 a b m n Ⅰ类 ϕ=0.4397ϕD+0.3717ϕN-0.0933ΔGR+0.9392 K=0.025e47.96ϕ 1 1 1.704 1.867 Ⅱ类 ϕ=1.0362ϕD-0.0617ϕN-0.0178 K=108ϕ8.509 1 1 1.704 1.867 Ⅲ类 ϕ=0.5823ϕD+0.4321ϕN-0.044ΔGR-0.03276 K=0.00006e63.62ϕ 1 1 1.515 2.013 Ⅳ类 ϕ=0.4496ϕD+0.4428ϕN-0.0105 1 1 1.828 1.940 注：ϕ—岩石孔隙度；K—岩石渗透率，mD；ΔGR—自然伽马相对值；ΔGR=（GR-GRmin）/（GRmax-GRmin）；GRmin—自然伽马最小值，API；GRmax—自然伽马最大值，API；a—岩性系数；b—系数；m—胶结指数；n—饱和度指数。

表 4 惠州凹陷珠江组、珠海组储层不同岩石物理相测井解释参数模型及参数表 Table 4 Statistical table of log interpretation models and parameters for different petrophysical facies of Zhujiang & Zhuhai Formations in Huizhou Sag

研究区的孔渗关系难以用统一的模型进行表述，岩石物理相的划分则为测井渗透率建模提供了框架模型^[20]。对图 2所示的孔渗数据进行岩石物理相分类，发现4类岩石物理相表现出较为明显的3种孔渗关系（图 5）。从图 5可以看出，Ⅲ类、Ⅳ类岩石物理相受泥质和钙质胶结的影响，孔隙结构明显变差，在相同的孔隙度条件下，其渗透率明显小于Ⅰ类、Ⅱ类岩石物理相。采用统计回归分析的方法分别建立4类岩石物理相对应的3种渗透率模型，如表 4所示。

图 5

图 5 惠州凹陷不同类型岩石物理相渗透率模型 Fig. 5 Permeability interpretation models for different types of petrophysical facies in Huizhou Sag

同样，储层岩电参数主要受控于岩性、孔隙结构、黏土矿物类型和分布形式^[21-24]。对于孔隙结构差异大的储层，如果采用相同的岩电参数，往往会造成含水饱和度的计算误差较大，影响流体类型判识和储量计算^[25-28]。通过分析研究区14口井137块岩心岩电实验数据，发现4类岩石物理相的地层因素与孔隙度、电阻率增大指数与含水饱和度在双对数坐标系中表现出3种关系（图 6）。从图 6可以看出，Ⅲ类岩石物理相受到黏土矿物含量较高和颗粒较细的影响，m值明显变小；Ⅳ类岩石物理相受钙质胶结的影响，岩石的导电路径更加复杂，m值明显变大。采用回归方法确定其岩电参数，如表 4所示。

图 6

图 6 惠州凹陷珠江组、珠海组储层不同岩石物理相岩电参数图版 Fig. 6 Plate of rock electric parameters of different petrophysical facies of Zhujiang & Zhuhai Formations in Huizhou Sag (a)—(c)依次为Ⅰ类 & Ⅱ类、Ⅲ类和Ⅳ类储层的地层因素与孔隙度关系图；(d)—(f)依次为Ⅰ类 & Ⅱ类、Ⅲ类和Ⅳ类储层的电阻率指数与含水饱和度关系图

图 7为惠州凹陷A井ZJ1-05和ZJ1-06层评价成果图。图中第11道为利用本文的方法进行岩石物理相分类的结果；第6、7、8道中的红色曲线分别为分类型计算的渗透率、含水饱和度和孔隙度曲线，蓝色曲线为原始单一模型的计算成果，离散点为密闭取心分析孔隙度、渗透率和含水饱和度数据。薄片分析数据表明该井ZJ1-05层1927~1932m井段以含泥中—细粒长石石英砂岩为主，孔隙空间主要为原生粒间孔，储层物性较好，以Ⅱ类岩石物理相为主；ZJ1-06层岩性以含泥极细粒长石石英砂岩为主，局部为含钙细粒石英砂岩，物性较差，以Ⅲ类岩石物理相为主，局部为Ⅳ类岩石物理相，与测井识别的岩石物理相类型一致。

图 7

图 7 惠州油区A井珠江组ZJ1-05层与ZJ1-06层储层分类评价成果图 Fig. 7 Reservoir classification and evaluation results graph of ZJ1-05 and ZJ1-06 layers of Well A

由图 7可知，采用岩石物理相分类计算储层参数的方法明显提高了孔隙度、渗透率、饱和度的计算精度。与岩心数据对比，分类建模后计算孔隙度的相对误差由单一模型的22.85%降低至5.58%，含水饱和度绝对误差由11.24%降低至4.29%，渗透率的相对误差减小至半个数量级以内。分类解释后将ZJ1-06层由单一模型解释的底水油藏调整为边水油藏，油水界面下移5.7m，DST测试获得日产纯油96m³，且稳定生产16个月，含水率保持在5%以下，与分类解释结论一致。

图 8为B井ZJ2-03层评价成果图，在2356.8~2357.9m井段，岩心分析渗透率小于0.5mD，按照区域标准该段应为干层。若采用单一的渗透率模型和岩电参数进行计算，其渗透率和换算后的含油饱和度被明显高估，从而解释为油层。从岩石物理相分类结果可以看出，该深度段为Ⅳ类钙质砂岩相。通过岩石物理相分类处理，有效消除了储层含钙对渗透性及导电性的影响，各项参数的计算精度明显提高。扣除非渗透性钙质夹层，ZJ2-03层的有效厚度减少了1.9m。

图 8

图 8 惠州油区B井珠江组ZJ2-03层储层分类评价成果图 Fig. 8 Reservoir classification and evaluation results graph of ZJ2-03 layer of Well B

依据上述成果对研究区65口老井进行重新处理，在其中20口井中新解释油层73层，累计有效厚度为187.5m，为老油田挖潜及开发评价提供了有效指导，生产实效显著。

（1）珠江口盆地惠州凹陷珠江组、珠海组砂岩储层特征差异主要受沉积和成岩作用影响，根据岩心资料，可将储层划分为4类岩石物理相，分类结果显示粒度均值和矿物成分是影响该区储层岩石物理响应的关键因素，通过构建正则化自然伽马、密度孔隙度与中子孔隙度差值及骨架岩性指数3项测井衍生参数可有效识别4种不同岩石物理相。

（2）不同的岩石物理相对储层物性及导电性有明显的控制作用，利用实验资料标定，针对不同岩石物理相分别构建孔隙度、渗透率、饱和度计算模型，相较于以往常用的单一解释模型，储层参数计算精度得到显著提高，从而为老油田挖潜及开发评价提供有效指导。

（3）本文采用的岩石物理相分类方法具有明确的地质意义，有效建立起地质特征与测井属性之间的“桥梁”，克服了常规分类方法地质意义模糊或与测井参数关联度低的不足，可推广应用于其他复杂碎屑岩储层测井评价。