间歇冷却结晶过程模型参数及操作敏感性分析

doi:10.11949/0438-1157.20240627

摘要/Abstract

摘要：

针对结晶模型中动力学模型参数和操作参数的不确定性，以间歇冷却结晶过程为研究对象，考虑了基于粒数衡算模型的三种简化结晶过程，评估了参数对产品晶体尺寸分布（CSD）和细晶形成的敏感性。在确保修正Morris筛选法和PAWN法收敛的情况下，分别定量分析了12个参数对晶体生长（CG）、晶体生长-成核（CGN）和晶体生长-成核-溶解（CGND）三种结晶过程模型有关产品指标的局部与全局敏感性，并进行了参数敏感性排序。结果表明，对于不同复杂度的模型情况，参数的敏感性有所不同。针对产品CSD，参数μ、C₀、b和m₀在CG和CGN中表现出显著敏感性，参数在CGND中不敏感；针对细晶，参数b、g和C₀在CGN中表现出显著敏感性，参数C₀、T₀、g、b和μ在CGND中表现出敏感性。通过敏感性分析法可以对不同结晶过程筛选出最重要的扰动参数，这将有助于指导鲁棒优化设计中选择需要考虑的模型误差和操作条件扰动。

关键词: 敏感性分析, 动力学模型, 粒度分布, 粒数衡算模型, 结晶

Abstract:

The study aimed to analyze the uncertainties associated with kinetic model parameters and operating parameters in the crystallization model. Three simplified crystallization process models, based on the population balance equations, were considered to analyze the batch cooling crystallization process. The sensitivity of the parameters to the product crystal size distribution (CSD) and the formation of fine crystals was evaluated. The local and global sensitivities of 12 parameters to the three crystallization process models of crystal growth (CG), crystal growth-nucleation (CGN) and crystal growth-nucleation-dissolution (CGND) were quantitatively analyzed, while the convergence of the modified Morris screening method and PAWN method was ensured. And the parameter sensitivities were sorted. The results show that the sensitivity of parameters is different for model situations of different complexity. For product CSD, the parameters μ, C₀, b, and m₀ exhibit significant sensitivity in CG and CGN, while all parameters are insensitive in CGND. For fine crystals, the parameters b, g, and C₀ exhibit significant sensitivity in CGN, and parameters C₀, T₀, g, b and μ show sensitivity in CGND. Sensitivity analysis can be used to identify the most crucial disturbance parameters for various crystallization processes. This can guide the selection of the consideration of model errors and operating condition perturbations in robust optimization design.

Key words: sensitivity analysis, kinetic model, size distribution, population balance model, crystallization

中图分类号:

TQ 028.5

连斌, 龙妍, 徐啟蕾, 单宝明, 王学重, 张方坤. 间歇冷却结晶过程模型参数及操作敏感性分析[J]. 化工学报, 2024, 75(12): 4587-4595.

Bin LIAN, Yan LONG, Qilei XU, Baoming SHAN, Xuezhong WANG, Fangkun ZHANG. Sensitivity analysis of model parameters and process operation for batch cooling crystallization process[J]. CIESC Journal, 2024, 75(12): 4587-4595.

图/表 8

参考文献 32

1	Zhang F K, Shan B M, Wang Y L, et al. Progress and opportunities for utilizing seeding techniques in crystallization processes[J]. Organic Process Research & Development, 2021, 25(7): 1496-1511.
2	He Y, Gao Z G, Zhang T, et al. Seeding techniques and optimization of solution crystallization processes[J]. Organic Process Research & Development, 2020, 24(10): 1839-1849.
3	Ferguson S, Morris G, Hao H X, et al. Automated self seeding of batch crystallizations via plug flow seed generation[J]. Chemical Engineering Research and Design, 2014, 92(11): 2534-2541.
4	Kalbasenka A N, Spierings L C P, Huesman A E M, et al. Application of seeding as a process actuator in a model predictive control framework for fed-batch crystallization of ammonium sulphate[J]. Particle & Particle Systems Characterization, 2007, 24(1): 40-48.
5	Unno J, Hirasawa I. Partial seeding policy for controlling size distribution of product crystal by batch cooling crystallization[J]. Journal of Chemical Engineering of Japan, 2019, 52(6): 501-507.
6	Aamir E, Nagy Z K, Rielly C D. Evaluation of the effect of seed preparation method on the product crystal size distribution for batch cooling crystallization processes[J]. Crystal Growth & Design, 2010, 10(11): 4728-4740.
7	龚俊波, 孙杰, 王静康. 面向智能制造的工业结晶研究进展[J]. 化工学报, 2018, 69(11): 4505-4517.
	Gong J B, Sun J, Wang J K. Research progress of industrial crystallization towards intelligent manufacturing[J]. CIESC Journal, 2018, 69(11): 4505-4517.
8	Nagy Z K, Braatz R D. Advances and new directions in crystallization control[J]. Annual Review of Chemical and Biomolecular Engineering, 2012, 3: 55-75.
9	Nyande B W, Thomas K M, Takarianto A A, et al. Control of crystal size distribution in batch protein crystallization by integrating a gapped Kenics static mixer to flexibly produce seed crystals[J]. Chemical Engineering Science, 2022, 263: 118085.
10	Mikami T, Kaizu R K. Influence of feed condition on crystal size distribution of potassium alum obtained by unseeded two-stage semi-batch cooling crystallization[J]. Journal of Chemical Engineering of Japan, 2019, 52(3): 317-324.
11	Trampuž M, Teslić D, Likozar B. Crystallization of fesoterodine fumarate active pharmaceutical ingredient: modelling of thermodynamic equilibrium, nucleation, growth, agglomeration and dissolution kinetics and temperature cycling[J]. Chemical Engineering Science, 2019, 201: 97-111.
12	Jordens J, Canini E, Gielen B, et al. Ultrasound assisted particle size control by continuous seed generation and batch growth[J]. Crystals, 2017, 7(7): 195.
13	Zhang F K, Liu T, Chen W X, et al. Seed recipe design for batch cooling crystallization with application to l-glutamic acid[J]. Industrial & Engineering Chemistry Research, 2019, 58(8): 3175-3187.
14	Zádor J, Zsély I G, Turányi T. Local and global uncertainty analysis of complex chemical kinetic systems[J]. Reliability Engineering & System Safety, 2006, 91(10/11): 1232-1240.
15	Wen T, Wang H, Huang W, et al. Research progress on controlling of crystal size distribution (CSD) in crystallization process[J]. Chemical Industry and Engineering, 2021, 38(4): 44-55.
16	Lenka M, Sarkar D. Improving crystal size distribution by internal seeding combined cooling/antisolvent crystallization with a cooling/heating cycle[J]. Journal of Crystal Growth, 2018, 486: 130-136.
17	Nagy Z K. Model based robust control approach for batch crystallization product design[J]. Computers & Chemical Engineering, 2009, 33(10): 1685-1691.
18	关润铎. 基于粒数衡算的结晶过程粒度分布预测与优化[D]. 大连: 大连理工大学, 2017.
	Guan R D. Prediction and optimization of particle size distribution in crystallization processe based on particle number balance[D]. Dalian: Dalian University of Technology, 2017.
19	赵绍磊, 王耀国, 张腾, 等. 制药结晶中的先进过程控制[J]. 化工学报, 2020, 71(2): 459-474.
	Zhao S L, Wang Y G, Zhang T, et al. Advanced process control of pharmaceutical crystallization[J]. CIESC Journal, 2020, 71(2): 459-474.
20	Seki H, Su Y. Robust optimal temperature swing operations for size control of seeded batch cooling crystallization[J]. Chemical Engineering Science, 2015, 133: 16-23.
21	Nagy Z K, Braatz R D. Open-loop and closed-loop robust optimal control of batch processes using distributional and worst-case analysis[J]. Journal of Process Control, 2004, 14(4): 411-422.
22	Long B W, Yang H T, Ding Y G. Impact of seed loading ratio on the growth kinetics of mono-ammonium phosphate under isothermal batch crystallization[J]. Korean Journal of Chemical Engineering, 2016, 33(2): 623-628.
23	Yang Y Y, Xie C X, Ma C, et al. Continuous crystallization kinetics of cefradine in a mixed suspension mixed product removal system[J]. Organic Process Research & Development, 2024. 28(7): 2534-2541.
24	Majumder A, Nagy Z K. Fines removal in a continuous plug flow crystallizer by optimal spatial temperature profiles with controlled dissolution[J]. AIChE Journal, 2013, 59(12): 4582-4594.
25	Aamir E, Nagy Z K, Rielly C D, et al. Combined quadrature method of moments and method of characteristics approach for efficient solution of population balance models for dynamic modeling and crystal size distribution control of crystallization processes[J]. Industrial & Engineering Chemistry Research, 2009, 48(18): 8575-8584.
26	Paleari L, Movedi E, Zoli M, et al. Sensitivity analysis using Morris: just screening or an effective ranking method?[J]. Ecological Modelling, 2021, 455: 109648.
27	Puy A, Lo Piano S, Saltelli A. A sensitivity analysis of the PAWN sensitivity index[J]. Environmental Modelling & Software, 2020, 127: 104679.
28	Al-Zoubi N, Malamataris S. Effects of initial concentration and seeding procedure on crystallisation of orthorhombic paracetamol from ethanolic solution[J]. International Journal of Pharmaceutics, 2003, 260(1): 123-135.
29	Doki N, Kubota N, Yokota M, et al. Production of sodium chloride crystals of uni-modal size distribution by batch dilution crystallization[J]. Journal of Chemical Engineering of Japan, 2002, 35(11): 1099-1104.
30	Öner M, Stocks S M, Sin G. Comprehensive sensitivity analysis and process risk assessment of large scale pharmaceutical crystallization processes[J]. Computers & Chemical Engineering, 2020, 135: 106746.
31	Mohammad K, Rahim S A A, Abu Bakar M R. Effect of seed loading and temperature of seeding on carbamazepine-saccharin co-crystal[J]. Indian Journal of Science and Technology, 2017, 10(6): 1-5.
32	康得军, 邱福杰, 温儒杰, 等. 基于SWMM的LID参数局部与全局敏感性分析[J]. 中国给水排水, 2023, 39(17): 131-138.
	Kang D J, Qiu F J, Wen R J, et al. Analysis on local and global sensitivity of LID parameters based on SWMM[J]. China Water & Wastewater, 2023, 39(17): 131-138.

编辑推荐 0

Metrics

阅读次数

全文

HTML			PDF

最新录用	在线预览	正式出版	最新录用	在线预览	正式出版
1	0	5	16	0	33

来源	本网站	其他网站

次数	38	17
比例	69%	31%

摘要

138

最新录用	在线预览	正式出版

51	0	87

来源	本网站	其他网站

次数	37	101
比例	27%	73%

类型	参数	初始值
操作参数	晶种质量m₀	1.470 g
	晶种分布均值μ	127.780 μm
	晶种分布标准差σ	18.650
	溶液初始温度T₀	40℃
	溶液初始浓度C₀	0.104 g/g
	批次时间t_f	4600 s
动力学参数	生长模型系数k_g	8.5708 μm/s
	生长模型系数g	1.0000
	成核模型系数k_b	0.0380 μm^-3·s^-1
	成核模型系数b	3.4174
	溶解模型系数k_d	10.7192 μm/s
	溶解模型系数d	0.5122

类型	参数	初始值
操作参数	晶种质量m₀	1.470 g
	晶种分布均值μ	127.780 μm
	晶种分布标准差σ	18.650
	溶液初始温度T₀	40℃
	溶液初始浓度C₀	0.104 g/g
	批次时间t_f	4600 s
动力学参数	生长模型系数k_g	8.5708 μm/s
	生长模型系数g	1.0000
	成核模型系数k_b	0.0380 μm^-3·s^-1
	成核模型系数b	3.4174
	溶解模型系数k_d	10.7192 μm/s
	溶解模型系数d	0.5122

参数	CG	CGN		CGND
参数	F₁	F₁	F₂	F₁	F₂
μ	1.392	1.363	0.013	0.049	0.085
σ	0.108	0.109	0.007	0.025	0.002
m₀	0.241	0.235	0.010	0.023	0.027
C₀	0.912	0.893	0.128	0.014	1.450
T₀	0.002	0.002	0.001	0.028	0.827
t_f	0.012	0.012	0.024	0.015	0.073
k_g	0.013	0.013	0.034	0.033	0.069
g	0.070	0.073	0.192	0.043	0.414
k_b	—	0.001	0.009	0.021	0.010
b	—	0.005	0.217	0.019	0.211
k_d	—	—	—	0.011	0.017
d	—	—	—	0.016	0.054

