黃土灌溉區域地面沉降InSAR監測結果分析

何 楊，趙超英,2*，張 勤,2（.長安大學地質工程與測繪學院，陜西?西安 70054；2.西部礦產資源與地質工程教育部重點實驗室，陜西?西安 70054）

?

黃土灌溉區域地面沉降InSAR監測結果分析

何 楊1，趙超英1,2*，張 勤1,2
（1.長安大學地質工程與測繪學院，陜西?西安 710054；2.西部礦產資源與地質工程教育部重點實驗室，陜西?西安 710054）

摘 要：合成孔徑雷達是近年發展起來的一種空間對地觀測技術，目前已被廣泛應用于獲取地表信息及地表形變。地面沉降形變緩慢且影響范圍大，已成為國內外城市主要地質災害。研究區域位于涇河南岸下游，距陜西省西安市約50km。通過選取覆蓋涇陽南塬黃土灌溉區的11景X波段TerraSAR數據，采用短基線集技術（SBAS-InSAR技術）進行地面沉降監測。結果表明涇陽南塬灌溉區具有兩個地面沉降區域，其在121天內沿剖線累計形變分別達到了6.4cm和2.7cm。經對地面沉降區域相關的水井分布情況的實地調查，分析了形變中心形變速率和附近水井地下水水位高程資料，發現地面沉降由灌溉引起合成孔徑雷達是近年發展起來的一種空間對地觀測技術，目前已被廣泛應用于獲取地表信息及地表形變。地面沉降形變緩慢且影響范圍大，已成為國內外城市主要地質災害。研究區域位于涇河南岸下游，距陜西省西安市約50km。通過選取覆蓋涇陽南塬黃土灌溉區的11景X波段TerraSAR數據，采用短基線集技術（SBAS-InSAR技術）進行地面沉降監測。結果表明涇陽南塬灌溉區具有兩個地面沉降區域，其在121天內沿剖線累計形變分別達到了6.4cm和2.7cm。經對地面沉降區域相關的水井分布情況的實地調查，分析了形變中心形變速率和附近水井地下水水位高程資料，發現地面沉降由灌溉引起。

關鍵詞：地面沉降；合成孔徑雷達干涉測量；小基線集技術；時間序列；黃土灌溉區

地面沉降是在自然因素和人為因素的共同作用下，由于地殼表層土體壓縮而導致的區域性地面標高緩慢下沉的一種環境地質現象[1]，是地質環境系統破壞所導致的惡果，嚴重時會發展成為一種地質災害[2]。據統計我國目前已先后有21個省90多個城市或地區發生不同程度的地面沉降，地面沉降總面積達到了近10萬km2，形成了華北平原、長江三角洲地區、汾渭盆地等地面沉降現象嚴重的典型區域。該類地面沉降主要是由于過量抽取地下水而引起，具有影響范圍大，緩變型且難以規避等特點，引起政府管理部門及科研學術界的高度重視，并開展了大量調查監測和機理研究工作[3～5]。

而另有一類規模相對較小的、由于農田灌溉引起地面沉降，其形變機理與城市地面沉降相反，還沒有引起足夠重視。如位于陜西涇陽南塬的黃土塬農田區域，由于季節性灌溉，地下水位呈現上升趨勢，但該區域仍存在明顯的持續的地面沉降現象，而且在1km范圍的黃土塬邊形成系列黃土滑坡區。對于農田灌溉引起的地面沉降的監測，對于研究地面沉降機理、保護農田以及研究其影響范圍內的黃土滑坡均具有重要理論意義和實際價值。

合成孔徑雷達干涉測量（InSAR）技術在相干性很好的城市區域得到成功的應用[6～10]，而在農田等相干性較低的區域進行形變監測需要采取一些特殊措施。本文研究了基于最短時間間隔的德國空間局TerraSAR數據，采用小基線集（SBAS）InSAR技術監測獲取了地面沉降的累計形變信息和時間序列結果，通過與區域歷史地下水變化結果比較，發現了灌溉型地面沉降的特殊形變機理。

