何亞東,徐許雄,李 林,孫 銘,黃 逢

(1. 自然資源部第三大地測量隊,四川 成都 610100; 2. 四川省冶勘設(shè)計集團有限公司,四川 成都 610051)

地基合成孔徑雷達干涉測量系統(tǒng)(ground based interferometric synthetic aperture radar, GBInSAR)能對滑坡體進行快速監(jiān)測,具有靈活的觀測姿態(tài)和連續(xù)觀測能力,其分辨率高、平臺穩(wěn)定、觀測周期短、成本低廉,目前已應(yīng)用于滑坡[1-2]、露天煤礦[3]、大壩[4]和橋梁[5-6]等領(lǐng)域的形變監(jiān)測中。部分學(xué)者聯(lián)合北斗和星載InSAR進行了相關(guān)研究。文獻[7]提出了北斗+GNSS+InSAR集成長周期基礎(chǔ)設(shè)施毫米級形變監(jiān)測方法。文獻[8]聯(lián)合時序InSAR和GNSS技術(shù)對深圳和香港海岸線的地表形變場進行監(jiān)測, 識別出填海區(qū)的地鐵建筑區(qū)和地基較淺的建筑物,存在明顯的沉降現(xiàn)象。文獻[9]提出融合星載InSAR和GNSS的三維形變監(jiān)測方法。文獻[10]探討了基于北斗聯(lián)合光學(xué)、雷達等多源衛(wèi)星遙感手段的地質(zhì)災(zāi)害監(jiān)測技術(shù)體系。

綜合上述研究,本文聯(lián)合北斗衛(wèi)星定位技術(shù)和地基InSAR技術(shù),對四川省茂縣的渭門鎮(zhèn)特大滑坡進行監(jiān)測,并結(jié)合裂縫計、雨量計傳感器進行數(shù)據(jù)分析。

1 地基InSAR原理

1.1 基本原理

地基InSAR沿雷達視線方向(line of sight,LOS)和距離方向形成像元面,該面被分割成許多二維像元小單元,將每個像元作為觀測對象(如圖1所示),進而以面為單位實現(xiàn)高空間分辨率形變監(jiān)測。其實質(zhì)為利用一個小的雷達天線在滑軌上運動,在運動方向上獲得一個等效的大孔徑天線(稱為合成孔徑)[11]。通過這種方式,雷達對同一場景重復(fù)掃描成像,利用重復(fù)圖像對進行干涉,提取目標(biāo)區(qū)域的相位變化,再通過相位變化差值(如圖2所示)反演監(jiān)測場景的形變信息。其中,LOS方向的形變量Δd與干涉相位φdef的關(guān)系[12]可表示為

圖1 空間分辨單元

圖2 差分干涉SAR形變監(jiān)測技術(shù)

(1)

式中,λ表示雷達信號的波長。

根據(jù)差分干涉SAR的基本理論,假設(shè)有N幅地基 SAR 圖像,形成N-1幅干涉圖像,目標(biāo)點的干涉相位可以表示為

φInt=φdef+Δφatmo+Δφnoise

(2)

式中,φdef表示由形變引入的相位;Δφatmo表示大氣相位;Δφnoise表示由時間去相關(guān)等因素引入的隨機相位。

1.2 PS選點算法

針對傳統(tǒng)選點方法的誤選率高、效率低的特點,本文采用基于動態(tài)選擇永久散射點(persistent scatterers,PS)的實時處理算法[13],分組篩選出PS點。其大致流程為:假設(shè)分組數(shù)為15,首先選擇其主圖像前30幅,并以PS點作為背景場,并將第30幅圖像作為主影像,第31～45幅影像開始反演形變;第45幅圖像作為主影像,第46～60幅影像開始反演形變,并以此類推。

2 研究區(qū)域與數(shù)據(jù)獲取

2.1 滑坡描述

結(jié)合氣象資料,四川茂縣全縣境內(nèi)多年平均降雨量為490.7 mm,降雨量分布不均,降水集中在5—9月,占全年降水量的72%,是降水高峰期[14]。地貌以高山峽谷為主,地勢由西北向東南傾斜,水流方位呈東南流勢,山峰海拔均約為 4000 m。

該滑坡位于茂縣渭門鎮(zhèn),距茂縣13 km?；聻闅v史滑坡,所在坡體上部、下部都為聚居區(qū)。滑坡主要形變集中于前緣,前緣坡體結(jié)構(gòu)松散,存在懸空情況和明顯的強烈形變區(qū)。根據(jù)現(xiàn)場測算,強烈形變的區(qū)域坡體長200 m、寬350 m,滑坡體堆積厚度大于15 m,整個形變區(qū)域的方量約為90×104m3。

2.2 設(shè)備布設(shè)情況

針對滑坡的強形變區(qū),在范圍內(nèi)布設(shè)裂縫計設(shè)備LF01、LF02、LF03、LF04、LF05、LF06,雨量計設(shè)備YL01,基于北斗定位的GNSS設(shè)備GP01、GP02、GP03、GP04,以及GBInSAR設(shè)備對該區(qū)域進行監(jiān)測。

