CN104111618A - 一种活性污泥污水处理曝气环节的建模、滑模控制器及其控制方法 - Google Patents
一种活性污泥污水处理曝气环节的建模、滑模控制器及其控制方法 Download PDFInfo
- Publication number
- CN104111618A CN104111618A CN201410231030.3A CN201410231030A CN104111618A CN 104111618 A CN104111618 A CN 104111618A CN 201410231030 A CN201410231030 A CN 201410231030A CN 104111618 A CN104111618 A CN 104111618A
- Authority
- CN
- China
- Prior art keywords
- mrow
- sliding mode
- msub
- mtd
- concentration
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000005273 aeration Methods 0.000 title claims abstract description 84
- 239000010865 sewage Substances 0.000 title claims abstract description 72
- 239000010802 sludge Substances 0.000 title claims abstract description 72
- 238000000034 method Methods 0.000 title claims abstract description 65
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 claims abstract description 64
- 239000001301 oxygen Substances 0.000 claims abstract description 64
- 229910052760 oxygen Inorganic materials 0.000 claims abstract description 64
- 230000008569 process Effects 0.000 claims abstract description 33
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 24
- 238000005842 biochemical reaction Methods 0.000 claims abstract description 7
- 238000013178 mathematical model Methods 0.000 claims description 65
- 244000005700 microbiome Species 0.000 claims description 16
- 238000004062 sedimentation Methods 0.000 claims description 16
- 239000005416 organic matter Substances 0.000 claims description 13
- 239000011159 matrix material Substances 0.000 claims description 12
- 239000002351 wastewater Substances 0.000 claims description 7
- 238000004422 calculation algorithm Methods 0.000 claims description 6
- 238000004891 communication Methods 0.000 claims description 6
- 239000007788 liquid Substances 0.000 claims description 6
- 239000000463 material Substances 0.000 claims description 6
- 239000000758 substrate Substances 0.000 claims description 6
- 238000012937 correction Methods 0.000 claims description 5
- 238000006065 biodegradation reaction Methods 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000000926 separation method Methods 0.000 claims description 4
- 239000006228 supernatant Substances 0.000 claims description 4
- 230000000813 microbial effect Effects 0.000 claims description 3
- 230000006641 stabilisation Effects 0.000 claims description 3
- 238000011105 stabilization Methods 0.000 claims description 3
- 230000009471 action Effects 0.000 claims description 2
- 238000004904 shortening Methods 0.000 claims 1
- 238000005202 decontamination Methods 0.000 abstract 1
- 230000003588 decontaminative effect Effects 0.000 abstract 1
- 238000004134 energy conservation Methods 0.000 abstract 1
- 238000013459 approach Methods 0.000 description 4
- 238000001514 detection method Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 238000011161 development Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 239000007787 solid Substances 0.000 description 3
- 230000003111 delayed effect Effects 0.000 description 2
- 238000009795 derivation Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000003912 environmental pollution Methods 0.000 description 2
- 239000010842 industrial wastewater Substances 0.000 description 2
- 239000002699 waste material Substances 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000002950 deficient Effects 0.000 description 1
- 239000013505 freshwater Substances 0.000 description 1
- 239000007789 gas Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 239000012263 liquid product Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 239000002957 persistent organic pollutant Substances 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 239000000047 product Substances 0.000 description 1
- 238000005067 remediation Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 239000012265 solid product Substances 0.000 description 1
- 238000005728 strengthening Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 238000003911 water pollution Methods 0.000 description 1
Landscapes
- Activated Sludge Processes (AREA)
Abstract
本发明属于自动控制领域,具体涉及一种活性污泥污水处理曝气环节的建模、滑模控制器及其控制方法。该建模和控制方法可求出污水处理曝气环节过程中,任一阶段曝气池中所需溶氧量的大小,从而实现对曝气池中溶氧量的自动化控制,方法包括曝气池生化反应过程的建模和溶氧浓度的滑模控制方法两个部分,本发明提高了活性污泥法污水处理曝气环节的去污能力,提升了处理后的水质,并且溶解氧浓度的实时控制避免了能源的浪费,更加节能。
Description
技术领域
本发明属自动控制技术领域,具体涉及一种活性污泥污水处理曝气环节的建模、滑模控制器及其控制方法。
背景技术
现如今,环境污染问题已经成为影响人类社会可持续发展的主要问题之一,而在环境污染问题当中污染最严重的就是水资源的污染,尤其是在发展中国家,水资源的污染最为严重。目前,对我国而言,淡水资源十分匮乏,但是由于我国工业化过程中某些技术的不成熟,进一步造成了水资源的污染。迄今为止,加强水资源的治理以及预防水资源的污染成为世界各国关注的焦点。我国由于人口众多,人均所占有的水资源更为短缺,在我国,造成水资源污染主要包括两方面原因,一方面是没有经过任何处理或者是经过处理但是没有达到排放标准的工业废水;另一方面是人们日常生活造成的生活污水。工业废水可以通过建立严格的监督制度以及加强治理使其排放量减少,但是由于我国生活污水处理技术远远落后与西方发达国家,生活污水成为造成我国水资源污染的主要因素。因此,提高治理生活污水的技术以及加强生活污水处理系统的智能控制成为我国改善水资源的污染首先要解决的问题。
发明内容
本发明为解决现在城市生活污水处理方法及控制系统的开发方面均比较落后的问题,提供一种活性污泥污水处理曝气环节的建模、滑模控制器及滑模控制方法。本发明方法简单易行,具有很好的动态品质和对外部干扰具有很强的适应能力,提高了控制系统对不确定参数、外界扰动及滞后性的鲁棒性能。
本发明通过以下技术方案予以解决:
一种活性污泥污水处理曝气环节的建模,其特征是包括以下步骤:
a.采用活性污泥污水处理曝气环节的工艺流程,即废水首先进入初次沉淀池沉淀,然后在曝气池中进行生物降解,然后在二次沉淀池进行固液分离,分离出上清水排出,底部污泥部分回流污泥至曝气池,剩余活性污泥处理后排出;
b.根据活性污泥1号模型(ASM1)和微生物的生长动力学方程、物料平衡定律,将步骤a中工艺流程建立活性污泥法污水处理曝气环节的数学模型,利用滑模控制器对其进行智能控制,其数学模型为:
其中,X是微生物浓度(grVSS/m3),S是可降解有机物浓度(grBOD5/m3),O是溶解氧浓度(grO2/m3),Xi是进水微生物浓度(grVSS/m3),Si是进水底物的浓度(grDOB5/m3),表示的是微生物的最大比增长速率(2/day),ks是半速度常数(grBOD5/m3),ko是半饱和系数(2grO2/m3),kd是子衰减系数(0.05/day),Q是进水流量(m3/day),V表示的是曝气池的体积(m3),Qw表示排出污泥的流量(m3/day),C是二次沉淀池中污泥的浓缩系数,f是联系有机物与需氧量的一个因子(0.67grBOD5/COD),fx是消耗因子(1.42grCOD/grVSS),u是曝气量(grO2/m3day),YNH是观察生长系数。
c.在废水进入曝气池初始阶段,可降解有机物浓度高于半速度常数,溶解氧浓度高于半饱和系数时,对数学模型进行推导简化得到:
其中,δ是对溶解氧设置的一个充量系数。
d.随着微生物生化反应的进行,可降解有机物逐渐被消耗,浓度在不断地降低,当可降解有机物浓度低于半速度常数,溶解氧浓度低于半饱和系数时,数学模型简化为:
e.以底物浓度的大小将活性污泥污水处理系统分为两个阶段,即可降解有机物浓度高于半速度常数和可降解有机物浓度低于半速度常数,改进的数学模型变换成状态空间表达式的形式:
其中,
B=[0 0 1]T,C=[0 0 1],D=0
优选的,所述的活性污泥污水处理曝气环节的建模,在加入不确定、时滞以及外界干扰因素时,数学模型为:
其中,ΔAσ、ΔAdσ为系统的不确定性,d为常时滞,w为外界干扰。
一种用于上述的活性污泥污水处理曝气环节的建模的滑模控制器,是将建立的数学模型应用滑模控制算法设计活性污泥污水处理曝气环节的滑模控制器,所述滑模控制器是根据系统所期望的动态特性来设计系统的切换超平面(滑模面),通过滑模控制器使系统状态从超平面之外向滑模面收束,系统一旦到达滑模面,控制作用将保证系统沿滑模面到达系统原点,由于系统的特性和参数只取决于设计的滑模面而与外界干扰没有关系,所以滑模变结构控制具有很强的鲁棒性。从而实现对切换数学模型的精确控制,使稳定时间缩短,控制精度提高,抗干扰能力和鲁棒性提高;
所述滑模控制器包括两部分,其中一部分是等效控制器ueq,用于保证从任意状态出发的系统始终满足滑模到达条件;另一部分是校正控制器uN,用于抵消系统中存在的各种不确定性和外界干扰。
所述滑模控制器的设计,具体步骤为:
所述滑模控制的滑模面的选取s(t)=BTPx(t);
其中B是在加入不确定、时滞以及外界干扰因素时,所建立的数学模型中的输入矩阵,P是正定对称矩阵,x(t)是权利要求2中所建立的数学模型中的状态向量。
所述滑模控制器的选取:
其中u(t)是滑模控制器,其作用是控制活性污泥污水处理曝气环节溶解氧浓度,uN=-(λ||x(t)||+k+μ)sign(s(t)),λ,k,μ是正常数,sign(s(t))是关于滑模面函数的符号函数。
所述滑模控制器的设计,具体步骤为:
第一步,滑模面的选取,
第二步,滑模控制器的选取
第三步,滑模控制器控制效果的稳定性证明。
1.针对第一步,滑模面选取,选取的依据是,保证存在某一区域,使得状态点趋近该区域时,就会被吸引到该区域内运动,该区域就是滑模面。按照滑模面必须能够吸引状态点的要求,滑模面必须满足以下条件,此条件称为滑模到达条件。
其中s代表滑模面。当s从大于0方向趋近于0时,满足的极限小于0,当s从小于0方向趋近于0时,满足的极限大于0,即。
此处针对在加入不确定、时滞以及外界干扰因素时所建立的数学模型,选 取滑模面s(t)=BTPx(t),保证滑模面以外的状态点都将在有限时间内到达滑模面。
其中B是在加入不确定、时滞以及外界干扰因素时所建立的数学模型中的输入矩阵,P是正定对称矩阵,x(t)是在加入不确定、时滞以及外界干扰因素时,所建立的数学模型中的状态向量。
2.针对第二步,滑模控制器的选取,选取的依据是,保证系统的状态点一旦到达滑模面,此状态点将只能沿滑模面运动,在有限时间内达到期望值。滑模控制器应该包括两部分,其中一部分是等效控制器ueq,作用是保证从任意状态出发的系统始终满足滑模到达条件;另一部分是校正控制器uN,其作用是抵消系统中存在的各种不确定性和外界干扰。
针对上一步中所选取的滑模面,滑模控制器选取如下,
其中u(t)是滑模控制器,其作用是控制活性污泥污水处理曝气环节溶解氧浓度,uN=-(λ||x(t)||+k+μ)sign(s(t)),λ,k,μ是正常数,sign(s(t))是关于滑模面函数的符号函数。
3.针对第三步,滑模控制器控制效果的稳定性证明。选取合适的李雅普诺夫函数V(t),使其满足的条件。证明如下:
其中Q是适当维数的常数矩阵。
将选取的李雅普诺夫函数对时间t求导,可以得到
将在加入不确定、时滞以及外界干扰因素时所建立数学模型中的代入上 式,得到
从上式结果中可以看出,由于,则此控制器u(t)能够保证活性污泥污水处理曝气环节溶解氧浓度渐近稳定到期望值。从而实现一种活性污泥污水处理曝气环节的建模、滑模控制器控制。
一种使用上述述滑模控制器对活性污泥污水处理曝气环节进行控制的方法,包括以下步骤:
步骤一,根据ASM1模型和微生物的生长动力学方程、物料平衡定律化简ASM1模型,得到化简数学模型;
步骤二,根据可降解有机物浓度的高低,将步骤一ASM1化简模型分为两个子系统,即可降解有机物浓度高于半速度常数子系统和可降解有机物浓度低于半速度常数子系统;
步骤三,利用切换系统理论将步骤二两个子系统结合为统一的数学模型,不考虑不确定性、时滞以及外界干扰因素,得到简化数学模型;
步骤四,根据污水处理实际情况,对步骤三简化数学模型,加入不确定性、时滞以及外界干扰因素,得到切换数学模型;
步骤五,根据步骤四切换数学模型,设计滑模控制器,编写MATLAB滑模控 制器子程序;
步骤六,开启工控机(为上位机)和PLC控制装置,利用组态王DDE通讯,定时在组态王中调用MATLAB滑模控制器子程序,实现滑模控制,实时在线调整溶解氧浓度;
步骤七,循环执行步骤六,即可得到溶解氧浓度实时控制曲线。
所述的一种对活性污泥污水处理曝气环节进行滑模控制的方法,包括以下具体步骤:
滑模控制采用上位机与下位机相结合的控制方式,所述上位机为安装有组态王(简称组态王)和MATLAB软件(简称MATLAB)的工业计算机,所述的下位机采用西门子公司的S7-315-2DP PLC,在本曝气环节自动控制过程中,溶氧仪采集的溶解氧浓度信号送给PLC,PLC以变量形式将溶解氧浓度信号送至上位机组态王,上位机组态王利用DDE通讯将其传给MATLAB,MATLAB运行根据所设计的滑模控制器编写的滑模控制程序,计算所需控制量u(t),计算完成后,将u(t)传送给组态王,再由组态王传给PLC,PLC将u(t)转换为0-10V电压信号,送至变频器MM440,调节输出频率,进而调节鼓风机转速,使溶解氧浓度根据可降解有机物浓度高低,分别控制在4mg/L和2.2mg/L。
本发明与现有技术相比具有的突出的实质性优点为:
(1)建立的数学模型满足污水处理曝气环节的工艺要求。
(2)设计的控制器能够使曝气池中溶解氧浓度根据曝气过程中的需求实现切换。
(3)控制精确,操作简单,控制系统的动态性能和稳态性能均能满足活性污泥污水处理的工艺要求。
新的建模方法从曝气环节的整体工艺出发,建立的数学模型符合曝气环节生化反应的整体流程;滑模控制方法的应用,提高了控制系统对不确定参数、外界扰动及滞后性的鲁棒性能。曝气环节中溶解氧浓度的实时控制,提高了活 性污泥法污水处理曝气环节的去污能力,提升了处理后的水质,并且溶解氧浓度的实时控制避免了能源的浪费,更加节能。
附图说明
图1为本发明中活性污泥污水处理过程基本流程图;
图2为活性污泥曝气环节滑模控制过程示意图;
图3为滑模控制曲线;
图4为PID控制曲线。
具体实施方式
下面结合附图对本发明作进一步的说明。
一.采用活性污泥污水处理曝气环节的工艺流程,即废水首先进入初次沉淀池沉淀,然后在曝气池中进行生物降解,然后在二次沉淀池进行固液分离,分离出上清水排出,底部污泥部分回流污泥至曝气池,剩余活性污泥处理后排出。
本发明为解决现在城市生活污水处理方法及控制系统的开发方面均比较落后的问题,提供一种活性污泥污水处理曝气环节的建模、滑模控制器及滑模控制方法。本发明方法简单易行,具有很好的动态品质和对外部干扰具有很强的适应能力,提高了控制系统对不确定参数、外界扰动及滞后性的鲁棒性能。
二.本发明的技术方案为下述内容:
1.一种活性污泥污水处理曝气环节的建模,包括以下步骤:
a.采用活性污泥污水处理曝气环节的工艺流程,即废水首先进入初次沉淀池沉淀,然后在曝气池中进行生物降解,然后在二次沉淀池进行固液分离,分离出上清水排出,底部污泥部分回流污泥至曝气池,剩余活性污泥处理后排出,如图1为本发明中活性污泥污水处理过程基本流程图。
b.根据活性污泥1号模型(ASM1)和微生物的生长动力学方程、物料平衡定律,将步骤a中工艺流程建立活性污泥法污水处理曝气环节的数学模型,利用 滑模控制器对其进行智能控制,其数学模型为:
其中,X是微生物浓度(grVSS/m3),S是可降解有机物浓度(grBOD5/m3),O是溶解氧浓度(grO2/m3),Xi是进水微生物浓度(grVSS/m3),Si是进水底物的浓度(grDOB5/m3),表示的是微生物的最大比增长速率(2/day),ks是半速度常数(grBOD5/m3),ko是半饱和系数(2grO2/m3),kd是子衰减系数(0.05/day),Q是进水流量(m3/day),V表示的是曝气池的体积(m3),Qw表示排出污泥的流量(m3/day),C是二次沉淀池中污泥的浓缩系数,f是联系有机物与需氧量的一个因子(0.67grBOD5/COD),fx是消耗因子(1.42grCOD/grVSS),u是曝气量(grO2/m3day),YNH是观察生长系数。
c.在废水进入曝气池初始阶段,可降解有机物浓度高于半速度常数,溶解氧浓度高于半饱和系数时,对数学模型进行推导简化得到:
其中,δ是对溶解氧设置的一个充量系数。
d.随着微生物生化反应的进行,可降解有机物逐渐被消耗,浓度在不断地 降低,当可降解有机物浓度低于半速度常数,溶解氧浓度低于半饱和系数时,数学模型简化为:
e.以底物浓度的大小将活性污泥污水处理系统分为两个阶段,即可降解有机物浓度高于半速度常数和可降解有机物浓度低于半速度常数,改进的数学模型变换成状态空间表达式的形式:
其中,
B=[0 0 1]T,C=[0 0 1],D=0
从ASM1模型中能够非常明显的观察出,利用活性污泥处理污水的过程是一个复杂的动态生物反应过程,在这一过程中,污水处理系统具有参数的不确定性、非线性、时变性以及各个因素之间的强耦合等特性。
本发明将非线性系统变换成两个线性系统之间的切换,从而将复杂的非线性过程变成一个简单的切换系统。
由于污水处理曝气过程是通过微生物的作用,将有机污染物转变成无害的 气体产物(如CO2,NO2,N2)、液体产物(如水)以及富含有机物的固体产物(生物污泥),其中生物污泥在沉淀池沉淀,从净化后的污水中除去。从控制的角度来看污水处处理曝气过程具有如下特性:一是非线性,污水处理过程中的各变量间为高度非线性,并且由于涉及到机理复杂的生化反应,这些非线性无法用传统的方法来描述;二是不确定性,污水处理过程的不确定性包括由于进水水量水质剧烈波动而引起的系统不确定性和传感器检测信息的不确定性;三是大滞后,进水水量水质的波动要滞后较长时间才与生化池的混合液悬浮固体浓度有关,生化池的混合液悬浮固体浓度、溶解氧浓度变化滞后一段时间后才影响到出水的悬浮物、生化需氧量、溶解氧等质量指标。并且滞后时间与进水量有关,因此滞后时间也是不确定的,另一方面传感器的检测也存在滞后。考虑到实际活性污泥污水处理的系统具有非线性、不确定性、大滞后的特点,为了不失一般性,故在上述状态方程中加入不确定、时滞以及外界干扰,上述状态方程变为:
其中,ΔAσ、ΔA1σ为系统的不确定性,d为常时滞,w为外界干扰。该方程组更适用于后边采用的滑模控制方法。
2.一种用于上述的活性污泥污水处理曝气环节的建模的滑模控制器,是将建立的数学模型应用滑模控制算法设计活性污泥污水处理曝气环节的滑模控制器,所述滑模控制器是根据系统所期望的动态特性来设计系统的切换超平面(滑模面),通过滑模控制器使系统状态从超平面之外向滑模面收束,系统一旦到达滑模面,控制作用将保证系统沿滑模面到达系统原点,由于系统的特性和参数只取决于设计的滑模面而与外界干扰没有关系,所以滑模变结构控制具有很强 的鲁棒性。从而实现对切换数学模型的精确控制,使稳定时间缩短,控制精度提高,抗干扰能力和鲁棒性提高;
所述滑模控制器包括两部分,其中一部分是等效控制器ueq,用于保证从任意状态出发的系统始终满足滑模到达条件;另一部分是校正控制器uN,用于抵消系统中存在的各种不确定性和外界干扰。
所述滑模控制器的设计,具体步骤为:
所述滑模控制的滑模面的选取s(t)=BTPx(t);
其中B是权利要求2中所建立的数学模型中的输入矩阵,P是正定对称矩阵,x(t)是权利要求2中所建立的数学模型中的状态向量。
所述滑模控制器的选取:
其中u(t)是滑模控制器,其作用是控制活性污泥污水处理曝气环节溶解氧浓度,uN=-(λ||x(t)||+k+μ)sign(s(t)),λ,k,μ是正常数,sign(s(t))是关于滑模面函数的符号函数。
所述滑模控制器的设计,具体步骤为:
第一步,滑模面的选取,
第二步,滑模控制器的选取
第三步,滑模控制器控制效果的稳定性证明。
1)针对第一步,滑模面选取,选取的依据是,保证存在某一区域,使得状态点趋近该区域时,就会被吸引到该区域内运动,该区域就是滑模面。按照滑模面必须能够吸引状态点的要求,滑模面必须满足以下条件,此条件称为滑模到达条件。
其中s代表滑模面。当s从大于0方向趋近于0时,满足的极限小于0,当s从小于0方向趋近于0时,满足的极限大于0,即。
此处针在加入不确定、时滞以及外界干扰因素时所建立的数学模型,选取滑模面s(t)=BTPx(t),保证滑模面以外的状态点都将在有限时间内到达滑模面。
其中B是在加入不确定、时滞以及外界干扰因素时所建立的数学模型中的输入矩阵,P是正定对称矩阵,x(t)是在加入不确定、时滞以及外界干扰因素时所建立的数学模型中的状态向量。
2)针对第二步,滑模控制器的选取,选取的依据是,保证系统的状态点一旦到达滑模面,此状态点将只能沿滑模面运动,在有限时间内达到期望值。滑模控制器应该包括两部分,其中一部分是等效控制器ueq,作用是保证从任意状态出发的系统始终满足滑模到达条件;另一部分是校正控制器uN,其作用是抵消系统中存在的各种不确定性和外界干扰。
针对上一步中所选取的滑模面,滑模控制器选取如下,
其中u(t)是滑模控制器,其作用是控制活性污泥污水处理曝气环节溶解氧浓度,uN=-(λ||x(t)||+k+μ)sign(s(t)),λ,k,μ是正常数,sign(s(t))是关于滑模面函数的符号函数。
3)针对第三步,滑模控制器控制效果的稳定性证明。选取合适的李雅普诺夫函数V(t),使其满足的条件。证明如下:
其中Q是适当维数的常数矩阵。
将选取的李雅普诺夫函数对时间t求导,可以得到
将权利要求2中所建立数学模型中的代入上式,得到
从上式结果中可以看出,由于则此控制器u(t)能够保证活性污泥污水处理曝气环节溶解氧浓度渐近稳定到期望值。从而实现一种活性污泥污水处理曝气环节的建模、滑模控制器控制。
三.一种使用上述述滑模控制器对活性污泥污水处理曝气环节进行控制的方法,包括以下步骤:
步骤一,根据ASM1模型和微生物的生长动力学方程、物料平衡定律化简ASM1模型,得到化简数学模型;
步骤二,根据可降解有机物浓度的高低,将步骤一ASM1化简数学模型分为两个子系统,即可降解有机物浓度高于半速度常数子系统和可降解有机物浓度低于半速度常数子系统;
步骤三,利用切换系统理论将步骤二两个子系统结合为统一的数学模型, 不考虑不确定性、时滞以及外界干扰因素,得到简化数学模型;
步骤四,根据污水处理实际情况,对步骤三简化数学模型,加入不确定性、时滞以及外界干扰因素,得到切换数学模型;
步骤五,根据步骤四切换数学模型,设计滑模控制器,编写MATLAB滑模控制器子程序;
步骤六,开启上位机(工控机)和下位机(PLC)控制装置,利用组态王DDE通讯,定时在组态王中调用MATLAB滑模控制器子程序,实现滑模控制,实时在线调整溶解氧浓度;
步骤七,循环执行步骤六,即可得到溶解氧浓度实时控制曲线。
所述的一种对活性污泥污水处理曝气环节进行滑模控制的方法,包括以下具体步骤:
滑模控制采用上位机与下位机相结合的控制方式,所述上位机为安装有组态王(简称组态王)和MATLAB软件(简称MATLAB)的工业计算机,所述的下位机采用西门子公司的S7-315-2DP PLC,在本曝气环节自动控制过程中,溶氧仪采集的溶解氧浓度信号送给PLC,PLC以变量形式将溶解氧浓度信号送至上位机组态王,上位机组态王利用DDE通讯将其传给MATLAB,MATLAB运行根据所设计的滑模控制器编写的滑模控制程序,计算所需控制量u(t),计算完成后,将u(t)传送给组态王,再由组态王传给PLC,PLC将u(t)转换为0-10V电压信号,送至变频器MM440,调节输出频率,进而调节鼓风机转速,使溶解氧浓度根据可降解有机物浓度高低,分别控制在4mg/L和2.2mg/L。
四、本实施例为污水处理厂的日处理能力为3万吨,其中动力鼓风机房采用2台200KW的罗茨风机对曝气池进行曝气。在本曝气环节自动控制系统中,现场采集的溶解氧浓度信号传送给下位机PLC,上位机通过组态王以变量形式读取PLC中的溶解氧数据,然后利用DDE通讯将其传给MATLAB,由MATLAB运行根据本发明中控制方法所编写的滑模控制计算程序,计算完成后,将输出值回传 给组态王,再由组态王传给PLC,PLC通过模拟量输出模块SM332输出0-10V电压信号,调节变频器MM430的输出频率,进而调节鼓风机转速,实现对溶解氧浓度的控制,图2为本发明活性污泥曝气环节滑模控制过程示意图。
新的建模方法从曝气环节的整体工艺出发,建立的数学模型符合曝气环节生化反应的整体流程;滑模控制方法的应用,提高了控制系统对不确定参数、外界扰动及滞后性的鲁棒性能。曝气环节中溶解氧浓度的实时控制,提高了活性污泥法污水处理曝气环节的去污能力,提升了处理后的水质,并且溶解氧浓度的实时控制避免了能源的浪费,更加节能。
组态王的DDE对话的内容是通过应用程序名、主题名和项目三个标识名来约定的。组态王运行系统的程序名是“VIEW”,主题名是“GNAME”,项目名则是在定义DDE变量时定义的项目名称,是具体的数据对象。MATLAB也支持DDE,当MATLAB以客户身份建立DDE通信时,需要使用MATLAB中的DDE客户端模块所提供的固有函数,这些函数包括:初始化函数,定义程序名和主题名的初始值;建立链接,本发明中采用热链接;请求数据;发送数据;释放链接;链接终止。
本发明中,滑模控制算法的编程实现,分为DDE初始化部分和控制程序两部分。在MATLAB中创建新的m文件,编写DDE初始化函数以及控制函数,针对本控制系统,初始化DDE程序按照DDE协议,定义上述6个函数,将组态王中的变量“S7300.DB.54”传送给MATLAB中的变量x3(i),此变量代表溶解氧浓度值。
DDE初始化成功后,则可运行调试主程序,现将程序运行过程描述如下:
1)定义函数,将溶解氧浓度检测值对应MATLAB中的变量。
2)请求数据,将组态王中溶解氧浓度变量传给MATLAB中变量。
3)调用滑模控制器子程序,实现滑模控制算法。
4)发送数据,将滑模控制器子程序的输出值回传给组态王中相应的变量,此变量用于控制变频器频率,进而控制鼓风机。
5)用ddeterm命令结束DDE通讯。
6)循环执行上述步骤,实现溶解氧浓度实时控制。
通过逐一完成以上所说的步骤,便可以实现活性污泥污水处理滑模控制,将得出的实时趋势曲线,与PID控制实时趋势曲线相比较,可明显看出滑模控制曲线波动小,超调量也小,响应时间快,能够更快速的达到设定值并保持稳定。如图3,图4。
本发明与现有技术相比具有的突出的实质性优点为:
1)建立的数学模型满足污水处理曝气环节的工艺要求
2)设计的控制器能够使曝气池中溶解氧浓度根据曝气过程中的需求实现切换。
3)控制精确,操作简单,控制系统的动态性能和稳态性能均能满足活性污泥污水处理的工艺要求。
Claims (5)
1.一种活性污泥污水处理曝气环节的建模,其特征是包括以下步骤:
a.采用活性污泥污水处理曝气环节的工艺流程,即废水首先进入初次沉淀池沉淀,然后在曝气池中进行生物降解,然后在二次沉淀池进行固液分离,分离出上清水排出,底部污泥部分回流污泥至曝气池,剩余活性污泥处理后排出;
b.根据活性污泥1号模型(ASM1)和微生物的生长动力学方程、物料平衡定律,将步骤a中工艺流程建立活性污泥法污水处理曝气环节的数学模型,利用滑模控制器对其进行智能控制,其数学模型为:
其中,X是微生物浓度(gr VSS/m3),S是可降解有机物浓度(gr BOD5/m3),O是溶解氧浓度(gr O2/m3),Xi是进水微生物浓度(gr VSS/m3),Si是进水底物的浓度(gr DOB5/m3),表示的是微生物的最大比增长速率(2/day),ks是半速度常数(gr BOD5/m3),ko是半饱和系数(2gr O2/m3),kd是子衰减系数(0.05/day),Q是进水流量(m3/day),V表示的是曝气池的体积(m3),Qw表示排出污泥的流量(m3/day),C是二次沉淀池中污泥的浓缩系数,f是联系有机物与需氧量的一个因子(0.67gr BOD5/COD),fx是消耗因子(1.42gr COD/grVSS),u是曝气量(grO2/m3day),YNH是观察生长系数。
c.在废水进入曝气池初始阶段,可降解有机物浓度高于半速度常数,溶解氧浓度高于半饱和系数时,对数学模型进行推导简化得到:
其中,δ是对溶解氧设置的一个充量系数;
d.随着微生物生化反应的进行,可降解有机物逐渐被消耗,浓度在不断地降低,当可降解有机物浓度低于半速度常数,溶解氧浓度低于半饱和系数时,数学模型简化为:
e.以底物浓度的大小将活性污泥污水处理系统分为两个阶段,即可降解有机物浓度高于半速度常数和可降解有机物浓度低于半速度常数,改进的数学模型变换成状态空间表达式的形式:
其中,
B=[0 0 1]T,C=[0 0 1],D=0 。
2.如权利要求1所述的活性污泥污水处理曝气环节的建模,其特征是,在加入不确定、时滞以及外界干扰因素时,数学模型为:
其中,ΔAσ、ΔAdσ为系统的不确定性,d为常时滞,w为外界干扰。
3.一种用于权利要求1所述的活性污泥污水处理曝气环节的建模的滑模控制器,其特征是,建立的数学模型应用滑模控制算法
将建立的数学模型应用滑模控制算法设计活性污泥污水处理曝气环节的滑模控制器,所述滑模控制器是根据系统所期望的动态特性来设计系统的切换超平面(滑模面),通过滑模控制器使系统状态从超平面之外向滑模面收敛,系统一旦到达滑模面,控制作用将保证系统沿滑模面到达系统原点,从而实现对切换数学模型的精确控制,使稳定时间缩短,控制精度提高,抗干扰能力和鲁棒性提高;
所述滑模控制器包括两部分,其中一部分是等效控制器ueq,用于保证从任意状态出发的系统始终满足滑模到达条件;另一部分是校正控制器uN,用于抵消系统中存在的各种不确定性和外界干扰。
所述滑模控制器的设计,具体步骤为:
所述滑模控制的滑模面的选取s(t)=BTPx(t);
其中B是权利要求2中所建立的数学模型中的输入矩阵,P是正定对称矩阵,x(t)是权利要求2中所建立的数学模型中的状态向量。
所述滑模控制器的选取:
其中u(t)是滑模控制器,其作用是控制活性污泥污水处理曝气环节溶解氧浓度,uN=-(λ||x(t)||+k+μ)sign(s(t)),λ,k,μ是正常数,sign(s(t))是关于滑模面函数的符号函数。
4.一种使用权利要求3所述滑模控制器对活性污泥污水处理曝气环节进行滑模控制的方法,其特征是包括以下步骤:
步骤一,根据ASM1模型和微生物的生长动力学方程、物料平衡定律化简ASM1模型,得到化简数学模型;
步骤二,根据可降解有机物浓度的高低,将步骤一ASM1化简数学模型分为两个子系统,即可降解有机物浓度高于半速度常数子系统和可降解有机物浓度低于半速度常数子系统;
步骤三,利用切换系统理论将步骤二两个子系统结合为统一的数学模型,不考虑不确定性、时滞以及外界干扰因素,得到简化数学模型;
步骤四,根据污水处理实际情况,对步骤三简化数学模型,加入不确定性、时滞以及外界干扰因素,得到切换数学模型;
步骤五,根据步骤四切换数学模型,设计滑模控制器,编写MATLAB滑模控制器子程序;
步骤六,开启上位机和下位机的控制装置,利用组态王DDE通讯,定时在组态王中调用MATLAB滑模控制器子程序,实现滑模控制,实时在线调整溶解氧浓度;
步骤七,循环执行步骤六,即可得到溶解氧浓度实时控制曲线。
5.如权利要求4所述的一种对活性污泥污水处理曝气环节进行滑模控制的 方法,其特征是包括以下具体步骤:
滑模控制采用上位机与下位机相结合的控制方式,所述上位机为安装有组态王(简称组态王)和MATLAB软件(简称MATLAB)的工业计算机,所述的下位机采用西门子公司的S7-315-2DP PLC,在本曝气环节自动控制过程中,溶氧仪采集的溶解氧浓度信号送给PLC,PLC以变量形式将溶解氧浓度信号送至上位机组态王,上位机组态王利用DDE通讯将其传给MATLAB,MATLAB运行根据所设计的滑模控制器编写的滑模控制程序,计算所需u(t)控制量,计算完成后,将u(t)传送给组态王,再由组态王传给PLC,PLC将u(t)转换为0-10V电压信号,送至变频器MM440,调节输出频率,进而调节鼓风机转速,使溶解氧浓度根据可降解有机物浓度高低,分别控制在4mg/L和2.2mg/L。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410231030.3A CN104111618A (zh) | 2014-05-28 | 2014-05-28 | 一种活性污泥污水处理曝气环节的建模、滑模控制器及其控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410231030.3A CN104111618A (zh) | 2014-05-28 | 2014-05-28 | 一种活性污泥污水处理曝气环节的建模、滑模控制器及其控制方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN104111618A true CN104111618A (zh) | 2014-10-22 |
Family
ID=51708449
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410231030.3A Pending CN104111618A (zh) | 2014-05-28 | 2014-05-28 | 一种活性污泥污水处理曝气环节的建模、滑模控制器及其控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104111618A (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105608976A (zh) * | 2015-12-24 | 2016-05-25 | 中国信息安全测评中心 | 一种虚实结合污水处理厂测试床及其构建方法 |
CN104238527B (zh) * | 2014-09-24 | 2017-02-15 | 天津创业环保集团股份有限公司 | 污水处理厂曝气总量的精确控制方法 |
CN106568913A (zh) * | 2016-11-04 | 2017-04-19 | 济南大学 | 一种基于不变观测器的活性污泥法污水处理进度估计方法 |
CN109534495A (zh) * | 2019-01-21 | 2019-03-29 | 清华大学深圳研究生院 | 一种基于微生物基因和活性酶的污水处理方法及系统 |
CN112327615A (zh) * | 2020-09-29 | 2021-02-05 | 北京工业大学 | 基于无抖振模糊滑模的污水处理脱氮过程控制方法 |
CN112987567A (zh) * | 2021-02-09 | 2021-06-18 | 河北科技大学 | 非线性系统的固定时间自适应神经网络滑模控制方法 |
CN112979058A (zh) * | 2021-02-03 | 2021-06-18 | 齐鲁工业大学 | 基于微生物燃料电池的啤酒厂废水处理装置及方法 |
CN113552837A (zh) * | 2021-07-26 | 2021-10-26 | 中冶华天南京工程技术有限公司 | 一种用于智能化料场的流程控制方法 |
CN114019797A (zh) * | 2021-10-29 | 2022-02-08 | 湖北工业大学 | 无时滞名义模型下污水处理设备臭氧浓度滑模抗扰控制方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100059437A1 (en) * | 2008-09-11 | 2010-03-11 | University Of Central Florida Research Foundation, Inc. | Passive Underground Drainfield For Septic Tank Nutrient Removal Using Functionalized Green Filtration Media |
CN201662710U (zh) * | 2010-08-19 | 2010-12-01 | 中持(北京)环保发展有限公司 | 用于污水处理系统的双参数曝气控制系统 |
CN102156432A (zh) * | 2011-02-22 | 2011-08-17 | 上海市城市建设设计研究院 | 一种实时控制生化反应池中曝气量的方法 |
CN102690015A (zh) * | 2011-03-24 | 2012-09-26 | 中国科学院沈阳自动化研究所 | 动态多级缺氧/好氧污水处理方法 |
CN103809436A (zh) * | 2012-11-06 | 2014-05-21 | 西安元朔科技有限公司 | 活性污泥法污水处理过程智能建模方法 |
-
2014
- 2014-05-28 CN CN201410231030.3A patent/CN104111618A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100059437A1 (en) * | 2008-09-11 | 2010-03-11 | University Of Central Florida Research Foundation, Inc. | Passive Underground Drainfield For Septic Tank Nutrient Removal Using Functionalized Green Filtration Media |
CN201662710U (zh) * | 2010-08-19 | 2010-12-01 | 中持(北京)环保发展有限公司 | 用于污水处理系统的双参数曝气控制系统 |
CN102156432A (zh) * | 2011-02-22 | 2011-08-17 | 上海市城市建设设计研究院 | 一种实时控制生化反应池中曝气量的方法 |
CN102690015A (zh) * | 2011-03-24 | 2012-09-26 | 中国科学院沈阳自动化研究所 | 动态多级缺氧/好氧污水处理方法 |
CN103809436A (zh) * | 2012-11-06 | 2014-05-21 | 西安元朔科技有限公司 | 活性污泥法污水处理过程智能建模方法 |
Non-Patent Citations (2)
Title |
---|
吕璇: "鲁棒控制方法在污水处理中的应用研究", 《中国优秀硕士学位论文全文数据库 工程科技I辑》 * |
张巧玲: "污水处理过程溶解氧浓度的滑模控制研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104238527B (zh) * | 2014-09-24 | 2017-02-15 | 天津创业环保集团股份有限公司 | 污水处理厂曝气总量的精确控制方法 |
CN105608976B (zh) * | 2015-12-24 | 2018-12-14 | 中国信息安全测评中心 | 一种虚实结合污水处理厂测试床及其构建方法 |
CN105608976A (zh) * | 2015-12-24 | 2016-05-25 | 中国信息安全测评中心 | 一种虚实结合污水处理厂测试床及其构建方法 |
CN106568913A (zh) * | 2016-11-04 | 2017-04-19 | 济南大学 | 一种基于不变观测器的活性污泥法污水处理进度估计方法 |
CN109534495A (zh) * | 2019-01-21 | 2019-03-29 | 清华大学深圳研究生院 | 一种基于微生物基因和活性酶的污水处理方法及系统 |
CN109534495B (zh) * | 2019-01-21 | 2021-08-24 | 清华大学深圳研究生院 | 一种基于微生物基因和活性酶的污水处理方法及系统 |
CN112327615B (zh) * | 2020-09-29 | 2022-03-15 | 北京工业大学 | 基于无抖振模糊滑模的污水处理脱氮过程控制方法 |
CN112327615A (zh) * | 2020-09-29 | 2021-02-05 | 北京工业大学 | 基于无抖振模糊滑模的污水处理脱氮过程控制方法 |
CN112979058A (zh) * | 2021-02-03 | 2021-06-18 | 齐鲁工业大学 | 基于微生物燃料电池的啤酒厂废水处理装置及方法 |
CN112987567A (zh) * | 2021-02-09 | 2021-06-18 | 河北科技大学 | 非线性系统的固定时间自适应神经网络滑模控制方法 |
CN112987567B (zh) * | 2021-02-09 | 2022-07-22 | 河北科技大学 | 非线性系统的固定时间自适应神经网络滑模控制方法 |
CN113552837A (zh) * | 2021-07-26 | 2021-10-26 | 中冶华天南京工程技术有限公司 | 一种用于智能化料场的流程控制方法 |
CN113552837B (zh) * | 2021-07-26 | 2022-11-01 | 中冶华天南京工程技术有限公司 | 一种用于智能化料场的流程控制方法 |
CN114019797A (zh) * | 2021-10-29 | 2022-02-08 | 湖北工业大学 | 无时滞名义模型下污水处理设备臭氧浓度滑模抗扰控制方法 |
CN114019797B (zh) * | 2021-10-29 | 2023-12-22 | 湖北工业大学 | 无时滞名义模型下污水处理设备臭氧浓度滑模抗扰控制方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104111618A (zh) | 一种活性污泥污水处理曝气环节的建模、滑模控制器及其控制方法 | |
CN103499982B (zh) | 一种污水处理过程自组织控制方法 | |
CN104865979A (zh) | 一种污水处理过程自适应广义预测控制方法及系统 | |
Lee et al. | Evaluation of aeration energy saving in two modified activated sludge processes | |
Piotrowski | Two‐Level Multivariable Control System of Dissolved Oxygen Tracking and Aeration System for Activated Sludge Processes | |
CN104478076A (zh) | 一种数字化微氧节能曝气装置及曝气工艺 | |
Åmand et al. | The optimal dissolved oxygen profile in a nitrifying activated sludge process–comparisons with ammonium feedback control | |
Hirsch et al. | Two-step model based adaptive controller for Dissolved Oxygen control in Sequencing Wastewater Batch Reactor | |
Simon-Várhelyi et al. | Dairy wastewater processing and automatic control for waste recovery at the municipal wastewater treatment plant based on modelling investigations | |
CN107728479A (zh) | 一种基于rbf神经网络的生物除磷精确控制方法 | |
CN105138716B (zh) | 硝化及亚硝化过程的运行优化方法 | |
Zhang et al. | Finite-Time Extended State Observer-based PI Control for Dissolved Oxygen | |
Swarnalatha | Simulation and Analysis of Ammonium-Based Aeration Control Strategies to Enhance Efficiency in Wastewater Treatment Plant | |
CN113087030B (zh) | 一种水资源双循环反馈控制系统及方法 | |
CN204369641U (zh) | 一种数字化微氧节能曝气装置 | |
CN102830720A (zh) | 基于模糊融合控制技术的污水溶解氧浓度控制方法 | |
Xu et al. | Application of fuzzy control on wastewater treatment plant for P-removal | |
Carp et al. | Robust control of an activated sludge wastewater treatment process | |
Nascu et al. | Multilevel predictive control system for an activated sludge wastewater treatment process | |
Zawadzki et al. | Nonlinear fuzzy control of the dissolved oxygen in activated sludge processes | |
Costello et al. | Control of anaerobic digesters using Generic Model Control | |
Angulo et al. | Control of an anaerobic upflow fixed bed bioreactor | |
EL-Zemity et al. | Wastewater treatment model with smart irrigation utilizing PID control | |
Islami et al. | Stability Analysis of Biological Wastewater Treatment in Facultative Stabilization Ponds: Mathematical Model | |
Ye et al. | Dissolved oxygen generic model control of wastewater treatment process based on immune optimization least squares support vector machine |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20141022 |
|
RJ01 | Rejection of invention patent application after publication |