厚松散層下條帶充填開采沉陷影響因素的數值模擬分析

張 磊 張兆民

(1. 大同大學，山西省大同市，037004； 2. 山東省煤田地質規劃勘察研究院，山東省濟南市，250100； 3. 山東省煤田資源數字化工程技術中心，山東省濟南市，250100)

隨著煤炭資源長期無節制的開采，中國東部老礦區的煤炭資源正面臨趨于枯竭的局面，并且尚未開采的煤炭資源大部分屬于“三下”壓煤。賦存有厚松散層的礦區，其地表沉陷規律具有一定的特殊性(松散層含水時，其下沉系數偏大)，通過采礦工作者歷年的實踐研究，對于普通地質條件下的煤炭開采引起的地表移動變形已經進行了大量深入的研究，然而特殊地質條件下開采，如厚松散層開采引起的地表下沉規律研究較少。

煤礦開采過程中地表沉陷造成環境破壞，房屋出現裂縫等災害現象，成為需要面臨的熱點問題。為此，部分學者提出了合理有效的開采方式來控制災害的發生，條帶充填開采成為一種可行的綠色開采方式。條帶充填開采采用部分充填的方法，在充填強度有保證的前提下，可以有效地控制地表大范圍下沉，減小開采對環境造成的破壞。

本文針對工作面上方賦存厚松散層的情況，運用FLAC3D三維數值模擬軟件，研究厚松散層下條帶充填開采影響地表沉陷的因素及其規律，為厚松散巖層下開采提供理論基礎，可以大大提高采出率、降低充填成本、將經濟效益最大化。

1 數值模型及模擬方案

1.1 數值計算模型

依據山東義能煤礦實測地質條件，確定模型參數，建立數值模擬基本模型，如圖1所示。為了簡

化數值模型，模型長度700 m，寬度220 m，因選取上覆巖層的組合結構不同，模型高度介于254～334 m之間。模型的上邊界是地表，下邊界是煤層底板。為了消除邊界效應的影響，走向方向上開采邊界以外設計留設200 m邊界條件，傾斜方向留設60 m邊界條件，即充填工作面的寬度為100 m，開采厚度3 m。在模擬計算過程中，采用摩爾—庫倫強度準則。巖體物理參數見表1。

此次研究地表沉陷影響因素的問題，可以看作平面移動應變問題，故假設工作面推進方向無限長，所研究的移動盆地的主斷面前后兩個面(XZ面)限制Y方向位移，左右兩個面(YZ面)限制X方向位移，底部限定垂直方向位移和水平方向位移。

圖1 計算模型網格劃分

巖層名稱體積模量/GPa切變模量/GPa粘聚力/MPa內摩擦角/(°)單軸抗拉強度/MPa容重/kg·m-3覆蓋土004012017190041880粘土2002009055180011980粘土1004014058190021980砂2005007027150031970砂1003005017180011970砂礫2004007019170021900砂礫1005006017180011900粗砂巖3655585658336452200中砂巖2590530552335792500粗砂巖2555385618316452200中砂巖1550550652326192500粗砂巖1555485718336452200泥巖510255320303552350煤層185130215292401500砂巖2560420953349352600砂巖158043098535995260020MPa充填體14105005528012220025MPa充填體25009010528017220030MPa充填體35212513728025220035MPa充填體45516717528032220040MPa充填體560210220280412200

1.2 數值模擬方案

本文共設計30個方案，其中20個方案采用單一控制變量法進行模擬研究，分別對不同松散層厚度(160 m、180 m、200 m、220 m、240 m)、不同充填率(51%、55.5%、60%、64.5%、69%)、不同充填體強度(2.0 MPa、2.5 MPa、3.0 MPa、3.5 MPa、4.0 MPa)以及不同關鍵層厚度(40 m、55 m、70 m、85 m、100 m)下對條帶充填開采進行地表下沉規律研究，模擬方案為充填率60%、充填體強度2.0 MPa、關鍵層厚度70 m、松散層厚度200 m。模擬方案見表2。

表2 模擬方案

圖2 工作面條帶充填示意圖

工作面條帶充填開采示意圖如圖2所示。由圖2可知，在推進方向上留設一定的充填條帶，充留間隔為60 m。充填率的計算是考慮煤炭采出后頂底板移近量和施工工藝(充填體的接頂率和泌水率)等因素影響后的實際充填率，其計算表達式為：

(1)

