石油与天然气化工  2019, Vol. 48 Issue (2): 111-115, 122
压力容器泄放过程建模与仿真分析
吕其宝1 , 张明洋2 , 陈丹霜2 , 宋满华1     
1. 中石化石油机械股份有限公司研究院;
2. 中石化石油机械股份有限公司三机分公司
摘要:在石油化工领域压力容器被广泛应用于存储和运输压缩气体。在特殊工况下,需要将压力容器气体进行泄放,以保证设备和管路的安全。基于热力学第一定律、质量守恒和AGA8状态方程,建立了压力容器泄放过程的动力学模型,分析了容器内气体压力、温度、压缩因子和绝热指数等参数的变化规律,研究了外界温度、气体组分和孔板面积对泄放时间的影响,给出了计算孔板直径的方法。结果表明,泄放过程中,压缩因子和绝热指数是变值,气体组分对泄放时间有明显的影响,在仿真计算时应该充分考虑。泄放时间随外界温度升高而增加。该研究结果可以为压力容器放空系统的设计和优化提供参考和帮助。
关键词压力容器    气体组分    物性参数    AGA8状态方程    动力学模型    
Modeling and simulation analysis of relief process of pressure vessel
Lyu Qibao1 , Zhang Mingyang2 , Chen Danshuang2 , Song Manhua1     
1. Institute of Sinopec Oilfield Equipment Corporation, Wuhan, Hubei, China;
2. San Ji Branch Company of Sinopec Oilfield Equipment Corporation, Wuhan, Hubei, China
Abstract: In petrochemical field, the pressure vessels are widely applied for storing and transporting compressed natural gas. In order to guarantee the safety of equipment and pipeline, the gas filled in the pressure vessels should be released under some special conditions. Based on the first law of thermodynamic, mass balance and AGA8 equation of state, a novel dynamic model was established to stimulate the relief process in this paper. According to this model, the variation of the compressor parameters, such as pressure, temperature, compressor factor and adiabatic index were predicted. Meanwhile, the influences of external temperature, gas component and orifice area on discharge time were investigated. Besides, a method for calculating orifice diameter was presented. The simulation results showed that the compressor factor and adiabatic index are variables and the gas component has significant effect on discharge time, which should not be ignored in calculations. Moreover, the discharge time increased with the elevation of the external temperature. All the above findings could provide reference and support for the design and optimization of the pressure vessel vent system.
Key words: pressure vessel    gas component    physics property    AGA8 equation of state    dynamic model    

随着天然气开采技术的发展和市场对天然气需求量的不断增加,高压存储系统在气田集输站场和石化机械制造领域受到越来越多的关注。压力容器作为高压存储系统的一种,由于其良好的存储性能和灵活的结构样式而被广泛应用。在一些特殊工况下,例如,天然气泄漏或者站场发生火灾等安全隐患,需要对压力容器进行安全及时的放空,以保证集输站场或机械设备的安全,防止损失和破坏的进一步扩大。国内外学者对压力容器的泄放进行了大量的研究[1-3], 研究方法主要有以下两种:一是按照相应标准与设计手册计算[4-5];二是利用商业软件仿真模拟[6-10]。黄翼翔等[6]比较了常用限流孔板的计算方法,并结合软件模拟和生产实际,提出了孔板优化的计算公式。陈俊文等[7]利用HYSYS软件计算了火灾工况下紧急放空系统限流孔板孔径。杨天宇等[9]利用HYSYS软件仿真模拟了火灾工况下泄放质量流量的变化规律。

压力容器中的高压气体组分和物性参数会影响容器泄放过程,成为研究的一个重要方面。AGA8状态方程是美国天然气协会推荐的计算天然气压缩因子的方程,该方程具有应用范围广和计算精度高等特点,是计算气体物性参数的常用方法[11-13]。目前,关于气体组分和物性参数对压力容器泄放过程影响的研究相对较少。因此,有必要进一步深入分析。

本研究将基于热力学第一定律、质量守恒定律和AGA8状态方程建立压力容器泄放过程动力学模型,仿真分析容器内气体压力、压缩因子等参数的变化规律以及气体组分和孔板面积等参数对泄放时间的影响,为压力容器放空系统的设计和优化提供理论指导。

1 动力学模型

将压力容器视为一个开放系统,在泄放过程中忽略泄漏影响。本节将利用热力学第一定律、质量守恒定律和AGA8状态方程建立压力容器泄放过程的动力学模型。

1.1 热力学第一定律方程

根据热力学第一定律,压力容器内能量的变化满足[14]

$ \frac{\mathrm{d} Q}{\mathrm{d} t}+h_{\rm s} \frac{\mathrm{d} m_{\rm s}}{\mathrm{d} t}=h_{\mathrm{e}} \frac{\mathrm{d} m_{\mathrm{e}}}{\mathrm{d} t}+\left(p_{\mathrm{c}} \frac{\mathrm{d} V_{\mathrm{c}}}{\mathrm{d} t}+\left|F_{\mathrm{f}} \frac{\mathrm{d} C}{\mathrm{d} t}\right|\right)+ \\ h_{\mathrm{c}} \frac{\mathrm{d} m_{\mathrm{c}}}{\mathrm{d} t}-Z_{\mathrm{c}} R \frac{\mathrm{d} T_{\mathrm{c}}}{\mathrm{d} t} $ (1)

式中: $\frac{\mathrm{d} Q}{\mathrm{d} t}$为容器与外界的热传递; $p_{\mathrm{c}} \frac{\mathrm{d} V_{\mathrm{c}}}{\mathrm{d} t}+\left|F_{\mathrm{f}} \frac{\mathrm{d} C}{\mathrm{d} t}\right|$为气体做功; $h_{\rm s} \frac{\mathrm{d} m_{\mathrm{s}}}{\mathrm{d} t}$为进入容器气体的能量; $h_{\mathrm{e}} \frac{\mathrm{d} m_{\mathrm{e}}}{\mathrm{d} t}$为流出容器气体的能量; $h_{\mathrm{c}} \frac{\mathrm{d} m_{\mathrm{c}}}{\mathrm{d} t}-Z_{\mathrm{c}} R \frac{\mathrm{d} T_{\mathrm{c}}}{\mathrm{d} t}$为容器内气体能量的变化; 下标s为进气,c为压力容器条件,e为排气。

压力容器泄放过程中,式(1)可以简化为:

$ \frac{\mathrm{d} Q}{\mathrm{d} t}=h_{\mathrm{e}} \frac{\mathrm{d} m_{\mathrm{e}}}{\mathrm{d} t}+h_{\mathrm{c}} \frac{\mathrm{d} m_{\mathrm{c}}}{\mathrm{d} t}-Z_{\mathrm{c}} R \frac{\mathrm{d} T_{\mathrm{c}}}{\mathrm{d} t} $ (2)
1.2 质量守恒方程

压力容器作为控制体,则质量守恒方程可以写为:

$ \frac{\mathrm{d} m_{\mathrm{c}}}{\mathrm{d} t}=\frac{\mathrm{d} m_{\mathrm{s}}}{\mathrm{d} t}-\frac{\mathrm{d} m_{\mathrm{e}}}{\mathrm{d} t} $ (3)

压力容器泄放时,有dms/dt=0。因此,式(3)可以简化为:

$ \frac{\mathrm{d} m_{\mathrm{c}}}{\mathrm{d} t}=-\frac{\mathrm{d} m_{\mathrm{e}}}{\mathrm{d} t} $ (4)

假定气体在孔板中为绝热流动,单位时间内流经孔板的气体质量流量可以表示为[15]:

$ \frac{\mathrm{d} m_{\mathrm{e}}}{\mathrm{d} t}=p_{\mathrm{c}} \gamma A \sqrt{\frac{2 k}{(k-1) Z_{\mathrm{c}} R T_{\mathrm{c}}}\left[\frac{p_{\mathrm{e}}}{p_{\mathrm{c}}}\right)^{\frac{2}{k}}-\left(\frac{p_{\mathrm{e}}}{p_{\mathrm{c}}}\right)^{\frac{k+1}{k}} ]} $ (5)

式中:γ为流量系数; A为孔板面积; k为绝热指数。

1.3 传热方程

压力容器内气体的热传递过程可以表示为[14]

