李春睿

(1.天地科技股份有限公司開采設(shè)計(jì)事業(yè)部,北京市朝陽(yáng)區(qū),100013;2.煤炭科學(xué)研究總院巖石力學(xué)與支護(hù)專業(yè)委員會(huì),北京市朝陽(yáng)區(qū),100013)

確定支架合理支護(hù)強(qiáng)度的數(shù)值模擬方法

李春睿1,2

(1.天地科技股份有限公司開采設(shè)計(jì)事業(yè)部,北京市朝陽(yáng)區(qū),100013;2.煤炭科學(xué)研究總院巖石力學(xué)與支護(hù)專業(yè)委員會(huì),北京市朝陽(yáng)區(qū),100013)

以支架和圍巖的相互作用關(guān)系為理論依據(jù),利用FLAC3D數(shù)值模擬軟件對(duì)液壓支架合理支護(hù)強(qiáng)度的確定方法進(jìn)行了研究,提出了兩種方法來(lái)確定液壓支架的支護(hù)強(qiáng)度,一種是以工作面控頂區(qū)巖(煤)體的下沉量來(lái)判斷支架的合理支護(hù)強(qiáng)度,另一種是以頂板(煤)對(duì)支架的直接作用力來(lái)判斷合理的支護(hù)強(qiáng)度,并結(jié)合現(xiàn)場(chǎng)應(yīng)用實(shí)例說(shuō)明了兩種方法的合理性。

液壓支架 支護(hù)強(qiáng)度 數(shù)值模擬 控頂區(qū)域

目前,對(duì)于確定液壓支架合理工作阻力的方法主要有基于巖重估算的倍重系數(shù)法、基于大量實(shí)測(cè)結(jié)論的統(tǒng)計(jì)類比法、基于頂煤體損傷力學(xué)的解析模型法和基于頂板分類的巖層結(jié)構(gòu)法等。近年來(lái),隨著計(jì)算機(jī)硬件及軟件技術(shù)的發(fā)展,基于計(jì)算機(jī)仿真技術(shù)的數(shù)值模擬法在解決采礦工程的問(wèn)題中得到了廣泛的應(yīng)用,然而在確定液壓支架的合理支護(hù)強(qiáng)度時(shí),數(shù)值模擬過(guò)程中參數(shù)的選擇、模擬方案的確定與結(jié)論密切相關(guān)。

2 數(shù)值模擬方法確定工作阻力

2.1 利用控頂區(qū)巖體的下沉量來(lái)判斷

圖1 理想結(jié)構(gòu)模型

利用控頂區(qū)煤(巖)體的下沉量來(lái)判斷合理的支護(hù)強(qiáng)度,即控頂區(qū)下沉法。理想結(jié)構(gòu)模型見圖1。在煤層開采時(shí),分別對(duì)工作面控頂區(qū)域上方的煤(巖)體施加不同支護(hù)強(qiáng)度載荷(一般從0～2.0 MPa),并監(jiān)測(cè)不同載荷時(shí)A、B、C、D、E 5個(gè)點(diǎn)處頂板(煤)的位移情況。為了避免連續(xù)介質(zhì)方法造成的邊界效應(yīng)影響,刪除掉采空區(qū)后方已破壞的煤巖體區(qū)域。

圖2 數(shù)值模擬模型

為了更詳盡地說(shuō)明該方法,以某礦工作面煤巖體的物理力學(xué)參數(shù)為依據(jù)建立模型,某一時(shí)態(tài)圍巖的垂直應(yīng)力云圖見圖2,考慮到支架的控頂區(qū)域在5 m左右,因此分別監(jiān)測(cè)了距離工作面煤壁1 m、2 m、3 m、4 m、5 m處頂板的A、B、C、D、E 5個(gè)點(diǎn),圖2中的箭頭為某時(shí)刻頂板的下沉位移矢量。

