基于常规测井的变质岩裂缝识别方法研究

廖海博，刘红岐，崔云江，汪瑞宏，李枝林，陈东，刘伟

1. 西南石油大学地球科学与技术学院，成都 610500

2. 中海石油（中国）有限公司天津分公司，天津 300459

3. 川庆钻探钻采工程技术研究院，广汉 618300

渤海湾含油气盆地是继大庆油田发现之后，在中国东部地区所发现的另一个重要的含油气盆地，可分为碎屑岩、古潜山变质岩、湖相碳酸盐岩、火山岩等4 大岩石类型。其中，古潜山变质岩的勘探程度较低，已探明储量仅占总储量的20%[1]。因此，古潜山油气藏是渤海湾盆地主要的勘探目标[2]。BZ19-6 气田是渤海湾盆地最大的凝析气田，探明储量超千亿方[3]，主要目的层为太古界潜山地层，岩性主要为花岗片麻岩，夹少量侵入岩[4-5]。受印支期、燕山期和喜马拉雅期等多期次构造运动的影响，太古界变质花岗岩潜山发育溶蚀孔隙-裂缝型和裂缝型储层，其中裂缝是主要的储集空间和渗流通道[6-7]。因此，有效的裂缝识别方法对该区块未探明单元的油气勘探至关重要。

目前，识别储层裂缝的方法有直接法和间接法。直接法是从露头、岩芯和成像测井中识别裂缝特征，可以直接反映储层裂缝，但其局限性强、价格昂贵、操作不易。间接法通过处理和分析常规测井资料、地震资料和动态生产资料识别裂缝，可以根据少量数据获得大面积内的裂缝特征，但其必须引入复杂的数学算法。近年来，随着计算机技术的迅速发展和数学理论基础的深入研究，许多学者将常规测井与数学方法相结合，间接地识别裂缝发育的位置和程度。例如，Xiao 等[8]利用重标极差（R/S）分析法对不同岩性进行了裂缝识别，表明运用R/S分析方法识别储层的裂缝具有可行性；陈义国等[9]采用R/S分析方法对常规测井资料进行裂缝识别，发现裂缝识别参数与裂缝线密度呈良好的线性关系；Li 等[10]通过R/S分析方法对北贵州地区下寒武统牛蹄塘组页岩储层的裂缝进行识别，表明R/S分析方法可以有效地区分裂缝与噪声信号；Zhang 等[11]基于R/S分析方法识别裂缝，表明可通过计算裂缝与基质之间的比例来确定裂缝线密度，进而预测天然裂缝的发育程度和空间分布规律；Xiao 等[12]利用R/S分析方法对致密砂岩储层的非均质性进行评估，认为R/S分析方法可定量评估储层的非均质性；Aghli 等[13]基于R/S分析方法识别碳酸盐岩储层裂缝，表明R/S分析方法能定量反映裂缝参数；Ge 等[14]探讨了火成岩多重分形参数与裂缝发育程度之间的关系，发现裂缝发育程度与多重分形特征呈负相关。研究成果显示，重标极差（R/S）分析与分形维数已广泛应用于天然裂缝的识别，但仍存在以下问题：①识别裂缝的尺度具有局限性；②对噪声和异常值较敏感，需要通过其他地质、地球物理等资料进行验证。

将此方法用于该区块的裂缝识别，仍面临着巨大的挑战。首先，BZ19-6 气田储层裂缝网络复杂，主要发育微裂缝，裂缝宽度集中分布于0.01～0.2 mm，且大多数裂缝被方解石、泥质等填充；其次，针对岩性更为复杂的花岗片麻岩储层，R/S分析方法识别裂缝的研究鲜见。鉴于此，本文基于R/S分析方法并结合测井、岩芯和薄片等资料，旨在：①研究R/S分析方法识别裂缝的最小尺度；②建立利用赫尔特指数识别花岗片麻岩储层裂缝发育程度分类标准；③分析岩性各向异性对R/S分析方法识别裂缝精度的影响。

1 地质背景

BZ19-6 气田位于渤中凹陷西南部，东南方向为渤南低凸起，西部为埕北低凸起，呈洼中隆的构造格局（图1）。如图2 所示（图1 中DD′剖面），太古界潜山地层在近南北向发育的构造脊上是古隆起的背景下发育的、经过多期次的构造作用改造的、被一系列断层复杂化的断层背斜构造，构造整体南高北低[15-16]，形态完整，圈闭面积较大，埋藏较深，次生断层、裂缝较为发育[17]。

图1 研究区位置图Fig.1 Location of study area

图2 太古界潜山气藏构造剖面图Fig.2 Archean buried hill gas reservoir structure profile

研究区太古界潜山与上部始新统孔店组呈不整合接触，潜山岩性主要为变质岩，局部可见闪长玢岩、辉绿岩等侵入岩，纵向上可划分为半风化带和潜山内部[18]。半风化带岩石颜色较浅，主要为石英及长石，少量暗色矿物，其中部分岩石呈灰黑色，发育孔隙-裂缝型和裂缝-孔隙型储层[19]；半风化带上部溶蚀作用较强，缝孔洞较为发育，风化比较严重，长石大多已经风化为高岭土；而在下部，风化程度逐渐变弱。潜山内部则主要发育裂缝型储层和孔隙-裂缝型，溶蚀孔、洞较少或者不发育。潜山储层主要发育三期裂缝[20]，第一期裂缝主要发育在印支期，受华北板块和扬子板块碰撞而产生的大量挤压构造裂缝；第二期裂缝主要发育在燕山期，郯庐断裂发生左旋走滑作用，导致岩石破碎作用强烈，形成大量碎裂岩和碎斑岩以及一系列动力破碎带，并派生出大量裂缝；第三期裂缝与新近纪太平洋俯冲和郯庐断裂发生的右旋走滑作用有关[21]。

2 裂缝特征

2.1 岩芯薄片裂缝特征

天然裂缝按成因可分为构造裂缝、岩溶裂缝、热胀冷缩裂缝、干缩裂缝等[22]，按充填情况可分为张开缝和充填缝两类[23]。通过对BZ19-6 区块的112 块岩芯的观察（图3 中a、b 是原始岩芯，c、d 为其岩芯切面），岩性为花岗片麻岩，裂缝较发育，多见微缝（＜0.1 mm）及小缝（0.1～0.5 mm），存在少量中缝（0.5～10 mm）。裂缝产状主要为中高角度，有少量的垂直缝和水平缝，裂缝多呈网状交叉，大多被充填-半充填。

图3 太古界天然裂缝岩芯图片Fig.3 Core pictures of Archean natural fractures

观察201 份岩芯薄片后，发现薄片裂缝较为发育，主要为微裂缝和小缝，大多处于充填-半充填状态。裂缝主要被云母化泥质、铁质、碳酸盐岩和黄铁矿充填，但也可见少量未充填的裂缝。其中，部分薄片显示斜长石绢云母化普遍，钾长石高岭土化，常见铁白云石交代现象及被铁质或铁泥质充填的岩石微裂缝（图4a）。另外，部分薄片显示黑云母见后期白云石交代现象，见微小裂缝，宽度为0.01～0.15 mm，其内被白云石、铁白云石和方解石所致密充填（图4b、c）。此外，部分岩石内裂缝和颗粒裂隙也较为发育，缝宽为0.01～0.20 mm，见被黄铁矿充填的微裂缝，但也可见少量未充填的裂缝（图4d）。另有部分岩石内发育数条裂缝，缝宽为0.01～0.17 mm，部分被碳酸盐矿物充填，也可见未充填裂缝，缝宽约为0.01～0.03 mm（图4e）。还有部分岩石见两期裂缝发育，缝宽为0.02～0.05 mm，裂缝先被石英、白云石、方解石完全充填，后被铁质或铁白云石充填，微裂隙较为发育，微裂隙宽度为0.01 mm 左右（图4f）。

图4 太古界岩石薄片Fig.4 Rock slices of Archean strata

2.2 随钻成像测井裂缝特征

测井资料处理解释过程中将裂缝分为天然裂缝和诱导缝。其中，天然裂缝可进一步分为高阻缝和高导缝，按裂缝的倾角又可分为高角度裂缝、中角度裂缝和低角度裂缝。为了进一步了解该区块的裂缝发育情况及特征，对该井进行了Microscope-HD 随钻成像测井资料处理分析。结果显示，该区块发育高导缝（图5a）、高阻缝（图5b）、诱导缝（图5c）及线理（图5d）。其中，高导缝较为发育，走向主要有NW-SE 和近E-W 向（图5e），以中高角度为主，倾角集中在40°～70°，均值为48°左右（图5f）。裂缝宽度集中分布在0.001～0.1 mm，其中＜0.1 mm 的裂缝占比91.12%（图5g）；裂缝密度为1～13 条/m，1～8 条/m 占所有的96.22%，平均为2.3～4.32 条/m（图5h）。区块高阻缝发育较少，走向为NNW-SSE，以中等角度为主。电成像测井只能识别到毫米级别的裂缝，埋深和充填情况都会降低其分辨率；而岩芯的分辨率可以达到亚毫米级别，可直接观察到岩石的构造和裂缝，但岩芯取样点不连续且难以覆盖整个井段。两者在识别裂缝的尺度上不统一，均难以识别研究区的微裂缝。因此，亟需一种更为精细、有效的裂缝评价方法。

3 研究方法

R/S分析，也称为重标极差分析，由Hurst 在1951 年提出，用于研究尼罗河水流量与库容之间的关系[24]，后被众多学者证明它是分析一维分形变量的有效技术[25-27]。对于单一岩性的地层，天然裂缝的发育给测井所测得的物理性质带来很大的变化，增强了地层的非均质性，裂缝发育程度不同，其物理性质的改变情况有所不同[28]。因此，可以根据测井曲线上的突变点大致判断裂缝的位置，根据突变点的强弱分析裂缝的发育程度。

通过分析某些测井参数的R(n)/S(n)（后文缩写为R/S）与测井采样点n的关系，预测裂缝的发育程度。如公式(1)、(2)、(3)所示，R(n)是最大累积偏差与最小累积偏差之差，代表测井序列的复杂性；S(n)是变化的平方根，代表时间序列的平均趋势；Z(i)是测井序列，每一个采样间隔点上的值；i、j、k是从0 到n变化的常数，其中k大于等于i，H表示赫尔特指数。

通过数值模拟，计算出R和S的值，做Lg(R/S)与Lg(n)的散点图，其斜率就代表赫尔特指数，H的范围是0 到2，可以通过H的大小判断裂缝的发育程度。再对Lg(R/S)的值进行牛顿差分运算，n 阶差分如公式(4)所示：

对该区块12 口井进行计算后发现，二阶差分模型对该区块的裂缝位置预测准确性较高，其模型如公式(5)所示：

4 结果与讨论

该区块随钻测井曲线包含井径曲线（CAL）、自然伽马曲线（GR）、冲洗带电阻率曲线（Rxo）、原状地层电阻率曲线（Rt）、密度曲线（RHOB）、中子曲线（TNPH）、声波时差曲线（DT）等，通过岩芯裂缝统计对比，选择对天然裂缝较为敏感的测井参数。TNPH 曲线受岩性、钻井液等因素干扰严重；GR 曲线主要反映岩性的变化，而泥质含量的变化和钻井液的变化会对GR 值影响明显，从而掩盖或放大裂缝存在的信息，而花岗片麻岩地层较为致密，密度变化较小。相比之下，CAL、Rxo、DT 三条测井曲线对裂缝更为敏感，而GR、RHOB、TNPH 对裂缝的敏感度较低（图6）。通过实际处理，发现在Lg(n)=3.5左右的深度为4 882 m 处，R/S曲线出现明显的下凹区间，表明测井响应在此处有较大的变化，结合电成像图显示，此段裂缝高度发育。因此，本文选择CAL、Rxo、DT 三条测井曲线进行处理分析，以进一步分类识别该井的裂缝发育程度。

图6 R/S 分析曲线Fig.6 R/S analysis curve

选取30 个层段共150 m，经过R/S分析发现，在花岗片麻岩中利用R/S分析方法能够有效识别裂缝的发育情况。该方法的准确性也得到了电成像资料的验证。针对这30 个层段的CAL、Rxo、DT 三条测井曲线的R/S分析结果，分别用HCAL、HRxo、HDT表示，并绘制了三维散点图（图7），发现分区性明显。因此，本文建立了针对花岗片麻岩储层的裂缝发育程度分类识别标准（表1）。

表1 花岗片麻岩储层裂缝发育程度分类识别标准Table 1 Classification and identification standard of fracture development degree of granite gneiss reservoir

图7 赫尔特指数三维散点图Fig.7 Three-dimensional scatter plot of Hurst index

对R/S分析处理后的值采用牛顿差分的方法处理，经过对比后发现二阶差分对裂缝的位置识别较为准确，其中CAL、Rxo、DT 曲线的二阶差分值分别用K-CAL、K-Rxo、K-DT 表示。如图8 所示，对电阻率Rxo 曲线进行R/S分析，进一步采用牛顿差分的方法处理，得出以下结论：①定义R/S曲线明显下凹的部分为下凹区间，其二阶差分的值为正且大于0.0001，主要反映测井曲线的变化情况，对应裂缝发育段；②定义R/S曲线明显上凸的部分为上凸区间，其二阶差分的值为负，主要反映基岩段，其裂缝不发育或发育程度低。由式(5)计算的二阶差分K值与测井曲线分析的R/S曲线的下凹区间吻合较好，证明K值对识别曲线的下凹区间是有效的，进一步证明了识别裂缝发育位置的准确性。

