崔增團(tuán)，郭世乾

(甘肅省農(nóng)業(yè)節(jié)水與土壤肥料管理總站，蘭州 730020)

耕地質(zhì)量是確保國(guó)家糧食安全和農(nóng)產(chǎn)品質(zhì)量安全的基礎(chǔ)，沒(méi)有可靠的耕地質(zhì)量，就不可能從根本上保障我國(guó)糧食安全和食品安全。耕地地力是耕地質(zhì)量的集中反映，包括耕地的立地條件、土壤條件和農(nóng)田基礎(chǔ)設(shè)施及培肥水平等［1-2］。近年來(lái)，隨著“3S”技術(shù)不斷發(fā)展，已在農(nóng)業(yè)領(lǐng)域得到廣泛應(yīng)用。如吳立忠以麥積區(qū)為例對(duì)黃土高原丘陵溝壑區(qū)耕地地力進(jìn)行了評(píng)價(jià)研究［3］，魯明星等以鄂州市為例選擇了13個(gè)評(píng)價(jià)因子開(kāi)展了基于GIS的華中丘陵區(qū)的耕地地力評(píng)價(jià)研究［4］。林碧珊等以土種為單元，采用限制因素法和綜合歸納法，對(duì)廣東省耕地地力因素進(jìn)行系統(tǒng)分析和評(píng)比，評(píng)價(jià)耕地 (土種)地力等級(jí)［5］。

河西走廊灌溉農(nóng)業(yè)區(qū)是我國(guó)著名的灌溉農(nóng)業(yè)區(qū)，也是我國(guó)主要的商品糧生產(chǎn)基地、玉米制種基地及經(jīng)濟(jì)作物集中生產(chǎn)基地。該研究結(jié)合農(nóng)業(yè)部耕地地力調(diào)查項(xiàng)目，借助于GIS技術(shù)，綜合運(yùn)用層次分析、模糊數(shù)學(xué)方法，以肅州區(qū)為例，探究其影響河西走廊灌溉農(nóng)業(yè)區(qū)耕地地力的主要因素，開(kāi)展耕地地力評(píng)價(jià)，摸清河西走廊灌溉農(nóng)業(yè)區(qū)耕地綜合生產(chǎn)能力。

1 研究區(qū)域情況與條件

1.1 研究區(qū)域基本情況

甘肅省肅州區(qū)位于河西走廊中部的西端，東接高臺(tái)縣，北連金塔縣，西鄰嘉峪關(guān)市，南與肅南縣接壤，地處北緯39°10＇～39°59＇，東經(jīng) 98°20＇～99°18＇范圍之內(nèi)，海拔 1 340 ～2 200m 之間。東西長(zhǎng)約104km，南北寬約84km，總面積3 449.5km2，耕地面積4.58萬(wàn)hm2。年平均氣溫7.3℃，年平均降水量85.3mm，年平均蒸發(fā)量為2 148.8mm，年平均日照時(shí)數(shù)為3 033.4小時(shí)。可利用的水資源總量為15.46億m3，其中可利用的地表水6.98億m3，可利用的地下水8.52億m3，是典型的內(nèi)陸干旱荒漠綠洲灌溉農(nóng)業(yè)區(qū)，沒(méi)有灌溉就沒(méi)有農(nóng)業(yè)。全區(qū)轄15個(gè)鄉(xiāng)鎮(zhèn)，10個(gè)國(guó)營(yíng)農(nóng)林場(chǎng)，124個(gè)村民委員會(huì)，790個(gè)村民小組。農(nóng)村總戶數(shù)5.87萬(wàn)戶，農(nóng)村總?cè)丝?1.78萬(wàn)人，農(nóng)村從業(yè)人員11.97萬(wàn)人。境內(nèi)地勢(shì)平坦，土壤肥沃，水源充足，是一個(gè)人口密集，人多地少，耕地資源十分緊張的區(qū) (縣)。耕種土壤主要為灌漠土和潮土。

1.2 研究條件

屬性數(shù)據(jù)庫(kù)的建立主要應(yīng)用ArcGIS 9.3及Microsoft Office Access 2007數(shù)據(jù)庫(kù)軟件，空間數(shù)據(jù)庫(kù)的建立主要應(yīng)用ArcGIS 9.3軟件，評(píng)價(jià)指標(biāo)的隸屬函數(shù)模型建立與層次分析主要應(yīng)用SPSS 13.0進(jìn)行分析。

2 研究方法

2.1 評(píng)價(jià)單元的確定

利用ArcGIS 9.3軟件對(duì)1∶5萬(wàn)土地利用現(xiàn)狀圖、土壤圖、行政區(qū)劃圖、水利分區(qū)圖、地貌類(lèi)型圖進(jìn)行疊加求交，共生成4 526個(gè)評(píng)價(jià)單元，使單元空間界線及行政隸屬關(guān)系明確，有準(zhǔn)確的面積，地貌類(lèi)型及土壤類(lèi)型一致，利用方式及耕作方法相同，便于農(nóng)業(yè)結(jié)構(gòu)調(diào)整規(guī)劃和指導(dǎo)實(shí)際的農(nóng)事操作為原則。

2.2 評(píng)價(jià)指標(biāo)的確定

2.2.1 確定原則

以選取的指標(biāo)對(duì)耕地生產(chǎn)力有較大的影響、在評(píng)價(jià)區(qū)域內(nèi)的變異較大、在時(shí)間序列上具有相對(duì)的穩(wěn)定性、與當(dāng)前生產(chǎn)密切相關(guān)等為原則。

2.2.2 確定方法

邀請(qǐng)18位有關(guān)專(zhuān)家在全國(guó)66個(gè)評(píng)價(jià)指標(biāo)中采用特爾斐法和模糊評(píng)價(jià)法進(jìn)行選取。

2.2.3 評(píng)價(jià)指標(biāo)

確定土壤侵蝕程度、海拔高度、有效土層厚度、質(zhì)地構(gòu)型、剖面構(gòu)型、有機(jī)質(zhì)、有效磷、速效鉀、障礙層厚度、1m土層含鹽量、障礙層出現(xiàn)位置、灌溉模數(shù)、灌溉保證率13項(xiàng)指標(biāo)耕地地力評(píng)價(jià)指標(biāo)。

2.3 數(shù)據(jù)庫(kù)的建立

海拔高度指標(biāo)的獲取，利用ArcGIS9.3軟件對(duì)肅州區(qū)地形圖進(jìn)行DEM分析以后，采用區(qū)域統(tǒng)計(jì)加權(quán)平均的方法將值賦給相應(yīng)的評(píng)價(jià)單元。

有機(jī)質(zhì)、有效磷、速效鉀、1m土層含鹽量等指標(biāo)的獲取，是以肅州區(qū)耕地地力調(diào)查項(xiàng)目中用GPS準(zhǔn)確定位采集的4 600個(gè)耕層調(diào)查樣，按照按常規(guī)方法進(jìn)行化驗(yàn)分析［6］，經(jīng)地統(tǒng)計(jì)分析符合空間插值條件，利用ArcGIS中普通Kriging插值方法進(jìn)行空間插值。

有效土層厚度、質(zhì)地、剖面構(gòu)型、障礙層厚度、障礙層出現(xiàn)位置、土壤侵蝕程度是通過(guò)土壤圖關(guān)聯(lián)屬性數(shù)據(jù)以后，生成評(píng)價(jià)單元過(guò)程中直接賦值。

灌溉模數(shù)、灌溉保證率是通過(guò)1∶5萬(wàn)灌溉分區(qū)圖與管理單元圖疊加以后，進(jìn)行屬性提取。

2.4 單因素評(píng)價(jià)

2.4.1 概念性指標(biāo)評(píng)價(jià)

概念型指標(biāo)評(píng)價(jià)采用專(zhuān)家經(jīng)驗(yàn)法，直接對(duì)評(píng)價(jià)單元進(jìn)行模糊評(píng)價(jià)確定隸屬度 (表1、表2、表3、表4，表5)。

表1 土壤侵蝕程度隸屬度

表2 質(zhì)地構(gòu)型隸屬度

表3 灌溉保證率隸屬度

表4 灌溉模數(shù)隸屬度

表5 剖面構(gòu)型隸屬度

2.4.2 數(shù)值性指標(biāo)評(píng)價(jià)

數(shù)值型指標(biāo)采用特爾斐法和模糊評(píng)價(jià)法，用SPSS 13.0軟件建立評(píng)價(jià)指標(biāo)實(shí)際值和代表生產(chǎn)力水平的隸屬度 (0～1值)間隸屬函數(shù)模型 (表6)。

表6 數(shù)值型評(píng)價(jià)指標(biāo)的隸屬函數(shù)類(lèi)型

2.5 計(jì)算評(píng)價(jià)指標(biāo)權(quán)重

表7 各個(gè)因素的組合權(quán)重計(jì)算結(jié)果

采用層次分析法確定每一個(gè)評(píng)價(jià)因素對(duì)耕地地力貢獻(xiàn)大小。依據(jù)14項(xiàng)指標(biāo)要素間的關(guān)系，構(gòu)造層次結(jié)構(gòu)圖，比較同一層次各因素對(duì)上層次的相對(duì)重要性，給出數(shù)字化的評(píng)估，形成判斷矩陣，計(jì)算層次分析結(jié)果 (表7)。

2.6 計(jì)算耕地地力綜合指數(shù) (IFI)

利用前面2.2節(jié)至2.5節(jié)中選擇的評(píng)價(jià)指標(biāo)及各指標(biāo)隸屬度、隸屬函數(shù)模型和組合權(quán)重，采用加法模型 (下式)計(jì)算耕地地力綜合指數(shù)IFI:

式中:IFI(Integrated Fertility Index)代表耕地地力指數(shù);

F i=第i個(gè)因素隸屬度;CiBi=第i個(gè)因素的組合權(quán)重。

2.7 確定耕地地力綜合指數(shù)分級(jí)方案

利用SPSS 13.0分析軟件評(píng)價(jià)耕地地力，采取累積曲線分級(jí)法劃分耕地地力等級(jí)，將肅州區(qū)耕地地力劃分為4級(jí) (表8)。

2.8 確定國(guó)家耕地地力等級(jí)

在每一個(gè)地力等級(jí)上選取10%的管理單元，獲得了453個(gè)管理單元上的IFI值，并獲取了所選取的管理單元上近3年實(shí)際的年平均產(chǎn)量，然后對(duì)所選管理單元上的IFI值和實(shí)際產(chǎn)量進(jìn)行回歸模擬，r=0.89＊＊，計(jì)算各相應(yīng)等級(jí)上的概念性產(chǎn)量，按照農(nóng)業(yè)部《全國(guó)耕地類(lèi)型區(qū)耕地地力等級(jí)劃分》標(biāo)準(zhǔn)歸入以概念性產(chǎn)量表示的地力等級(jí)體系 (表8)。

表8 肅州區(qū)耕地地力等級(jí)與國(guó)家地力等級(jí)劃分

3 結(jié)果與分析

3.1 各級(jí)等級(jí)基本情況

由表8可見(jiàn)，肅州區(qū)總耕地面積為4.58萬(wàn)hm2，耕地地力等級(jí)可分為4個(gè)等級(jí)，一、二、三等地占的比例較大，分別為36.1%、37.2%和21.6%，占總耕地的94.9%，四等地占的比例較小，為5.1%。根據(jù)《全國(guó)耕地類(lèi)型區(qū)耕地地力等級(jí)劃分》農(nóng)業(yè)行業(yè)標(biāo)準(zhǔn)，肅州區(qū)一至四等地，分別對(duì)應(yīng)國(guó)家三至六等地。

3.2 耕地地力等級(jí)空間分布特點(diǎn)及土壤屬性分析

由圖1、表9可知，肅州區(qū)一等地主要分布在肅州區(qū)北部地區(qū)，土種以灌淤土和潮土為主，地貌類(lèi)型以?xún)A斜沖洪積平原為主，耕層質(zhì)地構(gòu)型以輕壤質(zhì)和中壤質(zhì)為主，灌溉保證率均在75%以上，土層厚，有效土層厚度大于1m，地勢(shì)平坦，降水量少，無(wú)明顯障礙層，耕層土壤鹽分平均含量為0.26%，土壤理化性狀良好，可耕性強(qiáng)。耕層土壤pH值為8.56，有機(jī)質(zhì)14.86g/kg，有效鋅 1.19mg/kg，有效磷 18.39mg/kg，速效鉀 160.69mg/kg，堿解氮 72.45mg/kg。

圖1 肅州區(qū)耕地地力等級(jí)

二等地主要分布在肅州區(qū)中北部的各鄉(xiāng)鎮(zhèn)以及西南部的部分鄉(xiāng)鎮(zhèn)，土種以潮土和灰棕漠土為主，地貌類(lèi)型以淺山區(qū)和低洼地為主，耕層質(zhì)地構(gòu)型以均質(zhì)輕壤、均質(zhì)中壤和夾黏中壤為主，灌溉保證率在65%以上，有效土層厚度在29～100cm，降雨量少，有明顯障礙層，耕層土壤鹽分平均含量為0.22，有輕度鹽堿危害，可耕性較強(qiáng)。耕層土壤pH值為8.54，有機(jī)質(zhì)14.27g/kg，有效鋅1.22mg/kg，有效磷18.10mg/kg，速效鉀149.53mg/kg，堿解氮68.68mg/kg。主要影響因素是土壤障礙層位置和厚度。

表9 不同等級(jí)土壤屬性分析

三等地主要分布在肅州區(qū)中部西洞鎮(zhèn)、清水鎮(zhèn)、黃泥堡鄉(xiāng)、豐樂(lè)鄉(xiāng)、鏵尖鄉(xiāng)，土種以灰棕漠土為主，耕層質(zhì)地構(gòu)型以均質(zhì)輕壤、夾黏中壤和夾砂輕壤為主，灌溉保證率均在65%以下，有效土層厚度在29～100cm，降雨量少，海拔高，有明顯障礙層，1m土層鹽分含量平均值為0.26%，土壤理化性狀良好，可耕性較強(qiáng)。耕層土壤pH值為8.54，有機(jī)質(zhì)14.68g/kg，有效鋅1.20mg/kg，有效磷18.14mg/kg，速效鉀157.07mg/kg，堿解氮70.25mg/kg。主要影響因素是鹽堿危害，土層中存在障礙層。

四等地主要分布在肅州區(qū)祁連山一帶。土種以山地灰棕漠土為主，地貌類(lèi)型為緩坡地，耕層質(zhì)地構(gòu)型以均質(zhì)中壤和均質(zhì)輕壤為主，灌溉保證率均在65%以上，有效土層厚度22～48cm，土層薄，有明顯障礙層，耕層土壤鹽分含量平均值為0.25%;受鹽堿和風(fēng)蝕危害，可耕性較好。耕層土壤pH值為8.54，有機(jī)質(zhì) 13.97g/kg，有效鋅 1.19mg/kg，有效磷 17.78mg/kg，速效鉀 149.91mg/kg，堿解氮 69.53mg/kg。主要影響因素為灌溉保證率差。

4 結(jié)論

(1)以甘肅省肅州區(qū)為例，應(yīng)用GIS技術(shù)，結(jié)合層次分析法和模糊數(shù)學(xué)綜合評(píng)價(jià)方法等現(xiàn)代統(tǒng)計(jì)分析技術(shù)對(duì)河西走廊灌溉農(nóng)業(yè)區(qū)耕地地力評(píng)價(jià)指標(biāo)權(quán)重分析結(jié)果表明，剖面構(gòu)型、障礙因素、土壤管理是影響該區(qū)域耕地生產(chǎn)潛力最重要的因素，其重要程度依次為:剖面構(gòu)型＞障礙層位置＞質(zhì)地構(gòu)型＞灌溉保證率＞土體含鹽量＞有機(jī)質(zhì)＞有效土層厚＞海拔高度＞障礙層厚度＞有效磷＞灌溉模數(shù)＞土壤侵蝕程度＞速效鉀。

(2)肅州區(qū)耕地地力等級(jí)可分為4等，一、二、三等地占的比例較大，分別為36.1%、37.2%和21.6%，占總耕地的94.9%，四等地占的比例較小，為5.1%。與國(guó)家耕地地力等級(jí)對(duì)比，肅州區(qū)一至四等地，分別對(duì)應(yīng)國(guó)家三至六等地，該區(qū)域耕地綜合生產(chǎn)能力普遍較高。

(3)從肅州區(qū)耕地地力評(píng)價(jià)結(jié)果及屬性分析，一等地應(yīng)主要測(cè)土推廣配方施肥技術(shù)，加強(qiáng)耕地保護(hù)和建設(shè);二等地除推廣測(cè)土配方施肥外，還要注意深耕，加強(qiáng)水肥管理;三等地應(yīng)大力推廣防鹽堿技術(shù)，增施有機(jī)肥，推廣綠肥種植，還應(yīng)加強(qiáng)水利設(shè)施建設(shè)，提高灌溉保證水平;四等地應(yīng)大力發(fā)展農(nóng)業(yè)節(jié)水技術(shù)，如應(yīng)推廣地膜覆蓋防滲技術(shù)，渠道防滲技術(shù)，膜下滴灌技術(shù)，提高水分利用效率。

［1］王蓉芳，曹富有，彭世琪，等.中國(guó)耕地基礎(chǔ)地力與土壤改良.北京:中國(guó)農(nóng)業(yè)出版社，1996

［2］趙其國(guó)，孫波，張?zhí)伊?土壤質(zhì)量與持續(xù)環(huán)境:l.土壤質(zhì)量的定義及評(píng)價(jià)方法.土壤，1997，29(3):113～120

［3］吳立忠，郭世乾，張保田.基于GIS的黃土高原丘陵溝壑區(qū)的耕地地力評(píng)價(jià)研究——以天水市麥積區(qū)耕地地力評(píng)價(jià)為例.中國(guó)農(nóng)業(yè)資源與區(qū)劃，2009，30(6):35～40

［4］魯明星，賀立源，吳禮樹(shù)，等.基于GIS的華中丘陵溝壑區(qū)的耕地地力評(píng)價(jià)研究.農(nóng)業(yè)工程學(xué)報(bào)，2006，22(8):96～101

［5］林碧珊，湯建東，張滿紅.廣東省耕地地力等級(jí)研究與評(píng)價(jià).生態(tài)環(huán)境，2005，14(1):145～149

［6］鮑士旦.土壤農(nóng)化分析.北京:中國(guó)農(nóng)業(yè)出版社.2000