CN116607929A - 一种煤田钻探深孔钻探中的钻孔自动纠偏方法 - Google Patents
一种煤田钻探深孔钻探中的钻孔自动纠偏方法 Download PDFInfo
- Publication number
- CN116607929A CN116607929A CN202310657569.4A CN202310657569A CN116607929A CN 116607929 A CN116607929 A CN 116607929A CN 202310657569 A CN202310657569 A CN 202310657569A CN 116607929 A CN116607929 A CN 116607929A
- Authority
- CN
- China
- Prior art keywords
- drilling
- quaternion
- attitude
- state
- sensor
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000005553 drilling Methods 0.000 title claims abstract description 118
- 238000000034 method Methods 0.000 title claims abstract description 39
- 239000003245 coal Substances 0.000 title claims abstract description 29
- 238000012937 correction Methods 0.000 claims description 49
- 206010034719 Personality change Diseases 0.000 claims description 27
- 230000036544 posture Effects 0.000 claims description 19
- 238000004422 calculation algorithm Methods 0.000 claims description 16
- 238000004364 calculation method Methods 0.000 claims description 16
- 239000011159 matrix material Substances 0.000 claims description 12
- 230000005540 biological transmission Effects 0.000 claims description 9
- 230000008569 process Effects 0.000 claims description 9
- 238000007781 pre-processing Methods 0.000 claims description 7
- 238000005259 measurement Methods 0.000 claims description 6
- 238000001914 filtration Methods 0.000 claims description 5
- 238000012544 monitoring process Methods 0.000 claims description 5
- 238000012545 processing Methods 0.000 claims description 5
- 238000005516 engineering process Methods 0.000 claims description 4
- 238000006243 chemical reaction Methods 0.000 claims description 3
- QVFWZNCVPCJQOP-UHFFFAOYSA-N chloralodol Chemical compound CC(O)(C)CC(C)OC(O)C(Cl)(Cl)Cl QVFWZNCVPCJQOP-UHFFFAOYSA-N 0.000 claims description 3
- 238000005094 computer simulation Methods 0.000 claims description 3
- 238000012790 confirmation Methods 0.000 claims description 3
- 238000009825 accumulation Methods 0.000 description 8
- 238000001514 detection method Methods 0.000 description 4
- 238000005065 mining Methods 0.000 description 3
- 238000010276 construction Methods 0.000 description 2
- 230000006378 damage Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 239000011148 porous material Substances 0.000 description 2
- 239000011435 rock Substances 0.000 description 2
- 229910000831 Steel Inorganic materials 0.000 description 1
- 238000005452 bending Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004140 cleaning Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 239000002002 slurry Substances 0.000 description 1
- 230000000087 stabilizing effect Effects 0.000 description 1
- 239000010959 steel Substances 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B44/00—Automatic control systems specially adapted for drilling operations, i.e. self-operating systems which function to carry out or modify a drilling operation without intervention of a human operator, e.g. computer-controlled drilling systems; Systems specially adapted for monitoring a plurality of drilling variables or conditions
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
- E21B47/02—Determining slope or direction
- E21B47/024—Determining slope or direction of devices in the borehole
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B7/00—Special methods or apparatus for drilling
- E21B7/04—Directional drilling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Computational Mathematics (AREA)
- Fluid Mechanics (AREA)
- Environmental & Geological Engineering (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- Operations Research (AREA)
- Computing Systems (AREA)
- Geophysics (AREA)
- Earth Drilling (AREA)
Abstract
本发明公开了一种煤田钻探深孔钻探中的钻孔自动纠偏方法,涉及煤田钻探技术领域,本发明将通过四元姿态解算算法对纠偏姿态进行确定,同时通过基于卡尔曼滤波器建立的钻孔姿态估计问题建模避免出现姿态四元数的漂移和误差累积问题,有效抑制传感器数据中的漂移和误差累积,通过偏差计算将实际测得的钻孔姿态与预设方向进行比较,确定实际方向与预设方向之间的偏差,降低了操作人员在危险环境中的风险和受伤可能性,实现自动化的钻孔偏差检测和纠正,减少人工干预,提高钻探作业效率,降低人工操作误差。
Description
技术领域
本发明涉及煤田钻探技术领域,具体为一种煤田钻探深孔钻探中的钻孔自动纠偏方法。
背景技术
煤田钻探是指在煤矿地区进行的一种勘探和开采前的钻探活动。它主要用于获取有关煤矿地下地质情况的信息,包括煤层的分布、厚度、品质、倾角等参数。这些信息对于煤矿开采的规划和设计非常重要。煤田钻探的主要目的是确定煤田的储量和质量,以便进行矿井设计和开采方案的制定。通过钻探取样和测量,可以确定煤层的物理性质、化学成分以及其他地质参数,为后续的煤矿开采工作提供依据。在煤田钻探中,常用的钻探方法包括钻孔钻探和岩心钻探。钻孔钻探是通过旋转钻杆和钻头来打开地下岩石,然后用钻孔来获取地质样品和测量地下地质参数。岩心钻探则是使用特殊的岩心钻头来采集岩石样本,以获取更准确的地质信息。
但是在现有技术中,钻孔的实际钻取参数要与预设钻取参数相近,例如中国专利公开了一种煤田钻探深孔钻探中的钻孔自动纠偏方法,CN114482861A,包括以下步骤:S1、保证场地平整,使得储备泥浆高出施工水位,防止钻孔过程中由于孔内压力小于孔外压力而导致坍孔;S2、在确认好钻孔位置后,将钻机移动到施工位置处;S3、操作钻机起落钢丝绳中心应对准桩中心;S4、钻机就位后校正,底座平整稳固,确保在钻进中不发生倾斜和位移,保证钻进中钻具的平稳及钻孔质量;S5、在成孔前,安装垂直度监测系统和钻头纠偏系统,自动纠偏后再成孔操作;S6、清孔后,对钻孔进行阶段性的成孔质量检查;S7、对孔内进行注浆,完成钻孔。在成孔前,钻头纠偏系统将钻杆垂直度纠正回正常范围内,在钻孔的过程中自动纠偏,提高钻孔质量。
虽然上述方案具有如上的优势,但是针对煤田钻探的深孔钻探,在实际的钻探过程中,引起钻孔偏斜的因素是多方面的,主要分为力学弯曲因素和集合偏斜因素,而无论哪种因素导致的钻孔偏斜,传统人工钻探都无法进行及时调整,还可能出现煤田钻探事故,导致人员伤亡、设备损坏,而少有的带有纠偏功能的钻头系统如上述技术方案,采用钻头检测系统和钻头纠偏系统来对钻头钻孔时进行监测,可以在出现偏斜进行及时的纠偏,但是在实际对钻头偏移检测过程,需要由传感器来对钻孔各参数进行监控,由于无法有效降低或减少传感器的误差累计,即使进行了纠正,还会因为传感器的累积误差进行二次纠正,三次纠正,因此亟须一种煤田钻探深孔钻探中的钻孔自动纠偏方法来解决此类问题。
发明内容
本发明的目的是为了解决现有技术中存在的无法有效降低或减少传感器的误差累计,即使进行了纠正,还会因为传感器的累积误差进行二次纠正,三次纠正的问题。
为实现上述目的,本发明提供了如下技术方案:
本发明提供了一种煤田钻探深孔钻探中的钻孔自动纠偏方法,其特征在于,包括以下步骤:
(1)纠偏数据采集
在钻杆中部安装倾角传感器,在钻杆下部安装方位角传感器,在钻杆上部安装旋转速度传感器,在钻机臂上同样安装倾角传感器和方位角传感器,多个传感器持续对钻杆钻探时的倾角、方位角、旋转速度进行监测,并进行记录;
(2)纠偏数据处理
对钻探时的倾角、方位角、旋转速度数据进行预处理;
(3)纠偏姿态计算
采用四元姿态解算算法对纠偏姿态进行确定,具体的:初始化姿态解算、传感器数据更新,其中传感器更新具体步骤为:
3.1从传感器数据中获取角速度数据;
3.2由当前时刻的获取的旋转速度数据采用四阶龙格库塔法来估计姿态变化四元数,具体的:
3.2.1将旋转速度传感器更新频率定义为△t,即离散化的时间步长;
3.2.2由当前时刻的角速度数据,计算姿态变化四元数的增量△q1;
3.2.3将当前时刻的角速度数据乘以△t的一半,得到中间时刻的角速度:
3.2.4由中间时刻的角速度,计算中间时刻的姿态变化四元数的增量△q2;
3.2.5再次将中间时刻的角速度乘以△t的一半,得到另一个中间时刻的角速度:
3.2.6由另一个中间时刻的角速度,计算另一个中间时刻的姿态变化四元数的增量△q3;
3.2.7将三个姿态变化四元数增量并进行加权平均,得到最终的姿态变化四元数增量:△q=(△q1+2×△q2+2×△q3+△q4)/6,其中△q4为使用当前时刻的角速度计算的姿态变化四元数的增量;
3.2.8将姿态变化四元数增量与上一时刻的姿态四元数相乘,得到当前时刻的姿态四元数:q(t)=q(t-1)×△q;
3.3将步骤2.8中得到的姿态变化四元数与上一时刻的四元数相乘,得到当前时刻的姿态四元数;
3.4在步骤(3),对姿态四元数进行归一化,使其范数等于1;
(4)姿态结果输出转换
根据姿态四元数,计算出钻孔的倾角和方位角,具体的
4.1通过姿态四元数转换为欧拉角,对倾斜角度进行表示;
4.2通过姿态四元数转换为旋转矩阵,姿态四元数q的方向向量表示为:D=[2(q1q3-q2q0),2(q2q3-q1q0),q0 2-q1 2-q2 2+q3 2],根据方向向量D计算方位角,采用反正切函数计算方向向量的x分量与y分量的比值,即atan2(D[1],D[0]);
(5)钻孔姿态估计问题建模
基于卡尔曼滤波器建立钻孔姿态估计问题建模,具体的:
5.1将钻孔姿态估计问题建立为一个状态空间模型,建立状态方程和观测方程;
5.2初始化卡尔曼滤波器的状态向量和协方差矩阵;
5.3根据钻孔当前姿态估计问题建模中的状态方程,使用过程模型对当前状态进行预测,并计算预测状态的协方差矩阵,预测步骤利用先验信息预测下一时刻的状态;
5.4根据系统建模中的观测方程,将各传感器测量值与预测状态进行比较,计算预测状态的测量残差,并计算卡尔曼增益;
5.5利用卡尔曼增益修正预测状态和协方差矩阵,得到更新后的状态估计和协方差;
(6)实际纠偏值计算
将实际测得的钻孔姿态通过向量运算与预设方向进行比较,确定偏差;
(7)纠偏校正值确认
根据计算得到的偏差值,采用PID算法进行校正;
(8)纠偏校正指令输出
根据PID算法的输出,生成对应控制指令,调整钻机的工作参数;
(9)纠偏后告警
持续监测钻孔的姿态数据,根据纠偏结果提供实时反馈。
本发明进一步地设置为:在步骤(1)中,基于LoRa无线传输技术建立无线传输网络拓扑结构,以多个传感器节点组成网络,由中继点进行纠偏数据传输;
本发明进一步地设置为:在步骤(2)中,所述预处理包括校准和滤波;
本发明进一步地设置为:在步骤(3)中,所述初始化姿态解算将四元数进行初始化,并将其设为初始值[1,0,0,0]来表示钻杆非旋转状态即初始化状态;
本发明进一步地设置为:在步骤4.1中,所述欧拉角为俯仰角、横滚角和航向角,其中俯仰角表示钻孔的倾斜角度;
本发明进一步地设置为:在步骤5.1中,钻孔的真实姿态定义为状态变量,各传感器测量的姿态数据定义为观测变量,并定义状态变量和观测变量之间的关系;
本发明进一步地设置为:在步骤(9)中,所述实施反馈具体为:当偏差超过预设阈值时,进行告警;当纠偏后需要二次纠偏时,进行告警。
采用本发明提供的技术方案,与已知的公有技术相比,具有如下有益效果:
本发明通过四元姿态解算算法对纠偏姿态进行确定,同时通过基于卡尔曼滤波器建立的钻孔姿态估计问题建模避免出现姿态四元数的漂移和误差累积问题,其中卡尔曼增益将测量残差与预测步骤中的状态预测进行融合,得到最优估计的姿态,卡尔曼滤波不断地对状态进行预测和更新,将测量值和先验信息进行融合,有效抑制传感器数据中的漂移和误差累积,并提供更精确的姿态估计结果,再通过偏差计算将实际测得的钻孔姿态与预设方向进行比较,确定实际方向与预设方向之间的偏差,计算获取具体偏差值,采用PID算法进行校正,来调整钻机工作参数,从而实现自动校正,同时在调整后对钻杆状态进行监测,当偏差超过预设阈值时,进行告警,当纠偏后需要二次纠偏时,进行告警,从而及时提醒操作人员干预,避免出现钻卡事故,降低了操作人员在危险环境中的风险和受伤可能性,区别传统深孔钻探,传感器数据的实时采集和处理使得钻孔偏差的检测和纠正可以在钻进过程中及时进行,同时采用四元姿态解算算法对纠偏姿态进行确定并通过卡尔曼滤波器建立的钻孔姿态估计问题建模避免出现姿态四元数的漂移和误差累积问题,实现自动化的钻孔偏差检测和纠正,减少人工干预,提高钻探作业效率,降低人工操作误差,解决了现有技术中存在的无法有效降低或减少传感器的误差累计,即使进行了纠正,还会因为传感器的累积误差进行二次纠正,三次纠正的问题。
具体实施方式
下面将对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施条例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明提供一种技术方案:一种煤田钻探深孔钻探中的钻孔自动纠偏方法,其特征在于,包括以下步骤:
(1)纠偏数据采集
在钻杆中部安装倾角传感器,在钻杆下部安装方位角传感器,在钻杆上部安装旋转速度传感器,在钻机臂上同样安装倾角传感器和方位角传感器,多个传感器持续对钻杆钻探时的倾角、方位角、旋转速度进行监测,并进行记录;
(2)纠偏数据处理
对钻探时的倾角、方位角、旋转速度数据进行预处理;
(3)纠偏姿态计算
采用四元姿态解算算法对纠偏姿态进行确定,具体的:初始化姿态解算、传感器数据更新,其中传感器更新具体步骤为:
3.1从传感器数据中获取角速度数据;
3.2由当前时刻的获取的旋转速度数据采用四阶龙格库塔法来估计姿态变化四元数,具体的:
3.2.1将旋转速度传感器更新频率定义为△t,即离散化的时间步长;
3.2.2由当前时刻的角速度数据,计算姿态变化四元数的增量△q1;
3.2.3将当前时刻的角速度数据乘以△t的一半,得到中间时刻的角速度:
3.2.4由中间时刻的角速度,计算中间时刻的姿态变化四元数的增量△q2;
3.2.5再次将中间时刻的角速度乘以△t的一半,得到另一个中间时刻的角速度:
3.2.6由另一个中间时刻的角速度,计算另一个中间时刻的姿态变化四元数的增量△q3;
3.2.7将三个姿态变化四元数增量并进行加权平均,得到最终的姿态变化四元数增量:△q=(△q1+2×△q2+2×△q3+△q4)/6,其中△q4为使用当前时刻的角速度计算的姿态变化四元数的增量;
3.2.8将姿态变化四元数增量与上一时刻的姿态四元数相乘,得到当前时刻的姿态四元数:q(t)=q(t-1)×△q;
3.3将步骤2.8中得到的姿态变化四元数与上一时刻的四元数相乘,得到当前时刻的姿态四元数;
3.4在步骤(3),对姿态四元数进行归一化,使其范数等于1;
(4)姿态结果输出转换
根据姿态四元数,计算出钻孔的倾角和方位角,具体的
4.1通过姿态四元数转换为欧拉角,对倾斜角度进行表示;
4.2通过姿态四元数转换为旋转矩阵,姿态四元数q的方向向量表示为:D=[2(q1q3-q2q0),2(q2q3-q1q0),q0 2-q1 2-q2 2+q3 2],根据方向向量D计算方位角,采用反正切函数计算方向向量的x分量与y分量的比值,即atan2(D[1],D[0]);
(5)钻孔姿态估计问题建模
基于卡尔曼滤波器建立钻孔姿态估计问题建模,具体的:
5.1将钻孔姿态估计问题建立为一个状态空间模型,建立状态方程和观测方程;
5.2初始化卡尔曼滤波器的状态向量和协方差矩阵;
5.3根据钻孔当前姿态估计问题建模中的状态方程,使用过程模型对当前状态进行预测,并计算预测状态的协方差矩阵,预测步骤利用先验信息预测下一时刻的状态;
5.4根据系统建模中的观测方程,将各传感器测量值与预测状态进行比较,计算预测状态的测量残差,并计算卡尔曼增益;
5.5利用卡尔曼增益修正预测状态和协方差矩阵,得到更新后的状态估计和协方差;
(6)实际纠偏值计算
将实际测得的钻孔姿态通过向量运算与预设方向进行比较,确定偏差;
(7)纠偏校正值确认
根据计算得到的偏差值,采用PID算法进行校正;
(8)纠偏校正指令输出
根据PID算法的输出,生成对应控制指令,调整钻机的工作参数;
(9)纠偏后告警
持续监测钻孔的姿态数据,根据纠偏结果提供实时反馈。
进一步的设置为:在步骤(1)中,基于LoRa无线传输技术建立无线传输网络拓扑结构,以多个传感器节点组成网络,由中继点进行纠偏数据传输;
进一步的设置为:在步骤(2)中,所述预处理包括校准和滤波;
进一步的设置为:在步骤(3)中,所述初始化姿态解算将四元数进行初始化,并将其设为初始值[1,0,0,0]来表示钻杆非旋转状态即初始化状态;
进一步的设置为:在步骤4.1中,所述欧拉角为俯仰角、横滚角和航向角,其中俯仰角表示钻孔的倾斜角度;
进一步的设置为:在步骤5.1中,钻孔的真实姿态定义为状态变量,各传感器测量的姿态数据定义为观测变量,并定义状态变量和观测变量之间的关系;
进一步的设置为:在步骤(9)中,所述实施反馈具体为:当偏差超过预设阈值时,进行告警;当纠偏后需要二次纠偏时,进行告警。
工作原理
在钻杆上安装的倾角传感器,用于测量钻杆的倾角信息,方位角传感器用于测量钻杆的方位角信息,旋转速度传感器用于测量钻杆的旋转速度,在钻机臂上安装的倾角传感器和方位角传感器用于测量钻机结构的倾角信息和方位角信息,此处的中部、下部和上部以钻杆接触、伸入地面的一端为下进行定义,在纠偏数据采集中,对于数据采集后的记录也可以采用数据库建立将数据传输至数据库进行保存,再通过纠偏数据预处理消除传感器的误差和偏移,确保数据的准确性,其中滤波可以去除高频噪声和震荡,提高数据的稳定性。
通过四元姿态解算算法对纠偏姿态进行确定,在纠偏姿态计算中,传感器更新内的角速度数据用来表示当前钻杆的旋转速率,同时在姿态变化四元数的增量△q1时,对姿态四元数进行归一化,使其范数等于1,从而保持四元数的单位长度特性,同时通过基于卡尔曼滤波器建立的钻孔姿态估计问题建模避免出现姿态四元数的漂移和误差累积问题,其中卡尔曼增益将测量残差与预测步骤中的状态预测进行融合,得到最优估计的姿态,卡尔曼滤波不断地对状态进行预测和更新,将测量值和先验信息进行融合,有效抑制传感器数据中的漂移和误差累积,并提供更精确的姿态估计结果。
通过偏差计算将实际测得的钻孔姿态与预设方向进行比较,确定实际方向与预设方向之间的偏差,计算获取具体偏差值,采用PID算法进行校正,来调整钻机工作参数,从而实现自动校正,同时在调整后对钻杆状态进行监测,当偏差超过预设阈值时,进行告警,当纠偏后需要二次纠偏时,进行告警,从而及时提醒操作人员干预,避免出现钻卡事故,降低了操作人员在危险环境中的风险和受伤可能性。
区别传统深孔钻探,传感器数据的实时采集和处理使得钻孔偏差的检测和纠正可以在钻进过程中及时进行,同时采用四元姿态解算算法对纠偏姿态进行确定并通过卡尔曼滤波器建立的钻孔姿态估计问题建模避免出现姿态四元数的漂移和误差累积问题,实现自动化的钻孔偏差检测和纠正,减少人工干预,提高钻探作业效率,降低人工操作误差。
尽管参照前述实施例对本发明进行了详细的说明,对于本领域的技术人员来说,其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (7)
1.一种煤田钻探深孔钻探中的钻孔自动纠偏方法,其特征在于,包括以下步骤:
(1)纠偏数据采集
在钻杆中部安装倾角传感器,在钻杆下部安装方位角传感器,在钻杆上部安装旋转速度传感器,在钻机臂上同样安装倾角传感器和方位角传感器,多个传感器持续对钻杆钻探时的倾角、方位角、旋转速度进行监测,并进行记录;
(2)纠偏数据处理
对钻探时的倾角、方位角、旋转速度数据进行预处理;
(3)纠偏姿态计算
采用四元姿态解算算法对纠偏姿态进行确定,具体的:初始化姿态解算、传感器数据更新,其中传感器更新具体步骤为:
3.1从传感器数据中获取角速度数据;
3.2由当前时刻的获取的旋转速度数据采用四阶龙格库塔法来估计姿态变化四元数,具体的:
3.2.1将旋转速度传感器更新频率定义为△t,即离散化的时间步长;
3.2.2由当前时刻的角速度数据,计算姿态变化四元数的增量△q1;
3.2.3将当前时刻的角速度数据乘以△t的一半,得到中间时刻的角速度:
3.2.4由中间时刻的角速度,计算中间时刻的姿态变化四元数的增量△q2;
3.2.5再次将中间时刻的角速度乘以△t的一半,得到另一个中间时刻的角速度:
3.2.6由另一个中间时刻的角速度,计算另一个中间时刻的姿态变化四元数的增量△q3;
3.2.7将三个姿态变化四元数增量并进行加权平均,得到最终的姿态变化四元数增量:△q=(△q1+2×△q2+2×△q3+△q4)/6,其中△q4为使用当前时刻的角速度计算的姿态变化四元数的增量;
3.2.8将姿态变化四元数增量与上一时刻的姿态四元数相乘,得到当前时刻的姿态四元数:q(t)=q(t-1)×△q;
3.3将步骤2.8中得到的姿态变化四元数与上一时刻的四元数相乘,得到当前时刻的姿态四元数;
3.4在步骤(3),对姿态四元数进行归一化,使其范数等于1;
(4)姿态结果输出转换
根据姿态四元数,计算出钻孔的倾角和方位角,具体的:
4.1通过姿态四元数转换为欧拉角,对倾斜角度进行表示;
4.2通过姿态四元数转换为旋转矩阵,姿态四元数q的方向向量表示为:根据方向向量D计算方位角,采用反正切函数计算方向向量的x分量与y分量的比值,即atan2(D[1],D[0]);
(5)钻孔姿态估计问题建模
基于卡尔曼滤波器建立钻孔姿态估计问题建模,具体的:
5.1将钻孔姿态估计问题建立为一个状态空间模型,建立状态方程和观测方程;
5.2初始化卡尔曼滤波器的状态向量和协方差矩阵;
5.3根据钻孔当前姿态估计问题建模中的状态方程,使用过程模型对当前状态进行预测,并计算预测状态的协方差矩阵,预测步骤利用先验信息预测下一时刻的状态;
5.4根据系统建模中的观测方程,将各传感器测量值与预测状态进行比较,计算预测状态的测量残差,并计算卡尔曼增益;
5.5利用卡尔曼增益修正预测状态和协方差矩阵,得到更新后的状态估计和协方差;
(6)实际纠偏值计算
将实际测得的钻孔姿态通过向量运算与预设方向进行比较,确定偏差;
(7)纠偏校正值确认
根据计算得到的偏差值,采用PID算法进行校正;
(8)纠偏校正指令输出
根据PID算法的输出,生成对应控制指令,调整钻机的工作参数;
(9)纠偏后告警
持续监测钻孔的姿态数据,根据纠偏结果提供实时反馈。
2.根据权利要求1所述的一种煤田钻探深孔钻探中的钻孔自动纠偏方法,其特征在于:在步骤(1)中,基于LoRa无线传输技术建立无线传输网络拓扑结构,以多个传感器节点组成网络,由中继点进行纠偏数据传输。
3.根据权利要求1所述的一种煤田钻探深孔钻探中的钻孔自动纠偏方法,其特征在于:在步骤(2)中,所述预处理包括校准和滤波。
4.根据权利要求1所述的一种煤田钻探深孔钻探中的钻孔自动纠偏方法,其特征在于:在步骤(3)中,所述初始化姿态解算将四元数进行初始化,并将其设为初始值[1,0,0,0]来表示钻杆非旋转状态即初始化状态。
5.根据权利要求1所述的一种煤田钻探深孔钻探中的钻孔自动纠偏方法,其特征在于:在步骤4.1中,所述欧拉角为俯仰角、横滚角和航向角,其中俯仰角表示钻孔的倾斜角度。
6.根据权利要求1所述的一种煤田钻探深孔钻探中的钻孔自动纠偏方法,其特征在于:在步骤5.1中,钻孔的真实姿态定义为状态变量,各传感器测量的姿态数据定义为观测变量,并定义状态变量和观测变量之间的关系。
7.根据权利要求1所述的一种煤田钻探深孔钻探中的钻孔自动纠偏方法,其特征在于:在步骤(9)中,所述实施反馈具体为:当偏差超过预设阈值时,进行告警;当纠偏后需要二次纠偏时,进行告警。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310657569.4A CN116607929A (zh) | 2023-06-06 | 2023-06-06 | 一种煤田钻探深孔钻探中的钻孔自动纠偏方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310657569.4A CN116607929A (zh) | 2023-06-06 | 2023-06-06 | 一种煤田钻探深孔钻探中的钻孔自动纠偏方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN116607929A true CN116607929A (zh) | 2023-08-18 |
Family
ID=87678029
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310657569.4A Pending CN116607929A (zh) | 2023-06-06 | 2023-06-06 | 一种煤田钻探深孔钻探中的钻孔自动纠偏方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116607929A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116237818A (zh) * | 2022-12-29 | 2023-06-09 | 广东中海万泰技术有限公司 | 一种深孔加工的偏移量测量方法 |
-
2023
- 2023-06-06 CN CN202310657569.4A patent/CN116607929A/zh active Pending
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116237818A (zh) * | 2022-12-29 | 2023-06-09 | 广东中海万泰技术有限公司 | 一种深孔加工的偏移量测量方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11559149B2 (en) | Method and apparatus for transitioning between rotary drilling and slide drilling while maintaining a bit of a bottom hole assembly on a wellbore bottom | |
US11143011B2 (en) | Real-time modification of a slide drilling segment based on continuous downhole data | |
CN107407143B (zh) | 采用多个反馈回路的定向钻井方法和系统 | |
US10794169B2 (en) | Systems, devices, and methods for generating drilling windows | |
CN116607929A (zh) | 一种煤田钻探深孔钻探中的钻孔自动纠偏方法 | |
US11441411B2 (en) | Optimal drilling parameter machine learning system and methods | |
CN102892970A (zh) | 钻井方法和系统 | |
WO2009064656A2 (en) | Borehole survey method and apparatus | |
CN108643828B (zh) | 一种高坡锚固装备智能控制方法 | |
US11434742B2 (en) | Method and apparatus for identifying a potential problem with drilling equipment using a feedback control loop system | |
US11174718B2 (en) | Automatic steering instructions for directional motor drilling | |
CN111335879B (zh) | 一种钻孔轨迹测量装置 | |
US10760341B2 (en) | Automated steering of a drilling system using a smart bottom hole assembly | |
US10364666B2 (en) | Optimized directional drilling using MWD data | |
MX2011009048A (es) | Metodo para determinar la incertidumbre con la posicion y actitud del agujero del pozo proyectado. | |
WO2021066927A1 (en) | Probabilistic model calibration and estimation for drilling steering control | |
US11111770B2 (en) | Automated steering using operating constraints | |
CN111615582A (zh) | 用于钻井作业的方位角锁定的方法和系统 | |
US10934831B2 (en) | Downhole drilling using a network of drilling rigs | |
US11035971B2 (en) | Data quality monitoring and control systems and methods | |
US11408228B2 (en) | Methods and systems for improving confidence in automated steering guidance | |
US12031424B2 (en) | Methods and apparatus for optimizing downhole drilling conditions using a smart downhole system | |
US20230193741A1 (en) | Methods and apparatus for optimizing downhole drilling conditions using a smart downhole system | |
US10502043B2 (en) | Methods and devices to perform offset surveys | |
CN118166849A (zh) | 一种用于软土地层的沉井姿态实时监测系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination |