CN110852494A - 一种计及需求响应的电动汽车概率负荷预测方法 - Google Patents
一种计及需求响应的电动汽车概率负荷预测方法 Download PDFInfo
- Publication number
- CN110852494A CN110852494A CN201911022933.XA CN201911022933A CN110852494A CN 110852494 A CN110852494 A CN 110852494A CN 201911022933 A CN201911022933 A CN 201911022933A CN 110852494 A CN110852494 A CN 110852494A
- Authority
- CN
- China
- Prior art keywords
- time
- vehicle
- load
- charging
- strategy
- 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.)
- Pending
Links
- 230000004044 response Effects 0.000 title claims abstract description 75
- 238000000034 method Methods 0.000 title claims abstract description 36
- 238000007599 discharging Methods 0.000 claims abstract description 34
- 238000009826 distribution Methods 0.000 claims abstract description 21
- 230000005611 electricity Effects 0.000 claims abstract description 21
- 230000008901 benefit Effects 0.000 claims abstract description 13
- 230000007246 mechanism Effects 0.000 claims description 11
- 230000005284 excitation Effects 0.000 claims description 10
- 230000000087 stabilizing effect Effects 0.000 claims description 4
- 238000013519 translation Methods 0.000 claims description 4
- 238000012614 Monte-Carlo sampling Methods 0.000 claims description 3
- 238000011156 evaluation Methods 0.000 claims description 3
- FFBHFFJDDLITSX-UHFFFAOYSA-N benzyl N-[2-hydroxy-4-(3-oxomorpholin-4-yl)phenyl]carbamate Chemical compound OC1=C(NC(=O)OCC2=CC=CC=C2)C=CC(=C1)N1CCOCC1=O FFBHFFJDDLITSX-UHFFFAOYSA-N 0.000 claims description 2
- 238000005457 optimization Methods 0.000 claims description 2
- 239000000126 substance Substances 0.000 claims description 2
- 238000012546 transfer Methods 0.000 claims description 2
- 230000006399 behavior Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 4
- 238000013507 mapping Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 230000002457 bidirectional effect Effects 0.000 description 2
- 230000000739 chaotic effect Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000003912 environmental pollution Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000009191 jumping Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 230000008569 process 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/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression 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/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
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Economics (AREA)
- Theoretical Computer Science (AREA)
- Strategic Management (AREA)
- General Business, Economics & Management (AREA)
- Marketing (AREA)
- Data Mining & Analysis (AREA)
- Operations Research (AREA)
- Tourism & Hospitality (AREA)
- Entrepreneurship & Innovation (AREA)
- Mathematical Analysis (AREA)
- Mathematical Physics (AREA)
- Quality & Reliability (AREA)
- Game Theory and Decision Science (AREA)
- Development Economics (AREA)
- Computational Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Mathematical Optimization (AREA)
- General Health & Medical Sciences (AREA)
- Bioinformatics & Computational Biology (AREA)
- Public Health (AREA)
- Probability & Statistics with Applications (AREA)
- Evolutionary Biology (AREA)
- Primary Health Care (AREA)
- Algebra (AREA)
- Water Supply & Treatment (AREA)
- Bioinformatics & Cheminformatics (AREA)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Life Sciences & Earth Sciences (AREA)
- Databases & Information Systems (AREA)
- Educational Administration (AREA)
- Charge And Discharge Circuits For Batteries Or The Like (AREA)
Abstract
本发明涉及一种计及需求响应的电动汽车概率负荷预测方法,包括以下步骤:1)判断当前电动汽车的接入电网持续时间是否大于充电至期望荷电状态的时间时,若是,则进行步骤2),若否,则当前车辆采用无序充电策略,并进行步骤6);2)以车辆收益最大化为目标构建基于分时电价的电动汽车自动充放电策略;3)对电动汽车自动充放电策略进行修正;4)计算用户参与响应的收益度;5)计算当前电动汽车对修正后电动汽车自动充放电策略的接受度;6)计算在各种可能情况下的电动汽车负荷及其期望值;7)计算各时刻下电动汽车负荷的概率分布。与现有技术相比,本发明具有响应策略合理、考虑需求响应不确定性、预测信息全面等优点。
Description
技术领域
本发明涉及电动汽车负荷预测领域,尤其是涉及一种计及需求响应的电动汽车概率负荷预测方法。
背景技术
近年来,环境污染和能源紧缺日益严重,电动汽车(electric vehicle,EV)产业凭借其低能耗、高性能和无污染等特点在国家政策扶持下蓬勃发展,未来电动汽车保有量将会持续增加。大规模电动汽车接入电网势必会增加配电网负担,引发负荷峰谷差加剧、线路过载、电能质量下降等负面影。为了降低这些负面影响,一方面需要通过电力市场的手段引导其进行有序充放电,从而减少这些负面影响的产生;另一方面需要对电动汽车负荷进行准确的预测,从而安排合理的调度计划。
在减少EV负面影响方面,电动汽车通过V2G(vehicle to grid)技术以需求响应的形式参与到电网“双向互动”中,合理安排充放电过程,从而提高系统的安全性和可靠性。需求响应主要可以分为基于激励的需求响应和基于价格的需求响应两类。现有研究多是基于单一的响应方式制定电动汽车有序充电策略,缺少对两种响应机制的综合考虑。
大多数用户对收益因素是敏感的,但是由于诸多非收益因素的影响,用户的响应行为具有很大的不确定性,针对这种不确定性,现有研究大多通过各种不确定性建模方法,最终抽取得到某一定性的结果,缺少对用户响应不确定性的定量刻画。云模型是一种描述不确定性问题的定性定量双向转换模型,已成功应用于电力综合评价、信号识别、乘客行为预测等方面。然而现有研究基本都采用正态云模型,这显然不符合用户收益与响应意愿间的模糊映射关系,因此需要建立一种新型的云模型以适应这种映射关系。
在电动汽车负荷预测方面,目前已有多种方法可对电动汽车负荷进行有效的预测,但是现有预测方法均默认车辆一接入电网就开始充电,没有考虑电动汽车参与需求响应的情况,得出的预测结果往往也是一种确定性的结果,这显然不适用于考虑电动汽车参与需求响应的负荷预测这一不确定性问题。
因此,急需一种考虑电动汽车参与需求响应的负荷预测方法,且能够较好地刻画电动汽车参与需求响应的不确定性,得出各种情况下的概率性负荷预测结果,对电网调度更具参考意义。
发明内容
本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种计及需求响应的电动汽车概率负荷预测方法。
本发明的目的可以通过以下技术方案来实现:
一种计及需求响应的电动汽车概率负荷预测方法,包括以下步骤:
1)判断当前电动汽车的接入电网持续时间是否大于充电至期望荷电状态的时间时,若是,则进行步骤2),若否,则当前车辆采用无序充电策略u0,i,并进行步骤6);
2)基于分时电价的响应机制,以车辆收益最大化为目标构建基于分时电价的电动汽车自动充放电策略u1,i;
3)结合基于激励的续期响应机制,考虑供电侧利益,以平抑负荷波动为目标,对步骤2)构建的电动汽车自动充放电策略进行修正,得到修正后的电动汽车自动充放电策略u2,i;
4)计算用户参与响应的收益度;
5)基于用户收益度,通过Sigmoid云模型计算当前电动汽车对修正后电动汽车自动充放电策略的接受度;
6)计算在各种可能情况下的电动汽车负荷及其期望值;
7)计算各时刻下电动汽车负荷的概率分布。
所述的步骤1)具体包括以下步骤:
11)根据电动汽车行驶特性的概率分布,通过蒙特卡洛抽样得到每辆车的接入电网时刻、离开电网时刻及日行驶里程,计算电动汽车的初始SOC:
其中,S0为电动汽车的初始SOC,Se为充电结束时荷电状态期望值,C为电池容量,Ed100为行驶100km所需的电能,Rd为日行驶里程;
12)判断当前车辆是否满足参与需求响应的条件,即接入电网持续时间是否大于充电至期望SOC的时间,则判断式如下:
Ts,i>Tc,i+Tm
其中,Tm为留有的时间裕度,Ts,i为车辆i的接入电网持续时间,Tc,i为充电至期望SOC的时间。
所述的步骤2)中,基于分时电价的电动汽车自动充放电策略以车辆收益最大化为目标建立目标函数,则有:
Et,i(Ii(t-1),Ii(t))=ε1(Ii(t))2+ε2(Ii(t)-Ii(t-1))2+ε3
其中,u1,i为基于分时电价的电动汽车自动充放电策略,其形式为[Ii(1),Ii(2),…Ii(t)…Ii(T)],Ts,i为车辆i接入电网的持续时间,Ii(t)为车辆i在t时段的充放电状态,其取值+1为充电,-1为放电,0为闲置,Rp,i(u1,i)表示在策略u1,i下,相比于无序充电所减少的费用,F0,i为无序充电情形下的充电费用,Pc,i(t)、Pd,i(t)分别为车辆i在t时段的充、放电功率,p(t)为t时段电价,Ei(Ii(t-1),Ii(t))为t-1时段到t时段电池损耗成本,ε1,ε2,ε3为电池损耗系数。
基于分时电价的电动汽车自动充放电策略的约束条件包括:
荷电状态约束:
Pc,i(t)Pd,i(t)=0
Smin≤Si(t)≤Smax
其中,Si(t)为车辆i在t时段的荷电状态SOC,ξc、ξd分别为充放电效率,Ci为车辆i的电池容量,Smax、Smin为荷电状态允许的最大值和最小值;
充电需求约束:
其中,S0,i为车辆i接入电网时的初始SOC,Se,i为车辆i离开电网时的期望SOC。
响应时间约束:
Ts,i>Tc,i
Tc,i=(Se,i-S0,i)Ci/Pc,iξc;
配电变压器容量约束:全天各个时段的总负荷不能大于变压器的最大负载;
L(t)=LB(t)+LEV(t)
其中,L(t)为t时段配电网总负荷,T为时段总数,LB(t)为t时段变压器接入的基线负荷,kT为变压器效率,AT为变压器额定容量,LEV(t)为t时段电动汽车集群负荷,Vt为当前时段已完成充放电计划制定的车辆集合。
所述的步骤3)中,所述的步骤3)中,首先采用模糊聚类的方法将当前负荷曲线分为峰、平、谷三类时段,分别用Фp,i及Фv,i表示峰、谷时段集合,则峰、谷时段的标识变量分别记为:
由此得出峰谷时段的响应情况为:
此式成立时,说明峰时段执行充电决策的时刻过多,加重了电网负担。
此式成立时,说明谷时段执行充电决策的时刻过少,未能充分响应。
其中,sum(R,r)表示矩阵R中等于r的元素个数,δp∈[0,1]、δv∈[0,1]分别表示电动汽车在峰、谷时段的响应裕度,当上述两式成立时,则需要给予一定的激励补偿;
修正后的电动汽车自动充放电策略以平抑负荷波动为优化目标,评价指标为负荷波动方差最小建立目标函数,则有:
minD(L(t))=E((L(t))2)-E2(L(t)),t=1,2…,T
F2,i=F1,i+RI,i(u2,i)-ΔRi
其中,D(L(t))为负荷曲线波动方差,E(*)表示数列*的期望,F1,i为经策略u1,i下用户的收益,F2,i为经策略u2,i下用户的收益,RI,i(u2,i)为在策略u2,i下激励机制补偿的费用,QI,i为激励过程所产生的负荷转移量,K1、K2分别为激励金额的二次项系数和一次项系数,ΔRi为从策略u1,i转变到策略u2,i所产生的电费损失和电池损耗成本。
所述的修正后的电动汽车自动充放电策略的约束条件包括:
荷电状态约束:
Pc,i(t)Pd,i(t)=0
Smin≤Si(t)≤Smax
其中,Si(t)为车辆i在t时段的荷电状态SOC,ξc、ξd分别为充放电效率,Ci为车辆i的电池容量,Smax、Smin为荷电状态允许的最大值和最小值;
充电需求约束:
其中,S0,i为车辆i接入电网时的初始SOC,Se,i为车辆i离开电网时的期望SOC;
收益提升约束:
RI,i(u2,i)-ΔRi≥0
充放电时间约束:
所述的步骤4)具体包括以下步骤:
41)计算车主实际收益Fr,i,则有:
其中,δti为实际提车时间与预期提车时间的差值,δFi为提前提车损失的收益,b0、b1分别为提前提车和延迟提车的惩罚系数,a0、a1分别为允许提前提车和延迟提车的差额裕度,(x)+=max{0,x};
42)以拒绝响应的充电成本为基准,定义用户参与响应的收益因子λi,则有:
其中,F0,i为在无序充电策略u0,i下的充电成本;
43)定义收益度ηi为用户参与响应所得收益与其心理预期的比值,则其表达式为:
其中,λE,i为用户心理预期的收益因子,并且该用户心理预期收益因子符合区间(0~1.5)内的均匀分布。
所述的步骤5)具体包括以下步骤:
51)通过逆向云发生器获得Sigmoid云模型的数字特征(Eα,Eβ,En,He);
52)通过X条件云发生器计算当前第i辆车在收益度为ηi的情况下对响应策略的接受度μi,具体为:
首先生成以En为期望值、He为标准差的正态随机数En’,并生成以Eα为期望值、En’为标准差的正态随机数α,然后生成以Eβ为期望值、En’为标准差的正态随机数β,最后由下式计算出云滴(ηi,μi),则有:
其中,ξc、ξd分别为充放电效率,Pc,i、Pd,i分别为车辆i在t时段的充、放电功率,分别表示策略u0,i和策略u2,i下处于放电状态的时段集合,分别表示策略u0,i和策略u2,i下处于充电状态的时段集合,Va,j为事件Xj中采用基于分时电价的电动汽车自动充放电策略的电动汽车集合,Vr,j为事件Xj中采用无序充电策略的电动汽车集合,表示事件Xj发生的概率。
所述的步骤7)中,在某一事件中t时刻的电动汽车负荷为:
其中,Vt +为所有在t时刻进行充电的EV集合,Vt -为所有在t时刻进行放电的电动汽车集合,通过计算所有事件的电动汽车负荷,从而得到t时刻电动汽车负荷的概率分布。
与现有技术相比,本发明具有以下优点:
一、响应策略合理:本发明提出的需求响应策略综合协调了基于价格和基于激励两种需求响应方式,使供需双方的综合利益最大。
二、有效考虑了需求响应不确定性:针对电动汽车用户响应行为间与收益因素之间的不确定关系,提出一种新型的Sigmoid云模型来定量刻画这种不确定性。
三、预测信息更为全面:相比于传统的电动汽车负荷预测方法,本发明最终可以得到各个时刻电动汽车负荷的概率信息以及不同置信度下的电动汽车负荷曲线,预测信息更加全面,对电网调度具有更好的参考意义。
附图说明
图1为考虑需求响应不确定性的EV日前负荷预测流程图。
图2为n辆车参与需求响应的事件组合及发生概率。
图3为所制定ADR策略与无序冲电负荷曲线对比图。
图4为四种不同事件下的EV负荷预测曲线,其中,图(4a)为在事件X1下的负荷预测曲线,图(4b)为在事件X2下的负荷预测曲线,图(4c)为在事件X3下的负荷预测曲线,图(4d)为在事件X4下的负荷预测曲线。
图5为四个不同时刻的EV负荷概率分布图,其中,图(5a)为在0时刻下的EV负荷概率分布图,图(5b)为在6时刻下的EV负荷概率分布图,图(5c)为在12时刻下的EV负荷概率分布图,图(5d)为在18时刻下的EV负荷概率分布图。
图6为不同置信度下的EV负荷预测曲线。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。
实施例
如图1所示,本发明提出了一种考虑需求响应不确定性的EV日前负荷预测方法,首先结合价格型和激励型两种需求响应机制设计了EV有序充放电的自动需求响应策略,然后为了准确描述用户响应该策略的随机性,提出一种新型的Sigmoid云模型,刻画了受复杂因素影响下用户收益度与响应行为之间的不确定映射关系,进而预测了EV参与需求响应的概率以及响应前后的负荷变化,最终计算出各个时刻EV负荷的概率信息以及不同置信度下的EV负荷曲线。
根据典型的电池充电特性可以近似认为每一个计算时段内EV充放电功率恒定,通过蒙特卡洛抽样得到每辆EV的接入电网时刻、离开电网时刻及日行驶里程,再通过以下公式计算其初始SOC:
其中,S0为初始SOC,Se为充电结束时荷电状态期望值;C为电池容量;Ed100为行驶100km所需的电能;Rd为日行驶里程。
假定在无序充电条件下,每辆EV从接入电网的时刻就开始充电,直至达到期望SOC后停止充电,用u0,i表示,从而计算得到无序充电情况下EV各时段的充电功率,然后基于分时电价和激励补偿机制制定需求响应策略,再通过Sigmoid云模型判断对该策略的接受度,最后计算考虑需求响应后的EV负荷情况,具体步骤如下:
步骤1:判断第i辆EV是否满足响应条件,即接入电网持续时间是否大于充电至期望SOC的时间,表达式如下:
Ts,i>Tc,i+Tm
其中,Tm表示留有的时间裕度,若此式不成立,则视当前EV无法参与响应,采用无序充电策略,令u2,i=u0,i,μi=0,跳转步骤5。
步骤2:基于分时电价制定ADR策略。以用户侧利益最大化为目标设计基于分时电价的最优充放电策略u1,i,详见发明内容步骤2)。
步骤3:基于激励补偿机制修正ADR策略。在策略u1,i的基础上考虑供电侧利益,以平抑负荷波动为目标,对u1,i进行调整,得到u2,i,详见发明内容步骤3)。
步骤4:计算当前EV对所制定ADR策略的接受度。首先根据发明内容步骤4)计算用户通过策略u2,i可以获得的收益度ηi,再通过Sigmoid云模型得出用户对该策略的接受度μi。
步骤5:若所有车辆是否全部计算完成,进行步骤6,否则,重复步骤1-4。
步骤6:计算EV负荷。每辆EV都有两种选择,参与响应接受ADR策略和不参与响应进行无序充电,因此对于n辆车来说,就会产生2n种事件组合,如图2所示。算例以下列四个事件为例:
X1:所有车辆全部参与响应,概率为
每个事件Xj下的EV负荷为:
日前EV总负荷期望值为:
步骤7:计算t时刻EV负荷的概率分布。某一事件中t时刻的EV负荷为:
其中,Vt +表示所有在t时刻进行充电的EV集合,Vt -表示所有在t时刻进行放电的EV集合。计算所有事件的从而得到t时刻EV负荷的概率分布,0时刻、6时刻、12时刻、18时刻的负荷分布如图5所示。
步骤8:计算不同置信度Q下的EV负荷:
其中,Vμ>Q表示所有接受度大于Q的EV集合,Vμ<Q表示所有接受度小于Q的EV集合的EV集合,以Q取0.5、0.6、0.8、0.95为例,计算结果如图6所示。
本方法结合价格和激励两种需求响应机制制定了一种需求响应策略,在此基础上考虑了电动汽车用户接受这种策略的不确定性,进而对区域内日前电动汽车负荷进行预测。如图3所示,制定的需求响应策略综合考虑了用户侧和供电侧双方的利益,在减少用户充电成本的前提下能够有效平抑配电网负荷波动。在此基础上,通过Sigmoid云模型准确描述了用户收益与响应行为之间的不确定关系,对EV负荷进行了概率性预测,如图4-6所示,可以得到各种情况下EV负荷的概率性结果,预测信息更为全面,预测结果也更符合客观情况,对电网调度具有更好的参考意义。
Claims (10)
1.一种计及需求响应的电动汽车概率负荷预测方法,其特征在于,包括以下步骤:
1)判断当前电动汽车的接入电网持续时间是否大于充电至期望荷电状态的时间时,若是,则进行步骤2),若否,则当前车辆采用无序充电策略u0,i,并进行步骤6);
2)基于分时电价的响应机制,以车辆收益最大化为目标构建基于分时电价的电动汽车自动充放电策略u1,i;
3)结合基于激励的续期响应机制,考虑供电侧利益,以平抑负荷波动为目标,对步骤2)构建的电动汽车自动充放电策略进行修正,得到修正后的电动汽车自动充放电策略u2,i;
4)计算用户参与响应的收益度;
5)基于用户收益度,通过Sigmoid云模型计算当前电动汽车对修正后电动汽车自动充放电策略的接受度;
6)计算在各种可能情况下的电动汽车负荷及其期望值;
7)计算各时刻下电动汽车负荷的概率分布。
2.根据权利要求1所述的一种计及需求响应的电动汽车概率负荷预测方法,其特征在于,所述的步骤1)具体包括以下步骤:
11)根据电动汽车行驶特性的概率分布,通过蒙特卡洛抽样得到每辆车的接入电网时刻、离开电网时刻及日行驶里程,计算电动汽车的初始SOC:
其中,S0为电动汽车的初始SOC,Se为充电结束时荷电状态期望值,C为电池容量,Ed100为行驶100km所需的电能,Rd为日行驶里程;
12)判断当前车辆是否满足参与需求响应的条件,即接入电网持续时间是否大于充电至期望SOC的时间,则判断式如下:
Ts,i>Tc,i+Tm
其中,Tm为留有的时间裕度,Ts,i为车辆i的接入电网持续时间,Tc,i为充电至期望SOC的时间。
3.根据权利要求1所述的一种计及需求响应的电动汽车概率负荷预测方法,其特征在于,所述的步骤2)中,基于分时电价的电动汽车自动充放电策略以车辆收益最大化为目标建立目标函数,则有:
Et,i(Ii(t-1),Ii(t))=ε1(Ii(t))2+ε2(Ii(t)-Ii(t-1))2+ε3
其中,u1,i为基于分时电价的电动汽车自动充放电策略,其形式为[Ii(1),Ii(2),…Ii(t)…Ii(T)],Ts,i为车辆i接入电网的持续时间,Ii(t)为车辆i在t时段的充放电状态,其取值+1为充电,-1为放电,0为闲置,Rp,i(u1,i)表示在策略u1,i下,相比于无序充电所减少的费用,F0,i为无序充电情形下的充电费用,Pc,i(t)、Pd,i(t)分别为车辆i在t时段的充、放电功率,p(t)为t时段电价,Ei(Ii(t-1),Ii(t))为t-1时段到t时段电池损耗成本,ε1,ε2,ε3为电池损耗系数。
4.根据权利要求3所述的一种计及需求响应的电动汽车概率负荷预测方法,其特征在于,基于分时电价的电动汽车自动充放电策略的约束条件包括:
荷电状态约束:
Pc,i(t)Pd,i(t)=0
Smin≤Si(t)≤Smax
其中,Si(t)为车辆i在t时段的荷电状态SOC,ξc、ξd分别为充放电效率,Ci为车辆i的电池容量,Smax、Smin为荷电状态允许的最大值和最小值;
充电需求约束:
其中,S0,i为车辆i接入电网时的初始SOC,Se,i为车辆i离开电网时的期望SOC。
响应时间约束:
Ts,i>Tc,i
Tc,i=(Se,i-S0,i)Ci/Pc,iξc;
配电变压器容量约束:全天各个时段的总负荷不能大于变压器的最大负载;
L(t)=LB(t)+LEV(t)
其中,L(t)为t时段配电网总负荷,T为时段总数,LB(t)为t时段变压器接入的基线负荷,kT为变压器效率,AT为变压器额定容量,LEV(t)为t时段电动汽车集群负荷,Vt为当前时段已完成充放电计划制定的车辆集合。
5.根据权利要求1所述的一种计及需求响应的电动汽车概率负荷预测方法,其特征在于,所述的步骤3)中,修正后的电动汽车自动充放电策略以平抑负荷波动为优化目标,评价指标为负荷波动方差最小建立目标函数,则有:
minD(L(t))=E((L(t))2)-E2(L(t)),t=1,2…,T
F2,i=F1,i+RI,i(u2,i)-ΔRi
其中,D(L(t))为负荷曲线波动方差,E(*)表示数列*的期望,F1,i为经策略u1,i下用户的收益,F2,i为经策略u2,i下用户的收益,RI,i(u2,i)为在策略u2,i下激励机制补偿的费用,QI,i为激励过程所产生的负荷转移量,K1、K2分别为激励金额的二次项系数和一次项系数,ΔRi为从策略u1,i转变到策略u2,i所产生的电费损失和电池损耗成本。
6.根据权利要求5所述的一种计及需求响应的电动汽车概率负荷预测方法,其特征在于,所述的修正后的电动汽车自动充放电策略的约束条件包括:
荷电状态约束:
Pc,i(t)Pd,i(t)=0
Smin≤Si(t)≤Smax
其中,Si(t)为车辆i在t时段的荷电状态SOC,ξc、ξd分别为充放电效率,Ci为车辆i的电池容量,Smax、Smin为荷电状态允许的最大值和最小值;
充电需求约束:
其中,S0,i为车辆i接入电网时的初始SOC,Se,i为车辆i离开电网时的期望SOC;
收益提升约束:
RI,i(u2,i)-ΔRi≥0
充放电时间约束:
7.根据权利要求5所述的一种计及需求响应的电动汽车概率负荷预测方法,其特征在于,所述的步骤4)具体包括以下步骤:
41)计算车主实际收益Fr,i,则有:
其中,δti为实际提车时间与预期提车时间的差值,δFi为提前提车损失的收益,b0、b1分别为提前提车和延迟提车的惩罚系数,a0、a1分别为允许提前提车和延迟提车的差额裕度,(x)+=max{0,x};
42)以拒绝响应的充电成本为基准,定义用户参与响应的收益因子λi,则有:
其中,F0,i为在无序充电策略u0,i下的充电成本;
43)定义收益度ηi为用户参与响应所得收益与其心理预期的比值,则其表达式为:
其中,λE,i为用户心理预期的收益因子,并且该用户心理预期收益因子符合区间(0~1.5)内的均匀分布。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911022933.XA CN110852494A (zh) | 2019-10-25 | 2019-10-25 | 一种计及需求响应的电动汽车概率负荷预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911022933.XA CN110852494A (zh) | 2019-10-25 | 2019-10-25 | 一种计及需求响应的电动汽车概率负荷预测方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110852494A true CN110852494A (zh) | 2020-02-28 |
Family
ID=69596929
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911022933.XA Pending CN110852494A (zh) | 2019-10-25 | 2019-10-25 | 一种计及需求响应的电动汽车概率负荷预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110852494A (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111682521A (zh) * | 2020-05-08 | 2020-09-18 | 天津大学 | 物流中心多台ev参与电网需求响应的充放电控制方法 |
CN112332433A (zh) * | 2020-09-29 | 2021-02-05 | 国网电动汽车服务有限公司 | 一种电动汽车参与填谷辅助服务可转移负荷容量分析方法 |
CN113361892A (zh) * | 2021-05-31 | 2021-09-07 | 东南大学 | 基于时序云模型的电动汽车需求响应潜力估计方法 |
CN113949090A (zh) * | 2021-11-05 | 2022-01-18 | 国网江苏省电力有限公司无锡供电分公司 | 一种基于电动汽车集群的实时需求响应互动方法 |
CN114742118A (zh) * | 2020-12-23 | 2022-07-12 | 中国科学院广州能源研究所 | 电动汽车集群充放电负荷组合预测方法 |
CN114781061A (zh) * | 2022-04-20 | 2022-07-22 | 国网江苏省电力有限公司电力科学研究院 | 一种电动汽车集群响应能力评估方法及装置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20160126732A1 (en) * | 2014-11-05 | 2016-05-05 | Honda Motor Co., Ltd. | System and method for aggregating electric vehicle loads for demand response events |
CN105631553A (zh) * | 2016-02-01 | 2016-06-01 | 浙江工业大学 | 一种电动汽车充放电自动需求响应优化方法 |
CN107169273A (zh) * | 2017-05-05 | 2017-09-15 | 河海大学 | 计及延时和v2g充电模式的电动汽车充电功率预测方法 |
-
2019
- 2019-10-25 CN CN201911022933.XA patent/CN110852494A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20160126732A1 (en) * | 2014-11-05 | 2016-05-05 | Honda Motor Co., Ltd. | System and method for aggregating electric vehicle loads for demand response events |
CN105631553A (zh) * | 2016-02-01 | 2016-06-01 | 浙江工业大学 | 一种电动汽车充放电自动需求响应优化方法 |
CN107169273A (zh) * | 2017-05-05 | 2017-09-15 | 河海大学 | 计及延时和v2g充电模式的电动汽车充电功率预测方法 |
Non-Patent Citations (1)
Title |
---|
郭敏: "一种基于云模型的负荷频率控制方法", 《电子设计工程》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111682521A (zh) * | 2020-05-08 | 2020-09-18 | 天津大学 | 物流中心多台ev参与电网需求响应的充放电控制方法 |
CN111682521B (zh) * | 2020-05-08 | 2023-03-14 | 天津大学 | 物流中心多台ev参与电网需求响应的充放电控制方法 |
CN112332433A (zh) * | 2020-09-29 | 2021-02-05 | 国网电动汽车服务有限公司 | 一种电动汽车参与填谷辅助服务可转移负荷容量分析方法 |
CN114742118A (zh) * | 2020-12-23 | 2022-07-12 | 中国科学院广州能源研究所 | 电动汽车集群充放电负荷组合预测方法 |
CN114742118B (zh) * | 2020-12-23 | 2023-10-27 | 中国科学院广州能源研究所 | 电动汽车集群充放电负荷组合预测方法 |
CN113361892A (zh) * | 2021-05-31 | 2021-09-07 | 东南大学 | 基于时序云模型的电动汽车需求响应潜力估计方法 |
CN113949090A (zh) * | 2021-11-05 | 2022-01-18 | 国网江苏省电力有限公司无锡供电分公司 | 一种基于电动汽车集群的实时需求响应互动方法 |
CN113949090B (zh) * | 2021-11-05 | 2024-04-26 | 国网江苏省电力有限公司无锡供电分公司 | 一种基于电动汽车集群的实时需求响应互动方法 |
CN114781061A (zh) * | 2022-04-20 | 2022-07-22 | 国网江苏省电力有限公司电力科学研究院 | 一种电动汽车集群响应能力评估方法及装置 |
CN114781061B (zh) * | 2022-04-20 | 2024-01-23 | 国网江苏省电力有限公司电力科学研究院 | 一种电动汽车集群响应能力评估方法及装置 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110852494A (zh) | 一种计及需求响应的电动汽车概率负荷预测方法 | |
CN112134300B (zh) | 基于预约的电动汽车光储充电站滚动优化运行方法及系统 | |
CN105024432B (zh) | 一种基于虚拟电价的电动汽车充放电优化调度方法 | |
CN111431198B (zh) | 峰谷电价下考虑电动汽车v2g能力的配电网储能调峰方法 | |
CN103679299A (zh) | 兼顾车主满意度的电动汽车最优峰谷分时电价定价方法 | |
CN109849723B (zh) | 一种基于充电站功率裕度的电动汽车有序充电控制方法 | |
CN110406422B (zh) | 一种计及多方利益的电动公交车电池参与v2g控制方法 | |
CN111845453B (zh) | 考虑柔性控制的电动汽车充电站双层优化充放电策略 | |
CN109193718A (zh) | 一种适应于v2g的选择电动汽车入网调控方法 | |
CN113269372A (zh) | 一种考虑用户意愿的集群电动汽车可调度容量预测评估方法 | |
CN111738518B (zh) | 一种基于计及平均放电率的电动汽车充放电调度方法 | |
CN111681127B (zh) | 一种居民小区电动汽车有序充放电控制方法 | |
CN110390468B (zh) | 电动汽车充电的调度方法和系统 | |
CN112590601A (zh) | 基于边缘计算平台的新型v2g充电场站系统及其充电方法 | |
CN111619394A (zh) | 一种基于分时电价的电动公交车有序充放电方法 | |
CN113829934A (zh) | 一种电动汽车集群聚合响应能力确定方法及调度方法 | |
CN112348591A (zh) | 一种电动汽车的有序充电控制方法和装置 | |
CN114784838A (zh) | 计及供需两侧的充放电调度方法 | |
CN112332433B (zh) | 一种电动汽车参与填谷辅助服务可转移负荷容量分析方法 | |
CN117220281A (zh) | 一种电动汽车接入电网调节能力量化评估方法及系统 | |
CN110334903B (zh) | 基于背包算法的电动汽车充电调度方法 | |
CN107016504A (zh) | 一种计及电动汽车随机模糊需求响应的调度优化建模及算法 | |
CN115983441A (zh) | 一种电动汽车集群的灵活性概率预测与随机优化调控方法 | |
CN115626072A (zh) | 基于用户间博弈的网联电动汽车协同充放电调控方法 | |
CN115330062A (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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20200228 |
|
RJ01 | Rejection of invention patent application after publication |