污水处理控制方法

　　申请日2017.06.01

　　公开(公告)日2017.09.15

　　IPC分类号G05B13/02

　　摘要

　　本发明涉及一种差分自适应万有引力搜索算法和一种基于该算法的污水处理控制方法。一种差分自适应万有引力搜索算法，包括如下步骤：S1、初始化种群;S2、计算种群适应度值;S3、计算种群的惯性质量M及加速度a;S4、更新个体位置;S5、比对是否达到迭代次数，当达到迭代次数时，输出最优值;当未达到迭代次数时，重复所述步骤S2、S3和S4。本发明的有益效果是：在加快收敛速度，收敛速度的跟步长有关，现在最好的方法通过实时更新个体参数G0和加速度α，有效控制算法收敛速度，同时提高算法精度。

　　摘要附图

　　权利要求书

　　1.一种差分自适应万有引力搜索算法，其特征在于，包括如下步骤：

　　S1、初始化种群;

　　S2、计算种群适应度值;

　　S3、计算种群的惯性质量M及加速度a;

　　S4、更新个体位置;

　　S5、比对是否达到迭代次数，当达到迭代次数时，输出最优值;当未达到迭代次数时，重复所述步骤S2、S3和S4。

　　2.根据权利要求1所述一种差分自适应万有引力搜索算法，其特征在于，所述步骤S4包括将万有引力搜索算法的公式带入更新公式中;

　　所述万有引力搜索算法的公式包括：

　　粒子间的引力为：

　　 F = G M 1 M 2 R 2

　　式中，G表示引力常量，M1、M2分别表示两种粒子的惯性质量，R表示两种粒子中间的欧式距离，F表示引力;

　　依据牛顿第二定律，当某个粒子受到其他粒子作用力时，随即会产生F大小的合力，同时会产生向这个方向运动的加速度a，即

　　 a = F m

　　引力常数G跟随着宇宙实际年龄变化的趋势为

　　 G ( t ) = G ( t 0 ) ( t 0 t ) β , β ≤ 1

　　式中，G(t)表示在t时间时的引力常数G的数值，G(t0)表示在第一宇宙子间隔时刻t0的数值;

　　所述更新公式包括：

　　有吸引力就会产生相应的加速度，则物质i在第k维上的加速度为

　　 a i k ( t ) = F i k ( t ) M i i ( t )

　　每次迭代更新的公式为

　　 v i k ( t + 1 ) = rand i × v i k ( t ) + a i k ( t ) x i k ( t + 1 ) = x i k ( t ) + v i k ( t + 1 ) .

　　3.一种污水处理控制方法，其特征在于，包括如下步骤:

　　S100、设定该问题的搜索空间维度为2，设置如下的空间上下限范围：

　　 S NH 3 < 8 C O D < 60 0 < ρ D O < 4 0 < Q w < 2000 1

　　其中，是氨氮含量;COD是化学需氧量;ρDO溶氧量;Qw是污泥排放量;

　　S200、种群初始化;

　　S300、初始化全局最优值gbest;

　　S400、计算适应值和加速度;

　　S500、更新速度与位置;

　　S600、更新全局最优粒子与个体最优粒子，得到全局最优解pbesti和gbest;

　　S700、比对是否达到迭代次数，当达到迭代次数时，输出最优解，当未达到迭代次数时，重复所述步骤S300、S400、S500和S600;

　　S800、通过输出的最优值，找到函数值最小时方程所对应的解，即溶解氧浓度ρDO优化设定值和污泥排放量Qw的优化设定值。

　　4.根据权利要求3所述一种污水处理控制方法，其特征在于，所述步骤S200包括：设置万有引力场的粒子数目100，最大迭代步数M=200;根据电脑随机初始化种群中各个粒子的位置X、速度Y以及速度的上下限固定范围[Vmin,Vmax]，每个粒子个体历史最优位置pbesti赋值：对于每个粒子，初始时都有相应的赋值pbesti。

　　5.根据权利要求3所述一种污水处理控制方法，其特征在于，所述步骤S300包括：在规定的目标要求下，依次算出各个初始粒子X的目标函数值f1,f2，电脑随机取其中某个粒子做为初始的gbest。

　　6.根据权利要求3所述一种污水处理控制方法，其特征在于，所述步骤S400包括：

　　S401、计算其目标函数值，而后根据线性权重的方法，把原先设定好的目标适应度值统一成一个适应度值：

　　 f i t n e s s = f 1 × ( 1 - t M ) + f 2 × t M

　　式中的t是当前的迭代次数;

　　S402、通过下式计算加速度：

　　 F i j k ( t ) = G ( t ) M p t ( t ) M a t ( t ) R i j ( t ) + ϵ ( x j k ( t ) - x i k ( t ) )

　　 a i k ( t ) = F i k ( t ) M i i ( t ) .

　　7.根据权利要求6所述一种污水处理控制方法，其特征在于，所述步骤S401中目标函数方程为：

　　 J = a × m i n { &Integral; t t + T P d t + 1 T &Integral; t t + T [ c 1 ( C O D - 60 ) V + c 2 ( S NH 3 - 8 ) V ] d t }

　　目标函数中T为优化周期，P为鼓风机的功率，a为每度电的费用，V是污水处理厂体积，c1=0.455为每千克去除化学需氧量COD的电耗，c2=9.1为每千克去除氨氮含量NH3-H的电耗，系数c1V,c2V可以记为电费的奖惩因子。

　　8.根据权利要求7所述一种污水处理控制方法，其特征在于，所述步骤S401中的目标函数方程可以简化为：

　　 J = a × min { &Integral; t t + T P d t + 1 T &Integral; t t + T [ c 1 ( C O D - 60 ) V + c 2 ( S NH 3 - 8 ) V ] d t } = a × min { &Integral; t t + T P d t + 1 T &Integral; t t + T [ 0.455 × ( C O D - 60 ) × 6.1 + 9.1 ( S NH 3 - 8 ) × 6.1 ] d t } = a × min { &Integral; t t + T P d t + 1 T &Integral; t t + T [ 2.78 × ( C O D - 60 ) + 55.5 × ( S NH 3 - 8 ) ] d t } .

　　9.根据权利要求3所述一种污水处理控制方法，其特征在于，所述步骤S500包括，通过下式更新速度与位置：

　　 X i t + 1 = X i t + r a n d × V i t + G 0 × e - a t / T × Σ j = 1 N ( M j t R i j ( r a n d × ( X j t - X i t ) ) )

　　其中，是当前位置;前一时刻位置;Vit前一时刻速度;Rij粒子间的欧式距离;是惯性质量。

　　说明书

　　一种引力搜索算法及基于该算法的污水处理控制方法

　　技术领域

　　本发明涉及一种引力搜索算法及基于该算法的污水处理控制方法。

　　背景技术

　　只有在污水处理工艺中引入优质的计算机技术、当代优化算法和控制理论才能提升污水处理的效率，并减少处理成本。万有引力搜索算法(Gravitational searchalgorithm，GSA)是由伊朗克曼大学的Esmat Rashedi等人于2009年所分析的一种新的启蒙式优化算法。此种GSA算法的实质是模仿物理学中的万有引力现象，并将其演算成随机寻找最优解的过程。可是GSA目前使用较少，对自适应的研究更是寥寥无几。

　　虽然GSA算法是全局收敛的算法，能够寻找到最优值，但是如果算法在收敛之前没有搜索到全局(或者局部)最优解，将导致过早收敛而找不到最优值，也就是说，GSA跟其他优化算法一样，容易陷入局部最优。

　　现在我国已经建设了很多的污水处理厂，只要目的是减少水污染对环境形成的消极影响，当然处理过程的自动控制方法很前沿，但至今仍没有实现最优控制。只有污水处理已经到达《城镇污水处理厂污染物排放标准》(GB18918-2002)，然后才可以思量控制的目标。从环境保护的角度研究，出水水质指标越小越好，而从运行成本研究，出水水质越好意味着消耗更多的资源和更多的运营费用。宁夏某污水处理厂使用先进的污水处理设施，采纳改良型SBR技术。

　　改良型SBR工艺还存在两个主要问题：一是污水处理中的水质参数检测滞后，采用的是离线人工检测方式，一般是人工取水，然后人工在实验室化验，得出出水水质指标COD(化学需氧量)、氨氮HN3-N等的化验结果，根据水质情况调整鼓风曝气盘。这样的检测方式需要耗时半小时，明显的对水质控制不够及时，效率低下。二是只有在水质达标的前提下，较少工艺的运行费用，才能够增强其经济效益。该工艺的运行费用主要包括人工费、污泥处理费、鼓风曝气费这三部分，其中人工费和污泥处理费一般根据工厂规模可以固定，而鼓风曝气费消耗大，且难控制。鼓风曝气量需要根据污泥排放量和溶解氧的浓度进行实时快速的调整，才能是曝气量最少，从而费用最优。

　　发明内容

　　本发明所要解决的技术问题是提供一种有效控制算法收敛速度，同时提高算法精度的差分自适应万有引力搜索算法。

