本发明提供了一种基于改进SSA的燃煤锅炉减少NOx氮氧化物排放优化方法及装置，优化方法具体包括以下步骤：(1)采集燃煤锅炉的历史特征参量数据；(2)采用RReliefF算法和Pearson相关性分析获取最优影响特征；(3)确定优化目标函数以及可调特征的寻优区间；(4)采用Sin混沌产生个体初始位置；(5)选取适应度值优的个体作为发现者，更新发现者个体位置；(6)采用灵敏度‑信息素配合的方式改进跟随者选择发现者方式，更新跟随者个体位置；(7)更新边界个体位置，最终获得最优可调特征；(8)利用获得的最优可调特征，控制锅炉燃烧，以进一步降低NOx排放浓度。本发明提供的NOx排放优化方法及装置有效的降低了锅炉NOx的排放量，为电站锅炉实际运行过程中降低NOx的排放量提供决策支持，对电站节能减排，响应国家绿色环保政策具有重要
1.一种基于改进SSA的燃煤锅炉NOx排放优化方法，其特征在于，包括以下步骤：步骤1000：从电厂DCS系统中获取燃煤锅炉的特征参量历史数据与NOx排放量，形成第一数据集，所述第一数据集是一个n行、m列数据构成的二维矩阵X，其中，矩阵中行为采集到的n个样本数据，列为每个样本涉及的m‑1个特征及NOx排放量，将这n×m个数据构成矩阵X：可知xij(i＝1,2,…，n；j＝1,2,…,m)为第i个样本的第j个特征的取值；
步骤2000：通过归一化方法对步骤1000所述第一数据集中的数据进行归一化处理，形成第二数据集D，所述数据归一化方法采用Min‑Max归一化方法，归一化公式为：公式(1)中，MaxValue表示样本数据的最大值；MinValue表示样本数据的最小值；x表示样本原始数据；y表示归一化之后的数据；
步骤3000：确定锅炉燃烧过程中影响NOx排放量的若干特征，进一步包括步骤3100‑步骤3700：步骤3100：针对步骤2000所述的第二数据集，通过RReliefF算法计算各个特征参量与NOx之间的相对距离，通过相对距离对每个特征赋予权重，首先根据公式(2)计算相近样本中特征值A不同的概率：PdifA＝P(difvalue(A)|相近的样本) (2)公式(2)中，value(A)表示特征A，dif函数用于计算两个样本之间的距离并找到最近的邻近样本，PdifA表示相近样本中特征值A不同的概率，相近的样本表示两个样本在样本空间中相对距离最近；
然后根据公式(3)计算相近样本中NOx排放量不同的概率：
PdifC＝P(difNOx|相近的样本) (3)公式(3)中，PdifC表示相近样本中NOx排放量不同的概率；
最后根据条件概率得出公式(4)计算炉各个特征参量权重：
公式(4)中，W[A]表示锅炉各个特征参量权重，PdifC|difA表示在特征值相异的相近样本中NOx排放量不同的概率；
步骤3200：根据公式(5)计算各特征之间的Pearson相关系数：公式(5)中，i表示第i列特征，j表示第j列特征，为第i列特征样本均值， 为第j列样本特征均值，n为样本个数；
步骤3300：基于Bootstrap随机采样的思路，从步骤2000所述的第二数据集D中提取K个样本子集步骤3400：使用RReliefF算法对 的特征按照权重进行排序，并删除小于第一阈值的特征，得到K个不同的子集步骤3500：对 使用Pearson相关性分析，计算所有特征两两之间的Pearson相关系数，并取绝对值；
步骤3600：根据事先设定好的第二阈值，如果大于第二阈值，则删除在步骤3500所述的特征排序中靠后的一个特征，得到K个训练子集 通过这步来去除掉冗余数据；
步骤3700：所得到的结果进行汇总，输出出现次数最多的排序结果，得到对NOx排放量影响最大的若干个特征；
步骤4000：步骤3700所述的对NOx 排放量影响最大的若干个特征可分为可调特征和不可调特征，确定目标函数以及可调特征的寻优区间，进一步包括步骤4100‑步骤4300：步骤4100：确定将以锅炉NOx排放量最小为优化目标函数，优化目标函数见公式(6)：公式(6)中，X＝(x1,x2...xm)，i∈[1,m]， 为锅炉NOx的排放浓度，ai是第i个特征的可取值的最小值，bi是第i个特征可取值的最大值；
步骤4200：将步骤4100所述的锅炉运行过程中的可调特征作为优化过程的输入变量进行寻优，不可调特征在寻优过程中作为固定值保持不变；将锅炉运行可调特征定义为一个待优化向量x:x＝[x1,x2,x3,x4,x5…,xm] (7)公式(7)中，xi是各个可调参数，其中i＝1,2...m；
公式(8)中，xi是第i个可调参量，ai是第i个可调参量可取值范围的最小值，bi是第i个可调参量可取值范围的最大值，i∈[1,m]步骤5000：采用改进SSA算法对锅炉运行输入的可调特征进行寻优，以进一步降低NOx排放浓度，得到最优可调特征集，进一步包括步骤5100‑步骤5600：步骤5100：采用无限折叠Sin混沌初始化SSA种群以增加SSA的种群多样性，根据公式(9)产生N个可行解，每个可行解向量与所述可调特征一一对应：Xn＝sin(δ/xn),n＝0,1...,N (9)公式(9)中，xn表示第n个初始可行解，Xn表示第n个初始个体混沌映射后的可行解，δ∈(0,4]；‑1≤xn≤1且xn≠0；
步骤5200：根据发现者数量比例PD，选取适应度值优的作为发现者，根据式(10)新发现者个体位置：公式(10)中， 表示第t代的第i个发现者在第t代第j维中的位置信息，t代表当前迭代数次数，iterm表示最大的迭代次数，α∈(0,1]是一个随机数，Q是服从正态分布的随机数，R2∈(0,1]表示预警值，ST∈[0.5,1]表示安全值；
步骤5300：根据公式(11)(12)采用灵敏度‑信息素配合的方式改进跟随者选择发现者方式，跟随者选择发现者具体过程如下：信息素是与发现者的适应度值有比例关系的值，用于标记发现者，信息素计算如公式(11)所示：公式(11)中，P(i)是第i个发现者的信息素，i表示第i个发现者，f(i)是第i个发现者当前适应度值，fmin是数值最小的发现者适应度值，fmax是数值最大的发现者适应度值；
每个跟随者都有对信息素的灵敏度，灵敏度在优化过程中也有所不同，灵敏度计算如公式(12)所示：S(i)＝Smin+ΔSj (12)
公式(12)中，ΔSj＝(Smax‑Smin)·Rand(0,1)，Smax＝P(i)max，Smin＝P(i)min，P(i)max是当前数值最大的信息素，P(i)min是当前数值最小的信息素，Smax是当前数值最大的灵敏度，Smin是当前数值最小的灵敏度；
找出配合第j个跟随者灵敏度的发现者i：随机找出j，满足P(i)＝S(j)；
步骤5400：根据步骤4300对跟随者提出的约束，使用公式(13)更新跟随者个体位置：公式(13)中， 表示第t代的第i个跟随者在第j维中的位置信息，t代表当前迭代数次数， 是t+1代发现者所占据的最优位置， 则表示当前全局最差的发现者位置，d是维数，Q是满足标准正态分布的随机数，a∈(‑1,1)；
公式(14)中， 表示第t代的第i个边界个体在第j维中的位置信息，t代表当前迭代数次数， 是当前的全局最优位置；β是服从正态分布的随机数；K∈[‑1,1]是一个随机数，fi是当前麻雀个体的适应度值；fg和fw分别是当前全局最佳和最差的适应度值；ε是避免分母为零的最小常数；K表示麻雀移动的方向同时也是步长控制参数；
步骤5600：判断当前是否达到足够好的适应度值即最小NOx排放量或者达到最大迭代次数，如果满足，则程序终止，并输出此时最优的一组解x＝[x1,x2,x3,x4,x5…,xm]，从而获得最优可调特征；否则，迭代次数加1，跳转到步骤5200进行继续搜索。
2.一种采用如权利要求1所述的一种基于改进SSA的燃煤锅炉NOx氮氧化物排放优化方法的装置，其特征在于，所述装置包括：数据采集模块：用于采集步骤1000所述的燃煤锅炉的特征参量历史数据与NOx排放量；
数据处理模块：用于对第一数据集进行数据预处理，实现步骤2000所述的数据归一化以及步骤3000所述的数据降维处理；
优化模块：对锅炉运行中实时检测数据进行预处理后得到数据样本，通过步骤4000‑
5000采用改进SSA对锅炉的目标函数寻优获得最优解，利用最终的可调特征最优解，控制锅炉燃烧，以进一步降低NOx排放量。
一种基于改进SSA的燃煤锅炉NOx排放优化方法及装置
技术领域
[0001] 本发明涉及锅炉控制领域，特别涉及一种燃煤锅炉的NOx氮氧化物排放优化方法及装置。
背景技术
[0002] 火电厂大量的煤炭燃烧给环境带来了严重的破坏，主要污染源是煤炭在电厂锅炉内燃烧过程所产生的NOx等有害气体。它们与温室效应、酸雨光化学烟雾的形成以及地球臭氧层的破裂等都有直接的关系。目前，我国针对火电厂污染物的排放也颁布了相关的规定。这样一来，如果电厂锅炉的NOx排放过大，就会受到相关部门的严厉惩罚，直接影响到电厂的经济效益。因此，如何降低锅炉NOx排放量已是电厂急需解决的问题。
[0003] 锅炉燃烧系统是一个多变量、强耦合的非线性复杂系统，且燃烧过程存在着时变性和不确定性，致使依靠人工经验进行特征选择以及解析式优化的传统方法很难对锅炉进行有效的优化控制。近年来，由于人工智能技术的发展步伐加快，进而促进优化控制技术也进入一个崭新的阶段。基于人工智能的锅炉氮氧化物排放量优化技术不需要优化特征与优化目标之间精确的物理模型，可以很好地解决非线性复杂系统问题。因此，进一步深入研究基于人工智能的锅炉氮氧化物排放量优化方法，对于火电厂降低环境污染和经济损失具有重要的实际工程意义。
发明内容
[0004] 本发明解决的问题是提供一种基于改进SSA的燃煤锅炉NOx(氮氧化物)排放优化方法，采用RReliefF和Pearson相关性分析进行锅炉的特征选择，并采用改进的麻雀搜索算法(Sparrow Search Algorithm,SSA)对锅炉可调特征优化，能快速地获得影响锅炉氮氧化物排放的可调特征全局最优值。
[0005] 为实现上述目的，本发明提供了如下的技术方案：
[0006] 一种基于改进SSA的燃煤锅炉NOx氮氧化物排放优化方法，具体步骤如下：
[0007] 步骤1000：从电厂DCS系统中获取燃煤锅炉的特征参量历史数据与NOx排放量，形成第一数据集，所述第一数据集是一个n行、m列数据构成的二维矩阵X，其中，矩阵中行为采集到的n个样本数据，列为每个样本涉及的m‑1个特征及NOx排放量，将这n×m个数据构成矩阵X：
[0008]
[0009] 可知xij(i＝1,2,…，n；j＝1,2,…,m)为第i个样本的第j个特征的取值；
[0010] 步骤2000：通过归一化方法对步骤1000所述第一数据集中的数据进行归一化处理，形成第二数据集D，所述数据归一化方法采用Min‑Max归一化方法，归一化公式为：
[0011]
[0012] 公式(1)中，MaxValue表示样本数据的最大值；MinValue表示样本数据的最小值；x表示样本原始数据；y表示归一化之后的数据；
[0013] 步骤3000：确定锅炉燃烧过程中影响NOx排放量的若干特征，进一步包括步骤3100‑步骤3700：
[0014] 步骤3100：针对步骤2000所述的第二数据集，通过RReliefF算法计算各个特征参量与NOx之间的相对距离，通过相对距离对每个特征赋予权重，首先根据公式(2)计算相近样本中特征值A不同的概率：
[0015] PdifA＝P(difvalue(A)|相近的样本) (2)
[0016] 公式(2)中，value(A)表示特征A，dif函数用于计算两个样本之间的距离并找到最近的邻近样本，PdifA表示相近样本中特征值A不同的概率，相近的样本表示两个样本在样本空间中相对距离最近；
[0017] 然后根据公式(3)计算相近样本中NOx排放量不同的概率：
[0018] PdifC＝P(difNOx|相近的样本) (3)
[0019] 公式(3)中，PdifC表示相近样本中NOx排放量不同的概率；
[0020] 最后根据条件概率得出公式(4)计算燃煤锅炉各个特征参量权重：
[0021]
[0022] 公式(4)中，W[A]表示锅炉各个特征参量权重，PdifC|difA表示在特征值相异的相近样本中NOx排放量不同的概率；
[0023] 步骤3200：根据公式(5)计算各特征之间的Pearson相关系数：
[0024]
[0025] 公式(5)中，i表示第i列特征，j表示第j列特征，为第i列特征样本均值， 为第j列样本特征均值，n为样本个数；
[0026] 步骤3300：基于Bootstrap随机采样的思路，从步骤2000所述的第二数据集D中提取K个样本子集
[0027] 步骤3400：使用RReliefF算法对 的特征按照权重进行排序，并删除小于第一阈值的特征，得到K个不同的子集
[0028] 步骤3500：对 使用Pearson相关性分析，计算所有特征两两之间的Pearson相关系数，并取绝对值；
[0029] 步骤3600：根据事先设定好的第二阈值，如果大于第二阈值，则删除在步骤3500所述的特征排序中靠后的一个特征，得到K个训练子集 通过这步来去除掉冗余数据；
[0030] 步骤3700：所得到的结果进行汇总，输出出现次数最多的排序结果，得到对NOx排放量影响最大的若干个特征；
[0031] 步骤4000：步骤3700所述的对NOx排放量影响最大的若干个特征可分为可调特征和不可调特征，确定目标函数以及可调特征的寻优区间，进一步包括步骤4100‑步骤4300：
[0032] 步骤4100：确定将以锅炉NOx排放量最小为优化目标函数，优化目标函数见公式(6)：
[0033]
[0034] 公式(6)中，X＝(x1,x2...xm)，i∈[1,m]， 为锅炉NOx的排放浓度，ai是第i个特征的可取值的最小值，bi是第i个特征可取值的最大值；
[0035] 步骤4200：将步骤4100所述的锅炉运行过程中的可调特征作为优化过程的输入变量进行寻优，不可调特征在寻优过程中作为固定值保持不变；将锅炉运行可调特征定义为一个待优化向量x:
[0036] x＝[x1,x2,x3,x4,x5…,xm] (7)
[0037] 公式(7)中，xi是各个可调参数，其中i＝1,2...m；
[0038] 步骤4300：设置可调特征的寻优区间：
[0039]
[0040] 公式(8)中，xi是第i个可调参量，ai是第i个可调参量可取值范围的最小值，bi是第i个可调参量可取值范围的最大值，i∈[1,m]
[0041] 步骤5000：采用改进SSA算法对锅炉运行输入的可调特征进行寻优，以进一步降低NOx排放浓度，得到最优可调特征集，进一步包括步骤5100‑步骤5600：
[0042] 步骤5100：采用无限折叠Sin混沌初始化SSA种群以增加SSA的种群多样性，根据公式(9)产生N个可行解，每个可行解向量与所述可调特征一一对应：
[0043] Xn＝sin(δ/xn),n＝0,1…,N (9)
[0044] 公式(9)中，xn表示第n个初始可行解，Xn表示第n个初始个体混沌映射后的可行解，δ∈(0,4]；‑1≤xn≤1且xn≠0；
[0045] 步骤5200：根据发现者数量比例PD，选取适应度值优的作为发现者，根据式(10)新发现者个体位置：
[0046]
[0047] 公式(10)中， 表示第t代的第i个发现者在第t代第j维中的位置信息，t代表当前迭代数次数，iterm表示最大的迭代次数，α∈(0,1]是一个随机数，Q是服从正态分布的随机数，R2∈(0,1]表示预警值，ST∈[0.5,1]表示安全值；
[0048] 步骤5300：根据公式(11)(12)采用灵敏度‑信息素配合的方式改进跟随者选择发现者方式，跟随者选择发现者具体过程如下：信息素是与发现者的适应度值有比例关系的值，用于标记发现者，信息素计算如公式(11)所示：
[0049]
[0050] 公式(11)中，P(i)是第i个发现者的信息素，i表示第i个发现者，f(i)是第i个发现者当前适应度值，fmin是数值最小的发现者适应度值，fmax是数值最大的发现者适应度值；
[0051] 每个跟随者都有对信息素的灵敏度，灵敏度在优化过程中也有所不同，灵敏度计算如公式(12)所示：
[0052] S(i)＝Smin+ΔSj (12)
[0053] 公式(12)中，ΔSj＝(Smax‑Smin)·Rand(0,1)，Smax＝P(i)max，Smin＝P(i)min，P(i)max是当前数值最大的信息素，P(i)min是当前数值最小的信息素，Smax是当前数值最大的灵敏度，Smin是当前数值最小的灵敏度；
[0054] 找出配合第j个跟随者灵敏度的发现者i：随机找出j，满足P(i)＝S(j)；
[0055] 步骤5400：根据步骤4300对跟随者提出的约束，使用公式(13)更新跟随者个体位置：
[0056]
[0057] 公式(13)中， 表示第t代的第i个跟随者在第j维中的位置信息，t代表当前迭代数次数， 是t+1代发现者所占据的最优位置， 则表示当前全局最差的发现者位置，d是维数，Q是满足标准正态分布的随机数，a∈(‑1,1)；
[0058] 步骤5500：根据公式(14)更新边界个体位置：
[0059]
[0060] 公式(14)中， 表示第t代的第i个边界个体在第j维中的位置信息，t代表当前迭代数次数， 是当前的全局最优位置；β是服从正态分布的随机数；K∈[‑1,1]是一个随机数，fi是当前麻雀个体的适应度值；fg和fw分别是当前全局最佳和最差的适应度值；ε是避免分母为零的最小常数；K表示麻雀移动的方向同时也是步长控制参数；
[0061] 步骤5600：判断当前是否达到足够好的适应度值即最小NOx排放量或者达到最大迭代次数，如果满足，则程序终止，并输出此时最优的一组解x＝[x1,x2,x3,x4,x5…,xm]，从而获得最优可调特征；否则，迭代次数加1，跳转到步骤5200进行继续搜索。
[0062] 一种基于改进SSA的燃煤锅炉NOx氮氧化物排放优化方法的装置，其特征在于，所述装置包括：
[0063] 数据采集模块：用于采集步骤1000所述的燃煤锅炉的特征参量历史数据与NOx排放量；
[0064] 数据处理模块：用于对第一数据集进行数据预处理，实现步骤2000所述的数据归一化以及步骤3000所述的数据降维处理；
[0065] 优化模块：对锅炉运行中实时检测数据进行预处理后得到数据样本，通过步骤4000‑5000采用改进SSA对锅炉的目标函数寻优获得最优解，利用最终的可调特征最优解，控制锅炉燃烧，以进一步降低NOx排放量。
[0066] 与现有技术相比，本发明的有益效果是：
[0067] 本发明采用一种基于RReliefF和Pearson相关性分析的数据预处理方法，解决人工选择特征困难，采用一种基于Sin无限混沌策略和信息素灵敏度策略改进麻雀搜索算法，解决了因麻雀搜索算法全局寻找最优特征能力不强，而导致收敛速度慢、寻优精度不高的问题，为燃煤锅炉燃烧NOx排放量优化提供了一种基于人工智能的新方法，降低环境污染和经济损失。
附图说明
[0068] 图1是基于改进SSA的燃煤锅炉NOx氮氧化物排放优化方法流程图；
具体实施方式
[0069] 为了更清楚地理解本发明的上述方案，下面结合附图对本发明进行进一步的详细描述。需要说明的是，此处所描述的具体实施仅用于解释本申请，并不用于限定本申请。
[0070] 图1是本发明一种改进SSA燃煤锅炉NOx氮氧化物排放优化方法流程图，具体步骤如下：
[0071] 步骤1000：从电厂DCS系统中获取燃煤锅炉的特征参量历史数据与NOx排放量，形成第一数据集，所述第一数据集是一个n行、m列数据构成的二维矩阵X，其中，矩阵中行为采集到的n个样本数据，列为每个样本涉及的m‑1个特征及NOx排放量，将这n×m个数据构成矩阵X：
[0072]
[0073] 可知xij(i＝1,2,…，n；j＝1,2,…,m)为第i个样本的第j个特征的取值；
[0074] 步骤2000：通过归一化方法对步骤1000所述第一数据集中的数据进行归一化处理，形成第二数据集D，所述数据归一化方法采用Min‑Max归一化方法，归一化公式为：
[0075]
[0076] 公式(1)中，MaxValue表示样本数据的最大值；MinValue表示样本数据的最小值；x表示样本原始数据；y表示归一化之后的数据；
[0077] 步骤3000：确定燃煤锅炉燃烧过程中影响NOx排放量的若干特征，进一步包括步骤3100‑步骤3700：
[0078] 步骤3100：针对步骤2000所述的第二数据集，通过RReliefF算法计算各个特征参量与NOx之间的相对距离，通过相对距离对每个特征赋予权重，首先根据公式(2)计算相近样本中特征值A不同的概率：
[0079] PdifA＝P(difvalue(A)|相近的样本) (2)
[0080] 公式(2)中，value(A)表示特征A，dif函数用于计算两个样本之间的距离并找到最近的邻近样本，PdifA表示相近样本中特征值A不同的概率，相近的样本表示两个样本在样本空间中相对距离最近；
[0081] 然后根据公式(3)计算相近样本中NOx排放量不同的概率：
[0082] PdifC＝P(difNOx|相近的样本) (3)
[0083] 公式(3)中，PdifC表示相近样本中NOx排放量不同的概率；
[0084] 最后根据条件概率得出公式(4)计算炉各个特征参量权重：
[0085]
[0086] 公式(4)中，W[A]表示锅炉各个特征参量权重，PdifC|difA表示在特征值相异的相近样本中NOx排放量不同的概率；
[0087] 步骤3200：根据公式(5)计算各特征之间的Pearson相关系数：
[0088]
[0089] 公式(5)中，i表示第i列特征，j表示第j列特征，为第i列特征样本均值， 为第j列样本特征均值，n为样本个数；
[0090] 步骤3300：基于Bootstrap随机采样的思路，从步骤2000所述的第二数据集D中提取K个样本子集
[0091] 步骤3400：使用RReliefF算法对 的特征按照权重进行排序，并删除小于第一阈值的特征，得到K个不同的子集
[0092] 步骤3500：对 使用Pearson相关性分析，计算所有特征两两之间的Pearson相关系数，并取绝对值；
[0093] 步骤3600：根据事先设定好的第二阈值，如果大于第二阈值，则删除在步骤3500所述的特征排序中靠后的一个特征，得到K个训练子集 通过这步来去除掉冗余数据；
[0094] 步骤3700：所得到的结果进行汇总，输出出现次数最多的排序结果，得到对NOx排放量影响最大的若干个特征；
[0095] 步骤4000：步骤3700所述的对Nox排放量影响最大的若干个特征可分为可调特征和不可调特征，确定目标函数以及可调特征的寻优区间，进一步包括步骤4100‑步骤4300：
[0096] 步骤4100：确定将以锅炉NOx排放量最小为优化目标函数，优化目标函数见公式(6)：
[0097]
[0098] 公式(6)中，X＝(x1,x2...xm)，i∈[1,m]， 为锅炉NOx的排放浓度，ai是第i个特征的可取值的最小值，bi是第i个特征可取值的最大值；
[0099] 步骤4200：将步骤4100所述的锅炉运行过程中的可调特征作为优化过程的输入变量进行寻优，不可调特征在寻优过程中作为固定值保持不变；将锅炉运行可调特征定义为一个待优化向量x:
[0100] x＝[x1,x2,x3,x4,x5…,xm] (7)
[0101] 公式(7)中，xi是各个可调参数，其中i＝1,2...m；
[0102] 步骤4300：设置可调特征的寻优区间：
[0103]
[0104] 公式(8)中，xi是第i个可调参量，ai是第i个可调参量可取值范围的最小值，bi是第i个可调参量可取值范围的最大值，i∈[1,m]
[0105] 步骤5000：采用改进SSA算法对燃煤锅炉运行输入的可调特征进行寻优，以进一步降低NOx排放浓度，得到最优可调特征集，进一步包括步骤5100‑步骤5600：
[0106] 步骤5100：采用无限折叠Sin混沌初始化SSA种群以增加SSA的种群多样性，根据公式(9)产生N个可行解，每个可行解向量与所述可调特征一一对应：
[0107] Xn＝sin(δ/xn),n＝0,1...,N (9)
[0108] 公式(9)中，xn表示第n个初始可行解，Xn表示第n个初始个体混沌映射后的可行解，δ∈(0,4]；‑1≤xn≤1且xn≠0；
[0109] 步骤5200：根据发现者数量比例PD，选取适应度值优的作为发现者，根据公式(10)新发现者个体位置：
[0110]
[0111] 公式(10)中， 表示第t代的第i个发现者在第t代第j维中的位置信息，t代表当前迭代数次数，iterm表示最大的迭代次数，α∈(0,1]是一个随机数，Q是服从正态分布的随机数，R2∈(0,1]表示预警值，ST∈[0.5,1]表示安全值；
[0112] 步骤5300：根据公式(11)(12)采用灵敏度‑信息素配合的方式改进跟随者选择发现者方式，跟随者选择发现者具体过程如下：信息素是与发现者的适应度值有比例关系的值，用于标记发现者，信息素计算如公式(11)所示：
[0113]
[0114] 公式(11)中，P(i)是第i个发现者的信息素，i表示第i个发现者，f(i)是第i个发现者当前适应度值，fmin是数值最小的发现者适应度值，fmax是数值最大的发现者适应度值；
[0115] 每个跟随者都有对信息素的灵敏度，灵敏度在优化过程中也有所不同，灵敏度计算如公式(12)所示：
[0116] S(i)＝Smin+ΔSj (12)
[0117] 公式(12)中，ΔSj＝(Smax‑Smin)·Rand(0,1)，Smax＝P(i)max，Smin＝P(i)min，P(i)max是当前数值最大的信息素，P(i)min是当前数值最小的信息素，Smax是当前数值最大的灵敏度，Smin是当前数值最小的灵敏度；
[0118] 找出配合第j个跟随者灵敏度的发现者i：随机找出j，满足P(i)＝S(j)；
[0119] 步骤5400：根据步骤4300对跟随者提出的约束，使用公式(13)更新跟随者个体位置：
[0120]
[0121] 公式(13)中， 表示第t代的第i个跟随者在第j维中的位置信息，t代表当前迭代数次数， 是t+1代发现者所占据的最优位置， 则表示当前全局最差的发现者位置，d是维数，Q是满足标准正态分布的随机数，a∈(‑1,1)；
[0122] 步骤5500：根据公式(14)更新边界个体位置：
[0123]
[0124] 公式(14)中， 表示第t代的第i个边界个体在第j维中的位置信息，t代表当前迭代数次数， 是当前的全局最优位置；β是服从正态分布的随机数；K∈[‑1,1]是一个随机数，fi是当前麻雀个体的适应度值；fg和fw分别是当前全局最佳和最差的适应度值；ε是避免分母为零的最小常数；K表示麻雀移动的方向同时也是步长控制参数；
[0125] 步骤5600：判断当前是否达到足够好的适应度值即最小NOx排放量或者达到最大迭代次数，如果满足，则程序终止，并输出此时最优的一组解x＝[x1,x2,x3,x4,x5…,xm]，从而获得最优可调特征；否则，迭代次数加1，跳转到步骤5200进行继续搜索。
[0126] 一种基于改进SSA的燃煤锅炉NOx氮氧化物排放优化方法的装置，其特征在于，所述装置包括：
[0127] 数据采集模块：用于采集步骤1000所述的燃煤锅炉的特征参量历史数据与NOx排放量；
[0128] 数据处理模块：用于对第一数据集进行数据预处理，实现步骤2000所述的数据归一化以及步骤3000所述的数据降维处理；
[0129] 优化模块：对锅炉运行中实时检测数据进行预处理后得到数据样本，通过步骤4000‑5000采用改进SSA对锅炉的目标函数寻优获得最优解，利用最终的可调特征最优解，控制锅炉燃烧，以进一步降低NOx排放量。
[0130] 以上所述仅为本发明的具体实施例，并非因此限制本发明的专利范围，对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所做的任何修改、同等替换、改进等，均应包含在本发明的保护范围之内。
