王勝鵬，龔自明*，鄭鵬程，劉盼盼，高士偉，滕 靖，王雪萍，葉 飛，鄭 琳，桂安輝

（湖北省農(nóng)業(yè)科學(xué)院果樹茶葉研究所，湖北 武漢 430064）

青磚茶是我國獨有的一種黑茶，主產(chǎn)于湖北省赤壁市，距今已有200余年歷史。青磚茶加工時，渥堆是非常關(guān)鍵的一道工序，會發(fā)生一系列以多酚類物質(zhì)為主的生化反應(yīng)，形成青磚茶茶湯特有的湯色明亮、滋味醇和的品質(zhì)特點。由于青磚茶具有分解脂肪[1]、舒暢腸胃[2]、抗氧化[3]和降血糖[4]等功效，已成為少數(shù)民族群眾日常生活中的一種必需品。當(dāng)前，對青磚茶的研究主要集中在加工技術(shù)[5]、香氣成分分析[6]、菌類鑒別[7]和生理保健[8]上，而在青磚茶茶湯品質(zhì)評價方面的研究較少。青磚茶茶湯含有茶多酚、茶褐素、茶紅素和表沒食子兒茶素沒食子酸酯（epigallocatechin gallate，EGCG）等多達(dá)數(shù)十種內(nèi)含成分，這些成分與茶湯滋味品質(zhì)都具有一定的關(guān)聯(lián)關(guān)系，因此，如何對青磚茶茶湯滋味品質(zhì)進行評價非常重要。

通常應(yīng)用感官審評方法[9]對茶湯品質(zhì)進行評價。感官審評是專業(yè)人員運用自身視覺、味覺、嗅覺和觸覺等人體感官器官對評審物品質(zhì)進行綜合評判的一門實用技術(shù)，已實現(xiàn)對小麥[10]、水稻[11]、蜂蜜[12]和中藥[13]等進行感官審評。應(yīng)用感官審評方法已對綠茶[14]、紅茶[15]、烏龍茶[16]和普洱茶[17]等茶類進行品質(zhì)評價。傳統(tǒng)感官審評法雖然經(jīng)典，但專業(yè)性較強，易受審評員嗜好差異、身體狀況和環(huán)境等多種因素影響，結(jié)果主觀性強，重復(fù)性較差，且費時費力，因此，亟待研發(fā)一種科學(xué)、客觀、便捷的品質(zhì)評價方法。

雖然電子舌技術(shù)已用于酒類和茶湯滋味識別[18]，但常用線性方法建立模型，沒有考慮非線性方法，預(yù)測結(jié)果常與實際結(jié)果不符；測色技術(shù)已用于茶葉色澤感官審評方面研究[19]，但茶湯內(nèi)含成分復(fù)雜，應(yīng)用該方法建立的線性模型難以確保預(yù)測結(jié)果的穩(wěn)定性和一致性，還需進行深入研究。

青磚茶茶湯品質(zhì)高低與其內(nèi)含成分含量間具有密切關(guān)系，但茶湯內(nèi)含成分豐富，如何尋找一種或者多種與茶湯滋味品質(zhì)密切相關(guān)的內(nèi)含成分非常關(guān)鍵。因此，本實驗通過應(yīng)用國標(biāo)方法先測定茶湯內(nèi)含成分化學(xué)物質(zhì)含量，再應(yīng)用SPSS 19.0統(tǒng)計軟件對這些內(nèi)含成分進行主成分分析（principal component analysis，PCA）和相關(guān)關(guān)系分析，篩選代表性內(nèi)含成分，得出各成分間的相關(guān)關(guān)系，再應(yīng)用多種化學(xué)計量學(xué)方法嘗試建立茶湯滋味品質(zhì)得分與代表性內(nèi)含成分含量間的擬合方程，用以評價茶湯品質(zhì)，為青磚茶茶湯滋味品質(zhì)評價提供一種準(zhǔn)確、快捷的方法。

1 材料與方法

1.1 材料

青磚茶樣品共158 個，148 個樣品分別來自湖北省赤壁市羊樓洞茶業(yè)股份有限公司、趙李橋茶廠有限責(zé)任公司、洞莊茶業(yè)有限公司和思莊茶業(yè)有限公司，10 個樣品來源未知，用于檢驗回歸方程的實際預(yù)測效果；加工時間2015年5月—2015年9月，將樣品沖泡后進行感官審評，而后立刻測定內(nèi)含成分含量。按照3∶1比例隨機劃分樣品，其中111 個樣品用于建立校正擬合方程，37 個樣品用于驗證擬合方程的預(yù)測效果。

1.2 儀器與設(shè)備

