張英俊

【摘 要】傳統(tǒng)的地表變形監(jiān)測技術(shù)如GPS、水準(zhǔn)測量等，存在監(jiān)測范圍小、費用高和僅能獲取點目標(biāo)形變信息的局限性。D-InSAR技術(shù)具有全天時、全天候，能夠獲得面式數(shù)據(jù)等優(yōu)點。本文選取龍灘水電站附近36KM2的范圍為主研究區(qū)域，以Radarsat-2為數(shù)據(jù)源，采用GAMMA軟件對該地區(qū)5景數(shù)據(jù)進行了完整的差分干涉處理，提取變化量并對地表變形趨勢作出了分析，得出了幾處主要形變區(qū)域在時間上和空間上的形變趨勢，對庫區(qū)的安全監(jiān)測具有重要的意義。

【關(guān)鍵詞】D-InSAR；變形監(jiān)測；相位解纏

0 引言

變形在自然界中是普遍存在的，它是指變形體在各種載荷作用下，在時間上和空間上大小、形狀及其位置的變化。自然界的變形危害非常普遍，如地震、地表沉陷、滑坡、巖崩、潰壩橋梁與建筑物的倒塌等。它不僅嚴(yán)重破壞了自然環(huán)境和資源，而且對人民的生命財產(chǎn)安全造成嚴(yán)重威脅。近幾年，人類活動空間不斷擴大，大型交通、水利及資源開發(fā)工程項目大量實施，致使地表變形發(fā)生的頻率變得越來越高，造成經(jīng)濟損失和人員傷亡，必須要引起警惕。

龍灘水電站工程的施工建設(shè)期間以及蓄水運行以后，壩址基坑和庫區(qū)邊坡等都曾發(fā)生過巨大的變形，造成了嚴(yán)重危害，其穩(wěn)定性是大壩及水庫設(shè)計和建設(shè)的潛在問題之一，也是必須要考慮的重要方面。

1 地表變形常用監(jiān)測方法及其局限性

常用的地表變形監(jiān)測方法有很多種，主要包括高精度地面監(jiān)測技術(shù)、GPS監(jiān)測法及攝影測量方法等。下面對幾種重要的監(jiān)測方法作簡單的介紹：

1.1 GPS監(jiān)測法

GPS可以進行全天時全天候的監(jiān)測，不受氣候條件的影響，而且具有定位精度高、速度快、操作簡單的優(yōu)點。但它也有缺點：有比較嚴(yán)重的多路徑效應(yīng)，衛(wèi)星信號在復(fù)雜地區(qū)比較容易被遮擋，這些一定程度上影響測量的精度。

1.2 大地精密測量法

該方法是指采用高精度光電和光學(xué)測量儀器，如全站儀、精密水準(zhǔn)儀等儀器，通過測量角度、邊長和高程的變化來完成變形監(jiān)測任務(wù)。大地精密測量法是目前變形監(jiān)測的主要手段，具有如下優(yōu)點：能夠有效確定邊坡地表變形的范圍；能夠觀測到邊坡體的絕對位移量，靈活性大，精度高。缺點是地形和氣象條件對它的實施有很大的限制，而且周期長，進行連續(xù)觀測的能力較差，其工作量也比較大。

1.3 三維激光掃描技術(shù)

三維激光掃描技術(shù)它不需要接觸物體，白天和夜間都不影響外業(yè)工作，監(jiān)測速度快，測量精度高，能夠很好的反映坡體的總體變形趨勢。其缺點是，儀器價格非常昂貴，在山區(qū)儀器攜帶不太方便，測站需要有穩(wěn)定的架設(shè)場所，由于受到掃描距離的限制，必須要保證目標(biāo)和儀器之間沒有大的障礙。

1.4 光學(xué)遙感監(jiān)測法

利用遙感技術(shù)，不直接接觸地面目標(biāo)，根據(jù)其與電磁波的相互作用，收集其信息并且對其進行識別分類、判讀和分析，最終得到其變形信息。缺點是不能夠精確的提供地表變形動態(tài)信息，只是提供了宏觀的定性或者說是半定量的成果。

2 InSAR技術(shù)監(jiān)測地表變形的優(yōu)勢

合成孔徑雷達(dá)干涉測量技術(shù)（InSAR）是近幾年來國際遙感界的一個研究熱點，能以獲得某一地區(qū)厘米量級甚至更小尺度連續(xù)的地表形變信息，它們在地形制圖、土地利用分類、速度場探測、洪水監(jiān)測以及地表形變監(jiān)測等方面表現(xiàn)出極好的應(yīng)用前景，具體優(yōu)勢介紹如下：

2.1 全天時，全天候

因合成孔徑雷達(dá)（SAR）是主動式遙感，容易產(chǎn)生在時間序列上不缺失的雷達(dá)遙感圖像。

2.2 精度高

D-InSAR技術(shù)監(jiān)測地表微小形變的精度可以達(dá)到厘米級甚至更小尺度的量級。

2.3 時效性

隨著高精度衛(wèi)星星座的運營，衛(wèi)星之間的重訪間隔可以達(dá)到20秒，完全能達(dá)到實時監(jiān)測的要求，可實現(xiàn)對持續(xù)較慢地表變形活動的連續(xù)捕獲。

2.4 高空間分辨率

與傳統(tǒng)監(jiān)測方法得到的單點數(shù)據(jù)結(jié)果不同，SAR可以得到研究區(qū)域連續(xù)的面式地表形變信息。

2.5 大面積覆蓋，成本相對低廉。與其他傳統(tǒng)監(jiān)測方法相比整體成本較低。

3 三軌法D-INSAR技術(shù)的數(shù)據(jù)處理流程

D-InSAR技術(shù)原理是利用至少三次拍攝影像構(gòu)成兩個像對，對這兩個像對進行干涉處理而得到兩幅干涉條紋圖。對干涉圖進行二次差分處理，用來消除因地形起伏以及雷達(dá)成像的幾何而形成的相位差，則剩余的相位差即是在成像期間地表發(fā)生的變化所引起的。

三軌差分干涉需要三幅SAR圖像，將其中一幅圖像的坐標(biāo)系統(tǒng)作為參考坐標(biāo)系，那么在參考圖像和另外兩幅圖像之間將形成兩幅干涉圖。選擇相同圖像作為參考影像確保了生成的干涉產(chǎn)品都具有相同的坐標(biāo)系，而不需要進行額外的配準(zhǔn)和重采樣。三軌差分干涉只有一幅干涉圖包含了差分干涉的相位信息，另外一幅是作為參考消除由于地形所引起的相位變化，在精度和實效上要優(yōu)于二軌，其基本流程如圖1所示。

1）對選取的SLC數(shù)據(jù)進行精確的配準(zhǔn)。SLC數(shù)據(jù)即包括相位信息又包括強度信息，強度信息用于配準(zhǔn)，對相位信息進行處理就可以得到形變值。

2）主圖像與第二幅圖像經(jīng)過配準(zhǔn)后，經(jīng)干涉處理形成干涉像對1，然后復(fù)共軛相乘從而獲得干涉條紋圖。

3）干涉條紋圖中常常包含有噪聲干擾，而且平坦地表因沒有高程的變化，也會產(chǎn)生線性變化的干涉相位，稱之為平地效應(yīng)。平地效應(yīng)使干涉相位圖呈現(xiàn)為密集相見的條紋，一定程度上掩蓋了地形變化引起的干涉條紋變化。因此干涉條紋圖需要首先去除平地效應(yīng)，最后再把去平地效應(yīng)的干涉圖進行濾波處理。

4）對去平地效應(yīng)和濾波后的干涉條紋圖進行相位解纏。主圖像與第三幅圖像進行配準(zhǔn)，經(jīng)過干涉處理形成的干涉像對2，再去除平地效應(yīng)，進行濾波處理。

5）對解纏后得到的兩幅干涉條紋圖作差分處理，此時得到的結(jié)果是斜距坐標(biāo)，最后需經(jīng)過地理編碼轉(zhuǎn)換為地理坐標(biāo)系下的地表形變圖。

4 三軌法D-InSAR技術(shù)在龍灘水電站地表變形中的應(yīng)用

4.1 研究區(qū)域及數(shù)據(jù)介紹

研究區(qū)域處于龍灘水電站附近，采用的是Radarsat-2數(shù)據(jù)，綜合考慮雨季時間段，斜坡的傾斜角，分辨率，滑坡的朝向等因素，所選用數(shù)據(jù)是超精細(xì)HH極化產(chǎn)品模式、分辨率是3m、幅寬是6*6 km的降軌數(shù)據(jù) ，入射角為39.57°。

4.2 數(shù)據(jù)處理實例

采用三軌法對Radarsat-2數(shù)據(jù)進行差分干涉處理，影像1與影像2以影像1為參考形成干涉像對1，只包含地形信息，分別與影像3、4、5形成干涉像對2，包含了形變信息。

4.3 差分干涉圖的形成

①圖像配準(zhǔn)：利用上述影像，進行圖像配置流程如下：生成偏移值文件→偏移量的初始估計→生成精密偏差估計多項式→生成偏移多項式；最終配準(zhǔn)誤差為：距離向0.0847個像元 方位向0.0118個像元，配準(zhǔn)誤差小于0.1個像元，符合標(biāo)準(zhǔn)。

②生成干涉圖：利用GAMMA軟件生成干涉圖效果如圖3。

③去除平地效應(yīng)，然后經(jīng)過計算相干性→干涉圖濾波，獲得去平后的干涉圖3。

④相位解纏。利用GAMMA軟件，完成解纏掩膜的生成、自適應(yīng)降低掩膜的采樣、加權(quán)內(nèi)插填補解纏后圖像的缺口，用內(nèi)插后的解纏相位為模板解纏初始干涉圖3。

⑤差分干涉圖的生成：生成用于差分干涉的參數(shù)文件，然后相位轉(zhuǎn)換為形變圖。

4.4 變形區(qū)域的識別

在處理的過程中，很多環(huán)節(jié)都有可能產(chǎn)生誤差，導(dǎo)致獲得的SAR數(shù)據(jù)及生成的干涉圖產(chǎn)生偏差，從而影響分析結(jié)果。為了盡可能的去除誤差引起的錯誤值，把研究區(qū)域內(nèi)的處理結(jié)果根據(jù)相關(guān)系數(shù)的高低，位移值的大小等因素進行一定的過濾處理。本文選用了8月17日和10月25日的兩期數(shù)據(jù)來分析研究區(qū)域的形變情況，對比光學(xué)影像，找到變形的區(qū)域，并結(jié)合實地調(diào)查確定可信的變形區(qū)域，對變形作出分析。

區(qū)域1部分發(fā)生了沉降，并且形變比較明顯，呈現(xiàn)出一種時間上變化的連續(xù)性，而且變化的數(shù)值和范圍逐步變大。從空間上看，監(jiān)測結(jié)果中一直呈面狀分布，說明該處處于整體滑動狀態(tài)，后經(jīng)現(xiàn)場驗證此處為一處于活動狀態(tài)的古堆積體。區(qū)域2部分發(fā)生抬升，其在雨季的變動量較大，在此區(qū)域的上面可能發(fā)生了滑坡，松散土層、碎石土、風(fēng)化殼和半成巖土層的斜坡抗剪強度低，容易產(chǎn)生變形面下滑，這些土、巖下滑到此區(qū)域，導(dǎo)致出現(xiàn)地面抬升的現(xiàn)象。此移動區(qū)域離壩址較遠(yuǎn)，不會對庫區(qū)產(chǎn)生太大影響。

5 問題及其展望

D-InSAR處理過程中受到時間失相干、空間失相干和大氣影響等各種誤差的制約，雖本文在處理過程中結(jié)合了可信區(qū)域分析和過濾處理，這些處理能較好的減少一些無效的監(jiān)測值，但最后提取出變動區(qū)域的形變值還包括大氣效應(yīng)等引起的誤差值。這對區(qū)域地表位移的動態(tài)識別并不產(chǎn)生影響，卻會影響監(jiān)測位移值的精度。

【參考文獻(xiàn)】

[1]王桂杰，謝謨文，柴小慶，等.D-INSAR技術(shù)在庫區(qū)滑坡監(jiān)測上的實例分析[J].中國礦業(yè)，2011，20（3）：93-101.

[2]林子曄.GAMMA中ISP模塊SAR影像處理若干問題解析[J].北京測繪，2012，（3）：25-27.

[3]何秀鳳，何敏.InSAR對地觀測數(shù)據(jù)處理方法與綜合測量[M].北京：科學(xué)出版社，2011.

[4]Massonnet D，Briole P，Arnaud A.Deflation of Mount Etna monitored by spaceborne radar interferometry[J].Nature，1995，375：567-570.

[責(zé)任編輯：曹明明]