期刊佳文 | 混行环境下的交通碳排放测算与演化规律

目前，全球运输部门碳排放量占能源消耗产生碳排放的24%，其中道路交通约占运输部门碳排放量的74%［1，2］，道路交通碳排放已成为影响可持续发展的重大挑战之一。发展电动汽车已变为道路交通节能减排的重要手段，截至2022年中国新能源汽车保有量达到1149万辆，其中纯电动汽车占新能源汽车总量的80%以上［3，4］。出租车作为电动汽车推广的重要先驱，对推广新能源汽车和构建绿色发展公共交通体系具有重要作用。以出租车数据为基础，清楚把握电动汽车发展对碳排放的影响，分析其时空演化规律，可以为制定电动汽车减碳政策提供支撑［5～7］。

随着大数据时代的到来，多种数据被用于交通碳排放研究，比如使用手机信令数据量化城市拥堵造成的碳排放［8］、使用地铁IC卡数据估算地铁出行的碳减排量［9］和使用共享单车数据探索城市的减排路径［10］。在进行城市道路交通碳排放研究时，出租车GPS数据可以很好地反映城市出租车交通的动态变化，对基于GPS轨迹数据的出租车特性进行碳排放分析可以获得详细准确的结果［11～13］。传统的利用GPS数据计算碳排放的方法大部分只考虑燃油车或电动汽车并对某一区域进行碳排放计算，如Shang等［14］使用出租车轨迹数据及地理信息数据计算北京市燃油车辆碳排放［15］，Luo等［16］对上海市燃油出租车碳排放的时空分布特征进行分析，还有学者从全寿命周期角度出发研究电动出租车的碳排放［17，18］，但电动汽车发展是一个漫长的过程，之前的研究很少考虑在电动出租车与燃油出租车同时运作情况时的碳排放，以及探究其演化规律。

国际上碳排放测算常用的方法有实测法［19，20］、排放因子法［21］、质量平衡法［22］、生命周期法［23］和投入产出法［24］。其中，排放因子法是目前应用范围最为广泛、最为普遍的方法，也是交通运输行业碳排放测算的主要方法，包括基于燃料的“自上而下”方法［25～29］和基于行驶轨迹距离的“自下而上”方法［30～33］。“自上而下”的方法依赖于总体数据和统计模型，可以更快速地得出估计结果，但可能存在一定的误差［34］。“自下而上”的方法需要收集大量的具体数据，需要投入较大人力和物力资源，但能够更准确地估计碳排放量［35］。出租车作为城市常见交通工具，在各个地区车型比较固定，可以轻易获取各类车型的能耗数据，GPS数据中又包含地理位置信息、速度、载客状态和时间等信息，可精确得到出租车行驶距离与速度［36］，因此本文将两种方法进行结合，在自下而上计算时将距离速度等信息转换成与自上而下计算使用的燃料相关的参数，通过此方法可以将燃油车和电动出租车的碳排放计算结合起来，提高碳排放计算的精准度并减少计算的工作量。

综上所述，本文拟基于出租车GPS数据，探讨适用于燃油出租车和电动出租车混行环境下的交通碳排放测算方法，并在计算历年出租车碳排放水平的基础上，通过地图匹配将测算结果匹配至栅格，运用K-means算法对碳排放区域进行空间聚类与时空分布演化分析，解析电动出租车所占比例变化对碳排放的影响、探究碳排放与建成环境和各类用地之间的关系，以期为节能减排策略制定提供参考。

各年电动出租车占比数据来自西安市统计年鉴以及西安市人民政府公告，如图1所示。

可以看出西安市2016～2022年是一个完整的电动出租车发展过程，因此使用西安市2022年、2020年、2018年和2016年出租车的GPS数据与POI数据来进行碳排放测算与影响因素分析，GPS设备每隔30～35s进行一次记录，每天的数据量约为2500万条，出租车GPS数据主要参数如表1所示。出租车GPS数据可能存在因设备故障、信号干扰、司机操作失误等原因引起的数据缺失和数据错误，故需要首先对缺失数据和错误数据进行处理。处理过程如下。

