基于格網法的地面復墾模擬與分析

陳慧玲 肖 武 王 錚 李素萃 臺曉麗

（中國礦業大學（北京）土地復墾與生態重建研究所，北京市海淀區，100083）

基于格網法的地面復墾模擬與分析

陳慧玲 肖 武 王 錚 李素萃 臺曉麗

（中國礦業大學（北京）土地復墾與生態重建研究所，北京市海淀區，100083）

為實現地面復墾與井下開采的充分耦合，探索一種新的土地復墾方法，以我國東部高潛水位礦區為研究對象，采用格網法將地面劃分成方格單元，建立煤炭開采影響地面方格單元的煤層范圍的計算模型，基于該模型圈定影響方格單元A（i，j）的井下煤層范圍，進而采用概率積分法進行開采沉陷模擬，分析各開采階段方格單元A（i，j）的沉陷特點與規律，得出以下結論：方格單元A（i，j）應在開采第二階段采取剝離表土的措施，并通過地面—地下—地面的土地復墾方法，科學劃定復墾區域，實現井上井下的充分耦合。

采煤塌陷地 土地復墾 高潛水位 格網法

1　前言

土地復墾工程施工時，為了便于施工及工程量計算，多以塊段作為施工單元，但目前已開展的未穩沉塌陷地復墾時機的研究多是以地下煤層開采為前提的地面復墾時機的確定。因此，地面復墾施工單元何時何地與如何復墾仍是困擾復墾工作開展與現場復墾施工的難題。本文以地面為出發點，將復墾區域格網化后，開展基于格網法的地面方格單元的復墾模擬與分析的研究，將為指導復墾施工提供科學依據。

2　研究區概況

選擇安徽省淮北市某高潛水位礦區進行模擬研究，研究區屬于華北平原，地勢較平緩。地表自然標高+30～+32 m，平均標高約為+31 m，地面坡度為0°～2°，地下潛水位埋深約為1.5 m。礦井于1993年12月26日正式投產，生產能力為200萬t/a，礦區煤層結構簡單，以單一近水平煤層為主，試驗區位于礦井的東部，煤層平均埋藏深度為250 m，煤層厚度平均為2 m。井下分布有4個工作面，工作面尺寸為800 m×130 m，開采方式為順序開采。示范區氣候溫和，屬北方大陸性氣候與濕潤氣候之間的季風氣候，日照充足，四季分明，春秋季明顯短于冬夏季，冬季寒冷干燥、夏季炎熱多雨。

3　分析材料與方法

以往的開采沉陷研究與預測技術都是從地下出發，先研究地下單個開采單元對地面造成的影響程度及范圍，再通過積分疊加等方式研究地下任意開采尺寸對地面造成的影響。土地復墾以地面為主，人們較為關心的是地面點會受到地下哪些區域煤層開采的影響。基于此，首先由地面到地下，根據邊界角圈定地下煤炭開采可能影響地面方格單元的煤層范圍，再反過來由地下到地面，分析地表各方格單元的沉陷特點與規律，確定各方格單元適宜的復墾時間與復墾方向，指導地面復墾施工。

3.1地面方格網劃分

根據研究區實地情況，將沉陷區劃分成規則排列的正方形格網，每個格網作為復墾施工的最小單元。研究區面積為126 hm2，由于地勢較為平緩，便于復墾施工，因此將沉陷區劃分為8行7列大小為150 m×150 m的方格網，使用阿拉伯數字從上到下、從左到右依次編號，使任一方格單元編號為A（i，j），表示第i行j列的小方格網。

3.2地面方格單元的井下煤層影響范圍的確定

開采沉陷學中規定，在充分采動或接近充分采動的條件下，地表移動盆地主斷面上盆地邊界點至采空區邊界的連線與水平線在煤柱一側的夾角為邊界角。根據邊界角建立地面點受煤炭開采的煤層影響范圍的模型。

在走向方向上，設置地面某一方格單元A（i，j），幾何中心點坐標為O（xi，yi，zi），2R范圍內的煤層開采可能會對O點造成影響，R的計算公式為:

式中:R——地面點走向方向上煤層影響半徑，m；

H——地面點到煤層的鉛垂距離，m；

δ——走向邊界角，（°）。

定義角度θ為能影響地面點O的煤層范圍角，角度θ的值為走向邊界角δ的余角，如圖1（a）所示。

同理，在傾向方向上，R1與R2范圍內的煤層開采可能會對O點造成影響，R1與R2的計算公式為:

式中:R1——地面點傾向下山方向上煤層影響半徑，m；

R2——地面點傾向上山方向上煤層影響半徑，m；

β——下山邊界角，（°）；

γ——上山邊界角，（°）。

定義角度θ1為能影響地面點O的煤層范圍下山角，角度θ2為能影響地面點O的煤層范圍上山角，角度θ1的值為上山邊界角γ的余角，角度θ2的值為下山邊界角β的余角，如圖1（b）所示。

根據邊界角的余角可以確定影響地面點O的井下煤層范圍。以地面方格單元A（5，4）為例，確定影響地面各方格單元的煤層范圍。研究區煤層為單一近水平煤層，方格單元A（5，4）的幾何中心點為O，O點到煤層的鉛垂距離H取值為240 m，邊界角δ為62.4°，將H、δ值代入公式（1）計算R為125 m，煤層范圍角θ為δ的余角27.6°。

對方格單元A（5，4）的4個頂點a、b、c和d以27.6°角向煤層所在水平面進行投影，得到煤層所在水平面內以點O′為圓心，以過a′、b′、c′和d′4個頂點為內接四邊形的區域，該區域所包含的范圍即為影響地面方格單元A（5，4）的煤層區域，如圖2所示。

3.3地面方格單元的沉陷規律分析

將確定影響地面各方格單元的煤層范圍與井下工作面疊加，得到位于該煤層范圍內的工作面數量及空間分布情況。方格單元A（5，4）的井下煤層影響范圍內分布有4個工作面，左右呈對稱狀態分布，A（5，4）主要受2和3工作面煤炭開采的影響，如圖2所示。煤層平均埋藏深度為220～270 m，煤層平均厚度為2 m，屬中厚煤層，工作面的開采方式為順序開采，開采方向為1→2→3→4。按照一個工作面為一個階段，可將影響方格單元A（5，4）的煤層區域分為4個開采階段。

圖1　地面點受地下煤層開采影響圖

圖2　地面方格單元受地下煤層開采影響圖

根據礦區的自然條件、地質條件和采礦計劃等信息，基于概率積分法，采用沉陷預計軟件預測各方格單元在4個開采階段的影響，然后分析各方格單元的沉陷規律。第一階段開采結束后，方格單元A（5，4）受采動影響較小，地表略有下沉，無積水產生，只需采取簡單的平整措施即可恢復土地的使用；第二階段開采結束后，方格單元A（5，4）受采動影響較大，下沉值已達到1.61 m，大于地下潛水位埋深，地面產生季節性積水，此時若不采取復墾措施，肥沃的土壤將沉入水中，喪失其利用價值；第三階段開采結束后，方格單元A（5，4）的沉陷程度進一步加劇，整個方格單元形成常年積水區，陸生生態系統退變為水生生態系統，此時再復墾，不僅喪失珍貴的表土資源，而且增加復墾難度，復墾率也大大降低；第四階段開采結束后，方格單元A（5，4）達到最終的沉陷狀態，積水程度進一步加劇。因此，第二階段是方格單元A（5，4）最佳的復墾時間。方格單元A（5，4）各開采階段積水面積如圖3所示。

圖3　方格單元A（5，4）各開采階段積水面積圖

地面沉陷模擬分析時，通常將地下潛水位埋深作為判斷土地積水與否的分界線，將與潛水位埋深相等的下沉等值線作為臨界積水線，根據臨界積水線確定各方格單元的復墾方向，位于臨界積水線內側的方格單元應采取剝離表土措施；位于臨界積水線外側的方格單元應采取表土充填措施；位于臨界積水線上的方格單元應根據充填或剝離面積所占的比重確定復墾措施。方格單元A（5，4）位于臨界積水線內側，因此地面應在第二階段采取剝離表土措施。動態復墾時，復墾施工順序通常與工作面的開采順序相一致，各階段復墾范圍的劃分應以最大限度的保護表土為目的。將方格網劃分為3個復墾階段，復墾順序為從左到右依次剝離表土，各復墾階段表土剝離范圍如圖4所示。

圖4　方格網復墾階段與復墾方向圖

4　地面方格單元的井下工作面開采順序的優化

地面受開采影響最終的沉陷損毀狀態雖然相同，但不同的開采順序造成土地的損毀時序不一致，所以適當調整地下工作面的開采順序，可以延長土地的使用時間。分別采用順序開采（1-2-3 -4）、跳采（1-3-2-4）和兩端逼近式開采（1 -4-2-3）三種方式模擬地面方格單元A（5，4）的土地變化過程。對比上述三種開采方式的沉陷結果，如圖5所示，跳采與順序開采的沉陷過程基本一致，略有差異的是第二階段開采結束后地面的沉陷值為1.52 m，略小于順序開采的沉陷值1.61 m；兩端逼近式開采土地的沉陷過程不同于順序開采與跳采，第一、二階段開采結束后，地面的沉陷較小，第三階段地面才有明顯的下沉，地面的沉陷值達到順序開采第二階段的沉陷值，地面產生積水，因此，第三階段是兩端逼近式開采最佳的復墾時間。綜合上述分析結果，方格單元A（5，4）若調整工作面的開采順序，采用兩端逼近式開采（1-4-2-3）方式，可在第三開采階段再剝離表土，延長了土地的使用時間。

圖5　不同開采順序下方格單元A（5，4）沉陷路徑圖

5　結論

本文提出一種新的土地復墾確定方法，通過該方法能科學準確的分析地面各方格單元的復墾時間與復墾方向，指導復墾施工。

（1）從地面到地下，根據地面點受煤炭開采的煤層影響范圍的計算模型，以邊界角的余角27.6°圈定影響方格單元A（5，4）的井下煤層范圍，再由地下到地面，對方格單元A（5，4）進行沉陷模擬，分析各開采階段方格單元的沉陷規律。

（2）方格單元A（5，4）第二開采階段下沉值1.61 m大于地下潛水位埋深1.5 m，地面將有積水產生，所以第二階段是最佳的復墾時間，進而根據臨界積水線1.5 m確定方格單元A（5，4）應采取剝離表土的復墾措施。

（3）進一步探討調整工作面采序以延長土地的使用時間的可能性，通過對比分析方格單元A（5，4）三種不同開采順序下的地面的損毀時序，發現若采用兩端逼近式（1-4-2-3）開采地下煤層，可在第三開采階段采取復墾措施，延長土地的使用時間。

綜上所述，基于地面—地下—地面的三步驟分析方法，能夠快速簡便的確定地面任一方格單元受地下開采影響的煤層范圍，科學準確的分析地面方格單元在開采影響下的下沉規律并安排復墾，實現了地面復墾時序與井下煤層開采時序的充分耦合，為邊開采邊復墾提供技術支持。

［1］ 楊光華，胡振琪，趙艷玲.高潛水位采煤沉陷地復墾治理對策研究［J］.煤炭工程，2014（6）

［2］ 胡振琪，李晶，趙艷玲.礦產與糧食復合主產區環境質量與糧食安全的問題、成因與對策 ［J］.科技導報，2006（3）［3］ 肖武，胡振琪，李太啟等.采區地表動態沉陷模擬與復墾耕地率分析［J］.煤炭科學技術，2013（8）

［4］ 董祥林，陳銀翠，歐陽長敏.礦區塌陷地梯次動態復墾研究［J］.中國地質災害與防治學報，2002（3）

［5］ 張友明，代曉東.動態沉陷復墾在高潛水位采煤塌陷區中的應用［J］.礦山測量，2004（3）

［6］ 胡振琪，肖武，王培俊等.試論井工煤礦邊開采邊復墾技術［J］.煤炭學報，2013（2）

［7］ 肖武.井工煤礦區邊采邊復的復墾時機優選研究［D］.中國礦業大學（北京），2012

［8］ 馬占林，李積蘭，余筱蓉.土石方量計算方法比較［J］.地礦測繪，2010（1）

［9］ 許軍，李笑一，孫彩敏.我國礦區土地復墾的主要問題及其對策［J］.中國煤炭，2010（12）

［10］ 何廳廳，趙艷玲，石娟娟等.基于可拓法的高潛水位礦區復墾潛力評價［J］.中國煤炭，2014（4）

（責任編輯 孫英浩）

Simulation and analysis of land reclamation basing upon grid method

Chen Huiling，Xiao Wu，Wang Zheng，Li Sucui，Tai Xiaoli
（Institute of Land Reclamation and Ecological Restoration，China University of Mining and Technology，Beijing，Haidian，Beijing 100083，China）

For sufficient coupling between ground reclamation and underground mining，the authors explored a new method of land reclamation which took a mining area of China eastern high ground-water level area as subject of the study.The ground was divided into grid cells by grid method.The authors established computational model of coal seam region which influential the ground grid cells during mining.Basing upon the model，influential underground coal seam region of grid cell A（i，j）was delineated.Then，the authors conducted mining subsidence simulation by using probability integral method and analyzed subsidence characteristics and patterns of grid cell A（i，j）at each mining stage.The results showed that grid cell A（i，j）should take measures of topsoil stripping at second stage.According to land reclamation method of ground-underground-ground，reclamation region was scientifically delimited and ground and underground were achieved sufficient coupling between ground reclamation and underground mining.

coal mining subsidence area，land reclamation，high underground water level，grid method

TD88

A

陳慧玲（1989-），女，山東曹縣人，碩士研究生，研究方向為礦區土地復墾與生態重建。