烃储罐VOC排放建模：机理洞察与清单方法比较

《Computers & Chemical Engineering》：Modeling VOC emissions from hydrocarbon tanks: mechanistic insights vs. inventory approaches

　　

从储罐释放的挥发性有机化合物（VOCs）是炼油厂和石化设施中空气污染物的主要来源。本研究将Rota等人（2001）的唯象储罐顶部空间模型与美国环保局（U.S. EPA）的标准AP-42方法进行基准测试，使用了来自柴油储罐一年的小时级操作数据。前者公式通过数值稳定化（采用通量限制器）和冷凝子模型重新实施，提高了在陡峭梯度下的鲁棒性，并能够捕捉原始方案中缺失的可能蒸汽再冷凝。两种方法都再现了季节性排放循环，但在解析瞬态事件的能力上存在分歧：机理模型捕捉了填充/排放峰值和顶部空间动态，而AP-42主要仍受温度驱动。年度排放总量相当，AP-42为清单提供了一种快速、低输入的解决方案，而Rota等人（2001）为短期暴露提供了增强的时间分辨率。结果突显了计算需求与时间分辨率之间的权衡，对监管级与过程级应用具有影响。

**论文解读：烃储罐VOC排放建模——机理洞察与清单方法比较**

**研究背景**

挥发性有机化合物（VOCs）是造成空气污染的重要前体物，可导致地面臭氧和二次有机气溶胶（SOA）的形成，对人类健康及生态系统构成严重威胁。在炼油和石化设施中，烃类储罐是VOC的主要排放源之一，其排放机制包括静态呼吸损失（standing losses）和工作损失（working losses）。当前，量化储罐VOC排放的行业标准主要基于美国环保局（U.S. EPA）的AP-42方法，该方法采用半经验相关式，能够快速提供年度或月度平均排放清单。然而，AP-42的时间分辨率粗糙，且假设蒸汽空间均匀，难以捕捉由温度分层、产品切换或集中填充操作引起的短期排放峰值。为克服这些局限性，研究人员开发了基于瞬态质量和能量平衡的唯象模型，但此类模型计算成本高，且在实际工业条件下的验证仍有限。因此，开展这项研究的目的是将先进的唯象模型与标准的AP-42清单方法进行系统对比，评估两者在不同操作模式下的表现，以在计算需求与时间分辨率之间找到平衡。

**研究内容与意义**

研究人员重新实现了Rota等人（2001）提出的燃料储罐顶部空间模型，并对其进行了两项关键改进：引入总变差递减通量限制器（flux limiters）以提高数值稳定性，并增加了一个冷凝子模型以考虑蒸汽再冷凝现象。该改进模型与U.S. EPA的AP-42方法（第7章）进行了对比基准测试。研究使用了来自意大利南部一家炼油厂四座固定顶柴油储罐（TKA, TKB, TKC, TKD）一整年（8760小时）的小时级操作数据，包括液位、产品温度和周围环境温度。结论表明，两种方法均能再现排放的季节性温度调制，但唯象模型在解析瞬态事件（如填充/排放峰值）方面具有明显优势，而AP-42方法主要受温度驱动。年度排放总量基本相当，AP-42凭借其简单快速的特点适用于常规排放清单，而Rota等人的模型则更适合短期暴露评估和操作调度优化。这项研究的意义在于揭示了机理模型与清单方法之间的本质权衡：高时间分辨率以增加计算复杂度和输入数据需求为代价，这对监管层与过程层应用具有重要启示。该论文发表在《Computers》。

**关键技术方法**

研究人员主要采用了以下关键技术方法：
1. **数值稳定化**：引入总变差递减通量限制器（如Superbee限制器），在陡峭梯度处抑制数值振荡，保证长时间模拟的收敛性。
2. **冷凝子模型**：在每个时间步和空间节点，将计算出的蒸汽浓度与本地饱和浓度（Csat）比较，取两者较小值，确保物理一致性并避免非物理过饱和。
3. **坐标变换**：将物理垂直坐标无量纲化，将移动边界问题转化为固定计算域，简化数值实施。
4. **边界条件处理**：将微分边界条件转化为代数形式，避免高阶微分代数方程带来的数值困难。
5. **求解算法**：采用MATLAB R2024a编写有限差分代码，利用刚性求解器ode15s进行时间积分。样本队列来自意大利南部炼油厂的四座大型固定顶柴油储罐，数据通过炼油厂分布式控制系统（DCS）和当地气象站获取。

**研究结果**

**3.1 案例研究描述**
研究聚焦于四座操作模式截然不同的柴油储罐（TKA, TKB, TKC, TKD）。TKC和TKD展现几乎相反的行为：TKC液位剧烈变化（2米至17米），而TKD全年几乎恒定在10米。TKA和TKB操作模式中等，但TKA填充周期快速振荡，TKB则较为平滑。这些差异为评估不同模型的响应提供了基础。

**3.2 AP-42、Rota等人（2001）及冷凝模型的比较评估**
研究人员比较了四种方法：标准AP-42、带冷凝修正的AP-42、标准Rota模型、带冷凝修正的Rota模型。通过绘制全年屋顶蒸汽密度（g/m3）和计算年度排放总量进行分析。主要发现包括：
- 所有时间序列均呈现由环境温度驱动的季节性调制，夏季浓度高，冬季低。
- Rota模型对液位变化高度敏感，TKB与TKD之间的排放差异在Rota模型中显著，而AP-42因不依赖瞬时液位，对两者输出相似。
- 总体而言，Rota模型（含或不含冷凝）预测的排放通常低于AP-42。冷凝修正使AP-42的平均排放改变约11%，使Rota模型改变约14%；标准AP-42与带冷凝的Rota模型之间的平均差异约为17%。
- 在动态响应上，AP-42对温度变化反应迅速，而标准Rota模型表现出更强的惯性；添加冷凝后，Rota模型的瞬态峰值更尖锐但持续时间缩短。

**讨论与结论**

**讨论部分总结**：唯象模型对初始条件（如温度分布假设）敏感；冷凝处理能避免非物理过饱和并稳定数值结果。对于TKD储罐，由于持续微调导致的“液位抖动”引发高频扰动，Rota模型因其有限松弛时间能够捕捉众多部分消散的瞬态，累计年排放量反而高于AP-42。此外，蒸汽压相关性（PVA）的不确定性对排放估计影响极大，不同柴油模型之间的偏差可达53-55%，建议在清单或合规评估中使用校准的相关式或报告不确定性带。

**结论部分翻译**：本研究的目标是实现Rota等人（2001）的唯象储罐顶部空间模型，并在真实炼油厂环境中对其性能与U.S. EPA的AP-42进行基准测试。分析聚焦于四座固定顶柴油储罐，记录了一整年的小时级数据。计算框架（MATLAB 2024a）适应了现场条件，引入了通量限制器来处理陡峭梯度，并在屋顶处实施了简化的冷凝处理。模型输出随后以屋顶通风口浓度和相应排放（包括年度总量）进行评估。总之，两种方法都再现了排放的季节性温度调制。储罐间的差异主要由液位调度驱动，这在Rota等人（2001）中明确出现，而在AP-42中则不然。因此，Rota模型的输出中储罐间对比更明显，而AP-42时间序列主要受温度驱动。年度排放总量在方法间大致可比，Rota模型一般给出比AP-42更低的值，但存在储罐特定的例外。就既定目标而言，Rota等人（2001）模型在事件尺度分析上具有明显优势：它能在单一物理框架内解析填充/排放瞬态和顶部空间动态，并且一旦包含冷凝，可避免非物理过饱和并产生热力学一致的事件后松弛。这种瞬态解析能力使其特别适用于短期暴露评估、操作调度和峰值驱动事件的诊断分析。另一方面，对于清单类目标，该模型存在重要局限性：它适应操作变化缓慢，响应温度和液位波动比工厂数据观察到的更为迟缓。鉴于年度排放总量在方法间差异不大，瞬态解析求解器用于监管清单的性价比似乎有限；其增值在高等时间分辨率下最为清晰，而非年度汇总。此外，唯象方法对输入（初始和边界剖面、热特性表征和蒸汽压相关性选择）要求更高，如果现场特定数据稀缺，可能引入额外不确定性。总之，对于最初设定的目标，Rota等人（2001）模型在需要机械论的、瞬态解析的排放和峰值行为视角时最为有效，而U.S. EPA的AP-42因其简单、快速以及年度总量上的微小差异，仍然是年度清单开发的实用和适度工具。