CN104734171B - 一种用于配电网可靠性评估的电动汽车充电站建模方法及其应用 - Google Patents
一种用于配电网可靠性评估的电动汽车充电站建模方法及其应用 Download PDFInfo
- Publication number
- CN104734171B CN104734171B CN201510180855.1A CN201510180855A CN104734171B CN 104734171 B CN104734171 B CN 104734171B CN 201510180855 A CN201510180855 A CN 201510180855A CN 104734171 B CN104734171 B CN 104734171B
- Authority
- CN
- China
- Prior art keywords
- electric automobile
- sigma
- charging station
- formula
- distribution network
- 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
- 238000009826 distribution Methods 0.000 title claims abstract description 97
- 238000000034 method Methods 0.000 title claims abstract description 34
- 230000002123 temporal effect Effects 0.000 claims abstract description 28
- 238000000342 Monte Carlo simulation Methods 0.000 claims abstract description 12
- 238000005070 sampling Methods 0.000 claims description 39
- 230000007257 malfunction Effects 0.000 claims description 19
- 238000004088 simulation Methods 0.000 claims description 18
- 238000012546 transfer Methods 0.000 claims description 15
- 239000012141 concentrate Substances 0.000 claims description 10
- 230000005611 electricity Effects 0.000 claims description 10
- 238000010586 diagram Methods 0.000 claims description 7
- 239000011159 matrix material Substances 0.000 claims description 4
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 3
- 230000002085 persistent effect Effects 0.000 claims description 3
- 235000006508 Nelumbo nucifera Nutrition 0.000 claims 1
- 240000002853 Nelumbo nucifera Species 0.000 claims 1
- 235000006510 Nelumbo pentapetala Nutrition 0.000 claims 1
- 238000007599 discharging Methods 0.000 claims 1
- 238000004146 energy storage Methods 0.000 abstract description 9
- 238000011160 research Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000009977 dual effect Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000033001 locomotion Effects 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 101100426970 Caenorhabditis elegans ttr-1 gene Proteins 0.000 description 1
- 238000007476 Maximum Likelihood Methods 0.000 description 1
- 101100460495 Rattus norvegicus Nkx2-1 gene Proteins 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004140 cleaning Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005315 distribution function Methods 0.000 description 1
- 238000003912 environmental pollution Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000013277 forecasting method Methods 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 238000011144 upstream manufacturing Methods 0.000 description 1
Classifications
-
- H02J7/0027—
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Charge And Discharge Circuits For Batteries Or The Like (AREA)
- Secondary Cells (AREA)
- Electric Propulsion And Braking For Vehicles (AREA)
Abstract
本发明公开了一种用于配电网可靠性评估的电动汽车充电站建模方法及其应用,其特征是按如下步骤进行:1建立电动汽车运行状态的时序模型;2建立电动汽车时空分布模型;3基于dijkstra算法计算各充电站可集中充电的动力电池的总数;4建立基于可靠性评估的充电站出力时序模型;5采用蒙特卡洛模拟的方法对含电动汽车的配电网进行可靠性评估,统计可靠性指标。本发明能建立更加符合实际的电动汽车出力模型,从而提高电动汽车集中优化控制及其作为储能装置的利用率,有助于缓解配电网负荷高峰时缺电压力,进一步提高配电网的可靠性。
Description
技术领域
本发明涉及配电网可靠性评估领域,具体涉及一种用于配电网可靠性评估的电动汽车充电站建模方法及其应用。
背景技术
在环境污染和能源短缺的双重压力下,具有清洁、无污染和节能优势的电动汽车已经成为新能源领域的重要研究内容。从电力系统角度对电动汽车的研究主要集中在其动力电池属性上,大规模电动汽车接入后将使配电网的电能峰值需求增加,无序模式下的充电负荷势必会对配电系统的可靠性造成一定的影响。
随着电动汽车接入电网(Vehicle to Grid,V2G)技术的发展,电动汽车接入作为充电负荷的同时也可以作为移动的储能装置向电网放电,从而可以提高配电网的可靠性;此外可入网电动汽车的接入使得配电网在结构和运行方式上与传统配电网有较大差异,传统的可靠性评估方法存在一定的缺陷,因此研究含电动汽车的配电网可靠性评估对充分发挥电动汽车的双重角色具有重要的现实意义。
目前将电动汽车看作储能装置接入对配电网可靠性的影响是一个比较新颖的问题。当前作为储能装置的电动汽车研究主要集中在利用其调节电网负荷峰谷差,消纳新能源及在分布式电源和微电网建设中的应用,而在可靠性评估层面,现有的研究主要探讨了电动汽车接入的情况下发输电系统的可靠性评估,并没涉及配电网可靠性评估且电动汽车仅作为充电负荷,增加了配电网的峰值负荷压力且无法发挥电动汽车作为储能的作用。有少数虽然研究了电动汽车负荷对配电网可靠性影响的量化分析,但没有涉及到电动汽车出力建模,此外对于随机分散在配电网中的电动汽车简单的采用V2G技术恢复供电过于理想、操作难度较大且很难治理谐波污染。
发明内容
本发明是为避免上述现有技术所存在的不足之处,提供一种用于配电网可靠性评估的电动汽车建模方法及其应用,以期能建立更加符合实际的电动汽车出力模型,从而提高电动汽车集中优化控制及其作为储能装置的利用率,有助于缓解配电网负荷高峰时缺电压力,进一步提高配电网的可靠性。
为了达到上述目的,本发明所采用的技术方案为:
本发明一种用于配电网可靠性评估的电动汽车充电站建模方法及其应用,所述配电网含有电动汽车;所述电动汽车包括动力电池和电机并在所述配电网中通过充电站和集中充电模式进行充电;假设所述配电网中含有H辆电动汽车,记为V={V1,V2,…,Vm,…,VH};Vm表示第m辆电动汽车;1≤m≤H;其特点是,所述方法是按如下步骤进行:
步骤一、建立所述H辆电动汽车运行状态的时序模型:
步骤1.1、根据所述H辆电动汽车的动力电池和电机的运行状态,将所述第m辆电动汽车的运行状态离散为四个状态,分别为和表示所述第m辆电动汽车动力电池和电机的正常工作状态,表示所述第m辆电动汽车动力电池的故障状态和电机的正常工作状态、表示所述第m辆电动汽车动力电池的正常工作状态和电机的故障状态、表示所述第m辆电动汽车动力电池和电机的故障状态;
步骤1.2、建立第m辆电动汽车的运行状态空间图,并获得所述运行状态空间图的第m个状态转移矩阵Mm;
步骤1.3、根据马尔科夫过程逼近原理,获得所述第m辆电动汽车四个状态的概率,分别为和表示所述第m辆电动汽车动力电池和电机的正常工作状态的概率;表示所述第m辆电动汽车动力电池的故障状态和电机的正常工作状态的概率、表示所述第m辆电动汽车动力电池的正常工作状态和电机的故障状态的概率、表示所述第m辆电动汽车动力电池和电机的故障状态的概率;并有
步骤1.4、采用蒙特卡洛方法抽取服从均匀分布的随机数和并利用式(2)确定所述第m辆电动汽车的运行状态Sm:
利用式(3)抽取所述第m辆电动汽车运行状态Sm的实际持续时间
式(3)中:Lm表示所述第m辆电动汽车从运行状态Sm转移到其他三种状态的转移次数,是表示所述电动汽车从运行状态Sm转移到其他三种状态过程中第h个转移次数时的转移率;
步骤1.5、设所述蒙特卡洛方法的仿真时间为Y年,在所述仿真时间Y年以ΔT为时间步长获得采样时刻ti,并有N表示一年的总小时数;
所述第m辆电动汽车的运行状态在所述四个状态中互相转移,从而在仿真时间Y年内形成第m辆电动汽车运行状态的时序模型ψ(m)[t],若采样时刻ti下第m辆电动汽车的运行状态为所述车动力电池和电机的正常工作状态则ψ(m)[ti]=1,否则ψ(m)[ti]=0;从而获得所述H辆电动汽车运行状态在仿真时间Y年内的时序模型为:
{ψ(1)[t],ψ(2)[t],…,ψ(m)[t],…,ψ(H)[t]};
步骤二、建立所述H辆电动汽车的时空分布模型:
步骤2.1、令Am事件表示所述第m辆电动汽车第一次出行开始时刻处在集中充电时段te~tf之间,Bm事件表示所述第m辆电动汽车最后一次出行结束时刻落在采样时刻ti的一个邻域内,记作则利用式(4)和式(5)分别获得所述Am事件的概率P(Am)和所述Bm事件的概率P(Am):
式(4)和式(5)中,fstart(x)和fend(x)分别表示所述第m辆电动汽车在24小时中第一次出行开始时刻和最后一次出行结束时刻的概率密度函数;
步骤2.2、在所述集中充电模式下,在配电网中利用式(6)获得采样时刻ti下所述第m辆电动汽车的可调度概率
步骤2.3、假设配电网中有W负荷节点,记为集合G={G1,G2,…,Gj,…,GW},Gj表示第j个负荷节点,1≤j≤W;利用式(7)建立采样时刻ti下配电网中第j个负荷节点Gj对所述第m辆动汽车吸引度模型
式(7)中,表示第j个负荷节点Gj对所述第m辆电动汽车吸引度的权重系数;表示第j个负荷节点Gj在采样时刻ti下的有功负荷;
步骤2.3、利用式(8)计算采样时刻ti下配网中第j个负荷节点Gj所接入的电动汽车分布数量nj(τ):
步骤三、基于dijkstra算法计算各充电站可集中充电的动力电池的总数:
步骤3.1、设配电网中有N座充电站,将所有带有充电站的负荷节点记为集合E={E1,E2,…,Eσ,…EN};Eσ表示第σ个带有充电站的负荷节点;1≤σ≤N<W;
步骤3.2、计算各充电站的负荷节点到配电网中其它W-1个负荷节点的最短路径:
步骤3.2.1、令所述配电网的W个负荷节点中任意两个负荷节点Ga和Gb之间的距离为l(a,b),1≤a,b≤W;其中,将所述第σ个带有充电站的负荷节点Eσ到第j个负荷节点Gj的距离记为定义循环次数z;1≤z≤W-1;初始化σ=1;z=1;
步骤3.2.2、令集合则集合在集合G中的补集记作利用式(9)获得所述集合中元素Eσ到补集中任意元素的距离
步骤3.2.3、根据 获得负荷节点并利用 和获得更新集合和更新集合
步骤3.2.4、以负荷节点为中间节点,对任意利用式(10)获得更新的距离
步骤3.2.5、z+1赋值给z,判断z=W是否成立,若成立,则执行步骤3.2.6;否则,返回步骤3.2.3;
步骤3.2.6、σ+1赋值给σ,判断σ=N+1,是否成立,若成立,则执行步骤3.3;否则,返回步骤3.2.2;
步骤3.3、利用式(11)获得第j个负荷节点Gj内的所有电动汽车对第σ个带有充电站的负荷节点Eσ的选择系数为
式(11)中,
步骤3.4、利用式(12)计算采样时刻ti内的第σ个带有充电站的负荷节点Eσ可集中充电的动力电池数量
步骤3.5、利用式(13)获得第σ个带有充电站的负荷节点Eσ采样时刻ti内可集中充电的动力电池的总数
步骤四、建立基于可靠性评估的充电站出力时序模型:
步骤4.1、假设第m辆电动汽车的动力电池在采样时刻ti进入充电站,则利用式(14)和式(15)分别获得第m辆电动汽车动力电池在配电网无故障状态的荷电状态变化特征SOC(ti+ΔT)m和在配电网故障状态下的荷电状态变化特征SOC(ti+ΔT)m:
式(14)和式(15)中:和分别为第m辆电动汽车动力电池的充电效率和放电效率;和分别表示第m辆电动汽车动力电池在采样时刻ti充电功率和放电功率,表示第m辆电动汽车电池的额定容量,SOC(ti)m表示第m辆电动汽车动力电池在采样时刻ti下的荷电状态值;
步骤4.2、假设在任意时刻tc下配电网发生故障,故障时间为ttr;若存在第σ个带有充电站的负荷节点Eσ位于故障区域内,则利用式(16)获得第σ个带有充电站的负荷节点Eσ在采样时刻ti下的时序输出功率从而获得所述N座充电站在仿真时间Y年内的出力时序模型为:
式(16)中:当SOC(ti)ω>0.2时,ζω(ti)=1,否则ζω(ti)=0。
本发明所述的用于配电网可靠性评估的电动汽车充电站建模方法的应用的特点是:
在所述仿真时间Y年内,抽取配电网中电气元件的状态,所述电气元件的状态分为故障状态和正常状态;以所述出力时序模型作为孤岛恢复供电电源,采用所述蒙特卡洛方法对所述配电网进行可靠性评估,并统计可靠性指标。
与已有技术相比,本发明的有益效果体现在:
1、本发明基于马尔科夫链解析方法建立了电动汽车运行状态时序模型,使充电站出力模型具有更高的准确性、更加全面的反映实际情况。
2、本发明将电动汽车时空分布模型和dijkstra算法应用于建立集中型充电站出力时序模型,不仅能有效处理电动汽车在配电网中分布的随机性,还解决了配电网中各充电站可集中充电的电动汽车数量问题。
3、本发明建立的集中型充电站出力时序模型有效避免了电动汽车充放电在时间和空间上的随机性,有利于充放电的调度控制,也为谐波的集中治理创造了条件。
4、本发明建立的集中型充电站出力时序模型能用于配电网的可靠性评估,可以有效的改善配电网的可靠性指标。
5、本发明将充电站作为储能装置可以在一定程度上提高配电网的可靠性,为电动汽车接入配电网后的集中优化控制提供参考,有效避免了电动汽车充电的随机性对电力系统所造成的冲击。
附图说明
图1为本发明电动汽车运行状态空间图;
图2为本发明配电网中电动汽车的吸引度建模流程。
具体实施方式
本实施例中,配电网含有电动汽车;电动汽车的动力电池采取“集中充电、统一配送”的原则,电动汽车包括动力电池和电机并在配电网中通过充电站和集中充电模式进行充电;假设配电网中含有H辆电动汽车,记为V={V1,V2,…,Vm,…,VH};Vm表示第m辆电动汽车;1≤m≤H;一种用于配电网可靠性评估的电动汽车充电站建模方法及其应用是采用马尔科夫链的解析方法,建立电动汽车运行状态的时序模型;根据电动汽车出行时刻概率密度函数和居民负荷与电动汽车空间分布的关联性,建立电动汽车时空分布模型;采用dijkstra算法求解出充电站到配电网节点的最短路径进而确定充电站的服务范围;结合电动汽车的出力特点和运行状态序列,得到用于可靠性评估的集中型充电站出力时序模型;采用蒙特卡洛模拟的方法对含电动汽车的配电网进行可靠性评估,统计可靠性指标。具体地说,是按如下过程进行
步骤一、建立H辆电动汽车运行状态的时序模型:
从实际运行统计来看,电动汽车的电池系统和电机系统出现的故障次数较多,占总故障次数的80%左右,因此本发明主要基于电池和电机系统的运行状态建立电动汽车的运行状态时间序列。
步骤1.1、根据H辆电动汽车的动力电池和电机的运行状态,将第m辆电动汽车的运行状态离散为四个状态,分别为和表示第m辆电动汽车动力电池和电机的正常工作状态,表示第m辆电动汽车动力电池的故障状态和电机的正常工作状态、表示第m辆电动汽车动力电池的正常工作状态和电机的故障状态、表示第m辆电动汽车动力电池和电机的故障状态;
步骤1.2、建立第m辆电动汽车的运行状态空间图,如图1所示,并获得运行状态空间图的第m个状态转移矩阵Mm;
在图1中,分别为第m辆电动汽车动力电池和电机的故障率和修复率;
状态转移矩阵
由此可得第m辆电动汽车处于各个运行状态的概率
步骤1.3、根据马尔科夫过程逼近原理,获得第m辆电动汽车四个状态的概率,分别为和表示第m辆电动汽车动力电池和电机的正常工作状态的概率;表示第m辆电动汽车动力电池的故障状态和电机的正常工作状态的概率、表示第m辆电动汽车动力电池的正常工作状态和电机的故障状态的概率、表示第m辆电动汽车动力电池和电机的故障状态的概率;并有
步骤1.4、采用蒙特卡洛方法抽取服从均匀分布的随机数和并利用式(2)确定第m辆电动汽车的运行状态Sm:
利用式(3)抽取第m辆电动汽车运行状态Sm的实际持续时间
式(3)中:Lm表示第m辆电动汽车从运行状态Sm转移到其他三种状态的转移次数,是表示电动汽车从运行状态Sm转移到其他三种状态过程中第h个转移次数时的转移率;
步骤1.5、设蒙特卡洛方法的仿真时间为Y年,在仿真时间Y年以ΔT为时间步长获得采样时刻ti,并有ti=ΔT×i;N表示一年的总小时数;
第m辆电动汽车的运行状态在四个状态中互相转移,从而在仿真时间Y年内形成第m辆电动汽车运行状态的时序模型ψ(m)[t],若采样时刻ti下第m辆电动汽车的运行状态为车动力电池和电机的正常工作状态则ψ(m)[ti]=1,否则ψ(m)[ti]=0;从而获得H辆电动汽车运行状态在仿真时间Y年内的时序模型为:
{ψ(1)[t],ψ(2)[t],…,ψ(m)[t],…,ψ(H)[t]};
通过上面的方法可以建立起电动汽车运行状态模型,该模型在步骤3.4中可集中充电的动力电池数量计算中采用,能得到一个更精确的充电站出力模型。
步骤二、建立H辆电动汽车的时空分布模型:
步骤2.1、为了使模型更加负荷实际,可调度的电动汽车在集中充电时段内应处于闲置状态;令Am事件表示第m辆电动汽车第一次出行开始时刻处在集中充电时段te~tf之间,Bm事件表示第m辆电动汽车最后一次出行结束时刻落在采样时刻ti的一个邻域内,记作则利用式(4)和式(5)分别获得Am事件的概率P(Am)和Bm事件的概率P(Am):
式(4)和式(5)中,fstart(x)和fend(x)分别表示第m辆电动汽车在24小时中第一次出行开始时刻和最后一次出行结束时刻的概率密度函数;fstart(x)和fend(x)是通过美国交通部对全美家用车辆进行了调查,对统计数据进行归一化处理后,用极大似然估计的方法得到:
步骤2.2、在集中充电模式下,在配电网中利用式(6)获得采样时刻ti下第m辆电动汽车的可调度概率
步骤2.3、假设配电网中有W负荷节点,记为集合G={G1,G2,…,Gj,…,GW},Gj表示第j个负荷节点,1≤j≤W;如图2所示,针对家用电动汽车,基于空间负荷预测方法及居民用地和电动汽车空间分布的关联性,利用式(7)建立采样时刻ti下配电网中第j个负荷节点Gj对第m辆动汽车吸引度模型
式(7)中,表示第j个负荷节点Gj对第m辆电动汽车吸引度的权重系数;表示第j个负荷节点Gj在采样时刻ti下的有功负荷;
步骤2.3、利用式(8)计算采样时刻ti下配网中第j个负荷节点Gj所接入的电动汽车分布数量nj(τ):
步骤三、基于dijkstra算法计算各充电站可集中充电的动力电池的总数:
步骤3.1、设配电网中有N座充电站,将所有带有充电站的负荷节点记为集合E={E1,E2,…,Eσ,…EN};Eσ表示第σ个带有充电站的负荷节点;1≤σ≤N<W;
步骤3.2、计算各充电站的负荷节点到配电网中其它W-1个负荷节点的最短路径:
步骤3.2.1、令配电网的W个负荷节点中任意两个负荷节点Ga和Gb之间的距离为l(a,b),1≤a,b≤W;其中,将第σ个带有充电站的负荷节点Eσ到第j个负荷节点Gj的距离记为定义循环次数z;1≤z≤W-1;初始化σ=1;z=1;
其中,任意两点间的距离计算方法为:l(x,y)=ξ×lel,其中,lel为两点间的供电线路长度,ξ为曲折系数,目的是为了将两节点间供电线路长度折算为实际距离;
步骤3.2.2、令集合则集合在集合G中的补集记作利用式(9)获得集合中元素Eσ到补集中任意元素的距离
步骤3.2.3、根据 获得负荷节点并利用 和获得更新集合和更新集合
步骤3.2.4、以负荷节点为中间节点,对任意利用式(10)获得更新的距离
步骤3.2.5、z+1赋值给z,判断z=W是否成立,若成立,则执行步骤3.2.6;否则,返回步骤3.2.3;
步骤3.2.6、σ+1赋值给σ,判断σ=N+1,是否成立,若成立,则执行步骤3.3;否则,返回步骤3.2.2;
步骤3.3、利用式(11)获得第j个负荷节点Gj内的所有电动汽车对第σ个带有充电站的负荷节点Eσ的选择系数为
式(11)中,通过式(11)可以得到各充电站所要服务的配电网节点集合;
步骤3.4、利用式(12)计算采样时刻ti内的第σ个带有充电站的负荷节点Eσ可集中充电的动力电池数量
步骤3.5、利用式(13)获得第σ个带有充电站的负荷节点Eσ采样时刻ti内可集中充电的动力电池的总数
在式(13)中,当集中充电结束后充电站内动力电池将统一配送,从而不影响车主白天的正常使用,因此在ti=tf+ΔT时刻充电站内可集中充电动力电池数量将重新开始累加;
在一天时间周期内,根据步骤三可以得到各时段累计集中的动力电池数量序列,以此作为基础进行迭代更新得到蒙特卡洛仿真时间内各时刻的可集中动力电池数量;
步骤四、建立基于可靠性评估的充电站出力时序模型:
在集中充电模式下充电站可具有双重角色,作为负荷时在配电网中可用于削峰填谷,作为储能电源时可以向孤岛供电,提高配电网可靠性。在配电网无故障状态时,集中的动力电池只有在集中充电时段统一充电,在配电网故障状态下,电动汽车动力电池作为储能装置向配电网供电。
步骤4.1、假设第m辆电动汽车的动力电池在采样时刻ti进入充电站,则利用式(14)和式(15)分别获得第m辆电动汽车动力电池在配电网无故障状态的荷电状态变化特征SOC(ti+ΔT)m和在配电网故障状态下的荷电状态变化特征SOC(ti+ΔT)m:
式(14)和式(15)中:和分别为第m辆电动汽车动力电池的充电效率和放电效率;和分别表示第m辆电动汽车动力电池在采样时刻ti充电功率和放电功率,表示第m辆电动汽车电池的额定容量,SOC(ti)m表示第m辆电动汽车动力电池在采样时刻ti下的荷电状态值;
步骤4.2、假设在任意时刻tc下配电网发生故障,故障时间为ttr;若存在第σ个带有充电站的负荷节点Eσ位于故障区域内,则利用式(16)获得第σ个带有充电站的负荷节点Eσ在采样时刻ti下的时序输出功率从而获得N座充电站在仿真时间Y年内的出力时序模型为:
式(16)中:当SOC(ti)ω>0.2时,ζω(ti)=1,否则ζω(ti)=0;即只有电动汽车的SOC值大于0.2时,才允许参与放电。
集中型充电站作为储能装置时出力具有时变性,结合电动汽车动力电池的出力特点、充电站服务范围和配电网的故障状态,可以得到更为合理的用于可靠性评估的集中型充电站出力模型。
一种用于配电网可靠性评估的电动汽车充电站建模方法的应用是:
在仿真时间Y年内,抽取配电网中电气元件的状态,电气元件的状态分为故障状态和正常状态;以出力时序模型作为孤岛恢复供电电源,采用蒙特卡洛方法对配电网进行可靠性评估,并统计可靠性指标。具体的,用于配电网可靠性评估的电动汽车充电站建模方法的应用具体流程如下:
①输入网络拓扑结构,搜索网络的非电源元件总数为Q,在仿真设置时间内建立各充电站出力时序序列σ=1,2,…,N;
②设置模拟时序初始值为T=0,根据故障概率分布函数抽出系统中各非电源元件无故障工作时间和故障修复时间
③通过比较确定具有最小无故障工作时间的元件q(记为)及其位置,令
④当元件q发生故障时:
I、若元件q在分支线上,则直接统计该负荷点可靠性指标;
II、若元件q在子馈线上则该子馈线脱离主馈线,其它馈线正常运行。判断子馈线内有无充电站,若没有,则根据其停电时间(开关切换时间Tswitch或)直接计算可靠性指标;若有,则在其内部形成孤岛区和非孤岛区。在孤岛区,结合充电站出力时序序列和负荷分级消减策略,计算可靠性指标;在非孤岛区由于没有作为储能装置的充电站,则根据其停电时间(Tswitch或)直接计算可靠性指标;(负荷分级消减策略即依次考虑消减常规负荷中的三类、二类及一类负荷);
III、若元件q在主馈线上,则配电网中子馈线全部脱离主馈线形成孤岛区,位于元件q上游的孤岛运行时间TR=Tswitch,位于其下游的孤岛运行时间在各孤岛区内结合充电站出力时序序列和负荷分级消减策略,计算可靠性指标;主馈线的计算过程同(II);
⑤对元件q重新产生随机数并转化为新的无故障工作时间和修复时间令
⑥判断T是否小于仿真设置时间,如果小于仿真设置时间,返回步骤3继续执行;否则评估结束,统计配电网可靠性指标。
Claims (2)
1.一种用于配电网可靠性评估的电动汽车充电站建模方法,所述配电网含有电动汽车;所述电动汽车包括动力电池和电机并在所述配电网中通过充电站和集中充电模式进行充电;假设所述配电网中含有H辆电动汽车,记为V={V1,V2,…,Vm,…,VH};Vm表示第m辆电动汽车;1≤m≤H;其特征是,所述方法是按如下步骤进行:
步骤一、建立所述H辆电动汽车运行状态的时序模型:
步骤1.1、根据所述H辆电动汽车的动力电池和电机的运行状态,将所述第m辆电动汽车的运行状态离散为四个状态,分别为和 表示所述第m辆电动汽车动力电池和电机的正常工作状态,表示所述第m辆电动汽车动力电池的故障状态和电机的正常工作状态、表示所述第m辆电动汽车动力电池的正常工作状态和电机的故障状态、表示所述第m辆电动汽车动力电池和电机的故障状态;
步骤1.2、建立第m辆电动汽车的运行状态空间图,并获得所述运行状态空间图的第m个状态转移矩阵Mm;
步骤1.3、根据马尔科夫过程逼近原理,获得所述第m辆电动汽车四个状态的概率,分别为和 表示所述第m辆电动汽车动力电池和电机的正常工作状态的概率;表示所述第m辆电动汽车动力电池的故障状态和电机的正常工作状态的概率、表示所述第m辆电动汽车动力电池的正常工作状态和电机的故障状态的概率、表示所述第m辆电动汽车动力电池和电机的故障状态的概率;并有
步骤1.4、采用蒙特卡洛方法抽取服从均匀分布的随机数和并利用式(2)确定所述第m辆电动汽车的运行状态Sm:
利用式(3)抽取所述第m辆电动汽车运行状态Sm的实际持续时间
式(3)中:Lm表示所述第m辆电动汽车从运行状态Sm转移到其他三种状态的转移次数,是表示所述电动汽车从运行状态Sm转移到其他三种状态过程中第h个转移次数时的转移率;
步骤1.5、设所述蒙特卡洛方法的仿真时间为Y年,在所述仿真时间Y年以ΔT为时间步长获得采样时刻ti,并有ti=ΔT×i;N表示一年的总小时数;
所述第m辆电动汽车的运行状态在所述四个状态中互相转移,从而在仿真时间Y年内形成第m辆电动汽车运行状态的时序模型ψ(m)[t],若采样时刻ti下第m辆电动汽车的运行状态为所述电动汽车动力电池和电机的正常工作状态则ψ(m)[ti]=1,否则ψ(m)[ti]=0;从而获得所述H辆电动汽车运行状态在仿真时间Y年内的时序模型为:
{ψ(1)[t],ψ(2)[t],…,ψ(m)[t],…,ψ(H)[t]};
步骤二、建立所述H辆电动汽车的时空分布模型:
步骤2.1、令Am事件表示所述第m辆电动汽车第一次出行开始时刻处在集中充电时段te~tf之间,0<e,Bm事件表示所述第m辆电动汽车最后一次出行结束时刻落在采样时刻ti的一个邻域内,记作则利用式(4)和式(5)分别获得所述Am事件的概率P(Am)和所述Bm事件的概率P(Bm):
式(4)和式(5)中,fstart(x)和fend(x)分别表示所述第m辆电动汽车在24小时中第一次出行开始时刻和最后一次出行结束时刻的概率密度函数;
步骤2.2、在所述集中充电模式下,在配电网中利用式(6)获得采样时刻ti下所述第m辆电动汽车的可调度概率
步骤2.3、假设配电网中有W负荷节点,记为集合G={G1,G2,…,Gj,…,GW},Gj表示第j个负荷节点,1≤j≤W;利用式(7)建立采样时刻ti下配电网中第j个负荷节点Gj对所述第m辆电动汽车吸引度模型
式(7)中,表示第j个负荷节点Gj对所述第m辆电动汽车吸引度的权重系数;表示第j个负荷节点Gj在采样时刻ti下的有功负荷;
步骤2.4、利用式(8)计算采样时刻ti下配网中第j个负荷节点Gj所接入的电动汽车分布数量nj(ti):
步骤三、基于dijkstra算法计算各充电站可集中充电的动力电池的总数:
步骤3.1、设配电网中有N座充电站,将所有带有充电站的负荷节点记为集合E={E1,E2,…,Eσ,…EN};Eσ表示第σ个带有充电站的负荷节点;1≤σ≤N<W;
步骤3.2、计算各充电站的负荷节点到配电网中其它W-1个负荷节点的最短路径:
步骤3.2.1、令所述配电网的W个负荷节点中任意两个负荷节点Ga和Gb之间的距离为l(a,b),1≤a,b≤W;其中,将所述第σ个带有充电站的负荷节点Eσ到第j个负荷节点Gj的距离记为定义循环次数z;1≤z≤W-1;初始化σ=1;z=1;
步骤3.2.2、令集合则集合在集合G中的补集记作利用式(9)获得所述集合中元素Eσ到补集中任意元素的距离
步骤3.2.3、根据获得负荷节点并利用和获得更新集合和更新集合
步骤3.2.4、以负荷节点为中间节点,对任意利用式(10)获得更新的距离
步骤3.2.5、z+1赋值给z,判断z=W是否成立,若成立,则执行步骤3.2.6;否则,返回步骤3.2.3;
步骤3.2.6、σ+1赋值给σ,判断σ=N+1,是否成立,若成立,则执行步骤3.3;否则,返回步骤3.2.2;
步骤3.3、利用式(11)获得第j个负荷节点Gj内的所有电动汽车对第σ个带有充电站的负荷节点Eσ的选择系数为
式(11)中,
步骤3.4、利用式(12)计算采样时刻ti内的第σ个带有充电站的负荷节点Eσ可集中充电的动力电池数量
步骤3.5、利用式(13)获得第σ个带有充电站的负荷节点Eσ采样时刻ti内可集中充电的动力电池的总数
步骤四、建立基于可靠性评估的充电站出力时序模型:
步骤4.1、假设第m辆电动汽车的动力电池在采样时刻ti进入充电站,则利用式(14)和式(15)分别获得第m辆电动汽车动力电池在配电网无故障状态的荷电状态变化特征SOC(ti+ΔT)m和在配电网故障状态下的荷电状态变化特征SOC(ti+ΔT)m:
式(14)和式(15)中:和分别为第m辆电动汽车动力电池的充电效率和放电效率;和分别表示第m辆电动汽车动力电池在采样时刻ti充电功率和放电功率,表示第m辆电动汽车电池的额定容量,SOC(ti)m表示第m辆电动汽车动力电池在采样时刻ti下的荷电状态值;
步骤4.2、假设在任意时刻tc下配电网发生故障,故障时间为ttr;若存在第σ个带有充电站的负荷节点Eσ位于故障区域内,则利用式(16)获得第σ个带有充电站的负荷节点Eσ在采样时刻ti下的时序输出功率从而获得所述N座充电站在仿真时间Y年内的出力时序模型为:
式(16)中:当SOC(ti)ω>0.2时,ζω(ti)=1,否则ζω(ti)=0。
2.一种采用如权利要求1所述的用于配电网可靠性评估的电动汽车充电站建模方法,其特征是:
在所述仿真时间Y年内,抽取配电网中电气元件的状态,所述电气元件的状态分为故障状态和正常状态;以所述出力时序模型作为孤岛恢复供电电源,采用所述蒙特卡洛方法对所述配电网进行可靠性评估,并统计可靠性指标。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510180855.1A CN104734171B (zh) | 2015-04-16 | 2015-04-16 | 一种用于配电网可靠性评估的电动汽车充电站建模方法及其应用 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510180855.1A CN104734171B (zh) | 2015-04-16 | 2015-04-16 | 一种用于配电网可靠性评估的电动汽车充电站建模方法及其应用 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104734171A CN104734171A (zh) | 2015-06-24 |
CN104734171B true CN104734171B (zh) | 2017-01-18 |
Family
ID=53457765
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510180855.1A Active CN104734171B (zh) | 2015-04-16 | 2015-04-16 | 一种用于配电网可靠性评估的电动汽车充电站建模方法及其应用 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104734171B (zh) |
Families Citing this family (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105071508B (zh) * | 2015-09-25 | 2017-02-22 | 武汉电动汽车技术开发有限公司 | 一种电动汽车智能充电方法及其充电系统 |
CN105205502B (zh) * | 2015-10-30 | 2019-01-01 | 山东大学 | 一种基于马尔柯夫蒙特卡罗的负荷特性综合分类方法 |
CN106026076B (zh) * | 2016-05-17 | 2018-06-12 | 上海电力学院 | 一种计及电动汽车支撑能力的用户侧供电可靠性评估方法 |
CN106355290A (zh) * | 2016-09-21 | 2017-01-25 | 四川大学 | 一种基于马尔科夫链的电动汽车充电负荷预测方法及系统 |
CN107176046B (zh) * | 2017-05-10 | 2019-10-18 | 华南理工大学 | 基于充电失败风险排序的电动汽车充放电控制方法 |
CN107292533B (zh) * | 2017-07-12 | 2020-10-23 | 广东工业大学 | 一种四状态空间的一供一备配电网供电可靠性评估方法 |
CN107292532B (zh) * | 2017-07-12 | 2020-11-13 | 广东工业大学 | 一种五状态空间的一供一备配电网供电可靠性评估方法 |
CN109284891B (zh) * | 2018-08-01 | 2021-06-08 | 大连理工大学 | 基于时空索引的充电桩维修调度方法 |
CN109635854B (zh) * | 2018-11-26 | 2023-11-24 | 国网冀北电力有限公司 | 基于马尔科夫链的充电桩故障预测方法及装置 |
CN110210777B (zh) * | 2019-06-11 | 2021-05-04 | 上海电力学院 | 一种含微电网和电动汽车充电站的配电网可靠性评估方法 |
CN110429596B (zh) * | 2019-08-29 | 2021-05-28 | 重庆大学 | 计及电动汽车时空分布的配电网可靠性评估方法 |
CN111369159A (zh) * | 2020-03-11 | 2020-07-03 | 四川大学 | 一种电力交通耦合网络可靠性评估方法 |
CN111682550B (zh) * | 2020-06-18 | 2022-04-01 | 国网山东省电力公司德州供电公司 | 接入高渗透率电源及电动汽车的配电网可靠性评估方法 |
CN111985777A (zh) * | 2020-07-20 | 2020-11-24 | 中国农业大学 | 一种电动汽车负荷聚合体调控能力评估模型建立方法及系统 |
CN112329215B (zh) * | 2020-10-20 | 2024-02-27 | 华北电力大学 | 含电动汽车换电站的配电网的可靠性评估方法及计算设备 |
CN117175587B (zh) * | 2023-11-03 | 2024-03-15 | 国网山东省电力公司东营供电公司 | 考虑柔性负荷的配电网调度优化方法、系统、终端及介质 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102097829A (zh) * | 2009-12-11 | 2011-06-15 | 姚一凡 | 一种利用电动汽车电池储能和供电的分布式电力储能和供电的方法 |
CN102447294A (zh) * | 2010-10-08 | 2012-05-09 | 台达电子工业股份有限公司 | 具有充电效率控制及自适应充电服务功能的车辆充电系统 |
-
2015
- 2015-04-16 CN CN201510180855.1A patent/CN104734171B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN104734171A (zh) | 2015-06-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104734171B (zh) | 一种用于配电网可靠性评估的电动汽车充电站建模方法及其应用 | |
CN102761153B (zh) | 电力监视控制装置及电力监视控制系统 | |
CN103810539B (zh) | 考虑换电服务可用性的电动汽车换电站容量优化配置方法 | |
CN105160451A (zh) | 一种含电动汽车的微电网多目标优化调度方法 | |
CN107704947A (zh) | 一种考虑电动汽车随机接入的微电网多目标优化调度方法 | |
CN103499792B (zh) | 电动汽车动力电池集群可用容量的预测方法 | |
CN111626527B (zh) | 计及可调度电动汽车快/慢充放电形式的智能电网深度学习调度方法 | |
CN107169273A (zh) | 计及延时和v2g充电模式的电动汽车充电功率预测方法 | |
CN109948823B (zh) | 一种光储充电塔自适应鲁棒日前优化调度方法 | |
CN103605079B (zh) | 公交电动汽车及其梯次利用电池集群的v2g可用容量评估方法 | |
CN107634532A (zh) | 基于改进v2g和优先级调度的充放储一体站控制方法 | |
CN103473446A (zh) | 用于有源配电网可靠性评估的负荷削减模型及其实现方法 | |
CN107292516A (zh) | 一种计及负荷分级和能量调度的负荷可靠性评估方法 | |
CN103903090A (zh) | 基于用户意愿和出行规律的电动汽车充电负荷分配方法 | |
CN103246942A (zh) | 风电-电动汽车-火电联合运行模型的多目标调度方法 | |
CN108110790A (zh) | 基于多源数据阶梯式动态模糊系统的主动配电网运行状态评估方法 | |
CN106786977A (zh) | 一种电动汽车充电站的充电调度方法 | |
CN103456099A (zh) | 一种基于实时电价的插入式电动汽车充电控制方法 | |
CN106026076A (zh) | 一种计及电动汽车支撑能力的用户侧供电可靠性评估方法 | |
CN110429596A (zh) | 计及电动汽车时空分布的配电网可靠性评估方法 | |
CN107150601A (zh) | 一种基于多代理商的电动汽车有序充电方法 | |
CN114036459A (zh) | 基于v2g调度响应下的电动汽车用能绿色程度计算方法 | |
CN110571855A (zh) | 含储能设备和ev充电站的园区微电网联合功率响应控制方法 | |
Xue et al. | Adaptive dynamic programming method for optimal battery management of battery electric vehicle | |
CN107622332A (zh) | 一种基于静态安全性约束的电网侧储能容量优化配置方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |