寒區充填開采地表變形監測數據穩定性分析

郭慶彪 郭廣禮 王金濤

(1.中國礦業大學環境與測繪學院，江蘇 徐州 221116；2.國土環境與災害監測國家測繪地理信息中心重點實驗室，江蘇 徐州 221116；3.資源環境信息工程江蘇省重點實驗室，江蘇 徐州 221116)

寒區充填開采地表變形監測數據穩定性分析

郭慶彪1,2,3郭廣禮1,2,3王金濤1,2,3

(1.中國礦業大學環境與測繪學院，江蘇 徐州 221116；2.國土環境與災害監測國家測繪地理信息中心重點實驗室，江蘇 徐州 221116；3.資源環境信息工程江蘇省重點實驗室，江蘇 徐州 221116)

為了準確獲取充填面地表移動變形數據，有必要對各期監測數據的質量進行檢驗。利用F檢驗和t檢驗方法對寒區充填面觀測站中監測點進行穩定性分析，結果表明觀測期間控制點具有穩定性，進一步分析發現目標點存在3種形式的豎向移動：一是采煤引起的地表沉降，二是長期的汽車荷載作用下路基沉降現象，三是季節性凍土層的融沉現象。最后一種現象在我國高緯度地區較為普遍，該現象的產生將擴大地表移動盆地邊界，減小主要影響正切值，不利于該類地質采礦條件下地表移動變形預計參數的準確獲取。

充填開采 統計檢驗 穩定性分析 季節性凍土層

固體充填開采作為綠色開采理念的重要組成部分，已成功應用于我國10多個礦區[1]。結合現階段工程實踐經驗，固體充填開采充實率可達85%，有效地限制覆巖運動，大幅減小地表移動變形極值，很好地保護地面建(構)筑物[2]。但相比于傳統垮落法開采，固體充填開采覆巖及地表移動規律的研究尚不充分，急需建立不同地質采礦條件下地表移動觀測站，獲取充填開采地表移動變形預計參數，完善固體充填開采巖層及地表控制體系。

觀測站中的監測點分為控制點和目標點，控制點作為唯一起算數據決定著目標點的絕對位移，而目標點的絕對位移將直接反應充填開采的地表移動規律。因此，在實際觀測過程中，監測點的穩定性將間接地影響所獲取的充填開采地表移動規律可靠性。本研究以獲取內蒙古某礦充填開采地表移動規律為研究背景，建立了1801充填面地表移動變形觀測站，借助F檢驗和t檢驗對觀測站中監測點進行穩定性檢驗[3-5]。

1 穩定性檢驗數學模型

觀測站初期觀測中，點位出現微弱沉降現象，接近于測量誤差，無法判斷該變形是誤差所致還是真正的點位變化，重復觀測仍無法剔除上述點位，因此，利用F檢驗和t檢驗對點位穩定性進行檢驗。

1.1 F檢驗數學模型

(1)

式中，α為顯著水平(一般取0.05)，f為自由度。

(2)

(3)

式中，V為方差改正數；d=H2-H1(H為高程)；Pd、P為權重；h為間隙d協因數陣Q的秩。

(4)

如果上式成立，認為所有檢驗點位都是穩定的，反之認為檢測點位存在動點。

1.2 t檢驗數學模型

當式(4)不成立時，檢測點位存在動點，但不代表全部點位都是不穩定的，此時再利用t檢驗法進行單點位移的顯著性檢驗。假設2期觀測檢測點穩定，構造統計量

(5)

若上式成立，則接受原假設，該檢驗點穩定，反之該檢測點變化。

2 實例分析

2.1 地表移動觀測站布設

內蒙古某礦1801充填面設計面寬150 m,推進長度1 116 m，煤層傾角3°～ 12°，采高2 m，平均開采深度為450 m，有夾矸。8煤層位于二迭系山西組下部的第一巖段的中部，頂板為砂質泥巖，底板為砂泥巖，裂隙不發育，遇水不膨脹。地質構造屬于簡單型，直接、間接充水含水層的含水空間以裂隙為主，富水性中等。設計采用綜合機械化矸石充填采煤工藝，密實矸石充填法管理頂板。依據《煤礦測量規程》第259條的規定可取地表移動觀測點間距為25 m。布設控制點3個，走向觀測點26個。平面測量儀器可選用測角精度為2″、測距精度為2+2×10-6D(D為觀測距離， km)的全站儀，水準測量采用DS1級高精度水準儀，按三等水準測量精度要求進行所有高程測量工作，觀測站布設情況見圖1。并于2013年12月與2014年4月分別進行首次觀測與日常第一次觀測，利用上述統計檢驗方法分別對控制點和目標點的穩定性進行檢驗。

圖1 觀測站布設示意

2.2 控制點穩定性分析

地表移動觀測站的控制點不是嚴格意義的控制點，將控制點的高程作為待定點處理，利用秩虧自由網平差方法對控制點高程進行平差處理。取長度為1 m的水準測量路線觀測高差中誤差為單位權中誤差。首次觀測相關數據見表1，日常觀測第一次相關數據見表2。

表1 首次觀測相關數據

表2 日常第1次觀測相關數據

2.3 目標點穩定性分析

經過上述分析，2期觀測期間控制點保持穩定狀態，觀測站控制點高程可作為必要的起算數據，利用間接平差對目標點高程進行平差處理。定權方法與上相同。由于目標點較多，基于篇幅限制，2期原始高差觀測數據在此不再羅列。

得到首次觀測的參數改正數與單位權方差

0.0,0.5,1,1,0.5,1,1.5,1.5,2,2.5,2.0,2.0,

得到日常觀測第一次的各參數改正數與單位權方差

5.0,2.0,1.5,2.0,0.5,0.0,0.5,0.5,1.5,5.0,4.0,

對平差后的目標點高程進行F檢驗，F0=1.095F(0.05/26，52)=1.819，拒絕原假設，目標點中存在動點。利用t檢驗法尋求動點。求得間隙

取0.05的顯著水平，t(0.05/2,52)=2.007，得到各點的穩定情況見表3。

剔除上述所有動點后，再次進行F檢驗，第一步無顯著差異通過，取0.05顯著水平，F0=2.105

2.4 動點成因分析

結合2期目標點高程差曲線(如圖2)與目標點的穩定性分析結果，可知觀測線存在2個下沉區間，分別為Z4～Z11段和Z19～Z27段。以下將對上述2區間下沉現象進行分析。

表3 基于t檢驗的觀測點穩定情況

圖2 2期觀測點高程差值

Z4～Z11測點段位于1801充填面采空區上方，地表無水系相連，最大下沉點Z8位于采空區中央(見圖1)，下沉量為5 mm，以Z8為中心對稱，該區間的下沉值向兩側逐漸減小，在一定距離后收斂為零(Z2和Z12的2次高程相同)，與采煤引起的地表下沉曲線分布形態相似。符合開采沉陷中相關理論：地表移動變形初期變形值很小，且僅發生在采空區中心小范圍內[6]。因此可判斷該區間目標點的下沉為地下煤層開采所致，屬采煤引起地表沉陷現象。

Z19～Z27測點段距離當前工作面推進距離較遠，由開采沉陷學相關理論可知，1801充填面的開采影響不會波及此處，同時該區附近無任何采掘活動，可排除采煤因素。由圖2所示，可將該區間下沉曲線分為2部分：Z19～Z22和Z23～Z27子區間段，其中Z19～Z22下沉值基本相當，而Z23～Z27區間段下沉值基本符合線性遞增規律。分析圖1觀測點點位示意圖，發現Z19～Z22區間段橫跨棋石一級公路，且Z20與Z21號監測點布設在路面上，Z19與Z23號監測點位于公路路基上，考慮到該公路車流量大，多為運煤車輛，載重量大，在長時間汽車荷載的影響下，公路路基發生沉降現象，最終導致公路及路基上所有測點下沉，并且下沉值基本相當。而Z23～Z27區間段緊鄰河流，地層含水率較高，考慮到該礦區位于我國寒區，寒區地層中季節性凍土層具有凍脹與融沉特性，觀測站2期觀測分別為2013年12月和2014年4月，12月為該地區室外平均溫度為-25 ℃左右，低于水的凝固溫度，而4月該區室外平均溫度為18 ℃左右，高于水的凝固溫度，氣候的變暖引起凍土層發生融沉現象。已有研究表明，凍土層的含水率與地表融沉量呈正相關[7-9]，而本區Z23～Z27區間段含水率遞增，因此出現Z23～Z27區間段下沉值線性遞增，在Z27處達到最大的現象。

綜上所述，目標點在2期觀測期間存在3種形式的豎向移動，一是采煤引起的地表沉降，二是長期的汽車荷載作用下路基沉降現象，三是季節性凍土層的融沉現象。其中融沉現象導致地表下沉約8 mm，接近盆地邊界10 mm，將擴大地表移動盆地邊界，使得盆地邊界收斂減慢，縮小移動盆地邊界角，減小主要影響角正切，不利于該類地質采礦條件下地表移動變形預計參數的準確獲取。

4 結論與展望

(1)利用秩虧自由網平差方法對控制點高程進行平差處理，再用F檢驗與t檢驗對控制點和觀測點進行穩定性檢驗，驗證了控制點在觀測期間保持穩定狀態，可作為監測網的必要起算數據。

(2)利用間接平差方法對觀測點高程進行平差處理，并對點位穩定性進行檢驗，結果顯示觀測點大多為不穩定點，觀測點存在3種形式的豎向移動：一種是采煤引起的地表沉降；二是長期的汽車荷載作用下路基沉降現象；三是季節性凍土層的融沉現象。后一現象的產生將擴大地表移動盆地邊界，減小主要影響正切值tanβ，不利于該類地質采礦條件下地表移動變形預計參數的準確獲取。

(3)本研究僅闡述1801充填面存在季節性凍土層融沉現象，但并未提出基于融沉現象的寒區監測數據修正模型。后續工作將結合該區凍土層融沉系數分布規律，建立修正模型，為準確獲取寒區充填采煤地表移動規律、地表移動動態參數及角值參數提供科學依據。

[1] 繆協興,錢鳴高.中國煤炭資源綠色開采研究現狀與展望[J].采礦與安全工程學報,2009,26(1):1-14. Miao Xiexing，Qian Minggao.China coal resource green mining research status and prospect[J].Journal of mining and safety engineering，2009,26(1):1-14.

[2] 楊寶貴,彭楊皓,李 楊,等.充填開采地表移動變形規律數值模擬分析[J].金屬礦山,2014(12):169-174. Yang Baogui，Peng Yanghao，Li Yang，et al.Numerical simulation analysis of ground movement deformation law in backfill mining[J].Metal Mine,2014(12):169-174.

[3] 陳永奇.變形觀測數據處理[M].北京：測繪出版社，1988. Chen Yongqi.Deformation Observation Data Processing[M].Beijing:Surveying and Mapping Press，1988.

[4] 王 巖,岳建平，周保興,等.工程控制網點位穩定性分析方法的研究[J].測繪通報,2004(8):12-14. Wang Yan，Yue Jianping，Zhou Baoxing，et al.A stability analysis research of engineering control network[J].Bulletin of Surveying and Mapping,2004(8):12-14.

[5] 邱國會,李子申,岳 昊,等.間隙分塊法和t檢驗法對高程控制網點位穩定性的分析對比[J].工程地球物理學報,2009,6(1):247-250. Qiu Guohui，Li Zishen，Yue Hao，et al.A comparison between gap block method andttest analysis of the stability of the elevation control network[J].Journal of Engineering Geophysics,2009，6(1):247-250.

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

[7] 王 芳.博牙高速公路島狀多年凍土地區路基沉降儲值技術研究[D].西安：長安大學,2011. Wang Fang.Research on Treatment Technology for Subgrade Settlement in Segregated Permafrost Region of Boketu-Yakesshi Expressway[D].Xi′an:Chang ′an University,2011.

[8] 侯曙光.多年凍土路基融沉機理及路面結構對策研究[D].南京：東南大學,2005. Hou Shuguang.Study on Mechanism of Sub-grade Thawing and Counter-measures Against Pavement Structure Distress in Permafrost Region[D].Nanjing:Southeast University,2005.

[9] 金 龍,汪雙杰,陳建兵,等.基于變形分析的多年凍土地區路基高度效應研究[J].中外公路,2013,33(3):22-28. Jin Long，Wang Shuangjie，Chen Jianbing，et al.Based on the analysis of the deformation area of permafrost subgrade height effect research[J].Chinese and Foreign Road,2013(3):22-28.

(責任編輯 徐志宏)

Stability Analysis of Surface Deformation Monitoring Data of Backfill Mining in Cold Region

Guo Qingbiao1,2,3Guo Guangli1,2,3Wang Jintao1,2,3

(1.SchoolofEnvironmentScienceandSpatialInformatics，ChinaUniversityofMiningandTechnology，Xuzhou221116;2.NASGKeyLaboratoryofLandEnvironmentandDisasterMonitoring，Xuzhou221116，China;3.JiangsuKeyLaboratoryofResourcesandEnvironmentInformationEngineering，Xuzhou221116，China)

In order to accurately obtain filling mining surface deformation data，it is necessary to test the quality of monitoring data at each period.The results show that control points are stable during monitoring through F-test and t-test for the stability analysis of surface deformation monitoring data of filling in cold region.Then，further analysis indicates that there are three kinds of vertical movements：the first is the surface subsidence caused by mining；the second is roadbed subsidence under long-term car loads；the last one is thaw collapse of seasonal frozen soil.The last is more general in the high latitude regions of China，and this phenomenon will have some influence on the surface deformation monitoring data，which enlarges the surface movement boundary and reduces the main affecting tangent angle.It is not conducive for getting estimated parameter accurately for this kind of geological and mining conditions.

Backfill mining，Statistical tests，Stability analysis，Seasonal frozen soil

2015-02-02

“十二五”國家科技支撐計劃項目(編號：2012BAB13B03)，國家自然科學基金青年基金項目(編號：41104011)，江蘇省高校優勢學科建設工程項目(編號：SA1102)，中央高?；究蒲袠I務專項(編號：2012LWB32)。

郭慶彪(1990—)，男，博士研究生。

TD17

A

1001-1250(2015)-04-220-04