乙烯齐聚化工流程熵产率最小化_于亚杰.pdf

第 39 卷第 1 期 化学反应工程与工艺 Vol 39,No 1 2023 年2 月 Chemical Reaction Engineering and Technology Feb.2023 收稿日期收稿日期：2021-11-01；修订日期修订日期：2022-12-10。作者简介作者简介：于亚杰（1997），男，硕士研究生；夏少军（1986），男，副教授，通讯联系人。E-mail:shaojunxia_。基金项目基金项目：国家自然科学基金（51976235；51606218）；湖北省自然科学基金（2018CFB708）。文章编号：文章编号：10017631(2023)01007209 DOI:10.11730/j.issn.1001-7631.2023.01.0072.09 乙烯齐聚乙烯齐聚化工流程熵产率化工流程熵产率最小化最小化 于亚杰，夏少军，赵 明 海军工程大学动力工程学院，湖北 武汉 430033 摘要：摘要：乙烯齐聚反应合成液体燃油技术在军用和民用领域具有广阔的应用前景。基于有限时间热力学，对包含混合器、压缩机、换热器以及乙烯齐聚反应器的化工流程进行了建模与优化，计算给定癸烯产率条件下流程熵产率最小时的压缩机出口压力、换热器出口温度以及换热器与反应器的管外热源分布，并将优化结果与参考化工流程以及反应器熵产率最小的对应流程进行了对比。结果表明，换热器出口温度为 608.18 K 和压缩机出口压力为 2.68 MPa 时，流程的熵产率最小，优化后的流程熵产率相比参考化工流程下降了 3.12%，相较于反应器熵产率最小对应流程熵产率降低了 1.30%。研究结果对乙烯齐聚化工流程最优设计具有一定的指导意义。关键词：关键词：有限时间热力学有限时间热力学 乙烯齐乙烯齐聚聚 熵产率熵产率 流程流程优化优化 中图分类号：中图分类号：TK123 文献标识码：文献标识码：A 海基合成燃油技术1可从海水中获取原料（H2和 CO2）并合成燃油。乙烯齐聚反应作为海基合成燃油“短链烃聚合”合成路径中的一环，可用于将乙烯合成为液体燃油原料。目前针对海基合成燃油技术的研究还不够充分，存在产率低、耗能高的问题，距离实际工程应用具有较大差距，需进一步进行反应性能分析与优化研究。许多学者已经针对不同类型烯烃的齐聚反应进行了经典热力学与动力学研究。吴建民等2计算了 298700 K 时乙烯齐聚制-烯烃反应的反应热、吉布斯自由能以及反应平衡常数，研究了温度和压力对反应化学平衡的影响。刘平等3分析对比了不同条件下催化异戊烯二聚反应的分子筛 Beta、ZSM-5、丝光沸石（MOR）和 ITQ-1 的性能。顾雪萍等4研究了一种以 BF3/C2H5OH为催化剂的 1-癸烯连续齐聚反应，建立了其反应动力学模型，并对模型准确性进行了验证。郑征等5建立了 Ziegler-Natta 催化作用下的丙烯/1-丁烯共聚反应动力学模型，并分析了其催化特性。袁明6优化了乙烯齐聚反应工艺条件，并建立反应动力学模型来表征实验结果。Toch 等7-9对乙烯齐聚反应进行了动力学分析，建立了一维非均相工业反应器模型，并分析了模型参数对反应器热效应以及产率的影响。上述研究皆考察了反应过程，通过改变反应器参数实现产物产率的增加，而在实际工程过程中，需要通过混合器、压缩机、换热器等配套过程设备以达到反应所需条件，因此需要对整个工艺流程进行建模与优化研究。化工流程是一个包含换热、流动、化学反应等多个过程的耦合过程，只有综合考虑多个过程，建立准确的数学模型，才能对流程进行更加精确的性能描述与优化。有限时间热力学实现了热力学、传热学、流体力学等多个学科的交叉，能更好地对这种复杂的耦合过程进行描述，以计算其目标函数极第39 卷第1 期 于亚杰等.乙烯齐聚化工流程熵产率最小化 73 值并得到到达极值的路径1。Rsjorde 等10研究了换热器与反应器组合的流程模型，以流程的总熵产率最小为优化目标，以反应器的进口温度和管外热源温度为变量进行了优化。Kingston 等11-12建立了包含恒温的单元反应器、不同绝热效率的压缩机和恒定温差的换热器的二甲醚合成流程模型，求解了流程的最小熵产率及相应的最优运行参数11，并对两个反应器的耦合进行了研究12。本研究将在文献1的基础上，应用有限时间热力学，建立包含混合器、压缩机、换热器与反应器的乙烯齐聚过程的化工流程模型，在压缩机效率给定、换热器与反应器管外热源温度完全可控以及产物产率给定的条件下，以流程熵产率最小为优化目标，对压缩机出口压力、换热器出口温度以及换热器和反应器的管外热源温度分布进行优化设计，分析各参数对流程熵产率的影响，并将优化结果与管外热源温度恒定的参考反应器以及反应器熵产率最小的对应流程相比较，为乙烯齐聚过程的实际流程设计提供科学依据。1 模型建立模型建立 工业上原料气体常以压缩状态储存于储气罐中，而气体从低压到高压状态需消耗能量。为全面描述乙烯齐聚化工流程性能，本研究将以环境温度（298.15 K）与环境压力（0.101 MPa）为初始态对原料气体进行处理。乙烯齐聚过程的流程如图 1 所示。气体组分为C2H4，C4H8，C10H20和 N2，通过彭-罗宾森方程13计算出混合气体在 570 K，2.5 MPa 与 710 K，4.5 MPa 下的压缩因子分别为 1.003 4 与 1.012 3，与理想气体偏差分别为 0.34%与 1.23%，故可以假设气体为理想气体。将反应原料 C2H4、C4H8以及保护性气体 N2通入混合器充分混合，经过压缩机与换热器使混合气体达到所设定的温度与压力，最终进入反应器进行反应。乙烯齐聚反应需要在高温高压催化剂存在的环境下才能进行，因而设定该流程中反应仅在反应器中发生。1.1 混合器混合器 气体的混合过程为等温等压混合，是不可逆过程，即存在熵产。混合过程的熵产率可根据进出口处熵的变化量（SM）计算得到：MlniiiSRny=（1）式中：R 为摩尔气体常数，J/(molK)；ni为组分 i 物质的量，mol；yi为组分 i 的摩尔分数。1.2 压缩机压缩机 压缩机的压缩过程为不可逆过程，存在熵产。假设压缩过程为不可逆绝热压缩且不发生化学反应。压缩机的效率（C）可用来体现压缩过程的不可逆程度：()()()()isenoutoutininCoutoutinin,h Tph Tph Tph Tp=（2）式中：pin和 pout为压缩机入口和出口压力，MPa；isenoutT为可逆过程下压缩机的出口温度，K；Tin 图 1 化工流程模型 Fig.1 Schematic diagram of chemical process model N2C4H8混合器压缩机换热器管式反应器C10H20C2H474 化学反应工程与工艺 2023年2月 和 Tout为压缩机入口和出口温度，K；h 为摩尔焓，J/mol。通过式（2）计算得到 Tout，压缩机出口与入口的熵产率（SC）可根据如下公式计算得到：outin5,C,in,in1.013 10dlnTp iTiTiTiCSFyTFyTpR=|（3）式中：FT,in为进入混合器的总流量，mol/s；Cp,i为组分 i 的定压摩尔热容，J/(molK)；p 为工质压力，Pa；T 为工质温度，K。等式右侧第二项为压力熵。1.3 换热器换热器 换热过程为不可逆过程，存在熵产。对换热器作如下假设：（1）反应气体在换热管内换热的模型为一维稳态模型；（2）换热过程不存在化学反应的进行。考虑换热过程服从线性唯象传热规律，则：()HEHEHEHE,a=11qUTT （4）式中：qHE为换热热流密度，J/(m2s)；UHE为换热器传热系数，(WK)/m2；THE为换热器管内混合气体温度；THE,a为换热器管外热源温度，K。则局部熵产率（HE）可以用如下公式来表示14-15：()HEHE,iHEHEHE,a11dqTT=（5）式中：dHE,i为换热管内径，m。换热过程的总熵产率（SHE）即为单位管长内的熵产率在轴向上的积分：HEHEHEHE0dLSz=（6）式中：LHE为换热器长度，m；zHE为离换热器进口处的轴向距离，m。1.4 反应器反应器 本研究采用多管固定床反应器，因各管束内反应管运行条件相同，故以单管为例进行计算。假设反应物在反应器内为一维活塞流稳态模型1，不存在径向温度与浓度梯度，不存在轴向的流体混合。反应器存在传热、黏性流动以及化学反应三个部分的熵产。反应物在反应器内主要发生以下反应：12448r11C HC H 02rH （7）2241020r21C HC H 05rH （8）反应器中各组分的质量守恒方程为：()24C HCRcb1 12 2ddFzArr=+（9）48C HCRcb 1 11dd2FzAr=（10）1020C HCRcb 2 21dd5FzAr=（11）第39 卷第1 期 于亚杰等.乙烯齐聚化工流程熵产率最小化 75 2NCRdd0Fz=（12）式中：Ac为反应器横截面积，m2；zCR为离反应管进口处的轴向距离，m；b为催化床层堆密度，kg/m3；j为反应 j 的内扩散有效因子；rj为反应 j 的反应速率，mol/(kgs)。反应速率可根据经验公式16以及阿伦尼乌斯公式计算得到：()244811C12H1C HrppKk=（13）()2410201 5C HH22C2pprKk=（14）()()1 000expjjjkAERT=（15）式中：pi为组分 i 的分压，MPa；kj为反应速率常数，mol1-n/(L3-3ns)；Aj为指前因子，mol/(skgMPa)；Kj为反应平衡常数，molm/Ln；m 与 n 分别为生成物与反应物的化学反应计量系数；Ej为反应活化能，J/mol。Aj，Ej和 Kj的取值参见文献1。传热和化学反应过程遵循如下能量守恒方程：()CR,iCRCRCR,acbrCR,11ddj jjjkp kkdUTTArHTzF C=（16）式中：dCR,i为反应器内径，m；TCR和 TCR,a分别为反应器管内反应物温度与管外热源温度，K；UCR为总传热系数，(WK)/m2；rHj为反应焓，J/mol。黏性流动过程遵循如下动量守恒方程16：()2mm233CRpp1150d1.751dpGvzdd|=+|（17）式中：为催化床层空隙率；m为反应物黏性系数，kg/(ms)；dp为催化剂颗粒直径，m；G 为物料质量流率，kg/(sm2)；vm为物料平均流速，m/s。根据非平衡热力学理论，反应器的局部熵产率（CR）为14-15：2rCRCR,iCRc mcbCRCR,aCRCRCR111ddjj jjGpdUA vArTTTzT=+|（18）总熵产率（SCR）即为局部熵产率在轴向上的积分：CRCRCRCR0dLSz=（19）式中：LCR为反应器长度，m。1.5 目标函数目标函数 乙烯齐聚化工流程优化，即以流程总熵产率最小为优化目标，得到压缩机出口压力、换热器出口温度以及换热器与反应器最优管外温度分布，为最优构型问题。流程总熵产率（Stot）为：totMCHECRSSSSS=+（20）2 流程熵产率最小化求解流程熵产率最小化求解 求解管外热源温度完全可控的烯烃齐聚子流程熵产率最小最优构型问题，即计算给定约束下的烯76 化学反应工程与工艺 2023年2月 烃齐聚子流程最小熵产率以及对应的换热器和反应器热源温度分布，其中，该问题的约束条件为子流程产率一定以及各部件单元所满足的守恒方程。为研究压缩机出口压力与换热器出口温度对优化结果的影响，设定了不同边界条件进行对比分析，在各节点参数与参考流程相同的条件下，对换热器与反应器管外热源进行了优化设计。（1）对换热器与反应器管外热源恒定且各优化参数给定的参考流程及其性能进行求解。其中，给定换热器与反应器管外热源温度恒定为 637 K，压缩机出口压力为 3 MPa，该类情形用 Reference Process 表示。（2）在保证反应器进口温度与进口压力、C10H20产率以及流程入口摩尔流率与 Reference Process相同的情况下，将换热器与反应器的管外热源温度设为控制变量对乙烯齐聚子流程进行优化，该类情形用 Case 1 表示。（3）在保证 C10H20产率以及流程入口摩尔流率与 Reference Process 相同的情况下，以换热器与反应器的管外热源温度为控制变量，在反应器进口温度自由的条件下对乙烯齐聚子流程进行优化，该类情形用 Case 2 表示。（4）在保证 C10H20产率以及流程入口摩尔流率与 Reference Process 相同的情况下，以换热器与反应器的管外热源温度为控制变量，在反应器进口温度和进口压力均自由的条件下对乙烯齐聚子流程进行优化，该类情形用 Case 3 表示。参考流程各单元设备参数如表 1表 3 所示。表表 1 混合器参数混合器参数 Table 1 Mixer parameters Parameters TM,in/K PM,in/MPa 24C H,iny 48C H,iny 1020C H,iny 2N,iny FT,in/(mols-1)Values 298 0.101 0.399 0.100 0.001 0.500 1 表表 2 换热器参数换热器参数 Table 2 Heat exchanger parameters Parameters UHE/(WKm-2)dHE,i/m LHE/m Values 1.7107 0.08 5 表表 3 反应器参数反应器参数 Table 3 Reactor parameters Parameters UCR/(WKm-2)dCR,i/m dCR,o/m LCR/m b/(kgm-3)TCR,a/K dp/m Values 3107 0.08 0.084 5 800 637 0.005 0.45 3 优化结果及分析优化结果及分析 考虑反应过程工质温度变化较大，为充分考虑传热过程熵产率的减小，减小工质与热源间温差，设定换热器与反应器管外热源温度为 4001 000 K。在文献1的基础上，为保证最优反应器进口温度在优化范围内而非优化范围边界值，设定较大的反应器进口温度优化，为 570710 K，反应器进口压力为 2.54.5 MPa。在保证产物（C10H20）产率为 0.048 4 mol/s 的基础上，通过 Matlab 求解器“bvp4c”对反应器和换热器的最优构型、压缩机出口压力以及换热器出口温度进行了求解。3.1 不同边界条件不同边界条件下下熵产率熵产率对比分析对比分析 表 4 为不同边界条件下的熵产率以及压缩机出口压力与换热器进口温度取值。由表 4 可见，流程第39 卷第1 期 于亚杰等.乙烯齐聚化工流程熵产率最小化 77 的熵产率以混合器、压缩机和反应器部分为主，而换热器部分的熵产率相较于其他部分较小。Case 1相较于参考流程优化了换热器和反应器的管外热源分布，降低了换热器与反应器的熵产率，流程的总熵产率减小，其中，换热器的熵产率占比较小，优化以降低反应器的熵产率为主。Case 2 相较于 Case 1进一步降低了换热器的出口温度，即反应器的入口温度。换热器出口温度越接近入口温度，换热过程的换热量越少，换热器熵产率越小，因而流程总熵产率降低。Case 3 相较于 Case 2 优化了压缩机出口压力，即反应器入口压力。压缩机出口压力越小，压缩机熵产率越小，但反应器入口压力减小会导致反应器熵产率增大，因此 Case 3 对压缩机与换热器的熵产率进行了权衡，在增大反应器熵产率的同时降低压缩机的熵产率，最终使得流程总熵产率降低。表表 4 不同条件下不同条件下化工流程优化结果比较化工流程优化结果比较 Table 4 Comparisons of optimization results of the chemical process at different conditions Parameters Reference process Case 1 Case 2 Case 3 pout/MPa 3.00 3.00 3.00 2.68 THE,in/K 637.00 637.00 611.01 608.18 SM/(WK-1)7.901 6 7.901 6 7.901 6 7.901 6 SC/(WK-1)8.203 0 8.203 0 8.203 0 7.963 9 SHE/(WK-1)0.149 5 0.140 2 0.001 6 0.028 6 SCR/(WK-1)6.472 5 6.113 8 5.977 4 6.123 6 Stot/(WK-1)22.726 6 22.358 6 22.083 6 22.017 7 Reduce scale,%1.62 2.83 3.12 3.2 温度构型对比分析温度构型对比分析 换热器的熵产率为传热过程的熵产率。反应物在换热器中温度变化越小，熵产率越小。不同条件下最优流程的换热器管外热源温度分布以及反应器管外热源温度分布见图 2 和图 3。结果表明：随着优化的不断深入，换热器出口温度越来越小，换热器所需管外热源温度也不断降低，而 Case 3 压缩机出口压力的减小使得压缩机出口温度降低，换热器进出口温度差增大，换热过程熵产率小幅度升高。012345600620640660680 THE,a/KzHE/m reference process Case 1 Case 2 Case 3 012345600640680720760800 TCR,a/KzCR/m reference process Case 1 Case 2 Case 3 图 2 换热器管外热源温度对比 Fig.2 Comparison of heat source temperature outside the heat exchanger tube 图 3 反应器管外热源温度对比 Fig.3 Comparison of heat source temperature outside the reactor tube 图 4 为不同条件下反应物在换热器和反应器内的温度曲线。反应器的熵产率由传热过程、黏性流动以及化学反应三个部分组成，其中黏性流动的熵产率远小于其他部分的，可忽略不计。由于乙烯齐聚反应的两个主要反应为放热反应，当入口处反应物温度与管外热源温度差增大时，化学反应的反应速率降低，化学反应过程熵产率减小，但传热过程的熵产率增大；当温差减小时，放热反应的反应速 78 化学反应工程与工艺 2023年2月 率增大，但传热过程熵产率降低。因此，反应器最优构型是综合决策传热过程与化学反应两部分熵产率得到的。由图 3 可以看出，优化后的反应器管外热源温度为先降低再升高。在实际工程过程中，对于反应器管外热源降低段的反应器热源温度可考虑让工质通过顺流反应器实现；对于温度上升段，可大致分为斜率不同的三段，可考虑让工质通过三个不同的逆流反应器实现，如图 5 所示，其中逆流反应器 2热源工质流率逆流反应器 1 热源工质流率逆流反应器 3 热源工质流率。co-current reactor countercurrent reactor 1 countercurrent reactor 2 countercurrent reactor 3 图 5 最优反应路径示意 Fig.5 Diagram of optimal reaction paths 3.3 最优流程与反应器熵产率最小对应流程最优流程与反应器熵产率最小对应流程的的对比对比 优化流程 Case 3 与反应器熵产率最小的对应流程性能比较见表 5。由表 5 可知，反应器熵产率最小对应流程虽然使得反应器熵产率相较于 Case 3 得到优化，却是要求较高的压力条件，使得流程中压缩机熵产率变大，最终流程总熵产率大于最优流程熵产率。最优流程熵产率相较于反应器熵产率最小对应流程的熵产率降低了 1.30%，实际工程中采用多根反应管并联，使得总能耗较高，该优化得到的收益就尤为可观。因此，相比于单个反应器优化，考虑化工流程全局优化更具科学性和指导意义。表表 5 优化流程优化流程 Case 3 与反应器熵产率最小对应流程性能比较与反应器熵产率最小对应流程性能比较 Table 5 Performance comparisons between the optimized process for the Case 3 and the process with the minimum entropy generation rate of the reaction Parameters Case 3 Reactor optimization pout/MPa 2.68 3.34 THE,in/K 608.18 605.35 SM/(WK-1)7.901 6 7.901 6 SC/(WK-1)7.963 9 8.430 0 SHE/(WK-1)0.028 6 0.036 3 SCR/(WK-1)6.123 6 5.940 2 Stot/(WK-1)22.017 7 22.308 1 4 结结 论论 本研究应用有限时间热力学理论，建立了包含混合器、压缩机、换热器以及反应器的乙烯齐聚化工流程模型，并研究了在给定压缩机效率、换热器与反应器管外温度完全可控且产物产率一定的条件下流程熵产率最小化，计算了包括压缩机出口压力、换热器出口温度、换热器管外热源温度分布以及TR1,aT(z)co-current reactorTR2,aT(z)Countercurrent Reactor 1TR3,aT(z)Countercurrent Reactor 2TR4,aT(z)Countercurrent Reactor 30246810600640680720760800zCRzHE T/Kz/m reference process Case 1 Case 2 Case 3 图 4 反应物温度对比 Fig.4 Reactant temperature comparison 第39 卷第1 期 于亚杰等.乙烯齐聚化工流程熵产率最小化 79 反应器管外热源温度分布的最优构型。结果表明：a）与参考流程相比，对换热器与反应器管外热源分布进行优化使流程熵产率降低了 1.62%；若仅释放换热器出口温度约束，优化后的流程熵产率相较于参考流程熵产率降低了 2.83%；若同时释放换热器出口温度以及压缩机出口压力约束，优化后的流程熵产率相较于参考流程熵产率降低了3.12%。b）在本研究参数计算范围内，反应器熵产率最小对应流程的熵产率虽小于参考流程熵产率，但大于相同约束条件下的优化整体流程的流程熵产率；相较于优化单个反应器，优化整体流程对反应流程最优设计与运用更具科学性与指导意义。符号说明符号说明 A 指前因子，mol/(skgMPa)T 温度，K Ac 横截面积，m2 U 传热系数，WK/m2 Cp,i 组分 i 的定压摩尔热容，J/(molK)v 流速，m/s d 管内径，m yi 组分 i 的摩尔分数 do 管外径，m z 轴向坐标，m dp 催化剂颗粒直径，m 催化床层空隙率 E 反应活化能，J/mol 效率 FT,in 进入混合器的物质总流量，mol/s j 反应 j 的内扩散有效因子 G 质量流率，kg/(sm2)黏性系数，kg/(ms)h 摩尔焓，J/mol b 催化床层堆密度，kg/m3 K 平衡常数，molm/Ln 局部熵产率，W/(Km)k 化学反应速率常数，mol1-n/(L3-3ns)S 熵，J/(molK)L 管长，m kS 标准摩尔熵，J/(molK)n 反应物的化学反应系数 S 熵产率，W/K m 生成物的化学反应系数 rH 反应焓，J/mol ni 组分 i 物质的量，mol 下标 p 压力，MPa C 压缩机 pi 组分 i 的分压，MPa CR 反应器 q 热流密度，W/m2 HE 换热器 R 普适气体常数，J/(molK)M 混合器 r 化学反应速率，mol/(kgs)参考文献：参考文献：1 陈林根,王超,张磊,等.海基合成燃油系统关键单元的热力学分析与优化研究进展J.中国科学:技术科学,2021,51(2):137-175.CHEN Lingen,WANG Chao,ZHANG Lei,et al.Progress in thermodynamic analyses and optimizations for key component units in sea-based fuel synthesis systemJ.Scientia Sinica:Technologica,2021,51(2):137-175.2 吴建民,孙启文,张宗森,等.乙烯齐聚制-烯烃反应体系的热力学研究 J.分子催化,2013,27(4):385-391.WU Jianmin,SUN Qiwen,ZHANG Zongsen,et al.Thermodynamic study for ethylene oligomerization to-olefinJ.Journal of Molecular Catalysis,2013,27(4):385-391.3 刘平,任忠飞,周磊,等.ITQ-1 分子筛催化异戊烯高效二聚反应J.化学反应工程与工艺,2021,37(1):24-29,72.LIU Ping,REN Zhongfei,ZHOU Lei,et al.High efficient dimerization of isoamylene over ITQ-1 zeoliteJ.Chemical Reaction Engineering and Technology,2021,37(1):24-29,72.4 顾雪萍,袁诗谣,冯黄河,等.面向相对分子质量分布的 1-癸烯连续齐聚过程动力学模型J.化学反应工程与工艺,2021,37(1):30-38.80 化学反应工程与工艺 2023年2月 GU Xueping,YUAN Shiyao,FENG Huanghe,et al.Relative molecular mass distribution oriented kinetic model of 1-decene continuous oligomerization processJ.Chemical Reaction Engineering and Technology,2021,37(1):30-38.5 郑征,杨以琳,黄凯,等.Ziegler-Natta 催化剂催化的丙烯/1-丁烯共聚:催化特性和动力学模型J.高校化学工程学报,2019,33(2):355-363.ZHENG Zheng,YANG Yilin,HUANG Kai,et al.Copolymerization of propene and 1-butene using a Ziegler-Natta catalyst:catalytic characteristic and kinetic modelsJ.Journal of Chemical Engineering of Chinese Universities,2019,33(2):355-363.6 袁明.甲醇合成汽油工艺中 ZSM-5 催化乙烯齐聚关键步骤反应条件的优化研究D.武汉:武汉工程大学,2017.7 TOCH K,THYBAUT J W,ARRIBAS M,et al.Ethylene oligomerization catalyst optimization using fundamental kinetic modelingC.2011 AIChE Annual Meeting,Minneapoils,MN,2011:16-21.8 TOCH K,THYBAUT J W,MARIN G B.Ethene oligomerization on Ni-SiO2-Al2O3:experimental investigation and single-event microkinetic modelingJ.Applied Catalysis A:General,2015,489:292-304.9 TOCH K,THYBAUT J W,MARIN G B.An industrial reactor design for heterogeneously catalyzed ethene oligomerizationC.2013 AIChE Annual Meeting,San Francisco,CA,2013:3-8.10 RSJORDE A,KJELSTRUP S,JOHANNESSEN E,et al.Minimizing the entropy production in a chemical process for dehydrogenation of propaneJ.Energy,2007,32(4):335-343.11 KINGSTON D,RAZZITTE A.Entropy production in chemical reactorsJ.Journal of Non-Equilibrium Thermodynamics,2017,42(3):265-275.12 KINGSTON D,RAZZITTE A.Entropy generation minimization in dimethyl ether synthesis:a case studyJ.Journal of Non-Equilibrium Thermodynamics,2018,43(2):111-120.13 陈宏芳,杜建华.高等工程热力学M.北京:清华大学出版社,2002:70-85.14 KJELSTRUP S,BEDEAUX D,JOHANNESSEN E,et al.Non-equilibrium thermodynamics for engineersM.New Jersey:World Scientific,2010:235-246.15 BEDEAUX D,KJELSTRUP S,SENGERS J V.Experimental thermodynamics.Volume X:non-equilibrium thermodynamics with applications.cambridge M.UK:The Royal Society of Chemistry,2016:61-77.16 FROMENT G F,BISCHOFF K B.Chemical reactor analysis and designM.New York:Willey,1990:525-532.Minimization of Entropy Generation Rate in the Chemical Process of Olefins Oligomerization YU Yajie,XIA Shaojun,ZHAO Ming College of Power Engineering,Naval University of Engineering,Wuhan 430033,China Abstract:The technology of olefin oligomerization to synthesize liquid fuel has broad application prospects in military and civil fields.Based on the theory of finite-time thermodynamics,the chemical process model including a mixer,a compressor,a heat exchanger and an ethylene oligomerization reactor was established and optimized.The compressor outlet pressure,the heat exchanger outlet temperature and the optimal temperature configurations of the heat reservoirs outside the heat exchanger and the reactor for the minimum entropy generation rate of the chemical process were derived under the condition of a given yield of C10H20.The optimization results were also compared with those for the reference chemical process and the process with the minimum entropy generation rate of the reactor.The results showed that the entropy generation rate for the optimal process was reduced by 3.21%compared to that for the reference chemical process and by 1.30%compared to that for the process with the minimum entropy generation of the reactor when the heat exchanger outlet temperature was 608.18 K and compressor outlet pressure was 2.68 MPa.The obtained results in this paper could be a certain theoretical guidance for the optimal design of olefin oligomerization reaction process.Key words:finite-time thermodynamics;ethylene oligomerization;entropy generation rate;process optimization 栏目编辑 蒋斌波