一种基于实际采集传感器数据的钻铤姿态测量方法与流程

本发明涉及测量领域，特别涉及一种基于实际采集传感器数据的钻铤姿态测量方法。

背景技术：

旋转自动导向技术是上世纪九十年代发展起来的一项高新钻井技术。该技术可以实现对旋转钻井工具的几何导向控制，使实钻井眼轨迹与设计经验轨迹一致，可以提高钻井效率，降低钻井成本。国外各大石油服务公司大都研制了自己的旋转导向钻井系统，国内的很多单位也进行了相关研究，但总体和国外仍有较大差距。

旋转自动导向技术的关键在钻铤姿态的计算，钻铤工作在较为恶劣的环境下，各个传感器的输出数据随时间而积累可能会产生无法估计的误差，，捷联系统的姿态解算原理决定了传感器输出误差会严重影响钻铤姿态角的解算。因此提高钻铤姿态角的解算稳定性成为解决旋转导向技术的一个重要的手段。

技术实现要素：

本发明要解决的技术问题是提供一种基于实际采集传感器数据的钻铤姿态测量方法，能够提高钻铤磁工具面角、井斜角、井斜工具面角解算稳定性。

为了解决上述技术问题，本发明的技术方案为：

一种基于实际采集传感器数据的钻铤姿态解算方法，包括以下步骤：

(a)根据钻铤的运动过程及传感器数据特性，对加速度计数据进行低通滤波；

(b)使用钻铤静止时数据估计加速度计零偏；

(c)根据(b)结果补偿加速计零偏，对钻铤旋转运动阶段加速度计输出进行修正；

(d)计算相应时间序列下钻铤的磁工具面角mhs及重力在重力在钻铤的投影；

(e)计算相应时间序列下的井斜工具面角incmhs和井斜角inc。

优选地，在步骤(a)中，所述加速度计的有效数据集中于低频6hz之内。

优选地，在步骤(a)中，选用切比雪夫i型滤波器对所述数据进行滤波，其中通带上限频率为ωp＝15hz，阻带下限频率为ωs＝20hz，通带和阻带对应的最大衰减系数分别为：ap＝0.5db,as＝20db。

优选地，在步骤(a)中，其中所述切比雪夫i型滤波器的传递函数为

优选地，在步骤(b)中钻铤运动过程可知前20s内钻铤处于静止状态。

优选地，在步骤(b)中的零偏估算过程为：对于加速度计在钻铤运动过程中不可避免的存在振动加速度av，由传感器配置，x轴和y轴上都分别反正有两个加速度计，旋转中加速度计测量到的比力可以表示为：

其中，ax和ay分别是x轴和y轴真实要测量的加速度，av,x和av,y分别是x轴和y轴加速度计测得的振动加速度，和为加速计每次启动后常值零偏，随机零偏和为白噪声，联立不定方程组，并在线性空间内拟合求可对加速度计零偏进行估算。

优选地，在步骤(d)中，计算方式为：

优选地，在步骤(d)中，重力在钻铤的投影gp：

优选地，在步骤(e)中，计算方式为：

采用上述技术方案，使得本发明的有益效果在于：

(1)根据钻铤的运动过程及传感器数据特性，通过仿真分析，设计五阶切比雪夫低通滤波器，对传感器输出数据进行预处理，提高数据有效性的和精度；

(2)根据钻铤的运动过程，利用前20s内钻铤处于静止状态数据，对加速度计零偏进行估算，并根据加速度计零偏估算结果，对钻铤旋转运动阶段加速度计输出进行补偿；

(3)根据钻铤不同姿态及磁工具面角mhs与和关系，建立磁工具面角mhs真值解算表。

附图说明

图1为x轴磁通门局部数据曲线

图2为x正轴加速度计数据曲线

图3为数字低通滤波器的幅频响应曲线

图4为x正轴加速度计滤波效果曲线

图5为前20s磁工具面角mhs解算结果

图6为处理前后旋转阶段井斜角inc解算结果

图7为处理前后旋转阶段井斜工具面角incmhs解算结果

具体实施方式

下面结合附图对本发明的具体实施方式作进一步说明。在此需要说明的是，对于这些实施方式的说明用于帮助理解本发明，但并不构成对本发明的限定。此外，下面所描述的本发明各个实施方式中所涉及的技术特征只要彼此之间未构成冲突就可以相互组合。

本发明提供了一种测量系统，所述测量系统包括一实际采集传感器和一滤波器，所述实际采集传感器采集数据并处理，所述滤波器对所述实际采集传感器的输出数据进行预处理，提高了数据的有效性和精度。所述实际采集传感器包括多个加速度计和一磁通门，通过所述加速度计和所述磁通门的数据进行数据采集和数据处理。

根据所述测量系统所实施的测量过程包括以下步骤：

第一步、首先对加速度计数据低通滤波，提高数据有效性的和精度；

第二步、使用钻铤静止时数据估计加速度计零偏；

第三步、根据s2结果补偿加速度计零偏，对钻铤旋转运动阶段加速度计输出进行修正；

第四步、计算相应时间序列下钻铤的磁工具面角mhs及重力在重力在钻铤的投影；

第五步、计算相应时间序列下的井斜工具面角incmhs和井斜角inc。

根据上述过程，实际编排中在x轴、y轴与z轴的正半轴和负半轴各安装有一个加速度计，磁通门采用冗余设计。

采集信息如下：

(1)钻铤x、y、z轴加速度表分量(“+”和“－”分别表示坐标系的正轴和负轴)；

(2)磁通门x、y、z轴分量

其中x轴磁通门局部数据曲线、x正轴加速度计数据曲线分别如图2、图3所示。

具体地，在第一步中，根据所述测量系统应用背景，加速度计有效数据集中于低频6hz以内，存在频带比较宽的高频噪声，因此前期对加速度计数据进行低通滤波可以提高数据有效性的和精度，截止频率处要求较为陡降，而且对于数字滤波器来说，较低的阶数才能使得通带相位接近于线性，也便于计算，同时对于通带的纹波可以不必过多考虑。

具体地，在第二步中，根据上述的选定原则，本发明选定切比雪夫i型滤波器，选定通带上限频率为ωp＝15hz，阻带下限频率为ωs＝20hz，通带和阻带对应的最大衰减系数分别为：ap＝0.5db,as＝20db。

通过仿真分析，使用五阶滤波器可满足上述要求，设计的切比雪夫滤波器传递函数为：

其中数字低通滤波器的幅频响应如图4所示，x正轴加速度计滤波效果如图5所示。根据钻铤运动过程可知前20s内钻铤处于静止状态，再利用20s内采集到的传感器数据，计算磁工具面角mhs，则mhs平均值可作为此时钻铤真实磁工具面角。

在上述过程中，根据钻铤运动过程可知前20s内钻铤处于静止状态，再利用20s内采集到的传感器数据，计算磁工具面角mhs，则mhs平均值可作为此时钻铤真实磁工具面角。另一方面，对于加速度计在钻铤运动过程中不可避免的存在振动加速度，由传感器配置，x轴和y轴上都分别反正有两个加速度计，旋转中加速度计测量到的比力可以表示为：

其中，ax和ay分别是x轴和y轴真实要测量的加速度，av,x和av,y分别是x轴和y轴加速度计测得的振动加速度，和为加速计每次启动后常值零偏，可认为都是恒值，随机零偏和为白噪声，联立不定方程组，并在线性空间内拟合求可对加速度计零偏进行估算。

在第四步中，根据磁工具面角mhs与和关系，可得mhs计算公式如表1所示的mhs的真值解。其中前20s磁工具面叫mhs解算结果如图6所示。

具体地，根据磁工具面角mhs与和关系，可得mhs计算公式如表1所示：

表1mhs的真值解

进一步地，计算方式为：

其中，重力在钻铤的投影gp：

在第五步中，计算方式为：

利用本发明公开的方法对实际采集数据进行处理，通过对比分析，发现本发明可提高井斜角inc、井斜工具面角incmhs解算稳定性。处理前后旋转阶段井斜角inc、井斜工具面角incmhs解算结果如图7所示。

根据上述的实施内容，本发明进一步地提供了一种钻铤姿态解算方法，包括以下步骤：

(a)根据钻铤的运动过程及传感器数据特性，对加速度计数据进行低通滤波；

(b)使用钻铤静止时数据估计加速度计零偏；

(c)根据(b)结果补偿加速计零偏，对钻铤旋转运动阶段加速度计输出进行修正；

(d)计算相应时间序列下钻铤的磁工具面角mhs及重力在重力在钻铤的投影；

(e)计算相应时间序列下的井斜工具面角incmhs和井斜角inc。

根据上述钻铤姿态解算方法，在步骤(a)中，所述加速度计的有效数据集中于低频6hz之内。

根据上述钻铤姿态解算方法，在步骤(a)中，选用切比雪夫i型滤波器对所述数据进行滤波，其中通带上限频率为ωp＝15hz，阻带下限频率为ωs＝20hz，通带和阻带对应的最大衰减系数分别为：ap＝0.5db,as＝20db。

根据上述钻铤姿态解算方法，在步骤(a)中，其中所述切比雪夫i型滤波器的传递函数为

根据上述钻铤姿态解算方法，在步骤(b)中钻铤运动过程可知前20s内钻铤处于静止状态。

根据上述钻铤姿态解算方法，在步骤(b)中的零偏估算过程为：对于加速度计在钻铤运动过程中不可避免的存在振动加速度av，由传感器配置，x轴和y轴上都分别反正有两个加速度计，旋转中加速度计测量到的比力可以表示为：

其中，ax和ay分别是x轴和y轴真实要测量的加速度，av,x和av,y分别是x轴和y轴加速度计测得的振动加速度，和为加速计每次启动后常值零偏，随机零偏和为白噪声，联立不定方程组，并在线性空间内拟合求可对加速度计零偏进行估算。

根据上述钻铤姿态解算方法，在步骤(d)中，计算方式为：

根据上述钻铤姿态解算方法，在步骤(d)中，重力在钻铤的投影gp：

根据上述钻铤姿态解算方法，在步骤(e)中，计算方式为：

具体地，如图1和图2所示的x轴磁通门局部数据曲线、x正轴加速度计数据曲线。

如图3所示的数字低通滤波器的幅频响应。如图4所示的x正轴加速度计滤波效果；如图5所示的其中前20s磁工具面叫mhs解算结果。其中根据磁工具面角mhs与和关系，可得mhs计算公式如表1所示的mhs的真值解。

利用本发明公开的方法对实际采集数据进行处理，通过对比分析，发现本发明可提高井斜角inc、井斜工具面角incmhs解算稳定性。处理前后旋转阶段井斜角inc、井斜工具面角incmhs解算结果如图6、图7所示。

本发明未详细说明部分属本领域技术人员公认常识。

以上结合附图对本发明的实施方式作了详细说明，但本发明不限于所描述的实施方式。对于本领域的技术人员而言，在不脱离本发明原理和精神的情况下，对这些实施方式进行多种变化、修改、替换和变型，仍落入本发明的保护范围内。