楊澤鵬 文泰翔

摘 要：【目的】獲取桂林市興安縣羊角山區(qū)域的地面形變監(jiān)測結(jié)果，探究人類活動對地面形變的影響?！痉椒ā坎捎?017年3月17日—2021年12月9日覆蓋羊角山地區(qū)的90景升軌Sentinel-1A衛(wèi)星影像作為數(shù)據(jù)源，基于SBAS-InSAR方法對該區(qū)域進行地面形變監(jiān)測分析，SBAS-InSAR方法能夠有效地減弱時空失相干擾和大氣效應引起的誤差，提高地面形變監(jiān)測結(jié)果的準確性?！窘Y(jié)果】在長時間序列的小尺度地表開發(fā)形變監(jiān)測上使用SBAS-InSAR方法仍能得到較高精度的結(jié)果，研究區(qū)域地表平均形變速率處于-18 mm/a～10 mm/a，在研究區(qū)的西南方向累積形變量最大，達到了-60 mm。【結(jié)論】該區(qū)域由于開發(fā)時間長度和開發(fā)強度的不同表現(xiàn)出不同的形變結(jié)果，且隨著對該區(qū)域的持續(xù)開發(fā)，整體形變量有繼續(xù)增加的趨勢。

關(guān)鍵詞：SBAS-InSAR；地面形變；形變分析

中圖分類號：P237；TU196.1? ? ? ? ? ? ?文獻標志碼：A? ? ? ? ? ? ? 文章編號：1003-5168（2023）06-0113-04

DOI：10.19968/j.cnki.hnkj.1003-5168.2023.06.022

Analysis of Ground Deformation Monitoring in Yangjiaoshan Mining Area Based on SBAS-InSAR Technology

YANG Zepeng? ?WEN Taixiang

（Guilin University of Electronic Technology，Guilin 541000，China）

Abstract： [Purposes] To obtain the monitoring results of ground deformation in Yangjiao Mountain， Xing'an County， Guilin City， and show the impact of human activities on ground deformation. [Methods] In this paper， 90 satellite images of Satinel-1A， which covered the Yangjiao Mountain area from March 17， 2017 to December 9， 2021， are used as the data source to monitor and analyze the ground deformation in the area based on the SBAS-InSAR method， which can effectively reduce the errors caused by spatiotemporal incoherence and atmospheric effects， and improve the accuracy of the ground deformation monitoring results. [Findings] The SBAS-InSAR method can still obtain high precision results on the small-scale surface development deformation monitoring of long time series. The average surface deformation rate in the study area is between -18 mm/a and 10 mm/a， and the maximum cumulative deformation in the southwest of the study area is -60 mm. [Conclusions] The region shows different deformation results due to different development time length and development intensity， and with the continuous development of the region， the overall deformation variable continues to increase.

Keywords： SBAS-InSAR; ground deformation; deformation analysis

0 引言

近年來，出于經(jīng)濟效益人們加大了對山區(qū)資源的開發(fā)，在開發(fā)過程中容易改變山區(qū)原有的地形地貌，打破山區(qū)原有的土地受力平衡，引發(fā)地質(zhì)災害［1］。因此需要及時了解地面形變信息，從而降低地質(zhì)災害產(chǎn)生的破壞。

InSAR技術(shù)具有高分辨率和強穿透力等特性，不受氣候和地形因素的限制。然而，長時間的監(jiān)測過程容易受到時空失相干和大氣延遲效應等問題的影響［2］。小基線集干涉測量（Small Baseline Subset InSAR， SBAS-InSAR）以D-InSAR技術(shù)為基礎(chǔ)［3］，通過設(shè)置時間基線和空間基線的閾值，將多景SAR數(shù)據(jù)圖像分成若干個干涉子集，獲得受時空失相干較小影響的監(jiān)測點目標［4］，從而得到整個研究區(qū)的形變結(jié)果。

針對在小尺度區(qū)域的地表開發(fā)形變監(jiān)測中存在的相干性低、相干點位少的問題，本研究采用 SBAS－InSAR技術(shù)對桂林市興安縣羊角山2017—2021年的小尺度區(qū)域地面進行形變監(jiān)測，在時間序列上對其地面形變監(jiān)測進行精細化處理，結(jié)合該區(qū)域周圍工業(yè)信息的空間位置分布信息，分析了羊角山地面形變與人類活動（如工業(yè)化發(fā)展和土地負荷增加）之間的相互關(guān)系。

1 研究區(qū)概況

興安縣隸屬于廣西壯族自治區(qū)桂林市，位于廣西東北部，地處東經(jīng)110°14′～110°56′，北緯25°17′～25°55′之間，境內(nèi)地貌多樣且復雜，地勢向西南傾斜。研究區(qū)羊角山距興安縣城3 km，是一個石灰?guī)r開發(fā)山區(qū)，面積約為1.8 km2，周圍分布有大圓新村、丁家村、桂興村。在對該區(qū)域長期的開發(fā)過程中，使其地面產(chǎn)生了不同程度的形變。研究區(qū)域如圖1所示。

2 研究數(shù)據(jù)和方法

2.1 研究數(shù)據(jù)

本研究試驗數(shù)據(jù)選取覆蓋羊角山區(qū)域的90景升軌Sentinel-1A衛(wèi)星影像作為數(shù)據(jù)源，模式影像為IW，極化方式為VV，時間跨度從2017年3月17日—2021年12月9日。為了消除地形誤差對試驗結(jié)果的影響，以美國航空航天局提供的30 m分辨率的SRTM1數(shù)據(jù)作為外部數(shù)字參考高程模型（DEM）。使用歐洲航空局提供的衛(wèi)星精密軌道數(shù)據(jù)（POD）進行軌道糾正，提高影像配準度和基線估算的準確性［5］。

2.2 SBAS-InSAR方法原理

SBAS-InSAR方法是由Beradino等人提出的一種以InSAR技術(shù)為基礎(chǔ)的時間序列分析法［6］。通過設(shè)定適當?shù)臅r間基線與空間基線閾值，生成若干個小基線集合，采用奇異值分解（SingularValue Decomposition，SVD）［7］的方法求出最小二乘解，獲得目標點的平均形變速率。

將研究區(qū)域內(nèi)N+1幅SAR影像按t0，…，tN的時間進行排序，將其中的一幅作為主影像，其余的輔影像與主影像進行配準，配準得到M幅干涉圖，其中M需滿足式（1）條件［8］：

分別在tA和tB（tA>tB）時刻獲取的主影像和輔影像生成的第j幅差分干涉圖，該圖中任一像元的干涉相位可以表示為式（2）：

式中：j∈（1，…，M）；λ為雷達信號波長；φA（x，r）和φB（x，r）表示點（x，r）在tA和tB時刻相對于t0時刻的相位，dA（x，r）和dB（x，r）表示點（x，r）在tA和tB時刻相對于t0時刻的累積形變量；?φ[jtopo]（x，r）為該點的地形相位誤差；?φ[jatmo]（x，r）為大氣延遲相位誤差；?φ[jnoise]（x，r）為由噪聲引起的相位誤差。若不考慮上述誤差的影響，則式（2）可以簡化為式（3）。

為了表示沉降序列中的物理意義，可用式（3）中相位與時間的關(guān)系來表示相位的平均速度為式（4）。

則該第j幅差分干涉圖的相位值可表示為式（5）。

其矩陣表示為式（6）。

式（6）中，A是一個M×N矩陣，φ為SAR影像中相干點的相位值向量，δφ為解纏后的干涉圖中相位值向量，由于多個集合在基線組合時會出現(xiàn)矩陣A的秩虧，根據(jù)最小二乘法此時方程解的個數(shù)為無窮多?？刹捎闷娈愔捣纸夥ㄇ蠼饩仃嘇的最小二乘解，得出形變速率。

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

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

試驗主要流程如下。

①SAR影像干涉對連接圖生成。試驗以2020年3月1日的影像為主影像，將其他日期的輔影像與主影像配準進行處理，得到422對干涉對。

②干涉圖的處理。試驗中以30 m分辨率的SRTM1數(shù)據(jù)作為外部DEM，消除地形相位。多視系數(shù)比為4：1，采用Delaunay MCF的方法進行相位解纏，濾波方法為Goldstein，濾波窗口最小為16，濾波窗口最大Alpha值為2.5。

③軌道精煉和重去平。選擇相干性高、相位好的點作為地面控制點，這些點將作為無形變的點來修正影像數(shù)據(jù)，去除殘余的相位。經(jīng)過多次試驗，利用相干性閾值系數(shù)法共選取了34個地面控制點。

④反演計算生成結(jié)果。反演計算估計位移速率和殘余地形，并對生成的干涉圖去平，重做相位解纏，提高干涉圖的質(zhì)量。計算時間序列上的平均位移速率，利用濾波方法估算和去除大氣相位，獲得更加準確的研究區(qū)域形變速率結(jié)果。

3.2 結(jié)果分析

羊角山2017年3月—2021年12月的平均形變速度，如圖2所示。由圖2可知：該研究區(qū)域形變較為嚴重，區(qū)域范圍內(nèi)地面平均形變速率處于-18 mm/a～10 mm/a。經(jīng)查近年來對該區(qū)域的開發(fā)活動比較頻繁，該區(qū)域地面存在多處不同程度的形變。為了得到研究時間段內(nèi)該區(qū)域具體的累積形變時序，在形變較為嚴重的區(qū)域分別選擇了三個監(jiān)測點進行累計形變量的時序分析，分別為點 A1、A2和A3。

對三個監(jiān)測點分別進行時序分析得到如圖3所示的時序形變圖，從圖3中結(jié)果可得，在監(jiān)測點A1，該點在2017年3月—2018年12月期間形變量表現(xiàn)出連續(xù)下降的趨勢，但從2019年開始累積形變量基本維持在-40 mm上下。

在2017年3月17日—2021年12月9日期間，A2、A3監(jiān)測點的最大形變量分別為-30 mm和-60 mm，且隨著時間的變化，他們的累積形變量的變化趨勢基本表現(xiàn)為連續(xù)下降，由此判斷這兩個區(qū)域的地表未來依然會繼續(xù)下降。

4 結(jié)論

本研究采用2017—2021年90景Sentinel-1A數(shù)據(jù)，利用SBAS-InSAR技術(shù)獲取了桂林市興安縣羊角山區(qū)域的地面形變監(jiān)測結(jié)果，分析了形變區(qū)域監(jiān)測點的累積形變量，表明了人類活動對地面形變的影響，得出如下結(jié)論。

①研究區(qū)中形變現(xiàn)象最嚴重的區(qū)域在羊角山西南部，在形變現(xiàn)象最為突出的區(qū)域選擇監(jiān)測點進行研究，分析表明此點在研究時段內(nèi)的累積形變量最大為-60 mm，且判斷該監(jiān)測區(qū)域以后還存在持續(xù)下降的趨勢。

②桂興村處于研究區(qū)羊角山西南方，隨著對羊角山開發(fā)范圍的不斷擴大，已經(jīng)觸及了桂興村，由此而帶來的地表形變可能會導致泥石流等災害的發(fā)生，影響當?shù)卮迕竦某鲂邪踩腿粘Ｉ?。因此有關(guān)部門有必要加強對該區(qū)域的監(jiān)測并對后續(xù)的地表形變保持觀測，有關(guān)企業(yè)也需制定合適的開發(fā)計劃，防止因不合理的開發(fā)而產(chǎn)生的地表形變對周圍村民的生活造成影響。

參考文獻：

［1］師蕓，李杰，呂杰，等. 結(jié)合SBAS-InSAR與支持向量回歸的開采沉陷監(jiān)測與預測［J］. 遙感信息，2021，36（2）：6-12.

［2］何秀鳳，高壯，肖儒雅，等. 多時相Sentinel-1A InSAR的連鹽高鐵沉降監(jiān)測分析［J］.測繪學報，2021，50（5）：600-611.

［3］姜德才，張永紅，張繼賢，等.天津市地鐵線不均勻地表沉降InSAR監(jiān)測［J］.遙感信息，2017，32（6）：27-32.

［4］陸好健，李金平，邵九明，等.基于SBAS-InSAR的紅河斷裂帶南段形變特征研究［J］.測繪工程，2018，27（9）：16-20，25.

［5］LANARI R，MORA O，MANUNTA M，et al.A small-baseline approach for investigating? deformations on full-resolution differential SAR interferograms［J］. Geoscience & Remote Sensing IEEE Transactions，2004，42（7）：1377-1386.

［6］LANARI R，LUNDGREN P，MANZO M，et al.Satellite radar interferometry time series analysis of surface deformation for Los Angeles，California［J］.Geophysical Research Letters，2004，31（23）：345-357.

［7］莫瑩，江利明，孫奇石，等.杭州灣上虞地區(qū)SBAS-InSAR地面沉降監(jiān)測分析［J］. 測繪科學，2020，45 （10）：77-84.

［8］石固林，徐浪，張璇鈺，等. 西山村滑坡時序形變的SBAS-InSAR監(jiān)測［J］. 測繪科學，2021，46（2）：93-98，105.