间歇精馏过程优化问题

间歇精微过程优化问题

摘要：间歇精馈技术广泛应用于精细化工与制药工业中，用于提纯或回收高附加值的

产品。间歇精播是一个动态过程，其本身具有瞬时特性，这种性质为过程操作策略的

制定以及过程设施结构的选择供应很大的敏捷度，同时也为讨论者们通过优化综合的

方法去挖掘间歇精饰过程的经济潜能以及提高过程的分别效率供应了驱动力与宽阔的

空间。然而，间歇精储中各个决策变量的时变特性与其过程优化所包含的简单经济指

标权衡，也为间歇精储过程的优化综合问题带来了巨大的挑战。

本文旨在构建一套系统化解决间歇精微过程优化综合问题的框架。通过对间歇精

馈过程的物理模型和数学模型的分析，得到优化问题的约束条件，并依据实际状况，

提诞生产利润最大的优化目标。最终形成了混合整数动态优化模型（MIDO）,并形成

了相应的解法。文中最终依据上述建模过程，详细计算了一个苯•甲苯的二元简洁案例,

以验证优化模型的可行性和正确性。

关键词：过程优化；间歇精储；混合整数动态优化

1引言

间歇精溜，或许是世界上最古老的用于液体混合物分别的单元操作，几个世纪以

来直至今日，间歇精憎仍旧大量应用于精细化学品及特别化学品的生产当中，间歇精

储是间歇生产过程中使用最为广泛的分别技术1L间歇精储主要适用于小批量、多组分、

高浓度物系的分别，能够实现单塔分别多组分物系，并获得较高纯度的各组分产品。

并且允许进料组分浓度在很大范围内变化，还可以适用于不同分别要求的物料。其相

对于连续精饰来讲，设施简洁、敏捷性高、适应力强、适用于工况多变、一塔多用的

场合，因此在小批量、多元物系的分别过程中得到了广泛的应用⑵。

从目前来看，全世界由间歇精福过程得到的化工产品的种类和价值所占的比重还

是很大的，尤其是在精细化工方面。目前我们国家的精细化率在30%以上，精细化工

产品具有品种多、更新换代快、产量小、商品性强、产品质量要求高、设施投资小、

附加值高等特点，其中有很大一部分都是采用间歇精储的生产方式得到的同时精

细化工包含的范围很广，包括医药、农药、合成染料、有机颜料、涂料、香料与香精、

化妆品等4()多个行业和门类。这也使得间歇精镭的技术水平成为了精细化工进展的关

键因素。

目前我们国家大部分企业仍采纳传统的间歇精储技术，即设施采纳只具有精播段

的常规间歇精储塔，运行策略为恒定回流比变塔顶组成或恒定塔顶组成变回流比的常

规操作。然而，已有文献表明，传统的间歇精僧技术并不是最为经济有效的生产方式,

最佳的间歇精储运行方案是依据待处理的物料状态选择最相宜的间歇精馆塔结构，并

采纳优化的回流比操作策略，同时权衡设施投资及公用工程费用，以利润最大化为目

标而得到的。未经优化综合的常规生产模式也致使我们国家精细化工行业存在投资成

本高、生产力量小、能源消耗严峻以及产品质量差等问题，进而导致产品利润下降以

及环境污染等现象的发生。要解决以上问题，依靠计算机帮助设计、模拟，提高精细

化工过程的模型化、掌握和优化技术水平是一个关键，其中一个重要的方面便在于间

歇精偏过程的优化综合。

基于上述社会和讨论现状，本文将对间歇精储过程的优化综合问题绽开系统化的

讨论。通过文献阅读和已有基本化工学科和工程优化学科的基本学问，以间歇精馆系

统为对象，以操作时间最短、产品产量最大、过程利润最大为优化目标函数，约束条

件为已有的化工间歇精播计算模型和物性函数模型。通过对优化问题求解，获得最佳

设施参数等，由此实现间歇精储过程的优化综合。最终，本文将列举一个简洁的案例

来详细阐述这个优化问题，并进行求解，来验证方法的可行性和有效性。

2间歇精储塔的物理模型

2.1间歇精储塔概述

间歇精储是指间歇投料、间歇出料或连续出料的各种精福过程，是间歇生产过程

中一种特别重要的分别技术，经常被用来作为中间产品和最终产品分别和提纯的手段。

间歇精微的流程如图2/所示。

间歇精馀与连续精镭原理基本全都，都是采用混合物中各组分挥发力量的不同，

通过液相和气相的回流，使气、液两相逆向多级接触，在热能驱动和相平衡关系的约

束下，使得易挥发组分（轻组分）不断从液相向气相中转移，而难挥发组分（重组分）

由气相向液相中转移，从而达到使混合物得到不断分别纯化的目的。间歇精馔的操作

方式与简洁蒸播相类似，即每批次物料一次性加入釜中，然后进行加热精储。随着釜

液被加热汽化，蒸汽在塔内逐级提升至塔顶被冷凝器冷凝，其凝液一部分返回塔内作

为回流，另一部分作为塔顶产品采出。塔釜釜液量及组成随塔顶播出液不断采出而削

减和变化。当釜液组成达到设计目标时，则停止精镭，排出釜底残余液，完成一个生

产周期。

图2-1常规精循塔示意图

间歇精储是一个典型的动态过程，通过有效的时空间支配，可单塔分别多组分混

合物，而且设施简洁，能多体系共用。工艺通用性强和操作策略敏捷的优点，使间歇

精微在化学工业中具有很大的应用价值和市场前景，其内在规律的讨论也引起了越来

越多的关注与重视。对于间歇精储过程的讨论已有一百多年的历史，纵观讨论内容主

要包括：间歇精储过程的模拟计算，间歇精馀过程操作策略的优化，间歇精馔过程的

塔型设计以及同时考虑塔型设计与操作策略的过程优化综合。

2.2间歇精储过程的物理模型

为了深化探寻间歇精储过程，把握其内部规律，并精确猜测间歇精微

塔的运行状态，就必需将实际间歇精储过程抽象成为物理模型，然后运用数学手段对

模型进行求解分析。

间歇精饰的建模始于闻名的简洁Rayleigh蒸储模型，但是模型中有诸多假设，模

型的精确性有时并不令人满足。之后随着电子计算机的高速进展，人们开

头思索引入更多的限定方程：能量平衡、塔板持液、塔板水力学以及严格的物性方程

等，来使间歇精镭过程的模型更为精确。然而在很多建模过程中，由于描述间歇精微

简单模型的是刚性代数-微分方程组，且方程及变量个数随着塔板及组分数的增加而成

倍增加，导致求解耗时巨大，甚至会消失不收敛的状况，致使求解困难。因此，人们

在不同的设计阶段及不同场景中，开头关注并考虑间歇精锵过程数学模型的大小及简

单度、计算机的计算力量以及是否有可行的数学方法及求解策略。总而言之，就是需

要权衡对真实间歇精储运行场景猜测的精确性与求解数学模型所需要的计

算代价，从而去选择适当简单度的模型去满足当前的模拟要求。依据模型的简化程度

主要可以分为简捷模型、简洁模型以及严格模型⑵。

2.2.1Rayleigh模型

Rayleigh模型，是Rayleigh于1902年基于单平衡级的间歇蒸储而建立的，是最简

洁的间歇精微过程。其操作流程如图2-2所示。待处理的液体混合物原料被加入至釜

内，通过加热器加热，釜内液体开头汽化，气相由釜上方排出，经冷凝器冷凝作为产

品采出，过程中没有凝液作为液相回流打回釜中。对于任意时刻，釜液上方的气体都

与残留于釜内的液体成相平衡关系。随着精微的进行，釜液量不断削减，而两相组成

也不断变化，易挥发组分在釜液中不断下降，釜液温度随之提升。当釜液中易挥发组

分降至规定要求时，则停止精储。

假设釜中初始二元液相混合物的原料处理量为阳，液相组成为心o,则任意时刻釜

中残余的液体量3/及其组成.皿可以通过Rayleigh模型方程进行计算

ln^=Y_^_(2-1)

片/

式中，y表示任意时刻与釜液液相组成加成相平衡关系的气相组成。

任意时刻的馈出液量可以通过下式2-2计算：

HD=B()-(2-2)

气液相平衡关系则一般假设组分间相对挥发度。在肯定的操作温度范围内，近似为常

数，瞬时平衡气相组成)，可以通过式2-3进行计算：

y=——、(2-3)

图2-2Rayleigh精t留示意图

Rayleigh模型虽然简洁，只具有一个平衡级，却为今后的多平衡级间歇精储过程

模型建立打下了基础，此后的无持液模型主要都是基于Rayleigh方程(式2-1)建立的。

2.2.2间歇精微简捷模型

间歇精储过程可以被视为一系列的持续短时间的连续精播过程在不同起始点下的

拼接，也就是说，将间歇精饰塔假设为进料随时间变化的连续精饰塔。间歇精储简捷

模型正是基于这种假设而建立起来的，并且把基于连续精微过程的Fenske-Underwood-

Gillikmd(FUG)方法通过进一步的修改，成为间歇精储过程的简捷模型。

假设一座间歇精微塔在时刻初始进料状态为与、/。，运行一小段时间间隔

加=乙-0后，精饰塔在右和乙时刻下的状态可以用等价的连续精储过程来表示(如图

2-3)o

图2-3间歇精t留过程等价于一系列迂续精微过程的连接

简捷模型忽视了间歇精饰塔的持液量与各个平衡级的能量衡算方程，并且不包含

逐板计算，只关怀精微塔两端总的物料衡算，把模型中所包含的微分及代数方程数量

降至最低，所以计算简单度大大削减，求解较为简便，适用于间歇精馀过程的初步设

计、优化阶段。

2.2.3间歇精微半严格模型

间歇精微过程半严格模型是基于相对挥发度为在温度操作区间内近似为常数以及

恒摩尔流假设的基础上提出的，模型中包含了逐板计算。此外，依据模型要求，它还

包含了进一步的假设：

I）冷凝罐及塔板上为恒摩尔持液；

2）塔顶为全凝器，不考虑分凝的状况；

3）忽视塔内板间持气量；

4）气液两相在塔板上为抱负混合状态；

5）忽视热损失；

6）塔板为理论板，即塔板效率为1（）0%：

7）进料为泡点进料；

8）间歇精锚塔在绝热状态下操作；

9）忽视塔顶到塔底的压降。

相对于间歇精偏简捷模型，半严格模型开头考虑塔板持液量对整个间歇精微过程的影

响，能够体现出塔顶和塔身所持物料的动态“匕轮效应”和稳态“容量效应”，克服了

简捷模型模拟具有较大持液量的间歇精储塔与实际状况偏差较大的状况。

2.2.4间歇精储严格模型

相比半严格模型，严格模型更全面的考虑了塔板的水力学关系及压降等因素，更

为符合间歇精储过程的真实运行状态。有学者对严格模型与简化模型做了比较，结果

表明，严格模型可以捕获到因间歇精储塔内液相流股组成及流率变化引起的持液量波

动，和其对精播塔各个操作参数造成的影响，并能猜测更接近真实状态的决策变量掌

握轨迹。

然而间歇精储严格模型的缺点也特别明显，计算代价巨大，需要很大的计算量及

计算机存储空间，尤其是当把此类模型应用于过程优化综合问题中，需要多次反复调

用间歇精储过程模型，更是使得求解过程特别困难。因此间歇精福严格模型并不特别

适合于间歇精偏过程的优化综合方面，而计算代价适中，又能够满足大部分模拟精度

要求的半严格间歇精储模型则得到广泛的应用。

2.3间歇精储过程的优化目标

间歇精储的本质是一个动态过程，塔内各项指标及操作参数都始终随时间变化，

为了满足产品要求，必需通过决策变量的不断调，，形成一套完整的操作策略，来保证

分别任务的完成。回流比是间歇精微操作过程中极为重要的决策变量，大部分学者都

是通过对回流比的优化，来得到最佳的间歇精馆操作策略的。同时最佳的间歇精储操

作策略又剧烈依靠于设计者所盼望达到的最终目标，依据目标函数的不同，操作策略

的优化问题主要可以分为三类：

（1）操作时间最短。操作时间最短问题也被称为最优时间掌握问题，是指在最短

