林英松，姚勁松，劉 瑩，喬繼延，丁雁生

（1. 中國(guó)石油大學(xué)（華東）石油工程學(xué)院，山東 青島 266555；2. 中國(guó)科學(xué)院力學(xué)研究所，北京 100081）

近年來，頁(yè)巖氣的開發(fā)日益成為各國(guó)關(guān)注的焦點(diǎn)[1]。而頁(yè)巖氣開采的難點(diǎn)在于儲(chǔ)層具有低滲透率和低孔隙度的屬性，為產(chǎn)生工業(yè)油氣流通常需要采取水力壓裂等增產(chǎn)措施[2]。壓裂作業(yè)通常是在射孔作業(yè)基礎(chǔ)上進(jìn)行的，因此射孔孔眼處巖石的損傷，對(duì)壓裂作業(yè)有不可忽視的影響[3]。

關(guān)于射孔周圍巖石損傷情況，僅有少數(shù)學(xué)者研究了對(duì)砂巖靶射孔后孔道周圍的損傷量、應(yīng)力分布及滲透率變化等，胡柳青等[4]對(duì)不同沖擊載荷下裂紋響應(yīng)進(jìn)行數(shù)值模擬，得到了一系列隨時(shí)間變化的動(dòng)態(tài)應(yīng)力場(chǎng)及動(dòng)態(tài)應(yīng)變場(chǎng)。朱秀星等[5]、薛世峰等[6]利用LS-DYNA軟件模擬了射孔后砂巖骨架應(yīng)力、塑性應(yīng)變等參數(shù)變化。單清林等[7]建立了螺旋射孔損傷模型以研究巖體破壞后彈性參數(shù)、強(qiáng)度及滲透率變化規(guī)律。王成等[8]模擬了聚能射孔彈各參數(shù)變化時(shí)對(duì)漏斗坑直徑及深度、侵徹孔直徑及穿深等的影響。Nabipour等[9]對(duì)射孔后的砂巖進(jìn)行研究，發(fā)現(xiàn)巖石原始孔隙度越大，射孔損傷范圍越大。而針對(duì)射孔周圍巖石的微裂紋分布規(guī)律的研究較為缺乏。

本文采用FEPG有限元程序生成系統(tǒng)，根據(jù)射孔彈作用機(jī)理，考慮射孔作業(yè)時(shí)巖石上的載荷分布特征，建立模型進(jìn)行計(jì)算?？紤]巖石非均質(zhì)性，將材料細(xì)觀強(qiáng)度等參數(shù)設(shè)置為空間內(nèi)隨機(jī)分布，使其在載荷作用下能產(chǎn)生隨機(jī)裂紋。模擬計(jì)算后，將結(jié)果與物理實(shí)驗(yàn)結(jié)果進(jìn)行對(duì)比，以驗(yàn)證力學(xué)模型有效性，并分析不同條件下射孔孔道周圍巖石裂紋形態(tài)分布規(guī)律和寬度分布規(guī)律。

1 力學(xué)模型

射孔彈金屬射流侵徹巖石過程主要分3個(gè)階段[10]。

第一階段，射孔開始。射流以6～7 km/s的速度撞擊巖石，該速度高于巖石中的聲速，射流碰撞靶材時(shí)，在撞擊點(diǎn)處產(chǎn)生強(qiáng)沖擊波并向靶材內(nèi)部傳入，造成波后的巖石損傷，形成失效波。與此同時(shí)，高速射流在巖石表面砸出一個(gè)漏斗形坑，尺度比射流直徑大幾倍，僅占總孔深的很小一部分[11-12]。沖擊波、失效波與開坑，構(gòu)成射流侵徹的第一個(gè)階段。

第二階段，后續(xù)射流繼續(xù)在接觸面產(chǎn)生極高壓力，約為200 GPa，巖石產(chǎn)生塑性變形，侵徹孔持續(xù)加深，直徑為射流幾倍。此階段持續(xù)時(shí)間長(zhǎng)，對(duì)穿深貢獻(xiàn)大，因與波動(dòng)效應(yīng)無關(guān)，稱為準(zhǔn)定常階段[10]。

第三階段為終止階段，該階段射流速度低，碰撞點(diǎn)壓力大幅下降，逐漸趨近巖石強(qiáng)度，孔深停止增加。射流在孔底堆積，并在孔底前方產(chǎn)生壓實(shí)區(qū)，遠(yuǎn)處產(chǎn)生開裂區(qū)。

上述3個(gè)階段是通常的說法，針對(duì)的是孔深和孔徑。射孔周圍巖石的拉伸裂紋，作為壓裂問題的初始條件，常常不被注意，卻是本文關(guān)注的問題。

1.1 基本假設(shè)

射孔過程較為復(fù)雜，建模前需簡(jiǎn)化。

對(duì)于射流侵徹過程，第二階段持續(xù)時(shí)間最長(zhǎng)，射孔深度最深，產(chǎn)生裂紋的區(qū)域最大，因此選取此階段產(chǎn)生的侵徹孔鄰域內(nèi)壓實(shí)和遠(yuǎn)處的開裂作為研究對(duì)象。劃分網(wǎng)格時(shí)單元尺度應(yīng)與裂紋寬度保持一致，若進(jìn)行三維模擬，單元和節(jié)點(diǎn)數(shù)量太大。為減少計(jì)算量，選擇射孔圍巖某一切片作為數(shù)值模擬對(duì)象，從而將侵徹過程簡(jiǎn)化為二維問題。由于二維情形射流開孔時(shí)中心點(diǎn)處應(yīng)變率無窮大，較難模擬，初始時(shí)刻在巖石切片中心預(yù)制一個(gè)內(nèi)徑極小的孔，將開孔簡(jiǎn)化為擴(kuò)孔。最后，為保證模擬過程接近實(shí)際，需考慮巖石非均勻性。本模型中，為保留裂縫形態(tài)和方向的隨機(jī)性，不宜使用預(yù)制裂縫法，而將巖石視為細(xì)觀非均勻材料，巖石的力學(xué)參數(shù)在空間中隨機(jī)分布，較弱點(diǎn)會(huì)率先發(fā)生破壞，從而隨機(jī)產(chǎn)生裂紋。

綜上所述，做出如下假設(shè)：

(1)巖石切片無沿射孔軸向的位移，只有徑向和周向運(yùn)動(dòng)；

(2)忽略沖擊波及失效波對(duì)巖石損傷的貢獻(xiàn)；

(3)在巖石切片中心有一個(gè)小孔，按照給定徑向速度脈沖擴(kuò)孔，使孔周圍巖石破壞；

(4)巖石為細(xì)觀非均勻彈脆性材料，其彈性模量、抗剪強(qiáng)度和抗拉強(qiáng)度等在空間隨機(jī)分布，服從韋布爾分布。

綜上所述，將三維射孔過程簡(jiǎn)化為平面應(yīng)變問題；忽略沖擊波效應(yīng)的運(yùn)動(dòng)方程退化為平衡方程；細(xì)觀非均勻彈脆性材料的宏觀力學(xué)行為不限于彈脆性。

1.2 本構(gòu)關(guān)系、基本方程的弱形式和邊界條件

為便于描寫隨機(jī)裂紋，使用直角坐標(biāo)系。細(xì)觀彈性本構(gòu)關(guān)系為[13]：

式中：E為彈性模量，μ為泊松比。將式(1)簡(jiǎn)記為：

式中：σ為應(yīng)力張量，D為彈性矩陣，ε為應(yīng)變張量。

FEPG (finite element program generator)平臺(tái)適用許多偏微分方程描寫的問題，但需要給出偏微分方程的弱形式。針對(duì)該問題，運(yùn)用虛位移原理,由平衡方程得到虛功表達(dá)式：

式中：δu、δv分別為兩個(gè)方向上的虛位移，fx、fy分別為兩個(gè)方向的質(zhì)量力分量，S為巖石切片面積。利用分部積分公式將上式化為：

式中：Γ為巖石切片的邊界，Tx、Ty分別為邊界上兩個(gè)方向的作用力。式(4)改寫成向量形式：

式中：εT、uT分別為應(yīng)力邊界上的應(yīng)變與位移，f、T分別為質(zhì)量力矩陣和邊界力矩陣。代入本構(gòu)方程式(2)，得：

由外邊界上的邊界條件，可得：

式中：上標(biāo)中的o表示外邊界。這就是本問題的弱形式。本模型中預(yù)制孔內(nèi)徑為1.5 mm，內(nèi)邊界給定位移。內(nèi)邊界Γi的虛位移為零，使式(7)中內(nèi)邊界項(xiàng)為零。邊界條件表示為：

將隨侵徹過程變化的擴(kuò)孔速度簡(jiǎn)化為兩個(gè)階段，在起始時(shí)間步位移量很大，后續(xù)時(shí)間步位移量很小，如：

式中：ur（t）為第t步時(shí)內(nèi)邊界位移量，Δur為整個(gè)擴(kuò)孔過程中內(nèi)邊界總位移量。式(2)、(7)和(8)～(9)等，就是用FEPG軟件平臺(tái)生成本問題的有限元程序所需要的數(shù)學(xué)表達(dá)式。

1.3 細(xì)觀破壞準(zhǔn)則和細(xì)觀非均勻性

巖石的細(xì)觀壓剪開裂遵循Mohr-Coulomb破壞準(zhǔn)則，細(xì)觀拉伸開裂遵循拉伸破壞準(zhǔn)則。

巖石一旦開裂，基于連續(xù)介質(zhì)建立的偏微分方程便不再成立，因此進(jìn)行有限元計(jì)算時(shí)，需對(duì)裂紋進(jìn)行處理。本模型中使用“模量折減法”進(jìn)行處理。當(dāng)單元達(dá)到破壞條件時(shí)，仍保持材料連續(xù)，折減單元模量以擴(kuò)大其變形，從而反映宏觀斷裂。

經(jīng)過分析[14]，孔眼附近巖石將發(fā)生壓剪破壞。單元一旦發(fā)生剪切破壞，該單元的剪切彈性模量衰減到原值的若干分之一?？籽圻h(yuǎn)處巖石發(fā)生拉伸破壞，縫寬大于前者，因此破壞后單元的拉伸彈性模量衰減量大于前者的衰減量。模量衰減若干倍，即相同應(yīng)力下單元變形擴(kuò)大若干倍。模量的衰減量在計(jì)算中根據(jù)擬合圖像確定。

為了能夠使模擬結(jié)果產(chǎn)生近似真實(shí)的隨機(jī)裂紋，令彈性模量、剪切強(qiáng)度和拉伸強(qiáng)度等細(xì)觀參數(shù)在空間中隨機(jī)分布，分布模型選擇韋布爾分布。

為方便計(jì)算，使用一系列滿足均勻分布的偽隨機(jī)數(shù)通過變化得到滿足韋布爾分布的隨機(jī)數(shù)列A[15]：

式中：F為區(qū)間內(nèi)均勻分布的偽隨機(jī)數(shù)列，ηA和mA是關(guān)于韋布爾分布的兩個(gè)參數(shù)。則巖石的彈性模量為：

同理設(shè)置巖石抗拉強(qiáng)度與抗剪強(qiáng)度在切片中隨機(jī)分布。

1.4 裂縫寬度

拉伸破壞開裂后，需對(duì)縫寬進(jìn)行計(jì)算。由于模型中巖石切片只受徑向壓力作用，因此環(huán)向應(yīng)變?yōu)榈谝恢鲬?yīng)變，縫寬計(jì)算公式為：

式中：h為裂縫寬度，ε1為第一主應(yīng)變，c為縫寬系數(shù)。

2 數(shù)值模擬

2.1 幾何模型

侵徹巖石切片網(wǎng)格劃分如圖1所示。

圖1 幾何模型示意圖Fig. 1 Schematic diagram of the geometric model

2.2 參數(shù)設(shè)置

根據(jù)現(xiàn)場(chǎng)巖心資料及射孔資料，設(shè)置模型參數(shù)如表1。除此之外，為方便計(jì)算，韋布爾分布均勻度參數(shù)mA均取3，圍壓取1 MPa。

表1 模擬計(jì)算參數(shù)Table 1 Parameters of simulation calculation

2.3 計(jì)算結(jié)果

數(shù)值模擬計(jì)算結(jié)果如圖2所示。

圖2 數(shù)值模擬結(jié)果圖Fig. 2 Numerical simulation result

2.4 損傷規(guī)律分析

根據(jù)破壞形態(tài)，由內(nèi)而外可將切片分為四塊區(qū)域：首先是近孔眼處的壓剪破壞區(qū)，該區(qū)域的破壞幾乎全為壓剪破壞，分析認(rèn)為其原因是射孔時(shí)中心孔瞬間產(chǎn)生大變形，該區(qū)域首先達(dá)到壓剪破壞條件，產(chǎn)生環(huán)帶狀壓剪破壞區(qū)。第二層為拉伸破壞集中區(qū)，該區(qū)域單元徑向壓力不足以產(chǎn)生壓剪破壞，但能導(dǎo)致足以產(chǎn)生拉伸破壞的環(huán)向拉力。第三層為拉伸裂紋擴(kuò)展區(qū)，該區(qū)域主要是初期大變形后，在準(zhǔn)靜態(tài)載荷下稍遠(yuǎn)處單元逐漸拉伸破壞而形成的。最外層為未破壞區(qū)，隨著距孔眼距離逐漸增大，單元應(yīng)力逐漸降低，不足以產(chǎn)生破壞。切片破壞區(qū)域劃分如圖3所示。

圖3 巖石破壞分區(qū)圖Fig. 3 The distribution of rock damage zones

記切片上半徑與侵徹孔直徑之比r/d0為無量綱半徑，統(tǒng)計(jì)距射孔中心不同無量綱半徑處裂紋數(shù)量和縫寬，做裂紋數(shù)量分布圖和縫寬分布圖如圖4～5所示。

2.4.1 不同載荷下裂縫分布規(guī)律

使用不同型號(hào)射孔彈對(duì)巖石造成的破壞不同。由于本模型中，中心孔邊界條件為應(yīng)變邊界條件，故不同射孔彈作用載荷不同，反映在該模型中影響的是侵徹孔直徑。保持其他參數(shù)不變（圍壓為2 MPa），分別設(shè)置侵徹孔直徑d0為7、9、11、13、15 mm，分析裂紋分布的變化規(guī)律。圖6為其中侵徹孔直徑為7 mm時(shí)切片破壞形態(tài)：

圖4 數(shù)值模擬裂紋數(shù)分布圖Fig. 4 Crack number distribution of numerical simulation

圖5 數(shù)值模擬不同縫寬裂紋數(shù)分布圖Fig. 5 Simulated number distribution of different width crack

對(duì)比圖2，分析可知隨侵徹孔直徑增大，切片壓剪破壞區(qū)、拉伸破壞集中區(qū)、拉伸破壞擴(kuò)展區(qū)都有不同程度的增大。其中拉伸破壞擴(kuò)展區(qū)增大迅速，即對(duì)載荷大小最敏感，侵徹孔直徑為13 mm和15 mm時(shí)拉伸破壞擴(kuò)展區(qū)已擴(kuò)大到切片邊緣。裂紋數(shù)目分布如圖7所示。

圖6 侵徹孔直徑7 mm圍壓2 MPa時(shí)巖石切片破壞形態(tài)Fig. 6 The rock damage form with 7 mm perforating charge under 2 MPa confining pressure

圖7 不同侵徹孔直徑下巖石切片裂紋分布圖Fig. 7 Distribution of rock crack number with perforating charge of different diameters

記壓剪破壞區(qū)、拉伸破壞集中區(qū)和拉伸破壞擴(kuò)展區(qū)外邊界半徑分別為r1、r2和r3，如圖3所示。根據(jù)模擬結(jié)果，歸納侵徹孔直徑對(duì)巖石切片裂紋數(shù)目分布在圍壓不變時(shí)的影響如下：

(1)不同侵徹孔直徑d0條件下，裂紋數(shù)目分布形態(tài)類似。壓剪破壞區(qū)邊界無量綱半徑r1/d0≈1.5；無邊界影響時(shí)拉伸破壞擴(kuò)展區(qū)無量綱半徑r3/d0＜15。

(2)裂紋數(shù)目峰值半徑落在C1≤r/d0≤C2范圍，稱C2-C1為峰值帶寬。侵徹孔孔徑越大，裂紋數(shù)目越多，峰值帶寬越大。切片半徑記作rB。當(dāng)rB/d0＞20時(shí)，峰值帶寬C2-C1很?。划?dāng)rB/d0≈18時(shí)，峰值帶寬上升到稱C2-C1≈2；當(dāng)rB/d0≈15時(shí)，邊界開始出現(xiàn)裂紋；當(dāng)rB/d0≈14時(shí)，邊界出現(xiàn)大量裂紋，峰值上界C2迅速增長(zhǎng)。

