基于擴(kuò)展PFC2D-GBM模型的單齒切削花崗巖破碎機(jī)制

黃繼慶　胡?！》汲伞ゼ∽ＰA

摘要：常規(guī)鉆井工具在鉆進(jìn)深部油氣資源地層巖石的過(guò)程中，通常會(huì)遇到鉆速低、成本高等問(wèn)題。為提高鉆頭在深部硬地層的鉆進(jìn)效率、優(yōu)化破巖參數(shù)，利用擴(kuò)展PFC2D-GBM模型建立基于真實(shí)晶粒結(jié)構(gòu)的花崗巖模型，對(duì)PDC單齒切削作用下的破巖機(jī)制進(jìn)行深入分析，分別從宏細(xì)觀尺度討論圍壓、前傾角以及切削速度對(duì)破巖效果的影響。結(jié)果表明：在切削過(guò)程中，巖石中總會(huì)萌生晶內(nèi)拉伸裂紋、晶內(nèi)剪切裂紋、晶間拉伸裂紋以及晶間剪切裂紋4種微裂紋；無(wú)論在何種工況下，晶間拉伸裂紋的數(shù)量總是最多且分布范圍最廣；在無(wú)圍壓時(shí)，礦物晶粒在PDC齒作用下沿晶粒邊界脫離巖石形成巖屑，產(chǎn)生的裂紋基本是晶間拉伸裂紋，只存在極少的晶內(nèi)裂紋；隨著圍壓的增大，晶內(nèi)剪切裂紋的萌生會(huì)得到一定的促進(jìn)作用，并且圍壓越大，這種促進(jìn)作用更強(qiáng)。

關(guān)鍵詞：硬巖； 花崗巖； PFC2D； 細(xì)觀破碎； 數(shù)值模擬

中圖分類(lèi)號(hào)：TE 921 文獻(xiàn)標(biāo)志碼：A

引用格式：黃繼慶，胡海，樊思成，等.基于擴(kuò)展PFC2D-GBM模型的單齒切削花崗巖破碎機(jī)制［J］.中國(guó)石油大學(xué)學(xué)報(bào)（自然科學(xué)版），2023，47（2）：81-89.

HUANG Jiqing， HU Hai， FAN Sicheng， et al. Rock cutting mechanisms of granite using an extended PFC2D-GBM model［J］.Journal of China University of Petroleum（Edition of Natural Science），2023，47（2）：81-89.

Rock cutting mechanisms of granite using an extended PFC2D-GBM model

HUANG Jiqing1， HU Hai2， FAN Sicheng1， LIU Weiji2， ZHU Xiaohua2

（1.CNPC Bohai Petroleum Equipment Manufacturing Company Limited， Tianjin 300457，China; 2.School of Mechatronic Engineering， Southwest Petroleum University， Chengdu 610500， China）

Abstract： Low penetration rate and high cost are the main problems encountered in drilling of deep oil and gas formation using conventional drilling tools. In order to improve the drilling efficiency of drill bits in deep and hard rocks and to optimize the rock breaking parameters， a granite rock model based on real grain structures was established using an extended PFC2D-GBM model， and the rock breaking mechanisms with a PDC single-tooth were studied. The effects of confining pressure， rake angle and cutting speed on granite fragmentation were discussed at macro and micro scales. The results indicate that， in the rock cutting process， four types of micro cracks can be generated in the rock matrix，i.e. intragranular tensile crack， intragranular shear crack， intergranular tensile crack and intergranular shear crack. The number of the intergranular tensile cracks is always the largest and its distribution range is also the widest under any working conditions. When there is no confining pressure， mineral grains can be broken off from the rock matrix along the grain boundaries under the action of the PDC teeth， and the cracks generated are basically the intergranular tensile cracks， only with a few intragranular cracks. However， with the increase of the confining pressure， the initiation of intragranular shear cracks can be promoted， and the higher the confining pressure， the more of the intragranular shear cracks.

Keywords：hard rock; granite; PFC2D; meso-fragmentation; numerical simulation

