影响Scathard-Hildebrand统计热力学理论预测PVC中增塑剂迁移平衡分配系数的因素

张　艳,张钦发,肖少军,汪雪雁,蒋利珍

(华南农业大学食品学院,广东广州 510642)



影响Scathard-Hildebrand统计热力学理论预测PVC中增塑剂迁移平衡分配系数的因素

张艳,张钦发,肖少军,汪雪雁,蒋利珍

(华南农业大学食品学院,广东广州 510642)

为探究Scathard-Hildebrand热力学模型对PVC塑料中有害物迁移预测的适用性和准确性。搜集相关文献中有效迁移实验数据,将计算得出的分配系数实验值Fexp与模型预测分配系数模拟值Fsim进行对比,并探讨温度、迁移物分子量、迁移物醇-水分配系数、模拟液特性等因素对分配系数校正值(模拟值与实验值之比的对数值ln(Fsim/Fexp))的影响。结果表明,PVC材料中迁移物分配系数校正值ln(Fsim/Fexp)与以上几种因素均存在着一定的线性关系。影响因素的分析可为预测模型的修正提供有力的依据。

Scathard-Hildebrand热力学理论,PVC,迁移,分配系数

聚氯乙烯(PVC)是一种被广泛使用的食品包装材料,其在加工过程中常加入增塑剂、稳定剂、抗氧化剂等来改变其功能特性[1]。塑料包装与食品接触过程中,这些添加剂会迁移至食品中,影响食品风味,甚至危害消费者健康。单纯的迁移实验检测较复杂,通过迁移实验建立模型来预测和估算有害化合物迁移平衡量,可为安全食品包装材料的选择和加工提供依据[2]。

目前,塑料中有害物质迁移平衡时分配系数的估算方法主要基于Scatchard-Hildebrand统计热力学理论建立的数学模型[3-7]。这些数学模型能够很好的描述塑料包装材料中单体分子的迁移行为,从而可用来估算迁移平衡时的分配系数[8]。由于数学模型都是基于理想状态下建立的,忽略了实际迁移中的很多影响因素,导致分配系数的模拟值与实验值有一定的差距。随着实验条件的改进和先进仪器的发展,实际迁移量的测定也得到了很大的发展[9-13],得到了很多可靠的实际迁移数据。研究各影响因素与分配系数校正值(模拟值与实验值之比的对数值ln(Fsim/Fexp))之间的数量关系,是一种对所建模型进行修正的有效方法。

1　材料与方法

1.1材料与仪器

PVC,2型北京化工二厂;邻苯二甲酸二甲酯(DMP)、邻苯二甲酸二丁酯(DBP)邻苯二甲酸二(2-乙基己基)酯(DEHP)、邻苯二甲酸二异壬酯(DINP)均为分析纯北京市化学试剂公司;甲醇(色谱纯)、乙醇(分析纯)天津市北方天医化学试剂厂。

表2　PVC中增塑剂迁移分配系数实验值Fexp与模拟值lnFsim及其比值对数值ln(Fsim/Fexp)

注:a、b、c表示实验分配系数Fexp的文献来源,其中a(刘芃岩[17]等,2010);b(黄丽[9]等,2001);C(郑文芝[18]等,2006)。LC-6A高效液相色谱仪日本岛津;CS501-SP型超级恒温水浴槽、AS200电子分析天平、固相萃取装置美国安捷伦公司。

1.2实验方法

当包装材料(高聚物)与食品(或食品模拟物)相接触时,聚合物中低分子化合物(迁移物)会向食品(或模拟液)中迁移,根据Scatchard-Hildebrand统计热力学理论,当迁移达到平衡时,迁移物在聚合物和食品(或模拟液)中理想状态下分配系数为Fsim(模拟液中迁移物的量/聚合物中迁移物的量)[14-15]:

(1)

V迁移物为迁移物的摩尔体积,φ聚合物、φ迁移物、φ模拟液分别为塑料、迁移物、模拟液的体积分数,δ聚合物、δ迁移物、δ模拟液分别为塑料、迁移物、模拟液的溶解度参数;由于迁移物的体积相对于聚合物和模拟液都非常小,因此实际计算时,聚合物和模拟液的体积分数φ聚合物、φ模拟液一般都设为1[14],R为气体常数,T为温度。

根据文献资料查询或计算可得PVC、各迁移物和食品模拟物的相关参数如表1所示[16],把各相关参数代入式(1)中,可计算出PVC中各种增塑剂向不同食品模拟液中迁移平衡的模拟分配系数Fsim。同时根据PVC中化合物向食品模拟液迁移的相关文献中的研究结果[9,17-19]可计算出PVC中增塑剂在各种条件下迁移平衡时分配系数实验值Fexp,计算结果列于表2。

表1　物质特征参数

注:a为资料查询值(Hansen CM[16],2000),b位基团贡献法计算值,c,d均为软件程序计算值,其中c为(molinspiration),d为(Molecular Modeling Pro)。

2　结果与分析

Scatchard-Hildebrand统计热力学理论模拟分配系数与实验实际分配系数的比值对数值ln(Fsim/Fexp),即为Scatchard-Hildebrand统计热力学理论分配系数预测模型的修正系数K。表2表明,PVC中不同迁移物向不同模拟物及在不同的温度下迁移物分配系数的修正系数K都不相同。

2.1温度对修正系数K的影响

由图1可以看出,PVC中增塑剂DBP在3%乙酸、10%乙醇和正己烷中,DEHP在正己烷中的迁移修正系数K均随温度升高而降低。这是因为无论实验分配系数还是模拟分配系数均因温度升高而降低,但模拟分配系数并未考虑摩擦阻力、界面张力、溶胀性等因素对迁移造成的影响。随着温度的升高,PVC会产生膨胀现象,且两相之间的界面作用力会随温度升高而降低,从而导致迁移物更易迁出,实际迁移量增大,实际分配系数(Fexp)增大,修正系数K与温度呈负相关。

图1　修正系数K随温度变化情况Fig.1　Effect of temperature on correction coefficients K

2.2增塑剂分子量对修正系数K的影响

由图2可以看出,PVC中增塑剂DMP(M194)、DBP(M267)、DEHP(M397)、DINP(M430)在蒸馏水中迁移平衡时修正系数K随增塑剂分子量增大而增大。这是由于理论模型认为迁移物只是一质点,但实际迁移中增塑剂分子量越大,其分子结构越庞大,在PVC中所占据的空间也就越大,在同一条件下扩散和迁移受到阻力也就越大,增塑剂迁出也就越困难最终导致实际分配系数Fexp减小而K值增大。

图2　修正系数K随增塑剂分子量变化情况Fig.2　Effect of plasticizers molecular weight on correction coefficients K

2.3增塑剂醇-水分配系数修正系数K的影响

图3　修正系数K随增塑剂醇-水分配系数变化情况Fig.3　Effect of plasticizers alcohol-water partition coefficient on correction coefficients K

邻苯二甲酸酯类增塑剂的极性一般较弱,那么衡量物质疏水能力大小的参数醇-水分配系数logP就会对其迁移平衡时的分配系数造成影响。图3为PVC中四种邻苯二甲酸酯类增塑剂DMP(logP1.5)、DBP(logP4.4)、DEHP(logP7.9)、DINP(logP8.3)在蒸馏水中迁移各增塑剂醇-水分配系数与ln(Fsim/Fexp)的关系。

由图3可以看出,随着增塑剂醇-水分配系数的增大,PVC中增塑剂向水中迁移时分配系数校正系数ln(Fsim/Fexp)也不断增大。这是因为增塑剂的醇-水分配系数越大,说明其极性越弱,根据相似相容原理,迁移物在水中的溶解性也越差,导致PVC中增塑剂向水中迁移的也就更少,实际分配系数Fexp就会随之而减小。而理论模型计算时只考虑到迁移物的溶解度参数,并未涉及醇-水分配系数的影响。

2.4模拟液极性对修正系数K的影响

由图4可知,在各温度下,PVC中增塑剂DBP迁移平衡时分配系数校正系数K随着模拟液极性(正己烷Ps0.06,10%乙醇Ps15.28,3%乙酸Ps15.75)的增大而增大。这是因为模拟液极性越大,聚合物与模拟液的极性相差越大,两相之间的界面张力会随之增大而不利于迁移物的迁出,从而使实验分配系数Fexp降低。所建模型并未考虑模拟液与聚合物材料间的相互作用,从而使修正系数K与模拟液极性呈正相关。

图4　修正系数K随模拟液极性变化情况Fig.4　Effect of simulation liquid polarity on correction coefficients K

3　结论

PVC材料中迁移物分配系数由Scathard-Hildebrand统计热力学理论计算的模拟值与实验值有一定的差异,其比值的对数值ln(Fsim/Fexp)随着迁移物分子量、迁移物醇-水分配系数、模拟液极性的增大而增大,即分配系数实验值与模拟值相差越大;温度越高,校正系数ln(Fsim/Fexp)越小,即分配系数实验值与模拟值越接近。

