0 引 言
城市化进程和居民生活水平的提升促使小汽车成为城市居民日常出行的重要工具,然而其广泛使用引发了交通拥堵、环境污染和能源消耗等诸多问题。降低居民小汽车依赖并推动绿色出行的普及,已成为城市规划和交通管理的关键课题。城市建成环境对小汽车保有和使用的影响是交通地理学、城市规划学等领域的研究焦点。
国内外学者已围绕小汽车保有和使用开展广泛研究。Shao等
[1]基于梯度提升决策树模型探讨了建成环境对小汽车保有的非线性影响,结果表明,多中心化发展和建成环境密度提高可以有效抑制小汽车增长。袁晓瑞
[2]、Wu等
[3]、尹超英等
[4]从多尺度建成环境角度分析了其对居民小汽车保有量的影响差异,发现个人及家庭属性影响最为显著。Cao等
[5]基于梯度提升决策树模型揭示了建成环境与小汽车依赖之间存在阈值效应。王晓全等
[6-8]利用多级贝叶斯模型和条件自相关模型确定空间效应对建成环境与小汽车依赖之间的关系有影响。Yin等
[9]和Chowdhury等
[10]在建成环境方面的研究进一步证实了土地混合使用、紧凑型城市发展是抑制汽车出行的有效政策手段。
现有研究已深入探讨建成环境与小汽车保有和使用间的影响效应和机理,但未充分考虑不同时期社会经济环境变化对建成环境属性影响居民小汽车依赖行为的动态性
[10-12]。自2010年以来,杭州市经济飞速发展,城市交通进入快速发展阶段,地铁从无到有,形成包含12条线路的地铁网络,城市道路总里程实现了翻倍增长,城市拥堵状况得到了有效缓解,城市建成环境(尤其是交通设施状况)发生了巨大改变,这推动了居民出行行为的转变。因此,本文基于2010年和2023年杭州市居民出行调查数据,以小汽车保有量和出行占比表征居民的小汽车依赖行为
[13],运用轻量级梯度提升机(Light gradient boosting machine,LightGBM)方法构建非线性模型,旨在探究不同时期下随时期动态变化的建成环境对居民小汽车依赖行为的非线性影响,为城市规划与交通管理提供科学依据。
1 研究区域及数据描述
1.1 研究区域概况
本文聚焦于浙江省杭州市市区的核心区域,涵盖了西湖区、拱墅区、上城区、滨江区,以及余杭区、临平区、钱塘区、萧山区的部分区域,详细范围如
图1所示。截至2021年底,杭州市居民小汽车保有量已超过340万辆
[14]。尽管近年来杭州市地铁网络迅速扩张,公共交通设施不断完善,但交通拥堵问题仍然是城市发展中的一个重大挑战。因此,本文选取了杭州市的主城区这一交通问题最为突出的区域作为研究对象,旨在为改善交通状况提供更有效的交通需求管理政策建议。
1.2 数据描述
本文使用的是2010年和2023年杭州市居民出行调查数据。这些调查通过入户访问和网络调查的方式收集,在预发放和检验通过的基础上采用问卷调查的形式对居民进行分区域按比例抽样。具体而言,2010年的出行调查共回收了16 874份有效样本,而2023年的调查则回收了5 828份有效样本,总计获得了22 702个有效样本个体。社会经济属性包括受访者的性别、年龄、家庭收入和学历,其中,性别变量与2011年和2024年杭州市统计年鉴中变量比例差异均在3%以内,年龄变量差异均在5%以内。此外,问卷调查严格遵循学术伦理规范,所有被调查者均在充分知情的情况下签署了受试者知情同意书,且此调查已通过浙江大学科技伦理委员会的审批。
建成环境指的是城市中的建筑和空间布局,尤其是那些能够通过政策手段进行调整和优化的部分。本文主要关注居民居住地周边的建筑要素,并采用“5D”框架全面描述建成环境的特征,包括密度(Density)、多样性(Diversity)、设计(Design)、目的地可达性(Destination accessibility)和公交邻近度(Distance to transit)。具体定义和指标详见
表1。
兴趣点(Point of interest,POI)数据因其庞大的数据量、高精度的坐标位置以及丰富的建成环境信息而备受青睐。每条POI记录都包含名称、类型以及经纬度坐标等信息,能够精确地描绘出建筑环境和城市设施在地图上的具体位置和分布。因此,本文利用POI数据提取建成环境属性,以期更准确地反映实际建成环境的情况。
本文数据涵盖了13种不同类型的POI数据,包括餐饮服务、风景名胜、公共设施、购物服务、交通设施服务、金融保险服务、科技文化服务、商务住宅、生活服务、体育休闲服务、医疗保险服务、政府机构及社会团体、住宿服务。根据《城市综合交通体系规划标准》(GB/T 51328-2018),本文提取了交通设施服务POI中公交车站和地铁站的地理位置信息,计算居民居住地半径500 m缓冲区内的公交车站点数量和半径800 m缓冲区内的地铁站点数量,并以此表征公交邻近度指标。
密度维度由人口密度表征,这一指标依据2010年和2020年全国人口普查数据以及杭州市街道面积统计信息综合计算得到。多样性维度则通过土地利用混合度这一指标衡量。土地利用混合度是指居住地附近POI类型的多样性和丰富程度,用熵指数
[1]定量分析,具体计算公式如下:
式中:Li 为第i个缓冲区内的土地利用熵指数;Pij 为第i个缓冲区内第j种POI数量占总数的比例;Ni 为第i个缓冲区内包含POI的类型数量。
土地利用混合度的值在0~1之间,值越接近1,说明该区域的土地POI类型越丰富。
设计维度在本文中通过路网密度表征,包括快速路、主干路、次干路和支路的密度,根据路网级配比综合相加得到,路网数据来源于OpenStreetMap平台,并借助ArcGIS软件对缓冲区内各类型路网密度进行统计和分析。目的地可达性维度则用商圈易达性表征,即居住地到最近商圈的距离,通过百度地图开放平台进行最短路径计算,商圈共选用杭州市最大的10个商圈,包括武林商圈、湖滨商圈、吴山商圈、钱江新城商圈、滨江商圈、人民广场商圈、申花商圈、新天地商圈、下沙商圈、西溪商圈,取最近商圈距离的倒数作为指标。
2 非线性模型构建
梯度提升决策树(Gradient boosting decision tree,GBDT)模型是一种通过迭代训练弱分类器(决策树)构建最优模型的机器学习方法。它在工业界的多个领域得到了广泛应用,近年来在交通领域也获得了广泛的关注和应用
[1,3-5,13,15-19]。与传统的统计学模型相比,GBDT模型具有以下几个显著优势:①可以更好地拟合数据;②能够灵活处理各类型数据,无论是连续变量还是分类变量;③对异常值具有较强鲁棒性,抗干扰能力强;④可以揭示自变量与因变量之间的非线性关系,不需要自变量服从正态分布,解决了多重共线性问题。
此外,与其他机器学习模型(如支持向量机(Support vector machine,SVM)或神经网络(Neural network,NN))相比,GBDT模型也兼具树模型的优点,有更好的可解释性,可以通过输出特征相对重要性、非线性部分依赖图(Partial dependence plots,PDP)以及结合沙普利加性解释(SHapley additive explanation,SHAP)等方法解释自变量与因变量之间的非线性关系。
近年来,为了在保持模型准确率的同时提高GBDT模型的训练速度,并更好地处理大规模数据,微软亚洲研究院
[20]提出了LightGBM。相较于其他传统基于树的学习算法,LightGBM模型在训练速度、计算效率和预测准确性等方面均实现了显著提升,广泛应用于出行分类、交通方式预测、目的地识别等任务。LightGBM模型的核心在于其采用改进的直方图算法(见
图2)进行决策树学习。对于一维特征(如年龄),该算法通过离散化分箱生成梯度统计直方图,每个分箱存储梯度累加值和样本频数;对于多维特征,则采用基于梯度排序的特征并行化分割策略,通过互斥特征绑定技术实现高效多维空间划分。相较于传统GBDT算法,该方法通过直方图降维、基于梯度的单边采样(Gradient‑based one‑side sampling,GOSS)及直方图差加速技术,在保持分割精度的同时显著提升训练速度,并降低75%内存占用,有效解决了高维稀疏数据场景下的计算效率瓶颈问题。
此外,LightGBM模型还改进了分裂策略,传统GBDT采用的是深度分裂策略(见
图3(a))构建决策树,这种策略在每个层次上对所有节点进行分裂,导致计算效率不高。而LightGBM引入了叶分裂策略(见
图3(b)),在每一层选择分裂效益最大的节点从而更高效利用节点信息,减少无效的节点分裂,大幅提高了模型训练速度。
本文分别对4种树模型进行建模比较,结果发现LightGBM模型在解释方差、均方根误差、均方误差、相对绝对误差4个评价指标中均优于其他3个模型,如
图4所示。
3 结果及分析
本文构建了LightGBM模型,具体如
表2所示。总体而言,建成环境属性对小汽车保有和使用的影响显著,但从2010年到2023年建成环境的影响显著减弱,从2010年的30.75%、18.91%下降至2023年的16.19%、10.86%,这与其他学者
[6,7,15]的研究结果相符。其中,人口密度因素在2010年对居民小汽车依赖影响最大,而道路密度因素在2023年对居民小汽车依赖影响最大。
3.1 单因素非线性分析
为了深入分析各因素对居民小汽车依赖的具体影响,本文绘制了部分依赖图,如
图5和
图6所示。
在小汽车保有量分析中,
图5(a)~(c)反映了不同时期社会经济因素对居民小汽车保有量的不同影响,
图5(d)~(h)则反映了不同时期建成环境因素对小汽车保有量的影响。由
图5(a)可见,2010年和2023年的年龄因素对居民小汽车保有量影响不同,在2010年,40~50岁的居民小汽车保有量最多,且60岁以上居民明显更少保有小汽车;在2023年,30岁左右的居民小汽车保有量最多,且在各年龄段分布更均衡,可见随着时代变化,居民小汽车保有量呈均衡化和年轻化。
图5(b)(c)反映了相较于2010年,在2023年不同学历和家庭收入群体中小汽车保有量的差距缩小。由
图5(d)可知,在2010年人口密度因素与小汽车保有量呈负相关,居住地人口密度越大,居民的小汽车保有量越少;而在2023年,不同人口密度区域居民的小汽车保有量无明显区别,人口密度因素的影响显著降低。由
图5(e)可看出,2023年土地利用混合度更高,但对小汽车保有量的影响更小。由
图5(f)可见,2010年和2023年商圈易达性因素均与小汽车保有量呈负相关,但2023年商圈易达性因素的影响更小。由
图5(g)可见,不同时期的差别主要在高道路密度区域,2010年高道路密度区域居民小汽车保有量较低,而2023年则与其他区域居民相比无明显降低。在2010年,当出发地周围超过10个公交或地铁站点时,公交邻近度对小汽车保有量的影响显著降低,而在2023年公交邻近度却有一定升高,这说明公共交通的建设能降低居民小汽车保有量,但当整体公共交通设施建设完善后,继续建设公共交通设施对居民小汽车保有量的降低没有帮助。
在小汽车使用占比分析中,
图6(a)~(c)反映了不同时期社会经济因素对居民小汽车出行占比的影响;
图6(d)~(h)则反映了不同时期建成环境因素对小汽车出行占比的影响。年龄因素在不同时期的影响差别不大,在2010年,28~45岁的居民小汽车出行占比较高,而在2023年,则为25~50岁,范围略有扩大。不同于2010年,在2023年不同学历居民小汽车出行占比区别不大,而家庭收入因素则与居民小汽车出行占比呈正相关,家庭收入越高,居民小汽车出行占比越高。在
图6(d)中,低人口密度区域居民的小汽车出行占比较高,在2010年拐点为10 000人/km
2,居住在人口密度高于10 000人/km
2的居民的小汽车出行占比显著较低,而在2023年拐点为5 000人/km
2,且高于拐点的居民仍有较高的小汽车出行占比,这说明低人口密度区域居民由于位置偏远更依赖小汽车,而随着经济发展,在2023年高人口密度区域居民的小汽车出行占比也较高。2023年土地利用混合度更高更均衡,对小汽车出行占比有负影响,土地利用混合度越高,居民小汽车出行占比越低,这说明随着城市各功能建设完善,居民出行需求可以更多在家附近实现,小汽车出行依赖会降低。在商圈易达性因素分析中,2010年商圈易达性越高,小汽车出行占比越低,而到2023年,随着商圈分布更多更均衡,对小汽车出行占比影响变小。在
图6(g)中,道路密度与小汽车出行占比呈负相关,且随着道路密度进一步完善,该负影响会变得更大,说明更合理、密集的道路建设对小汽车出行减少有帮助。
3.2 双因素交互分析
为了深入探究多个特征属性与居民的小汽车依赖行为的关系,本文基于SHAP方法,引入双因素交互分析,揭示不同特征属性间的交互影响。如
图7所示(彩图参见电子版,后文同),依赖图中横轴表示人口密度因素的取值,纵轴为人口密度因素与家庭收入因素交互作用下对模型预测结果的SHAP value(SHAP贡献值)。颜色表示家庭收入水平高低,红色越深值越大,蓝色越深值越小。灰色直方图反映样本量随人口密度变化的分布情况。
上文分析中发现,人口密度因素对小汽车保有量的影响在2023年显著减弱,因此本文结合家庭收入变化进行交互分析。
图7表明,2023年家庭收入较2010年显著提升,导致不同人口密度区域居住居民的小汽车保有量差异缩小。当家庭收入普遍达到一定水平时,居住地建成环境变量对小汽车保有量的影响逐渐降低,而其他社会经济因素的作用更突出。
在小汽车出行分析中存在一个有趣的现象,在2010年和2023年中,人口密度和道路密度因素均对小汽车出行存在负向影响,这与Milakis和Ding
[21,22]等的研究结论一致,但在2023年的交互作用分析图中有显著阈值效应,右上角存在明显的反趋势点。因此,本文进一步分析了道路密度、人口密度与社会经济属性对小汽车出行占比的交互影响,结果如
图8所示。经交互分析可知,居住地人口密度超过20 000 人/km
2或道路密度超过0.04 km/缓冲区单位的群体中,中老年群体的小汽车出行占比较低,而中高收入青年群体的小汽车出行占比较高。一般而言,青年群体对出行成本的敏感性较高,更高人口密度、道路密度的居住地有着更便利的公共交通条件,青年群体会减少其小汽车出行占比。在徐浩文等
[23]的研究中发现,高收入群体选择小汽车出行概率更高,其出行成本敏感性更低,因此产生本文现象的主要原因为青年群体拥有中高收入(即月收入大于10 000元),其出行成本敏感性显著下降且对出行舒适度要求显著上升,因而更多选择小汽车出行。因此,在公共交通建设中应加强公共交通的舒适性、便捷性,以提升其对青年群体的吸引力,降低其小汽车出行占比。此外,这也反映了随着经济的发展,居民小汽车出行使用受社会经济属性影响更大,而受居住地建成环境因素影响显著减弱,说明在2010年可以通过改变建成环境因素来调节居民对小汽车的依赖,但随着建成环境的逐步完善,再通过改变建成环境来调节居民对小汽车的依赖效果已不明显。因此,建议相关政策决策注意阈值效应和群体区分,更合理、均衡地规划建设建成环境,将更多注意力转移到交通需求管理中。
4 结 论
(1)本文构建的LightGBM模型精准捕捉了建成环境属性与居民小汽车依赖之间的非线性关系。不同时期建成环境因素对居民小汽车依赖的影响存在显著差异。
(2)2010年,人口密度是影响小汽车保有量及出行占比的主导因素;而2023年,道路密度取代人口密度成为关键影响因素,且人口密度的影响显著降低。在2010年,低人口密度区域居民由于位置偏远更依赖小汽车出行,随着经济的发展,在2023年,高人口密度区域居民的小汽车出行占比也较高。道路密度与小汽车出行占比呈负相关,合理、密集的道路建设有助于减少小汽车出行。
(3)与2010年相比,2023年建成环境属性对居民小汽车依赖的影响程度明显减弱,居民小汽车保有量和出行占比也得到有效控制。当家庭收入普遍达到一定水平时,居住地建成环境变量对小汽车保有量的影响逐渐降低,而其他社会经济因素的作用更突出。随着建成环境的完善,尤其是地铁等公共交通的发展,居民小汽车依赖的增长得到了有效抑制。
(4)有趣的发现是,2010年高人口和道路密度区域居民的小汽车出行占比较低,而在2023年,该区域中高收入青年群体的小汽车出行占比却呈现出明显的阈值效应,即有较高的小汽车出行占比。这表明随着建成环境的逐渐完善,应从建成环境的规划建设逐步转向交通需求管理等,分析不同人群出行需求,准确把握阈值效应的影响,以避免建设资源浪费及过度建设带来的负向效应。
国家自然科学基金重点项目(52131202)
浙江省“尖兵”、“领雁”研发攻关计划项目(2023C01240)