烏日罕，白文龍，張?jiān)ｘP

(內(nèi)蒙古師范大學(xué)地理科學(xué)學(xué)院，內(nèi)蒙古 呼和浩特 010020)

草原承包經(jīng)營(yíng)制度，是我國(guó)農(nóng)業(yè)經(jīng)營(yíng)制度的重要組成部分，是黨在牧區(qū)政策的基石，是建設(shè)現(xiàn)代化草原畜牧業(yè)的根。2014年中央“一號(hào)文件”首次提出“在落實(shí)農(nóng)村土地集體所有權(quán)的基礎(chǔ)上，穩(wěn)定農(nóng)戶承包權(quán)，放活土地經(jīng)營(yíng)權(quán)，允許承包土地的經(jīng)營(yíng)權(quán)向金融機(jī)構(gòu)抵押融資”〔1〕，三權(quán)分置隨即被視為是現(xiàn)階段我國(guó)農(nóng)地產(chǎn)權(quán)制度變革的基本方向〔2〕。在鄉(xiāng)村振興、農(nóng)村土地制度改革等大背景下，刺激了牧草地流轉(zhuǎn)相關(guān)方面的研究，有關(guān)草地流轉(zhuǎn)的文獻(xiàn)越來(lái)越多，草原產(chǎn)權(quán)方面〔3-5〕、草地規(guī)模經(jīng)營(yíng)〔6-9〕、流轉(zhuǎn)現(xiàn)狀〔10〕、流轉(zhuǎn)意愿〔11,12〕、流轉(zhuǎn)價(jià)格影響因素〔13,14〕等方面居多，但均缺少對(duì)牧草地流轉(zhuǎn)價(jià)格和空間位置之間存在的相關(guān)性的檢測(cè),缺少對(duì)牧草地流轉(zhuǎn)價(jià)格分布的空間關(guān)聯(lián)性的描述。本文以2018～2019年入戶調(diào)研數(shù)據(jù)為基礎(chǔ)進(jìn)行西烏珠穆沁旗草地流轉(zhuǎn)價(jià)格研究，不僅可以為宏觀把控草地流轉(zhuǎn)價(jià)格總體情況提供參考，而且對(duì)于科學(xué)指導(dǎo)當(dāng)前牧草地流轉(zhuǎn)市場(chǎng)的規(guī)范化有重要意義。

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

1.1 研究區(qū)概況

西烏珠穆沁旗(下文簡(jiǎn)稱西烏旗)，是華北、東北結(jié)合部(北緯43°52′～45°23′，東經(jīng)116°21′～119°23′)，海拔835～1957米之間，位于內(nèi)蒙古自治區(qū)中部，錫林郭勒盟東部，草原總面積22434.5km2。植被類型以草甸草原、典型草原為主，草原生產(chǎn)力相對(duì)較高。無(wú)論從草地面積、草地類型和牧業(yè)生產(chǎn)的重要性方面，西烏旗都是錫林郭勒盟最為典型的草原牧區(qū)，因此，選擇西烏旗作為研究區(qū)域具有很強(qiáng)的代表性。

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

數(shù)據(jù)來(lái)自2018～2019年對(duì)西烏旗牧區(qū)的實(shí)地調(diào)研，調(diào)查選取分層隨機(jī)抽樣方式，選取6個(gè)鎮(zhèn)(蘇木)：巴彥花鎮(zhèn)、浩勒?qǐng)D高勒鎮(zhèn)、吉仁高勒鎮(zhèn)、烏蘭哈拉嘎蘇木、巴彥胡碩蘇木、高日罕鎮(zhèn)。從每個(gè)鎮(zhèn)(蘇木)抽取6～9個(gè)嘎查，每個(gè)嘎查隨機(jī)挑選4～10個(gè)牧戶進(jìn)行問(wèn)卷調(diào)查。所得適合的有效問(wèn)卷178份。通過(guò)ArcGIS軟件對(duì)樣本數(shù)據(jù)進(jìn)行錄入和數(shù)字化處理以及正態(tài)分布檢驗(yàn)，建立西烏旗牧草地流轉(zhuǎn)樣點(diǎn)的空間地理數(shù)據(jù)庫(kù)(見(jiàn)圖1)。

2 研究方法

2.1 空間自相關(guān)

該方法是檢驗(yàn)具有空間位置關(guān)系的某要素的觀測(cè)值是否顯著地與其相鄰空間點(diǎn)上的觀測(cè)值相關(guān)聯(lián)〔15〕。全局自相關(guān)表達(dá)式為：

式中：n為觀測(cè)次數(shù);xi和xj分別為位置i和j處的流轉(zhuǎn)價(jià)格;x為流轉(zhuǎn)價(jià)格的均值;wij為價(jià)格i與價(jià)格j之間的空間權(quán)重，其中，如果觀察到i與j是空間相鄰的wij，則j=i，否則wij=0。I為Moran指數(shù)，在〔-1，1〕取值，大于0、等于0和小于0分別表示正相關(guān)、不相關(guān)和負(fù)相關(guān)，絕對(duì)值越大則表示相關(guān)性越強(qiáng)。

局域空間自相關(guān)表達(dá)式為：

其各個(gè)指數(shù)均與公式(1)中的一樣，Z得分由Moran’I系數(shù)、期望值和方差三部分組成，計(jì)算公式為：

式中，Ii為Moran指數(shù)；E(Ii)為期望；Var(Ii)為方差。

2.2 克里金插值

該方法是一種依賴于測(cè)量誤差模型實(shí)現(xiàn)精確或平滑插值的插值法?？死锝鸱ㄊ褂玫氖墙y(tǒng)計(jì)模型，其公式為：

Z(x0)=∑n(i=1)λiZ(xi) (公式3)

式中：Z(x0)為未知樣點(diǎn)草地流轉(zhuǎn)價(jià)格；Z(xi)為已知樣點(diǎn)草地流轉(zhuǎn)價(jià)格；λi為第i個(gè)樣點(diǎn)草地流轉(zhuǎn)價(jià)格的權(quán)重，n為樣點(diǎn)草地流轉(zhuǎn)價(jià)格的樣本數(shù)。

3 結(jié)果與分析

3.1 牧草地流轉(zhuǎn)價(jià)格空間自相關(guān)分析

3.1.1 全局空間自相關(guān)分析

在進(jìn)行全局空間自相關(guān)分析之前，首先要構(gòu)建空間權(quán)重矩陣。筆者選擇生成4個(gè)鄰居的權(quán)重矩陣，西烏旗草地流轉(zhuǎn)價(jià)格樣點(diǎn)數(shù)據(jù)的全局自相關(guān)Moran’sⅠ的計(jì)算結(jié)果如表1所示：

表1 全局Moran’s Ⅰ值

西烏旗牧草地流轉(zhuǎn)價(jià)的Moran’sⅠ指數(shù)為0.5778，標(biāo)準(zhǔn)化Z統(tǒng)計(jì)量為12.5613，遠(yuǎn)大于99%置信區(qū)間的檢驗(yàn)閾值，說(shuō)明牧草地流轉(zhuǎn)價(jià)格在空間上具有較強(qiáng)的正相關(guān)性，呈現(xiàn)出各鎮(zhèn)(蘇木)牧草地流轉(zhuǎn)價(jià)格高的地區(qū)與低的地區(qū)相對(duì)聚集的狀態(tài)，并存在一定的差異。

3.1.2 局部空間自相關(guān)分析

通過(guò)Geoda軟件將研究樣本點(diǎn)構(gòu)成面，然后構(gòu)建面的Rook權(quán)重矩陣，計(jì)算局部Moran’sI數(shù)值為0.5299。Moran散點(diǎn)圖大部分樣點(diǎn)分布在一三象限，說(shuō)明西烏旗草地流轉(zhuǎn)價(jià)格在局域空間上相關(guān)性程度很高，空間差異較小。少量樣點(diǎn)分布在二三象限，體現(xiàn)了流轉(zhuǎn)價(jià)格樣點(diǎn)的局部差異性。

由圖2可知，上文所述的集聚和異質(zhì)性特征在牧草地流轉(zhuǎn)價(jià)格的LISA圖顯示，即在浩勒?qǐng)D高勒鎮(zhèn)中部、巴彥花鎮(zhèn)中部以及高日罕鎮(zhèn)的西北部一帶是牧草地流轉(zhuǎn)價(jià)格的“高-高”區(qū)域，集聚比較明顯，此部分高日罕鎮(zhèn)位于烏拉蓋盆地，內(nèi)有烏拉蓋水系并且全部集中于此而潛入地下，形成大小不等的湖泊(淖爾)和大片沼澤濕地，使得這片草原水草豐美，草地流轉(zhuǎn)價(jià)格相對(duì)較高。浩勒?qǐng)D高勒鎮(zhèn)的阿拉坦高勒嘎查人均草地少，因此流轉(zhuǎn)價(jià)格偏高。

與此相反，分布在巴彥花鎮(zhèn)的東北部、浩勒?qǐng)D高勒鎮(zhèn)的東部大部分地區(qū)是草地流轉(zhuǎn)價(jià)格“低-低”集聚地區(qū)，該地區(qū)屬于大興安嶺西麓低山丘陵〔16〕，低山丘陵地區(qū)不適于作為打草場(chǎng)，牧民一般將其作為夏季放牧場(chǎng)，因此草地流轉(zhuǎn)價(jià)格可能會(huì)降低；但同時(shí)山坡陽(yáng)坡的草地長(zhǎng)勢(shì)更好，個(gè)別山地內(nèi)的草場(chǎng)流轉(zhuǎn)價(jià)格有可能偏高。值得一提的是，巴彥胡碩蘇木南部有一小部分也呈“低-低”集聚，這是由于該地區(qū)處于嘎亥額勒蘇沙地，因此該地區(qū)草地質(zhì)量偏低，影響牧草地流轉(zhuǎn)價(jià)格偏低〔17〕。

