王宇童 唐卿雁 邵金良 林 奇 杜麗娟

（1. 云南省農(nóng)業(yè)科學(xué)院質(zhì)量標(biāo)準(zhǔn)與檢測(cè)技術(shù)研究所， 昆明 650205； 2. 云南省農(nóng)業(yè)大學(xué)食品科學(xué)技術(shù)學(xué)院，昆明 650500）

普洱茶是云南省的地理標(biāo)志性產(chǎn)品， 其是以云南大葉種曬青毛茶為原料， 經(jīng)過(guò)鮮葉攤放、 殺青、揉捻、 解塊、 日光干燥、 精制、 蒸壓成型等工藝制成的茶葉。 普洱生茶的品質(zhì)特征為外形色澤墨綠，香氣清純持久， 滋味濃厚回甘， 湯色綠黃清亮， 葉底肥厚黃綠[1]。 作為全球3 大天然飲料之一， 普洱茶不僅富含茶多酚、 咖啡堿、 氨基酸等多種營(yíng)養(yǎng)成分， 還含有多種人體所需的礦質(zhì)元素， 具有抗癌[2]、抗氧化[3]、 降血糖[4]、 降血脂[5]、 抗動(dòng)脈粥樣硬化[6]等多種功效。 普洱生茶具有更強(qiáng)的保護(hù)酒精性胃損傷[7]、 抗氧化活性、 降血糖血脂等功效[8]， 因而受到廣大消費(fèi)者的青睞。 但在普洱茶產(chǎn)業(yè)不斷發(fā)展的同時(shí)， 一些商家采用產(chǎn)地造假、 以次充好、 品牌冒用等方式用劣質(zhì)普洱茶冒充優(yōu)質(zhì)普洱茶進(jìn)行售賣，不僅侵犯了消費(fèi)者合法權(quán)益， 也阻礙了普洱茶產(chǎn)業(yè)的健康發(fā)展[9]。 因此， 解決普洱茶產(chǎn)地溯源問(wèn)題亟不可待。

近年來(lái)， 農(nóng)產(chǎn)品產(chǎn)地溯源發(fā)展迅速， 主要采用不同產(chǎn)地農(nóng)產(chǎn)品的礦質(zhì)元素、 有機(jī)元素、 穩(wěn)定同位素等成分的差異結(jié)合化學(xué)計(jì)量法來(lái)進(jìn)行農(nóng)產(chǎn)品產(chǎn)地溯源。 礦質(zhì)元素能很好地反映植物、 土壤間的遷移， ZHAO 和YANG[10]對(duì)稀土元素在茶葉地理溯源中的適宜性進(jìn)行研究， 證明了茶葉中稀土元素含量與種源土壤存在顯著相關(guān)性。 林昕等[11]基于稀土元素含量差異性對(duì)普洱茶古樹(shù)茶和臺(tái)地茶的產(chǎn)地和身份進(jìn)行判別， 建立的判別模型回判驗(yàn)證正確率超過(guò)90%， 表明利用稀土元素指紋技術(shù)來(lái)判別普洱古樹(shù)茶和臺(tái)地茶的分類方法是可行的。 石伊凡等[12]利用有機(jī)元素及礦質(zhì)元素對(duì)3 個(gè)產(chǎn)區(qū)的龍井茶進(jìn)行溯源， 判別準(zhǔn)確度達(dá)80%以上。 HUANG 等[13]通過(guò)測(cè)定羅漢果中14 種礦質(zhì)元素結(jié)合線性判別分析和支持向量機(jī)建立溯源模型， 判別成功率達(dá)到100%。

同位素指紋與作物生長(zhǎng)環(huán)境密切相關(guān)， 不同產(chǎn)地的作物穩(wěn)定同位素比值會(huì)有所差異， 成為作物天然的身份標(biāo)簽[14]。 聶晶等[15]通過(guò)測(cè)定輕穩(wěn)定同位素比值、 礦物元素含量對(duì)不同產(chǎn)地的綠茶進(jìn)行識(shí)別分析， 總體判別準(zhǔn)確度在95%以上。 BROMBIN 等[16]對(duì)意大利北部?jī)蓚€(gè)相距不到2 km 果園里蘋果的多穩(wěn)定同位素比值進(jìn)行分析， 結(jié)果表明， 即使產(chǎn)品種植的距離只有幾千米， 多穩(wěn)定同位素分析也可以對(duì)其進(jìn)行追溯。 張欣昕等[17]利用礦質(zhì)元素及穩(wěn)定同位素對(duì)馬鈴薯進(jìn)行產(chǎn)地辨別， 結(jié)果表明， 礦質(zhì)元素結(jié)合同位素比值進(jìn)行產(chǎn)地判別分析的正確率高于礦質(zhì)元素含量和同位素比值單獨(dú)判別正確率。 LIANG等[18]通過(guò)測(cè)定碳 （C）、 氮 （N）、 氫 （H）、 氧 （O）的穩(wěn)定同位素比值結(jié)合化學(xué)計(jì)量學(xué)方法， 對(duì)3 個(gè)三七產(chǎn)地的三七主根建立了產(chǎn)地判別模型， 不同計(jì)量方法對(duì)三七產(chǎn)地的正確判別率平均達(dá)到84.09%。

針對(duì)茶葉產(chǎn)地溯源研究， 現(xiàn)有研究多是通過(guò)茶葉本身所含的礦質(zhì)元素、 穩(wěn)定同位素、 化學(xué)成分等生化指標(biāo)進(jìn)行溯源分析， 而對(duì)土壤-茶葉中的穩(wěn)定同位素的相關(guān)關(guān)系， 以及對(duì)產(chǎn)地溯源影響等方面的研究較少。 鑒于此， 本研究在基于以往的研究模式上， 將土壤、 茶葉中的穩(wěn)定同位素進(jìn)行相關(guān)性分析[19]。 通過(guò)測(cè)定茶樣中穩(wěn)定同位素C、 N 比值和礦質(zhì)元素， 分析不同產(chǎn)地普洱茶中穩(wěn)定同位素比值、礦質(zhì)元素和有機(jī)成分的差異， 結(jié)合單因素方差分析、 正交偏最小二乘判別分析 （Orthogonal partial least squares-discriminant analysis， OPLS-DA）、 判別分析和神經(jīng)網(wǎng)絡(luò)-多層感知器分析構(gòu)建普洱茶在小尺度空間的產(chǎn)地溯源模型， 為云南高原特色地理標(biāo)志產(chǎn)品普洱茶產(chǎn)地溯源技術(shù)提供理論基礎(chǔ)。

一、 材料與方法

（一） 材料、 試劑與儀器

1.實(shí)驗(yàn)樣品。 以云南省臨滄市、 普洱市及西雙版納傣族自治州3 個(gè)地級(jí)市下的5 個(gè)縣（區(qū)） 為普洱茶生茶溯源產(chǎn)地， 共采集30 份普洱茶生茶茶樣及其對(duì)應(yīng)產(chǎn)地的土壤樣本 （見(jiàn)表1）。 普洱茶生茶樣品及土壤樣品均采集于2021 年11-12 月， 所有普洱茶生茶原料均為單一產(chǎn)地的古樹(shù)茶， 沒(méi)有多產(chǎn)地復(fù)配情況。 采集的土壤樣品深度均為60 cm。

表1 采集地相關(guān)信息和樣品數(shù)

2.實(shí)驗(yàn)試劑。 咖啡堿標(biāo)準(zhǔn)溶液、 茶氨酸標(biāo)準(zhǔn)溶液（純度＞99%， 國(guó)家標(biāo)準(zhǔn)物質(zhì)中心）； 硝酸（優(yōu)級(jí)純， 蘇州晶瑞化學(xué)股份有限公司）； GBW100494 小麥成分分析標(biāo)準(zhǔn)物質(zhì)（鋼研納克檢測(cè)技術(shù)股份有限公司） ； 多元素標(biāo)準(zhǔn)溶液SGB -YYA23001I、SGB-YYA040011 （100 mg/L， 上海市計(jì)量測(cè)試技術(shù)研究院）； 元素混合溶液標(biāo)準(zhǔn)物質(zhì)NIMRM5307（10 mg/L， 中國(guó)計(jì)量科學(xué)研究院）。

3. 實(shí)驗(yàn)儀器。 SPECORD210Plus 紫外可見(jiàn)分光光度計(jì) （德國(guó)耶拿分析儀器股份公司）； U1timate 3000-IAP RQ 電感耦合等離子體質(zhì)譜儀 （美國(guó)賽默飛世爾科技公司）； I-TOUCH 石墨爐消解儀（北京萊伯泰科公司）； OPTIMA8000 電感耦合等離子體發(fā)射光譜儀 （美國(guó)PE 公司）； EA-VisION 穩(wěn)定同位素質(zhì)譜儀 ［艾力蒙塔貿(mào)易 （上海） 有限公司］； DHG-9203A 電熱鼓風(fēng)干燥箱 （上海一恒科技有限公司）； TGL-10B 高速臺(tái)式離心機(jī) （上海安亭科學(xué)儀器廠）。

