張忠厚，左 彪，黃厚旭

(1.遼寧工程技術(shù)大學(xué)理學(xué)院，遼寧阜新 123000;2.遼寧工程技術(shù)大學(xué)力學(xué)與工程學(xué)院，遼寧阜新 123000;3.遼寧工程技術(shù)大學(xué)土木與交通學(xué)院，遼寧阜新 123000)

0 引言

錢鳴高提出了采場上覆巖層活動的關(guān)鍵層理論并用彈性組合梁模型確定了當(dāng)較軟巖層與堅硬關(guān)鍵層協(xié)調(diào)變形時關(guān)鍵層上的等效荷載［1］，李文秀等用彈性薄板法建立了開采時地表下沉的預(yù)測力學(xué)模型［2］，劉開云等利用彈性薄板法分析了采空區(qū)上覆巖層的組合運動特征并對關(guān)鍵層的位置進行了計算分析［3］。利用彈性板理論研究頂板撓曲變形已經(jīng)取得了許多有價值的成果［3-6］。但在頂板撓度計算中已有的研究大多是先建立薄板力學(xué)模型然后求解平衡方程，本文將試圖從另一個角度求解離層撓度問題。最小勢能原理表明:在合外力作用下真實存在的位移場使平衡體系所具有的勢能最小。將覆巖頂板簡化為四邊簡支的彈性薄板然后建立平衡方程求解彈性薄板在上部荷載及自重下的撓度值是目前研究覆巖離層及地表沉陷使用較多的一種方法。在上部荷載及自重作用下四邊簡支的彈性頂板的撓曲方程符合雙三角級數(shù)。注漿減沉的力學(xué)機理是:漿體充填離層空間，固結(jié)后對關(guān)鍵層提供支撐反力進而減小地表沉降。本文將以最小勢能原理為基礎(chǔ)計算關(guān)鍵層在注漿前后的最大撓度值，并用求得的結(jié)果評價注漿減沉的效果。

1 離層巖板彈性力學(xué)模型建立

在開采形成采空區(qū)后位于采空區(qū)上部的巖石失去下部支撐，在上部荷載及自重作用下將發(fā)生向下的撓曲變形。由于成巖年代、巖石成分及上覆荷載等存在差異，在下沉過程中不同巖板將產(chǎn)生不同的撓度。關(guān)鍵層對巖體的活動起全部的控制作用受關(guān)鍵層控制的各巖層將與關(guān)鍵層協(xié)調(diào)變形且撓曲形式相同。在離層空間上方關(guān)鍵層巖板的懸露面積遠大于其厚度，而關(guān)鍵層的硬脆性決定了在未破壞范圍內(nèi)關(guān)鍵層將發(fā)生彈性變形且其撓度小于其厚度，因此將關(guān)鍵層巖板簡化為彈性板是合理的。且根據(jù)已有的研究成果可將關(guān)鍵層巖板簡化為四邊簡支的矩形巖板。在應(yīng)用彈性薄板理論分析關(guān)鍵層的離層變形之前先對關(guān)鍵層做如下假設(shè)［7］:(1)垂直于薄板中平面的法線在變形后仍為垂直彈性曲面的直線且長度不變;(2)垂直于中平面方向的正應(yīng)力與其他兩個方向的正應(yīng)力相比可以忽略，即σz=0;(3)板的中平面無伸縮變形。

由彈性薄板理論知，薄板在合外力作用下的平衡方程為:

或表示為:

其中:

其中amn為關(guān)鍵層在彈性限度內(nèi)的最大撓度，m、n分別為巖板在x、y方向撓曲的半波數(shù)，其中m、n=1、2…，a、b分別為巖板在 x、y方向的長度，q為關(guān)鍵層上的等效均布荷載，D為關(guān)鍵層的抗彎剛度，E為關(guān)鍵層巖板的彈性模量，μ為關(guān)鍵層巖板的泊松比，d為巖板厚度。

2 關(guān)鍵層上等效均布荷載q的計算

將關(guān)鍵層當(dāng)作覆巖撓曲變形的分界點，與關(guān)鍵層協(xié)調(diào)變形的上覆巖層與關(guān)鍵層形成一個組合梁，依據(jù)巖層的自然分層將采場上覆巖層從下至上分別編號1、2、3…m … k，k ＞m 對應(yīng)的厚度分別為 d1、d2、d3… dm… dk，密度分別為 ρ1、ρ2、ρ3…ρm…ρk，并假設(shè)第 m 層為關(guān)鍵層(圖1)。則第m層至第k層將發(fā)生相同撓曲變形并且各層之間緊密接觸這種變形協(xié)調(diào)關(guān)系可表示為:

圖1 覆巖離層力學(xué)模型Fig.1 Mechanical model of cover rock seperation stata

由材料力學(xué)知:

即:

將(8)代入式(6)得:

由式(9)可解得:

式中Mi(i=m…k)為第i層巖板承受的彎矩，Ii第i層巖層的截面慣性矩。由組合梁理論知:

式中Mhe代表組合梁承受的總彎矩。

將式(10)代入式(12)可得:

由于組合梁的彎矩可視為由組合梁上的均布荷載組成，則有:

即:

由qisi=ρisidig得:

式中si為相應(yīng)巖板的有效計算面積。

將式(16)代入式(15)可得:

3 關(guān)鍵層的最大撓度計算

3.1 注漿前關(guān)鍵層的最大撓度計算

由于離層的產(chǎn)生是一個長期而緩慢的過程，因此可將關(guān)鍵層視為在合外力作用下的平衡體系.設(shè)在上部荷載及自重作用下關(guān)鍵層產(chǎn)生的最大撓度為amn。

當(dāng)關(guān)鍵層發(fā)生撓曲ω后所具有的彈性勢能為:

但由于關(guān)鍵層處于四邊簡支狀態(tài)，因此上式中第二項將不存在，所以此處關(guān)鍵層的勢能表達式為:

其中:

將式(20)代入式(19)中，則有:

關(guān)鍵層上等效均布荷載在撓曲方向上做的功為:

關(guān)鍵層巖板所具有的總勢能為:

將式(21)和(22)代入(23)中，則有:

最小勢能原理表明:對于給定的應(yīng)力場平衡體系將發(fā)生讓系統(tǒng)總勢能最小的實際位移，則此時的位移amn應(yīng)滿足:

將式(24)代入式(25)中，則有:

由式(21)可知:

若僅取級數(shù)的第一項，即m=1，n=1則得關(guān)鍵層的最大撓度為:

3.2 注漿后關(guān)鍵層的最大撓度計算

在離層空間中充填的粉煤灰漿體在固結(jié)后由于其硬度和彈性模量遠小于堅硬的關(guān)鍵層頂板，因此可將固結(jié)后的漿體視為處于關(guān)鍵層下的彈性墊層。假設(shè)注漿時關(guān)鍵層只發(fā)生了大小可以忽略的撓曲下沉且注漿量足夠，關(guān)鍵層的下沉和凝固漿體的變形是同步且兩種介質(zhì)之間接觸良好，由Winkler地基假設(shè)知固結(jié)的漿體對關(guān)鍵層的支撐反力為［7］:

式中F是關(guān)鍵層受到的支撐反力，K是Winkler地基系數(shù)與基礎(chǔ)本身的結(jié)構(gòu)及力學(xué)特性有關(guān)，彈性地基系數(shù)K=(E0/h0)-2，其中E0，h0分別為墊層的彈性模量和厚度。ωz為注漿后關(guān)鍵層的撓度，由于支撐反力方向與ωz的方向相反，故F的值取負號。

注漿后關(guān)鍵層仍處緩慢下沉狀態(tài)，但由于變形速率慢，因此可以將撓曲變形過程看作平衡狀態(tài)。當(dāng)關(guān)鍵層發(fā)生最大撓度amnz時具有的彈性勢能為:

關(guān)鍵層的撓曲方程可表示為:

由式(21)知:

此時合外力做功為:

將式(20)和(22)代入式子(33)則有:

注漿后關(guān)鍵層的總勢能為:

將式(32)和(34)代入(35)中則有:

由平衡體系滿足最小勢能原理則有:

將式(36)代入式(37)中則有:

經(jīng)整理后得:

因此充填注漿后關(guān)鍵層的撓曲方程可以表示為:

在式(39)中僅取級數(shù)第一項，當(dāng)m=n=1得:

4 注漿減沉效果分析

4.1 薄板所受最大拉力方向的確定

由于注漿前后巖板的變形都是處于小撓曲變形范圍內(nèi)，因此可認為中平面依然保持為平面［7］。則巖板沿x、y方向的應(yīng)力分別為:

式中μ為巖板材料的泊松比且μ＜1，z為受力點距中平面的距離，在中平面上部時取“+”號，在中平面下部時取“-”號。由于最大拉應(yīng)力出現(xiàn)在關(guān)鍵層最底端，即分別代入式(36)和(37)則得到最大拉應(yīng)力σx、σy的表達式為:

在式(46)和(47)中當(dāng)改變 a、b的相對大小時得:

上式表明四邊簡支矩形巖板的最大拉應(yīng)力發(fā)生在撓度最大處且沿頂板的短邊方向。

4.2 關(guān)鍵層巖板的極限跨距的計算

在計算薄巖板撓度時應(yīng)保證計算是在巖板的極限跨距內(nèi)進行且滿足，但實際開采中大多數(shù)情況下離層巖板的a、b值都大于巖板的極限跨距，因此必須對巖板極限撓度計算公式予以修正。彈性力學(xué)中將受到外部荷載作用時薄巖板最大撓度處正好發(fā)生拉裂時對應(yīng)的長度稱為薄巖板的極限跨距，分別用am、bm表示。薄巖板的極限跨距與巖板自身的力學(xué)性質(zhì)、外觀尺寸以及作用于巖板上的荷載的分布情況等都有關(guān)系，由板的抗拉強度可推導(dǎo)計算am、bm的公式為:［8］

式中k為薄巖板的形狀系數(shù)，由am、bm的表達式

可知薄巖板的抗拉強度是影響薄巖板極限跨距的主要因素。

將am、bm代入式(28)和(40)中可得關(guān)鍵層注漿前后的最大撓度分別為:

將式(51)代入式(29)和(41)中可得關(guān)鍵層注漿前后的撓度表達式為:

工程實際中常用注漿前后地表最大沉降值的變化作為衡量地表減沉的主要指標(biāo)，考慮注漿是將漿液注入導(dǎo)水裂隙帶之上的離層中，因此離層的下位巖層必須完整沒有破壞，即該巖層不處于冒落帶或斷裂帶范圍之內(nèi)。假設(shè)地表土層對因巖板撓曲而產(chǎn)生的沉陷的弱化作用可以忽略，關(guān)鍵層的撓曲變形可以最大限度的傳遞到地表，地表的最大沉降值和關(guān)鍵層的最大撓度相等，那么注漿的減沉效果就可以表示為:

將式(51)代入式(53)中可得到離層注漿減沉系數(shù)的計算公式為:

4.3 實例計算

某礦43下02工作面是四采區(qū)開采的第一個工作面，開采3下煤層。該工作面位于北翼3下回風(fēng)大巷東320m。工作面走向長800m，斜長180m，煤層厚度1.98～5.30m，平均4m，含夾矸一層，最厚達1.2m。煤層傾角3°～15°，平均8°開采深度 556～588m，平均572m，各巖層力學(xué)性質(zhì)見表1。

表1 巖層特性參數(shù)Table 1 Rock characteristic parameters

由上表巖石參數(shù)，通過覆巖關(guān)鍵層位置的快速判別［9］可確定巖層10為關(guān)鍵層。因為離層注漿是將漿液注入導(dǎo)水裂隙帶之上的離層中，所以必須確保離層的下位巖層完整，即下位巖層不發(fā)生破壞，否則該巖層將位于冒落帶或斷裂帶范圍內(nèi)［10］。故關(guān)鍵層的極限跨距尺寸取決于離層下位巖層的極限跨距尺寸。故此實例中關(guān)鍵層的極限跨距尺寸取決于巖層9，由工作面采礦條件分析知H=572m，a/b=4。將巖層9的相關(guān)參數(shù)代入式(49)后可以得到am=232m，bm=58m由式(17)及表1中數(shù)據(jù)可算得D10=7286.42GPa·m3。

表2 粉煤灰在封閉空間中的力學(xué)參數(shù)［11］Table 2 Mechanical parameters of fly ash

將注漿材料粉煤灰作取樣實驗，結(jié)果表明其彈性模量 E處于 0.03～0.05GPa之間，此處取 E=0.04GPa，經(jīng)過計算上覆巖層相鄰兩巖板間的最大離層空間高度為h0max=1.2m，而實際注漿范圍可等效為長和寬分別為800m和200m，面積S=160000m2的矩形，經(jīng)計算離層空間的大小為 V=68947m3［11］，假設(shè)注采比合適且漿液剛好注滿離層空間，將固結(jié)的注漿體看作一個矩形的薄板鋪在關(guān)鍵層之下，則薄板厚度等于注漿體的平均厚度為h0=V/S，計算得到h0=0.43m。取h0=0.43m，E=0.04GPa計算后可得等效地基系數(shù)為 K=0.025GPa/m，將am=232m，bm=58m，D10=7286.42GPa·m3代入式(54)得離層注漿的減沉率為 η=26.1%，這與注漿減沉的實際效果相符［12］。

如圖2所示，當(dāng)最接近地表的關(guān)鍵層和此關(guān)鍵層對應(yīng)的極限跨距為am=232m，bm=58m時，離層注漿的減沉效果隨注漿體固結(jié)后形成的下臥支撐體的等效地基系數(shù)的增大將逐漸提高。且等效地基系數(shù)K從0增加到0.3時減沉系數(shù)η快速增長并超過0.8，因此，當(dāng)?shù)刃У鼗禂?shù)K處于該范圍時增大K能明顯提高減沉效果;當(dāng)?shù)刃У鼗禂?shù)K大于0.8時減沉系數(shù)增長緩慢，此范圍內(nèi)K的增大對減沉效果的影響不明顯;而下臥支撐體的等效地基系數(shù)K主要與灌注材料的性質(zhì)、級配、含水率有關(guān)，在其他條件不變的情況下材料的彈性模量越高，地基系數(shù)K越大，材料級配越均勻密實，地基系數(shù)K越大，含水率達到最佳時，地基系數(shù)K最大。然而實際減沉效果ηd往往達不到理想效果η，實際工程中，由于不同的填充條件和物理環(huán)境，之間為了提高加固效果，獲得更好的減沉效果，注漿前應(yīng)先分析離層動態(tài)運動規(guī)律，合理設(shè)計布置注漿鉆孔，計算填充材料的注入量，注漿時可以適當(dāng)提高注漿壓力(充填材料反力大，也就是充填密實)。漿液比重最好為1.1～1.2t/m3，比重太大，不易輸送灌注，且易堵塞注漿管，比重太小，材料沉縮率大，影響加固減沉效果。

圖2 等效地基系數(shù)K和減沉效果η的關(guān)系曲線Fig.2 The relationship of the equivalent foundation coefficient K and the reduce sink effect η

5 結(jié)論

(1)建立了開采后關(guān)鍵層的力學(xué)模型。將緩慢下沉并最終達到穩(wěn)定的關(guān)鍵層巖板看作平衡體系，依據(jù)平衡體系具有最小勢能的原理求出了關(guān)鍵層注漿前后的最大撓度。

(2)頂板發(fā)生撓曲變形時，在巖板的最大撓度處沿其短邊方向?qū)⒊霈F(xiàn)最大拉應(yīng)力。

(3)將注漿前后關(guān)鍵層最大撓度的變化量間接描述地表的沉陷值，從而建立注漿減沉效果的評價模型，得出了注漿減沉評價的一般表達式。

(4)注漿前應(yīng)先分析離層動態(tài)運動規(guī)律，合理設(shè)計布置注漿鉆孔，計算填充材料的注入量，注漿時可以適當(dāng)提高注漿壓力(充填材料反力大，也就是充填密實)。

［1］錢鳴高，繆協(xié)興，許家林.巖層控制中的關(guān)鍵層理論研究［J］.煤炭學(xué)報，1996，21(3):225-230.QIAN Minggao，MIU Xiexing，XU Jialin.Theoretical study of key stratum in ground control［J］.Journal of China Coal Society，1996，21(3):225-230.

