基于AHP、負權重和模糊數學的土壤質量評價研究

摘要：以北京市房山區為研究區域，選用土壤肥力和土壤環境質量的8個指標，采用層次分析法（AHP）、負權重和模糊數學對北京房山區土壤質量進行綜合評價，發現房山區土壤綜合指數在0.247～0.996之間，其中燕山地區土壤綜合指數最低（0.247），蒲洼鄉的土壤綜合指數最高（0.996）。根據耕地地力等級的劃分范圍，北京市房山區Ⅰ級占9%，Ⅱ級占23%，Ⅲ級占55%，Ⅳ級占14%，說明土壤質量總體處于中等水平。通過對比分析Ⅰ級和Ⅳ級指標含量，發現土壤綜合指數高的區域具有土壤肥力指標高、土壤環境質量指標低、沒有指標超標的特點，而土壤綜合指數低的區域具有土壤肥力指標低、土壤環境質量指標高、至少有一個指標超過臨界值的特點。

關鍵詞：土壤質量；負權重；層次分析法（AHP）；模糊數學

中圖分類號：F224.0 文獻標識碼：A 文章編號：0439-8114（2016）17-4480-04

DOI：10.14088/j.cnki.issn0439-8114.2016.17.031

Abstract：The research area located at Beijing Fangshan district，soil fertility and soil environment quality was taken as evaluating indicators which included soil organic matter， alkaline N，available P，available K，Cu，Pb，Cr，Cd. Comprehensive evaluation of soil quality was carried out by analytic hierarchy process，negative weights and fuzzy in Fangshan district of Beijing city，the results showed that the soil quality index was between 0.247 and 0.996，in which soil quality index of Yanshan area was the lowest of 0.247，soil quality index of Puwa town was the higest of 0.996. According to the classification cultivated land fertility range，Ⅰgrade accounted for 9%，Ⅱ grade accounted for 23%，Ⅲ grade 55%，Ⅳ grade 14%. Through comparative analysis it was found that areas of high soil quality index had high soil fertility，low soil environmental quality indicators， no indicators exceeded the threshold. Areas with low soil quality index had low soil fertility， soil environmental quality index was higher， and there was at least a characteristic index exceeded the threshold.

Key words： soil quality； negative weights； analytic hierarchy process（AHP）； fuzzy math

土壤是重要的農業資源和自然環境要素之一，提供了農作物生長的大部分養分，又是農作物生長發育的重要載體。近年來，工業污染、過量施肥、污水灌溉等原因導致土壤質量下降，土壤質量評價越來越成為土壤科學研究的熱點之一[1]。

土壤質量是土壤肥力質量、環境質量和健康質量的綜合量度[2]，而土壤質量評價則是科學準確衡量土壤維持生產力、環境凈化能力以及保障動植物健康能力的集中體現[3]。目前，對土壤質量的評價主要以土壤質量3個部分（土壤肥力質量、環境質量和健康質量）的單獨評價為主。由于這3個部分對土壤質量的影響不同，土壤肥力的指標大多對土壤質量產生正向積極影響，而土壤環境的指標通常對土壤產生消極的影響。但土壤質量評價是一個綜合質量評價，需要綜合考慮土壤質量的正負向指標。選擇土壤質量的土壤肥力指標和土壤環境指標進行綜合評價，有助于更加準確地了解土壤質量狀況。

在土壤質量評價中，土壤肥力質量評價采用的方法主要有層次分析法（AHP）和模糊數學[4-7]、支持向量機方法[8]、灰色關聯模型[9]、神經網絡[10]等，而土壤環境質量評價以單項污染指數和內梅羅綜合污染指數[11-15]為主。在眾多評價方法中，利用層次分析法和模糊數學評價土壤質量最為常見和成熟。本研究以北京市房山區為研究對象，將土壤肥力指標和土壤環境指標結合起來，采用層次分析法為評價指標賦權重、模糊數學計算指標隸屬度，對土壤環境質量進行綜合評價。

1 研究區概況及樣品的采集和分析

1.1 研究區概況

