估计整体煤气化联合循环（IGCC)厂中的变量的方法和系统
[0002] 本
发明用政府支持、根据由美国
能源部授予的合同号DE-FC26-07NT43094做出。因 此，美国政府在本发明中具有某些权利。
[0004] 该申请与同此同时提交并且通过引用而全部合并于此的题为"Model Predictive Control System And Method For Integrated Gasification Combined Cycle Power Generation(
整体煤气化联合循环发电的模型预测控制系统和方法)"的美国
专利申请号 13/053731(律师案号 250557-1)有关。
技术领域
[0005] 本发明一般涉及整体煤气化联合循环（IGCC)发电，并且更特定地涉及用于估计 IGCC
发电厂中的变量的方法和系统。
背景技术
[0006] 整体煤气化联合循环（IGCC)技术作为清洁和高效发电的有吸引
力的技术而持续 发展，例如可从大量含
碳物质（例如煤或其他相对低成本的
燃料)生成。在IGCC的前端是称 为煤气化的工艺，其是将燃料(例如，煤)转变成可燃
合成气体(合成气)流的部分
氧化工艺。 因为污染引起的排放(例如，S0 X、N0X、汞、微粒，等)可大致上在发生燃烧之前从合成气流去 除，IGCC是环境友好的。尽管IGCC技术固有地保持清洁和高效发电的明显可能性，还要利用 机会以对于增强的可靠性、可用性、效率和灵活性来改进IGCC发电。
[0007] 已知用于操作IGCC电厂的本技术趋于基于过分简单的控制程序，如可通过严格且 繁琐的操作者指南而传达给操作者，而不必设计成实现任何有意义的优化策略，例如可由 于用于监测并且控制IGCC厂的有限的在线信息而引起。例如，IGCC厂的煤气化段可经受相 对恶劣的操作环境并且因此有限的在线
传感器可用于监测和控制。
[0008] 还已知基于模型的估计实现可有助于估计厂变量。然而，存在因为在这样的估计 实现中可通常未适当地解释建模和/或感测不确定性而可能出现的挑战。
[0009] 鉴于前面的考虑，制定其中可将厂测量(如可通过传感器套件获得的）与基于模型 的估计组合用于估计厂变量的估计策略，可对于建模和/或感测不确定性适当地校正这些 厂变量而不增加任何相当大的计算负担并且同时实现相当大的估计准确度，这将是可期望 的。
发明内容
[0010] -般地，本发明的至少一些方面可由用于估计整体煤气化联合循环（IGCC)发电厂 中的变量的系统实现，该系统可包括传感器套件，其耦合来感测指示相应的厂输入变量和 厂输出变量的
信号。该系统可进一步包括
扩展卡尔曼
滤波器(EKF)，其耦合来接收感测的厂 输入变量并且包括厂的动态模型，用于在离散时间处生成多个厂状态估计和对这些状态估 计的协方差矩阵。抢先约束处理器可配置成抢先地将状态估计和协方差矩阵约束为无约束 违规。测量-校正处理器可配置成基于感测的厂输出变量的处理来校正受约束的状态估计 和受约束的协方差矩阵。该测量-校正处理器耦合以用校正的状态估计和校正的协方差矩 阵来更新动态模型。更新的动态模型可配置成估计起初未被传感器套件感测的至少一个厂 变量的值。
[0011] 本发明的至少一些额外的方面可由用于估计整体煤气化联合循环（IGCC)发电厂 中的变量的方法实现。该方法可包括下面的动作:使传感器套件耦合来感测指示相应的厂 输入变量和厂输出变量的信号；向扩展卡尔曼滤波器(EKF)供应感测的厂输入变量，该扩展 卡尔曼滤波器包括厂的动态模型，用于在离散时间处生成多个状态估计和对这些状态估计 的协方差矩阵;抢先地将状态估计和协方差矩阵约束为无约束违规;基于感测的厂输出变 量的处理来校正受约束的状态估计和受约束的协方差矩阵；用校正的状态估计和校正的协 方差矩阵来更新动态模型；以及用更新的动态模型来估计起初未被传感器套件感测的至少 一个厂变量。
[0012] 本发明的再另外的方面可由具有计算机可执行指令的非暂时性有形计算机可读 介质实现，这些计算机可执行指令在由处理器执行时配置成执行下面的动作：测量来自传 感器套件的信号，这些信号指示相应的厂输入变量和厂输出变量；向扩展卡尔曼滤波器供 应感测的厂输入变量，该扩展卡尔曼滤波器包括厂的动态模型，用于在离散时间处生成多 个状态估计和对这些状态估计的协方差矩阵;抢先地将状态估计和协方差矩阵约束为无约 束违规;基于感测的厂输出变量的处理来校正受约束的状态估计和受约束的协方差矩阵； 用校正的状态估计和校正的协方差矩阵来更新动态模型；以及用更新的动态模型来估计起 初未被传感器套件感测的至少一个厂变量。
附图说明
[0013] 当下列详细描述参照附图（其中类似的符号代表整个图中类似的部件）阅读时，本 发明的这些和其他特征、方面和优势将变得更好理解，其中：
[0014] 图1是实施本发明的方面、如可用于在整体煤气化联合循环（IGCC)电厂中执行在 线估计的示例估计系统的
框图表示。
[0015] 图2是IGCC厂的示例煤气化段(如可由动态模型模拟的、如可以是在图1中示出的 估计系统的一部分)的简化示意模型表示。
[0016] 图3是图示关于抢先地受约束的扩展卡尔曼滤波器(EKF)(如可以是在图1中示出 的估计系统的一部分)的示例处理动作的
流程图。
具体实施方式
[0017] 本发明的方面可在用于整体煤气化联合循环（IGCC)发电厂的在线感测系统的上 下文中应用。在一个示例
实施例中，可将厂变量的测量(如可通过传感器套件获得的）与基 于模型的估计(例如，扩展卡尔曼滤波器(ΕΚΠ 估计)组合用于估计厂变量，其中的一些起初 可未被传感器套件测量。
[0018] 如将由本领域内技术人员意识到的，已经证明EKF估计是估计动态厂的状态的强 大的计算工具。然而，在传统的EKF实现中，存在因为在这样的EKF实现中可通常未适当地解 释模型或信号不确定性而可能出现的挑战。例如，对厂状态的约束可因为这样的约束可能 未容易地被传统EKF实现的处理架构所适应而被忽略。这些约束可基于厂的现实世界的物 理和/或操作考虑并且因此约束违规可导致不准确的EKF估计或甚至可以导致厂实际操作 条件的错误评估。
[0019] 本发明的
发明人提出创新性且简洁的策略或抢先地约束EKF使得可在不增加任何 大量计算负担的情况下可实现EKF估计的准确性方面的大大改进。如本文使用的，"抢先地 约束"指执行适于在执行(卡尔曼滤波器)测量更新之前避免约束违规的处理动作。
[0020] 图1是实施本发明的方面的示例估计系统10的框图表示。在一个示例实施例中，系 统10可用于估计整体煤气化联合循环（IGCC)发电厂12中的参数。传感器套件14可耦合来感 测指示例如燃料输入(例如，煤浆）、再循环C0 2、氧气馈送等相应的厂输入变量16的信号。传 感器套件14可进一步耦合来感测指示例如合成气组成和
温度、合成气流率、碳转换等相应 的厂输出变量18的信号。
[0021] 估计系统10可进一步包括扩展卡尔曼滤波器(EKF)20,其包括根据本发明的方面 的抢先地受约束的EKF。厂的动态模型22耦合来接收感测（即，测量)的厂输入变量16并且在 离散时间处传播多个厂状态估计和对这些状态估计的协方差矩阵。将意识到在该卡尔曼循 环阶段，厂状态估计(例如，先验估计)可经受感测不准确性，例如在来自传感器套件14的测 量中存在的偏置和/或噪声和/或建模不准确性，其可在动态模型12中存在。
[0022]抢先约束处理器24配置成抢先地将传播的状态估计和协方差矩阵约束为无约束 违规。在一个示例实施例中，抢先约束处理器24生成受约束的状态估计和受约束的协方差 矩阵。测量校正处理器26配置成基于感测的厂输出变量18来校正受约束的状态估计和受约 束的协方差矩阵。在一个示例实施例中，这样的校正可牵涉将至少一个感测的厂输出变量 与至少一个对应的模型输出变量比较。
[0023]测量校正处理器26耦合以用校正的状态估计和校正的协方差矩阵来更新动态模 型22。在一个示例实施例中，校正的状态估计可包括对厂的动态状态、输入(例如，控制输 入）、输出变量(例如可包括虚拟输出（即，未被传感器套件14测量的厂输出变量））、测量变 量的偏置和模型参数的相应校正。在一个示例实施例中，更新动态模型可包括限定先前未 知的模型参数和/或更新随时间变化的模型参数，例如气化器动力学、
辐射合成气冷却器 (RSC)
结垢等。在一个示例实施例中，更新的动态模型可配置成估计至少一个厂变量的值， 其起初可未被传感器套件14感测（即，虚拟感测）。尽管动态模型22、抢先约束处理器24和测 量校正处理器26在图1中图示为独立处理器(例如，处理模
块），本领域内技术人员将容易意 识到这样的处理器可集成在公共处理器中。
[0024]图2是如可通过动态模型22而建模的示例煤气化段32的简化示意模型表示。如将 由本领域内技术人员容易意识到的，煤气化段32的示例模型组件可包括气化器模型组件 34、辐射合成气冷却器(RSC)模型组件36、RSC淬火模型组件38、高压(HP)汽包模型组件40和 洗涤器模型组件42。由气化器模型组件34接收的示例输入变量可指示燃料输入(例如，煤 浆）、再循环CO 2和氧气馈送。由气化器模型组件34供应的示例模型输出变量可以是气化器 出口合成气组成和温度、灰流和碳转换。
[0025]由RSC模型组件36接收的示例输入可以是来自气化器模型组件34的合成气输出 流，和来自HP汽包模型组件40的
水/汽流以例如计算合成气温度中的瞬态变化、RSC管道应 力的水平以及RSC管道中的水流中的汽分数。淬火模型组件38可接收来自RSC出口的相应合 成气和/或灰/渣流作为示例输入以例如计算淬火出口合成气流流率、组成和
焓。RSC模型组 件36可耦合于HP汽包模型组件40,以例如计算HP汽包与RSC管件之间的水/汽循环流率。洗 涤器模型组件42可从淬火模型组件38接收淬火合成气输出并且然后计算洗涤的合成气组 成和温度。将意识到前面的煤气化段32的模型表示应该在示例意义上而非限制性意义上解 释，因为本发明的方面既不限于对于IGCC厂的煤气化段的任何特定建模实现也不限于任何 特定煤气化设计。
[0026]在一个示例实施例中，动态模型22可包括适合于在线模型预测和优化的煤气化段 的降阶动态模型。在一个示例实施例中，降阶模型可能够比全阶模型大致上更快执行模拟 (例如，快了近似至少1000倍)同时维持相对高水平的准确性。将意识到在实际实施例中，降 阶模型的更快的模拟速度对于实时模型和感测和/或控制设计是可期望的。
[0027]图3是图示关于实施本发明的方面的扩展卡尔曼滤波器20(图1)的示例处理动作 的流程图。框50代表初始化动作，例如可允许状态估计以及对状态估计的协方差矩阵的相 应初始化。框52允许处理动态模型22-个时间步，如可允许状态估计和状态协方差矩阵的 相应更新。框52代表抢先约束动作，例如可允许基于协方差矩阵的sigma-点处理的状态估 计和协方差矩阵的相应抢先约束。框56代表测量校正动作，例如可牵涉使用感测的厂输出 变量来校正受约束的状态和受约束的协方差矩阵。如将由本领域内技术人员意识到的，EKF 是递归过程并且在一个示例实施例中，框58允许在
迭代地行进到框52中的下一个时间步以 及随后的处理动作(例如约束动作54和测量校正动作56)之前执行动态模型22的更新动作。 [0028] 数学
基础[0029]下文的描述将聚焦在实施本发明的方面的扩展卡尔曼滤波器(EKF)的基本数学基 础上。下文的描述不意成为关于EKF的数学论文，因为这样的信息容易地在文献中可获得并 且我们的意图只意在适当地突出实施本发明的方面的概念。
[0030] Xk代表状态向量（即，整体状态向量），其在一个示例实施例中可包括动态厂状态 Xk、输入Uk和偏置bk以及模型参数pk，其中k指示离散
采样时间。
[0031]等式 1
[0032]
[0033] P代表状态向量的协方差矩阵。假定状态向量的示例大小是ηχΧ1，协方差矩阵P的 大小将是nxXnx。矩阵Q代表建模不确定性，例如可由于建模不准确性而引起。矩阵R代表测 量不确定性，例如可由于传感器噪声和/或偏置引起，UB和LB分别代表状态向量的上下界 (例如，约束极限），EQ代表状态向量的相等类型的约束。
[0034] EKF初始化
[0035]
[0036] 如将由本领域内技术人员意识到的，标记Aljhj.指示初始状态向量并且可以基于 多至并且包括时间0的测量而读为初始化时间（例如，时间0)处的初始状态估计。一般，可以 基于多至并且包括时间k的测量而读为"时间k处的X估计"（即，后验估计）。相似地，μ..丨 可以基于直到时间k-1的测量而读为"时间k处的X估计"（即，先验估计）。#〇| O指示状态向 量的初始协方差矩阵。
[0037] EKF 传播
[0038]等式 2
[0039]等式2指示用于生成当前时间步的先验估计的动态模型22(图1)的示例处理动作。
[0040]等式 3
[0041] 等式3指示当前时间步的预测输出。
[0042]等式 4
[0043] 图4表示用于获得相应状态-空间矩阵A和C的动态模型(22)的线性化。
[0044]等式 5
[0045] 图5指示使用模型的线性化形式并且包括模型噪声的协方差矩阵的更新处理动 作。应该意识到在一个示例实施例中，等式2-5可被框52(图3)所包含。
[0046] 根据本发明的方面，本发明的发明人已经开发高效计算
算法，用于在题为"执行测 量更新"的段下执行处理动作之前抢先地约束EFK。
[0047]注意传统的EKF估计技术未在下面的题为"EKF的抢先约束"的段下执行处理动作。 即，传统的EKF估计将在题为"执行测量更新"的段下继续于处理动作，其将被本领域内技术 人员识别为卡尔曼滤波器循环的部分，其中测量更新通过该时间处的实际测量调整先验估 计来生成后验估计。参见例如等式11的有机结构。将意识到测量更新牵涉计算卡尔曼增益 Kk，其进而受到协方差估计的影响。参见例如等式12的有机结构。另一个步骤是计算后验协 方差矩阵。参见例如等式13的有机结构。将意识到约束违规将不利地影响前面的计算的准 确性。从而，通过执行要在下文描述的在题为"EKF的抢先约束"的段下的处理动作，有利地 确保前面的等式的计算尊重任何可应用的约束而不引入任何大量的计算负担。
[0048] EKF的抢先约束
[0049] 确定 Sigma-点（SP)
[0050]等式 6 [00511等式6指示复
制动作，其中在传播步骤中获得的状态向量Xk将被复制2nx次。
[0052]等式 7
[0053] 等式7指示关于每个状态估计的均值的增加和扣除一个或多个标准偏差来生成 sigma点矩阵。
[0054] 应用约束
[0055]
[0056] 对于所有i=]_，…，]^和j=]_，…，2nx 等式8
[0057] UB和LB分别代表等式8中的ηχΧ1个上下界向量。
[0058] 确定具有受约束的均值和协方差元素的矩阵
[0059]等式 9
[0060]等式9允许通过计算它相应列上的sigma点矩阵的平均值来确定受约束的状态向量。
[0061]
[0062]等式10允许确定受约束的协方差矩阵(例如，与无迹卡尔曼滤波器(UKF)相似）。在 一个示例实施例中，等式6-10可被框54(图3)所包含。[0063I执行测量更新
[0064]等式 11
[0005]等式11允许基于预测Yk与测量的Ykm输出之间的差来校正状态向量。
[0066]等式 12
[0067] Kk代表用在等式10中确定的受约束的协方差矩阵而计算的卡尔曼增益。
[0068]
[0069] 等式13允许在考虑实际厂测量效应的情况下更新状态协方差矩阵。在一个示例实 施例中，等式11-13可被框56(图3)所包含。
[0070] 虚拟感测
[0071]等式 14
[0072] 如在前面的论述中指出并且通过等式14而在数学上表示的，起初未测量的厂变量 输出一由于多种示例目的（例如厂控制、监测和/或诊断）中的任何目的而可期望它们一可 由具有更新的状态向量的处理动态模型22来估计。
[0073] 预期继"执行测量更新"处理之后，根据给定应用的需要，可以通过使用对于本领 域内技术人员将可用的技术来对受约束优化(例如，二次规划）问题求解而可选地执行状态 估计和协方差矩阵的进一步约束。对于这样的受约束优化的解的示例，参见Dan Simon和 Donald L · Simon 的题为 "Aircraft Turbofan Engine Health Estimation Using Constrained Kalman Filtering!；使用受约束的卡尔曼滤波的
飞行器祸轮
风扇
发动机健康 估计）"（ASME Conf. Proc. 2〇〇3,485 (2〇〇3))的论文。
[0074] 将意识到本文公开的发明性估计系统和方法的方面可使用任何适合的编程语言 或编程技术由任何适合的处理系统实现。系统可以采取
硬件实施例、
软件实施例或包括硬 件和软件元件两者的实施例的形式。在一个实施例中，系统可通过软件(例如，抢先地受约 束的EKF)和硬件(例如，处理器、传感器)实现，其可包括但不限于
固件、常驻软件、
微码等。 此外，系统的部件可以采取从计算机可使用或计算机可读介质（其提供程序代码以供计算 机或任何指令执行系统使用或连同它们一起使用）
访问的
计算机程序产品的形式。计算机 可读介质的示例可包括非暂时性有形计算机可读介质，例如
半导体或固态
存储器、磁带、可 移动计算机
软盘、
随机存取存储器(RAM)、
只读存储器(ROM)、硬磁盘和光盘。光盘的当前示 例包括光盘-只读存储器(⑶-ROM)、光盘-读/写(CD-R/W)和DVD。
接口显示器可以是平板、平 面控制板显示器、PDA或诸如此类。
[0075]在一个示例实施例中，适合于存储和/或执行程序代码的处理系统在一个示例中 可包括直接或通过
系统总线间接耦合于存储器元件的至少一个处理器。这些存储器元件可 以包括在程序代码的实际执行期间采用的本地存储器、
大容量存储装置和高速缓存存储 器，其提供至少一些程序代码的暂时性存储以便使在执行期间必须从大容量存储装置检索 代码的次数减少。输入/输出或I/O设备(其包括但不限于
键盘、显示器、指点装置，等)可以 直接或通过介入的I/O
控制器而耦合于系统。
网络适配器也可耦合于系统来使
数据处理系 统能够通过介入的私有或公共网络而变成耦合于其他
数据处理系统或远程
打印机或存储 装置。
[0076]尽管已经在本文图示并且描述本发明的某些特征，本领域内技术人员将想起许多
修改和改变。因此，要理解附上的
权利要求意在涵盖所有这样的修改和改变，它们落入本发 明的真正精神内。
