基于多排污口核查与监测的中小流域污染源源解析方法与流程

本发明属环境监测技术领域，具体涉及一种基于多排污口核查与监测的中小流域污染源源解析方法。

背景技术：

目前水环境污染源解析主要依靠模型进行计算，模型主要有两种：

一种是苏丹等人于2009年在《水环境污染源解析研究进展》一文中所提出的模型，是以污染源为对象的正向扩散模型；另一种是以污染受体为对象的反向溯源模型，即受体模型。

扩散模型是一种预测式模型，它是通过输入各个污染源的排放数据和相关参数信息来预测污染物的时空变化情况。由于扩散模型参数复杂，解析繁琐，在水环境源解析中的应用还较少。

受体模型因不受污染源源强的限制，不依赖于距离、扩散系数等多种特性参数而得到广泛应用。

目前常见的水环境中污染物的受体模型源解析方法，主要包括因子分析法和主成分分析法；因子分析法和主成分分析法都被成功应用于水环境中污染物的源解析工作。这两种方法在应用中需要进行以下假设，如源成分谱在从源到受体这段距离没有显著变化、单个污染物通量的变化与浓度成比例、在给定时段污染物的总通量是所有已知源污染物通量的总和、源成分谱和贡献率都线性无关、所有采样点均主要受几个相同源的影响等。而实际水体往往要复杂的多，上述的两种方法往往在实际应用过程中并不能很好的适用于实地水体解析工作。

除此之外，这两种方法需要输入大量数据，对受体样品数量要求较高，一般超过50个，而且要想得到各类源的绝对贡献值，还需要再进行回归分析。并且，因子分析法不是对具体数值进行分析而是对偏差进行处理，如果某重要排放源比较恒定，而其他非重要源具有较大的排放强度变异，则可能会忽略排放强度较大的排放源。在实际使用中多元统计法一般鉴别出5～8个因子，如果重要排放源类型>10，那么这种方法不能提供较好的结果[gordonge.receptormodels[j].environmentalscienceandtechnology,1980,14(7):792-799]，即其结果具有不确定性。

因此，需要针对上述方法的缺陷进行改进，发明一种能适用于多排污口河流污染源源解析方法，以克服参数较多、结果不确定性高这些不足。

技术实现要素：

为了解决上述的技术问题，本发明建立了一种基于多排污口核查与监测的中小流域污染源源解析方法，该方法解决了现有方法所表现出的特定适用性、参数较多、结果不确定性较高的问题，提高了源解析结果客观性。

本发明主要是通过核查与监测的各入河排污口排放量，对流域污染源进行解析，得到每个排放源以及每类排放源的污染贡献；本发明尤其适用于多排污口、多支流的中小流域。

为了解决上述的技术问题，本发明所提供的基于多排污口的中小流域污染源源解析方法，具体如下：

(1)入河排污口核查

广泛搜集污染源及排污设施的相关资料，开展现场查勘，确定入河排污口分布状况和具体位置。采用gps法准确测量排污口地理位置，坐标系统为wgs-84坐标系统，高程系统为1956黄海高程系统。采用gps流动站的rtk模式准确测量出排污口的地理位置，绘制电子地图(jpg和shp数据格式)。

同时针对每个监测断面的具体情况和水流条件确定流量测验方法和测流方案。对所有连续或间歇向河流排污的入河排污口(不含入河的雨水排放口，农田沥水及涝水、退水排放口，未作为排污用的截洪沟和导洪渠入口)，包括通过涵洞、沟渠、管道进行排污的，要逐个核查。核查的内容包括入河排污口的位置、经纬度坐标、废污水性质、排放方式和水功能区等信息。

针对往年入河排污口监测成果中污染物含量过高的排污口，应根据其排污类型及特点着重调查；生活类排污口调查其污水来源，工业类排污口须了解企业污水处理能力，污水处理厂排污口须明确其处理级别以及中水用途等。原则上污染物含量过低的排污口不再监测。

根据核查结果，汇总核查表，包括入河排污口名称、入河排污口编码、详细地址、法人代表、单位性质、取用水量、服务面积和服务人口、污水性质、入河方式、排放方式、区县等；然后对核查成果进行编报。

排污口名称：按现有名称。

入河排污口编码：

干流：前面冠以gl，后面是编号，编号顺序从上游到下游，以两位阿拉伯数字进行记录，个位编号前面补充“0”；

支流：前面冠以zl，后面是编号，编号顺序从上游到下游，以两位阿拉伯数字进行记录，个位编号前面补充“0”；

排污口(沟)：前面冠以pw，后面是编号，编号顺序从上游到下游，编号方式：以两位阿拉伯数字进行记录，个位编号前面补充“0”；

溢流口：前面冠以yl，后面是编号，编号顺序从上游到下游，编号方式：以两位阿拉伯数字进行记录，个位编号前面补充“0”；

工业排污：前面冠以gp，后面是编号，编号顺序从上游到下游，编号方式：以两位阿拉伯数字进行记录，个位编号前面补充“0”；

污水处理厂：前面冠以szc，后面是编号，编号顺序从上游到下游，编号方式：以两位阿拉伯数字进行记录，个位编号前面补充“0”；

排污口位置电子地图的编制：

数据来源为gps(仪器型号：trimbler8，基站采用黄委会济南勘测局沿黄专用基站)。现场定位测量，坐标系统为：wgs-84坐标系统，高程系统为：1956黄海高程系统。

shp格式数据，导出后自动形成4个文件，分别是：.dbf、.prj、.shx和.shp，要查看该文件时，需要事先下载bigemap地图下载器方能打开但不能编辑，如要编辑需要预备加密程序。

编辑和查阅入河排污口的电子地图时：

1)标明支流、溢流口、直排企业等在河流上的位置；

2)标明流量测定采样分析的点位，就是现场核实确定的点位。一般支流分析点位大都往支流上游移动了些许。排入支流的企业标明其厂内排污口位置。

(2)入河口水量水质同步监测

在核查成果的基础上，制定水位、水量、水质监测方案，对入河排污口进行取样监测，开展水量水质同步监测。监测分不同水期开展，包括枯水期、丰水期、平水期等。

枯水期对每个排污口连续监测2天，每天间隔采样3次，取样时间分别为8：00、17：00和24：00点，为确保不受其他因素的影响，每个排污口的第一次取样时间务必为17：00，六次取样样品分别检测。监测采样时间应选择在前三天无降水日，避免受降水和河水顶托影响。入河排污口监测常规项目为：流量、流速、水位、化学需氧量(codcr)、氨氮(nh3-n)、总磷、总氮和浊度共8项。按照《污水综合排放标准》(gb8978-1996)、《地表水环境质量标准》(gb3838－2002)中水质分析方法进行，各项目分析方法按照各实验室资质认定许可的国家或行业执行。

丰水期对降雨后情况进行监测，在普遍降雨3天之后采样，时间分别为：8：00、17：00和24：00，对每个排污口连续监测2天。监测项目和监测方法同枯水期。

平水期对每个排污口连续监测2天，间隔采样3次，取样时间分别为8：00、17：00和24：00点。监测采样时间应选择在前三天无降水日，避免受降水和河水顶托影响。监测项目和监测方法同枯水期。

监测过程中严格各个采样和分析的质控工作，包括空白、平行样、质控样、加标回收率等。如有必要还需对重点排污单位进行金属和有机物特征污染物监测。最后编制监测报告。

对于重要的入河支流和重要的排污口监测断面要观测水位，按照《水位观测标准》(gb/t50138-2010)，首先在适合设立水尺的监测断面，设立临时水尺观测水位，水尺零点高程采用水准仪用不低于四等水准的要求，从附近国家三等以上水准点上引测其高程，然后人工读取水尺读数加上水尺零点高程得到相应监测断面的水位。

在不适合设立水尺的监测断面，可以采取如下方法：岸壁上沿、桥梁边沿等处设置固定点水准点，采用gps的rtk模式测定固定点的坐标和高程，用钢尺或测深杆量取国定水准点到水面的垂直距离，计算出支流或排污口监测断面的水位。

水量按照水文测流要求，可采用流速仪法、浮标法、溢流堰法、容量法等。

污水处理厂采用其在线监测设备的流量成果。根据监测断面水流条件和具体情况，首先采用流速仪--面积法进行流量测量：用流速仪测量每条垂线的流速，用测深仪或测深杆测量水深，用钢尺或全站仪测量出断面的起点距，然后按照流速-面积法计算出断面流量。

在不适合转子流速仪法的监测断面，采用超声波表面流速仪，测出每条垂线的表面流速，然后进行河道断面测量，用流速--面积法计算出表面虚流量，然后乘以系数(系数在0.8～0.9之间，一般取0.85)得实际流量。

不适合以上流量测验法的监测断面，可以采用双浮标法测流，水面浮标采用草把、木棒制作，下面坠以石块，上面插上易于观测的小旗。在夜间可以用电池、电珠或棉球火炬安装在浮标上作成夜明浮标。浮标下端坠以重物，上端加标志物，入水深度占全水深的9/10及以上，可以近似测得垂线平均流速，然后用监测断面的实测过水面积乘以浮标流速，得断面流量。

对于有溢流坝、涵闸结构的监测断面，采用测量闸、上游水位或上游行进河槽的水位推算出过水断面的水头，利用堰闸、溢流堰公式推算流量。

对于排污口用排污管道的排水泄流的监测断面，如果水流不大，同时容易被容器接纳的监测断面，可以采用有标准容积的容器，接收水流，从开始接纳计时到容器接满为止并计时，用容器的容积除以接纳时间得出流量。

对于闸控或用泵站抽排的入河排污口，必须通过调查，掌握其抽排规律，给出在无地表径流条件下的平均入河流量。

有地表径流影响的入河排污口或者是污染严重的沟渠，要设置并注明本底断面，计算排污量要扣除上游来量和本底污染物量。

所监测的各入河排污口排放量之和应占本河段入河排污总量的90％以上。

最后编制监测报告。

(3)污染源源解析

统计各水期各入河排污口水量水质情况，计算各水期水量和水质平均值，解析各排污口对流域的污染贡献。若对各类排污口进行分类汇总，可解析不同水期各类排污口对流域的污染贡献，从而提出针对该流域水质达标的控源截污方案。也可根据不同目的，进行细致源的源解析。

与现有源解析方法相比，本发明的有益效果在于：

(1)本发明所提出的一种基于多排污口核查与监测的中小流域污染源源解析方法，通过对排污口进行核查与监测，避免了源解析模型具有一定应用条件、参数较多、解析结果不确定性较高的缺点，将90％以上负荷的污染源均纳入流域污染源源解析，可得出每类排放源乃至每个排放源的污染贡献；

(2)本发明解析结果客观，接近实际排放情况，从而为精准控源提供了靶向，提高了流域污染治理的效率。

(3)与现有的因子分析法和主成分分析法相比，本发明的方法具有广泛的适用性，适用于所有中小流域，具有广泛的推广和应用价值。

附图说明

图1小清河入河排污口分布图(部分河段)；

图2不同时期各入河排污口水量占比(％)；

图3不同时期各类排污口污染物排污占比情况(％)；

图4溢流口封堵后小清河济南段出境断面氨氮浓度变化。

具体实施方式

下面结合具体实施例对本发明作更进一步的说明，以便本领域的技术人员更了解本发明，但并不因此限制本发明。

以济南的小清河为例，实施例1具体提供了“小清河济南段污染源源解析方法”，解析时间为2018年。

基于多排污口核查与监测的中小流域污染源源解析方法，包括以下的步骤：

(1)入河排污口核查

广泛搜集污染源及排污设施的相关资料，开展现场查勘，确定入河排污口分布状况和具体位置。采用gps法准确测量排污口地理位置，坐标系统为wgs-84坐标系统，高程系统为1956黄海高程系统。用gps流动站的rtk模式，准确测量出排污口的地理位置和并绘制了电子地图(jpg和shp数据格式)。

同时针对每个监测断面的具体情况和水流条件确定流量测验方法和测流方案。对所有连续或间歇向河流排污的入河排污口(不含入河的雨水排放口，农田沥水及涝水、退水排放口，未作为排污用的截洪沟和导洪渠入口)，包括通过涵洞、沟渠、管道进行排污的，要逐个核查。核查的内容包括入河排污口的位置、经纬度坐标、废污水性质、排放方式和水功能区等信息。

针对往年入河排污口监测成果中污染物含量过高的排污口，应根据其排污类型及特点着重调查；生活类排污口调查其污水来源，工业类排污口须了解企业污水处理能力，污水处理厂排污口须明确其处理级别以及中水用途等。原则上污染物含量过低的排污口不再监测。

根据核查结果，汇总核查表，包括入河排污口名称、入河排污口编码、详细地址、法人代表、单位性质、取用水量、服务面积和服务人口、污水性质、入河方式、排放方式、区县等；然后对核查成果进行编报。

排污口名称：按现有名称。

入河排污口编码：

干流：前面冠以gl，后面是编号，编号顺序从上游到下游，以两位阿拉伯数字进行记录，个位编号前面补充“0”；

支流：前面冠以zl，后面是编号，编号顺序从上游到下游，以两位阿拉伯数字进行记录，个位编号前面补充“0”；

排污口(沟)：前面冠以pw，后面是编号，编号顺序从上游到下游，编号方式：以两位阿拉伯数字进行记录，个位编号前面补充“0”；

溢流口：前面冠以yl，后面是编号，编号顺序从上游到下游，编号方式：以两位阿拉伯数字进行记录，个位编号前面补充“0”；

工业排污：前面冠以gp，后面是编号，编号顺序从上游到下游，编号方式：以两位阿拉伯数字进行记录，个位编号前面补充“0”；

污水处理厂：前面冠以szc，后面是编号，编号顺序从上游到下游，编号方式：以两位阿拉伯数字进行记录，个位编号前面补充“0”；

详细地址：根据现场调查并结合地图进行填写。

法人代表：只填写污水处理厂和企业的。

单位性质：按实际情况填写。

入河方式：直接排放、间接排放。

污水性质：分支流、生活污水、工业污水、污水厂排水等。

区县：按实际所在区县填写。

地理坐标：摘录现场gps测量数据。

以上内容，以电子表格形式记载录入，并形成电子文档。

排污口位置电子地图的编制：

数据来源为gps(仪器型号：trimbler8，基站采用黄委会济南勘测局沿黄专用基站)。现场定位测量，坐标系统为：wgs-84坐标系统，高程系统为：1956黄海高程系统。

shp格式数据，导出后自动形成4个文件，分别是：.dbf、.prj、.shx和.shp，要查看该文件时，需要事先下载bigemap地图下载器方能打开但不能编辑，如要编辑需要预备加密程序。

编辑和查阅入河排污口的电子地图时：

1)标明支流、溢流口、直排企业等在小清河干流上的位置；

2)标明流量测定采样分析的点位，就是现场核实确定的点位。一般支流分析点位大都往支流上游移动一定距离，避免回灌影响。排入支流的企业标明其厂内排污口位置。

对小清河济南段入河排污口进行现场核查，共核查排污口50个，其中支流25条，污水处理厂6个，排污口6个，工业排污口4个，溢流口5个，干流断面2个。

25条支流为：腊山河、南太平河、北太平河、兴济河、曹家圈虹吸干河、西工商河、前引河、后引河、西泺河、东泺河、柳行头河、全福河、华山大沟、大辛河、赵王河、小汉峪沟、冷水沟、龙脊河、石河、刘公河、杨家石河、巨野河、绣江河、章齐沟和漯河。其中4月份截流截污支流有9个：兴济河、北太平河、前引河、后引河、西泺河、全福河、华山大沟、赵王河和刘公河，河干支流有2个：绣江河和漯河(夏侯桥)；8月份截流截污河流为前引河；10月份截污截流河流有4个：前引河、华山大沟、赵王河、刘公河。

6个污水处理厂为：大金污水厂、光大四厂、光大二厂、光大一厂、光大三厂和临港污水厂。

6个排污口为：南太平河入小清河口东10米小清河北岸方形涵洞、王舍人北滩头桥东生活入河排污口、遥墙排污沟(大码头村)、机场路桥下西北排污口、董家镇排污沟和遥墙老工业园东排污口。其中南太平河入小清河口东10米小清河北岸方形涵洞排污，但是一直处于淹没状态，但是在排污。遥墙老工业园东排污口和董家排污沟，污水已被截污进入临港污水厂处理后排放。其余排污口除了汛期基本为河干或截污状态。

4个工业排污口为：济南炼油厂、山东绿霸(厂外、厂外)、蓝星石油(厂内、厂外)、齐鲁制药华山还乡店桥东侧入河排污口。其中济炼和绿霸正常排放；蓝星石油已经停产、且不再复工；齐鲁制药厂华山还乡店桥东侧入河排污口，不再排污，污水已被截污进入光大三厂处理。

5个溢流口为：黄岗桥溢流口、济洛路桥西溢流口、济洛路桥东溢流口、生产路桥溢流口和板桥溢流口。主要溢流口为黄岗桥溢流口、生产路桥溢流口和板桥溢流口。溢流口大部分处于淹没或半淹没状态。

2个干流断面为：源头断面睦里庄和出境断面辛丰庄。

已纳入污水厂不向小清河排放废水的排污口，此次不再监测，包括齐鲁制药厂华山还乡店桥东侧入河排污口、蓝星石油、遥墙老工业园东排污口、董家排污沟和完全截污时的刘公河。

对于河干、截流、截污、断流(统称为死水，即流量为零的监测断面)只拍照，不取样，并要在水样递送单上注明。位于水质厂内或企业内，只取样监测，不再观测水位和流速数据，流量摘抄厂内的在线自动监测数据。其他监测断面需要进行水位观测、流量测验和水质取样化验。

接纳污水厂排水的支流监测总流量和水质情况。

小清河入河排污口核查汇总见表1。小清河入河排污口分布图(部分河段)见图1。

(2)入河排污口水量水质同步监测

通过现场核查情况，最后确定排污口监测断面中，进行流量测验、水位观测和水质取样的站点45个。对这些入河排污口进行取样监测，开展水量水质同步监测。监测分不同水期开展，包括枯水期、丰水期、平水期等。

枯水期对每个排污口连续监测2天，每天间隔采样3次，取样时间分别为8：00、17：00和24：00点，为确保不受其他因素的影响，每个排污口的第一次取样时间务必为17：00，六次取样样品分别检测。监测采样时间应选择在前三天无降水日，避免受降水和河水顶托影响。入河排污口监测常规项目为：流量、流速、水位、化学需氧量(cod)、氨氮(nh3-n)、总磷、总氮和浊度共8项。按照《污水综合排放标准》(gb8978-1996)、《地表水环境质量标准》(gb3838－2002)中水质分析方法进行，各项目分析方法按照各实验室资质认定许可的国家或行业执行。

丰水期对降雨后情况进行监测，在普遍降雨3天之后采样，时间分别为：8：00、17：00和24：00，对每个排污口连续监测2天。监测项目和监测方法同枯水期。

平水期对每个排污口连续监测2天，间隔采样3次，取样时间分别为8：00、17：00和24：00点。监测采样时间应选择在前三天无降水日，避免受降水和河水顶托影响。监测项目和监测方法同枯水期。

监测过程中严格各个采样和分析的质控工作，包括空白、平行样、质控样、加标回收率等。如有必要还需对重点排污单位进行金属和有机物特征污染物监测。最后编制监测报告。

对于重要的入河支流和重要的排污口监测断面要观测水位，按照《水位观测标准》(gb/t50138-2010)，首先在适合设立水尺的监测断面，设立临时水尺观测水位，水尺零点高程采用水准仪用不低于四等水准的要求，从附近国家三等以上水准点上引测其高程，然后人工读取水尺读数加上水尺零点高程得到相应监测断面的水位。

在不适合设立水尺的监测断面，可以采取如下方法：岸壁上沿、桥梁边沿等处设置固定点水准点，采用gps的rtk模式测定固定点的坐标和高程，用钢尺或测深杆量取国定水准点到水面的垂直距离，计算出支流或排污口监测断面的水位。

水量按照水文测流要求，可采用流速仪法、浮标法、溢流堰法、容量法等。污水处理厂采用其在线监测设备的流量成果。

所监测的各入河排污口排放量之和应占本河段入河排污总量的90％以上。

1)各类入河排污口流量监测结果

对各入河排污口进行分类、统计，见表2。由表可知，小清河的主要来源是污水厂排水，其次是支流、溢流口和源头。污水厂排水为97.47～105.64万m³/d，支流流量为18.70～121.20万m³/d，溢流口流量为6.49～15.55万m³/d，睦里庄水量为0～28.08万m³/d，辛丰庄出境断面流量为127.01～271.30万m³/d。不同时期内，除了污水厂水量和企业水量变化较小，其余类别如源头、支流和溢流口，水量变化均较大。支流总量变化最大，其次是源头睦里庄。支流流量变化较大的是腊山河、南太平河、北太平河、兴济河、后引河、西泺河、东泺河、龙脊河、石河和刘公河。溢流口变化较大的是黄岗桥溢流口和生产路桥溢流口。

表2小清河各入河排污口流量监测结果

2)各类入河排污口水质监测结果

2018年小清河各入河排污口水质汇总见表3。

表3小清河入河排污口水质监结果

“/”表明：该排污口没水，所以没有相应的数据。

(3)小清河污染源源解析

统计各水期各入河排污口水量水质情况，计算各水期水量和水质平均值，解析各排污口对流域的污染贡献。若对各类排污口进行分类汇总，可解析不同水期各类排污口对流域的污染贡献，从而提出针对该流域水质达标的控源截污方案。

1)水量解析

不同时期各类别排污口流量占比见图2。该统计为不考虑水量损失的统计结果。由图可知，4月份水量占前三位的排污口排序为：污水厂(77.44％)>支流(14.86％)>溢流口(7.63％)；8月份水量占前四位的排序为：支流(44.73％)>污水厂(38.99％)>源头(10.36％)>溢流口(5.74％)；10月份水量占前四位的排序为：污水厂(55.75％)>支流(35.92％)>源头(4.55％)>溢流口(3.57％)。在枯水期，污水处理厂水量占比最高为77.4％，其次是支流为14.86％，然后是溢流口为7.63％。在汛期，支流水量占比最高，为44.73％，然后是污水厂为38.99％，其次是源头10.36％；在平水期，污水厂排水比例最高，为55.75％，其次是支流为35.92％，然后是源头为4.55％。

在各个时期，污水处理厂和支流的总和对小清河水量贡献最大，其中枯水期占总排水量的92.30％，汛期占83.73％，平水期占91.68％。随着污水处理厂的扩建，溢流口排污将得到解决，而支流流量将不容小觑。支流代表地表径流产生的排污，将是小清河水质改善下一步治理的重点。

2)流域污染源贡献解析

此实例中对各类排污口进行分类汇总，分为污水厂、溢流口、支流以及其他四类。结合流量统计出污染物排放量占比情况，重点解析了小清河codcr、氨氮、总氮、总磷的污染源排放量贡献率，见附图3。

codcr污染解析：4月份cod的污染贡献率为：污水厂(54.67％)>溢流口(31.44％)>支流(13.79％)，8月份cod的污染贡献率为：支流(40.12％)>溢流口(31.87％)>污水厂(16.31％)，10月份cod的污染贡献率为：支流(50.87％)>溢流口(23.27％)>污水厂(21.66％)。可见，污水厂对cod的贡献只有在枯水期占比最大，在汛期和平水期占比均下降。汛期和平水期时支流对cod贡献最大，其次是溢流口。支流的贡献主要体现在汛期和平水期，说明在点源得到有效控制之后，地表径流污染的控制是下一步的重点。

氨氮污染解析：4月份氨氮的污染贡献率为：溢流口(51.39％)>支流(26.92％)>污水厂(21.67％)，8月份氨氮的污染贡献率为：支流(49.88％)>溢流口(31.64％)>污水厂(8.24％)，10月份氨氮的污染贡献率为：溢流口(50.25％)>支流(31.94％)>污水厂(16.84％)。可见，枯水期和平水期时，溢流口对氨氮的贡献最大，其次是支流，然后是污水处理厂。在汛期，支流对氨氮的贡献最大，其次是溢流口，然后是污水处理厂。因此，为使小清河氨氮达标最迫切的是需要将溢流口污水进行处理后排放，然后是开展支流治理以及地表径流污染控制。

总氮：4月份总氮的污染贡献率为：污水厂(68.16％)>支流(28.49％)>溢流口(3.32％)，8月份总氮的污染贡献率为：支流(44.23％)>污水厂(37.98％)>源头(10.79％)，10月份总氮的污染贡献率为：支流(55.11％)>污水厂(34.25％)>溢流口(8.61％)。可见，枯水期时污水厂对总氮的贡献最大，其次是支流，然后是溢流口；在汛期，支流对总氮的贡献最大，其次是污水处理厂，然后是源头。在平水期，支流对总氮的贡献最大，其次是污水厂，然后是溢流口。因此，为使小清河总氮达标，首先需要对污水处理厂出水进行严格控制，然后是开展支流以及地表径流污染控制。整体来看污水厂和支流对小清河总氮贡献最大，两者几乎相当，其次是溢流口，雨季时源头的贡献率上升。

总磷：4月份总磷的污染贡献率为：污水厂(34.92％)>溢流口(34.63％)>支流(30.43％)，8月份总磷的污染贡献率为：支流(43.23％)>溢流口(38.25％)>污水厂(12.43％)，10月份总磷的污染贡献率为：支流(68.84％)>溢流口(16.96％)>污水厂(13.16％)。可见，枯水期时，污水厂对总磷的贡献最大，其次是溢流口，然后是支流。在汛期和平水期，支流对总磷的贡献最大，其次是溢流口，然后是污水处理厂。因此，为使小清河总磷达标，需要开展支流以及地表径流污染控制，然后将溢流口污水纳入污水厂处理。

目前，小清河济南段主要超标因子为氨氮，解析结果表明在无降雨时溢流口的贡献最大，超过50％。溢流口污水是经市政管网收集的生活污水，由于污水厂处理能力不足这部分污水无法处理，从而直接排入小清河。因此，为使小清河氨氮达标最迫切的是需要将溢流口污水处理后排放。2019年1月济南市光大水务一厂和二厂各扩建10万吨/天产能，将溢流口污水逐渐引入污水厂进行处理，然后再排入小清河，小清河水质出现明显好转，见图4。图4为溢流口封堵后小清河济南段出境断面氨氮浓度变化图。图4表明2019年济南市小清河出境断面氨氮监测数据明显优于2018年。自2019年2月至今，氨氮指标均达到《地表水环境质量标准》v类标准。说明解析结果准确、可靠，为解决流域问题提出了精准解决方案。