一种空分节能过程温度分布的非线性观测系统及方法

专利名称：一种空分节能过程温度分布的非线性观测系统及方法
技术领域：
本发明涉及空分动态建模领域，特别地，涉及空分节能过程温度分布的非线性观 测器系统设计及其方法。
背景技术：
空分是对空气进行分离，得到氧、氮、氩等高纯工业气体的国民经济重要行业，其 产品广泛用于石油、化工、冶金、电子、能源、航空航天、食品饮料、医疗保健等各种工业领 域。巨大的能量消耗一直为空分行业的瓶颈问题。在空分过程中，精馏塔是一个重要 的传质单元操作过程，也是主要的耗能单元。由于空分精馏过程的强非线性，耦合性等复杂 动态特性，传统的机理模型虽然具有较高的预测精度，但是在线运算效率太低，较难应用于 现今较流行的先进控制方案，而目前商业使用的基于数据的线性辨识模型，往往具有较低 的预测精度，对控制器的设计提出了较高的要求，从而建立空分精馏过程的高效率高精度 非线性模型实现空分精馏过程中的温度分布等状态预测已经成为一项关键的空分节能技 术。

发明内容
为了克服现有的空分节能过程温度分布的控制系统的预测精度低、精度较低、效 率较低的不足，本发明的提供一种预测精度高、高效率和高精度的空分节能过程温度分布 的非线性观测系统及方法。本发明解决其技术问题所采用的技术方案是一种空分节能过程温度分布的非线性观测系统，包括与空分塔直接连接的现场智 能仪表和DCS系统，所述DCS系统包括存储装置和上位机，所述现场智能仪表与存储装置、 上位机相连，所述上位机包括用以预测空分塔将来时刻的温度分布的观测器，所述观测器 包括温度分布函数拟合模块，用以采用历史数据库中的各块塔板温度Ti拟合温度分布 函数，并将拟合参数存储到历史数据库当中，通过式(1) (2)得到 其中，i为塔板编号，f为进料板编号，η为上塔塔底编号， 为第i块塔板处预估温 度，Tfflin,r> Tmax,r、kr, Tfflinjs, Tfflaxjs, ks为拟合参数，Sr、Ss分别为空分塔下塔、上塔液相温度分 布曲线的位置；将来时刻空分塔温度预测模块，用以根据当前温度数据以及温度分布函数预测将 来时刻的温度分布，该模块包括以下过程
其中，k为当前采样时刻，下角标i为塔板编号，1为下塔塔顶编号，η为上塔塔底 编号，下角标N、O分别代表氮和氧，上角标f代表进料，Xi,N、Xi, 0分别为空分塔中第i块塔 板处氮、氧的液相组分浓度，P为空分塔压强，Ti为第i块塔板处的温度，α N、C^分别为氮 和氧相对于氩的相对挥发度，aN、bN、cN、a。、b0, C0为安东尼常数，Fi (k)为k采样时刻第i ±夬 塔板进料流量，Ln为上塔塔底液相流率，V1为下塔塔顶汽相流率，xn,N(k)、χη, ο(k)分别为k 采样时刻上塔塔底氮、氧的液相组分浓度，Xl,N(k)、kdk)分别为k采样时刻下塔塔顶氮、 氧的液相组分浓度，yf,^k)分别为k采样时刻进料板氧汽相组分浓度，yuGO^uGO分别 为k采样时刻下塔塔顶氮、氧汽相组分浓度，Xfi,N(k)，/"(ΙΟ分别为第i块塔板的进料液
氮组分浓度和进料液氧组分浓度，M为塔板持液量，，分别为下塔上塔温度分
at at
布位置变化速度；②下一时刻波形位置以及下一时刻温度分布，通过式(10)_(13)得到
"(9)得至
①当前时刻温度分布曲线位置变化速度计算，通过式(3) 其中，δ为采样周期，Sr(k),Ss(k)分别为k采样时刻下塔、上塔温度分布曲线位 置，Sr(k+1)、Ss(k+l)分别为k+1采样时刻下塔、上塔温度分布曲线位置，+ 为k+Ι时刻 下塔或上塔第i块塔板温度的预测值；温度预测误差补偿模块，用以从历史数据库中获取当前温度测量值，以及预测值， 计算预测误差，修正温度预测值，通过式(14) (15)得到Ti(U) = TXk)-T1(U) i = 1,2, ......，η (14)f{^k + l) = TXk + !)+ Ti(Ji)I = 1,2, ......，η (15)其中，TiGO为k+Ι采样时刻第i块塔板的温度，r,(k)为k采样时刻第i块塔板 的温度预测误差，辩A+1)为k+Ι采样时刻第i块塔板的温度预测修正值。作为优选的一种方案所述的现场智能仪表与DCS系统中的存储装置、上位机通 过现场总线相互连接。再进一步，所述的上位机还包括人机界面模块，用以设定采样周期，设定预测时间 长度，显示观测器输出的将来时刻温度预测值和当前时刻温度测量值。一种空分节能过程温度分布的非线性观测方法，所述的观测方法包括以下步骤1)确定采样周期，并将氮和氧相对于氩的相对挥发度aN、a0，安东尼常数aN、bN、 cN、aQ、bQ、cQ保存在历史数据库当中；2)从现场智能仪表获取k采样时刻压强P，各塔板温度Ti,计算第i块塔板处氮、 氧的液相组分浓度值，通过式(3)-(4)得到 其中，Xi,Ν为空分塔第i块塔板处氮的液相组分浓度，Xy为空分塔中第i块塔板 处氧的液相组分浓度；检测k采样时刻下塔塔顶汽相流量和上塔塔底液相流量值，同上塔压强数据，各 塔板温度数据，组分浓度的测量值一起存储到观测器系统的历史数据库当中；3)从历史数据库当中获取k采样时刻温度的测量值，拟合k采样时刻的温度分布 函数如式(1) (2) 其中，1为下塔塔顶编号，f为进料板编号，η为上塔塔底编号，7；为第i块塔板处预 估温度，Tfflinjr, Tmax,r、kr, Tfflinjs, Tfflaxjs, ks为拟合参数，Sr、Ss分别为空分塔下塔、上塔液相温 度分布曲线的位置；4)计算k采样时刻的下塔上塔温度分布位置变化速度，通过式(5)-(9)得到
8 其中，k为当前采样时刻，下角标i为塔板编号，1为下塔塔顶编号，η为上塔塔底 的编号，下角标N、0分别代表氮和氧，上角标f代表进料，Xi,N、Xi, ο分别为空分塔中第i块 塔板处氮、氧的液相组分浓度，P为空分塔压强，Ti为第i块塔板处的温度，α N、C^分别为 氮和氧相对于氩的相对挥发度，aN、bN、cN、a。、b0, C0为安东尼常数，Fi (k)为k采样时刻第i 块塔板进料流量，Ln为上塔塔底液相流率，V1为下塔塔顶汽相流率，xn,N(k)、xn,0(k)分别为 k采样时刻上塔塔底氮、氧的液相组分浓度，Xl,N(k)、Xu(k)分别为k采样时刻下塔塔顶氮、 氧的液相组分浓度，yf,^k)分别为k采样时刻进料板氧汽相组分浓度，yuGO^uGO分别 为k采样时刻下塔塔顶氮、氧汽相组分浓度，Xfi,N(k)，/"(ΙΟ分别为第i块塔板的进料液
氮组分浓度和进料液氧组分浓度，M为塔板持液量，，分别为下塔上塔温度分
dt dt
布位置变化速度；5)下一时刻波形位置以及下一时刻温度分布，通过式(10)-(13)得到 其中，δ为采样周期，Sr(k),Ss(k)分别为k采样时刻下塔，上塔温度分布曲线位 置，Sr(k+1)、Ss(k+l)分别为k+1采样时刻下塔，上塔温度分布曲线位置，+ 为k+Ι时刻 下塔或上塔的第i块塔板温度的预测值；6)从历史数据库中获取当前温度测量值，以及预测值，计算预测误差，修正温度预 测值，通过式(14) (15)得到 r^k) = TXk)-TXk)! = 1,2,
·, η
(14)
f{{k + l) = TXk + l) + rXk)i = 1,2, ......，η (15)其中，TiGO为k+Ι采样时刻第i块塔板的温度，r,(k)为k采样时刻第i块塔板 的温度预测误差，辩A+1)为k+Ι采样时刻的温度预测修正值；进一步，所述历史数据库为DCS系统的存储装置，控制站读取历史数据库，显示观 测器工作过程状态。本发明的有益效果主要表现在1.基于空分机理建立高精度的简化预测模型，能 够高效率地准确地跟踪空分节能过程温度的复杂非线性变化特征；2.在线运行效率高，温 度状态变量预测比传统基于数据的辨识模型精确。


图1是本发明所提出的空分节能过程温度非线性观测系统的结构图。图2是上位机观测器实现方法的原理图。
具体实施例方式下面根据附图具体说明本发明。实施例1参照图1、图2，一种空分节能过程温度分布的非线性观测系统，包括与空分塔1直 接连接的现场智能仪表2和DCS系统13，所述DCS系统13包括存储装置4、控制站5及上 位机6，所述现场智能仪表2与存储装置4、控制站5和上位机6依次相连，上位机包括用以 预测空分塔将来时刻的温度分布的观测器，所述观测器包括温度分布函数拟合模块9，将来 时刻温度预测模块11，温度预测误差补偿模块10。所述的温度分布函数拟合模块9，采用历史数据库中的各块塔板温度Ti拟合温度 分布函数，并将拟合参数存储到历史数据库当中，通过式(1) (2)得到=……，f-1 (1)T1 = T^s + ^^ i = f, f+1,……，η (2)其中，1为下塔塔顶编号，f为进料板编号，η为上塔塔底编号， 为第i块塔板处预 估温度，Tmin,ρ ρ kr, Tfflinjs, Tfflaxjs, ks为拟合参数，民、Ss分别为空分塔下塔、上塔液相温 度分布曲线的位置。所述的将来时刻空分塔温度预测模块11，根据当前温度数据以及温度分布函数预 测将来时刻的温度分布，该模块包括以下几部分①当前时刻温度分布曲线位置变化速度计算，通过式(3)-(9)得到γ _ PxaNxW Ν(3)
Λλγ (众)=Λ0 ⑷=yf,0(k) =
其中，k为当前采样时刻，下角标i为塔板编号，1为下塔塔顶编号，η为上塔塔底 的编号，下角标N、0分别代表氮和氧，上角标f代表进料，Xi,N、Xi, ο分别为空分塔中第i 土夬 塔板处氮、氧的液相组分浓度，P为空分塔压强，Ti为第i块塔板处的温度，α N、C^分别为 氮和氧相对于氩的相对挥发度，aN、bN、cN、a。、bQ、C0为安东尼(Antonie)常数，Fi (k)为k采 样时刻第i块塔板进料流量，Ln为上塔塔底液相流率，V1为下塔塔顶汽相流率，xn,N(k)、xn, 。(k)分别为k采样时刻上塔塔底氮、氧的液相组分浓度，Xl,N(k),Xl,0(k)分别为k采样时刻 下塔塔顶氮、氧的液相组分浓度，yf,^k)分别为k采样时刻进料板氧汽相组分浓度，yi,N(k)、 yi,0(k)分别为k采样时刻下塔塔顶氮、氧汽相组分浓度，Xfi,N(k)，Xfi,0(k)分别为第i块塔
板的进料液氮组分浓度和进料液氧组分浓度，M为塔板持液量，，分别为下塔、
at at
上塔温度分布位置变化速度；
②下一时刻波形位置以及下一时刻温度分布，通过式(10)_(13)得到
Γ1Sr(k + \)^Sr(k) + ^-(k)x5(10)
dt
1Γ1Ss(k + \)^Ss(k) + ^(k)x5(11)
dt^ + = ^ +i = l,2,……，f-1 (12)T,(k + l) = Tmm:S +i = f,f+l,……，η (13)其中，δ为采样周期，Sr(k),Ss(k)分别为k采样时刻下塔、上塔温度分布曲线位 置，Sr(k+1)、Ss(k+l)分别为k+1采样时刻下塔，上塔温度分布曲线位置，+ 为k+Ι时亥Ij 下塔或上塔的第i块塔板温度的预测值；所述的温度预测误差补偿模块10，从历史数据库中获取当前温度测量值，以及预 测值，计算预测误差，修正温度预测值，通过式(14) (15)得到rXk) = ^)-^)! = 1,2,……，n (14)
其中，TiGO为k+Ι采样时刻第i块塔板的温度，r,(k)为k采样时刻第i块塔板 的温度预测误差，辩A+1)为k+Ι采样时刻第i块塔板的温度预测修正值。所述的现场智能仪表2、DCS系统中的存储装置4.、控制站5及上位机6通过现场 总线和数据接口 3依次相连；所述的智能仪表用以实现检测模块7的功能即检测得到空分 塔的温度、压强数据，所述的检测模块7包括温度检测元件和压力检测元件，所述的数据接 口用以实现I/O模块8的功能即数据的输入输出。所述的上位机包括人机界面模块12，用于设定采样周期，设定预测时间长度，显示 观测器输出的将来时刻温度预测值和当前时刻温度测量值。实施例2参照图1和图2，一种空分节能过程温度分布的非线性观测方法，包括以下步骤1)确定采样周期，并将氮和氧相对于氩的相对挥发度αΝ、α ^，安东尼(Antonie) 常数％、％、 、如、13。、(^保存在历史数据库当中；2)从现场智能仪表获取k采样时刻压强P，各塔板温度Ti,计算第i块塔板处氮、 氧的液相组分浓度值，通过式(3) _(4)得到 其中，Xi,Ν为空分塔第i块塔板处氮的液相组分浓度，Xu为空分塔中第i块塔板 处氧的液相组分浓度；检测k采样时刻下塔塔顶汽相流量和上塔塔底液相流量值，同上塔压强数据，各 塔板温度数据，组分浓度的测量值一起存储到观测器系统的历史数据库当中；3)从历史数据库当中获取k采样时刻温度的测量值，拟合k采样时刻的温度分布 函数如式(1) (2)7； 二 Tmm,+^Jifi = I, 2，……，f_l (1)i = f，f+1，……，η (2)其中，1为下塔塔顶编号，f为进料板编号，η为上塔塔底编号， 为第i块塔板处预 估温度，Tfflinjr, Tmax,r、kr, Tfflinjs, Tfflaxjs, ks为拟合参数，Sr、Ss分别为空分塔下塔、上塔液相温 度分布曲线的位置；4)计算k采样时刻的下塔上塔温度分布位置变化速度，通过式(5)_(9)得到
Γ π“、aNXl,N (^)/ C ΛΛλγ (众)=7-h "、丄,⑶AoW= 二 (6)
12 其中，k为当前采样时刻，下角标i为塔板编号，1为下塔塔顶编号，η为上塔塔底 的编号，下角标N、0分别代表氮和氧，上角标f代表进料，Xi,N、Xi, ο分别为空分塔中第i 土夬 塔板处氮、氧的液相组分浓度，P为空分塔压强，Ti为第i块塔板处的温度，α N、C^分别为 氮和氧相对于氩的相对挥发度，aN、bN、cN、a。、b0, C0为安东尼常数，Fi (k)为k采样时刻第i 块塔板进料流量，Ln为上塔塔底液相流率，V1为下塔塔顶汽相流率，xn,N(k)、xn,0(k)分别为 k采样时刻上塔塔底氮、氧的液相组分浓度，Xl,N(k)、Xu(k)分别为k采样时刻下塔塔顶氮、 氧的液相组分浓度，yf,^k)分别为k采样时刻进料板氧汽相组分浓度，yuGO^uGO分别 为k采样时刻下塔塔顶氮、氧汽相组分浓度，，Xfi,N(k)，Xf^0(k)分别为第i块塔板的进料液
氮组分浓度和进料液氧组分浓度，M为塔板持液量，^KQ，分别为下塔上塔温度分
at at
布位置变化速度；5)下一时刻波形位置以及下一时刻温度分布，通过式(10)-(13)得到
Γ1Sr{k + \)^Sr(k) + ^-(k)x5(10)
dt
Γ1Ss{k + \)^Ss(k) + ^(k)x5(11)
dtm + = += ……，f-1 (12)fiik + X) = T^s + ^^y)i = f, f+l，……，n (13)其中，δ为采样周期，Sr(k),Ss(k)分别为k采样时刻下塔，上塔温度分布曲线位 置，Sr(k+1)、Ss(k+l)分别为k+1采样时刻下塔，上塔温度分布曲线位置，iG + l)为k+Ι时亥Ij 下塔或上塔的第i块塔板温度的预测值；6)从历史数据库中获取当前温度测量值，以及预测值，计算预测误差，修正温度预 测值，通过式(14) (15)得到rt(k) = T^k)-T1(U)I = 1,2,……，η (14)
2, ......，η (15)其中，ri(k)为k采样时刻第i块塔板的温度预测误差，ff^ + l)为k+Ι采样时刻的 温度预测修正值。所述的历史数据库为DCS系统中的存储装置，所述的DCS系统包括数据接口，存储装置，控制站，其中控制站可以读取历史数据库，显示观测器工作过程状态。 上述实施例用来解释说明本发明，而不是对本发明进行限制，在本发明的精神和
权利要求的保护范围内，对本发明作出的任何修改和改变，都落入本发明的保护范围。
权利要求
一种空分节能过程温度分布的非线性观测系统，包括与空分塔直接连接的现场智能仪表和DCS系统，所述DCS系统包括存储装置、控制站和上位机，所述现场智能仪表与存储装置、控制站和上位机相连，其特征在于所述上位机包括用以预测空分塔将来时刻的温度分布的观测器，所述观测器包括温度分布函数拟合模块，用以采用历史数据库中的各块塔板温度Ti拟合温度分布函数，并将拟合参数存储到历史数据库当中，通过式(1)(2)得到 <mrow><msub> <mover><mi>T</mi><mo>^</mo> </mover> <mi>i</mi></msub><mo>=</mo><msub> <mi>T</mi> <mrow><mi>min</mi><mo>,</mo><mi>r</mi> </mrow></msub><mo>+</mo><mfrac> <mrow><msub> <mi>T</mi> <mrow><mi>max</mi><mo>,</mo><mi>r</mi> </mrow></msub><mo>-</mo><msub> <mi>T</mi> <mrow><mi>min</mi><mo>,</mo><mi>r</mi> </mrow></msub> </mrow> <msup><mrow> <mn>1</mn> <mo>+</mo> <mi>e</mi></mrow><mrow> <mo>-</mo> <msub><mi>k</mi><mi>r</mi> </msub> <mrow><mo>(</mo><mi>i</mi><mo>-</mo><msub> <mi>S</mi> <mi>r</mi></msub><mo>)</mo> </mrow></mrow> </msup></mfrac> </mrow>i＝1，2，……，f-1(1) <mrow><msub> <mover><mi>T</mi><mo>^</mo> </mover> <mi>i</mi></msub><mo>=</mo><msub> <mi>T</mi> <mrow><mi>min</mi><mo>,</mo><mi>s</mi> </mrow></msub><mo>+</mo><mfrac> <mrow><msub> <mi>T</mi> <mrow><mi>max</mi><mo>,</mo><mi>s</mi> </mrow></msub><mo>-</mo><msub> <mi>T</mi> <mrow><mi>min</mi><mo>,</mo><mi>s</mi> </mrow></msub> </mrow> <msup><mrow> <mn>1</mn> <mo>+</mo> <mi>e</mi></mrow><mrow> <mo>-</mo> <msub><mi>k</mi><mi>s</mi> </msub> <mrow><mo>(</mo><mi>i</mi><mo>-</mo><msub> <mi>S</mi> <mi>s</mi></msub><mo>)</mo> </mrow></mrow> </msup></mfrac> </mrow>i＝f，f+1，……，n(2)其中，i为塔板编号，f为进料板编号，n为上塔塔底编号，为第i块塔板处预估温度，Tmin，r、Tmax，r、kr、Tmin，s、Tmax，s、ks为拟合参数，Sr、Ss分别为空分塔下塔、上塔液相温度分布曲线的位置；将来时刻空分塔温度预测模块，用以根据当前温度数据以及温度分布函数预测将来时刻的温度分布，该模块包括以下过程①当前时刻温度分布曲线位置变化速度计算，通过式(3)-(9)得到 <mrow><msub> <mi>X</mi> <mrow><mi>i</mi><mo>,</mo><mi>N</mi> </mrow></msub><mo>=</mo><mfrac> <mrow><mi>P</mi><mo>&times;</mo><msub> <mi>&alpha;</mi> <mi>N</mi></msub><mo>&times;</mo><msup> <mn>10</mn> <mrow><mrow> <mo>(</mo> <mfrac><mrow> <msub><mi>T</mi><mi>i</mi> </msub> <mo>+</mo> <msub><mi>c</mi><mi>N</mi> </msub></mrow><msub> <mi>b</mi> <mi>N</mi></msub> </mfrac> <mo>-</mo> <msub><mi>a</mi><mi>N</mi> </msub> <mo>)</mo></mrow> </mrow></msup><mo>-</mo><mn>1</mn> </mrow> <mrow><msub> <mi>&alpha;</mi> <mi>N</mi></msub><mo>-</mo><mn>1</mn> </mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo></mrow> </mrow> <mrow><msub> <mi>X</mi> <mrow><mi>i</mi><mo>,</mo><mi>O</mi> </mrow></msub><mo>=</mo><mfrac> <mrow><mi>P</mi><mo>&times;</mo><msub> <mi>&alpha;</mi> <mi>O</mi></msub><mo>&times;</mo><msup> <mn>10</mn> <mrow><mrow> <mo>(</mo> <mfrac><mrow> <msub><mi>T</mi><mi>i</mi> </msub> <mo>+</mo> <msub><mi>c</mi><mi>O</mi> </msub></mrow><msub> <mi>b</mi> <mi>O</mi></msub> </mfrac> <mo>-</mo> <msub><mi>a</mi><mi>O</mi> </msub> <mo>)</mo></mrow> </mrow></msup><mo>-</mo><mn>1</mn> </mrow> <mrow><msub> <mi>&alpha;</mi> <mi>O</mi></msub><mo>-</mo><mn>1</mn> </mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo></mrow> </mrow> <mrow><msub> <mi>y</mi> <mrow><mn>1</mn><mo>,</mo><mi>N</mi> </mrow></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><mo>=</mo><mfrac> <mrow><msub> <mi>&alpha;</mi> <mi>N</mi></msub><msub> <mi>x</mi> <mrow><mn>1</mn><mo>,</mo><mi>N</mi> </mrow></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow> </mrow> <mrow><mrow> <mo>(</mo> <msub><mi>&alpha;</mi><mi>N</mi> </msub> <mo>-</mo> <mn>1</mn> <mo>)</mo></mrow><msub> <mi>x</mi> <mrow><mn>1</mn><mo>,</mo><mi>N</mi> </mrow></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><mo>+</mo><mn>1</mn> </mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo></mrow> </mrow> <mrow><msub> <mi>y</mi> <mrow><mn>1</mn><mo>,</mo><mi>O</mi> </mrow></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><mo>=</mo><mfrac> <mrow><msub> <mi>&alpha;</mi> <mi>O</mi></msub><msub> <mi>x</mi> <mrow><mn>1</mn><mo>,</mo><mi>O</mi> </mrow></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow> </mrow> <mrow><mrow> <mo>(</mo> <msub><mi>&alpha;</mi><mi>O</mi> </msub> <mo>-</mo> <mn>1</mn> <mo>)</mo></mrow><msub> <mi>x</mi> <mrow><mn>1</mn><mo>,</mo><mi>O</mi> </mrow></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><mo>+</mo><mn>1</mn> </mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo></mrow> </mrow> <mrow><msub> <mi>y</mi> <mrow><mi>f</mi><mo>,</mo><mi>O</mi> </mrow></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><mo>=</mo><mfrac> <mrow><msub> <mi>&alpha;</mi> <mi>O</mi></msub><msub> <mi>x</mi> <mrow><mi>f</mi><mo>,</mo><mi>O</mi> </mrow></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow> </mrow> <mrow><mrow> <mo>(</mo> <msub><mi>&alpha;</mi><mi>O</mi> </msub> <mo>-</mo> <mn>1</mn> <mo>)</mo></mrow><msub> <mi>x</mi> <mrow><mi>f</mi><mo>,</mo><mi>O</mi> </mrow></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><mo>+</mo><mn>1</mn> </mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo></mrow> </mrow> <mrow><mfrac> <msub><mi>dS</mi><mi>r</mi> </msub> <mi>dt</mi></mfrac><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><mo>=</mo><mfrac> <mrow><mo>[</mo><mo>-</mo><msub> <mi>V</mi> <mn>1</mn></msub><msub> <mi>y</mi> <mrow><mn>1</mn><mo>,</mo><mi>N</mi> </mrow></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><mo>-</mo><msub> <mi>L</mi> <mi>n</mi></msub><msub> <mi>x</mi> <mrow><mi>n</mi><mo>,</mo><mi>N</mi> </mrow></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><mo>+</mo><munderover> <mi>&Sigma;</mi> <mrow><mi>i</mi><mo>=</mo><mn>1</mn> </mrow> <mi>n</mi></munderover><msub> <mi>F</mi> <mi>i</mi></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><msub> <msup><mi>x</mi><mi>f</mi> </msup> <mrow><mi>i</mi><mo>,</mo><mi>N</mi> </mrow></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><mo>]</mo><mrow> <mo>(</mo> <msub><mi>&alpha;</mi><mi>N</mi> </msub> <mo>-</mo> <mn>1</mn> <mo>)</mo></mrow><msup> <mrow><mo>(</mo><mn>1</mn><mo>+</mo><msup> <mi>e</mi> <mrow><mo>-</mo><msub> <mi>k</mi> <mi>r</mi></msub><mrow> <mo>(</mo> <mi>i</mi> <mo>-</mo> <msub><mi>s</mi><mi>r</mi> </msub> <mo>)</mo></mrow> </mrow></msup><mo>)</mo> </mrow> <mn>2</mn></msup> </mrow> <mrow><mi>M</mi><mrow> <mo>(</mo> <msub><mi>x</mi><mrow> <mi>n</mi> <mo>,</mo> <mi>N</mi></mrow> </msub> <mrow><mo>(</mo><mi>k</mi><mo>)</mo> </mrow> <mo>-</mo> <msub><mi>x</mi><mrow> <mn>1</mn> <mo>,</mo> <mi>N</mi></mrow> </msub> <mrow><mo>(</mo><mi>k</mi><mo>)</mo> </mrow> <mo>)</mo></mrow><msub> <mi>P&alpha;</mi> <mi>N</mi></msub><mrow> <mo>(</mo> <msub><mi>T</mi><mrow> <mi>max</mi> <mo>,</mo> <mi>r</mi></mrow> </msub> <mo>-</mo> <msub><mi>T</mi><mrow> <mi>min</mi> <mo>,</mo> <mi>r</mi></mrow> </msub> <mo>)</mo></mrow><msub> <mi>k</mi> <mi>r</mi></msub><msup> <mi>e</mi> <mrow><mo>-</mo><msub> <mi>k</mi> <mi>r</mi></msub><mrow> <mo>(</mo> <mi>i</mi> <mo>-</mo> <msub><mi>s</mi><mi>r</mi> </msub> <mo>)</mo></mrow> </mrow></msup><munderover> <mi>&Sigma;</mi> <mrow><mi>i</mi><mo>=</mo><mn>1</mn> </mrow> <mrow><mi>f</mi><mo>-</mo><mn>1</mn> </mrow></munderover><mfrac> <mrow><mi>ln</mi><mn>10</mn> </mrow> <mrow><msub> <mi>T</mi> <mi>i</mi></msub><mo>+</mo><msub> <mi>c</mi> <mi>N</mi></msub><mo>-</mo><msub> <mi>a</mi> <mi>N</mi></msub><msub> <mi>b</mi> <mi>N</mi></msub> </mrow></mfrac> </mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow> <mo>(</mo> <mn>8</mn> <mo>)</mo></mrow> </mrow> <mrow><mfrac> <msub><mi>dS</mi><mi>s</mi> </msub> <mi>dt</mi></mfrac><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><mo>=</mo><mfrac> <mrow><mo>[</mo><mo>-</mo><msub> <mi>V</mi> <mi>f</mi></msub><msub> <mi>y</mi> <mrow><mi>f</mi><mo>,</mo><mi>O</mi> </mrow></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><mo>-</mo><msub> <mi>L</mi> <mi>n</mi></msub><msub> <mi>x</mi> <mrow><mi>n</mi><mo>,</mo><mi>O</mi> </mrow></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><mo>+</mo><munderover> <mi>&Sigma;</mi> <mrow><mi>i</mi><mo>=</mo><mn>1</mn> </mrow> <mi>n</mi></munderover><msub> <mi>F</mi> <mi>i</mi></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><msub> <msup><mi>x</mi><mi>f</mi> </msup> <mrow><mi>i</mi><mo>,</mo><mi>O</mi> </mrow></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><mo>]</mo><mrow> <mo>(</mo> <msub><mi>&alpha;</mi><mi>O</mi> </msub> <mo>-</mo> <mn>1</mn> <mo>)</mo></mrow><msup> <mrow><mo>(</mo><mn>1</mn><mo>+</mo><msup> <mi>e</mi> <mrow><mo>-</mo><msub> <mi>k</mi> <mi>s</mi></msub><mrow> <mo>(</mo> <mi>i</mi> <mo>-</mo> <msub><mi>s</mi><mi>s</mi> </msub> <mo>)</mo></mrow> </mrow></msup><mo>)</mo> </mrow> <mn>2</mn></msup> </mrow> <mrow><mi>M</mi><mrow> <mo>(</mo> <msub><mi>x</mi><mrow> <mi>n</mi> <mo>,</mo> <mi>O</mi></mrow> </msub> <mrow><mo>(</mo><mi>k</mi><mo>)</mo> </mrow> <mo>-</mo> <msub><mi>x</mi><mrow> <mn>1</mn> <mo>,</mo> <mi>O</mi></mrow> </msub> <mrow><mo>(</mo><mi>k</mi><mo>)</mo> </mrow> <mo>)</mo></mrow><msub> <mi>P&alpha;</mi> <mi>O</mi></msub><mrow> <mo>(</mo> <msub><mi>T</mi><mrow> <mi>max</mi> <mo>,</mo> <mi>s</mi></mrow> </msub> <mo>-</mo> <msub><mi>T</mi><mrow> <mi>min</mi> <mo>,</mo> <mi>s</mi></mrow> </msub> <mo>)</mo></mrow><msub> <mi>k</mi> <mi>s</mi></msub><msup> <mi>e</mi> <mrow><mo>-</mo><msub> <mi>k</mi> <mi>s</mi></msub><mrow> <mo>(</mo> <mi>i</mi> <mo>-</mo> <msub><mi>s</mi><mi>s</mi> </msub> <mo>)</mo></mrow> </mrow></msup><munderover> <mi>&Sigma;</mi> <mrow><mi>i</mi><mo>=</mo><mi>f</mi> </mrow> <mi>n</mi></munderover><mfrac> <mrow><mi>ln</mi><mn>10</mn> </mrow> <mrow><msub> <mi>T</mi> <mi>i</mi></msub><mo>+</mo><msub> <mi>c</mi> <mi>O</mi></msub><mo>-</mo><msub> <mi>a</mi> <mi>O</mi></msub><msub> <mi>b</mi> <mi>O</mi></msub> </mrow></mfrac> </mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow> <mo>(</mo> <mn>9</mn> <mo>)</mo></mrow> </mrow>其中，k为当前采样时刻，下角标i为塔板编号，1为下塔塔顶编号，n为上塔塔底编号，下角标N、O分别代表氮和氧，上角标f代表进料，Xi，N、Xi，O分别为空分塔中第i块塔板处氮、氧的液相组分浓度，P为空分塔压强，Ti为第i块塔板处的温度，αN、αO分别为氮和氧相对于氩的相对挥发度，aN、bN、cN、aO、bO、cO为安东尼常数，Fi(k)为k采样时刻第i块塔板进料流量，Ln为上塔塔底液相流率，V1为下塔塔顶汽相流率，xn，N(k)、xn，O(k)分别为k采样时刻上塔塔底氮、氧的液相组分浓度，xl，N(k)、xl，O(k)分别为k采样时刻下塔塔顶氮、氧的液相组分浓度，yf，O(k)分别为k采样时刻进料板氧汽相组分浓度，y1，N(k)、y1，O(k)分别为k采样时刻下塔塔顶氮、氧汽相组分浓度，xfi，N(k)，xfi，O(k)分别为第i块塔板的进料液氮组分浓度和进料液氧组分浓度，M为塔板持液量，分别为下塔、上塔温度分布位置变化速度；②下一时刻波形位置以及下一时刻温度分布，通过式(10)-(13)得到 <mrow><msub> <mi>S</mi> <mi>r</mi></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo></mrow><mo>=</mo><msub> <mi>S</mi> <mi>r</mi></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><mfrac> <msub><mi>dS</mi><mi>r</mi> </msub> <mi>dt</mi></mfrac><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><mo>&times;</mo><mi>&delta;</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow> <mo>(</mo> <mn>10</mn> <mo>)</mo></mrow> </mrow> <mrow><msub> <mi>S</mi> <mi>s</mi></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo></mrow><mo>=</mo><msub> <mi>S</mi> <mi>s</mi></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><mfrac> <msub><mi>dS</mi><mi>s</mi> </msub> <mi>dt</mi></mfrac><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><mo>&times;</mo><mi>&delta;</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow> <mo>(</mo> <mn>11</mn> <mo>)</mo></mrow> </mrow> <mrow><msub> <mover><mi>T</mi><mo>^</mo> </mover> <mi>i</mi></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo></mrow><mo>=</mo><msub> <mi>T</mi> <mrow><mi>min</mi><mo>,</mo><mi>r</mi> </mrow></msub><mo>+</mo><mfrac> <mrow><msub> <mi>T</mi> <mrow><mi>max</mi><mo>,</mo><mi>r</mi> </mrow></msub><mo>-</mo><msub> <mi>T</mi> <mrow><mi>min</mi><mo>,</mo><mi>r</mi> </mrow></msub> </mrow> <msup><mrow> <mn>1</mn> <mo>+</mo> <mi>e</mi></mrow><mrow> <mo>-</mo> <msub><mi>k</mi><mi>r</mi> </msub> <mrow><mo>(</mo><mi>i</mi><mo>-</mo><msub> <mi>S</mi> <mi>r</mi></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo></mrow><mo>)</mo> </mrow></mrow> </msup></mfrac> </mrow>i＝1，2，……，f-1(12) <mrow><msub> <mover><mi>T</mi><mo>^</mo> </mover> <mi>i</mi></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo></mrow><mo>=</mo><msub> <mi>T</mi> <mrow><mi>min</mi><mo>,</mo><mi>s</mi> </mrow></msub><mo>+</mo><mfrac> <mrow><msub> <mi>T</mi> <mrow><mi>max</mi><mo>,</mo><mi>s</mi> </mrow></msub><mo>-</mo><msub> <mi>T</mi> <mrow><mi>min</mi><mo>,</mo><mi>s</mi> </mrow></msub> </mrow> <msup><mrow> <mn>1</mn> <mo>+</mo> <mi>e</mi></mrow><mrow> <mo>-</mo> <msub><mi>k</mi><mi>s</mi> </msub> <mrow><mo>(</mo><mi>i</mi><mo>-</mo><msub> <mi>S</mi> <mi>s</mi></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo></mrow><mo>)</mo> </mrow></mrow> </msup></mfrac> </mrow>i＝f，f+1，……，n(13)其中，δ为采样周期，Sr(k)、Ss(k)分别为k采样时刻下塔、上塔温度分布曲线位置，Sr(k+1)、Ss(k+1)分别为k+1采样时刻下塔、上塔温度分布曲线位置，为k+1时刻下塔或上塔第i块塔板温度的预测值；温度预测误差补偿模块，用以从历史数据库中获取当前温度测量值，以及预测值，计算预测误差，修正温度预测值，通过式(14)(15)得到 <mrow><msub> <mi>r</mi> <mi>i</mi></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><mo>=</mo><msub> <mi>T</mi> <mi>i</mi></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><mo>-</mo><msub> <mover><mi>T</mi><mo>^</mo> </mover> <mi>i</mi></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow> </mrow>i＝1，2，……，n(14) <mrow><msubsup> <mi>X</mi> <mi>i</mi> <mo>%</mo></msubsup><mrow> <mo>(</mo> <mi>k</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo></mrow><mo>=</mo><msub> <mover><mi>T</mi><mo>^</mo> </mover> <mi>i</mi></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo></mrow><mo>+</mo><msub> <mi>r</mi> <mi>i</mi></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow> </mrow>i＝1，2，……，n(15)其中，Ti(k)为k+1采样时刻第i块塔板的温度，ri(k)为k采样时刻第i块塔板的温度预测误差，为k+1采样时刻第i块塔板的温度预测修正值。FDA0000022887840000013.tif,FDA0000022887840000023.tif,FDA0000022887840000028.tif,FDA0000022887840000033.tif
2.如权利要求1中所述的空分节能过程温度分布的非线性观测系统，其特征在于所 述的现场智能仪表与DCS系统中的存储装置、上位机通过现场总线相互连接。
3 如权利要求1或2所述的空分节能过程温度分布的非线性观测系统，其特征在于 所述的上位机还包括人机界面模块，用以设定采样周期，设定预测时间长度，显示观测器输 出的将来时刻温度预测值和当前时刻温度测量值。
4.一种如权利要求1所述的空分节能过程温度分布的非线性观测系统实现的非线性 观测方法，其特征在于所述的观测方法包括以下步骤1)确定采样周期，并将氮和氧相对于氩的相对挥发度αΝ、α<)，安东尼常数aN、bN、cN、a()、 h、q保存在历史数据库当中；2)从现场智能仪表获取k采样时刻压强P，各塔板温度Ti,计算第i块塔板处氮、氧的 液相组分浓度值，通过式(3) _(4)得到 其中，Xi, N为空分塔第i块塔板处氮的液相组分浓度，Xu为空分塔中第i块塔板处氧 的液相组分浓度；检测k采样时刻下塔塔顶汽相流量和上塔塔底液相流量值，同上塔压强数据，各塔板 温度数据，组分浓度的测量值一起存储到观测器系统的历史数据库当中；3)从历史数据库当中获取k采样时刻温度的测量值，拟合k采样时刻的温度分布函数 如式⑴⑵ 其中，1为下塔塔顶编号，f为进料板编号，η为上塔塔底编号， 为第i块塔板处预估温 度，Tfflinjr, Tmax,ρ kr, Tfflinjs, Tfflaxjs, ks为拟合参数，民、Ss分别为空分塔下塔、上塔液相温度分 布曲线的位置；4)计算k采样时刻的下塔上塔温度分布位置变化速度，通过式(5)-(9)得到 其中，k为当前采样时刻，下角标i为塔板编号，1为下塔塔顶编号，η为上塔塔底的编 号，下角标N、0分别代表氮和氧，上角标f代表进料，Xi, N、Xi, 0分别为空分塔中第i块塔板 处氮、氧的液相组分浓度，P为空分塔压强，Ti为第i块塔板处的温度，α N、C^分别为氮和 氧相对于氩的相对挥发度，aN、bN、cN、a。、b0, C0为安东尼常数，Fi (k)为k采样时刻第i块塔 板进料流量，Ln为上塔塔底液相流率，V1为下塔塔顶汽相流率，xn,N(k)、χη,ο (k)分别为k采 样时刻上塔塔底氮、氧的液相组分浓度，X1,N(k)、Xl, Jk)分别为k采样时刻下塔塔顶氮、氧的 液相组分浓度，yf,0(k)分别为k采样时刻进料板氧汽相组分浓度，yi,N(k).yi,0(k)分别为k 采样时刻下塔塔顶氮、氧汽相组分浓度，Xfi,N(k)，/"(ΙΟ分别为第i块塔板的进料液氮组分浓度和进料液氧组分浓度，M为塔板持液量，，分别为下塔上塔温度分布位at at置变化速度；5)下一时刻波形位置以及下一时刻温度分布，通过式(10)-(13)得到 其中，S为采样周期，民(k)、Ss(k)分别为k采样时刻下塔，上塔温度分布曲线位置， Sr(k+1)、Ss(k+1)分别为k+1采样时刻下塔，上塔温度分布曲线位置，f@ + l)为k+1时刻下 塔或上塔的第i块塔板温度的预测值；6)从历史数据库中获取当前温度测量值，以及预测值，计算预测误差，修正温度预测 值，通过式(14) (15)得到i = 1,2,……，n (14) ^ + = ^ + + = 1,2, ......，n (15)其中，，t (k)为k+1采样时刻第i块塔板的温度，r, (k)为k采样时刻第i块塔板的温 度预测误差，辩A + 1)为k+1采样时刻的温度预测修正值。
5.如权利要求4所述的非线性观测方法，其特征在于所述历史数据库为DCS系统的 存储装置，控制站读取历史数据库，显示观测器工作过程状态。
全文摘要
一种空分节能过程温度分布的非线性观测系统，包括与空分塔直接连接的现场智能仪表和DCS系统，所述DCS系统中包括存储装置、控制站和上位机，所述现场智能仪表与存储装置、控制站和上位机相连，所述的上位机包括用以预测空分塔将来时刻的温度分布的观测器，所述的观测器包括检测模块、温度分布函数拟合模块、将来时刻空分塔温度预测模块和温度预测误差补偿模块。本发明也提供了一种空分节能过程温度分布的非线性动态观测方法。本发明能够有效快速的在线预测空分节能过程的温度状态变化，准确地反应空分节能过程的强非线性性，耦合性等复杂非线性动态特性。
文档编号B01D3/14GK101884848SQ20101021388
公开日2010年11月17日 申请日期2010年6月30日 优先权日2010年6月30日
发明者刘兴高, 周叶翔 申请人:浙江大学