巷道圍巖爆破地震效應模擬分析

王驥俊，龍 賀

(1.山西汾西礦業集團 柳灣煤礦，山西 孝義 032300； 2.山東井亭實業有限公司 生產技術部,山東 棗莊 277525)

·試驗研究·

巷道圍巖爆破地震效應模擬分析

王驥俊1，龍 賀2

(1.山西汾西礦業集團 柳灣煤礦，山西 孝義 032300； 2.山東井亭實業有限公司 生產技術部,山東 棗莊 277525)

對某巷道圍巖爆破工程進行數值模擬，模擬出巷道圍巖中爆破地震效應，然后將數值模擬結果與理論上的速度解析解相對比，驗證了模擬的正確性。經過對比可知，軸向距離對巷道圍巖質點峰值速度的影響程度比徑向距離要小，所以爆破地震波對徑向圍巖的作用要大于對軸向圍巖的作用。因此在實際工程中要特別注意對巷道徑向圍巖的防護。

巷道圍巖；爆破地震；數值模擬分析；解析對比

隨著地下工程爆破技術的逐步應用，巷道圍巖爆破地震效應的研究愈發重要。當進行掘進爆破施工時，爆破地震波中的一部分能量向周圍巖體擴散，引起圍巖和支護結構震動[1,2]. 受爆破的復雜性和測試技術的有限性等的影響，巷道內部力學和變形性狀情況不容易觀測，從而影響施工安全。地震效應強度、圍巖的巖石特性、巷道斷面大小等是影響巷道圍巖及支護結構爆破地震效應的關鍵因素。國內外的專家學者對巷道圍巖爆破地震效應進行了很多研究，在數值模擬方面，通過模擬淺孔爆破得到質點速度，經對比，與實際檢測結果相吻合；通過編程對爆炸應力波進行計算，建立本構關系；此外還分別采用不同方法對巷道圍巖爆破效應進行模擬[3]. 但這些研究都缺乏實際論證，不適用于實際工程。在理論方面，通過分離變量，得到巷道圍巖中地震波的傳播速度解析解，從而得到圍巖質點峰值速度的變化規律。本文將通過運用ANSYS/LS-DYNA有限元分析軟件，對巷道圍巖爆破地震效應進行模擬，將模擬結果與理論結果相對比，研究爆破地震波在巷道圍巖中的傳播規律[4].

1 計算模型的建立

計算模型采用某巷道掘進爆破工程，其爆破施工技術采用微差爆破，炮孔半徑為2 cm，炮孔之間的距離即炮孔間距為45 cm. 巷道簡化為圓形，其直徑為5 m，循環進尺為1.5 m.

1.1 模擬材料選取

本計算模型采用的是乳化炸藥，其炸藥材料模型參數和狀態方程參數見表1.

表1 炸藥材料及狀態方程參數表

本文計算模型中巷道圍巖為風化的花崗巖，為了保證模擬的正確性，將巖體的本構模型加以簡化，定義其為連續介質線彈性模型，巖石材料的主要參數見表2.

表2 花崗巖主要參數表

1.2 爆破等效荷載的處理方法

由于巷道掘進爆破主要以鉆孔裝藥爆破為主，為了能全面模擬出整個爆破地震效應，經過大量的研究總結出，對于爆源附近巖體的動力響應問題的模擬采用直接模擬裝藥法較為理想，而對于爆破中遠區域爆破震動的模擬研究需要采用爆破荷載等效施加的方法，將爆破荷載進行簡化，簡化后的結果不僅要和采用直接模擬裝藥法得到的結果盡量接近，而且又要節省計算工作量。

目前，對爆破震動的傳播進行數值模擬，可采用的方法基本分為兩種：

1) 使用ANSYS/LS-DYNA程序中的高性能炸藥材料模型及其狀態方程來模擬炸藥的整個爆轟過程及爆轟產物與周圍巖石介質的相互作用過程。

2) 根據爆轟波理論和爆腔膨脹理論，結合工程經驗、實測結果、炸藥參數等得到半理論半經驗的爆破荷載壓力曲線，此曲線必須能最大程度地反映出炮孔壓力的變化歷程，然后將該荷載曲線直接施加在炮孔壁上。在實際計算應用中，爆破荷載壓力曲線的選取主要包括三角型脈沖荷載曲線和指數衰減型荷載曲線兩種。

第一種方法較第二種方法而言，模擬結果更具有真實性，但其只能較好地模擬出爆源近區的沖擊響應，對于模擬爆破中遠區的爆破動力響應，存在諸多問題。第二種方法雖然方便，但是對于炮孔周圍尤其是炮孔數很多情況下的網格劃分工作仍然十分復雜，且爆破荷載曲線是通過理論結合試驗得到的，對于解決實際問題仍存在很大的誤差。結合以上兩種方法，基于圣維南原理，提出了合理的爆破荷載等效施加方法，使其盡可能接近實際爆破荷載作用。即在不考慮炮孔形狀的前提下，提出合理的爆破荷載曲線，然后對其等效施加在對應的開挖輪廓面上，此開挖輪廓面是由同排炮孔連心線與炮孔軸線所確定的。

1.3 邊界條件的選取

采用平直的三維黏彈性人工邊界，即同一邊界面上的法向或切向具有相同的參數值[3].巷道圍巖爆破地震效應的計算模型簡圖見圖1，圖2.

圖1 巷道計算模型正視圖

圖2 巷道計算模型側視圖

2 模擬結果與分析

由于爆破具有瞬時性的特點，在本次模擬中，爆破等效荷載收縮速度較快，當時間達到0.002 s時，荷載已經趨近于零；同時爆源附近圍巖的受力變形也在短時間內完成，所以定義模擬的求解控制時間為0.01 s，以便觀察爆破地震效應的整個傳播過程[3].

2.1 圍巖質點速度云圖

對模型進行求解，為了更加直觀地反映出模擬的效果，現列出部分時間點的圍巖質點軸向速度云圖，見圖3～7.

圖3 t=0.003 977 2 s時圍巖質點軸向速度云圖

圖4 t=0.004 696 9 s時圍巖質點軸向速度云圖

圖5 t=0.005 07 s時圍巖質點軸向速度云圖

圖6 t=0.006 174 2 s時圍巖質點軸向速度云圖

圖7 t=0.008 863 7 s時圍巖質點軸向速度云圖

從圖3～7可以看出，當爆炸發生后，爆破地震波在圍巖中以爆源為原點開始沿著徑向傳播時，圍巖質點軸向速度隨著徑向距離的增大而有一定規律的遞減；且呈現正負交替狀態，對其速度取最大值和最小值的絕對值相比較，發現二者相差不大。此外，通過研究在爆破地震波作用下圍巖質點的徑向速度，發現其也有類似特點。這與實際爆破工程中所呈現的狀況相吻合，說明模擬結果是合理的。

2.2 模擬結果對比分析

通過查閱文獻，對于巷道圍巖在動力作用下爆破地震效應的傳播規律，速度的解析解方程為[4,5]：

(1)

(2)

其中：

(3)

式中：

m、n、A1、A2—實常數；

R—巷道半徑，m；

b—循環進尺，m；

cp—波速，m/s；

r—徑向距離，m；

z—軸向距離，m；其他參數含義同上。

將相關數據帶入式(1)，式(2)，得到巷道圍巖速度曲線。然后將模擬中得到的圍巖質點峰值速度整理并用曲線的形式表示出來，并對二者進行對比分析，見圖8～11.

通過觀察圖8～10，經對比可知從數值和曲線變化規律的角度上分析，模擬結果和解析結果是相似的，說明模擬結果是可靠的。

圖8 Z=3 m時巷道圍巖質點徑向峰值速度隨徑向距離的變化曲線圖

圖9 Z=3 m時巷道圍巖質點軸向峰值速度隨徑向距離的變化曲線圖

圖10 r=2.5 m時巷道圍巖表面質點徑向速度隨軸向距離的變化曲線圖

圖11 r=6 m時巷道圍巖質點軸向速度隨軸向距離的變化曲線圖

通過對比圖8中的兩組曲線，在軸向距離為3 m時，從數值上分析，兩種情況中的巷道圍巖質點徑向峰值速度的最大值是相近的，約為15 cm/s；從曲線變化規律上分析，兩者均呈現波浪式衰減，且衰減周期距離約為8 m.

通過對比圖9中的兩組曲線，發現均呈現上述特點。結果表明：在爆破地震波作用下，巷道圍巖質點的徑向和軸向峰值速度隨著徑向距離的增大呈現正負交替周期性變化且逐漸遞減狀態，這可能是導致深部巖體發生分區破裂現象最重要的原因，但其中的具體原因有待進一步研究。

對圖10中的兩組曲線進行對比分析，在徑向距離為2.5 m時，兩種情況下巷道圍巖質點的徑向峰值速度隨著軸向距離的增大都呈現負指數遞減狀態，且衰減速度很快，最終逐漸趨近于零。在圖11中，在徑向距離為6 m時，巷道圍巖質點的軸向峰值速度與軸向距離之間的關系同上。結果表明：在爆破地震波作用下，與徑向距離相比較，軸向距離對巷道圍巖質點徑向與軸向峰值速度的影響比較小。所以在巷道爆破工程中，爆破地震波對徑向圍巖的作用要更大一些，大于對軸向圍巖的作用。所以在巷道圍巖的爆破施工中要更加注意地震波對徑向圍巖的影響。

3 結 論

本文通過運用ANSYS/LS-DYNA有限元分析軟件，對巷道圍巖的爆破地震效應進行模擬。首先把問題簡化，然后選取巖石和炸藥材料模型，進行爆破效應等效荷載，選取邊界條件，從而建立計算和數值模擬模型，進行網格劃分，得到圍巖質點速度云圖和峰值速度曲線圖。將模擬結果與理論結果進行對比分析，發現兩種情況下的圍巖質點速度場在交替周期和衰減速度上都具有一致性，以此來研究爆破地震波在巷道圍巖中的傳播規律，得出以下結論：

1) 在爆破地震波的作用下，隨著徑向距離的增大，巷道圍巖質點峰值速度呈現波動且逐漸衰減狀態。而在深部巖土開挖工程中，如果爆破地震效應達到一定程度，使得巷道圍巖質點峰值速度達到其極限狀態時，圍巖會產生交替的破裂區和非破裂區，這與巷道圍巖質點峰值速度的波動狀態相類似，所以根據此規律有可能解釋圍巖的分區破裂現象。

2) 隨著軸向距離的增大，巷道圍巖在爆破地震波作用下的質點峰值速度會隨之減小，呈現負指數衰減，并且衰減速度比較快。并且經過對比可知，軸向距離對巷道圍巖質點峰值速度的影響程度比徑向距離要小，所以爆破地震波對徑向圍巖的作用要大于對軸向圍巖的作用。因此在實際工程中要特別注意對巷道徑向圍巖的防護。

3) 總體來說，結合不同的力學性質和不同的地質因素，對巷道圍巖中爆炸地震波的研究有助于增加巷道的穩定性和結構的安全性。

[1] 程久龍,李 飛,彭蘇萍. 礦井巷道地球物理方法超前探測研究進展與展望[J]. 煤炭學報,2014(08):82-84.

[2] 李術才,劉 斌,孫懷鳳. 隧道施工超前地質預報研究現狀及發展趨勢[J]. 巖石力學與工程學報,2014(06):632-634.

[3] 高 峰,周科平,周炳仁. 基于TRT技術的礦山井下地質超前預報[J]. 中國安全科學學報,2014(04):231-233.

[4] 相興華,鄧小鵬,王 鵬. TRT在隧道地質超前預報中的應用[J]. 地下空間與工程學報,2012(06):92-94.

[5] 王 琦,李術才,李為騰,等. 基于地質預報的煤巷頂板事故防治研究[J]. 采礦與安全工程學報,2012(01):92-94.

SimulationAnalysisofSeismicEffectforRoadwaySurroundingRock

WANGJijun,LONGHe

The numerical simulation of a rock blasting in a roadway is carried out to simulate the blasting seismic effect on the surrounding rock. The numerical simulation results are compared with the theoretical solution of the velocity analysis to verify the correctness of the simulation. The influence in the axial direction on the peak velocity of the surrounding rock is smaller than that of the radial direction measured in distance, so the effect of the blasting seismic wave on the radial rock is greater than that on the axial rock. Therefore, in the actual project to special attention should be paid to the protection of the roadway surrounding rock.

Roadway surrounding rock; Blasting seism; Numerical simulation analysis; Analytic comparison

TD235

：A

：1672-0652(2017)07-0031-05

2017-06-08

王驥俊(1985—)，男，山西平遙人，2010年畢業于中國礦業大學(北京)，工程師，主要從事煤炭開采技術研究工作(E-mail)zfg3259@163.com