作為世界能源結(jié)構(gòu)的主體，石油、天然氣的消耗量隨著世界經(jīng)濟(jì)不斷發(fā)展而日益增加。由于不斷的開(kāi)采，淺層容易被開(kāi)發(fā)的常規(guī)油氣資源越來(lái)越少，高效開(kāi)發(fā)深層、超深層油氣資源成為中國(guó)實(shí)現(xiàn)能源接替戰(zhàn)略的重大需求，是當(dāng)前和未來(lái)油氣勘探開(kāi)發(fā)的重點(diǎn)和熱點(diǎn)［1-2］。油氣資源的勘探開(kāi)發(fā)從淺層向深層、超深層發(fā)展的過(guò)程中，鉆遇硬度高、研磨性高、可鉆性級(jí)值高的巖石的概率越來(lái)越大［3］。牙輪鉆頭、金剛石鉆頭等常規(guī)破巖工具在深井、超深井的鉆進(jìn)過(guò)程中，通常會(huì)遇到鉆頭磨損嚴(yán)重、鉆井成本高等問(wèn)題［4］，這嚴(yán)重影響了深部硬地層的機(jī)械鉆速，從而導(dǎo)致勘探開(kāi)發(fā)成本大量增加。在常規(guī)破巖鉆頭中，PDC鉆頭相對(duì)于牙輪鉆頭等其他鉆頭具有較高硬度、耐磨性好、結(jié)構(gòu)簡(jiǎn)單和設(shè)計(jì)靈活的特點(diǎn)［5］，具有更大的改進(jìn)空間。為提高PDC鉆頭在硬巖地層中的鉆進(jìn)效率，許多學(xué)者和專(zhuān)業(yè)人員進(jìn)行了大量的嘗試和試驗(yàn)［6-7］，而眾多數(shù)值模擬試驗(yàn)在鉆頭的設(shè)計(jì)和優(yōu)化過(guò)程中功不可沒(méi)［8-11］。然而，許多模擬試驗(yàn)大多將巖石考慮成均質(zhì)結(jié)構(gòu)，忽略巖石內(nèi)部細(xì)觀結(jié)構(gòu)的影響，在一定程度上制約了對(duì)鉆頭結(jié)構(gòu)的進(jìn)一步有效優(yōu)化。為進(jìn)一步研究巖石在PDC齒切削作用下的破碎機(jī)制，筆者以典型的深部硬巖—花崗巖［12-13］為研究對(duì)象，引入巖石內(nèi)部細(xì)觀晶粒真實(shí)的結(jié)構(gòu)及其分布狀態(tài)，建立花崗巖擴(kuò)展PFC2D-GBM模型，討論在PDC單齒線(xiàn)性切削作用下，圍壓、前傾角以及切削速度對(duì)花崗巖破碎情況的影響，以期為PDC鉆頭的進(jìn)一步優(yōu)化以及深部油氣資源的高效開(kāi)采過(guò)程提供新思路、新方法。

1 花崗巖模型建立與標(biāo)定

1.1 擴(kuò)展PFC2D-GBM模型

破碎機(jī)制的研究必然與巖石的細(xì)觀力學(xué)行為等方面有關(guān)。然而，通過(guò)室內(nèi)試驗(yàn)只能觀測(cè)到巖石的宏觀力學(xué)行為，對(duì)于巖石內(nèi)部微裂紋的擴(kuò)展等現(xiàn)象卻難以觀測(cè)，并且室內(nèi)試驗(yàn)無(wú)法實(shí)現(xiàn)完全相同條件下的重復(fù)性操作。近年來(lái)，通過(guò)數(shù)值方法再現(xiàn)巖石破碎的完整過(guò)程成為解決上述問(wèn)題的新方法［14-16］。PFC（particle flow code）作為基于非連續(xù)介質(zhì)力學(xué)的代表性軟件［17］，不僅可以模擬巖石材料的宏觀特性和細(xì)觀力學(xué)特性，而且可以很方便地揭示細(xì)觀層面上巖石的動(dòng)態(tài)力學(xué)特性［18］。在PFC中，巖石由微小的、離散的顆粒組成，顆粒與顆粒之間通過(guò)特定的力學(xué)接觸而相互聯(lián)系在一起。當(dāng)有外力作用時(shí)，各個(gè)顆粒之間的接觸被壓縮、分離或滑動(dòng)，顆粒因此會(huì)發(fā)生平移、旋轉(zhuǎn)和變形，眾多微小顆粒的細(xì)觀力學(xué)行為匯集成巖石的宏觀力學(xué)行為。為使巖石模型貼近巖石的真實(shí)結(jié)構(gòu)，Potyondy［19］在PFC2D中建立了一種基于晶粒結(jié)構(gòu)的復(fù)雜巖石模型（grain based model， GBM）。GBM極大的促進(jìn)了對(duì)巖石細(xì)觀力學(xué)行為的研究［20-21］。

