某油田大型联合站多排放源VOCs扩散的数值模拟

周娟; 何战友; 徐自强; 朱恒; 王鑫雅; 张澄; 黄维秋

doi:10.12030/j.cjee.202210037

某油田大型联合站多排放源VOCs扩散的数值模拟

1.
中国石油长庆油田分公司油气工艺研究院，西安 710000
2.
低渗透油气田勘探开发国家工程实验室，西安 710000
3.
常州大学，江苏省油气储运技术重点实验室，常州 213164

作者简介: 周娟 (1983—) ，女，硕士，高级工程师，zhoujuan_cq@petrochina.com.cn

通讯作者: 黄维秋 (1965—)，男，博士，教授，hwq213@cczu.edu.cn

基金项目:
国家自然科学基金资助项目 (51574044)
中图分类号: X51;TE86

Numerical simulation of VOCs diffusion from multiple emission sources at a large combined station in an oilfield

1.
Oil & Gas Technology Research Institute of Changqing Oilfield Branch of CNPC, Xi’an 710000, China
2.
National Engineering Laboratory for Exploration and Development of Low Permeability Oil & Gas Fields, Xi’an 710000, China
3.
Changzhou University, Jiangsu Key Laboratory of Oil Vapor Storage and Transportation Technology, Changzhou 213164, China

Corresponding author: HUANG Weiqiu, hwq213@cczu.edu.cn

摘要: 油田联合站是油田的重要组成部分，存在油气 (VOCs) 排放的可能，有必要对其扩散规律进行研究，以制定相关污染防控及安全措施。以某典型联合站为例，建立1:1的实际模型，结合现场调研测试和数值模拟，重点分析了联合站正常工况下多排放源的VOCs扩散机理及储罐裂缝处VOCs泄漏扩散的叠加效应。结果表明：在风速影响下，罐间和背风侧由于出现绕流和回流，容易达到爆炸极限；在重力和涡流的作用下，联合站内的背风侧会出现VOCs聚集，但叠加效应不明显；当储罐发生罐壁破损时，在事故罐后方出现正压区，在涡流和强气流的影响下，叠加后的VOCs浓度会明显增强，并呈点射状向下风向扩散，油气爆炸危险区域加速扩展。本研究成果可为联合站设计、运行管理及制定安全环保措施提供参考。
- 油田联合站 /
- 油气(VOCs)扩散 /
- 储罐泄漏 /
- 数值模拟 /
- 叠加效应
Abstract: Oil field joint station is an important part of oil field, there is the possibility of oil and gas (VOCs) emission, it is necessary to carry out research on the oil vapor emission laws of oilfields. Therefore, taking a typical united station as an example, a 1:1 practical model was established. Combined with on-site investigation, testing and numerical simulation, an investigation of the oil vapor diffusion mechanism under normal conditions and the superposition effect of oil vapor leakage and diffusion at the crack of the storage tank was conducted emphatically. The results showed that: under the influence of wind speed, the oil vapor could easily reach the explosion limit due to the circumfluence and backflow between the tanks and the leeward side; under the action of gravity and eddy current, oil vapor would accumulate on the leeward side, but the superposition effect was not obvious; When the tank wall was damaged and leaked, a positive pressure area would appear behind the accident tank. Under the influence of eddy current and strong airflow, the superimposed oil vapor concentration would be significantly enhanced, and will spread to the downwind in a point-like manner, accelerating the extension of oil vapor explosion danger area. The results of this study can provide reference for the design, operation management and safety and environmental protection measures of the joint station.
- oilfield united station /
- oil vapor diffusion /
- tank leakage /
- numerical simulation /
- superposition effect
我国城镇污水处理厂常年在实际运行过程中会遭受到有机负荷冲击^[1-4]，进水污染物浓度的突然变化或有机负荷(OLR)冲击会影响生化系统中微生物的生命活动、营养物质的利用及微生物代谢途径，进而影响水中污染物的去除效果^[5-6]。我国城镇污水处理厂常用生物工艺为An/O、氧化沟、SBR和MBR等，An/O和SBR为城镇污水厂主要工艺^[7]，污水厂在运行的过程中常受高有机负荷冲击，生物系统遭受破坏且出水水质不达标，且高有机负荷冲击导致系统长时间处于低溶解氧状态，易造成污泥膨胀。因此如何采取有效的措施应对有机污染物冲击，保证出水水质稳定达标，对于实际工程应用价值更大。

董国日等^[8]通过研究连续有机负荷冲击对SBR反应器的影响，当周期有机冲击负荷为0.10 kg COD/(kg MLVSS·h)时，SBR系统能够耐受连续多个周期的高有机负荷冲击。金明姬等^[9]通过研究不同有机负荷条件下曝气时间对SBR反应器的影响，当搅拌/曝气时间分别为2.5 h/7.5 h时，COD和TP的去除效果最佳。本文通过改变进水有机负荷的大小，来了解不同有机负荷瞬时冲击对SBR反应器的影响，给出系统自然恢复时间。通过检测进出水COD、NH₄⁺-N、TN、TP和活性污泥胞外聚合物(EPS)，并采取提高曝气量的方法探索其对各项污染物的去除效果，以此为工程实践提供理论参数支持。

1.   材料与方法

1.1   试验装置

试验装置示意图，见图1。SBR反应器采用有机玻璃板制，300 mm×300 mm×400 mm，总有效容积为36 L。试验期间水温（21±1.5）℃。采用电磁式空气压缩机曝气，使用水泵进水，并用玻璃转子流量计进行气量控制。反应器运行方式：曝气8 h、沉淀3 h、排水10 min、闲置50 min，排水比50%。每个周期12 h，采用非限制性曝气方式。反应器活性污泥来自吉林省柏林水务集团污水厂。

图 1 试验装置

下载: 全尺寸图片幻灯片

1.2   实验用水

实验污水模拟城镇污水主要成分：淀粉、牛肉膏、蛋白胨、乙酸钠(CH₃COONa)、氯化铵(NH₄Cl)、磷酸二氢钾(KH₂PO₄)、碳酸氢钠(NaHCO₃)、硫酸镁(MgSO₄)、氯化钙(CaCl₂)、硫酸亚铁(FeSO₄)和微量元素1 mL。其中微量元素成分：氯化锌(ZnCl₂)、硫酸铜(CuSO₄)、硫酸锰(MnSO₄)和氯化铝(AlCl₃)。实验期间进水污染物浓度，见表1。

表 1 设计进水主要污染物浓度范围 mg·L⁻¹

污染物常负荷负荷1 负荷2 负荷3

COD 372.1～427.3 798.2～811.2 1 116.3～1 281.9 1 484.8～1 709.2
NH₄⁺-N 33.53～36.17
TN 38.53～41
TP 3.02～4.17

| Show Table

DownLoad: CSV

1.3   试验方法

冲击试验之前反应器已稳定运行30 d，运行期间曝气量1.2 L/min，混合液悬浮固体(MLSS)质量浓度3 500 mg/L，污泥龄30 d，运行期间平均出水COD、NH₄⁺-N、TN和TP分别为28.9、0.43、10.46和0.083 mg/L。反应器运行稳定之后进行单周期瞬时有机负荷冲击试验，实验通过调整淀粉、乙酸钠的量来控制进水有机负荷大小。冲击过程中观察不同负荷下污染物去除变化和污泥特性，给出系统自然恢复时间并结合反应器出水效果给出相应的解决对策。反应器运行参数，见表2。

表 2 运行参数

指标常负荷负荷1 负荷2 负荷3

曝气/L·min⁻¹ 1.20 1.20 1.20 1.50 1.80 2.10
污泥负荷/g·(L·d)⁻¹ 0.17 0.34 0.51 0.68

| Show Table

DownLoad: CSV

1.4   检测方法

实验过程中MLSS采用重量法测定^[10]；COD采用重铬酸钾法测定；TN采用碱性过硫酸钾消解紫外分光光度法测定；NH₄⁺-N采用纳氏试剂分光光度法测定；TP采用钼酸铵分光光度法测定；DO和温度用WTW(oxi3310)便携式DO仪检测测定；EPS采用甲醛-NaOH法提取^[11]；PN含量采用Folin-酚法测定^[12]；EPS中PS含量采用苯酚-硫酸法测定^[13]，PN和PS的含量以MLSS计，mg/g。

2.   结果与分析

2.1   有机负荷冲击对污染物去除效果的影响

2.1.1   不同负荷下污染物去除效果

有机负荷对COD的降解。反应初期去除COD主要以物理吸附为主，活性污泥吸附量随着负荷的增加而增大。0～4 h期间常负荷、负荷1和负荷2基本完全降解，4～8 h期间COD处于平稳状态，出水COD分别为31.6、42.14和31.6 mg/L。0～8 h时负荷3在降解的过程中，随着降解时间的增加，COD的降解速度变慢，出水COD为96.32 mg/L。

图2（b）可知，NH₄⁺-N完全降解的时间随着负荷的增加而增加，负荷3的NH₄⁺-N出水浓度为7.205 mg/L，末端溶解氧偏低，硝化反应无法充分进行。图3（c）可知，常负荷、负荷1和负荷2出水TN为11.045、8.302和6.59 mg/L ，结合图2（b）和图3可看出硝化反应充分进行，不同负荷冲击下为反硝化提供了充足的碳源，且随着有机负荷的增加而增加；DO偏低为反硝化提供了良好的缺氧环境，且充足的碳源有利于异养菌的增长，使得反硝化反应充分进行。图2（d）可知，TP降解历时中，PAOs的释放量随着负荷的增加而增加，负荷1、负荷2和负荷3情况下的TP检测最大值分别为16.096、37.4和43.6 mg/L，3种负荷条件下出水TP分别为0.056、0.104和19.3 mg/L。负荷3的TP出水超高的原因是进水负荷的增加，高浓度有机物能使聚磷菌(PAOs)更好的生成聚羟基烷酸，厌氧释磷不受影响可以很好的进行^[14]。PAOs在厌氧条件下得到充分的碳源，释磷时间也随之增加；导致好氧吸磷的时间较短，出水TP浓度过高，其次高负荷冲击下（图3），冲击过程中DO一直处在偏低状态，异养菌和硝化细菌争夺DO，硝化反应发生不完全。由于进水有机物浓度高、曝气时间短不利于硝化细菌的生长，而残留在反应器内的硝化反应产物硝态氮抑制PAOs厌氧释磷^[15]。另外，PAOs也需在好氧条件下吸收磷来合成ATP，无法为好氧吸磷和硝化反应提供足够的溶解氧。试验结果表明，系统在OLR达到0.68 g/(L·d)时反应器出水不满足《城镇污水处理厂污染物排放标准：GB 18918—2002》的一级A标准。通过上述研究发现DO对冲击过程影响较大，因此决定加大曝气量增加DO来改善氧环境，探索增加DO对活性污泥系统抗OLR冲击缓解作用大小。

图 2 不同负荷下的污染物降解历时

下载: 全尺寸图片幻灯片

2.1.2   自然恢复下污染物去除效果

图4（a）可知，反应器经过4个周期COD出水达到48.19 mg/L，又经过7个周期自然恢复，恢复速率呈下降的趋势。第10个周期到第11个周期恢复速率逐渐平缓，且第11个周期时反应器出水28.67 mg/L，满足GB 18918—2002排放标准。

图 3 不同负荷时的溶解氧变化

下载: 全尺寸图片幻灯片

图 4 自然恢复下系统出水的变化

下载: 全尺寸图片幻灯片

图4（b）可知，反应器经过3个周期NH₄⁺-N出水4.964 mg/L，又经过3个周期自然恢复，恢复速率呈上升的趋势，在第11个周期时反应器NH₄⁺-N出水达到0.682 mg/L。图4（b）和图4（c）可知，随着自然恢复期数的增长TN浓度也随之增长，NH₄⁺-N浓度随之下降。原因是高有机负荷冲击下反应器溶解氧过低，硝化反应无法充分进行，高浓度有机物冲击下有利于异养菌的生长，且在低溶解氧条件下更有利于反硝化的进行。系统自然恢复的过程中随着周期数的增加，微生物得到充分的溶解氧，硝化反应进行充分，反硝化受到抑制造成硝态氮和亚硝态氮堆积，所以TN浓度也随之增高。图4（d）可知，反应器经过9个周期TP出水达到0.346 mg/L，又经过3个周期自然恢复出水达到0.054 mg/L，且自然恢复速率逐渐下降。原因是高有机负荷冲下，反应器在高浓度有机物条件下好氧异养菌与硝化菌以及PAOs争夺DO,抑制了PAOs好氧吸磷过程，恢复常负荷后PAOs得到充足的溶解氧且维持生命活动所消耗的有机物的量充足，通过自然恢复11个周期反应器出水TP满足城镇污水处理厂污染物排放标准。

2.2   有机负荷冲击对EPS的影响

EPS是活性污泥在一定条件下分泌的、包裹在胞外的高分子聚合物，主要成分为PS、PN等^[16-18]。活性污泥在受到水质、水量的冲击条件下，EPS能对微生物起着保护作用。由于反应系统在负荷3的条件下系统出水不达标，因此通过检测冲击前后EPS含量来得出活性污泥系统自然恢复周期。负荷3冲击前PS和PN分别为16.62和14.62 mg/g，经过5个周期后PS达到最大值为30.17 mg/g，PN无较大变化，为16.48 mg/g，又经过5个周期后PS降至19.16 mg/g，EPS恢复至冲击前要需11个周期。负荷3时EPS恢复至起始时的含量，在恢复过程中，PS含量变化较大，PN含量变化较小，有机负荷较高的情况下，微生物无法利用全部碳源来维持新陈代谢，而是将过量的碳源转化为PS。系统自然恢复的过程中，高有机负荷冲击下促进系统产生大量的异养型微生物，PS又作为碳源被系统微生物利用而逐渐减少，见图5。

图 5 自然恢复下系统EPS的变化

下载: 全尺寸图片幻灯片

2.3   曝气量对有机负荷冲击的调控作用

根据以上分析得出反应器在负荷3时脱氮除磷效率下降，主要原因在于高有机负荷冲击下反应器DO的不足，异养菌与硝化菌以及聚磷菌争夺DO,抑制了硝化作用和好氧吸磷过程。因此通过提高曝气量来解决高有机负荷冲击下反应器DO偏低的问题，并给出应对该问题的解决参数。

图6（a）可知，COD降解速率随着曝气量的增大而变快，曝气量分别为1.5、1.8和2.1 L/min，平均出水COD分别为39.66、28.25和31.59 mg/L。当气量增至2.1 L/min时，COD的降解速率相比于1.8 L/min时无变化。图6（b）可知，NH₄⁺-N降解速率随着曝气量增加而增加，反应器得到充足的DO，硝化反应进行充分。在相同曝气时间时，出水NH₄⁺-N浓度降低，曝气量分别为1.5、1.8和2.1 L/min时，平均出水NH₄⁺-N分别为5.28、2.39和0.97 mg/L。图6（c）可知，气量越大TN降解速率越慢，原因是随着曝气量的增加反应器DO充足，厌氧时间缩短，不利于反硝化的充分进行。曝气量分别为1.5、1.8和2.1 L/min时平均出水TN分别为8.59、9.89和11.5 mg/L。图6（d）可知，随着曝气量的增加，PAOs得到充分的DO；厌氧释磷时间随曝气量增加而缩短。从图6中可以明显看到，不同的气量条件下磷的释放量也不同，气量越大磷的释放量相对越大，说明水中溶解氧越高，水力剪切力越大；把淀粉大分子物质分解成小分子糖类被PAOs储存，LIU et al ^[19]研究表明水力剪切力越大可促进PS分泌；因此PS也为厌氧释磷提供充足的碳源。当气量增至1.8 L/min时，反应器DO不足，硝化细菌和PAOs竞争反应器的溶解氧，此时DO不充足导致好氧吸磷无法充分进行。曝气量分别为1.5、1.8和2.1 L/min时平均出水TP分别为8.73、2.64和0.447 mg/L，提高曝气可使出水TP浓度降低。曝气量为2.1 L/min，污染物去除效率最高并且反应器出水达到城镇污水处理厂污染物排放标准。综合分析，负荷3下最适曝气量应为2.1 L/min。

图 6 负荷3下不同曝气量时的污染物降解历时

下载: 全尺寸图片幻灯片

3.   结论

（1）SBR活性污泥系统受单周期瞬时高浓度有机负荷冲击时，在0.68 g/(L·d)负荷条件下，反应8 h时的出水COD、NH₄⁺-N、TN和TP分别为96.32、7.205、7.34和19.3 mg/L，不满足城镇污水处理厂污染物排放标准。DO不充足影响有机物的降解、硝化反应和好氧吸磷，且充足的碳源有利于反硝化和厌氧释磷进行。

（2）在0.68 g/(L·d)负荷条件下，系统冲击前PS和 PN分别为16.62和14.62 mg/L，冲击后EPS组分中PN含量无明显变化，PS含量经过11个周期自然恢复后；在第5个周期PS达到最大值。经过11个周期自然恢复后反应器出水COD、NH₄⁺-N、TN和TP分别为28.67、0.682、11.12和0.054 mg/L满足GB 18918—2002排放标准，系统自然恢复至冲击前需要11个周期。

（3）瞬时有机负荷冲击导致反应器内DO不足，污染物去除效率降低的情况可采用提高曝气量的方式进行调控。在0.68 g/(L·d)负荷条件下，随着气量增加污染物去除效率提高，TN去除效率逐渐降低。通过调整曝气量的方式，在曝气量为2.1 L/min时，系统在0.68 g/(L·d)负荷条件下出水污染物满足污水处理厂污染物排放标准。

图 1 联合站建筑编号图

Figure 1. Building number of united station

下载: 全尺寸图片幻灯片

图 2 联合站三维几何模型图与计算域

Figure 2. 3D geometric model diagram and computational domain of united station

下载: 全尺寸图片幻灯片

图 3 联合站三维网格划分结果

Figure 3. Results of three-dimensional meshing of united station

下载: 全尺寸图片幻灯片

图 4 现场测得G3净化罐VOCs排放组成

Figure 4. VOCs emission composition of G3 purification tank measured at the united station

下载: 全尺寸图片幻灯片

图 5 各测点浓度对比

Figure 5. Concentration comparison of each monitoring point

下载: 全尺寸图片幻灯片

图 6 联合站建筑压力和XY截面压力分布云图 (v= 2 m·s⁻¹)

Figure 6. Cloud map of building pressure and XY section pressure distribution in united station

下载: 全尺寸图片幻灯片

图 7 不同风速下联合站内气流流动轨迹 (t=300 s，y= 0.5 m)

Figure 7. Air flow trajectory in united station at different wind speeds

下载: 全尺寸图片幻灯片

图 8 不同高度联合站内2 m·s^-1的速度分布云图

Figure 8. Cloud map of speed distribution of 2 m·s^-1 in the union station at different heights

下载: 全尺寸图片幻灯片

图 9 不同风速下联合站内的速度分布云图 (t=300 s，y= 0.5 m)

Figure 9. Cloud map of speed distribution in the union station under different wind speeds

下载: 全尺寸图片幻灯片

图 10 不同风速下联合站内VOCs质量分数分布云图 (t=300 s)

Figure 10. Cloud map of VOCs mass fraction distribution in the union station at different wind speeds（ t =300 s）

下载: 全尺寸图片幻灯片

图 11 不同高度下联合站内压力分布云图 (t=300 s)

Figure 11. Cloud map of pressure distribution in the union station at different heights (t=300 s)

下载: 全尺寸图片幻灯片

图 12 不同高度下联合站内VOCs质量分数分布云图

Figure 12. Cloud map of VOCs mass fraction distribution in the union station at different heights

下载: 全尺寸图片幻灯片

表 1 通用方程 (3) 中各参数的具体形式

Table 1. Detailed form of each parameter in general equation (3)

方程	广义变量Φ值	广义扩散系数Г值	广义源项S
质量守恒方程	1	0	0
动量守恒方程	${u_i}$	$\mu$	$- \partial p/\partial x+{S_i}$
能量守恒方程	$T$	$k/c$	${S_T}$
组分运输方程	${C_S}$	${D_S}\rho$	${S_S}$

方程	广义变量Φ值	广义扩散系数Г值	广义源项S
质量守恒方程	1	0	0
动量守恒方程	${u_i}$	$\mu$	$- \partial p/\partial x+{S_i}$
能量守恒方程	$T$	$k/c$	${S_T}$
组分运输方程	${C_S}$	${D_S}\rho$	${S_S}$

下载: 导出CSV

表 2 G区储罐的尺寸

Table 2. Dimension of oil storage tank in area G

储罐位号	储罐直径/m	储罐高度/m	公称容积/m³	最大液体高度/m	平均液体高度/m
G1	17	13.15	3 000	12	10.68
G2	17	13.2	3 000	12	9
G3	17	13.14	3 000	12	9

下载: 导出CSV

表 3 数值模拟网格无关性检验结果

Table 3. Results of grid independence test for numerical simulation

网格数量	沿x方向的速度/ (m·s⁻¹)
746 528	3.14
1 461 100	3.23
2 005 657	3.64
3 009 842	3.64

下载: 导出CSV

表 4 某油田联合站现场测试数据

Table 4. Data from field test of an united station

区域	采样点位	VOCs质量浓度/(mg·m⁻³)	采样时刻	风速/(m·s⁻¹)
I区	G3罐顶环境	11.40	14:45	2.0
	G2罐顶环境	9.95	14:40	2.0
	B环境 (停用)	24.50	15:15	0.5
	I区厂门口	0.95	18:05	1.5
Ⅱ区	D5罐顶环境	20.00	11:51	3.2
Ⅱ区	Ⅱ区罐区环境	11.00	16:15	3.0

下载: 导出CSV

[1]	祁甲民. 基于多源测量数据融合的三维建模技术研究[D]. 郑州: 河南工业大学, 2018.
[2]	LI T, ZHU J, ZHANG W. Cascade utilization of low temperature geothermal water in oilfield combined power generation, gathering heat tracing and oil recovery[J]. Applied Thermal Engineering, 2012, 40: 27-35. doi: 10.1016/j.applthermaleng.2012.01.049
[3]	LI Z, LIANG Y, WANG G, et al. A method for optimising pump configuration and operation in oilfield water injection network[J]. Chemical Engineering Transactions, 2021, 88: 1105-1110.
[4]	LI Z, GUO Y, XU N, et al. Integration of a novel distributed water and energy system in the oilfield[J]. Chemical Engineering Research and Design, 2022, 186: 350-361. doi: 10.1016/j.cherd.2022.08.008
[5]	张秀玲, 宋翠红, 刘春杨. 储油罐油气扩散规律的理论研究[J]. 中国石油和化工标准与质量, 2011, 31(7): 248-249. doi: 10.3969/j.issn.1673-4076.2011.07.213
[6]	TAMADDONI M, SOTUDEH-GHAREBAGH R, NARIO S, et al. Experimental study of the VOC emitted from crude oil tankers[J]. Process Safety and Environmental Protection, 2014, 92(6): 929-937. doi: 10.1016/j.psep.2013.10.005
[7]	黄维秋. 油气回收基础理论及其应用[M]. 北京: 中国石化出版社, 2011: 244-253.
[8]	ATKINSON G, COWPE E, HALLIDAY J, et al. A review of very large vapour cloud explosions: Cloud formation and explosion severity[J]. Journal of Loss Prevention in the Process Industries, 2017, 48: 367-375. doi: 10.1016/j.jlp.2017.03.021
[9]	MILAZZO M F, ANCIONE G, LISI R. Emissions of volatile organic compounds during the ship-loading of petroleum products: Dispersion modelling and environmental concerns[J]. Journal of environmental management, 2017, 204: 637-650. doi: 10.1016/j.jenvman.2017.09.045
[10]	CIRIMELLO P G, OTEGUI J L, RAMAJO D, et al. A major leak in a crude oil tank: Predictable and unexpected root causes[J]. Engineering Failure Analysis, 2019, 100: 456-469. doi: 10.1016/j.engfailanal.2019.02.005
[11]	SIMAYI M, SHI Y, XI Z, et al. Emission trends of industrial VOCs in China since the clean air action and future reduction perspectives[J]. Science of The Total Environment, 2022, 826: 153994. doi: 10.1016/j.scitotenv.2022.153994
[12]	WU T, CUI Y, LIAN A, et al. Vehicle emissions of primary air pollutants from 2009 to 2019 and projection for the 14th Five-Year Plan period in Beijing, China[J]. Journal of Environmental Sciences, 2023, 124: 513-521. doi: 10.1016/j.jes.2021.11.038
[13]	GOEURY C, HERVOUET J-M, BAUDIN-BIZIEN I, et al. A Lagrangian/Eulerian oil spill model for continental waters[J]. Journal of hydraulic Research, 2014, 52(1): 36-48. doi: 10.1080/00221686.2013.841778
[14]	刘泽阳. 联合站典型源VOCs排放规律研究及应用[D]. 青岛: 中国石油大学(华东), 2019.
[15]	HUANG W Q, HUANG F Y, FANG J, et al. A calculation method for the numerical simulation of oil products evaporation and vapor diffusion in an internal floating-roof tank under the unsteady operating state[J]. Journal Of Petroleum Science And Engineering, 2020, 188: 106867. doi: 10.1016/j.petrol.2019.106867
[16]	黄维秋, 陈风, 吕成, 等. 基于风洞平台实验的内浮顶罐油气泄漏扩散数值模拟[J]. 油气储运, 2020, 39(4): 425-433.
[17]	DEAVES D, GILHAM S, MITCHELL B, et al. Modelling of catastrophic flashing releases[J]. Journal of hazardous materials, 2001, 88(1): 1-32. doi: 10.1016/S0304-3894(01)00284-9
[18]	ZHANG T, LI G, YU Y, et al. Atmospheric diffusion profiles and health risks of typical VOC: Numerical modelling study[J]. Journal of Cleaner Production, 2020, 275: 122982. doi: 10.1016/j.jclepro.2020.122982
[19]	BELLEGONI M, OVIDI F, LANDUCCI G, et al. CFD analysis of the influence of a perimeter wall on the natural gas dispersion from an LNG pool[J]. Process Safety and Environmental Protection, 2021, 148: 751-764. doi: 10.1016/j.psep.2021.01.048
[20]	OKAMOTO K, ICHIKAWA T, FUJIMOTO J, et al. Prediction of evaporative diffusion behavior and explosion damage in gasoline leakage accidents[J]. Process Safety and Environmental Protection, 2021, 148: 893-902. doi: 10.1016/j.psep.2021.02.010
[21]	KOZO F. Progress and future prospects of CFD in aerospace—Wind tunnel and beyond[J]. Progress in Aerospace Sciences, 2005, 41(6): 455-470. doi: 10.1016/j.paerosci.2005.09.001
[22]	ALEXANDROS K, ALEIFERIS P G, CHARALAMBIDES A G. Numerical investigation of VOC levels in the area of petrol stations[J]. Science of the Total Environment, 2014, 470: 1205-1224.
[23]	黄维秋, 吕成, 郭淑婷, 等. 油气排放及回收的研究进展[J]. 石油学报(石油加工), 2019, 35(2): 421-432.
[24]	HUANG W Q, WANG S, JING H B, et al. A calculation method for simulation and evaluation of oil vapor diffusion and breathing loss in a dome roof tank subjected to the solar radiation[J]. Journal Of Petroleum Science And Engineering, 2020, 195: 107568. doi: 10.1016/j.petrol.2020.107568
[25]	许雪, 陈风, 黄维秋, 等. 基于风洞平台实验的大型罐区溢油事故后的油气扩散模拟[J]. 环境工程学报, 2021, 15(12): 3946-3956. doi: 10.12030/j.cjee.202101099
[26]	ALEXANDER M, LORENZO M, VAGESH D N. Comparison of k-ε models in gaseous release and dispersion simulations using the CFD code FLACS[J]. Process Safety and Environmental Protection, 2019, 130: 306-316. doi: 10.1016/j.psep.2019.08.016

期刊类型引用(5)

1.	胡昊杨，吴哲，侯成亮. 基于收尘管路设计优化和清灰系统的研究. 建筑机械. 2024(07): 13-15 . 百度学术
2.	李建龙，陈源正，林子捷，吴泉泉，吴代赦. 除尘滤筒脉喷清灰技术研究进展与展望. 金属矿山. 2022(11): 23-35 . 百度学术
3.	郗元，姜文文，代岩，王国际，闫志刚，任福良，牛凤娟. 基于CFD的锥形散射器强化清灰特性数值模拟及优化. 轻工机械. 2021(01): 98-103 . 百度学术
4.	邱俊，陈强，李建龙，吴代赦. 数值模拟滤筒中负载尘饼对脉冲喷吹清灰效果的影响. 环境工程学报. 2021(03): 1067-1074 . 本站查看
5.	莫布林. 基于设备组合控制方式的袋式除尘器脉冲喷吹控制系统. 东北电力技术. 2021(05): 28-31 . 百度学术

其他类型引用(14)

点击查看大图

图( 12) 表( 4)

计量

文章访问数: 2881
HTML全文浏览数: 2881
PDF下载数: 133
施引文献: 19

1. 材料与方法
1.1 试验装置
1.2 实验用水
1.3 试验方法
1.4 检测方法
2. 结果与分析
2.1 有机负荷冲击对污染物去除效果的影响
2.2 有机负荷冲击对EPS的影响
2.3 曝气量对有机负荷冲击的调控作用
3. 结论

全文HTML

石油工业是我国的重要支柱产业之一，油田联合站是油田原油集输和处理的中枢^[1]。在联合站中，各区块来油通过原油稳定系统将分离出的轻烃组分回收进入轻烃厂，处理后的原油组分通过联合站处理再进行外输^[2-4]。在联合站的日常操作运行中，各种拱顶罐呼吸阀、干化池等存在油气逸散的风险。另外，储罐在运行过程中可能因腐蚀、雷击、油罐失效、人为损坏等意外因素导致油罐罐壁破损^[5-6]，从而泄漏和蒸发出大量易燃易爆、有毒的油气混合物。这些混合物即挥发性有机物 (volatile organic compounds，VOCs) 。VOCs的泄漏不仅会污染环境，还会严重威胁企业安全和人群健康、安全和环境 (Health、Safety、Environment，HSE) ^[7-10]。“十四五”规划提出，要强化多污染物协同控制和区域协同治理，加快VOCs排放综合整治，到2025年，VOCs排放总量下降10%以上^[11-12]。经调研，油田联合站的多个环节 (如原油储罐及轻烃储罐、除油罐、净化罐、沉降罐、敞开液面) 存在VOCs排放问题，因此，有必要开展油田联合站VOCs排放的扩散规律研究。

目前，VOCs扩散规律分析方法主要有：现场测试^[13-14]、风洞实验^[15-16]和数值模拟^[17-20]。其中，基于计算流体力学 (computational fluid dynamics，CFD) 及其Fluent软件的数值模拟方法的可操作性强，已在航天、汽车、能源、化工、材料、生物医药等诸多领域广泛应用^[21]。KOUNTOURIOTIS等^[22]通过数值模拟，在风速、风向、温度、油气扩散源位置等多种影响因素下，研究了不同成分的汽油挥发出的VOCs扩散规律，并发现在扩散源附近的VOCs浓度远远高于爆炸极限。基于风洞平台实验验证和Fluent数值模拟，建立了基于单膜传质理论的油气蒸发过程当量膜厚数值模拟计算方法，以及基于Stefan-Fuchs方程、Clausius-Clapeyron方程及若干准则数的非稳态蒸发单相传质的数值模拟方法，揭示了在各操作条件下油罐非稳态石油蒸发和油罐排放气在大气环境中的扩散行为及其内在机理，以及影响因素间的关联性^[15,23-24]。

本课题组通过CFD 数值模拟和风洞实验平台研究了某石化企业的实体罐区发生溢油事故后油气蒸发的扩散规律，掌握了罐区VOCs浓度的变化特征^[25]。然而，针对油田联合站内多排放源的VOCs扩散规律研究仍未出现。本研究以某油田大型联合站多点排放源为对象，通过CFD数值模拟和现场调研数据相结合，以探究正常工况下VOCs扩散机理及储罐裂缝处VOCs泄漏扩散的叠加效应。本研究结果可与前期针对大型罐区VOCs扩散数值模拟^[25]成果相结合形成系列成果，为石油石化行业的运行管理及VOCs排放控制提供参考。

4. 结论

1) 在正常工况条件下，联合站内多排放源的受风速影响，内储罐迎风侧的压力达到最大，罐间和背风侧由于出现绕流和回流，出现大面积的负压区域，且有明显的反方向速度。另外，在重力和涡流的影响下，VOCs呈整体向下运动趋势，容易造成罐间和背风侧的VOCs聚集，且在漩涡引起的强气流作用下，Ⅱ区和某石化企业容易达到爆炸极限。但由于联合站内排放源之间的距离大于扩散范围，VOCs叠加效果不明显。

2) 通过对联合站Ⅰ区的排放源G2罐罐壁破损泄漏时进行模拟，发现在事故罐后方会出现红色正压区，且VOCs扩散会呈点射状向下风向扩散。加上“卡门涡街”效应和漩涡引起的强气流，多排放源的VOCs扩散会出现“1+1>2”的叠加效应。叠加后的VOCs质量浓度会明显增强，油气爆炸危险区域加速扩展。另外，Ⅱ区和某石化企业会处在爆炸极限范围内，应注意火灾爆炸的风险。

3) 本研究考虑条件为从联合站Ⅰ区到Ⅱ区的风向。但联合站常年风向不定，若联合站风向为Ⅱ区到Ⅰ区，按本模拟方法可推断处联合站内多排放源的VOCs扩散会对联合站Ⅰ区的影响较大，并严重影响作业区工作人员的健康与生产作业安全，且应注意火灾爆炸的风险。

参考文献 (26)

某油田大型联合站多排放源VOCs扩散的数值模拟

作者简介: 周娟 (1983—) ，女，硕士，高级工程师，zhoujuan_cq@petrochina.com.cn

通讯作者: 黄维秋 (1965—)，男，博士，教授，hwq213@cczu.edu.cn

Numerical simulation of VOCs diffusion from multiple emission sources at a large combined station in an oilfield

Corresponding author: HUANG Weiqiu, hwq213@cczu.edu.cn

1. 材料与方法

1.1 试验装置

1.2 实验用水

1.3 试验方法

1.4 检测方法

2. 结果与分析

2.1 有机负荷冲击对污染物去除效果的影响

2.1.1 不同负荷下污染物去除效果

2.1.2 自然恢复下污染物去除效果

2.2 有机负荷冲击对EPS的影响

2.3 曝气量对有机负荷冲击的调控作用