基于数据融合的近钻头井眼轨迹参数动态测量方法

宋晓健， 郑邦贤， 谭勇志， 黄秉亚， 马鸿彦， 董晨曦

（中国石油集团渤海钻探工程有限公司定向井技术服务分公司，天津300280）

现代地质导向钻井过程中，需采用近钻头随钻边界探测技术动态测量近钻头井眼轨迹参数（井斜角、方位角和工具面角），根据随钻成像预估地层和井眼变化趋势，据此调整钻头方向，提高油层钻遇率。由于振动、旋转、磁干扰等因素的影响，导致测量的近钻头井眼轨迹参数出现很大误差，严重影响了随钻测量精度，出现钻穿油层等情况。

目前旋转导向动态测量技术都掌握在国外油田技术服务公司手中，不对外公开，因此，未见到国外有井眼轨迹参数动态测量技术的报道。国内杨全进等人[1]建立了一种旋转导向系统有色噪声的改进无迹卡尔曼滤波方法，该算法限制条件较多。高怡等人[2]提出了采用多源动态姿态组合测量方法测量导向钻具的动态姿态，但加速度计采集数据的存储周期较长，造成测量误差较大。Xue Qilong 等人[3–4]提出了基于卡尔曼滤波状态空间模型的动态井眼轨迹测量方法，但其只适用于线性系统。徐宝昌等人[5]提出了基于无迹卡尔曼滤波的动态姿态测量方法，但没有解决振动及旋转对动态姿态测量的影响。为此，笔者提出了一种基于数据融合的近钻头井眼轨迹参数动态测量方法，该方法针对三轴加速度计、磁通门和速率陀螺的测量系统，建立了基于四元数井眼轨迹参数测量模型，推导出加速度计与井眼轨迹平滑预测的关系，运用3个捷联式无迹卡尔曼滤波器和磁干扰校正系对加速度计、磁通门进行消噪，校正，实时测量近钻头井眼轨迹参数，提高随钻近钻头探边能力，确定最佳储层位置。

1 近钻头动态井眼轨迹测量模型

近钻头动态测量系统由三轴加速度计、三轴磁通门和角速率陀螺仪组成，针对该系统，基于四元数方法建立井眼轨迹非线性数学模型。如图1所示，根据地理坐标系O-NED和钻具坐标系O-xyz的对应关系，建立欧拉角转换矩阵，并转换为四元数，k时刻姿态转换矩阵T表示为：

图1 井眼轨迹参数测量原理Fig.1 Measurement principle of borehole trajectory parameters

式中：q0，q1，q2和q3为四元数。

理论上由于x轴和y轴的加速度计、磁通门正交相隔90°，其信号为正弦或余弦曲线。由于钻头旋转、振动的干扰，近钻头x轴和y轴加速度计、磁通门的实测信号如图2所示。由图2可以看出，近钻头的旋转、振动对磁通门的影响相对于加速度计较小，所以采取磁通门读数校正加速度计读数。

图2 近钻头磁通门、加速度计实测数据Fig.2 Measured data of near-bit fluxgates and accelerometers

式中：g为重力加速度，m/s2；B为地磁场强度，μT；Bx(k)，By(k)和Bz(k)为k时刻三轴磁通门测量的磁场强度，μT；ax(k)，ay(k)和az(k)为k时刻三轴加速度计测量结果，m/s2；wx(k)，wy(k)和wz(k)为k时刻陀螺仪的角速度，rad/s。

2 数据融合近钻头井眼轨迹参数动态测量方法

基于数据融合算法的近钻头井眼轨迹参数动态测量方法的测量流程如图3所示，图中KF1和KF2为基于扩展卡尔曼滤波算法的滤波器，KF3为基于无迹卡尔曼滤波算法的滤波器。测量步骤：1）将加速度计、磁通门、转动角速度四元数带入KF1滤波器，进行扩展卡尔曼滤波，得出井斜角、方位角估计值；2）将加速度计四元数带入KF2滤波器，进行扩展卡尔曼滤波，得出测深增量Δhm；3）将测深增量Δhm、井斜角、方位角估计值带入KF3滤波器，进行无迹卡尔曼滤波，得出井斜角、方位角最终估计值；4）利用井斜角、方位角最终估计值计算磁性工具面角ωm与重力工具面角的差Δω；5）利用磁性工具面角和角差Δω求出重力工具面角ωg。

图3 数据融合测量过程Fig.3 Measurement process of data fusion

近钻头振动信号是一种幅值大、频率高、频带宽的噪声信号，可以近似等效为高斯白噪声。动态测量的动力学模型为典型的非线性模型。根据z轴陀螺仪测量的转速wz，运用转速补偿策略，对x轴和y轴加速度计测量结果进行转速补偿[5]，减小钻具旋转对加速度计测量结果的影响。根据转速补偿结果以及噪声特性，采用扩展、无迹卡尔曼滤波滤除振动干扰信号。

2.1 估计近钻头井斜角、方位角的扩展卡尔曼滤波算法

基于四元数的KF1的状态方程和量测方程：

式中：Q(k)为k时刻的状态值；I为单位矩阵；ts为采样周期；w(k)为k时刻系统高斯白噪声；v(k)为k时刻传感器观测噪声；A(k)为k时刻状态转移矩阵；F(x)为非线性函数；Z(k+1)为k+1时刻的观测值。

三轴加速度信号、三轴磁通门信号、角速率陀螺信号进行数据融合后，采用扩展卡尔曼滤波算法，得到最优姿态估计[6–7]，动态解算出钻井工具的实时姿态参数，确保钻具姿态测量计算的精度，减少计算量，对四元数Q进行更新，求出KF1滤波后的井斜角αKF1、方位角ϕKF1、高边工具面角ωg,KF1和磁性工具面角ωm,KF1。

式中：αKF1，ϕKF1，ωg,KF1和ωm,KF1分别为 KF1 滤波后的井斜角、方位角、重力工具面角和磁性工具面角，（°）；θ为地层倾角，（°）。

2.2 估计近钻头测深增量的扩展卡尔曼滤波算法

根据式（7）计算出四元数中的az，运用扩展卡尔曼滤波器计算系统经过ts后测深增量Δhm。z轴加速度计主要受到重力加速度和振动的干扰，由于采样时间ts为毫秒级，在单位采样周期内，重力加速度和振动的干扰可以视为近似相同[8–10]，可以忽略振动对加速度计测量结果的影响。k为当前采样点，z轴加速度增量Δaz：

为了提高对测深增量的估计，对Δhm进行二阶泰勒展开：

KF2的状态方程和量测方程为：

2.3 估计近钻头井眼轨迹参数的无迹卡尔曼滤波算法

如图4所示，在单位采样时间内，井眼轨迹趋于平滑曲线，可以根据前面2个测点的狗腿度和KF2输出测深增量对井眼轨迹进行递归式预测[11]。

图4 KF3预测井深参数的原理Fig.4 Prediction principle of well depth parameters by KF3

式中：γ为狗腿角，（°）；Δγ为系统经过ts的狗腿角增量，（°）；hm(k)为k时刻的测深，m。

根据 ∆hm(k)进行无迹卡尔曼滤波，KF3的状态方程和量测方程为：

式中：wα和wϕ分别为井斜角和方位角的系统高斯白噪声。

式中：αKF3和ϕKF3分别为KF3滤波后的井斜角和方位角，（°）；vα和vϕ分别为井斜角和方位角的系统观测噪声。

2.4 近钻头重力工具面角的估计

根据旋转测量原理，同一时刻的重力工具面角与磁工具面角的差与测量时刻的井斜角、方位角、地磁倾角呈现一定函数关系[12–13]。根据KF3求出的井眼井斜角和方位角计算磁性工具面角与重力工具面角的差Δω：

