陳向紅，　陶連金，　趙志榮

(1. 北京工業(yè)大學 建工學院， 北京 100124； 2. 中國礦業(yè)大學 安全學院， 北京 100083)

?

地表鉆孔注漿工藝下采空區(qū)覆巖力學性態(tài)的數(shù)值模擬

陳向紅1，陶連金1，趙志榮2

(1. 北京工業(yè)大學 建工學院， 北京 100124； 2. 中國礦業(yè)大學 安全學院， 北京 100083)

圍巖變形大、控制困難等問題一直以來都是采礦界著力研究的技術課題，尤其在薄巖體分布的特殊地質類工程中表現(xiàn)更加突出。以石圪臺煤礦22303回采工作面為例，采用數(shù)值模擬方法，研究在地表鉆孔全段注漿工藝下采空區(qū)覆巖位移和應力變化規(guī)律。結果表明：計入圍巖彈塑性變形特性，考慮作業(yè)面前方圍巖的空間約束效應，在距離作業(yè)面1～1.5倍采高范圍內覆巖存在較大的豎向變形和應力差，并隨距作業(yè)面距離應力逐步釋放，位移逐漸增大；當回采高度遠小于回采長度時，可不計采煤對側墻水平位移和應力的影響。利用切面等值線分析了覆巖內的力學性態(tài)，研究成果豐富了地表鉆孔注漿下的采空區(qū)覆巖力學理論，對同類條件下的覆巖變形控制具有重要的借鑒意義。

地表鉆孔注漿；采空區(qū)；覆巖；力學性態(tài)；數(shù)值模擬

0　引　言

作為國民經濟的支柱產業(yè)，煤炭在我國能源結構中仍占據(jù)著主導地位。隨著煤炭產量的逐年增長及城市環(huán)境的不斷惡化，安全和綠色開采受到學術界與工業(yè)領域的高度關注[1,2]。

目前，國內在控制采煤引起圍巖變形的研究，主要側重于充填采煤覆巖移動規(guī)律方面，如朱衛(wèi)兵等[3]提出覆巖離層分區(qū)隔離注漿充填技術，并采用數(shù)值模擬方法做了進一步證實；李興尚等[4]應用PFC顆粒流程序模擬條帶煤層整個動態(tài)變形發(fā)展過程；Xuan Dayang等[5]分析了深厚火成巖的離層注漿對覆巖穩(wěn)定性的影響；張吉雄[6]、陳杰等[7]運用彈性地基梁理論分析了充填材料彈性模量與覆巖撓度之間的影響關系；黃艷利[8]綜合理論建模、力學分析、實驗測試和數(shù)值模擬等多種方法系統(tǒng)研究了固體密實充填的巖層移動特征和控制機理；李猛等[9]基于彈性地基薄板理論建立了密實充填采煤充填體-基本頂力學模型，并推導了基本頂破斷臨界條件。上述研究，主要針對采空區(qū)和巖石離層的加固處理措施，對于上覆松散層全段注漿工藝[10]下采煤覆巖的變形和應力演化規(guī)律尚鮮見報道。鑒于此，筆者以石圪臺煤礦22303注漿工作面為例，通過數(shù)值模擬方法，研究地表深孔全段注漿采煤覆巖的力學性態(tài)，以期為注漿采煤覆巖的變形控制提供理論參考。

1　工程描述

石圪臺煤礦22303工作面具有上覆巖層厚度較薄松散層厚度較大的特點，回采至該區(qū)域時存在冒頂?shù)扔绊懓踩┕さ奈ｋU。為了控制采煤施工引起的覆巖變形過大，場區(qū)提出了“地表鉆孔全段注漿工藝”。針對采煤工作面基巖厚度小于30 m的上覆松散層進行注漿加固處理，設計注漿段長度為將基巖厚度加固至約40 m位置。

依據(jù)鉆孔提供資料，鉆孔設計間距為6.0 m，呈正方形6 m×6 m布置，見圖1。

注漿工程范圍內揭露地層由老到新分布為：細粒砂巖、煤巖(J1-2y)，細粒砂巖、泥巖(J2z)，淤泥質粉砂、粉土(Q4)。地層傾角平緩，近水平狀分布。

圖1　工作面注漿孔布置示意

2　采空區(qū)覆巖力學性態(tài)的數(shù)值模擬

2.1計算模型的建立

在建模中，綜合考慮模型幾何特征、收斂性和精度選擇單元類型和單元網(wǎng)格尺寸。模型坐標x、y、z分別對應回采掘進方向、采煤工作面寬度方向和地層深度，對應坐標方向模型尺寸36 m×72 m×88 m，采空區(qū)側邊界距離模型邊界距離取3倍采高，利用有限元軟件ANSYS建立實體單元進行模擬。模型邊界條件為：沿x和y方向邊界均限值水平位移，底部邊界約束豎向位移，上邊界施加豎向壓應力模擬未注漿松散層自重應力荷載。計算時按3 m增量模擬回采長度變化，利用軟件提供的單元“生死”功能實現(xiàn)采煤的動態(tài)施工過程。三維計算模型如圖2a所示。圖2b給出了用來分析不同采長下覆巖力學性態(tài)的四個觀測點。其中Point1～ Point4分別為采空區(qū)頂板軸線上距離始發(fā)端為6、12、18和24 m的四個點。

a

b

2.2材料參數(shù)

根據(jù)采場地質勘查報告，選用Drucker-Prager彈塑性模型作為計算分析時圍巖的本構模型，經過對材料物理力學性質相近地層的合并與等效，各地層參數(shù)如表1所示。

表1　地層劃分及參數(shù)

3　數(shù)值計算與結果分析

為了考察采空區(qū)覆巖力學性態(tài)的變化規(guī)律，文中選取覆巖厚度為20和30 m兩種情況下的計算結果進行分析。

3.1各觀測點的時空效應分析

經計算，四個觀測點豎向位移uz和σz應力隨回采長度l的變化規(guī)律如圖3所示。

a　豎向位移

b　豎向應力

由圖3a可知，不同的覆巖厚度下，各觀測點豎向位移均隨回采長度的增加而增大，其中，覆巖厚度為20 m、回采長度24 m時，觀測點1～4的豎向位移分別為20.5、24.8、27.4和28.3 mm；覆巖厚度為30 m時，各觀測點豎向位移分別為17、19.5、20.9和21.3 mm。薄巖層更容易產生較大的變形。

由圖3b可知，各觀測點均在工作面回采至該位置時出現(xiàn)應力極大值，之后隨回采長度逐漸衰減，最大應力幅值隨觀測點由Point1到Point4呈指數(shù)形式增長。

3.2沿回采方向頂板軸線豎向位移和應力

為了考察采空區(qū)覆巖的力學性態(tài)。圖4給出了回采長度lh為24 m時采空區(qū)頂板軸線豎向位移和應力變化規(guī)律。

a　豎向位移

b　豎向應力

由圖4可知，在回采工作面附近，由于受到前方圍巖的空間約束效應，覆巖應力逐步釋放，所以在作業(yè)面附近形成較大的應力差降；當離開作業(yè)面約1.0～1.5倍的采高時，覆巖不再受圍巖約束，豎向位移隨之增加。同時，由兩種覆巖厚度下的位移和應力曲線發(fā)現(xiàn)，覆巖厚度的變化對位移的影響程度明顯高于應力。

3.3沿回采方向側墻水平位移和應力

作為采空區(qū)覆巖支撐，側墻圍巖的力學性態(tài)也直接影響到覆巖的位移和應力發(fā)展。圖5給出了工作面回采在橫斷面內對采空區(qū)側墻水平位移uy和σy應力的擾動影響。

由圖5可知，采空區(qū)側墻水平位移和應力相比覆巖豎向位移要小，主要原因可歸結為采空區(qū)的空間高跨比較小，圍巖具備上覆地層載荷作用下的自承能力。受圍巖的約束影響，水平位移和應力差亦出現(xiàn)在距離作業(yè)面約1.0～1.5倍采高范圍內，之后逐漸趨于穩(wěn)定。

a　水平位移

b　水平應力

3.4覆巖內的力學性態(tài)模擬

為了更直觀的描述覆巖內的力學性態(tài)，圖6～8分別示出了沿模型頂板軸線截面的豎向位移、豎向應力和剪應力等值線圖。

a　h=30 mm

b　h=20 mm

a　h=30 mm

b　h=20 mm

a　h=30 mm

b　h=20 mm

由圖6、7可知，受回采工作面的空間約束影響，在采空區(qū)邊界存在較大的應力梯度，為應力集中區(qū)。之后在遠離邊界位置逐漸進入降低區(qū)，豎向位移增大，表現(xiàn)為梁板的結構撓曲變形特征。進而比較圖7和8兩種覆巖厚度下的應力等值線還可發(fā)現(xiàn)：當覆巖厚度為30 m時，松散層注漿對覆巖變形和應力的影響不明顯；而當覆巖厚度降至20 m時，覆巖內應力與注漿應力連成一體，注漿對覆巖的變形和應力發(fā)展已產生影響。該結論可進一步證實22303回采工作面所提出的薄覆巖注漿加固處理方案的可行性。

4　結　論

(1)煤炭回采工作面對采空區(qū)覆巖有空間約束效應，對于石圪臺22303工作面其影響范圍約為作業(yè)面前方1.0～1.5倍采高，在該范圍以外，覆巖豎向應力隨距作業(yè)面距離的增加而逐步釋放，豎向位移逐步增加。

(2)采空區(qū)高跨比較小時，煤炭回采對采空區(qū)橫斷面?zhèn)葔λ轿灰坪蛻Φ挠绊懴啾容^小，可協(xié)同覆巖共同形成空間承載環(huán)承擔上覆地層的自重載荷并傳遞給底板巖層。

(3)通過數(shù)值計算，當采用地表鉆孔注漿工藝加固地層時，需要在考慮不同厚度覆巖的力學性態(tài)后選擇注漿參數(shù)，以達到合理減沉目的。

[1]繆協(xié)興，錢鳴高．中國煤炭資源綠色開采研究現(xiàn)狀與展望[J]. 采礦與安全工程學報， 2009， 26(1)： 1-14.

[2]LIANG KAIWU, LIU JIKUN, WANG CUIXIA. The Coal mine accident causation model based on the hazard theory [J]. Procedia Engineering, 2011 (26)： 2199-2205.

[3]朱衛(wèi)兵， 許家林， 賴文奇， 等．覆巖離層分區(qū)隔離注漿充填減沉技術的理論研究[J]．煤炭學報，2007， 32(5)： 458-462．

[4]李興尚，許家林，朱衛(wèi)兵，等．條帶開采垮落區(qū)注漿充填技術的理論研究[J]． 煤炭學報， 2008， 33(11)： 1205-1210．

[5]XUAN DAYANG，XU JIALIN，ZHU WEIBING. Dynamic disaster control under a massive igneous-sill by grouting from surface boreholes [J]．International Journal of Rock Mechanics & Mining Sciences， 2014， 71(1)： 176-187．

[6]張吉雄，李 劍，安泰龍，等．矸石充填綜采覆巖關鍵層變形特征研究[J]．煤炭學報， 2010， 35(3)： 357-362．

[7]陳杰， 杜計平， 張衛(wèi)松， 等．矸石充填采煤覆巖移動的彈性地基梁模型分析[J]．中國礦業(yè)大學學報， 2012， 41(1)： 14-19．

[8]黃艷利．固體密實充填采煤的礦壓控制理論與應用研究[D]．北京： 中國礦業(yè)大學， 2012．

[9]李猛， 張吉雄， 姜海強， 等．固體密實充填采煤覆巖移動彈性地基薄板模型[J]．煤炭學報， 2014， 39(12)： 2369-2373．

[10]何朋朋．注漿段全段注漿的方法．中國， zh201310328853.3[P]. 2013-10-23．

(編輯徐巖)

Numerical simulation on mechanical properties of overlying strata above goaf under surface drilling grouting

CHENXianghong1，TAOLianjin1，ZHAOZhirong2

(1. College of Civil Engineering, Beijing University of Technology, Beijing 100124, China;2. College of Security, China University of Mining & Technology，Beijing 100083, China)

This paper is devoted to overcoming the notorious difficulties, such as a larger deformation in surrounding rocks and consequent difficult control, which have remained technical challenge facing scientific researchers of mining engineering, especially so in special geological engineering with deep and thin rock mass distribution. The study drawing on the example of the stone GeTai 22303 working face and using the numerical simulation method works toward the law underlying how the displacement and stress change occur in overlying strata above goaf in the case of the surface borehole subjected to an entire-section grouting process. The results demonstrate that the consideration of the elasto-plastic characteristics of surrounding rock and the space constraint effect of surrounding rock in front of the construction face justifies the occurrence of a large vertical differential deformation and stress in the overlying strata within the scope of 1-1.5 times mining height from the working face; a gradual stress release due to an increasing distance from the working face is accompanied by a gradually increasing displacement; and the effects of excavation on the horizontal displacement and stress along the side wall can be negligible when the coal mining height is far less than the length of extraction direction. The study ends with analyzing the mechanics characteristic in strata using cross-section contour, thus contributing to enriching the theory of mechanics of rock above goaf under the technology of surface drilling grouting, and providing an important reference for strata displacement control in the same condition.

surface drilling grouting; goaf; overlying rock; mechanical property; numerical simulation

2015-12-05

陳向紅(1977-)，男，河北省保定人，博士，研究方向：地下工程，E-mail:07872@bjut.edu.cn。

10.3969/j.issn.2095-7262.2016.01.002

TD32

2095-7262(2016)01-0005-05

A