1、茶叶化学分类编制说明茶叶化学分类国家标准编制说明二一六年八月三十日一、工作简介3 1.1目的意义3 1.2任务来源3 1.3 前期所做的工作3 1.4制标原则及标准化过程4 二、标准主要内容的确定5 2.1本标准方法的引入5 2.2样品制备5 2.3化学成分检测5 2.4第一次特征性因子及判别方法的选择6 2.4.1特征性因子第一次提取6 2.4.2 Bayes判别分析结果7 2.4.3新特征性因子提取8 2.4.4 再次利用Bayes判别分析9 2.4.5 非典型样品的剔除11 2.4.6判别方法的选则11 2.5逐步分类的基本思路11 2.6判别模型的建立及逐步判别12 2.6.1 完全发
2、酵茶判别12 2.6.2红茶黑茶判别13 2.6.3白茶鉴别14 2.6.4乌龙茶（青茶）判别15 2.6.5绿茶和黄茶判别16 2.7逐步判别结果16三、方法验证18四、意见稿讨论20茶叶化学分类国家标准编制说明一、工作简介1.1目的意义 目前，国际上公认的将茶叶分为红茶、绿茶、乌龙茶（青茶）、白茶、黄茶和黑茶六大茶类的方法，主要是依赖其特征性加工工艺的不同，侧重在工艺特征的表述上，品质上侧重于感官审评，而缺乏量化的分类方法和指标，加上不同茶类之间加工工艺的相互融合和借鉴，单纯依赖加工工艺和感官指标已经很难实现对茶类进行正确的鉴别，这不仅对消费者对茶叶品质的判别带来了困难，而且在国际贸易中，
3、由于茶类不同感官审评评判规则有异，也容易造成对茶样认识的差异而导致贸易纠纷。 随着植物化学分类学和天然植物化学成分分析检测技术、化学计量学的快速发展和应用，为量化的化学分类法在茶叶分类上的应用提供了可能。制定茶叶化学分类方法不仅符合植物化学学科发展趋势的要求，而且更为重要的是作为茶叶行业中茶叶是为基本的分类标准，对于推进茶产业的发展，促进茶叶国际间的公平、公正贸易具有重要作用。基于上述理由，提出茶叶化学分类方法标准化学分类的定义和基本要求。1.2任务来源 2013年7月22日以国标委综合201356号国家标准委关于下达2013年国家第一批制修订计划的通知中，茶叶化学分类（计划编号：201310
4、66-T-442）由安徽农业大学负责起草，全国供销总社杭州茶叶研究院、福建农林大学、中国农科院茶叶研究所、安徽省标准化研究院和华南农业大学共同参与。1.3前期工作 标准制定项目任务下达后，安徽农业大学组织专家成立了标准编制工作组，负责编制起草工作。编制组成立后，我们进行了广泛的调查和研究工作，收集茶样，茶样来源见表1。表1 茶样产地及数量茶类数量国内国外绿茶113 110：浙江20；云南1；四川3；陕西12；山东15；江西5；江苏7；湖南10；湖北5；河南1；贵州2；广西1；广东2；福建2；安徽 24 3：日本、越南、韩国各1个红茶277 136：浙江2；云南32；四川1；山东5；江西11；江
5、苏1；湖北5；河南5；广东3；福建22；安徽49 141：越南3；印尼21；印度18；乌拉维3；乌干达3；坦桑尼亚2；斯里兰卡26；莫桑比克15；缅甸3；肯尼亚40；布隆迪7 白茶87 86：福鼎48；政和38 1：斯里兰卡 1乌龙茶（青茶）94 94：闽北乌龙11；武夷岩茶15；大红袍3；广东乌龙4；闽南色种30；铁观音31 黄茶21 21：湖北4；浙江9；湖南2；安徽6 黑茶58 58：云南25；四川6；湖南12；湖北4；广西11 合计650 505 145 围绕茶叶的特性和标准适用的范围，重点考虑了以下问题：茶样收集的范围；茶样的等级确定。为了统一茶叶中儿茶素、咖啡碱和氨基酸的检测方法，
6、安徽农业大学、全国供销总社杭州茶叶研究院、福建农林大学、安徽省标准化研究院四个单位选派检测人员到全国供销总社杭州茶叶研究院进行统一培训。1.4制标原则及标准化过程 在广泛调查研究和必要的实验验证的工作基础上，根据GB/T 1.1-2009标准化工作导则 第1部分：标准的结构和编写给出的内容和格式编写完成茶叶化学分类。 本标准的制定编写过程，以提高茶类的判别率及减少判别因子为总的指导原则。使其能体现基本使用性能，确保标准既保持技术上先进性，又有经济上合理性。遵循了标准制定过程中的先进性和适用性。 在标准制定过程中严格遵循国家有关方针，严格执行国家标准，与同体系标准及相关的各种基础标准以及配套使用
7、的取样、检测等标准相衔接，遵循了政策和协调统一性原则。 在标准制定过程中力求做到：技术内容的叙述正确无误，文字表达准确、简明、易懂，标准的结构严禁合理，内容编排、层次划分等符合逻辑与规定。二、标准主要内容的确定2.1本标准方法的引入2.1.1 GB/T 8302 茶 取样2.1.2 GB/T 8303 茶 磨碎试样的制备及其干物质含量测定2.1.3 GB/T 8304 茶 水分测定2.1.4 GB/T 8312 茶 咖啡碱测定 2.1.5 GB/T 8313 茶叶中茶多酚和儿茶素类含量的检测方法 2.1.6 2.2样品制备 2.2.1 样品的收集充分考虑年份、产地、等级和季节等因素，所收集的样
8、品均匀分布，且能代表该茶类所覆盖的范围。对于不具备该茶类典型特点的茶样，通过感官审评的办法将其剔除。考虑到生茶是由晒青毛茶压制而成，其品质特点与晒青绿茶接近，在收集黑茶样品的时候，未将生茶统计在内。 2.2.2 样品收集由安徽农业大学、福建农林大学负责，样品由茶叶加工企业免费提供，或从茶叶加工企业采购，样品来源可信。样品收集后，从中称取一定量样品置于已知质量的铝盒中于电热烘箱中烘干，测定样品水分含量。 2.2.3样品制备由安徽农业大学统一负责制备。茶样经复火后用粉碎机粉碎，使用牢固的包装密封袋保存样品，并放在5冰柜中保存。按照GB/T 8302的方法，将样品6等份，一份留存，另外5份分别统一分
9、发到全国供销总社杭州茶叶研究院，中国农业科学院杭州茶叶研究所，福建农林大学，安徽省标准化研究院和安徽农业大学5家单位同时检测。2.3化学成分检测 茶样的化学成分如表2所示。表2 干茶化学成分含量化学成分含量（干物质%）绿茶红茶白茶乌龙茶（青茶）黄茶黑茶GAEGC+C咖啡碱EGCGECECG儿茶素总量简单儿茶素茶氨酸0.070.032.431.150.090.063.260.696.961.690.840.231.990.6712.302.583.361.261.050.440.240.090.400.420.210.283.050.820.350.540.240.270.550.451.751
10、.440.860.720.830.370.220.120.530.390.150.193.950.543.981.900.280.161.460.646.42.830.960.531.200.590.110.082.130.790.070.032.280.574.471.710.600.191.080.368.352.492.800.960.210.160.150.071.360.870.080.043.060.624.541.670.510.201.810.918.33.461.951.071.110.680.320.220.450.400.030.032.700.930.360.510.1
11、70.150.150.181.161.050.640.500.040.05注：表中值均为该类茶平均值2.4第一次特征性因子及判别方法的选择2.4.1特征性因子第一次提取将测定的8个有效化学成分及其构成的组合变量，运用SPSS数据处理软件进行主成分分析，得到的原始数据解释总方差如表3所示。从统计学上来说，累积方差贡献率达85%以上，选取的新变量就能基本反映原始变量的综合信息，即可作为抽提出来的主成分。从表3可以看出，前五个主成分的累积贡献率达85.08%，满足要求，可作为新变量来代替原始数据并作分析判别。表3 原始数据解释总方差表成份初始特征值提取平方和载入合计方差的 %累积 %合计方差的 %累
12、积 %17.44241.34641.3467.44241.34641.34624.37924.32765.6744.37924.32765.67431.5498.60574.2791.5498.60574.27941.0806.00280.2801.0806.00280.28050.8634.79585.0750.8634.79585.07560.7774.31889.39470.6843.79993.19280.4082.26695.45890.3011.67197.129100.1961.09098.220110.1540.85799.077120.0640.35599.432130.04
13、00.22399.655140.0300.16999.824150.0210.11899.942160.0110.058100.0017-9.397E-17-5.221E-16100.000根据上述主成分分析结果，我们对前五个主成分进行因子分析（如表4）。因子分析是主成分分析的推广，目的是选择每个主成分中的载荷贡献率最高的因子来代替该主成分，从而代表原始变量进行判别分析。因子对每个主成分的载荷贡献值越大，才能全面的反映原始信息，并作为主要成分因子来替代主成分。由表4可以得到五个主成分的代表因子，分别为儿茶素总量、茶氨酸*咖啡碱、EGC/儿茶素总量、咖啡碱/儿茶素总量、+C。表 4 主成分因子筛
14、选结果%成分12345GA-0.0710.143-0.2760.1980.111EGC0.758-0.3750.418-0.0240.069+C0.3700.025-0.2870.1560.827咖啡碱-0.0420.768-0.0360.4050.033EGCG0.9500.0450.0050.052-0.094EC0.866-0.1430.174-0.0200.148ECG0.8740.274-0.1350.081-0.019儿茶素总量0.975-0.0240.0860.041-0.019简单儿茶素0.824-0.3310.357-0.0150.138茶氨酸0.4560.8350.1360
15、.097-0.064茶氨酸20.4170.7820.0860.192-0.117EGC/儿茶素总量-0.425-0.4450.5750.2190.008EGCG/儿茶素总量0.8120.043-0.358-0.255-0.164咖啡碱/儿茶素总量-0.5890.0180.1790.529-0.020ECG/儿茶素总量-0.3170.613-0.260-0.3830.085茶氨酸/儿茶素总量-0.4920.5790.453-0.3600.169茶氨酸2/儿茶素总量-0.3610.6520.481-0.3180.159茶氨酸*咖啡碱0.3860.8800.0950.153-0.0682.4.2 B
16、ayes判别分析结果运用上述分析得到的主成分代表因子作为判别因子（儿茶素总量、茶氨酸*咖啡碱、EGC/儿茶素总量、咖啡碱/儿茶素总量、+C）来进行Bayes判别，判别分类结果如表5所示。我们从表中结果可以看出，六大类茶判别的判别率很低，绿茶的判别率仅为66.4%，黄茶的判别率仅为22.2%，不符合统计学判别分离要求；典则判别效果图显示分类效果较差，同类茶分布不集中，不同种茶类相互重叠，所以提取的成分因子无法实现六大茶类的判别。表5 六大茶类Bayes判别分类结果茶类绿茶红茶白茶乌龙茶（青茶）黄茶黑茶判别率66.4%76.3%70.0%94.7%22.2%69.0% 1：绿茶；2：红茶；3：白茶
17、；4：乌龙茶（青茶）；5：黄茶；6：黑茶图1六大茶类Bayes判别分析结果分布图 运用主成分分析和因子分析提取得到的成分因子并不能代替原始数据进行判 别分析，需要重新对测定的化学成分及其组合变量进行因子分析。从表3可见，第一主成分与第二主成分的特征值较大，抓住这个特点，重点分析前两个主成分。表3有前两个主成分的因子载荷，每一列代表一个主成分作为原来变量线性组合的系数。同一变量在这两个主成分上，有不同的载荷（即它同各主成分的相关系数），相关系数越大，主成分对它的代表性越高，也就是说该变量对这个主成分起主要作用。在第一主成分中，EGC、EGCG、EC、ECG、儿茶素总量、简单儿茶素、EGCG/儿茶
18、素总量的载荷较大，均起主要作用；第二主成分中，咖啡碱、茶氨酸、茶氨酸2、茶氨酸*咖啡碱的载荷较大，说明这4个因子起主要作用。仔细观察上述成分因子，不难发现，第一主成分中起主要作用的都是儿茶素类因子，在这里我们把它们命名为第一类因子；而第二主成分中起主要作用的是咖啡碱与茶氨酸相关因子，我们把这类因子命名为第二类因子。这两个主成分相互独立，所以我们接下来可以考虑把所有因子分成这样的两类因子，重新分析。2.4.3新特征性因子提取 根据上述分析，我们得到第一类因子，包括EGC、EGCG、EC、ECG、儿茶素总量、简单儿茶素和EGCG/儿茶素总量共7个因子，因子总数较多，因此需要对第一类因子做主成分分析
19、，从表6可以看出前两个主成分的累计方差贡献率已达到91.78%，说明前两个主成分可以代表第一类因子做判别分析。表6 原始数据解释总方差表成份初始特征值提取平方和载入合计方差的 %累积 %合计方差的 %累积 %15.50178.58378.5835.50178.58378.58320.92413.20091.7830.92413.20091.78330.3274.66596.44840.1642.34998.79750.0821.17699.97360.0020.027100.0071.464E-152.091E-14100.00 由表7可知，第一个主成分中儿茶素总量的因子载荷最大，所以提取儿茶
20、素总量作为第一主成分的代表因子；第二个主成分中EGCG/儿茶素总量的因子载荷最大，所以提取EGCG/儿茶素总量作为第二个主成分的代表因子。表7 主成分因子筛选结果成分12EGC0.842-0.497EGCG0.9490.224EC0.905-0.192ECG0.8480.378简单儿茶素0.895-0.436儿茶素总量0.9900.070EGCG/儿茶素总量0.7560.502 综上所述，我们找到了第一类因子的代表因子即儿茶素总量与EGCG/儿茶素总量。结合第二类因子（咖啡碱、茶氨酸、茶氨酸2、茶氨酸*咖啡碱），共选取6个因子作为原始变量的代表因子。2.4.4 再次利用Bayes判别分析 以儿
21、茶素总量、EGCG/儿茶素总量、咖啡碱、茶氨酸、茶氨酸2和茶氨酸*咖啡碱这6个因子作为新变量对六大茶类做判别分析。如表8可见，六大茶类的判别率达到了84.6%，较之前的判别率已经提高了很多，但仍不符合统计学上达到85%的要求。结合图2，我们发现红茶、黑茶、白茶和乌龙茶（青茶）的判别率达到85%以上，但是黄茶和绿茶的判别率很低，分别为75.0%、57.1%。表8 利用新因子Bayes交叉验证分类结果编号预测组成员合计123456初始计数187021017011622852125973114304150430091009453006120216000705158%175.00.01.78.614.
22、70.0100.022.187.62.11.02.15.2100.032.02.086.00.08.02.0100.043.20.00.096.80.00.0100.0514.30.00.028.657.10.0100.060.00.00.012.10.087.9100.0交叉验证a计数18602101801162184313597311400715043009010945300 6120216000705158%174.10.01.78.615.50.0100.021.086.63.11.03.15.2100.032.02.080.00.014.02.0100.043.20.00.095.7
23、1.10.0100.0514.30.00.028.657.10.0100.060.00.00.012.10.087.9100.0a. 仅对分析中的案例进行交叉验证。 在交叉验证中，每个案例都是按照从该案例以外的所有其他案例派生的函数来分类的。b. 已对初始分组案例中的 84.6% 个进行了正确分类。c. 已对交叉验证分组案例中的 83.3% 个进行了正确分类。d. 1：绿茶；2：红茶；3：白茶；4：乌龙茶（青茶）；5：黄茶；6：黑茶1：绿茶；2：红茶；3：白茶；4：乌龙茶（青茶）；5：黄茶；6：黑茶图2利用新因子Bayes判别六大茶类分布图2.4.5 非典型样品的剔除在数据分析结果中，我们发现
24、绿茶与黄茶的化学成分含量非常接近，建模过程中同样发现绿茶与黄茶样品较难区分。在这样的情况下，我们决定结合感官审评对Bayes判别分析中错判的绿茶与黄茶样品进行筛选。 邀请安徽农业大学的3位茶叶审评专家对误判的样品进行盲审。感官审评的内容是干茶的色泽、冲泡后茶样的香气、滋味、汤色，如果感官审评后茶叶的类别还是被错判了，说明这个茶样不具备它本身该有的品质特点，不具有典型性，更不具有代表性，所以我们予以剔除2.4.6判别方法的选择 从上述判别的结果来看，利用Bayes判别总体判别率不高，于是我们尝试采用Fisher逐步判别的方法。Fisher判别是以各未知茶样距离判别函数组质心处的距离来进行判别分类
25、，该方法在区分少数种类时具有直观、快速简便的优势，且最优为两类之间的判别，经过不断地分析与总结，发现逐步Fisher判别是区分六大茶类的最优思路。2.5逐步分类的基本思路将化学成分结合Fisher模式识别方法，建立判别方程，并确定质心点距离。由于fisher方法每次只能对两类进行判别，所以茶类的鉴别逐步实施，利用判别方程逐步地鉴别，直到鉴别出该茶类别，整体思路如下图3。图3 六类茶逐步判别思路2.6判别模型的建立及逐步判别根据主成分分析结果将咖啡碱含量（X1）、儿茶素总量（EGCG、ECG、 EGC、 EC、+C）（X2）、茶氨酸含量（X3）、EGCG含量/儿茶素总量（X4）、茶氨酸含量茶氨酸含量（X5）、茶氨酸含量咖啡碱含量（X6）这六个因子作为茶叶化学分类的特征性因子，采用Fisher方法，得出判别方程和质心点距离，判别方程及判别步骤如下： 2.6.1 完全发酵茶判别，将红茶和黑茶判为一个类别，绿茶、白茶、黄茶和乌龙茶（青茶）判为另外一个类别分开,判别方程为：Fisher1 = 0.150X1+0.254X2 -1.364X3+6.175X4+1.126X5-0.271