4.樣品前處理。 將采集到的普洱茶樣品用粉碎機(jī)粉碎均勻過(guò)篩， 75℃烘干至恒重， 放入密封袋，密封避光保存。

（二） 穩(wěn)定同位素比值的測(cè)定參考胡玲等[20]的方法， 取約0.8 mg 待測(cè)樣品包入錫紙杯， 放置在同位素質(zhì)譜儀自動(dòng)進(jìn)盤中， C、 N、 硫（S） 元素通過(guò)高溫燃燒后再經(jīng)過(guò)銅絲還原生成CO2、 N2和SO2， 之后載入質(zhì)譜檢測(cè)。 穩(wěn)定同位素質(zhì)譜儀具體參數(shù)： 燃燒管（Comb tube） 溫度為1150℃， 還原管 （Reduct tube） 溫度為850℃， TCD 流速為230 mL/min， 氦氣流速為230 mL/min。

（三） 礦質(zhì)元素的測(cè)定

1.樣品消解。 采用硝酸結(jié)合高氯酸對(duì)樣品、 國(guó)家標(biāo)準(zhǔn)物質(zhì)GBW100494 及試劑空白進(jìn)行消解。 消解過(guò)程為稱取0.3 g 樣品 （精確到0.001 g） 于聚四氟乙烯消解管， 加入10 mL 硝酸和1 mL 高氯酸，加蓋放置于石墨爐消解儀中， 按設(shè)定的溫度分階段消解。 消解結(jié)束后加入0.65 mL 硝酸 （1∶1，V∶V）， 倒入25 mL 容量瓶中， 使用超純水清洗3次消解管， 倒入容量瓶， 定容至25 mL。

2. 微量元素的測(cè)定。 參考胡玲等[20]的方法， 使用電感耦合等離子體質(zhì)譜儀進(jìn)行測(cè)定。 工作參數(shù)：分析模式為碰撞反應(yīng)池； 蠕動(dòng)泵清洗時(shí)間和提升時(shí)間分別為20 s 和30 s； 射頻功率1500 W； 等離子體氣流量、 載氣流量、 輔助氣流量、 氦氣流量分別為15、 0.8、 1.2 L/min 和4.5 mL/min； 采樣深度9.0 mm， 采樣模式為Spectrum； 重復(fù)采集數(shù)據(jù)3 次。

3. 常量元素的測(cè)定。 參考魏新軍等[21]的方法，使用電感耦合等離子體發(fā)射光譜儀進(jìn)行測(cè)定。 工作參數(shù)： 功率1300 W； 等離子體流量、 輔助氣流量和霧化氣氣體流量分別為20、 1 和0.5 L/min； 重復(fù)采集數(shù)據(jù)3 次。

（四） 有機(jī)成分的測(cè)定

1.游離氨基酸的測(cè)定。 游離氨基酸總量的測(cè)定參照GB/T 8314-2013 《茶 游離氨基酸總量的測(cè)定》。

2. 茶多酚的測(cè)定。 茶多酚含量測(cè)定參照GB/T 8313-2018 《茶葉中茶多酚和兒茶素類含量的檢測(cè)方法》。

3. 咖啡堿的測(cè)定。 咖啡堿含量測(cè)定參照GB/T 8312-2013 《茶 咖啡堿測(cè)定》。

（五） 數(shù)據(jù)處理數(shù)據(jù)均為3 次重復(fù)實(shí)驗(yàn)結(jié)果，采用SPSS 26.0 統(tǒng)計(jì)分析軟件進(jìn)行單因素方差分析（選用Duncan 進(jìn)行多重比較）、 Fisher 判別分析和多層感知器分析； 采用SIMCA 14.1 統(tǒng)計(jì)分析軟件進(jìn)行OPLS-DA； 作圖采用Origin 2022。

二、 結(jié)果與分析

（一） 穩(wěn)定同位素比值測(cè)定結(jié)果分析

1.穩(wěn)定同位素比值結(jié)果分析。 對(duì)5 個(gè)產(chǎn)地普洱茶生茶和對(duì)應(yīng)土壤的穩(wěn)定同位素比值進(jìn)行單因素方差分析， 結(jié)果顯示， 茶葉中δ15N、 δ13C 值均存在組間顯著差異 （p＜0.05）（見(jiàn)表2）。 原因可能是δ15N的組成差異不僅僅受到生態(tài)環(huán)境的影響， 很大程度上還受不同類型的肥料影響[20]， 而植物中的碳同位素會(huì)受到不同氣候和環(huán)境因子的影響[21]， 因此不同產(chǎn)地茶葉中的δ15N、 δ13C 值存在一定的差異性。土壤中的δ34S、 δ15N、 δ13C 也存在組間顯著差異性（p＜0.05）。

表2 不同產(chǎn)地普洱生茶、 土壤中的穩(wěn)定同位素比值（‰）

2.穩(wěn)定同位素比值相關(guān)性分析。 對(duì)茶葉、 土壤中的穩(wěn)定同位素比值進(jìn)行Pearson 相關(guān)性分析， 結(jié)果表明， 茶樣中δ13C 值與土壤中δ13C 值存在極顯著負(fù)相關(guān)性 （見(jiàn)圖1）， 與土壤中δ34S 值存在顯著的正相關(guān)。 茶葉中δ15N 值與土壤中δ15N 值存在極顯著正相關(guān)， 與土壤中δ13C 值存在顯著正相關(guān)，與土壤中δ34S 值存在極顯著的負(fù)相關(guān)性。 原因可能是茶葉中的氮主要來(lái)源于根部吸收的土壤養(yǎng)分， 不同外源氮的同位素會(huì)直接影響茶樹(shù)體內(nèi)氮同位素的變化[22]。 同時(shí)， 土壤中δ15N、 δ13C、 δ34S 值也存在一定相關(guān)性， 土壤中δ15N 值與土壤中δ13C 值存在極顯著正相關(guān)， 而土壤中δ34S 值與土壤中δ15N、δ13C 值存極顯著負(fù)相關(guān)。 這說(shuō)明茶葉與土壤間的同位素比值存在一定的相關(guān)性。

圖1 不同產(chǎn)地普洱生茶及對(duì)應(yīng)土壤中穩(wěn)定同位素比值的相關(guān)性熱圖

（二） 礦質(zhì)元素測(cè)定結(jié)果分析

1.常量元素結(jié)果分析。 比較各產(chǎn)區(qū)普洱生茶的6 種常量元素 （Ca、 K、 Mg、 P、 Na、 S）， 其中K、P、 Na 元素不存在顯著差異， Ca、 Mg、 S 元素均存在組間顯著性差異 （p＜0.05）（見(jiàn)表3）。 雙江自治縣普洱生茶的Ca、 Mg、 S 元素含量均低于5 地的平均值； 臨翔區(qū)產(chǎn)地普洱生茶的Ca 含量明顯低于其他4 地； 瀾滄自治縣普洱生茶的Ca、 Mg 含量均高于其余4 地， 勐?？h、 勐臘縣兩地生茶的S 含量均值分別為2149、 2011 mg/kg， 高于其他3 地。

表3 不同產(chǎn)地普洱生茶中的常量元素含量 （mg/kg）

2. 微量元素結(jié)果分析。 比較不同產(chǎn)區(qū)茶樣的10 種微量元素 （Se、 Mo、 Mn、 Sn、 Cu、 Fe、 Zn、Ba、 Al、 Kr） 含量， 除了Mo 不存在顯著差異性，其余9 種微量元素 （Se、 Mn、 Sn、 Cu、 Fe、 Zn、Ba、 Al、 Kr） 含量均具有組間顯著差異性（p＜0.05）（見(jiàn)表4）。 其中， 雙江自治縣產(chǎn)地茶樣中Mn、Zn 含量明顯低于其余4 地； 臨翔區(qū)產(chǎn)地茶樣中Mn、 Fe、 Al 含量明顯高于其余4 地； 勐?？h產(chǎn)地茶樣中的Sn、 Ba 明顯高于其余4 地； 勐臘縣產(chǎn)地茶樣中Se、 Mo、 Fe、 Ba、 Al 含量低于其余4 地。

表4 不同產(chǎn)地普洱生茶中的微量元素含量（mg/kg）