手動固相萃取進樣器、50/30 μm DVB/CAR/PDMS固相微萃取頭 美國Supeclo公司；HHS型恒溫水浴鍋上海博迅實業(yè)有限公司醫(yī)療設(shè)備廠；UV-2550紫外-可見光分光光度計 日本島津公司；Milli-RO PLUS 30純水機 法國Millipore公司；CM-5色差計 柯尼美能達(dá)投資有限公司；2695高效液相色譜（high performance liquid chromatography，HPLC）、2998二極管陣列檢測器、AccQ.Tag C18氨基酸專用分析柱 美國Waters公司；ODS C18反相柱（250 mm×4.6 mm，5 μm）美國Phenomenex公司。

1.3 方法

1.3.1 茶湯感官審評

根據(jù)GB/T 23776—2009《茶葉感官審評方法》，用四分法取得青磚茶樣品5.0 g，置于250 mL評審杯中，注滿沸水，加蓋浸泡5 min，按照沖泡次序依次等速將茶湯瀝入評茶碗中，茶湯以滋味醇厚、回味干爽為最佳。3 名感官審評專家對青磚茶茶湯品質(zhì)進行打分，滿分為100 分，茶湯品質(zhì)越好，分?jǐn)?shù)越高。審評過程中要求評審員身體健康，個人衛(wèi)生條件好，全程禁止吸煙。

1.3.2 內(nèi)含成分檢測

應(yīng)用國家標(biāo)準(zhǔn)方法檢測水浸出物含量[20]、茶多酚含量[21]、游離氨基酸含量[22]；蒽酮-硫酸法測定總糖含量[23]；三氯化鋁法測定黃酮含量[24]；系統(tǒng)分析法測定茶紅素、茶黃素和茶褐素含量[25]；HPLC法檢測沒食子酸、咖啡堿和兒茶素和其單體含量。

在英語的句法分析中，Ratnaparkhi最早引入了最大熵的方法，他利用上下文特征，通過最大熵的方法來預(yù)測下一步要執(zhí)行的操作。其上下文特征主要包括：成分的核心詞，核心詞的組合，非特定組合信息，以及部分已完成的子樹信息。

1.4 數(shù)據(jù)分析

應(yīng)用SPSS 19.0軟件對青磚茶茶湯內(nèi)含成分進行相關(guān)關(guān)系分析[26]；應(yīng)用Matlab 2012a軟件程序包對青磚茶茶湯內(nèi)含成分進行PCA[27]。應(yīng)用的化學(xué)計量學(xué)方法為主成分回歸方法和逐步回歸方法。所得結(jié)果用校正集決定系數(shù)（determination coefficient of cross validation，）、預(yù)測集決定系數(shù)（determination coefficient of prediction，Rp2）、交互驗證均方根誤差（root mean square error of cross validation，RMSECV）、預(yù)測均方根誤差（root mean square error of prediction，RMSEP）[28]和因子數(shù)表示。其中，R2越大、RMSEP越小時，方程預(yù)測效果越好。RMSECV、RMSEP、R2的計算見式（1）～（3）：

式中：n為校正集樣品數(shù)；yi為樣品i實測值；yi’為校正集樣品i預(yù)測值。

式中：n為驗證集樣品數(shù)；yi為樣品i實測值；yi’為驗證集樣品i預(yù)測值。

式中：yi和yi’分別為樣本中第i個樣本實測值和預(yù)測值；y為所有樣本實測值的平均值。

2 結(jié)果與分析

2.1 茶湯滋味評分及內(nèi)含成分測定結(jié)果

根據(jù)感官審評專家對茶湯滋味進行評分，然后應(yīng)用化學(xué)方法對茶湯內(nèi)含成分進行測定，結(jié)果見表1。

表1 茶湯滋味得分及其內(nèi)含成分測定結(jié)果Table 1 Chemical composition and sensory score of Qingzhuan tea infusion

從表1可以看出，茶湯滋味分?jǐn)?shù)最大值為93，最小值為85，平均值為88.37，標(biāo)準(zhǔn)偏差為2.21，樣品茶湯滋味得分間差距較小。19 種內(nèi)含成分中含量較高的分別為CAF、EGCG、兒茶素總量、茶多酚、水浸出物、茶紅素、茶褐素和總糖。含量較低的成分分別為GA、GC、EGC、C、ECG、CG和ECG。物質(zhì)成分的含量高低是相對的，這些成分對茶湯滋味的形成都起到一定作用，但對茶湯滋味貢獻(xiàn)度并不同。為求得對青磚茶茶湯滋味貢獻(xiàn)較大的內(nèi)含成分，舍棄貢獻(xiàn)較小的內(nèi)含成分，因此需要對上述內(nèi)含成分進行PCA和相關(guān)關(guān)系分析。

2.2 PCA和相關(guān)關(guān)系分析

表2 茶湯內(nèi)含成分PCA貢獻(xiàn)率Table 2 Contribution rates of first 6 principal components to total variance

從表2可以看出，PC1～PC6前6 個主成分的累計貢獻(xiàn)率達(dá)到99.98%，根據(jù)PCA原理，前6 個主成分就可以完全解釋全部內(nèi)含物的化學(xué)信息。PC1和PC2分別可以解釋84.65%和13.17%的青磚茶19 種主要內(nèi)含成分的變異，總解釋率為97.82%。從表3可以看出，青磚茶茶湯滋味得分與EGC（0.09%～0.46%）、EGCG（0.91%～2.38%）和ECG（0.1%～0.73%）呈顯著正相關(guān)（P＜0.05）。EGC、EGCG和ECG為兒茶素單體，均具有較強抗氧化活性，能夠清除體內(nèi)的自由基，具有抗衰老的功效。青磚茶原料在渥堆過程中，會發(fā)生多酚類物質(zhì)的氧化反應(yīng)，但由于翻堆不及時，無法跟外界空氣進行充分接觸，一些兒茶素物質(zhì)未被完全氧化而被保留了下來，是茶湯爽口和刺激性的成分，有利于茶湯的品質(zhì)提升。青磚茶茶湯滋味得分與C（1.99%～4.42%）和茶多酚（3.40%～11.11%）含量呈顯著正相關(guān)（P＜0.05）。茶多酚是一類存在于茶樹中的多元酚混合物，其主要成分是以兒茶素為主體的黃烷醇類物質(zhì)，在茶湯中主要呈現(xiàn)澀味，量多時不利于茶湯品質(zhì)的提高。由于茶多酚含有大量—OH化學(xué)鍵，青磚茶原料在渥堆時極易氧化變性成其他類物質(zhì)，導(dǎo)致茶多酚類物質(zhì)含量大量降低，可以起到降低茶湯澀味的功效，同時保留的茶多酚還具有抗氧化、抗病毒及殺菌作用[29]。從表3還可以看出，EGC、EGCG、ECG、兒茶素總量和茶多酚5 種成分間均呈極顯著正相關(guān)（P＜0.01），彼此間相互起正作用，共同影響著青磚茶茶湯品質(zhì)。除上述5 種成分外，茶褐素（3.40%～8.13%）與茶湯品質(zhì)得分的相關(guān)關(guān)系最大，為0.41，但二者間并不呈顯著相關(guān)（P＞0.05），雖然咖啡堿（1.18%～3.37%）和茶紅素（1.17%～4.62%）的貢獻(xiàn)值也較大，但它們與茶湯品質(zhì)得分間并無顯著相關(guān)（P＞0.05）?？梢?，茶湯內(nèi)含成分含量高低與茶湯滋味品質(zhì)得分間并無必然聯(lián)系，內(nèi)含成分含量高，不一定對茶湯品質(zhì)得分貢獻(xiàn)就大，內(nèi)含成分含量低，也不一定對茶湯品質(zhì)得分貢獻(xiàn)就小。從表3還可以看出，青磚茶茶湯滋味得分與EGC、EGCG、ECG、兒茶素總量和茶多酚的相關(guān)系數(shù)分別為0.49、0.57、0.50、0.56和0.51，與之對應(yīng)的P值分別為0.049 8、0.023、0.049 6、0.025 5和0.044，可見，這5 種成分對茶湯品質(zhì)得分的貢獻(xiàn)大小順序為EGCG、兒茶素總量、茶多酚、ECG和EGC。

表3 青磚茶茶湯組成成分間相關(guān)關(guān)系Table 3 Correlation between the components in Qingzhuan tea soup

2.3 茶湯滋味得分回歸方程建立

以EGC、EGCG、ECG、兒茶素總量和茶多酚成分為自變量，以茶湯滋味得分值為因變量，以SPSS 19.0軟件為平臺，分別應(yīng)用主成分回歸方法和逐步回歸方法，構(gòu)建內(nèi)含成分與滋味得分間的回歸方程。

當(dāng)取前4 個主成分時，建立的主成分回歸方程結(jié)果最佳，其方程為：

式中：y為滋味得分；x1為EGCG；x2為C；x3為茶多酚；x4為ECG；x5為EGC。

方程中樣品實測值和預(yù)測值間R2c為0.911 7，RMSECV為0.421 1；然后應(yīng)用主成分回歸方程對37 個樣品茶湯滋味分?jǐn)?shù)進行預(yù)測，得到為0.889 8，RMSEP為0.432 6。

2.3.2 逐步回歸擬合方程建立

在建立逐步回歸方程過程中，根據(jù)自變量間的組合以及對方程的貢獻(xiàn)大小，建立13 個回歸方程。

表4 逐步回歸方程結(jié)果Table 4 Results of step-wise regression equation