北京市房山區位于北京西南面，城區距市中心22 km，地處北緯39°30′～39°55′，東經115°25′～116°15′，面積201 900 hm2，東北與豐臺區相鄰，東與大興區一水相隔，南面和西面與河北涿州市、淶水縣相連，北與門頭溝以百花山為界。房山地形復雜多變，處于華北平原與太行山交界地帶，西部和北部是山地、丘陵，約占全區總面積2/3，人口數量83萬。境內有大小河流13條，拒馬河、大石河回旋曲折，永定河、小清河穿境而過。氣候屬于溫帶大陸性氣候，年平均氣溫為11.9 ℃，年平均降水量582.8 mm，年平均無霜期202 d。

1.2 樣品的采集與測定

土壤樣品采集時間為2006年3、5、8月，樣點布局考慮了土壤肥力指標和土壤環境指標的空間變異性、土地利用方式和土壤類型。采用梅花點法進行定點采樣，應用GPS確定樣點位置，采集的耕層為土壤的0～20 cm，采用多點混合和四分法制樣，共采集土壤樣品22個。采集的土壤樣品經風干后，采用四分法過20目、100目篩后風干保存備用。土壤有機質的測定采用重鉻酸鉀-外加熱氧化法，堿解氮的測定采用堿解擴散法，速效磷的測定采用 0.5 mol/L碳酸氫鈉浸提-鉬銻抗比色法，速效鉀的測定采用醋酸銨浸提-原子吸收法，Pb、Cd的測定為石墨爐原子吸收分光光度法，Cr和Cu的測定為火焰原子吸收分光光度法。本研究以鄉鎮作為評價單元，選取土壤肥力（有機質、堿解氮、速效磷、速效鉀）和重金屬污染元素（Cd、Pb、Cu、Cr）8項作為評價指標。

2 結果與分析

2.1 土壤質量綜合指數的計算

2.1.1 土壤質量評價指標的標準化處理 由于土壤質量的各個指標從不同角度反映土壤質量狀況，并且各指標實測值之間的量綱有很大區別，所以選用隸屬函數進行數據標準化。隸屬函數是根據一定范圍內評價指標與作物產量關系確定的，鑒于土壤環境質量指標與作物成負相關關系，選用降半梯形分布函數，而土壤肥力指標與農作物成正相關關系，選用S型分布函數。

降半梯形分布函數，分段函數界限值為相應的環境質量標準級數，函數如下：

土壤肥力指標與作物效應呈S型關系，所以隸屬函數也采用S型，為方便計算把曲線函數轉化為相應的折線函數，隸屬函數[16]如下：

結合研究區實際，參照國家土壤質量標準（GB15168-1995）和綠色食品中心發布的標準（NY/T391-2000），確定土壤肥力和土壤環境質量各指標隸屬度函數轉折點，如表1所示。

2.1.2 權重的設置 土壤質量評價中土壤肥力質量指標大多為正向因素，而土壤環境質量指標多為負向因素。確定正向因素和負向因素所占比例是土壤質量評價至關重要的一步。無論采用哪種方法，權重的加和必須為1，即正負權重之和也必須為1，對于土壤質量評價來說，土壤肥力直接作用于植物生長，理論上不應存在負值，而土壤污染危害植物生長，在土壤評價中發揮著負向作用。在類似的問題中，正向權重應大于負向權重，并且兩者差距較大。

陳彥清[17]利用X/Y的比值確定A值（X和Y的絕對值和）范圍，并詳細列出了X/Y在不同比值下A值的取值范圍表。根據研究區域的實際情況，本研究確定A=1.4，計算得X=1.2，Y=0.2，即正向權重總和為1.2，負向權重總和為-0.2。

2.1.3 確定評價指標的權重 目前，土壤質量評價指標權重確定的方法主要有主成分分析法、因子分析法、灰色關聯分析法、層次分析法等[18]，其中層次分析法把定性和定量分析結合起來，更符合客觀實際，應用較為普遍。根據土壤質量評價的特點，本研究采用層次分析法確定各項指標的權重。由于引入負權重，正向指標與負向指標在建立判斷矩陣和排序時相當于兩個獨立的體系，需分別建立判斷矩陣與排序，打分時不分正向影響與負向影響，遵循“最嚴重對應最高分”的規則，這樣在加入負權重計算時不利因素的“最嚴重”才可起到“減分最嚴重”的效果。

2.1.4 土壤質量綜合指數的計算 根據計算出的各指標的隸屬度和權重，計算出各土壤樣品的綜合評價指標值（Integrated quality index）采用下列公式計算：

