唐禮忠，范若楠，宋徉霖

(中南大學(xué)資源與安全工程學(xué)院，長(zhǎng)沙 410083)

深埋巷道與淺埋巷道在周邊圍巖的破壞形式上有著較大的不同。淺埋巷道周邊圍巖在破壞時(shí)有著良好的線性特征，從巷道壁沿徑向可依次形成松動(dòng)區(qū)、塑性區(qū)及彈性區(qū)，因此對(duì)淺埋巷道的研究主要是使用連續(xù)介質(zhì)中的彈塑性力學(xué)理論[1-3]。而深埋巷道圍巖破壞時(shí)，巖體在洞室徑向上的主裂縫帶交替存在，且各主裂縫帶之間存在巖體相對(duì)完整的區(qū)域，這種現(xiàn)象被稱為分區(qū)破裂化現(xiàn)象[4]。該現(xiàn)象在俄羅斯Taimyrskii的深層礦洞中被首次發(fā)現(xiàn)[5]，此后，在南非、烏克蘭及中國(guó)安徽等地的深層礦井中均觀測(cè)到類似的分區(qū)破裂化現(xiàn)象。

在實(shí)驗(yàn)研究方面，顧金才等[11]首次利用在內(nèi)含巷道的特制材料模型上施加圍壓的方式，得出分區(qū)破裂化現(xiàn)象的主要成因是高軸向壓力的結(jié)論。張強(qiáng)勇等[12]利用類似實(shí)驗(yàn)測(cè)量出模型內(nèi)部的應(yīng)變及位移的波浪形變化規(guī)律。

在數(shù)值模擬方面，錢七虎等[13]采用二維數(shù)值模擬，再現(xiàn)了中國(guó)錦屏II級(jí)水電站引水隧洞的破壞過(guò)程，發(fā)現(xiàn)在引水隧洞周圍產(chǎn)生了滑移帶狀地層裂縫。李樹(shù)忱等[14]在結(jié)合能量理論的基礎(chǔ)上，利用FLAC3D軟件對(duì)錦屏深部圍巖進(jìn)行了數(shù)值模擬。王紅英等[15]對(duì)應(yīng)變軟化模型巖體的變形機(jī)理進(jìn)行了分析，并用FLAC2D軟件模擬了地下洞室的應(yīng)力變形狀態(tài)。陳旭光等[16]利用ABAQUS軟件內(nèi)置擴(kuò)展有限元功能對(duì)分區(qū)破裂化現(xiàn)象進(jìn)行了模擬[16]。

在上述數(shù)值模擬中，分區(qū)破裂化現(xiàn)象的出現(xiàn)情況基本與實(shí)際情況吻合，但模擬時(shí)往往忽略了開(kāi)挖進(jìn)程對(duì)分區(qū)破裂化的影響，無(wú)法反映在不同開(kāi)挖階段時(shí)，分區(qū)破裂化現(xiàn)象的產(chǎn)生過(guò)程及破壞程度。而在實(shí)際工程施工時(shí)，在不同的開(kāi)挖階段都能保證施工的安全性是非常重要的。

故在上述研究的基礎(chǔ)上，結(jié)合安徽冬瓜山銅礦開(kāi)挖時(shí)的實(shí)際情況，利用ABAQUS有限元模擬軟件模擬了洞室在不同開(kāi)挖階段時(shí)的分區(qū)破裂化現(xiàn)象，對(duì)不同階段下的圍巖破壞形式進(jìn)行了研究。

1 實(shí)際工程中的分區(qū)破裂化現(xiàn)象

由于分區(qū)破裂化現(xiàn)象無(wú)法通過(guò)直觀方式觀察，通常存在兩種間接觀察方式：一是在隧道壁上鉆孔，通過(guò)分析鉆孔得到巖芯的圖像來(lái)推測(cè)破裂帶的存在；二是測(cè)量圍巖的物理力學(xué)參數(shù)(如聲速)并分析其變化規(guī)律來(lái)推測(cè)破裂帶的存在。因此，對(duì)冬瓜山銅礦巷道圍巖鉆孔得到的巖芯的完整度、采用率、缺失情況等信息進(jìn)行分析，從而推測(cè)出巷道周邊圍巖中是否存在分區(qū)破裂化現(xiàn)象。

