基于C波段和L波段的SBAS-InSAR技術監測廣州地面沉降

王會強,馮光財,喻永平,毛曉康,丁　超,陳晨月

(1.中南大學 雷達遙感研究室,湖南 長沙 410083；2.廣州市城市規劃勘測設計研究院,廣東 廣州 510000)

?

基于C波段和L波段的SBAS-InSAR技術監測廣州地面沉降

王會強1,馮光財1,喻永平2,毛曉康1,丁超1,陳晨月1

(1.中南大學 雷達遙感研究室,湖南 長沙 410083；2.廣州市城市規劃勘測設計研究院,廣東 廣州 510000)

近10年來廣州市頻繁發生地面沉降、塌陷等地質災害，造成巨大的生命財產損失，而且目前仍然在加劇。傳統的監測技術如GPS、水準測量等難以開展大范圍、高精度和高空間分辨率的地表沉降監測工作，而合成孔徑雷達干涉測量技術(InSAR)正逐漸成為城市地表沉降監測的有效手段。文中采用短基線集(SBAS-InSAR)技術，通過17景ENVISAT/ASAR數據和21景ALOS/PALSAR數據，探測廣州佛山地區2006—2011年的地表形變信息。將其與研究區內已有的水準測量數據進行比較，從而驗證InSAR技術監測結果的可靠性。最后圈定了研究區的重點沉降區域并對沉降成因進行分析。

廣州地區；地面沉降；SBAS-InSAR；形變場

珠江三角洲特別是廣州地區是我國經濟最發達的地區之一，但同時也是遭受地面沉降、塌陷等地質災害比較嚴重的區域。近10年來廣州及其周邊地區頻繁發生地面沉降、塌陷，給當地造成了巨大的生命財產損失，而且這種趨勢還在惡化，引起了相關部門的重視。傳統的監測技術如GPS、水準測量等很難實現大范圍、高精度和高空間分辨率的監測，而近些年發展起來的InSAR技術在城市地表沉降監測中得到廣泛應用，它不僅彌補了傳統測量手段的不足，而且兼具高精度、省時省力等監測優勢[1-2]。雖然目前已經有部分工作利用InSAR技術監測廣州及其周邊地區的沉降，獲得了覆蓋局部區域某時間段的平均形變速率圖[3-6]，但他們所采用的SAR數據量較少，平臺比較單一，也缺乏系統的比較驗證。這里同時采用C波段和L波段的SAR數據來監測廣州佛山地區的沉降，并利用水準數據對本文獲得的結果進行驗證。L波段較C波段穿透性更強，即使在植被覆蓋區也能具有很好的相干性[7-8]。

為了準確把握這一地區地面沉降動態和沉降區分布狀況及規律，建立完善的地面沉降地質災害監控體系，為政府制定具體防治措施和保證重要大型工程項目的建設安全提供參考資料，本文基于時序InSAR監測技術中的短基線集(SBAS-InSAR)技術[9]，通過17景ENVISAT/ASAR數據(時間跨度為2006—2009年)和21景ALOS/PALSAR數據(時間跨度為2007—2011年)，分別獲得了廣州、佛山地區在影像時間跨度內的形變場，兩組數據的結果可以相互對比驗證。然后，以研究區內已有的水準數據為基準，與InSAR監測結果進行比較，從而驗證了監測結果的可靠性。最后，圈定了研究區的重要沉降區域，并對監測結果的特征和成因進行分析。

1　研究區概況和數據

本文的研究區域涵蓋了廣州、佛山兩個市的城區，即從東經112°46′至113°36′，北緯22°48′至23°36′的范圍，該地區地勢大致呈現北部高、南部低，地面起伏不大。研究區屬于南亞熱帶典型的季風氣候，冬夏季風交替，雨量充沛，而且每年的3月至10月為雨季，降雨量大，而每年的11月至次年的2月，降雨量小，區內各地多年平均降雨量為1 665.0～1 678.0 mm。研究區內廣州番禺區和佛山北滘鎮、陳村一帶軟土層分布廣泛且厚度不均，廣花盆地地區地下隱伏巖溶(灰巖)分布多、發育強烈，特別是淺部巖溶[10-14]。研究區SAR影像覆蓋范圍如圖1所示。

圖1　研究區SAR影像覆蓋范圍ASAR數據的覆蓋范圍和PALSAR數據的覆蓋范圍，研究區則是兩種數據覆蓋相同的區域。

實驗采用ENVISAT/ASAR和ALOS/PALSAR影像，時間跨度分別為2006-10-30至2009-08-10和2007-01-08至2011-01-19。具體數據參數如表1所示。

表1　ASAR和PALSAR數據參數列表

續表1

2　InSAR數據處理

利用GAMMA軟件將SAR數據配準重采樣到相同的空間坐標系下，然后基于較短的時空基線閾值選取合適的干涉對進行DInSAR的數據處理，其中所選的干涉對和基線集情況如圖2所示。數據處理中，將ASAR數據進行多視處理為2×10(影像分辨率約為40 m×40 m)，將PALSAR數據進行多視處理為3×8(影像分辨率約為25 m×25 m)。采用的外部DEM是美國宇航局(NASA)提供的SRTM數據，分辨率是30 m×30 m，利用該DEM來去除干涉圖中的地形相位部分。而軌道誤差的去除則采用多項式曲面擬合的方法，作為常用的軌道誤差改正模型，該方法對于該地區的軌道趨勢去除效果比較明顯，特別是PALSAR數據[15-16]。該地區大氣影響嚴重，而且缺少相同時間段的用于改正大氣的外部數據源，所以根據干涉圖中大氣誤差的時空特性(空間域低通，時間域高通)，采用時空濾波的方法進行大氣延遲相位的去除。短基線集(Short BAseline Subsets，SBAS)作為一種先進的InSAR時間序列分析技術，解決了DInSAR的應用局限性(時空去相干，大氣延遲等)，為InSAR技術的拓展研究和應用提供了廣闊的空間。該技術將多景解纏后的多視差分干涉圖進行相位建模，通過奇異值分解(SVD)的方法最終生成平均形變速率圖和時間序列形變圖[2,9,16]。

圖2　基線集

3　監測結果與驗證

3.1監測結果

如圖3所示，根據相同覆蓋區域的ASAR數據和PALSAR數據，分別獲取了廣州地區2006年至2011年的地面沉降速率場。

圖3顯示，ALOS/PALSAR影像相干性較好，比ENVISAT/ASAR影像獲取的監測點更多，這與前者的穿透能力較強有關。另外，發現該地區存在局部小范圍的沉降，沉降速率大約在4～14 mm/a。圖3中結果是利用不同的SAR數據采用相同的InSAR形變監測方法獲得的，通過對比幾乎可以認為兩個數據探測到的形變場是一致的，這從SAR數據內部驗證了SBAS-InSAR技術用于監測廣州地區地表沉降結果的可靠性。

圖3　SBAS-InSAR技術監測結果

3.2形變結果驗證分析

為了進一步驗證監測結果的可靠性，將其與研究區已有水準數據進行比較。由于InSAR 監測的形變結果為雷達視線向的形變，而水準數據則是垂直向的形變，二者存在著差異，故需要利用式(1)把水準數據垂直方向的形變量轉換為雷達視線向。

(1)

式中：dLOS為轉換后的視線向的水準數據，d垂直向為轉換前垂直向的水準數據，θ為雷達入射角。

然后以水準結果為基準，即視為真值，將轉換后的InSAR結果與其進行比較。這里給出已有水準點的分布情況，如圖4所示。

(a)空間分布　　　　　　　　　　　　　　　　　　　(b)比較驗證圖4　水準點的空間分布和比較驗證

圖4(a)中顯示藍色方框為上述兩種SAR數據覆蓋的空間范圍，圖中紅色和白色符號代表用于比較的62個水準點，其中紅色三角符號共有14個，代表發生形變的水準點，而白色圓形符號則代表48個沒有發生形變的水準點。由于已有的水準數據只有一定時間段內的沉降量，所以將其轉換為年平均形變速率，利用InSAR獲取的平均形變速率分別與沒有發生形變的水準點和發生形變的水準點進行比較，結果如圖4(b)所示。圖4(b)顯示，對于穩定水準點ASAR數據監測的平均形變速率與水準數據之差的均值達到了-0.9 mm/a，其標準差(精度)為0.8 mm/a，而PALSAR數據的結果與水準數據之差的均值為-0.5 mm/a，與ASAR結果相差不大，其精度也達到了0.7 mm/a。對于形變水準點，我們發現SAR數據的精度較穩定點的監測精度低，ASAR數據較差的均值和標準差分別為0.3 mm/a和1.6 mm/a；而PALSAR數據較差的均值和標準差則為-0.6 mm/a和1.2 mm/a。總體而言，獲取的InSAR結果可靠性較高。

就ASAR和PALSAR數據而言，理論上對于同一區域，二者獨立監測的結果應該是一致的，結果中在形變大小和形變范圍方面也幾乎相同，與水準數據比較也說明ASAR數據和PALSAR數據都能達到很高的精度。然而兩種數據所采用的波段不同，ASAR數據采用C波段(波長約為5.6 cm)和PALSAR數據(波長約為23.6 cm)采用L波段，由于L波段穿透性更強，使得L波段數據較C波段即使在相干性較低的區域(如廣花盆地一帶的植被覆蓋區)也能獲得更多的監測信息，相干性更高。

3.3沉降成因分析

考慮到兩種數據共同監測的形變區域，在圖3中共圈定了3個主要形變區域，即佛山市北滘鎮陳村一帶(a區)、廣花盆地一帶(b區)和佛山丹灶鎮一帶(c區)。初步分析a、c區域沉降多與軟土層的分布和城市化的快速發展(如地鐵線和高樓大廈的興建)有關[10-12]。經查，該沉降區域多為工廠區和密集居民樓區。有可能是軟土層厚度分布不均及地基下土層承載力較低，隨著建筑區的荷載不斷增加和軟土地基受力不均等原因引起地面沉降。如2008-07-21，佛山市順德區北滘職業技術學校全校幾乎所有的建筑都出現了不同程度的沉降現象，并導致許多墻體開裂。

另外在廣州廣花盆地一帶(b區)也存在局部沉降區，如金沙洲、大坦沙地區以及龍歸鎮等，該區沉降范圍較小。廣花盆地是地面沉降塌陷的重點監測區，如金沙洲、大坦沙地區地下隱伏巖溶(灰巖)分布多、發育強烈，而且江村水源地和新華水源地是該地區地下水集中開采地帶，近幾十年來，由于這兩個水源地大量抽取地下水，已造成水源地附近的多處農田區的生產井附近發生塌陷沉降[12]。加之，近年來廣州城市總體規劃和發展戰略確定該地區為廣州發展的重點，地鐵、高速等重大工程項目不斷在此開展[12-14]。據調查，2006-08以來，受武廣高鐵金沙洲隧道施工抽排地下水的影響，金沙洲地區發生了數十起巖溶地面塌陷、地面沉降地質災害，導致多個單位，居民區的數十棟樓房相繼發生傾斜、變形。

正是因為這些復雜多變的地質環境和人類活動的綜合作用，致使該地區出現地面沉降、塌陷、裂縫等地質災害。

4　結束語

本文首先基于相同覆蓋區域的ASAR數據和PALSAR數據，分別獲取了廣州地區2006年至2011年的地面沉降速率場，兩種數據的監測結果一致性較高。然后利用研究區內已有的水準數據對InSAR技術監測的形變結果進行了比較分析，驗證了InSAR技術獲取結果的可靠性。最后，確定了廣州地區的主要沉降區域，即佛山北滘鎮、丹灶鎮一帶是研究區中沉降較為明顯的地區，沉降速率最大可達到14 mm/a左右，而廣花盆地一帶存在小范圍局部沉降。據相關資料分析，前者沉降可能與該地區軟土層分布有關，而后者沉降估計與該地區的隱伏巖溶分布有關。但是這些復雜的地質環境應該是地面沉降、塌陷發生的前提條件，沉降分析暗示人類活動與這些災害的發生有很大關系。

[1]MASSONNET D, ROSSI M, CARMONA C, et al. The displacement field of the Landers earthquake mapped by radar interferometry. Nature, 1993, 364(6433):138-142

[2]許文斌, 李志偉, 丁曉利，等. 利用InSAR短基線技術估計洛杉磯地區的地表時序形變和含水層參數[J]. 地球物理學報, 2012, 55(2):452-461.

[3]ZHAO Q, HUI L, JIANG L M, et al. A study of ground deformation in the guangzhou urban area with persistent scatterer interferometry. [J]. Sensors, 2009, 9(1):503-518.

[4]WANG, HUA, WRIGHT, et al. InSAR reveals coastal subsidence in the Pearl River Delta, China[J]. Geophysical Journal International, 2012, 191(3):1119-1128.

[5]CHEN F, LIN H, ZHANG Y, et al. Ground subsidence geo-hazards induced by rapid urbanization: implications from InSAR observation and geological analysis[J]. Natural Hazards & Earth System Sciences, 2012, 12(4):935-942.

[6]AO M, WANG C, XIE R, et al. Monitoring the land subsidence with persistent scatterer interferometry in Nansha District, Guangdong, China[J]. Natural Hazards, 2015，75(3)：2947-2964

[7]STROZZI T, WEGMULLER U, WERNER C L, et al. JERS SAR interferometry for land subsidence monitoring[J]. IEEE Transactions on Geoscience & Remote Sensing, 2002, 475(7):1702-1708.

[8]TAKEUCHI S, YAMADA S. Comparison of InSAR capability for land subsidence detection between C-band and L-band SAR[C]// Geoscience and Remote Sensing Symposium, 2002. IGARSS '02. 2002 IEEE International. IEEE, 2002:2379-2381 vol.4.

[9]BERARDINO P, FORNARO G, LANARI R, et al. A new algorithm for surface deformation monitoring based on small baseline differential SAR interferograms[J]. Geoscience & Remote Sensing IEEE Transactions on, 2002, 40(11):2375-2383.

[10] XIA F. Geological hazard and geological environment in guangdong province. Journal of Natural Disasters,1995.

[11] LIU H P, WANG Y L, LIU J L, et al. Cause mechanism and spatiotemporal distribution of major geological disasters in Guangzhou[J]. Journal of Natural Disasters, 2005,14(5):149-153.

[12] 易順民, 梁池生. 廣東省地質災害防治[M]. 北京:科學出版社, 2010.

[13] 鄭小戰. 廣花盆地巖溶地面塌陷災害形成機理及風險評估研究[D].長沙：中南大學, 2010.

[14] 黃健民, 鄧雄文, 胡讓全. 廣州金沙洲巖溶區地下水位變化與地面塌陷及地面沉降關系探討[J]. 中國地質, 2015,42(1)：300-308.

[16] 李珊珊,李志偉,胡俊,等. SBAS-InSAR技術監測青藏高原季節性凍土形變[J].地球物理學報, 2013, 56(5):1476-1486.

[責任編輯：劉文霞]

Ground subsidence monitoring in Guangzhou with SBAS-InSAR technique based on C- and L-band

WANG Huiqiang1, FENG Guangcai1, YU Yongping2, MAO Xiaokang1, DING Chao1, CHEN Chenyue1

(1.Laboratory of Radar Remote Sensing, School of Geoscience and Info-Physics, Central South University, Changsha 410083,Chaina; 2.Guangzhou Urban Planning & Design Survey Research Institure, Guangzhou 510000,China)

In recent 10 years, frequently reported geological disasters such as ground subsidence have occurred in Guangzhou area, causing huge damages to people’s life and property. More seriously, these disasters are continuing and further developing. However, traditional monitoring means, such as leveling and GPS surveying cannot meet the requirements of large-scale, high-resolution and high-precision monitoring for the surface subsidence. InSAR technique has gradually become an efficient way in monitoring surface subsidence. This paper uses the SBAS-InSAR technology based on high coherent and strong steady points to monitor the deformation field of Guangzhou and Foshan areas from 2006 to 2011 with 17 ENVISAT/ASAR images and 21 ALOS/PALSAR images respectively. Then the InSAR results are compared with the leveling data in the research area to validate InSAR technique’s reliability. Ultimately, analysis is made on the main subsidence areas and the reasons of ground subsidence.

Guangzhou area; ground subsidence; SBAS-InSAR; deformation field

10.19349/j.cnki.issn1006-7949.2016.11.012

2016-01-18

廣州市珠江科技新星項目(2013J2200072)；國家自然科學基金資助項目(41574005)

王會強(1990-)，男，碩士研究生.

P642.26

A

1006-7949(2016)11-0060-05