一種基于風(fēng)險(xiǎn)管控的稻田土壤重金屬污染分級方法

王琦，李芳柏*，黃小追，馬義兵，林玉鎖，劉曉文

1. 廣東省生態(tài)環(huán)境技術(shù)研究所/廣東省農(nóng)業(yè)環(huán)境綜合治理重點(diǎn)實(shí)驗(yàn)室，廣東 廣州 510650；2. 中國農(nóng)業(yè)科學(xué)院農(nóng)業(yè)資源與農(nóng)業(yè)區(qū)劃研究所，北京 100081；3. 生態(tài)環(huán)境部南京環(huán)境科學(xué)研究所，江蘇 南京 210042；4. 生態(tài)環(huán)境部華南環(huán)境科學(xué)研究所，廣東 廣州 510655

土壤是人類賴以生存和國家文明建設(shè)的基礎(chǔ)性自然資源（陳懷滿等，1999；張甘霖等，2018；朱永官等，2015）。土壤是所有陸地生態(tài)系統(tǒng)的基礎(chǔ)，土壤環(huán)境受到污染，會威脅生態(tài)系統(tǒng)的結(jié)構(gòu)和功能。近 30年來，隨著社會經(jīng)濟(jì)和農(nóng)業(yè)集約化的快速發(fā)展，中國土壤污染逐年加劇，尤其是農(nóng)產(chǎn)品產(chǎn)地土壤的污染，威脅食物鏈安全和人體健康（Huang et al.，2018；Li et al.，2012；鄧琳靜，2016；樊霆等，2013；顧繼光等，2003；邱孟龍等，2017）。2014年的《全國土壤污染狀況調(diào)查公報(bào)》顯示，全國土壤環(huán)境狀況總體不容樂觀，總超標(biāo)率為16.1%；其中，中度和重度污染點(diǎn)位比例分別為1.5%和1.1%（環(huán)境保護(hù)部等，2014）。為此，國務(wù)院印發(fā)了《土壤污染防治行動計(jì)劃》（簡稱“土十條”），作為當(dāng)前乃至今后一個時期內(nèi)全國土壤污染防治工作的行動綱領(lǐng)，其目標(biāo)是至 2020年，受污染耕地安全利用率達(dá)到90%左右（王國慶等，2016）。2018年8月31日，《中華人民共和國土壤污染防治法》（以下簡稱《土壤污染防治法》）正式審議通過，在農(nóng)用地土壤污染風(fēng)險(xiǎn)管控方面，《土壤污染防治法》規(guī)定由國家建立農(nóng)用地分類管理制度；按照土壤污染程度和相關(guān)標(biāo)準(zhǔn)，將農(nóng)用地劃分為優(yōu)先保護(hù)類、安全利用類和嚴(yán)格管控類（駱永明等，2018）。因此，系統(tǒng)認(rèn)識中國農(nóng)田土壤污染區(qū)域化等級化特征，探索基于風(fēng)險(xiǎn)管控的農(nóng)田土壤重金屬污染分級方法，實(shí)現(xiàn)土壤污染分類分級化治理和修復(fù)顯得重要且迫切（Liu et al.，2016；Zhang et al.，2018；蘇耀明等，2016）。

然而，土壤污染分類分級存在以下2個技術(shù)難點(diǎn)尚未解決，第一，如何劃分污染程度：優(yōu)先保護(hù)類、安全利用類和嚴(yán)格管控類？第二，如何實(shí)現(xiàn)土壤與農(nóng)產(chǎn)品重金屬污染風(fēng)險(xiǎn)的協(xié)同評價？這些難點(diǎn)嚴(yán)重阻礙了土壤分區(qū)分級治理方法與行動方案的建立，需要綜合考慮土壤和農(nóng)產(chǎn)品質(zhì)量及環(huán)境要素才能對土壤環(huán)境質(zhì)量進(jìn)行科學(xué)的風(fēng)險(xiǎn)評價（Li et al.，2017；Luo et al.，2016；Yang et al.，2018；劉蕊等，2014；秦普豐等，2010）。本研究旨在建立一套區(qū)域稻田土壤重金屬污染風(fēng)險(xiǎn)分級方法，針對珠三角稻田土壤重金屬Cd污染，采用土壤單因子指數(shù)、農(nóng)產(chǎn)品單因子指數(shù)和富集系數(shù)相結(jié)合的方法，基于區(qū)域自然地理要素和污染源等因素劃分評價的地理單元，進(jìn)一步結(jié)合點(diǎn)位土壤、農(nóng)產(chǎn)品和富集系數(shù)風(fēng)險(xiǎn)劃分等級，對區(qū)域稻田土壤重金屬污染風(fēng)險(xiǎn)進(jìn)行分級，分級分類結(jié)果將有利于針對性的土壤環(huán)境管理，符合當(dāng)前中國土壤污染防治工作的目標(biāo)和任務(wù)，是對現(xiàn)有土壤環(huán)境質(zhì)量分級方法的科學(xué)探索和有效補(bǔ)充。

1 材料和方法

1.1 研究區(qū)域概況

珠江三角洲行政區(qū)域包括廣州市、深圳市、珠海市、佛山市、惠州市、東莞市、中山市、江門市和肇慶市等9市，面積約4.17×104km2，約占廣東省陸域面積的23.2%（張文忠等，2003）。珠江三角洲屬亞熱帶季風(fēng)氣候，以珠江三角洲平原為主體，由西江、北江、東江、流溪河等諸河構(gòu)成的珠江水系在區(qū)內(nèi)縱橫交錯（黎夏，2004）。土地利用方式以林地為主，其次是耕地、水域和居民點(diǎn)。其中，耕地以灌溉水田、水澆地、旱地、菜地為主，主要分布于沖積平原區(qū)。主要農(nóng)作物以水稻和蔬菜為主，水稻主要分布在東江、西江、北江流域（曾輝等，1999）。據(jù)統(tǒng)計(jì)，2015年珠江三角洲耕地面積為6.1×105hm2。區(qū)內(nèi)土壤類型以赤紅壤和水稻土分布最廣，分別占陸域總面積的 44.8%和40.2%。其次為紅壤和潮土，所占比例分別為6.5%和4.8%，其他類型土壤所占比例均不足 1.0%（曹祺文等，2016；宗慶霞等，2017）。該區(qū)是中國經(jīng)濟(jì)最發(fā)達(dá)、增長速度最快的地區(qū)之一，工業(yè)化和城市化的迅速發(fā)展產(chǎn)生了大量工業(yè)廢物和生活垃圾，導(dǎo)致出現(xiàn)了嚴(yán)重的重金屬污染問題（陸大道，2017）。

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

農(nóng)田土壤、水稻重金屬點(diǎn)位數(shù)據(jù)來自 2011年以來的廣東省生態(tài)環(huán)境技術(shù)研究所珠三角土壤污染狀況調(diào)查和文獻(xiàn)數(shù)據(jù)，共1566個水稻產(chǎn)地土壤點(diǎn)位數(shù)據(jù)，其中土壤-水稻點(diǎn)位數(shù)據(jù)880個，土地利用現(xiàn)狀圖來自2010年廣東省第二次全國土地調(diào)查。采樣點(diǎn)、涉重金屬污染企業(yè)（340個國控和省控企業(yè)，數(shù)據(jù)來源：原國家環(huán)境保護(hù)部）和海拔空間分布如圖1所示。借助SPSS 13.0進(jìn)行數(shù)據(jù)描述性統(tǒng)計(jì)分析、正態(tài)分布檢驗(yàn)、方差分析后，觀察到數(shù)據(jù)的變異性較大，利用Microsoft Excel 2010采用4倍均方差法剔除了數(shù)據(jù)中的異常值，對土壤和稻米數(shù)據(jù)進(jìn)行對數(shù)轉(zhuǎn)換，使其服從正態(tài)分布（Eischeid et al.，1995）。

圖1 珠三角采樣點(diǎn)、基礎(chǔ)地理及涉重金屬企業(yè)分布Fig. 1 Sampling sites, basic geography and distribution of heavy metal enterprises in the Pearl River Delta

1.3 土壤重金屬污染分級方法

1.3.1 方法的背景

根據(jù)針對性、科學(xué)性、可操作性的原則建立污染分級方法。

首先，全面收集區(qū)域概況資料、土壤污染源資料、土壤環(huán)境和農(nóng)產(chǎn)品資料及其對應(yīng)的圖件資料等，并通過分析，研究有關(guān)信息劃分區(qū)域內(nèi)的地理單元，為污染風(fēng)險(xiǎn)評價提供基礎(chǔ)支撐；其次，通過評估土壤及農(nóng)產(chǎn)品調(diào)查點(diǎn)位的土壤重金屬污染情況、農(nóng)產(chǎn)品重金屬污染情況，反映該點(diǎn)位土壤和農(nóng)產(chǎn)品重金屬污染狀況，劃分點(diǎn)位土壤重金屬風(fēng)險(xiǎn)等級，為區(qū)域風(fēng)險(xiǎn)評價提供基礎(chǔ)依據(jù)；最后，通過評價單元內(nèi)土壤重金屬污染風(fēng)險(xiǎn)、農(nóng)產(chǎn)品重金屬危害風(fēng)險(xiǎn)、重金屬生物可利用性等因素的情況，反映區(qū)域內(nèi)耕地土壤和農(nóng)產(chǎn)品重金屬污染情況，劃分區(qū)域內(nèi)土壤重金屬風(fēng)險(xiǎn)等級，為重金屬污染耕地分級管理提供依據(jù)。

1.3.2 地理單元劃分方法

根據(jù)行政區(qū)域、流域水系、氣候氣象、地形地貌、水文地質(zhì)以及成土母質(zhì)、土壤類型等要素的相對一致性劃分地理單元。原則上，每個地理單元必須僅屬于同一行政邊界范圍內(nèi)，便于地塊的屬地管理。

1.3.3 點(diǎn)位風(fēng)險(xiǎn)評價方法

（1）等級劃分指標(biāo)

采用土壤單因子指數(shù)和農(nóng)產(chǎn)品單因子指數(shù)相結(jié)合的方法（Wang et al.，2018）。土壤指數(shù)（Pi）和農(nóng)產(chǎn)品指數(shù)（Ei）計(jì)算公式分別如下所示：

式中，Pi為土壤重金屬i的單因子指數(shù)；Ci為土壤重金屬i的實(shí)測濃度；Coi為土壤重金屬i的污染風(fēng)險(xiǎn)評估參比值，采用GB 15618（生態(tài)環(huán)境部等，2018）中規(guī)定的土壤污染風(fēng)險(xiǎn)篩選值。

式中，Ei為農(nóng)產(chǎn)品中重金屬i的單因子指數(shù)；Ai為農(nóng)產(chǎn)品中重金屬i的實(shí)測濃度；Si為農(nóng)產(chǎn)品重金屬i的限量標(biāo)準(zhǔn)值，采用GB 2762（中華人民共和國衛(wèi)生部，2012）中規(guī)定的食物中污染物限量標(biāo)準(zhǔn)。

（2）等級劃分

依據(jù)土壤指數(shù)和農(nóng)產(chǎn)品指數(shù)將點(diǎn)位稻田土壤重金屬風(fēng)險(xiǎn)等級劃分為6類7級，劃分依據(jù)見表1。

表1 稻田土壤重金屬Cd點(diǎn)位風(fēng)險(xiǎn)等級劃分依據(jù)Table 1 Risk classification based on heavy metal Cd points in paddy soil

1.3.4 區(qū)域風(fēng)險(xiǎn)評價方法

1.3.4.1 風(fēng)險(xiǎn)評價單元劃分

基于已有資料及圖件，采用地統(tǒng)計(jì)學(xué)、空間疊置等空間風(fēng)險(xiǎn)方法，依次疊置求交集劃定相關(guān)類別單元的邊界，最后劃定為相對均一的風(fēng)險(xiǎn)評價單元。在地理單元的基礎(chǔ)上，對受同一污染源影響、且污染程度相似的區(qū)域原則上劃為同一風(fēng)險(xiǎn)評價單元，具體需綜合考慮點(diǎn)位風(fēng)險(xiǎn)評價結(jié)果、區(qū)域污染源類型因素，結(jié)合土地利用方式、地塊邊界等因素綜合確定風(fēng)險(xiǎn)評價單元的區(qū)域和邊界（圖2）。

劃分方法具體如下：

圖2 基于點(diǎn)位評價的風(fēng)險(xiǎn)評價單元范圍加載示意圖Fig. 2 Diagram of risk assessment unit based on point-based risk assessment

（1）底圖搭建。工作底圖包括行政區(qū)劃、河流水系、交通（公路）、土壤類型分布、全國第二次土地利用現(xiàn)狀調(diào)查更新數(shù)據(jù)、地形地貌、植被類型及覆蓋率分布等基礎(chǔ)性圖件以及遙感影像圖。

（2）范圍加載。根據(jù)點(diǎn)位風(fēng)險(xiǎn)評價結(jié)果，加載點(diǎn)位風(fēng)險(xiǎn)評價結(jié)果，初步劃定風(fēng)險(xiǎn)評價單元。按照主導(dǎo)性原則，若每項(xiàng)污染物80%以上的土壤點(diǎn)位分類結(jié)果一致時，則采用該結(jié)果判定該項(xiàng)污染物所代表的評價單元類別。

（3）邊界劃定。以市（縣、區(qū)）為單位，根據(jù)點(diǎn)位評估情況、重點(diǎn)污染源，按照上述地理單元的對應(yīng)關(guān)系，結(jié)合實(shí)際地勢、地貌、農(nóng)用地分布等情況，劃定評價單元。

1.3.4.2 風(fēng)險(xiǎn)評價指標(biāo)

采用農(nóng)產(chǎn)品重金屬危害風(fēng)險(xiǎn)、土壤重金屬污染風(fēng)險(xiǎn)、重金屬生物可利用性和土壤重金屬污染源輸入風(fēng)險(xiǎn)相結(jié)合的評估方法，評估指標(biāo)包括農(nóng)產(chǎn)品重金屬污染指數(shù)（E）、土壤重金屬污染指數(shù)（P）、重金屬生物可利用風(fēng)險(xiǎn)商（QBCF）、土壤污染源輸入指數(shù)（Qs）。

（1）農(nóng)產(chǎn)品污染指數(shù)（E）計(jì)算公式（Moore，2010）如下：

式中，E為評價區(qū)域農(nóng)產(chǎn)品中重金屬的污染指數(shù)；A為研究區(qū)農(nóng)產(chǎn)品中重金屬的實(shí)測質(zhì)量分?jǐn)?shù)均值，mg·kg-1；SE為研究農(nóng)產(chǎn)品中重金屬的實(shí)測值的標(biāo)準(zhǔn)誤差；S為農(nóng)產(chǎn)品重金屬的限量標(biāo)準(zhǔn)值，采用 GB 2762中的農(nóng)產(chǎn)品中污染物限量標(biāo)準(zhǔn)，mg·kg-1；

（2）土壤重金屬污染指數(shù)（P）計(jì)算公式如下：

式中，P為評價區(qū)域土壤中重金屬的污染指數(shù)；C為評價單元土壤中重金屬的實(shí)測質(zhì)量分?jǐn)?shù)均值，mg·kg-1；Cs為土壤重金屬的限量值，采用 GB 15618中規(guī)定的土污染風(fēng)險(xiǎn)篩選值，mg·kg-1；

（3）重金屬生物可利用指數(shù)（QBCF）計(jì)算公式如下：

式中，QBCF為評價區(qū)域農(nóng)產(chǎn)品的重金屬生物可利用風(fēng)險(xiǎn)商；A為研究區(qū)農(nóng)產(chǎn)品中重金屬的實(shí)測質(zhì)量分?jǐn)?shù)均值，mg·kg-1；C為研究區(qū)土壤中重金屬的實(shí)測質(zhì)量分?jǐn)?shù)均值，mg·kg-1；CS為土壤重金屬的限量值，采用GB 15618中規(guī)定的土壤污染風(fēng)險(xiǎn)篩選值，mg·kg-1；S為農(nóng)產(chǎn)品重金屬的限量標(biāo)準(zhǔn)值，采用 GB 2762中規(guī)定的食物中污染物限量標(biāo)準(zhǔn)，mg·kg-1；SE為研究農(nóng)產(chǎn)品中重金屬i的實(shí)測值的標(biāo)準(zhǔn)誤差。

1.3.4.3 風(fēng)險(xiǎn)評價

根據(jù)土壤污染風(fēng)險(xiǎn)等級，將耕地劃分為3個類別，將無污染的耕地劃為優(yōu)先保護(hù)類，低風(fēng)險(xiǎn)和中度風(fēng)險(xiǎn)的耕地劃為安全利用類，高風(fēng)險(xiǎn)和極高風(fēng)險(xiǎn)的耕地劃為嚴(yán)格管控類。稻田土壤重金屬污染風(fēng)險(xiǎn)等級見表2。根據(jù)現(xiàn)階段Cd污染治理技術(shù)水平，以0.5為間隔劃分農(nóng)產(chǎn)品風(fēng)險(xiǎn)等級。

表2 稻田土壤重金屬污染風(fēng)險(xiǎn)等級劃分Table 2 Risk classification of heavy metal pollution in paddy soil

1.4 圖形制作與數(shù)據(jù)處理

數(shù)據(jù)分析與處理于Microsoft Excel 2010、R 3.4.4和ArcGIS 10.2空間平臺進(jìn)行，其中空間插值借助反距離權(quán)重插值（IDW）進(jìn)行。反距離加權(quán)法（IDW）以插值點(diǎn)與樣本點(diǎn)間的距離為權(quán)重進(jìn)行加權(quán)平均，離插值點(diǎn)越近的樣本賦予的權(quán)重越大，插值結(jié)果介于插值數(shù)據(jù)的最大值和最小值之間，是一種簡單易行且常用的空間插值方法。該方法在應(yīng)用中直觀高效，在已知點(diǎn)分布均勻的情況下插值效果好，經(jīng)剔除異常值減弱甚至消除極值的影響后，可信度較高。對比樣條插值法（Spline）和克里金法（Kring）空間插值方法的預(yù)測結(jié)果后發(fā)現(xiàn)，反距離加權(quán)法（IDW）對本研究樣本預(yù)測誤差較小，精度較高。

2 結(jié)果與分析

2.1 樣品描述性統(tǒng)計(jì)

水稻產(chǎn)地土壤樣品數(shù)據(jù)統(tǒng)計(jì)結(jié)果顯示（表3），Cd總質(zhì)量分?jǐn)?shù)為0.005~0.937 mg·kg-1（均值為0.14 mg·kg-1)，Cd 含量超過土壤背景值（0.110 mg·kg-1），變異系數(shù)為85.7%，Cd的超標(biāo)樣點(diǎn)有252個，占總樣點(diǎn)數(shù)的16.1%，超標(biāo)率較高，表明水稻產(chǎn)地土壤明顯受到Cd的污染。稻米樣品數(shù)據(jù)統(tǒng)計(jì)結(jié)果顯示，Cd總質(zhì)量分?jǐn)?shù)為 0~0.936 mg·kg-1（均值為 0.163 mg·kg-1)，超標(biāo)點(diǎn)位為294個，超標(biāo)率為33.4%。

表3 珠三角稻米產(chǎn)地土壤和稻米Cd含量描述統(tǒng)計(jì)分析Table 3 Descriptive statistics of Cd contents in soil and rice from rice production areas in the Pearl River Delta mg·kg-1

2.2 地理單元劃分

綜合分析已有的資料得知：珠三角地區(qū)農(nóng)業(yè)用地主要分布在受河流沖積和三角洲沖積的平原地帶，地帶性土壤以紅壤和赤紅壤為主；其成土母質(zhì)是由西江、北江、東江等河流帶來的泥沙在海灣內(nèi)沉積而形成的；沉積物受河、海共同作用影響；水系由西、北江思賢滘以下、東江石龍以下的網(wǎng)河水系和注入三角洲的其他河流組成，注入三角洲的河流主要有潭江、流溪河、增江、沙河、高明河；涉重金屬行業(yè)企業(yè)主要分布于以廣州市為中心的周邊50 km以內(nèi)；珠三角地區(qū)西面和北面以羅平山脈為界，東側(cè)羅浮山區(qū)則是三角洲的東界（圖2和圖3）。

