許哲明

(南昌大學 建筑工程學院， 南昌 330031)

基于卡爾曼濾波的地表移動變形預測*

許哲明

(南昌大學 建筑工程學院， 南昌 330031)

針對傳統(tǒng)卡爾曼濾波模型在地表變形預測中依賴于噪聲和數(shù)學模型而導致的運行發(fā)散問題，提出了一種基于方差補償自適應(yīng)卡爾曼濾波的建筑地表移動變形預測分析方法.提出的自適應(yīng)卡爾曼濾波主要是通過適當?shù)墓浪愫托拚到y(tǒng)模型的不確定參數(shù)以及噪聲的統(tǒng)計特性，來彌補濾波過程中噪聲方差的不足.建立了自適應(yīng)卡爾曼濾波的數(shù)學模型，并在Matlab軟件上對地表移動變形預測進行了仿真分析，其結(jié)果驗證了所提卡爾曼濾波方法在地表移動變形預測應(yīng)用中的可行性，表明該方法有效提高了實時預測的可靠性和預測精度.

地表變形； 卡爾曼濾波； 方差補償； 自適應(yīng)； 預測； 數(shù)學模型； Matlab仿真； 可靠性

建筑物的大量開發(fā)會破壞周圍巖體的原始應(yīng)力狀態(tài)，從而引起地表巖層的移動變形，為地表建筑帶來較大的安全隱患.如果能夠準確實時地對地表建筑物上任意點的移動變形數(shù)據(jù)進行預測，則可以進行有計劃的維護和防御，防止次生災害的發(fā)生，因此對地表變形進行監(jiān)測和數(shù)據(jù)處理顯得十分重要[1-3].

目前，常用的地表變形監(jiān)測數(shù)據(jù)處理主要包括：回歸分析法[4-5]、灰色關(guān)聯(lián)分析法[6-7]以及卡爾曼濾波法[8-10]，其中回歸分析法在引入自變量過程中需對方程中各自變量做顯著性實驗，會造成回歸方程因子的引入不夠準確；灰色關(guān)聯(lián)分析法中各指標的量綱或單位常常不同，指標值的量綱級相差較大，且在計算關(guān)聯(lián)度時存在客觀性不足的問題；而卡爾曼濾波方法具有較強的自適應(yīng)性，且對原數(shù)據(jù)不依賴，具有最小無偏方差，能夠在不存儲大量監(jiān)測數(shù)據(jù)的條件下進行參數(shù)估計，因此，被廣泛應(yīng)用于地表變形預測數(shù)據(jù)的動態(tài)處理中.卡爾曼濾波是通過狀態(tài)方程以及觀測方程來構(gòu)造系統(tǒng)的動態(tài)過程，通過改變?yōu)V波增益矩陣來實現(xiàn)監(jiān)測數(shù)據(jù)的識別和提取[11-13].

本文根據(jù)某建筑物地表變形數(shù)據(jù)，提出基于方差自適應(yīng)的卡爾曼濾波方法，對監(jiān)測數(shù)據(jù)進行計算處理而得到其預測值，該方法能夠有效提高預測精度，具有較強的可行性.

1 卡爾曼濾波模型

卡爾曼濾波模型可實現(xiàn)系統(tǒng)動態(tài)數(shù)據(jù)處理，它是利用觀測值來估算變化的狀態(tài)向量，通過對新的觀測值不斷修正和預測得到新的預測值，不需要存儲大量的歷史數(shù)據(jù)[14-15].卡爾曼濾波離散系統(tǒng)的狀態(tài)方程和觀測方程分別表示為

Xk=Φk，k-1Xk-1+Γk，k-1Ωk-1

(1)

Lk=BkXk+Δk

(2)

式中：Xk，Xk-1分別為第k時刻和第k-1時刻的狀態(tài)向量；Φk，k-1為第k次和第k-1次之間的狀態(tài)轉(zhuǎn)移矩陣；Γk，k-1為第k次和第k-1次的動態(tài)噪聲矩陣；Ωk-1為第k-1次動態(tài)噪聲向量；Lk為第k次的狀態(tài)觀測向量；Bk為第k次的狀態(tài)觀測向量系數(shù)矩陣；Δk為第k次的狀態(tài)觀測噪聲向量.

若Ωk-1和Δk是互不相關(guān)的零均值白噪聲序列，則兩者滿足

(3)

估算可得卡爾曼濾波的遞推方程式為

(4)

(5)

(6)

(7)

2 自適應(yīng)卡爾曼濾波模型

傳統(tǒng)卡爾曼濾波模型在很多情況下依賴于噪聲和數(shù)學模型，需要深入研究動態(tài)系統(tǒng)的噪聲特性以及數(shù)學模型特性，因此會產(chǎn)生各種誤差且出現(xiàn)運行發(fā)散的結(jié)果，達不到預期目的，很難滿足實際應(yīng)用場合.針對上述問題，自適應(yīng)卡爾曼濾波被提出并應(yīng)用，自適應(yīng)卡爾曼濾波主要是通過適當?shù)毓浪愫托拚到y(tǒng)模型的不確定參數(shù)以及噪聲的統(tǒng)計特性，有效減小誤差，使預測結(jié)果更精確.在地表變形監(jiān)測數(shù)據(jù)濾波中，可根據(jù)已知信息實時估算動態(tài)噪聲的方差來彌補濾波過程中噪聲方差的不足.

假設(shè)Ωk、Δk都是正態(tài)序列，則可定義第i步的預測殘差為

(8)

(9)

則預測殘差Vk+i的方差矩陣可表示為

(10)

令

Bk+iΦk+i，rΓr，r-1=A(k+i，r)=[a(k+i，r)]

(11)

式中：r=1，2，…，N；k=1，2，…，n.若假設(shè)DΩr-1Ωr-1在tk+1，tk+2，…，tk+N觀測時間段內(nèi)是一常數(shù)對角陣，即

(12)

(13)

(14)

則可以得到

F=AdiagDΩr-1Ωr-1+ξ

(15)

式(15)可以看作是diagDΩr-1Ωr-1的線性方程組形式，當n≥r時，該方程組存在唯一解，則diagDΩr-1Ωr-1的估計值可表示為

diagDΩr-1Ωr-1=(ATA)-1ATF

(16)

因此，可根據(jù)式(16)求得任意時間段內(nèi)的DΩr-1Ωr-1，并實時估算動態(tài)噪聲協(xié)方差.

3 實例分析

以某建筑物作為研究對象，分析了建筑物變形對該地表的影響.較高頻率的觀測對地表移動變形的規(guī)律研究十分有利，故對該段地表進行了21期的觀測，選擇觀測間隔為每天1次.本文以工作面走向線與建筑物走向線的交點作為觀測站，首先以前16期觀測數(shù)據(jù)為基礎(chǔ)，利用Matlab軟件編制自適應(yīng)卡爾曼濾波程序，對地表變形規(guī)律進行擬合分析，然后對后5期的數(shù)據(jù)進行預測，并將其與觀測值以及傳統(tǒng)的卡爾曼濾波預測方法進行了對比.

兩種預測值與實際觀測值的水平位移對比結(jié)果如表1所示，實測值與誤差值的對比曲線如圖1、2所示.由圖1、2可見，自適應(yīng)卡爾曼濾波方法對地表位移預測的結(jié)果更接近于原始的實測數(shù)據(jù)，最大誤差百分比僅為0.25%，而傳統(tǒng)方法最大誤差百分比為0.38%(均為第17期觀測值)，較傳統(tǒng)的卡爾曼濾波預測方法而言，其具有更高的預測精度，說明采用基于協(xié)方差補償自適應(yīng)卡爾曼濾波能夠有效模擬位移狀態(tài)變量的變化規(guī)律.

表1 觀測點水平位移預測結(jié)果Tab.1 Prediction results of horizontal displacement of observation points

圖1 監(jiān)測點水平位移的預測值和實測值對比曲線Fig.1 Contrast curves for predicted and measured valuesof horizontal displacement of monitoring points

兩種預測值與實際觀測值的沉降對比結(jié)果如表2所示，實測值與誤差值的沉降對比曲線如圖3、4所示.由表2可以看出，自適應(yīng)卡爾曼濾波方法對地表沉降預測結(jié)果的最大誤差百分比為0.28%，而傳統(tǒng)的卡爾曼濾波預測方法最大誤差百分比為0.64%(均為第21期觀測值)，證明自適應(yīng)卡爾曼濾波方法具有更高的預測精度，能夠有效模擬地表沉降狀態(tài)變量的變化規(guī)律.

圖2 監(jiān)測點水平位移預測誤差對比曲線Fig.2 Contrast curves for prediction error of horizontaldisplacement of monitoring points

表2 觀測點沉降預測結(jié)果Tab.2 Prediction results of settlement of monitoring point

圖3 監(jiān)測點沉降預測值和實測值對比曲線Fig.3 Contrast curves for predicted and measuredvalues of settlement of monitoring points

圖4 監(jiān)測點沉降預測誤差對比曲線Fig.4 Contrast curves for prediction errorof settlement of monitoring points

4 結(jié) 論

本文結(jié)合卡爾曼濾波原理和建筑物結(jié)構(gòu)特點，建立了基于方差補償?shù)淖赃m應(yīng)卡爾曼濾波預測模型，實現(xiàn)了對建筑物地表水平位移和沉降的預測.仿真結(jié)果表明：

1) 自適應(yīng)卡爾曼濾波通過遞推公式對濾波模型進行不斷修正，有利于對實時觀測數(shù)據(jù)進行處理，非常適合于地表變形監(jiān)測數(shù)據(jù)的處理和預測；

2) 自適應(yīng)卡爾曼濾波對地表變形數(shù)據(jù)的預測結(jié)果非常接近于實測值，比傳統(tǒng)的卡爾曼濾波預測具有更高的精度，能夠較好地預測地表狀態(tài)變量的變化規(guī)律.

[1] 張紅娟，祝文華，魏向輝，等.基于VRS技術(shù)的煤礦地表變形監(jiān)測 [J].地理空間信息，2015，13(2)：150-152.

(ZHANG Hong-juan，ZHU Wen-hua，WEI Xiang-hui，et al.Mining surface deformation monitoring based on VRS technology [J].Geospatial Information，2015，13(2)：150-152.)

[2] 楊理踐，李暉，靳鵬，等.管道地理位置測量系統(tǒng)的動態(tài)初始對準方法 [J].沈陽工業(yè)大學學報，2015，37(6)：656-661.

(YANG Li-jian，LI Hui，JIN Peng，et al.Dynamic initial alignment method for geographic position mea-surement system of pipeline [J].Journal of Shenyang University of Technology，2015，37(6)：656-661.)

[3] 夏開宗，陳從新，夏天游，等.結(jié)構(gòu)面對程潮鐵礦西區(qū)地表變形的影響分析 [J].巖土力學，2015，36(5)：1389-1396.

(XIA Kai-zong，CHEN Cong-xin，XIA Tian-you，et al.Analysis of influence of structural planes on surface deformation in western area of chengchao iron mine [J].Rock and Soil Mechanics，2015，36(5)：1389-1396.)

[4] 曲芳，明向蘭.基于回歸分析法的黑龍江省人口預測及應(yīng)用 [J].工業(yè)儀表與自動化裝置，2015(3)：66-68.

(QU Fang，MING Xiang-lan.Heilongjiang province population forecast and application based on regression analysis method [J].Industrial Instrumentation & Automation，2015(3)：66-68.)

[5] 張智韜，蘭玉彬，鄭永軍，等.影響大豆 NDVI 的氣象因素多元回歸分析 [J].農(nóng)業(yè)工程學報，2015，31(5)：188-193.

(ZHANG Zhi-tao，LAN Yu-bin，ZHENG Yong-jun，et al.Multiple regression analysis of soybean NDVI affected by meteorological [J].Transactions of the Chinese Society of Agricultural Engineering，2015，31(5)：188-193.)

[6] 董海彪，盧文喜，安永凱，等.基于對應(yīng)分析法的吉林西部地下水化學特征時空演化研究 [J].環(huán)境科學學報，2016，36(3)：820-826.

(DONG Hai-biao，LU Wen-xi，AN Yong-kai，et al.Spatio-temporal evolution of chemical characteristics of groundwater in western Jilin based on correspondence analysis method [J].Acta Scientiae Circumstantiae，2016，36(3)：820-826.)

[7] 羅德江，吳尚昆，郭科.基于組合權(quán)-灰色關(guān)聯(lián)分析法的礦產(chǎn)資源開發(fā)利用綜合評價 [J].金屬礦山，2015，44(2)：20-25.

(LUO De-jiang，WU Shang-kun，GUO Ke.The comprehensive evaluation of exploitation and utilization efficiency of mineral resources based on combination weight and gray relation analysis method [J].Metal Mine，2015，44(2)：20-25.)

[8] 陳鋒.卡爾曼濾波和卡爾曼預測方法的改進及其在結(jié)構(gòu)損傷識別中的應(yīng)用 [D].廈門：廈門大學，2014.

(CHEN Feng.Improvement of Kalman filter and Kalman estimator in the application of structural damage detection [D].Xiamen：Xiamen University，2014.)

