CN102809966A - 基于半闭环的软测量仪表及其软测量方法 - Google Patents
基于半闭环的软测量仪表及其软测量方法 Download PDFInfo
- Publication number
- CN102809966A CN102809966A CN2012102662765A CN201210266276A CN102809966A CN 102809966 A CN102809966 A CN 102809966A CN 2012102662765 A CN2012102662765 A CN 2012102662765A CN 201210266276 A CN201210266276 A CN 201210266276A CN 102809966 A CN102809966 A CN 102809966A
- Authority
- CN
- China
- Prior art keywords
- variable
- module
- soft
- measurement
- auxiliary 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.)
- Granted
Links
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Testing Or Calibration Of Command Recording Devices (AREA)
Abstract
一种工业过程控制技术领域的基于半闭环的软测量仪表,包括用于测量辅助变量的智能仪表,存放数据的DCS数据库以及上位机,现场智能仪表与工业生产过程连接,上位机与DCS数据库及软测量显示仪连接,所述的上位机为软测量智能处理器,包括数据处理模块、测量模块、预估模块、补偿模块和存储模块;其中,数据处理模块对辅助变量和主导变量进行预处理,测量模块产生主导变量的初始测量值,预估模块预测下一时刻的辅助变量,而补偿模块利用预估模块的输出和实际辅助变量测量值的差对主导变量的初始测量值进行补偿。本发明实现了在线校正,有效地提高软测量的精度,具有数确定方便、使用范围广、软测量效果好、精度高的特点。
Description
技术领域
本发明涉及工业过程中软测量技术,具体是一种基于半闭环的软测量仪表及其软测量方法。
背景技术
随着生产技术的发展和生产过程的日益复杂,在工业生产过程中,大量关键性过程状态、产品质量等参数由于技术和经济上的原因是很难或无法进行在线测量的,这已成为制约生产安全、产品质量、产量及生产效益进一步提升的瓶颈。近年来,研究者们开始利用生产过程中所测量到的大量数据建立测量模型,以实现对未知变量的间接测量,这类技术被称为软测量技术。其基本思想是通过测量辅助变量(易测变量)来估计主导变量的实际值。
过程设备实际上处于动态的***环境下,然而传统的软测量模型是开环结构并且缺乏在线校正机制,因此当***受外界干扰产生突然的扰动,原有的软测量模型的测量效果无法得到持续的保持,需要进行实时校正。而现有的软测量技术无法进行在线校正,只能定期的读取新数据并使用新的离线分析数据进行模型的更新和调整。显然,这种办法实现起来费时费力。
针对上述情况,设计一种半闭环的软测量仪表及软测量技术,对开环软测量的预估值进行补偿,以实现在线校正,可以获得更好的软测量效果。
发明内容
本发明的目的在于克服上述现有技术的不足,提供一种基于半闭环的软测量仪表及软测量方法,在实现在线校正的同时提高软测量的准确性,以满足工业生产的要求。
为实现上述目的,本发明采取以下技术方案实现:
一种基于半闭环的软测量仪表,包括用于测量辅助变量的智能仪表、存放数据的DCS数据库、用于显示主导变量输出值的软测量显示仪和上位机,所述的智能仪表与工业生产过程连接,所述的上位机分别与DCS数据库及软测量显示仪连接,其特点在于,所述的上位机为软测量智能处理器,该软测量智能处理器包括:
数据处理模块,用于对辅助变量进行归一化处理和对主导变量进行反归一化处理;
测量模块,用于产生主导变量的初始测量值;
预估模块,用于预测下一时刻的辅助变量;
补偿模块,利用预估模块输出的预测值和测量模块的实际辅助变量测量值的差对主导变量的初始测量值进行补偿;和
存储模块,用于存储辅助变量的预测值。
一种利用上述基于半闭环的软测量仪表的软测量方法,其特点在于,该软测量方法包括如下步骤:
步骤(1)设置测量辅助变量的智能仪表,将智能仪表测量到的数据传送到DCS实时数据库;
步骤(2)初始化k=1,在当前时刻k从DCS实时数据库中得到最新的变量数据x(1)作为辅助变量并送入上位机的数据处理模块,数据处理模块对辅助变量x(1)进行归一化处理,公式如下:
x*(k)=β+(x(k)-a)/(b-a)×(α-β) (1)
其中,x(k)表示辅助变量的实际值,x*(k)表示辅助变量归一化后的实际值,[a b]表示辅助变量的变化范围,[αβ]表示归一化范围,将归一化后的辅助变量x*(1)送入测量模块,得到主导变量的预测值y*'(1),同时当前时刻主导变量的软测量值y*(1)=y*'(1);
步骤(3)将当前时刻主导变量的软测量y*(1)送入数据处理模块,经过公式(2)反归一化后,得到主导变量的实际输出值y(1),
y(k)=(y*(k)-β)/(α-β)×(b-a)+a (2)
其中,y(k)表示主导变量的实际输出值,y*(k)表示主导变量的软测量值;
步骤(4)将归一化后的辅助变量x*(1)和主导变量的软测量值y*(1)送入预估模块得到下一时刻(k=2)辅助变量的预测值x*'(2)并在存储模块中保存,并置k=k+1;
步骤(5)在当前时刻k,从DCS实时数据库中得到最新的变量数据x(k)的在线实测数据送入上位机的数据处理模块,对其按(1)式进行归一化处理,将处理后的辅助变量x*(k)送入测量模块,得到主导变量的预测值y*'(k);
步骤(6)调用存储模块中辅助变量的预测值x*'(k),预测值x*'(k)与辅助变量归一化后的实际值x*(k)之间形成差值Δx*(k)=x*(k)-x*'(k),将Δx*(k)送入补偿模块;图1中+号表示Δx*(k)=x*(k)-x*'(k)中x*(k)的符号,-号为x*'(k)的符号;
步骤(7)通过补偿模块得到k时刻主导变量的补偿值Δy*(k),用k时刻主导变量的补偿值Δy*(k)补偿主导变量的预测值y*'(k)得到主导变量的软测量值y*(k)=y*'(k)+Δy*(k),
步骤(8)将主导变量的软测量值y*(k)送入数据处理模块,经过公式(2)反归一化后得到主导变量的实际输出值y(k)。
步骤(9)将主导变量的实际输出值传给软测量显示仪进行显示。
步骤(10)将当前时刻k归一化后的辅助变量x*(k)和主导变量的软测量值y*(k)分别送入预估模块预测下一时刻(k+1时刻)的辅助变量的预测值x*'(k+1)并在存储模块中保存同时覆盖上一时刻辅助变量的预测值,此时置k=k+1,返回步骤(5)继续进行后续的软测量。
本发明中的测量模块、预估模块、补偿模块的建立按照以下步骤完成:
步骤(1)从DCS历史数据库以及离线化验数据库中分析搜索相关的易测辅助变量与难测主导变量y(k)构成若干组数据样本用于半闭环软测量技术的训练,对数据进行归一化处理后,按采样时间递增的顺序形成下述相应的时间序列:
k=1,2…M j=1,2…N
步骤(2)建立测量模块:选取若干组样本,采用常用的映射技术(如,神经网络NN,支持向量机SVM等)并利用相关算法(BP算法,遗传算法GA)(参考《软测量技术原理与应用》,潘立登等编,中国电力出版社)训练得到映射组成测量模型,并送入到计算机中的测量模块中储存。
与现有技术相比,本发明利用既有数据构造半闭环框架,对软测量结果进行了校正和补偿,具有以下优点:
(1)本发明通过半闭环结构实现了在线校正,可以有效提高软测量的精确度。
(2)本发明提供的是一种框架结构,在3个子模块中,具体采用哪种映射技术和训练算法不受限制,因此具有很强的扩展性。
附图说明
图1是半闭环软测量技术上位机功能结构图;
图2是三层前向神经网络框图;
图3是开环软测量和半闭环软测量模型的仿真结果。
具体实施方式
下面结合附图对本发明的实施例作详细说明:本实施例在以本发明技术方案为前提下进行实施,给出了详细的实施方式和具体的操作过程,但本发明的保护范围不限于下述的实施例。
本实施例是对循环流化床锅炉中的烟气含氧量(FGOC)进行软测量。
由于三层以上结构的反向传播(BP,Back propagation)神经网络具有映射任意非线性函数的功能,这里,采用三层结构的BP神经网络来建立各子***的软测量模块,如图2所示;将得到的子模块按图1的框架组合,就得到FGOC的半闭环软测量框架。
FGOC的软测量技术包括以下几个步骤:
步骤(1)辅助变量的选取和数据的预处理
通过分析***流程,选取给煤量、一次风、二次风作为辅助变量来预测FGOC。在现场的DCS中,保存有大量的过程历史数据,从历史数据库中选取一天的数据用以建模,每十分钟计算一次平均值作为一个采样时刻的值。在经过滤波处理(参考《信号处理——模拟与数字信号***及滤波器》,迈耶尔编,机械工业出版社)后,选取500组数据作为训练样本,100组数据作为测试样本。为了消除不同量纲数据对模型的影响,需要利用式(1)对所有数据进行归一化处理。
步骤(2):构造FGOC的测量模块
FGOC的测量模块是三输入——单输出的,选取神经网络的隐层节点数nh=6。使用PSO算法对模型进行训练,PSO算法的迭代公式如下:
vid=ωvid+c1r1(pid-xid)+c2r2(pgd-xid) (3)
xid=xid+vid (4)
其中,每一个粒子代表神经网络的一个连接权值,因此将权值编码为粒子,粒子在D维搜索空间中位置表示为:Xi=(xi1,xi2,..,xiD),第i个粒子经历过的最好位置记为Pi=(pi1,pi2,...,piD),每个粒子的飞行速度为Vi=(vi1,vi2,...viD),i=1,2,...,m,在整个群体中,所有粒子经历过的最好位置为Pg=(pi1,pi2,...piD),每一代粒子根据式(3)、(4)来更新自己的速度和位置。PSO的参数设置:种群数目N为30,学习因子C1=C2=1.49,权重因子ω=0.73。
适应值函数选为:
训练误差err由下式给出:
测试误差err由下式给出:
最大适应值计算次数为2000,算法运行50次,选择最优适应值对应的神经网络作为映射关系,建立FGOC的测量模块,并送入到计算机中的测量模块中储存。
步骤(3):构造FGOC的预估模块
FGOC的预估模块是四输入——三输出的,选取神经网络的隐层节点数nh=8。适应值函数的选取和算法参数的设置同步骤(2)。使用PSO算法对模型进行训练,最大适应值计算次数为2000,算法运行50次,选择最优适应值对应的神经网络作为映射关系,建立FGOC的预估模块,并送入到计算机中的预估模块中储存。
步骤(4):构造FGOC的补偿模块
FGOC的补偿模块是三输入——单输出的,选取神经网络的隐层节点数nh=6。适应值函数的选取和算法参数的设置同步骤(2)。使用PSO算法对模型进行训练,最大适应值计算次数为2000,算法运行50次,选择最优适应值对应的神经网络作为映射关系,建立FGOC的补偿模块,并送入到计算机中的补偿模块中储存。
步骤(5):通过测量、预估、补偿3个模块来构造FGOC的半闭环软测量技术:
1)设置给煤量、一次风和二次风的数据采集设备,将采集到的数据传送到DCS实时数据库。在每个定时周期k从DCS的实时数据库中得到最新的变量数据x(k)作为辅助变量并送入上位机,所述上位机内部预设有数据处理模块、测量模块、预估模块、存储模块和补偿模块。
2)初始化:在k=1时,辅助变量x(1)的在线实测数据送入上位机的数据处理模块,对(1)式进行归一化处理。将处理后的辅助变量x*(1)送入测量模块,得到FGOC的预测值y*'(1),同时FGOC的软测量值y*(1)=y*'(1)。经过(2)式反归一化后得到FGOC的实际输出值y(1)。将x*(1)和y*(1)送入预估模块得到下一时刻(k=2)辅助变量的预测值x*'(2)并在存储模块中保存,置k=k+1。
3)当k≥2时,将辅助变量x(k)的在线实测数据送入上位机的数据处理模块,对其按(1)式进行归一化处理。将处理后的辅助变量x*(k)送入测量模块,得到FGOC在k时刻的预测值y*'(k)。
4)调用存储模块中辅助变量的预测值x*'(k),它与辅助变量归一化后的实际值x*(k)之间形成差值Δx*(k)=x*(k)-x*'(k),将Δx*(k)送入补偿模块。
5)Δx*(k)进入补偿模块,得到k时刻FGOC的补偿值Δy*(k),用Δy*(k)补偿预测值y*'(k)得到FGOC软测量值y*(k)=y*'(k)+Δy*(k)。将y*(k)送入数据处理模块,经过(2)式反归一化后得到FGOC的实际输出值y(k)。
6)将FGOC的实际输出值传给软测量显示仪进行显示。
7)将x*(k)和y*(k)送入计算机内部预设的预估模块预测下一时刻(k+1时刻)的辅助变量的预测值x*'(k+1)并在存储模块中保存同时覆盖上一时刻辅助变量的预测值。此时置k=k+1,返回步骤3)继续进行后续的软测量。
图3给出开环软测量和半闭环软测量模型的仿真结果,表1为两种方法的比较,可见,采用本发明进行在线软测量,其测量精度于鲁棒性与开环软测量相比有较大提高。
表1开环软测量和半闭环软测量模型的比较
Claims (6)
1.一种基于半闭环的软测量仪表,包括用于测量辅助变量的智能仪表、存放数据的DCS数据库、用于显示主导变量输出值的软测量显示仪和上位机,所述的智能仪表与工业生产过程连接,所述的上位机分别与DCS数据库及软测量显示仪连接,其特征在于,所述的上位机为软测量智能处理器,该软测量智能处理器包括:
数据处理模块,用于对辅助变量进行归一化处理和对主导变量进行反归一化处理;
测量模块,用于产生主导变量的初始测量值;
预估模块,用于预测下一时刻的辅助变量;
补偿模块,利用预估模块输出的预测值和测量模块的实际辅助变量测量值的差对主导变量的初始测量值进行补偿;和
存储模块,用于存储辅助变量的预测值。
6.一种利用权利要求1-5任一项所述的基于半闭环的软测量仪表的软测量方法,其特征在于,该软测量方法包括如下步骤:
步骤(1)设置测量辅助变量的智能仪表,将智能仪表测量到的数据传送到DCS实时数据库;
步骤(2)初始化k=1,在当前时刻k从DCS实时数据库中得到最新的变量数据x(1)作为辅助变量并送入上位机的数据处理模块,数据处理模块对辅助变量x(1)进行归一化处理,公式如下:
x*(k)=β+(x(k)-a)/(b-a)×(α-β) (1)
其中,x(k)表示辅助变量的实际值,x*(k)表示辅助变量归一化后的实际值,[a b]表示辅助变量的变化范围,[αβ]表示归一化范围,将归一化后的辅助变量x*(1)送入测量模块,得到主导变量的预测值y*'(1),同时当前时刻主导变量的软测量
y*(1)=y*'(1);
步骤(3)将当前时刻主导变量的软测量y*(1)送入数据处理模块,经过公式(2)反归一化后,得到主导变量的实际输出值y(1),
y(k)=(y*(k)-β)/(α-β)×(b-a)+a (2)
其中,y(k)表示主导变量的实际值,y*(k)表示主导变量的软测量值;
步骤(4)将归一化后的辅助变量x*(1)和主导变量的软测量值y*(1)送入预估模块得到下一时刻(k=2)辅助变量的预测值x*'(2)并在存储模块中保存,并置k=k+1;
步骤(5)在当前时刻k,从DCS实时数据库中得到最新的变量数据x(k)的在线实测数据送入上位机的数据处理模块,对其按(1)式进行归一化处理,将处理后的辅助变量x*(k)送入测量模块,得到主导变量的预测值y*'(k);
步骤(6)调用存储模块中辅助变量的预测值x*'(k),预测值x*'(k)与辅助变量归一化后的实际值x*(k)之间形成差值Δx*(k)=x*(k)-x*'(k),将Δx*(k)送入补偿模块;步骤(7)通过补偿模块得到k时刻主导变量的补偿值Δy*(k),用k时刻主导变量的补偿值Δy*(k)补偿主导变量的预测值y*'(k)得到主导变量的软测量值y*(k)=y*′(k)+Δy*(k),
步骤(8)将主导变量的软测量值y*(k)送入数据处理模块,经过公式(2)反归一化后得到主导变量的实际输出值y(k);
步骤(9)将主导变量的实际输出值传给软测量显示仪进行显示;
步骤(10)将当前时刻k归一化后的辅助变量x*(k)和主导变量的软测量值y*(k)分别送入预估模块预测下一时刻(k+1时刻)的辅助变量的预测值x*'(k+1)并在存储模块中保存同时覆盖上一时刻辅助变量的预测值,此时置k=k+1,返回步骤(5)继续进行后续的软测量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210266276.5A CN102809966B (zh) | 2012-07-30 | 2012-07-30 | 基于半闭环的软测量仪表及其软测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210266276.5A CN102809966B (zh) | 2012-07-30 | 2012-07-30 | 基于半闭环的软测量仪表及其软测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102809966A true CN102809966A (zh) | 2012-12-05 |
CN102809966B CN102809966B (zh) | 2014-10-15 |
Family
ID=47233680
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210266276.5A Active CN102809966B (zh) | 2012-07-30 | 2012-07-30 | 基于半闭环的软测量仪表及其软测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102809966B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103064286A (zh) * | 2013-01-17 | 2013-04-24 | 丽水学院 | 一种工业过程的控制方法与设备 |
WO2021093362A1 (zh) * | 2019-11-11 | 2021-05-20 | 东北大学 | 一种蒸发过程出料苛性碱浓度测量装置精度补偿方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH10244286A (ja) * | 1997-03-04 | 1998-09-14 | Mitsubishi Chem Corp | 活性汚泥法排水処理設備の栄養源供給量制御装置 |
US6639402B2 (en) * | 2001-01-31 | 2003-10-28 | University Of Kentucky Research Foundation | Temperature, stress, and corrosive sensing apparatus utilizing harmonic response of magnetically soft sensor element (s) |
CN101630376A (zh) * | 2009-08-12 | 2010-01-20 | 江苏大学 | 多模型神经网络的生物发酵过程软测量建模方法及软仪表 |
-
2012
- 2012-07-30 CN CN201210266276.5A patent/CN102809966B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH10244286A (ja) * | 1997-03-04 | 1998-09-14 | Mitsubishi Chem Corp | 活性汚泥法排水処理設備の栄養源供給量制御装置 |
US6639402B2 (en) * | 2001-01-31 | 2003-10-28 | University Of Kentucky Research Foundation | Temperature, stress, and corrosive sensing apparatus utilizing harmonic response of magnetically soft sensor element (s) |
CN101630376A (zh) * | 2009-08-12 | 2010-01-20 | 江苏大学 | 多模型神经网络的生物发酵过程软测量建模方法及软仪表 |
Non-Patent Citations (2)
Title |
---|
刘兴业等: "基于多个软测量模型集成的温度控制", 《计算机与数字工程》, vol. 39, no. 8, 31 August 2011 (2011-08-31), pages 169 - 2 * |
彭小奇等: "一种软测量模型性能监测评价及其自适应校正方法", 《化工学报》, vol. 63, no. 5, 31 May 2012 (2012-05-31) * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103064286A (zh) * | 2013-01-17 | 2013-04-24 | 丽水学院 | 一种工业过程的控制方法与设备 |
WO2021093362A1 (zh) * | 2019-11-11 | 2021-05-20 | 东北大学 | 一种蒸发过程出料苛性碱浓度测量装置精度补偿方法 |
US11893472B2 (en) | 2019-11-11 | 2024-02-06 | Northeastern University | Accuracy compensation method for discharge caustic alkali concentration measuring device in evaporation process |
Also Published As
Publication number | Publication date |
---|---|
CN102809966B (zh) | 2014-10-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102339040B (zh) | 精馏塔优化控制方法 | |
CN103205665B (zh) | 一种连续热镀锌线锌层厚度自动控制方法 | |
CN109669009A (zh) | 基于域校正的域自适应极限学习机电子鼻漂移抑制方法 | |
CN102411305A (zh) | 单旋翼直升机/涡轴发动机综合抗扰控制***设计方法 | |
CN102778538A (zh) | 一种基于改进支持向量机的锅炉飞灰含碳量软测量方法 | |
CN106684857A (zh) | 一种含统一潮流控制器的线性化最优潮流模型 | |
CN107966242B (zh) | 适用于深海浮标的一键式气压传感器现场校准***及方法 | |
CN110503268B (zh) | 一种模型数据联合驱动的综合能源***态势感知方法 | |
CN107704962A (zh) | 一种基于不完整时间序列数据集的冶金企业蒸汽流量区间预测方法 | |
CN109508818A (zh) | 一种基于LSSVM的在线NOx预测方法 | |
CN102880905A (zh) | 一种常顶油干点在线软测量方法 | |
CN103245759A (zh) | 基于自回归全潜结构投影模型的产品质量监测方法 | |
CN101879378B (zh) | 内部热耦合精馏塔温度非线性观测***及方法 | |
CN104063577A (zh) | 基于广义回归神经网络的变压器油中特征气体发展趋势预测方法 | |
CN106779139A (zh) | 基于小波分解和二阶灰色神经网络的短期风速预测方法 | |
CN107085371A (zh) | 基于数据驱动的原油蒸馏装置经济模型预测控制方法 | |
CN102809966B (zh) | 基于半闭环的软测量仪表及其软测量方法 | |
CN104865944B (zh) | 基于pca‑lssvm的气分装置控制***性能评估方法 | |
CN101285816A (zh) | 铜锍吹炼的过程参数软测量仪表及其软测量方法 | |
CN103344740B (zh) | 基于多输入维纳模型的谷氨酸产物浓度在线软测量方法 | |
CN103279030A (zh) | 基于贝叶斯框架的动态软测量建模方法及装置 | |
Si et al. | A new approach for function approximation in boiler combustion optimization based on modified structural AOSVR | |
CN108204997A (zh) | 常一线油闪点在线软测量方法 | |
CN109034497A (zh) | 多晶硅还原工序能耗值的预测方法、***、介质及设备 | |
CN102043380A (zh) | 基于二次多项式的非线性复合pid神经网络控制方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |