王雪蓉，羅瑞明，李亞蕾，馬夢斌

（寧夏大學(xué)食品與葡萄酒學(xué)院，寧夏 銀川 750021）

鹽池灘羊肉質(zhì)細(xì)嫩鮮美，不腥不膻，獨(dú)具特色。近年來，灘羊被其他地方引進(jìn)，但未能保持其原有品質(zhì)特性，冒用、濫用品牌現(xiàn)象阻礙灘羊產(chǎn)業(yè)的穩(wěn)定健康發(fā)展及地方特色產(chǎn)業(yè)保護(hù)[1-2]。因此，對(duì)標(biāo)識(shí)鹽池灘羊的副產(chǎn)品羊骨進(jìn)行原產(chǎn)地溯源鑒定，建立一種快速高效的肉骨鑒別技術(shù)，完善灘羊品種保護(hù)體系，是當(dāng)前亟需解決的問題。

原產(chǎn)地是影響骨產(chǎn)品商業(yè)價(jià)值的重要因素之一，比較常見的溯源方法有信息編碼技術(shù)[3]、氣相色譜法[4]、近紅外光譜法[5]、電感耦合等離子體質(zhì)譜（inductively coupled plasma-mass spectrometry，ICP-MS）法[6-8]、能量散射X射線熒光光譜方法[9]等。骨骼作為體內(nèi)礦質(zhì)元素主要的儲(chǔ)存部位，可以通過測定骨礦質(zhì)元素含量進(jìn)行區(qū)域鑒別。ICP-MS對(duì)礦質(zhì)元素檢出限低、選擇性好，被廣泛使用[10]。張欣昕[11]、駱仁軍[12]等用礦質(zhì)元素含量結(jié)合穩(wěn)定同位素比值法分別對(duì)馬鈴薯、中華絨螯蟹進(jìn)行產(chǎn)地鑒別，判別正確率均高于90.00%；Qi Jing等[13]利用機(jī)器學(xué)習(xí)輔助礦質(zhì)元素分析區(qū)分中國不同產(chǎn)地的豬肉，總體準(zhǔn)確率達(dá)95.71%；Zhang Hongru等[14]研究證明多元素分析結(jié)合線性判別模型是保證牦牛骨樣品信息可靠性的有效方法。此外，礦質(zhì)元素也被廣泛應(yīng)用于中國海參[15]、綠茶[16]、番茄[17]、葡萄酒[18]、牛羊肌肉組織[19-20]等產(chǎn)品的原產(chǎn)地溯源方面。這些研究充分說明基于礦質(zhì)元素的原產(chǎn)地鑒別技術(shù)對(duì)同一物種不同產(chǎn)區(qū)的鑒別相對(duì)較快、準(zhǔn)確性較高，具有很大的原產(chǎn)地區(qū)分潛力[21-22]。灘羊骨骼作為灘羊資源的一種副產(chǎn)品，國內(nèi)外研究相對(duì)較少，可通過建立灘羊骨骼礦質(zhì)元素溯源指紋圖譜對(duì)其進(jìn)行區(qū)域鑒別與保護(hù)。

本研究在已建立的灘羊肉原產(chǎn)地鑒別模型基礎(chǔ)上，利用ICP-MS測定灘羊骨骼中25 種（Ag、Ba、Be、Ca、Cd、Cu、Co、Cr、Cs、K、Fe、Mg、Mn、Mo、Na、Ni、P、Rb、Sb、Se、Sn、Sr、Te、V、Zn）礦質(zhì)元素含量，通過方差分析、主成分分析（principal components analysis，PCA）、聚類分析等多元統(tǒng)計(jì)方法，建立灘羊骨骼礦質(zhì)元素溯源指紋圖譜，采用Fisher判別驗(yàn)證模型有效性，為灘羊骨骼的原產(chǎn)地鑒別和檢測冒充、摻假灘羊骨產(chǎn)品提供理論依據(jù)。

1 材料與方法

1.1 材料與試劑

采集寧夏鹽池縣（麻黃山鄉(xiāng)、高沙窩鎮(zhèn)、馮記溝鄉(xiāng)）、內(nèi)蒙古自治區(qū)鄂托克前旗、陜西省定邊縣和甘肅省環(huán)縣4 地7 月齡灘羊頭骨、脊骨、四肢骨、肋骨樣品，共計(jì)240 個(gè)樣本，具體信息見表1。

表1 樣品采集信息Table 1 Information of sheep bone samples

氫氟酸（MOS級(jí)） 上海旭恒商貿(mào)有限公司；高氯酸（MOS級(jí)） 南京化學(xué)試劑有限公司；Ca、Fe、Ni、Mn、Cu、Zn、Se、Sr、Mo、Rb、Sn、P、Na、K、Mg、V、Cr、Co、Cs、Sb、Ba、Be、Cd、Te、Ag（1 000.00 μg/mL）元素標(biāo)準(zhǔn)溶液 北京儀化通標(biāo)科技有限公司；Dura series超純水處理系統(tǒng)（＞18.20 MΩ·cm） 上海和泰儀器有限公司；濃硝酸（65.00%）、濃鹽酸（37.00%） 國藥集團(tuán)化學(xué)試劑有限公司。

1.2 儀器與設(shè)備

FA-1004冷凍干燥機(jī) 上海舜宇恒平科學(xué)儀器有限公司；冷凍球磨儀 德國Restsch公司；WX-8000微波消解儀 上海新儀微波化學(xué)科技有限公司；iCAPQ ICPMS儀 美國Thermo公司。

1.3 方法

1.3.1 樣品前處理

將采集的骨骼樣品剔除表面肉、肌腱和其他雜質(zhì)，分別稱取100.00 g，冷凍干燥48 h，粉碎研磨成粉末狀，過200 目篩，備用。

1.3.2 樣品消解

稱取粉末狀樣品0.15 g至聚四氟乙烯消解罐中，加入2.00 mL雙氧水、6.00 mL濃硝酸預(yù)消解20 min，然后放入微波消解儀中消解。消解程序：0～10 min，由0 ℃以18 ℃/min速率升高到180 ℃，保持5 min；在10 min內(nèi)以6 ℃/min速率上升到240 ℃，保持25 min；最后，通風(fēng)趕酸10 min，得到消解溶液。取消解液1.00 mL定容至50 mL，開始測樣，同時(shí)做空白實(shí)驗(yàn)。

1.3.3 元素含量測定

使用ICP-MS儀測定灘羊骨骼樣品中Ag、Ba、Be、Ca、Cd、Cu、Co、Cr、Cs、K、Fe、Mg、Mn、Mo、Na、Ni、P、Rb、Sb、Se、Sn、Sr、Te、V、Zn 25 種礦質(zhì)元素含量，結(jié)果以每千克樣品干質(zhì)量計(jì)。氬用于等離子體產(chǎn)生，添加內(nèi)部標(biāo)準(zhǔn)（In、Ge、Sc元素）以補(bǔ)償任何酸效應(yīng)和儀器漂移。具體參數(shù)如下：掃描速率5 000 u/s；射頻發(fā)生器功率1 550 W；分辨率0.30～3.00 u；質(zhì)量范圍1～285 u；霧化氣流量1.12 L/min；載氣流量0.98 L/min；樣品提升速率1.20 L/min；冷卻氣流量14.00 L/min；駐留時(shí)間50.00 ms；離子透鏡電壓8.25 V；每質(zhì)量數(shù)采集數(shù)據(jù)點(diǎn)50.00 ms；積分時(shí)間500.00 ms。樣品中測定元素的含量按下式計(jì)算：

式中：x為試樣中待測元素含量/（mg/kg）；ρ1為試樣溶液中被測元素質(zhì)量濃度/（mg/L）；ρ0為試樣空白液中被測元素質(zhì)量濃度/（mg/L）；V為試樣消化液定容體積/mL；f為稀釋倍數(shù)；m為試樣稱取質(zhì)量/g。

1.4 數(shù)據(jù)處理

采用Office 2019、SPSS 23.0軟件對(duì)礦質(zhì)元素含量數(shù)據(jù)進(jìn)行處理并作方差分析、PCA、K-均值聚類等多元統(tǒng)計(jì)分析。

2 結(jié)果與分析

2.1 礦質(zhì)元素含量測定結(jié)果的準(zhǔn)確性分析

對(duì)所有樣本中25 種礦質(zhì)元素含量進(jìn)行檢測，各元素回歸方程、相關(guān)系數(shù)、相對(duì)標(biāo)準(zhǔn)偏差（relative standard deviation，RSD）和檢出限（limit of detection，LOD）見表2。

表2 各元素線性回歸方程、相關(guān)系數(shù)、RSD和LODTable 2 Linear regression equations,correlation coefficients,RSDs,and LODs for various metal elements

選擇擬合性較好的質(zhì)量濃度范圍繪制標(biāo)準(zhǔn)曲線，如表2所示，P、Se、Sn、Te 4 種礦質(zhì)元素線性方程相關(guān)系數(shù)大于0.990 0，其余20 種礦質(zhì)元素線性方程相關(guān)系數(shù)均大于0.999 0，RSD均小于7.00%。測定時(shí)用空白樣品對(duì)儀器進(jìn)樣8 次，計(jì)算得出各元素LOD在0.001 3～3.306 0 μg/L范圍內(nèi)，說明該方法可用于檢測灘羊骨骼礦質(zhì)元素含量。

2.2 不同產(chǎn)地灘羊骨骼礦質(zhì)元素分析

如表3所示，灘羊骨骼樣品中Ba、Ca、Cd、Cu、Cr、Cs、Fe、Mg、Mn、Mo、Na、Ni、P、Rb、Sb、Sn、Se、Sr、Te、V、Zn 21 種礦質(zhì)元素在4 個(gè)地區(qū)間存在顯著差異（P＜0.05），Be、Co、K元素在4 個(gè)地區(qū)間差異不顯著（P＞0.05）。

表3 不同地區(qū)間灘羊骨骼中礦物元素含量Table 3 Contents of mineral elements in Tan sheep bones from four county-level regions

內(nèi)蒙古鄂托克前旗樣品中Ba、Cr、Na、Ni、Sr、Te、V含量在4 個(gè)地區(qū)間最高，Cd、Mn、P含量最低；甘肅環(huán)縣樣品中Cd、Mg含量在4 個(gè)地區(qū)間最高，F(xiàn)e、Se 2 種元素含量最低；寧夏鹽池縣樣品中Ca、Cu、Fe、K、Mn、P、Rb和Se 8 種礦質(zhì)元素含量相對(duì)其他3 個(gè)地區(qū)最高，Ba、Ni、Te 3 種礦質(zhì)元素含量顯著低于其他3 個(gè)地區(qū)，高含量礦質(zhì)元素相對(duì)較多；陜西定邊縣樣品中Cs、Ni、Zn含量在4 個(gè)地區(qū)間最高，Sr、V 2 種元素含量最低。由此可見，不同產(chǎn)地灘羊骨骼元素組成分布特征各異，且相同礦質(zhì)元素在不同產(chǎn)地灘羊骨骼中含量也差異顯著。在此基礎(chǔ)上，對(duì)寧夏鹽池灘羊骨骼中礦質(zhì)元素進(jìn)一步分析，篩選顯著性差異礦質(zhì)元素，對(duì)其建立特有骨骼礦質(zhì)元素產(chǎn)地溯源指紋圖譜。

2.3 寧夏鹽池縣灘羊骨骼中礦質(zhì)元素含量差異分析

對(duì)寧夏鹽池縣馮記溝鄉(xiāng)、高沙窩鎮(zhèn)、麻黃山鄉(xiāng)灘羊骨骼中25 種礦質(zhì)元素含量作方差分析，如表4所示，P、Na、V、Mg、Ca、Cr、Be、Mn、Zn、Se、Sr、Cd、Cs、Ba、Rb、Te 16 種礦質(zhì)元素含量在寧夏鹽池縣3 個(gè)地區(qū)間存在顯著差異（P＜0.05）。

表4 寧夏鹽池縣馮記溝鄉(xiāng)、高沙窩鎮(zhèn)、麻黃山鄉(xiāng)灘羊骨骼中礦物元素含量Table 4 Contents of differential mineral elements among Tan sheep bones from Fengjigou,Gaoshawo and Mahuangshan in Yanchi,Ningxia

生物體中礦質(zhì)元素含量特征與地區(qū)密切相關(guān)[23]。馮記溝鄉(xiāng)灘羊骨骼中Ca、Cs、Mn、Rb、Se含量最高，其中僅Se含量與高沙窩鎮(zhèn)差異不顯著，Ba和Sr 2 種礦質(zhì)元素含量顯著低于其他兩地；高沙窩鎮(zhèn)灘羊骨骼中Ba、Be、Cd、P、Zn 5 種礦質(zhì)元素含量高于麻黃山鄉(xiāng)和馮記溝鄉(xiāng)，其中Ba、Be和Zn含量顯著高于麻黃山鄉(xiāng)和馮記溝鄉(xiāng)，而Rb、V含量顯著低于其他兩地；麻黃山鄉(xiāng)灘羊骨骼中Cr、Na、Mg、Sr、Te、V含量高于其他兩地，其中Cr、Mg和Sr含量顯著高于馮記溝鄉(xiāng)和高沙窩鎮(zhèn)，Cs、Cd、Se含量顯著低于馮記溝鄉(xiāng)和高沙窩鎮(zhèn)?？梢姡瑸┭蚬趋乐械V質(zhì)元素分布在小范圍相似地區(qū)間也存在差異，若將此特征作為產(chǎn)地鑒別標(biāo)識(shí)還存在不足，需將鹽池縣灘羊飼養(yǎng)代表性地區(qū)骨骼中具有顯著差異的礦質(zhì)元素進(jìn)一步作多元統(tǒng)計(jì)分析。

2.4 寧夏鹽池縣灘羊骨骼中礦質(zhì)元素的PCA

PCA通過降維思想研究實(shí)驗(yàn)各項(xiàng)指標(biāo)的內(nèi)在關(guān)系，在保留80%以上原有信息基礎(chǔ)上，將多項(xiàng)指標(biāo)減少為幾個(gè)相互獨(dú)立的綜合指標(biāo)[24]。相比其他統(tǒng)計(jì)學(xué)分析方法，PCA確定的權(quán)數(shù)是基于大量數(shù)據(jù)分析得到的指標(biāo)間內(nèi)在結(jié)構(gòu)關(guān)系，不受主觀因素影響，降維后得到的PC間彼此獨(dú)立，避免信息交叉，使分析結(jié)果更加客觀、可靠[25-26]。

對(duì)寧夏鹽池縣馮記溝鄉(xiāng)、高沙窩鎮(zhèn)、麻黃山鄉(xiāng)灘羊骨骼中16 種含量差異顯著的礦質(zhì)元素進(jìn)一步作PCA。由圖1可知，馮記溝鄉(xiāng)、高沙窩鎮(zhèn)、麻黃山鄉(xiāng)3 地灘羊骨骼礦質(zhì)元素中具有顯著差異的元素前4 個(gè)PC累計(jì)方差貢獻(xiàn)率分別達(dá)到85.90%、81.20%、85.50%，基本保留了原元素的變化信息。結(jié)合圖2可以非常直觀地看出，馮記溝鄉(xiāng)PC1主要由P、Na、Mg、Ca、Mn、Zn等構(gòu)成，PC2由Rb、Cd、Se、Ba 4 種元素構(gòu)成；PC3由Cs、Te構(gòu)成；高沙窩鎮(zhèn)PC1主要由P、Na、Ca、Mn、Zn、Te等元素構(gòu)成，PC2主要由Cs、V、Cd 3 種元素構(gòu)成，PC3由Rb構(gòu)成；麻黃山鄉(xiāng)PC1主要由P、Na、Mg、Ca、Zn、Cd等構(gòu)成，PC2由Cs、Ba、Rb 3 種元素構(gòu)成。PCA可以直觀地表現(xiàn)灘羊骨骼中16 種礦質(zhì)元素分布特征，找到不同地區(qū)間差異顯著的特征元素，為產(chǎn)地指紋圖譜構(gòu)建提供理論依據(jù)。

圖1 馮記溝鄉(xiāng)（a）、高沙窩鎮(zhèn)（b）、麻黃山鄉(xiāng)（c）灘羊骨骼中礦質(zhì)元素PC碎石圖Fig.1 PCA scree plot of mineral elements in Tan sheep bones from Fengjigou (a),Gaoshawo (b) and Mahuangshan (c)

圖2 馮記溝鄉(xiāng)（a）、高沙窩鎮(zhèn)（b）、麻黃山鄉(xiāng)（c）灘羊骨骼中礦質(zhì)元素PC載荷網(wǎng)絡(luò)圖Fig.2 PCA loading network plot of mineral elements in Tan sheep bones from Fengjigou (a),Gaoshawo (b) and Mahuangshan (c)

2.5 寧夏鹽池縣灘羊骨骼中礦物元素含量聚類分析

聚類即根據(jù)樣品間相似程度進(jìn)行簡化合并分組，相似性最大的元素優(yōu)先聚合在一起，最終按照類別的綜合性質(zhì)將多個(gè)樣品聚合，使結(jié)果具有最大的組內(nèi)相似性和最小的組間相似性[27-28]。采用K-均值聚類對(duì)鹽池灘羊骨骼中礦質(zhì)元素所取得的主要成分貢獻(xiàn)值進(jìn)行聚類分析，由表5可知，馮記溝鄉(xiāng)灘羊骨骼中16 種礦質(zhì)元素指標(biāo)被聚為3 類：第1類包括Ca、Mn、Na、P、Se、Zn；第2類包括Be、Cs、Mg、Rb、Te，第3類包括Ba、Cd、Cr、Sr、V，結(jié)合2.4節(jié)PCA結(jié)果，Ca、Mg、Mn、Na、P、Rb、Se、Zn這8 種元素可以作為馮記溝鄉(xiāng)灘羊骨骼的礦質(zhì)元素溯源指標(biāo)。高沙窩鎮(zhèn)灘羊骨骼中16 種礦質(zhì)元素指標(biāo)被聚為3 類：第1類包括Be、Ca、Mn、Na、P、V、Te、Zn；第2類包括Ba、Cd、Cs、Mg、Se，第3類包括Cr、Rb、Sr，結(jié)合2.4節(jié)PCA結(jié)果，Ca、Cd、Cs、Mn、Na、P、Te、Zn這8 種元素可以作為高沙窩鎮(zhèn)灘羊骨骼的礦質(zhì)元素溯源指標(biāo)。麻黃山鄉(xiāng)灘羊骨骼中16 種礦質(zhì)元素指標(biāo)被聚為3 類：第1類包括Ca、Cd、Mg、Mn、Na、P、Se、Zn；第2類包括Be、Cr、Cs、Rb、Sr、Te，第3類包括Ba、V，結(jié)合2.4節(jié)PCA結(jié)果，Ca、Cd、Cs、Mg、Na、P、Rb、Zn這8 種元素可作為麻黃山鄉(xiāng)灘羊骨骼的礦質(zhì)元素溯源指標(biāo)。綜上所述，Ca、Cd、Cs、Mg、Mn、Na、P、Rb、Te、Zn這10 種礦質(zhì)元素可以作為鹽池縣灘羊骨骼中礦質(zhì)元素的溯源指標(biāo)。

表5 馮記溝鄉(xiāng)、高沙窩鎮(zhèn)、麻黃山鄉(xiāng)灘羊骨骼中礦質(zhì)元素含量聚類分析結(jié)果Table 5 Cluster analysis of mineral elements in Tan sheep bones from three towns in Yanchi county

2.6 寧夏鹽池灘羊骨骼中礦物質(zhì)元素的判別分析

選取采集的陜西定邊縣、內(nèi)蒙古鄂托克前旗、寧夏鹽池縣和甘肅環(huán)縣4 地灘羊骨骼樣品，測定其礦質(zhì)元素含量，檢驗(yàn)已篩選出的溯源指標(biāo)對(duì)隨機(jī)灘羊骨骼樣品產(chǎn)地鑒別效果。利用2.5節(jié)中鹽池灘羊骨骼溯源指標(biāo)中的10 種元素（Ca、Cd、Cs、Mg、Mn、Na、P、Rb、Te、Zn）建立判別模型。在0.05顯著性水平上，將Ca、Cs、Mg、Mn、Se、P、Te 7 種對(duì)地域判別顯著的元素引入到Fisher判別模型中。具體判別模型如下：

y1＝71.07x1＋49.88x2－0.24x3－412.01x4＋11.05x5＋163.16x6－8.63x7－1 217.90

y2＝194.25x1－18.70x2＋0.14x3＋208.19x4－6.53x5＋107.54x6＋5.27x7－423.76

y3＝64.44x1－23.74x2＋0.10x3＋216.45x4＋5.59x5＋66.30x6＋4.17x7＋263.07

y4＝53.98x1＋14.52x2＋0.06x3＋321.12x4＋3.99x5－27.75x6＋1.98x7－122.76

式中：y1、y2、y3、y4分別為寧夏鹽池縣、內(nèi)蒙古鄂托克前旗、甘肅環(huán)縣、陜西定邊縣4 個(gè)地區(qū)模型計(jì)算值；x1、x2、x3、x4、x5、x6、x7分別為Ca、Cs、Mg、Mn、Se、P、Te含量/（mg/kg）。

利用此判別模型對(duì)樣品進(jìn)行歸類，將礦質(zhì)元素含量測定值代入上述模型中計(jì)算，比較所得y值，將樣品歸為y值最大的一類，并結(jié)合回代檢驗(yàn)和交叉檢驗(yàn)對(duì)模型有效性進(jìn)行驗(yàn)證。由表6可知，灘羊骨骼樣品來源整體判別效果較好。4 個(gè)地區(qū)的回代檢驗(yàn)和交叉檢驗(yàn)的整體正確率均在85.00%以上，其中對(duì)寧夏鹽池縣的判別效果最好，回代檢驗(yàn)和交叉檢驗(yàn)正確率均達(dá)90.00%。陜西定邊縣骨骼樣品被錯(cuò)判為內(nèi)蒙古鄂托克前旗和甘肅環(huán)縣，致使檢驗(yàn)正確率稍低，同時(shí)，其他地域被錯(cuò)判的樣品多被錯(cuò)判為陜西定邊縣樣品，可能是其中差異不顯著元素占比較大，區(qū)分性不高。此外，寧夏鹽池縣和甘肅環(huán)縣樣品也存在相互錯(cuò)判。

表6 4 個(gè)地區(qū)灘羊骨骼樣品的LDA結(jié)果Table 6 LDA results of Tan sheep bone samples from four county-level regions

對(duì)樣品前2 個(gè)判別函數(shù)得分作散點(diǎn)圖，如圖3所示，每個(gè)地區(qū)樣品均被較好的歸類，各具空間分布特征，這與PCA、聚類分析結(jié)果基本一致，Ca、Cd、Cs、Mg、Mn、Na、P、Rb、Te、Zn可作為鹽池灘羊骨骼礦質(zhì)元素指紋圖譜，表明篩選出用作灘羊骨骼產(chǎn)地鑒別的10 種礦質(zhì)元素是有效的，并且利用該10 種元素指標(biāo)建立的判別模型亦有效。

圖3 灘羊骨骼樣品的前兩個(gè)判別函數(shù)得分散點(diǎn)圖Fig.3 Scatter plot of Tan sheep bone samples on the first two discriminant functions

本研究中篩選的礦質(zhì)元素溯源指紋圖譜，整體正確判別率約為85.00%，與Sun Shumin[19]、Zhang Hongru[14]等的研究結(jié)果相比，判別正確率較低，可能是土壤、氣候、飼喂草料等諸多因素對(duì)礦質(zhì)元素含量有影響，進(jìn)而影響判別結(jié)果[29-30]。后期擬添加氣候、土壤等環(huán)境因素，采集大樣本量，結(jié)合具有地域代表性的穩(wěn)定同位素、化學(xué)計(jì)量方法進(jìn)行研究，繼續(xù)完善指紋圖譜信息，提高判別正確率，從而有效區(qū)分不同產(chǎn)地灘羊骨骼，為羊骨的產(chǎn)地鑒別提供理論依據(jù)。

3 結(jié)論

通過對(duì)寧夏鹽池縣、內(nèi)蒙古鄂托克前旗、甘肅環(huán)縣和陜西定邊縣4 個(gè)地區(qū)灘羊骨骼中25 種（Ag、Ba、Be、Ca、Cd、Cu、Co、Cr、Cs、K、Fe、Mg、Mn、Mo、Na、Ni、P、Rb、Sb、Se、Sn、Sr、Te、V、Zn）礦質(zhì)元素含量進(jìn)行測定，運(yùn)用PCA等多元統(tǒng)計(jì)分析方法得出以下結(jié)論：1）灘羊骨骼樣品中Ba、Ca、Cd、Cu、Cr、Cs、Fe、Mg、Mn、Mo、Na、Ni、P、Rb、Se、Sr、Te、V、Zn 19 種礦質(zhì)元素在4 個(gè)地區(qū)間存在顯著差異，內(nèi)蒙古鄂托克前旗樣品中Ba、Cr、Na、Ni、Sr、Te、V含量在4 個(gè)地區(qū)間最高，甘肅環(huán)縣樣品中Cd、Mg含量最高，寧夏鹽池縣樣品中Ca、Cu、Fe、K、Mn、P、Rb、Se含量最高，陜西定邊縣樣品中Cs、Ni、Zn含量最高，不同產(chǎn)地間灘羊骨骼的組成元素各具分布特征，同時(shí)相同礦質(zhì)元素在不同原產(chǎn)地間也差異顯著。2）篩選出鹽池灘羊骨骼礦質(zhì)元素指紋溯源圖譜為Ca、Cd、Cs、Mg、Mn、Na、P、Rb、Te、Zn，通過Fisher判別，Ca、Cs、Mg、Mn、Se、P、Te 7 種礦質(zhì)元素被引入判別模型中，回代檢驗(yàn)的整體正確判別率為87.50%，交叉檢驗(yàn)的整體正確判別率為86.87%，建立的礦質(zhì)元素溯源指紋圖譜可用于區(qū)分不同地區(qū)灘羊骨骼。