彭鵬

摘要：本文以安徽省淮南煤礦區(qū)為例，編程訂購10期RadarSAT-2雷達(dá)影像數(shù)據(jù)，采用D-InSAR技術(shù)，獲取了采空塌陷的范圍與沉降速率。結(jié)果表明：1.D-InSAR技術(shù)可以有效地識別出光學(xué)遙感未能識別的潛在或未穩(wěn)沉的塌陷區(qū)，經(jīng)實(shí)地調(diào)查，結(jié)果真實(shí)可靠，與野外具有一致性。2.可以定量估算出采空塌陷區(qū)的沉降速率，本區(qū)最大沉降速率為4.9cm/24天。3.InSAR技術(shù)在地下開采監(jiān)測方面與傳統(tǒng)的測量技術(shù)相比，表現(xiàn)出極大的優(yōu)勢，具有快速、準(zhǔn)確、大范圍、成本低的特點(diǎn)，應(yīng)用前景廣闊。

關(guān)鍵詞：采空塌陷；InSAR；遙感RadarSAT-2

前言

礦產(chǎn)資源的開采造成的采空塌陷對生態(tài)環(huán)境造成一定影響，成為嚴(yán)重制約礦區(qū)可持續(xù)發(fā)展的重要因素，它不僅破壞土地資源、導(dǎo)致生態(tài)環(huán)境惡化，而且影響人民的生產(chǎn)、房屋生活設(shè)施，進(jìn)而誘發(fā)一系列社會、經(jīng)濟(jì)問題（圖1、圖2）。當(dāng)前，生態(tài)文明建設(shè)已成為社會關(guān)注的熱點(diǎn)，如何對生態(tài)環(huán)境進(jìn)行恢復(fù)，如何快速的識別采空塌陷區(qū)的范圍和沉降速率，成為治理和監(jiān)測的首要目標(biāo)。常規(guī)的水準(zhǔn)測量、GPS測量監(jiān)測礦區(qū)采空塌陷的技術(shù)存在監(jiān)測周期長、成本高、無法達(dá)到區(qū)域全覆蓋監(jiān)測等問題。

近些年發(fā)展起來的合成孔徑雷達(dá)干涉測量技術(shù)（Inter-ferometric SAR，InSAR）具有全天時、全天候的特點(diǎn)，可以從空間直接獲取大范圍、高精度的形變信息、高效可重復(fù)的動態(tài)監(jiān)測，廣泛地應(yīng)用到地表形變、采空塌陷、地震形變、冰川移動、火山運(yùn)動以及山體滑坡等方面。InSAR技術(shù)可以快速、準(zhǔn)確、省時的監(jiān)測到塌陷區(qū)的動態(tài)變化，節(jié)省大量資金，和傳統(tǒng)的測量方式相比，有很大的優(yōu)越性（表1）。本文以淮南煤礦區(qū)為例，采用二軌法對研究區(qū)采空塌陷進(jìn)行信息提取與分析，為掌握和治理采空塌陷區(qū)提供基礎(chǔ)資料與技術(shù)支持。

1.研究區(qū)概況和以往形變概況

研究區(qū)位于安徽省淮南市，坐標(biāo)范圍為東經(jīng)114°50′～118°10′、北緯32°20′～34°40′之間。全市總面積2121km²，市轄田家庵、大通、謝家集、八公山、潘集五區(qū)和鳳臺縣以及毛集綜合實(shí)驗(yàn)區(qū)，總?cè)丝?10萬人?；茨鲜侵袊宕竺禾镏??；茨夏壳耙寻l(fā)展成為以煤炭、電力為主體，兼有紡織、電子等門類較為齊全的新型工業(yè)城市。淮南煤礦區(qū)長期礦產(chǎn)資源開采給社會帶來巨大效益的同時也嚴(yán)重地破壞了周邊的環(huán)境，采空塌陷就是直接表現(xiàn)，已成為影響和制約礦業(yè)城市可持續(xù)發(fā)展的重要因素。至2002年，全市范圍內(nèi)采煤塌陷區(qū)面積已近50km²。其中，形成水面的面積為13km²，占塌陷總面積的30%左右，主要分布在大通、謝家集、八公山、潘集4個區(qū)和鳳臺縣。

2.數(shù)據(jù)源

采用的雷達(dá)數(shù)據(jù)為2012年～2013年編程訂購加拿大MDA公司的RADARSAT-2（20～30m）作為主數(shù)據(jù)源，RA-DARSAT-2寬波束模式，分辨率為距離向30m，方位向26m，圖像為150^x150km，面積為22500km²，極化方式有HH、VV、HV和VH，本次選用VV，入射角為20°～45°，其基本信息見表2。

3.數(shù)據(jù)處理

對于采空塌陷快速形變類型，本次以GAMMA軟件為數(shù)據(jù)平臺，采用二軌法對獲取的數(shù)據(jù)進(jìn)行差分干涉處理，獲取了采空塌陷的范圍和沉降速率，主要步驟包括：（1）選取兩期影像，利用主輔圖像生成干涉圖，去除平地效應(yīng)，干涉條紋圖中包含地形信息和形變信息；（2）將DEM數(shù)據(jù)轉(zhuǎn)換為相位值，并獲得模擬干涉條紋圖；（3）從主輔圖像生成的干涉條紋圖中減去利用DEM模擬的干涉條紋圖；（4）對差分干涉相位圖進(jìn)行處理，獲得地表形變量；（5）將獲取的地表形變圖轉(zhuǎn)換到地理坐標(biāo)系中。通過以上處理，獲得了干涉強(qiáng)度圖、干涉圖等圖件（圖3、圖4），可以看出，獲得效果較好。

4.采空塌陷區(qū)分析

通過常規(guī)差分干涉處理，獲得了研究區(qū)的差分干涉圖（圖5），從圖上可以準(zhǔn)確的識別已發(fā)生的塌陷和潛在的塌陷區(qū)域，且通過干涉紋理測算出塌陷的沉降量（表3）。

從圖5可以看出：淮南煤礦區(qū)采空塌陷區(qū)域主要集中分布于大通區(qū)的大通——九龍崗一帶、謝家集區(qū)和八公山區(qū)的大山鎮(zhèn)——八公山——謝家集——李郢孜一帶、潘集區(qū)的董圩-泥河一帶、鳳臺縣的觀音——岳張集——新集——劉集一帶。根據(jù)塌陷區(qū)的分布，結(jié)合收集的礦山資料，采空塌陷區(qū)主要是由大型煤礦劉莊礦業(yè)區(qū)域、謝橋煤礦區(qū)域、張集煤礦區(qū)域、顧北煤礦區(qū)域、顧橋煤礦區(qū)域、丁集煤礦區(qū)域、潘集第三煤礦區(qū)域、潘集第一煤礦區(qū)域、新集二礦區(qū)域、新集三礦區(qū)域、李咀孜煤礦區(qū)域等礦山開采所形成。

5.野外實(shí)地驗(yàn)證

通過野外調(diào)查，結(jié)合2012年度的光學(xué)遙感影像和工作區(qū)礦山開采現(xiàn)狀分布圖等資料，進(jìn)行綜合分析，對采用常規(guī)差分干涉方法提取的采空塌陷區(qū)與實(shí)際情況做對比分析，其結(jié)果與實(shí)際情況吻合，結(jié)果真實(shí)可信。

如顧橋煤礦區(qū)域，通過遙感影像可以看出，在礦區(qū)內(nèi)地表已經(jīng)出現(xiàn)小范圍塌陷，積水（圖6），經(jīng)過雷達(dá)差分干涉處理，發(fā)現(xiàn)此位置在2012年2月4日與2012年2月28日期間出現(xiàn)沉降，沉降量達(dá)到4.5cm/24天，通過后續(xù)3、4月的差分干涉圖對比，發(fā)現(xiàn)采空區(qū)塌陷范圍在逐步加大，與實(shí)地調(diào)查的結(jié)果一致。

通過差分干涉處理，獲取圖7所示煤礦區(qū)域沉降量為4.2cm/24天，結(jié)合光學(xué)遙感數(shù)據(jù)和雷達(dá)數(shù)據(jù)提取的采空塌陷區(qū)域（圖8），進(jìn)行野外實(shí)地調(diào)查，發(fā)現(xiàn)此區(qū)域地勢明顯低于周圍，地上電線桿已發(fā)生傾斜，已開始形成塌陷（圖9）。

5.結(jié)論與展望

（1）InSAR技術(shù)可以有效地識別出光學(xué)遙感未能識別的或潛在的采空塌陷區(qū)；可以估算出沉降區(qū)域沉降速率，通過實(shí)地調(diào)查和已有資料的對比分析，結(jié)果與實(shí)際情況一致，真實(shí)可靠。

（2）通過大量數(shù)據(jù)積累，采用InSAR技術(shù)，可以做到長時序的區(qū)域地面塌陷進(jìn)行大范圍的監(jiān)測，為工礦企業(yè)和政府管理部門及時掌握塌陷分布和數(shù)量，也為地質(zhì)災(zāi)害預(yù)警提供新的監(jiān)測技術(shù)方法。