錐形PDC齒破巖機(jī)理仿真研究

豆寧輝，楊順輝

（中國石化石油工程技術(shù)研究院，北京100101）

錐形PDC齒破巖機(jī)理仿真研究

豆寧輝，楊順輝

（中國石化石油工程技術(shù)研究院，北京100101）

錐形PDC齒作為一種新型破巖方式，其相關(guān)理論在國內(nèi)研究還處于起步階段。采用DruckerPrager準(zhǔn)則，以等效塑性應(yīng)變作為破壞準(zhǔn)則，建立了基于ABAQUS的錐形PDC齒切削巖石有限元分析模型。通過數(shù)值仿真、理論分析闡述了錐形PDC齒切向上以剪切破壞為主、側(cè)向方向以擠壓拉伸斷裂破壞為主的破巖機(jī)理，以及“相鄰”效應(yīng)，并且數(shù)值仿真分析了切削深度、切削角度對錐形PDC齒破巖效率的影響規(guī)律，優(yōu)選出最優(yōu)切削角度20°。研究結(jié)果對錐形PDC齒為基礎(chǔ)的新型PDC鉆頭的設(shè)計(jì)提供了依據(jù)。

錐形PDC齒；破巖機(jī)理；數(shù)值仿真；切削角

在復(fù)雜地質(zhì)條件下的油氣鉆探工作中，經(jīng)常鉆遇碳酸鹽巖、火成巖等難鉆地層，其強(qiáng)度高，硬度大，可鉆性差，研磨性強(qiáng)?，F(xiàn)有破巖工具會存在鉆頭壽命短、鉆井周期長等問題。Simth鉆頭公司2014年推出了基于錐形PDC齒為基礎(chǔ)的StingBlade鉆頭，與常規(guī)PDC切削齒相比，其抗沖擊強(qiáng)度更大；隨著機(jī)械鉆速的增加，Stinger切削齒鉆頭的轉(zhuǎn)矩平均減少了26%；Stinger切削齒比PDC切削齒有著更加平衡的切削響應(yīng)，StingBlade鉆頭在全球300多口井中進(jìn)行了成功應(yīng)用，總作業(yè)進(jìn)尺已超過209 095 m（686 000英尺）。中國石油大學(xué)（華東）2014年進(jìn)行了錐形PDC齒鉆頭的研制及室內(nèi)試驗(yàn)評價(jià)［1-2］，但其破巖機(jī)理的理論研究尚不完善，筆者利用仿真研究方法分析了切削深度、切削角度對錐形PDC齒破巖效率的影響規(guī)律，以優(yōu)選錐形PDC齒最優(yōu)切削角度。通過對錐形PDC齒破巖機(jī)理與規(guī)律的研究，為新型PDC鉆頭的研制提供一定的理論依據(jù)。

1　錐形PDC齒破巖機(jī)理

錐形PDC齒作為一種新型石油鉆井破巖工具，既不同于牙輪鉆頭的沖擊破巖，也不同于PDC鉆頭的剪切破巖。其破巖方式是在軸向力和轉(zhuǎn)矩共同作用下連續(xù)旋轉(zhuǎn)破巖，特殊之處在于錐形PDC齒是以前傾角的角度將巖石“犁”出溝槽，在切向方向?qū)r石從母體上剝離，同時(shí)在側(cè)向方向?qū)r石產(chǎn)生擠壓破壞作用。

2　建立有限元模型

2.1優(yōu)選有限元分析軟件

ABAQUS作為最先進(jìn)的大型通用有限元分析軟件之一，能夠解決從相對簡單的線性分析到許多復(fù)雜的非線性問題，尤其是它的非線性力學(xué)（幾何、材料、接觸）分析功能具有世界領(lǐng)先水平，廣泛應(yīng)用于機(jī)械、化工、土木等各領(lǐng)域的工業(yè)設(shè)計(jì)中。

巖土介質(zhì)是具有非均質(zhì)、非線性以及幾何形狀的任意性、不連續(xù)性等因素的特殊介質(zhì)。ABAQUS能夠提供巖土工程分析中常用的本構(gòu)模型，如彈性和非線性彈性模型，MohrCoulomb模型，修正的DruckerPrager模型。

鑒于上述ABAQUS的強(qiáng)大功能，本文選擇ABAQUS作為錐形PDC齒破碎巖石有限元分析模型的平臺。

2.2模型假設(shè)

1） 將巖石看作各向同性介質(zhì)，彈塑性體，不考慮圍壓、溫度、鉆進(jìn)液等對巖石的影響。

2） 不考慮錐形PDC齒的磨損，錐形PDC齒與巖石的接觸部分為聚晶金剛石，硬度極大，視為剛性體。

3） 錐形PDC齒實(shí)際繞一定大小的旋轉(zhuǎn)半徑螺旋吃入地層，但螺旋角很小，所以將錐形PDC齒的運(yùn)動簡化為平面直線運(yùn)動，且在破巖過程中，切削速度和切削深度都固定不變。

2.3幾何模型建立與網(wǎng)格劃分

2.3.1巖石

1） 幾何模型。巖石建為長方體模型，尺寸為30 mm×20 mm×10 mm。

2） 網(wǎng)格劃分。綜合考慮計(jì)算量、計(jì)算時(shí)間以及計(jì)算精度，將巖石劃分為上下兩部分，與錐形PDC齒接觸的上部網(wǎng)格劃分較密，大小為0.5 mm，下部較疏，為1 mm，兩部分單元均采用C3D8R（8節(jié)點(diǎn)六面體線性減縮積分單元），使用結(jié)構(gòu)網(wǎng)格劃分技術(shù)（Structured），最終共劃分38 400個(gè)C3D8R單元，如圖1所示。

圖1　巖石和錐形PDC齒幾何模型

2.3.2錐形PDC齒

1） 幾何模型。錐形PDC齒設(shè)為錐頂半徑2.5 mm，錐角85°，錐齒直徑14 mm，總高20 mm。

2） 網(wǎng)格劃分。錐形PDC齒分為錐頂和其余部分，分別進(jìn)行網(wǎng)格劃分。錐頂部分，網(wǎng)格大小設(shè)為0.2mm，采用C3D4單元和自由網(wǎng)格劃分技術(shù)（Free）；其余部分網(wǎng)格大小設(shè)為0.5 mm，采用C3D8R單元和掃略網(wǎng)格劃分技術(shù)（Sweep），最終錐形PDC齒共劃分44 875個(gè)網(wǎng)格單元，如圖1所示。

2.4材料模型選取

本模型主要考慮巖石的損傷與破壞，所以將錐形PDC齒約束為剛性體，設(shè)置密度為7 800 kg／m3，彈性模量為200 GPa，泊松比為0.2。

本文采用典型砂巖的力學(xué)參數(shù)作為有限元模型中巖石的參數(shù)。巖石視為彈塑性體，彈性階段，采用線彈性本構(gòu)模型；塑性階段采用DruckerPrager準(zhǔn)則，并以等效塑性應(yīng)變作為巖石破壞準(zhǔn)則［3］，材料參數(shù)如表1所示。

表1　巖石材料屬性

2.5定義分析步

采用動態(tài)顯示分析模擬錐形PDC齒破巖過程。顯示分析不僅能夠模擬復(fù)雜的接觸，而且相比于隱式分析需要更小的時(shí)間增量步，不需要迭代和收斂準(zhǔn)則。

本文將模擬時(shí)間設(shè)置為0.05 s，并設(shè)置相關(guān)場變量和歷史變量輸出。

2.6邊界條件

對巖石底部進(jìn)行PINNED約束；對錐形PDC齒的參考點(diǎn)施加x方向的恒定速度，并約束其他5個(gè)自由度，如圖2所示。

圖2　模型邊界條件

參考實(shí)際鉆頭上切削齒的線速度，將錐形PDC齒的線速度定為0.5 m／s。

2.7定義接觸相互作用

接觸屬性中接觸對的法向行為法則采用Hard Contact，切向行為法則采用庫倫摩擦模型，摩擦因數(shù)設(shè)為0.3。

接觸采用ABAQUS／Explicit模塊中的General Contact，并在Input文件中將接觸雙方改為錐形PDC表面以及巖石內(nèi)部所有可能接觸的單元。

3　仿真結(jié)果分析

3.1巖石受力分析

本文利用錐形PDC齒破巖有限元模型，首先模擬了錐形PDC齒以15°前傾角、2mm切削深度下破壞巖石的過程。圖3～5分別為同一時(shí)刻巖石與錐形PDC齒接觸部分應(yīng)力云圖、巖石切向應(yīng)力云圖、巖石側(cè)向應(yīng)力云圖。

圖3　錐形PDC齒前巖石應(yīng)力云圖

圖4　巖石切向應(yīng)力云圖

圖5　巖石側(cè)向應(yīng)力云圖

從以上巖石應(yīng)力云圖可以看出：

1） 錐形PDC齒頂部正下方有一處應(yīng)力最大區(qū)，整體形狀為橢球形。

2） 沿錐形PDC齒運(yùn)動方向切向上，有兩處應(yīng)力最大區(qū)，分別位于錐形PDC齒正前方偏兩側(cè)，并且齒前方應(yīng)力等值線層拋物線形狀；巖石上表面破碎形狀為“∧”形，表明錐形PDC齒在切向上對巖石產(chǎn)生“犁削”作用。

3） 錐形PDC齒兩側(cè)各有一處應(yīng)力最大區(qū)域。

4） 錐形PDC齒產(chǎn)生的破碎帶側(cè)向截面呈“V”形?！癡”形底部寬度為3 mm，和錐形PDC齒正前方與巖石上表面接觸寬度相同；“V”形上部寬度為6 mm，是接觸寬度的兩倍，說明錐形PDC齒兩側(cè)對巖石產(chǎn)生了擠壓破壞作用。

根據(jù)以上數(shù)值仿真對錐形PDC齒破巖作用的分析，本文綜合彈性力學(xué)、斷裂力學(xué)，并借鑒盤形滾刀、鎬形截齒等破巖工具的破巖機(jī)理模型［4-8］，對錐形PDC齒的犁削破巖機(jī)理在切向、側(cè)向以及“相鄰效應(yīng)”作出如下闡釋。

3.2切向破巖機(jī)理分析

根據(jù)數(shù)值仿真結(jié)果，將錐形PDC齒的破巖過程分為變形及裂紋起裂、密實(shí)核形成、裂紋擴(kuò)展及塊體崩落4個(gè)階段。

1） 變形及裂紋起裂階段。

如圖6所示，當(dāng)錐形PDC齒作用于巖石時(shí)，錐形PDC齒與巖石之間會產(chǎn)生擠壓效應(yīng)，在錐頂前下方出現(xiàn)生壓縮區(qū)，同時(shí)沿切削力方向一定深度處（深色部分）產(chǎn)生最大拉應(yīng)力。在此階段，錐形PDC齒邊緣部分巖石在壓應(yīng)力作用下開始破壞，內(nèi)部拉應(yīng)力區(qū)的拉應(yīng)力達(dá)到巖石抗拉強(qiáng)度，產(chǎn)生局部破壞，形成初始裂紋。

圖6　變形及裂紋起裂階段

2） 密實(shí)核形成階段。

如圖7所示，隨著位移載荷的增加，齒前下方壓縮區(qū)巖石首先發(fā)生破碎，并被錐形PDC齒擠壓，形成密實(shí)核。

圖7　密實(shí)核形成及裂紋擴(kuò)展階段

3） 裂紋擴(kuò)展階段。

隨著位移載荷的增加，密實(shí)核作為傳力介質(zhì)向周圍巖石施加壓力，產(chǎn)生拉伸應(yīng)力。內(nèi)部裂紋沿最大拉應(yīng)力方向向自由面擴(kuò)展，如圖8所示。

圖8　裂紋擴(kuò)展階段

4） 宏觀破裂階段。

隨著載荷的進(jìn)一步增加，內(nèi)部裂紋與自由面迅速貫通，產(chǎn)生宏觀破裂，完成一次破壞過程。宏觀破裂形態(tài)如圖9所示。隨后，錐形PDC齒再次與巖石接觸，形成壓應(yīng)力區(qū)，開始下一個(gè)破壞過程［9-10］。

圖9　巖石宏觀破裂形態(tài)

3.3側(cè)向破巖機(jī)理分析

