[发明专利]基于储备池网络的污水处理出水总磷TP软测量方法

权利要求书

1.基于储备池网络的污水处理出水总磷TP软测量方法，其特征在于，包括以下步骤：

第一步：数据预处理及辅助变量的确定

采集污水处理厂实际水质参数数据，对数据进行预处理，通过主成分分析选取与出水总磷TP相关性强的辅助变量，最终选取进水总磷TP、出水温度T、好氧前段溶解氧DO、入水油类，出水油类，以及出水氨氮NH₄-N作为出水总磷TP测量的辅助变量，即软测量模型的输入变量；

第二步：设计基于贡献率的储备池网络结构优化算法

通过计算储备池神经元对网络输出的信息贡献程度，评价储备池神经元的贡献度；根据储备池神经元的贡献度优化网络结构，提高网络的泛化能力，所以，设计有效的储备池网络结构优化算法可以保证基于储备池网络建立的出水总磷TP软测量模型预测的准确性，提高对污水处理过程中复杂动态问题的适应能力；

1)构建一个储备池规模为N的储备池网络，网络通过储备池把输入数据转换到高维状态空间，且每一维对最终输出的贡献都不一样，为了分析储备池神经元的贡献度，首先定义两个神经元X和Y间信息连接强度函数：

IC(X；Y)＝H(X)+H(Y)-H(X,Y) (1)

其中，H(X)为神经元X的熵，H(Y)为神经元Y的熵，分别描述了神经元X，Y所包含的信息量；H(X,Y)为神经元X和Y的联合熵，描述了两个神经元间共有的信息量；

依据Shannon熵定义，H(X)、H(Y)和H(X,Y)可由式(2)-(4)计算：

H(X)＝-∑_x∈Xρ(x)logρ(x) (2)

H(Y)＝-∑_y∈Yρ(y)logρ(y) (3)

则公式(1)可以写成：

针对于第i个储备池神经元的贡献度进行评定，贡献度较小的输出连接权值将被修剪；第i个储备池神经元的贡献度C_i表示如下：

其中,ρ(x)、ρ(y)分别表示神经元X和神经元Y的边缘概率密度；ρ(x,y)表示神经元X和Y的联合概率密度；ρ_X(x)和ρ_Y(y)与ρ(x)和ρ(y)是相同定义，分别表示神经元X和神经元Y的边缘概率密度；m_i表示第i个储备池神经元对网络输出的信息贡献程度，N为储备池神经元总个数，具体计算公式如下：

其中，I(x_i,y)表示第i个储备池神经元x_i与神经元Y间信息连接强度；

2)将第i个储备池神经元贡献度C_i与设定的修剪阈值ρ进行对比，判断是否修剪第i个储备池神经元的输出连接权值，具体的修剪判别机制如下：

其中，w_i为第i个储备池神经元的输出连接权值，w_i′为修剪后的第i个储备池神经元的输出连接权值，ρ为修剪阈值，ρ为常数，0≤ρ1；通过设定修剪阈值，修剪贡献度较小的储备池神经元的输出连接权值，保留剩余输出连接权值；

第三步：建立出水总磷TP软测量的储备池网络模型

1)初始化储备池网络；储备池网络由输入层，储备池和输出层三部分组成，其中输入层神经元个数为K，储备池神经元个数为N，输出层神经元的个数为L，K，N和L为正整数；初始化网络输入矩阵Wⁱⁿ，反馈矩阵W^back和储备池内部连接权值矩阵W，为了保证储备池丰富的动态特性，储备池内部连接权值矩阵W的谱半径小于1；

2)将N_s个训练样本输入到储备池网络中，对网络进行训练；储备池网络的状态变量表示为x_l，l＝1,2,...100，根据公式(9)对状态变量进行更新

x(n)＝f(Wⁱⁿu(n)+Wx(n-1)+W^backy(n)) (9)

其中，f(·)为储备池状态变量的激活函数，u(n)表示输入变量，n＝1,2,...N_s,并收集储备池网络的状态变量，生成状态矩阵M：

计算储备池网络的输出变量y(n)，具体的计算公式如下：

y(n)＝f^out(W^outx(n)) (11)

其中，f^out(·)为网络输出的激活函数，n＝1,2,...N_s

3)分别计算网络的输出连接权值W^out和网络的检验误差E^v，计算公式如下：

W_out＝M⁺Y (12)

其中，M⁺表示状态矩阵M的Moore-Penrose逆矩阵；Y＝[y(1),y(2),…,y(N_s)]^T是储备池网络的输出，d(n)是目标输出变量，表示出水总磷TP真实值，y(n)是储备池网络的输出变量，表示储备池网络对出水总磷TP预测值；

4)优化储备池网络结构

根据公式(1)-(7)计算储备池内每一个神经元的贡献度，并根据公式(8)的修剪判别机制，修剪贡献度较小的储备池神经元的输出连接权值，得到优化后储备池网络的状态矩阵M_c＝[x₁,x₂,…,x_lprune],其中lprune为修剪后剩余的储备池神经元个数；

5)根据公式(12)对优化后储备池网络的输出连接权值进行调整，并根据公式(13)计算优化后网络的检验误差修剪终止条件如下：

其中，E_s表示优化前后网络检验误差的差值；若E_s≤0继续调整储备池网络结构，否则网络结构调整结束，得到最终的建立出水总磷TP软测量的储备池网络模型；

第四步：用测试数据对储备池网络进行测试

将测试样本中辅助变量数据作为训练好的基于贡献度优化后储备池网络的输入，网络的输出即为出水总磷TP的预测结果。