IQI=∑WiNi （3）

式中，Wi，Ni分別為第i種評價指標的權重和隸屬度值。綜合指標值是每個土壤采樣點的土壤質量水平的得分，是土壤質量狀況的綜合反映。該指數值的大小反映了土壤質量的高低，IQI取值在0和1之間，該值越高，表明土壤質量越好，IQI值越低， 表明土壤質量越差[19]。根據國家農業部標準（NY/T33092-1996）全國耕地類型區耕地地力等級劃分并結合房山區實際情況，本研究采用等間距法將IQI劃分為5個等級，如表3所示。

2.2 土壤質量綜合指數分析

通過土壤綜合指數計算發現，房山區土壤綜合指數區間為0.247～0.996，其中燕山地區土壤綜合指數最低（0.247），蒲洼鄉的土壤綜合指數最高（0.996）。根據土壤質量等級劃分，土壤綜合質量屬于Ⅰ級的有蒲洼鄉（0.996）、霞云嶺鄉（0.991）；土壤綜合質量屬于Ⅱ級的有5個，分別為十渡鎮（0.778）、長陽鎮（0.692）、青龍湖鎮（0.649）、張坊鎮（0.641）、韓村河鎮（0.64）；土壤綜合質量屬于Ⅲ級的有12個，分別為閻村鎮（0.590）、大安山鎮0.586）、長溝鎮（0.582）、南窯鎮（0.571）、史家營鎮（0.551）、石樓鎮（0.550）、竇店鎮（0.514）、周口店鎮（0.508）、良鄉鎮（0.483）、房山鎮（0.456）、佛子莊鄉（0.453）、河北鎮（0.411），占55%；Ⅳ級有3個，分別為大石窩鎮（0.397）、琉璃河鎮（0.363）、燕山地區（0.247）。統計發現，土壤綜合質量屬于Ⅰ級的鄉鎮占9%，Ⅱ級占32%，Ⅲ級占55%，Ⅳ級占14%，Ⅰ級、Ⅱ級、Ⅲ級之和占到86%，說明房山區土壤質量總體處于中等水平。分析土壤綜合質量屬于Ⅳ級的鄉鎮，發現大石窩鎮、琉璃河鎮和燕山地區的土壤存在污染情況，主要是鎘超標，鎘的含量分別為0.16 mg/kg、0.14 mg/kg、0.23 mg/kg，超過清潔標準0.13 mg/kg，但大石窩鎮和琉璃河鎮的鎘含量未超過尚清潔的界限0.20 mg/kg，而燕山地區的鎘含量超過了尚清潔的界限，但未超出輕污染的界限0.60 mg/kg。總體來看，大石窩鎮和琉璃河鎮的土壤還屬于安全狀態，但燕山地區需要查清鎘來源，防治進一步污染。

3 小結與討論

土壤質量評價是為了能夠準確衡量土壤肥力和土壤環境的情況。本研究采用模糊數學、層次分析法和負權重進行土壤質量評價，充分考慮了土壤環境質量指標對土壤質量評價的限制作用。研究發現，房山區的蒲洼鄉和霞云嶺鄉的土壤綜合指數高，而大石窩鎮、琉璃河鎮和燕山地區土壤綜合指數低。通過對比分析發現，土壤綜合指數高的區域具有土壤肥力指標高、土壤環境質量指標低、沒有指標超標的特點，而土壤綜合指數低的區域具有土壤肥力指標低、土壤環境質量指標高、至少有1個指標超過臨界值的特點。研究結果與實際調查情況相同，也驗證了采用負權重進行土壤質量評價的可行性，同時也給土壤質量評價提供了一個新的思路。

本研究以北京市房山區為研究對象，選擇影響土壤質量最核心的兩個綜合因素土壤肥力和土壤環境質量作為評價指標，由于這兩個綜合指標對土壤質量影響不同，一般認為土壤肥力指標對土壤質量產生積極作用，而土壤環境質量指標對土壤質量起消極作用，屬于控制因素。本研究嘗試對正向因素和負向因素統一賦權重，統一評價，充分考慮負向因素對土壤質量的影響，有助于提高評價的準確性。

