曝气控制方案综述

1、目录 TOC o 1-5 h z HYPERLINK l bookmark9 o Current Document 一、 方案1-直接控制法（多参数调节）2 HYPERLINK l bookmark36 o Current Document 二、 方案2-PID控制法（多参数调节，模型补偿）3 HYPERLINK l bookmark46 o Current Document 三、 方案3-无模控制法（单参数调节，模型补偿）5 HYPERLINK l bookmark63 o Current Document 四、 方案4-模式控制法（经验数据库，查表法）5 HYPERLINK l bookm

2、ark66 o Current Document 五、 方案5-PFC控制法（模型法，误差性能评价）6一、方案1-直接控制法（多参数调节）污水厂的曝气就是根据污水处理中的需氧量来调节鼓风机的风速。通常污水 处理中的需氧量由两部分组成，一是水中微生物的有氧活动所消耗的氧，二是一 定值的剩余氧，这是为了保证微生物高效有氧活动。那么大致可以得出，需氧量 等于微生物消耗的氧加上剩余的氧。需氧量表示为：（1-1）式中，0需为水处理中的需氧量，O耗为水中微生物的有氧活动所消耗的氧，0剩 为剩余氧。水中剩余的氧等于处理水体积乘以剩余氧浓度，即处理水体积乘以剩余溶解 氧。微生物消耗的氧等于微生物的数量乘以单个

3、所消耗的氧，微生物的数量等于 所需要处理水体积乘以单位体积微生物浓度，所需要处理水体积等于水流速度乘 以时间。表示如下，（1-2）（1-3）（1-4）（1-5）0 剩=V*D（1-2）（1-3）（1-4）（1-5）0耗=0单耗*NN = V*DmV 二T*F式中，0单耗为单个所消耗的氧，N为微生物的数量，V为处理水体积，Do，Dm 分别为水中剩余氧浓度和水中单位体积微生物浓度，F为水流速度，T为水流时 间。把（1-2）、（1-3）、（1-4）和（1-5）代入（1-1）中，表达如下，0需=0单耗 *T*F*D + T*F*D（1-6）曝气是指将空气中的氧强制向液体中转移的过程，溶解氧跟空气里氧的

4、分 压、大气压、水温和水质有密切的关系，因此在转移过程中会有一定氧的损失。 在此假设合适条件下可以完全溶于水中，氧损失可以忽略不计，那么可以认为曝 气量等于需氧量乘以氧气浓度。即，air需 = 0 需*D氧（1-7） 式中，air需为曝气量，D氧为氧气在空气中的浓度。污泥浓度就是单位体积污泥 含有的干固体重量，而水中的微生物来源于活性污泥中，那么可以认为水中单位 体积微生物浓度与污泥浓度有一定的比例关系，即（1-8）D = k*MLSS（1-8）把（1-6）和（1-8）代入（1-7）中，（1-9）air需=（0单耗 *k*MLSS + D）*T*F*D氧 令k广T*D氧，J。单耗*k，则（1-

5、9）（1-10）air =（k *MLSS + D）*k *F（1-10）需 2o 1即，（1-11）air需=% * Flow *（k2* MLSS + DO）（1-11）从式（1-11）可以看出，只要调整参数k1和k2既可以根据进水水量、污泥浓度 和剩余溶解氧调节曝气气量。为了减小误动作的影响，采用控制量增量法，即（1-11）u = % * AFlow * （七 * AMLSS + ADO） 为了保证控制达到期望值，对上述控制加上反馈调节，调节如下，（1-11）（1-12）Au = % * AFlow * （匕 * AMLSS + ADO）+ k3* e（1-12）二方案2-PID控制法（

6、多参数调节，模型补偿）PID作为工程中一种广泛应用的控制器，分为增量式和位置式。其区别有二： 一是增量型PID不需做累加，因误差产生的精度问题对控制量的影响较小；而位 置型PID则用到过去误差的累加，容易产生较大的累加误差，对控制量的影响较 大。二是增量型PID得出的是控制的增量，误动作影响小，必要时通过逻辑判断 限制或禁止本次输出，不会影响系统的工作；而位置型PID的输出是控制量的全 部输出，误动作影响大。根据曝气过程的特点，本为采用增量式PID算法。增量 式PID研究对象方程描述如下，Au(k) = Au(k-1) + k e(k) -e(k-1) + e(k) + 3(e(k) - 2e

7、(k -1) + e(k- 2) (2-1) pT T式中k为比例系数，T为积分时间常数，T为微分时间常数，T为采样周期。pid式(2-1)可以写成(2-2)式，Au(k) = Au(k-1) + k (e(k)-e(k-l)+k e(k) + k (e(k)-2e(k-1) + e(k-2) (2-2)pid式中kp、ki、kd分别为比例、积分和微分系数。式中的e(k)为实际值与期望值 的差值，如下，g(k) = e (k) = r(k)-y(k)( 2-3)式中r(k)为期望值，y(k)为当前输出，考虑到曝气过程的大滞后因素，本文对当 前输出值作补偿，从而使当前时刻的输出逼近因干扰和控制而

8、引起的真实输出 值。补偿方式如下，七(k) = y(k) + ym (k) - * (k-D)( 2-4)e(k) = r(k)-yp(k)(2-5)式中ym(k)为当前时刻模型值，D为滞后时间。通过曝气模型的预测输出值对当前时刻的输出值进行补偿，从而使当前时刻 的输出逼近因干扰和控制而引起的真实输出值。通过模型计算出当前时刻的输出 值和滞后时间前的输出值，用他们的差作为当前时刻输出值的补偿值。目前建立 曝气模型的方式有很多种，有非线性回归模型、BP神经网络模型和传递函数模 型等。回归模型建立最简单，它通过对数据的统计即可得到，但其逼近能力主要 取决于所选择结构的准确性。BP神经网络不需要确立

9、模型结构，其三层结构可 以逼近任意非线性，经过训练即可得到精度相当高的模型，但其结构复杂。传递 函数的建立则需要得到控制的单变量输出响应曲线，需要经过实验或数据的复杂 选择。根据曝气过程的特点，选择回归模型或BP神经网络模型最为恰当。通常污水厂中影响溶解氧的因素主要为进水水量、进水COD、曝气池的污泥 浓度和曝气气量，此外还有进水温度、进水PH、NH4-N、NO3-N、TP等。因此模 型的输入变量选择为进水水量、进水COD、曝气池的污泥浓度和曝气气量，输出 为曝气池的DO。三方案3无模控制法(单参数调节，模型补偿)无模型控制方法也是一个应用相当广泛的一个控制方法。其方法如下:$(k) = *k

10、-1) + 气&(k 1)(Ay(k)-i(k-1)Au(k-1)p + |Au(k-1)(3-1)(3-2)(3-3)u(k) = u(k-1) +X +Pg(k) y*(k+1)-y(k) $ (k)2i(k) (3-2)(3-3)u(k) = u(k-1) +X +Pg(k) y*(k+1)-y(k) $ (k)2其中，门上，PkG(0,2)，h，人是权重因子，是一个充分小的正数，i(1) 是i(k)的初值。无模型控制与PID控制一样，对于大时滞系统也会不稳定，那么 跟上述一样，也对其做补偿，补偿方式如下，七(k) = y (k) + y (k) -、(k-D)( 3-4), p $ (

11、k)u(k) = u(k-1) +kAy*(k+1)-y (k).+ $ (k)2p(3-5)式中ym(k)为当前时刻模型值，d为滞后时间。方案4-模式控制法(经验数据库，查表法)通常影响污水处理需氧量的主要因素有进水流量、进水COD、污泥浓度、水 温、PH和DO等。经过对椒江污水厂因素的分析，对需氧量有重大影响的变量是 进水流量、污泥浓度和DO值。通过经验查表法来进行曝气量的控制，表格结构 如表1，A1，A2，An为流量的区间，B1，B2，Bm为Mlss的区间。表2、 3、p+1(p=n*m)为流量在Ak区间和污泥浓度在Bq区间情况下溶解氧DO和曝气气量的关系对照表。表2、3、p+1(p=n

12、*m)的值可以取所在区间内 的平均值。根据对历史数据的统计分析，建立如下所示表格，供控制查询。表1进水流量和污泥浓度区间划分表XFlowMlssVA1A2AnB1B2Bm表2 A1-B1溶解氧DO和曝气气量的关系表DOa11a12a1nQb11b12b1n表3 A1-B2溶解氧DO和曝气气量的关系表DOa21a22a2nQb21b22b2n表p+1 An-Bm溶解氧DO和曝气气量的关系表DOap1ap2apnQbp1bp2bpn五.方案5-PFC控制法(模型法，误差性能评价)PFC控制方法对于曝气这种非线性、大时滞系统具有较好的效果，但是其模 型的精度相对要较高。预测函数控制研究对象方程描述如

13、下，y(k) =ay(k -1) + u (k -1)(5-1) TOC o 1-5 h z u (k) = r + ru(k) + r u2(k) +A + r un(k)(5-2)p012n式中n为多项式展开阶数。阶跃基加权预测函数控制的基函数结构形式如下，f(i)=七*l(t-iT) (i = 1,2,L ,H)(5-3)式中七(i = 1,2,A ,H)为权系数，H为预测时域长度。此时控制输入表达形式如下，u(k + i) = r * f (i) (i = 1,2,A , H)(5-4)式中R为组合系数。某一采样时刻其未来预测控制量用矩阵表达形式如下，U =x(k)* R(5-5)人(

14、k) = l1(k) 气(k) A、(k)(5-6)式中人(k)为k时刻权系数取值。在滚动优化过程中，可实时调节权系数人(k)的 值，保证参考轨迹与过程预测轨迹在预测时域内前部分的尽可能局部拟合。预测 函数控制中的控制量为up，其为基函数的组合，即式(5-4)中u的的计算公式。 用式(5-1)计算模型响应如下， TOC o 1-5 h z y (k + H) =ahy (k) + P(ah-1k +ah-2X +A +X )R (5-7) mm01H-1由预测函数控制原理知模型自由响应yt和模型函数响应yf结果如下，y (k + H) = a h y(k)(5-8)y (k + H) = y

15、(k + H) r(5-9)y (k + H) = P(ah-1k +ah-2k +A +k )(5-10)参考轨迹和过程预测轨迹结果如下，y (k + i) = c(k + i) + e- (c(k) - y (k)(5-11)y (k + i) = y (k + i) + e(k + i)(5-12)y (k + i) = y (k + i) + y (k + i)(5-13)e(k + i) = y (k) - y (k)(5-14)yp (k) = y (k) + ym (k) - y (k - D)(5-15)式中e(k + i)为过程补偿误差，D为过程时滞。优化计算就是用优化指标函数对过程预测轨迹与参考轨迹进行拟合，通常优化指