The Prediction of Hydrate Formation Regions in the Wellbore during Supercritical Carbon Dioxide Drilling
-
摘要: 为了保证超临界CO2钻井安全并快速钻进,需要解决井筒内水合物生成的问题。为此,在分析水合物形成机理的基础上,建立了超临界CO2钻井井筒水合物形成区域预测模型,并给出了模型的定解条件和数值求解方法。通过设计算例进行了计算分析,结果表明:环空内水合物形成区域的临界井深,随注入温度或水合物抑制剂加量增大呈二次多项式非线性下降;随井口回压增大先呈对数函数增大、后呈二次多项式关系增大,增大的幅度逐渐变小。分析结果可为超临界CO2钻井防治水合物形成提供理论参考。Abstract: In order to improve the safety and efficiency of SC-CO2 drilling (supercritical carbon dioxide drilling), it is necessary to deal with the formation of CO2 hydrate in the wellbores. Based on the hydrate formation mechanism analysis, a model for predicting the formation region of CO2 hydrate in the wellbore was built up, and the definite conditions and numerical resolution methods were proposed. Computational analysis was carried out by designing examples. It was shown that hydrate formation in the annulus decreased non-linearly in quadratic polynomial relations with the increasing of injection temperature or inhibitor dosage. And with the increasing of wellhead back pressure, hydrate formation increased with increasing rate reducing gradually (first in logarithmic relations and then in quadratic polynomial relations). The research achievements could provide a theoretical reference for hydrate prevention during SC-CO2 drilling.
-
油田长期注水开发导致地层的非均质性加剧,注入水沿高渗透层突入油井,导致油井含水率升高、产油量降低,增产稳产面临的难度很大。注水井深部调剖是油田稳产增产的重要技术措施之一,但对调剖剂的要求较高。冻胶型堵剂具有成胶时间可调、成胶强度高和价格便宜的特点,在油田调剖作业中得到了广泛的应用,但冻胶在注入地层的过程中受机械剪切、色谱分离和地层水稀释等多种因素的影响,其成冻时间、形成冻胶的强度和进入地层的深度难以控制,导致冻胶在地下的成胶效果变差,影响了调剖效果和有效期[1-7]。为解决以上问题,将水膨体发展为地面成胶体系,解决了地下成胶效果不可控的问题,但其初始粒径较大、且制备工艺复杂,影响了其在海上油田调剖调驱中的应用[8-12]。近年来发展起来聚合物微球调驱技术很好地解决了上述问题,但其采用乳液聚合方式进行制备,要求引发时间、聚合时间及聚合温度精准,对制备设备要求高,制备工艺相对复杂,不能在线生成注入,且合成原料中包含表面活性剂,增大了制备成本[13-16]。
针对目前冻胶、水膨体和聚合物微球调剖存在的问题,笔者采用机械剪切法制备了一种新的海上油田剖面调整用分散共聚物颗粒体系,不仅具有冻胶的特点,还不受地面剪切、稀释和色谱分离等因素的影响,能够变形进入地层深部,封堵地层深部的大孔道,调整渗流剖面,实现深部液流转向。同时,对制备工艺进行了探索和优化,为分散共聚物颗粒体系的现场应用奠定了基础。
1. 分散共聚物颗粒体系的制备
根据共聚物黏度的变化情况,制备过程分为共聚物制备和分散共聚物颗粒制备2个阶段。
1.1 共聚物的制备
化学材料:丙烯酰胺(AM),工业品;丙烯酸,工业品;N, N亚甲基双丙烯酰胺(MBA),分析纯;AMPS,分析纯;其余试剂均为分析纯。
以丙烯酰胺为主要原料,按比例加入其他添加剂,在水中溶解后置于65 ℃恒温水浴中,聚合物进行交联反应,形成黏度极高的高分子聚合物。其交联反应式为:
共聚物形成阶段分为引发阶段、快速交联阶段和稳定阶段。共聚物成胶后,形成三维凝胶网络结构,黏度不再增加。单体AM质量分数选用3%~6%,AM与MBA的质量比为250∶1~100∶1。不同AM/MBA质量比的共聚物样品合成时间及合成后的黏度如表1所示。
表 1 共聚物的制备参数Table 1. Preparation parameters of the copolymerAM质量分数,
%AM/MBA
质量比合成时间/
h共聚物黏度/
(mPa·s)5 167∶1 2.0 16 460 5 250∶1 2.5 15 320 1.2 分散共聚物颗粒的制备
分散共聚物颗粒形成阶段分为破碎阶段、研磨阶段和稳定阶段,研磨转速为1 000 r/min,研磨时间为3~15 min,制得不同粒径分布的均一分散共聚物颗粒水相溶液。为了更好地了解影响分散共聚物颗粒制备的因素,测试了制备过程中体系黏度、粒径分布的变化情况。
1.2.1 制备过程中黏度的变化
分散共聚物颗粒制备过程中,体系黏度的变化是表征其性能的主要指标。研磨时间为1,5,10和15 min时的黏度分别为9.8,3.9,3.8和3.7 mPa·s。由此可以看出,将共聚物溶液(AM质量分数为5%,AM/MBA质量比为250∶1)加入胶体磨中高速研磨,其黏度在5 min后就迅速降至5.0 mPa·s以下,得到分散共聚物颗粒。
1.2.2 粒径分布
使用马尔文3000激光粒度仪测试分散共聚物颗粒溶液的粒径分布情况,测量前使用蒸馏水将分散共聚物颗粒溶液稀释至400 mg/L,每个待测样品测3个平行样,得到研磨时间对体系粒径分布的影响结果,如图1和表2所示。
表 2 分散共聚物颗粒粒径随研磨时间的变化Table 2. Viscosity change of dispersed copolymer particles in grinding process剪切时间/min 平均粒径/μm 最大粒径/μm 最小粒径/μm 1 1 520.0 3 080.0 31.1 5 135.0 352.0 4.6 从图1可以看出,研磨1 min时,体系的粒径分布范围较广,平均粒径为1 520 μm;持续研磨5 min后,体系的粒径分布变窄,平均粒径降至135 μm,说明分散共聚物颗粒制备过程中,随着研磨时间增长,体系粒径分布更为均一。另外,从表2可以看出,最大粒径与最小粒径的比值随着研磨时间增长而降低,进一步验证了以上结果。
2. 分散共聚物颗粒溶液性能评价
2.1 流变性测试
采用双狭缝模型,剪切速率变化范围为1~1 000 s–1,测试质量分数为3%的分散共聚物颗粒溶液的黏度随剪切速率的变化情况。试验采用模拟垦利油田注入水,总矿化度4 485.0 mg/L,钙镁离子含量为71.0 mg/L,pH值为7.34。使用Anton Paar旋转流变仪,测定了65 ℃下分散共聚物颗粒溶液的流变性能,结果见图2。由图2可知,在低剪切速率下,分散共聚物颗粒溶液的黏度随剪切速率升高而降低;在中剪切速率下,其溶液黏度基本不受剪切速率的影响,呈现出牛顿流体的性质;在高剪切流动状态下,其溶液黏度随剪切速率升高而略有升高。分析认为,在高剪切速率下,分散共聚物不断重新排布、相互碰撞而形成网状结构,导致黏度略有升高。剪切速率为1~1 000 s–1时,3%分散共聚物颗粒溶液的黏度低于6.0 mPa·s;剪切速率为7.34 s–1时,其溶液黏度为1.1 mPa·s,略高于水。
2.2 分散稳定性测试
使用Turbiscan多重光稳定性分析仪,将质量分数为3%的分散共聚物颗粒溶液放入分析仪中,每间隔15 min,对测试样品从底部到顶部扫描一次,测试其背散射光强和透射光强,并据此直接计算得到测试样品的稳定性动力学指数。稳定性动力学指数累计了测试样品所有光强的变化,反映了测试样品的稳定程度。稳定性动力学指数越大,测试样品越不稳定;稳定性动力学指数小于3.0时,测试样品的稳定性较好。采用高级分析模块,计算分散共聚物颗粒溶液的稳定性动力学指数,得到其稳定性动力学指数随时间的变化,如图3所示。
从图3可以看出,分散共聚物颗粒溶液的稳定性动力学指数小于3.0,性能较为稳定,分散共聚物颗粒溶液从井口注入后运移至井筒及近井地带期间能够保持良好的分散稳定性能,不会发生沉降聚集。
3. 封堵性和运移性评价
3.1 封堵性评价试验
采用填砂管驱替装置,填砂管长度50.0 cm,渗透率5 000 mD左右,测试相同用量条件下分散共聚物颗粒体系和聚合物凝胶体系在多孔介质内的封堵情况。在模拟地层温度下,以2.0 mL/min的注入速度,分别注入1.0倍孔隙体积的分散共聚物颗粒体系和聚合物凝胶体系,记录注入压力的变化情况。后续水驱至压力平稳,记录后续水驱阶段压力的响应情况,结果如图4所示。
从图4可以看出,在用量相同条件下,注入1.0倍孔隙体积的1 500 mg/L聚合物凝胶体系时,注入压力升至0.1 MPa;注入1.0倍孔隙体积的分散共聚物颗粒体系时,注入压力升至0.5 MPa。后续水驱阶段,2种体系均显示出较高的残余阻力系数,其中,聚合物凝胶体系的注入压力升至0.9 MPa,分散共聚物颗粒体系的注入压力升至2.8 MPa。
因此,渗透率大于5 000 mD的高渗透储层进行剖面调整时,在调剖剂用量相同的条件下,分散共聚物颗粒体系具有更低的注入压力,即具有更优的注入性能和更强的封堵能力。
3.2 模拟地层条件的长填砂管运移性试验
填砂模型封堵性试验结果表明,分散共聚合颗粒体系具有良好的封堵性能,为了更真实地反映其在地层条件下的封堵性,开展了长填砂管条件下的封堵运移性试验。采用10.00 m填砂模型进行长距离运移性试验,模型渗透率为6 009 mD,沿填砂管均布5个测压点,评价其深部运移及封堵性能。
试验时首先以1.5 mL/min的注入速度向填砂模型内注入地层水,至模型内部压力平稳;然后再以1.5 mL/min的注入速度向填砂模型内注入1.0倍孔隙体积的分散共聚物颗粒体系;关闭注入端和采出端,在65 ℃恒温箱中放置10 d,再以1.5 mL/min的注入速度向填砂模型内注入后续水,记录驱替过程中压力的变化情况,结果见图5。
从图5可以看出,开始注水后,各测压点的压力迅速上升;上升至一定压力后,压力逐渐降低;后续水突破后,距离注入端越远,压力下降越缓慢。这说明分散共聚合颗粒被注入水突破后,封堵体系在填砂模型内仍保持良好的封堵性能。因此,在模拟地层条件下,分散共聚物颗粒体系对于渗透率大于5 000 mD的高渗透储层具有良好的运移性和封堵性。
4. 结论与建议
1)通过采用特殊的交联技术和分散技术,形成的高黏聚合物经研磨后,制得纳米–微米级的均一分散水溶液。
2)通过调节主剂AM的质量分数和AM/MBA的质量比,可得到不同粒径分布的分散共聚物颗粒溶液;其初始黏度可控制在10 mPa·s以内,具有良好的注入性、深部运移性和对高渗透储层的有效封堵性能;制备的分散共聚物颗粒溶液经高速剪切后,黏度和粒度变化较小,表现出良好的抗剪切性能。
3)分散共聚物颗粒体系的黏度可控,工艺上可实现在线注入,可进入地层深部,对高渗透地层具有较好的封堵效果,可实现深部调剖。
-
[1] Al-Adwani F A,Langlinais J,Hughes R G.Modeling of an underbalanced drilling operation utilizing supercritical carbon dioxide[R].SPE 114050,2008.
[2] 霍洪俊,王瑞和,倪红坚,等.超临界二氧化碳在水平井钻井中的携岩规律研究[J].石油钻探技术,2014,42(2):12-17. Huo Hongjun,Wang Ruihe,Ni Hongjian,et al.Cuttings carrying pattern of supercritical carbon dioxide in horizontal wells[J].Petroleum Drilling Techniques,2014,42(2):12-17. [3] Kolle J J.Coiled-tubing drilling with supercritical carbon dioxide[R].SPE 65534,2000.
[4] 邱正松,谢彬强,王在明,等.超临界二氧化碳钻井流体关键技术研究[J].石油钻探技术,2012,40(2):1-7. Qiu Zhengsong,Xie Binqiang,Wang Zaiming,et al.Key technology on drilling fliud of supercritical carbon dioxide[J].Petroleum Drilling Techniques,2012,40(2):1-7. [5] van der Waals J H,Platteeuw J C.Clathrate solutions[J].Advances in Chemical Physics,1959,2:1-57.
[6] Parrlsh W R,Prausnltz J M.Dissociation pressures of gas hydrates formed by gas mixtures[J].IEC Chemical Process Design and Development,1972,11(1):26-33.
[7] Ng Heng-Joo,Robinson D B.Hydrate formation in systems containing methane,ethane,propane,carbon dioxide or hydrogen sulfide in the presence of methanol[J].Fluid Phase Equilibria,1985,21(1/2):145-155.
[8] Holder G D,Angert P F.Simulation of gas production from a reservoir containing both gas hydrates and free natural gas[R].SPE 11105,1982.
[9] Englezos P.Computation of the incipient equilibrium carbon dioxide hydrate formation conditions in aqueous electrolyte solutions[J].Industrial Engineering Chemistry Research,1992,31(9):2232-2237.
[10] Bishnoi P R,Dholabhai P D.Equilibrium conditions for hydrate formation for a ternary mixture of methane,propane and carbon dioxide,and a natural gas mixture in the presence of electrolytes and methanol[J].Fluid Phase Equilibria,1999,158/159/160:821-827.
[11] Nasrifar K,Moshfeghian M.Computation of equilibrium hydrate formation temperature for CO2 and hydrocarbon gases containing CO2 in the presence of an alcohol,electrolytes and their mixtures[J].Journal of Petroleum Science and Engineering,2000,26(1/2/3/4):143-150.
[12] 宫智武,张亮,程海清,等.海底天然气水合物分解对海洋钻井安全的影响[J].石油钻探技术,2015,43(4):19-24. Gong Zhiwu,Zhang Liang,Cheng Haiqing,et al.The influence of subsea natural gas hydrate dissociation on the safety of offshore drilling[J].Petroleum Drilling Techniques,2015,43(4):19-24. [13] 宋中华,张士诚,王腾飞,等.塔里木油田高压气井井下节流防治水合物技术[J].石油钻探技术,2014,42(2):91-96. Song Zhonghua,Zhang Shicheng,Wang Tengfei,et al.Downhole throttling technology for gas hydrate prevention in deep gas wells of Tarim Oilfield[J].Petroleum Drilling Techniques,2014,42(2):91-96. [14] 王海柱,沈忠厚,李根生.超临界CO2钻井井筒压力温度耦合计算[J].石油勘探与开发,2011,38(1):97-102. Wang Haizhu,Shen Zhonghou,Li Gensheng.Wellbore temperature and pressure coupling calculation of drilling with supercritical carbon dioxide[J].Petroleum Exploration and Development,2011,38(1):97-102. [15] Hasan A R,Kabir C S.Heat transfer during two-phase flow in wellbores:part I:formation temperature[R].SPE 22866,1991.
[16] Wang Weihua,Wen Hao,He Xianfeng,et al.Design and development of database system on the physical properties and phase equilibria of gas hydrates[J].Computers and Applied Chemistry,2005,22(6):35-40.
[17] Patwardhan V S,Kumar A.A unified approach for prediction of thermodynamic properties of aqueous mixed-electrolyte solutions:part I:vapor pressure and heat of vaporization[J].AIChE Journal,1986,32(9):1419-1428.
[18] 查丽,梁德青,何松.CO2水合物在电解质溶液中的相平衡预测[J].工程热物理学报,2012,33(7):1105-1108. Zha Li,Liang Deqing,He Song.Prediction of CO2 hydrate phase equilibria in aqueous electrolytes solutions[J].Journal of Engineering Thermophysics,2012,33(7):1105-1108. [19] Adisasmito S,Frank R J,Sloan E D Jr.Hydrates of carbon dioxide and methane mixtures[J].Journal of Chemical and Engineering Data,1991,36(1):68-71.
[20] Vlahakis J G,Chen H S,Suawandi M S,et al.The growth rate of ice crystals:the properties of carbon dioxide hydrate,a review of properties of 51 gas hydrates[M].US Department of the Interior,1972.
[21] Nasrifar K,Moshfeghian M.A model for prediction of gas hydrate formation conditions in aqueous solutions containing electrolytes and/or alcohol[J].The Journal of Chemical Thermodynamics,2001,33(9):999-1014.
[22] Dholabhai P D,Kalogerakis N,Bishnoi P R.Equilibrium conditions for carbon dioxide hydrate formation in aqueous electrolyte solutions[J].Journal of Chemical and Engineering Data,1993,38(4):650-654.
[23] Bozzo A T,Hsiao-Sheng C,Kass J R,et al.The properties of the hydrates of chlorine and carbon dioxide[J].Desalination,1975,16(3):303-320.
[24] Daugherty R L,Franzini J B,Finnemore E J,et al.Fluid mechanics with engineering applications[M].New York:McGraw-Hill,1985.
[25] Span R,Wagner W.A new equation of state for carbon dioxide covering the fluid region from the triple-point temperature to 1 100 K at pressures up to 800 MPa[J].Journal of Physical and Chemical Reference Data,1996,25(6):1509-1596.
[26] Vesovic V,Wakeham W A,Olchowy G A,et al.The transport properties of carbon dioxide[J].Journal of Physical and Chemical Reference Data,1990,19(3):763-808.
[27] Wang Zhiyuan,Sun Baojiang,Wang Jintang,et al.Experimental study on the friction coefficient of supercritical carbon dioxide in pipes[J].International Journal of Greenhouse Gas Control,2014,25(6):151-161.
[28] Bruch A,Bontemps A,Colasson S.Experimental investigation of heat transfer of supercritical carbon dioxide flowing in a cooled vertical tube[J].International Journal of Heat and Mass Transfer,2009,52(11/12):2589-2598.
[29] 何淼,柳贡慧,李军,等.多相流全瞬态温度压力场耦合模型求解及分析[J].石油钻探技术,2015,43(2):25-32. He Miao,Liu Gonghui,Li Jun,et al.Solution and analysis of fully transient temperature and pressure coupling model for multiphase flow[J].Petroleum Drilling Techniques,2015,43(2):25-32. -
期刊类型引用(1)
1. 贾志伟,程长坤,朱秀雨,濮兰天,韩宇,扈福堂. 青海油田尕斯N_1–N_2~1超高盐油藏复合驱提高采收率技术. 石油钻探技术. 2021(05): 81-87 . 本站查看
其他类型引用(0)
计量
- 文章访问数: 2225
- HTML全文浏览量: 70
- PDF下载量: 2785
- 被引次数: 1