PNN测井清污混注水淹层剩余油饱和度计算方法

2021-06-08朱学娟单沙沙殷梓原孔雪

物探与化探 2021年3期

朱学娟，单沙沙，殷梓原，孔雪

(1.中国石油大学胜利学院油气工程学院，山东东营 257061； 2.中国石油集团测井有限公司油气评价中心，陕西西安 710077)

0 引言

传统水淹层定量评价主要依靠电法测井完成[1-2]，储层流体性质的识别难度大，测井解释符合率低，水淹层的定量评价及水淹等级划分存在很大困难[3-4]。过套管脉冲中子—中子测井(PNN)是建立在热中子俘获基础上的套管井剩余油饱和度测井方法，通过记录地层中未被地层俘获的热中子数量来确定储层的流体性质，计算油水饱和度。此方法不需要考虑伽马本底的影响，而且对于低孔隙度与低矿化度地层水的储层仍有较高的计数率，减小了统计误差[5-7]。但是测量条件、非地质因素影响等对地层热中子宏观俘获截面的提取增加了困难，且在计算饱和度过程中需要解释参数受岩性、流体性质、泥质、注入水等多种因素影响，所以，准确确定解释参数是利用PNN测井计算剩余油饱和度的前提[8]。

1 PNN测井原理

1.1 热中子俘获理论基础

原子核微观俘获截面指一个原子核俘获热中子的几率，单位为靶恩(b)，1 b=10-24cm2，岩石宏观俘获截面即1 cm3岩石中所有原子核微观俘获截面的总和单位为俘获单位(c.u.)，1 c.u.=10-3cm-1。岩石常见元素中氯的微观俘获截面最大，为31.6 b，而氢、碳、氧、镁、铝、硅、钙分别是0.329、0.004 5、0.001 6、0.4、0.215、0.13、0.43 b，石英、方解石、白云石和淡水的宏观俘获截面分别是0.003 4、0.007 1、0.004 6、0.022 cm-1，因此，岩石对热中子的俘获能力主要取决于含氯量，即地层含水量及地层水含盐量[9]。以热中子俘获为基础的测井方法通过测量地层的宏观俘获截面来研究地层性质，特别是储层流体性质，是在套管井中代替电阻率测井用来区分油气、水和研究开发动态的很好的测井方法[10]。

1.2 PNN测井原理及优势分析

中子寿命测井(热中子衰减时间测井)是最常用的在套管井中代替电阻率测井研究地层流体性质的测井方法。但中子寿命测井测量的是俘获伽马的计数率，为降低统计误差并消除地层的伽马本底，要求计数率越高越好，所以对中到高孔隙度、地层水矿化度较大的储集层效果最好。PNN测井通过发射较长时间间隔(75 ms)的中子脉冲与地层发生反应，并测量地层中未被俘获的热中子数量来计算地层对热中子的宏观俘获截面，其整个计数区间内的计数率都能反应热中子的衰减规律[11]。它的优势在于可有效避免自然伽马本底的影响，对于较低孔隙度、较低矿化度的储层仍有很好的应用效果[12-14]。

2 PNN测井计算水淹层剩余油饱和度

2.1 确定热中子宏观俘获截面

根据岩石体积物理模型，储层总的宏观俘获截面等于各部分的宏观俘获截面之和，对于水淹层来讲，地层的宏观俘获截面Σ表示为：

(1)

式中：Σma为岩石骨架的Σ；Σsh为泥质的Σ；Σh为剩余油气的Σ；Σw为混合地层水的Σ；Vsh为泥质含量；φ为孔隙度。则储层含水饱和度为：

(2)

式中：Sw为含水饱和度。为得到真实的地层宏观俘获截面Σ，消除流体性质、岩石骨架成分、储层孔隙度、仪器测量误差、统计误差等因素的影响，PNN仪器采用双源距(短源距425 mm，长源距745 mm)测量，36个时间道(每30 μs一道，共1 080 μs)计数，计算时取平均值的方法确定岩石宏观俘获截面[15]。

根据中子守恒定律，某一空间体积内中子密度随时间的变化率：

(3)

在脉冲中子发射的间隔，在源距适当的位置产生率与泄漏率为零，即在中子发射一段时间以后(快中子减速为热中子)，热中子密度变化率只与地层对热中子的俘获有关。

选择时刻t1和t2，此时热中子密度n1和n2分别为：

(4)

式中：τ为热中子寿命(μs)，热中子数量N1和N2与热中子密度n1和n2成正比，则有：

(5)

已知热中子寿命τ与宏观俘获截面的关系：

(6)

在常温条件下(25℃)热中子速度2 200 m/s，所以：