[9] 駱榮劍，李穎，錢廣華，等.機動目標跟蹤中一種改進的自適應(yīng)卡爾曼濾波算法 [J].重慶郵電大學學報(自然科學版)，2015，27(1)：31-36.

(LUO Rong-jian，LI Ying，QIAN Guang-hua，et al.Improved maneuvering target tracking adaptive Kalman filter algorithm [J].Journal of Chongqing University of Posts and Telecommunications (Natural Science Edition)，2015，27(1)：31-36.)

[10]王利，李亞林，劉萬紅.卡爾曼濾波在大壩動態(tài)變形監(jiān)測數(shù)據(jù)處理中的應(yīng)用 [J].西安科技大學學報，2006，26(3)：353-357.

(WANG Li，LI Ya-lin，LIU Wan-hong.Application of Kalman filtering in data processing of dam dynamic deformation monitoring [J].Journal of Xi’an University of Science and Technology，2006，26(3)：353-357.)

[11]胡振濤，胡玉梅，劉先省.量測提升卡爾曼濾波 [J].電子學報，2016，44(5)：1149-1155.

(HU Zhen-tao，HU Yu-mei，LIU Xian-sheng.Kalman filter based on measurement lifting strategy [J].Acta Electronica Sinica，2016，44(5)：1149-1155.)

[12]袁捷，唐龍，杜浩.機場道面使用性能的動態(tài)自回歸預測模型 [J].同濟大學學報(自然科學版)，2015，43(3)：399-404.

(YUAN Jie，TANG Long，DU Hao.Dynamic auto-regression prediction model of airport pavement performance [J].Journal of Tongji University(Natural Science Edition)，2015，43(3)：399-404.)

[13]林琳，史良勝，宋雪航.地下水參數(shù)反演的確定性集合卡爾曼濾波方法 [J].武漢大學學報(工學版)，2016，49(2)：161-167.

(LIN Lin，SHI Liang-sheng，SONG Xue-hang.A deterministic ensemble kalman filter method for inversing hydrogeological parameters [J].Engineering Journal of Wuhan University，2016，49(2)：161-167.)

[14]王琦，孫華，李偉華，等.卡爾曼濾波在變形監(jiān)測數(shù)據(jù)處理中的應(yīng)用[J].工程地球物理學報，2009，6(5)：658-661.

(WANG Qi，SUN Hua，LI Wei-hua，et al.Application of Kalman filter analysis method in the deformation monitoring data procession [J].Chinese Journal of Engineering Geophysics，2009，6(5)：658-661.)

[15]許國輝，余春林.卡爾曼濾波模型的建立及其在施工變形測量中的應(yīng)用 [J].測繪通報，2004(4)：22-24.

(XU Guo-hui，YU Chun-lin.Establishment of Kalman filtering model and its application in construction deformation measurement [J].Bulletin of Surveying and Mapping，2004(4)：22-24.)

(責任編輯：景 勇 英文審校：尹淑英)

PredictionofsurfacemovementdeformationbasedonKalmanfiltering

XU Zhe-ming

(School of Architectural Engineering, Nanchang University, Nanchang 330031, China)

Aiming at the operation divergence problem caused by the fact that the traditional Kalman filtering model in the prediction of surface deformation depends on the noise and mathematical models, a prediction analysis method of building surface movement deformation based on adaptive Kalman filtering with variance compensation was proposed. The proposed adaptive Kalman filtering could compensate the deficiency of noise variance in the filtering process mainly through appropriately estimating and correcting the uncertain parameters of the system model and the statistical properties of the noise. The mathematical model for adaptive Kalman filtering was established, and the simulation analysis for the prediction of surface movement deformation was performed with Matlab software. The results verify the feasibility of the proposed Kalman filtering method in the prediction of surface movement deformation. The proposed method effectively improves the reliability and forecasting precision of real-time prediction.

surface deformation; Kalman filtering; variance compensation; adaption; prediction; mathematical model; Matlab simulation; reliability

TM 135

： A

： 1000-1646(2017)05-0557-05

2016-07-04.

江西省青年科學基金資助項目(20151BAB206057).

許哲明(1974-)，男，江西南昌人，講師，碩士，主要從事建筑、橋梁工程測量等方面的研究.

* 本文已于2017-03-28 17∶08在中國知網(wǎng)優(yōu)先數(shù)字出版. 網(wǎng)絡(luò)出版地址： http：∥www.cnki.net/kcms/detail/21.1189.T.20170328.1708.022.html

10.7688/j.issn.1000-1646.2017.05.15