降压开采海域天然气水合物电潜泵排采生产优化

平晓琳 韩国庆 岑学齐 朱玮涛 白振强 王清华 彭龙

平晓琳,韩国庆,岑学齐,朱玮涛,白振强,王清华,彭龙. 降压开采海域天然气水合物电潜泵排采生产优化[J]. 石油钻采工艺,2022,44(2):225-232 doi:  10.13639/j.odpt.2022.02.014
引用本文: 平晓琳,韩国庆,岑学齐,朱玮涛,白振强,王清华,彭龙. 降压开采海域天然气水合物电潜泵排采生产优化[J]. 石油钻采工艺,2022,44(2):225-232 doi:  10.13639/j.odpt.2022.02.014
PING Xiaolin, HAN Guoqing, CEN Xueqi, ZHU Weitao, BAI Zhenqiang, WANG Qinghua, PENG Long. Drainage and production optimization of electric submersible pump for depressurization recovery of offshore natural gas hydrates[J]. Oil Drilling & Production Technology, 2022, 44(2): 225-232 doi:  10.13639/j.odpt.2022.02.014
Citation: PING Xiaolin, HAN Guoqing, CEN Xueqi, ZHU Weitao, BAI Zhenqiang, WANG Qinghua, PENG Long. Drainage and production optimization of electric submersible pump for depressurization recovery of offshore natural gas hydrates[J]. Oil Drilling & Production Technology, 2022, 44(2): 225-232 doi:  10.13639/j.odpt.2022.02.014

降压开采海域天然气水合物电潜泵排采生产优化

doi: 10.13639/j.odpt.2022.02.014
基金项目: 中石化优青项目“天然气水合物开采产能预测方法与有效开发模式探究”(编号:P20025);南方海洋科学与工程广东省实验室(广州)重大基金项目“琼东南盆地天然气水合物系统成藏机制研究”(编号:GML2019ZD0102)
详细信息
    作者简介:

    平晓琳 (1987-),在读博士研究生,从事天然气水合物相关方向的研究工作。通讯地址:(1002249)北京市昌平区府学路18号。E-mail: ogt_pingxl@vip.sina.com

    通讯作者:

    韩国庆 (1969-),从事采油理论与技术方面的研究工作,教授,博士生导师。通讯地址:(102249) 北京市昌平区府学路 18 号。电话:010-89734339。E-mail: hanguoqing@163.com

  • 中图分类号: TE53

Drainage and production optimization of electric submersible pump for depressurization recovery of offshore natural gas hydrates

  • 摘要: 天然气水合物试采实践表明,电潜泵是降压开采天然气水合物最适宜的人工举升工艺。综合考虑复杂的分采管柱、周围海水环境的温度场、电潜泵和气液分离器等多种因素的影响,建立了降压开采海域天然气水合物电潜泵排采的井筒气液两相流模型,分析了不同管线的传热过程,预测不同管线的流型分布、温度和压力剖面等,并利用该模型进行生产优化。研究结果表明,在天然气水合物降压开采过程中,通过增加井口回压或电潜泵频率能够降低采气管线中的动液面高度,从而降低采气管线连续出水的风险。该研究为降压开采海域天然气水合物电潜泵排水采气的生产优化提供参考依据。
  • 图  1  天然气水合物试采系统示意图

    Figure  1.  Schematic diagram of the pilot production system of natural gas hydrates

    图  2  降压开采天然气水合物井筒气液两相流模型求解流程

    Figure  2.  Workflow to solve the gas-liquid two-phase well flow model for depressurization recovery of natural gas hydrates

    图  3  不同管线的压力、温度和气体体积分数与井深的关系曲线

    Figure  3.  Pressure, temperature, and gas volume fractions vs. well depth for different pipes

    图  4  不同井口压力下的生产管线的压力

    Figure  4.  Pressure of production pipe at different wellhead pressures

    图  5  不同井口压力下生产管线的气体体积分数

    Figure  5.  Gas volume fraction of production pipe at different wellhead pressures

    图  6  不同电潜泵频率下生产管线的压力

    Figure  6.  Pressure of production pipe at different ESP frequencies

    图  7  不同电潜泵频率下生产管线的气体体积分数

    Figure  7.  Gas volume fraction of production pipe at different ESP frequencies

    图  8  CENTRILIFT GC1600电潜泵频率与泵效关系

    Figure  8.  Frequency vs. frequency of the CENTRILIFT GC1600 ESP

    表  1  全球天然气水合物试采情况

    Table  1.   Statistics of global pilot production of natural gas hydrates

    时间地点试采方法举升工艺
    2002年加拿大麦肯齐三角州注热法电潜泵
    2007年加拿大麦肯齐三角州降压法电潜泵
    2008年加拿大麦肯齐三角州降压法电潜泵
    2012年美国阿拉斯加北坡CO2置换法射流泵
    2013年日本南海海槽降压法电潜泵
    2017年日本南海海槽降压法电潜泵
    2017年中国南海神狐海域降压法电潜泵
    2017年中国南海荔湾海域固态流化法连续油管
    2020年中国南海神狐海域降压法电潜泵
    下载: 导出CSV

    表  2  天然气水合物试采井排水采气适应性分析

    Table  2.   Applicability analysis of dewatering gas recovery in pilot production wells of natural gas hydrates

    项目电潜泵电潜螺杆泵射流泵气举
    海上适应性良好良好良好良好
    最大下入
    深度/m
    27002 00028003000
    最大排液量/
    (m3·d−1)
    500100300400
    防砂能力较弱较强较弱
    应用情况日本、中国、加拿大美国
    下载: 导出CSV

    表  3  模型基础参数

    Table  3.   Basic model parameters

    参数数值参数数值
    井深/m1 350平均日产气量/(m3·d−1)20 000
    海水深度/m1 000平均日产水量/(m3·d−1)200
    电潜泵和EIH深度/m1 200气液分离系数0.95
    地温梯度/(℃·(100 m)−1)2电加热器功率/kW30
    井底流压/MPa5电潜泵频率/Hz60
    下载: 导出CSV

    表  4  模拟结果与现场监测数据对比

    Table  4.   Comparison between simulation result and field monitoring

    位置压力/MPa温度/℃
    模拟
    结果
    监测
    数据
    误差模拟
    结果
    监测
    数据
    误差
    泵吸入口4.14.70.618.318.50.2
    泵排出口14.114.80.710.311.51.2
    海底测试树
    (采气管线)
    3.33.80.513.515.21.7
    海底测试树
    (排水管线)
    12.212.50.315.3××
    采气管线井口
    (钻井平台)
    0.60.20.415.215.10.1
    排水管线井口
    (钻井平台)
    2.62.50.115.5××
      注:“×”表示无数据
    下载: 导出CSV
  • [1] 魏伟, 张金华, 于荣泽, 等. 2017年天然气水合物研发热点回眸[J]. 科技导报, 2018, 36(1):83-90. doi:  10.3981/j.issn.1000-7857.2018.01.009

    WEI Wei, ZHANG Jinhua, YU Rongze, et al. Review on natural gas hydrate in 2017[J]. Technology Review, 2018, 36(1): 83-90. doi:  10.3981/j.issn.1000-7857.2018.01.009
    [2] 李守定, 孙一鸣, 陈卫昌, 等. 天然气水合物开采方法及海域试采分析[J]. 工程地质学报, 2019, 27(1):55-68. doi:  10.13544/j.cnki.jeg.2019-065

    LI Shouding, SUN Yiming, CHEN Weichang, et al. Analyses of gas production methods and offshore production tests of natural gas hydrates[J]. Journal of Engineering Geology, 2019, 27(1): 55-68. doi:  10.13544/j.cnki.jeg.2019-065
    [3] DEMIRBAS A. Methane hydrates as potential energy resource: Part 2–Methane production processes from gas hydrates[J]. Energy Conversion and Management, 2010, 51(7): 1562-1571. doi:  10.1016/j.enconman.2010.02.014
    [4] YAMAMOTO K, WANG X X, TAMAKI M, et al. The second offshore production of methane hydrate in the Nankai Trough and gas production behavior from a heterogeneous methane hydrate reservoir[J]. RSC Advances, 2019, 9(45): 25987-26013. doi:  10.1039/c9ra00755e
    [5] WEI N, ZHAO J, SUN W, et al. Non-equilibrium multiphase wellbore flow characteristics in solid fluidization exploitation of marine gas hydrate reservoirs[J]. Natural Gas Industry B, 2019, 6(3): 282-292. doi:  10.1016/j.ngib.2018.10.008
    [6] SAKURAI S, NISHIOKA I, MATSUZAWA M, et al. Issues and challenges with controlling large drawdown in the first offshore Methane-Hydrate production test[J]. SPE Production & Operations, 2017, 32(4): 500-516. doi:  10.2118/182230-PA
    [7] LI J F, YE J L, QIN X W, et al. The first offshore natural gas hydrate production test in South China Sea[J]. China Geology, 2018, 1(1): 5-16. doi:  10.31035/cg2018003
    [8] 程心平, 沈琼, 马喜超, 等. 国内外天然气水合物试采概况[J]. 广东化工, 2020, 47(1):95-96.

    CHENG Xinping, SHEN Qiong, MA Xichao, et al. Overview of natural gas hydrate test mining at domestic and overseas[J]. Guangdong Chemical Industry, 2020, 47(1): 95-96.
    [9] 孙宝江, 王志远, 公培斌, 等. 深水井控的七组分多相流动模型[J]. 石油学报, 2011, 32(6):1042-1049. doi:  10.7623/syxb201106018

    SUN Baojiang, WANG Zhiyuan, GONG Peibin, et al. Application of a seven-component multiphase flow model to deepwater well control[J]. Acta Petrolei Sinica, 2011, 32(6): 1042-1049. doi:  10.7623/syxb201106018
    [10] WEI N, SUN W T, MENG Y F et al. Change mechanism of transient gas-liquid two-phase flow in wellbore during Marine natural gas hydrate reservoir drilling[J]. Thermal Science, 2019, 23(4): 2179-2187.
    [11] ZHANG J, WANG Z, TONG S, et al. Hydrate plugging prevention in deep water gas wells[C]// Paper presented at the SPE/AAPG/SEG Unconventional Resources Technology Conference, July 2020, Virtual: URTEC-2020-3938-MS.
    [12] DENG Z, WANG Z, ZHAO Y, et al. Flow assurance during gas hydrate production: Hydrate regeneration behavior and blockage risk analysis in wellbore[C]//Paper presented at the Abu Dhabi International Petroleum Exhibition & Conference, November 2017, Abu Dhabi, UAE: SPE-188684-MS.
    [13] LIAO Y, SUN X, SUN B, et al. Transient gas-liquid-solid flow model with heat and mass transfer for hydrate reservoir drilling[J]. International Journal of Heat and Mass Transfer, 2019, 141: 476-486. doi:  10.1016/j.ijheatmasstransfer.2019.06.097
    [14] HIROBAYASHI S, KIYONO F, MORITA H, et al. A study on gas-liquid two phase flow in methane hydrate production system[C]// Paper presented at the The 29th International Ocean and Polar Engineering Conference, June 2019, Honolulu, Hawaii, USA: ISOPE-I-19-504.
    [15] 曹光强, 姜晓华, 李楠, 等. 产水气田排水采气技术的国内外研究现状及发展方向[J]. 石油钻采工艺, 2019, 41(5):614-623. doi:  10.13639/j.odpt.2019.05.011

    CAO Guangqiang, JIANG Xiaohua, LI Nan, et al. Domestic and foreign research status and development direction of drainage gas recovery technologies in water-producing gas fields[J]. Oil Drilling & Production Technology, 2019, 41(5): 614-623. doi:  10.13639/j.odpt.2019.05.011
    [16] 薛蓥, 何秀萍, 丁英楠, 等. 东海高气液比电泵排采技术先导试验与推广应用[J]. 中国石油和化工标准与质量, 2020, 40(4):207-211,213. doi:  10.3969/j.issn.1673-4076.2020.04.102

    XUE Ying, HE Xiuping, DING Yingnan, et al. Pilot test and application of high gas-liquid ratio for electric pump drainage technology in East China Sea[J]. China Petroleum and Chemical Standard and Quality, 2020, 40(4): 207-211,213. doi:  10.3969/j.issn.1673-4076.2020.04.102
    [17] DU Q, LIU Y, HOU J, et al. Investigation on the secondary Generation of natural gas hydrates in horizontal wellbore caused by pressure jump during the depressurization development of hydrate bearing layers[J]. Geofluids, 2020, 2020: 8815455. doi:  10.1155/2020/8815455
  • [1] 叶正荣, 裘智超, 周祥, 伊然, 刘翔, 赵志宏.  克拉2气田排水采气工艺优选及效果分析 . 石油钻采工艺, 2022, 44(3): 1-6.
    [2] 宋来明, 王春秋, 卢川, 丁祖鹏, 李桂亮, 檀朝东, 程时清.  数据驱动的复杂油藏注采生产优化技术研究进展 . 石油钻采工艺, 2022, 44(2): 253-260. doi: 10.13639/j.odpt.2022.02.018
    [3] 赵克斌, 孙长青, 吴传芝.  天然气水合物开发技术研究进展 . 石油钻采工艺, 2021, 43(1): 7-14. doi: 10.13639/j.odpt.2021.01.002
    [4] 莫军, 刘勇, 张雅楠, 祝洪爽, 李濛.  葡北油田天然气重力混相驱注气井天然气水合物的形成及防治 . 石油钻采工艺, 2021, 43(5): 668-674. doi: 10.13639/j.odpt.2021.05.017
    [5] 郑成胜, 蓝强, 徐运波, 赵怀珍.  天然气水合物抑制剂YHHI-1的合成及评价 . 石油钻采工艺, 2020, 42(6): 708-713. doi: 10.13639/j.odpt.2020.06.007
    [6] 符东宇, 王雨生, 姜超, 严小勇, 鲁光亮, 黄文杰, 魏朝勇.  超音速雾化排水采气工艺数值模拟与现场应用 . 石油钻采工艺, 2020, 42(3): 322-328. doi: 10.13639/j.odpt.2020.03.012
    [7] 李牧.  下倾型页岩气水平井连续油管排水采气工艺 . 石油钻采工艺, 2020, 42(3): 329-333. doi: 10.13639/j.odpt.2020.03.013
    [8] 陈花, 关富佳, 肖启圣, 程亮.  基于温度震荡的天然气水合物增量合成实验 . 石油钻采工艺, 2019, 41(3): 399-404. doi: 10.13639/j.odpt.2019.03.022
    [9] 曹光强, 姜晓华, 李楠, 贾敏, 张义, 王浩宇.  产水气田排水采气技术的国内外研究现状及发展方向 . 石油钻采工艺, 2019, 41(5): 614-623. doi: 10.13639/j.odpt.2019.05.011
    [10] 王星, 孙子刚, 张自印, LawrenceLau, 粟玲, 王彬.  海域天然气水合物试采实践与技术分析 . 石油钻采工艺, 2017, 39(6): 744-750. doi: 10.13639/j.odpt.2017.06.015
    [11] 李锐, 蔡昌新, 李勇, 廖锐全.  多模式优化下的柱塞气举排水采气控制系统设计 . 石油钻采工艺, 2016, 38(5): 673-677. doi: 10.13639/j.odpt.2016.05.024
    [12] 李旭日, 田伟, 李耀德, 杨亚聪, 李丽, 惠艳妮.  柱塞气举排水采气远程控制系统 . 石油钻采工艺, 2015, 37(3): 76-79. doi: 10.13639/j.odpt.2015.03.017
    [13] 周忠城, 曹和平, 姜楚豪, 白晓弘, 李旭日.  小直径管悬挂+ 智能注剂泡沫排水采气技术 . 石油钻采工艺, 2014, 36(5): 103-105. doi: 10.13639/j.odpt.2014.05.025
    [14] 易劲, 顾友义, 潘云兵, 李秋.  油管穿孔气举排水采气技术 . 石油钻采工艺, 2014, 36(3): 103-105. doi: 10.13639/j.odpt.2014.03.0026
    [15] 李隽, 李楠, 李佳宜, 王云, 曹光强.  涡流排水采气技术数值模拟研究   . 石油钻采工艺, 2013, 35(6): 60-64.
    [16] 甘宝安, 聂暘, 路春明, 王鹏, 宁晓颖, 方红.  水力泵独立排水采气工艺 . 石油钻采工艺, 2012, 34(4): 75-76.
    [17] 张亮, 马认琦, 苏杰, 张向华, 苗典远.  天然气水合物形成机理及有效清除 . 石油钻采工艺, 2010, 32(3): 33-36.
    [18] 李淑霞, 姜兴兴, 姜汉桥, 李清平.  天然气水合物藏注热开采敏感参数分析 . 石油钻采工艺, 2010, 32(2): 54-57.
    [19] 胡海良, 唐海雄, 罗俊丰, 韦红术.  深水天然气水合物钻井及取心技术 . 石油钻采工艺, 2009, 31(1): 27-30.
    [20] 高贵民, 朱世民, 刘东, 崔金榜, 李向齐, 马凤芹, 孙卫华.  苏4-14井借气气举排水采气技术 . 石油钻采工艺, 2006, 28(4): 44-45,49. doi: 10.3969/j.issn.1000-7393.2006.04.014
  • 加载中
图(8) / 表 (4)
计量
  • 文章访问数:  0
  • HTML全文浏览量:  4
  • PDF下载量:  1
  • 被引次数: 0
出版历程
  • 修回日期:  2022-01-29
  • 网络出版日期:  2022-09-19
  • 刊出日期:  2022-09-19

降压开采海域天然气水合物电潜泵排采生产优化

doi: 10.13639/j.odpt.2022.02.014
    基金项目:  中石化优青项目“天然气水合物开采产能预测方法与有效开发模式探究”(编号:P20025);南方海洋科学与工程广东省实验室(广州)重大基金项目“琼东南盆地天然气水合物系统成藏机制研究”(编号:GML2019ZD0102)
    作者简介:

    平晓琳 (1987-),在读博士研究生,从事天然气水合物相关方向的研究工作。通讯地址:(1002249)北京市昌平区府学路18号。E-mail: ogt_pingxl@vip.sina.com

    通讯作者: 韩国庆 (1969-),从事采油理论与技术方面的研究工作,教授,博士生导师。通讯地址:(102249) 北京市昌平区府学路 18 号。电话:010-89734339。E-mail: hanguoqing@163.com
  • 中图分类号: TE53

摘要: 天然气水合物试采实践表明,电潜泵是降压开采天然气水合物最适宜的人工举升工艺。综合考虑复杂的分采管柱、周围海水环境的温度场、电潜泵和气液分离器等多种因素的影响,建立了降压开采海域天然气水合物电潜泵排采的井筒气液两相流模型,分析了不同管线的传热过程,预测不同管线的流型分布、温度和压力剖面等,并利用该模型进行生产优化。研究结果表明,在天然气水合物降压开采过程中,通过增加井口回压或电潜泵频率能够降低采气管线中的动液面高度,从而降低采气管线连续出水的风险。该研究为降压开采海域天然气水合物电潜泵排水采气的生产优化提供参考依据。

English Abstract

平晓琳,韩国庆,岑学齐,朱玮涛,白振强,王清华,彭龙. 降压开采海域天然气水合物电潜泵排采生产优化[J]. 石油钻采工艺,2022,44(2):225-232 doi:  10.13639/j.odpt.2022.02.014
引用本文: 平晓琳,韩国庆,岑学齐,朱玮涛,白振强,王清华,彭龙. 降压开采海域天然气水合物电潜泵排采生产优化[J]. 石油钻采工艺,2022,44(2):225-232 doi:  10.13639/j.odpt.2022.02.014
PING Xiaolin, HAN Guoqing, CEN Xueqi, ZHU Weitao, BAI Zhenqiang, WANG Qinghua, PENG Long. Drainage and production optimization of electric submersible pump for depressurization recovery of offshore natural gas hydrates[J]. Oil Drilling & Production Technology, 2022, 44(2): 225-232 doi:  10.13639/j.odpt.2022.02.014
Citation: PING Xiaolin, HAN Guoqing, CEN Xueqi, ZHU Weitao, BAI Zhenqiang, WANG Qinghua, PENG Long. Drainage and production optimization of electric submersible pump for depressurization recovery of offshore natural gas hydrates[J]. Oil Drilling & Production Technology, 2022, 44(2): 225-232 doi:  10.13639/j.odpt.2022.02.014
  • 天然气水合物又称可燃冰,资源潜力巨大。世界上天然气水合物的总资源量相当于全球已知煤、石油和天然气的2倍,99%以上的水合物资源赋存于海底,只有1%存在于陆地冻土层。1 m3天然气水合物可转化为164 m3的天然气和0.8 m3的水,具有非常高的使用价值[1-3]。近年来,国外多个国家开展了天然气水合物试采研究,我国南海神狐海域2次天然气水合物试采都取得了重大进展[2-8]

    目前,天然气水合物的开采还处于试采阶段,尚未实现大规模商业化开采。天然气水合物当前的研究方向主要集中于天然气水合物富集地的形成机理、勘探与识别方法、天然气水合物的开采方法与技术等。降压开采是一种经济有效的天然气水合物开采方式,而排水采气工艺是实现天然气水合物降压开采的重要手段。井筒多相流模拟是进行降压开采海域天然气水合物排采和生产优化的先决条件。许多学者对天然气水合物井筒多相流动规律进行了研究。孙宝江等[9]考虑深水钻井外部的多温度梯度环境和天然气水合物相变,建立了七组分井筒多相流控制方程;WEI等[10]研究了海洋天然气水合物储层钻井过程中井筒内瞬态气液两相流动变化机理;ZHANG、DENG等[11-13]研究了天然气水合物在深水气井测试过程中的井筒多相流规律,建立了井筒多相流模型,预测水合物在井筒内的生成风险;HIROBAYASHI等[14]通过实验的方法研究了甲烷水合物生产系统中的气液两相流动。

    目前针对降压开采海域天然气水合物排采工艺及生产优化的研究处于起步阶段。为实现天然气水合物储层长期试采的目标,需要进行天然气水合物降压开采人工举升适应性分析,明确压降控制机制并进行生产优化,避免井筒积液或产气管线连续产水等复杂工况的出现,实现天然气水合物的连续排采。笔者对降压开采海域天然气水合物电潜泵排采的可行性进行论证,并基于日本第1次试采过程中复杂的分采管柱和复杂的传热过程,建立井筒气液两相流模型,预测不同管线的温度和压力剖面,并利用该模型进行生产优化,为降压开采海域天然气水合物电潜泵排采生产优化的研究提供了理论支撑。

    • 开采方式是天然气水合物进行举升工艺适应性分析的基础。天然气水合物常见的开采方式包括降压法、注热法、CO2置换法、固态流化法等,其中降压法通过降低天然气水合物储层压力至相平衡压力之下而引起水合物相变分解,是一种较为简单且经济有效的开采方式[2]。根据全球天然气水合物的9次试采情况统计(见表1)[4-7],日本南海海槽的2次降压试采和中国南海神狐海域的2次降压开采全部采用电潜泵,利用泵吸作用抽汲地层与井中的水来降低水位和水合物储层的压力,打破相平衡条件,使水合物从固态变成气体和水,实现天然气水合物的成功试采。日本和中国的4次试采实践证明了电潜泵排水采气工艺适用于降压开采海域天然气水合物。

      表 1  全球天然气水合物试采情况

      Table 1.  Statistics of global pilot production of natural gas hydrates

      时间地点试采方法举升工艺
      2002年加拿大麦肯齐三角州注热法电潜泵
      2007年加拿大麦肯齐三角州降压法电潜泵
      2008年加拿大麦肯齐三角州降压法电潜泵
      2012年美国阿拉斯加北坡CO2置换法射流泵
      2013年日本南海海槽降压法电潜泵
      2017年日本南海海槽降压法电潜泵
      2017年中国南海神狐海域降压法电潜泵
      2017年中国南海荔湾海域固态流化法连续油管
      2020年中国南海神狐海域降压法电潜泵
    • 天然气水合物降压开采后,分解成游离气相和游离的水相,然后流动到海上平台。根据国内外天然气水合物试采井的资料,降压开采海域天然气水合物的井底流压较低,水气比较高,无法自喷到平台,必须采取有效的人工举升工艺防止井筒积液,以保证天然气水合物试采井能够正常生产。常用的海上气井排水采气方式包括电潜泵、电潜螺杆泵、射流泵、气举等,不同排水采气方式在海域天然气水合物排采中的适应性分析见表2[15-16]

      表 2  天然气水合物试采井排水采气适应性分析

      Table 2.  Applicability analysis of dewatering gas recovery in pilot production wells of natural gas hydrates

      项目电潜泵电潜螺杆泵射流泵气举
      海上适应性良好良好良好良好
      最大下入
      深度/m
      27002 00028003000
      最大排液量/
      (m3·d−1)
      500100300400
      防砂能力较弱较强较弱
      应用情况日本、中国、加拿大美国

      根据日本第1次天然气水合物试采资料,平均排液量约为200 m3/d,超过了电潜螺杆泵排水采气工艺允许的排液量上限,只有气举、射流泵、电潜泵排水采气工艺满足天然气水合物降压开采的排液量上限。电潜螺杆泵适用于各种井型的复杂工况,是成熟的排水采气技术,不会发生气锁,较适合于液气混输,但其排量变化范围小,常规电潜螺杆泵高造斜井段通过能力差,不适合用于水合物试采。

      射流泵的地面配套动力设备复杂,泵效率低,容易发生气蚀,要求较高的吸入压力和较高的沉没度。气举方式需要气源供气,即便使用自产天然气作为气源,也需要先建成一个气源启动气举过程,而且气举方式不利于天然气取样监测,制约条件多,装置设备多。因此,射流泵和气举不适合在水合物试采试验中采用。

      电潜泵排水采气工艺在日本、中国、加拿大的7次水合物试采中都发挥了重要作用,虽然射流泵在美国CO2置换试验中应用过,但泵效低,试验中产气量仅530 m3/d,而电潜螺杆泵、机抽和气举都没有成功应用的案例及可以借鉴的经验。虽然电潜泵防砂能力较弱,但利用电潜泵能够扩大生产压差,实现天然气水合物井的连续排采。试采期间如果出现砂堵,可使用连续油管冲砂,电潜泵无法运转时可使用气举进行助产。

      综上所述,对比各种排水采气工艺的技术适应性及在水合物试采中的应用情况,优选电潜泵作为水合物试采人工举升方式。

    • 日本海域第1次天然气水合物试采过程中采用了电潜泵系统来降低井筒内静液柱压力,管柱结构复杂(图1),包括电潜泵、分离器、加热器、电潜泵封隔器等。天然气水合物降压分解后,变成气液两相流,从井底向上流动。混合流体在向上过程中,流经加热器和气体分离器,分离器将气体和水进行分离,气体进入电潜泵封隔器下方的环空,与来自储层的混合水和气体相结合,在油套环空中向上流动然后流入封隔器上方的油管中到达海底测试树,最终通过隔水管中的钻杆进入钻井平台;水经过电潜泵加压后进入封隔器上方的油套环空,最终进入与防喷器相连的节流管线,流入到钻井平台。这种设计有利于通过调节电潜泵频率和井口回压来调节采气管线中的液位从而实现对压降的控制,防止水合物的二次生成[6]

      图  1  天然气水合物试采系统示意图

      Figure 1.  Schematic diagram of the pilot production system of natural gas hydrates

      降压开采海域天然气水合物电潜泵排水采气,与常规采油井、常规气井电潜泵排水采气相比,流动通道和流动规律更复杂,要考虑不同管线的流动规律分别建立模型。由于日本第1次天然气水合物试采中加热器发生电缆故障,没有发挥加热作用,因此,在此模型中忽略了电加热器的影响。基于日本第1次试采井资料,综合考虑复杂的海水环境和电潜泵影响,将流动通道简化处理成5段:(1)电潜泵以下的油管为气液混输管线;(2)电潜泵到海底测试树之间的油套环空为采气管线;(3)电潜泵到海底泥线之间的油管为排水管线;(4)海底泥线到钻井平台之间的隔水管内的钻杆为采气管线;(5)海底泥线到平台之间的隔水管外的节流管线为排水管线。

    • 综合考虑了复杂的分采管柱、周围海水环境的温度场、电潜泵、分离器等多种因素的影响,根据质量守恒原理,分别建立了混输管线、采气管线、排水管线的气液两相流的连续性方程。

      (1)混输管线。

      $$ 气相\quad\quad\quad\frac{\partial \left (A{E}_{\mathrm{g}}{\rho }_{\mathrm{g}}\right)}{\partial t}+\frac{\partial \left (A{E}_{\mathrm{g}}{\rho }_{\mathrm{g}}{\nu }_{\mathrm{g}}\right)}{\partial {\textit{z}}}={{\rho }_{\mathrm{g}}Q}_{\mathrm{g}\mathrm{T}} $$ (1)
      $$ 液相 \quad\quad \frac{\partial \left (A{E}_{\mathrm{w}}{\rho }_{\mathrm{w}}\right)}{\partial t}+\frac{\partial \left (A{E}_{\mathrm{w}}{\rho }_{\mathrm{w}}{\nu }_{\mathrm{w}}\right)}{\partial {\textit{z}}}={{\rho }_{\mathrm{w}}Q}_{\mathrm{w}\mathrm{T}} $$ (2)

      (2)采气管线。

      $$ 气相\; \frac{\partial \left (A{E}_{\mathrm{g}}{\rho }_{\mathrm{g}}\right)}{\partial t}+\frac{\partial \left (A{E}_{\mathrm{g}}{\rho }_{\mathrm{g}}{\nu }_{\mathrm{g}}\right)}{\partial {\textit{z}}} = {{{\rho }_{\mathrm{g}}Q}_{\mathrm{g}\mathrm{T}} - \rho }_{\mathrm{g}}{Q}_{\mathrm{p}}(1 - \lambda ) $$ (3)
      $$ 液相\;\;\; \frac{\partial \left (A{E}_{\mathrm{w}}{\rho }_{\mathrm{w}}\right)}{\partial t}+\frac{\partial \left (A{E}_{\mathrm{w}}{\rho }_{\mathrm{w}}{\nu }_{\mathrm{w}}\right)}{\partial {\textit{z}}}={{{\rho }_{\mathrm{w}}Q}_{\mathrm{w}\mathrm{T}}-\rho }_{\mathrm{w}}{Q}_{\mathrm{p}}\lambda $$ (4)

      (3)排水管线。

      $$ 气相 \quad\quad \frac{\partial \left (A{E}_{\mathrm{g}}{\rho }_{\mathrm{g}}\right)}{\partial t}+\frac{\partial \left (A{E}_{\mathrm{g}}{\rho }_{\mathrm{g}}{\nu }_{\mathrm{g}}\right)}{\partial {\textit{z}}}={\rho }_{\mathrm{g}}{Q}_{\mathrm{p}}(1-\lambda )$$ (5)
      $$ 液相 \quad\quad\quad \frac{\partial \left (A{E}_{\mathrm{w}}{\rho }_{\mathrm{w}}\right)}{\partial t}+\frac{\partial \left (A{E}_{\mathrm{w}}{\rho }_{\mathrm{w}}{\nu }_{\mathrm{w}}\right)}{\partial {\textit{z}}}={\rho }_{\mathrm{w}}{Q}_{\mathrm{p}}\lambda $$ (6)
    • 管线内的流体主要受重力、压力、摩擦力等,根据动量守恒原理,作用于单元体的合外力等于单元体内动量的变化率,建立动量守恒方程。

      $$ \begin{split} & \frac{\partial \left (A{E}_{\mathrm{g}}{\rho }_{\mathrm{g}}+A{E}_{\mathrm{w}}{\rho }_{\mathrm{w}}\right)}{\partial t}+\frac{\partial \left (A{E}_{\mathrm{g}}{\rho }_{\mathrm{g}}{\nu }_{\mathrm{g}}^{2}+A{E}_{\mathrm{w}}{\rho }_{\mathrm{w}}{\nu }_{\mathrm{w}}^{2}\right)}{\partial {\textit{z}}}+\\ & Ag\mathrm{c}\mathrm{o}\mathrm{s}\theta \left ({E}_{\mathrm{g}}{\rho }_{\mathrm{g}}+{E}_{\mathrm{w}}{\rho }_{\mathrm{w}}\right)+\frac{\mathrm{d}\left (Ap\right)}{\mathrm{d}{\textit{z}}}+A{f}_{\mathrm{r}}\frac{{\rho }_{\mathrm{m}}{\nu }_{\mathrm{m}}^{2}}{2d}=0 \end{split}$$ (7)
    • 不同管线的传热方式不同,其传热过程需要单独分析,5段流动通道的传热过程为:1)混输管线的换热过程主要包括套管内流体对流换热、套管壁导热、水泥环导热、地层非稳态传热;2)电潜泵到海底泥线之间的采气管线的换热过程主要包括排水管线内的流体与管线内壁的对流换热、排水管线的导热、排水管线与采气管线之间的环空的辐射换热与自然对流换热、套管壁导热、水泥环导热、地层非稳态传热;3)电潜泵到海底泥线之间的排水管线的换热过程与采气管线相同;4)海底泥线到平台之间的采气管线的换热过程主要包括采气管线内流体与采气管线之间的对流换热、采气管线的导热、隔水管与采气管线之间的流体对流换热与辐射换热、隔水管壁的导热、隔水管外壁与海水之间的对流换热;5)海底泥线到平台之间的排水管线的换热过程主要包括排水管线内流体与排水管线之间的对流换热、排水管线的导热、排水管线外壁和海水的对流换热。根据不同管线的传热过程建立对应的能量方程,并计算不同管线的综合传热系数。

      (1)混输管线(套管)。

      $$ \begin{split} & \mathrm{\pi }{r}_{\mathrm{m}\mathrm{i}}^{2}\left ({\rho }_{\mathrm{g}}{C}_{\mathrm{g}}{E}_{\mathrm{g}}+{\rho }_{\mathrm{w}}{C}_{\mathrm{w}}{E}_{\mathrm{w}}\right)\frac{{\partial T}_{\mathrm{m}}}{\partial t}+\\ &\left ({\rho }_{\mathrm{g}}{C}_{\mathrm{g}}{E}_{\mathrm{g}}+{\rho }_{\mathrm{w}}{C}_{\mathrm{w}}{E}_{\mathrm{w}}\right)\frac{{\partial T}_{\mathrm{m}}}{\partial {\textit{z}}}=\frac{2\mathrm{\pi }{r}_{\mathrm{m}\mathrm{o}}{U}_{\mathrm{m}}}{{{T}_{\mathrm{D}}U}_{\mathrm{m}}{r}_{\mathrm{m}\mathrm{o}}+{\lambda }_{\mathrm{e}}}\left ({T}_{\mathrm{e}\mathrm{i}}-{T}_{\mathrm{m}}\right) \end{split}$$ (8)
      $$ \frac{1}{{U}_{\mathrm{m}}}=\frac{{r}_{\mathrm{c}\mathrm{o}}}{{r}_{\mathrm{c}\mathrm{i}}{h}_{\mathrm{m}}}+\frac{\displaystyle\sum {r}_{\mathrm{c}\mathrm{o}}\mathrm{l}\mathrm{n}\frac{{r}_{\mathrm{c}\mathrm{o}}}{{r}_{\mathrm{c}\mathrm{i}}}}{{\lambda }_{\mathrm{c}\mathrm{a}\mathrm{s}}}+\frac{\displaystyle\sum {r}_{\mathrm{c}\mathrm{o}}\mathrm{l}\mathrm{n}\frac{{r}_{\mathrm{b}\mathrm{o}}}{{r}_{\mathrm{b}\mathrm{i}}}}{{\lambda }_{\mathrm{c}\mathrm{e}\mathrm{m}}} $$ (9)

      (2)电潜泵到海底泥线间采气管线(油套环空)。

      $$ \begin{split} & \mathrm{\pi }\left ({r}_{\mathrm{m}\mathrm{o}}^{2}-{r}_{\mathrm{w}\mathrm{o}}^{2}\right)\left ({\rho }_{\mathrm{g}}{C}_{\mathrm{g}}{E}_{\mathrm{g}}+{\rho }_{\mathrm{w}}{C}_{\mathrm{w}}{E}_{\mathrm{w}}\right)\frac{{\partial T}_{\mathrm{g}}}{\partial t}+\\ &\left ({\rho }_{\mathrm{g}}{C}_{\mathrm{g}}{E}_{\mathrm{g}}+{\rho }_{\mathrm{w}}{C}_{\mathrm{w}}{E}_{\mathrm{w}}\right)\frac{{\partial T}_{\mathrm{g}}}{\partial {\textit{z}}}=2\mathrm{\pi }{r}_{\mathrm{w}\mathrm{o}}{U}_{\mathrm{w}}\left ({T}_{\mathrm{w}}-{T}_{\mathrm{g}}\right)+\\ &\frac{2\mathrm{\pi }{\lambda }_{\mathrm{e}}{r}_{\mathrm{m}}{U}_{\mathrm{a}}}{{U}_{\mathrm{g}}{r}_{\mathrm{m}}{T}_{\mathrm{D}}+{\lambda }_{\mathrm{e}}}\left ({T}_{\mathrm{e}\mathrm{i}}-{T}_{\mathrm{g}}\right) \end{split}$$ (10)
      $$ \dfrac{1}{{U}_{\mathrm{a}}}=\dfrac{1}{{h}_{\mathrm{c}}+{h}_{\mathrm{r}}}+\dfrac{\displaystyle\sum {r}_{\mathrm{c}\mathrm{o}}\mathrm{l}\mathrm{n}\dfrac{{r}_{\mathrm{c}\mathrm{o}}}{{r}_{\mathrm{c}\mathrm{i}}}}{{\lambda }_{\mathrm{c}\mathrm{a}\mathrm{s}}}+\dfrac{\displaystyle\sum {r}_{\mathrm{c}\mathrm{o}}\mathrm{l}\mathrm{n}\dfrac{{r}_{\mathrm{b}\mathrm{o}}}{{r}_{\mathrm{b}\mathrm{i}}}}{{\lambda }_{\mathrm{c}\mathrm{e}\mathrm{m}}} $$ (11)
      $$ \frac{1}{{U}_{\mathrm{w}}}=\frac{{r}_{\mathrm{t}\mathrm{o}}}{{r}_{\mathrm{t}\mathrm{i}}{h}_{\mathrm{w}}}+\frac{{r}_{\mathrm{t}\mathrm{o}}\mathrm{l}\mathrm{n}\frac{{r}_{\mathrm{t}\mathrm{o}}}{{r}_{\mathrm{t}\mathrm{i}}}}{{\lambda }_{\mathrm{w}}} $$ (12)

      (3)电潜泵到海底泥线间的排水管线(油管)。

      $$\begin{split} & \mathrm{\pi }{r}_{\mathrm{w}\mathrm{i}}^{2}\left ({\rho }_{\mathrm{g}}{C}_{\mathrm{g}}{E}_{\mathrm{g}}+{\rho }_{\mathrm{w}}{C}_{\mathrm{w}}{E}_{\mathrm{w}}\right)\frac{{\partial T}_{\mathrm{w}}}{\partial t}+\\ &\left ({\rho }_{\mathrm{g}}{C}_{\mathrm{g}}{E}_{\mathrm{g}}+{\rho }_{\mathrm{w}}{C}_{\mathrm{w}}{E}_{\mathrm{w}}\right)\frac{{\partial T}_{\mathrm{w}}}{\partial {\textit{z}}}= \frac{2\mathrm{\pi }{r}_{\mathrm{w}\mathrm{o}}{U}_{\mathrm{w}}}{{{T}_{\mathrm{D}}U}_{\mathrm{w}}{r}_{\mathrm{w}\mathrm{o}}+{\lambda }_{\mathrm{e}}}\left ({T}_{\mathrm{e}\mathrm{i}}-{T}_{\mathrm{w}}\right) \end{split}$$ (13)
      $$ \dfrac{1}{{U}_{\mathrm{w}}}=\dfrac{{r}_{\mathrm{w}\mathrm{o}}}{{r}_{\mathrm{w}\mathrm{i}}{h}_{\mathrm{w}}} + \dfrac{{r}_{\mathrm{w}\mathrm{o}}\mathrm{l}\mathrm{n}\dfrac{{r}_{\mathrm{w}\mathrm{o}}}{{r}_{\mathrm{w}\mathrm{i}}}}{{\lambda }_{\mathrm{w}}} + \dfrac{1}{{h}_{\mathrm{c}} + {h}_{\mathrm{r}}} + \dfrac{\displaystyle\sum {r}_{\mathrm{w}\mathrm{o}}\mathrm{l}\mathrm{n}\dfrac{{r}_{\mathrm{c}\mathrm{o}}}{{r}_{\mathrm{c}\mathrm{i}}}}{{\mathrm{\lambda }}_{\mathrm{c}\mathrm{a}\mathrm{s}}} + \dfrac{\displaystyle\sum {r}_{\mathrm{w}\mathrm{o}}\mathrm{l}\mathrm{n}\dfrac{{r}_{\mathrm{b}\mathrm{o}}}{{r}_{\mathrm{b}\mathrm{i}}}}{{\lambda }_{\mathrm{c}\mathrm{e}\mathrm{m}}} $$ (14)

      (4)海底泥线到钻井平台间采气管线(钻杆)。

      $$ \begin{split} & \mathrm{\pi }{r}_{\mathrm{g}\mathrm{i}}^{2}\left ({\rho }_{\mathrm{g}}{C}_{\mathrm{g}}{E}_{\mathrm{g}}+{\rho }_{\mathrm{w}}{C}_{\mathrm{w}}{E}_{\mathrm{w}}\right)\frac{{\partial T}_{\mathrm{g}}}{\partial t}+\\ &\left ({\rho }_{\mathrm{g}}{C}_{\mathrm{g}}{E}_{\mathrm{g}}+{\rho }_{\mathrm{w}}{C}_{\mathrm{w}}{E}_{\mathrm{w}}\right)\frac{{\partial T}_{\mathrm{g}}}{\partial {\textit{z}}}= 2\mathrm{\pi }{r}_{\mathrm{g}\mathrm{o}}{U}_{\mathrm{g}}\left ({T}_{\mathrm{s}\mathrm{e}\mathrm{a}}-{T}_{\mathrm{g}}\right) \end{split} $$ (15)
      $$ \dfrac{1}{{U}_{\mathrm{g}}}=\dfrac{{r}_{\mathrm{g}\mathrm{o}}}{{r}_{\mathrm{g}\mathrm{i}}{h}_{\mathrm{g}}}+\dfrac{{r}_{\mathrm{g}\mathrm{o}}\mathrm{l}\mathrm{n}\dfrac{{r}_{\mathrm{g}\mathrm{o}}}{{r}_{\mathrm{g}\mathrm{i}}}}{{\lambda }_{\mathrm{g}}}+\dfrac{1}{{h}_{\mathrm{c}}+{h}_{\mathrm{r}}}+\dfrac{{r}_{\mathrm{g}\mathrm{o}}\mathrm{l}\mathrm{n}\dfrac{{r}_{\mathrm{r}\mathrm{o}}}{{r}_{\mathrm{r}\mathrm{i}}}}{{\lambda }_{\mathrm{r}\mathrm{i}\mathrm{s}}}+\dfrac{{r}_{\mathrm{g}\mathrm{o}}}{{r}_{\mathrm{g}\mathrm{i}}{h}_{\mathrm{r}\mathrm{s}}} $$ (16)

      (5)海底泥线到钻井平台间排水管线(节流管线)。

      $$\begin{split} & \pi {r}_{\mathrm{w}\mathrm{i}}^{2}\left ({\rho }_{\mathrm{g}}{C}_{\mathrm{g}}{E}_{\mathrm{g}}+{\rho }_{\mathrm{w}}{C}_{\mathrm{w}}{E}_{\mathrm{w}}\right)\frac{{\partial T}_{\mathrm{w}}}{\partial t}+\\ &\left ({\rho }_{\mathrm{g}}{C}_{\mathrm{g}}{E}_{\mathrm{g}}+{\rho }_{\mathrm{w}}{C}_{\mathrm{w}}{E}_{\mathrm{w}}\right)\frac{{\partial T}_{\mathrm{w}}}{\partial {\textit{z}}}=2\pi {r}_{\mathrm{w}\mathrm{o}}{U}_{\mathrm{w}}\left ({T}_{\mathrm{s}\mathrm{e}\mathrm{a}}-{T}_{\mathrm{w}}\right) \end{split}$$ (17)
      $$ \dfrac{1}{{U}_{\mathrm{w}}}=\dfrac{{r}_{\mathrm{w}\mathrm{o}}}{{r}_{\mathrm{w}\mathrm{i}}{h}_{\mathrm{w}}}+\dfrac{{r}_{\mathrm{w}\mathrm{o}}\mathrm{l}\mathrm{n}\dfrac{{r}_{\mathrm{w}\mathrm{o}}}{{r}_{\mathrm{w}\mathrm{i}}}}{{\lambda }_{\mathrm{w}}}+\dfrac{{r}_{\mathrm{w}\mathrm{o}}}{{r}_{\mathrm{w}\mathrm{i}}{h}_{\mathrm{w}\mathrm{s}}} $$ (18)
    • 天然气水合物降压开采井筒多相流计算过程选择从井底向上计算,井底温度等于周围地层的温度,井底流压假设保持不变。在电潜泵位置处,考虑了电潜泵的增压及电泵、电机产生的增温,忽略了电缆升温对井筒温度的影响。

      $$ {p}_{\mathrm{o}\mathrm{u}\mathrm{t}}={p}_{\mathrm{i}\mathrm{n}}+\Delta p $$ (19)
      $$ {T}_{\mathrm{o}\mathrm{u}\mathrm{t}}={T}_{\mathrm{i}\mathrm{n}}+\Delta T $$ (20)

      式中,$ A $为环空截面积,m2$ {E}_{\mathrm{g}} $$ {E}_{\mathrm{w}} $分别为持气率、持液率,无量纲;$ \;{\rho }_{\mathrm{g}} $$ {\rho }_{\mathrm{w}} $$ {\rho }_{\mathrm{m}} $分别为气体、液体、混合流体的密度,kg/m3$ {\nu }_{\mathrm{g}} $$ {\nu }_{\mathrm{l}} $$ {\nu }_{\mathrm{m}} $分别为气相、液相、混合流体的速度,m/s;$ {Q}_{\mathrm{g}\mathrm{T}} $$ {Q}_{\mathrm{w}\mathrm{T}} $分别为总产气量、总产水量,m3/s;$ {Q}_{\mathrm{P}} $为电潜泵排量,m3/s;$ \lambda $为气液分离效率;$ \theta $为井斜角,°;$ p $为压力,Pa;$ {f}_{\mathrm{r}} $为摩阻系数;$ d $为生产管柱的内径,m;$ {T}_{\mathrm{e}\mathrm{i}} $$ {T}_{\mathrm{s}\mathrm{e}\mathrm{a}} $为周围地层、海水温度,K;$ {U}_{\mathrm{m}} $$ {U}_{\mathrm{g}} $$ {U}_{\mathrm{w}} $为混输管线、采气管线、排水管线的综合传热系数,W/(m2·K);$ {U}_{\mathrm{a}} $为油管与套管之间的环空流体与地层的总传热系数W/(m2·K);$ {T}_{\mathrm{D}} $为无因次时间,无量纲;$ {\lambda }_{\mathrm{e}} $为地层导热系数,W/(m·K);$ {r}_{\mathrm{m}\mathrm{o}} $$ {r}_{\mathrm{m}\mathrm{i}} $$ {r}_{\mathrm{g}\mathrm{o}} $$ {r}_{\mathrm{g}\mathrm{i}} $$ {r}_{\mathrm{w}\mathrm{o}} $$ {r}_{\mathrm{w}\mathrm{i}} $$ {r}_{\mathrm{r}\mathrm{o}} $$ {r}_{\mathrm{r}\mathrm{i}} $分别为混输管线、采气管线、排水管线、隔水管的外半径和内半径,m;$ {r}_{\mathrm{c}\mathrm{o}} $$ {r}_{\mathrm{c}\mathrm{i}} $$ {r}_{\mathrm{b}\mathrm{o}} $$ {r}_{\mathrm{b}\mathrm{i}} $分别为套管、水泥环的外半径和内半径,m;$ {T}_{\mathrm{e}\mathrm{i}} $为地层温度,K;$ {T}_{\mathrm{m}} $$ {T}_{\mathrm{g}} $$ {T}_{\mathrm{w}} $分别为混输管线、采气管线、排水管线内的流体温度,K; $ {\lambda }_{\mathrm{g}} $$ {\lambda }_{\mathrm{w}} $$ {\lambda }_{\mathrm{r}\mathrm{i}\mathrm{s}} $分别为采气管线、排水管线、隔水管的导热系数,W/(m·K);$ {h}_{\mathrm{m}} $$ {h}_{\mathrm{g}} $$ {h}_{\mathrm{w}} $分别为流体与混输管线、采气管线、排水管线之间的对流换热系数,W/(m2·K);$ {h}_{\mathrm{c}} $$ {h}_{\mathrm{r}} $分别自然对流传热系数和辐射传热系数,W/(m2·K);$ {h}_{\mathrm{w}\mathrm{s}} $$ {h}_{rs} $为排水管线外壁、隔水管外壁与海水的对流换热系数,W/(m2·K);$ {p}_{\mathrm{o}\mathrm{u}\mathrm{t}} $$ {p}_{\mathrm{i}\mathrm{n}} $$ \Delta p $分别为电泵的出口、入口及增加的压力,Pa;$ {T}_{\mathrm{o}\mathrm{u}\mathrm{t}} $$ {T}_{\mathrm{i}\mathrm{n}} $$ \Delta T $分别为电泵的出口、入口及增加的温度,℃。

    • 根据日本第1次天然气水合物测试结果来看,试采初期,产量和压力都不稳定,这个时间段长达约18 h,而稳定生产阶段,井底流压在5 MPa左右,波动幅度小于0.5 MPa,波动较小,产气量和水气比保持稳定,总的产气量大约在20 000 m3/d,总的产水量大约在200 m3/d,稳定生产阶段约4 d[6]。模型基础参数如表3所示。

      表 3  模型基础参数

      Table 3.  Basic model parameters

      参数数值参数数值
      井深/m1 350平均日产气量/(m3·d−1)20 000
      海水深度/m1 000平均日产水量/(m3·d−1)200
      电潜泵和EIH深度/m1 200气液分离系数0.95
      地温梯度/(℃·(100 m)−1)2电加热器功率/kW30
      井底流压/MPa5电潜泵频率/Hz60

      日本第1次水合物试采时,稳定生产阶段时间最长,流压波动较小,井筒内没有水合物的二次生成,因此本文模型忽略了储层变化对井筒内水合物流动的影响和水合物在井筒中的生成和分解,重点解决天然气水合物降压分解成天然气和水后的井筒流动和生产优化问题,可简化为稳态模型[17],为天然气水合物排采优化提供理论基础。

    • 采用Hasan Kabir多相流计算模型,按照压力增量迭代的方法,从井底向上计算不同管线的温度和压力分布。气液两相流计算时首先要先判断流型,然后计算不同流型下相应的压降和综合传热系数等。不同管线的气液两相流求解流程见图2所示。

      图  2  降压开采天然气水合物井筒气液两相流模型求解流程

      Figure 2.  Workflow to solve the gas-liquid two-phase well flow model for depressurization recovery of natural gas hydrates

    • 利用上述模型可以得出电潜泵试采井的不同管线在稳定生产阶段的流型分布、温度和压力剖面、气体体积分数等(图3)。计算结果表明,混输管线为搅动流;采气管线在动液面以下为环状流,动液面以上为纯气相流;排水管线主要为泡状流,在井口附近出现段塞流。周围海水温度场对采气管线和排水管线的温度分布影响较大;电潜泵对采气管线和排水管线的流动规律影响很大,在电潜泵位置处温度、压力、气体体积分数等参数发生突变。

      图  3  不同管线的压力、温度和气体体积分数与井深的关系曲线

      Figure 3.  Pressure, temperature, and gas volume fractions vs. well depth for different pipes

    • 根据日本第1次天然气水合物试采结果对模型进行验证。将稳定生产阶段关键位置点的温度和压力模拟结果与现场监测数据进行对比(排水管线无温度监测点),压力误差在0.7 MPa以内,温度的误差在1.7 ℃以内,表明模型计算结果准确(表4)。

      表 4  模拟结果与现场监测数据对比

      Table 4.  Comparison between simulation result and field monitoring

      位置压力/MPa温度/℃
      模拟
      结果
      监测
      数据
      误差模拟
      结果
      监测
      数据
      误差
      泵吸入口4.14.70.618.318.50.2
      泵排出口14.114.80.710.311.51.2
      海底测试树
      (采气管线)
      3.33.80.513.515.21.7
      海底测试树
      (排水管线)
      12.212.50.315.3××
      采气管线井口
      (钻井平台)
      0.60.20.415.215.10.1
      排水管线井口
      (钻井平台)
      2.62.50.115.5××
        注:“×”表示无数据
    • 天然气水合物降压开采过程中,压降的控制是影响产气量的一个重要因素。井底流压过高,水合物分解速度降低,产气量降低,同时水合物在井底会二次生成堵塞井筒。因此,合理的压降控制有利于天然气水合物试采井的连续排采,防止采气管线连续出水,提高采收率。采气管线连续出水的原因主要包括两方面,一是采气管线内气体流速较高,携液能力较强;二是随着天然气水合物降压分解,采气管线内气体体积分数逐渐增高,气体的膨胀作用降低了液柱压力,动液面逐渐升高。通过控制井口压力和调节电潜泵的频率可有效控制动液面高度,防止采气管线连续出水,从而实现对压降的有效控制。

    • 影响天然气水合物电潜泵排采工艺的参数主要包括井底流压、井口压力、下泵深度及电潜泵的级数、气液分离效率等,其中生产过程中可以调节的参数主要为井口压力和电潜泵频率。本文利用降压开采天然气水合物电潜泵排采的井筒气液两相流模型,分析井口压力和电潜泵频率对天然气水合物井筒流动过程中的影响。

    • 当采气管线的井口压力增加后,动液面下降,能够降低采气管线连续出水的风险。从图4图5的模拟结果可以看出,当井口压力为0.6 MPa时,采气管线中的动液面位于井口,当井口压力继续升高为1 MPa时的动液面下降到274 m,1.5 MPa时的动液面下降到535 m。动液面的下降主要是由于井口压力的增加导致动液面上方气柱的压力增加,动液面处的压力增加。

      图  4  不同井口压力下的生产管线的压力

      Figure 4.  Pressure of production pipe at different wellhead pressures

      图  5  不同井口压力下生产管线的气体体积分数

      Figure 5.  Gas volume fraction of production pipe at different wellhead pressures

    • 图6图7的模拟结果可以看出,当电潜泵的频率增加后,排水管线的产水量增加,采气管线中的动液面下降,井底流压下降,但较低的转速能够降低气体聚集对电潜泵性能的影响。模拟结果表明,当电潜泵频率等于57 Hz时,采气管线中的动液面正好等于0,当电潜泵频率低于57 Hz时,生产压差减小,采气管线井口连续出水的风险增大。日本第1次天然气水合物试采过程中,电潜泵频率增大到70 Hz后,气体大量进入产水管线,电潜泵性能下降。利用商业化采油采气优化设计软件PROPSER模拟不同CENTRILIFT GC1600泵在不同频率下的效率(图8),证实了电潜泵频率大于70 Hz后,泵效急剧下降。因此,电潜泵的频率应该保持在57~70 Hz范围内。

      图  6  不同电潜泵频率下生产管线的压力

      Figure 6.  Pressure of production pipe at different ESP frequencies

      图  7  不同电潜泵频率下生产管线的气体体积分数

      Figure 7.  Gas volume fraction of production pipe at different ESP frequencies

      图  8  CENTRILIFT GC1600电潜泵频率与泵效关系

      Figure 8.  Frequency vs. frequency of the CENTRILIFT GC1600 ESP

    • (1)海域天然气水合物试采井具有低压、大液量等特点,电潜泵排液量范围大,占地面积小,操作简单,适用于降压开采海域天然气水合物排水采气,试采过程中如果出现砂堵,可以使用气举进行助产。

      (2)在降压开采海域天然气水合物电潜泵排采过程中,稳定生产阶段水气同采,混输管线和采气管线为搅动流,排水管线主要为泡状流。周围海水温度场对采气管线和排水管线的温度分布影响较大,电潜泵对不同管线的温度、压力和气体体积分数的分布影响很大。

      (3)天然气水合物试采过程中压降的控制可以通过调节井口压力和电潜泵频率实现,通过增大井口压力或在合理范围内增加电潜泵频率可以降低采气管线中的动液面高度,从而降低采气管线连续出水的风险,延长稳产时间。

参考文献 (17)

目录

    /

    返回文章
    返回