磨損狀況下聚晶金剛石鉆頭切削齒的切削參數(shù)求解

馬亞超 黃志強(qiáng) 李 琴 彭世金 謝 豆

1.西南石油大學(xué)機(jī)電工程學(xué)院，成都，6105002.西南石油大學(xué)實(shí)驗(yàn)室與設(shè)備管理處，成都，610500

磨損狀況下聚晶金剛石鉆頭切削齒的切削參數(shù)求解

馬亞超1黃志強(qiáng)1李 琴1彭世金2謝 豆1

1.西南石油大學(xué)機(jī)電工程學(xué)院，成都，6105002.西南石油大學(xué)實(shí)驗(yàn)室與設(shè)備管理處，成都，610500

針對(duì)磨損狀況下PDC鉆頭切削齒的切削參數(shù)預(yù)測(cè)方法尚未建立，無(wú)法很好地開展等磨損布齒設(shè)計(jì)的問(wèn)題，建立了磨損齒的幾何模型，并提出了求解切削參數(shù)的零點(diǎn)遍歷法?；贛ATLAB平臺(tái)對(duì)零點(diǎn)遍歷法進(jìn)行編譯，并對(duì)切削齒的切削弧長(zhǎng)、切削面積和切削體積進(jìn)行求解。結(jié)果表明，提出的零點(diǎn)遍歷法能夠用于磨損齒和未磨損齒的切削參數(shù)求解(誤差在2%以下)；結(jié)合切削齒的受力與磨損模型，能夠預(yù)測(cè)PDC鉆頭每顆齒的磨損趨勢(shì)。

聚晶金剛石鉆頭；切削齒；磨損；切削參數(shù)；零點(diǎn)遍歷法

0 引言

聚晶金剛石復(fù)合片(polycrystalline diamond compact，PDC) 鉆頭自誕生以來(lái)，在世界油氣鉆井總進(jìn)尺中的比例已由1988年的5%增長(zhǎng)到現(xiàn)在的90%，其鉆進(jìn)效率已是二三十年前的5～10倍[1]。

從提升PDC鉆頭性能的角度講，PDC鉆頭最理想的工作狀態(tài)是切削齒的均勻磨損，即符合等磨損布齒原則。然而該原則至今沒有真正實(shí)現(xiàn)，其主要原因是切削齒的均勻磨損設(shè)計(jì)需要精確計(jì)算鉆進(jìn)過(guò)程中每顆切削齒的切削載荷，而切削載荷是關(guān)于巖性、機(jī)械鉆速和鉆頭轉(zhuǎn)速等參數(shù)以及切削弧長(zhǎng)、切削面積等切削參數(shù)的函數(shù)。目前，國(guó)內(nèi)外的研究多集中在未磨損齒的切削參數(shù)求解，如采用微積分對(duì)切削參數(shù)進(jìn)行數(shù)值計(jì)算[2]；采用計(jì)算機(jī)模擬切削齒的真實(shí)破巖過(guò)程，通過(guò)切削齒與巖石間的布爾減運(yùn)算求解切削參數(shù)[3]；建立切削齒與巖石的三維模型并離散成3D網(wǎng)格，通過(guò)有限元仿真求得切削參數(shù)[4-5]。對(duì)磨損PDC鉆頭切削齒的切削參數(shù)研究還未見相關(guān)報(bào)道，這勢(shì)必導(dǎo)致無(wú)法準(zhǔn)確預(yù)測(cè)切削齒的受力變化和磨損趨勢(shì)，也就無(wú)法從根本上解決磨損不均問(wèn)題。

本文根據(jù)PDC鉆頭的實(shí)際鉆進(jìn)情況，考慮切削齒的空間位置以及先后切削順序，分析了切削齒的實(shí)際切削參數(shù)，提出了一種求解磨損齒切削參數(shù)的新方法。

1 磨損齒幾何模型的建立

圖1 切削齒磨損的主要形式Fig.1 The main wear type for PDC cutters

圖2 磨損齒幾何模型Fig.2 Geometry modeling for wear cutters

PDC切削齒最具代表性的磨損形態(tài)為弦截圓，如圖1所示[6]。結(jié)合IADC (international association of drilling contractors) 制定的固定切削齒磨損分級(jí)系統(tǒng)中的磨損模型，將磨損齒抽象為圖2中的幾何模型，其中，實(shí)線代表了磨損齒的齒刃邊線，齒刃邊線上的任意一點(diǎn)坐標(biāo)如下：圓弧段齒刃邊線的一點(diǎn)的坐標(biāo)(x4，y4，z4)=(rcsinφ，0，rccosφ)；直線段齒刃邊線(垂直于切削齒包絡(luò)線的法線，即圖2中的Z4軸)上的一點(diǎn)(x4，y4，z4)=((rc-Hd)tanφ，0，rc-Hd) ，其中，rc為切削齒半徑；φ為齒面輪廓線上任一點(diǎn)與齒心的連線和Z4軸的夾角；Hd為磨損高度。

2 磨損齒切削參數(shù)分析

為計(jì)算磨損齒的切削參數(shù)，將PDC鉆頭上的切削齒向鉆頭軸截面Mi內(nèi)投影，并取出5顆相鄰齒，從左至右的編號(hào)依次為i-2、i-1、i、i+1和i+2，如圖3所示。經(jīng)統(tǒng)計(jì)分析，第i顆磨損齒的實(shí)際切削參數(shù)與軸截面Mi內(nèi)的8條齒刃邊線投影相關(guān)。PDC鉆頭當(dāng)前鉆進(jìn)過(guò)程中，第i、i-1、i-2、i+1、i+2顆磨損齒的齒刃邊線投影的編號(hào)分別為p、p-1、p-2、p+1、p+2；前一圈鉆進(jìn)過(guò)程中，第i、i-1和i+1顆磨損齒的齒刃邊線投影的編號(hào)分別為p′、p′-1、p′+1。

圖3 PDC鉆頭磨損齒切削巖石示意圖Fig.3 Rock breaking of wear cutters in PDC bit

考慮齒的螺旋鉆進(jìn)，設(shè)第i顆齒上一圈鉆進(jìn)過(guò)程中，通過(guò)軸截面Mi時(shí)的標(biāo)高為h0，PDC鉆頭每轉(zhuǎn)的進(jìn)尺記為δ，對(duì)齒刃邊線上的任意一點(diǎn)坐標(biāo)進(jìn)行坐標(biāo)變換，求解得到上述8條齒刃邊線投影[7]：

(1)

(2)

(3)

式中，m為不同鉆進(jìn)狀態(tài)下δ的差值系數(shù)；n為齒刃邊線投影編號(hào)；r為齒刃邊線上任一點(diǎn)在鉆頭軸截面中的軸向坐標(biāo)；h為齒刃邊線上任一點(diǎn)在鉆頭軸截面中的徑向坐標(biāo)；θ為切削齒的周向角；α為切削齒的前傾角；β為切削齒的側(cè)轉(zhuǎn)角；γ為切削齒的法向角；下標(biāo)c代表關(guān)于齒心的參數(shù)；Rc為齒心在鉆頭軸截面中的軸向坐標(biāo)；Hc為齒心在鉆頭軸截面中的徑向坐標(biāo)。

當(dāng)n=p-2，p-1，p，p+1，p+2時(shí)，磨損齒處于當(dāng)前圈鉆進(jìn)狀態(tài)，m=1；當(dāng)n=p′-1，p′+1，p′時(shí)，磨損齒處于前一圈鉆進(jìn)狀態(tài)，m=0。

根據(jù)式(1)～式(3)，繪制出了磨損齒螺旋鉆進(jìn)情況下的8條齒面邊線投影，如圖 4所示，圖中陰影區(qū)域面積為第i顆磨損齒的切削截面面積，虛線為第i顆磨損齒的切削弧。切削截面繞鉆頭中心軸旋轉(zhuǎn)一周所形成的切削環(huán)體積為切削體積。

圖4 磨損齒的切削參數(shù)Fig.4 The cutting parameters of wear cutters

3 磨損齒切削參數(shù)求解

3.1 零點(diǎn)遍歷法的提出

通過(guò)對(duì)大量切削截面的研究分析得到：不管是磨損齒還是未磨損齒，齒刃邊線上φ=0的點(diǎn)均為關(guān)鍵點(diǎn)，如圖5中虛線上的圓點(diǎn)所示，稱這些圓點(diǎn)為“零點(diǎn)”。

