城市是一个具有开放性、动态性、自组织性的时空动态复杂系统, 城市空间扩展是在特定的地理环境和一定的社会经济发展阶段中, 不同时空尺度下的自然、人文、社会、经济、文化等多种因素的综合反映。随着城市空间的不断扩张, 对自然资源的侵占力度也不断加强, 由此带来一些区域生态环境问题, 如热岛效应、植被退化、生态系统服务功能下降、生物多样性降低等^[1-4]。因此, 对城市时空演化特征进行定量化描述和预测可为区域的可持续发展和土地利用管理决策提供科学依据, 城市空间动态模拟是城市生态环境研究中的重要方向^[5]。

近年来, 许多学者从不同角度建立了数学模型和空间模型开展城市空间扩展和动态模拟研究, 如Clarke等利用城市发展的历史数据, 根据城市的交通、地形条件等设定适当的参数, 建立SLEUTH模型对旧金山和华盛顿都市区进行了成功模拟和长期预测^[6-7]；黎夏和詹云军等将神经网络与元胞自动机(cellular automata model, CA)模型相结合, 挖掘城市扩展土地利用演变的规律^[8\9], 但该方法属于黑箱结构, 对模型参数的物理意义难以解释；刘小平和单玉红等利用智能体模型对城市空间扩张和土地利用格局演变进行模拟^[10-11], 智能体模型能模拟个体或群体的行为及决策模式对土地利用产生的影响, 但个体决策行为具有极大的复杂性, 增加了模型的不确定性, 为模型的验证和移植带来了困难^[12]。其中, SLEUTH模型是一种基于元胞自动机(CA)的城市扩展与土地利用变化模型, 其设计思路是基于交通、地形和城市化的约束条件计算每个元胞单元发展的可能性, 把城市化的元胞作为种子点, 通过其扩散带动整个区域的发展, 其中离种子点越近的元胞越容易被城市化^[6]。与其他模型相比, SLEUTH模型具有较强的普适性和可移植性, 最大的特点是根据模拟结果与历史训练数据的拟合度, 通过蒙特卡洛迭代自动寻找城市增长误差最小的组合方案, 从而得到城市扩展过程最佳的参数组合^[13], 解决了传统CA模型中转换规则不易确定的问题。近年来, SLEUTH模型被成功的应用于国内的城市空间研究, 如北京、广州、长沙等城市^[14-16]。

武汉是我国华中地区最大城市, 也是长江经济带国家战略中长江中游最重要的城市。近年来, 武汉市在“武汉8+1城市圈”的发展规划以及旧城改造的政策背景下, 城市空间向周围的快速扩展, 土地利用结构演变剧烈, 城市化进程日益加快, 由此也在一定程度上造成了武汉市水域湿地面积缩减、水质污染、耕地减少等生态问题, 城市发展、空间扩张与生态环境保护之间的矛盾成为区域可持续发展面临的重要难题。本文以武汉市武昌地区为例, 在已有的城市扩展机制和SLEUTH模型应用研究的基础上, 对武汉市城市空间进行动态模拟, 预测不同情景下城市空间扩张格局, 进而探讨城市土地扩张的生态环境影响及可能采取的对策及效果, 以期为武汉城市空间发展和可持续规划与管理提供参考。

1 研究区域和研究方法 1.1 研究区域

武汉市位于长江和汉水交汇处, 地处29°58′N到31°22′N, 113°41′E到115°05′E之间, 是华中地区的特大城市。本文以武汉市的武昌区、洪山区、青山区和江夏区的北部为研究区, 该区域大部分为平原地区, 湖泊水系众多, 空间格局上具有从城市核心区到郊区的过渡, 产业上由青山区的老工业区和洪山、江夏区的经济开发区并存, 研究区示意图如图 1所示。

1.2 方法与数据 1.2.1 SLEUTH模型

SLEUTH元胞自动机模型可以模拟和预测城市增长和土地利用演化, 模型名称来自6个输入图层的首字母缩写, 分别是坡度图层(Slope)、土地利用图层(Landuse)、排除图层(Exclusion)、城市化图层(Urban Extent)、道路图层(Transportation)和山体阴影图层(Hillshade)^{[6, 17]}。模型转换规则受扩散系数、繁衍系数、蔓延系数、坡度阻碍系数和道路引力系数5个参数控制, 它们决定了城市增长的4种类型, 即自发增长、新扩展中心增长、边缘增长和道路影响增长。本文的城市空间动态模拟具体涉及到居住用地的增长和其他用地类型之间的转换, 因此初始种子元胞是已经城市化的居住用地元胞, 自发增长指在满足坡度条件的可城市化的区域发展新的居住用地元胞的可能性；新的扩展中心增长表示自发增长产生的新的居住用地元胞成为新的扩展中心的可能性；边界增长是原有的扩展中心和新产生的扩展中心周围继续演变为居住用地的可能性；道路影响型增长表示城市交通网络对城市扩张的引力作用。

运用SLEUTH模型进行居住用地增长模拟前需进行校正处理, 即根据历史数据来进行样本训练, 这是一个蒙塔卡罗迭代处理过程^[6]。校正过程可分3个阶段, 即粗校正、精校正和终校正, 根据模型5个参数的初始范围, 通过反复计算不同的参数组合所产生的模拟结果与实际情况的拟合度, 逐步缩小参数范围, 最终确定最优的参数组合^[7], 代表了居住用地增长的历史轨迹方向, 并以此参数组合对城市未来空间扩展进行模拟和预测。

1.2.2 数据来源与预处理

本研究所用基础数据来源包括：2005年、2007年、2009年、2011年4期季相一致、高度30m的武汉市Google Earth截图, 从国际科学数据服务平台下载分辨率为30m的DEM数据。针对SLEUTH模型的运行要求, 对上述数据做了相应的校正、裁切处理, 然后结合遥感解译和目视解译, 在ARCGIS中制作了SLEUTH模型模型所需的居住用地范围图层、土地利用图层、交通图层和排除图层, 利用DEM数据制备了研究区的坡度图层及山体影图层, 输入数据图层如图 2所示, 由于土地利用图层在模拟结果验证中会展示, 在此部分不再列举。最后所有输入图层均按照统一的横轴墨卡托(UTM_WGS1984_49N)投影设定, 全部转换成无符8位灰度的GIF格式的栅格图像。

本文的城市空间动态演变主要关注居住用地与其他土地利用之间的转变, 因此对土地利用类型进行重新分类(表 1), 居住用地图层指各类住宅用地, 具体包括居住小区和城中村居住用地；商业服务业用地、公共基础设施用地和水域作为排除图层参与模型运行, 限制其向居住用地演化；此外, 考虑到除了交通网络良好的可达性对居住用地建设的导向作用外, 还将生态环境好的湖泊、教育用地和设施发达的商业服务业对居住用的吸引作用集成到模型中, 运用ARCGIS软件将这类吸引居住用地开发的区域从面转换成线, 成为“引力道路”纳入交通图层, 不同等级的道路设置不同的相对权重值来表示不同的可达性和吸引力, 一起参与模型的校准、模拟和预测过程。

1.2.3 情景设计

本文基于研究区城市空间变化的历史趋势和城市生态保护需求, 设置了现状趋势、基本保护、严格保护等3种发展情景。利用SLEUTH模型的参数获取结果, 通过修改排除图层和SLEUTH模型的输入参数, 以2011年居住用地范围为“种子点”, 预测其2012—2025年的土地利用变化, 实现3种情景下城市土地利用变化的模拟预测。(1) 现状趋势情景：将公共基础设施用地作为排除图层, 排除概率设为100%, 对其他用地不设发展限制, 同时限制发展坡度设为25°, 即地形超过25°的区域不能演化居住用地元胞；(2) 基本保护情景：从保护湖泊湿地和基本农田的角度, 适当控制城市发展的速度, 在现在趋势发展情景设计基础上, 将研究区域内的水域设置60%的保护概率, 农业用地设置50%的保护概率, 同时限制发展坡度设置为25°；(3) 严格保护情景：该情景加强了对农田、水域、林地等生态用地的保护力度, 将水域的保护概率设为100%, 农田和林地保护概率设为90%, 同时将限制发展坡度设置为15°。

2 结果与讨论 2.1 土地利用变化的历史动态分析

根据各用地类型对应的面积计算出2007年至2011年的土地利用转移概率(表 2), 分析土地利用转移概率矩阵可以得到, 农业用地和荒地的耗减度比较大, 耗减的农业用地向居住用地和荒地的转化占很大比例, 22.48%的荒地被开发成新的住宅用地。居住用地总体上的净增加量为10.34 km², 同时从表 2中可以得到2007年到2011年期间已有的居住用地的耗减度为13.18%, 原因是城区的一些破旧的城中村和老工业区的旧住宅被拆除, 转变成暂时的荒地(空地), 表明历史年份中居住用地元胞在增长的同时也有一小部分在衰败和死亡。研究区域内水域的整体面积变化不大, 可能是越来越注重生态宜居环境的住宅开发模式合理的保护了水域的边界, 也可能与政府采取的“还湖、还渔”等政策密切相关。

新增的居住用地呈散点填充在中心城区的待建用地中以及城郊区的工业园周围的闲置用地上, 部分新增居住用地元胞围绕湖泊沿岸分布, 这也是武汉特殊的自然环境优势之处, 住宅开发模式向靠近江河湖泊的生态宜居环境方向发展。2005年至2011年研究区土地利用变化的主要特征是城市空间的扩张及其对农业用地的占用, 这可以说是改革开放后武汉市土地利用变化的主要历史形态。

2.2 模型校正与模拟精度检验

本文用处理之后的4期研究区城市居住用地数据作为训练样本对模型进行校正, 以获取城市增长最佳参数。模型会将训练结果与真实数据进行对比, 并且用Lee-Sallee指数来评价校正结果精度, Lee-Sallee形状指数是各控制年份的模拟增长结果与实际城市居住用地范围的空间匹配度, 用公式表示为：

(1)

S为Lee-Sallee形状指数, A为模拟的城市居住用地面积, B为实际的城市居住用地面积。

执行校正模块, 本文中校正的3个阶段的参数设置和检验结果如表 3所示。

3个校正阶段的Lee-Sallee值均达到0.82以上, 说明模型校正结果非常理想。最后获得研究区的增长系数组合为：3(扩散系数)、4(繁衍系数)、1(蔓延系数)、95(坡度阻力系数)和45(道路引力系数)。运行SLEUTH模型的测试(Test)模块, 以2005年为初始年, 以校正阶段提取的最佳系数组合作为模型模拟的初始参数, 对研究区2006年至2011年的土地利用变化情况进行模拟, 得到逐年的土地利用模拟结果(图 3)。对土地利用模拟结果从空间匹配性和面积精确度两个方面进行评估, 选择Lee-Sallee形状指数为主要指标评估模拟结果与历史数据在空间位置上的匹配性, 各控制年份2007、2009和2011年的Lee-Sallee指数依次为0.87、0.94、0.67, 均在0.6以上, 空间匹配度比较高；利用面积比指数对模拟的土地利用面积进行精度验证(表 4), 计算结果表明各用地类型的演化在数量上的模拟精度也比较高, 精度均达到79%以上, 交通用地、公共基础设施用地、工业用地、农业用地、水域和林地的模拟精度均达到92%以上, 其中水域的模拟精度最高。相比较而言, 荒地的模拟精度稍低, 可能受“旧城改造”和产业上“退二进三”的政策影响, 一些污染较大的老工业区和城中村的拆除导致临时空地(待建用地)增多, 而SLEUTH模型模拟过程并没有捕捉突变式的社会经济政策影响。

2.3 不同情景下的城市空间预测结果

运用模型模拟, 得到研究区3种情景下2012—2025年土地利用变化结果。2025年研究区土地利用格局见图 4, 各土地利用类型面积变化见图 5。从面积变化看, 各情景下居住用地均明显增加, 城市居住用地对其他用地类型均有不同程度的侵占, 对农业用地和荒地的占用最多。从扩展模式看, 总体表现为从旧城中心区向外不断蔓延, 区域内没有产生新的大规模居住隔离组团, 新增的居住用地元胞主要是在原有的原有城市化基础上向外围郊区蔓延扩展, 空间格局较破碎化；同时也有内部填充式增长, 表现为中心城区内部的空地出现呈散点状分布的居住用地元胞。此外, 在城市交通的推动下, 居住用地沿交通便利区域蔓延速度比较迅速, 青山区和洪山区外围的高速路、洪山区的白沙洲大道、光谷一路以及江夏区的武昌大道周围出现了一定量的新增居住用地元胞。江河湖泊沿岸的新增居住用地元胞也很明显, 显示出“道路引力”的效果。

从图 5可以看出, 不同情景下土地利用变化趋势差异较大。现状趋势情景下城市扩展呈无序蔓延状态, 尤其是城郊结合区土地利用变化尤为剧烈。城镇用地扩张迅速, 占用了大量的农田和荒地, 林地和水域也有相当程度的损失。居住用地增加最为显著, 2011—2025年, 将增加5673.6hm², 年均增长378.24hm²。农田损失量约为2742.84hm², 林地减少667.08 hm², 水域将减少546.12 hm², 这会加剧区域的生境破碎、耕地功能下降、水资源匮乏、湖滨湿地萎缩等生态问题, 加大城市发展与生态保护之间的矛盾。

基本保护情景下, 城市扩张的速度得到了一定程度的控制。居住用地将增加4575.6 hm², 年均增长量为305.04hm²。农田损失较为显著下降, 损失量为1781.64hm²。林地和水域等自然生态类型仍然有显著损失, 林地减少了635.76 hm², 水域减少270.36 hm²。从空间格局来看, 此种情景下的城市空间扩展空间形态也比现状趋势情景下稍显紧凑。

严格保护情景下, 由于加大了对各类湖泊湿地和基本农田的保护力度, 居住用地的大肆扩张得到适当的控制, 土地利用变化速度明显降低。居住用地增幅显著减小, 增加2583.72 hm², 年均增长量为172.25hm²。林地、水域等生态用地损失明显降低, 林地减少203.4 hm², 水域减少129.96 hm², 相比基本保护情景, 损失量分别下降了68.01%、51.93%。这显示出生态保护措施具有显著的调控效果, 对于限制建设用地的肆意扩张和保护水域、农田等各类生态用地及其生态系统服务能力可以起到积极作用。

3种情景下, 工业用地和公共基础设施用地的面积变化均不明显, 主要是由于这些用地类型规模不大, 历史变动也较小, 从而根据历史外推趋势得到的预测面积变化也不显著。

3 结论

在GIS和RS技术的支持下, 利用SLEUTH元胞自动机模型对武汉市武昌、洪山和江夏部分地区的城市空间演变进行模拟和情景预测, 揭示不同土地保护强度下的城市扩展趋势, 主要结论如下：

(1) 从模型模拟精度来看, 2007年至2011年的城市空间模拟结果与历史数据的空间匹配度较高, 形状指数和面积比指数均较大, 模拟精度较高, 表明SLEUTH元胞自动机模型通过本地化校正后可以应用在城市空间动态模拟中, 具有较强的适用性和可移植性。

(2) 对武汉2025年城市空间扩展的模拟结果表明, 各情景下居住用地均明显增加, 对其他用地类型均有不同程度的侵占, 显示城市居住用地增长在未来土地利用变化过程中占据主导地位。从空间扩展模式来看, 新增的居住用地元胞主要是在原有的城市化基础上向外围郊区蔓延扩展, 空间格局较破碎化, 同时辅以内部空地填充增长；此外, 城市交通和自然环境优越的湖泊水系对城市空间的扩展有很强的吸引作用。

(3) 从不同情景下的模拟结果来看, 3种情景下农业用地、水域、林地和荒地的面积均有所减少, 现状趋势情景和基本保护情景下农田和水域减少的幅度较大, 会继续加剧区域的生境破碎、耕地功能下降、水资源匮乏等生态问题, 说明现状趋势和一般保护两种情景下, 城市重要生态组分及其生态服务能力将受到较为显著的损害, 进而威胁城市生态系统健康和可持续发展。严格保护情景下模拟结果显示, 城市居住用地扩张的程度得到了明显的控制, 水域和林地得到了有效的保护, 对于重要的自然生态系统组分保护及其服务能力维持可以起到显著作用, 保障城市系统健康和可持续发展。

本文的研究结果可以为城市空间动态相关研究和武汉城市发展调控提供参考。此外, 城市空间的演变受多种因素的综合影响, 存在许多不确定性, 进一步更为深入的研究, 应当研究将城市交通发展、新区建设、功能区转移、城市发展规划等多方面的影响因素纳入模型, 揭示其对城市空间动态的影响机制并在此基础上进行未来情景预测模拟, 为城市可持续管理提供更为精确和明确的参考依据。