内蒙古西乌珠穆沁旗草地流转价格的空间分布格局

乌日罕，白文龙，张裕凤

(内蒙古师范大学地理科学学院，内蒙古 呼和浩特 010020)

草原承包经营制度，是我国农业经营制度的重要组成部分，是党在牧区政策的基石，是建设现代化草原畜牧业的根。2014年中央“一号文件”首次提出“在落实农村土地集体所有权的基础上，稳定农户承包权，放活土地经营权，允许承包土地的经营权向金融机构抵押融资”〔1〕，三权分置随即被视为是现阶段我国农地产权制度变革的基本方向〔2〕。在乡村振兴、农村土地制度改革等大背景下，刺激了牧草地流转相关方面的研究，有关草地流转的文献越来越多，草原产权方面〔3-5〕、草地规模经营〔6-9〕、流转现状〔10〕、流转意愿〔11,12〕、流转价格影响因素〔13,14〕等方面居多，但均缺少对牧草地流转价格和空间位置之间存在的相关性的检测,缺少对牧草地流转价格分布的空间关联性的描述。本文以2018～2019年入户调研数据为基础进行西乌珠穆沁旗草地流转价格研究，不仅可以为宏观把控草地流转价格总体情况提供参考，而且对于科学指导当前牧草地流转市场的规范化有重要意义。

1 研究区概况与数据来源

1.1 研究区概况

西乌珠穆沁旗(下文简称西乌旗)，是华北、东北结合部(北纬43°52′～45°23′，东经116°21′～119°23′)，海拔835～1957米之间，位于内蒙古自治区中部，锡林郭勒盟东部，草原总面积22434.5km2。植被类型以草甸草原、典型草原为主，草原生产力相对较高。无论从草地面积、草地类型和牧业生产的重要性方面，西乌旗都是锡林郭勒盟最为典型的草原牧区，因此，选择西乌旗作为研究区域具有很强的代表性。

1.2 数据来源与处理

数据来自2018～2019年对西乌旗牧区的实地调研，调查选取分层随机抽样方式，选取6个镇(苏木)：巴彦花镇、浩勒图高勒镇、吉仁高勒镇、乌兰哈拉嘎苏木、巴彦胡硕苏木、高日罕镇。从每个镇(苏木)抽取6～9个嘎查，每个嘎查随机挑选4～10个牧户进行问卷调查。所得适合的有效问卷178份。通过ArcGIS软件对样本数据进行录入和数字化处理以及正态分布检验，建立西乌旗牧草地流转样点的空间地理数据库(见图1)。

2 研究方法

2.1 空间自相关

该方法是检验具有空间位置关系的某要素的观测值是否显著地与其相邻空间点上的观测值相关联〔15〕。全局自相关表达式为：

式中：n为观测次数;xi和xj分别为位置i和j处的流转价格;x为流转价格的均值;wij为价格i与价格j之间的空间权重，其中，如果观察到i与j是空间相邻的wij，则j=i，否则wij=0。I为Moran指数，在〔-1，1〕取值，大于0、等于0和小于0分别表示正相关、不相关和负相关，绝对值越大则表示相关性越强。

局域空间自相关表达式为：

其各个指数均与公式(1)中的一样，Z得分由Moran’I系数、期望值和方差三部分组成，计算公式为：

式中，Ii为Moran指数；E(Ii)为期望；Var(Ii)为方差。

2.2 克里金插值

该方法是一种依赖于测量误差模型实现精确或平滑插值的插值法。克里金法使用的是统计模型，其公式为：

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

式中：Z(x0)为未知样点草地流转价格；Z(xi)为已知样点草地流转价格；λi为第i个样点草地流转价格的权重，n为样点草地流转价格的样本数。

3 结果与分析

3.1 牧草地流转价格空间自相关分析

3.1.1 全局空间自相关分析

在进行全局空间自相关分析之前，首先要构建空间权重矩阵。笔者选择生成4个邻居的权重矩阵，西乌旗草地流转价格样点数据的全局自相关Moran’sⅠ的计算结果如表1所示：

表1 全局Moran’s Ⅰ值

西乌旗牧草地流转价的Moran’sⅠ指数为0.5778，标准化Z统计量为12.5613，远大于99%置信区间的检验阈值，说明牧草地流转价格在空间上具有较强的正相关性，呈现出各镇(苏木)牧草地流转价格高的地区与低的地区相对聚集的状态，并存在一定的差异。

3.1.2 局部空间自相关分析

通过Geoda软件将研究样本点构成面，然后构建面的Rook权重矩阵，计算局部Moran’sI数值为0.5299。Moran散点图大部分样点分布在一三象限，说明西乌旗草地流转价格在局域空间上相关性程度很高，空间差异较小。少量样点分布在二三象限，体现了流转价格样点的局部差异性。

由图2可知，上文所述的集聚和异质性特征在牧草地流转价格的LISA图显示，即在浩勒图高勒镇中部、巴彦花镇中部以及高日罕镇的西北部一带是牧草地流转价格的“高-高”区域，集聚比较明显，此部分高日罕镇位于乌拉盖盆地，内有乌拉盖水系并且全部集中于此而潜入地下，形成大小不等的湖泊(淖尔)和大片沼泽湿地，使得这片草原水草丰美，草地流转价格相对较高。浩勒图高勒镇的阿拉坦高勒嘎查人均草地少，因此流转价格偏高。

与此相反，分布在巴彦花镇的东北部、浩勒图高勒镇的东部大部分地区是草地流转价格“低-低”集聚地区，该地区属于大兴安岭西麓低山丘陵〔16〕，低山丘陵地区不适于作为打草场，牧民一般将其作为夏季放牧场，因此草地流转价格可能会降低；但同时山坡阳坡的草地长势更好，个别山地内的草场流转价格有可能偏高。值得一提的是，巴彦胡硕苏木南部有一小部分也呈“低-低”集聚，这是由于该地区处于嘎亥额勒苏沙地，因此该地区草地质量偏低，影响牧草地流转价格偏低〔17〕。

3.2 草地流转价格空间分异规律

3.2.1 空间趋势分析

为了研究牧草地流转数据样点的空间总体分布趋势，利用ArcGIS的三维透视工具对数据进行分析，如图3所示，研究区范围内牧草地流转价格从东西方向看，呈“倒U”型，流转价格在横向上由中心向两边逐渐降低，东部较西部价格较高。而南北方向上变化趋势呈现轻微的由高降低，整体呈下降趋势缓慢。草地流转价格总体上是北部偏高，南部偏低，且南北走向的曲线在边缘处出现较明显的不平衡性，草地流转价格呈现南高北低的不平衡特点。

3.2.2 空间分异特征分析

基于变异函数的理论模型，对样点数据进行了简单克里金插值，得到西乌旗草地流转价格空间分布图(见图4)。由图4可知，对插值结果进行交叉检验，标准平均值为0.0131，接近于0；标准均方根预测误差为0.9724，接近于1，平均标准误差为43.4815，接近于均方根42.5886，检验结果表明模型预测效果较为理想。

由图4得出西乌旗草地流转的空间分异规律如下：(1)总体上西乌旗草地流转价格的空间分布特征是由旗北部、东南部向外围逐渐递减，并呈现出东高西低的趋势。北部盆地和东南部山地丘陵地区的地形地貌使其成为牧草地流转价格的高值区。(2)草地流转价格的幅度存在区域差异。东南部草地流转价格等值线较密集，表明流转价格衰减速度快，随着距离的增大，流转价格等值线也逐渐变得稀疏。(3)西乌旗东南部及北部小部分区域出现流转价格次高峰，浩勒图高勒镇和巴彦花镇出现了明显的流转价格突起。

4 结论

西乌旗牧草地流转价格存在显著的全局空间自相关特征，全局空间自相关系数为0.5778，表明牧草地流转价格具有空间相关性与集聚特征。局部空间自相关分析表明：西乌旗牧草地流转价格空间关联性显著，空间异质性较弱，高值聚集地区集中在东南部低山丘陵区以及北部盆地区域，在东北小部分盐碱地以及旗中部沙地出现一定的空间异质性。

总体来看，西乌旗东部牧草地流转价格较西部偏高，北部较南部偏高，且南北走向的曲线在边缘处出现较明显的不平衡性，牧草地流转价格呈现南高北低的不平衡特点。

通过克里金插值方法建立了数字牧草地流转价格模型，结果表明：牧草地流转价格存在区域关联性与空间向心性，浩勒图高勒镇和巴彦花镇位于高价区；空间分布存在差异性和变异性，具有明显的区域性，东南部和北部区域高于西部和东北部区域。