1）運用層次分析法和負權重對北京市房山區的土壤肥力指標和土壤環境指標權重進行確定，得出W（有機質、堿解氮、速效磷、速效鉀、鎘、銅、鉛、鉻）=（0.374、0.336、0.154、0.336、-0.103、-0.038、-0.038、-0.021）。

2）利用模糊數學的隸屬函數和層次分析法獲得的權重，計算土壤質量綜合指數，房山區土壤綜合指數區間為0.247～0.996，其中燕山地區土壤綜合指數最低（0.247），蒲洼鄉的土壤綜合指數最高（0.996）。

3）根據耕地地力等級的劃分范圍，北京市房山區Ⅰ級占9%，Ⅱ級占23%，Ⅲ級占55%，Ⅳ級占14%。

4）通過對比分析發現，土壤綜合指數高的區域具有土壤肥力指標高、土壤環境質量指標低、沒有指標超標的特點，而土壤綜合指數低的區域具有土壤肥力指標低、土壤環境質量指標高、至少有一個指標超過臨界值的特點。

參考文獻：

[1] 黃 勇，楊忠芳.土壤質量評價國外研究進展[J].地質通報，2009，28（1）：137-143.

[2] 王媛媛，楊忠芳，余 濤.土壤質量評價研究進展[J].安徽農業科學，2011，39（36）：22617-22622.

[3] 桂東偉，穆桂金，雷加強，等.干旱區農田不同利用強度下土壤質量評價[J].應用生態學報，2009，20（4）：894-900.

[4] 王鳳春，李曉秀，岳子明.基于AHP—Fuzzy與GIS的綠色食品產地發展潛力評價研究[J].首都師范大學學報（自然科學版），2008，29（6）：98-104.

[5] 李 梅，張學雷.基于GIS的農田土壤肥力評價及其與土體構型的關系[J].應用生態學報，2011，22（1）：129-136.

[6] 葉回春，張世文，黃元仿，等.北京延慶盆地農田表層土壤肥力評價及其空間變異[J].中國農業科學，2013，46（15）：3151-3160.

[7] 韓 平，王紀華，陸安祥，等.北京順義區土壤重金屬分布與環境質量評價[J].農業環境科學學報，2012，31（1）：106-112.

[8] 趙汝東，王殿武，陳延華，等.應用支持向量機方法對北京平原良田區土壤養分肥力的評價研究[J].土壤通報，2009，40（3）：514-517.

[9] 楊奇勇，楊勁松，姚榮江，等.基于GIS和改進灰色關聯模型的土壤肥力評價[J].農業工程學報，2010，26（3）：100-105.

[10] 趙艷玲，何廳廳，侯占東，等.基于BFA-BP神經網絡的土壤肥力評價[J].江蘇農業科學，2013，41（8）：340-344.

[11] 李 科，丁晴晴，傅 珊，等.太原市土壤重金屬污染空間分布及評價[J].環境化學，2015，34（4）：772-778.

[12] 陳 江，張海燕，何小峰，等.湖州市土壤重金屬元素分布及潛在生態風險評價[J].土壤，2010，42（4）：595-599.

[13] 陳玉東，王火焰，周健民，等.黑龍江省海倫市農田土壤重金屬分布特征及污染評級[J].土壤，2012，44（4）：613-620.

[14] 韓 平，王紀華，潘立剛，等.北京郊區田塊尺度土壤質量評價[J].農業工程學報，2009，25（S2）：228-234.

[15] 陳瀟霖，楊 丹，胡迪青，等.北京土壤重金屬分布及評價——以五環以內為例[J].環境科學與技術，2012，35（S2）：78-81.

[16] 趙振亞，姬寶霖，宋小園，等.基于層次分析和模糊數學法的公烏素土壤質量評價[J].干旱區研究，2014，31（6）：1010-1016.

[17] 陳彥清，楊建宇，嚴泰來，等.基于正負權重的農用地自然質量分計算方法研究[J].沈陽農業大學學報，2013，44（3）：268-272.

[18] 岳西杰，葛璽祖，王旭東.改進AHP模型在黃土高原溝壑區土壤質量評價中的應用——以陜西長武縣耕地為例[J].農業現代化研究，2011，32（2）：239-242.

[19] 王建國，楊林章，單艷紅.模糊數學在土壤質量評價中應用研究[J].土壤學報，2001，38（2）：176-183.