GPS/INS組合系統(tǒng)的多重漸消魯棒容積卡爾曼濾波

張秋昭，張書畢，鄭南山，等

摘要：目的：在GPS/INS組合導(dǎo)航中一般采用卡爾曼濾波對系統(tǒng)的狀態(tài)進行估計，在實際應(yīng)用中，由于慣性導(dǎo)航系統(tǒng)模型的非線性及噪聲統(tǒng)計特性的不確定性，采用傳統(tǒng)卡爾曼濾波器，影響系統(tǒng)的定位精度。本文基于多重漸消濾波和非線性 H∞濾波原理，構(gòu)造另一種新的非線性魯棒濾波方法，以提高組合定位結(jié)果的精度和穩(wěn)定性。方法：在對適用于GPS/INS組合導(dǎo)航非線性模型的容積卡爾曼濾波進行深入研究的基礎(chǔ)上，提出一種改進的多重漸消 H∞魯棒容積卡爾曼濾波算法。基于系統(tǒng)狀態(tài)的可觀測性給出多重漸消因子矩陣求解過程，提高濾波算法的穩(wěn)定性，抑制濾波發(fā)散；引入H∞魯棒思想，構(gòu)造多重漸消 H∞魯棒容積卡爾曼濾波器；采用一種奇異值分解的矩陣分解策略代替標(biāo)準(zhǔn)容積卡爾曼濾波中的Cholesky分解，進一步提高算法的數(shù)值穩(wěn)定性。結(jié)果：采集了一組車載GPS/INS導(dǎo)航實驗數(shù)據(jù)，驗證論文提出的新濾波算法的有效性。論文采用5種算法進行組合導(dǎo)航數(shù)據(jù)處理，分別為：標(biāo)準(zhǔn)容積卡爾曼濾波（CKF）；多重漸消容積卡爾曼濾波（MF-CKF）；H∞魯棒容積卡爾曼濾波（R-CKF）；多重漸消魯棒容積卡爾曼濾波（MFRCKF）；基于 SVD 的多重漸消 H∞魯棒容積卡爾曼濾波（SVD-MFR-CKF）。從5種算法的位置誤差圖可以看出，標(biāo)準(zhǔn)容積卡爾曼濾波算法能正常收斂，這說明標(biāo)準(zhǔn)容積卡爾曼濾波算法能處理組合導(dǎo)航系統(tǒng)的非線性問題，但系統(tǒng)穩(wěn)定性和抗野值能力較弱。比較多重漸消容積卡爾曼濾波算法和標(biāo)準(zhǔn)容積卡爾曼濾波算法的誤差圖，可以發(fā)現(xiàn)多重漸消容積卡爾曼濾波能改善濾波的穩(wěn)定性，但由于觀測野值影響，部分歷元的誤差仍較大，算法的魯棒性還需要進一步提高。對比H∞魯棒容積卡爾曼濾波算法與標(biāo)準(zhǔn)容積卡爾曼濾波和多重漸消容積卡爾曼濾波的誤差圖結(jié)果，可以看出，H∞魯棒容積卡爾曼濾波算法能夠有效抑制觀測野值的影響，說明H∞魯棒容積卡爾曼濾波算法抗野值能力較強。對比多重漸消魯棒容積卡爾曼濾波算法與 H∞魯棒容積卡爾曼濾波和多重漸消容積卡爾曼濾波算法的位置誤差圖可知，多重漸消魯棒容積卡爾曼濾波算法的位置誤差精度明顯提高，但濾波結(jié)果整體波動性仍較明顯，濾波的穩(wěn)定性還需要進一步提高。對比基于SVD的多重漸消H∞魯棒容積卡爾曼濾波算法與其他 5種算法的位置誤差圖可知，基于 SVD的多重漸消 H∞魯棒容積卡爾曼濾波算法的位置誤差明顯減小，各個歷元的位置誤差數(shù)值也較小，且沒有明顯的野值，說明該算法不僅提高了濾波穩(wěn)定性，野值也得到了明顯的抑制。從5種不同數(shù)據(jù)處理方案的統(tǒng)計比較表中可以看出，多重漸消容積卡爾曼濾波算法和H∞魯棒容積卡爾曼濾波的位置解算精度比標(biāo)準(zhǔn)容積卡爾曼濾波算法都有所提高，而基于SVD的多重漸消H∞魯棒容積卡爾曼濾波效果最好，3個方向的均方根誤差（root mean square error，RMS）與標(biāo)準(zhǔn)CKF相比分別減小了55.8%，46.6%和39.7%，從數(shù)值統(tǒng)計角度證明了該算法的優(yōu)越性，與理論分析一致。從不同算法新息序列的偏態(tài)系數(shù)和峰態(tài)系數(shù)表可以看出，多重漸消容積卡爾曼濾波、H∞魯棒容積卡爾曼濾波、多重漸消魯棒容積卡爾曼濾波和基于SVD的多重漸消H∞魯棒容積卡爾曼濾波4種改進的算法的正態(tài)性較標(biāo)準(zhǔn)容積卡爾曼濾波的正態(tài)性均有所提高，其中基于SVD的多重漸消H∞魯棒容積卡爾曼濾波算法更能提高濾波解算的穩(wěn)定性。結(jié)論：標(biāo)準(zhǔn)容積卡爾曼濾波算法能處理 GPS/INS組合導(dǎo)航系統(tǒng)的非線性誤差模型，但濾波穩(wěn)定性和抗野值能力較差；多重漸消容積卡爾曼濾波算法，一定程度上能提高組合導(dǎo)航系統(tǒng)的濾波穩(wěn)定性，H∞魯棒容積卡爾曼濾波算法對奇異觀測值有一定的抑制作用；基于SVD的多重漸消H∞魯棒容積卡爾曼濾波算法，既能提高濾波器的穩(wěn)定性，也能抑制奇異值的影響，能有效提高組合導(dǎo)航的定位精度；基于SVD的多重漸消H∞魯棒容積卡爾曼濾波算法與標(biāo)準(zhǔn)容積卡爾曼濾波算法相比，XYZ這3個方向的位置精度分別提高了55.8%， 46.6%和39.7%。

來源出版物：中國礦業(yè)大學(xué)學(xué)報， 2014， 43(1)： 162-168

入選年份：2016

圓形巷道圍巖偏應(yīng)力場及塑性區(qū)分布規(guī)律研究

馬念杰，李季，趙志強

摘要：目的：在地下工程和采礦工程中，圓形巷道圍巖應(yīng)力計算一直以來都參照彈塑性力學(xué)中的圓形孔周圍應(yīng)力解。當(dāng)側(cè)壓系數(shù)不為1即非均勻圍壓條件下時，被普遍應(yīng)用的芬納公式或卡斯特奈公式無法計算圓形巷道圍巖塑性區(qū)半徑。同時在開挖巷道過程中，會使巷道圍巖產(chǎn)生偏應(yīng)力，而偏應(yīng)力控制著巖體的變形與破壞。因此，本文以圓孔應(yīng)力解和偏應(yīng)力理論為基礎(chǔ)，利用莫爾-庫侖強度準(zhǔn)則，研究非均勻圍壓條件下圍巖偏應(yīng)力場和塑性區(qū)分布規(guī)律。方法：利用彈塑性力學(xué)彈性平板中的圓孔周圍的二維應(yīng)力分布解，研究分析了均質(zhì)圍巖條件下圓形巷道圍巖的應(yīng)力狀態(tài)。并結(jié)合塑性力學(xué)中的偏應(yīng)力理論，推導(dǎo)得到了圓形巷道圍巖任意一點最大與最小主偏應(yīng)力的理論計算公式。利用該公式，研究了不同埋深、不同圍巖位置及不同側(cè)壓系數(shù)下的圓形巷道圍巖偏應(yīng)力場的分布規(guī)律。在此基礎(chǔ)之上，結(jié)合莫爾—庫侖強度理論，推到了非均勻圍巖條件下，圓形巷道圍巖塑性區(qū)半徑計算公式。利用該公式，計算了不同側(cè)壓系數(shù)、不同圍巖力學(xué)參數(shù)及不同巷道尺寸下的圍巖塑性區(qū)半徑，并繪制了塑性區(qū)分布圖。以此為基礎(chǔ)，利用理論分析研究了非均勻應(yīng)力場下圓形巷道塑性區(qū)的分布規(guī)律。結(jié)果：圍巖偏應(yīng)力場分布規(guī)律研究表明：在均勻應(yīng)力場，巷道埋深和巷道半徑?jīng)Q定了圓形巷道圍巖的偏應(yīng)力大小。隨著巷道埋深的增加，圍巖最大及最小主偏應(yīng)力大小都隨之呈線性增長的趨勢，且最大與最小主應(yīng)力的差值越來越大。圓形巷道圍巖偏應(yīng)力與極坐標(biāo)r的關(guān)系呈橫“八”字型分布，隨著極坐標(biāo)r的增加，最大與最小主偏應(yīng)力的差值越來越小。圓形巷道圍巖偏應(yīng)力與極坐標(biāo)θ的關(guān)系呈橫“筆尖”型分布，隨著極坐標(biāo)θ的增加，最大與最小主偏應(yīng)力的差值先越來越大，隨后逐漸減小。圓形巷道圍巖任意位置的偏應(yīng)力與側(cè)壓系數(shù)的關(guān)系呈“X”型分布，隨著側(cè)壓系數(shù)的增大，最大與最小主偏應(yīng)力的差值先減小后增大。巷道圍巖塑性區(qū)分布規(guī)律表明：隨著側(cè)壓系數(shù)的增大，水平方向的塑性區(qū)減小，而垂直方向的塑性區(qū)增大，在非均勻應(yīng)力場下，圓形巷道圍巖塑性區(qū)會出現(xiàn)蝶形分布。隨著圍巖內(nèi)摩擦角的不斷增加，圍巖塑性區(qū)的減小幅度卻逐漸減小，其對圍巖塑性區(qū)大小的影響作用將會逐漸減弱。隨著圍巖內(nèi)聚力的不斷增加，圍巖塑性區(qū)的減小幅度逐漸減小，圍巖內(nèi)聚力對圍巖塑性區(qū)大小的影響作用將會逐漸減弱。隨著巷道半徑的增加，巷道圍巖相同部位的塑性區(qū)大小增大幅度都一樣，這一增加趨勢與內(nèi)聚力、內(nèi)摩擦角變化時，塑性區(qū)的變化趨勢顯著不同。結(jié)論：（1）以彈塑性力學(xué)中的圓孔應(yīng)力解和塑性力學(xué)的偏應(yīng)力理論為基礎(chǔ)，利用莫爾-庫侖強度準(zhǔn)則，得到了非均勻應(yīng)力場下圓形巷道圍巖偏應(yīng)力的計算公式和塑性區(qū)半徑計算方法。（2）巷道埋深、極坐標(biāo)、極坐標(biāo)和側(cè)壓系數(shù)與最大和最小主偏應(yīng)力的大小曲線分別呈線性分布、“八”字型分布、“筆尖”型和“X”型分布；側(cè)壓系數(shù)小于和大于 1時，會出現(xiàn)不同形態(tài)的蝶形塑性區(qū)；隨著側(cè)壓系數(shù)的增大，塑性區(qū)邊界各位置的最大主偏應(yīng)力會逐漸增大，而最小主偏應(yīng)力卻逐漸減小，兩者之間的差值越來越大；內(nèi)摩擦角、內(nèi)聚力和圓形巷道半徑的增加，都會導(dǎo)致塑性區(qū)半徑的減小。

