基于礦物元素的西藏自治區(qū)牦牛肉產(chǎn)地溯源研究

宗萬里劉海金趙姍姍陳愛亮趙 燕

（1.中國(guó)農(nóng)業(yè)科學(xué)院農(nóng)業(yè)質(zhì)量標(biāo)準(zhǔn)與檢測(cè)技術(shù)研究所，農(nóng)業(yè)農(nóng)村部農(nóng)產(chǎn)品質(zhì)量安全重點(diǎn)實(shí)驗(yàn)室，北京 100081；2.威海市食品藥品檢驗(yàn)檢測(cè)研究院，山東威海 264210；3.西藏自治區(qū)農(nóng)畜產(chǎn)品質(zhì)量安全檢驗(yàn)檢測(cè)中心，拉薩 850211）

西藏自治區(qū)是我國(guó)牦牛的主產(chǎn)地之一，由于其特殊的地理?xiàng)l件和生態(tài)環(huán)境，形成了九龍牦牛、斯布牦牛、帕里牦牛、類烏齊牦牛、嘉黎牦牛等特有的類群，不同地域牦牛肉的品質(zhì)和口感存在很大的差異。為防止西藏自治區(qū)存在低質(zhì)量牦牛肉假冒高品質(zhì)牦牛肉的現(xiàn)象，針對(duì)不同產(chǎn)地牦牛肉的肉品性狀和顏色極為相似以及肉眼很難區(qū)分的問題，有必要對(duì)不同地域牦牛肉產(chǎn)地鑒別技術(shù)進(jìn)行研究。目前基于礦物元素分析的產(chǎn)地溯源技術(shù)已經(jīng)應(yīng)用于畜產(chǎn)品[1～12]、 乳及乳制品[13～17]、 水產(chǎn)品[18～19]、 水果[20～22]、谷物[23～29]、 中藥材[30～31]、 蔬菜[32～33]、 茶 葉[34]等產(chǎn)地溯源研究中。選擇礦物元素進(jìn)行溯源的依據(jù)是不同地域環(huán)境中礦物元素含量存在差異，從而導(dǎo)致植物中礦物元素含量存在差異，進(jìn)而以植物為飼料的動(dòng)物體內(nèi)礦物元素含量存在差異。本研究以西藏自治區(qū)4個(gè)不同地區(qū)的100個(gè)牦牛肉樣本為實(shí)驗(yàn)材料，利用電感耦合等離子體發(fā)射光譜儀（ICP－OES）、電感耦合等離子體質(zhì)譜儀（ICP－MS）分析4個(gè)不同地區(qū)牦牛肉樣品中的礦物元素組成差異特征。通過牦牛肉中礦物元素差異，利用化學(xué)計(jì)量學(xué)方法，建立區(qū)分西藏自治區(qū)不同產(chǎn)地牦牛肉的判別模型，對(duì)不同產(chǎn)地牦牛肉進(jìn)行判別，旨在為西藏自治區(qū)不同產(chǎn)地牦牛肉的溯源提供參考依據(jù)。

一、材料與方法

（一）材料與 試劑鉀 （K）、鈉 （Na）、鈣（Ca）、鎂 （Mg）、鋁 （Al）、鐵 （Fe）、銅 （Cu）、鋅（Zn）、釩 （V）、鉻 （Cr）、錳 （Mn）、鈷 （Co）、鎳（Ni）、鎵 （Ga）、砷 （As）、硒 （Se）、銣 （Rb）、鍶（Sr）、銀（Ag）、鎘（Cd）、鋇（Ba）、鉈（Tl）、鉛（Pb）、釔 （Y）、鑭 （La）、鈰 （Ce）、鐠 （Pr）、釹 （Nd）、釤（Sm）、銪（Eu）、釓（Gd）、鋱（Tb）、鏑（Dy）、鈥（Ho）、鉺（Er）、銩 （Tm）、鐿（Yb）元素標(biāo)準(zhǔn)溶液，100 μg/mL，國(guó)家有色金屬及電子材料分析測(cè)試中心；鈦（Ti）、錫（Sn）、銻（Sb）元素標(biāo)準(zhǔn)溶液，10 μg/mL，安捷倫科技（中國(guó)）有限公司；實(shí)驗(yàn)用水均為GB/T 6682規(guī)定的一級(jí)水；硝酸（色譜純），德國(guó)Merck公司；三氯甲烷（分析純），天津市科密歐化學(xué)試劑有限公司；甲醇（分析純），天津市科密歐化學(xué)試劑有限公司；雙氧水（優(yōu)級(jí)純），國(guó)藥集團(tuán)化學(xué)試劑有限公司；玻璃培養(yǎng)皿，直徑6 cm，上海五一玻璃儀器廠；刻度消解管，50 mL，美國(guó)EE公司。

（二）儀器與設(shè)備GM 300型刀式研磨儀，德國(guó)萊馳公司；MARS 6型微波消解儀，美國(guó)培安公司；7900型電感耦合等離子體質(zhì)譜儀，美國(guó)安捷倫公司；7400型電感耦合等離子體發(fā)射光譜儀，美國(guó)賽默飛世爾公司；Practum 1102－1 CN型電子天平，德國(guó)賽多利斯公司；SCIENTZ－48型高通量組織研磨器，寧波新芝生物科技股份有限公司；Comfort II型超純水機(jī)，德國(guó)賽多利斯公司；EHD－16型電熱趕酸儀，北京東航科儀儀器有限公司；FD260型干燥箱，德國(guó)賓得公司；UMV－2型多位渦旋振蕩器，煙臺(tái)盈聚儀器儀表有限公司；Z216MK型冷凍離心機(jī)，德國(guó)Hermle公司。

（三）方法

1.樣品采集。采集的牦牛肉樣品來自西藏自治區(qū)的拉薩達(dá)孜、阿里、昌都類烏齊和日喀則亞東，樣品詳細(xì)信息見表1。

表1 采集的牦牛肉樣品的地理及其環(huán)境信息

2.樣品前處理。取經(jīng)粉碎的牦牛肉樣品15.0 g左右平攤于平皿中，然后放到干燥箱中（70℃）干燥48 h，取干燥的樣品用高通量組織研磨器（參數(shù)設(shè)置為50 Hz、300 s）磨成粉。將磨成粉的牦牛肉樣品進(jìn)行脫脂，脫脂的步驟為稱取0.3 g牦牛肉粉至2 mL離心管中，加入1.5 mL三氯甲烷－甲醇（2∶1，體積比）溶液，置于離心管并蓋緊蓋子。將離心管放入多位渦旋振蕩器，以2 500 r/min的速度振蕩10 min。振蕩結(jié)束后，在12℃的溫度下以5 000 r/min的速度離心5 min。用移液槍吸取上清液至廢液瓶中。再重復(fù)脫脂兩次，將脫脂后的樣品置于干燥器中干燥。

稱取0.2 g干燥后的樣品至微波消解管中，加入7 mL硝酸、1 mL雙氧水，采用微波消解法消解。微波消解程序?yàn)閺氖覝厣?20℃（20 min），保持20 min，從120℃升至200℃（15 min），保持30 min，功率為1 400 W。

待消解完成后將消解管取出置于電熱趕酸儀上110℃趕酸2 h，冷卻后將消解液轉(zhuǎn)移至刻度消解管中，用純水少量多次潤(rùn)洗消解管內(nèi)壁并轉(zhuǎn)移至同一刻度消解管中，再用純水定容至25 mL，作為上機(jī)測(cè)試溶液。

元素K、Na、Ca、Mg、Al、Fe、Cu、Zn使用電感耦合等離子體發(fā)射光譜法檢測(cè)，其他元素V、Cr、Mn、Co、Ni、Ga、As、Se、Rb、Sr、Ag、Cd、Ba、Tl、Pb、Ti、Sn、Sb、Y、La、Ce、Pr、Nd、Sm、Eu、Gd、Tb、Dy、Ho、Er、Tm、Yb使用電感耦合等離子體質(zhì)譜法檢測(cè)[35～36]。

3.儀器條件。電感耦合等離子體發(fā)射光譜、電感耦合等離子體質(zhì)譜的儀器條件分別見表2和表3。

表2 電感耦合等離子體發(fā)射光譜條件

表3 電感耦合等離子體質(zhì)譜條件

4.檢出限及定量限。40種礦物元素的檢出限及定量限見表4。

表4 40種礦物元素檢出限及定量限

5.統(tǒng)計(jì)學(xué)分析。應(yīng)用統(tǒng)計(jì)學(xué)軟件SPSS 22.0（IBM公司）和SIMCA 14.1（Umetrics公司）對(duì)不同產(chǎn)地牦牛肉的元素含量進(jìn)行多重方差分析、聚類分析（HCA）、主成分分析（PCA）、偏最小二乘法判別分析（PLS－DA）和基于正交信號(hào)校正的偏最小二乘判別分析（OPLS－DA）。

二、結(jié)果與分析

（一）不同地區(qū)牦牛肉中礦物元素含量差異分析對(duì)西藏自治區(qū)4個(gè)不同產(chǎn)地耗牛肉樣品中的40種礦物元素含量進(jìn)行方差分析，分析結(jié)果見表5。由表5可知，牦牛肉樣品中Mg、K、Na、Ca、Zn、Cu、Al、Fe、Mn、Cr、Ni、Ga、Se、Rb、Sr、Ba、Pb、Ti、Sn含量較高，V、Co、As、Ag、Cd、Tl、Sb、Y、La、Ce、Pr、Nd、Sm、Eu、Gd、Tb、Dy、Ho、Er、Tm、Yb含量較低。

表5 西藏自治區(qū)4個(gè)不同產(chǎn)地牦牛肉中礦物元素方差分析

牦牛肉樣品中Er含量在4個(gè)產(chǎn)地間存在顯著差異（P<0.05），Na、Cu、Fe、Ga、As、Se、Rb、Cd、Ba、Y、La、Ce、Pr、Nd、Sm、Eu、Gd、Dy、Ho、Tm共20種礦物元素的含量在3個(gè)產(chǎn)地間存在顯著差異，Mg、K、Al、Zn、V、Cr、Mn、Co、Ni、Ag、Tl、Pb、Ti、Sn、Sb、Tb、Yb共17種礦物元素的含量在2個(gè)產(chǎn)地間存在顯著差異（P<0.05），元素Ca、Sr含量差異不顯著。拉薩達(dá)孜的牦牛肉樣品中K、Na、Cu、Rb、Cd、Ba、Pb、Sn、La、Ce、Pr、Nd、Sm、Eu、Tb、Dy、Ho、Er、Tm含量高于其他3個(gè)產(chǎn)地，F(xiàn)e、V、Cr、Co、Ni、Ga、Se、Ag低于其他3個(gè)產(chǎn)地；阿里的牦牛肉樣品中Zn、Fe、Mn、As、Sb含量高于其他3個(gè)產(chǎn)地，Rb、Cd、Tl、Tb含量低于其他3個(gè)產(chǎn)地；昌都類烏齊的牦牛肉樣品中Mg、Cr、Ni、Ga、Se含量高于其他3個(gè)產(chǎn)地，Cu、Al、Mn、As、Ba、Pb、Ti、Sn、Sb、Y、La、Ce、Pr、Nd、Sm、Eu、Gd、Dy、Ho、Er、Tm、Yb含量低于其他3個(gè)產(chǎn)地；日喀則亞東的牦牛肉樣品中Al、V、Co、Ag、Tl、Ti、Y、Gd含量高于其他3個(gè)產(chǎn)地，Mg、K、Na、Zn含量低于其他3個(gè)產(chǎn)地。綜上，在測(cè)定的40種元素中存在差異的有38種，不同產(chǎn)地的牦牛肉樣品中礦物元素含量有其各自的特征，說明基于礦物元素分析對(duì)牦牛肉產(chǎn)地進(jìn)行判別是可行的。

動(dòng)物體內(nèi)礦物元素含量受水、植物的影響[37]，水、植物中礦質(zhì)元素的積累與環(huán)境中礦物元素含量呈一定相關(guān)性[38]，動(dòng)物體內(nèi)的礦物元素含量與周圍的環(huán)境有一定相關(guān)性。西藏自治區(qū)的牦牛肉中稀土元素含量低，與地理環(huán)境有關(guān)，我國(guó)土壤中稀土元素的含量由南向北和由東向西呈逐漸降低的態(tài)勢(shì)[39]。從礦物元素含量差異分析結(jié)果可以看出，西藏自治區(qū)不同產(chǎn)地牦牛肉中礦物元素含量存在差異。牦牛肉中含量較高的元素Mg、K、Na、Zn、Cu、Al、Fe、Mn、Cr、Ni、Ga、Se、Rb、Ba、Pb、Ti、Sn在不同產(chǎn)地之間存在差異，這可能是由土壤中可利用的礦物元素含量差異引起的。但是，這些元素中沒有一種元素同時(shí)在4個(gè)產(chǎn)地之間存在差異，說明單純利用方差分析不能較好地對(duì)不同產(chǎn)地牦牛肉進(jìn)行區(qū)分，需要利用其他分析方法對(duì)數(shù)據(jù)進(jìn)行分析。

（二）聚類分析聚類分析是歸類劃分?jǐn)?shù)據(jù)中相似性數(shù)據(jù)的統(tǒng)計(jì)方法[33]，是將樣本根據(jù)其特征的相似性逐步聚合，把相似度最大的優(yōu)先聚合在一起，最后根據(jù)類別的綜合性對(duì)多個(gè)類別進(jìn)行聚合，從而完成聚類分析的過程[40]。近年來，聚類分析越來越多地應(yīng)用于農(nóng)產(chǎn)品產(chǎn)地溯源[41～45]。本研究對(duì)38種存在含量顯著差異的礦物元素進(jìn)行基于歐式距離的層次聚類分析，產(chǎn)生的聚類樹狀圖見圖1。圖1橫坐標(biāo)樣品編號(hào)中LS代表來自拉薩達(dá)孜的樣品、CD代表來自昌都類烏齊的樣品、RKZ代表來自日喀則亞東的樣品，AL代表來自阿里的樣品。將聚類距離200作為分割線，4個(gè)產(chǎn)地的牦牛肉樣品可被分為4類。其中阿里的牦牛肉樣品歸屬一類，有個(gè)別拉薩達(dá)孜牦牛肉樣品歸屬到這一類中；昌都類烏齊的牦牛肉樣品歸屬一類，有日喀則亞東的個(gè)別牦牛肉樣品分類至這類中；拉薩達(dá)孜和日喀則亞東的樣品各自歸為一類，但分別有少數(shù)樣品歸到其他類?？梢姡诘V物元素含量的聚類分析為牦牛肉產(chǎn)地鑒別提供了可能，但仍存在交叉現(xiàn)象。

圖1 基于礦物元素的4個(gè)產(chǎn)地牦牛肉HCA圖

