1970~2010年云南农田生产潜力时空演变特征及影响因素分析

Spatial-Temporal Evolution Characteristics and Influencing Factors Analysis of Farmland Production Potential in Yunnan from 1970 to 2010

摘要: 云南省区域内部地貌立体地形复杂多样，独特的地理条件造就了多样化自然条件组合，农业环境中不同组分的改变会对粮食生产系统产生影响，现有农田生产潜力计算模型估算结果总量变化经验证能够拟合实际粮食产量。本文以云南省的129个县作为基本研究单元，分析了1970年至2010年间每十年一个时间节点下的农田生产潜力空间格局的变化特征，为定义出农田生产潜力增产限制因素，选择年平均气温、降水、耕地面积、地形起伏度作为影响因素参与分析，增加县域人口总数、GDP总量作为社会经济要素分析是否与农田生产潜力存在相关关系。使用ArcGIS 10.2中栅格计算、重心分析、莫兰指数、热点分析、地理加权回归分析等工具进行实验操作，实验结果表明：1970至2010年50年间云南省农田生产潜力总体呈逐渐增长趋势，高值增长区县数量集中连片分布于云南省中部；在研究时段内农田生产潜力重心移动速度较慢波动幅度小，说明农田生产潜力空间分布格局较稳定；农田生产潜力高值区多集中于耕地面积资源丰富地形起伏度小的地区且人口数量与农田生产潜力也存在相关关系。

Abstract: The topography within the Yunnan Province is characterized by complex and diverse terrain, creating unique geographical conditions and a diverse combination of natural elements. Changes in different components of the agricultural environment can significantly impact the grain production system. The results of the existing farmland production potential calculation model have been validated to effectively estimate actual grain production. This study takes the 129 counties of Yunnan Province as the basic research units and analyzes the spatial-temporal evolution characteristics of farmland Production potential from 1970 to 2010, with each decade as a time node. In order to define the limiting factors of increasing production potential of farmland, annual average temperature, precipitation, cultivated land area and topographic relief were selected as influencing factors to participate in the analysis, and the total population and GDP of the county were also increased as social and economic factors to analyze whether there is a correlation with farmland production potential. The grid calculation, gravity center analysis, Moran index, hot spot analysis, Geographic Weighted Regression Analysis and other tools in ArcGIS10.2 were used for experimental operation. The experimental results showed that the farmland production potential of Yunnan Province showed a gradual growth trend in the 50 years from 1970 to 2010, the number of districts and counties with high value growth was concentrated in the central part of Yunnan Province. During the study period, the center of gravity of farmland production potential moves slowly and fluctuates less, indicating that the spatial distribution pattern of farmland production potential is stable. The high value areas of farmland production potential are mainly concentrated in the areas with rich arable land resources and small topographic relief, and there was also a correlation between population quantity and farmland production potential.

1. 引言

粮食生产是世界各国政府重点关注的发展战略问题，粮食生产系统是多因素共同作用形成的综合系统，受气候、土壤等自然因素以及科技、劳动力等社会经济要素的影响，粮食生产的时空特征存在差异 [1] 。农田生产潜力也被叫农田作物生产潜力，是指现有耕作技术和生产管理水平下农田能持续生产的最大生产能力，受光照、热量、水分等农业气候条件以及农田的数量和质量的影响 [2] 。研究表明，在现有条件下，农田生产潜力和实际粮食产量之间仍存在差距，生产潜力尚未得到充分挖掘 [3] ，对农田生产潜力进行研究可以快速获取研究区内部长时间尺度变化和区域分异规律，结合对应年份各县域范围内自然与社会经济统计数据分析农田生产潜力在不同区域产生差异的原因，定义出关键限制因素，从而为农业空间布局、以及种植规划提供决策依据。

国内外针对农业生产潜力时空格局及影响因素分析都进行了较为深入探讨。

国外大尺度范围研究区通常采用逐级限制法建立农田生产潜力计算模型，以验证分析农作物生长潜力。A. Verdoodt [4] 研究得出在南非半干旱地区玉米和向日葵试验种植区由于对降雨有严重依赖，干旱对作物生产产量产生毁灭性影响；M. V. K. Sivakumar [5] 使用气象数据作为约束条件，整合评估作物生产潜力，为可持续农业发展提供规划依据。外国学者认为气候、土壤等综合因素在不同区域不同组分条件下影响了农作物生育期的长短进而影响了农作物生产潜力的增减。A. Sattar [6] 研究了7个长期干旱易发区的降雨资料、潜在蒸散量、有效持水量在不同土质下对农作物生育期长短的影响。除了自然因素影响分析，也有学者研究了造成近年来农田生产潜力下降的因素与农业管理制度和由于快速城市化大量面积被占用有关。Gergely Toth [7] 分析土地利用变化数据得出结论：欧盟各成员国农田生产潜力下降的原因是伴随着人口增多、经济发展耕地不断转换为人工地表。J. Timsina [8] 等从水稻和小麦的播种管理制度、播种时间角度分析影响因素与产量增长之间的关系。

国内学者对农田生产潜力研究包括从全国范围按照省界或分区研究，以及针对县域小尺度空间单元进行差异性研究。以全国范围尺度研究可以宏观把握粮食生产格局，分析地区间时间和空间上的差异，对比潜在生产力与实际粮食产量，分析粮食潜力开发水平差异。贾琨 [9] 根据“一带一路”沿线包括中国在内的65个国家40年统计年鉴数据，分析粮食产量、面积、单产等数据研究从全区域到每个国别时序演变特征及空间分布格局；韩荣青 [10] 比较实际粮食产量和潜在粮食生产力的变化格局；张锦宗 [11] 从省际范围分析我国近30年以来粮食生产潜力格局呈南升北降分布格局是各省份播种面积之间存在差异而产生的结果。有学者认为对粮食生产时空格局和影响因子分析在省域范围研究不够充分，无法体现县域内部差异性，有研究从县域范围对粮食生产格局增减情况分析研究，使用不同工具分析具体县域范围内实际产量与计算潜力之间的差异造成因素。刘正兵 [12] 使用空间常系数回归、空间变系数回归分析山西省县域范围粮食生产格局影响因子进行实证分析，研究不同影响要素作用于粮食潜力生产格局的效应；杜国明 [13] 计算近五十年三江平原粮食生产潜力验证计算结果与实际各县粮食产量统计数据的拟合程度；屈艳辉 [14] 考虑到各市周边县域之间的距离，使用地理加权回归方法分析不同时间范围几种对县域粮食产量格局在不同县域位置的空间变异。国内学者多采用逐级修订法计算农田生产潜力，气候因子(气温、降水、辐射)对农业影响深刻，且其成分变化对区域内不同作物气候生产潜力影响有差异，各种农业限制因素相互协调程度影响农业资源优劣判断。庞艳梅 [15] 使用逐级订正模型使用计算四川盆地1961至2018年几种主要粮食作物气候生产潜力；冷疏影 [16] 使用地理信息系统计算农业生产潜力并进行等级分区；王情 [17] 选择七种淮河流域主要粮食计算粮食生产潜力转换率，将粮食生产潜力转换区高的地区划分为开发适宜性高地区。

2. 研究区与数据来源

2.1. 研究区概况

Figure 1. Location of study area

图1. 研究区地理位置

云南省地处低纬，北回归线从中部穿过，区域内海拔高差较大，西北部横断山区山高谷深，地形复杂、地貌立体是典型的喀斯特地貌区，气候垂直差异显著(图1)。云南省光热资源丰富，雨热同期，生物气候条件多种多样，适合发展林牧业和花卉、烟草等产业，海陆热力性质差异影响下的季风气候导致降水时空分布不均，地区有效利用的农业资源匮乏。合理利用农业气候资源，要做到因地制宜，避免和减轻不利农业生产环境条件的影响，完善种植结构和布局。

2.2. 数据来源

研究使用的数据包括：1970、1980、1990、2000、2010五个年份的中国农田生产潜力数据、年平均气温和降水、GDP、人口总数、中国90米DEM数据、云南省区县矢量数据。

中国农田生产潜力数据是基于中国耕地遥感监测空间分布、土壤和高程数据，采用GAEZ (Global Agro-Ecological Zones)模型，综合考虑气象、土壤等多个方面的因素，估算得到中国全国耕地的生产潜力。数据集主要考虑了小麦、玉米、水稻、大豆和甘薯五种中国主要种植作物，由于我国幅员辽阔，地形地貌复杂，气候类型多样，地域位置不同作物熟制不同，北方温带季风区、东北地区、华北地区农业制度多为一年一熟到两年三熟，南方地区水热条件较北方优越可达一年两熟到三熟，因此计算土地粮食生产潜力需考虑作物熟制。计算土地粮食生产潜力前先根据气温、降水等气候条件估算某种作物空间分布的气候适宜性，对适宜作物种植的地区使用逐级限制法对气温、降水、辐射三个气象条件限制得到气候生产力，依据土壤类型、肥力等农业属性添加限制，得到土地生产潜力，最后对农业耕作管理制度、土地利用方式添加限制得出土地粮食生产潜力，数据集单位为kg/hm²。5个数据集目前都已经公开出版，质量状况良好 [18] 。

实验使用中国科学院地理科学与资源研究所提供的2008年全国县级行政边界矢量数据。云南省数字高程模型(DEM)数据通过地理空间数据云端下载，随后进行了镶嵌和裁剪处理。年平均气温、降水量、人口数量和GDP数据均来自对应年份《云南省统计年鉴》。

3. 研究方法

Figure 2. Technical roadmap

图2. 技术路线图

在本研究中，我们以云南省129个区县为基本研究单元利用ArcGIS10.2中的重心分析、莫兰指数、热点分析、地理加权回归等工具每十年一个时间节点分析了1970至2010年间农田生产潜力的时空格局的变化特征，选择了年平均气温、降水、耕地面积、地形起伏度、县域人口总数和GDP六个影响因子分析各因子对农田生产潜力的作用，为云南省科学合理粮食安全政策构建提供理论参考。

实验所使用的农田生产潜力数据是基于耕地的分布、土壤、DEM高程数据，使用逐级限制法计算我国五种主要粮食作物的粮食生产潜力 [18] 。选取参与影响因素分析的气象数据包括年平均气温、降水两个县域记录较完整的数据；短时间内土壤质地、肥力、侵蚀强度等变化较小，因此未选取土壤因素参与分析；耕地数据为各县域耕地面积，地形起伏度是反映范围内最高海拔与最低海拔的差值，相较于DEM数据地形起伏度更能描述出一个区域地貌形态，基于DEM数据使用ArcGIS空间统计工具计算可得出地形起伏度结果；除以上自然条件外还增加县域总人口数与GDP总量作为社会经济要素参与分析，分析农田生产潜力与社会经济要素之间是否有相关性。由于1970年云南省县域统计资料有缺失，所以影响因素分析的年份为1980年、1990年、2000年、2010年四个时间节点(图2)。

3.1. 重心分析法

“重心”分析借用重心点及其移动方向、移动距离等指标不仅可以刻画区域地理现象的空间差异，还可以进一步探寻其动态过程及演化规律通过可视化表示出重心点位置以及移动方向、距离，从而描述区域内地理现象的空间分布差异和时空格局变化 [19] 。实验计算了1970、1980、1990、2000、2010五个时间节点云南省农田生产潜力重心，描述出迁移路线，总结该区域相关时间段内农田生产潜力重心移动过程以及其移动规律。

3.2. 全局空间自相关

全局莫兰指数(Global Morans’I) [20] 研究全局空间范围内某位置和周围邻接或临近的空间区域单元属性值是否存在聚集性，反映了观测变量在整个研究区域内空间相关性的总体趋势 [21] 。全局莫兰指数计算公式为：

$I=\frac{n}{s_0}\frac{{{\displaystyle \sum }}_{i=1}^n{{\displaystyle \sum }}_{j=1}^n{w}_{ij}\left(y_i-\bar{y}\right)\left(y_j-\bar{y}\right)}{{{\displaystyle \sum }}_{i=1}^n{\left(y_i-\bar{y}\right)}^2}$ (1)

$S_0={\displaystyle {\sum }_{i=1}^n{\displaystyle {\sum }_{j=1}^n{w}_{ij}}}$ (2)

上式中：I为莫兰指数；n为研究区内空间单元总数；w_ij为空间权重矩阵；y_i为区域i的观测值。根据莫兰指数结果可以将描述对象空间分布模式划分为聚类、离散、随机分布。莫兰指数值的大小还反映区域单元相关的程度。

P值，P-value表示概率，是研究对象空间分布模式呈随机分布的概率，P值越小说明越不太可能为随机分布模式。Z值得分，Z-score是标准差的倍数，Z值分正负，当P值小于0.1且Z得分小于−1.65或大于1.65，置信度为90%；当P值小于0.05且Z值小于−1.96或大于1.96置信度为95%；当P值小于0.01且Z值小于−2.58或大于2.58置信度为99%。在县域尺度下计算农田生产潜力及各相关因子莫兰指数，分析各相关因子与农田生产潜力在空间上分布的差异程度，分析各因子是否对农田生产潜力的格局分布有影响。

3.3. 局部空间自相关

聚类和异常值分析的指标则用于反映区(县)地区观测值与其相邻区(县)之间的差异程度和差异的显著性 [22] 。局部莫兰指数的计算公式为：

$I_i=\frac{y_i-\bar{y}}{s^2}{\displaystyle {\sum }_{i\ne j}^n{w}_{ij}\left(y_j-\bar{y}\right)}$ (3)

$S^2=\frac{1}{n}{{\displaystyle \sum \left(y_i-\bar{y}\right)}}^2$ (4)

上式中：I_i为莫兰指数，结果取值范围为(−1, 1)，值小于0表示分析单元相邻单元间分析属性值为较弱的集聚，即离散模式；值大于0表示相邻各单元的属性值存在比较显著地集聚现象；趋近于0表示各单元呈随机分布。局部莫兰指数结果分为H-H (高高聚集)、H-L (高低异常)、L-L (低低聚集)、L-H (低高异常)四种类型属性值，以及不显著聚集类型。

3.4. 热点分析法

Getis-Ord G_i^*用于检测局部小范围内局部空间显著热点、冷点。每个单元i的统计量G_i为：

$G_i{}^{*}=\frac{{\displaystyle {\sum }_{j=1}^n{w}_i{}_{,j}{x}_j-\bar{x}}{\displaystyle {\sum }_{j=1}^n{w}_{ij}}}{\sqrt{\frac{{\displaystyle {\sum }_{j=1}^n{x}_j{}^2}}{n}-{\left(\bar{x}\right)}^2}\sqrt{\frac{\left[n{{\displaystyle {\sum }_{j=1}^n{w}_i{{}_{,j}}^2-\left({\displaystyle {\sum }_{j=1}^n{w}_i{}_{,j}}\right)}}^2\right]}{n-1}}}$ (5)

上式中，x_j为第j个地理单元的密度值；x为所有单元的密度均值；w_i,j为单元i和j的空间权重系数；n为研究区域内所有单元的总数量。热点分析结果为显著正G_i值表示在观测值较高的区域邻近单元也具有较高观测值，显著的负G_i值表示在观测值较低的区域邻近单元也具有较低观测值，与全局莫兰指数相比，热点分析结果能反映出高值和低值聚类具体在空间中的分布 [23] 。县域尺度下对农田生产潜力和影响因素进行探测，分析是否存在统计上显著相关的高值集聚区和低值集聚区。

3.5. 地理加权回归分析法

传统的线性回归模型是估计参数的“平均”或“全局”意义，如果自变量是空间数据且存在空间自相关性，就不适用传统回归模型(Ordinary Least Squares, OLS)残差项独立的假设 [24] ，地理加权回归(Geographically Weighted Regression, GWR)容许一些不平稳的数据直接被模拟，利用非参数估计方法进行局部参数估计。通过线性回归模型中假定回归系数是观测点地理位置的位置函数，将数据空间特性纳入模型中，为分析回归关系的空间差异性创造了条件，能够反映参数在不同空间的空间非平稳性，使变量间的关系可以随空间位置的变化而变化 [25] 。GWR模型表达式为：

$y_i={\beta }_0\left(u_i,v_i\right)+{\displaystyle \underset{j=1}{\overset{k}{\sum }}{\beta }_j\left(u_i,v_i\right)}{x}_{ij}+{\epsilon }_i$ (6)

上式中，y_i表示地理位置 $\left(u_i,v_i\right)$ 的因变量， ${\beta }_0\left(u_i,v_i\right)$ 表示地理位置 $\left(u_i,v_i\right)$ 的常数项值， ${\beta }_j\left(u_i,v_i\right)$ 是在地理位置 $\left(u_i,v_i\right)$ 的未知参数，ε_i表示空间随机残差。本文将使用地理加权回归法分析各个因子对农田生产潜力的影响程度。

4. 结果与分析

4.1. 云南省农田生产潜力空间格局变化

