随着中国"双碳"目标的提出,以及政府对交通充电、补能行业的大量政策支持,中国已经成为全球新能源车辆推广和充换电基础设施建设的领先国家。中国2022年电动汽车销量占全球销量的60%,占国内所有新乘用车注册量的四分之一,成功实现了"十四五"节能减排工作计划中20%的电动汽车份额目标。
在出租车电动化领域,《北京市"十四五"时期交通发展建设规划》提出"十四五"时期巡游出租车100%新能源化的政策目标。截至2022年底,中国电动汽车公用充电桩保有量达到100万根,占全球充电桩总量的51%。
然而,由于充电基础设施建设初期缺乏电动车实际运营数据的支撑,部分充电站点在选址和容量规划上不尽合理。加之出租车充电负荷高、时空分布随机性强的特点,其充电供需时空错配现象明显,表现为充电热点地区"一桩难求",而其他地区利用率低下,造成资源浪费。
因此,仅依赖现有燃油车数据对其潜在充电需求进行精确估计,并挖掘其充电需求的时空耦合关系,可以为交通电动化转型起步地区的政府部门、充电运营企业和电力部门进行决策、规划、投资和调度提供关键支撑。
本文提出的出租车潜在充电需求估计及时空结构分析方法的技术路线包括:
图1 出租车潜在充电需求估计及时空特征分解技术路线
该概念模型以出租车是否停留为标准将轨迹划分为驾驶区段和停留区段。停留状态的判别需提前指定阈值δ₁=5分钟。当停留时间tⱼ>δ₁时该区段被识别为restⱼ,restⱼ和restⱼ₊₁间的区段被识别为driveⱼ。
根据已有文献对出租车司机充电行为、电动汽车电量消耗的研究,对充电需求识别与计算做出以下假设:
图2 出租车潜在充电需求估计概念模型示意图
算法将某车辆连续行驶轨迹初始化为根节点,判断当前轨迹长度是否小于一次充电行驶阈值threshold。使用基于C++语言开发的FMM在Ubuntu环境下对轨迹数据和北京市道路网进行匹配计算,将轨迹点匹配到路段后,累计经过的路段长度获得真实轨迹长度trajᵢ。
在获取到精确的轨迹长度trajᵢ后,若轨迹长度trajᵢ在续航里程threshold之内则无需充电,否则依据假设寻找车辆最大停留区段作为充电区段。
图3 基于出租车轨迹地图匹配的充电需求二叉树生成算法示意图
从2方面验证潜在充电需求估计的时间准确性:① 电量-时间曲线,用于反应区域总体的充电电量波动;② 定义各区段时间估计匹配率A,考察真实充电区段被其估计时间区段所覆盖的比率。
验证潜在充电需求估计的空间准确性,首先需要将充电地点空间位置的真实值和估计值进行集计,以网格作为空间分析单元进行误差度量。设计了Pₛₗₒₜ(E≤ε)和Pₚₒwₑᵣ(E≤ε)指标来度量出租车辆潜在充电需求的空间分布估计误差。
以北京市六环内道路网为空间对象,采用K-邻近算法将充电地点经纬度匹配至道路网,并将该点充电需求电量、充电区段起止时间赋值给道路网路段属性。
构建了出租车潜在充电需求时空分布矩阵Mₘ×ₙ对信息进行集计表示,矩阵的行数m表示一天内的24 h,列数n表示充电需求匹配到路段后有充电行为发生的路段数目。应用SVD将潜在充电需求时空矩阵分解为3个矩阵的乘积。
本文使用的数据包括北京市传统出租车轨迹数据、北京市小规模电动出租车轨迹数据、以及北京市道路网数据:
本文以北京市六环内道路网路段为空间分析的基本单元,标注出11个重点考察区域,包括回龙观、天通苑等超大型社区;望京、国贸等商业区域;北京城市副中心、亦庄等规划发展区域,以及北京奥运村、首钢工业园等文旅重点区域和北京首都国际机场。
图4 北京市六环内主要道路网络及部分重点区域
模型参数 | 充电区段数量/个 | 充电需求电量/(% SOC) | 平均每个区段需求电量/(% SOC/个) | 充电总时长/h |
---|---|---|---|---|
真实值 | 4 022 | 202 942 | 50.46 | 34 085 |
θ=50% | 7 822 [94.48%] | 205 761 [1.39%] | 26.31 [47.87%] | 33 986* [0.29%] |
θ=100% | 3 276 [18.55%] | 205 177* [1.10%] | 62.63 [24.12%] | 29 734 [12.77%] |
注:[ ]中的数字表示各估计指标同真实值间的相对误差ε,*表示各指标表现最好的估计值
不同θ参数使得本文所提模型和算法对提取不同指标产生不同适应性。较小θ参数值能够捕捉短时间充电需求,从而提高充电总时长估计准确率和各区段时间估计匹配率A;而较大θ参数值能够提升对充电区段数量、充电总能量等指标的准确率。
图5 不同空间尺度和θ值下的充电区段数量和充电电量估计误差累积分布
随着空间尺度增加,空间栅格对估计量的容错能力增大。当使用参数θ=90%时的模型估计充电区段数量的空间分布时,45.92%的500 m栅格误差小于0.2,而10 000 m栅格误差小于0.2的占比达到78.79%。此外,不同θ参数对充电需求估计量误差的影响较大,参数θ值越高,模型估计的空间分布误差越小。
选取中型车作为现阶段北京市出租车的代表车型,其百公里能耗δ和续驶里程分别为12.6 kWh/100 km和468 km。假设电动出租车在电池电量消耗θ为90%SOC时产生充电需求,即一次充电行驶阈值threshold为420 km。
应用基于轨迹地图匹配的潜在充电需求二叉树生成算法从2019年3月9日—11日连续3d的出租车轨迹数据中提取出北京市六环范围内发生充电请求的出租车1 913辆,有效充电区段4009段,充电总量达107 434 kWh,日均35 811 kWh。
图6 城市道路网尺度下的典型工作日北京市六环内出租车充电需求空间分布
图7 典型日北京市六环内出租车充电需求时空特征奇异值分解结果
通过奇异值分解对2019年3月10日北京市六环内出租车充电需求时空分布一阶矩阵进行分解降维,识别充电需求中存在的时空结构。前3项奇异值便能表示69%原矩阵的时空信息:
本文提出的充电需求估计模型和算法是一种基于燃油车轨迹大数据的自下而上的方法框架。所提的概念模型通过3种替代关系阐述了将以"驾驶-停驶"为循环的轨迹转换为以"充电-放电"为循环的轨迹的技术方法;而二叉树算法在Ubuntu环境下通过递归遍历对上述概念模型进行实现,将叶子节点输出为充电需求,包括充电电量、充电时间和充电空间位置分布。
本文使用890辆带有电池状态字段的电动出租车轨迹数据对模型和算法进行验证,估计模型的实验结果表明,其在不同电池电量消耗阈值θ和空间尺度下具有较高的时空分布估计精度。
以北京市六环内区域为实证案例,利用连续3d部分出租车轨迹数据,分析北京市出租车潜在充电需求分布的地理空间格局,研究发现出租车潜在充电需求量呈现明显的分区域、分通道的空间聚集特征。基于奇异值分解对潜在充电需求的时空结构进行分析,得到3种类型的充电需求:常态化充电需求、上下午异质性充电需求以及工作时段和非工作时段异质性充电需求。
随着出租车队电动化的推进、电池技术的进步以及居民活动强度的变化,出租车充电需求与供给的互动关系不断演化。相比于研究已经满足的充电需求,挖掘潜在的充电需求及其时空变化模式有助于解析出租车主真正的充电需求响应规律。
值得注意的是,使用本文所提模型和算法来估计出租车潜在充电需求的结果受电池容量、电量消耗速率和驾驶行为异质性影响,充电需求时空分布存在不确定性。此外,本文所提方法在充电需求估计精度、求解算法效率和模型迭代发展等方面仍有改进空间。
随着补能方式的多样化发展,还有待探索换电行为和驾驶轨迹特征的关联关系,不断迭代该模型用来协同估计充电需求和换电需求。最后,开发计算高效、强表征能力的充电需求时空结构挖掘方法也是未来研究方向之一。
政府规划部门以潜在充电需求估计为支撑制定出台近中长期的电动汽车充电基础设施规划,提供分批次、分阶段的需求基础数据,避免铺张建设导致的资源浪费。
充电桩服务建设、运营商等要以充电需求时空分布为参考有效部署充电基础设施,避免均质化与随机化,而是要以提升充电桩利用率为导向优化选址流程和定容策略。
电力部门要对电网负荷实施调度调节,提升电网的安全性、经济性。电力部门预计要为北京市待电动化的3万辆电车入网新增日均561 MWh的充电供给。
电力部门和交通运行管理部门要以高精度的充电需求估计为有力工具研究充电需求管理策略,包括被动管理策略如峰谷电价调节机制和主动管理策略如充电服务导航系统。