4. 一维对流弥散解析解与浓度曲线形态

在一维均匀地下水流场中，溶质沿水流方向（取 x 轴正方向）迁移，其浓度 C(x,t) 满足一维对流—弥散方程。若假定含水层均质、各向同性、流速 v 与弥散系数 D 为常数、溶质不发生吸附与降解，则方程可写为：

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

其中 D = α·v + D*，α 为纵向弥散度，D* 为有效分子扩散系数（通常远小于机械弥散项，工程计算中常忽略）。该方程的两个解析解分别对应两种典型源条件：连续注入（阶跃源）与瞬时投放（脉冲源）。理解这两个解的曲线形态，是判读实际监测浓度曲线的基础。

连续注入：阶跃源解与 S 形曲线

设 t=0 时刻起，在 x=0 处持续注入浓度为 C₀ 的污染水，边界条件为 C(0,t)=C₀（t>0），初始条件 C(x,0)=0。其解析解为：

C(x,t) = (C₀/2) · [ erfc( (x − v·t) / (2·√(D·t)) ) 
                    + exp(v·x/D) · erfc( (x + v·t) / (2·√(D·t)) ) ]

当佩克莱数 Pe = v·x/D 较大（即对流占优、距离较远）时，第二项可忽略，解简化为：

C(x,t) ≈ (C₀/2) · erfc( (x − v·t) / (2·√(D·t)) )

该曲线在空间上呈典型的 S 形（sigmoid）：在锋面位置 x ≈ v·t 附近，浓度由接近 C₀ 快速下降到接近 0。曲线的拐点恰好位于 C = C₀/2 处，即 x = v·t 的位置；拐点处斜率绝对值最大，对应浓度梯度最陡。拐点两侧的展宽程度由弥散项 2√(D·t) 控制：D 越大或时间越长，过渡带越宽，S 形越平缓。工程上常把 C/C₀ = 0.5 的位置视为“锋面中心”，把 C/C₀ 从 0.84 降到 0.16 的区间（对应 erfc 自变量 ±1）近似作为过渡带宽度，约为 4√(D·t)。

瞬时投放：脉冲源解与钟形曲线

设 t=0 时刻在 x=0 处一次性投放单位面积质量为 M 的溶质（面源强 M，量纲为 M/L²），初始条件为狄拉克函数 C(x,0) = (M/n)·δ(x)（n 为有效孔隙度）。其解析解为：

C(x,t) = (M / (n · √(4π·D·t))) · exp( −(x − v·t)² / (4·D·t) )

该解在空间上是一条高斯型钟形曲线：峰值位置位于 x = v·t，峰值浓度为 C_max = M / (n·√(4π·D·t))。峰值随时间按 1/√t 衰减，同时曲线不断展宽——这是质量守恒的必然结果：总质量 M 不变，展宽必然伴随峰值降低。曲线的标准差为 σ = √(2·D·t)，常用“峰值 ±σ”或“半高全宽”描述污染羽的展布范围。峰值两侧各有一个拐点，位于 x = v·t ± √(2·D·t) 处，拐点处浓度约为峰值的 0.607 倍，是浓度梯度变化最快的位置。

两种源条件的曲线对比与识别

把两种解放在同一坐标系下比较，可以提炼出实用的识别特征：连续注入曲线是“单侧阶跃”，浓度从源端向远端单调衰减，拐点位于 C₀/2 且随时间右移但形态趋于稳定；瞬时投放曲线是“双侧钟形”，存在明确的峰值，峰值位置随时间线性右移（x=v·t），峰值高度随时间按 1/√t 衰减。实际监测中，若某断面浓度长期稳定在某一水平并向下游递减，多为连续源；若观测到浓度先升后降、峰值逐期向下游推移且峰值降低，则更符合瞬时投放或有限体积源的特征。

示例计算：用参数判断曲线形态

设某均质砂层中 v = 0.5 m/d，纵向弥散度 α = 2 m，忽略分子扩散，则 D = α·v = 1.0 m²/d。考察 t = 100 d 时两种源条件的曲线特征（以下为教学示例参数，非实测数据）：

• 连续注入：锋面中心位于 x = v·t = 50 m，拐点浓度 C₀/2；过渡带宽度约 4√(D·t) = 4×√100 = 40 m，即从约 30 m 到 70 m 浓度由接近 C₀ 降至接近 0。
• 瞬时投放：峰值位于 x = 50 m，标准差 σ = √(2·D·t) = √200 ≈ 14.1 m，拐点位于 50 ± 14.1 m，即约 35.9 m 与 64.1 m 处；若取 M/n = 100 g/m²，则峰值浓度 C_max = 100/√(4π×100) ≈ 2.82 g/L（按一维面源公式估算）。

把 t 改为 400 d，连续注入的锋面中心移到 200 m，过渡带宽度增至约 80 m；瞬时投放的峰值同样移到 200 m，但峰值浓度降为约 1.41 g/L，σ 增至约 28.3 m。可见时间加倍时，连续源的过渡带按 √t 展宽，而脉冲源的峰值按 1/√t 衰减、宽度按 √t 增长，两者在“展宽”上规律一致，差别在于是否存在稳定平台与独立峰值。

曲线判读的实操要点

拿到一组不同时刻的浓度—距离监测数据后，可按以下顺序判读：先看曲线是单调下降还是存在峰值，区分连续源与脉冲源；再看过渡带或峰宽随时间的变化，若宽度按 √t 增长，说明弥散项主导、参数 D 估计合理；然后定位 C/C₀ = 0.5 或峰值位置，反算视流速 v ≈ x/t，与实测水位资料对比校验；最后检查曲线是否对称——脉冲源解在无反应条件下应近似对称，若明显拖尾，往往提示存在吸附、滞后或非均质影响，需在模型中补充相应过程。通过“形态—位置—宽度”三个特征与解析解对照，可以在缺少完整数值模拟时快速判断污染羽的迁移阶段与主控过程。