結(jié)合該礦區(qū)已有回采面的支護(hù)強(qiáng)度范圍,此次模擬的支護(hù)強(qiáng)度q分別取0 MPa、0.2 MPa、0.4 MPa、0.6 MPa、0.8 MPa、1.0 MPa、1.2 MPa、1.4 MPa和1.6 MPa,獲得前述5個(gè)點(diǎn)處的頂板下沉量數(shù)據(jù),不同支護(hù)強(qiáng)度作用下的頂板下沉量見表1,并分析了頂板下沉量與支架支護(hù)強(qiáng)度之間的關(guān)系曲線,見圖3。

圖3 控頂區(qū)域的頂板下沉量與支架支護(hù)強(qiáng)度關(guān)系曲線

表1 不同支護(hù)強(qiáng)度作用下的頂板下沉量

數(shù)值模擬表明,隨著支護(hù)強(qiáng)度的增大,頂板下沉位移量逐漸減小,當(dāng)支護(hù)強(qiáng)度大于1.2 MPa后,E點(diǎn)處巖體甚至向上移動(dòng)。進(jìn)一步分析可以得知,當(dāng)支架支護(hù)強(qiáng)度大于0.6 MPa時(shí),頂板位移量隨支架強(qiáng)度增加已不再發(fā)生明顯的變化,即當(dāng)支架支護(hù)強(qiáng)度大于0.6 MPa時(shí),能夠很好地控制該頂板,因此該工作面最后確定的合理支護(hù)強(qiáng)度為0.6 MPa。

2.2 利用支架頂梁的受力情況來(lái)判斷

利用支架頂梁的受力情況來(lái)判斷合理的支護(hù)強(qiáng)度,即支架頂梁受力法。需要建立液壓支架的數(shù)值模型,某四柱式液壓支架模型見圖4。在建立該數(shù)值模型時(shí),煤巖層采用應(yīng)變軟化模型(Strain-Hardenging/Softening Model),而支架部分采用了各向同性的彈性模型(Elastic Isotropic Model),這樣可以避免計(jì)算過(guò)程中支架被壓壞,計(jì)算時(shí)分別對(duì)液壓支架頂梁處A、B、C、D、E 5個(gè)點(diǎn)的受力情況進(jìn)行監(jiān)測(cè),根據(jù)監(jiān)測(cè)得到的頂梁受力變化情況來(lái)判斷液壓支架的合理工作阻力。

同樣以某礦工作面煤巖體的物理力學(xué)參數(shù)為依據(jù)建立模型,某一時(shí)態(tài)圍巖的垂直應(yīng)力云圖見圖5,圖5中的曲線為計(jì)算過(guò)程中支架頂梁的受力變化曲線,由于5個(gè)點(diǎn)處的曲線重合,因此可以判斷支架頂梁承受著均布的荷載。從頂梁的應(yīng)力變化曲線中可見,隨著計(jì)算時(shí)步的增加,該工作面液壓支架的頂梁受力狀態(tài)是逐漸增加的過(guò)程,直至增加到0.78 MPa時(shí)保持平穩(wěn)狀態(tài),據(jù)此,液壓支架選型時(shí)支護(hù)強(qiáng)度取0.8 MPa左右的結(jié)論。

圖4 某四柱式液壓支架理想結(jié)構(gòu)模型

圖5 某一時(shí)態(tài)圍巖的垂直應(yīng)力數(shù)值模擬模型

根據(jù)配套尺寸、支架頂梁長(zhǎng)度、空頂距算出支架工作阻力:

式中:P——支架工作阻力,kN;

q——支架的支護(hù)強(qiáng)度,取0.8 MPa;

LK——空頂距,取0.36 m;

LD——頂梁長(zhǎng)度,取4.53 m;

B——支架寬度,取1.5 m。

根據(jù)上面分析結(jié)果,P=5868 kN。作為4柱液壓支架,每個(gè)立柱額定阻力為1467 kN。因此安全閥開啟壓力為35.32 MPa。

圖6 運(yùn)算10000時(shí)步液壓支架立柱受力情況