图8 Lg(R/S)曲线的K 值和凹区间对比Fig.8 Comparison of K value and concave interval of Lg(R/S)curve

如表2 显示，通过牛顿二阶差分K值与随钻Microscope-HD 电成像图共识别出3 个有效的区域类别：①井壁崩落区：K-CAL 值较大，K-DT 和KRxo 值较小，因此K-CAL 值可以明显预测井壁崩落的位置；②裂缝发育区：K-DT 敏感程度强，从图像上也可以发现，裂缝发育的小层段，K-DT 值都比较突出；③基岩区：K-Rxo 值突出，成像图上显示特征为高阻亮色，裂缝不发育。对每一个区域类别的深入对比分析，可以得出以下结论：①K-CAL 值可以区分井壁崩落区和裂缝发育区，井壁崩落区KCAL＞0.00003，裂缝发育区K-CAL＜0.00003；②KDT 值对裂缝的发育位置识别较为准确，即裂缝越发育，K-DT 越大；③K-Rxo 值可以区分裂缝发育区和基岩区，裂缝发育区K-Rxo＜0.00001，基岩区KRxo＞0.00001。

表2 牛顿二阶差分K 值与Microscope-HD 电成像图对比Table 2 Comparison of Newton ' s second-order difference K value and Microscope-HD electrical imaging map

在裂缝发育区（K-Rxo＜0.00001），为进一步识别裂缝的发育程度，除了应用前文的赫尔特指数H值之外，电成像解释的裂缝线密度与K-Rxo 散点图 （图9）显示，裂缝线密度与K-Rxo 值呈正相关，即裂缝线密度越大，K-Rxo 值越大，其相关性R2=0.8895。因此，在裂缝发育区（K-Rxo＜0.00001），KRxo 值越大，裂缝线密度越大，裂缝越发育。

图9 裂缝线密度与K-Rxo 散点图Fig.9 Fracture line density and K-Rxo scatter plot

裂缝是油气井初周日产量的主要影响因素之一，初周的日产气和日产液可间接验证研究区裂缝的发育程度与分形维数的预测效果。选取研究区裂缝较为发育的层段，其垂直深度在4 650～4 700 m，在研究区井位平面图上绘制裂缝发育层段DT 曲线R/S分析曲线二阶差分K值的等值图，结合初周日产量与电成像图等资料，对BZ19-6 气田的裂缝相对发育情况进行预测（图10）。由图可知，A7H 井裂缝较发育，A3H 井裂缝发育程度低。在其他各种因素相近的条件下，牛顿差分K值与油气井初周日产液产量间存在正相关性，裂缝相对发育区的产量明显高于裂缝相对不发育区。

图10 BZ19-6 气田裂缝发育预测Fig.10 Prediction of Fracture Development in BZ19-6 Gas Field

计算结果表明，应用R/S分析结合二阶差分的方法进行裂缝识别和预测不仅具备坚实的理论基础，而且计算结果与成像测井所解释的裂缝信息具有较好的吻合度。客观地分析此方法识别裂缝精度存在误差的原因主要有以下3 个方面：一是测井信息不仅包括裂缝信息，还包括母岩本身的纵向非均质性信息（如岩性、物性、电性、含水饱和度等），其势必会给裂缝识别带来一定的影响；且岩性从花岗片麻岩变为闪长玢岩（图11），而又有裂缝发育，其中K-DT 值与K-Rxo 值较小，因此对非花岗片麻岩储层中的裂缝K值预测效果欠佳。二是研究区天然裂缝中充填方解石、白云石、石英、黄铁矿和泥质等矿物，当裂缝被矿物充填时，测井响应特征不明显，因此对完全充填缝的识别较难。三是该区块的裂缝宽度集中分布在0.001～0.1 mm，大量数据表明此方法对裂缝宽度小于0.005 mm 的微裂缝仍然有响应，但特征不明显，因此裂缝的宽度对识别效果有影响。

图11 非花岗片麻岩裂缝发育特征Fig.11 Fracture development characteristics of non-granitic gneiss

5 结论

（1）将R/S分析和牛顿差分法相结合改进的裂缝识别方法，在变质岩储层裂缝评价中具有可行性，并建立了利用赫尔特指数识别花岗片麻岩储层裂缝发育程度分类标准。

（2）此方法可识别宽度大于0.005 mm 的裂缝，提高了裂缝识别的精度。此外，研究还发现，K-Rxo值与裂缝线密度呈正相关关系，且相关性较好。

（3）岩性各向异性和裂缝充填情况影响常规测井曲线R/S分析方法识别裂缝的精度。