喬美英，王 波，肖學(xué)軍，許城寬，閆書豪

煤礦井下隨鉆測斜儀誤差聯(lián)合校正方法

喬美英1，王 波1，肖學(xué)軍2，許城寬1，閆書豪1

(1. 河南理工大學(xué) 電氣工程與自動化學(xué)院，河南 焦作 454000；2. 河北省煤田地質(zhì)局第四地質(zhì)隊，河北 張家口 075100)

為降低煤礦井下隨鉆測斜儀隨鉆測量鉆孔姿態(tài)角參數(shù)的誤差，建立基于遞推最小二乘的誤差補償數(shù)學(xué)模型，采用橢球擬合法、點積不變法和旋轉(zhuǎn)平面擬合法的聯(lián)合校正方法對隨鉆測斜儀進行誤差補償校正。對研發(fā)的隨鉆測斜儀進行數(shù)值仿真和雙軸轉(zhuǎn)臺實驗，用聯(lián)合校正方法對實驗數(shù)據(jù)進行誤差補償校正，校正前俯仰角和方位角的最大絕對誤差分別為4.7°和5.1°，校正后為0.8°和0.9°，校正后測量精度滿足隨鉆測斜儀要求。實驗結(jié)果表明，該方法可以對煤礦井下隨鉆測斜儀進行有效的誤差補償校正。

煤礦井下隨鉆測斜儀；聯(lián)合校正；誤差補償校正；俯仰角；方位角

煤礦井下鉆探技術(shù)在采集地層、構(gòu)造、水文地質(zhì)資料等地質(zhì)鉆孔以及探放水、抽采瓦斯等措施鉆孔中有著廣泛的應(yīng)用。為檢驗鉆孔軸跡，一定比例的鉆孔在鉆進時配備隨鉆測斜儀實時測量鉆孔的姿態(tài)角等參數(shù)，以掌握鉆孔軸跡，但是隨鉆測斜儀姿態(tài)測量模塊中的傳感器由于工藝及安裝技術(shù)等問題，導(dǎo)致測量姿態(tài)角參數(shù)存在誤差，無法精確地測繪鉆孔軸跡。

隨鉆測斜儀姿態(tài)測量模塊由三軸加速度計和三軸磁強計傳感器組合而成[1]，為了獲得精確的鉆孔測量參數(shù)，必須要對傳感器進行數(shù)據(jù)校正。研究傳感器校正方法較多，郝東等[2]采用赫姆霍茲線圈校正磁強計，代剛[3]采用離心機校正加速度計，這兩種方法的校正精度高，而且建模簡單，但是需要提供高精度的儀器，實驗環(huán)境也難以提供。郭衛(wèi)等[4]建立基于BP神經(jīng)網(wǎng)絡(luò)模型的方法校正加速度計；晁正正等[5]采用模擬退火、宋忠國等[6]采用粒子群優(yōu)化方法校正磁強計，通過訓(xùn)練優(yōu)化迭代得到校正需要的最優(yōu)擬合值，但訓(xùn)練需要大量數(shù)據(jù)而且難免出現(xiàn)泛化能力不足；陸欣等[7]采用最大似然校正法校正加速度計；龍禮等[8]采用自適應(yīng)校正方法校正磁強計，通過建立復(fù)雜數(shù)學(xué)模型可以校正加速度計或磁強計，但稍顯復(fù)雜。上述方法多難以應(yīng)用在隨鉆測斜儀工程實際中，并且多是對單一傳感器校正，隨鉆測斜儀姿態(tài)測量模塊需要考慮多傳感器之間校正對準問題、傳感器坐標系與載體坐標系的校正對準等問題。

在以上研究的基礎(chǔ)上，根據(jù)加速度計和磁強計誤差建立的測量模型，首先采用橢球擬合法校正加速度計，然后采用點積不變法校正磁強計以及傳感器之間的未對準誤差，最后采用旋轉(zhuǎn)平面法校正傳感器與載體坐標系的未對準誤差，并通過對研發(fā)的隨鉆測斜儀進行實驗來驗證該方法的有效性。

