水泥浆流变性对水平井固井顶替界面的影响——基于天河一号大规模集群计算平台的数值模拟

陶谦 周仕明 张晋凯 方春飞 薛亮 陈博

陶谦, 周仕明, 张晋凯, 方春飞, 薛亮, 陈博. 水泥浆流变性对水平井固井顶替界面的影响——基于天河一号大规模集群计算平台的数值模拟[J]. 石油钻采工艺, 2017, 39(2): 185-191. doi: 10.13639/j.odpt.2017.02.011
引用本文: 陶谦, 周仕明, 张晋凯, 方春飞, 薛亮, 陈博. 水泥浆流变性对水平井固井顶替界面的影响——基于天河一号大规模集群计算平台的数值模拟[J]. 石油钻采工艺, 2017, 39(2): 185-191. doi: 10.13639/j.odpt.2017.02.011
TAO Qian, ZHOU Shiming, ZHANG Jinkai, FANG Chunfei, XUE Liang, CHEN Bo. Effect of rheological property of slurry on cementing displacement interface of horizontal well: the numerical simulation based on large-scale cluster computing platform Tianhe-1[J]. Oil Drilling & Production Technology, 2017, 39(2): 185-191. doi: 10.13639/j.odpt.2017.02.011
Citation: TAO Qian, ZHOU Shiming, ZHANG Jinkai, FANG Chunfei, XUE Liang, CHEN Bo. Effect of rheological property of slurry on cementing displacement interface of horizontal well: the numerical simulation based on large-scale cluster computing platform Tianhe-1[J]. Oil Drilling & Production Technology, 2017, 39(2): 185-191. doi: 10.13639/j.odpt.2017.02.011

水泥浆流变性对水平井固井顶替界面的影响——基于天河一号大规模集群计算平台的数值模拟

doi: 10.13639/j.odpt.2017.02.011
基金项目: 国家科技重大专项“高压低渗油气藏固井完井技术”(编号:2016ZX05021-005)和国家科技重大专项“彭水地区高效钻井及压裂工程工艺优化技术”(编号:2016ZX05061-002)
详细信息
    作者简介:

    陶谦(1982-),2010年毕业于中国石油大学(北京)油气井工程专业,获工学博士学位,现从事固井新技术的研究工作,高级工程师。通讯地址:(100101)北京市朝阳区北辰东路8号北辰时代大厦520。电话:010-84988267。E-mail:taoqian.sripe@sinopec.com

  • 中图分类号: TE256

Effect of rheological property of slurry on cementing displacement interface of horizontal well: the numerical simulation based on large-scale cluster computing platform Tianhe-1

  • 摘要: 水平封固段的固井顶替效果是影响水平井固井质量的关键问题,套管居中度与水泥浆流变性是影响水平段固井顶替效果的重要参数。基于国家超级计算中心天河一号大规模集群计算平台,采用Fluent软件对300 m长水平段偏心环空进行1 000 s的固井顶替数值模拟,研究水泥浆流变性对固井顶替界面的影响规律。研究结果表明:在低居中度(居中度≤50.0%)条件下,水泥浆流性指数n对顶替的影响可以忽略;在中等居中度(50.0% < 居中度 < 85.0%)条件下,提高水泥浆n值可以大幅改善环空宽边隔离液滞留,降低水泥浆n值可以改善环空窄边隔离液滞留,设计时需要综合考虑水泥浆n值;在高居中度(85.0%≤居中度 < 100%)与理想居中度(居中度100%)条件下,设计水泥浆n值时要避免出现其值过小的情况,从而避免在环空出现水泥浆分层窜流现象。研究结果对固井顶替施工设计具有重要指导作用。
  • 图  1  水平段环空几何模型

    Figure  1.  Geometric model of annulus in horizontal section

    图  2  水平段环空网格模型

    Figure  2.  Grid model of annulus in horizontal section

    图  3  低居中度条件下不同水泥浆n值的顶替云图

    Figure  3.  Displacement cloud chart for different slurry n at low centralized degree

    图  4  低居中度条件下水泥浆n值对顶替界面长度与顶替效率的影响

    Figure  4.  Effect of slurry n on displacement interface length and displacement efficiency at low centralized degree

    图  5  中等居中度条件下不同水泥浆n值的顶替云图

    Figure  5.  Displacement cloud chart for different slurry n at moderate centralized degree

    图  6  中等居中度条件下水泥浆n值对顶替界面长度与顶替效率的影响

    Figure  6.  Effect of slurry n on displacement interface length and displacement efficiency at moderate centralized degree

    图  7  高居中度条件下不同水泥浆n值的顶替云图

    Figure  7.  Displacement cloud chart for different slurry n at high centralized degree

    图  8  高居中度条件下水泥浆n值对顶替界面长度与顶替效率的影响

    Figure  8.  Effect of slurry n on displacement interface length and displacement efficiency at high centralized degree

    图  9  理想居中度条件下不同水泥浆n值的顶替云图

    Figure  9.  Displacement cloud chart for different slurry n at ideal centralized degree

    图  10  理想居中度条件下水泥浆n值对顶替界面长度与顶替效率的影响

    Figure  10.  Effect of slurry n on displacement interface length and displacement efficiency at ideal centralized degree

    表  1  水泥浆顶替隔离液的其他参数

    Table  1.   Other parameter of slurry displacing spacer fluid

    下载: 导出CSV
  • [1] 陶文铨.数值传热学[M]. 2版.西安: 西安交通大学出版社, 2001: 351-352.

    TAO Wenquan. Numerical heat transfer[M]. 2nd. Xi'an: Xi'an Jiaotong University Press, 2001: 351-352.
    [2] 王福军.计算流体动力学[M].北京: 清华大学出版社, 2004: 39-43.

    WANG Fujun. Computational fluid dynamics[M]. Beijing: Tsinghua University Press, 2004: 39-43.
    [3] 郑永刚.非牛顿流体流动理论及其在石油工程中的应用[M].北京: 石油工业出版社, 1999: 39-88.

    ZHENG Yonggang. The theory of non Newtonian fluid flow and its application in petroleum engineering[M]. Beijing: Petroleum Industry Press, 1999: 39-88.
    [4] BITTLESTON S H, FERGUSON J, FRIGAARD I A. Mud removal and cement placement during primary cementing of an oil well—Laminar non-Newtonian displacements in an eccentric annular Hele-Shaw cell[J]. Journal of Engineering Mathematics, 2002, 43(2-4): 229-253.
    [5] PELIPENKO S, FRIGAARD I A. Mud removal and cement placement during primary cementing of an oil well—Part 2; steady-state displacements[J]. Journal of Engineering Mathematics, 2004, 48(1): 1-26.
    [6] 冯福平, 艾池, 杨丰宇, 孟翔, 钟鑫.偏心环空层流顶替滞留层边界位置研究[J].石油学报, 2010, 31(5): 858-862.

    FENG Fuping, AI Chi, YANG Fengyu, MENG Xiang, ZHONG Xin. A study on the replacement position of laminar flow in eccentric annulus at retention layer boundaries[J]. Acta Petrolei Sinica, 2010, 31(5): 858-862.
    [7] 王瑞和, 李明忠, 王成文, 方群.油气井注水泥顶替机理研究进展[J].天然气工业, 2013, 33(5): 69-76.

    WANG Ruihe, LI Mingzhong, WANG Chengwen, FANG Qun. Research progress in the cementing displacement mechanism[J]. Natural Gas Industry. 2013, 33(5): 69-76.
    [8] 李明忠, 王瑞和, 王成文, 方群.层流顶替钻井液滞留层厚度研究[J].水动力学研究与进展A辑, 2013, 28(5): 531-537.

    LI Mingzhong, WANG Ruihe, WANG Chengwen, FANG Qun. A study on drilling mud stagnant layer thickness under the condition of laminar displacement. Chinese Journal of Hydrodynamics, 2013, 28(5): 531-537.
    [9] 冯福平, 艾池, 于法浩, 崔志华, 徐海粟, 张嵇南.基于Hele-shaw模型的斜井偏心环空顶替流体密度差优化[J].数学的实践与认识, 2014, 44(22): 125-133.

    FENG Fuping, AI Chi, YU Fahao, CUI Zhihua, XU Haisu, ZHANG Jinan. The optimization of density difference in the eccentric annulus of inclined well on the basis of Heleshaw model[J]. Mathematics in Practice and Theory, 2014, 44(22): 125-133.
    [10] ADRIAN R J. Particle-imaging techniques for experimental fluid mechanics[J]. Annual Review of Fluid Mechanics, 1991, 23(1): 261-304. doi:  10.1146/annurev.fl.23.010191.001401
    [11] 郑永刚, 方铎, 郝俊芳.水平井注水泥的理论与实验研究[J].水动力学研究与进展(A辑), 1996, 11(1): 19-23.

    ZHENG Yonggang, FANG Duo, HAO Junfang. A theoretical and experimental study of cementing for horizontal well[J]. Journal of Hydrodynamics(edition A), 1996, 11(1): 19-23.
    [12] MALEKMOHAMMADI S, CARRASCO-TEJA M, STOREY S, FRIGAARD I A, MARTINEZ D M. An experimental study of laminar displacement flows in narrow vertical eccentric annuli[J]. Journal of Fluid Mechanics, 2010, 649(6): 371-398.
    [13] 朱亮, 徐璧华, 谭元铭, 冯青豪.水平井注水泥顶替效率模拟方法研究[J].西部探矿工程, 2015, 27(11): 32-34.

    ZHU Liang, XU Bihua, TAN Yuanming, FENG Qinghao. Study on simulation method of cement displacement efficiency in horizontal wells[J]. West China Exploration Engineering, 2015, 27(11): 32-42.
    [14] 张松杰, 薛亮, 汪志明, 张文华.密度差对水平井固井顶替影响规律数值模拟研究[J].钻采工艺, 2012, 35(6): 15-17.

    ZHANG Songjie, XUE Liang, WANG Zhiming, ZHANG Wenhua. Numerical simulation study on the influence law of density difference on cementing displacement interface in the horizontal well[J]. Drilling & Production Technology, 2012, 35(6): 15-17.
    [15] 冯福平, 艾池, 崔志华, 于法浩, 陈顶峰, 徐海粟.水平井偏心环空顶替流体密度差优化[J].石油钻采工艺, 2014, 36(1): 61-65.

    FENG Fuping, AI Chi, CUI Zhihua, YU Fahao, CHEN Dingfeng, XU Haisu. Density difference optimization of displacing fluids in eccentric annulus of horizontal wells[J]. Oil Drilling & Production Technology, 2014, 36(1): 61-65.
    [16] 罗恒荣, 张晋凯, 周仕明, 陶谦, 方春飞.偏心度和密度差耦合条件下水平井顶替界面特征研究[J].石油钻探技术, 2016, 44(4): 65-71.

    LUO Hengrong, ZHANG Jinkai, ZHOU Shiming, TAO Qian, FANG Chunfei. Horizontal well cement displacement interface features under the coupling of eccentricity and density difference[J]. Petroleum Drilling Techniques, 2016, 44(4): 65-71.
    [17] 步玉环, 郭胜来, 马明新, 李红.复杂井眼条件下旋转套管速度对固井质量的影响[J].石油学报, 2011, 32(3): 529-533.

    BU Yuhuan, GUO Shenglai, MA Mingxin, LI Hong. Influence of rotating casing speed on well-cementing quality under the bad hole conditions[J]. Acta Petrolei Sinica, 2011, 32(3): 529-533.
    [18] 高永海, 孙宝江, 赵欣欣, 刘东清.前置液流变性对顶替界面稳定性影响的数值模拟[J].中国石油大学学报(自然科学版), 2007, 31(6): 51-54.

    GAO Yonghai, SUN Baojiang, ZHAO Xinxin, LIU Dongqing. Numerical simulation on influence of rheological property of front fluid on stability of cement displacement interface[J]. Journal of China University of Petroleum(Edition of Natural Science), 2007, 31(6): 51-54.
  • [1] 甘泉泉, 刘少胡, 管锋, 李明红.  水平井段工具内脉冲射流冲砂技术 . 石油钻采工艺, 2019, 41(4): 549-554. doi: 10.13639/j.odpt.2019.04.023
    [2] 郑忠茂, 刘天恩, 周宝义, 张铁.  浮力作用对大斜度井套管居中度的影响 . 石油钻采工艺, 2017, 39(3): 313-318. doi: 10.13639/j.odpt.2017.03.011
    [3] 冯福平, 艾池, 崔志华, 于法浩, 陈顶峰, 徐海粟.  水平井偏心环空顶替流体密度差优化 . 石油钻采工艺, 2014, 36(1): 61-65.
    [4] 徐鑫, 魏新芳, 王书峰.  水平井固井液压扶正器的研究与应用 . 石油钻采工艺, 2012, 34(2): 102-103.
    [5] 朱礼平, 廖忠会, 刁素, 严焱诚.  川西水平井水泥浆防气窜评价方法 . 石油钻采工艺, 2012, 34(1): 60-62.
    [6] 帅春岗, 段永刚, 冉林, 彭建云, 方全堂.  水平井ICD控水方法研究 . 石油钻采工艺, 2012, 34(1): 85-88.
    [7] 苏博勇, 姜宏图, 马淑梅, 霍维维.  胶乳水泥浆体系在深井水平井的应用 . 石油钻采工艺, 2010, 32(5): 53-55.
    [8] 王合林, 钟福海.  华北油田水平井完井固井技术 . 石油钻采工艺, 2009, 31(4): 113-117.
    [9] 汪志明, 洪玉奎, 郭晓乐, 李明.  泡沫质量和环空返速对泡沫携岩能力影响的数值模拟 . 石油钻采工艺, 2008, 30(3): 34-37.
    [10] 李兆敏, 李松岩, 尚朝辉, 李勇.  流速和环空偏心对泡沫携砂能力的影响数值模拟 . 石油钻采工艺, 2007, 29(3): 97-100. doi: 10.3969/j.issn.1000-7393.2007.03.031
    [11] 周长静, 李华斌, 苏国辉, 闫志强, 吕海燕.  渤海绥中油田J3井疏水缔合聚合物驱数值模拟研究 . 石油钻采工艺, 2006, 28(5): 31-33. doi: 10.3969/j.issn.1000-7393.2006.05.011
    [12] 王炳印, 卢晓峰, 曲从峰, 邓金根, 邹灵战.  疏松砂岩储层水平井安全钻井周期数值模拟计算 . 石油钻采工艺, 2005, 27(3): 11-12. doi: 10.3969/j.issn.1000-7393.2005.03.004
    [13] 董长银, 张琪, 曲占庆.  水平井砾石充填参数优化设计的正交数值试验法 . 石油钻采工艺, 2004, 26(2): 42-44. doi: 10.3969/j.issn.1000-7393.2004.02.013
    [14] 丁士东.  塔河油田紊流、塞流复合顶替固井技术 . 石油钻采工艺, 2002, 24(1): 20-22. doi: 10.3969/j.issn.1000-7393.2002.01.007
    [15] 杨启贞.  在水平井、大位移井中选用合适的套管扶正器 . 石油钻采工艺, 2000, 22(4): 37-39. doi: 10.3969/j.issn.1000-7393.2000.04.011
    [16] 黄学军, 崔桂香, 周海兵.  水平井变质量管内流动损失的数值研究(二):多支管流动 . 石油钻采工艺, 2000, 22(3): 50-53. doi: 10.3969/j.issn.1000-7393.2000.03.015
    [17] 周海兵, 崔桂香, 黄学军.  水平井变质量管内流动损失的数值研究(一):单支孔流动 . 石油钻采工艺, 2000, 22(2): 42-45. doi: 10.3969/j.issn.1000-7393.2000.02.012
    [18] 吴晓东, 李兆文, 王杰祥, 冯建华, 周赤峰.  定向井有杆抽油系统的数值模拟 . 石油钻采工艺, 1998, 20(4): 54-58. doi: 10.3969/j.issn.1000-7393.1998.04.010
    [19] 赵达壮, 张允真, 张红.  水平井井眼轨迹及套管初弯曲的计算 . 石油钻采工艺, 1994, 16(3): 5-11. doi: 10.3969/j.issn.1000-7393.1994.03.002
    [20] 杨全盛, 陈道元, 俞平方.  卫2—25水平井固井技术 . 石油钻采工艺, 1992, 14(4): 31-34. doi: 10.3969/j.issn.1000-7393.1992.04.005
  • 加载中
图(10) / 表 (1)
计量
  • 文章访问数:  466
  • HTML全文浏览量:  50
  • PDF下载量:  694
  • 被引次数: 0
出版历程
  • 修回日期:  2017-02-19
  • 刊出日期:  2017-03-20

水泥浆流变性对水平井固井顶替界面的影响——基于天河一号大规模集群计算平台的数值模拟

doi: 10.13639/j.odpt.2017.02.011
    基金项目:  国家科技重大专项“高压低渗油气藏固井完井技术”(编号:2016ZX05021-005)和国家科技重大专项“彭水地区高效钻井及压裂工程工艺优化技术”(编号:2016ZX05061-002)
    作者简介:

    陶谦(1982-),2010年毕业于中国石油大学(北京)油气井工程专业,获工学博士学位,现从事固井新技术的研究工作,高级工程师。通讯地址:(100101)北京市朝阳区北辰东路8号北辰时代大厦520。电话:010-84988267。E-mail:taoqian.sripe@sinopec.com

  • 中图分类号: TE256

摘要: 水平封固段的固井顶替效果是影响水平井固井质量的关键问题,套管居中度与水泥浆流变性是影响水平段固井顶替效果的重要参数。基于国家超级计算中心天河一号大规模集群计算平台,采用Fluent软件对300 m长水平段偏心环空进行1 000 s的固井顶替数值模拟,研究水泥浆流变性对固井顶替界面的影响规律。研究结果表明:在低居中度(居中度≤50.0%)条件下,水泥浆流性指数n对顶替的影响可以忽略;在中等居中度(50.0% < 居中度 < 85.0%)条件下,提高水泥浆n值可以大幅改善环空宽边隔离液滞留,降低水泥浆n值可以改善环空窄边隔离液滞留,设计时需要综合考虑水泥浆n值;在高居中度(85.0%≤居中度 < 100%)与理想居中度(居中度100%)条件下,设计水泥浆n值时要避免出现其值过小的情况,从而避免在环空出现水泥浆分层窜流现象。研究结果对固井顶替施工设计具有重要指导作用。

English Abstract

陶谦, 周仕明, 张晋凯, 方春飞, 薛亮, 陈博. 水泥浆流变性对水平井固井顶替界面的影响——基于天河一号大规模集群计算平台的数值模拟[J]. 石油钻采工艺, 2017, 39(2): 185-191. doi: 10.13639/j.odpt.2017.02.011
引用本文: 陶谦, 周仕明, 张晋凯, 方春飞, 薛亮, 陈博. 水泥浆流变性对水平井固井顶替界面的影响——基于天河一号大规模集群计算平台的数值模拟[J]. 石油钻采工艺, 2017, 39(2): 185-191. doi: 10.13639/j.odpt.2017.02.011
TAO Qian, ZHOU Shiming, ZHANG Jinkai, FANG Chunfei, XUE Liang, CHEN Bo. Effect of rheological property of slurry on cementing displacement interface of horizontal well: the numerical simulation based on large-scale cluster computing platform Tianhe-1[J]. Oil Drilling & Production Technology, 2017, 39(2): 185-191. doi: 10.13639/j.odpt.2017.02.011
Citation: TAO Qian, ZHOU Shiming, ZHANG Jinkai, FANG Chunfei, XUE Liang, CHEN Bo. Effect of rheological property of slurry on cementing displacement interface of horizontal well: the numerical simulation based on large-scale cluster computing platform Tianhe-1[J]. Oil Drilling & Production Technology, 2017, 39(2): 185-191. doi: 10.13639/j.odpt.2017.02.011
  • 水平井固井过程中, 水泥浆、隔离液、钻井液顺序顶替, 存在多个顶替界面, 每个顶替界面都会对固井质量产生重要影响, 深入研究水平井水平段顶替界面的形成过程与规律具有重要意义。目前固井顶替界面的研究方向主要分为3类:顶替界面的理论分析, 顶替界面的室内实验, 顶替界面的数值实验。顶替界面的理论分析以一维层流顶替模型[1-2]和二维Hele-Shaw模型为代表[3-9], 特别是二维Hele-Shaw模型可以很好地描述固井顶替流体界面的发展动态, 具有比较好的研究前景, 但建模和求解都具有相当大的难度。顶替界面的室内实验研究[10-13]在室内实验系统建设、实验相似液配置、顶替界面的显示与观测等方面都取得了很大的进步, 但井筒模拟长度限制了这一研究手段的使用。CFD数值模拟方法[14-18]以三维非定常组分多相流方程为控制方程, 通过计算机数值模拟固井流体的顶替界面发展过程, 能够真实再现三维顶替界面的形成与发展特征, 但需要采用大型机来进行。笔者基于国家超级计算中心天河一号大规模集群计算平台, 采用Fluent软件对300 m长水平段偏心环空进行1 000 s的固井顶替数值模拟, 研究了套管居中度与水泥浆流变性对固井顶替界面的影响规律, 为水平井固井顶替施工方案的合理设计提供理论基础。

    • 针对水平井水平段固井顶替问题, 建立水平井水平段偏心环空固井顶替物理模型, 基于国家超级计算中心天河一号大规模集群计算平台, 采用Fluent软件对流场进行求解, 开展水平井偏心环空固井顶替数值模拟实验。

    • 水平井井筒直径为215.9 mm, 考虑8%井径扩大率, 直径取233 mm; 套管外径取139.7 mm, 水平段环空长度为300 m。水平井由于套管不居中, 无法使用轴对称条件, 但可以采用左右对称边界条件。根据现场扶正器使用情况, 居中度设置为30%、50%、66.7%、85%、100%, 建立5个不同居中度的水平段偏心环空流场计算网格模型。图 1为300 m水平段偏心环空流场一段长度几何模型。因为采用左右对称条件, 所以取环空的一半进行计算。

      图  1  水平段环空几何模型

      Figure 1.  Geometric model of annulus in horizontal section

      采用数值模拟方法, 需要把连续解析的空间离散为一组有限个离散的点。不考虑扶正器与井眼的不规则性, 可以采用结构网格来布置流场空间离散点。300 m环空模型网格数量为600万网格点, 进行1 000 s的固井顶替计算。另外为了更好地捕捉顶替界面, 采用了网格局部加密技术进行流场计算。图 2为300 m水平段局部长度偏心环空网格分布图。图中截面分别为宽边截面、窄边截面、套管外壁面, 也就是环空的内壁面。壁面附近的网格进行了加密处理。

      图  2  水平段环空网格模型

      Figure 2.  Grid model of annulus in horizontal section

    • 水平井固井顶替过程中, 无论是水泥浆顶替隔离液还是隔离液顶替钻井液, 都属于长距离环空间隙中的液、液两相流动过程, 满足流体力学基本方程组。考虑到固井顶替流体的黏度等水力参数, 环空流态一般为非牛顿流体层流流动。流体力学基本方程组详见参考文献[1-2]。顶替界面采用组分模型进行捕捉。水泥浆、隔离液采用幂律非牛顿流变模式。

      1.2.1  顶替界面的组分方程  采用物质输运模型, 第i种组分的质量守恒方程用于计算第i种物质的质量分数Yi, 可用于捕捉两相界面, 以及两相界面附近的质量扩散。

      $$ \frac{\partial}{\partial t}\left(\rho Y_{i}\right)+\nabla \cdot\left(\rho V Y_{i}\right)=-\nabla J_{i}+R_{i}+S_{i} $$ (1)

      式中, t为时间, s; ρ为密度, kg/m3; V为矢量速度, m/s; Yi为组分i的质量分数; Ji为组分i的扩散通量, m3/s; Ri为系统内部单位时间内单位体积通过化学反应消耗或生成该种组分的净生成率, mol/(L·s); Si表示通过其他方式, 如异相反应、相变所生成该种组分的净生产率, 以及自定义的其他质量源项, mol/(L·s)。

      Yi=0时, 为一种固井流体; 当Yi=1.0时, 为另外一种固井流体; 当Yi为0~1.0时, 为两种流体顶替的混浆界面。

      1.2.2  初、边界条件  水平井偏心环空顶替流场包括入口边界、出口边界、环空内外壁面边界以及对称边界。

      (1) 初始条件:顶替流体初始时刻的顶替界面位置设置在300 m环空正中间150 m处, 初始时刻的速度均设为0。

      (2) 运动壁面边界条件:基于相对运动的概念, 环空壁面采用运动壁面边界条件。壁面运动速度设置为顶替流体平均速度, 运动方向与顶替流体流动方向相反。壁面上的流体随壁面一起反方向运动, 满足相对壁面无滑移的黏性流体边界条件。

      (3) 入口边界条件:与运动壁面边界条件配合, 环空入口平均速度设置为0。

      (4) 出口边界:利用相临内部流场节点数据进行外推插值计算, 满足质量守恒条件。

      (5) 对称边界:在半环空的宽边纵界面与窄边纵界面, 采用平面对称边界条件。

      基于相对运动的原理, 运动壁面边界条件与0速度入口边界条件配合使用, 2种顶替流体的界面可以一直停留流场中间, 顶替界面与混浆过程不断向两侧发展, 而不会很快到达出口边界。这样就可以采用较短的环空网格模型(300 m)计算很长时间(1000 s)的顶替时间, 而且还方便捕捉顶替界面随顶替时间的发展过程。

    • 水平井固井顶替中一般存在水泥浆顶替隔离液、隔离液顶替钻井液2个界面, 以水泥浆顶替隔离液为例进行顶替界面研究, 研究结果对隔离液顶替钻井液也具有一定的参考价值。

      根据固井顶替实践可知, 环空居中度、水泥浆流变性、隔离液流变性等对固井顶替界面具有重要影响。水平段套管无法严格居中会产生严重的宽窄边效应, 流体在环空宽边速度大于窄边速度, 导致水泥浆在环空宽边指进。

      水泥浆可根据需要来调整密度与流变性, 隔离液配置过程中一般对流性指数n值做重点考虑, 所以本文重点研究水泥浆n值对顶替界面的影响规律。居中度设置30%、50%、66.7%、85%、100%, 水泥浆n值设置0.5、0.6、0.7、0.8、0.9, 其他参数见表 1

      表 1  水泥浆顶替隔离液的其他参数

      Table 1.  Other parameter of slurry displacing spacer fluid

    • 居中度与水泥浆n值对固井顶替具有重要影响, 有必要对其进行耦合分析。在试算过程中, 发现不同居中度条件下水泥浆n值的影响存在差异, 所以为了全面研究水泥浆n值对固井顶替的影响, 在每个居中度条件下分别进行了5组水泥浆n值的数值模拟实验, 共25个算例。

      数值模拟过程保存了每秒的顶替界面特征, 限于篇幅仅以第100 s为例给出了不同居中度条件下的顶替界面体积分数云图。环空顶替界面云图长度方向按5:1进行缩比, 图示环空横截面间隔为5 m。红色区域代表水泥浆, 蓝色区域代表隔离液, 红色逐渐转为蓝色的区域为混浆界面。每个居中度条件下, 分别对比分析不同水泥浆n值条件下顶替界面体积分数云图, 可以直观获得居中度和水泥浆n值耦合对顶替界面的影响规律。

      2.2.1  低居中度  低居中度条件, 是指居中度30%和50%。图 3为居中度50%条件下的不同水泥浆n值第100 s顶替界面体积分数云图。在低居中度条件下, 水泥浆顶替界面前缘在环空宽边快速指进, 在环空窄边形成严重的隔离液滞留, 在环空宽边也会有很长的隔离液滞留, 质量扩散混浆现象明显, 顶替界面很长, 顶替效率很低。对比不同水泥浆n值在低居中度条件下顶替界面特征可以发现, 提高水泥浆n值可以小幅改善环空宽边隔离液滞留, 降低水泥浆n值只能微幅改善环空窄边隔离液滞留。或者说, 水泥浆n值对顶替界面发展影响很小, 可以忽略。

      图  3  低居中度条件下不同水泥浆n值的顶替云图

      Figure 3.  Displacement cloud chart for different slurry n at low centralized degree

      图 4为低居中度条件下水泥浆n值对顶替界面长度与顶替效率影响规律。水泥浆n值在50%居中度条件下对顶替界面长度与顶替效率影响很小。降低水泥浆n值, 顶替界面长度微幅减小, 顶替效率略有提高。n=0.8和n=0.9两条曲线在800 s以后顶替界面逐渐超出300 m计算域范围, 计算数据逐渐失真。

      图  4  低居中度条件下水泥浆n值对顶替界面长度与顶替效率的影响

      Figure 4.  Effect of slurry n on displacement interface length and displacement efficiency at low centralized degree

      综上, 在低居中度条件下, 水泥浆n值对顶替界面发展过程、顶替界面长度与顶替效率影响很小, 可以忽略, 即在低居中度条件下无法通过配置水泥浆n值来改善顶替效果。

      2.2.2  中等居中度  图 5为中等居中度(66.7%)条件下的不同水泥浆n值第100 s顶替界面体积分数云图, 水泥浆顶替界面前缘位于环空中部, 在环空窄边隔离液滞留不明显, 在环空宽边隔离液滞留也不明显, 质量扩散混浆现象明显, 顶替界面长度中等。对比不同水泥浆n值在中等居中度条件下顶替界面特征可以发现, 提高水泥浆n值可以改善环空宽边隔离液滞留, 降低水泥浆n值可以大幅改善环空窄边隔离液滞留, 并降低顶替界面长度。

      图  5  中等居中度条件下不同水泥浆n值的顶替云图

      Figure 5.  Displacement cloud chart for different slurry n at moderate centralized degree

      图 6为中等居中度条件下水泥浆n值对顶替界面长度与顶替效率影响规律。水泥浆n值在66.7%居中度条件下对顶替界面长度有一定影响, 顶替界面长度随水泥浆n值降低而减小, 顶替效率随水泥浆n值降低而提高。图 6中顶替界面长度与顶替效率随水泥浆n值变化曲线, 与图 5中的顶替界面云图特征一致, 进一步验证了水泥浆n值对顶替界面的影响规律。

      图  6  中等居中度条件下水泥浆n值对顶替界面长度与顶替效率的影响

      Figure 6.  Effect of slurry n on displacement interface length and displacement efficiency at moderate centralized degree

      综上, 在中等居中度条件下, 降低水泥浆n值有利于改善环空窄边隔离液滞留, 进一步降低顶替界面长度, 提高顶替效率。但提高水泥浆n值可以大幅改善环空宽边隔离液滞留, 所以在中等居中度条件下, 水泥浆n值要综合考虑。

      2.2.3  高居中度  图 7为高居中度(85%)条件下水泥浆n值对顶替界面的影响规律, 水泥浆顶替界面前缘位于环空中下部, 在环空窄边隔离液滞留很不明显, 混浆现象明显, 顶替界面环空宽边隔离液滞留也很不明显, 顶替界面长度较短。

      图  7  高居中度条件下不同水泥浆n值的顶替云图

      Figure 7.  Displacement cloud chart for different slurry n at high centralized degree

      图 7可知, 高居中度条件下提高水泥浆n值对顶替界面影响与其他居中度有重要区别。在水泥浆n值为0.5时, 由于水泥浆流动性好, 宽窄边效应不明显, 在正密度差(240 kg/m3)作用下, 水泥浆在环空出现了指进式窜流, 水泥浆与隔离液顶替界面出现了上下分层流动的现象, 导致顶替界面长度大幅提高。在水泥浆n值相对较大时(n=0.6, 0.7, 0.8, 0.9), 水泥浆n值对顶替界面影响与中等居中度条件下类似, 提高水泥浆n值有利于改善环空宽边隔离液滞留, 有利于减小顶替界面长度。

      图 8为高居中度条件下水泥浆n值对顶替界面长度与顶替效率影响规律, 与图 7中的顶替界面云图规律一致, 当水泥浆n值较小时(n=0.5), 由于水泥浆在环空出现窜流, 导致顶替界面很长, 顶替效率比较低。当水泥浆n值比较大时(n=0.6, 0.7, 0.8, 0.9), 顶替界面长度随水泥浆n值增加而增加, 顶替效率随水泥浆n值增加而下降。

      图  8  高居中度条件下水泥浆n值对顶替界面长度与顶替效率的影响

      Figure 8.  Effect of slurry n on displacement interface length and displacement efficiency at high centralized degree

      综上, 高居中度条件下, 设计水泥浆n值时要避免出现n值过小的情况, 避免在环空出现水泥浆窜流现象。在不出现水泥浆窜流的情况下, 减小水泥浆n值可以进一步降低顶替界面长度, 提高顶替效率。但由于高居中度条件下本身就具有良好的顶替效果, 所以调整水泥浆n值的意义有限。

      2.2.4  理想居中度图 9是理想居中度(100%)条件下水泥浆n值对顶替界面的影响规律, 水泥浆顶替界面前缘位于环空下部, 在环空窄边没有隔离液滞留, 在环空宽边隔离液滞留明显, 顶替界面在环空宽边存在质量扩散混浆现象, 顶替界面很短, 具有良好的顶替效率。

      图  9  理想居中度条件下不同水泥浆n值的顶替云图

      Figure 9.  Displacement cloud chart for different slurry n at ideal centralized degree

      对比高居中度条件下环空水泥浆窜流现象可以发现, 在理想居中度条件下水泥浆出现窜流的n值范围扩大, 由高居中度条件下的0.5扩大到0.6。在理想居中度条件下, 当水泥浆n值为0.5和0.6时, 由于水泥浆流动性相对较好, 不存在宽窄边效应, 在正密度差(240 kg/m3)等几个条件的综合作用下, 水泥浆在环空出现了指进式窜流, 水泥浆与隔离液顶替界面出现了上下分层流动的现象, 导致顶替界面长度大幅提高, 顶替效率大幅降低。

      图 10为理想居中度条件下水泥浆n值对顶替界面长度与顶替效率影响规律, 与图 9中的顶替界面云图规律一致, 当水泥浆n值较小时(n=0.5, 0.6), 由于水泥浆在环空出现窜流, 导致顶替界面很长, 顶替效率比较低。当水泥浆n值比较大时(n=0.7, 0.8, 0.9), 顶替界面长度随水泥浆n值降低而减小, 顶替效率随水泥浆n值减小而提高。

      图  10  理想居中度条件下水泥浆n值对顶替界面长度与顶替效率的影响

      Figure 10.  Effect of slurry n on displacement interface length and displacement efficiency at ideal centralized degree

      综上, 在理想居中度条件下, 设计水泥浆n值时要避免出现n值过小的情况, 避免在环空出现水泥浆窜流现象。在不出现水泥浆窜流的情况下, 由于理想居中度条件下本身就具有良好的顶替效果, 所以调整水泥浆n值的意义有限。

    • (1) 根据水平井水平段固井顶替数值模拟结果, 相对全面地研究了居中度与水泥浆n值对顶替界面的影响规律。

      (2) 在套管居中度较差时, 由于窄边滞留严重, 通过调整水泥浆流变性无法改善水泥浆顶替效率; 在居中度较高条件下, 通过降低n值改善水泥浆流动性, 有利于提高水泥浆顶替效率。

      (3) 结合对不同居中度条件下开展不同流变性水泥浆模拟分析可以发现, 提高顶替效率的关键是保证套管居中度, 同时通过调整水泥浆流变性, 合理控制水泥浆n值, 以提高顶替效率。

      (4) 本文主要是针对居中度和水泥浆流变性开展顶替效率的分析, 后续将进一步结合实际固井施工参数和浆体性能, 综合开展替浆排量、固井流体密度差等因素开展耦合分析。

参考文献 (18)

目录

    /

    返回文章
    返回