3.2 草地流轉(zhuǎn)價(jià)格空間分異規(guī)律

3.2.1 空間趨勢(shì)分析

為了研究牧草地流轉(zhuǎn)數(shù)據(jù)樣點(diǎn)的空間總體分布趨勢(shì)，利用ArcGIS的三維透視工具對(duì)數(shù)據(jù)進(jìn)行分析，如圖3所示，研究區(qū)范圍內(nèi)牧草地流轉(zhuǎn)價(jià)格從東西方向看，呈“倒U”型，流轉(zhuǎn)價(jià)格在橫向上由中心向兩邊逐漸降低，東部較西部?jī)r(jià)格較高。而南北方向上變化趨勢(shì)呈現(xiàn)輕微的由高降低，整體呈下降趨勢(shì)緩慢。草地流轉(zhuǎn)價(jià)格總體上是北部偏高，南部偏低，且南北走向的曲線在邊緣處出現(xiàn)較明顯的不平衡性，草地流轉(zhuǎn)價(jià)格呈現(xiàn)南高北低的不平衡特點(diǎn)。

3.2.2 空間分異特征分析

基于變異函數(shù)的理論模型，對(duì)樣點(diǎn)數(shù)據(jù)進(jìn)行了簡(jiǎn)單克里金插值，得到西烏旗草地流轉(zhuǎn)價(jià)格空間分布圖(見(jiàn)圖4)。由圖4可知，對(duì)插值結(jié)果進(jìn)行交叉檢驗(yàn)，標(biāo)準(zhǔn)平均值為0.0131，接近于0；標(biāo)準(zhǔn)均方根預(yù)測(cè)誤差為0.9724，接近于1，平均標(biāo)準(zhǔn)誤差為43.4815，接近于均方根42.5886，檢驗(yàn)結(jié)果表明模型預(yù)測(cè)效果較為理想。

由圖4得出西烏旗草地流轉(zhuǎn)的空間分異規(guī)律如下：(1)總體上西烏旗草地流轉(zhuǎn)價(jià)格的空間分布特征是由旗北部、東南部向外圍逐漸遞減，并呈現(xiàn)出東高西低的趨勢(shì)。北部盆地和東南部山地丘陵地區(qū)的地形地貌使其成為牧草地流轉(zhuǎn)價(jià)格的高值區(qū)。(2)草地流轉(zhuǎn)價(jià)格的幅度存在區(qū)域差異。東南部草地流轉(zhuǎn)價(jià)格等值線較密集，表明流轉(zhuǎn)價(jià)格衰減速度快，隨著距離的增大，流轉(zhuǎn)價(jià)格等值線也逐漸變得稀疏。(3)西烏旗東南部及北部小部分區(qū)域出現(xiàn)流轉(zhuǎn)價(jià)格次高峰，浩勒?qǐng)D高勒鎮(zhèn)和巴彥花鎮(zhèn)出現(xiàn)了明顯的流轉(zhuǎn)價(jià)格突起。

4 結(jié)論

西烏旗牧草地流轉(zhuǎn)價(jià)格存在顯著的全局空間自相關(guān)特征，全局空間自相關(guān)系數(shù)為0.5778，表明牧草地流轉(zhuǎn)價(jià)格具有空間相關(guān)性與集聚特征。局部空間自相關(guān)分析表明：西烏旗牧草地流轉(zhuǎn)價(jià)格空間關(guān)聯(lián)性顯著，空間異質(zhì)性較弱，高值聚集地區(qū)集中在東南部低山丘陵區(qū)以及北部盆地區(qū)域，在東北小部分鹽堿地以及旗中部沙地出現(xiàn)一定的空間異質(zhì)性。

總體來(lái)看，西烏旗東部牧草地流轉(zhuǎn)價(jià)格較西部偏高，北部較南部偏高，且南北走向的曲線在邊緣處出現(xiàn)較明顯的不平衡性，牧草地流轉(zhuǎn)價(jià)格呈現(xiàn)南高北低的不平衡特點(diǎn)。

通過(guò)克里金插值方法建立了數(shù)字牧草地流轉(zhuǎn)價(jià)格模型，結(jié)果表明：牧草地流轉(zhuǎn)價(jià)格存在區(qū)域關(guān)聯(lián)性與空間向心性，浩勒?qǐng)D高勒鎮(zhèn)和巴彥花鎮(zhèn)位于高價(jià)區(qū)；空間分布存在差異性和變異性，具有明顯的區(qū)域性，東南部和北部區(qū)域高于西部和東北部區(qū)域。