1　SBAS InSAR技術原理

SBAS InSAR技術通過設置時間和垂直基線閾值，來減弱時間和基線引起的失相干[11]。增加SAR影像數據，可以有效提高InSAR形變監測精度。假設第j幅干涉圖是由tA和tB時刻獲得的兩幅SAR圖像產生的，并已去除了地形相位部分，假設tB>tA，在方位距離像素坐標系(x,r)中，(x,r)處的差分解纏相位可以表示為：

其中λ為波長，M為干涉對個數。假設有N景影像，于是可得M滿足下列關系式（假設N為奇數）：

如果考慮到不同時間間隔的非線性形變，tA和tB時刻間的累積形變，可以表示成如下形式：

k為tA和tB時刻間SAR數據序列。最終得到了M個解纏干涉對，不同SAR數據間的累積形變將通過最小二乘（LS）或奇異值分解（SVD）獲得。

2　SBAS數據處理

2.1 研究區域概況

陜西涇陽南塬位于涇河下游南岸，橫跨太平、蔣劉及高莊3個行政鄉，東西長達27.1km，塬區面積約70 km2[12]。由于臺塬抬升和涇河不斷側蝕，右岸階地逐漸缺失，塬邊形成了高陡的斜坡帶，陡坡高為40～90m，坡度為40～70°[13]。

2.2 實驗數據及處理

本文共獲取2011年11月至2012年3月近121天的TerraSAR數據。該時間段數據覆蓋兩個地面沉降區域，總數據量11景。

首先，將所有從影像統一與主影像進行粗、精配準，所選擇的主影像拍攝時間為2012-01-01；通過設置時間、空間基線閾值來削弱時空失相關對干涉對的影響，本文所選的時間基線閾值為66天，空間基線閾值為±300m，共得到14對符合條件的干涉對，如圖1所示。以其中一個干涉對為例進行處理過程說明，先將精確配準后主從影像進行干涉處理，為抑制噪聲，本文進行了方位向距離向為2:2的多視處理，得到的干涉對的空間分辨率約為3m×3m；由等高線插值成10m分辨率的DEM，去除地形相位，從而獲得差分干涉圖；利用改進的自適應濾波方法對所有差分干涉圖進行頻域濾波；通過掩模方式削弱或避免低相干區域（海域或者孤島）對相位解纏的影響，使其不參與相位解纏，同時通過設置一定的相干性閾值（>0.2）來選擇相對穩定的點進行最小費用流（MCF）解纏。其中，某些解纏后的相位圖中仍殘留部分的非線性軌道相位，因此進一步采用二次多項式對其進行去除。最后通過SBAS-InSAR技術得到研

圖1　基線分布圖Fig.1 Interferometric SAR configurations for time series deformation calculation

3　結果及分析

采用SBAS技術對上述數據進行解算，得到研究區域的地面沉降速率圖。圖2為研究區域121天地面沉降累積形變圖。編號Ⅰ位于蔣劉村正東邊，編號Ⅱ位于余家堡與蔣劉村之間。

圖2　121天累積形變Fig.2 Cumulative deformation maps in 121 days

兩個地面沉降區域Ⅰ和Ⅱ分別距離塬邊約620m和330m，其剖面線AB和CD如圖3所示。由圖3(a)可得，地面沉降Ⅰ呈條帶狀展布，長和寬分別為570m和130m，沿剖面線AB最大累積形變為6.39cm。由圖3(b)可得，地面沉降Ⅱ長和寬分別為360m和230m，沿剖面線CD最大累積形變為2.65cm。

圖3　地面沉降區域剖面線形變時間序列Fig.3 The time series of deformation along profile in land subsidence area

資料顯示涇陽南塬分布著許多水井，該監測區域分布有四口水井（如圖2）。地面沉降Ⅱ位于井3與井4之間，地面沉降Ⅰ與井2相近，而井1位于地面沉降Ⅰ的東南方向。井1、井3、井4為灌溉水井，而井2為飲用水井。

