顏芳 孫芳 黃鑫 徐銳 岳雷



摘? ?要:GNSS接收機因其高精度的點位監測在地表位移形變監測領域的市場逐年增加。目前,已成為研究熱點的InSAR技術在地面形變監測方面具有高空間分辨率和連續覆蓋的特點。文章將這兩種技術進行組合設計,重點描述了組合系統的組成、原理、現場設備的設計、云平臺處理軟件設計,并在涿州園區進行試驗驗證,后期可推廣應用于形變監測領域。
關鍵詞:全球導航衛星系統;接收機;合成孔徑雷達干涉測量;
1? ? 項目背景
監測各類地質災害現象(如滑坡、地面沉降、采空塌陷)所引發的三維地表位移,對于準確認識地質災害現象的發生機制、運動過程和發展趨勢,具有重要的意義,近年來各類地質災害監測項目中越來越多地采用了全球導航衛星系統(Global Navigation Satellite System,GNSS)接收機進行三維地表位移的監測。但GNSS監測存在面定位精度不高的問題,將結合合成孔徑雷達干涉測量技術(Interferometric Synthetic Aperture Radar,InSAR)[1-2]高空間分辨率和連續覆蓋的特點,形成點面結合的三維地表位移監測系統,既突破單一監測技術的應用局限,又極大地提高了三維形變監測結果的時空分辨率[3]。
2? ? 總體方案
2.1? 系統組成
GNSS和InSAR組合監測系統由現場監測設備和云平臺數據處理軟件組成,現場監測設備安裝在災害監測點附近,主要包括GNSS接收機和角反射器。云平臺數據處理軟件可對GNSS接收機和角反射器的數據進行處理。其組成框架如圖1所示。
2.2? 原理介紹
GNSS接收機用來測量地表的三維形變,在一個變形區域內布設一個GNSS基準站和若干個GNSS觀測站,接收視場內的所有衛星數據(包括GPS、北斗、GLONASS),并通過4G無線通信網絡將采集的衛星觀測數據上傳至云平臺。
角反射器是InSAR中的輔助測量工具,可增強星載雷達的反射,作為監測區域在InSAR影像中的“地面控制點”[4],用于精確配準InSAR影像[4],以提高區域形變監測的準確性。角反射器最好布設在GNSS觀測站5 m以內的范圍。人工按照衛星的觀測周期(一般為7~14天),從云端下載衛星數據,然后通過InSAR數據處理軟件進行處理。
GNSS數據處理軟件負責接收GNSS接收機的衛星觀測數據,進行解算得到觀測站的坐標值,并與初始測量值進行對比,得到觀測站X,Y,Z 3個方向上的位移變化值。
InSAR數據處理軟件利用GNSS接收機得到的三維坐標值對衛星InSAR影像進行精確配準,然后進行解譯處理,得到角反射器安裝區域的高程變化。原理架構如圖2所示。
3? ? 現場設備設計
3.1? GNSS接收機設計
GNSS接收機主要由太陽能充電管理電路、電源管理系統、MCU主控模塊、GNSS板卡、GNSS天線、4G模塊、數據存儲模塊和調試接口等組成。其組成框架如圖3所示。
為了提高接收機的處理速度并降低設備的功耗,最終準備采用ST公司的STM32L系列單片機實現,它是超低功耗高性能微控制單元(Microcontroller Unit,MCU),采用Cortex-M4內核,工作頻率為80 MHz,有3個串行外設接口(Serial Peripheral Interface,SPI)接口、4個通用異步收發傳輸器(Universal Asynchronous Receiver/Transmitter,UART)接口和1個低功耗的UART接口;帶有日歷功能的32位RTC小于1 μA的實時時鐘,1 s精度;多路12位模擬/數字轉換器(Analog to Digital Converter,ADC)通道等,主要用于各種低功耗應用如物聯網(Internet of Things,IoT)等。
(1)太陽能充電管理電路可對太陽能電池的最大功率點進行跟蹤,自動調整充電電流,對深度放電的電池進行涓流充電,具備充電狀態、充電結束狀態指示和電池溫度檢測功能。
(2)電源管理單元主要用于將電池電壓轉換為各模塊需要的工作電壓,并通過單片機控制各模塊的供電,以降低系統功耗。
(3)GNSS板卡選用三系統全頻板卡,靜態精度:水平不大于(±2.5±1)×10-6,垂直不大于(±5±1)×10-6;采樣頻率可根據現場的要求進行設定,默認值設為10 s。
(4)GNSS天線為全頻天線,采用抗多路徑干擾設計。
(5)數據存儲模塊采用SD卡存儲的方式,可根據客戶的要求配置SD卡的容量。
(6)調試接口采用RS232接口,方便維護。
(7)4G模塊選用移遠的EC20 R2.1 Mini PCIe來實現。該模塊支持標準的AT指令集3GPP TS 27.007,27.005及移遠通信增強型AT命令;其工作電壓:3.0~3.6 V;支持USIM和SIM:1.8 V和3.3 V;尺寸:(51.0±0.15)mm×(30.0±0.15)mm×(4.9±0.2)mm;正常工作溫度:﹣35~﹢75 ℃;擴展工作溫度:﹣40~﹢80 ℃。該模塊在網絡注冊時的功耗約為70 mA,網絡注冊成功但無數據交互的功耗為15~19 mA,睡眠階段的功耗約為3 mA,發送數據階段的功耗約為100 mA。
3.2? 角反射器設計
本方案中擬選用三面角的角反射器[6-7],材料為鋁板和鍍鋅鐵皮雙層結構,鋁板厚3 mm,外加鍍鋅鐵皮(1 mm厚)以保護反射面(鋁板),邊側加三角角鋼加固;3塊金屬板之間的相互垂直角度加工公差不超過±1;在棱邊設置了3個活動關節,通過伸縮桿來調節人工角反射器(Corner Reflector,CR)的仰角;在頂底處設置了一個漏水孔,使積水不至于影響其反射路線;為獲得雷達散射截面(Radar Cross Section,RCS),可根據雷達數據軌道信息來調整角反射器的底邊方位角,并使角反射器的底邊與衛星飛行方向平行。角反射器的示意和實物如圖4—5所示。
4? ? 云平臺數據處理軟件設計
云平臺數據處理軟件包括GNSS數據處理軟件和InSAR數據處理軟件,GNSS數據處理軟件可對GNSS接收機的數據進行靜態實時解算,并按照設定的周期輸出高精度的三維坐標值。InSAR數據處理軟件采用sarmap公司的SARscape軟件,它是國際知名的雷達圖像處理軟件。該軟件架構于專業的ENVI遙感圖像處理軟件之上,提供完整的SAR數據處理功能,全面支持雷達干涉測量(InSAR/DInSAR)數據處理和分析,輸出SAR圖像產品、數字高程模型(Digital Elevation Model,DEM)和地表形變圖等信息。
本系統將GNSS數據處理軟件解算的GNSS觀測站的高精度三維坐標值(X,Y,H)帶入到SARscape軟件中,不僅能對InSAR的影像進行插值,還能快速識別該GNSS觀測站附近的角反射器,并利用角反射器的三維坐標對衛星的干涉基線進行重新估算,從而提高觀測精度。
5? ? 試驗驗證
5.1? 試驗設備
在涿州工業園區安裝了兩套GNSS接收機和一套CR,兩套GNSS接收機分別配置為基準站和觀測站。在靠近GNSS觀測站5 m的距離范圍內布設了一臺反射器,角反射器針對哨兵1號衛星布設。設備如圖6所示。
5.2? 試驗結果
5.2.1? GNSS觀測站的3個方向的變化曲線
GNSS觀測站和基準站的原始數據采集周期設定為10 s,星歷周期為60 s,數據解算周期為1 h。為驗證GNSS是否能測試到高程的變化,設計了一個15 mm的墊塊,在2018年7月3日將GNSS觀測站的高程調高了15 mm。圖7—9選取了GNSS觀測站從2018年6月24日到2018年7月5號的數據,x,y,h方向分別為高斯投影北、東、高3個方向的坐標值,測試結果如下。
5.2.2? 通過InSAR監測涿州區域的形變曲線
本次試驗獲取了從2018年1月—11月的20景哨兵影像。可能是CR的設計和安裝問題,導致其在SAR影像上的十字星特征不明顯。
經過IN-SAR形變計算,獲得監測區域的平均形變如圖10所示,最小變化區間為-10.569~-2.221,最大形變區間在-89.877~-60.062,標記為2的點代表角反射器安裝的區域,其變形區間在-33.228 633 65~-21.302 676 48。可以看出,大部分形變變化在-10.569~-2.221之間。
5.2.3? 數據分析
可以看出,GNSS觀測站的x方向最大變化值約為4 mm,y方向最大變化值約為5 mm,h方向在2018年7月3號調整了15 mm,未調整高度前h方向的最大變化值為5 mm,可以明顯地看出調整前后數據有15 mm的形變趨勢。圖10中標記為2的點代表角反射器的安裝位置,因角反射器設計和安裝的誤差,從SAR影像中未能獲取到角反射器明顯的十字星,故通過SAR影像獲取的涿州監測區域的平均形變也是存在一定誤差的。
6? ? 結語
本系統簡要描述了一套集GNSS技術和InSAR技術的組合監測系統,通過在涿州園區進行測試,驗證了GNSS技術的mm級形變監測能力。盡管角反射器的設計和安裝位置不夠理想,給后期SAR影像的處理帶來了一定的困難,但也從側面反映了該組合系統不僅可以監測一些高精度點位上的形變量,而且可以測試一定區域的形變。考慮到InSAR技術對角反射器的設計、安裝、選點及后期數據處理有很高的要求,要進行大面積推廣,需有專業的機構進行指導。后期會對這套裝置進行設計完善,確保這套點面測量相結合的監測系統盡快推向市場。
[參考文獻]
[1]李德仁.衛星雷達干涉測量原理與應用[J].測繪科學,2000(1):9-12.
[2]何敏.合成孔徑雷達干涉測量技術及其在形變災害監測中的應用[J].水電自動化與大壩監測,2005(2):45-48.
[3]金雙根,朱文耀.GPS觀測數據提高InSAR干涉測量精度的分析[J].遙感信息,2001(4):8-11.
[4]諶華.CR-InSAR大氣校正模型研究及其初步應用[D].北京:中國地震局地質研究所,2006.
[5]諶華,甘衛軍.利用GPS與InSAR融合提高形變監測精度方法研究[J].大地測量與地球動力學,2010(3):59-62.
[6]姜山.三角形三面角反射器加工公差對其單站RCS影響研究[J].航空兵器,2006(4):24-27.
[7]袁禮海.SAR輻射定標精度設計與分析[J].雷達科學與技術,2009(1):35-39.
GNSS and InSAR integrated monitoring system design
Yan Fang, Sun Fang, Huang Xin, Xu Rui, Yue Lei
(Aerospace Science & Industry Inertial Technology Co., Ltd., Beijing 100070, China)
Abstract:The market of GNSS receiver in the field of surface displacement deformation monitoring is increasing year by year because of its high precision point monitoring. At present, the InSAR technology, which has become a research hotspot, has the characteristics of high spatial resolution and continuous coverage in ground deformation monitoring. In this paper, the combination design of these two technologies is carried out, focusing on the composition, principle, field equipment design and cloud platform processing software design, and the test verification is carried out in Zhuozhou Park, which can be popularized and applied in the field of deformation monitoring in the later stage.
Key words:global navigation satellite system; receiver; interferometric synthetic aperture radar; angle reflector