3 个产地不同等级庐山云雾茶挥发性成分主成分分析

刘 晔，葛丽琴，王远兴*

（南昌大学 食品科学与技术国家重点实验室，江西 南昌 330047）

庐山云雾茶是中国十大名茶之一，叶厚毫多、醇香甘润、富含营养，以“味醇、色秀、香馨、汤清”闻名[1]，原地产、品质等级、生产年份和相似品种是市场中茶叶掺伪常见的几种类型。常规的化学成分鉴定方法也不足以进行正确的鉴别，因为这些相似茶叶品种有着非常接近的品种亲缘性，其主要品质成分含量差异并不大。快速准确鉴别茶叶品质是当前亟待解决的一项课题[2]。

茶叶的等级主要由其内在的化学成分和冲泡过程中表现出的感官性状所决定。茶叶香气是其所含众多芳香类物质综合表现的结果[3]。挥发性物质的测定多采用气相色谱-质谱（gas chromatography-mass spectrometry，GCMS）[4]和GC-MS-嗅闻[5]的方式，富集常采用固相微萃取（solid phase microextraction，SPME）技术[6]。保留指数（retention index，RI）测定常被用来区分不同的同分异构化合物，是一种重现性和准确度都很好的定性分析参数[7]。主成分分析常采用减少数据集的维数，同时保持对数据集方差贡献最大的特征，通过保留低阶主成分，忽略高阶主成分[8]，可以更加形象地表现原样本的信息如等级分类、产地的鉴别等[9]。

本实验采用SPME法对不同等级庐山云雾茶样本中的挥发性物质进行萃取，并采用GC-MS对样本中挥发性物质进行定性分析。利用主成分分析法对不同等级庐山云雾茶进行等级分类，建立模型，为今后庐山云雾茶等级区分提供一种方便快速的方法。

1 材料与方法

1.1 材料与试剂

供试庐山云雾茶来源于通远镇（5 个等级，每种6 个批次） 庐山茶科所茶场；赛阳镇（4 个等级，每种6 个批次） 庐山露语茶业有限公司；星子县（3 个等级，每种6 个批次） 星子县庐山七尖云雾茶有限公司；正构烷烃（C6～C10、C8～C20）的混合标准品（色谱级） 美国AccuStandard公司；正己烷（色谱纯）德国Merck公司。

等级评定由茶叶采摘日期而定，通常越早采摘的茶叶，产量少等级高。供试品均采于2016年5月份，样品放置于4 ℃冰箱中保存。

1.2 仪器与设备

7890-7000 GC-QQQ-MS气相色谱-三重串联四极杆质谱联用仪（配有MassHunter Workstation色谱工作站和NIST MS Search 2.0质谱检索数据库） 美国Agilent公司；57330-U SPME手柄、50/30 μm二乙烯基苯/碳分子筛/聚二甲基硅氧烷（divinylbenzene/carboxen/polydimethylsiloxane，DVB/CAR/PDMS）萃取纤维头美国Supelco公司；AL 104电子天平 瑞士Mettler-Toledo公司；HH-2数显恒温水浴锅 国华电器有限公司。

1.3 方法

1.3.1 前处理条件

清明前后采摘鲜叶，随海拔高度的增加，开采期亦相应延迟到谷雨前后。以一芽一叶初展为标准，长度在3 cm左右。采回后，薄摊于洁净簸箕上，置于阴凉通风处，保持鲜叶纯净。炒制分杀青、抖散、抖捻、炒二青（初干），理条、搓条、拣剔、提毫、除茶末、烘干等工序。结合宋萌萌[10]、陈熠敏[11]等的前处理方法，将通远镇、赛阳镇、星子县3 个地区的庐山云雾茶，粉碎过筛，用电子天平称取2.00 g茶样于20 mL顶空萃取瓶中，立即用聚四氟乙烯隔垫盖子密封；将顶空萃取瓶置于80 ℃水浴锅中水浴平衡30 min，每隔10 min晃动瓶身；再将固相萃取头插入顶空萃取瓶中，萃取吸附30 min。萃取结束后将萃取头放入GC-MS进样口中，250 ℃解吸5 min，使得挥发性成分得到分离。

1.3.2 GC条件

色谱柱：HP-5MS石英毛细柱（30 m×0.25 mm，0.25 µm）；升温程序：起始柱温50 ℃，保持5 min，以5 ℃/min升至80 ℃，然后以2 ℃/min升至140 ℃，最后以5 ℃/min升至210 ℃保持5 min；载气为高纯氦气（纯度≥99.999%）；载气流速1 mL/min；进样口温度250 ℃；分流比10∶1。

1.3.3 MS条件

电子电离源；电子能量70 eV；传输线温度250 ℃；离子源温度230 ℃；四极杆温度150 ℃；溶剂延迟时间1 min；质量扫描方式为全扫描；质量扫描范围m/z 40～400。

1.3.4 庐山云雾茶挥发性成分RI的测定

将正构烷烃标准品以液体进样的形式，按照1.3.2节升温程序进行分离。进样量为1 µL，分流比为20∶1。记录每个正构烷烃标准品出峰的保留时间，各组分的RI按下式计算，采用Kovats保留指数计算[12]。

RI＝100n＋（tR（x）－tR（n））/（tR（n＋1）－tR（n））×100

式中：tR（x）为待测组分的保留时间/min；n和n＋1分别表示正构烷烃的碳原子数；其中tR（n）＜tR（x）＜tR（n＋1）。

1.4 数据处理

根据NIST MS Search 2.0质谱图库检索结果，并结合RI文献值对庐山云雾茶挥发性成分进行了初步筛查鉴定，并运用峰面积归一化法测得各挥发性成分的相对含量。以样本挥发性物质的相对峰面积为变量，将数据导入SIMCA-P 13.0版软件，将数据标准化后进行无监督的主成分分析，并得到特征值，贡献率及预测能力相关数据。

2 结果与分析

2.1 庐山云雾茶挥发性成分的定性分析

图1 3 个产地庐山云雾茶挥发性成分纵向叠加总离子流色谱图Fig. 1 Overlapped total ion current chromatograms of volatile components in Lu Mountain Clouds-Mist Tea

表1 庐山云雾茶挥发性成分定性分析结果Table 1 Qualitative analysis of volatile components in Lu Mountain Clouds-Mist Tea samples

续表1

如图1、表1所示，样品经过GC-MS检测后，得到较为典型的庐山云雾茶总离子流色谱图。可以看出本实验所采用的GC-MS条件可以较好分离这3 个产地庐山云雾茶的挥发性成分。

2.2 庐山云雾茶挥发性物质的主成分分析

通常来说累计方差贡献率大于80%，认为数据是可信的。如表2、3所示，通远镇庐山云雾茶前4 个成分累计方差贡献率为80.4%并拥有0.673的累计预测能力。赛阳镇庐山云雾茶前3 个成分的累计方差贡献率已经达到87.9%，累计预测能力达到0.807，可见前3 个主成分足以说明该数据的变化趋势。星子县庐山云雾茶前3 个成分的累计方差贡献率已经达到89.8%，累计预测能力达到0.704。累计方差贡献率代表着模型拟合后Y变量（因变量）的变化，累计预测能力代表着通过模型可以被预测的Y变量的变化。累计方差贡献率和累计预测能力值接近1.0说明建立了一个很好的模型，以上可知3 个模型均建立良好。

续表2

表3 3 个产地庐山云雾茶主成分的特征值、贡献率及预测能力Table 3 Eigenvalues, cumulative contribution and predication performance of principal components in Lu Mountain Clouds-Mist Tea from three regions

2.2.1 通远镇庐山云雾茶挥发性物质的主成分分析

由图2可知，通远镇5 个等级庐山云雾茶样品总离子流色谱图相似度很高，除了三等约在保留时间17.5、27.5 min时出现硅峰，疑似少许由萃取头带来的硅氧烷类杂峰[44]，进行样品峰面积计算时应扣除杂峰。结合图2和表2可知，通远镇庐山云雾茶通过富集的挥发性成分主要有：C20苄醇、C30苯乙醇、C71雪松醇、C27十一烷、C44凤梨醛、C42十二烷。挥发性成分包括醇类（14 种）、醛类（8 种）、酮类（7 种）、酯类（5 种）。醇类：三等（47.87%）＞二等（47.74%）＞一等（43.86%）＞普通（40.68%）＞特等（37.99%）；醛类：普通（13.05%）＞特等（12.49%）＞一等（11.57%）＞三等（10.76%）＞二等（8.67%）；酮类：普通（5.41%）＞三等（4.82%）＞二等（4.41%）＞一等（4.37%）＞特等（3.84%）；酯类：特等（1.47%）＝三等（1.47%）＞二等（1.44%）＞普通（1.24%）＞一等（1.17%）。

图2 通远镇5 个等级庐山云雾茶挥发性成分纵向叠加总离子流色谱图Fig. 2 Overlapped total ion current chromatograms of volatile components in five grades of Lu Mountain Clouds-Mist Tea samples from Tongyuan

图3 通远镇5 个等级庐山云雾茶样本的主成分得分图（a）和载荷图（b）Fig. 3 PCA score plots (a) and loading plots (b) of five grades of Lu Mountain Clouds-Mist Tea samples from Tongyuan