錐形PDC齒前部在剪切破巖的同時(shí)，兩側(cè)對巖體也會產(chǎn)生擠壓作用，尤其是密實(shí)核形成后，擠壓作用增加，兩側(cè)出現(xiàn)拉伸裂紋，隨著錐形PDC齒載荷的增加，拉伸裂紋主要向兩側(cè)自由面和水平方向擴(kuò)展，最終產(chǎn)生破壞［11-13］，如圖10所示。

圖10　錐形PDC齒側(cè)向破巖模型

3.4“相鄰”效應(yīng)分析

多個(gè)錐形PDC齒同時(shí)破巖過程中將會產(chǎn)生“相鄰”效應(yīng)［15-16］：一方面侵入相同深度所需的切削力比單齒要變小一些，另一方面?zhèn)认蛄鸭y的發(fā)展將沿著相鄰兩切槽間最短的直線距離發(fā)展。裂紋是否貫通與切削齒間距、切削深度或切削力、切削齒頂部本身結(jié)構(gòu)參數(shù)有關(guān)。即使裂紋不能貫通產(chǎn)生巖塊，相鄰切削齒產(chǎn)生的巖脊再次被切削時(shí)也會很容易破碎。如圖11。

圖11　“相鄰效應(yīng)”示意

3.5切削角對破巖效率影響分析

為研究切削角度在不同切削深度下對錐形PDC齒破巖效率的影響規(guī)律，本文制定如下模擬方案，如表2所示。

根據(jù)數(shù)值仿真方案，共模擬6種切削深度、6種切削角度下36種切削狀態(tài)。破巖過程如圖12所示。

本文通過計(jì)算破碎單位體積巖石所需要的能量，即破巖比能，來比較不同切削深度、切削角度下錐形PDC齒的破巖效率。破巖比能越小，破巖效率越高，反之，破巖效率越低。

表2　數(shù)值仿真方案

圖12　破巖過程

破巖能量可以通過破巖過程中的能量曲線讀取，如圖13所示，或者在XY Data選項(xiàng)中提取最后時(shí)刻破巖能量值；破巖體積通過計(jì)算失效單元的體積和得到，最終計(jì)算兩者的比值得到每種切削狀態(tài)下的破巖比能，如圖14～15所示。

從圖14～15可以得出如下結(jié)論：

1） 破巖比能隨著切削角度的增加，先減小后增加，從小到大排列依次為20、15、10、25、5、30°，從而可以看出切削角度為20°時(shí)破巖效率最高，5°和30°時(shí)破巖效率最低。

2） 破巖比能隨著切削深度的增加而減小，切削深度越大，破巖效率越高，而且當(dāng)切削深度大于1.5 mm時(shí)，破巖效率隨切削深度增加趨勢大幅減緩。

圖13　破巖過程能量曲線

圖14　切削角度與破巖比能關(guān)系曲線

圖15　切削深度與破巖比能關(guān)系曲線

4　結(jié)論

1） 建立了基于ABAQUS的錐形PDC齒切削巖石有限元分析模型，塑性階段采用DruckerPrager準(zhǔn)則，以等效塑性應(yīng)變作為巖石的破壞準(zhǔn)則。

2） 根據(jù)數(shù)值仿真、理論分析闡述了錐形PDC齒切向上以剪切破壞為主、側(cè)向上以擠壓拉伸斷裂破壞為主的破巖機(jī)理。

3） 模擬結(jié)果表明，切削角度為20°時(shí)錐形PDC齒破巖效率最高，5°和30°時(shí)破巖效率最低；破巖效率隨切削深度的增加而增加。

4） 建議繼續(xù)研究建立包含圍壓、溫度等影響因素在內(nèi)的巖石彈塑性本構(gòu)模型，并通過室內(nèi)試驗(yàn)修正模型參數(shù)。

5） 建議室內(nèi)試驗(yàn)研究不同切削角度下，錐形PDC齒的破巖效率，并與數(shù)值仿真結(jié)果進(jìn)行比較分析。

［1］宋洵成，鄒德永，管志川.PDC鉆頭等切削體積布齒優(yōu)化設(shè)計(jì)［J］.石油礦場機(jī)械，2006，35（4）：61-64.

［2］魯桂榮，肖文生，董維彬，等.PDC鉆頭CAD系統(tǒng)研究與開發(fā)［J］.石油礦場機(jī)械，2008，37（5）：22-27.

［3］楊順輝.錐形PDC齒鉆頭的研制及室內(nèi)試驗(yàn)評價(jià)［J］.石油機(jī)械，2015，43（2）：14-17.

［4］王福修.混合切削結(jié)構(gòu)PDC鉆頭簡論［J］.石油礦場機(jī)械，2005，34（3）：30-32.

［5］孫源秀，鄒德永.錐形PDC齒犁切破巖受力試驗(yàn)研究［J］.石油機(jī)械，2014，42（9）：23-26.

［6］金國棟.全斷面巖石隧道掘進(jìn)機(jī)盤刀破巖機(jī)理分析（下）［J］.電站施工機(jī)械，1981（4）：21-27.

［7］虞和濟(jì)，金國棟，蘇禾.盤形滾刀破巖機(jī)理的有限元分析［J］.東北工學(xué)院學(xué)報(bào)，1981（3）：79-86.

［8］王春華.截齒截割作用下煤體變形破壞規(guī)律研究［D］.阜新：遼寧工程技術(shù)大學(xué)，2004.

［9］王崢榮.采煤機(jī)截齒截割過程的動力學(xué)仿真［D］.太原：太原理工大學(xué)，2009.

［10］張時(shí)忠.在截刀作用下巖石破碎機(jī)理的實(shí)驗(yàn)研究［J］.巖石力學(xué)與工程學(xué)報(bào)，1990，9（3）：202-208.

［11］賴海輝.巖石破碎機(jī)理的顯微研究［J］.金屬礦山，1981（7）：2-5.

［12］東北工學(xué)院采礦系巖石破碎研究組.牙輪鉆單齒破碎巖石的初步研究［J］.東北工學(xué)院學(xué)報(bào)，1974（2）：23-45.

［13］于躍.盤刀破巖機(jī)理的細(xì)觀數(shù)值模擬研究［D］.大連：大連理工大學(xué)，2010.

［14］徐小荷，余靜.巖石破碎學(xué)［M］.北京：煤炭工業(yè)出版社，1984：191-192.

［15］林家驥.研究盤形滾刀破巖機(jī)理與提高刀具使用壽命［J］.福州大學(xué)學(xué)報(bào)（自然科學(xué)版），1990，18（1）：108-113.

［16］劉德順.滾刀破巖拉槽間距研究［J］.焦作礦業(yè)學(xué)院學(xué)報(bào)，1993（3）：86-91.

Simulation Research of Conical PDC Cutter Rock-Breaking Mechanism

DOU Ninghui，YANG Shunhui
（Sinopec Research Institute of Petroleum Engineering，Beijing 100101，China）

Conical PDC cutter as a new break the rock，the theory is still in initial phase in China. By choosing DruckerPrager criterion as the rock constitutive model，equivalent plastic strain as the failure criterion，the FEA model of rockcutting by conical PDC cutter was established，which was based on ABAQUS.Through numerical simulation and theoretical analysis，the rockbreaking mechanism of conical PDC cutter was elaborated：shear failure in tangential and tensile fracture failure in lateral direction，“adjacent”effect was also elaborated.The influence laws between rockbreaking efficiency and cutting depth as well as cutting angle were simulated，and the optimal cutting angle 20°was selected.The rock breaking model of the conical PDC cutter is plowbreaking of the PDC bit.

conical PDC cutter；rockbreaking mechanism；numerical simulation；cutting angle

TE921.102

A

10.3969／j.issn.1001.3482.2015.11.003

1001-3482（2015）11-0012-06

2015-05-03

國家科技重大專項(xiàng)“中東富油氣區(qū)復(fù)雜地層井筒關(guān)鍵技術(shù)”（2011ZX05031-004）

豆寧輝（1984），男，甘肅慶陽人，工程師，2010年畢業(yè)于中國石油大學(xué)（北京）油氣井工程專業(yè)，主要從事鉆頭、深水鉆井、井壁穩(wěn)定等方面研究，Email：dounh.sripe＠sinopec.com。