基于分子量分布指标的聚酯生产过程模拟方法

doi:10.11949/0438-1157.20190499

摘要/Abstract

摘要：

提出一种以分子量分布为目标的聚酯工业生产过程流程模拟方法，针对通用流程模拟软件Aspen Plus在逐步聚合（step-growth）反应机理下不能计算分子量分布指数和分子量分布曲线的不足，通过外置接口程序从后台运行的Aspen源文件（*.bkp/*.apw）中读取数据，实现分子量分布曲线的可视化和其他质量指标的计算。以聚对苯二甲酸丁二醇酯（PBT）工业生产中的预缩聚过程为建模实例，在Aspen Polymers Plus软件平台上建立基于严格过程机理的聚酯工业预缩过程框架模型。根据分子量及端羧基浓度随动力学参数变化的灵敏度，利用工业操作参数及分析数据整定动力学参数，进一步通过外置接口程序计算PBT的分子量分布，重均分子量、分子量分布指数误差分别为4.2%、5.6%，为聚酯生产过程的分子量分布可视化提供基础。

关键词: 数学模拟, 分子量分布, 聚酯, 反应动力学, 建模方法

Abstract:

A polyester industrial production process simulation method targeting molecular weight distribution was proposed. The general process simulation software Aspen Plus could not calculate the molecular weight distribution index and molecular weight distribution curve under the step-growth reaction mechanism._{To overcome this disadvantage, an external interface program was developed to calculate MWD and PDI with major data from the *.bkp/*.apw file running in the background. This strategy was validated by the pre-condensation process simulation of PBT industrial plant. The pre-condensation process model framework with rigorous process mechanism including phase equilibrium and reaction kinetics was established. The sensitivity analysis of the kinetic parameters to molecular weight and concentration of end carboxyl group was applied for the kinetic parameters tuning. Furthermore, the model was validated by process data and analysis data. The calculation errors of weight average molecular weight and PDI were 4.2% and 5.6%. The simulation strategy can be applied for molecular weight distribution in polyester industrial process.}

Key words: mathematical modeling, molecular weight distribution, polyester, reaction kinetics, modeling strategy

中图分类号:

TQ 316.3

常诚, 冯连芳, 顾雪萍, 陈曦, 张才亮. 基于分子量分布指标的聚酯生产过程模拟方法[J]. 化工学报, 2020, 71(2): 708-714.

Cheng CHANG, Liangfang FENG, Xueping GU, Xi CHEN, Cailiang ZHANG. Development of simulation strategy based on molecular weight distribution in polyester production process[J]. CIESC Journal, 2020, 71(2): 708-714.

图/表 8

图1 以聚合物分子量分布为目标的聚酯工业流程模拟策略

Fig.1 Simulation strategy for polyester industrial process based on molecular weight distribution

表1 缩聚过程反应机理

Table 1 Kinetic reactions for polymerization process

No.	Reaction equations
main reactions
1	$T P A + B D O ? k 1 k 1 / K 1 t T P A + t B D O + H 2 O$
2	$t T P A + B D O ? k 2 k 2 / K 2 b T P A + t B D O + H 2 O$
3	$T P A + t B D O ? k 3 k 3 / K 3 t T P A + b B D O + H 2 O$
4	$t T P A + t B D O ? k 4 k 4 / K 4 b T P A + b B D O + H 2 O$
5	$t B D O ? + t T P A + t B D O ? k 5 k 5 / K 5 b B D O + t T P A + H 2 O$ $t B D O ? + b T P A + t B D O ? k 5 k 5 / K 5 b B D O + b T P A + H 2 O$
side reactions
6	$B D O → k s 1 T H F + H 2 O$
7	$t T P A + t B D O → k s 2 T H F + T P A$ $b T P A + t B D O → k s 2 T H F + t T P A$
8	$t T P A + b B D O → k s 3 T P A + t B U E N$ $b T P A + b B D O → k s 3 t T P A + t B U E N$
9	$b B D O - 1 + b T P A : t B D O - 2 → k s 4 b B D O - 2 + b T P A : t B D O - 1$

表1 缩聚过程反应机理

Table 1 Kinetic reactions for polymerization process

No.	Reaction equations
main reactions
1	$T P A + B D O ? k 1 k 1 / K 1 t T P A + t B D O + H 2 O$
2	$t T P A + B D O ? k 2 k 2 / K 2 b T P A + t B D O + H 2 O$
3	$T P A + t B D O ? k 3 k 3 / K 3 t T P A + b B D O + H 2 O$
4	$t T P A + t B D O ? k 4 k 4 / K 4 b T P A + b B D O + H 2 O$
5	$t B D O ? + t T P A + t B D O ? k 5 k 5 / K 5 b B D O + t T P A + H 2 O$ $t B D O ? + b T P A + t B D O ? k 5 k 5 / K 5 b B D O + b T P A + H 2 O$
side reactions
6	$B D O → k s 1 T H F + H 2 O$
7	$t T P A + t B D O → k s 2 T H F + T P A$ $b T P A + t B D O → k s 2 T H F + t T P A$
8	$t T P A + b B D O → k s 3 T P A + t B U E N$ $b T P A + b B D O → k s 3 t T P A + t B U E N$
9	$b B D O - 1 + b T P A : t B D O - 2 → k s 4 b B D O - 2 + b T P A : t B D O - 1$

图2 预缩过程Aspen建模示意图

Fig.2 Flowsheet for pre-condensation process

图3 动力学参数对数均分子量的影响

Fig.3 Influence of kinetic parameters on MWN

图4 动力学参数对端羧基浓度影响

Fig.4 Influence of kinetic parameters on tTPA1concentration

表2 Aspen模型计算值与分析值比较

Table 2 Comparison of data between analysis and simulation

Data	Concentration of carboxyl end group/（mol/t）	Number average molecular weight
analysis data	22.002	6292
simulation data	21.973	6302
relative error/%	-0.132	0.165

图5 预缩产品分子量分布曲线比较

Fig.5 Comparison of MWD curve between calculated and analysis curve of pre-condensation product

表3 分子量分布模拟值与离线分析值比较

Table 3 Comparison of MWD between simulation and analysis data

Data	Weight average molecular weight	Polymer dispersity index
analysis data	10346	1.644
simulation data	10778	1.736
relative error /%	4.2	5.6

参考文献 31

1	王梅, 张健. 聚酯切片的流变性与可纺性研究[J]. 合成纤维工业, 1995, 18(3): 1-4.
	Wang M, Zhang J. Study on the rheological behaviour and spinnability of pet chip[J]. China Synthetic Fiber Industry, 1995, 18(3): 1-4.
2	Kilkson H. Effect of reaction path and initial distribution on molecular weight distribution of irreversible condensation polymers[J]. Industrial & Engineering Chemistry Fundamentals, 1964, 3(4): 95-105.
3	张美琳, 姜行知, 魏晓波, 等. PET分子量分布的测定及对其性能影响的研究[J]. 聚酯工业, 1992, 5(3): 7-11.
	Zhang M L, Jiang X Z, Wei X B, et al. Research of determination of molecular weight distribution of PET and the effects to properties[J]. Polymer Industry, 1992, 5(3): 7-11.
4	蔡夫柳, 姜胶东, 张明宏, 等. 分子量分布对聚酯性能的影响[J]. 合成纤维工业, 1991, 14(4): 23-29.
	Cai F L, Jiang J D, Zhang M H, et al. Effect of the molecular weight distribution on pet property[J]. China Synthetic Fiber Industry, 1991, 14(4): 23-29.
5	张满, 郝双红, 卜晓军. 辽化聚酯切片的可纺性研究[J]. 聚酯工业, 1992, 5(2): 5-7.
	Zhang M, Hao S H, Bu X J. Study on spinnability of Liaoning Petrochemical Company[J]. Polymer Industry, 1992, 5(2): 5-7.
6	Laus M, Angeloni A S, Galli G, et al. Molecular weight and molecular weight distribution effects on the liquid-crystalline and biphasic behavior of a thermotropic polyester[J]. Macromolecules, 1992, 25(22): 5901-5906.
7	沈瀛坪. TPA和EG连续生产PET的第一酯化反应器操作模拟[J]. 化学反应工程与工艺, 1993, 9(4): 423-431.
	Shen Y P. A mathematical model for reactor simulation of a continuous esterification process between TPA and EG[J]. Chemical Reaction Engineering and Technology, 1993, 9(4): 423-431.
8	沈瀛坪, 赵玲. 聚酯(PET)反应过程研究(Ⅵ): 预缩聚反应过程的数学模拟[J]. 化学反应工程与工艺, 1997, 13(4): 383-388.
	Shen Y P, Zhao L. Studies on reaction process for poly(ethylene terephthalate) (Ⅵ): The continuous prepolymerization process analysis and the study for reaction kinetics of polymerization[J]. Chemical Reaction Engineering and Technology, 1997, 13(4): 383-388.
9	沈瀛坪, 赵玲. 聚酯(PET)反应过程研究(Ⅶ): 后缩聚反应过程研究和数学模拟[J]. 化学反应工程与工艺, 1997, 13(4): 389-394.
	Shen Y P, Zhao L. Studies on reaction process for poly(ethylene terephthalate)(Ⅶ): An experimental research and mathematical model for computer simulation of a continuous final stage polymerization process[J]. Chemical Reaction Engineering and Technology, 1997, 13(4): 389-394.
10	王金堂, 周兴贵, 赵玲, 等. 五釜工艺聚酯工业装置的稳态模拟[J]. 合成纤维工业, 2012, 35(1): 59-63.
	Wang J T, Zhou X G, Zhao L, et al. A steady model for five-reactor PET industrial process[J]. China Synthetic Fiber Industry, 2012, 35(1): 59-63.
11	胡耀辉. 聚酯直接酯化过程数学模型[D]. 北京: 北京化工大学, 2003.
	Hu Y H. A Mathematical model for direct esterification process of PET[D]. Beijing: Beijing University of Chemical Technology. 2003.
12	曾思. 工业聚酯装置酯化过程的建模与流程模拟[D]. 杭州: 浙江大学, 2010.
	Zeng S. Modeling and simulation of esterification section in industrial poly(ethylene terephthalate) process[D]. Hangzhou: Zhejiang University, 2010.
13	李文艳. 工业聚酯装置生产过程建模与分析[D]. 杭州: 浙江大学, 2011.
	Li W Y. Modeling and analysis of the industrial poly(ethylene terephthalate) plant[D]. Hangzhou: Zhejiang University, 2011.
14	Woo B, Choi K Y, Song K H, et al. Melt polymerization of bisphenol-A and diphenyl carbonate in a semibatch reactor[J]. Journal of Applied Polymer Science, 2001, 80(8): 1253-1266.
15	Ye Y, Machado B, Choi K Y. Modeling of solid-state polymerization of bisphenol A polycarbonate[J]. Industrial & Engineering Chemistry Research, 2005, 44(8): 2494-2505.
16	Bokis C P, Orbey H, Chen C C. Properly model polymer processes[J]. Chemical Engineering Progress, 1999, 95(4): 39-52.
17	Fredenslund A, Jones R L, Prausnitz J M. Group-contribution estimation of activity coefficients in non-ideal liquid mixtures[J]. AIChE Journal, 1975, 21(6): 1086-1099.
18	Soave G. Equilibrium constants from a modified Redlich-Kwong equation of state[J]. Chemical Engineering Science, 1972, 27(6): 1197-1203.
19	Mark J E. Polymer Data Handbook [M]. 2nd ed. ACS Publications, 2009: 349-351.
20	黄志恭. 聚酯缩聚工艺及反应器优化分析[J]. 合成纤维工业, 2003, 26(4): 38-41
	Huang Z G. Optimization of polyester condensation process and reactor design[J]. China Synthetic Fiber Industry, 2003, 26(4): 38-41.
21	黄志恭. 聚酯装置反应器剖析[J]. 聚酯工业, 2003, 16(2): 1-5.
	Huang Z G. Analysis of PET reactors[J]. Polyester Industry, 2003, 26(4): 38-41.
22	Pilati F, Manaresi P, Fortunato B, et al. Formation of poly(butylene terephthalate): growing reactions studied by model molecules[J]. Polymer, 1981, 22(6): 799-803.
23	Pilati F, Manaresi P, Fortunato B, et al. Formation of poly(butylene terephthalate): secondary reactions studied by model molecules[J]. Polymer, 1981, 22(11): 1566-1570.
24	Suganuma H, Iida H, Harada H, et al. Kinetic analysis of the formation of poly(1, 4-butyleneterephthalate) prepolymer[J]. Sen'i Gakkaishi, 1987, 43(4): 186-191.
25	Suganuma H, Ikeuchi H, Yamamoto Y, et al. Solid state polycondensation of poly(1,4-butyleneterephthalate)[J]. Sen'i Gakkaishi, 1987, 43(10): 538-543.
26	Iida H, Ikeuchi H, Nakazawa S, et al. Polycondensation of poly(1,4-butyleneterephthalate) prepolymer[J]. Sen'i Gakkaishi, 1987, 43(5): 251-256.
27	黄海南. 对苯二甲酸双羟丁酯的直接酯化合成工艺及动力学研究[D]. 上海: 华东理工大学, 2009.
	Huang H N. Study on synthesis and kinetics of bishydroxybutyl terephthalates through direct esterification[D]. Shanghai: East China University of Science and Technology, 2009.
28	陈星枚, 刘先密, 刘尔嘉, 等. 直接酯化法合成BHBT过程中副产物四氢呋喃生成动力学研究[J]. 聚酯工业, 1996, (4): 8-12.
	Chen X M, Liu X M, Liu E J, et al. Study on kinetics of side product tetrahydrofuran formation in BHBT synthesis by direct esterification[J]. Polyester Industry, 1996, (4): 8-12.
29	Passalacqua V, Pilati F, Zamboni V, et al. Thermal degradation of poly(butylene terephthalate)[J]. Polymer, 1976, 17(12): 1044-1048.
30	Darda P J, Hogendoorn J A, Versteeg G F, et al. Reaction kinetics of polybutylene terephthalate polycondensation reaction[J]. AIChE Journal, 2005, 51(2): 622-630.
31	Flory P J. Molecular size distribution in linear condensation polymers[J]. Journal of the American Chemical Society, 1936, 58(10): 1877-1885.

[1]	程成, 段钟弟, 孙浩然, 胡海涛, 薛鸿祥. 表面微结构对析晶沉积特性影响的格子Boltzmann模拟[J]. 化工学报, 2023, 74(S1): 74-86.
[2]	汪林正, 陆俞冰, 张睿智, 罗永浩. 基于分子动力学模拟的VOCs热氧化特性分析[J]. 化工学报, 2023, 74(8): 3242-3255.
[3]	张蒙蒙, 颜冬, 沈永峰, 李文翠. 电解液类型对双离子电池阴阳离子储存行为的影响[J]. 化工学报, 2023, 74(7): 3116-3126.
[4]	陈朝光, 贾玉香, 汪锰. 以低浓度废酸驱动中和渗析脱盐的模拟与验证[J]. 化工学报, 2023, 74(6): 2486-2494.
[5]	罗来明, 张劲, 郭志斌, 王海宁, 卢善富, 相艳. 1~5 kW高温聚合物电解质膜燃料电池堆的理论模拟与组装测试[J]. 化工学报, 2023, 74(4): 1724-1734.
[6]	禹进, 余彬彬, 蒋新生. 一种基于虚拟组分的燃烧调控化学作用量化及分析方法研究[J]. 化工学报, 2023, 74(3): 1303-1312.
[7]	赵亚静, 胡激江, 介素云, 李伯耿. HTPB引入方式对不饱和树脂改性效果的影响[J]. 化工学报, 2023, 74(2): 883-892.
[8]	陈晨, 杨倩, 陈云, 张睿, 刘冬. 不同氧浓度下煤挥发分燃烧的化学动力学研究[J]. 化工学报, 2022, 73(9): 4133-4146.
[9]	陈玉弓, 陈昊, 黄耀松. 基于分子反应动力学模拟的六甲基二硅氧烷热解机理研究[J]. 化工学报, 2022, 73(7): 2844-2857.
[10]	何聪, 钟文琪, 周冠文, 陈曦. 高海拔地区水泥生料悬浮炉分解特性研究[J]. 化工学报, 2022, 73(5): 2120-2129.
[11]	王敏, 程金兰, 李鑫, 陆晶晶, 尹崇鑫, 戴红旗. 酸性助水溶剂脱除木质素机理分析[J]. 化工学报, 2022, 73(5): 2206-2221.
[12]	杨霄, 丁锐, 李墨含, 宋正昶. 氧浓度对微通道内甲烷均相/非均相耦合反应特性的影响[J]. 化工学报, 2022, 73(12): 5427-5437.
[13]	迟子怡, 刘成伟, 张欲凌, 李学刚, 肖文德. CO氧化偶联反应器模拟与优化[J]. 化工学报, 2022, 73(11): 4974-4986.
[14]	龚翔, 李林森, 姜召. PdCo/SiO₂双金属催化剂用于杂环储氢载体的高效脱氢[J]. 化工学报, 2022, 73(10): 4448-4460.
[15]	张利合, 张凡, 李昌伦, 许德平, 徐振刚, 王永刚. BGL煤气化动力学模型构建与验证[J]. 化工学报, 2022, 73(10): 4668-4678.