王夢(mèng)銳，畢如田

（山西農(nóng)業(yè)大學(xué)資源環(huán)境學(xué)院，山西太谷030801）

0 引言

耕地細(xì)碎化不僅包括耕地自然斑塊的細(xì)碎化，還包括耕地經(jīng)營(yíng)權(quán)屬的細(xì)碎化[1-2]。耕地自然斑塊的細(xì)碎化又叫耕地自然細(xì)碎化，它是由于耕地受到農(nóng)業(yè)資源、氣候?yàn)?zāi)害、地形條件等因素的影響，使得耕地出現(xiàn)地塊大小不一，質(zhì)量不等，數(shù)目較多，地塊分散的問(wèn)題[3]。耕地經(jīng)營(yíng)權(quán)屬的細(xì)碎化又叫耕地經(jīng)營(yíng)細(xì)碎化，是社會(huì)經(jīng)濟(jì)、人文地理、科學(xué)技術(shù)等各種因素綜合影響的結(jié)果，耕地受農(nóng)戶(hù)耕作權(quán)屬影響，將一塊自然斑塊耕地人為分割成了面積不同的幾個(gè)田塊，一般情況下，每個(gè)農(nóng)戶(hù)的耕地田塊分散，但分布在一定空間范圍內(nèi)[4-6]。尤其中國(guó)實(shí)行家庭聯(lián)產(chǎn)承包責(zé)任制以來(lái)，極大地促進(jìn)了中國(guó)農(nóng)業(yè)的發(fā)展，激發(fā)了農(nóng)民生產(chǎn)熱情，然而該制度下耕地平均分配原則使耕地細(xì)碎化進(jìn)一步加深[7]，造成耕地浪費(fèi)、規(guī)?；a(chǎn)受阻[10-11]和生產(chǎn)效率低下[12-14]等問(wèn)題。19世紀(jì)人們才開(kāi)始研究耕地細(xì)碎化對(duì)農(nóng)業(yè)的影響，直到20世紀(jì)初，解決耕地細(xì)碎化問(wèn)題的研究正式起步[13]。國(guó)外對(duì)耕地細(xì)碎化的定義雖然沒(méi)有中國(guó)的權(quán)威，但國(guó)外在這方面的研究早于中國(guó)很多年，其中在耕地細(xì)碎化定義方面，Binns[14]認(rèn)為耕地細(xì)碎化是在農(nóng)場(chǎng)等較大范圍內(nèi)，很多單個(gè)面積較小且不相連的地塊空間分布呈現(xiàn)出來(lái)的狀態(tài)；Papageorgiou和Tin Nguyen等[15-16]則強(qiáng)調(diào)地塊間的距離和多于一塊離散地塊的經(jīng)營(yíng)產(chǎn)權(quán)；Dovring[17]把同一塊自然斑塊耕地上有不同的經(jīng)營(yíng)產(chǎn)權(quán)和相同產(chǎn)權(quán)地塊不在同一個(gè)自然斑塊的地塊定義為土地細(xì)碎化；Bellon[18]強(qiáng)調(diào)耕地細(xì)碎化是耕地經(jīng)營(yíng)產(chǎn)權(quán)的分散。在農(nóng)業(yè)生產(chǎn)方面，Schultz[19]認(rèn)為土地的細(xì)碎化是農(nóng)場(chǎng)中土地低效率利用、管理不合理的體現(xiàn)；Jabarin等和Hazarika等[20-21]的研究表明地塊破碎度越高生產(chǎn)成本越高，地塊越大對(duì)農(nóng)業(yè)生產(chǎn)效率的提高越明顯。在中國(guó)，譚淑豪[22]教授分別從耕作地塊大小、經(jīng)營(yíng)產(chǎn)權(quán)分散、地塊間距離定義耕地細(xì)碎化；李功奎[2]認(rèn)為耕地細(xì)碎化能夠提供多元化種植，有效發(fā)揮剩余勞動(dòng)力；萬(wàn)廣華、蘇旭霞等[8-9]調(diào)查指出耕地細(xì)碎化限制了農(nóng)業(yè)的規(guī)模化生產(chǎn)；呂挺、紀(jì)月清等[23]的研究顯示地塊面積對(duì)水稻的單產(chǎn)無(wú)差異，但面積較大的地塊生產(chǎn)過(guò)程中生產(chǎn)成本尤其是勞動(dòng)成本較低。

在中國(guó)偏遠(yuǎn)山區(qū)，受自然條件、社會(huì)基礎(chǔ)設(shè)施等影響較大，加上家庭聯(lián)產(chǎn)承包責(zé)任制的影響，大多只能實(shí)行傳統(tǒng)農(nóng)業(yè)耕作方式，生產(chǎn)水平較低。因此明晰耕地空間分布規(guī)律，對(duì)耕地質(zhì)量重新分類(lèi)，調(diào)整村落內(nèi)部耕地承包經(jīng)營(yíng)權(quán)，改善小面積、低效率的生產(chǎn)方式，實(shí)現(xiàn)農(nóng)業(yè)規(guī)?；鹘?jīng)營(yíng)，改善農(nóng)業(yè)生產(chǎn)效率等成為丘陵地區(qū)耕地流轉(zhuǎn)問(wèn)題的關(guān)鍵。目前耕地質(zhì)量評(píng)價(jià)分級(jí)方面的研究有很多，但基于農(nóng)村田塊特征的耕地分類(lèi)及實(shí)現(xiàn)同類(lèi)地塊間承包經(jīng)營(yíng)權(quán)調(diào)整的研究較少。筆者以山西省太谷縣溫莊村為例，基于研究區(qū)田塊特征，對(duì)耕地進(jìn)行重新分類(lèi)，通過(guò)調(diào)整村落內(nèi)部耕地承包經(jīng)營(yíng)權(quán)，為解決耕地細(xì)碎化問(wèn)題，實(shí)現(xiàn)農(nóng)業(yè)規(guī)?；a(chǎn)，提高農(nóng)業(yè)生產(chǎn)率提供參考依據(jù)。

1 研究區(qū)概況與數(shù)據(jù)來(lái)源

1.1 研究區(qū)概況

溫莊村地處山西省太谷縣范村鎮(zhèn)東北部，距離太谷縣縣城32 km，北與坪上村相接，西鄰上安村，東接下寨村。研究區(qū)位于丘陵地帶，屬暖溫帶大陸性氣候，全村約有耕地32 hm2，戶(hù)數(shù)31戶(hù)，人均年收入2000元。全年降水季節(jié)分布不均勻，作物生長(zhǎng)靠人工灌溉來(lái)促進(jìn)；土壤呈微堿性，并且侵蝕較為嚴(yán)重；植被種類(lèi)豐富，玉米和棗樹(shù)為主要的農(nóng)作物，經(jīng)濟(jì)來(lái)源以農(nóng)業(yè)為主，占農(nóng)民收入的70%。研究區(qū)為典型的丘陵地貌，常住人口稀少，經(jīng)濟(jì)發(fā)展速度緩慢，全村耕地分布在山區(qū)和平地的地塊數(shù)目均衡，且地塊形狀各異，面積大小和便捷程度不同，相同承包經(jīng)營(yíng)權(quán)的耕地空間分布散亂、耕地細(xì)碎化現(xiàn)象嚴(yán)重，該村地塊特征對(duì)研究欠發(fā)達(dá)地區(qū)村內(nèi)耕地經(jīng)營(yíng)權(quán)調(diào)整、解決耕地細(xì)碎化問(wèn)題具有一定的代表性。

1.2 數(shù)據(jù)來(lái)源

矢量數(shù)據(jù)來(lái)源于2017年太谷縣農(nóng)村土地承包經(jīng)營(yíng)權(quán)數(shù)據(jù)庫(kù)，經(jīng)ArcGIS提取、投影轉(zhuǎn)換獲取耕地?cái)?shù)據(jù)；收集DEM數(shù)據(jù)，經(jīng)裁剪、投影轉(zhuǎn)換為與耕地?cái)?shù)據(jù)一致坐標(biāo)系的柵格數(shù)據(jù)；道路矢量數(shù)據(jù)根據(jù)太谷縣1:2000正射影像，并結(jié)合田間實(shí)地調(diào)查繪制獲得；田塊屬性數(shù)據(jù)及產(chǎn)量數(shù)據(jù)通過(guò)實(shí)地調(diào)查訪問(wèn)獲得。

2 研究方法

2.1 村落尺度耕地質(zhì)量影響因子的選取與計(jì)算

本研究選取了地塊面積、地塊坡度、地塊形狀指數(shù)、地塊連片性、地塊規(guī)整度、耕作通達(dá)度、耕作遠(yuǎn)距離及地塊產(chǎn)量等8個(gè)丘陵區(qū)耕地細(xì)碎化的影響因子，并在在各因素的影響下分別對(duì)研究區(qū)地塊進(jìn)行分類(lèi)，分析地塊的空間分布特征。

（1）地塊面積。地塊面積是衡量耕地細(xì)碎化的一項(xiàng)重要指標(biāo)，地塊面積大小與農(nóng)戶(hù)對(duì)地塊的成本投入、經(jīng)營(yíng)方式及產(chǎn)出利益有一定的關(guān)系[24]。研究區(qū)地塊屬性數(shù)據(jù)提取自2017年太谷縣農(nóng)村土地承包經(jīng)營(yíng)權(quán)確權(quán)數(shù)據(jù)庫(kù)，其中包括地塊面積，地塊周長(zhǎng)，地塊經(jīng)營(yíng)權(quán)屬以及農(nóng)戶(hù)的家庭情況等，通過(guò)地塊面積的提取發(fā)現(xiàn)該研究區(qū)地塊面積最大為0.56 hm2，最小為0.01 hm2，反映出耕地的細(xì)碎化較為嚴(yán)重，可見(jiàn)地塊面積可作為影響耕地細(xì)碎化的一個(gè)重要影響因素。

（2）地塊坡度。地形地貌是影響耕地細(xì)碎化的重要因子，反映在地塊屬性上即為地塊的坡度，地塊坡度影響著耕地的綜合價(jià)值，坡度越平緩耕作越便利，坡度較大的耕地養(yǎng)分流失會(huì)比較嚴(yán)重故不適合耕種。根據(jù)研究區(qū)的DEM數(shù)據(jù)提取坡度值，并通過(guò)計(jì)算地塊的平均坡度值確定地塊坡度大小。

（3）地塊形狀指數(shù)。地塊形狀指數(shù)用來(lái)衡量地塊的復(fù)雜程度，是反映地塊耕作形狀的重要指標(biāo)，往往地塊形狀越接近規(guī)則的矩形，農(nóng)戶(hù)耕作起來(lái)會(huì)更加方便，農(nóng)戶(hù)也更加傾向于規(guī)則的地塊。計(jì)算方法見(jiàn)公式(1)。

式中S為地塊的形狀指數(shù)，Li為第i個(gè)地塊的周長(zhǎng)，ai為第i個(gè)地塊的面積，0.25為正方形矯正常數(shù)，S值為1表示地塊的形狀為正方形，S值越大，地塊形狀越不規(guī)則。

（4）地塊連片性。地塊的連片性反應(yīng)的是耕地在空間分布上的細(xì)碎化，指同一類(lèi)型或同一耕作等級(jí)的地塊在空間分布上的連接程度[25]，地塊連片性由公式(2)計(jì)算可得。

式中P為地塊連片性，Si為連片面積現(xiàn)狀值，Smin為連片面積中地塊的最小值，Smax為連片面積中地塊的最大值，其中地塊連片度P的取值范圍為[0，1]。

（5）地塊規(guī)整度。地塊規(guī)整度是評(píng)價(jià)地塊質(zhì)量的一個(gè)重要因素,更是農(nóng)戶(hù)對(duì)地塊價(jià)值進(jìn)行評(píng)價(jià)時(shí)參考的一個(gè)重要的影響因素。地塊規(guī)整度計(jì)算見(jiàn)公式(3)。

