李文艷，顧雪萍，王嘉駿，馮連芳

(浙江大學 化學工程與生物工程學系 化學工程聯(lián)合國家重點實驗室，浙江 杭州 310027)

工業(yè)技術

規(guī)?；蹖Ρ蕉姿嵋叶减ソK縮聚的流程模擬

李文艷，顧雪萍，王嘉駿，馮連芳

(浙江大學 化學工程與生物工程學系 化學工程聯(lián)合國家重點實驗室，浙江 杭州 310027)

歸納了目前已報道的關于聚對苯二甲酸乙二醇酯(PET)終縮聚的動力學機理及參數(shù)、傳質(zhì)機理及參數(shù)，建立了相應的數(shù)學模型;考察了縮聚、醇解、酯化、水解、熱裂解、二甘醇和乙醛生成反應的速率常數(shù)和傳質(zhì)因子對PET終縮聚的影響;采用工廠的實際數(shù)據(jù)對模型進行了修正。結果表明，縮聚、醇解、熱裂解反應的速率常數(shù)對二甘醇、乙醛和PET的質(zhì)量流量及PET的相對分子質(zhì)量影響較大;熱裂解反應的速率常數(shù)、酯化反應速率常數(shù)對對苯二甲酸的質(zhì)量流量影響較大;傳質(zhì)因子對PET的數(shù)均相對分子質(zhì)量(Mn)影響較大。根據(jù)修正后的PET終縮聚的機理模型，得到羧基轉化率、PET的Mn、端羧基值等模型數(shù)據(jù)與實測數(shù)據(jù)相對偏差在2%以內(nèi)。

聚對苯二甲酸乙二醇酯;動力學;傳質(zhì);流程模擬;縮聚

聚對苯二甲酸乙二醇酯(PET)是重要的纖維原料和工程塑料，在塑料領域中占有重要地位。終縮聚是PET生產(chǎn)的關鍵環(huán)節(jié)之一，其模型化對穩(wěn)定生產(chǎn)過程、優(yōu)化工藝條件等均具有重要的指導意義。

PET終縮聚時，其反應體系的物性計算、動力學參數(shù)的修正和傳質(zhì)模型的正確建立是實現(xiàn)PET工業(yè)裝置模型化的3個重要部分。Ravindranath等［1－2］建立了終縮聚的二階數(shù)學模型、PET相對分子質(zhì)量分布模型及傳質(zhì)過程的數(shù)學模型［3］，得出了影響PET生產(chǎn)的主要因素，但模擬的結果只體現(xiàn)了各種因素對聚合度的影響。Laubriet等［4］應用兩相模型模擬了穩(wěn)態(tài)下聚合條件對PET產(chǎn)品質(zhì)量指標的影響，但模擬結果未與實驗數(shù)據(jù)或工廠數(shù)據(jù)進行對比。Yoon等［5］認為縮聚反應發(fā)生在膜相，縮聚反應受小分子擴散速率的控制，但模型參數(shù)的數(shù)據(jù)僅來源于實驗室。趙玲等［6］通過靜止膜實驗研究了PET的縮聚規(guī)律，并考察了在PET縮聚過程中不同厚度的靜止膜對傳質(zhì)的影響［7］，提出了傳質(zhì)速率表達式，但未考慮終縮聚的副反應?？梢?，目前PET終縮聚的數(shù)學模型大多基于實驗室規(guī)?；蛑性囈?guī)模。基于反應機理和規(guī)?；b置的PET終縮聚的數(shù)學模型更適用于模擬工業(yè)反應器。

本工作歸納了已報道的PET終縮聚的動力學機理及參數(shù)、傳質(zhì)機理及參數(shù)，考慮了主副反應和傳質(zhì)對PET終縮聚的影響，建立了相應的數(shù)學模型，并采用工業(yè)裝置的實際數(shù)據(jù)對模型進行了驗證。

1 PET終縮聚工藝與反應機理

1.1 PET終縮聚工藝

圓盤反應器是PET生產(chǎn)中常用的終縮聚釜。特性黏數(shù)(［η］)為0.17～0.19 dL/g的PET預聚物進入圓盤反應器后以近似平推流的方式流動，通過抽真空的方法將反應過程中生成的乙二醇、少量水、少量的乙醛和低聚合度PET形成的氣態(tài)混合物移出反應器。終縮聚的反應條件及結果見表1。

1.2 PET終縮聚反應機理

PET終縮聚反應主要包括縮聚/醇解和酯化/水解反應，同時還存在生成乙烯基的熱裂解及生成乙醛和二甘醇的熱降解副反應，具體反應見表2。

