第33卷第8期

　2008年8月煤　　炭　　學(xué)　　報(bào)JOURNALOFCHINACOALSOCIETYVol.33　No.8　Aug.　2008　　　文章編號(hào):0253-9993(2008)08-0850-06

采空區(qū)堅(jiān)硬頂板破斷機(jī)理與災(zāi)變塌陷研究

王金安,尚新春,劉　紅,侯志鷹

司,山西大同　037003)1123(1.北京科技大學(xué)土木與環(huán)境工程學(xué)院,北京　100083;2.煤炭科學(xué)研究總院重慶研究院,重慶　400037;3.大同煤業(yè)集團(tuán)有限責(zé)任公

摘　要:以堅(jiān)硬頂板下的房柱式和條帶式采礦工程為背景,建立了表征采空區(qū)內(nèi)礦柱支撐頂板的

彈性基礎(chǔ)板力學(xué)模型,研究頂板不同階段的破斷模式與突變失穩(wěn)的力學(xué)過(guò)程.研究表明:當(dāng)煤柱

的有效承載面積逐漸減小到臨界值時(shí),非線性控制參數(shù)即可穿越分岔點(diǎn)集,頂板位移突跳產(chǎn)生極

限點(diǎn)失穩(wěn),煤柱-頂板系統(tǒng)出現(xiàn)突然塌陷失穩(wěn).

關(guān)鍵詞:堅(jiān)硬頂板;破斷;災(zāi)變;塌陷

中圖分類(lèi)號(hào):TD823.85　　　文獻(xiàn)標(biāo)識(shí)碼:A

Studyonfracturemechanismandcatastrophiccollapseofstrong

roofstrataabovetheminedarea

WANGJin-an,SHANGXin-chun,LIUHong,HOUZhi-ying

(1.SchoolofCivilandEnvironmentEngineering,UniversityofScienceandTechnologyBeijing,Beijing　100083,China;2.ChongqingResearchInsti-tute,ChinaCoalResearchInstitute,Chongqing　400037,China;3.DatongCoalMineGroupCo.Ltd.,Datong　037003,China)1123

Abstract:Basedonthepillar-roofandstrippillarminingengineeringbeneaththestrongroofstrata,themechani-calmodelofaplateseatedontheelasticfoundationwasestablishedtoanalyzethefailureprocessoftheroofstrataanditscatastrophiccollapse.Thestudyshowsthatthenon-linearcontrolparametercanpassthroughthepoint-setofbifurcationwhentheeffectiveloadbearingareaofpillarsdecreasestoitscriticalvalue,resultingintheroofstratatoproduceasuddenjumpofsubsidence,thepillar-roofstratasystembecomescatastrophiccollapse.

Keywords:strongroofstrata;fracture;catastrophe;collapse

　　采空區(qū)頂板穩(wěn)定性狀況由礦柱和頂板2個(gè)基本要素共同決定[1].各個(gè)巖層在地下開(kāi)采中所發(fā)揮的作用是不同的.有些厚度較大且較為堅(jiān)硬的巖層,在整個(gè)巖體運(yùn)動(dòng)過(guò)程中起控制作用(支撐載體作用),而有些較為軟弱、厚度較大的巖層在整個(gè)巖體運(yùn)動(dòng)中只起加載作用.覆巖的大部分自重載荷主要由堅(jiān)硬的厚巖層承擔(dān),它的穩(wěn)定性狀況主導(dǎo)著整個(gè)地下開(kāi)采上覆巖層運(yùn)動(dòng)的模式和規(guī)模.因此,通常將巖體活動(dòng)中起主要控制作用的巖層稱為關(guān)鍵層

礦區(qū)為例[3-4][2].頂板關(guān)鍵層對(duì)覆巖運(yùn)動(dòng)的全局起著至關(guān)重要和決定性的作用.以大同,礦區(qū)屬侏羅紀(jì)煤系頂板巖體結(jié)構(gòu),大多數(shù)采場(chǎng)頂板屬于堅(jiān)硬或極堅(jiān)硬的頂板巖體.堅(jiān)硬難冒頂板巖體分層厚度大、整體性強(qiáng).1961年挖金灣煤礦(房柱式)和1975年馬脊梁礦(刀柱式)曾先后發(fā)生采空區(qū)大面積頂板瞬時(shí)一次垮落,造成井下設(shè)施嚴(yán)重?fù)p毀.2005年邢臺(tái)縣尚汪莊石膏礦區(qū)發(fā)生特別重大坍塌事故,造成重大人員傷亡和財(cái)產(chǎn)損失.采空區(qū)塌陷是由多種原因造成支承煤柱和頂板破裂失穩(wěn)的

