張旭飛 郭運鴻 王興亞 郭奇峰

收稿日期：2023-12-28； 修回日期：2024-01-30

基金項目：國家重點研發(fā)計劃項目（2021YFC3001301）

作者簡介：張旭飛（1992—），男，工程師，從事采礦技術(shù)研究工作；E-mail：1253031969@qq.com

*通信作者：郭運鴻（2002—），男，碩士研究生，研究方向為巖石力學；E-mail：793936117@qq.com

摘要：為解決深部高應力下巷道變形嚴重及支護困難的問題，以焦家金礦-710 m水平巷道為研究對象，建立Flac3D數(shù)值分析模型，模擬分析巷道掘進過程中圍巖的擾動響應規(guī)律。結(jié)果顯示：因巷道開采擾動，巷道周邊應力得到釋放，應力進行重新分布，巷道周邊應力集中，應力極值主要在頂板、拱肩及邊幫上下兩側(cè)，應加強頂板、拱肩及邊幫上下兩側(cè)支護；距離巷道兩倍寬度范圍內(nèi)是支護的重點關(guān)注區(qū)域；巷道掘進后，周邊圍巖會出現(xiàn)剪切破壞和張拉破壞。通過對巷道變形和應力分布演化分析，揭示了掘進過程中巷道圍巖的變形失穩(wěn)過程，為支護參數(shù)設(shè)計和支護時機選擇提供了科學依據(jù)。

關(guān)鍵詞：深部開采；高應力；圍巖變形；應力分布；巷道掘進

中圖分類號：TD32????????? 文章編號：1001-1277（2024）06-0001-05

文獻標志碼：Adoi：10.11792/hj20240601

引? 言

隨著中國工業(yè)的蓬勃發(fā)展和社會需求的不斷增長，經(jīng)濟發(fā)展對資源的需求日益迫切［1-2］。為滿足社會生產(chǎn)的需求，礦產(chǎn)資源開發(fā)逐漸向深部轉(zhuǎn)移。深部地下開采工程具有“三高一擾動”的特點［3-7］。深部采掘過程中，巷道圍巖容易發(fā)生巨大變形，同時還存在高溫和突水等問題。作為礦山井巷工程的重要組成部分，井下各類巷道的安全性直接影響礦產(chǎn)資源開采的安全性，作為運輸通道的巷道更是直接制約礦山生產(chǎn)。深部高應力作用下巷道變形特征呈現(xiàn)出的擾動響應較淺部有所變化，掌握采掘過程中巷道圍巖變形和應力分布規(guī)律，是做好巷道穩(wěn)定性控制的關(guān)鍵。

眾多學者采用數(shù)值模擬分析方法對巷道圍巖穩(wěn)定性進行超前模擬分析，為工程建設(shè)和支護控制提供指導。孫宇超［8］以Flac3D數(shù)值模擬為技術(shù)手段，對等腰梯形巷道和直角梯形巷道周邊的水平應力、垂直應力和最大主應力展開了研究，發(fā)現(xiàn)無論采用等腰梯形巷道還是直角梯形巷道，巷道圍巖穩(wěn)定性的主控應力都為垂直應力。穆磊等［9-13］采用理論分析和Flac3D數(shù)值模擬相結(jié)合的方法對巷道圍巖穩(wěn)定性進行分析，研究成果為類似地質(zhì)條件下巷道的支護設(shè)計提供了很好的借鑒。任賽等［14］采用Flac3D數(shù)值模擬方法系統(tǒng)地評估了某礦山基地深井圍巖的巖爆風險程度，完成了巷道方位、掘進過程、巷道交叉、巷道與采空區(qū)相對關(guān)系等對巖爆風險等級程度、巖爆出現(xiàn)的潛在位置、開采過程中潛在巖爆發(fā)生時機的影響分析，研究結(jié)果為現(xiàn)場施工提供了指導，其建模方法和分析巖爆傾向性的方法對其他賦存復雜的礦山具有一定的借鑒意義。許永山等［15］運用Flac3D軟件，對某礦山進行了多礦房同采條件下的圍巖穩(wěn)定性分析，取得了良好的模擬效果。郭子源等［16］為解決深部高應力軟巖巷道圍巖的支護困難問題，采用Flac3D軟件有限差分程序并結(jié)合現(xiàn)場實測的方法為深層軟巖巷道的掘進和支護設(shè)計提供了技術(shù)指導。孫闖等［17］以Flac3D軟件為工具，分析應變軟化模型、剪脹角及不同巷道斷面類型對圍巖-支護相互作用的影響，并對支護系統(tǒng)穩(wěn)定性進行評價，為高應力軟巖巷道的支護結(jié)構(gòu)設(shè)計及施工提供理論支撐。

本文針對山東黃金礦業(yè)（萊州）有限公司焦家金礦（下稱“焦家金礦”）進入深部開采面臨的巷道圍巖變形和支護問題，充分考慮工況、礦巖力學特性、開采擾動等因素，結(jié)合巖石力學數(shù)據(jù)開展數(shù)值模擬分析，研究深部巷道掘進過程中的應力、變形演化規(guī)律，為支護設(shè)計和優(yōu)化提供依據(jù)。

1? 礦區(qū)概況

焦家金礦位于山東省黃縣弧形大斷裂帶，總面積25 km2，下轄3個礦區(qū)：焦家礦區(qū)、望兒山礦區(qū)和寺莊礦區(qū)。礦區(qū)內(nèi)礦體賦礦圍巖主要為玲瓏花崗巖，位于礦體下盤及深部，淺肉紅色，中?；◢徑Y(jié)構(gòu)，塊狀構(gòu)造、片麻狀構(gòu)造。主要礦物成分有斜長石（35 %～40 %）、鉀長石（25 %～35 %）、石英（25 %～30 %）、黑云母（5 %），以及少量榍石、鋯石等副礦物。礦化類型可分為浸染狀（Ⅰ號礦體）、細脈浸染狀（Ⅱ號礦體）及細脈狀（Ⅲ號礦體）。圍巖蝕變類型多樣，可分為綠泥石化、紅化、鉀長石化、絹云母化、硅化等，且有良好的分帶性。焦家金礦采用上向水平分層和上向水平進路尾砂膠結(jié)充填采礦法進行開采，深部巷道掘進過程中圍巖穩(wěn)定性控制難度大，支護結(jié)構(gòu)易破壞，返修率高。隨著開采范圍和深度的增加，焦家金礦已進入深部開采，巖體在高應力環(huán)境下的采掘擾動響應更為明顯，易出現(xiàn)破碎巖體變形大、完整硬巖動力破壞等潛在的地壓災害。2024年第6期/第45卷? 礦業(yè)工程礦業(yè)工程? 黃? 金

2? 模型建立

以焦家礦區(qū)-710 m水平巷道為工程背景，建立Flac3D數(shù)值分析模型，研究巷道掘進過程中圍巖變形及應力分布規(guī)律。數(shù)值分析模型采用的巖體物理力學參數(shù)見表1。

計算模型四周邊界與巷道中心的距離約為5倍的巷道寬度，減少邊界條件對模擬分析結(jié)果的影響。開挖計算采用位移邊界條件，模型根據(jù)區(qū)域地應力場分布特征施加梯度應力邊界，并對模型整體施加重力場，最大水平主應力與巷道軸向平行。巷道采用分步開挖，開挖總長度為16 m，每次開挖2 m。數(shù)值計算模型及網(wǎng)格劃分見圖1。

3? 巷道應力分布規(guī)律分析

應力狀態(tài)是影響巷道穩(wěn)定性的主要因素，為實現(xiàn)安全開采必須分析巷道應力分布特征。巷道開挖引起應力重新分布。巷道分步開挖完成后的圍巖應力分布云圖見圖2。由圖2可知：巷道周邊最大主應力集中于巷道拱肩和邊幫與底板拐角處，應力集中區(qū)最大應力達到80.55 MPa，極易產(chǎn)生破壞；最小主應力云圖顯示，在巷道表面產(chǎn)生一定的拉應力，極值約為1.32 MPa，易破壞剝落；水平主應力極值集中在巷道頂板處及兩幫墻底角，極值約為52.77 MPa；垂直主應力極值與前三者有所不同，主要集中在邊幫上下兩側(cè)，極值約為57.97 MPa。模擬結(jié)果顯示：隨著巷道開挖擾動，巷道周邊應力得到一定的釋放，但在局部區(qū)域產(chǎn)生嚴重的應力集中，主要在頂板、拱肩及邊幫上下兩側(cè)，施工和支護過程中應重點關(guān)注上述區(qū)域，及時排險。

在巷道邊幫及頂板中設(shè)置監(jiān)測點，分析掘進過程中應力演化規(guī)律。隨著掘進工作面的推進，巷道左幫、頂板監(jiān)測點的三向應力變化曲線見圖3。其中，σxx、σyy、σzz分別為x、y、z方向的應力。掘進過程中，垂直巷道走向的邊幫水平應力逐漸減小，隨著掘進工作面的遠離逐漸趨于穩(wěn)定，主要是開挖臨空面形成后兩幫向巷道變形，水平應力得到釋放；垂直應力隨著掘進面的推進，呈現(xiàn)逐步增加的趨勢，最后趨于穩(wěn)定，其主要原因是巷道頂部圍巖失去豎向支撐，垂直荷載向兩幫擴散轉(zhuǎn)移，導致邊幫垂直應力增大；巷道軸向應力出現(xiàn)先增大后降低的趨勢，但總體變化較小。頂板巷道軸向應力隨進尺增加先增大后降低，與邊幫一致；然而，水平應力和垂直應力變化趨勢與邊幫相反，垂直巷道走向的水平應力逐漸增加后趨于穩(wěn)定；垂直應力則先下降后緩慢增加至穩(wěn)定值，源于頂板垂直應力向兩幫的擴散及巷道豎向變形得到一定程度的應力釋放。

至巷道邊幫和頂板不同距離巖體的應力演化曲線見圖4。巷道邊幫垂直應力隨進尺增加而增大，增大幅度越來越小，最后趨于穩(wěn)定。由圖4-a）可知：越靠近巷道，垂直應力值變化越大；與之相反，遠離巷道時，圍巖垂直應力受開采擾動的影響極小。頂板垂直應力隨進尺增加出現(xiàn)先極速減小，后緩慢增大至穩(wěn)定的趨勢；距離巷道頂板越遠，垂直應力變化越不明顯，這與邊幫應力變化相同。邊幫和頂板的水平應力變化與垂直應力變化截然相反，邊幫水平應力出現(xiàn)先減小后增大的趨勢，而頂板水平應力則與邊幫垂直應力變化一致，隨著工作面的推進逐漸增大到穩(wěn)定。

4? 巷道變形規(guī)律

巷道開挖完成后的位移變化云圖見圖5。由圖5可知，開挖后巷道出現(xiàn)較為明顯的邊幫收斂現(xiàn)象。至掘進工作面距離不同時的巷道表面位移變化曲線見圖6。由圖6可知：巷道頂板和邊幫表面位移均隨掘進工作面推進逐漸增大，掘進工作面持續(xù)推進過監(jiān)測點7 m左右時，邊幫變形達到最大值并趨于穩(wěn)定，頂板變形雖持續(xù)增加，但增速明顯減緩，說明巷道表面變形在掘進工作面持續(xù)推進約巷道兩倍寬度時趨于穩(wěn)定。這為柔性支護和讓壓支護設(shè)計提供了依據(jù)和優(yōu)化思路。

at different distances to excavation face

距巷道表面不同距離巖體的位移變化見圖7。邊幫水平位移和頂板位移均隨著進尺的增加而增加，近巷道表面1 m范圍內(nèi)的變形量較為突出，一定程度上能夠反映出巷道圍巖在采掘擾動后的松動損傷范圍，為支護結(jié)構(gòu)的設(shè)計和參數(shù)確定提供參考。

5? 結(jié)? 論

采用Flac3D軟件對焦家金礦焦家礦區(qū)-710 m水平巷道圍巖開挖擾動過程進行了模擬分析，得到了巷道變形和應力的時間和空間演化規(guī)律，主要結(jié)論如下：

1）因巷道開采擾動，巷道周邊應力得到釋放，應力進行重新分布，應力集中主要聚集于頂板、拱肩及邊幫上下兩側(cè)，掘進施工過程中應加強安全管理和支護。

2）巷道頂板和邊幫的位移隨掘進工作面的推進而增加，掘進工作面后方約7 m處的巷道邊幫變形趨于穩(wěn)定，頂板下沉速度明顯降低，趨于穩(wěn)定，對柔性讓壓支護和支護時機選擇具有指導意義。

［參 考 文 獻］

［1］? 郭良銀，蔣萬飛，宋召法，等.新城金礦階段空場嗣后充填法開采大直徑深孔切槽爆破方法［J］.黃金科學技術(shù)，2022，30（4）：585-593.

［2］? 梁旭.井下采礦技術(shù)及井下采礦的發(fā)展趨勢探究［J］.能源與節(jié)能，2021（12）：149-150.

［3］? 鄭要偉.我國金屬礦山采礦技術(shù)現(xiàn)狀與發(fā)展趨勢［J］.工程建設(shè)與設(shè)計，2021（16）：33-35.

［4］? 張永政.井下采礦技術(shù)和井下采礦的發(fā)展形勢［J］.世界有色金屬，2019（8）：65，67.

［5］? 趙鉞.金屬礦深部開采現(xiàn)狀與發(fā)展戰(zhàn)略［J］.中國金屬通報，2020（9）：3-4.

［6］? 何滿潮，謝和平，彭蘇萍，等.深部開采巖體力學研究［J］.巖石力學與工程學報，2005，24（16）：2 803-2 813.

［7］? 徐冬東.金屬礦深部開采現(xiàn)狀與發(fā)展探討［J］.世界有色金屬，2020（22）：51-52.

［8］? 孫宇超.不同梯形巷道圍巖應力特征及變形規(guī)律的數(shù)值模擬研究［J］.采礦技術(shù)，2021，21（5）：113-116，132.

［9］? 穆磊，趙龍剛.基于FLAC3D數(shù)值模擬的巷道圍巖穩(wěn)定性及支護參數(shù)設(shè)計研究［J］.能源與環(huán)保，2018，40（12）：170-173.

［10］? 何富連，張廣超.深部高水平構(gòu)造應力巷道圍巖穩(wěn)定性分析及控制［J］.中國礦業(yè)大學學報，2015，44（3）：466-476.

［11］? 代碧波，李懷賓，張姝婧，等.深部采動影響下硬巖巷道圍巖穩(wěn)定性數(shù)值模擬研究［J］.金屬礦山，2021（4）：70-75.

［12］? 夏德威，冷坤鵬，張鎮(zhèn)鑠，等.基于FLAC3D數(shù)值模擬的軟巖巷道圍巖穩(wěn)定性及布置方案研究［J］.礦業(yè)研究與開發(fā)，2020，40（12）：92-99.

［13］? 范順剛，姜永恒，趙橋波，等.鎮(zhèn)沅金礦松軟破碎巖體巷道穩(wěn)定性分析及支護技術(shù)［J］.黃金，2021，42（2）：41-47.

［14］? 任賽，程光華，于寶民.基于FLAC3D的某鈾礦深井開采巖爆傾向性研究［J］.現(xiàn)代礦業(yè)，2017，33（8）：220-223.

［15］? 許永山，孫忠文，曲玉峰，等.基于FLAC3D多礦房同采圍巖的穩(wěn)定性研究［J］.現(xiàn)代礦業(yè)，2015，31（9）：7-8，11.

［16］? 郭子源，趙國彥，彭康.深部高應力軟巖巷道開挖與支護圍巖變形的FLAC3D模擬［J］.礦冶工程，2012，32（2）：18-22，29.

［17］? 孫闖，張向東，李永靖.高應力軟巖巷道圍巖與支護結(jié)構(gòu)相互作用分析［J］.巖土力學，2013，34（9）：2 601-2 607，2 614.

Deformation and stress distribution evolution of surrounding

rock during the roadway excavation in Jiaojia Gold Mine

Zhang Xufei1，Guo Yunhong2，Wang Xingya1，Guo Qifeng2

（1.Jiaojia Gold Mine，Shandong Gold Mining Industry （Laizhou） Co.，Ltd.;

2.School of Civil and Resources Engineering，University of Science and Technology Beijing）

Abstract：To solve the problem of serious deformation and difficult support of the roadway under deep high stress，the research centers on the -710 m horizontal roadway in Jiaojia Gold Mine，establishes Flac3D numerical model to simulate and analyze the disturbance response law of the surrounding rock during the roadway excavation process.The results show that due to the mining disturbance of the roadway，the stress around the roadway is released.Thus，the stress is redistributed，with the stress around the roadway concentrated，and the stress peaks mainly occurring in the roof，the arch shoulder，and the upper and lower sides of the side walls，so the support should be strengthened in the roof，the arch shoulder，and the upper and lower sides of the side walls;within a range twice the width of the roadway is the focus area for support;after excavation，surrounding rock mass will experience shear and tensile failure.By analyzing the evolution of deformation and stress distribution in the roadway，the process of deformation failure of surrounding rock during excavation is revealed，providing a scientific basis for support parameter design and support timing.

Keywords：deep mining;high stress;surrounding rock deformation;stress distribution;roadway excavation