CN112966902A - 一种兼顾区域公平和污染源治理差异的水污染负荷分配方法 - Google Patents

一种兼顾区域公平和污染源治理差异的水污染负荷分配方法 Download PDF

Info

Publication number
CN112966902A
CN112966902A CN202110156017.6A CN202110156017A CN112966902A CN 112966902 A CN112966902 A CN 112966902A CN 202110156017 A CN202110156017 A CN 202110156017A CN 112966902 A CN112966902 A CN 112966902A
Authority
CN
China
Prior art keywords
water
pollutant
treatment cost
pollution
ton
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202110156017.6A
Other languages
English (en)
Other versions
CN112966902B (zh
Inventor
黄晓敏
谢红忠
徐成剑
张仲伟
马俊超
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Changjiang Institute of Survey Planning Design and Research Co Ltd
Original Assignee
Changjiang Institute of Survey Planning Design and Research Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Changjiang Institute of Survey Planning Design and Research Co Ltd filed Critical Changjiang Institute of Survey Planning Design and Research Co Ltd
Priority to CN202110156017.6A priority Critical patent/CN112966902B/zh
Publication of CN112966902A publication Critical patent/CN112966902A/zh
Application granted granted Critical
Publication of CN112966902B publication Critical patent/CN112966902B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063—Operations research, analysis or management
    • G06Q10/0637—Strategic management or analysis, e.g. setting a goal or target of an organisation; Planning actions based on goals; Analysis or evaluation of effectiveness of goals
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063—Operations research, analysis or management
    • G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
    • G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/10—Services
    • G06Q50/26—Government or public services
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A20/00—Water conservation; Efficient water supply; Efficient water use
    • Y02A20/152—Water filtration

Landscapes

  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Engineering & Computer Science (AREA)
  • Strategic Management (AREA)
  • Educational Administration (AREA)
  • Economics (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Tourism & Hospitality (AREA)
  • Development Economics (AREA)
  • Marketing (AREA)
  • Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Game Theory and Decision Science (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Primary Health Care (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Sewage (AREA)

Abstract

本发明提供一种兼顾区域公平和污染源治理差异的水污染负荷分配方法，包括选取评价指标，计算各指标基尼系数及其权重；计算现状综合基尼系数；设置总量控制约束；设置公平性约束；设置削减率约束；确定各行政区涉及工业废水的吨水治理费用；确定各行政区城镇生活污水吨水治理费用；确定市政提标改造的吨水治理费用；确定各行政区工业废水、城镇生活污水、市政提标改造的单位污染负荷治理费用；计算综合治理费用；水污染负荷分配问题进行求解。本发明能够科学合理地指导各区域和各类污染源中污染物削减顺序及削减量方案制定。

Description

一种兼顾区域公平和污染源治理差异的水污染负荷分配方法
技术领域
本发明涉及水环境治理领域，具体地指一种兼顾区域公平和污染源治理差异的水污染负荷分配方法。
背景技术
对于涉及多区域、多类型污染源的水体，如何确定各区域和各污染源的削减率是水污染防治工作亟待解决的重要内容之一。近年来众多学者陆续尝试提出兼顾公平性和经济最优性的水污染负荷分配方法。李如忠等基于以综合基尼系数作为公平性目标，以污染物削减费用最小作为经济最优性目标，构建了基于双目标的水污染负荷分配方法；该方法将每个子区域“虚拟”为一座污水处理厂，每个子区域内水污染负荷的削减均由“虚拟”污水处理厂承担，削减费用即为“虚拟”污水处理厂的处理费用；该方法提出了一种可行的概化污染物削减费用的思路，但在污染物削减费用函数中并未考虑不同污染源废污水特点和治理费用的差异。罗军刚等提出了一种双层多目标优化水污染负荷分配方法，其中，上层基于环境基尼系数最小和单位污染物排放成本最低对水污染物在区域间进行分配，下层基于工业产值最大和削减率不均最小对水污染物在排污单位间进行分配；该方法较好地考虑了区域及排污口两个等级的水污染物分配，以污水处理成本和水环境保护税税额来表征污染物排放成本，同样未细化不同污染源治理成本差异。
发明内容
为了解决上述背景技术提出的技术问题，本发明提供一种兼顾区域公平和污染源治理差异的水污染负荷分配方法。
本发明的技术方案：
一种兼顾区域公平和污染源治理差异的水污染负荷分配方法，包括以下步骤：
步骤一.选取评价指标，基于各行政区评价指标及现状污染负荷，并计算各指标基尼系数及其权重；
步骤二.根据各指标基尼系数及其权重，计算现状综合基尼系数；
步骤三.设置总量控制约束，根据各水功能区限制排污总量作为总量控制目标，并据此确定每种污染物的区域和污染源削减率；
步骤四.设置公平性约束，候选分配方案的综合基尼系数不大于现状综合基尼系数，公平性不会变差；
步骤五.设置削减率约束，根据各水功能区的总量控制目标，首先确定水功能区污染物总削减率，同时结合各污染源治理水平的差异，在技术可行和经济合理条件下，为各污染源设定污染物入河量削减率上限和下限；
步骤六.确定各行政区涉及工业废水的吨水治理费用；
步骤七.确定各行政区城镇生活污水吨水治理费用；
步骤八.基于各行政区市政提标需求，确定市政提标改造的吨水治理费用；
步骤九.确定各行政区工业废水、城镇生活污水、市政提标改造的单位污染负荷治理费用；
步骤十.根据各行政区现状污染负荷和最大允许排放负荷，计算综合治理费用；
步骤十一.基于多目标优化算法对以综合基尼系数最小和综合治理成本最低为目标函数、以总量控制、公平性不变差、削减率技术可行经济合理为约束条件的水污染负荷分配问题进行求解。
所述步骤一中计算各指标基尼系数及其权重的具体方法为：
Figure BDA0002934722260000021
Figure BDA0002934722260000022
Figure BDA0002934722260000031
式中，Gj为基于某个指标j的基尼系数；Xj(i)为指标j的累计百分比；Mj(i)为第i个分区内j指标值；Yj(i)为指标j的污染物入河量累计百分比；Wj(i)为第i个分区内j污染物的入河量；n为分配区域的个数；当i＝1时，(Xj-1，Yj-1)视为(0，0)，
yij＝xi/zij
Figure BDA0002934722260000032
Figure BDA0002934722260000033
Figure BDA0002934722260000034
式中：xi为第j个分区内分配的最大允许污染物入河量；zij表示第i个分区内第j个指标的实际值；yij为第i个分区内第j个指标的单位负荷污染物量；pij为第j个指标、第i个分区的单位负荷污染物量所占权重；ej为第j个指标的单位负荷污染物量的信息熵；wj为第j个指标的权重；n为分区的个数；m为评价指标的个数。
所述步骤二中计算现状综合基尼系数具体方法为：
Figure BDA0002934722260000035
式中，G0为现状综合基尼系数；Gj0为基于某个指标j的现状基尼系数；wj0为第j个指标的权重；m为指标的个数。
所述步骤三中确定每种污染物的区域和污染源削减率的方法为：
Figure BDA0002934722260000036
Figure BDA0002934722260000041
WDi＝(1-qdi)×PDi
Figure BDA0002934722260000042
式中：n为分区数量；np为污染源种类；W为某污染物的限制排污总量；WDi为第i个分区内污染物最高允许入河量；
Figure BDA0002934722260000043
为第i个分区内第k种污染源的污染物最高允许入河量；PDi为第i个分区内污染物现状入河量；
Figure BDA0002934722260000044
为第i个分区内第k种污染源的现状污染物入河量；qdi为第i个分区内污染物削减率；
Figure BDA0002934722260000045
第i个分区内第k种污染源的污染物削减率。
所述步骤四中公平性约束为：
G≤G0
式中：G为候选分配方案综合基尼系数；G0为现状综合基尼系数；
所述步骤五中削减率约束为：
Figure BDA0002934722260000046
Figure BDA0002934722260000047
式中：P为水功能区现状入河负荷；W为水功能区限制排污总量；r为水功能区污染负荷削减率；i为分区编号；
Figure BDA0002934722260000048
为第i个分区内第k种污染源的现状污染物入河量；
Figure BDA0002934722260000049
为第i个分区内第k种污染源的污染物最高允许入河量；
Figure BDA00029347222600000410
分别为第i个分区内第k种污染源污染物削减比例的可行上下限。
所述步骤六中确定各行政区涉及工业废水的吨水治理费用方法为：
NC工业＝(Y2+M×Y3)×(1+i′)n′/10000
式中：NC工业为工业废水吨水治理费用(万元/t)；Y2为运行费用(元/m3水)，包括能耗费、折旧费、原材料费、维修费、工资、管理费和其他费用；Y3为电耗费用(度/m3水)，即设施运行过程中的用电量；Q为设计处理水量(m3/d)；M为当地电费单价(元/度)；i′为物价上涨率5％；n′为治理费用调研的基准年至今年限。
所述步骤七中确定各行政区城镇生活污水吨水治理费用方法为：
Figure BDA0002934722260000051
NC城镇＝Y×(1+i′)n′/Q×24×365×nd
式中：Y为城镇污水处理总费用(万元)；Q为处理水量(t/h)；η为处理效率(％)；k1、k2、k3、k4为参数；NC城镇为城镇污水的吨水治理费用(万元/t)；i′为物价上涨率5％；n′为治理费用调研的基准年至今年限；Q为处理水量(t/h)；nd为折算年数。
所述步骤八中确定市政提标改造的吨水治理费用方法为：
Figure BDA0002934722260000052
Figure BDA0002934722260000053
Figure BDA0002934722260000054
式中：
Figure BDA0002934722260000055
为城镇污水处理厂污染物排放标准中三级标准提标为二级标准的吨水治理费用(万元/t)；
Figure BDA0002934722260000056
为城镇污水处理厂污染物排放标准中二级标准提标为一级标准的吨水治理费用(万元/t)；
Figure BDA0002934722260000057
为城镇污水处理厂污染物排放标准中三级标准提标为一级标准的吨水治理费用(万元/t)；W为实际处理规模(万t)；D为设计规模(万t)；a、b、c、d、l、g、k、h为参数。
所述步骤九中确定各行政区工业废水、城镇生活污水、市政提标改造的单位污染负荷治理费用计算方法为：
Figure BDA0002934722260000061
Figure BDA0002934722260000062
Figure BDA0002934722260000063
式中：ηl废污水处理设施对第l种污染物的处理效益；Il废污水处理设施第l种污染物的进水浓度(mg/L)；El废污水处理设施第l种污染物的出水浓度(mg/L)；Sl第l种污染物的排放标准(mg/L)；βl废污水处理设施第l种污染物的处理费用系数；Q为设计处理水量(m3/d)；Cl废污水第l种污染物的单位负荷处理费用，(元/t)；NC废污水的吨水处理费用，(元/t)；np污染源种类。
所述步骤十中计算综合治理费用的方法为：
Figure BDA0002934722260000064
式中：C为污染物治理费用(万元)；n为分区个数；s为污染源数量；
Figure BDA0002934722260000065
为第i个分区第k种污染源的污染物现状入河量(t)；
Figure BDA0002934722260000066
为第i个分区第k种污染源的污染物最高允许入河量(t)；
Figure BDA0002934722260000067
为第i个分区第k种污染源的单位污染负荷治理费用(元/t)。
与现有技术相比，本发明的有益效果是：针对多行政区、多污染源的超标水体，提供兼顾区域公平和削减效率的一种水污染负荷分配方法，能够科学合理地指导各区域和各类污染源中污染物削减顺序及削减量方案制定。
附图说明
图1为本发明的方法流程图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
一种兼顾区域公平和污染源治理差异的水污染负荷分配方法，包括以下步骤：
步骤一.选取评价指标，基于各行政区评价指标及现状污染负荷，并计算各指标基尼系数及其权重；
步骤二.根据各指标基尼系数及其权重，计算现状综合基尼系数；
步骤三.设置总量控制约束，根据各水功能区限制排污总量作为总量控制目标，并据此确定每种污染物的区域和污染源削减率；
步骤四.设置公平性约束，候选分配方案的综合基尼系数不大于现状综合基尼系数，公平性不会变差；
步骤五.设置削减率约束，根据各水功能区的总量控制目标，首先确定水功能区污染物总削减率，同时结合各污染源治理水平的差异，在技术可行和经济合理条件下，为各污染源设定污染物入河量削减率上限和下限；
步骤六.确定各行政区涉及工业废水的吨水治理费用；
步骤七.确定各行政区城镇生活污水吨水治理费用；
步骤八.基于各行政区市政提标需求，确定市政提标改造的吨水治理费用；
步骤九.确定各行政区工业废水、城镇生活污水、市政提标改造的单位污染负荷治理费用；
步骤十.根据各行政区现状污染负荷和最大允许排放负荷，计算综合治理费用；
步骤十一.基于多目标优化算法对以综合基尼系数最小和综合治理成本最低为目标函数、以总量控制、公平性不变差、削减率技术可行经济合理为约束条件的水污染负荷分配问题进行求解。
<实施例>
以该方法求解某流域水污染负荷分配问题为例进行详细说明：
1)选取常住人口和生产总值作为评价指标，流域内常住人口和生产总值分别为202813人和22167万元，COD现状入河总量为580.72t，流域内涉及三个行政区，常住人口为
Figure BDA0002934722260000081
人，生产总值为
Figure BDA0002934722260000082
万元，COD入河量为
Figure BDA0002934722260000083
t，污染源类型涉及工业废水、城镇生活污水和市政尾水，对应COD入河量为
Figure BDA0002934722260000084
t，据此计算常住人口现状基尼系数及对应权重为0.3252和0.3039，生产总值现状基尼系数及对应权重为0.3528和0.6961。
2)根据常住人口、生产总值的现状基尼系数和对应权重，计算现状综合基尼系数为0.3444。
G0＝0.3252×0.3039+0.3528×0.6961＝0.3444
3)根据水功能区COD限制排污总量，确定总量控制目标为213.83t。设置总量约束条件为：
213.83＝(1-qd1)×199.83+(1-qd2)×274.33+(1-qd3)×106.57
对于行政区一：
Figure BDA0002934722260000085
对于行政区二：
Figure BDA0002934722260000086
对于行政区三：
Figure BDA0002934722260000087
4)候选分配方案的综合基尼系数不大于现状综合基尼系数，公平性不会变差。设置公平性约束为：
G≤0.3444
5)根据各水功能区的总量控制目标，首先确定水功能区污染物总削减率。同时结合各污染源治理水平的差异，在技术可行和经济合理条件下，为各污染源设定污染物入河量削减率上限和下限。设置削减率约束为：
Figure BDA0002934722260000091
对于行政区一：
城镇生活污水削减率
Figure BDA0002934722260000092
市政提标改造削减率
Figure BDA0002934722260000093
对于行政区二：
城镇生活污水削减率
Figure BDA0002934722260000094
市政提标改造削减率
Figure BDA0002934722260000095
对于行政区三：
工业废水治理削减率为
Figure BDA0002934722260000096
城镇生活污水削减率
Figure BDA0002934722260000097
6)确定工业废水的吨水治理费用。水功能区汇水范围内仅行政区三设有食品行业排口，其吨水治理费用为：
NC工业＝(0.17+0.74×0.04)×(1+0.05)33′/10000＝1.00(元/t)
7)确定城镇生活污水吨水治理费用。按照《城镇污水处理厂污染物排放标准》(GB18918-2002)一级A标准，其吨水治理费用为：
行政区一：
Y1＝9×112.660.657+22×112.660.657×0.91.7＝610.44(万元)
NC城镇1＝610.44×(1+0.05)33/112.66×24×365×20＝1.55(元/t)
行政区二：
Y2＝9×197.660.657+22×197.660.657×0.91.7＝883.17(万元)
NC城镇2＝883.17×(1+0.05)33/197.66×24×365×20＝1.28(元/t)
行政区三：
Y3＝9×39.550.657+22×39.550.657×0.91.7＝306.89(万元)
NC城镇3＝306.89×(1+0.05)33/39.55×24×365×20＝2.22(元/t)
8)确定市政提标改造的吨水治理费用。行政区一和行政区二的城镇污水处理厂尾水排放标准由《城镇污水处理厂污染物排放标准》(GB18918-2002)二级标准提升为《城镇污水处理厂污染物排放标准》(GB18918-2002)一级A标准，其吨水治理费用为：
对于行政区一：
NC市政提标1＝(e1.093×4500.830-1-e1.093-0.315×4500.830-1)+(e-0.305×8000.889-1-e-0.305-0.234×8000.889-1)/365×20＝0.39(元/t)
对于行政区二：
NC市政提标2＝(e1.093×5000.830-1-e1.093-0.315×5000.830-1)+(e-0.305×5000.889-1-e-0.305-0.234×5000.889-1)/365×20＝0.39(元/t)
9)确定各行政区工业废水、城镇生活污水、市政提标改造的单位污染负荷治理费用。污染物的处理费用系数主要考虑COD和NH3-N。
对于行政区一：
Figure BDA0002934722260000101
Figure BDA0002934722260000102
Figure BDA0002934722260000103
Figure BDA0002934722260000104
Figure BDA0002934722260000111
Figure BDA0002934722260000112
Figure BDA0002934722260000113
Figure BDA00029347222600001118
对于行政区二：
Figure BDA0002934722260000115
Figure BDA0002934722260000116
Figure BDA0002934722260000117
Figure BDA00029347222600001119
Figure BDA0002934722260000119
Figure BDA00029347222600001110
Figure BDA00029347222600001111
Figure BDA00029347222600001120
对于行政区三：
Figure BDA00029347222600001113
Figure BDA00029347222600001114
Figure BDA00029347222600001115
Figure BDA00029347222600001121
Figure BDA00029347222600001117
Figure BDA0002934722260000121
Figure BDA0002934722260000122
Figure BDA00029347222600001213
10)根据各行政区现状污染负荷和最大允许排放负荷，计算综合治理费用。
对于行政区一：
Figure BDA00029347222600001210
对于行政区二：
Figure BDA00029347222600001211
对于行政区三：
Figure BDA00029347222600001212
Figure BDA0002934722260000127
11)基于MOPSO算法对以综合基尼系数最小和综合治理成本最低为目标函数、以总量控制、公平性不变差、削减率技术可行经济合理为约束条件的水污染负荷分配问题进行求解。设置MOPSO初始种群为100，进化代数为50，变异概率为1。兼顾区域公平与污染源治理差异的COD多目标pareto解集详见表1,用户可根据实际需求选取意向分配方案。优化后分配方案的综合基尼系数均小于现状综合基尼系数，各污染源COD削减比例均符合削减率上限和下限要求，三个行政区COD入河量符合总量控制目标要求，故该发明提出的兼顾区域公平与污染源治理差异的水污染负荷分配方法是可行的。
表1兼顾区域公平与污染源治理差异的水污染负荷分配解集
Figure BDA0002934722260000129
Figure BDA0002934722260000131
以上实施例仅仅是针对本发明技术方案所做的举例说明。本发明所涉及的水文模型参数时变形式构造方法并不限定于在以上实施例中所描述的内容，而是以权利要求所限定的范围为准。本发明所述领域技术人员在该实施例的基础上所做的任何修改或补充或等效替换，均在本发明所要求保护范围内。

Claims (10)

1.一种兼顾区域公平和污染源治理差异的水污染负荷分配方法，其特征在于，包括以下步骤：
步骤一.选取评价指标，基于各行政区评价指标及现状污染负荷，并计算各指标基尼系数及其权重；
步骤二.根据各指标基尼系数及其权重，计算现状综合基尼系数；
步骤三.设置总量控制约束，根据各水功能区限制排污总量作为总量控制目标，并据此确定每种污染物的区域和污染源削减率；
步骤四.设置公平性约束，候选分配方案的综合基尼系数不大于现状综合基尼系数，公平性不会变差；
步骤五.设置削减率约束，根据各水功能区的总量控制目标，首先确定水功能区污染物总削减率，同时结合各污染源治理水平的差异，在技术可行和经济合理条件下，为各污染源设定污染物入河量削减率上限和下限；
步骤六.确定各行政区涉及工业废水的吨水治理费用；
步骤七.确定各行政区城镇生活污水吨水治理费用；
步骤八.基于各行政区市政提标需求，确定市政提标改造的吨水治理费用；
步骤九.确定各行政区工业废水、城镇生活污水、市政提标改造的单位污染负荷治理费用；
步骤十.根据各行政区现状污染负荷和最大允许排放负荷，计算综合治理费用；
步骤十一.基于多目标优化算法对以综合基尼系数最小和综合治理成本最低为目标函数、以总量控制、公平性不变差、削减率技术可行经济合理为约束条件的水污染负荷分配问题进行求解。
2.根据权利要求1所述的一种兼顾区域公平和污染源治理差异的水污染负荷分配方法，其特征在于，所述步骤一中计算各指标基尼系数及其权重的具体方法为：
Figure FDA0002934722250000021
Figure FDA0002934722250000022
Figure FDA0002934722250000023
式中，Gj为基于某个指标j的基尼系数；Xj(i)为指标j的累计百分比；Mj(i)为第i个分区内j指标值；Yj(i)为指标j的污染物入河量累计百分比；Wj(i)为第i个分区内j污染物的入河量；n为分配区域的个数；当i＝1时，(Xj-1，Yj-1)视为(0，0)，
yij＝xi/zij
Figure FDA0002934722250000024
Figure FDA0002934722250000025
Figure FDA0002934722250000026
式中：xi为第j个分区内分配的最大允许污染物入河量；zij表示第i个分区内第j个指标的实际值；yij为第i个分区内第j个指标的单位负荷污染物量；pij为第j个指标、第i个分区的单位负荷污染物量所占权重；ej为第j个指标的单位负荷污染物量的信息熵；wj为第j个指标的权重；n为分区的个数；m为评价指标的个数。
3.根据权利要求2所述的一种兼顾区域公平和污染源治理差异的水污染负荷分配方法，其特征在于，所述步骤二中计算现状综合基尼系数具体方法为：
Figure FDA0002934722250000027
式中，G0为现状综合基尼系数；Gj0为基于某个指标j的现状基尼系数；wj0为第j个指标的权重；m为指标的个数。
4.根据权利要求3所述的一种兼顾区域公平和污染源治理差异的水污染负荷分配方法，其特征在于，所述步骤三中确定每种污染物的区域和污染源削减率的方法为：
Figure FDA0002934722250000031
Figure FDA0002934722250000032
WDi＝(1-qdi)×PDi
Figure FDA0002934722250000033
式中：n为分区数量；np为污染源种类；W为某污染物的限制排污总量；WDi为第i个分区内污染物最高允许入河量；
Figure FDA0002934722250000034
为第i个分区内第k种污染源的污染物最高允许入河量；PDi为第i个分区内污染物现状入河量；
Figure FDA0002934722250000035
为第i个分区内第k种污染源的现状污染物入河量；qdi为第i个分区内污染物削减率；
Figure FDA0002934722250000036
第i个分区内第k种污染源的污染物削减率。
5.根据权利要求4所述的一种兼顾区域公平和污染源治理差异的水污染负荷分配方法，其特征在于，所述步骤四中公平性约束为：
G≤G0
式中：G为候选分配方案综合基尼系数；G0为现状综合基尼系数；
所述步骤五中削减率约束为：
Figure FDA0002934722250000037
Figure FDA0002934722250000038
式中：P为水功能区现状入河负荷；W为水功能区限制排污总量；r为水功能区污染负荷削减率；i为分区编号；
Figure FDA0002934722250000041
为第i个分区内第k种污染源的现状污染物入河量；
Figure FDA0002934722250000042
为第i个分区内第k种污染源的污染物最高允许入河量；
Figure FDA0002934722250000043
分别为第i个分区内第k种污染源污染物削减比例的可行上下限。
6.根据权利要求5所述的一种兼顾区域公平和污染源治理差异的水污染负荷分配方法，其特征在于，所述步骤六中确定各行政区涉及工业废水的吨水治理费用方法为：
NC工业＝(Y2+M×Y3)×(1+i′)n′/10000
式中：NC工业为工业废水吨水治理费用，单位为万元/t，即万元每吨；Y2为运行费用，单位为元/m3水，包括能耗费、折旧费、原材料费、维修费、工资、管理费和其他费用；Y3为电耗费用，单位为度/m3水，即设施运行过程中的用电量；Q为设计处理水量，单位为m3/d即立方米每天；M为当地电费单价，单位为元/度；i′为物价上涨率5％；n′为治理费用调研的基准年至今年限。
7.根据权利要求6所述的一种兼顾区域公平和污染源治理差异的水污染负荷分配方法，其特征在于，所述步骤七中确定各行政区城镇生活污水吨水治理费用方法为：
Figure FDA0002934722250000044
NC城镇＝Y×(1+i′)n′/Q×24×365×nd
式中：Y为城镇污水处理总费用，单位为万元；Q为处理水量，单位为t/h，即吨每小时；η为处理效率；k1、k2、k3、k4为参数；NC城镇为城镇污水的吨水治理费用，单位为万元/t，即万元每吨；i′为物价上涨率5％；n′为治理费用调研的基准年至今年限；nd为折算年数。
8.根据权利要求7所述的一种兼顾区域公平和污染源治理差异的水污染负荷分配方法，其特征在于，所述步骤八中确定市政提标改造的吨水治理费用方法为：
Figure FDA0002934722250000051
Figure FDA0002934722250000052
Figure FDA0002934722250000053
式中：
Figure FDA0002934722250000054
为城镇污水处理厂污染物排放标准中三级标准提标为二级标准的吨水治理费用，单位为万元/t，即万元每吨；
Figure FDA0002934722250000055
为城镇污水处理厂污染物排放标准中二级标准提标为一级标准的吨水治理费用，单位为万元/t，即万元每吨；
Figure FDA0002934722250000056
为城镇污水处理厂污染物排放标准中三级标准提标为一级标准的吨水治理费用，单位为万元/t，即万元每吨；W为实际处理规模，单位为万t，即万吨；D为设计规模，单位为万t，即万吨；a、b、c、d、l、g、k、h为参数。
9.根据权利要求8所述的一种兼顾区域公平和污染源治理差异的水污染负荷分配方法，其特征在于，所述步骤九中确定各行政区工业废水、城镇生活污水、市政提标改造的单位污染负荷治理费用计算方法为：
Figure FDA0002934722250000057
Figure FDA0002934722250000058
Figure FDA0002934722250000059
式中：ηl废污水处理设施对第l种污染物的处理效益；Il废污水处理设施第l种污染物的进水浓度，单位为mg/L；El废污水处理设施第l种污染物的出水浓度，单位为mg/L；Sl第l种污染物的排放标准，单位为mg/L；βl废污水处理设施第l种污染物的处理费用系数；Q为设计处理水量，单位为m3/d，即立方米每天；Cl废污水第l种污染物的单位负荷处理费用，单位为元/t，即元每吨；NC废污水的吨水处理费用，单位为元/t，即元每吨；np污染源种类。
10.根据权利要求9所述的一种兼顾区域公平和污染源治理差异的水污染负荷分配方法，其特征在于，所述步骤十中计算综合治理费用的方法为：
Figure FDA0002934722250000061
式中：C为污染物治理费用单位为万元；n为分区个数；s为污染源数量；
Figure FDA0002934722250000062
为第i个分区第k种污染源的污染物现状入河量，单位为吨；
Figure FDA0002934722250000063
为第i个分区第k种污染源的污染物最高允许入河量，单位为吨；
Figure FDA0002934722250000064
为第i个分区第k种污染源的单位污染负荷治理费用，单位为元/t，即元每吨。
CN202110156017.6A 2021-02-04 2021-02-04 一种兼顾区域公平和污染源治理差异的水污染负荷分配方法 Active CN112966902B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110156017.6A CN112966902B (zh) 2021-02-04 2021-02-04 一种兼顾区域公平和污染源治理差异的水污染负荷分配方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110156017.6A CN112966902B (zh) 2021-02-04 2021-02-04 一种兼顾区域公平和污染源治理差异的水污染负荷分配方法

Publications (2)

Publication Number Publication Date
CN112966902A true CN112966902A (zh) 2021-06-15
CN112966902B CN112966902B (zh) 2022-05-17

Family

ID=76273921

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110156017.6A Active CN112966902B (zh) 2021-02-04 2021-02-04 一种兼顾区域公平和污染源治理差异的水污染负荷分配方法

Country Status (1)

Country Link
CN (1) CN112966902B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115578232A (zh) * 2022-11-11 2023-01-06 北京华科仪科技股份有限公司 一种基于水质分析的水污染治理方法及系统

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005030912A (ja) * 2003-07-14 2005-02-03 Penta Ocean Constr Co Ltd 水質環境の監視方法及び監視システム
JP2007065966A (ja) * 2005-08-31 2007-03-15 Yamatake Corp 河川汚濁負荷源推定装置、方法、およびプログラム
CN110210710A (zh) * 2019-05-06 2019-09-06 河海大学 一种基于负荷均衡的水资源承载能力量化方法
CN110619435A (zh) * 2019-09-18 2019-12-27 西安理工大学 一种双层多目标优化水污染负荷分配方法
CN110969346A (zh) * 2019-11-22 2020-04-07 江苏省环境科学研究院 基于指标筛选的流域水生态功能分区治理需求评价方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005030912A (ja) * 2003-07-14 2005-02-03 Penta Ocean Constr Co Ltd 水質環境の監視方法及び監視システム
JP2007065966A (ja) * 2005-08-31 2007-03-15 Yamatake Corp 河川汚濁負荷源推定装置、方法、およびプログラム
CN110210710A (zh) * 2019-05-06 2019-09-06 河海大学 一种基于负荷均衡的水资源承载能力量化方法
CN110619435A (zh) * 2019-09-18 2019-12-27 西安理工大学 一种双层多目标优化水污染负荷分配方法
CN110969346A (zh) * 2019-11-22 2020-04-07 江苏省环境科学研究院 基于指标筛选的流域水生态功能分区治理需求评价方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
HAI, YE ETAL.: "Polluted river water treatment using multi-soil-layering system (MSL) in different surface loads", 《2013 INTERNATIONAL CONFERENCE ON APPLICATION OF MATERIALS SCIENCE AND ENVIRONMENTAL MATERIALS, AMSEM 2013》 *
庞莹莹: "湖泊流域纳污能力及污染负荷分配研究", 《中国优秀博硕士学位论文全文数据库(硕士)工程科技Ⅰ辑》 *
李如忠等: "基于基尼系数的水污染负荷分配模糊优化决策模型", 《环境科学学报》 *
王媛等: "基尼系数法在水污染物总量区域分配中的应用", 《中国人口.资源与环境》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115578232A (zh) * 2022-11-11 2023-01-06 北京华科仪科技股份有限公司 一种基于水质分析的水污染治理方法及系统

