CN103472729A - 群智的农药废液焚烧炉有害物排放达标控制系统及方法 - Google Patents
群智的农药废液焚烧炉有害物排放达标控制系统及方法 Download PDFInfo
- Publication number
- CN103472729A CN103472729A CN2013104338648A CN201310433864A CN103472729A CN 103472729 A CN103472729 A CN 103472729A CN 2013104338648 A CN2013104338648 A CN 2013104338648A CN 201310433864 A CN201310433864 A CN 201310433864A CN 103472729 A CN103472729 A CN 103472729A
- Authority
- CN
- China
- Prior art keywords
- training sample
- cod
- particle
- iter
- value
- 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.)
- Granted
Links
Images
Landscapes
- Feedback Control In General (AREA)
Abstract
本发明公开了一种群智的农药废液焚烧炉有害物排放达标控制系统及方法。它包括焚烧炉、智能仪表、DCS系统、数据接口以及上位机;DCS系统包括控制站和数据库;用于测量易测变量的智能仪表与DCS系统连接,DCS系统通过数据接口与上位机连接。上位机首先对训练样本进行预处理和模糊化,获得新的输入矩阵,然后采用BP神经网络对新的输入矩阵建立回归模型,获得预测输出,再对预测输出进行反模糊化,获得系统的输出,最后引入粒子群算法对BP神经网络的参数进行优化;上位机还具有判别模型更新和结果显示的功能。本发明具有在线测量COD、有效监测COD是否超标、控制COD排放达标、抗噪声和自学习能力强、在线自动优化等优点。
Description
技术领域
本发明涉及农药生产领域,尤其涉及一种群智的农药废液焚烧炉有害物排放达标控制系统及方法。
背景技术
我国是农药生产和使用大国,农药生产企业已达4100家左右,其中原药生产企业为500多家,国家农业部统计数据显示2008年1~11月农药总产量达171.1万吨。我国农药品种结构的不合理性更加大了环境治理的难度。据不完全统计,全国农药工业每年排放的废水约为15亿吨。其中,处理达标的仅占已处理的1%。焚烧法是目前处理农药残液和废渣最有效、彻底、应用最普遍的方法。焚烧后废水的化学耗氧量(COD)是农药废液焚烧有害物排放的最重要指标,但是其无法在线测量,离线测量一次需要四五个小时,无法及时反映工况变化和指导实际生产。因此,在实际焚烧过程中,COD严重超标。
发明内容
为了克服已有的焚烧炉过程COD无法在线测量、COD严重超标的不足,本发明提供一种群智的农药废液焚烧炉有害物排放达标控制系统及方法,其具有在线测量COD、有效监测COD是否超标、控制COD排放达标、抗噪声和自学习能力强、在线自动优化等优点。
本发明解决其技术问题所采用的技术方案是:
群智的农药废液焚烧炉有害物排放达标控制系统,包括焚烧炉、智能仪表、DCS系统、数据接口以及上位机,所述的DCS系统包括控制站和数据库;所述现场智能仪表与DCS系统连接,所述DCS系统与上位机连接,所述的上位机包括:
数据预处理模块,用于将从DCS数据库输入的模型训练样本进行预处理,对训练样本中心化,即减去样本的平均值,然后对其进行标准化:
计算均值:
计算方差:
标准化:
其中,TXi为第i个训练样本,是从DCS数据库中采集的生产正常时的关键变量、化学耗氧量(COD)和相应的使COD排放达标时的操作变量的数据,N为训练样本数,为训练样本的均值,X为标准化后的训练样本。σx表示训练样本的标准差,σ2 x表示训练样本的方差。模糊方程模块,对从数据预处理模块传过来的标准化后的训练样本X,进行模糊化。设模糊方程系统中有c*个模糊群,模糊群k、j的中心分别为vk、vj,则标准化后的第i个训练样本Xi对于模糊群k的隶属度μik为:
式中,n为模糊分类过程中需要的分块矩阵指数,通常取作2,||·||为范数表达式。
使用以上隶属度值或者它的变形以获得新的输入矩阵,对于模糊群k,其输入矩阵变形为:
式中,wlk是隐层神经元l与输出层神经元之间的连接权,f(·)为神经元的输出函数,通常取为Sigmoid函数,表示为:
式中,h为输出层神经元的阈值,θ0为陡度参数,用以调节Sigmoid函数的陡度;最后,由反模糊方法中的重心法得到最后的模糊方程系统的输出:
粒子群算法优化模块,用于采用粒子群算法对模糊方程中BP神经网络局部方程的wlk进行优化,具体步骤如下:
①确定粒子群的优化参数为BP神经网络局部方程的wlk、粒子群个体数目popsize、最大循环寻优次数itermax、第p个粒子的初始位置rp、初始速度vp、局部最优值Lbestp以及整个粒子群的全局最优值Gbest。
②设定优化目标函数,将其转换为适应度,对每个局部模糊方程进行评价;通过相应的误差函数计算适应度函数,并认为误差大的粒子适应度小,粒子p的适应度函数表示为:
fp=1/(Ep+1) (10)
式中,Ep是模糊方程系统的误差函数,表示为:
③按照如下公式,循环更新每个粒子的速度和位置,
vp(iter+1)=ω×vp(iter)+m1a1(Lbestp-rp(iter))+m2a2(Gbest-rp(iter))
(12)
rp(iter+1)=rp(iter)+vp(iter+1) (13)
式中,vp表示更新粒子p的速度,rp表示更新粒子p的位置,Lbestp表示更新粒子p的个体最优值,Gbest即为对应于标准化后的第i个训练样本Xi的COD预报值和使COD排放达标的操作变量值,iter表示循环次数,ω是粒子群算法中的惯性权重,m1、m2是对应的加速系数,a1、a2是[0,1]之间的随机数;
④对于粒子p,如果新的适应度大于原来的个体最优值,更新粒子的个体最优值:
Lbestp=fp (14)
⑤如果粒子p的个体最优值Lbestp大于原来的粒子群全局最优值Gbest,更新原来的粒子群全局最优值Gbest:
Gbest=Lbestp (15)
⑥判断是否满足性能要求,若是,结束寻优,得到一组优化的模糊方程的局部方程参数;否则返回步骤③,继续迭代寻优,直至达到最大迭代次数itermax。
迭代终止时的Gbest即为对应于标准化后的第i个训练样本Xi的COD预报值和使COD排放达标的操作变量值。
作为优选的一种方案:所述的上位机还包括:判别模型更新模块,用于按设定的采样时间间隔,采集现场智能仪表信号,将得到的实测化学耗氧量与函数预报值比较,如果相对误差大于10%或实测COD数据不达标,则将DCS数据库中生产正常时的达标的新数据加入训练样本数据。
进一步,所述的上位机还包括:结果显示模块,用于将COD预报值和使COD排放达标的操作变量值传给DCS系统,并在DCS的控制站显示过程状态,同时通过DCS系统和现场总线将过程状态信息传递到现场操作站进行显示;同时,DCS系统将所得到的使COD排放达标的操作变量值作为新的操作变量设定值,自动执行COD排放达标操作。信号采集模块,用于依照设定的每次采样的时间间隔,从数据库中采集数据。
再进一步,所述的关键变量包括进入焚烧炉的废液流量、进入焚烧炉的空气流量和进入焚烧炉的燃料流量;所述的操作变量包括进入焚烧炉的空气流量和进入焚烧炉的燃料流量。
所述的农药废液焚烧炉有害物排放达标控制系统实现的有害物排放达标控制方法,所述的控制方法具体实现步骤如下:
1)、对农药废液焚烧炉有害物排放过程对象,根据工艺分析和操作分析,确定所用的关键变量,从DCS数据库中采集生产正常时所述变量的数据作为训练样本TX的输入矩阵,采集对应的COD和使COD排放达标的操作变量数据作为输出矩阵Y;
2)、将从DCS数据库输入的模型训练样本进行预处理,对训练样本中心化,即减去样本的平均值,然后对其进行标准化,使得其均值为0,方差为1。该处理采用以下算式过程来完成:
2.1)计算均值:
2.2)计算方差:
2.3)标准化:
其中,TXi为第i个训练样本,是从DCS数据库中采集的生产正常时的关键变量、COD和相应的使COD排放达标时的操作变量的数据,N为训练样本数,为训练样本的均值,X为标准化后的训练样本。σx表示训练样本的标准差,σ2 x表示训练样本的方差。
3)对从数据预处理模块传过来的训练样本,进行模糊化。设模糊方程系统中有c*个模糊群,模糊群k、j的中心分别为vk、vj,则标准化后的第i个训练样本Xi对于模糊群k的隶属度μik为:
式中,n为模糊分类过程中需要的分块矩阵指数,通常取作2,||·||为范数表达式。
使用以上隶属度值或者它的变形以获得新的输入矩阵,对于模糊群k,其输入矩阵变形为:
BP神经网络作为模糊方程系统的局部方程,设第k个BP神经网络模糊方程输出层的预测输出为输入为net,与此层相邻的隐含层中任一神经元l的输出为sl,则有:
式中,wlk是隐层神经元l与输出层神经元之间的连接权,f(·)为神经元的输出函数,通常取为Sigmoid函数,表示为:
式中,h为输出层神经元的阈值,θ0为陡度参数,用以调节Sigmoid函数的陡度;最后,由反模糊方法中的重心法得到最后的模糊方程系统的输出:
4)、采用粒子群算法对模糊方程中BP神经网络局部方程的wlk进行优化,具体步骤如下:①确定粒子群的优化参数为BP神经网络局部方程的wlk、粒子群个体数目popsize、最大循环寻优次数itermax、第p个粒子的初始位置rp、初始速度vp、局部最优值Lbestp以及整个粒子群的全局最优值Gbest。
②设定优化目标函数,将其转换为适应度,对每个局部模糊方程进行评价;通过相应的误差函数计算适应度函数,并认为误差大的粒子适应度小,粒子p的适应度函数表示为:
fp=1/(Ep+1) (10)
式中,Ep是模糊方程系统的误差函数,表示为:
式中,是模糊方程系统的预测输出,Oi为模糊方程系统的目标输出;
③按照如下公式,循环更新每个粒子的速度和位置,
vp(iter+1)=ω×vp(iter)+m1a1(Lbestp-rp(iter))+m2a2(Gbest-rp(iter))
(12)
rp(iter+1)=rp(iter)+vp(iter+1)(13)
式中,vp表示更新粒子p的速度,rp表示更新粒子p的位置,Lbestp表示更新粒子p的个体最优值,Gbest即为对应于标准化后的第i个训练样本Xi的COD预报值和使COD排放达标的操作变量值,iter表示循环次数,ω是粒子群算法中的惯性权重,m1、m2是对应的加速系数,a1、a2是[0,1]之间的随机数;
④对于粒子p,如果新的适应度大于原来的个体最优值,更新粒子的个体最优值:
Lbestp=fp (14)
⑤如果粒子p的个体最优值Lbestp大于原来的粒子群全局最优值Gbest,更新原来的粒子群全局最优值Gbest:
Gbest=Lbestp (15)
⑥判断是否满足性能要求,若是,结束寻优,得到一组优化的模糊方程的局部方程参数;否则返回步骤③,继续迭代寻优,直至达到最大迭代次数itermax。
迭代终止时得到Gbest即为对应于标准化后的第i个训练样本Xi的COD预报值和使COD排放达标的操作变量值。
作为优选的一种方案:所述的方法还包括:5)按设定的采样时间间隔,采集现场智能仪表信号,将得到的实测化学耗氧量与函数预报值比较,如果相对误差大于10%或实测COD数据不达标,则将DCS数据库中生产正常时的达标的新数据加入训练样本数据。
进一步,在所述的步骤4)中,将COD预报值和使COD排放达标的操作变量值传给DCS系统,并在DCS的控制站显示过程状态,同时通过DCS系统和现场总线将过程状态信息传递到现场操作站进行显示;同时,DCS系统将所得到的使COD排放达标的操作变量值作为新的操作变量设定值,自动执行COD排放达标操作。
再进一步,所述的关键变量包括进入焚烧炉的废液流量,进入焚烧炉的空气流量,进入焚烧炉的燃料流量;所述的操作变量包括进入焚烧炉的空气流量和进入焚烧炉的燃料流量。
本发明的技术构思为:本发明提供一种群智的农药废液焚烧炉有害物排放达标控制系统及方法,寻找到使得有害物排放达标的操作变量值。
本发明的有益效果主要表现在:1、建立了系统关键变量和化学耗氧量排放之间定量关系的在线软测量模型;2、迅速找到使得化学耗氧量排放达标的操作条件。
附图说明
图1是本发明所提出的系统的硬件结构图;
图2是本发明所提出的上位机的功能结构图。
具体实施方式
下面结合附图对本发明作进一步描述。本发明实施例用来解释说明本发明,而不是对本发明进行限制,在本发明的精神和权利要求的保护范围内,对本发明作出的任何修改和改变,都落入本发明的保护范围。
实施例1
参照图1、图2,群智的农药废液焚烧炉有害物排放达标控制系统,包括与焚烧炉1连接的现场智能仪表2、DCS系统以及上位机6,所述DCS系统包括数据接口3、控制站5和数据库4,所述现场智能仪表2与数据接口3连接,所述数据接口与控制站5、数据库4和上位机6连接,所述的上位机6包括:
数据预处理模块,用于将从DCS数据库输入的模型训练样本进行预处理,对训练样本中心化,即减去样本的平均值,然后对其进行标准化:
计算均值:
计算方差:
标准化:
其中,TXi为第i个训练样本,是从DCS数据库中采集的生产正常时的关键变量、COD和相应的使COD排放达标时的操作变量的数据,N为训练样本数,为训练样本的均值,X为标准化后的训练样本。σx表示训练样本的标准差,σ2 x表示训练样本的方差。
模糊方程模块,对从数据预处理模块传过来的标准化后的训练样本X,进行模糊化。设模糊方程系统中有c*个模糊群,模糊群k、j的中心分别为vk、vj,则标准化后的第i个训练样本Xi对于模糊群k的隶属度μik为:
式中,n为模糊分类过程中需要的分块矩阵指数,通常取作2,||·||为范数表达式。
使用以上隶属度值或者它的变形以获得新的输入矩阵,对于模糊群k,其输入矩阵变形为:
式中,wlk是隐层神经元l与输出层神经元之间的连接权,f(·)为神经元的输出函数,通常取为Sigmoid函数,表示为:
式中,h为输出层神经元的阈值,θ0为陡度参数,用以调节Sigmoid函数的陡度;最后,由反模糊方法中的重心法得到最后的模糊方程系统的输出:
粒子群算法优化模块,用于采用粒子群算法对模糊方程中BP神经网络局部方程的wlk进行优化,具体步骤如下:
①确定粒子群的优化参数为BP神经网络局部方程的wlk、粒子群个体数目popsize、最大循环寻优次数itermax、第p个粒子的初始位置rp、初始速度vp、局部最优值Lbestp以及整个粒子群的全局最优值Gbest。
②设定优化目标函数,将其转换为适应度,对每个局部模糊方程进行评价;通过相应的误差函数计算适应度函数,并认为误差大的粒子适应度小,粒子p的适应度函数表示为:
fp=1/(Ep+1) (10)
式中,Ep是模糊方程系统的误差函数,表示为:
③按照如下公式,循环更新每个粒子的速度和位置,
vp(iter+1)=ω×vp(iter)+m1a1(Lbestp-rp(iter))+m2a2(Gbest-rp(iter))
(12)
rp(iter+1)=rp(iter)+vp(iter+1) (13)
式中,vp表示更新粒子p的速度,rp表示更新粒子p的位置,Lbestp表示更新粒子p的个体最优值,Gbest即为对应于标准化后的第i个训练样本Xi的COD预报值和使COD排放达标的操作变量值,iter表示循环次数,ω是粒子群算法中的惯性权重,m1、m2是对应的加速系数,a1、a2是[0,1]之间的随机数;
④对于粒子p,如果新的适应度大于原来的个体最优值,更新粒子的个体最优值:
Lbestp=fp (14)
⑤如果粒子p的个体最优值Lbestp大于原来的粒子群全局最优值Gbest,更新原来的粒子群全局最优值Gbest:
Gbest=Lbestp (15)
⑥判断是否满足性能要求,若是,结束寻优,得到一组优化的模糊方程的局部方程参数;否则返回步骤③,继续迭代寻优,直至达到最大迭代次数itermax。
迭代终止时的Gbest即为对应于标准化后的第i个训练样本Xi的COD预报值和使COD排放达标的操作变量值。
所述上位机6还包括:信号采集模块11,用于依照设定的每次采样的时间间隔,从数据库中采集数据;
所述的上位机6还包括:判别模型更新模块12,按设定的采样时间间隔,采集现场智能仪表信号,将得到的实测COD与函数计算值比较,如果相对误差大于10%,则将新数据加入训练样本数据。
所述的关键变量包括进入焚烧炉的废液流量、进入焚烧炉的空气流量和进入焚烧炉的燃料流量;所述的操作变量包括进入焚烧炉的空气流量和进入焚烧炉的燃料流量。
所述系统还包括DCS系统,所述的DCS系统由数据接口3、控制站5、数据库4构成;智能仪表2、DCS系统、上位机6通过现场总线依次相连;上位机6还包括结果显示模块10,用于将COD预报值和使COD排放达标的操作变量值传给DCS系统,并在DCS的控制站显示过程状态,同时通过DCS系统和现场总线将过程状态信息传递到现场操作站进行显示;同时,DCS系统将所得到的使COD排放达标的操作变量值作为新的操作变量设定值,自动执行COD排放达标操作。
当废液焚烧过程已配有DCS系统时,样本实时动态数据的检测、存储利用DCS系统的实时和历史数据库,得到COD预报值和使COD排放达标的操作变量值的功能主要在上位机上完成。
当废液焚烧过程没有配备DCS系统时,采用数据存储器来替代DCS系统的实时和历史数据库的数据存储功能,并将得到COD预报值和使COD排放达标的操作变量值的功能系统制造成包括I/O元件、数据存储器、程序存储器、运算器、显示模块几大构件的不依赖于DCS系统的一个独立的完整的片上系统,在不管焚烧过程是否配备DCS的情况下,都能够独立使用,更有益于推广使用。
实施例2
参照图1、图2,群智的农药废液焚烧炉有害物排放达标控制方法,所述控制方法具体实现步骤如下:
1)、对农药废液焚烧炉有害物排放过程对象,根据工艺分析和操作分析,确定所用的关键变量,从DCS数据库中采集生产正常时所述变量的数据作为训练样本TX的输入矩阵,采集对应的COD和使COD排放达标的操作变量数据作为输出矩阵Y;
2)、将从DCS数据库输入的模型训练样本进行预处理,对训练样本中心化,即减去样本的平均值,然后对其进行标准化,使得其均值为0,方差为1。该处理采用以下算式过程来完成:
2.1)计算均值:
2.2)计算方差:
2.3)标准化:
其中,TXi为第i个训练样本,是从DCS数据库中采集的生产正常时的关键变量、COD和相应的使COD排放达标时的操作变量的数据,N为训练样本数,为训练样本的均值,X为标准化后的训练样本。σx表示训练样本的标准差,σ2 x表示训练样本的方差。
3)对从数据预处理模块传过来的经过标准化后的训练样本,进行模糊化。设模糊方程系统中有c*个模糊群,模糊群k、j的中心分别为vk、vj,则标准化后的第i个训练样本Xi对于模糊群k的隶属度μik为:
式中,n为模糊分类过程中需要的分块矩阵指数,通常取作2,||·||为范数表达式。
使用以上隶属度值或者它的变形以获得新的输入矩阵,对于模糊群k,其输入矩阵变形为:
式中,wlk是隐层神经元l与输出层神经元之间的连接权,f(·)为神经元的输出函数,通常取为Sigmoid函数,表示为:
式中,h为输出层神经元的阈值,θ0为陡度参数,用以调节Sigmoid函数的陡度;最后,由反模糊方法中的重心法得到最后的模糊方程系统的输出:
4)采用粒子群算法对模糊方程中BP神经网络局部方程的wlk进行优化,具体步骤如下:①确定粒子群的优化参数为BP神经网络局部方程的wlk、粒子群个体数目popsize、最大循环寻优次数itermax、第p个粒子的初始位置rp、初始速度vp、局部最优值Lbestp以及整个粒子群的全局最优值Gbest。
②设定优化目标函数,将其转换为适应度,对每个局部模糊方程进行评价;通过相应的误差函数计算适应度函数,并认为误差大的粒子适应度小,粒子p的适应度函数表示为:
fp=1/(Ep+1) (10)
式中,Ep是模糊方程系统的误差函数,表示为:
③按照如下公式,循环更新每个粒子的速度和位置,
vp(iter+1)=ω×vp(iter)+m1a1(Lbestp-rp(iter))+m2a2(Gbest-rp(iter))
(12)
rp(iter+1)=rp(iter)+vp(iter+1) (13)
式中,vp表示更新粒子p的速度,rp表示更新粒子p的位置,Lbestp表示更新粒子p的个体最优值,Gbest即为对应于标准化后的第i个训练样本Xi的COD预报值和使COD排放达标的操作变量值,iter表示循环次数,ω是粒子群算法中的惯性权重,m1、m2是对应的加速系数,a1、a2是[0,1]之间的随机数;
④对于粒子p,如果新的适应度大于原来的个体最优值,更新粒子的个体最优值:
Lbestp=fp (14)
⑤如果粒子p的个体最优值Lbestp大于原来的粒子群全局最优值Gbest,更新原来的粒子群全局最优值Gbest:
Gbest=Lbestp (15)
⑥判断是否满足性能要求,若是,结束寻优,得到一组优化的模糊方程的局部方程参数;否则返回步骤③,继续迭代寻优,直至达到最大迭代次数itermax。
迭代终止时的Gbest即为对应于标准化后的第i个训练样本Xi的COD预报值和使COD排放达标的操作变量值。
所述的方法还包括:5)、按设定的采样时间间隔,采集现场智能仪表信号,将得到的实测COD与函数预报值比较,如果相对误差大于10%或实测COD数据不达标,则将DCS数据库中生产正常时的达标的新数据加入训练样本数据。
6)、在所述的步骤4)中,将COD预报值和使COD排放达标的操作变量值传给DCS系统,并在DCS的控制站显示过程状态,同时通过DCS系统和现场总线将过程状态信息传递到现场操作站进行显示;同时,DCS系统将所得到的使COD排放达标的操作变量值作为新的操作变量设定值,自动执行COD排放达标操作。
所述的关键变量包括进入焚烧炉的废液流量,进入焚烧炉的空气流量,进入焚烧炉的燃料流量;所述的操作变量包括进入焚烧炉的空气流量和进入焚烧炉的燃料流量。
Claims (2)
1.一种群智的农药废液焚烧炉有害物排放达标控制系统,包括焚烧炉、智能仪表、DCS系统、数据接口以及上位机,所述的DCS系统包括控制站和数据库;所述现场智能仪表与DCS系统连接,所述DCS系统与上位机连接,其特征在于:所述的上位机包括:
数据预处理模块,用于将从DCS数据库输入的模型训练样本进行预处理,对训练样本中心化,即减去样本的平均值,然后对其进行标准化:
计算均值:
计算方差:
标准化:
其中,TXi为第i个训练样本,是从DCS数据库中采集的生产正常时的关键变量、化学耗氧量(COD)和相应的使COD排放达标时的操作变量的数据,N为训练样本数,为训练样本的均值,X为标准化后的训练样本。σx表示训练样本的标准差,σ2 x表示训练样本的方差。模糊方程模块,对从数据预处理模块传过来的标准化后的训练样本X,进行模糊化。设模糊方程系统中有c*个模糊群,模糊群k、j的中心分别为vk、vj,则标准化后的第i个训练样本Xi对于模糊群k的隶属度μik为:
式中,n为模糊分类过程中需要的分块矩阵指数,通常取作2,||·||为范数表达式。
使用以上隶属度值或者它的变形以获得新的输入矩阵,对于模糊群k,其输入矩阵变形为:
式中,wlk是隐层神经元l与输出层神经元之间的连接权,f(·)为神经元的输出函数,通常取为Sigmoid函数,表示为:
式中,h为输出层神经元的阈值,θ0为陡度参数,用以调节Sigmoid函数的陡度;最后,由反模糊方法中的重心法得到最后的模糊方程系统的输出:
粒子群算法优化模块,用于采用粒子群算法对模糊方程中BP神经网络局部方程的wlk进行优化,具体步骤如下:
①确定粒子群的优化参数为BP神经网络局部方程的wlk、粒子群个体数目popsize、最大循环寻优次数itermax、第p个粒子的初始位置rp、初始速度vp、局部最优值Lbestp以及整个粒子群的全局最优值Gbest。
②设定优化目标函数,将其转换为适应度,对每个局部模糊方程进行评价;通过相应的误差函数计算适应度函数,并认为误差大的粒子适应度小,粒子p的适应度函数表示为:
fp=1/(Ep+1)(10)
式中,Ep是模糊方程系统的误差函数,表示为:
式中,是模糊方程系统的预测输出,Oi为模糊方程系统的目标输出;
③按照如下公式,循环更新每个粒子的速度和位置,
vp(iter+1)=ω×vp(iter)+m1a1(Lbestp-rp(iter))+m2a2(Gbest-rp(iter))
(12)
rp(iter+1)=rp(iter)+vp(iter+1) (13)
式中,vp表示更新粒子p的速度,rp表示更新粒子p的位置,Lbestp表示更新粒子p的个体最优值,Gbest即为对应于标准化后的第i个训练样本Xi的COD预报值和使COD排放达标的操作变量值,iter表示循环次数,ω是粒子群算法中的惯性权重,m1、m2是对应的加速系数,a1、a2是[0,1]之间的随机数;
④对于粒子p,如果新的适应度大于原来的个体最优值,更新粒子的个体最优值:
Lbestp=fp (14)
⑤如果粒子p的个体最优值Lbestp大于原来的粒子群全局最优值Gbest,更新原来的粒子群全局最优值Gbest:
Gbest=Lbestp (15)
⑥判断是否满足性能要求,若是,结束寻优,得到一组优化的模糊方程的局部方程参数;否则返回步骤③,继续迭代寻优,直至达到最大迭代次数itermax。
迭代终止时的Gbest即为对应于标准化后的第i个训练样本Xi的COD预报值和使COD排放达标的操作变量值。
所述的上位机还包括:判别模型更新模块,用于按设定的采样时间间隔,采集现场智能仪表信号,将得到的实测化学耗氧量与函数预报值比较,如果相对误差大于10%或实测COD数据不达标,则将DCS数据库中生产正常时的达标的新数据加入训练样本数据,更新模型。
结果显示模块,用于将COD预报值和使COD排放达标的操作变量值传给DCS系统,在DCS的控制站显示,并通过DCS系统和现场总线传递到现场操作站进行显示;同时,DCS系统将所得到的使COD排放达标的操作变量值作为新的操作变量设定值,自动执行COD排放达标操作。信号采集模块,用于依照设定的每次采样的时间间隔,从数据库中采集数据。
所述的关键变量包括进入焚烧炉的废液流量、进入焚烧炉的空气流量和进入焚烧炉的燃料流量;所述的操作变量包括进入焚烧炉的空气流量和进入焚烧炉的燃料流量。
2.一种用如权利要求1所述的群智的农药废液焚烧炉有害物排放达标控制系统实现的控制方法,其特征在于:所述控制方法具体实现步骤如下:
1)、对农药废液焚烧炉有害物排放过程对象,根据工艺分析和操作分析,确定所用的关键变量,从DCS数据库中采集生产正常时所述变量的数据作为训练样本TX的输入矩阵,采集对应的COD和使COD排放达标的操作变量数据作为输出矩阵Y;
2)、将从DCS数据库输入的模型训练样本进行预处理,对训练样本中心化,即减去样本的平均值,然后对其进行标准化,使得其均值为0,方差为1。该处理采用以下算式过程来完成:
2.1)计算均值:
2.2)计算方差:
2.3)标准化:
其中,TXi为第i个训练样本,是从DCS数据库中采集的生产正常时的关键变量、化学耗氧量(COD)和相应的使COD排放达标时的操作变量的数据,N为训练样本数,为训练样本的均值,X为标准化后的训练样本。σx表示训练样本的标准差,σ2 x表示训练样本的方差。
3)对从数据预处理模块传过来的训练样本,进行模糊化。设模糊方程系统中有c*个模糊群,模糊群k、j的中心分别为vk、vj,则标准化后的第i个训练样本Xi对于模糊群k的隶属度μik为:
式中,n为模糊分类过程中需要的分块矩阵指数,通常取作2,||·||为范数表达式。
使用以上隶属度值或者它的变形以获得新的输入矩阵,对于模糊群k,其输入矩阵变形为:
其中func(μik)为隶属度值μik的变形函数,一般取exp(μik)等,Φik(Xi,μik)表示第i个输入变量Xi及其模糊群k的隶属度μik所对应的新的输入矩阵。
BP神经网络作为模糊方程系统的局部方程,设第k个BP神经网络模糊方程输出层的预测输出为输入为net,与此层相邻的隐含层中任一神经元l的输出为yl,则有:
式中,wlk是隐层神经元l与输出层神经元之间的连接权,f(·)为神经元的输出函数,通常取为Sigmoid函数,表示为:
式中,h为输出层神经元的阈值,θ0为陡度参数,用以调节Sigmoid函数的陡度;最后,由反模糊方法中的重心法得到最后的模糊方程系统的输出:
4)采用粒子群算法对模糊方程中BP神经网络局部方程的wlk进行优化,具体步骤如下:①确定粒子群的优化参数为BP神经网络局部方程的wlk、粒子群个体数目popsize、最大循环寻优次数itermax、第p个粒子的初始位置rp、初始速度vp、局部最优值Lbestp以及整个粒子群的全局最优值Gbest。
②设定优化目标函数,将其转换为适应度,对每个局部模糊方程进行评价;通过相应的误差函数计算适应度函数,并认为误差大的粒子适应度小,粒子p的适应度函数表示为:
fp=1/(Ep+1) (10)
式中,Ep是模糊方程系统的误差函数,表示为:
式中,是模糊方程系统的预测输出,Oi为模糊方程系统的目标输出;
③按照如下公式,循环更新每个粒子的速度和位置,
vp(iter+1)=ω×vp(iter)+m1a1(Lbestp-rp(iter))+m2a2(Gbest-rp(iter)) (12)
rp(iter+1)=rp(iter)+vp(iter+1) (13)
式中,vp表示更新粒子p的速度,rp表示更新粒子p的位置,Lbestp表示更新粒子p的个体最优值,Gbest即为对应于标准化后的第i个训练样本Xi的COD预报值和使COD排放达标的操作变量值,iter表示循环次数,ω是粒子群算法中的惯性权重,m1、m2是对应的加速系数,a1、a2是[0,1]之间的随机数;
④对于粒子p,如果新的适应度大于原来的个体最优值,更新粒子的个体最优值:
Lbestp=fp (14)
⑤如果粒子p的个体最优值Lbestp大于原来的粒子群全局最优值Gbest,更新原来的粒子群全局最优值Gbest:
Gbest=Lbestp (15)
⑥判断是否满足性能要求,若是,结束寻优,得到一组优化的模糊方程的局部方程参数;否则返回步骤③,继续迭代寻优,直至达到最大迭代次数itermax。
迭代终止时的Gbest即为对应于标准化后的第i个训练样本Xi的COD预报值和使COD排放达标的操作变量值。
所述方法还包括:5)、判别模型更新模块,用于按设定的采样时间间隔,采集现场智能仪表信号,将得到的实测化学耗氧量与函数预报值比较,如果相对误差大于10%或实测COD数据不达标,则将DCS数据库中生产正常时的达标的新数据加入训练样本数据,更新模型。
6)、在所述的步骤4)中得到的COD预报值和使COD排放达标的操作变量值,将结果传给DCS系统,在DCS的控制站显示,并通过DCS系统和现场总线传递到现场操作站进行显示;同时,DCS系统将所得到的使COD排放达标的操作变量值作为新的操作变量设定值,自动执行COD排放达标操作。
所述的关键变量包括进入焚烧炉的废液流量、进入焚烧炉的空气流量和进入焚烧炉的燃料流量;所述的操作变量包括进入焚烧炉的空气流量和进入焚烧炉的燃料流量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310433864.8A CN103472729B (zh) | 2013-09-22 | 2013-09-22 | 群智的农药废液焚烧炉有害物排放达标控制系统及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310433864.8A CN103472729B (zh) | 2013-09-22 | 2013-09-22 | 群智的农药废液焚烧炉有害物排放达标控制系统及方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103472729A true CN103472729A (zh) | 2013-12-25 |
CN103472729B CN103472729B (zh) | 2016-01-06 |
Family
ID=49797621
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310433864.8A Expired - Fee Related CN103472729B (zh) | 2013-09-22 | 2013-09-22 | 群智的农药废液焚烧炉有害物排放达标控制系统及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103472729B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108319146A (zh) * | 2018-03-09 | 2018-07-24 | 西安西热控制技术有限公司 | 一种径向基神经网络基于离散粒子群训练的方法 |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106931453B (zh) * | 2017-02-27 | 2018-02-16 | 浙江大学 | 循环流化床生活垃圾焚烧锅炉NOx排放的预测系统及方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040143149A1 (en) * | 2002-08-02 | 2004-07-22 | Decourcy Michael Stanley | Method for reducing nitrogen oxide emissions in industrial process |
CN101457264A (zh) * | 2008-12-29 | 2009-06-17 | 杭州电子科技大学 | 高炉炉温优化控制方法 |
CN101763084A (zh) * | 2009-12-29 | 2010-06-30 | 浙江大学 | 农药生产废液焚烧炉化学耗氧量排放最小化系统及方法 |
CN101763085A (zh) * | 2009-12-29 | 2010-06-30 | 浙江大学 | 农药生产废液焚烧炉炉温最佳化系统及方法 |
-
2013
- 2013-09-22 CN CN201310433864.8A patent/CN103472729B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040143149A1 (en) * | 2002-08-02 | 2004-07-22 | Decourcy Michael Stanley | Method for reducing nitrogen oxide emissions in industrial process |
CN101457264A (zh) * | 2008-12-29 | 2009-06-17 | 杭州电子科技大学 | 高炉炉温优化控制方法 |
CN101763084A (zh) * | 2009-12-29 | 2010-06-30 | 浙江大学 | 农药生产废液焚烧炉化学耗氧量排放最小化系统及方法 |
CN101763085A (zh) * | 2009-12-29 | 2010-06-30 | 浙江大学 | 农药生产废液焚烧炉炉温最佳化系统及方法 |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108319146A (zh) * | 2018-03-09 | 2018-07-24 | 西安西热控制技术有限公司 | 一种径向基神经网络基于离散粒子群训练的方法 |
CN108319146B (zh) * | 2018-03-09 | 2020-08-11 | 西安西热控制技术有限公司 | 一种径向基神经网络基于离散粒子群训练的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103472729B (zh) | 2016-01-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101763084B (zh) | 农药生产废液焚烧炉化学耗氧量排放最小化系统及方法 | |
CN101782771B (zh) | 基于免疫神经网络的纺丝工艺智能优化设计方法 | |
CN101763085B (zh) | 农药生产废液焚烧炉炉温最佳化系统及方法 | |
CN103195727B (zh) | 一种风机在线状态监测评估装置 | |
CN103472865A (zh) | 智能最小二乘的农药废液焚烧炉炉温最佳化系统及方法 | |
CN106682316A (zh) | 一种基于尖峰径向基神经网络的出水总磷实时监测系统 | |
CN106199174A (zh) | 基于迁移学习的挤压机能耗异常预测方法 | |
CN110824923A (zh) | 一种基于深度学习与云计算的污水处理控制方法及其系统 | |
CN103674778B (zh) | Rbf粒子群优化的工业熔融指数软测量仪表和方法 | |
CN113591374B (zh) | 基于卷积简单循环单元网络的水泥成品比表面积预测方法 | |
CN111260117A (zh) | 一种基于气象因子的ca-narx水质预测方法 | |
CN103675006A (zh) | 最小二乘的工业熔融指数软测量仪表和方法 | |
CN103472729B (zh) | 群智的农药废液焚烧炉有害物排放达标控制系统及方法 | |
CN103488145B (zh) | 群智模糊网络的焚烧炉有害物排放达标控制系统及方法 | |
CN103488089B (zh) | 自适应的农药废液焚烧炉有害物排放达标控制系统及方法 | |
CN103472727A (zh) | 群智加权的农药焚烧炉有害物排放达标控制系统及方法 | |
CN115700324A (zh) | 一种瓦温预警的方法、装置和电子设备 | |
CN103472867A (zh) | 支持向量机的农药生产废液焚烧炉炉温最佳化系统及方法 | |
CN103488090B (zh) | 群智机器学习的焚烧炉有害物排放达标控制系统及方法 | |
CN103472728A (zh) | 机器学习的农药焚烧炉有害物排放达标控制系统及方法 | |
CN103472721A (zh) | 自适应机器学习的农药废液焚烧炉炉温最佳化系统及方法 | |
CN103488088B (zh) | 误差反向传播的焚烧炉有害物排放达标控制系统及方法 | |
CN103675012B (zh) | Bp粒子群优化的工业熔融指数软测量仪表和方法 | |
CN103488087B (zh) | 最优的农药废液焚烧炉有害物排放达标控制系统及方法 | |
CN103488084B (zh) | 模糊网络的农药焚烧炉有害物排放达标控制系统及方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160106 Termination date: 20180922 |
|
CF01 | Termination of patent right due to non-payment of annual fee |