表1 PET終縮聚的反應條件(100 t/d)Table1 Polycondensation conditions for polyethylene terephthalate(PET)(100 t/d)

表2 PET終縮聚階段的反應類型Table 2 Reaction types for PET polycondensation process

2 PET終縮聚的傳質(zhì)模型

PET終縮聚的主反應(縮聚反應)是可逆的，小分子產(chǎn)物的移除有利于反應的正向進行，因此對小分子傳質(zhì)方面的研究有助于PET產(chǎn)品質(zhì)量的控制和反應器的設計。

圓盤反應器中的物料分為兩相，即平推流的熔融相和完全混合的氣相，采用文獻［4］報道的傳質(zhì)模型模擬PET終縮聚的傳質(zhì)過程，將傳質(zhì)阻力因素歸屬于一個傳質(zhì)系數(shù)kLa。根據(jù)膜理論，在平衡物料中對揮發(fā)性物料乙二醇、水和二甘醇僅使用同一個kLa，用線性濃度項的形式加入到揮發(fā)性的平衡物料中。式(1)是揮發(fā)性物料的傳質(zhì)通量(J)的表達式：

kLa的關聯(lián)式為：

式中，A為模型中設定的參考傳質(zhì)因子，是可調(diào)參數(shù)，A確定于特定的工業(yè)裝置;Mn為模型計算所得的PET數(shù)均相對分子質(zhì)量;ω為終縮聚釜攪拌轉速。

3 靈敏度分析

文獻［8］報道的反應條件和催化體系與本工作的建模條件相近。因此本工作以文獻［8］報道的反應速率常數(shù)為基準，對PET終縮聚反應器的出口物料進行靈敏度分析。反應速率常數(shù)的合理范圍參見文獻［9］。

3.1 反應速率常數(shù)對對苯二甲酸質(zhì)量流量的影響

反應速率常數(shù)對對苯二甲酸質(zhì)量流量的影響見圖1。

圖1 反應速率常數(shù)對對苯二甲酸質(zhì)量流量的影響Fig.1 The effects of the reaction rate constants on themass flow of TPA.

從圖1可知，酯化反應速率常數(shù)(k2和k3)對對苯二甲酸的質(zhì)量流量影響最大，對苯二甲酸的質(zhì)量流量隨k2和k3的增大而迅速降低。與k2相比，k3對對苯二甲酸質(zhì)量流量的影響更大，這是因為在終縮聚階段乙二醇端基含量比乙二醇含量高，故酯化反應(3)(見表2)對對苯二甲酸質(zhì)量流量的影響更顯著;隨熱裂解速率常數(shù)(k7)的增大，對苯二甲酸的質(zhì)量流量增大，這是因為熱裂解使對苯二甲酸端基含量增加;其他反應對體系中對苯二甲酸的質(zhì)量流量和轉化率的影響較小。

3.2 反應速率常數(shù)對二甘醇質(zhì)量流量的影響

終縮聚階段的反應(5)和反應(6)(見表2)生成二甘醇端基和二甘醇重復鏈段，聚合物鏈的形成減少了小分子二甘醇的生成。反應速率常數(shù)對二甘醇質(zhì)量流量的影響見圖2。從圖2可見，k7的增大會使乙二醇端基的含量減小，從而間接使二甘醇生成量減小，故二甘醇的質(zhì)量流量隨k1和k7的增大而降低。其他反應對體系中二甘醇的質(zhì)量流量和轉化率的影響較小。

圖2 反應速率常數(shù)對二甘醇質(zhì)量流量的影響Fig.2 The effects of the reaction rate constants on DEG mass flow.

3.3 反應速率常數(shù)對PET質(zhì)量流量的影響

反應速率常數(shù)對PET質(zhì)量流量的影響見圖3。從圖3可見，隨醇解反應速率常數(shù)(k'1)和 k7的增大，低聚度PET質(zhì)量流量減小，這是因為醇解反應和乙醛生成的反應不利于PET聚合物鏈的形成。其他反應對PET的質(zhì)量流量影響較小。

3.4 反應速率常數(shù)對乙醛質(zhì)量流量的影響

反應速率常數(shù)對乙醛質(zhì)量流量的影響見圖4。從圖4可見，對乙醛質(zhì)量流量影響最大的是k7，其次是k'1。乙醛由熱降解反應生成，參與熱降解反應(8)(見表2)的反應物為乙二醇端基、對苯二甲酸重復單元，其中對苯二甲酸重復單元的濃度變化很小，乙二醇端基的濃度主要由k'1決定，故k'1對乙醛質(zhì)量流量影響較大。參與熱降解反應(9)(見表2)的反應物為乙烯端基、乙二醇端基，乙烯端基來源于熱裂解反應(7)(見表2)，且熱降解反應速率常數(shù)k9較熱降解反應速率常數(shù)k8大4個數(shù)量級，故k7對乙醛質(zhì)量流量的影響最大。

圖3 反應速率常數(shù)對PET質(zhì)量流量的影響Fig.3 The effects of the reaction rate constants on PET mass flow.

圖4 反應速率常數(shù)對乙醛質(zhì)量流量的影響Fig.4 The effect of the reaction rate constants on acetaldehyde(AA)mass flow.

3.5 反應速率常數(shù)對PET的M n的影響

反應速率常數(shù)對PET的Mn的影響見圖5。從圖5可見，Mn/Mn0(Mn0：反應速率常數(shù)為k0時的PET的Mn)隨k1的增大而增大，即PET的Mn隨k1的增大而增大;Mn/Mn0隨k'1的增大而減小，即PET的Mn隨k'1的增大而減小。這是因為縮聚反應有利于聚合物鏈的形成，而醇解反應不利于聚合物鏈的形成;PET的Mn隨k7的增大先增大后減小;其他反應速率常數(shù)對PET的Mn影響較小。

圖5 反應速率常數(shù)對PET的M n/M n0的影響Fig.5 The effect of the reaction rate constants on M n/M n0 of PET.

3.6 傳質(zhì)因子對PET的M n的影響

傳質(zhì)因子對PET的Mn的影響見圖6。從圖6可見，隨傳質(zhì)因子的增大，產(chǎn)物PET的Mn也逐漸增大。因此，PET終縮聚的脫揮是保證產(chǎn)品質(zhì)量的關鍵因素，這也是反應器設計重點考慮的因素。

圖6 傳質(zhì)因子對PET的M n的影響Fig.6 The effect of themass transfer factor(A)on M n of PET.

靈敏度分析結果表明，為準確模擬出口物料質(zhì)量指標如二甘醇、PET、乙醛的質(zhì)量流量和PET的Mn，需對反應速率常數(shù)k1，k'1，k7進行修正;為準確模擬對苯二甲酸的質(zhì)量流量，需修正反應速率常數(shù)k2，k3，k7。

4 PET終縮聚的過程模擬

對酯化過程［9］建立了基于反應機理的模型，并經(jīng)過了3組工業(yè)數(shù)據(jù)的驗證，修正結果表明，所建立的模型具有代表性。過程建模分為3個部分：物性模型、動力學模型和反應器模型(傳質(zhì)模型)。

用于含聚合物體系的物性模型通常有兩種，即活度系數(shù)模型和狀態(tài)方程模型，其中活度系數(shù)模型適用于極性體系，狀態(tài)方程模型較適合非極性體系?；疃认禂?shù)模型比較適合PET生成的多相反應體系，因此采用UNIFAC模型［10］進行相平衡和物性計算。

終縮聚反應器由旋轉的多圓盤片串聯(lián)組成，反應物料流接近平推流。反應器內(nèi)物料的黏度較高，小分子脫除受傳質(zhì)速率的限制。因此，終縮聚釜作為平推流反應器考慮，并耦合傳質(zhì)模型。

根據(jù)靈敏度分析的結果修正式(3)，得到反應速率常數(shù)(k)的指數(shù)前因子(k0)和傳質(zhì)因子見表3。

利用Mark－Houw ink方程式(見式(4))，可根據(jù)工廠實際測得的 PET的 ［η］模擬得到 PET的Mn：

式中，常數(shù)b=7.55×10－4，c=0.685，適用于Mn在3 000～30 000之間、［η］≥0.15 dL/g的PET。

生產(chǎn)負荷為100 t/d的PET終縮聚的模擬結果與工廠的數(shù)據(jù)比較見表4。由表4可見，出口物料中重要的質(zhì)量指標如羧基轉化率、PET的Mn模擬數(shù)據(jù)均與工業(yè)生產(chǎn)的實際數(shù)據(jù)相符，羧基轉化率、PET的Mn、羧基濃度、二甘醇的質(zhì)量分數(shù)的相對偏差分別為0.01%，1.46%，0.37%，0.23%。

表3 PET終縮聚反應k0和傳質(zhì)因子的確定Table The determ ination of the pre-exponential factors(k0)and A in themodel

表4 PET終縮聚釜模擬結果與工廠數(shù)據(jù)比較(100 t/d)Table 4 Comparison between the plant data and themodel results in polycondensation reactor for PET(100 t/d)

5 結論

(1)對PET終縮聚的主副反應機理和動力學參數(shù)進行了系統(tǒng)的分析。

(2)反應速率常數(shù)k1，k'1，k7對出口物料質(zhì)量指標如二甘醇、乙醛和PET的質(zhì)量流量及PET的Mn影響較大;反應速率常數(shù)k2，k3，k7對對苯二甲酸的質(zhì)量流量影響較大;傳質(zhì)因子A對PET的Mn影響較大。

(3)根據(jù)工廠的實際數(shù)據(jù)分別修正了PET終縮聚過程中的各反應速率常數(shù)，建立了規(guī)?；疨ET裝置終縮聚的機理模型。羧基轉化率、PET的Mn、端羧基值等模型數(shù)據(jù)與實測數(shù)據(jù)相對偏差在2%以內(nèi)。

符 號 說 明

［1］ Ravindranath K，Mashelkar R A.Modeling of Poly(Ethylene Terephthalate)Reactors：6.A Continuous Process for Final Stages of Polycondensation［J］.Polym Eng Sci，1982，22(10)：628－636.

［2］ Ravindranath K，Mashelkar R A.Modeling of Polyethylene Terephthalate Reactors：7.MWD Considerations［J］.Polym Eng Sci，1984，24(1)：30－41.

［3］ Ravindranath K，Mashelkar R A.Finishing Stages of PET Synthesis：A Comprehensive Model［J］.AIChE J，1984，30 (3)：415－422.

［4］ Laubriet C，Corre B，Choi K Y.Two-Phase Model for Continuous Final Stage Melt Polycondensation of Poly(Ethylene Terephthalate)：1.Steady-State Analysis［J］.Ind Eng Chem Res，1991，30(1)：2－12.

［5］ Yoon K H，Park O O.Analysis of a Reactor w ith Surface Renewal for Poly(Ethylene Terephthalate)Synthesis［J］.Polym Eng Sci，1994，34(3)：190－200.

［6］ 趙玲，朱中南，戴干策.PET縮聚過程反應與傳質(zhì)：Ⅰ.反應動力學研究［J］.化學反應工程與工藝，2000，16(2)：159－163.

［7］ 趙玲，朱中南，戴干策.PET縮聚過程反應與傳質(zhì)：Ⅱ.傳質(zhì)規(guī)律研究［J］.化學反應工程與工藝，2000，16(2)：164－168.

［8］ Yamada T，Imamura Y.Simulation of Continuous Direct Esterification Process Between Terephthalic Acid and Ethylene Glycol［J］.Polym Plast Technol Eng，1989，28(7/8)：811－876.

［9］ 李文艷，顧雪萍，曾思，等.規(guī)?；疨ET裝置酯化過程模型化［J］.化學工程，2010，38(10)：158－162.

［10］ 曾思，顧雪萍，王嘉駿，等.聚對苯二甲酸乙二醇酯酯化過程組分物性及相平衡計算［J］.石油化工，2009，38(6)：645－650.

Modeling of an Industrial Polycondensation Process for Polyethylene Terephthalate

LiWenyan，Gu Xueping，Wang Jiajun，F(xiàn)eng Lianfang

(State Key Laboratory of Chem ical Engineering，Department of Chem ical and Biochem ical Engineering，Zhejiang University，Hangzhou Zhejiang 310027，China)

The reported kinetics mechanism，mass-transfer mechanism and their parameters for polyethylene terephthalate(PET)polycondensation were summarized，and the corresponding model was established.The effects of the polycondensation rate constant，the esterification rate constant，the hydrolysis rate constant，the thermal degradation rate constant，diethylene glycol formation，acetaldehyde formation and mass-transfer factor on PET polycondensation were studied.The model was revised by using plant data.The results showed that the rate constants for the polycondensation，the hydrolysis and the thermal degradation could evidently affect the mass flows of diglycol，acetaldehyde and PET，and PET relativemolecularmass.The rate constants of the thermal degradation and the esterification had significant effects on the terephthalic acid mass flow.The effect of masstransfer factor on the number average relativemolecularmass(Mn)of PET was obvious.Based on the revised model of PET polycondensation，all the errors between plant data and simulation data for carboxyl conversion，Mnof PET and massmole concentration of the terminal carboxyl group were less than 2%.

polyethylene terephthalate;kinetics;mass-transfer;modeling;polycondensation

1000－8144(2011)10－1094－06

TQ 015.9

A

2011－04－28;［修改稿日期］2011－07－22。

李文艷(1985—)，女，遼寧省盤錦市人，碩士生，電郵liwy_85@yahoo.com.cn。聯(lián)系人：顧雪萍，電話0571－87951307，電郵guxueping@zju.edu.cn。

國家科技支撐項目(2007BAE28B04)。

(編輯 王小蘭)