來源出版物：中國礦業(yè)大學(xué)學(xué)報， 2015， 44(2)： 206-213

入選年份：2016

深部軟巖巷道錨注支護機理數(shù)值模擬研究

孟慶彬，韓立軍，喬衛(wèi)國，等

摘要：目的：隨著礦產(chǎn)資源進入深部開采階段，深部軟巖巷道圍巖大變形、難支護等難題越凸顯，錨注支護技術(shù)可改善巷道圍巖的受力狀態(tài)，提高破裂圍巖的自承能力與穩(wěn)定性，是解決深部軟巖巷道支護難題的有效技術(shù)途徑。本文針對破裂巖體注漿加固力學(xué)特性試驗與巷道錨注支護機理研究難題，采用破裂巖體注漿加固力學(xué)特征試驗、數(shù)值模擬及理論分析等綜合技術(shù)手段，力圖揭示深部軟巖巷道錨注支護機理。方法：首先，為研究破裂巖體注漿加固效應(yīng)，自主研制了破裂巖樣承壓注漿試驗設(shè)備（由圍壓注漿系統(tǒng)、活塞加壓系統(tǒng)、手動注漿泵、高壓閥門、高壓軟管及壓力表等組成），采用該注漿試驗裝置與 TATW-2000型微機控制電液伺服巖石三軸試驗機，開展了兩類破裂巖體注漿加固力學(xué)特性試驗，對比分析了破裂巖樣注漿加固前后的力學(xué)特性；并采用KH-3000型視頻儀對破裂巖樣裂紋及其注漿后的狀態(tài)進行了探測，得出水泥注漿加固對巖石中的大孔徑孔隙的充填密實作用明顯。然后，結(jié)合巷道錨注支護機理的理論分析，且考慮注漿后破裂巖體的內(nèi)聚力、內(nèi)摩擦角及彈性模量等參數(shù)的提高在數(shù)值模擬分析中如何體現(xiàn)，提出了“錨注加固體等效層”的概念，采用FLAC3D模擬揭示了“錨注加固體等效層”厚度、彈性模量、內(nèi)聚力、內(nèi)摩擦角等參數(shù)變化對巷道圍巖位移及塑性區(qū)的影響規(guī)律。結(jié)果：（1）采用破裂巖樣承壓注漿試驗設(shè)備，開展了兩類破裂巖體注漿加固力學(xué)特性試驗，破裂巖樣注漿后，所形成的注漿固結(jié)體與破裂巖體殘余強度相比有較大幅度的提高，約為1.5～2.0倍。注漿量隨著圍壓的提高而減小，注漿加固體強度隨圍壓的提高而增大，即說明圍壓對破裂巖體起了較大的約束作用。（2）錨注支護是在軟巖巷道原有支護的基礎(chǔ)上進行壁后注漿，對原有支護結(jié)構(gòu)有強化作用。錨注支護可配合錨網(wǎng)噴支護，形成一個多層有效組合拱，即錨桿錨固壓縮區(qū)組合拱、噴網(wǎng)組合拱和漿液擴散加固拱，形成的多層有效組合拱提高了支護結(jié)構(gòu)的整體性和承載能力，擴大了支護結(jié)構(gòu)的承載范圍。（3）采用注漿加固后，破裂巖樣的裂隙被漿液充填，使破裂的巖石黏結(jié)成整體，提高了破裂巖體的強度，使得莫爾圓遠離強度包絡(luò)線，有利于巷道圍巖的穩(wěn)定。將漿液在巖石裂隙中擴散范圍（徑向擴散半徑）的巖層定義為“錨注加固體等效層”。在數(shù)值計算時，將“錨注加固體等效層”范圍內(nèi)的圍巖體參數(shù)按破裂巖體注漿加固試驗結(jié)果進行調(diào)整，以體現(xiàn)錨注加固后圍巖體參數(shù)的改變影響。（4）隨著“錨注加固體等效層”參數(shù)的提高，巷道圍巖收斂變形量不斷減小，頂板下沉量最大降低幅度為18.36%～28.72%，底板底臌量最大降低幅度為15.26%～24.45%，兩幫內(nèi)擠量最大降低幅度為19.15%～29.1%。巷道圍巖塑性區(qū)范圍也逐漸減少，圍巖塑性區(qū)范圍最大降低幅度為44.0%～83.33%，且圍巖塑性區(qū)分布趨于均勻，防止因圍巖局部破壞而導(dǎo)致整體失穩(wěn)。但當(dāng)內(nèi)聚力提高到2.5倍之后，巷道圍巖變形降低幅度較小、曲線斜率趨于平緩。結(jié)論：本文自主研制了破裂巖樣承壓注漿試驗設(shè)備，開展了破裂巖體注漿加固力學(xué)特性試驗，分析了破裂巖樣注漿加固前后的力學(xué)特性與微觀結(jié)構(gòu)。提出了“錨注加固體等效層”概念，采用FLAC3D模擬揭示了“錨注加固體等效層”的厚度、彈性模量、內(nèi)聚力、內(nèi)摩擦角等參數(shù)變化對巷道圍巖位移及塑性區(qū)的影響規(guī)律，體現(xiàn)了注漿對破裂巖體的完整性提高與強度強化效應(yīng)，為錨注支護參數(shù)優(yōu)化設(shè)計奠定了基礎(chǔ)。

