金 顥, 張俊斌, 何玉發
天然氣水合物降壓開采過程中儲層應力規律分析
金 顥, 張俊斌, 何玉發
(中海石油(中國)有限公司深圳分公司深水工程技術中心, 廣東 深圳 518000)
為了研究天然氣水合物降壓開采過程的儲層應力及其穩定性, 運用線性多孔彈性力學和巖石力學知識, 考慮水合物儲層原始應力、孔隙壓力、滲流附加應力及降壓開采水合物過程中水合物飽和度的變化, 建立了降壓開采天然氣水合物儲層的力學模型, 結合墨西哥灣某處水合物藏的基本參數, 對降壓開采水合物儲層應力變化和開采過程的儲層穩定性進行研究。結果表明: 井底壓力是影響水合物儲層應力變化的關鍵因素之一; 滲流附加應力在一定程度上減小了儲層的應力; 水合物分解儲層應力發生變化, 儲層應力在井壁處的波動最大, 井壁處是整個儲層所受軸向偏應力最大的位置, 因此井壁處是優先發生剪切破壞的位置; 為了儲層的穩定性, 降壓開采水合物生產壓差應小于2.19 MPa。
降壓開采; 天然氣水合物; 滲流附加應力; 水合物飽和度; 臨界生產壓差; 儲層穩定性
目前, 天然氣被認為是21世紀的一種未來優質、高效、潔凈的能源, 天然氣的蘊藏量為現有地球化石燃料總碳量的2倍左右[1]。然而, 天然氣能否真正成為未來的可用能源, 關鍵在于能否有效地對天然氣水合物進行合理的開采[2]。
天然氣水合物廣泛分布于海底沉積層和極地凍土帶。通常天然氣水合物以固態形式充填于多孔介質儲層中, 天然氣水合物成藏相比于常規儲層條件要復雜, 并且不具備流動性。因此, 根據天然氣水合物的存在條件和基本性質決定了天然氣水合物的開采方式[3]。天然氣水合物理論上的開采思路是: 首先通過一定的方式打破天然氣水合物儲層的平衡性, 促使水合物在儲層中發生分解, 進而開采出分解的天然氣。目前, 天然氣水合物廣泛的開采方法有: 注化學劑法、加熱法和降壓法等三種[4]。注化學劑法較加熱法作用緩慢, 且二者使用的費用昂貴; 降壓法是通過降低儲層壓力使其低于水合物相平衡壓力引起儲層中水合物分解的方法, 降壓法開采天然氣水合物由于其不需要昂貴的連續激發, 所以被認為是最具有商業前景的開采方法[5]。
天然氣水合物在降壓開采過程中, 由于降低井底壓力, 打破了水合物儲層的穩定性, 促使水合物發生分解。受水合物分解的影響, 分解區水合物儲層膠結性變差、力學強度大幅度降低, 進而近井區儲層呈現出弱膠結、低強度、高孔滲的特點, 類似于常規疏松砂巖[6]。由于水合物降壓分解引起水合物儲層物性參數及巖石力學性質等發生一系列變化, 如水合物飽和度減小、儲層內聚力降低、儲層孔隙度增大、滲透率增大等, 導致水合物儲層應力發生改變, 致使水合物儲層穩定性降低, 達到或超過本身固有的剪切強度, 使近井區水合物儲層發生剪切破壞, 從而影響天然氣水合物安全有效的開采。
由于天然生成的天然氣水合物巖樣取芯困難, 目前, 國內外學者對天然氣水合物的研究大多是在實驗室內完成的。因此, 基于前人的研究, 通過分析水合物儲層的穩定性, 本文在采用降壓法開采天然氣水合物過程中考慮儲層滲流附加應力的影響, 同時利用疊加原理建立了天然氣水合物降壓開采過程中儲層穩定性的力學模型, 進行了降壓開采水合物過程中儲層應力的模擬計算, 分析天然氣水合物在降壓開采過程中儲層應力變化特征及儲層的穩定性, 為有效開采天然氣水合物提供一定的理論基礎。
通常天然氣水合物存在于疏松砂巖沉積層中, 埋藏深度淺。在降壓開采天然氣水合物過程中, 可將整個水合物儲層劃分為兩個區域, 即水合物分解區、儲層穩定區, 如圖1所示[4, 7]。降壓開采天然氣水合物時, 當水合物儲層內壓力低于水合物相平衡時, 打破了水合物儲層的平衡條件, 促使水合物發生分解。天然氣水合物的分解將會導致儲層的穩定性發生改變, 而水合物儲層的應力在一定程度上影響儲層的穩定性。因此, 水合物儲層應力的變化成為影響水合物儲層穩定性的重要因素, 而井底壓力是分析水合物儲層應力變化的關鍵因素。

圖1 水合物儲層區域劃分示意圖
在研究降壓開采天然氣水合物儲層應力變化時, 結合國內外學者建立的有關水合物飽和度、儲層孔隙度、內聚力及儲層應力等相關模型之間的關系, 構建水合物儲層力學模型。由于降壓開采天然氣水合物的過程是非等溫滲流過程, 因此該模型在研究過程中做了如下簡化假設: (1) 天然氣水合物儲層為均質、各向同性的線彈性體。(2) 水合物儲層僅考慮三相三組分。三相: 氣、水、水合物; 三組分: 氣、水、水合物。其中, 水合物為固相, 氣體僅含甲烷, 且不重新生成新的水合物。(3) 降壓開采水合物過程中, 水合物儲層的溫度變化忽略不計, 本文注重考慮井底壓力的變化影響水合物儲層應力改變的問題。(4) 僅考慮氣、水兩相流, 且兩相滲流符合廣義Darcy定律。(5) 井為垂直井。
在降壓開采水合物前, 水合物儲層受上覆巖層壓力、最大水平主應力和最小水平主應力, 水合物儲層處于應力平衡狀態。當水合物儲層被鉆開時, 近井區水合物儲層應力平衡被打破, 儲層應力發生變化, 引起水合物儲層應力的重新分布, 水合物儲層的穩定性發生改變。因此, 井周圍的應力分布可表示為[8-9]:


劉玉石等[10]通過運用有限元方法分析了井周儲層穩態達西滲流時儲層井壁圍巖應力的變化, 從而得出井壁處的徑向應力和周向應力與線彈性理論計算的結果存在很大的差別。水合物儲層為多孔介質儲層, 考慮到水合物儲層的滲透性, 同時又由于儲層流體在向孔隙中的徑向流動過程中會對儲層周圍產生應力, 因此滲流附加應力會對儲層產生一定的影響, 其相關表達式為[11]:


應用線性疊加原理, 由公式(1)和(2)得到天然氣水合物儲層近井區任意一點(≥1)處的應力分布式如下:

(1) 初始條件
天然氣水合物開采前儲層壓力滿足:

(2) 近井區水合物儲層初始分量



(1) 水合物儲層孔隙壓力的分布
開采前, 整個水合物儲層處于穩定的狀態, 儲層中各點壓力均等于原始地層應力。假設天然氣水合物儲層為一圓形區域, 在降壓開采水合物過程中, 圓形封閉儲層中任意一點時刻的壓力滿足彈性不穩定滲流邊界條件, 其控制方程可表示為[12-13]:




整理方程(8), 得到降壓過程中水合物儲層孔隙壓力的分布函數:

當=1時, 得到井底壓力w隨時間變化的關系式為:

式中:為水合物儲層厚度, m;為通過井孔的砂巖流量, cm3/s。
通過式(9)及式(10)可得:

式(11)可化為:

(2) 水合物飽和度隨時間變化
降壓開采水合物過程中水合物分解動力學方程可表示為:

根據降壓法開采中Kim-Bishnoi模型, 天然氣水合物分解速率可表示為:

水興百業旺。結合飲水安全工程建設,各地的鄉村振興建設也搞得有聲有色。有了水,農戶就把改廁、改灶、改路結合起來,很多農村的面貌煥然一新。



分解區水合物儲層內氣、水、水合物之間的關系式:
H+W+g= 1, (16)

其中,H=HH/H, 式中:H為分解水合物的體積;H為水合物物質的量, mol;H為水合物密度, 0.918 g/cm3。
多孔介質孔隙體積為:

其中,s= 2πd。

(3) 水合物儲層孔隙度
隨著水合物儲層埋深、狀態分布等特點, 水合物儲層物性參數發生變化。隨著埋深的增加, 地層沉積壓實作用增大, 水合物儲層孔隙度減小。多孔介質儲層孔隙度與有效應力關系式為[15-16]:

Nazridoust等[17]假設天然氣水合物儲層孔隙度均勻分布且忽略了溫度的影響, 得出儲層孔隙度與水合物飽和度之間的關系式為:

由式(20)和式(21), 可得到孔隙度變化動態模型:

式中:為儲層有效應力, MPa。
(4) 水合物儲層內聚力及有效應力
在降壓開采水合物過程中, 由于水合物的分解導致儲層膠結性變差, 儲層內聚力減小。隨著水合物的分解, 水合物儲層的內聚力降低, 儲層內聚力與孔隙度之間的數學表達式如下[18]:

式中:為水合物分解后儲層內聚力, MPa;0為初始水合物儲層內聚力, MPa;為假設為1.2;為水合物儲層孔隙度增量。
應用Terzaghi有效應力原理, 設儲層巖石應力受壓為負, 受拉為正[19], 其表達式為:

式中:為有效應力, MPa;v為儲層原始垂向應力, MPa;為Biot系數;為Biot有效應力參數。
(5) 天然氣水合物相平衡
天然氣水合物廣泛分布于極地凍土帶和海底沉積層中, 水合物在儲層中存在穩定帶即水合物在儲層中特定溫度和壓力條件下的一個范圍, 在這個范圍內水合物達到相平衡, 甲烷水合物相溫度和壓力平衡采用Makogon[20]模型計算, 得出表達式為:
lge=(–0) +(–0)2+, (25)
式中:e、0為分別是水合物分解的平衡壓力(MPa)和平衡溫度(K);= 0.034 2K–1,= 0.000 5K–2,= 6.480 4。
由于降壓開采水合物導致的變形幾乎是注熱的10倍[21], 因此本文在研究降壓開采水合物過程中, 假設水合物儲層內的溫度變化忽略不計。
(6) 儲層穩定性判據準則
在研究水合物儲層的穩定性時, 運用Drucker- Prager準則進行判斷:

即,



降壓開采水合物過程中生產壓差的計算公式表示為:
Δ=0–w, (28)
式中: Δ為水合物儲層臨界生產壓差, MPa;0為原始水合物儲層的孔隙壓力, MPa;w為井底壓力, MPa。
將式(3)、(29)代入(28)化簡, 最終得到降壓開采水合物過程中儲層穩定性指數與生產壓差之間的關系式為:

模型計算所選用的參數采用墨西哥灣某處甲烷水合物藏的部分基本參數[18]來分析降壓開采水合物過程中近井區儲層應力的變化特征。
墨西哥灣某處水合物藏的部分基本參數: 水合物儲層溫度288 K; 儲層初始壓力16.9 MPa; 井孔半徑0.15 m; 儲層垂向原始主應力21.8 MPa; 儲層最大水平原始主應力20.45 MPa; 儲層最小水平原始主應力19.70 MPa; 儲層厚度30 m; Biot有效應力參數0.6; Biot系數1.0; 水合物飽和度0.5; 儲層孔隙度0.4; 泊松比0.4; 儲層初始滲透率0.2 μm2; 多孔介質壓縮系數10–9Pa–1; 水的壓縮系數4.6×10–10Pa–2; 儲層內聚力1.8 MPa; 內摩擦角30°。其中, 泄油半徑100 m。
圖2、圖3分別給出了生產壓差3 MPa, 不同生產時間下水合物飽和度分布、儲層滲透率分布曲線圖。降壓開采水合物過程中, 初始水合物儲層飽和度受已存自由氣體的影響。

圖2 水合物飽和度分布曲線
從圖2中可以看出, 由降壓開采引起的水合物飽和度變化主要發生在近井區, 在降壓開采水合物初期(20 d), 靠近井壁處的水合物飽和度減小的相對緩慢, 這主要是因為原水合物儲層內存在的自由氣體使得儲層內的壓力升高, 在一定程度上阻礙了水合物的分解速率, 使得水合物的分解速率減慢; 在較遠處, 水合物飽和度趨于穩定。隨著時間不斷增加(60 d), 儲層內的自由氣體逐漸減少, 水合物飽和度減小的速度逐漸加快, 這與唐良廣等[22]研究的水合物飽和度變化規律相符, 進一步證明了該水合物飽和度變化模型的準確性。

圖3 水合物儲層滲透率分布曲線
分析圖3可以看出, 井壁處的儲層滲透率增大的最快, 遠離井壁儲層滲透率逐漸恢復到原始狀態。在水合物分解區, 特別是井壁處水合物分解速率最快, 水合物飽和度減小的最快, 根據式(19)可知, 水合物飽和度和儲層滲透率呈反比。因此, 生產時間越長, 水合物飽和度變得越小, 儲層滲透率就越大。
水合物儲層為多孔介質儲層, 通過模擬計算得到在生產壓差3 MPa, 開井生產100 d后, 研究儲層應力過程中考慮滲流附加應力與不考慮滲流附加應力兩種情況下對水合物儲層徑向應力、周向應力進行了對比, 如圖4所示。

圖4 存在滲流附加應力與無滲流附加應力兩種情況下儲層應力的對比
從圖4中可以得出, 考慮滲流附加應力時的儲層徑向、周向應力小于不考慮滲流附加應力時的儲層徑向、周向應力。通過對式(2)分析可知, 水合物降壓分解過程中, 飽和度逐漸減小, 儲層有效孔隙度逐漸增大, 這說明滲流附加應力降低了儲層的穩定性。
在井底壓力13.9 MPa, 井眼半徑0.15 m的條件下, 分析考慮近井區水合物儲層滲流附加應力和不考慮滲流附加應力時, 儲層所受最大軸向偏應力分布曲線如圖5所示。結合圖4, 由圖5不難得出, 考慮滲流附加應力時儲層所受最大軸向偏應力小于不考慮滲流附加應力時儲層所受最大軸向偏應力。因此, 在建立水合物儲層的力學模型時, 必須考慮滲流附加應力對水合物儲層的影響, 以減小預測結果的誤差。

圖5 儲層最大軸向偏應力分布
通過降低井底壓力, 使水合物儲層壓力低于其平衡壓力, 促使水合物發生分解生成天然氣和水。井底壓力變化會改變一定區域的原始儲層應力, 造成分解區水合物儲層發生應力集中, 從而降低水合物儲層的穩定性。在降壓開采水合物過程中水合物儲層徑向、周向、垂向應力表現出不同的變化。圖6顯示了在生產壓差3 MPa, 開采時間100 d后, 塑性區水合物儲層應力的分布曲線。
通過式(3)不難看出, 水合物儲層應力是角度的函數, 儲層周向、垂向應力隨角度呈周期性變化, 儲層徑向應力不發生變化。近井區水合物儲層中最大、最小應力之間的差值隨著距井孔距離的增大逐漸減小, 在井周上某點與水平最大主應力方向的夾角90°或270°方位處, 水合物儲層的最大、最小應力之間的差值最大, 即近井區水合物儲層所受到的擠壓應力最大。根據Drucker-Prager巖石破壞準則, 儲層周向、徑向應力之間的差值越大, 儲層越容易發生剪切破壞。因此, 這兩個方位處也是儲層最容易發生剪切破壞的位置, 儲層發生剪切破壞一般是從井壁處開始的。

圖6 儲層應力隨角度的變化曲線
Drucker-Prager準則認為儲層發生的破壞主要是剪切破壞, 根據式(29)進一步說明儲層破壞的臨界條件。保持其他參數不變, 水合物儲層穩定性隨生產壓差變化的關系圖如圖7所示。由圖7可知, 水合物儲層穩定性指數隨生產壓差的增大而減小。這是因為井底壓力越小, 即生產壓差越大, 儲層有效應力增大越快, 造成儲層失穩的可能性就越大; 隨著生產壓差的增大, 水合物分解速率加快, 儲層物性變化越顯著。

圖7 儲層穩定性隨生產壓差變化關系圖
通過對圖7曲線進行線性回歸得到儲層穩定性指數與生產壓差之間的關系式(30)所示:
= 0.247 3Δ2– 2.803 7Δ+ 5.092 9, (30)
式中:為儲層穩定性指數, MPa; Δ為生產壓差, MPa。
通過解式(30), 得到水合物儲層臨界生產壓差為2.19 MPa。當生產壓差小于2.19 MPa, 水合物儲層保持穩定; 相反, 當生產壓差大于2.19 MPa, 水合物儲層發生剪切強度破壞, 有可能會引起儲層出砂。因此, 降壓開采水合物過程中, 應保持井底壓力大于14.71 MPa。
綜合以上分析可知, 在降壓開采水合物過程中, 水合物分解區是整個儲層穩定性較差的區域。降低井底壓力: 一方面, 促進了水合物的分解, 儲層內聚力減小, 儲層各應力發生改變; 另一方面, 不合理的降壓會造成儲層發生剪切破壞, 進而影響水合物的有效開采。因此, 在降壓開采水合物過程中應保證合理的井底壓力。
1) 降壓開采水合物過程中, 由于水合物發生分解, 分解區儲層呈現出弱膠結、低強度、高孔滲的特點, 在不同生產壓差下, 近井區儲層應力發生了明顯的改變, 距井孔較遠處儲層應力不再發生變化而趨于原始儲層應力。說明近井區儲層出現應力集中的現象, 井底壓力是影響儲層穩定性的重要因素。
2) 生產壓差越大, 塑性區水合物儲層所受軸向偏應力就越大, 井壁處所受軸向偏應力最大, 是整個水合物儲層優先破壞的位置。其中, 生產壓差2.19 MPa作為降壓開采水合物過程中儲層發生剪切破壞的臨界壓差。
3) 有待進一步開展對天然氣水合物降壓開采過程中儲層應力變化分析, 確定水合物儲層發生剪切破壞的生產壓差的合理臨界點, 為今后有效開采天然氣水合物奠定一定的基礎。
[1] Makogon Y F, Holditch S A, Makogon T Y. Natural gas- hydrates: A potential energy source for the 21st Century[J]. Journal of Petroleum Science & Engineering, 2007, 56(1-3): 14-31.
[2] Moridis G J, Collett T S, Dallimore S R, et al. Numerical studies of gas production from several CH4, hydrate zones at the Mallik site, Mackenzie Delta, Canada[J]. Journal of Petroleum Science & Engineering, 2004, 43(3-4): 219-238.
[3] Sloan E D. Clathrate hydrates of natural gases[J]. Journal of the American Chemical Society, 1998, 120(43): 11212.
[4] Ji C, Ahmadi G, Smith D H. Natural gas production from hydrate decomposition by depressurization[J]. Chemical Engineering Science, 2001, 56(20): 5801-5814.
[5] Wang Yi, Li Xiaosen, Xu Wenyue, et al. Experimental investigation into factors influencing methane hydrate formation and a novel method for hydrate formation in porous media[J]. Energy and Fuel, 2013, 27(7): 3751- 3757.
[6] 喻西崇, 吳應湘, 安維杰, 等. 開采地層中的天然氣水合物的數學模型[J]. 天然氣工業, 2004, 24(1): 63-67. Yu Xichong, Wu Yingxiang, An Weijie, et al. Mathema-tical model to recover gas hydrate from formations[J]. Natural Gas Industry, 2004, 24(1): 63-67.
[7] Ahmadi G, Ji C, Smith D H. Numerical solution for natural gas production from methane hydrate dissociation[J]. Journal of Petroleum Science & Engineering, 2004, 41(4): 269-285.
[8] Aadnoy B S, Chenevert M E. Stability of highly inclined boreholes[J]. SPE Drilling Engineering, 1987, 2: 364-374.
[9] Fjar E. Petroleum Related Rock Mechanics[M]. Amsterdam: Elsevier, 2008.
[10] 劉玉石, 黃克累. 孔隙流體對井眼穩定的影響[J]. 石油鉆探技術, 1995, 23(3): 4-6. Liu Yushi, Huang Kelei. Effects of pore fluids on borehole stability[J]. Petroleum Drilling Techniques, 1995, 23(3): 4-6.
[11] 閻鐵, 李士斌. 深部井眼巖石力學理論與應用[M]. 北京: 石油工業出版社, 2002. Yan Tie, Li Shibin. Deep Borehole Rock Mechanics Theory and Applications[M]. Beijing: Petroleum Industry Press, 2002.
[12] 張建國, 杜殿發, 侯建, 等. 油氣層滲流力學[M]. 東營: 中國石油大學出版社, 2010. Zhang Jianguo, Du Dianfa, Hou Jian, et al. The Oil and Gas Flow through Porous Media[M]. Dongying: China University of Petroleum Press, 2010.
[13] Jacob B. Dynamics of Fluids in Porous Media[M]. New York: American Elsevier Pub. Co., 1972.
[14] Masuda Y, Kurihara M. A filed scale simulation study on gas productivity of formation containing gas hydrates[J]. Proceedings of the Fourth International Conference on Gas Hydrates, 2002, 15(6): 622-627.
[15] 張新紅, 秦積舜. 低滲巖心物性參數與應力關系的試驗研究[J]. 中國石油大學學報: 自然科學版, 2001, 25(4): 56-57. Zhang Xinhong, Qin Jishun. Relation of Ground Stress with the Properties of Low Permeable Core[J]. Journal of China University of Petroleum (Edition of Natural Science), 2001, 25(4): 56-57.
[16] 沈海超. 天然氣水合物藏降壓開采流固耦合數值模擬研究[D]. 青島: 中國石油大學, 2009. Shen Haichao. Fluid-solid coupling numerical simulation on natural gas production from hydrate reservoirs by depressurization[D]. Qingdao: China University of Petroleum, 2009.
[17] Nazridoust K, Ahmadi G. Computational modeling of methane hydrate dissociation in a sandstone core[J]. Che-mical Engineering Science, 2007, 62(22): 6155-6177.
[18] Freij-Ayoub R, Tan C, Clennell B, et al. A wellbore stability model for hydrate bearing sediments[J]. Journal of Petroleum Science & Engineering, 2007, 57(1): 209-220.
[19] 冉啟全, 顧小蕓. 油藏滲流與應力耦合分析[J]. 巖土工程學報, 1998, 20(2): 69-73. Ran Qiquan, Gu Xiaoyun. Coupling analysis of multiphase flow and stress for oil reservoir[J]. Chinese Journal of Geotechnical Engineering, 1998, 20(2): 69-73.
[20] Makogon Y F. Hydrates of hydrocarbons[J]. Penn Well Books, 1997, 610(52): 50-62.
[21] Kimoto S, Oka F, Fushita T, et al. A chemo-thermo- mechanically coupled numerical simulation of the subsurface ground deformations due to methane hydrate dissociation[J]. Computers & Geotechnics, 2007, 34(4): 216-228.
[22] 唐良廣, 馮自平, 李小森, 等. 天然氣水合物降壓開采的數值模擬研究[C]//中國科學技術協會.中國工程熱物理學會多相流學術會議. 北京: 中國科學技術協會, 2006: 691-698. Tang Lianguang, Feng Ziping, Li Xiaoshen, et al. Numerical Simulation on Natural Gas Production from Hydrate Reservoirs by Depressurization[C]// Science and Technology Conference of Chinese. Multiphase Flow Academic Conference of Chinese Society of Engineering Thermophysics. Beijing: Science and Technology Conference of Chinese, 2006: 691-698.
Geomechanical stress law analysis of hydrate reservoir for gas hydrate production by depressurization
JIN Hao, ZHANG Jun-bin, HE Yu-fa
(Deepwater Engineering Technology Center of CNOOC (China) Co., Ltd., Shenzhen 518000, China)
A model combining linear porous elasticity and rock mechanics knowledge is developed to analyze the geomechanical stress of the hydrate reservoir and assess the wellbore stability during depressurization. The proposed model considers the original reservoir stress, pore pressure, seepage additional stress, and depressurization of hydrate mining during the process of hydrate saturation change. This model is applied to the case of the Gulf of Mexico, where the hydrate reservoir basic parameters have been publicly published and depressurization of the hydrate mining reservoir has been simulated and analyzed. Results show that the effect of bottom hole pressure is the dominant factor affecting the formation of reservoir stress. Moreover, reservoir stress significantly affects the wall. Therefore, the wall is the location that is most prone to shear failure. To assess the wellbore stability during the process of depressurization, the production pressure difference of the hydrate reservoir must be <2.19 MPa.
depressurization; gas hydrate reservoir; seepage additional stress; hydrate saturation; critical production pressure differential; wellbore stability
Oct. 31, 2019
TE319
A
1000-3096(2020)11-0001-09
10.11759/hykx20191031002
2019-10-31;
2019-12-10
國家重點研發計劃項目 (2016YFC0304000);國家重點研發計劃課題 (2016YFC03044007);國家重點研發計劃課題 (2018YFC0310205)
[The National Key R&D Program of China, No. 2016YFC0304000; The National Key R&D Program of China, No. 2016YFC03044007; The National Key R&D Program of China, No. 2018YFC0310205]
金顥(1984-),男,廣東深圳人,高級工程師,主要從事鉆井、完井、測試、水合物等領域相關內容研究,電話:15899755760,E-mail:jinhao@cnooc.com.cn
(本文編輯: 劉珊珊)