劉斌， 葛大慶， 李曼， 張玲， 王艷， 郭小方， 張曉博

(1.中國(guó)國(guó)土資源航空物探遙感中心，北京 100083；2.中國(guó)礦業(yè)大學(xué)(北京)地球科學(xué)與測(cè)繪工程學(xué)院，北京 100083)

地基合成孔徑雷達(dá)干涉測(cè)量技術(shù)及其應(yīng)用

劉斌1， 葛大慶1， 李曼1， 張玲1， 王艷1， 郭小方1， 張曉博2

(1.中國(guó)國(guó)土資源航空物探遙感中心，北京 100083；2.中國(guó)礦業(yè)大學(xué)(北京)地球科學(xué)與測(cè)繪工程學(xué)院，北京 100083)

在簡(jiǎn)要討論地基合成孔徑雷達(dá)干涉測(cè)量(ground-based interferometric synthetic aperture radar，GB InSAR)技術(shù)原理的基礎(chǔ)上，歸納了常規(guī)地基InSAR形變測(cè)量的主要處理流程。全面總結(jié)了地基干涉測(cè)量雷達(dá)系統(tǒng)的主要類(lèi)型及其發(fā)展趨勢(shì)，地基InSAR技術(shù)的主要應(yīng)用領(lǐng)域以及目前存在的問(wèn)題，對(duì)比分析了地基和星載InSAR在形變測(cè)量上的優(yōu)勢(shì)和不足。

地基合成孔徑雷達(dá)干涉測(cè)量； 步進(jìn)頻率連續(xù)波； 調(diào)頻連續(xù)波； 噪聲雷達(dá)； 多入多出技術(shù)

0 引言

星載干涉測(cè)量雷達(dá)系統(tǒng)的優(yōu)勢(shì)為覆蓋范圍廣、監(jiān)測(cè)精度可靠，能以較高頻率和密度獲取觀(guān)測(cè)區(qū)的變形狀況，適用于具有一定空間尺度的形變場(chǎng)連續(xù)監(jiān)測(cè)。針對(duì)地面多類(lèi)型、多尺度形變測(cè)量的監(jiān)測(cè)技術(shù)需求，星載雷達(dá)系統(tǒng)仍然存在如下缺點(diǎn)： ①觀(guān)測(cè)實(shí)時(shí)性受制于衛(wèi)星重訪(fǎng)周期，難以滿(mǎn)足動(dòng)態(tài)和應(yīng)急監(jiān)測(cè)需要； ②雷達(dá)入射角由于衛(wèi)星飛行方向固定(升軌或降軌)，且單一，邊坡觀(guān)測(cè)存在疊掩、陰影和頂?shù)椎怪玫痊F(xiàn)象，難以滿(mǎn)足邊坡監(jiān)測(cè)等觀(guān)測(cè)需要； ③進(jìn)行山體滑坡、礦山邊坡失穩(wěn)、冰川運(yùn)動(dòng)等局部變形監(jiān)測(cè)時(shí)，星載干涉測(cè)量雷達(dá)系統(tǒng)的空間分辨率、覆蓋范圍往往難以達(dá)到最佳匹配。

地基干涉測(cè)量雷達(dá)系統(tǒng)是利用布設(shè)在監(jiān)測(cè)對(duì)象附近的地面雷達(dá)測(cè)量系統(tǒng)，實(shí)現(xiàn)對(duì)被觀(guān)測(cè)對(duì)象的近實(shí)時(shí)動(dòng)態(tài)連續(xù)監(jiān)測(cè)。由于采用的雷達(dá)波長(zhǎng)更短，因此測(cè)量精度能達(dá)到毫米―亞毫米級(jí)，能夠根據(jù)觀(guān)測(cè)目標(biāo)形變場(chǎng)的演化特征靈活安置，彌補(bǔ)星載雷達(dá)重訪(fǎng)周期固定、入射角度單一等不足。

本文在詳細(xì)討論地基合成孔徑雷達(dá)干涉測(cè)量(ground-based interferometric synthetic aperture radar，GB InSAR)技術(shù)原理的基礎(chǔ)上，歸納常規(guī)GB InSAR形變測(cè)量的主要處理流程，并通過(guò)總結(jié)地基干涉測(cè)量雷達(dá)系統(tǒng)的發(fā)展現(xiàn)狀，以及GB InSAR技術(shù)的應(yīng)用領(lǐng)域，對(duì)比分析地基和星載InSAR在形變測(cè)量上的優(yōu)勢(shì)和不足，全面掌握地基干涉測(cè)量雷達(dá)系統(tǒng)的實(shí)用能力。

1 GB InSAR技術(shù)原理

一般來(lái)說(shuō)，地基雷達(dá)觀(guān)測(cè)期間的空間基線(xiàn)為0，因此地形相位為0。為了確定地形相位，2次掃描過(guò)程中，垂直移動(dòng)傳感器的位置可以產(chǎn)生空間基線(xiàn)。如圖1(據(jù)文獻(xiàn)[1]修改)所示，輕微平移滑軌，使得2次觀(guān)測(cè)過(guò)程中雷達(dá)傳感器的位置發(fā)生變化(M1→M2)，將2次觀(guān)測(cè)獲得的復(fù)數(shù)影像共軛相乘，那么得到的干涉相位中就會(huì)包含觀(guān)測(cè)區(qū)域的地形相位。觀(guān)測(cè)相位與地形的關(guān)系可以表示為

圖1 GB InSAR示意圖(據(jù)文獻(xiàn)[1]修改)

(1)

根據(jù)圖1中的幾何關(guān)系，式(1)可以寫(xiě)為

(2)

式中：R1和R2分別為第1次和第2次測(cè)量距離；λ為波長(zhǎng)；b為M1和M2之間的距離(空間基線(xiàn))；h為目標(biāo)點(diǎn)的高程。假定(b/R)<<1，且R1-R2<

(3)

(4)

按照基線(xiàn)的類(lèi)型，干涉相位△φ21主要分為4部分，即

(5)

式中：φtopo為空間基線(xiàn)引起的地形相位；φdisp和φatm為形變和大氣影響產(chǎn)生的時(shí)間相位變化；φnoise為噪聲；n為整周模糊數(shù)。式(5)左側(cè)為觀(guān)測(cè)干涉相位，而右側(cè)含有未知參數(shù)。為了求解地形相位，式(5)可寫(xiě)成

φtopo=△φ21-φdisp-φatm-φnoise+2πn。

(6)

2次掃描期間垂直移動(dòng)傳感器的位置可以產(chǎn)生空間基線(xiàn)Bs，由于存在時(shí)間基線(xiàn)Bt，這樣時(shí)間基線(xiàn)相關(guān)分量仍是方程的一部分。采用2個(gè)具有一定距離的垂向分布的天線(xiàn)同時(shí)接收信號(hào)，可以大大簡(jiǎn)化方程的計(jì)算，然后估計(jì)函數(shù)模型中的時(shí)間基線(xiàn)相關(guān)分量φdisp和φatm。

對(duì)于地基雷達(dá)形變監(jiān)測(cè)，假定空間基線(xiàn)為0，φdisp可以寫(xiě)為

φdisp=△φ21-φatm-φnoise+2πn。

(7)

2 常規(guī)GB InSAR形變測(cè)量處理流程

GBInSAR形變測(cè)量數(shù)據(jù)處理與星載InSAR類(lèi)似，但其處理流程相對(duì)更為復(fù)雜。GBInSAR形變測(cè)量處理主要流程如圖2所示(據(jù)文獻(xiàn)[2]修改)。

圖2 地基InSAR形變測(cè)量數(shù)據(jù)處理流程(據(jù)文獻(xiàn)[2]修改)

Fig.2 Displacement process flow of GB InSAR

1)GB InSAR影像配準(zhǔn)。影像配準(zhǔn)是GB InSAR形變測(cè)量的關(guān)鍵步驟之一，是將2景影像中代表相同地物的像元匹配到同一位置。影像配準(zhǔn)的步驟主要分為影像匹配、仿射變換和影像重采樣等3步。如果影像間的像元位置沒(méi)有發(fā)生較大變動(dòng)，一般來(lái)說(shuō)，配準(zhǔn)過(guò)程能夠高度自動(dòng)化完成。地基雷達(dá)對(duì)同一目標(biāo)區(qū)域進(jìn)行重復(fù)觀(guān)測(cè)時(shí)，由于觀(guān)測(cè)平臺(tái)的微小偏移會(huì)使地基雷達(dá)觀(guān)測(cè)軌道和觀(guān)測(cè)視角發(fā)生微小改變，造成影像在距離向和方位向發(fā)生一定的錯(cuò)位和扭曲，因此，非連續(xù)模式觀(guān)測(cè)必須進(jìn)行影像配準(zhǔn)。而連續(xù)模式觀(guān)測(cè)時(shí)由于軌道固定不動(dòng)，GB InSAR通常無(wú)需進(jìn)行影像配準(zhǔn)。但對(duì)于長(zhǎng)距離觀(guān)測(cè)并且觀(guān)測(cè)時(shí)間較長(zhǎng)的情況，由于某些影像獲取期間大氣的變化使得影像發(fā)生畸變，因此部分影像也應(yīng)進(jìn)行校正補(bǔ)償[3]。

(8)

由于連續(xù)模式觀(guān)測(cè)的空間基線(xiàn)為0，因此無(wú)需進(jìn)行地形相位補(bǔ)償； 非連續(xù)模式觀(guān)測(cè)的空間基線(xiàn)不為0，引入額外的地形相位時(shí)需要將干涉圖中的地形相位消除。

3)相位解纏。相位解纏是干涉相位主值恢復(fù)至真實(shí)相位值的過(guò)程，是GBInSAR形變測(cè)量的重要環(huán)節(jié)之一。目前的相位解纏算法通常需滿(mǎn)足

(9)

4)大氣相位校正。解纏相位中包括形變相位、大氣相位和噪聲。噪聲可以通過(guò)空間低通濾波削弱，大氣相位主要有3種方法予以消除： ①利用氣象觀(guān)測(cè)數(shù)據(jù)(溫度、濕度和壓力)確定大氣相位； ②通過(guò)穩(wěn)定點(diǎn)(地面控制點(diǎn))估計(jì)大氣相位； ③空間插值或?yàn)V波。

5)地理編碼。為了解譯觀(guān)測(cè)數(shù)據(jù)以及與其他觀(guān)測(cè)數(shù)據(jù)進(jìn)行對(duì)比分析，需要將雷達(dá)坐標(biāo)系中的結(jié)果投影到地理坐標(biāo)系下，即獲得主影像每一像元的三維位置矢量。

3 地基干涉測(cè)量雷達(dá)系統(tǒng)發(fā)展和應(yīng)用

3.1 發(fā)展現(xiàn)狀

地基干涉測(cè)量雷達(dá)系統(tǒng)目前多采用Ku波段雷達(dá)波，主要包括雷達(dá)傳感器、能量供應(yīng)單元、數(shù)據(jù)采集和存儲(chǔ)單元以及數(shù)據(jù)處理模塊等。按照成像方式的不同可分為2大類(lèi)： 真實(shí)孔徑雷達(dá)(real aperture radar，RAR)和合成孔徑雷達(dá)(synthetic aperture radar，SAR)。商用型地基干涉測(cè)量雷達(dá)系統(tǒng)主要有： IBIS系列地基雷達(dá)系統(tǒng)(IDS，意大利)[4]、GPRI便攜式雷達(dá)干涉儀(GAMMA，瑞士)[5]、FastGBSAR系統(tǒng)(MetaSensing，荷蘭)[6]、ARAMIS系統(tǒng)(ARELECS，意大利)[7]、SSR系統(tǒng)(GroundProbe，澳大利亞)[8]和S-SAR系統(tǒng)(中國(guó)安全生產(chǎn)科學(xué)研究院)[9]； 典型的科研型包括： LiSA[10]和MELISSA[11]系統(tǒng)(JRC，歐盟委員會(huì))、RiskSAR系統(tǒng)[12](UPC，西班牙)、GB NW-SAR系統(tǒng)[13](IRPHE，烏克蘭)、ARCSAR系統(tǒng)[14](KNU，韓國(guó))、ASTRO系統(tǒng)[15](中國(guó)科學(xué)院電子學(xué)研究所，中國(guó))和超寬帶軌道SAR系統(tǒng)[16](國(guó)防科技大學(xué)，中國(guó))。其中，GPRI便攜式雷達(dá)干涉儀和SSR地基干涉測(cè)量雷達(dá)系統(tǒng)屬于真實(shí)孔徑雷達(dá)范疇。

按照雷達(dá)的工作體制劃分主要包括： 步進(jìn)頻率連續(xù)波(stepped frequency continuous wave，SFCW)[4,9,10,12,16]、調(diào)頻連續(xù)波(frequency modulation continuous wave，F(xiàn)MCW)[5-7]、噪聲雷達(dá)[13-14]以及多入多出技術(shù)(multiple input multiple output，MIMO)[11,15]等雷達(dá)系統(tǒng)。地基干涉測(cè)量雷達(dá)系統(tǒng)研制初期，由于矢量網(wǎng)絡(luò)分析儀(vectorial network analyzer，VNA)無(wú)需復(fù)雜的電子元器件，便可產(chǎn)生類(lèi)似于SFCW的電磁波，一些地基雷達(dá)傳感器(如IBIS和LiSA)的核心通常采用VNA產(chǎn)生帶通信號(hào)、檢波和模擬數(shù)字轉(zhuǎn)換?；赩NA的地基雷達(dá)系統(tǒng)是以連續(xù)波的形式發(fā)射離散的頻率分量，然后在時(shí)間域內(nèi)重建電磁脈沖的波形合成一個(gè)等效脈沖，因此單景影像的采集時(shí)間較長(zhǎng)(最短幾min)。較長(zhǎng)的掃描時(shí)間增加了大氣擾動(dòng)的可能性，同時(shí)由于掃描期間觀(guān)測(cè)目標(biāo)的不穩(wěn)定也可能會(huì)導(dǎo)致影像幅度和相位的失真[17]。近幾年來(lái)，科研人員開(kāi)始利用FMCW雷達(dá)(如FastGBSAR、ARAMIS)取代VNA，F(xiàn)MCW雷達(dá)，這不但可以提高設(shè)備的可操作性，而且能夠減少數(shù)據(jù)采集時(shí)間(幾s—幾十s)，降低觀(guān)測(cè)成本。此外，研究人員還將噪聲雷達(dá)(如GB NW-SAR)和多入多出技術(shù)(如MELISSA)應(yīng)用于地基雷達(dá)系統(tǒng)以獲得更快的數(shù)據(jù)采集時(shí)間，如MELISSA系統(tǒng)最短采集時(shí)間僅為0.26 s。相對(duì)于真實(shí)孔徑雷達(dá)，合成孔徑雷達(dá)通過(guò)沿較長(zhǎng)的軌道運(yùn)動(dòng)以獲得較大的孔徑，目前除了幾種最新的系統(tǒng)采用弧形掃描(如GB NW-SAR)以及無(wú)運(yùn)動(dòng)掃描(如MELISSA)外，雷達(dá)傳感器沿線(xiàn)性軌道滑動(dòng)掃描是最常用的掃描方式。

3.2 應(yīng)用領(lǐng)域

GB InSAR技術(shù)按照數(shù)據(jù)獲取方式的不同分為連續(xù)模式(continuous ground-based interferometric synthetic aperture radar，C-GB InSAR)和非連續(xù)模式(dis-continuous ground-based interferometric synthetic aperture radar，D-GB InSAR)2類(lèi)[18]。連續(xù)模式觀(guān)測(cè)通常將設(shè)備安裝在觀(guān)測(cè)現(xiàn)場(chǎng)，根據(jù)需要合理配置雷達(dá)參數(shù)，長(zhǎng)時(shí)間持續(xù)采集數(shù)據(jù)(如每隔幾分鐘采集一次)，這種模式允許用戶(hù)采用“近似實(shí)時(shí)”監(jiān)測(cè)目標(biāo)區(qū)，適宜觀(guān)測(cè)快速形變現(xiàn)象(如每天幾cm)。由于設(shè)備安裝在固定觀(guān)測(cè)站，因此空間基線(xiàn)和地形相位皆為0。非連續(xù)模式觀(guān)測(cè)通常需在觀(guān)測(cè)區(qū)建立固定觀(guān)測(cè)站，根據(jù)形變的運(yùn)動(dòng)學(xué)特征設(shè)置一定的觀(guān)測(cè)周期，定期觀(guān)測(cè)獲取數(shù)據(jù)(如以周、月或年為周期)，非連續(xù)模式一般用來(lái)觀(guān)測(cè)緩慢變形現(xiàn)象。連續(xù)模式和非連續(xù)模式2種類(lèi)型的GB InSAR數(shù)據(jù)處理策略基本相同，但非連續(xù)模式往往面臨2個(gè)難點(diǎn)： 時(shí)間失相干和相位混疊現(xiàn)象，觀(guān)測(cè)點(diǎn)密度、監(jiān)測(cè)精度以及形變測(cè)量的可靠性有所降低。

目前GB InSAR已在多個(gè)領(lǐng)域得到了廣泛應(yīng)用，如露天采礦邊坡[19]、滑坡[20]、崩塌[21]、火山[22]、單體建(構(gòu))筑物[23]、水壩/堤壩[24]、冰川[25]等地物的變形測(cè)量。其中，GB InSAR在露天采礦邊坡穩(wěn)定性、水壩/堤壩形變、單體建(構(gòu))筑物結(jié)構(gòu)變形等領(lǐng)域的應(yīng)用最為廣泛和成熟。由于這些目標(biāo)物適宜于雷達(dá)波反射，GB InSAR獲得的相干目標(biāo)數(shù)量多，其信噪比、相干系數(shù)和相位穩(wěn)定性都較高，能夠探測(cè)到觀(guān)測(cè)目標(biāo)的局部微小變形，為工程施工提供早期預(yù)警。對(duì)于單體建(構(gòu))筑物的結(jié)構(gòu)變形監(jiān)測(cè)，關(guān)鍵是對(duì)每一觀(guān)測(cè)像元進(jìn)行精確的地理編碼，而地理編碼精度取決于所用的DEM質(zhì)量。因此，通常采用激光雷達(dá)獲取精確的地形數(shù)據(jù)，或者利用GB InSAR得到同一區(qū)域的DEM，但后者得到的DEM一般來(lái)說(shuō)遠(yuǎn)不如激光雷達(dá)結(jié)果可靠。與露天采礦邊坡監(jiān)測(cè)相比，崩塌、山體滑坡和火山等地質(zhì)災(zāi)害涉及到邊坡穩(wěn)定性監(jiān)測(cè)，多數(shù)情況下觀(guān)測(cè)區(qū)域植被較為茂密，不利于雷達(dá)回波信號(hào)的接收，這些觀(guān)測(cè)區(qū)域的相干目標(biāo)相對(duì)稀少，有的區(qū)域甚至無(wú)法獲得相干目標(biāo)，但提取的相干目標(biāo)的測(cè)量精度依然可以達(dá)到毫米―亞毫米級(jí)。GB InSAR遠(yuǎn)距離、高時(shí)空采樣能力，也為冰川運(yùn)動(dòng)提供了一種新型監(jiān)測(cè)工具，利用GB InSAR觀(guān)測(cè)冰川的應(yīng)用限制主要?dú)w結(jié)于失相干問(wèn)題： 一方面由于冰川存在積雪覆蓋，不利于雷達(dá)波反射，觀(guān)測(cè)期間如遇降雪也會(huì)影響觀(guān)測(cè)效果； 另一方面對(duì)于SFCW 地基雷達(dá)系統(tǒng)的單景數(shù)據(jù)采集時(shí)間較長(zhǎng)，相對(duì)于冰川局部的移動(dòng)速度，幾min采集時(shí)間內(nèi)冰川局部已發(fā)生明顯變化。因此，對(duì)于移動(dòng)較快的冰川監(jiān)測(cè)適宜采用連續(xù)模式觀(guān)測(cè)。

從數(shù)據(jù)處理層面上講，大多數(shù)GB InSAR系統(tǒng)采用Ku波段雷達(dá)系統(tǒng)，這類(lèi)系統(tǒng)的硬件性能易于實(shí)現(xiàn)，可以獲得高空間分辨率影像，而且Ku波段對(duì)地物的微小變化比較敏感。但溫度、濕度等外界觀(guān)測(cè)條件的變化都可能會(huì)引起觀(guān)測(cè)目標(biāo)的微小形變，Ku波段的這種高形變敏感度也會(huì)帶來(lái)數(shù)據(jù)處理上的額外負(fù)擔(dān)，在數(shù)據(jù)處理過(guò)程中，如何消除環(huán)境變化引起的目標(biāo)微小形變的影響至關(guān)重要[16]。GB InSAR在2次數(shù)據(jù)采集間隔大氣條件的差異會(huì)產(chǎn)生額外的相位貢獻(xiàn)也不容忽視，研究人員已提出多種大氣模型改進(jìn)數(shù)據(jù)處理流程[12,16-17]，然而這些大氣校正模型均是基于均一大氣介質(zhì)模型的假設(shè)，對(duì)于非均一大氣的影響，目前還沒(méi)有一種可靠的方案消除大氣相位影響。另外，波長(zhǎng)越短失相干越嚴(yán)重，影像失相干可能會(huì)使相位產(chǎn)生混疊現(xiàn)象，增加了相位解纏的難度。在觀(guān)測(cè)區(qū)域的某些部位，采用連續(xù)觀(guān)測(cè)模式觀(guān)測(cè)失相干現(xiàn)象相對(duì)容易控制，如果采用非連續(xù)觀(guān)測(cè)模式，這種失相干現(xiàn)象更為明顯。為了提高相干性，一種較為有效的方法就是采用較短的滑軌，但這種做法勢(shì)必會(huì)降低觀(guān)測(cè)精度和空間分辨率。

3.3 地基與星載 InSAR形變測(cè)量比較

1)零空間基線(xiàn)。地基雷達(dá)系統(tǒng)可以完全控制空間基線(xiàn)。對(duì)于連續(xù)觀(guān)測(cè)模式形變監(jiān)測(cè)來(lái)說(shuō)，空間基線(xiàn)為0，零空間基線(xiàn)干涉測(cè)量無(wú)需DEM參與，因此它是最理想的干涉測(cè)量形變獲取方式。

2)軌道控制。由于觀(guān)測(cè)設(shè)備位置固定，因此干涉測(cè)量處理沒(méi)有軌道誤差。

3)時(shí)間采樣率。地基雷達(dá)系統(tǒng)的最短數(shù)據(jù)采集間隔在亞s到幾min之間，而星載雷達(dá)的重訪(fǎng)周期通常為幾d—幾十d。高時(shí)間采樣率在很大程度上簡(jiǎn)化了相位解纏，甚至無(wú)需進(jìn)行相位解纏。

4)入射角。根據(jù)監(jiān)測(cè)目標(biāo)的形變特征差異，星載雷達(dá)和地基雷達(dá)各有自己的優(yōu)勢(shì)： 對(duì)于地面沉降觀(guān)測(cè)，星載雷達(dá)具有最佳視角； 而地基雷達(dá)更適合陡坡形變觀(guān)測(cè)。

5)觀(guān)測(cè)范圍。星載雷達(dá)可以監(jiān)測(cè)大范圍區(qū)域，而地基雷達(dá)只能監(jiān)測(cè)幾km2的區(qū)域。

6)測(cè)量精度。地基和星載InSAR的觀(guān)測(cè)精度都取決于目標(biāo)的散射特征以及與參考點(diǎn)的距離。GB InSAR的觀(guān)測(cè)精度往往與設(shè)備和目標(biāo)間的距離有關(guān)，由于波長(zhǎng)更短(如Ku波段)，監(jiān)測(cè)精度一般在亞mm到幾個(gè)mm之間。

4 結(jié)論

本文簡(jiǎn)要討論了地基合成孔徑雷達(dá)干涉測(cè)量技術(shù)原理，歸納了常規(guī)地基 InSAR形變測(cè)量的主要處理流程； 系統(tǒng)總結(jié)了地基干涉測(cè)量雷達(dá)系統(tǒng)的主要類(lèi)型、發(fā)展趨勢(shì)及地基合成孔徑雷達(dá)干涉測(cè)量技術(shù)的主要應(yīng)用領(lǐng)域，比較了地基和星載InSAR在形變測(cè)量上的優(yōu)勢(shì)和不足。

1)地基干涉測(cè)量雷達(dá)系統(tǒng)經(jīng)過(guò)十幾a的研發(fā)，核心傳感器從原來(lái)的SFCW雷達(dá)發(fā)展到目前的FMCW雷達(dá)、噪聲雷達(dá)、MIMO技術(shù)與FMCW雷達(dá)耦合技術(shù)，數(shù)據(jù)采集頻率由原來(lái)的幾分鐘縮短為幾秒鐘，甚至亞秒級(jí)，二維成像觀(guān)測(cè)由“近似動(dòng)態(tài)”發(fā)展到了動(dòng)態(tài)觀(guān)測(cè)目標(biāo)區(qū)的變形過(guò)程。設(shè)備可操作性和穩(wěn)定性越來(lái)越高，能在惡劣工作情況下持續(xù)觀(guān)測(cè)； 觀(guān)測(cè)方式更加靈活，可根據(jù)研究區(qū)形變特征自由設(shè)計(jì)最佳觀(guān)測(cè)方案，能從多個(gè)角度觀(guān)測(cè)和分析目標(biāo)區(qū)的運(yùn)動(dòng)特征； 數(shù)據(jù)處理簡(jiǎn)單快速，能在十幾分鐘內(nèi)處理完幾天的觀(guān)測(cè)數(shù)據(jù)； 監(jiān)測(cè)精度更高，觀(guān)測(cè)目標(biāo)的微小變化在地基雷達(dá)測(cè)量結(jié)果中均有所反映。

2)地基合成孔徑雷達(dá)干涉測(cè)量技術(shù)應(yīng)用涵蓋了露天采礦邊坡、崩塌、山體滑坡、火山、橋梁等結(jié)構(gòu)變形、壩體變形、冰川運(yùn)動(dòng)、區(qū)域數(shù)字高程模型(DEM)等十幾個(gè)領(lǐng)域，尤其是在露天邊坡穩(wěn)定性監(jiān)測(cè)領(lǐng)域技術(shù)成熟，能探測(cè)到邊坡的局部極其微小的變形，為工程施工提供早期預(yù)警。

3)相位解纏仍然是制約地基合成孔徑雷達(dá)干涉測(cè)量技術(shù)的主要因素。由于采用短波長(zhǎng)雷達(dá)波，增加了相位解纏的難度。新型雷達(dá)傳感器的應(yīng)用能夠有效提高數(shù)據(jù)采樣率，通過(guò)縮短數(shù)據(jù)采集時(shí)間間隔以降低形變梯度，從而可以提高相位解纏的準(zhǔn)確度。另外，大氣相位校正也是數(shù)據(jù)處理的一個(gè)難點(diǎn)，目前尚無(wú)可靠的大氣校正模型消除大氣相位影響。

[1] Monserrat Hernández O.Deformation Measurement and Monitoring with Ground-based SAR[D].Barcelona:Universitat Politècnica de Catalunya,2012.

[2] R?delsperger S.Real-time Processing of Ground Based Synthetic Aperture Radar(GB-SAR) Measurements[D].Hessen:Technische Universit?t Darmstadt,2011.

[3] Martínez-Vazquez A,Fortuny-Guasch J.A GB-SAR processor for snow avalanche identification[J].IEEE Transactions on Geoscience and Remote Sensing,2008,46(11):3948-3956.

[4] IDS Corporation.Interferometric radar(IBIS) technical papers[EB/OL].https://www.idscorporation.com/georadar/more-information/technical-papers.

[5] GAMMA Corporation.GAMMA Portable Radar Interferometer[EB/OL].http://www.gamma-rs.ch/rud/microwave-hardware/gpri.html.

[6] R?delsperger S,Coccia A,Vicente D,et al.Introduction to the new metasensing ground-based SAR:Technical description and data analysis[C]//Proceedings of 2012 IEEE International Geoscience and Remote Sensing Symposium(IGARSS).Munich,Germany:IEEE,2012:4790-4792.

[7] ARELECS Corporation.ARAMIS applications[EB/OL].http://www.arelecs.com/wordpress/?page_id=166.

[8] Groundprobe Corporation.Products and services[EB/OL].http://www.groundprobe.com/products-and-services.

[9] 趙東寅,申其鴻,馬海濤,等.國(guó)產(chǎn)地基合成孔徑雷達(dá)監(jiān)測(cè)預(yù)警系統(tǒng)在紫金山金銅礦露天采場(chǎng)邊坡位移監(jiān)測(cè)的應(yīng)用[J].中國(guó)安全生產(chǎn)科學(xué)技術(shù),2015,11(4):54-58. Zhao D Y,Shen Q H,Ma H T,et al.Application of domestic ground based synthetic aperture radar monitoring and warning system in slope displacement monitoring on open pit of Zijinshan gold/copper mine[J].Journal of Safety Science and Technology,2015,11(4):54-58.

[10]Martínez-Vázquez A.Snow Cover Monitoring Techniques with GB-SAR[D].Barcelona:Universitat Politècnica de Catalunya,2008.

[11]Broussolle J,Kyovtorov V,Basso M,et al.MELISSA,a new class of ground based InSAR system.An example of application in support to the Costa Concordia emergency[J].ISPRS Journal of Photogrammetry and Remote Sensing,2014,91:50-58.

[12]Iglesias González R.Orbital and Ground-based Differential SAR Interferometry for the Monitoring of Slow-Moving Landslides[D].Barcelona:Universitat Politècnica de Catalunya,2013.

[13]Lukin K,Mogyla A,Palamarchuk V,et al.Monitoring of St. Sophia Cathedral interior using Ka-band ground based noise waveform SAR[C]//Proceedings of the European 6th Radar Conference.Rome,Italy:IEEE,2009:215-217.

[14]Lee H,Lee J H,Kim K E,et al.Development of a truck-mounted arc-scanning synthetic aperture radar[J].IEEE Transactions on Geoscience and Remote Sensing,2014,52(5):2773-2779.

[15]曲世勃,王彥平,譚維賢,等.地基SAR形變監(jiān)測(cè)誤差分析與實(shí)驗(yàn)[J].電子與信息學(xué)報(bào),2011,33(1):1-7. Qu S B,Wang Y P,Tan W X,et al.Deformation detection error analysis and experiment using ground based SAR[J].Journal of Electronics & Information Technology,2011,33(1):1-7.

[16]張祥,陸必應(yīng),宋千.地基SAR差分干涉測(cè)量大氣擾動(dòng)誤差校正[J].雷達(dá)科學(xué)與技術(shù),2011,9(6):502-506,512. Zhang X,Lu B Y,Song Q.Atmospheric disturbance correction in ground-based SAR differential interferometry[J].Radar Science and Technology,2011,9(6):502-506,512.

[17]Lglesias R,Fabregas X,Aguasca A,et al.Atmospheric phase screen compensation in ground-based sar with a multiple-regression model over mountainous regions[J].IEEE Transactions on Geoscience and Remote Sensing,2014,52(5):2436-2449.

[18]Monserrat O,Crosetto M,Luzi G.A review of ground-based SAR interferometry for deformation measurement[J].ISPRS Journal of Photogrammetry and Remote Sensing,2014,93:40-48.

[19]楊紅磊,彭軍還,崔洪曜.GB-InSAR監(jiān)測(cè)大型露天礦邊坡形變[J].地球物理學(xué)進(jìn)展,2012,27(4):1804-1811. Yang H L,Peng J H,Cui H Y.Slope of large-scale open-pit mine monitoring deformations by using ground-based interferometry[J].Progress in Geophysics,2012,27(4):1804-1811.

[20]Del Ventisette C,Casagli N,Fortuny-Guasch J,et al.Ruinon landslide(Valfurva,Italy) activity in relation to rainfall by means of GB InSAR monitoring[J].Landslides,2012,9(4):497-509.

[21]Gischig V,Loew S,Kos A,et al.Identification of active release planes using ground-based differential InSAR at the Randa rock slope instability,Switzerland[J].Natural Hazards and Earth System Sciences,2009,9(6):2027-2038.

[22]Nolesini T,Di Traglia F,Del Ventisette C,et al.Deformations and slope instability on Stromboli volcano:Integration of GBInSAR data and analog modeling[J].Geomorphology,2013,180-181:242-254.

[23]Pieraccini M.Monitoring of civil infrastructures by interferometric radar:A review[J].The Scientific World Journal,2013,2013:786961.

[24]Mao W J,Chang W L.Deformation monitoring by ground-based SAR interferometry(GB-InSAR):A field test in dam[J].Advances in Information Sciences and Service Sciences,2015,7(2):133-140.

[25]Strozzi T,Werner C,Wiesmann A,et al.Topography mapping with a portable real-aperture radar interferometer[J].IEEE Geoscience and Remote Sensing Letters,2012,9(2):277-281.

(責(zé)任編輯： 李瑜)

Ground-based interferometric synthetic aperture radar and its applications

LIU Bin1, GE Daqing1, LI Man1, ZHANG Ling1, WANG Yan1, GUO Xiaofang1, ZHANG Xiaobo2

(1.ChinaAeroGeophysicalSurveyandRemoteSensingCenterforLandandResources,Beijing100083,China;2.CollegeofGeoscienceandSurveyingEngineering,ChinaUniversityofMiningandTechnology(Beijing),Beijing100083,China)

On the basis of a detailed discussion on the principle of GB InSAR, the main data processing and analysis stages for estimating deformations starting with the GB InSAR observations are described. This paper gives a review of the main types and development trend of ground-based radar system, the main application domain and some existent problems of GB InSAR, and then summarizes the pros and cons of ground-based and space-borne InSAR for deformation monitoring.

ground-based interferometric synthetic aperture radar; stepped frequency continuous wave(SFCW); frequency modulation continuous wave(FMCW); noise radar; multiple input multiple output(MIMO) technology

10.6046/gtzyyg.2017.01.01

劉斌，葛大慶，李曼,等.地基合成孔徑雷達(dá)干涉測(cè)量技術(shù)及其應(yīng)用[J].國(guó)土資源遙感,2017,29(1):1-6.(Liu B,Ge D Q,Li M,et al.Ground-based interferometric synthetic aperture radar and its applications[J].Remote Sensing for Land and Resources,2017,29(1):1-6.)

2015-07-22；

2015-10-28

國(guó)家自然科學(xué)基金項(xiàng)目“基于改進(jìn)的高分辨率時(shí)序InSAR技術(shù)研究Khash Mw7.7地震震后形變機(jī)制”(編號(hào)： 41504048)和中國(guó)地質(zhì)調(diào)查局計(jì)劃項(xiàng)目“地基雷達(dá)干涉測(cè)量系統(tǒng)與邊坡穩(wěn)定性監(jiān)測(cè)示范”(編號(hào)： 12120114036401)共同資助。

劉斌(1982-)，男，博士，主要從事星載和地基InSAR技術(shù)理論與應(yīng)用研究。Email： lbin0226@163.com。

TP 79

A

1001-070X(2017)01-0001-06