為了進(jìn)一步分析上述由頂梁間接得到的支護(hù)阻力是否準(zhǔn)確,又監(jiān)測(cè)了液壓支架立柱上所受到的壓力變化情況,運(yùn)算10000時(shí)步液壓支架立柱受力情況見圖6,從圖6曲線可見,立柱所受壓力由0 MPa逐漸增加,當(dāng)增加至36 MPa時(shí),立柱所受力已達(dá)最大值,此后繼續(xù)運(yùn)算發(fā)現(xiàn),由于上覆巖體斷裂的影響,液壓支架立柱所承受的壓力有所波動(dòng),但最大值仍然沒有超過(guò)36 MPa,說(shuō)明安全閥開啟壓力為35.32 MPa基本能夠滿足要求。

3 現(xiàn)場(chǎng)應(yīng)用實(shí)例與驗(yàn)證

3.1 控頂區(qū)下沉法的應(yīng)用與驗(yàn)證

水簾洞煤礦煤層平均埋深370 m,煤層平均厚度為12.2 m,煤層傾角為5～6°,煤層的硬度系數(shù)f為3.02;直接頂厚度為3.51 m,f為4.03;老頂厚度為7.9 m,f為5.6。根據(jù)煤巖層物理力學(xué)條件,建立數(shù)值模型并應(yīng)用控頂區(qū)下沉法進(jìn)行了數(shù)值模擬,得到綜放支架支護(hù)強(qiáng)度為0.9 MPa。根據(jù)支護(hù)強(qiáng)度(q)、支架控頂長(zhǎng)度(LK)、端面距(LD)和支架寬度(B)計(jì)算得出支架工作阻力為6377 kN,選用ZF6400/17/32型低位放頂煤液壓支架。

在現(xiàn)場(chǎng)對(duì)水簾洞礦ZF2801工作面進(jìn)行了礦壓觀測(cè),觀測(cè)了初期來(lái)壓和6個(gè)周期來(lái)壓。初次來(lái)壓時(shí)最大工作阻力為6586.8 kN,時(shí)間加權(quán)平均工作阻力為3927.6 kN;周期來(lái)壓時(shí)最大工作阻力為5940.4 kN,時(shí)間加權(quán)平均工作阻力為3424.6 kN;從監(jiān)測(cè)結(jié)果及現(xiàn)場(chǎng)應(yīng)用情況可以看出支架工作阻力選擇得比較合適,能夠很好地適應(yīng)工作面的正常生產(chǎn)。

3.2 支架頂梁受力法的應(yīng)用與驗(yàn)證

照金煤礦煤層平均埋深350 m,煤層平均厚度為9 m,采煤機(jī)割煤高度2.6 m,放煤高度6.4 m,煤的硬度系數(shù)f為2.5。直接頂一般厚3 m,粉細(xì)砂巖,局部為泥巖,屬不堅(jiān)實(shí)至中等堅(jiān)實(shí)頂板?；卷敒?.03～19.81 m的粗砂巖,較為堅(jiān)硬。根據(jù)煤巖層物理力學(xué)條件,建立數(shù)值模型并應(yīng)用支架頂梁受力法進(jìn)行了數(shù)值模擬,得到綜放支架支護(hù)強(qiáng)度為0.8 MPa左右。該煤層選擇5868 kN的四柱式放頂煤液壓支架,安全閥開啟壓力35.32 MPa,能夠滿足頂板要求。

工業(yè)試驗(yàn)階段,對(duì)照金礦101綜放工作面進(jìn)行了礦壓觀測(cè),觀測(cè)初期來(lái)壓和3個(gè)周期來(lái)壓。最大工作阻力為5142 kN(0.781 MPa),時(shí)間加權(quán)平均工作阻力為3597 kN(0.735 MPa)。通過(guò)模擬可知,與數(shù)值模擬結(jié)論的比值分別為87.62%、61.29%。根據(jù)現(xiàn)場(chǎng)反饋的情況,在該工作面開始使用的ZF5200/18/28型號(hào)支架的工作阻力略有偏小,來(lái)壓期間安全閥開啟頻繁,可見,此次選擇的5868 kN的液壓支架較為合理,能夠滿足安全生產(chǎn)要求。

4 結(jié)語(yǔ)

控頂區(qū)下沉法是主動(dòng)向頂板(煤)施加作用力,觀測(cè)頂板(煤)的下沉量,當(dāng)所施加的作用力能夠使頂板(煤)控制在一定范圍內(nèi)時(shí),則所施加的力為合理的支護(hù)強(qiáng)度;支架頂梁受力法是建立了支架模型,由支架被動(dòng)地承受頂板(煤)的作用力,當(dāng)支架的頂梁所承受的力達(dá)到平穩(wěn)時(shí),此時(shí)的壓力值視為合理的支護(hù)強(qiáng)度。

如果側(cè)重于研究支架初撐力的變化對(duì)頂板(煤)的破壞作用效果、分析頂板(煤)的下沉量,選擇控頂區(qū)下沉法會(huì)更好一些;如果側(cè)重于研究支架立柱的受力狀態(tài)、支架與圍巖的整體性受力關(guān)系、研究架前空頂距與片幫的關(guān)系,可選擇支架頂梁受力法。

[1] 史元偉.液壓支架初撐力及工作阻力的確定[J].煤炭學(xué)報(bào),1985(3)

[2] 史元偉.綜放工作面圍巖動(dòng)態(tài)及液壓支架載荷力學(xué)模型[J].煤炭學(xué)報(bào),1997(3)

[3] 于海勇,吳鳳東,孟金鎖等.放頂煤工作面支架受力特征及分析[J].煤炭學(xué)報(bào),1995(5)

[4] 范韶剛,閆少宏,毛德兵,富強(qiáng).確定綜放支架支護(hù)強(qiáng)度的統(tǒng)計(jì)類比法[J].礦山壓力與頂板管理,2000(4)

[5] 閆少宏.放頂煤開采支架工作阻力的確定[J].煤炭學(xué)報(bào),1997(1)

[6] 閆少宏.放頂煤開采頂煤與頂板活動(dòng)規(guī)律研究[D].北京:中國(guó)礦業(yè)大學(xué)(北京),1995

[7] 康立軍.長(zhǎng)壁綜放開采支架與頂煤相互作用關(guān)系研究[D].北京:煤炭科學(xué)研究總院,2000

[8] Ian Farmer.Coal Mine Structure[M].Chapman and Hall.Ltd.,1985

[9] 毛德兵,王延峰.數(shù)值模擬方法確定綜放工作面支架阻力[J].煤礦開采,2005(2)

[10] 王永秀.綜放開采中頂煤對(duì)支架合理支護(hù)強(qiáng)度影響的數(shù)值分析[J].煤炭學(xué)報(bào),2006(2)

Numerical simulation to determine rational supporting strength of roof supports

Li Chunrui1,2
(1.Coal Mining&Designing Branch,Tiandi Coal Co.,Ltd,Chaoyang,Beijing 100013,China;2.Rock Mechanics and Support Professional Committee,Coal Research Institute Chaoyang,Beijing 100013,China)

Research is conducted by taking the interaction between roof supports and the wall rock as the theoretical basis and using the FLAC3D to determine the rational supporting strength of hydraulic roof supports,which results in two methods-one using the level of subsidence of rock(or coal)mass in the roof control zone at the coalface and the other using the direct force of roof on the supports for determining the rational supporting strength.The rationality of the two methods is demonstrated in field application.

hydraulic roof support,supporting strength,numerical simulation,roof control zone

TD355.41

A

＊“十一五”國(guó)家科技支撐計(jì)劃重點(diǎn)項(xiàng)目課題(2008BAB36B01)

李春睿(1980-),男,漢,吉林九臺(tái)人,博士。2003年畢業(yè)于遼寧工程技術(shù)大學(xué),現(xiàn)就職于天地公司開采設(shè)計(jì)事業(yè)部工作。

(責(zé)任編輯 張艷華)