從表4可以看出，根據(jù)自變量對茶湯滋味得分的貢獻(xiàn)率而建立的13 個回歸方程中，隨著方程中自變量個數(shù)增加，回歸方程的預(yù)測效果逐漸提高，其中效果最差的回歸方程只應(yīng)用x1單個自變量（R2c=0.494 4，RMSECV=1.632 2），用此回歸方程對驗證樣品茶湯滋味分?jǐn)?shù)進行預(yù)測時，Rp2=0.490 1，RMSEP=1.654 5；當(dāng)應(yīng)用2 個自變量建立回歸方程時，以應(yīng)用x1和x3建立的回歸方程預(yù)測較佳（=0.628 0，RMSECV=1.450 3），對驗證樣品茶湯滋味分?jǐn)?shù)的預(yù)測結(jié)果Rp2=0.600 5，RMSEP=1.472 4，這是由于x3對茶湯滋味分?jǐn)?shù)的貢獻(xiàn)要大于x5，所以回歸方程預(yù)測結(jié)果較佳，但此時建立回歸方程的信息較少，方程預(yù)測結(jié)果欠佳；同理，以應(yīng)用3 個自變量x1、x3和x2建立的回歸方程預(yù)測結(jié)果（R2c=0.838 5，RMSECV=1.324 0；R2p=0.811 7，RMSEP=1.376 2）要優(yōu)于其他3 個自變量建立的回歸方程得到的預(yù)測結(jié)果，也優(yōu)于全部應(yīng)用2 個自變量或者1 個自變量建立回歸方程的預(yù)測結(jié)果。這主要是由于回歸方程的自變量信息不充分造成的。因此，當(dāng)以x1、x2、x3、x4、x5全部5 個自變量建立回歸方程時預(yù)測結(jié)果最佳（Rc2=0.973 5，RMSECV=0.380 7），當(dāng)用此回歸方程對驗證樣品茶湯滋味分?jǐn)?shù)進行預(yù)測時，R2p=0.968 1，RMSEP=0.400 0（圖1）。

建立的最佳逐步回歸方程為：

式中：y為滋味得分；x1為EGCG；x2為C；x3為茶多酚；x4為ECG；x5為EGC。

圖1 最佳逐步回歸方程預(yù)測結(jié)果比較Fig. 1 True values versus predicted values from the best stepwise regression equation

2.3.3 回歸方程預(yù)測效果檢驗

為檢驗回歸方程的實際預(yù)測效果，應(yīng)用建立的主成分回歸方程和逐步多元回歸方程對10 個未知樣品茶湯滋味品質(zhì)得分進行預(yù)測，所得結(jié)果見表5。

表5 未知樣品預(yù)測效果Table 5 Prediction efficiency of unknown samples by two equations

從表5可以看出，主成分回歸方程對10 個未知樣品的茶湯滋味品質(zhì)得分進行預(yù)測，所的結(jié)果為0.958 6，RMSEP為0.401 1；當(dāng)用逐步多元回歸方程對10 個未知樣品的茶湯滋味品質(zhì)得分進行預(yù)測時，為0.974 6，RMSEP為0.391 5。逐步多元回歸方程預(yù)測結(jié)果要優(yōu)于主成分回歸方程結(jié)果，所得結(jié)論與建立回歸方程結(jié)果相符。

3 結(jié)論與討論

本實驗嘗試建立一種青磚茶茶湯品質(zhì)的化學(xué)評價方法。在感官審評評價基礎(chǔ)上，分析檢測茶湯19 種內(nèi)含成分含量，通過PCA和相關(guān)關(guān)系分析，篩選5 種與茶湯滋味品質(zhì)得分關(guān)系密切的成分，分別為EGCG、兒茶素總量、茶多酚、ECG和EGC。然后應(yīng)用主成分回歸方法和逐步回歸方法分別建立青磚茶茶湯滋味得分的擬合回歸方程，經(jīng)比較2 種方法的預(yù)測結(jié)果，以應(yīng)用5 種內(nèi)含成分建立的逐步回歸方程預(yù)測結(jié)果最佳（=0.973 5，RMSECV=0.380 7；0.968 1，RMSEP=0.400 0）。結(jié)果表明，應(yīng)用化學(xué)分析方法結(jié)合化學(xué)計量學(xué)方法預(yù)測青磚茶茶湯滋味品質(zhì)是可行的，實現(xiàn)青磚茶茶湯品質(zhì)的快速、準(zhǔn)確評價。但本實驗也存在一些不足之處，樣品采集的多樣性不夠，在實際應(yīng)用過程中，建議進一步補充采集不同渥堆程度青磚茶樣品數(shù)據(jù)，尤其是品質(zhì)較差青磚茶樣品，擴大樣品的品質(zhì)得分范圍，使得實驗數(shù)據(jù)更具有代表性和全面性，建立的擬合方程有較大的使用范圍，更加準(zhǔn)確地預(yù)測樣品品質(zhì)得分。