基于上述模型，將傳統(tǒng)GBM中表征晶粒形狀的多邊形替換為真實(shí)的花崗巖晶粒輪廓圖，使模型的微觀結(jié)構(gòu)更貼近真實(shí)的巖石。采用的巖樣來(lái)自湖北隨州的花崗巖，主要礦物組分為斜長(zhǎng)石（21.9%）、鈉長(zhǎng)石（47.5%）、石英（19.3%）以及黑云母（9.3%），其晶粒輪廓明顯、礦物成分易區(qū)分。按照劉偉吉等［22］提出的方法建立基于真實(shí)晶粒輪廓的花崗巖擴(kuò)展PFC2D-GBM模型。模型建立過(guò)程如圖1所示。

1.2 細(xì)觀參數(shù)標(biāo)定

PFC中的數(shù)值模型能用于探究巖石破碎機(jī)制的前提，是必須能夠再現(xiàn)特定的宏觀力學(xué)行為。根據(jù)圖1所示的過(guò)程建立非均質(zhì)花崗巖模型，并按礦物種類(lèi)在合理范圍內(nèi)賦予顆粒力學(xué)參數(shù)。通過(guò)單軸壓縮試驗(yàn)和巴西劈裂試驗(yàn)獲得巖石模型的抗壓強(qiáng)度、彈性模量以及抗拉強(qiáng)度。用于單軸壓縮試驗(yàn)的模型是一個(gè)矩形，其高度為50 mm、寬度為25 mm，包含35 078個(gè)顆粒；用于巴西劈裂試驗(yàn)的模型則是一直徑為25 mm的圓，包含13 644個(gè)顆粒。標(biāo)定試驗(yàn)的結(jié)果如圖2、3所示。

兩個(gè)標(biāo)定試驗(yàn)的抗壓強(qiáng)度為84.8 MPa、彈性模量為10 451.4 MPa、抗拉強(qiáng)度為5.31 MPa。通過(guò)室內(nèi)試驗(yàn)得到抗壓強(qiáng)度為86.37 MPa、彈性模量為10 707.6 MPa、抗拉強(qiáng)度為4.85 MPa?？箟簭?qiáng)度、彈性模量、抗拉強(qiáng)度標(biāo)定結(jié)果與試驗(yàn)結(jié)果間誤差分別為1.82%、2.39%、9.4%，誤差均小于10%，證明相關(guān)參數(shù)能用于后續(xù)試驗(yàn)。標(biāo)定模型中顆粒的最小半徑為7.5×10-5 m、最大半徑為1.25×10-4 m，顆粒密度為2 700 kg/m3，摩擦系數(shù)為0.5，平行鍵半徑倍數(shù)為1，平行鍵摩擦角為30°。模型中的其他微觀參數(shù)［22］如表1所示。

2 切削破巖數(shù)值模擬

2.1 切削破巖模型

基于花崗巖模型，采用表1中的各項(xiàng)細(xì)觀參數(shù)并導(dǎo)入PDC齒，建立如圖4所示的切削破巖模型。該模型包含7 716個(gè)顆粒，長(zhǎng)度為20 mm，高度為12 mm。PDC齒的直徑為12 mm。討論花崗巖地層處于平衡壓力鉆井時(shí)的破碎情況，該狀態(tài)下井筒內(nèi)的液柱壓力與地層壓力相等。將花崗巖模型的底邊固定，通過(guò)兩側(cè)的剛性壁及上方的柔性邊界顆粒鏈向其施加圍壓pc。兩側(cè)的剛性壁上施加的圍壓代表地層壓力，上方柔性邊界顆粒鏈?zhǔn)┘拥膲毫Υ磴@井過(guò)程中的液柱壓力。PDC齒與巖石之間的接觸不會(huì)受到柔性邊界顆粒鏈的影響，因?yàn)镻DC齒能夠穿過(guò)柔性邊界顆粒鏈與花崗巖直接接觸。

2.2 結(jié)果分析

