发明名称 空分塔生产潜力优化系统及方法
摘要 一种空分塔生产潜力优化系统,包括与空分塔连接的现场智能仪表以及控制站、数据库和上位机,上位机包括优化计算模块,采用以下过程:设定塔的结构参数和操作参数,指定进料空气流量初值;假定各塔板液相组成;对每一个塔板,分别计算其平衡温度和汽相组成;对每一个塔板,分别计算汽液相的焓值;计算各塔板的汽液相流量;判断式(3)是否成立,如果成立,则继续,否则,更新各塔板液相组成;判断产品氮气、氧气的纯度是否满足约束,如果不满足则结束迭代,输出结果,如果满足则将空气进料流量增加一个迭代步长Δ,返回继续迭代。以及提出了一种生产潜力优化方法。本发明获取使空分塔产品纯度满足生产要求且生产能力最大、并提高节能性。
申请公布号 CN101776893B 申请公布日期 2011.12.28
申请号 CN200910157175.2 申请日期 2009.12.23
申请人 浙江大学 发明人 刘兴高
分类号 G05B19/418(2006.01)I 主分类号 G05B19/418(2006.01)I
代理机构 杭州天正专利事务所有限公司 33201 代理人 王兵;王利强
主权项 1.一种空分塔生产潜力优化系统,包括与空分塔连接的现场智能仪表、以及控制站、数据库和上位机,智能仪表与控制站、数据库和上位机连接,其特征在于:所述的上位机包括:优化计算模块,用以优化计算,采用以下过程来完成:1)设定塔的结构参数和操作参数,指定进料空气流量初值;2)假定各塔板液相组成;3)对每一个塔板,分别由泡点法计算其平衡温度和汽相组成;4)对每一个塔板,分别计算汽液相的焓值;5)由式(1)(2)计算各塔板的汽液相流量:<img file="RE-FSB00000586991400011.GIF" wi="1447" he="78" /><img file="RE-FSB00000586991400012.GIF" wi="1447" he="82" />其中,V表示汽相流量,U表示液相流量,F表示进料流量,H<sup>F</sup>表示进料焓值,S表示侧提流量,H<sup>G</sup>和H<sup>L</sup>分别是汽液相焓值,下标j-1、j、j+1分别表示第j-1、j、j+1块板,上标L表示液相,上标G表示汽相,Q表示塔板传出的能量;6)判断下式(3)是否成立,如果成立,则继续7),否则,更新各塔板液相组成,返回3)迭代;<img file="RE-FSB00000586991400013.GIF" wi="1447" he="73" />其中,x是液相组成,y是汽相组成,z是进料组成,下标i=1、2、3表示组分,依次对应氮、氩、氧;7)判断产品氮气、氧气的纯度是否满足约束,如果不满足则结束迭代,输出结果,前一步的进料空气流量即为最大空气进料量,如果满足则将空气进料流量增加一个迭代步长Δ,返回2)继续迭代;所述上位机还包括:泡点法模块,用以由泡点法计算其平衡温度和汽相组成,其过程如下:3.1)假定塔板平衡温度;3.2)计算汽液平衡常数,采用以下过程完成:<img file="RE-FSB00000586991400014.GIF" wi="1553" he="164" /><img file="RE-FSB00000586991400015.GIF" wi="1478" he="165" /><img file="RE-FSB00000586991400016.GIF" wi="1467" he="63" />y<sub>i</sub>=K<sub>i</sub>x<sub>i</sub>                               (7)其中,Φ表示逸度系数,上标L表示液相,上标G表示汽相,R是气体常数,T是温度,P是塔板压强,下标m=1、2、3表示组分,依次对应氮、氩、氧,摩尔体积v、物性参数b<sup>G</sup>、b<sup>L</sup>、b<sub>i</sub>、a<sup>G</sup>、a<sup>L</sup>、a<sub>i,m</sub>、ξ<sup>G</sup>、ξ<sup>L</sup>、汽相压缩因子Z<sup>G</sup>、液相压缩因子Z<sup>L</sup>由物性模块计算; 3.3)检验<img file="RE-FSB00000586991400021.GIF" wi="345" he="132" />是否成立,成立则结束迭代,返回计算结果,否则,更新塔板平衡温度,返回3.2)继续迭代;所述上位机还包括:焓模块,用以计算汽液相混合焓,其过程如下:<img file="RE-FSB00000586991400022.GIF" wi="1438" he="78" /><img file="RE-FSB00000586991400023.GIF" wi="1439" he="113" /><img file="RE-FSB00000586991400024.GIF" wi="1463" he="125" /><img file="RE-FSB00000586991400025.GIF" wi="1461" he="125" />其中<img file="RE-FSB00000586991400026.GIF" wi="56" he="59" />表示第i个纯组分理想气体的焓值,H<sup>*</sup>是混合物理想气体焓值,c、d、e、f、h为常数;所述上位机还包括:物性模块,用以计算物性参数,其过程如下:<img file="RE-FSB00000586991400027.GIF" wi="1537" he="90" />b<sub>i</sub>=Ω<sub>b</sub>RT<sub>ci</sub>/P<sub>ci</sub>a                     (13)<img file="RE-FSB00000586991400028.GIF" wi="1543" he="74" /><img file="RE-FSB00000586991400029.GIF" wi="1832" he="86" />Z<sub>ci,m</sub>=0.5(Z<sub>ci</sub>+Z<sub>cm</sub>)                 (16)P<sub>ci,m</sub>=RT<sub>ci,m</sub>Z<sub>ci,m</sub>/V<sub>ci,m</sub>         (17)Ω<sub>ai,m</sub>=0.5(Ω<sub>ai</sub>+Ω<sub>am</sub>)              (18)对汽相:<img file="RE-FSB000005869914000210.GIF" wi="1769" he="104" /><img file="RE-FSB000005869914000211.GIF" wi="1817" he="109" />令A<sup>G</sup>=a<sup>G</sup>P/R<sup>2</sup>T<sup>2</sup>                         (21)B<sup>G</sup>=b<sup>G</sup>P/RT                           (22)α<sup>G</sup>=2B<sup>G</sup>-1                           (23)<img file="RE-FSB000005869914000212.GIF" wi="1542" he="64" /><img file="RE-FSB000005869914000213.GIF" wi="1542" he="64" />取初值为1-0.6P<sub>r</sub>,用牛顿法解如下方程,即得到汽相压缩因子Z<sup>G</sup><img file="RE-FSB000005869914000214.GIF" wi="1542" he="62" />则,v<sup>G</sup>=RT/PZ<sup>G</sup>                           (27)<img file="RE-FSB000005869914000215.GIF" wi="1543" he="127" />对液相:<img file="RE-FSB00000586991400031.GIF" wi="1541" he="93" /><img file="RE-FSB00000586991400032.GIF" wi="1541" he="94" />令A<sup>L</sup>=a<sup>L</sup>P/R<sup>2</sup>T<sup>2</sup>                    (31)B<sup>L</sup>=b<sup>L</sup>P/RT                      (32)α<sup>L</sup>=2B<sup>L</sup>-1                      (33)<img file="RE-FSB00000586991400033.GIF" wi="1541" he="63" /><img file="RE-FSB00000586991400034.GIF" wi="1541" he="73" />取初值为P<sub>r</sub>(0.106+0.078P<sub>r</sub>),用牛顿法解如下方程,即得到液相压缩因子Z<sup>L</sup><img file="RE-FSB00000586991400035.GIF" wi="1540" he="64" />则,v<sup>L</sup>=RT/PZ<sup>L</sup>                      (37)<img file="RE-FSB00000586991400036.GIF" wi="1541" he="118" />Ω<sub>ai</sub>=C<sub>i</sub>-D<sub>i</sub>τ+E<sub>i</sub>τ<sup>2</sup>-W<sub>i</sub>τ<sup>3</sup>       (39)Ω<sub>b</sub>=0.070721                   (40)τ=0.01T                       (41)其中,A、B、α、β、γ、τ是中间变量,C、D、E、W是常数,T<sub>c</sub>、P<sub>c</sub>、V<sub>c</sub>、Z<sub>c</sub>分别是临界温度、压力、体积和压缩因子,P<sub>r</sub>是对比压力,R是气体常数,k<sub>i,m</sub>表示第i组分和第m组分的二元交互系数,k<sub>i,m</sub>是常数,下标c表示临界点的性质,下标r表示对比态,下标i,m表示第i组分和第m组分的二元混合物,Ω<sub>a</sub>、Ω<sub>b</sub>是中间变量。
地址 310027 浙江省杭州市西湖区浙大路38号
您可能感兴趣的专利