3.稀土元素結(jié)果分析。 對(duì)不同產(chǎn)地普洱生茶的16 種稀土元素含量進(jìn)行比較， 結(jié)果顯示， 不同產(chǎn)地茶樣中16 種稀土元素含量存在組間顯著差異 （p＜0.05）（見(jiàn)表5）， 表明稀土元素與不同產(chǎn)地存在一定的相關(guān)性。

表5 不同產(chǎn)地普洱生茶中的稀土元素含量 （μg/kg）

（三） 有機(jī)成分測(cè)定結(jié)果分析不同產(chǎn)地茶樣中的茶多酚、 游離氨基酸、 咖啡堿含量均呈現(xiàn)組間顯著性差異（p＜0.05）（見(jiàn)表6）。 其中， 勐臘縣產(chǎn)地普洱生茶茶多酚含量低于其他4 個(gè)產(chǎn)地； 臨翔區(qū)產(chǎn)地普洱生茶游離氨基酸含量明顯均低于其他4 個(gè)產(chǎn)地； 5 個(gè)產(chǎn)地的茶樣中， 勐?？h的咖啡堿平均含量最高， 瀾滄自治縣的游離氨基酸平均含量最高。 這說(shuō)明不同產(chǎn)地有機(jī)成分含量具有一定差異， 對(duì)產(chǎn)地溯源也具有一定的作用。

表6 不同產(chǎn)地普洱生茶中的有機(jī)成分含量（%）

（四）普洱生茶產(chǎn)地溯源分析

1.OPLS-DA。 將穩(wěn)定同位素及礦質(zhì)元素在單因素方差分析結(jié)果中存在顯著差異的30 項(xiàng)指標(biāo)進(jìn)行OPLS-DA， 結(jié)果顯示， 自變量擬合指數(shù)R2x＝0.894， 因變量擬合指數(shù)R2y＝0.842， 模型預(yù)測(cè)指數(shù)Q2＝0.747， 即對(duì)5 個(gè)不同產(chǎn)地茶樣的預(yù)測(cè)率為74.7%，R2和Q2均大于0.5 表示模型擬合的結(jié)果可接受[23]。 OPLS-DA 得分圖 （見(jiàn)圖2） 中可以看出不同產(chǎn)地的樣品分離效果一般， 臨翔區(qū)和勐臘縣茶樣相對(duì)聚集， 可被明顯區(qū)分開(kāi)； 雙江自治縣與瀾滄自治縣樣品雖然能被區(qū)分開(kāi)， 但相距很短； 勐?？h茶樣較為分散， 有個(gè)別茶樣靠近雙江自治縣茶樣。 對(duì)OPLS-DA 模型進(jìn)行200 次置換檢驗(yàn)， 結(jié)果顯示，Q2回歸線與縱軸的相交點(diǎn)小于0 （見(jiàn)圖3）， 說(shuō)明構(gòu)建的OPLS-DA 模型沒(méi)有出現(xiàn)過(guò)擬合，模型的驗(yàn)證有效。 變量重要性投影 （Variable importance in projection， VIP） 值見(jiàn)圖4， VIP 值越大代表對(duì)于區(qū)分組別的貢獻(xiàn)越大， 因?yàn)樗凶兞縑IP值的平方和為1， 所以常將VIP＞1 的指標(biāo)作為判定區(qū)分組別的主要差異指標(biāo)[24]。 在該模型中共有11 個(gè)指標(biāo) （Cu、 Fe、 Mn、 Ca、 Mg、 Zn、 Al、 S、 Kr、Ce、 δ13C） 的VIP 值大于1， 表示在此模型中這11種成分是區(qū)分不同產(chǎn)地普洱生茶的差異物質(zhì)。

圖2 不同產(chǎn)地普洱生茶樣品的OPLS－DA 得分圖

圖3 OPLS－DA 模型200 次迭代計(jì)算驗(yàn)證圖

圖4 存在顯著組間差異的30 項(xiàng)指標(biāo)的VIP 圖

2.判別分析。 根據(jù)指標(biāo)（p＜0.05、 VIP 值＞1 的標(biāo)準(zhǔn)， 將上述篩選出的11 種差異物質(zhì)使用SPSS進(jìn)行判別分析， 得到4 個(gè)判別方程式［見(jiàn)公式（1）～（4）]和5 個(gè)組質(zhì)心函數(shù)（見(jiàn)表7）。

表7 組質(zhì)心函數(shù)

以第1、 第2 判別函數(shù)為橫縱坐標(biāo)， 根據(jù)組質(zhì)心處的函數(shù)繪制散點(diǎn)圖 （見(jiàn)圖5）， 可以看出5 個(gè)產(chǎn)地的茶樣聚集性效果較好， 只有雙江自治縣和勐?？h出現(xiàn)輕微重疊的現(xiàn)象， 但雙江自治縣、 瀾滄自治縣、 勐?？h3 個(gè)產(chǎn)地組質(zhì)心距離較為相近， 勐臘縣和臨翔區(qū)能夠與其他產(chǎn)區(qū)完全分開(kāi)。 利用11 個(gè)指標(biāo)所建立的普洱生茶判別模型對(duì)5 個(gè)產(chǎn)區(qū)的茶樣進(jìn)行歸類， 并結(jié)合回判驗(yàn)證和留一交叉驗(yàn)證， 結(jié)果見(jiàn)表8。 從初始分類結(jié)果來(lái)看， 所建模型對(duì)90.0%的原始已分組生茶樣品進(jìn)行了正確分類， 并且對(duì)90.0%的交叉驗(yàn)證已分組個(gè)案進(jìn)行了正確分類。

圖5 分類散點(diǎn)圖

表8 判別分析檢驗(yàn)

3. 神經(jīng)網(wǎng)絡(luò)-多層感知器分析。 根據(jù)指標(biāo)p＜0.05、 VIP 值＞1 的標(biāo)準(zhǔn)， 將篩選出的11 種差異物質(zhì)進(jìn)行神經(jīng)網(wǎng)絡(luò)-多層感知器溯源建模。 將11 項(xiàng)指標(biāo)設(shè)置為模型的輸入層， 輸出層由5 個(gè)產(chǎn)地構(gòu)成， 以70%隨機(jī)樣品作為訓(xùn)練樣本， 30%隨機(jī)樣本作為檢驗(yàn)樣本， 結(jié)果見(jiàn)表9。 結(jié)果顯示， 訓(xùn)練集與測(cè)試集正確率總和為200.0%， 訓(xùn)練集的分類正確率為100.0%， 測(cè)試集分類正確率為100.0%。

表9 多層感知器分類正確率

三、 結(jié)論

本研究通過(guò)單因素方差分析5 個(gè)產(chǎn)地普洱生茶樣品中N、 C 穩(wěn)定同位素比值、 32 種礦質(zhì)元素、 3種有機(jī)成分含量共37 項(xiàng)指標(biāo)， 發(fā)現(xiàn)不同產(chǎn)地茶樣間的穩(wěn)定同位素、 礦質(zhì)元素、 有機(jī)成分存在一定的差異， 篩選出茶樣中具有產(chǎn)地差異且VIP 值大于1的11 種差異物質(zhì) （Cu、 Fe、 Mn、 Ca、 Mg、 Zn、Al、 S、 Kr、 Ce、 δ13C）， 構(gòu)成普洱生茶產(chǎn)地溯源信息指標(biāo)。 結(jié)合判別分析、 神經(jīng)網(wǎng)絡(luò)-多層感知器分析建立了雙江自治縣、 臨翔區(qū)、 瀾滄自治縣、 勐海縣、 勐臘縣5 個(gè)產(chǎn)區(qū)的普洱生茶溯源模型， 其中基于神經(jīng)網(wǎng)絡(luò)-多層感知器的溯源模型溯源結(jié)果最優(yōu)， 訓(xùn)練集和測(cè)試集分類正確率均達(dá)到100%； 而判別分析構(gòu)建的模型判別正確率及交叉驗(yàn)證率均為90%。 本研究為普洱生茶在小尺度空間溯源提供了一定的理論依據(jù)。

同時(shí)， 本研究在土壤-茶葉的穩(wěn)定同位素相關(guān)性分析中， 證實(shí)了不同產(chǎn)地的茶樣和土壤穩(wěn)定同位素比值存在一定的相關(guān)性。 本研究收集的樣品量不足以代表當(dāng)?shù)厮衅斩琛?由于不同產(chǎn)地普洱茶生茶的穩(wěn)定同位素、 礦質(zhì)元素、 有機(jī)元素可能受到水源、 空氣等因素影響， 還需進(jìn)一步控制變量并增加實(shí)驗(yàn)樣本， 這樣實(shí)驗(yàn)的結(jié)果將更加準(zhǔn)確和完善。