綜合珠三角地區(qū)行政區(qū)域、流域水系、地形、水文地質(zhì)以及成土母質(zhì)、土壤類型等要素的相對一致性，最終將風(fēng)險(xiǎn)評價區(qū)域劃分為 64個基礎(chǔ)地理單元（圖3）。

2.3 珠三角土壤Cd風(fēng)險(xiǎn)等級評價

2.3.1 土壤Cd風(fēng)險(xiǎn)等級評價

采用調(diào)查點(diǎn)位土壤樣品數(shù)據(jù)，借助反距離權(quán)重插值（IDW）法進(jìn)行空間插值，并對插值結(jié)果進(jìn)行土壤Cd風(fēng)險(xiǎn)等級評價，結(jié)果顯示珠三角地區(qū)土壤Cd有風(fēng)險(xiǎn)區(qū)域集中分布在北江流域兩岸，包括肇慶東南部、佛山中北部、中山市、廣州市南端及珠海市東北部（圖4）。

2.3.2 稻米Cd風(fēng)險(xiǎn)等級評價

圖4 基于點(diǎn)位數(shù)據(jù)的珠三角稻田土壤和稻米重金屬Cd污染風(fēng)險(xiǎn)等級空間分布Fig. 4 Spatial distribution of heavy metal Cd pollution risks in paddy soil and rice in the Pearl River Delta based on samplings

采用調(diào)查點(diǎn)位稻米樣品數(shù)據(jù)，借助反距離權(quán)重插值（IDW）法進(jìn)行空間插值，并對插值結(jié)果進(jìn)行稻米Cd風(fēng)險(xiǎn)等級評價，結(jié)果顯示珠三角地區(qū)稻米Cd有風(fēng)險(xiǎn)區(qū)域集中分布在北江流域兩岸，包括肇慶東部、佛山中北部、中山市東北部、廣州市南端、珠海市東北部和東莞東北部（圖4）。

2.3.3 土壤-稻米-富集系數(shù)協(xié)同Cd風(fēng)險(xiǎn)等級評價

采用調(diào)查點(diǎn)位土壤、稻米和富集系數(shù)數(shù)據(jù)，借助反距離權(quán)重插值（IDW）法進(jìn)行空間插值，并對插值結(jié)果進(jìn)行土壤-稻米-富集系數(shù) Cd協(xié)同風(fēng)險(xiǎn)等級評價，結(jié)果表明珠三角地區(qū)55.0%的稻田為無風(fēng)險(xiǎn)，低風(fēng)險(xiǎn)和潛在風(fēng)險(xiǎn)稻田的占比分別為19.0%和16.0%，中、高風(fēng)險(xiǎn)稻田的占比分別為 9.00%和1.00%（圖5）。

2.3.4 風(fēng)險(xiǎn)評價單元及等級

圖5 珠三角稻田土壤重金屬評價單元劃分及土壤重金屬Cd污染分級Fig. 5 Division of risk assessment units of heavy metals and classification of Cd pollution in paddy soils of the Pearl River Delta

根據(jù)上述基礎(chǔ)地理單元劃分、點(diǎn)位風(fēng)險(xiǎn)評價、水稻產(chǎn)地土壤風(fēng)險(xiǎn)評價、稻米風(fēng)險(xiǎn)評價、水稻富集系數(shù)及土壤-稻米協(xié)同風(fēng)險(xiǎn)評價結(jié)果劃分風(fēng)險(xiǎn)評價單元，將研究區(qū)域內(nèi)稻田劃分成125個風(fēng)險(xiǎn)評價單元，其中，優(yōu)先保護(hù)類、安全利用類和嚴(yán)格管控類的評價單元分別有85個、38個和2個。按風(fēng)險(xiǎn)等級劃分，區(qū)域內(nèi)無、低、中、高風(fēng)險(xiǎn)單元分別有85個、30個、8個和2個，未出現(xiàn)極高風(fēng)險(xiǎn)單元（圖5）。

3 討論

本研究首次建立了基于土壤、農(nóng)產(chǎn)品和富集系數(shù)相結(jié)合的土壤重金屬 Cd污染風(fēng)險(xiǎn)分類分級方法，它是基于土壤、農(nóng)產(chǎn)品以及富集系數(shù)的風(fēng)險(xiǎn)等級評價，評估參比值均選用了中國現(xiàn)行標(biāo)準(zhǔn)，而現(xiàn)行標(biāo)準(zhǔn)是根據(jù)生態(tài)毒理學(xué)數(shù)據(jù)、土壤理化性質(zhì)、土壤及農(nóng)產(chǎn)品種類等，基于不同暴露情景和暴露途徑下的暴露模型，參考相關(guān)毒理學(xué)模型模擬得出的生態(tài)閾值，因此本方法的評價基礎(chǔ)是科學(xué)可靠的。從劃分結(jié)果上看，評價過程綜合考慮了土壤和農(nóng)產(chǎn)品點(diǎn)位調(diào)查數(shù)據(jù)結(jié)果，因此該方法能夠有效識別珠三角稻田土壤Cd污染風(fēng)險(xiǎn)，針對性強(qiáng)；在進(jìn)行風(fēng)險(xiǎn)評價前劃分了地理評價單元，再結(jié)合地統(tǒng)計(jì)分析結(jié)果，綜合考慮了行政邊界、地形、河流水系等要素來劃分風(fēng)險(xiǎn)評價單元，有效地降低了土壤空間變異性對區(qū)域土壤污染風(fēng)險(xiǎn)等級評價的影響，從而在保證科學(xué)性的前提下，提高了方法的可操作性。風(fēng)險(xiǎn)評價過程嚴(yán)謹(jǐn)，風(fēng)險(xiǎn)分級思路明確，分級方法簡單易行。此外，本研究建立的稻田土壤重金屬污染分級方法對中國現(xiàn)有土壤污染防治相關(guān)政策及標(biāo)準(zhǔn)進(jìn)行了科學(xué)的完善和補(bǔ)充，適應(yīng)性強(qiáng)，易推廣應(yīng)用，適宜作為今后土壤分級分類管理的可靠的技術(shù)方法。

4 結(jié)論

土壤污染關(guān)系到區(qū)域可持續(xù)發(fā)展、國民健康和生態(tài)文明建設(shè)。中國的土壤污染防治修復(fù)工作，應(yīng)以經(jīng)濟(jì)與環(huán)境協(xié)調(diào)發(fā)展為指導(dǎo)思想，以全面落實(shí)“土十條”為主要工作目標(biāo)，提出土壤污染綜合分級防治新對策，建立健全區(qū)域土壤環(huán)境質(zhì)量綜合監(jiān)管與改善技術(shù)新體系。本研究基于風(fēng)險(xiǎn)管控的思路，通過構(gòu)建土壤、農(nóng)產(chǎn)品和富集系數(shù)相結(jié)合的稻田土壤重金屬污染分級方法，將珠三角典型重金屬污染物 Cd污染風(fēng)險(xiǎn)等級劃分為無風(fēng)險(xiǎn)、低風(fēng)險(xiǎn)、中度風(fēng)險(xiǎn)、高風(fēng)險(xiǎn) 4個等級。同時綜合分析縣級行政邊界、地形、河流水系、點(diǎn)位風(fēng)險(xiǎn)評價結(jié)果，將珠三角區(qū)域內(nèi)稻田土壤劃分為125個風(fēng)險(xiǎn)評價單元。本研究建立的稻田土壤重金屬污染分級方法可為中國土壤污染等級劃分提供有效、有力的參考，為維護(hù)中國土壤資源永續(xù)利用和生態(tài)安全，保護(hù)國民健康和實(shí)現(xiàn)區(qū)域經(jīng)濟(jì)社會可持續(xù)發(fā)展提供有力保障。