圖5 切削齒上的零點(diǎn)Fig.5 Zero points of cutters

由圖2可得，零點(diǎn)與齒心的連線為切削齒包絡(luò)線的法線(圖2中的Z4軸)，而法線是垂直于巖石表面的，因此零點(diǎn)必定參與破碎巖石，同時(shí)也必被切削弧長(zhǎng)所包含?；诖?，提出一種從零點(diǎn)開始并在齒刃上進(jìn)行雙向遍歷[8-9]求解切削參數(shù)的方法——零點(diǎn)遍歷法。零點(diǎn)遍歷法的遍歷路徑如圖 6所示，從零點(diǎn)開始，沿齒刃向兩邊遍歷，直到遍歷路徑圍成環(huán)后終止，再根據(jù)所得的遍歷路徑計(jì)算切削參數(shù)。這種遍歷過(guò)程與切削參數(shù)的求解均通過(guò)MATLAB編碼實(shí)現(xiàn)。

圖6 磨損齒的遍歷路徑Fig.6 Traversal path of a wear cutter

3.2 切削參數(shù)具體求解步驟

切削參數(shù)的具體求解步驟如下：

(1)由零點(diǎn)X0開始，沿第p條齒刃邊線投影向左遍歷，求出與第p條齒刃邊線投影最先相交的點(diǎn)X1與齒刃邊線N1。

(2)由零點(diǎn)X0開始，沿第p條齒刃邊線投影向右遍歷，求出與第p條齒刃邊線投影最先相交的點(diǎn)X2與齒刃邊線N2。

(3)點(diǎn)X1沿N1號(hào)齒刃邊線投影向右遍歷，求出與N1號(hào)齒刃邊線投影最先相交的點(diǎn)X3和齒刃邊線N3。

(4)點(diǎn)X2沿N2號(hào)齒刃邊線投影向左遍歷，求出與N2號(hào)齒刃邊線投影最先相交的點(diǎn)X4和齒刃邊線N4。

(5)若N3與N4不是同一齒刃邊線的投影，則令X1=X3，N1=N3，X2=X4，N2=N4，返回到步驟(3)繼續(xù)計(jì)算；若N3與N4是同一齒刃邊線的投影，則終止遍歷，并記錄遍歷路徑。

(6)將遍歷路徑離散成K個(gè)點(diǎn)，構(gòu)建K邊形，并求出K個(gè)點(diǎn)在鉆頭軸截面內(nèi)的坐標(biāo)(xi，yi)，i=1，2，…，K。

(7)求解切削截面積A、切削體積V，分別為

(4)

V=2πRA

(5)

(6)

式中，R為K邊形的形心到鉆頭軸心的距離。

步驟(1)～步驟(5)與圖6中的序號(hào)對(duì)應(yīng)，為遍歷過(guò)程；步驟(6)～步驟(9)是對(duì)遍歷過(guò)程的數(shù)據(jù)處理，以用于求解切削參數(shù)。

在上述求解過(guò)程中，N1、N2、N3和N4代表了7條齒刃邊線投影p-1、p-2、p+2、p+1、p′、p′-1和p′+1中的某一條。K邊形的邊數(shù)影響了切削參數(shù)的求解精度，K邊形的邊越多，求解越精確。在下一節(jié)的計(jì)算中，K均取1000。

3.3 求解結(jié)果與分析

目前已有的切削參數(shù)求解算法的對(duì)象均是未磨損齒。為了便于將零點(diǎn)遍歷法與已有算法進(jìn)行對(duì)比分析，特將零點(diǎn)遍歷法的求解范圍拓展至未磨損齒切削參數(shù)的求解領(lǐng)域，此時(shí)式(3)中的Hdn設(shè)置為0，切削參數(shù)的具體求解步驟不變。

3.3.1 未磨損齒切削參數(shù)求解結(jié)果

