超临界CO2水平环空携砂试验研究

2022-06-24王海柱李英杰

石油钻探技术 2022年3期

孙晓，王海柱，李英杰，郑永，陆群

（1.陕西延长石油（集团）有限责任公司研究院，陕西西安 710075；2.油气资源与探测国家重点实验室（中国石油大学（北京）），北京 102249；3.中国石油集团海洋工程有限公司，北京 100176）

超临界CO2具有类似气体的低黏度和接近液体的高密度，同时具有高扩散系数和良好的溶解性，因此，在非常规油气开发中具有巨大的优势[1-11]，有望形成集钻井、压裂、采气和CO2埋存为一体的高效油气开采技术。采用超临界CO2作为钻井流体钻井具有破岩门限压力低、机械钻速快的特点，能够缩短钻井周期[12-16]。此外，超临界CO2用于无水压裂时易形成复杂裂缝网络，并且能够在提高油气采收率的同时实现 CO2的地下埋存[2,17-18]。

超临界CO2水平环空携砂是超临界CO2钻井及压裂技术应用中的基础问题，涉及井眼清洁及压裂支撑剂泵注的水力参数优化设计[19-20]。与常规工作流体相比，超临界CO2物性参数对温度和压力敏感，工程参数对其在井筒内流动和传热影响的规律复杂。沈忠厚等人[20]通过模拟发现，超临界CO2水平环空携岩存在一个临界密度，可根据最小动能准则来确定不同密度超临界CO2携岩的最小流动速度。李良川等人[21]通过试验研究了超临界CO2在井筒内的携岩能力，发现试验数据与Angel理论模型的计算结果基本吻合。霍洪俊等人[22-23]通过数值模拟发现超临界CO2携带岩屑在水平环空中运移时会出现分层流动的现象，并分析了超临界CO2在斜井临界环空上返速度的变化规律。宋维强等人[24-25]模拟了超临界CO2在直井段和水平段环空中携带岩屑的运移规律，发现温度和压力对岩屑粒径的分布有影响。

由于进行超临界CO2携砂试验时压力、温度比较难控制，因此当前多采用模拟手段研究水平段环空中超临界CO2的携砂性能，并且已有试验研究由于试验装置尺寸和性能参数的限制，试验参数选取范围较窄，仅对超临界CO2的井筒携砂性能进行了初步探索，关键参数的影响规律尚不明确。基于此，笔者运用相似原理，设计了超临界CO2水平环空携砂试验装置，并以环空中砂粒运移速度和砂床高度为标准，利用该试验装置研究了注入质量流量、砂比、出口压力和流体温度对超临界CO2水平环空携砂能力的影响。研究结果可为超临界CO2钻井和压裂过程中水力参数的设计提供理论基础。

1 超临界 CO2 水平环空携砂试验装置

试验过程涉及CO2存储、超临界CO2生成及泵注、水平环空井筒、安全控制及循环管路等，其中水平环空井筒是携砂试验的关键。为了模拟现场水平环空流体携砂运移过程，笔者结合现场实际和室内试验条件，运用相似原理，设计了水平环空携砂模拟试验装置，其主要由水平环空井筒、加砂装置和除砂除水装置组成，如图1所示。

图1 水平环空携砂模拟试验装置Fig.1 Simulation device for sand-carrying tests in horizontal annulus

1）水平环空井筒。井筒外管为1Cr18Ni9Ti薄壁管，标准承压 30 MPa，总长 3.50 m，环空内径 38.0 mm、外径62.0 mm；内管主要为有机玻璃管，为防止砂粒碰撞摩擦而使有机玻璃管表面变粗糙导致透光性变差，正对可视窗口处的内管使用有机耐磨玻璃管。井筒上分布有压力传感器、温度传感器，并设有直径20.0 mm的对开耐高压石英玻璃观察窗，距井筒出口0.32 m。试验中可通过数据采集系统自动采集压力、流体温度等数据，并利用高速摄影装置监测记录砂粒在环空中的运移过程、流体相态变化、砂粒运动速度和砂床高度等相关试验数据。

水平井筒的设计满足几何相似和动力相似：模拟环空井筒的尺寸与现场油套环空尺寸满足几何相似；动力相似选择雷诺数，用于研究封闭系统内流体的运动，并考虑了流体的质量惯性力和黏性力。

a）几何相似。选择现场油套环空尺寸作为原型，套管外径dco为 139.7 mm、内径dci为 121.5 mm，油管外径dto为 74.5 mm、内径ddi为 62.0 mm，水力直径dcf为 46.9 mm；模拟井筒外管内径doi为 62.0 mm，内管外径dio为 38.0 mm，水力直径dmf为 24.0 mm，满足几何相似条件：

b）雷诺数相似。雷诺数相似的关系式为：

式中：ρf为实际流体密度，kg/m3；vf为实际环空流体流速，m/s；lf为实际流场特征长度，m；μf为实际流体黏度，Pa·s；ρm为模拟流体密度，kg/m3；vm为模拟环空流体流速，m/s；lm为模拟流场特征长度，m；μm为模拟流体黏度，Pa·s。

假设现场施工和室内模拟试验用超临界CO2的密度和黏度相同，则式（2）可变为：

由于试验所用泵排量的限制，试验过程中注入超临界CO2的质量流量较小，相当于模拟现场低排量下超临界CO2携带支撑剂运移的情况。

2）加砂装置。加砂装置由加砂罐、加砂套筒、加砂电机和搅笼等组成。搅笼与加砂电机连接，并置于加砂套筒内部，可通过加砂电机控制加砂速度，实现对砂比的精确控制，最高加砂速度350 L/h。试验中的砂比为加砂排量与携砂液的总排量之比。

3）除砂除水装置。除砂除水装置主要由砂粒收集装置和净化装置组成，砂粒收集装置用来分离CO2和砂粒，收集的砂粒可取出称重；净化装置主要用来对除砂后的CO2进行干燥净化，以满足循环使用的要求。

此外，试验所用泵注系统的最大排量为30 L/min，额定压力为 100 MPa，限定安全压力为 50 MPa。超临界CO2生成系统采用循环水浴加热，水温可在0～100 ℃ 间任意调节，误差±0.5 ℃。

2 试验方案及试验步骤

2.1 试验材料及方案

试验选用纯度为 99.9%的瓶装CO2，使用的模拟砂为20/40目的陶粒，平均粒径为0.6 mm，使用前进行洗净干燥，以避免杂质对试验结果和循环系统造成影响。

水平环空携砂试验主要研究不同泵注参数对砂粒在水平环空中的运移速度及砂床高度的影响。具体方案为：1）保持井筒流体温度为50 ℃，井筒流体出口压力为12 MPa，砂比为9%，观察超临界CO2质量流量分别为 0.135，0.160，0.178，0.223 和 0.268 kg/s时水平环空中砂粒的运移情况，以分析超临界CO2质量流量对砂粒运移的影响；2）保持超临界CO2泵注质量流量为 0.26 kg/s，井筒流体出口压力为 12 MPa，井筒流体温度为50 ℃，观察砂比分别为3.5%、4.5%、5.5%、6.5%和7.5%时水平环空中砂粒的运移情况，以分析砂比对砂粒运移的影响；3）保持超临界CO2泵注质量流量为0.26 kg/s，砂比为9%，井筒流体温度为 50 ℃，观察出口压力分别为 9，10，12，14 和 16 MPa时水平环空中砂粒的运移情况，以分析出口压力对砂粒运移的影响；4）保持井筒流体出口压力为12 MPa，超临界 CO2泵注质量流量为 0.26 kg/s，砂比为 9 %，观察井筒流体温度分别为40，45，50，55和 60 ℃ 时水平环空中砂粒的运移情况，以分析井筒流体温度对砂粒运移的影响。

2.2 试验步骤

1）打开气源，利用制冷机组将CO2液化存储在液态CO2储罐中，启动柱塞泵，将液态CO2泵至缓冲罐中进行加热，生成超临界CO2；

2）将超临界CO2缓慢注入井筒，同时打开井筒循环水浴保温装置，将柱塞泵排量调至设计流量；

3）关闭气源，开启井筒出口与储液罐间的节流阀，使系统内部流体达到循环状态，通过调节井筒出口回压控制井筒内部流体压力；

4）井筒内流体循环达到稳定状态后，开启加砂电机，通过调节电机转速调整加砂速度；

5）监测记录水平环空内颗粒的运动；

6）针对不同的试验参数重复以上各步骤。

3 砂粒运移影响因素分析

3.1 超临界CO2质量流量

流速是影响环空内流体携岩输砂效果的主要参数，不同超临界CO2质量流量下水平环空内砂粒的运移速度如图2所示。从图2可以看出，随着超临界CO2注入质量流量增大，砂粒在环空内的运移速度明显提高。此外，随着砂粒不断注入，不同质量流量下砂粒的运移速度均表现为逐渐增大并趋于稳定的趋势。这是因为环空中的砂粒以悬浮、跃移的形式输送，随着砂粒持续注入，砂粒会逐渐沉降在底部形成砂床（见图3），导致环空过流面积逐渐减小，使得超临界CO2流速不断增大，直至达到砂粒能够克服沉降阻力的临界流速并保持稳定。

图2 超临界CO2质量流量对水平环空中砂粒运移速度的影响Fig.2 Effect of supercritical CO2 mass flow on sand transport velocity in horizontal annulus

图3 水平环空中超临界CO2输送砂粒形态Fig.3 Shape of sand transported by supercritical CO2 in horizontal annulus

不同超临界CO2质量流量下水平环空内砂床的平衡高度如图4所示。从图4可以看出，砂床平衡高度随超临界CO2质量流量增大整体上呈线性减小趋势。由于超临界CO2在环空上部的流动速度比下部快，且高速区的范围会随超临界CO2质量流量增大而扩大[26]，导致环空内紊流区的范围和强度随超临界CO2质量流量增大而增大，作用在砂粒上的曳力和升力随之增大，使更多的砂粒进入悬浮层，并且运移距离也随之增大，最终使砂床平衡高度降低。因此，施工过程中在设备能力允许的情况下，应尽量提高超临界CO2的注入质量流量，以在环空内产生范围较大的能悬浮运移砂粒的高速流动区，避免井下砂埋和砂卡等的发生。

图4 不同超临界CO2质量流量下水平环空中砂床的平衡高度Fig.4 Equilibrium height of sand bed in horizontal annulus with different mass flow of supercritical CO2

3.2 砂比

由上文分析可知，砂粒运移速度趋于稳定时形成了高度稳定的砂床。而不同砂比下水平环空中砂粒的运移速度如图5所示。

图5 砂比对水平环空砂粒粒运移速度的影响Fig.5 Effect of sand concentration on sand transport velocity in horizontal annulus

由图5可知：环空内形成高度稳定的砂床前，砂比对砂粒在环空中运移速度的影响较弱；环空内形成高度稳定的砂床后，随着砂比增大，砂粒在环空中的运移速度加快，这说明砂比增大，形成稳定砂床的高度增大，使得过流面积减小，导致砂粒运移速度随砂比增大而加快。

不同砂比条件下，水平环空砂床的平衡高度如图6所示。

图6 不同砂比下水平环空砂床的平衡高度Fig.6 Equilibrium height of sand bed in horizontal annulus at different sand concentrations

由图6可知，砂床的平衡高度随砂比增大呈线性增大。这是因为环空内砂比增大，会使砂粒与砂粒、砂粒与管壁之间的相互作用增大，导致砂粒运移需要更大的驱动力；同时，砂比增大，超临界CO2和砂粒之间的动量交换相对减少[24]，砂粒更容易沉积在环空底部，最终导致砂床平衡高度增大。

3.3 出口压力

超临界CO2具有较强的压缩性，其物性参数对温度和压力的变化比较敏感。不同出口压力下水平环空中砂粒运移速度的变化如图7所示。

图7 出口压力对水平环空砂粒运移速度的影响Fig.7 Effect of outlet pressure on sand transport velocity in horizontal annulus

由图7可知，砂粒运移速度随出口升高而降低，但出口压力升至10 MPa后，砂粒运移速度的降低幅度逐渐减小。这是因为出口压力升高，导致环空内超临界CO2的密度增大，而超临界CO2的质量流量不变，因此根据质量守恒定律，可知其流速减小，使得砂粒的运移速度降低；另外，在高压条件下超临界CO2的密度和黏度对压力变化的敏感性较低，使压力升高对砂粒运移速度的影响减弱。

不同出口压力下砂床的平衡高度如图8所示。由图8可知，砂床平衡高度整体上呈现出随出口压力升高而升高的趋势，这主要与不同出口压力下环空内流体的流速相关。现场应用中，应合理调节环空内压力，以使超临界CO2达到所需要的密度，同时优化设计其他参数，尽量降低压力升高对环空内流体流速的影响，避免砂粒过多堆积。

图8 不同出口压力下水平环空砂床的平衡高度Fig.8 Equilibrium height of sand bed in horizontal annulus under different outlet pressure

3.4 流体温度

流体温度是利用超临界CO2开发非常规油气资源的关键参数，涉及相态调控等多个方面。不同流体温度下水平环空中砂粒运移速度的变化情况如图9所示。

图9 流体温度对水平环空砂粒运移速度的影响Fig.9 Effect of fluid temperature on sand transport velocity in horizontal annulus

由图9可知，在相同超临界CO2质量流量下，砂粒的运移速度随着流体温度升高而增大。这是因为随着流体温度升高，超临界CO2密度降低，在质量流量相同的情况下会使环空流速增大，砂粒获得了更大的驱动力，进而使其运移速度增大。利用王海柱编制的超临界CO2物性参数计算软件[26]，计算得到温度 40，45，50，55 和 60 ℃ 下井筒内超临界CO2流体的密度分别为 717.49，657.25，583.99，503.47和 433.15 kg/m3，即随着温度升高，环空中流体的密度分别降低了60.24，73.26，80.52和70.32 kg/m3。由此可知，水平环空中超临界CO2的密度在温度由50 ℃升至55 ℃时降低幅度最大，达80.52 kg/m3，导致砂粒的运移速度在该温度区间上有较大的增大幅度。

不同温度下水平环空砂床的平衡高度如图10所示。