根据Δω，计算旋近钻头动态重力工具面角估计值ωdg,e：

式中：ωdg,e为旋近钻头动态重力工具面角估计值，（°）；ωm,KF3为 KF3 滤波后的磁性工具面角，（°）。

2.5 磁干扰情况下的磁性工具面角

近钻头重力工具面角需根据x轴和z轴磁通门传感器的读数求出[14–15]，而磁通门测量结果不可避免地会受到周围电磁场的影响，为了降低井下钻具周围电磁场的影响程度，采用径向磁干扰的方法校正x轴和z轴磁通门的测量结果[16]。

磁场的干扰导致磁通门测量的磁场强度发生偏移和变形。磁干扰下的测量结果如图5所示。

图5 磁干扰校正示意Fig.5 Magnetic interference correction

在实际钻井过程中，井下仪器旋转一圈时，钻深可以忽略不计，可以看作仪器在原地旋转了一圈。z轴磁通门的测量结果可以认为没有发生变化，而x轴和y轴磁通门的测量值不断发生变化，如图2（a）所示。三轴磁通门传感器的测量数据记为（Bx，By，Bx），地球磁场可以看成一个固定值，即：

式中 ：C为常数。

根据椭圆校正原理，对短时间内采集的Bx，By进行磁干扰校正，得出排除磁干扰的Bxm和Bym：

式中：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°。

图6 测量井斜角、方位角与实际井斜角、方位角的对比Fig.6 Comparison of measured and actual deviation angles and azimuths

图7为采用上文测量方法（以下称为动态测量方法）测量工具面角、NWD探管直接测量工具面角与实际工具面角的差值。由图7可以看出，NWD探管直接测得工具面角与实际工具面角的差最大接近40°，已经不能满足测量近钻头伽马重力工具面角的要求，而采用动态测量方法测得的工具面角与实际工具面角的差小于10°，满足测量近钻头工具面角的要求。

图7 动态测量和直接测量工具面角与实际工具面角的差值Fig.7 Difference between actual tool face angles and tool face angles measured by dynamic and direct methods

在井斜角为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.934 2°。图9 为动态测得方位角与静态测得方位角的对比。由图9可以看出，动态测得方位角与静态测得方位角的平均误差为 1.19°，方差为 1.943 642°。

图8 动态测量井斜角与静态测量井斜角的对比Fig.8 Comparison of deviation angles measured by dynamic and static methods

图9 动态测量方位角与静态测量方位角的对比Fig.9 Comparison of azimuths measured by dynamic and static methods

图10为动态测量工具面角与静态测量工具面角的对比。由图10可以看出，动态测得工具面角与静态测得工具面角的误差小于5°。

图10 动态测量工具面角与静态测量工具面角的对比Fig.10 Comparison of tool face angles measured by dynamic and static methods

以上分析可以得出，基于数据融合的井眼轨迹参数测量方法能够有效消除旋转、振动、磁干扰的影响，井斜角、方位角和工具面角的测量精度得到明显提高，满足了随钻地质导向对井眼轨迹参数测量的要求。

4 结 论

1）针对旋转、振动、磁干扰对近钻头测量仪器的影响，提出了基于数据融合的井眼轨迹参数动态测量方法。该方法采用捷联式卡尔曼滤波器和磁干扰校正系统对测量信号进行滤波、校正，利用四元数法处理测量数据，求出近钻头井眼轨迹参数。

2）实验室模拟试验和现场试验均表明，基于数据融合的井眼轨迹参数动态测量方法可以消除旋转、振动、磁干扰的影响，使测得井眼轨迹参数的精度明显提高，满足了随钻地质导向对井眼轨迹参数测量的要求。

3）建议在现在研究成果基础上，进一步优化数据融合算法，比如基于机器学习的井眼轨迹预测模型，基于遗传算法的卡尔曼滤波器的参数调优。