论文部分内容阅读
摘 要：采用空间自相关方法，通过网络爬虫技术获取到的网络平台土地流转数据，探讨浙江省各地级市的农用地流转时空模式和地区差异。实验结果表明：随着网络平台建设推广和农用地流转市场的发展，2014年以后浙江省各地级市的流转规模有明显空间聚集特性，表现出以衢州市-金华市-绍兴市为分界带，呈现出西北角为热点区域，东南角为冷点区域的特征。本文的贡献在于：①提供基于互联网数据的农地流转研究新视角，突破以往农户调查方式受限于调查精力、技术水平和体制机制问题无法得到多年份、跨地区数据的限制；②引入空间自相关分析方法，对农用地流转的地区差异和时空模式的定量分析提供参考。
关键词：农用地流转；网络平台；空间自相关；时间序列分析
1 引言
农用地是人类不可替代的自然资源基础，其数量和质量特征决定着一个国家或地区社会经济的可持续发展[1-4]。为盘活农村土地市场，我国于十八届三中全会将发展和引导农用地流转作为农村土地改革的重点之一[5]。在国家政策和农业发展需求双重刺激下，我国农用地流转规模不断加大，流转市场稳步发展。特别是近年来，随着网络集成商等主体纷纷试水土地流转业务[6]，土地流转网络平台的建设打破了过去村集体组织、用地单位、地方政府容易形成一个利益集团的局面[4]，大大缓解了因信息不对称造成的中央政府调控失灵、农民只能被动接受土地流转的情况[7，8]，土地流转网络平台逐渐在农用地流转市场中发挥着重要的角色。
通过文献阅读[1，3，4，7-12]和分析发现，以往农用地流转研究存在两个缺陷：一方面，研究内容主要来源于抽样问卷调查，样本容量的限制研究区域多局限于均质性的小区域，研究多基于现状时段，缺少对大范围区域和长时间序列的农用地流转的研究。另一方面，研究方法多采用统计方法和计量模型，较少考虑到土地流转的区域性，更无法从微观层面考察时空差异的机理和影响因素。
针对以上问题，本研究尝试利用土地流转网络平台累计的数据，结合空间自相关理论，对网络平台的农用地流转的时空差异特征和区域特征展开探究。空间自相关分析是研究人文和自然地理现象的空间差异性主要方法，本文引入该方法，以流转数量表征流转规模，结合时间变化、尺度变化，分析农用地流转的时空差异与特征。
2 数据抓取与预处理
目前土地流轉市场中较大的代理网站包括土流网、聚土网、地合网和土地资源网等，从研究的区域特征、时间跨度、数据质量等角度综合分析后，本研究使用Python的Requests库从土地资源网（http：//www.tdzyw.com）爬取浙江省农用地流转数据。爬取数据总量为1406条，发布时间从2010-07-29至2017-3-30，数据内容包括流转方式、用地类别、流转面积、使用期限、流转价格、发布时间和土地坐落信息。
经过数据清洗（关键字提取、单位换算、地址解析和判读），去除无法解析市场区域的交易数据，剩余1382条规范的农用地流转信息（见表1）。通过关键字信息提取，将流转方式分转让、转包、出租、入股和其他五类，用地类别分为耕地、园地、山地、养殖场等。流转面积统一以亩作为计量单位，使用期限以年作为计量单位。流转价格根据信息源分为面议和参考价格，将参考价格换算为单位价格，即每亩每年流转价格。数据中的市场区域信息定位至县区级行政单位，以2016年浙江省最新行政区划为基准进行统计，包括90个县区级行政单元。
3 网络平台农用地流转规模时空特征
我们使用空间自相关理论对网络平台农用地的流转规模进行时空特征分析。由于舟山市与其他地级行政区相隔离，选用邻接权重矩阵将产生无邻居的情况，因此在进行农用地流转空间自相关分析时采用距离权重矩阵。选取不同距离阈值的空间权重矩阵，对2010-2016年各地级市进行土地流转统计汇总量计算全局空间自相关指数（Moran's I）。比较发现当距离阈值选取176.3公里时，全省农地流转规模的空间聚集特性最显著（P=0.029<0.05，且Moran's I=0.128>0）。选用其作为本研究空间自相关分析的基础，此时各地级市的邻居数量分布较为均匀，舟山市具有嘉兴市、宁波市、绍兴市、台州市四个近邻。
计算每年的土地流转数量的全局Moran's I指数，进行时间序列分析（图2）。2014年以前，全局Moran's I指数有正有负、变化跳跃大，且显著性水平P值较大（>0.05），反映出2010年～2014年浙江省农用地流转不具备明显的空间自相关性，体现了起步阶段的网络流转平台的农用地流转呈现出空间无序性。随着时间推移，2014年后农用地流转呈现出显著的空间自相关性（P<0.05），Moran's I指数为正并在2016年达到了最大值（0.236），呈现出空间聚集特征。Moran's I指数的增强反映出地级市与周围的地级市间的差异在减少，也反映出网络流转平台的农用地流转市场逐渐推广发展的时间特征。
全局自相关能从整体水平上考量浙江省各地级市的空间自相关情况，但无法聚焦于具体分异情况。因此，对各个年份的农地交易数据进行局部空间自相关分析，在0.05的置信水平上通过GeoDa绘制局部自相关指数的LISA聚集图（见图3）。
与全局自相关的分析趋势相对一致，2014年以前全省的农地流转数量的空间结构特征并不明显，未呈现出一定规律：2010年杭州市属于高高类型，台州市属于低低类型。2011年和2012年仅嘉兴市表现出低高类型。2013年，金华市为高高类型，衢州市为低高类型。2014年开始，全省的农地流转交易以衢州市-金华市-绍兴市-宁波市形成分界带，呈现出西北角为热点区域，东南角为冷点区域的空间结构特征。且杭州市为长期热点城市，台州市为长期冷点城市。
进一步，对县区级尺度进行流转规模的局部空间自相关分析（2010～2016县级），可得到更细粒度的冷热点信息：热点地区包括杭州市下属的余杭区和临安市，湖州市下属的长兴县，嘉兴市下属的嘉善縣，宁波市下属的余姚市，以及金华市下属的兰溪市；冷点地区包括台州市下属的临海市、椒江区、黄岩区和玉环县，温州市下属的乐清市、龙湾区、洞头县、瓯海区、文成县、苍南县、瑞安市。相较于地级市尺度的分析，县区级流转规模的分析结果仍然保持了西北角为热点区域、东南角为冷点区域的结构特征，并且可以发现冷点地区连接成簇，集聚性明显高于热点地区。 4 结论与讨论
（1）土地流转网络平台，作为农用地流转的新形式。一方面能够打破信息不对称，帮助建立透明、公开的流转市场，另一方面也为农用地流转研究提供了真实、可靠的流转信息（流转面积、流转方式、流转价格、土地坐落等）。通过网络爬虫技术抓取到的网络平台数据作为传统问卷调查数据的补充和对照，能对农用地的流转空间差异、时间序列分析和影响因素分析等研究提供新研究视角。
（2）农用地流转空间差异的研究方法上，借助空间自相关的探索性空间分析方法能从空间结构和区域差异角度对流转规模定量描述。GIS制图技术与农用地市场发育测度模型的结合，能对农地流转研究中农地市场的发展程度研究提供有力支撑。
（3）通过对流转规模的时空特征进行分析，发现浙江省流转规模在2014年以前不具备空间自相关特征，2014年以后有明显空间聚集特性。在空间结构上，以衢州市-金华市-绍兴市为分界带，呈现出西北角为热点区域，东南角为冷点区域的特征。进一步，在县区级尺度上进行空间自相关研究，发现：热点包括杭州市下属的余杭区、临安市，湖州市下属的长兴县，嘉兴市下属的嘉善县，宁波市下属的余姚市和金华市下属的兰溪市；冷点包括台州市下属的临海市、椒江区、黄岩区和玉环县，温州市下属的乐清市、龙湾区、洞头县、瓯海区、文成县、苍南县、瑞安市等。
本研究从网络平台农用地流转的视角出发，对流转规模的时空特征分析，揭示浙江省农地流转的空间结构特征和时间演变规律，对进一步引导建立更活跃、规范的农用地流转市场有一定的参考意义。
参考文献
[1]何欣，蒋涛，郭良燕，等.中国农地流转市场的发展与农户流转农地行为研究——基于2013～2015年29省的农户调查数据[J].管理世界.2016（06）：79-89.
[2]张乐天，陆洋.乡土秩序与土地流转的非均衡实践——浙江省Z村一带土地流转研究[J].南京社会科学.2011（07）：55-64.
[3]許恒周，金晶.农地流转市场发育对农民养老保障模式选择的影响分析——基于有序Probit模型的估计[J].资源科学.2011（08）：1578-1583.
[4]李庚.农村土地流转的空间差异及成因分析——以陕西关中三县为例[J].人文地理.2012（01）：102-106.
[5]张四梅.经济结构变化对农用地流转的影响研究[J].湖南师范大学社会科学学报.2015（03）：124-130.
[6]张祚本.2015年网络平台土地流转价格监测分析[J].农村经营管理.2016（07）：32-33.
[7]宋伟，倪九派.基于农户调查的农用地流转影响因素研究——以重庆市典型区（县）为例[J].中国人口.资源与环境.2011（S1）：172-175.
[8]邵景安，魏朝富，谢德体.家庭承包制下土地流轉的农户解释：对重庆不同经济类型区七个村的调查分析[J].地理研究.2007（02）：275-286.
[9]杜文星，翟文侠，黄贤金.土地流转的空间分异研究——基于泰州、南京、上海的农户实地调查[Z].中国乌鲁木齐：20059.
[10]陈立双.农村土地流转空间差异分析——基于辽宁省4城市农村土地流转数据的分析[J].农村经济与科技.2014（06）：31-33.
[11]卢新海，望萌.农用地流转的武汉模式研究——基于武汉农村综合产权交易所的启示[J].农林经济管理学报.2014（03）：244-251.
[12]刘成武，黄利民.农地边际化过程中农户土地利用行为变化及其对粮食生产的影响[J].地理研究.2015（12）：2268-2282.
关键词：农用地流转；网络平台；空间自相关；时间序列分析
1 引言
农用地是人类不可替代的自然资源基础，其数量和质量特征决定着一个国家或地区社会经济的可持续发展[1-4]。为盘活农村土地市场，我国于十八届三中全会将发展和引导农用地流转作为农村土地改革的重点之一[5]。在国家政策和农业发展需求双重刺激下，我国农用地流转规模不断加大，流转市场稳步发展。特别是近年来，随着网络集成商等主体纷纷试水土地流转业务[6]，土地流转网络平台的建设打破了过去村集体组织、用地单位、地方政府容易形成一个利益集团的局面[4]，大大缓解了因信息不对称造成的中央政府调控失灵、农民只能被动接受土地流转的情况[7，8]，土地流转网络平台逐渐在农用地流转市场中发挥着重要的角色。
通过文献阅读[1，3，4，7-12]和分析发现，以往农用地流转研究存在两个缺陷：一方面，研究内容主要来源于抽样问卷调查，样本容量的限制研究区域多局限于均质性的小区域，研究多基于现状时段，缺少对大范围区域和长时间序列的农用地流转的研究。另一方面，研究方法多采用统计方法和计量模型，较少考虑到土地流转的区域性，更无法从微观层面考察时空差异的机理和影响因素。
针对以上问题，本研究尝试利用土地流转网络平台累计的数据，结合空间自相关理论，对网络平台的农用地流转的时空差异特征和区域特征展开探究。空间自相关分析是研究人文和自然地理现象的空间差异性主要方法，本文引入该方法，以流转数量表征流转规模，结合时间变化、尺度变化，分析农用地流转的时空差异与特征。
2 数据抓取与预处理
目前土地流轉市场中较大的代理网站包括土流网、聚土网、地合网和土地资源网等，从研究的区域特征、时间跨度、数据质量等角度综合分析后，本研究使用Python的Requests库从土地资源网（http：//www.tdzyw.com）爬取浙江省农用地流转数据。爬取数据总量为1406条，发布时间从2010-07-29至2017-3-30，数据内容包括流转方式、用地类别、流转面积、使用期限、流转价格、发布时间和土地坐落信息。
经过数据清洗（关键字提取、单位换算、地址解析和判读），去除无法解析市场区域的交易数据，剩余1382条规范的农用地流转信息（见表1）。通过关键字信息提取，将流转方式分转让、转包、出租、入股和其他五类，用地类别分为耕地、园地、山地、养殖场等。流转面积统一以亩作为计量单位，使用期限以年作为计量单位。流转价格根据信息源分为面议和参考价格，将参考价格换算为单位价格，即每亩每年流转价格。数据中的市场区域信息定位至县区级行政单位，以2016年浙江省最新行政区划为基准进行统计，包括90个县区级行政单元。
3 网络平台农用地流转规模时空特征
我们使用空间自相关理论对网络平台农用地的流转规模进行时空特征分析。由于舟山市与其他地级行政区相隔离，选用邻接权重矩阵将产生无邻居的情况，因此在进行农用地流转空间自相关分析时采用距离权重矩阵。选取不同距离阈值的空间权重矩阵，对2010-2016年各地级市进行土地流转统计汇总量计算全局空间自相关指数（Moran's I）。比较发现当距离阈值选取176.3公里时，全省农地流转规模的空间聚集特性最显著（P=0.029<0.05，且Moran's I=0.128>0）。选用其作为本研究空间自相关分析的基础，此时各地级市的邻居数量分布较为均匀，舟山市具有嘉兴市、宁波市、绍兴市、台州市四个近邻。
计算每年的土地流转数量的全局Moran's I指数，进行时间序列分析（图2）。2014年以前，全局Moran's I指数有正有负、变化跳跃大，且显著性水平P值较大（>0.05），反映出2010年～2014年浙江省农用地流转不具备明显的空间自相关性，体现了起步阶段的网络流转平台的农用地流转呈现出空间无序性。随着时间推移，2014年后农用地流转呈现出显著的空间自相关性（P<0.05），Moran's I指数为正并在2016年达到了最大值（0.236），呈现出空间聚集特征。Moran's I指数的增强反映出地级市与周围的地级市间的差异在减少，也反映出网络流转平台的农用地流转市场逐渐推广发展的时间特征。
全局自相关能从整体水平上考量浙江省各地级市的空间自相关情况，但无法聚焦于具体分异情况。因此，对各个年份的农地交易数据进行局部空间自相关分析，在0.05的置信水平上通过GeoDa绘制局部自相关指数的LISA聚集图（见图3）。
与全局自相关的分析趋势相对一致，2014年以前全省的农地流转数量的空间结构特征并不明显，未呈现出一定规律：2010年杭州市属于高高类型，台州市属于低低类型。2011年和2012年仅嘉兴市表现出低高类型。2013年，金华市为高高类型，衢州市为低高类型。2014年开始，全省的农地流转交易以衢州市-金华市-绍兴市-宁波市形成分界带，呈现出西北角为热点区域，东南角为冷点区域的空间结构特征。且杭州市为长期热点城市，台州市为长期冷点城市。
进一步，对县区级尺度进行流转规模的局部空间自相关分析（2010～2016县级），可得到更细粒度的冷热点信息：热点地区包括杭州市下属的余杭区和临安市，湖州市下属的长兴县，嘉兴市下属的嘉善縣，宁波市下属的余姚市，以及金华市下属的兰溪市；冷点地区包括台州市下属的临海市、椒江区、黄岩区和玉环县，温州市下属的乐清市、龙湾区、洞头县、瓯海区、文成县、苍南县、瑞安市。相较于地级市尺度的分析，县区级流转规模的分析结果仍然保持了西北角为热点区域、东南角为冷点区域的结构特征，并且可以发现冷点地区连接成簇，集聚性明显高于热点地区。 4 结论与讨论
（1）土地流转网络平台，作为农用地流转的新形式。一方面能够打破信息不对称，帮助建立透明、公开的流转市场，另一方面也为农用地流转研究提供了真实、可靠的流转信息（流转面积、流转方式、流转价格、土地坐落等）。通过网络爬虫技术抓取到的网络平台数据作为传统问卷调查数据的补充和对照，能对农用地的流转空间差异、时间序列分析和影响因素分析等研究提供新研究视角。
（2）农用地流转空间差异的研究方法上，借助空间自相关的探索性空间分析方法能从空间结构和区域差异角度对流转规模定量描述。GIS制图技术与农用地市场发育测度模型的结合，能对农地流转研究中农地市场的发展程度研究提供有力支撑。
（3）通过对流转规模的时空特征进行分析，发现浙江省流转规模在2014年以前不具备空间自相关特征，2014年以后有明显空间聚集特性。在空间结构上，以衢州市-金华市-绍兴市为分界带，呈现出西北角为热点区域，东南角为冷点区域的特征。进一步，在县区级尺度上进行空间自相关研究，发现：热点包括杭州市下属的余杭区、临安市，湖州市下属的长兴县，嘉兴市下属的嘉善县，宁波市下属的余姚市和金华市下属的兰溪市；冷点包括台州市下属的临海市、椒江区、黄岩区和玉环县，温州市下属的乐清市、龙湾区、洞头县、瓯海区、文成县、苍南县、瑞安市等。
本研究从网络平台农用地流转的视角出发，对流转规模的时空特征分析，揭示浙江省农地流转的空间结构特征和时间演变规律，对进一步引导建立更活跃、规范的农用地流转市场有一定的参考意义。
参考文献
[1]何欣，蒋涛，郭良燕，等.中国农地流转市场的发展与农户流转农地行为研究——基于2013～2015年29省的农户调查数据[J].管理世界.2016（06）：79-89.
[2]张乐天，陆洋.乡土秩序与土地流转的非均衡实践——浙江省Z村一带土地流转研究[J].南京社会科学.2011（07）：55-64.
[3]許恒周，金晶.农地流转市场发育对农民养老保障模式选择的影响分析——基于有序Probit模型的估计[J].资源科学.2011（08）：1578-1583.
[4]李庚.农村土地流转的空间差异及成因分析——以陕西关中三县为例[J].人文地理.2012（01）：102-106.
[5]张四梅.经济结构变化对农用地流转的影响研究[J].湖南师范大学社会科学学报.2015（03）：124-130.
[6]张祚本.2015年网络平台土地流转价格监测分析[J].农村经营管理.2016（07）：32-33.
[7]宋伟，倪九派.基于农户调查的农用地流转影响因素研究——以重庆市典型区（县）为例[J].中国人口.资源与环境.2011（S1）：172-175.
[8]邵景安，魏朝富，谢德体.家庭承包制下土地流轉的农户解释：对重庆不同经济类型区七个村的调查分析[J].地理研究.2007（02）：275-286.
[9]杜文星，翟文侠，黄贤金.土地流转的空间分异研究——基于泰州、南京、上海的农户实地调查[Z].中国乌鲁木齐：20059.
[10]陈立双.农村土地流转空间差异分析——基于辽宁省4城市农村土地流转数据的分析[J].农村经济与科技.2014（06）：31-33.
[11]卢新海，望萌.农用地流转的武汉模式研究——基于武汉农村综合产权交易所的启示[J].农林经济管理学报.2014（03）：244-251.
[12]刘成武，黄利民.农地边际化过程中农户土地利用行为变化及其对粮食生产的影响[J].地理研究.2015（12）：2268-2282.