一种基于有监督的线性动态系统模型的软测量建模方法
【专利摘要】本发明公开了一种基于有监督的线性动态系统模型的软测量建模方法，用于工业生产在噪声环境下的动态过程的软测量建模，并实现对于难以直接测量的质量变量的预测。本发明基于有监督的线性动态系统模型，建立了一个有效的软测量建模，并克服了工业生产中过程的动态性和采集数据的随机性特征。相比目前的其它方法，本发明建立的模型更加精确，模型的预测更加准确，使产品质量更加稳定；而且改善了软测量建模对过程知识的依赖性，更加有利于工业过程的自动化实施。
【专利说明】
一种基于有监督的线性动态系统模型的软测量建模方法
技术领域
[0001] 本发明属于工业生产过程软测量建模和应用领域，尤其涉及一种基于有监督的线 性动态系统模型的软测量建模方法。
【背景技术】
[0002] 随着科技的发展，工业生产过程越来越大型化、复杂化。现代工业过程中存在着许 多难以甚至无法用传感器直接测量的重要变量，比如产品的反应速率、产品的成分含量等 等。但是这些重要变量对于保证产品质量和提高生产效益有极其重要的作用，是工业生产 过程中必须严格监视和控制的参数。将传感器可以直接测量的或容易测得的变量称为过程 变量，将难以测量或无法测量的重要变量称为质量变量。那么软测量建模方法是指通过建 立工业生产中过程变量和质量变量之间的数学模型，实现利用过程变量预测质量变量的方 法。近年来，工业生产过程的软测量建模越来越得到工业界和学术界的广泛重视，成为研究 热点。
[0003] 传统的工业过程软测量建模方法主要是基于机理模型的方法。但是现代工业生产 过程的机理模型越来越难以获得，此时基于数据驱动的多元统计分析方法已经成为工业过 程软测量建模的主流方法，比如偏最小二乘(PLS)，概率主元回归(PPCR)。但是，传统的多元 统计分析方法大多没有考虑过程的动态性和过程变量的随机性，比如PLS。虽然PPCR考虑了 过程变量的随机性，但没有考虑过程的动态性。这些因素没有考虑会大大影响了软测量建 模的精确性和模型预测的准确性。相比之下，有监督的线性动态系统模型既考虑了动态性 又考虑了随机性，本发明采用该方法替代原有的多元统计分析方法实现工业过程的软测量 建模。
【发明内容】
[0004] 本发明的目的在于针对现有技术的不足，提供一种基于有监督的线性动态系统模 型的软测量建模方法。
[0005] 本发明的目的是通过以下技术方案来实现的：一种基于有监督的线性动态系统模 型的软测量建模方法，包括以下步骤：
[0006] (1)利用集散控制系统以及离线检测方法收集工业生产中的过程变量和质量变量 的数据，组成建模用的训练样本集 :X=[Xl，X2，~，XN]eRvx%PY=[ yi，y2，~，yN]eRLXNj 中，R为实数集，RVXN表示X满足V X N的二维分布，RlxN表示Y满足L X N的二维分布，V为过程变 量个数，L为质量变量个数，N为采样数据点数，将数据存入历史数据库。
[0007] (2)从历史数据库中调用训练样本集X和Y，分别对训练样本集X中的各个样本和Y中 的各个样本按照时间点方向进行排序，得到 y'2,…t，…n] jPy't分别为t时刻采集到的过程变量的训练样本和质量变 量的训练样本，t = 1，2，…，N。对每一个训练样本进行标准化处理，即使得各个过程变量和 质量变量的均值为〇,方差为1，得到的新数据矩阵分别为和
。冥为Y t经标准化处理后得到的样本，叉为y't经标准化处理 后得到的样本，t = l，2,…，N。
[0008] (3)根据训练样本集X和f采用期望最大化方法建立有监督的线性动态系统模型， 得到模型参数Θ。
[0009] (4)将建模数据和模型参数Θ存入历史数据库中备用。
[0010] (5)收集新的过程变量的在线数据： 为当前t时刻的在线过程变量数据。对其标准化处理，得到文_。
[0011] (6)采用基于有监督的线性动态系统模型的软测量建模方法，根据历史数据库中 的数据和过程变量的在线数据X ww ,预测对应时刻难以直接测量的质量变量
，实现工业生产上对质量变量的监控。
[0012] 本发明的有益效果是：
[0013] 本发明通过离线采集过程变量数据和质量变量数据，建立有监督的线性动态系统 模型，然后在线根据容易测量的过程变量数据预测对应的难以测量的质量变量数据。基于 有监督的线性动态系统模型的学习方法，对于工业过程的随机性和动态性等数据特性，具 有非常强的建模能力。相比目前的其它软测量建模方法，本发明把握了更多的工业过程的 数据特性，所以建立的模型更加精确，模型的预测更加准确，最终提高了工业过程对质量变 量的监控效果，使生产的产品质量更加稳定。
【附图说明】
[0014] 图1是基于有监督的线性动态系统模型的脱丁烷塔丁烷含量的在线软测量结果示 意图；
[0015] 图2是基于偏最小二乘(PLS)模型的脱丁烷塔丁烷含量的在线软测量结果示意图； [0016]图3是基于概率主元回归(PPCR)模型的脱丁烷塔丁烷含量的在线软测量结果示意 图。
【具体实施方式】
[0017] 下面结合附图和具体实施例对本发明进行详细说明。
[0018] 本发明给出一种基于有监督的线性动态系统模型的软测量建模方法，该方法针对 工业过程的软测量建模问题，首先利用集散控制系统以及离线检测方法收集工业生产中的 过程变量和质量变量的数据，然后建立有监督的线性动态系统模型，把所有建模数据和模 型参数存入数据库中备用。预测在线质量变量数据的时候，首先利用前向滤波方法计算出 对应的隐变量数据，然后根据模型参数预测难以直接测量的质量变量数据。
[0019] 本发明采用的技术方案的主要步骤如下：
[0020] 第一步:利用集散控制系统以及离线检测方法收集工业生产中的过程变量和质量 变量的数据，组成建模用的训练样本集:X=[ X1，X2r_，XN]eRvxqPY=[yi， y2r_，yN]eRL XN，其中，R为实数集，RVXN表示X满足V X N的二维分布，RlxN表示Y满足L X N的二维分布，V为过 程变量个数，L为质量变量个数，N为采样数据点数，将数据存入历史数据库；
[0021] 第二步:从历史数据库中调用训练样本集X和Y，分别对训练样本集X中的各个样本 和γ中的各个样本按照时间点方向进行排序，得到
[/! y 2，…y t，…，y'n] *和/1分别为t时刻采集到的过程变量的训练样本和质量 变量的训练样本，t = l, 2，···，N。对每一个训练样本进行标准化处理，即使得各个过程变量 和质量变量的均值为〇，方差为1，得到的新数据矩阵分别为！ = ??,…Λ7；，…e fxV 和? = ,γ2，…，?,,…,?#] e i?ixAr。I为V t经标准化处理后得到的样本，瓦为/ t经标准化处 理后得到的样本，t = l，2,…，N;
[0022] 在历史数据库中对采集到的过程数据进行预处理，剔除野值点和明显的粗糙误差 数据。
[0023] 第三步:根据训练样本集又和?采用期望最大化方法建立有监督的线性动态系统 模型，得到模型参数Θ;
[0024] 对于训练样本集X和f,采用期望最大化方法求出有监督的线性动态系统模型参 数0 = {A，P，C，Ση，Σβ, Σ?，μπ，Σπ}，其中AERHXH为传递矩阵，Η为隐变量个数;PeR VXH为映 射矩阵，CeRW为回归矩阵；ΣηαΗΧΗ为隐空间噪声的方差，Σ#Κ νχν为过程变量噪声的方 差，为质量变量噪声的方差;μπΕΚΗΧ4ΡΣ πΕΚΗΧΗ分别为服从高斯分布的初始时刻 隐变量heRa1的均值和方差。用期望最大化方法建模的具体实现步骤如下所示：
[0025] (1)设置初始的模型参数Θ为0 = {A，P，C，Ση，Σβ，Σ?，μπ，Σπ}
[0026] (2)求期望:在当前模型参数θ下，根据训练样本集X和f计算每个时刻隐变量的平 滑均值g teRHX1、方差GteRHXH，t = 1，2，…，N;相邻时刻隐变量的协方差 <h,h丨:p,?)e沪xi^V=t+l，且t，=l，2，…，N。其中E(·)表示括号中变量的均值，h t为t 时刻Xt，yt对应的隐变量。具体实现步骤如下所示：
[0027] (2.1)通过前向滤波方法可以得到每个时刻隐变量的滤波均值ft e RHX1和方差Ft ERHXH，t = l，2，."，N 如下：
[0032] (2.2)通过后向平滑方法可以得到每个时刻隐变量的平滑均值gt eRHX1，方差Gt e RHXH，t = l，2，...，N;相邻时刻隐变量的协方差，= / + l,且t，=1， 2，…，N。如下：
[0033] gt = FtAT (AFtAT+ Σ n) (gt^ -Aft) +f t (18)
[0034] Gt = FtAT(AFtAT+ Σn)[FtAT(AFtAT+ Σn)]T+Ft-FtAT(AF tAT+ Σn)^AFt (19)
[0035] ￡(h,h,r |Χ,?) = Γ;ΑΓ (AF,Ar + Σ" )_1 G, +g,g^ (20)
[0036] 其中当时刻 t = T 时，gT=fr，GT=FT。
[0037] (3)最大化:根据训练样本集f和f,以及隐变量的期望最大化来重新估计模型参 数
[0048] (4)按照步骤3.2和3.3进行反复迭代，直至满足收敛条件，所述收敛条件为：
[0050] 其中ε为收敛因子，Μ · | |2表示二范数。
[0051] 第四步:将建模数据和模型参数Θ存入历史数据库中备用；
[0052] 第五步：收集新的过程变量的在线数据：χ_=[χΓ，ΧΓ, ···，<'/,其中 X厂e 为当前t时刻的在线过程变量数据。对其标准化处理，得到X"?'；
[0053] 第六步:采用基于有监督的线性动态系统模型的软测量建模方法，根据历史数据 库中的数据和过程变量的在线数据，预测对应时刻难以直接测量的质量变量 Y-=tyr，yr，一，y：T，yr]，实现工业生产上对质量变量的监控。
[0054] 首先，根据有监督的线性动态系统模型参数Θ和过程变量的在线数据充-'，通过无监 督的线性动态系统模型的前向滤波方法计算出对应的隐变量的值[rvr，…，dr*]， 其中为当前t时刻隐变量的值，（这里是用隐变量的滤波均值来作为隐变量的 值），如下所示：
[0058]然后，根据有监督的线性动态系统模型参数Θ和历史数据库中的数据，预测出难以
[0060] 其中，σ!，σ2，…，oL为训练样本集γ中每个质量变量的方差，μι，μ 2，…，yL为训练样本 集Y中每个质量变量的均值
[0061] 以下结合一个具体的脱丁烷塔例子来说明本发明的有效性。脱丁烷塔是炼油厂脱 硫和石脑油分馏装置的重要组成部分。脱丁烷塔用于将石脑油中的丙烷和丁烷除去，质量 控制要求为:最小化脱丁烷塔底部丁烷的含量。该过程收集了 2394组数据，选取数据集中的 奇数样本作为训练集，偶数样本作为测试集。针对该过程，一共选取了 7个容易测量的变量 作为过程变量，如表1所示。质量变量为丁烷的含量。
[0062] 表1:监控变量说明
[0064] 接下来结合该具体过程对本发明的实施步骤进行详细地阐述：
[0065] 1.采集表1中的7个过程变量的数据，以及对其对应的丁烷含量值进行离线分析和 记 D
[0066]然后对所有数据进行预处理和标准化。
[0067]采集到的过程变量数据为XeR7xl197,对应的丁烷含量数据为YeRlxl197。剔除数据 中的野值点和粗糙误差点。然后标准化处理，即使得各个过程变量和质量变量的均值为〇， 方差为1，得到的新数据矩阵分别为X e i?M197和f e i?M197。
[0068] 2.针对训练数据，建立基于有监督的线性动态系统模型的软测量建模。
[0069] 根据训练数据1和?,按照实施步骤中给出的详细方法，建立基于有监督的线性动 态系统模型的软测量建模。
[0070] 3.获取在线过程变量的数据，并对其进行预处理和归一化。
[0071] 为了测试新方法的有效性，对在线测试样本集Xne?eR7xl197进行测试，并利用建模 时的标准化参数对其进行处理。
[0072] 4.在线预测丁烧含量
[0073] 根据有监督的线性动态系统模型的参数，对在线数据Xn~进行在线软测量，获得相 应的丁烷含量的在线预测值Y ne3WeRlxl197。图1给出了本发明方法针对在线测试数据的在线 预测结果。图2和图3分别给出了基于偏最小二乘模型和基于概率主元回归模型的在线预测 结果。从图中可以看出，考虑了过程动态性和数据随机性的本发明方法的在线预测效果要 好于偏最小二乘和概率主元回归方法，预测更加准确。
[0074] 上述实施例用来解释说明本发明，而不是对本发明进行限制，在本发明的精神和 权利要求的保护范围内，对本发明做出的任何修改和改变，都落入本发明的保护范围。
【主权项】
1. 一种基于有监督的线性动态系统模型的软测量建模方法，其特征在于，包括以下步 骤： (1) 利用集散控制系统以及离线检测方法收集工业生产中的过程变量和质量变量的数 据，组成建模用的训练样本集:X=[ X1，X2，."，XN]eRvx^Y=[yi， y2，."，yN]eRLXNj*，R 为实数集，RVXN表示X满足V X N的二维分布，RlxN表示Y满足L X N的二维分布，V为过程变量个 数，L为质量变量个数，N为采样数据点数，将数据存入历史数据库。 (2) 从历史数据库中调用训练样本集X和Y，分别对训练样本集X中的各个样本和Y中的 各个样本按照时间点方向进行排序，得到 y'2,…t，…n] eRaW jPy't分别为t时刻采集到的过程变量的训练样本和质量变 量的训练样本，t = 1，2，…，N。对每一个训练样本进行标准化处理，即使得各个过程变量和 质量变量的均值为〇，方差为1，得到的新数据矩阵分别为又=??，…A，…e fχΛ和 歹=民，5^-，5^-,^>，#。其中，写为匕经标准化处理后得到的样本，充为5^经标准 化处理后得到的样本，t = l，2,…，Ν。 (3) 根据训练样本集X和?采用期望最大化方法建立有监督的线性动态系统模型，得到 模型参数Θ。 (4) 将建模数据和模型参数Θ存入历史数据库中备用。 (5) 收集新的过程变量的在线数据〔中1￡厂eJTxl为当 前t时刻的在线过程变量数据。对其标准化处理，得到文 (6) 采用基于有监督的线性动态系统模型的软测量建模方法，根据历史数据库中的数 据和过程变量的在线数据文_，预测对应时刻难以直接测量的质量变量，实现工业生产上对质量变量的监控。2. 根据权利要求1所述基于有监督的线性动态系统模型的软测量建模方法，其特征在 于，所述步骤(3)具体为:对于训练样本集X和?，采用期望最大化方法求出有监督的线性 动态系统模型参数9 = {A，P，C，Ση，Σβ，Σ?，μπ，Σπ}，其中AER HXH为传递矩阵，Η为隐变量个 数;P e RVXH为映射矩阵，C e RLXH为回归矩阵；Σ n e RHXH为隐空间噪声的方差，Σ e e Rvxν为过 程变量噪声的方差，为质量变量噪声的方差;μπΕΚΗΧ4ΡΣ πΕΚΗΧΗ分别为服从高斯 分布的初始时刻隐变量heRa1的均值和方差。用期望最大化方法建模的具体实现步骤如 下所示： (3.1) 设置初始的模型参数0为0={六少，(：，2^2(3,2^11，211} (3.2) 求期望:在当前模型参数Θ下，根据训练样本集又和?计算每个时刻隐变量的平滑 均值gtERHX1、方差G tERHXH，t = l，2,…，N;相邻时刻隐变量的协方差文 V =t+l，且V =1，2,…，N。其中E( ·)表示括号中变量的均值，ht为t时刻Xt，yt对应的隐变 量。具体实现步骤如下所示： (3.2.1)通过前向滤波方法可以得到每个时刻隐变量的滤波均值f t e Rhx1和方差Ft e Rh XH，t=l，2，'"，N 如下：其中当时亥Ijt = 1时，Fj =:+ 0?}? + Σ:1)--1.,(3.2.2)通过后向平滑方法可以得到每个时刻隐变量的平滑均值gt e Rhx1，方差Gt e Rh 抑八=1，2，...札相邻时刻隐变量的协方差五^,4|文，交^俨'?4+1，且？=1，2，···， N。如下： gt = FtAT (AFtAT+ Σ η)-1 (gt' -Aft) +ft (3) Gt = FtAT(AFtAT+ Σn)-bt' [FtAT(AFtAT+ Σn)-1]T+Ft-FtA T(AFtAT+ Σn)-iFt (4) ,.,. (5) 其中当时亥lJt = T时，gT = fT，GT=FT。 (3.3) 最大化:根据训练样本集！和?，以及隐变量的期望最大化来重新估计模型参数(3.4) 按照步骤3.2和3.3进行反复迭代，直至满足收敛条件，所述收敛条件为：其中ε为收敛因子，| | · | |2表示二范数。3.根据权利要求1所述基于有监督的线性动态系统模型的软测量建模方法，其特征在 于，所述步骤(6)具体为： (6.1) 根据有监督的线性动态系统模型参数Θ和过程变量的在线数据又_，通过无监督 的线性动态系统模型的前向滤波方法计算出对应的隐变量的值[dv-cr'l，其 中e 为当前t时刻隐变量的值，（这里是用隐变量的滤波均值来作为隐变量的值）。 (6.2) 根据有监督的线性动态系统模型参数Θ和历史数据库中的数据，预测出难以直接 测量的质量变量￥_ = [〇Γ、.··ΟΓΜ?],如下所示：(15) 其中，σ:, σ2,…，为训练样本集Υ中每个质量变量的方差，μι, μ2,…，μΑ训练样本集Υ 中每个质量变量的均值。
【文档编号】G06F17/18GK105868164SQ201610165416
【公开日】2016年8月17日
【申请日】2016年3月19日
【发明人】葛志强, 陈新如
【申请人】浙江大学