1. 引言
耕地非农化是指耕地转为非农建设用地过程,其转化过程涉及政府、农户、土地开发企业等复杂主体 [1] 。目前耕地非农化成为全球十分关注的一个重大问题,已经受到国内外众多学者的关注,研究领域主要集中在流转成因(驱动力) [2] - [7] 、影响结果(经济增长贡献、粮食安全) [8] [9] [10] [11] [12] 、调控措施 [13] [14] [15] 等方面。目前对于中国耕地非农化的特征、规律以及驱动力等机制的研究已经成熟,然而对耕地如何合理非农化流转的关键科学问题还没搞清楚,对耕地非农化合理性的相关研究较少,对耕地非农化流转决策的关注不足。其主要原因是耕地非农化流转的科学决策过程属于多主体、多目标决策问题,不同主体之间的竞争性和复杂性使得对其优化决策变得非常困难。
耕地非农化流转决策研究目前仍处于起步阶段,定性研究多,定量研究少,静态方法多,空间的动态研究方法少。由于影响耕地非农化流转的影响因素和决策主体不断地在空间和时间上进行协调与变化,这就要求耕地非农化流转决策模型能够更加灵活实用。另外,在空间决策与优化问题中,空间信息的处理与空间搜索能力是关键,常规的GIS空间分析功能无法满足其需求,一般的穷尽搜索方法由于耗时长,搜索速度慢,难以胜任复杂空间搜索问题的求解,这些问题都极大的限制了耕地非农化流转决策优化技术的应用与推广。近年来,基于多智能体系统的土地利用优化配置与空间模拟已经得到广泛的研究和应用 [16] [17] ,这为耕地非农化决策研究提供了研究思路与方法。
本文以广东省广州市黄埔区为例,通过探讨智能体与外部环境间的作用和反馈关系、智能体之间互动关系等问题,在多智能体决策框架的基础上,建立耕地非农化流转的时空决策优化模型,以提升耕地非农化流转决策的科学性。
2. 研究区域概况
本文以我国经济快速发展地区广州市黄埔区为例(图1),黄埔区由2014年的旧黄埔区和萝岗区组成,其位于东经113˚23'29''~113。36'2'',北纬23˚1'57''~23˚24'57'',地处广州市中部,珠江口内,地形南北狭长,地势北向南逐渐降低。包括1个国家级经济技术开发区(广州科学城)和1个省级经济技术开发区(中新知识城),土地总面积约48417 hm2,下辖14个街道、1个镇,其中社区居委会91个,街道15个,行政村48个,村民委员会542个。
3. 数据来源与模型构建
3.1. 数据来源
研究基础数据包括:2000年、2015年1:10000的土地利用变化调查数据,将每个图班作为数据转换的基本单元,2015年10 m分辨率SPOT遥感影像数据、统计年鉴数据、土壤调查数据等,所需基础地理信息数据如交通、基础设施数据等来自百度地图。
3.2. 模型构建
在进行耕地非农化流转空间决策过程中,不同的决策主体有其各自的活动目标,如对耕地承包经营者而言,他们考虑的是如何从农业生产中获取最大的经济收入,土地开发企业考虑的是如何以最小的投入获取最大的经济报酬,而政府追求的是经济、社会和生态效益三者的最大化。决策主体之间不同的目标追求决定了他们有不同的行为特征,这些行为可以是不同个体之间的直接关系,也可以是许多个体在共同环境下的间接关系。如何将不同利益相关者的决策行为特征以定量化的形式确定下来,并形成不同的决策规则,是复杂空间决策模型建立之前要解决的关键问题。
3.2.1. 模型智能体及其决策规则
建立基于多智能系统的耕地非农化流转决策优化模型,关键在于如何对各主体agent进行合理的定量描述,明确说明各类型Agent是如何对外部信息进行有效探测,并对相对应的外部信息做出反应并最终做出的行为选择。在本文研究中,主要考虑三种类型Agent及其行为决策,分别为:政府Agent、土地开发商Agent、农户Agent,其他Agent暂不在考虑范围内。
1) 政府Agent及其决策规则
在耕地非农化流转决策优化模型中,政府作为一类比较特殊的Agent,它与其他Agent不同的是它不具备空间属性,而是通过在制定土地利用和城市规划,从而在时空上对耕地资源进行合理配置。
政府的规划与调控对耕地非农化过程是一个不可忽略的因素 [18] 。当土地开发企业向政府申请耕地流转时,政府会根据该地块目前的耕地利用状况和未来规划情况进行对比,给出不同的接受概率。在此分两种情况,第一种情况是该地块符合规划为建设用地时,说明政府接受的可能性较大,在此定为1;第二种情况是该地块不符合相关规划要求时,一个地块被申请的次数越多,它被政府接受的概率就会增加,这充分体现政府在宏观规划的同时也要考虑公众及开发商的要求。借鉴微粒群优化算法与多智能体的优势,结合GIS技术,从满足耕地动态适宜性和周围土地利用约束对耕地非农化影响的角度,计算耕地单元用地非农化概率 [19] [20] 。其公式表示为:
(1)
其中:
为该地块
被政府接受的概率;
为政府原始接受的概率;g是该地块被申请的次数;
是每申请一次,政府所增加的接受概率;
是以
为中心邻域窗口内已经被政府接受的网络数,
是邻域窗口内每增加一个被政府接受的网格,政府所增加的接受概率。
2) 土地开发商Agent及其决策规则
土地开发商Agent是耕地非农化的主要助力之一,选择土地开发成本最低为出发点,主要考虑因素包括绿地面积(
)、交通通达度(
)、商服水平(
)、教育水平(
)和医疗水平(
)等方面,土地开发商Agent进行选址的公式为:
(2)
式中,
为最低成本开发地块,a、b、c、d、e为上述影响因素的权重,
。其中:
(3)
式中,
表示绿地
状态值,
表示邻域内土地利用类型绿地网格数量,
表示领域内建设用地网格数量。
(4)
式中,
表示交通运输网
状态值,
、
、
表示到地铁、主干道以及普通公路的平均距离;
、
、
表示因子权重系数,
;
、
、
表示空间加强系数;
、
、
表示空间减弱系数。
(5)
式中,
表示商服
状态值,
、
、
表示到超市、百货商店和集市的平均距离。
(6)
式中,
表示教育水平
状态值,
、
、
表示到大学、高中和小学的平均距离。
(7)
式中,
表示医疗水平
状态值,
、
、
表示医院、卫生站和药店的平均距离。
3) 农户Agent及其决策规则
农户Agent及决策行为。影响农户决策的因素有家庭人口、经济状况、耕地质量、耕种习惯以及外部的制度与政策等因素,其中农户从耕地所获取的经济效益是影响农户进行耕地非农流转意愿的主导因素。在农户经济中我们可以构建家庭收入与支出函数,并通过效益与时间模拟的耦合,可以解释农户在何时进行耕地非农化流转时的效益最大,具体的公式如下:
(8)
(9)
式中:
为在时间t地块j在状态r下未来回报的现值;
为在时间t内地块j从状态u转化到状态r的费用;
为耕地未改变用途时的回报;
为贴现率。该公式表明农户在时间T + 1下得到的回报大于在T时段得到的回报,那么农户将会在时间T + 1时存在进行耕地非农化流转的意愿。
4) 耕地非农化流转多智能体系统综合决策
就某个耕地单元而言,耕地非农化是政府、开发商、农户等多个Agent共同决策的结果,是在宏观智能体与微观智能体的共同作用下实现的,采用加权求和模型,该耕地单元在二者影响下的非农化概率可以表示为:
(10)
式中,
为权重参数;P表示耕地单元
在宏观智能体与微观智能体的共同作用下的耕地非农化概率;
表示耕地单元
在宏观智能体作用下的非农化流转概率;
是在微观智能体共同作用下耕地非农化概率,借鉴离散选择模型,得到其求解公式为:
(11)
式中,j为微观智能体编号;L为微观智能体个数,如核心利益相关者为农户、土地开发企业和政府,则L = 3;
为微观智能体j对单元
的耕地非农化流转期望概率;
为微观智能体j所占的权重。在模型中,P值的高低决定了耕地单元非农化的顺序。
依据上述设定Agent决策规则,得到耕地非农化流转多智能体综合决策规则:
(12)
式中,
表示
时点
单元耕地非农化概率;
表示
时点参与运算的耕地非农化适宜性的最大值;
表示离散参数;
表示
单元在政府、开发商和农户多个Agent共同作用下的耕地非农化流转适宜性分值,计算公式如下:
(13)
式中,
表示Agent的编号;
表示Agent的个数;即政府Agent、开发商Agent和农户Agent;
表示t时点x单元期望流转概率;
表示Agent
在的权重;k表示耕地非农化限制因子,取值规则如公式(12):
(14)
式中,0表示如果x单元的耕地类型为水田、旱地、水浇地;1表示如果x单元的耕地非农化为其他类型建设用地等。
3.2.2. 目标函数及约束条件
选择4个耕地非农化流转过程中的优化目标,即最小化地块开发费用、土地利用的适宜性、最大配置环境兼容性、资源节约型和环境友好型 [18] ,模型的目标函数与约束性条件如公式(13)~(16):
1) 最小化地块开发费用(
)。主要通过实现耕地非农化后在经济方面以低成本高收益为主的优化配置。
(13)
式中,I、R和B分别表示耕地
单元转变为工业、居住和商业类
的最小成本。
2) 宜耕后备土地资源的适宜性(
)。是指其他类型的土地转化为耕地的适宜性。
(14)
式中,
表示其他类型土地转变为耕地的经济产出系数,
表示其他类型土地流转为耕地的面积。
3) 环境友好型(
)。指耕地非农化对生态环境友好度和生态服务功能的影响。因而该目标以对生态环境的兼容性的总价值为优化对象。
(15)
式中,
表示单位面积耕地非农化流转提供的生态服务功能价值系数。
4) 资源节约型(
)。
(16)
式中,
为耕地可流转最大面积,
为转换系数,
为区域优化过程中需保障的耕地面积。
3.2.3. 模型算法
为提高耕地非农化流转决策的精度与效率,本文采用多智能体系统(MAS)和微粒群优化算法(PSO)进行运算,选择该模型的主要原因是:微粒群优化算法作为多目标智能优化算法,具有开放性大、效率高、问题优化求解能力强等特点,但缺乏微观主体决策过程,难以反映微观空间演变规律,通过与多智能体系统算法的结合可以有效克服优化算法的缺陷,并获得最优的方案选择和结果。
以2015年土地利用变更数据作为基础数据源,以2020年为目标年。将该区域土地利用现状图,栅格化为10 × 10 m的网格,其空间布局包括耕地、园地、林地、城镇建设用地、交通用地、水域和其他未利用地。耕地流转为其他类型建设用地的环境兼容性为0,保留耕地属性的环境兼容性为1。耕地网格数为325,743、园地网格数为586,392、林地网格数为1,613,117、城镇建设用地网格数为1,644,208、交通用地257,713、水域383,275、其他未利用地31,252。测试要求为根据本文中提出的决策规则及约束条件,对耕地非农化流转配置方案进行空间决策优化,优化后耕地网格数不得低于280,000、城镇建设用地网格数最高达到1,676,136和林地覆盖不得低于1,676,136,依据测试要求,共有82,589个网格需要重新分配。
4. 结果与分析
模型参数取值如下:邻域合作概率和竞争概率均为0.85,最小误差值为0.0000001,学习因子取值1.22,常数因子为0.77,迭代次数为10次,搜索半径为0.2,研究为各个Agent设置相同的优势度,取值为1,对耕地非农化流转过程的多个目标进行优化,计算获得全局最优解作为判断依据。模型运算初期,全局最优解主要集中在四个目标值的较低水平,随着模型迭代次数的增加,非支配解的各个目标逐步提升,最终达到收敛状态。整个运算过程是在黄埔区《基于多智能体的耕地非农化流转决策支持系统》中运行,经最后检验,kappa系数达到0.8035,超过0.75的一致性检验水平,表明此次运行结果良好。最终运算结果如表1所示:
通过2015年土地利用现状图与优化配置结果图进行对比(图2),耕地由原来2256个斑块通过优化配置后减少为1963个斑块,城镇建设用地斑块从3031个斑块减少到2744个斑块,两者的空间紧凑度在优化后获得了不同程度的提高。
将2015年土地利用现状图与优化结果进行叠加对比,进一步分析耕地流向及空间分布格局。表2中显示发生变化的用地主要为耕地、城镇建设用地和交通用地,分别占变化总量的50.54% (47,104)、32.71% (30,485)和13.99% (13,036),林地、水域和其他用地未发生转变。具体来看,耕地流转过程中,约22.60%由于退耕还林政策调整为林地,约58.94%由于城市扩张转变为城镇建设用地,约27.67%由于完善城市交通网而转变为交通用地。同时,为落实耕地占补平衡政策,通过建设用地整理复垦,约占城镇建设用地的35.25%的10,747个单元网格调整为耕地。
![](Images/Table_Tmp.jpg)
Table 1. Operation optimization results of different models of cultivated land conversion
表1. 耕地非农化流转运行优化结果
![](//html.hanspub.org/file/26-2391040x114_hanspub.png)
Figure 2. Actual land of area in 2015 and optimized configuration results in 2020
图2. 2015年现状图与2020年耕地流转结果图
![](Images/Table_Tmp.jpg)
Table 2. Transition matrix of farmland
表2. 耕地转换矩阵
在耕地转换过程中,位于中部和北部坡地或者土地质量较差的耕地被转为其他用途,而土地质量较好的单元则优先被补充为耕地。耕地总面积虽有所减少了,但是耕地质量有明显的提高(见表3)。在资源节约和环境友好目标的指导下,新增耕地主要集中在黄埔区九龙镇北部地区,该地区为广州市主要高标准基本农田建设区,考虑到高标准基本农田的集中连片性,该区域范围内新增耕地占85.63%。经过优化
![](Images/Table_Tmp.jpg)
Table 3. Changes of arable land quality
表3. 耕地质量变化分析
后耕地等级质量上升了两个等别,从i-V等变为i-iii等,且i等的耕地数量占到所有耕地面积的86.01%,大幅度提高了优质耕地数量;而新增城镇建设用地主要集中在黄埔区中部、九龙镇南部地区,主要原因在于中部地区和九龙镇南部地区相对地势平坦,非城镇化区域范围随着政府“北进政策”,城镇建设用地逐渐趋于集中布局,提高了土地资源的集约利用,避免了优质农田的过渡占用。
5. 结论与讨论
1) 耕地资源非农化流转决策优化是以实现社会、经济和生态效益最大化为目标,合理安排耕地非农效率,最终实现生态系统平衡和土地资源可持续利用。本文以区域耕地非农化流转决策为研究目标,在宏观层面上考虑到影响耕地非农化的主要制定主体和执行主体,在微观层面上涉及农户主体行为特征,最终在地理空间上实现耕地资源优化配置,证明了本文的实用性。
2) 将2020年与2015年耕地非农化优化前后空间格局进行对比,发现MASPSO模型运行的结果中建设用地呈现出集中连片的空间布局,零星的土地利用斑块大幅度减少,土地利用的空间集聚度明显提高,土地利用空间格局得以优化,耕地质量得到有效提升,从而进一步证明了MASPSO模型进行耕地非农化流转决策的可行性与优越性。
3) 耕地非农化流转决策优化是一个多目标的复杂决策过程,本文结合了“资源节约”和“环境友好”政策对耕地非农化流转决策进行了常规化指标优化,侧重考虑了政府、开发商和农户Agent的空间决策行为,但与基期有较大的差异,这表明研究区的耕地非农化流转存在明显的不合理性,需要相关部门实际工作中进一步改进。
基金项目
国家自然科学基金项目(批准号:41001054);教育部人文社会科学(批准号:13YJCZH016、13YJA840009、14YJA630053);国土资源部城市土地资源监测与仿真重点实验室开放基金资助课题(批准号:KF-2016-02-003)项目资助。
NOTES
*通讯作者。