1 隨鉆測斜儀隨鉆測量原理

煤礦井下隨鉆測斜儀姿態(tài)測量模塊由三軸加速度計和三軸磁強計組成，其原理如圖1所示。

圖1 姿態(tài)測量原理示意圖

在鉆孔作業(yè)時，加速度計實時測量地球的重力場分量G、G、G，磁強計實時測量地球的磁場分量B、B、B。測量值經(jīng)過數(shù)值轉(zhuǎn)換根據(jù)式(1)計算出俯仰角，工具面角和方位角，可以確定鉆孔的實時姿態(tài)[9-10](其中0為重力矢量和)。

2 加速度計的單獨校正

隨鉆測斜儀中的磁強計容易受到周圍環(huán)境磁干擾，測量數(shù)據(jù)不精確，校正后的結(jié)果對第二步的點積不變法產(chǎn)生影響，因此，第一步選擇單獨校正加速度計。

2.1 加速度計誤差校正模型

加速度計由于工藝及安裝技術(shù)等的影響存在各種誤差[11]。根據(jù)文獻[4,7,12]得到加速度計測量模型：

式中：c為加速度計測量值；z為加速度計真實值；b為比例系數(shù)誤差；f為非正交誤差；w為未對準誤差；0為零位偏差。

引入可逆矩陣=bfw和=–1并代入式(2)，整理得到加速度計誤差補償校正模型：

2.2 誤差補償校正與橢球擬合模型的轉(zhuǎn)換

理想情況下，加速度計以任意姿態(tài)旋轉(zhuǎn)，穩(wěn)定下加速度計測量重力矢量不變，輸出軌跡是一個正球體，但是受誤差影響軌跡會畸變?yōu)闄E球。

對加速度真實值z求模的平方：

式(4)整理為橢球標準式：

橢球的一般式為：

式(6)整理為橢球矩陣式：

根據(jù)式(7)可以得到：

—，—是橢球的系數(shù)，其中0表示球心。根據(jù)式(5)、式(7)、式(8)得到：

將式(10)代入式(3)，加速度計誤差補償校正模型可以表示為：

加速度計校正，關(guān)鍵是擬合出橢球系數(shù)。

2.3 基于遞推最小二乘的橢球系數(shù)

遞推最小二乘算法是一種在線算法，利用原有參數(shù)進行下一步估計，可以做到運算量小，實時進行估計，逐步得到最優(yōu)解，而且遞推過程不需要批量處理，隨著迭代次數(shù)的增加估計精度也越來越高[13-14]。

根據(jù)2.2節(jié)可知加速度計校正模型可以用橢球模型來表示。

帶有噪聲的橢球方程為：

式中：()為噪聲；、、為加速度計測量分量。橢球系數(shù)定義為：

加速度計的實時測量數(shù)據(jù)為[,,]T，按照橢球參數(shù)可定義：

由式(12)—式(14)得到遞推最小二乘式：

根據(jù)式(16)遞推擬合橢球系數(shù)，其中，各項的初值設(shè)為0，初值設(shè)為106(為單位矩陣)，設(shè)為常值1。

3 加速度計和磁強計之間的相互校正

采用點積不變法校正磁強計，校正過程中加速度計和磁強計之間的未對準誤差也可以消除[15-16]。點積不變法是利用了相對位置不變的兩條矢量數(shù)量積為常數(shù)的性質(zhì)。假設(shè)地磁矢量是大小和方向都不變的常數(shù)矢量，重力矢量與地磁矢量在任意姿態(tài)下的點積都為常數(shù)(const)。設(shè)重力矢量為z，地磁矢量為z，矢量公式為：

以上研究為基層國家機關(guān)“誰執(zhí)法誰普法”責(zé)任制的建立和實施提供了理論依據(jù)，對法律在全社會的普及大有裨益。

根據(jù)文獻[5-6,8,16]得到磁強計測量模型：

式中：c為磁強計測量值；z為磁強計真實值；is為比例系數(shù)誤差；n為非正交誤差；w為未對準誤差；s為軟鐵誤差；h為硬鐵誤差；為零位偏差。引入=isnws、0=isnwh+(0為總偏差)和m=–1代入式(18)整理得：

將z和式(19)代入矢量公式(17)得：

定義m=，m0=代入式(20)得：

式中：zT為經(jīng)過校正的加速度計真實值；c為磁強計測量值。多次測量建立超定方程組，利用遞推最小二乘擬合得到和的值，代入m=，m0=計算得到磁強計的誤差補償矩陣m和總偏差0。

4 傳感器與載體坐標系的未對準誤差校正

橢球擬合法不能校正傳感器坐標系與載體坐標系之間的未對準誤差[17]。這里引入坐標系旋轉(zhuǎn)理論，未對準誤差可以看作傳感器坐標系()與載體坐標系()之間的多次旋轉(zhuǎn)，假設(shè)兩個坐標系之間的未對準角度為、、，按、、順序旋轉(zhuǎn)，旋轉(zhuǎn)矩陣由式(22)來表示[11]：

旋轉(zhuǎn)平面擬合，是指在均勻的矢量場中，在固定的位置下，傳感器繞任意軸旋轉(zhuǎn)時，在旋轉(zhuǎn)軸輸出是一個常數(shù)，另外兩軸的輸出矢量將會形成一個圓環(huán)平面[18]。假設(shè)A、A、A是加速度計經(jīng)過橢球擬合校正后的值，同時引入：

當繞軸轉(zhuǎn)動，在載體坐標軸輸出為常數(shù)。

式(24)、式(25)為兩個空間平面方程，按照++=1代入橢球擬合校正后的值建立超定方程組，求出、、和、、，代入式(23)計算出未對準小角度；代入式(22)計算出結(jié)果并用w表示，代入式(3)、式(19)求出加速度計誤差補償矩陣w–1，磁強計誤差補償矩陣w–1m。

5 實驗驗證

實驗分為數(shù)值仿真和雙軸轉(zhuǎn)臺(實測數(shù)據(jù))實驗。數(shù)值仿真是為了驗證方法的準確性，雙軸轉(zhuǎn)臺用于驗證方法實際有效性。

5.1 數(shù)值仿真

根據(jù)式(2)加速度計的測量模型做出如下假設(shè)：比例誤差設(shè)為：b=diag(1.2,1,0.8)；偏差設(shè)為：0= [0.1,0.2,0.3]T；三軸非正交誤差角設(shè)為=1°，=2°，=3°；坐標系未對準誤差角設(shè)為=2°，=1°，=1.5°；三軸分量設(shè)為[0,0,1]T[7]。根據(jù)設(shè)置的誤差參數(shù)計算誤差補償矩陣和偏差的參考值。

根據(jù)式(18)磁強計的測量模型做出如下假設(shè)：比例誤差設(shè)為：is=diag(1.2,0.8,1)；偏差設(shè)為：=[900,860,–780]T；三軸非正交誤差角設(shè)為=1°，=0.5°，=1°；坐標系未對準誤差角設(shè)為=1°，=2°，=3°；硬鐵誤差設(shè)為h=[–75,–70,85]T；軟鐵誤差設(shè)為s=[1.2,0.2,0.3;0.3,1.4,0.4;0.2,0.4,1.5]；三軸分量設(shè)為[0,0,48 152]T[16]，根據(jù)設(shè)置的誤差參數(shù)計算誤差補償矩陣和偏差的參考值。

圖2 橢球擬合參數(shù)估計

假設(shè)載體坐標系三軸為東北天指向即軸指向東，軸指向北，軸指向天。根據(jù)坐標系旋轉(zhuǎn)理論做旋轉(zhuǎn)運動，旋轉(zhuǎn)角度分別設(shè)為繞軸、軸、軸旋轉(zhuǎn)360°，繞軸旋轉(zhuǎn)每次間隔120°，繞軸旋轉(zhuǎn)每次間隔60°，繞軸旋轉(zhuǎn)每次間隔30°，共計216個測點。根據(jù)文獻[11]中的旋轉(zhuǎn)矩陣公式求出每次旋轉(zhuǎn)后的加速度計和磁強計的三軸分量，作為加速度計和磁強計無誤差輸出數(shù)據(jù)，按照上文求出的參考值計算加速度計和磁強計誤差輸出數(shù)據(jù)。

