徐凱帆 張書畢 鮑 國(guó)

(1.廣東省重工建筑設(shè)計(jì)院有限公司，廣東 廣州510034;2.江蘇省資源環(huán)境信息工程重點(diǎn)實(shí)驗(yàn)室，江蘇徐州221116;3.中國(guó)礦業(yè)大學(xué)環(huán)境與測(cè)繪學(xué)院，江蘇徐州221116;4.中國(guó)人民解放軍空軍勤務(wù)學(xué)院機(jī)場(chǎng)工程與保障系，江蘇徐州221008)

井下控制網(wǎng)測(cè)量不僅為井下提供準(zhǔn)確的空間基準(zhǔn)，也為礦山安全生產(chǎn)提供所需的各項(xiàng)數(shù)據(jù)，在礦業(yè)開發(fā)過程中發(fā)揮著重要作用［1-8］。由于受井下巷道條件的限制，控制網(wǎng)一般以導(dǎo)線的形式沿巷道分期布設(shè)［1，9］。隨著巷道的掘進(jìn)，需要加測(cè)一定數(shù)量的陀螺定向邊來提高井下控制網(wǎng)的精度。在控制網(wǎng)內(nèi)業(yè)計(jì)算時(shí)，傳統(tǒng)方法是對(duì)各期觀測(cè)數(shù)據(jù)進(jìn)行獨(dú)立的平差計(jì)算［1］，但由于受到系統(tǒng)誤差的影響，導(dǎo)致控制網(wǎng)局部誤差累積過大。此外，每期控制網(wǎng)平差計(jì)算的起算點(diǎn)位置不同，在一定程度上影響了點(diǎn)位精度的評(píng)定。若采用傳統(tǒng)方法對(duì)各期的觀測(cè)值進(jìn)行整體平差，可能會(huì)由于前幾觀測(cè)數(shù)據(jù)的部分丟失給整體平差帶來困難。為此，利用前期控制網(wǎng)點(diǎn)的平差值和方差、協(xié)方差矩陣，結(jié)合新的觀測(cè)值推導(dǎo)變參數(shù)序貫平差模型，并利用無限權(quán)理論將陀螺方位角條件轉(zhuǎn)化為變參數(shù)序貫平差模型的觀測(cè)方程。

1 理論分析

1.1 變參數(shù)序貫平差

變參數(shù)序貫平差是指在平差過程中的未知參數(shù)的數(shù)量有增有減［6-8］，主要包括3 種情況:①第2 次平差增加新的參數(shù);②第2 次平差的參數(shù)僅為第1 次平差參數(shù)的一部分;③上述2 種情況的綜合。為了便于推導(dǎo)，給出井下控制網(wǎng)2 期觀測(cè)值的誤差方程

式中，V1、V2分別為2 期觀測(cè)值的改正數(shù)矩陣; 為第1 期出現(xiàn)的參數(shù)；為2 期共有的參數(shù);為第2期新增加的參數(shù);A、B1、B2和C 為系數(shù)矩陣;l1和l2為常數(shù)向量。

對(duì)式(1)進(jìn)行第1 次平差，其法方程為［10-12］

式中，P1為第1 期觀測(cè)值的權(quán)陣;和為第1 次平差參數(shù)的改正數(shù)。

Ea、Eb為單位矩陣，分別為參數(shù)改正數(shù)的系數(shù)矩陣。令

根據(jù)最小二乘原理，令

將式(4)代入式(5)，整理后得到法方程為

式中，P1和P2分別為2 期觀測(cè)值的權(quán)陣。

由(6)式可得參數(shù)平差值為

精度評(píng)定時(shí)，單位權(quán)中誤差估值可表示為［13-17］

式中，t 為獨(dú)立參數(shù)的個(gè)數(shù);n 為觀測(cè)值個(gè)數(shù);P 為2期觀測(cè)值的權(quán)陣;VTPV 可進(jìn)一步地表示為

1.2 附加約束條件

井下布設(shè)基本控制導(dǎo)線時(shí)，一般每間隔1. 5 ～2.0 km 加測(cè)1 條陀螺定向邊。對(duì)于已經(jīng)建立井下控制網(wǎng)的礦井，為了提高井下控制網(wǎng)的平面精度，有時(shí)需要加測(cè)一定數(shù)量的陀螺定向邊(見圖1)。在數(shù)據(jù)處理時(shí)，通常將陀螺定向邊作為堅(jiān)強(qiáng)邊處理［1］，即將實(shí)測(cè)的陀螺定向邊方位角作為約束條件參與控制網(wǎng)的平差計(jì)算。

圖1 陀螺定向邊Fig.1 Gyro directional edge

圖1 中，gf 為一條陀螺定向邊，當(dāng)以陀螺定向邊方位角為約束條件時(shí)，觀測(cè)值為真值應(yīng)滿足

式中，Lgf為陀螺定向邊方位角觀測(cè)值，(°);αgf為坐標(biāo)平差值反算的坐標(biāo)方位角，(°)。

將式(10)線性化后可得待求參數(shù)的限制條件方程，該方程與式(1)、式(2)為2 種不同類型的方程，不便于利用變參數(shù)序貫平差模型進(jìn)行計(jì)算。為此，根據(jù)無限權(quán)理論［6］，將陀螺定向邊方位角作為一般觀測(cè)值處理。此時(shí)，該觀測(cè)值對(duì)應(yīng)的權(quán)值為

式(11)表明陀螺方位角值的權(quán)重?zé)o窮大，即對(duì)應(yīng)觀測(cè)值的改正數(shù)為0。

由近似坐標(biāo)改正數(shù)引起的近似坐標(biāo)方位角的改正數(shù)為δαfg，即

式中，Vgyr為陀螺方位角改正數(shù);Lgyr為陀螺方位角觀測(cè)值。

將式(13)整理后，得

式(14)為陀螺定向邊的觀測(cè)值誤差方程，可與水平角、邊長(zhǎng)觀測(cè)值的誤差方程統(tǒng)一表示為式(1)、式(2)的形式，可利用變參數(shù)序貫平差模型進(jìn)行求解。

1.3 權(quán)陣的確定

2 實(shí)例分析

2.1 井下控制網(wǎng)概況

井下控制網(wǎng)往往隨井下巷道的道掘進(jìn)而逐步布設(shè)。圖2 為某礦井7″級(jí)局部導(dǎo)線網(wǎng)，共由119 個(gè)點(diǎn)組成，導(dǎo)線全長(zhǎng)約12 km，其中JB1和JB2為已知點(diǎn)。第1 期形成的網(wǎng)形僅為一條與JB1點(diǎn)相連的閉合導(dǎo)線。虛線框中的導(dǎo)線及WE1－WE2、YS7－YS8和XF22－XF23為第2 期施測(cè)的3 條陀螺定向邊，其中陀螺定向邊被視為堅(jiān)強(qiáng)邊，且3 條陀螺定向邊距起算點(diǎn)JB1分別為1.5，3，6 km。現(xiàn)需要實(shí)現(xiàn)2 期控制網(wǎng)的整體平差。

2.2 平差結(jié)果分析

2.2.1 定權(quán)系數(shù)確定

陀螺定向邊作為堅(jiān)強(qiáng)邊處理時(shí)，雖然其權(quán)值Pg→+ ∞，但在實(shí)際應(yīng)用中為了避免法方程病態(tài)［18-20］，陀螺定向邊觀測(cè)值的權(quán)值(Pg)往往選取一個(gè)有限值，根據(jù)觀測(cè)值權(quán)的定義

圖3 不同定權(quán)系數(shù)對(duì)應(yīng)的陀螺定向邊方位角差值Fig.3 Azimuth variation of Gyro directional edges with different weight coefficient

由圖3 可知，當(dāng)定權(quán)系數(shù)m = 1 時(shí)，最大的差值達(dá)到5.1″，但隨著m 值增大，陀螺邊方位角的差值迅速減小，當(dāng)m = 3 時(shí)，各陀螺邊方位角的平差值與已知值之差均小于1″;當(dāng)m = 10 時(shí)，其差值均小于0.1″;當(dāng)m = 20 時(shí)，各陀螺定向邊的平差值與已知值之差均小于0.01″。

為了全面分析定權(quán)系數(shù)的取值對(duì)平差結(jié)果的影響，表1 給出了m =1 ～20 時(shí)，陀螺定向邊各端點(diǎn)坐標(biāo)的變化量。

表1 陀螺定向邊各端點(diǎn)坐標(biāo)變化量Table 1 Coordinate variation of each gyro directional edge endpoint

由表1 可知，陀螺定向邊距起算點(diǎn)越近，其端點(diǎn)坐標(biāo)的變化量則越大。

2.2.2 平差結(jié)果驗(yàn)證

井下控制網(wǎng)平差計(jì)算分別采用2 種方法進(jìn)行:①利用控制網(wǎng)2 期完整的觀測(cè)數(shù)據(jù)，采用附有限制條件的間接平差方法進(jìn)行整體平差，該結(jié)果視為真值;②根據(jù)控制網(wǎng)第1 期平差結(jié)果，采用附加約束條件的變參數(shù)序貫平差模型進(jìn)行控制網(wǎng)整體平差，其中，定權(quán)系數(shù)m 取20。為了全面比較2 種方法平差結(jié)果的差異，圖4 給出了控制網(wǎng)各點(diǎn)的坐標(biāo)差。

圖4 各控制點(diǎn)坐標(biāo)殘差值Fig.4 Coordinate residuals of each control point

由圖4可知，坐標(biāo)差值均小于0. 03mm，超出了結(jié)果需要保留的有效位數(shù)，其差值可以忽略不計(jì)。也就是說，如果坐標(biāo)殘差值以毫米為單位，結(jié)果保留至0.1 mm，單位權(quán)中誤差以秒為單位，保留至0.1″時(shí)，則2 種方法平差結(jié)果的精度完全一致。但文中提出的平差模型無需存儲(chǔ)控制網(wǎng)的前期觀測(cè)值，即使前期數(shù)據(jù)丟失，仍可根據(jù)前期的平差結(jié)果進(jìn)行計(jì)算，以達(dá)到整體平差的效果。

3 結(jié) 論

(1)變參數(shù)序貫平差模型僅需前期計(jì)算結(jié)果，便可實(shí)現(xiàn)井下控制網(wǎng)多期觀測(cè)整體平差，可有效克服控制網(wǎng)前期觀測(cè)數(shù)據(jù)丟失給整體平差帶來的困難。

(2)在井下控制網(wǎng)整體平差時(shí)，可利用無限權(quán)理論將陀螺堅(jiān)強(qiáng)邊轉(zhuǎn)化為變參數(shù)序貫平差模型的觀測(cè)方程。實(shí)例分析結(jié)果表明，陀螺邊的定權(quán)系數(shù)m 取20 時(shí)，可滿足井下控制網(wǎng)計(jì)算的精度要求。

(3)附加限制條件的變參數(shù)序貫平差模型思路簡(jiǎn)單、易于編程、無需前期觀測(cè)數(shù)據(jù)，適用于井下分期布設(shè)的控制網(wǎng)整體平差。

［1］ 張國(guó)良，朱家鈺，顧和和.礦山測(cè)量學(xué)［M］. 徐州:中國(guó)礦業(yè)大學(xué)出版社，2008.

Zhang Guoliang，Zhu Jiayu，Gu Hehe.Mine Surveying［M］.Xuzhou:China University of Mining and Technology Press，2008.

［2］ 崔 清，陳 曦，賈旭光，等. 基于測(cè)距技術(shù)的邊坡角測(cè)定誤差及敏感性分析［J］.金屬礦山，2014(5):121-124.

Cui Qing，Chen Xi，Jia Xuguang，et al. Slope angle measuring error and sensitivity analysis based on the distance measurement technology［J］.Metal Mine，2014(5):121-124.

［3］ 賈寶新，李國(guó)臻. 礦山地震監(jiān)測(cè)臺(tái)站的空間分布研究與應(yīng)用［J］.煤炭學(xué)報(bào)，2010，35(12):2045-2048.

Jia Baoxin，Li Guozhen.Research and application for spatial distribution of mines seismic monitoring stations［J］. Journal of China Coal Society，2010，35(12):2045-2048.

［4］ 周志輝.三維數(shù)字化礦山信息系統(tǒng)的動(dòng)態(tài)更新與應(yīng)用［J］.現(xiàn)代礦業(yè)，2014(5):11-15.

Zhou Zhihui. Dynamic updating and application of the 3D digital mine information system［J］.Modern Mining，2014(5):11-15.

［5］ 蔣 晨，張書畢.基于VB 的TM30 測(cè)量機(jī)器人變形監(jiān)測(cè)系統(tǒng)開發(fā)［J］.金屬礦山，2015(1):104-107.

Jiang Chen，Zhang Shubi.Development of the deformation monitoring system of TM30 measuring robot based on VB［J］.Metal Mine，2015(1):104-107.

［6］ 吳立新，汪云甲，丁恩杰，等. 三論數(shù)字礦山——借力物聯(lián)網(wǎng)保障礦山安全與智能采礦［J］.煤炭學(xué)報(bào)，2012，37(3):357-365.

Wu Lixin，Wang Yunjia，Ding Enjie，et al. Thirdly study on digital mine:serve for mine safety and intellimine with support from internet of tings［J］.Journal of China Coal Society，2012，37(3):357-365.

［7］ 張廣偉.山區(qū)地表移動(dòng)觀測(cè)站布設(shè)新方法探討［J］. 金屬礦山，2015(4):182-185.

Zhang Guangwei.A new method of layout mode of surface movement observation station in mountain area［J］.Metal Mine，2015(4):182-185.

［8］ 武繼軍，胡圣武. 測(cè)量平差的基準(zhǔn)方程與相對(duì)形變及其精度分析［J］，煤炭學(xué)報(bào)，2007，32(8):808-812.

Wu Jijun，Hu Shengwu.Surveying adjustment datum and relative deformation accuracy analysis［J］. Journal of China Coal Society，2007，32(8):808-812.

［9］ 高俊強(qiáng)，胡 燦.Helmert 方差分量估計(jì)在隧道貫通中控制網(wǎng)平差的研究［J］.中國(guó)礦業(yè)大學(xué)學(xué)報(bào)，2006，35(1):125-129.

GaoO Junqiang，Hu Can.Study of Helmert variance component estimation in tunnel run-through control traverse difference［J］. Journal of China University of Mining＆Technology，2006，35(1):125-129.

［10］ 黃維彬. 近代平差理論及其應(yīng)用［M］. 北京:解放軍出版社，1992.

Huang Weibin. Modern Adjustment Theory and Its Application［M］.Beijing:People＇s Liberation Army Press，1992.

［11］ 張書畢.測(cè)量平差［M］.徐州:中國(guó)礦業(yè)大學(xué)出版社，2012.

Zhang Shubi.Surveying Adjustment［M］.Xuzhou:China University of Mining and Technology Press，2012.

［12］ Strang G，Borre K. Linear Algebra Geodesy and GPS［M］. Massachusetts:Wellesley-Cambridge Press，1997.

［13］ Frank Neitzel. Generalization of total least-squares on example of unweighted and weighted 2D similarity transformation［J］. Journal of Geodesy，2010，84(9):751-762.

［14］ 卞和方，張書畢，李益斌，等.概括平差模型的通用解法［J］.測(cè)繪科學(xué)，2009，34(6):35-36.

Bian Hefang，Zhang Shubi，Li Yibin，et al.A universal solution for generalization adjustment model［J］. Science of Surveying and Mapping，2009，34(6):35-36.

［15］ 朱建軍，田玉淼，陶肖靜.帶準(zhǔn)則參數(shù)的平差準(zhǔn)則及其統(tǒng)一與解算［J］.測(cè)繪學(xué)報(bào)，2012，41(1):8-13.

Zhu Jianjun，Tian Yumiao，Tao Xiaojing.United expression and solution of adjustment criteria with parameters［J］.Acta Geodaetica et Cartographic Sinica，2012，41(1):8-13.

［16］ 李昌貴，呂志平.分組平差網(wǎng)格計(jì)算研究［J］.測(cè)繪科學(xué)，2011，36(5):25-27.

Li Changgui，Lu Zhiping. Group adjustment based on grid computing［J］.Science of Surveying and Mapping，2011，36(5):25-27.

［17］ Teunissen P J G，Amiri-Simkooei A R.Least-squares variance component estimation［J］.Journal of Geodesy，2008，82(10):65-82.

［18］ 馬 洋，歐吉坤，袁運(yùn)斌.采用聯(lián)合平差法處理附有病態(tài)等式約束的反演問題［J］.武漢大學(xué)學(xué)報(bào):信息科學(xué)版，2011，36(7):816-819.

Ma Yang，Ou Jikun，Yuan Yunbin. Solving equality constraint inversion with ill-posed constraint matrix using united method［J］.Geomatics and Information Science of Wuhan University，2011，36(7):816-819.

［19］ 葛旭明，伍吉倉(cāng).病態(tài)總體最小二乘問題的廣義正則化［J］.測(cè)繪學(xué)報(bào)，2012，41(3):372-377.

Ge Xuming，Wu Jicang.Generalized regularization to ill-posed total least squares problem［J］.Acta Geodaetica et Cartographic Sinica，2012，41(3):372-377.

［20］ 李 斐，郝衛(wèi)峰，王文睿，等.非線性病態(tài)問題解算的擾動(dòng)分析［J］.測(cè)繪學(xué)報(bào)，2011，40(1):1-5.

Li Fei，Hao Weifeng，Wang Wenrui，et al.The perturbation analysis of nonlinear ill-conditioned solution［J］.Acta Geodaetica et Cartographic Sinica，2011，40(1):1-5.