式中K為地塊規(guī)整度，Si為第i個(gè)地塊的面積，Li為第i個(gè)地塊的周長(zhǎng)。K值越小表示地塊的形狀越規(guī)則，K值越大則地塊的形狀越復(fù)雜。

（6）耕作通達(dá)度。耕作通達(dá)度是在基礎(chǔ)設(shè)施方面影響耕地細(xì)碎化的一個(gè)重要因子，是指各地塊到道路的最短距離，耕作通達(dá)度在很大程度上決定了耕地的綜合價(jià)值，在實(shí)際耕作中農(nóng)戶(hù)更加傾向于通達(dá)度較好的耕地，所以較好通達(dá)度的自然地塊被經(jīng)營(yíng)權(quán)屬分割的更加細(xì)碎化。耕地通達(dá)度通過(guò)ArcGIS軟件工具箱里的near工具計(jì)算各地塊的幾何中心到可通農(nóng)用車(chē)輛道路的最近距離。

（7）耕作距離。耕作距離是指耕地到居民區(qū)的距離，距離越近，農(nóng)戶(hù)耕作越方便，耕作距離也是衡量耕地綜合質(zhì)量的重要因素。耕作距離根據(jù)ArcGIS軟件的鄰近分析得出，將研究區(qū)地塊和居民區(qū)的幾何中心進(jìn)行鄰近分析得出每塊耕地到居民區(qū)的距離。

（8）地塊產(chǎn)量。家庭聯(lián)產(chǎn)承包責(zé)任制下的耕地平均分配不僅包括田塊數(shù)量上的平均分配，還包括地塊肥田瘦地的平均分配，即在農(nóng)村土地承包責(zé)任制度下，村集體結(jié)合耕地肥瘦按人口或勞動(dòng)力平均分配給每個(gè)農(nóng)戶(hù)。而地塊產(chǎn)量是對(duì)耕地肥瘦情況最直觀的反映，因此，地塊產(chǎn)量是研究影響耕地細(xì)碎化的一項(xiàng)重要指標(biāo)，在耕地承包經(jīng)營(yíng)權(quán)調(diào)整的過(guò)程中農(nóng)戶(hù)將傾向于地塊產(chǎn)量較高的田塊。

2.2 聚類(lèi)分析

利用統(tǒng)計(jì)軟件SPSS中的系統(tǒng)聚類(lèi)分析法，采用歐氏距離和組平均法進(jìn)行耕地分類(lèi)。根據(jù)分類(lèi)結(jié)果，將地塊的空間分布情況在ArcGIS軟件中顯示。

3 結(jié)果與分析

3.1 單因子影響的地塊分布

本研究選取太谷縣溫莊村，該村經(jīng)濟(jì)來(lái)源以農(nóng)業(yè)為主，經(jīng)濟(jì)條件相對(duì)較差。區(qū)域內(nèi)的耕地地塊較均勻的分布在平原和山區(qū)，且細(xì)碎化問(wèn)題較為嚴(yán)重，對(duì)研究欠發(fā)達(dá)地區(qū)耕地細(xì)碎化問(wèn)題以及通過(guò)調(diào)整地塊承包經(jīng)營(yíng)權(quán)實(shí)現(xiàn)村集體內(nèi)部耕地規(guī)?；a(chǎn)具有特殊意義。

分析8項(xiàng)田塊空間數(shù)據(jù)對(duì)耕地分類(lèi)影響，其中地塊面積分類(lèi)根據(jù)地塊面積由小到大分成了5個(gè)區(qū)間；地塊坡度根據(jù)國(guó)土第二次調(diào)查耕地坡度分級(jí)標(biāo)準(zhǔn)將地塊坡度分為四類(lèi)：≤2°、2°～6°、6°～15°、15°～25°和＞25°；地塊形狀指數(shù)、地塊連片性、地塊規(guī)整度、耕作通達(dá)度及耕作距離等數(shù)據(jù)通過(guò)地塊數(shù)據(jù)矢量化，利用ArcGIS軟件分析計(jì)算得到，將計(jì)算結(jié)果值由小到大分為5個(gè)區(qū)間；將實(shí)際訪問(wèn)調(diào)查獲得的產(chǎn)量數(shù)據(jù)由高產(chǎn)到低產(chǎn)分為四類(lèi)。結(jié)果見(jiàn)圖1。

圖1 單因子影響下地塊分類(lèi)示意圖

由圖1和表1可知：（1）研究區(qū)地塊面積在第Ⅰ類(lèi)的地塊最多為124塊，面積區(qū)間為0.01～0.12 hm2，該區(qū)間田塊面積為7.10 hm2；田塊面積在第Ⅱ類(lèi)的面積區(qū)間為 0.12～0.23 hm2，地塊數(shù)目為 61 塊，總面積為9.91 hm2；隨著地塊面積的增大，地塊的數(shù)目在逐漸減小，可以看出研究區(qū)的地塊面積普遍較小，耕地經(jīng)營(yíng)細(xì)碎化比較嚴(yán)重。（2）研究區(qū)地塊坡度在6°～15°之間的地塊最多為111塊，面積也最大為14.75 hm2，占到該研究區(qū)總面積的60.99%，其次地塊坡度在2°～6°數(shù)目較多有83塊，面積達(dá)到8.35 hm2，占該研究區(qū)總面積的34.51%?？梢?jiàn)研究區(qū)耕地坡度較緩集中在2°～15°之間。（3）研究區(qū)地塊形狀指數(shù)在第Ⅰ類(lèi)的地塊較多有121塊，其面積占到研究區(qū)面積的69.35%，其次為第Ⅱ類(lèi)，地塊面積占到總面積26.44%，研究區(qū)地塊形狀較為規(guī)則。（4）研究區(qū)耕地地塊的連片程度較低，連片性在第Ⅱ類(lèi)的地塊數(shù)目最多，統(tǒng)計(jì)后為125塊，面積占到了研究區(qū)耕地面積的60.08%，研究區(qū)耕地連片性較差。（5）地塊的規(guī)整度在第Ⅳ類(lèi)的地塊較多為64塊，面積占到研究區(qū)耕地面積的36.24%，第Ⅲ類(lèi)中地塊僅次于第Ⅳ類(lèi)有58塊，面積達(dá)到總面積的30.98%，研究區(qū)耕作情況較為復(fù)雜。（6）研究區(qū)耕作通達(dá)度在第Ⅰ類(lèi)的耕地?cái)?shù)量最多為103塊，面積為10.70 hm2，占總面積的44.25%，第Ⅱ類(lèi)的地塊數(shù)相對(duì)較多為73塊，面積為9.38 hm2，占全部耕地面積的38.79%，由此可見(jiàn)該地區(qū)的耕地通達(dá)度較好。其中耕地到道路的最小距離為2.98 m，最大距離為214.26 m。（7）耕作距離在第Ⅰ類(lèi)的地塊數(shù)目較多為92塊，面積為7.75 hm2，占到研究區(qū)耕地總面積的32.02%，耕作距離在第Ⅱ類(lèi)的耕地?cái)?shù)目雖然不多，但地塊面積最大，達(dá)到9.50 hm2，占到研究區(qū)耕地總面積的39.26%，隨著耕作距離的漸遠(yuǎn)，耕作便利度下降，地塊的數(shù)目也隨著減少。（8）地塊產(chǎn)量在第Ⅰ類(lèi)的地塊數(shù)目較多，通過(guò)計(jì)算有97塊，面積達(dá)到9.49 hm2，占總面積的39.22%，產(chǎn)量較高的地塊分布在居民區(qū)周?chē)钠降?，產(chǎn)量較低的耕地分布距離居民區(qū)較遠(yuǎn)。

3.2 系統(tǒng)聚類(lèi)分析結(jié)果