收稿日期:2007-09-05　　責(zé)任編輯:柴海濤

　　基金項(xiàng)目:教育部高等學(xué)校博士學(xué)科點(diǎn)專(zhuān)項(xiàng)科研基金資助項(xiàng)目(20040008025)

　　作者簡(jiǎn)介:王金安(1958—),男,河北昌黎人,教授,博士生導(dǎo)師.E-mail:wjarock@ces.ustb.edu.cn

第8期王金安等:采空區(qū)堅(jiān)硬頂板破斷機(jī)理與災(zāi)變塌陷研究

[1,5-7]851交互過(guò)程.本文在彈性基礎(chǔ)板理論基礎(chǔ)上,引入非線性分析,研究頂板和礦柱漸進(jìn)破壞導(dǎo)致采空區(qū)堅(jiān)硬頂板突然失穩(wěn)和塌陷的力學(xué)機(jī)理與失穩(wěn)過(guò)程.

1　

采空區(qū)頂板破壞的力學(xué)分析

1.1　力學(xué)模型

房柱式或條帶式采煤方法在采空區(qū)內(nèi)留下大量的煤柱,

這些煤柱支撐著堅(jiān)硬難冒頂板巖體.當(dāng)頂板和煤柱接近或

超過(guò)極限荷載狀態(tài),將導(dǎo)致頂板大面積突然垮落,造成礦

山地質(zhì)災(zāi)害.研究頂板破壞機(jī)理可在深層次上認(rèn)識(shí)堅(jiān)硬難

冒頂板大面積瞬時(shí)垮落這一礦山壓力現(xiàn)象的本質(zhì).

采空區(qū)垂直剖面如圖1所示,可將采空區(qū)頂板巖體視

為邊界固定的彈性矩形平板(圖2).建立坐標(biāo)系oxyz,設(shè)圖1　采空區(qū)垂直剖面Fig.1　Crosssectionalprofileofroom-pillarminingsystem

彈性矩形平板長(zhǎng)度為2a、寬度為2b(b≤a)、厚度為h,頂板巖體的彈性模量為E、泊松比為ν、體密度為ρ、抗拉強(qiáng)度極限為[σ].設(shè)上層巖土介質(zhì)對(duì)頂板上表面的壓力為均布載荷qs0

.

圖2　頂板巖體簡(jiǎn)化為四邊固支的等效彈性基礎(chǔ)上的平板

Fig.2　Roofstratumissimplifiedasanequivalentelasticplatewithfoursidesfixedonelasticfoundation

　　基本假設(shè):①采空區(qū)尺寸已達(dá)到或超過(guò)頂板關(guān)鍵層厚度的5～8倍以上,采空區(qū)上方堅(jiān)硬頂板近似按薄板處理;②采空區(qū)上方堅(jiān)硬頂板在破斷前為小變形,其上覆巖載荷在頂板變形過(guò)程中不發(fā)生明顯變化,故作用在堅(jiān)硬頂板上的覆巖荷載按恒載荷對(duì)待;③煤層堅(jiān)硬,采空區(qū)四周煤壁的微小塑性變形對(duì)堅(jiān)硬頂板的支撐影響作用忽略不計(jì);④將采空區(qū)內(nèi)每個(gè)煤柱視為相同的受壓彈性直桿,其初始小變形時(shí)的彈性彈性模量為E,高度為H.假設(shè)煤柱是等距分布的,其總數(shù)目為n.將彈性直桿近似1,平均橫截面積為A

地等效成連續(xù)分布的溫克爾彈性基礎(chǔ)

