反应精馏塔是将反应过程和精馏过程综合在一起的复杂化工装置。由于能使可逆反应向目标产物的方向进行,反应精馏塔极大程度地提高了反应转化率,节约了能源消耗与设备投资成本,已经成为现代化工过程研究与开发的一个重要研究对象。针对包含两步连续可逆反应的复杂混合物的分离问题,单反应段精馏塔(RDC-SRS)由于难以充分协调两步反应之间的相互耦合,常常导致非常高的能源消耗和设备投资成本。双反应段精馏塔(RDC-DRS)具有较多的自由度,能够有效协调两步反应之间的相互耦合,因而在节省能源消耗和设备投资成本方面具有较大的潜力。Yu等[1]首先对RDC-DRS进行了较为系统的研究,主要对理想两步连续可逆反应的分离和碳酸二甲酯与乙醇两段连续酯交换反应的分离进行了RDC-DRS的设计,并与RDC-SRS进行了比较,结果表明在节约能源和设备投资费用方面RDC-DRS具有良好的效果。夏春英等[2]以合成碳酸二乙酯为例对RDC-DRS的闭环浓度控制进行了研究,所得结果显示RDC-DRS具有良好的抗扰动性能。然而迄今为止,对RDC-DRS温度控制的研究还寥寥无几。Kaymak等[3-4]基于灵敏度分析方法对RDC-DRS的温度控制进行了较为系统的研究,他们发现灵敏度分析方法虽然有时能够实现对RDC-DRS产品质量的严格控制,但有时却表现出极其宽松的控制效果。这种现象不仅充分表明RDC-DRS具有非常复杂的稳态与动态特性,还体现了研究RDC-DRS温度控制的必要性。
本文主要对RDC-DRS温度控制系统的设计方法进行研究。基于稳态偏差最小原理,提出一种有效的温度控制系统设计方法,并以理想两步连续可逆反应的分离为例进行了RDC-DRS温度控制系统的设计与分析。通过与灵敏度分析方法设计的温度控制系统的控制效果相比较,表明了本文提出的温度控制系统综合与设计方法的有效性与可行性。
1 RDC-DRS的稳态设计本文在Yu等[1]所建立的RDC-DRS稳态模型的基础上进行温度控制系统的研究。该反应物系包含A、B、C、D、E 5种物质,其相对挥发度的比值为αD:αB:αC:αA:αE = 16:8:4:2:1。反应动力学与热力学模型如式(1)、(2)所示
$ {\rm{A}} + {\rm{B}} \leftrightarrow {\rm{C}} + {\rm{D}}\;\;\Delta {H_{\rm{r}}} = 0\;{\rm{kJ}}/{\rm{kmol}} $ | (1) |
$ {\rm{C}} + {\rm{B}} \leftrightarrow {\rm{E}} + {\rm{D}}\;\;\Delta {H_{\rm{r}}} = 0\;{\rm{kJ}}/{\rm{kmol}} $ | (2) |
两个反应段的反应速率为
$ \begin{array}{*{20}{l}} {{R_{1,i}} = {\nu _{1,i}}\;H({k_{ + 1}}{x_{\rm{A}}}{x_{\rm{B}}} - {k_{ - 1}}{x_{\rm{C}}}{x_{\rm{D}}})\;\;\;i = {\rm{A}},\;{\rm{B}},}\\ {{\rm{C}},\;{\rm{D}}} \end{array} $ | (3) |
$ \begin{array}{*{20}{l}} {{R_{2,i}} = {\nu _{2,i}}\;H({k_{ + 2}}{x_{\rm{C}}}{x_{\rm{B}}} - {k_{ - 2}}{x_{\rm{E}}}{x_{\rm{D}}})\;\;\;i = {\rm{B}},{\rm{C}},}\\ {{\rm{D}},\;{\rm{E}}} \end{array} $ | (4) |
式(3)、(4)中,ν1、ν2分别为第一步、第二步反应计量比;H为催化剂用量;正向反应速率k+m=a+ me-E+m/RT,反向反应速率k-m=a-me-E-m/RT,其中,m=1、2,a为指前系数,E为活化能,R为理想气体常数,T为温度。
RDC-DRS的结构如图 1所示。较重的反应物A从上反应段的顶部进料,而较轻的反应物B分别从上、下两个反应段的底部进料。最重的产物E从塔底出料,最轻的产物D从塔顶出料,产物出料浓度(摩尔分数)均设定为0.95。其中,RR表示回流量,Qreb表示再沸器热负荷。
目前已经有很多灵敏板的选择方法用于精馏塔的控制,但灵敏度分析方法仍是应用最为广泛的一种[5-8]。它是以操作变量所引起的塔板温度的最大变化为判断依据,具有原理简单、计算强度小等特点。图 2给出了RDC-DRS的灵敏度分析结果,其中4个操作变量的变化幅度均为-1‰。
RDC-DRS具有4个操作变量,即回流量(RR)、再沸器热负荷(Qreb)和两个B组分进料流量(FB1,FB2)。4个操作变量所对应的最大温度变化均发生在第26块板。另外,回流量和两个B组分进料流量的次最大温度变化均发生在第20块板,其再次最大温度变化发生在第10块板。基于就近配对原则,温度控制系统的结构应为RR控制第10块板的温度,Qreb控制第26块板的温度,FB1控制第20块板的温度。对于FB2控制回路,选择最近的第24块板作为被控塔板。虽然就近选取第24块板不一定是最佳的选择,但它并不会显著影响温度控制系统的动态品质。FB2回路与Qreb回路所控塔板距离很近,导致二者之间存在较强的相互作用,但这有助于比较和分析灵敏度分析方法与本文所提方法的性能。由此得到如图 3所示的温度控制系统(CS1),其中两个液位控制器用来调节回流罐和再沸器的液位,一个压力控制器用来调节塔压,一个流量控制器用来控制A进料的进料量FA。
虽然温度控制系统的性能应根据被控产品质量的瞬态偏差和稳态偏差来评估,但在选择温度灵敏板时如果同时考虑这两种偏差,温度控制系统的综合与设计问题就成为一个混合整数非线性动态优化问题。其求解不但原理非常复杂,而且计算强度大,极不适应于控制结构与控制器参数的有效搜索。由于在严格的产品质量控制条件下较小的稳态偏差对应着较小的瞬态偏差,所以以稳态偏差最小为原则来搜索灵敏板位置,这样温度控制系统的综合与设计问题就成为一个较简单的混合整数非线性稳态优化问题,既可以采用单变量优化方法进行求解,又能大大降低计算强度。
灵敏板搜索采用如式(5)所示的目标函数
$ J({L_1},{\rm{ }}{L_2},{\rm{ }}{L_3},{\rm{ }}{L_4}) = \sum\limits_i {(|\Delta {x_{\rm{D}}}{|_i} + |\Delta {x_{\rm{B}}}{|_i})} $ | (5) |
式中,i表示FA+20%、FA-20%以及在FA中添加5% B的3种扰动情况;L1、L2、L3、L4分别表示RR、FB1、FB2、Qreb所对应的灵敏板的位置;|ΔxD|和|ΔxB|分别表示塔顶和塔底产品出料浓度的稳态偏差值。
为了简化计算,采用单变量优化方法对上述问题进行求解。具体操作步骤如下:
(1) 以CS1为基础,只改变RR所对应的被控塔板L1的位置,保持其他被控塔板不变,搜索使目标函数J最小的塔板并作为新的灵敏板L1;
(2) 在步骤(1)的基础上,再改变FB1所对应被控塔板L2的位置,搜索使J最小的塔板并作为新的灵敏板L2;
(3) L3,L4的搜索过程同上;
(4) 在完成一轮搜索后,再从头进行新一轮搜索,直至目标函数J达到最小值。
具体搜索结果如表 1所示。在搜索过程中,L1、L2、L3、L4搜索的先后顺序对得到的最终结果没有影响。通过上述搜索流程,得到如图 4所示的温度控制系统(CS2)。其中,RR被用来控制第2块板的温度,Qreb被用来控制第29块板的温度,FB1被用来控制第11块板的温度以及FB2被用来控制第26块板的温度。不难看出,CS2具有完全不同于CS1的灵敏板位置,这是由RDC-DRS非常复杂的稳态与动态特性所决定的。虽然利用回流量与再沸器分别控制最顶端与最底端塔板的温度有时会因为远离进料位置而一定程度地恶化控制系统的动态品质,但通过其他被控塔板的有效选择,本文所提方法仍能够使得被控塔板温度之间具有较弱的相互作用,且与产品浓度之间保持较紧密的对应关系,并最终改善温度控制的效果。
RDC-DRS的稳态与动态模型均采用平衡级模型,包括物料衡算方程、能量衡算方程、归一化方程、相平衡方程以及反应动力学方程,本文利用Mathematica对其进行求解。扰动发生后,两个温度控制系统均可以在较短的时间内使被控温度达到设定值。由于本文的主要评价指标是产品浓度,因此在这里对温度控制的效果不再赘述。
图 5展现了A进料流量发生±20%扰动后,CS1和CS2的动态响应。图 6给出了在A进料中添加5% B组分的扰动后CS1和CS2的动态响应。从图 6中可以明显地看出,CS2具有与CS1相近的调节时间,但最大偏差很小。对于CS1,A进料流量发生±20%变化所导致的塔顶产品D的出料浓度分别为0.941 7和0.958,塔底产品E的出料浓度分别为0.952和0.947 4,在FA中添加5% B组分的扰动所导致的塔顶、塔底出料浓度分别为0.939 8和0.949 6。对于CS2,A进料流量发生±20%变化所导致的塔顶产品D的出料浓度分别为0.949 5和0.950 7,塔底产品E的出料浓度分别为0.949 8和0.950 1,在FA中添加5% B组分的扰动所导致的塔顶、塔底出料浓度分别为0.949 6和0.949 7。相较于CS1,CS2具有更好的控制效果,其产品质量远比CS1的更接近于0.95。CS1的控制效果与Kaymak等[3]设计的温度控制系统的控制效果非常相似。这也说明了在设计CS1时采用FB2控制第24块塔板温度不是影响控制系统性能的主要因素。CS1,CS2具体的稳态偏差数据见表 2。可以看出,CS2的偏差值均比CS1的小了1到2个数量级,且与设定值0.95极其接近,改进显著。这些结果有力地表明了本文所提的温度控制系统综合与设计方法的有效性与可行性。
虽然双反应段精馏塔复杂的拓扑结构以及两步连续可逆反应的内在关联使得基于灵敏度分析方法的温度控制系统综合与设计难以得到良好的效果,但这并不意味温度控制系统不能用于RDC-DRS的操作。本文提出的基于稳态偏差最小原理的温度控制系统的综合与设计方法能够体现温度与产品浓度之间的对应关系,并且同时对4个被控塔板进行搜索可以充分调节4个变量之间的相互作用。CS1和CS2的闭环响应比较结果显示,CS2可以对产品浓度进行更为严格的控制,证明本文所提方法有效且可行。
[1] |
YU C H, YAO X H, HUANG K J, et al. A reactive distillation column with double reactive sections for the separations of two-stage consecutive reversible reactions[J]. Chemical Engineering and Processing:Process Intensification, 2014, 79: 56-68. DOI:10.1016/j.cep.2014.03.003 |
[2] |
夏春英, 黄克谨, 王韶锋. 应用于两步连续可逆反应的新型反应精馏塔的控制[J]. 现代化工, 2015, 35(8): 173-176. XIA C Y, HUANG K J, WANG S F. Control of a novel reactive distillation column for the two-stage consecutive reversible reactions[J]. Modern Chemical Industry, 2015, 35(8): 173-176. (in Chinese) |
[3] |
KAYMAK D B, VNLV H, ÖFKELI T. Control of a reactive distillation column with double reactive sections for two-stage consecutive reactions[J]. Chemical Engineering and Processing:Process Intensification, 2017, 113: 86-93. DOI:10.1016/j.cep.2016.09.010 |
[4] |
OKSAL I N, KAYMAK D B. Dynamic controllability comparison of reactive distillation columns with single and double reactive sections for two-stage consecutive reactions[J]. Chemical Engineering Research and Design, 2018, 129: 391-402. DOI:10.1016/j.cherd.2017.11.025 |
[5] |
LUYBEN W L. Evaluation of criteria for selecting temperature control trays in distillation columns[J]. Journal of Process Control, 2006, 16(2): 115-134. DOI:10.1016/j.jprocont.2005.05.004 |
[6] |
RADEMAKER O, RIJNSDORP J E, MAARLEVELD A. Dynamics and control of continuous distillation units[M]. Amsterdam: Elsevier, 1975.
|
[7] |
HEATH J A, KOOKOS I K, PERKINS J D. Process control structure selection based on economics[J]. AIChE Journal, 2000, 46(10): 1998-2016. DOI:10.1002/(ISSN)1547-5905 |
[8] |
KOOKOS I K. Real-time regulatory control structure selection based on economics[J]. Industrial & Engineering Chemistry Research, 2005, 44(11): 3993-4000. |