李　瑋， 紀(jì)照生, 董智煜， 李卓倫， 李　悅

(1.東北石油大學(xué) 石油工程學(xué)院，黑龍江 大慶　163318; 2.大慶油田有限責(zé)任公司第三采油廠,黑龍江 大慶　163318)

?

基于重整化方法的沖擊載荷下巖石振動分析

李瑋1， 紀(jì)照生1, 董智煜2， 李卓倫1， 李悅1

(1.東北石油大學(xué) 石油工程學(xué)院，黑龍江 大慶163318; 2.大慶油田有限責(zé)任公司第三采油廠,黑龍江 大慶163318)

旋轉(zhuǎn)沖擊鉆井是當(dāng)前鉆進(jìn)深井硬地層應(yīng)用較多且效果顯著的一種高效破巖技術(shù)?；谡駝訉W(xué)理論，在考慮巖石重力的情況下，建立鉆頭沖擊載荷下巖石振動響應(yīng)的數(shù)學(xué)模型，采用重整化方法對其求解，利用MATLAB軟件分析各參數(shù)對鉆頭高頻沖擊下巖石振動的影響，并進(jìn)行了室內(nèi)實驗對理論分析結(jié)果進(jìn)行了驗證。結(jié)果顯示：沖擊載荷下巖石的運動為兩個簡諧振動的互擾振動；巖石的密度越小，剛度越小，沖擊力越大，巖石的振動幅度越大，振動速度也越快；巖石存在一個固有頻率，沖擊載荷沖擊頻率越接近這個固有頻率，巖石的振動越劇烈。當(dāng)鉆頭高頻振動頻率和巖石固有頻率相等時，巖石振動幅度和振動速度最大，即達(dá)到了共振。沖擊載荷下巖石響應(yīng)的分析對于揭示動載作用下巖石的破碎機(jī)理，指導(dǎo)沖擊工具的設(shè)計具有重要意義。

旋轉(zhuǎn)沖擊鉆井；巖石振動；幅頻特性；共振

隨著石油勘探開發(fā)的不斷深入，深井超深井的比例逐年增加，復(fù)雜地層與難鉆地層鉆遇率升高，旋沖鉆井作為解決當(dāng)前鉆井作業(yè)難題較為有效的高效破巖技術(shù)越來越受到科研人員的重視[1-6]。

沖擊載荷下巖石的響應(yīng)機(jī)制和破碎機(jī)理作為沖擊鉆井技術(shù)的重中之重已經(jīng)被部分學(xué)者所研究[7-10]。金解放等[11]利用巖石動靜組合加載SHPB試驗裝置對不同靜載砂巖試件進(jìn)行循環(huán)沖擊試驗，對不同靜載作用下巖石試件的應(yīng)力狀態(tài)進(jìn)行分析，探索具有一定靜載的巖石在循環(huán)沖擊作用下的破裂機(jī)理。宋義敏等[12]以可調(diào)速落錘沖擊試驗機(jī)作為試驗加載裝置，通過搭建高速數(shù)據(jù)采集系統(tǒng)，對巖石Ⅰ型裂紋在沖擊載荷作用下的動態(tài)斷裂進(jìn)行試驗研究。李瑋等[13]根據(jù)機(jī)械振動原理建立了巖石穩(wěn)態(tài)振動響應(yīng)計算模型，根據(jù)旋轉(zhuǎn)剪切和高頻振動沖擊的破巖試驗結(jié)果，分析了機(jī)械參數(shù)和振動參數(shù)對破巖效果的影響。趙伏軍等[14]根據(jù)壓頭侵入巖石的斷裂特征，確定了裂紋擴(kuò)展計算公式，分析了沖擊-靜壓切削組合破巖模式下沖擊間距(頻率)對破巖效果的影響。

作者基于振動學(xué)理論，考慮巖石重力，建立鉆頭高頻沖擊下巖石振動響應(yīng)的數(shù)學(xué)模型，采用重整化方法對其求解，利用MATLAB軟件分析各參數(shù)對鉆頭高頻沖擊下巖石振動的影響，并通過室內(nèi)實驗對其進(jìn)行了驗證。沖擊載荷下巖石響應(yīng)的分析對于揭示動載作用下巖石的破碎機(jī)理，指導(dǎo)沖擊工具的設(shè)計具有重要意義。

1　旋轉(zhuǎn)沖擊鉆井技術(shù)

旋轉(zhuǎn)沖擊鉆井技術(shù)實質(zhì)上就是將沖擊鉆進(jìn)的單次破碎巖石作用和旋轉(zhuǎn)鉆進(jìn)的連續(xù)破碎巖石作用相結(jié)合的一種鉆井方法，其主要方式是在鉆頭上聯(lián)接一個專用的沖擊器。鉆進(jìn)時，鉆頭在常規(guī)鉆壓和扭矩作用的基礎(chǔ)上同時通過沖擊器給鉆頭施加一定頻率的沖擊載荷，井底鉆頭就在沖擊和旋轉(zhuǎn)共同作用下破碎巖石，進(jìn)行鉆進(jìn)(見圖1)。交變的沖擊動載，使鉆頭上的比載荷增加，加載速度提高，巖石產(chǎn)生形變所需時間縮短，變形速度增大，被沖擊點還來不及對作用力重新分配，應(yīng)力便很快接近或超過強(qiáng)度極限，使巖石脆性增加，塑性下降。被沖擊巖體容易破碎形成坑穴和產(chǎn)生剪切體，有利于體積破碎的發(fā)生，從而可以提高鉆井速度[15]。

圖1　旋轉(zhuǎn)沖擊破巖示意圖Fig.1 Rotary impact rock breaking schematic

2　巖石振動方程的建立

雖然在現(xiàn)場實踐中，井底鉆頭是在沖擊和旋轉(zhuǎn)共同作用下破碎巖石。但為了使獲得的高頻沖擊下巖石響應(yīng)機(jī)制更具有普遍性，以巖石僅受沖擊載荷為例進(jìn)行分析。

2.1物理模型

2.2數(shù)值模型

將巖石區(qū)域分割，單獨取出一個彈簧微元進(jìn)行分析?，F(xiàn)以u1，u2……uN-1，uN表示各個小球的坐標(biāo)(見圖4)。

圖4　等效模型微元Fig.4 Micro unit of equivalent model

易得：

(1)

(2)

則式(1)變?yōu)椋?/p>

(3)

(4)

式中：ω即為彈簧固有振動頻率。

對xn(t)求二階導(dǎo)得

(5)

將式(5)代入式(3)并整理得：

(6)

2.3利用重整化方法進(jìn)行求解

在量子場論發(fā)展的早期，人們發(fā)現(xiàn)許多微擾展開的高階項的計算結(jié)果含有發(fā)散項，這些發(fā)散項和量子場論所描述的系統(tǒng)具有無窮自由度有關(guān)。于是他們在拉氏量中引進(jìn)有限數(shù)目的項來抵消這些發(fā)散項，這種方法即為重整化方法。重整化以后的理論不發(fā)散，而且能與實驗進(jìn)行比較，從而可以得到相應(yīng)的正確觀測量[16]。

由式(6)可知，本文所建立模型的自由度為2N，當(dāng)N趨向于無窮大時，采用常規(guī)解法不能求解。利用重整化方法，可引入有限數(shù)目的項，來逐漸減少系統(tǒng)的自由度，最終獲得模型的解。

(7)

上述過程重復(fù)q次得：

(8)

令bq=2yq得：yq=2yq-12-1

對比cos(2α)=2cos2α-1的形式，若令y1=cosα則易得：

(9)

此時運動方程變?yōu)椋?/p>

(10)

結(jié)合式(8)～式(10)，可解得:

(11)

(12)

計算可得其通解為：

(13)

式中：F0為施加的沖擊載荷的幅值，ω1為施加沖擊載荷的沖擊頻率，φ1為施加沖擊載荷的相位角，A，φ為初始振幅和相位，反應(yīng)的是沖擊載荷下巖石初始能量和初始位置。

式(13)即為鉆頭高頻振動沖擊下巖石響應(yīng)的運動方程。觀察其特點可知，巖石響應(yīng)為兩個簡諧振動的互擾振動。

3　巖石振動響應(yīng)的影響因素分析

鉆頭高頻沖擊下巖石的振動情況直接關(guān)系到破巖效果的好壞，因此應(yīng)用MATLAB軟件對鉆頭高頻振動沖擊下巖石的振動方程進(jìn)行分析，確定鉆頭下巖石的密度、剛度、沖擊力的大小、沖擊的頻率等參數(shù)對巖石振動的影響規(guī)律，為現(xiàn)場沖擊鉆井提供有益的指導(dǎo)。

(1) 沖擊載荷大小對巖石振動的影響

圖5給出了不同大小鉆頭沖擊力作用下巖石表面振動位移、振動速度隨時間變化的關(guān)系曲線。可知，鉆頭高頻振動沖擊下，巖石呈現(xiàn)較明顯的擾動特征。在相同激勵頻率下，巖石振動的最大位移和振動速度隨動態(tài)激振力的增加而增加。鉆頭的沖擊力越大，巖石能夠達(dá)到的最大振動位移和最大振動速度越大。

圖5　不同鉆頭沖擊力下巖石振動位移、振動速度隨時間的變化Fig.5 Vibration displacement and vibration velocity of rock under different bit impact force

(2) 沖擊頻率對巖石振動的影響及共振頻率影響因素分析

圖6中不同鉆頭沖擊頻率作用下巖石表面振動位移、振動速度隨時間變化的關(guān)系曲線表明，沖擊力大小相同時，巖石振動的最大位移和振動速度均隨著沖擊頻率的增大先增大后變小。經(jīng)驗可知，從增大到減小的轉(zhuǎn)折點，即為巖石的固有頻率。當(dāng)鉆頭的沖擊頻率與巖石的固有頻率相等時，鉆頭便與巖石產(chǎn)生共振。

圖6　不同鉆頭沖擊頻率下巖石振動位移、振動速度隨時間的變化Fig.6 Vibration displacement and vibration velocity of rock under different bit impact frequency

共振破巖是高頻沖擊破巖的最佳狀態(tài)，此時沖擊載荷的利用率最高，破巖效果最好，鉆進(jìn)速度最快。共振頻率即巖石的固有頻率，是巖石的固有屬性。由巖石振動方程知，巖石固有頻率與巖石的密度和剛度有關(guān)，關(guān)系曲線(見圖7、圖8)。由圖知，巖石固有頻率與巖石剛度的開方成正比關(guān)系，與巖石密度的開方成反比關(guān)系。

圖7　巖石固有頻率隨巖石剛度的變化Fig.7 Changes of natural frequency with the rock rigidity

圖8　巖石固有頻率隨巖石密度的變化Fig.8 Changes of natural frequency with the rock density

(3) 巖石剛度對巖石振動的影響

圖9給出了不同剛度下巖石表面振動位移、振動速度隨時間變化的關(guān)系曲線。由圖9可知，巖石的剛度越小，巖石表面在鉆頭振動沖擊下能夠達(dá)到的最大振動位移和最大振動速度就越大。這是由于在沖擊力大小相同和巖石密度相同的條件下，巖石剛度越小，巖石越易產(chǎn)生變形，也越易產(chǎn)生振動。

圖9　不同巖石剛度下巖石振動位移、振動速度隨時間的變化Fig.9 Vibration displacement and vibration velocity of rock of different rigidity

(4) 巖石密度對巖石振動的影響

圖10給出了不同密度下巖石表面振動位移、振動速度隨時間變化的關(guān)系曲線?？芍?，巖石密度越小，越容易產(chǎn)生振動，相應(yīng)的振動位移也越大，越容易被破碎。

