姚丹丹 吳 侃 何 強

(中國礦業(yè)大學環(huán)境與測繪學院，江蘇 徐州 221116)

基于D-InSAR技術(shù)的煤礦區(qū)沉陷監(jiān)測

姚丹丹 吳 侃 何 強

(中國礦業(yè)大學環(huán)境與測繪學院，江蘇 徐州 221116)

煤炭開采引起的地表移動變形時間集中、形變量大，使得基于D-InSAR技術(shù)的煤礦區(qū)沉陷監(jiān)測受到了諸多的限制，很難獲取準確的地表沉陷信息。為使D-InSAR技術(shù)能夠更好地應(yīng)用于煤礦區(qū)沉陷監(jiān)測，采用“兩軌法”差分干涉技術(shù)處理了多景RadarSat-2影像數(shù)據(jù)，通過實驗分析了植被因素對影像相干性的影響，探究了監(jiān)測量級對地表沉陷信息解譯結(jié)果的影響，并以峰峰礦區(qū)的九龍礦為試驗區(qū)對D-InSAR技術(shù)的監(jiān)測精度進行了分析。結(jié)果表明：地表植被覆蓋率是影響影像間相干性的決定性因素，在地表植被稀疏的冬季相比植被茂盛的夏季更適合使用D-InSAR技術(shù)進行地表形變監(jiān)測；由于地質(zhì)采礦條件復雜，導致地表下沉速度超過衛(wèi)星可監(jiān)測范圍，通常只能準確監(jiān)測到盆地邊緣的形變；當沉降量不超過監(jiān)測量級時，D-InSAR技術(shù)的監(jiān)測結(jié)果與常規(guī)水準測量的結(jié)果大致吻合；隨著相關(guān)理論研究的不斷深入與處理軟件的日趨成熟，D-InSAR技術(shù)在煤礦區(qū)沉陷監(jiān)測中的應(yīng)用前景將更為廣闊。

D-InSAR 礦區(qū)沉陷 形變監(jiān)測 植被 監(jiān)測量級 監(jiān)測精度

長期以來，煤炭資源的開發(fā)與利用為我國經(jīng)濟發(fā)展和社會進步做出了巨大貢獻，但也不可避免地嚴重破壞了礦區(qū)生態(tài)環(huán)境，使得礦區(qū)及其附近產(chǎn)生了不同程度的地面沉降。與傳統(tǒng)的工作量大、觀測時間長、成本高、空間分辨率低相比[1]，D-InSAR技術(shù)具有快捷方便、低成本、高精度、對微小變形的高敏感度等優(yōu)點，已被廣泛應(yīng)用于地面沉降方面的監(jiān)測。但由于煤礦區(qū)地表條件復雜且變化迅速，地表沉降量大、形變時間集中，使得D-InSAR技術(shù)在監(jiān)測煤炭開采引起的礦區(qū)地表沉陷應(yīng)用上受到了諸多的限制，亟需對基于D-InSAR技術(shù)的煤礦區(qū)地面沉降監(jiān)測技術(shù)體系、算法進行研究改進，由此發(fā)展出了PS技術(shù)、CPT技術(shù)、相干目標分析、小基線集算法等新技術(shù)[1-7]。為此，本研究利用GAMMA軟件，分析了植被和監(jiān)測量級對D-InSAR礦區(qū)沉降監(jiān)測技術(shù)的影響，通過實驗探究了植被因素對影像相干性的影響，根據(jù)實例驗證了在監(jiān)測量級范圍內(nèi)將D-InSAR技術(shù)應(yīng)用于煤礦區(qū)地面沉降監(jiān)測可行性，并在此基礎(chǔ)上分析了沉降監(jiān)測的精度，為改進基于D-InSAR的煤礦區(qū)地面沉降監(jiān)測技術(shù)提供解決思路。

1 D-InSAR變形監(jiān)測限制因素

D-InSAR技術(shù)在監(jiān)測煤炭開采引起的礦區(qū)地表沉陷應(yīng)用上受到了諸多的限制：①影像失相干使得D-InSAR技術(shù)無法正常反演出地表特征，影響去相干的因素有空間、時間、熱噪聲及數(shù)據(jù)處理過程的失相干等；②大氣效應(yīng)是干涉處理過程中非常重要的影響因素，大氣延遲相位可能掩埋掉地形相位或形變相位，導致2幅影像的失相干；③當?shù)乇硇》秶霈F(xiàn)較大的移動變形會造成干涉圖相位的急劇變化而使干涉條紋變得模糊，使得地表變形造成的沉降量很難應(yīng)用D-InSAR技術(shù)監(jiān)測到[8-11]。

2 地表植被因素的影響

影像對的相干性直接影響了差分干涉測量結(jié)果的質(zhì)量，關(guān)系到能否獲得可用的地表形變。導致時間失相干的因素主要有：①季節(jié)交替或農(nóng)作物收割導致的地表植被覆蓋變化，②季節(jié)交替導致的冰雪覆蓋變化，③天氣原因?qū)е碌牡乇頋穸茸兓芡寥莱煞值淖兓取?/p>

為了研究植被變化對D-InSAR技術(shù)的影響，實驗選取了同一區(qū)域不同時間段的多景Radarsat-2影像，影像位于平原區(qū)域，建筑物、山區(qū)、河流變化因素可以忽略不計，實驗數(shù)據(jù)如表1所示。不同時序的相干系數(shù)圖見圖1所示，統(tǒng)計直方圖如圖2所示。

表1 實驗設(shè)計的影像對

由圖2可知，實驗二影像獲取時間處于冬季，地表植被覆蓋率很低，相干性最高，平均相干系數(shù)值為0.932 3，獲得的干涉相位圖條紋也最為清晰；實驗一和實驗三影像獲取時間分別處于秋、春2季，植被覆蓋率較高，相干性較好，實驗一、實驗三的平均相干系數(shù)值分別為0.792 3、0.800 8，也可形成干涉圖，但條紋較密；實驗四影像獲取時間處于夏季，植被生長茂盛、覆蓋率高，導致相干性大大降低，平均相干系數(shù)值為0.320 4，無法獲得干涉相位圖。

圖1 相干系數(shù)

圖2 相干系數(shù)統(tǒng)計直方圖

影像監(jiān)測受到植被因素影響，產(chǎn)生季節(jié)性失相干，冬季相干性最佳，春、秋2季相干性較好，均可獲取沉降圖，但夏季相干性很低，該季節(jié)的數(shù)據(jù)無法使用。但是，礦區(qū)沉陷變形持續(xù)時間一般為1 a甚至數(shù)年，時間跨度大、季節(jié)交替變化，受植被因素影響，處于夏季時間段內(nèi)的監(jiān)測影像數(shù)據(jù)相干性很低甚至失相干，導致該季節(jié)期間的沉陷情況均無法監(jiān)測到，從而很難獲取礦區(qū)長時間緩慢累積的地表形變。

3 沉陷監(jiān)測量級的影響

在單位像元內(nèi)能夠監(jiān)測到的臨界形變梯度是一個條紋周期，當?shù)乇淼男巫兞考壋^這個臨界形變梯度時就會造成影像之間的空間失相干，從而不能夠解算出真實的形變量[2]。D-InSAR技術(shù)可監(jiān)測到的臨界形變量可表示為[12]：

(1)

式中，Wmax為地表可監(jiān)測的最大沉降量；r為沉降盆地的主要影響半徑；g為雷達影像的地面分辨率；λ為傳感器波長。

根據(jù)各衛(wèi)星系統(tǒng)所采用的波長和地面分辨率的不同，由式(1)可計算得到其能監(jiān)測的地表最大下沉值。當假設(shè)沉降盆地的主要影響半徑r為100 m時，目前常用衛(wèi)星系統(tǒng)可監(jiān)測的最大地表沉降量如表2。

以RadarSat-2衛(wèi)星為例，由表2可知，該衛(wèi)星可監(jiān)測最大沉降量為28 cm，重訪周期為24 d，平均每天可監(jiān)測到的下沉量為11.7 mm。當下沉速度超過11.7 mm/d時，超出該衛(wèi)星可監(jiān)測量級，導致大的形變量監(jiān)測不到。

表2 不同衛(wèi)星系統(tǒng)可探測的最大地表沉降量

注：不同衛(wèi)星系統(tǒng)具有多種不同的產(chǎn)品模式，其地表分辨率不同，此表中僅列出了各衛(wèi)星系統(tǒng)常見分辨率的產(chǎn)品。

[13]，地表最大下沉速度V為

(2)

式中，K為下沉速度系數(shù)，一般取為1.8；C為工作面推進速度，m/d；Wfm為地表最大下沉值，m；H0為平均開采深度，m；M為煤層厚度，m；q為下沉系數(shù)；α為煤層傾角,(°)。

假設(shè)某一工作面的推進速度為6m/d，下沉速度系數(shù)取1.8、下沉系數(shù)取0.8、煤層傾角取10°，只有最大下沉速度小于11.7mm/d，才可以監(jiān)測到地表沉降量。若開采深度為300m，根據(jù)式(2)計算可知，煤層厚度要小于等于0.4m。反之，若煤層厚度已知為1m，則開采深度要大于等于727m，否則無法監(jiān)測到正確的地面沉降。

