曝气控制方案综述.doc

目录一、方案 1-直接控制法（多参数调节） ............................................................................................1二、方案 2-PID 控制法（多参数调节，模型补偿） ..........................................................................3三、方案 3-无模控制法（单参数调节，模型补偿） .........................................................................4四、方案 4-模式控制法（经验数据库， 查表法） .............................................................................5五、方案 5-PFC 控制法（模型法，误差性能评价） .........................................................................5一、方案 1-直接控制法（多参数调节）污水厂的曝气就是根据污水处理中的需氧量来调节鼓风机的风速。

通常污水处理中的需氧量由两部分组成，一是水中微生物的有氧活动所消耗的氧，二是一定值的剩余氧，这是为了保证微生物高效有氧活动那么大致可以得出，需氧量等于微生物消耗的氧加上剩余的氧需氧量表示为：（1-剩耗需 O1）式中， 为水处理中的需氧量， 为水中微生物的有氧活动所消耗的氧，'O需 耗为剩余氧剩水中剩余的氧等于处理水体积乘以剩余氧浓度，即处理水体积乘以剩余溶解氧微生物消耗的氧等于微生物的数量乘以单个所消耗的氧，微生物的数量等于所需要处理水体积乘以单位体积微生物浓度，所需要处理水体积等于水流速度乘以时间表示如下，（1-oD*VO剩2）（1-N单 耗耗3）（1-mD*V4）（ 1-FT5）式中， 为单个所消耗的氧， 为微生物的数量， 为处理水体积， ，单 耗ONVoD分别为水中剩余氧浓度和水中单位体积微生物浓度， 为水流速度， 为水mD FT流时间把（1-2） 、 （1-3） 、 （1-4）和（1-5）代入（1-1）中，表达如下，（1-om*TD*FO单 耗需6）曝气是指将空气中的氧强制向液体中转移的过程，溶解氧跟空气里氧的分压、大气压、水温和水质有密切的关系，因此在转移过程中会有一定氧的损失。

在此假设合适条件下可以完全溶于水中，氧损失可以忽略不计，那么可以认为曝气量等于需氧量乘以氧气浓度即，（1-氧需需 D*Oair7）式中， 为曝气量， 为氧气在空气中的浓度污泥浓度就是单位体积污泥需air氧D含有的干固体重量，而水中的微生物来源于活性污泥中，那么可以认为水中单位体积微生物浓度与污泥浓度有一定的比例关系，即（1-MLS*km8）把（1-6）和（1-8）代入（1-7）中，（1-氧单 耗需 ）（ D*FTS*kOair o9）令 ， ，则氧D*Tk1k2单 耗（1-F*kDMLS*air 1o）（需 10）即，（1-）（需 OS*kFlow*kair21 11）从式（1-11）可以看出，只要调整参数 k1 和 k2 既可以根据进水水量、污泥浓度和剩余溶解氧调节曝气气量为了减小误动作的影响，采用控制量增量法，即（1-）（ DOMLS*kFlow*k21 u11）为了保证控制达到期望值，对上述控制加上反馈调节，调节如下， （1-eku *DOMLS*kFlow*k 321 ）（12）二、方案 2-PID 控制法（多参数调节，模型补偿）PID 作为工程中一种广泛应用的控制器，分为增量式和位置式。

其区别有二：一是增量型 PID 不需做累加，因误差产生的精度问题对控制量的影响较小；而位置型 PID 则用到过去误差的累加，容易产生较大的累加误差，对控制量的影响较大二是增量型 PID 得出的是控制的增量，误动作影响小，必要时通过逻辑判断限制或禁止本次输出，不会影响系统的工作；而位置型 PID 的输出是控制量的全部输出，误动作影响大根据曝气过程的特点，本为采用增量式PID 算法增量式 PID 研究对象方程描述如下，（2-2)]-e(k1)-2((ek)T)1)-e(k[1)-u(k) dip 1）式中 为比例系数， 为积分时间常数， 为微分时间常数， 为采样周期piTdT式（2-1）可以写成（2-2）式，（2-2)-e(k1)-2(ek)()k1)-e(k1)-u(k) dip 2）式中 、 、 分别为比例、积分和微分系数式中的 为实际值与期望pid ()值的差值，如下，（2-y(k)-r(k)e3）式中 为期望值， 为当前输出，考虑到曝气过程的大滞后因素，本文对r(k)y(k)当前输出值作补偿，从而使当前时刻的输出逼近因干扰和控制而引起的真实输出值补偿方式如下，（2-D)-(k)(k)mpyy4）（2-(k)y-r(k)pe5）式中 为当前时刻模型值， 为滞后时间。

k)ymD通过曝气模型的预测输出值对当前时刻的输出值进行补偿，从而使当前时刻的输出逼近因干扰和控制而引起的真实输出值通过模型计算出当前时刻的输出值和滞后时间前的输出值，用他们的差作为当前时刻输出值的补偿值目前建立曝气模型的方式有很多种，有非线性回归模型、BP 神经网络模型和传递函数模型等回归模型建立最简单，它通过对数据的统计即可得到，但其逼近能力主要取决于所选择结构的准确性BP 神经网络不需要确立模型结构，其三层结构可以逼近任意非线性，经过训练即可得到精度相当高的模型，但其结构复杂传递函数的建立则需要得到控制的单变量输出响应曲线，需要经过实验或数据的复杂选择根据曝气过程的特点，选择回归模型或 BP 神经网络模型最为恰当通常污水厂中影响溶解氧的因素主要为进水水量、进水 COD、曝气池的污泥浓度和曝气气量，此外还有进水温度、进水 PH、NH4-N、NO3-N、TP 等因此模型的输入变量选择为进水水量、进水 COD、曝气池的污泥浓度和曝气气量，输出为曝气池的 DO 三、方案 3-无模控制法（单参数调节，模型补偿）无模型控制方法也是一个应用相当广泛的一个控制方法其方法如下：（3-1)-u(k)-(yk)(1)-(k1)-(k) 2u1）（3- 1)-(k(k)(),k21 u或若 2） （3-y(k)]-1([(k)1)-(k)*2u3）其中， ， （0,2 ） ， ， 是权重因子， 是一个充分小的正数，k是 的初值。

无模型控制与 PID 控制一样，对于大时滞系统也会不稳定，(1)那么跟上述一样，也对其做补偿，补偿方式如下，（3-D)-(k)(k)mpyy4）（3-()]-1([(k)1)-(k) p*2u5）式中 为当前时刻模型值， 为滞后时间k)ymD四、方案 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*m）的值可以取所在区间内的平均值根据对历史数据的统计分析，建立如下所示表格，供控制查询表 1 进水流量和污泥浓度区间划分表FlowMlssA1 A2 … AnB1B2…Bm表 2 A1-B1 溶解氧 DO 和曝气气量的关系表DO a11 a12 … a1nQ b11 b12 … b1n表 3 A1-B2 溶解氧 DO 和曝气气量的关系表DO a21 a22 … a2nQ b21 b22 … b2n...表 p+1 An-Bm 溶解氧 DO 和曝气气量的关系表DO ap1 ap2 … apnQ bp1 bp2 … bpn五、方案 5-PFC 控制法（模型法，误差性能评价）PFC 控制方法对于曝气这种非线性、大时滞系统具有较好的效果，但是其模型的精度相对要较高。

预测函数控制研究对象方程描述如下，（5-1）)1()1()kukyp（5-2）((20 rrrkunp 式中 为多项式展开阶数阶跃基加权预测函数控制的基函数结构形式如下，n（5-3）()*1()(1,)iftiTiH式中 为权系数， 为预测时域长度此时控制输入表达形式如,21(Hi下，（5-4）),21()*)(ifiku式中 为组合系数某一采样时刻其未来预测控制量用矩阵表达形式如下，（5-5）*)(kU（5-6）)()(21kkH式中 为 时刻权系数取值在滚动优化过程中，可实时调节权系数 的)(k )(k值，保证参考轨迹与过程预测轨迹在预测时域内前部分的尽可能局部拟合预测函数控制中的控制量为 ，其为基函数的组合，即式（5-4）中 的的计算公pu u式用式（5-1）计算模型响应如下，（5-7） )())( 11201 HHmHmkyky 由预测函数控制原理知模型自由响应 和模型函数响应 结果如下，l fy（5-8）)()(kykyHl（5-9）of （5-10）)()( 112010 HHky参考轨迹和过程预测轨迹结果如下，（5-11）)(()()( kyceikci pir（5-12） iymp（5-13）)()()(kiikyflm（5-14）yemp（5-15）)()()Dkkyp式中 为过程补偿误差， 为过程时滞。

优化计算就是用优化指标函数对)(ikeD过程预测轨迹与参考轨迹进行拟合，通常优化指标函数取为误差平方和，公式如下，（5-16）Hi prikyiJ12)]()([令 ，则优化指标函数形式改为，)()()( keyikidlr （5-17）Hi oikyidJ12])()([极小化优化指标函数，则组合系数计算公式如下，（5-18）)()(11 ikyiikyidoHiooHi 那么预测时域内第一步预测控制量结果如下，（5-19）1)(u。