21世纪以来, 我国城镇化发展水平较快, 常住人口城市化率从2000年的36.2%增加到2018年的59.58%, 年均增长超过1%。城镇化快速增长给区域生态环境承载力带来巨大压力, 人口的剧增、工业化进程的加速、经济总量的飞跃式发展、经济增长方式的粗放型特征引发了严峻的生态环境问题, 同时区域生态环境恶化也制约城镇化的进一步发展, 阻碍城镇化发展质量的提升, 二者是相互影响、交互胁迫的耦合关系。面对快速城镇化引发的生态环境问题, 2014年出台《国家新型城镇化规划(2014—2020)》, 明确提出建设“以人为本”的中国特色新型城镇化道路, 党的十八大和十九大报告也高度重视“生态文明建设”发展战略, 因此城镇化与生态环境的耦合协调研究已经成为国家进行生态文明建设的重要内容之一。城镇化如何在生态环境承载力范围内实现合理、有效、优质发展, 生态环境如何更好的为城镇化发展提供良好空间和基础载体, 成为区域可持续发展的关键问题。开展城镇化与生态环境交互耦合效应的研究, 已经成为国际上未来10年地球系统科学与可持续性科学研究的热点与前沿领域^[1-5]。

国外学者David J.Rapport和Tony Friend 1979年提出PSR模型, 即从“压力-状态-响应”角度研究经济活动与生态环境之间的相互关系^[6]; Pearce认为在城市发展不同阶段所采取的环境政策不一样, 提出“城市发展阶段环境对策模型”^[7]; Grossman和Krueger通过研究人均收入与环境指标的关系认为, 城镇化与生态环境呈现出倒“U”字型的环境库兹涅次假设(EKC)^[8], 并判定不同地区的耦合状态^[9]。国内相关研究发展较晚, 主要集中在城镇化与生态环境耦合理论^[10-11]、机制^[12-14]、相互关系^[15-16]、耦合评价^[17-18]、规律与协同路径^[19-20]等方面展开, 研究方法主要运用双指数曲线^[21]、综合响应度模型^[22]、灰色关联模型^[23]、耦合协调度模型^[24-25]、BP模型^[26]、SD模型^[27]、VAR模型^[28]。研究尺度集中在国家^[29-30]、省级^[31-34]、城市^[35-36]、县域^[37], 马海涛^[38]对不同尺度下城镇化与生态环境关系进行了研究评述, 刘海猛^[39]构建了城镇化与生态环境耦合机理的分析框架-耦合魔方, 提出城镇化与生态环境耦合的4个分析维度——时间、空间、表象、组织, 拓展了近程耦合、远程耦合、近期耦合、远期耦合、组内耦合、组间耦合、显性耦合、隐形耦合8种耦合机制。崔学刚^[40]总结了4种城镇化与生态环境耦合动态模拟模型：系统动力学模型、人工智能算法模型(ANN和BN)、土地利用变化模型(CLUE、CA和MAS)、多模型集成复合模型。但是国内以城市群尺度进行研究的相对较少, 部分成果也多集中在对国家级重点城市群进行有效的探索, 如京津冀城市群、长三角城市群、珠三角城市群, 针对国内其他城市群的关注较少。方创琳^[41-42]对中国城市群取得的重大进展进行了梳理, 提出开展特大城市群地区城镇化与生态环境交互耦合效应的研究, 是未来10年地球系统科学研究的前沿领域和高优先研究主题, 并从理论上解析了特大城市群地区城镇化与生态环境的交互耦合效应, 创建了特大城市群地区城镇化与生态环境的交互耦合效应解析的理论框架及技术路径。王少剑^[43]对京津冀城市群地区城镇化与生态环境耦合关系和动态变化进行定量测度, 梁龙武、王振波^[44]进一步分析了京津冀城市群内部城镇化与生态环境耦合关系的空间分异和协同发展格局, 陈肖飞^[45]从利奥波德大地伦理观思想出发, 对长三角城市群新型城镇化与生态环境承载力耦合协调研究, 崔木花^[46]以中原城市群9个城市为例研究城镇化与生态环境耦合协调关系, 还有部分学者对特大城市群的生态安全^[47]、大气污染^[48]、生态效率等特征进行分析, 均取得了一定的研究成果。中国有19个城市群, 由于历史发展、资源禀赋、产业基础、经济结构等条件各不相同, 每一个城市群的城镇化发展水平存在较大差异, 生态环境条件和生态问题更是复杂多样, 导致二者出现不同的耦合特征, 因此对于每一个城市群都要结合自身特点进行研究。

本文在借鉴相关研究成果的基础之上, 以山东半岛城市群为例, 探讨城镇化与生态环境的时空演变及耦合格局。山东半岛城市群位于京津冀城市群和长三角城市群中间地带, 是国家主体功能区规划中的优化开发区, 是国家重点建设的沿海城市群之一。最初有济南、青岛、淄博、潍坊、威海、烟台、东营、日照8个沿海城市, 2014年规模扩大到13个城市, 2017年出台《山东半岛城市群发展规划(2016—2030)》, 扩大到17个城市。国家在“十三五”规划中明确提出“优化提升东部地区城市群, 提升山东半岛城市群开放竞争水平”, 打造北方重要开放门户, 打破行政壁垒, 构建“两圈四区”的城市群空间格局(图 1)。截止到2018年底, 山东半岛城市群常住人口10047万人, 地区生产总值76469亿元, 常住人口城市化率达到61.18%, 比全国平均水平59.58%提高近2个百分点。快速城镇化对区域人口容量、资源承载力、环境承载力产生巨大的压力, 传统资源产业的快速发展给生态环境带来一定的挑战, 城镇化与生态环境协调发展问题越来越受到关注。基于以上分析, 本文以山东半岛城市群为例, 研究2000—2016年城镇化与生态环境发展特征和耦合时空格局, 探索影响二者耦合协调发展存在的问题, 以期对山东半岛城市群城镇化与生态环境协调发展提供一定的建议, 为其他同类城市群城镇化与生态环境协调发展提供借鉴, 推动城市群健康绿色发展和区域经济高质量发展, 带动全国城镇化水平和生态水平的进一步提升, 对于丰富城市群的可持续发展研究具有一定的理论意义与实践意义。

1 研究方法与指标体系 1.1 研究方法 1.1.1 耦合协调度模型

耦合度源于物理学, 指多个系统或要素之间彼此相互作用、相互影响的程度。本文研究城镇化与生态环境两个系统之间的耦合协调关系, 借鉴耦合概念模型得出两个系统相互作用的耦合度模型, 公式如下：

(1)

(2)

式中, C为系统耦合度, f(a)是城镇化发展水平综合评价指数, f(b)是生态环境综合评价指数, f(a)和 f(b)可以通过多指标综合评价模型计算得出, T为城镇化发展水平与生态环境水平系统总得分, 反映城镇化与生态环境的整体协同效应。n为调节系数(n≥2), 本文取n=2。由于城镇化与生态环境的相互关系中, 二者同等重要, 因此取α=β=0.5。C∈[0, 1], 当C = 1时, 表明城镇化和生态环境系统位于最佳耦合协调状态; 当C=0时, 表明城镇化和生态环境系统基本处于不耦合状态。

耦合度的高低并不能反映两个子系统的发展水平。当生态环境与城镇化都比较低时, 也能有较高的耦合度, 但这种高耦合数值与高水平耦合的内涵不同, 因此构建能反映生态环境与城镇化协调发展水平的耦合协调度模型:

(3)

D为系统耦合协调度, D的数值介于0和1之间, 耦合协调度越接近1, 说明系统之间的耦合协调程度越高, 越接近0, 说明系统耦合协调度越低。

1.1.2 耦合协调类型

结合已有研究成果^[43-44], 综合耦合协调度和子系统综合评价指数得出城镇化与生态环境耦合协调发展类型及特征, 如表 1所示。

1.2 指标体系

城镇化是农村人口转化为城镇人口、农业用地转化为城市用地的过程, 不仅仅是人口和土地的空间变化, 还有人口空间移动所带来的城市经济结构、社会服务、生活方式、城市文化等全方位的城镇化。中国新型城镇化明确提出建立“以人为本”的城镇化, 人口迁移到城市是否能够享受到和城市居民同样的福利待遇、生活水平、生存空间, 人口城镇化与就业岗位是否匹配成为制约新型城镇化发展的关键。生态环境是人类赖以生存的自然资源与自然环境系统, 主要包括水资源、大气环境、能源系统、土地资源、生态绿地系统等多方面。生态环境的优劣取决于区域某一段时期内人类的生产与生活消费对资源环境产生的压力、资源禀赋的生态生产力以及生态环境治理和保护力度。结合已有研究成果, 从人口城镇化、经济城镇化、社会城镇化、空间城镇化4个方面构建城镇化系统指标体系, 从生态环境压力、生态环境状态、生态环境保护3个方面构建生态环境系统指标体系, 具体指标如表 2所示。

指标权重的确定方法可以分为主观赋权评价法和客观赋权评价法, 常用主观赋权评价法有AHP法、Delphi法、模糊评价法等, 常用客观赋权评价法主要有熵值法、主成分分析法、多目标规划法、TOPSIS法、变异系数法等。不同的评价方法都有一定的科学性和局限性, 本文在确定指标权重时采用了两种方法：主观赋值层次分析法(AHP)和客观赋值熵值法, 分别计算二者权重求平均值, 有利于消除单一方法所带来的弊端, 通过数据处理之后的结果如表 2所示。

1.3 数据来源

本文时间序列以2000—2016年为研究阶段, 空间格局以2000年、2005年、2010年、2016年为截面, 数据主要来源于《中国城市统计年鉴(2001—2017)》、山东省国民经济和社会发展统计公报(2001—2017)、山东省城市国民经济和社会发展统计公报(2001—2017)、《山东省城市统计年鉴(2001—2017)》、《中国统计年鉴(2001—2017)》、《山东省统计年鉴(2001—2017)》、《中国环境统计年鉴(2001—2017)》、《中国农村统计年鉴(2001—2017)》、《中国国土资源公报(2001—2017)》, 个别缺失数据通过移动平均法补齐。

2 实证结果分析 2.1 山东半岛城市群城镇化与生态环境耦合协调的时间演变特征 2.1.1 城镇化水平

2000—2016年, 山东半岛城市群城镇化综合水平一直保持着持续上升的趋势(图 2), 城镇化水平得到有效发展。2003年出台山东半岛城市群规划, 2007年提出“一体两翼”区域发展战略, 北翼黄河三角洲作为新的区域增长极, 同时作为欠发达地区的南翼鲁南经济带区域发展规划也正式实施, 2009年黄河三角洲高效生态经济区发展规划正式上升为国家战略, 2011年山东半岛蓝色经济区发展规划上升为国家战略, 正式提出“一蓝一黄”发展战略, 2013年提出“一圈一带”规划, 即省会城市圈发展规划和西部经济隆起带发展规划, 所有这些规划的颁布与实施加快了山东半岛城市群城镇化发展速度, 全面提升了城镇化发展水平。

通过图 2可以进一步分析城镇化子系统的变化趋势情况, 多年平均值表现为空间城镇化>人口城镇化>经济城镇化>社会城镇化。2000年到2005年, 经济城镇化略大于人口城镇化和空间城镇化, 但是幅度较小。从2006年到2008年, 人口城镇化开始超过空间城镇化, 空间城镇化超过经济城镇化, 人口增长进入高速发展阶段, 城市地域面积扩大, 大量人口从乡村涌入到城市, 导致人口城市化率快速提高, 同时由于经济增长开始从增长速度向增长质量转变, 导致经济增长率下降, 特别是2008年、2009年受全球金融危机的影响使经济城镇化处于最低点。2009年到2014年, 城市土地扩张进入快速发展期, 城市新区、高新产业园区、城市工业园区、经济技术开发区的建设导致城市用地面积不断扩大, 人均城市建设用地面积不断增加, 造城运动盛行, 空间城镇化逐渐超过人口城镇化。经济城镇化在这一时期随着经济增长率的持续下降增长缓慢, 2008年GDP增长率为20%, 2009年下降到9.6%, 2014年继续下降到7.6%, 初步改变了长期以来依靠资源型产业、重工业带动经济增长的局面, 为实现经济增长方式的变革而放慢了经济增长的速度。2015年到2016年, 人口城镇化超过空间城镇化, 2014年出台了《国家新型城镇化规划(2014—2020)》, 提出建设“以人为本”的新型城镇化, 严格控制城市建设用地扩张规模, 避免城市“摊大饼”式的空间增长, 抑制了土地城镇化的发展速度。但是随着我国2015年陆续放开城市人口“二胎”政策, 使城市人口自然增长率增加, 人口总量不断提高, 从而超过空间城镇化。社会城镇化子系统的贡献度最低, 但是也呈现出不断上升的趋势, 提高幅度是四个子系统中最大的, 说明山东半岛城市群在不断完善医疗、卫生、文化、公共交通等各项基础设施和服务功能方面做了很大努力, 社会城镇化水平也不断提高。

2.1.2 生态环境水平

2000—2016年山东半岛城市群综合生态水平呈现出两次波动中不断上升的趋势(图 3), 这说明山东半岛城市群生态环境出现过紧张状态但总体上不断走向优质。第一次波动是在2003年综合生态指数开始上升, 明确提出“生态省”建设, 生态环境发展水平上了一个台阶, 2006年指数下降是由于2006年人均水资源降到214m³, 降水量出现极低值。第二次波动上升是在2009年, 综合生态环境水平出现骤增, 主要由于森林城市创建取得一定进展, 2009年威海创建成为区域第一个国家森林城市, 森林覆盖率从2008年的13%提高到2009年的16%, 提高了区域生态资源容量, 临沂、淄博、枣庄、济南、青岛、泰安也陆续创建成为国家森林城市, 缓解了生态压力带来的负面影响, 随后生态环境水平缓慢上升。

进一步分析生态环境子系统的变化趋势如图 2所示, 2000年—2016年, 生态压力曲线表现为波动中下降的趋势, 说明在这16年间随着城镇化发展水平的快速增长和人口总量的快速增加, 对资源的挖掘和消耗日益严重, 污染物的排放和资源的消耗逐渐增大, 给区域生态环境造成一定的压力, 生态环境问题比较严峻。从2012年开始, 生态压力曲线开始提升, 生态压力现状才开始有所缓解。生态状态和生态保护曲线都一直呈现出不断上升的趋势, 主要由于山东半岛城市群提出“生态山东”建设, 推进绿色发展、循环发展、低碳发展, 生态资源开发、生态环境治理、生态环境保护以及产业结构转型升级大大改善了区域生态环境的总体水平, 同时居民对生态环境的保护意识逐渐加强以及政府相关部门对生态环境问题关注度不断提高, 环保投入增加、资源生态生产力增强, 逐渐形成山东半岛城市群生态建设的大格局。工业污染物的排放和技术处理、城市绿地面积和森林面积是制约山东半岛城市群城市生态环境优劣的主要影响因素, 政府要重点考虑环境污染和环境资源、环境治理之间的关系。

2.1.3 耦合协调度

根据2000—2016年山东半岛城市群城镇化综合评价指数和生态环境综合评价指数, 计算得出城镇化与生态环境耦合协调度D, 其变化趋势如图 4所示：2000—2016年山东半岛城市群耦合度和耦合协调度都不断提高, C值从0.4提高到0.9以上, 耦合协调度由基本不协调向高级协调发展转变, 子系统协调特征从城镇化滞后型向生态环境滞后型转变。

2000—2016年, 山东半岛城市群城镇化和生态环境耦合度表现为“快速增长-保持稳定”的发展特征, 耦合协调度由0.196缓慢上升至0.914, 经历了严重不协调、基本不协调、基本协调、中度协调、高级协调5个阶段。2003年之前耦合协调度不高, D值处于0.2—0.4之间, 2003年进入基本协调阶段, 2005年进入中度协调阶段, 2010年之后基本稳定在0.8—0.9之间, 处于高水平协调阶段, 说明山东半岛城市群城镇化与生态环境耦合协调越来越密切。具体来说可以分为以下几个时期：

第一阶段为不协调阶段(2000年—2001年), 耦合协调特征为不协调-城镇化滞后阶段。2000年D值为0.196, 处于严重不协调阶段, 但是已经接近基本不协调, 2001年D值是0.381, 完全进入到基本不协调阶段。生态环境系统综合水平要远远高于城镇化系统综合水平, 21世纪初城镇化发展速度还没有进入到高速增长阶段, 対生态环境系统还没有构成严重的威胁, 城镇化发展和生态环境的耦合关系并不明显。

第二阶段为转型阶段(2002年—2005年), 耦合协调度不断提高, D值由2002年0.434增长到2005年0.591, 由不协调-城镇化滞后阶段进入到基本协调-同步阶段。城镇化系统发展水平开始迅速提高, 逐渐追平生态环境系统发展水平, 在城镇化发展过程中开始注意生态环境的保护与治理, 生态环境有所改善。

第三阶段为中级协调阶段(2006年—2009年), 耦合协调度继续提高, D值由2006年0.636增长到2009年0.772, 由基本协调-同步阶段进入到中度协调-生态环境滞后阶段。城镇化进入高速发展期, 大规模的人口城镇化导致人口短时间内大量集聚到城市, 给城市生态环境带来一定的压力和冲击, 导致生态环境系统发展水平出现暂时的回落。但是随着对生态环境污染治理投资的加大以及对生态资源的有效利用与开发, 生态环境综合发展水平也不断提高, 二者的耦合关系越来越密切, 但是生态环境综合发展水平始终落后于城镇化综合发展水平, 因此二者耦合协调类型从基本协调-同步型转向中度协调-生态环境滞后型转化。

第四阶段为高级协调阶段(2010—2016年), 耦合协调度继续提高, 上升速度减慢, 耦合协调度曲线逐渐放缓, D值稳定在0.8以上, 由中度协调-生态环境滞后阶段进入到高度协调-生态环境滞后阶段。城镇化发展水平继续提高, 生态环境综合发展水平在波动中也不断提高, 二者协调程度不断优化, 主要由于贯彻执行国家十八大、十九大陆续提出和深化“生态文明”建设以及“新型城镇化建设”的结果。但是后期城镇化发展与生态环境发展有拉大趋势, 生态环境出现一个下降, 说明高水平城镇化发展给生态环境带来的问题仍然十分严峻。

2.2 山东半岛城市群城镇化与生态环境耦合协调空间格局

以2000年、2005年、2010年和2016年为时间点, 对山东半岛城市群城镇化综合评价指数与生态环境综合评价指数进行截面耦合分析, 得到17个地级市城镇化发展水平、生态环境发展水平、城镇化与生态环境耦合协调度的时空格局演化特征。由于城镇化发展速度和生态环境发展速度不相同, 因此在进行发展水平测度时的标准略有不同, 城镇化发展速度较快, 划分为五级, 分别是高水平城镇化区[0.6, 1]、较高水平城镇化区[0.4, 0.6)、中等城镇化区[0.3—0.4)、较低水平城镇化区[0.2, 0.3)、低水平城镇化区[0, 0.2), 如图 5所示; 生态环境发展速度较慢, 划分为四级, 分别是生态环境优质区[0.6, 1]、生态环境中等区[0.5, 0.6)、生态环境较差区[0.4—0.5)、生态环境恶劣区[0, 0.4), 如图 6所示。

2.2.1 城镇化发展水平

从图 5可以看出, 首先, 绝大多数城市的城镇化发展水平不断提高, 但是城镇化发展速度不同, 城镇化发展水平空间差异有拉大的趋势, 即高水平城镇化地区城镇化发展速度快于低水平城镇化地区。济南城镇化水平从2000年的0.581提高到2016年的0.762, 增加了将进0.2, 而聊城从0.172提高到0.282, 只提高了0.1, 发展较为缓慢。第二, 高水平城镇化地区“双核”特征明显, 中等城镇化地区范围不断扩大, 较低城镇化地区逐渐缩小, 低水平城镇化区消失。高水平城镇化区分布范围逐渐扩大, 从2000年的东营不断蔓延到2016年的核心城市济南、青岛, 济南和青岛的“双核”特征明显, 形成了高城镇化的第一梯队。东营人口较少, 石油产业发达, 人均GDP在区域排名第一位, 2009年成为黄河三角洲高效生态经济区的核心城市, 城镇化发展水平高, 但2014年沾化县撤县换区使东营建成区面积扩大, 导致整体城镇化水平有所下滑。济南作为省会城市, 人口集聚能力、产业集聚能力、政策吸引能力强, 社会基础服务设施完善, 技术人才、城市建设和发展、资金支持等都走在了前列。青岛作为一座边缘港口城市, 承接山东半岛城市群对外发展功能, 交通区位优势和高端科技人才优势明显, 产业基础雄厚, 创新能力强。围绕济南都市圈的淄博和莱芜以及沿海的威海和烟台形成城镇化发展的第二梯队, 依靠“双核”的涓滴效应和自身的资源优势成为城镇化发展水平较高的地区。淄博人口城市化率达到70%, 是全国重要的建材基地和石化基地, 城市化发展水平较高。莱芜从中等城镇化水平上升到较高城镇化水平区, 钢铁工业在产业结构中占有较大比重, 但是单一的工业结构类型也给未来莱芜抗风险压力有所降低。威海和烟台属于东部滨海城市, 自然环境优越, 旅游产业、机械制造、化工产业高度发达, 是山东半岛蓝色经济区的重要核心城市。第三梯队主要集中在鲁南地区的枣庄、济宁、临沂以及日照、潍坊。济宁从较低城市化地区上升到中等城市化地区, 临沂从低水平城市化地区上升到中等水平城市化地区。临沂医学制药业、食品工业、建材工业、现代物流业高度发达, 但是临沂总人口1000万左右, 大部分分布在县级行政区划单位, 从人均水平上降低了该市综合城镇化水平。日照从低水平城镇化地区上升到中等水平地区, 主要受到青岛的辐射带动作用。城镇化水平较低的有聊城、德州、菏泽以及滨州, 从低水平城镇化区上升为较低水平区, 农业人口比重较大, 城市产业不发达, 今后应继续挖掘本地资源优势, 提高基础服务功能, 探索适合本地发展的新型城镇化道路。

2.2.2 生态环境发展水平

从图 6可以看出, 研究时间范围内整体生态环境水平逐渐提高, 高水平生态环境区主要集中在东部沿海地区, 中等水平生态环境区从胶济沿线向南北扩散, 低水平生态环境区主要集中在鲁西南、鲁北部分地区。东部沿海烟威都市区、青岛都市圈核心城市青岛、临日都市区日照生态环境优质, 具有先天的自然条件优势, 滨海风景优美, 第三产业高度发达, 居民的环保意识强烈, 政府对生态环境的监管和处置力度到位, 特别是威海空气质量、环境优良度在山东半岛城市群都处于首位, 构成了山东半岛城市群的东部绿色沿海防线, 泰安成为中部生态环境最为优质的城市。济南都市圈中的济南、莱芜、淄博、德州、潍坊、济枣菏都市区的济宁、枣庄、临日都市区的临沂、东滨都市区的东营组成了中等水平生态环境区。济南虽然是省会城市, 但是长期钢铁、机械、化工产业的发展给城市生态环境带来了一定的问题, 工业烟尘和固体废弃物的排放较多, 政府对环境要给与充分的重视, 今后要不断增强省会中心的人口凝聚力和高新技术产业凝聚力。莱芜产业结构以重工业和煤能源产业为主, 钢铁产业占整个工业结构的40%, 导致城市的工业烟尘、工业固体排放物、工业SO₂排放量严重超标, 但近年来生态环境质量有所改善。临沂和枣庄长期以来一直处于生态环境低水平状态, 临沂陶瓷建材工业和枣庄煤炭工业给当地生态环境造成巨大的破坏, 但是随着枣庄煤炭工业向煤化工产业的转型、台儿庄古城的旅游开发以及临沂现代物流商贸业的发展, 区域产业结构和经济增长方式不断优化升级, 生态环境逐渐好转, 形成生态环境发展水平的第二梯队。菏泽、聊城、滨州的生态环境质量最差, 产业结构以工业为主, 政府相关机构环保长效机制不健全, 监督和奖惩措施不到位, 污染物科学治理、排污、处置能力不够, 沿袭重经济发展、弱环境保护的传统老路。菏泽化工园区环境问题突出, 城市基础服务设施不完善。聊城工业占国民经济结构的60%, 是全国重要铝和铜有色金属基地, 化石燃料消耗严重, 热电燃煤所需要的火电释放大量的污染物。滨州工业废水和工业SO₂排放量严重超标, 人均排放量分别达到51.62t/人和404/人, 排污治污设施标准低。

2.2.3 耦合协调度

基于GIS分析软件平台, 根据耦合协调度区间划分标准, 得出研究期内山东半岛城市群城镇化发展与生态环境耦合协调度空间格局, 如图 7所示。

山东半岛城市群城镇化与生态环境整体耦合协调性不断向前推进, 空间差异明显, 呈现出东部沿海地区、鲁中地区向鲁南、鲁西地区逐渐降低的分布趋势, 东部大于西部, 沿海大于内陆。山东半岛城市群经济发展基础和自然环境条件不同, 导致耦合协调度出现空间差异。首先, “双核”耦合协调度高, 以青岛为核心的山东半岛蓝色经济区是城镇化与生态环境耦合协调度整体水平最高的地区, 青岛和威海从2000年中度协调走向2016年高级协调, 城市化水平和生态环境水平发展程度都很高, 日照从基本协调上升到中度协调, 烟台一直保持在中度协调状态。“双核”中的另一核——以济南为核心的济南都市圈耦合协调度较高, 济南从中级协调进入到高度协调, 周边区域泰安、莱芜也从基本协调阶段进入到中度协调阶段。第二, 中度协调城市较多, 包括东营、淄博、烟台、潍坊、济宁、泰安、日照、莱芜、临沂、德州。东营作为黄河三角洲高效生态经济区的核心城市, 城市化发展水平高, 但是在发展初期一度环境污染非常严重, 后期虽在“黄河三角洲高效生态经济区”规划不断实施的过程中, 生态环境质量有所改善, 但是距离高效生态区的标准还有一定距离。枣庄、临沂、德州、莱芜从2000年基本协调进入到2016年中度协调, 资源枯竭型城市枣庄目前正处于城市旅游转型关键时期, 临沂市2013年入选首批“国家生态文明先行示范区”建设, 对于城市可持续发展和绿色赶超提供了生态文明创新之路。第三, 聊城、菏泽、滨州是城镇化与生态环境耦合协调度最低的地区。虽然摆脱了不协调的状态, 但在整个城市群仍处于最低的位置。菏泽和聊城产业结构以农业活动为主, 城镇化发展速度和城镇化发展质量都不高, 生态环境治理和保护意识较弱。滨州市的生态环境压力问题较为突出, 今后应提高工业技术, 从源头上降低污染物的排放。

2.2.4 耦合协调类型分析

根据表 1得出的山东半岛城市群城镇化和生态环境耦合协调类型的分类, 得出山东半岛城市群城镇化与生态环境耦合类型变化特征, 如表 3所示：

从表 3可以发现, 山东半岛城市群城镇化与生态环境耦合协调类型从城镇化滞后型向同步型、生态环境滞后型开始转化, 但落后地区的城镇化质量和环境问题仍不容乐观。2000年, 威海等13个城市为城镇化滞后, 济南、淄博、莱芜为均衡发展, 东营为生态环境滞后; 2005年, 烟台等10个城市为城镇化滞后, 济南、青岛、淄博、东营、莱芜5个城市为生态环境滞后, 枣庄、日照为均衡发展, 城镇化快速发展对生态环境胁迫作用增强; 2010年, 生态环境的紧张状况有所缓解, 济南仍然维持生态环境滞后, 青岛等4个城市为均衡发展, 其他城市为城镇化滞后; 2016年, 济南、东营、滨州为生态环境滞后, 青岛等6个城市为均衡发展, 其他7个城市为城镇化滞后。从数量上来说, 在这16年的发展过程中, 城镇化滞后的数量由2000年的12个减少到7个, 均衡发展的数量从2个提高到6个, 生态环境滞后型由1个增加到3个, 因此总体上山东半岛城市群城镇化发展水平有所提高。生态环境滞后型虽然只有东营、济南、滨州3个城市, 但是其环境污染程度和生态环境压力是巨大的, 滨州2016年人均SO²排放量404.77, 比青岛14.02的28倍。落后地区的城镇化水平和生态环境问题都不容乐观, 特别是聊城、菏泽的耦合协调度最低, 城镇化发展水平和生态环境发展水平都处于最低状态。菏泽、聊城等欠发达地区的人口城镇化发展水平较快, 从2000年的17%提高到2016年的45%左右, 但是经济城镇化和社会城镇化发展速度慢、水平低, 城镇化发展质量落后于城市化发展速度, 抑制了城镇居民生活水平和享有社会基础设施服务水平的提升。落后的城镇化水平刺激了区域的大面积经济开发, 但由于经济基础薄弱, 技术创新能力低端落后, 导致整体生态环境不断恶化, 因此如何共同提高落后地区的城镇化与生态环境综合发展水平是政府相关部门今后要开展工作的方向。

发达城市的核心作用和集聚作用明显, 对周围地区形成强大的吸引力成为区域增长极。济南的人口城镇化和社会城镇化高度集聚, 2016年人口城镇化指数0.843, 社会城镇化指数0.921, 比青岛还多0.2个指数点, 东营的经济城镇化水平高, 经济城镇化指数为0.818, 青岛的人口、经济、社会、空间城镇化水平比较均匀, 综合城镇化水平较高, 发达城市成为带动山东半岛城市群城镇发展的主要区域。因此城镇化发展水平的空间差异较大, 落后地区与发达地区之间的城镇化发展指数差距较大, 这也很好的解释了即使部分城市如聊城、菏泽、临沂是城镇化发展滞后类型, 但是整个山东半岛城市群的综合城镇化水平仍然高于综合生态环境发展水平。生态环境发展水平同样如此, 威海的生态环境在每一年都是最优的, 滨州、菏泽、聊城、莱芜的生态环境最恶劣, 2016年也是刚刚脱离了基本不协调类型。济南虽然城镇化发展水平高, 但是一直处于生态环境滞后型, 主要由于济南人口高度集聚, 资源消耗大, 重工业产业结构导致城市环境问题突出。

3 结论与讨论

根据2000—2016年的面板数据, 借助耦合协调度模型, 研究山东半岛城市群城镇化与生态环境耦合协调发展的时空格局演变, 得出以下结论：

(1) 2000—2016年山东半岛城市群城镇化综合水平高速发展, 子系统发展特征表现为多年平均值空间城镇化>人口城镇化>经济城镇化>社会城镇化的演化趋势, 生态环境综合发展水平呈现出波动中不断上升的趋势, 生态压力子系统曲线先降后升, 生态压力有所缓解, 生态状态和生态保护曲线一直处在上升状态, 总体生态环境不断走向优质。

(2) 2000—2016年山东半岛城市群城镇化与生态环境耦合度和耦合协调度都不断提高, 逐渐走向优质协调, 依次经历了严重不协调、基本不协调、基本协调、中度协调、高级协调五个阶段, 实现了耦合协调性从严重不协调向高级协调的转变, 子系统耦合协调特征经历了四个阶段：不协调-城镇化滞后、基本协调-同步型、中度协调-生态环境滞后、高级协调-生态环境滞后, 从城镇化滞后型向生态环境滞后型转变。

(3) 分析2000年、2005年、2010年和2016年山东半岛城市群城镇化发展水平与生态环境发展水平的时空演化, 得出：城镇化水平不断提高, 但空间差异逐渐拉大, 高水平城镇化地区“双核”特征明显, 中等城镇化地区范围不断扩大, 较低城镇化地区逐渐缩小, 低水平城镇化区消失; 高水平生态环境区主要集中在东部沿海的威海、烟台和鲁中的泰安, 中等水平生态环境区从胶济沿线向南北扩散, 低水平生态环境区主要集中在鲁西南、鲁北和鲁中部分地区, 生态环境恶劣区不存在。

(4) 分析2000年、2005年、2010年和2016年山东半岛城市群城镇化与生态环境系统耦合协调度的时空演化, 得出山东半岛城市群城镇化与生态环境耦合协调度空间差异明显, 整体耦合协调性不断向前推进, 耦合协调度呈现出由东部沿海地区、鲁中地区向鲁南、鲁西地区逐渐降低的分布趋势。青岛和济南由中度协调转化为高级协调; “双核”周边的城市以中度协调为主, 鲁西南地区的聊城和菏泽是城镇化与生态环境耦合协调度最低的地区。子系统耦合协调类型由城镇化滞后型向同步型、生态环境滞后型开始转化, 发达城市的集聚作用明显。

为了更好的推动山东半岛城市群城镇化与生态环境之间的协调发展, 首先要树立与生态文明建设相适应的新型城镇化发展理念, 强化绿色、低碳、环保发展意识, 维护和扩大绿色生态空间, 增强区域生态环境承载力, 出台相关政策支持节能减排、严把环保标准、淘汰落后产能, 转变城市经济增长方式, 促进产业结构升级, 探索绿色城镇化的转型之路。同时在国家新旧动能转换背景下, 积极培育特色小镇, 推进农民就地城镇化, 为促进山东半岛城市群绿色城镇化水平的提高提供了空间和载体, 从而形成功能定位科学合理的生态战略格局, 促进山东半岛城市群可持续发展。

城镇化和生态环境耦合协调发展是目前国内外学术界研究的热点, 也是我国进行生态文明建设和新型城镇化建设的核心问题之一, 但是目前国内以城市群尺度进行城镇化与生态环境协调发展的研究尚处于探索阶段, 城市群和城市、区域相比是一个更复杂的系统, 城市群城镇化系统和生态环境系统的庞大性、复杂性、特殊性更为突出, 因此本研究在系统指标体系的构建上难免存在不足和缺陷, 如城市群研究关注城市之间的经济联系、交通联系、社会联系等关联因素受到指标获取性的影响考虑不全面, 今后要不断延伸和丰富城市群城镇化和生态环境指标体系。但是本论文在进行耦合时空特征及相关原因分析时尝试从城市群关联的角度进行分析, 并且对城镇化子系统和生态环境子系统的变化特征也进行了充分的论证, 对于深入剖析二者之间的耦合协调关系具有一定的意义, 也对其他类似城市群的发展提供一定的思路。今后可继续围绕城市群城镇化与生态环境耦合协调的驱动机制、发展趋势和优化路径展开下一步的深入研究, 不断丰富城市群经济、社会、环境的可持续发展。