膏體充填參數(shù)對上覆巖層破壞規(guī)律的數(shù)值模擬研究

張誠

（忻州市煤田地質(zhì)勘探隊,山西忻州 034000）

膏體充填參數(shù)對上覆巖層破壞規(guī)律的數(shù)值模擬研究

張誠

（忻州市煤田地質(zhì)勘探隊,山西忻州 034000）

傳統(tǒng)的條帶開采中，巷道開挖后擾動了巷道上方原巖應(yīng)力的分布狀態(tài)，巷道上方的頂板在二次應(yīng)力場的作用下會發(fā)生變形、破壞，嚴重時會導(dǎo)致結(jié)構(gòu)失穩(wěn)、垮塌和破壞。為研究不同的膏體充填參數(shù)對覆巖應(yīng)力分布規(guī)律的影響，本文運用數(shù)值模擬的手段以小屯礦條帶充填開采為背景，選取采空區(qū)頂板及煤柱塑性區(qū)破壞范圍、巷道頂板和煤柱的下沉量作為充填體的充填效果，得出不同充填體的充填率、充填體強度對上覆巖層破壞規(guī)律的影響。

膏體充填參數(shù)；巖層；數(shù)值模擬

1 數(shù)值模擬方法的確定

在巖土工程與地質(zhì)材料的力學(xué)行為數(shù)值模擬計算中以拉格朗日法為代表的FLAC軟件與離散單元法所代表的UDEC軟件是兩種行之有效的解決手段[1-3]。

FLAC3D是一種基于三維有限差分法以巖石力學(xué)為基礎(chǔ)理論，根據(jù)巖土或其他材料的物理力學(xué)參數(shù)及地質(zhì)構(gòu)造特性的數(shù)值分析方法。由于混合了單元離散模型，采用有限差分格式可以準確模擬多種材料與非規(guī)則區(qū)域邊界條件的力學(xué)動態(tài)行為，在外力載荷及邊界條件的約束下，同時可以準確的模擬材料的屈服、塑性流動、甚至大變形的應(yīng)力-應(yīng)變關(guān)系?？迓鋷Ш蛷澢鲁翈Ф粫霈F(xiàn)垮落破壞，頂板的上覆巖層可以按照連續(xù)介質(zhì)的問題研究。因此，選用FLAC3D數(shù)值模擬軟件進行膏體條帶充填的模擬計算，為膏體條帶充填提供一定的理論依據(jù)。

2 FLCA3D數(shù)值模擬方案的建立

2.1 基本假定

小屯煤礦屬于近水平煤層，因此在數(shù)值模擬計算時，煤層按水平考慮。為提高運算速度與精度，精確的分析采空區(qū)頂板及煤柱的受力狀態(tài)及應(yīng)力分布規(guī)律，將模型分為頂?shù)装迮c煤層共三層。模型的上邊界以上覆巖層的地應(yīng)力作為邊界加載條件，根據(jù)公式（1）得出模型上邊界的應(yīng)力：

式中：q模型上邊界的應(yīng)力，MPa；γ為為上覆巖層的容重，N/m3；g為重力加速度，g/cm3；h為煤層埋深，m。

煤層頂板為粗粉砂巖，底板為細粒砂巖。由于采空區(qū)頂板巷道使用錨噴網(wǎng)聯(lián)合支護，所以在模型力學(xué)參數(shù)上進行加強。由于巷道圍巖采用了錨網(wǎng)支護，所以在力學(xué)參數(shù)上進行了加強。模型網(wǎng)格剖分圖，見圖1。

模型在X方向上取110 m，在Y方向上取100 m。模型高40 m，其中頂板與底板為便于計算各取18 m，煤厚為4 m，位于Z=18～20 m之間。在煤層18～22 m的區(qū)間內(nèi)，布置4條充填巷道，同時留設(shè)3個隔離煤柱，同時為減輕邊界效應(yīng)的影響在模型的左右邊界在X方向各留設(shè)28 m的煤柱。

2.2 膏體材料與巖石的物理力學(xué)參數(shù)的選取

為提高數(shù)值模擬的準確性，使模擬結(jié)果更真實的采空區(qū)上覆巖層與煤柱的變形與破壞以及充填體的充填效果，模擬采用的煤巖層物理力學(xué)參數(shù)主要根據(jù)實驗室實測數(shù)據(jù)與相關(guān)理論研究成果得出，頂?shù)装宓牟此杀圈倘?.47；充填體的泊松比μ取0.25；煤體的泊松比μ取0.3。煤巖層的物理力學(xué)參數(shù)，如表1所示。

2.3 邊界條件的確定

小屯煤礦膏體條帶充填計算模型邊界條件確定如下：

1）水平與初始位移：模型底部邊界的水平與垂直的初始位移為零；而模型前后與左右的兩端的邊界水平初始位移均為零。

2）水平與垂直應(yīng)力：根據(jù)公式（2）可知，模型的底部施加上負巖層自重應(yīng)力的的等效載荷。式中：γ為上覆巖層的容重，N/m3；H為模型底部邊界的埋深，m；σs為垂直應(yīng)力，MPa。

由公式（3）得出，模型邊界的前后左右在水平方向上施加由自重應(yīng)力產(chǎn)生的側(cè)向應(yīng)力。

式中：λ為側(cè)壓系數(shù)；μ為泊松比；σs為垂直應(yīng)力，MPa。

3 不同充填率對上覆巖層的數(shù)值模擬分析