選取地面沉降Ⅰ的形變中心，計算該點的地面沉降速率（如圖4），該速率呈現加速趨勢。2012年3月與2012年5月野外調查顯示井2處地下水位上升了5.7m[14]。由此可見，該地區的沉降非飲用取水所致，究其原因是由于農田灌溉引起的。大量的農田灌溉使得黃土的含水量不斷上升，導致黃土本身的結構被破壞，黃土顆粒之間的空隙減小，在地表表現為地面沉降。

圖4　地面沉降I形變速率時間序列與井2地下水位高程Fig.4 The land subsidence rate and groundwater level elevation

4　結論

本文通過2011年11月至2012年3月11景TerraSAR數據，采用SBAS-InSAR技術成功獲取研究區域地面沉降地理位置、空間分布、累積形變和形變時間序列結果。野外實地調查發現涇陽南塬地面沉降與灌溉引起的地下水上升有密切關系。

致謝

感謝德國空間局DLR提供的TerraSAR數據（No.LAN2471）。本文由國家重點基礎研究計劃項目（973計劃）（No.2014CB744703）和中國自然科學基金（No.41372375）聯合資助。

參考文獻(References)

[1] Galloway D L.國際地面沉降研究現狀與發展趨勢評述[J].上海國土資源,2014,35(2):1-8.Galloway D L.The current situations and developing trend of international land subsidence research[J].Shanghai Land &Resources,2014,35(2):1-8.

[2] 鄭銑鑫,武強,侯艷聲,等.關于城市地面沉降研究的幾個前沿問題[J].地球學報,2002,23(3):279-282.Zheng X X,Wu Q,Hou Y S,et al.Some frontier problems on land subsidence research[J].Journal of earth,2002,23(3):279-282.

[3] 方正,武健強,趙建康.長江三角洲地區地面沉降監測[J].上海地質,2003,24(2):1-4.Fang Z,Wu J Q,Zhao J K.Land subsidence monitoring of Yangtze delta region[J].Shanghai Geology,2003,24(2):1-4.

[4] 趙超英,張勤,丁曉利,等.基于InSAR的西安地面沉降與地裂縫發育特征研究[J].工程地質學報,2009,17(3):390-392.Zhao C Y,Zhang Q,Ding X L,et al.InSAR based evaluation of land subsidence and ground fissure evolution at Xi’an[J].Journal of Engineering Geology,2009,17(3):390-392.

[5] 葛大慶,郭小方,張玲,等.華北平原地面沉降區InSAR監測[R].北京:中國國土資源航空物探遙感中心,2010.Ge D Q,Guo X F,Zhang L,et al.Subsidence monitoring with InSAR in the North China Plain[R].Beijing:China Aero Geophysical Surveying and Remote Sensing Center for Land and Resources,2010.

[6] 李德仁.InSAR技術進步與地面沉降監測應用[J].上海國土資源,2013,34(4):1-6.Li D R.InSAR:Technological progress and its application to land subsidence monitoring[J].Shanghai Land &Resources,2013,34(4):1-6.

[7] 廖明生,裴媛媛,王寒梅,等.永久散射體雷達干涉技術監測上海地面沉降[J].上海國土資源,2012,33(3):5-10.Liao M S,Pei Y Y,Wang H M,et al.Subsidence monitoring in Shanghai using the PS-InSAR technique[J].Shanghai Land &Resources,2012,33(3):5-10.

[8] 張靜,張勤,曲菲霏.運城市地面沉降SBAS-InSAR監測和敏感性GIS分析[J].上海國土資源,2012,33(1):58-61.Zhang J,Zhang Q,Qu F F.Monitoring land subsistence using SBAS-InSAR and GIS to conduct sensitivity analysis in Yuncheng[J].Shanghai Land &Resources,2012,33(1):58-61.

