2. 污染羽迁移的驱动机制

要判断污染羽会往哪里走、走多远、浓度如何变化，必须把地下水中的污染物迁移拆解为几个同时发生、相互叠加的物理与化学过程。对流决定“质心往哪走”，弥散与扩散决定“羽状体如何摊开变淡”，吸附决定“污染物是否被滞后”，降解决定“总量是否被削减”。本节围绕这五种机制，说明它们各自在迁移趋势中的作用，以及如何用可计算的参数把它们组合起来判断污染羽的走向。

一、对流：污染羽质心的搬运者

对流（advection）是污染物随地下水整体流动而迁移的过程，是污染羽迁移的主导机制。其核心参数是达西流速与孔隙流速的换算：

达西流速   q = K · i
孔隙流速   v = q / n_e = K · i / n_e
对流迁移距离  L = v · t

其中 K 为渗透系数（m/d），i 为水力梯度（无量纲），n_e 为有效孔隙度。举例说明（示例参数，非实测）：某砂层 K = 5 m/d，i = 0.002，n_e = 0.25，则 q = 0.01 m/d，v = 0.04 m/d。若污染已持续 10 年（约 3650 天），仅考虑对流时污染质心约下移 146 m。这个量级说明：对流参数直接决定污染羽的“前锋位置”和受体是否在到达时间窗口内被暴露。实际工作中，K 的空间非均质性往往比对流计算本身更关键，因此应优先用钻孔抽水试验或压水试验分层获取 K，而不是全场地取一个平均值。

二、弥散与扩散：让羽状体摊开、峰值降低

机械弥散由孔隙尺度流速差异与流线弯曲引起，沿水流方向（纵向弥散度 α_L）和垂直方向（横向弥散度 α_T）表现不同；分子扩散由浓度梯度驱动，在低速黏性土层中不可忽略。两者合称水动力弥散，常用弥散度描述：

D_L = α_L · v + D*      D_T = α_T · v + D*

弥散的作用不是搬运质心，而是使浓度峰变平、羽状体边缘变宽。工程经验上 α_L 常取迁移距离的 1/10～1/100 量级，α_T 通常为 α_L 的 1/10～1/3；这些是经验取值范围，具体场地应通过示踪试验或浓度剖面反演校核。弥散对“受体暴露识别”的意义在于：即使质心尚未到达取水井，弥散前锋也可能使井中浓度提前出现可检出水平的上升，因此判断暴露不能只看质心距离。

三、吸附：让污染物“走慢一点”

吸附使部分污染物暂时附着于固相，表现为迁移滞后。对线性等温吸附，用阻滞因子 R 描述：

R = 1 + (ρ_b / n_e) · K_d
实际迁移速度  v_c = v / R

其中 ρ_b 为干密度（g/cm³），K_d 为分配系数（mL/g）。仍用上例：若 ρ_b = 1.6 g/cm³、n_e = 0.25、K_d = 0.5 mL/g，则 R = 1 + (1.6/0.25)×0.5 = 4.2，污染物的实际迁移速度降为对流的约 1/4.2，10 年质心距离由 146 m 缩短到约 35 m。需要注意：R 只适用于线性、可逆、平衡吸附；对重金属等可能出现非线性或竞争吸附的污染物，应改用 Freundlich/Langmuir 等温线或数值模型，不能简单套用单一 R 值。

四、降解：削减总量，改变趋势方向

降解（含生物降解、水解、氧化还原转化）使污染物总量减少，常用一阶衰减描述：

C(t) = C₀ · exp(−λ · t)
半衰期  t₁/₂ = ln2 / λ ≈ 0.693 / λ

例如 λ = 0.0005 /d 时，半衰期约 1386 天（约 3.8 年），10 年后剩余浓度约为初始的 16%。降解与对流、弥散叠加后，污染羽可能形成“稳态长度”——前锋推进与衰减削减达到平衡，羽状体不再明显扩展。这正是判断污染羽“是否仍在扩展”的关键：若监测显示羽状体长度多年基本稳定，往往说明衰减与迁移已接近平衡，而非污染源已消失。

五、机制叠加：如何判断迁移趋势

五种机制并非独立，而是共同作用于同一条迁移路径。下图给出各机制在污染羽迁移链中的位置与作用方向，可作为分析迁移趋势时的对照框架。

把上述机制组合起来，可以形成一条实用的判断链：先用 K、i、n_e 估算对流前锋到达受体的时间；再用 α_L、α_T 判断弥散前锋是否提前造成可检出浓度；用 R 修正实际迁移速度；最后用 λ 判断浓度是否在到达前已被显著削减。若对流前锋到达时间远大于受体可接受的时间窗口，且降解半衰期较短，则暴露风险相对可控；反之，即使质心尚远，也应警惕弥散前锋与优先流通道造成的提前暴露。

综上，对流回答“往哪走、多快到”，弥散与扩散回答“摊多宽、峰多高”，吸附回答“慢多少”，降解回答“少多少”。只有把这四类作用放在同一时间尺度上比较，才能对污染羽迁移趋势和受体暴露识别作出可靠判断。