胡 波 汪漢勝 賈路路 王 平 胡 俊 樊文靜

(1)中國科學院測量與地球物理研究所動力大地測量學重點實驗室,武漢 430077 2)中國科學院研究生院,北京 100049 3)中南大學信息物理工程學院,長沙 410083 4)長安大學地質工程與測繪工程學院,西安 710054)

DInSAR技術監(jiān)測青藏高原凍土形變的試驗研究＊

胡 波1,2)汪漢勝1)賈路路1,2)王 平3)胡 俊3)樊文靜4)

(1)中國科學院測量與地球物理研究所動力大地測量學重點實驗室,武漢 430077 2)中國科學院研究生院,北京 100049 3)中南大學信息物理工程學院,長沙 410083 4)長安大學地質工程與測繪工程學院,西安 710054)

以青藏鐵路拉薩段為試驗區(qū),利用D InSAR技術和 Palsar數(shù)據(jù)監(jiān)測得到了該地區(qū)的凍土形變,對形變圖的分析表明,其形變非常吻合凍土凍脹融沉的物理變化規(guī)律,說明D InSAR技術適合在青藏高原大范圍內監(jiān)測凍土形變。

青藏高原;凍土;青藏鐵路;凍脹融沉;D InSAR

1 引言

合成孔徑雷達差分干涉 (D InSAR)是以合成孔徑雷達復數(shù)影像的相位信息獲取地表變化的技術,利用遙感衛(wèi)星多時相的單視復數(shù)雷達圖像的相干信息進行差分或利用外部DEM數(shù)據(jù)消除地形效應后,提取地面信息,最終達到探測地表微小形變的目的,視線向形變監(jiān)測精度可達到毫米級[1]。D InSAR已發(fā)展成為一種專門監(jiān)測地表形變的新技術,由于其具有高覆蓋、全天候和高精度等優(yōu)勢,在監(jiān)測大面積微小地面形變如地震、滑坡、泥石流、火山、構造活動、冰川消融、城市沉降等監(jiān)測和預報方面有著廣泛的應用。因此,D InSAR已成為監(jiān)測環(huán)境災害最有效的空間對地觀測手段之一。

青藏高原凍土面積占該地區(qū)總面積的 57%。凍土是一種對溫度極為敏感的土體介質,溫度的變化所引起的凍脹和融沉對凍土地區(qū)的工程安全有很大的影響。青藏高原大多屬高溫凍土,高原大氣溫度和地表溫度隨季節(jié)性變化較大,導致了青藏高原凍土的不穩(wěn)定[2]。

高原凍土的凍脹和融沉歸根結底是地表形變問題,我們嘗試采用D InSAR技術進行監(jiān)測,并對凍融區(qū)進行精確劃分,詳細給出凍土的凍脹融沉大小和深度。關于D InSAR監(jiān)測凍土,國內外學者已做了大量研究。我國于 2001年在青藏高原北麓河附近建立了凍土工程研究基地,李震等[3]基于重復軌道的 ERS1/2 SAR圖像,研究了利用干涉 SAR技術探測凍土形變的方法,得到了研究區(qū)凍土形變結果,通過與實測數(shù)據(jù)比較,表明重軌差分干涉測量可以精確地探測凍土表面形變。謝酬等[4]研究了利用PALSAR數(shù)據(jù)進行凍土形變檢測的方法,并對其檢測結果進行了分析,其結果也表明了該方法與水準測量方法有較好的一致性。本文將在前人工作的基礎上,選擇青藏鐵路經(jīng)過的拉薩及附近的某一段為試驗區(qū),研究青藏高原季節(jié)性的和長期變化的凍脹融沉規(guī)律。

2 研究區(qū)概況及數(shù)據(jù)來源

研究區(qū)內主要土地類型有冰川、積雪、裸巖、灌叢等,植被以高寒草原為主,混生有墊狀植物。該研究區(qū)主要是高溫多年凍土區(qū)。高溫多年凍土由于冷蓄量比較小,對外界環(huán)境變化比較敏感,穩(wěn)定性差,從而易引發(fā)凍脹融沉現(xiàn)象。凍土活動層的年變化過程一般分 4個階段：夏季融化、秋季凍結、冬季降溫、春季升溫[5]。凍土的凍結和融化隨溫度和季節(jié)的變化交替出現(xiàn),在這兩種因素的作用下,地表就會反復的隆升和下降,引起地表不均勻沉降,從而引發(fā)一系列的地質災害如滑坡、泥石流、建筑物塌陷等。

青藏高原地區(qū)地形、地貌、生態(tài)結構都比較復雜。我們知道波長越長,穿透性越強,故相對于 ERS和 ENV ISAT數(shù)據(jù)來說,日本ALOS衛(wèi)星采用了先進的對地觀測技術,可獲得更高分辨率的雷達影像。在高植被覆蓋和地形復雜多變地區(qū),L波段的 PALSAR數(shù)據(jù)可以保持很高的相干性,尤其對于形變梯度大的區(qū)域,不易引起干涉條紋的混疊,因此,PALSAR數(shù)據(jù)比較適合地形復雜區(qū)域和凍土形變監(jiān)測。故本試驗選擇 3景 PALSAR數(shù)據(jù),其成像時間分別為2008年 12月 17日、2009年2月 1日和 2009年 3月 19日。圖幅的 4個角坐標及中心坐標見表 1。

表1 PALSAR數(shù)據(jù)4角及中心坐標Tab.1 Coordinates of the four corners and the center of the PALSAR data

地形數(shù)據(jù)為美國宇航局提供的 90 m×90 m分辨率的DEM。圖 1為青藏高原試驗區(qū)的 SRT M數(shù)據(jù),圖 2為 SRT M數(shù)據(jù)的三維顯示。

圖1 青藏高原試驗區(qū)的 SRT M數(shù)據(jù)Fig.1 SRT M data of test area in the Qinghai-Tibetan plateau

圖2 SRT M數(shù)據(jù)的三維顯示Fig.2 Three-dimensional display of the SRT M

3 數(shù)據(jù)處理

鑒于衛(wèi)星的重返周期及復雜地形對地形相位的影響,外部DEM的輔助正好可以克服計算地形相位貢獻值的缺陷[6]。因此,選用雙軌法進行差分干涉提取地表形變。

以 2008年 12月 17日為主影像,2009年 2月 1日、2009年 3月 19日為從影像,分別進行干涉,形成 2個干涉對,分別命名為 20081217—20090201, 20081217—20090319,數(shù)據(jù)處理主要包括配準、干涉、濾波、解纏、地理編碼等 5個步驟[7-9]。相干圖、差分干涉圖、形變圖見圖 3～8。

圖 3 2008-12-17—2009-02-01相干圖Fig.3 Coherence map of 2008-12-17—2009-02-01

圖 4 2008-12-17—2009-03-19相干圖Fig.4 Coherence map of 2008-12-17—2009-03-19

圖 5 2008-12-17—2009-02-01差分干涉圖Fig.5 Differential interferogram of 2008-12-17—2009-02-01

圖 6 2008-12-17—2009-03-19差分干涉圖Fig.6 Differential interferogram of 2008-12-17—2009-03-19

圖 7 2008-12-17—2009-02-01形變圖Fig.7 Defor mation map of 2008-12-17—2009-02-01

圖 8 2008-12-17—2009-03-19形變圖Fig.8 Defor mation map of 2008-12-17—2009-02-01

從圖 7可以看出,研究區(qū)域發(fā)生了整體抬升,最大抬升量達 0.16 m。這是因為在 2008-12-17日—2009-02-01日時間段內,試驗區(qū)經(jīng)歷了冬季降溫過程從而引起凍土凍脹,故研究區(qū)域普遍發(fā)生抬升,同時,凍脹形變達到最大。從圖 8可以看出,研究區(qū)域發(fā)生了大范圍的沉降,最大沉降量達 0.28 m。這是因為從 2008-12-17日—2009-03-19日,天氣轉暖,凍土融沉而引起大范圍的沉降。從形變圖中,我們還可以分析出在沒有出現(xiàn)人為活動的區(qū)域,由于熱平衡保持完好,形變量相對較小,而在青藏鐵路所經(jīng)過的區(qū)域,地表發(fā)生了明顯的沉降,這是因為大型工程活動改變了多年凍土的熱量平衡,使得季節(jié)融化深度加深和次年回凍深度減薄。這一結果與馬小杰[10]給出的利用水準測量的結果比較穩(wěn)合。

4 結論與討論

研究驗證了 PALSAR數(shù)據(jù)在凍土形變監(jiān)測方面具有巨大的優(yōu)勢,為青藏鐵路的監(jiān)測提供了一種便捷手段。同時也可指導凍土地區(qū)的工程建設活動,可以通過 PALSAR數(shù)據(jù)的監(jiān)測結果來判斷凍土地區(qū)的形變大小是否超出可承受范圍,預防病害的發(fā)生。

下一步的研究將通過處理研究區(qū)域不同時段的數(shù)據(jù),得到研究區(qū)域的凍土形變規(guī)律,擬合出形變隨時間變化的關系。建立凍脹融沉預測模型,對各種特征區(qū)的時間序列形變觀測結果,進行最小二乘回歸分析,回歸模型包括季節(jié)性變化、半年變化和長期線性變化等成分,建立各種特征區(qū)特別是鐵路和公路的變形預測模型,可對凍土地區(qū)的凍融進行預測和預警。

致謝 感謝NASA與日本航空局提供 SRT M數(shù)據(jù)和PALSAR數(shù)據(jù)!

1 Xia Ye.Bam earthquake：Surface defor mation measurement using radar interferometry[J].Acta Seis mologica Sinica, 2005,18(4)：451-459.

2 答治華,童長江.建立青藏鐵路多年凍土區(qū)工程的監(jiān)測體系的必要性[J].鐵道工程學報,2002,(2)：68-70.(Da Zhihua and Tong Changjiang.Necessity of the establishment monitoring system in the permafrost regions of the Qinghai-Tibet Railway[J].Railway Engineering,2002,(2)：68-70)

3 李震,等.差分干涉 SAR凍土形變檢測方法研究[J].冰川凍土,2004,26(8)：389-396.(Li Zhen,et al.Differential interferometric SAR detection of permafrost deformation [J].glacial soil,2004,26(8)：389-396)

4 謝酬,等.基于 PALSAR數(shù)據(jù)的青藏高原凍土形變檢測方法研究[J].國土資源遙感,2008,9(3)：15-19.(Xie Chou,et al.A study of deformation in permafrost regions of Qinghai-Tibet Plateau based on AloS/PALSAR D InSAR Interferometry[J].Journal of Remote Sensing,2008,9(3)：15-19)

5 趙林,等,青藏高原五道梁附近多年凍土活動層凍結和融化過程[J].科學通報,2000,45(11)：1 205-1 211. (Zhao Lin,et al.Qinghai-Tibet Plateau permafrost Wudaoliang active layer near the freezing and thawing process [J].Science Bulletin,2000,45(11)：1 205-1 211)

6 LiaoM S,et al.Reconstruction ofDEMS From ERS1/2 Tandem Data in Mountainous Area Facilitated by SRT M Data [J].IEEE Transactions on Geoscience and Remote Sensing, 2007,45(7)：2 325-2 335.

7 胡波,汪漢勝.二軌法 D InSAR技術監(jiān)測城市地表沉降[J].測繪工程,2010,19(2)：37-41.(Hu Bo and Wang Hansheng.Two-pass method D InSAR technology to monitor ground subsidence[J].Engineering of Surveying and Mapping,2010,19(2)：37-41)

8 胡波,汪漢勝.D InSAR技術對地震同震形變場的研究[J].測繪工程,2010,19(1)：9-12.(Hu Bo and Wang Hansheng.Research on extract earthquake co-seis mic deformation by Differential Interferometry Synthetic Aperture Radar[J].Engineering of Surveying andMapping,2010,19 (1)：9-12)

9 胡波,汪漢勝,PSInSAR技術監(jiān)測地表形變的研究[J].大地測量與地球動力學,2010,(1)：34-39.(Hu Bo and Wan Hansheng.Monitoring ground subsidence with permanent scatterers interferometry[J].Geodesy and Geodynamics,2010,(1)：34-39)

10 馬小杰,等.青藏鐵路路基沉降變形的灰色預測模型研究[J].路基工程,2006(5)：1-3.(Ma Xiaojie,et al. The Qinghai-Tibet railway on the settlement of the gray prediction model[J].Foundation Engineering,2006,(5)：1 -3)

USING D InSAR TO MONITOR DEFORMATION OF FROZEN GROUND IN TIBETAN PLATEAU

Hu Bo1,2),Wang Hansheng1),Jia Lulu1,2),Wang Ping3),Hu Jun3)and FanWenjing4)

(1)Key Laboratory of Dynam ic Geodesy,Institute of Geodesyamp;Geophysics,CAS,W uhan 430077 2)Graduate Universty of Chinese Academ y of Sciences,Beijing 100049 3)Departm ent of Surveying and Land Infor m ation Engineering,Central South University,Changsha 410083 4)Shool of Geology Engineering and Geom atics of Chang’an University,Xi’an 710054)

The Lasa section of theQinghai-Tibet railwaywas selected as a test region.The defor mation of frozen ground in the test region were successfully monitored with D InSAR.Through the analysis of the deformation,it is found that the deformation conformes to the law of physical changes of frozen ground,frost heave and thawing settlement.It is shown that D InSAR technique has good prospects in defor mation monitoring for the frozen ground in Qinghai-Tibetan plateau.

Tibetan Plateau;frozen ground;Qinghai-Tibet rail way;frost heave and thawing settlement;D InSAR

1671-5942(2010)05-0053-04

2010-05-10

國家杰出青年科學基金(40825012)

胡波,男,1983年生,博士研究生,主要從事 InSAR研究.E-mail：hubo@asch.whigg.ac.cn

P227

A