内部热耦合空分塔平衡级节能控制系统及方法.pdf

一种内部热耦合空分塔平衡级节能控制系统，包括与内部热耦合空分塔连接的现场智能仪表、以及控制站、数据库和上位机，上位机包括信号采集模块，节能控制模块，过程成：设定塔的结构参数和操作参数，指定进料空气流量初值；假定各塔板液相组成；对每一个塔板，分别计算其平衡温度和汽相组成、汽液相的焓值和汽液相流量；判断下式(4)是否成立，如果成立则继续，否则，更新各塔板液相组成；判断产品氮气、氧气的纯度是否满足约束，如果不满足则结束迭代，输出结果，如果满足则将空气进料流量增加一个迭代步长。以及提出一种内部热耦合空分塔平衡级节能控制方法。本发明在当前生产工况条件下使得内部热耦合空分塔单位能耗最小，并提高节能性。

1.  一种内部热耦合空分塔平衡级节能控制系统，包括与内部热耦合空分塔连接的现场智能仪表、以及控制站、数据库和上位机，智能仪表与控制站、数据库、上位机连接，其特征在于：所述的上位机包括：
信号采集模块，用以采集当前生产工况数据；
节能控制模块，用以节能控制，采用以下过程来完成：
1)设定塔的结构参数和操作参数，指定进料空气流量初值；
2)假定各塔板液相组成；
3)对每一个塔板，分别由泡点法计算其平衡温度和汽相组成；
4)对每一个塔板，分别计算汽液相的焓值；
5)由式(1)(2)计算各塔板的汽液相流量：
Vj+1Hj+1G+Uj-1Hj-1L+FjHjF-(Vj+SjG)HjG-(Uj+SjL)HjL-Qj=0---(1)]]>
Vj+1+Uj-1+FjG+FjL-(Vj+SjG)-(Uj+SjL)=0---(2)]]>
其中，V表示汽相流量，U表示液相流量，F表示进料流量，H^F表示进料焓值，S表示侧提流量，H^G和H^L分别是汽液相焓值，下标j-1、j、j+1分别表示第j-1、j、j+1块板，上标L表示液相，上标G表示汽相，Q表示热耦合量，由下式计算：
Q＝UAΔT                            (3)
其中，UA表示热耦合系数，ΔT表示耦合塔板间的温差；
6)判断下式(4)是否成立，如果成立，则继续7)，否则，更新各塔板液相组成，返回3)迭代；
Vj+1yi,j+1+Uj-1xi,j-1+Fjzi,j-(Vj+SjG)yi,j-(Uj+SjL)xi,j<0.0001---(4)]]>
其中，x是液相组成，y是汽相组成，z是进料组成，下标i＝1、2、3表示组分，依次对应氮、氩、氧；
7)判断产品氮气、氧气的纯度是否满足约束，如果不满足则结束迭代，输出结果，前一步的进料空气流量即为最大空气进料量，如果满足则将空气进料流量增加一个迭代步长Δ，返回2)继续迭代。

2.  如权利要求1所述的内部热耦合空分塔平衡级节能控制系统，其特征在于：所述上位机还包括：泡点法模块，用以由泡点法计算其平衡温度和汽相组成，其过程如下：

3.  1)假定塔板平衡温度；

3.  2)计算汽液平衡常数，采用以下过程完成：
lnΦiL=lnRTP(vL-bL)-bibL(1-ZL)+ξLaL(bibL-2Σmxmai,maL)/bLRT---(5)]]>
lnΦiG=lnRTP(vG-bG)-bibG(1-ZG)+ξGaG(bibG-2Σmxmai,maG)/bGRT---(6)]]>
Ki=ΦiL/ΦiG---(7)]]>
y_i＝K_ix_i                            (8)
其中，Φ表示逸度系数，上标L表示液相，上标G表示汽相，R是气体常数，T是温度，P是塔板压强，下标m＝1、2、3表示组分，依次对应氮、氩、氧，摩尔体积v、物性参数b^G、b^L、b_i、a^G、a^L、a_i，m、ξ^G、ξ^L、汽相压缩因子Z^G、液相压缩因子Z^L由物性模块计算；

3.  3)检验是否成立，成立则结束迭代，返回计算结果，否则，更新塔板平衡温度，返回3.2)继续迭代。

3.  如权利要求1或2所述的内部热耦合空分塔平衡级节能控制系统，其特征在于：所述上位机还包括：焓模块，用以计算汽液相混合焓，其过程如下：
Hi*=ci+diT+eiT2+fiT3+hiT4---(9)]]>
H*=ΣiyiHi*---(10)]]>
HG=H*-RT(1-ZG)-ξG(aG-TdaTdT)/bG---(11)]]>
HL=H*-RT(1-ZL)-ξL(aL-TdaLdT)/bL---(12)]]>
其中H_i^*表示第i个纯组分理想气体的焓值，H^*是混合物理想气体焓值，c、d、e、f、h为常数。

4.  如权利要求2所述的内部热耦合空分塔平衡级节能控制系统，其特征在于：所述上位机还包括：物性模块，用以计算物性参数，其过程如下：
ai,m=Ωai,mR2Tci,m2/Pci,m---(13)]]>
b_i＝Ω_bRT_ci/P_cia                            (14)
Tci,m=TciTcj(1-ki,m)---(15)]]>
Vci,m=0.125(Vci,m1/3+Vci,m1/3)3---(16)]]>
Z_ci，m＝0.5(Z_ci+Z_cm)                        (17)
P_ci，m＝RT_ci，mZ_ci，m/V_ci，m                (18)
Ω_ai，m＝0.5(Ω_ai+Ω_am)                    (19)
对汽相：
aG=ΣiΣmyiymai,m---(20)]]>
bG=Σiyibi---(21)]]>
令
A^G＝a^GP/R²T²                                (22)
B^G＝b^GP/RT                                  (23)
α^G＝2B^G-1                                  (24)
βG=AG-3BG-5BG2---(25)]]>
γG=2(BG3+BG2)-AGBG---(26)]]>
取初值为1-0.6P_r，用牛顿法解如下方程，即得到汽相压缩因子Z^G
ZG3+αGZG2+βZG+γG=0---(27)]]>
则，
v^G＝RT/PZ^G                                (28)
ξG=0.242536lnvG+3.561553bGvG-0.561553bG---(29)]]>
对液相：
aL=ΣiΣmxixmai,m---(30)]]>
bL=Σixibi---(31)]]>
令
A^L＝a^LP/R²T²                              (32)
B^L＝b^LP/RT                                (33)
α^L＝2B^L-1                                (34)
βL=AL-3BL-5BL2---(35)]]>
γL=2(BL3+BL2)-ALBL---(36)]]>
取初值为P_r(0.106+0.078P_r)，用牛顿法解如下方程，即得到液相压缩因子Z^L
ZL3+αLZL2+βZL+γL=0---(37)]]>
则，
v^L＝RT/PZ^L                                (38)
ξL=0.242536lnvL+3.561553bLvL-0.561553bL---(39)]]>
Ω_ai＝C_i-D_iτ+E_iτ²-W_iτ³                (40)
Ω_b＝0.070721                             (41)
τ＝0.01T                                 (42)
其中，A、B、a、β、γ、τ是中间变量，C、D、E、W是常数，T_c、P_c、V_c、Z_c分别是临界温度、压力、体积和压缩因子，P_r是对比压力，R是气体常数，k_i，m表示第i组分和第m组分的二元交互系数，k_i，m是常数，下标c表示临界点的性质，下标r表示对比态，下标i，m表示第i组分和第m组分的二元混合物，Ω_a、Ω_b是中间变量。

5.  如权利要求1或2所述的内部热耦合空分塔平衡级节能控制系统，其特征在于：所述的上位机还包括：
结果显示模块，用于将节能控制结果传给控制站进行显示，并通过现场总线将节能控制结果传递到现场操作站进行显示。

6.  一种用如权利要求1所述的内部热耦合空分塔平衡级节能控制系统实现的节能控制方法，其特征在于：所述的节能控制方法包括以下步骤：
1)设定塔的结构参数，采集塔的生产工况数据，采集当前进料空气流量作为初值；
2)假定各塔板液相组成；
3)对每一个塔板，分别由泡点法计算其平衡温度和汽相组成；
4)对每一个塔板，计算其汽液相的焓值；
5)联立式(1)(2)计算各塔板的汽液相流量：
Vj+1Hj+1G+Uj-1Hj-1L+FjHjF-(Vj+SjG)HjG-(Uj+SjL)HjL-Qj=0---(1)]]>
Vj+1+Uj-1+FjG+FjL-(Vj+SjG)-(Uj+SjL)=0---(2)]]>
其中，V表示汽相流量，U表示液相流量，F表示进料流量，H^F表示进料焓值，S表示侧提流量，H^G和H^L分别是汽液相焓值，下标j-1、j、j+1分别表示第j-1、j、j+1块板，上标L表示液相，上标G表示汽相，Q表示热耦合量，由下式计算：
Q＝UAΔT                    (3)；
6)判断式(3)是否成立，如果成立，则继续步骤7)，否则，更新液相组成，返回3)迭代；
Vj+1yi,j+1+Uj-1xi,j-1+FjGzi,jG+FjLzi,jL-(Vj+SjG)yi,j-(Uj+SjL)xi,j<0.0001---(4)]]>
其中，x是液相组成，y是汽相组成，z是进料组成，下标i＝1、2、3表示组分，依次对应氮、氩、氧；
7)判断产品氮气、氧气的纯度是否满足约束，如果不满足则结束迭代，输出结果，前一步的进料空气流量即为最大空气进料量，如果满足则将空气进料流量增加一个迭代步长Δ，返回步骤2)继续迭代。

7.  如权利要求6所述的节能控制方法，其特征在于：所述的步骤3)中，泡点法计算其平衡温度和汽相，采用以下过程完成：

3.  1)假定塔板平衡温度；

3.  2)计算汽液平衡常数，采用以下过程完成：
lnΦiL=lnRTP(vL-bL)-bibL(1-ZL)+ξLaL(bibL-2Σmxmai,maL)/bLRT---(5)]]>
lnΦiG=lnRTP(vG-bG)-bibG(1-ZG)+ξGaG(bibG-2Σmxmai,maG)/bGRT---(6)]]>
Ki=ΦiL/ΦiG---(7)]]>
y_i＝K_ix_i                            (8)
其中，Φ表示逸度系数，上标L表示液相，上标G表示汽相，R是气体常数，T是温度，P是塔板压强，下标m＝1、2、3表示组分，依次对应氮、氩、氧，摩尔体积v、物性参数b^G、b^L、b_i、a^G、a^L、a_i，m、ξ^G、ξ^L、汽相压缩因子Z^G、液相压缩因子Z^L由物性参数计算方法计算；

3.  3)检验是否成立，成立则结束迭代，返回计算结果，否则，更新塔板平衡温度，返回3.2)继续迭代。

8.  如权利要求6或7所述的节能控制方法，其特征在于：所述步骤4)中，所述焓计算过程如下：
Hi*=ci+diT+eiT2+fiT3+hiT4---(9)]]>
H*=ΣiyiHi*---(10)]]>
HG=H*-RT(1-ZG)-ξG(aG-TdaTdT)/bG---(11)]]>
HL=H*-RT(1-ZL)-ξL(aL-TdaLdT)/bL---(12)]]>
其中H_i^*表示第i个纯组分理想气体的焓值，H^*是混合物理想气体焓值，c、d、e、f、h为常数。

9.  如权利要求7所述的节能控制方法，其特征在于：所述物性参数计算方法的过程如下：
ai,m=Ωai,mR2Tci,m2/Pci,m---(13)]]>
b_i＝Ω_bRT_ci/P_cia                            (14)
Tci,m=TciTcj(1-ki,m)---(15)]]>
Vci,m=0.125(Vci,m1/3+Vci,m1/3)3---(16)]]>
Z_ci，m＝0.5(Z_ci+Z_cm)                        (17)
P_ci，m＝RT_ci，mZ_ci，m/V_ci，m                (18)
Ω_ai，m＝0.5(Ω_ai+Ω_am)                     (19)
对汽相：
aG=ΣiΣmyiymai,m---(20)]]>
bG=Σiyibi---(21)]]>
令
A^G＝a^GP/R²T²                                (22)
B^G＝b^GP/RT                                  (23)
α^G＝2B^G-1                                  (24)
βG=AG-3BG-5BG2---(25)]]>
γG=2(BG3+BG2)-AGBG---(26)]]>
取初值为1-0.6P_r，用牛顿法解如下方程，即得到汽相压缩因子Z^G
ZG3+αGZG2+βZG+γG=0---(27)]]>
则，
v^G＝RT/PZ^G                                  (28)
ξG=0.242536lnvG+3.561553bGvG-0.561553bG---(29)]]>
对液相：
aL=ΣiΣmxixmai,m---(30)]]>
bL=Σixibi---(31)]]>
令
A^L＝a^LP/R²T²                              (32)
B^L＝b^LP/RT                                (33)
α^L＝2B^L-1                                (34)
βL=AL-3BL-5BL2---(35)]]>
γL=2(BL3+BL2)-ALBL---(36)]]>
取初值为P_r(0.106+0.078P_r)，用牛顿法解如下方程，即得到液相压缩因子Z^L
ZL3+αLZL2+βZL+γL=0---(37)]]>
则，
v^L＝RT/PZ^L                                (38)
ξL=0.242536lnvL+3.561553bLvL-0.561553bL---(39)]]>
Ω_ai＝C_i-D_iτ+E_iτ²-W_iτ³                 (40)
Ω_b＝0.070721                             (41)
τ＝0.01T                                 (42)
其中，A、B、α、β、γ、τ是中间变量，C、D、E、W是常数，T_c、P_c、V_c、Z_c分别是临界温度、压力、体积和压缩因子，P_r是对比压力，R是气体常数，k_i，m表示第i组分和第m组分的二元交互系数，k_i，m是常数，下标c表示临界点的性质，下标r表示对比态，下标i，m表示第i组分和第m组分的二元混合物，Ω_a、Ω_b是中间变量。

10.  如权利要求6或7所述的节能控制方法，其特征在于：在所述的步骤7)中，上位机将节能控制结果传给控制站进行显示，并通过现场总线将节能控制结果传递到现场操作站进行显示。

内部热耦合空分塔平衡级节能控制系统及方法
技术领域
本发明涉及空分领域，特别地，涉及一种内部热耦合空分塔平衡级节能控制系统及方法。
背景技术
内部热耦合精馏技术是至今为止的四大节能精馏技术中节能效能最高的一项节能技术，在世界范围内得到了巨大的重视。
在低温空分过程中，空分塔是一个重要的操作单元，也是最重要的耗能单元。通过对内部热耦合技术的研究发现：(1)该技术的实现，需要过程装置中具有不同高低压强的区域，而空分装置本身的化工工艺流程就给出了这样区域，不需要像其他过程装置需要人为制造出不同高低压强的区域。(2)研究表明内部热耦合技术在高纯度产品生产过程上将具有更大的经济效益，而空分过程恰好是一个产品纯度要求极高的过程。(3)内部热耦合精馏技术在低操作温度情况下更加有利于发挥它的节能潜力，而本身就是低温精馏的空分装置较之其他精馏过程，无疑具有这个优势。(4)内部热耦合精馏技术在处理沸点接近的物系时，能体现出很高的节能效率，而空气中氮一氩一氧三元物系又恰好满足了这个要求，也为热耦合技术的应用开发提高了非常有利条件。因此，内部热耦合技术可以被应用到空气分离过程，改变传统空分塔结构，达到良好的节能效果。
过程在线优化是指综合应用过程建模技术、优化技术、先进控制技术以及计算机技术，在满足生产安全要求以及产品质量约束等条件下，不断计算并改变过程的操作条件，使得生产过程始终运行在“最优状态”。研究表明，在线优化所带来的经济效益，相当于DCS和各种先进控制手段所带来的经济效益的总和，其投资却只有先进控制的三分之一。可见在线优化的成功实施是企业获取最小经济效益的必要保证。内部热耦合空分塔平衡级节能控制是指基于内部热耦合空分塔平衡级模型，在满足产品纯度要求的前提下，不断计算并改变空分塔的进料流量，使得空分过程始终运行在产量最大即单位能耗最小的状态。
发明内容
为了克服现有的空分塔内部热耦合精馏过程的单位能耗较大、节能性较差的不足，本发明提供一种在当前生产工况条件下使得内部热耦合空分塔单位能耗最小，并提高节能性的内部热耦合空分塔平衡级节能控制系统及方法。
本发明解决其技术问题所采用的技术方案是：
一种内部热耦合空分塔平衡级节能控制系统，包括与内部热耦合空分塔连接的现场智能仪表、以及控制站、数据库和上位机，智能仪表与控制站、数据库、上位机连接，所述的上位机包括：
信号采集模块，用以采集当前生产工况数据；
节能控制模块，用以节能控制，采用以下过程来完成：
1)设定塔的结构参数和操作参数，指定进料空气流量初值；
2)假定各塔板液相组成；
3)对每一个塔板，分别由泡点法计算其平衡温度和汽相组成；
4)对每一个塔板，分别计算汽液相的焓值；
5)由式(1)(2)计算各塔板的汽液相流量：
Vj+1Hj+1G+Uj-1Hj-1L+FjHjF-(Vj+SjG)HjG-(Uj+SjL)HjL-Qj=0---(1)]]>
Vj+1+Uj-1+FjG+FjL-(Vj+SjG)-(Uj+SjL)=0---(2)]]>
其中，V表示汽相流量，U表示液相流量，F表示进料流量，H^F表示进料焓值，S表示侧提流量，H^G和H^L分别是汽液相焓值，下标j-1、j、j+1分别表示第j-1、j、j+1块板，上标L表示液相，上标G表示汽相，Q表示热耦合量，由下式计算：
Q＝UAΔT                        (3)
其中，UA表示热耦合系数，ΔT表示耦合塔板间的温差；
6)判断下式(4)是否成立，如果成立，则继续7)，否则，更新各塔板液相组成，返回3)迭代；
Vj+1yi,j+1+Uj-1xi,j-1+Fjzi,j-(Vj+SjG)yi,j-(Uj+SjL)xi,j<0.0001---(4)]]>
其中，x是液相组成，y是汽相组成，z是进料组成，下标i＝1、2、3表示组分，依次对应氮、氩、氧；
7)判断产品氮气、氧气的纯度是否满足约束，如果不满足则结束迭代，输出结果，前一步的进料空气流量即为最大空气进料量，如果满足则将空气进料流量增加一个迭代步长Δ，返回2)继续迭代。
作为优选的一种方案：所述上位机还包括：泡点法模块，用以由泡点法计算其平衡温度和汽相组成，其过程如下：
3.1)假定塔板平衡温度；
3.2)计算汽液平衡常数，采用以下过程完成：
lnΦiL=lnRTP(vL-bL)-bibL(1-ZL)+ξLaL(bibL-2Σmxmai,maL)/bLRT---(5)]]>
lnΦiG=lnRTP(vG-bG)-bibG(1-ZG)+ξGaG(bibG-2Σmxmai,maG)/bGRT---(6)]]>
Ki=ΦiL/ΦiG---(7)]]>
y_i＝K_ix_i                            (8)
其中，Φ表示逸度系数，上标L表示液相，上标G表示汽相，R是气体常数，T是温度，P是塔板压强，下标m＝1、2、3表示组分，依次对应氮、氩、氧，摩尔体积v、物性参数b^G、b^L、b_i、a^G、a^L、a_i，m、ξ^G、ξ^L、汽相压缩因子Z^G、液相压缩因子Z^L由物性模块计算；
3.3)检验是否成立，成立则结束迭代，返回计算结果，否则，更新塔板平衡温度，返回3.2)继续迭代。
作为优选的另一种方案：所述上位机还包括：焓模块，用以计算汽液相混合焓，其过程如下：
Hi*=ci+diT+eiT2+fiT3+hiT4---(9)]]>
H*=ΣiyiHi*---(10)]]>
HG=H*-RT(1-ZG)-ξG(aG-TdaGdT)/bG---(11)]]>
HL=H*-RT(1-ZL)-ξL(aL-TdaLdT)/bL---(12)]]>
其中H_i^*表示第i个纯组分理想气体的焓值，H^*是混合物理想气体焓值，c、d、e、f、h为常数。
作为优选的再一种方案：所述上位机还包括：物性模块，用以计算物性参数，其过程如下：
ai,m=Ωai,mR2Tci,m2/Pci,m---(13)]]>
b_i＝Ω_bRT_ci/P_cia                (14)
Tci,m=TciTcj(1-ki,m)---(15)]]>
Vci,m=0.125(Vci,m1/3+Vci,m1/3)3---(16)]]>
Z_ci，m＝0.5(Z_ci+Z_cm)            (17)
P_ci，m＝RT_ci，mZ_ci，m/V_ci，m    (18)
Ω_ai，m＝0.5(Ω_ai+Ω_am)         (19)
对汽相：
aG=ΣiΣmyiymai,m---(20)]]>
bG=Σiyibi---(21)]]>
令
A^G＝a^GP/R²T²                    (22)
B^G＝b^GP/RT                      (23)
α^G＝2B^G-1                      (24)
βG=AG-3BG-5BG2---(25)]]>
γG=2(BG3+BG2)-AGBG---(26)]]>
取初值为1-0.6P_r，用牛顿法解如下方程，即得到汽相压缩因子Z^G
ZG3+αGZG2+βZG+γG=0---(27)]]>
则，
v^G＝RT/PZ^G                    (28)
ξG=0.242536lnvG+3.561553bGvG-0.561553bG---(29)]]>
对液相：
aL=ΣiΣmxixmai,m---(30)]]>
bL=Σixibi---(31)]]>
令
A^L＝a^LP/R²T²                  (32)
B^L＝b^LP/RT                    (33)
α^L＝2B^L-1                    (34)
βL=AL-3BL-5BL2---(35)]]>
γL=2(BL3+BL2)-ALBL---(36)]]>
取初值为P_r(0.106+0.078P_r)，用牛顿法解如下方程，即得到液相压缩因子Z^L
ZL3+αLZL2+βZL+γL=0---(37)]]>
则，
v^L＝RT/PZ^L                    (38)
ξL=0.242536lnvL+3.561553bLvL-0.561553bL---(39)]]>
Ω_ai＝C_i-D_iτ+E_iτ²-W_iτ³     (40)
Ω_b＝0.070721                 (41)
τ＝0.01T                     (42)
其中，A、B、α、β、γ、τ是中间变量，C、D、E、W是常数，T_c、P_c、V_c、Z_c分别是临界温度、压力、体积和压缩因子，P_r是对比压力，R是气体常数，k_i，m表示第i组分和第m组分的二元交互系数，k_i，m是常数，下标c表示临界点的性质，下标r表示对比态，下标i，m表示第i组分和第m组分的二元混合物，Ω_a、Ω_b是中间变量。
进一步，所述的上位机还包括：结果显示模块，用于将节能控制结果传给控制站进行显示，并通过现场总线将节能控制结果传递到现场操作站进行显示。
一种内部热耦合空分塔平衡级节能控制方法，所述的节能控制方法包括以下步骤：
1)设定塔的结构参数，采集塔的生产工况数据，采集当前进料空气流量作为初值；
2)假定各塔板液相组成；
3)对每一个塔板，分别由泡点法计算其平衡温度和汽相组成；
4)对每一个塔板，计算其汽液相的焓值；
5)联立式(1)(2)计算各塔板的汽液相流量：
Vj+1Hj+1G+Uj-1Hj-1L+FjHjF-(Vj+SjG)HjG-(Uj+SjL)HjL-Qj=0---(1)]]>
Vj+1+Uj-1+FjG+FjL-(Vj+SjG)-(Uj+SjL)=0---(2)]]>
其中，V表示汽相流量，U表示液相流量，F表示进料流量，H^F表示进料焓值，S表示侧提流量，H^G和H^L分别是汽液相焓值，下标j-1、j、j+1分别表示第j-1、j、j+1块板，上标L表示液相，上标G表示汽相，Q表示热耦合量，由下式计算：
Q＝UAΔT                        (3)；
6)判断式(3)是否成立，如果成立，则继续步骤7)，否则，更新液相组成，返回3)迭代；
Vj+1yi,j+1+Uj-1xi,j-1+FjGzi,jG+FjLzi,jL-(Vj+SjG)yi,j-(Uj+SjL)xi,j<0.0001---(4)]]>
其中，x是液相组成，y是汽相组成，z是进料组成，下标i＝1、2、3表示组分，依次对应氮、氩、氧；
7)判断产品氮气、氧气的纯度是否满足约束，如果不满足则结束迭代，输出结果，前一步的进料空气流量即为最大空气进料量，如果满足则将空气进料流量增加一个迭代步长Δ，返回步骤2)继续迭代。
作为优选的一种方案：所述的步骤3)中，泡点法计算其平衡温度和汽相，采用以下过程完成：
3.1)假定塔板平衡温度；
3.2)计算汽液平衡常数，采用以下过程完成：
lnΦiL=lnRTP(vL-bL)-bibL(1-ZL)+ξLaL(bibL-2Σmxmai,maL)/bLRT---(5)]]>
lnΦiG=lnRTP(vG-bG)-bibG(1-ZG)+ξGaG(bibG-2Σmxmai,maG)/bGRT---(6)]]>
Ki=ΦiL/ΦiG---(7)]]>
y_i＝K_ix_i                                  (8)
其中，Φ表示逸度系数，上标L表示液相，上标G表示汽相，R是气体常数，T是温度，P是塔板压强，下标m＝1、2、3表示组分，依次对应氮、氩、氧，摩尔体积v、物性参数b^G、b^L、b_i、a^G、a^L、a_i，m、ξ^G、ξ^L、汽相压缩因子Z^G、液相压缩因子Z^L由物性参数计算方法计算；
3.3)检验是否成立，成立则结束迭代，返回计算结果，否则，更新塔板平衡温度，返回3.2)继续迭代。
作为优选的另一种方案：所述步骤4)中，所述焓计算过程如下：
Hi*=ci+diT+eiT2+fiT3+hiT4---(9)]]>
H*=ΣiyiHi*---(10)]]>
HG=H*-RT(1-ZG)-ξG(aG-TdaGdT)/bG---(11)]]>
HL=H*-RT(1-ZL)-ξL(aL-TdaLdT)/bL---(12)]]>
其中H_i^*表示第i个纯组分理想气体的焓值，H^*是混合物理想气体焓值，c、d、e、f、h为常数。
作为优选的再一种方案：所述物性参数计算方法的过程如下：
ai,m=Ωai,mR2Tci,m2/Pci,m---(13)]]>
b_i＝Ω_bRT_ci/P_cia                    (14)
Tci,m=TciTcj(1-ki,m)---(15)]]>
Vci,m=0.125(Vci,m1/3+Vci,m1/3)3---(16)]]>
Z_ci，m＝0.5(Z_ci+Z_cm)                (17)
P_ci，m＝RT_ci，mZ_ci，m/V_ci，m        (18)
Ω_ai，m＝0.5(Ω_ai+Ω_am)             (19)
对汽相：
aG=ΣiΣmyiymai,m---(20)]]>
bG=Σiyibi---(21)]]>
令
A^G＝a^GP/R²T²                        (22)
B^G＝b^GP/RT                          (23)
α^G＝2B^G-1                          (24)
βG=AG-3BG-5BG2---(25)]]>
γG=2(BG3+BG2)-AGBG---(26)]]>
取初值为1-0.6P_r，用牛顿法解如下方程，即得到汽相压缩因子Z^G
ZG3+αGZG2+βZG+γG=0---(27)]]>
则，
v^G＝RT/PZ^G                        (28)
ξG=0.242536lnvG+3.561553bGvG-0.561553bG---(29)]]>
对液相：
aL=ΣiΣmxixmai,m---(30)]]>
bL=Σixibi---(31)]]>
令
A^L＝a^LP/R²T²                      (32)
B^L＝b^LP/RT                        (33)
α^L＝2B^L-1                        (34)
βL=AL-3BL-5BL2---(35)]]>
γL=2(BL3+BL2)-ALBL---(36)]]>
取初值为P_r(0.106+0.078P_r)，用牛顿法解如下方程，即得到液相压缩因子Z^L
ZL3+αLZL2+βZL+γL=0---(37)]]>
则，
v^L＝RT/PZ^L                        (38)
ξL=0.242536lnvL+3.561553bLvL-0.561553bL---(39)]]>
Ωai=Ci-Diτ+Eiτ2-Wiτ3---(40)]]>
Ω_b＝0.070721                     (41)
τ＝0.01T                         (42)
其中，A、B、α、β、γ、τ是中间变量，C、D、E、W是常数，T_c、P_c、V_c、Z_c分别是临界温度、压力、体积和压缩因子，P_r是对比压力，R是气体常数，k_i，m表示第i组分和第m组分的二元交互系数，k_i，m是常数，下标c表示临界点的性质，下标r表示对比态，下标i，m表示第i组分和第m组分的二元混合物，Ω_a、Ω_b是中间变量。
进一步，在所述的步骤7)中，上位机将节能控制结果传给控制站进行显示，并通过现场总线将节能控制结果传递到现场操作站进行显示。
本发明的有益效果主要表现在：对内部热耦合空分塔进行平衡级节能控制，可以用于指导生产，在满足当前生产工况要求的前提下提高产量，降低单位产品能耗，从而提高生产效益。
附图说明
图1是本发明所提出的平衡级节能控制系统的硬件结构图。
图2是本发明所述内部热耦合空分塔结构示意图。
图3是本发明上位机的功能模块图。
具体实施方式
下面结合附图对本发明作进一步描述。
实施例1
参照图1、图2、图3，一种内部热耦合空分塔平衡级节能控制系统，包括内部热耦合空分塔1连接的现场智能仪表2、数据接口3、控制站4、数据库5以及上位机6，智能仪表2与现场总线连接、所述现场总线与数据接口3连接，所述数据接口3与控制站4、数据库5和上位机6连接，所述的上位机6包括：
信号采集模块7，用以采集当前生产工况数据；
节能控制模块8，用以节能控制，采用以下过程来完成：
1)设定塔的结构参数和操作参数，指定进料空气流量初值；
2)假定各塔板液相组成；
3)对每一个塔板，分别由泡点法计算其平衡温度和汽相组成；
4)对每一个塔板，分别计算汽液相的焓值；
5)由式(1)(2)计算各塔板的汽液相流量：
Vj+1Hj+1G+Uj-1Hj-1L+FjHjF-(Vj+SjG)HjG-(Uj+SjL)HjL-Qj=0---(1)]]>
Vj+1+Uj-1+FjG+FjL-(Vj+SjG)-(Uj+SjL)=0---(2)]]>
其中，V表示汽相流量，U表示液相流量，F表示进料流量，H^F表示进料焓值，S表示侧提流量，H^G和H^L分别是汽液相焓值，下标j-1、j、j+1分别表示第j-1、j、j+1块板，上标L表示液相，上标G表示汽相，Q表示热耦合量，由下式计算：
Q＝UAΔT                            (3)
其中，UA表示热耦合系数，ΔT表示耦合塔板间的温差；
6)判断下式(4)是否成立，如果成立，则继续7)，否则，更新各塔板液相组成，返回3)迭代；
Vj+1yi,j+1+Uj-1xi,j-1+Fjzi,j-(Vj+SjG)yi,j-(Uj+SjL)xi,j<0.0001---(4)]]>
其中，x是液相组成，y是汽相组成，z是进料组成，下标i＝1、2、3表示组分，依次对应氮、氩、氧；
7)判断产品氮气、氧气的纯度是否满足约束，如果不满足则结束迭代，输出结果，前一步的进料空气流量即为最大空气进料量，如果满足则将空气进料流量增加一个迭代步长Δ，返回2)继续迭代。
所述上位机6还包括：泡点法模块9，用以由泡点法计算其平衡温度和汽相组成，其过程如下：
3.1)假定塔板平衡温度；
3.2)计算汽液平衡常数，采用以下过程完成：
lnΦiL=lnRTP(vL-bL)-bibL(1-ZL)+ξLaL(bibL-2Σmxmai,maL)/bLRT---(5)]]>
lnΦiG=lnRTP(vG-bG)-bibG(1-ZG)+ξGaG(bibG-2Σmxmai,maG)/bGRT---(6)]]>
Ki=ΦiL/ΦiG---(7)]]>
y_i＝K_ix_i                            (8)
其中，Φ表示逸度系数，上标L表示液相，上标G表示汽相，R是气体常数，T是温度，P是塔板压强，下标m＝1、2、3表示组分，依次对应氮、氩、氧，摩尔体积v、物性参数b^G、b^L、b_i、a^G、a^L、a_i，m、ξ^G、ξ^L、汽相压缩因子Z^G、液相压缩因子Z^L由物性模块计算；
3.3)检验是否成立，成立则结束迭代，返回计算结果，否则，更新塔板平衡温度，返回3.2)继续迭代。
所述上位机6还包括：焓模块10，用以计算汽液相混合焓，其过程如下：
Hi*=ci+diT+eiT2+fiT3+hiT4---(9)]]>
H*=ΣiyiHi*---(10)]]>
HG=H*-RT(1-ZG)-ξG(aG-TdaGdT)/bG---(11)]]>
HL=H*-RT(1-ZL)-ξL(aL-TdaLdT)/bL---(12)]]>
其中H_i^*表示第i个纯组分理想气体的焓值，H^*是混合物理想气体焓值，c、d、e、f、h为常数；
所述上位机6还包括：物性模块11，用以计算物性参数，其过程如下：
ai,m=Ωai,mR2Tci,m2/Pci,m---(13)]]>
b_i＝Ω_bRT_ci/P_cia                    (14)
Tci,m=TciTcj(1-ki,m)---(15)]]>
Vci,m=0.125(Vci,m1/3+Vci,m1/3)3---(16)]]>
Z_ci，m＝0.5(Z_ci+Z_cm)                (17)
P_ci，m＝RT_ci，mZ_ci，m/V_ci，m        (18)
Ω_ai，m＝0.5(Ω_ai+Ω_am)             (19)
对汽相：
aG=ΣiΣmyiymai,m---(20)]]>
bG=Σiyibi---(21)]]>
令
A^G＝a^GP/R²T²                    (22)
B^G＝b^GP/RT                      (23)
α^G＝2B^G-1                      (24)
βG=AG-3BG-5BG2---(25)]]>
γG=2(BG3+BG2)-AGBG---(26)]]>
取初值为1-0.6P_r，用牛顿法解如下方程，即得到汽相压缩因子Z^G
ZG3+αGZG2+βZG+γG=0---(27)]]>
则，
v^G＝RT/PZ^G                      (28)
ξG=0.242536lnvG+3.561553bGvG-0.561553bG---(29)]]>
对液相：
aL=ΣiΣmxixmai,m---(30)]]>
bL=Σixibi---(31)]]>
令
A^L＝a^LP/R²T²                    (32)
B^L＝b^LP/RT                      (33)
α^L＝2B^L-1                      (34)
βL=AL-3BL-5BL2---(35)]]>
γL=2(BL3+BL2)-ALBL---(36)]]>
取初值为P_r(0.106+0.078P_r)，用牛顿法解如下方程，即得到液相压缩因子Z^L
ZL3+αLZL2+βZL+γL=0---(37)]]>
则，
v^L＝RT/PZ^L                        (38)
ξL=0.242536lnvL+3.561553bLvL-0.561553bL---(39)]]>
Ω_ai＝C_i-D_iτ+E_iτ²-W_iτ³         (40)
Ω_b＝0.070721                     (41)
τ＝0.01T                         (42)
其中，A、B、α、β、γ、τ是中间变量，C、D、E、W是常数，T_c、P_c、V_c、Z_c分别是临界温度、压力、体积和压缩因子，P_r是对比压力，R是气体常数，k_i，m表示第i组分和第m组分的二元交互系数，k_i，m是常数，下标c表示临界点的性质，下标r表示对比态，下标i，m表示第i组分和第m组分的二元混合物，Ω_a、Ω_b是中间变量。
所述的上位机6还包括：结果显示模块12，用于将节能控制结果传给控制站进行显示，并通过现场总线将节能控制结果传递到现场操作站进行显示。
本实施例的内部热耦合空分塔节能潜力优化系统的硬件结构图如附图1所示，所述的优化系统核心由包括信号采集模块7、节能控制模块8、泡点法模块9、焓模块10、物性模块11、结果显示模块12和人机界面的上位机6构成，此外还包括：现场智能仪表2，数据接口3、控制站4、数据库5和现场总线。内部热耦合空分塔1、智能仪表2、数据接口3、控制站4、数据库5、上位机6通过现场总线依次相连，实现信息流的上传和下达。节能控制系统在上位机6上运行，可以方便地与底层系统进行信息交换。
本实施例的优化系统的功能模块图如附图3所示，主要包括信号采集模块7、节能控制模块8、泡点法模块9、焓模块10、物性模块11、结果显示模块12等。
所述的平衡级节能控制方法按照如下步骤进行实施：
1)设定塔的结构参数，采集塔的生产工况数据，采集当前进料空气流量作为初值；
2)假定各塔板液相组成；
3)对每一个塔板，分别由泡点法计算其平衡温度和汽相组成；
4)对每一个塔板，计算其汽液相的焓值；
5)联立式(1)(2)计算各塔板的汽液相流量：
Vj+1Hj+1G+Uj-1Hj-1L+FjHjF-(Vj+SjG)HjG-(Uj+SjL)HjL-Qj=0---(1)]]>
Vj+1+Uj-1+FjG+FjL-(Vj+SjG)-(Uj+SjL)=0---(2)]]>
其中，V表示汽相流量，U表示液相流量，F表示进料流量，H^F表示进料焓值，S表示侧提流量，H^G和H^L分别是汽液相焓值，下标j-1、j、j+1分别表示第j-1、j、j+1块板，上标L表示液相，上标G表示汽相，Q表示热耦合量，由下式计算：
Q＝UAΔT                        (3)；
6)判断式(4)是否成立，如果成立，则继续步骤7)，否则，更新液相组成，返回3)迭代；
Vj+1yi,j+1+Uj-1xi,j-1+FjGzi,jG+FjLzi,jL-(Vj+SjG)yi,j-(Uj+SjL)xi,j<0.0001---(4)]]>
其中，x是液相组成，y是汽相组成，z是进料组成，下标i＝1、2、3表示组分，依次对应氮、氩、氧；
7)判断产品氮气、氧气的纯度是否满足约束，如果不满足则结束迭代，输出结果，前一步的进料空气流量即为最大空气进料量，如果满足则将空气进料流量增加一个迭代步长Δ，返回步骤2)继续迭代：
实施例2
参照图1、图2、图3，一种内部热耦合空分塔平衡级节能控制方法，所述的平衡级节能控制方法包括以下步骤：
1)设定塔的结构参数，采集塔的生产工况数据，采集当前进料空气流量作为初值；
2)假定各塔板液相组成；
3)对每一个塔板，分别由泡点法计算其平衡温度和汽相组成；
4)对每一个塔板，计算其汽液相的焓值；
5)联立式(1)(2计算各塔板的汽液相流量：
Vj+1Hj+1G+Uj-1Hj-1L+FjHjF-(Vj+SjG)HjG-(Uj+SjL)HjL-Qj=0---(1)]]>
Vj+1+Uj-1+FjG+FjL-(Vj+SjG)-(Uj+SjL)=0---(2)]]>
其中，V表示汽相流量，U表示液相流量，F表示进料流量，H^F表示进料焓值，S表示侧提流量，H^G和H^L分别是汽液相焓值，下标j-1、j、j+1分别表示第j-1、j、j+1块板，上标L表示液相，上标G表示汽相，Q表示热耦合量，由下式计算：
Q＝UAΔT                            (3)；
6)判断式(4)是否成立，如果成立，则继续步骤7)，否则，更新液相组成，返回步骤3)迭代；
Vj+1yi,j+1+Uj-1xi,j-1+FjGzi,jG+FjLzi,jL-(Vj+SjG)yi,j-(Uj+SjL)xi,j<0.0001---(4)]]>
其中，x是液相组成，y是汽相组成，z是进料组成，下标i＝1、2、3表示组分，依次对应氮、氩、氧；
7)判断产品氮气、氧气的纯度是否满足约束，如果不满足则结束迭代，输出结果，前一步的进料空气流量即为最大空气进料量，如果满足则将空气进料流量增加一个迭代步长Δ，返回步骤2)继续迭代。
所述的步骤3)中，泡点法计算其平衡温度和汽相，采用以下过程完成：
3.1)假定塔板平衡温度；
3.2)计算汽液平衡常数，采用以下过程完成：
lnΦiL=lnRTP(vL-bL)-bibL(1-ZL)+ξLaL(bibL-2Σmxmai,maL)/bLRT---(5)]]>
lnΦiG=lnRTP(vG-bG)-bibG(1-ZG)+ξGaG(bibG-2Σmxmai,maG)/bGRT---(6)]]>
Ki=ΦiL/ΦiG---(7)]]>
y_i＝K_ix_i                    (8)
其中，Φ表示逸度系数，上标L表示液相，上标G表示汽相，R是气体常数，T是温度，P是塔板压强，下标m＝1、2、3表示组分，依次对应氮、氩、氧，摩尔体积v、物性参数b^G、b^L、b_i、a^G、a^L、a_i，m、ξ^G、ξ^L、汽相压缩因子Z^G、液相压缩因子Z^L由物性参数计算方法计算；
3.3)检验是否成立，成立则结束迭代，返回计算结果，否则，更新塔板平衡温度，返回3.2)继续迭代。
所述步骤4)中，所述焓计算过程如下：
Hi*=ci+diT+eiT2+fiT3+hiT4---(9)]]>
H*=ΣiyiHi*---(10)]]>
HG=H*-RT(1-ZG)-ξG(aG-TdaGdT)/bG---(11)]]>
HL=H*-RT(1-ZL)-ξL(aL-TdaLdT)/bL---(12)]]>
其中H_i^*表示第i个纯组分理想气体的焓值，H^*是混合物理想气体焓值，c、d、e、f、h为常数。
所述物性参数计算方法的过程如下：
ai,m=Ωai,mR2Tci,m2/Pci,m---(13)]]>
b_i＝Ω_bRT_ci/P_cia                    (14)
Tci,m=TciTcj(1-ki,m)---(15)]]>
Vci,m=0.125(Vci,m1/3+Vci,m1/3)3---(16)]]>
Z_ci，m＝0.5(Z_ci+Z_cm)                (17)
P_ci，m＝RT_ci，mZ_ci，m/V_ci，m        (18)
Ω_ai，m＝0.5(Ω_ai+Ω_am)             (19)
对汽相：
aG=ΣiΣmyiymai,m---(20)]]>
bG=Σiyibi---(21)]]>
令
A^G＝a^GP/R²T²                    (22)
B^G＝b^GP/RT                      (23)
α^G＝2B^G-1                      (24)
βG=AG-3BG-5BG2---(25)]]>
γG=2(BG3+BG2)-AGBG---(26)]]>
取初值为1-0.6P_r，用牛顿法解如下方程，即得到汽相压缩因子Z^G
ZG3+αGZG2+βZG+γG=0---(27)]]>
则，
v^G＝RT/PZ^G                      (28)
ξG=0.242536lnvG+3.561553bGvG-0.561553bG---(29)]]>
对液相：
aL=ΣiΣmxixmai,m---(30)]]>
bL=Σixibi---(31)]]>
令
A^L＝a^LP/R²T²                    (32)
B^L＝b^LP/RT                      (33)
α^L＝2B^L-1                      (34)
βL=AL-3BL-5BL2---(35)]]>
γL=2(BL3+BL2)-ALBL---(36)]]>
取初值为P_r(0.106+0.078P_r)，用牛顿法解如下方程，即得到液相压缩因子Z^L
ZL3+αLZL2+βZL+γL=0---(37)]]>
则，
v^L＝RT/PZ^L                        (38)
ξL=0.242536lnvL+3.561553bLvL-0.561553bL---(39)]]>
Ω_ai＝C_i-D_iτ+E_iτ²-W_iτ³         (40)
Ω_b＝0.070721                     (41)
τ＝0.01T                         (42)
其中，A、B、α、β、γ、τ是中间变量，C、D、E、W是常数，T_c、P_c、V_c、Z_c分别是临界温度、压力、体积和压缩因子，P_r是对比压力，R是气体常数，k_i，m表示第i组分和第m组分的二元交互系数，k_i，m是常数，下标c表示临界点的性质，下标r表示对比态，下标i，m表示第i组分和第m组分的二元混合物，Ω_a、Ω_b是中间变量。
在所述的步骤7)中，上位机将节能控制结果传给控制站进行显示，并通过现场总线将计算结果传递到现场操作站进行显示。
本发明所提出的内部热耦合空分塔平衡级节能控制系统及方法，已通过上述具体实施步骤进行了描述，相关技术人员明显能在不脱离本发明内容、精神和范围内对本文所述的装置和操作方法进行改动或适当变更与组合，来实现本发明技术。特别需要指出的是，所有相类似的替换和改动对本领域的技术人员是显而易见的，它们都会被视为包括在本发明精神、范围和内容中。