（1）数据剔除先对数据中字段缺失部分进行剔除，之后通过经纬度和运营状态信息对定位异常与非运营状态的数据进行删除，对所有速度大于180km·h⁻¹的数据与城市内部行驶速度大于80km·h⁻¹的数据进行剔除。

（2）坐标转换由于年份不同数据采用的坐标系可能存在差异，需要对所有数据的坐标系进行统一，以国际上通用的WGS-84坐标系作为基准坐标系对所有数据进行统一。

（3）OD提取并处理异常出行根据载客状态中0与1之间的转变点进行车辆OD提取，对完整的OD段进行筛选进而确定出行产生的热点区域，运用数据中携带的经纬度信息将轨迹点匹配至地图上，对偏离路网的轨迹点进行修正。

为展示出租车空间分布特征，将出租车轨迹栅格化处理后进行地图匹配，其中栅格大小为250m×250m，生成出租车轨迹点分布热力图，如图2所示。

选取出行量较大并且出行状况复杂的地区，可以更明显地反映特征变化情况［37］。由图2可知，出租车轨迹点主要分布在西安市各大主干道上，热点区域主要集中在主城区内，因此本研究选取西安市主城区作为研究范围，包括雁塔区、碑林区、新城区、莲湖区和未央区。

建成环境因素在选择的时候一般遵循5D原则，即密度（density）、设计（design）、多样性（diversity）、距离（distance）和目的可达性（destination accessibility）［38，39］。基于此，围绕土地利用，兴趣点分布，交通特征等方面，在研究范围内选取9类变量，其中土地利用分为6类，分别为居住用地、商业用地、工业用地、公共设施用地、休闲娱乐用地和其他用地；除用地外，公共交通站点与出租车的替代关系会对交通碳排放产生一定的影响［40］；土地利用混合度由土地利用熵指数表征，如表2所示。

为了直观展示碳排放的演化情况，利用Python将出租车各个轨迹点连接成轨迹段，创建250m×250m的栅格，对研究区域进行栅格化，轨迹段与栅格相交，通过相交点将轨迹段进行打断，一个网格内可能会有多个轨迹段；之后通过对栅格中每一段轨迹的速度进行油耗率的求取进而对出租车碳排放进行测算后求和，得出各栅格中碳排放大小，按照各年出租车的比例来确定最终的碳排放；最后通过K-means聚类与梯度提升决策树（gradient boosting iterative decision tree model，GBDT）模型进行影响因素分析，具体流程见图3。

车碳排放测算方法提出一种将自上而下和自下而上相结合的出租，在自下而上计算时通过将距离速度等信息转换成自上而下计算时要用到的油耗率等燃料相关的参数，计算出各轨迹段的燃油消耗进而将其转化为碳排放，如式（1）所示。

根据IPCC国家温室气体清单指南，得到各种燃料对应的CO₂排放系数，通过式（2）计算燃油车的碳排放：

电动汽车的碳排放根据电动汽车能耗折算方法计算［42］，通过将电动汽车的能耗转化为同等条件下燃油车的能耗进行计算，转化公式如下：

（1）K-means聚类算法常用于处理大规模数据集，能在保持较高效率的同时很好地表示聚类特征。采用K-means聚类算法对碳排放区域进行聚类分析，步骤如下：选择k个聚类中心。

（2）计算各个栅格中心到聚类中心的间距，以距离和碳排放大小作为条件进行分配。

（3）更新每个簇的聚类中心，如果变化则重新进行分配，直到收敛。

（4）输出聚类结果，在开始之前，需要选择合适的距离度量方法和目标函数来计算聚类质心。

计算簇内误方差：

GBDT模型是一种基于迭代的集成学习算法［43～45］。它通过不断迭代来训练多个决策树模型，并利用前一个树的预测结果与实际值之间的残差来训练下一个树，逐步提升模型的准确性；通过信息增益考虑了特征在决策树中的贡献，信息增益是基于特征在决策树的分裂点上的影响来计算的，能够更好地反映特征对模型的贡献程度。因此，采用GBDT模型计算建成环境对碳排放的相对重要度影响，将部分碳排放数据作为训练集输入与建成环境指标一起进行模型训练，训练后将得到的特征重要度与测试集一同输入进行得分归一化与参数调整，最终得到各指标特征重要度，具体流程如图4所示。

将西安市各年电动出租车的占比代入电动出租车碳排放计算公式求出具体转化率；通过Python将每个栅格内的轨迹段长度与速度代入公式进行遍历计算，显示出每个栅格内具体的碳排放大小。得到各年碳排放分布如图5所示。

由图5可知，2016～2018年的碳排放整体较高，碳排放量呈下降趋势但下降幅度较小，约12%；2018～2020年碳排放量急剧下降，在2020年到达最低点，这是由新冠疫情刚暴发人们的出行量大幅度减少，并且西安电动出租车在此年数量急剧上升造成的；2022年碳排放量相对回升但相较于疫情前其他年份有明显的下降，此时疫情逐渐结束，出行量明显增加，与此同时燃油出租车逐渐被淘汰，电动出租车成为西安市出租车的主要组成部分，这使碳排放在出行量增大的情况下仅出现小幅度回升。

对所有栅格中的碳排放进行求和得到2016～2022年各年工作日与周末碳排放总量，并将其与电动出租车发展规律进行对比（图6）。可以看出，电动出租车的占比与碳排放大小呈负相关，2022年工作日碳排放相比于2016年降低了41%，周末碳排放降低了56%，并且随着电动出租车的占比增加，周末与工作日的碳排放差异在逐年减小，这是因为在出行量增加的情况下每辆出租车的平均碳排放减小了，符合实际情况。

在考虑空间位置碳排放量大小的情况下，通过K-means聚类进行碳排放空间演化分析。当聚类结果为6类时SSE值下降速度明显变缓，聚类效果较好.根据聚类结果将各个区域划分成3类：即低碳排放区域（1类与2类区域）、中碳排放区域（3类与4类区域）、高碳排放区域（5类与6类区域）.聚类后的碳排放区域分类效果如图7所示。

由图7可知，2016～2022年低碳排放区域总面积没有明显的变化，形状也始终保持环状，但低碳排放区域中的2类区域在逐年缩小并转化为1类区域，此现象在图中西北位置（未央区西部）表现比较明显。造成这一现象的原因是电动汽车的发展使车辆平均碳排放降低进而使区域平均碳排放有所降低。

中碳排放区域出现了分布上的明显变化，2016～2018年中碳排放区域面积缩小了约1/4，表现为整体向内收缩的趋势，2020～2022年中碳排放区域面积没有明显的缩减，但存在4类区域逐渐转化为3类区域并向市中心位置靠拢，部分区域碳排放聚集在路网上的现象。

与其他区域相比，高碳排放区域变化最为明显，面积减少52%，2016～2018年地图北部（未央区）高碳排放区域开始明显减小，此时其分布仍聚集在一起保持环状；2020年和2022年6类区域面积减少了70%以上，高碳排放区域不再聚成一团而是逐渐分配到各主干路上，从环状演化为网状。2022年碳排放聚集的各大路段主要有长安南路、西大街、北大街、南大街、高新路和西影路等，为西安出行量最大且经常出现拥堵状况的几条道路；几个碳排放突出的地点有钟楼、小寨、西安站、高新区（钟楼和小寨位于市中心位置并且是热门旅游地点与商圈聚集地，出行量巨大；西安站作为主要的交通枢纽全天交通需求量都比较大；高新区聚集着大量企业，通勤出行多），其碳排放维持在较高水平。

对各区域栅格中碳排放最大值、最小值和均值进行统计见表3。

从表3可以看出，各碳排放区域的碳排放最大值与均值都在逐年减小，1类区域的碳排放均值下降45%，2类区域下降39%，3类区域下降38%，4类区域下降40%，5类区域下降38%，6类区域下降27%，表明随着电动出租车的发展各区域平均碳排放均明显降低。

总体而言，在空间演化方面，碳排放区域随着电动出租车的不断增加由最开始的全区域环状逐渐变为外围低碳排放区域环状，中碳排放区域部分环状部分网状，高碳排放区域网状，且各碳排放区域均值都有明显的下降，整体趋势为低碳排放区域不断扩大，高碳排放区域逐渐向各出行量巨大的主干道以及人流量大的热点区域聚集。

将各年碳排放作为因变量，建成环境作为自变量运用GBDT模型求取出租车碳排放与各建成环境指标之间的相对重要度。

2016～2022年建成环境与碳排放之间的相对重要度关系如表4所示，从中可知，影响各年碳排放的建成环境变量重要度存在显著差异，其中居住用地与人口密度在各年对于碳排放重要度较高，平均贡献度达到了27.8%与27.55%，说明人口对于碳排放的增长有显著影响，这与Tan等［18］的研究结论一致。有几个因素促成了这一结果的合理性，首先，更高的人口密度意味着该地区有更多的交通需求，从而导致碳排量的增加.其次，居住用地密集的城市区域常伴有交通拥堵，在降低出租车运行效率的同时增加了碳排放.此外，居住用地的重要度表现为逐年降低的趋势，这主要由电动出租车的发展造成人均碳排放量的降低导致的。

各年的公共设施用地在工作日的重要度均大于周末，并且在2020年重要度较低，与其他学者的研究结果类似［46］，进一步说明日常活动（如上班和上学），会导致碳排放的时间异质性。特别是学校等教育公共设施，工作日其所引发的学生接送活动增加了出行需求，促进了碳排放。2020年受疫情影响，许多学校转为线上教学导致这部分出行量大大降低，导致该用地对碳排放重要度随之降低。

与公共设施用地相反，休闲娱乐用地的重要度在周末大于工作日，表明人们更加偏好在周末进行休闲娱乐，同时因为休闲娱乐用地中包含旅游用地，西安市作为旅游城市随着疫情结束出行量大幅度增加，并且游客在没有疫情时选择出租车的倾向会有所加强［47］，以上导致在2022年休闲娱乐用地重要度出现巨大提升。

人口密度与居住用地在2022年的相对重要度出现明显降低，通过以上分析可以发现，产生这一现象的主要原因是在该年公共设施与休闲娱乐用地重要度有较大的提升，同时人均碳排放的降低也对其有一定的影响。

剩余因素中，工业用地、公共交通站点与商业用地对于出租车的碳排放的重要度在各年均较低，说明相较于出租车在通勤出行时人们更倾向于选择公共交通或其他交通方式，这与人们以通勤为目的出行时的交通方式选择一致［48］。与吕莹等［40］在成都的研究结果不同，后者表明公共交通站点对碳排放有显著影响，这种差异可能归因于西安和成都两个城市不同的结构特征。

（1）各年周末碳排放均大于工作日碳排放量，且周末与工作日碳排放差异随着电动出租车比例的增加在逐年减小。

（2）随着时间推移总体碳排放演化趋势呈现为不断下降，到2022年时周末总体碳排放下降了约56%，工作日总体碳排放下降了约40%。

（3）各碳排放区域的碳排放最大值与均值都在逐年减小，中碳排放区域在不断向低碳排放区域转化，高碳排放区域向各大主干道转变；且随着电动出租车的发展西安市的碳排放区域由最开始的全区域环状逐渐变为外围低碳排放区域环状，中碳排放区域部分环状部分网状，高碳排放区域网状。

（4）居住用地与人口密度对于碳排放重要度在各年均较高，公共设施用地在工作日重要度均大于周末，休闲娱乐用地周末重要度均大于工作日，工业用地、公共交通站点与商业用地的重要度在各年均较低。

（5）基于以上研究，道路交通减碳应重点做好以下工作：大力发展电动汽车，增加电动汽车的使用比例，提高汽车的能量利用效率；从各大主干道入手降低交通高碳排放区域的碳排放是城市道路交通减碳的有效途径，通过有效的交通管理和控制来减少交通拥堵从而减小高碳排放区域的碳排；本研究对电动汽车碳排放计算以火力发电为基础，大力发展太阳能、水力、风力等各种清洁能源发电可以进一步减少碳排放。