1. 课程导入与问题界定

地下水污染羽流（contaminant plume）是指污染物随地下水流动，从污染源向下游迁移扩散所形成的三维污染范围。它不是一个静止的"污染团块"，而是一个随水流方向不断被拉长、稀释、弥散，并可能因吸附、降解、密度差异而发生形态变化的动态体。在污染调查与风险管控中，能否准确刻画羽流的空间范围与浓度分布，直接决定了监测井布设是否合理、修复范围是否被低估、以及健康风险评估结论是否可靠。

理解羽流，需要同时抓住两个维度：空间维度——污染物在水流方向上迁移多远、在横向和垂向上扩散多宽；时间维度——浓度随时间的衰减或稳定趋势。前者依赖对流—弥散方程与水文地质参数的刻画，后者依赖长期监测数据形成的浓度—时间曲线。二者结合，才能回答"污染到哪了、还会到哪、什么时候能达标"这三个调查核心问题。

为什么迁移计算与浓度曲线判读缺一不可

迁移计算提供的是"预测能力"。以最常用的一维对流—弥散方程为例，在稳态连续源条件下，下游某点浓度可近似表达为：

C(x,t) = (C0/2) · erfc[ (x - v·t) / (2·√(D·t)) ]

其中 C0 为源浓度，v 为 seepage 流速（达西流速除以有效孔隙度），D 为水动力弥散系数，erfc 为余误差函数。这个式子告诉我们：羽流前锋位置大致由 v·t 决定，而前锋的"陡缓"由弥散系数 D 控制。若某场地实测 v≈0.05 m/d、D≈0.5 m²/d，则一年后前锋约在 18 m 处，但浓度过渡带可能宽达十余米——这正是仅凭"污染源位置"判断范围会严重出错的根本原因。

浓度曲线判读提供的是"验证与诊断能力"。把同一口井不同时间的浓度连成曲线，可以识别出三类典型形态：持续上升型（羽流尚未到达峰值或源仍在释放）、平台稳定型（源强与衰减达到平衡）、指数衰减型（源已切断或自然衰减主导）。曲线形态的差异，直接对应不同的调查结论与管控策略。若只做计算不读曲线，模型无法校准；若只读曲线不做计算，又无法外推未来趋势。

本课程的任务框架

基于上述认识，本课程围绕一条主线展开：参数 → 计算 → 曲线 → 判读 → 决策。具体包括四个能力层次。第一层是概念与参数层，明确有效孔隙度、水力梯度、弥散度、阻滞因子、衰减常数等参数的地下水文学含义及其获取途径（抽水试验、示踪试验、土柱试验或经验取值）。第二层是计算层，掌握一维、二维对流—弥散解析解的适用条件与手算方法，理解量纲与边界条件对结果的影响。第三层是数据层，学会把监测井浓度序列整理为可判读的曲线，识别异常点、检出限效应与季节性波动。第四层是判读与决策层，把计算结果与实测曲线对照，判断羽流是否稳定、是否到达敏感受体、修复目标是否可达。

需要特别强调参照约定：本课程统一规定地下水流向为 x 正方向，污染源位于 x=0 处，下游监测井 x>0；浓度以质量浓度 mg/L 表示，时间以天（d）为单位。所有公式、图示与案例均遵循这一约定，若某案例采用相反方向，将明确说明。

一个贯穿全课程的示例场景

为使抽象概念落地，本课程设置一个贯穿案例（示例参数，非实测）：某场地浅层含水层为粉细砂，水力坡度 i=0.002，渗透系数 K=5 m/d，有效孔隙度 nₑ=0.15，则 seepage 流速 v=K·i/nₑ≈0.067 m/d。污染源为历史储罐泄漏，源区浓度 C0=50 mg/L，纵向弥散度 αL=5 m。据此可估算一年后羽流前锋约在 24 m 处，而浓度降至 5 mg/L 的位置约在 40 m 之外。后续章节将用这组参数演示完整计算、曲线绘制与判读过程，读者可自行替换为实际场地数据。

最后需要界定本课程不涉及的内容：不展开三维数值模拟软件的建模操作，不讨论污染物化学反应机理的分子层面细节，也不替代场地-specific 的法规符合性判定。课程聚焦于"用可手算、可解释的方法把羽流迁移与浓度曲线讲清楚"，为后续深入数值模拟与风险评估打下判断基础。