基于污染税的流域污染负荷分配二层决策优化方法
基于污染税的流域污染负荷分配二层决策优化方法
【技术领域】
[0001] 本发明涉及流域水质目标管理，特别涉及基于污染税的流域污染负荷分配决策方 法。
【背景技术】
[0002] 二层决策优化方法是一类针对具有二层递阶结构问题的优化方法。其主要研宄具 有两个层次系统的规划与管理问题。分为上层和下层，上下两层问题都具有各自的决策变 量、约束条件和目标函数。这种方法区别于其他优化方法的一个主要方面是不再是一个决 策者进行规划与决策，而是有上下两层决策者，上层决策者只是通过自己的决策去指导下 层决策者，并不直接干涉下层的决策；而下层的决策将上层的决策作为一个影响条件，其在 自己的可行域进行自由决策，这就导致上层决策者在做任何决策的时候必须考虑到下层决 策者可能做出的反馈，以规避下层决策者的决策可能带来的不利影响。
[0003] 在水污染严重的今天，合理的流域水质管理显得至关重要。目前，为进一步治理各 流域水环境，我国环保部制定了一系列的环境经济政策，包括环境税费、排污权交易等。
[0004] 其中，污染税费为主要方式。其旨在国家以立法的形式明确规定污染者要按其排 放的污染物种类、浓度和数量缴纳一定的费用。通过征收污染税，使经济活动所造成的环境 损失由其本身承担，使外部经济内在化；此外，税金的收入可用作改善环境项目的拨款。污 染税费制定在限制各流域排放者的排污行为的同时，也一定程度的增加了地方环保部门的 水污染治理工作难度。具体表现为：①如何设计一套科学的收费体系，可以使某一流域和地 区各排污者实现特定水质标准所需的成本相同；②如何制定一套合理的流域污染负荷控制 技术，能够指导地方环保部门在污染税费政策下进行水污染防治工作。因此，本发明考虑了 基于污染税的流域污染负荷分配问题，旨在提出一种新的水环境污染分配方案。
[0005] 在过去的几十年中，决策优化技术在水资源管理方面起到了至关重要的作用，主 要相关的技术包括多目标决策优化，多属性决策优化以及多阶段决策优化等，这些决策的 一个主要特点是由一个主体进行决策，没有考虑其他利益相关者的意见及反馈。
【发明内容】
[0006] 本发明所要解决的技术问题，就是提供一种基于污染税的流域污染负荷分配二层 决策优化方法以解决基于污染税的流域污染负荷分配问题。
[0007] 本发明解决所述技术问题，采用的技术方案是，基于污染税的流域污染负荷分配 二层决策优化方法，包括以下步骤：
[0008] 步骤1、系统获取水流量，并将其转换为水流量的确定性变量；
[0009] 步骤2、系统获取河流中溶解氧浓度关系及生化需氧量浓度关系；
[0010] 步骤3、系统根据水流量的确定性变量、溶解氧浓度关系及生化需氧量浓度关系构 建下层模型；
[0011] 步骤4、系统根据污染税收标准制定上层模型；
[0012] 步骤5、系统根据下层模型及上层模型构建基于污染税的流域污染负荷分配二层 决策模型；
[0013] 步骤6、系统对基于污染税的流域污染负荷分配二层决策模型进行求解，得到最优 解。
[0014] 具体的，所述步骤1中，首先，系统采用模糊随机变量对未来时间水流量进行量 化，系统将水流量分为三个不同的级别：低水位（L)，中水位（M)和高水位（H)，服从概率分 布（a，b，c)，三种水位的概率分别可以表示为九，pjPpH，综上，水流量匕为：
[0016] 其次，将水流量匕转换为模糊随机变量，并计算其模糊期望值，具体包括：
[0017] ①假设水流量可以表述成为任意一个模糊随机变量X，含有以下数值：
[0019] 其中，4，4,，表示三个水流量水平，因此其概率分布如下：
[0021] 其中，D表示三种水流水平及其模糊逻辑；
[0022] 假设：0(x)是模糊随机变量x的模糊参数，如果0(x)= Ef [X | P]，V? e [0,1J 存在，则有：
[0025] 那么，模糊随机变量x的模糊期望值可以表示为E'虎]，则有：
[0027] 其中，aE=inf( 0 (x)) 0，bE=inf( 0 (x))i=sup( 0 (x)) "cE=sup( 0 (x)) 0， ElGl表示模糊随机变量x的模糊期望值；
[0028] ②模糊随机水流量被转化成为三角模糊数，系统采用的带乐观-悲观指数的期望 值算子将其转化成为确定性值；
[0029] 首先，定义一个乐观-悲观指数（0<e< 1)用于量化决策者的决策偏好，然后 使用期望值算子计算该三角模糊数的期望值，并将其作为水流量的确定性变量，则有：
[0031] 其中，ElEli]]表示水流量的确定性变量。
[0032] 具体的，所述步骤2中，系统根据Streeter-Phelps模型的结论，河流中生化需氧 量浓度的随时间变化量与其自身浓度成正比，即：
[0034] 其中，Cb表示生化需氧量浓度，kd表示脱氧系数；
[0035] 因此，在一个河段t，用kd(t)表示脱氧系数，d(t)表示该河段长度，v(t)表示该河 段评价水流流速，Cb(t)表示流入该河段水流的生化需氧量浓度，则流出该河段水流的生化 需氧量浓度Cb (t+1)为
该生化需氧量浓度等价于：
[0036]Cb(t+1) =a(t)Cb(t)；
[0037] 系统根据Streeter-Phelps模型的结论，河流中的溶解氧从空气中补充，其补充 浓度的随时间变化量与其氧亏浓度成正比，即：
[0039] 其中，Cb表示生化需氧量浓度，D。表示氧亏浓度，kd表示脱氧系数，ka表示充氧系 数。与生化需氧量情况相似，在一个河段1用ka(t)表示充氧系数，Cb(t)，Cjt)分别表示 流入该河段水流的生化需氧量浓度，Cd。表示该河段水流的饱和溶解氧浓度，则流出该河段 水流的溶解氧浓度C。(t+1)为：
[0041] 定义：
和S⑴=(^。[1-|3(1：)]，该溶解氧浓 度等价于：
[0042]C0(t+1) = |3 (t)C0(t)-y(t)Cb(t)-8 (t)；
[0043]其中，a(t)、Mt)分别可以表示河段t生化需氧量浓度（BOD)、溶氧浓度（DO)的 转换系数，Y(t)、S(t)分别表示影响河段t溶氧浓度的生化需氧量浓度、饱溶氧浓度转换 系数。
[0044] 具体的，所述步骤3中，系统构下层模型包括：
[0045] (1)河流系统状态方程组构建，具体包括：
[0046]①流量状态转移方程：基于质量守恒定律以及流量连续性条件，用Qjt-l)，Qjt) 分别表示在监测点t_l和t的流量，用Qw(t)，Qd(t)分别表示排污控制单元t的取水量和排 污量，则可构建流量状态转移方程如下：
[0047]Qr(t) =Qr(t-1)-Qw(t)+Qd(t)；
[0048] ②BOD浓度状态转移方程：基于质量守恒定律、Streeter-Phelps水质模拟模型以 及流量状态转移方程，用Cb(t-1)，Cb(t)表示监测点t-1和t的B0D浓度，表示排污控 制单元t未处理的污水B0D浓度，以及n(t)为排污控制单元t进行削减污水B0D浓度的 比例，其为下层排污者t的决策变量，则可构建B0D浓度状态转移方程如下：
[0050] 定义变量
，则该状态转移方程等价于：
[0052] ③DO浓度状态转移方程：相似于BOD情形，基于质量守恒定律、Streeter-Phelps 水质模拟模型以及流量状态转移方程，用（；(t-l)，Cjt)表示监测点t-1和t的DO浓度， 表示排污控制单元t未处理的污水D0浓度，则可构建D0浓度状态转移方程如下：
[0054] 同样的定义变量
，则该状态转移方程等价于：
[0056] 另外，考虑河流中饱和溶解氧浓度Cd。的限制，以上D0浓度状态转移方程可改写 为：
[0058] (2)构建初始模型河流状态条件：
[0059] 对于位于整个河流系统中上游第一个排污控制单元，用q，C:分别表示上游水流 B0D和D0浓度，则该控制单元的模型考虑初始河流状态条件如下：
[0061] (3)构建下层模型约束条件：
[0062] ①允许的B0D和D0浓度标准：
[0063] 排污控制单元t排放污水到河流中，不能使得在相应检测点河流的B0D浓度 (Cb(t))超过环保部门所允许的最大B0D浓度标准（Bu(t))，即：
[0064] 0 ^Cb(t) ^Bu(t)；
[0065] 同样，不能使相应检测点河流的DO浓度（C；(t))低于环保部门所允许的最小DO浓 度标准汍⑴），即：
[0066] Dx(t) ^C0(t)；
[0067] ②B0D削减能力约束：
[0068] 每个排污控制单元具有污水B0D浓度削减能力的限制，通常采用二级污水处理工 艺。对于排污控制单元t，第一级处理后的污水B0D浓度削减率为，第二级处理后的污 水B0D浓度削减率为^厂，在实际的污水处理运行中，环保部门要求所以排污者的污水必须 进行第一级处理；因此，排污控制单元t进行污水B0D浓度削减决策，其B0D削减率（n(t)) 不能够超过nr，同时不能低于f，即：
[0070] (4)构建下层模型目标函数：
[0071] 下层决策者的目标为最小的单位污染成本目标函数，污染成本由污水处理成本和 污染税费组成；
[0072] ①污水处理成本：
[0073] 基于B0D削减和污水处理成本的经验曲线，本发明对于排污控制单元的污水处理 费用采用幂函数形式，用kel (t)和ke2 (t)表示排污控制单元t的成本参数1和成本参数2，kra(t)和表示排污控制单元t的污水和BOD削减参数，则污水处理成本C(t)可以表 示为以下方程：
[0075] 其次，参数kel (t)，ke2 (t)，kee⑴和keb⑴数值以污水处理构筑物建设和运营所得 数据，采用最小二乘法进行曲线拟合所得；
[0076] ②污染税费：
[0077] 在每个监测点，地方环保部门制定相应的污水税收标准（Alb(t)，上层模型的 决策变量）和税率（Tb)，如果检 测点t所对应的河流B0D浓度（Cb(t))低于该税收标准 (Cb(t) <Alb(t))，贝1J排污控制单元t不需要征收污染税（P(t));如果检测点t所对应的 河流如0浓度（(；(〇)超过该税收标准（(；(〇>4113(〇)，则排污者〖需征收污染税？(〇 = Qjt)Tb [Cb (t)-Alb (t)]。因此，排放者t的污染税费模型定义为：
[0079] 综上所述，排污控制单元t的目标函数，即最小化的单位污染成本如下：
[0081] 具体的，所述步骤4中，系统构上层模型包括：
[0082] (1)上层模型约束条件构建：
[0083] 对于上层模型的决策者-地方环保部门，在进行相应污染税收标准的制定过程 中，需满足允许的B0D浓度标准约束条件，环保部门所制定的污染税收标准（Alb(t))不能 超过允许的最大B0D浓度标准（Bu(t))，即：
[0084] 0 ^Alb(t) ^Bu(t)；
[0085] (2)上层模型目标函数构建：
[0086] 上层决策者需保证公平的污染负荷分配，因此，本发明采用各排污控制单元的单 位污染成本的方差来量化水污染负荷分配的公平性，用
表示各排污控制单元平 均的单位污染成本，则其最小的方差定义为：
[0088] 具体的，所述步骤5中，系统根据下层模型及上层模型构建的基于污染税的流域 污染负荷分配二层决策模型为：
[0090]
[0091] 具体的，所述步骤6中，假设决策排污者为T个，则在特定上层决策变量下地下层 模型为T个一元的非线性规划模型；首先，系统利用极值算法计算出单一的非线性规划模 型的极值，再通过下层约束条件计算出该单一的非线性规划模型的最优解，进而获得下层T 个一元的非线性规划模型的全部最优解集；利用粒子群算法对上述最优解集进行计算，从 而获得基于污染税的流域污染负荷分配二层决策模型的最优解。
[0092] 具体的，粒子群算法的粒子更新方程如下：
[0093]
[0094] Px(T+1)=Vx(T+1)+?!(T);
[0095] 其中，
，对o -1，2, ? ??，T
，对 1=1，2,...，L;0(T)表示惯性权重，计算过程为e(T) =e(i)+g「e(i)-e(i)i，其中0⑴和 0 (I)表示第1和第I代的惯性权重；Cp和c 别表示个体和全局最优学习因子；rp和r1 表示0到1之间的随机数。
[0096] 具体的，系统采用Boltzmann选择算子代替传统全局最优粒子的更新机制。
[0097] 本发明的有益效果是：本发明针对河流流量这一不确定参数，采用模糊随机数进 行描述，并采用两步法进行处理，即首先将其转化成为模糊数，然后采用模糊期望值算子将 模糊数转化为确定值。针对二层决策模型的求解，根据在污染税政策下水污染负荷分配的 二层决策结构特征以及河流系统中污染物的传递特性，提出了混合嵌套基于多智能体的动 态极值算法的改进粒子群算法（PS0)的方法进行求解，从而得到合理的流域污染税费方案 和污染物削减方案。
【附图说明】
[0098] 图1为本发明基于污染税的流域污染负荷分配二层决策优化方法实施例中基于 污染税的流域污染分配问题的二层决策模型结构；
[0099] 图2为本发明基于污染税的流域污染负荷分配二层决策优化方法实施例中B0D和 D0转化机制的概念图；
[0100] 图3为本发明基于污染税的流域污染负荷分配二层决策优化方法实施例中河流 系统转态转态的概念图；
[0101] 图4为本发明基于污染税的流域污染负荷分配二层决策优化方法实施例中排污 控制单元污水处理分析；
[0102] 图5为本发明基于污染税的流域污染负荷分配二层决策优化方法实施例中针对 基于污染税的流域污染负荷分配问题的二层决策建模技术流程；
[0103] 图6为本发明基于污染税的流域污染负荷分配二层决策优化方法实施例中粒子 结构表示；
[0104] 图7为本发明基于污染税的流域污染负荷分配二层决策优化方法实施例中粒子 性能的混合嵌套计算流程图；
[0105] 图8为本发明基于污染税的流域污染负荷分配二层决策优化方法实施例中更新 全局最优粒子的Boltzmann选择算子流程图；
[0106] 图9为本发明基于污染税的流域污染负荷分配二层决策优化方法实施例中河流 系统图。
【具体实施方式】
[0107] 下面结合附图及实施例详细描述本发明的技术方案：
[0108] 本发明针对基于污染税的流域污染负荷分配问题，提出了基于污染税的流域污染 负荷分配二层决策优化方法，首先，系统获取水流量，并将其转换为水流量的确定性变量； 其次，系统获取河流中溶解氧浓度关系及生化需氧量浓度关系；再其次，系统根据水流量的 确定性变量、溶解氧浓度关系及生化需氧量浓度关系构建下层模型；然后，系统根据污染税 收标准制定上层模型；再然后，系统根据下层模型及上层模型构建基于污染税的流域污染 负荷分配二层决策模型；最后，系统对基于污染税的流域污染负荷分配二层决策模型进行 求解，得到最优解。针对河流流量这一不确定参数，采用模糊随机数进行描述，并采用两步 法进行处理，即首先将其转化成为模糊数，然后采用模糊期望值算子将模糊数转化为确定 值。针对二层决策模型的求解，根据在污染税政策下水污染负荷分配的二层决策结构特征 以及河流系统中污染物的传递特性，提出了混合嵌套基于多智能体的动态极值算法的改进 粒子群算法（PS0)的方法进行求解。从而得到合理的流域污染税费方案和污染物削减方 案。
[0109] 实施例
[0110] 本发明主要研宄单一河流系统中污染税政策下的水污染负荷分配设计。单一河流 系统由地方环保部门（EPA)以及一系列坐落于同一流域的排污控制单元（Dischargers)组 成。环保部门为了衡量各排污者污染水排放情况，在其下游安置水质监测点以监测河流水 质情况，并根据流域污染物对人类和生态系统健康的影响，以及流域环境生态功能分区情 况，在各监测点制定相应的水质标准。在满足水质标准条件下，每个排污者被允许将部分污 染物削减后的污水排放入河流中。在污染税政策下，环保部门为化解上下游排污者之间水 污染负荷分配的不公平以及进一步限制各排污者的污水排放，在各个监测点制定相关污染 物污染税费标准，如果某监测点监测相应污染物浓度大于该标准，环保部门将向相应的排 污控制单元征收相应的污染税收。每个排污控制单元需要平衡相应的污水处理成本和污染 税费比例，决策污水处理的削减比例以实现最小成本达到相应水质标准。
[0111] 一般情况下，即使地方环保局和各排污控制单元是独立的决策主体，地方环保局 作为领导者具有优先决策权，所以各排污控制单元是从属者。因此，基于污染税的流域污染 负荷分配问题可以抽象成为一个合作型二层多从属者规划问题，整个过程可以分为两个环 节：污染税制定环节和污染物削减环节：在污染税制度环节，地方环保部门（领导者）制定 T个监测点的污染税收标准以获得最小的所有排污控制单元的单位排污成本的偏差。在污 染物削减环节：各排污控制单元（从属者）依据环保部门制定的税收标准以及上游水质情 况依次决策自身污水削减比例以获得最小的单位排污成本，如图1所示。
[0112] 为解决上述基于污染税的流域污染负荷分配问题，提出如下合作型二层多从属者 决策建模技术，具体步骤如下：
[0113] 1.不确定变量的量化与处理。
[0114] 在水质量管理问题中需要考虑不确定问题已经得到广泛的认知，因为不确定存在 于各种不同的系统组成部分和各种预决策环境当中。因此，现实水质量管理决策中的内部 复杂以及随机性使其不能单单的用传统的确定性优化方法来考虑与解决。因此，如要进一 步对决策中的不确定变量进行量化和处理。
[0115] (1)量化过程：在污染负荷分配决策中，对上游水流量的量化十分重要，然而，上 游水量是一个典型的不确定性，其大小受到多种因素影响，如温度、降雨量、蒸发量等。因 此，本发明采用模糊随机变量对未来时间水流量进行量化，即可以采取相应的水资源配置 决策由于河流水流量的变化，未来的水供应是不确定的，水流量一般可以分为三个不同的 级别：低水位（L)，中水位（M)和高水位（H)，服从概率分布（a，b，c)，三种水位的概率分别 可以表示为&，pjPpH。然而，对于每一个水位，区域水资源管理局很难给定一个确切的值， 因此需要用模糊概念来描述它。为了搜集每一个水位的数据，针对水位k，找相应的水文专 家进行数据采集，他们通常用口头表述"水流量在ak和ck之间，并且有一个最有可能的值 bk。"其中，k表示水位（t=L表示低水位，t=M表示中水位，t=H表示高水位）。值得 一提的是，这些专家给定的口头描述是通过多次的观察的数据获得。因此在每一个水位上 的水流量可以用三角模糊数描述，综上，水流量（可供的分配水量）可以考虑成为一个模糊 随机变量。
[0117] (2)处理过程：一般地，由于模糊随机变量的存在我们很难获得模型的解。本发明 提出一种两步转化方法，将模糊随机变量转化成为确定型变量。第一步：将模糊随机变量转 化成为模糊数；第二步：运用带乐观-悲观指数的期望值算子将模糊数转化成为确定型变 量。关于上游水流量^ ,具体的转换步骤如下：
[0118] ①假设水流量可以表述成为任意一个模糊随机变量x，含有以下数值：
[0120] 其中，A，0、,,，0"表示三个水流量水平，因此其概率分布如下：
[0122] 其中，Q表示三种水流水平及其模糊逻辑；
[0123] 假设：0 (x)是模糊随机变量x的模糊参数，如果0 (x) =Ef[X|P]，e[0,1] 存在，则有：
[0126] 那么，模糊随机变量x的模糊期望值可以表示为f兩]，则有：
[0128]其中，aE=inf( 0 (x)) 0，bE=inf( 0 (x))i=sup( 0 (x)) "cE=sup( 0 (x)) 0， E'[§]表示模糊随机变量x的模糊期望值。
[0129] ②模糊随机水流量被转化成为三角模糊数，系统采用的带乐观-悲观指数的期望 值算子将其转化成为确定性值；
[0130] 首先，定义一个乐观-悲观指数（〇<e< 1)用于量化决 策者的决策偏好，然后 使用期望值算子计算该三角模糊数的期望值，并将其作为水流量的确定性变量，则有：
[0132] 其中，EIE1^]]表示水流量的确定性变量。
[0133] 2.河流同化能力模拟。
[0134] 河流同化能力表示一定的水域、水量条件下，基于维持水域一定的含氧水平，该水 域对有机污染物的自净量。其与河段的流量、流速、水温、允许的最低溶解氧含量及污染物 本身的特性有关。确定水体对污染物的同化能力是确定地区有机污染物排放标准以及制度 合理的污染负荷分配方案的基础。本发明使用Streeter-Phelps模型来模拟河流对有机污 染物的同化过程。其主要原理为：当河流受纳有机物后，沿水流方向产生的输移有机物远大 于扩散稀释量，当河水流量与污水流量稳定，河水温度不变时，有机物生化降解的好氧量与 该时期河水中存在的有机物量成正比。因此，有两个微分耦合模型被定义来描述河流中溶 解氧0〇)和生化需氧量（B0D)浓度变化规律。
[0135] (1)生化需氧量（B0D)转换机制：生化需氧量表示水中有机物等需氧污染物质含 量的一个综合指标。说明水中有机物由于微生物的生化作用进行氧化分解，使之无机化或 气体化时所消耗水中溶解氧的总数量。其值越高说明水中有机污染物质越多，污染也就越 严重。根据Streeter-Phelps模型的结论，河流中B0D浓度的随时间变化量与其自身浓度成 正比，即
，其中Cb表示B0D浓度，kd表示脱氧系数。因此，在一个河段t，用kd(t) 表示脱氧系数，d(t)表示该河段长度，v(t)表示该河段平均水流流速，Cb(t)表示流入该河 段水流的B0D浓度，基于以上微分模型的解析解，则流出该河段水流的B0D浓度Cb(t+1)为
，该浓度等价于：
[0136]Cb(t+1) =a(t)Cb(t) 〇
[0137] (2)溶解氧（DO)转换机制：水中溶解氧的多少是衡量水体自净能力的另一个 指标。根据Streeter-Phelps模型的结论，河流中的溶解氧从空气中补充，其补充浓度 的随时间变化量与其氧亏（即水体中饱和溶解氧和现存溶解氧的差）浓度成正比，即
1其中Cb表示B0D浓度，D。表示氧亏浓度，kd，匕表示分别脱氧和充氧系 数。与B0D情况相似，在一个河段1用ka(t)表示充氧系数，Cb(t)，Cjt)分别表示流入该 河段水流的B0D浓度，Cd。表示该河段水流的饱和溶解氧浓度，则流出该河段水流的D0浓度 C0(t+1)为：
[0139] 具体如图2所示。
[0140] 定义
和S⑴=CdJl-|3⑴]，该浓度等价于：
[0141]C0(t+1)=|3 (t)C0(t)-y(t)Cb(t)-8 (t)〇
[0142]其中，0⑴表示河段t的溶氧浓度的转换系数，Y(t)、S⑴分别表示影响河段 t溶氧浓度的生化需氧量浓度、饱溶氧浓度的转换系数。
[0143] 3.下层模型构建。
[0144] 对于下层模型，从河流上游到下游，每个排污控制单元为了达到最小的污染成本， 根据地方环保部门制定的污染税收标准以及上游的水质情况依次进行独立的污水B0D削 减决策。
[0145] (1)河流系统状态方程组构建。
[0146] 河流系统中的水以及水中的各类化学物质沿着河道顺流而下，在流经过各个的排 污控制单元过程中，控制单元取水、排污等活动会一定程度的影响河流状态（水量、各类污 染物浓度等）的改变，并且该状态会随水的流动转移到下一个控制单元，如图3所示。因 此，在污染负荷分配问题中，需关注河流状态在各个排污之间的转移情况，即流量状态转移 方程、B0D浓度状态转移方程、D0浓度状态转移方程被构建。
[0147] ①流量状态转移方程。基于质量守恒定律以及流量连续性条件，用Qjt-1)，Qjt) 分别表示在监测点t-1和t的流量，用Qw(t)，Qd(t)分别表示排污控制单元t的取水量和排 污量，则可构建流量状态转移方程如下：
[0148]Qr (t) =Qr (t-1) -Qw (t) +Qd (t) 〇
[0149] ②BOD浓度状态转移方程。基于质量守恒定律、Streeter-Phelps水质模拟模型 以及流量状态转移方程，用Cb(t_l)，Cb(t)表示监测点t-1和t的B0D浓度，q⑴表示排污 控制单元t未处理的污水B0D浓度，以及n(t)为排污控制单元t进行削减污水B0D浓度 的比例，其为下层排污者t的决策变量，则可构建B0D浓度状态转移方程如下：
[0151]定义变量
，则该状态转移方程等价于：
[0153] ③D0浓度状态转移方程。相似于B0D情形，基于质量守恒定律、Streeter-Phelps 水质模拟模型以及流量状态转移方程，用（；(t-l)，Cjt)表示监测点t-1和t的DO浓度，C.,⑴表示排污控制单元t未处理的污水D0浓度，则可构建D0浓度状态转移方程如下：
[0155] 同样的定义变量
，则该状态转移方程等价于：
[0157] 另外，考虑河流中饱和溶解氧浓度Cd。的限制，以上D0浓度状态转移方程可改写 为：
[0159] (2)初始模型河流状态条件构建。
[0160] 对于位于整个河流系统中上游第一个排污控制单元，用C:，C:分别表示上游水流 B0D和D0浓度，则该控制单元的模型考虑初始河流状态条件如下：
[0162] (3)下层模型约束条件构建。对于下层模型的决策者-排污控制单元，在进行相应 污水B0D浓度削减的决策过程中，需满足一下两个约束条件：允许的B0D和D0浓度标准以 及B0D削减能力约束。
[0163] ①允许的B0D和D0浓度标准。排污控制单元t排放污水到河流中，不能使得在相 应检测点河流的B0D浓度（Cb(t))超过环保部门所允许的最大B0D浓度标准（Bu(t))，即：
[0164] 0 ^Cb(t) ^Bu(t) 〇
[0165] 同样，不能使相应检测点河流的DO浓度（C；(t))低于环保部门所允许的最小DO浓 度标准汍⑴），即：
[0166] Dx(t) ^C0(t)〇
[0167] ②BOD削减能力约束。每个排污控制单元具有污水BOD浓度削减能力的限制，通 常采用二级污水处理工艺，如图4所示。对于排污控制单元t，第一级处理后的污水B0D浓 度削减率为，第二级处理后的污水bod浓度削减率为nr。在实际的污水处理运行中，环 保部门要求所以排污者的污水必须进行第一级处理。因此，排污控制单元t进行污水B0D 浓度削减决策，其B0D削减率（n(t))不能够超过nr，同时不能低于<in，BP:
[0169] (4)下层模型目标函数构建。下层决策者的目标为最小的单位污染成本目标函数， 污染成本由污水处理成本和污染税费组成。
[0170] ①污水处理成本。基于B0D削减和污水处理成本的经验曲线，如图4所示，本发明 对于排污控制单元的污水处理费用采用幂函数形式，用kel(t)和ke2(t)表示排污控制单元 t的成本参数1和2,Mt)和keb⑴表示排污控制单元t的污水和B0D削减参数，则污水 处理成本C(t)可以表示为以下方程：
[0172] 其次，参数1^(〇，1^2(〇，1^(〇和1^(〇数值以污水处理构筑物建设和运营所得 数据，采用最小二乘法进行曲线拟合所得。
[0173] ②污染税费。在每个监测点，地方环保部门制定相应的污水税收标准（Alb(t)，上 层模型的决策变量）和税率（Tb)。如果检测点t所对应的河流B0D浓度（Cb(t))低于该税 收标准（Cb(t)彡Alb(t))，则排污控制单元t不需要征收污染税（P(t));如果检测点t所对 应的河流B0D浓度（Cb(t))超过该税收标准（(；(〇>八113(〇)，则排污者七需征收污染税？(〇 =(t)Tb [Cb (t)-Alb⑴]。因此，排放者t的污染税费模型定义为：
[0175] 综上所述，排污控制单元t的目标函数，即最小化的单位污染成本如下：
[0177] 4?上层模型（地方环保部门）构建。
[0178] 对于上层模型，地方环保部门为实现公平的流域污染负荷分配，于河流系统中每 个监测点制度相应的污染税收标准。
[0179] (1)上层模型约束条件构建。对于上层模型的决策者-地方环保部门，在进行相应 污染税收标准的制定过程中，需满足允许的BOD浓度标准约束条件。
[0180] 环保部门所制定的污染税收标准（Alb(t))不能超过允许的最大B0D浓度标准 (Bu(t))，即：
[0181] 0 ^Alb(t) ^Bu(t) 〇
[0182] (2)上层模型目标函数构建。上层决策者需保证公平的污染负荷分配，因此，本 发明采用各排污控制单元的单位污染成本的方差来量化水污染负荷分配的公平性。用
表示各排污控制单元平均的单位污染成本，则其最小的方差定义为：
[0184] 5.全局模型
[0185] 集成上下层目标函数和上下层约束条件，我们可以得到以下基于污染税的流域污 染负荷分配问题的合作型二层多从属者规划全局模型。
[0187]
[0188] 如图5所示，表示了针对基于污染税的流域污染负荷分配问题的合作型二层多从 属者决策优化建模技术流程图。
[0189] 针对上述合作型二层多从属者决策模型，本发明提出了混合嵌套基于多智能体的 动态极值算法的改进粒子群算法进行求解。由于二层规划问题具有较为特殊的决策规则， 即上层决策者只是通过自己的决策去指导下层决策者，并不直接干涉下层的决策；而下层 决策者只需要把上层的决策作为参数，他可以在自己的可能范围内自由决策。这种决策机 制使得上层决策者在选择策略以优化自己的目标达成时，必须考虑到下层决策者可能采取 的策略对自己的不利影响。其求解问题已经为证明为NP-难题，即使是最简单的情况。然 而，上述模型中存在非线性的目标函数和约束条件，增加模型计算的难度；另一方面，该模 型相比于传统的二层模型具有更为特殊的动态决策结构，即下层模型需要通过河流系统状 态转移方程进行依次计算，显然，使得该模型的求解问题变的更加困难。因此，本发明针对 基于污染税的流域污染负荷分配二层决策模型，将传统精确算法与现代启发式智能算法相 结合，设计了一种混合 嵌套粒子群算法，具体包括：①针对上述下层模型特殊的多元动态决 策结构，设计一种基于多智能体的动态极值算法为精确快速的获得基于特定上层决策变量 的下层模型的最优解；②为计算上层模型，将该算法嵌套至粒子群算法的各个粒子适应值 计算的过程中；③为提高整体算法的全局搜索能力，采用Boltzmann选择算子代替传统全 局最优粒子的更新机制。下面介绍其主要设计思路。
[0190] 1.粒子结构设计。为简化描述，以决策排污者为T为例介绍算法的粒子结构设计。 在混合嵌套粒子群优化算法中，每个粒子的位置编码为T维的向量，以表示解空间中的一 点。每个粒子位子代表上层模型一个决策变量（Alb (t)，t= 1，2,. . .，T)，如下
[0192] 其中，t表示迭代次数指标，t= 1，2, ...，I，I表示最大迭代次数，1表示粒子指 标，1 = 1，2, ...，L，L表示种群指标。图8表示了一个包括T个排污控制单元的河流系统 的粒子结构表达。
[0193] 2.粒子性能的混合嵌套计算。基于二层规划特殊的决策规则，在使用粒子群算法 进行上层模型的计算的过程中，对各代粒子性能的计算时，需要获得基于特定上层决策变 量的下层模型的最优解。因此，针对基于污染税的流域污染负荷分配问题的合作型二层多 从属者决策模型中下层模型的多元动态特征，本发明为求得其最优反应函数，设计了一种 基于多智能体的动态极值算法。
[0194] (1)基于多智能体的动态极值算法。为求解在特定上层决策变量下的上述模型的 下层模型，需求解T个一元的非线性规划模型。因此，针对单一的非线性规划模型，本发明 提出了一种两阶段极值算法，即第一阶段通过极值算法求解该问题极值，第二阶段通过约 束处理方法从第一次获得极值中找到该问题的最优解。在解决单一下层模型求解问题后， 本发明基于下层各个模型间的动态传递关系，构建了一种多智能体动态框架，以获得下层T 个模型的全部最优解。
[0195] ①极值算法。以下层第t个单一模型为例，介绍极值算法求解下层单一模型 的极值。令下层模型目标函数中的Cb(t) =Alb(t)，我们能得到该目标函数的断点为
因此，该模型的极值存现两种情况：
[0196] 情况I当Cb(t)彡Alb(t)时，单一下层模型的目标函数如下：
[0198]
该目标函数等价于
[0199] 巧=atri(t)p'+ct
[0200] 因此，在该种情况下，目标函数的极值点为nel(t) =np(t)。
[0201] 情况II当Cb(t)>Alb(t)时，单一下层模型的目标函数如下：
[0202]
[0203] 令1^ =TBCUt)，该目标函数等价于：
[0205]因此，在该种情况下，该目标函数的极值点为
>考虑极值点的范 围，对该种情况下目标函数的极值点做以下修正：当ne2(t)彡np(t)时，下层模型目标 函数的极值点为nel(t) =np(t);当o〈ru2(t)〈np(t)，下层模型目标函数的极值点为
。该极值算法的具体步骤如表1所示：
[0206] 表1:极值算法的步骤：
[0208]
[0209] ②约束处理方法。在计算下层单一模型目标函数的极值后，需对下层目标进行约 束处理，通过对下层模型相应的约束条件的等价转换，允许的B0D浓度标准约束变为：
[0211]
该约束可以等价于：
[0213] 允许的BOD浓度标准约束变为：
[0215] 令gt等于上述不等式的右边，该不等式二次转换为：
[0217]
该约束可以等价于：
[0219] 因此，下层单一模型约束处理方法的具体步骤如表2所示。
[0220] 表2,下层单一模型约束处理方法的具体步骤：
[0223] ③多智能体动态框架。基于下层模型之间河流系统动态转移方程，本发明设计了 一种基于多智能体的动态框架，结合上述提出的极值算法，以求得在特定上层决策下的全 部下层模型的最优反应集。其具体步骤如表3所示。
[0224] 表3基于多智能体的动态极值算法的求解过程；
[0225]
[0227] (2)适应值函数。通过上述所提出的算法能够计算获得在特定上层决策变量下 的下层模型的最优解，根据双层模型的决策规则，可以得到量化例子性能的适应值函数，即 为：
[0229] 其中，
。因此，整个粒子性能的混合嵌套计算过程可以描述为图7所 示的计算方法。
[0230] 3.粒子更新机制。基于粒子群算法的运行程序，在量化粒子性能后需进行粒子更 新，为获得更好性能的粒子。粒子更新方程如下：
[0232] 卩丄（t+1) = % (t+1)+Pi(t);
[0233]其中
表不当前个体最优粒子，即Fitness^13)SFitness^。)），对o =1，2,…，T;Pe=[PlQ(l)，p「(2)，...，pf(T)]表示当前粒子群中最优粒子，即FitneSS(P°) 莳腦8(杧），对 1 = 1，2，...，L; 0 (t)表示惯性权重，计算过程为
，其中0 (1)和 0 (I)表示第1和第I代的惯性权重；cp和c 别表示个体和全局最优学习因子；rp和r1 表示0到1之间的随机数。为提供粒子的全局搜索能力，本发明采用Boltzmann选择算子 代替传统全局最优粒子的更新机制。
[0234] (1)Boltzmann选择算子。采用Boltzmann选择算子更新全局最优粒子，其基本思 路：每代全局最优粒子1><3=[1 )「(1)，?「(2)，...,1)「(1')]，通过运用13〇1〖21]^1111选择算子在各粒子的个 体最优粒子矿=[pf(l)，pf(2)，...，pf(T)]中进行选择。为计算各个粒子的个体最优粒子的选择概 率，其评价函数定义如下：
[0236] 其各粒子的个体最优粒子的选择概率为
[0238] 其中TejPTe:表示初代和第t代温度。因此，运用Boltzmann选择算子进行各 粒子的全局最优粒子的具体步骤如表4及图8所示。
[0239] 表4运用Boltzmann选择算子进行各粒子的全局最优粒子的具体步骤；
[0240]
[0242] 4.混合嵌套粒子群算法。该算法的整体步骤如表5所示：
[0243] 表5混合嵌套粒子群算法步骤；
[0245]
[0246] 为了能进一步了解本发明的
【发明内容】
、特点及用法，兹举以如下实例详细说明。
[0247] 实例：四川省沱江流域某河段为例，其具有5个排污控制单元= 〇, 1，2, 3, 4, 5)，整个区域的考虑基于污染税政策下的污染负荷分配决策，该河流系统见图 9。通过对该流域水位数据访问调研，该流域饱和溶解氧浓度为10. 08mg/L，其他的水位数据 见表6。对于各排污控制单元，其污水处理成本参数（kel⑴，ke2 (t)，Mt)和(t))数值 以污水处理构筑物建设和运营所得数据，采用最小二乘法进行曲线拟合所得，其他相关数 据见表7。对于河流系统的上游情况，其B0D和D0浓度为3. 2mg/L和6. 2mg/L，其水流量为 模糊随机变量，其分布数据见8。假设上层决策的偏好乐观-悲观指数设置为e=0.5,通 过计算上游水流量的期望值为12. 10m3/s。该河流系统各监测点的环境水质标准见表9,另 外基本污染税率为〇? 005yuan?L/mg。
[0248] 表6河段水位数据；
[0250]
[0251] 表7排污控单元数据；
[0253] 表8河流水流量分别数据；
[0254]
[0255]表9河流系统各监测点的环境水质标准
[0257] 计算结果。应用本发明提出的基于污染税的问题的多阶段建模技术对该实例问题 进行建模，并采用本发明提出的两步转化方法对模型中的三角模糊随机数进行处理，最后 应用本发明提出的混合嵌套粒子群算法对模型进行求解计算，结果如表10所示。通过将算 法优化的结果反映，各排污控制单元的污染成本都为1. 72yuan?s/m3。因此，地方环保部 门通过污染税收的制定，实现了公平的污染负荷分配。这种分配结果可以能够实现流域健 康可持续的发展。尽管算法优化的结果是基于模型假设的条件得出的，但在一定程度上仍 能反映出实际决策方案尚有提升的空间。在实际决策过程中，各决策者仅仅根据经验判断 决定决策方案具有一定风险。运用模型和算法计算的结果可对决策者进行方案选择提供参 考，以更好地指导水质管理实践。
[0258] 表10模型求解结果
[0259]
【主权项】
1. 基于污染税的流域污染负荷分配二层决策优化方法，其特征在于，包括以下步骤： 步骤1、系统获取水流量，并将其转换为水流量的确定性变量； 步骤2、系统获取河流中溶解氧浓度关系及生化需氧量浓度关系； 步骤3、系统根据水流量的确定性变量、溶解氧浓度关系及生化需氧量浓度关系构建下 层模型； 步骤4、系统根据污染税收标准制定上层模型； 步骤5、系统根据下层模型及上层模型构建基于污染税的流域污染负荷分配二层决策 丰旲型； 步骤6、系统对基于污染税的流域污染负荷分配二层决策模型进行求解，得到最优解。2. 根据权利要求1所述的基于污染税的流域污染负荷分配二层决策优化方法，其特征 在于，所述步骤1中，首先，系统采用模糊随机变量对未来时间水流量进行量化，系统将水 流量分为三个不同的级别：低水位（L)，中水位（M)和高水位（H)，服从概率分布（a，b，c)， 三种水位的概率分别可以表示为Py Pm和P H，综上，水流量$为：其次，将水流量^转换为模糊随机变量，并计算其模糊期望值，具体包括： ① 假设水流量可以表述成为任意一个模糊随机变量X，含有以下数值： Ol _ (aL'bL'CL)，Qm _ (aM'bjvI，CM)，Qh _ (aH，bH'CH). 其中，I，，O11表示三个水流量水平，因此其概率分布如下： Ρ({ω e Ω I χ(ω) = QL}) = pL, Ρ({ω e Ω | χ(ω) = Qm }) = ρΜ, Ρ({ω ε Ω | χ(ω) = Qh }) = ρΗ 9 其中，Ω表示三种水流水平及 其模糊逻辑； 假设：Θ (X)是模糊随机变量X的模糊参数，如果Θ (X) =Ef [X I P]，Vae [0，1]存在， 则有：那么，模糊随机变量X的模糊期望值可以表示为1^造1，则有： E'[Q] = (aE,bE,cE) _ 其中，aE= inf ( θ ( x )) Q, bE= inf ( θ ( x )) i = sup ( θ ( x )) Ce= sup ( θ ( x )) 〇， 表示模糊随机变量x的模糊期望值； ② 模糊随机水流量被转化成为三角模糊数，系统采用的带乐观-悲观指数的期望值算 子将其转化成为确定性值； 首先，定义一个乐观-悲观指数（〇< ε <1)用于量化决策者的决策偏好，然后使用 期望值算子计算该三角模糊数的期望值，并将其作为水流量的确定性变量，则有：其中，Ε^Ε?^]]表示水流量的确定性变量。3. 根据权利要求1所述的基于污染税的流域污染负荷分配二层决策优化方法，其特征 在于，所述步骤2中，系统根据Streeter-Phelps模型的结论，河流中生化需氧量浓度的随 时间变化暈与其自身浓度成正比，即：其中，Cb表示生化需氧量浓度，k 示脱氧系数； 因此，在一个河段表示脱氧系数，d(t)表示该河段长度，v(t)表示该河段 评价水流流速，Cb(t)表示流入该河段水流的生化需氧量浓度，则流出该河段水流的生化需 氧量浓度Cb (t+Ι)为，该生化需氧量浓度等价于： Cb(t+1) = a (t)Cb(t)； 系统根据Streeter-Phelps模型的结论，河流中的溶解氧从空气中补充，其补充浓度 的随时间变化量与其氧亏浓度成正比，即：其中，Cb表示生化需氧量浓度，D。表示氧亏浓度，1^表示脱氧系数，1表示充氧系数。 与生化需氧量情况相似，在一个河段1用ka(t)表示充氧系数，Cb (t)，CJt)分别表示流入 该河段水流的生化需氧量浓度，Cd。表示该河段水流的饱和溶解氧浓度，则流出该河段水流 的溶解氧浓度cjt+l)为：定义：-和δ (t) = CdJl-β⑴]，该溶解氧浓度等 价于： C0(t+1) = β (t)C0(t)-y (t)Cb(t)-5 (t)； 其中，a (t)、β (t)分别可以表示河段t的生化需氧量浓度、溶氧浓度的转换系数， γ (t)、δ (t)分别表示影响河段t溶氧浓度的生化需氧量浓度、饱溶氧浓度转换系数。4. 根据权利要求1所述的基于污染税的流域污染负荷分配二层决策优化方法，其特征 在于，所述步骤3中，系统构下层模型包括： (1)河流系统状态方程组构建，具体包括： ① 流量状态转移方程：基于质量守恒定律以及流量连续性条件，用A(t-l)，Qjt)分别 表示在监测点t-Ι和t的流量，用Qw (t)，Qd (t)分别表示排污控制单元t的取水量和排污 量，则可构建流量状态转移方程如下： Qr (t) = Qr(t-1)-Qw(t)+Qd(t)； ② BOD浓度状态转移方程：基于质量守恒定律、Streeter-Phelps水质模拟模型以及流 量状态转移方程，用Cb (t-1)，Cb (t)表示监测点t-Ι和t的BOD浓度，（'Ut)表示排污控制单 元t未处理的污水BOD浓度，以及n (t)为排污控制单元t进行削减污水BOD浓度的比例， 其为下层排污者t的决策变量，则可构建BOD浓度状态转移方程如下： Cb (t)Qr(t) = a(t- l)Cb (t - l)[Qr (t -1) - Qw (t)]+(t)[l - n(t)]Qd(t); 定义变I?则该状态转移方程等价于： Cb (t) = LI - r(t)ja(t - l)Cb (t -1) + r(t)C；1 (t)[l - η(?)]; ③DO浓度状态转移方程：相似于BOD情形，基于质量守恒定律、Streeter-Phelps水质 模拟模型以及流量状态转移方程，用C。(t-1)，C。⑴表示监测点t-Ι和t的DO浓度，C〖(t)表 示排污控制单元t未处理的污水DO浓度，则可构建DO浓度状态转移方程如下： C0(t)Qr(〇 = LP(t-l)Q(t-l)-T(t-l)Cb(l-l)-6(l-l)JLQI(t-l)-Q w(i)J+C· (I)Qd(I)： 同样的定义变量I则该状态转移方程等价于：另外，考虑河流中饱和溶解氧浓度Cd。的限制，以上DO浓度状态转移方程可改写为： C^t) = inin{[l-r(t)][p(l-l)C；At-l)-y(t-l)Cb(t-l)-6(t-l)]+r(t)CUi).C dJ ； (2) 构建初始模型河流状态条件： 对于位于整个河流系统中上游第一个排污控制单元，用C；；，C:分别表示上游水流BOD 和DO浓度，则该控制单元的模型考虑初始河流状态条件如下： Qr(O) = Εε[Ε'[4]], Cb(O) = Q, C0(O) = C^ ； (3) 构建下层模型约束条件： ① 允许的BOD和DO浓度标准： 排污控制单元t排放污水到河流中，不能使得在相应检测点河流的BOD浓度（Cb(t))超 过环保部门所允许的最大BOD浓度标准（Bu (t))，即： 0. Cb (t) ^ Bu (t)； 同样，不能使相应检测点河流的DO浓度（C；(t))低于环保部门所允许的最小DO浓度 标准（D1⑴），即： D1 (t) ^ C0 (t)； ② BOD削减能力约束： 每个排污控制单元具有污水BOD浓度削减能力的限制，通常采用二级污水处理工艺。 对于排污控制单元t，第一级处理后的污水BOD浓度削减率为η;111π，第二级处理后的污水BOD 浓度削减率为ηΓ，在实际的污水处理运行中，环保部门要求所以排污者的污水必须进行第 一级处理；因此，排污控制单元t进行污水BOD浓度削减决策，其BOD削减率（τι (t))不能 够超过ηΓ，同时不能低于ηΓ，即： OWWt; (4) 构建下层模型目标函数： 下层决策者的目标为最小的单位污染成本目标函数，污染成本由污水处理成本和污染 税费组成； ① 污水处理成本： 基于BOD削减和污水处理成本的经验曲线，本发明对于排污控制单元的污水处理费用 采用幂函数形式，用L⑴和U⑴表示排污控制单元t的成本参数1和成本参数2, Mt) 和Ι^α)表示排污控制单元t的污水和BOD削减参数，则污水处理成本C(t)可以表示为以 下方程：其次，参数ha)，u(t)，Mt)和u⑴数值以污水处理构筑物建设和运营所得数 据，采用最小二乘法进行曲线拟合所得； ② 污染税费： 在每个监测点，地方环保部门制定相应的污水税收标准（Alb(t)，上层模型的决 策变量）和税率（Tb)，如果检测点t所对应的河流BOD浓度（Cb(t))低于该税收标准 (Cb(t) <Alb(t))，贝IJ排污控制单元t不需要征收污染税（P(t));如果检测点t所对应的 河流如0浓度（(；(〇)超过该税收标准（(；(〇>4113(〇)，则排污者〖需征收污染税？(〇 = Qi (t) Tb [Cb (t)-Alb⑴]，因此，排放者t的污染税费模型定义为：综上所述，排污控制单元t的目标函数，即最小化的单位污染成本如下：5. 根据权利要求1所述的基于污染税的流域污染负荷分配二层决策优化方法，其特征 在于，所述步骤4中，系统构上层模型包括： (1) 上层模型约束条件构建： 对于上层模型的决策者-地方环保部门，在进行相应污染税收标准的制定过程中，需 满足允许的BOD浓度标准约束条件，环保部门所制定的污染税收标准（Alb(t))不能超过允 许的最大BOD浓度标准（Bu⑴），即： 0 ^Alb (t) ^ Bu (t)； (2) 上层模型目标函数构建： 上层决策者需保证公平的污染负荷分配，因此，本发明采用各排污控制单元的单位污 染成本的方差来量化水污染负荷分配的公平性，用表示各排污控制单元平均的 单位污染成本，则其最小的方差定义为：6. 根据权利要求1所述的基于污染税的流域污染负荷分配二层决策优化方法，其特征 在于，所述步骤5中，系统根据下层模型及上层模型构建的基于污染税的流域污染负荷分 配二层决策模型为：7. 根据权利要求1所述的基于污染税的流域污染负荷分配二层决策优化方法，其特征 在于，所述步骤6中，假设决策排污者为T个，则在特定上层决策变量下地下层模型为T个 一元的非线性规划模型；首先，系统利用极值算法计算出单一的非线性规划模型的极值，再 通过下层约束条件计算出该单一的非线性规划模型的最优解，进而获得下层T个一元的非 线性规划模型的全部最优解集；利用粒子群算法对上述最优解集进行计算，从而获得基于 污染税的流域污染负荷分配二层决策模型的最优解。8. 根据权利要求7所述的基于污染税的流域污染负荷分配二层决策优化方法，其特征 在于，粒子群算法的粒子更新方程如下： Ρ^τ+Ι) = V^t+D+P^t)；其中，杧=[?!5(1)，1^(2)，...￥(丁)]表示当前个体最优粒子，即府11 (^叱)<]%^(?1(〇))，对0 = 1，2，···，τ ; Pc=丨pf (1),4(2),..., pf (T)]表示当前粒子群中最优粒子，SPFitness(Pc)i; Fitness(P1B),对 I =1，2，...，L; θ (τ)表示惯性权重，计算过程为，其中Θ⑴和 Θ⑴表不第1和第I代的惯性权重；Cp和％分别表不个体和全局最优学习因子；rp和％表 示O到1之间的随机数。9.根据权利要求7所述的基于污染税的流域污染负荷分配二层决策优化方法，其特征 在于，系统采用Boltzmann选择算子代替传统全局最优粒子的更新机制。
【专利摘要】本发明提出了基于污染税的流域污染负荷分配二层决策优化方法，首先，系统获取水流量，并将其转换为水流量的确定性变量；其次，系统获取河流中溶解氧浓度关系及生化需氧量浓度关系；再其次，系统根据水流量的确定性变量、溶解氧浓度关系及生化需氧量浓度关系构建下层模型；然后，系统根据污染税收标准制定上层模型；再然后，系统根据下层模型及上层模型构建基于污染税的流域污染负荷分配二层决策模型；最后，系统对基于污染税的流域污染负荷分配二层决策模型进行求解，得到最优解。适用于基于污染税的流域污染负荷分配问题。
【IPC分类】G06Q10/04
【公开号】CN104899668
【申请号】CN201510390525
【发明人】徐玖平, 张梦翔, 曾自强
【申请人】四川大学
【公开日】2015年9月9日
【申请日】2015年7月3日
【技术领域】
[0001] 本发明涉及流域水质目标管理，特别涉及基于污染税的流域污染负荷分配决策方 法。
【背景技术】
[0002] 二层决策优化方法是一类针对具有二层递阶结构问题的优化方法。其主要研宄具 有两个层次系统的规划与管理问题。分为上层和下层，上下两层问题都具有各自的决策变 量、约束条件和目标函数。这种方法区别于其他优化方法的一个主要方面是不再是一个决 策者进行规划与决策，而是有上下两层决策者，上层决策者只是通过自己的决策去指导下 层决策者，并不直接干涉下层的决策；而下层的决策将上层的决策作为一个影响条件，其在 自己的可行域进行自由决策，这就导致上层决策者在做任何决策的时候必须考虑到下层决 策者可能做出的反馈，以规避下层决策者的决策可能带来的不利影响。
[0003] 在水污染严重的今天，合理的流域水质管理显得至关重要。目前，为进一步治理各 流域水环境，我国环保部制定了一系列的环境经济政策，包括环境税费、排污权交易等。
[0004] 其中，污染税费为主要方式。其旨在国家以立法的形式明确规定污染者要按其排 放的污染物种类、浓度和数量缴纳一定的费用。通过征收污染税，使经济活动所造成的环境 损失由其本身承担，使外部经济内在化；此外，税金的收入可用作改善环境项目的拨款。污 染税费制定在限制各流域排放者的排污行为的同时，也一定程度的增加了地方环保部门的 水污染治理工作难度。具体表现为：①如何设计一套科学的收费体系，可以使某一流域和地 区各排污者实现特定水质标准所需的成本相同；②如何制定一套合理的流域污染负荷控制 技术，能够指导地方环保部门在污染税费政策下进行水污染防治工作。因此，本发明考虑了 基于污染税的流域污染负荷分配问题，旨在提出一种新的水环境污染分配方案。
[0005] 在过去的几十年中，决策优化技术在水资源管理方面起到了至关重要的作用，主 要相关的技术包括多目标决策优化，多属性决策优化以及多阶段决策优化等，这些决策的 一个主要特点是由一个主体进行决策，没有考虑其他利益相关者的意见及反馈。
【发明内容】
[0006] 本发明所要解决的技术问题，就是提供一种基于污染税的流域污染负荷分配二层 决策优化方法以解决基于污染税的流域污染负荷分配问题。
[0007] 本发明解决所述技术问题，采用的技术方案是，基于污染税的流域污染负荷分配 二层决策优化方法，包括以下步骤：
[0008] 步骤1、系统获取水流量，并将其转换为水流量的确定性变量；
[0009] 步骤2、系统获取河流中溶解氧浓度关系及生化需氧量浓度关系；
[0010] 步骤3、系统根据水流量的确定性变量、溶解氧浓度关系及生化需氧量浓度关系构 建下层模型；
[0011] 步骤4、系统根据污染税收标准制定上层模型；
[0012] 步骤5、系统根据下层模型及上层模型构建基于污染税的流域污染负荷分配二层 决策模型；
[0013] 步骤6、系统对基于污染税的流域污染负荷分配二层决策模型进行求解，得到最优 解。
[0014] 具体的，所述步骤1中，首先，系统采用模糊随机变量对未来时间水流量进行量 化，系统将水流量分为三个不同的级别：低水位（L)，中水位（M)和高水位（H)，服从概率分 布（a，b，c)，三种水位的概率分别可以表示为九，pjPpH，综上，水流量匕为：
[0016] 其次，将水流量匕转换为模糊随机变量，并计算其模糊期望值，具体包括：
[0017] ①假设水流量可以表述成为任意一个模糊随机变量X，含有以下数值：
[0019] 其中，4，4,，表示三个水流量水平，因此其概率分布如下：
[0021] 其中，D表示三种水流水平及其模糊逻辑；
[0022] 假设：0(x)是模糊随机变量x的模糊参数，如果0(x)= Ef [X | P]，V? e [0,1J 存在，则有：
[0025] 那么，模糊随机变量x的模糊期望值可以表示为E'虎]，则有：
[0027] 其中，aE=inf( 0 (x)) 0，bE=inf( 0 (x))i=sup( 0 (x)) "cE=sup( 0 (x)) 0， ElGl表示模糊随机变量x的模糊期望值；
[0028] ②模糊随机水流量被转化成为三角模糊数，系统采用的带乐观-悲观指数的期望 值算子将其转化成为确定性值；
[0029] 首先，定义一个乐观-悲观指数（0<e< 1)用于量化决策者的决策偏好，然后 使用期望值算子计算该三角模糊数的期望值，并将其作为水流量的确定性变量，则有：
[0031] 其中，ElEli]]表示水流量的确定性变量。
[0032] 具体的，所述步骤2中，系统根据Streeter-Phelps模型的结论，河流中生化需氧 量浓度的随时间变化量与其自身浓度成正比，即：
[0034] 其中，Cb表示生化需氧量浓度，kd表示脱氧系数；
[0035] 因此，在一个河段t，用kd(t)表示脱氧系数，d(t)表示该河段长度，v(t)表示该河 段评价水流流速，Cb(t)表示流入该河段水流的生化需氧量浓度，则流出该河段水流的生化 需氧量浓度Cb (t+1)为
该生化需氧量浓度等价于：
[0036]Cb(t+1) =a(t)Cb(t)；
[0037] 系统根据Streeter-Phelps模型的结论，河流中的溶解氧从空气中补充，其补充 浓度的随时间变化量与其氧亏浓度成正比，即：
[0039] 其中，Cb表示生化需氧量浓度，D。表示氧亏浓度，kd表示脱氧系数，ka表示充氧系 数。与生化需氧量情况相似，在一个河段1用ka(t)表示充氧系数，Cb(t)，Cjt)分别表示 流入该河段水流的生化需氧量浓度，Cd。表示该河段水流的饱和溶解氧浓度，则流出该河段 水流的溶解氧浓度C。(t+1)为：
[0041] 定义：
和S⑴=(^。[1-|3(1：)]，该溶解氧浓 度等价于：
[0042]C0(t+1) = |3 (t)C0(t)-y(t)Cb(t)-8 (t)；
[0043]其中，a(t)、Mt)分别可以表示河段t生化需氧量浓度（BOD)、溶氧浓度（DO)的 转换系数，Y(t)、S(t)分别表示影响河段t溶氧浓度的生化需氧量浓度、饱溶氧浓度转换 系数。
[0044] 具体的，所述步骤3中，系统构下层模型包括：
[0045] (1)河流系统状态方程组构建，具体包括：
[0046]①流量状态转移方程：基于质量守恒定律以及流量连续性条件，用Qjt-l)，Qjt) 分别表示在监测点t_l和t的流量，用Qw(t)，Qd(t)分别表示排污控制单元t的取水量和排 污量，则可构建流量状态转移方程如下：
[0047]Qr(t) =Qr(t-1)-Qw(t)+Qd(t)；
[0048] ②BOD浓度状态转移方程：基于质量守恒定律、Streeter-Phelps水质模拟模型以 及流量状态转移方程，用Cb(t-1)，Cb(t)表示监测点t-1和t的B0D浓度，表示排污控 制单元t未处理的污水B0D浓度，以及n(t)为排污控制单元t进行削减污水B0D浓度的 比例，其为下层排污者t的决策变量，则可构建B0D浓度状态转移方程如下：
[0050] 定义变量
，则该状态转移方程等价于：
[0052] ③DO浓度状态转移方程：相似于BOD情形，基于质量守恒定律、Streeter-Phelps 水质模拟模型以及流量状态转移方程，用（；(t-l)，Cjt)表示监测点t-1和t的DO浓度， 表示排污控制单元t未处理的污水D0浓度，则可构建D0浓度状态转移方程如下：
[0054] 同样的定义变量
，则该状态转移方程等价于：
[0056] 另外，考虑河流中饱和溶解氧浓度Cd。的限制，以上D0浓度状态转移方程可改写 为：
[0058] (2)构建初始模型河流状态条件：
[0059] 对于位于整个河流系统中上游第一个排污控制单元，用q，C:分别表示上游水流 B0D和D0浓度，则该控制单元的模型考虑初始河流状态条件如下：
[0061] (3)构建下层模型约束条件：
[0062] ①允许的B0D和D0浓度标准：
[0063] 排污控制单元t排放污水到河流中，不能使得在相应检测点河流的B0D浓度 (Cb(t))超过环保部门所允许的最大B0D浓度标准（Bu(t))，即：
[0064] 0 ^Cb(t) ^Bu(t)；
[0065] 同样，不能使相应检测点河流的DO浓度（C；(t))低于环保部门所允许的最小DO浓 度标准汍⑴），即：
[0066] Dx(t) ^C0(t)；
[0067] ②B0D削减能力约束：
[0068] 每个排污控制单元具有污水B0D浓度削减能力的限制，通常采用二级污水处理工 艺。对于排污控制单元t，第一级处理后的污水B0D浓度削减率为，第二级处理后的污 水B0D浓度削减率为^厂，在实际的污水处理运行中，环保部门要求所以排污者的污水必须 进行第一级处理；因此，排污控制单元t进行污水B0D浓度削减决策，其B0D削减率（n(t)) 不能够超过nr，同时不能低于f，即：
[0070] (4)构建下层模型目标函数：
[0071] 下层决策者的目标为最小的单位污染成本目标函数，污染成本由污水处理成本和 污染税费组成；
[0072] ①污水处理成本：
[0073] 基于B0D削减和污水处理成本的经验曲线，本发明对于排污控制单元的污水处理 费用采用幂函数形式，用kel (t)和ke2 (t)表示排污控制单元t的成本参数1和成本参数2，kra(t)和表示排污控制单元t的污水和BOD削减参数，则污水处理成本C(t)可以表 示为以下方程：
[0075] 其次，参数kel (t)，ke2 (t)，kee⑴和keb⑴数值以污水处理构筑物建设和运营所得 数据，采用最小二乘法进行曲线拟合所得；
[0076] ②污染税费：
[0077] 在每个监测点，地方环保部门制定相应的污水税收标准（Alb(t)，上层模型的 决策变量）和税率（Tb)，如果检 测点t所对应的河流B0D浓度（Cb(t))低于该税收标准 (Cb(t) <Alb(t))，贝1J排污控制单元t不需要征收污染税（P(t));如果检测点t所对应的 河流如0浓度（(；(〇)超过该税收标准（(；(〇>4113(〇)，则排污者〖需征收污染税？(〇 = Qjt)Tb [Cb (t)-Alb (t)]。因此，排放者t的污染税费模型定义为：
[0079] 综上所述，排污控制单元t的目标函数，即最小化的单位污染成本如下：
[0081] 具体的，所述步骤4中，系统构上层模型包括：
[0082] (1)上层模型约束条件构建：
[0083] 对于上层模型的决策者-地方环保部门，在进行相应污染税收标准的制定过程 中，需满足允许的B0D浓度标准约束条件，环保部门所制定的污染税收标准（Alb(t))不能 超过允许的最大B0D浓度标准（Bu(t))，即：
[0084] 0 ^Alb(t) ^Bu(t)；
[0085] (2)上层模型目标函数构建：
[0086] 上层决策者需保证公平的污染负荷分配，因此，本发明采用各排污控制单元的单 位污染成本的方差来量化水污染负荷分配的公平性，用
表示各排污控制单元平 均的单位污染成本，则其最小的方差定义为：
[0088] 具体的，所述步骤5中，系统根据下层模型及上层模型构建的基于污染税的流域 污染负荷分配二层决策模型为：
[0090]
[0091] 具体的，所述步骤6中，假设决策排污者为T个，则在特定上层决策变量下地下层 模型为T个一元的非线性规划模型；首先，系统利用极值算法计算出单一的非线性规划模 型的极值，再通过下层约束条件计算出该单一的非线性规划模型的最优解，进而获得下层T 个一元的非线性规划模型的全部最优解集；利用粒子群算法对上述最优解集进行计算，从 而获得基于污染税的流域污染负荷分配二层决策模型的最优解。
[0092] 具体的，粒子群算法的粒子更新方程如下：
[0093]
[0094] Px(T+1)=Vx(T+1)+?!(T);
[0095] 其中，
，对o -1，2, ? ??，T
，对 1=1，2,...，L;0(T)表示惯性权重，计算过程为e(T) =e(i)+g「e(i)-e(i)i，其中0⑴和 0 (I)表示第1和第I代的惯性权重；Cp和c 别表示个体和全局最优学习因子；rp和r1 表示0到1之间的随机数。
[0096] 具体的，系统采用Boltzmann选择算子代替传统全局最优粒子的更新机制。
[0097] 本发明的有益效果是：本发明针对河流流量这一不确定参数，采用模糊随机数进 行描述，并采用两步法进行处理，即首先将其转化成为模糊数，然后采用模糊期望值算子将 模糊数转化为确定值。针对二层决策模型的求解，根据在污染税政策下水污染负荷分配的 二层决策结构特征以及河流系统中污染物的传递特性，提出了混合嵌套基于多智能体的动 态极值算法的改进粒子群算法（PS0)的方法进行求解，从而得到合理的流域污染税费方案 和污染物削减方案。
【附图说明】
[0098] 图1为本发明基于污染税的流域污染负荷分配二层决策优化方法实施例中基于 污染税的流域污染分配问题的二层决策模型结构；
[0099] 图2为本发明基于污染税的流域污染负荷分配二层决策优化方法实施例中B0D和 D0转化机制的概念图；
[0100] 图3为本发明基于污染税的流域污染负荷分配二层决策优化方法实施例中河流 系统转态转态的概念图；
[0101] 图4为本发明基于污染税的流域污染负荷分配二层决策优化方法实施例中排污 控制单元污水处理分析；
[0102] 图5为本发明基于污染税的流域污染负荷分配二层决策优化方法实施例中针对 基于污染税的流域污染负荷分配问题的二层决策建模技术流程；
[0103] 图6为本发明基于污染税的流域污染负荷分配二层决策优化方法实施例中粒子 结构表示；
[0104] 图7为本发明基于污染税的流域污染负荷分配二层决策优化方法实施例中粒子 性能的混合嵌套计算流程图；
[0105] 图8为本发明基于污染税的流域污染负荷分配二层决策优化方法实施例中更新 全局最优粒子的Boltzmann选择算子流程图；
[0106] 图9为本发明基于污染税的流域污染负荷分配二层决策优化方法实施例中河流 系统图。
【具体实施方式】
[0107] 下面结合附图及实施例详细描述本发明的技术方案：
[0108] 本发明针对基于污染税的流域污染负荷分配问题，提出了基于污染税的流域污染 负荷分配二层决策优化方法，首先，系统获取水流量，并将其转换为水流量的确定性变量； 其次，系统获取河流中溶解氧浓度关系及生化需氧量浓度关系；再其次，系统根据水流量的 确定性变量、溶解氧浓度关系及生化需氧量浓度关系构建下层模型；然后，系统根据污染税 收标准制定上层模型；再然后，系统根据下层模型及上层模型构建基于污染税的流域污染 负荷分配二层决策模型；最后，系统对基于污染税的流域污染负荷分配二层决策模型进行 求解，得到最优解。针对河流流量这一不确定参数，采用模糊随机数进行描述，并采用两步 法进行处理，即首先将其转化成为模糊数，然后采用模糊期望值算子将模糊数转化为确定 值。针对二层决策模型的求解，根据在污染税政策下水污染负荷分配的二层决策结构特征 以及河流系统中污染物的传递特性，提出了混合嵌套基于多智能体的动态极值算法的改进 粒子群算法（PS0)的方法进行求解。从而得到合理的流域污染税费方案和污染物削减方 案。
[0109] 实施例
[0110] 本发明主要研宄单一河流系统中污染税政策下的水污染负荷分配设计。单一河流 系统由地方环保部门（EPA)以及一系列坐落于同一流域的排污控制单元（Dischargers)组 成。环保部门为了衡量各排污者污染水排放情况，在其下游安置水质监测点以监测河流水 质情况，并根据流域污染物对人类和生态系统健康的影响，以及流域环境生态功能分区情 况，在各监测点制定相应的水质标准。在满足水质标准条件下，每个排污者被允许将部分污 染物削减后的污水排放入河流中。在污染税政策下，环保部门为化解上下游排污者之间水 污染负荷分配的不公平以及进一步限制各排污者的污水排放，在各个监测点制定相关污染 物污染税费标准，如果某监测点监测相应污染物浓度大于该标准，环保部门将向相应的排 污控制单元征收相应的污染税收。每个排污控制单元需要平衡相应的污水处理成本和污染 税费比例，决策污水处理的削减比例以实现最小成本达到相应水质标准。
[0111] 一般情况下，即使地方环保局和各排污控制单元是独立的决策主体，地方环保局 作为领导者具有优先决策权，所以各排污控制单元是从属者。因此，基于污染税的流域污染 负荷分配问题可以抽象成为一个合作型二层多从属者规划问题，整个过程可以分为两个环 节：污染税制定环节和污染物削减环节：在污染税制度环节，地方环保部门（领导者）制定 T个监测点的污染税收标准以获得最小的所有排污控制单元的单位排污成本的偏差。在污 染物削减环节：各排污控制单元（从属者）依据环保部门制定的税收标准以及上游水质情 况依次决策自身污水削减比例以获得最小的单位排污成本，如图1所示。
[0112] 为解决上述基于污染税的流域污染负荷分配问题，提出如下合作型二层多从属者 决策建模技术，具体步骤如下：
[0113] 1.不确定变量的量化与处理。
[0114] 在水质量管理问题中需要考虑不确定问题已经得到广泛的认知，因为不确定存在 于各种不同的系统组成部分和各种预决策环境当中。因此，现实水质量管理决策中的内部 复杂以及随机性使其不能单单的用传统的确定性优化方法来考虑与解决。因此，如要进一 步对决策中的不确定变量进行量化和处理。
[0115] (1)量化过程：在污染负荷分配决策中，对上游水流量的量化十分重要，然而，上 游水量是一个典型的不确定性，其大小受到多种因素影响，如温度、降雨量、蒸发量等。因 此，本发明采用模糊随机变量对未来时间水流量进行量化，即可以采取相应的水资源配置 决策由于河流水流量的变化，未来的水供应是不确定的，水流量一般可以分为三个不同的 级别：低水位（L)，中水位（M)和高水位（H)，服从概率分布（a，b，c)，三种水位的概率分别 可以表示为&，pjPpH。然而，对于每一个水位，区域水资源管理局很难给定一个确切的值， 因此需要用模糊概念来描述它。为了搜集每一个水位的数据，针对水位k，找相应的水文专 家进行数据采集，他们通常用口头表述"水流量在ak和ck之间，并且有一个最有可能的值 bk。"其中，k表示水位（t=L表示低水位，t=M表示中水位，t=H表示高水位）。值得 一提的是，这些专家给定的口头描述是通过多次的观察的数据获得。因此在每一个水位上 的水流量可以用三角模糊数描述，综上，水流量（可供的分配水量）可以考虑成为一个模糊 随机变量。
[0117] (2)处理过程：一般地，由于模糊随机变量的存在我们很难获得模型的解。本发明 提出一种两步转化方法，将模糊随机变量转化成为确定型变量。第一步：将模糊随机变量转 化成为模糊数；第二步：运用带乐观-悲观指数的期望值算子将模糊数转化成为确定型变 量。关于上游水流量^ ,具体的转换步骤如下：
[0118] ①假设水流量可以表述成为任意一个模糊随机变量x，含有以下数值：
[0120] 其中，A，0、,,，0"表示三个水流量水平，因此其概率分布如下：
[0122] 其中，Q表示三种水流水平及其模糊逻辑；
[0123] 假设：0 (x)是模糊随机变量x的模糊参数，如果0 (x) =Ef[X|P]，e[0,1] 存在，则有：
[0126] 那么，模糊随机变量x的模糊期望值可以表示为f兩]，则有：
[0128]其中，aE=inf( 0 (x)) 0，bE=inf( 0 (x))i=sup( 0 (x)) "cE=sup( 0 (x)) 0， E'[§]表示模糊随机变量x的模糊期望值。
[0129] ②模糊随机水流量被转化成为三角模糊数，系统采用的带乐观-悲观指数的期望 值算子将其转化成为确定性值；
[0130] 首先，定义一个乐观-悲观指数（〇<e< 1)用于量化决 策者的决策偏好，然后 使用期望值算子计算该三角模糊数的期望值，并将其作为水流量的确定性变量，则有：
[0132] 其中，EIE1^]]表示水流量的确定性变量。
[0133] 2.河流同化能力模拟。
[0134] 河流同化能力表示一定的水域、水量条件下，基于维持水域一定的含氧水平，该水 域对有机污染物的自净量。其与河段的流量、流速、水温、允许的最低溶解氧含量及污染物 本身的特性有关。确定水体对污染物的同化能力是确定地区有机污染物排放标准以及制度 合理的污染负荷分配方案的基础。本发明使用Streeter-Phelps模型来模拟河流对有机污 染物的同化过程。其主要原理为：当河流受纳有机物后，沿水流方向产生的输移有机物远大 于扩散稀释量，当河水流量与污水流量稳定，河水温度不变时，有机物生化降解的好氧量与 该时期河水中存在的有机物量成正比。因此，有两个微分耦合模型被定义来描述河流中溶 解氧0〇)和生化需氧量（B0D)浓度变化规律。
[0135] (1)生化需氧量（B0D)转换机制：生化需氧量表示水中有机物等需氧污染物质含 量的一个综合指标。说明水中有机物由于微生物的生化作用进行氧化分解，使之无机化或 气体化时所消耗水中溶解氧的总数量。其值越高说明水中有机污染物质越多，污染也就越 严重。根据Streeter-Phelps模型的结论，河流中B0D浓度的随时间变化量与其自身浓度成 正比，即
，其中Cb表示B0D浓度，kd表示脱氧系数。因此，在一个河段t，用kd(t) 表示脱氧系数，d(t)表示该河段长度，v(t)表示该河段平均水流流速，Cb(t)表示流入该河 段水流的B0D浓度，基于以上微分模型的解析解，则流出该河段水流的B0D浓度Cb(t+1)为
，该浓度等价于：
[0136]Cb(t+1) =a(t)Cb(t) 〇
[0137] (2)溶解氧（DO)转换机制：水中溶解氧的多少是衡量水体自净能力的另一个 指标。根据Streeter-Phelps模型的结论，河流中的溶解氧从空气中补充，其补充浓度 的随时间变化量与其氧亏（即水体中饱和溶解氧和现存溶解氧的差）浓度成正比，即
1其中Cb表示B0D浓度，D。表示氧亏浓度，kd，匕表示分别脱氧和充氧系 数。与B0D情况相似，在一个河段1用ka(t)表示充氧系数，Cb(t)，Cjt)分别表示流入该 河段水流的B0D浓度，Cd。表示该河段水流的饱和溶解氧浓度，则流出该河段水流的D0浓度 C0(t+1)为：
[0139] 具体如图2所示。
[0140] 定义
和S⑴=CdJl-|3⑴]，该浓度等价于：
[0141]C0(t+1)=|3 (t)C0(t)-y(t)Cb(t)-8 (t)〇
[0142]其中，0⑴表示河段t的溶氧浓度的转换系数，Y(t)、S⑴分别表示影响河段 t溶氧浓度的生化需氧量浓度、饱溶氧浓度的转换系数。
[0143] 3.下层模型构建。
[0144] 对于下层模型，从河流上游到下游，每个排污控制单元为了达到最小的污染成本， 根据地方环保部门制定的污染税收标准以及上游的水质情况依次进行独立的污水B0D削 减决策。
[0145] (1)河流系统状态方程组构建。
[0146] 河流系统中的水以及水中的各类化学物质沿着河道顺流而下，在流经过各个的排 污控制单元过程中，控制单元取水、排污等活动会一定程度的影响河流状态（水量、各类污 染物浓度等）的改变，并且该状态会随水的流动转移到下一个控制单元，如图3所示。因 此，在污染负荷分配问题中，需关注河流状态在各个排污之间的转移情况，即流量状态转移 方程、B0D浓度状态转移方程、D0浓度状态转移方程被构建。
[0147] ①流量状态转移方程。基于质量守恒定律以及流量连续性条件，用Qjt-1)，Qjt) 分别表示在监测点t-1和t的流量，用Qw(t)，Qd(t)分别表示排污控制单元t的取水量和排 污量，则可构建流量状态转移方程如下：
[0148]Qr (t) =Qr (t-1) -Qw (t) +Qd (t) 〇
[0149] ②BOD浓度状态转移方程。基于质量守恒定律、Streeter-Phelps水质模拟模型 以及流量状态转移方程，用Cb(t_l)，Cb(t)表示监测点t-1和t的B0D浓度，q⑴表示排污 控制单元t未处理的污水B0D浓度，以及n(t)为排污控制单元t进行削减污水B0D浓度 的比例，其为下层排污者t的决策变量，则可构建B0D浓度状态转移方程如下：
[0151]定义变量
，则该状态转移方程等价于：
[0153] ③D0浓度状态转移方程。相似于B0D情形，基于质量守恒定律、Streeter-Phelps 水质模拟模型以及流量状态转移方程，用（；(t-l)，Cjt)表示监测点t-1和t的DO浓度，C.,⑴表示排污控制单元t未处理的污水D0浓度，则可构建D0浓度状态转移方程如下：
[0155] 同样的定义变量
，则该状态转移方程等价于：
[0157] 另外，考虑河流中饱和溶解氧浓度Cd。的限制，以上D0浓度状态转移方程可改写 为：
[0159] (2)初始模型河流状态条件构建。
[0160] 对于位于整个河流系统中上游第一个排污控制单元，用C:，C:分别表示上游水流 B0D和D0浓度，则该控制单元的模型考虑初始河流状态条件如下：
[0162] (3)下层模型约束条件构建。对于下层模型的决策者-排污控制单元，在进行相应 污水B0D浓度削减的决策过程中，需满足一下两个约束条件：允许的B0D和D0浓度标准以 及B0D削减能力约束。
[0163] ①允许的B0D和D0浓度标准。排污控制单元t排放污水到河流中，不能使得在相 应检测点河流的B0D浓度（Cb(t))超过环保部门所允许的最大B0D浓度标准（Bu(t))，即：
[0164] 0 ^Cb(t) ^Bu(t) 〇
[0165] 同样，不能使相应检测点河流的DO浓度（C；(t))低于环保部门所允许的最小DO浓 度标准汍⑴），即：
[0166] Dx(t) ^C0(t)〇
[0167] ②BOD削减能力约束。每个排污控制单元具有污水BOD浓度削减能力的限制，通 常采用二级污水处理工艺，如图4所示。对于排污控制单元t，第一级处理后的污水B0D浓 度削减率为，第二级处理后的污水bod浓度削减率为nr。在实际的污水处理运行中，环 保部门要求所以排污者的污水必须进行第一级处理。因此，排污控制单元t进行污水B0D 浓度削减决策，其B0D削减率（n(t))不能够超过nr，同时不能低于<in，BP:
[0169] (4)下层模型目标函数构建。下层决策者的目标为最小的单位污染成本目标函数， 污染成本由污水处理成本和污染税费组成。
[0170] ①污水处理成本。基于B0D削减和污水处理成本的经验曲线，如图4所示，本发明 对于排污控制单元的污水处理费用采用幂函数形式，用kel(t)和ke2(t)表示排污控制单元 t的成本参数1和2,Mt)和keb⑴表示排污控制单元t的污水和B0D削减参数，则污水 处理成本C(t)可以表示为以下方程：
[0172] 其次，参数1^(〇，1^2(〇，1^(〇和1^(〇数值以污水处理构筑物建设和运营所得 数据，采用最小二乘法进行曲线拟合所得。
[0173] ②污染税费。在每个监测点，地方环保部门制定相应的污水税收标准（Alb(t)，上 层模型的决策变量）和税率（Tb)。如果检测点t所对应的河流B0D浓度（Cb(t))低于该税 收标准（Cb(t)彡Alb(t))，则排污控制单元t不需要征收污染税（P(t));如果检测点t所对 应的河流B0D浓度（Cb(t))超过该税收标准（(；(〇>八113(〇)，则排污者七需征收污染税？(〇 =(t)Tb [Cb (t)-Alb⑴]。因此，排放者t的污染税费模型定义为：
[0175] 综上所述，排污控制单元t的目标函数，即最小化的单位污染成本如下：
[0177] 4?上层模型（地方环保部门）构建。
[0178] 对于上层模型，地方环保部门为实现公平的流域污染负荷分配，于河流系统中每 个监测点制度相应的污染税收标准。
[0179] (1)上层模型约束条件构建。对于上层模型的决策者-地方环保部门，在进行相应 污染税收标准的制定过程中，需满足允许的BOD浓度标准约束条件。
[0180] 环保部门所制定的污染税收标准（Alb(t))不能超过允许的最大B0D浓度标准 (Bu(t))，即：
[0181] 0 ^Alb(t) ^Bu(t) 〇
[0182] (2)上层模型目标函数构建。上层决策者需保证公平的污染负荷分配，因此，本 发明采用各排污控制单元的单位污染成本的方差来量化水污染负荷分配的公平性。用
表示各排污控制单元平均的单位污染成本，则其最小的方差定义为：
[0184] 5.全局模型
[0185] 集成上下层目标函数和上下层约束条件，我们可以得到以下基于污染税的流域污 染负荷分配问题的合作型二层多从属者规划全局模型。
[0187]
[0188] 如图5所示，表示了针对基于污染税的流域污染负荷分配问题的合作型二层多从 属者决策优化建模技术流程图。
[0189] 针对上述合作型二层多从属者决策模型，本发明提出了混合嵌套基于多智能体的 动态极值算法的改进粒子群算法进行求解。由于二层规划问题具有较为特殊的决策规则， 即上层决策者只是通过自己的决策去指导下层决策者，并不直接干涉下层的决策；而下层 决策者只需要把上层的决策作为参数，他可以在自己的可能范围内自由决策。这种决策机 制使得上层决策者在选择策略以优化自己的目标达成时，必须考虑到下层决策者可能采取 的策略对自己的不利影响。其求解问题已经为证明为NP-难题，即使是最简单的情况。然 而，上述模型中存在非线性的目标函数和约束条件，增加模型计算的难度；另一方面，该模 型相比于传统的二层模型具有更为特殊的动态决策结构，即下层模型需要通过河流系统状 态转移方程进行依次计算，显然，使得该模型的求解问题变的更加困难。因此，本发明针对 基于污染税的流域污染负荷分配二层决策模型，将传统精确算法与现代启发式智能算法相 结合，设计了一种混合 嵌套粒子群算法，具体包括：①针对上述下层模型特殊的多元动态决 策结构，设计一种基于多智能体的动态极值算法为精确快速的获得基于特定上层决策变量 的下层模型的最优解；②为计算上层模型，将该算法嵌套至粒子群算法的各个粒子适应值 计算的过程中；③为提高整体算法的全局搜索能力，采用Boltzmann选择算子代替传统全 局最优粒子的更新机制。下面介绍其主要设计思路。
[0190] 1.粒子结构设计。为简化描述，以决策排污者为T为例介绍算法的粒子结构设计。 在混合嵌套粒子群优化算法中，每个粒子的位置编码为T维的向量，以表示解空间中的一 点。每个粒子位子代表上层模型一个决策变量（Alb (t)，t= 1，2,. . .，T)，如下
[0192] 其中，t表示迭代次数指标，t= 1，2, ...，I，I表示最大迭代次数，1表示粒子指 标，1 = 1，2, ...，L，L表示种群指标。图8表示了一个包括T个排污控制单元的河流系统 的粒子结构表达。
[0193] 2.粒子性能的混合嵌套计算。基于二层规划特殊的决策规则，在使用粒子群算法 进行上层模型的计算的过程中，对各代粒子性能的计算时，需要获得基于特定上层决策变 量的下层模型的最优解。因此，针对基于污染税的流域污染负荷分配问题的合作型二层多 从属者决策模型中下层模型的多元动态特征，本发明为求得其最优反应函数，设计了一种 基于多智能体的动态极值算法。
[0194] (1)基于多智能体的动态极值算法。为求解在特定上层决策变量下的上述模型的 下层模型，需求解T个一元的非线性规划模型。因此，针对单一的非线性规划模型，本发明 提出了一种两阶段极值算法，即第一阶段通过极值算法求解该问题极值，第二阶段通过约 束处理方法从第一次获得极值中找到该问题的最优解。在解决单一下层模型求解问题后， 本发明基于下层各个模型间的动态传递关系，构建了一种多智能体动态框架，以获得下层T 个模型的全部最优解。
[0195] ①极值算法。以下层第t个单一模型为例，介绍极值算法求解下层单一模型 的极值。令下层模型目标函数中的Cb(t) =Alb(t)，我们能得到该目标函数的断点为
因此，该模型的极值存现两种情况：
[0196] 情况I当Cb(t)彡Alb(t)时，单一下层模型的目标函数如下：
[0198]
该目标函数等价于
[0199] 巧=atri(t)p'+ct
[0200] 因此，在该种情况下，目标函数的极值点为nel(t) =np(t)。
[0201] 情况II当Cb(t)>Alb(t)时，单一下层模型的目标函数如下：
[0202]
[0203] 令1^ =TBCUt)，该目标函数等价于：
[0205]因此，在该种情况下，该目标函数的极值点为
>考虑极值点的范 围，对该种情况下目标函数的极值点做以下修正：当ne2(t)彡np(t)时，下层模型目标 函数的极值点为nel(t) =np(t);当o〈ru2(t)〈np(t)，下层模型目标函数的极值点为
。该极值算法的具体步骤如表1所示：
[0206] 表1:极值算法的步骤：
[0208]
[0209] ②约束处理方法。在计算下层单一模型目标函数的极值后，需对下层目标进行约 束处理，通过对下层模型相应的约束条件的等价转换，允许的B0D浓度标准约束变为：
[0211]
该约束可以等价于：
[0213] 允许的BOD浓度标准约束变为：
[0215] 令gt等于上述不等式的右边，该不等式二次转换为：
[0217]
该约束可以等价于：
[0219] 因此，下层单一模型约束处理方法的具体步骤如表2所示。
[0220] 表2,下层单一模型约束处理方法的具体步骤：
[0223] ③多智能体动态框架。基于下层模型之间河流系统动态转移方程，本发明设计了 一种基于多智能体的动态框架，结合上述提出的极值算法，以求得在特定上层决策下的全 部下层模型的最优反应集。其具体步骤如表3所示。
[0224] 表3基于多智能体的动态极值算法的求解过程；
[0225]
[0227] (2)适应值函数。通过上述所提出的算法能够计算获得在特定上层决策变量下 的下层模型的最优解，根据双层模型的决策规则，可以得到量化例子性能的适应值函数，即 为：
[0229] 其中，
。因此，整个粒子性能的混合嵌套计算过程可以描述为图7所 示的计算方法。
[0230] 3.粒子更新机制。基于粒子群算法的运行程序，在量化粒子性能后需进行粒子更 新，为获得更好性能的粒子。粒子更新方程如下：
[0232] 卩丄（t+1) = % (t+1)+Pi(t);
[0233]其中
表不当前个体最优粒子，即Fitness^13)SFitness^。)），对o =1，2,…，T;Pe=[PlQ(l)，p「(2)，...，pf(T)]表示当前粒子群中最优粒子，即FitneSS(P°) 莳腦8(杧），对 1 = 1，2，...，L; 0 (t)表示惯性权重，计算过程为
，其中0 (1)和 0 (I)表示第1和第I代的惯性权重；cp和c 别表示个体和全局最优学习因子；rp和r1 表示0到1之间的随机数。为提供粒子的全局搜索能力，本发明采用Boltzmann选择算子 代替传统全局最优粒子的更新机制。
[0234] (1)Boltzmann选择算子。采用Boltzmann选择算子更新全局最优粒子，其基本思 路：每代全局最优粒子1><3=[1 )「(1)，?「(2)，...,1)「(1')]，通过运用13〇1〖21]^1111选择算子在各粒子的个 体最优粒子矿=[pf(l)，pf(2)，...，pf(T)]中进行选择。为计算各个粒子的个体最优粒子的选择概 率，其评价函数定义如下：
[0236] 其各粒子的个体最优粒子的选择概率为
[0238] 其中TejPTe:表示初代和第t代温度。因此，运用Boltzmann选择算子进行各 粒子的全局最优粒子的具体步骤如表4及图8所示。
[0239] 表4运用Boltzmann选择算子进行各粒子的全局最优粒子的具体步骤；
[0240]
[0242] 4.混合嵌套粒子群算法。该算法的整体步骤如表5所示：
[0243] 表5混合嵌套粒子群算法步骤；
[0245]
[0246] 为了能进一步了解本发明的
【发明内容】
、特点及用法，兹举以如下实例详细说明。
[0247] 实例：四川省沱江流域某河段为例，其具有5个排污控制单元= 〇, 1，2, 3, 4, 5)，整个区域的考虑基于污染税政策下的污染负荷分配决策，该河流系统见图 9。通过对该流域水位数据访问调研，该流域饱和溶解氧浓度为10. 08mg/L，其他的水位数据 见表6。对于各排污控制单元，其污水处理成本参数（kel⑴，ke2 (t)，Mt)和(t))数值 以污水处理构筑物建设和运营所得数据，采用最小二乘法进行曲线拟合所得，其他相关数 据见表7。对于河流系统的上游情况，其B0D和D0浓度为3. 2mg/L和6. 2mg/L，其水流量为 模糊随机变量，其分布数据见8。假设上层决策的偏好乐观-悲观指数设置为e=0.5,通 过计算上游水流量的期望值为12. 10m3/s。该河流系统各监测点的环境水质标准见表9,另 外基本污染税率为〇? 005yuan?L/mg。
[0248] 表6河段水位数据；
[0250]
[0251] 表7排污控单元数据；
[0253] 表8河流水流量分别数据；
[0254]
[0255]表9河流系统各监测点的环境水质标准
[0257] 计算结果。应用本发明提出的基于污染税的问题的多阶段建模技术对该实例问题 进行建模，并采用本发明提出的两步转化方法对模型中的三角模糊随机数进行处理，最后 应用本发明提出的混合嵌套粒子群算法对模型进行求解计算，结果如表10所示。通过将算 法优化的结果反映，各排污控制单元的污染成本都为1. 72yuan?s/m3。因此，地方环保部 门通过污染税收的制定，实现了公平的污染负荷分配。这种分配结果可以能够实现流域健 康可持续的发展。尽管算法优化的结果是基于模型假设的条件得出的，但在一定程度上仍 能反映出实际决策方案尚有提升的空间。在实际决策过程中，各决策者仅仅根据经验判断 决定决策方案具有一定风险。运用模型和算法计算的结果可对决策者进行方案选择提供参 考，以更好地指导水质管理实践。
[0258] 表10模型求解结果
[0259]
【主权项】
1. 基于污染税的流域污染负荷分配二层决策优化方法，其特征在于，包括以下步骤： 步骤1、系统获取水流量，并将其转换为水流量的确定性变量； 步骤2、系统获取河流中溶解氧浓度关系及生化需氧量浓度关系； 步骤3、系统根据水流量的确定性变量、溶解氧浓度关系及生化需氧量浓度关系构建下 层模型； 步骤4、系统根据污染税收标准制定上层模型； 步骤5、系统根据下层模型及上层模型构建基于污染税的流域污染负荷分配二层决策 丰旲型； 步骤6、系统对基于污染税的流域污染负荷分配二层决策模型进行求解，得到最优解。2. 根据权利要求1所述的基于污染税的流域污染负荷分配二层决策优化方法，其特征 在于，所述步骤1中，首先，系统采用模糊随机变量对未来时间水流量进行量化，系统将水 流量分为三个不同的级别：低水位（L)，中水位（M)和高水位（H)，服从概率分布（a，b，c)， 三种水位的概率分别可以表示为Py Pm和P H，综上，水流量$为：其次，将水流量^转换为模糊随机变量，并计算其模糊期望值，具体包括： ① 假设水流量可以表述成为任意一个模糊随机变量X，含有以下数值： Ol _ (aL'bL'CL)，Qm _ (aM'bjvI，CM)，Qh _ (aH，bH'CH). 其中，I，，O11表示三个水流量水平，因此其概率分布如下： Ρ({ω e Ω I χ(ω) = QL}) = pL, Ρ({ω e Ω | χ(ω) = Qm }) = ρΜ, Ρ({ω ε Ω | χ(ω) = Qh }) = ρΗ 9 其中，Ω表示三种水流水平及 其模糊逻辑； 假设：Θ (X)是模糊随机变量X的模糊参数，如果Θ (X) =Ef [X I P]，Vae [0，1]存在， 则有：那么，模糊随机变量X的模糊期望值可以表示为1^造1，则有： E'[Q] = (aE,bE,cE) _ 其中，aE= inf ( θ ( x )) Q, bE= inf ( θ ( x )) i = sup ( θ ( x )) Ce= sup ( θ ( x )) 〇， 表示模糊随机变量x的模糊期望值； ② 模糊随机水流量被转化成为三角模糊数，系统采用的带乐观-悲观指数的期望值算 子将其转化成为确定性值； 首先，定义一个乐观-悲观指数（〇< ε <1)用于量化决策者的决策偏好，然后使用 期望值算子计算该三角模糊数的期望值，并将其作为水流量的确定性变量，则有：其中，Ε^Ε?^]]表示水流量的确定性变量。3. 根据权利要求1所述的基于污染税的流域污染负荷分配二层决策优化方法，其特征 在于，所述步骤2中，系统根据Streeter-Phelps模型的结论，河流中生化需氧量浓度的随 时间变化暈与其自身浓度成正比，即：其中，Cb表示生化需氧量浓度，k 示脱氧系数； 因此，在一个河段表示脱氧系数，d(t)表示该河段长度，v(t)表示该河段 评价水流流速，Cb(t)表示流入该河段水流的生化需氧量浓度，则流出该河段水流的生化需 氧量浓度Cb (t+Ι)为，该生化需氧量浓度等价于： Cb(t+1) = a (t)Cb(t)； 系统根据Streeter-Phelps模型的结论，河流中的溶解氧从空气中补充，其补充浓度 的随时间变化量与其氧亏浓度成正比，即：其中，Cb表示生化需氧量浓度，D。表示氧亏浓度，1^表示脱氧系数，1表示充氧系数。 与生化需氧量情况相似，在一个河段1用ka(t)表示充氧系数，Cb (t)，CJt)分别表示流入 该河段水流的生化需氧量浓度，Cd。表示该河段水流的饱和溶解氧浓度，则流出该河段水流 的溶解氧浓度cjt+l)为：定义：-和δ (t) = CdJl-β⑴]，该溶解氧浓度等 价于： C0(t+1) = β (t)C0(t)-y (t)Cb(t)-5 (t)； 其中，a (t)、β (t)分别可以表示河段t的生化需氧量浓度、溶氧浓度的转换系数， γ (t)、δ (t)分别表示影响河段t溶氧浓度的生化需氧量浓度、饱溶氧浓度转换系数。4. 根据权利要求1所述的基于污染税的流域污染负荷分配二层决策优化方法，其特征 在于，所述步骤3中，系统构下层模型包括： (1)河流系统状态方程组构建，具体包括： ① 流量状态转移方程：基于质量守恒定律以及流量连续性条件，用A(t-l)，Qjt)分别 表示在监测点t-Ι和t的流量，用Qw (t)，Qd (t)分别表示排污控制单元t的取水量和排污 量，则可构建流量状态转移方程如下： Qr (t) = Qr(t-1)-Qw(t)+Qd(t)； ② BOD浓度状态转移方程：基于质量守恒定律、Streeter-Phelps水质模拟模型以及流 量状态转移方程，用Cb (t-1)，Cb (t)表示监测点t-Ι和t的BOD浓度，（'Ut)表示排污控制单 元t未处理的污水BOD浓度，以及n (t)为排污控制单元t进行削减污水BOD浓度的比例， 其为下层排污者t的决策变量，则可构建BOD浓度状态转移方程如下： Cb (t)Qr(t) = a(t- l)Cb (t - l)[Qr (t -1) - Qw (t)]+(t)[l - n(t)]Qd(t); 定义变I?则该状态转移方程等价于： Cb (t) = LI - r(t)ja(t - l)Cb (t -1) + r(t)C；1 (t)[l - η(?)]; ③DO浓度状态转移方程：相似于BOD情形，基于质量守恒定律、Streeter-Phelps水质 模拟模型以及流量状态转移方程，用C。(t-1)，C。⑴表示监测点t-Ι和t的DO浓度，C〖(t)表 示排污控制单元t未处理的污水DO浓度，则可构建DO浓度状态转移方程如下： C0(t)Qr(〇 = LP(t-l)Q(t-l)-T(t-l)Cb(l-l)-6(l-l)JLQI(t-l)-Q w(i)J+C· (I)Qd(I)： 同样的定义变量I则该状态转移方程等价于：另外，考虑河流中饱和溶解氧浓度Cd。的限制，以上DO浓度状态转移方程可改写为： C^t) = inin{[l-r(t)][p(l-l)C；At-l)-y(t-l)Cb(t-l)-6(t-l)]+r(t)CUi).C dJ ； (2) 构建初始模型河流状态条件： 对于位于整个河流系统中上游第一个排污控制单元，用C；；，C:分别表示上游水流BOD 和DO浓度，则该控制单元的模型考虑初始河流状态条件如下： Qr(O) = Εε[Ε'[4]], Cb(O) = Q, C0(O) = C^ ； (3) 构建下层模型约束条件： ① 允许的BOD和DO浓度标准： 排污控制单元t排放污水到河流中，不能使得在相应检测点河流的BOD浓度（Cb(t))超 过环保部门所允许的最大BOD浓度标准（Bu (t))，即： 0. Cb (t) ^ Bu (t)； 同样，不能使相应检测点河流的DO浓度（C；(t))低于环保部门所允许的最小DO浓度 标准（D1⑴），即： D1 (t) ^ C0 (t)； ② BOD削减能力约束： 每个排污控制单元具有污水BOD浓度削减能力的限制，通常采用二级污水处理工艺。 对于排污控制单元t，第一级处理后的污水BOD浓度削减率为η;111π，第二级处理后的污水BOD 浓度削减率为ηΓ，在实际的污水处理运行中，环保部门要求所以排污者的污水必须进行第 一级处理；因此，排污控制单元t进行污水BOD浓度削减决策，其BOD削减率（τι (t))不能 够超过ηΓ，同时不能低于ηΓ，即： OWWt; (4) 构建下层模型目标函数： 下层决策者的目标为最小的单位污染成本目标函数，污染成本由污水处理成本和污染 税费组成； ① 污水处理成本： 基于BOD削减和污水处理成本的经验曲线，本发明对于排污控制单元的污水处理费用 采用幂函数形式，用L⑴和U⑴表示排污控制单元t的成本参数1和成本参数2, Mt) 和Ι^α)表示排污控制单元t的污水和BOD削减参数，则污水处理成本C(t)可以表示为以 下方程：其次，参数ha)，u(t)，Mt)和u⑴数值以污水处理构筑物建设和运营所得数 据，采用最小二乘法进行曲线拟合所得； ② 污染税费： 在每个监测点，地方环保部门制定相应的污水税收标准（Alb(t)，上层模型的决 策变量）和税率（Tb)，如果检测点t所对应的河流BOD浓度（Cb(t))低于该税收标准 (Cb(t) <Alb(t))，贝IJ排污控制单元t不需要征收污染税（P(t));如果检测点t所对应的 河流如0浓度（(；(〇)超过该税收标准（(；(〇>4113(〇)，则排污者〖需征收污染税？(〇 = Qi (t) Tb [Cb (t)-Alb⑴]，因此，排放者t的污染税费模型定义为：综上所述，排污控制单元t的目标函数，即最小化的单位污染成本如下：5. 根据权利要求1所述的基于污染税的流域污染负荷分配二层决策优化方法，其特征 在于，所述步骤4中，系统构上层模型包括： (1) 上层模型约束条件构建： 对于上层模型的决策者-地方环保部门，在进行相应污染税收标准的制定过程中，需 满足允许的BOD浓度标准约束条件，环保部门所制定的污染税收标准（Alb(t))不能超过允 许的最大BOD浓度标准（Bu⑴），即： 0 ^Alb (t) ^ Bu (t)； (2) 上层模型目标函数构建： 上层决策者需保证公平的污染负荷分配，因此，本发明采用各排污控制单元的单位污 染成本的方差来量化水污染负荷分配的公平性，用表示各排污控制单元平均的 单位污染成本，则其最小的方差定义为：6. 根据权利要求1所述的基于污染税的流域污染负荷分配二层决策优化方法，其特征 在于，所述步骤5中，系统根据下层模型及上层模型构建的基于污染税的流域污染负荷分 配二层决策模型为：7. 根据权利要求1所述的基于污染税的流域污染负荷分配二层决策优化方法，其特征 在于，所述步骤6中，假设决策排污者为T个，则在特定上层决策变量下地下层模型为T个 一元的非线性规划模型；首先，系统利用极值算法计算出单一的非线性规划模型的极值，再 通过下层约束条件计算出该单一的非线性规划模型的最优解，进而获得下层T个一元的非 线性规划模型的全部最优解集；利用粒子群算法对上述最优解集进行计算，从而获得基于 污染税的流域污染负荷分配二层决策模型的最优解。8. 根据权利要求7所述的基于污染税的流域污染负荷分配二层决策优化方法，其特征 在于，粒子群算法的粒子更新方程如下： Ρ^τ+Ι) = V^t+D+P^t)；其中，杧=[?!5(1)，1^(2)，...￥(丁)]表示当前个体最优粒子，即府11 (^叱)<]%^(?1(〇))，对0 = 1，2，···，τ ; Pc=丨pf (1),4(2),..., pf (T)]表示当前粒子群中最优粒子，SPFitness(Pc)i; Fitness(P1B),对 I =1，2，...，L; θ (τ)表示惯性权重，计算过程为，其中Θ⑴和 Θ⑴表不第1和第I代的惯性权重；Cp和％分别表不个体和全局最优学习因子；rp和％表 示O到1之间的随机数。9.根据权利要求7所述的基于污染税的流域污染负荷分配二层决策优化方法，其特征 在于，系统采用Boltzmann选择算子代替传统全局最优粒子的更新机制。
【专利摘要】本发明提出了基于污染税的流域污染负荷分配二层决策优化方法，首先，系统获取水流量，并将其转换为水流量的确定性变量；其次，系统获取河流中溶解氧浓度关系及生化需氧量浓度关系；再其次，系统根据水流量的确定性变量、溶解氧浓度关系及生化需氧量浓度关系构建下层模型；然后，系统根据污染税收标准制定上层模型；再然后，系统根据下层模型及上层模型构建基于污染税的流域污染负荷分配二层决策模型；最后，系统对基于污染税的流域污染负荷分配二层决策模型进行求解，得到最优解。适用于基于污染税的流域污染负荷分配问题。
【IPC分类】G06Q10/04
【公开号】CN104899668
【申请号】CN201510390525
【发明人】徐玖平, 张梦翔, 曾自强
【申请人】四川大学
【公开日】2015年9月9日
【申请日】2015年7月3日
最新回复(0)