（三）主成分分析主成分分析是利用降維的思想，在損失很少信息的前提下，把多個(gè)指標(biāo)轉(zhuǎn)化為幾個(gè)綜合指標(biāo)的多元統(tǒng)計(jì)方法[46]。通過主成分分析結(jié)果可以判斷篩選的特征元素的差異情況。本研究對(duì)在牦牛肉中具有產(chǎn)地間顯著差異的38種元素進(jìn)行主成分分析，共提取了38個(gè)主成分，其中初始特征值大于1的主成分有8個(gè)，累積方差百分比為79.525%（見表6），表明8個(gè)主成分包含了原始變量中79%以上的信息，可代表大部分?jǐn)?shù)據(jù)特征。第1主成分主要綜合了牦牛肉樣品的Al、V、Ba、Ti、Y、La、Ce、Pr、Nd、Sm、Eu、Gd、Dy、Ho、Er、Tm和Yb共17種 元素 含 量信息；第2主成分主要綜合了牦牛肉樣品的Fe、V、Cr、Mn、Co、Ni、Ga、Sb共8種元素含量信息；第3主成分主要綜合了牦牛肉樣品的V、Cr、Co、Ni、Ga、Tl共6種元素含量信息；第4主成分主要綜合了牦牛肉樣品的Cu、Rb、Pb、Sn共4種元素含量信息；第5主成分主要綜合了牦牛肉樣品的Mg、Rb和Pb共3種元素含量信息；第6主成分主要綜合了牦牛肉樣品的Na、Ag和Cd共3種元素含量信息；第7主成分主要綜合了牦牛肉樣品的Mg元素含量信息；第8主成分主要綜合了牦牛肉樣品的K元素含量信息。

表6 總方差解釋

在主成分分析中，利用第1主成分、第2主成分得到得分圖（見圖2）。由圖2可以看出，4個(gè)產(chǎn)地牦牛肉樣品區(qū)分不夠明顯，日喀則亞東的牦牛肉樣品與拉薩達(dá)孜、阿里的牦牛肉樣品存在交叉的情況；拉薩達(dá)孜有3個(gè)樣品落在置信橢圓外，其中有1個(gè)樣品與同產(chǎn)地其他樣品相距較遠(yuǎn)；日喀則亞東有兩個(gè)樣品落在置信橢圓外，與同產(chǎn)地其他樣品相距較遠(yuǎn)；阿里有3個(gè)樣品落在置信橢圓外，與同地區(qū)其他樣品相距較遠(yuǎn)。從地理位置上看，日喀則亞東與阿里和拉薩達(dá)孜的距離較近，與昌都類烏齊的距離遠(yuǎn)，在圖2中可以看到日喀則亞東的樣品與昌都類烏齊的樣品區(qū)分明顯，而與拉薩達(dá)孜和阿里的樣品區(qū)分效果較差，其分散在拉薩達(dá)孜和阿里的樣品之間。

圖2 基于礦物元素的4個(gè)產(chǎn)地牦牛肉PCA得分圖

（四）判別分析

1.Fisher判別模型?；陉笈Ｈ猱a(chǎn)地進(jìn)行元素差異分析、聚類分析和主成分分析后可以看出，利用礦物元素的含量對(duì)牦牛肉產(chǎn)地溯源具有一定可行性。為進(jìn)一步了解各礦物元素含量指標(biāo)對(duì)牦牛肉產(chǎn)地判別的作用及效果，對(duì)在產(chǎn)地間存在差異的38種礦物元素，用Fisher函數(shù)建立判別模型，進(jìn)行判別分析。判別分析是一種有監(jiān)督的數(shù)據(jù)分類方式。其可利用各類樣品之間的相關(guān)性建立判別模型而對(duì)未知樣品進(jìn)行分類，還可利用交叉驗(yàn)證把除某一個(gè)樣品之外的其余樣品作為訓(xùn)練集，進(jìn)而建立判別函數(shù)對(duì)該樣品進(jìn)行判別，從而檢驗(yàn)判別模型的穩(wěn)定性。

本研究在Fisher判別分析中，對(duì)樣品進(jìn)行了原始判別分析和交叉驗(yàn)證分析。結(jié)果顯示，來自拉薩達(dá)孜、阿里、昌都類烏齊、日喀則亞東的牦牛肉的原始正確判別率均為100%，總體原始正確判別率為100%；來自拉薩達(dá)孜的牦牛肉的交叉驗(yàn)證正確判別率為100%，來自阿里的牦牛肉的交叉驗(yàn)證正確判別率為96.7%，來自昌都類烏齊的牦牛肉的交叉驗(yàn)證正確判別率為97.7%，來自日喀則亞東的牦牛肉的交叉驗(yàn)證正確判別率為94.1%，總體交叉驗(yàn)證正確判別率為97%。可見西藏自治區(qū)不同產(chǎn)地的牦牛肉，借助于礦物元素能夠進(jìn)行有效的產(chǎn)地判別。

2.偏最小二乘法判別分析。偏最小二乘法判別分析是一種用于判別分析的多變量統(tǒng)計(jì)分析方法，是特征投影顯示方法的一種，可揭示自變量（X）和因變量（Y）之間的模型關(guān)系[47]。為綜合分析牦牛肉不同產(chǎn)地間的差異性，去掉0值較多的Tb元素，將剩余的39個(gè)元素含量輸入SIMCA 14.1分析軟件，建立有監(jiān)督模式的PLS－DA模型。提取2個(gè)主成分，模型解釋率參數(shù)R2Y為0.559，模型預(yù)測(cè)率參數(shù)Q2為0.537。經(jīng)PLS－DA模型分析分別生成得分圖、載荷圖和變量重要性投影值圖（VIP圖），分別見圖3～圖5。

圖3 基于礦物元素的4個(gè)產(chǎn)地牦牛肉PLS－DA得分圖

圖3中的每個(gè)點(diǎn)代表1個(gè)樣品，聚合程度反映同一產(chǎn)地牦牛肉樣品間相似度。由圖3可知，PLS－DA模型中數(shù)據(jù)點(diǎn)均在95%置信區(qū)間內(nèi)，聚合程度較好，但日喀則亞東牦牛肉樣品與拉薩達(dá)孜以及昌都類烏齊牦牛肉樣品存在交叉現(xiàn)象。

圖4中每1個(gè)點(diǎn)代表1種元素，距離原點(diǎn)越遠(yuǎn)的點(diǎn)權(quán)重值越大，決定樣本差異的作用越大。圖5中的VIP值可反映變量（元素）對(duì)模型分類的整體貢獻(xiàn)度，PLS－DA模型中變量VIP值＞1.0說明該變量對(duì)整體模型的貢獻(xiàn)度高于平均水平。由圖5可知，As、Tl、Rb、Ag、Se、Sb、Mn、Cd、Cr、Fe、Na、Ni、Y、Ga、Co貢獻(xiàn)度較大，這些元素可能是不同產(chǎn)地牦牛肉的主要標(biāo)志性元素。

圖4 基于礦物元素的4個(gè)產(chǎn)地牦牛肉PLS-DA載荷圖

圖5 PLS-DA模型的元素VIP值分布情況

3.基于正交信號(hào)校正的偏最小二乘判別分析。按和PLS－DA同樣的方法建立有監(jiān)督模式的OPLS－DA模型。提取2個(gè)主成分，模型解釋率參數(shù)R2Y為0.709，高于0.5，說明所建模型具有較強(qiáng)的解釋率；模型預(yù)測(cè)率參數(shù)Q2為0.596，大于0.5。經(jīng)OPLS－DA模型分析分別生成得分圖、載荷圖和變量重要性投影值圖（VIP圖），分別見圖6～圖8。

圖6 基于礦物元素的4個(gè)產(chǎn)地牦牛肉OPLS-DA得分圖

圖8 OPLS-DA模型的元素VIP值分布情況

從圖6中可以看出，基于牦牛肉中礦物元素含量構(gòu)建的OPLS－DA模型可以明顯區(qū)分拉薩達(dá)孜、阿里、昌都類烏齊、日喀則亞東4個(gè)產(chǎn)地的牦牛肉樣品。OPLS－DA模型中數(shù)據(jù)點(diǎn)均在95%置信區(qū)間內(nèi)，且聚合程度與區(qū)分性均優(yōu)于PLS－DA模型，說明在區(qū)分4個(gè)不同產(chǎn)地牦牛肉樣品的元素含量差異性上，OPLS－DA構(gòu)建的判別模型效果更好。

與PLS－DA類似，圖7中每個(gè)點(diǎn)代表1種元素，距離原點(diǎn)越遠(yuǎn)的點(diǎn)權(quán)重值越大，決定樣本差異的作用越大。圖8中變量（元素）VIP值＞1.0說明該變量（元素）對(duì)OPLS－DA模型的貢獻(xiàn)度高于平均水平。由圖8可知，Tl、As、Ag、Se、Cr、Rb、Ni、Mn、Sb、Co、Ce、Al、Y、Pr、Fe、La、Ga、Cd、V、Na、Er貢 獻(xiàn) 度 較 大，這 些 元 素 具 有 在OPLS－DA模型中作為判別產(chǎn)地差異元素的潛力。

圖7 基于礦物元素的4個(gè)產(chǎn)地牦牛肉OPLS-DA載荷圖

在判別分析中，F(xiàn)isher判別分析、PLS－DA、OPLS－DA均能較好實(shí)現(xiàn)對(duì)4個(gè)產(chǎn)地牦牛肉樣品的區(qū)分，其中Fisher判別分析和OPLS－DA比PLSDA的區(qū)分效果好。在牦牛的產(chǎn)地溯源中，項(xiàng)洋等[48]利用穩(wěn)定同位素與礦物元素相結(jié)合的方法實(shí)現(xiàn)了對(duì)青海省大通縣、四川省紅原縣、四川省九龍縣、青海省河南縣牦牛肉的區(qū)分，選用了δ13C、δ2H、Ti，總體正確分類率為85%，交叉驗(yàn)證率為82.5%。本研究進(jìn)一步證實(shí)了利用礦物元素分析技術(shù)結(jié)合化學(xué)計(jì)量學(xué)的方法能夠?qū)崿F(xiàn)對(duì)同一省內(nèi)不同產(chǎn)地牦牛肉的產(chǎn)地溯源。

三、結(jié)論

本研究采用基于電感耦合等離子質(zhì)譜法和電感耦合等離子體發(fā)射光譜法的礦物元素分析技術(shù)，對(duì)拉薩達(dá)孜、阿里、昌都類烏齊、日喀則亞東4個(gè)產(chǎn)地100個(gè)牦牛肉樣品中的40種礦物元素進(jìn)行了測(cè)定，并進(jìn)行多重方差分析、聚類分析、主成分分析和判別分析等化學(xué)計(jì)量學(xué)分析，實(shí)現(xiàn)了對(duì)西藏自治區(qū)不同產(chǎn)地牦牛肉的區(qū)分，并成功建立了判別模型。由分析結(jié)果可知，多重方差分析、聚類分析、主成分分析能夠?qū)?個(gè)產(chǎn)地的牦牛肉進(jìn)行一定程度的區(qū)分。Fisher判別分析可以有效判別牦牛肉產(chǎn)地，其總體原始正確判別率為100%，交叉驗(yàn)證正確判別率為97%。建立的PLS－DA數(shù)據(jù)模型和OPLS－DA數(shù)據(jù)模型可對(duì)牦牛肉產(chǎn)地實(shí)現(xiàn)有效的區(qū)分，其中OPLS－DA模型的區(qū)分效果優(yōu)于PLSDA模型。本研究結(jié)果證明了利用礦物元素分析技術(shù)建立的產(chǎn)地溯源模型可對(duì)來自我國(guó)西藏自治區(qū)拉薩達(dá)孜、阿里、昌都類烏齊、日喀則亞東4個(gè)產(chǎn)地的牦牛肉進(jìn)行有效的產(chǎn)地判別，礦物元素分析技術(shù)在西藏自治區(qū)牦牛肉產(chǎn)地溯源應(yīng)用方面具有可行性。