由于礦區(qū)地質(zhì)條件復雜，采煤工作面的布設(shè)常受到煤層埋藏深度、煤層厚度等因素的影響，開采深度和煤層厚度無法同時滿足臨界值的要求，導致開采產(chǎn)生的最大沉降超過監(jiān)測量級，通常只能獲取盆地邊緣小形變量，而不能監(jiān)測到完整的沉陷盆地，限制了D-InSAR技術(shù)在煤礦區(qū)沉陷監(jiān)測中的應(yīng)用范圍。

4 D-InSAR技術(shù)沉陷監(jiān)測精度分析

目前，普遍采用常規(guī)觀測站的水準測量結(jié)果對D-InSAR技術(shù)的監(jiān)測結(jié)果進行評價。峰峰礦區(qū)九龍礦在152下31S工作面上方布設(shè)有地表移動觀測站，選擇覆蓋九龍礦的RadarSat-2影像數(shù)據(jù)，將2011-01-09和2011-02-02的影像形成干涉對，影像信息見表3，此期間正在開采的工作面如圖3所示粗線條框格部分，與影像監(jiān)測到的塌陷區(qū)進行對比，結(jié)果如圖4所示。

表3 SAR影像信息

圖3 九龍礦152下31S工作面示意

圖4 九龍礦監(jiān)測塌陷區(qū)與工作面對應(yīng)圖

將應(yīng)用D-InSAR方法監(jiān)測的地表沉降量與觀測站監(jiān)測的沉降量進行對比如圖5所示，從圖中可知，兩者的最小絕對誤差為0 mm，最大絕對誤差為30 mm，18個點誤差在10 mm以內(nèi)，11個點誤差在10 mm及以上，平均誤差為8.97 mm，標準偏差為9.68 mm。在監(jiān)測時間段內(nèi)，影像監(jiān)測到礦區(qū)沉降在100 mm以內(nèi)，不超過監(jiān)測量級280 mm，與地表觀測站實測數(shù)據(jù)進行對比，可知兩者誤差較小，監(jiān)測結(jié)果準確且接近實際情況。

圖5 D-InSAR與水準測量監(jiān)測沉降量對比

另外，利用淮南礦區(qū)的ALOS PALSAR數(shù)據(jù)，提取出2007年12月10日到2008年1月25日這段時間內(nèi)礦區(qū)的地面沉陷，其中潘一礦、潘二礦、潘三礦、顧橋礦北部、顧北礦東部、張集礦、新集礦、花園湖礦、謝橋礦的地面沉陷明顯，利用D-InSAR技術(shù)監(jiān)測獲取的沉降圖，可以清晰地標示沉陷區(qū)域和下沉程度，監(jiān)測到的最大地面沉降量值為508.6 mm，在該衛(wèi)星最大監(jiān)測量級(580 mm)范圍內(nèi)，滿足沉陷監(jiān)測量級的要求。

D-InSAR技術(shù)可應(yīng)用于礦區(qū)沉陷監(jiān)測，與實測結(jié)果大致吻合，并具有一定的可靠性和精度。但是，當?shù)V區(qū)沉降量超過衛(wèi)星監(jiān)測量級時，則無法獲得準確的沉陷信息。

對于礦區(qū)地表沉陷監(jiān)測而言，一方面應(yīng)該盡可能地設(shè)置觀測站、獲取地面實測數(shù)據(jù)，并與D-InSAR技術(shù)監(jiān)測結(jié)果進行對比分析，在一定程度上定量分析D-InSAR監(jiān)測精度的可靠性；另一方面，針對D-InSAR技術(shù)監(jiān)測中的誤差源和干涉圖像質(zhì)量，建立合理的定性或定量分析的可靠性評價體系，是D-InSAR技術(shù)精度評定的技術(shù)關(guān)鍵[8]。

5 結(jié) 論

D-InSAR技術(shù)具有常規(guī)監(jiān)測方法不可比擬的優(yōu)勢，在礦區(qū)沉陷監(jiān)測方面的應(yīng)用前景非常廣闊。

(1)地表植被覆蓋率是影響影像間相干性的決定性因素，植被稀疏區(qū)域即可應(yīng)用D-InSAR技術(shù)進行地表形變監(jiān)測，不受影像成像時間和地形條件的限制，但是夏季植被覆蓋率高、降低了影像相干性，無法成功獲取該時間段內(nèi)的沉降量。

(2)不同衛(wèi)星系統(tǒng)可探測最大地表沉降量不同，以RadarSat-2衛(wèi)星監(jiān)測量級為例，只有當工作面開采條件達到監(jiān)測量級允許范圍內(nèi)，才能夠正確獲取地面沉降。

(3)通過對九龍礦以及淮南礦沉陷監(jiān)測的實例說明，在影像監(jiān)測量級范圍內(nèi)，D-InSAR技術(shù)可以很好地應(yīng)用于礦區(qū)監(jiān)測，并具有較高的精度。

(4)隨著高分辨率合成孔徑雷達的發(fā)射、干涉處理中關(guān)鍵算法與處理軟件的日趨完善成熟、D-InSAR相關(guān)技術(shù)如永久散射體、短基線的發(fā)展以及與其他數(shù)據(jù)源融合技術(shù)的研究，使得差分干涉測量技術(shù)在礦區(qū)的應(yīng)用前景更為廣闊。

參 考 文 獻

[1] 吳 侃,汪云甲,王歲權(quán),等.礦山開采沉陷監(jiān)測及預測新技術(shù)[M].北京:中國環(huán)境科學出版社,2012:83-151. Wu Kan，Wang Yunjia，Wang Suiquan，et al.New Technologies of Monitoring and Forecasting Mining Subsidence[M].Beijing:China Environmental Science Press,2012:83-151.

[2] 刁鑫鵬.D-InSAR技術(shù)在峰峰礦區(qū)開采沉陷監(jiān)測中的應(yīng)用研究[D].徐州：中國礦業(yè)大學,2013. Diao Xinpeng.The Study of D-InSAR Technology on Mining Subsidence Monitoring in Fengfeng[D].Xuzhou:China University of Mining and Technology,2013.

[3] 董玉森,葛玲玲,張新春，等.基于差分雷達干涉測量的礦區(qū)地面沉降監(jiān)測研究[J].武漢大學學報:信息科學版,2007,32(10):888-892. Dong Yusen，Ge Linlin，Zhang Xinchun，et al.Mine subsidence monitoring by differential InSAR[J].Geomatics and Information Science of Wuhan University:Information Science,2007,32(10):888-892.

[4] 傅文學,田慶久,郭小方,等.PS技術(shù)及其在地表形變監(jiān)測中的應(yīng)用現(xiàn)狀與發(fā)展[J].地球科學進展,2006,21(11):1193-1198. Fu Wenxue，Tian Qingjiu，Guo Xiaofang，et al.The permanent scatterers technique and its application to monitoring ground deformation[J].Advances in Earth Science,2006,21(11):1193-1198.

[5] 葛大慶,王 艷,郭小方,等.利用短基線差分干涉紋圖集監(jiān)測地表形變場[J].大地測量與地球動力學,2008,28(2):61-66. Ge Daqing，Wang Yan，Guo Xiaofang，et al.Surface deformation field monitoring by use of small-baseline differential interferograms stack[J].Journal of Geodesy and Geodynamics,2008,28(2):61-66.

[6] 尹宏杰,朱建軍,李志偉,等.基于SBAS的礦區(qū)形變監(jiān)測研究[J].測繪學報,2011,40(l):52-58. Yin Hongjie，Zhu Jianjun，Li Zhiwei，et al.Ground Subsidence monitoring in mining area using SBAS algorithm[J].Acta Geodaetica et Cartographica Sinica,2011,40(l):52-58.

[7] 湯益先.基于相干目標的 DInSAR 方法及其地表沉降應(yīng)用研究[D].北京：中國科學院研究生院,2006. Tang Yixian.Research of Coherent Target Interferometry and Its Applications to Surface Subsidence Fields[D].Beijing：Graduate University of Chinese Academy Sciences.,2006.

[8] 吳立新,高均海,葛大慶,等.工礦區(qū)地表沉陷D-InSAR監(jiān)測實驗研究[J].東北大學學報:自然科學版,2005(8):778-781. Wu Lixin，Gao Junhai，Ge Daqing，et al.Experimental study on surface subsidence monitoring with D-InSAR in mining area[J].Journal of Northeastern University：Natural Science,2005(8):778-781.

[9] 成 樞,劉 元,張洪勝.基于D-InSAR的開采沉陷監(jiān)測技術(shù)[J].測繪與空間地理信息,2013,36(6):7-9. Cheng Shu，Liu Yuan，Zhang Hongsheng.Exploration and settlement monitoring technology based on D-InSAR[J].Geomatics & Spatial Information Technology,2013,36(6):7-9.

