天然氣水合物地層PCS 取樣擾動有限元仿真分析

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

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

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

圖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)時其應(yīng)力應(yīng)變場的分布。根據(jù)相關(guān)文獻(xiàn)及ODP 航次報告所給出的鉆進(jìn)參數(shù)[16]，進(jìn)而選定3 個不同的鉆壓(20，30 和40 kN)以及3 個不同的轉(zhuǎn)速(80，100 和120 r/min)。鉆壓施加于整個鉆頭模型，方向為沿鉆進(jìn)方向Z 軸向下，施加類型為恒定載荷；初始轉(zhuǎn)速同樣施加于整個鉆頭模型，方向為順時針方向，類型為勻速轉(zhuǎn)動。

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

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

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

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

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

當(dāng)鉆壓為20 kN，轉(zhuǎn)速為80 r/min，鉆進(jìn)深度為0.41 m 時，于鉆頭附近井底地層別選取5 個距超前鉆頭不同距離的單元節(jié)點，得到了這5 個單元節(jié)點所在區(qū)域的最大等效應(yīng)力與距超前鉆頭距離的關(guān)系曲線，如圖6 所示。由圖6 可知：距離超前鉆頭越近，其最大等效應(yīng)力越大，擾動程度越大，距離越遠(yuǎn)，擾動程度越小。其中距離超前鉆頭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 時，得到了鉆進(jìn)過程中鉆頭附近井周地層(取距鉆頭0.02 m處)的最大等效應(yīng)力(最大擾動程度)與鉆壓的對應(yīng)關(guān)系曲線，如圖7 所示。由圖7 可知：當(dāng)轉(zhuǎn)速不變、鉆壓增大時，井周地層的最大等效應(yīng)力整體呈線性增大趨勢，其中鉆壓由20 kN增加到40 kN 時，井周地層的最大等效應(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)速對地層最大等效應(yīng)力的影響

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

當(dāng)鉆壓為20 kN 時，得到了當(dāng)鉆進(jìn)過程中，鉆頭附近井周地層(取距PCS 鉆頭0.02 m 處)的最大等效應(yīng)力(最大擾動程度)與轉(zhuǎn)速的對應(yīng)關(guān)系曲線，如圖9 所示。由圖9 可知：當(dāng)鉆壓不變，轉(zhuǎn)速增大，井周地層的最大等效應(yīng)力整體呈線性增大趨勢，其中轉(zhuǎn)速由80 r/min 增加到120 r/min 時，井周地層的最大等效應(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 時，鉆進(jìn)深度為0.41 m 時，于鉆頭附近井周地層別選取5 個距鉆頭不同距離的單元節(jié)點，得到了這5 個單元節(jié)點所在區(qū)域的最大等效應(yīng)力與距鉆頭體距離的關(guān)系曲線，如圖10 所示。由圖10 可知：距離鉆頭體越近，其最大等效應(yīng)力越大，擾動程度越大，距離越遠(yuǎn)，擾動程度越小。其中距離鉆頭體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)速一定、鉆壓增大時，超前鉆頭附近井底地層及鉆頭附近井周地層擾動程度呈增大趨勢；當(dāng)鉆壓一定、轉(zhuǎn)速增大時，超前鉆頭附近井底地層最大擾動程度基本不受轉(zhuǎn)速改變的影響，而鉆頭附近井周地層最大擾動程度整體呈線性增大趨勢。

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

[1] 周文軍. 天然氣水合物鉆探的井控工藝和參數(shù)計算[D]. 青島:中國石油大學(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] 沈珠江. 天然軟粘土的壓縮曲線究竟是什么樣子?——答鄒越強的質(zhì)疑[J]. 巖土工程學(xué)報, 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]. 巖土工程學(xué)報, 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ì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] 黃志強, 范永濤, 魏振強, 等. 沖旋鉆頭破巖機制仿真研究[J].西南石油大學(xué)學(xué)報(自然科學(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] 黃志強, 宋嘉寧, 李琴, 等. 基于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]. 巖石力學(xué)與工程學(xué)報, 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]. 山東大學(xué)學(xué)報(工學(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] 蔣國盛, 王榮璟, 黎忠文. 天然氣水合物的鉆進(jìn)過程控制和取樣技術(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.