[发明专利]一种电站锅炉NOx排放动态软测量方法

权利要求书

1.一种电站锅炉NO_x排放动态软测量方法，其特征在于：包括以下步骤：

(1)数据采集及预处理：确定模型输入变量和输出变量，从DCS系统中采集得到原始样本对原始样本采用归一化预处理方法，并将各工况下测得的NO_x浓度折算到6％氧量下的值；其中x_i∈R^m表示第i组模型输入样本，m为输入变量的个数，y_i∈R表示第i组模型输出样本，n为样本数量；

(2)初始化自适应粒子群算法：每个粒子p(i)包含m+3个变量，其中包括m个输入变量x_i对应的阶次d_i、1个输出变量y对应的阶次d_y以及支持向量回归机SVMR模型的惩罚因子C和核参数σ；

(3)将输入数据和输出数据整理成非线性自回归滑动平均模型框架NARMAX模型结构：以粒子p(i)包含的输入变量x_i的阶次d_i和输出变量y的阶次d_y，将输入数据和输出数据整理成NARMAX模型结构，然后计算得到软测量模型的输入X(t)和输出Y(t)如下：

X(t)=[x1(t-T),K,x1(t-d1T),K,xm(t-T),K,xm(t-dmT),y(t-T),K,y(t-dyT)]Y(t)=y(t)]]>

上式中，T指的是收集测点数据的采样周期；X(t)为采样时刻t对应的模型输入量，包含m个输入变量x_i在采样时刻t之前的连续d_i个采样周期的状态量，以及被测量y在采样时刻t之前的连续d_y个采样周期的状态量；Y(t)为采样时刻t对应的模型输出量，为被测量y在采样时刻t时的状态量；

(4)SVMR模型初始化：以粒子p(i)包含的惩罚因子C和核参数σ采用高斯核函数即RBF来设置SVMR模型：

Y(t)=Σi=1l(αi*-αi)exp{-||X(t)-Xisv||22σ2}+b]]>

上式中，α_i与b为支持向量回归机参数，为支持向量，l为支持向量的个数，σ为RBF核参数；

(5)进行SVMR训练和拟合，计算粒子个体适应度f_i，保留最优个体：以X(t)作为SVMR模型输入，Y(t)作为SVMR模型输出，将整理后的样本集的前3/4作为训练样本，后1/4作为预测样本，进行模型的训练和拟合；以预测样本偏差指标为适应度函数，计算粒子个体适应度f_i，并得到群体最优适应度为f_m以及保留其粒子位置Pm，那么粒子p(i)的适应度函数为：

f(pi)=1nΣj=1n(yj-y^j)2]]>

上式中，n为预测时间序列长度，y_j和分别为第j个时刻的实测值与模型预测值；

(6)寻优收敛判断：若f_m<f_e，其中f_e为期望收敛阈值，当k＝K_max时，则寻优终止，执行步骤(9)；否则，令k＝k+1，执行步骤(7)；其中,k为当前迭代次数，K_max为最大迭代次数，阶次的寻优范围为[d_min,d_max]；

(7)计算粒子群体平均适应度，将粒子群体分为三个子群：第k次迭代中粒子p(i)的适应度值为f_i^k，群体最优粒子的适应度值为f_m；粒子群平均适应度为将适应度值优于的粒子适应度求平均得到根据粒子个体适应度f_i^k将其分为局部寻优、平衡寻优、全局寻优三个子群，分别计算其对应的惯性权重w和学习因子c₁、c₂；

其中，根据粒子个体适应度f_i^k将粒子群分为三个子群的过程为：

(7.1)计算粒子群体平均适应度，然后设定第k次迭代中粒子p(i)的适应度值为f_i^k，群体最优粒子的适应度值为f_m，粒子群体平均适应度值为

(7.2)将适应度值优于的粒子适应度求平均值得到定义用来评价粒子群的早熟收敛程度，△越小说明粒子群体越趋于早熟收敛；

(7.3)最后根据粒子个体适应度f_i^k将其分为以下三类：

①若则表明该粒子为群体中较好的粒子，分为局部寻优子群，按下式取较小的惯性权重w，使其在小范围内变动，同时取较大的c₁和较小的c₂，以加强粒子个体自身的学习能力并强化粒子的局部寻优能力：

w=w0-(w0-wmin)·|fik-favgbkfm-favgbk|]]>

c1=c1min+(c1max-c1min)·|fik-favgbkfm-favgbk|]]>

c2=c2max-(c2max-c2min)·|fik-favgbkfm-favgbk|]]>

②若则表明该粒子为群体中一般的粒子，分为平衡寻优子群，同时具有良好的全局寻优能力和局部寻优能力，维持惯性权重w不变，并取c₁＝c₂＝1，平衡粒子自身认知和社会认知的学习能力；

③若则表明该粒子为群体中较差的粒子，分为全局寻优子群，取较大的w，以增大变化幅度，强化粒子群体的全局搜索能力，同时取较小的c₁和较大的c₂，强化粒子的社会学习能力，调整策略照下式进行：

w=1.5-11+k1·exp(-k2·Δ)]]>

c1=c1max-12(c1max-c1min)·|fik-fmfavgk-fm|]]>

c2=c2min+12(c2max-c2min)·|fik-fmfavgk-fm|]]>

其中，c_1max和c_2max分别为认知学习因子c₁和社会学习因子c₂设定的取值范围的最大值，c_1min和c_2min为其取值范围的最小值；w₀为惯性权重的初始设定值，w_min为惯性权重的设定最小值，k₁、k₂为惯性权重w的调节因子，k₁决定惯性权重w的搜索上限，k₂用于控制惯性权重w的收敛速度；

(8)粒子个体速度和位置更新：根据步骤(7)得到的惯性权重w和学习因子c₁、c₂，对粒子个体速度V_i和位置P_i进行计算更新；

(9)输出最优结果，完成动态软测量模型的建立：输出f_m对应的最优个体所包含的d_i、d_y以及惩罚因子C和核参数σ，并以其参数作为设定参数建立最终的动态软测量SVMR模型；

(10)动态软测量的实现：对于任意新采集的样本，首先对其进行预处理后，输入步骤(9)得到的动态软测量SVMR模型，便可得到相应的软测量输出。