基于实际工况提取的电动公交车能耗估计方法

本发明属于城市公共交通运营管理，具体涉及基于实际工况提取的电动公交车能耗估计方法。

背景技术：

1、城市公共交通系统是传统化石能源消耗和温室气体排放的重要领域，在我国实现双碳战略目标的进程中具有巨大潜能。公共交通系统的电动化转型是城市交通系统可持续性低碳转型中最有效的措施之一。

2、然而，相比于传统的燃油公交车，电动公交存在续航里程短、充电时间长等明显缺陷，公交系统的电动化转型、电动公交系统的资源配置和运营(如用于替换传统车辆的电动公交车数量、充电基础设施配置等)，高度依赖可靠和易于计算的公交车能耗估计算法。

3、电动公交车的能耗受到多种因素的影响，按影响因素的不同大致可以分为两大类——与车辆本身相关的内部因素和与驾驶条件相关的外部因素。前者包括车辆部件的特性，车辆设计参数，效率和惯性，辅助设备使用等；后者包括道路类型和条件，环境和交通状况，驾驶行为等。

4、现有电动汽车能耗估计模型主要包含分析模型、统计模型和计算模型。其中，分析模型以车辆纵向动力学为基础，由可用效率图预估电机损耗。车辆纵向动力学从车辆动力学理论出发建模，计算在车轮处克服对立力所需的功率。后两种模型虽然精度更高，但往往需要提前收集更加全面和高精度微观汽车行驶数据，且由于模型构建往往脱离汽车动力学原理，无法有效估计制动或下坡行驶时回收的能量，难以具有普遍适用性。

5、现有技术中，别一鸣等公开了一种基于多元非线性回归的电动公交车行程能耗估计方法(cn 113361792a)，能够一定程度上降低能耗预测方法对输入较高数据采集频率和精度的依赖，但采用的能耗估计算法中包含了有限的影响因素，未考虑道路坡度、车辆速度等造成的额外电量消耗。

6、肖大伟等公开了一种纯电动汽车能耗计算方法(cn 105426672a)，基于电动汽车行车记录仪采集的数据提取车辆工况速度信息，进而构建了动力功耗模型和辅助系统功耗模型来计算汽车能耗，采用的算法多针对单车能耗计算，基于行车记录仪数据提取车速信息的方法在大规模电动汽车能耗估算场景中计算负担较大。

7、现有电动汽车能耗估计方法对收集的数据质量和数据采集设备均有较高的要求。然而实际运营过程中，不同公交公司对车辆数据的收集精细度和质量参差不齐；在电动化转型过程中，一些线路甚至无法收集电动公交车行驶数据，因此提出的方法也无法估计电动公交车能耗。因此，基于汽车动力学理论和基础行驶工况数据的分析模型，应用能够普遍获得的数据格式，能更好地满足能耗估计模型普遍适用性的要求。

技术实现思路

1、为了克服上述电动公交能耗估计方法的缺陷和不足，本发明提出了一种电动公交车能耗估计方法，能够基于全球广泛使用和发布的通用公交提要规范数据，考虑了道路类型和交通状况等外在因素，提取公交线路的实际行驶工况，基于汽车动力学原理构建包括制动回收能量在内的公交车能耗估计模型。该方法极大地降低了能耗估计模型对数据精细度的要求，具有广泛的适用性。

2、基于实际工况提取的电动公交车能耗估计方法，具体步骤如下：

3、步骤一、采集电动公交车基础运行数据；

4、步骤二、基于步骤一的数据估计电动公交车牵引力；

5、步骤三、采集电动公交车辅助负载功率数据；

6、步骤四、基于路线段的电动公交车能耗计算；

7、步骤五、基于步骤二、步骤三和步骤四的数据估算电动公交车行程能耗。

8、进一步，步骤一中采集电动公交车基础运行数据，具体内容为：

9、建模的标准公交车的相关参数及其符号如下：

10、参数：公交车质量；符号：m；

11、参数：正面面积；符号：af；

12、参数：空气密度；符号：ρ；

13、参数：空气阻力系数；符号：cd；

14、参数：固定滚动阻力系数；符号：cr0；

15、参数：可变滚动阻力系数；符号：cr2；

16、参数：车轮半径；符号：rw；

17、参数：传动系效率；符号：ηt；

18、参数：变速器转动比；符号：gr；

19、参数：转换器效率；符号：ηc；

20、参数：辅助负载；符号：paux；

21、步骤11、提取出公交车的质量、迎风面积和轮胎的可变滚动阻力；

22、步骤12、搜集资料，包括：公交车在线路不同位置对应的速度、高程信息(道路坡度)、空气密度、经纬度和地球曲率半径；

23、步骤13、计算公交路线总长度，沿着原始静态路线点位置进行插值，并以di为间隔对线路位置进行重新采样；

24、通过半正矢方法来计算公交线路中第i个控制点到i+1控制点间的距离d；

25、

26、

27、式中，为位置i的纬度，λi为位置i的经度，re为地球的曲率半径；

28、步骤14、道路坡度计算：

29、

30、式中，mi为道路坡度，hi为某点的高程。

31、步骤15、当多个实时公交位置与单个静态路线位置匹配时，选取该位置的平均值速度为该点的速度，当平均值数值明显高于平均行驶速度时，选取速度中值为该点对应的速度。

32、步骤16、采用二阶sg(savitzky-golay)滤波法过滤生成的速度曲线的噪音，对公交路线中的控制点速度数据进行平滑和连续性处理；

33、优选的，步骤16中采用最小二乘法，以公交路线中控制点i为中心的w个数据点为窗口，对离散的数据点进行平滑处理；

34、优选的，步骤16中在拟合公交路线中的点i+1为中心的新的二次曲线之前，使用该二次曲线来获得点i的过滤值。

35、步骤17、计算得到每个线路控制点位置的加速度ai和控制点间的时间跨度ti：

36、

37、

38、式中，vi是公交线路中控制点i的公交车速度，di为路线段的长度。

39、进一步，步骤二中电动公交车牵引力的计算方法为：

40、步骤21、空气阻力faer表达式为：

41、

42、式中，ρ是空气密度，af是公交车的迎风面积，cd是空气阻力系数，v是公交车速度。

43、步骤22、爬坡阻力fgra的表达式为：

44、fgra＝m*g*sin(θ)

45、式中，m为含乘客的公交车的质量，g为重力加速度常数，θ为道路的坡度。

46、步骤23、滚动阻力frol表达式为：

47、frol＝(cr0+cr2*v)*m*g

48、式中，cr0为固定阻力系数，cr2为公交车轮胎的可变滚动阻力系数。

49、步骤24、加速阻力facc表达式为：

50、facc＝a*m

51、式中，a是公交车加速度。

52、步骤25、步骤2.1-2.4中力的总和与公交车车轮在运行期间产生的牵引力fpro相平衡。

53、fpro＝facc+faer+frol+fgra

54、进一步，步骤三中采集电动公交车辅助负载功率，辅助负载包括为公交车提供加热、冷却或照明等功能的设备运行时所产生的负载。

55、进一步，步骤四中基于路线段的能耗计算，具体内容为：

56、电动机转速ωm可通过下式计算：

57、

58、式中，v为公交车速度、g为变速器的传动比，rw为车轮半径；

59、车轮所需扭矩tw表达式为：

60、tw＝fpro*rw

61、电机扭矩tm表达式为：

62、

63、式中，ηt传动系统的传递效率；

64、电机转速和扭矩在电机效率图中找到对应的电机效率ηm。

65、基于电动公交车传动系统能量流，电动公交车在一个路线段(每两个控制点间的路段长度)所消耗的电能，使用下式计算：

66、

67、式中，paux为辅助负载，ηm为电机效率，ηc为转换器效率，s为再生制动系统回收的能量比例。当推进力为负值时表示公共汽车正在制动，此时使用上式计算返回到电池系统的能量。

68、步骤五中电动公交车的行程能耗计算，包括：

69、步骤51、通过提取不同阶段的电池功率和时间长度，计算出整条线路的能源消耗e。计算方法如下：

70、

71、式中，i是公交线路的控制点数量，起点编号为0，终点编号为i。

72、能源消耗可根据能量强度进行描述，具体计算公式如下式：

73、

74、式中，drou是公交线路总长度，单位为km。

75、本发明的有益效果：

76、本发明提出了一种电动公交车能耗估计方法，基于汽车动力学理论和通用的公交运行数据格式，融合公交线路静态数据，能够同时考虑影响汽车能耗的动力学特征和交通状况，极大地降低能耗预测方法对高频率、高精度数据采集技术的依赖，减小模型计算负担，具有广泛的适用性。