將加速度計誤差數(shù)據(jù)代入MATLAB軟件編寫的遞推最小二乘程序，其中設(shè)置為106，橢球系數(shù)—、—初始值都設(shè)置為0，經(jīng)過24次迭代，參數(shù)趨于穩(wěn)定，迭代輸出結(jié)果如圖2所示，橢球系數(shù)=[0.846 2，1.220 1，1.908 0，–0.016 4，–0.040 9，–0.083 4，–0.069 5，–0.217 3，–0.551 8]T。

根據(jù)遞推的結(jié)果計算經(jīng)橢球擬合的加速度計誤差補償矩陣和偏差；計算經(jīng)橢球+旋轉(zhuǎn)平面擬合的加速度計誤差補償矩陣和偏差，結(jié)果見表1。

表1 加速度計參考值以及各校正方法結(jié)果

將磁強計誤差數(shù)據(jù)按照點積不變法計算磁強計誤差補償矩陣和偏差；計算點積不變法+旋轉(zhuǎn)平面擬合的磁強計誤差補償矩陣和偏差，結(jié)果見表2。

表2 磁強計參考值以及各校正方法結(jié)果

MATLAB軟件繪制加速度計校正前和經(jīng)橢球擬合校正后測點空間分布圖以及磁強計校正前和經(jīng)點積不變法校正后測點空間分布圖，如圖3所示，兩者校正前都表現(xiàn)為橢圓，并且分布在球體外，校正后表現(xiàn)為圓形分布在球體表面。

橢球擬合校正后的加速度計與點積校正后磁強計的未對準誤差消除，即傳感器坐標系對準。將磁強計數(shù)值縮小48 152倍(根據(jù)設(shè)的分量大小縮小，只改變大小，不改變方向)，目的是將兩個傳感器的坐標系繪制在同一個坐標系。繪制橢球擬合校正后的加速度計與點積法校正后的磁強計測點俯視圖，如圖4a所示。

圖3 校正前后測點空間分布

經(jīng)旋轉(zhuǎn)平面擬合的傳感器坐標系和載體坐標系未對準誤差消除。數(shù)值仿真中表現(xiàn)為傳感器所在坐標系和載體坐標系坐標軸方向重合，如圖4b所示。

5.2 雙軸轉(zhuǎn)臺實驗

為進一步說明方法的有效性，采用自主研發(fā)的隨鉆測斜儀進行實驗。該隨鉆測斜儀主要由單片機(STC15F2k 60S2)、姿態(tài)測量單元(MMA8451加速度計和HMC5883L磁強計)、存儲器(16兆字節(jié)Flash芯片)、時鐘(DS1302)和通信接口構(gòu)成。該主控芯片單片機速度是傳統(tǒng)單片機的8~12倍，對傳感器數(shù)據(jù)的采集和處理會更加高效與準確，實驗過程中姿態(tài)測量模塊測量的數(shù)據(jù)，通過UART協(xié)議傳輸?shù)絾纹瑱C，并存儲到Flash中。

5.2.1 擬合校正參數(shù)

將隨鉆測斜儀用固定板固定在雙軸轉(zhuǎn)臺上，分別繞、、軸以任意姿態(tài)轉(zhuǎn)動若干次，為避免非重力加速度干擾，獲取轉(zhuǎn)動結(jié)束時靜止狀態(tài)下數(shù)據(jù)，對每個姿態(tài)測量10次，并求均值，共獲得100組測量數(shù)據(jù)(加速度計和磁強計讀數(shù))。

表3是利用本文方法求出加速度計和磁強計的校正參數(shù)，計算時需要進行數(shù)值轉(zhuǎn)換。

圖4 未對準誤差校正

表3 加速度計和磁強計校正參數(shù)

5.2.2 姿態(tài)角補償校正

補償校正最終目的是降低隨鉆測斜儀測量的姿態(tài)角誤差，即俯仰角和方位角誤差。

