含杂质CO2管道泄漏扩散模拟分析

陈兵 郭焕焕 崔维刚 肖红亮

西安石油大学

对CO2的捕集、利用和封存(CCUS)技术是目前减少CO2排放最经济有效的办法,其主要分为CO2的捕集、运输、封存和利用4个部分，而CO2的运输是CCUS技术中最为重要的一环[1-3]。对于长距离大规模CO2输送，管道输送有其独特的优势,即:输送量大、耗能低、经济性好，因此成为运输CO2最可行的方式[4-6]。CO2可以以液态、固态、气态、超临界/密相4种相态进行输送。超临界相态输送时，流体的密度和它的液态密度相近，黏度和气态黏度相近[7-8]，因此超临界相态是CO2管道运输最有效的状态。在CO2管道输送过程中，不可避免地会出现管道腐蚀、地表下陷、材料失效、第三方破坏等事故。管道破裂事故是极其危险的，其泄放的大量CO2及杂质会对周边动植物安全构成严重威胁[9-13]。因此，对管道气体泄漏扩散规律进行研究，对事故发生后的应急处理以及环境分析都有重要意义。

当CO2管道发生泄漏扩散时，杂质的潜在危害不可忽视，无论是天然气管道输送，还是CO2管道输送，杂质H2S作为一种剧毒酸性气体均普遍存在。因此，本研究取加拿大韦本地区CO2管道中的CO2气流，以H2S杂质为研究对象，进行泄漏扩散模拟计算，其气流组成见表1。

1 模型建立

PHAST(process hazard analysis software tool)是一种工艺危险源分析软件工具，是由DNV公司开发的一款专门用于天然气领域和石油石化危险分析和安全计算预测的软件。该软件对于气体的泄漏、放空扩散后果计算的针对性较强，其内部数据基于大量现场实验，尤其最新版专门对CO2泄漏、扩散模块进行了较大改善，可输出与气体浓度、范围等参数相关的曲线图，视觉效果和可读性较强[14-19]。通常情况下，在泄漏或放空过程中截止阀是工作的，泄放流体属于无源状态，本研究采用PHAST建立含杂质CO2管道泄漏扩散模型，选用软件中的短管(short pipeline)模型，假设管道发生全口径破裂，泄漏口将管道分成了上游分支和下游分支两个部分，如图1所示。

2 控制方程

CO2在泄漏或放空过程中，其流体流动仍遵循三大守恒定律：动量、能量、质量守恒定律，由此推导出CO2气体适用的最基本的连续方程、运动方程和能量方程。

(1)连续方程

(1)

式中：ρ为流体密度，kg/m3;t为时间，s；u、v、w分别为x、y、z轴的分量，该方程本身守恒。

(2)运动方程

(2)

(3)能量方程

(3)

式中：e为单位质量的内能，J;V为速度，m/s；ρ为流体密度，kg/m3;为拉普拉斯算子；为x轴的压力分量；为x轴的速度分量；γ为二阶单位张量；T为温度，K；qk为辐射加热，J。

(4)状态方程

本研究在计算CO2流体物性参数和CO2流体扩散过程中均选用P-R方程：

(4)

式中：μ、ω为常数，μ=2，ω=-1；a、b为物质特有常数，与临界参数有关，关系式如下：

(5)

(6)

式中：Tc为临界温度，K；Pc为临界压力，MPa；R为常数；T为温度,K；V为摩尔体积,L/mol。

3 模拟结果验证

大连理工大学与欧盟联合设计和搭建了世界上首套工业规模的CO2管道泄放装置，该装置开展了气相和超临界两种相态的纯CO2水平泄放实验，管道规格273 mm×20 mm，长258 m，主要实验条件见表2。

表2 实验条件Table 2 Experimental conditions条件初始温度/℃输送压力/MPa孔径/mm质量/t大气温度/℃湿度/%风速/(m·s-1)气相34.541001.528.00.9003.5超临界36.081004.031.80.5541.0

实验主要结论：对CO2浓度进行分析时，在气相下，CO2扩散长度至少为30 m，宽度至少3 m；在超临界相态下，长度至少50 m，宽度至少5 m。

本研究利用PHAST软件进行验证，选择短管和UDM模型，设置相同的实验条件，得出模拟结果，将实验数据与模拟数据进行对比(见表3)。

表3 实验、模拟结果对比mTable 3 Comparison of experiment and simulation results相态气相超临界长度宽度长度宽度实验30.03.050.05.0模拟31.02.259.04.2

由表3可以看出，用PHAST软件模拟的结果为：气相相态下CO2扩散的最大长度为31 m，最大宽度为2.2 m；超临界相态下最大长度为59 m，最大宽度4.2 m。总体上，模拟值略大于其实验值。但因为该实验属于工业级别，现场操作较复杂，而且实验所采用的双爆破装置会使管线压力具有不均匀性，因此实验值与模拟值存在一定误差。但该模拟结果符合喻健良等[14]的实验值预测范围，证明该模拟方法可行。

4 模拟结果与分析

本研究根据韦本气源组成进行模拟分析，其中CO2和H2S纯度分别为96%(摩尔分数)、0.9%(摩尔分数)。美国国家职业安全卫生研究所(NIOSH)认为，CO2浓度职业接触标准为5000×10-6(摩尔分数，下同)，危及健康的标准为20 000×10-6(最低危险浓度值)，一旦超过40 000×10-6，就会立即危及生命和健康。在较高的浓度下，杂质H2S也会对身体造成严重伤害，其职业接触标准值也即最低危险浓度值在10×10-6以下。因此，这里设置的CO2追踪浓度为20 000×10-6、40 000×10-6、60 000×10-6；H2S的追踪浓度为10×10-6、20×10-6、50×10-6。风速为1 m/s 时为无风状态。在不同孔径、不同风速、不同输送压力下，对管道泄漏扩散进行模拟分析，模拟结果见表4～表5，以及图2～图8。

表4 CO2泄漏扩散模拟结果Table 4 Simulation results of impurity CO2 leakage diffusion项目CO2扩散范围/m高度下风高度下风高度下风孔径/mm25024.94.417.31.913.31.130030.95.521.22.316.41.335037.06.625.22.719.51.640043.37.529.53.122.61.845049.88.233.93.525.72.1风速/(m·s-1)1107.111.769.54.552.32.6370.613.051.05.640.93.4555.314.841.76.234.53.9746.715.936.17.030.14.01038.616.630.47.325.84.5压力/MPa121.86.216.02.612.81.5230.56.321.32.616.61.5337.06.625.22.719.51.6442.76.728.92.822.01.7548.07.332.32.924.51.7

表5 H2S泄漏扩散模拟结果Table 5 Simulation results of impurity H2S leakage diffusion项目H2S扩散范围/m高度下风高度下风高度下风孔径/mm25062.1141.553.154.239.117.230077.2187.266.571.648.619.635092.1226.380.285.360.225.4400107.2285.794.6103.670.630.7450121.8327.6109.5119.880.734.2风速/(m·s-1)101 456.80295.6111.865.9301 225.20469.5103.975.8525.0705.773.8306.680.370.6750.8546.574.1251.266.264.91076.5463.362.4217.453.670.5压力/MPa144.1201.541.080.932.524.4272.7210.964.184.247.325.1392.5219.280.785.960.325.44108.9225.393.888.370.527.65121.3230.1107.692.580.629.3

经调研可知，CO2和H2S的管道泄漏浓度是对周围动植物造成严重威胁的最低浓度。因此，本研究选用最低浓度对危险范围进行分析。

4.1 孔径对泄漏扩散的影响

由图2和图3可以发现，在风速稳定时，随着孔径的增加，H2S和CO2扩散高度和顺风扩散距离显著增加，且两者都近似正比增加，孔径越大，扩散范围越大，但H2S的扩散高度约为CO2的两倍以上，H2S的扩散距离远远大于CO2的扩散距离。所以，同等浓度下，H2S扩散得更高、更远。因此，当CO2气流中含有较高浓度的H2S时，H2S扩散危险范围更大，安全范围的划定应以H2S为依据。

由图4可知，风速稳定时，随着孔径的增加，扩散到一定浓度时扩散的距离和扩散高度随着孔径的增加而增加，但增幅都较小，常用孔径300 mm到达浓度的时间为1.3 s，相邻孔径时间差均小于0.5 s。因此，孔径的变化对达到某一浓度的时间影响可忽略不计。

4.2 风速对泄漏扩散的影响

(1)由图5可知，在H2S、CO2扩散浓度下，风速越大，扩散高度的变化越小，这是由于风速越大，大气湍流强度越强，同时扩散初始动能被逐渐弱化。由图6可知，风速为1 m/s时，H2S的顺风距离达到了1500 m左右。显然，在相同的工况下，与CO2扩散的最低危险浓度相比，H2S扩散的职业接触浓度的范围更大。因此，在控制好危险范围的同时，更要做好气源地的除杂工作，严格限制H2S杂质的含量。

(2)同一条件下，风速越大，H2S最低危险浓度范围越小，随着风速的变化，低浓度变化范围要远大于高浓度，如风速为1～10 m/s，H2S下风向扩散距离的变化量为993.5 m，CO2对应的变化量为4.6 m。显然，较大的风速是有利于气体扩散的。

4.3 压力对泄漏扩散的影响

由图7和图8可以看出，输送压力对H2S扩散高度的影响与对CO2类似，输送压力沿射流方向的影响更显著。同等条件下，以3 MPa为例，考察最低危险浓度，H2S、CO2的扩散高度分别为93 m、37 m，下风向距离为219 m、6.6 m，显然H2S扩散范围更广。

同一扩散浓度下，随着输送压力的增加，CO2扩散高度和顺风扩散距离逐渐增加，但高度增加量明显更大，1～5 MPa过程增量为26.2 m，而下风向距离的增量仅为1.1 m，几乎不变，说明压力主要影响射流方向，射流方向与主扩散方向越趋近，扩散距离就越远。

5 不同截止阀间距结果分析

通过模拟计算分析可知，影响气体扩散的主要因素有风速、孔径大小、出口压力。通过分析三大因素对扩散范围的影响，并结合我国多数地区人口密度较高，以最大危险范围为计算目标，选取最常见的大气环境，采用韦本CO2长输管道常用的超临界/密相为模拟对象，将截止阀安装距离设置为8 km、10 km、12 km、15 km。

本研究关注的CO2最低浓度值为职业接触最大浓度值5000×10-6，H2S最低检测浓度为10×10-6，因此，5000×10-6为CO2职业接触浓度值，10×10-6为H2S最低危险浓度值。不同截止阀间距下CO2、H2S扩散持续到最低检测浓度的数据列于表6。

表6 不同截止阀间距扩散数据Table 6 Diffusion data at different stop valve spacing截止阀间距/kmCO2扩散H2S扩散最大距离/m云团宽度/m持续时间/s最大距离/m云团宽度/m持续时间/s8635.262.6976.61 100.1201.21 144.210662.466.51 302.51 156.4211.91 463.312691.370.01 666.41 208.4225.71 820.715728.173.92 274.81 277.4240.52 418.8

(1)由图9和图10可知，CO2和H2S扩散至最大距离和云团高度具有相同的变化趋势，持续时间随着阀间距增加呈线性增加，且增加幅度较大。间距为15 km时，最长CO2持续时间2 274.8 s，增量为(2 278.4 s-976.6 s)/(15 km-8 km)=185.5 s/km，最长H2S持续时间为2 418.8 s，增量为(2 418.8 s-1 144.2 s)/(15 km-8 km)=182.1 s/km，两者每千米时间增量相近，均在3 min左右。因此，阀间距越大，对人、动物和环境的影响越大，故在设置阀间距时持续时间是主要考虑对象。考虑到阀室安装成本及现场施工时间，可将职业接触浓度的持续时间降低至30 min，而阀间距为10 km时，CO2扩散距离为662 m，持续时间为1302 s，H2S扩散距离为1156 m，持续时间为1463 s，两相指标均不太高，因此可选择阀间距10 km作为参考标准。

(2)阀室间距的安装需要考虑当地的人口密度、经济性、放空安全范围、扩散至职业接触浓度值的持续时间等。综合前面的分析，我国不同地区人口密度变化较大，阀间距的安装可“因地制宜”，一般阀间距的安装以10 km为主。

6 结论

利用PHAST软件对地面架空管道泄漏扩散进行了模拟，研究了不同泄漏孔径、不同风速以及不同输送压力对泄漏扩散的危险范围影响。结果表明：

(1)随着孔径的增加，H2S和CO2扩散高度和顺风扩散距离显著增加，且两者都近似正比增加，但不同泄放孔径对管道泄漏扩散的影响不是很大。

(2)在相同的工况下，与CO2扩散的最低危险浓度相比，H2S扩散的职业接触浓度的范围更大。因此，在控制好危险范围的同时，更要做好气源地的除杂工作，同时随着风速的变化，CO2、H2S在低浓度情况下变化范围要远大于高浓度情况。

(3)同等条件下，输送压力对H2S扩散高度的影响与CO2类似，输送压力沿射流方向的影响更显著，压力主要影响射流方向，射流方向与主扩散方向越趋近，扩散距离越远。

(4)我国不同地区人口密度变化较大，CO2管道阀间距的安装可“因地制宜”，一般阀间距的安装以10 km为主。