李小慶

摘 要：該研究采用模糊綜合評(píng)判法，對(duì)山南市乃東縣、貢嘎縣、扎朗縣、隆子縣、桑日縣、瓊結(jié)6個(gè)主要的農(nóng)業(yè)區(qū)縣856個(gè)土壤樣品進(jìn)行了土壤養(yǎng)分等級(jí)評(píng)價(jià)。結(jié)果表明：山南市農(nóng)耕地土壤的整體養(yǎng)分等級(jí)為中等水平，其中中等偏上的占23.01%、中等水平占47.08%、中等偏下的占21.96%。山南市農(nóng)耕地土壤pH整體偏中性，有機(jī)質(zhì)、全氮含量適宜；堿解氮、全磷、有效磷、全鉀和速效鉀含量較低；有效硼含量豐富。影響土壤肥力的主要因素是有機(jī)質(zhì)、全磷、速效鉀、全氮、堿解氮、有效磷，各主要養(yǎng)分元素對(duì)土壤肥力的貢獻(xiàn)率為有機(jī)質(zhì)>全磷>速效鉀>全氮>堿解氮>有效磷。

關(guān)鍵詞：土壤肥力；土壤養(yǎng)分綜合指數(shù)；模糊綜合評(píng)判

中圖分類號(hào) S158.2 文獻(xiàn)標(biāo)識(shí)碼 A 文章編號(hào) 1007-7731（2017）12-0054-04

Astract：The nutrient content of 856 soils sampled from Naidong，Kongga，Chanang，Longzi，Sangri and Qiongjie counties of Shannan were measured.and the corresponding soil nutrient status was evaluated by fussy synthetic estimation.The results showed that the soil fertility in Shanan was at secondary level，among which 23.01% was better，14.08% was middling and 21.96% was worse.Most of the soil pHs were neutral.The content of available bon was at rich level and the contents of organic matter and total nitrogen were at intermediate level.Unfortunately，the contents of available nitrogen，total phosphorus，available phosphorus and available potassium were low.The main factors that decide the soil fertility were organic matter，total phosphorus，available phosphorus，total nitrogen，available nitrogen and available potassium.And the contribution rates of main factors were organic matter>total phosphorus>available potassium>total nitrogen>available nitrogen>available phosphorus.

Key words：Soil fertility；Soil nutrient synthesis index；Fuzzy synthetic evaluation

土壤是人們賴以生存和發(fā)展的物質(zhì)基礎(chǔ)，是一切物質(zhì)生產(chǎn)最基本的源泉，是農(nóng)業(yè)生產(chǎn)最重要的資源。土地資源對(duì)農(nóng)業(yè)生產(chǎn)的發(fā)展，對(duì)人們物質(zhì)生活水平的提高以及對(duì)國(guó)家經(jīng)濟(jì)的發(fā)展都起著十分重要的作用。了解土地資源的質(zhì)量狀況，可以為保持和培肥土壤、合理利用土壤提供依據(jù)。在土壤質(zhì)量的確定方面，土壤養(yǎng)分是其中的重要指標(biāo)之一[1]。影響土壤肥力質(zhì)量的土壤屬性很多，土壤的每個(gè)理化性狀、生物屬性都可能是最終的評(píng)價(jià)指標(biāo)，但對(duì)于特定區(qū)域，由于考慮到土壤屬性的時(shí)空變異性、數(shù)據(jù)獲取的成本及因子間的相關(guān)性等因素，顯然不可能獲取所有因子的數(shù)據(jù)。參評(píng)因子的選定與分級(jí)指標(biāo)的劃分是土壤養(yǎng)分評(píng)價(jià)的基礎(chǔ)和核心工作，根據(jù)農(nóng)業(yè)研究者，在使用頻率最高的且具有穩(wěn)定性的評(píng)價(jià)耕作土壤肥力的因子中選取pH、有機(jī)質(zhì)、全氮、堿解氮、全磷、有效磷、全鉀、速效鉀、有效硼共9項(xiàng)作為本次土壤肥力評(píng)價(jià)指標(biāo)[9]。為此，本研究對(duì)各評(píng)價(jià)指標(biāo)進(jìn)行了分析和測(cè)定，應(yīng)用相關(guān)分析和模糊數(shù)學(xué)的原理，采用模糊綜合評(píng)價(jià)模型，對(duì)各評(píng)價(jià)指標(biāo)進(jìn)行綜合評(píng)價(jià)，確定土壤肥力等級(jí)，為山南市主要農(nóng)耕地的土壤養(yǎng)分診斷和平衡施肥提供科學(xué)依據(jù)。

1 材料與方法

1.1 研究區(qū)概況 研究區(qū)位于崗底斯山脈和喜瑪拉雅山脈之間的河谷地帶，雅魯藏布江中游，北緯26°51′10″～29°49′27″，東經(jīng)90°04′04″～94°21′40″。區(qū)域平均海拔3700m左右，總國(guó)土面積7.97萬(wàn)km2，邊境線長(zhǎng)934km，戰(zhàn)略位置十分重要。研究區(qū)具有豐富的氣候資源，有經(jīng)向和緯向的地帶性變化，根據(jù)地帶性差異，自南而北全地區(qū)分為山地?zé)釒А⑸降貋啛釒Ш透咴瓬貛?個(gè)氣候帶。全年干濕季分明，多夜雨，年降水量大約在200～500mm，降水多集中在6—9月份，占全年降水量70%～85%。太陽(yáng)輻射強(qiáng)，年平均太陽(yáng)輻射量為7322MJ/m2，年平均日照時(shí)數(shù)為2600～3300h，年日照百分率66%，氣溫年差較小，日差較大，全地區(qū)年平均氣溫5.6℃左右，極端最高氣溫31℃，極端最低氣溫-27℃。冬季降水少，氣候干燥，多大風(fēng)，年平均大風(fēng)日數(shù)為70d左右（瞬時(shí)風(fēng)速達(dá)到17m/s以上成為大風(fēng)），風(fēng)期主要分布在12月至翌年3月，全年無(wú)霜期120d左右。研究區(qū)盛產(chǎn)青稞、小麥、豌豆、蠶豆、玉米、蕎麥、胡豆、油菜等作物，是西藏糧油主要產(chǎn)地之一。

1.2 土壤樣品采集 土壤樣品采自山南市乃東、貢嘎、扎朗、隆子、桑日和瓊結(jié)6個(gè)主要農(nóng)業(yè)區(qū)縣，也是測(cè)土配方區(qū)縣。土樣采集在9月15日至10月15日進(jìn)行。采樣工具用不銹鋼土鉆。采樣深度為0～20㎝耕層，田塊采樣路線用“S”法，均勻、等量、隨機(jī)采集15～20個(gè)土樣，充分混合后，四分法留取1kg土樣。剔去石礫和植物殘?jiān)冉塍w，在室內(nèi)自然風(fēng)干后磨碎，分別過(guò)2.00mm和0.25mm的尼龍篩，裝入塑料袋中，存放在陰涼干燥的地方備用。

1.3 土樣樣品化驗(yàn)分析 土壤樣品各理化性質(zhì)均按《測(cè)土配方施肥土壤分析技術(shù)規(guī)范》操作步驟完成，各分析方法見(jiàn)表1。

2 結(jié)果與分析

2.1 研究區(qū)土壤各養(yǎng)分因子含量水平分析 土壤各養(yǎng)分因子測(cè)定值見(jiàn)表2，由表2可知，研究區(qū)土壤pH值處在5.5～8.8，土壤中性居多。土壤有機(jī)質(zhì)、全氮、堿解氮、全磷、有效磷、全鉀、速效鉀含量均值分別為23.2g/kg、2.09g/kg、110mg/kg、0.551g/kg、10.5mg/kg、11.9g/kg、78mg/kg，均處在缺乏和潛在缺乏狀態(tài)，有效硼含量均值為1.89mg/kg，含量比較豐富。從變異系數(shù)看：pH、全氮、堿解氮、全磷變異系數(shù)相對(duì)較低，有機(jī)質(zhì)、有效磷、全鉀、速效鉀、有效硼變異系數(shù)都在50%以上，說(shuō)明這部分養(yǎng)分指標(biāo)在研究區(qū)的分布很不均衡。

2.2 土壤養(yǎng)分的綜合評(píng)價(jià)

2.2.1 土壤養(yǎng)分等級(jí)評(píng)價(jià)模型的建立 進(jìn)行土壤肥力評(píng)價(jià)首先要對(duì)各評(píng)價(jià)因子的優(yōu)劣進(jìn)行分析研究，目前農(nóng)業(yè)工作者建立了許多評(píng)價(jià)土壤肥力的方法和模型，列如田間描述評(píng)價(jià)法、指數(shù)法評(píng)價(jià)模型、環(huán)境指數(shù)評(píng)價(jià)模型、模糊綜合評(píng)價(jià)模型地統(tǒng)計(jì)學(xué)評(píng)價(jià)、系統(tǒng)評(píng)價(jià)方法、動(dòng)力學(xué)評(píng)價(jià)、決策樹(shù)方法等[3，9，14-15]。而從土壤養(yǎng)分狀況出發(fā)進(jìn)行土壤肥力評(píng)價(jià)，應(yīng)用最多的還是模糊綜合評(píng)價(jià)模型。首先確定隸屬度函數(shù)計(jì)算隸屬度值Ni，隸屬度函數(shù)實(shí)際上是所要評(píng)價(jià)的肥力指標(biāo)與作物生長(zhǎng)效應(yīng)曲線之間的數(shù)學(xué)表達(dá)式，它可以將肥力評(píng)價(jià)指標(biāo)標(biāo)準(zhǔn)化，轉(zhuǎn)變?yōu)闊o(wú)量綱值的隸屬度。由于土壤因子變化具有連續(xù)性，故各評(píng)價(jià)指標(biāo)采用連續(xù)性質(zhì)的隸屬度函數(shù)。常用的隸屬度函數(shù)類型有兩類，即拋物線型隸屬函數(shù)和S型隸屬函數(shù)[3]。確定各肥力指標(biāo)適用的函數(shù)類型及函數(shù)上限（X2）、下限（X1）等參數(shù)，將各指標(biāo)實(shí)測(cè)數(shù)據(jù)代入函數(shù)公式得到隸屬度值。如作物的效應(yīng)隨著養(yǎng)分值在一定范圍內(nèi)增加而增加，達(dá)到最高點(diǎn)后又隨養(yǎng)分值得增加而下降，呈現(xiàn)拋物線型的，如土壤pH，其函數(shù)公式見(jiàn)公式）[4，10]（2）；如作物的效應(yīng)隨著養(yǎng)分值的增加而增加，達(dá)到最高值后，隨著養(yǎng)分因子的增加，作物效應(yīng)值則穩(wěn)定在一定水平上，呈現(xiàn)為S型隸屬函數(shù)，如有機(jī)質(zhì)、氮、磷、鉀等含量，相應(yīng)的隸屬函數(shù)公式見(jiàn)公式[4，10]（3）。參考相關(guān)文獻(xiàn)和專家建議[7-8]，根據(jù)研究區(qū)土壤pH統(tǒng)計(jì)分析結(jié)果，研究區(qū)土壤pH下限為5.5，上限為8.5，最優(yōu)值為6.5～7.5。

各單項(xiàng)評(píng)價(jià)指標(biāo)豐缺標(biāo)準(zhǔn)因不同的作物和土壤而有差異，因此，依據(jù)本文研究區(qū)土壤的理化特性，已有的研究結(jié)果表明，不同土壤肥力指標(biāo)的豐缺狀況因作物類型而異。筆者結(jié)合生產(chǎn)實(shí)踐經(jīng)驗(yàn)，并匯總專家建議，確定了各肥力指標(biāo)適合的隸屬度函數(shù)類型和相應(yīng)取值[7，8]（表3）。將各指標(biāo)實(shí)測(cè)值代入相應(yīng)表達(dá)式，便可得到各指標(biāo)的隸屬度值。

