丁 楊

(中煤科工集團(tuán) 北京華宇工程有限公司，北京 100120)

礦山壓力和圍巖活動(dòng)密切相關(guān)，因此，探究礦山壓力的顯現(xiàn)規(guī)律必須和采場(chǎng)圍巖活動(dòng)規(guī)律聯(lián)系起來(lái)，通過(guò)現(xiàn)場(chǎng)礦壓觀測(cè)和對(duì)采場(chǎng)頂板結(jié)構(gòu)的假定，尋求不同條件下礦山壓力的現(xiàn)象規(guī)律。前人對(duì)采場(chǎng)圍巖結(jié)構(gòu)提出了許多假說(shuō)和理論，諸如壓力拱假說(shuō)、懸臂梁假說(shuō)、預(yù)成裂隙假說(shuō)、鉸接巖塊假說(shuō)，砌體梁理論[1]、關(guān)鍵層理論[2]和采場(chǎng)薄板理論。其中，梁理論忽略了頂板的空間影響[3]，而整體薄板理論忽略了結(jié)構(gòu)面的影響。一些學(xué)者提出鉸接沿板組假說(shuō)[3-4]，此假說(shuō)考慮了結(jié)構(gòu)面的影響，但巖板斷裂的隨機(jī)性無(wú)法預(yù)測(cè)。對(duì)于采場(chǎng)頂板來(lái)說(shuō)，除結(jié)構(gòu)面之外，巖體層間作用對(duì)頂板的斷裂位置及應(yīng)力分布也有較大的影響，本文從薄板理論和巖體層間作用的角度分析頂板斷裂規(guī)律，細(xì)化了頂板破斷的過(guò)程，為進(jìn)一步精確分析和預(yù)測(cè)頂板來(lái)壓提供了新的思路。

1 力學(xué)模型的分析

巖體的力學(xué)性質(zhì)，是巖塊和結(jié)構(gòu)面力學(xué)性質(zhì)的綜合反映，節(jié)理發(fā)育的巖體抗拉強(qiáng)度極弱，實(shí)際可認(rèn)為不可抗拉。當(dāng)采場(chǎng)頂板懸露至一定面積后，結(jié)構(gòu)面很大程度上影響了頂板的斷裂形式。對(duì)于層狀頂板巖層，當(dāng)處于懸伸狀態(tài)時(shí)，將會(huì)發(fā)生彎曲變形甚至斷裂回轉(zhuǎn)[5]。圖1為頂板斷裂結(jié)構(gòu)示意圖，頂板斷裂遵循“X-O”破斷形態(tài)。

圖1 頂板斷裂示意

Krichhoff在將材料視為均勻、連續(xù)、各向同性和線彈性的基礎(chǔ)上建立了薄板小撓度彎曲理論[6]，即：薄板發(fā)生彎曲變形后，垂直于薄板中面的法線長(zhǎng)度不變且仍垂直于薄板曲面，這忽略了剪應(yīng)變，薄板彎曲時(shí)，各點(diǎn)撓度與中面撓度相同。

頂板破斷前邊界條件如圖2所示。

圖2 頂板破斷前邊界條件

初次斷裂前，頂板嵌固在煤層中，可視為薄板四端固定的情況，根據(jù)采場(chǎng)薄板礦壓理論，其撓度及應(yīng)力表達(dá)式為[7]：

式中，ω為撓度；a，b為矩形板長(zhǎng)短邊長(zhǎng)；D為薄板抗彎強(qiáng)度；q為頂板均布載荷；E為彈性模量；μ為泊松比；h為頂板厚度，m；A為層間垂直作用力分布線傾角。

在初步來(lái)壓前采場(chǎng)長(zhǎng)邊與短邊中心先破壞，其受力變?yōu)樗亩撕?jiǎn)支，其撓度與應(yīng)力表達(dá)式為：

采場(chǎng)中部頂板斷裂之后，頂板兩端的受力變?yōu)槿斯讨В欢俗杂?，其撓度與應(yīng)力表達(dá)式為：

初始條件下，采場(chǎng)基本頂遵循“X-O”型破斷，但在采場(chǎng)中部巖體斷裂之后，兩端頂板則受層間作用影響，部分處于懸伸狀態(tài)，靠近端部的頂板與上層頂板粘結(jié)，如圖3所示。

圖3 巖層層間作用模型

當(dāng)采場(chǎng)中部巖層斷裂之后，分析采場(chǎng)上部與下部頂板，考慮巖層層間作用，其作用力主要包括巖層自重應(yīng)力、層間垂直作用力、層間摩擦作用力和層間粘結(jié)滑動(dòng)阻力[8]。利用彈性力學(xué)平面理論可解得其應(yīng)力分布。

自重應(yīng)力下的應(yīng)力分布：

式中，h為頂板厚度，m；γ為頂板巖層容重，MN/m3。

層間垂直作用力應(yīng)力分布:

σx=0,σy=0 (x

式中，a為頂板懸露長(zhǎng)度，m。

層間摩擦力應(yīng)力分布：

σx=0,σy=0 (x

式中，B為層間摩擦角。

層間粘結(jié)滑動(dòng)阻力應(yīng)力分布：

σx=0,σy=0 (x

式中，C為層間黏結(jié)力，MPa。

由以上分析可知，頂板懸伸狀態(tài)下只受自重應(yīng)力影響；而粘結(jié)狀態(tài)下，受到自重應(yīng)力、垂直應(yīng)力、摩擦應(yīng)力和粘結(jié)滑動(dòng)阻力的合力影響，所以在頂板懸伸范圍內(nèi)和在粘結(jié)范圍內(nèi)頂板垂直應(yīng)力分布為：

可見(jiàn)，在頂板懸伸范圍內(nèi)巖體垂直應(yīng)力小于粘結(jié)狀態(tài)下的頂板垂直應(yīng)力。

巖體抗拉強(qiáng)度遠(yuǎn)小于其抗壓或抗剪強(qiáng)度，采用最大拉應(yīng)力準(zhǔn)則判斷頂板斷裂位置。當(dāng)采場(chǎng)基本頂與上覆巖層處于粘結(jié)狀態(tài)時(shí)，頂板拉應(yīng)力為：

通過(guò)對(duì)基本頂?shù)牧W(xué)分析可知，基本頂遵循“X-O”型破斷，在中部頂板破斷之后，兩端頂板滯后斷裂，在此過(guò)程中，兩端頂板狀態(tài)如圖3所示，部分處于懸伸狀態(tài)，其余部分處于粘結(jié)狀態(tài)。因懸伸狀態(tài)下基本頂只受自重應(yīng)力影響，而粘結(jié)狀態(tài)下的頂板受層間作用影響，所以?xún)煞N狀態(tài)下基本頂?shù)拇怪睉?yīng)力不同。

2 工程背景

三溝鑫都煤業(yè)2101工作面長(zhǎng)度150m，埋深200m左右，煤層厚度為1.8m。直接頂板為砂質(zhì)泥巖，局部相變?yōu)榉凵皫r，厚度2.0m左右，基本頂為深灰色粉砂巖及砂巖，厚度3.0m左右，屬較堅(jiān)硬巖石，一般不易冒落。安裝液壓支架103架，在3號(hào)架、4號(hào)架、30號(hào)架、31號(hào)架、32號(hào)架、60號(hào)架、61號(hào)架、62號(hào)架、90號(hào)架、91號(hào)架、92號(hào)架、100號(hào)架、101號(hào)架分別安裝自動(dòng)記錄儀，并按順序編號(hào)。

3 力學(xué)模型的計(jì)算分析

在頂板斷裂的數(shù)值計(jì)算中，用ANSYS有限元數(shù)值模擬軟件對(duì)其進(jìn)行計(jì)算，以三溝鑫都煤業(yè)2101工作面地質(zhì)情況為背景，工作面長(zhǎng)度150m，基本頂厚度為3m，初始條件下基本頂受均布?jí)毫Γ卷敒樯罨疑凵皫r，彈性模量為25GPa，泊松比為0.28，巖體抗拉強(qiáng)度為5MPa，采用最大拉應(yīng)力破斷準(zhǔn)則來(lái)判斷頂板破壞情況[9]。當(dāng)工作面推進(jìn)50m時(shí)，頂板拉應(yīng)力如圖4所示。

圖4 基本頂拉應(yīng)力

分析以上計(jì)算結(jié)果可知，最大拉應(yīng)力出現(xiàn)在薄板長(zhǎng)邊的中部，且已超出巖體極限抗拉強(qiáng)度，達(dá)到7MPa，所以，頂板長(zhǎng)邊的中部首先破裂[10]，之后頂板短邊裂隙擴(kuò)展，頂板的應(yīng)力狀態(tài)發(fā)生改變。隨著工作面的推進(jìn)，由四邊簡(jiǎn)支應(yīng)力狀態(tài)下薄板的破斷方式可知，頂板中心裂隙沿x方向擴(kuò)展，工作面中部頂板先斷裂，產(chǎn)生初次來(lái)壓。

