4. 污染源与释放历史分析

污染羽的现状不是凭空出现的，它是源区在时间轴上持续或间歇释放的结果。要判断污染羽未来往哪走、会不会到达某个受体，第一步必须回答三个问题：源在哪里、释放有多强、释放历史与今天的羽状分布是否自洽。本节围绕这三问展开，重点讲清楚如何把历史生产信息转化为可用的源强参数，并用它解释污染羽的现状。

一、源区识别：从“浓度最高点”到“释放单元”

初学者常把监测井中浓度最高的位置直接当作源区，这在工程上往往不准确。源区应理解为一个释放单元：具有相对明确的边界、持续的污染物质量储备、以及与下游污染羽存在可追溯的水力与地球化学联系。识别源区通常需要三类证据交叉验证。

第一是空间浓度结构。在源区附近，浓度等值线往往呈现闭合高值核，且垂向上可能出现分层富集，例如DNAPL类污染物在含水层底部形成残留饱和带。第二是地球化学与同位素指纹。特定污染物或降解产物的比值、稳定同位素组成可以区分不同释放单元，判断污染羽是否来自单一源还是多源叠加。第三是历史土地利用与生产布局。原料储罐区、装卸栈台、排污沟、事故泄漏点、废渣堆场，往往对应不同的释放机制与释放强度。

实际操作中，可先圈定“疑似源区候选清单”，再逐项用证据打分。例如某化工场地历史平面图显示，1990年代存在一个溶剂储罐区和一个废液暂存池，两者相距约80米。监测数据显示储罐区下游浅层含水层中三氯乙烯浓度达数十毫克每升，而暂存池下游以苯系物为主。同位素与组分指纹表明两者属于不同释放单元，因此应作为两个独立源区分别估算源强，而不是合并为一个“场地源”。

二、释放速率估算：从历史产量反推质量通量

释放速率（源强）是污染羽模拟与受体暴露评估的核心输入。工程上常用三类方法，精度与数据需求各不相同。

方法一：历史物料平衡法。依据生产记录、原料采购量、产品产量与工艺收率，估算污染物的“输入量减去向量”。例如某车间1995—2005年累计使用三氯乙烯约120吨，工艺回收与产品带走合计约95吨，则潜在损失约25吨。再结合泄漏年限与泄漏面积，可粗算平均释放速率。该方法适合有较完整生产台账的场地，但需注意区分正常损耗、废气排放与进入土壤地下水的部分。

方法二：源区质量估算与释放时间反演。通过源区土壤和地下水浓度、含水层厚度、孔隙度与分配系数，估算源区当前残留质量，再结合释放持续时间反推历史平均通量。示例参数：源区面积约600平方米，污染层厚4米，孔隙度0.3，实测土壤中污染物平均浓度按合理假设取80毫克每千克，则可估算残留质量量级，再与历史损失量对照，判断释放是否已基本衰减。

方法三：污染羽通量法。在源区下游设置垂直于地下水流向的监测断面，利用浓度与达西通量估算穿过断面的质量通量。该方法反映的是“当前”释放强度，适合判断源区是否仍在活跃释放。若断面通量远低于历史平均释放速率，说明源区已进入衰减阶段，污染羽未来扩展速度将放缓。

三、历史生产信息与污染羽现状的关联

历史生产信息不是背景资料，而是解释污染羽形态的关键证据。生产起止时间决定释放持续时间，工艺变更决定污染物种类与浓度变化，泄漏事故决定脉冲式释放，厂区改造决定源区位置迁移。把这些时间节点与污染羽的空间分布叠加，可以判断污染羽处于扩展期、稳定期还是衰减期。

例如，某场地1998年投产、2008年停产，主要污染物为氯代烃。若下游污染羽前锋距离源区约300米，含水层地下水流速按合理假设为每年15米，则粗略推断污染羽前锋迁移时间约20年，与“停产至今约十余年、生产期约十年”的时间框架基本吻合，说明污染羽仍可能处于缓慢扩展或趋于稳定的过渡阶段。若前锋距离明显小于按流速推算的距离，则需考虑吸附、降解或源区实际释放时间较短等因素。

下图展示了源区识别、释放速率估算与历史信息关联三者的逻辑关系，以及它们如何共同支撑污染羽现状解释与受体暴露判断。

四、把分析落到可执行的判断上

完成上述分析后，应形成一份简明的源区与释放历史结论，至少包含：源区数量与位置、各源区污染物种类、释放起止时间、历史平均释放速率与当前释放趋势、以及与污染羽前锋位置的时间一致性判断。若历史物料平衡估算的累计释放量远大于源区残留量与已迁移通量之和，说明可能存在未识别的释放路径或源区，需要补充调查；若当前断面通量仍较高而污染羽前锋已接近受体，则应优先将该源区列为风险管控对象。

需要强调的是，历史生产信息往往存在缺失或口径不一，释放速率估算结果应给出合理区间而非单一数值。在后续污染羽迁移趋势预测中，建议采用保守情景与基准情景并行，观察受体到达时间对源强参数的敏感程度，从而判断哪些源区参数值得进一步投入调查精度。