2. 污染羽流迁移的物理过程解析

地下水中的污染物一旦进入含水层，其空间分布并不是静止不变的，而是被四类基本机制共同驱动和改造：对流决定羽流整体往哪里走、走多快；弥散（含机械弥散与分子扩散）决定羽流边缘如何展宽、峰值如何降低；吸附决定污染物有多少被固相滞留、迁移被延迟多少；降解决定污染物总量如何随时间衰减。理解这四者的叠加关系，是后续把实测浓度曲线"翻译"成物理过程的前提。

一、对流：羽流的"骨架"与时间尺度

对流由地下水达西流速驱动，污染物随水流整体平移。工程上常用渗流速度（达西流速）与实际孔隙流速两个量：

v = K · i / n_e        （实际孔隙流速，m/d）
K   — 渗透系数（m/d）
i   — 水力梯度（无量纲）
n_e — 有效孔隙度（无量纲）

示例（合理假设）：某砂层 K = 8 m/d，i = 0.002，n_e = 0.25，则 v = 8×0.002/0.25 = 0.064 m/d ≈ 23 m/a。若污染源距下游监测井 200 m，仅考虑对流时污染物约 8.7 年到达该井。这个"到达时间"就是浓度曲线开始抬升的时间尺度，也是判断曲线首现时间是否合理的基准。

对流本身不改变羽流的形状，只平移它。若实测曲线首现时间远早于对流估算值，应怀疑存在优先流通道（裂隙、透镜体、废弃井）；若远晚于估算值，则要考虑吸附延迟或源区释放受限。

二、弥散：边缘展宽与峰值衰减

弥散包括机械弥散（孔隙尺度流速差异与路径弯曲引起）和分子扩散（浓度梯度驱动，通常远小于机械弥散）。一维条件下常用弥散系数描述：

D = α · v + D*        （α 为弥散度，m；D* 为有效扩散系数）

弥散度 α 是尺度相关的经验参数，砂性介质常见量级为 0.1～10 m，粘性土更小。弥散的作用在浓度曲线上表现为：锋面不再是阶跃，而是呈 S 形过渡；峰值随距离降低；拖尾逐渐拉长。纵向弥散控制沿流向的展宽，横向弥散控制羽流宽度——横向弥散度通常只有纵向的 1/10～1/3，这正是污染羽往往"细长"的原因。

三、吸附：迁移的"延迟因子"

吸附使部分污染物附着于固相，表现为迁移速度慢于水流。对线性等温吸附，延迟因子为：

R = 1 + (ρ_b / n_e) · K_d
ρ_b — 干容重（kg/L 或 g/cm³）
K_d — 分配系数（L/kg）

示例：ρ_b = 1.6 g/cm³，n_e = 0.25，K_d = 0.5 L/kg，则 R = 1 + (1.6/0.25)×0.5 = 4.2。污染物实际迁移速度约为孔隙流速的 1/4.2，前述 8.7 年的到达时间将延长至约 36 年。吸附不改变稳态羽流的总量关系，但显著推迟曲线抬升时刻，并使穿透曲线更平缓。对非线性（Freundlich/Langmuir）吸附，高浓度段延迟小、低浓度拖尾段延迟大，曲线会出现明显的不对称拖尾。

四、降解：总量衰减与稳态长度

降解（生物降解、水解、放射性衰变等）使污染物浓度按一级动力学衰减：

C(t) = C₀ · exp(−λt)，  λ = ln2 / t₁/₂

降解与对流、弥散耦合后，羽流存在一个稳态长度：当对流输运与衰减达到平衡，羽流不再向前延伸。一维稳态近似下特征衰减长度为 L = v/(R·λ)。示例：v = 0.064 m/d，R = 4.2，t₁/₂ = 2 年（λ ≈ 9.5×10⁻⁴ /d），则 L ≈ 0.064/(4.2×9.5×10⁻⁴) ≈ 16 m。这意味着即使源区持续释放，污染羽也大致稳定在源下游数十米范围内，浓度沿程呈指数下降——这与无降解时"持续拉长"的形态有本质区别。

五、四机制叠加与浓度曲线的对应关系

把四类机制放到同一张"过程—曲线"图上，可以建立如下判读逻辑：对流定首现时间与整体走向；弥散定过渡带宽度与峰值高低；吸附定延迟倍数与拖尾形态；降解定稳态长度与沿程衰减斜率。下图以单点监测井的穿透曲线为主线，展示各机制对曲线不同部位的控制作用。

需要强调的是，上述参数均为教学示例中的合理假设值，实际场地必须通过抽水试验、示踪试验、批试验或土柱试验获取 K、n_e、α、K_d 与降解速率，不能直接套用。此外，四类机制并非简单线性叠加：吸附会降低有效迁移速度，从而同时改变弥散的时间尺度和降解的暴露时间；降解会削减高浓度段，使非线性吸附的拖尾特征更明显。因此在实际判读中，应先用首现时间与峰值位置约束对流与吸附，再用过渡带宽度约束弥散度，最后用曲线回落段与沿程浓度斜率约束降解速率，形成"由粗到细"的参数识别顺序。