作者简介:徐婕 工程师,1988 年生,2010 年毕业于长江大学资源勘查专业,现在中国石油渤海钻探第二录井公司从事工程地质一体化工作。通信地址:062552 河北省任丘市渤海钻探第二录井公司地质研究评价中心。电话:(0317)270427。E-mail:xu_jie@cnpc.com.cn
在同一岩性地层中,由于钻速、转盘转速、钻压、钻井液密度和岩石密度等录井参数存在不确定性,导致基于录井参数的随钻地层压力预测存在一定的不确定性。针对该问题,分析了钻速、转盘转速、钻压、钻井液密度和岩石密度等录井参数的不确定性特征,进而采用改进的Rosenbluthe方法推导得出 dc指数、上覆地层垂向应力及随钻地层压力的不确定性分析表达式,同时提出了基于录井参数的随钻地层压力不确定性分析方法,最后得到了地层压力的概率分布。以某区块HT 101井为例,分别应用不确定性分析方法和蒙特卡罗方法对地层压力的概率分布进行预测,结果表明,不确定性分析方法与蒙特卡罗方法的压力预测一致性达95%以上,其压力预测精度能够满足工程需要。
In the strata of the same lithology, due to the uncertainties of mud logging parameters such as drilling speed, rotary speed, weight on bit, drilling fluid density and rock density, some uncertainties exist in predicting formation pressure while drilling based on mud logging parameters.To solve this problem, the uncertainty characteristics of mud logging parameters including drilling speed, rotary speed, weight on bit, drilling fluid density and rock density were analyzed. Then, the analytical expressions of uncertainties in dc index, vertical stress of overlying strata and formation pressure while drilling were derived by using the improved Rosenbluthe method. At the same time, a method for analyzing the uncertainties of formation pressure while drilling based on mud logging parameters was proposed, and the probability distribution of formation pressure was finally obtained. Taking well HT 101 in a block as an example, uncertainty analysis method and Monte Carlo method were used to predict the probability distribution of formation pressure respectively. The results showed that the consistency of pressure prediction between uncertainty analysis method and Monte Carlo method was more than 95%, and the accuracy of pressure prediction can meet the engineering needs.
地层压力是指地层中流体所承受的压力[1, 2], 其对于维持钻井过程中的井壁稳定性具有重要作用[3, 4]。地层压力预测尤其是随钻地层压力预测一直是钻井工程的核心问题, 目前随钻地层压力预测中大多应用录井数据[5], 然而即使对于同一地层, 钻速、转盘转速、钻压、钻井液密度、岩石密度等录井数据均存在一定的不确定性, 这导致基于上述参数计算得出的dc指数和上覆地层垂向应力存在不确定性, 进而使得随钻地层压力预测存在不确定性, 其预测结果也不是确切数值, 而是一个概率分布区间。因此, 通过对随钻地层压力进行不确定性分析, 获取地层压力的概率分布区间, 对实际钻井工程施工具有重要指导作用。
目前, 通过录井参数和dc指数建立随钻地层压力确定性模型已取得进展:Hottman和Johnson通过建立确定性模型预测随钻地层压力[6]; Penebaker利用地震资料确定性预测随钻地层压力, 并确认准确获取层速度是用地震资料预测地层压力的关键; Eaton总结前人研究成果, 提出了基于电阻率、声波时差和dc指数的随钻地层压力关系式, 即伊顿法[7]; Fillippone提出一个不依赖于正常压实趋势线, 而是利用层速度来预测地层压力的经验模型[8]; 刘震等对Fillippone确定性模型进行了修正, 提高了地层压力预测精度[9]。
此外, 部分学者引入不确定性分析方法预测地层压力:胜亚楠等[10, 11]结合蒙特卡罗法模拟和概率分布理论, 建立了含可信度的地层压力求取方法, 得到了地层压力和钻井液密度的概率分布; 李莹莹[12]通过处理已钻地层地震或测井资料, 利用概率统计理论对地层压力进行了不确定性分析, 定量描述含可信度区间的地层压力纵向剖面; 于浩[13]使用基于机器学习的多变量地层压力预测方法获取多组预测地层压力并使用多组模拟对预测压力不确定性进行量化分析; 黄熠等[14]分析了地层压力的不确定性, 并基于不确定度分析理论对钻井液密度的不确定性进行了概率分布分析; 魏凯等[15]通过录井参数不确定性的概率分布和不确定度两种描述方法, 建立地质参数不确定性分析方法, 确定钻井地质参数的概率分布函数。
虽然前人在地层压力预测方面已取得很多研究成果, 但仍然存在以下不足:一是由于地层压力的不确定性, 预测的结果应该是一个区间, 而非一个固定的数值, 所以确定性模型很难对地层压力进行精确评价; 二是虽然已开展了地层压力的不确定性分析, 但是几乎所有的不确定性分析都是基于蒙特卡罗法, 而蒙特卡罗法计算效率较低, 制约了该方法在实际工程中的应用。针对上述问题, 本文通过引入改进的Rosenbluthe不确定性方法建立了基于录井参数的随钻地层压力预测模型, 获取地层压力的概率分布, 并与蒙特卡罗法对比验证其准确性。
在同一岩性地层中, 钻速、转盘转速、钻压、钻井液密度、岩石密度等录井参数都存在一定的不确定性, 这导致同一地层各深度处的随钻地层压力预测存在不确定性。本文以某区块HT 101井5 010~6 357 m井段的录井资料为例, 分别对其钻速、转盘转速、钻压、钻井液密度及岩石密度等录井参数的不确定性数学特征进行统计。
由图1a可知, HT 101井5 010~6 357 m井段的钻速在0~0.4 m/min范围内, 通过统计每0.04 m/min区间内的钻速数据个数, 拟合出钻速概率分布如图1b所示, 钻速具有不确定性, 其分布形式基本上满足正态分布。
由图2a可知, HT 101井5 010~6 357 m井段的转盘转速在50~100 r/min范围内, 通过统计每10 r/min区间内的转盘转速数据个数, 拟合出转盘转速概率分
布如图2b所示, 转盘转速具有不确定性, 其分布形式基本上满足正态分布。
同理, 采用上述分析方法对HT 101井5 010~6 357 m井段钻压、钻井液密度、岩石密度等其他录井参数进行了不确定性分析, 结果表明上述录井参数同钻速、转盘转速一样, 均具有不确定性, 其分布形式基本上均满足正态分布。钻速、转盘转速、钻压、钻井液密度及岩石密度等录井参数作为随机变量, 其不确定性特征如表1所示。
![]() | 表1 HT 101井5 010~6 357 m井段录井参数不确定性特征 |
Rosenbluthe于1985年提出了一种计算由不同随机变量X构成的函数Y的简便方法, 称为Rosenbluthe法[16]。崔维成等[17]针对Rosenbluthe法存在的计算数值奇异和计算精度不稳等问题进行了改进, 提出了改进的Rosenbluthe法, 其计算过程如下。
假设随机变量为X1, X2, …, Xn, 其均值和标准差分别为μ X1, μ X2, …, μ Xn和σ X1, σ X2, …, σ Xn, 由此计算各随机变量的函数值y* 、yi+、yi-, 公式如下。
y* =G(μ X1, …, μ Xi, …, μ Xn) (1)
yi+=G(μ X1, …, μ Xi+σ Xi, …, μ Xn) (2)
yi-=G(μ X1, …, μ Xi-σ Xi, …, μ Xn) (3)
根据公式(2)和公式(3)可计算两个系数μ i和σ i, 公式如下。
通过以下公式可确定
yij+=G(μ X1, …, μ Xi+σ Xi, …, μ Xj+σ Xj, …, μ Xn) (6)
yij-=G(μ X1, …, μ Xi-σ Xi, …, μ Xj-σ Xj, …, μ Xn) (7)
通过以下公式计算μ ij、Δ μ i、Δ μ ij的函数值:
Δ μ i=μ i-y* (9)
Δ μ ij=μ ij-y* -Δ μ i-Δ μ j (10)
最后, 将上述函数值代入公式(11), 得到由不同随机变量X构成的函数Y的均值和标准差:
由于钻速、转盘转速、钻压等录井工程参数, 以及钻井液密度、岩石密度参数均为具有不确定性的随机变量, 导致基于录井参数得出的dc指数和上覆地层垂向应力也属于具有不确定性的随机变量。
本文针对录井参数的不确定性, 采用改进的Rosenbluthe法推导得出了dc指数[16]这一随机变量的不确定性分析表达式, 其均值和标准差公式如下:
其中:
式中:μ dc为dc指数的均值, 无量纲; Db为钻头直径, mm; Gn为静水压力梯度, MPa/m; μ X1为钻速的均值, m/min; μ X2为转盘转速的均值, r/min; μ X3为钻压的均值, kN; μ X4为钻井液密度的均值, g/cm3; σ dc为dc指数的标准差, 无量纲; yi+、yi-、yij+、yij-(i, j=1、2、3、4, 且i< j)为函数值, 可由μ X1、σ X1、μ X2、σ X2、μ X3、σ X3、μ X4、σ X4代入公式(2)、(3)、(6)、(7)获得(σ X1为钻速的标准差, m/min; σ X2为转盘转速的标准差, r/min; σ X3为钻压的标准差, kN; σ X4为钻井液密度的标准差, g/cm3)。
采用改进的Rosenbluthe法推导得出了上覆地层垂向应力[17]这一随机变量的不确定性分析表达式, 其均值和标准差见公式(13):
其中:
式中:μ σ V为上覆地层垂向应力的均值, MPa; H为深度, m; σ σ V为上覆地层垂向应力的标准差, MPa; y5+、y5-为函数值, 由μ X5、σ X5代入公式(2)、(3)获得(μ X5为岩石密度的均值, g/cm3; σ X5为岩石密度的标准差, g/cm3)。
针对dc指数和上覆地层垂向应力的不确定性, 基于Eaton法[18]和改进的Rosenbluthe法推导得出了随钻地层压力的不确定性分析表达式, 其均值和标准差见公式(14):
其中:
式中:
计算求得随钻地层压力的均值和标准差后, 再结合地层压力作为随机变量的最合理分布形式, 即可得到随钻地层压力的概率分布。
以HT 101井5 010~6 357 m井段为例, 进行基于录井参数不确定性的随钻地层压力测算应用分析, 并采用行业内公认的测算精度较高的蒙特卡罗法对dc指数和随钻地层压力进行不确定性分析所获得的预测结果来验证本文研究方法的精准度。
根据公式(1)-公式(13)和表1中的数据, 可以得到HT 101井5 010~6 357 m井段dc指数和上覆地层垂向应力的均值和标准差(表2)。由于表1中的随机变量满足正态分布, 由表1中的随机变量构成的dc指数和上覆地层垂向应力函数也满足正态分布。图3为HT 101井5 010~6 357 m井段dc指数和上覆地层垂向应力的概率分布。
![]() | 表2 dc指数和上覆地层垂向应力的均值和标准差 |
根据公式(14)和表2中的数据, 可以计算得到HT 101井5 010~6 357 m井段随钻地层压力的均值为69.78 MPa, 标准差为2.23 MPa。由于随钻地层压力公式中的随机变量dc指数和上覆地层垂向应力满足正态分布, 随钻地层压力也应满足正态分布, 其概率分布如图4所示。
为了验证本文所研究方法的精度, 利用蒙特卡罗法[19, 20]对HT 101井5 010~6 357 m井段dc指数和随钻地层压力进行了不确定性分析, 模拟次数为20万次, 得到了其dc指数和随钻地层压力的蒙特卡罗法模拟结果(图5)。将本文所研究的不确定性分析方法得到的dc指数和随钻地层压力概率分布曲线与蒙特卡罗法模拟得到的结果进行比较可知(图6), 两种方法的概率分布曲线具有很高的一致性, 吻合程度超过95%。同时, 对两种方法得到的dc指数和随钻地层压力的最大值、最小值、均值、标准差等特征属性进行对比(图7)可知, dc指数和随钻地层压力的最大值、最小值、均值、标准差之间的差异非常小, 吻合程度高达98%。综合以上分析结果, 从概率分布曲线和随机变量特征属性的角度来看, 本文研究方法所得到的结果与蒙特卡罗法得到的结果非常接近, 可以满足实际工程需要。
基于HT 101井的密度测井和声波测井曲线及表1数据, 并结合前文所提方法, 计算得出HT 101井5 010~6 357 m井段地层压力分布区间。图8为该井段地层压力分布区间和实测值对比。计算结果表明, 所有地层压力实测值均位于地层压力分布区间的中值线附近, 这在一定程度上验证了本文研究方法所测算得出的地层压力分布区间的合理性。以HT 101井5 704 m处的地层压力为例, 该深度处的地层压力分布区间为57.87~67.87 MPa, 分布区间中值为62.87 MPa, 地层压力实测值为59.15 MPa, 分布区间中值与实测值之间的误差仅为6.3%。另外, 其他深度处的地层压力分布区间中值和实测值误差也均低于10%, 表明本文所研究方法的计算精度可以满足现场要求。
(1)针对钻速、转盘转速、钻压、钻井液密度、岩石密度不确定性以及由此引起的dc指数和上覆地层垂向应力不确定性对随钻地层压力的影响, 提出了基于改进的Rosenbluthe方法和Eaton法的随钻地层压力不确定性分析方法。
(2)建立了随钻地层压力不确定性分析的完整计算流程。首先统计钻速、转盘转速、钻压以及钻井液密度、岩石密度的不确定性特征; 然后由推导得出的dc指数和上覆地层垂向应力的不确定性表达式, 求出其均值和标准差; 最后将dc指数和上覆地层垂向应力的均值和标准差代入随钻地层压力的不确定性表达式, 得到地层压力的概率分布。
(3)将本文所研究的不确定性分析方法得到的dc指数和随钻地层压力的概率分布与蒙特卡罗法模拟结果进行了比较, 两者的吻合程度超过95%; 将不确定性分析方法所得到的地层压力与密度测井和声波测井曲线得到的结果进行比较, 误差均低于10%, 验证了本文所提方法的准确性。
(编辑 唐艳军)
[1] |
|
[2] |
|
[3] |
|
[4] |
|
[5] |
|
[6] |
|
[7] |
|
[8] |
|
[9] |
|
[10] |
|
[11] |
|
[12] |
|
[13] |
|
[14] |
|
[15] |
|
[16] |
|
[17] |
|
[18] |
|
[19] |
|
[20] |
|