a.俯仰角 在雙軸轉(zhuǎn)臺上固定好自主研發(fā)的隨鉆測斜儀。以水平0°為基準，控制轉(zhuǎn)臺以間隔10°順序上仰到80°，接著回到水平基準再控制轉(zhuǎn)臺以間隔10°順序下俯到–80°，為避免非重力加速度干擾，每次轉(zhuǎn)動結(jié)束時在靜止狀態(tài)下采集測量數(shù)據(jù)并存儲在Flash中，每個角度測量10次求均值，最終獲得17個測量值。對采集數(shù)據(jù)按照表3中加速度計參數(shù)校正，繪制校正前與校正后的示意圖，如圖5所示，校正前分布在球體外，表現(xiàn)為橢球，校正后分布在球體表面。

圖5 加速度計校正前后示意圖

根據(jù)式(1)計算校正前和校正后的俯仰角，并分別計算校正前和校正后與參考值(–80∶0∶80)的差。繪制校正前和校正后的俯仰角誤差曲線，如圖6所示。上仰和下俯的誤差最大絕對值分別為3.9°和4.7°，校正后為0.7°和0.8°，校正前后，測量精度提高了一個等級。

圖6 俯仰角誤差曲線

b.方位角 在同樣的實驗環(huán)境下以水平0°為基準，精密羅盤確定地理北向為0°，將隨鉆測斜儀指北向0°固定。控制轉(zhuǎn)臺以0°為起始隔10°測量，順時針旋轉(zhuǎn)360°，測量數(shù)據(jù)存儲在Flash中，每個方位角測量10次，求其均值，最終獲得36個測量點。對采集數(shù)據(jù)按照表3中磁強計參數(shù)校正，繪制校正前與校正后的示意圖，如圖7所示，校正前分布在球體外，表現(xiàn)為橢球，偏離球心，校正后分布在球體表面，指向球心。

根據(jù)式(1)計算校正前和校正后的方位角，并分別計算校正前和校正后與參考值(0∶10∶350)的差。繪制校正前和校正后的方位角誤差曲線如圖8所示。校正前，從0°轉(zhuǎn)至180°和180°轉(zhuǎn)至350°之間誤差最大絕對值分別為5.1°和4.9°，校正后分別為0.8°和0.9°，校正前后，測量精度提高了一個等級。

圖7 磁強計校正前后示意圖

圖8 方位角誤差曲線

6 結(jié)論

a. 對煤礦井下隨鉆測斜儀測量模塊中的加速度計和磁強計進行誤差源分析并基于遞推最小二乘，采用橢球擬合法、點積不變法和旋轉(zhuǎn)平面擬合法的聯(lián)合校正方法完成對隨鉆測斜儀的校正，俯仰角和方位角的測量精度提高了一個數(shù)量級。

b. 聯(lián)合校正方法是一種低成本，計算方法簡單，精度可靠，方便應(yīng)用的校正，適用于煤礦井下隨鉆測斜儀校正。

請聽作者語音介紹創(chuàng)新技術(shù)成果等信息，歡迎與作者進行交流

[1] GIZAWY E L. Wellbore surveying while drilling based on Kalman filtering[J]. American Journal of Engineering and Applied Sciences，2010，3(2)：240–259.

[2] 郝東，繩濤，陳小前. 三軸磁強計測量誤差修正方法[J]. 航天器環(huán)境工程，2011，28(5)：463–466. HAO Dong，SHENG Tao，CHEN Xiaoqian. The error correction of three-axis magnetometer measurement[J]. Spacecraft Environment Engineering，2011，28(5)：463–466.

[3] 代剛. MEMS–IMU誤差分析補償與實驗研究[D]. 北京：清華大學(xué)，2011. DAI Gang. MEMS-IMU error analysis compensation method and experiment research[D]. Beijing：Tsinghua University，2011.

[4] 郭衛(wèi)，楊鵬飛，張武剛，等. 基于BP神經(jīng)網(wǎng)絡(luò)模型的MEMS加速度計誤差補償方法[J]. 中國測試，2018，44(3)：109–113. GUO Wei，YANG Pengfei，ZHANG Wugang，et al. Error compensation method for MEMS accelerometer based on BP neural network model[J]. China Measurement & Test，2018，44(3)：109–113.

[5] 晁正正，張曉明，馬喜宏，等. 一種新的電子羅盤校準算法研究[J]. 傳感技術(shù)學(xué)報，2019，32(1)：106–110.CHAO Zhengzheng，ZHANG Xiaoming，MA Xihong，et al. Research on a new electronic compass calibration algorithm[J]. Chinese Journal of Sensors and Actuators，2019，32(1)：106–110.

[6] 宋忠國，鄭家歡，張金生，等. 基于粒子群優(yōu)化的三軸磁強計非線性誤差校正[J]. 傳感器與微系統(tǒng)，2017，36(12)：40–42. SONG Zhongguo，ZHENG Jiahuan，ZHANG Jinsheng，et al. Nonlinear error calibration of three-axis magnetometer based on PSO algorithm[J]. Transducer and Microsystem Technologies，2017，36(12)：40–42.

[7] 陸欣，劉忠，張宏欣，等. 三軸MEMS加速度計的最大似然校正算法[J]. 國防科技大學(xué)學(xué)報，2017，39(5)：185–191. LU Xin，LIU Zhong，ZHANG Hongxin，et al. Maximum likelihood calibration for MEMS triaxial accelerometer[J]. Journal of National University of Defense Technology，2017，39(5)：185–191.

[8] 龍禮，張合. 三軸地磁傳感器誤差的自適應(yīng)校正方法[J]. 儀器儀表學(xué)報，2013，34(1)：161–165. LONG Li，ZHANG He. Automatic and adaptive calibration method of tri-axial magnetometer[J]. Chinese Journal of Scientific Instrument，2013，34(1)：161–165.

[9] 江澤宇，謝洪波，文廣超，等. 煤礦井下電磁波無線隨鉆軌跡測量系統(tǒng)設(shè)計與應(yīng)用[J]. 煤田地質(zhì)與勘探，2017，45(3)：156–161. JIANG Zeyu，XIE Hongbo，WEN Guangchao，et al. Design and application of electromagnetic radio MWD system of drilling track in coal mine[J]. Coal Geology & Exploration，2017，45(3)：156–161.

[10] 王成立，CHIKHOTKIN V F，盧春華. 隨鉆測斜儀誤差補償實驗研究[J]. 煤田地質(zhì)與勘探，2016，44(4)：147–152. WANG Chengli，CHIKHOTKIN V F，LU Chunhua. Error compensation of measurement while drilling inclinometer[J]. Coal Geology & Exploration，2016，44(4)：147–152.

[11] 鄧志紅. 慣性器件與慣性導(dǎo)航系統(tǒng)[M]. 北京：科學(xué)出版社，2012. DENG Zhihong. Inert components and inert navigation system[M]. Beijing：Science Press，2012.

[12] 劉艷霞，方建軍，楊清梅. 基于橢球假設(shè)的三軸加速度計誤差標定與補償[J]. 傳感器與微系統(tǒng)，2014，33(6)：52–54. LIU Yanxia，F(xiàn)ANG Jianjun，YANG Qingmei. Three-axis accelerometer error calibration and compensation based on ellipsoid hypothesis[J]. Transducer and Microsystem Technologies，2014，33(6)：52–54.

[13] 郭鵬飛，華春紅，任章，等. 基于遞推最小二乘的航姿系統(tǒng)羅差校正[J]. 中國慣性技術(shù)學(xué)報，2008，16(1)：24–27. GUO Pengfei，HUA Chunhong，REN Zhang，et al. Magneticdeviation compensation using recursive least square for AHRS[J]. Journal of Chinese Inertial Technology，2008，16(1)：24–27.

[14] 龍禮，黃家才. 基于遞推最小二乘法的地磁測量誤差校正方法[J]. 儀器儀表學(xué)報，2017，38(6)：1440–1446. LONG Li，HUANG Jiacai. Recursive least square based online error calibration method in geomagnetic detection[J]. Chinese Journal of Scientific Instrument，2017，38(6)：1440–1446.