[9] 孫赫,張勤,楊成生,等.PS-InSAR技術監測分析遼寧盤錦地區地面沉降[J].上海國土資源,2014,35(4):68-71.Sun H,Zhang Q,Yang C S,et al.PS-InSAR monitoring of land subsidence in the Panjin region of Liaoning province[J].Shanghai Land &Resources,2014,35(4):68-71.

[10] 劉媛媛,張勤,趙超英,等.基于多源SAR數據的InSAR地表形變監測[J].上海國土資源,2014,35(4):71.Liu Y Y,Zhang Q,Zhao C Y,et al.The surface deformation monitoring based on multi-source SAR data of InSAR[J].Shanghai Land &Resources,2014,35(4):71.

[11] Zhao C Y,Zhang Q,Yang C S,et al.Integration of MODIS data and short baseline subset (SBAS) technique for land subsidence monitoring in Datong,China[J].Journal of Geodynamics,2011,52:16-23.

[12] 許領,戴福初.涇陽南塬黃土滑坡特征參數統計分析[J].水文地質工程地質,2008,35(5):28-32.Xu L,Dai F C.Statistical analysis of the characteristic parameters of loess landslides at the South Jingyang Plateau[J].Hydrogeology &Engineering Geology,2008,35(5):28-32.

[13] 李文,龍建輝,李同錄,等.陜西涇陽南塬黃土滑坡滑帶土殘余剪切強度特性[J].地球科學與環境學報,2012,34(4):95-101.Li W,Long J H,Li T L,et al.Residual shear strength characteristics of landslide soil in loess landslide of Jingyang southern highland,Shaanxi[J].Journal of Earth Sciences and Environment,2012,34(4):95-101.

[14] 段釗.黃土滑坡觸發機理研究—以涇河下游南岸黃土塬區滑坡為例[D].西安:長安大學博士學位論文,2013:28-29.Duan Z.Study on the trigger mechanism of loess landslide:A case study of the loess landslide in the south bank of lower Jing River[D].Doctoral dissertation,Xi'an:Chang'an University,2013:28-29.

InSAR monitoring result analysis of land subsidence in loess irrigation area

HE Yang1,ZHAO Chao-Ying1,2,ZHANG Qin1,2
(1.College of Geological Engineering and Surveying,Chang'an University,Shaanxi Xi'an 710054,China;2.Key Laboratory of Western China Mineral Resources and Geological Engineering,Ministry of Education of China,Shaanxi Xi’an 710054,China)

Abstract:Synthetic aperture radar interferometry is a newly developed technology with potential for terrestrial observation from space,and is now being used to obtain information about the earth’s surface.This technology can be used to monitor land subsidence deformation,which is a slow process that can affect large areas and a major geological hazard for cities.We monitored the characteristics of loess irrigated areas using short baseline set technology (SBAS-InSAR technique) with 11 X-band TerraSAR images of the South Jingyang platform.The region of study is located downstream of the Jinghe river’s south bank,about 50 km from Xi’an,Shaanxi.Results demonstrated that the South Jingyang platform has two ground subsidence areas and the accumulated deformation along the profile reached 6.4 cm and 2.7 cm,respectively,in 121 days.We investigated the distribution of water wells,which related to ground subsidence areas.Based on the deformation rate and water level elevation data analysis,we found that land subsidence is caused by irrigation.

Key words:land subsidence;InSAR;small baseline subsets (SBAS) InSAR;time series;loess irrigation area

基金項目:國家重點基礎研究計劃(973計劃)項目(2014CB744703);中國自然科學基金項目(41372375)*通訊作者:趙超英(博士/教授):zhaochaoying@163.com

作者簡介:何楊(1990-),男,碩士生,主要從事地質災害InSAR監測分析.

修訂日期:2015-12-20

收稿日期:2015-11-18

doi:10.3969/j.issn.2095-1329.2016.01.015

中圖分類號：P642.26

文獻標志碼：A

文章編號：2095-1329(2016)01-0066-03

電子郵箱:838913114@qq.com

聯系電話:029-82339251