CN111598754B - 基于水质模型的河网区水质断面污染贡献率核算方法 - Google Patents

Abstract

本发明提供一种基于水质模型的河网区水质断面污染贡献率核算方法，包括：S1、根据源质空间关系构建影响各水质断面的备选污染源清单；S2、根据河网区的特征建立河网水动力水质模型，并对所述河网水动力水质模型进行验证率定；S3、根据备选污染源清单，将污染源相关数据作为输入条件输入至河网水动力水质模型；S4、依次对目标水质断面的每个备选污染源按照排除法进行情景模拟，计算每个备选污染源对目标水质断面的贡献率并排序。本发明的有益效果：将现有成熟的水质模型与大量实测统计数据，通过数值模拟得到各污染源对水质断面的影响贡献率大小，能为水环境管理领域对污染源调控方案作出提前预测，提供理论支持。

Description

基于水质模型的河网区水质断面污染贡献率核算方法

技术领域

本发明涉及水环境管理领域，尤其涉及一种基于水质模型的河网区水质断面污染贡献率核算方法。

背景技术

在水环境模型领域，河流水质模型从20世纪30年代开始起步，截至目前，发展成熟的模型已有QUAL系列、WASP5、MIKE11、EFDC等几十个，逐步从简单单一化发展成能从不同维度解析、分析多种变量，可结合多种技术方法进行综合运用的工具，河流水质模型强大的功能能够针对不同的问题进行运用分析。但目前水质模型的分析主要结合ANN、灰色理论等方法运用在分析预测、评估、治理等方面，解决水功能区水质解析、水质动态模拟等问题，极少结合排除法及贡献率计算方法解析污染源对水质断面的影响问题。

目前针对污染源对水质断面的污染贡献影响研究较少，水质断面贡献率的研究还没有相关的方法和设备，大部分关于贡献率的研究都与污染源排放贡献率的计算方法有关，比如输出系数法、排污系数法、监测法等，很少有方法针对污染源对水质断面的污染贡献进行研究。且排除法一般应用于理论知识学习、计算机科学领域，极少涉及水环境精细化管理领域。

相关科研结果表明，通过大量有效的实测统计数据对水质模型进行验证和率定后，水质模型的预测能力可达到尽可能还原实际水质状况的要求，从而使得水环境管理领域不仅仅停留在假象阶段，更具实际借鉴意义。

发明内容

有鉴于此，本发明提供了一种基于水质模型的河网区水质断面污染贡献率核算方法，利用实测统计的相关数据，通过构建适用于相关河网区的水质模型，结合污染源全排和依次停排的排除法，获得污染源在不同排放情况下相应水质断面的主要控制污染物的数值模拟值，以主要控制污染物浓度之和代表污染源，再利用水质断面贡献率计算方法得到不同排放条件下各污染源对水质断面贡献率，并根据贡献率大小获得排名，能为水环境管理领域对污染源调控方案作出提前预测，提供理论支持。

本发明提供一种基于水质模型的河网区水质断面污染贡献率核算方法，包括以下步骤：

S1、根据河网区的水流方向以及水质断面和污染源之间的上下游关系建立水质断面与污染源清单之间的空间关系，确定能够影响到水质断面的备选污染源，从而构建影响各水质断面的备选污染源清单；

S2、根据河网区的特征建立河网水动力水质模型，并对所述河网水动力水质模型进行验证率定；

S3、根据所述备选污染源清单，将污染源相关数据作为河网水动力水质模型的输入条件，模拟各控制污染物沿河段距离动态变化的浓度值；

S4、根据备选污染源清单，明确影响目标水质断面的污染源集合，基于排除法，在河网水动力水质模型中设定污染源全排和依次停排的情况，模拟得到不同情况下控制污染物的浓度及目标水质断面的水流量，进一步计算每个备选污染源对目标水质断面的贡献率，并根据贡献率大小进行排序，从而确定水质断面的主要污染源。

进一步地，所述步骤S1的具体地过程为：

S11、根据河流方向，将河流重要干流支流进行概化，确定干流、一级支流和二级支流，并以水力学特征为基础，在污染源排放口下游污水与河水基本混匀处以及特殊要求地区的河段上设置水质断面，以水质断面、支流汇入干流口和特殊要求地区为节点，对河段进行划分并排序；

S12、确定河网区所属的水质控制功能区，明确各水质断面的水质控制目标，并将所述水质控制目标作为水质评价指标，确定各水质断面的控制污染物，构建包含水质断面所处河段序号、距离河段起始点的距离、控制污染物类型的水质断面清单；

S13、基于实测统计数据，考虑点源和非点源污染源，根据污染源对水质断面的影响大小，筛除部分污染源，并将相邻密集的点源污染源概化为集中的污染源，非点源污染源则记录下所影响河段的起始点与终点距离河段上边界断面的位置，参照由水质控制目标所确定的各水质断面的控制污染物，构建包含污染源所处河段序号、污染源类型、控制污染物的排放浓度、排放流量、排放口距离河段上边界的距离、属于出流还是入流污染源的污染源排放清单；

S14、比对所述水质断面清单与所述污染源清单中所记录的水质断面和污染源所处河段及距离，确定污染源排污口和水质断面之间的上下游关系，确定每个水质断面的备选污染源，构建水质断面的备选污染源清单。

进一步地，所述步骤S2的具体过程为：

根据所要研究的河网区的特征，选用合适的水质模型，将水质断面概化成底坡水平、断面矩形或梯形的简单几何形状，通过对实测水质断面的水深与过水面积关系进行回归分析，得到概化断面的下底宽和边坡系数，使得概化后的水质断面与实测水质断面保持一致的水位关系；

确定河段的河道坡度、河段水力参数、河段水质参数、流域水文资料、河道上游、下游边界条件和初始条件、断面节点的情况、流域气象数据以及水质监测站点所测污染物数据，按照所选用的水质模型的操作规则，将数据输入模型中，选择合适的计算步长进行计算，将所模拟的结果数据与水质监测站点的数据曲线进行比较，分析相对误差，若所述相对误差在允许范围之内，则说明建模成功，否则，根据实测数据修改模型参数，重新验证，直到构建得到适用于所要研究的河网区的水质模型。

进一步地，所述河段水力参数根据河段河道实际情况，选择性使用水力特征值的计算方法；所述河段水质参数根据河网区所处环境、污染物及污染物变化程度，通过实测统计数据或经验公式法获得河段的沉降、水解速率。

进一步地，所述步骤S3中，所述污染源相关数据包括每个污染源的流量、主要控制污染物的浓度、所处河段序号以及距离河段上边界距离。

进一步地，所述步骤S4中计算每个备选污染源对目标水质断面的贡献率的具体过程为：

模拟污染源全部排放时，向河网水动力水质模型中输入所有污染源的数据，由此得到所有污染源的各主要控制污染物在目标水质断面处的浓度模拟值，以及目标水质断面处的水流量模拟值；模拟污染源依次停排时，不向河网水动力水质模型中输入被停排的污染源的相关数据，由此得到未被停排的污染源的各主要控制污染物在目标水质断面处的浓度模拟值，以及此时目标水质断面处的水流量模拟值；多次进行上述模拟试验，最后计算污染源对于水质断面的贡献率：

式中，p_i表示被停排的污染源i对目标水质断面的贡献率，i表示被停排的污染源的编号，t表示按时间进行的试验次数，C_0t表示第t次试验时污染源全排的情况下所有污染源的主要控制污染物在目标水质断面处的浓度之和，Q_0t表示第t次试验时污染源全排的情况下目标水质断面处的水流量，C_st表示第t次试验时关停污染源i的情况下未被停排的污染源的主要控制污染物在目标水质断面处的浓度之和；Q_st表示第t次试验时关停污染源i的情况下目标水质断面处的水流量。

本发明提供的技术方案带来的有益效果是：将现有成熟的水质模型与排除法结合，设定污染源全排和依次进行停排的情况，根据现有的大量实测统计数据，通过水质模型的数值模拟，结合水质断面贡献率计算方法，得到各污染源对水质断面的影响贡献率大小，能为水环境管理领域对污染源调控方案作出提前预测，提供理论支持。

附图说明

图1是本发明实施例提供的基于水质模型的河网区水质断面污染贡献率核算方法的流程图；

图2是本发明实施例提供的对河网区河段进行划分并排序的结果示意图；

图3是本发明实施例提供的水质断面清单的示意图；

图4是本发明实施例提供的污染源排放清单的示意图；

图5是本发明实施例提供的备选污染源清单的示意图。

具体实施方式

为使本发明的目的、技术方案和优点更加清楚，下面将结合附图对本发明实施方式作进一步地描述。

请参考图1，本发明的实施例提供了一种基于水质模型的河网区水质断面污染贡献率核算方法，包括以下步骤：

S1、根据河网区的水流方向以及水质断面和污染源之间的上下游关系建立水质断面与污染源清单之间的空间关系，确定能够影响到水质断面的备选污染源，从而构建影响各水质断面的备选污染源清单。

步骤S1的具体地过程为：

S11、根据河流方向，将河流重要干流支流进行概化，确定干流、一级和二级支流，并以水力学特征为基础，在污染源排放口下游污水与河水基本混匀处以及特殊要求地区(饮用水源地、风景游览区等)的河段上设置水质断面，以水质断面、支流汇入干流口和特殊要求地区(行政边界、水库、闸泵站等)为节点，对河段进行划分，以保证同一河段具有大致相同的水力特性和参数；请参考图2，对河段进行排序，从主干河流源头至下游进行升序编号，当到达支流时，从支流源头开始编号；

S12、根据《地表示环境质量标准》确定河网区的用水主要属于五大地表水中的哪一类，获取其所属的水质控制功能区，从而明确各水质断面的水质控制目标，将其作为水质评价指标，确认各水质断面的主要控制污染物，构建包含水质断面所处河段序号、距离河段起始点的距离、水质控制目标、主要控制污染物类型的水质断面清单，如图3所示，其中所属河段的水质断面是指河段中出流一侧的水质断面；

S13、基于实测统计数据，考虑入河的多种污染源，包括点源和非点源(比如城市生活废水、工业废水、畜禽养殖、农田中的农药化肥等)污染源，根据污染源对水质断面的影响大小，筛除对河网区影响效果可以忽略的污染源，并将相邻密集的点源污染源概化为集中的污染源，非点源污染源则记录下其所影响河段的起始点与终点距离河段上边界断面的位置；参照由水质控制目标所确定的各水质断面的主要控制污染物，构建包含所处河段序号、污染源类型、主要控制污染物的排放浓度、排放流量、排放口距离河段上边界的距离、属于出流还是入流污染源的污染源排放清单，如图4所示；

S14、比对构建的水质断面清单与污染源清单中所记录的水质断面和污染源所处河段以及距离，确定污染源排污口与水质断面之间的上下游关系，确定每个水质断面的备选污染源，构建水质断面的备选污染源清单，如图5所示。

S2、根据河网区的特征建立河网水动力水质模型，在河网区中选择水质监测站点分布合理、实测数据较为完备的流域进行参数率定。

具体地，根据所要研究的河网区的特征，选用合适的水质模型，将水质断面概化成底坡水平、断面矩形或梯形的简单几何形状，通过对实测断面的水深与过水面积关系进行回归分析，得到概化断面的下底宽和边坡系数两个重要参数，使得概化后的水质断面与实测断面保持一致的水位关系；

同时确定河段的河道坡度、河段水质水力参数、流域水文资料、河道上游和下游边界条件及初始条件、断面节点的情况、流域气象数据以及水质监测站点所测污染物数据，按照所选用的水质模型的操作规则，将数据输入模型中，选择合适的计算步长进行计算，将所模拟的结果数据与水质监测站点的数据曲线进行比较，分析相对误差，以此开展水质模型的率定与验证，构建能够适用于此河网区水文水质模拟的河网水动力水质模型，若所述相对误差在允许范围之内，则说明建模成功，否则，应根据实测数据修改参数，重新验证。其中，河段水力参数根据河段河道实际情况，选择性使用水力特征值(流量、流速、水深)的计算方法，例如，若选用水位-流量关系曲线，可运用最小二乘法确定与流量流速有关的经验系数，若自行输入纵向弥散系数，则根据河底物质组成确定曼宁粗糙系数，再利用平均水深、平均流速和曼宁粗糙系数通过经验公式确定纵向弥散系数；河段水质参数根据河网区所处环境、污染物及污染物变化程度，通过实测统计数据或经验公式法获得沉降、水解速率等。

S3、根据备选污染源清单明确每个污染源的流量、主要控制污染物的浓度、所处河段序号以及距离河段上边界距离等数据，按照河流分段顺序依次将数据填入河网水动力水质模型的相应板块中进行数值计算，从而模拟得到各主要控制污染物沿河道距离动态变化的浓度值。

S4、根据备选污染源清单明确影响目标水质断面的污染源集合，基于排除法，在河网水动力水质模型中设定污染源全排和依次停排的情况，模拟得到不同情况下主要控制污染物的浓度以及目标水质断面的水流量；根据情景模拟结果和贡献率计算公式，计算每个备选污染源对目标水质断面的贡献率，并根据贡献率大小进行排序，从而确定水质断面的主要污染源，达到精确治理的目标要求。

其中，计算每个备选污染源对目标水质断面的贡献率的具体过程为：模拟污染源全部排放时，向河网水动力水质模型中输入所有污染源的数据，由此得到所有污染源的各主要控制污染物在目标水质断面处的浓度模拟值，以及目标水质断面处的水流量模拟值；模拟污染源依次停排时，不向河网水动力水质模型中输入被停排的污染源的相关数据，由此得到未被停排的污染源的各主要控制污染物在目标水质断面处的浓度模拟值，以及此时目标水质断面处的水流量模拟值；多次进行上述模拟试验，最后计算污染源对于水质断面的贡献率：

式中，p_i表示被停排的污染源i对目标水质断面的贡献率，i表示被停排的污染源的编号，t表示按时间进行的试验次数，C_0t表示第t次试验时污染源全排的情况下所有污染源的主要控制污染物在目标水质断面处的浓度模拟值之和(mg/L)，Q_0t表示第t次试验时污染源全排的情况下目标水质断面处的水流量模拟值(m³/s)，C_st表示第t次试验时关停污染源i的情况下未被停排的污染源的主要控制污染物在目标水质断面处的浓度模拟值之和(mg/L)；Q_st表示第t次试验时关停污染源i的情况下目标水质断面处的水流量模拟值(m³/s)。

在本文中，所涉及的前、后、上、下等方位词是以附图中零部件位于图中以及零部件相互之间的位置来定义的，只是为了表达技术方案的清楚及方便。应当理解，所述方位词的使用不应限制本申请请求保护的范围。

在不冲突的情况下，本文中上述实施例及实施例中的特征可以相互结合。

以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (1)

1.基于水质模型的河网区水质断面污染贡献率核算方法，其特征在于，包括以下步骤：

S1、根据河网区的水流方向以及水质断面和污染源之间的上下游关系建立水质断面与污染源清单之间的空间关系，确定能够影响到水质断面的备选污染源，从而构建影响各水质断面的备选污染源清单；

S2、根据河网区的特征建立河网水动力水质模型，并对所述河网水动力水质模型进行验证率定；

S3、根据所述备选污染源清单，将污染源相关数据作为河网水动力水质模型的输入条件，模拟各控制污染物沿河段距离动态变化的浓度值；

S4、根据备选污染源清单，明确影响目标水质断面的污染源集合，基于排除法，在河网水动力水质模型中设定污染源全排和依次停排的情况，模拟得到不同情况下控制污染物的浓度及目标水质断面的水流量，进一步计算每个备选污染源对目标水质断面的贡献率，并根据贡献率大小进行排序，从而确定水质断面的主要污染源；

所述步骤S4中计算每个备选污染源对目标水质断面的贡献率的具体过程为：

模拟污染源全部排放时，向河网水动力水质模型中输入所有污染源的数据，由此得到所有污染源的各主要控制污染物在目标水质断面处的浓度模拟值，以及目标水质断面处的水流量模拟值；模拟污染源依次停排时，不向河网水动力水质模型中输入被停排的污染源的相关数据，由此得到未被停排的污染源的各主要控制污染物在目标水质断面处的浓度模拟值，以及此时目标水质断面处的水流量模拟值；多次进行上述模拟试验，最后计算污染源对于水质断面的贡献率：

式中，p_i表示被停排的污染源i对目标水质断面的贡献率，i表示被停排的污染源的编号，t表示按时间进行的试验次数，C_0t表示第t次试验时污染源全排的情况下所有污染源的主要控制污染物在目标水质断面处的浓度之和，Q_0t表示第t次试验时污染源全排的情况下目标水质断面处的水流量，C_st表示第t次试验时关停污染源i的情况下未被停排的污染源的主要控制污染物在目标水质断面处的浓度之和；Q_st表示第t次试验时关停污染源i的情况下目标水质断面处的水流量；

所述步骤S1的具体地过程为：

S11、根据河流方向，将河流重要干流支流进行概化，确定干流、一级支流和二级支流，并以水力学特征为基础，在污染源排放口下游污水与河水基本混匀处以及特殊要求地区的河段上设置水质断面，以水质断面、支流汇入干流口和特殊要求地区为节点，对河段进行划分并排序；

S12、确定河网区所属的水质控制功能区，明确各水质断面的水质控制目标，并将所述水质控制目标作为水质评价指标，确定各水质断面的控制污染物，构建包含水质断面所处河段序号、距离河段起始点的距离、控制污染物类型的水质断面清单；

S13、基于实测统计数据，考虑点源和非点源污染源，根据污染源对水质断面的影响大小，筛除部分污染源，并将相邻密集的点源污染源概化为集中的污染源，非点源污染源则记录下所影响河段的起始点与终点距离河段上边界断面的位置，参照由水质控制目标所确定的各水质断面的控制污染物，构建包含污染源所处河段序号、污染源类型、控制污染物的排放浓度、排放流量、排放口距离河段上边界的距离、属于出流还是入流污染源的污染源排放清单；

S14、比对所述水质断面清单与所述污染源清单中所记录的水质断面和污染源所处河段及距离，确定污染源排污口和水质断面之间的上下游关系，确定每个水质断面的备选污染源，构建水质断面的备选污染源清单；

所述步骤S2的具体过程为：

根据所要研究的河网区的特征，选用合适的水质模型，将水质断面概化成底坡水平、断面矩形或梯形的简单几何形状，通过对实测水质断面的水深与过水面积关系进行回归分析，得到概化断面的下底宽和边坡系数，使得概化后的水质断面与实测水质断面保持一致的水位关系；

确定河段的河道坡度、河段水力参数、河段水质参数、流域水文资料、河道上游、下游边界条件和初始条件、断面节点的情况、流域气象数据以及水质监测站点所测污染物数据，按照所选用的水质模型的操作规则，将数据输入模型中，选择合适的计算步长进行计算，将所模拟的结果数据与水质监测站点的数据曲线进行比较，分析相对误差，若所述相对误差在允许范围之内，则说明建模成功，否则，根据实测数据修改模型参数，重新验证，进而得到适用于所要研究的河网区的水质模型；

所述河段水力参数根据河段河道实际情况，选择性使用水力特征值的计算方法；所述河段水质参数根据河网区所处环境、污染物及污染物变化程度，通过实测统计数据或经验公式法获得河段的沉降、水解速率；

所述步骤S3中，所述污染源相关数据包括每个污染源的流量、主要控制污染物的浓度、所处河段序号以及距离河段上边界距离。

Images