SLAB 用户手册:模拟重气体泄漏的空气扩散模型中文简要用户使用手册环境保护部环境工程评估中心国家环境保护环境影响评价数值模拟重点实验室手册说明本用户手册基于《USER’S MANUAL FOR SLAB: AN ATMOSPHERIC DISPERSION MODEL FOR DENSER-THAN-AIR RELEASES》(1990.06)编写，仅对美国EPA网站所提供的模拟重气体泄漏的空气扩散模型SLAB的使用方法提供中文版简要说明，更详细的程序使用说明请查阅相关的软件手册及文档，或采用带图形界面版的商业软件。
本手册由环境保护部环境工程评估中心国家环境保护环境影响评价数值模拟重点实验室负责编写，参与人员包括：易爱华、陈陆霞、胡翠娟、梁昊、杨晔、丁峰等。
本手册版权所有，转载及印刷请与环境保护部环境工程评估中心联系。
本手册所涉及的模型系统及本手册电子版本下载地址：一、SLAB简介SLAB是用于模拟重气体泄漏的空气扩散模型。
该模型最初基于Zeman于1982年提出的关于重气体云的空气卷吸和重力扩散的理念而开发。
SLAB早期相关工作由美国能源署支持。
SLAB的进一步开发由USAF工程和服务中心（1986年开始）和美国石油学会（1987开始）共同提供支持。
现行的SLAB版本可以模拟连续的、限时的和瞬时的物质泄漏，泄漏源包括以下4种：地面液池蒸发、高于地面的水平射流、烟囱或高于地面的垂直射流，以及瞬时释放。
SLAB除可以用于模拟重气体的扩散，还可以模拟中性浮力气体的烟云扩散，以及烟云轻于空气时的上升过程。
泄漏时的空气扩散过程可以通过求解质量、动量、能力和物质的守恒方程来计算，如图1所示。
为了简化守恒方程的求解过程，方程可以通过将烟云作为稳态烟羽或瞬时烟团在空间上进行平均。
连续排放（持续时间非常长的排放源）可以作为稳态烟羽。
有限时间的排放采用稳态烟羽模式描述最初烟云的扩散，而且在该排放源持续泄漏的时间段内，可以一直使用稳态烟羽模式。
释放一旦终止，烟团被视为瞬时烟团，之后的扩散采用瞬时烟团模式来计算。
对于瞬时泄漏的排放源，整个过程都均使用瞬时烟团扩散模式。
二、理论介绍2.1重气体扩散模型简介重气体泄漏的空气扩散模型受到多种物理现象的影响，这些物理现象在中性或浮力气体泄漏中可能不会发生或者即便发生也不是很重要。
这些物理现象包括：重气体烟云的稳定密度分层导致的湍流衰减；由于重力流和初始排放源动量导致的环境速率场的改变；由于液滴形成和挥发以及在过热或低温液体排放情况下的地面加热对烟云温度、浮力和湍流的热力学效应；此外，我们所关注的某种特定的重气体的浓度可能和典型大气污染物关注的累积浓度差别很大。
例入，对于易燃气体，关注的是瞬间浓度；而对于有毒气体，关注的则是几分钟到几小时的浓度，以及累积浓度。
因此，为了能够更好的预测出重气体泄漏时有毒浓度区的大小和持续时间，所有重要的物力现象都需要进行考虑，而且预测过程中要使用最合理的浓度平均时间。
为了满足重气体泄漏情形的要求，SLAB模型以质量、动量、能量和物质守恒方程的平均形式为起点，在该理论框架的基础上进行开发（如图1所示）。
这些方程用于计算扩散烟团的空间平均性质，并且以两种方式来代表两种不同的扩散模式：稳态烟羽扩散模式和瞬时烟团扩散模式。
2.2模型组成SLAB模型的计算流如图2所示。
一个典型的模拟过程主要包括三个层次：源类型识别和初始化，烟云扩散计算和时间平均浓度计算。
通过守恒方程求解，获得瞬时空间平均气云特性，从而计算得到时间平均体积浓度。
2.3控制扩散方程2.3.1稳态烟羽模式SLAB的稳态烟羽模式基于稳态侧风向平均的质量、动量、能量和物质守恒方程，并且使用空气卷吸概念来考虑气体云与环境大气的湍流混合，如图3所示。
图3 SLAB 烟羽扩散模型预测的重气云扩散2.3.2瞬时烟团模式Slab的瞬时烟团模式基于体积平均浓度的质量、动量、能量和物质守恒方程，采用空气卷吸理念来考虑云和周围空气的湍流混合。
图4 SLAB 烟团扩散模型预测的重气体云扩散2.3.3烟羽向烟团转化模式瞬时排放源或短时间持续的液池蒸发可以采用烟团模式进行模拟，同时泄漏结束或者稳态期结束后的泄漏源也可以采用烟团模式进行模拟。
在后一种情况下，当扩散方程由稳态烟羽方程向瞬时烟团方程转变时，存在空间平均气云性质在计算上的转变。
稳态烟羽模式下，在气云的侧风向上进行空间平均；瞬时烟团模式下，在气云的整个体积上进行空间平均。
因此，为了开始烟团模式的计算，需要定义这一转变时刻的传输时间，该时间的云长度，以及此时的质量中心。
三、用户指南3.1输入文件运行SLAB需要约30个输入参数。
这些参数包括源类型、源性质、泄漏性质、场地性质、气象参数等。
表1 输入变量的定义源类型和数值子步骤参数—IDSPL泄漏源类型1——液池泄漏2——水平喷射泄漏3——垂直喷射或烟囱泄漏4——瞬时或短持续时间液池挥发—NCALC数值子步骤参数源性质—WMS泄漏物质的分子质量（KG）—CPS恒压下的蒸汽热容（J/KG K）—TBP沸点温度（K）—CMEDO初始液体质量比—DHE汽化热（J/KG）—CPSL液体热容（J/KG K）—RHOSL源物质的液体密度（KG/M3）—SPB饱和气压常数（默认：SPB=-1.0）—SPC饱和气压常数（默认：SPC=0.0）泄漏参数—TS初始泄漏温度（K）—QS源泄漏速率（KG/S）—AS源初始扩散面积（M2）—TSD连续源持续时间（S）—QTIS瞬时源泄漏量（KG）—HS源高度（M）场地参数—TAV浓度平均时间（S）—XFFM最大下风向距离（M）—ZP(I)浓度计算高度（M）；I=1，4气象参数—ZO地表粗糙度（M）—ZA环境测量高度（M）—UA环境风速（M/S）—TA环境温度（K）—RH相对湿度（百分娩比）—STAB稳定度等级等级值描述A-F 1.0-6.0不稳定—稳定默认0.0输入稳定度参数“ALA”—ALA M-O长度的倒数（1/M）（ALA是仅当STAB=0.0成立时的输入参数）3.1.1源类型和数值子步骤参数3.1.2.1IDSPL-源类型SLAB可处理四种类型的源：1-液池泄漏，2-水平喷射泄漏，3-垂直喷射或烟囱泄漏，4-瞬时或短持续时间的液池泄漏。
3.1.2.2NCALC-数值子步骤参数参数NCALC是一个整数子步骤乘数，它指定了在守恒方程积分期间执行的计算子步骤的数量。
通常推荐NCALC=1。
3.1.2源性质3.1.2.1WMS——源物质的分子质量（kg）3.1.2.2CPS——恒压下的蒸汽热容（J/kg K）3.1.2.3TBP——源物质的沸点温度（K）3.1.2.4CMEDO——初始液体质量比其中CPSL是物质在液相中的比热，DHE是在沸点温度TSP下的挥发热。
3.1.2.5DHE——沸点温度下的汽化热（J/kg）3.1.2.6CPSL——源物质的液体比热（J/kg K）3.1.2.7RHOSL——源物质的液体密度（kg/m3）3.1.2.8SPB，SPC——饱和蒸气压常数SLAB中使用的关于饱和蒸气压的饱和蒸气压常数表达式如下：其中PA是环境气压（在SLAB中经常是PA=101325 N/m3=1 atm），T是当地云温度，以K为单位。
当饱和气压常数未知时，可用默认的选项，指定SPB值为“-1.0”及SPC值为“0.0”。
此时将使用Clapeyron方程来定义SPB值。
3.1.3泄漏参数3.1.3.1TS——初始泄漏温度（。
K）初始泄漏温度的定义（TS）取决于泄漏的类型。
当泄漏是挥发池（IDSPL=1或4），TS 是沸点温度TBP。
当瞬时泄漏（IDSPL=4）时，TS是物质在排放瞬间的温度，或者当源是爆炸的结果时，TS是物质在完全膨胀并降至一个大气压时的温度。
对于喷射泄漏（IDSPL=2或3），源相关参数为物质在其完全膨胀后的性质。
其性质为可根据公式进行计算。
当源物质作为蒸汽加压储存时，它们作为蒸汽排放（CMEDO=0.0），此时推荐将膨胀作为绝热过程处理。
此时TS由下式给出：其中γ=Cp/Cv是比热的比值，Pa是环境大气压，Pst和Tst分别是储存压强和温度。
如果对于蒸汽排放计算的初始泄漏温度（TS）在沸点温度（TSP）以下，那么TS应该设置等于沸点温度。
类似地，当源物质作为液体加压储存并作为两相液滴—蒸汽混合物排放时，源温度是沸点温度TBP。
初始泄漏温度（TS）必须等于或大于沸点温度TBP，因为源或者完全在蒸汽相（TS≥TBP）或者是平衡的液滴和蒸汽相的混合物（TS=TBP）。
3.1.3.2QS—源泄漏速率（kg/s）是指连续泄漏的排放速率，即液池排放（IDSPL=1）、喷射排放中的一种（IDSPL=2或3）以及短持续时间液池挥发（IDSPL=4）。
对于瞬时排放（IDSPL=4），源排放速率应该设置为等于0（QS=0.0）。
3.1.3.3AS—源初始扩散面积（m2）源初始扩散面积对于不同泄漏类型有不同定义。
对于液池挥发（IDSPL=1或4），AS是挥发池的面积。
如果AS未知，它可以通过有效挥发速度WS使用以下公式来进行计算：其中QS是输入源排放速率，RHOS是源物质在沸点温度TBP下的蒸汽密度，WS是以速度表达的已知挥发速度（m/s）。
蒸汽密度RHOS由理想气体定律计算：其中WMS是输入的源物质的分子质量，Pa是环境大气压（Pa=101325 N/m2），Rc是气体常数（Rc=8.31431 J/mol K），TBP是输入沸点温度。
当源是压力状态下的水平或垂直喷射排放（IDSPL=2或3）时，AS是完全膨胀后，气压降至环境大气压水平后的源面积。
如果源物质作为纯蒸汽储存和排放（CMEDO=0.0），建议将膨胀视为绝热过程。
源面积可以表达为：11 / 13其中Pst 是储存压强，Pa 是环境大气压，TS 是输入源温度，Tst 是储存温度，Ap 是破裂或开口的实际面积。
当源物质作为液体加压储存并作为两相喷射排放时，AS 是闪蒸并且形成纯物质液滴—蒸汽混合物后的源面积。
在这种情况下，AS 值由下面公式给出：其中RHOSL 是输入源物质的液体密度，Ap 是破裂或开口的实际面积，ρm 值可以由状态方程计算并且将其中一些项进行重组其中CMEDO 是输入初始液体质量比，RHOSL 是输入源物质的液体密度，RHOS 是源物质在沸点温度下的蒸汽密度，如本节中前面部分所述。
在瞬时泄漏（IDSPL=4）的情况下，AS 是中心在点X=Y=Z=0.0的地面上的体源的面积。
源面积定义为：其中Va 是瞬时排放的体积，HS 是体积的高度，QTIS 是输入排放质量，ρsi 是初始排放密度。
当源是纯蒸汽排放时，ρsi 是纯物质在源温度TS 下的蒸汽密度，由下式给出：当源是液体—蒸汽混合物时，ρsi 是在沸点温度TSP 下的混合物密度ρm ，其液体质量比为CMEDO 。