由于主成分分析是在没有进行预先分类的情况下，通过每个样本点在得分图上的位置来显示样本的分类信息，因此分类结果更加客观[45]。以各个等级通远镇庐山云雾茶在第1、2主成分上的得分作图，获得所有茶样的二维散点分布图，如图3所示。5 个不同等级的庐山云雾茶在主成分空间上的分布较为分散，各等级的样本均处于相对独立的空间。特等位于图3a右侧，一等位于图3a的中部，二等位于图3a上部，三等位于图3a左侧，普通位于图3a左下侧。由此可见，主成分分析方法对通远镇不同等级云雾茶有较好的区分效果，说明不同等级茶样挥发性成分的相对含量有着一定程度的区别。

表4 通远镇庐山云雾茶挥发性成分中潜在标志物（VIP＞1）Table 4 Potential markers of volatile components in Lu Mountain Clouds-Mist Tea from Tongyuan (VIP ＞ 1)

载荷图表明挥发性成分的分布情况，其分布情况与得分图中样本点的分布和位置对应。重要性因子（variable importance in the projection，VIP）值可以量化每个变量对分类的贡献，VIP值越大，变量在庐山云雾茶不同级别间的差异越显著[46]。通过分析计算出的VIP值，发现有37 个挥发性成分VIP值大于1，由表4所示。可以为通远镇云雾茶的级别分类提供一定的理论支持。并在图3b中，将表4中VIP值大于1标注出来。图3b中右边三角形标记的物质，与图3a中特等茶叶样本位置对应，C3己醛（1.67%）、C4顺-3-己烯-1-醇（2.24%）、C5正己醇（0.77%）、C7庚醛（3.81%）、C8己酸甲酯（0.16%）、C9 α-蒎烯（0.36%）、C15正癸烷（1.09%）、C16辛醛（0.84%）、C28芳樟醇（3.99%），说明这些物质是通远镇特等茶叶的特征物质，并且相对含量较其他等级高。具有香气风味的物质[5,47]如：己醛呈生油脂和青草气及苹果味，浓度低时有水果样的特有香气；顺-3-己烯-1-醇有青叶香气；正己醇有淡青的嫩枝叶香气，微带酒香、果香和脂肪香气；庚醛具有强烈和不愉快的粗糙刺鼻的油脂气味、稀释后具有类似甜杏、坚果香气；α-蒎烯有松节油香气[41]；辛醛有醛香、蜡香、青果皮香及明显的脂肪和水果气味；芳樟醇有柑橘、萜香、铃兰香、玉兰花香、玫瑰花香、甜的青香及木香，香气柔和、持久。说明这些香气成分是通远镇特等庐山云雾茶的重要气味。图3a中二等庐山云雾茶样本和图3b中C21苯乙醛（0.46%）、C40水杨酸甲酯（0.86%）、C47橙花醇甲醚（4.08%）、C69 δ-毕澄茄烯（1.29%）位置类似，相对含量相比于其他等级较高。其中有香气成分如苯乙醛呈强烈风信子香气、玫瑰、巧克力香气，低浓度时有杏仁、樱桃香味；水杨酸甲酯呈冬青油草药香气、薄荷味。说明这些香气成分是体现通远镇二等庐山云雾茶特征香气。图3a中普通等级庐山云雾茶样本和图3b中C13 1-辛烯-3-醇（0.58%）、C24异佛尔酮（0.87%）、C37萘（0.57%）、C38 3-甲基十一烷（0.52%）、C41马鞭草烯酮（0.54%）、C44凤梨醛（5.09%）、C48 2-甲基十二烷（0.09%）、C61表圆线藻烯（0.21%）、C62长叶烯（0.20%）、C63姜黄烯（0.97%）等位置相似，说明这些挥发性成分相对含量较高。其中有香气风味物质：1-辛烯-3-醇有泥土的清香、蔬菜、真菌香气；异佛尔酮有凉香、木香、清甜香气、果香、樟脑及霉香。这些香气是通远镇普通等级庐山云雾茶的特征气味。综合可以看出，通远镇特等庐山云雾茶宜人香气成分种类很多且相对含量很高，随着等级的下降，提供香气的挥发性物质如己醛、顺-3-己烯-1-醇、正己醇、庚醛、α-蒎烯、辛醛、芳樟醇、苯乙醛、水杨酸甲酯等相对含量有所降低。

2.2.2 赛阳镇庐山云雾茶挥发性物质的主成分分析

图4 赛阳镇4 个等级庐山云雾茶挥发性成分纵向叠加总离子流色谱图Fig. 4 Overlapped total ion current chromatograms of volatile components in four grades of Lu Mountain Clouds-Mist Tea samples from Saiyang

由图4可知，赛阳镇4 个等级庐山云雾茶样品总离子流色谱图极为相似。结合图4和表2可知，赛阳镇庐山云雾茶通过本方式富集的主要挥发性成分相对含量由高到低为：C27十一烷、C30苯乙醇、C20苄醇、C48 2-甲基十二烷、C51 4,6-二甲基十二烷、C28芳樟醇等。挥发性成分包括醇类（14 种）、醛类（8 种）、酮类（7 种）、酯类（5 种）。醇类：二等（31.97%）＞三等（30.70%）＞特等（29.33%）＞一等（28.04%）；醛类：一等（12.01%）＞二等（11.89%）＞三等（10.50%）＞特等（10.06%）；酮类：三等（6.29%）＞一等（5.97%）＞特等（4.96%）＞二等（4.92%）；酯类：二等（1.45%）＞三等（1.14%）＝一等（1.14%）＞特等（1.05%）。

以各个茶样在第1、2主成分上的得分作图，获得所有茶样的二维散点分布图，如图5所示。由图5a可以看出，特等位于得分图左侧，距离其他等级样本较远，一等位于图5a右下侧，距离其他等级样本较远。二等和三等茶叶样本距离很近，二等在右侧，三等在左侧，由于等级是按采摘时间分类的，二等和三等采摘时间均在4月下旬，时间较近，以至二等和三等的庐山云雾茶在主成分空间上的分布较为靠拢，而与特等和一等区分明显。可见，主成分分析方法对赛阳镇庐山云雾茶不同等级有较好的区分效果，说明不同等级茶样挥发性成分的相对含量有着一定程度的区别。

图5 赛阳镇4 个等级庐山云雾茶样本的主成分得分图（a）和载荷图（b）Fig. 5 PCA score plots (a) and loading plots (b) of four grades of Lu Mountain Clouds-Mist Tea samples from Saiyang

表5 赛阳镇庐山云雾茶挥发性成分中潜在标志物（VIP＞1）Table 5 Potential markers of volatile components in Lu Mountain Clouds-Mist Tea from Saiyang (VIP ＞ 1)

在图5b中，已将表5中VIP值大于1标注出来。图5a中赛阳镇特等庐山云雾茶与图5b中C21苯乙醛（0.32%）、C48 2-甲基十二烷（6.03%）、C49 1-甲基萘（0.30%）、C51 4,6-二甲基十二烷（4.49%）、C54茉莉酮（1.27%）、C60异丁子香烯（1.75%）、C63姜黄烯（0.06%）、C65 α-广藿香烯（0.49%）、C66 α-衣兰油烯（3.05%）等挥发性物质位置相似，说明其相对于其他等级庐山云雾茶相对含量较高。其中有香气成分如苯乙醛呈强烈风信子香气、玫瑰、巧克力香气，低浓度时有杏仁、樱桃香味；茉莉酮具有花香；异丁子香烯呈紫丁香香气，这些宜人的香气是赛阳镇特等庐山云雾茶代表性的气味。图5a中一等庐山云雾茶与图5b中C2甲苯（1.60%）、C3己醛（3.13%）、C5正己醇（1.57%）、C11苯甲醛（1.87%）、C14 2-正戊基呋喃（3.72%）、C18邻异丙基甲苯（0.79%）、C19柠檬烯（1.19%）、C24异佛尔酮（4.69%）等物质位置相似，相对含量较高。其中有香气风味的物质如己醛呈生油脂和青草气及苹果味，浓度低时有水果样的特有香气；正己醇有淡青的嫩枝叶香气，微带酒香、果香和脂肪香气；苯甲醛具有特殊的苦杏仁气味；2-正戊基呋喃有果香、青香、泥土味和霉臭味；柠檬烯有柠檬、柑橘的果皮香气；异佛尔酮有凉香、木香、清甜香气、果香、樟脑及霉香。图5a中二等云雾茶与图5b中C34橙花醚（0.23%）、C36 3-癸炔-1-醇（2.63%）、C47橙花醇甲醚（4.30%）、C53可巴烯（0.76%）、C69 δ-毕澄茄烯（1.83%）等挥发性物质位置类似，说明这些物质相对含量较高。三等庐山云雾茶在图5a的中间位置，而在图5b中类似位置VIP值大于1的物质很少，说明三等庐山云雾茶特征物质较少。综上可以看出，赛阳镇特等和一等庐山云雾茶具有香气的挥发性物质相对含量较高，随着等级的下降，己醛、正己醇、苯甲醛、2-正戊基呋喃、柠檬烯、异佛尔酮等挥发性物质相对含量降低。

2.2.3 星子县庐山云雾茶挥发性物质的主成分分析

图6 星子县3 个等级庐山云雾茶挥发性成分纵向叠加总离子流色谱图Fig. 6 Overlapped total ion current chromatograms of volatile components in three grades of Lu Mountain Clouds-Mist Tea samples from Xingzi

由图6可知，星子县3 个等级庐山云雾茶样品总离子流色谱图极为相似。结合图6和表2可知，星子县庐山云雾茶通过本方法富集的主要挥发性成分相对含量由高到低有：C20苄醇、C30苯乙醇、C27十一烷、C31 3-亚甲基-1,1-二甲基-2-乙烯基环己烷、C42十二烷、C7庚醛等。挥发性成分包括醇类（13 种）、醛类（8 种）、酮类（7 种）、酯类（5 种）。醇类：一等（48.26%）＞特等（43.72%）＞二等（42.00%）；醛类：一等（9.70%）＞特等（9.30%）＞二等（8.21%）；酮类：一等（2.7%）＞特等（2.67%）＞二等（2.3%）；酯类：二等（2.65%）＞特等（2.59%）＞一等（2.11%）。

图7 星子县3 个等级庐山云雾茶样本的主成分得分图（a）和载荷图（b）Fig. 7 PCA score plots (a) and loading plots (b) of three grades of Lu Mountain Clouds-Mist Tea samples from Xingzi

表6 星子县庐山云雾茶挥发性成分中潜在标志物（VIP＞1）Table 6 Potential markers of volatile components in Lu Mountain Clouds-Mist Tea from Xingzi (VIP ＞ 1)

以各个茶样在第1、2主成分上的得分作图，获得所有茶样的二维散点分布图，如图7所示。由图7a可以看出，特等位于得分图右侧，一等位于得分图左侧，二等位于得分图下部，距离另外2 个等级样本较远。特等和一等庐山云雾茶样本在主成分空间上的分布较为靠拢，而与二等区分明显。可见，主成分分析方法对星子县庐山云雾茶不同等级有较好的区分效果，说明星子县不同等级茶样挥发性成分的相对含量有着一定程度的差别。

通过分析计算出VIP值，发现有29 个挥发性成分VIP值大于1，由表6所示。结合图7a中星子县特等庐山云雾茶所处位置和图7b中C47橙花醇甲醚（3.11%）、C52 α-荜澄茄油烯（0.53%）、C53 可巴烯（0.22%）、C56 α-柏木烯（0.46%）、C61表圆线藻烯（0.17%）、C66 α-衣兰油烯（0.30%）、C67 γ-雪松烯（0.15%）、C69 δ-毕澄茄烯（1.30%）、C73刺柏脑（0.13%）等挥发性物质位置类似，相对含量较另2 个等级高或持平。其中γ-雪松烯具有木香、柏木、类似檀香木并带花香、香水的韵调。图7a中一等庐山云雾茶和图7b中C10 3-乙烯基-环己酮（0.13%）、C11苯甲醛（1.28%）、C18邻异丙基甲苯（0.45%）、C31 3-亚甲基-1,1-二甲基-2-乙烯基环己烷（4.96%）等挥发性物质位置类似，说明这些物质是一等庐山云雾茶的特征物质，相对含量较高。图7a中二等庐山云雾茶与图7b中C21苯乙醛（0.20%）、C23 2,6-二甲基辛烷（0.35%）、C34橙花醚（0.53%）、C35四氢薰衣草醇（1.21%）、C38 3-甲基十一烷（0.22%）、C41马鞭草烯酮（0.26%）、C42十二烷（4.07%）等挥发性物质位置相似，这些物质相对含量较高。C14 2-正戊基呋喃（1.26%）、C44凤梨醛（0.75%）、C60异丁子香烯（0.06%）与二等茶叶样品位置相反，其中2-正戊基呋喃具有果香、青香、泥土味和霉臭味；凤梨醛具有水果香气；异丁子香烯具有紫丁香香气。说明较另外2 个等级茶叶而言，二等庐山云雾茶这些香气成分含量较低。综合可以看出，星子县庐山云雾茶宜人香气成分相对含量随着等级的下降，提供香气的挥发性物质如γ-雪松烯、2-正戊基呋喃、凤梨醛、异丁子香烯等相对含量降低。

3 结 论

本实验通过主成分分析结合GC-MS方法对3 个产地不同等级庐山云雾茶挥发性物质进行了研究。结果显示，在庐山云雾茶样本中共鉴定出76 种挥发性成分，3 个产地不同等级茶叶样本进行主成分分析。其中通远镇庐山云雾茶提取了4 个主成分累计方差贡献率达到80.4%；赛阳镇庐山云雾茶提取了3 个主成分累计方差贡献率达到87.9%；星子县庐山云雾茶提取了3 个主成分累计方差贡献率达到89.8%。从主成分得分图中可直观看出不同等级庐山云雾茶间的差异。该方法能够实现不同产地庐山云雾茶的区分，所建立起来的模型和庐山云雾茶的挥发性物质相对含量信息较吻合。主成分分析不仅能够减少分析的变量，提取有效信息，还能利用建立的模型来区分不同等级的庐山云雾茶。分析结果表明，3 个产地不同等级庐山云雾茶都呈现这样的规律：随着等级的降低，具有香气特征的挥发性成分如己醛、顺-3-己烯-1-醇、正己醇、庚醛、苯甲醛、2-正戊基呋喃、柠檬烯、异佛尔酮、芳樟醇、茉莉酮、异丁子香烯等，相对含量也随之减少。本实验的数据及分析可为今后庐山云雾茶的等级区分提供参考。

参考文献：

[1]陈剑, 刘频健, 吴晓言. 庐山云雾茶多糖提取工艺研究[J].怀化学院学报, 2014, 33(11)： 10-12. DOI：10.3969/j.issn.1671-9743.2014.11.003

[2]林杰. 茶叶香气的图谱分析及在茶叶品质真实性鉴定中的应用[D].杭州： 浙江大学, 2013.

[3]程权, 杨方, 李捷, 等. 全二维气相色谱-飞行时间质谱结合聚类分析与Fisher判别分析对铁观音品质等级的评价研究[J]. 分析测试学报,2015, 34(5)： 525-531. DOI：10.3969/j.issn.1004-4957.2015.05.005.

[4]FENG Y, CAI Y, SUN-WATERHOUSE D, et al. Approaches of aroma extraction dilution analysis (AEDA) for headspace solid phase microextraction and gas chromatography-olfactometry (HS-SPMEGC-O)： altering sample amount, diluting the sample or adjusting split ratio?[J]. Food Chemistry, 2015, 187： 44-52. DOI：10.1016/j.foodchem.2015.03.138.

[5]吕世懂, 孟庆雄, 徐咏全, 等. 普洱茶香气分析方法及香气活性物质研究进展[J]. 食品科学, 2014, 35(11)： 292-298. DOI：10.7506/spkx1002-6630-201411058.

[6]SANG Z Y, WANG Y T, TSOI Y K, et al. CODEX-compliant eleven organophosphorus pesticides screening in multiple commodities using headspace-solid phase microextraction-gas chromatographymass spectrometry[J]. Food Chemistry, 2013, 136(2)： 710-717.DOI：10.1016/j.foodchem.2012.08.060.

[7]尹团章, 熊国玺, 王娜, 等. GC-MS结合保留指数分析烘焙咖啡挥发性成分[J]. 香料香精化妆品, 2016, 6(2)： 9-12. DOI：10.3969/j.issn.1000-4475.2016.02.003.

[8]LILAND K H. Multivariate methods in metabolomics-from preprocessing to dimension reduction and statistical analysis[J]. Trends in Analytical Chemistry, 2011, 30(6)： 827-841.

[9]康受姈. GC-MS结合化学计量学对茶叶种类及品质判断研究[D].杭州： 浙江大学, 2012.

[10]宋萌萌, 王远兴. SPME-GC-QQQ-MS测定两种婺源名优绿茶香气成分[J]. 南昌大学学报(理科版), 2016, 40(4)： 334-340. DOI：10.3969/j.issn.1006-0464.2016.04.007.

[11]陈熠敏, 何洪, 王远兴. 靖安白茶与安吉白茶特征香气成分的比较[J]. 南昌大学学报(理科版), 2015, 39(6)： 573-578. DOI：10.3969/j.issn.1006-0464.2015.06.012.

[12]WANG Y, YANG C, LI S, et al. Volatile characteristics of 50 peaches and nectarines evaluated by HP-SPME with GC-MS[J]. Food Chemistry,2009, 116(1)： 356-364. DOI：10.1016/j.foodchem.2009.02.004.

[13]BONAITI C, IRLINGER F, SPINNLER H E. An iterative sensory procedure to select odor-active associations in complex consortia of microorganisms: application to the construction of a cheese model[J].Journal of Dairy Science, 2005, 88(5)： 1671-1684. DOI:10.3168/jds.S0022-0302(05)72839-3.

[14]XIE J, SUN B, ZHENG F, et al. Volatile flavor constituents in roasted pork of Mini-pig[J]. Food Chemistry, 2008, 109(3)： 506-514.DOI：10.1016/j.foodchem.2007.12.074.

[15]ROUX M L, CRONJE J C, BURGER B V, et al. Characterization of volatiles and aroma-active compounds in honeybush (Cyclopia subternata) by GC-MS and GC-O analysis[J]. Journal of Agricultural and Food Chemistry, 2012, 60(10)： 2657-2664. DOI：10.1021/jf2048383.

[16]IBRAHIM M A, EGIGU M C, KASURINEN A, et al. Diversity of volatile organic compound emissions from flowering and vegetative branches of Yeheb, Cordeauxia edulis (Caesalpiniaceae), a threatened evergreen desert shrub[J]. Flavour and Fragrance Journal, 2010, 25(2)：83-92. DOI：10.1002/ffj.1971.

[17]VARARU F, MORENO-GARCIA J, ZAMFIR C I, et al. Selection of aroma compounds for the differentiation of wines obtained by fermenting musts with starter cultures of commercial yeast strains[J]. Food chemistry, 2016, 197： 373-381. DOI：10.1016/j.foodchem.2015.10.111.

[18]陈玮玲, 钟培培, 范琳琳, 等. 固相微萃取-气相色谱-质谱分析青钱柳叶挥发性成分[J]. 食品工业科技, 2016, 37(22)： 52-58.DOI：10.13386/j.issn1002-0306.2016.22.002.

[19]PINO J A, QUERIS O. Analysis of volatile compounds of mango wine[J]. Food Chemistry, 2011, 125(4)： 1141-1146. DOI：10.1016/j.foodchem.2010.09.056.

[20]MAIETTI S, ROSSI D, GUERRINI A, et al. A multivariate analysis approach to the study of chemical and functional properties of chemodiverse plant derivatives： lavender essential oils[J]. Flavour and Fragrance Journal, 2013, 28(3)： 144-154.

[21]FORERO M D, QUIJANO C E, PINO J A. Volatile compounds of chile pepper (Capsicum annuum L. var. glabriusculum) at two ripening stages[J]. Flavour and Fragrance Journal, 2009, 24(1)： 25-30.

[22]JAVIDNIA K, MIRI R, JAMALIAN A. Composition of the essential oil of Salvia macrosiphon Boiss. from Iran[J]. Flavour and Fragrance Journal, 2005, 20(5)： 542-543. DOI：10.1002/ffj.1468.

[23]FERNANDO L N, GRUN I U. Headspace-SPME analysis of volatiles of the ridge gourd (Luffa acutangula) and bitter gourd (Momordica charantia) flowers[J]. Flavour and Fragrance Journal, 2001, 16(4)：289-293. DOI：10.1002/ffj.999.

[24]WU Y, LÜ S, LIAN M, et al. Study of characteristic aroma components of baked Wujiatai green tea by HS-SPME/GC-MS combined with principal component analysis[J]. CyTA-Journal of Food, 2016, 14(3)： 423-432. DOI：10.1080/19476337.2015.1123298.

[25]申明月, 刘玲玲, 聂少平, 等. 顶空-气相色谱-四极杆质谱结合保留指数法测定普洱茶香气成分[J]. 食品科学, 2014, 35(6)： 103-106.DOI：10.7506/spkx1002-6630-201406021.

[26]ZENG Y X, ZHAO C X, LIANG Y Z, et al. Comparative analysis of volatile components from Clematis species growing in China[J].Analytica Chimica Acta, 2007, 595(1)： 328-339. DOI：10.1016/j.aca.2006.12.022.

[27]吕世懂, 吴远双, 姜玉芳, 等. 不同产区乌龙茶香气特征及差异分析[J].食品科学, 2014, 35(2)： 146-153. DOI：10.7506/spkx1002-6630-201402027.

[28]CAPRIOLI G, IANNARELLI R, CIANFAGLIONE K, et al. Volatile profile, nutritional value and secretory structures of the berry-like fruits of Hypericum androsaemum L.[J]. Food Research International,2016, 79： 1-10. DOI：10.1016/j.foodres.2015.11.021.

[29]CARDILE V, RUSSO A, FORMISANO C, et al. Essential oils of Salvia bracteata and Salvia rubifolia from Lebanon： chemical composition, antimicrobial activity and inhibitory effect on human melanoma cells[J]. Journal of Ethnopharmacology, 2009, 126(2)： 265-272. DOI：10.1016/j.jep.2009.08.034.

[30]MUHAIDAT R, AL-QUDAH M A, SAMIR O, et al. Phytochemical investigation and in vitro antibacterial activity of essential oils from Cleome droserifolia (Forssk.) Delile and C. trinervia Fresen.(Cleomaceae)[J]. South African Journal of Botany, 2015, 99： 21-28.DOI：10.1016/j.sajb.2015.03.184.

[31]JAVIDNIA K, MIRI R, KAMALINEJAD M, et al. Chemical composition of the essential oils of Anthemis altissima L. grown in Iran[J]. Flavour & Fragrance Journal, 2004, 19(3)： 213-216.DOI：10.1002/ffj.1277.

[32]ALISSANDRAKIS E, TARANTILIS P A, HARIZANIS P C, et al.Comparison of the volatile composition in thyme honeys from several origins in Greece[J]. Journal of Agricultural & Food Chemistry, 2007,55(20)： 8152-8157. DOI：10.1021/jf071442y.

[33]SAROGLOU V, MARIN P D, RANCIC A, et al. Composition and antimicrobial activity of the essential oil of six Hypericum species from Serbia[J]. Biochemical Systematics & Ecology, 2007, 35(3)：146-152. DOI：10.1016/j.bse.2006.09.009.

[34]LIU Z, CHEN Z, HAN F, et al. Microwave-assisted method for simultaneous hydrolysis and extraction in obtaining ellagic acid, gallic acid and essential oil from Eucalyptus globulus leaves using Brönsted acidic ionic liquid [HO3S(CH2)4mim]HSO4[J]. Industrial Crops and Products, 2016, 81： 152-161.

[35]张晓珊, 陈图锋, 陈耀文, 等. 固相微萃取-气质联用分析沙茶酱香气物质的研究[J]. 食品与药品, 2014, 16(6)： 410-413. DOI：10.3969/j.issn.1672-979X.2014.06.010

[36]ROUT P K, NAIK S N, RAO Y R. Composition of the concrete,absolute, headspace and essential oil of the flowers of Michelia champaca Linn[J]. Flavour and Fragrance Journal, 2006, 21(6)： 906-911. DOI：10.1002/ffj.1742.

[37]BENKACI-ALI F, BAALIOUAMER A, MEKLATI B Y, et al.Chemical composition of seed essential oils from Algerian Nigella sativa extracted by microwave and hydrodistillation[J]. Flavour and Fragrance Journal, 2007, 22(2)： 148-153. DOI：10.1002/ffj.1773

[38]CUI L, WANG Z Y, ZHOU X H. Volatile constituents in the roots and rhizomes oils of Valeriana amurensis[J]. Journal of Essential Oil Bearing Plants, 2010, 13(1)： 130-134. DOI：10.1080/0972060X.2010.10643801.

[39]XIE Y, WANG J, YANG F, et al. Comparative analysis of essential oil components of two Cryptomeria species from China[J]. Industrial Crops and Products, 2011, 34(1)： 1226-1230. DOI：10.1016/j.indcrop.2011.04.016.

[40]DU L, LI J, LI W, et al. Characterization of volatile compounds of Puerh tea using solid-phase microextraction and simultaneous distillationextraction coupled with gas chromatography-mass spectrometry[J].Food Research International, 2014, 57： 61-70.

[41]TEIXEIRA B, MARQUES A, RAMOS C, et al. Chemical composition and antibacterial and antioxidant properties of commercial essential oils[J]. Industrial Crops and Products, 2013, 43： 587-595.DOI：10.1016/j.indcrop.2012.07.069.

[42]PAPA F, MAGGI F, CIANFAGLIONE K, et al. Volatile profiles of flavedo, pulp and seeds in Poncirus trifoliata fruits[J]. Journal of the Science of Food and Agriculture, 2014, 94(14)： 2874-2887.

[43]DE SOUZA A, LOPES E M C, DA SILVA M C, et al. Chemical composition and acetylcholinesterase inhibitory activity of essential oils of Myrceugenia myrciolides (Cambess.) O. Berg and Eugenia riedelliana O.Berg, Myrtaceae, Brazs[J]. Pharmacognosy, 2010, 20(2)： 175-179.

[44]申东, 沈强, 潘科, 等. 基于HS-SPME/GC-MS定性半定量分析藤茶特征香气[J]. 贵州农业科学, 2010(11)： 67-70. DOI：10.3969/j.issn.1001-3601.2010.11.020.

[45]刘彬球, 陈孝权, 吴晓刚, 等. PCA和PLS-DA用于晒青毛茶级别分类研究[J]. 茶叶科学, 2015, 35(2)： 179-184. DOI：10.3969/j.issn.1000-369X.2015.02.014.

[46]朱荫, 杨停, 施江, 等. 西湖龙井茶香气成分的全二维气相色谱-飞行时间质谱分析[J]. 中国农业科学, 2015, 48(20)： 4120-4146.DOI：10.3864/j.issn.0578-1752.2015.20.013.

[47]GOLDENBERG L, YANIV Y, DORON-FAIGENBOIM A, et al.Diversity among mandarin varieties and natural sub-groups inaroma volatiles compositions[J]. Journal of the Science of Food and Agriculture, 2015, 96(1)： 57-65. DOI：10.1002/jsfa.7191.