其中，實際充填高度取2.7 m。在模型頂部(地表)垂直于工作面推進方向的中央(X=0～700，Y=120，Z=地表值)，每隔10 m設置一個監測點，共設計70個監測點，用來監測開采過程中地表下沉的垂直位移(Z方向)和橫向水平位移(X方向)，因縱向水平位移(Y方向)不在下沉盆地主斷面內，故暫不研究。

2 數值模擬分析結果

2.1 充填率對地表下沉的影響

充填率不同時，其下沉盆地主斷面的地表下沉曲線及最大下沉值的變化趨勢分別如圖3和圖4所示。

圖3 地表下沉曲線

圖4 地表最大下沉值變化趨勢圖

由圖3可以看出，各方案下沉盆地主斷面的地表下沉曲線關于地表最大下沉值點呈近似對稱分布形態，最大下沉點位于采空區下沉盆地中央，下沉盆地邊界相較于采空區邊界向外擴展約100 m。下沉盆地的最大下沉值分別為581 mm、520 mm、461 mm、404 mm和298 mm。在保持充填體強度的前提下，設計方案對下沉系數進行了控制，控制在0.1～0.19之間。由模擬結果可知，地表最大下沉量較小，表明條帶充填體—覆巖主關鍵層組合結構對控制覆巖運動以及減小地表下沉起到重要作用。

由圖4可以看出，隨著充填率的不斷增加，地表最大下沉值逐漸減小，其中充填率從64.5%增大到69%，最大值減小較明顯。主要原因在于充填率增加后，相鄰條帶之間的間隔減小，采空區上方頂板的懸跨度減小，且充填體內部形成的彈性核區域增加，穩定性和強度增加。條帶充填體自身變形量較小，減少了上覆巖層運動的空間，能有效控制關鍵層彎曲下沉，因此，地表下沉量減小。

2.2 松散層厚度對地表下沉的影響

對于不同的松散層厚度條帶充填開采方案，其下沉盆地主斷面的地表下沉曲線及最大下沉值的變化趨勢如圖5和圖6所示。

圖5 地表下沉曲線圖

圖6 地表最大下沉值變化趨勢圖

由圖5可以看出，各方案下沉盆地主斷面的地表下沉曲線關于地表最大下沉值點呈近似對稱分布形態，最大下沉點位于采空區下沉盆地中央，下沉盆地邊界相較于采空區邊界向外擴展了100～180 m左右，隨著松散層厚度的增加地表下沉范圍逐漸增大，同時地表下沉盆地邊緣收斂逐漸變緩。下沉盆地的最大下沉值分別為538 mm、492 mm、461 mm、386 mm和401 mm。在保持充填體強度的前提下，設計方案對下沉系數進行了控制，控制在0.12～0.18之間，地表整體下沉值較小，表明條帶充填開采能有效控制上覆巖層運動減緩地表下沉。

由圖6可知，當松散層厚度在160～180 m范圍內時，隨著松散層厚度的增加地表最大下沉值小幅度減小。當松散層厚度在180～240 m范圍內時，地表最大下沉值變為大幅度增大趨勢。主要原因在于當松散層厚度在160～180 m范圍內時，厚松散層對地表下沉有一定的緩解作用；在180～240 m范圍內，厚松散層對地表下沉的緩解作用弱于松散層增加的自重對條帶充填體壓縮作用以及覆巖主關鍵層的彎曲變形，因此地表下沉表現為隨松散層厚度的增加而增加。

2.3 充填體強度對地表下沉的影響

對于不同的充填體強度條帶充填開采方案，其下沉盆地主斷面的地表下沉曲線及最大下沉值的變化趨勢如圖7和圖8所示。

圖7 地表下沉曲線

由圖7可以看出，各方案下沉盆地主斷面的地表下沉曲線關于地表最大下沉值點呈近似對稱分布形態，最大下沉點位于采空區下沉盆地中央，下沉盆地邊界相較于采空區邊界向外擴展超過100 m。下沉盆地的最大下沉值分別560 mm、488 mm、459 mm、431 mm和403 mm。在保持充填體強度的前提下，設計方案對下沉系數進行了控制，控制在0.13～0.19之間，由模擬結果可知，地表最大下沉量較小，表明條帶充填體—覆巖主關鍵層組合結構對控制覆巖運動以及減小地表下沉起到重要作用。

圖8 地表最大下沉值變化趨勢圖

由圖8可以看出，當充填體強度在2.5～4.0 MPa時，充填體強度變化趨勢和地表下沉相反，充填體強度增加，下沉值減小，主要原因在于在充填體強度能滿足支承上覆巖體重量時，充填體強度越大抵抗變形的能力就越強，當受到上覆巖層壓縮時，自身變形量較小，關鍵層不能充分運動，地表發生下沉就小。充填體強度由2.5 MPa減小到2.0 MPa時，地表下沉值表現出奇異性，下沉值增加幅度較大，因為充填強度較低時，在上覆巖層的重力作用下充填體狀態發生變化，由彈性狀態變為塑性狀態，使充填體自身變形量增加，因此，地表下沉量隨之增加。

2.4 關鍵層厚度對地表下沉的影響

對于不同的關鍵層厚度條帶充填開采方案，其下沉盆地主斷面的地表下沉曲線及最大下沉值的變化趨勢如圖9和圖10所示。

由圖9可以看出，各方案下沉盆地主斷面的地表下沉曲線關于地表最大下沉值點呈近似對稱分布形態，最大下沉點位于采空區下沉盆地中央，下沉盆地邊界相較于采空區邊界向外擴展超過100 m。下沉盆地的最大下沉值分別為705 mm、509 mm、461 mm、271 mm和325 mm。在保證充填體穩定的前提下，本組條帶充填開采試驗方案將主斷面地表下沉系數控制在0.19～0.24之間，地表整體下沉值較小，表明條帶充填開采能有效控制上覆巖層運動減緩地表下沉。

圖9 地表下沉曲線

圖10 地表最大下沉值變化趨勢圖

由圖10可以看出，當關鍵層厚度小于55 m時，地表最大下沉值隨關鍵層厚度的減小顯著增大；當關鍵層厚度為55～70 m時，隨關鍵層厚度的增加地表最大下沉值減小幅度較小，而當關鍵層厚度為70～85 m時，其減小幅度較大；當關鍵層厚度為85～100 m時，地表最大下沉值隨厚度的增加而增大。主要原因在于當關鍵層厚度在一定范圍內，并且保證充填體穩定性的前提下，關鍵層越厚覆巖主關鍵層強度越大，其自身抵抗形變的能力越強，能承擔上覆巖層及其自身的重量，起到關鍵層作用，下沉量較小。當關鍵層厚度逐漸增加時，上覆巖層的自重也開始變大，條帶充填體承受的實際載荷不斷增大，當載荷增大到一定程度時，條帶充填體將會由彈性狀態向塑性狀態轉變，隨著塑性區不斷擴大，充填體自身的壓縮量也隨之增加，當條帶充填體的壓縮量超過關鍵層對上覆松散層下沉的控制量時，地表下沉量表現出增加的趨勢。

3結論

(1)充填率變化規律與地表最大下沉值呈負相關，充填率越大，下沉量越小，且當充填率由64.5%增長至69%時，最大下沉值減小幅度較大。

(2)隨著松散層厚度增加，地表移動變形的下沉值先減小后增大，當松散層厚度為160～180 m時，隨著松散層厚度的增加地表最大下沉值逐漸減??；當松散層厚度為180～240 m時，隨著松散層厚度的增加地表最大下沉值逐漸增大。

(3)充填體強度變化規律與地表最大下沉值呈負相關，充填強度越大，下沉量越小。

(4)當關鍵層厚度在40～85 m范圍內，并且保證充填體穩定性的前提下，隨著關鍵層厚度的增加，地表移動變形的下沉值逐漸減??；當關鍵層厚度超過85 m時，地表移動變形的下沉最大值由減小趨勢變為增大趨勢。

[1] 賈民，徐法奎. “三下”壓煤開采方案優化設計[J]. 煤礦開采，2012(4)

[2] 許延春. 礦區深厚復合含水松散層的工程、力學特性及其應用 [D]. 煤炭科學研究總院，1999

[3] 王同旭，馬磊，馬文強. 煤柱群穩定性綜合指數法與模糊評價法綜合分析 [J]. 山東科技大學學報(自然科學版)，2015(4)

[4] 劉貴，張華興，徐乃忠. 深部厚煤層條帶開采煤柱的穩定性 [J]. 煤炭學報，2008(10)

[5] 劉曉明，趙同彬，王明強等. 固體充填工作面支架工作特性及頂板控制分析 [J]. 山東科技大學學報(自然科學版)，2017(2)

[6] 陳俊杰，鄒友峰，郭文兵．厚松散層下下沉系數與采動程度關系研究[J]. 采礦與安全工程學報，2012(2)