一种板式精馏过程的混合整型非线性优化设计方法

文档序号:8528331阅读:226来源:国知局
一种板式精馏过程的混合整型非线性优化设计方法
【技术领域】
[0001] 本发明涉及一种化工过程的优化设计方法,具体地说,设计一种板式精馏过程的 混合整型非线性优化设计方法。
【背景技术】
[0002] 精馏过程的优化是实现精馏过程最优经济效益的主要手段,对于提高产品的质 量、降低能耗有着重要的意义。对精馏塔的优化设计包含塔板数目、进料位置和操作条件 的选取,最小化总投资成本和操作成本。该优化过程包含与操作条件、能量等相关的连续变 量以及塔板数、进料位置等离散变量。所以精馏过程的优化问题为非线性混合整型规划问 题。通常,混合整型非线性规划(MINLP)算法可以求解这类问题。但混合整型非线性规划 (MINLP)算法求解精馏模型所需的时间依赖于模型的复杂程度。

【发明内容】

[0003] 本发明的目的是提供板式精馏过程的优化设计方法。本发明分别通过建立板式精 馏过程的广义析取规划(GDP)模型、预处理过程、模型进行求解得到最佳进料位置和最优 塔板数,给出了一种有效的苯-甲苯精馏过程的优化设计方法,解决了精馏过程同时优化 离散变量和连续变量的难题。同时本发明给出的预处理过程,为精馏模型的求解提供了恰 当的初始值和边界值,能够保证精馏模型求解的全局收敛性和全局最优解。
[0004] 本发明的目的可以通过以下技术方案实现:
[0005] -种苯-甲苯精馏过程的优化设计方法,包括以下步骤:
[0006] (1)建立精馏过程的⑶P模型:
[0007] 精馏塔的⑶P模型结构如图1所示:热量的交换发生在精馏塔两端的塔板上,这 两块塔板和进料板都被称为"固定塔板"。其余的塔板则被称为"条件塔板",其存在条件由 相应布尔变量定义。塔板超结构如图2所示。根据精馏过程体系的物系特征分别建立固定 塔板和条件塔板由物料平衡方程、气液平衡方程、摩尔分数归一化方程和能量平衡方程等 数学表达式组成的机理模型。优化的目标函数为最小化投资费用和操作费用之和的年总成 本。
[0008] (2)建立预处理过程:
[0009] 在预处理阶段,初始值和边界值可以通过系统地解决优化问题来获取。主要通过 获取一个合理的设计作为经济优化问题的初始估计。通过理论上可逆模型的计算可以得到 严格精馏模型的初始值和边界值,能够保证精馏模型求解的全局收敛性和全局最优解。给 定进料混合物的性质和组成、气化率和操作压力,用闪蒸模型求解精馏过程的气液平衡方 程。该闪蒸模型中,气相压力由Wagner's方程求解得到。气相和液相焓值由温度的多项式 求解得到。预处理过程的原理如图3所示,该过程包含了两个基本的NLP子问题。
[0010] 第一个NLP问题:根据可逆精馏段和提馏段的物料平衡方程、能量平衡方程和归 一化方程、进料段的物料平衡方程、总的气液平衡方程建立可逆精馏过程的数学模型。
[0011] 第二个NLP问题:根据精馏过程夹点模型的塔顶、塔底、重组分完全采出的部分精 馏段和轻组分完全采出的部分提馏段的物料平衡方程、能量平衡方程和归一化方程,建立 精馏过程夹点的数学模型。
[0012] (3)精馏过程⑶P模型求解过程:
[0013] ⑶P模型运用了LB0A算法,该算法通过简化的NLP子问题和MILP主问题的迭代进 行求解。GDP模型的求解步骤(图5)为:1)预处理过程的求解;2)求解GDP相关的预处理 项:求解所有塔板存在的NLP子问题和部分塔板存在的NLP子问题;3)求解运用了GDP算 法的系统综合问题。
[0014] 本发明的有益效果:
[0015] 本发明针对板式精馏过程,建立该过程的⑶P模型。并在模型求解前进行预处理, 该过程通过求解一系列的NLP问题得到了GDP模型的初始值和边界值,降低了求解的难度。 GDP模型在MINLP模型的基础上引入了逻辑变量和析取表达式,使得建模过程简单直观,同 时也使得NLP子问题的求解难度降低。可以为苯-甲苯精馏过程的优化控制提供结构参数 和操作参数的初始值,对实际生产过程具有指导意义。
【附图说明】
[0016] 图1存在固定塔板和条件塔板的GDP精馏塔结构示意图
[0017] 图2单个塔板GDP超结构:(a)固定塔板(b)条件塔板
[0018] 图3预处理过程原理图
[0019] 图4可逆产物计算原理图
[0020] 图5鞍夹点平衡控制容积的计算原理图
[0021] 图6⑶P求解步骤示意图
【具体实施方式】
[0022] 苯(C6H6)和甲苯(C7H8)都是芳烃类碳氢化合物,都是重要的有机化工原料,在化学 工业中占有极其重要的地位。苯的沸点为80. 2°C,甲苯的沸点为110. 6°C。在常温下,苯是 一种无色、有芳香气味的透明液体,易挥发;甲苯则呈无色液体状态。
[0023] 苯、甲苯等芳烃的主要生产原料为石油催化重整油和石油裂解加氢汽油。分离芳 烃的主要手段有液液萃取和萃取精馏。萃取精馏提纯苯、甲苯的工艺流程中包含苯-甲苯 精馏塔,从该塔的塔顶和塔底分别得到纯度较高的苯和甲苯。
[0024] 步骤1 :苯-甲苯精馏过程⑶P模型的建立:
[0025] 苯、甲苯、邻二甲苯的三元混合物来测试上述所提到的方法。进料中苯、甲苯、邻二 甲苯的比值为:〇. 33、0. 33、0. 34,精馏塔的压力为1.Olbar,进料流量为10mol/s产物为饱 和液体,塔板数下限11_为5,塔板数上限n_为60,产品的纯度要求为98 %。苯-甲苯精 馏塔的GDP模型结构如图1所示:热量的交换发生在精馏塔两端的塔板上,这两块塔板和进 料板都被称为"固定塔板"。其余的塔板则被称为"条件塔板",条件塔板的存在与布尔变量 有关。
[0026] C表示进料成分i的集合;N表示塔板n的集合,n= 1,2, . . .,Nmax;top表示顶部 塔板,top= 1 ;bot表示底部塔板,bot=Nmax;prod=topUbot。目标函数为最小化投资 费用和操作费用之和的年总成本1八(:。投资费用通过函数£(阶',0。。1,4,^给定。这里目 标函数表示为:
[0027]minTAC=f(NT,Dcol,Ar,Ac)+TWQC+TSQH (1)
[0028] 将固定塔板定义为pt={n:n= 1,n=Nmax},中间的塔板定义为it={n: 2 彡n彡(N-_l)}〇
[0029] 对精馏塔的直径、冷凝器传热面积、再沸器传热面积和塔板数的约束方程如下:
【主权项】
1. 一种板式精馏过程的混合整型非线性优化设计方法,其特征在于包括以下步骤: (1) 建立板式精馏过程的广义析取规划(GDP)模型: 首先根据给定的板式精馏塔的进料成分和流率,以及最终分离得到的产物的纯度要 求,建立板式精馏塔的GDP模型:热量的交换发生在精馏塔两端的塔板上,这两块塔板和进 料板都被称为"固定塔板",其余的塔板则被称为"条件塔板",条件塔板的存在与布尔变量 有关。根据精馏过程体系的物系特征,分别建立固定塔板和条件塔板由物料平衡方程、气液 平衡方程、摩尔分数归一化方程和能量平衡方程等数学表达式组成的机理模型。优化的目 标函数为最小化投资费用和操作费用之和的年总成本。 (2) 建立预处理过程: 预处理过程通过理论上可逆模型的计算可以得到严格精馏模型的初始值和边界值, 能够保证精馏模型求解的全局收敛性和全局最优解。给定进料混合物的性质和组成、气 化率和操作压力,用闪蒸模型求解精馏过程的气液平衡方程。该闪蒸模型中,气相压力由 Wager's方程求解得到。气相和液相焓值由温度的多项式求解得到。预处理过程包含了两 个基本的NLP子问题。 第一个NLP问题:根据可逆精馏段和提馏段物料平衡方程、能量平衡方程和归一化方 程、进料段的物料平衡方程、总的气液平衡方程建立可逆精馏过程的数学模型。 第二个NLP问题:根据精馏过程夹点模型的塔顶、塔底、重组分完全采出的部分精馏段 和轻组分完全采出的部分提馏段的物料平衡方程、能量平衡方程和归一化方程,建立精馏 过程夹点的数学模型。 (3) 精馏过程GDP模型求解过程: ⑶P模型的求解运用LBOA算法,该算法通过简化的NLP子问题和MILP主问题的迭代进 行求解。GDP模型的求解步骤为:1)预处理过程的求解;2)求解GDP相关的预处理项:求解 所有塔板存在的NLP子问题和部分塔板存在的NLP子问题;3)求解运用了 GDP算法的系统 综合问题。
2. 根据权利要求1所述的一种板式精馏过程的混合整型非线性优化设计方法,其特征 在于将广义析取规划(GDP)算法运用到了精馏塔的建模和优化求解过程,是以塔板和进料 板的是否存在为逻辑变量,以年费用最小为目标函数,利用GAMS进行优化求解。
3. 根据权利要求1所述的一种板式精馏过程的混合整型非线性优化设计方法,其特征 在于采用了一般预处理过程,为GDP模型的求解提供了恰当的初始值和边界值,降低了求 解的难度。
4. 根据权利要求1所述的一种板式精馏过程的混合整型非线性优化设计方法,其特征 在于所述的苯、甲苯、邻二甲苯的三元混合物精馏分离过程优化设计。
【专利摘要】本发明公开一种板式精馏过程的优化设计方法,包括以下步骤:(1)建立苯-甲苯精馏过程的广义析取规划(GDP)模型;(2)建立苯-甲苯求解预处理过程;(3)求解苯-甲苯精馏过程的广义析取规划(GDP)模型。本发明给出了苯-甲苯精馏过程的优化设计方法,在模型求解前进行了预处理,该过程通过求解一系列的非线性问题得到了广义析取规划(GDP)模型的初始值和边界值,降低了求解的难度。广义析取规划(GDP)引入了逻辑变量和析取表达式,使得建模过程简单直观,同时也使得非线性子问题的求解难度降低。
【IPC分类】G06F17-50
【公开号】CN104850719
【申请号】CN201510297410
【发明人】薄翠梅, 黄庆庆, 李俊, 陆兵, 李芳芳
【申请人】南京工业大学
【公开日】2015年8月19日
【申请日】2015年6月1日
网友询问留言 已有0条留言
  • 还没有人留言评论。精彩留言会获得点赞!
1