CN107169273B - 计及延时和v2g充电模式的电动汽车充电功率预测方法 - Google Patents
计及延时和v2g充电模式的电动汽车充电功率预测方法 Download PDFInfo
- Publication number
- CN107169273B CN107169273B CN201710310372.8A CN201710310372A CN107169273B CN 107169273 B CN107169273 B CN 107169273B CN 201710310372 A CN201710310372 A CN 201710310372A CN 107169273 B CN107169273 B CN 107169273B
- Authority
- CN
- China
- Prior art keywords
- charging
- time
- power
- electric vehicle
- electric
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 50
- 238000007599 discharging Methods 0.000 claims abstract description 48
- 238000004088 simulation Methods 0.000 claims abstract description 20
- 230000004044 response Effects 0.000 claims abstract description 16
- 230000008901 benefit Effects 0.000 claims abstract description 6
- 238000012360 testing method Methods 0.000 claims abstract description 5
- 238000000342 Monte Carlo simulation Methods 0.000 claims description 23
- 230000006399 behavior Effects 0.000 claims description 18
- 230000006870 function Effects 0.000 claims description 18
- 230000005611 electricity Effects 0.000 claims description 16
- 230000003111 delayed effect Effects 0.000 claims description 12
- 230000009471 action Effects 0.000 claims description 9
- 238000005315 distribution function Methods 0.000 claims description 9
- 230000000694 effects Effects 0.000 claims description 7
- 230000002035 prolonged effect Effects 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 5
- 230000008859 change Effects 0.000 claims description 3
- 230000005284 excitation Effects 0.000 claims description 3
- 238000009499 grossing Methods 0.000 claims description 3
- 239000004576 sand Substances 0.000 claims description 3
- 230000005540 biological transmission Effects 0.000 description 5
- 238000011161 development Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 2
- 230000000638 stimulation Effects 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- 229910052799 carbon Inorganic materials 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000004134 energy conservation Methods 0.000 description 1
- 238000004146 energy storage Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000006855 networking Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
- G06Q10/06315—Needs-based resource requirements planning or analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0637—Strategic management or analysis, e.g. setting a goal or target of an organisation; Planning actions based on goals; Analysis or evaluation of effectiveness of goals
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Marketing (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Educational Administration (AREA)
- Quality & Reliability (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Development Economics (AREA)
- Health & Medical Sciences (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Charge And Discharge Circuits For Batteries Or The Like (AREA)
- Secondary Cells (AREA)
- Electric Propulsion And Braking For Vehicles (AREA)
Abstract
本发明公布了一种计及延时和V2G充电模式的电动汽车充电功率预测方法,本发明用于研究电动汽车充放电随机模型,并通过模拟得到电动汽车的充放电功率。本发明首先针对电动汽车的时空不确定性进行分析,建立电动汽车充电功率的动态概率模型。其次,选择合适的充电方式,针对普通充电模式、延时充电模式和V2G充电模式对电动汽车充电功率特性进行建模。最后,在3种充电模式下对单台电动汽车进行仿真测试,计算各时段电动汽车负荷的期望、标准差。本发明能够获取需求响应背景下3种充电方式的电动汽车充放电功率期望值,具有原理简单,实现方便的优点,所得概率分布特性可以使调度人员对电动汽车充放电集群负荷有更为清晰的认识。
Description
技术领域
本发明属于电力系统运行分析和负荷预测领域,涉及一种计及延时和V2G充电模式的电动汽车充电功率预测方法,适用于研究需求响应背景下同时计及电动汽车在时间、空间以及驾驶行为等随机性因素影响下的充放电随机模型,并通过模拟得到电动汽车在延时和V2G充电模式下的充放电功率。
背景技术
电动汽车作为一种节能减排的有效手段,以电力为驱动,是一种新型的交通工具,为实现低碳生活、缓解能源危机提供新的契机。电动汽车因其基本无尾气排放得到了高度重视,也成为未来汽车发展的潮流,因而对其发展提出了迫切的要求。近年来,电力储能、充电、电机控制等技术的快速发展,充电设施的合理建设,都推进了电动汽车的普及进程。未来电动汽车将大规模接入电网进行充电,对电力系统的调度和运行产生深远影响。由于电动汽车在配电网接入,电动汽车的充电行为对配电网有直接影响,然而随着电动汽车数量的增加,电动汽车充电所产生的集群效应将对输电网产生影响。因此,在输电网的调度运行中,有必要考虑电动汽车集群效应的影响。如果大规模电动汽车在电网负荷高峰时段进行充电,将导致负荷“峰上加峰”,给电网运行带来一定负担。
作为一类特殊负荷,电动汽车在时间、空间及行为上均存在大量的不确定性,其是否可双向送电、司机的驾驶行为的差异都将产生巨大随机特性,大规模电动汽车的无序充电将引发负荷的快速增长,电动汽车在负荷高峰时段的充电行为将造成“峰上加峰”,给电力系统的运行带来挑战。解决上述问题的关键在于引导电动汽车的有序充电,即协调充电时间。针对此,国内外已有学者进行了相关研究,电动汽车入网(vehicle to grid,V2G)的概念也应运而生。有研究指出,如果将电动汽车以一个移动分布式储能单元的形式接入电网,并根据系统需要进行充放电方式转换,可以达到削峰填谷的作用,即电网负荷低谷时充电,高峰时放电,从而有效减轻电网的运行负担。在需求响应视角下,研究同时计及电动汽车在时间、空间以及驾驶行为等随机性因素影响下的充放电随机模型以及如何得到电动汽车在不同充电方式下充放电功率具有现实意义。
发明内容
发明目的:本发明的目的在于针对现有技术的不足,提出了一种计及延时和V2G充电模式的电动汽车充电功率预测方法,适用于研究需求响应背景下同时计及电动汽车在时间、空间以及驾驶行为等随机性因素影响下的充放电随机模型,并通过模拟得到电动汽车在延时和V2G充电模式下的充放电功率。
技术方案:本发明提供了一种计及延时和V2G充电模式的电动汽车充电功率预测方法,包括以下步骤:
步骤1:针对电动汽车的时空不确定性进行分析,建立电动汽车充电功率的动态概率模型;
步骤2:选择合适的充电方式,针对普通充电模式、延时充电模式和V2G充电模式对电动汽车充电功率特性进行建模,具体建模过程如下:
201)普通充电方式指电动汽车充电行为不受电价的影响,即电动汽车不改变其开始充电时间也不参与V2G放电,在普通充电方式下,对于单台电动汽车模拟其需要充电电量,结合电动汽车的充电功率和充电效率,得到其充电持续时间:
式中:TC表示充电持续时间,ηC表示充电效率,D表示日行驶里程数,PC为电动汽车充电功率,W100表示百公里耗电量;
202)采用蒙特卡罗模拟法,随机产生满足分布的最后一次返回时刻和日行驶里程数样本,模拟得到各时段单台电动汽车充电功率的期望,由于主要考虑大量电动汽车接入后,进行充放电的集群效应,因此,通过蒙特卡罗模拟法得到单台电动汽车的充电功率后,需要将充电功率进行累加,得到大量电动汽车的充电功率:
式中PEV[t]表示t时段母线上电动汽车的总充放电功率,nEV表示电动汽车总数量,PCi[t]表示第i辆电动汽车在t时段的充电功率;
203)通过蒙特卡罗模拟法,得到电动汽车充电功率的样本数据,为了进一步分析充电功率的概率特性,需要得到其概率分布曲线和概率密度曲线,为避免主观假设分布与实际分布存在偏差时导致所得结果产生较大的误差,采用非参数核密度估计,该模型不需要任何概率分布形式的假设,仅仅假设充电功率的概率分布函数和概率密度函数存在且连续可微,更客观地反映充电功率的概率特性;
概率X1,X2,…,XN是蒙特卡罗模拟法得到的N个充电功率的样本,假设充电功率的概率密度函数为f(x),则此概率密度函数的核估计为:
式中:h表示带宽,即窗宽或平滑系数,N表示样本容量,K(·)表示核函数;
204)延时充电方式指电动汽车车主在电价激励的作用下,改变其开始充电时刻,即将开始充电时刻转移至低谷时刻,延长响应的充电时间,节省充电费用,但不参与V2G放电,将一天24h分为高峰时段和低谷时段,高峰时段从Tpeak1时刻开始到Tpeak2时刻结束,其余时刻为低谷时刻,当电动汽车开始充电时刻处于高峰时段时,电动汽车不进行充电,等低谷时段再进行充电,通过延时充电,使一部分高峰时段的电动汽车负荷转移至低谷时段,实现削峰填谷,平抑电网总负荷量;
采用蒙特卡罗模拟法,随机产生满足分布的最后一次返回时刻和日行驶里程数样本,当最后一次返回时刻处于高峰时段时,将其开始充电时间转移至低谷时段,由此,模拟得到各时段单台电动汽车充电功率的期望、标准差和概率分布,将单台电动汽车的充电功率进行累加,得到大量电动汽车的充电功率;
205)V2G充电方式指电动汽车车主在电价激励的作用下,在高峰时刻作为电源向电网放电,并在低谷时刻作为负荷由电网向其充电,由于峰谷时段存在电价差,车主利用该电价差,节省总充电费用,甚至实现盈利,但相应的也会延长充电时间;
考虑到对电动汽车电池的保护以及车主的感受,电动汽车电池的最大放电状态不小于其容量的20%,即电动汽车在高峰时段接入电网,且入网的荷电状态大于电池容量的20%时,才可对电网进行放电,因而,结合电动汽车的日行驶里程数其持续放电至电池容量的20%所需的时间:
式中Tf0.2表示电动汽车放电至电池容量的20%所需时间,Smax表示电动汽车电池的总电量,ηdisc表示电动汽车的放电效率,Pdisc表示电动汽车的放电功率;
206)当电动汽车接入电网的荷电状态大于电池容量的20%时,其放电时间由开始充电时刻、放电至电池容量的20%所需的时间、高峰时段共同决定,放电时间如下:
式中Tf表示电动汽车的放电持续时间,Tpeak2为高峰时段右侧端点,TS为最后一次返回时刻;
根据电动汽车的放电时间,结合电动汽车的放电参数可以得到其充电持续时间:
式中TC表示电动汽车的充电持续时间;
207)采用蒙特卡罗模拟法,随机产生满足分布的最后一次返回时刻和日行驶里程数样本,对于高峰时段接入电网的电动汽车,结合其入网荷电状态对其是否进行放电进行判断,并计算其放电时间,到低谷时段,对电动汽车进行充电,由此,模拟得到各时段单台电动汽车充电功率的期望、标准差和概率分布;
208)由于输电网相比配电网具有更高的电压等级,而电动汽车往往不会直接接入输电网,电动汽车对输电网的影响要相对间接,因而对电动汽车建模的目的是得到大量电动汽车接入后充放电行为的集群特征,即母线上所有电动汽车在t时段的充放电功率之和PEV[t]:
步骤3:在步骤2的基础上,在3种充电模式下对单台电动汽车进行仿真测试,计算各时段电动汽车负荷的期望、标准差。
所述步骤1包括以下步骤:
步骤101:电动汽车的充电功率在时间和空间上均有较强的不确定性,针对电动汽车的时空不确定性进行分析是建立充电功率特性模型的基础,其中用户行为是影响电动汽车充电功率的关键,电动汽车的充电特性受电动汽车接入电网的时刻、需要充电电量和充电功率3方面因素的影响,本发明考虑私用调度汽车作为电网可调度资源,有多种充电方式,并以V2G的方式参与电网调度;
步骤102:电动汽车接入电网时刻与用户使用行为有密切关系,假设用户最后一次行驶返回时刻即为电动汽车的入网时刻,同时根据美国能源部2008年对电动汽车用户的调查统计报告,用户的最后一次返回时刻近似服从正态分布,最后一次返回时刻x的概率密度函数如下:
式中,TS表示最后一次返回时刻,μS和σS分别表示正态分布的期望和标准差;
步骤103:假设电动汽车一次充电即全部充满,则电动汽车当日的充电电量与其当日的行驶里程有密切关系,即电动汽车当日从电网获取的电能全部转化为当日行驶里程所消耗的能量,根据调查报告,电动汽车的日行驶里程数服从对数正态分布,电动汽车的日行驶里程数的概率密度函数如下:
式中:D表示日行驶里程数,μD和σD分别表示lnD的期望和标准差;
得到日行驶里程数后,结合电动汽车百公里电耗,得到需要充电电量如下:
式中:WC表示电动汽车需要充电电量,W100表示百公里耗电量;
步骤104:假设考虑电动汽车以常规慢充方式进行充电,由于起始充电阶段和结束充电阶段相较整个充电过程的时间较短,可以忽略,因此,假设电动汽车充电和放电均采用恒功率充放电方式进行,其充电功率为PC,放电功率为Pdisc;由于电动汽车模型是在峰谷电价的背景下建立的,电动汽车车主可以根据自身需要,采用不同的充电方式,以实现自身利益的最大化,即电动汽车对电价产生不同的响应方式,电动汽车的充电方式包括普通充电方式、延时充电方式和V2G充电方式3种。
所述步骤3包括以下步骤:
步骤301:单台电动汽车接入系统进行充放电时,假设高峰时段的开始时刻为Tpeak1=9,结束时刻Tpeak2=22,在步骤2的基础上进行蒙特卡洛模拟,得到单台电动汽车的充放电负荷期望值以及不同充电方式下单台电动汽车各时段充放电负荷期望曲线;
步骤302:在步骤301的基础上,对模拟得到的N个充电功率样本采用非参数核密度估计,求得各充电方式下100台电动汽车在各时段的充放电集群负荷的概率分布函数,通过这些概率分布函数的曲线图可以对电动汽车充放电集群负荷的概率特性有清晰的认识。
有益效果:与现有技术相比,本发明具有如下优点和技术效果:
(1)同时计及电动汽车在时间、空间以及驾驶行为等随机性因素影响下的充放电功率的随机特性,重点考虑用户的出行习惯以及驾驶习惯;
(2)在需求响应背景下,针对普通充电模式、延时充电模式和V2G充电模式3种不同充电模式对电动汽车充电功率特性进行建模,通过蒙特卡洛模拟法得到大量电动汽车的充电功率;
(3)获得单台电动汽车负荷的期望、标准差以及各时段充放电负荷期望曲线,拟合100台电动汽车在各时段充放电负荷概率分布曲线,所得概率分布特性可以使调度人员对电动汽车充放电集群负荷有更为清晰的认识。
附图说明
图1为计及延时和V2G充电模式的电动汽车充电功率预测方法流程图;
图2为普通充电方式下电动汽车功率特性模拟流程图;
图3为延时充电方式下电动汽车功率特性模拟流程图;
图4为V2G充电方式下电动汽车功率特性模拟流程图;
图5为单台电动汽车充放电负荷期望曲线图;
图6为普通充电方式下电动汽车充放电负荷概率分布曲线;
图7为延时充电方式下电动汽车充放电负荷概率分布曲线;
图8为V2G充电方式下电动汽车充放电负荷概率分布曲线。
具体实施方式
以下结合附图和实例对本发明的实施作进一步说明,但本发明的实施和包含不限于此。
一种计及延时和V2G充电模式的电动汽车充电功率预测方法,包括以下如下步骤:
步骤1:针对电动汽车的时空不确定性进行分析,建立电动汽车充电功率的动态概率模型;
步骤2:选择合适的充电方式,针对普通充电模式、延时充电模式和V2G充电模式对电动汽车充电功率特性进行建模;
步骤3:在步骤2的基础上,在3种充电模式下对单台电动汽车进行仿真测试,计算各时段电动汽车负荷的期望、标准差;
进一步,所述步骤1包括以下步骤:
步骤101:电动汽车的充电功率在时间和空间上均有较强的不确定性,针对电动汽车的时空不确定性进行分析是建立充电功率特性模型的基础,其中用户行为是影响电动汽车充电功率的关键。电动汽车的充电特性主要受电动汽车接入电网的时刻、需要充电电量和充电功率等3方面因素的影响。本发明主要考虑私用调度汽车作为电网可调度资源,有多种充电方式,并以V2G的方式参与电网调度。
步骤102:电动汽车接入电网时刻与用户使用行为有密切关系,假设用户最后一次行驶返回时刻即为电动汽车的入网时刻。同时根据美国能源部2008年对电动汽车用户的调查统计报告,用户的最后一次返回时刻近似服从正态分布。最后一次返回时刻的概率密度函数如下:
式中,TS表示最后一次返回时刻,μS和σS分别表示正态分布的期望和标准差。
步骤103:假设电动汽车一次充电即全部充满,则电动汽车当日的充电电量与其当日的行驶里程有密切关系,即电动汽车当日从电网获取的电能全部转化为当日行驶里程所消耗的能量。根据调查报告,电动汽车的日行驶里程数服从对数正态分布。电动汽车的日行驶里程数的概率密度函数如下:
式中:D表示日行驶里程数,μD和σD分别表示lnD的期望和标准差;
得到日行驶里程数后,结合电动汽车百公里电耗,得到需要充电电量如下:
式中:WC表示电动汽车需要充电电量,W100表示百公里耗电量。
步骤104:假设本发明考虑电动汽车以常规慢充方式进行充电。由于起始充电阶段和结束充电阶段相较整个充电过程的时间较短,可以忽略。因此,假设电动汽车充电和放电均采用恒功率充放电方式进行,其充电功率为PC,放电功率为Pdisc。由于本发明的电动汽车模型是在峰谷电价的背景下建立的,电动汽车车主可以根据自身需要,采用不同的充电方式,以实现自身利益的最大化,即电动汽车对电价产生不同的响应方式。电动汽车的充电方式包括普通充电方式、延时充电方式和V2G充电方式等3种。
进一步:所述步骤2包括以下步骤:
步骤201:普通充电方式指电动汽车充电行为不受电价的影响,即电动汽车不改变其开始充电时间也不参与V2G放电。在普通充电方式下,对于单台电动汽车模拟其需要充电电量,结合电动汽车的充电功率和充电效率,得到其充电持续时间:
式中:TC表示充电持续时间,ηC表示充电效率。
步骤202:采用蒙特卡罗模拟法,随机产生满足分布的最后一次返回时刻和日行驶里程数样本,模拟得到各时段单台电动汽车充电功率的期望,由于本发明主要考虑大量电动汽车接入后,进行充放电的集群效应,因此,通过蒙特卡罗模拟法得到单台电动汽车的充电功率后,需要将充电功率进行累加,得到大量电动汽车的充电功率:
式中:PEV[t]表示t时段母线上电动汽车的总充放电功率,nEV表示电动汽车总数量,PCi[t]表示第i辆电动汽车在t时段的充电功率。
步骤203:通过蒙特卡罗模拟,得到电动汽车充电功率的样本数据。为了进一步分析充电功率的概率特性,需要得到其概率分布曲线和概率密度曲线。为避免主观假设分布与实际分布存在偏差时导致所得结果产生较大的误差,本发明采用非参数核密度估计,该模型不需要任何概率分布形式的假设,仅仅假设充电功率的概率分布函数和概率密度函数存在且连续可微,更客观地反映充电功率的概率特性。
概率X1,X2,…,XN是蒙特卡罗法模拟得到的N个充电功率的样本,假设充电功率的概率密度函数为f(x),则此概率密度函数的核估计为:
式中:h表示带宽(窗宽或平滑系数),N表示样本容量,K(·)表示核函数。
步骤204:延时充电方式指电动汽车车主在电价激励的作用下,改变其开始充电时刻,即将开始充电时刻转移至低谷时刻,延长响应的充电时间,节省充电费用,但不参与V2G放电。将一天24h分为高峰时段和低谷时段,高峰时段从Tpeak1时刻开始到Tpeak2时刻结束,其余时刻为低谷时刻。当电动汽车开始充电时刻处于高峰时段时,电动汽车不进行充电,等低谷时段再进行充电。通过延时充电,使一部分高峰时段的电动汽车负荷转移至低谷时段,实现削峰填谷,平抑电网总负荷量。
采用蒙特卡罗模拟法,随机产生满足分布的最后一次返回时刻和日行驶里程数样本。当最后一次返回时刻处于高峰时段时,将其开始充电时间转移至低谷时段。由此,模拟得到各时段单台电动汽车充电功率的期望、标准差和概率分布。将单台电动汽车的充电功率进行累加,得到大量电动汽车的充电功率。
步骤205:V2G充电方式指电动汽车车主在电价激励的作用下,在高峰时刻作为电源向电网放电,并在低谷时刻作为负荷由电网向其充电。由于峰谷时段存在电价差,车主利用该电价差,节省总充电费用,甚至实现盈利,但相应的也会延长充电时间。
考虑到对电动汽车电池的保护以及车主的感受,电动汽车电池的最大放电状态不小于其容量的20%,即电动汽车在高峰时段接入电网,且入网的荷电状态(state ofcharge,SOC)大于电池容量的20%时,才可对电网进行放电。因而,结合电动汽车的日行驶里程数其持续放电至电池容量的20%所需的时间:
式中:Tf0.2表示电动汽车放电至电池容量的20%所需时间,Smax表示电动汽车电池的总电量,ηdisc表示电动汽车的放电效率,Pdisc表示电动汽车的放电功率。
步骤206:当电动汽车接入电网的SOC大于电池容量的20%时,其放电时间由开始充电时刻、放电至电池容量的20%所需的时间、高峰时段等共同决定。放电时间如下:
式中:Tf表示电动汽车的放电持续时间。
根据电动汽车的放电时间,结合电动汽车的放电参数可以得到其充电持续时间:
式中:TC表示电动汽车的充电持续时间。
步骤207:采用蒙特卡罗模拟法,随机产生满足分布的最后一次返回时刻和日行驶里程数样本。对于高峰时段接入电网的电动汽车,结合其入网SOC对其是否进行放电进行判断,并计算其放电时间。到低谷时段,对电动汽车进行充电。由此,模拟得到各时段单台电动汽车充电功率的期望、标准差和概率分布。
步骤208:由于输电网相比配电网具有更高的电压等级,而电动汽车往往不会直接接入输电网,电动汽车对输电网的影响要相对间接,因而本发明对电动汽车建模的目的是得到大量电动汽车接入后充放电行为的集群特征,即母线上所有电动汽车在t时段的充放电功率之和PEV[t]:
进一步:所述步骤3包括以下步骤:
步骤301:单台电动汽车接入系统进行充放电时,假设高峰时段的开始时刻为Tpeak1=9,结束时刻Tpeak2=22。在步骤2的基础上进行蒙特卡洛模拟,得到单台电动汽车的充放电负荷期望值以及不同充电方式下单台电动汽车各时段充放电负荷期望曲线。
步骤302:在步骤301的基础上,对模拟得到的N个充电功率样本采用非参数核密度估计,求得各充电方式下100台电动汽车在各时段的充放电集群负荷的概率分布函数,通过这些概率分布曲线可以对电动汽车充放电集群负荷的概率特性有清晰的认识。
工作原理:本发明首先针对电动汽车的时空不确定性进行分析,重点考虑用户的出行习惯以及驾驶习惯,在此基础上建立电动汽车充电功率的动态概率模型。其次,在需求响应背景下,根据用户与电力公司签订的相关协议及其他因素,选择合适的充电方式,包括普通充电模式、延时充电模式和V2G充电模式,针对3种不同的充电模式对电动汽车充电功率特性进行建模,通过蒙特卡洛模拟法得到大量电动汽车的充电功率。最后,在3种充电模式下对单台电动汽车进行仿真测试,计算各时段电动汽车负荷的期望、标准差并且最终得到单台电动汽车在各时段的充放电负荷期望曲线以及100台电动汽车充放电集群负荷概率分布曲线。
算例分析
单台电动汽车接入系统进行充放电时,假设高峰时段的开始时刻为Tpeak1=9,结束时刻Tpeak2=22,其中单台电动汽车充放电参数数值:μS取17.6h,σS取3.4h;μD取3.2km,σD取0.88km;W100取25kWh/100km,PC取4.5kW,ηC取90%,Smax取50kWh,Pdisc取4.5Kw,ηdisc取90%;
按照图2—图4所示流程图进行模拟,得到单台电动汽车的充放电负荷期望,充放电负荷期望曲线如图5所示。采用非参数核密度估计,分别求得普通充电模式、延时充电模式和V2G充电模式3种不同充电模式下100台电动汽车的充放电集群负荷的概率分布函数,如图6-图8所示,通过这些概率分布曲线可以对电动汽车充放电集群负荷的概率特性有清晰的认识,在需求响应背景下,由于峰谷电价的实施,电动汽车车主可以根据自身情况,选择合适的充电方式。不同的峰谷电价激励作用下,车主会产生不同的响应情况,电网公司可以据此对峰谷电价进行调节,以达到不同的响应目标。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明保护的范围之内。
Claims (3)
1.一种计及延时和V2G充电模式的电动汽车充电功率预测方法,其特征在于,包括以下步骤:
步骤1:针对电动汽车的时空不确定性进行分析,建立电动汽车充电功率的动态概率模型;
步骤2:选择合适的充电方式,针对普通充电模式、延时充电模式和V2G充电模式对电动汽车充电功率特性进行建模,具体建模过程如下:
201)普通充电方式指电动汽车充电行为不受电价的影响,即电动汽车不改变其开始充电时间也不参与V2G放电,在普通充电方式下,对于单台电动汽车模拟其需要充电电量,结合电动汽车的充电功率和充电效率,得到其充电持续时间:
式中:TC表示充电持续时间,ηC表示充电效率,D表示日行驶里程数,PC为电动汽车充电功率,W100表示百公里耗电量;
202)采用蒙特卡罗模拟法,随机产生满足分布的最后一次返回时刻和日行驶里程数样本,模拟得到各时段单台电动汽车充电功率的期望,由于主要考虑大量电动汽车接入后,进行充放电的集群效应,因此,通过蒙特卡罗模拟法得到单台电动汽车的充电功率后,需要将充电功率进行累加,得到大量电动汽车的充电功率:
式中PEV[t]表示t时段母线上电动汽车的总充放电功率,nEV表示电动汽车总数量,PCi[t]表示第i辆电动汽车在t时段的充电功率;
203)通过蒙特卡罗模拟法,得到电动汽车充电功率的样本数据,为了进一步分析充电功率的概率特性,需要得到其概率分布曲线和概率密度曲线,为避免主观假设分布与实际分布存在偏差时导致所得结果产生较大的误差,采用非参数核密度估计,该模型不需要任何概率分布形式的假设,仅仅假设充电功率的概率分布函数和概率密度函数存在且连续可微,更客观地反映充电功率的概率特性;
概率X1,X2,…,XN是蒙特卡罗模拟法得到的N个充电功率的样本,假设充电功率的概率密度函数为f(x),则此概率密度函数的核估计为:
式中:h表示带宽,即窗宽或平滑系数,N表示样本容量,K(·)表示核函数;
204)延时充电方式指电动汽车车主在电价激励的作用下,改变其开始充电时刻,即将开始充电时刻转移至低谷时刻,延长响应的充电时间,节省充电费用,但不参与V2G放电,将一天24h分为高峰时段和低谷时段,高峰时段从Tpeak1时刻开始到Tpeak2时刻结束,其余时刻为低谷时刻,当电动汽车开始充电时刻处于高峰时段时,电动汽车不进行充电,等低谷时段再进行充电,通过延时充电,使一部分高峰时段的电动汽车负荷转移至低谷时段,实现削峰填谷,平抑电网总负荷量;
采用蒙特卡罗模拟法,随机产生满足分布的最后一次返回时刻和日行驶里程数样本,当最后一次返回时刻处于高峰时段时,将其开始充电时间转移至低谷时段,由此,模拟得到各时段单台电动汽车充电功率的期望、标准差和概率分布,将单台电动汽车的充电功率进行累加,得到大量电动汽车的充电功率;
205)V2G充电方式指电动汽车车主在电价激励的作用下,在高峰时刻作为电源向电网放电,并在低谷时刻作为负荷由电网向其充电,由于峰谷时段存在电价差,车主利用该电价差,节省总充电费用,甚至实现盈利,但相应的也会延长充电时间;
考虑到对电动汽车电池的保护以及车主的感受,电动汽车电池的最大放电状态不小于其容量的20%,即电动汽车在高峰时段接入电网,且入网的荷电状态大于电池容量的20%时,才可对电网进行放电,因而,结合电动汽车的日行驶里程数其持续放电至电池容量的20%所需的时间:
式中Tf0.2表示电动汽车放电至电池容量的20%所需时间,Smax表示电动汽车电池的总电量,ηdisc表示电动汽车的放电效率,Pdisc表示电动汽车的放电功率;
206)当电动汽车接入电网的荷电状态大于电池容量的20%时,其放电时间由开始充电时刻、放电至电池容量的20%所需的时间、高峰时段共同决定,放电时间如下:
式中Tf表示电动汽车的放电持续时间,Tpeak2为高峰时段右侧端点,TS为最后一次返回时刻;
根据电动汽车的放电时间,结合电动汽车的放电参数可以得到其充电持续时间:
式中TC表示电动汽车的充电持续时间;
207)采用蒙特卡罗模拟法,随机产生满足分布的最后一次返回时刻和日行驶里程数样本,对于高峰时段接入电网的电动汽车,结合其入网荷电状态对其是否进行放电进行判断,并计算其放电时间,到低谷时段,对电动汽车进行充电,由此,模拟得到各时段单台电动汽车充电功率的期望、标准差和概率分布;
208)由于输电网相比配电网具有更高的电压等级,而电动汽车往往不会直接接入输电网,电动汽车对输电网的影响要相对间接,因而对电动汽车建模的目的是得到大量电动汽车接入后充放电行为的集群特征,即母线上所有电动汽车在t时段的充放电功率之和PEV[t]:
步骤3:在步骤2的基础上,在3种充电模式下对单台电动汽车进行仿真测试,计算各时段电动汽车负荷的期望、标准差。
2.根据权利要求1所述的计及延时和V2G充电模式的电动汽车充电功率预测方法,其特征在于:所述步骤1包括以下步骤:
步骤101:电动汽车的充电功率在时间和空间上均有较强的不确定性,针对电动汽车的时空不确定性进行分析是建立充电功率特性模型的基础,其中用户行为是影响电动汽车充电功率的关键,电动汽车的充电特性受电动汽车接入电网的时刻、需要充电电量和充电功率3方面因素的影响,本发明考虑私用调度汽车作为电网可调度资源,有多种充电方式,并以V2G的方式参与电网调度;
步骤102:电动汽车接入电网时刻与用户使用行为有密切关系,假设用户最后一次行驶返回时刻即为电动汽车的入网时刻,同时根据美国能源部2008年对电动汽车用户的调查统计报告,用户的最后一次返回时刻近似服从正态分布,最后一次返回时刻x的概率密度函数如下:
式中,TS表示最后一次返回时刻,μS和σS分别表示正态分布的期望和标准差;
步骤103:假设电动汽车一次充电即全部充满,则电动汽车当日的充电电量与其当日的行驶里程有密切关系,即电动汽车当日从电网获取的电能全部转化为当日行驶里程所消耗的能量,根据调查报告,电动汽车的日行驶里程数服从对数正态分布,电动汽车的日行驶里程数的概率密度函数如下:
式中:D表示日行驶里程数,μD和σD分别表示lnD的期望和标准差;
得到日行驶里程数后,结合电动汽车百公里电耗,得到需要充电电量如下:
式中:WC表示电动汽车需要充电电量,W100表示百公里耗电量;
步骤104:假设考虑电动汽车以常规慢充方式进行充电,由于起始充电阶段和结束充电阶段相较整个充电过程的时间较短,可以忽略,因此,假设电动汽车充电和放电均采用恒功率充放电方式进行,其充电功率为PC,放电功率为Pdisc;由于电动汽车模型是在峰谷电价的背景下建立的,电动汽车车主可以根据自身需要,采用不同的充电方式,以实现自身利益的最大化,即电动汽车对电价产生不同的响应方式,电动汽车的充电方式包括普通充电方式、延时充电方式和V2G充电方式3种。
3.根据权利要求1所述的计及延时和V2G充电模式的电动汽车充电功率预测方法,其特征在于:所述步骤3包括以下步骤:
步骤301:单台电动汽车接入系统进行充放电时,假设高峰时段的开始时刻为Tpeak1=9,结束时刻Tpeak2=22,在步骤2的基础上进行蒙特卡洛模拟,得到单台电动汽车的充放电负荷期望值以及不同充电方式下单台电动汽车各时段充放电负荷期望曲线;
步骤302:在步骤301的基础上,对模拟得到的N个充电功率样本采用非参数核密度估计,求得各充电方式下100台电动汽车在各时段的充放电集群负荷的概率分布函数,通过这些概率分布函数的曲线图可以对电动汽车充放电集群负荷的概率特性有清晰的认识。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710310372.8A CN107169273B (zh) | 2017-05-05 | 2017-05-05 | 计及延时和v2g充电模式的电动汽车充电功率预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710310372.8A CN107169273B (zh) | 2017-05-05 | 2017-05-05 | 计及延时和v2g充电模式的电动汽车充电功率预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107169273A CN107169273A (zh) | 2017-09-15 |
CN107169273B true CN107169273B (zh) | 2021-01-05 |
Family
ID=59813122
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710310372.8A Active CN107169273B (zh) | 2017-05-05 | 2017-05-05 | 计及延时和v2g充电模式的电动汽车充电功率预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107169273B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111486573A (zh) * | 2020-04-16 | 2020-08-04 | 南方电网科学研究院有限责任公司 | 一种中央空调集群调控方法、系统以及设备 |
Families Citing this family (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108400626A (zh) * | 2017-11-23 | 2018-08-14 | 蔚来汽车有限公司 | 电动型移动充电车的交直流充电方法 |
CN108229733A (zh) * | 2017-12-21 | 2018-06-29 | 中南大学 | 一种考虑用户满意度评价的电动汽车充电调度方法 |
CN108376985B (zh) * | 2018-02-02 | 2020-06-02 | 国家电网公司 | 一种电动汽车充电负荷模型建立及节点电压计算方法 |
CN109149561A (zh) * | 2018-08-13 | 2019-01-04 | 国网江苏省电力有限公司南京供电分公司 | 一种基于光储充电塔接入的配电网静态优化方法 |
CN109672199B (zh) * | 2018-11-29 | 2022-03-18 | 杭州电子科技大学 | 一种基于能量平衡的电动汽车削峰填谷能力估计方法 |
CN109624781B (zh) * | 2018-12-11 | 2021-02-12 | 华北电力大学 | 一种基于电动汽车能量响应能力预测的储能容量配置方法 |
CN109583136B (zh) * | 2018-12-28 | 2023-05-05 | 上海电力学院 | 基于可调度潜力的电动汽车充换储一体站模型建立方法 |
CN109866628B (zh) * | 2019-01-18 | 2022-07-15 | 国网上海市电力公司 | 一种主动配电网电动汽车有序充电控制方法 |
CN109886578A (zh) * | 2019-02-25 | 2019-06-14 | 电子科技大学 | 一种停车场内电动汽车充电调度方法 |
CN109910671B (zh) * | 2019-03-11 | 2019-10-08 | 三峡大学 | 基于智能充电桩的电动汽车v2g控制方法 |
CN110303931B (zh) * | 2019-05-29 | 2021-03-23 | 北京航盛新能科技有限公司 | 一种适用于v2g的电动汽车智能充电方法 |
CN110492509B (zh) * | 2019-08-19 | 2022-03-01 | 西南交通大学 | 一种基于电动汽车放电行为的分布式电源等效方法 |
CN110852494A (zh) * | 2019-10-25 | 2020-02-28 | 上海电力大学 | 一种计及需求响应的电动汽车概率负荷预测方法 |
CN112825102A (zh) * | 2019-11-20 | 2021-05-21 | 河海大学 | 一种电动汽车充电负荷特性模拟方法 |
CN111985777A (zh) * | 2020-07-20 | 2020-11-24 | 中国农业大学 | 一种电动汽车负荷聚合体调控能力评估模型建立方法及系统 |
CN114742118B (zh) * | 2020-12-23 | 2023-10-27 | 中国科学院广州能源研究所 | 电动汽车集群充放电负荷组合预测方法 |
CN113036753B (zh) * | 2021-02-04 | 2023-05-02 | 南昌工学院 | 一种模拟电动汽车用户群实现提高充电稳定性的叠加模型 |
CN113420953B (zh) * | 2021-05-24 | 2022-08-02 | 国网上海市电力公司电力科学研究院 | 一种基于容量可信度的柔性负荷可调节能力分析方法 |
CN114643892A (zh) * | 2022-04-11 | 2022-06-21 | 广州万城万充新能源科技有限公司 | 一种基于多模态数据感知的电动车充电功率预测系统 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8027812B2 (en) * | 2007-02-08 | 2011-09-27 | Kabushiki Kaisha Toshiba | Charge trajectory calculating method, system, and program |
CN102222928A (zh) * | 2011-06-16 | 2011-10-19 | 北京许继电气有限公司 | 电动汽车动力电池大型集中储能智能型充放电系统 |
CN104123598A (zh) * | 2014-08-07 | 2014-10-29 | 山东大学 | 一种基于多目标优化的电动私家车充电模式选择方法 |
CN106026076A (zh) * | 2016-05-17 | 2016-10-12 | 上海电力学院 | 一种计及电动汽车支撑能力的用户侧供电可靠性评估方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103679299B (zh) * | 2013-12-30 | 2016-09-28 | 华北电力大学(保定) | 兼顾车主满意度的电动汽车最优峰谷分时电价定价方法 |
CN106295860A (zh) * | 2016-07-29 | 2017-01-04 | 国网山东省电力公司经济技术研究院 | 一种基于蒙特卡洛模拟法的电动汽车规模化充电需求预测方法 |
-
2017
- 2017-05-05 CN CN201710310372.8A patent/CN107169273B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8027812B2 (en) * | 2007-02-08 | 2011-09-27 | Kabushiki Kaisha Toshiba | Charge trajectory calculating method, system, and program |
CN102222928A (zh) * | 2011-06-16 | 2011-10-19 | 北京许继电气有限公司 | 电动汽车动力电池大型集中储能智能型充放电系统 |
CN104123598A (zh) * | 2014-08-07 | 2014-10-29 | 山东大学 | 一种基于多目标优化的电动私家车充电模式选择方法 |
CN106026076A (zh) * | 2016-05-17 | 2016-10-12 | 上海电力学院 | 一种计及电动汽车支撑能力的用户侧供电可靠性评估方法 |
Non-Patent Citations (1)
Title |
---|
微网下V2G技术经济性的分析与研究;郭晓丽等;《电网与清洁能源》;20170331;第33卷(第3期);第43-48页、第55页 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111486573A (zh) * | 2020-04-16 | 2020-08-04 | 南方电网科学研究院有限责任公司 | 一种中央空调集群调控方法、系统以及设备 |
CN111486573B (zh) * | 2020-04-16 | 2021-09-14 | 南方电网科学研究院有限责任公司 | 一种中央空调集群调控方法、系统以及设备 |
Also Published As
Publication number | Publication date |
---|---|
CN107169273A (zh) | 2017-09-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107169273B (zh) | 计及延时和v2g充电模式的电动汽车充电功率预测方法 | |
Gan et al. | A probabilistic evaluation method of household EVs dispatching potential considering users’ multiple travel needs | |
CN103269107B (zh) | 一种电动汽车充换电站充换电控制方法 | |
CN103499792B (zh) | 电动汽车动力电池集群可用容量的预测方法 | |
CN111626527B (zh) | 计及可调度电动汽车快/慢充放电形式的智能电网深度学习调度方法 | |
CN103679299A (zh) | 兼顾车主满意度的电动汽车最优峰谷分时电价定价方法 | |
CN108596667B (zh) | 一种基于车联网的电动汽车实时充电电价计算方法 | |
CN109353242B (zh) | 一种智能充电桩系统实现双向有序充放电的充电算法 | |
CN111310966A (zh) | 含电动汽车充电站的微电网选址及优化配置方法 | |
CN109523087A (zh) | 基于深度学习的电动汽车快充站储能监测系统及其方法 | |
CN110232219B (zh) | 一种基于数据挖掘的电动汽车可调度容量核定方法 | |
Guner et al. | Impact of car arrival/departure patterns on EV parking lot energy storage capacity | |
CN110533222B (zh) | 基于峰平谷电价的电动汽车充电负荷预测方法及装置 | |
CN113060036A (zh) | 光储充电站的充电控制方法、装置、服务器及存储介质 | |
Guner et al. | Distributed storage capacity modelling of EV parking lots | |
CN108493959A (zh) | 基于Markov模型的电动汽车参与电网调频控制与分析方法 | |
CN102938094A (zh) | 参与调频服务的电动汽车充电时序确定与服务车辆选择方法 | |
CN113067370A (zh) | V2g充电站的充电控制方法、装置、服务器及存储介质 | |
Ma et al. | Real-time plug-in electric vehicles charging control for V2G frequency regulation | |
CN109672199B (zh) | 一种基于能量平衡的电动汽车削峰填谷能力估计方法 | |
CN109353243B (zh) | 一种智能充电桩系统实现双向有序充放电的放电算法 | |
CN117458477A (zh) | 考虑分群优化模式下负荷聚合商参与的电动汽车调度方法 | |
CN110861508B (zh) | 居民区直流充电机共享充电控制方法、系统及存储介质 | |
CN115330062A (zh) | 一种社区场景下的新能源汽车有序充电服务调度优化方法 | |
CN112109580B (zh) | 一种电量自分配的微电网电动汽车充放电控制系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |