张建昌 王立涛 周晓宇 张嘉安 张 妮 于 丹 丁 妮

1. 中国石油长庆油田分公司第三输油处, 宁夏 银川 750000;2. 西安石大派普特科技工程有限公司, 陕西 西安 710065

0 前言

含蜡原油输送过程中,当油品温度低于析蜡点温度时,蜡晶析出并在管壁沉积,导致管道内有效流通面积减小,管道输送能力下降,严重时甚至会引发蜡堵事故,造成重大经济损失[1-3]。

针对原油管道蜡沉积问题,国内外学者大多采用室内环道实验或软件仿真实验的方式开展研究[4-5]。室内环道实验通过控制变量法定量研究各变量对蜡沉积规律的影响,基于实验数据建立蜡沉积生成速率预测模型[6-7];软件仿真实验通过仿真软件建立管道动态模型,模拟蜡沉积厚度变化规律[8-10]。最后将室内环道实验或软件仿真实验的研究成果应用于现场,制定相应的清蜡防蜡方案。

然而,上述研究方法在现场管输应用中存在一定局限性,主要因为：1)生产现场管道规模庞大、工艺复杂,室内环道实验规模较小,难以还原现场真实工况;2)生产现场工况多变(增压、反输、倒罐、清管等),油品进出站压力、温度、流量均处于动态变化过程,且埋地温度受季节与天气影响而变化频繁;3)室内环道实验与软件仿真实验难以模拟现场复杂多变的工况,导致实验成果在生产现场的适用性较差。

针对上述问题,本文充分利用SCADA系统所采集的运行数据进行数据分析和挖掘,进而评价实际生产管道的蜡沉积规律和状态,并论证了该方案的可行性;即采用压差法计算原油管道当量蜡沉积厚度,分析油品输量、平均温度、管道埋地温度对当量蜡沉积厚度的影响,明确油品温度与油流冲刷对当量蜡沉积厚度的影响,并建立“蜡沉积动力学”预测模型,为生产现场评估管道蜡沉积情况提供新思路。

1 现场管输蜡沉积分析

1.1 管道概况

西北地区JF原油管道全长113.5 km,管道规格Φ377 mm×6 mm,任务输量250×104t/a,最小起输量117×104t/a,起终点高程差173 m。油品物性经过化验后,密度与黏度随温度变化情况见图1,平均温度下原油密度与黏度分别按式(1)与式(2)计算;原油析蜡点温度为34.7 ℃,凝点温度为21 ℃。

ρPJ=0.875 3-7.498×10-4TPJ

(1)

lgμPJ=-1.271-0.015 83TPJ

(2)

式中:ρPJ为油品平均温度下的密度,g/cm3;μPJ为油品平均温度下的黏度,Pa·s;TPJ为油品平均温度,℃。

图1 油品密温关系与黏温关系图Fig.1 Oil product density-temperature relationship and viscosity-temperature relationship

1.2 计算方法

准确测量蜡沉积厚度是开展蜡沉积研究的前提。室内环道实验中,管道蜡沉积厚度测量方法主要有直接测量法、传热测量法、压差测量法等。其中,直接测量法将蜡沉积物清理出管道并测量其质量或体积,进而计算出蜡沉积层厚度[11];传热测量法根据油品与环境之间热阻的变化计算出沉积层厚度[12-13];压差测量法通过监测管道两端压差变化计算出沉积层厚度[14-15]。

考虑到在生产现场蜡沉积物取样困难,且蜡沉积量及油品与环境换热量难以实时计量,直接测量法与传热测量法在现场应用不便,故本次研究选择压差测量法计算蜡沉积厚度,该方法仅需要油品进出站压力、温度、流量等数据,根据管道压差变化反算出管道有效内径,见式(3)。

(3)

式中:Δpf为蜡沉积管段压降,Pa;L为蜡沉积管段长度,m;d为管段有效内径,m;Q为体积流量,m3/s;ρ为流体密度,kg/m3;f为范宁摩阻系数。

f可由式(4)计算:

(4)

式中:c、n为常数,当流态为层流时,c=16,n=1;当流态为紊流时,c=0.046,n=0.2；μ为流体黏度，Pa·s。

根据管道摩擦压降间接计算管道有效内径后,可计算出管道当量蜡沉积厚度:

(5)

式中:y为管道当量蜡沉积厚度，m;d0为管道原始管内径,m。

考虑到管道沿程油品温度不断变化,原油出站温度与进站温度均不能概括管道整体油品温度情况,故选用油品平均温度计算油品密度与黏度。

(6)

式中:TPJ为油品沿程平均温度,℃;TR为油品起点温度,℃;TZ为油品终点温度,℃。

在原油输送过程中,仅当原油温度低于析蜡点温度时蜡晶才会析出,原油管道蜡沉积分布并不均匀[16];而随着管道沿线埋地温度变化以及油品输送工艺的调整,管道内沉积物厚度变化规律会更加复杂[17];因此式(5)计算结果并不是管道实际蜡沉积厚度,而是根据管路压降反算出的当量蜡沉积厚度,体现了局部沉积物厚度变化对管道整体水力学特性的影响。本文1.3节选用2020年5—6月管道9次清管数据验证计算结果准确性。

1.3 准确性验证

清管作业时,蜡沉积物被剥离并被油流裹挟流向下游,清管球前会形成一段蜡沉积物与油品混合的区间;收球时,大部分蜡油混合物会继续流向下游,少部分蜡油混合物跟随清管球进入收球筒。因此,可通过分析收球筒中蜡沉积物与原油占比判断清管效果。

取清管前后的站点压力、温度、流量数据计算当量蜡沉积厚度;记录每次清管后收球筒中蜡油混合物质量及蜡沉积物质量,计算结果见表1。

表1 9次清管数据与计算结果表

将收球筒中蜡沉积物质量与蜡油混合物质量之比定义为蜡沉积物质量分数,见式(7);蜡沉积物质量分数越大,表明蜡油混合物中蜡沉积物含量越多,即该次清管作业清出蜡沉积物越多。

(7)

式中:Mwax为收球筒中蜡沉积物质量,kg;Mwax-oil为收球筒中蜡油混合物质量,kg;η为蜡沉积物质量分数。清管前后当量蜡沉积厚度变化量与蜡沉积物质量分数关系见图2。

图2 当量蜡沉积厚度变化量与蜡沉积物质量分数关系图Fig.2 The relationship between the change of equivalent wax deposition thickness and the mass fraction of wax deposition

由图2可知,当量蜡沉积厚度变化量与蜡沉积物占蜡油混合物中质量分数呈线性正比关系;清管前后当量蜡沉积厚度变化量较大时,该次清管清出蜡沉积物较多,蜡沉积物占蜡油混合物中质量分数较大;当量蜡沉积厚度变化量较小时,该次清管清出蜡沉积物较少,蜡沉积物占蜡油混合物中质量分数较小。

当量蜡沉积厚度的变化规律与管道内蜡沉积物质量分数的变化规律一致，证明本文提出的当量蜡沉积厚度计算方法准确性较高。

2 蜡沉积厚度规律分析

影响蜡沉积的因素诸多,包括油品组成、流型、流态、油温、流速、油壁温差、沉积时间等[18-19];而生产现场工况多变,无法采用室内实验的方式(控制变量法)分析蜡沉积规律;本节基于JF原油管道2018—2019年712组生产数据,分析当量蜡沉积厚度变化规律。

考虑到影响蜡沉积的因素众多,而原油输送过程中油品组成、流型、流态变化情况较少,粗糙度难以准确测量,故本文主要研究油品流量、油品平均温度以及管道埋地温度3个参数对蜡沉积的影响情况。

2.1 油品流量

JF原油管道2018年1月—2019年12月油品流量在190～340 m3/h范围内变化频繁,无法如室内实验般保持油流稳定,油品流量与当量蜡沉积厚度关系见图3。

图3 油品流量与当量蜡沉积厚度关系图Fig.3 The relationship between pipeline equivalent wax deposition thickness and flow rate

由图3可知,当量蜡沉积厚度并未随时间一味增大,而是与管道流量呈负相关,当流量减小时,管道当量蜡沉积厚度增大;当流量增大时,管道当量蜡沉积厚度随之减小。这是因为流量变化会影响剪切剥离作用[20-21]，当流量增大时,由于剪切剥离作用增强,油流对蜡沉积层的冲刷作用随着增大,当量蜡沉积厚度随着减小。研究结果表明,在生产现场,剪切剥离作用对蜡沉积影响较大。

生产现场输量变化频繁,而室内环道实验与仿真实验时,难以模拟如此复杂的输量变化,通常会保持输量恒定,这必然导致蜡沉积厚度预测结果与实际相差较大。

2.2 油品平均温度

油品沿程温降变化受埋地温度、输量等因素影响,油品出站温度与进站温度均不能概括油品温度变化情况,因此分析油品平均温度与当量蜡沉积厚度的关系影响见图4。

图4 当量蜡沉积厚度与油品平均温度关系图Fig.4 The relationship between the equivalent wax deposition thickness of the pipeline and the average temperature of the oil

由图4可知,当量蜡沉积厚度与油品平均温度呈负相关,即当油品平均温度较高时,当量蜡沉积厚度处于较低水平,这是因为油品平均温度较高时,油品温度低于析蜡点管段的长度会缩短,蜡沉积对油品流动产生的影响降低;但在局部时段(如2019年4—6月),随着油品平均温度升高,当量蜡沉积厚度逐渐升高,这是因为该时段输量较低(图3),剪切剥离作用影响较小,油品平均温度增大会增大油壁温差,增强分子扩散作用的影响。

2.3 管道埋地温度

管道埋地温度与当地气候及季节变化相关,JF原油管道所处地区夏秋两季管道埋地温度较高,春冬两季埋地温度较低。管道埋地温度与当量蜡沉积厚度关系见图5。

图5 当量蜡沉积厚度与管道埋地温度关系图Fig.5 Relationship between pipeline equivalent wax deposition thickness and pipeline buried temperature

由图5可知,管道埋地温度与当量蜡沉积厚度呈负相关,当管道埋地温度增大时,当量蜡沉积厚度随之减小;管道埋地温度减小时,当量蜡沉积厚度随之增大。这是因为管道埋地温度变化会导致原油与管壁之间温度梯度改变,进而影响分子扩散作用[22];当管道埋地温度升高时,原油与管壁之间温度梯度降低,分子扩散作用减弱,蜡沉积速率降低;反之,当管道埋地温度降低时,原油与管壁之间温度梯度升高,分子扩散作用增强,蜡沉积速率随着升高。

由上述分析可知,分子扩散与剪切剥离作用对管道蜡沉积的影响均较大,而实际生产过程中,油品流量、油品平均温度与管道埋地温度均处于动态变化过程,分析蜡沉积厚度变化规律较困难,因此需要建立蜡沉积动力学预测模型,根据实时数据预测蜡沉积厚度变化情况。

3 蜡沉积动力学预测模型

考虑到原油管道蜡沉积同时受分子扩散作用与剪切剥离作用影响,本次研究选用Hsu模型[23-24]建立蜡沉积速率预测模型,见式(8)。

(8)

由于管道沿程温降较大,本文选用油品平均温度与管道埋地温度计算径向温度梯度,见式(9)。

(9)

式中:T0为管道埋地温度,℃;r为有效半径,m。

根据管道实测数据拟合蜡沉积速率公式,回归系数为k与fs，拟合结果见式(10),模型预测结果与当量蜡沉积厚度变化量见图6。

图6 模型预测结果与实际值对比图Fig.6 Comparison of model prediction results with actual values

(10)

711组数据计算结果中,有518组数据计算结果与实际变化趋势相同,占比72.9%,最大误差为1.74 mm/d,说明通过蜡沉积动力学预测模型计算当量蜡沉积厚度具有可行性。计算结果存在一定误差,原因如下。

1)现场提供生产数据不够精确,其中压力单位为MPa,输量单位为m3/h,难以准确计算管道当量蜡沉积厚度、范宁摩阻系数等参数。

2)本次研究选用第i-1天生产数据预测第i天当量蜡沉积厚度变化量,且第i-1天生产数据为各参数计量值的平均值,蜡沉积动力学预测模型未能考虑1天内工况变化(如停泵、增压、倒罐等)的影响。

3)现场管道跨度较长,但只能提供油品进站与出站数据,故研究中选用沿线油品平均温度计算密度、黏度与温度梯度等,不能准确反映温度变化对蜡沉积的影响。

4)JF管道投产至今已运行17年,管道内壁腐蚀与保温层失效情况严重,管段检修、更换频繁,各段管道粗糙度、传热系数以及内径不一,当量蜡沉积厚度计算时未能考虑上述情况。

受现场条件约束,蜡沉积动力学预测模型计算结果虽不能精准地根据第i-1天生产数据预测第i天蜡沉积厚度变化量,但计算结果与实际当量蜡沉积厚度变化趋势相同;因此,上述Hsu模型可用于评估工况(出站油温与输量)变化与地温变化对蜡沉积的影响。

4 结论与建议

1)当量蜡沉积厚度虽然不是真实蜡沉积厚度,但其能够反映蜡沉积对管道整体压降的影响;清管前后,当量蜡沉积厚度变化量与蜡沉积物质量分数呈线性正比,证明当量蜡沉积厚度计算结果准确性较高。

2)西北地区JF原油管道蜡沉积同时受分子扩散与剪切剥离作用影响,提高油品输量与出站温度能够有效降低管道当量蜡沉积厚度;选取Hsu模型预测蜡沉积速率;预测模型最大误差为1.74 mm/d,712组数据计算结果中,有518组数据计算结果与实际变化趋势相同,占比72.9%。

3)在本次研究中,现场数据质量较差导致当量蜡沉积厚度计算不够精确;蜡沉积动力学预测模型也未能考虑工况变化与内壁腐蚀等情况,因此以后的研究应该围绕提高计算精度与模型适应性两方面开展;而随着管网仿真技术的成熟与应用,蜡沉积速率动态预测、沉积厚度实时监控将成为可能。