多梯度下降污水处理过程优化控制方法

发布时间：2018-6-29 16:52:04  中国污水处理工程网

　　申请日2014.11.02

　　公开(公告)日2015.02.18

　　IPC分类号G06F19/00; G05B13/04

　　摘要

　　针对污水处理过程高度非线性、强耦合性和不确定性严重等特点，本发明提出一种基于多梯度下降的污水处理过程优化控制方法，实现对污水处理过程中溶解氧DO和硝态氮SNO浓度的同时控制;该控制方法通过构建控制系统的多目标函数，基于多梯度下降的优化方法求解污水处理过程的多目标问题，通过控制优化后的曝气量和内循环回流量达到控制溶解氧DO和硝态氮SNO浓度的目的;解决了污水处理过程多目标优化控制的问题，提高了污水处理过程中溶解氧DO和硝态氮SNO的控制精度，保障污水处理过程正常运行，促进污水处理厂高效稳定运行。

　　权利要求书

　　1.一种基于多梯度下降的污水处理过程优化控制方法，

　　针对序批式间歇活性污泥系统中溶解氧DO和硝态氮SNO浓度进 行控制，以曝气量和内循环回流量为控制量，溶解氧DO和硝态氮 SNO浓度为被控量;

　　其特征在于，包括以下步骤：

　　(1)设计用于污水处理过程中溶解氧DO和硝态氮SNO浓度预测 控制方法的多目标函数：

　　$J_1(U,X)=\frac{1}{2}{U}^T{W}_{1}U+\frac{1}{2}{X}^T{P}_{1}X;$

　　$J_2(U,X)=\frac{1}{2}{U}^T{W}_{2}U+\frac{1}{2}{X}^T{P}_{2}X;---\left(1\right)$

　　其中W1和W2是目标函数J1和J2的Hessian矩阵，P1和P2是目标函 数J1和J2的权值矩阵，T为矩阵的转置，溶解氧DO和硝态氮SNO的状态量X和控制量U为：

　　U(t)=[u(t)，u(t+1)，L，u(t+H-1)]T;

　　u(t)=[u1(t)，u2(t)];

　　X(t)=[x(t)，x(t+1)，L，x(t+H-1)]T;

　　x(t)=[x1(t)，x2(t)]; (2)

　　u1是曝气量，u2是内循环回流量，x1是溶解氧DO浓度值，x2是硝态 氮SNO浓度值，H为变化时域，H∈[1，5];限制条件：

　　Δu1(t)=u1(t)-u1(t-1);

　　Δu2(t)=u2(t)-u2(t-1);

　　|Δu1(t)|≤Δu1，max;|Δu2(t)|≤Δu2，max;

　　u1，min≤u1(t)≤u1，max;u2，min≤u2(t)≤u2，max;

　　x1，min≤x1(t)≤x1，max;x2，min≤x2(t)≤x2，max; (3)

　　其中，Δu1，max是控制器允许的最大曝气调整量，u1，min是控制器允许的 最小曝气量，u1，max是控制器允许的最大曝气量，Δu2，max是控制器允许 的最大内循环回流调整量，u2，min是控制器允许的最小内循环回流量， u2，max是控制器允许的最大内循环回流量，x1，min是控制对象允许的最 小溶解氧浓度，x1，max是控制对象允许的最大溶解氧浓度，x2，min是控 制对象允许的最小硝态氮浓度，x2，max是控制对象允许的最大硝态氮 浓度，Δu1，max、u1，min、u1，max、Δu2，max、u2，min、u2，max、x1，min、x1，max、x2，min和x2，max根据控制系统设备设置;

　　(2)设计用于污水处理过程中溶解氧DO和硝态氮SNO浓度预测 控制的多目标优化方法，具体为：

　　①根据状态量X(t)和上一时刻的控制量U(t-1)，计算每个目标函 数的下降梯度：

　　$\begin{array}{c}{\bar{u}}_i\left(t\right)=-\frac{{▿j}_i(U(t-1),X\left(t\right))}{\left|\right|{▿j}_i(U(t-1),X\left(t\right))\left|\right|}\\ =-\frac{W_{i}U(t-1)+{P_i}^{T}X\left(t\right)}{\left|\right|W_{i}U(t-1)+{P_i}^{T}X\left(t\right)\left|\right|},i=\mathrm{1,2}\end{array};---\left(5\right)$

　　②根据每个目标函数的下降梯度计算多个目标函数的共同下降 梯度：

　　$▿J(U(t-1),X\left(t\right))=\underset{i=1}{\overset{2}{\Sigma }}{\beta }_i\left(t\right){\bar{u}}_i\left(t\right);---\left(6\right)$

　　其中，目标函数间的系数βi(t)表示为：

　　${\beta }_i\left(t\right)=\frac{{\left|\right|{\bar{u}}_i\left(t\right)\left|\right|}^2}{{\left|\right|{\bar{U}}^T\left(t\right)\bar{U}\left(t\right)\left|\right|}^2},i=\mathrm{1,2};---\left(7\right)$

　　并且$\bar{U}\left(t\right)={[{\bar{u}}_1^T\left(t\right),{\bar{u}}_2^T\left(t\right)]}^T;$

　　③判断共同下降梯度的大小，如果转到步骤④; 如果转到步骤⑤;σ是设定的共同下降梯度阀值，σ∈(0， 0.01];

　　④计算当前时刻的控制量U(t)

　　$U\left(t\right)=U(t-1)+h▿J(U(t-1),x\left(t\right));---\left(8\right)$

　　其中，h为梯度下降步长，h∈(0，0.1];转到步骤⑥;

　　⑤计算当前时刻的控制量U(t)

　　U(t)=U(t-1); (9)

　　⑥抽取当前时刻控制量U(t)中的u(t)，u(t)为当前时刻污水处理 系统预测控制系统的实际输入量;

　　(3)利用求解出的u(t)对溶解氧DO和硝态氮SNO进行控制，u1(t) 即为曝气量，u2(t)即为内循环回流量，整个控制系统的输出为实际溶 解氧DO和硝态氮SNO浓度值。

　　说明书

　　一种基于多梯度下降的污水处理过程优化控制方法

　　技术领域

　　本发明利用基于多梯度下降的模型预测控制方法实现污水处理过程中 溶解氧DO和硝态氮SNO浓度的控制，溶解氧DO和硝态氮SNO的浓度直接决 定了污水处理的效果，对出水水质和污水处理过程的能耗有着重要影响。 将基于多梯度下降的模型预测控制方法应用于污水处理系统，对溶解氧DO 和硝态氮SNO进行优化控制，既可节约投资和运行成本，又能促使污水处理 厂高效稳定运行，是先进制造技术领域的重要分支，既属于控制领域，又 属于水处理领域。

　　背景技术

　　污水生化处理过程的实质是利用人工培养的微生物群体来吸附、分解、 氧化污水中可生物降解的有机物，通过生物化学反应，将这些有机物从污 水中分离出来，使污水得到净化。城市污水处理厂的有效运行对于削减污 染物排放量、改善水环境、促进水资源有效利用发挥了重要作用，是缓解 水资源匮乏的重要举措;国家中长期科技发展规划中提出要研究并推广高 效、低能耗的污水处理新技术。因此，本发明的研究成果具有广阔的应用 前景。

　　污水处理过程中硝化反应主要是在有氧条件下发生的，溶解氧DO浓度 的大小直接影响了硝化反应进程，溶解氧DO浓度变大时，系统中出水氨氮 和总氮的浓度就会呈下降趋势，但是当溶解氧DO浓度达到一定值时，出水 中的氨氮的变化幅度就减弱了，而且总氮也受硝态氮的影响，硝态氮增加 时，总氮浓度也会升高。然而，污水处理过程中反硝化反应主要是在缺氧 环境下进行的，缺氧区的硝态氮SNO浓度是衡量脱氮效果的重要指标，它反 映了反硝化反应过程的进程，将硝态氮SNO浓度控制在一个合适的范围内， 能够提高反硝化反应的潜力。因此，对曝气池中溶解氧DO和硝态氮SNO的 控制非常重要，需要将溶解氧DO和硝态氮SNO浓度控制在一定范围内，才 能高效利用硝化反应的潜力。传统的开关控制或者PID控制，虽然是当前应 用较为广泛的控制方法，但是由于氧气的溶解过程受入水水质、温度和pH 值等方面的影响，具有高度非线性、强耦合性、时变、大滞后和不确定性 等特点。采用传统的开关控制或者PID控制方法自适应能力较差，往往不能 取得理想的控制效果。因此，必须寻求新的控制方法，以满足污水处理过 程控制的需要。优良的控制可以节省污水处理运行费用，同时也是减少和 应对异常工况发生、保障污水处理过程正常运行的关键。此外，通过提高 污水处理过程自动化水平，还可以有效地减少运行管理和操作人员，降低 运行费用。

　　本发明设计了一种基于多梯度下降的模型预测控制方法，主要通过多 梯度下降方法对控制目标进行优化，实现溶解氧DO和硝态氮SNO的在线控 制。

　　发明内容

　　本发明获得了一种基于多梯度下降的污水处理模型预测控制方法，该 控制器基于多梯度下降的优化方法求解污水处理过程的多目标问题，通过 控制优化后的曝气量和内循环回流量达到控制溶解氧DO和硝态氮SNO浓度 的目的;解决了污水处理过程多目标优化控制的问题，提高了污水处理过 程中溶解氧DO和硝态氮SNO浓度的控制精度，保障污水处理过程正常运行;

　　本发明采用了如下的技术方案及实现步骤：

　　1.一种基于多梯度下降的污水处理过程优化控制方法包括以下步骤：

　　针对序批式间歇活性污泥系统中溶解氧DO和硝态氮SNO进行控制，以曝 气量和内循环回流量为控制量，溶解氧DO和硝态氮SNO浓度为被控量;

　　(1)设计用于污水处理过程中溶解氧DO和硝态氮SNO浓度预测控制方 法的多目标函数：

　　$\begin{array}{c}{J}_1(U,X)=\frac{1}{2}{U}^T{W}_{1}U+\frac{1}{2}{X}^T{P}_{1}X;\\ J_2(U,X)=\frac{1}{2}{U}^T{W}_{2}U+\frac{1}{2}{X}^T{P}_{2}X;\end{array}---\left(1\right)$

　　其中W1和W2是目标函数J1和J2的Hessian矩阵，P1和P2是目标函数J1和J2的权值矩阵，T为矩阵的转置，溶解氧DO和硝态氮SNO的状态量X和控制量 U为：

　　U(t)=[u(t),u(t+1),…,u(t+H-1)]T;

　　u(t)=[u1(t),u2(t)];

　　X(t)=[x(t),x(t+1),…,x(t+H-1)]T;

　　x(t)=[x1(t),x2(t)]; (2)

　　u1是曝气量，u2是内循环回流量，x1是溶解氧DO浓度值，x2是硝态氮SNO浓度值，H为变化时域，H∈[1,5];限制条件：

　　△u1(t)=u1(t)-u1(t-1);

　　△u2(t)=u2(t)-u2(t-1);

　　|△u1(t)|≤△u1,max;|△u2(t)|≤△u2,max;

　　u1,min≤u1(t)≤u1,max;u2,min≤u2(t)≤u2,max;

　　x1,min≤x1(t)≤x1,max;x2,min≤x2(t)≤x2,max; (3)

　　其中，Δu1,max是控制器允许的最大曝气调整量，u1,min是控制器允许的最小 曝气量，u1,max是控制器允许的最大曝气量，Δu2,max是控制器允许的最大内 循环回流调整量，u2,min是控制器允许的最小内循环回流量，u2,max是控制器 允许的最大内循环回流量，x1,min是控制对象允许的最小溶解氧浓度，x1,max是控制对象允许的最大溶解氧浓度，x2,min是控制对象允许的最小硝态氮浓 度，x2,max是控制对象允许的最大硝态氮浓度，Δu1,max、u1,min、u1,max、Δu2,max、 u2,min、u2,max、x1,min、x1,max、x2,min和x2,max根据控制系统设备设置;

　　(2)设计用于污水处理过程中溶解氧DO和硝态氮SNO浓度预测控制的 多目标优化方法，具体为：

　　①根据状态量X(t)和上一时刻的控制量U(t-1)，计算每个目标函数的 下降梯度：

　　${\bar{u}}_i\left(t\right)=-\frac{▿J_i(U(t-1),X\left(t\right))}{\left|\right|▿J_i(U(t-1),X\left(t\right))\left|\right|}=-\frac{W_{i}U(t-1)+P_i^{T}X\left(t\right)}{\left|\right|W_{i}U(t-1)+P_i^{T}X\left(t\right)\left|\right|},i=\mathrm{1,2};---\left(5\right)$

　　②根据每个目标函数的下降梯度计算多个目标函数的共同下降梯度：

　　$▿J(U(t-1),X\left(t\right))=\underset{i=1}{\overset{2}{\Sigma }}{\beta }_i\left(t\right){\bar{u}}_i\left(t\right);---\left(6\right)$

　　其中，目标函数间的系数βi(t)表示为：

　　${\beta }_i\left(t\right)=\frac{{\left|\right|{\bar{u}}_i\left(t\right)\left|\right|}^2}{{\left|\right|{\bar{U}}^T\left(t\right)\bar{U}\left(t\right)\left|\right|}^2},i=\mathrm{1,2};---\left(7\right)$

　　并且$\bar{U}\left(t\right)={[{\bar{u}}_1^T\left(t\right),{\bar{u}}_2^T\left(t\right)]}^T;$

　　③判断共同下降梯度的大小，如果转到步骤④;如果 转到步骤⑤;σ是设定的共同下降梯度阀值，σ∈(0,0.01];

　　④计算当前时刻的控制量U(t)

　　$U\left(t\right)=U(t-1)+h▿J(U(t-1),x\left(t\right));---\left(8\right)$

　　其中，h为梯度下降步长，h∈(0,0.1];转到步骤⑥;

　　⑤计算当前时刻的控制量U(t)

　　U(t)=U(t-1); (9)

　　⑥抽取当前时刻控制量U(t)中的u(t)，u(t)为当前时刻污水处理系统预 测控制系统的实际输入量;

　　(3)利用求解出的u(t)对溶解氧和硝态氮进行控制，u1(t)即为曝气量， u2(t)即为内循环回流量，整个控制系统的输出为实际溶解氧DO和硝态氮 SNO浓度值;

　　本发明的创造性主要体现在：

　　(1)本发明针对当前活性污泥法污水处理过程是一个复杂的、动态的生 物反应过程，需要溶解氧DO和硝态氮SNO浓度同时达标，然而污水处理过 程不仅具有非线性、大时变等特点，而且溶解氧DO和硝态氮SNO浓度之间 存在强耦合关系，难以实现溶解氧DO和硝态氮SNO浓度的多目标控制;根 据多目标优化的特点，采用了基于多梯度下降的模型预测控制方法实现溶 解氧DO和硝态氮SNO浓度的控制，具有控制精度高、稳定性好等特点;

　　(2)本发明采用基于多梯度下降的模型预测控制方法对污水处理过程溶 解氧DO和硝态氮SNO浓度进行控制，该优化控制方法解决了多个目标函数 的优化求解问题，实现了溶解氧DO和硝态氮SNO浓度实时闭环精确控制， 避免了当前污水处理厂需要设计多个控制器进行控制的复杂过程，具有结 构简单、实时性好等特点;

　　特别要注意：本发明只是为了描述方便，采用的是对溶解氧DO和硝态 氮SNO浓度的控制，同样该发明也可适用污水处理过程氨氮的控制等，只要 采用了本发明的原理进行控制都应该属于本发明的范围。