2.3 數(shù)據(jù)采集

基于北斗的GNSS設(shè)備為滑坡體安裝的固定設(shè)施,分別為GP01、GP02、GP03、GP04,該設(shè)備采用太陽能供電,24 h連續(xù)監(jiān)測。GBInSAR數(shù)據(jù)采集于2021年8月19日18:05—24日17:37,共獲得199景有效影像,采樣間隔為35 min。采樣數(shù)據(jù)距離向分辨率為0.375 m,方位向分辨率為4.0 mrad。

3 數(shù)據(jù)處理結(jié)果

3.1 傳感器數(shù)據(jù)結(jié)果分析

圖3為雨量計設(shè)備2021年8月1日—8月31日的實際降水量情況?？梢钥闯?該地區(qū)整個8月份有1/2的時間存在降雨過程。

圖3 雨量計降雨量

圖4表示LF01、LF02、LF03、LF04、LF05、LF06共6臺設(shè)備的拉力值,為2021年8月1日—8月31日的拉力值變化情況。可以看出LF02、LF03、LF04、LF06在不同時刻存在拉力值突變的情況,范圍為5～40 mm,表明滑坡體表面存在一定程度的形變。

圖4 裂縫計拉力值

3.2 北斗數(shù)據(jù)結(jié)果分析

圖5表示GP01、GP02、GP03、GP04于8月19日18:00—8月24日18:00采集的監(jiān)測數(shù)據(jù)。GP01、GP02、GP03設(shè)備累計形變量在20 mm以內(nèi),GP04設(shè)備累計形變達46 mm。

圖5 監(jiān)測點形變量

3.3 GBInSAR結(jié)果分析

GBInSAR結(jié)果如圖6所示,橢圓圈出的5處明顯形變區(qū)域為A1、A2、A3、A4、A5,取該區(qū)域內(nèi)形變最大的5個點位P001、P002、P003、P004、P005,其形變量監(jiān)測結(jié)果如圖7所示,其中P003位置累計形變量達272 mm。

圖6 案例形變累計量

圖7 地GBInSAR數(shù)據(jù)LOS向累計形變量

4 結(jié)果分析與討論

4.1 北斗與地基InSAR監(jiān)測數(shù)據(jù)結(jié)果分析

圖8為選擇相同觀測時間和觀測位置的監(jiān)測對比。其中,★表示北斗監(jiān)測GP04形變量變化,▲表示GBInSAR監(jiān)測GP04同一位置形變量變化,實線表示形變量趨勢線。

圖8 監(jiān)測對比

整個監(jiān)測周期內(nèi),北斗/GNSS設(shè)備監(jiān)測的GP04最大形變量為46.2 mm,平均值為22.3 mm,平均形變速度為0.35 mm/h;GBInSAR設(shè)備監(jiān)測的GP04最大形變量為49.7 mm,平均值為18.2 mm,平均形變速度為0.34 mm/h。對比分析發(fā)現(xiàn),兩種設(shè)備監(jiān)測結(jié)果基本一致。

4.2 監(jiān)測能力分析

北斗設(shè)備是對固定安裝位置進行監(jiān)測,可以安裝在滑坡體不同位置,達到對滑坡體的長期實時監(jiān)測,但缺點是只能做到點監(jiān)測。

GBInSAR設(shè)備可以實現(xiàn)對監(jiān)測面整體情況的監(jiān)測,既能監(jiān)測滑坡表面變化較大的區(qū)域塊,又能很好地監(jiān)控不同風(fēng)險等級,實現(xiàn)面監(jiān)測。其具有架設(shè)位置靈活、安裝簡單、快速出成果、實時監(jiān)測的優(yōu)點,能方便地應(yīng)用于應(yīng)急災(zāi)害監(jiān)測和次生災(zāi)害監(jiān)測。

5 結(jié) 語

北斗和GBInSAR聯(lián)合分析的方法用于滑坡監(jiān)測,能充分發(fā)揮兩種技術(shù)手段各自的優(yōu)勢,在滑坡形變的不同階段,從不同位置上分析滑坡表面形變信息,對于滑坡的早期防范和局部實時監(jiān)測有重大意義。本文通過四川茂縣渭門鎮(zhèn)滑坡這一案例,探討了北斗聯(lián)合GBInSAR進行形變監(jiān)測的模式。首先利用裂縫計、雨量計傳感器設(shè)備,證實滑坡存在下滑的現(xiàn)象;然后布設(shè)北斗區(qū)域網(wǎng)獲取實時監(jiān)測數(shù)據(jù),結(jié)合GBInSAR 24 h不間斷的監(jiān)測結(jié)果,識別出滑坡體存在多處不同程度的形變區(qū)域。北斗和GBInSAR聯(lián)合監(jiān)測,能多方位地識別出形變區(qū)域,將在今后滑坡穩(wěn)定性監(jiān)測和次生災(zāi)害監(jiān)測中發(fā)揮重要作用。