(3)隨侵徹孔直徑增大，峰值前裂紋數(shù)目增加速度及峰值后裂紋數(shù)目衰減速度加快。

(4)侵徹孔直徑較大，rB/d0接近15和小于15時(shí)，峰值后裂紋數(shù)目受邊界效應(yīng)的影響出現(xiàn)波動(dòng)。

2.4.2 不同圍壓下裂縫分布規(guī)律

保持其他參數(shù)不變（侵徹孔直徑為11 mm），分別將模擬圍壓設(shè)為0、0.5、1.0、1.5、2.0 MPa，分析裂紋變化規(guī)律。圖8為其中零圍壓下巖石切片破壞情況。

對(duì)比圖2，總結(jié)模擬結(jié)果可得，隨著圍壓增大，壓剪破壞區(qū)范圍基本不變，拉伸破壞集中區(qū)和拉伸破壞擴(kuò)展區(qū)范圍有不同程度的減小，其中前者減小緩慢，后者對(duì)圍壓變化較敏感。統(tǒng)計(jì)5種圍壓下裂縫數(shù)目分布如圖9所示。

圖8 侵徹孔徑9 mm圍壓為0時(shí)巖石切片裂縫形態(tài)Fig. 8 The rock damage form with 9 mm charge under no confining pressure

圖9 不同圍壓下巖石切片裂紋分布圖Fig. 9 The distribution of rock crack number under different confining pressures

歸納圍壓對(duì)巖石切片裂紋分布規(guī)律影響如下：

(1)圍壓變化基本不影響壓剪破壞區(qū)和拉伸破壞集中區(qū)，峰值前裂紋數(shù)目的增長(zhǎng)幾乎相同；

(2)圍壓減小時(shí)，裂紋數(shù)目的峰值有所增大，峰值帶寬C2-C1逐漸變大，拉伸破壞擴(kuò)展區(qū)無量綱半徑逐漸增長(zhǎng)，直至達(dá)到切片邊緣；

(3)圍壓越小，裂紋數(shù)目越多，且峰值后裂紋數(shù)目衰減得越慢；

(4)圍壓較小時(shí)，峰值后裂紋數(shù)目受邊界效應(yīng)的影響出現(xiàn)波動(dòng)。

3 物理模擬實(shí)驗(yàn)

3.1 實(shí)驗(yàn)方法

因天然巖心力學(xué)性質(zhì)難以重復(fù)，采用水泥石代替天然巖石進(jìn)行模擬實(shí)驗(yàn)。為保證實(shí)驗(yàn)可靠性，根據(jù)天然巖心的滲透率、孔隙度等參數(shù)選擇性能相近的水泥型號(hào)。經(jīng)對(duì)比選擇勝維G級(jí)水泥（水灰比0.44）作為試樣，兩者參數(shù)對(duì)比如表2。

經(jīng)實(shí)驗(yàn)測(cè)定，水泥試樣抗拉強(qiáng)度為1.7 MPa，抗壓強(qiáng)度為36 MPa。

使用勝利油田測(cè)井公司高壓射孔設(shè)備，按照API RP 19B標(biāo)準(zhǔn)[16]流程進(jìn)行實(shí)驗(yàn)。實(shí)驗(yàn)裝置如圖10所示。由于實(shí)驗(yàn)設(shè)備尺寸受限，水泥試樣直徑最大為200 mm，為盡量避免邊界效應(yīng)，試樣直徑應(yīng)盡可能大[17]，因此試樣直徑設(shè)為200 mm。為觀察橫剖面裂紋分布，采用預(yù)制剖面法，將3段分別長(zhǎng)200、300和500 mm的試樣堆疊放置，如圖11所示。

表2 天然頁(yè)巖與水泥試樣數(shù)據(jù)對(duì)比表Table 2 Comparison of performance parameters betweennatural shale and cement sample

圖10 實(shí)驗(yàn)裝置示意圖Fig. 10 Schematic diagram of experimental device

此外，為方便后續(xù)工作，將試樣預(yù)制剖面編號(hào)，各剖面編號(hào)名稱見表3。

圖11 水泥試樣安裝圖Fig. 11 Diagram of cement sample installation

表3 水泥試樣尺寸及編號(hào)表Table 3 Size and number of cement samples

3.2 物理實(shí)驗(yàn)結(jié)果

實(shí)驗(yàn)后[17]發(fā)現(xiàn)X1、Y1試樣被射穿，X2、Y2僅侵徹一半，未被完全射穿。X1、Y1下端面裂紋分布如圖12所示。

圖12 射孔后水泥靶材破壞形態(tài)Fig. 12 Damaged form of cement sample surface after perforation

用水潤(rùn)濕試樣表面以方便肉眼觀察宏觀裂紋。使用便攜式數(shù)碼顯微鏡（最高放大倍數(shù)為150），結(jié)合Nano Measurer軟件觀察微觀裂紋。以射孔孔道中心為圓心，圓心至靶材邊緣最近處為半徑，將試樣半徑16等分，分別以孔道中心為圓心做圓，孔道邊緣處標(biāo)為0。統(tǒng)計(jì)以孔眼中心，與各同心圓相交徑向裂紋數(shù)量及X1-B面裂縫縫寬。統(tǒng)計(jì)結(jié)果如圖13～14所示。

3.3 物理模型與數(shù)值模擬結(jié)果對(duì)比分析

物理模型與數(shù)值模擬巖石剖面裂紋形態(tài)對(duì)比及裂紋分布如圖15～16所示。

對(duì)比數(shù)值模型與物理實(shí)驗(yàn)參數(shù)設(shè)置，發(fā)現(xiàn)前者圍壓相較于后者小一個(gè)數(shù)量級(jí)時(shí)，裂縫數(shù)目分布才較相似。初步分析知，徑向裂紋始于拉伸破壞，它達(dá)到一定長(zhǎng)度后的擴(kuò)展取決于斷裂機(jī)制，后者在較小拉伸應(yīng)力作用下可以繼續(xù)延伸。本文的力學(xué)模型尚未引入斷裂機(jī)制，因此導(dǎo)致這一偏差。

分析可知物理實(shí)驗(yàn)與數(shù)值模擬裂紋分布較相似，中心一圈裂紋較少，外圍出現(xiàn)大量裂紋并呈輻射狀擴(kuò)展，后逐漸衰減。由于尺寸有限，物理實(shí)驗(yàn)靶材剖面不存在未破壞區(qū)，且靶材受反射波影響，邊界效應(yīng)明顯，剖面中存在幾條的環(huán)向裂紋，而模擬計(jì)算時(shí)忽略爆炸沖擊波且尺寸較大，沒有此類問題。統(tǒng)計(jì)兩結(jié)果裂紋數(shù)目作圖16?？梢钥闯隽藘山Y(jié)果中裂縫數(shù)目分布規(guī)律趨勢(shì)類似。

圖13 靶材剖面裂縫數(shù)目分布圖Fig. 13 The distribution of crack number on cement samples

圖14 X1-B面不同縫寬裂縫數(shù)目分布圖Fig. 14 The distribution of crack with different width on X1-B surface

圖15 數(shù)值模擬與物理實(shí)驗(yàn)靶材破壞形態(tài)對(duì)比圖Fig. 15 Comparison of damage forms between physical experiment and numerical simulation

圖16 數(shù)值模擬與物理實(shí)驗(yàn)裂縫數(shù)目分布圖Fig. 16 Crack number distribution by physical experiment and numerical simulation

4 結(jié) 論

（1）實(shí)驗(yàn)結(jié)果初步驗(yàn)證了所建立的射流侵徹巖石二維簡(jiǎn)化力學(xué)模型的有效性，并提示應(yīng)當(dāng)增加斷裂機(jī)制。

（2）根據(jù)數(shù)值模擬計(jì)算結(jié)果，以破壞類型及裂縫形態(tài)數(shù)量為依據(jù)，可將靶材切片由內(nèi)而外分為4個(gè)區(qū)域：壓剪破壞區(qū)、拉伸破壞集中區(qū)、拉伸破壞擴(kuò)展區(qū)和未破壞區(qū)。

（3）隨著侵徹孔孔徑的增大及圍壓的降低，壓剪破壞區(qū)、拉伸破壞集中區(qū)及拉伸破壞擴(kuò)展區(qū)范圍均有不同程度的增大，其中拉伸破壞擴(kuò)展區(qū)最為敏感，增長(zhǎng)最快；且侵徹孔孔徑改變時(shí)，裂紋數(shù)目無量綱半徑上的分布趨勢(shì)大致相同。