研究圍壓、切削速度以及前傾角對(duì)花崗巖在PDC單齒線(xiàn)性切削作用下破碎情況的影響。圍壓pc為0～6 MPa，增量為2 MPa；前傾角θ為5°～20°，增量為5°；切削速度υ為0.5～2.0 m/s，增量為0.5 m/s。PDC齒切削花崗巖的整個(gè)過(guò)程中切削深度為1 mm，切削行程為5 mm。

2.2.1 巖屑與裂紋

圖5為υ=1.0 m/s時(shí)，不同pc和θ作用下花崗巖的破碎情況，圖6為θ=5°時(shí)，不同pc和υ作用下花崗巖的破碎情況。圖5、6展示了巖屑的生成以及不同種類(lèi)裂紋的萌生和分布情況。為更好地展示巖石的破碎情況，其中未顯示柔性邊界顆粒鏈。

由圖5可以看出：pc=0 MPa時(shí)產(chǎn)生的巖屑尺寸都比較大，并且基本沿PDC齒的切削面分布；當(dāng)pc≠0 MPa時(shí)，巖屑的尺寸都比較小。在整個(gè)切削過(guò)程中，巖石中萌生晶內(nèi)拉伸裂紋、晶內(nèi)剪切裂紋、晶間拉伸裂紋以及晶間剪切裂紋4種微裂紋。

由圖5、6可知，無(wú)論在何種工況下，總是會(huì)萌生大量的晶間拉伸裂紋，其數(shù)量最多且分布最廣，晶間拉伸裂紋不僅在直接損傷區(qū)大量萌生，也廣泛分布在遠(yuǎn)離直接損傷區(qū)的更深處的各晶粒邊界上。pc=0 MPa時(shí)，PDC齒主要造成晶粒與晶粒之間的破壞，各晶粒在PDC齒作用下沿晶粒邊界脫離巖石形成巖屑，巖屑的尺寸普遍較大。此時(shí)在直接損傷區(qū)產(chǎn)生的裂紋基本是晶間拉伸裂紋，只存在非常少的晶內(nèi)裂紋；而隨著圍壓增大，直接損傷區(qū)內(nèi)晶內(nèi)裂紋的萌生得到促進(jìn)，并且圍壓越大，這種促進(jìn)作用越強(qiáng)。這是由于在PDC齒將礦物晶粒從巖石上剝離后，柔性邊界顆粒鏈向下施加的壓力使PDC齒再次與巖屑相互作用，導(dǎo)致巖屑發(fā)生二次破碎，從而產(chǎn)生大量晶內(nèi)裂紋。并且當(dāng)pc=0 MPa、θ=5°時(shí)，切削過(guò)程產(chǎn)生的晶內(nèi)裂紋數(shù)量最少。pc≠0 MPa時(shí)，在同一圍壓條件下，不同前傾角對(duì)花崗巖造成的破碎相差并不是很明顯。

與圖5展示的情況相似，圖6中pc=0 MPa時(shí)的巖屑尺寸也普遍較大。而在同一圍壓條件下，不同切削速度對(duì)花崗巖造成的破壞相差并不明顯。但隨著圍壓的增大，同一切削速度下花崗巖的破碎難度增加，因?yàn)閲鷫涸龃髸?huì)使巖石內(nèi)部的微孔隙逐漸閉合，各礦物晶粒之間變得更加致密，PDC齒對(duì)花崗巖造成破壞的難度逐漸加大。在pc沒(méi)有達(dá)到極限值以前這種趨勢(shì)不會(huì)發(fā)生變化。

2.2.2 切削力與裂紋變化

切削花崗巖的過(guò)程中PDC齒會(huì)受到巖石對(duì)其產(chǎn)生的作用力，即切削力。通過(guò)對(duì)PDC齒和巖石模型的動(dòng)態(tài)監(jiān)測(cè)，可以得到切削力以及各類(lèi)型裂紋的變化情況。圖7為υ=1.0 m/s時(shí)，不同pc和θ作用下切削力與裂紋數(shù)量的變化曲線(xiàn)。圖8為切削模型局部展示。圖9為θ=5°時(shí)，不同pc和υ作用下切削力與裂紋數(shù)量的變化曲線(xiàn)。

在pc=0 MPa的條件下，切削力變化曲線(xiàn)波峰的數(shù)量最少，并且波峰的峰值小，即無(wú)圍壓狀態(tài)下切削力的頻率和幅度都比較小。這是因?yàn)榫Я＿吔绲牧W(xué)參數(shù)比晶粒內(nèi)部各粒子之間的力學(xué)參數(shù)弱，PDC齒對(duì)花崗巖造成的破壞主要發(fā)生在各晶粒邊界處，產(chǎn)生的巖屑通常是整塊的礦物晶粒，這與巖屑尺寸在無(wú)圍壓狀態(tài)下普遍偏大的現(xiàn)象相吻合。從圖9可知，pc=0 MPa時(shí)，切削力變化曲線(xiàn)的波峰數(shù)量會(huì)隨著切削速度的增大而增多，但每個(gè)波峰的峰值并不大，這意味著PDC齒的切削速度越大，無(wú)圍壓狀態(tài)下切削力的變化頻率越高。然而，圖7顯示前傾角的增大又會(huì)減小切削力變化曲線(xiàn)的波峰數(shù)量，即前傾角的增大可以降低切削力的變化頻率。另外，由圖7和9可以看出，在pc≠0 MPa時(shí)，在切削行程超過(guò)3 mm后，切削力變化曲線(xiàn)在大多數(shù)情況下會(huì)出現(xiàn)明顯的下降。由圖8可知，當(dāng)切削行程為3 mm時(shí)，PDC切削齒直接作用在3號(hào)晶粒上，而由圖7與9中pc=0 MPa時(shí)的破碎情況可知，3號(hào)晶粒幾乎是整塊從花崗巖上脫離，因此3號(hào)與4號(hào)晶粒之間的邊界比較薄弱， PDC齒極易對(duì)其造成破壞，因此造成了上述切削力變化曲線(xiàn)的下降。

從圖7中可以看出，隨著θ的不斷增大，相同圍壓條件下切削力變化曲線(xiàn)中大峰值波峰的數(shù)量減小，相應(yīng)地，峰值平均值減小。在θ為5°和20°時(shí)，相同圍壓條件下的切削力平均值最小，且波峰的分布更加均勻，這表明PDC齒受到的力最小且受力的頻率更穩(wěn)定，PDC齒受到的損傷最小。由圖9可知，當(dāng)υ=0.5 m/s時(shí)，有圍壓狀態(tài)下切削力曲線(xiàn)在切削過(guò)程的初期總是會(huì)出現(xiàn)峰值很高的大波峰；在υ=1.5 m/s時(shí)，切削力曲線(xiàn)在有圍壓狀態(tài)下的波谷在整個(gè)切削過(guò)程中始終維持在一個(gè)較高水平，PDC齒處在持續(xù)受力的狀態(tài)；而與υ=1.0 m/s時(shí)的情況相比，υ=2.0 m/s時(shí)切削力曲線(xiàn)的波峰的峰值更大。綜合切削力曲線(xiàn)在4個(gè)切削速度與不同圍壓條件下的變化情況，PDC齒在υ=1.0 m/s時(shí)受力最小。

在所有工況下，裂紋總數(shù)在切削過(guò)程初期總是會(huì)出現(xiàn)驟增，此時(shí)PDC齒與花崗巖之間發(fā)生相互作用，切削力出現(xiàn)第一個(gè)波峰。此后，裂紋總數(shù)的變化與切削力曲線(xiàn)波峰的變化情況直接相關(guān)：當(dāng)切削力曲線(xiàn)出現(xiàn)波峰時(shí)，裂紋總數(shù)增加。在大多數(shù)情況下，波峰的峰值越大且大峰值波峰越密集時(shí)裂紋總數(shù)的增加越明顯。而當(dāng)切削力曲線(xiàn)不出現(xiàn)波峰或波峰的峰值處于較低水平時(shí)，各種類(lèi)型的裂紋數(shù)量幾乎不發(fā)生變化。晶內(nèi)剪切裂紋與晶間拉伸裂紋的數(shù)量與裂紋總數(shù)變化趨勢(shì)一致，

不同工況下晶間拉伸裂紋的數(shù)量與裂紋總數(shù)變化趨勢(shì)一致。

4種微裂紋中，晶間拉伸裂紋的數(shù)量從始至終都處在最高水平，晶內(nèi)剪切裂紋的數(shù)量次之，晶內(nèi)拉伸裂紋數(shù)量始終處于較低水平，而晶間剪切裂紋數(shù)量基本為0且?guī)缀鯖](méi)有變化。當(dāng)其他條件相同時(shí)，晶內(nèi)剪切裂紋的數(shù)量總是會(huì)隨圍壓的增大而增大。晶內(nèi)拉伸裂紋與晶間拉伸裂紋隨圍壓的變化規(guī)律不太明顯。

