浙江省某尾礦庫周邊農(nóng)田土壤重金屬
污染特征及來源解析
呂玉娟1,，王秋月1,2，孫雪梅1,2,，張志偉1，張毅敏1,，高月香1*
1.生態(tài)環(huán)境部南京環(huán)境科學研究所
2.河海大學環(huán)境學院
摘要：尾礦庫周邊生態(tài)環(huán)境安全受到高度關注,，土壤重金屬污染是農(nóng)田治理和保護的風險源之一。以地處浙東丘陵山地的浙江省某銅礦尾礦庫周邊農(nóng)田為研究對象,，測定了農(nóng)田土壤中8種重金屬元素Cd,、Hg、As,、Pb,、Zn、Cu,、Cr,、Ni的濃度，運用地累積指數(shù)法,、污染指數(shù)法,、潛在生態(tài)風險指數(shù)法和生態(tài)風險預警指數(shù)法對農(nóng)田土壤重金屬污染程度以及生態(tài)風險進行評價，結(jié)合正定矩陣受體模型（PMF）,，定量解析農(nóng)田土壤重金屬的來源,。結(jié)果表明：1）（略）農(nóng)田土壤中Cd、Hg,、Cu,、Zn濃度分別是土壤元素背景值的5.36、2.06,、8.19,、5.36倍，具有高度變異性,；污染指數(shù)評價結(jié)果表明,，Cu、Zn,、Cd重度污染占比均達到10.5%,，中度污染占比為5.26%，靠近尾礦庫（<300m）的15.8%的點位處于重度污染等級,；地累積指數(shù)評價結(jié)果表明,，Cd、Cu,、Zn和Hg可能具有累積風險,。2）潛在生態(tài)風險評價結(jié)果表明,，Cd為很強生態(tài)風險，Hg為較強生態(tài)風險,，Cu為中等生態(tài)風險,，其余重金屬均為輕微風險；綜合潛在生態(tài)風險指數(shù)（RI）為308.91,，綜合潛在風險為較強風險,。生態(tài)風險預警評估結(jié)果表明，Cu為重警,，Cd和Zn為中警,，Hg為輕警，As為預警,，Pb,、Cr和Ni為無警；綜合生態(tài)風險預警指數(shù)（IER）為16.06,，綜合生態(tài)風險預警為重警,。RI和IER空間分布基本一致，主要受Cd,、Cu,、Zn和Hg的影響。3）PMF解析出3個源,，Cd,、Zn、Cu主要受銅礦尾礦庫尾砂和壩下滲水的混合源影響,，貢獻率分別為94.4%,、94.3%和67.1%；Hg可能是以肥料,、農(nóng)藥施用等農(nóng)業(yè)活動源為主,，貢獻率為61.5%；Cr,、Ni,、Pb和As主要受成土母質(zhì)和交通運輸活動混合源的影響，貢獻率分別為89.7%,、82.7%,、75.0%和68.3%。
關鍵詞：尾礦庫,；土壤重金屬,；污染評價；潛在生態(tài)風險,；正定矩陣受體模型（PMF）,；源解析
礦產(chǎn)資源是人類生存和發(fā)展的重要物質(zhì)基礎，也是人類生產(chǎn)資料和生活資料的重要來源,。我國是世界上為數(shù)不多的礦產(chǎn)資源總量豐富且礦種齊全的大國之一,，現(xiàn)已發(fā)現(xiàn)礦產(chǎn)173種。隨著礦產(chǎn)資源不斷被開發(fā)利用,，國民經(jīng)濟得到極大發(fā)展[1-2],。然而，在礦山開采過程中,，必然會產(chǎn)生Hg,、As、Pb,、Cd等成礦或伴礦元素污染,，其通過廢液排放、廢渣堆積,、雨水和地表水淋濾等途徑污染周邊土壤[3],。土壤重金屬具有持久性、累積性,、隱蔽性和循環(huán)性等特點,，因難以被降解而長期穩(wěn)定存在于土壤中，可通過食物鏈在人體聚集,，（略）功能和結(jié)構,，威脅人體健康安全[4-7]。
開（略）周邊土壤重金屬污染評價與污染源識別對（略）周邊土壤修復與復墾具有重要意義,。地累積指數(shù)法,、污染指數(shù)法、潛在生態(tài)風險指數(shù)法等方法是當前應用較多的重金屬污染評價方法[8-9],。除（略）域背景值作為參比外,，以GB15618—2018《土壤環(huán)境質(zhì)量農(nóng)用地土壤污染風險管控標準（試行）》風險篩選值作為參比值[10-13]成為國內(nèi)主流評價方法。土壤重金屬主要來自于自然本底（與地殼中微量元素有關）和人類活動（工業(yè)生產(chǎn),、農(nóng)業(yè)活動等）[14],，常用的土壤重金屬分布與源識別方法包括地統(tǒng)計與空間插值方法（如普通克里格法、反距離加權法）,，多元統(tǒng)計分析（如主成分分析,、相關分析和聚類分析），定量方法〔如穩(wěn)定同位素法,、Unmix模型和正定矩陣受體模型（PMF）〕等,。與其他受體模型相比，PMF模型可以同時將測量不確定性和非負性約束納入計算過程[15],，且可以自動處理缺失和錯誤的數(shù)據(jù),，被廣泛應用于土壤污染物源解析中,。Zhang等[16]采用PMF模型對（略）鉛（略）土壤重金屬來源進行解析，發(fā)現(xiàn)人為礦業(yè)活動（黃鐵礦和鉛鋅礦）,、土壤母質(zhì)（玄武巖和碳酸鹽巖）和農(nóng)業(yè)生產(chǎn)活動是周邊土壤重金屬的主要來源,。朱曉麗等[17]基于PMF模型對寶雞鉛鋅尾礦庫周邊農(nóng)田進行土壤重金屬源解析發(fā)現(xiàn)，周邊土壤重金屬主要來源于農(nóng)業(yè)活動源,、自然源,、化石燃燒釋放源、鉛鋅選礦源,、鉛鋅工業(yè)活動源和交通源等,。對于已閉庫（略），人為采礦活動基本終止,，對周邊土壤的重金屬污染程度因后期管理水平差異而不同,，而在沒有人為采礦活動源時，對重金屬其他來源的解析成為后（略）土壤修復或復墾的重要依據(jù),。
因此,，筆者以浙江省某閉庫尾礦庫下游農(nóng)田土壤為研究對象，對農(nóng)田土壤重金屬濃度和空間分布特征進行分析,，使用地累積指數(shù)法,、污染指數(shù)法和潛在生態(tài)風險指數(shù)法等對農(nóng)田土壤環(huán)境質(zhì)量進行評價，并結(jié)合PMF模型,，定量解（略）農(nóng)田土壤重金屬潛在來源,，以期為閉庫尾礦庫周邊農(nóng)田土壤復墾及重金屬修復提供一定的理論依據(jù)。
1.（略）與研究方法
1.1（略）概況
該銅礦尾礦庫地處浙江省中北部,、杭州灣南岸的（略）,，屬于亞熱帶季風氣候，年平均氣溫16.4℃,，年均降水量1400mm,，盛行東北偏東風。礦區(qū)所在地屬浙西山地丘陵,、浙東丘陵山地和浙北平原三大地貌單元的交接地帶,，地勢南高北低，形成群山環(huán)繞,、盆地內(nèi)含,、平原集中的地貌特征，最高海拔為1194.6m,，最低海拔僅為3.1m,。礦區(qū)主要盛產(chǎn)銅精礦、鋅精礦、硫精礦,，伴生的金,、銀冶煉回收。尾礦庫于1982年設計建設,，1987年11月投入使用,，服務年限36年，為傳統(tǒng)上游式濕排山谷型尾礦庫,，設計總壩高85m、總庫容343萬m3,，為三等庫,。1992年該尾礦庫因尾砂輸送管道爆管，尾砂溢出導致壩下52000m2農(nóng)田受到不同程度污染,，2020年7月起尾礦庫閉庫,。目前，該尾礦庫下游52000m2農(nóng)田正進行生態(tài)修復工作,，距離尾礦庫較遠的農(nóng)田上有農(nóng)戶自發(fā)種植的零散蔬菜,，其他大（略）域閑置有雜草生長。
通過實地調(diào)查,，選取該閉庫尾礦庫周邊農(nóng)田作為研究對象,。所選的尾礦庫由三面山坡自然形成的近南北走向山谷組成，南面筑壩圍庫,，在壩底修建有沉降池,，尾礦庫周圍多為巖石裸露，中下部植被較好,，以灌木林為主,，高大喬木較少，庫邊滲水渠有尾礦庫滲漏水流經(jīng),，岸邊為種植有零散蔬菜的農(nóng)田（圖1）,。
1.2樣品采集處理
（略）樣品采集
在尾礦庫下游約52000m2的農(nóng)田范圍內(nèi)，距離尾礦庫由遠及近共采集農(nóng)田土壤表層（0~20cm）樣品19個,。（略）戶外助手APP中GPS（WGS-84坐標）對每個樣點進行定位,，使用五點采樣法采集耕層土壤樣品，先去掉表面雜物和礫石,，使用鐵鍬采集0~20cm土壤（用鐵鏟挖掘后,，用木片刮去與金屬采樣器接觸的部分，再用木片采取樣品）,，充分混合后取約1kg（略）帶回實驗室,。樣品采集點位見圖2。
（略）測定方法
將采集的樣品放置于陰涼的通風處自然風干,，剔除根莖,、石頭等雜物,，采用瑪瑙研缽研磨后過100目篩，隨后裝入統(tǒng)一規(guī)格的密封袋中保存待測,。土壤中重金屬（Cd,、Hg、As,、Pb,、Zn、Cu,、Cr和Ni）濃度測定參照HJ/T166—2004《土壤環(huán)境監(jiān)測技術規(guī)范》,，其中Cd、Pb,、Zn,、Cu、Cr和Ni元素提取采用王水消解法〔V（濃鹽酸）∶V（濃硝酸）=3∶1〕,，并采用電感耦合等離子體質(zhì)譜儀（STS-188,，中國）測定其濃度，Hg和As濃度采用王水消解法和原子熒光光度計（STS-041,，中國）測定,。所有試劑均為優(yōu)級純，采用國家標準土壤樣品（GSS-25）進行質(zhì)量控制,，元素回收率控制在90%~105%,，測量結(jié)果的相對誤差在±5%以內(nèi)。
1.3研究方法（略）重金屬污染評價方法
（略）地累積指數(shù)法
地累積指數(shù)法最早由德國學者Müller提出并被用于沉積物中重金屬污染評價,。該指數(shù)不僅考慮了人為污染因素和環(huán)境地球化學背景值,，還特別考慮了自然成巖作用對背景值的影響，被廣泛用于大氣沉降,、土壤和現(xiàn)代沉積物中重金屬的污染評價,。其公式如下：
（略）單因子污染指數(shù)法
單因子污染指數(shù)是土壤表層重金屬污染物濃度與污染物評價標準的比值，能夠準確反映某一位點單項重金屬污染狀況,。其計算公式如下：
（略）內(nèi)梅羅綜合指數(shù)法
內(nèi)梅羅綜合指數(shù)法是當前國內(nèi)外進行綜合污染指數(shù)計算最常用的方法之一,，涵蓋各單項污染指數(shù)，并突出高濃度污染物在評價結(jié)果中的權重,，用于評價土壤復合污染狀況,。其計算公式如下：
（略）重金屬污染風險評價法
（略）潛在生態(tài)風險指數(shù)法
潛在生態(tài)風險指數(shù)法由Hakanson提出，主要應用于土壤和沉積物中重金屬的污染生態(tài)風險評價,，該方法不僅考慮了重金屬濃度,，還將重金屬的生態(tài)效應、環(huán)境效應和毒理學效應聯(lián)系起來。其計算公式如下：
（略）生態(tài)風險預警指數(shù)
生態(tài)風險預警是指對自然資源或生態(tài)風險可能出現(xiàn)的衰竭或危機而建立的報警,，其計算公式如下：
（略）重金屬污染源解析方法
PMF模型是一種基于受體模型對污染源進行矩陣分析的定量方法,，該模型在求解過程中對因子載荷和因子得分均做非負約束，確保因子矩陣分解的結(jié)果中不出現(xiàn)負值,，使得到的源成分譜與源貢獻率具有更實際的物理意義[24-26],。該模型將原始矩陣Xab因子化，分解為2個因子矩陣fkb和gak以及1個殘差矩陣eab,，計算公式如下：
PMF模型需要通過多次迭代計算分解原始矩陣,，最終得出最優(yōu)化的污染源貢獻因子矩陣G和污染源因子矩陣F，使得目標函數(shù)Q最小化,。目標函數(shù)Q定義如下：
當各個重金屬元素含量小于或等于相應的方法檢出限（MDL）時,，不確定度uab的計算公式為：
當各個重金屬元素含量大于相應的MDL時，不確定度uab的計算公式如下：
式中EF為相對標準偏差,。
1.4數(shù)據(jù)處理采用Excel2019和Origin2022軟件進行數(shù)據(jù)處理分析，采用SPSS16.0軟件進行相關性分析和主成分分析,，采用ArcGIS10.3軟件進行空間插值制圖,，源解析基于EPAPMF5.0模型完成。
2.結(jié)果與討論
2.1農(nóng)田土壤重金屬濃度特征
由描述統(tǒng)計結(jié)果（表3）可知,，尾礦庫周邊農(nóng)田土壤中8種重金屬元素偏度為0.01~3.44,，屬右偏分布。Cd,、Hg,、Cu、Zn濃度均值分別是0.75,、0.31,、172、419mg/kg,，分別是土壤元素背景值的5.36,、2.06、8.19,、5.36倍,。其中，Cd,、Cu,、Zn濃度均值高于風險篩選值，超標率分別為25.0%,、72.7%,、39.7%。與土壤元素背景值和農(nóng)用地土壤污染風險篩選值相比，Cd,、Hg,、As、Pb,、Cu,、Zn6種元素濃度均存在不同程度的超標，其中Cu元素超標率最高,，只有1個采樣點（S15）濃度低于背景值,，其余點位均高于背景值；而所有點位的Cr,、Ni濃度均低于背景值,，說明（略）基本上不存在Ni、Cr污染風險,。
變異系數(shù)是衡量各參數(shù)在數(shù)據(jù)上的變異程度,，空間變異性可以說明自然因子或外界因子的干擾程度。尾礦庫壩下農(nóng)田土壤中各重金屬濃度表現(xiàn)出很大的變異性,。根據(jù)Wilding[27]對變異系數(shù)的分類,，Pb和Ni的變異系數(shù)介于15%~36%，屬中等變異,；Cd,、Hg、As,、Cu和Zn的變異系數(shù)為37%~211%,，屬高度變異。重金屬元素的空間分布具有自然變異性和外在變異性,，自然變異性主要是由于自然界母質(zhì)巖風化釋放導致的,，而外在變異性則是由于（略）域人類活動影響致使重金屬進入土壤中導致的。（略）Cd,、Hg,、As、Cu和Zn的濃度均值都明顯高于（略）土壤元素背景值,，這些重金屬元素在土壤中高度富集,，同時具有較高的變異性，說明它們在（略）域的空間分布極其不均勻,。這有可能是（略）域人為活動所導致的結(jié)果,。與之相反，Cr元素的濃度和變異程度均較低,，說明Cr在土壤中主要來自于自然源,。Pb和Ni元素在（略）域土壤中的濃度均值低于背景值,，卻表現(xiàn)出中度變異性，說明Pb和Ni有可能受到自然因素和人為活動雙重因子的影響,。
空間插值法是采用觀測值預測（略）域的一種空間分析方法,，目前被廣泛應用于生態(tài)環(huán)境治理領域[28]。對（略）19個采樣點Cd,、Hg,、As、Pb,、Cu和Zn濃度進行反距離權重插值并進行掩膜提取,，結(jié)果（圖3）顯示，Cd,、As,、Pb、Cu,、Zn5種重金屬濃度最大值點位均分布在距尾礦庫200m以內(nèi),，呈現(xiàn)靠近尾礦（略）域濃度較高的現(xiàn)象。Cd,、Cu,、Zn大部分點位濃度超過背景值，靠近尾礦庫的3個點位（小于200m）均超過篩選值,，而As、Pb（略）域濃度均超背景值,，但未超篩選值,，As在（略）域濃度也處于背景值和篩選值之間；Hg空間分布呈現(xiàn)距尾礦庫400~650m區(qū)域濃度高的特點,，且處于背景值和篩選值之間,，其中距尾礦庫500~550m區(qū)域濃度最高，這可能與（略）地塊長期零散種植蔬菜有關,。
2.2農(nóng)田土壤重金屬污染評價（略）地累積指數(shù)評價
尾礦庫周邊農(nóng)田土壤中重金屬Igeo評價結(jié)果如圖4所示,。Igeo均值由高到低依次為Cd>Cu>Zn>Hg>As>Pb>Cr>Ni。其中,，As,、Pb、Cr和Ni的Igeo均值均小于0,，總體上處于無污染狀態(tài),，但仍有10.53%的采樣點As為輕度污染狀態(tài)；Hg的Igeo均值小于0,，分別有36.84%,、5.26%和5.26%的采樣點達到了輕度污染,、中度污染和偏重污染；Cd,、Cu和Zn的Igeo均值為0~1,，總體為輕度污染，但有10.53%的采樣點為重度污染,，5.26%的采樣點為中度污染,。
（略）污染指數(shù)評價
使用單因子污染指數(shù)和內(nèi)梅羅綜合指數(shù)對尾礦庫周邊農(nóng)田土壤重金屬污染狀況進行評價（表4）。農(nóng)田土壤重金屬單因子指數(shù)污染程度排列順序為Cu>Zn>Cd>As>Pb>Cr>Hg>Ni,。Cu,、Zn、Cd單因子污染指數(shù)分別為0.19~15.36,、0.12~9.42,、0.12~8.07，重度污染率均為10.5%,，中度污染率均為5.26%,。As、Pb,、Cr,、Hg、Ni單因子指數(shù)均小于1.0,，處于無污染等級,。內(nèi)梅羅綜合污染指數(shù)為0.22~11.2，變異系數(shù)達到197%,，屬于高度變異,，15.8%的點位（S17、S18,、S19）處于重度污染等級,，其他點位無污染，靠近尾礦庫的部分農(nóng)田具有較高的污染風險,。
重金屬Hg單因子污染指數(shù)判斷為無污染,，而（略）域背景值的Igeo評價結(jié)果卻存在中度污染和偏重度污染。綜上可知,，（略）農(nóng)田土壤重金屬Cd,、Cu和Zn存在一定程度的污染，需同時（略）域重金屬Hg的累積污染風險問題,。
2.3農(nóng)田土壤重金屬潛在生態(tài)風險評價（略）潛在生態(tài)風險評價
根據(jù)（略）的土壤元素背景值計算農(nóng)田土壤重金屬潛在生態(tài)風險指數(shù),，結(jié)果見圖5。土壤中8種重金屬元素的單項生態(tài)風險指數(shù)平均值表現(xiàn)為Cd>Cu>Hg>As>Zn>Pb>Ni>Cr,。其中,，As,、Pb、Cr,、Ni和Zn的均值小于40,，屬于輕微生態(tài)風險等級；Cd屬于很強生態(tài)風險,，其中輕微,、中等、較強和極強生態(tài)風險等級樣點占比分別為47.37%,、26.32%,、10.53%和15.79%；Hg屬于較強生態(tài)風險,，其中輕微,、中等、較強,、很強和極強生態(tài)風險等級樣點占比分別為36.84%,、36.84%、15.79%,、5.26%和5.26%,；Cu的生態(tài)風險等級以輕微為主，占比為84.2%,，較強,、很強和極強占比均為5.26%。農(nóng)田土壤重金屬的RI為79.1~1302,，平均值為308.9,，屬于較強生態(tài)風險。其中輕微,、中等、較強和極強生態(tài)風險樣點占比分別為15.79%,、57.89%,、5.26%和21.05%。Cd,、Hg,、Cu元素潛在生態(tài)風險和RI的空間分布顯示，西北部極強（略）主要影響元素是Cd和Cu,，（略）域主要受Hg元素的影響,。
（略）生態(tài)風險預警評估
根據(jù)（略）的土壤元素背景值計算生態(tài)風險預警指數(shù)，結(jié)果見圖6,。IER的平均值表現(xiàn)為Cu>Cd>Zn>Hg>As>Pb>Cr>Ni,，其中Cu為重警,，Cd和Zn為中警，Hg為輕警,，As為預警,，Pb、Cr和Ni為無警,。8種重金屬的IER為?2.56~124.59,，平均值為16.06，屬于重警,。IER處于無警,、預警、輕警,、中警,、重警樣點占比分別為36.84%、10.53%,、21.05%,、5.26%、26.3%,?？臻g分布顯示，（略）生態(tài)風險重警等級分布在西北部和東部,，這與潛在生態(tài)風險指數(shù)的空間分布基本一致,，相比潛在生態(tài)風險指數(shù)，單項生態(tài)風險預警評估發(fā)現(xiàn)Cu在西北部有重警等級,，其單項生態(tài)風險預警評估方面更加敏感,，這與杜貫新等[29]的研究結(jié)果一致。綜上分析,，（略）內(nèi)RI,、IER的空間分布主要受Cu、Zn,、Cd和Hg4項元素的影響,，這與Igeo評價結(jié)果（主要影響因素是Cu、Zn和Cd,，其次是Hg）一致,。
2.4農(nóng)田土壤重金屬污染源解析尾礦庫周邊土壤中重金屬來源途徑比較廣泛，包括人為采選礦活動,、農(nóng)業(yè)生產(chǎn),、交通運輸、大氣沉降和土壤母質(zhì)等,。根據(jù)元素地球化學理論,，相同來源的重金屬元素之間通常存在較顯著的相關關系,，通過對重金屬元素間的統(tǒng)計分析可以闡釋元素間的同源關系，為辨析重金屬來源提供依據(jù)[30],。由（略）農(nóng)田土壤中8種重金屬濃度Pearson相關性分析熱點圖（圖7）可知,，Cd、As,、Pb,、Zn、Ni兩兩之間存在極顯著（P<0.01）正相關性,，說明這5種元素具有較大的同源性和復合污染的可能,。Cu與Cd、Pb,、Zn之間存在極顯著相關性,，推斷Cu與Cd、Pb,、Zn來自同一種污染源的可能性較大,；Cu與As、Cu與Ni及Cr與Hg具有顯著（P<0.05）相關性,，推斷兩兩之間可能具有相似的污染源,；Cr、Hg與其他6種金屬呈現(xiàn)出兩兩不相關性,，說明Cr,、Hg的來源途徑不同。
采用PMF5.0模型對（略）農(nóng)田土壤中8種重金屬來源進行解析,，得到3個重金屬污染源（因子F1,、F2、F3）對8種重金屬的貢獻率（圖8）,。結(jié)果顯示,，因子F1對農(nóng)田土壤重金屬貢獻率為77.62%，且對Cd,、Zn和Cu的貢獻率較高,，分別為94.4%、94.3%和67.1%,。由相關性分析可知，Cd,、Zn和Cu兩兩之間相關性達到極顯著水平（P<0.01）,，說明這3種元素可能來自同一污染源。該銅礦尾礦庫采用傳統(tǒng)濕排法堆存,，導致尾礦中大量的重金屬殘留伴隨酸性礦山廢水向周圍擴散,，造（略）周邊土壤中重金屬元素Cd,、Zn和Cu的復合污染。陳雪等[31]在寧鎮(zhèn)（略）土壤重金屬污染來源解析中發(fā)現(xiàn),，該區(qū)域存在Cd,、Zn和Cu復合污染，與本研究結(jié)果相似,。另外,，陳航等[32]對銅（略）的周圍土壤重金屬來源解析發(fā)現(xiàn)，富含大量Cu元素尾礦和殘渣堆的淋濾擴散（略）周邊土壤Cu污染的主要來源,；杜立宇等[33]發(fā)現(xiàn)重金屬元素Zn在銅尾（略）存在一定的富集,；馬杰等[34]發(fā)現(xiàn)鉛鋅尾礦庫的Cd元素污染最嚴重。根據(jù)走訪調(diào)查,，該尾礦庫在90年代曾發(fā)生過一次潰壩事故,，導致含重金屬的尾礦砂在壩下農(nóng)田土壤中歷史遺留。同時該尾礦庫壩下滲水通常與尾礦庫上游山塘水庫流水混合,，被用（略）域農(nóng)田的灌溉,，長期的滲水灌溉也會導致農(nóng)田土壤重金屬累積。因此,，可以推斷重金屬元素Cd,、Zn和Cu可能存在伴生關系，（略）農(nóng)田土壤中重金屬元素Cd,、Zn和Cu的主要來源是銅礦尾礦庫尾砂和壩下滲水,。
因子F2對農(nóng)田土壤重金屬貢獻率為5.19%，對Hg貢獻率最高,，達到61.5%,。Hg變異系數(shù)大于1，高度變異,，Igeo評價結(jié)果顯示47.36%的樣品處于輕度污染至偏重污染水平,，說明Hg受人為影響程度較大。相關性分析結(jié)果顯示,，Hg僅與Cr為顯著相關（P<0.05）,，而與其他重金屬無相關性。實地調(diào)查與人員訪談發(fā)現(xiàn),，Hg濃度較（略）域為種植有油菜,、葉用萵苣和包菜等綠葉類蔬菜的地塊，這些土地長期施用農(nóng)藥,、農(nóng)家肥和化肥,。一般認為農(nóng)業(yè)活動的污水灌溉、肥料、農(nóng)藥和抗菌劑施用等是土壤中Hg的重要來源[35],。魏迎輝等[26]發(fā)現(xiàn)鉛鋅礦周邊農(nóng)田Hg濃度（略）除了礦山附近外,，種植蔬（略）域也是Hg濃度（略）,。陳雪等[31]發(fā)現(xiàn)農(nóng)事活動對Hg富集影響大,，是Hg的主要來源,。以上結(jié)果均與本研究結(jié)果相似,。因此,，可以推斷（略）Hg的富集可能與局部肥料,、農(nóng)藥施用的農(nóng)業(yè)活動有關,。
因子F3對農(nóng)田土壤重金屬貢獻率為17.19%,，其中Cr,、Ni,、Pb和As在因子F3上具有較高載荷,，貢獻率分別為89.7%、82.7%,、75.0%和68.3%,。相關研究[31,36]表明Cr和Ni主要來源于成土母質(zhì)。Ni與Pb之間呈極顯著（P<0.01）相關,，Cr,、Ni和Pb濃度均值低于（略）土壤元素背景值,，變異系數(shù)小,，表明Cr,、Ni,、和Pb主要受地質(zhì)背景影響，Igeo評價結(jié)果也說明這一點,。另外,，本研究中As與Hg沒有相關性，這與通常認為的農(nóng)業(yè)活動（肥料農(nóng)藥施用）會導致As和Hg存在伴生污染[35]的結(jié)果不一致。但As與Pb具有極顯著相關性,，且二者具有中高度變異性，排除尾礦庫人為活動和農(nóng)業(yè)活動因素，分析Pb和As富集的原因可能是由于該尾礦庫正進行閉庫后生態(tài)修復的平整覆土等工作,，在主導風向為東風的氣象條件下,，（略）載土車的頻繁來往等交通運輸活動導致的尾礦砂土掉落、揚塵沉降等,。因此,，推斷（略）土壤Cr,、Ni,、Pb和As一定程度上受成土母質(zhì)的控制,，同（略）域交通運輸活動的影響,。
土壤重金屬源解析是土壤重金屬污染精準防控的基礎,，該方面研究逐漸集中在將數(shù)理統(tǒng)計分析模型與土壤重金屬的信息相結(jié)合，基于受體模型來分析重金屬的來源。與其他受體模型相比,，PMF模型對因子載荷和因子得分均做非負約束，且可以處理缺失,、錯誤的數(shù)據(jù),，因此通過每個可解釋的源成分譜,，定量識別污染物的來源。然而,，由于污染源排放復雜，疊加伴生現(xiàn)象等,，PMF模型難以有效分割共線性強的污染源,，識別出的因子往往為幾種污染源的混合,。因此，后續(xù)研究中,，應通過加強對潛在污染源的實際監(jiān)測,，并將其作為約束條件帶入模型計算，以提高結(jié)果的可靠性,，為重金屬污染有效防治提供參考依據(jù)。
3.結(jié)論
（1）（略）農(nóng)田土壤除Cr,、Ni,、Pb、As以外,，Cd,、Hg、Cu,、Zn濃度分別為土壤元素背景值的5.36,、2.06、8.19,、5.36倍,，具有高度變異性；單因子污染指數(shù)顯示,，As,、Pb、Cr,、Hg,、Ni單因子指數(shù)均小于1.0，處于無污染等級,，Cu,、Zn、Cd重度污染占比均達到10.5%,，中度污染占比為5.26%,；內(nèi)梅羅綜合指數(shù)評價顯示,，15.8%的點位（S17、S18,、S19）處于重度污染等級,，其他點位無污染，靠近尾礦庫的部分農(nóng)田土壤具有較高的污染等級,；地累積指數(shù)評價結(jié)果表明,，Cd、Cu,、Zn和Hg可能具有累積風險,。
（2）潛在生態(tài)風險評估結(jié)果表明，Cd為很強生態(tài)風險,，Hg為較強生態(tài)風險,，Cu為中等生態(tài)風險，其余重金屬均為輕微風險,，綜合潛在風險為較強風險（RI=308.91）,。生態(tài)風險預警評估結(jié)果表明，Cu為重警,，Cd和Zn為中警,，Hg為輕警，As為預警,，Pb,、Cr和Ni為無警，綜合生態(tài)風險預警為重警（IER=16.06）,。RI和IER空間分布基本一致,，主要受來Cu、Zn,、Cd和Hg的影響,。
（3）基于PMF解析出（略）土壤重金屬3個源，Cd,、Zn,、Cu主要受銅礦尾礦庫尾砂和壩下滲水的混合源影響，貢獻率分別為94.4%,、94.3%和67.1%,；Hg可能是受肥料農(nóng)藥施用等農(nóng)業(yè)活動源的影響，貢獻率為61.5%,；Cr,、Ni、Pb和As主要受成土母質(zhì)和交通運輸活動混合源的影響,，貢獻率分別為89.7%,、82.7%,、75.0%和68.3%。
向上滑動閱覽
參考文獻
[1]李海東,馬偉波,胡國長.礦區(qū)修復生態(tài)學理論與實踐[M].北京:中國環(huán)境出版集團,2022.
[2]LICF,WANGAJ,CHENXJ,etal.RegionaldistributionandsustainabledevelopmentstrategyofmineralresourcesinChina[J].ChineseGeographicalScience,2013,23(4):470-481.DOI:10.1007/s（略）11-z
[3]朱點鈺,楊倩琪.中（略）重金屬污染現(xiàn)狀及生態(tài)風險研究[J].礦產(chǎn)勘查,2018,9(4):747-750.
ZHUDY,YANGQQ.StudyonheavymetalcontaminationandecologicalriskinminingareasinChina[J].MineralExploration,2018,9(4):747-750.
[4]李傳飛,劉登璐,趙平,等.某區(qū)域（略）土壤重金屬污染與生態(tài)風險評價[J].四川環(huán)境,2021,40(2):141-148.
LICF,LIUDL,ZHAOP,etal.Evaluationofsoilheavymetalpollutionandecologicalriskinaminingarea[J].SichuanEnvironment,2021,40(2):141-148.
[5]吳洋,楊軍,周小勇,等.廣西（略）耕地土壤重金屬污染風險評價[J].環(huán)境科學,2015,36(8):（略）.
WUY,YANGJ,ZHOUXY,etal.RiskassessmentofheavymetalcontaminationinfarmlandsoilinDu'anautonomousCountyofGuangxiZhuangAutonomousregion,China[J].EnvironmentalScience,2015,36(8):（略）.
[6]戴彬,呂建樹,戰(zhàn)金成,等.山東省典型工（略）土壤重金屬來源,、空間分布及潛在生態(tài)風險評價[J].環(huán)境科學,2015,36(2):507-515.
DAIB,LüJS,ZHANJC,etal.Assessmentofsources,spatialdistributionandecologicalriskofheavymetalsinsoilsinatypicalindustry-basedcityofShandongProvince,EasternChina[J].EnvironmentalScience,2015,36(2):507-515.
[7]LüJS,LIUY,ZHANGZL,etal.IdentifyingtheoriginsandspatialdistributionsofheavymetalsinsoilsofJuCountry(EasternChina)usingmultivariateandgeostatisticalapproach[J].JournalofSoilsandSediments,2015,15(1):163-178.DOI:10.1007/s（略）37-x
[8]王恩瑞,楊光,昌盛,等.乾務水庫表層沉積物氮磷和重金屬時空分布特征與生態(tài)風險評價[J].環(huán)境工程技術學報,2023,13(3):（略）.DOI:10.12153/（略）-991X.（略）
WANGER,YANGG,CHANGS,etal.Occurrence,spatiotemporaldistributionandecologicalriskassessmentofnitrogen,phosphorusandheavymetalsinthesurfacesedimentsofQianwuReservoir[J].JournalofEnvironmentalEngineeringTechnology,2023,13(3):（略）.DOI:10.12153/（略）-991X.（略）
[9]鄧紅芮,李積普,羅琳,等.礦山公園土壤重金屬的生態(tài)與健康風險評價:以湖南M礦為例[J].環(huán)境工程技術學報,2023,13(2):760-768.DOI:10.12153/（略）-991X.（略）
DENGHR,LIJP,LUOL,etal.Ecologicalandhealthriskassessmentofsoilheavymetalsinaminingpark:acasestudyofMMineinHunanProvince[J].JournalofEnvironmentalEngineeringTechnology,2023,13(2):760-768.DOI:10.12153/（略）-991X.（略）
[10]韓娟娟,吳大鵬,張濤.北京某（略）土壤重金屬污染特征及風險評價[J].城市地質(zhì),2021,16(4):424-431.
HANJJ,WUDP,ZHANGT.PollutioncharacteristicsandriskassessmentofsoilheavymetalsinacertaingoldmineareainBeijing[J].UrbanGeology,2021,16(4):424-431.
[11]施建飛,靳正忠,周智彬,等.額爾齊斯河流域典型尾（略）周邊土壤重金屬污染評價[J].生態(tài)環(huán)境學報,2022,31(5):（略）.
SHIJF,JINZZ,ZHOUZB,etal.EvaluationofheavymetalpollutioninthesoilaroundatypicaltailingreservoirinIrtyshRiverBasin[J].EcologyandEnvironmentalSciences,2022,31(5):（略）.
[12]孫德堯,薛忠財,韓興,等.冀（略）（略）周邊耕地土壤重金屬污染特征及生態(tài)風險評價[J].生態(tài)與農(nóng)村環(huán)境學報,2020,36(2):242-249.
SUNDY,XUEZC,HANX,etal.PollutingcharacteristicsandecologicalriskassessmentofheavymetalsincultivatedlandaroundaminingareainnorthernHebeiProvince[J].JournalofEcologyandRuralEnvironment,2020,36(2):242-249.
[13]廖文靜,李波,陳杰,等.云南淘金溝錫礦山廢棄地重金屬來源分布及污染評價[J].地質(zhì)災害與環(huán)境保護,2022,33(1):107-113.
LIAOWJ,LIB,CHENJ,etal.SourcedistributionandpollutionevaluationofheavymetalsinthewastelandofTaojingoutinmine,Yunnan[J].JournalofGeologicalHazardsandEnvironmentPreservation,2022,33(1):107-113.
[14]趙家印,楊地,楊湘智,等.云南省某煤礦開采遺址周邊農(nóng)用地土壤重金屬污染評價及源解析研究[J].生態(tài)與農(nóng)村環(huán)境學報,2022,38(11):（略）.
ZHAOJY,YANGD,YANGXZ,etal.Pollutionassessmentandsourceidenti:cationofheavymetalsinfarmlandsoilsaroundacoalmineareainYunnanProvince[J].JournalofEcologyandRuralEnvironment,2022,38(11):（略）.
[15]LüJS.Multivariatereceptormodelsandrobustgeostatisticstoestimatesourceapportionmentofheavymetalsinsoils[J].EnvironmentalPollution,2019,244:72-83.DOI:10.1016/（略）.09.147
[16]ZHANGFG,WANGCW,CHENGXM,etal.Ecologicalassessment,spatialanalysis,andpotentialsourcesofheavymetals(HMs)insoilswithhighbackgroundvaluesinthelead-zincmine,HezhangCounty,southwesternChina[J].Water,2022,14(5):783-783.DOI:10.3390/w（略）
[17]朱曉麗,薛博倩,李雪,等.基于PMF模型的寶雞鉛鋅尾礦庫周邊農(nóng)田土壤重金屬源解析[J].西北大學學報(自然科學版),2021,51(1):43-53.
ZHUXL,XUEBQ,LIX,etal.Sourcesapportionmentofheavymetalsinfarmlandsoilaroundlead-zinctailingsreservoirbasedonPMFmodel[J].JournalofNorthwestUniversity(NaturalScienceEdition),2021,51(1):43-53.
[18]范允慧,王艷青.浙江省四大（略）土壤元素背景值特征[J].物探與化探,2009,33(2):132-134.
FANYH,WANGYQ.BackgroundcharacteristicsofsoilelementsinfourplainsofZhejiangProvince[J].GeophysicalandGeochemicalExploration,2009,33(2):132-134.
[19]李軍,李開明,王曉槐,等.（略）黃河風情線地表積塵及周邊綠地土壤重金屬污染特征及風險評價[J/OL].環(huán)境科學,2022.[（略）].https://（略）.（略）
LIJ,LIKM,WANGXH,etal.PollutioncharacteristicsandriskassessmentofheavymetalsinsurfacedustsanditssurroundinggreenlandsoilsfromYellowRivercustomtouristlineinLanzhou[J/OL].EnvironmentalScience,2022.[（略）].https://（略）.（略）
[20]劉洋,劉明慶,王磊,等.云南某廢棄硅廠周邊農(nóng)田土壤重金屬污染評價[J].農(nóng)業(yè)環(huán)境科學學報,2022,41(4):785-793.
LIUY,LIUMQ,WANGL,etal.EvaluationofheavymetalpollutioninfarmlandsoilaroundanabandonedsiliconplantinYunnan[J].JournalofAgro-EnvironmentScience,2022,41(4):785-793.
[21]國家環(huán)境保護總局.土壤環(huán)境監(jiān)測技術規(guī)范:HJ/T166—2004[S].北京:中國環(huán)境出版社,2004.
[22]徐爭啟,倪師軍,庹先國,等.潛在生態(tài)危害指數(shù)法評價中重金屬毒性系數(shù)計算[J].環(huán)境科學與技術,2008,31(2):112-115.
XUZQ,NISJ,TUOXG,etal.Calculationofheavymetals'toxicitycoefficientintheevaluationofpotentialecologicalriskindex[J].EnvironmentalScience&Technology,2008,31(2):112-115.
[23]YUANGL,SUNTH,HANP,etal.Sourceidentificationandecologicalriskassessmentofheavymetalsintopsoilusingenvironmentalgeochemicalmapping:typicalurbanrenewalareainBeijing,China[J].JournalofGeochemicalExploration,2014,136:40-47.DOI:10.1016/（略）.10.002
[24]李嬌,滕彥國,吳勁,等.基于PMF模型及地統(tǒng)計法的樂安河中上（略）土壤重金屬來源解析[J].環(huán)境科學研究,2019,32(6):984-992.
LIJ,TENGYG,WUJ,etal.SourceapportionmentofsoilheavymetalinthemiddleandupperreachesofLe'anRiverbasedonPMFmodelandgeostatistics[J].ResearchofEnvironmentalSciences,2019,32(6):984-992.
[25]陳林,馬琨,馬建軍,等.寧夏引（略）農(nóng)田土壤重金屬生態(tài)風險評價及來源解析[J/OL].環(huán)境科學,2022:1-19.[（略）].https://（略）.（略）
CHENL,MAK,MAJJ,etal.RiskassessmentandsourcesofheavymetalsinfarmlandsoilsofYellowRiverirrigationareaofNingxia[J/OL].EnvironmentalScience,2022:1-19.[（略）].https://（略）.（略）
[26]魏迎輝,李國琛,王顏紅,等.PMF模型的影響因素考察:以某鉛鋅礦周邊農(nóng)田土壤重金屬源解析為例[J].農(nóng)業(yè)環(huán)境科學學報,2018,37(11):（略）.
WEIYH,LIGC,WANGYH,etal.InvestigatingfactorsinfluencingthePMFmodel:acasestudyofsourceapportionmentofheavymetalsinfarmlandsoilsnearalead-zincore[J].JournalofAgro-EnvironmentScience,2018,37(11):（略）.
[27]WILDINGL.Spatialvariability:itsdocumentation,accommodationandimplicationtosoilsurveys[C]//Soilspatialvariability.LasVegas,1985:166-194.
[28]馬宏宏,余濤,楊忠芳,等.（略）土壤重金屬空間插值方法與污染評價[J].環(huán)境科學,2018,39(10):（略）.
MAHH,YUT,YANGZF,etal.Spatialinterpolationmethodsandpollutionassessmentofheavymetalsofsoilintypicalareas[J].EnvironmentalScience,2018,39(10):（略）.
[29]杜貫新,閆百泉,孫雨,等.松嫩平原（略）西北部阿榮旗黑土重金屬分布特征,、評價及預警[J].地質(zhì)科學,2022,57(2):606-621.
DUGX,YANBQ,SUNY,etal.Distributioncharacteristics,evaluationandearlywarningofheavymetalsinArongqiblacksoilinnorthwestofSongnenblacksoilregion[J].ChineseJournalofGeology(ScientiaGeologicaSinica),2022,57(2):606-621.
[30]米曉軍,任雯,雒瓊,等.新疆準噶爾盆地未開墾鹽堿地土壤重金屬評價及其來源[J].（略）研究,2019,36(4):824-834.
MIXJ,RENW,LUOQ,etal.Evaluationandtheirsourcesofheavymetalsinuncultivatedsaline-alkalinesoilintheJunggarBasin,Xinjiang[J].AridZoneResearch,2019,36(4):824-834.
[31]陳雪,汪小祥,劉敬青.小流域（略）土壤重金屬污染評價與來源綜合解析[J].有色金屬(礦山部分),2022,74(6):74-81.
CHENX,WANGXX,LIUJQ.Assessmentandsourcescomprehensiveidentificationofsoilheavymetalpollutioninsmall-scaledrainagecatchmentofminingareas[J].NonferrousMetals(MiningSection),2022,74(6):74-81.
[32]陳航,王穎,王澍.銅（略）周邊農(nóng)田土壤重金屬來源解析及污染評價[J].環(huán)境科學,2022,43(5):（略）.
CHENH,WANGY,WANGS.SourceanalysisandpollutionassessmentofheavymetalsinfarmlandsoilaroundTongshanminingarea[J].EnvironmentalScience,2022,43(5):（略）.
[33]杜立宇,梁成華,劉桂琴.紅透山銅尾礦重金屬分布及其對土壤重金屬污染的影響[J].土壤通報,2008,39(4):938-941.DOI:10.19336/（略）.（略）
DULY,LIANGCH,LIUGQ.ThedistributioncharacteristicsofheavymtealsinCuminetailingandeffectofheavymetalsonpollutionofsoilsintheareasofHongToumountians[J].ChineseJournalofSoilScience,2008,39(4):938-941.DOI:10.19336/（略）.（略）
[34]馬杰,張秀,劉今朝,等.某鉛鋅尾礦庫周邊土壤重金屬污染特征及其來源分析[J].有色金屬(冶煉部分),2022(6):101-109.
MAJ,ZHANGX,LIUJZ,etal.Pollutioncharacteristicsandsourceanalysisofheavymetalsinsoilsofalead-zinctailingspond[J].NonferrousMetals(ExtractiveMetallurgy),2022(6):101-109.
[35]孫晶,李偉,呂學斌.土壤中汞和砷的環(huán)境效應、來源及影響因素[J].綠色科技,2020(14):148-151.
SUNJ,LIW,LüXB.Environmentaleffect,sourcesandinfluencingfactorsofmercuryandarsenicinsoil[J].JournalofGreenScienceandTechnology,2020(14):148-151.
[36]黃華斌,林承奇,胡恭任,等.基于PMF模型的九龍江流域農(nóng)田土壤重金屬來源解析[J].環(huán)境科學,2020,41(1):430-437.
HUANGHB,LINCQ,HUGR,etal.SourceappointmentofheavymetalsinagriculturalsoilsoftheJiulongRiverBasinbasedonpositivematrixfactorization[J].EnvironmentalScience,2020,41(1):430-437.□
引用本文:呂玉娟,，王秋月,，孫雪梅，等.浙江省某尾礦庫周邊農(nóng)田土壤重金屬污染特征及來源解析[J].環(huán)境工程技術學報,，2023,，13（4）：（略）.DOI:10.12153/（略）-991X.（略）文章來源：環(huán)境工程技術學報
免責聲明：（略）絡轉(zhuǎn)載，版權歸原作者和原出處所有,。部分文章因轉(zhuǎn)載眾多,，僅標明轉(zhuǎn)載來源，若涉及版權問題,，請聯(lián)系我們,。