1. 前言
如何实现土地资源优化配置已经成为土地科学研究的热点与重点内容 [1] ,其研究的主要目的是通过一定的科学技术手段,在不同尺度上合理安排各类土地利用布局与结构,以提高土地利用效益 [2] 。上世纪90年代以来,国内外众多学者从不同侧面对土地资源优化配置开展了深入研究,如Wang Xinhao等基于GIS技术实现了流域尺度上的土地资源优化配置 [3] ;Fuhu Ren等耦合了GIWIN-LRA模型,在土地适宜性评价的基础上对土地资源优化配置进行研究 [4] ;Jerry M. Faris构建ULAM模型、Chuvieco运用线性规划模型实现了土地利用结构的优化配置 [5] [6] ;周宗丽等人基于线性规划模型实现了农业用地的优化配置 [7] ;柯新利等人建立了基于粮食安全模型的多种土地资源优化配置方案 [8] ;魏伟等人以干旱区内陆河流域为例,建立了基于CLUE-S模型的土地利用优化配置分区方法 [9] ;黎夏、刘小平、杨青生等人将CA模型与神经网络、多智能体系统、蚁群算法等相结合模拟预测城市用地扩展过程,以实现土地资源的时空配置 [10] - [15] 。从相关研究看,多数学者将人工智能、数理统计与GIS相结合以构建土地资源时空配置优化决策模型,取得了一定的研究成果,但较少将土地利用优化配置与区域土地利用总体规划相结合,难以获得土地利用规划约束指标下的最优布局方案,从而确保土地资源可持续利用。本文以经济快速发展地区广州市黄埔区为例,探讨不同限制条件下土地资源时空配置的技术与方法,为区域土地资源合理利用提供政策建议。
2. 思路与方法
2.1. 典型样区选择——广东省广州市黄埔区为例
广州市新黄埔区地处广州市中部,位于珠三角洲中心地区,介于东经113˚23'29~113˚36'2,北纬23˚1'57~23˚24'57,东西宽约17.4 km,南北长13.5 km,其中陆地面积占全区总面积的77%,水域面积约为23%。它是新一轮对外开放“一带一路”重点发展区域,“三大战略枢纽”关键枢纽,“一江两岸三带”宜居城市建设重叠区域,是粤中向粤东、粤北辐射的过渡地域。
2.2. 研究思路
以土地利用总体规划主要参数为参考标准,以经济、社会和生态效益最大化为目标,以参与土地利用优化配置的主要主体为研究对象,定义多智能体土地利用优化配置过程中行为主体决策,最后将多智能体系统(Multi-Agent System,简称MAS)与微粒群算法(Particle Swarm Optimization,简称PSO)相结合,以帕累托验证寻求模型最优解,从而实现土地利用时间与空间上的优化配置。
2.3. 土地利用优化模型
2.3.1. 变量设定
模型中,首先将黄埔区的2015年土地利用现状图合并为k种土地利用类型,并分割为
个10 × 10 m大小的单元,各土地利用类型面积用
表示,k为土地利用类型代码。
2.3.2. 目标函数
1) 经济效益。主要通过实现各类土地利用在经济方面以低成本高收益为主的优化配置。因此,选取土地利用经济总产出最大化区域和最小化区域所转化的成本作为经济效益优化目标结果。
(1)
式中:
代表区域经济总产出,
为各种土地利用类型经济产出系数。
(2)
式中,
表示土地利用转换总成本;
为土地利用单元
位置上土地利用类型k的成本输出;
为二值变量,即空间单元
的土地利用类型为k时,其值为1或0。
2) 社会效益。是指为社会发展所提供的基本社会保障服务,本文以社会基本保障功能作为社会效益目标优化的结果。
(3)
式中:
为区域社会功能保障的总值;
代表土地利用类型k单位面积提供社会功能保障的价值系数。
3) 生态效益。指土地利用优化配置对生态环境友好度和生态服务功能的影响。因而,该目标以对生态环境的兼容性和对生态系统的服务功能的总价值为优化对象。
(4)
式中,
为土地利用对生态系统功能总值;
为某一土地利用类型k在单位面积上的生态服务功能价值系数值。
(5)
式中,
为区域生态环境与土地利用兼容性的总和;
为二值变量,即空间单元与
之间存在邻接,其值为1或0。
为土地利用类型k空间单元与
空间单元
之间的生态环境兼容性。
2.3.3. 约束条件
约束条件为实现设定目标函数所需的限制因素,本文主要从以下方面进行设定:
1) 人口总量是土地承载总人口不应该超过国家规定预期人口数量:
(6)
式中,
为农用地的人口密度;
为城市建设用地人口密度;
表示农用地面积;P为未来规划区域人口总数。
2) 土地面积为各种土地利用类型面积之和为研究区域总面积:
(7)
式中,S为研究区域总面积。
3) 劳动力资源为不同土地利用类型从事相关职业的劳动力数量:
(8)
(9)
式中,
表示第k种农业用地单位面积所承载劳动力数量;
表示第k种非农业用地单位面积所承载劳动力数量;
表示从事农业的劳动力数量;
从事非农业的劳动力数量。
4) 粮食安全
(10)
式中,
表示可以确保期限内粮食安全的耕地需求面积;
表示规定人均食用粮食的标准;
表示研究区域粮食自给率;
表示研究区域粮食单位面积产量;
表示研究区域粮食作物占全部种植农作物比重;
表示研究区域复种指数。
5) 生态环境
(11)
式中,
表示土地利用规划期林地面积;c表示林地面积转换的系数;r表示规划期内确保的森林覆盖面积。
6) 市场经济
(12)
(13)
式中,
表示城镇建设用地、工矿用地、交通设施用地等在土地利用规划期末的面积;
表示上述类型土地在土地利用规划基期的面积;
表示研究区域土地利用规划末期未利用土地面积;
表示研究区域土地利用规划基期未利用土地面积。
7) 土地利用公平性
a) 土地利用规划期末的各类土地生态能力总和不得低于基期:
(14)
式中,
表示农用地类型k在土地利用规划基期的生态生产力;
表示农用地类型k在土地利用规划期内最大能承受的生态产出损失;
表示农用地类型k的生态产出系数。
b) 耕地非农化时序约束条件:
(15)
式中,
表示在土地利用规划期内t时间点耕地非农化面积;Q表示在土地利用规划期限内允许耕地非农化面积。
8) 空间布局
a) 空间单元所含至少一种土地类型:
(16)
b) 新增建设用地在空间扩展时的紧凑性:
(17)
式中,
表示为二值变量,即空间单元
为建设用地时,其值为1或0;
空间单元
存在邻接的集合,确保新增建设用地空间单元
相邻的单元建设用地空间单元不少于2,从而避免建设用地的无限扩展 [16] 。
2.4. 多智能体及决策行为
土地资源优化配置的核心主体具有明显的层次性,第一层次是对土地配置起宏观主导作用的上级政府,第二层次是对土地配置进行合理引导的地方政府,第三层次为各土地利用实际主体,三者之间的决策行为具有承上启下的作用。
2.4.1. 上级政府Agent (Agent1)决策
Agent1依据国家政策和区域总体发展目标,凝练区域土地利用战略,并在土地利用规划中分解各项规划指标的限制性因子。本文采用层次分析法对各因子依次赋予权重,并进行多因素综合分析模型,确定112个行政村(社区)占黄埔区的比例,并依据(公式19)确定各行政村的指标。
(18)
式中,
表示第j个行政村(社区)指标分值;
和
表示每个因子的权重;
表示每个因子的标准化数值。
(19)
式中,
表示第j个行政村(社区)指标;
表示行政村(社区)指标分值;
表示全区的总指标。
2.4.2. 地方政府Agent (Agent2)决策
借鉴美国学者约翰逊提出的“生态位”的思想,将Agent2所组织的黄埔区不同土地利用类型对各种资源的需求构成n维资源空间,公式如下:
(20)
式中,ES表示不同类型土地利用在多维资源结构中的生态适宜水平;
表示资源中某一种资源生态适宜水平;n表示资源空间的维数。
通过地方政府Agent2的发展目标可以计算不同维度的资源需求类型,而后根据不同的发展目标,确定生态因子适宜水平
,最后依照生态适宜水平单元评价结果为基础,以适宜水平最高值所对应土地类型为标准,结合黄埔区实际划分的土地利用功能区,并划定出不同功能区中相应的土地主导类型和空间管制规则。
2.4.3. 利用主体Agent(Agent3)决策
该决策从土地利用主体之间的合作与竞争为基础,采用微粒群优化算法(PSO)设定使用主体Agentu的领域,但不能与其他领域中的Agentu产生合作与竞争的操作,需要与适应能力最强的Agentu在群体中实现信息共享,因而需设计Agentu三种算子实现空间决策,其中包括合作与竞争、自学性和适应性算子,而所有的Agentu均分布在
的土地利用网格上。
1) Agentu合作与竞争算子
为确保Agentu在模型中的交互性和自治性以及空间检索能力,实现Agentu顺畅的信息交流。本文结合PSO算法的全局空间检索能力,设计Agentu合作与竞争算子,以完成Agentu在空间上的决策。
区域土地利用优化配置模型将每个解都看成为一个n维空间中粒子。
的空间解表示的土地空间优化配置的位置向量表示为
,其中,
代表
对目标的决策偏好,因Agentu具有适应性,因而随着模型的变化,产生迭代变化,从而推进模型趋向于最优化。
粒子的飞行速度
。在众多适应度中
为Agentu所拥有的最大适应度,记作。
表示领域间竞争概率,
表示领域间合作概率,
表示均匀分布的随机常数。当
,则模型执行竞争算子;当
,则模型执行合作算子。
a) 竞争算子。当
满足公式21时,其空间解的位置不变,否则
的空间解的位置将按照公式22进行修正,且当
所在网格的土地类型将调整的与
保持一致。
(21)
(22)
式中,
表示为
之间的随机常数量。当
,则
;当
,则。
表示优化解的可行空间下限,
表示可行解上限。
b) 合作算子。在
和
选择随机交叉点,并将交叉点后半部分进行交换,同时在交叉点处进行交叉,从而获得新的Agent。假设
和
在k点进行交叉,而后所产生的两个Agent分支则为:
(23)
(24)
式中,
,
,
和
为
的取值范围,β为
之间的随机常数。
在迭代过程中,但用公式25和公式26对Agentu的空间解和速度的位置进行不断更新:
(25)
(26)
式中,d为模型迭代次数,w表示为惯性常数,
和
表示为学习因子,用于调节个体极值点和值点的飞行步长;
为(0,1)之间的随机常数。
表示
第d次迭代时的n维分量;
表示
第d次迭代速度n维分量;
表示
自身寻求最优Agent的n维分量;
表示模型中最优Agent的n维分量,而后通过全局最有交换,可加快环境中的信息快速传递,提高POS的收敛性。
2) Agentu自学性算子
Agentu自学性算子主要用于模型中自适应性最强的Agentu,学习检索的半径为小于5-10个空间单元(T),并以此为中心,检索产生
的单元合集,进而对Agentu执行竞争和合作算子,经过不断迭代,该Agentu的自学性将会进一步提高。
3) Agentu适应性算子
区域土地利用优化配置模型中,通过适应性调整Agentu的适应能力,而Agentu又可用适应度的数值表示大小,用
,表示自适应度函数,n为目标个数;对于每个不同的目标x,Agentu
会通过优劣生成可行解排序Y,Agentu初始自适应度公式如下:
(27)
(28)
式中,
表示
所有Agentu对目标
的优劣排序序号;
表示
对目标x的适应水平;
表示
目标x的决策偏好;
表示
对所有目标初始适应水平总和。
在优化配置过程中,
的适应水平会受时间t的空间单元对K土地类型的适宜性
的影响,周围空间单元对
在时间t转化为K土地类型的影响
,Agentu在时间对
实行土地引导
以及随机干扰项V的互相影响,综合考虑上述影响后的
自适应性函数如下:
(29)
式中,
为多因素影响Agentu
的自适应性总和。
4) Pareto寻优策略
在Agentu的空间决策中,一方面为了提高模型在时序上的动态调节性以确保土地利用状态达到最大功能,另一方面则为可以提供多种优化配置方案,满足不同需求下的土地利用状态,达到因地制宜的效果。本文采用动态Pareto寻求最优解以获得最优策略,同时排除非Pareto最优解,从而保证Pareto库中均为最优解。
3. 结果与分析
3.1. 模型简化
3.1.1. 上级政府Agent1的决策行为选择
根据广州市在规划期内下达的土地利用约束性与预期性指标,确定黄埔区土地利用优化配置的主要控制指标(表1),并将其落实到具体的行政村。
3.1.2. 地方政府Agent2的决策行为选择
根据Agent1制定的黄埔区土地利用规划战略目标,从黄埔区实际情况出发,Agent2当下需要组织的黄埔区土地利用功能主要包含四种类型:城镇发展结合部、新兴产业园区拓展部、高产农业发展区、原始生态保护区。针对不同类型的发展战略目标,城镇发展的生态适宜水平、工业类的生态适宜水平、农业发展生态适宜水平以及生态保护的适宜水平评价因子体系如表2。
以10 m * 10 m的栅格作为黄埔区生态位适宜水平评价单元,先将数据进行统一标准化,然后针对不同土地利用功能所对应的土地需求类型的方法可以得到多个生态位适宜水平值
,而后将不同发展战略下的
代入公式20,从而获得不同土地利用功能土地在不同评价单元生态适宜水平ES。初步规划黄埔区土地利用功能区,再根据黄埔区经济社会发展战略目标、土地资源现状、土地利用政策等以及Agent2下达给各行政村(社区)及各部门的指标进行局部调控,最终将黄埔区分为4种类型的土地利用功能区:城镇发展区(引导建设区)、新兴产业园区(引导建设区)、高产农业发展区(限制建设区)、原始生态保护区(禁止建设区) (图1),分别对应的土地利用类型位:城镇建设用地、独立工矿用地、农业用地(耕地和园地)、生态用地(林地、水域和未利用地)。
3.1.3. 利用主体Agent3的决策行为选择
根据2007年的土地利用分类标准,将Agent3分为耕地、园地、林地、牧草地、城镇建设用地、农村居民点用地、独立工矿用地、交通运输用地、水利设施用地等9大类Agent,其中未利用土地没有相对应的土地利用行为,因而暂不考虑。每一个网格上只能分配唯一土地利用类型,可容纳的Agent3也只能为一个。在这里Agent3仅能反映的是简单的比例关系,在本文的实际含义为10 m * 10 m的网格上可容纳的实际使用单位或个人的平均值,将不同类型的Agent3配置在2015年土地利用网格上,随着Agent3之
![](Images/Table_Tmp.jpg)
Table 1. Main indicators of land use optimization allocation in Huangpu (hm2)
表1. 黄埔区土地利用优化配置主要控制指标
![](Images/Table_Tmp.jpg)
Table 2. Factors of niche-fitness evaluation for land use function zoning
表2. 土地利用功能的生态适宜水平因子评价体系
间的合作、竞争和自学性行为的迭代,将会引起Agent3的不断变化。
Agent3通过合作、竞争和自学性行为运算进行的空间决策须按照公式29计算Agent3的自适应水平。确定Agentu模型设定目标的观测值以及、
、
、
。
以
空间单元拟被转化的k土地类型与
所在功能区兼容性程度越强,
数值越大,采用Delphi Method求得各种土地利用类型的与功能区的兼容性如表3。Agent3对目标(1)~(5)的观测值所在的经济产出、转换成本、生态功能服务功能价值、生态环境兼容性和社会保障功能总和得到。
![](//html.hanspub.org/file/12-1120696x170_hanspub.png)
Figure 1. Land use function zones of Huangpu
图1. 黄埔区土地利用功能分区图
![](Images/Table_Tmp.jpg)
Table 3. Compatibility between land use types and function zones
表3. 土地利用类型与功能区的兼容性
3.2. 结果与分析
经过最终求解,土地利用优化配置模型共获得6个帕累托最优解。因而得到黄埔区土地利用优化配置方案共计6个,如图2和表4所示。
为更好的辅助决策规划,根据所建立模型的目标函数设置,严格遵循可操作性、科学性和整体性的原则,选取经济产出、转换成本、生态功能服务功能价值、生态环境兼容性、社会保障功能、人均耕地面积、人均建设用地面积、景观聚集水平以及景观平均斑块等9个评价指标体系,采用灰色关联投影法
![](//html.hanspub.org/file/12-1120696x171_hanspub.png)
Figure 2. Comparison of spatial patterns of alternative sustainable land use optimization
图2. 土地利用优化配置方案空间对比
![](Images/Table_Tmp.jpg)
Table 4. Comparison of structures of alternative land use optimization allocation (hm2)
表4. 土地利用优化配置方案数量结构对比(hm2)
![](//html.hanspub.org/file/12-1120696x172_hanspub.png)
Figure 3. Preferred result of sustainable land use optimization allocation in Huangpu
图3. 黄埔区土地利用优化配置方案择优结果
优化方案选择,投影值越大,说明模型所得方案与理想方案越接近,在社会、经济和生态效益上的契合度以及产值上就越高。其基本思路是:结合灰色关联投影的方法和原理,构造对比矩阵并对所选指标进行统一标准化,然后采用SPSS进行关联度计算与判别矩阵的构建,对指标再次进行正交变换最终构造出理想方案,最后对图2和表4中的方案进行投影确定最优方案,结果如图3所示,从中可以看出,最终优方案排序为:方案2 > 方案4 > 方案6 > 方案3 > 方案5 > 方案1。方案1和方案5需要配置的建设用地较少,不能满足黄埔区未来经济和社会发展过程中所需用地,虽然耕地和生态用地较多,满足了区域的生态效益,但忽略了区域经济效益和社会效益;方案3绝大多数未利用地被全部转换为耕地且部分低效工业用地通过复垦转变为农用地,城镇建设用地和独立工况以及交通设施用地配比较小,该方案能够有效保障黄埔区粮食安全问题,短期内提高了社会效益,但从长远区域经济发展角度来看,耕地过多,对支持经济快速发展的力度较小;方案6中,建设用地配置相比较其他方案来说最多,但耕地和林地面积相对较少,能够有效满足社会、经济效益,但忽略了区域生态效益;方案2和方案4为最接近模型所构建的黄埔区土地利用规划功能分区。方案4中的林地配置较为充沛,能够满足区域生态效益的发展,但是耕地保有量未能达到规划中土地利用规划控制标准;方案2中农用地和建设用地数量配置最为合理,土地利用规划的控制指标达到规划要求,且空间分布紧凑,是众多方案中社会、经济和生态效益最大的方案。
方案2、4、6的灰色关联的投影均在0.85以上,其中方案2的投影值达到0.951,非常接近理想值1,这说明模型的建立能够实现社会发展可接受、经济发展可满足、环境发展可承载的土地利用格局。为了进一步验证各方案的可行性和科学性,在公式1模型目标和约束条件下结合元胞自动机模型和线性规划对3个方案数量和空间进行第二阶段的优化配置,采用上述灰色关联投影法对上述3个较优方案进行最终评价,方案2、4、6的灰色关联投影值依次为0.697、0.654、0.631,从而再次印证了方案2为最佳方案。
4. 结论与讨论
1) 本文以我国经济快速发展地区广州市黄埔区为例,以土地利用经济、社会和生态效益最大化为目标,定义了三个不同类型的智能体,通过嵌入帕累托最优算法与POS优化算法相结合,构建了不同的土地资源优化配置的决策方案,结果表明该研究思路是基本可行的,为区域土地资源优化配置提供了一种新思路,研究结果可为区域土地资源合理利用提供参考建议。
2) 多智能体系统(MAS)与微粒群算法(PSO)“自下而上”的集群模型在其它领域的应用已日趋成熟,该模型的关键是如何准确地获得各Agents的决策行为及规则,并通过定量化方法进行表达,如果准确确定公式中的参数非常困难,需要开展更多的实验论证。
3) 土地资源的优化配置是一项系统工程,其配置的合理与否关键是要实现社会、经济和生态效益等综合效益的最大化,本文的研究思路中,对社会与经济因素考虑的较多,对环境因素,如地形、河流等考虑较少,对模型运行的多种方案的综合比较可能还不够合理,今后还需要对研究思路与方法作进一步的改进以提高其运行精度。
基金项目
国家自然科学基金项目(批准号:41001054);教育部人文社会科学(批准号:13YJCZH016、13YJA840009、14YJA630053);国土资源部城市土地资源监测与仿真重点实验室开放基金资助课题(批准号:KF-2016-02-003)项目资助。
NOTES
*通讯作者。