于林松，宋蕊蕊, 康桂玲, 范海印，胡尊芳，邵麗娟

(1.山東省物化探勘查院，山東省土壤地球化學(xué)工程實(shí)驗(yàn)室，山東省地質(zhì)勘查工程技術(shù)研究中心，山東 濟(jì)南 250013；2.成都理工大學(xué) 地球科學(xué)學(xué)院，四川 成都 610059)

0 引言

土壤是位于地球陸地表面和淺水域底部具有生命力、生產(chǎn)力的疏松而不均勻的聚積層，是地球系統(tǒng)的組成部分和調(diào)控環(huán)境質(zhì)量的中心要素[1]。隨著工、農(nóng)業(yè)生產(chǎn)現(xiàn)代化，人類的各種生產(chǎn)活動(dòng)已經(jīng)對(duì)地表耕地土壤生態(tài)環(huán)境產(chǎn)生了巨大的破壞。耕地土壤重金屬污染受到各學(xué)科的普遍重視[2-3]。據(jù)統(tǒng)計(jì)(1)中國(guó)地質(zhì)調(diào)查局，中國(guó)耕地地球化學(xué)調(diào)查報(bào)告，2015年。，我國(guó)耕地土壤重金屬中—重度污染或超標(biāo)的點(diǎn)位比例占2.5%，覆蓋面積232.53萬hm2，輕微—輕度污染或超標(biāo)的點(diǎn)位比例占5.7%，覆蓋面積526.8萬hm2。國(guó)家“十三五”規(guī)劃綱要中明確提出，要深入實(shí)施土壤污染防治行動(dòng)計(jì)劃，圍繞最嚴(yán)格的耕地保護(hù)制度，實(shí)施土壤污染分類分級(jí)防治，優(yōu)先保護(hù)農(nóng)用地土壤環(huán)境質(zhì)量安全[4]。有關(guān)耕地重金屬污染問題的評(píng)價(jià)、探討再次引起社會(huì)的廣泛關(guān)注。近年來，我國(guó)學(xué)者在農(nóng)用地分等定級(jí)[5-6]和生態(tài)地球化學(xué)調(diào)查評(píng)價(jià)[7-8]基礎(chǔ)上，逐步形成了耕地質(zhì)量評(píng)價(jià)[9-10]、土壤環(huán)境質(zhì)量[11-12]與土地質(zhì)量地球化學(xué)[13-14]評(píng)價(jià)等綜合方法體系，并形成了以水土重金屬生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)為主的多種評(píng)價(jià)方法[15-21]，主要有污染負(fù)荷指數(shù)法、單因子污染指數(shù)法、內(nèi)梅羅污染指數(shù)法、地累積指數(shù)法、潛在風(fēng)險(xiǎn)評(píng)價(jià)法等。然而，隨著當(dāng)前大比例尺土地質(zhì)量調(diào)查工作的持續(xù)推進(jìn)，面對(duì)調(diào)查精度高、重金屬元素?cái)?shù)據(jù)量大的特點(diǎn)，尚未形成一套快捷有效的系統(tǒng)分區(qū)評(píng)價(jià)方式。該次以郯城縣某鎮(zhèn)為例，嘗試在土壤地球化學(xué)分區(qū)基礎(chǔ)上，進(jìn)行耕地土壤重金屬污染評(píng)價(jià)，以期為鄉(xiāng)鎮(zhèn)級(jí)土地質(zhì)量地球化學(xué)安全性評(píng)價(jià)提供參考，也為區(qū)域土地資源污染分級(jí)防控提供科學(xué)依據(jù)。

1 材料與方法

1.1 研究區(qū)概況

研究區(qū)位于山東省臨沂市郯城縣某鄉(xiāng)鎮(zhèn)，下轄29個(gè)行政村委會(huì)，面積為60.04km2，其中，耕地面積約46.94km2。按地形地貌特征看，研究區(qū)屬于沂沭河沖積平原。屬東亞暖溫帶季風(fēng)氣候，四季分明，雨熱同季。研究區(qū)大地構(gòu)造單元位于華北板塊魯中隆起區(qū)的馬頭凹陷和沂沭斷裂帶的郯城凹陷附近。區(qū)內(nèi)未見基巖出露，第四系覆蓋廣泛，厚度大。研究區(qū)耕地類型主要為水澆地、水田、旱地和園地，稻麥輪作；土壤類型以水稻土為主，少量潮土；成土母質(zhì)主要為第四紀(jì)臨沂組沖積相沉積物和黑土湖組湖沼相沉積物。研究區(qū)是郯城縣傳統(tǒng)的農(nóng)業(yè)鄉(xiāng)鎮(zhèn)，水稻種植歷史悠久，素有“魯南米倉(cāng)”之稱；工業(yè)不甚發(fā)達(dá)，以輕工業(yè)為主，次為零散的鄉(xiāng)鎮(zhèn)企業(yè)。

1.2 樣品采集與測(cè)試

按照相關(guān)地質(zhì)調(diào)查規(guī)范[22]要求，以網(wǎng)格加地塊、均勻控制的采樣思路，按1∶10 000調(diào)查精度采集樣品，采樣區(qū)域?yàn)檑俺强h某鎮(zhèn)全域。研究區(qū)采集表層土壤(0～20cm)樣品1279件，其中耕地土壤樣品1150件，采集密度為24點(diǎn)/km2。野外用GPS精確定位采樣點(diǎn)，現(xiàn)場(chǎng)填寫“土壤樣品采樣記錄卡”。采樣過程中，避開外來土和新近擾動(dòng)過的土層，去掉表面雜物和土壤中的碎石及動(dòng)植物殘?bào)w，按五點(diǎn)等量組合采樣，四分法保留1.5kg新鮮土樣。

土壤樣品經(jīng)自然風(fēng)干，先過10目(孔徑2mm)篩，然后縮分法取10g樣品研磨過100目(孔徑0.149mm)篩，用于重金屬元素分析。土壤樣品消解方法：稱取0.2 g土壤樣于聚四氟乙烯坩堝中，加入混合酸(HF-HClO4-HNO3)進(jìn)行消解。重金屬元素分析方法為：①Pb，Ni，Cr，Zn，Cu采用X-射線熒光光譜法測(cè)定(Axiosmax,PANalytical B.V,荷蘭)；②Cd采用石墨爐原子吸收光譜法測(cè)定(GF-AAS,PE600, thermo Elemental,美國(guó))；③As，Hg采用氫化物發(fā)生-原子熒光光譜法測(cè)定(AFS9750，北京海光儀器,中國(guó))。試驗(yàn)所用玻璃器皿均用10%硝酸溶液浸泡24h，所有試劑均為優(yōu)級(jí)純，分析用水均為超純水。測(cè)定過程中，每50件樣品插入一件重復(fù)樣和4件土壤國(guó)家標(biāo)準(zhǔn)物質(zhì)(GBW—07403)進(jìn)行質(zhì)量控制，樣品重復(fù)性檢驗(yàn)合格率91.25%～98.75%，異常點(diǎn)重復(fù)檢查合格率96.00%～100%，標(biāo)樣測(cè)定結(jié)果均在允許誤差范圍內(nèi)。

1.3 評(píng)價(jià)方法

1.3.1 土壤地球化學(xué)分區(qū)

元素地球化學(xué)分區(qū)在20世紀(jì)八九十年代的區(qū)域地球化學(xué)調(diào)查研究中被廣泛應(yīng)用[23]，為直接或間接快速、便捷指示地質(zhì)礦產(chǎn)及環(huán)境信息提供了翔實(shí)的依據(jù)[24]。后來，地球化學(xué)分區(qū)也被表生地球化學(xué)研究、農(nóng)業(yè)及環(huán)境領(lǐng)域的元素地球化學(xué)研究所應(yīng)用[23]。土壤地球化學(xué)調(diào)查數(shù)據(jù)是土壤地球化學(xué)分區(qū)的關(guān)鍵。因此，土壤地球化學(xué)分區(qū)前提是必須要掌握一個(gè)地區(qū)系統(tǒng)的土壤元素含量分布等調(diào)查數(shù)據(jù)[25]。進(jìn)行土壤地球化學(xué)分區(qū)，其實(shí)質(zhì)就是依據(jù)土壤元素地球化學(xué)特征的差異所反映的土壤化學(xué)元素組合，用統(tǒng)一的地球化學(xué)標(biāo)準(zhǔn)劃分不同的土壤地球化學(xué)范圍。土壤地球化學(xué)分區(qū)內(nèi)的元素間彼此存在相關(guān)性，而重金屬元素的含量往往決定了該土壤地球化學(xué)分區(qū)的生態(tài)風(fēng)險(xiǎn)狀況。因此，可以通過重金屬元素的含量狀況作出快捷的生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)。

1.3.2 土壤重金屬污染指數(shù)[26]

單項(xiàng)污染指數(shù)法是評(píng)價(jià)土壤中單個(gè)重金屬污染元素的污染指數(shù)。其計(jì)算公式為：

(1)

式中：Pi為土壤中污染物i的單個(gè)重金屬污染元素的污染指數(shù)；Ci為污染物i的實(shí)測(cè)濃度；Si為污染物i的評(píng)價(jià)標(biāo)準(zhǔn)。指數(shù)(Pi)越大表明土壤重金屬累積程度越高。

當(dāng)區(qū)域內(nèi)土壤環(huán)境質(zhì)量作為一個(gè)整體與外區(qū)域進(jìn)行比較時(shí)，除用單項(xiàng)污染指數(shù)外，還常用綜合污染指數(shù)，而其中最常用的為內(nèi)梅羅污染指數(shù)。土壤由于地區(qū)背景差異較大，用土壤污染累積指數(shù)更能反映土壤的人為污染程度。其計(jì)算公式為：

(2)

式中：P為綜合污染指數(shù)；Pmean為各單因子污染指數(shù)的平均值；Pimax為各單因子污染指數(shù)中的最大值?；诰C合污染指數(shù)的土壤污染等級(jí)劃分標(biāo)準(zhǔn)見表1。

表1 土壤內(nèi)梅羅污染指數(shù)評(píng)價(jià)標(biāo)準(zhǔn)

1.3.3 地累積指數(shù)法

地累積指數(shù)由德國(guó)海德堡大學(xué)的Muller等(1969)提出，通常稱為Muller指數(shù)。它是一種研究土壤、沉積物中重金屬污染程度的定量指標(biāo)，已被國(guó)內(nèi)外學(xué)者[15，18-19]廣泛應(yīng)用于人為活動(dòng)產(chǎn)生的重金屬對(duì)土壤污染的評(píng)價(jià)。其公式為：

(3)

式中：Ci是樣品中元素i的實(shí)測(cè)濃度；Bi是土壤中元素i的地球化學(xué)背景值。重金屬地累積指數(shù)(Igeo)分級(jí)與污染程度的關(guān)系列于表2。

表2 地累積指數(shù)的污染級(jí)別標(biāo)準(zhǔn)

2 結(jié)果與討論

2.1 土壤地球化學(xué)分區(qū)

研究區(qū)位于山東省南部與江蘇新沂市接壤(圖1)，成土母質(zhì)屬區(qū)域性“新沂-贛榆北部殘坡積-沖積亞區(qū)[25]”之沂沭河沖積體。以此次土壤地球化學(xué)調(diào)查數(shù)據(jù)為基礎(chǔ)，以SPSS 19.0為處理工具，采用主成分分析方法，對(duì)研究區(qū)1150件表層土壤樣品的22種測(cè)試指標(biāo)的檢測(cè)數(shù)據(jù)進(jìn)行R型因子分析。首先，對(duì)變量樣本進(jìn)行數(shù)據(jù)標(biāo)準(zhǔn)化處理，同時(shí)進(jìn)行取樣足夠度的度量和顯著性檢驗(yàn)，KMO值為0.790(KMO>0.50，其值愈逼近1，表明對(duì)這些變量進(jìn)行因子分析的效果愈好)，Bartlett球形檢驗(yàn)顯著(Sig.=0.000<0.05)，認(rèn)為變量樣本數(shù)據(jù)因子分析的效果非常好，其分析結(jié)果能較好地反應(yīng)元素之間成因上的聯(lián)系；其次，對(duì)因子載荷矩陣進(jìn)行正交旋轉(zhuǎn)(表3)，在正交旋轉(zhuǎn)載荷矩陣中，將取值大于0.5的變量作為該因子的主要載荷元素組，即認(rèn)為該組元素是具有成因聯(lián)系的多種元素，而初始特征值大于1的前6個(gè)因子(F1～F6)認(rèn)為是“有意義”的主要影響因子，計(jì)算其累計(jì)方差貢獻(xiàn)率達(dá)69.401%。

從研究區(qū)土壤類型及微地貌景觀分區(qū)等表生環(huán)境出發(fā)，綜合考慮pH分等(研究區(qū)耕地土壤pH值在6.15～6.69間，為偏弱酸性，值域變化不大，認(rèn)為是適合小麥、水稻生長(zhǎng)的pH范圍)，疊加6個(gè)主成分因子(F1～F6)得分累頻地球化學(xué)圖。

1—薛家寨子-玉皇廟弱酸性土壤地球化學(xué)分區(qū)；2—墨河小流域中至弱酸性土壤地球化學(xué)分區(qū)；3—興旺村-東樊村弱酸性土壤地球化學(xué)分區(qū)；4—土壤地球化學(xué)分區(qū)界線；5—土壤地球化學(xué)亞區(qū)界線；6—土壤地球化學(xué)亞區(qū)編號(hào)圖1 研究區(qū)位置、微地貌景觀及土壤地球化學(xué)分區(qū)圖 注：圖中白色區(qū)域?yàn)槟有×饔?土壤類型為黏質(zhì)潮土)，灰色區(qū)域?yàn)槲⑿逼皆瓍^(qū)(土壤類型為水稻土)

表3 因子分析旋轉(zhuǎn)成份矩陣

元素主成分因子F1F2F3F4F5F6元素主成分因子F1F2F3F4F5F6Fe2O30.9410.127-0.0400.089-0.060-0.027F-0.0970.959-0.001-0.010-0.0350.038V0.8900.0440.0420.161-0.315-0.028Mn0.1630.0510.8490.118-0.1230.032As0.8290.0400.2030.0010.112-0.005I0.021-0.0040.794-0.0280.140-0.034N0.7120.049-0.4590.332-0.1650.022P0.112-0.011-0.1690.885-0.101-0.021Se0.6990.060-0.0710.2560.011-0.080Zn0.4380.0920.2200.7370.132-0.004Cr0.6830.6480.0370.124-0.053-0.010Cu0.5670.0350.1490.6080.1360.035Corg0.6780.058-0.5180.262-0.1410.043Cd0.3290.0350.0130.3880.2020.305Co0.6390.3080.3870.114-0.1360.004K2O-0.228-0.1790.027-0.0420.732-0.068Pb0.635-0.007-0.0470.0640.1320.170Mo0.0850.018-0.0980.1640.5710.012Ge0.529-0.0610.1730.180-0.136-0.004Hg-0.062-0.0550.0220.136-0.1430.857Ni0.1980.9670.0470.0320.0370.025B-0.054-0.1490.0380.196-0.249-0.400

注：提取方法—主成分；旋轉(zhuǎn)法—具有Kaiser標(biāo)準(zhǔn)化的正交旋轉(zhuǎn)法，旋轉(zhuǎn)在6次迭代后收斂。

表4 土壤地球化學(xué)分區(qū)及其特征描述

具體為：以提取的6主成分因子累頻75%以上的中高值域?yàn)榀B加線，進(jìn)行分區(qū)邊界描繪；同時(shí)根據(jù)pH與地貌等綜合特征對(duì)分區(qū)界線進(jìn)行適當(dāng)完善。比如，圖2展示了F1，F(xiàn)2，F(xiàn)3三個(gè)主成分因子得分計(jì)量圖，展現(xiàn)了不同元素組合在空間分布上呈現(xiàn)顯著的差異；F1因子代表的Fe2O3，V，As，N，Se等10個(gè)元素沿墨河小流域呈現(xiàn)明顯的低值區(qū)，向兩側(cè)逐步過度為背景及中高值分布區(qū)，相關(guān)性顯著以外圍疊加線描繪分區(qū)邊界，可對(duì)應(yīng)典型分區(qū)為Ⅱ-1亞區(qū)和Ⅲ-4亞區(qū)。根據(jù)以上分區(qū)原則，將研究區(qū)土壤劃分為3個(gè)土壤地球化學(xué)分區(qū)、10個(gè)土壤地球化學(xué)亞區(qū)(圖1、表4)。每個(gè)土壤地球化學(xué)分區(qū)或亞區(qū)是相關(guān)因子土壤化學(xué)元素組的反映(圖2)，元素間彼此存在地球化學(xué)相關(guān)性，而重金屬元素的含量往往決定了該土壤地球化學(xué)分區(qū)或亞區(qū)的污染等級(jí)。

2.2 重金屬元素含量特征

從研究區(qū)土壤地球化學(xué)分區(qū)重金屬元素平均含量來看(表5)，該鎮(zhèn)耕地土壤分區(qū)中重金屬含量均值均低于標(biāo)準(zhǔn)[25]規(guī)定的農(nóng)用地風(fēng)險(xiǎn)管控值[27]，土壤環(huán)境質(zhì)量重金屬污染風(fēng)險(xiǎn)低。其中，I-1，III-1，III-2，III-3，III-4土壤分區(qū)樣品均值多高于全區(qū)均值；I-2，I-4，II-1，II-2土壤分區(qū)樣品均值多低于全區(qū)均值；Hg，Cr，Pb及Cu，Ni，Zn等元素含量均值均高于或不同程度高于山東省背景值。從研究區(qū)土壤重金屬元素含量變化系數(shù)分布情況來看，8項(xiàng)重金屬元素變化系數(shù)絕大多數(shù)值域分布在9.8%～35.5%間，屬弱至低等變異程度，表明數(shù)據(jù)總體離散度不大，均一化程度較高；而圖3中明顯反應(yīng)I-3亞區(qū)Ni和II-1亞區(qū)Hg的強(qiáng)變異特征，其變化系數(shù)分別達(dá)241.7%和136.1%，這可能代表外界干擾成為影響元素含量分異的主導(dǎo)因素；同時(shí)，就II-1，I-3兩個(gè)土壤亞區(qū)對(duì)應(yīng)土壤類型分布面積和重金屬污染影響程度來看，黏質(zhì)潮土重金屬風(fēng)險(xiǎn)面積大于水稻土分布面積(表4)。此外，就全區(qū)Hg含量分布來看，介于0.028～0.361mg/kg間，高于背景值的采樣點(diǎn)超標(biāo)率達(dá)到73.3%，這表明表層土壤中的Hg元素受后期人為活動(dòng)的干擾影響較大。

圖2 研究區(qū)典型因子得分計(jì)量圖

表5 土壤地球化學(xué)分區(qū)重金屬元素平均含量(mg/kg)

元素均值統(tǒng)計(jì)范圍I-1I-2I-3I-4II-1II-2III-1III-2III-3III-4全區(qū)山東省表層土壤地球化學(xué)背景值[28]農(nóng)用地風(fēng)險(xiǎn)管控值*As7.777.867.677.876.386.497.446.598.87.917.1448.625Hg0.040.030.040.030.050.040.030.040.040.040.0410.0310.5Cd0.130.110.130.130.130.120.120.110.110.140.1280.1320.3Cr70.7164.9779.7465.8362.0767.6670.7570.3874.7478.7967.95162150Cu24.1921.4522.8720.6721.9721.923.9623.825.0326.323.17122.650Ni27.0325.0450.1424.0222.9324.9226.6726.6227.7131.4226.03627.170Pb28.0125.928.726.2425.6326.2327.327.3528.8430.3327.04223.690Zn62.553.4862.1151.056259.3764.265.3163.2870.0862.4863.3200

注：“*”為標(biāo)準(zhǔn)[27]“表1 農(nóng)用地土壤污染風(fēng)險(xiǎn)篩選值”中“5.5

圖3 不同土壤地球化學(xué)分區(qū)重金屬元素變異系數(shù)對(duì)比圖

2.3 評(píng)價(jià)結(jié)果

2.3.1 土壤重金屬污染指數(shù)評(píng)價(jià)

根據(jù)研究區(qū)不同土壤地球化學(xué)分區(qū)重金屬元素含量，采用表5中“農(nóng)用地風(fēng)險(xiǎn)管控值”作為該區(qū)土壤環(huán)境質(zhì)量評(píng)價(jià)標(biāo)準(zhǔn)。將土壤亞區(qū)重金屬元素平均含量帶入公式(1)(2)中，計(jì)算結(jié)果見表6。內(nèi)梅羅污染指數(shù)中只有I-3亞區(qū)的值超過0.5，其余分區(qū)均在0.5以下，均處于土壤清潔等級(jí)(I級(jí))。從內(nèi)梅羅污染指數(shù)回看單項(xiàng)污染指數(shù)可以看出，在各土壤亞區(qū)中8項(xiàng)重金屬的單項(xiàng)污染指數(shù)最高的為I-3亞區(qū)中Ni的0.716和Cr的0.532。經(jīng)查證，這是I-3亞區(qū)內(nèi)的單點(diǎn)異常引起，Cr，Ni處于同點(diǎn)位，表明存在單個(gè)因子人為污染的現(xiàn)象，這可能與農(nóng)田施肥及噴灑農(nóng)藥等相關(guān)。

表6 土壤地球化學(xué)分區(qū)重金屬污染指數(shù)

表7 土壤地球化學(xué)分區(qū)重金屬地累積指數(shù)

2.3.2 地累積指數(shù)評(píng)價(jià)

以山東省表層土壤地球化學(xué)背景值[27]作為土壤背景參比值進(jìn)行計(jì)算。將土壤亞區(qū)重金屬元素平均含量帶入公式(3)，計(jì)算結(jié)果見表7。按照木桶原理，取各種元素污染的最高級(jí)別表示該土壤亞區(qū)的污染程度。比如，I-3亞區(qū)重金屬元素Igeo排序?yàn)椋篘i(0.303)>Hg(-0.131)>Cr(-0.222)>Pb(-0.302)>Cu(-0.568)>Zn(-0.612)>Cd(-0.664)>As(-0.750)，其中，Ni元素的Igeo>0，區(qū)域富集污染程度為輕微污染，其他元素均屬未污染(Igeo<0)，按最高級(jí)污染級(jí)別該亞區(qū)為輕微污染級(jí)；同理，II-1亞區(qū)中Hg元素的Igeo>0，亦為輕微污染，重金屬元素Igeo排序?yàn)椋篐g(0.019)>Pb(-0.466)>Cr(-0.538)>Zn(-0.615)>Cd(-0.618)>Cu(-0.626)>Ni(-0.826)>As(-1.016)。在其余的I-1、I-2、Ⅱ-2、Ⅲ-1、Ⅲ-2、Ⅲ-3、Ⅲ-4亞區(qū)中重金屬元素均屬未污染級(jí)別。

2.4 評(píng)價(jià)結(jié)果一致性討論

利用金維地學(xué)軟件Geo IPAS 4.0軟件分別繪制耕地土壤重金屬累加地球化學(xué)圖(圖4)和耕地土壤環(huán)境質(zhì)量綜合等級(jí)圖(圖5)。從重金屬累加圖上可以看出，累加值350 mg/kg以上的區(qū)域主要集中在I-3、II-1土壤地球化學(xué)分區(qū)，其中I-3分區(qū)累積值濃集中心顯著，極值為1028.866 mg/kg，II-1分區(qū)沿墨河小流域呈現(xiàn)串珠狀的多個(gè)濃集點(diǎn)位，極值為1023.664 mg/kg。從耕地土壤環(huán)境質(zhì)量綜合分等圖上可以看出，玉皇廟村附近I-3分區(qū)耕地土壤存在中—重度污染點(diǎn)位，對(duì)應(yīng)重金屬累加地球化學(xué)的單點(diǎn)異常點(diǎn)位；沿墨河小流域II-1分區(qū)存在輕微—輕度污染等級(jí)地塊，對(duì)應(yīng)墨河附近的多個(gè)濃集點(diǎn)位。總體來看，土壤地球化學(xué)亞區(qū)潛在生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)信息與圖件反映信息的一致性較好，均指示了墨河小流域和玉皇廟村附近耕地存在一定的生態(tài)風(fēng)險(xiǎn)。

圖中“+”即重金屬累加高值點(diǎn)圖4 土壤重金屬累加地球化學(xué)圖

圖5 土壤環(huán)境質(zhì)量綜合分等圖

3 結(jié)語

通過該鄉(xiāng)鎮(zhèn)此次大比例尺土地質(zhì)量調(diào)查區(qū)耕地土壤中重金屬評(píng)價(jià)得到以下結(jié)論：

(1)綜合考慮研究區(qū)土壤類型及微地貌景觀，利用多因子疊加的方式，劃分了3個(gè)土壤地球化學(xué)分區(qū)，10個(gè)土壤地球化學(xué)亞區(qū)；耕地土壤重金屬含量分布結(jié)果顯示，8項(xiàng)重金屬元素含量屬弱至低等變異程度，元素含量均一化程度高，環(huán)境質(zhì)量總體處于清潔；但個(gè)別土壤分區(qū)存在Hg，Ni，Cr含量分異程度大，點(diǎn)位高于背景比例偏高的現(xiàn)象，黏質(zhì)潮土重金屬污染風(fēng)險(xiǎn)面積大于水稻土分布面積。

(2)各土壤亞區(qū)重金屬平均含量均低于當(dāng)前農(nóng)用地土壤污染風(fēng)險(xiǎn)篩選值，總體處于土壤清潔等級(jí)；而單項(xiàng)污染指數(shù)中Ni、Cr的值相對(duì)偏高，表明存在單個(gè)因子人為污染的現(xiàn)象；內(nèi)梅羅污染指數(shù)指示I-3亞區(qū)為該鎮(zhèn)最高，存在一定的污染隱患；地累積指數(shù)指示了I-3亞區(qū)中Ni、II-1亞區(qū)中Hg>1，土壤亞區(qū)為輕微污染等級(jí)。

(3)基于土壤地球化學(xué)分區(qū)的綜合評(píng)價(jià)結(jié)果表明，基于土壤地球化學(xué)分區(qū)的耕地土壤重金屬污染評(píng)價(jià)結(jié)果與土壤環(huán)境質(zhì)量評(píng)價(jià)圖件反映的信息較為吻合。在后期開展的大比例尺土地質(zhì)量調(diào)查工作中，可以采用基于土壤地球化學(xué)分區(qū)的形式進(jìn)行耕地土壤污染風(fēng)險(xiǎn)評(píng)價(jià)，并可據(jù)此次評(píng)價(jià)結(jié)果對(duì)土壤地球化學(xué)亞區(qū)重金屬污染情況進(jìn)行分類分級(jí)防控，對(duì)單項(xiàng)指數(shù)偏高的亞區(qū)進(jìn)行采樣點(diǎn)位的分級(jí)評(píng)價(jià)，以了解實(shí)際的污染風(fēng)險(xiǎn)水平。