綜上可知，圍壓的增大會(huì)促進(jìn)晶內(nèi)剪切裂紋的萌生，并且當(dāng)υ=1.0 m/s時(shí)，在θ=5°與θ=20°的情況下，切削力的最大峰值相對(duì)其他情況更小，且波峰分布更加平均，此時(shí)PDC齒受力最小且受力頻率更穩(wěn)定。

2.2.3 破碎比

為進(jìn)一步描述花崗巖的破碎情況，引入破碎比K來(lái)表示花崗巖在切削過(guò)程受到的破壞，K的計(jì)算公式為

K=Sd/St.

式中，Sd為發(fā)生破壞的顆?？偯娣e，m2；St為整個(gè)模型包含的所有顆粒的面積總和，m2。

圖10為υ=1.0 m/s時(shí)不同pc和θ作用下破碎比的變化情況。由圖10（a）可以看出，不同切削角度下的破碎比隨著圍壓的增大而減小。這是由于隨著圍壓的增大，巖石內(nèi)部的微孔隙逐漸閉合，各礦物晶粒之間變得更加致密，PDC齒對(duì)花崗巖造成破壞的難度逐漸加大。而在pc沒(méi)有達(dá)到極限值以前，這種趨勢(shì)不會(huì)發(fā)生變化。由圖10（b）可以看出，相同圍壓條件下破碎比隨θ的增大呈波浪式變化，即K隨θ的增大先減小后增大，然后再減小。在不同工況下，θ=5°時(shí)的破巖效率最高。

圖11為θ=5°時(shí)不同pc和υ作用下K的變化情況。由圖11（a）可以看出，同一圍壓條件下，不同切削速度的破碎比隨著圍壓的增大而有不斷減小的趨勢(shì)，這與圖10中情況相似，都是由于巖石內(nèi)部的微孔隙隨著圍壓的增大而逐漸閉合造成的。相對(duì)于其他3種切削速度，υ=1.0 m/s時(shí)的K大多數(shù)情況下都是最大的。由圖11（b）可以看出，同一圍壓條件下切削速度的改變對(duì)破碎比的改變并不明顯。在沒(méi)有圍壓時(shí)，切削速度的增加會(huì)導(dǎo)致破碎比逐漸減小。

3 結(jié) 論

（1）在整個(gè)切削過(guò)程中，巖石中會(huì)萌生晶內(nèi)拉伸裂紋、晶內(nèi)剪切裂紋、晶間拉伸裂紋以及晶間剪切裂紋4種微裂紋。其中晶間拉伸裂紋數(shù)量最多且分布最廣，不僅在直接損傷區(qū)大量萌生，也廣泛分布在遠(yuǎn)離直接損傷區(qū)的更深處的各晶粒邊界上。

（2）不存在圍壓時(shí)，PDC齒主要造成晶粒之間的破壞，礦物晶粒在PDC齒作用下沿晶粒邊界脫離巖石形成巖屑，此時(shí)產(chǎn)生的裂紋基本是晶間拉伸裂紋，只存在很少的晶內(nèi)裂紋。而隨著圍壓的增大，晶內(nèi)剪切裂紋的萌生會(huì)得到一定的促進(jìn)作用，并且圍壓越大，這種促進(jìn)作用越強(qiáng)。

（3）當(dāng)υ=1.0 m/s、θ=5°時(shí)，PDC齒受力最小且受力的頻率更穩(wěn)定，并能保證較高的破巖效率。

參考文獻(xiàn)：

［1］陳浩，李邦潤(rùn)，蔡燦，等.深井高溫巖石破巖機(jī)理及生熱分析［J］.石油機(jī)械，2021，49（1）：1-10.

CHEN Hao， LI Bangrun， CAI Can， et al. Rock breaking mechanisms and heat generation analysis on high temperature rocks in deep wells［J］. China Petroleum Machinery， 2021，49（1）：1-10.

［2］胡文瑞，鮑敬偉.石油行業(yè)發(fā)展趨勢(shì)及中國(guó)對(duì)策研究［J］.中國(guó)石油大學(xué)學(xué)報(bào)（自然科學(xué)版），2018，42（4）：1-10.

