GBInSAR在大壩變形監(jiān)測中的應(yīng)用

郭樂萍，岳建平，岳 順，邱志偉，汪學(xué)琴

(河海大學(xué) 地球科學(xué)與工程學(xué)院，江蘇 南京 211100)

?

GBInSAR在大壩變形監(jiān)測中的應(yīng)用

郭樂萍，岳建平，岳 順，邱志偉，汪學(xué)琴

(河海大學(xué) 地球科學(xué)與工程學(xué)院，江蘇 南京 211100)

評判水工建筑物的工作性態(tài)，需高精度的變形、應(yīng)力等信息。根據(jù)大壩安全監(jiān)測特點和要求，介紹GBSAR監(jiān)測系統(tǒng)構(gòu)成和測量基本原理，分析該系統(tǒng)測量誤差來源及特征，探討該技術(shù)在大壩變形監(jiān)測中的可行性。以隔河巖大壩3 d的實測數(shù)據(jù)為例，分析該技術(shù)實際測量精度，研究顯示成果達亞毫米級，能夠滿足大壩變形監(jiān)測的精度要求。

變形監(jiān)測；GBInSAR；大壩；視線向位移；時間序列

我國自進入21世紀(jì)以來重點發(fā)展水利工程建設(shè)，興建一批以小浪底、紫坪鋪等水利樞紐為代表的水利重點工程項目。為了保證這些水利工程建筑物的穩(wěn)定性，須進行變形監(jiān)測[1-2]。傳統(tǒng)的變形監(jiān)測方法因其野外工作量大、受觀測環(huán)境影響、定位精度低以及測程有限等缺點，難以滿足長距離大范圍變形監(jiān)測的需求。合成孔徑雷達干涉技術(shù)(InSAR)因其全天時全天候連續(xù)觀測等優(yōu)點而得到迅速發(fā)展，成為變形監(jiān)測領(lǐng)域的新技術(shù)和研究熱點[3]。

GBInSAR技術(shù)由于測量范圍大、精度高、非接觸等優(yōu)點[4]，在變形監(jiān)測領(lǐng)域的應(yīng)用具有明顯的現(xiàn)實意義和實用價值。地基合成孔徑雷達(Ground-Based Synthetic Aperture Radar，GBSAR)將步進頻率連續(xù)波技術(shù)(Stepped-Frequency Continuous Wave)、雷達差分干涉技術(shù)和合成孔徑雷達技術(shù)(SAR)相結(jié)合，系統(tǒng)根據(jù)監(jiān)測需要選擇最佳觀測視角和相應(yīng)的觀測平臺，以非接觸的測量方式在安全的距離內(nèi)獲取被監(jiān)測危險區(qū)域的高精度變形數(shù)據(jù)[5]。目前，GBInSAR技術(shù)在滑坡[6-7]、冰川[8]和大壩[9]的變形監(jiān)測中得到廣泛應(yīng)用。

1 GBInSAR的基本原理

GBInSAR技術(shù)源于空基InSAR技術(shù)(Airborne-Based InSAR，ABInSAR)，它是ABInSAR在地面的實現(xiàn)。GBInSAR的原理包括合成孔徑雷達技術(shù)(SAR)和差分雷達干涉技術(shù)(D-InSAR)。

1.1 合成孔徑雷達技術(shù)

合成孔徑雷達技術(shù)[10-11]采用步進頻率連續(xù)波(SFCW)信號體制實現(xiàn)距離向高分辨率，同時保證雷達波的長距離傳輸。合成孔徑技術(shù)是通過雷達在地面上的固定軌道上滑動來實現(xiàn)的，方位向分辨率

(1)

式中:λ為雷達發(fā)射波長，L為雷達發(fā)射點相對于監(jiān)測目標(biāo)的相對移動距離，即合成后的天線孔徑。通過將合成孔徑雷達技術(shù)(SAR)和步進頻率連續(xù)波技術(shù)(SFCW)進行結(jié)合，可以在較大程度上同時提高距離向和方位向分辨率。

1.2 雷達差分干涉技術(shù)

假設(shè)GBSAR在不同時間獲取了同一目標(biāo)區(qū)域的兩幅SAR復(fù)影像，由于兩次采集圖像中間具有時間間隔，因此圖像之間因為回波信號的不同而產(chǎn)生了相應(yīng)的相位差。由解纏后的相位差可反演出目標(biāo)LOS方向的真實形變值：

(2)

由式(2)可知，目標(biāo)形變值的測量精度同時受系統(tǒng)工作頻率fc和相位測量精度的影響。

2 GBInSAR測量精度分析

GBInSAR應(yīng)用于變形監(jiān)測時的精度要達到毫米或亞毫米級，要提高變形監(jiān)測的精度就要掌握影響測量精度的因素。本文從GBInSAR的數(shù)據(jù)獲取和處理兩個過程著手，將影響GBInSAR測量精度的因素分為干涉相位誤差和形變真值解算誤差兩部分。

2.1 干涉相位誤差

干涉相位誤差是GBInSAR形變測量誤差的主要部分，受觀測期間大氣環(huán)境的改變、噪聲以及相位解纏誤差等因素的影響，

(3)

式中:φdisp為目標(biāo)形變造成的相位差，φatm為大氣擾動引起的相位差，φnoise為噪聲的影響，εφ為相位解纏誤差。

大氣效應(yīng)的影響是干涉相位誤差的重要來源，現(xiàn)在還沒有系統(tǒng)的解決方案。噪聲的影響在濾波和后續(xù)信號處理中會得到抑制。相位解纏是GBInSAR數(shù)據(jù)處理中的難點，當(dāng)進行短時間間隔或連續(xù)監(jiān)測時，由于觀測時間間隔短，獲得的相鄰圖像間相干性較好，一般不存在相位解纏問題。當(dāng)觀測時間間隔較長時，受時間、空間和噪聲相關(guān)等影響，難以獲得整個觀測區(qū)域精確的解纏相位。

2.2 形變真值解算誤差

形變真值是由視向形變值根據(jù)系統(tǒng)工作幾何關(guān)系計算出來的，觀測距離不同，采用的計算模型也有所區(qū)別。GBSAR系統(tǒng)的幾何關(guān)系有兩種假設(shè)模型[12]，分別是平行近似幾何關(guān)系和精確幾何關(guān)系，如圖1所示。

(a)平行近似幾何關(guān)系

(b)精確幾何關(guān)系圖1 形變真值假設(shè)模型

圖1中R1，R2為前后兩次觀測時觀測點到目標(biāo)的斜距長度；ΔR為視向形變值；Δr為形變真值；θ為雷達入射角。經(jīng)過推算可得兩種近似所產(chǎn)生的誤差為

(4)

當(dāng)取斜距R1=100 m，目標(biāo)形變值ΔR=0.1 m時，由式(4)可以計算出雷達入射角θ和形變真值近似計算誤差δr之間的關(guān)系，如圖2所示。

圖2 雷達入射角θ與形變真值近似計算誤差δr

通過圖2可以看出當(dāng)雷達入射角θ=30°時，近似計算誤差仍然可達亞毫米級，但當(dāng)θ趨于0，也即雷達入射波與目標(biāo)形變方向垂直時，近似誤差為無窮大，此時GBSAR系統(tǒng)失去形變監(jiān)測的能力。θ越接近90°時，計算誤差越小，在測量時可適當(dāng)增大雷達入射角，使得測量精度更高。

3 IBIS-L在隔河巖大壩變形監(jiān)測中的應(yīng)用

本次形變測量實驗以隔河巖大壩為例，實驗所用數(shù)據(jù)為從2013-07-30 T11 h 20 min至2013-08-02 T11 h 7 min觀測的數(shù)據(jù)，觀測時長2 d23 h47 min，每隔約5 min采集一幅圖像，共獲取SAR影像685幅。圖3展示實驗時IBIS-L儀器與大壩的相對位置，設(shè)備可達到的最大觀測距離為4 km，而實驗時所需觀測的最大距離約為1.3 km，完全滿足

實驗需求。實驗時設(shè)備安放在比較穩(wěn)定的地方，與觀測區(qū)域通視，中間無任何遮擋物，可直觀方便的觀測到整個壩體區(qū)域以及周邊坡岸。圖4為大壩信號反射強度圖。

圖3 IBIS-L設(shè)備裝置及觀測區(qū)域

圖4 大壩真實圖與信號反射強度

在圖4中用矩形框和橢圓形標(biāo)注出現(xiàn)場中物體在雷達信號反射圖中的位置。其中框代表壩體，橢圓代表旁邊山坡。通過圖4的對比可以看出IBIS-L可準(zhǔn)確的獲取壩體表面的雷達反射信息。

對獲取的影像進行干涉處理后，選取質(zhì)量較好的地面控制點，利用固定點法濾除大氣擾動誤差[13-14]，并進行中值濾波和控制點改正，得到大壩壩體觀測數(shù)據(jù)的累計視線向位移圖(如圖5所示)。從圖5中可以看出系統(tǒng)可以獲取監(jiān)測區(qū)域較高空間分辨率的變形信息，大壩壩體在觀測時間內(nèi)累計變形在2 mm之內(nèi)。為了對變形監(jiān)測結(jié)果進行驗證，在大壩壩體選擇了3個像元點進行時間序列分析。圖6為P1-P3雷達視線向形變時間序列。

圖5 累計視線向位移圖

圖6分別為所選的3個像元點在整個觀測過程中的位移時間關(guān)系圖。綜合圖5和圖6看出單點位

移隨時間呈現(xiàn)不規(guī)律的變化，在某些時刻波動較大，可能是受大氣延遲校正不完善的影響，但整體具有一定的規(guī)律性。從圖6看出像元點的變形值精度為亞毫米級，最終的變形保持在2 mm以內(nèi)，這與累計視線向位移圖結(jié)果一致。

(a)P1雷達視線向時間序列

(b)P2雷達視線向時間序列

(c)P3雷達視線向時間序列圖6 P1～P3雷達視線向時間序列

4 結(jié) 論

地基合成孔徑雷達采用步進頻率連續(xù)波信號體制實現(xiàn)距離向的高分辨率，同時保證雷達波的長距離傳輸，在獲取目標(biāo)微小形變方面具有很大的優(yōu)勢，現(xiàn)已在大壩、橋梁、邊坡等工程的變形監(jiān)測中得到廣泛應(yīng)用。通過對隔河巖大壩2 d 23 h 47 min的連續(xù)觀測，得到大壩觀測時間內(nèi)的累計視線向位移。實驗結(jié)果表明GBInSAR技術(shù)可獲得變形區(qū)域高空間密度的變形信息；GBInSAR技術(shù)全天時全天候的連續(xù)觀測能力，在大壩變形監(jiān)測方面具有廣闊的應(yīng)用前景。

[1] 邱志偉,岳建平,汪學(xué)琴. 地基雷達系統(tǒng)IBIS-L在大壩變形監(jiān)測中的應(yīng)用[J]. 長江科學(xué)院院報,2014(10):104-107.

[2] 劉學(xué)敏，黃其歡，田林亞. IBIS-L系統(tǒng)及其在大壩變形監(jiān)測中的應(yīng)用[J].測繪與空間地理信息，2015,38(7):34-36.

[3] 何儒云,王耀南,毛建旭. 合成孔徑雷達干涉測量(InSAR)關(guān)鍵技術(shù)研究[J]. 測繪工程,2007,16(5):53-56,60.

[4] 鄒進貴,田徑,陳艷華,等. 地基SAR與三維激光掃描數(shù)據(jù)融合方法研究[J]. 測繪地理信息,2015(3):26-30.

[5] 黃其歡,張理想. 基于GBInSAR技術(shù)的微變形監(jiān)測系統(tǒng)及其在大壩變形監(jiān)測中的應(yīng)用[J]. 水利水電科技進展,2011(3):84-87.

[6] LUZI G,NOFERINI L,MECATTI D, et al. Using a ground-based SAR interferometer and a terrestrial laser scanner to monitor a snow-covered slope: results from an experimental data collection in Tyrol(Austria)[J]. IEEE Trasactions on Geoscience and Remote Sensing,2009,47(2):382-393.

[7] HERRERAG,FERN NDEZ-MERODO J A,MULAS J, et al. A landslide forecasting model using ground based SAR data: the portalet case study[J].Engineering Geology,2009,105(3/4):220-230.

[8] LUZI G,PIERACCINI M,MECATTI D, et al. Monitoring of an alpine glacier by means of ground based SAR interferometry[J].Geoscience and Remote Sensing Letters,IEEE,2007,4 (3):495-499.

[9] MARIO A,GIULIA B,ALBERTO G. Measurement of dam deformations by terrestrial interferometric techniques[ C]// CHEN Jun, JIANG Jie, ALAIN B. The XXI Congress of the International Society for Photogrammetry and Remote Sensing. Beijing:ISPRS,2008:133-139.

[10] 侯建國,初禹. 合成孔徑雷達差分干涉測量技術(shù)在城市地面沉降監(jiān)測中的應(yīng)用[J]. 測繪工程,2014,23(8):40-44.

[11] 周校,王鵬,邢誠. 基于GB-SAR的建筑物微變形測量研究[J]. 測繪地理信息,2012(5):40-43.

[12] 曲世勃,王彥平,譚維賢,等.地基SAR形變監(jiān)測誤差分析與實驗[J].電子與信息學(xué)報,2011(1):1-7.

[13] 劉振宇，高雪峰，袁欣華，等.GPS技術(shù)在變形監(jiān)測中的應(yīng)用和發(fā)展趨勢[J].測繪與空間地理信息，2015,38(6):140-143.

[14] 董杰,董妍. 基于氣象數(shù)據(jù)的地基雷達大氣擾動校正方法研究[J]. 測繪工程,2014,23(10):72-75.

[責(zé)任編輯：李銘娜]

Application research of GBInSAR in dam deformation monitoring

GUO Leping, YUE Jianping, YUE Shun, QIU Zhiwei, WANG Xueqin

(School of Earth Sciences and Engineering, Hohai University, Nanjing 211100, China)

In order to evaluate the working behaviour of hydraulic structure, the deformation and stress information concerning its high precision are needed. In this paper, according to the characteristics and requirements of the dam safety monitoring, the constitution and the basic principle of measurement of the GBSAR monitoring system have been introduced. The measurement error sources and characteristics of the system have been analyzed. The feasibility of GBInSAR technology used in dam deformation monitoring has been discussed. By taking the measured data during three days of Geheyan dam as an example,it analyzes the actual measurement precision of the technology can reach millimeter level, and the technology can meet the accuracy requirement of the dam deformation monitoring.

deformation monitoring; GBInSAR; dam; LOS displacement; time series

10.19349/j.cnki.issn1006-7949.2017.02.014

2015-11-24

國家自然科學(xué)基金資助項目(41174002)；湖南省重點實驗室開放研究基金資助項目(PKLHD201311)；中央高?；究蒲袠I(yè)務(wù)費項目(2013/B14020383) .

郭樂萍(1991-)，女，碩士研究生.

P231

A

1006-7949(2017)02-0062-04

引用著錄：郭樂萍，岳建平，岳順，等.GBInSAR在大壩變形監(jiān)測中的應(yīng)用[J].測繪工程，2017，26(2)：62-65.