使用区(县)矢量面数据对农田生产潜力数据集进行分区统计，分析结果得出1970年以来，云南省农田生产潜力总量整体为增长趋势(图3)，其中1980年到1990年呈下降趋势，1970年到2010年40年间云南省农田生产潜力总量增长了13.76%，从2994.36万t增加到3406.60万t；县域农田生产潜力均值从4395.37 kg/hm²增加到5252.53 kg/hm²，增长19.50%；最大县域农田生产潜力生产均值增长趋势与总量一致，1980年至1990年呈下降趋势，其余年份均为上升趋势。

县域农田生产潜力均值为使用分区统计各县域农田生产潜力均值与区县耕地面积之积。耕地是主要用于种植水稻、小麦等主要粮食作物以及蔬菜的土地，包括水田、望天田(雨育水田)、旱地、水浇地、几种常见类型。实验选择使用耕地面积数据计算各县域内农田生产潜力均值，使用耕地面积数据为国家青藏高原科学数据中心多年土地利用现状遥感监测数据集，通过面积制表工具计算出各区县耕地面积，经验证，使用该数据集计算所得出的2000年耕地面积是同一年份下《云南省统计年鉴》中统计耕地面积1.45倍，以云南省县域为单位进行数据比较如图4所示，虚线为云南省各区县农田生产潜力数据集计算得出的耕地面积与《统计年鉴》中实际耕地面积的相关趋势线，可从两者拟合情况中得出各区县计算得出的耕地面积与统计实际耕地面积相关系数为0.64，拟合情况较好。因此使用统计面积计算县域农田生产潜力均值能够反映实际情况变化趋势。

Figure 3. Statistics of potential production of farmland in Yunnan Province from 1970 to 2010

图3. 1970年~2010年云南省农田生产潜力统计图

Figure 4. Comparison of calculated and actual cultivated land area in county of Yunnan Province in 2000

图4. 2000年云南省县域计算耕地面积与实际耕地面积对比

使用栅格计算器计算每10年一个时间节点之间的县域农田生产潜力差值，分析1970至2010年40年间农田生产潜力数量变化的空间格局特征。

1970年至1980年(图5(a))，农田生产潜力增加的区(县)多于减少的区(县)，农田生产潜力增加超过50,000 t的区(县)共8个，主要集中在昭通市，包括永善县、镇雄县、宜良县、盐津县、威信县、和绥江县，另有昆明市下辖禄劝彝族苗族自治县，和楚雄彝族自治州禄丰县；农田生产潜力减少超过50,000 t的区(县)仅曲靖市罗平县，减少量在50,000~100,000 t的区(县)共9个，主要集中于德宏傣族景颇族自治州下辖的盈江县、芒市，保山市隆阳区、昌宁县、腾冲县，临沧市凤庆县、云县，文山壮族苗族自治州砚山县、丘北县。

1980年至1990年(图5(b))，农田生产潜力增加区(县)多于减少区(县)，农田生产潜力增加超过50,000 t的区(县)集中于文山市，包括砚山县、文山市2个区(县)；农田生产潜力减少超过50,000 t的区(县)仅楚雄彝族自治州元谋县，减少量在50,000~100,000 t的区(县)共14个，主要集中于昆明市包括西山区、晋宁县、禄劝彝族苗族自治县、安宁市，曲靖市陆良县、沾益县，楚雄彝族自治州楚雄市、武定县、禄丰县，昭通市永善县，丽江市永胜县、迪庆藏族自治州及香格里拉市，玉溪市下辖新平彝族傣族自治县。

1999年至2000年(图5(c))，农田生产潜力增加超过50,000 t的区(县)共22个，主要集中于昆明市，包括晋宁县、安宁市，玉溪市华宁县、易门县、新平县、元江哈尼族彝族傣族自治县，普洱市下辖思茅区、镇沅彝族哈尼族拉祜族自治县、墨江哈尼族自治县、景谷傣族彝族自治县、宁洱哈尼族彝族自治县，楚雄彝族自治州楚雄市、双柏县、元谋县、禄丰县，红河哈尼族彝族自治州石屏县、红河县，文山壮族苗族自治州马关县、丘北县、广南县，曲靖市罗平县，西双版纳州景洪市；农田生产潜力减少的区(县)仅有昭通市镇雄县。

Figure 5. Difference of county farmland production potential in Yunnan Province

图5. 云南省县域农田生产潜力差值

2000年至2010年(图5(d))，农田生产潜力增加的区(县)比减少的区(县)多，农田生产潜力增加超过50,000 t的区(县)共36个，主要集中于楚雄彝族自治州包括永仁县、双柏县、大姚县、楚雄市、元谋县、武定县、禄丰县，昆明市石林彝族自治县、禄劝县，曲靖市麒麟区、陆良县、会泽县、沾益县，玉溪市通海县、华宁县、新平彝族傣族自治县，临沧市临翔区、双江县、云县、凤庆县、耿马县，丽江市永胜县、宁蒗彝族自治县，红河哈尼族彝族自治州个旧市、蒙自市、建水县、弥勒市、泸西县，普洱市景洪市、勐海县，普洱市景谷傣族彝族自治县、镇沅县、景东彝族自治县、澜沧拉祜族自治县，大理白族自治州祥云县、宾川县；农田生产潜力减少超过50,000 t的区(县)仅文山壮族苗族自治州广南县，减少量在50,000~100,000 t的区(县) 5个，主要集中于昭通市，包括昭阳区、镇雄县、彝良县，和文山壮族苗族自治州砚山县。

1970年至2010年(图5(e))，农田生产潜力增加的区(县)多于减少的区(县)，农田生产潜力增加超过50,000 t的区(县)共33个，集中成片分布在云南省中部；农田生产潜力减少超过50,000 t的区(县)仅保山市隆阳区，减少量在50,000~100,000 t的区(县)共6个，主要集中于曲靖市，包括罗平县、富源县，保山市腾冲县，迪庆藏族自治州香格里拉市，德宏州下辖盈江县，空间位置分布无显著特征。

4.2. 云南省农田生产潜力重心时空格局演变

4.2.1. 云南省农田生产潜力重心区位及空间移动路径

基于区域空间重心计算1970年云南省农田生产潜力重心为(102˚04'29''E, 24˚59'04''N)，位于楚雄彝族自治州禄丰县；1980年农田生产潜力重心转变为(102˚12'27''E, 25˚02'35''N)，空间位置从禄丰县转移至昆明市下辖安宁市辖区西部；1990年农田生产潜力重心为(102˚12'54''E, 25˚00'00''N)，位于安宁市与禄丰县交界点；2000年重心位置移动至(102˚13'38''E, 24˚54'24''N)，位于玉溪市易门县境内；2010年重心位置移动至(102˚6'11''E, 24˚51'18''N)，仍位于玉溪市易门县境内。在研究时间段内，2010年重心终点位置位于1970年重心起始点位置的南偏东方向，移动路径在空间中近似于平行四边形如图6所示。

Figure 6. Transfer path of gravity center of farmland production potential in Yunnan Province

图6. 云南省农田生产潜力重心转移路径

4.2.2. 云南省农田生产潜力重心移动方向和距离

1970年至1980年重心迁移路径为西北走向，该阶段移动速度相对较快，年均移动1.45 km，自1980年后重心逐渐南移，1980年至2000年均为东南方向，1980年至1990年重心年均移动0.44 km，移动速度相对最慢，1990年至2000年移动速度增加，重心年均移动距离增加为0.94 km，2000年至2010年重心移动路径转变为西南走向且年均移动距离增加至1.36 km，云南省农田生产潜力重心移动距离较短速度较慢，农田生产潜力空间格局变化较稳定(表1)。

Table 1. Changes in gravity center of potential farmland production in Yunnan Province from 1970 to 2010

表1. 1970至2010年云南省农田生产潜力重心变动

4.3. 云南省农田生产潜力及影响因素空间自相关分析

4.3.1. 全局莫兰指数

结果如表2所示，1980年至2010年农田生产潜力及相关因素Morans’I值均为正值，除1980至2000年县域GDP外其余年份农田生产潜力以及各影响因素在99%的置信度下Z得分均大于2.58，因此可以拒绝零假设，认为农田生产潜力以及各影响因素在空间中随机分布的可能小于0.01，则在空间分布上表现出显著的集聚性，农田生产潜力以及各影响因素值较高的区域都趋于与值较高的区域相邻，值较低的值趋于与值较低的区域相邻，而1980年至2000年县域GDP则在99%的置信水平下在全局空间分布为随机分布模式的概率超过10%，不能拒绝0假设，在全局空间范围内无聚集并不一定在局部空间也没有集聚现象，所以GDP仍然参与后续的局部莫兰指数以及热点探测分析。

Table 2. Global Morans’I index of county farmland production potential and its influencing factors in Yunnan Province from 1980 to 2010

表2. 云南省1980~2010年县域农田生产潜力及各影响因素全局Morans’I指数

4.3.2. 热点探测

Figure 7. Hotspot analysis results of potential production and influencing factors of county farmland in Yunnan Province in 1980

图7. 1980年云南省县域农田生产潜力及各影响因素热点分析结果

冷热点探测结果显示，1980至2010年云南省县域农田生产潜力及各因素热点探测结果分布格局基本一致，县域范围内农田生产潜力热点分布区与耕地面积热点分布有重合部分，重合部分位于云南省东部曲靖市南部区县，农田生产潜力热点分布区也与地形起伏度冷点有重合部分，重合部分位于昆明市西部和曲靖市中部及南部；冷点分布区域与地形起伏度热点、耕地面积冷点、降水热点、人口总数冷点、GDP热点位置有重合部分(图7~10)。

Figure 8. Hotspot analysis results of potential production and influencing factors of county farmland in Yunnan Province in 1990

图8. 1990年云南省县域农田生产潜力及各影响因素热点分析结果

Figure 9. Hotspot analysis results of potential production and influencing factors of county farmland in Yunnan Province in 2000

图9. 2000年云南省县域农田生产潜力及各影响因素热点分析结果

Figure 10. Hotspot analysis results of potential production and influencing factors of county farmland in Yunnan Province in 2010

图10. 2010年云南省县域农田生产潜力及各影响因素热点分析结果

4.3.3. 局部莫兰指数

Figure 11. LISA cluster diagram of farmland production potential and its influencing factors in counties of Yunnan Province in 1980

图11. 1980年云南省县域农田生产潜力及各影响因素LISA集聚图

根据LISA集聚图进一步分析农田生产潜力以及各因子的局部空间关系。1980年农田生产潜力“高–高”聚集区共14个区县，主要聚集于云南省东部曲靖市中部下辖区县，其中文山壮族苗族自治州丘北县和曲靖市罗平县两个县所处区域也是耕地面积“高–高”聚集区，农田生产潜力“高–高”聚集区中曲靖市麒麟区也是GDP和区(县)人口总数的高聚集区，楚雄彝族自治州元谋县也是年平均降水“高–高”聚集区；农田生产潜力“低–低”聚集区为怒江傈僳族自治州福贡县及其邻县贡山独龙族自治县，贡山县是年平均降水、地形起伏度“高–高”聚集区；耕地面积、气温、人口“低–低”聚集区福贡县所处区域是地形起伏、年平均降水“高–高”聚集区，年平均气温、县域总人口的“低–低”聚集区(图11)。

1990年农田生产潜力“高–高”聚集区共12个区县，空间分布与1980年分布一致，文山壮族自治州丘北县为“低–高”聚集区，“低–高”聚集表示该县域周围存在农田生产潜力较低的县域，丘北县和罗平县两个县所处区域是耕地面积“高–高”聚集区，“高–高”聚集区中曲靖市麒麟区、宣威市、罗平县也是GDP和区(县)人口总数的高值聚集区；“低–低”聚集区相较于1980年新增加兰坪县共三个区县，三个区县所处区域也是地形起伏“高–高”聚集区，年平均气温“低–低”聚集区，贡山县、福贡县区范围也是县域总人口的“低–低”聚集区(图12)。

Figure 12. LISA cluster diagram of farmland production potential and its influencing factors in counties of Yunnan Province in 1990

图12. 1990年云南省县域农田生产潜力及各影响因素LISA集聚图

2000年农田生产潜力“高–高”聚集区共12个区县，与1980年分布基本一致，南部新增红河州开远市，“高–高”聚集区中曲靖市麒麟区、宣威市、富源县、罗平县也是区(县)人口总数的高值聚集区；“低–低”聚集区相较于1990年新增加维西傈僳族自治县、泸水县、大理白族自治州云龙县共六个区县，所处区域也是地形起伏度“高–高”聚集区，年平均气温低值聚集区，贡山县、福贡县、维西县区域也是县域总人口的低值聚集区(图13)。

2010年农田生产潜力“高–高”聚集区共12个区县，新增楚雄州元谋县及思茅市孟连傣族拉祜族佤族自治县，“高–高”聚集区与区(县)人口总数的高聚集区位置一致，其中元谋县是耕地面积数量统计高值集聚区；“低–低”聚集区相较于2000年新增加德钦县共七个区县，德钦县县域位置也是气温与人口分布的“低–低”集聚区、地形起伏变化的高值区(图14)。

Figure 13. LISA cluster diagram of farmland production potential and its influencing factors in counties of Yunnan Province in 2000

图13. 2000年云南省县域农田生产潜力及各影响因素LISA集聚图

Figure 14. LISA cluster diagram of farmland production potential and its influencing factors in counties of Yunnan Province in 2010

图14. 2010年云南省县域农田生产潜力及各影响因素LISA集聚图

4.3.4. 地理加权回归分析

如表3所示将耕地面积、地形起伏度、年平均气温、年平均降水、人口、县域GDP几个影响因子作为自变量z-score标准化处理后运用回归模型进行计算，构建OLS模型分析自变量对因变量的影响程度， 每个因子都通过了显著性检验，其中耕地面积与人口在研究时段内系数都为正向系数，说明两个因子对农田生产潜力有正向拉动力；地形起伏度与年平均降水在研究时段内系数都为负向系数，说明两个因子对农田生产潜力有负向拉动力，年平均气温OLS模型分析结果则是在1980年至1990年为负向系数，2000至2010年为正向系数；县域GDP的OLS模型分析结果则是在1980为正向系数，1990至2010年为负向系数。OLS模型虽然在一定程度上解释了影响因子的正负极性和影响力，但是其回归系数是建立在全局均一化的基础上的，不能辨析农田生产潜力影响因素空间异质性。

Table 3. The regression coefficient of OLS model

表3. OLS模型回归系数

如表4所示研究时段五个时间节点云南省区县内耕地面积与对应区域的农田生产潜力GWR回归系数为正值的区县均超过123个县，表明区县内随着耕地面积增加农田生产潜力也随之增加；人口与农田生产潜力回归系数为正值的区县均超过117个县，表明区县内农田生产潜力随人口数量而增加；地形起伏度与农田生产潜力回归系数为负值的区县均超过125个县，负值系数表明区县内随着地形起伏度增加农田生产潜力减少。

从1980年至2010年，云南省各区县年平均气温与对应的农田生产潜力回归系数为正值的区县从74个增加至118个区县；但随着县域GDP与农田生产潜力回归系数为负值的区县从57个增加至98个区县，研究时间段内各区县GDP呈上升趋势，伴随着花卉种植、烟草生产、生态旅游等农业人口城镇化进程调整了产业结构，导致区县内农田生产潜力并未随GDP增加而增加。

Table 4. The regression coefficient of GWR model

表4. GWR模型回归系数

5. 结果与分析

本文通过莫兰指数、热点分析、地理加权回归几种空间分析的方法研究了1970、1980、1990、2000、2010五个时间节点农田生产潜力与各影响因素之间的关系。结果表明：

1) 1970年至2010年云南省农田生产潜力总体呈增长趋势，农田生产潜力增量超过5000 t的县域持续增加，且50年间农田生产潜力增加的区县主要分布在云南省中部，呈现接连成片分布格局。

2) 1970年至2010年云南省农田生产潜力重心在研究时间范围内农田生产潜力中心位置从楚雄州禄丰县到玉溪市易门县，云南省农田生产潜力重心移动距离较短速度较慢，空间格局变化较稳定。

3) 云南省农田生产潜力空间分布具有显著的空间集聚特征，云南省西北部横断山区位置受地形影响地形起伏变化大耕地面积少、年均温较低，尽管降水丰富，但农田生产潜力仍为较低水平。云南省中部包括曲靖市下辖区县和昆明市西部区县耕地面积资源较丰富，地形起伏度较小区域农田生产潜力值较高。

4) 耕地面积、人口对农田生产潜力为正向拉动力，地形起伏度对农田生产潜力为负向拉动力。

1970年至2010年云南省农田生产潜力呈增长趋势且空间格局变化稳定，因此，对于农田生产潜力较高的中部地区包括曲靖市以及昆明市西部应当切实保护耕地、推进农业机械化，用科学技术、生产管理方式缩小农田生产潜力与实际粮食生产力之间的差距；对于自然环境条件恶劣的西北地区要注重农业生产与资源开发的协调度，顺应自然条件因地制宜规划农业生产。