安徽省銅陵有色集團(tuán)下屬冬瓜山銅礦是我國(guó)最早深入巖層內(nèi)部的金屬礦山之一，其平均深度達(dá)800 m，是我國(guó)目前采掘深度最深的銅礦床，其巷道圍巖內(nèi)最大主應(yīng)力已達(dá)38 MPa[17-18]，并伴隨著高應(yīng)變、高變形的特點(diǎn)。其中較為典型斷面(46-4#采場(chǎng)B斷面)的圍巖探測(cè)鉆孔布置如圖1所示。

注：1～5表示鉆孔編號(hào)圖1 冬瓜山銅礦46-4#采場(chǎng)B斷面探測(cè)鉆孔布置Fig.1 Drilling layout of the section B 46-4# stope in Dongguashan Copper Mine

通過(guò)對(duì)鉆孔取得的巖芯樣本數(shù)據(jù)的分析，將巖芯尺寸按一定數(shù)量級(jí)劃分為破裂狀態(tài)和相對(duì)完整狀態(tài)，再將不同孔中破裂狀態(tài)區(qū)域間用弧線相連，即可大致預(yù)測(cè)圍巖中的分區(qū)破裂化現(xiàn)象的情況。對(duì)冬瓜山巷道圍巖的多次調(diào)查表明，其深部巷道圍巖疑似發(fā)生分區(qū)破裂化現(xiàn)象。其中46-4#采場(chǎng)B斷面分區(qū)破裂化情況如圖2所示。

圖2 46-4#采場(chǎng)B斷面預(yù)測(cè)巷道圍巖分區(qū)破壞情況Fig.2 The prediction of surrounding rock zonal disintegration of roadway in section B 46-4 stop

由圖2知，圍巖基本被分為4層破碎帶：第1層緊貼巷道壁，以巷道壁為原點(diǎn)，徑向?yàn)檎较颍瑒t第1層破碎帶范圍在0～1.2 m，其破碎程度嚴(yán)重，可被看作傳統(tǒng)上的巷道圍巖松動(dòng)區(qū)；第2層破碎帶范圍在1.6～2.1 m；第3層破碎帶范圍在2.6～3.4 m；第4層破碎帶范圍在3.7～4.6 m。由于該斷面左側(cè)存在礦區(qū)采場(chǎng)對(duì)該斷面的持續(xù)擾動(dòng)，使已有裂縫持續(xù)擴(kuò)張發(fā)展，導(dǎo)致該斷面左側(cè)部分的裂縫帶比右側(cè)部分的裂縫帶普遍更寬。

1.選擇合適的后備豬，一般5～6月齡60 kg左右體重為宜，奶頭7對(duì)以上，大小勻稱，外陰發(fā)育良好，四肢強(qiáng)健，性情溫順的后備豬中選擇。

原巖測(cè)試實(shí)驗(yàn)得到該現(xiàn)場(chǎng)的地應(yīng)力參數(shù)如表1所示。

表1 巷道現(xiàn)場(chǎng)地應(yīng)力參數(shù)

2 數(shù)值模擬

2.1 有限元模型

鑒于所模擬巷道的埋深較深，故與其橫斷面方向平行的地應(yīng)力可看作相同。為提高有限元計(jì)算效率，本模擬僅截取巷道及其較近處的圍巖進(jìn)行計(jì)算。模擬所采用的有限元模型為立方體，其尺寸(長(zhǎng)×寬×高)為40 m×14 m×40 m(見(jiàn)圖3)。

圖3 巷道的有限元模型Fig.3 Finite element model of the roadway

冬瓜山礦洞開(kāi)挖時(shí)采用了毫秒延時(shí)爆破方法，以毫秒級(jí)時(shí)差順序起爆同一斷面上的各組炸藥包。該爆破模式可實(shí)現(xiàn)全斷面一次爆破，提高掘進(jìn)效率，并降低對(duì)洞內(nèi)人員的健康損害。雖然各組炸藥的起爆順序存在時(shí)間間隔，但時(shí)間間隔為毫秒級(jí)，可被看作瞬時(shí)卸載，故在本有限元模擬中采用刪除單元法的開(kāi)挖方式，實(shí)現(xiàn)全斷面一次性開(kāi)挖，更加擬合實(shí)際中的掘進(jìn)模式。而在實(shí)際開(kāi)挖中，一處礦洞需分多次爆破才能完成，且礦洞長(zhǎng)度較長(zhǎng)，故模擬中分2次開(kāi)挖，第一次開(kāi)挖長(zhǎng)度為7 m(第一階段)，之后加載平衡再開(kāi)挖7 m(第二階段)。整個(gè)模擬采用“先加載，后開(kāi)挖”的模式，使模擬開(kāi)挖時(shí)的圍巖應(yīng)力也更加符合實(shí)際情況。

實(shí)際冬瓜山銅礦46-4#采場(chǎng)B斷面的界面尺寸如圖4所示。

圖4 巷道的斷面尺寸Fig.4 Section size of the roadway

模型上加載垂直及水平方向地應(yīng)力采用表1中的數(shù)據(jù)，軸向應(yīng)力實(shí)際開(kāi)挖時(shí)遠(yuǎn)大于現(xiàn)場(chǎng)測(cè)試時(shí)原巖應(yīng)力，故采用較大軸向應(yīng)力。模擬采用Mohr-Coulomb本構(gòu)模型，其材料參數(shù)如表2所示。

表2 Mohr-Column模型參數(shù)

該有限元模型共有21 350個(gè)單元，單元類型采用C3D20。軟件分析時(shí)，在開(kāi)挖后共設(shè)置60個(gè)時(shí)步令模型回到平衡狀態(tài)。分析時(shí)采用非對(duì)稱求解器進(jìn)行計(jì)算，開(kāi)啟允許大變形影響選項(xiàng)。該模擬主要研究三軸受壓對(duì)分區(qū)破裂化的影響，因此真實(shí)圍巖中的不連續(xù)性特征如節(jié)理和裂縫等不在模擬范圍之內(nèi)。

2.2 分區(qū)破裂化的產(chǎn)生過(guò)程

ABAQUS軟件可反映模型中發(fā)生的等效塑性應(yīng)變(Equivalent Plastic Strain)情況。在云圖中，等效應(yīng)變值越高的區(qū)域，表明模型中該區(qū)域較其他區(qū)域所受的破壞更嚴(yán)重。各開(kāi)挖階段結(jié)束后的模擬結(jié)果見(jiàn)圖5所示。

圖5各開(kāi)挖階段巷道圍巖等效塑性應(yīng)變Fig.5 Equivalent plastic strain of surrounding rock of roadway in each excavation stage

由圖5a知，隧道底部和兩側(cè)向內(nèi)隆起，說(shuō)明其等效塑性應(yīng)變最大處為開(kāi)挖工作面與隧道的交界處。也說(shuō)明在此階段分區(qū)破裂化現(xiàn)象已經(jīng)產(chǎn)生，其第1層破裂帶可看作傳統(tǒng)破壞區(qū)域。在傳統(tǒng)破壞區(qū)域與第2層破裂帶之間存在相對(duì)完整區(qū)域。從模型的側(cè)面來(lái)看，破裂帶從開(kāi)挖工作面附近向已開(kāi)挖區(qū)域延伸出來(lái)，呈漏斗形狀。

由圖5b知，模型中新產(chǎn)生了4層清晰的破裂帶，而在開(kāi)挖第一階段中產(chǎn)生的第2層破裂帶已經(jīng)并入了傳統(tǒng)開(kāi)挖破壞區(qū)(第1層破裂帶)，說(shuō)明隨著時(shí)間的推移，破壞程度逐漸加深。在外側(cè)的3層非傳統(tǒng)破裂帶中，第3層的破裂帶相對(duì)破壞程度最嚴(yán)重，呈明顯的閉合環(huán)狀。第4層破裂帶并未形成閉合環(huán)狀，但其與第3層破裂帶之間的相對(duì)完整區(qū)域比較明顯。從側(cè)面來(lái)看，分區(qū)破裂帶的發(fā)展路徑與第一階段類似，分區(qū)破裂帶呈漏斗狀由巷道壁附近向外延伸出來(lái)，但由于巷道開(kāi)挖是貫穿整個(gè)模型的緣故，各破裂帶的發(fā)展路徑出現(xiàn)了交叉的情況。

由模擬結(jié)果可知，分區(qū)破裂化的形成是伴隨著隧道的開(kāi)挖進(jìn)程，其破壞的嚴(yán)重程度與開(kāi)挖的深度及開(kāi)挖后經(jīng)過(guò)的時(shí)間呈正相關(guān)。

各開(kāi)挖階段結(jié)束后，在巷道圍巖中沿平行巷道縱軸方向切應(yīng)力(τyz)云圖如圖6所示。

比較圖6a和圖5a，在第一開(kāi)挖階段結(jié)束時(shí)，τyz在第2層破壞帶的內(nèi)邊緣處的值高于其他位置的值，推測(cè)該處較高的應(yīng)力可能是引起第2層破裂帶的主要原因。比較圖6b和圖5b，τyz在各主破裂帶位置與其兩側(cè)圍巖的值相同但方向相反，推測(cè)這是引發(fā)破裂帶但保持其兩側(cè)存在相對(duì)完整區(qū)域的原因。

3 結(jié)果對(duì)比

在冬瓜山銅礦46-4#采場(chǎng)B斷面中，共形成了4層破裂帶。在模擬結(jié)果中，最終形成了4層清晰的破裂帶，在模擬的第一階段中形成的第2層破裂帶已并入隧道的傳統(tǒng)破壞區(qū)。將模擬中各階段分區(qū)破裂數(shù)據(jù)與實(shí)際測(cè)量數(shù)據(jù)進(jìn)行對(duì)比(見(jiàn)表3)。

表3 數(shù)值模擬結(jié)果與實(shí)際測(cè)量結(jié)果對(duì)比

從模擬數(shù)據(jù)來(lái)看，除第2階段的第1層破裂帶的寬度外，其余破裂帶的寬度基本與實(shí)際測(cè)量類似，但模擬結(jié)果中的破裂帶范圍普遍較實(shí)測(cè)結(jié)果更靠外側(cè)，推測(cè)可能是因?yàn)槟M忽略了真實(shí)巖石中的不連續(xù)性特征如既有裂縫和節(jié)理等，導(dǎo)致模擬中的主裂縫帶更不易發(fā)展。此外，由于實(shí)際情況中，巷道西側(cè)(圖1中巷道左側(cè))存在其他采場(chǎng)對(duì)巷道的持續(xù)擾動(dòng)，導(dǎo)致實(shí)際巷道較本模擬更容易發(fā)生分區(qū)破裂化現(xiàn)象。

4 結(jié)語(yǔ)

1)數(shù)值模擬的結(jié)果顯示了在深層巷道的開(kāi)挖過(guò)程中，其周邊圍巖產(chǎn)生4層分區(qū)破裂帶，前3層破壞較為均勻，且破壞具有方向性，距掘進(jìn)面越近破壞越嚴(yán)重。

2)分區(qū)破裂化現(xiàn)象產(chǎn)生的嚴(yán)重程度不僅與巖石的力學(xué)屬性及地應(yīng)力大小有關(guān)，也與開(kāi)挖隧道的長(zhǎng)度以及開(kāi)挖后經(jīng)過(guò)的時(shí)間有關(guān)。開(kāi)挖隧道的長(zhǎng)度越長(zhǎng)，開(kāi)挖后經(jīng)過(guò)的時(shí)間越久，分區(qū)破裂化現(xiàn)象越嚴(yán)重；且分區(qū)破裂化現(xiàn)象是伴隨著開(kāi)挖進(jìn)程不斷加深的。

3)數(shù)值模擬中可能會(huì)忽略真實(shí)巖石中的不連續(xù)特征，導(dǎo)致模擬結(jié)果較真實(shí)情況出現(xiàn)偏差，分區(qū)破裂帶范圍普遍距巷道更遠(yuǎn)。在未來(lái)的數(shù)值模擬中可以對(duì)此加以改進(jìn)。

4)數(shù)值模擬和現(xiàn)場(chǎng)實(shí)測(cè)數(shù)據(jù)進(jìn)行對(duì)比可知，該模擬的擬合效果較好，對(duì)深層巷道的模擬和分區(qū)破裂化現(xiàn)象的預(yù)測(cè)具有一定的參考意義。