1、第第7 7講講n土壤質量指標土壤質量指標n土壤環境質量評價標準土壤環境質量評價標準n土壤環境質量評價方法土壤環境質量評價方法n污染權重的確定污染權重的確定n土壤環境質量綜合評價方法土壤環境質量綜合評價方法n土壤環境質量評價實例土壤環境質量評價實例土壤質量指標土壤質量指標 肥肥 力力 質質 量量 環環 境境 質質 量量 健健 康康 質質 量量 內內 涵涵 安全生產安全生產土土壤壤質質量量 保證最大生產力保證最大生產力 最佳生物學質量最佳生物學質量 不引發二次環境污染不引發二次環境污染n土壤質量是指土壤肥力質量、土壤環境質量及土壤土壤質量是指土壤肥力質量、土壤環境質量及土壤健康質量健康質量3 3方

2、面的綜合量度，既土壤在生態系統的范方面的綜合量度，既土壤在生態系統的范圍內，維持生物的生產能力、保護環境質量及促進圍內，維持生物的生產能力、保護環境質量及促進動植物健康的能力。動植物健康的能力。n土壤肥力質量是土壤提供植物養分和生產生物物質土壤肥力質量是土壤提供植物養分和生產生物物質能力的，是保證糧食生產的根本。能力的，是保證糧食生產的根本。n土壤環境質量是土壤容納、吸收和降解各種環境污土壤環境質量是土壤容納、吸收和降解各種環境污染物的能力。染物的能力。n土壤健康質量是土壤影響或促進人類和動植物健康土壤健康質量是土壤影響或促進人類和動植物健康的能力，是土壤環境質量在人類和動植物體上的反的能力，

3、是土壤環境質量在人類和動植物體上的反應。應。n描述性指標描述性指標n分析性指標分析性指標nHarrisHarris等提供了以等提供了以解釋框圖和訪談指南解釋框圖和訪談指南為基礎，包括通用調為基礎，包括通用調查表、特定地點調查表、相關報告卡組成的一套較為完整的查表、特定地點調查表、相關報告卡組成的一套較為完整的土壤質量評價信息收集工具土壤質量評價信息收集工具. .nRomigRomig等等 基于農民的土地評價方法中給出的基于農民的土地評價方法中給出的WisconsinWisconsin土壤土壤健康評分卡健康評分卡中，包括了中，包括了24 24 項土壤指標項土壤指標 14 14 項植物指標和項植物

4、指標和 3 3 項動物指標及項動物指標及 2 2 項水環境指標，根據農民對這些指標的評分項水環境指標，根據農民對這些指標的評分可以得到土壤的健康狀況可以得到土壤的健康狀況nUSDA-NRCSUSDA-NRCS設計的設計的 MarylandMaryland土壤質量評價手冊將土壤動物、土壤質量評價手冊將土壤動物、有機質顏色、根系和殘留物、表土緊實性、土壤耕性、侵蝕有機質顏色、根系和殘留物、表土緊實性、土壤耕性、侵蝕狀況、持水性、滲透性、作物長勢、狀況、持水性、滲透性、作物長勢、pH pH 狀況和保肥性分為差、狀況和保肥性分為差、中、好中、好 3 3 個等級，對每個項目各個等級的特征有詳細描述，個等

5、級，對每個項目各個等級的特征有詳細描述，在對這些項目等級評定的基礎上得到土壤質量的定性狀況在對這些項目等級評定的基礎上得到土壤質量的定性狀況 土壤質量評土壤質量評價中常用的價中常用的指標介紹指標介紹土壤環境質量標準土壤環境質量標準土壤環境質量標準土壤環境質量標準土壤環境質量標準 -GB 15618-1995-GB 15618-1995展覽會用地土壤環境質量評價標準（暫行）展覽會用地土壤環境質量評價標準（暫行） -HJ 3502007 HJ 3502007 溫室蔬菜產地環境質量評價標準溫室蔬菜產地環境質量評價標準 -HJ 333-HJ 333-2006 2006 食用農產品產地環境質量評價標準食

6、用農產品產地環境質量評價標準 -HJ 332-HJ 332-2006 2006 土壤環境質量標準值土壤環境質量標準值 mg/kgmg/kg級別級別 土壤土壤 pH值值一級一級 二級二級 三級三級 項目項目自然背景自然背景 7.5 6.5 鎘鎘 0.20 0.30 0.60 1.0 汞汞 0.15 0.30 0.50 1.0 1.5 砷砷 水田水田 旱地旱地 15 30 25 20 30 15 40 30 25 40 銅銅 農田等農田等 果園果園 35 50 100 100 400 150 200 200 400 鉛鉛 35 250 300 350 500 鉻鉻 水田水田 旱地旱地 90 250

7、 300 350 400 90 150 200 250 300 鋅鋅 100 200 250 300 500 鎳鎳 40 40 50 60 200 六六六六六六 0.05 0.50 1.0 滴滴娣滴滴娣 0.05 0.50 1.0 注：重金屬（鉻主要是三價）和砷均按元素量計，適用于陽離子交換量5cmol（+）/kg的土壤，若5cmol（+）/kg，其標準值為表內數值的半數。六六六為四種異構體總量，滴滴涕為四種衍生物總量。水旱輪作地的土壤環境質量標準，砷采用水田值，鉻采用旱地值。標準分級n類類為國家規定的自然保護區、集中式生活飲為國家規定的自然保護區、集中式生活飲用水源地和其他保護地區土壤標準用

8、水源地和其他保護地區土壤標準 , ,旨在保護旨在保護土壤環境質量基本保持自然背景水平土壤環境質量基本保持自然背景水平n類類為一般農田、蔬菜地、茶園、果園、牧場為一般農田、蔬菜地、茶園、果園、牧場等土壤等土壤 , ,土壤環境質量基本不對植物和環境造土壤環境質量基本不對植物和環境造成危害和污染成危害和污染n類類為林地土壤及污染物容量較大的高背景值為林地土壤及污染物容量較大的高背景值土壤和礦區附近的農田土壤土壤和礦區附近的農田土壤 （蔬菜地除外（蔬菜地除外) )標準制定的方法一級標準一級標準二級和三級二級和三級標準標準生態環境效應法生態環境效應法土壤背景值即地球土壤背景值即地球化學法化學法國外經驗n

9、內容體系內容體系n方法論上方法論上n其它其它保護土壤資源自身的“土壤自然質量指導值/標準”保護土壤生態和人體健康的“土壤環境和健康質量指導值/標準”類標準的制定上盡可能考慮不同土壤的母質和性質，基于區域內土壤的背景值、類標準的制定考慮人和生態受體對土壤污染物的取食攝入、皮膚接觸和呼吸攝入等引起的直接暴露風險環境介質：考慮了某些土壤污染物因侵蝕、地表徑流或淋溶進入地表水體或地下水體;揮發性有機污染物進入地下水或室內空氣、土壤污染物的遷移等帶來的潛在風險，并進行量化表征土地利用方式：考慮農業、居住區、工業和商業、飲用水水源地等標準存在的問題與建議l一級標準的制定過分強調一級標準的制定過分強調統一統

10、一l二級標準難以操作二級標準難以操作l標準中鉛的臨界含量定值標準中鉛的臨界含量定值偏高偏高l標準中有機污染物種類過標準中有機污染物種類過少少l重金屬形態選用單一重金屬形態選用單一l一類標準建議采用當地背一類標準建議采用當地背景值景值l二類標準可以按照土類分二類標準可以按照土類分類類l應降低土壤鉛的標準值應降低土壤鉛的標準值l增加標準中污染物的種類增加標準中污染物的種類l各標準值應增加有效態各標準值應增加有效態l加強土壤學研究加強土壤學研究l盡快制定地方標準盡快制定地方標準問題問題建議建議土壤環境質量評價方法土壤環境質量評價方法 (土壤環境監測技術規范土壤環境監測技術規范HJ/T 166 -20

11、04)HJ/T 166 -2004)評價方法n污染指數、超標率（倍數）評價法污染指數、超標率（倍數）評價法n內梅羅污染指數評價法內梅羅污染指數評價法n背景值及標準偏差評價法背景值及標準偏差評價法 n綜合污染指數法綜合污染指數法 污染指數、超標率（倍數）評價法土壤單項污染指數土壤單項污染指數= =土壤污染物實測值土壤污染物實測值/ /土壤污染物質土壤污染物質量標準量標準土壤污染累積指數土壤污染累積指數= =土壤污染物實測值土壤污染物實測值/ /污染物背景值污染物背景值土壤污染物分擔率（土壤污染物分擔率（% %）= =（土壤某項污染指數（土壤某項污染指數/ /各項污各項污染指數之和）染指數之和）1

12、00%100%土壤污染超標倍數土壤污染超標倍數= =（土壤某污染物實測值某污染物（土壤某污染物實測值某污染物質量標準）質量標準）/ /某污染物質量標準某污染物質量標準土壤污染樣本超標率（土壤污染樣本超標率（% %）= =（土壤樣本超標總數（土壤樣本超標總數/ /監測監測樣本總數）樣本總數）100%100%內梅羅污染指數評價法n內梅羅污染指數內梅羅污染指數（PNPN）= = （PIPI均）均）2 2+ + （PIPI最大最大2 2/2/21/21/2n污染等級分級污染等級分級 PN0.7PN0.7 清潔（安全）清潔（安全） 0.70.7PN1.0PN1.0 尚清潔（警戒限）尚清潔（警戒限） 1.

13、01.0PN2.0PN2.0 輕度污染輕度污染 2.02.0PN3.0PN3.0 中度污染中度污染 PNPN3.03.0 重污染重污染背景值及標準偏差評價法用區域土壤環境背景值（用區域土壤環境背景值（x x）95%95%置信度的范圍置信度的范圍（x x2s2s）來評價：）來評價：若土壤某元素監測值若土壤某元素監測值xIxIx x2s2s，則該元素缺乏或屬于，則該元素缺乏或屬于低背景土壤低背景土壤若土壤某元素監測值在若土壤某元素監測值在x x2s2s，則該元素含量正常，則該元素含量正常若土壤某元素監測值若土壤某元素監測值xIxIx x2s2s，則土壤已受該元素污，則土壤已受該元素污染，或屬于高背

14、景土壤染，或屬于高背景土壤綜合污染指數綜合污染指數法法綜合污染指數（綜合污染指數（CPICPI）包含了土壤元素背景值、土壤元素標準尺）包含了土壤元素背景值、土壤元素標準尺度因素和價態效應綜合影響。其表達式：度因素和價態效應綜合影響。其表達式： CPI=X(1+REP)+YDDMB/(ZDDSB)CPI=X(1+REP)+YDDMB/(ZDDSB) X X、Y Y分別為測量值超過標準值和背景值的數目，分別為測量值超過標準值和背景值的數目，RPE-RPE-相對污相對污染當量，染當量，DDMB-DDMB-元素測定濃度偏離背景值的程度，元素測定濃度偏離背景值的程度，DDSB-DDSB-土壤標準土壤標準

15、偏離背景值的程度，偏離背景值的程度，Z Z為用作標準元素的數目為用作標準元素的數目 X X Y Y CPI CPI 評價評價0 0 0 0 0 0 背景狀態背景狀態0 0 11 0 0CPICPI1 1 未污染狀態，表未污染狀態，表示偏離背景值相對程度示偏離背景值相對程度1 1 1 1 11 污染狀態，數值越大表示污染程度相對污染狀態，數值越大表示污染程度相對 越嚴重越嚴重污染權重的確定污染權重的確定賦權方法賦權方法主觀賦權法主觀賦權法客觀賦權法客觀賦權法由專家根據經驗進行由專家根據經驗進行主觀判斷而得到權數，主觀判斷而得到權數，然后再對指標進行綜然后再對指標進行綜合評價。合評價。層次分析法層

16、次分析法功效系數法功效系數法模糊評價法模糊評價法專家調查法專家調查法根據指標之間的相關根據指標之間的相關關系或各項指標的變關系或各項指標的變異系數來確定權數進異系數來確定權數進行綜合評價行綜合評價因子分析法因子分析法主成分分析法主成分分析法熵值法聚類分析熵值法聚類分析法法變異系數法灰色變異系數法灰色關聯分析法、關聯分析法、TOPSISTOPSIS法法神經網絡分析法神經網絡分析法判別分析法判別分析法 污染權重的確定方法污染權重的確定方法污染權重的確定方法權重確定方法權重確定方法以專家咨詢值為判以專家咨詢值為判定依據的賦權方法定依據的賦權方法以因子標準值為判定依以因子標準值為判定依據的賦權方法據的

17、賦權方法以因子實測值以因子實測值為判定依據的為判定依據的賦權方法賦權方法以因子實測值與標準值以因子實測值與標準值為雙重判定依據的賦權為雙重判定依據的賦權方法方法以專家咨詢值為判定依據的賦權方法DelphiDelphi法法德爾菲法德爾菲法(Delphi(Delphi Method)Method)系以一系列系以一系列問卷向各類專家征詢意見，依據所有專家對原問問卷向各類專家征詢意見，依據所有專家對原問卷的答復再擬定下一份問卷，再次向各類專家征卷的答復再擬定下一份問卷，再次向各類專家征詢意見，直到大多數專家的意見看法趨于一致才詢意見，直到大多數專家的意見看法趨于一致才獲至結論，獲至結論，此法通常較適合

18、作長期預測此法通常較適合作長期預測。層次分析法（層次分析法（AHPAHP）是將決策總是有關的元素分是將決策總是有關的元素分解成目標、準則、方案等層次，在此基礎之上進解成目標、準則、方案等層次，在此基礎之上進行定性和定量分析的決策方法。行定性和定量分析的決策方法。這種方法的特點是在對復雜的決策問題的本質、影這種方法的特點是在對復雜的決策問題的本質、影響因素及其內在關系等進行深入分析的基礎上，響因素及其內在關系等進行深入分析的基礎上，利用較少的定量信息使決策的思維過程數學化，利用較少的定量信息使決策的思維過程數學化，從而為多目標、多準則或無結構特性的復雜決策從而為多目標、多準則或無結構特性的復雜決

19、策問題提供簡便的決策方法。問題提供簡便的決策方法。尤其適合于對決策結尤其適合于對決策結果難于直接準確計量的場合。果難于直接準確計量的場合。灰色關聯法灰色關聯法是對于兩個系統之間的因素，其是對于兩個系統之間的因素，其隨時間或不同對象而變化的關聯性大小的量度，隨時間或不同對象而變化的關聯性大小的量度，稱為關聯度。在系統發展過程中，若兩個因素稱為關聯度。在系統發展過程中，若兩個因素變化的趨勢具有一致性，即同步變化程度較高，變化的趨勢具有一致性，即同步變化程度較高，即可謂二者關聯程度較高；反之，則較低。因即可謂二者關聯程度較高；反之，則較低。因此，灰色關聯分析方法，是根據因素之間發展此，灰色關聯分析方

20、法，是根據因素之間發展趨勢的相似或相異程度，亦即趨勢的相似或相異程度，亦即“灰色關聯度灰色關聯度”，作為衡量因素間關聯程度的一種方法。因此，作為衡量因素間關聯程度的一種方法。因此，灰色關聯度分析對于一個系統發展變化態勢提灰色關聯度分析對于一個系統發展變化態勢提供了量化的度量，非常適合動態歷程分析。供了量化的度量，非常適合動態歷程分析。 以因子標準值為判定依據的賦權方法n閾值賦權法閾值賦權法 此法是根據環境質量標準中因子在各級別間的平均相此法是根據環境質量標準中因子在各級別間的平均相對差值的大小來賦權的。標準差值越大，說明指標危對差值的大小來賦權的。標準差值越大，說明指標危害性較小，故權值也應較

21、小。害性較小，故權值也應較小。 fifi= 1/(m-1)= 1/(m-1) (S(Si,j+1i,j+1-S-Sijij) WiWi=(1/Si)/ =(1/Si)/ (1/Si)(1/Si)n簡單賦權法簡單賦權法 以某因子的標準值為以某因子的標準值為SiSi, ,則相應的權重為：則相應的權重為： WiWi=(1/Si)/ =(1/Si)/ (1/Si)(1/Si)m-1j=1ni=1i=1n以因子實測值為判定依據的賦權方法熵賦權法熵賦權法在信息論中，熵是對不確定性的一種度量。信在信息論中，熵是對不確定性的一種度量。信息量越大，不確定性就越小，熵也就越小；信息量越小，不息量越大，不確定性就越

22、小，熵也就越小；信息量越小，不確定性越大，熵也越大。根據熵的特性，可以通過計算熵值確定性越大，熵也越大。根據熵的特性，可以通過計算熵值來判斷一個事件的隨機性及無序程度，也可以用熵值來判斷來判斷一個事件的隨機性及無序程度，也可以用熵值來判斷某個指標的離散程度，指標的離散程度越大，該指標對綜合某個指標的離散程度，指標的離散程度越大，該指標對綜合評價的影響越大。評價的影響越大。主分量分析法主分量分析法主成分分析也稱主分量分析，旨在利用降主成分分析也稱主分量分析，旨在利用降維的思想，把多指標轉化為少數幾個綜合指標。在統計學中，維的思想，把多指標轉化為少數幾個綜合指標。在統計學中，主成分分析（主成分分析

23、（principal components principal components analysis,PCAanalysis,PCA）是一種）是一種簡化數據集的技術。它是一個線性變換。這個變換把數據變簡化數據集的技術。它是一個線性變換。這個變換把數據變換到一個新的坐標系統中，使得任何數據投影的第一大方差換到一個新的坐標系統中，使得任何數據投影的第一大方差在第一個坐標在第一個坐標( (稱為第一主成分稱為第一主成分) )上，第二大方差在第二個坐上，第二大方差在第二個坐標標( (第二主成分第二主成分) )上，依次類推。主成分分析經常用減少數據上，依次類推。主成分分析經常用減少數據集的維數，同時保持

24、數據集的對方差貢獻最大的特征集的維數，同時保持數據集的對方差貢獻最大的特征. .以因子實測值與標準值為雙重判定依據的賦權方法n污染貢獻率法污染貢獻率法 I Ii i= =Ci/SiCi/Si WiWi=I=Ii i/ / I Ii i(i(i=1,2n)=1,2n)n污染分擔率法污染分擔率法 I Iijij=(=(C Cijij-Bi)/Si-Bi)/Si uiui= = I Iijij WiWi= =u ui i/ / u ui in統計概率法統計概率法 假定污染因子的實測值在時間和空間假定污染因子的實測值在時間和空間范圍內的變化近似正態分布，然后按統計范圍內的變化近似正態分布，然后按統計概

25、率求的權重概率求的權重 mj=1ni=1土壤環境質量綜合評價方法土壤環境質量綜合評價方法模糊數學的方法模糊數學的方法模糊綜合評價的優點土壤是一個復雜的多相有機土壤是一個復雜的多相有機無機復合體無機復合體 , ,涉及涉及到各種因素的相互作用到各種因素的相互作用 , , 而且評價中存在大量的而且評價中存在大量的模糊現象和模糊概念。目前廣泛使用的綜合指數模糊現象和模糊概念。目前廣泛使用的綜合指數評價方法將這種模糊性用一種比較明確的界限加評價方法將這種模糊性用一種比較明確的界限加以區分以區分 , , 結論與實際情況存在一定差異。模糊綜結論與實際情況存在一定差異。模糊綜合評價法評價土壤環境質量狀況合評價

26、法評價土壤環境質量狀況 , , 可以較好地克可以較好地克服這種不科學性服這種不科學性 , , 有效地解決評價標準邊界模糊有效地解決評價標準邊界模糊和監測誤差對評價結果的影響。和監測誤差對評價結果的影響。模糊綜合評價的方法原理模糊綜合評價的方法原理建立因素集建立因素集U = u U = u 1 1, u , u 2 2, , u, , un n , , 即有即有 n n 個評價指標個評價指標 確定評價確定評價（評語）集（評語）集V= v v1 1 , v , v 2 2, , , ,v vn n ,即代表評價等級、分類的集合即代表評價等級、分類的集合 建立隸屬函數，構造模糊關系矩陣。建立隸屬函數

27、，構造模糊關系矩陣。確定模糊確定模糊權重權重向量向量 A = aA = a1 1 , a , a 2 2, , a, , an n , , 模糊復合運算模糊復合運算模糊數學模型突出單個決定因素作突出單個決定因素作用的模型用的模型體現各個參評因素綜體現各個參評因素綜合作用的模型合作用的模型單因素決定模型單因素決定模型主導因素模型主導因素模型加權平均模型加權平均模型幾何平均模型幾何平均模型權重的計算權重的計算Wi為第為第i因子的權重因子的權重Si為為W= 0.134,0.141，0.152，0.118，0.040，0.109，0.131，0.175單因素決定模型單因素決定模型加權平均模型兩種模型的

28、出發點不同n單因素決定模型單因素決定模型 只考慮了最突出的因素只考慮了最突出的因素 , ,其它因素的作用被其它因素的作用被弱化，主要是體現了相對含量較高、污染狀況弱化，主要是體現了相對含量較高、污染狀況較嚴重的評價因子對整個評價結果的影響較嚴重的評價因子對整個評價結果的影響n加權平均模型加權平均模型 充分考慮了每個因素對綜合評價的貢獻充分考慮了每個因素對綜合評價的貢獻 , ,并并把貢獻按權重進行分配把貢獻按權重進行分配 , ,其評價結果是各個參其評價結果是各個參評指標綜合作用的產物評指標綜合作用的產物對模糊綜合評價的改進將污染物濃度和毒性級別指數加權疊加將污染物濃度和毒性級別指數加權疊加, ,

29、作歸一作歸一化處理化處理 , ,得到某污染組分的權重公式得到某污染組分的權重公式 土壤環境質量評價實例土壤環境質量評價實例杭州市郊蔬菜基地土壤和蔬菜中杭州市郊蔬菜基地土壤和蔬菜中 PbPb、Zn Zn 和和 CuCu含量的環境質量評價含量的環境質量評價研究區概況杭州市屬亞熱帶季風性氣候杭州市屬亞熱帶季風性氣候, ,四季分明四季分明 , ,溫暖濕潤溫暖濕潤 , ,光照充足光照充足 雨量充沛雨量充沛. .杭州市區年平均氣溫為杭州市區年平均氣溫為16.4 16.4 ,年光照時數為年光照時數為1752h,1752h,無霜期無霜期 235d 235d 左右左右 , ,年平年平均降雨量均降雨量1439mm

30、,1439mm,年相對濕度為年相對濕度為77%,77%,全年各季水、全年各季水、熱差異大熱差異大. .杭州市地處錢塘江水系入海口杭州市地處錢塘江水系入海口 , ,位于水位于水網平原與濱海平原之間網平原與濱海平原之間 , ,地勢平坦地勢平坦 , ,海拔海拔 3 35m,5m,境內地貌類別多樣境內地貌類別多樣 , ,江河縱橫江河縱橫 , ,湖泊密布湖泊密布 , ,成土母成土母質主要為河湖相或淺海相沉積物質主要為河湖相或淺海相沉積物 . .本研究主要涉及本研究主要涉及杭州市的江干、下城、拱墅、西湖杭州市的江干、下城、拱墅、西湖 4 4 個區個區 , ,是杭是杭州市的主要蔬菜基地。州市的主要蔬菜基地。

31、土壤樣品的采集土壤樣品的分析土壤中土壤中Pb,ZnPb,Zn，CuCu土壤中水溶土壤中水溶性性Pb,ZnPb,Zn，CuCu用用HF-HClO4-HNO3 消化，原消化，原子吸收分光光度子吸收分光光度法法用用 15mL 去離子水在室溫下去離子水在室溫下振蕩振蕩 2h 后后 ,離心離心 30min, 分分離定容后上清液離定容后上清液 ,用原子吸收用原子吸收分光光度計測定分光光度計測定評價方法n單因子指數法單因子指數法n內梅羅綜合污染指數法內梅羅綜合污染指數法評價標準參考文獻參考文獻1、謝正苗等、謝正苗等.杭州市郊蔬菜基地土壤和蔬菜中杭州市郊蔬菜基地土壤和蔬菜中 PbPb、Zn Zn 和和 CuC

32、u含量的環境質量評價（含量的環境質量評價（J J）. .環環境科學境科學.2006,27.2006,27（4 4）：）：742-747.742-747.2 2、郭巖等、郭巖等. .汕頭市典型區域土壤重金屬污染特征及評價（汕頭市典型區域土壤重金屬污染特征及評價（J J）. .環境科學環境科學.2007,28(5).2007,28(5)：1067-1074.1067-1074.3 3、朱青等、朱青等. .兩種模糊數學模型在土壤重金屬綜合污染評價中的應用與比較兩種模糊數學模型在土壤重金屬綜合污染評價中的應用與比較（J J）. .環境評環境評價價.2004.2004：3030（123123）：）：53

33、-57.53-57.4 4、石曉翠等、石曉翠等. .模糊數學模型在土壤重金屬污染評價中的應用模糊數學模型在土壤重金屬污染評價中的應用（J J）. .土壤通報土壤通報.2006.2006，37(2)37(2)：332-337.332-337.5 5、竇磊等、竇磊等. .針對土壤重金屬污染評價的模糊數學模型的改進及應用（針對土壤重金屬污染評價的模糊數學模型的改進及應用（J J）. .土壤通土壤通報報.2007,38(1):101-104.2007,38(1):101-104.6 6、武偉等、武偉等. .土壤養分的模糊綜合評價土壤養分的模糊綜合評價（J J）. .西南農業大學學報西南農業大學學報.2000,22.2000,22（3 3）:270-272.:270-272.7 7、劉曉冰等、劉曉冰等. .土壤質量及其評價指標（土壤質量及其評價指標（J J）. .農業系統科學與綜合研究農業系統科學與綜合研究.2002.2002；1818（2 2）:109-112.:109-112.8 8、張華，張甘霖、張華，張甘霖. .土壤質量指標和評價方法（土壤質量指標和評價方法（J J）. .土壤土壤.2001,(6):326-330.2001,(6):32