　　本发明解决上述技术问题的技术方案如下：一种差分自适应万有引力搜索算法，包括如下步骤：

　　S1、初始化种群;

　　S2、计算种群适应度值;

　　S3、计算种群的惯性质量M及加速度a;

　　S4、更新个体位置;

　　S5、比对是否达到迭代次数，当达到迭代次数时，输出最优值;当未达到迭代次数时，重复所述步骤S2、S3和S4。

　　本发明的有益效果是：在加快收敛速度，收敛速度的跟步长有关，现在最好的方法通过实时更新个体参数G0和加速度α，有效控制算法收敛速度，同时提高算法精度。

　　在上述技术方案的基础上，本发明还可以做如下改进。

　　进一步，所述步骤S4包括将万有引力搜索算法的公式带入更新公式中;

　　所述万有引力搜索算法的公式包括：

　　粒子间的引力为：

　　式中，G表示引力常量，M1、M2分别表示两种粒子的惯性质量，R表示两种粒子中间的欧式距离，F表示引力;

　　依据牛顿第二定律，当某个粒子受到其他粒子作用力时，随即会产生F大小的合力，同时会产生向这个方向运动的加速度a，即

　　引力常数G跟随着宇宙实际年龄变化的趋势为

　　式中，G(t)表示在t时间时的引力常数G的数值，G(t0)表示在第一宇宙子间隔时刻t0的数值;

　　所述更新公式包括：

　　有吸引力就会产生相应的加速度，则物质i在第k维上的加速度为

　　每次迭代更新的公式为

　　本发明所要解决的技术问题是提供一种能预测出水水质的污水处理控制方法。

　　本发明解决上述技术问题的技术方案如下：一种污水处理控制方法，包括如下步骤:

　　S100、设定该问题的搜索空间维度为2，设置如下的空间上下限范围：

　　其中，是氨氮含量;COD是化学需氧量;ρDO溶氧量;Qw是污泥排放量。

　　S200、种群初始化;

　　S300、初始化全局最优值gbest;

　　S400、计算适应值和加速度;

　　S500、更新速度与位置;

　　S600、更新全局最优粒子与个体最优粒子，得到全局最优解pbesti和gbest;

　　S700、比对是否达到迭代次数，当达到迭代次数时，输出最优解，当未达到迭代次数时，重复所述步骤S300、S400、S500和S600;

　　S800、通过输出的最优值，找到函数值最小时方程所对应的解，即溶解氧浓度ρDO优化设定值和污泥排放量Qw的优化设定值。

　　本发明的有益效果是，通过输出的最优值，找到函数值最小时候，此时方程所对应的解，这就是得到了优AGSA算法寻求得到的溶解氧浓度ρDO优化设定值和污泥排放量Qw的优化设定值。该方法结合差分自适应引力搜索算法，使得鼓风曝气量根据进水口量和溶解氧的浓度进行实时快速的调整，使曝气量最少，从而费用最优。

　　所述步骤S200包括：设置万有引力场的粒子数目100，最大迭代步数M=200;根据电脑随机初始化种群中各个粒子的位置X、速度Y以及速度的上下限固定范围[Vmin,Vmax]，每个粒子个体历史最优位置pbesti赋值：对于每个粒子，初始时都有相应的赋值pbesti。

　　所述步骤S300包括：在规定的目标要求下，依次算出各个初始粒子X的目标函数值f1,f2，电脑随机取其中某个粒子做为初始的gbest。

　　所述步骤S400包括：

　　S401、计算其目标函数值，而后根据线性权重的方法，把原先设定好的目标适应度值统一成一个适应度值：

　　式中的t是当前的迭代次数;

　　S402、通过下式计算加速度：

　　所述步骤S401中目标函数为：

　　目标函数中T为优化周期，P为鼓风机的功率，a为每度电的费用，V是污水处理厂体积，c1=0.455为每千克去除COD(化学需氧量)的电耗，c2=9.1为每千克去除NH3-H(氨氮含量)的电耗，系数c1V,c2V可以记为电费的奖惩因子。

　　所述步骤S401中的目标方程可以简化为：

　　所述步骤S500包括，通过下式更新速度与位置：

　　其中，是当前位置;前一时刻位置;前一时刻速度;Rij粒子间的欧式距离;是惯性质量。

	搜索
	污水宝	污水处理风向标	投运污水厂电话查询


首页	行业新闻	招标项目	审批公告	拟建项目	中标公示	十三五规划	业主委托	污水处理技术	水质标准	政策法规