CN105576684B - 一种含高渗透率光电微电网中的电动汽车优化调度方法 - Google Patents
一种含高渗透率光电微电网中的电动汽车优化调度方法 Download PDFInfo
- Publication number
- CN105576684B CN105576684B CN201610071449.6A CN201610071449A CN105576684B CN 105576684 B CN105576684 B CN 105576684B CN 201610071449 A CN201610071449 A CN 201610071449A CN 105576684 B CN105576684 B CN 105576684B
- Authority
- CN
- China
- Prior art keywords
- game
- microgrid
- charging
- electric vehicle
- grid
- 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
- 230000035699 permeability Effects 0.000 title claims abstract description 15
- 238000005457 optimization Methods 0.000 title claims abstract description 10
- 230000005622 photoelectricity Effects 0.000 title abstract description 5
- 238000004146 energy storage Methods 0.000 claims abstract description 74
- 238000005286 illumination Methods 0.000 claims abstract description 6
- 238000007599 discharging Methods 0.000 claims description 72
- 230000008901 benefit Effects 0.000 claims description 40
- 238000010248 power generation Methods 0.000 claims description 37
- 238000000034 method Methods 0.000 claims description 28
- 230000005611 electricity Effects 0.000 claims description 12
- 238000007726 management method Methods 0.000 claims description 7
- 230000008569 process Effects 0.000 claims description 5
- 238000013439 planning Methods 0.000 claims description 4
- 239000000126 substance Substances 0.000 claims description 4
- 230000005540 biological transmission Effects 0.000 claims description 3
- 230000005251 gamma ray Effects 0.000 claims description 3
- 230000003993 interaction Effects 0.000 claims description 3
- 230000000694 effects Effects 0.000 claims description 2
- 150000001875 compounds Chemical class 0.000 claims 1
- 238000011160 research Methods 0.000 description 3
- 230000002195 synergetic effect Effects 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000004134 energy conservation Methods 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 238000012384 transportation and delivery Methods 0.000 description 2
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 239000003245 coal Substances 0.000 description 1
- 238000011217 control strategy Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 239000000295 fuel oil Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000000087 stabilizing effect Effects 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/28—Arrangements for balancing of the load in a network by storage of energy
- H02J3/32—Arrangements for balancing of the load in a network by storage of energy using batteries with converting means
-
- 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
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/50—Photovoltaic [PV] energy
- Y02E10/56—Power conversion systems, e.g. maximum power point trackers
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E70/00—Other energy conversion or management systems reducing GHG emissions
- Y02E70/30—Systems combining energy storage with energy generation of non-fossil origin
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Economics (AREA)
- Marketing (AREA)
- Tourism & Hospitality (AREA)
- General Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Public Health (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Water Supply & Treatment (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Power Engineering (AREA)
- Charge And Discharge Circuits For Batteries Or The Like (AREA)
Abstract
一种含高渗透率光电微电网中的电动汽车优化调度方法,包括以下步骤:S1:将全天24h均分为J个时段,对于任意第k时段,有k∈{1,2,...,J},且第k时段的时长为Δt,获取光伏发电的光照强度和环境温度信息、储能系统的蓄电池容量信息以及微电网的常规负荷信息;S2:记录当前接入微电网的电动汽车的电池信息和客户充电需求信息;S3:根据步骤S1、S2获取的信息基于虚拟费用构建电动汽车集群的NCD博弈模型;S4:求解所构建的NCD博弈模型,以最大化利用光伏发电为目的,获取电动汽车的最优充放电计划,以实现微电网内电动汽车的优化调度。
Description
技术领域
本发明涉及一种含高渗透率光电微电网中的电动汽车优化调度方法。
背景技术
电力行业是节能减排的重点领域,作为节能减排的最佳途径之一,电动汽车(electric vehicle,EV)产业获得了迅速的发展,逐渐成为智能电网领域不可或缺的重要组成部分。电动汽车接入电网(Vehicle to Grid,V2G)概念的提出,使得大规模电动汽车入网将对电力系统的规划、运行以及电力市场的运营产生深刻影响。然而,在发电侧一次能源仍以煤炭为主的情形下,EV直接接入电网充电,会带来大量的间接碳排放,相比传统燃油汽车,难以体现其节能减排优势。现阶段,以光电为代表的新能源的战略地位开始由目前的补充能源上升至替代能源,装机容量迅速提升。
从当前发展情况来看,通过调整电网侧一次能源结构的途径实现对大规模光电的消纳会大幅增加电网改造成本,而以微电网为平台实现EV对光电的集成利用,将成为未来大规模可再生能源并网的主要方式之一,相对于传统电网,其主要优势体现在:实现高渗透率的光电接入,具有较好的经济性和环保性。EV与光伏发电的本地化协调控制更有利于EV备用功能的发挥以及系统电能质量和稳定性的提高。
光伏发电的出力特性导致其规模化利用将给系统带来巨大挑战,为了满足系统的安全、经济运行以及供电电能质量要求,已出现严重的弃光现象。含高渗透率光电并网型微电网中,有效的调度和控制为手段充分发挥EV集群的储能作用成为最大限度提高光伏发电利用率、体现EV与大规模光电协同增效利用的关键。文献“含电动汽车换电站的微电网孤岛运行优化”,电力自动化设备,2012,32(5)的研究表明,EV换电站作为储能可以提高微电网的可再生能源接纳能力和可靠性,并更具经济性,但较少涉及含高渗透率光电并网型微电网环境下EV与光伏的协同增效模型;文献“Economic impacts of small-scale owngenerating and storage units,and electric vehicles under different demandresponse strategies for smart households”,Applied Energy,126(2014)针对智能家居,以需求响应视角进行了EV、光伏和储能集成利用的经济性分析,为EV和可再生能源的协同利用研究提供了借鉴,但属于集中式优化,优化效率仍存在提升空间;文献“Adecentralized charging control strategy for plug-in electric vehicles tomitigate wind farm intermittency and enhance frequency regulation”,Journal ofPower Sources,248(2014)提出EV的分散式控制方法,并据此侧重研究了EV的平抑风电出力波动以及调频作用,其并未着眼于提高光伏并网效率。
发明内容
本发明提供一种含高渗透率光电微电网中的电动汽车优化调度方法,本发明针对含EV集群、高渗透率光电的并网型微电网,研究了EV与大规模光电协同增效利用的方法,提出的电动汽车优化调度方法,更加符合光伏发电情况,从而提升光电利用率、提高微电网的经济效益;基于光伏出力信息的虚拟费用理论,在此基础上构造了电动汽车集群的非合作充放电(Non-cooperative Charging and Discharging,NCD)博弈模型,推导出博弈模型纳什均衡与集中式优化模型最优解的一致性。
本发明的技术方案是:
一种含高渗透率光电微电网中的电动汽车优化调度方法,所述的微电网由高渗透率光伏发电系统、储能系统、常规负荷和电动汽车集群构成,常规负荷和电动汽车集群的充电负荷构成微电网总负荷,电动汽车集群由若干电动汽车构成,电动汽车既可作为充电负荷又可作为移动储能单元向微电网供电,微电网内的高渗透率的光伏发电系统、储能系统和作为移动储能单元时的电动汽车优先给微电网供电,富余的电量送到地区大电网,当微电网供电不足时,再从地区大电网购电,以满足微电网总负荷的需求;包括以下步骤:
S1:将全天24h均分为J个时段,对于任意第k时段,有k∈{1,2,...,J},且第k时段的时长为Δt,获取光伏发电的光照强度和环境温度信息、储能系统的蓄电池容量信息以及微电网的常规负荷信息;
S2:记录当前接入微电网的电动汽车的电池信息和客户充电需求信息;
S3:根据步骤S1、S2获取的信息基于虚拟费用构建电动汽车集群的NCD博弈模型;
S4:求解所构建的NCD博弈模型,以最大化利用光伏发电为目的,获取电动汽车的最优充放电计划,以实现微电网内电动汽车的优化调度。
进一步,本发明步骤S2中用一个七维行向量记录电动汽车的电池信息和客户充电需求信息,假设接入微电网的电动汽车的总数为N,对于任一电动汽车l∈{1,2,...,N},则有:
Xl=[Tin,l,Tout,l,S0,l,SE,l,Qs,l,Pc,l,Pd,l] (1)
式中:
Tin,l、Tout,l分别表示电动汽车l接入微电网的时间和预期离开微电网的时间;
S0,l、SE,l分别表示电动汽车l的电池的起始SOC(State of Charge,荷电状态,电池剩余能量与电池容量的比值)和离开微电网时客户期望电池达到的SOC,0≤S0,l≤1、0≤SE,l≤1;
Qs,l表示电动汽车l的电池容量;
Pc,l、Pd,l分别表示电动汽车l的额定充、放电功率。
进一步,本发明所述步骤S3具体包括以下步骤:
S31.以所有接入微电网的电动汽车和微电网内储能系统的蓄电池为博弈参与者,即决策主体U={1,2,...,N,N+1},构建任一博弈参与者i∈U的充放电计划模型;
对于博弈参与者之一的电动汽车l,其在全天J个时段内的充放电计划可表示如下:
Pl=(Pl(1),…,Pl(k),…,Pl(J)) (2)
其中,Pl为电动汽车l的充放电计划;Pl(k)表示第k时段,电动汽车l的充放电功率;
视电动汽车l的电池以恒功率充放电,忽略电池自身放电的影响,建立电动汽车l的电池模型和约束条件分别如下:
Sl(k)=Sl(k-1)+Pl(k)η(Pl(k))Δt/Qs,l (3)
SEV,min≤Sl(k)≤SEV,max (4)
-Pd,l≤Pl(k)≤Pc,l (5)
其中:
式(3)是电动汽车l的电池模型,Sl(k-1)、Sl(k)分别表示电动汽车l的电池在第k-1个和第k个时段的SOC,η(Pl(k))表示第k时段的电动汽车l的充放电效率;
SEV,max、SEV,min分别为电动汽车l的电池的SOC的最大值和最小值;
式(4)、(5)分别为电动汽车l的电池的SOC的约束和充放电功率Pl(k)的约束;
式(6)表示客户充电需求的约束,即电动汽车l如约离开时,电动汽车l的电池的SOC需满足客户期望;
微电网内储能系统的蓄电池是另一个博弈参与者,可视为一类无充电负荷需求的特殊“电动汽车”,其充放电计划和约束条件可类比于式(2)、式(4)~(5);此外,储能系统的蓄电池需满足的等式约束为:
即全天J个时段内储能系统的蓄电池的总充、放电电量需相等以满足充放电的循环,式中分别为储能系统的蓄电池的充、放电功率;ηBc、ηBd分别表示储能系统的蓄电池的充、放电效率;
S32.生成所有博弈参与者的策略空间,即生成所有博弈参与者的充放电计划集合[P1,P2,…,Pi,…,PN,PN+1],其中Pi=(Pi(1),Pi(2),...Pi(k),...,Pi(J))表示博弈参与者i在J个时段内充放电计划的集合,Pi(k)表示博弈参与者i在第k个时段的充放功率;在生成策略空间的过程中,除了考虑步骤S31中如式(4)~(7)所述的博弈参与者的运行约束外,还应考虑微电网与地区大电网之间的关系约束,主要包括:
Pgridout≤Pgridout,max (9)
式(8)表示微电网与地区大电网之间的功率传输平衡约束,其中,PPV(k)表示第k时段,高渗透率光伏发电系统的光伏发电输出功率,Pgrid(k)表示第k时段微电网与地区大电网的交互功率,Pgrid(k)>0表示微电网从地区大电网购电,Pgrid(k)<0表示微电网向地区大电网倒送输电,LB(k)表示第k时段微电网的常规负荷;
式(9)表示微电网向地区大电网倒送输电的倒送功率约束,其中,Pgridout表示微电网向地区大电网倒送输电的倒送功率,Pgridout,max为倒送功率允许的最大值;
综上,对于任一博弈参与者i的策略空间ρi可表示为:
ρi={Pi|constraint(2),(4),(5),(6),(7),(8),(9)} (10)
式(10)表示博弈参与者i的策略空间ρi包括博弈参与者i在J个时段内的充放电计划集合以及所要满足的约束(Constraint)条件:式(2)、式(4)~式(9);
S33.电动汽车和储能系统的蓄电池接入微电网后,计算博弈参与者i在J个时段内获得的虚拟利益;首先计算所有博弈参与者i的总虚拟利益:
式中:
γa表示总虚拟利益;γi表示博弈参与者i的虚拟利益;
γ(LB(k),k)=σ(LB(k)-PPV(k))2表示第k时段的基本虚拟费用;γ(La(k),k)=σ(La(k)-PPV(k))2表示第k时段的总虚拟费用;其中,表示第k时段微电网总负荷;
σ表示虚拟系数,且σ>0;
与实际电费不同,虚拟费用仅作为规划博弈参与者充放电计划的依据,用以制定各博弈参与者的最佳充放电计划;对于任一博弈参与者i的虚拟利益γi与博弈参与者i的充电需求以及负荷转移能力密切相关,假设博弈参与者i的负荷转移能力与其电池容量正相关,则博弈参与者i分享总虚拟利益的比例为:
式中,εi为博弈参与者i分享总虚拟利益的比例;Eneed,i表示博弈参与者i的充电需求量;κ表示正相关系数;特别地,由于接入微电网的电动汽车总数为N,因此当i∈{1,2,...,N}时,博弈参与者i指电动汽车,Eneed,i≠0;当i=N+1时博弈参与者i指储能系统的蓄电池,且有Eneed,i=0,Qs,i表示储能系统的蓄电池的配置容量;
根据式(11)和式(12)便可计算博弈参与者i的虚拟利益γi:
γi=εiγa (13)
S34.建立用于度量博弈参与者i参与博弈后的成本或效益的效用函数;定义博弈参与者i的虚拟利益γi为博弈参与者的效用函数,则博弈参与者i在充放电计划集合[P1,P2,…,Pi,…,PN,PN+1]下的效用函数Fi(Pi,Pi -)表示为:
式中,Pi -=[P1,P2,…,Pi-1,Pi+1,…,PN,PN+1],表示除博弈参与者i之外的其他所有博弈参与者的充放电计划集合;表示第k时段微电网内除博弈参与者i的负荷之外的其他负荷之和,其中τ表示除博弈参与者i之外的任一博弈参与者;
S35.构建决策主体的NCD博弈模型,博弈模型中,所有博弈参与者都是独立的决策者,各博弈参与者通过观察微电网的常规负荷、光伏发电输出功率以及其他博弈参与者的充放电计划来确定自身的最优充放电计划,以最大化虚拟利益;根据步骤S31~S34建立NCD博弈模型如下:
GNCD={U,{ρi}i∈U,{Fi}i∈U} (15)
基于上述NCD博弈模型,所有博弈参与者均会在假定其余博弈参与者都已选最优充放电计划的前提下选择能使其虚拟利益最大的最优充放电计划,即
其中,Pi *为博弈参与者i的最优充放电计划;
在NCD博弈模型中,如果对于博弈参与者充放电计划的其中一个集合满足如下条件:
Fi(Pi *,Pi *-)≥Fi(Pi,Pi *-),Pi∈ρi,i∈U (17)
则称P*为该NCD博弈模型的纳什均衡,即P*为各博弈参与者的最优充放电计划集合;其中Pi *-表示除博弈参与者i之外的其他所有博弈参与者的最优充放电计划集合。
进一步,本发明所述步骤S4中,基于NCD博弈模型的电动汽车的最优充放电计划求解过程具体如下:
(1)初始化电动汽车和储能系统的蓄电池的参数、第k时段配电网常规负荷LB(k)信息和光伏发电输出功率PPV(k);
(2)能量管理系统实时更新该时段的微电网总负荷La(k),并分享给所有博弈参与者;
(3)对于任一博弈参与者i∈U,其收到微电网总负荷La(k)信息后,计算出并以最大化式(14)为目标计算得到Pi *;
(4)如果计算得到的Pi *与当前最优值比较发生变化,更新Pi *,并将Pi *反馈给能量管理系统,转至(2);
(5)重复步骤(2)~(4),直至结果收敛(即Pi *不再发生变化),输出最优充放电计划Pi *;
为提高NCD博弈模型收敛至纳什均衡点的速度,步骤(4)中的各参与者可采用异步更新方式,即所有参与者不在同一时刻更新其充放电计划。
本发明的有益效果是:
1)本发明不依托于能量管理系统的集中式优化,基于NCD博弈模型发展EV及蓄电池的分散自治,大大降低决策难度,适于大规模、分散式EV的优化调度。
2)采用本发明,能够在满足EV用户的用电需求前提下,有效改善负荷特性,使其更加符合光伏发电情况,从而提升光电利用率、提高微电网的经济效益。
3)以微电网为平台集成利用具有一定出行规律的EV与高渗透率光电,能够有效发挥EV集群的分布式储能作用、降低储能配置需求,提高微电网对光伏的消纳能力,是一种理想的协同增效利用模式。
附图说明
图1为年光照辐照度(即光照强度)曲线;
图2为微电网年常规负荷曲线;
图3为微电网日常规负荷和三种典型天气下的日光伏出力(即光伏发电功率)曲线;
图4为微电网总成本和光伏发电利用率φPPUR与储能配置量(即储能系统蓄电池的容量)的关系;
图5为不同电动汽车总数下微电网总成本与储能配置量的关系;
图6为本发明的不同电动汽车总数下φPPUR与储能配置量的关系。
具体实施方式
一种含高渗透率光电微电网中的电动汽车优化调度方法,所述的微电网由高渗透率光伏发电系统、储能系统、常规负荷和电动汽车集群构成,常规负荷和电动汽车集群的充电负构成荷微电网总负荷,电动汽车集群由若干电动汽构成,电动汽车既可作为充电负荷又可作为移动储能单元向微电网供电,微电网内的高渗透率的光伏发电系统、储能系统和作为移动储能单元时的电动汽车优先给微电网供电,富余的电量送到地区大电网,当微电网供电不足时,再从地区大电网购电,以满足微电网总负荷的需求;包括以下步骤:
S1:将全天24h均分为J个时段,对于任意第k时段,有k∈{1,2,...,J},且第k时段的时长为Δt,获取光伏发电的光照强度和环境温度信息、储能系统的蓄电池容量信息以及微电网的常规负荷信息;
S2:记录当前接入微电网的电动汽车的电池信息和客户充电需求信息;
S3:根据步骤S1、S2获取的信息基于虚拟费用构建电动汽车集群的NCD博弈模型;
S4:求解所构建的NCD博弈模型,以最大化利用光伏发电为目的,获取电动汽车的最优充放电计划,以实现微电网内电动汽车的优化调度。
进一步,本发明步骤S2中用一个七维行向量记录电动汽车的电池信息和客户充电需求信息,假设接入微电网的电动汽车的总数为N,对于任一电动汽车l∈{1,2,...,N},则有:
Xl=[Tin,l,Tout,l,S0,l,SE,l,Qs,l,Pc,l,Pd,l] (1)
式中:
Tin,l、Tout,l分别表示电动汽车l接入微电网的时间和预期离开微电网的时间;
S0,l、SE,l分别表示电动汽车l的电池的起始SOC(State of Charge,荷电状态,电池剩余能量与电池容量的比值)和离开微电网时客户期望电池达到的SOC,0≤S0,l≤1、0≤SE,l≤1;
Qs,l表示电动汽车l的电池容量;
Pc,l、Pd,l分别表示电动汽车l的额定充、放电功率。
进一步,本发明所述步骤S3具体包括以下步骤:
S31.以所有接入微电网的电动汽车和微电网内储能系统的蓄电池为博弈参与者,即决策主体U={1,2,...,N,N+1},构建任一博弈参与者i∈U的充放电计划模型;
对于博弈参与者之一的电动汽车l,其在全天J个时段内的充放电计划可表示如下:
Pl=(Pl(1),…,Pl(k),…,Pl(J)) (2)
其中,Pl为电动汽车l的充放电计划;Pl(k)表示第k时段,电动汽车l的充放电功率;
视电动汽车l的电池以恒功率充放电,忽略电池自身放电的影响,建立电动汽车l的电池模型和约束条件分别如下:
Sl(k)=Sl(k-1)+Pl(k)η(Pl(k))Δt/Qs,l (3)
SEV,min≤Sl(k)≤SEV,max (4)
-Pd,l≤Pl(k)≤Pc,l (5)
其中:
式(3)是电动汽车l的电池模型,Sl(k-1)、Sl(k)分别表示电动汽车l的电池在第k-1个和第k个时段的SOC,η(Pl(k))表示第k时段的电动汽车l的充放电效率;
SEV,max、SEV,min分别为电动汽车l的电池的SOC的最大值和最小值;
式(4)、(5)分别为电动汽车l的电池的SOC的约束和充放电功率Pl(k)的约束;
式(6)表示客户充电需求的约束,即电动汽车l如约离开时,电动汽车l的电池的SOC需满足客户期望;
微电网内储能系统的蓄电池是另一个博弈参与者,可视为一类无充电负荷需求的特殊“电动汽车”,其充放电计划和约束条件可类比于式(2)、式(4)~(5);此外,储能系统的蓄电池需满足的等式约束为:
即全天J个时段内储能系统的蓄电池的总充、放电电量需相等以满足充放电的循环,式中分别为储能系统的蓄电池的充、放电功率;ηBc、ηBd分别表示储能系统的蓄电池的充、放电效率;
S32.生成所有博弈参与者的策略空间,即生成所有博弈参与者的充放电计划集合[P1,P2,…,Pi,…,PN,PN+1],其中Pi=(Pi(1),Pi(2),...Pi(k),...,Pi(J))表示博弈参与者i在J个时段内充放电计划的集合,Pi(k)表示博弈参与者i在第k个时段的充放功率;在生成策略空间的过程中,除了考虑步骤S31中如式(4)~(7)所述的博弈参与者的运行约束外,还应考虑微电网与地区大电网之间的关系约束,主要包括:
Pgridout≤Pgridout,max (9)
式(8)表示微电网与地区大电网之间的功率传输平衡约束,其中,PPV(k)表示第k时段,高渗透率光伏发电系统的光伏发电输出功率,Pgrid(k)表示第k时段微电网与地区大电网的交互功率,Pgrid(k)>0表示微电网从地区大电网购电,Pgrid(k)<0表示微电网向地区大电网倒送输电,LB(k)表示第k时段微电网的常规负荷;
式(9)表示微电网向地区大电网倒送输电的倒送功率约束,其中,Pgridout表示微电网向地区大电网倒送输电的倒送功率,Pgridout,max为倒送功率允许的最大值;
综上,对于任一博弈参与者i的策略空间ρi可表示为:
ρi={Pi|constraint(2),(4),(5),(6),(7),(8),(9)} (10)
式(10)表示博弈参与者i的策略空间ρi包括博弈参与者i在J个时段内的充放电计划集合以及所要满足的约束(Constraint)条件:式(2)、式(4)~式(9);
S33.电动汽车和储能系统的蓄电池接入微电网后,计算博弈参与者i在J个时段内获得的虚拟利益;首先计算所有博弈参与者i的总虚拟利益:
式中:
γa表示总虚拟利益;γi表示博弈参与者i的虚拟利益;
γ(LB(k),k)=σ(LB(k)-PPV(k))2表示第k时段的基本虚拟费用;γ(La(k),k)=σ(La(k)-PPV(k))2表示第k时段的总虚拟费用;其中,表示第k时段微电网总负荷;
σ表示虚拟系数,且σ>0;
与实际电费不同,虚拟费用仅作为规划博弈参与者充放电计划的依据,用以制定各博弈参与者的最佳充放电计划;对于任一博弈参与者i的虚拟利益γi与博弈参与者i的充电需求以及负荷转移能力密切相关,假设博弈参与者i的负荷转移能力与其电池容量正相关,则博弈参与者i分享总虚拟利益的比例为:
式中,εi为博弈参与者i分享总虚拟利益的比例;Eneed,i表示博弈参与者i的充电需求量;κ表示正相关系数;特别地,由于接入微电网的电动汽车总数为N,因此当i∈{1,2,...,N}时,博弈参与者i指电动汽车,Eneed,i≠0;当i=N+1时博弈参与者i指储能系统的蓄电池,且有Eneed,i=0,Qs,i表示储能系统的蓄电池的配置容量;
根据式(11)和式(12)便可计算博弈参与者i的虚拟利益γi:
γi=εiγa (13)
S34.建立用于度量博弈参与者i参与博弈后的成本或效益的效用函数;定义博弈参与者i的虚拟利益γi为博弈参与者的效用函数,则博弈参与者i在充放电计划集合[P1,P2,…,Pi,…,PN,PN+1]下的效用函数Fi(Pi,Pi -)表示为:
式中,Pi -=[P1,P2,…,Pi-1,Pi+1,…,PN,PN+1],表示除博弈参与者i之外的其他所有博弈参与者的充放电计划集合;表示第k时段微电网内除博弈参与者i的负荷之外的其他负荷之和,其中τ表示除博弈参与者i之外的任一博弈参与者;
S35.构建决策主体的NCD博弈模型,博弈模型中,所有博弈参与者都是独立的决策者,各博弈参与者通过观察微电网的常规负荷、光伏发电输出功率以及其他博弈参与者的充放电计划来确定自身的最优充放电计划,以最大化虚拟利益;根据步骤S31~S34建立NCD博弈模型如下:
GNCD={U,{ρi}i∈U,{Fi}i∈U} (15)
基于上述NCD博弈模型,所有博弈参与者均会在假定其余博弈参与者都已选最优充放电计划的前提下选择能使其虚拟利益最大的最优充放电计划,即
其中,Pi *为博弈参与者i的最优充放电计划;
在NCD博弈模型中,如果对于博弈参与者充放电计划的其中一个集合满足如下条件:
Fi(Pi *,Pi *-)≥Fi(Pi,Pi *-),Pi∈ρi,i∈U (17)
则称P*为该NCD博弈模型的纳什均衡,即P*为各博弈参与者的最优充放电计划集合;其中Pi *-表示除博弈参与者i之外的其他所有博弈参与者的最优充放电计划集合。
进一步,本发明所述步骤S4中,基于NCD博弈模型的电动汽车的最优充放电计划求解过程具体如下:
(1)初始化电动汽车和储能系统的蓄电池的参数、第k时段配电网常规负荷LB(k)信息和光伏发电输出功率PPV(k);
(2)能量管理系统实时更新该时段的微电网总负荷La(k),并分享给所有博弈参与者;
(3)对于任一博弈参与者i∈U,其收到微电网总负荷La(k)信息后,计算出并以最大化式(14)为目标计算得到Pi *;
(4)如果计算得到的Pi *与当前最优值比较发生变化,更新Pi *,并将Pi *反馈给能量管理系统,转至(2);
(5)重复步骤(2)~(4),直至结果收敛(即Pi *不再发生变化),输出最优充放电计划Pi *;
为提高NCD博弈模型收敛至纳什均衡点的速度,步骤(4)中的各参与者可采用异步更新方式,即所有参与者不在同一时刻更新其充放电计划。
以某办公楼所在微电网为例验证本发明所述的电动汽车优化调度方法的有效性。该微电网的光伏发电系统的装机容量为1500kW,该地区年辐照度曲线如附图1所示,微电网年常规负荷曲线如附图2所示;微电网日常规负荷曲线和三种典型天气下的日光伏出力曲线如附图3所示。
相关参数设置如表1所示,由于本实施例针对办公楼区域微电网进行能量调度,车辆接入时段为高峰电价时段,车主并不希望在该种情况下为电动汽车补充能量,因此本实施例设SE,l=S0,l,微电网相关设备规格及价格参数如表2所示;办公楼区域峰谷电价如表3所示。
表1 参数设置
表2 设备参数
表3 办公楼区域峰谷电价
基于蒙特卡洛仿真方法模拟电动汽车的充电行为,根据所述的客户日充电需求数据分布情况,随机抽取电动汽车的电池的初始SOC、充电起止时间,根据本发明提出的优化调度方法,在Matlab中进行建模、求解,得到几种典型配置方案下的光伏电利用率φPPUR、微电网总等年值成本等微电网运行性能仿真结果,如表4所示。
表4 几种典型配置方案下的微电网系统运行性能对比
表4中,倒送量和购电量分别指微电网向地区大电网的倒送量和购电量;PV成本是光伏发电总等年值成本;储能成本指储能系统的总等年值成本;Grid成本指微电网向地区大电网年购电成本,EV补偿指为EV辅助服务年补偿费用;PV补偿指光伏年发电补贴;总成本指微电网总等年值成本。
微电网总成本和φPPUR与储能系统蓄电池的配置量的关系如附图4所示;从图中可以看出,微电网总成本曲线呈现“U”型,即随着储能系统储能容量增加,微电网总成本先减少后增加,在储能容量达到2200kW·h时,微电网经济性最优,总成本最低为186.19万元;另一方面,随着储能容量增加,微电网内的光电利用率持续增加,呈现单调向好趋势。
对比表4中情况1和情况2可以看出,尽管增加储能成本32.87万元,情况2仍具有较好的经济性,总成本相比于情况1减少了11.28万元;主要原因在于添加储能系统组成光储微电网后,弃光量减少,φPPUR从74.67%增加到88.06%,从而使得微电网净出力减少,导致向地区大电网购电成本减少25.95万元以及光伏发电补贴增加16.82万元。
可见单一的分布式光伏发电系统增加储能系统组成光储微电网后可以有效提高φPPUR和系统经济性,弃光问题得到有效的改善。
对比表4中情况2和情况3可以看出,在储能系统蓄电池配置量相同情况下,接入有序控制下的EV集群,可以使微电网净出力和光伏弃电量有所减少,使得向地区大电网购电成本减少11.11万元及光伏发电补贴增加11.22万元;φPPUR由88.06%提高到97.00%;进一步分析,可见EV集群改善微电网运行效果明显,在本发明设定的负荷转移补偿价格下EV集群的补偿成本14.25万元,但微电网总成本减少达10.72万元,可见EV集群的接入,能够有效提高φPPUR、降低微电网总成本。
对比表4中情况2和情况4可以看出,在接入总数为20辆的EV集群后,最优储能容量配置需求从2200kW·h降低至1200kW·h,使储能成本减少14.86万元,情况4微电网总成本比情况2总成本减少12.83万元;可见从经济性角度出发,V2G模式下EV的优化调度可以有效减少微电网储能配置容量,进而获得更好的经济效益。
附图5为不同EV接入规模下系统总成本与储能配置容量的关系;从图中可以看出,EV接入规模增加,最优储能配置需求减少,系统总成本会进一步降低。附图6为不同EV接入规模下φPPUR与储能配置容量的关系,同一电动汽车总数下,φPPUR随储能配置量的增加而持续增大,直至光电被完全消纳时稳定;在光电饱和利用之前,相同储能配置下,接入车辆总数越大,φPPUR越大。本说明书实施例所述的内容仅仅是对发明构思的实现形式的列举,本发明的保护范围不应当被视为仅限于实施例所陈述的具体形式,本发明的保护范围也包括本领域技术人员根据本发明构思所能够想到的等同技术手段。
Claims (3)
1.一种含高渗透率光电微电网中的电动汽车优化调度方法,所述的微电网由高渗透率光伏发电系统、储能系统、常规负荷和电动汽车集群构成,常规负荷和电动汽车集群的充电负荷构成微电网总负荷,电动汽车集群由若干电动汽车构成,电动汽车既可作为充电负荷又可作为移动储能单元向微电网供电,微电网内的高渗透率的光伏发电系统、储能系统和作为移动储能单元时的电动汽车优先给微电网供电,富余的电量送到地区大电网,当微电网供电不足时,再从地区大电网购电,以满足微电网总负荷的需求;其特征在于,包括以下步骤:
S1:将全天24h均分为J个时段,对于任意第k时段,有k∈{1,2,...,J},且第k时段的时长为Δt,获取光伏发电的光照强度和环境温度信息、储能系统的蓄电池容量信息以及微电网的常规负荷信息;
S2:记录当前接入微电网的电动汽车的电池信息和客户充电需求信息;
S3:根据步骤S1、S2获取的信息基于虚拟费用构建电动汽车集群的NCD博弈模型,NCD表示非合作充放电;
S4:求解所构建的NCD博弈模型,以最大化利用光伏发电为目的,获取电动汽车的最优充放电计划,以实现微电网内电动汽车的优化调度;
所述步骤S3具体包括以下步骤:
S31.以所有接入微电网的电动汽车和微电网内储能系统的蓄电池为博弈参与者,即决策主体U={1,2,...,N,N+1},
构建任一博弈参与者i∈U的充放电计划模型;
对于博弈参与者之一的电动汽车l,其在全天J个时段内的充放电计划可表示如下:
Pl=(Pl(1),…,Pl(k),…,Pl(J)) (2)
其中,Pl为电动汽车l的充放电计划;Pl(k)表示第k时段,电动汽车l的充放电功率;
视电动汽车l的电池以恒功率充放电,忽略电池自身放电的影响,建立电动汽车l的电池模型和约束条件分别如下:
Sl(k)=Sl(k-1)+Pl(k)η(Pl(k))Δt/Qs,l (3)
SEV,min≤Sl(k)≤SEV,max (4)
-Pd,l≤Pl(k)≤Pc,l (5)
<mrow>
<msub>
<mi>S</mi>
<mrow>
<mn>0</mn>
<mo>,</mo>
<mi>l</mi>
</mrow>
</msub>
<mo>+</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>J</mi>
</munderover>
<msub>
<mi>P</mi>
<mi>l</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mi>&eta;</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>P</mi>
<mi>l</mi>
</msub>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
<msub>
<mi>Q</mi>
<mrow>
<mi>s</mi>
<mo>,</mo>
<mi>l</mi>
</mrow>
</msub>
</mfrac>
<mo>&GreaterEqual;</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>E</mi>
<mo>,</mo>
<mi>l</mi>
</mrow>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
其中:
式(3)是电动汽车l的电池模型,Sl(k-1)、Sl(k)分别表示电动汽车l的电池在第k-1个和第k个时段的SOC,η(Pl(k))表示第k时段的电动汽车l的充放电效率;
SEV,max、SEV,min分别为电动汽车l的电池的SOC的最大值和最小值;
式(4)、(5)分别为电动汽车l的电池的SOC的约束和充放电功率Pl(k)的约束;
式(6)表示客户充电需求的约束,即电动汽车l如约离开时,电动汽车l的电池的SOC需满足客户期望;
微电网内储能系统的蓄电池是另一个博弈参与者,可视为一类无充电负荷需求的特殊“电动汽车”,其充放电计划和约束条件可类比于式(2)、式(4)~(5);此外,储能系统的蓄电池需满足的等式约束为:
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>J</mi>
</munderover>
<mrow>
<mo>(</mo>
<msubsup>
<mi>P</mi>
<mi>B</mi>
<mo>+</mo>
</msubsup>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
<msub>
<mi>&eta;</mi>
<mrow>
<mi>B</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>+</mo>
<mfrac>
<mrow>
<msubsup>
<mi>P</mi>
<mi>B</mi>
<mo>-</mo>
</msubsup>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
<msub>
<mi>&eta;</mi>
<mrow>
<mi>B</mi>
<mi>d</mi>
</mrow>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mi>&Delta;</mi>
<mi>t</mi>
<mo>=</mo>
<mn>0</mn>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
即全天J个时段内储能系统的蓄电池的总充、放电电量需相等以满足充放电的循环,式中分别为储能系统的蓄电池的充、放电功率;ηBc、ηBd分别表示储能系统的蓄电池的充、放电效率;
S32.生成所有博弈参与者的策略空间,即生成所有博弈参与者的充放电计划集合[P1,P2,…,Pi,…,PN,PN+1],其中Pi=(Pi(1),Pi(2),...Pi(k),...,Pi(J))表示博弈参与者i在J个时段内充放电计划的集合,Pi(k)表示博弈参与者i在第k个时段的充放功率;在生成策略空间的过程中,除了考虑步骤S31中如式(4)~(7)所述的博弈参与者的运行约束外,还应考虑微电网与地区大电网之间的关系约束,主要包括:
<mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>P</mi>
<mi>V</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>g</mi>
<mi>r</mi>
<mi>i</mi>
<mi>d</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>L</mi>
<mi>B</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>&Element;</mo>
<mi>U</mi>
</mrow>
</munder>
<msub>
<mi>P</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
Pgridout≤Pgridout,max (9)
式(8)表示微电网与地区大电网之间的功率传输平衡约束,其中,PPV(k)表示第k时段,高渗透率光伏发电系统的光伏发电输出功率,Pgrid(k)表示第k时段微电网与地区大电网的交互功率,Pgrid(k)>0表示微电网从地区大电网购电,Pgrid(k)<0表示微电网向地区大电网倒送输电,LB(k)表示第k时段微电网的常规负荷;
式(9)表示微电网向地区大电网倒送输电的倒送功率约束,其中,Pgridout表示微电网向地区大电网倒送输电的倒送功率,Pgridout,max为倒送功率允许的最大值;
综上,对于任一博弈参与者i的策略空间ρi可表示为:
ρi={Pi|constraint(2),(4),(5),(6),(7),(8),(9)} (10)
式(10)表示博弈参与者i的策略空间ρi包括博弈参与者i在J个时段内的充放电计划集合以及所要满足的约束(Constraint)条件:式(2)、式(4)~式(9);
S33.电动汽车和储能系统的蓄电池接入微电网后,计算博弈参与者i在J个时段内获得的虚拟利益;首先计算所有博弈参与者i的总虚拟利益:
<mrow>
<msub>
<mi>&gamma;</mi>
<mi>a</mi>
</msub>
<mo>=</mo>
<munder>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>&Element;</mo>
<mi>U</mi>
</mrow>
</munder>
<msub>
<mi>&gamma;</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>J</mi>
</munderover>
<mo>&lsqb;</mo>
<mi>&gamma;</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>L</mi>
<mi>B</mi>
</msub>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
<mo>,</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>&gamma;</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>L</mi>
<mi>a</mi>
</msub>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
<mo>,</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>11</mn>
<mo>)</mo>
</mrow>
</mrow>
式中:
γa表示总虚拟利益;γi表示博弈参与者i的虚拟利益;
γ(LB(k),k)=σ(LB(k)-PPV(k))2表示第k时段的基本虚拟费用;γ(La(k),k)=σ(La(k)-PPV(k))2表示第k时段的总虚拟费用;其中,表示第k时段微电网总负荷;
σ表示虚拟系数,且σ>0;
与实际电费不同,虚拟费用仅作为规划博弈参与者充放电计划的依据,用以制定各博弈参与者的最佳充放电计划;对于任一博弈参与者i的虚拟利益γi与博弈参与者i的充电需求以及负荷转移能力密切相关,假设博弈参与者i的负荷转移能力与其电池容量正相关,则博弈参与者i分享总虚拟利益的比例为:
<mrow>
<msub>
<mi>&epsiv;</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>&kappa;Q</mi>
<mrow>
<mi>s</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>E</mi>
<mrow>
<mi>n</mi>
<mi>e</mi>
<mi>e</mi>
<mi>d</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
</mrow>
<mrow>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>&Element;</mo>
<mi>U</mi>
</mrow>
</munder>
<mrow>
<mo>(</mo>
<msub>
<mi>&kappa;Q</mi>
<mrow>
<mi>s</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>E</mi>
<mrow>
<mi>n</mi>
<mi>e</mi>
<mi>e</mi>
<mi>d</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>12</mn>
<mo>)</mo>
</mrow>
</mrow>
式中,εi为博弈参与者i分享总虚拟利益的比例;Eneed,i表示博弈参与者i的充电需求量;κ表示正相关系数;特别地,由于接入微电网的电动汽车总数为N,因此当i∈{1,2,...,N}时,博弈参与者i指电动汽车,Eneed,i≠0;当i=N+1时博弈参与者i指储能系统的蓄电池,且有Eneed,i=0,Qs,i表示储能系统的蓄电池的配置容量;
根据式(11)和式(12)便可计算博弈参与者i的虚拟利益γi:
γi=εiγa (13)
S34.建立用于度量博弈参与者i参与博弈后的成本或效益的效用函数;定义博弈参与者i的虚拟利益γi为博弈参与者的效用函数,则博弈参与者i在充放电计划集合[P1,P2,…,Pi,…,PN,PN+1]下的效用函数表示为:
<mrow>
<msub>
<mi>F</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>P</mi>
<mi>i</mi>
</msub>
<mo>,</mo>
<msubsup>
<mi>P</mi>
<mi>i</mi>
<mo>-</mo>
</msubsup>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>&gamma;</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<msub>
<mi>&epsiv;</mi>
<mi>i</mi>
</msub>
<mi>&sigma;</mi>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>J</mi>
</munderover>
<mo>{</mo>
<msup>
<mrow>
<mo>&lsqb;</mo>
<msub>
<mi>L</mi>
<mi>B</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>P</mi>
<mi>V</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>-</mo>
<msup>
<mrow>
<mo>&lsqb;</mo>
<msubsup>
<mi>L</mi>
<mi>a</mi>
<mrow>
<mo>-</mo>
<mi>i</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>P</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>P</mi>
<mi>V</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>}</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>14</mn>
<mo>)</mo>
</mrow>
</mrow>
式中,表示除博弈参与者i之外的其他所有博弈参与者的充放电计划集合;表示第k时段微电网内除博弈参与者i的负荷之外的其他负荷之和,其中τ表示除博弈参与者i之外的任一博弈参与者;
S35.构建决策主体的NCD博弈模型,博弈模型中,所有博弈参与者都是独立的决策者,各博弈参与者通过观察微电网的常规负荷、光伏发电输出功率以及其他博弈参与者的充放电计划来确定自身的最优充放电计划,以最大化虚拟利益;根据步骤S31~S34建立NCD博弈模型如下:
GNCD={U,{ρi}i∈U,{Fi}i∈U} (15)
基于上述NCD博弈模型,所有博弈参与者均会在假定其余博弈参与者都已选最优充放电计划的前提下选择能使其虚拟利益最大的最优充放电计划,即
<mrow>
<msubsup>
<mi>P</mi>
<mi>i</mi>
<mo>*</mo>
</msubsup>
<mo>&Element;</mo>
<munder>
<mi>argmax</mi>
<mrow>
<msub>
<mi>P</mi>
<mi>i</mi>
</msub>
<mo>&Element;</mo>
<msub>
<mi>&rho;</mi>
<mi>i</mi>
</msub>
</mrow>
</munder>
<msub>
<mi>F</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>16</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,为博弈参与者i的最优充放电计划;
在NCD博弈模型中,如果对于博弈参与者充放电计划的其中一个集合满足如下条件:
<mrow>
<msub>
<mi>F</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<msubsup>
<mi>P</mi>
<mi>i</mi>
<mo>*</mo>
</msubsup>
<mo>,</mo>
<msubsup>
<mi>P</mi>
<mi>i</mi>
<mrow>
<mo>*</mo>
<mo>-</mo>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mo>&GreaterEqual;</mo>
<msub>
<mi>F</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>P</mi>
<mi>i</mi>
</msub>
<mo>,</mo>
<msubsup>
<mi>P</mi>
<mi>i</mi>
<mrow>
<mo>*</mo>
<mo>-</mo>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>P</mi>
<mi>i</mi>
</msub>
<mo>&Element;</mo>
<msub>
<mi>&rho;</mi>
<mi>i</mi>
</msub>
<mo>,</mo>
<mi>i</mi>
<mo>&Element;</mo>
<mi>U</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>17</mn>
<mo>)</mo>
</mrow>
</mrow>
则称P*为该NCD博弈模型的纳什均衡,即P*为各博弈参与者的最优充放电计划集合;其中表示除博弈参与者i之外的其他所有博弈参与者的最优充放电计划集合。
2.根据权利要求1所述的一种含高渗透率光电微电网中的电动汽车优化调度方法,其特征在于,步骤S2中用一个七维行向量记录电动汽车的电池信息和客户充电需求信息,假设接入微电网的电动汽车的总数为N,对于任一电动汽车l∈{1,2,...,N},则有:
Xl=[Tin,l,Tout,l,S0,l,SE,l,Qs,l,Pc,l,Pd,l] (1)
式中:
Tin,l、Tout,l分别表示电动汽车l接入微电网的时间和预期离开微电网的时间;
S0,l、SE,l分别表示电动汽车l的电池的起始SOC和离开微电网时客户期望电池达到的SOC,SOC是荷电状态,表示电池剩余能量与电池容量的比值,且0≤S0,l≤1、0≤SE,l≤1;
Qs,l表示电动汽车l的电池容量;
Pc,l、Pd,l分别表示电动汽车l的额定充、放电功率。
3.根据权利要求2所述的一种含高渗透率光电微电网中的电动汽车优化调度方法,其特征在于,所述步骤S4中,基于NCD博弈模型的电动汽车的最优充放电计划求解过程具体如下:
(1)初始化电动汽车和储能系统的蓄电池的参数、第k时段配电网常规负荷LB(k)信息和光伏发电输出功率PPV(k);
(2)能量管理系统实时更新该时段的微电网总负荷La(k),并分享给所有博弈参与者;
(3)对于任一博弈参与者i∈U,其收到微电网总负荷La(k)信息后,计算出并以最大化式(14)为目标计算得到
(4)如果计算得到的与当前最优值比较发生变化,更新并将反馈给能量管理系统,转至(2);
(5)重复(2)~(4),直至结果收敛(即不再发生变化),输出最优充放电计划
为提高NCD博弈模型收敛至纳什均衡点的速度,步骤(4)中的各参与者可采用异步更新方式,即所有参与者不在同一时刻更新其充放电计划。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610071449.6A CN105576684B (zh) | 2016-02-01 | 2016-02-01 | 一种含高渗透率光电微电网中的电动汽车优化调度方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610071449.6A CN105576684B (zh) | 2016-02-01 | 2016-02-01 | 一种含高渗透率光电微电网中的电动汽车优化调度方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105576684A CN105576684A (zh) | 2016-05-11 |
CN105576684B true CN105576684B (zh) | 2018-06-01 |
Family
ID=55886518
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610071449.6A Active CN105576684B (zh) | 2016-02-01 | 2016-02-01 | 一种含高渗透率光电微电网中的电动汽车优化调度方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105576684B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107017619B (zh) * | 2017-03-29 | 2019-09-27 | 华北电力大学 | 非合作博弈视角的光伏充电站网络分布式能量管理方法 |
CN108805325B (zh) * | 2018-04-11 | 2022-03-01 | 杭州电子科技大学 | 一种生产计划与调度集成优化方法 |
CN108923536B (zh) * | 2018-07-12 | 2020-09-08 | 中国南方电网有限责任公司 | 可调度潜力分析方法、系统、计算机设备及存储介质 |
CN110289611A (zh) * | 2019-03-29 | 2019-09-27 | 国网浙江海宁市供电有限公司 | 一种提升分布式电源消纳能力的主动配电网规划方法 |
CN111799786B (zh) * | 2020-07-07 | 2021-03-16 | 中原工学院 | 一种基于博弈论的新能源混合电力系统容量配置方法 |
CN117485147B (zh) * | 2024-01-02 | 2024-03-26 | 国网湖北省电力有限公司营销服务中心(计量中心) | 一种电动汽车虚拟储能优化方法和系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104218597A (zh) * | 2014-08-25 | 2014-12-17 | 浙江工业大学 | 微网下v2g的智能充放电管理方法 |
JP5680222B2 (ja) * | 2011-12-27 | 2015-03-04 | 三菱電機株式会社 | エネルギーマネジメントシステム |
-
2016
- 2016-02-01 CN CN201610071449.6A patent/CN105576684B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5680222B2 (ja) * | 2011-12-27 | 2015-03-04 | 三菱電機株式会社 | エネルギーマネジメントシステム |
CN104218597A (zh) * | 2014-08-25 | 2014-12-17 | 浙江工业大学 | 微网下v2g的智能充放电管理方法 |
Non-Patent Citations (1)
Title |
---|
含光伏发电系统的电动汽车充电站多目标容量优化配置方法;陈征等;《电工技术学报》;20130630;第28卷(第6期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN105576684A (zh) | 2016-05-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105576684B (zh) | 一种含高渗透率光电微电网中的电动汽车优化调度方法 | |
Hou et al. | Multi-objective economic dispatch of a microgrid considering electric vehicle and transferable load | |
CN106169108B (zh) | 含电池储能系统的主动配电网短期有功优化方法 | |
CN107017619B (zh) | 非合作博弈视角的光伏充电站网络分布式能量管理方法 | |
CN110350523B (zh) | 基于需求响应的多能源互补优化调度方法 | |
Liu et al. | Online energy sharing for nanogrid clusters: A Lyapunov optimization approach | |
CN109599856B (zh) | 一种微网多楼宇中电动汽车充放电管理优化方法及装置 | |
Bayod-Rujula et al. | Sizing criteria of hybrid photovoltaic–wind systems with battery storage and self-consumption considering interaction with the grid | |
CN109657993B (zh) | 一种基于非合作博弈的能源局域网储能系统自动需求响应方法 | |
Tao et al. | Pricing strategy and charging management for PV-assisted electric vehicle charging station | |
Kheradmand-Khanekehdani et al. | Well-being analysis of distribution network in the presence of electric vehicles | |
CN113240210B (zh) | 一种基于多主体主从博弈的负荷聚合商经济优化方法 | |
Arab et al. | Suitable various-goal energy management system for smart home based on photovoltaic generator and electric vehicles | |
Ren et al. | Multitime scale coordinated scheduling for electric vehicles considering photovoltaic/wind/battery generation in microgrid | |
CN116402307A (zh) | 考虑可调度柔性资源运行特性的电网规划容量分析方法 | |
Yuanyuan et al. | Optimization scheduling method of power grid energy-saving based on fuzzy decision | |
Fu et al. | Research on the stackelberg game method of building micro-grid with electric vehicles | |
Arab et al. | A priority-based seven-layer strategy for energy management cooperation in a smart city integrated green technology | |
CN110991928B (zh) | 一种多微能源网的综合能源系统能量管理方法及系统 | |
Phonphan et al. | Home Energy Management System Based on The Photovoltaic–Battery Hybrid Power System | |
CN111987717A (zh) | 车联网充放电管理方法及系统、电力系统及存储介质 | |
CN108736496B (zh) | 一种分布式储能系统的补充性规划方法及系统 | |
CN106253356B (zh) | 基于储能电量价值评估的交直流混合微网实时经济调度法 | |
CN111952996B (zh) | 一种基于经济效益评估的含储能的分布式光伏调峰控制方法 | |
Wu et al. | Optimal operation of residential microgrids in the Harbin area |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |