[發(fā)明專(zhuān)利]基于動(dòng)態(tài)數(shù)學(xué)模型的精餾塔求解方法

在審

• 申請(qǐng)?zhí)枺?/td>: 202010889422.4
• 申請(qǐng)日: 2020-08-28
• 公開(kāi)（公告）號(hào): CN112052580A
• 公開(kāi)（公告）日: 2020-12-08
• 發(fā)明（設(shè)計(jì)）人: 黃彩霞;沈井學(xué)
• 申請(qǐng)（專(zhuān)利權(quán)）人: 杭州百子尖科技股份有限公司
• 主分類(lèi)號(hào): G06F30/20
• 分類(lèi)號(hào): G06F30/20
• 代理公司: 杭州九洲專(zhuān)利事務(wù)所有限公司 33101
• 代理人: 沈敏強(qiáng)

本發(fā)明涉及一種基于動(dòng)態(tài)數(shù)學(xué)模型的精餾塔求解方法。本發(fā)明的目的是提供一種基于動(dòng)態(tài)數(shù)學(xué)模型的精餾塔求解方法，該精餾塔具有N塊塔板，由上而下依次為第1至第N塊塔板，S1、獲取用戶給定的塔頂塔底壓力線性插值得到精餾塔內(nèi)各塔板的塔壓，第j塊塔板的塔壓P_j，j＝1...N；S2、獲取假定的各塔板的T_j、V_j、L_j、x_j,i和y_j,i；T_j為第j塊塔板的溫度，V_j為第j塊塔板的氣相流量，L_j第j塊塔板的液相流量，x_j,i為第j塊塔板上組分i的液相組分，yj,i為第j塊塔板上組分i的氣相組分，i＝1...m，i代表組分編號(hào)，m為組分的數(shù)目；S3、每塊塔板調(diào)用PH閃蒸，計(jì)算各塔板新的T_j、V_j、L_j、x_j,i和y_j,i；S4、如新的Tj和Vj滿足∑(T_i,j+1?T_i,j)²＜tol并且∑(V_i,j+1?V_i,j)²＜tol，則退出迭代，計(jì)算結(jié)束。本發(fā)明適用于精餾塔仿真模型領(lǐng)域。

技術(shù)領(lǐng)域

本發(fā)明涉及一種基于動(dòng)態(tài)數(shù)學(xué)模型的精餾塔求解方法。適用于精餾塔仿真模型領(lǐng)域。

背景技術(shù)

常見(jiàn)的精餾塔算法有：流量加和法、泡露點(diǎn)法、內(nèi)外圈法等等，此類(lèi)算法都是通過(guò)MESH模型方程來(lái)求解各塔板溫度、流量及組成。但是不同的求解算法適用的體系是不同的，如流量加和法適用于沸點(diǎn)相差較大的體系；泡露點(diǎn)法適用于沸點(diǎn)相差較小的體系；內(nèi)外圈算法適用的范圍較為廣泛，但是對(duì)于一些沸程較寬的體系，內(nèi)外圈法仍然無(wú)法較好的收斂。

發(fā)明內(nèi)容

本發(fā)明要解決的技術(shù)問(wèn)題是：針對(duì)上述存在的問(wèn)題，提供一種基于動(dòng)態(tài)數(shù)學(xué)模型的精餾塔求解方法。

本發(fā)明所采用的技術(shù)方案是：一種基于動(dòng)態(tài)數(shù)學(xué)模型的精餾塔求解方法，精餾塔具有N塊塔板，由上而下依次為第1至第N塊塔板，其特征在于：

S1、獲取用戶給定的塔頂塔底壓力線性插值得到精餾塔內(nèi)各塔板的塔壓，其中第j塊塔板的塔壓P_j，j＝1...N；

S2、獲取假定的各塔板的T_j、V_j、L_j、x_j,i和y_j,i；其中T_j為第j塊塔板的溫度，V_j為第j塊塔板的氣相流量，L_j第j塊塔板的液相流量，x_j,i為第j塊塔板上組分i的液相組分，yj,i為第j塊塔板上組分i的氣相組分，i＝1...m，i代表組分編號(hào)，m為組分的數(shù)目；

S3、每塊塔板調(diào)用PH閃蒸，計(jì)算各塔板新的T_j、V_j、L_j、x_j,i和y_j,i；

S4、如新的Tj和Vj滿足∑(T_i,j+1-T_i,j)²＜tol并且∑(V_i,j+1-V_i,j)²＜tol，則退出迭代，計(jì)算結(jié)束；反之，將PH閃蒸計(jì)算出新的T_j、V_j、L_j、x_j,i和y_j,i作為下一輪的值，返回到步驟S3，迭代計(jì)算。

步驟S2中各塔板的初始溫度T_j的獲取方法，包括：