基于PLS局部拟合的SCR反应器信息快速获取方法技术

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

Rapid Information Acquisition Method for SCR Reactor Based on PLS Local Fitting

【技术实现步骤摘要】
基于PLS局部拟合的SCR反应器信息快速获取方法
本专利技术属于工业脱销
，具体涉及一种基于PLS局部拟合的SCR反应器信息快速获取方法。
技术介绍
由于国家新的排放政策的出台，对污染物的排放标准提出了新的要求。为了能够达到新的排放要求，各大电厂已基本完成脱硝系统改造任务。目前，电站锅炉中一般联合使用低NOx燃烧技术和烟气脱硝技术来达到超低排放。选择性催化还原法(selectivecatalyticreduction，SCR）烟气后处理技术以其脱硝效率高、技术成熟等优点，在各大电站脱硝系统中得到大规模的应用。SCR脱硝技术是指在催化剂和氧气存在的条件下，在一定温度范围内，通过氨作为还原剂，有选择的与烟气中的NOx反应生成无害的氮和水，从而除去烟气中的NOx污染物，保证烟气达标排放甚至污染物零排放。影响SCR脱硝效率的因素包括催化剂活性、浓度场分布、反应温度、空速、氨氮比等，还和喷氨格栅（ammoniainjectiongrid，AIG）处及反应器内流场、浓度场有关。现如今为了评估流场以及浓度场，许多学者采用CFD（计算流体动力学）对SCR反应器进行建模分析，通过数值计算的方式得到SCR催化剂层出口的NOx浓度分布和氨逃逸分布。该方法一般能够较好的得出一般性结论，但该种方法存在着计算时间过长、变工况运行时误差较大、不便于实时在线监测等缺点。
技术实现思路
：专利技术目的：为了克服现有技术中的不足，本专利技术提供一种基于PLS局部拟合的SCR反应器信息快速获取方法，在针对某一特定的SCR脱硝系统基础之上，通过CFD软件和试验台得出宽负荷工况下的SCR催化剂层的仿真数据和实验数据，并将其录入到SQLsever数据库中，当输入特定的入口参数之后，通过SQL语句在方差最小化原则的基础上提取得到距离输入的入口参数最近的局部拟合数据集，并通过微软.net平台上采用偏最小二乘法（partialleast-squaremethod，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的计算公式如下：偏最小二乘法的回归分析建模的具体步骤如下：(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=w1TX，u1是因变量集Y=(y1,...yp)T的线性组合；u1=v11y1+…+v1pyp=v1TY。为了回归分析的需要，需求t1和u1各自尽可能多地提取出所在变量组的变异信息且两者相关程度最大。(3.2.3)由两组变量集的标准化观测数据阵E0和F0，可以计算第一对成分的得分向量，记为和第一对成分t1和u1的协方差Cov(t1,u1)可用第一对成分的得分向量和的内积来计算，故而，以上两个要求可化为数学上的条件极值问题：w1Tw＝||w1||2＝1，v1Tv1＝||v1||2＝1公式六利用Lagrange乘数法，问题化为求单位向量w1和v1，使得最大。问题的求解只需通过计算m×m矩阵M＝E0TF0F0TE0的特征值与特征向量，且M的最大特征值为θ12，相应的单位特征向量就是所求的解w1，而v1可由w1计算得到(3.2.4)建立y1,y2,…yp，对t1的回归及x1,x2,…,xm对t1的回归。假定回归模型为：其中，a1＝(a11,…,a1m)T，β1＝(β11,...β1p)T分别是多对一的回归模型中的参数向量，E1和F1是残差阵。回归系数向量a1，β1的最小二乘估计为：称a1，β1为模型的效应负荷量。(3.2.5)用残差阵E1和F1代替E0和F0重复以上步骤。记则残差阵如果残差阵F1中的元素的绝对值近似为0，则认为用第一个成分建立的回归式精度已满足需要了，可以停止抽取成分，否则用残差阵E1和F1代替E0和F0重复以上步骤即得：w2＝(w21,…w2m)T；v2＝(v21,…,v2m)T分别为第二对成分的权数，而为第二对成分的得分向量。分别为X，Y的第二对成分的负荷量，这时有(3.2.6)设n×m数据阵E0的秩为r<＝min(n-1,m)，则存在r个成分t1,t2,…,tr，使得：把tk＝wk1x1+…+wkmxm(k＝1,2,…,r)，代入Y＝t1β1+...+trβr，即得P个因变量的偏最小二乘法回归方程式：yj＝aj1x1+…+ajmxm，(j＝1,2…,m)。有益效果：本专利技术本文档来自技高网...

【技术保护点】
1.基于PLS局部拟合的SCR反应器信息快速获取方法，其特征在于：所述步骤为：（1）建立自变量数据库和因变量数据库：根据SCR新鲜、劣化催化剂化学反应实验，得出与化学反应速率常数相关参数：温度、指前因子以及活化能；建立化学反应速率常数与相关参数数据库，并把数据库带入到fluent中，得到关于数值模拟的数据库，并最终获取自变量数据库和因变量数据库；（2）在任意输入一条参数后，通过SQL数据库中select语句快速查找并提取与所输入参数在方差最小化原则基础之上的距离输入参数最近的若干个样本数据，从而形成一个局部拟合数据集dataset；（3）利用偏最小二乘法对局部拟合标准化数据集进行建模，得到该局部拟合标准化数据集的偏最小二乘回归方程，并对输入参数进行回归，得到该输入参数下对应的输出值，将该输入参数与其对应的输出值分别加入到自变量数据库和因变量数据库。

【技术特征摘要】
1.基于PLS局部拟合的SCR反应器信息快速获取方法，其特征在于：所述步骤为：（1）建立自变量数据库和因变量数据库：根据SCR新鲜、劣化催化剂化学反应实验，得出与化学反应速率常数相关参数：温度、指前因子以及活化能；建立化学反应速率常数与相关参数数据库，并把数据库带入到fluent中，得到关于数值模拟的数据库，并最终获取自变量数据库和因变量数据库；（2）在任意输入一条参数后，通过SQL数据库中select语句快速查找并提取与所输入参数在方差最小化原则基础之上的距离输入参数最近的若干个样本数据，从而形成一个局部拟合数据集dataset；（3）利用偏最小二乘法对局部拟合标准化数据集进行建模，得到该局部拟合标准化数据集的偏最小二乘回归方程，并对输入参数进行回归，得到该输入参数下对应的输出值，将该输入参数与其对应的输出值分别加入到自变量数据库和因变量数据库。2.根据权利要求1所述的基于PLS局部拟合的SCR反应器信息快速获取方法，其特征在于：步骤（2）中，任意输入一条参数后，通过SQLsever数据库中select语句快速查找并提取与所输入参数在方差最小化原则基础之上的距离输入参数最近的8个样本数据；步骤（2）的步骤为：（2.1）获取标准化数据库：将步骤（2）中提取出每个队列数据进行标准化处理，减去自身的均值并除以其标准差，从而得到标准化数据库；（2.2）将所输入的参数根据标准化数据库进行同样的标准化，得到标准化输入参数；（2.3）建立局部拟合数据集dataset：将步骤（2.1）中快速提取出的数据建立一个局部拟合标准化数据集，利用偏最小二乘法（PLS）对局部拟合标准化数据集进行建模，进行步骤（3）。3.根据权利要求1所述的基于PLS局部拟合的SCR反应器信息快速获取方法，其特征在于：步骤（3）的具体步骤为：（3.1）设因变量组和自变量组的n次标准化观测数据阵分别记为F0、E0，F0是由输出参数出口NO浓度和氨逃逸率组成的因变量矩阵，E0是由输入参数入口烟气速度、氨氮摩尔比、入口烟气温度以及入口NO浓度组成的自变量矩阵，每一列代表一个样本，F0、E0的计算公式如下：（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与原自变量的回归方程式，即偏最小二乘回归方程式。4.根据权利要求3所述的基于PLS局部拟合的SCR反应器信息快速获取方法，其特征在于：步骤(3.2)中，偏最小二乘回归方程式的展开如下：(3.2.1)分别提取两变量组的第一对成分，并使之相关性达最大；(3.2.2)从两组变量分别提取第...

【专利技术属性】
技术研发人员：王虎，常峥峰，黄力，纵宇浩，张鑫，陈立志，司风琪，
申请(专利权)人：大唐南京环保科技有限责任公司，东南大学，
类型：发明
国别省市：江苏,32