电动汽车能耗模型的校准:一种结合上下文的转移贝叶斯优化方法
《TRANSPORTATION RESEARCH PART E-LOGISTICS AND TRANSPORTATION REVIEW》:Calibration of electric vehicle energy consumption models: A context-Integrated transfer bayesian optimization approach
【字体:
大
中
小
】
时间:2026年09月02日
来源:TRANSPORTATION RESEARCH PART E-LOGISTICS AND TRANSPORTATION REVIEW 9.3
编辑推荐:
• 建立了一个具有多档动力系统的高保真电动汽车能量模型。
• 一个CBO框架能够捕捉不同条件下的上下文感知参数变化。
• TBO实现了跨不同车辆平台的快速知识迁移。
• 该框架使校准收敛速度提高了五倍。
• 实际数据验证了旅行级能量估计误差小于10%。
1.
• 建立了一个具有多档动力系统的高保真电动汽车能量模型。
• 一个CBO框架能够捕捉不同条件下的上下文感知参数变化。
• TBO实现了跨不同车辆平台的快速知识迁移。
• 该框架使校准收敛速度提高了五倍。
• 实际数据验证了旅行级能量估计误差小于10%。
1. 引言
1.1. 背景与动机
随着全球交通向电气化转型,电动汽车(EVs)已成为减少碳排放和推进可持续出行的关键技术(Chen和Liu,2023年;Liu等人,2024年;Wang等人,2024年)。在这一背景下,准确的能量消耗估计成为现代智能交通系统的核心要求,包括生态路线规划、充电计划、车队运营和网络级能源管理(Li等人,2026年;Liu等人,2023年;Liu等人,2022b年;Zhang等人,2020年)。因此,能够准确计算各种运行条件下电动汽车能量消耗的高保真能量模型对于解决这一问题至关重要(Liu等人,2025年)。准确的能量消耗估计可以缓解驾驶员的续航焦虑,提高路线规划可靠性,并支持电动汽车车队与能源基础设施之间的协调互动(He和Zhen,2025年;Liu等人,2022a年)。当前基于物理的能量消耗模型可以从纵向车辆动力学中得出,在不同的驾驶循环下提供可解释且可外推的预测(Koch等人,2021年;Luin等人,2019年;Xie等人,2020年)。这些模型的准确性在很大程度上取决于关键参数的正确识别,如车辆质量、滚动阻力、传动系统效率和支持功率。参数校准是从观测到的车辆运行数据中推断参数值的逆过程,因此对能量消耗模型的有效性至关重要。已经付出了大量努力来使用离线或低频校准策略来匹配观测数据,这有效地解决了不变参数中的系统偏差(Fiori等人,2016年;Miri等人,2021年)。
尽管取得了这些进展,但在动态和异构环境中实际应用高保真能量建模仍面临三个关键挑战。首先,校准精度的上限往往受到底层物理模型简化结构的限制。许多现有公式依赖于理想化的假设,例如恒定的再生制动效率或恒定的支持功率负载(Koch等人,2021年)。这种结构简化引入了固有的建模偏差,仅通过参数优化无法完全补偿,可能导致在极端运行条件下参数估计不准确或拟合效果不佳。其次,实时适应性对于实际能源管理至关重要。现代电动汽车配备了能够提供速度、充电状态和功率即时测量的车载传感器。但这些原始测量数据本质上是反应性的,而不是预测性的。相比之下,在线能源感知路线规划和规划等应用需要一个校准后的前向模型,以便模拟未来驾驶段的能量消耗(Fotouhi等人,2020年;Wang等人,2020年;Zhang等人,2020年;Zhen等人,2025年)。因此,模型应在早期阶段识别或精细化这些潜在参数(Kim等人,2023年)。这产生了对一个在线高效校准范式的实际需求,该范式在物理上可解释且在操作上具有适应性。第三,在多样化的运行条件和车辆平台上实际应用能量模型存在可扩展性障碍。一方面,电动汽车作为一个开放系统与其外部环境互动,其效率受驾驶行为、交通条件、环境条件和道路特征的复杂相互作用影响(Hao等人,2020年)。例如,空气动力学损失和滚动阻力都受到外部因素(如空气密度、温度和道路表面状况)的影响,这些因素在每次行驶中可能会有很大变化(Liu等人,2018年)。然而,现有的校准策略通常假设所有运行条件下的物理参数都是固定的,这限制了在不同行驶剖面和环境条件下的泛化能力(Fiori等人,2016年)。这种环境非稳定性激发了需要根据当前上下文条件快速调整内部模型参数的校准机制的需求。另一方面,即使为特定来源车辆使用大量实际数据仔细校准了高质量的能量模型,将其直接转移到新的目标车辆上也具有挑战性,因为车辆配置、电机效率和支持系统存在差异(Fukushima等人,2018年;Pan等人,2023年)。实际上,为每一种新引入的车辆模型收集精确数据既耗时又成本高昂。为了减轻这种数据依赖性,迁移学习作为一种强大的范式应运而生,它被定义为通过利用从一个或多个相关来源任务中预先获得的知识来提高学习者在目标任务上的性能(Pan和Yang,2010年;Zhuang等人,2021年)。尽管迁移学习在许多问题和案例中显示出潜力,但对于跨异构车辆平台的可转移校准策略的关注相对较少(Rahman等人,2026年)。
为了解决这些挑战,本研究旨在开发一个统一的两阶段校准框架,平衡物理基础的建模与数据驱动的灵活性。具体而言,我们提出了一种集成方法,首先通过上下文贝叶斯优化(CBO)学习环境条件与潜在物理参数之间的关系,然后通过迁移贝叶斯优化(TBO)实现跨车辆的知识迁移。通过结合上下文感知的参数适应和可转移的先验知识,所提出的框架旨在在不同驾驶条件和车辆平台上保持高能量估计的准确性和效率。
本研究希望为电动汽车行业提供一种在线校准工具。通过将能量模型与我们的算法集成,可以在最少数据需求的情况下实现不同类型车辆和驾驶条件下的准确能量估计。这种范式的转变可以大幅减少对劳动密集型现场实验和工程师手动参数调整的依赖,从而降低大规模电动汽车车队的成本和部署周期。
1.2. 研究贡献
本文的主要贡献总结如下:
• 改进的高保真物理模型:我们建立了一个包含复杂车辆动力系统的结构物理架构,包括多档变速比、动态再生制动效率和速度依赖的支持功率消耗。这种全面的公式减少了结构建模误差,同时使校准后的参数表现出与物理预期大体一致的趋势。
• 跨条件在线校准:我们开发了一种CBO机制,该机制具有基于上下文局部信息的专门获取功能。这使得系统能够通过利用区域环境相似性进行“零样本”或“少样本”参数适应。通过在搜索过程中显式建模非静态的上下文效应,所提出的方法防止了搜索停滞,并确保在不同未见过的驾驶行程中实现稳健的能量估计。
• 跨车辆知识迁移:我们提出了一种新的迁移学习范式,该范式基于来自源域的优化CBO映射构建目标车辆的替代模型。这种方法将学到的参数-性能相关性迁移到新的车辆平台上,有效地将高级物理启发式注入目标优化过程,降低了异构电动汽车高精度建模所需的数据阈值。
本文的其余部分组织如下:第2节回顾相关文献。第3节介绍基于物理的能量模型。第4节和第5节分别详细描述了所提出的CBO和TBO算法。第6节提供数值实验和分析,最后第7节给出结论。
2. 相关工作
过去十年中,电动汽车能量消耗模型的研究方向发生了变化。现有研究在如何平衡物理可解释性、数据依赖性和计算效率方面有所不同。总体而言,当前方法可以大致分为基于物理的公式、数据驱动的方法以及两者的组合,其中参数校准在连接理论与实际应用中起着核心作用。
2.1. 能量消耗建模方法
大部分文献基于车辆动力学的物理建模。在这些研究中,能量消耗通常是通过纵向力分析得出的,其中牵引力需求通过考虑空气动力学阻力、滚动阻力和道路坡度来计算(Fiori等人,2016年)。这一类别还包括将车辆动力学集成到微观交通模拟器中的基于模拟的方法(Koch等人,2021年;Luin等人,2019年;Xie等人,2020年)。这些模型在物理上是可解释的,但其性能在很大程度上取决于参数的准确性。当外部条件(如交通、温度或道路属性)发生变化时,固定的参数集往往无法保持令人满意的准确性(Liu等人,2018年)。
随着实际驾驶轨迹数据的日益可用,数据驱动的方法引起了越来越多的关注。这些方法包括基于现场数据训练的回归模型(Chen等人,2021年;Wu等人,2023年)以及先进的机器学习技术,后者在捕捉复杂非线性关系方面具有灵活性。先前的研究探索了多种技术,从集成学习方法(Ullah等人,2021年)到深度神经网络(Adedeji和Kabir,2023年;Pan等人,2023年;Topi?等人,2019年)。尽管这些模型通常能实现强大的预测性能,但它们通常需要大量数据集,并且比基于物理的模型更难以解释。此外,当应用于训练领域之外时,它们的预测准确性往往会显著下降。
鉴于基于物理的方法和数据驱动方法的互补优势,混合建模已成为一种主流趋势。现有的混合方法通常可以分为三种类型。第一种类型结合了基于物理的模型和机器学习模型的输出,通常通过加权融合或集成策略,确保预测受益于物理约束和数据驱动的灵活性(Pan等人,2023年;Ullah等人,2021年)。第二种类型使用基于物理的模型作为基线,并应用机器学习方法来估计由简化物理假设引起的残差误差(Fiori等人,2016年;Hjelkrem等人,2021年)。第三种类型通过将控制关系引入训练目标,将物理知识直接嵌入学习过程中,通常称为基于物理的学习(Karniadakis等人,2021年;Raissi等人,2019年)。与纯数据驱动的预测模型相比,这种混合策略保留了更强的可解释性,因为底层模型仍然是车辆动力学的确定性表示。然而,为了在不同条件下实现稳健的泛化,确定这些模型的最佳物理参数仍然是一个关键前提,这推动了先进参数校准技术的需求。
2.2. 参数校准和可转移性
参数校准被定义为从观测到的车辆运行数据中识别物理常数的逆过程,以确保能量消耗模型能够准确再现现实世界的行为(Rakha等人,2004年)。传统上,这被形式化为一个优化问题,其中使用历史观测值通过最小二乘拟合或贝叶斯估计来最小化残差(Frazier,2018年;Wang等人,2020年)。然而,这些数值策略通常计算效率低且缺乏可扩展性,通常在引入新的车辆平台或运行条件时需要重新运行(Kim等人,2023年)。为了缓解这一问题,特别是贝叶斯优化(BO)等辅助方法由于其能够通过平衡探索和利用在高维参数空间中导航而受到关注(Huo等人,2023年)。
尽管标准BO效率较高,但由于它通常假设参数和输出之间的静态映射,其在动态环境中的性能往往下降。最近的进展通过引入上下文贝叶斯优化(CBO)来解决了这一问题,该算法将环境或操作变量(上下文)集成到获取函数中。例如,Ayala-Romero等人(2021年)展示了贝叶斯在线学习在资源协调方面的有效性,突出了其在实时能源感知适应中的潜力。为了解决高维上下文的挑战,Martinelli等人(2023年)提出了识别哪些上下文变量真正影响目标函数的方法。这些研究表明,集成上下文对于实现跨不同驾驶条件的在线适用性至关重要。
虽然CBO处理了域内变异性,但在不同域之间的适应问题仍然是一个瓶颈。迁移学习旨在通过重用从相关源域学习的表示或先验来提高目标任务性能(Pan和Yang,2010年;Zhuang等人,2021年)。尽管迁移学习在深度强化学习和一般预测任务中取得了成功(Tan等人,2021年;Zhu等人,2023年),但其在基于优化的校准中的应用仍在发展中。Yang等人(2025年)回顾了智慧城市和可持续交通的跨城市迁移学习,强调了在异构域之间转移知识的机会和挑战。最近的研究将边界推向了极端样本效率;例如,Meindl等人(2025年)引入了ZeroShotOpt的概念,利用预训练模型以零样本方式实现高效的黑盒优化。然而,一种能够将情境感知与跨车辆可转移性协同结合在一起,以实现高保真能量建模的系统架构尚未完全建立。2.3. 讨论与研究空白总结来说,现有文献明显展示了从基于硬性物理模型的发展,向灵活的混合模型和辅助校准框架的转变。然而,尽管取得了这些进展,对于现实世界中的大规模电动汽车能源管理来说,仍有两个关键的研究空白尚未得到解决。一方面,大多数现有的校准方法将参数识别视为一个静态的、离线的过程。它们无法捕捉到由于不同运行情境导致的参数瞬态变化,这限制了它们的在线适应能力。另一方面,当前的基于优化的校准框架在很大程度上是特定于车辆的。当在新的车辆平台上部署时,校准过程必须完全从头开始,忽略了不同电动汽车之间共享的物理和结构相似性。为了弥合这些空白,迫切需要一个统一的框架,该框架不仅能够适应动态驾驶情境,还能够跨不同车辆领域传递学到的先验知识。受此启发,我们的研究提出了一种结合情境贝叶斯优化和转移贝叶斯优化的集成架构,旨在实现复杂交叉条件和跨车辆场景下的高效、少样本参数校准。3. 基于物理的能量模型在本节中,我们提出了用于估算车辆能耗的物理模型。我们基于Koch等人(2021年)的研究,定制了一种自下而上的建模方法,整合了纵向动力学、传动损失、电机效率图和电池内阻,以确保在不同运行条件和车辆类型下的适应性。3.1. 模型概述所提出模型的计算过程如图1所示。通过考虑每个仿真时间步骤中的瞬时车辆运行状态(包括速度v、加速度a和道路坡度α),该模型通过考虑空气阻力、滚动阻力、坡度阻力和加速度需求来计算车轮所需的驱动力。随后,该模型通过传动效率映射、电机功率损失图和电池内阻模型来估算电机侧功率和电池化学功率。最后,总能耗通过时间步长累积得到。下载:下载高分辨率图像(671KB)下载:下载全尺寸图像图1. 所提出的物理模型概述。3.2. 驱动力模型和纵向动力学车轮所需的驱动力FW由四个主要阻力组件和加速度需求的总和决定。根据车辆纵向动力学的原理,FW表示为:(1)????=????????????????+?????????????+?????????????+??????????????????????各个组件的定义如下:(1) 坡度阻力:当车辆在具有纵向坡度的道路上行驶时,沿道路方向的重力分量构成了坡度阻力:(2)????????????????=??????sin???其中道路坡度α来自高程剖面。(2) 滚动阻力:这来源于轮胎与路面接触区域的弹性变形:(3)?????????????=???????????cos???其中fr是反映轮胎和路面条件的滚动阻力系数。(3) 空气阻力:空气阻力是能耗的一个相对重要的来源,尤其是在高速行驶时:(4)?????????????=12????????????????2其中Cw是表征车辆空气动力形状的空气阻力系数。(4) 惯性力:为了考虑车轮和电机转子等组件的平移惯性和旋转惯性,引入了一个等效惯性因子ei:(5)??????????????????????=??????????(6)????=1+Θ?????2??其中Θ是转换为轮轴的等效旋转惯性,rw是有效滚动半径。然后推导出所需的轮端扭矩TW:(7)????=????·????正数值(TW?>?0)表示电机必须输出正扭矩的驾驶条件。负数值(TW?0)表示制动或滑行,在这种情况下,可以通过再生制动回收部分能量。3.3. 传动系统模型传动系统配备了一个减速器,其变速比iT(可以是固定值或从齿轮比集合中选择),将高速电机输出转换为车轮所需的低速高扭矩。对于具有nm个相同电机的车辆配置,电机速度nEM和每个电机所需的扭矩TEM,s计算如下:(8)???????=??·????2????????·60(9)???????=2???????????60(10)???????,??=1????·{????????·????,????≥0????·????????,????<0其中ηT是传动效率。单个电机的扭矩TEM,s还受到电机物理能力的限制:(11)???????????,????????≤???????,??≤????????(12)???????????,????????≤???????,??·???????≤????????其中Tmax, Pmax和Treg,max, Preg,max分别代表每个电机的推进和再生峰值限制。3.4. 电机效率和损失模型电机系统PEM的总电能需求el考虑了所有电机的机械输出、它们的内部损失以及再生制动效率。总电能表示为:(13)???????,??????????=????·(???????,??·???????)(14)???????,???????????=????·???(???????,???????,??)(15)???????,?????={???????,??????????+???????,???????????,???????,??????????≥0???????,??????????·????????????????+???????,???????????,???????,??????????<0其中ηregen是再生效率系数,表示可以有效地回收到直流母线的电机侧负功率的比例;电机损失PEM,loss是速度nEM和扭矩TEM的二元函数,由损失图表征:(16)???????,???????????=???(???????,???????)损失图作为一个二维查找表存储,水平轴表示电机速度,垂直轴表示电机扭矩。每个网格点代表该特定运行点的可测量功率损失。根据此图,电机效率ηEM定义为???????=???????,??????????/???????,?????。如图2所示,该图展示了宝马i3电机的效率特性,系统在中高扭矩和速度区域表现出高效,而在低负载条件下效率降低。下载:下载高分辨率图像(424KB)下载:下载全尺寸图像图2. 所提出的基于物理模型中使用的宝马i3电机效率图。3.5. 电池和辅助功率模型电池PBat的输出总电能包括电机需求和辅助负载:(17)??????????,?????=???????,?????+??????????在这个模型中,辅助功率消耗Paux不再是一个简单的常数。相反,它考虑了一个与速度相关的残差项和气候控制负载:(18)??????????=????????????????+????????????????_????????·|??|+???????+????????????其中Pconst是基础恒定负载,????????????????_????????是速度相关损失的系数,PAC/Pheat代表空调和制热功率。为了估算内部能量损失,我们采用了内阻模型(U0, Ri)。终端电压U和PBat之间的关系为??????????,?????=?????=(??0????????)???。求解电池化学功率PBat,chem得到:(19)??????????,??????????=??0???=??202????????02??????√??20?4????????????????,?????3.6. 能量消耗通过积分瞬时化学功率,可以得到整个行程期间的总能耗E:(20)??=∑????????????,???????????(??)·Δ???这种自下而上的逐步计算框架确保了模型准确反映了能量转换链,为参数校准和优化提供了可靠的物理基础。4. 情境贝叶斯优化在本节中,我们介绍了为车辆能量模型的自适应校准设计的情境贝叶斯优化(CBO)框架。如第3章所示,物理能量模型中的几个关键参数并非车辆硬件的固有常数,而是对外部运行条件(情境)非常敏感。这些参数包括滚动阻力系数fr、空气阻力系数Cd、电池内阻Ri、恒定辅助功率Pconst、速度相关的功率系数kv和再生效率ηregen。公式19中的开路电压U0被视为一个常量参数,固定在标称电池电压。这是一个合理的近似,因为在校准行程期间的SOC变化相对较小,内阻Ri被校准以吸收二阶SOC依赖的效应。这些参数的最优值根据不同的行程特性(如环境温度、天气条件和道路类型)而显著变化。传统的校准方法通常在遇到新环境时需要重新进行彻底的优化,这对于实时应用来说是计算上不可行的。为了使能量模型能够快速适应未见过的情况,我们提出了基于CBO的方法,如图3所示。通过将环境情境作为额外的输入维度,CBO使模型能够跨不同行程进行知识转移,有效地学习从运行条件到最优参数配置的通用映射。以下小节详细介绍了我们提出方法的数学公式和实现。下载:下载高分辨率图像(783KB)下载:下载全尺寸图像图3. 所提出的CBO模型概述。4.1. 基础贝叶斯优化贝叶斯优化是一种强大的全局优化工具,适用于计算密集型目标函数,特别是当函数的解析形式不可用或其梯度难以计算时。在车辆能量模型校准的背景下,目标函数f(x)将一组物理参数??∈??映射到模拟能耗与实际测量值之间的差异。由于高保真动力总成仿真的复杂性,BO通过构建概率替代模型来指导搜索过程,从而成为一种样本效率高的策略。4.1.1. 使用高斯过程的替代建模为了近似昂贵的物理仿真,采用了高斯过程(GP)作为替代模型。GP是一组随机变量的集合,任意有限数量的这些变量具有联合高斯分布。它完全由均值函数m(x)和协方差函数(核)κ(x, x′)指定,表示为:(21)???(??)~?????(???(??),???(??,??′))其中核函数κ(x, x′)编码了参数空间的平滑度和相关结构。给定一组t个观测值????={(????,????)}????=1,函数在新点x处的后验分布仍然是高斯的,提供了预测均值μt(x)和预测方差??2???(??):(22)?????(??)=???????1???1:??(23)??2???(??)=???(??,??)????????1???其中k是新点与先前观测值之间的协方差向量,K是观测数据的协方差矩阵。这种不确定性量化使算法能够区分高置信区域和未探索区域。4.1.2. 获取函数获取函数α(x)利用后验均值和方差来确定评估候选参数集的效用。它有效地平衡了利用(查询预测误差低的区域)和探索(查询不确定性高的区域)之间的权衡。在本研究中,我们使用期望改进(EI)标准,它量化了相对于当前最佳观测值??*=min?{??1,?,????}的预期收益。对于最小化问题,EI定义为:(24)????????(??)=???[max?(0,??*????(??))]通过最大化αEI(x),优化器选择下一个最有可能显著提高校准精度的点。4.1.3. 迭代优化过程参数校准的BO迭代循环总结如下:1. 初始化:使用拉丁超立方抽样(LHS)生成N个初始点,形成起始数据集??0={(????,????)}????=1,通过执行物理仿真。2. 替代模型拟合:将GP模型拟合到当前数据集??,并通过最大似然估计(MLE)优化核超参数(例如,长度尺度)。3. 候选建议:解决内部优化问题,找到最大化获取函数的点:(25)????+1=arg?max??∈??????(??∣????)4. 仿真评估:在????+1处评估昂贵的物理模型,以获得实际观测值????+1。5. 迭代:更新数据集????+1=????∪{(????+1,????+1)},并重复步骤2–4,直到达到最大迭代次数或满足收敛条件。4.2. 适应性校准的情境公式标准BO将每个校准任务视为一个独立的优化问题。然而,本研究的关键观察结果是,最优物理参数取决于每次行程的环境和运行情境。因此,我们不是重复从头开始解决孤立的优化问题,而是将校准过程表述为一个情境优化问题,在其中情境描述为了确保在低负载区间或空闲情况下的数值稳定性,分母被下限定为1:(28)?WMAPE=∑????=1∣??????(??;??)??????????????(??)∣∑????=1max?(1,?????????????(??))•累积能量误差(?Total):总行程级别能量消耗的绝对百分比误差。•RMSE百分比(?RMSE):通过平均观测消耗量标准化的功率轮廓的均方根误差,用于捕捉时间上的不匹配。•物理有效性惩罚(?Penalty):一个约束项,用于量化物理模型无法计算出可行扭矩输出的时间步长百分比,确保建议的参数保持在物理上一致的范围内。4.2.2. 联合上下文–参数高斯过程替代模型为了实现具有不同运行条件的行程之间的转移,我们构建了一个在校准参数和上下文联合输入空间上的GP替代模型。具体来说,定义扩展输入向量(29)??=[???,???]?∈??=??×??,并将上下文目标建模为(30)???(??,??)=???(??)~?????(???(??),???(??,??′)),其中m(z)是先验均值,κ(z, z′)是协方差核。为了反映物理参数和上下文变量的不同作用,实现的核使用了两个子空间的单独自动相关性确定(ARD)长度尺度。等效地,该核可以表示为两个平方指数分量的乘积:(31)???(??,??′)=??2???exp??? ???12?????∑??=1(?????????′??)2?2??,???12?????∑??=1(?????????′??)2?2??,???? ??=?????(??,??′)·?????(??,??′),其中?x,j和?c,k分别是参数和上下文维度的维度长度尺度,??2??是信号方差。此外,为了数值稳定性,在格拉姆矩阵的对角线上添加了一个观测噪声项??2??:(32)????=[???(????,????)]????,??=1+??2?????给定累积的上下文数据集(33)????={(????,????,????)}????=1,在候选对(x, c)处的后验预测分布仍然是高斯的:(34)?????(??,??)=?????????1?????1:??(35)??2???(??,??)=???(??,??)??????????1???????其中kz是候选输入z与所有历史观测值之间的协方差向量。4.2.3. 超参数估计和定期重新调整而不是预先固定核的超参数,该框架定期从累积的上下文数据集中重新估计它们,以平衡适应质量和在线计算成本。设(36)??={???,1,…,???,????,???,1,…,???,????,??2??,??2??}收集所有GP超参数。这些是通过最大化对数边际似然来估计的:(37)log????(??1:??∣??1:??,??)=?12????1:??????1?????1:???12?log?|????|???2?log?(2???)其中??1:??={??1,…,????},Kt通过方程(32)依赖于θ。在实现中,超参数更新不是在每次迭代时都进行。相反,它们在初始化时触发,然后根据可配置的时间表在在线校准期间定期进行,这减少了重复优化的开销,同时保持了足够的适应性。为了避免当历史上下文数据集变大时出现过高的立方复杂性,超参数拟合步骤可以进一步在累积观测值的一个 capped-size 子样本上进行。4.2.4. 基于上下文的获取与局部期望改进在部署时,假设新行程的上下文c(new)在校准开始前是可观察或可推断的。因此,优化问题简化为在固定上下文下选择最佳物理参数向量:(38)????+1=arg?max??∈??????(??∣??(????????),????)其中获取是从基于上下文的后验中评估的:(39)???(??,??(????????))~???(?????(??,??(????????)),??2???(??,??(????????)))重要的是,为了避免由不相关的运行条件引起的过于乐观的基线,我们用上下文局部当前值替换了全局历史最佳值。设???(??(????????))表示与其上下文足够相似的历史观测值的邻域,可以通过基于半径的规则构建,或者当样本太少时,通过在标准化上下文空间中使用k最近邻来构建。然后定义上下文局部当前值为:(40)????(??(????????))=min??∈???(??????????))?????而不是全局最佳值min?iyi。因此,实现的EI为:(41)????????(??∣??(????????),????)=???[max?(0,????(??(????????))????(??,??(????????))]这种基于上下文的EI是标准CBO循环的关键实际改进。它保留了EI的探索-利用权衡,同时使获取基线更好地与新行程的目标运行状态对齐。注意,为了确保数值稳健性并避免昂贵的连续内部优化,获取函数是在通过LHS生成的离散候选池?????????????上评估的。4.3. 算法工作流程为了与实际实现保持一致,所提出的框架包括两个紧密耦合的阶段:(i) 一个离线上下文训练阶段,在处理历史行程时逐步训练联合上下文–参数GP替代模型;(ii) 一个在线上下文校准阶段,该阶段重新使用这个逐步预训练的GP对新遇到的行程进行零样本推荐或少量样本自适应细化。4.3.1. 离线上下文数据集构建离线阶段是一个完整的替代模型预训练阶段:每个新评估的历史样本立即追加到共享的上下文数据集中,并且在BO迭代过程中重复拟合/更新(x, c)上的联合GP。因此,离线阶段的输出是上下文数据集和已经为在线部署预训练的GP。完整的离线上下文数据集构建过程在算法1中总结。下载:下载高分辨率图像(285KB)下载:下载全尺寸图像算法1. 离线上下文数据集构建。在离线阶段之后,共享的上下文数据集定义为:(42)??????????????????????????????={(????,????,????)}????????????=1,相应的联合GP已经在这个数据集上逐步预训练完毕。4.3.2. 在线上下文校准对于新遇到的具有观测上下文c(new)的行程,在线阶段将上下文数据集初始化为:(43)??????????0=??????????????????????????????,并从预训练的离线模型状态初始化联合GP替代模型,然后进行可选的在线更新。根据可用的在线评估预算,支持两种操作模式。(1) 零样本上下文推荐。如果没有在线评估预算,框架直接使用预训练的联合GP后验对有限的候选池进行排序,并返回预测后验均值最小的参数向量:(44)?????????????=arg?min??∈?????????????0?(??,??(????????))其中??????????是通过LHS生成的。这种模式提供了即时的上下文感知参数推荐,无需任何特定于行程的评估。(2) 少量样本在线上下文校准。如果在线评估预算有限,框架在固定上下文c(new)下进行少量行程特定的更新。在每次迭代中,生成一个候选池,并通过最大化上下文局部EI获取函数选择下一个查询点。新评估的样本被追加到联合数据集中,然后从离线预训练的状态中启动GP,并根据更新计划定期重新估计。完整的在线上下文校准过程在算法2中总结。下载:下载高分辨率图像(243KB)下载:下载全尺寸图像算法2. 新行程的在线上下文校准。5. 迁移学习5.1. 问题设定我们刚刚建立了一个用于单一车辆平台内跨条件参数校准的上下文贝叶斯优化框架。然而,仅在同一车辆类型内的跨条件适应仍然不足以满足实际应用的可扩展性要求。当模型从源车辆转移到目标车辆时,源域中学到的知识不能直接在目标域中复用。为了解决这个限制,我们进一步研究了跨车辆转移校准问题。具体来说,在保持上下文建模原则的同时,它将源域知识引入目标域优化中,以提高在目标域样本有限的情况下的参数校准效率,如图4所示。下载:下载高分辨率图像(358KB)下载:下载全尺寸图像图4. 提出的TBO模型概述。为了形式化这个问题,设S表示源域,T表示目标域。对于任何给定的上下文c,让(45)?????(??,??),?????(??,??),??∈??,??∈??,分别表示源域和目标域中真实但未知的校准目标函数。因此,目标域中的校准目标可以写为(46)??*???(??)=arg?min??∈????????(??,??)。这表明本节的核心问题不再是源域内的上下文泛化,而是如何在目标域评估有限的预算下获得新的目标域样本的上下文特定最优解??*???(??)。需要强调的是,本研究考虑的迁移并不是直接复用在源域获得的最优参数,而是转移关于上下文、参数配置和校准损失之间关系的结构知识。为了实现这种跨域迁移,本研究进一步假设源域和目标域采用语义上一致的上下文表示。在此基础上,本节的核心任务是构建一个基于替代模型的优化机制,该机制可以继承源域的先验知识,并通过目标域的观测逐步进行纠正,以便在有限的评估预算下在目标域中高效搜索近似最优参数解。本质上,这里考虑的迁移问题是在信息丰富但不完美的跨车辆先验知识下的样本高效目标域校准问题。在提出的迁移公式中,源域知识通过训练后的上下文替代模型的后验均值传递,作为目标域搜索的固定源信息基线。源后验方差不会传播到目标替代模型;目标域的预测不确定性由下面引入的残差GP表征。单独的源信息基线是不够的,因为车辆特定的偏差仍然不可避免。因此,关键的方法学要求不仅是复用从源车辆学到的可迁移结构,还包括建立一种可以使用少量目标域观测逐步纠正这种先验偏差的机制。5.2. 基于源先验的迁移替代模型与残差校正所提出的方法将这种源域知识纳入目标域建模中,以便目标域优化可以从一个有信息的先验开始,而不是从完全无信息的搜索状态开始。根据第4节中的联合高斯过程公式,源域中的上下文目标函数写为(47)?????(??,??)=?????(??)~?????(?????(??),?????(??,??′)),??=[???,???]?。这里,fS(x, c)表示源域中真实但未知的校准目标函数,而mS(z)和κS(z, z′)分别是源域高斯过程的先验均值函数和协方差核。给定源域历史数据集????={(????,??????)}??????=1,任何输入(x, c)处源域目标的后验预测分布可以写为(48)?????(??,??)∣????~???(?????(??,??),??2???(??,??))。这里,μS(x, c)和??2???(??,??)分别表示源域替代模型的后验均值和后验方差。由于这个后验已经编码了校准参数的变化如何影响源域中不同上下文条件下的校准损失,μS(x, c)在提出的迁移框架中用作目标域建模的先验基线。然而,源域和目标域通常在车辆质量、动力系统特性、电机效率和电池属性以及辅助系统配置方面存在差异。因此,源后验均值不能完全代表目标域中的真实响应表面。如果直接将μS(x, c)作为目标域响应的近似值,那么在跨车辆迁移中常见的系统差异将无法捕捉到。为了解决这个问题,引入了一个目标域残差项,目标域目标模型表示为(49)?????(??,??)=?????(??,??)+?????(??,??),其中δT(x, c)表示相对于源域基线的目标域残差过程。这种分解意味着源后验均值提供了整体结构先验,而目标域残差项逐步补偿了源预测和真实目标域响应之间的差异。残差过程进一步建模为在相同联合输入空间上定义的高斯过程:(50)?????(??,??)~?????(0,?????((??,??),??′,??′)))。由于残差建模仍然在参数和上下文的联合空间中进行,这个过程同样可以捕捉参数接近性和上下文相似性对目标域差异分布的联合影响。在迭代t时,设目标域观测数据集为??(??)??={(????,??????)}??????=1,其中??????=?????(????,????)表示目标域样本点(xi, ci)处的真实评估,nt是当前累积据此,转移后的目标替代模型的后验均值定义为(54)??(??)???(??,??)=?????(??,??)+???????(??)???(??,??)。这个表达式表明,当仅有少量目标域样本可用时,模型主要依赖于源先验来提供基线预测;随着收集到更多目标域样本,目标残差校正的作用变得越来越重要。对于预测不确定性,由于源随机过程没有被传播,源后验均值被视为固定的基线,因此目标域的预测方差仅由残差GP表征:(55)??2,(??)???(??,??)=??2,(??)???(??,??)。这种改进的公式去除了之前模糊的加性方差表达式。目标替代模型不再包含源后验方差或样本依赖系数λt在方差项中,因为λt仅缩放残差均值校正,而不缩放潜在的残差随机过程本身。因此,目标域中的转移替代模型可以写为(56)???(??)???(??,??)~???(??(??)???(??,??),??2,(??)???(??,??))。5.3. 算法工作流程提出的转移校准过程通过一个包含冷启动初始化和基于残差的迭代搜索的两阶段程序来实现。详细步骤总结在算法3中。下载:下载高分辨率图片(192KB)下载:下载完整尺寸图片算法3. 用于跨车辆校准的转移贝叶斯优化。在冷启动阶段,我们利用源域后验均值?????(??,?藏(????????)??)作为零样本预测器来识别高质量的初始点。具体来说,在参数空间生成一个候选池?????????????,并选择具有最低先验预测损失的m个点进行初始物理模拟。这种策略确保了起始观测集??(0)??位于在源域已被证明有效的区域,减少了早期探索的浪费。在迭代搜索阶段,使用累积的观测数据和残差过程逐步细化目标替代模型。在每次迭代t时,构建一个新的候选池??(??)???????????。下一个评估点通过最大化转移获取函数??(??)??????(??)来确定,该函数平衡了源域的结构先验和目标域的残差不确定性。信任系数λt确保了从先验引导的搜索到数据驱动细化的无缝过渡,随着目标域样本量的增加,从而能够精确适应目标车辆的具体特征。从优化的角度来看,这种工作流程允许目标域的搜索从一个更有信息量的起始区域开始,而不是进行纯粹的无信息探索。冷启动阶段主要提高了初始目标观测的质量,而随后的基于残差的迭代逐步增强了対目标车辆的局部适应。因此,所提出的算法在保持从源域继承的结构指导的同时,保留了在有限的评估预算下适应车辆特定目标域特征的足够灵活性。6. 数值实验本节从准确性、收敛效率和上下文依赖校准参数的物理可解释性等方面评估所提出的框架。6.1. 数据集和车辆初始化本研究使用了两个不同的真实世界数据集,代表不同的车辆平台和运行环境。特别是,BMW i3数据集用于物理模型和CBO机制的主要验证,而Tesla Model 3数据集用于TBO模块的迁移学习分析。我们物理模型中使用的电机损失图来自Koch等人(2021年)和MotorXP(2020年)的公开可用基准研究。(1) BMW i3数据集。第一个数据集包含来自BMW i3(60 Ah)的大约70个真实世界行程,由慕尼黑工业大学(Steinstraeter等人,2020年)发布。每个行程包含四个类别的同步高分辨率测量数据:(i) 环境状态(例如,环境温度和海拔),(ii) 车辆动态(例如,速度和油门位置),(iii) 电池状态(例如,电压、电流和充电状态),以及(iv) 客舱/供暖系统信号(例如,室内温度和供暖功率)。这个数据集的一个关键特性是其广泛的环境温度范围(-10°C至35°C),这为评估上下文信息是否能在变化较大的运行条件下改善校准提供了一个合适的测试平台。(2) Tesla Model 3数据集。第二个数据集是通过阿贡国家实验室(Livewire Data Platform,2024年)提供的Tesla Model 3自动驾驶在路上行驶的数据集。它包含大约60次行驶,覆盖了超过1000英里的 distance,主要是在郊区和高速公路条件下。与BMW i3数据集相比,这个数据集反映了更现代的动力系统和感知堆栈,结合了车辆CAN总线信号和来自摄像头和LiDAR的后处理感知数据。CBO框架中使用的上下文向量c包括五个变量:环境温度、天气条件、路线类型、行驶距离和前缀平均速度。对于BMW i3数据集,前四个变量在记录的数据中是自然存在的;对于Tesla Model 3数据集,环境温度和天气条件是根据行驶的日期和时间从历史天气记录中查询的。前缀平均速度是根据两个数据集的前10%的观测轨迹计算得出的。故意排除了需要知道完整未来轨迹的变量。分类变量被 one-hot 编码,所有连续变量都使用训练集统计信息进行了标准化。这两个数据集中这些上下文变量的协调在表1中总结。表1. 数据集中上下文变量的协调。上下文变量BMW i3Tesla Model 3环境温度数据集中自然存在从天气数据中查询天气条件数据集中自然存在从天气数据中查询路线类型数据集中自然存在数据集中自然存在行驶距离数据集中自然存在前缀平均速度根据前10%计算根据前10%计算对于两种车辆,初始物理参数都是从公开可用的技术规格和经验工程估计中得出的。这些值定义了可行的搜索范围,并作为物理能量模型的先验初始化。主要规格在表2中总结。对于参数空间,两种车辆都校准了相同的六个物理参数(fr、Cd、Ri、Pconst、kv、ηregen)。搜索范围是车辆特定的,但在相同的物理尺度上定义,并且在GP训练期间所有参数都标准化到[0,1]以确保兼容性。表2. BMW i3和Tesla Model 3的比较物理规格。参数符号BMW i3Tesla Model 3整车质量(kg)m14171847车轮半径(m)rw0.34980.33前面积(m2)Af2.382.22传动比iG9.6659.04传动效率ηG0.960.98最大电机功率(kW)Pmax125188最大电机扭矩(Nm)Mmax250430最大再生扭矩(Nm)Mrecup7780最大再生功率(kW)Precup5060电池容量(Wh)Ecap22,20075,000内阻(Ω)Ri0.07680.086.2. 物理模型基准测试和CBO评估本节重点关注BMW i3数据集,并评估所提出的物理模型以及基于CBO的校准策略。目标有两个:首先,量化通过改进模型结构引入的准确性提升;其次,验证上下文先验是否能够显著提高校准效率,与非上下文搜索相比。6.2.1. 基线和实验设置为了建立一个严格的基准,我们将所提出的框架与模型结构和优化器基线进行比较:MMPEVEM + BO:作为主要的模型结构基线,我们采用了亚琛工业大学(RWTH Aachen University)机械电子教学与研究领域(MMP)开发的MMPEVEM(MMP的电动汽车排放模型)(Koch等人,2021年)。作为集成到SUMO交通微模拟器中的最广泛使用的EV能量模型之一,MMPEVEM校准了四个核心参数:滚动阻力系数fr、空气动力阻力系数Cd、电池内阻Ri和一个常数辅助功率项Paux。这个基线是针对每次行程独立使用标准BO进行校准的。提出的 + 随机:所提出的模型使用来自可行搜索空间的单个随机参数样本进行评估。所提出的模型通过额外引入动态再生效率项ηregen和速度依赖的辅助功率系数kv来扩展MMPEVEM的公式。提出的 + LHS:所提出的模型带有LHS,其中最佳候选项是从空间填充设计中选择的。提出的 + BO:所提出的模型使用标准贝叶斯优化独立校准。提出的 + CBO:完整的框架,其中所提出的模型使用CBO策略进行校准。MLP:作为像Adedeji和Kabir(2023年)以及Topi?等人(2019年)这样的纯数据驱动的黑盒模型的代表,我们实现了一个多层感知器(MLP)回归器。该模型包含两个隐藏层,每个层各有512个神经元,使用层归一化和ReLU激活函数。输入向量包括瞬时速度、加速度和道路坡度,而输出是瞬时能量消耗Pinst。对于这种逐次行程的基线,每个10次测试行程分别独立划分为训练和测试子集,比例为80:20,以便在每次单独的行程中训练和评估MLP。c-MLP:为了评估纯数据驱动模型的跨行程泛化能力,我们还实现了一个上下文MLP(c-MLP)。在60次历史训练行程上训练一个c-MLP,并直接在10次保留的测试行程上进行评估,遵循与CBO相同的60/10行程级别划分。其输入包括三个时间序列特征(速度、加速度和道路坡度)以及CBO框架中使用的五个行程级别上下文特征:环境温度、天气条件、路线类型、行驶距离和前缀平均速度。实验设置旨在确保评估结果在不同时间和上下文尺度上的可靠性和泛化能力。对于BMW i3数据集,随机选择了60次行程来构建历史上下文库????????????,而剩余的10次行程构成测试集?????????????。对于特斯拉,我们考虑的是少样本迁移,因此所有行程都独立地被视为在线样本,没有保留单独的历史上下文集。为了严格评估框架的性能,对历史和测试阶段应用了不同的数据分割策略:•离线阶段:对于????????????中的每个行程,时间序列样本被分为80%用于离线参数识别和20%用于保留验证。•在线阶段:为了模拟车辆仅短暂行驶的真实世界在线适应情况,每个测试行程的最初10%被用作少样本校准数据,而行程的剩余90%作为未见的评估集。除非另有说明,所有迭代优化方法都受到最多??=50次评估的限制。为了减轻来自优化器初始化和获取过程的随机方差,每个实验使用五个不同的随机种子重复五次进行BO/LHS初始化和候选池生成。因此,报告的标准差反映了随机优化过程引入的方差。6.2.2. 准确性比较首先检查所提出的物理模型和CBO方法的准确性。使用五个指标来评估性能:加权平均绝对百分比误差(WMAPE)、平均绝对误差(MAE)、均方根误差(RMSE)、总能量误差和有效率。其中,MAE和RMSE衡量瞬时功率偏差,而总能量误差评估行程级别的累积能量一致性。有效率量化了保持在物理可接受操作范围内的模拟样本的百分比。表3总结了10次行程测试套件的预测性能。表3. BMW i3测试套件中能量模型和校准框架的性能比较。配置WMAPE (%)MAE (Wh)RMSE (Wh)总误差 (%)有效率 (%)MMPEVEM + BO19.52 ± 10.430.42 ± 0.180.56 ± 0.2414.02 ± 24.61100.00 ± 0.00提出的 + 随机32.73 ± 5.910.66 ± 0.231.00 ± 0.458.63 ± 14.7599.90 ± 0.17提出的 + LHS19.06 ± 6.340.37 ± 0.130.58 ± 0.185.76 ± 5.5399.96 ± 0.12提出的 + BO18.52 ± 5.770.38 ± 0.130.63 ± 0.246.94 ± 9.62100.00 ± 0.00提出的 + CBO17.27 ± 6.310.34 ± 0.130.53 ± 0.184.03 ± 6.3699.96 ± 0.13MLP34.05 ± 13.160.68 ± 0.241.13 ± 0.408.14 ± 10.71—c-MLP47.52 ± 19.880.94 ± 0.431.40 ± 0.5631.74 ± 37.11—从表3中可以得出几个观察结果。首先,提出的 + 随机设置的表现明显不如所有优化变体,这证实了物理模型本身在没有信息校准的情况下无法保证良好的准确性。其次,当使用相同的优化器时,提出的 + BO在行程级别的能量误差和一致性上始终优于MMPEVEM + BO,表明改进的物理模型具有更好的表示保真度。第三,提出的 + CBO配置实现了最佳的整体性能,具有最低的WMAPE(17.27%)、MAE(0.34 Wh)和RMSE(0.53 Wh)。最后,两种神经网络基线在两种评估设置下表现出显著不同的性能。逐次行程的MLP在每个10次测试行程中使用80:20的比例进行训练和测试,实现了34.05%的WMAPE。相比之下,c-MLP在60次历史行程上训练,并直接在10次保留的行程上进行评估,导致WMAPE为47.52%。c-MLP的相对较大误差主要是由于一些在训练数据中未能充分表示的运行条件导致的。例如,TripA14、TripB12和TripA15的MAPE值分别为67.77%、91.37%和66.85%。这些较大的误差表明,纯数据驱动的模型需要足够多样和具有代表性的训练数据才能可靠地泛化到以前未见或代表性不足的运行条件。相比之下,CBO在相同的跨行程设置下实现了17.27%的WMAPE,显示出更大的鲁棒性。请注意,MMPEVEM + BO在总能量误差(24.61%)上表现出非常大的标准差。MMPEVEM基线的不稳定性可以归因于两个因素。首先,其对再生制动的简化处理倾向于在负功率段引入系统偏差。其次,在短距离或低消耗的行程中(例如,TripA13),会出现小分母效应,其中较小的绝对偏差可能会放大基于百分比的能量误差。相比之下,所提出的公式结合上下文感知的初始化方法,在不同的操作条件下产生了更稳定的参数识别。图5通过测量能量消耗和模拟能量消耗之间的均衡图进一步说明了即时预测质量。基线MMPEVEM + BO在负功率区域显示出明显的分散,表明在能量回收事件期间存在系统性的不匹配。在提出的+随机设置中,高能量消耗下的散布显著增加,表明未经信息的参数选择无法捕捉到在激进加速或高负载条件下的动力系统行为。相比之下,提出的+CBO配置在整个操作范围内围绕??=??线形成了最紧密的集中,包括再生和以牵引为主的模式。所得的??2=0.969表明,在不同的驾驶操作下,校准后的参数集在物理上是一致的,而不仅仅适合主要的操作区域。
下载:下载高清晰度图像(304KB)
下载:下载全尺寸图像
图5. 不同校准框架下的即时功率一致性。
6.2.3. 优化效率
另一个重要的实际要求是校准策略是否能够足够快速地收敛,以便在线或接近在线部署。表4比较了所提出的CBO框架与在不同评估预算下的非上下文基线。
表4. 在不同评估预算下的效率和收敛性分析。
配置 迭代次数 WMAPE (%) MAE (Wh) RMSE (Wh) 总误差 (%) 有效百分比 (%)
空白单元 提出的+随机 13 2.73 ± 5.91 10.66 ± 0.23 1.00 ± 0.45 8.63 ± 14.75 99.90 ± 0.17
提出的+LHS 50 19.06 ± 6.34 0.37 ± 0.13 0.58 ± 0.18 5.76 ± 5.53 99.96 ± 0.12
提出的+LHS 100 17.37 ± 6.28 0.34 ± 0.13 0.54 ± 0.21 3.93 ± 6.69 99.92 ± 0.17
提出的+CBO(零次学习) 12 1.87 ± 9.08 0.45 ± 0.22 0.69 ± 0.34 8.77 ± 9.94 99.96 ± 0.13
提出的+CBO(少量学习) 10 17.53 ± 6.09 0.34 ± 0.12 0.53 ± 0.17 6.44 ± 7.07 100.00 ± 0.00
提出的+CBO(完全学习) 50 17.27 ± 6.31 0.34 ± 0.13 0.53 ± 0.18 4.03 ± 6.36 99.96 ± 0.13
结果清楚地显示了所提出的CBO策略的 warm-start 优势。在第一次评估(??=1)时,CBO(零次学习)的WMAPE为21.87%,远优于纯粹随机的基线(32.73%)。这表明历史上下文经验即使在没有任何实质性的特定行程适应之前也能提供高质量的先验估计。在少量学习的情况下,效率增益更为明显。仅通过10次评估,CBO(少量学习)的WMAPE就达到了17.53%,这已经接近于LHS在1000次评估时的17.37%。根据配对t检验(??=.87),CBO在10次评估(WMAPE = 17.53%)与CBO在50次评估(WMAPE = 17.27%)之间的差异并不显著。这个结果表明,上下文先验使得在前几次迭代中就能快速收敛,而额外的评估带来的改进可以忽略不计。
关于详细的时间消耗,所有实验都是在配备Intel Core i5-9400F CPU(6核,2.90GHz)和32GB RAM的工作站上进行的,没有使用GPU加速。每次物理模型评估的平均墙钟时间约为0.099秒。离线训练阶段需要大约207秒的模拟时间,加上大约36秒的GP超参数优化时间(约35次触发,每次约1.03秒),从而使得端到端的离线训练时间约为4.1分钟。这种离线成本是一次性投资,将在所有未来的行程中分摊。相比之下,在线少量学习校准只需要10次评估(约0.99秒),而冷启动的BO需要50次评估(约4.94秒),或者LHS需要1000次评估(约98.7秒)。我们强调,CBO报告的效率增益特指在校准后每次新行程所需的在线物理模型评估次数的减少。这里,物理模型评估是指在给定参数配置下的前向模拟,与收集新的车辆轨迹不同。LHS在1000次评估的结果是作为穷举搜索的参考,而不是作为竞争性的在线基线。
图6中的收敛轨迹提供了这种行为的补充视角。标准BO最终会有所改善,但由于缺乏信息丰富的先验和获取过程的随机性,会显示出可见的暂停。LHS的改进更为稳定,但只是通过蛮力的样本积累实现的。相比之下,CBO在前几次迭代中迅速下降,并迅速达到一个良好的性能区域,表明大部分优化负担已经被上下文条件化的先验所吸收。这正是在线适应所期望的行为:快速early收敛,而长时间的搜索带来的边际收益有限。
下载:下载高清晰度图像(641KB)
下载:下载全尺寸图像
图6. 不同校准策略的收敛性比较。
6.2.4. 上下文感知校准的物理可解释性
除了预测性能之外,所提出的框架还保留了一种物理上基于现实的可参数化方式,其中校准后的量可以直接与车辆能量模型的可解释组件相关联。由于优化器校准的是明确的模型参数而不是潜在的神经表示,因此可以检查得到的上下文-参数关系,以评估有效参数趋势是否与物理预期大体一致。
为了检验校准后的参数模式是否受到限制性或任意参数范围的影响,表5报告了六个校准参数的搜索范围以及在优化过程中实际探索的范围和最优值的中位数。搜索范围是基于车辆规格、车辆动力学和EV建模研究中报告的经验范围以及基本的物理可行性约束来确定的。探索的范围表明,优化器积极地在可行的参数空间内进行搜索,而不是局限于一个狭窄的区域。
表5. 校准参数的搜索范围、探索范围和最优值的中位数。
参数 搜索下限 搜索上限 探索最小值 探索最大值 最优中位数
fr 0.001 0 0.02 0 0.02 0 0.02 0 0.008
Cd 0.20 0 0.50 0 0.50 0 0.21 84
Ri [Ω] 0.005 0 0.50 0 0.50 0 0.21 50
Pconst [W] 20.00 1 0 0.10 99.94 50.0
kv 0 8 0.00 0.01 22 0 8 0.00 23.92
ηregen 0.30 0 1.00 0 0.30 1 1.00 0.96 33
为了评估校准参数的可识别性,我们使用了60次校准行程中收集的所有样本计算了成对的皮尔逊相关系数。如图7所示,任意两个参数之间的绝对相关性都低于0.15,表明这些参数没有严重混淆,可以根据驾驶循环中包含的信息合理区分。
图8总结了行程上下文与校准参数之间的关系。对于连续上下文(图8a),出现了一些与物理预期大体一致的趋势。首先,滚动阻力系数fr随着环境温度的降低而降低,这与轮胎物理学一致。较暖和的条件通常通过改变橡胶的柔顺性和充气压力来减少滞后损失。其次,辅助功率Paux与平均速度呈正相关。尽管在可用时HVAC负载是分开建模的,但更高的持续速度通常伴随着更大的热管理需求和增加的低电压系统活动,这可以纳入有效的辅助负载中。相比之下,校准后的电池内阻对环境温度的依赖性较弱。一个合理的解释是,车辆的热管理系统稳定了电池核心温度,从而减弱了环境条件对电化学电阻的直接影响。
下载:下载高清晰度图像(1MB)
下载:下载全尺寸图像
图8. 校准后的物理模型参数对不同行程上下文的敏感性。
对于分类上下文(图8b),天气和路线属性以与物理直觉大体一致的方式改变了参数分布,尽管效应大小在不同类别间有所不同。例如,雨天条件改变了有效的道路-轮胎相互作用以及与维持稳定运行相关的能量需求,这可能表现为滚动阻力和等效电池侧损失的改变。同样,不同的路线类型影响了牵引和再生机会之间的平衡,从而改变了摩擦和再生相关参数的最优值。一些类别仍然显示出相对较高的方差,这是由于在少数子集中样本数量有限所致。
尽管如此,我们将这些校准后的值解释为有效参数,它们共同捕捉了车辆在给定操作条件下的能量响应,而不是作为唯一可识别的物理常数。这项工作的主要贡献在于所提出框架的预测准确性和样本效率高的校准能力,物理一致性作为一个辅助的观察结果。
6.3. 跨车辆传递评估TBO
为了评估所提出的TBO方法在跨车辆设置中的传递校准准确性,我们将其与标准BO、LHS和Random Select进行了比较。所有方法都使用每个Tesla行程的前30%来进行参数识别,然后对剩余的70%未见部分和完整的100%行程进行评估。保留的结果是对参数预测泛化的严格测试,而完整行程的结果反映了校准的总体稳定性。
表6报告了在Tesla行程的保留的70%部分上的预测性能。由于这些部分完全排除了目标域参数校准,它们提供了评估跨车辆泛化的主要样本外证据。
表6. 在Tesla行程的保留的70%部分上的性能比较。
配置 WMAPE (%) MAE (Wh) RMSE (Wh) 总误差 (%) 有效百分比 (%)
随机 27.36 ± 3.78 0.96 ± 0.27 1.38 ± 0.37 16.58 ± 10.71 100.00 ± 0.00
LHS 27.04 ± 4.31 0.95 ± 0.27 1.32 ± 0.36 14.04 ± 10.02 100.00 ± 0.00
BO 27.16 ± 3.73 0.94 ± 0.24 1.36 ± 0.34 12.16 ± 8.38 100.00 ± 0.00
TBO 24.82 ± 3.95 0.86 ± 0.24 1.23 ± 0.34 12.65 ± 8.53 100.00 ± 0.00
TBO在所有四种方法中实现了最低的点预测误差。其WMAPE为24.82%,而标准BO为27.16%,相对减少了大约8.6%。MAE和RMSE也呈现出相同的模式:TBO将MAE从0.94 Wh降低到0.86 Wh,RMSE从1.36 Wh降低到1.23 Wh,分别对应大约8.5%和9.6%的减少。这些结果表明,传递的源域信息,加上目标域的残差校正,使得校准后的参数能够更准确地再现目标车辆未见部分的点能量消耗动态。
累积总能量指标呈现出略有不同的模式。BO实现了最低的保留总能量误差12.16%,而TBO获得了12.65%。相对于这种累积指标的大幅行程间变化,TBO在WMAPE、MAE和RMSE方面仍然具有明显优势。因此,保留的结果并不表明TBO在每个指标上都占主导地位;相反,它的主要优势在于在保持与标准BO相当的累积能量一致性的同时,提高了点预测的精确度。所有四种方法的有效率均为100.00%,表明校准后的参数集仍在物理模型的可行操作范围内。
表7提供了对完整Tesla行程的补充评估。因为每次行程的前30%用于目标域校准,这些结果应该被解释为对校准后总体稳定性的评估,而不是完全的样本外泛化测试。TBO再次实现了最低的点误差指标,WMAPE为23.89%,MAE为0.82 Wh,RMSE为1.19 Wh。与标准BO相比,这些分别相当于大约7.8%、7.5%和9.1%的减少。
表7. 在完整100% Tesla行程上的性能比较。
配置 WMAPE (%) MAE (Wh) RMSE (Wh) 总误差 (%) 有效百分比 (%)
随机 26.10 ± 3.46 0.90 ± 0.19 1.31 ± 0.29 12.87 ± 8.64 100.00 ± 0.00
LHS 25.80 ± 3.46 0.89 ± 0.19 1.27 ± 0.27 10.93 ± 7.55 100.00 ± 0.00
BO 25.90 ± 3.14 0.89 ± 0.17 1.31 ± 0.28 8.69 ± 6.10 100.00 ± 0.00
TBO 23.89 ± 3.23 0.82 ± 0.17 1.19 ± 0.26 8.69 ± 6.10 100.00 ± 0.00
对于累积能量一致性,TBO和BO在完整行程上基本没有区别:它们的总能量误差分别为8.69%和8.69%。这一结果很重要,因为它表明点预测精确度的提高并没有以牺牲行程级别的累积能量估计为代价。总的来说,保留和完整行程的评估表明,TBO的主要跨车辆好处在于一致降低了点预测误差,而累积能量准确性仍与标准BO相当。特别是,保留的70%结果提供了更强的证据,表明传递的替代模型可以识别出超出目标域校准段的参数配置。
6.3.2. 收敛性和样本效率分析
为了进一步比较在不同目标域预算下不同方法的优化效率,我们评估了LHS、标准BO和TBO在四个代表性阶段(即在选定的Tesla测试行程上进行5次、10次、20次和30次迭代后)所达到的最佳目标值分布。较低的目标值表明了更好的校准质量和更快的收敛速度。
如图9所示,在相同的评估预算下,TBO的收敛速度始终比LHS和标准BO快。在早期阶段,TBO已经实现了更低的中位数目标值,这主要得益于其冷启动设计。在中期阶段,TBO的收敛优势更加明显。其目标值分布下降得比两个基线都快,表明TBO可以用更少的目标域评估次数识别出低损失区域。这证实了传递的先验不仅改善了初始化,还加速了后续的顺序优化。到第30次迭代时,三种方法之间的差距变小了。这是预期之中的,因为随着优化预算的增加,标准BO也逐渐学会了目标域的响应结构,并接近了模型和搜索空间施加的性能限制。尽管如此,TBO仍然保持了更好或相当的最佳目标值水平。
下载:下载高清晰度图像(910KB)
下载:下载全尺寸图像箱线图比较了LHS、BO和TBO在选定的特斯拉测试行程中分别经过(a) 5次、(b) 10次、(c) 20次和(d) 30次迭代后所取得的最佳目标值。总体而言,这些结果表明TBO的主要优势在于收敛速度快和样本效率高。在目标域预算有限的情况下,TBO能够比LHS和标准BO更快地获得高质量的校准解,这对于实际的跨车辆校准任务来说尤其有价值。
6.3.3 信任增长超参数κ的敏感性分析
为了检验TBO是否强烈依赖于方程(53)中的信任增长超参数κ,我们在特斯拉目标车辆行程上评估了??∈{1,3,5,10,20,50}的不同值。在整个实验过程中,源域替代模型、目标域数据划分、参数搜索空间、评估预算和特定行程的随机种子保持不变,只有κ进行了变化。如表8所示,保留的70%数据的平均WMAPE范围为24.82%至26.63%,而整个行程的平均WMAPE范围为23.89%至25.86%。根据上述统计报告惯例,默认设置??=5在两个评估区间内均实现了最低的WMAPE。与??=5相比,??=5使保留的70%数据的WMAPE和整个行程的WMAPE分别增加了1.81个百分点和1.97个百分点。
表8. TBO对信任增长超参数κ的敏感性
| ?? | 保留的70%数据WMAPE (%) | 整个行程WMAPE (%) | 第30次评估时的最佳目标值 |
|------|-------------|----------------|-------------|
| 5 | 25.94?±?4.16 | 26.55 | 26.56 |
| 10 | 25.14?±?3.76 | 26.01 | 26.59 |
| 20 | 26.01?±?4.26 | 25.19?±?3.89 | 26.89 |
| 50 | 26.36?±?4.06 | 25.56?±?3.69 | 27.18 |
| ...... | ...... | ...... | ...... |
在第30次评估时,最佳目标值从??=5时的26.59提高到了??=50时的27.71,增加了1.12。这种模式与信任系数的增长行为一致。当??=5时,λt在5次目标域观测后从0.63增加到30次观测后的0.99;而当??=1时,残差校正几乎在5次观测后就被完全激活,此时????=0.99;当??=50时,即使在30次观测后,残差校正仍然相对较弱,此时????=0.45。因此,??=5在早期源域指导和后续目标域适应之间提供了一个适度的过渡。我们将它作为一个稳健且实用的默认值,因为它位于一个相对稳定的范围内,而不是将其解释为一个理论上推导出的或唯一最优的超参数。
7. 结论
本文提出了一个综合框架,通过整合精细的高保真物理模型与CBO和TBO架构来实现电动车能量估算。实验结果表明,CBO框架在瞬时水平上实现了17.27%的较低WMAPE,每次行程的总能量误差为4.03%。此外,在保留的70%特斯拉行程上的主要样本外评估显示WMAPE为24.82%,总能量误差为12.65%。这些高保真结果以及高效的零样本或小样本能力为实时导航和车队管理系统中的能量感知决策提供了坚实支持。尽管取得了这些进展,未来仍有许多研究方向可以探索。未来的工作应该采用更严格的上下文保留验证方法,即在训练中不包含特定的操作条件(如温度范围或天气类别),以更系统地评估模型对未见过的上下文环境的泛化能力。跨车辆评估也有限,个别目标行程观察到了负 transfer 效应。首先,通过数据集扩展和异构转移进行进一步验证是必要的。具体来说,使用更广泛的车辆架构和多样化的地理数据集来测试该框架,可以揭示所提方法在多车辆能量模型同时可用的场景中的可扩展性。其次,解决纯物理模型固有的精度上限仍然是一个优先事项;结合残差学习或物理信息神经网络是一个有前景的方向,可以捕捉未建模的非线性动态。第三,探索先进的上下文学习方法对于提高系统在遇到前所未有的驾驶条件或极端环境下的适应性具有重要意义。这些发展将有助于确保校准框架在现实世界应用中的长期可靠性和稳健性。未来的工作还将研究参数的可识别性以及跨更多数据集和车辆平台的更广泛统计验证。
CRediT作者贡献声明
刘成奇:撰写 - 原始稿、验证、方法论、数据整理、概念化
史欣瑜:撰写 - 审阅与编辑、验证、方法论、概念化
霍金宝:验证、方法论
张宇豪:可视化、调查
孙音可:撰写 - 审阅与编辑、调查
艾子博:验证
胡洪涛:方法论
刘志远:撰写 - 审阅与编辑、监督、资金获取