馬立龍 蔡瑞慶 張帆 李倩文

摘 要：【目的】煤炭資源的開采引發(fā)一系列的生態(tài)環(huán)境問題，地面塌陷問題尤為嚴(yán)重，而地表形變監(jiān)測是修復(fù)生態(tài)環(huán)境問題的首要工作?！痉椒ā縄nSAR技術(shù)作為一種面監(jiān)測技術(shù)，已成功應(yīng)用到煤礦地面塌陷監(jiān)測中，本研究利用雷達(dá)差分干涉測量（D-InSAR）技術(shù)，計(jì)算2020年8月至2022年10月間共31景Sentinel-1A數(shù)據(jù)。【結(jié)果】在監(jiān)測期內(nèi)，開采工作面的地表共形成3個(gè)形變中心，分別位于301和303盤區(qū)的工作面，地表在煤礦開采后進(jìn)入快速形變期，之后開始緩慢形變，直至趨于穩(wěn)定狀態(tài)?！窘Y(jié)論】D-InSAR技術(shù)對于此類地面塌陷問題，相較于點(diǎn)監(jiān)測技術(shù)，具有先天的優(yōu)勢，能夠快速、準(zhǔn)確地獲取整個(gè)工作面的形變情況，對煤礦生態(tài)環(huán)境問題的改善起到了指導(dǎo)作用。

關(guān)鍵詞：楊伙盤煤礦；D-InSAR；Sentinel-1A；地表形變監(jiān)測；工作面

中圖分類號：TD327；P237；P642.26? 文獻(xiàn)標(biāo)志碼：A? 文章編號：1003-5168（2024）03-0037-05

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

Monitoring and Analysis of Surface Deformation in Yanghuopan

Mining Area Based on D-InSAR Technology

MA Lilong1 CAI Ruiqing1 ZHANG Fan2 LI Qianwen2

（1.Xinjiang Bureau of Geo-exploration & Mineral Development the Third Geological Group，Korla 841000，China;2. Xi'an Meihang Remote Sensing Information Co.， Ltd.， Xi'an 710100，China）

Abstract： [Purposes] The mining of coal resources has caused a series of ecological and environmental problems， and among them， the problem of ground collapse is particularly serious， and the monitoring of surface deformation is the primary work to repair ecological environment. [Methods] InSAR technology has been successfully applied to the monitoring of coal mine surface collapse as a surface monitoring technology， and 31 Sentinel-1A data fromAugust 2020 to October 2022 were calculated by using radar differential interferometry （D-InSAR） technology to monitor surface deformation in Yanghuopan Coal mine in this paper. [Findings] During the monitoring period， according to the monthly shape variables obtained by D-InSAR technology， three deformation centers were formed， which were located in the working face of the 301 and 303 panels respectively， and the surface of the coal mine entered the rapid deformation period after mining， and then began to deform slowly until it reached a stable state. [Conclusions] D-InSAR technology has a congenital advantage over point monitoring technology for this kind of ground collapse problem， which can quickly and conveniently obtain the deformation of the whole working face， and plays a reference role in improving the ecological environment monitoring..

Keywords： Yanghuopan Coal Mine; D-InSAR; Sentinel-1A; surface deformation monitoring; working face

0 引言

我國的煤炭等礦產(chǎn)資源較為豐富，但開采的同時(shí)，也會引起較多的生態(tài)環(huán)境問題。例如地面塌陷及其引發(fā)的一系列次生災(zāi)害，導(dǎo)致農(nóng)田被破壞［1］，發(fā)生崩塌或滑坡等［2］，給人民的生產(chǎn)生活帶來了嚴(yán)重的經(jīng)濟(jì)損失，甚至?xí)＜叭嗣裆?cái)產(chǎn)安全［3-4］。因此，在地下煤炭開采過程中對礦區(qū)進(jìn)行地表形變監(jiān)測、分析及預(yù)測顯得尤為重要［5-6］。

目前，在礦區(qū)地表形變監(jiān)測方面，主要還是依賴于傳統(tǒng)的水準(zhǔn)及GPS測量，即在生產(chǎn)工作面上每隔一段距離布設(shè)一個(gè)水準(zhǔn)點(diǎn)或GPS點(diǎn)進(jìn)行測量。此方法一方面人工消耗成本高，周期長；另一方面兩者都是基于點(diǎn)觀測，對于礦區(qū)大范圍、大區(qū)域的監(jiān)測需要大量加密點(diǎn)通過擬合獲取最終結(jié)果，無法保證精度［7-9］。合成孔徑雷達(dá)干涉測量（InSAR）技術(shù)作為一種高效觀測技術(shù)，不受天氣、云霧的影響，具有大范圍面監(jiān)測的優(yōu)勢，很大程度上降低了人力、物力、財(cái)力的消耗［10-12］。

本研究基于D-InSAR技術(shù)，利用中分辨率Sentinel-1A數(shù)據(jù)對楊伙盤煤礦礦區(qū)地表形變進(jìn)行監(jiān)測分析，結(jié)合礦區(qū)的工作面分布等資料，提取礦區(qū)工作面開采沉陷形變場，分析沉陷變化規(guī)律。

1 研究區(qū)概況

楊伙盤煤礦位于陜西省神木市城北約30 km的府店一級公路北側(cè)，行政區(qū)劃屬神木市店塔鎮(zhèn)。工作區(qū)交通便利，西（安）—包（頭）鐵路、210國道西（安）—包（頭）公路從楊伙盤煤礦西側(cè)通過，神（木北）—朔（州）鐵路及府店一級公路均從煤礦南緣黃羊城溝內(nèi)通過，緊鄰煤礦工業(yè)場地。

根據(jù)現(xiàn)有巷道布置情況，結(jié)合各煤層開采范圍及煤層賦存特點(diǎn)，一水平3-1煤劃分為4個(gè)盤區(qū)，分別為301、302、303、304盤區(qū)。3-1號煤層開采工作面分布情況如圖1所示，301盤區(qū)工作面為30101～30120，302盤區(qū)工作面為30201～30208，303盤區(qū)工作面為30301～30304，304盤區(qū)工作面為30401～30405。2020年對30113、30115、30117、30112和30114工作面進(jìn)行回采；2021年對30117、30114、30116、30301和30302工作面進(jìn)行回采；2022年對30302、30303、30116、30118和30119工作面進(jìn)行回采。

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

2.1 SAR數(shù)據(jù)

本研究使用的SAR數(shù)據(jù)為歐洲太空局Sentinel-1A衛(wèi)星在TOPS成像模式下的寬幅干涉SLC數(shù)據(jù)，為目前最常用的中分辨率SAR數(shù)據(jù)，數(shù)據(jù)詳細(xì)參數(shù)見表1。本次共收集的2020年8月至2022年11月間的31景升軌數(shù)據(jù)見表2，軌道號為113/126。

DEM數(shù)據(jù)為ALOS Global Digital Surface Model“ALOS World 3D-30 m”（AW3D30），水平分辨率為30 m（1弧秒），高程精度5 m。其是目前世界上最精確的3D地圖，覆蓋全球所有的土地尺度。精密軌道數(shù)據(jù)采用成像21 d之后發(fā)布的POD精密軌道數(shù)據(jù)。

2.2 D-InSAR技術(shù)方法

D-InSAR技術(shù)是以合成孔徑雷達(dá)兩次及以上數(shù)據(jù)提取的相位信息為信息源，獲取地表三維信息和變化信息的一項(xiàng)技術(shù)。雙軌D-InSAR是對同一研究區(qū)在形變發(fā)生前后獲取的兩幅SAR影像，進(jìn)行干涉得到干涉圖，干涉圖中的相位值即為兩次成像的相位差測量值。然后基于外部的DEM模擬地形相位，與干涉圖進(jìn)行差分處理，得到差分干涉圖，最終提取目標(biāo)點(diǎn)處地表沿視線向形變量，流程如圖2所示。

D-InSAR技術(shù)主要是利用相位差信息來獲取地形及形變信息［13］，相位差信息計(jì)算見式（1）。

[δ?=?1??2=4πλ（ρ1?ρ2）] （1）

式中：[δ?]表示相位差；[?]表示相位；[λ]表示波長；[ρ]表示斜距。

將衛(wèi)星分別在兩個(gè)軌道處獲取的SAR圖像做相干運(yùn)算，獲取形變前的地形差分相位信息，同理可以獲取形變發(fā)生后的地形差分相位信息。將形變前后的地形差分相位信息作差分，就可以獲取地表形變的干涉紋圖。

3 結(jié)果與分析

通過D-InSAR技術(shù)獲取工作面的30幅月度形變量圖。由于工作面采煤的時(shí)間和強(qiáng)度情況不同，每個(gè)時(shí)間段的形變分布范圍和量級均有差異，總體來說，形變主要分布在①、②、③三個(gè)區(qū)域，分別位于301和303盤區(qū)。例如，20210226—20210403時(shí)間段的形變量分布如圖3所示。由圖3可知，①、②區(qū)域均有形變，而③區(qū)域的工作面尚未進(jìn)行開采，因此無形變分布。在①、②、③三個(gè)區(qū)域的工作面分別做一條剖線，形變分析如下。

①區(qū)域的月度形變分布如圖4所示，工作面開采及形變時(shí)間間隔基本在2020年8月至2022年10月之間。涉及的工作面為30112、30114、30116、30118。本研究對30114工作面的開采情況進(jìn)行說明。在2020年10月左右［圖4（4）］，30114工作面開始開采，到2021年9月末［圖4（14）］，開采導(dǎo)致的形變基本結(jié)束。在30114工作面上畫一條剖線，分析其在開采階段的形變情況，如圖5所示。在開采過后，地表沿著工作面隨著時(shí)間推移逐段出現(xiàn)突變，最大形變量約為-40 mm，處于20210521—20210801時(shí)間段，主要是由于時(shí)間間隔較長。其余時(shí)間段的形變量級基本都在-30 mm范圍內(nèi)。

②區(qū)域的月度形變分布情況如圖6所示，工作面開采及形變時(shí)間間隔為2020年8月至2021年5月，之后該區(qū)域的工作面地表無形變。開采的工作面主要有30115和30117工作面，而形變覆蓋的工作面為30113、30115、30117和30119。本次對30117工作面的開采情況進(jìn)行說明。大概在2020年11月左右［圖6（5）］，30117工作面開始開采，到2021年5月末［圖6（10）］，形變基本結(jié)束。在30117

工作面上畫一條剖線分析其在開采階段的形變情況，如圖7所示。在開采過后，地表沿著工作面逐段出現(xiàn)突變，最大形變量約為-35 mm。

③區(qū)域的月度形變分布情況如圖8所示，工作面開采及形變時(shí)間間隔為2021年5月至2022年10月。開采和形變涉及的工作面主要有30301、30302和30303工作面。本研究對30302工作面的開采情況進(jìn)行說明。2021年10月［圖8（6）］，30302工作面開始開采，到2022年5月［圖8（13）］，形變結(jié)束。在30302工作面上畫一條剖線，分析其在開采階段的形變情況，如圖9所示。在開采過后，地表沿著工作面逐段出現(xiàn)突變，最大形變量約為-45 mm。

4 結(jié)論

本研究利用傳統(tǒng)的D-InSAR技術(shù)對楊伙盤煤礦地表進(jìn)行了月度形變量監(jiān)測，并結(jié)合工作面的開采情況進(jìn)行了具體分析，得出以下結(jié)論。

①獲取到2020年8月至2022年11月之間的31幅形變量分布圖，形變隨著開采工作面的采煤推進(jìn)進(jìn)行展布，可探測的最大形變量約為-63 mm。

②在監(jiān)測期內(nèi)，煤礦開采的工作面為30112、30114、30116、30118、30115、30117、30301、30302、30303。而形變覆蓋的工作面比開采工作面范圍廣，30113和30119工作面也有形變分布。

參考文獻(xiàn)：

［1］黃忠倫.煤礦開采對土地資源的破壞及對策研究［J］. 內(nèi)蒙古煤炭經(jīng)濟(jì)， 2018（1）：92，94.

［2］王志平.煤礦地質(zhì)災(zāi)害特征及防治措施［J］.中國石油和化工標(biāo)準(zhǔn)與質(zhì)量， 2019， 39（20）：140-141.

［3］劉貴.采煤沉陷區(qū)地面塌陷原因及塌陷坑附近房屋居住安全分析［J］.煤炭工程，2019，51（8）：145-148.

［4］賀衛(wèi)中，向茂西，劉海南，等.榆神府礦區(qū)地面塌陷特征及環(huán)境問題［J］.煤田地質(zhì)與勘探，2016，44（5）：131-135.

［5］劉曉玲，魏奧林，王毅，等.淺析陜北煤礦礦區(qū)地質(zhì)災(zāi)害發(fā)育特征及其成災(zāi)過程［J］.中國地質(zhì)災(zāi)害與防治學(xué)報(bào)，2016，27（4）：70-73.

［6］劉善軍，吳立新，毛亞純，等.天-空-地協(xié)同的露天礦邊坡智能監(jiān)測技術(shù)及典型應(yīng)用［J］.煤炭學(xué)報(bào)，2020，45（6）： 2265-2276.

［7］王志勇，張繼賢，黃國滿.基于InSAR的濟(jì)寧礦區(qū)沉降精細(xì)化監(jiān)測與分析［J］.中國礦業(yè)大學(xué)學(xué)報(bào)，2014，43（1）：169-174.

［8］黃寶偉.基于D-InSAR和GIS技術(shù)的煤礦區(qū)地面沉降監(jiān)測研究［D］.青島：中國石油大學(xué)，2011.

［9］麻源源， 陳云波， 左小清， 等. Sentinel-1A數(shù)據(jù)及短基線集的昆明地面沉降分析［J］.測繪科學(xué)，2019，44（11）：59-66，95.

［10］孫曉鵬，魯小丫，文學(xué)虎， 等.基于SBAS-InSAR的成都平原地面沉降監(jiān)測［J］.國土資源遙感，2016， 28（3）： 123-129.

［11］RABUS B，EINEDER M，ROTH A ，et al.The shuttle radar topography mission：a new class of digital elevation models acquired by spaceborne radar［J］.ISPRS Journal of Photogrammetry and Remote Sensing， 2003， 57（ 4）：241-262.

［12］ROSSI C，GONZALEZ F R，F(xiàn)RITZ T，et al.TanDEM-X calibrated Raw DEM generation［J］.Isprs Journal of Photogrammetry & Remote Sensing， 2012， 73（SEP.）：12-20.

［13］查顯杰，傅容珊，戴志陽.DInSAR技術(shù)對不同方位形變的敏感性研究［J］.測繪學(xué)報(bào)，2006（2）：133-137.