充填率與充填體強度是膏體充填的關(guān)鍵參數(shù)。為研究不同充填率及充填體強度對采空區(qū)頂板及煤柱的應(yīng)力分析，數(shù)值模擬不同充填參數(shù)對上覆巖層的應(yīng)力分布及塑性區(qū)破壞影響范圍分兩種情況：一是保持充填體強度不變，模擬不同膏體材料充填率（70%、80%、90%、100%）對工作面頂板與隔離煤柱的塑形區(qū)破壞范圍，以及頂板的垂直應(yīng)力于煤柱的垂直應(yīng)力。二是保持膏體材料的充填率不變，通過分析充填區(qū)域不同充填體強度的塑性區(qū)破壞范圍，得出不同充填體強度對采空區(qū)頂板及煤柱的充填效果。

隨著原煤不斷采出，采空區(qū)頂板在上覆巖層垂直應(yīng)力的作用下會發(fā)生沉降。

由于充填體的充填率是相對于采空區(qū)頂板沉降后的實際空間，因而充填體的體積要小于原煤采出的體積。根據(jù)配比試驗得出的最優(yōu)配比，其充填體強度為3.57 MPa，模擬在充填體強度確定的情況下，選擇不同的膏體材料充填率（70%、80%、90%、100%）對工作面頂板與隔離煤柱的塑形區(qū)破壞范圍，以及頂板的垂直應(yīng)力與垂直位移。不同充填率的塑性區(qū)破壞范圍及頂板的垂直應(yīng)力與垂直位移，見圖2-圖6。

由圖2-圖6可以看出，當(dāng)采空區(qū)完全不充填時，采空區(qū)的頂板與底板板會發(fā)生剪切與拉伸破壞，造成頂板垮落破壞與底鼓。隨著充填體的充填率由70%增加至100%，上覆巖層的垂直應(yīng)力對塑性區(qū)破壞影響幅度逐漸變小。對比充填率為90% 和100%的采空區(qū)塑性破壞范圍及頂板的垂直位移可知，兩種不同充填率的破壞范圍基本相同，前者的采空區(qū)頂板垂直位移比后者增加2.14%，而配比膏體材料的膨脹率為2.25%，完全可以滿足主動接頂要求。

由圖7-圖11可以看出,采空區(qū)頂板垂直位移隨充填率的增加呈線性減小。隨著充填率的不斷增加，充填體控制頂板的能力不斷增強，從而約束頂板與隔離煤柱的位移。由于在模擬的膏體材料不是隨著巷道的推進而隨采隨充，頂板在未充入膏體材料時已有部分下沉量。因此，在實際充填作業(yè)時應(yīng)最大限度的減小頂板下沉量。

4 結(jié)論

充填率與充填體強度是膏體充填的關(guān)鍵參數(shù)。通過研究不同充填率及充填體強度對采空區(qū)頂板及煤柱的應(yīng)力分析得出當(dāng)采空區(qū)完全不充填時，采空區(qū)的頂板與底板板會發(fā)生剪切與拉伸破壞，造成頂板垮落破壞與底鼓。采空區(qū)頂板垂直位移隨充填率的增加呈線性減小。隨著充填率的不斷增加，充填體控制頂板的能力不斷增強。

參考文獻：

[1]李俊，寇云鵬，楊震，等.“三下”壓煤充填開采研究現(xiàn)狀綜述[J].黑龍江科技信息,2012（12）:43-44.

[2]余偉健，馮濤，王衛(wèi)軍，等.充填開采的協(xié)作支撐系統(tǒng)及其力學(xué)特征[J].巖石力學(xué)與工程學(xué)報,2012（S1）:2803-2813.

[3]繆協(xié)興.綜合機械化固體充填采煤技術(shù)研究進展[J].煤炭學(xué)報,2012（8）:1247-1255.

[4]王家臣，楊勝利，楊寶貴，等.長壁矸石充填開采上覆巖層移動特征模擬實驗[J].煤炭學(xué)報,2012（8）:1256-1262.

[5]胡炳南.我國煤礦充填開采技術(shù)及其發(fā)展趨勢[J].煤炭科學(xué)技術(shù),2012（11）:1-5，18.

[6]張吉雄，繆協(xié)興，郭廣禮.矸石(固體廢物)直接充填采煤技術(shù)發(fā)展現(xiàn)狀[J].采礦與安全工程學(xué)報,2009（4）:395-401.

[7]繆協(xié)興，張吉雄，郭廣禮.綜合機械化固體充填采煤方法與技術(shù)研究[J].煤炭學(xué)報,2010（1）:1-6.

Numerical Simulation on Failure Law of Overlying Strata by Paste Filling Parameters

ZHANG Cheng
(Xinzhou Geological Prospecting Team of Coal Field,Xinzhou 034000,China)

In the traditional belt mining,roadway digging disturbs the stress distribution of upper rock.Under the action of secondary stress field,the roof of roadways will be deformed and damaged.When the problem gets more serious,structural instability,collapse,and failure will occur.To study the impact of different paste filling parameters on the stress distribution of the overlying rock,numerical simulation was used to achieve the influence of different filling radio and paste strength on the failure law of belt-filling mining in Xiaotun mine,taking the plastic failure range and the subsidence of the roof and coal pillars in goaf as filling effect parameters.

paste filling parameters;strata;numerical simulation

TD823.7

A

1672-5050（2015）02-0028-04

10.3969/j.cnki.issn1672-5050sxmt.2015.02.010

（編輯：劉新光）

2015-01-10

張誠(1968-)，男，山西偏關(guān)人，大學(xué)本科，工程師，從事礦山地質(zhì)測量工作。