[10] 田 馨,廖明生.InSAR技術(shù)在監(jiān)測形變中的干涉條件分析[J].地球物理學報,2013,56(3):812-823. Tian Xin，Liao Mingsheng.Analysis of conditions for InSAR in the field of deformation monitoring[J].Chinese Journal of Geophysics，2013,56(3):812-823.

[11] 馬海濤,李 輝,劉勇峰,等D-InSAR技術(shù)在礦區(qū)地表沉降監(jiān)測中的應(yīng)用[J].金屬礦山,2011(2):95-98. Ma Haitao，Li Hui，Liu Yongfeng，et al.Application of D-InSAR technique to the land subsidence monitoring in mining area[J].Metal Mine,2011(2):95-98.

[12] 武繼峰,楊志強.InSAR技術(shù)應(yīng)用于礦區(qū)開采沉陷監(jiān)測的數(shù)據(jù)優(yōu)選研究[J].測繪通報,2013(增刊):92-94. Wu Jifeng，Yang Zhiqiang.Data optimization of mining subsidence monitoring by InSAR technology[J].Bulletin of Surveying and Mapping,2013(S):92-94.

[13] 何國清,楊 倫,凌賡娣,等.礦山開采沉陷學[M].徐州:中國礦業(yè)大學出版社,1991：89-130. He Guoqing，Yang Lun，Ling Gengdi，et al.Study of Mining Subsidence[M].Xuzhou:China University of Mining and Technology Press,1991：89-130.

[14] 鄧喀中,姚 寧,盧 正,等.D-InSAR監(jiān)測開采沉陷的實驗研究[J].金屬礦山,2009(12):25-27. Deng Kazhong，Yao Ning，Lu Zheng，et al.Experimental research on monitoring mining subsidence by D-InSAR technique[J].Metal Mine,2009(12):25-27.

[15] 王亞男.InSAR技術(shù)用于礦區(qū)大量級塌陷監(jiān)測研究[D].西安：長安大學,2011. Wang Yanan.Research on Large-scale Mining Collapse Monitoring with InSAR Technology[D].Xi'an:Chang'an University,2011.

[16] 蔣 彌,李志偉,丁曉利,等.InSAR可檢測的最大最小變形梯度的函數(shù)模型研究[J].地球物理學報,2009,52(7):1715-1724. Jiang Mi，Li Zhiwei，Ding Xiaoli，et al.Study on the maximum and minimum detectable deformation gradients resolved by InSAR[J].Chinese Journal of Geophysics,2009,52(7):1715-1724.

[17] 閆建偉,汪云甲,陳國良,等.錢營孜煤礦地表沉陷的D-InSAR監(jiān)測[J].金屬礦山,2011(3):105-107. Yan Jianwei，Wang Yunjia，Chen Guoliang，et al.Ground subsidence monitoring with D-InSAR in Qianyingzi Coal Mine[J].Metal Mine,2011(3):105-107.

(責任編輯 徐志宏)

Coal Mine Subsidence Monitoring Based on D-InSAR Technology

Yao Dandan Wu Kan He Qiang

(School of Environmental Science and Spatial Informatics，China University of Mining and Technology，Xuzhou 221116，China)

Due to time concentration and large deformation in ground movement and deformation caused by coal mining，there are various limitations in monitoring coal mining subsidence with D-InSAR，and it is difficult to obtain accurate surface subsidence information.In order to apply D-InSAR technology to monitoring land subsidence in coal-mining areas，RadarSat-2 data are processed by “two-pass” differential interference technology.The effects of vegetation factor on the coherence between images were analyzed through experiments，and the influence of monitoring magnitude on the interpretation results of surface subsidence was explored.With Jiulong mine in Fengfeng coal mining area as an experimental zone，the monitoring precision based on D-InSAR technology was analyzed.The results show that vegetation coverage rate is the decisive factor of affecting the coherence between images，and D-InSAR technology can be better used in less vegetation in winter than more vegetation in summer.Subsidence velocity exceeds the range that can be monitored by the satellite and only the deformation of the basin edge can be accurately monitored due to the complexity of geological and mining conditions；when mining subsidence does not exceed monitoring magnitude，the monitoring results by D-InSAR technology are consistent with the usual leveling monitoring data.With the depth of relevant theory and the maturity of processing software，the D-InSAR technology will have a broader application prospect in coal-mining areas.

D-InSAR，Mining subsidence，Deformation monitoring，Vegetation，Monitoring magnitude，Monitoring precision

2014-09-01

江蘇高校優(yōu)勢學科建設(shè)工程項目(編號：SZBF2011-6-B35)。

姚丹丹(1989—)，女，碩士研究生。

TD325.2

A

1001-1250(2014)-11-151-05