[1]蔡晶,柴丽月,胡秋辉.食品中邻苯二甲酸酯的检测及安全评价[J].食品科学,2005,26(1):242-245.

[2]E. Helmroth,R Rijk,M Dekker,et al. Predictive modeling of migration from packaging materials into food products for regulatory purposes[J]. Trends in Food Science & Technology,2002(13):102-109.

[3]Cheng F L. Research of mathematical model for plastic packaging materials migration[J].Jilin Inst Chem Tech,2011,28(7):77-78.

[4]D Chung,S E Papadakis,K L Yam.Simple Models for Assessing Migration from Food-packaging Films[J].Food Additives & Contaminants,2002,19(6):611-617.

[5]Welle F. Research of mathematical model for plastic packaging materials migration[J].Package Technol Sci,2013,26(4):56-59

[6]刘志刚,王志伟.塑料包装材料化学物向食品迁移的模型研究进展[J].高分子材料科学与工程,2007,23(5):19-23.

[7]Frank Welle. A New Method for the Prediction of Diffusion Coefficients in Poly(ethylene terephthalate)[J].Journal of Applied Polymer Science,2013,129(4):63-70.

[8]Lau Wong. Mathematical model for the migration of plasticisers from food contact materials into solid food[J]. Analytica Chimica Acta,1997,347(3):249-256.

[9]黄丽,姜志国,张金生.PVC薄膜中增塑剂在水环境中迁移规律研究[J].中国塑料,2001,15(4):51-53.

[10]朱勇,王志伟.食品包装用PVC膜增塑剂迁移的研究[J].包装工程,2006,27(l):40-41.

[11]郭春海,薄海波,贾海涛,等.食品接触材料PVC中32种增塑剂在4种食品模拟物中的迁移规律研究[J].包装工程,2011,32(7):9-14.

[12]王君,许超,杨学军,等. 食品包装用PVC中的物质在5种介质中总迁移规律的研究[J].包装工程,2012,33(19):79-84.

[13]李波平,林勤保. PVC包装材料中化学物向食品迁移的研究[J].食品工业科技,2007,28(1):197-200.

[14]肖少军,张钦发,向红,等.热力学数学模型对塑料中化学物向食品模拟液中迁移平衡时分配系数的预测研究[J].食品安全质量检测学报,2014,5(6):1746-1750.

[15]李敏雯,肖少军,张钦发,等. Scatchard-Hildebrand热力学模型在预测PE中有害物质迁移分配系数的应用[J].食品工业科技,2016(5):86-88.

[16]CharlesM.Hansen. Hansen solubility parameters-A User’s Handbook[M]. CRC,1999,12(17):280-287.

[17]刘芃岩,王有旺,李睿,等.不同食品模拟物对塑料中酞酸酯溶出效果研究[J].安徽农业科学,2010,38(27):15078-15080.

[18]郑文芝,周勇强,张霖霖.PVC塑料制品中增塑剂PAEs在水环境中迁移规律的研究[J].广东化工,2006,33(3):29-30.

[19]李喜宏,袁军伟,马骏,等.PVC保鲜膜中增塑剂DOP在水中迁移规律研究[J].食品科技,2007,4(6):219-221.

Influencing factors of predicting migration partion coefficient of plasticizers in PVC by Scatchard-Hildebrand thermodynamic model

ZHANG Yan,ZHANG Qin-fa,XIAO Shao-jun,WANG Xue-yan,JIANG Li-zhen

(College of Food,South China Agricultural University,Guangzhou 510642,China)

In order to explore the applicability and accuracy of the Scatchard-Hildeb-rand thermodynamic model in predicting the partition coefficient of harmful compounds in PVC migrant to food. The relevant number of effective migration experiment data were collected,then the partition coefficients of experiment data were calculated and compared with the modeling values. The relationship of correction partition coefficient(logarithm the ratio of simulation value and experiment value)between temperature,molecular weight of migrants,alcohol-water distribution coefficient of migrants,polarity of simulants was studied. The results showed that the correction value ln(Fsim/Fexp)of migrants in PVC was linear correlated with the described factors.The analysis of influencing factors provided the powerful basis for correcting predicting model.

Scatchard-Hildebrand thermodynamic theory;PVC;migrant;partition coefficient

2016-03-14

张艳(1993-),女,硕士研究生,研究方向:食品包装材料,E-mail:15039077004@163.com。

张钦发(1963-),男,博士,教授,研究方向:包装新工艺和新材料,E-mail:zqfzgn@163.com。

TS201.2

A

1002-0306(2016)17-0138-04

10.13386/j.issn1002-0306.2016.17.018