來源出版物：采礦與安全工程學(xué)報， 2016， 33(1)： 27-34

入選年份：2016

深部軟巖巷道深—淺耦合全斷面錨注支護研究

王連國，陸銀龍，黃耀光，等

摘要：目的：深部巖土體所處地質(zhì)環(huán)境和所受地質(zhì)構(gòu)造影響的差異性和多樣性，導(dǎo)致軟巖巷道變形和破壞特征呈不確定性。特別是深部地質(zhì)構(gòu)造發(fā)育、礦壓顯現(xiàn)劇烈區(qū)內(nèi)的軟巖巷道，由于圍巖極為破碎、節(jié)理裂隙發(fā)育、強度極低，在高自重應(yīng)力、高構(gòu)造應(yīng)力和強烈開挖擾動等作用下，圍巖應(yīng)力急劇向巷道深部轉(zhuǎn)移，發(fā)生應(yīng)力重新分布，很大程度上促進了巷道深淺部圍巖的孔隙和裂隙的再發(fā)育，降低了圍巖的強度。這不但使巷道淺部圍巖變形和破壞異常嚴重，并且使得巷道深部圍巖的遭到整體性破壞，其承載力和強度都進一步降低。因此，進一步深入研究深部軟巖巷道圍巖穩(wěn)定性控制技術(shù)，對于復(fù)雜地質(zhì)構(gòu)造發(fā)育區(qū)的巷道支護和維護具有重要的理論意義和工程實用價值。方法：本文以安徽淮北某深部礦井的巷道支護工程為背景，提出了針對破碎區(qū)軟巖巷道的深—淺耦合全斷面錨注支護方法，分析了深—淺耦合錨注支護作用機理。利用滲流力學(xué)理論推得圍巖擾動應(yīng)力和注漿壓力耦合作用下的漿液非穩(wěn)態(tài)滲透擴散基本方程，并運用多場耦合軟件 COMSOL建立了錨注漿液在圍巖中滲透擴散的數(shù)值計算模型。在相同注漿壓力和注漿時間的條件下，研究注漿錨索和注漿錨桿單獨作用下的錨注漿液在巷道圍巖中的擴散規(guī)律和擴散范圍。在此基礎(chǔ)上，進一步研究注漿錨桿和注漿錨索耦合作用下巷道全斷面圍巖內(nèi)的漿液擴散規(guī)律和深—淺耦合加固作用機理，進而完善了深—淺耦合錨注支護理論及技術(shù)體系。為了評價深—淺耦合錨注聯(lián)合支護的加固效果，將該支護技術(shù)體系應(yīng)用于某采區(qū)軌道大巷的支護設(shè)計中，主要對深—淺耦合錨注加固后的巷道深部和表面圍巖的變形進行了現(xiàn)場監(jiān)測，驗證其合理性。結(jié)果：提出了針對深部軟巖巷道圍巖的以注漿錨索和注漿錨桿為核心的“深—淺耦合全斷面錨注支護體系”，建立了深—淺耦合錨注作用下的漿液滲流基本方程，并借助COMSOL多場耦合軟件獲得了深—淺耦合全斷面錨注支護下漿液的擴散規(guī)律以及深—淺耦合加固作用機理，（1）錨注漿液以各自鉆孔為中心，呈橢圓形狀向外擴散，且隨著注漿錨桿和注漿錨索的增加而發(fā)生孔間相互滲透、交叉擴散的現(xiàn)象；（2）僅采用注漿錨桿對巷道進行錨注支護時，錨注漿液只能在巷道淺部圍巖裂隙和孔隙內(nèi)進行滲透擴散，最終在巷道淺部圍巖內(nèi)形成小范圍的錨注加固圈；（3）當(dāng)對巷道圍巖采用深—淺耦合錨注支護時，錨注漿液不僅能在巷道淺部圍巖內(nèi)形成完整錨注加固圈，而且注漿錨索所注漿液能同時在巷道淺部、深部圍巖內(nèi)進行滲透擴散，使得深部破碎圍巖內(nèi)的裂隙和孔隙結(jié)構(gòu)等能被充填、密實，進而形成局部錨注加固體；(4）隨著注漿錨索和注漿錨桿密度的增加，巷道淺部圍巖內(nèi)的錨注加固圈厚度逐漸增大，分別達到1.15，2.35和3.05 m。現(xiàn)場監(jiān)測證明該支護體系有效地控制了巷道圍巖的變形與破壞，驗證其合理性。結(jié)論：（1）深—淺耦合全斷面錨注時，巷道深部圍巖內(nèi)形成的局部錨注加固體與淺部圍巖內(nèi)的錨注加固圈協(xié)調(diào)作用，形成互為支撐的承載體，從而達到深—淺耦合錨注加固圍巖的目的；（2）采用以“中空注漿錨索和高強注漿錨桿”為核心的新型全斷面錨注支護技術(shù)體系能夠有效地控制了深部破碎區(qū)軟巖巷道圍巖的變形與破壞，為地質(zhì)構(gòu)造發(fā)育、礦壓顯現(xiàn)劇烈區(qū)軟巖巷道的支護設(shè)計提供了一種可借鑒的新方法。

來源出版物：中國礦業(yè)大學(xué)學(xué)報， 2016， 45(1)： 11-18

入選年份：2016