王維，彭枧明，王紅巖，陳晨，王馨靚，付偉根

(1. 吉林大學(xué) 建設(shè)工程學(xué)院，吉林 長(zhǎng)春，130026；2. 中國(guó)石油西部鉆探工程有限公司 國(guó)際合作處，新疆 烏魯木齊，830011；3. 大慶油田有限責(zé)任公司 水務(wù)公司，黑龍江 大慶，163000；4. 中水北方勘測(cè)設(shè)計(jì)研究有限責(zé)任公司，天津，300202)

天然氣水合物主要賦存于海底以下0～1 500 m、年代新、強(qiáng)度低的非固結(jié)松散沉積地層和高緯度大陸地區(qū)永凍土帶及水深100～250 m 以下極地陸架海[1]。三軸實(shí)驗(yàn)結(jié)果表明：實(shí)鉆獲得的天然氣水合物巖心強(qiáng)度較低，與天然軟黏土和凍土的巖心強(qiáng)度接近[2-4]。由此可以知道，在鉆井過(guò)程中，井底鉆頭在承受一定鉆壓及扭矩的條件下切削破碎區(qū)以下一定范圍內(nèi)的天然氣水合物地層，將類似受一定外力作用下的土，經(jīng)歷骨架顆粒的擠密和彎曲、骨架顆粒間連接形狀的扭曲和骨架顆粒間的錯(cuò)動(dòng)過(guò)程，部分骨架顆粒落入原來(lái)的孔隙中，發(fā)生塑性體積變形，從而導(dǎo)致孔隙度、飽和度、滲透率等重要儲(chǔ)層參數(shù)的變化[5]。在取心鉆井過(guò)程中，經(jīng)鉆頭擾動(dòng)的這部分地層物質(zhì)，一部分進(jìn)入巖心管形成巖心，一部分構(gòu)成了井周力學(xué)擾動(dòng)帶，其余被鉆井液帶走。巖心和井周地層受嚴(yán)重?cái)_動(dòng)后，由于孔隙度、飽和度等參數(shù)的變化導(dǎo)致測(cè)井波速等用于儲(chǔ)量計(jì)算的主要依據(jù)數(shù)據(jù)的有效性顯著較低，因此，天然氣水合物鉆井過(guò)程中鉆頭擾動(dòng)將不利于獲取較為精確的巖心和井周地層的波速。由天然氣水合物地層資源量計(jì)算公式可以得知，巖心孔隙度和聲波、密度、中子測(cè)井孔隙度等較小的相對(duì)誤差將引起計(jì)算的資源量出現(xiàn)很大的絕對(duì)誤差[6]。由此可知，井底擾動(dòng)區(qū)域及擾動(dòng)程度對(duì)井周地層的波速變化產(chǎn)生很大的影響，因此，研究一定鉆壓及轉(zhuǎn)速條件下井底及井周地層應(yīng)力分布狀態(tài)對(duì)于獲取天然氣水合物地層準(zhǔn)確信息具有重要意義。國(guó)際大洋鉆探計(jì)劃(ODP，ocean drilling project)采用的保壓取心器PCS(pressure core sampler)是目前國(guó)內(nèi)外使用最廣泛的保壓取樣器之一，該取心器的研制成功被認(rèn)為在天然氣水合物取心器方面取得了重大突破[7]。PCS 在ODP Leg124，Leg139，Leg141，Leg146，Leg164，Leg196，Leg201 和Leg204 等多個(gè)航次中進(jìn)行了應(yīng)用并成功取得天然氣水合物樣品[8]。在鉆井過(guò)程中，常規(guī)實(shí)驗(yàn)方法難以獲得井底及井周鉆頭附近地層應(yīng)力場(chǎng)及應(yīng)變場(chǎng)，所以，要想了解其擾動(dòng)區(qū)域及擾動(dòng)程度十分困難。而通過(guò)數(shù)值模擬方法來(lái)進(jìn)行鉆井過(guò)程中鉆頭的破巖仿真分析較常規(guī)實(shí)驗(yàn)方法容易實(shí)現(xiàn)，且能清楚地了解鉆頭擾動(dòng)區(qū)域和擾動(dòng)程度。借助于有限元軟件ANSYS/LS-DYNA，黃志強(qiáng)等[9-10]進(jìn)行了沖旋鉆頭及潛孔鉆頭的破巖機(jī)理仿真研究；祝效華等[11]進(jìn)行了空氣沖旋鉆井破巖的數(shù)值模擬研究；卜長(zhǎng)根等[12]進(jìn)行了潛孔錘鑿巖瞬態(tài)沖擊過(guò)程的仿真研究；譚青等[13]進(jìn)行了TBM 刀具三維破巖仿真研究。由此可知，借助于有限元軟件ANSYS/LS-DYNA 來(lái)進(jìn)行鉆頭破巖的仿真研究是可行的。本文作者借助有限元軟件ANSYS/LS-DYNA 來(lái)模擬PCS 取樣器鉆進(jìn)天然氣水合物地層過(guò)程，給出了鉆進(jìn)過(guò)程中鉆頭附近井底及井周地層不同時(shí)刻應(yīng)力場(chǎng)及應(yīng)變場(chǎng)從而得出其擾動(dòng)程度和擾動(dòng)區(qū)域以及鉆壓和轉(zhuǎn)速對(duì)擾動(dòng)程度的影響。由于實(shí)際鉆進(jìn)過(guò)程中井下條件復(fù)雜，井底PCS 鉆頭受鉆井液、巖屑等多方面因素影響，而數(shù)值模擬難以完全模擬實(shí)際情況，因此，本文在進(jìn)行數(shù)值模擬時(shí)進(jìn)行了以下簡(jiǎn)化：不考慮井底鉆井液對(duì)鉆具和井壁的影響；不考慮井底巖屑對(duì)鉆具本身的摩擦阻力影響。雖然簡(jiǎn)化后的數(shù)值模擬結(jié)果與真實(shí)擾動(dòng)情況會(huì)有所差別，但是所得結(jié)論對(duì)于實(shí)際工程仍具有一定的參考意義。

1 模型建立

1.1 PCS 鉆頭模型

PCS(圖1)是一種自由下落式展開、液壓驅(qū)動(dòng)、鋼絲繩提取的取心工具，它既采用了目前油田壓力取心技術(shù)，又采用了DSDP 中的取心技術(shù)。PCS 鉆頭由三牙輪鉆頭和超前鉆頭[14](cutting shoe)兩部分組成。

超前鉆頭有多種形式，圖2 所示為201 航次中使用的3 種超前鉆頭。本文建立的PCS 鉆頭模型中的超前鉆頭為圖2 中左下所示結(jié)構(gòu)。由于PCS 鉆頭模型較復(fù)雜，為了簡(jiǎn)化模型，本文建立的模型將不考慮牙輪的自轉(zhuǎn)而只考慮其公轉(zhuǎn)。由于牙輪鉆頭牙齒附近容易產(chǎn)生應(yīng)力集中，而本文研究的重點(diǎn)是井底及井周地層的擾動(dòng)應(yīng)力大小，因此，將牙輪表面簡(jiǎn)化為光滑結(jié)構(gòu)，同時(shí)不考慮鉆頭各部件之間的相對(duì)運(yùn)動(dòng)將鉆頭視為一個(gè)整體。鉆頭模型采用剛體材料，具體材料參數(shù)如表1 所示。

圖1 ODP-PCS 的基本結(jié)構(gòu)Fig.1 Schematic diagram of ODP-PCS

圖2 ODP202 航次中使用的超前鉆頭Fig.2 Cutting shoes used on ODP Leg 201

表1 模型物理力學(xué)性能參數(shù)Table 1 Physico-mechanical properties of models

1.2 水合物地層模型

本文天然氣水合物地層的本構(gòu)模型選取Drucker-Prager(DP)彈塑性模型，服從Mohr-Coulomb模型的近似修正Drucker-Prager 屈服準(zhǔn)則，模型選取的水合物其具體材料參數(shù)如表1 所示。為了達(dá)到減小計(jì)算量而又能同時(shí)獲得PCS鉆頭完全鉆入對(duì)地層的擾動(dòng)程度，建立的地層模型高度為0.6 m。為了保證計(jì)算精度，將PCS 鉆頭與水合物地層接觸處即假定擾動(dòng)嚴(yán)重區(qū)域進(jìn)行網(wǎng)格加密，假定中心擾動(dòng)嚴(yán)重區(qū)半徑為0.2 m，整個(gè)地層模型的半徑為0.4 m。圖3 所示為PCS鉆頭鉆進(jìn)水合物地層的整體模型，共劃分網(wǎng)格單元729 928 個(gè)，其中PCS 鉆頭模型89 928 個(gè)，水合物地層640 000 個(gè)。

圖3 PCS 鉆頭鉆進(jìn)水合物地層的有限元模型Fig.3 Finite element model of PCS bit drilling in gas hydrate formation

1.3 邊界條件及加載方式

添加PCS 鉆頭與水合物地層之間的接觸為面-面侵蝕接觸算法。水合物地層上表面為自由界面，下底面限制豎直方向的位移，水合物地層模型外圍限制所有方向的位移。井底鉆頭受力情況復(fù)雜，本文只考慮PCS 鉆頭在承受一定的鉆壓并以一定的轉(zhuǎn)速進(jìn)行鉆進(jìn)時(shí)其應(yīng)力應(yīng)變場(chǎng)的分布。根據(jù)相關(guān)文獻(xiàn)及ODP 航次報(bào)告所給出的鉆進(jìn)參數(shù)[16]，進(jìn)而選定3 個(gè)不同的鉆壓(20，30 和40 kN)以及3 個(gè)不同的轉(zhuǎn)速(80，100 和120 r/min)。鉆壓施加于整個(gè)鉆頭模型，方向?yàn)檠劂@進(jìn)方向Z 軸向下，施加類型為恒定載荷；初始轉(zhuǎn)速同樣施加于整個(gè)鉆頭模型，方向?yàn)轫槙r(shí)針方向，類型為勻速轉(zhuǎn)動(dòng)。

2 模擬結(jié)果及分析

圖4 PCS 鉆頭鉆進(jìn)水合物地層的模擬過(guò)程Fig.4 Numerical simulation of PCS bit drilling in gas hydrate formation

采用 ANSYS/LS-DYN A 軟件進(jìn)行計(jì)算，LS-PrePost 軟件對(duì)計(jì)算結(jié)果進(jìn)行后處理分析。圖4 所示為鉆壓為20 kN、轉(zhuǎn)速為80 r/min 時(shí)，PCS 鉆頭鉆進(jìn)水合物地層的模擬過(guò)程，分別給出了時(shí)間為0.020，0.065，0.110 和0.190 s 時(shí)水合物地層的等效應(yīng)力云圖。

2.1 鉆壓對(duì)地層最大等效應(yīng)力的影響

對(duì)轉(zhuǎn)速為80 r/min，鉆壓分別為20，30 和40 kN時(shí)的鉆進(jìn)過(guò)程進(jìn)行模擬，得到當(dāng)鉆進(jìn)深度為0.41 m時(shí)，井底地層(取距超前鉆頭0.04 m處)的最大等效應(yīng)力(最大擾動(dòng)程度)與鉆壓的對(duì)應(yīng)關(guān)系曲線，如圖5 所示。由圖5 可知：當(dāng)轉(zhuǎn)速不變，鉆壓增大，井底地層的最大等效應(yīng)力呈增大趨勢(shì)，其中鉆壓由20 kN 增加到30 kN時(shí)，井底地層的最大等效應(yīng)力增加了1.7 MPa；由30 kN 增加到40 kN 時(shí)，井底地層的最大等效應(yīng)力增加了0.8 MPa，增大趨勢(shì)減慢。

當(dāng)鉆壓為20 kN，轉(zhuǎn)速為80 r/min，鉆進(jìn)深度為0.41 m 時(shí)，于鉆頭附近井底地層別選取5 個(gè)距超前鉆頭不同距離的單元節(jié)點(diǎn)，得到了這5 個(gè)單元節(jié)點(diǎn)所在區(qū)域的最大等效應(yīng)力與距超前鉆頭距離的關(guān)系曲線，如圖6 所示。由圖6 可知：距離超前鉆頭越近，其最大等效應(yīng)力越大，擾動(dòng)程度越大，距離越遠(yuǎn)，擾動(dòng)程度越小。其中距離超前鉆頭0.09 m 以內(nèi)的地層其最大等效應(yīng)力隨距離的增大而急劇減小，而距離超前鉆頭0.09 m 以外的地層其最大等效應(yīng)力隨距離的增大而緩慢減小。

圖5 超前鉆頭附近井底地層最大等效應(yīng)力與鉆壓關(guān)系曲線Fig.5 Maximum effective stress of undrilled sub bottom hole neared cutting shoes versus bit weigh

當(dāng)轉(zhuǎn)速為80 r/min 時(shí)，得到了鉆進(jìn)過(guò)程中鉆頭附近井周地層(取距鉆頭0.02 m處)的最大等效應(yīng)力(最大擾動(dòng)程度)與鉆壓的對(duì)應(yīng)關(guān)系曲線，如圖7 所示。由圖7 可知：當(dāng)轉(zhuǎn)速不變、鉆壓增大時(shí)，井周地層的最大等效應(yīng)力整體呈線性增大趨勢(shì)，其中鉆壓由20 kN增加到40 kN 時(shí)，井周地層的最大等效應(yīng)力增加了3 MPa。

圖6 鉆頭附近井底地層最大等效應(yīng)力與距超前鉆頭距離關(guān)系曲線Fig.6 Maximum effective stress of undrilled sub bottom hole versus distance to cutting shoe

圖7 鉆頭附近井周地層最大等效應(yīng)力與鉆壓關(guān)系曲線Fig.7 Maximum effective stress of underground layer around well versus bit weight

2.2 轉(zhuǎn)速對(duì)地層最大等效應(yīng)力的影響

對(duì)鉆壓為20 kN，轉(zhuǎn)速分別為80，100 和120 r/min時(shí)的鉆進(jìn)過(guò)程進(jìn)行了模擬，得到了當(dāng)鉆進(jìn)深度為0.41 m 時(shí)，井底地層(取距超前鉆頭0.04 m 處)的最大等效應(yīng)力(最大擾動(dòng)程度)與轉(zhuǎn)速的對(duì)應(yīng)關(guān)系曲線，如圖8所示。由圖8 可知：當(dāng)鉆壓不變，轉(zhuǎn)速增大時(shí)，井底地層的最大等效應(yīng)力基本保持不變，即井底地層的最大等效應(yīng)力不受轉(zhuǎn)速的影響。

當(dāng)鉆壓為20 kN 時(shí)，得到了當(dāng)鉆進(jìn)過(guò)程中，鉆頭附近井周地層(取距PCS 鉆頭0.02 m 處)的最大等效應(yīng)力(最大擾動(dòng)程度)與轉(zhuǎn)速的對(duì)應(yīng)關(guān)系曲線，如圖9 所示。由圖9 可知：當(dāng)鉆壓不變，轉(zhuǎn)速增大，井周地層的最大等效應(yīng)力整體呈線性增大趨勢(shì)，其中轉(zhuǎn)速由80 r/min 增加到120 r/min 時(shí)，井周地層的最大等效應(yīng)力增加了1.9 MPa。

圖8 鉆頭附近井底地層最大等效應(yīng)力與轉(zhuǎn)速關(guān)系曲線Fig.8 Maximum effective stress of undrilled sub bottom hole versus rotary speed

圖9 鉆頭附近井底地層最大等效應(yīng)力與轉(zhuǎn)速關(guān)系曲線Fig.9 Maximum effective stress of underground layer around well versus rotary speed

當(dāng)鉆壓為20 kN，轉(zhuǎn)速為120 r/min 時(shí)，鉆進(jìn)深度為0.41 m 時(shí)，于鉆頭附近井周地層別選取5 個(gè)距鉆頭不同距離的單元節(jié)點(diǎn)，得到了這5 個(gè)單元節(jié)點(diǎn)所在區(qū)域的最大等效應(yīng)力與距鉆頭體距離的關(guān)系曲線，如圖10 所示。由圖10 可知：距離鉆頭體越近，其最大等效應(yīng)力越大，擾動(dòng)程度越大，距離越遠(yuǎn)，擾動(dòng)程度越小。其中距離鉆頭體0.055 m 以內(nèi)的井周地層其最大等效應(yīng)力隨距離的增大而急劇減小，而距離超前鉆頭0.055 m 以外的井周地層其最大等效應(yīng)力隨距離的增大而緩慢減小。

圖10 鉆頭附近井周地層最大等效應(yīng)力與距鉆頭體距離關(guān)系曲線Fig.10 Maximum effective stress of underground layer around well versus distance to bit

3 結(jié)論

(1) 當(dāng)轉(zhuǎn)速一定、鉆壓增大時(shí)，超前鉆頭附近井底地層及鉆頭附近井周地層擾動(dòng)程度呈增大趨勢(shì)；當(dāng)鉆壓一定、轉(zhuǎn)速增大時(shí)，超前鉆頭附近井底地層最大擾動(dòng)程度基本不受轉(zhuǎn)速改變的影響，而鉆頭附近井周地層最大擾動(dòng)程度整體呈線性增大趨勢(shì)。

(2) 鉆壓和轉(zhuǎn)速一定時(shí)，鉆進(jìn)方向上井底地層距離超前鉆頭越近，其水合物地層擾動(dòng)程度越強(qiáng)，反之則越弱，距離超前鉆頭0.09 m 以內(nèi)的井底地層擾動(dòng)程度較大，在此范圍之外的井底地層則擾動(dòng)程度較?。煌瑫r(shí)，徑向上距離鉆頭體越近，其水合物地層擾動(dòng)程度越強(qiáng)，反之則越弱；距離鉆頭體0.055 m 以內(nèi)的井周地層擾動(dòng)程度較大，在此范圍之外的井周地層則擾動(dòng)程度較小。

[1] 周文軍. 天然氣水合物鉆探的井控工藝和參數(shù)計(jì)算[D]. 青島:中國(guó)石油大學(xué)(華東)石油工程學(xué)院, 2008: 1-2.ZHOU Wenjun. Study on well control technology and hydraulic parameters calculation of the gas hydrate drilling[D]. Qingdao:China University of Petroleum (East China). School of Petroleum Engineering, 2008: 1-2.

[2] Winters W J, Pecher I A, Waite W F, et al. Physical properties and rock physics models of sediment containing natural and laboratory-formed methane gas hydrate[J]. American Mineralogist, 2004, 89: 1221-1227.

[3] 沈珠江. 天然軟粘土的壓縮曲線究竟是什么樣子?——答鄒越強(qiáng)的質(zhì)疑[J]. 巖土工程學(xué)報(bào), 1997, 19(3): 118.SHEN Zhujiang. What is the compression curve of natural soft clays[J]. Chinese Journal of Geotechnical Engineering, 1997,19(3): 118.

[4] 馬巍, 吳紫汪, 盛煌. 圍壓對(duì)凍土強(qiáng)度特性的影響[J]. 巖土工程學(xué)報(bào), 1995, 17(5): 7-11.MA Wei, WU Ziwang, SHENG Huang. effect of confining pressure on strength behaviour of frozen soil[J]. Chinese Journal of Geotechnical Engineering, 1995, 17(5): 7-11.

[5] 殷宗澤. 土體沉降與固結(jié)[M]. 北京: 中國(guó)電力出版社, 1998:1-73.YIN Zongze. Soil settlement and consolidation[M]. Beijing:China Power Press, 1998: 1-73.

[6] 楊通佑. 石油及天然氣儲(chǔ)量計(jì)算方法[M]. 2 版. 北京: 石油工業(yè)出版社, 1998: 71-72.YANG Tongyou, Calculation method of oil and natural gas reserves[M]. 2nd ed. Beijing: Petroleum Industry Press, 1998:71-72.

[7] Dickens G R, Paull C K, Wallace P, et al. Direct measurement of in situ methane quantities in a large gas hydrate reservoir[J].Nature, 1997, 385: 426-428.

[8] 白玉湖, 李清平. 天然氣水合物取樣技術(shù)及裝置進(jìn)展[J]. 石油鉆探技術(shù), 2010, 38(6): 116-123.BAI Yuhu, LI Qingping. Progress on natural gas hydrate sampling techniques and tools[J]. Petroleum Drilling Techniques,2010, 38(6):116-123.

[9] 黃志強(qiáng), 范永濤, 魏振強(qiáng), 等. 沖旋鉆頭破巖機(jī)制仿真研究[J].西南石油大學(xué)學(xué)報(bào)(自然科學(xué)版), 2010, 32(1): 148-150.HUANG Zhiqiang, FAN Yongtao, WEI Zhenqiang, et al.Emulation study of rock breaking mechanism of percussion rotary bit[J]. Chinese Journal of Southwest Petroleum University(Science and Technology), 2010, 32(1): 148-150.

[10] 黃志強(qiáng), 宋嘉寧, 李琴, 等. 基于LS-DYNA 的潛孔鉆頭仿真分析與布齒優(yōu)化設(shè)計(jì)[J]. 礦山機(jī)械, 2009, 37(5): 14-17.HUANG Zhiqiang, SONG Jianing, LI Qing, et al. Simulation analysis on the downhole bit based on LS-DYNA and optimization of tooth distribution[J]. Mining & Processing Equipment, 2009, 37(5): 14-17.

[11] 祝效華, 羅衡, 賈彥杰. 考慮巖石疲勞損傷的空氣沖旋鉆井破巖數(shù)值模擬研究[J]. 巖石力學(xué)與工程學(xué)報(bào), 2012, 31(4):754-761.ZHU Xiaohua, LUO Heng, JIA Yanjie. Numerical analysis of air hammer bit drilling based on rock fatigue model[J]. Chinese Journal of Rock Mechanics and Engineering, 2012, 31(4):754-761.

[12] 卜長(zhǎng)根, 龔漢松, 夏柏如. 空隙率對(duì)潛孔錘鑿巖瞬態(tài)沖擊過(guò)程的影響[J]. 探礦工程(巖土鉆掘工程), 2011, 38(4): 6-9.BU Changgen, GONG Hansong, XIA Boru. Simulation study of rock voids rate affecting on impact force of DTH hammer[J].Exploration Engineering (Rock & Soil Drilling and Tunneling),2011, 38(4): 6-9.

[13] 譚青, 張魁, 夏毅敏. TBM 刀具三維破巖仿真[J]. 山東大學(xué)學(xué)報(bào)(工學(xué)版), 2009, 39(6): 72-77.TAN Qing, ZHANG Kui, XIA Yimin, et al. Three-dimensional simulation of rock breaking by TBM cutter[J]. Journal of Shandong University (Engineering Science), 2009, 39(6): 72-77.

[14] 蔣國(guó)盛, 王榮璟, 黎忠文. 天然氣水合物的鉆進(jìn)過(guò)程控制和取樣技術(shù)[J]. 探礦工程(巖土鉆掘工程), 2001(3): 33-35.JIANG Guosheng, WANG Rongjing, LI Zhongwen. Drilling procedure control and sampling of gas hydrates[J]. Exploration Engineering (Rock & Soil Drilling and Tunneling), 2001(3):33-35.

[15] Manoochehr Salehabadi, JIN Min, YANG Jinhai, et al. Finite element modelling of casing in gas hydrate bearing sediments[J].SPE Drilling & Completion, 2008, 24(4): 545-552.

[16] Shipboard Scientific Party Leg 204. Leg 204 preliminary report[EB/OL]. [2009-12-20]. http://www.odp.tamu.edu/publications/prelim/204_prel/204PREL.PDF.