$ \frac{\mathrm{d} Q}{\mathrm{d} t}=\alpha A_{0}\left(T_{\mathrm{c}}-T_{\mathrm{w}}\right) $ (6)

式中:αA0TcTw分别为热传递系数、气体的接触面积、压力容器内气体温度与外界温度,由于外界温度变化较小,在仿真时看作常数。

1.4 AGA8状态方程

压力容器泄放是一个动态过程,随着压力和温度的变化,气体的物性参数也随之变化。根据上面的热力学第一定律方程和质量守恒方程可知,仿真时需要计算压缩因子、绝热指数等参数。AGA8方程是国际标准化组织推荐的一种计算工作状态下天然气物性参数的方程,具有计算精度高,使用范围广的特点。本研究采用AGA8状态方程对相关物性参数进行计算。

AGA8状态方程的一般形式为[11]

$ p=\rho_{\mathrm{m}} Z R T $ (7)
$ Z=1+B \rho_{\mathrm{m}}-\rho_{r} \sum\limits_{n=13}^{18} C_{n}^{*}+ \\ \qquad\sum\limits_{n=13}^{58} C_{n}^{*}\left(b_{n}-c_{n} k_{n} \rho_{r}^{k_{n}}\right) \rho_{r}^{b_{n}} \exp \left(-c_{n} \rho_{r}^{k_{n}}\right) $ (8)
$ \rho_{r}=K^{3} \rho_{\mathrm{m}} $ (9)
$ \begin{aligned} C_{n}^{*}=& a_{n}\left(G+1-g_{n}\right)^{g_{n}}\left(Q^{2}+1-q_{n}\right)^{q_{n}} \cdot \\ &\left(F+1-f_{n}\right)^{f_{n}} U^{u_{n}} T^{-u_{n}} \end{aligned} $ (10)
$ B=\sum\limits_{n=1}^{18} a_{n} T^{-u_{n}} \sum\limits_{i=1}^{N} \sum\limits_{j=1}^{N} x_{i} x_{j} B_{n i j}^{*} E_{i j}^{u_{n}}\left(K_{i} K_{i}\right)^{3 / 2} $ (11)
$ B_{n i j}^{*}=\left(G_{i j}+1-g_{n}\right)^{g_{n}}\left(Q_{i} Q_{j}+1-q_{n}\right)^{q_{n}} \cdot \\ \qquad \left(F_{i}^{1 / 2} F_{j}^{1 / 2}+1-f_{n}\right)^{f_{n}}\left(S_{i} S_{j}+1-s_{n} \right)^{s_{n}} \cdot \\ \qquad \left(W_{i} W_{j}+1-w_{n}\right)^{w_{n}} $ (12)

式中:Z为压缩因子;B为第二维利系数;xi为组分i在气体混合物中的摩尔分数;anbncnknungnqnfnsnwn为常数;KiQiFiWi为组分i的组合参数;EijGij为二元交互作用参数。

2 数值仿真与结果分析
2.1 仿真流程

动力学模型数值仿真流程如图 1所示。首先,根据输入的初始压力p、初始温度T和气体组分x,利用AGA8状态方程(7)和(8)计算压缩因子Z、气体密度ρm和绝热指数k。其次,将计算得到的物性参数代入热力学第一定律方程(2)和质量守恒方程(4),得到下一时刻压力容器的温度T和密度ρm。利用AGA8状态方程计算下一时刻压力容器压力p,并与排气压力pd进行比较,如果ppd,重复上面的循环,直到ppd,计算结束。最后,输出压力容器内压力、温度以及物性参数数据。

图 1     压力容器动力学模型仿真流程 Figure 1     Simulation flow of the pressure vessel dynamic model

2.2 仿真结果与分析
2.2.1 仿真参数设置

本节将利用前面建立的动力学模型分析泄放过程中压力容器内气体状态参数、物性参数,质量流量等变化趋势。另外,研究不同节流面积、外部温度和气体类型对泄放时间的影响。压力容器内的气体组分如表 1所列,结构参数和其他仿真参数如表 2所列。

表 1    气体1各组分的摩尔分数 Table 1    Mole fraction of gas 1

表 2    仿真参数 Table 2    Parameters for simulation

2.2.2 仿真结果分析

为了验证AGA8状态方程结果的准确性,分别计算了温度300 K和压力2 MPa条件下,气体密度和压缩因子仿真结果与专业物性计算软件之间的误差值,如图 2所示。当温度为300 K,压力为1~10 MPa时,误差值小于0.25%;当压力为2 MPa,温度为300~340 K时,误差值小于0.08%。通过误差分析可知,本研究使用AGA8状态方程计算压力容器内气体物性参数变化是合理可行的。

图 2     AGA8状态方程误差分析 Figure 2     Error analysis of AGA8 equation of state

图 3所示为泄放过程中压力容器内气体的压力、温度与时间的关系。从图 3可知,容器内气体泄放时间为5.25 min。随着时间的增加,气体压力逐渐减小,压力越大,单位时间内压力下降得越快。另外,通过观察气体温度的变化趋势可以发现,当时间小于2.5 min时,气体的温度随时间的增加而逐渐降低,但是,当大于2.5 min时,气体的温度随着时间的增加出现小幅的上升。造成这种现象的原因为:一方面,在气体泄放过程中,容器内气体的内能减小,温度降低;另一方面,气体与周围环境进行热传递,在开始阶段,气体温度高于外界温度,热量由气体传递给周围环境,气体温度降低, 当气体温度低于环境温度时,气体从环境中吸热,吸收的热量大于泄放减少的内能时,容器内气体温度上升。

图 3     气体压力与温度变化曲线 Figure 3     Curves of gas pressure and temperature

压力容器内气体的物性参数变化曲线如图 4所示。从图 4(a)可以看出,气体压缩因子和绝热指数并不是一个定值,而是随着泄放时间不断变化。压缩因子随着时间的增加先减小后增大,而绝热指数的变化趋势与压缩因子正好相反。因此,在仿真时,不应该忽略压缩因子和绝热指数的影响。

图 4     气体物性参数变化曲线 Figure 4     Curves of gas physical parameters

气体的质量流量和流速变化如图 4(b)所示。根据式(5)可知,气体流经节流孔板的质量流量与容器内气体压力、压缩因子、绝热指数等因素有关。其中,气体压力的影响最大,气体压力随时间的增加而逐渐减小,因此气体质量流量随时间增加而减小。另外,通过观察发现,质量流量与时间近似成线性关系。气体流速的变化趋势与质量流量的变化趋势相同。

由于压力容器内气体与外界存在温度差,在泄放过程中存在热传递现象。根据热力学第一定律可知,容器与外界的热传递会引起系统的压力变化,从而影响气体的流速、排气温度以及实际泄放时间等参数。图 5所示为不同外界温度条件下,压力容器泄放时间的变化规律。随着外界温度的升高,泄放时间逐渐增加。产生这种现象是因为外界温度升高,气体吸收的热量增加,容器内气体单位时间总内能变化减小,导致压力降低速率变慢,因此气体泄放时间随温度的升高而增加。

图 5     泄放时间与外界温度变化曲线 Figure 5     Variation of discharge time for ambient temperature

根据质量守恒方程和AGA8状态方程可知,孔板面积和容器内气体组分是影响泄放时间的重要参数。例如,相同压力和温度条件下,气体组分不同,由状态方程可知,气体的密度、压缩因子等物性参数不同。因此,在相同结构参数下,当容器内气体组分变化时,泄放时间也会发生变化。图 6所示为孔板面积、气体组分与泄放时间的关系曲线。表 3表 4所列分别为气体2和气体3的摩尔分数。通过比较可以发现,在相同条件下,不同的气体组分泄放时间差值最大可达0.5 min。其中,气体1的泄放时间最长,在泄放时间一定的情况下,对于气体1应该适当增加孔板面积。通过图 6还发现,容器泄放时间随孔板面积的增大而减小,当孔板面积较小时,增加孔板面积可以明显缩短泄放时间。但是,当孔板面积超过一定数值后,采用增加孔板面积来缩短泄放时间的效果减弱。另外,当气体组分确定和泄放时间确定时,可以通过图 6中的曲线确定需要的最小孔板面积。

图 6     泄放时间与孔板面积、气体组分变化曲线 Figure 6     Variation of discharge time and gas components for orifice area

表 3    气体2各组分的摩尔分数 Table 3    Mole fraction of gas 2

表 4    气体3各组分的摩尔分数 Table 4    Mole fraction of gas 3

3 结论

根据热力学第一定律、质量守恒定律和AGA8状态方程,建立了压力容器泄放过程的动力学模型。通过该模型,可以对泄放过程中压力容器内气体状态参数和物性参数的变化进行实时仿真。结论如下:

(1) 模型考虑了泄放过程中压缩因子、绝热指数、泄放速度和质量流量的动态变化,更接近真实情况。

(2) 模型分析了外界温度和气体组分对泄放时间的影响。相同结构参数条件下,外界温度越高,泄放时间越长;不同气体组分,泄放时间不同。在结构设计时,不能忽略这些因素的影响。

(3) 模型计算了不同孔板面积对泄放时间的影响,得到了孔板面积与泄放时间的关系,可以为孔板设计和优化提供参考与指导。

参考文献
[1]
WON K W, SMITH A R, ZEININGER G A. Thermodynamic methods for pressure relief system design parameters[J]. Fluid Phase Equilibria, 2006, 241(1/2): 41-50.
[2]
HUFF J E. Multiphase flashing flow in pressure relief systems[J]. Plant/Operations Progress, 1985, 4(4): 191-199. DOI:10.1002/(ISSN)1549-4632
[3]
刘熔, 李峰, 盛炳林. 天然气净化厂超压自动放空系统设置[J]. 天然气与石油, 2006, 24(5): 37-40. DOI:10.3969/j.issn.1006-5539.2006.05.013
[4]
API. Pressure-relieving and depressuring systems: ANSI/API STD 521-2007[S]. UK: API, 2007.
[5]
《炼油装置工艺管线安装设计手册》编写组. 炼油装置工艺管线安装设计手册[M]. 北京: 石油工业出版社, 1978: 7-8.
[6]
黄翼翔, 陈俊文, 陈庆. 限流孔板孔径计算优化[J]. 质量管理, 2013(9): 179-181.
[7]
陈俊文, 陈庆, 汤晓勇, 等. 紧急放空系统限流孔板孔径HYSYS模拟计算[J]. 石油与天然气化工, 2016, 45(6): 93-98. DOI:10.3969/j.issn.1007-3426.2016.06.019
[8]
陈文峰, 刘培林, 郭洲, 等. 复杂物系压力容器安全阀泄放过程的HYSYS动态模拟[J]. 天然气与石油, 2010, 28(6): 55-58. DOI:10.3969/j.issn.1006-5539.2010.06.015
[9]
杨天宇, 朱海山, 郝蕴, 等. 压力容器火灾工况安全泄放质量流量的动态研究[J]. 石油与天然气化工, 2014, 43(2): 208-212. DOI:10.3969/j.issn.1007-3426.2014.02.022
[10]
李涛, 李昱江, 苏敏, 等. HYSYS动态模拟方法进行安全阀尺寸分析[J]. 油气田地面工程, 2015, 34(5): 22-24. DOI:10.3969/j.issn.1006-6896.2015.5.011
[11]
中国石油天然气总公司.天然气压缩因子的计算第2部分: 用摩尔组成进行计算: GB/T 17747.2-1999[S].北京: 中国标准出版社, 2004.
[12]
梁光川, 左果. 基于气体组成的天然气压缩因子计算方法[J]. 石油与天然气化工, 2014, 43(4): 395-400. DOI:10.3969/j.issn.1007-3426.2014.04.011
[13]
李丹华, 姜东琪. AGA8-92DC计算方法天然气压缩因子计算[J]. 煤气与热力, 2011, 31(3): A35-A38.
[14]
JIN J M, HONG W R. Valve dynamic and thermal cycle model in stepless capacity regulation for reciprocating compressor[J]. Chinese Journal of Mechanical Engineering, 2012, 25(6): 1151-1160. DOI:10.3901/CJME.2012.06.1151
[15]
吴沛宜, 马元. 变质量系统热力学及其应用[M]. 北京: 高等教育出版社, 1983.