烏日罕，白文龍，張裕鳳

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

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

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

1.1 研究區(qū)概況

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

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

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

2 研究方法

2.1 空間自相關(guān)

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

式中：n為觀測次數(shù);xi和xj分別為位置i和j處的流轉(zhuǎn)價格;x為流轉(zhuǎn)價格的均值;wij為價格i與價格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)，絕對值越大則表示相關(guān)性越強(qiáng)。

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

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

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

2.2 克里金插值

該方法是一種依賴于測量誤差模型實(shí)現(xiàn)精確或平滑插值的插值法。克里金法使用的是統(tǒng)計模型，其公式為：

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

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

3 結(jié)果與分析

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

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

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

表1 全局Moran’s Ⅰ值

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

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

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

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

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

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

3.2.1 空間趨勢分析

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

3.2.2 空間分異特征分析

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

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

4 結(jié)論

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

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

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