CN109766577B - 基于pls局部拟合的scr反应器信息快速获取方法 - Google Patents

基于pls局部拟合的scr反应器信息快速获取方法 Download PDF

Info

Publication number
CN109766577B
CN109766577B CN201811492297.2A CN201811492297A CN109766577B CN 109766577 B CN109766577 B CN 109766577B CN 201811492297 A CN201811492297 A CN 201811492297A CN 109766577 B CN109766577 B CN 109766577B
Authority
CN
China
Prior art keywords
database
regression
partial
data
variable
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.)
Active
Application number
CN201811492297.2A
Other languages
English (en)
Other versions
CN109766577A (zh
Inventor
王虎
常峥峰
黄力
纵宇浩
张鑫
陈立志
司风琪
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Southeast University
Datang Nanjing Environmental Protection Technology Co Ltd
Original Assignee
Southeast University
Datang Nanjing Environmental Protection Technology Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Southeast University, Datang Nanjing Environmental Protection Technology Co Ltd filed Critical Southeast University
Priority to CN201811492297.2A priority Critical patent/CN109766577B/zh
Publication of CN109766577A publication Critical patent/CN109766577A/zh
Application granted granted Critical
Publication of CN109766577B publication Critical patent/CN109766577B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种基于PLS局部拟合的SCR反应器信息快速获取方法，步骤包括：针对某一特定的SCR脱硝系统，通过CFD软件和试验台得出宽负荷工况下的SCR催化剂层的仿真数据和实验数据，并将其录入到SQLsever数据库中；当输入特定的入口参数之后，通过SQL语句在方差最小化原则的基础之上提取出距离输入参数最近的局部拟合数据集；并采用偏最小二乘法对其快速的进行局部拟合；并将其结果返回添加到原始数据库中，丰富并发展数据集。本发明提供的基于PLS局部拟合的SCR反应器信息快速获取方法，能够得到出口NO浓度和氨逃逸率，可快速得到催化剂层出口浓度场和流场的分布情况，给电厂运行人员提供合理的运行指导。

Description

基于PLS局部拟合的SCR反应器信息快速获取方法
技术领域
本发明属于工业脱销技术领域，具体涉及一种基于PLS局部拟合的SCR反应器信息快速获取方法。
背景技术
由于国家新的排放政策的出台，对污染物的排放标准提出了新的要求。为了能够达到新的排放要求，各大电厂已基本完成脱硝系统改造任务。目前，电站锅炉中一般联合使用低NOx燃烧技术和烟气脱硝技术来达到超低排放。选择性催化还原法(selectivecatalytic reduction，SCR)烟气后处理技术以其脱硝效率高、技术成熟等优点，在各大电站脱硝系统中得到大规模的应用。SCR脱硝技术是指在催化剂和氧气存在的条件下，在一定温度范围内，通过氨作为还原剂，有选择的与烟气中的NOx反应生成无害的氮和水，从而除去烟气中的NOx污染物，保证烟气达标排放甚至污染物零排放。影响SCR脱硝效率的因素包括催化剂活性、浓度场分布、反应温度、空速、氨氮比等，还和喷氨格栅(ammoniainjectiongrid，AIG)处及反应器内流场、浓度场有关。
现如今为了评估流场以及浓度场，许多学者采用CFD(计算流体动力学)对SCR反应器进行建模分析，通过数值计算的方式得到SCR催化剂层出口的NOx浓度分布和氨逃逸分布。该方法一般能够较好的得出一般性结论，但该种方法存在着计算时间过长、变工况运行时误差较大、不便于实时在线监测等缺点。
发明内容：
发明目的：为了克服现有技术中的不足，本发明提供一种基于PLS局部拟合的SCR反应器信息快速获取方法，在针对某一特定的SCR脱硝系统基础之上，通过CFD软件和试验台得出宽负荷工况下的SCR催化剂层的仿真数据和实验数据，并将其录入到SQLsever数据库中，当输入特定的入口参数之后，通过SQL语句在方差最小化原则的基础上提取得到距离输入的入口参数最近的局部拟合数据集，并通过微软.net平台上采用偏最小二乘法(partial least-square method，PLS)对其快速的进行局部拟合，并将其结果返回添加到原始数据库中，丰富并发展数据集，通过拟合能够得到催化剂层出口流场的分布情况，从而给电站脱硝系统运行人员提供运行指导意见。
技术方案：为解决上述技术问题，本发明提供一种基于PLS局部拟合的SCR反应器信息快速获取方法，包括如下步骤：
(1)建立自变量数据库和因变量数据库：根据SCR新鲜、劣化催化剂化学反应实验，得出与化学反应速率常数相关参数：温度、指前因子以及活化能；建立化学反应速率常数与相关参数数据库，并把数据库以及SCR脱硝反应模型带入到fluent中，得到关于数值模拟的数据库，并最终获取自变量数据库和因变量数据库，测量的点包括入口烟气速度、氨氮摩尔比NH3/NO、入口烟气温度、入口NO浓度、出口NO浓度、氨逃逸率等测点数据，其中，自变量数据库包括入口烟气速度、氨氮摩尔比NH3/NO、入口烟气温度和入口NO浓度，因变量数据库包括出口NO浓度和氨逃逸率，并将其录入到SQL数据库中。
(2)任意输入一条参数后，通过SQLsever数据库中select语句快速查找并提取与所输入参数在方差最小化原则基础之上的距离输入参数最近的8个样本数据；
(2.1)获取标准化数据库：将步骤(2)中提取出每个队列数据进行标准化处理，减去自身的均值并除以其标准差，从而得到标准化数据库；
(2.2)将所输入的参数根据标准化数据库进行同样的标准化，得到标准化输入参数；
(2.3)建立局部拟合数据集dataset：将步骤(2.1)中快速提取出的数据建立一个局部拟合标准化数据集，并利用微软的.net技术对其进行进一步的处理；
(2.4)利用微软.net技术中的偏最小二乘法(PLS)对局部拟合标准化数据集进行建模，并对输入参数进行回归，得到该输入参数下对应的输出值，将该输入参数与其对应的输出值分别加入到自变量数据库和因变量数据库，即进行步骤(3)。
(3)采用偏最小二乘法(PLS)对所提取出的局部拟合数据集进行建模：
(3.1)设因变量组和自变量组的n次标准化观测数据阵分别记为F0、E0，F0是由输出参数出口NO浓度和氨逃逸率组成的因变量矩阵，E0是由输入参数入口烟气速度、氨氮摩尔比、入口烟气温度以及入口NO浓度组成的自变量矩阵(其中每一列代表一个样本)，F0、E0的计算公式如下：
Figure GDA0001968719530000031
Figure GDA0001968719530000032
偏最小二乘法的回归分析建模的具体步骤如下：
(3.2)偏最小二乘法的基本做法是首先在自变量集中提出第一成分t1，t1是x1,x2,…xm的线性组合，且尽可能多地提取原自变量集中的变异信息；同时在因变量集中也提取第一成分u1，并要求t1和u1相关程度达到最大；然后建立因变量y1,y2,…yp与t1的回归，如果回归方程已达满意的精度，则算法中止；若最终对自变量集提取r个成分t1,t2,…tr，偏最小二乘回归将通过建立y1,y2,…yp与t1,t2,…tr的回归式，然后再表示为y1,y2,…yp与原自变量的回归方程式，具体展开如下所示：
(3.2.1)分别提取两变量组的第一对成分，并使之相关性达最大。
(3.2.2)假设从两组变量分别提取第一对t1和u1，t1是自变量集X＝(x1,…xm)T的线性组合；t1＝w11x1+…+w1mxm＝w1 TX，u1是因变量集Y＝(y1,...yp)T的线性组合；u1＝v11y1+…+v1pyp＝v1 TY。为了回归分析的需要，需求t1和u1各自尽可能多地提取出所在变量组的变异信息且两者相关程度最大。
(3.2.3)由两组变量集的标准化观测数据阵E0和F0，可以计算第一对成分的得分向量，记为
Figure GDA0001968719530000033
和/>
Figure GDA0001968719530000037
Figure GDA0001968719530000035
Figure GDA0001968719530000036
第一对成分t1和u1的协方差Cov(t1,u1)可用第一对成分的得分向量
Figure GDA0001968719530000041
和/>
Figure GDA0001968719530000042
的内积来计算，故而，以上两个要求可化为数学上的条件极值问题：
Figure GDA0001968719530000043
w1 Tw＝||w1||2＝1，v1 Tv1＝||v1||2＝1 公式六
利用Lagrange乘数法，问题化为求单位向量w1和v1，使得
Figure GDA0001968719530000044
最大。问题的求解只需通过计算m×m矩阵M＝E0 TF0F0 TE0的特征值与特征向量，且M的最大特征值为θ1 2，相应的单位特征向量就是所求的解w1，而v1可由w1计算得到/>
Figure GDA0001968719530000045
(3.2.4)建立y1,y2,…yp，对t1的回归及x1,x2,…,xm对t1的回归。
假定回归模型为：
Figure GDA0001968719530000046
Figure GDA0001968719530000047
其中，a1＝(a11,…,a1m)T，β1＝(β11,...β1p)T分别是多对一的回归模型中的参数向量，E1和F1是残差阵。回归系数向量a1，β1的最小二乘估计为：
Figure GDA0001968719530000048
Figure GDA0001968719530000049
称a1，β1为模型的效应负荷量。
(3.2.5)用残差阵E1和F1代替E0和F0重复以上步骤。
记
Figure GDA00019687195300000410
则残差阵/>
Figure GDA00019687195300000411
如果残差阵F1中的元素的绝对值近似为0，则认为用第一个成分建立的回归式精度已满足需要了，可以停止抽取成分，否则用残差阵E1和F1代替E0和F0重复以上步骤即得：w2＝(w21,…w2m)T；v2＝(v21,…,v2m)T分别为第二对成分的权数，而/>
Figure GDA00019687195300000412
为第二对成分的得分向量。
Figure GDA0001968719530000051
分别为X，Y的第二对成分的负荷量，这时有
Figure GDA0001968719530000052
Figure GDA0001968719530000053
(3.2.6)设n×m数据阵E0的秩为r<＝min(n-1,m)，则存在r个成分t1,t2,…,tr，使得：
Figure GDA0001968719530000054
Figure GDA0001968719530000055
把tk＝wk1x1+…+wkmxm(k＝1,2,…,r)，代入Y＝t1β1+...+trβr，即得P个因变量的偏最小二乘法回归方程式：yj＝aj1x1+…+ajmxm，(j＝1,2…,m)。
有益效果：本发明相对于现有技术而言，具有以下优点：
(1)依据CFD数值模拟的仿真数据与搭建试验台实验得出的测量数据建立起原始数据库模型，并将其分为自变量数据库、因变量数据库，能够建立起任何工况下的数据库，当随机输入任意工况参数时都能快速而准确的找到与其匹配的数据集。
(2)快速的从SQL数据库中提取出对应的自变量数据和因变量数据，并通过基于微软的.net技术采用偏最小二乘法PLS进行准确的建模，得出每一层催化剂出口的实时参数，同时可将所预测出的数据分别的添加到自变量数据库和因变量数据库中，丰富和发展了数据库的样本信息。
(3)通过实时预测催化剂层出口NO浓度和氨逃逸率，可以快速的得到催化剂层出口浓度场和流场的分布情况，从而给电厂运行人员提供合理的运行指导。
附图说明
图1是本发明的流程图；
图2是本发明的PLS算法流程图。
具体实施方式
下面结合附图对本发明作更进一步的说明。
如图1-2所示，一种基于PLS局部拟合的SCR反应器信息快速获取方法，包括以下步骤：
(1)建立自变量数据库和因变量数据库：根据SCR新鲜、劣化催化剂化学反应实验，得出与化学反应速率常数相关参数：温度、指前因子以及活化能。建立化学反应速率常数与相关参数数据库。并把数据库以及SCR脱硝反应模型带入到fluent中，得到关于数值模拟的数据库，并最终获取自变量数据库和因变量数据库；
(2)在任意输入一条参数后，通过SQL数据库中select语句快速查找并提取与所输入参数在方差最小化原则基础之上的距离输入参数最近的若干个样本数据，从而形成一个局部拟合数据集dataset；
(3)利用偏最小二乘法对局部拟合标准化数据集进行PLS建模，得到该局部拟合标准化数据集的偏最小二乘回归方程，并对输入参数进行回归，得到该输入参数下对应的出口参数输出值，将该输入参数与其对应的输出值分别加入到自变量数据库和因变量数据库。
进一步的，步骤(2)的具体步骤为：
(2.1)在输入参数后，通过SQL数据库中select语句提取得到距离输入样本最近的八个样本数据集；
(2.2)获取标准化数据库：将步骤(2.1)中提取的每个队列数据进行标准化处理，减去自身的均值并除以其标准差，从而得到标准化数据库；
(2.3)将所输入的参数根据标准化数据库进行同样的标准化，得到标准化输入参数；
(2.4)建立局部拟合数据集：将步骤(2.2)中获得的数据建立成一个局部拟合标准化数据集。
进一步的，步骤(3)的具体步骤为：
(3.1)设因变量组和自变量组的n次标准化观测数据阵分别记为F0、E0，F0是由输出参数出口NO浓度和氨逃逸率组成的因变量矩阵，E0是由输入参数入口烟气速度、氨氮摩尔比、入口烟气温度以及入口NO浓度组成的自变量矩阵，每一列代表一个样本，F0、E0的计算公式如下：
Figure GDA0001968719530000061
Figure GDA0001968719530000062
(3.2)偏最小二乘法的回归分析建模的具体步骤如下：
偏最小二乘法的基本做法是首先在自变量集中提出第一成分t1，t1是x1,x2,…xm的线性组合，且尽可能多地提取原自变量集中的变异信息；同时在因变量集中也提取第一成分u1，并要求t1和u1相关程度达到最大，然后建立因变量y1，y2，...yp与t1的回归，如果回归方程已达满意的精度，则算法中止，若最终对自变量集提取r个成分t1，t2，...tr，偏最小二乘回归将通过建立y1，y2，...yp与t1，t2，...tr的回归式，然后再表示为y1，y2，...yp与原自变量的回归方程式，即偏最小二乘回归方程式。
进一步的，步骤(3.2)中，偏最小二乘回归方程式的展开如下：
(3.2.1)分别提取两变量组的第一对成分，并使之相关性达最大；
(3.2.2)从两组变量分别提取第一对t1和u1，t1是自变量集X＝(x1，...xm)T的线性组合，t1＝w11x1+...+w1mxm＝w1 TX，u1是因变量集Y＝(y1，...yp)T的线性组合，u1＝v11y1+...+v1pyp＝v1 TY，为了回归分析的需要，需求t1和u1各自尽可能多地提取出所在变量组的变异信息且两者相关程度最大；
(3.2.3)由两组变量集的标准化观测数据阵E0和F0，可以计算第一对成分的得分向量，记为
Figure GDA0001968719530000071
和/>
Figure GDA0001968719530000072
Figure GDA0001968719530000073
Figure GDA0001968719530000074
第一对成分t1和u1的协方差Cov(t1，u1)可用第一对成分的得分向量
Figure GDA0001968719530000075
和/>
Figure GDA0001968719530000076
的内积来计算，以上两个要求可化为数学上的条件极值问题：
Figure GDA0001968719530000077
w1 Tw＝||w1||2＝1，v1 Tv1＝||v1||2＝1 公式六
利用Lagrange乘数法，问题化为求单位向量w1和v1，使得
Figure GDA0001968719530000078
最大，问题的求解只需通过计算m×m矩阵M＝E0 TF0F0 TE0的特征值与特征向量，且M的最大特征值为θ1 2，相应的单位特征向量就是所求的解w1，而v1可由w1计算得到
Figure GDA0001968719530000081
(3.2.4)建立y1,y2,…yp，对t1的回归及x1,x2,…,xm对t1的回归；
假定回归模型为：
Figure GDA0001968719530000082
Figure GDA0001968719530000083
其中，a1＝(a11,…,a1m)T、β1＝(β11,...β1p)T分别是多对一的回归模型中的参数向量，E1和F1是残差阵；回归系数向量a1、β1的最小二乘估计分别为：
Figure GDA0001968719530000084
Figure GDA0001968719530000085
称a1，β1为模型的效应负荷量；
(3.2.5)用残差阵E1和F1代替E0和F0重复以上步骤；
记
Figure GDA0001968719530000086
则残差阵/>
Figure GDA0001968719530000087
残差阵F1中的元素的绝对值近似为0，则认为用第一个成分建立的回归式精度已满足需要，可停止抽取成分，否则用残差阵E1和F1代替E0和F0重复以上步骤。
进一步的，步骤(3.2.5)中，用残差阵E1和F1代替E0和F0重复以上步骤即得到：w2＝(w21,…w2m)T；v2＝(v21,…,v2m)T分别为第二对成分的权数，
Figure GDA0001968719530000088
分别为第二对成分的得分向量；/>
Figure GDA0001968719530000089
分别为X，Y的第二对成分的负荷量，这时有如下公式：
Figure GDA00019687195300000810
Figure GDA00019687195300000811
设n×m数据阵E0的秩为r<＝min(n-1,m)，则存在r个成分t1,t2,…,tr，使得：
Figure GDA00019687195300000812
Figure GDA00019687195300000813
把tk＝wk1x1+…+wkmxm(k＝1，2，…，r)代入Y＝t1β1+...+trβr，即得P个因变量的偏最小二乘法回归方程式：yj＝aj1x1+…+ajmxm，(j＝1,2…,m)。
进一步的，步骤(1)中，根据SCR新鲜、劣化催化剂化学反应实验，得出与化学反应速率常数相关参数：温度、指前因子以及活化能。建立化学反应速率常数与相关参数数据库。并把数据库以及SCR脱硝反应模型带入到fluent中，得到关于数值模拟的数据库，并最终获取自变量数据库和因变量数据库，测量的点包括入口烟气速度、氨氮摩尔比NH3/NO、入口烟气温度、入口NO浓度、出口NO浓度、氨逃逸率等测点数据，其中，自变量数据库包括入口烟气速度、氨氮摩尔比NH3/NO、入口烟气温度和入口NO浓度，因变量数据库包括出口NO浓度和氨逃逸率，并将其录入到SQL数据库中。
实施例1
如图1-2所示，以某600MW燃煤机组的脱硝系统为例，建立其CFD数学模型并得出数值模拟数据，并通过试验台得到该机组催化剂层的单通道实验数据，且将其录入到SQL数据库中，步骤如下：
(1)建立自变量数据库和因变量数据库：根据SCR新鲜、劣化催化剂化学反应实验，得出与化学反应速率常数相关参数：温度、指前因子以及活化能；建立化学反应速率常数与相关参数数据库，并把数据库以及SCR脱硝反应模型带入到fluent中，得到关于数值模拟的数据库，并最终获取自变量数据库和因变量数据库，测量的数据包括入口烟气速度、氨氮摩尔比NH3/NO、入口烟气温度、入口NO浓度、出口NO浓度、氨逃逸率等测点数据，其中，自变量数据库包括入口烟气速度、氨氮摩尔比NH3/NO、入口烟气温度和入口NO浓度，因变量数据库包括出口NO浓度和氨逃逸率，并将其录入到SQL数据库(sqlsever数据库)中，总共1250条数据，测量得到的部分原始数据如表1所示。
表1
Figure GDA0001968719530000091
Figure GDA0001968719530000101
(2)在输入参数后，通过SQL数据库中select语句快速的查找与输入参数在方差最小化原则基础之上的距离输入参数最近的8个样本数据，如输入参数为入口烟气速度5.5m/s、氨氮摩尔比0.86、入口烟气温度312℃、入口NO浓度330mg/m3，则提取出数据为，表2为局部数据集dataset：
表2
Figure GDA0001968719530000102
(3)对每个队列数据进行标准化处理，即：减去自身的均值并除以其标准差，从而得到标准化数据库，即得到表3。
表3
Figure GDA0001968719530000103
(4)将输入参数进行如上的标准化步骤得到入口烟气速度0.25819、氨氮摩尔比-1.11837、入口烟气温度-0.77754、入口NO浓度0.7610592。
(5)借助微软的.net技术对其进行进一步的处理，利用偏最小二乘法(PLS)对局部拟合标准化数据集进行建模，该方法的建模思路如图2所示，并对输入参数进行回归预测，得到该输入参数下对应的输出值：出口NO浓度为321.45mg/m3，氨逃逸率为121.956。
本发明实现了在通过CFD数值模拟计算与试验台实验测量得出的数据后并建立数据库，通过SQL语句实现快速的在最小化方差的原则下提取出距离输入参数最短的局部拟合标准化数据集，并分别从因变量数据库和自变量数据库中提取出对应的数据，基于微软的.net平台并采用PLS模型进行预测出在特定输入参数的条件下能够快速的得出催化剂出口NO浓度以及氨逃逸率指标，从而能够快速而准确的得出催化剂层出口处的流场和浓度场分布情况，能够给电站脱硝系统运行人员提供运行意见。
以上所述仅是本发明的优选实施方式，应当指出：对于本技术领域的普通技术人员来说，在不脱离本发明的原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。

Claims (5)

1.基于PLS局部拟合的SCR反应器信息快速获取方法，其特征在于：所述步骤为：
(1)建立自变量数据库和因变量数据库：根据SCR新鲜、劣化催化剂化学反应实验，得出与化学反应速率常数相关参数：温度、指前因子以及活化能；建立化学反应速率常数与相关参数数据库，并把数据库以及SCR脱硝反应模型带入到fluent中，得到关于数值模拟的数据库，并最终获取自变量数据库和因变量数据库；
(2)在任意输入一条参数后，通过SQL数据库中select语句快速查找并提取与所输入参数在方差最小化原则基础之上的距离输入参数最近的若干个样本数据，从而形成一个局部拟合数据集dataset，具体步骤为：
(2.1)在输入参数后，通过SQL数据库中select语句提取得到距离输入样本最近的八个样本数据集；
(2.2)获取标准化数据库：将(2.1)提取出的每个队列数据进行标准化处理，减去自身的均值并除以其标准差，从而得到标准化数据库；
(2.3)将所输入的参数根据标准化数据库进行同样的标准化，得到标准化输入参数；
(2.4)建立局部拟合数据集dataset：将步骤(2.2)中获得的数据建立成一个局部拟合标准化数据集；
(3)利用偏最小二乘法对局部拟合标准化数据集进行建模，得到该局部拟合标准化数据集的偏最小二乘回归方程，并对输入参数进行回归，得到该输入参数下对应的输出值，将该输入参数与其对应的输出值分别加入到自变量数据库和因变量数据库。
2.根据权利要求1所述的基于PLS局部拟合的SCR反应器信息快速获取方法，其特征在于：步骤(3)的具体步骤为：
(3.1)设因变量组和自变量组的n次标准化观测数据阵分别记为F0、E0，F0是由输出参数出口NO浓度和氨逃逸率组成的因变量矩阵，E0是由输入参数入口烟气速度、氨氮摩尔比、入口烟气温度以及入口NO浓度组成的自变量矩阵，每一列代表一个样本，F0、E0的计算公式如下：
Figure FDA0004263485980000021
Figure FDA0004263485980000022
(3.2)偏最小二乘法的回归分析建模的具体步骤如下：
偏最小二乘法的基本做法是首先在自变量集中提出第一成分t1，t1是x1,x2,…xm的线性组合，且尽可能多地提取原自变量集中的变异信息；同时在因变量集中也提取第一成分u1，并要求t1和u1相关程度达到最大，然后建立因变量y1,y2,…yp与t1的回归，如果回归方程已达满意的精度，则算法中止，若最终对自变量集提取r个成分t1,t2,…tr，偏最小二乘回归将通过建立y1,y2,…yp与t1,t2,…tr的回归式，然后再表示为y1,y2,…yp与原自变量的回归方程式，即偏最小二乘回归方程式。
3.根据权利要求2所述的基于PLS局部拟合的SCR反应器信息快速获取方法，其特征在于：步骤(3.2)中，偏最小二乘回归方程式的展开如下：
(3.2.1)分别提取两变量组的第一对成分，并使之相关性达最大；
(3.2.2)从两组变量分别提取第一对t1和u1，t1是自变量集X＝(x1,…xm)T的线性组合，t1＝w11x1+…+w1mxm＝w1 TX，u1是因变量集Y＝(y1,...yp)T的线性组合，u1＝v11y1+…+v1pyp＝v1 TY，为了回归分析的需要，需求t1和u1各自尽可能多地提取出所在变量组的变异信息且两者相关程度最大；
(3.2.3)由两组变量集的标准化观测数据阵E0和F0，可以计算第一对成分的得分向量，记为
Figure FDA0004263485980000031
和/>
Figure FDA0004263485980000032
Figure FDA0004263485980000033
Figure FDA0004263485980000034
第一对成分t1和u1的协方差Cov(t1,u1)可用第一对成分的得分向量
Figure FDA0004263485980000035
和/>
Figure FDA0004263485980000036
的内积来计算，以上两个要求可化为数学上的条件极值问题：
Figure FDA0004263485980000037
w1 Tw＝||w1||2＝1,v1 Tv1＝||v1||2＝1 公式六
利用Lagrange乘数法，问题化为求单位向量w1和v1，使得
Figure FDA0004263485980000038
最大，问题的求解只需通过计算m×m矩阵M＝E0 TF0F0 TE0的特征值与特征向量，且M的最大特征值为θ1 2，相应的单位特征向量就是所求的解w1，而v1可由w1计算得到/>
Figure FDA0004263485980000041
(3.2.4)建立y1,y2,…yp，对t1的回归及x1,x2,…,xm对t1的回归；
假定回归模型为：
Figure FDA0004263485980000042
Figure FDA0004263485980000043
其中，a1＝(a11,…,a1m)T、β1＝(β11,...β1p)T分别是多对一的回归模型中的参数向量，E1和F1是残差阵；回归系数向量a1、β1的最小二乘估计分别为：
Figure FDA0004263485980000044
Figure FDA0004263485980000045
称a1，β1为模型的效应负荷量；
(3.2.5)用残差阵E1和F1代替E0和F0重复以上步骤；
记
Figure FDA0004263485980000046
则残差阵/>
Figure FDA0004263485980000047
残差阵F1中的元素的绝对值近似为0，则认为用第一个成分建立的回归式精度已满足需要，可停止抽取成分，否则用残差阵E1和F1代替E0和F0重复以上步骤。
4.根据权利要求2所述的基于PLS局部拟合的SCR反应器信息快速获取方法，其特征在于：步骤(3.2.5)中，用残差阵E1和F1代替E0和F0重复以上步骤即得到：
w2＝(w21,…w2m)T；v2＝(v21,…,v2m)T分别为第二对成分的权数，
Figure FDA0004263485980000048
Figure FDA00042634859800000410
分别为第二对成分的得分向量；/>
Figure FDA0004263485980000049
分别为X，Y的第二对成分的负荷量，这时有如下公式：
Figure FDA0004263485980000051
Figure FDA0004263485980000052
设n×m数据阵E0的秩为r<＝min(n-1,m)，则存在r个成分t1,t2,…,tr，使得：
Figure FDA0004263485980000053
Figure FDA0004263485980000054
把tk＝wk1x1+…+wkmxm，k＝1,2,…,r，代入Y＝t1β1+...+trβr，即得P个因变量的偏最小二乘法回归方程式：yj＝aj1x1+…+ajmxm，j＝1,2…,m。
5.根据权利要求2所述的基于PLS局部拟合的SCR反应器信息快速获取方法，其特征在于：步骤(1)中，根据SCR新鲜、劣化催化剂化学反应实验，得出与化学反应速率常数相关参数：温度、指前因子以及活化能；建立化学反应速率常数与相关参数数据库，并把数据库以及SCR脱硝反应模型带入到fluent中，得到关于数值模拟的数据库，并最终获取自变量数据库和因变量数据库，测量的点包括入口烟气速度、氨氮摩尔比NH3/NO、入口烟气温度、入口NO浓度、出口NO浓度、氨逃逸率等测点数据，其中，自变量数据库包括入口烟气速度、氨氮摩尔比NH3/NO、入口烟气温度和入口NO浓度，因变量数据库包括出口NO浓度和氨逃逸率，并将其录入到SQL数据库中。
CN201811492297.2A 2018-12-07 2018-12-07 基于pls局部拟合的scr反应器信息快速获取方法 Active CN109766577B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811492297.2A CN109766577B (zh) 2018-12-07 2018-12-07 基于pls局部拟合的scr反应器信息快速获取方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811492297.2A CN109766577B (zh) 2018-12-07 2018-12-07 基于pls局部拟合的scr反应器信息快速获取方法

Publications (2)

Publication Number Publication Date
CN109766577A CN109766577A (zh) 2019-05-17
CN109766577B true CN109766577B (zh) 2023-07-14

Family

ID=66451315

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811492297.2A Active CN109766577B (zh) 2018-12-07 2018-12-07 基于pls局部拟合的scr反应器信息快速获取方法

Country Status (1)

Country Link
CN (1) CN109766577B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112461995A (zh) * 2020-11-03 2021-03-09 西安热工研究院有限公司 一种预测火电厂氨逃逸的方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108804758A (zh) * 2018-05-01 2018-11-13 南京信息工程大学 基于数据拟合的稀土矿尾中稀土元素含量的规律研究方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108804758A (zh) * 2018-05-01 2018-11-13 南京信息工程大学 基于数据拟合的稀土矿尾中稀土元素含量的规律研究方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
一种简化递推偏最小二乘建模算法及其应用;刘强等;《北京航空航天大学学报》;20030830(第07期);全文 *
基于LabVIEW-化学计量学的化学分析系统;徐昊 等;《微型机与应用》;20121210(第23期);第4-7页 *

Also Published As

Publication number Publication date
CN109766577A (zh) 2019-05-17

Similar Documents

Publication Publication Date Title
CN111804146B (zh) 一种智能喷氨控制方法以及智能喷氨控制装置
CN101180590B (zh) 预测性排放物监测系统及方法
CN112085277B (zh) 一种基于机器学习的scr脱硝系统预测模型优化方法
CN102565274A (zh) 一种电厂scr脱硝系统中催化剂活性的修正方法
CN105971738A (zh) 一种燃气轮机组烟气流量在线监测的方法
Campbell Technology innovation & advancements for shell cansolv CO2 capture solvents
CN109499364B (zh) 一种基于数字镜像的催化剂辅助设计方法
CN110189800B (zh) 基于多粒度级联循环神经网络的炉氧含量软测量建模方法
CN113689917A (zh) 一种基于SCR出口NOx浓度的可视化喷氨优化方法及装置
CN116822380B (zh) 一种基于数字孪生在铜冶炼过程中尾气回收利用的协同优化方法
CN109766577B (zh) 基于pls局部拟合的scr反应器信息快速获取方法
CN114225662A (zh) 一种基于滞后模型的烟气脱硫脱硝优化控制方法
CN116189791A (zh) 一种scr系统氨氮双控喷氨优化方法、装置及存储介质
Hao et al. Prediction of nitrogen oxide emission concentration in cement production process: a method of deep belief network with clustering and time series
CN103529161A (zh) 选择性催化还原法脱硝系统催化剂活性的检测装置和方法
Shi et al. Modeling and simulation of marine SCR system based on Modelica
Dvoretskii et al. A new approach to the optimal design of industrial chemical-engineering apparatuses
Paul et al. Sensor network design for maximizing process efficiency: an algorithm and its application
CN105241823B (zh) 基于稀疏表示的火电厂烟气光谱定量分析方法
CN116128136A (zh) 一种基于LSO-CatBoost的燃煤电厂锅炉NOX排放预测方法
Li et al. Fault separation and detection algorithm based on Mason Young Tracy decomposition and Gaussian mixture models
Wang Model-based control of selective catalytic reduction systems
Han et al. A simulation study on NOx reduction efficiency in SCR catalysts utilizing a modern C3-CNN algorithm
Hadlocon et al. Statistical modelling of ammonia absorption in an acid spray scrubber
Tingting et al. Modeling on SCR process of a coal-fired boiler using LSSVM method

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant