基于改变相比法结合HS-GC-MS测定烟丝基质中酯类单体香料的固-气分配系数.pdf

2024 年 3 月第 57 卷第 3 期烟草科技Tobacco Science&TechnologyMar.2024Vol.57 No.3摘要：为了借助固-气分配系数对烟草基质中外加风味成分的释放行为进行深入剖析，基于改变相比法结合 HS-GC-MS 分析方法，建立了一种适用于烟丝基质中单体香料分配系数测定的PRV-HS-GC-MS方法，采用PRV-HS-GC-MS法测定了6种酯类香料在烟丝基质（固相）和气相中的分配系数（Kcut/air）。结果表明：在90110 条件下6种酯类香料Kcut/air的相对标准偏差（RSD）为0.92%9.14%，所建立的方法精密度高，可重复性好，便捷准确。烟丝中6种酯类单体香料的Kcut/air值与温度呈负相关关系；以低温区域内Kcut/air值所建立的预测方程可对难以测定的高温区域下的Kcut/air值进行预测，90110 时Kcut/air实际测定值与预测值的相对误差在1.04%5.32%之间，预测效果较好。Kcut/air值可以作为烟丝基质中单体香料释放性能参数，量化表征单体香料释放行为的差异，为烟用单体香料释放行为提供新的研究途径。关键词：酯类化合物；单体香料；烟草基质；分配系数；改变相比法；静态顶空气相色谱中图分类号：TS411.1文献标志码：A文章编号：1002-0861（2024）03-0001-10收稿日期：2023-10-13录用日期：2023-12-18基金项目：国家烟草专卖局创新平台科研活动稳定支持专项经费资助；山东中烟工业有限责任公司科技项目“烟用香料香气物质释放规律研究”（432023DG0240）。第一作者：赵佳明（1999），男，硕士研究生，从事烟用香精香料研究。E-mail：*通信作者：李军（1986），男，高级工程师，从事卷烟调香技术研究。E-mail：；孙世豪（1980），男，研究员，从事风味化学研究。E-mail：引用本文：赵佳明，付英杰，罗曾源，等.基于改变相比法结合HS-GC-MS测定烟丝基质中酯类单体香料的固-气分配系数 J.烟草科技，2024，57（3）：1-10.（ZHAO Jiaming，FUYingjie，LUO Zengyuan，et al.Determing solid-gas partition coefficientsof monomer ester flavors in cut tobacco matrix by changing comparison method with HS-GC-MS J.Tobacco Science&Technology，2024，57（3）：1-10.DOI：10.16135/j.issn1002-0861.2023.0644）基于改变相比法结合HS-GC-MS测定烟丝基质中酯类单体香料的固-气分配系数赵佳明1，付英杰1，罗曾源1，席辉1，卢斌斌1，张帅2，陈帅伟2，王增瑜2，宗永立1，李军*2，孙世豪*11.中国烟草总公司郑州烟草研究院，郑州高新技术产业开发区枫杨街2号4500012.山东中烟工业有限责任公司技术中心，济南市历下区经十路11888号250014Determing solid-gas partition coefficients of monomer ester flavors in cuttobacco matrix by changing comparison method with HS-GC-MSZHAO Jiaming1,FU Yingjie1,LUO Zengyuan1,XI Hui1,LU Binbin1,ZHANG Shuai2,CHEN Shuaiwei2,WANG Zengyu2,ZONG Yongli1,LI Jun*2,SUN Shihao*11.Zhengzhou Tobacco Research Institute of CNTC,Zhengzhou 450001,China2.China Tobacco Shandong Industrial Co.,Ltd.,Jinan 250014,ChinaAbstract：In order to analyze in depth the release behavior of added flavor components in tobaccomatrixusingthesolid-gaspartitioncoefficient,aPRV-HS-GC-MSmethodwasdevelopedtodetermine the partition coefficients of monomeric flavors in cut tobacco matrix based on the changingcomparison method combined with HS-GC-MS analysis.The partition coefficients（Kcut/air）of six esterflavors in cut tobacco matrix（solid phase）and gas phase were determined by the PRV-HS-GC-MSmethod.The results showed that:1）The relative standard deviation（RSD）of Kcut/airof the six esterflavors at 90-110 ranged from 0.92%to 9.14%.The developed method was convenient andaccurate with good repeatability.2）The Kcut/airvalues of the six monomer ester flavors in cut tobacco2024 年烟草科技were negatively correlated with temperature.The prediction equations established on the basis of Kcut/airvalues in low-temperature regions could predict Kcut/airvalues that were difficult to determine inhigh-temperature regions.The relative errors between the actually determined and predicted Kcut/airvalues at 90-110 were between 1.04%and 5.32%,indicating good prediction performance.3）TheKcut/airvalues could be used as a parameter for the release performance of monomer flavors in cuttobacco matrix to quantitatively characterize the differences in monomer flavor release behavior,andit provided a new strategy to study the release behavior of tobacco monomer flavors.Keywords：Ester compound;Monomer flavor;Tobacco matrix;Partition coefficient;Changingcomparison method;Static headspace gas chromatography消费风味产品时，风味成分从产品（固体、液体或凝胶）中释放出来，为消费者所感知1。消费者感受到风味产品的风味特征并非产品中所含风味成分的简单加和，而是消费过程中风味成分释放特征的综合体现，这不仅受产品中风味成分浓度和理化性质的影响，也取决于构成产品主体的其他化学物质（如水、糖、有机酸和蛋白质等），即基质成分2。产品中风味成分的种类差异及其从基质中释放的行为差异共同塑造了差异化和特征化的产品风味。因此，在特定的基质背景下，阐述风味成分的释放行为特征对于产品风味的调控和风格塑造具有理论支撑意义和实际应用价值。卷烟是一种特殊的嗜好型风味产品，为了使消费者得到更加舒适的感官体验，通常会在卷烟中添加一些风味物质以改善和丰富卷烟香味和吃味特征，掩盖杂气，减少刺激性，或者赋予烟气某些特征香气，达到塑造产品风格特征的目的。因此，阐明风味成分在卷烟释放过程中的行为差异，可对今后更加科学合理地使用单体香料提供技术支撑。为了探明单体香料在卷烟中的释放规律，研究人员对不同种类烟用单体香料在烟草基质（卷烟、加热卷烟）中的转移行为进行了深入的研究3-4，通过考察酯类单体香料在烟支内部（烟丝、滤嘴）的持留率、迁移率5-8以及酯类单体香料在烟气中的转移率9-11尝试揭示其在卷烟中的转移规律。风味成分的迁移率/转移率反映的是烟气中风味成分的最终释放结果，对于卷烟加香有很大的指导意义；但无论是单体香料在烟气中的总转移率还是逐口转移率，均是单体香料由烟草基质中释放到气相并转移进入烟气，再经历冷凝吸附和截留12-13等作用，最终在烟气中呈现，其体现的是风味成分最终释放的累积结果，对风味成分在卷烟中的释放行为差异解释较为有限。固-气分配系数（K）定义为分析物在固相中或固相上的浓度与其在气相中质量浓度的平衡比14。固-气分配系数能很好地反映分析物在固相或气相中的分配释放能力，是许多应用中的重要参数，例如，表征食品中风味成分释放能力大小的参数，用于香气释放研究15-16；评估包装材料的安全性能，以尽量减少挥发性有机化合物（VOCs）的释放17。此前，分配系数的测定大多处于液体体系当中，其他固相体系中所使用的微量点滴法18也不适用于烟草基质中单体香料分配系数的测定。因此，开发准确可靠的方法测定烟草基质中目标化合物的固-气分配系数很重要。基于此，本研究中以6种酯类单体香料为研究对象，建立了适用于烟草基质中单体香料固-气分配系数的检测方法，旨在为烟草基质中单体香料释放行为研究提供技术支持。1材料与方法1.1材料、试剂和仪器云南玉溪产C3F等级烟叶切后烟丝样品。乙酸丁酯、乙酸戊酯（99%，阿拉丁生化科技股份有限公司）；乙酸异戊酯、丁酸乙酯、戊酸乙酯、异戊酸乙酯（99%，梯希爱化成工业发展有限公司）；无水乙醇（色谱纯，德国Merck公司）。其理化性质见表1。O37-医用喉头喷雾器（江苏平安医疗器械有限公司）；CP224S 型电子天平（感量 0.000 1 g，德国Sartious公司）；Lab Dancer涡旋振荡器（德国IKA公司）；8890/5977B 气相色谱-质谱联用仪（美国Agilent 公司）；MPS robotic 多功能进样平台（德国Gerstel公司）。表1香气化合物的理化性质Tab.1Physicochemical properties of the studiedaroma compounds酯类单体香料乙酸丁酯乙酸戊酯乙酸异戊酯丁酸乙酯戊酸乙酯异戊酸乙酯相对分子量116.16130.18130.18116.16130.18130.18沸点/126149142120145133 2第 57 卷第 3 期赵佳明，等：基于改变相比法结合HS-GC-MS测定烟丝基质中酯类单体香料的固-气分配系数1.2方法1.2.1改变相比法测定理论分配系数（K）表示一种挥发性化合物在热力学平衡时在固相（Cs）与气相（Cg）中的质量浓度之比：K=CsCg（1）式中：K为热力学平衡状态下香气化合物的分配系数，Cs为平衡时香气化合物在固相（烟丝）中的质量浓度，g/cm3；Cg为平衡时香气化合物在气相中的质量浓度，g/cm3。两相间的质量浓度也可以描述为Cs=msVs（2）Cg=mgVg（3）式中：ms为平衡时香气化合物在固相中的质量，g；mg为平衡时香气化合物在气相的质量，g。因此，联立公式（1）、（2）和（3）可得公式（4）。K=CsCg=msmgVgVs=msmg（4）式中：为顶空瓶中气相（Vg）与固相（Vs）的体积之比。定义顶空瓶内烟丝中香气化合物的初始质量为ms*，ms*可以表示为平衡时两相间的质量之和，即：ms*=ms+mg（5）ms*Vs=msVs+mgVs（6）ms*Vs=msVs+mgVg（7）Cs*=Cs+Cg（8）式中：Cs*为顶空瓶内烟丝中香气化合物的初始质量浓度，g/cm3。联立公式（3），可以表示为Cs*=K Cg+Cg=Cg（K+）（9）Cg=Cs*K+（10）公式（10）最初由Ettre等19建立，等式两边取倒数，得公式（11）：1Cg=KCs*+1Cs*（11）在顶空分析时，色谱峰面积与顶空气相中的平衡浓度成比例，见公式（12）：A=fi Cg（12）式中：A为顶空瓶内待测化合物在给定温度下达到热力学平衡时色谱峰的面积，fi为检测器的响应系数。联立公式（11）与（12）可得：1A=Kfi Cs*+1fi Cs*（13）1A=a+b（14）式中：a=K/（fi Cs*），b=1/（fi Cs*）。K=ab（15）1.2.2加香烟丝及顶空分析样品制备加香烟丝样品：将未添加任何香料的云南C3F切后烟丝过孔径0.250 mm（60目）筛除杂，过滤掉烟丝废渣、烟末和梗签，放置在（221），相对湿度60%3%的恒温恒湿箱中平衡48 h以上。称取50 g平衡后的烟丝，按其质量的0.01%称取1.1节中所选用的6种酯类单体香料，以8 mL无水乙醇稀释后，用医用喉头喷雾器均匀地喷洒在烟丝表面，加香过程中保持翻动以保证加香均匀，将制备好的加香烟丝样品置密封袋中保存48 h。顶空分析样品：将不同质量的加香烟丝分别放入20 mL顶空瓶中，密封后放入MPS多功能进样平台，由全自动进样器置于孵化器中，孵化平衡条件：60，240 min；65，240 min；70，120 min；75，120 min；80，90 min；90，90 min；100，45 min；110，45 min。将制备好的样品进行HS-GC-MS分析。分析条件：MPS多功能进样平台条件：孵化器分别在60、65、70、75、80、90、100、110 下孵化；自动进样平台气密针注射温度：70120（气密针温度设定高于孵化温度10）；气密针填充速率：250 L/s；气体注射速率：250 L/s；进样量：1 000 L。色谱条件：色谱柱：DB-5MS（60 m0.25 mm0.25 m）；载气：氦气；进样模式：分流进样；分流比：101；进样口温度：250；程序升温：50（7 min）10/min250 （1 min）；传输线温度：280；EI离子源温度：230；电离能量：70 eV；四极杆温度：150；扫描模式：选择离子监测；质量扫描范围：30350 amu；溶剂延迟：5 min。1.2.3加香烟丝堆积密度测定20与计算称取20 mL顶空瓶的质量，记为m0。使用称量纸和镊子将加香后烟丝置于20 mL顶空瓶中，并手持顶空瓶在桌面上振动10次，辅助使用涡旋振荡器振荡顶空瓶，使烟丝自由堆积密实，直至烟丝在顶空瓶瓶口处的高度不再变化。称取20 mL顶空瓶与填充 32024 年烟草科技后烟丝的总质量，记为m1，烟丝的堆积密度（，g/mL）按照公式（16）计算：=m1-m0VL（16）=VL-VsVs（17）式中：VL表示顶空瓶体积，cm-3；VS表示顶空瓶中烟丝基质的体积，cm-3。测定烟丝堆积密度，从而确定顶空瓶中所加入不同质量比烟丝对应的气-固两相体积比。使用同一批次加香后烟丝重复5次实验，测得其平均堆积密度为 0.15 g/cm3，RSD 为 2.61%。在体积为 20mL的5个顶空瓶中分别加入0.15、0.30、0.45、0.60和0.75 g的烟丝（对应的 值分别为19.00、9.00、5.67、4.00和3.00）。1.2.4数据处理使用 Origin 2019 软件对 与色谱峰面积倒数（1/A）进行线性回归分析。2结果与讨论2.1基于改变相比法的分配系数测定方法建立将不同体积的待测烟丝分别放入5个顶空瓶中密封，用HS-GC-MS法分别测定每个样品，因为每个样品中待测成分是相同的，且处于同一温度下，所以其分配系数相同。同时，因为每个相热力学系统是相同的，所以联立几次独立的顶空检测结果就可以计算出Kcut/air值。此外，分配系数的测定受温度、平衡时间和基质体系作用力等多方面的影响；其中，与平衡时间对方法精确度的影响很大，因此需要对 与平衡时间两个实验条件进行优化，提高测定方法的准确度。2.1.1顶空固相体积优化采用改变相比法测定分配系数时，需要待测化合物处于热力学平衡状态下，测定每一个 值下所对应的化合物色谱峰面积。的选择对方法的精确度有很大影响，当体积比过小或过大时，均无法得到准确的分配系数Kcut/air值21-22。如图1所示，在添加1.00 g烟丝时，烟丝占据了顶空瓶较大的体积，气密针在扎取采样时存在与烟丝接触的现象，导致气密针吸取体积小于设定体积，色谱峰面积（对应的值为2.00）小于添加0.80 g（对应的 值为2.75）烟丝时的色谱峰面积。本研究中通过优化实验得到了结论：在测定烟丝基质中酯类单体香料在固-气两相间的分配系数时，应避免出现固相体积过大导致气密针与烟丝基质接触的情况，选择值在3.0019.00范围内较为合适。a.优化前b.优化后图160 条件下 优化前后异戊酸乙酯所对应的色谱峰面积Fig.1Chromatographic peak areas corresponding to ethyl isovalerate before and after optimized values at 60 2.1.2顶空平衡温度选择及平衡孵化时间优化顶空平衡温度可选择范围为30200，但是当平衡温度过高时，会破坏原有基质结构，从而影响加香成分在固-气两相间的分配。低碳酯类单体香料热稳定性相对较低，处于烟支冷凝区较低温度时，就已经开始释放。此外，温度较高时其分子活动剧烈，不同酯类单体香料间释放差异较不明显。因此，选择处于冷凝区温度区间的60110 平衡温度作为实验温度，其中，90、100、110 作为分配系数测定的外推温度。为了保证加香成分在设定温度下达到固-气两相间分配平衡，在60110 内考察平衡孵化时间对烟丝中6种低碳酯类加香成分丰度的影响23，结果如图2所示。由图2a可知，随着平衡孵化时间增加，6种酯类单体香料在60 条件下的色谱峰面积不断增加；当平衡孵化时间大于240 min时，色谱峰面积趋近于稳定。图2b中6种酯类单体香料在70 条件下于120 min时基本达到平衡状态。图2c、2d表 4第 57 卷第 3 期明，在80、100 条件下，6种酯类单体香料分别在90、45 min时达到平衡。因此，选择240 min作为60、65 条件下的平衡孵化时间；选择120 min作为70、75 的平衡孵化时间；选择90 min作为80、90 的平衡孵化时间；选择45 min作为100、110 的平衡孵化时间。2.1.3分配系数Kcut/air测定方法评价采用HS-GC-MS法逐个测定每个待测样品的顶空平衡色谱峰面积，结果见图3。通过绘制色谱峰面积倒数（1/A）与 的曲线得到线性回归方程，结果见图4。以异戊酸乙酯为例，将线性方程中的截距（a）与斜率（b），按照k=a/b计算得到其在烟丝基质a.60 b.70 c.80 d.100 图2不同温度下顶空平衡孵化时间对目标物色谱峰面积的影响Fig.2Effects of headspace equilibrium incubation time at different temperatureson chromatographic peak areas of the target compoundsa.60 孵化条件b.70 孵化条件图3不同 值所对应的异戊酸乙酯色谱峰面积Fig.3Peak areas of ethyl isovalerate corresponding to different values赵佳明，等：基于改变相比法结合HS-GC-MS测定烟丝基质中酯类单体香料的固-气分配系数 52024 年烟草科技中的固-气分配系数Kcut/air。以 为自变量（X），色谱峰面积倒数（1/A）为因变量（Y），进行线性回归分析，得线性回归方程和相关系数（R2），并计算5次重复实验的相对标准偏差（RSD），结果见表2。根据文献 24-25 报道，基于改变相比（Phasea.60 孵化条件b.70 孵化条件图4值与对应异戊酸乙酯色谱峰面积倒数（1/A）之间的关系Fig.4Relationships between values and the reciprocals of the corresponding chromatographic peak areas表26种酯类单体香料在不同温度条件下的线性回归方程、相关系数（R2）Tab.2Linear regression equations and correlation coefficients（R2）of six monomer ester flavors at different temperatures酯类单体香料乙酸丁酯乙酸戊酯乙酸异戊酯丁酸乙酯戊酸乙酯异戊酸乙酯酯类单体香料乙酸丁酯乙酸戊酯乙酸异戊酯丁酸乙酯戊酸乙酯异戊酸乙酯酯类单体香料乙酸丁酯乙酸戊酯乙酸异戊酯丁酸乙酯戊酸乙酯异戊酸乙酯酯类单体香料乙酸丁酯乙酸戊酯乙酸异戊酯丁酸乙酯戊酸乙酯异戊酸乙酯60 线性回归方程Y=2.8710-7X+3.7710-6Y=1.1910-7X+2.2910-6Y=9.4310-8X+1.2910-6Y=1.9710-7X+1.8110-6Y=7.4510-8X+1.1010-6Y=1.0610-7X+1.2710-670 线性回归方程Y=1.1610-7X+1.2510-6Y=5.5710-8X+8.8510-7Y=5.4810-8X+6.3110-7Y=8.5610-8X+6.7310-7Y=3.9510-8X+4.5110-7Y=6.2210-8X+5.6710-780 线性回归方程Y=1.2810-7X+9.7010-7Y=5.5710-8X+5.7310-7Y=5.5310-8X+4.1510-7Y=9.0510-8X+4.9010-7Y=3.9410-8X+2.9610-7Y=5.9010-8X+3.9810-7100 线性回归方程Y=1.0610-7X+5.3210-7Y=2.6710-8X+1.7310-7Y=4.2810-8X+2.0910-7Y=1.3510-7X+4.5710-7Y=2.8510-8X+1.5410-7Y=5.8610-8X+2.1810-7R20.9910.9980.9990.9950.9950.996R20.9980.9970.9980.9960.9960.998R20.9980.9900.9990.9980.9980.998R20.9950.9950.9950.9980.9990.996RSD/%7.636.039.143.936.563.54RSD/%5.952.137.381.793.507.92RSD/%1.628.578.764.765.331.20RSD/%5.501.683.652.474.033.1265 线性回归方程Y=1.5410-7X+1.8310-6Y=6.7110-8X+1.1210-6Y=6.0810-8X+7.2110-7Y=1.0010-7X+8.0510-7Y=4.1410-8X+5.1810-7Y=5.6710-8X+6.4910-775 线性回归方程Y=1.0210-7X+8.9210-7Y=4.5610-8X+6.0810-7Y=4.7110-8X+4.3710-7Y=7.8410-8X+5.0010-7Y=3.2910-8X+3.2210-7Y=5.5010-8X+4.1910-790 线性回归方程Y=1.7110-7X+1.1010-6Y=3.8510-8X+3.1510-7Y=6.0310-8X+3.7110-7Y=1.8310-7X+7.7010-7Y=3.9010-8X+3.6210-7Y=7.4110-8X+3.8210-7110 线性回归方程Y=2.7010-7X+1.0610-6Y=5.5010-8X+2.9810-7Y=9.5510-8X+3.8810-7Y=3.0510-7X+8.0010-7Y=5.9710-8X+2.6010-7Y=1.2110-7X+4.2610-7R20.9940.9920.9910.9930.9940.996R20.9950.9960.9980.9980.9960.992R20.9960.9940.9970.9960.9970.999R20.9950.9990.9980.9960.9980.999RSD/%7.613.484.262.662.664.20RSD/%2.693.681.953.830.922.37RSD/%4.073.154.093.151.273.14RSD/%3.402.847.273.453.813.18 6第 57 卷第 3 期ratio variation，PRV）方法，利用PRV-HS-GS-MS检测溶液基质中化合物的分配系数时，色谱峰面积倒数与 的线性回归方程的相关系数（R2）均大于等于0.990。由表2可知，在60110 条件下，6种酯类单体香料方程线性关系良好，R20.991，RSD为0.92%9.14%。表明该方法重复性较好，适用于烟丝等固相基质中酯类香料分配系数的测定。2.2烟丝基质中单体香料的分配系数分析2.2.1酯类单体香料的分配系数与温度的关系温度对烟丝中酯类单体香料释放的影响如表3所示。可知，总的趋势是，随着温度的升高，Kcut/air值逐渐减小。这表明随着温度的升高，单体香料与烟丝基质的吸附结合作用变弱，香料的释放行为增强，在气相中的分配增加。进一步分析互为同系物、同分异构体的香料单体的分配系数可知，互为碳链异构的异戊酸乙酯与乙酸异戊酯，乙酸异戊酯分子的对称性更好，色散力较大，具有较高的结构稳定性，其分配系数也更大。此外，分配系数的大小不仅与化合物自身的结构有关，也与基质有关26-27。在 65 时乙酸丁酯、乙酸异戊酯与异戊酸乙酯的分配系数以及100 时丁酸乙酯与异戊酸乙酯的分配系数较为接近，此时，酯类单体香料与烟丝基质之间的作用力可能是主要的决定因素。表36种酯类单体香料的Kcut/air值Tab.3Kcut/airvalues of six monomer ester flavors酯类单体香料乙酸丁酯乙酸戊酯乙酸异戊酯丁酸乙酯戊酸乙酯异戊酸乙酯60 13.1619.2013.699.1714.7711.9065 11.8616.7311.868.0412.5111.4470 10.8015.8911.507.8711.429.1075 8.7113.349.286.379.797.6280 7.5610.297.505.417.516.7590 6.418.166.164.226.705.16100 5.006.474.873.375.393.72110 3.915.414.062.624.363.53绘制ln K和热力学温度倒数（1/Tabs）之间的回归曲线，如图5所示，可以发现ln K和1/Tabs间基本呈线性关系。因此，可以描述为公式（18），在公式（18）中a和b是回归系数，斜率（b）是反映分析物分配系数与温度关系的直接系数。此外，ln K与1/Tabs之间的线性关系也可以检验分配系数数据的准确性28，即ln K-1/Tabs回归曲线的R20.9时，可证明分配系数测定结果的准确性。虽然从乙酸丁酯和异戊酸乙酯的回归曲线上能观察到存在轻微曲率，但在研究范围内的线性关系（如相关系数所示）仍然良好，这也意味着，在本研究中的温度条件下测定的6种酯类单体香料的固-气分配系数，在合理范围内可以外推至 其 他 温 度 范 围，尤 其 是 高 于 顶 空 最 高 温 度（200）的温度范围25。ln K=a+b/Tabs（18）式中：a和b为回归方程回归系数，a为截距，b为斜率。2.2.2酯类单体香料的分配系数与释放预测酯类单体香料分配系数的大小对应了在该温度下酯类单体香料由烟丝基质释放到气相过程中能力的大小，将ln K 和1/Tabs的回归曲线作为预测方程，以90、100和110 （对应图中温度倒数绝对值2.7510-3、2.6810-3、2.6110-3K-1）作为外推温度进行预测，结果如表4、图6所示。从图6可以看出，6种酯类单体香料分配系数的实际测定值与预测值相对误差不大，预测拟合程度较好。6种酯类单体香料在90、100和110 下分配系数的实际值与预测值的平均相对误差（AE）在1.04%5.32%范围内，预测效果较好，可为实验温度区域之外的分配系数测定提供一定的数据支持。此外，有研究表明，疏水性香气化合物在油溶液中的油-气分配系数在设定温度下可以准确地预测香气化合物的释放速率（g/s）29。目前，关于分配系图56种酯类单体香料Kcut/air的对数值与热力学温度倒数的线性回归分析结果Fig.5Results of linear regression analysis of Kcut/airlogarithmic values and thermodynamic temperaturereciprocals for six monomer ester flavors赵佳明，等：基于改变相比法结合HS-GC-MS测定烟丝基质中酯类单体香料的固-气分配系数 72024 年烟草科技数的研究多是基于恒定的低温度场体系下进行的30，而烟草基质中单体香料的释放往往涉及温度的变化，Kcut/air值反映的是对应温度下该单体香料从烟草基质释放到气相中的能力，对于温度变化过程中单体香料的动态释放的预测存在一定局限性。若假设烟支中烟丝升温是与时间（t）相关的函数，如公式（19）所示：T=（t）+f（19）式中：t 为时间，s；（t）表示时间的函数；f 为常数参数。则ln K和1/T线性回归方程的斜率（b）从一定程度上可解释烟丝基质中酯类单体香料在温度变化过程中的释放行为差异，以及预测其在动态释放过程的释放行为差异。3结论（1）建立了烟草基质中单体香料分配系数测定的改变相比-顶空-气质联用法（PRV-HS-GC-MS）法。该方法操作简便、重复性好，可实现对烟草基质中单体香料分配系数（Kcut/air）值的快速和准确检测。（2）在60110 条件下，6种酯类单体香料Kcut/air值的RSD值为0.92%9.14%，Kcut/air值与温度呈负相关关系。（3）建立了6种酯类单体香料的Kcut/air预测方程，在90110 条件下Kcut/air的实际测定值与预测值平均相对误差在1.04%5.32%范围内，预测效果较好，可实现直接测量困难或在难以测量的较高温度下预测分配系数。参考文献1 刘源，陈艳萍.食品风味研究动态 J.食品安全质量检测学报，2019，10（15）：4830-4831.LIU Yuan，CHEN Yanping.Research trend of foodflavor in China J.Journal of Food Safety and Quality，2019，10（15）：4830-4831.2 Liu Y，Xi H，Fu Y J，et al.Effects of organic acids onthe release of fruity esters in water：an insight at themolecular level J.Molecules，2022，27（9）：2942.3 朱梦薇，胡志忠，务文涛，等.杂环类单体香料在电加热卷烟中的转移行为 J.烟草科技，2022，55（12）：59-68.ZHU Mengwei，HU Zhizhong，WU Wentao，et al.Transfer behaviors of single heterocyclic flavors inelectricallyheatedtobaccoproducts J.TobaccoScience&Technology，2022，55（12）：59-68.4 罗海涛，张峻松.5 种烯酮类香料单体在卷烟烟气中的转移行为分析 J.食品工业，2012，33（3）：136-138.LUO Haitao，ZHANG Junsong.The transfer behavior ofketone flavors in“Jinsheng”cigarettes J.Food Industry，2012，33（3）：136-138.5 黄勇兵.代表性香料单体在中、低焦油卷烟中的转移行为研究 D.长沙：湖南师范大学，2014.HUANG Yongbing.Research on the transfer behavior oftypical flavor monomers in moderate/low tar cigarettesD.Changsha：Hunan Normal University，2014.6 徐如彦，朱梦薇，饶先立，等.醇类香料单体在加热卷烟中的转移行为研究 J.化学试剂，2022，44（6）：894-903.XU Ruyan，ZHU Mengwei，RAO Xianli，et al.Studyon transfer behavior of alcohol flavor monomers in图66种酯类单体香料Kcut/air对数值的预测方程Fig.6Prediction equations for Kcut/airlogarithmic values ofsix monomer ester flavors酯类单体香料乙酸丁酯乙酸戊酯乙酸异戊酯丁酸乙酯戊酸乙酯异戊酸乙酯90 Kcut/air实际值6.418.166.164.226.705.16Kcut/air预测值5.998.206.304.106.755.25AE/%4.972.593.732.991.042.06100 Kcut/air实际值5.006.474.873.375.393.72Kcut/air预测值4.806.615.073.285.444.20AE/%5.221.974.333.003.173.42110 Kcut/air实际值3.915.414.062.624.363.53Kcut/air预测值3.835.344.082.664.393.37AE/%3.022.295.322.483.322.53表46种酯类单体香料在不同温度下的Kcut/air实际值和Kcut/air预测值及其平均相对误差（AE）Tab.4Actual and predicted Kcut/airvalues and average relative errors of six monomer ester flavors at different temperatures 8第 57 卷第 3 期heated cigarettes J.Chemical Reagents，2022，44（6）：894-903.7 沈宏林.十六种烟用香料卷烟迁移行为研究 D.昆明：昆明理工大学，2009.SHEN Honglin.Study on the migration behavior ofsixteencigaretteflavors D.Kunming：KunmingUniversity of Science and Technology，2009.8 张杰，宗永立，周会舜，等.一些醛酮类