基于PCA-MLR的礦區(qū)耕地土壤重金屬污染來源分析

摘要：以鉛鋅礦區(qū)的X縣1 078個耕地樣本點為研究對象，通過多元統(tǒng)計分析法得出污染源溯源的定性結(jié)果，利用PCA-MLR模型計算各污染源的貢獻率。結(jié)果表明，X縣土壤重金屬污染的總體情況良好，但存在輕度污染、中度污染的風(fēng)險；相關(guān)性分析和聚類分析表明，Cr、As、Pb三者存在明顯相關(guān)性，Cu、Zn二者也存在明顯相關(guān)性，可能存在相似的污染源；Pb、Cr、As的污染源主要為礦業(yè)排放，Cu、Zn、Cd的污染源主要為農(nóng)業(yè)面源污染，Ni的污染源主要為自然來源。在對污染源定性分析的基礎(chǔ)上，通過PCA-MLR模型定量確認(rèn)不同污染源的貢獻率，X縣及A鎮(zhèn)的主要污染源均為農(nóng)業(yè)面源污染。相較于傳統(tǒng)方法，該方法在污染源解析方面具有較高的準(zhǔn)確性和實用性。

關(guān)鍵詞：PCA-MLR模型；礦區(qū)耕地；土壤；重金屬；污染來源

中圖分類號：S151.9" " " " "文獻標(biāo)識碼：A

文章編號：0439-8114（2024）11-0185-06

DOI：10.14088/j.cnki.issn0439-8114.2024.11.031 開放科學(xué)（資源服務(wù)）標(biāo)識碼（OSID）：

Source analysis of heavy metal pollution in farmland soil in mining areas

based on PCA-MLR

SONG Ke1， WANG Hao1， YANG Wei2， XU Jia-jun1， CHU Zhao-yuan1， DU Ming-ming1， CAO Zhi-yong1

（1. College of Big Data， Yunnan Agricultural University， Kunming" 650500，China； 2. Kunming Agricultural Product Quality and Safety Center， Kunming" 650118，China）

Abstract： Taking the 1 078 farmland sampling points in X County， a lead-zinc mining area， as the research object， the qualitative results of pollution source tracing were obtained through multivariate statistical analysis， and the contribution rate of each pollution source was calculated using the PCA-MLR model. The results indicated that the overall situation of heavy metal pollution in the soil of X County was good， but there was a risk of mild to moderate pollution；correlation analysis and cluster analysis indicated that there was a significant correlation among Cr， As， and Pb， as well as among Cu， Zn， and Cd， suggesting the possibility of similar pollution sources；the main sources of pollution for Pb， Cr， and As were mining emissions，the main sources of pollution for Cu and Zn were agricultural non-point source pollution，and the main source of Ni pollution was natural sources. On the basis of qualitative analysis of pollution sources， the contribution rates of different pollution sources were quantitatively confirmed through the PCA-MLR model. The main pollution sources in X County and A Town were agricultural non-point source pollution. Compared with traditional methods， this method had higher accuracy and practicality in pollution source analysis.

Key words： PCA-MLR model； farmland in mining areas； soil； heavy metals； source of pollution

土壤重金屬污染是指土壤中重金屬元素含量超過環(huán)境的承載能力，從而對生態(tài)環(huán)境和人類健康造成危害的現(xiàn)象。土壤重金屬污染的成因較復(fù)雜，有自然因素、人為因素以及兩者疊加的因素。區(qū)域大尺度上自然因素是主要成因，在特定區(qū)域等局部范圍內(nèi)，人為因素是主要成因。耕地土壤重金屬成分分析是耕地土壤污染治理的前提和基礎(chǔ)。

于旦洋等［1］闡述了中國土壤重金屬污染源解析的研究進展，主要方法包括主成分分析法（PCA）、正定矩陣因子分解法（PMF）、絕對因子得分-多元線性回歸法（APCS-MLR）以及單一和多參數(shù)同位素示蹤法等。劉宏波等［2］的研究指出，定量源解析技術(shù)主要包括源清單法、擴散模型法、化學(xué)質(zhì)量平衡模型法、正定矩陣因子分解法、UNMIX模型法和同位素法等。地理探測器（Geo detector）被廣泛用于分析土壤重金屬空間分異及其影響因素，龔倉等［3］利用地理探測器分析土壤重金屬的空間分異及其影響因素。趙嘉麗等［4］基于AWRSEI模型利用地理探測器分析岱海流域生態(tài)環(huán)境質(zhì)量的時空演變及其驅(qū)動因子。本研究采用PCA-MLR（主成分-多元線性回歸）模型分析礦區(qū)耕地土壤的重金屬來源。

1 數(shù)據(jù)來源與方法

1.1 數(shù)據(jù)來源

采樣點位于鉛鋅礦區(qū)附近的X縣，共設(shè)有1 078個耕地采樣點，每個采樣點均進行3次獨立的土壤樣本采集。

1.2 污染評價標(biāo)準(zhǔn)

單因子指數(shù)法只能反映單一污染物的污染程度，當(dāng)土壤樣本量大、污染元素較多時，單因子指數(shù)法無法準(zhǔn)確反映土壤污染情況。因此本研究采取內(nèi)梅羅綜合污染指數(shù)法來對污染狀況進行評價。

[P綜=（Pimax）2+i=1nPin22] （1）

式中，[P綜]為土壤重金屬的內(nèi)梅羅綜合污染指數(shù)；[Pimax]為單因子污染指數(shù)的最大值；[Pi]為土壤中重金屬元素i的污染指數(shù)；[n]表示監(jiān)測的污染物種類數(shù)量。

內(nèi)梅羅綜合污染指數(shù)的評價標(biāo)準(zhǔn)如表1所示。

1.3 主成分分析

主成分分析法是將樣本間差距較大的指標(biāo)抽取出來，通過對變量進行轉(zhuǎn)換的方法實現(xiàn)數(shù)據(jù)降維，是一種可行的重金屬污染源解析方法，可以提取出主要因子并保留原始數(shù)據(jù)信息［5］。主成分分析需要考慮因子之間的相關(guān)性［6］，因此需要進行相關(guān)性和聚類分析。

1.4 多元線性回歸

主成分分析的結(jié)果與污染源的貢獻相關(guān)，但是只是定性推斷潛在的污染源，直接用于源解析效果不好，PCA-MLR模型是基于主成分分析的結(jié)果，用多元線性回歸定量確認(rèn)不同污染源的貢獻率。

[G=QnQn×100%] （2）

式中，G為貢獻率；[Qn]為回歸系數(shù)。

2 結(jié)果與分析

2.1 土壤重金屬分析

由圖1可知， X縣土壤重金屬污染的總體情況良好，但存在輕度污染、中度污染的風(fēng)險。

本次采樣共涉及21個鄉(xiāng)鎮(zhèn)，存在污染采樣點的鄉(xiāng)鎮(zhèn)有13個，其中輕度污染采樣點的鄉(xiāng)鎮(zhèn)有10個，中度污染采樣點的鄉(xiāng)鎮(zhèn)有3個，從中隨機選取6個鄉(xiāng)鎮(zhèn)進行分析，如表2所示。在X縣的土壤重金屬污染檢測中，A鎮(zhèn)和E鄉(xiāng)的重金屬含量較高，且A鎮(zhèn)同時具有中度污染和輕度污染的采樣點。因此，本研究以A鎮(zhèn)為例進一步分析可能的污染來源。

2.2 土壤重金屬來源分析

2.2.1 土壤重金屬相關(guān)性分析 為了深入理解X縣土壤中不同重金屬之間的潛在聯(lián)系，本研究進行了皮爾森相關(guān)性分析。皮爾森相關(guān)系數(shù)的熱力圖顏色越深，數(shù)值越大，表示土壤各重金屬之間相關(guān)性越強。由圖2可知，Cd和Cu、Zu之間的相關(guān)系數(shù)分別為0.222、0.186，相關(guān)性較弱。Zn和Cu之間有明顯的相關(guān)性，相關(guān)系數(shù)為0.578。As與Cr、Pb之間有明顯的相關(guān)性，相關(guān)系數(shù)分別為0.588、0.418，而Cr與Pb之間的相關(guān)系數(shù)也達0.484，Pb、As、Cr可能存在相似的污染來源。Ni與Cu、Zn之間的相關(guān)系數(shù)分別為0.127、0.235，相關(guān)性較弱。

2.2.2 土壤重金屬聚類分析 使用組間連接法對A鎮(zhèn)土壤數(shù)據(jù)進行系統(tǒng)聚類，進一步對土壤中重金屬元素進行分析。由圖3可知，Pb、Cr、As聚為一類，Cu、Zn、Cd聚為一類，Ni單獨聚為一類，聚類結(jié)果和相關(guān)性分析結(jié)果基本吻合。

2.2.3 KMO檢驗和Bartlett球形檢驗 根據(jù)相關(guān)性分析和聚類分析的結(jié)果可知，X縣土壤重金屬含量之間存在一定的信息交叉。將數(shù)據(jù)進行KMO檢驗和Bartlett球形檢驗，結(jié)果如表3所示，KMO為0.613，通過KMO檢驗，說明各變量之間存在相關(guān)性，符合主成分分析要求。Bartlett球形檢驗的結(jié)果顯示，Plt;0.01，通過Bartlett檢驗，可進行主成分分析。

2.2.4 土壤重金屬主成分分析 對特征值大于1.000的3個主成分進行提取，3個主成分累計貢獻率為68.896%，如表4所示，使用因子載荷法進一步分析，得到土壤重金屬的因子載荷矩陣熱力圖，如圖4所示。

根據(jù)各主成分的特征值以及主成分得分系數(shù)矩陣（表5），得到3個主成分與變量指標(biāo)間的線性關(guān)系。

[F1=0.002X1+0.431X2+0.376X3-0.092X4-]

[0.052X5-0.014X6+0.399X7]

[F2=0.255X1-0.024X2+0.070X3+0.045X4+0.469X5+0.219X6+0.131X7]

[F3=0.556X1-0.021X2+0.176X3-0.128X4+0.068X5+0.781X6-0.133X7]

表5 主成分得分系數(shù)矩陣

[指標(biāo) 主成分 1 2 3 Cd（X1） 0.002 0.255 0.556 Cr（X2） 0.431 -0.024 -0.021 Pb（X3） 0.376 0.070 0.176 Cu（X4） -0.092 0.045 -0.128 Zn（X5） -0.052 0.469 0.068 Ni（X6） -0.014 0.219 0.781 As（X7） 0.399 0.131 -0.133 ]

主成分1貢獻了28.676%的信息。Cr、Pb、As的因子載荷分別為0.866、0.755、0.801，遠(yuǎn)高于其他重金屬元素的因子載荷。Pb、As的異常受人為影響嚴(yán)重［7，8］，Pb的常見污染源為燃煤、金屬冶煉排放、尾氣排放等［9-11］。As主要受地質(zhì)背景影響，馮新斌等［12］指出Pb和As經(jīng)金屬冶煉排放后，可通過大氣沉降污染土壤。Cr含量受成土母質(zhì)控制［13］，采礦區(qū)的礦冶活動對周圍土壤Cr含量有一定的影響［14］?？紤]到Pb、Cr、As具有較高相關(guān)性，很可能存在相似的來源，而A鎮(zhèn)有大型的鉛鋅礦區(qū)，因此，主成分1的來源是礦業(yè)排放。

主成分2貢獻了25.775%的信息。Cu、Zn和Cd的因子載荷系數(shù)分別為0.804、0.847、0.460，說明這3種元素可能存在相似的污染源。通過查閱文獻可知，含Zn、Cu的殺菌劑在農(nóng)業(yè)生產(chǎn)中被廣泛使用［15］，而含銅農(nóng)藥是當(dāng)前農(nóng)業(yè)銅污染的主要來源［16，17］。Cd常作為磷礦石雜質(zhì)隨磷肥施于土壤中［18］，同時也應(yīng)用于地膜生產(chǎn)。有研究表明，Cd常被視為農(nóng)業(yè)行為的標(biāo)志元素［19］。王美等［20，21］認(rèn)為過磷酸鈣中Zn、Cu、Cd含量明顯高于氮肥、鉀肥；有機肥會提高Zn、Cu、Cd在土壤中的含量；Zn、Cu在高肥、高復(fù)種率的種植模式下，容易隨著禽畜糞便等有機肥進入耕地［22］。A縣是農(nóng)業(yè)縣，有地膜、農(nóng)藥、化肥、有機肥等使用背景，因此主成分2的來源是農(nóng)業(yè)面源污染。

主成分3貢獻了14.445%的信息。Ni的因子載荷最高，為0.790。有文獻表明，Ni主要來源為成土母質(zhì)［23-25］，受人為影響較小。Ni與其他重金屬的相關(guān)性較低，有孤立存在的特點。因此將主成分3的來源解釋為自然來源。

3 基于PCA-MLR模型的污染來源解析

3.1 多元線性回歸

將7種重金屬元素的濃度標(biāo)準(zhǔn)化，自變量選取" 3個主成分的公因子（F1、F2和F3），構(gòu)建PCA-MLR模型。模型的相關(guān)系數(shù)（R2）為0.761，說明模型的擬合度較好。模型的顯著性為0.000，小于0.05，假設(shè)成立。

3.2 污染來源解析結(jié)果

圖5為X縣各污染源的貢獻率，X縣礦業(yè)（礦區(qū)、礦業(yè)）排放污染貢獻率為1.76%，自然來源污染的貢獻率為25.69%；農(nóng)業(yè)面源污染貢獻率最高，為72.55%，說明長期的農(nóng)業(yè)生產(chǎn)活動，如有機肥、化肥、農(nóng)藥和地膜的不合理使用等，可能對該區(qū)域造成一定的污染風(fēng)險。

3.3 A鎮(zhèn)各采樣點農(nóng)業(yè)面源污染的貢獻率

A鎮(zhèn)有鉛鋅礦區(qū)，根據(jù)前人針對鉛鋅礦區(qū)周邊土壤重金屬污染的相關(guān)研究［26-28］，礦區(qū)周邊土壤Cd、Zn、Pb容易存在不同程度的富集。為了排除礦區(qū)、礦業(yè)等因素影響，單獨計算農(nóng)業(yè)面源污染的貢獻率，如圖6所示。采樣點tr11、tr15、tr16、tr21、tr298的農(nóng)業(yè)面源污染貢獻率高于其他采樣點，且這5個采樣點都遠(yuǎn)離礦區(qū)、礦業(yè)。

土壤重金屬農(nóng)業(yè)面源污染中，X縣處于輕度污染的采樣點有20個，中度污染的采樣點有3個，其中A鎮(zhèn)處于輕度污染的采樣點有4個，中度污染的采樣點有1個，如表6所示。

鉛鋅礦區(qū)、工廠/冶煉廠、居民區(qū)、公路距離5個污染點最近的距離如表7所示。5個污染點距離鉛鋅礦區(qū)和工廠/冶煉廠均較遠(yuǎn)，且位于礦區(qū)和工廠/冶煉廠的上風(fēng)口，因此認(rèn)為受礦區(qū)、工廠/冶煉廠影響的可能性很小。5個污染點距離居民區(qū)和公路較近，研究顯示［29］，交通運輸導(dǎo)致的污染物為Cr、Cu、Pb、Zn、Cd等，根據(jù)全國土壤污染調(diào)查公報［30］，這些污染物集中在公路兩側(cè)150 m的范圍內(nèi)。羅婭君等［31］的研究發(fā)現(xiàn)，距公路的距離越遠(yuǎn)，受污染土壤的重金屬含量越低，因此交通運輸導(dǎo)致土壤重金屬污染的可能性較小。通過查詢流經(jīng)A鎮(zhèn)的支流水質(zhì)狀況監(jiān)測月報，A鎮(zhèn)周邊的水域水質(zhì)良好，不太可能受到灌溉污染。陳雅麗等［32］對全國土壤重金屬污染特征與來源研究做了梳理分析，發(fā)現(xiàn)土壤中的Cd、Cu、Zn主要源于施肥等農(nóng)業(yè)活動。A鎮(zhèn)耕地存在重金屬污染風(fēng)險，重金屬污染風(fēng)險的主要來源為農(nóng)業(yè)面源污染，并且可能是化肥、有機肥、農(nóng)藥殺菌劑污染。PCA-MLR模型解析結(jié)果可靠，為當(dāng)?shù)氐母赝寥牢廴痉乐翁峁┝祟A(yù)警。

本研究利用內(nèi)梅羅綜合污染指數(shù)評價X縣的污染狀況，發(fā)現(xiàn)X縣土壤總體情況良好，但存在輕度污染、中度污染的風(fēng)險；本次采樣共涉及21個鄉(xiāng)鎮(zhèn)，存在污染采樣點的鄉(xiāng)鎮(zhèn)有13個。

4 小結(jié)與討論

本研究利用內(nèi)梅羅綜合污染指數(shù)評價X縣的污染狀況，發(fā)現(xiàn)土壤總體情況良好，但存在輕度污染、中度污染的風(fēng)險；本次采樣共涉及21個鄉(xiāng)鎮(zhèn)，存在污染采樣點的鄉(xiāng)鎮(zhèn)有13個。

相關(guān)性分析和聚類分析表明，Cr、As、Pb三者存在明顯相關(guān)性，Cu、Zn二者存在明顯相關(guān)性，可能存在相似的污染源。Pb、Cr、As的污染源主要為礦業(yè)排放，Cu、Zn、Cd的污染源主要為農(nóng)業(yè)面源污染，Ni的污染源為主要為自然來源。

在對污染源定性分析的基礎(chǔ)上，通過PCA-MLR模型定量確認(rèn)不同污染源的貢獻率。X縣礦業(yè)（礦區(qū)、礦業(yè)）排放污染貢獻率為1.76%，自然來源污染的貢獻率為25.69%，農(nóng)業(yè)面源污染貢獻率最高，為72.55%。X縣化肥使用量為854.5 kg/hm2，遠(yuǎn)高于國內(nèi)化肥的施量標(biāo)準(zhǔn)（不超過250 kg/hm2）［33］；農(nóng)藥使用量為10.2 kg/hm2，存在明顯污染風(fēng)險；地膜殘留量較低，但回收率不高，存在污染風(fēng)險；秸稈利用率僅70%左右。綜上，X縣的主要污染源為農(nóng)業(yè)面源污染。

在X縣的土壤重金屬污染檢測中，A鎮(zhèn)的重金屬含量較高，且A鎮(zhèn)同時具有中度污染和輕度污染的采樣點。因此，本研究以A鎮(zhèn)為例進一步分析可能的污染來源。A鎮(zhèn)處于輕度污染的采樣點有4個，中度污染的采樣點有1個。通過排除交通運輸、灌溉等污染源，發(fā)現(xiàn)A鎮(zhèn)土地重金屬污染風(fēng)險的主要來源為農(nóng)業(yè)面源污染，并且可能是由于化肥、有機肥、農(nóng)藥殺菌劑污染導(dǎo)致。PCA-MLR模型解析結(jié)果可靠，為當(dāng)?shù)氐耐寥牢廴痉乐翁峁﹨⒖肌?/p>

參考文獻：

［1］ 于旦洋，王顏紅，丁 茯，等.近十年來我國土壤重金屬污染源解析方法比較［J］.土壤通報，2021，52（4）：1000-1008.

［2］ 劉宏波，瞿明凱，張健琳，等.土壤污染物源解析技術(shù)研究進展［J］.環(huán)境監(jiān)控與預(yù)警，2021，13（1）：1-6，19.

［3］ 龔 倉，王順祥，陸海川，等.基于地理探測器的土壤重金屬空間分異及其影響因素分析研究進展［J］.環(huán)境科學(xué)，2023，44（5）：2799-2816.

