2. 污染羽概念模型构建
污染羽概念模型是对“污染物从哪里来、何时来、经过什么介质、以什么方式扩散、最终形成什么形状”这一完整链条的结构化表达。它不是一张漂亮的三维示意图，而是后续数值模型、监测方案和修复决策的共同基础。概念模型错了，后面所有计算都只是在精确地犯错。本节围绕四个相互咬合的环节展开：污染源识别、释放历史反演、水文地质概念模型要素，以及污染羽形态的控制因素。
一、污染源识别：先回答“源在哪里、是什么”
污染源识别要同时确定源的位置、类型、强度与空间尺度。工程上通常分三步推进。第一步是资料溯源：查阅企业生产工艺、原辅材料清单、储罐与管线台账、历史航拍与厂区平面图变更记录，锁定潜在源区。第二步是现场筛查：在疑似源区布设土壤气、表层土和浅层地下水快速筛查点，用便携式光离子化检测器（PID）、X射线荧光光谱（XRF）或现场试剂盒判断异常中心。第三步是确认性采样：对异常点送实验室分析特征污染物，并结合地层岩性判断污染物是滞留于包气带还是已进入饱和带。
源的类型直接决定后续建模方式。以某电镀厂为例（示例假设）：历史槽液泄漏形成的源可近似为面状持续源，源区约 40 m×25 m，特征污染物为六价铬与镍；而废液罐区泄漏则更接近点状源，泄漏点集中、初始浓度高。若源区存在非水相液体（DNAPL/NAPL），还必须判断其是否已到达弱透水层顶板并形成池状聚集，因为这决定了源区是“可清除的有限质量”还是“长期缓慢释放的持久源”。
二、释放历史反演：把“现在看到的”还原成“过去发生的”
监测井中今天测到的浓度分布，是几十年释放与迁移叠加的结果。释放历史反演的目标是估计释放起始时间、释放持续时间与释放强度随时间的变化。常用手段有三类：一是特征污染物比值法，利用母体与降解产物的浓度比推断污染年龄；二是历史资料与生产记录对照，把产量、原料消耗与泄漏事件对应到时间轴；三是解析解或简化数值反演，用当前羽状分布反推源强。
一个可操作的简化做法是“一维对流弥散反演”。设地下水流速 v、纵向弥散度 αL，污染物前锋到达某监测井的距离为 L，则特征到达时间约为 t ≈ L/v（忽略弥散时）。若该井实测污染起始于 1998 年，而 L=320 m、v=8 m/a，则推算释放起点约在 1998−40=1958 年，与厂区建厂年代对照即可校验合理性。若考虑弥散，前锋到达会更早、峰值更晚，需用解析解修正。反演结果通常表达为“释放情景”：如“1970—1990 年高强度持续释放，1990 年后衰减为残余释放”。这一情景将作为数值模型的源项输入。
三、水文地质概念模型要素
水文地质概念模型要回答“污染物在什么样的介质中、以什么驱动力迁移”。核心要素包括：含水层结构与分层（潜水含水层、承压含水层、弱透水层的位置与厚度）；岩性组合与渗透性（砂、粉土、黏土的垂向组合决定污染物是垂向入渗还是侧向扩散）；地下水流场（流向、水力梯度、流速）；边界条件（河流、抽水井、隔水边界）；以及污染物的地球化学环境（氧化还原电位、pH、吸附与降解能力）。
以某冲积平原场地为例（示例假设）：0～6 m 为粉质黏土包气带，6～18 m 为细砂潜水含水层，18～22 m 为淤泥质弱透水层，22 m 以下为中砂承压含水层。地下水流向自西北向东南，水力梯度约 0.003，实测流速 6～10 m/a。这样的结构意味着：污染物在包气带中垂向迁移慢、吸附强；进入潜水含水层后以侧向扩散为主；弱透水层是重要的垂向阻隔，但一旦被击穿，承压含水层将形成新的迁移通道。概念模型必须把这些“结构控制点”明确标出。
四、污染羽形态的控制因素
污染羽的最终形态是源、介质、流场和污染物性质共同作用的结果。源强与释放时长决定羽的长度与质量；含水层渗透性的非均质性造成指状推进和绕流，使羽缘不规则；地下水流向决定羽的延伸方向，流速决定延伸速度；弥散度控制羽的宽度与浓度梯度；吸附与降解则使羽的前缘衰减、浓度峰值后移。对于密度大于水的 DNAPL，垂向迁移会形成“下沉型”羽；对于轻质 LNAPL，则易在潜水面附近形成“漂浮型”羽。
判断羽形态是否“合理”，可用几个简单指标交叉校验：羽的长宽比是否与流场方向一致；浓度等值线是否沿流向拉长；羽的前缘是否与某条地层界线或抽水井位置吻合。若模型给出的羽形态与这些地质、水动力证据矛盾，应优先怀疑源项或渗透性参数，而不是急于调弥散度。
把上述四部分整合成一份可交付的概念模型，通常以“源项情景表 + 水文地质剖面图 + 羽形态判读说明”三件套呈现。源项情景表列出释放起止年份、源强区间与特征污染物；剖面图标明分层、流场与监测井位置；判读说明则用长宽比、等值线走向和地层界线交叉验证羽形态的合理性。概念模型不是一次成型的，随着新监测数据加入，应定期复核源项与渗透性参数，使模型始终与证据保持一致。