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

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

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

1 模型建立

1.1 PCS 鉆頭模型

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

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

1.2 水合物地層模型

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

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

1.3 邊界條件及加載方式

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

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

圖4 PCS 鉆頭鉆進水合物地層的模擬過程Fig.4 Numerical simulation of PCS bit drilling in gas hydrate formation

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

2.1 鉆壓對地層最大等效應力的影響

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

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

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

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

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

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

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

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

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

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

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

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

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

3 結(jié)論

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

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

[1] 周文軍. 天然氣水合物鉆探的井控工藝和參數(shù)計算[D]. 青島:中國石油大學(華東)石油工程學院, 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] 沈珠江. 天然軟粘土的壓縮曲線究竟是什么樣子?——答鄒越強的質(zhì)疑[J]. 巖土工程學報, 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] 馬巍, 吳紫汪, 盛煌. 圍壓對凍土強度特性的影響[J]. 巖土工程學報, 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]. 北京: 中國電力出版社, 1998:1-73.YIN Zongze. Soil settlement and consolidation[M]. Beijing:China Power Press, 1998: 1-73.

[6] 楊通佑. 石油及天然氣儲量計算方法[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]. 石油鉆探技術(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] 黃志強, 范永濤, 魏振強, 等. 沖旋鉆頭破巖機制仿真研究[J].西南石油大學學報(自然科學版), 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] 黃志強, 宋嘉寧, 李琴, 等. 基于LS-DYNA 的潛孔鉆頭仿真分析與布齒優(yōu)化設(shè)計[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]. 巖石力學與工程學報, 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] 卜長根, 龔漢松, 夏柏如. 空隙率對潛孔錘鑿巖瞬態(tài)沖擊過程的影響[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]. 山東大學學報(工學版), 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] 蔣國盛, 王榮璟, 黎忠文. 天然氣水合物的鉆進過程控制和取樣技術(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.