申请日2011.05.27
公开(公告)日2012.01.11
IPC分类号G01N33/18; G06N3/08
摘要
一种污水处理出水生化需氧量软测量方法属于污水处理领域。污水处理过程机理复杂，为使污水处理系统处于良好的运行工况，获得更好的出水水质，需要对污水处理系统中的过程参数及水质参数进行检测。本发明针对生化需氧量不易在线测得的问题，应用一种动态结构自组织神经网络建立软测量模型。根据隐含层神经元的输出对整个网络输出的贡献大小，动态地增长-删减神经元，并能自动地调整节点的增长-修剪训练步数，实现了神经网络网络结构的动态变化，对出水生化需氧量进行预测。根据预测结果，能够对污水处理过程中的相关控制环节及生化反应中的物料进行调节，提高污水处理的出水水质质量，为污水处理过程优化控制和节能运行提供技术保障。
权利要求书 [支持框选翻译]
1.一种污水处理生化需氧量软测量方法，其特征在于，包括以下步骤：
(1).对污水处理过程机理进行分析，初步确定与生化需氧量相关的辅助变量;
(2).在污水处理工业现场采集k组步骤(1)中确定的辅助变量及生化需氧量的样本数据;
(3).运用主元分析方法对步骤(2)的数据再次进行分析，精选出跟出水生化需氧量相关性高的辅助变量;
(4).用神经网络结构自组织设计方法对污水处理过程现场的进出水数据建立模型，预测输出下一时刻的出水生化需氧量。
2.根据权利要求1所述的一种污水处理生化需氧量软测量方法，其特征在于：用神经网络结构自组织设计方法对污水处理过程现场的进出水数据建立模型，预测输出下一时刻的出水生化需氧量结果，具体包括以下步骤：
(1).建立一个神经网络模型，神经网络包括输入层、隐含层及输出层，神经网络模型的输入为权利要求1中精选出的辅助变量，输出为生化需氧量;
(2).对神经网络进行μ步训练，每训练一步计算一次神经网络的误差，误差函数定义如下：
其中：
y(t)表示将采集到的精选出的辅助变量的第t组样本数据输入到神经网络所产生的实际输出值，
yd(t)表示神经网络的期望输出，即与精选出的辅助变量的第t组样本数据对应的已经采集到的生化需氧量的样本数据，
k表示采集到的样本数;
将μ步训练结束后的神经网络的误差计为E0p;
(3).判断神经网络训练至最后一步时的误差是否满足设定目标值，若满足，跳转至(9)，若不满足，继续向下执行;
(4).计算隐含层每个神经元的全局显著性指数，全局显著性指数计算公式如下：
其中，
GSIj为隐含层第j个神经元的全局显著性指数，
K为隐含层神经元个数，
SIj为隐含层第j个神经元的显著性指数，计算公式如下;
其中，
M为输入层神经元个数，
wijF为神经网络训练终止时输入层第i个神经元与隐含层第j个神经元之间的连接权值，
wijI为神经网络训练初始时输入层第i个神经元与隐含层第j个神经元之间的连接权值，
σj为隐含层第j个神经元输出的方差;
(5).根据隐含层神经元变化条件，对每个隐含层神经元进行增长或删减;
如果隐含层第j个神经元的全局显著性指数满足
GSIj≥Gth1 (4)
则由原来的一个神经元增长为三个神经元，其中，Gth1为设定的增长阈值;相应的隐含层新增神经元的权值调整为：
∑αm=1 (6)
其中，
wijm表示新增隐含层神经元与输入层第i个神经元之间的连接权值，
wjm表示新增隐含层神经元与输出层神经元之间的连接权值，
wij表示输入层第i个神经元和隐含层第j个神经元之间的权值，
wj表示隐含层第j个神经元和输出层神经元之间的权值。
αm的选择服从均值为0，方差为1的高斯分布;
如果隐含层第j个神经元的全局显著性指数满足
Gth2≤GSIj 则由原来的一个神经元增长为两个神经元，其中，Gth2为设定的增长阈值：相应的隐含层新增节点的权值调整为： 如果隐含层第j个神经元全局显著性较低，满足 GSIj≤Pth (9) 则将其删除;其中，Pth为设定的删减阈值， (6).对神经网络再次进行ε步训练，每训练一步计算一次神经网络的误差，训练结束后计算平均误差Emean，Ei为神经网络训练至第i步时的误差;将第一步训练结束时神经网络的误差计为E1P; (7).确定自适应训练步长; 当Emean<ΔEmax·E0p时，ΔEmax=|E0p-E1P|，训练步长调整为： 当Emean>ΔEmax·E0p时，训练步长调整为 ε=γ·ε (12) 其中，γ>1; (8)根据(7)得到的训练步长对神经网络再次进行训练并跳转至(3); (9)神经网络训练结束，利用训练好的神经网络对出水生化需氧量进行预测。 说明书 [支持框选翻译] 一种污水处理生化需氧量软测量方法 技术领域 本发明涉及一种污水处理出水生化需氧量(BOD)的软测量方法，属于污水处理领域。 背景技术 污水生物处理过程机理复杂，为使污水处理系统处于良好的运行工况，获得更好的出水水质，必须及时取得污水处理系统当中重要的过程参数及水质参数。目前通常采用两种方法对各种参数进行在线检测：一种是采用新型传感器，以硬件的方式实现难测参数的检测。但由于污水中有机物千差万别，若根据每种有机物的特性研发出相对应的传感器，将是一个耗资大、历时长的工程，而且这种传感器具有仪器造价高、寿命短，测量范围窄、稳定性差等缺点;另一种是软测量的方法。现有的软测量方法主要有基于工艺机理分析的软测量方法和基于神经网络的软测量方法。基于工艺机理分析的软测量方法需要对工艺过程机理深刻认识，否则不能准确预测出不易测量的参数值。基于神经网络的软测量方法适用于非线性和不确定性系统，但传统的神经网络软测量方法不能根据所预测参数的复杂程度自动地调整神经网络结构，具有预测精度较低，时间较长的缺点。 发明内容 为了解决污水处理过程中难以实时检测的水质参数在较短的时间内准确测量的问题，本发明提供了一种利用结构自组织神经网络建立软测量模型，对难测数据进行准确预测的技术方法。具体方案如下： (1).对污水处理过程机理进行分析，根据污水处理流程需经过的沉沙池、初沉池、曝气池、二沉池及接触池，确定出水生化需氧量与污水流经各处理池时所检测到的其他水质参数的变化情况，初步确定与出水生化需氧量相关的辅助变量; (2).在污水处理工业现场采集k组步骤(1)中确定的辅助变量及生化需氧量的样本数据; (3).运用主元分析方法对步骤(2)的数据再次进行分析，精选出跟出水生化需氧量相关性高的辅助变量; 将样本数据用零均值标准化方法进行归一化处理，然后计算辅助变量的累积方差贡献率，若累积方差贡献率大于85%，则认为该辅助变量与所需预测的出水生化需氧量相关性高，可将其精选出，进行下一步的工作;其中具体方法步骤参见《软测量技术原理与应用》一书。 (4).用神经网络结构自组织设计方法对污水处理过程现场的进出水数据建立模型，预测输出下一时刻的出水生化需氧量。具体步骤为： ①建立一个神经网络模型，神经网络包括输入层、隐含层及输出层，神经网络模型的输入为步骤(3)中精选出的M个辅助变量，输出为生化需氧量; 神经网络的输入层有M个神经元，隐含层有N个神经元，输出层有1个神经元。输入层的输入为V=[v1，v2，...vM]，对于输入层第i个神经元，其输出为： i=1，2，Λ，M (3) 其中，fI为输入层的激励函数。 隐含层第j个神经元的输入为 其中 为输入层第i个神经元与隐含层第j个神经元之间的连接权值;M为输入层神经元个数。 隐含层第j个神经元的输出为 其中，fH为隐含层的激励函数， 输出层神经元的输入为： 其中， 为隐含层第j个神经元与输出层神经元之间的连接权值，N为隐含层的神经元数; 输出层神经元的输出，即实际网络的输出为 yO=fO(xO) (7) 其中，fO为输出层神经元的激励函数函数 ②对神经网络进行μ步训练，每训练一步计算一次神经网络的误差，误差函数定义如下： 其中： y(t)表示将采集到的精选出的辅助变量的第t组样本数据输入到神经网络所产生的实际输出值， yd(t)表示神经网络的期望输出，即与精选出的辅助变量的第t组样本数据对应的已经采集到的生化需氧量的样本数据， k表示采集到的样本数; 将μ步训练结束后的神经网络的误差计为E0p; ③判断神经网络训练至最后一步时的误差是否满足设定目标值，若满足，跳转至(9)，若不满足，继续向下执行; ④计算隐含层每个神经元的全局显著性指数，全局显著性指数计算公式如下： 其中， GSIj为隐含层第j个神经元的全局显著性指数， K为隐含层神经元个数， SIj为隐含层第j个神经元的显著性指数，计算公式如下; 其中， M为输入层神经元个数， wijF为神经网络训练终止时输入层第i个神经元与隐含层第j个神经元之间的连接权值， wijI为神经网络训练初始时输入层第i个神经元与隐含层第j个神经元之间的连接权值， σj为隐含层第j个神经元输出的方差; ⑤根据隐含层神经元变化条件，对每个隐含层神经元进行增长或删减; 如果隐含层第j个神经元的全局显著性指数满足 GSIj≥Gth1 (11) 则由原来的一个神经元增长为三个神经元，其中，Gth1为设定的增长阈值;相应的隐含层新增神经元的权值调整为： ∑αm=1 (13) 其中， wijm表示新增隐含层神经元与输入层第i个神经元之间的连接权值， wjm表示新增隐含层神经元与输出层神经元之间的连接权值， wij表示输入层第i个神经元和隐含层第j个神经元之间的权值。 wj表示隐含层第j个神经元和输出层神经元之间的权值。 αm的选择服从均值为0，方差为1的高斯分布; 如果隐含层第j个神经元的全局显著性指数满足 Gth2≤GSIj 则由原来的一个神经元增长为两个神经元，其中，Gth2为设定的增长阈值;相应的隐含层新增节点的权值调整为： 如果隐含层第j个神经元全局显著性较低，满足 GSIj≤Pth (16) 则将其删除;其中，Pth为设定的删减阈值， ⑥对神经网络再次进行ε步训练，每训练一步计算一次神经网络的误差，训练结束后计算平均误差Emean， Ei为神经网络训练至第i步时的误差;将第一步训练结束时神经网络的误差计为E1P; ⑦确定自适应训练步长; 当Emean<ΔEmax·E0p时，ΔEmax=|E0p-E1P|，训练步长调整为： 当Emean>ΔEmax·E0p时，训练步长调整为 ε=γ·ε (18) 其中，γ>1; ⑧根据⑦得到的训练步长对神经网络再次进行训练并跳转至(3); ⑨神经网络训练结束。 本发明的优势主要体现在： 本发明针对污水处理生化需氧量难以实时测量的问题，采用了基于结构自组织神经网络建立出水生化需氧量的软测量模型方法，通过神经网络结构自组织设计能够在较短时间内预测出精度较高的生化需氧量值，避免了参数检测的大滞后问题，并且具有操作简单，易于维护等特点。