由于單因子影響村落耕地經(jīng)營(yíng)細(xì)碎化的程度存在差異，使耕地重新分類(lèi)結(jié)果各不相同，為了減少耕地分類(lèi)結(jié)果的差異，實(shí)現(xiàn)村落內(nèi)部耕地經(jīng)營(yíng)權(quán)調(diào)整，研究8項(xiàng)單因子對(duì)耕地經(jīng)營(yíng)細(xì)碎化的綜合影響，采用系統(tǒng)聚類(lèi)法對(duì)所選因子數(shù)據(jù)在統(tǒng)計(jì)軟件中進(jìn)行分析，并對(duì)研究區(qū)地塊進(jìn)行分類(lèi)，聚類(lèi)結(jié)果如圖2所示。

表1 單因子影響下地塊分類(lèi)結(jié)果

圖2 聚類(lèi)分析結(jié)果示意圖

從圖2可以看出，聚類(lèi)分析下地塊分類(lèi)結(jié)果以居民區(qū)為中心成環(huán)形分成了5類(lèi)，其中第Ⅰ類(lèi)和第Ⅱ類(lèi)地塊面積相對(duì)較大，分別為9.21 hm2和9.81 hm2，兩類(lèi)地塊總面積占到全村耕地面積的78.61%，且從圖中可以看出第Ⅰ類(lèi)和第Ⅱ類(lèi)地塊距離居民區(qū)相對(duì)較近，距離居民區(qū)較遠(yuǎn)的第Ⅲ、Ⅳ、Ⅴ類(lèi)地塊分散且數(shù)目較少，三類(lèi)地塊面積為5.17 hm2，占全村耕地面積21.40%。通過(guò)系統(tǒng)聚類(lèi)法后耕地分類(lèi)結(jié)果，同類(lèi)地塊可以根據(jù)地塊產(chǎn)量系數(shù)進(jìn)行田塊承包經(jīng)營(yíng)權(quán)的調(diào)整，不同類(lèi)地塊依據(jù)農(nóng)戶(hù)意愿參考田塊產(chǎn)量系數(shù)進(jìn)行互換，將擁有相同承包經(jīng)營(yíng)權(quán)的耕地在空間分布上盡量集中在一起，使農(nóng)戶(hù)田塊數(shù)量由多變少且保持原有人均耕地質(zhì)量均衡或不變，降低耕地經(jīng)營(yíng)細(xì)碎化，實(shí)現(xiàn)耕地集中連片經(jīng)營(yíng)，便于機(jī)械化耕作，同時(shí)降低運(yùn)輸成本和減少勞作時(shí)間，提高生產(chǎn)效率。

表2 聚類(lèi)分析后地塊數(shù)目與面積

4 結(jié)論與討論

本文基于田塊特征選取了地塊面積、地塊坡度、地塊形狀指數(shù)、地塊連片性、地塊規(guī)整度、耕作通達(dá)度、耕作距離和地塊產(chǎn)量等8個(gè)因素來(lái)對(duì)研究區(qū)耕地的經(jīng)營(yíng)細(xì)碎化進(jìn)行研究，并將8項(xiàng)影響因子數(shù)據(jù)進(jìn)行系統(tǒng)聚類(lèi)分析。對(duì)系統(tǒng)聚類(lèi)后的分類(lèi)結(jié)果進(jìn)行分析得出同類(lèi)地塊以居民區(qū)為中心呈環(huán)形的方式分布在居民區(qū)周?chē)?，相近區(qū)間的耕作距離和地塊坡度的耕地地塊屬于同一個(gè)類(lèi)中，且同一個(gè)類(lèi)中的地塊具有相似的連片性和地塊產(chǎn)量，而同一類(lèi)中地塊具有不同的地塊面積、耕作通達(dá)度、地塊形狀指數(shù)和地塊規(guī)整度。基于多因素系統(tǒng)聚類(lèi)分析法對(duì)耕地進(jìn)行重新分類(lèi)的方法適用于農(nóng)村內(nèi)部耕地承包經(jīng)營(yíng)權(quán)的調(diào)整，為農(nóng)村內(nèi)部耕地互換，解決耕地細(xì)碎化，提高農(nóng)業(yè)生產(chǎn)效率提供了一種新思路。

本研究選取的研究區(qū)域較小，地貌類(lèi)型僅為丘陵山區(qū)，不同地貌下的耕地經(jīng)營(yíng)細(xì)碎化影響因素還有待研究。