• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    有機(jī)成分與元素分析相結(jié)合鑒別6 種中國(guó)名茶

    2014-02-27 08:39:46靳保輝涂小珂王丙濤趙瓊暉王宏菊謝麗琪
    食品科學(xué) 2014年18期
    關(guān)鍵詞:分析模型

    陳 波,靳保輝,顏 治,涂小珂,王丙濤,趙瓊暉,王宏菊,謝麗琪

    (深圳出入境檢驗(yàn)檢疫局食品檢驗(yàn)檢疫技術(shù)中心,深圳市食品安全檢測(cè)技術(shù)研發(fā)重點(diǎn)實(shí)驗(yàn)室,廣東 深圳 518067)

    有機(jī)成分與元素分析相結(jié)合鑒別6 種中國(guó)名茶

    陳 波,靳保輝,顏 治,涂小珂,王丙濤,趙瓊暉,王宏菊,謝麗琪

    (深圳出入境檢驗(yàn)檢疫局食品檢驗(yàn)檢疫技術(shù)中心,深圳市食品安全檢測(cè)技術(shù)研發(fā)重點(diǎn)實(shí)驗(yàn)室,廣東 深圳 518067)

    利用有機(jī)成分與元素分析相結(jié)合,建立6 種中國(guó)名茶的鑒別方法。通過(guò)原產(chǎn)地購(gòu)買67 份不同產(chǎn)地的茶葉,包括綠茶(杭州龍井、蘇州碧螺春、黃山毛峰、廬山云霧)和烏龍茶(武夷巖茶、安溪鐵觀音),利用電感耦合等離子體-原子發(fā)射光譜測(cè)定茶葉中K、Ca、Mg、Al、Mn、Fe、Zn、Rb、Ti 9 種元素,利用電感耦合等離子體-質(zhì)譜法測(cè)定Pb、As、Cd、Cr、V、Ni、Co、Cu、Sr、Cs、Ba、Se 12 種元素,利用超高效液相色譜測(cè)定茶葉中10 種有機(jī)成分。通過(guò)偏最小二乘變量篩選法分別建立6 種中國(guó)名茶的元素判別模型和有機(jī)成分判別模型,并設(shè)定已知茶葉樣本的確證及未知茶葉樣本的溯源流程。實(shí)際樣品檢測(cè)結(jié)果表明,6 種中國(guó)名茶分類效果較好,通過(guò)2 種模型交叉驗(yàn)證可有效減少結(jié)果誤判。

    茶;電感耦合等離子體-原子發(fā)射光譜法;電感耦合等離子體-質(zhì)譜法;超高效液相色譜法;產(chǎn)地鑒別;交叉驗(yàn)證

    茶葉是我國(guó)重要的經(jīng)濟(jì)作物,但目前我國(guó)市場(chǎng)上茶葉商品質(zhì)量混亂,產(chǎn)地信息不明,以假充真、以次充好的現(xiàn)象普遍存在。因此,研究不同產(chǎn)地茶葉之間的差異,發(fā)展茶葉的真實(shí)屬性表征技術(shù)對(duì)于茶葉的質(zhì)量控制,保護(hù)原產(chǎn)地域產(chǎn)品和消費(fèi)者權(quán)益,市場(chǎng)的規(guī)范運(yùn)作,幫助我國(guó)茶葉走向國(guó)際市場(chǎng)都具有積極意義[1]。

    傳統(tǒng)的茶葉識(shí)別方法是感官評(píng)定法和化學(xué)方法。感官評(píng)定的結(jié)果受人為因素和外界環(huán)境的干擾很大,影響到結(jié)果的客觀性。因此,通過(guò)現(xiàn)代分析儀器尋求茶葉鑒別方法的研究一直是茶葉研究的熱點(diǎn),如核磁共振[2]、近紅外光譜[1,3-5]、高效液相色譜[6-7]、電感耦合等離子體-質(zhì)譜(inductively coupled plasma-mass spectrometry,ICP-MS)[8-10]、電感耦合等離子體-原子發(fā)射光譜(inductively coupled plasma-atomic emission spectrometry,ICPAES)[11-13]、穩(wěn)定同位素質(zhì)譜[9,14]、離子色譜[15]、氣相色譜-質(zhì)譜[16-17]、電子鼻技術(shù)[18-19]及激光誘導(dǎo)熒光[20]等?,F(xiàn)代分析技術(shù)的一個(gè)共同特點(diǎn)就是對(duì)每一個(gè)樣品測(cè)量所得到的分析結(jié)果都是多參量數(shù)據(jù),對(duì)多參量數(shù)據(jù)進(jìn)行分析,僅用肉眼判斷或簡(jiǎn)

    單的數(shù)據(jù)比較來(lái)處理是難以奏效的。因此,通常使用模式識(shí)別方法來(lái)分析數(shù)據(jù),以便最大限度提取信息。常見(jiàn)的模式識(shí)別方法包括主成分分析[2,10,17,21]、聚類分析[1-2,21]、偏最小二乘法[4]、人工神經(jīng)網(wǎng)絡(luò)[5]及多元線性回歸[22]等。

    常見(jiàn)的茶葉溯源研究是將檢測(cè)獲得數(shù)據(jù)利用一種或幾種不同的模式識(shí)別方法進(jìn)行建模分析。然而,還沒(méi)有見(jiàn)到將多種模型結(jié)合起來(lái)實(shí)現(xiàn)判別的研究報(bào)道。因此,本研究分別利用ICP-MS與ICP-AES測(cè)定茶葉中的21 種元素含量,利用超高效液相色譜(ultra high performance liquid chromatography,UPLC)法測(cè)定茶葉中10 種有機(jī)成分,將元素成分與茶葉有機(jī)成分結(jié)合在一起,利用偏最小二乘變量篩選法分別建立6 種中國(guó)名茶的元素判別模型和有機(jī)成分判別模型,將2 個(gè)模型進(jìn)行交叉驗(yàn)證,降低了誤判率,取得很好的區(qū)分效果。

    1 材料與方法

    1.1 材料與試劑

    67 份茶葉包括綠茶(杭州龍井12 份、蘇州碧螺春10 份、黃山毛峰15 份、廬山云霧15 份)和烏龍茶(武夷巖茶7 份、安溪鐵觀音8 份)均在原產(chǎn)地購(gòu)買,茶葉樣品經(jīng)高速粉碎機(jī)粉碎后室溫保存。

    甲醇(色譜純)、硝酸(質(zhì)量分?jǐn)?shù)65%,優(yōu)級(jí)純)德國(guó)Merck公司;過(guò)氧化氫(質(zhì)量分?jǐn)?shù)30%,優(yōu)級(jí)純)天津市科密歐化學(xué)試劑有限公司;三氟乙酸 美國(guó)Sigma公司;超純水 美國(guó)Millipore公司。

    標(biāo)準(zhǔn)品:元素標(biāo)準(zhǔn)儲(chǔ)備液(1 000 mg/L,使用體積分?jǐn)?shù)5%硝酸(以濃硝酸為基準(zhǔn))稀釋至適當(dāng)質(zhì)量濃度)國(guó)家標(biāo)準(zhǔn)物質(zhì)中心;咖啡因(caffeine,CAFF) 中國(guó)食品藥品檢定研究院;沒(méi)食子兒茶素(gallocatechin,GC)、表沒(méi)食子兒茶素(epigallocatechin,EGC)、表兒茶素沒(méi)食子酸酯(epigallocatechin gallate,ECG)、表兒茶素(epicatechin,EC)、表沒(méi)食子兒茶素沒(méi)食子酸酯(epigallocatechin-3-gallate,EGCG) 日本和光純藥工業(yè)株式會(huì)社;沒(méi)食子酸(gallic acid,GA) 美國(guó)Alfa Aesar公司;沒(méi)食子兒茶素沒(méi)食子酸酯(gallocatechin gallate,GCG)、兒茶素(catechin,C)、茶葉堿(theophylline,THEO) 美國(guó)Sigma公司。

    1.2 儀器與設(shè)備

    710-ES電感耦合等離子體-原子發(fā)射光譜儀 美國(guó)瓦里安公司;XSERIES 2電感耦合等離子體-質(zhì)譜儀美國(guó)熱電公司;ACQUITY超高效液相色譜儀(二極管陣列檢測(cè)器) 美國(guó)Waters公司。

    1.3 方法

    1.3.1 樣品前處理

    1.3.1.1 元素分析

    稱取0.5 g左右試樣(精確至0.001 g)于微波消解罐中,加2 mL水潤(rùn)濕樣品,加5 mL硝酸,靜置反應(yīng)2 h左右,加1 mL過(guò)氧化氫,蓋好消解罐放入微波消解儀中,根據(jù)儀器條件設(shè)定最優(yōu)程序進(jìn)行消解。冷卻后轉(zhuǎn)移定容至25 mL,待測(cè)。若樣液元素質(zhì)量濃度過(guò)高,可適當(dāng)稀釋或適當(dāng)減少取樣量。同時(shí)做試劑空白實(shí)驗(yàn)。

    1.3.1.2 有機(jī)成分分析

    參照涂小珂等[23]所建立的方法進(jìn)行操作。準(zhǔn)確稱取0.25 g茶葉樣品于50 mL離心管中,加入25 mL 90 ℃熱水,于90 ℃水浴中振蕩提取25 min后,取上清液5 mL置于50 mL容量瓶中,用二次水定容后,0.45 μm微孔濾膜過(guò)濾待測(cè)。

    1.3.2 儀器條件

    ICP-MS測(cè)量參數(shù):射頻功率1 400 W;駐留時(shí)間10 ms;測(cè)量通道3;掃描次數(shù)120 次;泵速30 r/min;冷卻氣流量14 L/min;輔助氣流量0.9 L/min;霧化氣流量0.88 L/min;其他儀器條件按儀器調(diào)諧最佳狀態(tài)時(shí)設(shè)置。

    ICP-AES測(cè)定參數(shù):功率1.20 kW;等離子氣流量15.0 L/min;輔助氣流量1.50 L/min;霧化氣壓力220 kPa;讀數(shù)時(shí)間5 s;泵速15 r/min。各元素參考分析譜線波長(zhǎng)如下:K 769.897、Fe 234.350、Al 237.312、Zn 206.200、Ca 315.887、Mn 259.372、Mg 279.800、Rb 780.026、Ti 336.122 nm。

    UPLC分析條件,參照涂小珂等[23]所建色譜方法。

    2 結(jié)果與分析

    2.1 檢測(cè)結(jié)果

    2.1.1 元素檢測(cè)

    表1 茶葉中21 種元素檢測(cè)結(jié)果Table 1 Analytical results for 21 elements in teas mg/kg

    采用ICP-AES測(cè)定所有茶葉樣本中K、Ca、Mg、Al、Mn、Fe、Zn、Rb、Ti 9 種元素,采用ICP-MS測(cè)定所有茶葉樣本Pb、As、Cd、Cr、V、Ni、Co、Cu、Sr、Cs、Ba、Se 12 種元素,共計(jì)21 種元素。結(jié)果顯示,不同元素之間含量相差很大,K含量很高,都在10 000 mg/kg以上。Mg、Ca含量在1 000~6 100 mg/kg之間,Mn、Al、Fe、Rb、Zn等金屬含量稍低,在18~2 600 mg/kg之間,其他元素含量更低,最低的As只有0.013 mg/kg。要從如此大量的數(shù)據(jù)中分析出各來(lái)源茶葉的特征信息,僅靠肉眼顯然無(wú)法實(shí)現(xiàn)。

    2.1.2 有機(jī)成分檢測(cè)

    利用UPLC測(cè)定67 份茶葉中GA、GC、THEO、CAFF、C、EC、EGC、EGCG、GCG、ECG 10 種有機(jī)成分。結(jié)果表明,CAFF與EGCG的含量最高,在300~2 600 mg/kg之間;其余有機(jī)成分含量相對(duì)低,為十幾到幾百mg/kg,其中GA的含量最低,平均含量為53.6 mg/kg。與元素檢測(cè)結(jié)果相同,要想僅靠肉眼從這么多有機(jī)成分的數(shù)據(jù)中找出規(guī)律,十分困難。

    表2 茶葉中10 種有機(jī)成分檢測(cè)結(jié)果Table 2 Analytical results for 10 organic components in teas mg/kg

    2. 2 偏最小二乘法(partial least squares,PLS)分析

    在處理自變量數(shù)很大的實(shí)際問(wèn)題時(shí),PLS得到的模型由于回歸系數(shù)個(gè)數(shù)很多而不夠簡(jiǎn)單。如果能根據(jù)PLS建模過(guò)程中的一些信息,找到一種刪除冗余變量或影響不顯著的變量的判據(jù),在保證不損失模型的預(yù)報(bào)穩(wěn)定性情況下簡(jiǎn)化模型,將對(duì)分析和處理實(shí)際問(wèn)題有很大的意義。本實(shí)驗(yàn)采用廈門大學(xué)朱爾一研制的偏最小二乘變量篩選法[24-25]建立模型,實(shí)現(xiàn)很高的判別準(zhǔn)確率。變量篩選主要通過(guò)誤差方程計(jì)算出每刪除特定變量引起的誤差變化的大小來(lái)選擇變量,重新建立預(yù)報(bào)模型。方法在刪除那些對(duì)所研究問(wèn)題影響不大的變量的同時(shí),保證模型的總誤差值的增長(zhǎng)最小。

    式中:ΔEi表示刪除第i個(gè)變量時(shí),模型總誤差的增加值,通過(guò)不斷刪除ΔE(i)值為最小的變量達(dá)到變量篩選優(yōu)化的效果;bi為線性模型第i個(gè)變量的回歸系數(shù)矢量;1i為第i個(gè)分量為1,其余分量為0的一種特殊矢量;T為正交矩陣;矩陣(TTT)-1為對(duì)角矩陣;R是PLS正交分解得到的矩陣。

    2.2.1 不同產(chǎn)地茶葉元素判別模型分析

    本實(shí)驗(yàn)一共測(cè)定67 份茶葉中的21 種元素。以67 個(gè)樣本的21 種元素含量組成(67×21)的自變量矩陣X,以6 個(gè)產(chǎn)地構(gòu)成目標(biāo)變量矩陣Y,其中矩陣Y有6 列矢量,第1列矢量y1對(duì)應(yīng)鐵觀音,如果樣本為鐵觀音,則取1,否則取0。以此類推:第2列矢量y2對(duì)應(yīng)武夷巖茶,第3列矢量y3對(duì)應(yīng)黃山毛峰,第4列矢量y4對(duì)應(yīng)西湖龍井,第5列矢量y5對(duì)應(yīng)廬山云霧,第6列矢量y6對(duì)應(yīng)碧螺春。用PLS直接對(duì)21 個(gè)元素變量數(shù)據(jù)進(jìn)行處理,根據(jù)Y的信息分別對(duì)變量矩陣X作正交分解,并采用使預(yù)測(cè)殘差平方和(predicted residual sum of squares,PRESS)值最低或接近最低的隱變量個(gè)數(shù)建立回歸模型,得到對(duì)應(yīng)于y1、y2、y3、y4、y5、y6模型PRESS值為最低值時(shí)的隱變量數(shù)分別為8、8、9、6、7、5,最終所得模型的交叉檢驗(yàn)相關(guān)系數(shù)(cross-validation correlation coefficient,CR)[24]分別為0.878、0.922、0.883、0.837、0.958、0.760。CR值越接近1則說(shuō)明模型效果越理想,通常CR值大于0.9時(shí)就可以得到十分理想的區(qū)分效果。直接PLS分析結(jié)果中除武夷巖茶與廬山云霧茶CR值大于0.9,其余4 種茶的CR值均小于0.9,碧螺春的CR值則小于0.8,因此所建模型分類效果還不能令人滿意。

    圖1 PLS篩選過(guò)程中模型CR值隨入選變量數(shù)變化關(guān)系圖(ys)Fig.1 Relationship between model CR and the number of variables in the model (y5)

    對(duì)上述21 個(gè)元素進(jìn)行變量擴(kuò)維,將各元素含量相除,引入各元素間的非線性項(xiàng),這樣自變量X矩陣中變量個(gè)數(shù)從21 個(gè)增加至274 個(gè),Y矩陣不變,對(duì)X與Y矩陣進(jìn)行PLS正交分解,采用一般PLS方法建模,所得CR值分別為0.994、0.993、0.988、0.979、0.991、0.954,說(shuō)明變量擴(kuò)維后模型的可靠性比擴(kuò)維前顯著提高。再用PLS變量篩選方法對(duì)擴(kuò)維后的大量數(shù)據(jù)進(jìn)行處理,據(jù)刪除冗余變量或影響不顯著的變量,使模型得到簡(jiǎn)化,并在該過(guò)程中不損失模型的預(yù)報(bào)穩(wěn)定性。以廬山云霧為例,目標(biāo)變量y5篩選過(guò)程中入選自變量數(shù)與模型CR值的關(guān)系見(jiàn)圖1。PLS模型信息見(jiàn)表3,CR值較變量篩選前幾 乎沒(méi)有變化。從表3可見(jiàn),對(duì)判別影響最大的參數(shù),即模型系數(shù)絕對(duì)值

    最大的5 個(gè)參數(shù)均為元素比值,如對(duì)鐵觀音判別影響最大的參數(shù)為Cd/Mg比值,其次是Se/Mn、As/Ca、Co/Mg和Co/Fe,對(duì)武夷巖茶判別影響最大的是Cd/Mg比值,對(duì)黃山毛峰判別影響最大的是Ba/Ca比值,對(duì)西湖龍井判別影響最大的是Co/Ca比值,對(duì)廬山云霧判別影響最大的是Cd/Mg比值,對(duì)碧螺春判別影響最大的是Co/Mn比值。

    表3 茶葉元素模型基本參數(shù)表Table 3 Basic parameters of element model for teas

    分別采用矢量模型y1、y2和y4、y5的擬合值作圖,如圖2所示,鐵觀音、武夷巖茶、黃山毛峰、碧螺春、廬山云霧和黃山毛峰均顯示出很好的分類效果。通常,對(duì)于產(chǎn)地模型yi計(jì)算出來(lái)的該產(chǎn)地樣本點(diǎn)估計(jì)值大于0.5,則判斷正確,否則判錯(cuò),同時(shí)模型計(jì)算出來(lái)的其他產(chǎn)地樣本的估計(jì)值小于0.5,則判斷正確,否則判錯(cuò),通過(guò)樣本總數(shù)和判別正確數(shù)來(lái)計(jì)算全部樣品實(shí)際判別的準(zhǔn)確率[25-26]。本實(shí)驗(yàn)的元素PLS模型幾乎所有茶葉對(duì)應(yīng)的模型計(jì)算結(jié)果均在0.8~1.2之間,實(shí)際判別準(zhǔn)確率為100%。

    圖2 PLS模式識(shí)別分類圖Fig.2 PLS-based classification models for 6 kinds of tea

    2.2.2 有機(jī)成分判別模型

    有機(jī)成分判別模型的建立與元素判別模型的建立相似,將10 種有機(jī)成分含量?jī)蓛上喑敫饔袡C(jī)成分間的非線性項(xiàng),這樣自變量X矩陣中變量個(gè)數(shù)從10 個(gè)增加至90 個(gè),對(duì)X與Y進(jìn)行PLS正交分解,采用PLS方法建模并進(jìn)行變量篩選,PLS模型參數(shù)見(jiàn)表4。有機(jī)成分建立判別模型的總體效果較好,但不及元素判別模型,其中西湖龍井和碧螺春的相關(guān)系數(shù)CR值小于0.9。從PLS模型系數(shù)可知,GCG含量對(duì)于鐵觀音判別的影響最大,對(duì)武夷巖茶判別影響最大的是CAFF含量,對(duì)黃山毛峰和西湖龍井判別影響最大是GC含量,廬山云霧判別影響最大是EC和GA含量,碧螺春判別影響最大是C的含量。

    表4 茶葉有機(jī)成分判別模型基本參數(shù)表Table 4 Basic parameters of organic component model for teas

    2.3 模型交叉驗(yàn)證

    本實(shí)驗(yàn)一共建立2 個(gè)模型:元素判別模型和有機(jī)成分判別模型。根據(jù)茶葉判別的需求分別采用2 種流程進(jìn)行分析,如圖3所示,將2 個(gè)判別模型結(jié)合在一起判定,取得很好的結(jié)果,減少了溯源分析的誤判率。對(duì)于在模型范圍內(nèi)的已知茶葉樣本類別的確證,只需要采用元素判別模型進(jìn)行分析,若計(jì)算結(jié)果為0.8~1.2間,則可確證。若是針對(duì)未知茶葉樣品進(jìn)行溯源分析,單獨(dú)使用元素判別模型可能會(huì)導(dǎo)致誤判,此時(shí)需結(jié)合有機(jī)成分判別模型,只有2 個(gè)模型均判別為某類茶葉時(shí)才可確認(rèn)溯源結(jié)果,這可有效減少誤判率。

    圖3 茶葉判別流程圖Fig.3 Flow chart for the discrimination of 6 kinds of tea

    為進(jìn)一步驗(yàn)證茶葉分類模型的有效性,本實(shí)驗(yàn)還搜集了16 份茶葉樣品,測(cè)定了21 種元素含量以及10 種有

    機(jī)成分含量。這16 份茶葉樣品,有8 份屬于茶葉數(shù)據(jù)庫(kù)范圍內(nèi)的品種,8 種茶葉在所建立的譜庫(kù)之外。對(duì)這16 種茶葉樣品按圖3采用2 種不同的判別方式來(lái)實(shí)現(xiàn)判別。結(jié)果表明,8 種譜庫(kù)內(nèi)的茶葉樣品通過(guò)元素判別模型準(zhǔn)確確證歸屬,8 種譜庫(kù)外茶葉樣品使用元素判別模型分析,大部分無(wú)誤判,但其中有一云南普洱茶判別顯示為黃山毛峰,經(jīng)過(guò)有機(jī)成分判別模型分析,該普洱茶在黃山毛峰的計(jì)算值為-1 288.47(表5),顯然并非黃山毛峰。普洱茶屬于全發(fā)酵茶,且越陳越好,茶多酚在發(fā)酵及存放過(guò)程中損失嚴(yán)重,然而黃山毛峰屬于綠茶,茶多酚含量較高,因此通過(guò)有機(jī)成分判別模型能夠有效的避免誤判。

    表5 普洱茶樣品的元素判別模型及有機(jī)成分判別模型分析結(jié)果Table 5 Analytical results from the element model and organic component model against an Pu’er tea sample

    3 結(jié) 論

    不同產(chǎn)地的茶葉在21 種元素含量和10 種有機(jī)成分含量上存在明顯差異。利用PLS變量篩選法,通過(guò)變量擴(kuò)維,引入變量?jī)蓛上喑?xiàng),同時(shí)通過(guò)變量篩選刪除了一些冗余變量和影響不顯著的變量,實(shí)現(xiàn)了很好的分類效果。所建立的6 種中國(guó)名茶的元素判別模型和有機(jī)成分判別模型能夠?qū)崿F(xiàn)不同茶葉的區(qū)分。設(shè)定了已知茶葉樣本的確證及未知茶葉樣本的溯源流程,通過(guò)實(shí)際樣品的驗(yàn)證,本實(shí)驗(yàn)顯示出很好的分類效果,同時(shí)有效減少了誤判率。

    [1] 陳波, 靳 保輝, 梁宏, 等. 近紅外光譜技術(shù)與聚類分析鑒別中國(guó)名茶[J].光譜實(shí)驗(yàn)室, 2012, 29(5): 2901-2905.

    [2] 陳波, 張巍, 康海寧, 等. 茶葉的1H NMR指紋圖譜研究[J]. 波譜學(xué)雜志, 2006, 23(2): 169-180.

    [3] 張龍, 潘家榮, 朱誠(chéng), 等. 基于近紅外光譜的不同發(fā)酵類型茶葉判別[J].食品科學(xué), 2012, 33(20): 149-152.

    [4] 周健, 成浩, 曾建明, 等. 基于近紅外的多相偏最小二乘模型組合分析實(shí)現(xiàn)茶葉原料品種鑒定與溯源的研究[J]. 光譜學(xué)與光譜分析, 2010, 30(10): 2650-2653.

    [5] 寧井銘, 宛曉春, 張正竹, 等. 近紅外光譜技術(shù)結(jié)合人工神經(jīng)網(wǎng)絡(luò)判別普洱茶發(fā)酵程度[J]. 農(nóng)業(yè)工程學(xué)報(bào), 2013, 29(11): 255-260.

    [6] 石明明, 唐欣, 李曉, 等. 陜西省不同產(chǎn)地綠茶中6 種活性成分含量的比較[J]. 食品科學(xué), 2013, 34(8): 232-235.

    [7] LI Yuanwang, KANG Wei, HAO Cheng, et al. Geographical tracing of Xihu Longjing tea using high performance liquid chromatography[J]. Food Chemistry, 2014, 146(1): 98-103.

    [8] 康海寧, 楊妙峰, 陳波, 等. 利用礦質(zhì)元素的測(cè)定數(shù)據(jù)判別茶葉的產(chǎn)地和品種[J]. 巖礦測(cè)試, 2006, 25(1): 22-26.

    [9] TAMARA S P, JOHN R W, KLITI G. Application of trace element and stable isotope signatures to determine the provenance of tea (Camellia sinensis) samples[J]. Food Chemistry, 2010, 118(4): 921-926.

    [10] 袁玉偉, 張永志, 付海燕. 茶葉中同位素與多元素特征及其原產(chǎn)地PCA-LDA判別研究[J]. 核農(nóng)學(xué)報(bào), 2013, 27(1): 47-55.

    [11] PEDRO L F C, MARIA J M, FERNANDO P, et al. Differentiation of tea (Camellia sinensis) varieties and their geographical origin according to their metal content[J]. Journal of Agricultural and Food Chemistry, 2001, 49(10): 4775-4779.

    [12] ANGELES M H, GUSTAVO A G. Pattern recognition procedures for differentiation of Green, Black and Oolong teas according to their metal content from inductively coupled plasma atomic emission spectrometry[J]. Talanta, 2001, 53(6): 1249-1257.

    [13] 龔自明, 王雪萍, 高士偉, 等. 礦物元素分析判別綠茶產(chǎn)地來(lái)源研究[J].四川農(nóng)業(yè)大學(xué)學(xué)報(bào), 2012, 30(4): 429-433.

    [14] RUPALI A L, DEVANATHAN A, AMZAD H L, et al. Isotope signature study of the tea samples produced at four different regions in India[J]. Analytical Methods, 2013, 5(6): 1604-1611.

    [15] 余衛(wèi)娟, 張杰, 殷月芬, 等. 離子色譜法測(cè)定茶葉中無(wú)機(jī)陰離子及其用于茶葉的分類[J]. 理化檢驗(yàn): 化學(xué)分冊(cè), 2012, 48(3): 272-279.

    [16] LIN Jie, ZHANG Pan, PAN Zhiqiang, et al. Discrimination of oolong tea (Camellia sinensis) varieties based on feature extraction and selection from aromatic profiles analysed by HS-SPME/GC-MS[J]. Food Chemistry, 2013, 144(1): 259-265.

    [17] 張雪波, 肖世青, 杜先鋒, 等. 基于主成分分析法的安溪鐵觀音香氣質(zhì)量評(píng)價(jià)模型的構(gòu)建[J]. 食品科學(xué), 2012, 33(22): 225-230.

    [18] 史波林, 趙鐳, 支瑞聰, 等. 應(yīng)用電子鼻判別西湖龍井茶香氣品質(zhì)[J].農(nóng)業(yè)工程學(xué)報(bào), 2011, 27(2): 302-306.

    [19] 張紅梅, 田輝, 何玉靜, 等. 茶葉中茶多酚含量電子鼻技術(shù)檢測(cè)模型研究[J]. 河南農(nóng)業(yè)大學(xué)學(xué)報(bào), 2012, 46(3): 302-306.

    [20] LIANG Mei, PATRIC L, MIKKEL B, et al. Tea classification and quality assessment using laser-induced fluorescence and chemometric evaluation[J]. Applied Optics, 2012, 51(7): 803-811.

    [21] 張瑞蓮, 袁海波, 尹軍峰, 等. 主成分分析與聚類分析在茶飲料湯色穩(wěn)定性評(píng)價(jià)中的應(yīng)用[J]. 茶葉科學(xué), 2010, 30(4): 287-294.

    [22] 馬燕紅, 張生萬(wàn), 李美萍, 等. 清香型白酒酒齡鑒別的方法研究[J].食品科學(xué), 2012, 33(10): 184-189.

    [23] 涂小珂, 梁宏, 吳鳳琪, 等. 超高效液相色譜法同時(shí)測(cè)定茶葉中兒茶素、酚酸及生物堿[J]. 中國(guó)衛(wèi)生檢驗(yàn)雜志, 2013, 23(16): 3167-3169.

    [24] 朱爾一, 林燕. 利用偏最小二乘法的一種變量篩選法[J]. 計(jì)算機(jī)與應(yīng)用化學(xué), 2007, 24(6): 741-745.

    [25] 朱爾一, 林燕, 莊贊勇. 偏最小二乘變量篩選法在毒品來(lái)源分析中的應(yīng)用[J]. 分析化學(xué), 2007, 35(7): 973-977.

    [26] 王丙濤, 陳波, 涂小珂, 等. 葡萄酒中元素分布與其原產(chǎn)地關(guān)系的分類模型[J]. 食品科學(xué), 2014, 35(2): 61-64.

    Discrimination of 6 Kinds of Chinese Tea by Combination of Organic Components and Multielement Analysis

    CHEN Bo, JIN Bao-hui, YAN Zhi, TU Xiao-ke, WANG Bing-tao, ZHAO Qiong-hui, WANG Hong-ju, XIE Li-qi
    (Shenzhen Key Laboratory of Detection Technology for Food Safety, Food Inspection Center of Shenzhen Entry-Exit Inspection and Quarantine Bureau, Shenzhen 518067, China)

    A discrimination method for 6 kinds of Chine se famous tea by combination of organic components and multielement analysis was established. Totally 67 different tea samples from their places of origin were collected, including green tea (Huangshan Maofeng, West Lake Longjing, Biluochun and Lushan Fragrant) and oolong tea (Wuyi Rock Tea, Ti Kuan Yin). The concentrations of 9 elements (K, Ca, Mg, Al, Mn, Fe, Zn, Rb and Ti) were determined by ICP-AES, 12 elements (Pb, As, Cd, Cr, V, Ni, Co, Cu, Sr, Cs, Ba and Se) by ICP-MS, and 10 organic compounds by UPLC. Partial least squares variable selection method was applied in the establishment of the discrimination models for elements and organic components, and a tracing procedure was set up for the confirmation of the known and unknown tea samples. The analysis of real samples showed that the 6 kinds of Chinese famous tea were well discriminated and misjudgment was reduced effectively by cross validation of the two models.

    tea; inductively coupled plasma-atomic emission spectrometry (ICP-AES); inductively coupled plasma-mass spectrometry (ICP-MS); ultra high performance liquid chromatography (UPLC); discrimination; cross-validation

    O657.33

    A

    1002-6630(2014)18-0119-05

    10.7506/spkx1002-6630-201418024

    2013-12-16

    國(guó)家質(zhì)檢總局科技計(jì)劃項(xiàng)目(2013IK180)

    陳波(1981—),男,工程師,碩士,研究方向?yàn)槭称钒踩珯z測(cè)。E-mail:15207783@qq.com

    猜你喜歡
    分析模型
    一半模型
    隱蔽失效適航要求符合性驗(yàn)證分析
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計(jì)的漸近分布
    電力系統(tǒng)不平衡分析
    電子制作(2018年18期)2018-11-14 01:48:24
    電力系統(tǒng)及其自動(dòng)化發(fā)展趨勢(shì)分析
    3D打印中的模型分割與打包
    FLUKA幾何模型到CAD幾何模型轉(zhuǎn)換方法初步研究
    中西醫(yī)結(jié)合治療抑郁癥100例分析
    在線教育與MOOC的比較分析
    黄色视频在线播放观看不卡| 免费大片黄手机在线观看| 内地一区二区视频在线| 欧美亚洲日本最大视频资源| 国产精品99久久99久久久不卡 | 欧美精品亚洲一区二区| a级毛片黄视频| 熟女人妻精品中文字幕| 日韩亚洲欧美综合| 亚洲四区av| 色婷婷久久久亚洲欧美| 99久久精品国产国产毛片| 一区二区三区免费毛片| 久久青草综合色| 久久精品久久精品一区二区三区| 在线观看美女被高潮喷水网站| 在线观看国产h片| 99九九线精品视频在线观看视频| 18在线观看网站| 在线天堂最新版资源| 欧美三级亚洲精品| 老司机影院毛片| 中文字幕人妻丝袜制服| 夫妻性生交免费视频一级片| 在线精品无人区一区二区三| 免费看不卡的av| 久久av网站| 大又大粗又爽又黄少妇毛片口| 一级爰片在线观看| 中文字幕人妻熟人妻熟丝袜美| 女的被弄到高潮叫床怎么办| 蜜桃国产av成人99| 狂野欧美白嫩少妇大欣赏| 成人无遮挡网站| 精品卡一卡二卡四卡免费| 最近手机中文字幕大全| 中文乱码字字幕精品一区二区三区| 五月天丁香电影| 国内精品宾馆在线| 亚洲av.av天堂| 26uuu在线亚洲综合色| 26uuu在线亚洲综合色| 男的添女的下面高潮视频| 亚洲欧美精品自产自拍| 大又大粗又爽又黄少妇毛片口| tube8黄色片| h视频一区二区三区| 久久久欧美国产精品| 一区二区日韩欧美中文字幕 | 卡戴珊不雅视频在线播放| 日本av免费视频播放| 日韩亚洲欧美综合| 久久 成人 亚洲| 丝袜美足系列| 日韩欧美精品免费久久| 午夜福利视频精品| 国模一区二区三区四区视频| 亚洲欧美日韩卡通动漫| 91精品伊人久久大香线蕉| videossex国产| 亚洲高清免费不卡视频| 王馨瑶露胸无遮挡在线观看| av卡一久久| 国产免费视频播放在线视频| 一级二级三级毛片免费看| 男女高潮啪啪啪动态图| av一本久久久久| 视频中文字幕在线观看| 9色porny在线观看| 久久久久人妻精品一区果冻| 亚洲国产精品一区三区| 一级二级三级毛片免费看| 国产高清有码在线观看视频| 精品国产国语对白av| 曰老女人黄片| 97精品久久久久久久久久精品| 999精品在线视频| 嘟嘟电影网在线观看| 高清毛片免费看| 天天影视国产精品| 亚洲人成网站在线观看播放| 国产探花极品一区二区| 黄色配什么色好看| 夫妻性生交免费视频一级片| 寂寞人妻少妇视频99o| 婷婷色麻豆天堂久久| 国产 一区精品| 99热网站在线观看| 国产片内射在线| 国产综合精华液| 亚洲国产最新在线播放| a 毛片基地| 丝袜在线中文字幕| 久久久久久久久久久丰满| 精品亚洲成国产av| 两个人免费观看高清视频| 亚洲,欧美,日韩| 99视频精品全部免费 在线| 亚洲激情五月婷婷啪啪| 一区在线观看完整版| 久久久久视频综合| 久久ye,这里只有精品| 精品国产国语对白av| 五月天丁香电影| 晚上一个人看的免费电影| 国产一区有黄有色的免费视频| 伦精品一区二区三区| 全区人妻精品视频| 一级,二级,三级黄色视频| 观看av在线不卡| 五月玫瑰六月丁香| 久久ye,这里只有精品| 天美传媒精品一区二区| 国内精品宾馆在线| 亚洲成人一二三区av| 超碰97精品在线观看| 亚洲av在线观看美女高潮| 欧美97在线视频| 国产欧美亚洲国产| 成人亚洲精品一区在线观看| 在线看a的网站| 国产精品成人在线| av播播在线观看一区| 九九久久精品国产亚洲av麻豆| 夫妻午夜视频| 国产成人av激情在线播放 | 成人二区视频| 国产色爽女视频免费观看| 中国美白少妇内射xxxbb| 中文字幕人妻丝袜制服| 超色免费av| 爱豆传媒免费全集在线观看| 欧美另类一区| 欧美日韩国产mv在线观看视频| 人妻少妇偷人精品九色| 老司机亚洲免费影院| 美女视频免费永久观看网站| 欧美成人午夜免费资源| 亚洲成人一二三区av| 一级爰片在线观看| 亚洲,欧美,日韩| 成人毛片60女人毛片免费| 哪个播放器可以免费观看大片| 日本猛色少妇xxxxx猛交久久| 亚洲av电影在线观看一区二区三区| 99久久综合免费| 国产亚洲一区二区精品| 亚洲第一av免费看| 亚洲精品美女久久av网站| 日本黄大片高清| 视频区图区小说| 夫妻午夜视频| 国产免费一区二区三区四区乱码| 午夜激情福利司机影院| 老熟女久久久| 一区二区日韩欧美中文字幕 | 美女脱内裤让男人舔精品视频| 伊人亚洲综合成人网| av又黄又爽大尺度在线免费看| 少妇猛男粗大的猛烈进出视频| 国产日韩欧美视频二区| av在线老鸭窝| 成人免费观看视频高清| 免费观看a级毛片全部| √禁漫天堂资源中文www| 大片免费播放器 马上看| 麻豆乱淫一区二区| 亚洲精品乱码久久久v下载方式| 赤兔流量卡办理| 欧美精品一区二区大全| xxx大片免费视频| 午夜91福利影院| 美女视频免费永久观看网站| 十八禁网站网址无遮挡| 高清黄色对白视频在线免费看| av福利片在线| 久久综合国产亚洲精品| 丝袜脚勾引网站| 国产乱人偷精品视频| 国产精品国产av在线观看| av.在线天堂| 国产精品秋霞免费鲁丝片| 99九九线精品视频在线观看视频| 亚洲精品第二区| 久久精品熟女亚洲av麻豆精品| videosex国产| 亚洲欧美精品自产自拍| 七月丁香在线播放| 日日摸夜夜添夜夜爱| 三级国产精品片| 日韩电影二区| 人人妻人人澡人人爽人人夜夜| 视频在线观看一区二区三区| 搡女人真爽免费视频火全软件| 黑人猛操日本美女一级片| 王馨瑶露胸无遮挡在线观看| 91精品国产九色| av视频免费观看在线观看| 制服丝袜香蕉在线| 亚洲国产精品专区欧美| 久久人人爽av亚洲精品天堂| 黄色视频在线播放观看不卡| 一个人免费看片子| 日韩熟女老妇一区二区性免费视频| 伊人亚洲综合成人网| 我的女老师完整版在线观看| 免费观看无遮挡的男女| 麻豆精品久久久久久蜜桃| 亚洲第一av免费看| 国产在线免费精品| 国产黄频视频在线观看| 免费av中文字幕在线| 2021少妇久久久久久久久久久| 亚洲美女搞黄在线观看| 欧美老熟妇乱子伦牲交| 寂寞人妻少妇视频99o| 久久久久久久久大av| 黄色怎么调成土黄色| 亚洲五月色婷婷综合| 精品一区二区三卡| 欧美日韩一区二区视频在线观看视频在线| av免费观看日本| 天堂8中文在线网| 99热全是精品| 大又大粗又爽又黄少妇毛片口| 夫妻性生交免费视频一级片| 99热国产这里只有精品6| 亚洲av电影在线观看一区二区三区| 简卡轻食公司| 精品一区二区三卡| 在线看a的网站| 精品午夜福利在线看| 18禁在线无遮挡免费观看视频| 99热这里只有精品一区| 日韩,欧美,国产一区二区三区| 婷婷色麻豆天堂久久| 国产成人精品婷婷| 日本爱情动作片www.在线观看| 少妇熟女欧美另类| 午夜视频国产福利| 你懂的网址亚洲精品在线观看| 最黄视频免费看| 亚洲色图综合在线观看| 国产精品久久久久久久久免| 亚洲国产日韩一区二区| 国产精品一国产av| 中文字幕最新亚洲高清| 国产精品久久久久久精品电影小说| av黄色大香蕉| 波野结衣二区三区在线| 国产欧美日韩综合在线一区二区| 欧美性感艳星| 丰满饥渴人妻一区二区三| 80岁老熟妇乱子伦牲交| 欧美xxⅹ黑人| 亚州av有码| 最近手机中文字幕大全| 国产亚洲精品第一综合不卡 | 国产毛片在线视频| 国产午夜精品久久久久久一区二区三区| 亚洲综合色网址| 热99国产精品久久久久久7| 中文字幕久久专区| 永久网站在线| 一本—道久久a久久精品蜜桃钙片| 最近手机中文字幕大全| 精品亚洲成国产av| 亚洲第一区二区三区不卡| 久久精品国产自在天天线| 日韩成人av中文字幕在线观看| 嘟嘟电影网在线观看| 日韩,欧美,国产一区二区三区| 最黄视频免费看| 黄色怎么调成土黄色| 麻豆精品久久久久久蜜桃| 999精品在线视频| 色哟哟·www| 亚洲国产精品999| 亚洲国产日韩一区二区| 青青草视频在线视频观看| videossex国产| 国产无遮挡羞羞视频在线观看| 日韩,欧美,国产一区二区三区| 一区二区av电影网| 黄色怎么调成土黄色| 国产一级毛片在线| 黄色配什么色好看| 久久久久精品久久久久真实原创| 国产成人免费无遮挡视频| 精品亚洲成a人片在线观看| 亚洲国产精品一区二区三区在线| 一边亲一边摸免费视频| 热99国产精品久久久久久7| 日韩一区二区视频免费看| 亚洲精品久久久久久婷婷小说| 免费av中文字幕在线| 蜜桃国产av成人99| 免费大片黄手机在线观看| 极品人妻少妇av视频| 日韩亚洲欧美综合| 大又大粗又爽又黄少妇毛片口| 日韩一本色道免费dvd| 一级a做视频免费观看| 97精品久久久久久久久久精品| 国产精品国产av在线观看| 中文欧美无线码| 久久97久久精品| 亚洲精品国产av蜜桃| 大陆偷拍与自拍| 女人精品久久久久毛片| 亚洲精品国产av成人精品| 青春草视频在线免费观看| 欧美亚洲 丝袜 人妻 在线| 黄片播放在线免费| 久久久久久久国产电影| 亚洲成色77777| 中国三级夫妇交换| 热99国产精品久久久久久7| av视频免费观看在线观看| av免费观看日本| 人妻人人澡人人爽人人| 欧美97在线视频| 蜜臀久久99精品久久宅男| 免费观看无遮挡的男女| 涩涩av久久男人的天堂| 亚洲欧美一区二区三区黑人 | 另类亚洲欧美激情| 免费人妻精品一区二区三区视频| 人人澡人人妻人| 老熟女久久久| 午夜福利在线观看免费完整高清在| 男女边吃奶边做爰视频| 人妻系列 视频| 亚洲在久久综合| 亚洲怡红院男人天堂| 久久国产亚洲av麻豆专区| 国产深夜福利视频在线观看| 亚洲欧美精品自产自拍| 国产男人的电影天堂91| 人妻系列 视频| 九草在线视频观看| 国产精品国产三级国产av玫瑰| 丝袜脚勾引网站| 日本黄色日本黄色录像| 色吧在线观看| 狂野欧美白嫩少妇大欣赏| 亚洲精品美女久久av网站| 美女国产视频在线观看| 美女国产高潮福利片在线看| 国产成人午夜福利电影在线观看| 美女主播在线视频| 一本久久精品| 亚洲精品一二三| 国产欧美日韩综合在线一区二区| 春色校园在线视频观看| 久久亚洲国产成人精品v| 国产精品秋霞免费鲁丝片| 中文字幕精品免费在线观看视频 | 国产日韩欧美亚洲二区| 交换朋友夫妻互换小说| 最近中文字幕高清免费大全6| 极品少妇高潮喷水抽搐| 国产黄色免费在线视频| 午夜福利视频精品| 欧美 日韩 精品 国产| 欧美日韩综合久久久久久| 久久久国产一区二区| 飞空精品影院首页| 精品国产露脸久久av麻豆| 91久久精品国产一区二区三区| 内地一区二区视频在线| av网站免费在线观看视频| 黄片播放在线免费| 精品亚洲成a人片在线观看| 日韩人妻高清精品专区| 亚洲色图综合在线观看| 久久狼人影院| 国产精品 国内视频| 日日摸夜夜添夜夜爱| 黄色欧美视频在线观看| 国产黄频视频在线观看| 人妻人人澡人人爽人人| 秋霞伦理黄片| 少妇被粗大猛烈的视频| 我要看黄色一级片免费的| 亚洲欧美成人精品一区二区| 亚洲怡红院男人天堂| 最近的中文字幕免费完整| 色94色欧美一区二区| 国产男女内射视频| 三上悠亚av全集在线观看| 国产免费一级a男人的天堂| 啦啦啦中文免费视频观看日本| videos熟女内射| 国产黄频视频在线观看| 日韩欧美精品免费久久| 日韩在线高清观看一区二区三区| 免费播放大片免费观看视频在线观看| a级毛片黄视频| 精品酒店卫生间| 日本猛色少妇xxxxx猛交久久| 午夜精品国产一区二区电影| 午夜福利网站1000一区二区三区| 亚洲,一卡二卡三卡| 欧美97在线视频| 欧美人与性动交α欧美精品济南到 | 爱豆传媒免费全集在线观看| 三级国产精品片| 中文字幕人妻丝袜制服| 黑人欧美特级aaaaaa片| 亚洲欧美色中文字幕在线| 日本色播在线视频| 久久久久精品性色| 91精品三级在线观看| 一级片'在线观看视频| 成年人午夜在线观看视频| 亚洲精华国产精华液的使用体验| 极品人妻少妇av视频| 国产乱来视频区| 日本免费在线观看一区| 黄片播放在线免费| 久久久久久久亚洲中文字幕| 哪个播放器可以免费观看大片| 国产亚洲精品久久久com| 久久久久国产网址| 国产亚洲午夜精品一区二区久久| 亚洲精品乱码久久久v下载方式| 亚洲av国产av综合av卡| 婷婷色综合www| 国产精品99久久久久久久久| 中文精品一卡2卡3卡4更新| 91精品国产九色| 国产精品 国内视频| 大片免费播放器 马上看| 国产成人午夜福利电影在线观看| 我的女老师完整版在线观看| 免费看av在线观看网站| 成人毛片60女人毛片免费| av国产久精品久网站免费入址| 欧美激情国产日韩精品一区| 色哟哟·www| av在线app专区| 少妇被粗大猛烈的视频| 亚洲欧美一区二区三区国产| 黑丝袜美女国产一区| 久久久久久久精品精品| 久久精品熟女亚洲av麻豆精品| 成人国产麻豆网| 韩国av在线不卡| 你懂的网址亚洲精品在线观看| 少妇猛男粗大的猛烈进出视频| 精品人妻在线不人妻| 草草在线视频免费看| 欧美亚洲日本最大视频资源| 精品人妻偷拍中文字幕| 国模一区二区三区四区视频| 亚洲久久久国产精品| 国产精品免费大片| 伊人久久国产一区二区| 一区在线观看完整版| 在线免费观看不下载黄p国产| 国产精品99久久久久久久久| 欧美日韩在线观看h| 高清欧美精品videossex| 色94色欧美一区二区| 中文字幕最新亚洲高清| 国产免费又黄又爽又色| 成人黄色视频免费在线看| 少妇人妻精品综合一区二区| 精品熟女少妇av免费看| 日日摸夜夜添夜夜爱| 欧美激情 高清一区二区三区| 午夜激情久久久久久久| 蜜桃久久精品国产亚洲av| 国产精品一国产av| 亚洲欧美日韩另类电影网站| 亚洲欧美精品自产自拍| 好男人视频免费观看在线| 中国美白少妇内射xxxbb| xxx大片免费视频| 精品亚洲成国产av| 久久久久精品久久久久真实原创| 久久久精品区二区三区| 91午夜精品亚洲一区二区三区| 色婷婷久久久亚洲欧美| 高清午夜精品一区二区三区| 91精品国产国语对白视频| 91精品三级在线观看| 9色porny在线观看| 欧美变态另类bdsm刘玥| 一边亲一边摸免费视频| 亚洲欧美色中文字幕在线| 日日撸夜夜添| 免费黄网站久久成人精品| 91aial.com中文字幕在线观看| 国产精品一区www在线观看| 一级,二级,三级黄色视频| 婷婷色麻豆天堂久久| 最近2019中文字幕mv第一页| 我的老师免费观看完整版| 2021少妇久久久久久久久久久| 国产高清不卡午夜福利| 午夜福利视频精品| 人妻制服诱惑在线中文字幕| 天堂8中文在线网| 在线精品无人区一区二区三| 国产视频首页在线观看| 久久久欧美国产精品| 国产精品欧美亚洲77777| 国产亚洲一区二区精品| 亚洲一区二区三区欧美精品| 欧美精品亚洲一区二区| 91久久精品电影网| 成年人午夜在线观看视频| 永久免费av网站大全| 国产黄色视频一区二区在线观看| 久久久午夜欧美精品| 性色av一级| 亚洲色图 男人天堂 中文字幕 | 各种免费的搞黄视频| 久久狼人影院| 伦理电影免费视频| 亚洲国产精品999| 久久综合国产亚洲精品| 亚洲欧美成人精品一区二区| 搡老乐熟女国产| 亚洲精华国产精华液的使用体验| 国产片内射在线| 免费看不卡的av| 丰满迷人的少妇在线观看| 久久婷婷青草| 欧美另类一区| 少妇的逼水好多| 免费av中文字幕在线| 久久国产精品男人的天堂亚洲 | 日本vs欧美在线观看视频| 一级黄片播放器| 三级国产精品欧美在线观看| 一级毛片我不卡| 欧美老熟妇乱子伦牲交| xxx大片免费视频| 久久狼人影院| 午夜久久久在线观看| av专区在线播放| 中文精品一卡2卡3卡4更新| 精品一区二区三卡| 九色成人免费人妻av| 中文字幕久久专区| 日本猛色少妇xxxxx猛交久久| 人妻 亚洲 视频| 男女国产视频网站| 熟女av电影| 18禁在线播放成人免费| 午夜福利,免费看| 日本爱情动作片www.在线观看| 亚洲国产日韩一区二区| 久久国产亚洲av麻豆专区| 国产男女超爽视频在线观看| 免费观看的影片在线观看| 免费观看在线日韩| 久久久久精品性色| 国产永久视频网站| 最近2019中文字幕mv第一页| 亚洲一区二区三区欧美精品| 国产视频内射| 国产 一区精品| 久久精品国产亚洲av天美| 亚洲欧美成人综合另类久久久| 黑人巨大精品欧美一区二区蜜桃 | 亚洲精品乱久久久久久| 亚洲少妇的诱惑av| 精品午夜福利在线看| 80岁老熟妇乱子伦牲交| 国产精品嫩草影院av在线观看| 一级毛片黄色毛片免费观看视频| 2021少妇久久久久久久久久久| 天天躁夜夜躁狠狠久久av| 在线免费观看不下载黄p国产| 黄片无遮挡物在线观看| 女的被弄到高潮叫床怎么办| 3wmmmm亚洲av在线观看| 丝瓜视频免费看黄片| 我的女老师完整版在线观看| 精品亚洲成国产av| 国产国拍精品亚洲av在线观看| 少妇被粗大的猛进出69影院 | 国产一区二区在线观看av| 亚洲国产毛片av蜜桃av| 国产国拍精品亚洲av在线观看| 在线播放无遮挡| 国精品久久久久久国模美| 国产男女超爽视频在线观看| 99久国产av精品国产电影| 久久热精品热| 免费看av在线观看网站| 久久久久精品久久久久真实原创| 国产亚洲精品第一综合不卡 | 国产精品三级大全| 亚洲国产精品999| 久久韩国三级中文字幕| 国产有黄有色有爽视频| 曰老女人黄片| 制服诱惑二区| 自线自在国产av| 欧美精品国产亚洲| 中国美白少妇内射xxxbb| 亚洲精品av麻豆狂野| 老司机影院毛片| 国产 精品1| 国产熟女午夜一区二区三区 | 国产成人91sexporn| 亚洲丝袜综合中文字幕| 免费大片18禁| 成人亚洲欧美一区二区av| 下体分泌物呈黄色|