3. 对流弥散方程（ADE）的建立与参数含义

要判读一条地下水污染羽的浓度曲线，先得知道这条曲线是被什么方程“塑造”出来的。对流弥散方程（Advection–Dispersion Equation，ADE）就是描述溶质在地下水中迁移的基本方程。它并不神秘，本质上是“质量守恒”在含水层中的具体写法：一个微小控制体内溶质质量的变化率，等于流入与流出之差，再加上源汇项与反应项。

从质量守恒到一维 ADE

取沿地下水流方向的一维控制体，厚度为 Δx、单位面积。溶质通量由两部分组成：一是随水流整体搬运的对流通量 q_c = n·v·C（n 为有效孔隙度，v 为孔隙流速，C 为浓度）；二是由浓度梯度驱动的弥散通量，按 Fick 定律写作 q_D = −n·D·∂C/∂x，D 为水动力弥散系数。对通量做空间差分、再除以孔隙体积，并计入吸附阻滞与一阶衰减，得到一维 ADE：

R·∂C/∂t = D·∂²C/∂x² − v·∂C/∂x − λ·C

其中 R 为阻滞因子，λ 为一阶衰减系数。若写成通量守恒形式，则左侧为 ∂(nC)/∂t，右侧为 −∂(nvc)/∂x + ∂(nD∂C/∂x)/∂x 再减去反应项，两种写法在 n 为常数时等价。这里约定 x 正方向与地下水流向一致，v 取正值；若反向取坐标，则对流项符号随之改变，全文与图示必须统一。

推广到二维：把方向性写进弥散张量

实际羽流不会只沿一个方向铺开。把控制体扩展到 x–y 平面，对流通量变为矢量 (v_x, v_y)，弥散则要用二阶张量 D 描述，因为纵向弥散度 α_L 与横向弥散度 α_T 不同，且张量主轴随水流方向旋转。二维 ADE 写作：

R·∂C/∂t = ∂/∂x(D_xx·∂C/∂x + D_xy·∂C/∂y)
         + ∂/∂y(D_yx·∂C/∂x + D_yy·∂C/∂y)
         − v_x·∂C/∂x − v_y·∂C/∂y − λ·C

当水流沿 x 轴时，张量退化为对角形式：D_xx = α_L·v，D_yy = α_T·v，交叉项为零。这一步很关键——它说明弥散系数不是独立参数，而是“弥散度 × 流速”的乘积，流速越大，同样的弥散度会产生更强的机械弥散。

四个参数各自“管什么”

孔隙流速 v 决定羽流质心的推进速度。在忽略阻滞与衰减时，浓度峰值大致以 v 前进，因此“峰值到达时间 ≈ 距离 / v”是判读曲线时最直接的估算。注意 v 是孔隙流速，与达西流速 q 之间差一个有效孔隙度：v = q/n，n 通常取 0.1～0.35，这个换算在实测中极易出错。

弥散度 α_L、α_T 决定曲线的“胖瘦”。纵向弥散度控制沿流向的展宽，横向弥散度控制侧向铺开。工程经验上 α_L 常取运移距离的 1/10 量级（如运移 100 m 时 α_L 约 10 m），α_T 通常为 α_L 的 1/10～1/3，但这些只是量级参考，必须以场地示踪或浓度拟合为准。弥散越强，峰值越低、拖尾越长，曲线越“平缓”。

阻滞因子 R 来自吸附，R = 1 + (ρ_b/n)·K_d，ρ_b 为干容重，K_d 为分配系数。R > 1 时溶质被固相“拖住”，曲线整体右移，到达时间约为无吸附时的 R 倍，但峰形宽度基本不变。R = 1 表示保守溶质（如 Cl⁻、Br⁻），常被用作示踪剂来反推流速与弥散度。

衰减系数 λ 描述一阶降解或放射性衰变，半衰期 t₁/₂ = ln2/λ。它使整条曲线按 exp(−λt) 整体压低，而不改变峰的位置与宽度。若观察到曲线“右移且变矮”，往往是 R 与 λ 共同作用；若“只变矮不移位”，则更可能是 λ 主导。

一个可复算的教学示例

设某砂层 n = 0.30，达西流速 q = 0.15 m/d，则 v = 0.5 m/d。取 α_L = 5 m，则 D = α_L·v = 2.5 m²/d。若污染物为保守溶质（R = 1、λ = 0），在 x = 50 m 处，峰值到达时间约 50/0.5 = 100 d。若该溶质 K_d 使 R = 3，则峰值推迟到约 300 d；若同时 λ = 0.005 /d，则 300 d 时浓度已衰减为初始的 exp(−1.5) ≈ 0.22 倍。这组数字说明：判读曲线时先看峰位定 v 与 R，再看峰宽定 α_L，最后看整体高度定 λ，三者可以相对独立地反推。