煤礦周邊土壤重金屬空間分布研究—以H煤礦為例摘要本研究分析了煤礦周邊土壤重金屬（Cr、Ni、Cu、Zn、As、Cd、Pb）的含量特征，空間分布特征，垂向遷移規(guī)律和來源。結果表明，研究區(qū)土壤7種重金屬中，Cr、Pb、Cu、As、Cd的平均含量均高于安徽省表層土壤背景值，但土壤生態(tài)環(huán)境的風險低；研究區(qū)域內(nèi)，Cr的高值區(qū)域出現(xiàn)在北部，Ni的高值區(qū)域出現(xiàn)在南部煤矸石山附近和東北部，Cu和Zn的高值區(qū)域出現(xiàn)在東南部煤矸石山和燃煤電廠之間，As和Cd的高值區(qū)域分別出現(xiàn)在西北部和東部農(nóng)田區(qū)；土壤重金屬垂向遷移難易程度由易到難依次為：Cr＞Cd＞Zn、Ni＞Cu&Zn＞Pb；Cu、Zn和Pb與汽車尾氣排放及交通運輸揚塵有關，Cd和Ni主要是自然地質作用以及煤矸石山淋濾、工業(yè)污染物的排放所致，As主要由農(nóng)業(yè)活動引起，Cr則與煤炭燃燒關系密切。關鍵詞：重金屬；克里金插值法；垂向遷移；多元統(tǒng)計分析目錄TOC\o"1-2"\h\u4250引言 198871材料與方法 38081.1研究區(qū)概況 3304541.2樣品采集與處理 440321.3樣品測試 4305812結果與討論 5158832.1表層土壤重金屬含量特征 5300462.2表層土壤重金屬的空間分布特征 6249772.3土壤重金屬垂向遷移特征 835702.4相關性分析 9223262.5聚類分析和主成分分析 10219173結論 1428553參考文獻 15引言隨著煤炭資源的開發(fā)與利用，會產(chǎn)生大量的廢棄物，如煤矸石、礦井水、粉煤灰等，對礦區(qū)周邊環(huán)境造成污染[1-3]。特別是露天堆放的煤矸石和粉煤灰，在自燃、雨水淋濾、風力吹揚等作用下，其中的微量重金屬元素向周邊環(huán)境介質擴散、遷移，污染土壤，導致土壤重金屬含量增加[4-6]。煤礦電廠在燃煤過程中更會釋放出有毒有害物質，特別是重金屬元素，對礦區(qū)環(huán)境也會產(chǎn)生嚴重影響。土壤是地球上一切生命賴以生存和發(fā)育的最基本的資源，重金屬在土壤中隱蔽性強、潛伏周期長、毒副作用強烈、降解能力低[7]，土壤一旦受到污染，將很難修復。因此，對于煤礦礦業(yè)工程影響區(qū)土壤重金屬污染的研究顯得尤為重要。袁新田[8]等對宿州市煤礦區(qū)周邊農(nóng)田表層土壤中的Cu、Zn、Cr、Cd、Pb、Hg、As7種重金屬污染特征進行研究，結果表明：煤礦區(qū)周邊500m內(nèi)農(nóng)田表層土壤除Cu、Zn、Pb外，Cr、Cd、Hg、As均超出國家土壤環(huán)境質量標準(I級)。袁建梅[9]等以重慶市的6個煤礦區(qū)為研究對象，測定了重金屬Cr、Ni、Cu、Zn、Cd、Pb、As、Hg含量，分析結果表明：研究區(qū)部分土樣Cr、Ni、Cu、Zn、Cd超標外，其余各土樣均滿足國家土壤環(huán)境質量相關標準，綜合污染程度為清潔，潛在生態(tài)風險為輕度或中度。熊佳[10]等對貴州省獨山縣某銻礦冶煉廠周邊土壤重金屬進行了空間分布特征分析，得出新廠區(qū)附近300m范圍內(nèi)土壤中Sb含量較高，呈現(xiàn)出隨距離增加而降低的分布特征。樊志穎[11]等分析了色季拉山森林土壤重金屬含量的空間分布特征，得出重金屬含量在空間分布上受坡向海拔與土壤深度的影響，影響最大的是坡向。蘇海民[12]等采用Hakanson潛在生態(tài)風險指數(shù)法，研究了宿州市礦區(qū)周邊土壤重金屬污染潛在的生態(tài)危害性，認為除Cd和Hg外，重金屬潛在生態(tài)風險系數(shù)均屬于輕度風險危害。AblizAbdugheni[13]等研究了西北采煤區(qū)土壤重金屬的來源，Zn、Cu來源于母體物質，Cr、As、Hg來源于燃煤、化工、交通等人類活動，Pb受自然因素和人類活動的雙重影響。煤礦的開采歷史、開發(fā)規(guī)模以及煤炭資源的利用方式不同，對礦區(qū)周邊的環(huán)境影響截然不同。本文選擇宿南礦區(qū)開采歷史最長、開發(fā)規(guī)模最大且有坑口電廠的蘆嶺煤礦為研究對象，分析研究周邊土壤中的Cr、Ni、Cu、Zn、As、Cd和Pb的含量特征、空間分布特征及遷移規(guī)律，探討重金屬污染可能的來源，以期為煤礦周邊土壤環(huán)境保護和修復治理提供科學依據(jù)。1材料與方法1.1研究區(qū)概況蘆嶺煤礦地處淮北煤田宿南礦區(qū)，面積23km2，1969年12月建成投產(chǎn)，有50多年的開采歷史。該礦煤炭產(chǎn)量220萬噸/a。2015年7月，該礦電廠裝機30000kW，日均發(fā)電46000kWh。多年的煤炭開采與利用，礦區(qū)表生環(huán)境發(fā)生了巨大變化，特別是煤矸石山、塌陷湖等構成了礦區(qū)周邊土壤污染的潛在威脅。研究區(qū)發(fā)育的沱河為季節(jié)性河流；地勢平坦，海拔在+24m左右；根據(jù)CST土壤分類，判斷該研究區(qū)土壤主要為砂姜黑土；主要農(nóng)作物為玉米、小麥；夏季主導風向為東南風，冬季主導風向為東北風。研究區(qū)地理位置見圖1。圖1研究區(qū)地理位置及采樣點分布圖1.2樣品采集與處理2019年7月，在考慮了研究區(qū)的主導風向和污染源后，以煤矸石山和坑口電廠為中心，沿東南順風方向，采用扇形布點法，采集19個表層土壤樣品，另外在煤矸石堆和坑口電廠之間采集2個表層土壤樣品，共21個表層土壤樣品（B1~B21）；5個垂直分層土壤樣品（P2~P5），垂直分層采樣以10cm的深度為間隔，分為六層，共60cm。采樣點分布見圖1。在20cm×20cm方形范圍內(nèi)取土壤樣品，采集深度0～10cm，樣品裝入干凈的密封袋中，貼好標簽，并用GPS對采樣點進行定位，同時記錄采樣點周邊環(huán)境情況。樣品運回實驗室后，經(jīng)過自然風干—研磨—過200目篩—裝入密封袋并貼好標簽，備測。1.3樣品測試（1）壓片用熒光專用硼酸磨具和手動粉末壓片機壓片。取適量過篩之后的土壤樣品放到漏斗中間，用壓桿將樣品壓平；再將硼酸粉末放到漏斗外側，之后迅速將漏斗抽出，依次放入墊片和上壓頭；將組裝好的模具放到壓片機中，加壓到10Mpa左右，保壓1分鐘；將模具倒置，去掉模具底，然后將模具放到壓片機中，用絲杠頂出已壓好的樣品。（2）測試用X射線熒光光譜儀對樣品的Cr、Ni、Cu、Zn、As、Cd和Pb元素進行含量測定，所用標準物為GBW07430（GSS-16）。2結果與討論2.1表層土壤重金屬含量特征將研究區(qū)表層土壤重金屬的含量測試結果，與安徽省表層土壤重金屬含量背景值、國家土壤環(huán)境質量土壤污染風險篩選值[14]和歐洲標準指南土壤重金屬最大允許限值[15]進行對比，其統(tǒng)計特征如表1所示。表1表層土壤重金屬含統(tǒng)計特征分析元素平均值/(mg·kg-1)標準差/(mg·kg-1)最大值/(mg·kg-1)最小值/(mg·kg-1)背景值/(mg·kg-1)中國安全限值/(mg·kg-1)歐洲安全限值/(mg·kg-1)加拿大安全限值/(%)變異系數(shù)/(%)Cr77.326.3288.7062.8066.5200100648Zn62.779.7493.9852.976225030025016Ni21.578.2835.5411.0529.8100504538Pb36.893.4049.2233.8026.6120100709Cu24.237.0039.3412.0620.41001006329As17.973.6925.2111.20930201221Cd0.220.020.270.160.0970.3031.411可見，除Ni的均值含量低于背景值、Zn的均值含量基本接近背景值外，5種重金屬Cr、Pb、Cu、As、Cd的均值含量分別為77.32、36.89、24.23、17.97、0.22mg·kg-1，分別是安徽省表層土壤背景值的1.16、1.39、1.19、2.00、2.27倍（圖2），說明研究區(qū)內(nèi)土壤中這5種重金屬已經(jīng)產(chǎn)生了不同程度的污染累積，其中Cr、Cu含量相對較低，Cd含量相對較高。與國家土壤環(huán)境質量土壤污染風險篩選值和歐洲標準指南土壤重金屬最大允許限值相比較，研究區(qū)表層土壤中Cr、Zn、Ni、Pb、Cu、As、Cd7種重金屬的均值含量均低于安全限值，說明該區(qū)域土壤生態(tài)環(huán)境的風險低，一般情況下可以忽略[14]。圖2研究區(qū)重金屬含量與安徽省表層土壤背景值柱狀圖變異系數(shù)為標準差與平均值之比，可以反映出數(shù)據(jù)的離散程度，從而間接反映出外源因素對重金屬含量的影響[16]。一般認為，CV＜10%為弱變異，10%～100%為中等強度變異，CV≥100%為強變異[17]?？芍狽i、Zn、Cu、As、Cd（38％、16％、29％、21％、11％）屬于中等強度變異，Cr、Pb（8％、9%）屬于弱變異，其中Ni、Cu的變異系數(shù)最高，離散程度最大，可能存在外源因素的干擾。2.2表層土壤重金屬的空間分布特征利用克里金插值分析方法，所得出的重金屬的空間分布如圖3所示。Cr（圖3a）的高值區(qū)域出現(xiàn)在北部，此處多為居民區(qū)，可能與居民生活污水排放有關[18]。重金屬的擴散受到盛行風向的影響，含量隨距離的增加而逐漸下降[19]。以燃煤電廠為起點，在研究區(qū)主導風向（夏季東南風）影響區(qū)域內(nèi)，土壤中Cr含量呈現(xiàn)出隨距離增加逐漸降低的趨勢。Ni（圖3b）的高值區(qū)域出現(xiàn)在南部煤矸石山附近，位于下風向處的Ni含量要明顯高于上風向，且在順風方向上Ni的含量隨著距離的增加有逐漸降低的趨勢。另外，Ni的高值區(qū)域還出現(xiàn)在東北部，可能與此處居民區(qū)[18]和農(nóng)田的分布有關，受到人為攪動較大。Cu（圖3c）和Zn（圖3d）的高值區(qū)域出現(xiàn)在東南部煤矸石山和燃煤電廠之間，此處分布有煤礦以及礦區(qū)主干道路，過往車輛頻繁，可能與道路揚塵以及汽車尾氣排放有關[20]；Cu含量在一定程度上呈現(xiàn)出從東南向西北逐漸遞減的趨勢，而Zn含量的分布除高值區(qū)域外，整體分布比較均勻。As（圖3e）的高值區(qū)域出現(xiàn)在西北部農(nóng)田區(qū)，Cd（圖3f）的高值區(qū)域出現(xiàn)在東部農(nóng)田區(qū)，可能與農(nóng)業(yè)活動有關[21]；Cd的含量整體呈現(xiàn)出由東向西逐漸遞減的趨勢。圖3研究區(qū)土壤重金屬空間分布圖2.3土壤重金屬垂向遷移特征為了研究礦區(qū)土壤重金屬的垂向遷移特征，對研究區(qū)土壤進行垂直分層采樣，分析結果見圖4、5。由淺入深，Cr、Zn、Ni、Pb、Cu、As和Cd的含量范圍依次為61.15～79.71、58.89～69.31、11.21～39.84、36.80～37.75、19.74～32.73、10.35～15.69和0.17～0.23mg·kg-1，6個垂直分層樣的7種重金屬元素均已達到中國土壤環(huán)境質量土壤污染風險篩選值和歐洲標準指南土壤重金屬最大允許限值。圖4土壤重金屬剖面分布特征柱狀圖圖5土壤重金屬剖面分布特征散點圖由圖4和圖5可知，除Cd和Pb的含量在0~60cm深度范圍內(nèi)基本保持不變外，其余5種重金屬元素（Cr、Ni、Cu、Zn、As）的含量在0~30cm范圍內(nèi)都具有較為明顯的波動，說明這5種重金屬元素在0~30cm范圍內(nèi)存在較為明顯的遷移。從0至30cm,，As含量的變化趨勢為先降低后升高，最高含量出現(xiàn)在20~30cm處，可能與天然降水和引水灌溉有關；Ni和Cu含量的變化趨勢是先升高后降低，且都是在10~20cm處急劇升高，在20~30cm處急劇降低，濃度最大值均出現(xiàn)在10~20cm處，最低值均出現(xiàn)在20~30cm處，可能是因為10~20cm處屬于施肥深度和植物根系發(fā)育深度，有機肥的施用和植物根系的發(fā)育均會導致重金屬難以向下遷移[22]，從而出現(xiàn)重金屬濃度在10～20cm急劇升高的現(xiàn)象；Cr和Zn的含量在0~30cm范圍內(nèi)呈現(xiàn)逐漸降低的趨勢，且最高含量都出現(xiàn)在表層土壤，說明這2種元素受到人為因素干擾較大。從30至60cm，Ni、Cu、Cr的含量都表現(xiàn)出逐漸增加的趨勢，這可能與淀積作用有關[23]，天然降雨或灌溉引起的淋溶作用使重金屬元素向下遷移，在30~60cm范圍內(nèi)發(fā)生淀積，導致重金屬元素含量升高；As和Zn的含量呈現(xiàn)逐漸降低的趨勢，說明這兩種元素也具有垂向遷移的能力，但遷移量較小。由以上分析可得，研究區(qū)土壤7種重金屬元素的垂向遷移特征有所差異，其難易程度由易到難依次為Cr＞Ni&Cu＞As＞Zn＞Cd&Pb(&連接的2種元素代表遷移能力相近)。造成這種差異的原因與土壤中有機質和石灰的含量、黏土礦物的種類和數(shù)量、重金屬對某些負電荷表面的吸附能力、植物根系對重金屬的吸附能力以及土壤中水的蒸發(fā)等因素皆有關[24-25]。2.4相關性分析研究表明，重金屬元素之間的相關性是反應各元素是否具有相同來源的重要指標[26]。一般認為元素的相關性系數(shù)大于0.3時[27]，表現(xiàn)為較強相關性，說明元素可能具有相同的來源。重金屬含量間的Pearson相關性分析結果見表2，Cd—Ni、Zn—Cu、Pb—Cu、Pb—Zn、Pb—Cd的Pearson相關系數(shù)的絕對值分別為0.563、0.469、0.518、0.375、0.397，說明Cd和Ni，Cd和Pb，Zn、Cu、Pb兩兩之間存在較強的相關性，可能具有相同的來源。其它元素之間的相關性系數(shù)均小于0.3，相關性較弱，這些元素是否來自同一來源可以進一步通過聚類分析和主成分分析進行判斷。表2重金屬含量相關性分析結果元素CrNiCuZnAsCdPbCr1.000Ni0.0321.000Cu0.0170.0991.000Zn0.186-0.1040.469*1.000As-0.086-0.1690.072-0.1931.000Cd-0.1120.563**-0.007-0.168-0.0301.000Pb-0.142-0.1450.518*0.375-0.294-0.3971.000**相關性在0.01水平（雙側）上顯著*相關性在0.05（雙側）水平上顯著2.5聚類分析和主成分分析聚類分析和主成分分析也可以反映出各元素之間是否具有相同的來源[26]。采樣點的聚類分析結果見圖6，當組間距離標準為10時，21個采樣點可以劃分為差異明顯的4個組，但86％的采樣點集中于前兩個組，說明采樣點主要受到了兩個污染源的影響，而第一個組的采樣點大多數(shù)以煤矸石山為中心而分布，第二個組的采樣點大多以電廠為中心而分布，說明煤矸石山和電廠是兩個不同的污染源。接下來采用主成分分析法對土壤中7種重金屬元素的來源做進一步判斷。圖6重金屬元素含量聚類分析結果主成分分析法利用降維的思想，把多個評價指標轉化為少數(shù)幾個具有代表性的綜合性指標[28-29]。研究區(qū)土壤重金屬主成分因子的載荷量如表3和圖7所示，7種重金屬可由4個公因子反映其83.662％的信息量，丟失信息較少，而且旋轉前后總的方差累積貢獻率相同。表3土壤重金屬主成分分析結果元素旋轉前旋轉后PC1PC2PC3PC4PC1PC2PC3PC4Cr0.0890.050-0.8470.4340.039-0.032-0.0440.955Ni-0.3710.805-0.0010.0010.0460.871-0.1580.029Cu0.6310.4500.3330.4140.9120.1320.174-0.050Zn0.7110.278-0.1900.2300.730-0.124-0.1630.312As-0.232-0.4030.4470.732-0.038-0.0940.968-0.063Cd-0.5870.6600.1470.090-0.1320.8830.078-0.071Pb0.8270.1280.229-0.2900.693-0.332-0.390-0.309特征值2.1321.5441.1381.0411.8681.6911.1791.118方差貢獻率/%30.46122.06216.26314.87426.68724.16016.84915.965總方差貢獻率/%30.46152.52568.78783.66226.68750.84867.69683.662圖7土壤重金屬因子載荷三維散點圖從旋轉后因子的載荷量可知，主成分1（PC1）的方差貢獻率為26.687％，Cu、Zn、Pb都具有較大的正載荷，則Cu、Zn、Pb可能來自同一污染源，這與上述相關性分析結果一致。由表1可知，Cu的變異系數(shù)較大，且少數(shù)樣品含量較高，Pb雖屬于低等變異，但其平均含量為安徽省土壤背景值的1.39倍，已經(jīng)產(chǎn)生了污染累積。在大多數(shù)情況下，人為源對微量元素輸入的貢獻大于自然源[30]，因此Cu、Zn和Pb可以認為受人為因子的控制。由于采樣點位于煤矸石山與電廠附近農(nóng)田區(qū)，以及礦區(qū)主干道路沿線，過往車輛頻繁，汽車尾氣排放嚴重，而Cu、Pb為機動車尾氣的標識元素[31]，Zn主要來源于輪胎和剎車片的磨損[32]。因此，PC1可解析為汽車尾氣排放和交通運輸揚塵。主成分2（PC2）的貢獻率占總貢獻率的24.160％，Cd和Ni的載荷較大，與上述相關性分析結果一致。本研究區(qū)中Cd的累積程度和生態(tài)危害最大，原因是長期堆放的煤矸石山經(jīng)過雨水淋濾后，重金屬Cd隨著流水遷移擴散到周邊土壤中，形成累積。土壤中Ni主要來源于土壤、巖石的形成過程以及巖石的風化等[33]。但研究區(qū)Ni的變異系數(shù)較大，說明還受到人為因素的干擾。Ni被廣泛應用于生產(chǎn)工業(yè)機械和精密電子儀器、冶金和電鍍等領域，Ni的氧化物和氫氧化物還可用于充電電池[33]。距研究區(qū)正南方向2500m處有一電子廠，由于該地區(qū)夏季盛行東南風，會對該研究區(qū)土壤中的重金屬Ni造成累積。因此，PC2可以解析為地質作用（即自然資源）和煤矸石山淋濾，其次是工業(yè)污染物的排放。主成分3（PC3）解釋了總方差貢獻率的16.849％，As的載荷較大。因為采樣點位于農(nóng)田區(qū)，種植作物需大量使用化肥，化肥中存在的有害重金屬元素（如As、Cd、Pb等）可能會以化肥的形態(tài)而保存在土壤里，造成As的大量積累[34]。因此，PC3代表農(nóng)業(yè)活動。主成分4（PC4）的貢獻率占總貢獻率的15.965％，Cr的載荷較大。Cr在燃煤電廠運行過程中可通過污水排放、大氣沉降等方式進入周邊土壤，造成一定的累積。因此，PC4可解釋為煤炭燃燒?？傮w來說，土壤污染主要有三個途徑，一是自然來源，二是人為來源，三是自然和人為來源的綜合。人為來源主要有交通污染、工農(nóng)業(yè)活動、采礦活動等。AbdugheniAbliz[13]等對西北采煤區(qū)土壤中重金屬來源進行解析，表明Zn和Cu主要來源于稀土元素，Cr、As、Hg主要來源于人類活動，如燃煤、化工、交通，Pb的產(chǎn)生是自然因素和人類活動共同作用的結果。ShuangXiong[5]等對臨渙礦區(qū)土壤重金屬的研究表明，Pb、Cr、Cd的污染來源于農(nóng)藥和化肥，F(xiàn)e和Ni來源于母體物質和工業(yè)活動，Cu和Mn受采煤活動影響，Zn主要來源于交通污染。KuangjiaLi[35]等對河南省某煤礦周邊土壤重金屬進行源解析，結論是，Cd、Pb、Cu、Zn元素來源于煤矸石堆。劉軍[36]等通過對草原區(qū)煤礦周邊旱作農(nóng)田土壤重金屬的研究發(fā)現(xiàn)，Zn主要來源于交通運輸，Ni主要來源于農(nóng)業(yè)生產(chǎn)活動，Cd與Hg來源于自然成土母質與人為活動。根據(jù)前人的研究可發(fā)現(xiàn)，Cu、Zn、Pb主要來源于交通污染，Cd、Cr主要來源于農(nóng)業(yè)活動，Ni、As主要與成土母質有關。本研究中，Cu、Zn、Pb的來源與前人研究基本一致，而對Cd、Cr、Ni、As的來源解析與前人研究略有差異，主要原因是蘆嶺煤礦開采歷史悠久，環(huán)境累積效應較強。此外本研究區(qū)不僅有采礦活動，還有電廠、電子企業(yè)等，工礦活動及周邊環(huán)境復雜，土壤污染源的判斷具有多解性。3結論通過上述研究，得出以下主要結論：（1）研究區(qū)土壤7種重金屬，除Ni和Zn外，其余5種重金屬（Cr、Pb、Cu、As、Cd）的平均含量均高于安徽省表層土壤背景值，但7種重金屬的平均含量均低于國家土壤環(huán)境質量土壤污染風險篩選值，土壤生態(tài)環(huán)境的風險低。（2）研究區(qū)域內(nèi)，Cr的高值區(qū)域出現(xiàn)在北部，且以燃煤電廠為起點，呈現(xiàn)出隨距離增加而降低的分布特征；Ni的高值區(qū)域出現(xiàn)在南部煤矸石山附近和東北部，受到風向的影響明顯；Cu和Zn的高值區(qū)域出現(xiàn)在東南部煤矸石山和燃煤電廠之間，Cu含量在一定程度上呈現(xiàn)出從東南向西北逐漸遞減的趨勢，而Zn含量的分布除高值區(qū)域外，整體分布比較均勻；As的高值區(qū)域出現(xiàn)在西北部農(nóng)田區(qū)，整體趨勢為由東向西逐漸遞減；Cd的高值區(qū)域出現(xiàn)在東部農(nóng)田區(qū)。（3）土壤重金屬剖面特征表明，該研究區(qū)土壤中的重金屬含量并沒有隨著深度增加而逐漸降低，0~30cm范圍內(nèi)5種重金屬（Cr、Ni、Cu、Zn、As）含量出現(xiàn)較為明顯的波動，說明這5種重金屬元素在0~30cm范圍內(nèi)存在較為明顯的遷移。研究區(qū)土壤重金屬垂向遷移難易程度由易到難依次為Cr＞Ni&Cu＞As＞Zn＞Cd&Pb，存在差異的原因與土壤中黏土礦物的種類和數(shù)量、重金屬對某些負電荷表面的吸附能力等皆有關。（4）Pearson相關性分析表明，研究區(qū)重金屬Cd和Ni，Cd和Pb，Zn、Cu和Pb兩兩之間存在較強的相關性，體現(xiàn)煤矸石山和電廠兩個主要的污染源特征。主成分分析結果顯示，Cu、Zn和Pb的來源主要為汽車尾氣排放及交通運輸揚塵，Cd和Ni主要來源是自然地質作用以及煤矸石山淋濾、工業(yè)污染物的排放，As主要是農(nóng)業(yè)活動引起的，Cr則與煤炭燃燒關系密切。參考文獻HuiliQiu,HerongGui,LinCui,etal.Hydrogeochemicalcharacteristicsandwaterqualityassessmentofshallowgroundwater:acasestudyfromLinhuancoal-miningdistrictinnorthernAnhuiProvince,China[J].WaterScience&Technology-WaterSupply,2019,19(5):1572-1578.ManliLin,WeihuaPeng,HerongGui.Heavymetalsindeepgroundwaterwithincoalminingarea,northernAnhuiprovince,China:concentration,relationship,andsourceapportionment[J]ADDINNE.Rep.ArabianJournalofGeosciences,2016,9(4):319.SiyanZeng,JingMa,YongjunYang.SpatialassessmentoffarmlandsoilpollutionanditspotentialhumanhealthrisksinChina[J].ScienceoftheTotalEnvironment,2019,687:642-653.DaweiHuang,HerongGui.Distributionfeaturesandinternalrelationsofheavymetalsinsoil–maizesystemofminingarea,AnhuiProvince,EasternChina[J].HumanandEcologicalRiskAssessment:AnInternationalJournal,2019,25

(4):863-881.ShuangXiong,HerongGui,ManliLin,etal.Contentsandpollutioncharacteristicsofheavymetalsinsoilfromcoalminingarea:acasestudyinLinhuanminingdistrict,northernAnhuiProvince[J].FreseniusEnvironmentalBulletin,2017,26(3):1989-1997.XianghongLiu,LinhuaSun,SongChen.Heavymetalpollutionofsoilaroundtheganguehill:acasestudyfromZhuxiangzhuangcoalmine,northernAnhuiprovince,China[J].AdvancedMaterialsResearch,2012,356-360:114-118.李靜,常勇,潘淑穎.土壤重金屬污染評價方法的研究[J].農(nóng)業(yè)災害研究,2012,2(4):50-52.袁新田,張春麗,孫倩,等.宿州市煤礦區(qū)農(nóng)田土壤重金屬含量特征[J].環(huán)境化學,2011,30(8):1451-1455.袁建梅,楊德敏,郭呈宇,等.重慶市煤礦區(qū)土壤重金屬污染評價研析[J].環(huán)境影響評價,2019,41(2):58-62.熊佳,韓志偉,吳攀,等.獨山銻冶煉廠周邊土壤銻砷空間分布特征、污染評價及健康風險評估[J].環(huán)境科學學報,2020,40(2):655-664.樊志穎,李江榮,高郯,等.色季拉山森林土壤重金屬空間分布特征及污染評價[J].西北農(nóng)林科技大學學報(自然科學版),2020(8):1-8.蘇海民,汪兆國,黃園.宿州市煤礦區(qū)農(nóng)田土壤重金屬污染及生態(tài)風險評價[J].山東師范大學學報(自然科學版),2017,32(02):99-104.AblizAbdugheni,QingdongShi,KeyimuMaierdang.Spatialdistribution,source,andriskassessmentofsoiltoxicmetalsinthecoal-miningregionofnorthwesternChina[J].ArabianJournalofGeosciences,2018,11(24):793.GB15618—2018.土壤環(huán)境質量農(nóng)用地土壤污染風險管控標準[S].生態(tài)環(huán)境部,2018.ChiromaTM,EbeweleRO,HymoreFK.Comparativeassessmentofheavymetallevelsinsoil,vegetablesandurbangreywaterusedforirrigationinYolaandKano.2014.劉玉紅.龍口市北部平原土壤重金屬空間分布特征分析[J].安徽農(nóng)學通報,2020,26(1):98-100.謝龍濤,潘劍君,白浩然,等.基于GIS的農(nóng)田土壤重金屬空間分布及污染評價——以南京市江寧區(qū)某鄉(xiāng)鎮(zhèn)為例[J].土壤學報:1-10.焦海林.城市生活污水和污泥中的重金屬[J].化學工程與裝備,2019(8):319-321.CsavinaJ,FieldJ,TaylorMP,etal.Areviewontheimportanceofmetalsandmetalloidsinatmosphericdustandaerosolfromminingoperations[J].ScienceoftheTotalEnvironment,2012,433:58-73.黃大偉,桂和榮.宿南礦區(qū)土壤重金屬含量特征及其來源解析[J].地球與環(huán)境,2017,45(5):546-554.王鳳英.化肥中的有害元素[J].環(huán)境科學動態(tài),1983(15):8.李昊哲.不同地膜條件下重金屬在農(nóng)田土壤的垂向分布特性研究[D].太原理工大學,2018:25-27.王萍,劉靜,朱健,等.巖溶山區(qū)磷石膏堆場重金屬遷移對耕地質量的影響及污染風險管控[J].水土保持通報,2019,39(04):294-299.BlahaU,AppelE,StanjekH.Determinationofanthropogenicboundarydepthinindustriallypollutedsoilandsemi-quantificationofheavymetalloadsusingmagneticsusceptibility[J].EnvironmentalPollution,2008,156(2):278-289.CuiH,MaK,FanY,etal.Stabilityandheavymetaldistributionofsoilaggregatesaffectedbyapplicationofapatite,lime,andcharcoal[J].EnvironSciPollutResInt,2016,23(11):10808-10817.MaereT,VillezK,Marsili-LibelliS,etal.Membranebioreactorfoulingbehaviourassessmentthroughprincipalcomponentanalysisandfuzzyclustering[J].WaterResearch,2012,46(18):6132-6142.劉斌.基于主成分分析法的中國超大城