[15] 李翔，李智. 航姿參考系統(tǒng)三軸磁強計校正的點積不變法[J]. 儀器儀表學(xué)報，2012，33(8)：1813–1818. LI Xiang，LI Zhi. Dot product invariance method for the calibration of three-axis magnetometer in attitude and heading reference system[J]. Chinese Journal of Scientific Instrument，2012，33(8)：1813–1818.

[16] 李翔，王勇軍，李智. 航姿系統(tǒng)矢量傳感器非對準誤差及其校正[J]. 傳感技術(shù)學(xué)報，2017，30(2)：266–271. LI Xiang，WANG Yongjun，LI Zhi. Inter-triad misalignment of vector field sensors in attitude andheading reference systems and its calibration[J]. Chinese Journal of Sensors and Actuators，2017，30(2)：266–271.

[17] 李智，李翔. 基于橢球假設(shè)的三軸電子羅盤羅差補償研究[J]. 儀器儀表學(xué)報，2011，32(10)：2210–2215. LI Zhi，LI Xiang. Research on magnetic deviation compensation of three-axis electronic compass based on ellipsoid hypothesis[J]. Chinese Journal of Scientific Instrument，2011，32(10)：2210–2215.

[18] 梁偉嘉，吝伶艷，劉宗偉. 基于橢球和平面擬合的隨鉆測量系統(tǒng)校正方法[J]. 煤炭技術(shù)，2016，35(9)：250–251. LIANG Weijia，LIN Lingyan，LIU Zongwei. Calibration method of measurement while drilling system based on ellipsoid and plane fitting[J]. Coal Technology，2016，35(9)：250–251.

Joint correction method of errors of MWD inclinometer in underground coal mine

QIAO Meiying1, WANG Bo1, XIAO Xuejun2, XU Chengkuan1, YAN Shuhao1

(1. School of Electrical Engineering and Automation, Henan Polytechnic University, Jiaozuo 454000, China; 2. Fourth Geological Team, Hebei Provincial Bureau of Coal Geology, Zhangjiakou 075100, China)

In order to reduce the error of measuring the drilling attitude angle parameter of MWD inclinometer in coal mine, a mathematical model of error compensation based on recursive least squares was established, the joint correction method of ellipsoid fitting method, dot product invariance method and rotation plane fitting method were used in the error compensation correction for MWD inclinometer. Numerical simulation experiment and double-axis turntable experiment were carried out on the developed MWD inclinometer, error compensation correction of experimental data was conducted with joint correction method, the maximum absolute errors of pitch angle and azimuth angle before correction are 4.7° and 5.1°, after correction are 0.8° and 0.9°, the corrected measurement accuracy meets the requirements of MWD inclinometer. Experimental results show that this method can be used to effectively compensate the error correction of MWD inclinometer in underground coal mine.

MWDinclinometer in underground coal mine; joint correction; the error compensation correction; pitch angle; azimuth angle

TE243

A

10.3969/j.issn.1001-1986.2020.02.030

1001-1986(2020)02-0202-07

2019-07-16；

2019-08-22

國家自然科學(xué)基金項目(61573129)；河南省教育廳重點科研項目(16A120004)

National Natural Science Foundation of China(61573129)；Key Research Project of Henan Provincial Educational Department(16A120004)

喬美英，1976年生，女，內(nèi)蒙古烏拉特前旗人，博士，副教授，從事機器學(xué)習(xí)及時間序列數(shù)據(jù)分析的研究. E-mail：qiaomy@hpu.edu.cn

喬美英，王波，肖學(xué)軍，等. 煤礦井下隨鉆測斜儀誤差聯(lián)合校正方法[J]. 煤田地質(zhì)與勘探，2020，48(2)：202–208.

QIAO Meiying，WANG Bo，XIAO Xuejun，et al. Joint correction method of errors of MWD inclinometer in underground coal mine[J]. Coal Geology & Exploration，2020，48(2)：202–208.

(責(zé)任編輯 聶愛蘭)