Data-knowledge Driven Multiobjective Optimal Control for Municipal Wastewater Treatment Process
-
摘要: 城市污水处理过程优化控制是降低能耗的有效手段, 然而, 如何提高出水水质的同时降低能耗依然是当前城市污水处理过程面临的挑战. 围绕上述挑战, 文中提出了一种数据和知识驱动的多目标优化控制(Data-knowledge driven multiobjective optimal control, DK-MOC)方法. 首先, 建立了出水水质、能耗以及系统运行状态的表达关系, 获得了运行过程优化目标模型. 其次, 提出了一种基于知识迁徙学习的动态多目标粒子群优化算法, 实现了控制变量优化设定值的自适应求解. 最后, 将提出的DK-MOC应用于城市污水处理过程基准仿真模型1 (Benchmark simulation model No. 1, BSM1). 结果表明该方法能够实时获取控制变量的优化设定值, 提高了出水水质, 并且有效降低了运行能耗.
-
关键词:
- 城市污水处理过程 /
- 数据和知识驱动方法, 多目标优化控制 /
- 知识迁徙学习 /
- 动态多目标粒子群优化
Abstract: The optimal control is an effective method to reduce energy consumption for municipal wastewater treatment process. However, it is still a challenge to improve the effluent qualities and reduce energy consumption simultaneously for the municipal wastewater treatment process. To solve this problem, a data-knowledge driven multiobjective optimal control (DK-MOC) method is proposed in this paper. First, the expression relationship among effluent qualities, energy consumption and system operation state is established to obtain the operational optimal objective model. Second, a dynamic multiobjective particle swarm optimization algorithm, based on knowledge transfer learning method, is proposed to obtain the optimal set-points of control variables adaptively. Finally, the proposed DK-MOC method is applied to the benchmark simulation model No. 1 (BSM1) of the municipal wastewater treatment process. The results demonstrate that this proposed method can obtain the optimal set-points of control variables online, which not only improve the effluent qualities, but also reduce the operation energy consumption effectively. -
城市污水处理过程通过微生物吸附、分解以及氧化等过程, 将污水中的污染物去除, 是维持水资源可持续利用的重要手段[1-3]. 近年来, 随着城市污水排放量不断增加和国家排放标准的日趋严格, 城市污水处理厂规模不断增容, 导致城市污水处理过程运行能耗逐渐增大, 且出水水质不稳定[4-6]. 因此, 如何降低运行能耗和提高出水水质是城市污水处理过程亟待解决的问题.
为了降低城市污水处理过程运行能耗, Sun等提出了基于活性污泥2d模型的优化控制策略, 该优化控制策略建立了能耗与化学需氧量(Chemical oxygen demand, COD)、生化需氧量(Biochemical oxygen demand, BOD)、氨氮、总氮等状态变量的关系, 利用氨氮偏差计算溶解氧浓度的优化设定值[7]. 应用结果表明, 该优化控制策略能够有效降低能耗. De Gussem等提出了一种能耗实时优化控制策略, 该优化控制策略采用活性污泥能耗模型实现能耗与城市污水处理过程运行状态之间的关系描述, 并基于计算流体力学完成曝气池中叶轮开闭数量优化[8]. 应用效果表明, 该优化控制策略可以在不影响有机物去除效率的情况下降低能耗. 上述优化控制策略构建了基于机理分析的能耗目标模型描述城市污水处理过程运行能耗与运行状态变量的关系. 然而, 由于城市污水处理过程入水流量、入水水质等动态变化, 模型精度难以保证, 影响优化控制效果[9-11]. 因此, 基于数据驱动优化目标模型的优化控制方法获得了一些学者的关注, Zeng等提出了基于混合整数非线性规划模型的优化控制方法, 该优化控制方法完成了泵送能耗与泵送速度、操作时间的动态关系表达, 并采用贪婪类电磁机制算法获取最优的泵送速度和操作时间[12]. 仿真结果表明, 该优化控制方法可以显著地降低城市污水处理过程泵送能耗. Filipe等提出了能耗优化控制策略, 该优化控制策略采用梯度增强树建立了泵送能耗与运行状态变量之间的动态关系, 并基于深度强化学习实现控制变量优化设定值的求解[13]. 仿真结果表明, 该优化控制策略可以有效节约城市污水处理过程泵送能耗. 然而, 城市污水处理过程出水水质也是一个重要的运行指标, 如何降低城市污水处理过程能耗的同时提高出水水质是亟待解决的问题[14-16].
为了同时优化能耗和出水水质, Vega等提出了实时优化控制策略, 该优化控制策略采用实时优化算法实现溶解氧和硝态氮浓度优化设定值的求解[17]. 仿真结果表明, 该优化控制策略能够在降低能耗成本的同时提高出水水质. Zhang等提出了基于城市污水处理过程动态特征的优化控制方法, 该优化控制方法采用动态遗传算法求解水力停留时间、内循环比等控制变量的优化设定值, 并利用反向传播算法获取能耗和出水水质的权重系数[18]. 仿真结果表明, 该优化控制方法能够有效降低城市污水处理过程能耗并提高出水水质. 上述两种优化控制策略主要通过引入权重系数, 将多目标优化问题转化为单目标优化问题对控制变量设定值进行求解. 然而, 由于城市污水处理过程运行状态难以确定, 导致权重系数难以准确获取[19-20]. 因此, De Faria等提出了基于档案多目标进化算法的优化控制策略, 该优化控制策略能够通过Pareto支配关系平衡能耗成本、出水水质和生命周期影响, 获得合适的溶解氧浓度、外部碳源等控制变量优化设定值[21]. 仿真结果表明, 该优化控制策略有效降低了能耗成本, 提高了出水水质. Qiao等提出了基于自适应多目标差分进化算法的优化控制策略, 该优化控制策略能够获取合适的溶解氧和硝态氮浓度的优化设定值, 实现了城市污水处理过程能耗和出水水质同时优化[22]. 另外, 一些学者提出了基于多目标优化算法的优化控制方法同时优化了能耗和出水水质[23-25]. 然而, 城市污水处理过程优化目标模型是时变的, 如何获得有效的控制变量优化设定值依然是城市污水处理过程面临的挑战性问题[26].
城市污水处理过程优化控制方法主要包括以下三个过程: 首先, 建立运行过程目标模型, 描述运行目标与系统状态变量的动态特征; 其次, 设计优化方法, 获取控制变量优化设定值; 最后, 利用控制器跟踪设定值, 实现城市污水处理过程优化控制[27-29].文中提出了一种数据和知识驱动的多目标优化控制(Data-knowledge driven multiobjective optimal control, DK-MOC)方法, 主要贡献如下:
1) 设计基于自适应模糊神经网络(Adaptive fuzzy neural network, AFNN)的数据驱动优化目标模型, 描述能耗、出水水质以及相关状态变量的动态关系.
2) 提出基于知识迁徙学习的动态多目标粒子群优化算法, 将目标变化前具有较好分布性的解应用于目标变化后控制变量设定值的求解过程, 并检测种群进化状态调整种群进化方向, 提高种群的搜索性能, 获取有效的控制变量优化设定值.
1. 城市污水处理过程能耗分析
城市污水处理过程能耗主要包括泵送能耗和曝气能耗, 其中, 泵送能耗主要指内回流、外回流和污泥排放所产生的能耗, 与泵送能耗有关的主要状态变量有入水流量、硝态氮浓度、混合固体悬浮物浓度; 曝气能耗主要指控制生化反应池中的曝气过程产生的能耗, 与曝气能耗有关的主要状态变量有入水流量、溶解氧浓度、硝态氮浓度、氨氮浓度、混合固体悬浮物浓度. 为了实现对能耗的定量描述, 文中给出基于城市污水处理过程基准仿真模型1 (Benchmark simulation model No. 1, BSM1)的曝气能耗计算方式[30]:
$${\rm{AE(}}t{\rm{) = }}\frac{{{S_{O,sat}}}}{{1\;800\tau }}\int_t^{t + \tau } {\sum\limits_{i = 1}^{ 5} {{V_i} \times {K_L}a_i(t){\rm{d}}t} } $$ (1) 式中, AE(t)为t时刻曝气能耗,
$\tau $ 为计算间隔,${{{S_{O,sat}}}} $ 为溶解氧饱和浓度, Vi为第i个反应池的体积,${K_L}a_i(t)$ 为t时刻氧传递系数, 溶解氧浓度可通过调整氧传递系数进行调节. 泵送能耗的计算方式为$$\begin{split} {\rm{PE}}(t)=\;&\frac{1}{\tau }\int_t^{t + \tau } ( 0.004 \times {Q_a}(t) + 0.05 \times {Q_w}(t)+\\ & 0.008 \times {Q_r}(t) ) {\rm{d}}t \\[-10pt]\end{split}$$ (2) 式中, PE(t)为t时刻泵送能耗, Qa(t)为t时刻的内回流量, 硝态氮浓度可通过调整内回流量进行调节, Qr(t)为t时刻的外回流量, Qw(t)为t时刻的污泥排放量.
城市污水处理过程的能耗和出水水质是相互冲突的目标, 同时受多个变量共同影响. 为了提高城市污水处理过程出水水质, 城市污水处理厂通常采用增加曝气量, 延长曝气时间, 调节内外回流等方式调节生化反应过程. 然而, 曝气量和回流流量的增加会导致能耗提高. 目前, 高能耗已成为制约城市污水处理厂发展的瓶颈问题, 给国家能源和资源节约产生了严峻挑战. 为此, 城市污水处理多目标优化控制得到了广泛的研究, 成为提高出水水质, 降低能耗, 促进可持续发展的重要技术手段. 然而, 在城市污水处理过程中, 能耗和出水水质受入水水质、入水流量等多个变量的影响, 具有动态时变特性. 在动态时变环境下, 如何获取合适的控制变量设定值是城市污水处理过程优化控制面临的挑战性问题.
2. 数据和知识驱动的城市污水处理过程多目标优化控制
文中提出了DK-MOC方法, 主要包括以下部分: 首先, 建立数据驱动的优化目标模型, 实现能耗、出水水质和运行状态变量之间动态关系的准确描述. 其次, 设计基于知识迁徙学习的动态多目标粒子群优化算法求解控制变量优化设定值. 最后, 采用比例积分微分(Proportional integral differential, PID)控制器跟踪优化设定值. DK-MOC流程如图1所示.
2.1 运行过程目标模型构建
为了描述城市污水处理过程优化目标与系统运行状态变量之间的动态关系, 文中建立了基于AFNN的出水水质和能耗优化目标模型:
$$\min {{F}}(t) = {[{{EQ}}(t),{{EC}}(t)]^{\rm{T}}}$$ (3) 式中, F(t)表示运行目标函数, EQ(t)为出水水质模型, EC(t)为能耗模型.
$${{EQ}}(t) = {\boldsymbol{\varphi }}(t){{\boldsymbol{w}}^{1{\rm{T}}}}(t)$$ (4) $${{EC}}(t) = {\boldsymbol{\varphi }}(t){{\boldsymbol{w}}^{2{\rm{T}}}}(t)$$ (5) 式中,
${\boldsymbol{\varphi }}(t)=[{{\varphi }}_1(t),\cdots,{{\varphi }}_{10}(t) ]$ 为模糊神经网络归一化层输出向量,${{\boldsymbol{w}}^q}(t)=[{{{w_1}}^q}(t),\cdots,{{{w}_{10}}^q}(t)]$ 为模糊神经网路归一化层与第q个输出连接的权重向量, 当q = 1时, 输出为出水水质, 当q = 2时, 输出为能耗. 归一化层的输出为$${\varphi _l}(t) = {{{\phi _l}(t)} \bigg/ {\sum\limits_{j = 1}^{10} {{\phi _j}} }}(t)$$ (6) 式中,
$l=1,\cdots,10 $ 为归一化层神经元数量.${\phi _j} (t)$ 为模糊规则层第j个神经元的输出.$${\phi _j}(t) = \prod\limits_{i = 1}^4 {{{\rm{e}}^{ - \tfrac{{{{({s_i}(t) - {\mu _{ij}}(t))}^2}}}{{2{{({{\rm{\sigma}} _{ij}}(t))}^2}}}}}} $$ (7) 式中, si(t)为输入层第i个输入,
${\boldsymbol{\mu}}_j(t) =[\mu_{1j}(t),\cdots, $ $ \;\mu_{ij}(t),\cdots,\;\mu_{4j}(t)]$ 为模糊规则层的中心向量,${\boldsymbol{\sigma}}_j(t) = $ $ [\sigma_{1j}(t),\cdots, \;\sigma_{ij}(t),\cdots,\;\sigma_{4j}(t)]$ 为模糊规则层宽度向量, 输入向量通过主成分分析和指标可操作性确定, 表示为$${\boldsymbol{s}}(t) = [{S_O}(t),{S_{NO}}(t),{S_{NH}}(t),MLSS(t)]$$ (8) 式中,
${\boldsymbol{s}}(t) $ 为t时刻模糊神经网络的输入向量,${S_O}(t) $ 为t时刻溶解氧浓度,${S_{NO}}(t) $ 为t时刻硝态氮浓度,${S_{NH}}(t) $ 为t时刻氨氮浓度,$MLSS(t) $ 为t时刻混合固体悬浮物浓度.为了提高模糊神经网络的自适应能力, 其参数更新规则为
$${{{\boldsymbol{\Theta}} }}(t + 1) = {{{\boldsymbol{\Theta}} }}(t) + {({{{\boldsymbol{\Phi}} }}(t) + \chi (t) \times {\boldsymbol{I}})^{ - 1}} \times {{{\boldsymbol{\Lambda}} }}(t)$$ (9) 式中,
${\boldsymbol{\Theta }}(t)\;=\;[{\boldsymbol{w}}^1(t),\;\;{\boldsymbol{w}}^2(t),\;\;{\boldsymbol{\mu}}_1(t),\cdots,\;\;{\boldsymbol{\mu}}_{10}(t), $ $ \;{\boldsymbol{\sigma}}_1(t),\cdots,\; {\boldsymbol{\sigma}}_{10}(t)]^{\rm{T}}$ 为参数矩阵,${\boldsymbol{\Phi }}(t)$ 为拟海瑟矩阵,${\boldsymbol{\Lambda }}(t)$ 为梯度向量, I为单位矩阵,$ \chi(t) $ 为自适应参数, 更新方式为$$\chi (t) = \chi (t - 1) \cdot { \frac{{\delta ^{\min }}(t) + \chi (t - 1)}{{\delta ^{\max }}(t) + 1}}$$ (10) ${\delta ^{\min }}(t) $ 和${\delta ^{\max }}(t) $ 分别为拟海瑟矩阵${\boldsymbol{\Phi }}(t)$ 的最小和最大特征值, 自适应参数$ \chi(t) $ 满足$0 < \chi(t)<1$ . 拟海瑟矩阵${\boldsymbol{\Phi }}(t)$ 和梯度向量${\boldsymbol{\Lambda }}(t)$ 的计算方式为$${\boldsymbol{\Phi }}(t) = \sum\limits_{q = 1}^2 {{{{\boldsymbol{J}}}}_q^{\rm{T}}(t){{\boldsymbol{J}}_q}(t)} $$ (11) $${\boldsymbol{\Lambda }}(t) = \sum\limits_{q = 1}^2 {{\boldsymbol{J}}_q^{\rm{T}}(t){{{{{e}}}}_q}(t)} $$ (12) 式中,
${{{e}}}_{{q}}(t)$ 为输出层第q个神经元的误差,${{\boldsymbol{J}}_q}(t) $ 为雅克比矩阵.AFNN充分考虑了城市污水处理过程动态以及易受扰动特性, 描述了出水水质、能耗以及运行状态变量之间的动态关系, 基于模糊规则增加了模型的抗干扰能力. 另外, 模型的参数更新规则能够提高模型的自适应能力, 提高模型精度.
2.2 控制变量优化设定值求解
为了降低能耗, 提高出水水质, 文中提出了基于知识迁徙学习的动态多目标优化算法, 获取控制变量溶解氧和硝态氮浓度的优化设定值. 在优化过程中, 知识的表示形式为
$${\boldsymbol{A}}{\rm{(}}t{\rm{) = [}}{\boldsymbol{a}}_1^*{\rm{(}}t{\rm{ - 1),}}\;{\boldsymbol{a}}_2^*{\rm{(}}t{\rm{ - 1),}} \cdots {\rm{,}}\;{\boldsymbol{a}}_H^*{\rm{(}}t{\rm{ - 1)]}}$$ (13) 式中, A(t)为t时刻种群初始化过程获得的知识,
$ {\boldsymbol{a}}_h^*{\rm{(}}t{\rm{ - 1)}}$ 为在t – 1 时刻转移的作为知识的第h个粒子位置信息, 即控制变量溶解氧和硝态氮浓度的优化值,$h=1,\;2,\cdots,\;H ,$ H为所转移的知识量.在新优化周期开始, 种群中H个粒子采用A(t)进行初始化, A(t)包括t − 1代中Pareto前沿上的边缘粒子和具有较强收敛性的粒子Knee points的位置信息, 剩余粒子采用如下公式分别生成
$${{\boldsymbol{x}}_i}{\rm{(}}0{\rm{) = }}{{\boldsymbol{x}}_{{\rm{min}}}}{\rm{ + (}}{{\boldsymbol{x}}_{{\rm{max}}}}{\rm{ - }}{{\boldsymbol{x}}_{{\rm{min}}}}{\rm{)}}\circ{{\boldsymbol{r}}_{1 \times 2}}$$ (14) 式中,
${\boldsymbol{x}}_i(0)$ 为第i个粒子的初始位置向量,$i\in [H+ $ $ 1,\;N]$ , N为种群规模,${{\boldsymbol{x}}_{{\rm{max}}}}=[S_{O,\;\max}, \;S_{NO,\;\max}]$ 和${{\boldsymbol{x}}_{{\rm{min}}}}=[S_{O,\;\min},\;S_{NO,\;\min}]$ 分别为溶解氧和硝态氮浓度的上界和下界,${{\boldsymbol{r}}_{1 \times 2}} $ 为分布于[0, 1]的1×2维随机向量.为了提高多目标粒子群优化算法的搜索能力, 文中基于历史信息的种群进化状态检测机制, 设计自适应进化方向选择策略, 平衡种群的探索和开发能力. 其中, 种群进化状态通过种群收敛性和种群多样性检测, 种群收敛性表示为
$${P_C}{\rm{(}}k{\rm{) = }}\sum\limits_{n = 1}^N {\sum\limits_{u = k - {k_0}}^k {{{\rm{e}}^{ \tfrac{{{{\tilde C}_n}{\rm{(}}k{\rm{)}}}}{{k - u + 1}}}}} } $$ (15) 式中,
${P_C}(k)$ 为第k代种群收敛性指标, k0表示历史迭代次数,$ {\tilde C_n}(k)$ 为种群中第n个粒子在第k代的收敛性, 计算方式为$${\tilde C_n}{\rm{(}}k{\rm{) = }}{{{E{{Q}}}}_n}{\rm{(}}k{\rm{ - }}1{\rm{) - }}{{{E{{Q}}}}_n}{\rm{(}}k{\rm{)}} + {{E}}{{C}_n}{\rm{(}}k{\rm{ - }}1{{) - E}}{{{C}}_n}{\rm{(}}k{\rm{)}}$$ (16) ${{{E{{Q}}}}_n}(k)$ 和${{{E{{C}}}}_n}(k)$ 分别为第n个粒子在第k代的出水水质和能耗值. 种群多样性表示为$${P_D}{\rm{(}}k{\rm{) = }}\sum\limits_{n = 1}^N {\sum\limits_{u = k - {k_0}}^k {{{\rm{e}}^{ \tfrac{{{{\tilde D}_n}{\rm{(}}k{\rm{)}}}}{{k - u + 1}}}}} } $$ (17) 式中,
${P_D}(k)$ 为第k代种群多样性指标,$ {\tilde D_n}(k)$ 表示第n个粒子在第k代的多样性分布信息, 计算方式为$$\begin{split} {\tilde D_n}{\rm{(}}k{\rm{)}} =\;& \left( \sum\limits_{i = 1}^N {|({{E}}{{{Q}}_n}(k) - {{E}}{{{Q}}_i}(k)|}\right. \\ & \left.+ \sum\limits_{i = 1}^N {|({{E}}{{{C}}_n}(k) - {{E}}{{{C}}_i}(k)|} \right)\bigg/N\end{split}$$ (18) 根据种群的收敛性和多样性, 文中将进化状态划分为5个阶段, 并设计自适应进化方向选择策略平衡种群的探索和开发能力, 自适应进化方向选择策略如下所示:
阶段1. 当PC(k) − PC(k − 1) > 0且PD(k) − PD(k − 1)>0时, 种群处于进化阶段, 速度和位置的更新公式为
$$\begin{split} v_{i,d}^{}(k + 1) =\;& \omega {v_{i,d}}(k) + {c_1}{r_1}({p_{i,d}}(k) - {x_{i,d}}(k))+\\ &{c_2}{r_2}({g_d}{\rm{(}}k{\rm{)}} - {x_{i,d}}(k))\\[-10pt] \end{split}$$ (19) $$x_{i,d}^{}(k + 1) = {x_{i,d}}(k) + v_{i,d}^{}(k + 1)$$ (20) 式中,
${x_{i,d}}(k) $ 为第i个粒子在第k代的第d维位置,${v_{i,d}}(k) $ 为第i个粒子在第k代的第d维速度,${p_{i,d}}(k) $ 为第i个粒子在第k代个体最优解的第d个维度,${g_{d}}(k) $ 为种群在第k代全局最优解的第d个维度, w为惯性权重, r1和r2为区间[0, 1]的随机数, c1和c2分别为自我认知参数和社会认知参数.阶段2. 当PC(k) − PC(k − 1) < 0且PD(k) −PD(k − 1) > 0时, 种群处于探索阶段, 速度和位置的更新公式为
$$\begin{split} {v_{i,d}}(k{\rm{ + }}1) = \;&\omega {v_{i,d}}(k){\rm{ + }}{c_1}{r_1}({p_{i,d}}(k){\rm{ - }}{x_{i,d}}(k))+\\ &{c_2}{r_2}({g_d}{\rm{(}}k{\rm{) - }}{x_{i,d}}(k)){\rm{ + }}{c_{\rm{3}}}{r_3}{C_d}(k) \end{split}$$ (21) $$x_{i,d}^{}(k + 1) = {x_{i,d}}(k) + v_{i,d}^{}(k + 1)$$ (22) 式中, r3为分布于区间[0, 1]的随机数, c3为收敛性方向飞行参数, Cd(k)为第k代收敛性方向的第d个维度, 表示为
$${C_d}(k){\rm{ = }}\frac{1}{{{k_0}{\rm{ + }}1}}({x_{l(k),d}}(k){\rm{ - }}{x_{l(k),d}}(k - {k_0}))$$ (23) $$l(k)=\mathop {{\rm{argmax}}}\limits_{n \in [1,N]} \left(\sum\limits_{u = k - {k_0}}^k {{{{{\rm{e}}}}^{ \tfrac{{{{\tilde C}_n}(k)}}{{k - u + 1}}}}} \right)$$ (24) 式中, l(k)表示在第k次迭代中具有最大收敛性的粒子编号.
阶段3. 当PC(k) − PC(k − 1) > 0且PD(k) − PD(k − 1) < 0时, 种群处于开发阶段, 速度和位置的更新公式为
$$\begin{split} v_{i,d}^{}(k+1) =\;& \omega {v_{i,d}}(k)+{c_1}{r_1}({p_{i,d}}(k)-{x_{i,d}}(k))+\\ &{c_2}{r_2}({g_d}(k) - {x_{i,d}}(k))+{c_4}{r_4}{D_d}(k) \end{split}$$ (25) $$x_{i,d}^{}(k +1) = {x_{i,d}}(k) + v_{i,d}^{}(k +1)$$ (26) 式中, r4为分布于区间[0, 1]的随机数, c4为多样性方向飞行参数, Dd(k)表示第k代多样性方向的第d个维度, 计算方式为
$${D_d}(k){\rm{ = }}\frac{1}{{{k_0} + 1}}({x_{p(k),d}}(k) - {x_{p(k),d}}(k - {k_0}))$$ (27) $$p(k)=\mathop {{\rm{argmax}}}\limits_{n \in [1,N]} \left(\sum\limits_{u = k - {k_0}}^k {{{\rm{e}}^{ \tfrac{{{{\tilde D}_n}(k)}}{{k - u + 1}}}}} \right)$$ (28) 式中, p(k)表示第k次迭代中种群多样性最大的粒子编号.
阶段4. 当PC(k) − PC(k − 1) < 0且PD(k) − PD(k − 1) < 0时, 种群处于退化阶段, 速度和位置的更新公式为
$$\begin{split} v_{i,d}^{}(k+1)=\;&\omega {v_{i,d}}(k)+{c_1}{r_1}({p_{i,d}}(k)-{x_{i,d}}(k))+\\ & {c_2}{r_2}({g_d}{\rm{(}}k{\rm{) - }}{x_{i,d}}(k))+\\ & 0{\rm{.5}}\times({c_3}{r_3}{C_d}(k)+{c_4}{r_4}{D_d}(k))\\[-10pt] \end{split}$$ (29) $$x_{i,d}^{}(k{\rm{ + }}1) = {x_{i,d}}(k){\rm{ + }}v_{i,d}^{}(k{\rm{ + }}1)$$ (30) 阶段5. 当PC(k) − PC(k − 1) = 0或PD(k) − PD(k −1) = 0时, 种群处于停滞阶段, 速度和位置的更新公式为
$$\begin{split} v_{i,d}^{}(k+1) =\;& \omega {v_{i,d}}(k)+{c_1}{r_1}({p_{i,d}}(k)-{x_{i,d}}(k))+\\ &{c_2}{r_2}({g_d}(k) -{x_{i,d}}(k))\\[-10pt] \end{split}$$ (31) $$x_{i,d}^{}(k{\rm{ + }}1) = \left\{ {\begin{aligned} &{x_{d,{\rm{min}}}}{\rm{ + }}({x_{d,{\rm{max}}}}{\rm{ - }}{x_{d,{\rm{min}}}}) \times G(0,1),\\ &\qquad\qquad\qquad\;\,{{r_{\rm{3}}} \le {p_b}(k{\rm{ + }}1)}\\ & {x_{i,d}}(k),\;\;\;\;\qquad{{r_{\rm{3}}} > {p_b}(k{\rm{ + }}1)} \end{aligned}} \right.$$ (32) 式中, G(0, 1)是一个服从均匀分布的随机值,
${x_{d,{\rm{min}}}} $ 和${x_{d,{\rm{max}}}} $ 是第d维的最小值和最大值, pb(k)为突变率:$${p_b}(k) = 0.5 - 0.5\times\frac{k}{K}$$ (33) 式中, K为种群最大迭代次数.
基于Pareto支配关系建立的外部档案采用拥挤距离维护并选择全局最优粒子. 城市污水处理过程水质达标即可排放, 因此, 在第K代外部档案中选择溶解氧和硝态氮浓度优化设定值时, 增加决策者偏好信息, 控制变量优化设定值选择对应能耗值最低的非支配解, t时刻溶解氧和硝态氮的优化设定值为SO*(K)和SNO*(K).
基于知识迁徙学习的动态多目标粒子群优化算法采用t时刻的最优粒子对t + 1时刻种群的位置进行初始化, 并利用历史信息来检测种群进化状态, 为种群提供不同的进化方向, 提高种群的搜索能力, 求解溶解氧和硝态氮浓度的优化设定值. 其中, Tmax为终止时间. 控制变量优化设定值的获取流程如图2所示.
基于知识迁徙学习的动态多目标粒子群优化算法的计算复杂度分析如下: 假设粒子群规模为N, 计算种群目标值的时间复杂度为O(MN), 计算每个粒子收敛性的时间复杂度为O(N), 计算每个粒子的多样性指标的时间复杂度为O(MN2), 计算种群收敛性和多样性指标的时间复杂度均为O(N), 粒子速度和位置更新的计算复杂度为O(N), 外部档案更新的计算复杂度为O(MN2), 选择全局最优解的计算复杂度为O(N), 因此, 基于知识迁徙学习的动态多目标粒子群优化算法的计算复杂度为O(MN2). 在城市污水处理过程中, 优化控制的周期一般为2小时, 因此, 多目标粒子群优化算法的计算复杂度是可以满足实际应用需求的.
3. 实验结果及分析
文中从优化目标模型和优化控制性能两方面验证DK-MOC策略的有效性, 采用BSM1仿真平台作为公共验证平台, 该平台是国际水协会专家提出的基于生物脱氮策略的城市污水处理过程仿真平台, 主要用于城市污水处理过程控制策略的测试, 能够为不同优化控制策略提供相同的测试环境, 可以公平地比较不同策略的性能. 而且, 该平台可设置不同的实验环境, 文中设置了干旱和阴雨天气两种实验环境. 实验在Windows10、CPU 1.80 GHz、MATLAB 2018上运行. 为了保证评价结果的有效性, 所有的实验运行30次, 取平均值作为结果.
3.1 优化目标模型性能
AFNN与自适应核函数(Adaptive kernel function, AKF)[22, 25]、遗传算法−人工神经网络(Genetic algorithm-artificial neural network, GA-ANN)[23] 以及最小二乘支持向量机(Least squares support vector machine, LSSVM)[24] 进行比较, 所有模型的更新频率为15 min, 训练过程中的迭代次数设置为10次. 均方根误差(Root mean square error, RMSE)和预测精度(Prediction accuracy, PA)用于评价模型的性能
$${\rm{RMSE}}=\sqrt{\frac{\sum\limits_{i = 1 }^L({y_i-{\hat{y}}_i})} {L}}$$ (34) $${\rm{PA}} = \left[ {1 - \left| {\frac{{\sum\limits_{i = 1}^L {({y_i} - {{\hat y}_i})} }}{{\sum\limits_{i = 1}^L {{y_i}} }}} \right|} \right] \times 100\;{\text {\%}} $$ (35) 式中, yi为第i次测试的实际值,
$ {\hat y_i}$ 为第i次测试时模型输出值, L为测试总次数.在干旱和阴雨天气下, RMSE和PA如表1和表2所示, 可以看出, 在能耗和出水水质指标预测中, 与AKF、GA-ANN和LSSVM相比, AFNN在干旱天气和阴雨天气下具有较低的RMSE和较高PA, 说明AFNN能够更精确地构建能耗和出水水质模型. 另外, 图3 ~ 4给出了干旱天气下基于AFNN的能耗和出水水质建模结果. 由于每日出水水质和能耗值为每日采样多次的累加和, 因此, 图3和图4的出水水质和能耗在数值上相对图5和图6较小. 图3 ~ 4的建模结果表明, 所采用的AFNN可以准确描述能耗和出水水质以及运行状态变量之间的动态关系.
表 1 干旱天气下能耗和出水水质的测量精度Table 1 The measurement accuracy of energy consumption and effluent quality in dry weather模型 EC EQ RMSE PA RMSE PA AKF 0.0082 98.71 % 0.0079 97.06 % GA-ANN 0.0123 93.36 % 0.0124 95.52 % LSSVM 0.0115 96.21 % 0.0139 94.31 % AFNN 0.0079 98.99 % 0.0070 98.73 % 表 2 阴雨天气下能耗和出水水质的测量精度Table 2 The measurement accuracy of energy consumption and effluent quality in rainy weather模型 EC EQ RMSE PA RMSE PA AKF 0.0093 96.60 % 0.0833 97.37 % GA-ANN 0.0146 95.72 % 0.1002 94.32 % LSSVM 0.0138 95.74 % 0.0940 96.76 % AFNN 0.0089 98.88 % 0.0776 98.89 % 3.2 优化控制性能
DK-MOC主要包括AFNN、基于知识迁徙学习的动态多目标粒子群优化算法和PID控制器. 为了验证基于知识迁徙学习的动态多目标粒子群优化算法与PID控制器共同构成的优化控制器的有效性, 文中利用文献[23-25]中提出的优化方法与PID控制器共同构建优化控制器, 建模方法都选择具有较高精度的AFNN. 对比的优化方法主要包括非支配排序遗传算法II (NSGAII)[23]、基于聚类的多目标粒子群优化(Cluster-based multiobjective particle swarm optimization, CMOPSO)[24]和自适应多目标粒子群优化(Adaptive multiobjective particle swarm optimization, AMOPSO)[25]. 另外, 为了验证基于多目标优化方法和基于单目标优化方法的优化控制性能, 文中根据基于知识迁徙学习的动态多目标粒子群优化算法设计基于知识迁徙的动态单目标粒子群优化算法(Knowledge transfer-based dynamic single objective particle swarm optimization, KT-DSOPSO), 其中, 出水水质作为约束, 能耗作为优化目标. 同时, 在所有的优化控制策略中, PID控制器的参数设置: 溶解氧比例系数为20, 硝态氮比例系数为10000, 溶解氧积分系数为5, 硝态氮积分系数为3000, 溶解氧微分系数为1, 硝态氮微分系数为100. 在所有的优化算法中, 设置种群规模和外部档案最大尺寸为40, 种群迭代次数均为50. BSM1中设置的总采样时间为14天, 优化周期为2小时, 在动态多目标粒子群优化算法中, 惯性权重ω为0.7, c1、c2、c3和c4分别为2.5、2.5、0.5和0.5, k0为4, H为8.
文中采用能耗和出水水质两个指标对优化控制策略的性能进行验证, 能耗值低说明优化控制过程中产生较少的能耗. 根据国际污水排放检测标准, 出水水质采用出水中氨氮、总氮、悬浮物、COD、BOD的浓度计算, 浓度越低, 水质越好. 干旱天气下平均每天的出水水质和能耗表示如图5 ~ 6所示. 为了描述DK-MOC的优化控制性能, 图7 ~ 8给出了干旱天气下溶解氧和硝态氮浓度的优化设定值和跟踪值, 可以看出, 所提出的优化算法可以获得动态的控制变量优化设定值, 使用PID可以跟踪变化的优化设定值. 根据图9所示, 干旱天气下溶解氧和硝态氮的平均跟踪误差为SO < 0.5 mg·L−1, SNO < 0.2 mg·L−1, 结果表明, 在DK-MOC策略中, PID能较好地跟踪SO和SNO的优化设定值, 保证了跟踪控制的有效性.
为了精确描述DK-MOC性能, 表3给出了DK-MOC的出水水质的组分结果, SNH表示氨氮浓度, Ntot表示总氮浓度, TSS表示悬浮物浓度, 表4给出了不同优化控制策略在干旱和阴雨天气下的能耗和出水水质结果, 由表3 ~ 4可以看出, 在干旱天气下, DK-MOC的TSS、COD和BOD浓度均远小于排放标准. 虽然SNH、Ntot接近排放标准, 但所有的水质指标都达标, 且其总出水水质指标优于其他优化控制策略; 在阴雨天气下, DK-MOC的SNH浓度为3.89 mg·L−1, Ntot浓度为17.01 mg·L−1, 水质指标均达标, 且运行能耗具有最小值. 另外, DK-MOC与基于KT-DSOPSO的优化控制器对比, 可以看出, 所提出的基于知识迁徙学习的动态多目标粒子群优化方法可以权衡能耗和出水水质的关系. 结果表明, DK-MOC可以通过获取动态优化设定值, 有效降低能耗, 提高出水水质.
表 3 DK-MOC的出水水质结果Table 3 Effluent quality results obtained by DK-MOC水质指标 干旱天气 阴雨天气 排放标准 SNH (mg·L−1) 3.44 3.89 <4 Ntot (mg·L−1) 17.41 17.01 <18 TSS (mg·L−1) 12.57 13.51 <30 COD (mg·L−1) 47.75 46.53 <100 BOD (mg·L−1) 2.71 2.93 <10 表 4 不同优化控制策略的比较结果Table 4 Comparison results of different optimal control strategies优化控制策略 干旱天气 阴雨天气 EC
(kW·h)EQ
(kg poll.
units)EC
(kW·h)EQ
(kg poll.
units)AFNN+NSGAII+PID 4015 8867 4969 9049 AFNN+CMOPSO+PID 4770 8915 4869 9293 AFNN+AMOPSO+PID 3848 7695 3870 7 944 AFNN+KT-DSOPSO+PID 3727 9062 3892 9479 DK-MOC 3 641 7 179 3 821 8106 4. 结论
文中提出的DK-MOC方法主要包括以下优点:
1) 设计数据驱动的城市污水处理过程多目标优化模型, 实现运行状态变量与出水水质、能耗之间的动态关系的准确描述.
2) 提出了基于知识迁徙学习的动态多目标粒子群优化算法, 通过知识迁徙学习和自适应进化方向选择策略提高种群搜索性能, 从而获取有效的优化设定值.
BSM1平台的测试效果显示: 与其他建模方式相比, AFNN可以提高出水水质和能耗模型精度; 与一些优化控制策略相比, 所提出的DK-MOC方法可以获得有效的控制变量优化设定值, 提高了出水水质, 降低了能耗. 结果证明了DK-MOC方法的有效性.
-
表 1 干旱天气下能耗和出水水质的测量精度
Table 1 The measurement accuracy of energy consumption and effluent quality in dry weather
模型 EC EQ RMSE PA RMSE PA AKF 0.0082 98.71 % 0.0079 97.06 % GA-ANN 0.0123 93.36 % 0.0124 95.52 % LSSVM 0.0115 96.21 % 0.0139 94.31 % AFNN 0.0079 98.99 % 0.0070 98.73 % 
下载: 导出CSV
表 2 阴雨天气下能耗和出水水质的测量精度
Table 2 The measurement accuracy of energy consumption and effluent quality in rainy weather
模型 EC EQ RMSE PA RMSE PA AKF 0.0093 96.60 % 0.0833 97.37 % GA-ANN 0.0146 95.72 % 0.1002 94.32 % LSSVM 0.0138 95.74 % 0.0940 96.76 % AFNN 0.0089 98.88 % 0.0776 98.89 %
下载: 导出CSV
表 3 DK-MOC的出水水质结果
Table 3 Effluent quality results obtained by DK-MOC
水质指标 干旱天气 阴雨天气 排放标准 SNH (mg·L−1) 3.44 3.89 <4 Ntot (mg·L−1) 17.41 17.01 <18 TSS (mg·L−1) 12.57 13.51 <30 COD (mg·L−1) 47.75 46.53 <100 BOD (mg·L−1) 2.71 2.93 <10
下载: 导出CSV
表 4 不同优化控制策略的比较结果
Table 4 Comparison results of different optimal control strategies
优化控制策略 干旱天气 阴雨天气 EC
(kW·h)EQ
(kg poll.
units)EC
(kW·h)EQ
(kg poll.
units)AFNN+NSGAII+PID 4015 8867 4969 9049 AFNN+CMOPSO+PID 4770 8915 4869 9293 AFNN+AMOPSO+PID 3848 7695 3870 7 944 AFNN+KT-DSOPSO+PID 3727 9062 3892 9479 DK-MOC 3 641 7 179 3 821 8106
下载: 导出CSV
-
[1] Tang R, Wang Y, Yuan S, Wang W, Yue Z, Zhan X, et al. Organoarsenic feed additives in biological wastewater treatment processes: Removal, biotransformation, and associated impacts. Journal of Hazardous Materials, 2021, 406: 124789. doi: 10.1016/j.jhazmat.2020.124789 [2] Iratni A, Chang N. Advances in control technologies for wastewater treatment processes: Status, challenges, and perspectives. IEEE/CAA Journal of Automatica Sinica, 2019, 6(2): 337-363. doi: 10.1109/JAS.2019.1911372 [3] 权利敏, 杨翠丽, 乔俊飞. 数据驱动的溶解氧浓度在线自组织控制方法. 自动化学报, DOI: 10.16383/j.aas.c210041Quan Lin-Min, Yang Cui-Li, Qiao Jun-Fei. Data-driven online self-organizing control for dissolved oxygen concentration, Acta Automatica Sinica, DOI: 10.16383/j.aas.c210041 [4] Chistiakova T, Wigren T, Carlsson B. Combined L2 -stable feedback and feedforward aeration control in a wastewater treatment plant. IEEE Transactions on Control Systems Technology, 2020, 28(3): 1017-1024. doi: 10.1109/TCST.2019.2891410 [5] Ben-David E A, Habibi M, Haddad E, Hasanin M, Angel D L, Booth A M, et al. Microplastic distributions in a domestic wastewater treatment plant: Removal efficiency, seasonal variation and influence of sampling technique. Science of The Total Environment, 2021, 752: 141880. doi: 10.1016/j.scitotenv.2020.141880 [6] Borzooei S, Amerlinck Y, Panepinto D, Abolfathi S, Nopens I, Nopens G, et al. Energy optimization of a wastewater treatment plant based on energy audit data: small investment with high return. Environmental Science and Pollution Research, 2020, 27: 17972–17985. doi: 10.1007/s11356-020-08277-3 [7] Sun J Y, Liang P, Yan X X, Zuo K C, Xiao K, Xia J L, et al. Reducing aeration energy consumption in a large-scale membrane bioreactor: Process simulation and engineering application. Water Research, 2016, 93: 205−213. doi: 10.1016/j.watres.2016.02.026 [8] De Gussem K, Fenu A, Wambecq T, Weemaes M. Energy saving on wastewater treatment plants through improved online control: Case study wastewater treatment plant antwerp-south. Water Science &Technology, 2014, 69(5): 1074−1079. [9] Zonta Z J, Kocijan J, Flotats X, Vrečko D. Multi-criteria analyses of wastewater treatment bio-processes under an uncertainty and a multiplicity of steady states. Water Research, 2012, 46(18): 6121-6131. doi: 10.1016/j.watres.2012.08.035 [10] Godini K, Azarian G, Kimiaei A, Dragoi E N, Curteanu S. Modeling of a real industrial wastewater treatment plant based on aerated lagoon using a neuro-evolutive technique. Process Safety and Environmental Protection, 2021, 148: 114-124. doi: 10.1016/j.psep.2020.09.057 [11] Qiu Y, Li J, Huang X, Shi H. A feasible data-driven mining system to optimize wastewater treatment process design and operation. Water, 2018, 10(10): 1342. doi: 10.3390/w10101342 [12] Zeng Y, Zhang Z, Kusiak A, Tang F, Wei X. Optimizing wastewater pumping system with data-driven models and a greedy electromagnetism-like algorithm. Stochastic En- vironmental Research and Risk Assessment, 2016, 30: 1263−1275. doi: 10.1007/s00477-015-1115-4 [13] Filipe J, Bessa R J, Reis M, Alves R, Póvoa P. Data-driven predictive energy optimization in a wastewater pumping station. Applied Energy, 2019, 252: 113423. doi: 10.1016/j.apenergy.2019.113423 [14] Canete J F, Saz-Orozco P, Gabriel J G, Baratti R, Ruano A, Rivas-Blanco I. Control and soft sensing strategies for a wastewater treatment plant using a neuro-genetic approach. Computers & Chemical Engineering, 2021, 144: 107146. [15] Han H G, Zhang L, Liu H X, Qiao J F. Multiobjective design of fuzzy neural network controller for wastewater treatment process. Applied Soft Computing, 2018, 67: 467–478. doi: 10.1016/j.asoc.2018.03.020 [16] 栗三一, 乔俊飞, 李文静, 顾锞. 污水处理决策优化控制. 自动化学报, 2018, 44(12): 2198-2209.Li San-Yi, Qiao Jun-Fei, Li Wen-Jing, Gu Ke. Advanced decision and optimization control for wastewater treatment plants. Acta Automatica Sinica, 2018, 44(12): 2198-2209. [17] Vega P, Revollar S, Francisco M, Martín J M. Integration of set point optimization techniques into nonlinear MPC for improving the operation of WWTPs. Computers and Chemical Engineering, 2014, 68: 78-95. doi: 10.1016/j.compchemeng.2014.03.027 [18] Zhang R, Xie W M, Yu H Q, Li W W. Optimizing municipal wastewater treatment plants using an improved multi-objective optimization method. Bioresource Technology, 2014, 157: 161-165. doi: 10.1016/j.biortech.2014.01.103 [19] Mooselu M G, Nikoo M R, Latifi M, Sadegh M, Al-Wardy M, Al-Rawas G A. A multi-objective optimal allocation of treated wastewater in urban areas using leader-follower game. Journal of Cleaner Production, 2020, 267: 122189. doi: 10.1016/j.jclepro.2020.122189 [20] Tejaswini E S S, Panjwani S, Gara U B B, Ambati S R. Multi-objective optimization based controller design for improved wastewater treatment plant operation. Environmental Technology & Innovation, 2021, 23: 101591 [21] De Faria A B B, Ahmadi A, Tiruta-Barna L, Spérandio M. Feasibility of rigorous multi-objective optimization of wastewater management and treatment plants. Chemical Engineering Research and Design, 2016, 155(Part B): 394-406. [22] Qiao J F, Hou Y, Zhang L, Han H G. Adaptive fuzzy neural network control of wastewater treatment process with multiobjective operation. Neurocomputing, 2018, 275: 383-393. doi: 10.1016/j.neucom.2017.08.059 [23] Huang M, Han W, Wan J, Ma Y, Chen X. Multi-objective optimisation for design and operation of anaerobic digestion using GA-ANN and NSGA-II. Journal of Chemical Technology & Biotechnology, 2016, 91(1): 226-233. [24] 李永明, 史旭东, 熊伟丽. 基于工况识别的污水处理过程多目标优化控制. 化工学报, 2019, 70(11): 4325-4336.Li Yong-Ming, Shi Xu-Dong, Xiong Wei-Li. Condition recognition based intelligent multi-objective optimal control for wastewater treatment. CIESC Journal, 2019, 70(11): 4325-4336. [25] 韩红桂, 张璐, 卢薇, 乔俊飞. 城市污水处理过程动态多目标智能优化控制研究. 自动化学报, 2021, 47(3): 620-629.Han Hong-Gui, Zhang Lu, Lu Wei, Qiao Jun-Fei. Research on dynamic multiobjective intelligent optimal control for municipal wastewater treatment process. Acta Automatica Sinica, 2021, 47(3): 620-629. [26] 杨翠丽, 武战红, 韩红桂, 乔俊飞. 城市污水处理过程优化设定方法研究进展. 自动化学报, 2020, 46(10): 2092-2108.Yang Cui-Li, Wu Zhan-Hong, Han Hong-Gui, Qiao Jun-Fei. Perspectives on optimal setting methods for municipal wastewater treatment processes. Acta Automatica Sinica, 2020, 46(10): 2092-2108. [27] Newhart K B, Holloway R W, Hering A S, Cath T Y. Data-driven performance analyses of wastewater treatment plants: A review. Water research, 2019, 157: 498-513. doi: 10.1016/j.watres.2019.03.030 [28] Han H G. Liu Z, Lu W, Hou Y, Qiao J F. Dynamic MOPSO-based optimal control for wastewater treatment process. IEEE Transactions on Cybernetics, 2021, 51(5): 2518−2528. [29] Christine S, Fu G, Butler D. Multi-objective optimization of wastewater treatment plant control to reduce greenhouse gas emissions. Water Research, 2014, 55: 52-62. doi: 10.1016/j.watres.2014.02.018 [30] Jeppsson U, Pons M N. The COST benchmark simulation model—current state and future perspective. Control Engineering Practice, 2004, 12(3): 299-304. doi: 10.1016/j.conengprac.2003.07.001 期刊类型引用(12)
1. 赵霞，石诗义，毛杨，武桢寓. 城市污水处理厂的“水-能-碳”协调调度模型. 电网技术. 2024(05): 1918-1928 . 
百度学术2. 韩红桂，王玉爽，刘峥，孙浩源，乔俊飞. 知识和数据驱动的污水处理反硝化脱氮过程协同优化控制. 自动化学报. 2024(06): 1221-1233 . 本站查看3. 王勇，苏晓亮，刘佳，焦峰，季子禾. 基于AI的污水处理监控与自动调节技术. 中国新技术新产品. 2024(13): 115-117+127 . 百度学术4. 韩红桂，白星，侯莹. 基于领域自适应的城市污水处理运行过程多工况优化控制. 中国科学:技术科学. 2024(09): 1652-1664 . 百度学术5. 郭鑫，李文静，乔俊飞. 基于自组织模块化神经网络的污水处理过程出水参数预测. 化工学报. 2024(09): 3242-3254 . 百度学术6. 姚清圆. 电气自动化技术在污水处理厂节能中的应用及发展趋势. 中国资源综合利用. 2024(12): 262-264 . 百度学术7. 杨静，王晓，王雨桐，刘忠民，李小双，王飞跃. 平行智能与CPSS:三十年发展的回顾与展望. 自动化学报. 2023(03): 614-634 . 本站查看8. 韩红桂，秦晨辉，孙浩源，乔俊飞. 城市污水处理过程自适应滑模控制. 自动化学报. 2023(05): 1010-1018 . 本站查看9. 魏华. 污水处理中泵站的优化控制研究. 中国设备工程. 2023(18): 120-122 . 百度学术10. 王鼎，赵明明，哈明鸣，乔俊飞. 基于折扣广义值迭代的智能最优跟踪及应用验证. 自动化学报. 2022(01): 182-193 . 本站查看11. 刘莹，王司阳. 基于改进遗传算法的污水处理优化控制方法. 清洗世界. 2022(10): 36-38 . 百度学术12. 张永合. 城市环境污水处理过程节能优化控制方法探索. 科技创新与应用. 2022(32): 142-145 . 百度学术其他类型引用(19)
-
下载:
计量
- 文章访问数: 2464
- HTML全文浏览量: 665
- PDF下载量: 866
- 被引次数: 31