以MATLAB程序設(shè)計(jì)平臺(tái)為基礎(chǔ)，將零點(diǎn)遍歷法編譯為MATLAB語(yǔ)言，并對(duì)某五刀翼PDC鉆頭上未磨損齒的切削弧長(zhǎng)和切削體積進(jìn)行求解。為便于比較，同時(shí)采用文獻(xiàn)[4]中的PDC鉆頭鉆進(jìn)仿真數(shù)字化技術(shù)對(duì)上述參數(shù)進(jìn)行了求解。如表 1所示，2種求解方法所得的切削弧長(zhǎng)相對(duì)誤差為1.13%，切削體積相對(duì)誤差為2.00%，誤差較小，驗(yàn)證了零點(diǎn)遍歷法求解未磨損齒切削參數(shù)的正確性。

3.3.2 磨損齒切削參數(shù)求解結(jié)果

采用零點(diǎn)遍歷法及Pro/E建模法對(duì)鉆頭上不同部位磨損齒的切削截面形狀及面積進(jìn)行分析與求解。為保證Pro/E建模法的可行性，將切削齒的前傾角和側(cè)轉(zhuǎn)角均設(shè)為0，在Pro/E的草繪模塊輸入相應(yīng)的尺寸位置參數(shù)，繪制切削截面并進(jìn)行測(cè)量，所得結(jié)果如表 2所示，2種算法求得的切削截面形狀一致，切削截面面積的相對(duì)誤差為0.035%。但Pro/E建模法只能對(duì)單齒進(jìn)行手動(dòng)操作計(jì)算，不能實(shí)現(xiàn)快速科學(xué)的計(jì)算。切削截面求出后，切削弧長(zhǎng)、切削體積等其他切削參數(shù)的求解就迎刃而解，因此，零點(diǎn)遍歷法可用于磨損齒切削參數(shù)的求解。

3.3.3 求解效率和精度分析

采用零點(diǎn)遍歷法對(duì)PDC鉆頭上切削齒的切削參數(shù)進(jìn)行求解，所需的求解時(shí)間約為5min。UG/Open二次開發(fā)技術(shù)和PDC鉆頭鉆進(jìn)仿真技術(shù)的求解時(shí)間分別為0.36 min和72 min[5]。零點(diǎn)遍歷法中，離散的K邊形邊數(shù)影響計(jì)算精度，采用直徑為19 mm的圓進(jìn)行測(cè)試，當(dāng)離散點(diǎn)數(shù)為628，即用628邊形面積代表圓面積時(shí)，計(jì)算精度達(dá)到1×10-8m，而零點(diǎn)遍歷法中每條齒刃邊線投影的離散點(diǎn)數(shù)為1000。PDC鉆頭鉆進(jìn)仿真技術(shù)采用的是有限元仿真技術(shù)，網(wǎng)格越大、精度越低，實(shí)驗(yàn)采用的網(wǎng)格尺寸在1mm左右；UG/Open二次開發(fā)技術(shù)的理論精度等級(jí)為1×10-6m[5]。因此，零點(diǎn)遍歷法與上述兩種方法相比，計(jì)算精度成倍增加，同時(shí)保證了計(jì)算時(shí)間不至于過(guò)長(zhǎng)。此外，零點(diǎn)遍歷法的精度能夠人為調(diào)節(jié)，可以很好地平衡精度與效率兩者之間的關(guān)系，適應(yīng)不同的計(jì)算需求。

表1 PDC切削齒的切削弧長(zhǎng)和切削體積

表2 磨損齒的切削截面形狀及面積

4 應(yīng)用實(shí)例

4.1 鉆頭受力分析

根據(jù)某A型鉆頭在在渤海油田KL20-1-1井鉆穿泥質(zhì)砂巖層段的實(shí)際鉆井記錄，鉆頭的每轉(zhuǎn)平均切深為4 mm，地層可鉆性極值為3.2。結(jié)合PDC鉆頭受力模型[10]，對(duì)鉆頭進(jìn)行受力分析，并將與實(shí)際鉆進(jìn)中的PDC鉆頭受力狀況進(jìn)行對(duì)比。如表 3所示，鉆壓W、扭矩T的計(jì)算誤差分別為0.95%和4.67%，驗(yàn)證了該方法可用于鉆頭的受力分析。

表3 A型鉆頭的力學(xué)參數(shù)比較

4.2 磨損趨勢(shì)預(yù)測(cè)

結(jié)合PDC鉆頭磨損模型[11]，可對(duì)切削齒的磨損趨勢(shì)做出預(yù)測(cè)。圖 7所示為某PDC鉆頭上的32顆切削齒在不同時(shí)間(第7號(hào)齒磨損4 mm時(shí)的時(shí)間記為1)下的磨損趨勢(shì)預(yù)測(cè)?？梢钥闯觯S著磨損時(shí)間的延長(zhǎng)，磨損量不斷增大；隨著齒號(hào)的增大，磨損量呈先增大、后減小的趨勢(shì)，且處于PDC鉆頭冠頂部位的切削齒(7～12號(hào))的磨損量最大，這與PDC鉆頭在真實(shí)鉆進(jìn)中冠頂部位齒磨損嚴(yán)重[6]的現(xiàn)象相符(圖8)。因此，基于零點(diǎn)遍歷法的切削齒磨損趨勢(shì)預(yù)測(cè)是可行的，且能夠量化不同時(shí)間段內(nèi)的磨損量，為等磨損布齒的設(shè)計(jì)與評(píng)價(jià)提供依據(jù)。

圖7 鉆頭上各個(gè)齒的磨損趨勢(shì)預(yù)測(cè)Fig.7 Wear trend predicting for each cutter in PDC bit

圖8 PDC鉆頭磨損嚴(yán)重的冠頂部位齒Fig.8 High wear cutters on the nose of PDC profile

5 結(jié)論

(1)構(gòu)建的齒刃邊線投影方程能很好地表述磨損齒的切削弧和切削截面等切削參數(shù)的特征。

(2)提出的零點(diǎn)遍歷法能夠用于磨損前后切削齒切削參數(shù)的求解。零點(diǎn)遍歷法與Pro/E建模法的相對(duì)誤差為0.035%，與鉆進(jìn)仿真數(shù)字化技術(shù)的相對(duì)誤差小于2%。

(3)零點(diǎn)遍歷法的求解精度為1×10-8m，計(jì)算精度能夠人為調(diào)節(jié)，適應(yīng)不同的計(jì)算需求；采用零點(diǎn)遍歷法，并結(jié)合切削齒受力和磨損模型，能預(yù)測(cè)鉆頭上每顆切削齒的磨損趨勢(shì)，為等磨損布齒的設(shè)計(jì)與評(píng)價(jià)提供依據(jù)。

[1] 左汝強(qiáng)．國(guó)際油氣井鉆頭進(jìn)展概述(三)PDC鉆頭發(fā)展進(jìn)程及當(dāng)今態(tài)勢(shì)(上)[J]．探礦工程(巖土鉆掘工程)，2016, 43(3)：1-8． ZUO Ruqiang. International Advancement of Drilling Bits for Oil and Gas Well(3)—PDC Bits Progress and Present Trend(Ⅰ)[J]. Exploration Engineering(Rock & Soil Drilling and Tunneling), 2016, 43(3):1-8.

[2] 鄒德永，王家駿，于金平，等．定向鉆井PDC鉆頭切削參數(shù)計(jì)算方法[J]．石油鉆采工藝，2014, 36(5)：5-9． ZOU Deyong, WANG Jiajun, YU Jinping, et al. Computing Method for Cutting Parameters of PDC Bit in Directional Drilling[J]. Oil Drilling & Production Technology, 2014, 36(5):5-9.

[3] 況雨春，陳玉中，屠俊文，等．基于 UG/OPEN 的 PDC鉆頭切削參數(shù)仿真方法[J]．石油鉆探技術(shù)，2014, 42(4)：111-115． KUANG Yuchun, CHEN Yuzhong, TU Junwen, et al. Simulation of Cutting Parameters of PDC Bit Based on UG/OPEN[J]. Petroleum Drilling Techniques, 2014, 42(4):111-115.

[4] 任海濤，楊迎新，陳煉，等．PDC鉆頭鉆進(jìn)仿真系統(tǒng)數(shù)字化方法研究[J]．西南石油大學(xué)學(xué)報(bào)(自然科學(xué)版)，2010, 32(5)：150-154． REN Haitao, YANG Yingxin, CHEN Lian, et al. Research on the Digital Method of the PDC Bit Drilling Progress Simulation System[J]. Journal of Southwest Petroleum University( Science & Technology Edition), 2010, 32(5):150-154.

[5] QU J L, LIU X D. Optimization Research on the Geometrical Parameters of PDC Core Bits Based on ANSYS[J]. Applied Mechanics & Materials, 2013, 275/277:803-807.

[6] 張富曉，黃志強(qiáng)，周已．PDC鉆頭切削齒失效分析[J]．石油礦場(chǎng)機(jī)械，2015, 44(9)：44-49． ZHANG Fuxiao, HUANG Zhijiang, ZHOU Yi. Failure Analysis of PDC Bit Cutter[J]. Oil Field Equipment, 2015, 44(9):44-49.

[7] 楊云志，王小平，張維中，等．基于投影法的過(guò)渡曲面構(gòu)造方法[J]．中國(guó)機(jī)械工程，2014, 25(6)：780-783． YANG Yunzhi, WANG Xiaoping, ZHANG Weizhong, et al. Construction Transition Surface Based on Curve Projection[J]. China Mechanical Engineering, 2014, 25(6):780-783.

[8] 張赤斌，王興松．基于蟻群算法的完全遍歷路徑規(guī)劃研究[J]．中國(guó)機(jī)械工程，2008, 19(16)：1945-1949． ZHANG Chibin, WANG Xingsong. Complete Coverage Path Planning Based on Ant Colony Algorithm[J]. China Mechanical Engineering, 2008, 19(16):1945-1949.

[9] ZHANG H, WANG P, GU X, et al. Least Operation Traversal Method Applied in Optimization of Logic Circuits[C]//2009 IEEE 8th International Conference on ASIC. Changsha, 2009:887-890.

[10] 梁爾國(guó)，李子豐，鄒德永．PDC鉆頭綜合受力模型的試驗(yàn)研究[J]．巖土力學(xué)，2009, 30(4)：938-942． LIANG Erguo, LI Zifeng, ZOU Deyong. Experimental Research on Integrated Mechanical Model of PDC Bit[J]. Rock and Soil Mechanics, 2009, 30(4):938-942.

[11] 曲冠政，曲大偉，郭瑞，等．PDC鉆頭復(fù)合片磨損速度模型研究[J]．復(fù)雜油氣藏，2013, 6(1)：62-64． QU Guanzheng, QU Dawei, GUO Rui, et al. Study on a Wear Rate Model of PDC Composite Cutters[J]. Complex Hydrocarbon Reservoirs, 2013, 6(1):62-64.

(編輯 張 洋)

Calculation of Cutting Parameters of Cutters of PDC Bit under Wear Conditions

MA Yachao1HUANG Zhiqiang1LI Qin1PENG Shijin2XIE Dou1

1.Mechatronic Engineering College,Southwest Petroleum University,Chengdu,610500 2.Laboratory and Equipment Management,Southwest Petroleum University,Chengdu,610500

For the problems of cutting parameters calculation and cutter layout design under wearing conditions might not meet the requirements of equal wear principles very well, a geometric model of waring cutters was built. Zero point traversing method used to calculate the cutting parameters was proposed. Then, cutting arc length, cutting area and cutting volume were calculated by zero point traversing method based on MATLAB software. The results reveal that zero point traversing method may be used to calculate the cutting parameters of wearing cutters and unworn cutters, and the calculation error is less than 2%. Combined with the force model and the wearing model, the wear tendency of the cutters in PDC bit may be predicted.

polycrystalline diamond compact(PDC) bit; cutting tooth; wear; cutting parameter; zero point traversing method

趙林峰，男，1979年生。合肥工業(yè)大學(xué)汽車與交通工程學(xué)院副教授。主要研究方向?yàn)檐囕v動(dòng)力學(xué)與控制。發(fā)表論文20余篇。E-mail:13856952050@163.com。徐 磊，男，1992年生。合肥工業(yè)大學(xué)汽車與交通工程學(xué)院碩士研究生。陳無(wú)畏，男，1951年生。合肥工業(yè)大學(xué)汽車與交通工程學(xué)院教授、博士研究生導(dǎo)師。

2016-06-02

TE921

10.3969/j.issn.1004-132X.2017.08.016