的运行周期内生产出规定的产品总量及纯度。

（2）产品产量最大。产品产量最大操作策略是指在规定的操作时间下，生产出最

多的满足纯度要求的产品。这种操作策略适合于产品唯一且运行周期严格限定的生产

工艺。

（3）过程利润最大。将操作时间最短和产品产量最大问题加以结合，得出了过程

利润最大化的目标函数形式（式2-4）

P=:品、嘈嘈-操作费用"CQY”_Q(2-4)

4二产周期，尸+G

式中，尸为利润（$/阮），q表示产品价格（$/kmol）,G表示原料价格（$/kmol）,

G表示操作费用（$小「），包括公用工程费用、人力成本、折旧费等。小表示运行周

期（hr）,4表示预备时间（hr）,包括操作周期间隔中的进料、卸料、清洗等时间消耗。

由此，产品产量与运行周期都被包含在目标函数中进行优化权衡。

3间歇精储过程的优化问题

通过上一章分析，可以知道，间歇精储过程的最优化综合问题，可以概括如下：

目标函数：

单位时间生产的利润最大化。

约束条件：

原料处理量和组成；

产品产量和纯度要求；

年度生产时间。

同时需要求解：

（1）塔型选择：包括进料位置、采出位置、塔板、储罐连接方式；

（2）塔板结构设计：塔板数目和储罐数目、塔板尺寸；

（3）操作策略：回流比（对于单块塔板为液气比）、塔釜再沸器加热功率（或塔

釜气相流率）、运行周期〔或年度生产批次数量）。

3.1间歇精储塔的数学模型

本文中，我们采纳详细的间歇精微动态模型，针对各个平衡级和全部的混合器、

安排器考虑质量衡算方程、能量衡算方程、气液相平衡方程以及归一方程（MESH方程）,

进行数学建模网。首先，定义3个变量，分别为塔板数集合、储罐数集合以及蛆分数集

合。其中塔板数集合设定为NT=[1,2,・・・，〃},表示共有〃块塔板可用于间歇精情的

分别过程，并用下角标i来表示与第i块塔板相关联的各个变量；储罐数集合为N8二{1,

2,・・・，什，表示共有攵个储罐可用于间歇精储过程的进料或持液，并用下角标b来表

示与第b个储罐相关联的各个变量；组分数集合则用C={1,2,・••卜来表示，并用下

角标）来表示与第，个组分相关联的各个变量。需要补充的是，如前所述，由于各个气

液流股的安排港与混合器数量也都和与其具有连接关系的各个平衡级数量相对应，故

NT={\,2,・・・，〃}与NB二{1,2,・・・，。也分别作为相应混合器与安排器的数量集合。

3.1.1安排网络模块

在安排网络模块中，主要包含各个气、液相流股的安排器与混合器。在间歇精储

过程优化综合问题中，我们所要确定的就是它们之间的连接方式，这种连接方式也正

对应着各个平衡级之间的连接结构，表达了间歇精储塔的类型与结构设计方案。

各个安排器与混合器之间的连接方式则需要通过相应的连接流股状态进行表达，

主要包括该流股的流率及组成。在分别过程中，各个流股分为两种相态，即气相与液

相，故我们分别采纳匕L来表达气、液相流股的流率°相应组分的流股流量则用下标

，来表达。由于流股依附于其所在的单元而存在，下标i或b则表达了与其相连的操作

单元，而每个流股又都具有方向性，故通过上标％或。，〃来说明流股的流进或流出。

详细而言，表达了从第力个储罐中流出的气相流率，4。衅，则表达了流入第i块

塔板的/组分液相流股。

对于安排网络上侧的安排器，在任一时刻，物料衡算方程⑺如下：

%rbsNB(3-1)

Z>=l

nk

Fi,i'wNTZbsNB(3-2)

r=ib=i

y-j=之(％x九)+乞(％xyhJ)〃，£NTZbeNBHjeC(3-3)

r-ife=i

%=z(Lpx+/XxhJ)V/,i'eNTRbwNBZjwC(3-4)

r=ib=\

甲X//")=£(%x//(L))+£(%x/『(％))V7,rGNT,"bGNB(3-5)

r=ib=i

(伉))+Z(%x//(%))V/,fGNT,VbGNB(3-6)

r=ih=\

对于安排网络下侧的流股安排器，在任一时刻，其质量及能量衡算方程如下：

匕产=力匕,WsNTRbwNB(3-7)

t=I

=2"+乃TieNTRbeNB(3-8)

V"汽(％x券JXfiGNTRbeNB、\/jeC(3-9)

r=i

%•=f(4,3NTRbsNBRjeC(3-10)

/=!

v

Yxhv(vr)=之伍JXh(%))V/eNTRbwNB(3-11)

r=i

以"x".(以"卜t&jx/(凡))+O///(a)riwNTRbeNB(3-12)

3.1.2平衡级操作算子模块

所谓平衡级算子模块，其主要表示了各个储罐与塔板上的气液相流股传质、传热

过程。对于任一平衡级，其中包括有持液，并分别含有两条气液相流股的流入流出，

如图3-1所示，分别为储罐与塔板的平衡级示意图网。

图3-1储罐(左)和塔板(右)平衡级示意图

首先，储罐的质量平衡与能量衡算方程可表达如下：

四」琛+V；f}YPbwNB(3-13)

dt

d(Bjj)=琛j+嚷一V："；-琛jV力wNBRjwC(3-14)

"叫⑻)=琛心(琛)+吸力(啜)-啜W(啜")一以W(以")+QKJ-V。eNB,刃£C

(3-15)

式中。依、。1.力分别表示第〃个储罐中再沸器或冷凝器的热负荷。

对于塔板的质量和能量衡算方程可表示如下：

鬻二窗+匕加一匕加一以VieNT(3-16)

d(H：J=以j+V：：-V：；-琮;V/wNTRjwC(3-17)

"(〃/〃(”,))=/濯(")_VinhL”)_VoU，hL)_《"#(邛")+QR(+0clV/cNT,V/>C

(3-18)

同理，式中Q/、Qc,i分别表示第，块塔板上中间再沸器或中间冷凝器的热负荷。

当然，各个平衡级及混合器和安排器还具有气液相平衡方程和归一方程的约束，

表达如下：

%=-%^ieNTyjeC(3-19)

见1=入.PbwNBRjwC(3-20)

力

目x.y=1ViwNfeC(3-21)

c

zyj=1FiwNTRjwC(3-22)

:

z

xhj=1BbeNBRjeC(3-23)

;=l

力券j=lYbeNBRjeC(3-24)

j=i

从上述数学模型可以发觉，变量匕方、%、匕八4〃表达了间歇精馈过程中塔板

与储罐之间的全部连接方式，匕•»<「、L•»•"、匕•»,••、•则表达了各个塔板之间全部连接关

系。通过对这些连接流股变量在不同时间间隔内的分析，可以得到间歇精镭过程优化

综合问题中的相关信息。

3.1.3目标函数

在上•章中对间歇精储过程优化综合问题的目标分析中，我们知道要想得到间歇

精储过程的最优化流程要在设施投资、公用工程消耗、运行周期与产品产量之间做出

最佳的权衡，转化为经济指标即为单位生产时间或年度生产时间尺度下的设旅费用、

操作费用以及产品收益。由于这几者之间存在着相互影响、相互制约的关系，所以要

在目标函数中将这几项全部包含进去。则目标函数形式可表示为下式3-25,以单位时

间为基准⑼：

。匕一综C0T.A.C

min”=——--------------------(3-25)

Hrs

式中表示一个运行周期内的产品利润，丸。0为原料成本，Hrs为年度可

7=1

用生产时间，r.AC则为年度总费用，其计算式可表达如下：

上式中，表示设施费用，精储塔部分与塔板数N,、塔板直

径。相关联，换热器部分则与塔板或储罐上所需再沸器或冷凝器换热面积

心―,4出相关联。c(XQc+Z&j+G(Z-+E。/则代表了公用

\/eATIwNT7\/€/VTtwNT/

工程费用，其中&..〃、。却、。心、。心.分别表示第i块塔板和第8个储罐所需的冷凝

器或再沸器的热负荷，C,、Q则表示单位冷、热公用工程的价格。

3.1.4求解策略

间歇精/离过程的优化综合问题经数学模型构建，形成混合整数动态优化问题

(MixedIntegerDynamicOptimization,MIDO),这种数学规划模型常被用于化学工业领

域具有时变动态特性的设计变量与操作变量的同时优化问题当中。对于混合整数动态

优化问题，通常的求解策略是通过离散方法(包括差分法、正交配置分以及有限元法等

等)将微分方程分解为代数方程，由此形成混合整数非线性规划(Mixedinteger

Non-linearProgramming,MINLP)问题，进而求解。然而混合整数非线性规划(MINLP)

的求解仍旧特别困难，其中包含了两个子问题，分别为混合整数规划(MixedInteger

Programming,MIP)与非线性规划(Non-LinearProgramming,NLP)。混合整数规划具有

多种组合优化的特性，而非线性规划则面临着求解空间非凸以及目标函数与约束条件

非线性的挑战C当二者的求解难度交织在一-起.使得混合整数非线性规划的求解难度

大大增加，并无法保证命中全局最优解。因此，我们考虑能否将数学模型中的整数变

量剔除，通过连续的变量对其进行表达，由此便可将混合整数动态优化问题(MIDO)转

化为动态优化问题，再通过离散方法，转化为非线性规划问题(NLP),如此便可在不牺

牲求解精度的前提下，大大降低求解难度，提高求解效率。

4苯-甲苯二元混合模型验证

为确保本文所提出数学模型与求解策略的正确性，在本节中，我们将分别采用

AspenPlus工业流程模拟软件以及文献中的算例对所提数学模型及求解策略进行验证。

首先我们将采用AspenBatchDistillation工业流程模拟软件对常规间歇精僧塔茉.

甲苯二组分混合物分别过程进行模拟计算，再采用所提数学模型及求解策略对相同工

况进行求解，之后再将二者分别得到的结果进行比照"3。案例在常压下操作，以全回

流状态为过程起点，其初始条件如表4-1：

表4-1案例初始条件

项目条件

塔板数（包括塔釜和冷凝器）12

冷凝器持液量（kmol）5

塔板持液量（kmol）3

塔釜进料量（kmol）65

1,1,0.999,0,999,0.996,0.991,0.976,

各塔板起点轻组分浓度分布

0.940,0.359,0.704,0.487,0.279

热力学模型抱负体系

塔釜加热功率（GJ/Hr）2

塔顶冷凝器类型全凝

恒定回流比5

产品纯度要求塔釜中甲苯摩尔浓度大于95%

经过分别使用AspenPlus工业流程模拟软件以及上述数学规划建模方法模拟计算,

我们可以得到如下对比结果：

图4/塔釜持液中各组分液相组成变化趋势

图4-2塔釜持液随时间的变化趋势

图4-1与图4-2分别表达了通过两种模拟途径得到的塔釜持液中各组分液相组成随

时间变化趋势，以及塔釜再沸器持液罐中持液量随时间变化趋势。从以上两幅图中可

以看出，通过GAMS求解所得结果与经AspenBatchDistillation软件模拟所得结果特别

接近，经数值分析，误差范围掌握在（）.7%之内。比照结果验证了上述数学模型和求解

算法的可行性和正确性。

5总结

对间歇精储过程做出整体分析，提出过程综合问题的优化变量，主要包括间歇精

播系统塔型选择、结构设计方案以及操作策略的制定。通过变量分析得出，设计变量

与操作变量之间存在简单的经济指标权衡，必需在同一框架内同时考虑。基于以上规

律，文章详细分析了间歇精饰过程的物理模型和数学模型，形成了模型方程，并提出

了方程的解法。

间歇精储过程的求解过程是一个典型的优化问题，在已知的给定条件的约束下，以生

产利润最大化为优化目标，通过求解这样的优化问题得到全局最优解。在这个解的基

础上，设计塔结构，回流比等参数。文章中着重介绍了混合整数动态优化模型（MIDO）,

采纳综合目标函数权衡操作时间、流程空间，包括过程设施、工艺物料和公用工程等

要素，同时获得最佳间歇精储塔类型选择，结构设计方案和过程操作策略，由此实现

了间歇精播过程的优化综合。