含杂质CO2 管道输送水合物形成规律研究

2021-04-21陈兵张曙旋郭焕焕

石油与天然气化工 2021年2期

陈兵张曙旋郭焕焕

西安石油大学机械工程学院

为了减轻全球温室效应,我国目前高度重视CCUS技术的相关研究工作,多条CO2长输管线开始设计筹建。但长输管线工况复杂,在高温低压条件下,气体在输送过程中很容易产生水合物,引起如冲击、阻塞、液位波动、管道堵塞甚至爆管等问题[1-6]。CO2水合物一旦形成,轻则影响流体流量,重则将管线堵塞甚至报废。管道中CO2水合物堵塞由两个过程引起:①达到水合物平衡的压力、温度,水合物开始成核、生长;②依靠流体颗粒之间的黏附力致使水合物聚集[7-12]。因此,在长输管线中防治水合物的生成是管道安全输送的关键,研究CO2水合物在管道输送过程中的形成规律是防治水合物形成的基础[13-14]。

研究CO2水合物的形成规律对防治和利用水合物均有重要意义,目前,国内外关于CO2水合物的生成多是在静态条件下,或研究在高压反应釜内搅拌条件下生成水合物,对管道内流动条件下水合物生长的研究相对较少。Davies等[15]开发了一套多相流模拟软件可以预测多相流管线中的水合物堵塞,能够预测多相流管路中水合物形成速率等多种流动参数,但是该软件不够成熟;Turner等[16]提出了一种水合物生成动力模型,他认为水合物形成和生长要经历成簇反应、临界尺寸(水合物晶核)的形成以及晶体的生长3个阶段;唐翠萍团队在长度50 m 的高压循环线路里,对CO2水合物形成过程和流动特性进行了研究,为后续水合物形成规律及堵塞提供理论和实验基础[17-19]。本研究根据PVTSIM 软件生成的含杂质CO2水合物的相平衡曲线,利用OLGA 软件建立水平、弯管管输模型对CO2的形成规律进行研究,从而加快CO2长输管线安全性研究进展,对实际工况的安全运行具有重要的现实意义。

1 适用于含杂质CO2 水合物的状态方程研究

掌握含杂质CO2的相态特征和物理特性是精确计算质量平衡、能量平衡以及传质传热和压力损耗的基础,含杂质CO2的输送模型建立在用状态方程来估算热力学性能与物理性能的基础上。研究CO2水合物在管道输送的相平衡条件,必须要了解整个CO2流的相态特征、压力、温度、黏度等其他相关热力学参数,再选择合适的状态方程,然后用PVTSIM 生成相应的CO2水合物相平衡曲线。

目前,关于用于含杂质CO2的状态方程还没有达到共识,常被用来预测含杂质CO2相态图和物理性质的状态方程有:Peng-Robinson(PR)方程、Soave-Redlich-Kwong(SRK)方程、Benedict-Webb-Rubin-Starling(BWRS)方程。用于模拟计算含杂质CO2物理性质的状态方程都是在综合相关文献的基础上,对可用实验数据进行对比评估而获得,其相关结论如下[20-24]:在7.0 MPa＜p＜15.0 MPa、0.0℃＜t＜12.0℃(即273.15 K＜T＜285.15 K)条件下,RK 方程对于CO2水合物的相平衡条件准确性比其他状态方程高,RK 方程见式(1)、式(2):

式中:p为气体压强,MPa;V为气体的摩尔体积,m3/mol;T为温度,K;R为气体常数,R=8.3143 J/(mol·K);a为常数,用于修正分子间引力;b为常数,用于修正体积;Tc为临界温度参数,K;pc为临界压力参数,MPa;w a=0.42748,w b=0.08664,与物质种类无关。

2 含杂质CO2 物理性质的研究

2.1 某油田捕集的CO2 组分与相平衡研究

本研究选用某油田一期36×104t/a 管道输送路线图,榆能化首站经过增压后的通过水平管道L11、L12输送至杏子川末站地区,输送距离总长105000 m。该油田的气源组成见表1所列。

表1 气源组成 %

选用表1某油田气源组成,利用PVTSIM 软件生成CO2水合物的相平衡曲线模拟值,如图1所示。研究认为,达到水合物生成的相平衡条件即有水合物生成[25-27]。图1中由水合物相平衡曲线将整体分为两个区,分别是水合物生成区和无水合物区。在水合物生成区,管道中有CO2水合物生成,且从图1可以看出:在高压低温条件下易形成水合物,即进入水合物生成区;当温度高于某临界值时,不再有水合物生成,即进入无水合物区。

2.2 含杂质CO2 流体的黏度特征研究

由于CO2流体黏度的大小决定了流体流动的传递性能,进而会影响水合物的生成速率,所以选用HYSYS软件对表1中含杂质CO2流体的黏度进行了数值模拟计算,其黏度随温度和压力的变化曲线如图2所示。

由图2可知,含杂质CO2流体的黏度的整体变化规律为:在同一压力下,随着温度的增加,CO2流体的黏度降低,在某一处出现大的突变,黏度逐渐平缓减小,这是因为随着温度的升高,CO2流体发生了相态的变化;在同一温度下,随着压力的增加,CO2流体的黏度增加。因此,CO2流体的黏度与温度和压力存在线性关系。

3 CO2 管道输送工艺水合物形成的模拟与分析

3.1 水平管道输送水合物形成的模拟与分析

用OLGA 建立水平管道模型,管道分L11、L12两段,管道的入口温度为45.0 ℃,入口压力为13.0 MPa,出口温度为6.6℃,出口压力为8.4 MPa。管道输送基本参数见表2。

表2 管道输送基本参数

参照图1 水合物生成模拟值,在距离入口47000 m 处已开始形成水合物。因此,将L11、L12总管道分为无水合物区和水合物生成区两个区块。如图3所示,在无水合物区0～20000 m 管段,温降比较明显。在管道输送过程中,L11段管道的入口压力、温度均比较高,且管道温度与土壤环境温度差距较大,因此,管道与土壤迅速换热,管道内部温度迅速降温。在无水合物区0～20000 m 管段,压降比较小,当流体刚进入管道L11时,流速较低,因此流阻损失小,压降较缓慢。

当CO2流体流经水合物生成区时,水合物开始形成。管道与土壤换热速率变缓,温降较缓慢,达到一定温度值后不再与环境进行换热;随着流速的增大,流阻损失也越大。因此,压降较无水合物生成区明显增大。

图4为不同出口压力下CO2流体黏度的变化曲线。由图4可知:在0～20000 m 管段内,由于温降非常明显,因此流体黏度骤增,黏度比接近3 倍;在20000～47000 m 管段内,流体的黏度增幅也在不断增加,但黏度较0～20000 m 管段变小;参考水合物生成模拟曲线,在47000 m 处已开始形成水合物;在47000 m 以后,随着CO2水合物的不断生成,流体黏度依旧在不断增加;在60000～105000 m 管段,参考图3压力温度变化曲线,虽然温度变化率降低了,但压力变化率增大,所以CO2流体黏度几乎不发生变化,这就导致水合物的生长驱动力降低,进而水合物停止生长。

根据不同的管道出口压力,流体黏度曲线的趋势是一致的。但随着管道出口压力从8.5 MPa增加到11.0 MPa,同一管道位置处,出口压力越大,形成CO2流体的黏度越大,更容易造成管道堵塞。而在出口压力较小的工况下,CO2流体黏度相对较小,对管道堵塞的可能性也较小。因此,在管道输送设计过程中,需要注意对管道出口压力的设计。

3.2 弯管管道输送水合物形成的模拟与分析

本研究选用了实际工况中比较常见的一段不利管路,即90°的垂直弯管建模,如图5所示:弯管总共分L13、L14两段,其中L13为水平管道,L14为斜管段。流体的入口温度为15.0 ℃,出口温度为8.5 ℃,出口压力为9.5 MPa。其管道输送基本参数见表3。

表3 管道输送基本参数

参考图1水合物生成模拟值,在距离入口9500 m处已开始形成水合物。因此,将L13、L14总弯管分为无水合物区和水合物生成区两个区块。由图6可知,在无水合物区,管道的入口压力较高、温度较低。因此,管道与土壤换热效率较低,管道内部温度降温较缓慢。且流体刚流进管道时,流速较低,流阻损失小,所以压降比较小。

当流体快要流经弯管段时,温度和压力达到了水合物生成区条件。此时,要防止生成的水合物在弯管处造成堵塞,影响流体在管道的流动性。流体流经弯道进入斜管段后,因为斜管与水平有一定的高度差,所以压降明显增加;但管道与土壤温差较小,因此温降依旧较缓慢。

如图7所示,在0～9500 m 管段流体的黏度在不断增加,在9500 m 处靠近弯管处已进入水合物生成区,流体的黏度降低后又升高,水合物浆在该处聚集有部分沉降。过了弯管处,水合物浆黏性依旧在不断增加,说明水合物在弯管处并未因沉降而堵塞管道,且过了弯管段后水合物在不断生长。在17500 m 管段后,流体黏性变化幅度不大,说明水合物在相对水平的管道中,晶核不易积聚,水合物生长缓慢。而弯管段易积聚,造成水合物堵塞的几率大。

由图7可见,随着管道出口压力从8.5 MPa增加到11.0 MPa,CO2流体黏度变化趋势不变,但同一位置处水合物的黏度在不断增加。压力驱动增加CO2流体的黏度,有助于水合物晶核的长大,从而加速水合物生成过程。因此,控制出口压力的大小,是控制水合物堵塞管道的有效手段。

3.3 CO2 管道输送应用实证

由延长油田产生的CO2经数千米的管道输送到气站,在输送过程中,现场管线时常会发生冰堵现象。相比于直管段,水合物多集中在弯管处,轻则影响气体流量,重则将管线堵死。尤其在寒冷的冬季,水合物堵塞现象几乎每天都在发生,严重制约了生产的正常进行。结合本研究内容,进一步验证了含杂质CO2水合物的形成规律,为今后采取相应预防CO2水合物生成的措施提供一定的依据。