［4］ 趙嘉麗，李 興，孫 冰.基于AWRSEI的岱海流域生態(tài)環(huán)境質(zhì)量時空演變及驅(qū)動因子分析［J］.環(huán)境科學(xué)，2024，45（3）：1598-1614.

［5］ FACCHINELLI A，SACCHI E，MALLEN L.Multivariate statistical and GIS-based approach to identify heavy metal sources in soils［J］.Environmental pollution，2001，114（3）：313-324.

［6］ 袁 宏，鐘紅梅，趙 利，等.基于PCA/APCS受體模型的崇州市典型農(nóng)田土壤重金屬污染源解析［J］.四川環(huán)境，2019，38（6）：35-43.

［7］ 韓志軒，王學(xué)求，遲清華，等.珠江三角洲沖積平原土壤重金屬元素含量和來源解析［J］.中國環(huán)境科學(xué)，2018，38（9）：3455-3463.

［8］ 李小平，高 瑜，張 蒙，等.城市土壤重金屬空間分布、污染與來源［J］.環(huán)境科學(xué)與技術(shù)，2018，41（6）：138-146.

［9］ 鄭宇娜，劉 鵬，劉金河，等.臺州市典型電子垃圾拆解場地周邊農(nóng)田土壤重金屬污染特征和來源解析［J］.農(nóng)業(yè)環(huán)境科學(xué)學(xué)報，2022，41（7）： 1442-1451.

［10］ YANG S，HE M，ZHI Y，et al.An integrated analysis on source-exposure risk of heavy metals in agricultural soils near intense electronic waste recycling activities［J］.Environment international，2019，133：105239.

［11］ BI X Y，LI Z G，WANG S X，et al.Lead isotopic compositions of selected coals， Pb/Zn ores and fuels in China and the application for source tracing［J］.Environmental science technology，2017，51（22）：13502-13508.

［12］ 馮新斌，倪建宇，洪業(yè)湯，等.貴州省煤中揮發(fā)性和半揮發(fā)性微量元素分布規(guī)律的初步研究［J］.環(huán)境化學(xué)，1998（2）：148-153.

［13］ LIU L L，LIU Q Y，MA J，et al.Heavy metal（loid）s in the topsoil of urban parks in Beijing，China： Concentrations， potential sources，and risk assessment［J］.Environmental pollution，2020，260：114083.

［14］ 黃 赫，周 勇，劉宇杰，等.基于多源環(huán)境變量和隨機森林的農(nóng)用地土壤重金屬源解析——以襄陽市襄州區(qū)為例［J］.環(huán)境科學(xué)學(xué)報，2020，40（12）：4548-4558.

［15］ 石寧寧. 不同類型農(nóng)藥園區(qū)周邊農(nóng)田土壤—作物重金屬分布特征及評價研究［D］.南京：南京農(nóng)業(yè)大學(xué)，2010.

［16］ 王 萌，李杉杉，李曉越，等.我國土壤中銅的污染現(xiàn)狀與修復(fù)研究進展［J］.地學(xué)前緣，2018，25（5）：305-313.

［17］ 陳 航，王 穎，王 澍.銅山礦區(qū)周邊農(nóng)田土壤重金屬來源解析及污染評價［J］.環(huán)境科學(xué)，2022，43（5）：2719-1372.

［18］ LV J S，LIU Y，ZHANG Z L，et al.Identifying the origins and spatial distributions of heavy metals in soils of Ju country （Eastern China） using multivariate and geostatistical approach［J］.Journal of soils and sediments，2015，15（1）：163-178.

［19］ 穆德苗，陳艷秋，胡 濤，等.基于田塊尺度的農(nóng)田土壤重金屬污染評價及來源解析［J］.農(nóng)業(yè)環(huán)境科學(xué)學(xué)報，2022， 41（6）：1271-1283.

［20］ 王 美，李書田，馬義兵，等.長期不同施肥措施對土壤和作物重金屬累積的影響［J］.農(nóng)業(yè)環(huán)境科學(xué)學(xué)報，2014，33（1）：63-74.

［21］ 王 美，李書田.肥料重金屬含量狀況及施肥對土壤和作物重金屬富集的影響［J］.植物營養(yǎng)與肥料學(xué)報，2014，20（2）：466-480.

［22］ 董達誠，胡夢淩，羅 昱，等.陸良縣菜地土壤重金屬的污染評價及來源解析［J］.環(huán)境污染與防治，2021，43（6）：732-737，790.

［23］ 湛潤生，胡冬梅，甄莉娜，等.山西省天鎮(zhèn)縣設(shè)施菜地土壤重金屬污染評價及其源解析［J］.環(huán)境污染與防治，2021，43（12）：1573-1577.

［24］ 張 佳. 地質(zhì)高背景區(qū)農(nóng)田土壤重金屬污染來源辨析與防治分區(qū)［D］.北京：中國地質(zhì)大學(xué)（北京），2019.

［25］ 陳丹丹，譚 璐，聶紫萌，等.湖南典型金屬冶煉與采選行業(yè)企業(yè)周邊土壤重金屬污染評價及源解析［J］.環(huán)境化學(xué)，2021，40（9）：2667-2679.

［26］ 顧 會，趙 濤，高 月，等.貴州省典型鉛鋅礦區(qū)土壤重金屬污染特征及來源解析［J］.地球與環(huán)境，2022，50（4）：506-515.

［27］ 袁鑫奇，俞乃琪，郭兆來，等.會澤鉛鋅礦區(qū)廢棄地優(yōu)勢草本植物的重金屬富集特征［J］.生態(tài)與農(nóng)村環(huán)境學(xué)報，2022，38（3）：399-408.

［28］ 楊 寧，李東海，楊小波，等.鉛鋅礦區(qū)周邊土壤重金屬污染及植物富集特征［J］.熱帶生物學(xué)報，2021，12（4）：500-507.

［29］ 于向楠. 北京市道路塵中重金屬元素地球化學(xué)特征的時空變化規(guī)律［D］.北京：中國地質(zhì)大學(xué)（北京），2021.

［30］ 龐榮麗，王瑞萍，謝漢忠，等.農(nóng)業(yè)土壤中鎘污染現(xiàn)狀及污染途徑分析［J］.天津農(nóng)業(yè)科學(xué)，2016，22（12）：87-91.

［31］ 羅婭君，王照麗，張 露，等.成綿高速公路兩側(cè)土壤中4種重金屬的污染特征及分布規(guī)律［J］.安全與環(huán)境學(xué)報，2014，14（3）：283-287.

［32］ 陳雅麗，翁莉萍，馬 杰，等.近十年中國土壤重金屬污染源解析研究進展［J］.農(nóng)業(yè)環(huán)境科學(xué)學(xué)報，2019，38（10）：2219-2238.

［33］ 李明飛，張海瑩，劉向華，等.黃河流域化肥施用強度區(qū)域分布及環(huán)境風(fēng)險評價［J］.河南大學(xué)學(xué)報（自然科學(xué)版），2022，52（2）：146-153.

宋 科，王 皓，楊 偉，等. 基于PCA-MLR的礦區(qū)耕地土壤重金屬污染來源分析［J］. 湖北農(nóng)業(yè)科學(xué)，2024，63（11）：185-190．

收稿日期：2023-06-02

基金項目：云南省重大科技專項（202202AE090066）

作者簡介：宋 科（1999-），男，湖南衡陽人，碩士，主要從事農(nóng)業(yè)工程與信息技術(shù)研究，（電話）17810795056（電子信箱）sk_9906@foxmail.com；共同第一作者，王 皓（1999-），男，重慶人，碩士，主要從事農(nóng)業(yè)工程與信息技術(shù)研究，（電話）15608819594（電子信箱）renhaowang2@gmail.com；通信作者，曹志勇（1976-），男，云南會澤人，教授，主要從事農(nóng)業(yè)環(huán)境信息處理、圖像采集與處理研究，（電話）13708718797（電子信箱）hi.czy@qq.com。