专利名称：一种高纯精馏过程的浓度非线性观测系统及方法
技术领域：
本发明涉及精馏过程的动态非线性建模领域，特别提出高纯精馏过程的浓度非线 性观测器系统的设计和实现方法。
背景技术：
精馏过程是化工过程中的一种核心过程，精馏塔是其中的一个关键单元。长久以 来，精馏塔因为高耗能，低能效的问题成为国际精馏领域研究的焦点。精馏塔具有十分复杂 的强非线性特征。导致高纯精馏过程的简化非线性建模，动态过程观测器设计成为一种世 界性难题。传统的复杂机理模型尽管能够精确的观测精馏塔的各种状态变量，但复杂的机理 模型结构导致模型在线求解效率太低，不能直接应用于控制器设计当中。而基于阶跃，脉冲 响应辨识的近似线性模型尽管具有较快的求解速度但是状态变量观测精度太低，甚至显著 偏离真实值。实现高纯精馏过程的高效率高精确度的非线性观测器对于该塔的控制设计节 能研究等方面起着至关重要的作用。
发明内容
为了克服现有的高纯精馏过程的浓度观测系统的状态变量观测精度太低、效率低 的不足，本发明提供一种高纯精馏过程的高效率，高精度的浓度观测的高纯精馏过程的浓 度非线性观测系统及方法。本发明解决其技术问题所采用的技术方案是一种高纯精馏过程的浓度非线性观测系统，包括与高纯精馏过程直接连接的现场 智能仪表和DCS系统，所述DCS系统包括存储装置。控制站和上位机，所述现场智能仪表与 存储装置、控制站和上位机相连，所述的上位机包括用以预测高纯精馏过程将来时刻的组 分浓度分布的观测器，所述观测器包括组分推断模块，用以从现场智能仪表获取温度，压强数据，计算高纯精馏过程各块 塔板的组分浓度，并将组分浓度计算结果存储在历史数据库当中，采用式(1)得到
其中，下脚标i为塔板编号，1为塔顶编号，f为进料板编号，η为塔底编号，Xi为第 i块塔板的液相组分浓度，P分别为精馏塔压强，Ti为第i块塔板的的温度，α为相对挥发 度，a、b、c为安东尼常数；组分分布函数拟合模块，用以采用历史数据库中组分推断模块计 算出的组分浓度数据，拟合液相组分浓度分布函数，并将拟合参数存储到历史数据库当中， 采用式⑵⑶得到
其中，友为第i块塔板处液相组分浓度预估浓度，XmiM、Xmau、kr、Xmin,s、Xmax,s、ks为 拟合参数，sr, Ss分别为高纯精馏过程精馏段、提馏段液相组分浓度分布曲线的位置；将来时刻高纯精馏过程状态预测模块，用以根据当前组分浓度数据以及组分浓度 分布函数预测将来时刻的液相组分浓度分布，该模块包括以下①从历史数据库中读取组分浓度测量值，计算当前时刻液相轻组分浓度分布曲线 位置变化速度，并将曲线位置变化速度存储到历史数据库当中，采用式(4)-(6)得到 其中，k为当前采样时刻，Xi (k)、Yi (k)分别为k采样时刻第i块塔板液相轻组分 浓度和汽相轻组分浓度，q为进料热状况，F为进料流率，V、L分别为气、液相摩尔流率，H为 持液量，Zf为进料组分浓度，Xh (k)、Xn(k)分别为k采样时刻第f-Ι块塔板和塔底的液相
轻组分浓度，Y1GihYfGO分别为k采样时刻塔顶和进料板处的汽相轻组分浓度，幻、幻分别为精馏段和提馏段的轻组分浓度分布曲线位置变化速度；②从历史数据库中获取曲线位置变化速度数据，组分浓度分布函数拟合参数，采 样周期并计算下一时刻波形位置以及下一时刻组分浓度分布，采用式(7)_(10)得到 其中，δ为采样周期，Sr(k),Ss(k)分别为k采样时刻精馏段、提馏段液相轻组分 浓度分布位置，sr(k+l)、Ss(k+1)分别为k+1采样时刻精馏段、提馏段液相轻组分浓度分布 位置，i,,0 + l)、i,JA: + l)分别为k+Ι时刻精馏段、提馏段第i块塔板液相轻组分浓度分布 的预测值；③从历史数据库中获取当前组分浓度测量值，以及预测值，修正组分浓度预测值， 采用式(11) (12)得到 其中，$(幻为k时刻精馏段或提馏段的第i块塔板液相轻组分浓度分布的预测值， ri(k)为k采样时刻第i块塔板的预测误差，！灿+1)为k+i采样时刻的液相轻组分浓度预 测修正值即液相轻组分浓度的观测值。所述的上位机还包括人机界面模块，用以设定采样周期，显示观测器输出的将来 时刻组分浓度预测值和当前时刻组分浓度测量值。一种高纯精馏过程的浓度非线性过程观测方法，所述的观测方法包括以下步骤1)确定采样周期，并将相对挥发度α，安东尼常数a、b、C、保存在历史数据库当 中；2)检测k采样时刻精馏塔压强P，以及各塔板温度Ti,计算液相轻组分浓度值，采 用式⑴得到
Tj+c a 其中，下脚标i为塔板编号，1为塔顶编号，f为进料板编号，η为塔底编号，Xi为第 i块塔板的液相组分浓度，P分别为精馏塔压强，Ti为第i块塔板的的温度，α为相对挥发 度，a、b、c为安东尼常数；3)从历史数据库当中获取k采样时刻组分浓度的测量值，拟合k采样时刻的组分 浓度分布函数如式(2) (3)Ii = Zmm, + :(》Γ i = 1'2'……，f"l (2)Ii = Xmm, +土 = f，f+1，……，η (3)其中，友为第i块塔板处液相组分浓度预估浓度，Xfflin,r>xfflax,r>krafflin,safflax,s>ks为 拟合参数，sr, Ss分别为高纯精馏过程精馏段、提馏段液相组分浓度分布曲线的位置；4)计算k采样时刻的组分浓度分布曲线位置变化速度，并将曲线位置变化速度存 储到历史数据库当中，采用式(4)_(6)得到
其中，k为当前采样时刻，Xi (k)、Yi (k)分别为k采样时刻第i块塔板液相轻组分 浓度和汽相轻组分浓度，q为进料热状况，F为进料流率，V、L分别为气、液相摩尔流率，H为 持液量，Zf为进料组分浓度，Xh (k)、Xn(k)分别为k采样时刻第f-Ι块塔板和塔底的液相
轻组分浓度，Y1GO、Yf (k)分别为k时刻塔顶和进料板处的汽相轻组分浓度，^(k)
at at
分别为精馏段和提馏段的轻组分浓度分布曲线位置变化速度； 5)预测k+Ι采样时刻的组分浓度分布，从历史数据库中获取曲线位置变化速度数 据，组分浓度分布函数拟合参数，采样周期值，采用式(7)_(10)得到
H ⑶、(k+ l) = x^s + 7^Jk +1)) 土 = f，f +1，……，η (10)其中，δ为采样周期，Sr(k),Ss(k)分别为k采样时刻精馏段、提馏段液相轻组分 浓度分布位置，sr(k+l)、Ss(k+1)分别为k+1采样时刻精馏段、提馏段液相轻组分浓度分布 位置，++ 分别为k+Ι采样时刻精馏段、提馏段段第i块塔板液相轻组分浓度 分布的预测值；6)从历史数据库中获取k采样时刻组分浓度测量值，以及预测值，计算k采样时刻 的组分浓度预测误差过程如式(11)ri{k) = Xi{k)-Xi{k)i = 1,2,……，η (11)其中， (幻为k时刻精馏段或提馏段的第i块塔板液相轻组分浓度分布的预测值， A (k)为k采样时刻第i块塔板的预测误差；7)修正组分浓度预测值，采用式(12)得到f^k + l) = X^k+ += 1,2, ......，η (12)其中，！ΚΑ: + !)为k+1采样时刻的液相轻组分浓度预测修正值即液相轻组分浓度的 观测值。进一步，所述历史数据库为DCS系统的存储装置，控制站读取历史数据库，显示观 测器工作过程状态。本发明的有益效果主要表现在1.基于高纯精馏过程的机理模型以及传热、传质 原理建立了高精度的高纯精馏过程的简化非线性动态模型，能够准确跟踪反应高纯精馏过 程的浓度非线性特征，状态预测比传统辨识模型更为精确；2.在确保浓度观测器跟踪高精 度的情况下，同时具有较快的在线求解速度，效率远远高于机理模型
图1是本发明所提出的高纯精馏过程非线性观测系统的结构图。图2是上位机观测器实现方法的原理图。图3是本发明所提出的高纯精馏过程非线性观测方法流程图。
具体实施例方式
下面根据附图具体说明本发明。实施例1参照图1、图2和图3，一种高纯精馏过程的浓度非线性观测系统，包括与高纯精馏 过程1直接连接的现场智能仪表2和DCS系统13，所述DCS系统中的存储装置4、控制站5和上位机6，其中所述现场智能仪表2与存储装置4、控制站5和上位机6通过现场总线和 数据接口 3依次相连；所述的智能仪表用以实现检测模块7的功能即检测得到高纯精馏过 程的温度、压强数据，所述的数据接口用以实现I/O模块8的功能即数据的输入输出所述的上位机6包括用以预测高纯精馏过程将来时刻的组分浓度分布的观测器， 所述的观测器包括组分推断模块9、组分分布函数拟合模块10和将来时刻高纯精馏过程状 态预测模块11。所述的组分推断模块9，上位机从现场智能仪表获取温度，压强数据，计算高纯精 馏过程各块塔板的组分浓度，并将组分浓度计算结果存储在历史数据库当中，采用式(1) 得到
其中，下脚标i为塔板编号，1为塔顶编号，f为进料板编号，η为塔底编号，Xi为第 i块塔板的液相组分浓度，P分别为精馏塔压强，Ti为第i块塔板的的温度，α为相对挥发 度，a、b、c为安东尼(Antonie)常数；所述的组分分布函数拟合模块10，采用历史数据库中组分推断模块计算出的组分 浓度数据，拟合液相组分浓度分布函数，并将拟合参数存储到历史数据库当中，采用式(2) ⑶得到 其中，友为第i块塔板处液相组分浓度预估浓度，Xfflin,r>xfflax,r>krafflin,safflax,s>ks为 拟合参数，sr, Ss分别为高纯精馏过程精馏段、提馏段液相组分浓度分布曲线的位置。所述的将来时刻高纯精馏过程状态预测模块11，根据当前组分浓度数据以及组分 分布函数预测将来时刻的液相组分浓度分布，该模块包括以下三部分①从历史数据库中读取组分浓度测量值，计算当前时刻组分浓度分布曲线位置变 化速度，并将曲线位置变化速度存储到历史数据库当中，采用式(4)_(6)得到
其中，k为当前采样时刻，Xi (k)、Y1 (k)分别为k采样时刻第i块塔板轻组分液相 轻组分浓度和汽相轻组分浓度，q为进料热状况，F为进料流率，V、L分别为气、液相摩尔流 率，H为持液量，Zf为进料组分浓度，Xh (k)、Xn(k)分别为k采样时刻第f-Ι块塔板和塔底
的液相轻组分浓度，Y1 (k)、Yf (k)分别为k时刻塔顶和进料板处的汽相轻组分浓度，、幻分别为精馏段和提馏段的轻组分浓度分布曲线位置变化速度；
at②从历史数据库中获取曲线位置变化速度数据，组分浓度分布函数拟合参数，采 样周期并计算下一时刻波形位置以及下一时刻组分浓度分布，采用式(7)_(10)得到 其中，δ为采样周期，Sr(k),Ss(k)分别为k采样时刻精馏段、提馏段液相轻组分 浓度分布位置，sr(k+l)、Ss(k+1)分别为k+1采样时刻精馏段、提馏段液相轻组分浓度分布 位置，i,,@ + l)、i,JA: + l)分别为k+Ι时刻精馏段、提馏段第i块塔板液相轻组分浓度分布 的预测值；③从历史数据库中获取当前组分浓度测量值，以及预测值，修正组分浓度预测值， 采用式(11) (12)得到 其中， (幻为k时刻精馏段或提馏段的第i块塔板液相轻组分浓度分布的预测值， ri(k)为k采样时刻第i块塔板的预测误差，！灿+1)为k+i采样时刻的液相轻组分浓度预 测修正值即液相轻组分浓度的观测值。所述的上位机还包括人机界面模块13，用于设定采样周期，显示观测器输出的将 来时刻组分浓度预测值和当前时刻组分浓度测量值.实施例2参照图1、图2和图3，一种高纯精馏过程的浓度非线性过程观测方法，如图3所 示，所述的观测方法包括以下步骤1)如图3模块14所示，确定采样周期，并将相对挥发度α，安东尼(Antonie)常 数a、b、C、保存在历史数据库当中；2)如图3模块15所示，检测k采样时刻精馏塔压强P，以及各塔板温度Ti,计算液 相轻组分浓度值，采用式(1)得到 其中，下脚标i为塔板编号，1为塔顶编号，f为进料板编号，η为塔底编号，Xi为第 i块塔板的液相组分浓度，P分别为精馏塔压强，Ti为第i块塔板的的温度，α为相对挥发 度，a、b、c为安东尼(Antonie)常数；；3)如图3模块16所示，从历史数据库当中获取k采样时刻组分浓度的测量值，拟合k采样时刻的组分浓度分布函数如式(2) (3) 其中，戈为第i块塔板处液相组分浓度预估浓度，Xfflin,r>xfflax,r>krafflin,safflax,s>ks为 拟合参数，sr, Ss分别为高纯精馏过程精馏段、提馏段液相组分浓度分布曲线的位置；4)如图3模块17所示，计算k采样时刻的组分浓度分布曲线位置变化速度，并将
曲线位置变化速度存储到历史数据库当中，采用式(4)-(6)得到 _0] + = ……，n ⑷
其中，k为当前采样时刻，Xi (k) ,Yi (k)分别为k采样时刻第i块塔板轻组分液相轻 组分浓度和汽相轻组分浓度，q为进料热状况，F为进料流率，V、L分别为气、液相摩尔流率， H为持液量，Zf为进料组分浓度，Xf_i (k)、Xn(k)分别为k采样时刻第f-Ι块塔板和塔底的液
相轻组分浓度,Y1GO、Yf (k)分别为k采样时刻塔顶和进料板处的汽相轻组分浓度，、
at
幻分别为精馏段和提馏段的轻组分浓度分布曲线位置变化速度；
at5)如图3模块18所示，预测k+1采样时刻的组分浓度分布。从历史数据库中获取 曲线位置变化速度数据，组分浓度分布函数拟合参数，采样周期值，采用式(7)-(10)得到 其中，δ为采样周期，Sr(k),Ss(k)分别为k采样时刻精馏段、提馏段液相轻组分 浓度分布位置，Sr (k+1)、Ss(k+1)分别为k+1采样时刻精馏段、提馏段液相轻组分浓度分布 位置，+ 、友+ 分别为k+Ι采样时刻精馏段、提馏段段第i块塔板液相轻组分浓 度分布的预测值；6)如图3模块19所示，从历史数据库中获取k采样时刻组分浓度测量值，以及预 测值，计算k采样时刻的组分浓度预测误差过程如式(11)ri(k) = Xi(k)-Xi(k)i = 1,2,……，η (11)
11
其中， (幻为k采样时刻精馏段或提馏段的第i块塔板液相轻组分浓度分布的预 测值，r, (k)为k采样时刻第i块塔板的预测误差7)如图3模块20所示，修正组分浓度预测值，采用式(12)得到f^k + \) = XXk + l) + rXk)i = 1,2, .........，η (12)其中；^^ + 1)为k+Ι采样时刻的液相轻组分浓度预测修正值即液相轻组分浓度的 观测值。所述的历史数据库为DCS系统中的存储装置4，所述的DCS系统包括数据接口 3， 存储装置4和控制站5，其中控制站可以读取历史数据库，显示观测器工作过程状态。上述实施例用来解释说明本发明，而不是对本发明进行限制，在本发明的精神和 权利要求的保护范围内，对本发明作出的任何修改和改变，都落入本发明的保护范围。
权利要求
一种高纯精馏过程的浓度非线性观测系统，包括与高纯精馏过程直接连接的现场智能仪表和DCS系统，所述DCS系统包括存储装置。控制站和上位机，所述现场智能仪表与存储装置、控制站和上位机相连，其特征在于所述的上位机包括用以预测高纯精馏过程将来时刻的组分浓度分布的观测器，所述观测器包括组分推断模块，用以从现场智能仪表获取温度，压强数据，计算高纯精馏过程各块塔板的组分浓度，并将组分浓度计算结果存储在历史数据库当中，采用式(1)得到 <mrow><msub> <mi>X</mi> <mi>i</mi></msub><mo>=</mo><mfrac> <mrow><mi>P</mi><mo>×</mo><mi>α</mi><mo>×</mo><msup> <mn>10</mn> <mrow><mfrac> <mrow><msub> <mi>T</mi> <mi>i</mi></msub><mo>+</mo><mi>c</mi> </mrow> <mi>b</mi></mfrac><mo>-</mo><mi>a</mi> </mrow></msup><mo>-</mo><mn>1</mn> </mrow> <mrow><mi>α</mi><mo>-</mo><mn>1</mn> </mrow></mfrac> </mrow>i＝1，2，……，n(1)其中，下脚标i为塔板编号，1为塔顶编号，f为进料板编号，n为塔底编号，Xi为第i块塔板的液相组分浓度，P分别为精馏塔压强，Ti为第i块塔板的的温度，α为相对挥发度，a、b、c为安东尼常数；组分分布函数拟合模块，用以采用历史数据库中组分推断模块计算出的组分浓度数据，拟合液相组分浓度分布函数，并将拟合参数存储到历史数据库当中，采用式(2)(3)得到 <mrow><msub> <mover><mi>X</mi><mo>^</mo> </mover> <mi>i</mi></msub><mo>=</mo><msub> <mi>X</mi> <mrow><mi>min</mi><mo>,</mo><mi>r</mi> </mrow></msub><mo>+</mo><mfrac> <mrow><msub> <mi>X</mi> <mrow><mi>max</mi><mo>,</mo><mi>r</mi> </mrow></msub><mo>-</mo><msub> <mi>X</mi> <mrow><mi>min</mi><mo>,</mo><mi>r</mi> </mrow></msub> </mrow> <msup><mrow> <mi>l</mi> <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(2) <mrow><msub> <mover><mi>X</mi><mo>^</mo> </mover> <mi>i</mi></msub><mo>=</mo><msub> <mi>X</mi> <mrow><mi>min</mi><mo>,</mo><mi>s</mi> </mrow></msub><mo>+</mo><mfrac> <mrow><msub> <mi>X</mi> <mrow><mi>max</mi><mo>,</mo><mi>s</mi> </mrow></msub><mo>-</mo><msub> <mi>X</mi> <mrow><mi>min</mi><mo>,</mo><mi>s</mi> </mrow></msub> </mrow> <mrow><mn>1</mn><msup> <mrow><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> </mrow></mfrac> </mrow>i＝f，f+1，……，n(3)其中，为第i块塔板处液相组分浓度预估浓度，Xmin，r、Xmax，r、kr、Xmin，s、Xmax，s、ks为拟合参数，Sr、Ss分别为高纯精馏过程精馏段、提馏段液相组分浓度分布曲线的位置；将来时刻高纯精馏过程状态预测模块，用以根据当前组分浓度数据以及组分浓度分布函数预测将来时刻的液相组分浓度分布，该模块包括以下①从历史数据库中读取组分浓度测量值，计算当前时刻液相轻组分浓度分布曲线位置变化速度，并将曲线位置变化速度存储到历史数据库当中，采用式(4)-(6)得到 <mrow><msub> <mi>Y</mi> <mi>i</mi></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><mo>=</mo><mfrac> <mrow><mi>α</mi><msub> <mi>X</mi> <mi>i</mi></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow> </mrow> <mrow><mrow> <mo>(</mo> <mi>α</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo></mrow><msub> <mi>X</mi> <mi>i</mi></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><mo>+</mo><mn>1</mn> </mrow></mfrac> </mrow>i＝1，2，……，n(4) <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><msub> <mi>VY</mi> <mi>f</mi></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><mo>-</mo><msub> <mi>LX</mi> <mrow><mi>f</mi><mo>-</mo><mn>1</mn> </mrow></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><mo>-</mo><msub> <mi>VY</mi> <mn>1</mn></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow> </mrow> <mrow><msub> <mi>HX</mi> <mrow><mi>f</mi><mo>-</mo><mn>1</mn> </mrow></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow> </mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow> <mo>(</mo> <mn>5</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><msub> <mi>VY</mi> <mi>f</mi></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><mo>-</mo><msub> <mi>LX</mi> <mi>n</mi></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><mo>+</mo><msub> <mi>LX</mi> <mrow><mi>f</mi><mo>-</mo><mn>1</mn> </mrow></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><mo>+</mo><msub> <mi>FZ</mi> <mi>f</mi></msub> </mrow> <mrow><mi>H</mi><mrow> <mo>(</mo> <msub><mi>X</mi><mrow> <mi>f</mi> <mo>-</mo> <mn>1</mn></mrow> </msub> <mrow><mo>(</mo><mi>k</mi><mo>)</mo> </mrow> <mo>-</mo> <msub><mi>X</mi><mi>n</mi> </msub> <mrow><mo>(</mo><mi>k</mi><mo>)</mo> </mrow> <mo>)</mo></mrow> </mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo></mrow> </mrow>其中，k为当前采样时刻，Xi(k)、Yi(k)分别为k采样时刻第i块塔板液相轻组分浓度和汽相轻组分浓度，q为进料热状况，F为进料流率，V、L分别为气、液相摩尔流率，H为持液量，Zf为进料组分浓度，Xf-1(k)、Xn(k)分别为k采样时刻第f-1块塔板和塔底的液相轻组分浓度，Y1(k)、Yf(k)分别为k时刻塔顶和进料板处的汽相轻组分浓度，分别为精馏段和提馏段的轻组分浓度分布曲线位置变化速度；②从历史数据库中获取曲线位置变化速度数据，组分浓度分布函数拟合参数，采样周期并计算下一时刻波形位置以及下一时刻组分浓度分布，采用式(7)-(10)得到 <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>×</mo><mi>δ</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow> <mo>(</mo> <mn>7</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>×</mo><mi>δ</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow> <mo>(</mo> <mn>8</mn> <mo>)</mo></mrow> </mrow> <mrow><msub> <mover><mi>X</mi><mo>^</mo> </mover> <mrow><mi>i</mi><mo>,</mo><mi>r</mi> </mrow></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo></mrow><mo>=</mo><msub> <mi>X</mi> <mrow><mi>min</mi><mo>,</mo><mi>r</mi> </mrow></msub><mo>+</mo><mfrac> <mrow><msub> <mi>X</mi> <mrow><mi>max</mi><mo>,</mo><mi>r</mi> </mrow></msub><mo>-</mo><msub> <mi>X</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(9) <mrow><msub> <mover><mi>X</mi><mo>^</mo> </mover> <mrow><mi>i</mi><mo>,</mo><mi>s</mi> </mrow></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo></mrow><mo>=</mo><msub> <mi>X</mi> <mrow><mi>min</mi><mo>,</mo><mi>s</mi> </mrow></msub><mo>+</mo><mfrac> <mrow><msub> <mi>X</mi> <mrow><mi>max</mi><mo>,</mo><mi>s</mi> </mrow></msub><mo>-</mo><msub> <mi>X</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(10)其中，δ为采样周期，Sr(k)、Ss(k)分别为k采样时刻精馏段、提馏段液相轻组分浓度分布位置，Sr(k+1)、Ss(k+1)分别为k+1采样时刻精馏段、提馏段液相轻组分浓度分布位置，分别为k+1采样时刻精馏段、提馏段第i块塔板液相轻组分浓度分布的预测值；③从历史数据库中获取当前组分浓度测量值，以及预测值，修正组分浓度预测值，采用式(11)(12)得到 <mrow><msub> <mi>r</mi> <mi>i</mi></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><mo>=</mo><msub> <mi>X</mi> <mi>i</mi></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><mo>-</mo><msub> <mover><mi>X</mi><mo>^</mo> </mover> <mi>i</mi></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow> </mrow>i＝1，2，……，n(11) <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>X</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(12)其中，为k采样时刻精馏段或提馏段的第i块塔板液相轻组分浓度分布的预测值，ri(k)为k采样时刻第i块塔板的预测误差，为k+1采样时刻的液相轻组分浓度预测修正值即液相轻组分浓度的观测值。FDA0000022884930000014.tif,FDA0000022884930000023.tif,FDA0000022884930000024.tif,FDA0000022884930000029.tif,FDA00000228849300000210.tif,FDA00000228849300000213.tif,FDA00000228849300000214.tif
2.如权利要求1所述的高纯精馏过程的浓度非线性观测器系统，其特征在于所述的 上位机还包括人机界面模块，用以设定采样周期，显示观测器输出的将来时刻组分浓度预 测值和当前时刻组分浓度测量值。
3.—种如权利要求1所述的高纯精馏过程的浓度非线性过程观测系统实现的浓度非 线性过程观测方法，其特征在于所述的观测方法包括以下步骤1)确定采样周期，并将相对挥发度a，安东尼常数a、b、c、保存在历史数据库当中；2)检测k采样时刻精馏塔压强P，以及各塔板温度！\，计算液相轻组分浓度值，采用式 ⑴得到 其中，下脚标i为塔板编号，1为塔顶编号，f 为进料板编号，n为塔底编号，&为第i块 塔板的液相组分浓度，P分别为精馏塔压强，T,为第i块塔板的的温度，a为相对挥发度， a、b、c为安东尼常数；3)从历史数据库当中获取k采样时刻组分浓度的测量值，拟合k采 样时刻的组分浓度分布函数如式(2) (3)友=(m’々::ixr i = i，2，……，f—1 ⑵ 友i = f，f+1，……，n ⑶其中，友为第i块塔板处液相组分浓度预估浓度，xfflin,rafflax,r>kr.xfflin,safflax,s>ks为拟合 参数，sr,ss分别为高纯精馏过程精馏段、提馏段液相组分浓度分布曲线的位置；4)计算k采样时刻的组分浓度分布曲线位置变化速度，并将曲线位置变化速度存储到历史数据库当中，采用式(4)-(6)得到 其中，k为当前采样时刻，X, (k)、t (k)分别为k采样时刻第i块塔板液相轻组分浓度 和汽相轻组分浓度，q为进料热状况，F为进料流率，V、L分别为气、液相摩尔流率，H为持液 量，Zf为进料组分浓度，Xf_i (k)、Xn(k)分别为k采样时刻第f_l块塔板和塔底的液相轻组分浓度，Yi(k)、Yf (k)分别为k时刻塔顶和进料板处的汽相轻组分浓度，、^分别at at为精馏段和提馏段的轻组分浓度分布曲线位置变化速度；5)预测k+1采样时刻的组分浓度分布，从历史数据库中获取曲线位置变化速度数据， 组分浓度分布函数拟合参数，采样周期值，采用式(7)_(10)得到 其中，S为采样周期，民(k)、Ss(k)分别为k采样时刻精馏段、提馏段液相轻组分浓度 分布位置，S, (k+1)、Ss(k+l)分别为k+1采样时刻精馏段、提馏段液相轻组分浓度分布位置， Xhr(k + l) + 分别为k+1采样时刻精馏段、提馏段段第i块塔板液相轻组分浓度分布 的预测值；6)从历史数据库中获取k采样时刻组分浓度测量值，以及预测值，计算k采样时刻的组 分浓度预测误差过程如式(11)ri(k) = Xi(k)-Xi(k) i = 1,2,……，n (11)其中，f (幻为k时刻精馏段或提馏段的第i块塔板液相轻组分浓度分布的预测值，^ (k) 为k采样时刻第i块塔板的预测误差；7)修正组分浓度预测值，采用式(12)得到+ = XXk + Vf + rXk)! = 1,2, ......, n (12)其中，为k+i采样时刻的液相轻组分浓度预测修正值即液相轻组分浓度的观测值。
4.如权利要求3所述的浓度非线性过程观测方法，其特征在于所述历史数据库为DCS 系统的存储装置，控制站读取历史数据库，显示观测器工作过程状态。
全文摘要
一种高纯精馏过程的浓度非线性观测系统，包括与高纯精馏过程直接连接的现场智能仪表和DCS系统，所述DCS系统包括存储装置、控制站和上位机，所述现场智能仪表与存储装置、控制站和上位机相连，所述的上位机包括用以预测高纯精馏过程将来时刻的组分浓度分布的观测器，所述的观测器包括组分推断模块、组分分布函数拟合模块和将来时刻高纯精馏过程浓度预测模块。本发明也提供了一种高纯精馏过程的浓度非线性观测方法。本发明能够有效快速的在线预测高纯精馏过程将来时刻的浓度变化，准确地反应高纯精馏过程的强非线性动态特性。
文档编号B01D3/42GK101884849SQ20101021398
公开日2010年11月17日 申请日期2010年6月30日 优先权日2010年6月30日
发明者刘兴高, 周叶翔 申请人:浙江大学