Also Published As

Publication number Publication date
CN112966902B (zh) 2022-05-17

Similar Documents

Publication Publication Date Title
CN105297827A (zh) 一种考虑多用户需水和多水源供水的水资源配置方法
Ahluwalia Challenges of urbanisation in India
JP5918672B2 (ja) 上下水道管理システム及び方法
CN113011662B (zh) 一种区域企厂网河一体化联合优化调度的系统及方法
WO2021238170A1 (zh) 二元水权核算交易方法及系统
CN113159387B (zh) 一种工业园区污水近零排放的水资源优化配置系统及方法
CN109389251A (zh) 基于控制断面水质达标的污染物总量优化分配方法
Zhang et al. Effect of the tap water supply system on China's economy and energy consumption, and its emissions’ impact
CN105913146A (zh) 南方湿润地区的水资源优化配置系统
Zeng et al. A developed fuzzy-stochastic optimization for coordinating human activity and eco-environmental protection in a regional wetland ecosystem under uncertainties
KR20100095819A (ko) 유역통합관리를 위한 수질오염 총량관리방법
CN112966902B (zh) 一种兼顾区域公平和污染源治理差异的水污染负荷分配方法
Ai et al. Optimization of ecological reservoir operation rules for a northern river in China: Balancing ecological and socio-economic water use
Baş et al. Environmental and techno-economic analysis of the integration of biogas and solar power systems into urban wastewater treatment plants
CN105906032A (zh) 污水处理厂拟人化经验管理控制系统及方法
CN117114918A (zh) 一种基于区域生态需水与多水源供水的水资源调配方法
Zhao et al. Allocation of SO2 emission rights in city agglomerations considering cross-border transmission of pollutants: A new network DEA model
CN116484647B (zh) 面向统筹协调的分布式水资源配置方法及系统
Liang et al. Cost benefit analysis of centralized wastewater reuse systems
Zeng et al. Strategic zoning for urban wastewater reuse in China
Dong et al. Analysis of the ultimate water resources carrying capacity in Yancheng, China
CN114035619B (zh) 一种家庭用水控制方法和系统
Ma et al. Study on water resources optimal regulation of inter-basin water transfer project
Zhang Calculation model and management system of WRCC based on grey water footprint calculation
CN117829559B (zh) 一种多水源到多用户的水资源一体化调配方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant