徐良驥，劉 哲，龐 會，王振兵

(安徽理工大學測繪學院，安徽 淮南 232001)

多源數(shù)據(jù)下老采空區(qū)上方地表殘余變形規(guī)律分析

徐良驥，劉 哲，龐 會，王振兵

(安徽理工大學測繪學院，安徽 淮南 232001)

為動態(tài)分析老采空區(qū)上方地表殘余變形規(guī)律，本文以淮南礦區(qū)新莊孜礦新淮工廠為例，綜合運用D-InSAR遙感解譯反演與地表移動變形監(jiān)測兩種方法分析其殘余變形規(guī)律，并以遙感解譯數(shù)據(jù)與實測數(shù)據(jù)為基礎資料，采用GM(1,1)模型進行地表殘余變形規(guī)律的動態(tài)預測。D-InSAR存檔數(shù)據(jù)反演結果顯示2009年1月—2011年12月期間地表最大下沉量為100 mm，最大下沉速率為2.5 mm/d；地表移動變形監(jiān)測結果顯示2012年1月—2013年12月期間地表最大下沉量為55 mm，最大下沉速率為0.45 mm/d；GM(1,1)模型預測結果顯示地表監(jiān)測點下沉速度小于1.67 mm/d，地表處于殘余衰退階段；同時模型相關系數(shù)R2>0.95，預測結果和曲線擬合精度可靠。研究成果可為煤礦老采空區(qū)上方建筑地基殘余變形監(jiān)測工作提供參考。

老采空區(qū)；地表殘余變形；D-InSAR；移動變形監(jiān)測；GM(1,1)

我國大部分以煤炭經(jīng)濟為主的城市都是“因煤而建、緣煤而興”，隨著礦業(yè)經(jīng)濟的發(fā)展和城市的擴張，建筑用地日益缺乏，很多重要建筑物不得不建在老采空上方。目前國內外學者主要采用殘余變形監(jiān)測、相似材料模擬試驗、數(shù)值模擬、開采沉陷預計反演等方法綜合分析老采空區(qū)上方地表殘余變形規(guī)律[1-5]，而對于采用高分辨率遙感存檔數(shù)據(jù)反演和現(xiàn)時實測沉降數(shù)據(jù)進行綜合評價研究的較少。本文以淮南礦區(qū)新莊孜礦新淮工業(yè)廣場為例，采用D-InSAR遙感解譯與地表殘余變形監(jiān)測結合的方法，對老采空區(qū)上方地表殘余變形規(guī)律進行動態(tài)分析。

1 研究區(qū)概況

新莊孜礦新淮工業(yè)廣場于2006年建于50年前開采的C13、B11b煤層老采空區(qū)上方的地表，建成投入使用后不久，地表出現(xiàn)了不均勻沉降，部分重要建(構)筑物先后出現(xiàn)了程度不等的損壞變形情況。新淮工廣井上下對照圖及建(構)筑物分布示意圖如圖1、圖2所示。

2 基于D-InSAR的殘余變形規(guī)律反演分析

2.1 數(shù)據(jù)源與解譯方法

采用ALOS PALSAR雷達影像數(shù)據(jù)為新淮工業(yè)廣場的影像數(shù)據(jù)源，同時將4景影像數(shù)據(jù)源組合成3組干涉像對，具體參數(shù)見表1。

采用D-InSAR雙軌差分干涉處理方法對表1中的干涉像對進行處理，經(jīng)配準、重采樣、濾波、平地效應消除、干涉條紋濾波和相位解纏、相高轉換等數(shù)據(jù)處理后[6-8]，獲得地表沉降分布圖。

圖1 B11b、C13槽井上下對照示意圖

2.2 D-InSAR反演結果分析

將研究區(qū)的高分辨率影像圖、D-InSAR解譯圖

像、下沉等值線圖進行疊加分析,如圖3所示，反演結果顯示新淮工業(yè)廣場地表殘余變形規(guī)律為：2009年1月—2010年1月期間最大下沉量為100 mm，最大下沉速率為2.5 mm/d；2010年1月—2011年2月期間，最大下沉量為27 mm，最大下沉速率為0.07 mm/d。

圖2 新淮工業(yè)廣場內建(構)筑物分布示意圖

干涉像對主圖像時間輔圖像時間時間基線/d空間基線/m臨界基線/m12009-01-272009-12-1532214801306822009-12-152010-01-30466731307432010-01-302011-02-02368165913070

圖3 新淮工業(yè)廣場地表殘余形變分析

3 基于地表移動觀測站的殘余變形規(guī)律監(jiān)測分析

3.1 地表移動觀測站設置

為分析新淮工業(yè)廣場殘余變形動態(tài)演變規(guī)律，在新淮工業(yè)廣場內設立地表移動觀測站并進行實時監(jiān)測，如圖4所示。觀測站由T線和N線組成，其中T線長510 m，觀測點16個；N線長716 m，觀測點10個，2012—2013年間共進行了7次日常觀測和4次全面觀測。

圖4 新淮工業(yè)廣場地表移動觀測站布設示意圖

3.2 殘余變形監(jiān)測結果分析

圖5(a)、(b)分別為觀測站T、N線的下沉曲線圖。觀測期間，T線中T1監(jiān)測點的下沉量最大，為55.6 mm，T11監(jiān)測點的下沉量最小，為21.6 mm；N線中N8監(jiān)測點的下沉量最大，為52.5 mm；N15監(jiān)測點下沉量最小，為23.5 mm。從不同時期下沉曲線可以發(fā)現(xiàn),T、N兩條監(jiān)測線的監(jiān)測點在初始時期下沉量逐漸變大,達到峰值后趨于穩(wěn)定的規(guī)律；T、N兩條監(jiān)測線上采空區(qū)邊緣點的下沉量逐漸增大且下沉速度非常緩慢。

圖5 新淮工業(yè)廣場殘余變形實測曲線

圖5(c)、(d)分別為T、N線的下沉速率曲線圖，T、N線的下沉速率均小于1.67 mm/d，處于下沉非活躍期[9]，對建筑地基不會產生危害。

3.3 殘余變形規(guī)律預測

根據(jù)地表移動觀測站監(jiān)測點設置情況，在采空區(qū)邊緣分別選取了N8、N13、N16、T1、T10、T15 6個特征點，采用GM(1,1)模型結合D-InSAR歷史數(shù)據(jù)和實時監(jiān)測數(shù)據(jù)對地表殘余變形動態(tài)發(fā)展規(guī)律進行預測[10-12]，預測結果如圖6所示，預測值與實測值的殘差最大值為2.4 mm，相關系數(shù)R2>0.95表明曲線擬合精度很高。

圖6 各監(jiān)測點殘余變形反演值、實測值與預測值對比分析

由圖6可知，選擇的監(jiān)測點均表現(xiàn)出基本一致的殘余變形規(guī)律，監(jiān)測點下沉量達到最大值后逐漸趨于穩(wěn)定。試驗結果表明，將2009—2011年間的雷達影像存檔數(shù)據(jù)反演結果和2012—2013年間觀測站實測數(shù)據(jù)相結合，再采用GM(1,1)模型預測提高了預測結果的精度和時效性。

4 結 論

(1) 采用2009—2011年間的ALOS衛(wèi)星PALSAR數(shù)據(jù)，通過雙軌差分干涉處理方法獲得新淮工業(yè)廣場的沉降分布圖和下沉等值線圖，反演結果顯示研究區(qū)殘余變形最大下沉量為100 mm，最大下沉速率為2.5 mm/d。

(2) 地表移動觀測站監(jiān)測結果顯示，監(jiān)測點和老采空區(qū)的相對位置影響其殘余下沉量，保護煤柱上方的地表監(jiān)測點的殘余下沉量較小，靠近老采空區(qū)的地表監(jiān)測點殘余下沉量較大。

(3) 采用GM(1,1)模型結合D-InSAR反演結果與觀測站實測數(shù)據(jù)進行地表殘余變形動態(tài)預計能夠提高預計結果的精度與時效性，實現(xiàn)了老采空區(qū)殘余變形的歷史反演和未來發(fā)展預測的有機結合。

[1] 王磊,郭廣禮，查劍鋒，等.老采空區(qū)地表殘余沉降預計與應用[J].采礦與安全工程學報,2011,28(2):283-287.

[2] 李斌，鄧喀中，牛海鵬.老采空區(qū)上方地基穩(wěn)定性評價[J].煤炭工程，2012，7(7)：64-67.

[3] 徐良驥,郭輝,朱楠,等.建筑荷載作用下老采空區(qū)上方覆巖與地表殘余變形規(guī)律研究[J].測繪通報,2015(5):103-108.

[4] 方磊,李春苗,裴宗平.煤礦采空區(qū)穩(wěn)定性評價[J].地下空間與工程學報,2013, 9(S2):2034-2039.

[5] 王亞軍,賀江洪,閆鵬陽.老采空區(qū)地基穩(wěn)定性附加應力-概率積分法評價[J].金屬礦山,2013,42(2):135-139.

[6] 劉曉菲,鄧喀中,范洪東,等.D-InSAR監(jiān)測老采空區(qū)殘余變形的實驗[J].煤炭學報,2014,39(3):467-472.

[7] 胡波,汪汗勝.二軌法DInSAR技術監(jiān)測城市地表沉降[J].測繪工程,2010,19(2):42-46.

[8] DING X L,LI Z W,ZHU J J,et al.Atmospheric Effects on InSAR Measurements and Their Mitigation[J].Sensors,2008(8):5426-5448.

[9] 何國清，楊倫，凌賡娣,等.礦山開采沉陷學[M].徐州：中國礦業(yè)大學出版社，1991.

[10] 王正帥, 鄧喀中.老采空區(qū)殘余沉降的離散灰色預測模型[J].煤炭學報,2010,35(7):1084-1088.

[11] 李斌,朱健.非等間隔灰色GM(1,1)模型在沉降數(shù)據(jù)分析中的應用[J].測繪科學,2007,32(4):52-55.

[12] 佘娣,謝劭峰,彭家頔,等. 非等時距GM(1,1)在變形預報中的應用[J]. 測繪科學,2013,38(3):110-111.

Analysis of Old Goaf Surface Residual Deformation Regulation Based on Multi-source Data

XU Liangji,LIU Zhe,PANG Hui,WANG Zhenbing

(Faculty of Surveying and Mapping,Anhui University of Science and Technology, Huainan 232001, China)

For analyzing the regulation of surface residual deformation above old goaf, the surface residual deformation regulation above Xinhuai factory by remote sensing interpretation method and surface residual deformation monitoring is studied, and the dynamic residual deformation regulation with historical archive data and real-time data monitoring based on GM(1,1) model is predicted. According to the inversion of historical archive data. The result shows that the maximum sinkage is 100 mm from January 2009 to December 2011; the maximum residual sinkage is 55 mm from January 2012 to December 2013; the residual subsidence of surface monitoring point is stable. The research result could provide reference for the building foundation residual deformation monitoring work above old goaf.

old goaf; surface residual deformation; D-InSAR; monitoring of surface movement and deformation; GM (1,1)

徐良驥，劉哲，龐會，等.多源數(shù)據(jù)下老采空區(qū)上方地表殘余變形規(guī)律分析[J].測繪通報，2017(2)：45-48.

10.13474/j.cnki.11-2246.2017.0046.

2016-10-10

國家自然科學基金(41472323)；安徽省對外科技合作計劃項目(1503062020)

徐良驥(1978—)，男，教授，研究方向為礦山開采形變監(jiān)測。E-mail:ljxu@aust.edu.cn

P237

A

0494-0911(2017)02-0045-04