稳态追踪自矫正的联合循环机组排烟含氧量在线测量方法
【专利摘要】本发明公开了一种稳态追踪及自矫正的燃气?蒸汽联合循环机组排烟含氧量在线测量方法，包括如下步骤：利用网络交换机获取机组实时运行数据，将采集的各个测点数据按时间顺序进行队列存储；利用样本熵从各队列数据中筛选最靠近当前时刻的稳态时间段，并同时将该段时间内各测点数据进行平均处理，进行故障诊断，若筛选的测点数据判定没有问题，则进行下一步的排烟含氧量计算，同时导入历史数据库，若筛选的测点数据判定出现问题，则利用该诊断模型进行测点数据重构，再进行下一步烟气氧含量计算；烟气氧含量经迭代计算得到。本发明提供一种燃气?蒸汽联合循环机组排烟含氧量在线测量方法，可实现稳态追踪及自适应矫正计算。
【专利说明】
稳态追踪自矫正的联合循环机组排烟含氧量在线测量方法
技术领域
[0001] 本发明涉及一种稳态追踪自矫正的联合循环机组排烟含氧量在线测量方法，兼有 样本熵稳态判定和鲁棒输入神经网络的故障诊断模型属于数据挖掘和机器学习建模领域。
【背景技术】：
[0002] 在工业过程中，存在一些无法直接测量或测量有很大滞后的变量，需要通过软测 量技术建立模型进行估测。对于燃气机组或燃气-蒸汽联合循环机组而言，由于燃烧室内温 度较高(超过1〇〇〇 °C)，无法利用仪器直接测量燃烧室内温度或燃气透平初温来监测燃烧状 况，但燃机排烟含氧量是监测燃烧室燃烧状况以优化燃烧的重要指标。燃机排烟中氧气含 量与过量空气系数有确定(单值）的函数关系，且燃气-蒸汽连环循环机组中燃料特性相对 稳定，因此排烟含氧量实际反映的是过量空气系数或空燃比，若空燃比过大，则多余的空气 吸收能量导致排烟损失增大，而当空燃比过小时，则会发生燃料不完全燃焼，因此通过监测 燃机排烟含氧量使其保持在一定的范围内，可以保证燃烧室的燃烧效率，促使机组的正常 运行，其次考虑到燃机燃烧室存在材料温度极限，排烟含氧量监测的安全意义更为重大。
[0003] 目前计算排烟含氧量的方法主要有直接测量法和分别基于工艺机理分析、过程对 象动态数学模型、神经网络、支持向量机等软测量方法。直接测量法是利用氧气传感器直接 测量，氧气传感器主要有热磁式氧气传感器和氧化锆传感器两种，但直接测量存在测量滞 后大、误差大，投资大、寿命短及传感器易发生故障等缺点。
[0004] 相对于燃煤机组，燃气机组或燃气-蒸汽联合循环机组燃料为特性稳定的洁净能 源天然气，燃烧过程相对简单，因此可以利用基于工艺机理分析的软测量方法计算排烟含 氧量。本发明根据现场燃机功率、天然气体积流量、压气机进气温度和压损、燃机排气温度 和压损等已有测点数据计算得到燃机排烟含氧量。进一步，考虑到测点传感器潜在的漂移 问题，建立了燃机侧鲁棒输入神经网络模型，以此及时诊断传感器故障并重构故障测点数 据，确保排烟含氧量计算的输入数据正常。
【发明内容】
：
[0005] 发明目的：为了克服现有技术中的不足，本发明提供一种稳态追踪自矫正的联合 循环机组排烟含氧量在线测量方法。
[0006] 技术方案:为解决上述技术问题，本发明的一种稳态追踪自矫正的联合循环机组 排烟含氧量在线测量方法，所述步骤包括如下：
[0007] (1)现场DCS数据经由网络交换机输出至数据接口，间隔5s采集燃机功率Wgt、天然 气体积流量q f、压气机进口气温度h和压力P1、压气机出口温度t2和压力p2、燃机排气温度t 4 和压力P4、IGV开度α、大气压力pQ测点数据。将采集的各测点数据按时间顺序进行队列存储， 连续存储lh内的测点数据，其中间隔采样周期和连续存储时间可自定义。
[0008] (2)利用样本熵从各队列数据中筛选最靠近当前时刻的稳态时间段，并同时将该 段时间内各测点数据进行平均处理，经RITNN测点数据诊断模型进行故障诊断。RITNN测点 数据诊断模型是由测点稳态历史数据库内的数据样本训练得到的。
[0009] (2.1)将步骤（1)中每个队列数据进行标准化处理（即为归一化处理，是对原始数 据的线性变换，使结果值映射到[0，1]之间。），设定滑动窗体参数(包括滑动窗口数据串长 度1及连续两滑动窗间的数据间隔t)，按照参数将队列数据划分为若干组滑动窗体(即在序 列中放置一个固定长度为1的小窗口，随着时间的推移，小窗口以数据间隔t不断向后移动， 同时求解窗体内数据串的样本熵，最终得到一组样本熵序列。），分段计算每个滑动内各测 点数据的样本熵。
[0010] (2.3)设定稳态容忍度，当各测点计算出的样本熵值大于稳态容忍度时，判定该工 况下各测点数据为稳态数据，并在机组运行时进行长期筛选和保留稳态数据工作，建立测 点稳态历史数据库。
[0011] (2.4)训练测点稳态历史数据库内的数据样本，设定隐藏层数目和各层神经元个 数，建立鲁棒输入神经网络模型(RITNN)，作为测点数据诊断模型。
[0012] (3)若出现测点故障，则将诊断模型给出的故障测点重构数据作为该测点数值;若 无测点故障，则将该组样本导入历史数据库中，更新测点数据诊断模型建模的数据样本。
[0013] (4)从重构数据或无故障数据中提取燃机功率Wgt、天然气体积流量qf、压气机进口 温度tl、燃机排气温度t4和大气压力P0。
[0014] (5)进入排烟含氧量计算模块，根据质量、能量及元素守恒定律，经迭代计算得到 稳定工况下的排烟含氧量。
[0015] (5.1)初始化空气燃料质量比
[0016]首先初始化一个空燃比u，u = mair/mf，其中mail·和mf分别为空气和天然气的质量流 量，假设u的初始值为10。
[0017] (5.2)计算天然气质量流量。
[0018] mf = qf · pf
[0019] 式中，qf为天然气体积流量，天然气密度。
[0020] (5.3)计算压气机进口空气焓值。
[0021 ] Aj = + &l " ij + ??.λ "ti + d*, 'if + ?Ζ4 + 〇-s · ij5 kJ./ kg
[0022] 式中，ti为进口空气温度，匕为。的多项式函数，各项系数为&〇 = -4.77;&1=1.05; a2 = -2.07;a3 = 3.39;a4 = -l .51 ;a5=l .79。
[0023] (5.4)计算燃机出口烟气焓值。
[0024]由天然气燃烧的化学方程式：
[0026]计算排烟中各气体成分体积分数分别为：
[0029] 式中，％2、： XN2、％32、、XAr分别为空气中各组分(^、^、(^、^?、△^?的体积分 数;Mail·和Mf分别为空气和天然气的摩尔质量。
[0030] 则燃机排气焓值可表示为：
[0031 ] h_t = Vo, -h〇i (/,) +VN, 'K2 (^4) + -^0, 'Κ?2 (^) + -^ H.O * ^7H,0 〇4)+-VAr - ,，Α々4)
[0032]式中，'?〇2 (G ) ' 心:（b) ' 々CO: (G) ' 力H」0 (Q) 'hAr ( t4 )分别为 t4温度下〇2、N2、C〇2、 H20、Ar的焓值，均可表示为t4的多项式函数，即：
[0034]式中，各多项式的系数矩阵为：
[0036] (5.5)计算燃机功率。
[0037] 将整个燃气侧看作一个整体，根据能量守恒方程，得到燃机功率的计算值：
[0038] Wgtc=mair · hi-mgas · h4+mf · LHV · qcc
[0039] 式中，mair = u · mf ;mgas = mair+mf=(l+u) · mf ;LHV为天然气的低位热值;n。。为燃烧 效率。
[0040] (5.6)计算燃机功率测点值与计算值的误差。
[0041 ]燃机功率测点值Wgt与程序计算值WZ的相对误差为：
[0042] σ= |ffgtc-ffgt|/ffgt
[0043] 若σ彡0.001，则u = u+0.01，进入步骤(3.4)不断迭代计算，直至σ<〇·〇〇1，得到最 终结果u。
[0045] (6)输出并记录最靠近当前时刻的稳态时间及排烟含氧量值。
[0046] 本发明实现了过程数据中稳态数据的提取和稳态工况下排烟含氧量的测算，进一 步，考虑到测点传感器潜在的漂移问题，建立了燃机侧鲁棒输入神经网络模型，以此诊断传 感器故障并重构故障测点数据，可实现排烟含氧量稳态下自适应矫正计算，目的是为了获 得各指标间的静态耦合关系以此监测优化燃烧。
[0047] 有益效果:本发明相对于现有技术而言，具有以下优点：
[0048] (1)依据质量、能量及元素守恒定律来进行排烟含氧量测算，不仅准确率高，而且 滞后性较小，可较为真实地再现燃机排烟含氧量。
[0049] (2)计算并记录稳定工况下的排烟含氧量，便于分析各指标间的静态耦合关系，达 到监测和优化燃烧的效果。
[0050] (3)建立了燃机侧鲁棒输入神经网络模型，进行测点传感器的故障诊断及数据重 构，同时通过不断的更新样本数据库使得模型更加接近实际。
【附图说明】
[0051]图1是本发明的流程图。
[0052]图2是本实施例中燃机功率曲线图。
【具体实施方式】：
[0053]下面结合附图对本发明作更进一步的说明。
[0054]以某200MW燃气-蒸汽联合循环机组为例，采集SIS系统中2015年05月12日20: 30: 00至24:00:00间的数据，采集间隔为5s。本发明框架主要由数据采集、数据筛选、数据校验 和排烟含氧量计算等核心模块，详细流程图如图1所示：
[0055] 1)现场DCS数据经由网络交换机输出至数据接口，间隔5s采集燃机功率Wgt、天然气 体积流量qf、压气机进口气温度和压力P1、压气机出口温度t2和压力p 2、燃机排气温度t4和 压力P4、IGV开度α、大气压力p〇测点数据。
[0056] 2)将采集的各测点数据按时间顺序进行队列存储，连续存储lh内的测点数据。利 用样本熵从各队列数据中筛选最靠近当前时刻的稳态时间段，并同时将该段时间内各测点 数据进行平均处理，经RITNN测点数据诊断模型进行故障诊断。
[0057] 3)若出现测点故障，则将诊断模型给出的故障测点重构数据作为该测点数值;若 无测点故障，则将该组样本导入历史数据库中，更新测点数据诊断模型建模的数据样本。
[0058] 4)从重构数据或无故障数据中提取燃机功率、天然气体积流量、压气机进口温度、 燃机排气温度和大气压力。
[0059] 5)进入排烟含氧量计算模块，计算天然气质量流量mf和压气机进口空气焓值lu，初 始化空燃比u为10。
[0060] 6 )由天然气燃烧的化学方程式，计算排烟中各气体成分体积分数 .V〇2、_FCQ 2.、.VH:。和yAr，通过体积分数加权计算排烟洽值h4。根据整个燃机侧的能量守 恒，得到燃机功率的计算值wy。
[0061 ] 7)计算燃机功率测点值与计算值的相对误差σ，若σ多0.001，则u = u+0.01，进入步 骤6)迭代计算，直至〇<0.001，得到最终结果u(
算得排烟含氧量&2 °
[0062] 8)输出并记录最靠近当前时刻的稳态时间及排烟含氧量值。
[0063]计算实例：
[0064] 图2为2015年05月12日20:30 :00至24:00:00间的燃机功率，可以看出在此期间燃 机降负荷运行，四个不同时刻记录并计算的排烟含氧量如表1所示。通过本发明可追踪到上 一小时内机组最近的稳态时间，同时计算该段时间内燃机排烟含氧量。在23:35:45至24: 〇〇:〇〇间发现天然气体积流量测点传感器发生漂移，修正前后的结果如表2所示，体现了本 发明自适应矫正的计算能力。
[0065]表1不同时刻排烟含氧量
[0067]表2传感器滑移后自适应矫正计算
[0070]以上所述仅是本发明的优选实施方式，应当指出：对于本技术领域的普通技术人 员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应 视为本发明的保护范围。
【主权项】
1. 一种稳态追踪自矫正的联合循环机组排烟含氧量在线测量方法，其特征在于:所述 测量步骤为： (1) 现场DCS数据经由网络交换机输出至数据接口，间隔5s采集燃机功率Wgt、天然气体 积流量qf、压气机进口气温度和压力 P1、压气机出口温度t2和压力p2、燃机排气温度t4和压 力P4、IGV开度α、大气压力po测点数据，将采集的各测点数据按时间顺序进行队列存储，连续 存储lh内的测点数据，其中间隔采样周期和连续存储时间可自定义； (2) 利用样本熵从各队列数据中筛选最靠近当前时刻的稳态时间段，并同时将该段时 间内各测点数据进行平均处理，经RITNN测点数据诊断模型进行故障诊断，RITNN测点数据 诊断模型是由测点稳态历史数据库内的数据样本训练得到的； (3) 若出现测点故障，则将诊断模型给出的故障测点重构数据作为该测点数值;若无测 点故障，则将该组样本导入历史数据库中，更新测点数据诊断模型建模的数据样本； (4) 从重构数据或无故障数据中提取燃机功率Wgt、天然气体积流量qf、压气机进口温度 tl、燃机排气温度t4和大气压力P0; (5) 进入排烟含氧量计算模块，根据质量、能量及元素守恒定律，经迭代计算得到稳定 工况下的排烟含氧量； (6) 输出并记录最靠近当前时刻的稳态时间及排烟含氧量值。2. 根据权利要求1所述的稳态追踪自矫正的联合循环机组排烟含氧量在线测量方法， 其特征在于:所述步骤(2)具体步骤为： (2.1) 将步骤（1)中存储的每个队列数据进行标准化处理，设定滑动窗体参数，包括滑 动窗口数据串长度1及连续两滑动窗间的数据间隔t，按照滑动窗体参数将队列数据划分为 若干组滑动窗体，分段计算每个滑动内各测点数据的样本熵； (2.3) 设定稳态容忍度，当各测点计算出的样本熵值大于稳态容忍度时，判定该工况下 各测点数据为稳态数据，并在机组运行时进行长期筛选和保留稳态数据工作，建立测点稳 态历史数据库； (2.4) 训练测点稳态历史数据库内的数据样本，设定隐藏层数目和各层神经元个数，建 立鲁棒输入神经网络模型，作为测点数据诊断模型。3. 根据权利要求1所述的稳态追踪自矫正的联合循环机组排烟含氧量在线测量方法， 其特征在于:所述步骤(5)具体步骤为： (5.1) 初始化空气燃料质量比 首先初始化一个空燃比U，u=mair/mf，其中mair和mf分别为空气和天然气的质量流量，假 设u的初始值为10; (5.2) 计算天然气质量流量； mf = qf · Pf 式中，qf为天然气体积流量，Pf为天然气密度； (5 · 3)计算压气机讲口空气焓倌式中，1：1为进口空气温度，111为1：1的多项式函数，各项系数为3() = -4.77;31 = 1.05;32 =-2.07 ；a3 = 3.39 ；a4 = ~l. 51 ；a5 = 1.79〇 (5.4) 计算燃机出口烟气焓值 由天然气燃烧的化学方程式：计算排烟中各气体成分体积分数分别为：> 式中，λ'0.2、％:、XC02、λΗ 20、ΧΑΓ分别为空气中各组分〇2、他、(1)2、112〇、41'的体 积分数;Mail·和Mf分别为空气和天然气的摩尔质量； 则燃机排气焓值可表示为：式中，/7。:（匕），/7卜、:（0，、，（，4)，1^(七4)分别为七4温度下〇2、~工〇2、112〇、八『 的焓值，均可表示为t4的多项式函数，BP:式中，各多项式的系数矩阵为：(5.5) 计算燃机功率 将整个燃气侧看作一个整体，根据能量守恒方程，得到燃机功率的计算值： ffgtc=mair · hi-mgas · h4+mf · LHV · η〇〇 式中，mair = u · mf;mgas=mair+mf=(l+u) · mf;LHV为天然气的低位热值;η。。为燃烧效率； (5.6) 计算燃机功率测点值与计算值的误差 燃机功率测点值Wgt与程序计算值的相对误差为： 〇=|ffgtc-ffgt|/ffgt 若σ彡0.001，则u = u+0.01，进入步骤(3.4)不断迭代计算，直至σ<〇. 〇〇1，得到最终结 果u; (5.7) 根据公式计算得到排烟含氧量y〇2。
【文档编号】G05B13/04GK105867132SQ201610247071
【公开日】2016年8月17日
【申请日】2016年4月19日
【发明人】司风琪, 孔羽, 黄志军, 周建新
【申请人】东南大学, 大唐苏州热电有限责任公司