［2］李文秀，梁旭黎，趙勝濤，等.地下開采引起地表沉陷預(yù)測的彈性模版法［J］.工程力學(xué)，2006，23(8):177-181.LI Wenxiu，LIANG Xuli，ZHAO Shengtao，et al.The elastic-thin-plate bending method for predicting ground subsidence due to undergroundmining［J］.Engineering Mechanics，2006，23(8):177-181.

［3］劉開云，喬春生，周輝.覆巖組合運動特征及關(guān)鍵層位置研究［J］.巖石力學(xué)與工程學(xué)報，2004，23(8):1301-1306.LIU Kaiyun，QIAO Chunsheng，ZHOU Hui.Research on combined motion characteristics of overlying rock stratum and position of key strata［J］.Chinese Journal of Rock Mechanics and Engineering，2004，23(8):1301-1306.

［4］郝延錦，吳立新，戴華陽.用彈性板理論建立地表沉陷預(yù)計模型［J］.巖石力學(xué)與工程學(xué)報，2006，25(增1):2958—2962.HAO Yanjin，WU Lixin，DAI Huayang.Establing a ground settlement prediction model with elastic slab theory［J］.Chinese Journal of Rock Mechanics and Engineering，2006，25(S1):2958-2962.

［5］王紅衛(wèi)，陳忠輝，杜澤超，等.彈性薄板理論在地下采場頂板變化規(guī)律研究中的應(yīng)用［J］.巖石力學(xué)與工程學(xué)報，2006，25(S2):3769—3774.WANG Hongwei，CHEN Zhonghui，DU Zechao，et al.Application of elastic thin plate theory to changerule of roof in underground stope［J］.Chinese Journal of Rock Mechanics and Engineering，2006，25(S2):3769-3774.

［6］翟所業(yè)，張開智.用彈性板理論分析采場覆巖中的關(guān)鍵層［J］.巖石力學(xué)與工程學(xué)報，2004，23(11):1856-1860.ZHAI Suoye，ZHANG Kaizhi.Analysis on key layer of overlying strata with elastic slab theory［J］.Chinese Journal of Rock Mechanics and Engineering，2004，23(11):1856-1860.

［7］曲慶璋，章權(quán)，季求和，等.彈性板理論［M］.北京:人民交通出版社，1999.QU Qingzhang，ZHANG Quan，JI Qiuhe，et al.Theory of elastic plates［M］.Beijin:China Communications Press，1999.

［8］李小強，白世偉，李鈾.最小勢能原理在二維邊坡穩(wěn)定分析中的應(yīng)用［J］.巖土力學(xué)，2004，25(6):909-912.LI Xiaoqiang，BAI Shiwei，LI You.2D slope stability analysis using principle of minimum potential energy［J］.Rock and Soil mechanics，2004，25(6):909-912.

［9］楊本生，張磊.覆巖關(guān)鍵層位置的快速判別及來壓預(yù)測［J］.煤炭工程，2013(7):79-81.YANG Bensheng，ZHANG Lei.The distinguish and forecast of the place of the key cover stratum［T］Coal Engineering，2013(7):79-81.

［10］蘇仲杰，于廣明，楊倫.覆巖離層變形力學(xué)模型及應(yīng)用［J］.巖土工程學(xué)報，2002，24(6):778-781 SU Zhongjie，YU Guangming，YANG Lun.Application of mechanicalmodel to deformation of covered rock separation strata［J］.Chinese Journal of Geotechnical Engineering，2002，24(6):778-781.

［11］蘇仲杰.采動覆巖離層變形機理研究［D］，阜新:遼寧工程技術(shù)大學(xué)，2001.SU Zhongjie.Research of deformation mechanism of mining overburden separated strata［D］.Fuxin:Liaoning Technical University，2001.

［12］Palaski J.The experiment and practical results of applying backfill［C］∥Innovations in Mining Backfill Technology，Proceedings of 4th Intemational Symposium 011 Mining with Backfill，Montreal:［s.n］，1989:33-37.