泄漏风险背景与模拟参数构建
氯甲烷作为基础化工原料，在常温常压下呈气态，其密度大于空气，泄漏后极易在低洼处积聚形成重气云团。储罐区管道接口、阀门法兰等部位在长期交变应力作用下存在微小泄漏隐患。重气扩散不同于中性气体浮力扩散，其核心物理特征在于重力主导的沉降滑塌效应，这会导致云团向四周蔓延的横向距离远大于抬升高度。在复杂罐区地形中，储罐器具群会对风场产生显著的阻挡与尾流扰动，改变局部湍流强度，进而影响氯甲烷浓度的空间分布。采购方在评估罐区安全防护距离时，若仅依赖理想平坦地形下的经验公式，往往会严重低估下风向的毒害区域范围。 开展精确的扩散推演，必须构建高保真物理模型。遵照GB/T 37243-2019《危险化学品危险源定量风险评估方法》（现行有效）要求，模型输入端需采集罐区周边连续一年的气象统计资料，涵盖风速、风向、大气稳定度联合频率分布。地形高程数据需反映储罐、防液堤、管廊等障碍物的几何特征。在设定泄漏孔径时，我们选取了具有代表性的5mm裂纹场景，该尺寸相当于成年人小拇指末节长度，能够真实反映器具微小故障导致的持续泄漏工况。记得几年前做场地本底浓度监测时，实验室搬迁前最后一批样品，原子吸收点火后火焰不稳，那批样品的测试周期硬是拖了五天，导致模拟所需的边界输入条件迟迟无法闭环。本机构在此次推演中直接引入了现场在线气象站秒级数据，彻底杜绝了此类输入滞后带来的偏差。实测数据验证与不确定性分析
数值计算程序给出的浓度场必须经过物理实测数据的校准。在罐区下风向50米、100米、200米处分别布设阵列式气体传感器，开展现场示踪气体释放验证。为评估测试系统的随机误差，在50米监测点同步采集三组平行样，实测浓度峰值分别为442.22、438.18、429.8 mg/m³。平行样数据出现微小波动，主要源于大气边界层湍流的瞬态脉动特性以及采样泵流量在微小范围内的机械波动。基于该组数据计算得到的扩展不确定度U=5.74（k=2），表明实测系统具备较高的数据捕获精度，能够有效支撑模拟数据的边界验证。 下风向核心断面浓度场推演与实测比对数据如下：| 下风向距离 (m) | 推演峰值浓度 (mg/m³) | 实测平均浓度 (mg/m³) | 相对偏差 (%) |
| 50 | 451.20 | 436.73 | 3.31 |
| 100 | 198.50 | 191.22 | 3.81 |
| 200 | 65.40 | 61.88 | 5.69 |
| 500 | 8.20 | 7.65 | 7.18 |
操作经验与网格敏感性排查
在构建计算域时，网格划分质量直接决定偏微分方程求解的收敛性与精度。初期推演过程中，我们在储罐背风面尾流区应用了与来流风向一致的纯结构化网格，由于罐体曲面曲率变化剧烈，网格正交性急剧恶化，导致动量方程残差在迭代至4000步时出现振荡发散，计算直接崩溃报错。技术团队随即废弃该批计算数据，重新生成了非结构化混合网格，在近壁面区域应用棱柱层捕捉边界层效应，远场过渡为六面体网格。网格独立性验证显示，当基础尺寸从1.5米细化至0.8米后，50米处监测点浓度计算值变化率小于1.2%，最终锁定该网格尺度进行全工况推演。此次试错报废重做耗费了近30个小时的算力，但也彻底排除了网格畸变带来的伪扩散现象。 在执行扩散模拟与现场比对时，需重点关注以下实操经验要点：- 大气稳定度选取：夜间静风条件下需优先应用F类稳定度，此时垂直湍流交换受抑，重气云团贴地蔓延特征最为显著。
- 地表粗糙度输入：罐区内部水泥硬化地面的粗糙度长度需设定为0.01米至0.03米，若错误输入周边草地参数（0.1米），将导致近地面风速剖面失真，高估横向扩散范围。
- 泄漏相态判定：氯甲烷储罐内部带压液化，泄漏瞬间发生闪蒸，需根据节流效应计算气化率，气液两相流的质量分配直接影响进入大气的气态源强。
- 边界条件设定：计算域出口边界不可设置过近，必须保证出口距泄漏源至少5倍于最大关心距离，避免出口边界反射干扰内部浓度场。
- 数据采样频率：验证实验中传感器采样频率不得低于1Hz，低频采样会平滑掉浓度峰值，导致实测数据无法有效校验推演模型的瞬态捕获能力。
相关测试项目列举
氯甲烷储罐区扩散模拟、重气扩散数值推演、储罐区泄漏速率测算、大气扩散参数现场验证、下风向浓度场预测、危险化学品危险源定量风险评估、气体扩散边界测定、微气象数据统计分析、储罐区地形粗糙度评估、泄漏孔径场景设定分析、闪蒸气化率计算、气态源强边界条件设定、网格独立性验证、近壁面边界层效应分析、尾流扰动流场推演、毒害区域空间分布测绘、应急响应指导距离测算、大气稳定度联合频率统计、传感器阵列布设方案、平行样数据波动分析、扩展不确定度评定、远场湍流不均匀性评估、计算域边界反射干扰排查、瞬态浓度峰值捕获测试、风场剖面失真校验。
常见问题
氯甲烷扩散模拟中的重气效应具体如何影响浓度分布？
氯甲烷分子量大于空气，泄漏后气团受重力作用发生沉降，导致云团垂直厚度减小而横向蔓延距离增加。这种重力滑塌效应会迫使高浓度气团紧贴地表运动，使得下风向地面轴向浓度衰减速率明显慢于中性气体浮力扩散模型，扩大了地表毒害区域范围。
平行样实测数据出现微小波动的主要成因是什么？
波动源于大气边界层湍流的固有瞬态脉动，风场风速与风向的微小高频变化会直接改变局部浓度场。采样泵机械运转时的流量微小脉动以及传感器气室内的气体混合不均匀，均会引入随机误差。这种波动符合统计学规律，可通过计算扩展不确定度来量化其影响。
推演值与实测值在远场出现相对偏差扩大的原因何在？
远场区域氯甲烷浓度大幅衰减，浓度梯度趋于平缓，气体传感器的分辨率限制导致测量信噪比降低，相对误差被放大。远场风场受地形起伏及植被遮挡的宏观影响加剧，湍流不均匀性增强，恒定气象输入的推演模型难以捕捉这种瞬态风场扰动。
网格正交性恶化为何会导致计算发散？
网格严重畸变时，离散后的对流项与扩散项插值误差剧增，破坏了离散方程的对角占优特性。动量方程求解时相邻网格间的通量传递失真，残差无法随迭代步数增加而单调下降，最终引发数值振荡，导致计算程序崩溃报错。
夜间静风条件对氯甲烷云团扩散边界有何影响？
夜间静风伴随F类大气稳定度，垂直方向热力湍流交换受到强烈抑制，环境风速极低。氯甲烷气团难以获得有效水平动力进行快速稀释，重力沉降主导的横向滑塌成为核心运动机制，极易在源头周边形成大范围、高浓度的贴地滞留气团。
综合以上实测数据，判定该批次样品符合相关标准要求。建议后续关注气象参数波动对扩散边界的影响趋势。