1 基本概念
1.1 定义与研究对象
污染扩散模拟是指借助数学方程、计算算法和环境观测数据，对污染物在空气、水体、土壤及地下介质中的迁移、扩散、沉降、稀释和转化过程进行定量分析与预测的方法体系。其研究对象通常包括颗粒物、气态污染物、溶解性化学物质、热污染以及某些具有反应性的复合污染组分。
这类模拟既关注污染物从源头释放后的空间分布，也关注其随时间变化的浓度演化。实际应用中，模型输出常用于判断污染影响范围、估计暴露水平，并为治理、预警和管理提供依据。
1.2 发展背景
污染扩散模拟的发展，与环境监测技术进步、计算机性能提升以及环境管理需求增长密切相关。早期研究多以简化理论为主，强调均匀介质中的扩散规律；随后，随着数值方法和观测数据的丰富，模型逐渐能够描述更复杂的地形、气象、水动力和介质异质性。
在工业化、城市化加速背景下，污染事件的类型更复杂，时空变化更迅速，对定量预测的需求也更强烈。由此，污染扩散模拟从单一学术研究工具，逐步发展为环境评估和应急决策中的常用手段。
1.3 学科归属与应用范围
污染扩散模拟通常位于环境科学、环境工程、流体力学、计算科学和地球信息科学的交叉区域。其核心任务是将物理、化学和生物过程以可计算形式表达，并通过观测数据进行校准与验证。
其应用范围较广，既可用于大气、水体和土壤污染，也可用于热量、盐分或其他保守性物质的输运分析。除环境领域外，在城市规划、工业安全、资源管理和公共健康评估中，也常可见其应用。
2 污染扩散过程
2.1 迁移机制
污染物进入环境后，通常不会停留在单一位置，而是随介质运动发生迁移。迁移过程决定了污染物到达何处、传播多远，以及浓度在空间上的分布格局。
2.1.1 平流作用
平流是指污染物随流体整体运动而被输送的过程。在大气中，平流主要受风场驱动；在水体中，则由水流速度和流向决定。该过程往往决定污染羽的主传播方向，是扩散模拟中的基础项。
2.1.2 扩散作用
扩散来源于分子热运动和湍流混合作用，会使污染物由高浓度区向低浓度区逐渐分散。与平流相比，扩散更强调浓度梯度的消减与污染云团的“摊开”。在湍流强烈或介质混合充分的环境中，这一作用尤为明显。
2.1.3 沉降与沉积
对于颗粒态污染物，重力作用会使其从气相或水相中逐步沉降到地表、河床或底泥中。沉降和沉积不仅影响污染物的空间分布，还会改变后续再悬浮或再释放的可能性，因此常被视为重要的质量交换环节。
2.2 转化机制
污染物在迁移过程中，可能发生化学结构、形态或生物可利用性的变化。此类转化会影响其毒性、持续时间和扩散行为。
2.2.1 化学反应
部分污染物会在环境中发生氧化、还原、水解或光化学反应，从而生成新的化合物。反应速率通常受温度、光照、pH值和反应物浓度影响。对反应性污染物而言，化学转化往往是模型必须考虑的重要过程。
2.2.2 生物降解
在土壤和水体中，微生物可将某些有机污染物分解为更简单的物质。生物降解通常与含氧条件、营养状况和微生物群落有关。该过程在污染自然衰减和生态恢复中具有重要意义。
2.2.3 吸附与解吸
污染物可附着在颗粒物、土壤表面或有机质上，形成吸附；在条件变化时，又可能重新释放，发生解吸。吸附会降低污染物的瞬时迁移能力，而解吸则可能造成滞后性再污染，因此是土壤和沉积环境模拟中的常见参数项。
2.3 影响因素
污染扩散并非仅由源强决定，还受到环境状态和介质性质的共同控制。不同因素之间往往相互耦合，使扩散结果具有明显的时空差异。
2.3.1 风速与流场
风速和流场结构直接影响污染物的输送速度与传播方向。风向的变化会使污染羽发生偏转，湍流强度增大则会加速混合。对于大气扩散来说，风场是决定性因素之一。
2.3.2 温度与湿度
温度会影响气体密度、反应速率和湍流稳定度；湿度则会改变颗粒吸湿、生长和沉降行为。在某些条件下，温度逆温层还可能抑制垂直扩散，使污染物在近地层积聚。
2.3.3 地形与介质特性
地形起伏、建筑物分布、河道形态以及土壤孔隙结构，都会改变污染物的传播路径和混合方式。粗糙度、孔隙率、渗透性等介质属性，则会进一步影响污染的迁移效率与停留时间。
3 模型类型
3.1 解析模型
解析模型通常基于简化假设，将污染扩散过程表示为可直接求解的公式。其优点是结构清晰、计算量小，适合快速估算与理论分析。
3.1.1 高斯烟羽模型
高斯烟羽模型是大气污染扩散中最经典的解析模型之一，常用于连续点源在稳定风场下的浓度估算。该模型假定污染物在横向和垂向上的分布服从高斯函数，适合远场近似和工程快速评估。
3.1.2 一维与二维简化模型
当研究对象具有明显主导流向，或空间变化可在较低维度内近似时，常采用一维或二维模型。这类方法能显著降低计算复杂度，适用于河道输运、浅层扩散和线性污染带分析。
3.2 数值模型
数值模型通过离散化控制方程，对复杂边界和非线性过程进行计算，更适合真实环境中的精细模拟。它们通常需要较多参数与较高计算资源，但适应性也更强。
3.2.1 有限差分法
有限差分法通过网格节点上的差分近似替代连续导数，实现对扩散方程的求解。该方法结构简单，便于编程实现，在规则网格问题中应用广泛。
3.2.2 有限元法
有限元法将研究区域划分为若干单元，在单元上构建近似函数并进行整体组装。它对不规则边界和复杂地形具有较强适应性，适合需要较高几何表达能力的场景。
3.2.3 有限体积法
有限体积法以控制体为基础，强调质量守恒，常用于流体与输运耦合问题。由于其守恒性较好，在大气、水动力与污染输运联立计算中具有较高实用价值。
3.3 统计与经验模型
统计与经验模型主要依赖历史观测数据和经验关系建立污染扩散预测公式，适合数据较多但机理描述有限的情况。此类模型通常计算快捷，便于工程应用。
3.3.1 回归模型
回归模型通过分析污染浓度与气象、排放或地形变量之间的对应关系，建立预测方程。其优点是形式直观，缺点是对样本质量和适用范围较为敏感。
3.3.2 概率扩散模型
概率扩散模型将污染行为视为随机过程，用概率分布描述污染物到达某一区域的可能性。它常用于风险评估和不确定条件下的情景推演，能够反映污染扩散的随机性特征。
3.4 耦合模型
耦合模型将多个子系统联立起来，以更完整地表达污染扩散过程中的相互作用。随着研究对象复杂化，这类模型的重要性不断提高。
3.4.1 大气-地表耦合
大气-地表耦合模型用于分析污染物在空气与地表之间的交换过程，例如干湿沉降、地表再释放和热量交换。它能够更真实地反映近地层污染演化。
3.4.2 水文-水质耦合
水文-水质耦合模型将水流运动与污染物输运、反应过程联合模拟，适用于河流、湖库和地下水系统。该类模型尤其适合研究降雨、径流和污染响应之间的关系。
4 模拟对象
4.1 大气污染扩散
大气污染扩散模拟主要关注污染物在边界层中的传播、稀释与沉降，是应用最广的扩散模拟类型之一。其结果常用于空气质量评估和事故预警。
4.1.1 工业排放
工业排放通常具有连续、稳定或阶段性释放的特点，污染源位置相对明确。模拟时需要考虑烟囱高度、排放速率、废气温度以及周边建筑对流场的扰动。
4.1.2 交通源排放
交通源排放多沿道路或交通走廊分布，具有线源特征。城市道路两侧的浓度梯度、车流密度和街谷效应，往往是此类模拟关注的重点。
4.1.3 事故性泄漏
事故性泄漏具有突发性、释放量大和时变性强等特点，常需要快速预测扩散范围。此类模拟多用于应急响应，强调短时间内的风险判断。
4.2 水体污染扩散
水体污染扩散涉及河流、湖泊、水库和海洋等多种环境，其运动形式受水动力条件显著影响。污染物在水中既可随流输运，也可能发生沉积、再悬浮和反应转化。
4.2.1 河流扩散
河流中的污染扩散通常表现为沿程输运与横向混合并存。流速、河宽、弯道结构和支流汇入，都会影响污染带的形态与到达时间。
4.2.2 湖泊与水库扩散
湖泊与水库水体相对滞缓，分层现象和局部环流对污染扩散影响较大。模拟时常需考虑风致环流、热分层及入湖口附近的局部混合特征。
4.2.3 海洋与近岸扩散
海洋与近岸环境受潮汐、波浪、海流和岸线形态共同控制，扩散过程更为复杂。近岸区域的污染模拟常涉及潮汐往复、海陆交换以及沉积物再悬浮等因素。
4.3 土壤与地下水污染扩散
土壤和地下水污染扩散通常具有隐蔽性强、持续时间长、治理难度大的特点。其模拟重点在于污染物在孔隙介质中的迁移、滞留和衰减。
4.3.1 土层迁移
污染物进入土层后，会沿土壤孔隙水或颗粒表面缓慢迁移。土壤质地、有机质含量和含水状态，会显著影响其运动速度和分布深度。
4.3.2 渗滤与入渗
降雨或地表积液可将污染物带入更深层土壤，这一过程称为渗滤或入渗。它往往决定污染是否进一步进入包气带并威胁地下水系统。
4.3.3 地下水运移
地下水中的污染扩散主要随含水层流动而移动，同时受弥散、吸附和降解控制。由于地下环境难以直接观测，相关模拟往往依赖较细致的参数估计与长期监测数据。
5 数据与参数
5.1 气象与水文数据
模拟污染扩散通常离不开基础环境数据。气象和水文数据决定了介质运动状态，是模型驱动的重要输入。
5.1.1 风场数据
风速、风向、湍流强度和边界层稳定度等风场数据，对大气扩散预测尤为关键。它们不仅影响污染输送速度，也影响垂直混合程度。
5.1.2 降雨与径流数据
降雨量、降雨历时和径流过程会影响污染物的冲刷、迁移与汇集。在地表径流显著的区域，这类数据对水体和土壤污染模拟具有基础意义。
5.1.3 流速与流向数据
流速和流向用于描述河流、湖库和海域中的水动力条件。它们直接决定污染羽的传播路径，并影响混合和稀释效率。
5.2 污染源参数
污染源参数描述污染物从何时、何地、以何种方式释放，是构建模拟场景的核心信息之一。
5.2.1 排放强度
排放强度通常表示单位时间内释放的污染量。该参数直接决定初始污染负荷，并对预测浓度具有显著影响。
5.2.2 排放时间与位置
污染释放的时间和位置决定了污染扩散的起点以及与环境条件的耦合关系。不同时间段的大气稳定度或水文状态，可能造成完全不同的扩散结果。
5.2.3 释放形态与粒径
污染物以气体、液滴、颗粒或溶液形式释放时，其扩散行为会明显不同。颗粒粒径尤其影响沉降、悬浮和再分布过程。
5.3 环境介质参数
环境介质参数反映介质本身对污染迁移的影响，是模型能够逼近真实环境的重要条件。
5.3.1 扩散系数
扩散系数用于表征污染物在介质中的分散能力。其大小受温度、介质结构和湍流条件影响，通常是模型中的关键参数。
5.3.2 孔隙率与渗透性
孔隙率决定介质中可供流体通过的空间比例，渗透性则反映流体通过介质的难易程度。二者在土壤、沉积物和含水层模拟中尤为重要。
5.3.3 湍流与混合系数
湍流与混合系数用于描述不规则流动引起的额外输运效果。它们能够弥补纯分子扩散不足以解释的现实混合现象。
6 模拟方法与流程
6.1 研究区域构建
在正式模拟前，需要对研究区域进行空间表达与结构化处理，以便模型能够读取地形、水系和边界信息。
6.1.1 地形地貌建模
地形地貌建模通过数字高程、岸线、河道或建筑轮廓等信息构建空间基础。对于复杂区域，地形表达越精细，模型结果通常越接近实际。
6.1.2 网格划分
网格划分是数值模拟中的关键步骤，决定计算精度与效率。网格可以是规则或非规则形式，通常需要在污染源附近适当加密。
6.2 边界条件设置
边界条件决定污染物在研究域边缘的行为，是模型稳定运行的重要前提。设置不当往往会导致误差放大或结果失真。
6.2.1 初始条件
初始条件描述模拟开始时污染物和环境状态的分布情况，包括初始浓度、流场和温度等。它是非稳态模拟的起点。
6.2.2 开边界与封闭边界
开边界允许物质和能量进入或离开研究区域，适用于开放流域和大气边界；封闭边界则限制通量穿越，适合封闭容器、局地系统或理想化分析。
6.3 模型求解
模型求解是将控制方程转化为可计算结果的核心环节，通常依赖数值迭代和时间步推进。
6.3.1 稳态求解
稳态求解适用于污染场随时间变化较小或可近似不变的情况。该方法关注长期平衡状态，计算相对简洁。
6.3.2 非稳态求解
非稳态求解用于描述污染扩散随时间连续演化的过程，是更符合实际的模拟方式。它能够反映释放、迁移和衰减的动态变化。
6.3.3 迭代与收敛控制
复杂模型往往需要通过迭代逐步逼近解，并设置收敛准则以判断计算是否稳定。收敛控制直接影响结果可靠性与计算成本。
6.4 结果输出与可视化
模型结果通常以图像、曲线和空间分布图的形式展示，以便于理解扩散特征和识别重点区域。
6.4.1 浓度场展示
浓度场展示可直观反映污染物在空间中的高低分布，是最常见的输出形式之一。它有助于识别污染核心区和影响边界。
6.4.2 时空演化动画
时空演化动画通过连续帧展示污染扩散过程，可清晰呈现污染羽随时间移动、扩张或衰减的过程。它在教学、汇报和应急分析中较为常用。
6.4.3 风险分区图
风险分区图按照浓度阈值或暴露标准，将区域划分为不同等级。该类图件便于与管理措施、预警等级或管控边界对应。
7 模型校验与评估
7.1 实测数据对比
模型是否可信，通常需要通过实测数据进行检验。对比过程可帮助判断模型是否有效再现了真实扩散行为。
7.1.1 现场采样
现场采样是获取污染浓度和环境状态数据的重要方式。采样点布设合理与否，往往会影响校验结论的代表性。
7.1.2 监测站数据
监测站数据具有连续性和规范性，适合用于长期比较和趋势分析。与模型结果对照时，可用于检验时序变化和空间分布一致性。
7.2 精度指标
精度指标用于衡量模拟结果与观测值之间的接近程度，是模型评价的量化依据。
7.2.1 均方误差
均方误差反映预测值与实测值偏离的平均程度，数值越小通常表示拟合越好。它对较大误差较为敏感。
7.2.2 决定系数
决定系数用于衡量模型对观测变化的解释能力。该指标越接近1，通常表示模型对数据趋势的再现能力越强。
7.2.3 偏差分析
偏差分析关注模型是否系统性高估或低估污染浓度。它有助于发现模型结构、参数设定或边界条件中的潜在问题。
7.3 不确定性分析
不确定性分析用于识别影响结果稳定性的来源，并评估模型输出的可信范围。对于复杂环境系统，这一步尤为重要。
7.3.1 参数敏感性
参数敏感性分析用于判断哪些参数对结果影响最大。通过该方法，可以优先改进关键参数的测定或估计精度。
7.3.2 误差传播
误差传播分析关注输入数据误差如何在计算过程中放大或传递到输出结果。它能够帮助理解不确定性的来源与累积方式。
7.3.3 场景假设检验
场景假设检验通过比较不同边界条件、排放情景或气象条件下的结果，评估模型的稳健性。此类检验常用于应急预案和方案比选。
8 应用领域
8.1 环境影响评价
在环境影响评价中，污染扩散模拟常用于预测建设项目或工程活动对周边环境的潜在影响。它能够帮助判断是否会出现超标区域，并为缓解措施设计提供依据。
8.2 事故应急响应
当污染事件突然发生时，扩散模拟可迅速提供污染范围、传播方向和到达时间等信息，是应急决策的重要支撑工具。
8.2.1 泄漏扩散预警
泄漏扩散预警侧重于在事故初期判断污染是否会影响敏感区域。通过快速模拟，可为提前通报和处置争取时间。
8.2.2 疏散与管控辅助
在需要临时管控的场景中，模拟结果可为人员疏散路线、封控边界和防护范围提供参考。其目标是尽可能降低暴露风险。
8.3 污染治理方案设计
污染扩散模拟可以比较不同治理措施的效果，例如截流、稀释、覆盖、通风或净化方案。通过情景对比，能够辅助选择更高效的治理路径。
8.4 风险评估与预警系统
在风险评估与预警系统中，模拟结果常与阈值规则、监测数据和预报信息结合使用。这样可实现对潜在污染事件的提前识别和分级响应。
8.5 城市与区域环境管理
在城市和区域尺度上，污染扩散模拟可用于优化空间布局、道路管理、排放控制和环境容量评估。它也能帮助管理者识别重点治理区和长期累积风险区。
9 典型软件与工具
9.1 通用仿真平台
通用仿真平台通常集成流体、传输和反应模块，便于开展多场耦合分析。此类平台往往适合科研和工程应用的统一建模。
9.1.1 环境扩散模块
环境扩散模块主要负责污染物输运、稀释和沉降计算，通常内置多种扩散方程和源项形式。它是污染模拟软件中的核心功能之一。
9.1.2 流体计算模块
流体计算模块用于求解风场、水流和湍流结构，为污染扩散提供动力背景。没有可靠的流体场，扩散结果往往难以准确。
9.2 地理信息系统集成
地理信息系统为污染扩散模拟提供空间数据管理、地图表达和区域分析能力。它常与模型联用，以提高结果的空间解释性。
9.2.1 空间数据处理
空间数据处理包括坐标转换、栅格化、矢量叠加和地形提取等操作。它能将不同来源的数据统一到同一空间框架中。
9.2.2 结果叠加分析
结果叠加分析可将污染浓度图与人口分布、土地利用或敏感目标位置进行组合，便于识别重点防护区域。该方法在风险管理中十分实用。
9.3 开源与商业工具
污染扩散模拟工具既有开放源代码框架，也有成熟的商业软件。不同工具在功能完整性、可扩展性、使用门槛和技术支持方面各有侧重。
9.3.1 开源建模框架
开源建模框架便于定制算法、修改参数和二次开发，适合研究型用户。其优势在于透明度高，但通常需要较强技术基础。
9.3.2 商业数值软件
商业数值软件通常提供较完善的界面、求解器和后处理能力，适合工程项目中的快速部署。它们在稳定性和支持服务方面往往更成熟。
10 发展趋势
10.1 多尺度模拟
未来的污染扩散模拟将更加重视从局地到区域、从瞬时到长期的多尺度统一描述。这样既能保留细节，又能兼顾整体过程。
10.2 高分辨率与实时预测
随着监测网络和计算能力提升，模型正向更高时空分辨率与更快更新频率发展。实时预测能力将使污染预警和应急响应更具时效性。
10.3 人工智能辅助建模
人工智能可用于参数反演、缺测补全、模型校准和快速代理计算。它并不替代机理模型，而是常作为辅助工具提升建模效率。
10.4 传感网络与数字孪生
传感网络提供连续观测数据，数字孪生则将真实环境与虚拟模型同步更新。二者结合后，可实现更动态的污染监测、模拟与决策支持。