CO₂相變致裂巖石的損傷范圍的計算及影響因素研究

［摘 要］在CO₂相變致裂巖石技術(shù)中，為探究巖石損傷范圍與CO₂相變致裂參數(shù)的關(guān)系，根據(jù)Mises準(zhǔn)則和損傷力學(xué)理論建立了CO₂相變沖擊巖石的致裂損傷半徑計算公式。利用3DEC軟件模擬了不同泄能片厚度下巖石的致裂損傷半徑，并與理論計算相比較，驗證了所建模型的準(zhǔn)確性。分析了泄能片厚度、巖石強度、致裂孔孔徑、CO₂充裝質(zhì)量等因素對巖石損傷的影響規(guī)律。結(jié)果表明：巖石的致裂損傷半徑與泄能片厚度、CO₂充裝質(zhì)量呈對數(shù)增長關(guān)系，與巖石強度、致裂孔孔徑呈指數(shù)衰減關(guān)系；增加泄能片厚度是增大巖石致裂損傷半徑的最有效手段；致裂孔孔徑的增大會使巖石的損傷范圍迅速減小。可為CO₂相變致裂巖石的損傷機制、技術(shù)參數(shù)選取等提供參考。

［關(guān)鍵詞］液態(tài)CO₂；相變致裂；損傷半徑；泄能片厚度；CO₂充裝質(zhì)量

［分類號］TD235

Damage Range Calculation and Influencing Factors of Rock Fractures Caused

by CO₂ Phase Transition Blasting

ABI Erdi^①②， HU Zunrong^①②， FU Lin^①②， HAN Yafeng^①②， LIU Mingwei^①②

① Nationl Engineering Research Center for Inland Waterway Regulation， Chongqing Jiaotong University

（Chongqing， 400074）

② Key Laboratory of Geological Hazards Mitigation for Mountainous Highway and Waterway， Chongqing

Municipal Education Commission， Chongqing Jiaotong University （Chongqing， 400074）

［ABSTRACT］In the technology of rock fractures caused by CO₂ phase transition， in order to study the relationship between the damage range of rocks and the CO₂ phase transition fracturing parameters， a calculation formula for the damage radius of rock fractures caused by CO₂ phase transition impact was established based on Mises criterion and damage mechanics theory. The damage radius of rocks with different thicknesses of energy release sheet was simulated using 3DEC software， and compared with theoretical calculations to verify the accuracy of the established model. The influence of factors such as thickness of the energy release sheet， rock strength， diameter of borehole， and CO₂ filling mass on rock damage was analyzed. The results show that the damage radius has a logarithmic growth relationship with the thickness of the energy release sheet and CO₂ filling mass， and an exponential decay relationship with the rock strength and the diameter of borehole. Increasing the thickness of the energy release sheet is the most effective way to increase the damage radius. The increase in the diameter of the borehole will rapidly reduce the damage range of rock. It can provide references for the damage mechanism and technical parameter selection in CO₂ phase trasition blasting.

［KEYWORDS］liquid CO₂ ; phase transition blasting; damage radius; thickness of energy release sheet; CO₂ filling mass

0 引言

作為一種經(jīng)濟、便捷、高效的破巖方法，爆破技術(shù)一直都是煤炭、礦山資源開采的重要手段。傳統(tǒng)的炸藥爆破在高效破巖的同時也會帶來巨大灰塵，產(chǎn)生有毒氣體、高振動沖擊波等危害居民和作業(yè)人員安全的問題。因此，一些非炸藥的綠色爆破技術(shù)相繼出現(xiàn)。其中，利用超臨界CO₂與氣態(tài)CO₂之間的壓力差作為破巖動力的CO₂相變致裂巖石技術(shù)，具有吸熱、無火花、產(chǎn)生的沖擊振動較小、產(chǎn)物為無害阻燃的CO₂氣體等優(yōu)點，在眾多復(fù)雜情況的巖體破碎開挖中得以廣泛應(yīng)用^［1-3］。

CO₂相變致裂巖石的原理主要包括2個方面^［4-6］。一方面是應(yīng)力波對巖石的沖擊作用。這是一個高溫、高壓、高速的瞬態(tài)過程，會在較短時間內(nèi)對巖石產(chǎn)生損傷破壞。也正是由于破壞的時間極為短暫，加之巖體本身的不均勻性和復(fù)雜性，使得對巖石的破碎機理的研究非常困難。另一方面是高壓CO₂氣體的氣楔作用。巖石在沖擊應(yīng)力波作用下會產(chǎn)生破碎及裂紋擴展，隨后的高壓CO₂氣體會貫入初始裂紋，并進一步促使裂紋擴展。一些學(xué)者利用液態(tài)CO₂與炸藥等當(dāng)量代換的方式，采用數(shù)值仿真來研究破巖機理，將液態(tài)CO₂相變等當(dāng)量換算為乳化炸藥^［7-8］。陳少波等^［9］使用CO₂爆破產(chǎn)生的應(yīng)力波加載時程曲線，在LS-DYNA軟件中模擬CO₂膨脹致裂巖石的過程；Zhang等^［10］分析巖石各階段的破壞規(guī)律，建立了煤層壓裂壓力和壓裂影響范圍的數(shù)學(xué)模型。還有一些學(xué)者利用實驗對CO₂相變致裂參數(shù)進行了研究，通過對比爆破應(yīng)力波的傳播和衰減規(guī)律，分析了液態(tài)CO₂爆炸中裂紋的擴展范圍、爆炸壓力、釋放孔的數(shù)量和半徑對液態(tài)CO₂爆破效果的影響^［11-13］。研究認(rèn)為，CO₂的充裝質(zhì)量和剪切片強度是決定破巖效果的2個主要因素。CO₂的充裝質(zhì)量越大，裂紋發(fā)育越多，巖塊破碎程度就越高。

綜上可知，關(guān)于CO₂相變致裂巖石的損傷范圍，當(dāng)前常以等當(dāng)量炸藥代替的方式進行分析，但二者振動差異較大，且容易受到所選炸藥參數(shù)的影響。而直接研究CO₂相變致裂巖石的損傷范圍的文獻(xiàn)較少。

本文中，通過理論分析建立了CO₂相變致裂巖石的損傷范圍的計算公式，運用3DEC離散元程序模擬了不同泄能片厚度下巖石的損傷范圍，分析了CO₂相變沖擊波下巖石損傷的演化規(guī)律，討論了巖石損傷范圍的影響因素。可為CO₂相變致裂巖石施工的參數(shù)選取等提供參考。

1 CO₂相變致裂巖石損傷的計算模型

1.1 致裂系統(tǒng)介紹

在CO₂相變致裂巖石技術(shù)中，泄能片的厚度、CO₂的充裝質(zhì)量等致裂系統(tǒng)的結(jié)構(gòu)差異，均會產(chǎn)生不同的沖擊壓力，進而影響致裂巖石的效果。CO₂相變致裂系統(tǒng)由發(fā)熱管、儲液管、泄能片、釋放管以及泄放孔等部分組成，主要結(jié)構(gòu)如圖1所示。

在工程運用中，通過充裝閥充填液態(tài)CO₂，通過發(fā)熱管內(nèi)的加熱裝置提供液態(tài)CO₂相變所需的能量，使得液態(tài)CO₂在短時間內(nèi)迅速汽化，相變壓力在10～50 ms內(nèi)達(dá)到峰值；當(dāng)儲液管內(nèi)的壓力達(dá)到泄能片的最大臨界壓力時，高壓氣態(tài)CO₂會沖破泄能片，從泄放孔向外釋放，迅速作用于周邊的巖石，從而導(dǎo)致巖石破壞。

1.2 CO₂相變致裂巖石的沖擊壓力計算

CO₂相變致裂巖石時，液態(tài)CO₂在發(fā)熱管的作用下，快速轉(zhuǎn)化為高壓CO₂氣體，高壓脈沖氣體沖壞泄能片，進入巖石致裂孔，并壓縮孔內(nèi)空氣，然后與孔壁巖石發(fā)生碰撞，并在巖石中形成沖擊波。針對正拱普通型致裂管，剪切破壞公式可較好地反映致裂管的爆破峰值壓力^［14-15］。致裂管的爆破壓力

p_k=2σ_sδR。（1）

式中：R為泄放孔半徑；σ_s為致裂管泄能片的剪切強度，取0.6～0.8倍的屈服強度；δ為泄能片厚度。

選取實際工程中常用的51型、85型以及100型致裂管，參數(shù)如表1所示。通過式（1）可以計算得到不同型號的致裂管的爆破壓力p_k。

致裂管的爆破壓力傳遞至致裂孔表面，與孔壁巖石發(fā)生碰撞，產(chǎn)生放大效應(yīng)，并在巖石中形成沖擊波。CO₂致裂管多數(shù)為不耦合裝藥，高壓氣體碰撞孔壁激起的巖石沖擊波初始壓力

p_R=p_kd_gd₀⁶n。（2）

式中：n為CO₂高壓氣體和孔壁發(fā)生碰撞后的壓力放大倍數(shù)，受入射空氣沖擊波的壓力影響，取n=8；d_g為致裂管的外徑；d₀為致裂孔的直徑。

致裂沖擊波在巖石中傳播時不斷衰減，最后變成應(yīng)力波。巖石中任一點引起的徑向應(yīng)力和切向應(yīng)力可表示為^［15］

σ_r=p_Rr_br^α;（3）

σ_θ=-bσ_r。（4）

式中：σ_r、σ_θ分別為巖石中的徑向應(yīng)力和切向應(yīng)力；r_b為致裂孔半徑；r為巖石距致裂孔中心的距離；α為載荷傳播衰減指數(shù)，α=2±μ_d/（1-μ_d），正、負(fù)號分別對應(yīng)沖擊波和應(yīng)力波區(qū)；μ_d為巖石的動態(tài)泊松比，μ_d=0.8μ， μ是巖石的靜態(tài)泊松比；b為側(cè)向應(yīng)力系數(shù)，b=μ_d/（1-μ_d） 。

利用式（3），可得到在致裂孔孔壁處的峰值應(yīng)力。CO₂致裂管在致裂孔孔壁處的峰值應(yīng)力隨泄能片厚度變化的曲線如圖2所示。在巖石距致裂孔中心的距離r=0.1 m處測得的峰值應(yīng)力隨CO₂充裝質(zhì)量變化的曲線如圖3所示。

從圖2可知，對于同一型號的致裂管，致裂孔孔壁處的峰值應(yīng)力隨著泄能片厚度的增加呈直線上升趨勢。對于不同型號的致裂管，在同一泄能片厚度情況下，爆破壓力隨著致裂管外徑變大而逐漸增大。

圖3為在不同初始壓力的條件下，3種型號的致裂管中分別充填0.9、1.4 kg和3.5 kg的液態(tài)CO₂時，相變致裂巖石所得到的峰值壓力。圖3中，在初始壓力一定的條件下，離致裂孔中心0.1 m處，巖石受到的峰值應(yīng)力隨致裂管中CO₂充裝質(zhì)量的增加呈上升趨勢。CO₂充裝質(zhì)量從0.9 kg增加到1.4 kg時，峰值應(yīng)力上升較為明顯;而CO₂充裝質(zhì)量從1.4 kg增加到3.5 kg時的峰值應(yīng)力增加較為緩慢。

1.3 CO₂相變致裂巖石的損傷半徑

液態(tài)CO₂相變致裂巖石過程中，原始巖石將以致裂孔為中心，依次形成壓碎區(qū)、裂隙區(qū)和振動區(qū)^［12］。大量實驗證明，壓碎區(qū)和裂隙區(qū)形成了巖石損傷的有效范圍。

巖石中任一點的應(yīng)力強度可表示為

σ_i=σ_r2［（1+b）²-2μ_d（1-b）²（1-μ_d）+（1+b²）］¹²。（5）

由于該爆破所產(chǎn)生的爆炸載荷在內(nèi)部是沿軸向均勻分布的^［16］，為簡化計算，可以將該問題視為平面應(yīng)變問題，令

B＝［（1+b）²-2μ_d（1-b）²（1-μ_d）+（1+b²）］¹²。（6）

根據(jù)Mises準(zhǔn)則，對巖石的破壞可以表示為

σ_i≥σ₀；（7）

σ₀=σ_cd（壓碎區(qū)）;

σ_td（裂隙區(qū)）。（8）

式中：σ₀為巖石的單軸受力條件下的破壞強度；σ_cd為巖石的單軸動態(tài)抗壓強度；σ_td為巖石的單軸動態(tài)抗拉強度。

CO₂爆破后，若巖石中的應(yīng)力超過σ_cd，則將在巖石致裂孔孔壁周圍形成壓碎區(qū)，利用式（1）、式（2）、式（5）～式（7），可得到壓碎區(qū)半徑

R_s=Bp_R2σ_cd^1αr_b。（9）

在壓碎區(qū)之外，巖石中的應(yīng)力小于σ_cd，壓碎區(qū)范圍不擴展，但巖石受拉應(yīng)力作用產(chǎn)生張拉裂隙，形成裂隙區(qū)。根據(jù)式（5）～式（7），可得在壓碎區(qū)與裂隙區(qū)分界面上徑向應(yīng)力為

σ_R=σ_r|_r=Rs=2σ_cdB。（10）

在壓碎區(qū)之外，載荷以應(yīng)力波的形式繼續(xù)向外傳播，衰減指數(shù)為

β=2-μ_d1-μ_d。（11）

于是，整理式（5）～式（7）及式（10），可得到巖石中裂隙區(qū)半徑為

R_l=σ_cdσ_td^1βR_s。（12）

由式（2）、式（9）、式（12）可得，非耦合情況下，當(dāng)非耦合系數(shù)較小時，CO₂相變致裂巖石的損傷半徑為

R=Bnp_k2σ_cdd_gd₀^61ασ_cdσ_td^1βr_b。（13）

2 CO₂相變致裂巖石損傷的模型驗證

2.1 模型建立

采用3DEC離散元程序?qū)O₂相變致裂過程進行模擬，將致裂作用簡化為CO₂相變產(chǎn)生的沖擊壓力直接作用在巖壁上?？紤]實際有效范圍，確定模型及模型的邊界條件如圖4所示。模型長4 m、寬4 m、高0.3 m，中部致裂孔半徑為0.1 m。

在結(jié)構(gòu)面間距發(fā)生變化的條件下，需對塊體進行網(wǎng)格化，從而滿足頻率和波長控制的網(wǎng)格尺寸限制。為了精確地描述沖擊應(yīng)力波的傳播過程，通常網(wǎng)格單元的尺寸Δl應(yīng)不大于最高頻率波長的1/8～1/10。給定單元網(wǎng)格尺寸所允許的輸入波最高頻率

f_max=cλ=c10Δl。（14）

式中：λ為最高頻率所對應(yīng)的波長；Δl為網(wǎng)格單元的最大長度，取2 cm；c為波在巖石中的傳播速度，取縱波速度和橫波速度中的較小者。

對模型采用犀牛軟件劃分塊體，大多數(shù)塊體劃分尺寸范圍為0.5～1.5 cm，整個模型共有12 569個節(jié)點和58 380個塊體。為了減少動載荷作用下邊界反射波的影響，采用黏性邊界。

數(shù)值模擬中，巖石介質(zhì)采用理想線彈性模型，節(jié)理采用庫倫滑移模型。對巖石節(jié)理的法向剛度k_n和剪切剛度k_s進行估算^［17-18］：

k_n=k_s=10 maxK+4G/3Δz_min。（15）

式中：K為體積模量；G為剪切模量；Δz_min為接觸面法向方向上連接區(qū)域的最小尺寸。

考慮到?jīng)_擊載荷的動態(tài)硬化效應(yīng)，巖石的動彈性模量取靜態(tài)模量的1.2倍，其余參數(shù)與靜態(tài)參數(shù)相同，巖石的力學(xué)特性參數(shù)見表2～表3。

2.2 CO₂相變沖擊載荷

結(jié)合文獻(xiàn)［19］的相關(guān)研究，CO₂壓力曲線接近于高斯曲線，壓力從0到峰值、再降為0的時間變化范圍約為30 ms。故CO₂相變致裂巖石過程中，高壓氣體對孔壁內(nèi)側(cè)的初始沖擊壓力p_R的時間衰減函數(shù)為

p（t）=p_Rae^-［b（t-t_m^）］22。（16）

式中：a、b為致裂衰減系數(shù)，與泄能片厚度有關(guān)；t_m為峰值載荷出現(xiàn)的時間。

CO₂相變致裂壓力從劇增點到峰值的持續(xù)時間約為10～15 ms。取增壓持續(xù)時間t_m為15 ms。

2.3 數(shù)值模型驗證

2.3.1 CO₂相變致裂巖石應(yīng)力場對比

選用特征最為明顯的85型致裂管，在泄能片厚度為4.19 mm的情況下進行CO₂相變沖擊載荷作用于巖石的損傷數(shù)值模擬。在模型中部距邊界1.9、1.8、1.4、1.1、0.8、0.5、0.2 m和0 m位置處布置測點，巖石中部不同測點處的水平應(yīng)力隨時間的變化如圖5所示。根據(jù)不同位置與時間上水平應(yīng)力的分布特征，可以看出：CO₂相變產(chǎn)生的沖擊壓力直接作用于致裂孔孔壁巖石上，孔壁周圍巖石所受水平峰值應(yīng)力最大；隨著致裂時間的增加，巖石所受水平應(yīng)力呈先增加、后衰減的趨勢。相變產(chǎn)生的沖擊波不斷向遠(yuǎn)處傳播，由于巖石的阻尼作用導(dǎo)致沖擊波能量消耗變?yōu)閼?yīng)力波，在巖石中產(chǎn)生的應(yīng)力逐漸減弱。距離致裂孔越遠(yuǎn)，巖石所受的峰值應(yīng)力越小；同時，升壓時間也不斷滯后。

利用式（2）、式（3）得到與數(shù)值模擬測點相同位置處的理論峰值應(yīng)力，與數(shù)值模擬得到的結(jié)果進行比較。如圖6所示，理論峰值應(yīng)力曲線和數(shù)值計算結(jié)果較為接近，近域（0.1 m）處的計算誤差為10%，遠(yuǎn)域（2.0 m）處的計算誤差為21%。該模型可以較好地模擬CO₂相變沖擊應(yīng)力場的衰減規(guī)律。

2.3.2 CO₂相變致裂巖石的損傷半徑對比

泄能片厚度δ分別為3.36、4.19、5.03、5.87 mm和6.71 mm時，進行85型致裂管致裂巖石的損傷半徑的數(shù)值計算，獲得不同泄能片厚度下巖石損傷范圍的模擬結(jié)果，如圖7所示。

由圖7可知，不同泄能片厚度下的致裂損傷半徑分別約為0.53、0.59、0.71、0.79 m和0.86 m。可見，泄能片厚度越大，產(chǎn)生的沖擊壓力越大，巖石致裂損傷半徑也逐漸增大。同理，由式（13）獲得理論致裂損傷半徑分別為0.59、0.66、0.71、0.76 m和0.81 m。將兩者進行比較，可以得到圖8。

由圖8可以看出，隨著泄能片厚度的增加，致裂損傷半徑逐漸增大，數(shù)值模擬與理論計算大致相近，最大誤差10%。原因在于：液態(tài)CO₂相變所產(chǎn)生的高壓氣體剪壞泄能片后，形成沖擊波沖擊巖石，致使出氣孔周圍產(chǎn)生主裂紋；后續(xù)高壓氣體嵌入由主裂紋所衍生的裂縫中；裂縫繼續(xù)延伸擴展，直至高壓氣體全部耗散。沖擊波能夠引起巖石產(chǎn)生徑向裂紋，而致裂管泄能片的厚度與產(chǎn)生沖擊波的大小成正比，這將促進巖石裂紋的產(chǎn)生，增大裂紋的擴展范圍。因此，隨著泄能片厚度的增大，巖石損傷范圍亦逐漸增大。可為實際施工過程中合理布置巖石的致裂孔間距和位置提供理論依據(jù)，以達(dá)到最佳的致裂效果。

3 致裂參數(shù)對巖石損傷半徑的影響

根據(jù)上述模擬結(jié)果可知，理論計算與數(shù)值模擬的計算結(jié)果接近?？梢哉J(rèn)為，基于Mises準(zhǔn)則，利用損傷半徑來描述巖石損傷的壓碎區(qū)和裂隙區(qū)是合理的。因此，選取巖石損傷半徑作為母因素，子因素分別為泄能片厚度、巖石強度、致裂孔孔徑、CO₂充裝質(zhì)量。通過理論計算公式分析這4個因素與損傷半徑之間的關(guān)系。

給定不同的泄能片厚度、巖石強度和致裂孔孔徑，計算85型致裂管爆破在巖石內(nèi)產(chǎn)生的損傷半徑，如圖9～圖11所示。

3.1 泄能片厚度

從圖9可以看出，在不同巖石強度和致裂孔孔徑下，損傷半徑均隨著泄能片厚度的增加而不斷增加，且在較低強度和較小孔徑下增長幅度較為明顯。泄能片厚度從3 mm增加到6 mm，當(dāng)巖石強度為20 MPa時，損傷半徑增大0.38 m，增幅為36%；當(dāng)致裂孔孔徑為90 mm時，損傷半徑增大0.25 m，增幅為38%?？芍?，當(dāng)泄能片厚度增加一倍時，損傷半徑的最大增幅可達(dá)38%左右。

3.2 巖石強度

從圖10可以看出，隨著巖石強度的增大，損傷半徑逐漸減小。巖石強度越小，受致裂孔孔徑和泄能片厚度的影響越大；且隨著巖石強度的增大，這種影響越來越小。巖石強度從20 MPa增加到40 MPa，泄能片厚度為2 mm時，巖石損傷半徑減小0.25 m，減幅為29%；致裂孔孔徑為90 mm時，巖石損傷半徑減小0.40 m，減幅為28%?？梢姡瑤r石強度增大一倍，損傷半徑的最大減幅為29%。

3.3 致裂孔孔徑

由圖11可知，損傷半徑隨致裂孔孔徑的增大而迅速減小，且呈指數(shù)衰減。當(dāng)致裂孔孔徑增加為致裂管外徑的2倍、泄能片厚度為4 mm時，損傷半徑減小到0.3 m，減幅為68%；巖石強度為20 MPa時，損傷半徑減小到0.45 m，減幅為66%。可以看出，致裂孔孔徑增大1倍時，損傷半徑的最大減幅為68%。

3.4 CO₂充裝質(zhì)量

選取表1中3種不同CO₂充裝質(zhì)量的致裂管，計算在不同泄能片厚度、不同巖石強度下致裂巖石所產(chǎn)生的損傷半徑，如圖12所示。

從圖12可以看出，隨著CO₂充裝質(zhì)量的增加，損傷半徑的增加幅度逐漸減小，且呈對數(shù)增長。當(dāng)CO₂充裝質(zhì)量從0.9 kg增長到2.0 kg，泄能片厚度為6 mm時，損傷半徑增長0.27 m，增幅為39%；巖石強度為20 MPa時，損傷半徑增長0.31 m，增幅為40%。當(dāng)CO₂充裝質(zhì)量超過2.0 kg，損傷半徑幾乎不再增長?？梢钥闯?，當(dāng)充裝的CO₂達(dá)到致裂管的容積之后，再增加CO₂質(zhì)量，對損傷半徑的影響微乎其微。

從以上分析可以看出，泄能片厚度和CO₂充裝質(zhì)量的增加在一定范圍內(nèi)能有效增大損傷半徑，而巖石強度和致裂孔孔徑的增大會減小損傷半徑。且當(dāng)致裂管型號一定、CO₂充裝質(zhì)量達(dá)到致裂管的最大充裝量時，增加泄能片厚度會顯著增加巖體的損傷半徑，此時的損傷半徑增長率最大；同樣的，增大致裂孔孔徑會使巖石損傷半徑迅速減小。

4 現(xiàn)場實驗分析

為更好地分析CO₂相變致裂巖石的損傷范圍，在重慶南桐綠色礦山進行超臨界CO₂相變致裂巖石現(xiàn)場實驗，如圖13所示。采用RSM-SY6樁基聲波探測儀進行聲波測試。依據(jù)SL—94《水工建筑物巖石基礎(chǔ)開挖工程施工技術(shù)規(guī)范》規(guī)定，以聲波變化率η＞10%作為判定巖石損傷的標(biāo)準(zhǔn)，可得出損傷門檻D_cr=0.19。當(dāng)損傷變量超過損傷門檻時，巖石進入破壞狀態(tài)。

根據(jù)圖14可以看出，所建立的理論模型比現(xiàn)場實驗測得的數(shù)值稍大，但規(guī)律性相同。這是因為，理論分析中假定巖石為連續(xù)、均值的材料，但實際工程中的巖石存在節(jié)理裂隙以及結(jié)構(gòu)面等形態(tài)。

現(xiàn)場實驗時，超臨界CO₂相變致裂巖石的過程中無明顯飛石，致裂孔有部分灰塵噴出，致裂孔外產(chǎn)生大量裂紋，引起的周邊巖石振動很小。致裂引起的巖塊形成率較高，顆粒級配適中，可以有效提高巖石開挖規(guī)格石料的獲得率。

5 結(jié)論

通過對CO₂相變致裂巖石的損傷范圍的理論分析和基于3DEC軟件的巖石致裂損傷的數(shù)值模擬研究，得出以下結(jié)論：

1）通過CO₂相變致裂巖石損傷范圍理論分析，給出了CO₂相變產(chǎn)生的沖擊壓力的計算方法，并根據(jù)Mises準(zhǔn)則和損傷力學(xué)理論，建立了CO₂相變致裂巖石的損傷半徑的計算模型。

2）采用數(shù)值模擬分析了CO₂相變產(chǎn)生的沖擊波在巖石中的衰減規(guī)律，得到了不同泄能片厚度時巖石的致裂損傷半徑。結(jié)果表明，對于85型致裂管，在泄能片的厚度為3.36～6.71 mm時，產(chǎn)生的沖擊波大小與泄能片厚度成正相關(guān)；隨著泄能片厚度的增加，產(chǎn)生的沖擊壓力增大，引起的巖石致裂損傷半徑也越大。

3）泄能片厚度和CO₂充裝質(zhì)量的增加能有效增大損傷半徑，巖石強度和致裂孔孔徑的增大會減小損傷半徑。且增大泄能片厚度是增大損傷半徑最有效的手段。巖石的損傷范圍隨著致裂孔孔徑的增大呈指數(shù)衰減。

參考文獻(xiàn)

［1］周盛濤， 羅學(xué)東， 蔣楠， 等. 二氧化碳相變致裂技術(shù)研究進展與展望［J］. 工程科學(xué)學(xué)報， 2021， 43（7）： 883-893.

ZHOU S T， LUO X D， JIANG N， et al. A review on fracturing technique with carbon dioxide phase transition［J］. Chinese Journal of Engineering， 2021， 43（7）： 883-893.

［2］周西華， 門金龍， 宋東平， 等. 液態(tài)CO₂爆破煤層增透最優(yōu)鉆孔參數(shù)研究［J］. 巖石力學(xué)與工程學(xué)報， 2016， 35（3）： 524-529.

ZHOU X H， MEN J L， SONG D P， et al. Research on optimal borehole parameters of antireflection in coal seam by liquid CO₂ blasting［J］. Chinese Journal of Rock Mechanics and Engineering， 2016， 35（3）： 524-529.

［3］孫小明.液態(tài)二氧化碳相變致裂掏槽破巖試驗研究［J］.煤炭科學(xué)技術(shù)， 2021， 49（8）： 81-87.

SUN X M. Experimental study on cutting and rock brea-king by liquid CO₂ phase transition fracturing technology［J］. Coal Science and Technology， 2021， 49（8）：81-87.

［4］孫可明， 辛利偉， 吳迪. 超臨界CO₂氣爆煤體致裂機理實驗研究［J］.爆炸與沖擊， 2018， 38（2）： 302-308.

SUN K M， XIN L W， WU D． Experimental study on fracture mechanism of coal caused by supercritical CO₂ explosion ［J］． Explosion and Shock Waves， 2018， 38（2）： 302-308.

［5］孫可明， 辛利偉， 吳迪， 等. 初應(yīng)力條件下超臨界CO₂氣爆致裂規(guī)律研究［J］.固體力學(xué)學(xué)報， 2017， 38（5）： 473-482.

SUN K M， XIN L W， WU D， et al． Mechanism of fracture caused by supercritical CO₂ explosion under the impact of initial stress［J］． Chinese Journal of Solid Mechanics， 2017， 38（5）： 473-482.

［6］WANG H， CHENG Z， ZOU Q， et al. Elimination of coal and gas outburst risk of an outburst-prone coal seam using controllable liquid CO₂ phase transition fracturing［J］. Fuel， 2021， 284： 119091.

［7］王杰， 張云鵬， 葛曉東， 等.液態(tài)CO₂相變破巖參數(shù)及數(shù)值模擬研究［J］.爆破器材， 2022， 51（4）： 16-22.

WANG J， ZHANG Y P， GE X D， et al. Study on rock breaking parameters and numerical simulation of phase transition of liquid CO₂ ［J］. Explosive Materials， 2022， 51（4）： 16-22.

［8］郭云龍， 李紀(jì)寶， 孫崔源， 等. 孔間距對二氧化碳相變爆破效果的影響模擬分析［J］. 鐵道建筑， 2022， 62（6）： 125-129.

GUO Y L， LI J B， SUN C Y， et al. Simulation analysis of influence of borehole spacing on effect of CO₂ phase change blasting［J］. Railway Engineering， 2022， 62（6）： 125-129.

［9］陳少波， 盧玉斌， 袁海梁， 等. 基于ANSYS/LS-DYNA的液態(tài)CO₂相變破巖裂紋擴展特征研究［J］.中國安全生產(chǎn)科學(xué)技術(shù)， 2023， 19（6）： 74-82.

CHEN S B， LU Y B， YUAN H L， et al. Study on crack propagation characteristics of liquid CO₂ phase transition rock breaking based on ANSYS/LS-DYNA ［J］. China Production Safety Science and Technology， 2023， 19（6）： 74-82.

［10］ZHANG W， XU K， LEI Y， et al. Evolutionary features in damage and destruction of gas-rich coal seam by CO₂ phase-transition blasting ［J］. Ekoloji， 2018， 27（106）： 1605-1613.

［11］郭志興. 液態(tài)二氧化碳爆破筒及現(xiàn)場試爆［J］.爆破， 1994（3）： 72-74.

［12］田澤礎(chǔ). 液態(tài)二氧化碳相變致裂裂縫形態(tài)及影響因素研究［D］. 徐州： 中國礦業(yè)大學(xué)， 2018.

TIAN Z C. Crack form and influencing factors of liquid carbon dioxide phase transition fracturing ［D］. Xuzhou： China University of Mining and Technology， 2018.

［13］程選生， 劉華東， 張尚龍， 等.巖質(zhì)路塹液態(tài)CO₂相變爆破充裝量的試驗研究［J］.土木工程學(xué)報， 2023， 56（7）： 137-146.

CHENG X S， LIU H D， ZHANG S L， et al. Test on filling of liquid CO₂ phase transition blasting in rock cutting ［J］. China Civil Engineering Journal， 2023， 56 （7）： 137-146.

［14］曾永慶. 二氧化碳致裂爆破振動效應(yīng)研究［D］. 武漢： 中國科學(xué)院大學(xué)， 2019.

ZENG Y Q. Study on blasting vibration effect for carbon dioxide blasting ［D］. Wuhan： University of Chinese Academy of Sciences， 2019.

［15］王文龍. 鉆眼爆破［M］. 北京： 煤炭工業(yè)出版社， 1984.

［16］冷振東， 盧文波， 陳明， 等.巖石鉆孔爆破粉碎區(qū)計算模型的改進［J］.爆炸與沖擊， 2015， 35（1）： 101-107.

LENG Z D， LU W B， CHEN M， et al. Improvement calculation model for the size of crushed zone around blasthole ［J］. Explosion and Shock Waves， 2015， 35（1）： 101-107.

［17］HART R D. An introduction to distinct element mode-ling for rock engineering ［J］. Analysis and Design Methods， 1993： 245-261.

［18］DENG X F， CHEN S G， ZHU J B， et al. UDEC-AUTODYN hybrid modeling of a large-scale underground explosion test ［J］. Rock Mechanics and Rock Enginee-ring， 2015， 48（2）： 737-747.

［19］ZHU W C， GAI D， WEI C H， et al. High-pressure air blasting experiments on concrete and implications for enhanced coal gas drainage［J］. Journal of Natural Gas Science and Engineering， 2016， 36： 1253-1263.