参数	CG	CGN		CGND
参数	F₁	F₁	F₂	F₁	F₂
μ	1.392	1.363	0.013	0.049	0.085
σ	0.108	0.109	0.007	0.025	0.002
m₀	0.241	0.235	0.010	0.023	0.027
C₀	0.912	0.893	0.128	0.014	1.450
T₀	0.002	0.002	0.001	0.028	0.827
t_f	0.012	0.012	0.024	0.015	0.073
k_g	0.013	0.013	0.034	0.033	0.069
g	0.070	0.073	0.192	0.043	0.414
k_b	—	0.001	0.009	0.021	0.010
b	—	0.005	0.217	0.019	0.211
k_d	—	—	—	0.011	0.017
d	—	—	—	0.016	0.054

[1]	黄俊豪, 庞克亮, 孙方远, 刘福军, 谷致远, 韩龙, 段衍泉, 冯妍卉. 干熄炉料钟结构对焦炭布料粒径均匀度影响的模拟研究[J]. 化工学报, 2024, 75(S1): 158-169.
[2]	齐琪, 郭利平, 石李明, 郑映, 潘鹏举. 山梨醇类成核剂改性聚丙烯及其共聚物的结晶行为与性能[J]. 化工学报, 2024, 75(7): 2688-2699.
[3]	马君霞, 李林涛, 熊伟丽. 基于Tri-training GPR的半监督软测量建模方法[J]. 化工学报, 2024, 75(7): 2613-2623.
[4]	陈彦伶, 袁炳志, 王丽伟, 张宸, 朱涵玉. 非平衡条件下金属氯化物-氨工质对的吸附动力学研究[J]. 化工学报, 2024, 75(6): 2252-2261.
[5]	师毓辉, 邢继远, 姜雪晗, 叶爽, 黄伟光. 基于PBM的离心式叶轮内气泡破碎合并数值模拟[J]. 化工学报, 2024, 75(5): 1816-1829.
[6]	李昂, 赵振宇, 李洪, 高鑫. 微波诱导高分散Pd/FeP催化剂构筑及其电催化性能研究[J]. 化工学报, 2024, 75(4): 1594-1606.
[7]	刘静, 杨文博, 吕英迪, 陶胜洋. 喷雾-反溶剂结晶法制备掺杂铝粉的复合微球[J]. 化工学报, 2024, 75(4): 1724-1734.
[8]	曾玉娇, 肖炘, 杨刚, 张意博, 郑光明, 李防, 汪凤玲. 基于机理与数据混合驱动的湿法磷酸生产过程代理建模与优化[J]. 化工学报, 2024, 75(3): 936-944.
[9]	刘昌会, 肖桐, 刘庆祎, 耿龙, 赵佳腾. 多孔二氧化钛强化的相变材料储热机理研究[J]. 化工学报, 2024, 75(2): 706-714.
[10]	李川, 洪振取, 单宝明, 徐啟蕾, 张方坤. 求解多维粒数衡算方程的高阶紧致差分方法[J]. 化工学报, 2024, 75(12): 4513-4522.
[11]	王茂先, 孙启典, 付哲, 华放, 纪晔, 程易. 分子水平动力学模型和机器学习方法相结合研究废弃塑料热解[J]. 化工学报, 2024, 75(11): 4320-4332.
[12]	余留洋, 刘书博, 贾晟哲, 马航, 万邦隆, 苏琦雯, 王静康, 汤伟伟, 贺豫娟, 龚俊波. 电子级磷酸的纯化精制技术发展现状与研究进展[J]. 化工学报, 2024, 75(1): 1-19.
[13]	麻雪怡, 刘克勤, 胡激江, 姚臻. POE溶液聚合反应器内混合与反应过程的CFD研究[J]. 化工学报, 2024, 75(1): 322-337.
[14]	于宏鑫, 邵双全. 水结晶过程的分子动力学模拟分析[J]. 化工学报, 2023, 74(S1): 250-258.
[15]	金正浩, 封立杰, 李舒宏. 氨水溶液交叉型再吸收式热泵的能量及分析[J]. 化工学报, 2023, 74(S1): 53-63.