ANSYS中單元“殺死”并不是將單元?jiǎng)h除，而是將其剛度矩陣乘以一個(gè)很小的因子。死單元的單元載荷變?yōu)榱鉡11]。將中部斷裂薄板單元?dú)⑺篮?，采?chǎng)上部與下部薄板拉應(yīng)力如圖5所示。

圖5 單元?dú)⑺篮箜敯謇瓚?yīng)力

分析以上計(jì)算結(jié)果可以看出，工作面中部頂板斷裂之后，上部和下部頂板應(yīng)力狀態(tài)發(fā)生改變。將工作面中部頂板單元?dú)⑺乐?，由圖6可看出，工作面端部處頂板并非是在煤壁處達(dá)到抗拉強(qiáng)度，而且在頂板懸伸范圍內(nèi)垂直應(yīng)力小于端部處，且板自由邊中部區(qū)域出現(xiàn)壓應(yīng)力。

圖6 單元?dú)⑺篮箜敯宕怪睉?yīng)力

通過(guò)數(shù)值模擬計(jì)算可以看出，初始狀態(tài)下，工作面推進(jìn)至一定距離后，頂板的四邊最先發(fā)生破裂。在中部頂板斷裂后，即將中部單元?dú)⑺乐螅卷敳⒎鞘窃诿罕谔庍_(dá)到抗拉強(qiáng)度，而且懸伸范圍內(nèi)的頂板的垂直應(yīng)力小于端部處，模擬結(jié)果與力學(xué)分析結(jié)果相吻合。

4 現(xiàn)場(chǎng)實(shí)測(cè)數(shù)據(jù)

工作面方向上初次來(lái)壓垮落步距的變化曲線如圖7所示，可見(jiàn)，工作面中部頂板先來(lái)壓，上部與下部來(lái)壓滯后于中部。

圖7 工作面初次來(lái)壓步距曲線

由圖7可以看出，工作面中部初次來(lái)壓步距為26m左右，而工作面上部和下部初次來(lái)壓步距為46m左右。結(jié)合力學(xué)模型分析，工作面中部先斷裂，其邊界條件改變，故而工作面中部先來(lái)壓，之后采場(chǎng)兩端頂板斷裂，隨之來(lái)壓。從頂板初次斷裂的步距分析，頂板斷裂形態(tài)基本遵循“X-O”型。

由初次來(lái)壓強(qiáng)度曲線圖(圖8)看出，工作面中部來(lái)壓強(qiáng)度明顯高于中部與下部，而工作面端部并沒(méi)有明顯來(lái)壓。來(lái)壓強(qiáng)度與頂板的巖性以及厚度有關(guān)[12]，除此之外，頂板斷裂時(shí)的力學(xué)狀態(tài)也影響著來(lái)壓的強(qiáng)度。通過(guò)在層間作用下頂板垂直應(yīng)力的分析，可知懸伸范圍內(nèi)巖體垂直應(yīng)力小于與上層頂板粘結(jié)狀態(tài)下的頂板垂直應(yīng)力，因此，在工作面方向初次來(lái)壓強(qiáng)度出現(xiàn)波動(dòng)。

圖8 工作面初次來(lái)壓強(qiáng)度曲線

在層間作用影響下，其拉應(yīng)力為：

圖9 端部來(lái)壓強(qiáng)度

5 結(jié) 論

(1)層間作用對(duì)采場(chǎng)基本頂?shù)臄嗔岩?guī)律有著顯著影響，通過(guò)對(duì)基本頂在工作面方向上不同位置的垂直應(yīng)力分析，揭示了沿工作面方向上來(lái)壓強(qiáng)度存在波動(dòng)的機(jī)理。

(2)通過(guò)對(duì)層間作用影響下采場(chǎng)端部頂板的拉應(yīng)力分析，得出在特定參數(shù)下判斷頂板斷裂位置的方法，即x=a+4.4m，其中a為頂板的懸伸長(zhǎng)度。

(3)因頂板端部斷裂位置不同，所以端部頂板的懸伸與粘結(jié)長(zhǎng)度不同，導(dǎo)致采場(chǎng)端部頂板在工作面推進(jìn)方向上來(lái)壓強(qiáng)度存在波動(dòng)。