污水处理过程动态多目标优化控制方法
【技术领域】
[0001] 本发明提出一种基于神经网络和非占优精英遗传算法NSGA-II的污水处理过程 动态多目标优化控制方法，对污水处理过程能耗与水质两个重要且相互冲突的性能指标同 时优化，实现好氧区溶解氧浓度和厌氧区硝态氮浓度设定值的动态寻优。污水处理优化控 制是水处理领域重要的研究课题，属于水研究领域，也属于智能控制领域。
【背景技术】
[0002] 污水处理过程是高能耗的复杂流程工业系统，其操作连续运行并需要保证出水水 质达标。随着国家对污水处理过程节能降耗的要求和人类对环境保护意识的增强，污水处 理过程的优化控制引起了广泛关注。污水处理过程优化控制是提高污水处理过程性能、实 现节能降耗的有效途径。因此，本发明的研究成果具有广阔的应用前景。
[0003] 污水处理优化控制的目标是在满足出水水质达标的前提下，获得能耗和成本的降 低。在污水处理过程优化控制研究初期，单变量优化取得了较为丰硕的成果，并以溶解氧浓 度优化为主。为了获得更好的优化性能，污水处理过程多变量优化控制近年来引起了众多 学者关注。前置反硝化活性污泥法处理的污水过程，第五分区的溶解氧浓度和第二分区的 硝态氮浓度强烈影响污水处理过程的硝化和反硝化水平，直接决定出水水质，同时也是曝 气能耗AE和栗送能耗PE高低的重要影响参数。因此，对溶解氧浓度S。和硝态氮浓度SN。的 设定值进行优化已经成为提高污水处理过程优化性能的重要手段。
[0004] 由于活性污泥处理过程中，存在复杂的物理、生化反应，且各种反应过程相互影 响，使得整个系统呈现较强的耦合特性，甚至表现为性能指标的相互冲突，如，能耗与出水 水质指标间便具有这种冲突特性。因此，将污水处理优化控制视为本质多目标优化问题更 符合污水处理过程优化问题的实质，但这方面的研究工作才刚起步。国内外现有污水处理 多目标优化控制方法中，存在的主要问题有：考虑了污水过程具有多目标优化特性，但优化 算法依然将其处理成单目标优化问题；以优化处理厂设计参数为主的多目标优化方法，采 用离线优化模式，不适于在线优化控制；在模型预测控制框架下考虑控制性能指标的冲突 特性，但未考虑水质、能耗等经济指标优化。此外，污水处理过程能耗、水质的多目标优化控 制问题，难点还在于能耗-水质-优化设定值模型较难建立。
【发明内容】
[0005] 本发明获得一种基于神经网络和非占优精英遗传算法NSGA-II的污水处理过程 多目标优化控制方法，实现污水处理过程溶解氧浓度和硝态氮浓度设定值优化，保证出水 水质达标前提下，降低能耗和处理成本；本发明采用如下技术方案及实现步骤，其特征在 于，包括以下步骤：
[0006] 1.以出水水质达标为约束条件，同时优化能耗和出水水质两个性能指标，构造污 水处理过程的约束多目标优化模型；
[0007] 1. 1多目标优化模型构造
[0008] 建立优化设定值与性能指标间的函数关系，其实质是建立溶解氧浓度和硝态氮浓 度的优化设定值与能耗指标、出水水质指标间的代理模型，实现对性能指标的预测和评价； 令xjk)为溶解氧浓度设定值，x2(k)为硝态氮浓度设定值，x(k) = [XiGOAGO]为设定 值组成的优化向量，k为时刻，为方便表述，记x(k)为x;fAE(x)为优化变量与曝气能耗间的 函数表达，fPE(x)为优化变量与栗送能耗间的函数表达，能耗指标定义为曝气能耗与栗送能 耗之和，即能耗指标fEC；(X) =fAE(X)+fpE(X)，fEQ(X)为优化变量与出水水质指标间的函数表 达；
[0009] 构建如下多目标优化模型，
[0012] 式中，xi,蚌分别为溶解氧浓度优化设定值的下限和上限值，4 <分别硝态氮浓 度优化设定值的下限和上限值；gl(x)为出水氨氮浓度SNH与优化设定值间的函数关系，gl(X)-4<0表示出水氨氮需满足的约束不等式，4(mg/L)为其约束上限值，g2(x)为出水总 氮浓度Ntot与优化设定值间的函数关系，g2(x)_18 < 0表示出水总氮需满足的约束不等式， 18(mg/L)为其约束上限值；本发明选择出水总氮浓度和出水氨氮浓度做为水质约束条件， 原因在于这两个出水浓度指标最为重要，且易出现超标现象；考虑其它出水浓度参数作为 约束时，如，固体悬浮物浓度SS低于其约束上限值30 (mg/L)，化学需氧量COD低于其约束上 限值100 (mg/L)，生物需氧量B0D5低于其约束上限值10 (mg/L)，本发明方法仍然适用，处理 方法如同出水总氮和出水氨氮，即建立相应的优化设定值与各出水浓度参数的函数关系， 并与约束上限值一起构成多目标优化模型中的约束条件；各出水浓度参数的约束上限值基 于国际水质协会制定的基准给出；本发明中，溶解氧浓度为污水处理第五分区溶解氧浓度， 硝态氮浓度为污水处理第二分区硝态氮浓度；
[0013] 1. 2模糊神经网络在线建模
[0014] 优化设定值与曝气能耗、栗送能耗、出水水质的函数关系，采用TS模糊神经网络 在线建模方法得到，模糊神经网络结构如图1所示；神经网络的输入向量由Z(k)表示，对于 曝气能耗模型、栗送能耗模型，网络的输入为k时刻溶解氧浓度和k时刻硝态氮浓度；对于 出水氨氮、出水总氮、水质性能指标的函数模型，网络输入为k时刻溶解氧浓度、k时刻硝态 氮浓度，以及k时刻入水流量；为方便表述，z(k)简记为z;
[0015] 模糊神经网络的输出表示为
[0016] >·(Α〇 =Φ(ζ)·Η(ζ,θ) (3)
[0017] 其中，Φ(ζ) =f(A(z)) (4)
[0018]
[0019]H(z，θ) = [1，ζτ] ·θ(5)式中，Α(ζ)为网络输入经模糊化后得到的隶属度矩 阵，f为每条规则隶属度乘积及规一化运算，隶属度函数取高斯函数；Θ为Θ(k)的简记，表 示k时刻神经网络的后件参数矩阵；上标T表示向量的转置；.?·(幻为k时刻网络输出；
[0020] 设k时刻网络学习的目标函数为：
[0022] 采用梯度下降算法更新网络参数，调整公式为：
[0024] 式中，a(k)为k时刻网络的学习参数向量，包括网络的后件参数、隶属度函数的 中心和宽度；α(k+Ι)为k+Ι时刻网络学习参数向量；η为网络的学习率，一般在（〇, 1)之 间取值，较大的学习率能提高学习速度，但不利于网络收敛，较小的学习率利于网络收敛， 但学习速度较慢，基于实验运行情况，取定学习率η为〇. 1 ;y(k)为k时刻系统实际物理 量输出，基于污水处理过程基准仿真模型BSM1产生，BSM1由国际水质协会与欧盟科学技术 与合作组联合提出，是当前国际上公认的污水处理过程控制性能评价的基准，其布局如图2 所示；
[0025] 1. 3模型中约束条件处理
[0026] 采用惩罚函数法将建立的约束优化问题转化成无约束优化问题；定义惩罚项
[0027]fpenalty(X)=max{g!(X)-4, 0}+max{g2(X)-18, 0} (8)
[0028] 加入惩ill项的能耗指标与水质性能指标表达为
[0030] 其中，C为惩罚因子，一般取较大正实数，其值越大对不满足约束条件的解的惩罚 力度越大，但易产生不可行域附近可行解丢失现象，根据实验运行结果，取定惩罚因子为 100000 ；
[0031] 2.多目标优化算法
[0032] 对建立的污水处理过程多目标优化模型采用NSGA-II进化算法进行求解，在溶 解氧浓度和硝态氮浓度的优化设定值范围内，获得多目标优化问题的一组等同优秀的解， 即Pareto最优解；NSGA-II进化算法中待优化性能指标为（9)式带有惩罚项的能耗指标 fEe(x)和水质性能指标fEQ(x);根据本发明中的实际优化问题及实验运行结果，取定 搜索空间D为2维，种群规模N为40,最大进化代数Μ为30;按照标准NSGA-II进化算法 运行程序，直至达到最大进化代数M，算法结束，给出污水处理过程多目标优化问题的一组 Pareto最优解{Xp，1彡p彡N}，即一组等同优秀的溶解氧浓度和硝态氮浓度的优化设定 值；
[0033] 3.智能决策
[0034] 为了实现污水处理过程的闭环控制，在多目标优化算法获得的一组Pareto最优 解中，确