圖10　不同巖石密度下巖石振動位移、振動速度隨時間的變化Fig.10 Vibration displacement and vibration velocity of rock of different density

4　室內(nèi)軸向沖擊破巖試驗

為了證明理論模型及影響因素分析的正確性，進(jìn)行室內(nèi)軸向沖擊試驗。獲得不同沖擊力、沖擊頻率下不同巖性巖樣表面的反應(yīng)，驗證各種因素對鉆頭沖擊作用下巖石振動的影響規(guī)律。

4.1試驗裝置

試驗裝置為東北石油大學(xué)“高效鉆井破巖技術(shù)”研究室自主研發(fā)的“巖石振動激勵模擬測試系統(tǒng)”，見圖11。測試系統(tǒng)由主體機(jī)架、巖樣夾持機(jī)構(gòu)、鉆頭高頻振動發(fā)生裝置、鉆進(jìn)系統(tǒng)、旋轉(zhuǎn)裝置、鉆壓模擬控制系統(tǒng)、數(shù)據(jù)采集控制系統(tǒng)等部分構(gòu)成。其功能主要是在普通鉆進(jìn)過程中產(chǎn)生各種振動激勵，如低頻大振幅和鉆頭高頻小振幅的振動激勵，通過鉆進(jìn)參數(shù)對試驗樣品進(jìn)行振動特性采集及各參數(shù)分析。

4.2試驗方法

試驗選取2塊灰質(zhì)紅砂巖，1塊泥質(zhì)粉砂巖進(jìn)行實驗，其物性參數(shù)見表1。設(shè)置系統(tǒng)不旋轉(zhuǎn)，只振動。振動系統(tǒng)對鉆頭的沖擊頻率為800～2 100 Hz，沖擊力為2 000～8 000 N。

表1　巖石試樣的物性參數(shù)

圖11　巖石振動激勵模擬測試系統(tǒng)Fig.11 Rock vibration excitation simulation test system

4.3實驗結(jié)果分析

分別對不同彈性模量、相同密度的1號灰質(zhì)粉砂巖和2號灰質(zhì)粉砂巖進(jìn)行2組沖擊試驗。由于材料的剛度通常用彈性模量來衡量，于是得到剛度對巖石振動的影響規(guī)律。由圖12可知，巖石的振動位移隨著剛度的增大而減小。這是由于巖石剛度越小，越容易產(chǎn)生變形，在振動沖擊下，巖石表面位移越大。

由圖13可知，在其他參數(shù)相同的情況下，巖石的最大振動位移隨著沖擊力的增大而增大。但是隨著沖擊力的增大，最大振動位移的增長速度變慢。這是因為當(dāng)沖擊力比較小時，巖石空隙壓實所產(chǎn)生的位移所占的比例比較大。隨著沖擊力的增大，巖石逐漸壓實，最大位移增長速度逐漸變慢。

由圖14可知,隨著沖擊激勵頻率增大，巖石最大振動位移先增大后減小，在沖擊激勵頻率為1 500 Hz左右時達(dá)到最大，即實現(xiàn)了共振。實驗結(jié)果與理論分析結(jié)果一致。

5　結(jié)　論

(1) 考慮巖石重力，建立了鉆頭高頻沖擊下巖石模型，應(yīng)用重整化方法，求解出了巖石運動方程。鉆頭高頻沖擊載荷下巖石的響應(yīng)為兩個簡諧振動的互擾振動。

(2) 巖石的密度、巖石的剛度、鉆頭沖擊力的大小、沖擊的頻率等對鉆頭下巖石振動有影響。巖石的密度越小，剛度越小，鉆頭高頻沖擊的沖擊力越大，巖石的振動幅度越大，振動速度也越快。

(3) 巖石存在一個固有頻率，其與巖石剛度正相關(guān)，與巖石密度負(fù)相關(guān)。鉆頭高頻振動沖擊頻率越接近巖石的固有頻率，巖石的振動越劇烈。當(dāng)鉆頭高頻振動頻率和巖石固有頻率相等時，巖石振動幅度和振動速度最大，即達(dá)到了共振。

[1] 常德玉，李根生，沈忠厚，等. 深井超深井井底應(yīng)力場[J]. 石油學(xué)報，2011，32(4)：697-703.

CHANG Deyu, LI Gensheng, SHEN Zhonghou, et al. The stress field of bottom hole in deep and ultra-deepwells[J]. Acta Petrolei Sinica, 2011, 32(4): 697-703.

[2] 祝效華，湯歷平，童華. 高頻扭轉(zhuǎn)沖擊鉆進(jìn)的減振與提速機(jī)理研究[J]. 振動與沖擊，2012，31(20)：75-78.

ZHU Xiaohua, TANG Liping, TONG Hua. Rock breaking mechanism of a high frequency torsional impact drilling[J]. Journal of Vibration and Shock, 2012, 31(21): 75-78.

[3] LI Wei, YAN Tie, LI Siqi, et al. Rock fragmentation mechanisms and an experimental study of drilling tools during high-frequency harmonic vibration[J]. Petroleum Science, 2013, 10(2): 205-211.

[4] AGUIAR R R, WEBER H I. Development of a vibroimpact device for the resonance hammer drilling[C]//Proceedings of the XII International Symposium on Dynamic Problems of Mechanics. Ilhabela,Brazil, 2007.

[5] 熊繼有，錢聲華，嚴(yán)仁俊，等. 鉆井高效破巖新進(jìn)展[J]. 天然氣工業(yè)，2004，24(4)：27-29.

XIONG Jiyou, QIAN Shenghua, YAN Renjun, et al. New progress of high efficiency of drilling[J]. Natural Gas Industry, 2004, 24(4): 27-29.

[6] 雷鵬，倪紅堅，王瑞和，等．自激振蕩式旋轉(zhuǎn)沖擊鉆井工具水力元件性能分析與優(yōu)化[J]．振動與沖擊，2014，33(19)：175-180.

LEI Peng, NI Hongjian, WANG Ruihe, et al.Performance analysis and optimization for hydraulic components of self-oscillating rotary impact drilling tool[J]. Journal of Vibration and Shock, 2014, 33(19): 175-180.

[7] LI Siqi, YAN Tie, LI Wei, et al. Modeling of vibration response of rock by harmonic impact[J]. Journal of Natural Gas Science and Engineering, 2015, 23: 90-96.

[8] 文平，陳波，雷巨鵬，等．液動沖擊旋轉(zhuǎn)鉆井技術(shù)在玉門青西油田的應(yīng)用[J]．天然氣工業(yè)，2004，24(9)：64-67.

WEN Ping, CHEN Bo, LEI Jupeng, et al. Application of hydraulic impact rotary drilling technology in Qingxi Oilfield[J]. Natural Gas Industry, 2004, 24(9): 64-67.

[9] 姚曄，徐寶富，陸敏恂. 低頻大振幅液壓沖擊器的建模和分析方法[J]. 同濟(jì)大學(xué)學(xué)報，1999，27(6)：741-744.

YAO Ye, XU Baofu, LU Minxun. Modeling and analytic method of low frequency and large amplitude hydraulic shocker[J]. Journal of Tongji University,1999,27(6):741-744.

[10] WIERCIGROCH M, KRIVTSOV A M, WOJEWODA J. Vibrational energy transfer via modulated impacts for percussive drilling[J]. Journal of Theoretical and Applied Mechanics, 2008, 46(3): 715-726.

[11] 金解放，李夕兵，王觀石，等．循環(huán)沖擊載荷作用下砂巖破壞模式及其機(jī)理[J]．中南大學(xué)學(xué)報:自然科學(xué)版，2012，43(4)：1453-1460.

JIN Jiefang, LI Xibing, WANG Guanshi, et al. Failure modes and mechanisms of sandstone under cyclic impact loadings[J]. Journal of Central South University:Science and Technology,2012, 43(4): 1453-1460.

[12] 宋義敏，楊小彬，金璐，等．沖擊載荷作用下巖石Ⅰ型裂紋動態(tài)斷裂試驗研究[J]．振動與沖擊，2014，33(11):49-53.

SONG Yimin, YANG Xiaobin, JIN Lu, et al. Dynamic fracture test for rock I-type crack under impact load[J]. Journal of Vibration and Shock, 2014, 33(11): 49-53.

[13] 李瑋，閆鐵，張志超，等. 高頻振動鉆具沖擊下巖石響應(yīng)機(jī)理及破巖試驗分析[J]. 石油鉆探技術(shù)，2013，41(6)：25-28.

LI Wei, YAN Tie, ZHANG Zhichao, et al. Rock response mechanism and rock breaking test analysis for impact of high frequency vibration drilling tool[J]. Petroleum Drilling Techniques, 2013, 41(6): 25-28.

[14] 趙伏軍，李夕兵，馮濤，等. 動靜載荷耦合作用下巖石破碎理論分析及試驗研究[J]. 巖石力學(xué)與工程學(xué)報，2005，24(8)：1315-1320.

ZHAO Fujun, LI Xibing, FENG Tao, et al. Theoretical analysis and experimental study on the fracture of the rock under static and dynamic loads[J]. Chinese Journal of Rock Mechanics and Engineering, 2005, 24(8): 1315-1320.

[15] 陶興華，張建龍，曾義金. 石油旋沖鉆井技術(shù)研究及應(yīng)用[J]. 石油鉆采工藝，1998，20(2)：27-30.

TAO Xinghua, ZHANG Jianlong, ZENG Yijin. Study and application of rotary percussion drilling technology[J]. Oil Drilling & Production Technology, 1998, 20(2): 27-30.

[16] 柴立和. 重整化方法及其應(yīng)用的研究進(jìn)展[J]. 現(xiàn)代物理知識，2005，17(2)：32-35.

CHAI Lihe. Progress in reforming method and its applications[J]. Modern Physics, 2005, 17(2): 32-35.

Vibration analysis of rock under impact loads based on the renormalization method

LI Wei1, JI Zhaosheng1, DONG Zhiyu2, LI Zhuolun1, LI Yue1

(1. Petroleum Engineering School, Northeast Petroleum University, Daqing 163318, China;2. No.3 Oil Production Plant of Daqing Oilfield Company, Daqing 163318, China)

Currently rotary percussion drilling is the most used and efficient rock breaking technology in deep hard formations drilling and the effect is remarkable. Based on vibration theory and the consideration of the gravity of rocks, a mathematical model of rock vibration under impact loads was established. Renormalization method was used to get the solution and MATLAB software was used to analyze the influence of various parameters on the vibration of the rock. Finally, laboratory experiments were conducted to verify the analytical results. It is shown that the movement of rock under drill bit’s impact is mutual interference vibration between two harmonic vibrations. The vibration amplitude of the rock as well as the vibration velocity becomes larger with the decrease of rock density, the decrease of rock stiffness and the increase of the impact force. There is an inherent frequency of the rock. Vibration of the rock becomes more intense with the decrease of the difference between the impact frequency and the inherent frequency. When the high impact frequency equals to the inherent frequency, the vibration amplitude and vibration velocity reach the largest, which is the resonance. For revealing the crushing mechanism of the rock under dynamic loads and guiding the design of impact tools, analysis of the response of the rock under impact loads is of great significance.

rotary percussion drilling; rock vibration; amplitude-frequency characteristics; resonance

中國博士后科學(xué)基金第6次特別資助“微進(jìn)尺諧振沖擊器設(shè)計及其破巖實驗研究”(2013T60343);東北石油大學(xué)研究生創(chuàng)新科研項目“軸向沖擊破巖技術(shù)研究”(YJSCX2016-015NEPU)

2015-11-18修改稿收到日期：2016-01-19

李瑋 男，教授，博士生導(dǎo)師，1979年生E-mail: cyyping@sina.com

TE242

A

10.13465/j.cnki.jvs.2016.16.009