基于分段副产物瞬时选择性的反应器网络综合❋

胡仰栋， 周 勇， 伍联营

(中国海洋大学化学化工学院，山东 青岛 266100)

反应器网络综合是化工过程综合的重要组成部分，对化工过程开发和反应器工艺的优化设计等具有重要意义。依据动力学特征，Chitra 等[1-2]、赵文[3]以及本文作者[4]都将反应体系分为三类。虽然分类方法不同，但对第三类反应体系的描述都是相同的，即允许反应物到目的产物的过程中含多个中间产品的反应体系。针对第一二类简单反应体系已发展了多种综合方法，如依据反应过程化工知识和经验进行反应器网络综合的直接推断法[5-7]；可得区法由Horn[8]最先提出，Hildebrandt等[9-10]采用几何法更清楚的阐述了可得区的基本思想，赵文等[3,11]将可得区法用于求解反应器网络废料最少的问题，Glasser等[12-13]在可得区内引入温度变量，深入研究了可逆放热反应在不同换热操作下的综合问题，但对复杂反应体系可得区法的应用受维数限制。Kokossis 等[14-15]提出了针对等温和非等温问题的超结构优化法，Marquardt等[16]提出流程和分散控制系统同时优化的超结构法，及耿宁[17]提出基于功能性催化剂排列的反应路径和反应器网络同步优化的超结构法，超结构法虽不受维数限制，但开发合适的数学算法困难。胡仰栋等[18-19]针对简单反应体系提出了用目标函数瞬时值的导数与反应器选型之间的关系进行反应器网络综合的导数分析法，并应用于几个实例证明了导数分析法的有效性和简捷性。

本文发展了适用于复杂反应体系反应器网络综合的导数分析法，提出以反应步数为依据的新分段策略，对分段副产物瞬时选择性进行导数分析选择各阶段合适的反应器，以目的产物总选择性为目标函数优化反应器网络。

1 导数分析法

导数分析法是用目标函数瞬时值的导数与反应器选型之间的关系进行反应器网络综合的方法[18]。

2 第三类反应体系的反应器网络综合

2.1 反应分段策略

Denbigh反应式：

串联反应过程为A→B→D，A为关键反应物，D为最终的目的产物。从反应物A到目的产物D的反应步数为两步，故将Denbigh反应分为A→B和B→D两个阶段。将反应物相同的副反应归为同一阶段，则最终第一阶段为A→B和A→C，第二阶段为B→D和B→E。

当反应的中间产品增多，如A→B→C→D→E的串联反应，其中A为反应物，E为目的产物。从反应物到目的产物的反应步数为四步，则将反应分为A→B，B→C，C→D和D→E四个阶段。

2.2 过程参数分析与反应器网络综合

由于复杂反应的中间产品和副产物种类繁多，原简单反应体系速率瞬时值或目的产物瞬时选择性不再适用，须寻找新目标函数瞬时值描述复杂反应的过程特性。本文所研究的反应均为恒温、恒容过程。

通过Denbigh反应为例介绍本文设计的四种目标函数瞬时值S，分别是反应在不分段时目的产物对关键反应物的瞬时选择性S1，中间产品和目的产物共同组成目的产物对关键反应物的瞬时选择性S2，以及反应分段时两个阶段目的产物瞬时选择性S31和S32，两个阶段副产物瞬时选择性S41和S42。对这四种目标函数瞬时值导数分析后选择的反应器进行比较，选择最合适的参数进行反应器网络综合。四种不同瞬时选择性S及S～XA曲线如下所示：

(1)

(2)

(3)

(4)

图1 不分段Denbigh反应S～XA的曲线Fig.1 S～XA curve for Denbigh reaction

图2 分段Denbigh反应S～XA的曲线Fig.2 S～XA curve for segmented Denbigh reaction

由S41和S31均选择PFR可知S1选择单个CSTR反应器不是最优的，优化结果也证明了该点，S2存在负值，故不分段时的S1和S2均无法准确描述反应过程。S41和S42选择了各阶段合适的反应器，但S31存在负值且S32在一阶段内确定了两个反应器，故分段副产物瞬时选择性优于分段目的产物瞬时选择性。

分析分段目的产物瞬时选择性不如分段副产物瞬时选择性准确的原因，是中间产品参与两个阶段的反应，将其作为某一阶段的过程参数是不准确的。如S31中dCB的变化与A的反应速率r1和D、E的生成速率r3、r4有关，仅将其作为第一阶段参数会出现负值则无法准确描述第一阶段的反应过程。而分段副产物瞬时选择性无此类问题，因为副产物仅由各阶段的反应物生成，分段副产物瞬时选择性仅与该阶段反应物有关。如第一阶段副产物C仅由A反应生成，S41仅与CA有关，其导数分析选择的是该阶段最优的反应器。故本文以分段副产物瞬时选择性作为参数，对其采用导数分析进行反应器网络综合。

采用分段副产物瞬时选择性导数分析法得到Denbigh反应的最优反应器网络为PFR+CSTR。当反应更复杂时，本方法能确保各阶段为最优的反应器，串联得到最优反应器网络。故本方法得到的最优反应器网络是CSTR和PFR组成的串联结构，最优反应器网络和工艺参数由这些串联结构依据不同目标优化得到。

2.3 工艺参数优化

由分段副产物瞬时选择性的导数分析得出反应器网络结构后，需要依据不同的工业目标优化工艺参数。如本文优化模型的目标函数为最终目的产物对关键反应物的总选择性St最大，决策变量为各反应器出口处关键反应物的转化率XAn。

目标函数：

(5)

式中：CA表示关键反应物A的浓度，CD表示目的产物D的浓度；CA0和CD0表示进料时A和D的浓度；CAexit和CDexit表示反应器网络出口处A和D的浓度；φ与计量系数有关，本文取为1。

约束条件：

(6)

PFR:

(7)

CSTR:

(8)

本文采用遗传算法优化反应器网络[21]。优化时种群数为50，交叉和变异概率分别为0.8和0.01，采用轮盘赌选择父辈并进行交叉和变异操作，产生的子代个体不断重复以上操作直至St最优。其中第一个转化率的边界为[0,1)，第二个转化率的边界条件为(XA1,1)。

3 算例研究

3.1 Denbigh反应

其反应动力学参数如下：

k1= 1l/(mol·s),k2= 0.6 s-1,k3= 0.6 s-1,k4= 0.1l/(mol·s)

采用遗传算法对不同反应器网络进行优化，并对结果进行比较。优化结果列于表1，其中XA1和XA2为各反应器出口处A的总转化率。图3、4分别为PFR+DSR和循环PFR+CSTR结构图。

表1 不同反应器网络下Denbigh反应优化结果Table 1 Denbigh reaction optimization results under different reactor networks

从表1可知，PFR+CSTR结构的总选择性大于其他结构，故PFR+CSTR为Denbigh反应的最优反应器网络，且本文的总选择性比文献值更高。

图3 PFR+DSR结构Fig.3 Structure of PFR+DSR

图4 循环PFR+CSTR结构Fig.4 Structure of circulation PFR+CSTR

3.2 九组分反应

九组分反应式和动力学参数如下：

k1=5l/(mol·s)，k2=0.5 s-1，

k3=5 s-1，k4=0.5l/(mol·s)，

k5=5l/(mol·s)，k2=0.5 s-1，

k3=5 s-1，k4=0.5l/(mol·s)。

其中：A为关键反应物；E为最终目的产物；B、C、D为中间产品；F、G、H、I为副产物。反应进料为纯A，CA0=10 mol/L。

先将反应按新分段策略分为四个阶段，第一阶段为A→B和A→F，第二阶段为B→C和B→G，第三阶段为C→D和C→H，第四阶段为D→E和D→I。

计算各阶段的副产物瞬时选择性，依据导数分析确定最终的反应器网络结构为：PFR+CSTR+PFR+CSTR。各阶段瞬时选择性结果如图5所示：

图5 九组分反应各阶段S～XA曲线Fig.5 S Profiles in each stage of the nine-component reaction

用遗传算法求解不同反应器网络进行求解，并对结果进行比较，不同结构的最优结果列于表2，其中XA1,XA2,XA3,XA4,XA5为对应各反应器出口处A的总转化率。

表2 不同反应器网络下九组分反应优化结果Table 2 Nine component system optimization results under different reactor networks

从表2可知，PFR+CSTR+PFR+CSTR结构的总选择性优于其他结构，且比文献结果更优，说明本方法得到的PFR+CSTR+PFR+CSTR结构是最优反应器网络。第四阶段的反应器内A变化很小，舍去CSTR后的PFR+CSTR+PFR结构也能获得较高的总选择性。

3.3 Van de Vusse反应

Van de Vusse反应式是第二类反应体系中的典型反应，其反应式如下：

k1=10 s-1,k2=1 s-1,k3=1′l/(mol·s)

目的产物为B，反应进料为纯A，CA=5.8mol/L。由于反应的中间产品B即为反应的目的产物，故本例将反应分为一个阶段，即A→B和2A→D。其副产物瞬时选择性的导数如式(9)所示：

(9)

由导数分析结果可知van de Vusse反应的最优反应器网络为单个CSTR，对其总选择性进行优化并与PFR和CSTR+PFR结构的结果比较。各结构的优化结果列于表3，其中XA1和XA2为各反应器出口处A的总转化率。

从表3可知CSTR和CSTR+PFR结构的选择性均优于单个PFR，且前两者选择性基本相等，即仅CSTR就可使总选择性最大；图6是传统以CB最大为目标得到的CSTR+PFR结构在最优反应器网络中目的产物总选择性与转化率XA的关系图，可知总选择性在PFR中呈下降趋势，表明CSTR+PFR结构中的PFR对提高总选择性不利，故最优反应器网络为CSTR。

表3 不同反应器网络下van de Vusse反应优化结果Table 3 van de Vusse reaction optimization results under different reactor networks

图6 CSTR+PFR结构最优转化率下目的产物总选择性St～XA曲线Fig.6 St profiles under optimal conversion rate of CSTR+PFR structure

Van de Vusse反应是第二类反应体系中的典型反应，其结果证实本方法选择的过程参数和目标函数是准确的，且本方法同样适用于第二类反应体系。

4 结论

针对复杂反应体系中间产品较多的特点，发展了适用于复杂反应体系反应器网络综合的导数分析法。

(1) 本文提出将串联反应以反应物到目的产物的反应步数为依据的分段策略，通过对分段副产物瞬时选择性进行导数分析确定反应器类型，运用遗传算法进行求解确定最优的反应器网络。

(2) 采用三个算例验证了以分段副产物瞬时选择性的导数分析进行反应器网络综合的有效性和简捷性，且得到的最优反应器网络是由CSTR和PFR的串联而成。其中，van de Vusse反应的综合结果说明本方法同样适用于第二类反应体系。

符号说明：

Ci组分i的浓度/mol·L-1τ物料在反应器中的停留时间/sCA0 组分A的进料浓度/mol·L-1下标S 瞬时选择性i反应器入口k 反应速率常数exit反应器出口r反应速率/mol·L-1·s-1m组分mSt目的产物总选择性上标XA关键组分A的总转化率C 全混流反应器φ计量系数有关参数,本文取lP平推流反应器