HU Wenrui， BAO Jingwei. Development trends of oil industry and Chinas countermeasures［J］.Journal of China University of Petroleum（Edition of Natural Science），2018，42（4）：1-10.

［3］黃家根，汪海閣，紀(jì)國(guó)棟，等.超聲波高頻旋沖鉆井技術(shù)破巖機(jī)理研究［J］.石油鉆探技術(shù)，2018，46（4）：23-29.

HUANG Jiagen， WANG Haige， JI Guodong， et al. The rock breaking mechanism of ultrasonic high frequency rotary-percussive drilling technology［J］. Petroleum Drilling Techniques， 2018，46（4）：23-29.

［4］王濱，鄒德永，李軍，等.深部及復(fù)雜地層中PDC鉆頭綜合改進(jìn)方法［J］.石油鉆采工藝，2018，40（1）：44-51.

WANG Bin， ZOU Deyong， LI Jun， et al. A comprehensive method to improve the performance of PDC bits in deep and complex formations［J］. Oil Drilling & Production Technology， 2018，40（1）：44-51.

［5］唐勝利，曹小軍，高欣，等.PDC切削齒后傾角與破巖關(guān)系的數(shù)值模擬［J］.煤炭技術(shù)，2019，38（4）：154-157.

TANG Shengli， CAO Xiaojun， GAO Xin， et al. Numerical simulation of relationship between back rake of PDC cutter and rock breaking［J］. Coal Technology， 2019，38（4）：154-157.

［6］鄒德永，蔡環(huán).布齒參數(shù)對(duì)PDC鉆頭破巖效率影響的試驗(yàn)［J］.中國(guó)石油大學(xué)學(xué)報(bào)（自然科學(xué)版），2009，33（5）：76-79.

ZOU Deyong， CAI Huan. Experiment on effect of cutter parameters of PDC bit on rate of penetration［J］.Journal of China University of Petroleum（Edition of Natural Science）， 2009，33（5）：76-79.

［7］王家駿，鄒德永，楊光，等.PDC切削齒與巖石相互作用模型［J］.中國(guó)石油大學(xué)學(xué)報(bào)（自然科學(xué)版），2014，38（4）：104-109.

WANG Jiajun， ZOU Deyong， YANG Guang， et al. Interaction model of PDC cutter and rock［J］.Journal of China University of Petroleum（Edition of Natural Science）， 2014，38（4）：104-109.

［8］祝效華，但昭旺.PDC切削齒破碎干熱巖數(shù)值模擬［J］.天然氣工業(yè)，2019，39（4）：125-134.

ZHU Xiaohua， DAN Zhaowang. Numerical simulation of rock breaking by PDC cutters in hot dry rocks［J］. Natural Gas Industry， 2019，39（4）：125-134.

［9］祝效華，李海.PDC切削齒破巖效率數(shù)值模擬研究［J］.應(yīng)用基礎(chǔ)與工程科學(xué)學(xué)報(bào)，2015，23（1）：182-191.

ZHU Xiaohua， LI Hai. Numerical simulation on mechanical special energy of PDC cutter rock-cutting［J］. Journal of Basic Science and Engineering， 2015，23（1）：182-191.

［10］馬清明，王瑞和.PDC切削齒破巖受力的試驗(yàn)研究［J］.中國(guó)石油大學(xué)學(xué)報(bào)（自然科學(xué)版），2006，30（2）：45-47.

MA Qingming， WANG Ruihe. Experimental study on force of PDC cutter breaking rock［J］. Journal of China University of Petroleum（Edition of Natural Science）， 2006，30（2）：45-47.

［11］景英華，袁鑫偉，姜磊，等.高速旋轉(zhuǎn)沖擊鉆井破巖數(shù)值模擬及現(xiàn)場(chǎng)實(shí)驗(yàn)［J］.中國(guó)石油大學(xué)學(xué)報(bào)（自然科學(xué)版），2019，43（1）：75-80.

JING Yinghua， YUAN Xinwei， JIANG Lei， et al. Numerical simulation and field experimental study on rock breaking in high speed rotating percussion drilling［J］. Journal of China University of Petroleum（Edition of Natural Science），2019，43（1）：75-80.

［12］宮鳳強(qiáng)，伍武星，李天斌，等.深部硬巖矩形隧洞圍巖板裂破壞的試驗(yàn)?zāi)M研究［J］.巖土力學(xué)，2019，40（6）：2085-2098.

GONG Fengqiang， WU Wuxing， LI Tianbin， et al. Simulation experimental study of spalling failure of surrounding rock of rectangular tunnel of deep hard rock［J］. Rock and Soil Mechanics，2019，40（6）：2085-2098.

［13］李地元，胡楚維，朱泉企.預(yù)制裂隙花崗巖動(dòng)靜組合加載力學(xué)特性和破壞規(guī)律試驗(yàn)研究［J］.巖石力學(xué)與工程學(xué)報(bào)，2020，39（6）：1081-1093.

LI Diyuan， HU Chuwei， ZHU Quanqi. Experimental study on mechanical properties and failure laws of granite with an artificial flaw under coupled static and dynamic loads［J］.Chinese Journal of Rock Mechanics and Engineering， 2020，39（6）：1081-1093.

［14］祝效華，劉偉吉.楔形鉆齒作用下巖石破碎及微觀劣化過(guò)程［J］.中國(guó)石油大學(xué)學(xué)報(bào)（自然科學(xué)版），2017，41（1）：90-95.

ZHU Xiaohua， LIU Weiji. Numerical simulation of rock breaking and microscopic deterioration under wedge-shaped drill bit teeth［J］. Journal of China University of Petroleum （Edition of Natural Science）， 2017，41（1）：90-95.

［15］FAIRHURST C. Rock engineering： where is the laboratory？［J］. Rock Mechanics and Rock Engineering， 2019，52（12）：4865-4888.

［16］趙吉坤，張子明，仝興華，等.巖石宏細(xì)觀彈塑性損傷破壞對(duì)比研究［J］.中國(guó)石油大學(xué)學(xué)報(bào)（自然科學(xué)版），2007，31（4）：78-84.

ZHAO Jikun， ZHANG Ziming， TONG Xinghua， et al. Comparison of macro-scale and meso-scale elasto-plastic damage failure in rocks［J］. Journal of China University of Petroleum （Edition of Natural Science）， 2007，31（4）：78-84.

［17］CUNDALL P A， STRACK O D L. A discrete numerical mode for granular assemblies［J］. Géotechnique， 1979，29（1）：

47-65.

［18］HUANG Haiying， DETOURNAY E. Discrete element modeling of tool-rock interaction II： rockindentation［J］. International Journal for Numerical and Analytical Methods in Geomechanics， 2013，37（13）：1930-1947.

［19］POTYONDY D O. A grain-based model for rock： approaching the true microstructure： Proceedings of Bergmekanikki Norden 2010—Rock Mechanics in the Nordic Countries 2010［C］. Kongsberg： Norwegian Group for Rock Mechanics， 2010.

［20］ZHOU Jian， LAN Hengxing， ZHANG Luqing， et al. Novel grain-based model for simulation of brittle failure of alxa porphyritic granite［J］. Engineering Geology， 2019，251：100-114.

［21］SAADAT M， TAHERI A. A numerical approach to investigate the effects of rock texture on the damage and crack propagation of a pre-cracked granite［J］. Computers & Geotechnics， 2019，111：89-111.

［22］劉偉吉，胡海，祝效華，等.圍壓條件下電脈沖破碎干熱巖機(jī)理［J］.天然氣工業(yè)，2022，42（4）：117-129.

LIU Weiji， HU Hai， ZHU Xiaohua， et al. Breaking mechanisms of hot dry rock by electric pulse under confining pressure［J］. Natural Gas Industry， 2022，42（4）：117-129.

（編輯 李志芬）

收稿日期：2022-09-20

基金項(xiàng)目：國(guó)家能源局能源領(lǐng)域“科技助力經(jīng)濟(jì)2020”重點(diǎn)專(zhuān)項(xiàng)（9）;國(guó)家自然科學(xué)基金項(xiàng)目（52004229）

第一作者：黃繼慶（1969-），男，教授級(jí)高級(jí)工程師，碩士，研究方向?yàn)槭蜋C(jī)械技術(shù)開(kāi)發(fā)和管理。E-mail：huangjiqing@cnpc.com.cn。

通信作者：劉偉吉（1989-），男，副教授，博士，研究方向?yàn)楦咝茙r理論與方法、鉆井提速提效技術(shù)。E-mial;lwj2017_swpu@163.com。