建好隸屬度函數(shù)后，用Excel軟件計(jì)算隸屬度值，隸屬度值大小在0.1～1.0，最大值1.0表示土壤肥力指標(biāo)完全適宜作物生長(zhǎng)，最低值0.1表示土壤肥力指標(biāo)嚴(yán)重缺乏。最小值0.1表示作物生長(zhǎng)所需求的土壤屬性作用分值的最低閾值，結(jié)合生產(chǎn)實(shí)際將最小值定為0.1而非0[3，7]。

2.2.2 單項(xiàng)指標(biāo)權(quán)重系數(shù)的確定 權(quán)重是指各指標(biāo)對(duì)土壤肥力的影響程度或貢獻(xiàn)率，表示各指標(biāo)在土壤肥力中作用和地位的不同，對(duì)肥力質(zhì)量的影響程度也就不一樣，如何確定各項(xiàng)肥力指標(biāo)的權(quán)重是土壤肥力綜合評(píng)價(jià)的一個(gè)關(guān)鍵問(wèn)題和難點(diǎn)問(wèn)題[4]。以往的評(píng)價(jià)大多由專家打分方法確定權(quán)重系數(shù)，主觀性比較強(qiáng)。除專家經(jīng)驗(yàn)打分法外，當(dāng)今一些數(shù)學(xué)和統(tǒng)計(jì)學(xué)的方法逐漸用于確定權(quán)重系數(shù)，諸如主成分分析法，層次分析法，等值法，等差法，回歸分析法，相關(guān)系數(shù)分析法[3，13]等。筆者采用相關(guān)系數(shù)法，相關(guān)系數(shù)法采用各指標(biāo)間相關(guān)系數(shù)來(lái)確定權(quán)重系數(shù)，應(yīng)用統(tǒng)計(jì)軟件SPSS 22求出各養(yǎng)分指標(biāo)的相關(guān)系數(shù)，以某項(xiàng)養(yǎng)分指標(biāo)與其它養(yǎng)分指標(biāo)之間相關(guān)系數(shù)的平均值占所有養(yǎng)分指標(biāo)相關(guān)系數(shù)平均值總和的比值，作為單項(xiàng)養(yǎng)分指標(biāo)在表征土壤養(yǎng)分中的貢獻(xiàn)[7]。各土壤養(yǎng)分因子的相關(guān)系數(shù)見(jiàn)表4。從表4可看出，pH與有機(jī)質(zhì)、全氮、全磷、全鉀、速效鉀和有效硼均達(dá)到極顯著的正相關(guān)關(guān)系。有機(jī)質(zhì)與全氮、堿解氮、全磷、有效磷、全鉀、速效鉀、有效硼之間均達(dá)到極顯著正相關(guān)關(guān)系。全氮與堿解氮、全磷與有效磷達(dá)到極顯著正相關(guān)關(guān)系。堿解氮與全磷、有效磷、全鉀、速效鉀、有效硼均達(dá)到極顯著正相關(guān)關(guān)系。全磷與有效磷、全鉀、速效鉀、有效硼也達(dá)到極顯著正相關(guān)關(guān)系。有效磷與速效鉀、有效硼也達(dá)到極顯著正相關(guān)關(guān)系。各養(yǎng)分因子相關(guān)性系數(shù)平均值和權(quán)重系數(shù)見(jiàn)表5，從表5中可看出，有機(jī)質(zhì)權(quán)重系數(shù)最大，達(dá)到0.1656，全鉀權(quán)重系數(shù)最小，為0.0697。

2.2.3 土壤養(yǎng)分綜合指數(shù)計(jì)算 綜合肥力指標(biāo)值IFI（Integrated Fertility Index）是一個(gè)全面反映土壤養(yǎng)分狀況的指標(biāo)值，其值大小表示土壤養(yǎng)分綜合肥力等級(jí)，用加乘法則求得，根據(jù)上述求得的各養(yǎng)分評(píng)價(jià)指標(biāo)的隸屬度值和相應(yīng)的權(quán)重，進(jìn)一步求得各土壤樣品的IFI值。山南市土壤養(yǎng)分綜合指數(shù)變幅在0.146～0.993，平均值0.527，標(biāo)準(zhǔn)差0.162，變異系數(shù)32.82%。山南市土壤養(yǎng)分綜合指數(shù)的頻數(shù)分布情況見(jiàn)圖1。

2.2.4 土壤肥力水平等級(jí)的劃分 計(jì)算出的土壤肥力綜合指標(biāo)值，按一定的方法即可劃分土壤肥力等級(jí)。常用的劃分方法有等間距法、依據(jù)曲線拐點(diǎn)的非等間距法、系統(tǒng)聚類分析法等，本研究采用等間距法[5]。本研究采用等間距法將山南市土壤肥力劃分為5個(gè)等級(jí)：Ⅰ級(jí)（IFI值0.8～1），Ⅱ級(jí)（IFI值0.6～0.8），Ⅲ級(jí)（IFI值0.4～0.6），Ⅳ級(jí)（IFI值0.2～0.4），Ⅴ級(jí)（IFI值0～0.2）。其中Ⅴ級(jí)土壤肥力的樣品最少占該區(qū)總樣品的1.29%，次之是屬于Ⅰ級(jí)的樣品占總樣品的6.66%。屬于Ⅲ級(jí)的土壤樣品最多，占總樣品的47.08%（表6）。

3 結(jié)論與討論

3.1 結(jié)論 山南市農(nóng)耕地養(yǎng)分等級(jí)在中等水平，其中一等地57個(gè)樣本，所占比例為6.66%，二等地197個(gè)樣本，占23.01%；三等地403個(gè)樣本，其所占比例為47.08%；四等地188個(gè)樣本，占21.96%；五等地11個(gè)樣本，所占比例1.29%。由分析結(jié)果可見(jiàn)，山南市主要農(nóng)耕地土壤肥力等級(jí)以三等地為主，處于中等水平。土壤養(yǎng)分指標(biāo)中的有機(jī)質(zhì)、全磷、有效磷、全鉀、速效鉀含量水平比較低，76.6%的土樣有機(jī)質(zhì)含量小于30g/kg，處于缺乏狀態(tài)；69.04%土樣堿解氮含量小于120mg/kg。59%土樣有效磷含量小于10mg/kg；92.52%的土樣全磷含量小于0.8g/kg；79.32%的土樣速效鉀含量小于100mg/kg；68.12%的土樣全鉀含量小于15g/kg；76.75%的土樣有效硼含量水平豐富。土壤肥力綜合模型充分考慮了評(píng)價(jià)因素指標(biāo)值，評(píng)價(jià)因素權(quán)重和評(píng)價(jià)因素之間交互作用對(duì)土壤肥力質(zhì)量的共同作用，與加權(quán)指數(shù)模型相比能夠綜合反映土壤肥力質(zhì)量的狀況。本次研究區(qū)土壤肥力綜合模型評(píng)價(jià)結(jié)果表明，全鉀和有效磷是影響土壤肥力綜合指標(biāo)的限制因子。

3.2 討論 研究區(qū)處于青藏高原，地理、氣候條件特殊，若再考慮土壤的某些物理性狀指標(biāo)（如土壤質(zhì)地等）、灌溉能力、海拔、成土母質(zhì)等立地條件，結(jié)果會(huì)更加客觀、全面。應(yīng)當(dāng)注意的是，土壤肥力綜合評(píng)價(jià)的結(jié)果僅代表一種潛在的生產(chǎn)能力，在實(shí)際應(yīng)用中，還要結(jié)合作物品種特性、栽培制度、氣候特點(diǎn)、作物需肥規(guī)律以及土壤養(yǎng)分的有效性等統(tǒng)籌考慮。研究區(qū)的整體土壤肥力水平中等，氣候特殊，總體上應(yīng)注意增加有機(jī)肥料的施入量，前期增施基肥應(yīng)與后期追肥相結(jié)合；在配方施肥策略上，應(yīng)根據(jù)大田作物營(yíng)養(yǎng)診斷實(shí)施平衡施肥技術(shù)，注意“有機(jī)與無(wú)機(jī)”、“增氮與控氮”相結(jié)合、適當(dāng)增施磷肥和鉀肥。

參考文獻(xiàn)

[1]武偉，唐明華，劉洪斌.土壤養(yǎng)分的模糊綜合評(píng)價(jià)[J].西南農(nóng)業(yè)大學(xué)學(xué)報(bào)，2000，22（3）：270-272.

[2]吳玉紅，田霄鴻，同延安，等.基于主成分分析土壤綜合指標(biāo)評(píng)價(jià)[J].生態(tài)學(xué)雜志，2010，29（1）：173-180.

[3]徐建明，張甘霖，謝正苗，等.土壤質(zhì)量指標(biāo)與評(píng)價(jià)[M].北京：科學(xué)出版社，2010.

[4]王建國(guó)，楊林章，單艷紅.模糊數(shù)學(xué)在土壤質(zhì)量評(píng)價(jià)中的應(yīng)用研究[J].土壤學(xué)報(bào)，2001，38（2）：176-183.

[5]邵慧芳，鄭聰，許自成，等.三門峽煙區(qū)土壤養(yǎng)分狀況的綜合評(píng)價(jià)[J].西南農(nóng)業(yè)學(xué)報(bào)，2009，22（4）：1011-1015.

[6]楊劍虹，王成秋，代亨林，等.土壤農(nóng)化分析與環(huán)境監(jiān)測(cè)[M].北京：中國(guó)大地出版社，2008.

[7]呂蘇丹，汪光宇，鄔亞浪，等.東陽(yáng)萬(wàn)畝園區(qū)土壤養(yǎng)分綜合評(píng)價(jià)研究[J].浙江大學(xué)學(xué)報(bào)（農(nóng)業(yè)與生命科學(xué)出版），2002，28（3）：272-276.

[8]陸允甫，呂曉男.中國(guó)測(cè)土施肥工作的進(jìn)展和展望[J].土壤學(xué)報(bào)，1995，32（3）：241-252.

[9]駱東奇，白潔，謝德體.論土壤肥力評(píng)價(jià)指標(biāo)和方法[J].土壤與環(huán)境，2002，11（2）：202-205.

[10]李忠意，程永毅，宋偉，等.重慶紫色土壤養(yǎng)分的模糊綜合評(píng)判[J].中國(guó)農(nóng)學(xué)通報(bào)，2011，27（23）：238-247.

[11]李大皓，崔守富，唐海波.黑龍江省齊齊哈爾市區(qū)耕地地力評(píng)價(jià)[M].北京：中國(guó)農(nóng)業(yè)出版社，2016.

[12]郭亞軍，于兆吉.綜合評(píng)價(jià)的合理性問(wèn)題[J].東北大學(xué)學(xué)報(bào)（自然科學(xué)版），2002，23（9）：844-847.

[13]孫波，張?zhí)伊?，趙其國(guó).我國(guó)東南丘陵山區(qū)土壤肥力的綜合評(píng)價(jià)[J].土壤學(xué)報(bào)，1995，324，362-369.

[14]呂曉男，陸允甫，王人朝.土壤肥力綜合評(píng)價(jià)初步研究[J].浙江大學(xué)學(xué)報(bào)農(nóng)業(yè)與生命科學(xué)，1999，25（4）：378-382.

[15]黎孟波，張先婉.土壤肥沃度理化評(píng)價(jià)：土壤肥力研究之三[A]//土壤肥力研究進(jìn)展[C].北京：中國(guó)科學(xué)技術(shù)出版社，1991：221-234. （責(zé)編：張宏民）