以溶解氧浓度作为SBR法模糊控制参数.doc

以溶解氧浓度作为SBR法模糊控制参数模糊控制在水处理领域已得到成功的应用，尤其是活性污泥法污水处理系统，它本身属于复杂的动态工程系统，无法找到精确的模型来描述，目前研究较多的是对出水悬浮物进行预测和控制的动态活性污泥法模糊控制1；在高纯氧活性污泥法中采取模糊控制，既能防止能量浪费，又能避免DO不足2。对硝酸态氮污染水脱氮处理的新方法生物电极法采用模糊控制，有利于避免过量地投加有机物，节省运行费用3、4。80年代以来，序批式活性污泥法（SBR法）用于处理间歇排放的水质水量变化很大的工业废水取得了很大成功并被广泛应用。SBR法的主要缺点是运行管理复杂，只有实现SBR法的自动控制，才能发挥其优势。传统的控制方法是时间程序控制和流量程序控制，即确定SBR法五个阶段所需要的时间后，实现其计算机控制。而工业废水的排放不仅是变化的或间歇的，而且其有机物浓度也随时间变化很大，往往相差几倍或十几倍。如果按相同的反应时间控制SBR的运行，当进水浓度高时出水不达标，当进水浓度低时曝气时间过长，浪费能源还易发生污泥膨胀。为了实现SBR法更高层次的计算机在线控制，必须寻找一个参数既能反映进水COD浓度的变化及反应过程中的降解情况，又能作为计算机控制参数，这也是深化SBR自动控制的重要研究课题。溶解氧浓度由于其能够在线检测、响应时间短、精确度高，人们在活性污泥法中围绕溶解氧浓度已做了大量研究，包括用溶解氧浓度作为SBR法过程控制和反应时间控制参数5，在脱氮反应过程中以在线检测的DO值模糊控制曝气量6等。但未见到采用SBR法处理工业废水时以在线检测的DO值作为模糊控制参数的研究。1 试验设备与方法SBR法试验装置如图1所示。反应器高70cm，直径30cm，总有效容积38L，采用鼓风曝气，转子流量计调节曝气量。在反应过程中在线检测DO值，并根据DO值的变化在一定的时间间隔内取样测定COD及混合液污泥浓度（MLSS）等指标。试验用石化废水中含有乙酸、偏苯三酸及苯酐等多种有机化合物，加适量的自来水稀释，使COD浓度在4001900mg/L范围内变化。由于原水中氮、磷元素含量极少，本试验按BOD5NP=10051的比例投加氯化铵（NH4Cl）和磷酸二氢钾（KH2PO4）配制的营养液，并调整pH值为66.5，水温控制在20 左右。2 试验结果与分析为了实现以DO作为SBR法运行过程中曝气量及反应时间的模糊控制参数，必须寻找DO与有机物降解之间的规律性。因此，本试验在两种不同条件下研究DO浓度对反应过程的影响。2.1 同一进水浓度不同曝气量试验配制原水COD浓度为1800mg/L，进水混合后反应初始COD浓度为1600mg/L。反应过程平均MLSS=2000 mg/L，曝气量分别为0.6m3/h、1.0m3/h、1.2m3/h时试验结果如图2所示。从图中可以看出，反应开始810min左右，不同曝气量下DO值显著不同。曝气量越小，此时的DO值越低，致使整个反应过程DO处于较低水平，大大延长了反应时间。随着曝气量的增大，初始DO值也在增大，带来整个反应过程DO的提高，缩短了有机物达到难降解程度的时间，即减少了反应时间。但是，过分地增大曝气量，初始DO的过高会造成整个过程DO处于过高的水平，却不能有效地缩短反应时间，还造成能量的巨大浪费。产生这种现象可以从耗氧速率与DO的相关关系的研究中得到解释，因为在低DO浓度下，DO浓度对生化反应速率的影响较大。当DO在12mg/L范围内，随着DO的提高，耗氧速率大幅度提高，标志着有机物降解速率的加快，从而缩短反应时间。当DO超过2mg/L后，继续增大DO值，由于受污泥浓度（MLSS）的限制，有机物降解速率的增加幅度较小。取进水COD浓度分别为400、600、800、1000、1300、1600、1900mg/L，每一进水浓度又在不同曝气量下进行试验，比较DO与有机物降解情况。试验中发现控制反应过程中DO为2.0 mg/L左右时最经济合理，其相应的曝气量为合适的曝气量。通过试验找出每一进水浓度下合适的曝气量，并应随着进水COD浓度的不同，根据检测的初始DO（810min左右）值调节曝气量，使整个反应过程的DO处于适当的水平，既能保证出水水质，又能节省运行费用。2.2 不同进水浓度相同曝气量试验配制原水COD浓度分别为850、1200、1500mg/L，进水混合后COD浓度分别为650、980、1300mg/L，反应过程平均MLSS为2000mg/L，曝气量为0.6 m3/h，试验结果如图3所示。从图可以看出，不同的进水COD浓度，反应10 min左右时DO值有很大差别，COD浓度越高，DO值越低，二者有很好的相关性。在COD浓度为650mg/L时，反应10 min左右DO值就升到4.5 mg/L；而COD浓度为1300 mg/L时，反应10 min时的DO值仅为1.3 mg/L。因此，在反应开始后较短的时间内就可以根据检测的DO值的大小预测出相应的进水COD浓度。用SBR法处理石化废水，以上述试验研究结果为基础，设定每一周期初始的曝气量均为0.6m3/h，在不知进水COD浓度的情况下，以在线检测反应10min左右时DO值的大小为依据，预测出该进水COD浓度，再找出在该进水COD浓度下适宜的曝气量，将其归纳总结如表1所示。与此同时，发现在上述每一试验过程中，当有机物达到难降解程度时，DO都有迅速大幅度升高的现象发生，并且在较短的时间内上升到56mg/L。根据反应期间DO的变化，实现对SBR供气量和反应时间的模糊控制。表1 进水COD浓度和反应时间预测初始DO值（mg/L） 预测进水COD浓度（mh/L） 适宜曝气量（m3/h） 反应时间（min） 5.554006000.20.39053.706008000.30.41003.702.2080010000.40.61002.201.50100013000.60.81201.50181.01201.151.0160019001.01.21303 DO作为模糊控制参数的基本思想SBR法的模糊控制目的有三：一是实现计算机自动控制；二是在保证出水水质前提下尽可能节省运行费用；三是避免曝气量不足或反应时间过长而引起的污泥膨胀。3.1 以初始DO作为模糊控制参数设定每一周期反应初始的曝气量均为0.6m3/h，以在线检测的反应初始阶段（810 min左右）的DO值（DOoff）作为被控制变量，以曝气量（u）为控制变量。在模糊控制系统的设计时，以在线检测的DOoff与设定的DO标准浓度（DOs）的偏差Ei作为模糊控制器的一个输入变量，1min为一个采样周期，一个采样周期后该偏差Ei的变化量CEi作为模糊控制器的另一个输入变量。Ei=DOoff-DOs i=1，2，3(1)CEi=Ei-Ei-1i=1，2，3(2)式中i-第i次采样的相应数据Ei-1-第i-1次采样处理水DOoff的偏差。根据这两个输入变量，经过模糊控制器的计算、判断与决策，作为模糊控制系统输出变量的则是控制变量的变化量ui，即曝气量的变化量。确定了模糊控制器的输入与输出变量后，根据模糊控制理论按照以下步骤实现模糊控制系统。3.1.1 精确量的模糊化根据表1所得的试验结果，将DO的偏差及偏差的变化量用模糊变量来表示，即将被控制变量进行模糊化处理，得到模糊集合向量。对误差E,误差变化CE及控制量u的模糊集及其论域定义如下：CE和u的模糊集均为：NB,NM,NS,O,PS,PM,PBE的模糊集为：NB,NM,NS,NO,PO,PS,PM,PB上述模糊语言变量的意义：NB=Negative Big(负大)PB=Positive Big(正大)NM=Negative Medium（负中）PM=Positive Medium(正中)NS=Negative Small(负小)PS=Positive Small(正小)NO=Negative Zero(负零)PO=Positive Zero(正零)O=ZeroE和CE的论域均为：-6，-5，-4，-3，-2，-1，0，+1，+2，+3，+4，+5，+6u的论域为：-7，-6，-5，-4，-3，-2，-1，0，+1，+2，+3，+4，+5，+6，+7将Ei、CEi、ui化为各自论域上离散的整型变量，即整型化处理。在此设定DOs=2.0mg/L,是相对于COD=1000 mg/L、曝气量为0.6m3/h、反应9min时的DO浓度。E、CE及u的模糊集和论域确定后，须对模糊语言变量确定隶属函数，即所谓对模糊变量赋值，就是确定论域内元素对模糊语言变量的隶属度。上述论域E、CE、u上的模糊变量均假定为正态型模糊变量，其正态函数为：F(x)=exp-(x-a)/)2(3)此函数确定了模糊隶属函数曲线的形状。将确定的隶属函数曲线离散化，就得到了有限个点上的隶属度，便构成了一个相应的模糊变量的模糊子集。3.1.2 建立模糊控制规则在试验基础上，分析DO与有机物降解及曝气量之间的关系，建立以模糊语言表示的模糊控制推理的合成规则和模糊控制规则。根据操作过程中可能遇到的各种情况和系统的运行数据，将相应的控制策略归纳为表2，这是一组根据系统输出的误差及误差的变化趋势来消除误差的模糊控制规则。表2 SBR法曝气过程模糊控制规则表ECENBNMNSOPSPMPBuNBPBPBPBPBPMPSONMPBPBPBPMPSONSNSPBPMPMPSONSNMNOPMPMPSONSNSNMPOPMPSPSONSNMNMPSPMPSONSNMNMNBPMPSONSNMNBNBNBPBONSNMNBNBNBNB上述这些模糊控制规则可以用模糊条件语句来描述，例如：if E=NB or NM and CE=NB or NM then u=PB orif E=NB or NM and CE=NS then u=PB or 上述选取控制量变化的原则是：当误差大或较大时，选择控制量以尽快消除误差为主；而当误差较小时，选择控制量要注意防止超调，以系统的稳定性为主要出发点。例如，当DO误差及误差变化均为负大时，就是SBR反应器内DO浓度很低，而且有进一步降低的趋势，如果不加以调整，势必造成反应时间过长或引起污泥膨胀，为尽快提高DO浓度，消除误差，必须增大曝气量，所以u取正大。当误差为负小，误差变化为正小时，系统本身具有消除误差的能力，可以不调整曝气量。3.1.3 模糊推理及其模糊量的非模糊化在模糊控制规则的指导下，经过模糊决策后，得到模糊控制变量ui：为了对被控制对象SBR处理系统施加精确的控制，还需要将模糊控制变量ui转化为可执行的精确量，即曝气量变化量的准确量，这就是非模糊化处理过程。上述模糊控制规则所确定的每一条模糊条件语句都可以计算出相应的模糊控制量ui。例如，由第1条语句所确定的模糊关系可用式（4）表示：R=(NBE+NME)PBu(NBCE+NMCE)PBu(4)如果此刻采样所得到的实际误差模糊变量为e，误差变化的模糊变量为ce,根据推理的合成规则进行模糊决策，得到模糊控制量为u1:u1=eo(NBE+NME)PBuceo(NBCE+NMCE)PBu(5)同理，可由其余模糊条件语句计算出相应的模糊控制量u2,u3,由于各条件语句之间是或的关系，则控制量的模糊集合u表示为：u=u1+u2+un(6)由式（6）所计算出的控制量是一个模糊子集，不能直接应用于被控对象，必须经过非模糊化处理转化为精确量。本文采用加权平均法作为非模糊化处理方法，计算式如下：这种方法可以充分利用模糊推理结果、模糊子集提供的有用信息量，得到SBR法曝气过程的模糊控制表，储存在计算机中。上述过程是反应进行到89min，根据在线检测的DO值预测进水COD浓度，并在第10min对曝气量进行调整。为了保证较好的控制效果，可在后续反应过程中继续检测DO值，并对曝气量进行再一次调整。调整次数的多少应视预测进水COD浓度的大小及相应反应时间的长短而定。在有机物达到难降解程度时，DO迅速大幅度升高，这是停止曝气的信号。如果在反应过程中，频繁地调整曝气量，势必使DO始终维持在2.0mg/L左右，而不会发生DO迅速大幅度升高的现象，进而影响反应时间的控制。因此，试验中根据表1的试验结果，在反应810 min预测到进水COD浓度后，便相应地预知反应时间。如果预测的反应时间在120min之内，进行两次调整，第二次是在反应进行到30 min，DOs仍取2.0 mg/L，具体的控制方法同前。如果预测的反应时间超过120 min，则进行三次调整，第三次是在反应进行到60 min，DOs仍取2.0 mg/L，具体的控制方法同前。经过这样的方法控制，就可以弥补第一次调整曝气量（810 min）所造成的误差，而且不会影响对反应时间的控制。3.2 作为反应时间的控制参数在反应初期（810 min），以上述的模糊控制方法对曝气量进行调整后，使后续反应过程中DO处于合适的水平。由于SBR法间歇运行的特点，当有机物达到难降解程度时，DO迅速大幅度升高，在10 min内DO便可上升到3.5 mg/L,反应结束时DO可升到6 mg/L左右。这一变化特点可用模糊语言变量加以描述，以此作为停止曝气的信号，故引入了用DO误差的大小和误差变化的快慢作为模糊控制器的两个输入变量。SBR法处理石油化工废水的试验结果表明：当DO(升高的幅度)3.5 mg/L，且DO的变化速率0.30 mg/(Lmin)时，反应已结束，应停止曝气，这是建立模糊控制规则的重要依据。由于两种现象中误差及误差的变化均为正，可以简化对误差E、误差变化CE及控制量u的模糊集和论域的定：E、CE的模糊集均为：PS，PM，PBE、CE的论域为：1，2，3，4，5，6而对控制量u（即曝气量）而言，只有两种选择：一种是维持原来的曝气量不变，继续等待；另一种是立即停止曝气，不存在改变曝气量大小的问题。对输入变量E和CE进行模糊化处理，此时DO的设定值DOs取2.0 mg/L。模糊化方法与3.1.1节中阐述的方法相同，不再重复。在此仅对模糊规则的建立进行说明，模糊规则控制如表3所示。表3 SBR反应时间模糊控制表ECEPSPMPBuPS000PM000PB0注 表示停止曝气当DO的误差为PS和PM时，无论DO的变化速率如何，均维持原有的曝气量，避免因曝气时间不够而使出水达不到排放标准。只有当E达到PB,且CE达到PM或PB，才认为有机物不再被降解，应该停止曝气。将3.1、3.2阐述的内容结合起来，便可以DO作为控制参数，实现对SBR法曝气过程和反应时间的模糊控制。当然，任何一种控制方法都无法保证每一次预测和控制都是准确无误的，有时难免会出现异常现象。例如：根据810min的DO值预知反应时间为120min左右，可是在90min就出现了DO迅速大幅度升高的现象，说明COD的预测值大于实际值。如果这种情况发生，此时不应该再调整曝气量，为防止因较早停止曝气而使出水不达标，只能改变反应时间控制规则中的某些参数，使控制规则更为严格。本试验中DO升高的幅度由正常时的3.5 mg/L提高到4.5 mg/L，DO变化速率也由0.30 mg/(Lmin)提高到0.35 mg/(Lmin)，即只有Ei4.5 mg/L,CEi0.35 mg/(Lmin)，才可以停止曝气。还有另外一种相反的情况，就是预测的反应时间为90 min，而实际运行至120 min DO仍没有迅速大幅度升高，说明COD的预测值小于实际值。如果这种情况发生，则应放宽控制规则，将Ei由3.5mg/L降到2.5 mg/L,CEi由0.30 mg/(Lmin)降到0.20 mg/(Lmin)，即Ei2.5 mg/L,CEi0.20mg/(Lmin)，就可以停止曝气。因此，在对反应时间进行控制时，既要根据在线检测的DO变化情况，也要参考前面的预测结果，如果有上述的异常情况发生，就要采取不同的控制规则，保证控制效果和系统的正常运行。4 结论 采用SBR法处理石油化工废水，有机物降解初期（810 min）的溶解氧浓度对整个反应过程有重要影响。如果初始DO值过高必将导致反应过程DO值的普遍过高，不但不会有效地缩短反应时间，反而会增大运行费用。相反，如果 DO值过低，则延长了反应时间，而且容易引起污泥膨胀。 可以用反应初始阶段溶解氧浓度作为SBR法曝气量的模糊控制参数。假定每一反应周期初始曝气量相同，可根据反应开始后较短时间内（810 min）DO值的变化情况预测进水COD浓度，进而调整到该浓度下适宜的曝气量。在反应过程中，也可以根据DO变化情况对曝气量进行再一次调整。 SBR反应器内，当有机物达到难降解程度时，DO迅速大幅度升高，这一变化特点可用模糊语言变量加以描述，实现对SBR反应时间的模糊控制。银杏酮生产废水的治理 银杏树嫩叶可生产出我国特产的软化心脑血管，降血脂的良药银杏酮。一般每生产 1 kg银杏酮，需消耗 3 吨银杏嫩叶。大量的残渣夹杂着生产过程中加入的助剂排出，污染最重的银杏酮废水 CODcr 质量浓度达 8.54 104 mgL以上。l 废水来源和排放标准1.l 废水来源银杏酮生产废水属于中成药中单一成分醇提废水。它来源于设备清洗、柱再生和酸碱冲洗、蒸发冷凝液等。废水中含有乙醇、氨基酸、植物纤维、脂肪质、叶绿素以及 NaOH 等。主要废水来源及污染物含量见表1。表1 银杏国废水主裹来源的污染物检测结果废水来源（CODcr）/（mgL-1）（SS）（mgL-1）pH值黄酮上往8.54l04 4.25水洗去杂3.92l0430.5l043.76酸再生柱6.80l04485.51042.64碱再生柱3.32l0493510411.74综合废水3.02l04483.71047.63经检测综合废水 （CODcr）可达 3.02 104mg/L，m（BOD5）m（CODcr）0.330.38，生化性一般。1.2 设计水质水量废水水量：45 m3d，处理能力 5 m3h。废水水质：（CODcr）3.02 104 mgL，（BOD5）= 11 000 mgL，（SS）= 483.7 mgL，pH值 7.63。1.3 排放标准根据江苏省及邓州市环境保护局规定，排放标准执行国家污水综合排放标准（GB89781996）中新建企业二级标准：（CODcr）150 mgL，（BOD5）60 mgL，（SS） 100 mgL，pH值 69。2 处理工艺银杏酮生产废水中含有银杏叶的植物纤维和脂肪质。生物稳定性强，不易生物降解；同时由于是全厂综合废水含有锅炉房的酸性排污水。因此采用pH值调节，并结合加药混凝沉池工艺，首先去除废水中的植物纤维和脂肪质，兼具化学水解效果，使UASB进水浓度减少一半以上，CODcr质量浓度可降低到11 114 mgL以下。废水中溶解性成分主要是乙醇、氨基酸、叶绿素等，易于生化分解。因此选用目前制药行业高浓度废水普遍使用的 UASB艺；并在 UASB 池的中上部放置 2m 高的组合填料，构成厌氧复合床。厌氧出水采用推流式延迟接触氧化池，作为达标出水的末端处理工艺。银杏酮废水处理工艺流程见图1。3 主要设备和构（建）筑物主要设备见表2，主要构建（筑）物见表3。表2 主要设备序号设备名称型号、规格特性参数数量配套电机功率/kw备注1格栅长 500 mm 宽 400mm 1台 2污水泵32ZW1020流量 10 m3/h 扬程 20 m3台2.2 3泥浆泵l1B3（b）流量 10 m3/h 扬程 60 m1台7.5 4罗茨风机SSR80风量 3.36 m3/min风压 53.9 kPa2台5.5 5絮凝反应器1600 mm 3150 mm 1台0.37减速机L-D0 反应30min6板框压滤机BMS20630/25过滤面积20 m21台 7气液分离器 600 mm 600 mm 1台 8酸贮槽 800 mm l200 mm 1台 9组合填料l50 mm 100 mm 200m3 表3 主要构（建）筑物序号构（建）筑物名称规格有效容积/m3停留时间/h其它工艺参数1调节池5100 mm 5100 mm 4650 mm5024 2中间水池5100 mm 5100 mm 4650 mm504 3UASB池 7 910 mm 7 280 mm20040容积负荷（以 BOD计）2.5 kg/（m3d）组合填料 80 m34厌氧出水池5100 mm 5100 mm 4650 mm504 5初沉池2500 mm 2500 mm 5500 mm302表面负荷 1.6 m3/（m2h）6生物选择池2250 mm 2500 mm 5500 mm282.5 7接触氧化池6000 mm 5000 mm 4500 mm 2250 mm 2500 mm 4500 mm142.512填料负荷（以 BOD计）012 kd/（m3d）,填料 120 m3，水气比 1：408二沉池2500 mm 2500 mm 5100 mm222表面负荷l6m3/（m2h）9物化处理间6740 mm 6520 mm 5800