(7)

(8)

实际计算时，一般选取36个时间道的计数点，即时间间隔为1 080 μs，把数据分为6组，每组6个点的计数分别表示为N0、N1、N2、N3、N4、N5，按照式(7)计算每组内的热中子寿命。为减小涨落误差，第j组数据实际处理为：

(9)

第j组数据求得的地层宏观俘获截面为：

(10)

则地层实际宏观俘获截面为：

(11)

2.2 解释参数确定

2.2.1 参数取值范围及常规确定方法

式(1)和式(2)给出了准确计算含水饱和度及剩余油饱和度的算法，需在详细获取地层岩性参数及地层流体俘获特性的条件下才能进行准确计算。常用的参数确定方法如下。

1)Vsh、φ：由裸眼井测井解释得到，一般认为油层水淹后注水冲刷会导致泥质含量发生变化，所以Vsh也可以根据PNN带测的伽马曲线重新计算得到。而水淹过程对孔隙度的影响不大，可以用裸眼井的孔隙度解释结果。

2)Σma：主要由岩石矿物成分及其含量决定，不同岩性的Σma变化范围见表1。

3)Σsh：主要取决于黏土的矿物成分、分布形式，需要岩石物性分析资料，一般取解释井段内纯泥岩地层的Σ作为Σsh。

4)Σh：可根据油的密度、溶解油气比查相关的图版，也可以根据理论计算公式计算，如天然气的宏观俘获截面按下式计算：

(12)

油的宏观俘获截面为：

Σo=22.3/(1+GOR/22 000)0.715，

(13)

式中：γg为天然气对空气的相对密度；p为地层压力(lb/in2，1 lb/in2=0.332 5 Pa)；T为地层温度(℃)；GOR为油气比(ft3/bbl，1 ft3=0.028 3 m3，1 bbl=0.159 m3)。大多数油田的Σo接近21 c.u.。

5)Σw：主要与水的矿化度有关，根据等效的NaCl浓度查相应的图版，也可按公式计算：

(14)

式中：22是淡水的宏观俘获截面，C是地层水等效NaCl矿化度(mg/L)。地层水中盐离子主要俘获核素是氯，其他常见离子俘获能力比氯小得多，一般可忽略，但如果地层水中含硼和锂，因其俘获能力又比氯大得多，所以要计算等效NaCl矿化度。按式(14)计算，Σw的范围一般在22～120 c.u.。

表1 常见岩石骨架的宏观俘获截面Table 1 Macroscopic capture sections of common rock skeletons c.u.

2.2.2 遗传算法确定解释参数

从式(2)可知，利用宏观俘获截面计算Sw需要准确输入Σma、Σh、Σw、Σsh等参数，而由2.2.1节的分析可知，对于某个区块而言，以上参数是在一定范围内变化的，且受很多因素影响，这必然导致计算公式具有多解性。为此，提出基于遗传算法的PNN参数优化技术，即在标准层中选取曲线无异常，完井解释精度高的点，利用裸眼井解释的Vsh、Φ、Sw等参数，采用遗传算法求取Σma、Σh、Σsh。

1)遗传算法流程。遗传算法是用于解决最优化问题的一种搜索算法，其主要步骤包括设置初始群体，计算群体中每个个体的适应度，经过选择、交叉、变异运算之后得到下一代群体，然后进行终止条件判断，若满足条件则以进化过程中所得到的具有最大适应度个体作为最优解输出，终止计算，否则继续迭代。其流程图如图1[16-17]。

图1 遗传算法流程Fig.1 Genetic algorithm flow

2)选择、交叉、变异概率设定[18]。将某个体被选中的概率与其适应度函数建立如下关系式：

(15)

式中：F为该个体的适应度。

首先将P和F作归一化处理，处理后记为p和f，当F≥1时，p≥f。为避免算法过快收敛于局部最优解，可利用p代替P作为选择概率，以增加适应度较小的个体被选中的概率。

(16)

(17)

变异概率pm设定如下:

(18)

(19)

3)目标函数及样本点选择。遗传算法的计算过程即目标函数最大值(或最小值)的优化问题[19]，其数学通用模型表示为：

(20)

其中:f(X)为目标函数，X∈R、R⊂U为约束条件，U是基本空间，R是U的子集，满足约束条件的解X为可行解。本例中要求极小化的目标函数，可表示为：

(21)

2.3 逐层计算水淹层混合地层水宏观俘获截面

由2.2.1节可知，Σw主要与地层水的矿化度有关且受盐离子类型的影响，而长期的注水开发导致水淹层混合地层水矿化度变化大，同时受注采期次及采收率影响，即使相邻油层组其混合地层水矿化度也差距明显。实际计算结果也显示，根据标准层数据计算出的Σma、Σh、Σsh较为准确，而Σw的结果却无法与各个小层完全匹配，导致根据2.2.2节中的方法利用标准层的参数来计算相邻目的层的含水饱和度结果不准确。所以，分别计算每个小层的Σw，对利用PNN测井计算剩余油饱和度十分必要[20]。

根据式(1)可知，若将孔隙内的油气和水看成一个整体，其宏观俘获截面记为孔隙宏观俘获截面Σφ，则：

(22)

为消除孔隙度大小对宏观俘获截面的影响，设Σf=Σφ/φ为流体宏观俘获截面，只受流体内油、水相对含量的影响，即：

(23)

则剩余油饱和度可表示为：

(24)

在要确定Σw的目的层内，选取水淹程度较一致且原始含油饱和度不同的两点M、N，为保证满足上述条件，M、N点最好在一个厚层的中部且自然伽马、中子孔隙度、密度、声波时差曲线特征相同，即二者的岩性、物性、泥质含量等基本无差别，最大限度地保证在注水开发过程中驱油效率一致，即水淹程度相同，而二者的电阻率不同，即含油饱和度不同。

根据裸眼井解释和遗传算法迭代的结果，原始含水(含油)饱和度、孔隙度和泥质含量及Σma、Σh、Σsh、Σh等参数已知，Σw未知，则根据式(23)可分别计算M、N点的流体宏观俘获截面ΣfM和ΣfN。已知两点驱油效率一致，则：

SoM-SoN=SooM-SooN，

(25)

式中：SoM、SoN为M、N点的剩余油饱和度，未知；SooM、SooN为M、N点的原始含油饱和度，已知。

将式(24)代入(25)，得到：

(26)

可以直接根据式(26)得到本层的混合地层水宏观俘获截面Σw，M和N点的Σw相同。

3 应用效果分析

3.1 应用工区概况

3.2 PNN测井解释结果分析

根据研究区注采情况及产出水矿化度变化(图2)，可以推断同一解释井段各个小层的混合地层水宏观俘获截面应该不同，单个小层混合地层水宏观俘获截面的变化对PNN解释结果影响较大。所以，利用逐层计算混合地层水宏观俘获截面的方法，将每小层的Σw代入饱和度计算公式，最终计算得到的含油饱和度更接近真实地层的剩余油饱和度。

如图3所示，研究区某井于2007年11月投产，初期日产油量4.4 t，日产水3 m3，含水40.41%。2008年4月补孔转注，转注后由于反洗井不通换封。转注前累计产油566 t，累计产水529 m3。2009年4月换封。2011年5月因套管错断换封失败，下混注完井。目前该井套管错断停注。截止到2011年9月底累计注水5.234 4 万m3。2012年2月PNN测井在套管井条件下完成，测量井段内包含砂岩、泥岩以及泥质砂岩储层序列。

图2 研究区块产出水矿化度变化Fig.2 Changes in salinity of produced water in the study area

图3 PNN测井解释结果Fig.3 Interpretation results of PNN logging

裸眼井解释结果显示，Ⅴ2-19、Ⅴ3-25、Ⅴ3-26、Ⅴ3-27四个小层皆为油层，而根据PNN测井资料，Ⅴ2-19、Ⅴ3-25、Ⅴ3-26三小层地层宏观俘获截面值分别为16.1 c.u.、18.2 c.u.、22.4 c.u.，分别为低俘获值、中等俘获值、较高俘获值，PNN解释含油饱和度分别为66%、62.6%、45.38%，综合分析分别为油层、油层、油水层。水淹解释分别为未水淹、低水淹、高水淹。Ⅴ3-27小层地层宏观俘获截面值20.4 c.u.，为中高俘获值，但该层较Ⅴ3-26小层孔隙度高，岩石骨架对地层宏观俘获截面的贡献更小，因此PNN解释含油饱和度为39.35%。该层温度呈异常反应，该层下面围岩段地层宏观俘获截面值呈异常高值，因该井1 670 m处套管错段，疑为该段水泥环缺失，注入水从错段处进入上窜至该层，综合分析为含油水层，水淹解释为高水淹层。通过PNN测试检查各小层剩余油分布情况，根据定性分析和定量解释计算结果，并结合裸眼井资料、邻井动态资料综合分析，建议射开Ⅴ-19、Ⅴ-25、Ⅴ-26三小层进行生产，投产后产水率分别为23.47%、27.34%、45.96%，半年后稳定产量日产油4 t左右，含水较少。