基于D-InSAR監測的黃土礦區地表沉陷特征分析

王宗良 湯伏全

摘 要：合成孔徑雷達差分干涉測量（D-InSAR）技術是20世紀90年代發展起來的微波遙感技術，具有高精度地監測地表微小變形的能力，已經成為礦區地表沉陷監測的有效技術手段。以陜西彬長礦區為研究區，利用多期SAR影像數據，應用D-InSAR技術方法開展地表沉陷監測實驗，并結合礦區的實際情況，對礦區下沉等值線圖進行簡要分析。實驗發現3種礦區工作面均發生了不同程度的沉陷變形，其中沉降量最大值出現在亭南煤礦304工作面內，達到17.5 cm.結果表明：D-InSAR技術可以反映出礦區沉陷位置和下沉程度，展示出D-InSAR技術在黃土礦區地表沉陷監測方面具有良好實用性和廣泛應用前景。

關鍵詞：差分干涉測量;黃土礦區;沉陷變形;形變監測

中圖分類號：X 936?? 文獻標識碼：A?? 文章編號：1672-7312（2019）03-0394-04

Abstract：Synthetic Aperture Radar Differential Interferometry （D-InSAR） technology is a microwave remote sensing technology developed in the 1990s. It has the ability to monitor surface micro-deformation with high precision and has become an effective technical means for monitoring surface subsidence in mining areas. Taking the Binchang mining area in Shaanxi as the research area， using the multi-period SAR image data， the D-InSAR technology method was used to carry out the surface subsidence monitoring experiment， and combined with the actual situation of the mining area， the sinking contour map of the mining area was briefly analyzed. The experiment found that the working faces of the three mining areas have experienced different degrees of subsidence deformation， and the maximum settlement amount appears in the 304 working face of the Tingnan Coal Mine， reaching 17.5 cm.The results show that D-InSAR technology can reflect the location and subsidence of mining area， and demonstrates that D-InSAR technology has good practicability and wide application prospects in surface subsidence monitoring in loess mining area.

Key words：differential interferometry;loess mining area;subsidence deformation;deformation monitoring

0 引 言

煤礦資源的開發利用在社會和生活中發揮越來越重要的作用[1]。以陜西彬縣為例，彬長煤田是陜西省第二大煤田，具有煤層厚、儲量大、易開采的特點。礦區地表沉降監測主要采用建立觀測站，采用GPS測量、全站儀測量、精密水準測量等傳統測量方法進行定期觀測，存在著費時費力、效率低下等缺陷。因此，應用新的監測技術已經成為礦山開采沉陷監測的必然趨勢。差分干涉測量（D-InSAR）技術作為一種新型地面觀測技術，它是通過2個回波信號的相位差與衛星軌道參數來反演出地球表面的微小變形，并且能夠在短周期內提供空間連續表面變形信息，可以獲得厘米甚至毫米的高精度表面變形信息[2]，可以與常規地面測量技術起到互補作用。

國內外研究學者相繼利用雷達遙感技術監測地表沉降。1989年，Garbriel等人利用D-InSAR技術實現了地表形變監測，監測結果達到厘米級精度[3];2009年，Z.Perski等人利用D-InSAR技術監測波蘭的Wieliczka鹽礦，成功獲得到該區域的地表形變[4];2012年，郭炳躍等人利用D-InSAR技術提取徐州礦區地表形變[5];洪卓眾等人采用D-InSAR技術獲得陜北神木柳塔礦區的地表沉降[6];魏長靖等人采用D-InSAR技術獲取安徽省錢營孜煤礦的地表沉降圖[7]。這些研究表明礦區周圍通常受到開采的影響，礦區地表周圍均會發生不同程度的沉陷。因此，有效監測并獲取礦區地表形變信息具有十分重要的意義。

1 D-InSAR的基本原理

差分干涉測量（D-InSAR）是指使用2幅干涉圖像相同的區域，其中一幅是變形前獲取的干涉圖像，另一幅是變形后獲取的干涉圖像，通過差分處理來獲取地表形變的測量技術[8]。根據采用干涉數據的多少和處理方法的差異，一般可將D-InSAR分為2-pass型、3-pass型和4-pass型[9-10]。實驗采用的D-InSAR方法處理流程圖如圖1所示。

2 實驗區與影像數據選取

2.1 礦區概況

彬長礦區位于陜西省關中西北部的長武和彬縣境內，是國家規劃的13個煤炭基地之一的黃隴基地的主要礦區，屬于黃隴侏羅紀煤田的中斷[11]。礦區的南部和東部主要以無煤區作為邊界線，北部和西部以陜甘省界作為邊界線，礦區東西長約為46 km，南北寬約為36.5 km，規劃面積約為978 km2，煤炭資源地質儲量約為89億噸，可采儲量約54萬噸，主采煤層平均厚度10.6 m.礦區地處涇河流域，地表水與地下水資料較豐富，地表植被較少，自然村落大部分處在開采沉陷區域以外的保護煤柱范圍，這一特征有利于SAR影像的數據處理。

2.2 SAR影像的選取

對于SAR數據，不僅包含了在接收信號時的幅度信息，同時也包含了它的相位信息[12]。D-InSAR正是利用相位信息獲得目標的高程信息，進而獲得形變信息。在進行差分干涉處理之前，就要滿足影像之間具備良好的相干性，其中，合理的基線估計是在處理數據前必須要完成的一項內容。研究區域彬長礦區處于黃土高原過渡地帶，西北黃土地區具有土壤含水量小、氣候條件干燥、植被較為稀疏等特征，為后續處理工作提供了良好的數據基礎。因此，利用彬縣行政區多景Sentinel-1A衛星的TOPSAR雷達數據（C-band）作為實驗的干涉影像數據，采用的4個具體干涉對影像的數據參數見表1.

2.3 干涉處理軟件簡介

在SAR數據的處理方面，世界各國研究學者普遍使用的軟件，如ERDAS Imagine，ENVI，PCL Geomatica，ER Mapper等，實驗選用美國研發的ENVI（Version 5.3.1）軟件。ENVI中的遙感圖像處理工具ENVI SARscape由核心模塊及5個擴展模塊構成，支持一系列機載和星載雷達系統的數據，包括Sentinel-1A，ERS-1/2，JERS-1，ENVISAT ASAR， RADARSAT-1/2等，通過對SAR數據在ENVI SARscape中進行處理后，能夠便捷地提取到高精度地形數據（DEM）和地表形變圖，還可以在ArcGis中將提取的這些有用的數據和形變圖與光學遙感數據等集成起來并疊加分析，全面提升SAR數據應用價值。實驗采用SARscape模塊中的D-InSAR工作流對哨兵Sentinel-1A衛星的TOPSAR雷達數據（C-band）進行處理，確定礦區位置并得到下沉等值線圖。

3 D-InSAR數據處理

根據D-InSAR的基本原理，對SAR數據進行處理，流程主要包括：基線估算、干涉圖生成及去平、自適應濾波、相位解纏、相位轉形變與地理編碼等諸多步驟。

在進行差分干涉處理前，基線估算的結果是判斷SAR干涉測量成像能否進行下去的基礎條件，同時是解釋SAR影像對失相干的有力佐證。實驗數據的基線估算結果表明影像干涉對可以進行后續的差分干涉處理。干涉圖的生成是通過對主、幅影像對應像元的復數值進行共軛相乘得到的，干涉圖中的相位值就是相位差的主值[13]。在生成干涉圖的過程中，存在數據處理噪聲、系統熱噪聲等影響，為了能順利完成相位解纏等后續步驟，就要對噪聲加以抑制，從而降低雷達噪聲的影響，這就需要對干涉圖進行濾波處理[14]，實驗選擇的Goldstein濾波是最常用的濾波方法，能在一定程度上提高干涉條紋的清晰度以及減少由于時間基線或空間基線造成的失相干問題。

其中較為關鍵的是相位解纏，它是獲取地表微小形變量的關鍵環節，解纏結果的好壞直接影響著最終高程數據的精度[15]。考慮到研究區周邊地質地貌條件較為復雜的情況下，文中選用最小費用流法進行相位解纏，該算法的好處在于當影像上存在有大面積的低相干部分時，或者是有其他限制增長因素導致解纏困難時，可以得到比其他解纏算法更加可靠、準確的結果。

地理編碼是D-InSAR工作流處理中的最后一步，它是將雷達的斜距坐標轉變為地理坐標的過程，經過相位轉形變和地理編碼，得到精度較高的數字高程模型。

4 結果分析

4.1 SAR沉降區與礦區坐標系的統一

SAR圖像做干涉處理是在雷達影像的坐標系下完成的，經過地理編碼獲取最終的DEM以及研究區形變圖是WGS-84坐標系下的坐標;而礦區生產及地表沉陷監測使用的是以X，Y表示的平面坐標。因此，通過七參數轉化法，將亭南煤礦304工作面、燕家河煤礦8211工作面、下溝煤礦ZF2404工作面由1980西安坐標系轉化到WGS-84坐標系下，投影方式設為高斯投影3°帶，根據彬長礦區地理位置，中央子午線設為108°.

4.2 不同煤礦沉降分析

根據D-InSAR數據處理流程，對獲取的形變圖與實地工作面進行坐標之間的統一，找到了不同煤礦的工作面（圖2中用黑色線框表示）在形變圖中不同沉降區域的具體位置。通過使用ArcMap 10.1軟件，對形變圖進行疊加，并生成不同煤礦的等值線圖（均以1 cm沉降量為邊界）時間跨度為2017年12月10日至2018年03月28日，圖2中的（a）、（b）、（c）分別代表亭南煤礦、燕家河煤礦和下溝煤礦。

從圖2中可以看出，在監測時段內，亭南煤礦304工作面（即（a）中的黑色邊框）以及附近地表均有不同程度的沉降，隨著時間推進，地表沉降量在逐漸增大，范圍由工作面向四周不斷延伸，并且越接近工作面中間，沉降量越大，最大約達到17.5 cm;燕家河煤礦8211工作面（即（b）中的黑色邊框）以及附近地表均有不同程度的沉降，隨著開采的進行，地表沉降量在逐漸增大，范圍主要沿開采工作面方向，即東西方向不斷延伸，并且越接近工作面，沉降量越大，最大約達到13.8 cm;下溝煤礦ZF2404工作面（即（c）中的黑色邊框）以及附近地表均有不同程度的沉降，隨著煤礦沿南北方向開采的不斷進行，地表沉降量在逐漸增大，范圍也在南北方向有著明顯的延伸，并且越接近工作面，沉降量越大，最大約達到10.3 cm.

5 結 論

1）以陜西彬長礦區的開采沉陷監測為研究對象，基于D-InSAR技術進行了地表沉陷監測數據處理的實驗研究，確定了適用于西部黃土礦區D-InSAR監測的影像數據處理方法與實用技術流程。

2）采用ENVI軟件針對實驗區域的干涉SAR影像對進行差分干涉處理。結果表明，監測時間段內，彬長礦區多個區域發生了明顯的沉陷變形。其中有3個變形區分別位于亭南煤礦、燕家河煤礦、下溝煤礦相應開采工作面上方地表，將上述3個變形區域進行下沉疊加并繪出等值線圖，在同一坐標系統下與3個煤礦的開采工作面位置進行分析。結果表明：隨著開采的進行，3個不同煤礦的工作面均發生了沉陷變形，其變形發展趨勢以及最大沉陷量與實際資料基本一致，從而驗證了D-InSAR技術在監測礦區沉陷變形方面的有效性。

綜上，D-InSAR技術為黃土礦區地表沉陷的監測提供了極具前景的技術手段，未來將吸引更多的研究學者開展廣泛而又深入的研究。

參考文獻：

[1] 姜 巖，高均海.合成孔徑雷達干涉測量技術在礦山開采地表沉陷監測中的應用[J].礦山測量，2003（01）：5-7.

[2]Massonnet D，Holzer T，Vadon H.Correction to《land subsidence caused by the east mesa geothermal gield，California，observed using SAR interferometry》[J].Geophysical Research Letters，2013，24（08）：3213-3213.

[3]Gabriel A K，Goldstein R M，Zebker H A.Mapping small elevation changes over large areas：Differential radar interferometry[J].Journal of Geophysical Research Solid Earth，1989，94（B7）：9183-9191.

[4]Perski Z，Hanssen R，Wojcik A，et al.InSAR analyses of terrain deformation near the Wieliczka Salt Mine，Poland[J].Engineering Geology，2009，106（01）：58-67.

[5]郭炳躍，何 敏，劉建東.利用InSAR技術監測徐州市礦區地表變形[J].地質學刊，2012，36（01）：99-103.

[6]洪卓眾.基于D-InSAR的礦區沉降監測[J].山西建筑，2012，38（22）：73-75.

[7]魏長靖，閏建偉.D-InSAR技術二軌法監測礦區地表沉陷的方法研究[J].煤炭技術，2012，31（07）：129-130.

[8]成 樞，劉 元，張洪勝.基于D-InSAR的開采沉陷監測技術[J].測繪與空間地理信息，2013，36（06）：7-9.

[9]王 超，張 紅，于 勇，等.雷達差分干涉測量[J].地理學與國土研究，2002，18（03）：13-17.

[10]田 馨，廖明生.InSAR技術在監測形變中的干涉條件分析[J].地球物理學報，2013，56（03）：812-823.

[11]姚頑強，馬 飛，龔 云，等.基于D-InSAR的彬長礦區沉陷變形監測[J].西安科技大學學報，2011，31（05）：563-568.

[12]王 超，張 紅，劉 智.星載合成孔徑雷達干涉測量[M].北京：科學出版社，2002.

[13]董玉森，葛玲玲，張新春，等.基于差分雷達干涉測量的礦區地面沉降監測研究[J].武漢大學學報：信息科學版，2007，32（10）：888-892.

[14]何秀風，何敏.InSAR對地觀測數據處理方法與綜合測量[M].北京：科學出版社，2012.

[15]舒 寧.雷達影像干涉測量原理[M].武漢：武漢大學出版社，2003.

（責任編輯：張 江）