CN101890246A - 一种精馏塔温度非线性观测系统及方法 - Google Patents
一种精馏塔温度非线性观测系统及方法 Download PDFInfo
- Publication number
- CN101890246A CN101890246A CN 201010213753 CN201010213753A CN101890246A CN 101890246 A CN101890246 A CN 101890246A CN 201010213753 CN201010213753 CN 201010213753 CN 201010213753 A CN201010213753 A CN 201010213753A CN 101890246 A CN101890246 A CN 101890246A
- Authority
- CN
- China
- Prior art keywords
- temperature
- sampling instant
- max
- alpha
- rectifying
- 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
- Feedback Control In General (AREA)
Abstract
一种精馏塔温度非线性观测系统,包括与精馏塔直接连接的现场智能仪表和DCS系统,所述DCS系统包括存储装置、控制站和上位机,所述现场智能仪表与存储装置、控制站和上位机相连,所述的上位机包括用以预测精馏塔将来时刻的塔板温度分布的观测器,所述的观测器包括温度分布函数拟合模块、将来时刻精馏塔塔板温度预测模块和温度预测误差补偿模块。本发明也提供了一种精馏塔的非线性温度观测方法。本发明能够有效快速的在线预测精馏塔将来时刻的温度变化,准确地反映精馏塔的强非线性动态特性,实现高效率、高精度的温度观测。
Description
技术领域
本发明涉及精馏塔的动态非线性建模领域,特别提出精馏塔温度非线性观测器系统的设计和实现方法。
背景技术
精馏过程是化工过程中的一种核心过程,精馏塔是其中的一个关键单元。长久以来,精馏塔因为高耗能,低能效的问题成为国际精馏领域研究的焦点。精馏塔具有十分复杂的强非线性特征。导致精馏塔的简化非线性建模,非线性动态过程观测器设计成为一种世界性难题。
传统的复杂机理模型尽管能够精确的观测精馏塔的各种状态变量,但复杂的机理模型结构导致模型在线求解效率太低,不能直接应用于控制器设计当中。而基于阶跃、脉冲响应辨识的近似线性模型尽管具有较快的求解速度但是状态变量观测精度太低,甚至显著偏离真实值。可见实现精馏塔的高效率高精确度的非线性观测器对于该塔的控制设计节能研究等方面起着至关重要的作用。
发明内容
为了克服现有的精馏塔温度观测系统的状态变量观测精度太低、效率低的不足,本发明提供一种精馏塔的高效率、高精度的温度观测的精馏塔温度非线性观测系统及方法、
本发明解决其技术问题所采用的技术方案是:
一种精馏塔温度非线性观测系统,包括与精馏塔直接连接的现场智能仪表和DCS系统,所述DCS系统包括存储装置、控制站和上位机,所述现场智能仪表与存储装置、控制站和上位机相连,所述的上位机包括用以预测精馏塔将来时刻的温度分布的观测器,所述观测器包括:
温度分布函数拟合模块,用以采用历史数据库中的各块塔板温度Ti拟合温度分布函数,并将拟合参数存储到历史数据库当中,采用式(1)(2)得到:
其中,1为塔顶编号,f为进料板编号,n为塔底编号,为第i块塔板处预估温度,Tmin,r、Tmax,r、kr、Tmin,s、Tmax,s、ks为拟合参数,Sr、Ss分别为精馏塔精馏段、提馏段液温度分布曲线的位置;
将来时刻精馏塔温度预测模块,用以根据当前温度数据以及温度分布函数预测将来时刻的温度变化趋势,该模块包括以下两部分:
①从历史数据库中读取k采样时刻第i块塔板温度测量值Ti(k),计算当前时刻温度分布曲线位置变化速度,并将曲线位置变化速度存储到历史数据库当中,采用式(3)-(6)得到:
其中,k为当前采样时刻,P为精馏塔压强,α为相对挥发度,a、b、c为安东尼常数,Ti(k)、Xi(k)、Yi(k)分别为k采样时刻第i块塔板温度测量值、液相轻组分浓度和汽相轻组分浓度,F为进料流率,V、L分别为气、液相摩尔流率,H为持液量,Zf为进料组分浓度,Xf-1(k)、Xn(k)分别为k采样时刻第f-1块塔板和塔底的液相轻组分浓度,Y1(k)、Yf(k)分别为k采样时刻塔顶和进料板处的汽相轻组分浓度, 分别为k采样时刻精馏段和提馏段的温度分布曲线位置变化速度;
②从历史数据库中获取温度分布曲线位置变化速度数据,温度分布函数拟合参数,采样周期并计算下一时刻波形位置以及下一时刻温度分布,采用式(7)-(10)得到:
其中,δ为采样周期,Sr(k)、Ss(k)分别为k采样时刻精馏段、提馏段温度分布曲线位置,Sr(k+1)、Ss(k+1)分别为k+1采样时刻精馏段、提馏段温度分布曲线位置,为k+1时刻精馏段或提馏段的第i块塔板温度的预测值;
温度预测误差补偿模块,用以从历史数据库中获取当前温度测量值,以及预测值,计算预测误差,修正温度预测值,采用式(11)(12)得到
作为优选的一种方案:所述的上位机还包括人机界面模块,用以设定采样周期,显示观测器输出的将来时刻温度预测值和当前时刻温度测量值。
一种精馏塔的温度非线性过程观测方法,所述的观测方法包括以下步骤:
1)确定采样周期,并将相对挥发度α,安东尼常数a、b、c、保存在历史数据库当中;
2)检测k采样时刻精馏塔压强P,以及各塔板温度Ti,计算液相轻组分浓度值:
其中,下脚标i为塔板编号,1为塔顶编号,f为进料板编号,n为塔底编号,Xi为第i块塔板的液相组分浓度,P为精馏塔压强,Ti为第i块塔板的温度,α为相对挥发度,a、b、c为安东尼常数;
3)从历史数据库当中获取k采样时刻温度的测量值,拟合k采样时刻的温度分布函数如式(1)(2):
其中,为第i块塔板处预估温度,Tmin,r、Tmax,r、kr、Tmin,s、Tmax,s、ks为拟合参数,Sr、Ss分别为精馏塔精馏段、提馏段液温度分布曲线的位置;
4)计算k采样时刻的温度分布曲线位置变化速度,并将曲线位置变化速度存储到历史数据库当中,采用式(4)-(6)得到:
其中,k为当前采样时刻,P为精馏塔压强,α为相对挥发度,a、b、c为安东尼常数,Ti(k)、Xi(k)、Yi(k)分别为k采样时刻第i块塔板温度测量值、液相轻组分浓度和汽相轻组分浓度,F为进料流率,V、L分别为气、液相摩尔流率,H为持液量,Zf为进料组分浓度,Xf-1(k),Xn(k)分别为k采样时刻第f-1块塔板和塔底的液相轻组分浓度,Y1(k)、Yf(k)分别为k采样时刻塔顶和进料板处的汽相轻组分浓度,分别为k采样时刻精馏段和提馏段的温度分布曲线位置变化速度;
5)预测k+1采样时刻的温度分布,从历史数据库中获取曲线位置变化速度数据,温度分布函数拟合参数,采样周期值,采用式(7)-(10)得到:
其中,δ为采样周期,Sr(k)、Ss(k)分别为k采样时刻精馏段、提馏段温度分布曲线位置,Sr(k+1)、Ss(k+1)分别为k+1采样时刻精馏段、提馏段温度分布曲线位置,为k+1时刻精馏段、提馏段段第i块塔板温度的预测值;
6)从历史数据库中获取当前温度测量值,以及预测值,计算k采样时刻的温度预测误差过程如式(11):
7)修正温度预测值,采用式(12)得到:
进一步,所述的历史数据库为DCS系统的存储装置,控制站读取历史数据库,显示观测器工作过程状态。
本发明的有益效果主要表现在:1.基于精馏塔的机理模型以及传热传质原理建立了高精度的精馏塔的简化非线性动态模型,能够准确跟踪反映精馏塔温度分布的各种非线性特征,状态预测比传统辨识模型更为精确;2.在确保观测器跟踪高精度效能的情况下,能够同时保持较快的在线求解速度,效率远远高于机理模型
附图说明
图1是本发明所提出的精馏塔温度非线性观测系统的结构图。
图2是上位机观测器实现方法的原理图。
图3是本发明所提出的精馏塔非线性观测方法流程图。
具体实施方式
下面根据附图具体说明本发明。
实施例1
参照图1,图2和图3,一种精馏塔温度非线性观测系统,包括与精馏塔1直接连接的现场智能仪表2和DCS系统13,所述DCS系统13包括存储装置4、控制站5及上位机6,所述现场智能仪表2与存储装置4、控制站5和上位机6通过现场总线和数据接口3依次相连;所述的智能仪表用以实现检测模块7的功能即检测得到精馏塔的温度、压强数据,所述的数据接口用以实现I/O模块8的功能即数据的输入输出。
所述的上位机6包括用以预测精馏塔将来时刻的温度分布的观测器,所述的观测器包括温度分布函数拟合模块9,将来时刻精馏塔温度预测模块11,温度预测误差补偿模块10。
所述的温度分布函数拟合模块9,采用历史数据库中的各块塔板温度Ti拟合温度分布函数,并将拟合参数存储到历史数据库当中,采用式(1)(2)得到:
其中,1为塔顶编号,f为进料板编号,n为塔底编号,为第i块塔板处预估温度,Tmin,r、Tmax,r、kr、Tmin,s、Tmax,s、ks为拟合参数,Sr、Ss分别为精馏塔精馏段、提馏段液温度分布曲线的位置。
所述的将来时刻精馏塔温度预测模块11,根据当前温度数据以及温度分布函数预测将来时刻的温度变化趋势,该模块包括以下两部分:
①从历史数据库中读取温度测量值Ti(k),计算当前时刻温度分布曲线位置变化速度,并将曲线位置变化速度存储到历史数据库当中,采用式(3)-(6)得到:
其中,k为当前采样时刻,P为精馏塔压强,α为相对挥发度,a、b、c为安东尼(Antonie)常数,Ti(k)、Xi(k)、Yi(k)分别为k采样时刻第i块塔板温度测量值、液相轻组分浓度和汽相轻组分浓度,F为进料流率,V、L分别为气、液相摩尔流率,H为持液量,Zf为进料组分浓度,Xf-1(k),Xn(k)分别为k采样时刻第f-1块塔板和塔底的液相轻组分浓度,Y1(k)、Yf(k)分别为k采样时刻塔顶和进料板处的汽相轻组分浓度,分别为k采样时刻精馏段和提馏段的温度分布曲线位置变化速度。
②从历史数据库中获取温度分布曲线位置变化速度数据,温度分布函数拟合参数,采样周期并计算下一时刻波形位置以及下一时刻温度分布,采用式(7)-(10)得到:
其中,δ为采样周期,Sr(k)、Ss(k)分别为k采样时刻精馏段、提馏段温度分布曲线位置,Sr(k+1)、Ss(k+1)分别为k+1采样时刻精馏段、提馏段温度分布曲线位置,为k+1时刻精馏段或提馏段的第i块塔板温度的预测值。
所述的温度预测误差补偿模块10,从历史数据库中获取当前温度测量值,以及预测值,计算预测误差,修正温度预测值,采用式(11)(12)得到
所述的上位机包括人机界面模块12,用于设定采样周期,显示观测器输出的将来时刻温度预测值和当前时刻温度测量值。
实施例2
参照图1、图2和图3,一种精馏塔的温度非线性过程观测方法,所述的观测方法包括以下步骤:
1)如图3模块14确定采样周期,并将相对挥发度α,安东尼(Antonie)常数a、b、c、保存在历史数据库当中;
2)如图3模块15检测k采样时刻精馏塔压强P,以及各塔板温度Ti,计算液相轻组分浓度值:
其中,下脚标i为塔板编号,1为塔顶编号,f为进料板编号,n为塔底编号,Xi为第i块塔板的液相组分浓度,P为精馏塔压强,Ti为第i块塔板的温度,α为相对挥发度,a、b、c为安东尼(Antonie)常数;
(3)如图3模块16从历史数据库当中获取k采样时刻温度的测量值,拟合k采样时刻的温度分布函数如式(1)(2):
4)如图3模块17计算k采样时刻的温度分布曲线位置变化速度,并将曲线位置变化速度存储到历史数据库当中,采用式(4)-(6)得到
其中,k为当前采样时刻,P为精馏塔压强,α为相对挥发度,a、b、c为安东尼常数,Ti(k)、Xi(k)、Yi(k)分别为k采样时刻第i块塔板温度测量值、液相轻组分浓度和汽相轻组分浓度,F为进料流率,V、L分别为气、液相摩尔流率,H为持液量,Zf为进料组分浓度,Xf-1(k),Xn(k)分别为k采样时刻第f-1块塔板和塔底的液相轻组分浓度,Y1(k)、Yf(k)分别为k采样时刻塔顶和进料板处的汽相轻组分浓度,分别为k采样时刻精馏段和提馏段的温度分布曲线位置变化速度;
5)如图3模块18预测k+1采样时刻的温度分布,从历史数据库中获取曲线位置变化速度数据,温度分布函数拟合参数,采样周期值,采用式(7)-(10)得到:
其中,δ为采样周期,Sr(k)、Ss(k)分别为k采样时刻精馏段、提馏段温度分布曲线位置,Sr(k+1)、Ss(k+1)分别为k+1采样时刻精馏段、提馏段温度分布曲线位置,为k+1时刻精馏段或提馏段的第i块塔板温度的预测值;
6)如图3模块19从历史数据库中获取当前温度测量值,以及预测值,计算k采样时刻的温度预测误差过程如式(11):
其中,ri(k)为k采样时刻第i块塔板的温度预测误差;
7)如图3模块20修正温度预测值,采用式(12)得到:
所述的历史数据库为DCS系统中上午存储装置4,所述的DCS系统13包括数据接口3,存储装置4,控制站5,其中控制站5可以读取历史数据库,显示观测器工作过程状态。
上述实施例用来解释说明本发明,而不是对本发明进行限制,在本发明的精神和权利要求的保护范围内,对本发明作出的任何修改和改变,都落入本发明的保护范围。
Claims (4)
1.一种精馏塔温度非线性观测系统,包括与精馏塔直接连接的现场智能仪表和DCS系统,所述DCS系统包括存储装置、控制站和上位机,所述现场智能仪表与存储装置、控制站和上位机相连,其特征在于:所述的上位机包括用以预测精馏塔将来时刻的温度分布的观测器,所述观测器包括:
温度分布函数拟合模块,用以采用历史数据库中的各块塔板温度Ti拟合温度分布函数,并将拟合参数存储到历史数据库当中,采用式(1)(2)得到:
其中,1为塔顶编号,f为进料板编号,n为塔底编号,为第i块塔板处预估温度,Tmin,r、Tmax,r、kr、Tmin,s、Tmax,s、ks为拟合参数,Sr、Ss分别为精馏塔精馏段、提馏段液温度分布曲线的位置;
将来时刻精馏塔温度预测模块,用以根据当前温度数据以及温度分布函数预测将来时刻的温度变化趋势,该模块包括以下两部分:
①从历史数据库中读取k采样时刻第i块塔板温度测量值Ti(k),计算当前时刻温度分布曲线位置变化速度,并将曲线位置变化速度存储到历史数据库当中,采用式(3)-(6)得到:
其中,k为当前采样时刻,P为精馏塔压强,α为相对挥发度,a、b、c为安东尼常数,Ti(k)、Xi(k)、Yi(k)分别为k采样时刻第i块塔板温度测量值、液相轻组分浓度和汽相轻组分浓度,F为进料流率,V、L分别为气、液相摩尔流率,H为持液量,Zf为进料组分浓度,Xf-1(k)、Xn(k)分别为k采样时刻第f-1块塔板和塔底的液相轻组分浓度,Y1(k)、Yf(k)分别为k采样时刻塔顶和进料板处的汽相轻组分浓度, 分别为k采样时刻精馏段和提馏段的温度分布曲线位置变化速度;
②从历史数据库中获取温度分布曲线位置变化速度数据,温度分布函数拟合参数,采样周期并计算下一时刻波形位置以及下一时刻温度分布,采用式(7)-(10)得到:
其中,δ为采样周期,Sr(k)、Ss(k)分别为k采样时刻精馏段、提馏段温度分布曲线位置,Sr(k+1)、Ss(k+1)分别为k+1采样时刻精馏段、提馏段温度分布曲线位置,为k+1采样时刻精馏段或提馏段的第i块塔板温度的预测值;
温度预测误差补偿模块,用以从历史数据库中获取当前温度测量值,以及预测值,计算预测误差,修正温度预测值,采用式(11)(12)得到
2.如权利要求1所述的精馏塔温度非线性观测器系统,其特征在于:所述的上位机还包括人机界面模块,用以设定采样周期,显示观测器输出的将来时刻温度预测值和当前时刻温度测量值。
3.一种如权利要求1所述的精馏塔的温度非线性过程观测系统实现的温度非线性过程观测方法,其特征在于:所述的观测方法包括以下步骤:
1)确定采样周期,并将相对挥发度α,安东尼常数a、b、c、保存在历史数据库当中;
2)检测k采样时刻精馏塔压强P,以及各塔板温度Ti,计算液相轻组分浓度值:
其中,下脚标i为塔板编号,1为塔顶编号,f为进料板编号,n为塔底编号,Xi为第i块塔板的液相组分浓度,P为精馏塔压强,Ti为第i块塔板的温度,α为相对挥发度,a、b、c为安东尼常数;
3)从历史数据库当中获取k采样时刻温度的测量值,拟合k采样时刻的温度分布函数如式(1)(2):
4)计算k采样时刻的温度分布曲线位置变化速度,并将曲线位置变化速度存储到历史数据库当中,采用式(4)-(6)得到:
其中,k为当前采样时刻,P为精馏塔压强,α为相对挥发度,a、b、c为安东尼常数,Ti(k)、Xi(k)、Yi(k)分别为k采样时刻第i块塔板温度测量值、液相轻组分浓度和汽相轻组分浓度,F为进料流率,V、L分别为气、液相摩尔流率,H为持液量,Zf为进料组分浓度,Xf-1(k),Xn(k)分别为k采样时刻第f-1块塔板和塔底的液相轻组分浓度,Y1(k)、Yf(k)分别为k采样时刻塔顶和进料板处的汽相轻组分浓度,分别为k采样时刻精馏段和提馏段的温度分布曲线位置变化速度;
5)预测k+1采样时刻的温度分布,从历史数据库中获取曲线位置变化速度数据,温度分布函数拟合参数,采样周期值,采用式(7)-(10)得到:
其中,δ为采样周期,Sr(k)、Ss(k)分别为k采样时刻精馏段、提馏段温度分布曲线位置,Sr(k+1)、Ss(k+1)分别为k+1采样时刻精馏段、提馏段温度分布曲线位置,为K+1时刻精馏段、提馏段段第i块塔板温度的预测值;
6)从历史数据库中获取当前温度测量值,以及预测值,计算k采样时刻的温度预测误差过程如式(11):
7)修正温度预测值,采用式(12)得到:
4.如权利要求3温度非线性过程观测方法,其特征在于:所述的历史数据库为DCS系统的存储装置,控制站读取历史数据库,显示观测器工作过程状态。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010102137532A CN101890246B (zh) | 2010-06-30 | 2010-06-30 | 一种精馏塔温度非线性观测系统及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010102137532A CN101890246B (zh) | 2010-06-30 | 2010-06-30 | 一种精馏塔温度非线性观测系统及方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101890246A true CN101890246A (zh) | 2010-11-24 |
CN101890246B CN101890246B (zh) | 2012-06-20 |
Family
ID=43099687
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2010102137532A Expired - Fee Related CN101890246B (zh) | 2010-06-30 | 2010-06-30 | 一种精馏塔温度非线性观测系统及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101890246B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107885081A (zh) * | 2017-11-13 | 2018-04-06 | 浙江大学 | 一种针对产品浓度曲线的内部热耦合空分塔在线观测器 |
CN108345953A (zh) * | 2017-01-24 | 2018-07-31 | 上海温尔信息科技有限公司 | 温度预测方法及装置 |
CN111537679A (zh) * | 2020-05-25 | 2020-08-14 | 北京化工大学 | 一种基于非线性滤波的隔离壁精馏塔浓度软测量方法 |
CN114154706A (zh) * | 2021-11-29 | 2022-03-08 | 江南大学 | 精馏过程单吨能耗估计精度评估方法及系统 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4419188A (en) * | 1980-06-02 | 1983-12-06 | Mccall Thomas F | Thermally coupled extractive distillation process |
US6347533B1 (en) * | 2000-03-04 | 2002-02-19 | Peter Tung | Hydraulically balanced fully thermally coupled system |
CN200955890Y (zh) * | 2006-09-08 | 2007-10-03 | 浙江大学 | 一种内部热耦合空分装置 |
CN101364355A (zh) * | 2008-09-26 | 2009-02-11 | 东南大学 | 整体煤气化联合循环空气分离系统的建模与虚拟仿真方法 |
-
2010
- 2010-06-30 CN CN2010102137532A patent/CN101890246B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4419188A (en) * | 1980-06-02 | 1983-12-06 | Mccall Thomas F | Thermally coupled extractive distillation process |
US6347533B1 (en) * | 2000-03-04 | 2002-02-19 | Peter Tung | Hydraulically balanced fully thermally coupled system |
CN200955890Y (zh) * | 2006-09-08 | 2007-10-03 | 浙江大学 | 一种内部热耦合空分装置 |
CN101364355A (zh) * | 2008-09-26 | 2009-02-11 | 东南大学 | 整体煤气化联合循环空气分离系统的建模与虚拟仿真方法 |
Non-Patent Citations (1)
Title |
---|
《化工学报》 20000630 刘兴高等 刘兴高等 内部热耦合精馏塔的初步设计(I) 第421-424页 1-4 , 第03期 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108345953A (zh) * | 2017-01-24 | 2018-07-31 | 上海温尔信息科技有限公司 | 温度预测方法及装置 |
CN108345953B (zh) * | 2017-01-24 | 2021-12-17 | 上海温尔信息科技有限公司 | 温度预测方法及装置 |
CN107885081A (zh) * | 2017-11-13 | 2018-04-06 | 浙江大学 | 一种针对产品浓度曲线的内部热耦合空分塔在线观测器 |
CN111537679A (zh) * | 2020-05-25 | 2020-08-14 | 北京化工大学 | 一种基于非线性滤波的隔离壁精馏塔浓度软测量方法 |
CN111537679B (zh) * | 2020-05-25 | 2021-06-08 | 北京化工大学 | 一种基于非线性滤波的隔离壁精馏塔浓度软测量方法 |
CN114154706A (zh) * | 2021-11-29 | 2022-03-08 | 江南大学 | 精馏过程单吨能耗估计精度评估方法及系统 |
CN114154706B (zh) * | 2021-11-29 | 2022-07-15 | 江南大学 | 精馏过程单吨能耗估计精度评估方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN101890246B (zh) | 2012-06-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101879378B (zh) | 内部热耦合精馏塔温度非线性观测系统及方法 | |
CN101890246B (zh) | 一种精馏塔温度非线性观测系统及方法 | |
CN101901005B (zh) | 空分节能过程的比例积分控制系统及方法 | |
CN102466412A (zh) | 一种多温区炉综合控制系统及方法 | |
JP6003247B2 (ja) | 情報処理装置、情報処理装置の制御方法及びプログラム | |
CN104635493A (zh) | 基于温度波模型预测控制的内部热耦合精馏控制装置 | |
CN101887262B (zh) | 内部热耦合精馏塔的非线性模型预测控制系统及方法 | |
CN101763037B (zh) | 一种空分节能过程的非线性预测控制系统及方法 | |
CN104587695A (zh) | 基于温度波特性的内部热耦合精馏塔控制装置 | |
CN101708373B (zh) | 一种内部热耦合精馏塔高纯非线性观测系统及方法 | |
CN101788810B (zh) | 一种内部热耦合精馏过程的非线性预测控制系统及方法 | |
CN101884849B (zh) | 一种高纯精馏过程的浓度非线性观测系统及方法 | |
CN102435334B (zh) | 一种变压器绕组温度测量方法 | |
CN1962015B (zh) | 高纯精馏的动态矩阵控制系统和方法 | |
CN101763081B (zh) | 一种空分节能过程的非线性过程观测系统及方法 | |
CN105549533A (zh) | 数控机床温度场在线监测及智能预警系统 | |
CN101884848B (zh) | 一种空分节能过程温度分布的非线性观测系统及方法 | |
CN101890247B (zh) | 一种内部热耦合精馏塔的高纯非线性控制系统及方法 | |
CN101887261B (zh) | 一种内部热耦合精馏塔的自适应非线性控制系统及方法 | |
CN101708379A (zh) | 一种高效节能精馏过程的高纯控制系统及方法 | |
CN104606912B (zh) | 基于温度波特性的内部热耦合精馏在线观测器 | |
CN101881962B (zh) | 一种内部热耦合精馏塔的超高纯非线性控制系统及方法 | |
CN103413029B (zh) | 具有多速率采样连续搅拌釜式反应器的滚动时域估计方法 | |
CN101881964B (zh) | 一种内部热耦合精馏塔的高纯自适应非线性控制系统及方法 | |
Fu et al. | Nonlinear wave modeling and dynamic analysis of high-purity heat integrated air separation column |
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 | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20120620 Termination date: 20130630 |