污水处理过程动态多目标优化控制方法

　　申请日2015.12.12

　　公开(公告)日2016.03.16

　　IPC分类号G05B13/04

　　摘要

　　污水处理过程动态多目标优化控制方法属于水研究领域，也属于智能控制领域。本发明在保证出水水质达标的前提下，同时优化能耗和出水水质指标。首先，通过神经网络在线建模方法构造污水处理过程多目标优化模型，解决优化变量与性能指标间没有精确数学描述的问题;其次，基于建立的优化模型，由多目标优化算法得到问题的Pareto最优解，并根据决策信息，从Pareto解集中确定一个满意优化解，即溶解氧浓度和硝态氮浓度的优化设定值;最后，由底层控制器实现优化设定值的跟踪任务。该方法在保证出水水质达标前提下，能够有效降低污水处理过程能耗和运行成本。

　　摘要附图

 

　　权利要求书

　　1.污水处理过程动态多目标优化控制方法，其特征在于，包括以下步骤：

　　1).以出水水质达标为约束条件，同时优化能耗和出水水质两个性能指标，构造污水处理过程的约束多目标优化模型;

　　1.1多目标优化模型构造

　　建立优化设定值与性能指标间的函数关系，其实质是建立溶解氧浓度和硝态氮浓度的优化设定值与能耗指标、出水水质指标间的代理模型，实现对性能指标的预测和评价;令x1(k)为溶解氧浓度设定值，x2(k)为硝态氮浓度设定值，x(k)=[x1(k),x2(k)]为设定值组成的优化向量，k为时刻，为方便表述，记x(k)为x;fAE(x)为优化变量与曝气能耗间的函数表达，fPE(x)为优化变量与泵送能耗间的函数表达，能耗指标定义为曝气能耗与泵送能耗之和，即能耗指标fEC(x)=fAE(x)+fPE(x)，fEQ(x)为优化变量与出水水质指标间的函数表达;

　　构建如下多目标优化模型，

　　minF(x)={fEC(x)，fEQ(x)}(1)

　　$s.t.\left\{\begin{array}{c}{g}_1\left(x\right)-4\le 0\\ g_2\left(x\right)-18\le 0\\ x_1^l<x_1\left(k\right)<x_1^u\\ x_2^l<x_2\left(k\right)<x_2^u\end{array}\right.---\left(2\right)$

　　式中，分别为溶解氧浓度优化设定值的下限和上限值，分别硝态氮浓度优化设定值的下限和上限值;g1(x)为出水氨氮浓度SNH与优化设定值间的函数关系，g1(x)-4≤0表示出水氨氮需满足的约束不等式，4(mg/L)为其约束上限值，g2(x)为出水总氮浓度Ntot与优化设定值间的函数关系，g2(x)-18≤0表示出水总氮需满足的约束不等式，18(mg/L)为其约束上限值;溶解氧浓度为污水处理第五分区溶解氧浓度，硝态氮浓度为污水处理第二分区硝态氮浓度;

　　1.2模糊神经网络在线建模

　　优化设定值与曝气能耗、泵送能耗、出水水质的函数关系，采用TS模糊神经网络在线建模方法得到;神经网络的输入向量由z(k)表示，对于曝气能耗模型、泵送能耗模型，网络的输入为k时刻溶解氧浓度和k时刻硝态氮浓度;对于出水氨氮、出水总氮、水质性能指标的函数模型，网络输入为k时刻溶解氧浓度、k时刻硝态氮浓度，以及k时刻入水流量;为方便表述，z(k)简记为z;

　　模糊神经网络的输出表示为

　　$\hat{y}\left(k\right)=\Phi \left(z\right)\&CenterDot;H(z,\theta )---\left(3\right)$

　　其中，Φ(z)=f(A(z))

　　(4)

　　H(z,θ)=[1,zT]·θ(5)

　　式中，A(z)为网络输入经模糊化后得到的隶属度矩阵，f为每条规则隶属度乘积及规一化运算，隶属度函数取高斯函数;θ为θ(k)的简记，表示k时刻神经网络的后件参数矩阵;上标T表示向量的转置;为k时刻网络输出;

　　设k时刻网络学习的目标函数为：

　　$J\left(k\right)=\frac{1}{2}e{\left(k\right)}^2=\frac{1}{2}{\left(\hat{y}\right(k)-y(k\left)\right)}^2---\left(6\right)$

　　采用梯度下降算法更新网络参数，调整公式为：

　　$\alpha (k+1)=\alpha \left(k\right)-\eta \frac{\partial J\left(k\right)}{\partial \alpha \left(k\right)}---\left(7\right)$

　　式中，α(k)为k时刻网络的学习参数向量，包括网络的后件参数、隶属度函数的中心和宽度;α(k+1)为k+1时刻网络学习参数向量;η为网络的学习率，在(0,1)之间取值，y(k)为k时刻系统实际物理量输出，基于污水处理过程基准仿真模型BSM1产生;

　　1.3模型中约束条件处理

　　采用惩罚函数法将建立的约束优化问题转化成无约束优化问题;定义惩罚项

　　fpenalty(x)=max{g1(x)-4,0}+max{g2(x)-18,0}(8)

　　加入惩罚项的能耗指标与水质性能指标表达为

　　$\left\{\begin{array}{c}{f}_{EC}^{\prime }\left(x\right)=f_{EC}\left(x\right)+C\&CenterDot;f_{penalty}\left(x\right)\\ f_{EQ}^{\prime }\left(x\right)=f_{EQ}\left(x\right)+C\&CenterDot;f_{penalty}\left(x\right)\end{array}\right.---\left(9\right)$

　　其中，C为惩罚因子，取定惩罚因子为100000;

　　2).多目标优化算法

　　对建立的污水处理过程多目标优化模型采用NSGA-II进化算法进行求解，在溶解氧浓度和硝态氮浓度的优化设定值范围内，获得多目标优化问题的一组等同优秀的解，即Pareto最优解;NSGA-II进化算法中待优化性能指标为(9)式带有惩罚项的能耗指标f′EC(x)和水质性能指标f′EQ(x);取定搜索空间D为2维，种群规模N为40，最大进化代数M为30;按照标准NSGA-II进化算法运行程序，直至达到最大进化代数M，算法结束，给出污水处理过程多目标优化问题的一组Pareto最优解{Xp,1≤p≤N}，即一组等同优秀的溶解氧浓度和硝态氮浓度的优化设定值;

　　3).智能决策

　　为了实现污水处理过程的闭环控制，在多目标优化算法获得的一组Pareto最优解中，确定当前状态下的一个最优满意解，作为溶解氧浓度和硝态氮浓度的优化设定值;计算Pareto解集中所有解的效用函数值

　　dutility(Xp)=ω1fAE(Xp)+ω2fPE(Xp)(10)

　　选取效用函数最小值所对应的解为最优满意解，即得到k时刻由优化算法确定的溶解氧浓度优化设定值x1(k)和硝态氮浓度优化设定值x2(k);式中，ω1、ω2为决策者确定的性能指标权重，且ω1+ω2=1，;

　　4).底层PID控制

　　将溶解氧浓度设定值x1(k)和硝态氮浓度设定值x2(k)送至底层PID控制器，实现污水处理过程设定值跟踪控制;溶解氧浓度通过第五分区氧气转换系数KLa5进行控制，硝态氮浓度通过内回流量Qa进行控制;按基准定义，氧气转换系数KLa5变化范围取为0

　　定义k时刻溶解氧浓度和硝态氮浓度的系统误差：

　　$\left\{\begin{array}{c}{e}_1\left(k\right)=x_1\left(k\right)-y_1\left(k\right)\\ e_2\left(k\right)=x_2\left(k\right)-y_2\left(k\right)\end{array}\right.---\left(11\right)$

　　式中y1(k)，y2(k)分别为k时刻溶解氧浓度和硝态氮浓度的实际测量值;

　　采用增量式PID控制，k时刻控制增量表达为

　　$\left\{\begin{array}{c}{\mathrm{\Delta u}}_1\left(k\right)=K_{P,1}\&CenterDot;\&lsqb;e_1\left(k\right)-e_1(k-1)\&rsqb;+K_{I,1}\&CenterDot;e_1\left(k\right)+K_{D,1}\&CenterDot;\&lsqb;e_1\left(k\right)-2\&CenterDot;e_1(k-1)+e_1(k-2)\&rsqb;\\ {\mathrm{\Delta u}}_2\left(k\right)=K_{P,2}\&CenterDot;\&lsqb;e_2\left(k\right)-e_2(k-1)\&rsqb;+K_{I,2}\&CenterDot;e_2\left(k\right)+K_{D,2}\&CenterDot;\&lsqb;e_2\left(k\right)-2\&CenterDot;e_2(k-1)+e_2(k-2)\&rsqb;\end{array}\right.---\left(12\right)$

　　k时刻系统控制量表达为

　　$\left\{\begin{array}{c}{u}_1\left(k\right)=u_1(k-1)+{\mathrm{\Delta u}}_1\left(k\right)\\ u_2\left(k\right)=u_2(k-1)+{\mathrm{\Delta u}}_2\left(k\right)\end{array}\right.---\left(13\right)$

　　式中KP,1，KI,1，KD,1分别为溶解氧浓度控制器的比例、积分和微分系数;KP,2，KI,2，KD,2分别为硝态氮浓度控制器的比例、积分和微分系数;取定：KP,1=200，KI,1=15，KD,1=2，KP,2=20000，KI,2=5000，KD,2=400;

　　5).每隔2小时优化周期，进行多目标优化算法和智能决策计算，更新一次溶解氧浓度优化设定值x1(k)和硝态氮浓度优化设定值x2(k)，由底层PID执行跟踪控制任务。

　　2.根据权利要求1所述的污水处理过程动态多目标优化控制方法，其特征在于，考虑固体悬浮物浓度SS低于其约束上限值30(mg/L)，化学需氧量COD低于其约束上限值100(mg/L)或/和生物需氧量BOD5低于其约束上限值10(mg/L)，并与约束上限值一起构成多目标优化模型中的约束条件。

　　说明书

　　污水处理过程动态多目标优化控制方法

　　技术领域

　　本发明提出一种基于神经网络和非占优精英遗传算法NSGA-II的污水处理过程动态多目标优化控制方法，对污水处理过程能耗与水质两个重要且相互冲突的性能指标同时优化，实现好氧区溶解氧浓度和厌氧区硝态氮浓度设定值的动态寻优。污水处理优化控制是水处理领域重要的研究课题，属于水研究领域，也属于智能控制领域。

　　背景技术

　　污水处理过程是高能耗的复杂流程工业系统，其操作连续运行并需要保证出水水质达标。随着国家对污水处理过程节能降耗的要求和人类对环境保护意识的增强，污水处理过程的优化控制引起了广泛关注。污水处理过程优化控制是提高污水处理过程性能、实现节能降耗的有效途径。因此，本发明的研究成果具有广阔的应用前景。

　　污水处理优化控制的目标是在满足出水水质达标的前提下，获得能耗和成本的降低。在污水处理过程优化控制研究初期，单变量优化取得了较为丰硕的成果，并以溶解氧浓度优化为主。为了获得更好的优化性能，污水处理过程多变量优化控制近年来引起了众多学者关注。前置反硝化活性污泥法处理的污水过程，第五分区的溶解氧浓度和第二分区的硝态氮浓度强烈影响污水处理过程的硝化和反硝化水平，直接决定出水水质，同时也是曝气能耗AE和泵送能耗PE高低的重要影响参数。因此，对溶解氧浓度So和硝态氮浓度SNo的设定值进行优化已经成为提高污水处理过程优化性能的重要手段。

　　由于活性污泥处理过程中，存在复杂的物理、生化反应，且各种反应过程相互影响，使得整个系统呈现较强的耦合特性，甚至表现为性能指标的相互冲突，如，能耗与出水水质指标间便具有这种冲突特性。因此，将污水处理优化控制视为本质多目标优化问题更符合污水处理过程优化问题的实质，但这方面的研究工作才刚起步。国内外现有污水处理多目标优化控制方法中，存在的主要问题有：考虑了污水过程具有多目标优化特性，但优化算法依然将其处理成单目标优化问题;以优化处理厂设计参数为主的多目标优化方法，采用离线优化模式，不适于在线优化控制;在模型预测控制框架下考虑控制性能指标的冲突特性，但未考虑水质、能耗等经济指标优化。此外，污水处理过程能耗、水质的多目标优化控制问题，难点还在于能耗-水质-优化设定值模型较难建立。

　　发明内容

　　本发明获得一种基于神经网络和非占优精英遗传算法NSGA-II的污水处理过程多目标优化控制方法，实现污水处理过程溶解氧浓度和硝态氮浓度设定值优化，保证出水水质达标前提下，降低能耗和处理成本;本发明采用如下技术方案及实现步骤，其特征在于，包括以下步骤：

　　1.以出水水质达标为约束条件，同时优化能耗和出水水质两个性能指标，构造污水处理过程的约束多目标优化模型;

　　1.1多目标优化模型构造

　　建立优化设定值与性能指标间的函数关系，其实质是建立溶解氧浓度和硝态氮浓度的优化设定值与能耗指标、出水水质指标间的代理模型，实现对性能指标的预测和评价;令x1(k)为溶解氧浓度设定值，x2(k)为硝态氮浓度设定值，x(k)=[x1(k),x2(k)]为设定值组成的优化向量，k为时刻，为方便表述，记x(k)为x;fAE(x)为优化变量与曝气能耗间的函数表达，fPE(x)为优化变量与泵送能耗间的函数表达，能耗指标定义为曝气能耗与泵送能耗之和，即能耗指标fEC(x)=fAE(x)+fPE(x)，fEQ(x)为优化变量与出水水质指标间的函数表达;

　　构建如下多目标优化模型，

　　minF(x)={fEC(x),fEQ(x)}(1)

　　$\begin{array}{cc}s.t.& \left\{\begin{array}{c}{g}_1\left(x\right)-4\le 0\\ g_2\left(x\right)-18\le 0\\ x_1^l<x_1\left(k\right)<x_1^u\\ x_2^l<x_2\left(k\right)<x_2^u\end{array}\right.\end{array}---\left(2\right)$

　　式中，分别为溶解氧浓度优化设定值的下限和上限值，分别硝态氮浓度优化设定值的下限和上限值;g1(x)为出水氨氮浓度SNH与优化设定值间的函数关系，g1(x)-4≤0表示出水氨氮需满足的约束不等式，4(mg/L)为其约束上限值，g2(x)为出水总氮浓度Ntot与优化设定值间的函数关系，g2(x)-18≤0表示出水总氮需满足的约束不等式，18(mg/L)为其约束上限值;本发明选择出水总氮浓度和出水氨氮浓度做为水质约束条件，原因在于这两个出水浓度指标最为重要，且易出现超标现象;考虑其它出水浓度参数作为约束时，如，固体悬浮物浓度SS低于其约束上限值30(mg/L)，化学需氧量COD低于其约束上限值100(mg/L)，生物需氧量BOD5低于其约束上限值10(mg/L)，本发明方法仍然适用，处理方法如同出水总氮和出水氨氮，即建立相应的优化设定值与各出水浓度参数的函数关系，并与约束上限值一起构成多目标优化模型中的约束条件;各出水浓度参数的约束上限值基于国际水质协会制定的基准给出;本发明中，溶解氧浓度为污水处理第五分区溶解氧浓度，硝态氮浓度为污水处理第二分区硝态氮浓度;

　　1.2模糊神经网络在线建模

　　优化设定值与曝气能耗、泵送能耗、出水水质的函数关系，采用TS模糊神经网络在线建模方法得到，模糊神经网络结构如图1所示;神经网络的输入向量由z(k)表示，对于曝气能耗模型、泵送能耗模型，网络的输入为k时刻溶解氧浓度和k时刻硝态氮浓度;对于出水氨氮、出水总氮、水质性能指标的函数模型，网络输入为k时刻溶解氧浓度、k时刻硝态氮浓度，以及k时刻入水流量;为方便表述，z(k)简记为z;

　　模糊神经网络的输出表示为

　　$\hat{y}\left(k\right)=\Phi \left(z\right)\&CenterDot;H(z,\theta )---\left(3\right)$

　　其中，Φ(z)=f(A(z))

　　(4)

　　H(z,θ)=[1,zT]·θ(5)式中，A(z)为网络输入经模糊化后得到的隶属度矩阵，f为每条规则隶属度乘积及规一化运算，隶属度函数取高斯函数;θ为θ(k)的简记，表示k时刻神经网络的后件参数矩阵;上标T表示向量的转置;为k时刻网络输出;

　　设k时刻网络学习的目标函数为：

　　$J\left(k\right)=\frac{1}{2}e{\left(k\right)}^2=\frac{1}{2}{\left(\hat{y}\right(k)-y(k\left)\right)}^2---\left(6\right)$

　　采用梯度下降算法更新网络参数，调整公式为：

　　$\alpha (k+1)=\alpha \left(k\right)-\eta \frac{\partial J\left(k\right)}{\partial \alpha \left(k\right)}---\left(7\right)$

　　式中，α(k)为k时刻网络的学习参数向量，包括网络的后件参数、隶属度函数的中心和宽度;α(k+1)为k+1时刻网络学习参数向量;η为网络的学习率，一般在(0,1)之间取值，较大的学习率能提高学习速度，但不利于网络收敛，较小的学习率利于网络收敛，但学习速度较慢，基于实验运行情况，取定学习率η为0.1;y(k)为k时刻系统实际物理量输出，基于污水处理过程基准仿真模型BSM1产生，BSM1由国际水质协会与欧盟科学技术与合作组联合提出，是当前国际上公认的污水处理过程控制性能评价的基准，其布局如图2所示;

　　1.3模型中约束条件处理

　　采用惩罚函数法将建立的约束优化问题转化成无约束优化问题;定义惩罚项

　　fpenalty(x)=max{g1(x)-4,0}+max{g2(x)-18,0}(8)

　　加入惩罚项的能耗指标与水质性能指标表达为

　　$\left\{\begin{array}{c}{f}_{EC}^{\prime }\left(x\right)=f_{EC}\left(x\right)+C\&CenterDot;f_{penalty}\left(x\right)\\ f_{EQ}^{\prime }\left(x\right)=f_{EQ}\left(x\right)+C\&CenterDot;f_{penalty}\left(x\right)\end{array}\right.---\left(9\right)$

　　其中，C为惩罚因子，一般取较大正实数，其值越大对不满足约束条件的解的惩罚力度越大，但易产生不可行域附近可行解丢失现象，根据实验运行结果，取定惩罚因子为100000;

　　2.多目标优化算法

　　对建立的污水处理过程多目标优化模型采用NSGA-II进化算法进行求解，在溶解氧浓度和硝态氮浓度的优化设定值范围内，获得多目标优化问题的一组等同优秀的解，即Pareto最优解;NSGA-II进化算法中待优化性能指标为(9)式带有惩罚项的能耗指标f′EC(x)和水质性能指标f′EQ(x);根据本发明中的实际优化问题及实验运行结果，取定搜索空间D为2维，种群规模N为40，最大进化代数M为30;按照标准NSGA-II进化算法运行程序，直至达到最大进化代数M，算法结束，给出污水处理过程多目标优化问题的一组Pareto最优解{Xp,1≤p≤N}，即一组等同优秀的溶解氧浓度和硝态氮浓度的优化设定值;

　　3.智能决策

　　为了实现污水处理过程的闭环控制，在多目标优化算法获得的一组Pareto最优解中，确定当前状态下的一个最优满意解，作为溶解氧浓度和硝态氮浓度的优化设定值;计算Pareto解集中所有解的效用函数值

　　dutility(Xp)=ω1fAE(Xp)+ω2fPE(Xp)(10)选取效用函数最小值所对应的解为最优满意解，即得到k时刻由优化算法确定的溶解氧浓度优化设定值x1(k)和硝态氮浓度的优化设定值x2(k);式中，ω1、ω2为决策者确定的性能指标权重，且ω1+ω2=1，权重大小体现了决策者对某一性能指标的偏重情况，本发明同时强调能耗指标和水质指标，确定ω1=0.5，ω2=0.5;

　　4.底层PID控制

　　将溶解氧浓度设定值x1(k)和硝态氮浓度设定值x2(k)送至底层PID控制器，实现污水处理过程设定值跟踪控制;溶解氧浓度通过第五分区氧气转换系数KLa5进行控制，硝态氮浓度通过内回流量Qa进行控制;按基准定义，氧气转换系数KLa5变化范围取为0

　　定义k时刻溶解氧浓度和硝态氮浓度的系统误差：

　　$\left\{\begin{array}{c}{e}_1\left(k\right)=x_1\left(k\right)-y_1\left(k\right)\\ e_2\left(k\right)=x_2\left(k\right)-y_2\left(k\right)\end{array}\right.---\left(11\right)$

　　式中y1(k)，y2(k)分别为k时刻溶解氧浓度和硝态氮浓度的实际测量值;

　　采用增量式PID控制，k时刻控制增量表达为

　　$\left\{\begin{array}{c}{\mathrm{\Delta u}}_1\left(k\right)=K_{P,1}\&CenterDot;\&lsqb;e_1\left(k\right)-e_1(k-1)\&rsqb;+K_{I,1}\&CenterDot;e_1\left(k\right)+K_{D,1}\&CenterDot;\&lsqb;e_1\left(k\right)-2\&CenterDot;e_1(k-1)+e_1(k-2)\&rsqb;\\ {\mathrm{\Delta u}}_2\left(k\right)=K_{P,2}\&CenterDot;\&lsqb;e_2\left(k\right)-e_2(k-1)\&rsqb;+K_{I,2}\&CenterDot;e_2\left(k\right)+K_{D,2}\&CenterDot;\&lsqb;e_2\left(k\right)-2\&CenterDot;e_2(k-1)+e_2(k-2)\&rsqb;\end{array}\right.---\left(12\right)$

　　k时刻系统控制量表达为

　　$\left\{\begin{array}{c}{u}_1\left(k\right)=u_1(k-1)+{\mathrm{\Delta u}}_1\left(k\right)\\ u_2\left(k\right)=u_2(k-1)+{\mathrm{\Delta u}}_2\left(k\right)\end{array}\right.---\left(13\right)$

　　式中KP,1，KI,1，KD,1分别为溶解氧浓度控制器的比例、积分和微分系数;KP,2，KI,2，KD,2分别为硝态氮浓度控制器的比例、积分和微分系数;根据实验实际运行结果，取定：KP,1=200，KI,1=15，KD,1=2，KP,2=20000，KI,2=5000，KD,2=400;

　　5.每隔2小时优化周期，进行多目标优化算法和智能决策计算，更新一次溶解氧浓度优化设定值x1(k)和硝态氮浓度优化设定值x2(k)，由底层PID执行跟踪控制任务;

　　本发明获得一种污水处理过程动态多目标优化控制方法，其创造性主要体现在：(1)将污水处理过程的能耗指标和水质指标当成本质多目标优化问题进行处理，考虑了目标间具有相互冲突的矛盾特性，采用NSGA-II算法同时优化能耗和水质指标，实现溶解氧浓度和硝酸氮浓度设定值优化;(2)由模糊神经网络建立污水处理过程带有约束的动态多目标优化模型，为污水处理过程能耗与水质的多目标优化问题提供了模型基础;以上两部分内容为本发明的保护范围。

　　本发明提出的污水处理过程动态多目标优化控制方法，利用神经网络建模方法在线建立污水处理过程的约束多目标优化模型，考虑能耗指标和水质指标的同时优化，在保证污水出水水质达标的情况下，可以实现污水处理过程能耗降低和出水水质指标优化。