专利名称：：一种高效节能精馏过程的高纯控制系统及方法
技术领域：
：本发明涉及精馏控制领域，尤其提出了一种高效节能精馏过程的高纯控制系统及方法。
背景技术：
：精馏过程是化工过程中的一种核心过程，精馏塔是其中的一个关键单元。长久以来，精馏塔因为高耗能，低能效的问题成为国际精馏领域研究的焦点。目前针对精馏过程的能耗问题主要有两方面的解决方案一种设计新型结构，利用热量耦合实现能量重复利用达到节能目的，一种设计高效精馏过程控制策略，提高产品生产质量减少废料从而达到节能目的。尽管有较多实验研究证明高效节能精馏塔能够显著提高能源利用率，但是由于高效节能精馏塔的精馏段与提馏段之间存在极强的耦合性且该塔具有十分复杂的强非线性，该塔的控制策略设计显得尤为困难。传统的PID、内模控制方案等已经不能满足高效节能精馏过程的控制品质要求，尤其在高纯精馏过程控制当中，这些方案已经很难使精馏过程稳定。而基于线性辨识模型的预测控制方案只能工作在稳态工作点附近，稍微增大干扰幅度，或者设定值阶跃变化系统控制质量则出现明显下降。事实表明设计有效的高效精馏过程的有效高纯控制方案是高效精馏过程平稳操作的关键环节。
发明内容为了克服目前高效精馏过程的控制效果不好、控制品质不理想的不足，本发明提供一种控制效果良好、控制品质理想的高效节能精馏过程的高纯控制系统及方法。本发明解决其技术问题所采用的技术方案是—种高效节能精馏过程的高纯控制系统，包括与高效节能精馏塔直接连接的现场智能仪表和DCS系统，所述DCS系统包括存储装置、控制站和上位机，现场智能仪表与存储装置、控制站和上位机连接，所述的上位机包括用以计算输出高效节能精馏塔控制变量值的高纯控制器，所述高纯控制器包括组分推断模块，用以从智能仪表获取温度，压强数据，计算高效节能精馏塔各块塔板的组分浓度，并将组分浓度计算结果存储在历史数据库当中，计算式为(1)(2):<formula>formulaseeoriginaldocumentpage5</formula><formula>formulaseeoriginaldocumentpage5</formula>其中k为当前采样时刻，下脚标i为塔板编号，l为塔顶编号，f为进料板编号，n为塔底编号，Xi(k)为k采样时刻液相轻组分浓度，Pr(k)为k采样时刻精馏段压强、Ps提馏段压强，Ti(k)为k采样时刻塔内各块塔板的温度，a为相对挥发度，a、b、c为安东尼常数。模型参数自适应校正拟合模块，用以采用历史数据库中组分推断模块计算出的组分浓度数据，在线拟合模型函数，并将拟合参数存储到历史数据库当中，拟合函数如式(3)(4):<formula>formulaseeoriginaldocumentpage6</formula><formula>formulaseeoriginaldocumentpage6</formula>其中义,为第i块塔板处液相组分浓度预估浓度，Xmin,r，XMX,r，kr，Xmin,s，Xmax,s，ks为拟合参数，Sr，Ss分别为高效节能精馏塔精馏段、提馏段液相组分浓度分布的位置；高效节能精馏过程的高纯控制律求解模块，用以根据当前组分浓度数据，模型函数和当前时刻操作变量值求取当前的控制变量的理想改变值，求解控制律代数方程组如式(5)至式(12)柳=,=、，f=",……,"(5)<formula>formulaseeoriginaldocumentpage6</formula>其中k为当前采样时刻，T为采样周期，&(k)、1(k)分别为k采样时刻第i块塔板轻组分液相轻组分浓度和汽相轻组分浓度，Qi(k)为第i块塔板之间的热耦合量，UA为传热速率，Xi+f—Jk)为k采样时刻第i+f-1块塔板液相轻组分浓度，q(k)为k采样时刻进料热状况，Pr(k)为当k采样时刻精馏段压强，F为进料流率，Zf为进料组分浓度，K(k)，Vf(k)分别为k采样时刻塔顶和进料板处的汽相流率，Lf—工(k)，Ln(k)分别为k采样时刻第f-1块塔板和塔底的液相流率，H为持液量，Xf—工(k)，Xn(k)分别为k采样时刻第f-1块塔板和塔底的液相轻组分浓度，l(k)，Yf(k)分别为k采样时刻塔顶和进料板处的汽相轻组分浓度，Kl，K2，K3，K4为控制律参数，X二Xn*分别为塔顶塔底液相轻组分浓度设定值，&(k)，Xn(k)为k时刻塔顶塔底的液相轻组分浓度值Aq(k)，APjk)，分别为当前时刻高效节能精馏塔控制变量即进料热状况和精馏段压强的当前理想改变值。所述的上位机还包括人机界面模块，用以设定采样周期T，控制律参数K2，K3，K4和塔顶塔底液相轻组分浓度设定值&*，Xn*，并显示控制器的输出曲线和被控变量即塔顶塔底液相轻组分浓度的记录曲线。—种高效节能精馏过程的高纯控制方法，所述的控制方法包括以下步骤1)确定采样周期T，并将T值，相对挥发度a，提馏段压强Ps，安东尼常数a、b、c、保存在历史数据库当中；2)设定控制律参数&，K2，K3，K4和塔顶塔底液相轻组分浓度设定值X,，Xn*3)从智能仪表获取k采样时刻精馏段压强&提馏段压强Ps，以及各塔板温度I,计算液相轻组分浓度值，计算式为(1)(2):，+糊=，)=—1a-1-1a-1(1)(2)其中k为当前采样时刻，下脚标i为塔板编号，1为塔顶编号，f为进料板编号，n为塔底编号，Xi(k)为k采样时刻液相轻组分浓度，Pr(k)为k采样时刻精馏段压强、Ps提馏段压强，1\(k)为k采样时刻塔内各块塔板的温度，a为相对挥发度，a、b、c为安东尼常数；4)用历史数据库中组分推断模块计算出的组分浓度数据，在线拟合模型函数，并将拟合参数存储到历史数据库当中，拟合函数如式(3)(4):+<formula>formulaseeoriginaldocumentpage7</formula>其中l,为第i块塔板处液相组分浓度预估浓度，Ss为拟合参数；5)根据当前组分浓度数据，模型函数和当前时刻操作变量值求取当前的控制变量的理想改变值，求解控制律代数方程组如式(5)至式(12)<formula>formulaseeoriginaldocumentpage7</formula><formula>formulaseeoriginaldocumentpage7</formula><formula>formulaseeoriginaldocumentpage8</formula>其中k为当前采样时刻，T为采样周期，&(k)、1(k)分别为k采样时刻第i块塔板轻组分液相轻组分浓度和汽相轻组分浓度，Qi(k)为第i块塔板之间的热耦合量，UA为传热速率，Xi+f—Jk)为k采样时刻第i+f-1块塔板液相轻组分浓度，q(k)为k采样时刻进料热状况，Pr(k)为当k采样时刻精馏段压强，F为进料流率，Zf为进料组分浓度，K(k)，Vf(k)分别为k采样时刻塔顶和进料板处的汽相流率，Lf—工(k)，Ln(k)分别为k采样时刻第f-1块塔板和塔底的液相流率，H为持液量，Xf—工(k)，Xn(k)分别为k采样时刻第f-1块塔板和塔底的液相轻组分浓度，l(k)，Yf(k)分别为k采样时刻塔顶和进料板处的汽相轻组分浓度，Kl，K2，K3，K4为控制律参数，X二Xn*分别为塔顶塔底液相轻组分浓度设定值，&(k)，Xn(k)为k时刻塔顶塔底的液相轻组分浓度值Aq(k)，APjk)，分别为当前时刻高效节能精馏塔控制变量即进料热状况和精馏段压强的当前理想改变值；6)将当前时刻高效节能精馏塔控制变量即进料热状况和精馏段压强的当前理想改变值Aq(k)，APr(k)输送给DCS系统中的控制站，调整高效节能精馏塔的进料热状况值和精馏段压强值。所述的历史数据库为DCS系统中存储装置，控制站读取历史数据库，显示高效节能精馏塔工作过程状态。本发明的有益效果主要表现在1.高纯控制方案建立在高精度非线性动态模型基础上，能够及时抑制干扰作用；2.控制方案较好地处理了耦合问题，能够快速准确地跟踪设定值变化；3.控制效果好、控制品质高。图1是本发明所提出的高效节能精馏过程的高纯控制系统的结构图。图2是上位机控制器实现方法的功能模块结构图。具体实施例方式下面根据附图具体说明本发明。实施例1参照图1，图2，一种高效节能精馏过程的高纯控制系统，包括与高效节能精馏塔1直接连接的现场智能仪表2和DCS系统，所述DCS系统包括存储装置4，控制站5及上位机6，其中现场智能仪表2与数据接口3连接，所述数据接口3与现场总线连接，所述现场总线与存储装置4、控制站5和上位机6连接；所述的上位机6包括用以计算输出高效节能精馏塔控制变量值的高纯控制器，所述高纯控制器包括组分推断模块9，模型参数自适应校正拟合模块IO，高效节能精馏过程的高纯控制律求解模块11;组分推断模块9，用以上位机从智能仪表2获取温度，压强数据，计算高效节能精馏塔各块塔板的组分浓度，并将组分浓度计算结果存储在历史数据库当中，计算式为(1)(2):<formula>formulaseeoriginaldocumentpage9</formula>其中k为当前采样时刻，下脚标i为塔板编号，1为塔顶编号，f为进料板编号，n为塔底编号，Xi(k)为k采样时刻液相轻组分浓度，Pr(k)为k采样时刻精馏段压强、Ps提馏段压强，Ti(k)为k采样时刻塔内各块塔板的温度，a为相对挥发度，a、b、c为安东尼常数。所述的模型参数自适应校正拟合模块IO采用历史数据库中组分推断模块计算出的组分浓度数据，在线拟合模型函数，并将拟合参数存储到历史数据库当中，拟合函数如式(3)(4):<formula>formulaseeoriginaldocumentpage9</formula><formula>formulaseeoriginaldocumentpage9</formula>其中i,为第i块塔板处液相组分浓度预估浓度，xmin,r，xMX,r，kr，xmin,s，xmax,s，ks为拟合参数，Sr，ss分别为高效节能精馏塔精馏段、提馏段液相组分浓度分布的位置。所述的高效节能精馏过程的高纯控制律求解模块11根据当前组分浓度数据，模型函数和当前时刻操作变量值求取当前的控制变量的理想改变值，求解控制律代数方程组如式(5)至式(12)<formula>formulaseeoriginaldocumentpage9</formula><formula>formulaseeoriginaldocumentpage9</formula>其中k为当前采样时刻，T为采样周期，&(k)、1(k)分别为k采样时刻第i块塔板轻组分液相轻组分浓度和汽相轻组分浓度，Qi(k)为第i块塔板之间的热耦合量，UA为传热速率，Xi+f—Jk)为k采样时刻第i+f-1块塔板液相轻组分浓度，q(k)为k采样时刻进料热状况，Pr(k)为当k采样时刻精馏段压强，F为进料流率，Zf为进料组分浓度，K(k)，Vf(k)分别为k采样时刻塔顶和进料板处的汽相流率，Lf—工(k)，Ln(k)分别为k采样时刻第f-1块塔板和塔底的液相流率，H为持液量，Xf—工(k)，Xn(k)分别为k采样时刻第f-1块塔板和塔底的液相轻组分浓度，l(k)，Yf(k)分别为k采样时刻塔顶和进料板处的汽相轻组分浓度，Kl，K2，K3，K4为控制律参数，X二Xn*分别为塔顶塔底液相轻组分浓度设定值，&(k)，Xn(k)为k时刻塔顶塔底的液相轻组分浓度值Aq(k)，APjk)，分别为当前时刻高效节能精馏塔控制变量及进料热状况和精馏段压强的当前理想改变值。所述的上位机6还包括人机界面模块12，用于设定采样周期T，控制律参数K"K2，K3，K4和塔顶塔底液相轻组分浓度设定值X二X二并显示控制器的输出曲线和被控变量即塔顶塔底液相轻组分浓度的记录曲线。实施例2参照图1和图2，一种高效节能精馏过程的高纯控制方法，所述的控制方法包括以下步骤1)确定采样周期T，并将T值，相对挥发度a，提馏段压强Ps，安东尼常数a、b、c、保存在历史数据库当中；2)设定控制律参数&，K2，K3，K4和塔顶塔底液相轻组分浓度设定值X,，Xn*3)从智能仪表获取k采样时刻精馏段压强&提馏段压强Ps，以及各塔板温度I,计算液相轻组分浓度值，计算式为(1)(2):<formula>formulaseeoriginaldocumentpage10</formula>(2)其中k为当前采样时刻，下脚标i为塔板编号，1为塔顶编号，f为进料板编号，n为塔底编号，Xi(k)为k采样时刻液相轻组分浓度，Pr(k)为k采样时刻精馏段压强、Ps提馏段压强，Ti(k)为k采样时刻塔内各块塔板的温度，a为相对挥发度，a、b、c为安东尼常数。4)用历史数据库中组分推断模块计算出的组分浓度数据，在线拟合模型函数，并将拟合参数存储到历史数据库当中，拟合函数如式(3)(4):(<formula>formulaseeoriginaldocumentpage10</formula>(4)其中l,为第i块塔板处液相组分浓度预估浓度，Xi一—一—.Ss为拟合参数。5)根据当前组分浓度数据，模型函数和当前时刻操作变量值求取当前的控制变〗的理想改变值，求解控制律代数方程组如式(5)至式(12)<formula>formulaseeoriginaldocumentpage11</formula>其中k为当前采样时刻，T为采样周期，&(k)、&(k)分别为k采样时刻第i块塔板轻组分液相轻组分浓度和汽相轻组分浓度，Qi(k)为第i块塔板之间的热耦合量，UA为传热速率，Xi+f—Jk)为k采样时刻第i+f-1块塔板液相轻组分浓度，q(k)为k采样时刻进料热状况，Pr(k)为当k采样时刻精馏段压强，F为进料流率，Zf为进料组分浓度，K(k)，Vf(k)分别为k采样时刻塔顶和进料板处的汽相流率，Lf—工(k)，Ln(k)分别为k采样时刻第f-1块塔板和塔底的液相流率，H为持液量，Xf—工(k)，Xn(k)分别为k采样时刻第f-1块塔板和塔底的液相轻组分浓度，l(k)，Yf(k)分别为k采样时刻塔顶和进料板处的汽相轻组分浓度，Kl，K2，K3，K4为控制律参数，X二Xn*分别为塔顶塔底液相轻组分浓度设定值，&(k)，Xn(k)为k时刻塔顶塔底的液相轻组分浓度值Aq(k)，APjk)，分别为当前时刻高效节能精馏塔控制变量即进料热状况和精馏段压强的当前理想改变值；6)将当前时刻高效节能精馏塔控制变量即进料热状况和精馏段压强的当前理想改变值Aq(k)，APr(k)输送给DCS系统中的控制站，调整高效节能精馏塔的进料热状况值和精馏段压强值。所述的历史数据库为DCS系统中存储装置4，所述的DCS系统包括数据接口3、存储装置4和控制站5，其中现场智能仪表2与数据接口3连接，所述数据接口3与现场总线连接，所述现场总线与存储装置4、控制站5和上位机6连接，其中控制站可以读取历史数据库，显示高效节能精馏塔工作过程状态。权利要求一种高效节能精馏过程的高纯控制系统，包括与高效节能精馏塔直接连接的现场智能仪表和DCS系统，所述DCS系统包括存储装置、控制站和上位机，现场智能仪表与存储装置、控制站和上位机连接，所述的上位机包括用以计算输出高效节能精馏塔控制变量值的高纯控制器，所述高纯控制器包括组分推断模块，用以从智能仪表获取温度，压强数据，计算高效节能精馏塔各块塔板的组分浓度，并将组分浓度计算结果存储在历史数据库当中，计算式为(1)(2)<mrow><msub><mi>X</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>P</mi><mi>r</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>×</mo><mi>α</mi><mo>×</mo><msup><mn>10</mn><mrow><mfrac><mrow><msub><mi>T</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><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><mo>,</mo><mi>i</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>·</mo><mo>·</mo><mo>·</mo><mo>·</mo><mo>·</mo><mo>·</mo><mo>,</mo><mi>f</mi><mo>-</mo><mn>1</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow><mrow><msub><mi>X</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>P</mi><mi>s</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>×</mo><mi>α</mi><mo>×</mo><msup><mn>10</mn><mrow><mfrac><mrow><msub><mi>T</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><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><mo>,</mo><mi>i</mi><mo>=</mo><mi>f</mi><mo>,</mo><mi>f</mi><mo>+</mo><mn>1</mn><mo>,</mo><mo>·</mo><mo>·</mo><mo>·</mo><mo>·</mo><mo>·</mo><mo>·</mo><mo>,</mo><mi>n</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>其中k为当前采样时刻，下脚标i为塔板编号，l为塔顶编号，f为进料板编号，n为塔底编号，Xi(k)为k采样时刻液相轻组分浓度，Pr(k)为k采样时刻精馏段压强、Ps提馏段压强，Ti(k)为k采样时刻塔内各块塔板的温度，α为相对挥发度，a、b、c为安东尼常数。模型参数自适应校正拟合模块，用以采用历史数据库中组分推断模块计算出的组分浓度数据，在线拟合模型函数，并将拟合参数存储到历史数据库当中，拟合函数如式(3)(4)<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><mrow><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></mrow></mfrac><mo>,</mo><mi>i</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>·</mo><mo>·</mo><mo>·</mo><mo>·</mo><mo>·</mo><mo>·</mo><mo>,</mo><mi>f</mi><mo>-</mo><mn>1</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow><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><mrow><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></mrow></mfrac><mo>,</mo><mi>i</mi><mo>=</mo><mi>f</mi><mo>,</mo><mi>f</mi><mo>+</mo><mn>1</mn><mo>,</mo><mo>·</mo><mo>·</mo><mo>·</mo><mo>·</mo><mo>·</mo><mo>·</mo><mo>,</mo><mi>n</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>其中为第i块塔板处液相组分浓度预估浓度，Xmin，r，Xmax，r，kr，Xmin，s，Xmax，s，ks为拟合参数，Sr，Ss分别为高效节能精馏塔精馏段、提馏段液相组分浓度分布的位置；高效节能精馏过程的高纯控制律求解模块，用以根据当前组分浓度数据，模型函数和当前时刻操作变量值求取当前的控制变量的理想改变值，求解控制律代数方程组如式(5)至式(12)<mrow><msub><mi>Y</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mi>α</mi><msub><mi>X</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>t</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>t</mi><mo>)</mo></mrow><mo>+</mo><mn>1</mn></mrow></mfrac><mo>,</mo><mi>i</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>·</mo><mo>·</mo><mo>·</mo><mo>·</mo><mo>·</mo><mo>·</mo><mo>,</mo><mi>n</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow><mfencedopen=''close=''><mtable><mtr><mtd><msub><mi>Q</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mi>UA</mi><mo>×</mo><mi>b</mi><mrow><mo>(</mo><mfrac><mn>1</mn><mrow><mi>a</mi><mo>-</mo><mi>ln</mi><mo>{</mo><msub><mi>P</mi><mi>r</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><mi>ΔPr</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mi>T</mi><mo>/</mo><mo>[</mo><msub><mi>X</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>X</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>/</mo><mi>α</mi><mo>]</mo><mo>}</mo></mrow></mfrac><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mo>-</mo><mfrac><mn>1</mn><mrow><mi>a</mi><mo>-</mo><mi>ln</mi><mo>{</mo><mi>Ps</mi><mo>/</mo><mo>[</mo><msub><mi>X</mi><mrow><mi>i</mi><mo>+</mo><mi>f</mi><mo>-</mo><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>X</mi><mrow><mi>i</mi><mo>+</mo><mi>f</mi><mo>-</mo><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>/</mo><mi>α</mi><mo>)</mo></mrow></mrow></mfrac></mtd></mtr></mtable><mrow><mi>i</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>·</mo><mo>·</mo><mo>·</mo><mo>·</mo><mo>·</mo><mo>·</mo><mo>,</mo><mi>f</mi><mo>-</mo><mn>1</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow></mfenced>Vl(k)＝F(1-q(k)-Δq(k)T)(7)Ln(k)＝F(q(k)+Δq(k)T)(8)<mrow><msub><mi>L</mi><mrow><mi>f</mi><mo>-</mo><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>Σ</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>f</mi><mo>-</mo><mn>1</mn></mrow></munderover><mfrac><mrow><msub><mi>Q</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow><mi>λ</mi></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow>Vf(k)＝Vl(k)+Lf-1(k)(10)<mrow><mfrac><mrow><msub><mi>V</mi><mi>f</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><msub><mi>Y</mi><mi>f</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>L</mi><mrow><mi>f</mi><mo>-</mo><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><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>V</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><msub><mi>Y</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow><mrow><mi>H</mi><msub><mi>X</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>11</mn><mo>)</mo></mrow></mrow><mrow><mo>=</mo><msub><mi>K</mi><mn>1</mn></msub><mrow><mo>(</mo><msubsup><mi>X</mi><mn>1</mn><mo>*</mo></msubsup><mo>-</mo><msub><mi>X</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>+</mo><msub><mi>K</mi><mn>2</mn></msub><munderover><mi>Σ</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>k</mi></munderover><mrow><mo>(</mo><msubsup><mi>X</mi><mn>1</mn><mo>*</mo></msubsup><mo>-</mo><msub><mi>X</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>)</mo></mrow><mi>T</mi></mrow><mrow><mfrac><mrow><mo>-</mo><msub><mi>V</mi><mi>f</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><msub><mi>Y</mi><mi>f</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>L</mi><mi>n</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><msub><mi>X</mi><mi>n</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>L</mi><mrow><mi>f</mi><mo>-</mo><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><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><mi>F</mi><msub><mi>Z</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>12</mn><mo>)</mo></mrow></mrow><mrow><mo>=</mo><msub><mi>K</mi><mn>3</mn></msub><mrow><mo>(</mo><msubsup><mi>X</mi><mi>n</mi><mo>*</mo></msubsup><mo>-</mo><msub><mi>X</mi><mi>n</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>+</mo><msub><mi>K</mi><mn>4</mn></msub><munderover><mi>Σ</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>k</mi></munderover><mrow><mo>(</mo><msubsup><mi>X</mi><mi>n</mi><mo>*</mo></msubsup><mo>-</mo><msub><mi>X</mi><mi>n</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>)</mo></mrow><mi>T</mi></mrow>其中k为当前采样时刻，T为采样周期，Xi(k)、Yi(k)分别为k采样时刻第i块塔板轻组分液相轻组分浓度和汽相轻组分浓度，Qi(k)为第i块塔板之间的热耦合量，UA为传热速率，Xi+f-1(k)为k采样时刻第i+f-1块塔板液相轻组分浓度，q(k)为k采样时刻进料热状况，Pr(k)为当k采样时刻精馏段压强，F为进料流率，Zf为进料组分浓度，Vl(k)，Vf(k)分别为k采样时刻塔顶和进料板处的汽相流率，Lf-1(k)，Ln(k)分别为k采样时刻第f-1块塔板和塔底的液相流率，H为持液量，Xf-1(k)，Xn(k)分别为k采样时刻第f-1块塔板和塔底的液相轻组分浓度，Yl(k)，Yf(k)分别为k采样时刻塔顶和进料板处的汽相轻组分浓度，K1，K2，K3，K4为控制律参数，Xl*，Xn*分别为塔顶塔底液相轻组分浓度设定值，Xl(k)，Xn(k)k时刻塔顶塔底的液相轻组分浓度值Δq(k)，ΔPr(k)，分别为当前时刻高效节能精馏塔控制变量即进料热状况和精馏段压强的当前理想改变值。F2009101556663C00015.tif2.如权利要求1所述的高效节能精馏过程的高纯控制系统，其特征在于所述的上位机还包括人机界面模块，用以设定采样周期T，控制律参数KpK2，K3，1(4和塔顶塔底液相轻组分浓度设定值&*，Xn*，并显示控制器的输出曲线和被控变量即塔顶塔底液相轻组分浓度的记录曲线。3.—种用如权利要求1所述的高效节能精馏过程的高纯控制系统实现的高纯控制方法，其特征在于所述的控制方法包括以下步骤1)确定采样周期T，并将T值，相对挥发度a，提馏段压强Ps，安东尼常数a、b、c、保存在历史数据库当中；2)设定控制律参数KnK2，K3，K4和塔顶塔底液相轻组分浓度设定值X,，Xn*3)从智能仪表获取k采样时刻精馏段压强&提馏段压强Ps，以及各塔板温度Ti，计算液相轻组分浓度值，计算式为(1)(2):<formula>formulaseeoriginaldocumentpage3</formula>其中k为当前采样时刻，下脚标i为塔板编号，1为塔顶编号，f为进料板编号，n为塔底编号，Xi(k)为k采样时刻液相轻组分浓度，Pr(k)为k采样时刻精馏段压强、Ps提馏段压强，Ti(k)为k采样时刻塔内各块塔板的温度，a为相对挥发度，a、b、c为安东尼常数；4)用历史数据库中组分推断模块计算出的组分浓度数据，在线拟合模型函数，并将拟合参数存储到历史数据库当中，拟合函数如式(3)(4):<formula>formulaseeoriginaldocumentpage3</formula>其中《为第i块塔板处液相组分浓度预估浓度，XmiM，Xmaxj，Xmin,s，Xmax,s，ks，&，Ss为拟合参数；5)根据当前组分浓度数据，模型函数和当前时刻操作变量值求取当前的控制变量的理想改变值，求解控制律代数方程组如式(5)至式(12)<formula>formulaseeoriginaldocumentpage4</formula>其中k为当前采样时刻，T为采样周期，Xi(k)、1(k)分别为k采样时刻第i块塔板轻组分液相轻组分浓度和汽相轻组分浓度，Qi(k)为第i块塔板之间的热耦合量，UA为传热速率，Xi+f—Jk)为k采样时刻第i+f-l块塔板液相轻组分浓度，q(k)为k采样时刻进料热状况，Pr(k)为当k采样时刻精馏段压强，F为进料流率，Zf为进料组分浓度，Vjk)，Vf(k)分别为k采样时刻塔顶和进料板处的汽相流率，Lf—工(k)，Ln(k)分别为k采样时刻第f-l块塔板和塔底的液相流率，H为持液量，Xf—工(k)，Xn(k)分别为k采样时刻第f-l块塔板和塔底的液相轻组分浓度，Yjk)，Yf(k)分别为k采样时刻塔顶和进料板处的汽相轻组分浓度，Kl，K2，K3，K4为控制律参数，Xj，Xn*分别为塔顶塔底液相轻组分浓度设定值，&(k)，Xn(k)为k时刻塔顶塔底的液相轻组分浓度值Aq(k)，APjk)，分别为当前时刻高效节能精馏塔控制变量即进料热状况和精馏段压强的当前理想改变值；6)将当前时刻高效节能精馏塔控制变量即进料热状况和精馏段压强的当前理想改变值Aq(k)，APr(k)输送给DCS系统中的控制站，调整高效节能精馏塔的进料热状况值和精馏段压强值。4.如权利要求3所述的高纯控制方法，其特征在于所述的历史数据库为DCS系统中存储装置，控制站读取历史数据库，显示高效节能精馏塔工作过程状态。全文摘要一种高效节能精馏过程的高纯控制系统，包括与高效节能精馏塔直接连接的现场智能仪表和DCS系统，DCS系统包括存储装置、控制站和上位机，上位机包括用以计算输出高效节能精馏塔控制变量值的高纯控制器，高纯控制器包括组分推断模块，用以从智能仪表获取温度，压强数据，计算高效节能精馏塔各块塔板的组分浓度；模型参数自适应校正拟合模块，用以采用历史数据库中组分推断模块计算出的组分浓度数据，在线拟合模型函数；高效节能精馏过程的高纯控制律求解模块，用以根据当前组分浓度数据，模型函数和当前时刻操作变量值求取当前的控制变量的理想改变值。以及提供了一种高纯控制方法。本发明控制效果良好、控制品质理想。文档编号B01D3/14GK101708379SQ20091015566公开日2010年5月19日申请日期2009年12月29日优先权日2009年12月29日发明者刘兴高,周叶翔申请人:浙江大学