基于蒙特卡洛模型的工业颗粒物风致PM排放精细化与定量评估

《Atmosphere》：Refinement and Quantitative Evaluation of a Monte Carlo Model for Wind-Driven PM Emissions from Industrial Granular Materials

　　

论文摘要翻译 先前已开发了一个物理-数学模型，用于估算暴露于风蚀的颗粒表面产生的粉尘排放。该模型基于控制风致粉尘排放的主要物理机制，其中细颗粒物的释放由盐跳（saltation）及其伴随的喷砂（sandblasting）过程控制。采用概率性蒙特卡洛（Monte Carlo）方法来模拟盐跳体对可蚀表面的撞击，并根据每次碰撞过程中释放的基本粒子质量估算颗粒物（PM）排放。虽然最初的研究仅提供了定性评估，但本研究引入了计算改进，并从数值性能以及预测的PM排放行为的物理一致性两方面，对该模型进行了首次定量评估。算法通过引入固定数量的模拟撞击次数进行了修改，从而降低了计算成本。修订后的模型被应用于意大利撒丁岛一家工业工厂的铅和锌硫化物精矿，以评估所提出模型改进的效果。评估重点关注：(i) 模拟排放量对模拟撞击次数的敏感性；(ii) 修订后模型再现PM排放量级、喷砂效率及其对风摩擦速度依赖性的能力。结果表明，将模拟撞击次数从10,000次减少到500次，与所研究的最高采样配置相比，模拟PM排放量的中位数相对差异为3.35%，同时计算时间减少了约95%。此外，该模型复现的排放量级和关键排放参数，总体上与文献中针对具有相似物理性质的材料所报道的结果一致。总体而言，修订后的模型为估算工业颗粒材料的PM排放提供了一个高效、基于物理的工具。

论文解读文章

1. 研究背景与问题

在采矿、矿物加工和冶金设施中，大量颗粒材料以料堆或尾矿库的形式储存。这些表面受到风蚀会产生显著的颗粒物（Particulate Matter, PM）及有毒化合物排放，对空气质量产生不利影响，尤其在干旱多风地区，可能导致周边区域空气质量标准超标，对人类健康和生物多样性构成严重威胁。欧盟《关键原材料法案》也强调了减少采矿废物和临时储存区环境影响的必要性。因此，可靠的风蚀PM排放估算工具对于大气影响评估和制定有效缓解策略至关重要。

然而，精确估算面临挑战，因为控制风蚀的主要变量（如粒径分布、含水量、颗粒间内聚力和矿物组成）之间存在高度非线性的相互作用。现有方法存在局限性：美国环保署等机构开发的标准化方法虽简便，但忽略了场地特异性因素及时空变异性，可能导致不准确的估算；而使用环境风洞等实验设备进行现场测量则需要大量时间、资源和实验工作，且测量结果本身也存在不确定性。因此，亟需一种既能真实描述风蚀和PM排放过程，又具备足够计算效率的建模工具。为此，Dentoni等人曾开发了一个基于概率性蒙特卡洛方法的物理-数学模型，但其高计算成本和缺乏系统性数值性能评估限制了其广泛应用。

本研究旨在巩固该前期模型，通过提高计算效率和系统评估其数值可靠性，提供首次定量评价。研究聚焦于两个具体目标：一是展示为简化模型使用和降低计算成本而引入的更新，并批判性地评估其效果和量化相关不确定性；二是通过与实验观测和文献数据对比，评估模型预测的物理一致性。

2. 关键技术方法

研究人员采用了以下主要技术方法：

1. 1.
   物理-数学模型：基于概率性蒙特卡洛（Monte Carlo）方法，模拟盐跳体（saltator）对颗粒表面的撞击，并通过两种替代的“损伤比”（Damage Ratio, DR）公式（Thornton等人提出的公式以及Deng和Davé提出的公式）来量化撞击导致的团聚体破碎及细颗粒释放。
2. 2.
   实验表征：对来自意大利撒丁岛某冶金厂的铅硫化物（PbS）和锌硫化物（ZnS）精矿进行物理特性测定，包括利用激光衍射法测定颗粒粒径分布（Particle Size Distribution, PSD）、干筛法测定团聚体粒径分布（Aggregate Size Distribution, ASD）、以及在DICAAR环境风洞（Environmental Wind Tunnel, EWT）中测定阈值摩擦速度（threshold friction velocity, u_*t）。
3. 3.
   数值模拟与性能评估：应用改进后的模型，对不同风摩擦速度、不同DR公式和不同模拟撞击次数（n = 500至10,000）的组合进行仿真。通过比较计算时间、相对标准偏差（Relative Standard Deviation, RSD）和95%置信区间半宽等指标，评估模型的数值性能和统计稳定性。

3. 研究结果

4.1. 模拟撞击次数对模型数值性能的影响

研究表明，模拟撞击次数主要影响蒙特卡洛解的统计稳健性，而非预测的PM排放值本身。将撞击次数从10,000次降至500次，计算时间大幅减少（例如，使用Thornton DR公式时，硫化锌精矿平均模拟时间从9.01小时降至0.44小时）。同时，统计变异性随撞击次数增加而系统性降低，例如，硫化锌精矿的平均相对标准偏差（RSD）在使用Thornton公式时从10.45%降至2.09%。最低（500次）和最高（10,000次）采样水平的直接比较显示，中位数相对差异仅为3.35%，表明少量撞击即可获得可靠的均值估计。

4.2.1. PM排放量

模型预测的PM10排放范围为0.04至14.20 mg·m^-2·s^-1，PM2.5排放范围为0.01至2.71 mg·m^-2·s^-1。两种DR公式均显示相同趋势：在较低摩擦速度下，硫化锌精矿排放更高；在高摩擦速度下，铅硫化物精矿排放更高。两种DR公式预测的排放量存在显著差异（约1-2个数量级），这直接源于它们预测的平均DR值相差约两个数量级。预测的PM2.5/PM10排放比（约0.18-0.19）基本恒定，与美国环保署推荐值一致。

4.2.2. PM排放对风摩擦速度的依赖性

对于两种材料，PM排放随摩擦速度的增加均可用幂律关系（F_D = a·u*^b）精确描述（R² = 0.98-0.99）。铅硫化物精矿的指数b（约6.7）高于硫化锌精矿（约5.1），表明其对摩擦速度增加的敏感性更强。

4.2.3. 喷砂效率

预测的喷砂效率（α）在所研究的摩擦速度范围内基本与风速无关。Deng和Davé公式预测的α值（约10^-2 m^-1）与文献上限范围一致，而Thornton公式预测的值（约10^-1 m^-1）则高出约一个数量级。铅硫化物精矿的α值约为硫化锌精矿的1.5倍，这与前者更低的表面内聚能相一致。

4. 讨论与结论总结

讨论部分总结：

讨论指出，模拟撞击次数主要影响结果的统计稳定性，而非均值本身，研究者可根据需求在计算成本和稳健性之间权衡。两种DR公式导致排放量级存在巨大差异，这是模型形式不确定性的主要来源，因为这两个公式并非专为风致盐跳条件开发，无法确定孰优孰劣。预测的PM排放量级与文献中高侵蚀性材料（如粉砂/黏土丰富的土壤、尾矿）的数据一致。预测的PM2.5/PM10比值恒定，源于模型中假设盐跳体冲击速度与风速无关。预测的幂律指数b（5.1-6.7）高于多数文献值，这可能是因为模型假设了未耗竭的表面，而实验中表面会逐渐耗竭，从而抑制了排放随风速的增长。喷砂效率α与风速无关，同样源于冲击速度恒定的假设。模型的主要局限包括：DR公式的不确定性、假设无限且不变的侵蚀物质可用性（未考虑表面耗竭），以及仅限于干燥材料。

研究结论部分翻译：

本研究描述了用于模拟工业颗粒材料风蚀产生PM排放的蒙特卡洛数值模型的更新与应用。提出的改进主要涉及两个方面：盐跳体粒径分布的自动确定，这减少了所需输入数据量并简化了模型使用；以及撞击模拟过程的优化，这大幅降低了计算成本。计算性能分析表明，模拟撞击次数可以大幅减少，而对估计的平均PM排放量影响有限。收敛分析进一步表明，随着模拟撞击次数的增加，模拟的统计不确定性和连续配置之间的差异均系统性减小，这表明采用的10,000次撞击配置提供了令人满意的收敛数值解近似。这一改进将该模型的潜在应用扩展到需要大量模拟的参数研究、敏感性分析和其他应用中。该模型再现了风致PM排放的主要特征。模拟排放的量级、其对风摩擦速度的依赖性以及PM2.5/PM10比值，总体上与文献中针对具有相似特性的材料所报告的值一致。预测的喷砂效率对采用的DR公式表现出更强的依赖性，使用Deng和Davé公式获得的值大致与文献上限范围一致，而Thornton公式则系统性产生更高的值。总体而言，结果支持了所提出的建模框架用于研究材料特性对PM排放影响的能力，同时也突显了采用的DR公式对预测排放水平的重要性。从本研究中得出了未来模型发展的两个优先事项：第一是开发专门设计并经实验校准的DR公式，以描述盐跳过程中颗粒团聚体的破碎；第二是引入侵蚀表面的演化表示，以考虑风蚀过程中易蚀物质的逐步耗竭。解决这些问题将改善风蚀过程的物理表示，并将模型的适用性扩展到更广泛的环境条件。