　　根據(jù)彈性基礎(chǔ)上的平板彎曲理論

式中,D為板的抗彎剛度,D=E3[7],等效彈性系數(shù)為k,于是有nE/H=4abk.得到等效彈性系數(shù)1Ak=nE/(4abH).(1)1A,頂板下沉位移(撓度)w(x,y)滿足4D w+kw=q,2[8](2)[12(1-ν)].

作用在頂板上總的均布載荷q由qgh疊加而成,即q=qρgh.0與頂板的自重載荷ρ0+

頂板在破壞前的邊界條件為固定:x=±a=0,y=±b=0,=0,=0.x=±a y=±b

1.2　頂板邊緣破裂的強(qiáng)度分析

根據(jù)頂板的下沉位移(撓度)和邊界約束條件情況,近似假設(shè)撓度具有如下解析解的形式

w0222222w44(x-a)(y-b).ab

　　將其代入式(2)的伽遼金弱形式方程,

(中心撓度)為[8],即(3)

852煤　　炭　　學(xué)　　報(bào)2008年第33卷　　板的彎矩表達(dá)式為

224Dw0 w 2222222222M=-+=(3x-a)(y-b)+ν(x-a)(3y-b,x22ab x 224Dw0 w 2222222222M-2+2=(x-a)(3y-b)+ν

(3x-a)(y-b).y= y ab

　　在板的邊界中點(diǎn)(±a,0)和(±b,0)彎矩分別達(dá)到最大值: MwWx max=8D0/a和 y max=8Dw,且彎曲應(yīng)力最大.根據(jù)強(qiáng)度理論,頂板是在邊界中點(diǎn)處首先進(jìn)入塑性狀態(tài)而發(fā)生初始破壞的,0/b

其初始破壞的條件為

σxmax6 M48Dw6 M48Dwx max0y max0≥[σ],σ≥[σ].symaxshahhbh(4)22

　　將式(1),(3)代入式(4),可得如下等價(jià)的初始破壞條件,即

9Haq4422E≤[E]1(7a+7b+4ab,1A1An4bh[σs]ab

4422E≤[E]21A1A(7a+7b+4ab.n4ah[σs]ab(5)(6)

　　隨著煤柱流變和表面的風(fēng)化,彈性模量E減小,致使煤柱的有效抗壓剛度E減1和平均橫截面積A1A小,它的支承能力逐漸減小.由于b≤a,比較式(5),(6),有[E]2≤[E]1.當(dāng)E減小到臨界1A1A1A值[E]2時(shí),頂板的長(zhǎng)邊中點(diǎn)先進(jìn)入塑性狀態(tài),而后沿長(zhǎng)邊擴(kuò)展形成塑性鉸;當(dāng)E繼續(xù)減小到臨界值1A

1A

[E]1時(shí),頂板的短邊中點(diǎn)也進(jìn)入塑性狀態(tài),也沿邊界形成塑性鉸,最后邊界條件由固定支承變成了簡(jiǎn)1A

支(圖3).此時(shí),頂板只是邊界約束減弱了,但整體并沒(méi)有破壞.這一過(guò)程為頂板破壞的第1階段.

圖3　頂板破裂過(guò)程

Fig.3　Thefractureprocessandthepositionofyieldinghingesinroofstratum

1.3　頂板內(nèi)部破裂的極限分析

四邊鉸支的邊界條件為

x=±a=0,y=±=0,=0,22ba xx=± y

πxπy　　假設(shè)撓度的近似解析解為w=wco,其中0c2a2b

416w+0k16aπ22y=±b=0..(7)

　　彎矩的絕對(duì)值和彎曲應(yīng)力在頂板的中心點(diǎn)(0,0)處達(dá)到最大值,因此,破壞條件為

221ν1σ],σ].xmax0≥[σsymax0≥[σs22w22w2ha2ha　　將式(1

,即

433,(8)(9)

第8期王金安等:采空區(qū)堅(jiān)硬頂板破斷機(jī)理與災(zāi)變塌陷研究

44H222b(a+b.24nah[σ]16abs853E≤[E]4

1A1A(10)

　　比較式(9),(10),有[E]4≤[E]3.當(dāng)煤柱實(shí)際的抗壓強(qiáng)度E降低到臨界值[E]3時(shí),1A1A1A1A頂板的中心下表面處σ當(dāng)煤柱抗壓強(qiáng)度E小于臨界值[E]4時(shí),頂板中y較σx先達(dá)到最大拉應(yīng)力值.1A1A

心點(diǎn)處會(huì)產(chǎn)生沿x軸的塑性鉸,并向兩側(cè)擴(kuò)展形成內(nèi)部塑性鉸線(破裂線),塑性鉸線擴(kuò)展到一定長(zhǎng)度后,會(huì)進(jìn)一步分岔成4條向角點(diǎn)擴(kuò)展的塑性鉸線,最后頂板成為可繞塑性鉸線轉(zhuǎn)動(dòng)的機(jī)構(gòu),即頂板內(nèi)部發(fā)生“X”型破裂破壞(圖3(c)),這與礦山頂板的破斷過(guò)程的相似材料模擬試驗(yàn)結(jié)果相一致

[9][2,7].根據(jù)結(jié)構(gòu)極限分析的破裂線理論(又稱機(jī)動(dòng)法),外部載荷所做的外功等于內(nèi)部塑性鉸線上所消耗

2的內(nèi)功與彈性基礎(chǔ)內(nèi)儲(chǔ)存的彈性勢(shì)能之和.作用在塑性鉸線上單位長(zhǎng)度的彎矩為M]h/4,頂板s=[σs

極限狀態(tài)下的最大(中心)撓度為w.假設(shè)與短邊相聯(lián)的板塊三角形是直角三角形,則塑性鉸線上所消s

耗的內(nèi)功A分別為n和外部載荷所做的外功Aw

A(a/b-1)w]h(a/b-1)w,A(a-b/3)bw.n=4Mss=[σssw=2qs

　　彈性基礎(chǔ)內(nèi)儲(chǔ)存的彈性勢(shì)能為

U=2k(0[w(1-y/b)]dy)dx+4k(-(a[w(1-y/b)]dy)dx+ss-(a-b)a-bx-b)

ba2

0y+a-bs2s2[9]∫∫∫∫2kb(2a-b)w.4k([w(a-x)/b]dx)d∫∫32(a-b)bab2

　　因此,由AAU可得到頂板在極限狀態(tài)時(shí)煤柱抗壓強(qiáng)度的臨界值為w=n+

6Ha2[E]s{2q(a-b/3)b-[σ]h(a/b-1)}.1Asn(2a-b)ws

　　當(dāng)煤柱的抗壓強(qiáng)度E降低到臨界值[E]s時(shí),頂板達(dá)到極限狀態(tài)而產(chǎn)生內(nèi)部破裂破壞.從邊緣形1A1A

成塑性鉸線到內(nèi)部形成塑性鉸線而成為機(jī)構(gòu),這是頂板破壞的第2階段.盡管在第2階段結(jié)束(極限狀態(tài))時(shí),頂板內(nèi)部已破裂成多個(gè)板塊組成的機(jī)構(gòu)已完全破壞.但是,由于煤柱仍有抗壓支承能力,頂板并沒(méi)有全局塌陷而失穩(wěn).

2　頂板-煤柱系統(tǒng)整體失穩(wěn)分析

頂板在第2破壞階段結(jié)束時(shí),破裂的頂板中部區(qū)域的下陷會(huì)使其邊緣沿水平向內(nèi)滑移.這里考慮一種極端危險(xiǎn)的情形,即頂板邊緣從采空區(qū)邊界幫部脫落下來(lái),從而失去邊界幫部對(duì)頂板的支承.破裂后的頂板板塊的內(nèi)部變形能被全部釋放,同時(shí)頂板板塊的載荷和自重全由煤柱支撐.此時(shí),頂板的下沉量w僅與煤柱的支撐情況(剛度)有關(guān),而與頂板剛度D無(wú)關(guān),這樣作為一種簡(jiǎn)化的分析,式(2)退化為