Key Drilling Fluid Technology in the Ultra Deep Section of Well Ying-1 in the Shunbei Oil and Gas Field
-
摘要:
鹰1井是顺北油气田的一口超深重点风险预探井,设计井深9 016.85 m(垂深8 603.00 m)。该井超深井段志留系柯坪塔格组与奥陶系桑塔木组等硬脆性泥岩地层、志留系裂缝性地层和奥陶系破碎性地层,在钻进过程中易出现井眼失稳、井漏、坍塌掉块等井下故障。为此,通过室内试验研究,分析了该井超深井段硬脆性泥岩地层井眼失稳机理、强压力敏感性裂缝性地层漏失原因及破碎性碳酸盐岩地层井眼失稳原因,应用“多元协同”井壁稳定基本理论,构建了SMHP–1强抑制强封堵钻井液,并制定了针对性强的防塌防漏技术措施。该井顺利钻穿大段硬脆性泥岩、裂缝性地层和破碎性地层,未发生井眼失稳及钻井液漏失,顺利钻至井深8 588.00 m完钻,创亚洲陆上井深最深纪录。现场应用表明,SMHP–1强抑制强封堵钻井液能够解决深部地层大段泥岩及破碎性地层的井眼失稳与漏失难题,为国内外深井超深井安全钻进提供了技术借鉴。
Abstract:Well Ying-1 is an ultra-deep key investigative risk management pre-exploration well deployed by Sinopec in the Shunbei Oil and Gas Field with designed well depth of 9 016.85 m (TVD 8 603.00 m). Drilling the welln encountered downhole problems including wellbore instability, well leakage and borehole wall sloughing. They also easily occur in drilling hard brittle mudstone formations, such as the Silurian Kepingtage Formation and the Ordovician Sangtamu Formation. For this reason, laboratory studies have been carried out to analyze mechanisms contributing to the instability of large section of hard brittle mudstone shale. As such, the reasons for leakage in the Silurian high pressure sensitive fissured formation, and the reasons for instability in Ordovician fractured formations. Based on the basic theory of " multivariate synergistic” wellbore stability, the drilling fluid system SMHP-1 with strong inhibition and sealing capacity was constructed, and the technical measures of mud loss and borehole wall collapse prevention were worked out. The well successfully drilled through a large section of hard brittle mudstone and broken formation to a total depth of 8 588 m without borehole wall instability and drilling fluid loss, setting a record of the deepest onshore well depth in Asia. The field application showed that the drilling fluid system SMHP-1 could effectively solve the problems of wellbore stability and leakage in deep mudstone and broken formations by strong inhibition and sealing capacity, and provide best practices for safe drilling of deep or ultra-deep wells at home and abroad.
-
现代地质导向钻井过程中,需采用近钻头随钻边界探测技术动态测量近钻头井眼轨迹参数(井斜角、方位角和工具面角),根据随钻成像预估地层和井眼变化趋势,据此调整钻头方向,提高油层钻遇率。由于振动、旋转、磁干扰等因素的影响,导致测量的近钻头井眼轨迹参数出现很大误差,严重影响了随钻测量精度,出现钻穿油层等情况。
目前旋转导向动态测量技术都掌握在国外油田技术服务公司手中,不对外公开,因此,未见到国外有井眼轨迹参数动态测量技术的报道。国内杨全进等人[1]建立了一种旋转导向系统有色噪声的改进无迹卡尔曼滤波方法,该算法限制条件较多。高怡等人[2]提出了采用多源动态姿态组合测量方法测量导向钻具的动态姿态,但加速度计采集数据的存储周期较长,造成测量误差较大。Xue Qilong等人[3-4]提出了基于卡尔曼滤波状态空间模型的动态井眼轨迹测量方法,但其只适用于线性系统。徐宝昌等人[5]提出了基于无迹卡尔曼滤波的动态姿态测量方法,但没有解决振动及旋转对动态姿态测量的影响。为此,笔者提出了一种基于数据融合的近钻头井眼轨迹参数动态测量方法,该方法针对三轴加速度计、磁通门和速率陀螺的测量系统,建立了基于四元数井眼轨迹参数测量模型,推导出加速度计与井眼轨迹平滑预测的关系,运用3个捷联式无迹卡尔曼滤波器和磁干扰校正系对加速度计、磁通门进行消噪,校正,实时测量近钻头井眼轨迹参数,提高随钻近钻头探边能力,确定最佳储层位置。
1. 近钻头动态井眼轨迹测量模型
近钻头动态测量系统由三轴加速度计、三轴磁通门和角速率陀螺仪组成,针对该系统,基于四元数方法建立井眼轨迹非线性数学模型。如图1所示,根据地理坐标系O-NED和钻具坐标系O-xyz的对应关系,建立欧拉角转换矩阵,并转换为四元数,k时刻姿态转换矩阵T表示为:
{\boldsymbol{T}}\left( k \right){{ = }} \left[ \begin{array}{l} {q_{\text{0}}^{\text{2}} + q_{\text{1}}^{\text{2}} - q_{\text{2}}^{\text{2}} - q_{\text{3}}^{\text{2}}}\\ {{\text{2}}\left( {{q_{\text{1}}}{q_{\text{2}}} + {q_{\text{0}}}{q_{\text{3}}}} \right)}\\ {{\text{2}}\left( {{q_{\text{1}}}{q_{\text{3}}} - {q_{\text{0}}}{q_{\text{2}}}} \right)}\end{array} \begin{array}{l} {{\text{2}}\left( {{q_{\text{1}}}{q_{\text{2}}} - {q_{\text{0}}}{q_{\text{3}}}} \right)}\\ {q_{\text{0}}^{\text{2}} - q_{\text{1}}^{\text{2}} + q_{\text{2}}^{\text{2}} - q_{\text{3}}^{\text{2}}}\\ {{\text{2}}\left( {{q_{\text{2}}}{q_{\text{3}}} + {q_{\text{1}}}{q_{\text{0}}}} \right)}\end{array} \begin{array}{l} {{\text{2}}\left( {{q_{\text{1}}}{q_{\text{3}}} + {q_{\text{0}}}{q_{\text{2}}}} \right)} \\ {{\text{2}}\left( {{q_{\text{2}}}{q_{\text{3}}} - {q_{\text{0}}}{q_{\text{1}}}} \right)} \\ {q_{\text{0}}^{\text{2}} - q_{\text{1}}^{\text{2}} - q_{\text{2}}^{\text{2}} + q_{\text{3}}^{\text{2}}} \end{array} \right] (1) 式中:q0,q1,q2和 q3为四元数。
理论上由于x轴和y轴的加速度计、磁通门正交相隔90°,其信号为正弦或余弦曲线。由于钻头旋转、振动的干扰,近钻头x轴和y轴加速度计、磁通门的实测信号如图2所示。由图2可以看出,近钻头的旋转、振动对磁通门的影响相对于加速度计较小,所以采取磁通门读数校正加速度计读数。
\frac{{{B_x}(k + 1)}}{{{B_x}(k)}} = \frac{{{a_x}(k + 1)}}{{{a_x}(k)}} = \frac{{B\sin ({w_x}(k + 1))}}{{B\sin ({w_x}(k))}} = \frac{{g\sin ({w_x}(k + 1))}}{{g\sin ({w_x}(k + 1))}} (2) \frac{{B\sin ({w_y}(k + 1))}}{{B\sin ({w_y}(k))}} = \frac{{g\sin ({w_y}(k + 1))}}{{g\sin ({w_y}(k + 1))}} = \frac{{{B_y}(k + 1)}}{{{B_y}(k)}} = \frac{{{a_y}(k + 1)}}{{{a_y}(k)}} (3) \left[ \begin{gathered} {a_x}(k + 1) \hfill \\ {a_y}(k + 1) \hfill \end{gathered}\right] = \left[ \begin{gathered} \frac{{{B_x}(k + 1)}}{{{B_x}(k)}} \hfill \\ 0 \hfill \\ \end{gathered} \begin{gathered} 0 \hfill \\ \frac{{{B_y}(k + 1)}}{{{B_y}(k)}} \hfill \\ \end{gathered}\right] \left[ \begin{gathered} {a_x}(k) \hfill \\ {a_y}(k) \hfill \\ \end{gathered}\right] (4) 式中:g为重力加速度,m/s2;B为地磁场强度,μT;Bx
\left(k\right) ,By\left(k\right) 和Bz\left(k\right) 为k时刻三轴磁通门测量的磁场强度,μT;ax\left(k\right) ,ay\left(k\right) 和az\left(k\right) 为k时刻三轴加速度计测量结果,m/s2;wx\left(k\right) ,wy\left(k\right) 和wz\left(k\right) 为k时刻陀螺仪的角速度,rad/s。2. 数据融合近钻头井眼轨迹参数动态测量方法
基于数据融合算法的近钻头井眼轨迹参数动态测量方法的测量流程如图3所示,图中KF1和KF2为基于扩展卡尔曼滤波算法的滤波器,KF3为基于无迹卡尔曼滤波算法的滤波器。测量步骤:1)将加速度计、磁通门、转动角速度四元数带入KF1滤波器,进行扩展卡尔曼滤波,得出井斜角、方位角估计值;2)将加速度计四元数带入KF2滤波器,进行扩展卡尔曼滤波,得出测深增量Δhm;3)将测深增量Δhm、井斜角、方位角估计值带入KF3滤波器,进行无迹卡尔曼滤波,得出井斜角、方位角最终估计值;4)利用井斜角、方位角最终估计值计算磁性工具面角ωm与重力工具面角的差Δω;5)利用磁性工具面角和角差Δω求出重力工具面角ωg。
近钻头振动信号是一种幅值大、频率高、频带宽的噪声信号,可以近似等效为高斯白噪声。动态测量的动力学模型为典型的非线性模型。根据z轴陀螺仪测量的转速wz,运用转速补偿策略,对x轴和y轴加速度计测量结果进行转速补偿[5],减小钻具旋转对加速度计测量结果的影响。根据转速补偿结果以及噪声特性,采用扩展、无迹卡尔曼滤波滤除振动干扰信号。
2.1 估计近钻头井斜角、方位角的扩展卡尔曼滤波算法
基于四元数的KF1的状态方程和量测方程:
Q(k + 1) = ({\boldsymbol{I}}{\text{ + }}{t_{\text{s}}}{\boldsymbol{A}}(k))Q(k) + w(k) (5) Z(k + 1) = F(Q(k)) + v(k) (6) 式中:Q(k)为k时刻的状态值;I为单位矩阵;ts为采样周期;w(k)为k时刻系统高斯白噪声;v(k)为k时刻传感器观测噪声;A(k)为k时刻状态转移矩阵;F(x)为非线性函数;Z(k+1)为k+1时刻的观测值。
Z(k + 1) = \left[ \begin{gathered} {B_x} \hfill \\ {B_y} \hfill \\ {B_z} \hfill \\ {a_x} \hfill \\ {a_y} \hfill \\ {a_z} \hfill \\ \end{gathered} \right] = \left[ \begin{array}{l} T(k) \left[ \begin{array}{c} B\cos \theta \hfill \\ 0 \hfill \\ B\sin \theta \hfill \end{array} \right]\\ T(k) \left[\begin{gathered} 0 \hfill \\ 0 \hfill \\ g \hfill \\ \end{gathered} \right] \end{array}\right] + v(k) (7) \begin{split} & Q(k + 1) = \left(I + {t_{\text{s}}} \left[ \begin{array}{c} 0 \hfill \\ {w_x}(k) \hfill \\ {w_y}(k) \hfill \\ {w_z}(k) \hfill \\ \end{array} \begin{array}{c} - {w_x}(k) \hfill \\ 0 \hfill \\ - {w_z}(k) \hfill \\ {w_y}(k) \hfill \\ \end{array} \begin{array}{c} - {w_y}(k) \hfill \\ {w_z}(k) \hfill \\ 0 \hfill \\ - {w_x}(k) \hfill \\ \end{array} \begin{array}{c} - {w_z}(k) \hfill \\ - {w_y}(k) \hfill \\ {w_x}(k) \hfill \\ 0 \hfill \\ \end{array} \right] \right) \cdot \\ & \qquad\qquad Q(k) + w(k)\end{split} (8) 三轴加速度信号、三轴磁通门信号、角速率陀螺信号进行数据融合后,采用扩展卡尔曼滤波算法,得到最优姿态估计[6-7],动态解算出钻井工具的实时姿态参数,确保钻具姿态测量计算的精度,减少计算量,对四元数Q进行更新,求出KF1滤波后的井斜角αKF1、方位角ϕKF1、高边工具面角ωg,KF1和磁性工具面角ωm,KF1。
{\alpha _{{\text{KF}}1}}{\text{ = }}\arctan \frac{{2({q_0}{q_1} + {q_2}{q_3})}}{{1 - 2(q_1^2 + q_2^2)}} (9) {\phi _{{\text{KF}}1}}{\text{ = }}\arctan \frac{{2({q_0}{q_3} + {q_1}{q_2})}}{{1 - 2(q_0^2 + q_3^2)}} (10) {\omega _{{\text{g}},{\text{KF}}1}}{\text{ = }}\arctan \frac{{({q_0}{q_2} + {q_1}{q_3})}}{{({q_0}{q_1} - {q_2}{q_3})}} (11) {\omega _{{\text{m}},{\text{KF}}1}}{{ = }}\arctan \frac{{({q_1}{q_2} + {q_0}{q_3})\cos \theta + ({q_1}{q_2} + {q_0}{q_3})\sin \theta }}{{(q_0^2 - q_1^2 - q_2^2 + q_3^2)\cos \theta + ({q_1}{q_3} - {q_0}{q_2})\sin \theta }} (12) 式中:αKF1,ϕKF1,ωg,KF1和ωm,KF1分别为KF1滤波后的井斜角、方位角、重力工具面角和磁性工具面角,(°);θ为地层倾角,(°)。
2.2 估计近钻头测深增量的扩展卡尔曼滤波算法
根据式(7)计算出四元数中的az,运用扩展卡尔曼滤波器计算系统经过ts后测深增量Δhm。z轴加速度计主要受到重力加速度和振动的干扰,由于采样时间ts为毫秒级,在单位采样周期内,重力加速度和振动的干扰可以视为近似相同[8-10],可以忽略振动对加速度计测量结果的影响。k为当前采样点,z轴加速度增量Δaz:
\Delta {a_z} = {a_z}(k + 1) - g\cos ({\alpha _{{\text{KF1}}}}(k)) (13) \Delta {a_z} = \Delta h_{\text{m}}^{''} (14) 为了提高对测深增量的估计,对Δhm进行二阶泰勒展开:
\Delta {h_{\text{m}}}(k + 1) = \Delta {h_{\text{m}}}(k) + \Delta {h_{\text{m}}}{(k)^{'}}{t_{\text{s}}} + 0.5\Delta {h_{\text{m}}}{(k)^{''}}t_{\text{s}}^{\text{2}} (15) KF2的状态方程和量测方程为:
\left[ \begin{gathered} \Delta {h_{\text{m}}}(k + 1) \\ \Delta {h_{\text{m}}}{(k + 1)^{'}} \\ \Delta {h_{\text{m}}}{(k + 1)^{''}} \\ \end{gathered} \right] = \left[ \begin{gathered} 1 \;\; {t_{\rm{s}}} \;\; t_{\rm{s}}^2 \hfill \\ 0 \;\;1 \;\;0 \hfill \\ 0\;\; 0 \;\;1 \hfill \end{gathered}\right] \left[ \begin{gathered} \Delta {h_{\text{m}}}(k) \\ \Delta {h_{\text{m}}}{(k)^{'}} \\ \Delta {h_{\text{m}}}{(k)^{''}} \\ \end{gathered} \right] + w(k) (16) \Delta {a_z} = \left[0 \; 0 \; 1\right] \left[ \begin{gathered} \Delta {h_{\text{m}}}(k + 1) \hfill \\ \Delta {h_{\text{m}}}{(k + 1)^{'}} \hfill \\ \Delta {h_{\text{m}}}{(k + 1)^{''}} \hfill \\ \end{gathered} \right] + v(k) (17) 2.3 估计近钻头井眼轨迹参数的无迹卡尔曼滤波算法
如图4所示,在单位采样时间内,井眼轨迹趋于平滑曲线,可以根据前面2个测点的狗腿度和KF2输出测深增量对井眼轨迹进行递归式预测[11]。
{h_{\text{m}}}(k + 1) = {h_{\text{m}}}(k) + \Delta {h_{\text{m}}}(k + 1) (18) \begin{split} & \;\;\;\;\;\;\; \gamma = \arccos \left[ {\cos \alpha (k)\cos \alpha (k - 1) + }\right. \\ & \left.{ \sin \alpha (k)\sin \alpha (k - 1)\cos (\phi (k) - \phi (k - 1))} \right] \end{split} (19) \Delta \gamma (k + 1) = \frac{{\gamma \Delta {h_{\text{m}}}(k + 1)}}{{\Delta {h_{\text{m}}}(k)}} (20) 式中:γ为狗腿角,(°);Δγ为系统经过ts的狗腿角增量,(°);
{h_{\text{m}}}(k) 为k时刻的测深,m。根据
\Delta {h_{\text{m}}}(k) 进行无迹卡尔曼滤波,KF3的状态方程和量测方程为:\begin{array}{l} \quad \alpha (k + 1) = \arccos [\cos (\gamma + \Delta \gamma (k + 1))\cos \alpha (k - 1) +\\ \dfrac{{\sin (\gamma + \Delta \gamma (k + 1))}}{{\sin \gamma }}(\cos \alpha (k - 1)\cos \gamma - \cos \alpha (k))] + {w_{\alpha (k)}} \end{array} (21) \begin{array}{l} \qquad\; \phi (k+1)=\phi (k)+\mathrm{sgn}(\phi (k)-\phi (k-1) \cdot \\ {\rm{arccos}}\dfrac{\mathrm{cos}(\Delta \gamma (k+1)-\mathrm{cos}\alpha (k)\mathrm{cos}\alpha (k{-1})}{\mathrm{sin}\alpha (k)\mathrm{sin}\alpha (k{-1})}+{w}_{\phi (k)} \end{array} (22) 式中:wα和wϕ分别为井斜角和方位角的系统高斯白噪声。
{\alpha _{{\text{KF3}}}} = \alpha {\text{(}}k + 1) + \nu {}_{\alpha (k)} (23) {\phi _{{\text{KF3}}}} = \phi {\text{(}}k + 1) + {\nu _{\phi (k)}} (24) 式中:αKF3和ϕKF3分别为KF3滤波后的井斜角和方位角,(°);vα和vϕ分别为井斜角和方位角的系统观测噪声。
2.4 近钻头重力工具面角的估计
根据旋转测量原理,同一时刻的重力工具面角与磁工具面角的差与测量时刻的井斜角、方位角、地磁倾角呈现一定函数关系[12-13]。根据KF3求出的井眼井斜角和方位角计算磁性工具面角与重力工具面角的差Δω:
\Delta \omega = - 90 + \arctan \frac{{\sin {\phi _{{\text{KF}}3}}}}{{\cos {\alpha _{{\text{KF}}3}}\cos {\phi _{{\text{KF}}3}} - \tan \theta \sin {\alpha _{{\text{KF}}3}}}} (25) 根据Δω,计算旋近钻头动态重力工具面角估计值ωdg,e:
{\omega _{{\text{dg,e}}}} = {\omega _{{\text{m,KF3}}}} + \Delta \omega (26) 式中:ωdg,e为旋近钻头动态重力工具面角估计值,(°);ωm,KF3为KF3滤波后的磁性工具面角,(°)。
2.5 磁干扰情况下的磁性工具面角
近钻头重力工具面角需根据x轴和z轴磁通门传感器的读数求出[14-15],而磁通门测量结果不可避免地会受到周围电磁场的影响,为了降低井下钻具周围电磁场的影响程度,采用径向磁干扰的方法校正x轴和z轴磁通门的测量结果[16]。
磁场的干扰导致磁通门测量的磁场强度发生偏移和变形。磁干扰下的测量结果如图5所示。
在实际钻井过程中,井下仪器旋转一圈时,钻深可以忽略不计,可以看作仪器在原地旋转了一圈。z轴磁通门的测量结果可以认为没有发生变化,而x轴和y轴磁通门的测量值不断发生变化,如图2(a)所示。三轴磁通门传感器的测量数据记为(Bx,By,Bx),地球磁场可以看成一个固定值,即:
B_x^2 + B_y^2 + B_z^2 = {C^2} (27) 式中 :C为常数。
根据椭圆校正原理,对短时间内采集的Bx,By进行磁干扰校正,得出排除磁干扰的Bxm和Bym:
{a_{x{\text{f}}}} = \max \left(1,\frac{{{B_{y\max }} - {B_{y\min }}}}{{{B_{x\max }} - {B_{x\min }}}}\right) (28) b{}_{x{\text{f}}} = \frac{{{B_{x\max }} - {B_{x\min }}}}{2} - {a_{x{\text{f}}}}{B_{x\max }} (29) B{}_{x{\text{m}}} = {B_x}{a_{x{\text{f}}}} + {b_{x{\text{f}}}} (30) {a_{y{\text{f}}}} = \max \left(1,\frac{{{B_{x\max }} - {B_{x\min }}}}{{{B_{y\max }} - {B_{y\min }}}}\right) (31) b{}_{y\text{f}}={a}_{y\text{f}} \left(\frac{{B}_{y\mathrm{max}}-{B}_{y\mathrm{min}}}{2}-{B}_{y\mathrm{max}}\right) (32) B{}_{y{\text{m}}} = {B_y}{a_{y{\text{f}}}} + {b_{y{\text{f}}}} (33) 式中:axf和ayf为x轴磁通门测量磁场强度的校正系数;bxf和byf为 y轴磁通门测量磁场强度的校正值,μT;Bxm和Bym为x轴和y轴排除磁干扰后的磁场强度,μT。
3. 试验与分析
在实验室对NWD施加20~50 r/min转速、 1g~3g振幅的高斯白噪声振动信号,模拟井下钻进过程。采用上述测量方法测量井斜角和方位角,并与实际井斜角和方位角进行对比,结果见图6。由图6可以看出:测量井斜角与实际井斜角的最大误差为1.30°,最小误差为0.02°,平均误差为0.12°,方差为1.26°;测量方位角与实际方位角的最大误差为1.95°,最小误差为0.15°,平均误差为0.57°,方差为2.46°。
图7为采用上文测量方法(以下称为动态测量方法)测量工具面角、NWD探管直接测量工具面角与实际工具面角的差值。由图7可以看出,NWD探管直接测得工具面角与实际工具面角的差最大接近40°,已经不能满足测量近钻头伽马重力工具面角的要求,而采用动态测量方法测得的工具面角与实际工具面角的差小于10°,满足测量近钻头工具面角的要求。
在井斜角为0.5°、方位角度为142.818°,地磁场强度为53.74 μT、磁倾角为55.8°、磁性工具面角30°的测点,对NWD探管x、y方向分别施加0.1和0.2 μT的磁干扰,未进行磁干扰校正前测得磁性工具面角为23.56°,经过磁干扰校正后,磁性工具面角为29.82°,与真实工具面角相差0.18°。
为了验证动态测量方法的可行性和有效性,在某井3 800.00~4 100.00 m井段进行验证,首先在钻进过程利用动态测量方法测量井斜角、方位角和工具面角,然后再停钻测量井斜角、方位角和工具面角(下文称之为静态测量),然后对比静态和动态测量结果。因为钻具停止时没有旋转、振动和磁干扰,因此以静态测量结果为真实值。
图8为动态测得井斜角与静态测得井斜角的对比。由图8可以看出,动态测得井斜角与静态测得井斜角的平均误差为0.19°,方差为0.9342°。图9为动态测得方位角与静态测得方位角的对比。由图9可以看出,动态测得方位角与静态测得方位角的平均误差为1.19°,方差为1.943 642°。
图10为动态测量工具面角与静态测量工具面角的对比。由图10可以看出,动态测得工具面角与静态测得工具面角的误差小于5°。
以上分析可以得出,基于数据融合的井眼轨迹参数测量方法能够有效消除旋转、振动、磁干扰的影响,井斜角、方位角和工具面角的测量精度得到明显提高,满足了随钻地质导向对井眼轨迹参数测量的要求。
4. 结 论
1)针对旋转、振动、磁干扰对近钻头测量仪器的影响,提出了基于数据融合的井眼轨迹参数动态测量方法。该方法采用捷联式卡尔曼滤波器和磁干扰校正系统对测量信号进行滤波、校正,利用四元数法处理测量数据,求出近钻头井眼轨迹参数。
2)实验室模拟试验和现场试验均表明,基于数据融合的井眼轨迹参数动态测量方法可以消除旋转、振动、磁干扰的影响,使测得井眼轨迹参数的精度明显提高,满足了随钻地质导向对井眼轨迹参数测量的要求。
3)建议在现在研究成果基础上,进一步优化数据融合算法,比如基于机器学习的井眼轨迹预测模型,基于遗传算法的卡尔曼滤波器的参数调优。
-
表 1 鹰1井深部泥岩的矿物组成
Table 1 Mineral composition of deep mudstone in Well Ying–1
编号 全岩矿物组成,% 黏土矿物组成,% 石英 长石 方解石 铁白云石 黄铁矿 黏土矿物 伊利石 蒙脱石 伊/蒙混层 绿泥石 高岭石 1 37.10 11.20 6.80 5.40 4.90 30.20 57.10 27.30 3.20 6.50 12.80 2 38.90 12.60 5.20 1.30 5.70 38.40 44.60 24.40 11.90 7.00 13.00 3 45.10 10.40 3.90 4.70 4.20 32.40 58.20 16.80 11.30 4.90 9.30 4 29.90 14.10 7.30 7.60 3.60 30.90 52.60 19.80 4.30 5.50 10.20 5 33.70 13.90 5.20 7.20 4.00 38.50 49.10 25.60 9.50 6.20 14.20 6 38.40 12.80 4.90 5.70 3.80 36.70 43.20 22.80 11.20 7.80 10.20 7 40.30 11.90 3.60 6.30 3.20 39.80 47.10 22.20 8.20 4.60 13.50 8 41.40 12.20 3.20 8.60 2.40 27.50 52.00 19.90 10.70 6.30 8.50 9 39.60 12.60 4.10 3.80 7.80 22.80 44.10 26.40 7.70 6.00 9.70 10 37.90 13.60 4.70 3.90 4.10 40.10 47.80 25.20 8.20 8.20 11.10 平均 38.23 12.53 4.89 5.45 4.37 33.73 49.58 23.04 8.62 6.30 11.25 表 2 鹰1井柯坪塔格组和桑塔木组地层泥岩理化性能
Table 2 Physical and chemical properties of mudstone in Kepingtage Formation and Sangtamu Formation of Well Ying–1
岩样编号 地层 比表面积/(m2·g–1) 总吸水量/(g·g–1) 比亲水量/(mg·m–2) 清水回收率,% 清水膨胀率,% 1 柯坪塔格组 54.12 0.50 9.21 76.35 11.6 2 62.53 0.62 9.93 72.32 13.8 3 50.80 0.51 10.13 80.98 13.2 4 桑塔木组 57.72 0.62 10.72 75.53 14.7 5 51.91 0.49 9.36 73.95 12.8 6 68.92 0.62 9.07 81.73 13.2 注:蒙脱石比亲水量为9.91 mg/m2,伊利石比亲水量为11.62 mg/m2。 表 3 随钻堵漏剂中封堵材料的加量配比
Table 3 Concentration ratio of plugging materials in lost circulation additive while drilling
配方 刚性架桥及充填材料加量,% 弹性可变形封堵材料加量,% 惰性纤维材料加量,% 软化封堵材料加量,% 1 2.0 1.0 1.5 0.5 2 2.0 2.0 0.5 0.5 3 2.0 1.5 1.0 0.5 4 3.0 1.0 0.5 0.5 5 3.0 0.5 1.0 0.5 6 3.0 0.5 0.5 1.0 -
[1] 赵志国,白彬珍,何世明,等. 顺北油田超深井优快钻井技术[J]. 石油钻探技术, 2017, 45(6): 8–13. ZHAO Zhiguo, BAI Binzhen, HE Shiming, et al. Optimization of fast drilling technology for ultra-deep wells in the Shunbei Oilfield[J]. Petroleum Drilling Techniques, 2017, 45(6): 8–13.
[2] 张平. 顺北蓬1井ϕ444.5 mm长裸眼井筒强化钻井液技术[J]. 石油钻探技术, 2018, 46(3): 27–33. ZHANG Ping. Wellbore enhancing technology for ϕ444.5 mm openhole section in Well SHBP1 by means of drilling fluid optimization[J]. Petroleum Drilling Techniques, 2018, 46(3): 27–33.
[3] JARVIE D M, HILL R J, RUBLE T E, et al. Unconventional shale-gas systems: the Mississippian Barnett shale of north-central Texas as one model for thermogenic shale-gas assessment[J]. AAPG Bulletin, 2007, 91(4): 475–499. doi: 10.1306/12190606068
[4] SIMPSON J P, WALKER T O, JIANG G Z. Environmentally acceptable water-base mud can prevent shale hydration and maintain borehole stability[R]. SPE 27496, 1995.
[5] 李佳,邱正松,宋丁丁,等. 井壁强化作用影响因素的数值模拟[J]. 钻井液与完井液, 2017, 34(2): 1–8. doi: 10.3969/j.issn.1001-5620.2017.02.001 LI Jia, QIU Zhengsong, SONG Dingding, et al. Numeric simulation of factors affecting the strengthening of borehole wall[J]. Drilling Fluid & Completion Fluid, 2017, 34(2): 1–8. doi: 10.3969/j.issn.1001-5620.2017.02.001
[6] 于雷,张敬辉,刘宝锋,等. 微裂缝发育泥页岩地层井壁稳定技术研究与应用[J]. 石油钻探技术, 2017, 45(3): 27–31. YU Lei, ZHANG Jinghui, LIU Baofeng, et al. Study and application of borehole stabilization technology in shale strata containing micro-fractures[J]. Petroleum Drilling Techniques, 2017, 45(3): 27–31.
[7] 俞杨烽. 富有机质页岩多尺度结构描述及失稳机理[D]. 成都: 西南石油大学, 2013. YU Yangfeng. Multi-scale structure description and borehole instability mechanism of organic rich shale[D]. Chengdu: Southwest Petroleum University, 2013.
[8] ZHAO Tianyi, LI Xiangfang, ZHAO Huawei, et al. Molecular simulation of adsorption and thermodynamic properties on type Ⅱ kerogen: influence of maturity and moisture content[J]. Fuel, 2016, 190: 198–207.
[9] 牛晓,潘丽娟,甄玉辉,等. SHB1-6H井长裸眼钻井液技术[J]. 钻井液与完井液, 2016, 33(5): 30–34. NIU Xiao, PAN Lijuan, ZHEN Yuhui, et al. Drilling fluid technology for long open hole section of Well SHB1-6H[J]. Drilling Fluid & Completion Fluid, 2016, 33(5): 30–34.
[10] SIGAL R F. A note on the intrinsic porosity of organic material in shale gas reservoir rocks[J]. Petrophysics, 2013, 54(3): 236–239.
[11] 金军斌. 塔里木盆地顺北区块超深井火成岩钻井液技术[J]. 石油钻探技术, 2016, 44(6): 17–23. JIN Junbin. Drilling fluid technology for igneous rocks in ultra-deep wells in the Shunbei Area, Tarim Basin[J]. Petroleum Drilling Techniques, 2016, 44(6): 17–23.
[12] 赵佩,李贤庆,田兴旺,等. 川南地区龙马溪组页岩气储层微孔隙结构特征[J]. 天然气地球科学, 2014, 25(6): 947–956. ZHAO Pei, LI Xianqing, TIAN Xingwang, et al. Study on micropore structure characteristics of Longmaxi Formation shale gas reservoirs in the Southern Sichuan Basin[J]. Natural Gas Geoscience, 2014, 25(6): 947–956.
[13] 薛世峰,马国顺,葛洪魁,等. 液-固-水化耦合形式的井眼稳定性模型研究[J]. 石油钻探技术, 2007, 35(1): 41–44. doi: 10.3969/j.issn.1001-0890.2007.01.012 XUE Shifeng, MA Guoshun, GE Hongkui, et al. Study of a fluid-solid-wetting coupling wellbore stability model[J]. Petroleum Drilling Techniques, 2007, 35(1): 41–44. doi: 10.3969/j.issn.1001-0890.2007.01.012
[14] 宋世超. 泥页岩井壁稳定的力学与化学协同作用研究与应用[D]. 武汉: 长江大学, 2013. SONG Shichao. Mechanical and chemical shale stability research and application of collaborative action[D]. Wuhan: Yangtze University, 2013.
[15] 李芷,贾长贵,杨春和,等. 页岩水力压裂水力裂缝与层理面扩展规律研究[J]. 岩石力学与工程学报, 2015, 34(1): 12–20. LI Zhi, JIA Changgui, YANG Chunhe, et al. Propagation of hydraulic fissures and bedding planes in hydraulic fracturing of shale[J]. Chinese Journal of Rock Mechanics and Engineering, 2015, 34(1): 12–20.
[16] SAVARI S, WHITFILL D L, JAMISON D E, et al. A method to evaluate lost circulation materials-investigation of effective wellbore strengthening applications[R]. SPE 167977, 2014.
[17] 卢运虎,陈勉,安生. 页岩气井脆性页岩井壁裂缝扩展机理[J]. 石油钻探技术, 2012, 40(4): 13–16. doi: 10.3969/j.issn.1001-0890.2012.04.003 LU Yunhu, CHEN Mian, AN Sheng. Brittle shale wellbore fracture propagation mechanism[J]. Petroleum Drilling Techniques, 2012, 40(4): 13–16. doi: 10.3969/j.issn.1001-0890.2012.04.003
[18] 李益寿. 柯193井井壁稳定钻井液技术应用研究[J]. 新疆石油天然气, 2018, 14(2): 37–41. doi: 10.3969/j.issn.1673-2677.2018.02.008 LI Yishou. Study on application for hole stability of Ke 193 Well[J]. Xinjiang Oil & Gas, 2018, 14(2): 37–41. doi: 10.3969/j.issn.1673-2677.2018.02.008
[19] 张金波,鄢捷年,赵海燕. 优选暂堵剂粒度分布的新方法[J]. 钻井液与完井液, 2004, 21(5): 4–7. doi: 10.3969/j.issn.1001-5620.2004.05.002 ZHANG Jinbo, YAN Jienian, ZHAO Haiyan. Optimization of bridging particle size distribution of drilling fluid for formation protection[J]. Drilling Fluid & Completion Fluid, 2004, 21(5): 4–7. doi: 10.3969/j.issn.1001-5620.2004.05.002
[20] 蒋官澄,鄢捷年,王富华,等. 新型屏蔽暂堵技术在大宛齐地区的应用[J]. 石油钻探技术, 1999, 27(6): 21–23. doi: 10.3969/j.issn.1001-0890.1999.06.008 JIANG Guancheng, YAN Jienian, WANG Fuhua, et al. Applications of temporary plugging techniques in Dawanqi Area[J]. Petroleum Drilling Techniques, 1999, 27(6): 21–23. doi: 10.3969/j.issn.1001-0890.1999.06.008
[21] 张金波,鄢捷年. 钻井液暂堵剂颗粒粒径分布的最优化选择[J]. 油田化学, 2005, 22(1): 1–5. doi: 10.3969/j.issn.1000-4092.2005.01.001 ZHANG Jinbo, YAN Jienian. Optimization of particle size distribution for temporarily plugging/shielding agents in water base reservoir drilling fluids[J]. Oilfield Chemistry, 2005, 22(1): 1–5. doi: 10.3969/j.issn.1000-4092.2005.01.001
[22] 舒勇,鄢捷年,宋付英,等. 暂堵剂图解优化新方法在钻井液设计中的应用[J]. 石油钻探技术, 2008, 36(6): 48–51. doi: 10.3969/j.issn.1001-0890.2008.06.011 SHU Yong, YAN Jienian, SONG Fuying, et al. The application of new graphical optimizing method of temporary plugging particle size distribution in drilling fluids design[J]. Petroleum Drilling Techniques, 2008, 36(6): 48–51. doi: 10.3969/j.issn.1001-0890.2008.06.011
-
期刊类型引用(5)
1. 李中,吴怡,李辉,张天玮,焦金刚,邓月,满曰南,袁亚圣,李帅. 海洋智能钻井井下工程参数实时测量技术及工具研发. 钻采工艺. 2025(01): 63-69 . 百度学术
2. 蔡峰,朱美静. 基于EKF和UKF的随钻姿态解算方法研究. 安徽理工大学学报(自然科学版). 2024(01): 12-20 . 百度学术
3. 郝世俊,褚志伟,李泉新,方俊,陈龙,刘建林. 煤矿井下近钻头随钻测量技术研究现状和发展趋势. 煤田地质与勘探. 2023(09): 10-19 . 百度学术
4. 潘勇,李国光,柳敬云. 基于国产可编程逻辑门阵列的近钻头随钻测量系统设计. 钻探工程. 2023(S1): 410-413 . 百度学术
5. 蒲文学,范光第,朱建建,赵国山. 磁性随钻测斜仪所需无磁钻具长度及影响因素研究. 石油钻探技术. 2022(04): 129-134 . 本站查看
其他类型引用(5)