CN110826841A - 考虑用户充电体验和配电网运行风险的充电站规划方法 - Google Patents
考虑用户充电体验和配电网运行风险的充电站规划方法 Download PDFInfo
- Publication number
- CN110826841A CN110826841A CN201910820177.9A CN201910820177A CN110826841A CN 110826841 A CN110826841 A CN 110826841A CN 201910820177 A CN201910820177 A CN 201910820177A CN 110826841 A CN110826841 A CN 110826841A
- Authority
- CN
- China
- Prior art keywords
- charging
- charging station
- point
- formula
- electric automobile
- 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
- 238000007600 charging Methods 0.000 title claims abstract description 583
- 238000013439 planning Methods 0.000 title claims abstract description 88
- 238000000034 method Methods 0.000 title claims abstract description 68
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 57
- 238000005457 optimization Methods 0.000 claims abstract description 42
- 238000010586 diagram Methods 0.000 claims abstract description 41
- 238000010276 construction Methods 0.000 claims abstract description 39
- 238000005192 partition Methods 0.000 claims abstract description 11
- 230000002068 genetic effect Effects 0.000 claims abstract description 10
- 125000004122 cyclic group Chemical group 0.000 claims abstract description 6
- 238000004513 sizing Methods 0.000 claims abstract description 4
- 210000000349 chromosome Anatomy 0.000 claims description 55
- 230000008569 process Effects 0.000 claims description 33
- 238000004364 calculation method Methods 0.000 claims description 13
- 238000012614 Monte-Carlo sampling Methods 0.000 claims description 12
- 108090000623 proteins and genes Proteins 0.000 claims description 12
- 238000005070 sampling Methods 0.000 claims description 8
- FGUUSXIOTUKUDN-IBGZPJMESA-N C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 Chemical compound C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 FGUUSXIOTUKUDN-IBGZPJMESA-N 0.000 claims description 6
- 210000004027 cell Anatomy 0.000 claims description 6
- 238000010277 constant-current charging Methods 0.000 claims description 6
- 238000004088 simulation Methods 0.000 claims description 6
- 230000000694 effects Effects 0.000 claims description 5
- 230000035772 mutation Effects 0.000 claims description 5
- 238000005215 recombination Methods 0.000 claims description 5
- 238000011160 research Methods 0.000 claims description 5
- 108700028369 Alleles Proteins 0.000 claims description 4
- 230000003321 amplification Effects 0.000 claims description 4
- 230000015572 biosynthetic process Effects 0.000 claims description 4
- 238000003199 nucleic acid amplification method Methods 0.000 claims description 4
- 230000035945 sensitivity Effects 0.000 claims description 4
- 238000000638 solvent extraction Methods 0.000 claims description 4
- 238000003786 synthesis reaction Methods 0.000 claims description 4
- 230000002238 attenuated effect Effects 0.000 claims description 3
- 230000008859 change Effects 0.000 claims description 3
- 230000007547 defect Effects 0.000 claims description 3
- 239000000284 extract Substances 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 3
- 230000010287 polarization Effects 0.000 claims description 3
- 230000001172 regenerating effect Effects 0.000 claims description 3
- 238000012502 risk assessment Methods 0.000 claims description 3
- 230000029305 taxis Effects 0.000 claims description 3
- 239000002699 waste material Substances 0.000 claims description 3
- 230000009194 climbing Effects 0.000 description 4
- 230000006386 memory function Effects 0.000 description 4
- 230000006399 behavior Effects 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 230000002452 interceptive effect Effects 0.000 description 2
- 230000006798 recombination Effects 0.000 description 2
- WHXSMMKQMYFTQS-UHFFFAOYSA-N Lithium Chemical compound [Li] WHXSMMKQMYFTQS-UHFFFAOYSA-N 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000007599 discharging Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 229910052744 lithium Inorganic materials 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000010845 search algorithm Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 238000013316 zoning Methods 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/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
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/12—Computing arrangements based on biological models using genetic models
- G06N3/126—Evolutionary algorithms, e.g. genetic algorithms or genetic programming
-
- 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
- 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/0635—Risk analysis of enterprise or organisation activities
-
- 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)
- Economics (AREA)
- Physics & Mathematics (AREA)
- Strategic Management (AREA)
- Theoretical Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Entrepreneurship & Innovation (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Marketing (AREA)
- Tourism & Hospitality (AREA)
- Quality & Reliability (AREA)
- Game Theory and Decision Science (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- Operations Research (AREA)
- Biophysics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Artificial Intelligence (AREA)
- Computational Linguistics (AREA)
- Water Supply & Treatment (AREA)
- Physiology (AREA)
- Genetics & Genomics (AREA)
- Primary Health Care (AREA)
- Biomedical Technology (AREA)
- Public Health (AREA)
- Data Mining & Analysis (AREA)
- Evolutionary Computation (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Charge And Discharge Circuits For Batteries Or The Like (AREA)
Abstract
本发明提供考虑用户充电体验和配电网运行风险的充电站规划方法,包括以下步骤:步骤1、建立电动汽车负荷模型;步骤2、建立充电体验指数,以充电体验指数、充电站建设成本和用户便捷性为充电站选址定容优化目标建立上层综合目标,并提出上层规划约束条件;步骤3、提出计及空间权重的元‑点加权Voronoi图进行服务范围分区;步骤4、基于上层优化目标采用遗传禁忌混合算法建立上层规划流程;步骤5、以新建馈线成本和配电网的运行风险为充电站负荷接入节点优化目标建立下层综合目标;步骤6、基于下层综合目标采用循环替换法建立下层规划流程;本发明能综合用户充电体验、电网运行情况以及充电站建设成本等需求来进行充电站规划。
Description
技术领域
本申请涉及充电站规划领域,特别是涉及一种考虑用户充电体验和配电网运行风险的充电站规划方法。
背景技术
随着电动汽车的迅猛发展,电动汽车充电需求将会迅速增大,建设满足其充电需求的充电站是实现电动汽车快速普及的基础保证。
在实现过程中,发明人发现传统电动汽车充电站规划中至少存在以下问题:
在互联网时代,人们越来越重视生活体验满足感;充电站的规划建设,与电动汽车用户的充电体验密切相关;目前已有研究从充电路上花费的时间、充电前等待时间以及用户行驶到充电站消耗的电能等方面来考虑用户充电体验;此外,充电站规划对用户充电体验的影响,还体现在用户在充电前是否需要排队、能否在电池深度放电前到达充电站等方面;另一方面,不同功能区域的人们,由于时间的紧迫程度不同,对前往充电过程中所花时间的要求不同,而且电动汽车密集程度也不同,因此在评估用户充电体验时应充分考虑上述因素;
在充电站服务范围划分方面,很多研究采用了加权Voronoi图;采用加权Voronoi图时,只考虑充电站的权重,没有考虑不同功能区域对充电站需求的差异;另一方面,在地价低的地方建设容量更大的充电站,建设成本更小,因此充电站权重应考虑地价;再者,以往的研究多以直线距离进行分区,但由于多数道路都是垂直或平行的田字型,采用直角距离可提高分区的准确性;
为了保证在电网安全运行的前提下进行充电站规划,多数研究将电网的运行状态参数(如电压、电流和负荷功率等)作为约束条件;但是,当电动汽车规模较大时,运行越限问题难以避免;同时,由于电动汽车充电的随机性,配电网运行状态不确定性增大,作为风险的基本特征,这使得配电网发生严重越限的可能性增加;此外,目前的研究侧重于充电设施的选址定容优化,很少关注充电站负荷接入节点的优化问题,简单假设充电站负荷就近接入相应负荷节点;但实际上,由于节点的脆弱性不同,充电站负荷接入不同节点带来的运行风险是不同的;因此,通过充电负荷接入节点优化来减小各种运行风险,具有重要的工程意义。
发明内容
基于此,本发明提出计及空间权重的元-点加权Voronoi图进行服务范围分区,建立充电体验指数新指标来衡量充电过程中用户的充电体验要求得到满足的程度,并以充电站的建设成本、用户便捷性和和充电体验指数为目标进行充电站选址定容优化;以新建馈线成本和配电网的运行风险为目标进行充电站负荷接入节点优化;并以此为基础提出考虑用户充电体验和配电网运行风险的充电站互动规划方法。
本发明至少通过如下技术方案之一实现。
考虑用户充电体验和配电网运行风险的充电站规划方法,包括以下步骤:
步骤1、根据电动汽车历史数据,建立电动汽车负荷模型;
步骤2、建立充电体验指数,以充电体验指数、充电站建设成本和用户便捷性为充电站选址定容优化目标建立上层综合目标,并提出上层规划约束条件;
步骤3、提出计及空间权重的元-点加权Voronoi图进行服务范围分区;
步骤4、基于上层优化目标采用遗传禁忌混合算法建立上层规划流程;
步骤5、以新建馈线成本和配电网的运行风险为充电站负荷接入节点优化目标建立下层综合目标,并提出下层约束条件;
步骤6、基于下层综合目标采用循环替换法建立下层规划流程。
进一步地,步骤1所述的根据电动汽车历史数据,建立电动汽车负荷模型包括以下步骤:
步骤1-1、采用电动汽车容量来描述区域内电动汽车的规模,所述电动汽车容量为区域内所有电动汽车额定充电功率的总和,包括处于充电状态和非充电状态的全部电动汽车;将单位面积电动汽车容量定义为电动汽车容量密度;
步骤1-2、根据电动汽车电池实际充电过程,构建电动汽车恒流-恒压二阶段变功率充电模型:
电动汽车在充电过程中,采用一种恒流-恒压二阶段变功率充电模式,当电池端电压小于最大电池端电压Ubmax时,电池是以恒电流模式进行充电;当电池端电压达到最大电池端电压Ubmax时,电池进入恒压阶段,电流迅速衰减为零,充电功率为零;描述其充电过程的方程为:
式中,Ub(t)为t时刻电池两端的电压,ic为充电时的电流,C0是电池的标称容量,R为电池内阻,K0是电池的极化常数,E0是电池的恒定电势,A0和B0分别为电池充电指数段的乘积系数和指数系数,SOC(t0)代表电池在开始充电时的剩余电量,h(t)为t时刻电池的可充电容量,即已放电容量;
在恒流充电过程中,电流i为常数,电动汽车的充电功率表示为:
P(t)=Ub(t)In (2)
式中,In为恒流充电模式下的充电电流,P(t)为t时刻的充电功率;
步骤1-3、根据电动汽车历史出行数据,构建电动汽车充电负荷时序概率模型:
电动汽车的充电负荷模拟须根据电动汽车的类型来选择各抽样环节的概率模型,抽取其行驶距离、到达时间、起始荷电状态和离开时间,从而确定每个时刻电动汽车的行为状态,进而计算相应时刻电动汽车的SOC值,得到其一天的SOC曲线,计算公式为:
(1)充电时
式中,SOC(t)是电动汽车的起始荷电状态;
(2)行驶时
SOC(t+Δt)=SOC(t)-dr/Dr·Δt/Tr (4)
式中,dr是电动汽车的行驶距离,Dr是电动汽车的最大续航里程,Tr是电动汽车的行驶时长,Δt为计算时间间隔;
然后,借助公式(1)和公式(2)算出其一天的充电功率;再将各类电动汽车的负荷逐辆叠加,得到电动汽车充电总负荷;t时刻电动汽车充电总负荷有功功率的计算公式为:
式中:N1、N2和N3分别代表t时刻该地区处于充电状态的电动私家车、公交车和出租车的数量;P1nt、P2nt和P3nt分别代表t时刻第n辆处于充电状态的电动私家车、公交车和出租车的实时充电功率;
步骤1-4、建立最大充电同时率:
为了清楚地了解某个区域电动汽车一天中同时充电的最大比例,定义电动汽车最大充电同时率为一天中电动汽车最大负荷与电动汽车容量的比值,即:
式中,Pt为t时刻区域内电动汽车充电总负荷,Kz为该区域的电动汽车容量。
进一步地,步骤2所述的建立充电体验指数,以充电体验指数、充电站建设成本和用户便捷性为充电站选址定容优化目标建立上层综合目标,并提出上层规划约束条件包括以下步骤:
步骤2-1、综合充电路程用时、充前SOC值和是否需要排队等要求满足程度建立充电体验指数:
定义充电体验指数为充电过程中用户的充电体验要求得到满足的程度,并综合充电路程用时、充前SOC值和是否需要排队等要求的满足程度来反映用户的充电体验指数;
(1)充电路程用时满足率
不同区域的人们由于时间的紧迫程度不同,对充电路上花费时间的要求不同,即充电路程用时要求不同;定义充电路程用时限值为用户对充电路上花费时间的接受极限值;充电过程中花费的时间大于充电路程用时限值时,则不满足充电体验要求;基于此,为考察所有用户的充电路程用时要求的满足程度,可定义充电路程用时满足率为:
式中,ηt_sat为充电路程用时满足率;εd(i)为第i个点的电动汽车容量密度;va为充电过程中电动汽车的平均行驶速度;Td(i)为第i个点的充电路程用时限值,即其所属区域的充电路程用时限值;Mt_sat为满足充电路程用时要求的点的电动汽车容量密度的集合;N0为点的个数;Ld(i)为规划区第i个点到所属充电站的直角距离;
(2)充前SOC值满足率
电动汽车在前往充电前的SOC值满足一定概率分布,当SOC值低于一定值时,将会损害电动汽车电池的寿命,因此可定义充前SOC值满足率为用户到达充电站时SOC值满足不损害电动汽车电池寿命的最低SOC值的概率。通过M次蒙特卡罗抽样得到M个SOC值,则充前SOC值满足率的算法为:
式中,ηs_sat为充前SOC值满足率;pd(i)为第i个点满足充前SOC值要求的概率;εd(i)为第i个点的电动汽车容量密度;Srange为相应电动汽车的最大续航里程;SOC(j)为蒙特卡罗抽样的第j个SOC值;Ωd(i)为第i个点M个SOC值中满足充前SOC值要求的SOC值的集合;card()为集合中元素的个数,SOCmin为不损害电池寿命的最低SOC值;M为蒙特卡罗抽样次数;
(3)无需排队满足率
若电动汽车充电需求大于充电站的配置容量,则用户充电前需要排队,不满足无需排队要求。定义无需排队满足率为用户不需要排队的概率。通过蒙特卡罗抽样得到各充电站M次一天24个时刻的充电站负荷数据,则无需排队满足率的算法为:
式中,ηp_sat为无需排队满足率;pCD(i)为充电站i满足无需排队要求的平均概率;П(j)为充电站i第j次蒙特卡罗抽样一天中充电站负荷满足无需排队要求的时刻的集合;J(i)为充电站i的充电机配置数量;PCD(j,k)为充电站i第j次蒙特卡罗抽样第k个时刻的负荷;Pc为单台充电机的额定功率;αc为充电效率;
(4)用户充电体验指数计算
综合充电路程用时满足率、充前SOC值满足率和无需排队满足率来反映充电体验指数,算法为:
ηg_sat=ηt_satηs_satηp_sat (10)
式中,ηg_sat为充电体验指数;
步骤2-2、建立充电站建设成本:
充电站建设成本包括变压器、充电机和其他设备的建设成本以及土地费用;因此充电站i的建设成本的算法为:
Cz_t(i)=Cs(i)+Cd(i) (11)
式中,Cs(i)为第i个充电站的变压器、充电机和其他设备的建设成本,由充电站的充电机数量决定;Cd(i)为第i个充电站的土地费用,其算法为:
Cd(i)=Sd(i)Dd (12)
式中,Sd(i)为第i个充电站的占用土地面积,由充电站的充电机数量决定,Dd为单位面积土地费用;
则规划区的充电站建设总成本为:
Cg_t=∑Cz_t(i) (13)
步骤2-3、建立用户便捷性指标:
用户便捷性以充电平均距离来表示;在规划区内平均取点,每个点的电动汽车容量密度为其所属区域的电动汽车容量密度,以每个点的电动汽车容量密度为权重对充电距离进行加权平均得到充电平均距离;则规划区的充电平均距离的算法为:
式中,εd(i)和Ld(i)分别为规划区第i个点的电动汽车容量密度和其到所属充电站的直角距离;
步骤2-4、建立上层综合目标:
充电站的上层规划以充电站建设成本、用户便捷性和充电体验指数为目标,对充电站进行选址定容优化;因此,上层综合目标的算法为:
式中,α1、α2、α3为权重系数,其作用为各指标的数量级相同,同时用于消掉各指标的单位,使得上层综合目标为无量纲,B1为上层综合目标;
步骤2-5、建立上层规划约束条件:
(1)相邻站间距离约束
为保证用户的充电需求,充电站的服务范围不能太大,同时为了避免资源浪费,应避免充电站分布过于集中;因此相邻充电站间的直角距离Lcc要满足约束条件:
Lmin<Lcc<Lmax (16)
式中,Lmin和Lmax分别为相邻站间距离约束的上下限;
(2)充电机数约束
为了满足用户的充电需求,同时出于对社会资源的充分利用,保证充电机有较高的利用率,充电站的充电机配置数量应满足充电机数约束:
Jmin≤JCD≤Jmax (17)
式中,Jmin和Jmax分别为最小充电机配置数和最大充电机配置数,JCD为充电站充电机的配置数量,算法为:
式中,ZCD是充电站服务范围内的最大充电需求;KCD为充电站服务范围内的电动汽车容量。
进一步地,步骤3提出计及空间权重的元-点加权Voronoi图进行服务范围分区包括以下步骤:
步骤3-1、提出元-点加权Voronoi图:
针对传统的加权Voronoi图不能反映不同区域权重差异的缺陷,本发明提出元-点加权Voronoi图来进行服务范围划分,以空间权重描述不同区域权重的大小;此时生成元生成Voronoi图的速度不仅与生成元的权重有关,还与各区域的空间权重有关;由于考虑了空间权重,不同的行驶路径两点之间的加权距离可能不同,因此本发明定义点A到点B的有向路径为A向B先纵向行驶,再横向行驶的路径,并定义加权有向距离为有向路径的加权距离;
基于此,元-点加权Voronoi图定义为:
V(Ri)={P∈V(Ri)|D(Ri,P)≤D(Rj,P),j=1,2,…,m1;j≠i} (19)
式中,D(Ri,P)为点Ri到点P的加权有向距离,m1为生成元的个数,D(Rj,P)为点Rj到点P的加权有向距离,若Ri的权重为vi,Ri到点P的有向路径经过区域Q1,Q2,…,Qn,n为Ri到点P的有向路径经过的区域个数,相应的空间权重为θ1,θ2,…,θn,Ri到点P的有向路径经过的边界点分别为P1,P2,…,Pn-1,d(Pn-1,P)为点Pn-1到点P的有向距离,则Ri到点P的加权有向距离为:
步骤3-2、建立充电站权重:
在地价低的地方建设服务范围更大、容量更大的充电站,能减小充电站建设总成本;因此在服务范围划分时,可设充电站权重的平方与其占用土地单价的倒数成反比,则充电站i的权重为:
式中,Di为充电站i的占用土地单价;
步骤3-3、建立空间权重:
本发明综合不同类型区域的充电路程用时要求和电动汽车容量密度来建立区域的空间权重,则区域类型i的空间权重的算法为:
式中,θq(i)、εq(i)、Tq(i)分别为区域类型i的空间权重、电动汽车容量密度和充电路程用时限值,β1和β2为权重系数;
步骤3-4、充电站服务范围划分
为便于计算,在利用元-点加权Voronoi图进行充电站规划时,规划区内均匀生成点,邻近点的移动步长为Δx=Δy=Δl,用矩阵分别记录每个点的坐标、电动汽车容量密度、土地单价、充电路程用时限值等信息,建立点数据库;在画Voronoi图时,每个点都用一个方块来代表,方块的边长为移动步长,点位于方块的中心;方块的电动汽车容量为方块的面积与点的电动汽车容量密度的乘积;
若权重为vi的充电站Ri到点O的有向路径依次经过O1,O2,…,On-1,各点相应的空间权重为θq1,θq2,…,θqn,则Ri到点O的加权有向距离为:
将各生成元到某点的加权有向距离进行比较,加权有向距离最小对应的生成元将该点划进其服务范围;
充电站的电动汽车容量为其服务范围内所有点的电动汽车容量之和,算法为:
式中,εd(i)为该充电站服务范围内第i个点的电动汽车容量密度,ΔS为单个方块的面积。
进一步地,步骤4所述的基于上层优化目标采用遗传禁忌混合算法建立上层规划流程包括以下内容:
本发明采用遗传禁忌混合算法(GATS)进行充电站上层规划;GATS综合了GA(遗传算法)具有多出发点和TS(禁忌搜索算法)的记忆功能和爬山能力强的特点,把TS独有的记忆功能引入到GA进化搜索过程中,构造新的重组算子(TSR);使用TS改进GA的爬山能力,即把TS作为GA的变异算子(TSM);充电站上层规划主要步骤如下:
Ⅰ.采集点数据库信息,作为优化的基础数据;
Ⅱ.利用点数据库的电动汽车容量密度信息,求出规划区内的电动汽车容量,即Kz,进而计算规划区的最大充电需求Zz,算法为:
Ⅲ.根据规划区内的最大充电需求Zz,以充电站的最小充电机配置数Jmin和最大充电机配置数Jmax为约束,计算充电站的最小座数Nmin和最大座数Nmax,即:
Ⅳ.对于不同的充电站数方案,基于上层综合目标采用GATS和元-点加权Voronoi图进行选址定容优化,包括以下步骤:
(1)记充电站数为g,设置GA的最大迭代数Ngen、禁忌表长度Tlength1、群体规模Npop、重组概率pch和变异概率pm;TSM中禁忌表长度Tlength2、邻域解个数Nh和最大迭代数Nt等参数;
(2)tr表示迭代次数,令tr=0,将规划区平均分为g个区,分别从各区中随机抽取一个地址,得到染色体的地址序列[x1,y1;…;xi,yi;…;xg,yg],(xi,yi)为染色体的地址序列中的第i个地址,地址之间不满足站间距离约束的重新生成;
(3)计算当代群体的上层综合目标;
(4)按滚轮盘方式选出Npop个染色体;
(5)生成0和1之间的随机数ri,i=1,2,…,Npop,若ri<pch,则第i个染色体为父代;对每对染色体进行交叉操作,随机互换两个基因,判断子代是否满足站间距离约束,不满足则选最优父代进下一代;若满足,则将子代的上层综合目标同渴望水平相比,若更好,则破禁进下一代;若更差,但不属于禁忌,则接受这个子代,若属于禁忌,则选最优父代进下一代;
(6)生成0和1之间的随机数qi,i=1,2,…,Npop;若qi<pm,则对第i个染色体进行变异操作:
a.以当前染色体作为最优染色体,将其上层综合目标作为最优上层综合目标,置TSM的禁忌表为空;
b.以当前地址基因为中心生成邻域地址基因(x+Δx,y)、(x-Δx,y)、(x,y+Δy)和(x,y-Δy),加上当前地址基因构成等位基因序列,进行全排列得到5^g个染色体并从中随机抽出Nh个;
c.将这Nh个染色体中不满足站间距离约束的剔除;剩余的不在禁忌表中的染色体为候选染色体,计算其中分区后满足充电机数约束的上层综合目标;
d.以上层综合目标最小选出最优候选染色体,并用它替代当前染色体和最早进入不满足的染色体;
e.判断最优候选染色体的上层综合目标是否小于最优上层综合目标?是则用它替换最优染色体,并将其上层综合目标替换最优上层综合目标,然后转步骤f;否则直接转到步骤f;
f.判断是否达到最大迭代数Nt,或最优上层综合目标已连续30次保持不变?是则输出最优染色体;否则返回步骤b;
(7)tr=tr+1,若tr<Ngen,转步骤(3);否则输出最优解,终止算法。
进一步地,步骤5所述的以新建馈线成本和配电网的运行风险为充电站负荷接入节点优化目标建立下层综合目标,并提出下层约束条件包括以下步骤:
步骤5-1、计算新建馈线成本:
充电站负荷要接入配电网,要通过建设新的馈线来接入配电网的节点;假定新建馈线是沿直线铺设的,则充电站的新建馈线成本算法为:
Ck=Dk∑Lk(i) (27)
式中,Dk为馈线的单位长度造价;Lk(i)为第i个充电站的新建馈线长度,即第i个充电站到接入节点的直角距离;
步骤5-2、计算配电网运行风险:
(1)元件风险指标
基于加权熵原理,以加权电压分布熵表征电压越限及电压分布不确定性引起的电压风险,以加权潮流分布熵表征潮流越限及潮流分布不确定性引起的潮流风险,算法为:
式中,Hv_i(t)和Hl_i(t)分别为t时刻节点i的加权电压分布熵和支路i的加权潮流分布熵;nv_i(t)和nl_i(t)分别为t时刻节点i的电压状态数和支路i的潮流状态数;Sv_j(t)和Pv_j(t)分别为t时刻节点i第j个电压状态的电压损失严重度和概率;Sl_j(t)和Pl_j(t)分别为t时刻支路i第j个潮流状态的潮流损失严重度和概率;
损失严重度S的算法为:
式中,μ为损失严重度的放大系数,以调节损失严重度的灵敏性;w为损失量;
(2)网络风险指标
本发明综合电压风险和潮流风险对含大规模电动汽车的配电网进行运行风险评估;为考察配电网运行风险一天的变化情况,可将t时刻网络运行风险定义为:
Hw(t)=γ1Hv(t)+γ2Hl(t) (30)
式中,Hw(t)为t时刻的网络运行风险;γ1和γ2为权重系数;Hv(t)和Hl(t)分别为t时刻的网络电压风险和网络潮流风险,算法为:
式中,z1和z2分别为配电网的节点数和支路数;
为研究配电网安全运行面对的最严重的风险情况,定义网络最大运行风险为网络运行风险一天的最大值,即:
Hwmax=max(Hw(t))t=1,2,…,24 (32)
步骤5-3、提出下层综合目标:
充电站下层规划综合考虑新建馈线成本和配电网的运行风险,对充电站负荷接入节点进行优化;定义下层综合目标为:
B2=ρ1Ck+ρ2Hwmax (33)
式中,Ck和Hwmax分别为新建馈线成本和网络最大运行风险;ρ1和ρ2为权重系数;
步骤5-4、提出下层规划约束条件:
(1)馈线长度约束
当新建馈线太长时,不仅馈线成本过大,同时馈线的阻抗过大,会增大功率损耗和电压损耗。因此新建的馈线应满足馈线长度约束:
Lk(i)≤Lkmax (34)
式中,Lkmax为新建馈线的长度限值;
(2)充电站负荷约束
由于充电站的充电负荷不可能大于其安装容量,因此要满足约束:
PCD(i)≤αcJ(i)Pc (35)
式中,PCD(i)为充电站i的充电负荷。
进一步地,步骤6所述的基于下层综合目标采用循环替换法建立下层规划流程包括以下内容:
(1)输入充电站地址及各分区电动汽车容量、充电过程参数、配电网络参数数据;
(2)通过蒙特卡罗抽样计算各个充电站分区的Nm次日充电负荷数据;
(3)将充电站负荷就近接入配电网节点作为当前方案,计算下层综合目标;
(4)记充电站的个数为g,确定各分区的优化顺序,bm=0;
(5)bm=bm+1,计算充电站bm到各节点的距离,将满足约束的节点作为候选节点;
(6)保持其他接入节点不变,将充电站bm的负荷分别接入各候选节点,计算下层综合目标,以下层综合目标最小选最优候选方案;
(7)判断最优候选方案与当前方案的优劣,若最优候选方案比当前方案更优,则以最优候选方案替换当前方案,否则保留当前方案不变;
(8)判断是否有bm<g,是则返回步骤(5),否则转到步骤(9);
(9)不断循环优化,直至充电站负荷接入方案不再改变为止。
与现有技术相比,本发明考虑用户充电体验和配电网运行风险的充电站互动规划方法具有如下显著的效果:本发明提出的计及空间权重的元-点加权Voronoi图能够综合考虑占用土地价格、充电路程用时要求和电动汽车容量密度,使充电站尽量选择低价低的位置作为站址,同时能够反映不同区域对充电站需求的差异,使充电站更加靠近商业区和办公区等对充电要求更高的区域,更能反映用户对充电站的实际要求;考虑用户的充电体验,将充电体验指数作为充电站选址定容的优化目标,有效地改善用户的充电体验,从而更好地满足用户的充电体验需求;将配电网的运行风险作为充电站负荷接入节点优化目标,以较小的新建馈线成本代价有效减小大规模充电站负荷接入配电网造成的运行风险,保证配电网能够更加安全稳定地运行,具有实际意义。
附图说明
图1为本发明实施例中的基于循环替换法的下层规划流程图;
图2为本发明实施例中规划区的区域分类及配网图;
图3为本发明实施例中各类型常规负荷的曲线图;
图4为本发明实施例中充电站数为6时最优规划方案的充电站服务范围分区图。
具体实施方式
为了使本申请的目的、技术方案及优点更加清楚明白,以下通过实施例,对本申请进行进一步详细说明。应当理解,此处描述的具体实施例仅仅用以解释本申请,并不用于限定本申请。
本发明实施例将给出考虑用户充电体验和配电网运行风险的充电站互动规划方法的具体应用实例,包括以下步骤:
步骤1:根据电动汽车历史数据,建立电动汽车负荷模型包括以下步骤:
步骤1-1、采用电动汽车容量来描述区域内电动汽车的规模,所述电动汽车容量为区域内所有电动汽车额定充电功率的总和,包括处于充电状态和非充电状态的全部电动汽车;将单位面积电动汽车容量定义为电动汽车容量密度;
步骤1-2、根据电动汽车电池实际充电过程,构建电动汽车恒流-恒压二阶段变功率充电模型:
电动汽车在充电过程中,采用一种恒流-恒压二阶段变功率充电模式,当电池端电压小于最大电池端电压Ubmax时,电池是以恒电流模式进行充电;当电池端电压达到最大电池端电压Ubmax时,电池进入恒压阶段,电流迅速衰减为零,充电功率为零;描述其充电过程的方程为:
式中,Ub(t)为t时刻电池两端的电压,ic为充电时的电流,C0是电池的标称容量,R为电池内阻,K0是电池的极化常数,E0是电池的恒定电势,A0和B0分别为电池充电指数段的乘积系数和指数系数,SOC(t0)代表电池在开始充电时的剩余电量,h(t)为t时刻电池的可充电容量,即已放电容量;
在恒流充电过程中,电流i为常数,电动汽车的充电功率表示为:
P(t)=Ub(t)In (2)
式中,In为恒流充电模式下的充电电流,P(t)为t时刻的充电功率;
步骤1-3、根据电动汽车历史出行数据,构建电动汽车充电负荷时序概率模型:
电动汽车的充电负荷模拟须根据电动汽车的类型来选择各抽样环节的概率模型,抽取其行驶距离、到达时间、起始荷电状态和离开时间,从而确定每个时刻电动汽车的行为状态,进而计算相应时刻电动汽车的SOC值,得到其一天的SOC曲线,计算公式为:
(1)充电时
式中,SOC(t)是电动汽车的起始荷电状态;
(2)行驶时
SOC(t+Δt)=SOC(t)-dr/Dr·Δt/Tr (4)
式中,dr是电动汽车的行驶距离,Dr是电动汽车的最大续航里程,Tr是电动汽车的行驶时长,Δt为时间间隔;
然后,借助公式(1)和公式(2)算出其一天的充电功率;再将各类电动汽车的负荷逐辆叠加,得到电动汽车充电总负荷;t时刻电动汽车充电总负荷有功功率的计算公式为:
式中:N1、N2和N3分别代表t时刻该地区处于充电状态的电动私家车、公交车和出租车的数量;P1nt、P2nt和P3nt分别代表t时刻第n辆处于充电状态的电动私家车、公交车和出租车的实时充电功率;
步骤1-4、建立最大充电同时率:
为了清楚地了解某个区域电动汽车一天中同时充电的最大比例,定义电动汽车最大充电同时率为一天中电动汽车最大负荷与电动汽车容量的比值,即:
式中,Pt为t时刻区域内电动汽车充电总负荷,Kz为该区域的电动汽车容量。
步骤2:建立充电体验指数,以充电体验指数、充电站建设成本和用户便捷性为充电站选址定容优化目标建立上层综合目标,并提出上层规划约束条件包括以下步骤:
步骤2-1、综合充电路程用时、充前SOC值和是否需要排队等要求满足程度建立充电体验指数:
定义充电体验指数为充电过程中用户的充电体验要求得到满足的程度,并综合充电路程用时、充前SOC值和是否需要排队等要求的满足程度来反映用户的充电体验指数;
(1)充电路程用时满足率
不同区域的人们由于时间的紧迫程度不同,对充电路上花费时间的要求不同,即充电路程用时要求不同;定义充电路程用时限值为用户对充电路上花费时间的接受极限值;充电过程中花费的时间大于充电路程用时限值时,则不满足充电体验要求;基于此,为考察所有用户的充电路程用时要求的满足程度,可定义充电路程用时满足率为:
式中,ηt_sat为充电路程用时满足率;εd(i)为第i个点的电动汽车容量密度;va为充电过程中电动汽车的平均行驶速度;Td(i)为第i个点的充电路程用时限值,即其所属区域的充电路程用时限值;Mt_sat为满足充电路程用时要求的点的电动汽车容量密度的集合;N0为点的个数;Ld(i)为规划区第i个点到所属充电站的直角距离;
(2)充前SOC值满足率
电动汽车在前往充电前的SOC值满足一定概率分布,当SOC值低于一定值时,将会损害电动汽车电池的寿命,因此可定义充前SOC值满足率为用户到达充电站时SOC值满足不损害电动汽车电池寿命的最低SOC值的概率。通过M次蒙特卡罗抽样得到M个SOC值,则充前SOC值满足率的算法为:
式中,ηs_sat为充前SOC值满足率;pd(i)为第i个点满足充前SOC值要求的概率;εd(i)为第i个点的电动汽车容量密度;Srange为相应电动汽车的最大续航里程;SOC(j)为蒙特卡罗抽样的第j个SOC值;Ωd(i)为第i个点M个SOC值中满足充前SOC值要求的SOC值的集合;card()为集合中元素的个数,SOCmin为不损害电池寿命的最低SOC值;M为蒙特卡罗抽样次数;
(3)无需排队满足率
若电动汽车充电需求大于充电站的配置容量,则用户充电前需要排队,不满足无需排队要求。定义无需排队满足率为用户不需要排队的概率。通过蒙特卡罗抽样得到各充电站M次一天24个时刻的充电站负荷数据,则无需排队满足率的算法为:
式中,ηp_sat为无需排队满足率;pCD(i)为充电站i满足无需排队要求的平均概率;П(j)为充电站i第j次蒙特卡罗抽样一天中充电站负荷满足无需排队要求的时刻的集合;J(i)为充电站i的充电机配置数量;PCD(j,k)为充电站i第j次蒙特卡罗抽样第k个时刻的负荷;Pc为单台充电机的额定功率;αc为充电效率;
(4)用户充电体验指数计算
综合充电路程用时满足率、充前SOC值满足率和无需排队满足率来反映充电体验指数,算法为:
ηg_sat=ηt_satηs_satηp_sat (10)
式中,ηg_sat为充电体验指数;
步骤2-2、建立充电站建设成本:
充电站建设成本包括变压器、充电机和其他设备的建设成本以及土地费用;因此充电站i的建设成本的算法为:
Cz_t(i)=Cs(i)+Cd(i) (11)
式中,Cs(i)为第i个充电站的变压器、充电机和其他设备的建设成本,由充电站的充电机数量决定;Cd(i)为第i个充电站的土地费用,其算法为:
Cd(i)=Sd(i)Dd (12)
式中,Sd(i)为第i个充电站的占用土地面积,由充电站的充电机数量决定,Dd为单位面积土地费用;
则规划区的充电站建设总成本为:
Cg_t=∑Cz_t(i) (13)
步骤2-3、建立用户便捷性指标:
用户便捷性以充电平均距离来表示;在规划区内平均取点,每个点的电动汽车容量密度为其所属区域的电动汽车容量密度,以每个点的电动汽车容量密度为权重对充电距离进行加权平均得到充电平均距离;则规划区的充电平均距离的算法为:
式中,εd(i)和Ld(i)分别为规划区第i个点的电动汽车容量密度和其到所属充电站的直角距离;
步骤2-4、建立上层综合目标:
充电站的上层规划以充电站建设成本、用户便捷性和充电体验指数为目标,对充电站进行选址定容优化;因此,上层综合目标的算法为:
式中,α1、α2、α3为权重系数,其作用为各指标的数量级相同,同时用于消掉各指标的单位,使得上层综合目标为无量纲,B1为上层综合目标;
步骤2-5、建立上层规划约束条件:
(1)相邻站间距离约束
为保证用户的充电需求,充电站的服务范围不能太大,同时为了避免资源浪费,应避免充电站分布过于集中;因此相邻充电站间的直角距离Lcc要满足约束条件:
Lmin<Lcc<Lmax (16)
式中,Lmin和Lmax分别为相邻站间距离约束的上下限;
(2)充电机数约束
为了满足用户的充电需求,同时出于对社会资源的充分利用,保证充电机有较高的利用率,充电站的充电机配置数量应满足充电机数约束:
Jmin≤JCD≤Jmax (17)
式中,Jmin和Jmax分别为最小充电机配置数和最大充电机配置数,JCD为充电站充电机的配置数量,算法为:
式中,ZCD是充电站服务范围内的最大充电需求;KCD为充电站服务范围内的电动汽车容量。
步骤3:提出计及空间权重的元-点加权Voronoi图进行服务范围分区包括以下步骤:
步骤3-1、提出元-点加权Voronoi图:
针对传统的加权Voronoi图不能反映不同区域权重差异的缺陷,本发明提出元-点加权Voronoi图来进行服务范围划分,以空间权重描述不同区域权重的大小;此时生成元生成Voronoi图的速度不仅与生成元的权重有关,还与各区域的空间权重有关;由于考虑了空间权重,不同的行驶路径两点之间的加权距离可能不同,因此本发明定义点A到点B的有向路径为A向B先纵向行驶,再横向行驶的路径,并定义加权有向距离为有向路径的加权距离;
基于此,元-点加权Voronoi图定义为:
V(Ri)={P∈V(Ri)|D(Ri,P)≤D(Rj,P),j=1,2,…,m1;j≠i} (19)
式中,D(Ri,P)为点Ri到点P的加权有向距离,m1为生成元的个数,D(Rj,P)为点Rj到点P的加权有向距离,若Ri的权重为vi,Ri到点P的有向路径经过区域Q1,Q2,…,Qn,n为Ri到点P的有向路径经过的区域个数,相应的空间权重为θ1,θ2,…,θn,Ri到点P的有向路径经过的边界点分别为P1,P2,…,Pn-1,d(Pn-1,P)为点Pn-1到点P的有向距离,则Ri到点P的加权有向距离为:
步骤3-2、建立充电站权重:
在地价低的地方建设服务范围更大、容量更大的充电站,能减小充电站建设总成本;因此在服务范围划分时,可设充电站权重的平方与其占用土地单价的倒数成反比,则充电站i的权重为:
式中,Di为充电站i的占用土地单价;
步骤3-3、建立空间权重:
本发明综合不同类型区域的充电路程用时要求和电动汽车容量密度来建立区域的空间权重,则区域类型i的空间权重的算法为:
式中,θq(i)、εq(i)、Tq(i)分别为区域类型i的空间权重、电动汽车容量密度和充电路程用时限值,β1和β2为权重系数;
步骤3-4、充电站服务范围划分
为便于计算,在利用元-点加权Voronoi图进行充电站规划时,规划区内均匀生成点,邻近点的移动步长为Δx=Δy=Δl,用矩阵分别记录每个点的坐标、电动汽车容量密度、土地单价、充电路程用时限值等信息,建立点数据库;在画Voronoi图时,每个点都用一个方块来代表,方块的边长为移动步长,点位于方块的中心;方块的电动汽车容量为方块的面积与点的电动汽车容量密度的乘积;
若权重为vi的充电站Ri到点O的有向路径依次经过O1,O2,…,On-1,各点相应的空间权重为θq1,θq2,…,θqn,则Ri到点O的加权有向距离为:
将各生成元到某点的加权有向距离进行比较,加权有向距离最小对应的生成元将该点划进其服务范围;
充电站的电动汽车容量为其服务范围内所有点的电动汽车容量之和,算法为:
式中,εd(i)为该充电站服务范围内第i个点的电动汽车容量密度,ΔS为单个方块的面积。
步骤4:基于上层优化目标采用遗传禁忌混合算法建立上层规划流程包括以下内容:
本发明采用遗传禁忌混合算法(GATS)进行充电站上层规划;GATS综合了GA具有多出发点和TS的记忆功能和爬山能力强的特点,把TS独有的记忆功能引入到GA进化搜索过程中,构造新的重组算子(TSR);使用TS改进GA的爬山能力,即把TS作为GA的变异算子(TSM);充电站上层规划主要步骤如下:
Ⅰ.采集点数据库信息,作为优化的基础数据;
Ⅱ.利用点数据库的电动汽车容量密度信息,求出规划区内的电动汽车容量,即Kz,进而计算规划区的最大充电需求Zz,算法为:
Ⅲ.根据规划区内的最大充电需求Zz,以充电站的最小充电机配置数Jmin和最大充电机配置数Jmax为约束,计算充电站的最小座数Nmin和最大座数Nmax,即:
Ⅳ.对于不同的充电站数方案,基于上层综合目标采用GATS和元-点加权Voronoi图进行选址定容优化,包括以下步骤:
(1)记充电站数为g,设置GA的最大迭代数Ngen、禁忌表长度Tlength1、群体规模Npop、重组概率pch和变异概率pm;TSM中禁忌表长度Tlength2、邻域解个数Nh和最大迭代数Nt等参数;
(2)tr表示迭代次数,令tr=0,将规划区平均分为g个区,分别从各区中随机抽取一个地址,得到染色体的地址序列[x1,y1;…;xi,yi;…;xg,yg],(xi,yi)为染色体的地址序列中的第i个地址,地址之间不满足站间距离约束的重新生成;
(3)计算当代群体的上层综合目标;
(4)按滚轮盘方式选出Npop个染色体;
(5)生成0,1之间的随机数ri,i=1,2,…,Npop,若ri<pch,则第i个染色体为父代;对每对染色体进行交叉操作,随机互换两个基因,判断子代是否满足站间距离约束,不满足则选最优父代进下一代;若满足,则将子代的上层综合目标同渴望水平相比,若更好,则破禁进下一代;若更差,但不属于禁忌,则接受这个子代,若属于禁忌,则选最优父代进下一代;
(6)生成0,1之间的随机数qi,i=1,2,…,Npop;若qi<pm,则对第i个染色体进行变异操作:
a.以当前染色体作为最优染色体,将其上层综合目标作为最优上层综合目标,置TSM的禁忌表为空;
b.以当前地址基因为中心生成邻域地址基因(x+Δx,y)、(x-Δx,y)、(x,y+Δy)和(x,y-Δy),加上当前地址基因构成等位基因序列,进行全排列得到5^g个染色体并从中随机抽出Nh个;
c.将这Nh个染色体中不满足站间距离约束的剔除;剩余的不在禁忌表中的染色体为候选染色体,计算其中分区后满足充电机数约束的上层综合目标;
d.以上层综合目标最小选出最优候选染色体,并用它替代当前染色体和最早进入不满足的染色体;
e.判断最优候选染色体的上层综合目标是否小于最优上层综合目标?是则用它替换最优染色体,并将其上层综合目标替换最优上层综合目标,然后转步骤f;否则直接转到步骤f;
f.判断是否达到最大迭代数Nt,或最优上层综合目标已连续30次保持不变?是则输出最优染色体;否则返回步骤b;
(7)tr=tr+1,若tr<Ngen,转步骤(3);否则输出最优解,终止算法。
步骤5:以新建馈线成本和配电网的运行风险为充电站负荷接入节点优化目标建立下层综合目标,并提出下层约束条件包括以下步骤:
步骤5-1、计算新建馈线成本:
充电站负荷要接入配电网,要通过建设新的馈线来接入配电网的节点;假定新建馈线是沿直线铺设的,则充电站的新建馈线成本算法为:
Ck=Dk∑Lk(i) (27)
式中,Dk为馈线的单位长度造价;Lk(i)为第i个充电站的新建馈线长度,即第i个充电站到接入节点的直角距离;
步骤5-2、计算配电网运行风险:
(1)元件风险指标
基于加权熵原理,以加权电压分布熵表征电压越限及电压分布不确定性引起的电压风险,以加权潮流分布熵表征潮流越限及潮流分布不确定性引起的潮流风险,算法为:
式中,Hv_i(t)和Hl_i(t)分别为t时刻节点i的加权电压分布熵和支路i的加权潮流分布熵;nv_i(t)和nl_i(t)分别为t时刻节点i的电压状态数和支路i的潮流状态数;Sv_j(t)和Pv_j(t)分别为t时刻节点i第j个电压状态的电压损失严重度和概率;Sl_j(t)和Pl_j(t)分别为t时刻支路i第j个潮流状态的潮流损失严重度和概率;
损失严重度S的算法为:
式中,μ为损失严重度的放大系数,以调节损失严重度的灵敏性;w为损失量;
(2)网络风险指标
本发明综合电压风险和潮流风险对含大规模电动汽车的配电网进行运行风险评估;为考察配电网运行风险一天的变化情况,可将t时刻网络运行风险定义为:
Hw(t)=γ1Hv(t)+γ2Hl(t) (30)
式中,Hw(t)为t时刻的网络运行风险;γ1和γ2为权重系数;Hv(t)和Hl(t)分别为t时刻的网络电压风险和网络潮流风险,算法为:
式中,z1和z2分别为配电网的节点数和支路数;
为研究配电网安全运行面对的最严重的风险情况,定义网络最大运行风险为网络运行风险一天的最大值,即:
Hwmax=max(Hw(t))t=1,2,…,24 (32)
步骤5-3、提出下层综合目标:
充电站下层规划综合考虑新建馈线成本和配电网的运行风险,对充电站负荷接入节点进行优化;定义下层综合目标为:
B2=ρ1Ck+ρ2Hwmax (33)
式中,Ck和Hwmax分别为新建馈线成本和网络最大运行风险;ρ1和ρ2为权重系数;
步骤5-4、提出下层规划约束条件:
(1)馈线长度约束
当新建馈线太长时,不仅馈线成本过大,同时馈线的阻抗过大,会增大功率损耗和电压损耗。因此新建的馈线应满足馈线长度约束:
Lk(i)≤Lkmax (34)
式中,Lkmax为新建馈线的长度限值;
(2)充电站负荷约束
由于充电站的充电负荷不可能大于其安装容量,因此要满足约束:
PCD(i)≤αcJ(i)Pc (35)
式中,PCD(i)为充电站i的充电负荷。
步骤6:基于下层综合目标采用循环替换法建立下层规划流程如图1所示,包括以下内容:
(1)输入充电站地址及各分区电动汽车容量、充电过程参数、配电网络参数等数据;
(2)通过蒙特卡罗抽样计算各个充电站分区的Nm次日充电负荷数据;
(3)将充电站负荷就近接入配电网节点作为当前方案,计算下层综合目标;
(4)记充电站的个数为g,确定各分区的优化顺序,bm=0;
(5)bm=bm+1,计算充电站bm到各节点的距离,将满足约束的节点作为候选节点;
(6)保持其他接入节点不变,将充电站bm的负荷分别接入各候选节点,计算下层综合目标,以下层综合目标最小选最优候选方案;
(7)判断最优候选方案与当前方案的优劣,若最优候选方案比当前方案更优,则以最优候选方案替换当前方案,否则保留当前方案不变;
(8)判断是否有bm<g,是则返回步骤(5),否则转到步骤(9);
(9)不断循环优化,直至充电站负荷接入方案不再改变为止。
在实施例中,规划区如图2所示,由10kV配电系统供电。节点es为平衡节点,电压设为1.05p.u.,系统等值电抗Xs为2Ω。配电网线路线型为LGJ-185,线路参数如表1所示。e0为不接入负荷的联络节点。图2中1-30为配电网的负荷节点,负荷节点位置及负荷峰值如表2所示。电动汽车负荷功率因数为0.95。本发明设单台充电机功率为120kW,最大充电同时率λmax取0.2,SOCmin取0.15,va取20km/h。各类型常规负荷曲线如图3所示。站间距离约束上下限Lmax和Lmin分别取6km和3km。新建馈线的长度限值Lkmax取4km。
为使损失严重度具有较大的灵敏性,令损失量增大0.001时损失严重度最小的增量为0.05,此时放大系数μ的值为82.422。蒙特卡罗抽样次数M设为500次。权重系数α1、α2、α3分别取0.02、0.49和0.49;β1、β2、γ1、γ2皆取0.5。ρ1和ρ2分别取0.9和0.1。
GATS参数设为:GA中Npop为200,Ngen为50,pch为0.5,Tlength1为15,pm为0.2;TS中Nh为30,Tlength2为15,Nt为50。
其他仿真参数设置如表1-表6。
表1线路参数
表2负荷节点位置
表3各区域参数设置
表4电动汽车参数
表5锂电池充电过程参数
表6充电站参数
为证明本发明提出的充电站规划方法的有效性,本文进行了以下仿真。
采用GATS进行充电站规划,只有在充电站数为4、5和6时存在满足约束条件的规划方案,不同充电站数下的最优规划方案的仿真结果如表7所示。
表7不同充电站数最优方案
由表7可知,随着建设充电站数目的增多,充电站建设总成本增大,但是充电平均距离减小,用户充电体验指数增大。因此,用户便捷性和充电体验指数的提高,是以建设成本的提高为代价的。当充电站数为6时,上层综合目标最小,能够同时较好地满足建设成本、用户便捷性和充电体验的需求,相应的最优规划方案为充电站选址定容优化的最优方案。
充电站数为6时的最优规划方案的服务范围分区如图4所示,各指标数据如表8所示。
表8充电站站址与成本
由图4和表8可知,充电站选址时会优先选择工业区或休闲区这些地价较低的地方作为站址,以降低建设成本。位于工业区的充电站3和6由于地价低,充电站权重大,因此不仅建设成本较小,而且服务范围较大。因此验证了本发明在充电站服务范围分区时以占用土地价格来衡量充电站权重的合理性。
充电站数为6时,分别采用加权Voronoi图分区(情景1)和本发明提出的元-点加权Voronoi图分区(情景2)时的最优规划方案中,各区域的点到所属充电站的充电平均距离如表9所示。
表9各区域充电平均距离
由表9可知,相比于加权Voronoi图分区,采用考虑空间权重的元-点加权Voronoi图进行分区后,商业区和办公区的充电平均距离明显减小了,而工业区和休闲区的充电平均距离增大了。这是由于商业区和办公区的充电路程用时要求较高,电动汽车容量密度较大,其空间权重较大,区域内有向加权距离相对更小,充电站的分区速度更快,因此总体上与充电站的距离更近。相反,工业区和休闲区的空间权重较小,因此总体上与充电站的距离更远。充电站的位置更靠近商业区等充电路程用时要求更高和电动汽车容量密度更大即对充电站需求更大的区域,证明了发明提出的元-点加权Voronoi图考虑空间权重的合理性以及元-点加权Voronoi图相对于加权Voronoi图在范围划分方面的优越性。
考虑用户充电体验指数前后,最优规划方案的充电站数分别为4座和6座,用户充电体验指标如表10所示。
表10充电体验指标
由表10可知,考虑充电体验指数后,充电路程用时满足率、充前SOC值满足率和无需排队满足率都得到了明显的提高,充电体验指数提高了23.6%,用户充电体验得到了更大的满足。因此以充电体验指数为优化目标,能够更好地满足用户的充电体验要求。
以充电站数为6时的最优规划方案为基础,进行充电站负荷接入节点优化,优化结果如表11和表12所示,优化前后越限最严重的几个元件在20:00时的越限概率及其状态数分别如表13和表14所示。。
表11优化前后接入节点
表12优化前后目标值
表13就近接入时风险最严重元件
表14优化后风险最严重元件
由表11和表12可知,优化后新建馈线成本增大为原来的3.4倍,配电网的运行风险减小为原来的4.83%。由表13和14可知,与充电站负荷就近接入时相比,优化后各元件的越限概率极大地减小了,电压和潮流状态数减小即元件的状态分布更加稳定了。由此可见,节点优化效果非常明显,能够以较小的新建馈线成本为代价,有效地减小配电网的运行风险,使得大规模充电负荷接入配电网时能够保证配电网的安全运行,证明了本发明提出的进行充电负荷接入节点优化的下层规划模型的有效性。
上述实施例为本发明较佳的实施方式,但本发明的实施方式并不受上述实施例的限制,其他任何未背离本发明的精神实质和原理下所作的修改、修饰、替代、组合、简化,均应为等效的置换方式,都应包含在本发明的保护范围之内。
Claims (7)
1.考虑用户充电体验和配电网运行风险的充电站规划方法,其特征在于,包括步骤:
步骤1、根据电动汽车历史数据,建立电动汽车负荷模型;
步骤2、建立充电体验指数,以充电体验指数、充电站建设成本和用户便捷性为充电站选址定容优化目标建立上层综合目标,并提出上层规划约束条件;
步骤3、提出计及空间权重的元-点加权Voronoi图进行服务范围分区;
步骤4、基于上层优化目标采用遗传禁忌混合算法建立上层规划流程;
步骤5、以新建馈线成本和配电网的运行风险为充电站负荷接入节点优化目标建立下层综合目标,并提出下层约束条件;
步骤6、基于下层综合目标采用循环替换法建立下层规划流程。
2.根据权利要求1所述的考虑用户充电体验和配电网运行风险的充电站规划方法,其特征在于,步骤1所述的根据电动汽车历史数据,建立电动汽车负荷模型包括以下步骤:
步骤1-1、采用电动汽车容量来描述区域内电动汽车的规模,所述电动汽车容量为区域内所有电动汽车额定充电功率的总和,包括处于充电状态和非充电状态的全部电动汽车;将单位面积电动汽车容量定义为电动汽车容量密度;
步骤1-2、根据电动汽车电池实际充电过程,构建电动汽车恒流-恒压二阶段变功率充电模型:
电动汽车在充电过程中,采用一种恒流-恒压二阶段变功率充电模式,当电池端电压小于最大电池端电压Ubmax时,电池是以恒电流模式进行充电;当电池端电压达到最大电池端电压Ubmax时,电池进入恒压阶段,电流迅速衰减为零,充电功率为零;描述其充电过程的方程为:
式中,Ub(t)为t时刻电池两端的电压,ic为充电时的电流,C0是电池的标称容量,R为电池内阻,K0是电池的极化常数,E0是电池的恒定电势,A0和B0分别为电池充电指数段的乘积系数和指数系数,SOC(t0)代表电池在开始充电时的剩余电量,h(t)为t时刻电池的可充电容量,即已放电容量;
在恒流充电过程中,电流i为常数,电动汽车的充电功率表示为:
P(t)=Ub(t)In (2)
式中,In为恒流充电模式下的充电电流,P(t)为t时刻的充电功率;
步骤1-3、根据电动汽车历史出行数据,构建电动汽车充电负荷时序概率模型:
电动汽车的充电负荷模拟须根据电动汽车的类型来选择各抽样环节的概率模型,抽取其行驶距离、到达时间、起始荷电状态和离开时间,从而确定每个时刻电动汽车的行为状态,进而计算相应时刻电动汽车的SOC值,得到其一天的SOC曲线,计算公式为:
(1)充电时
式中,SOC(t)是电动汽车的起始荷电状态;
(2)行驶时
SOC(t+Δt)=SOC(t)-dr/Dr·Δt/Tr (4)
式中,dr是电动汽车的行驶距离,Dr是电动汽车的最大续航里程,Tr是电动汽车的行驶时长,Δt为时间间隔;
然后,借助公式(1)和公式(2)算出其一天的充电功率;再将各类电动汽车的负荷逐辆叠加,得到电动汽车充电总负荷;t时刻电动汽车充电总负荷有功功率的计算公式为:
式中:N1、N2和N3分别代表t时刻该地区处于充电状态的电动私家车、公交车和出租车的数量;P1nt、P2nt和P3nt分别代表t时刻第n辆处于充电状态的电动私家车、公交车和出租车的实时充电功率;
步骤1-4、建立最大充电同时率:
为了清楚地了解某个区域电动汽车一天中同时充电的最大比例,定义电动汽车最大充电同时率为一天中电动汽车最大负荷与电动汽车容量的比值,即:
式中,Pt为t时刻区域内电动汽车充电总负荷,Kz为该区域的电动汽车容量。
3.根据权利要求1所述的考虑用户充电体验和配电网运行风险的充电站规划方法,其特征在于,步骤2建立充电体验指数,以充电体验指数、充电站建设成本和用户便捷性为充电站选址定容优化目标建立上层综合目标,并提出上层规划约束条件包括以下步骤:
步骤2-1、综合充电路程用时、充前SOC值和是否需要排队要求满足程度建立充电体验指数:
定义充电体验指数为充电过程中用户的充电体验要求得到满足的程度,并综合充电路程用时、充前SOC值和是否需要排队要求的满足程度来反映用户的充电体验指数;
(1)充电路程用时满足率
不同区域的人们由于时间的紧迫程度不同,对充电路上花费时间的要求不同,即充电路程用时要求不同;定义充电路程用时限值为用户对充电路上花费时间的接受极限值;充电过程中花费的时间大于充电路程用时限值时,则不满足充电体验要求;基于此,为考察所有用户的充电路程用时要求的满足程度,定义充电路程用时满足率为:
式中,ηt_sat为充电路程用时满足率;εd(i)为第i个点的电动汽车容量密度;va为充电过程中电动汽车的平均行驶速度;Td(i)为第i个点的充电路程用时限值,即其所属区域的充电路程用时限值;Mt_sat为满足充电路程用时要求的点的电动汽车容量密度的集合;N0为点的个数;Ld(i)为规划区第i个点到所属充电站的直角距离;
(2)充前SOC值满足率
电动汽车在前往充电前的SOC值满足一定概率分布,当SOC值低于一定值时,将会损害电动汽车电池的寿命,因此可定义充前SOC值满足率为用户到达充电站时SOC值满足不损害电动汽车电池寿命的最低SOC值的概率。通过M次蒙特卡罗抽样得到M个SOC值,则充前SOC值满足率的算法为:
式中,ηs_sat为充前SOC值满足率;pd(i)为第i个点满足充前SOC值要求的概率;εd(i)为第i个点的电动汽车容量密度;Srange为相应电动汽车的最大续航里程;SOC(j)为蒙特卡罗抽样的第j个SOC值;Ωd(i)为第i个点M个SOC值中满足充前SOC值要求的SOC值的集合;card()为集合中元素的个数,SOCmin为不损害电池寿命的最低SOC值;M为蒙特卡罗抽样次数;
(3)无需排队满足率
若电动汽车充电需求大于充电站的配置容量,则用户充电前需要排队,不满足无需排队要求。定义无需排队满足率为用户不需要排队的概率。通过蒙特卡罗抽样得到各充电站M次一天24个时刻的充电站负荷数据,则无需排队满足率的算法为:
式中,ηp_sat为无需排队满足率;pCD(i)为充电站i满足无需排队要求的平均概率;П(j)为充电站i第j次蒙特卡罗抽样一天中充电站负荷满足无需排队要求的时刻的集合;J(i)为充电站i的充电机配置数量;PCD(j,k)为充电站i第j次蒙特卡罗抽样第k个时刻的负荷;Pc为单台充电机的额定功率;αc为充电效率;KCD(i)表示充电站i的电动汽车容量,为其服务范围内所有点的电动汽车容量之和;
(4)用户充电体验指数计算方法
综合充电路程用时满足率、充前SOC值满足率和无需排队满足率来反映充电体验指数,算法为:
ηg_sat=ηt_satηs_satηp_sat (10)
式中,ηg_sat为充电体验指数;
步骤2-2、建立充电站建设成本:
充电站建设成本包括变压器、充电机和其他设备的建设成本以及土地费用;因此充电站i的建设成本的算法为:
Cz_t(i)=Cs(i)+Cd(i) (11)
式中,Cs(i)为第i个充电站的变压器、充电机和设备的建设成本,由充电站的充电机数量决定;Cd(i)为第i个充电站的土地费用,其算法为:
Cd(i)=Sd(i)Dd (12)
式中,Sd(i)为第i个充电站的占用土地面积,由充电站的充电机数量决定,Dd为单位面积土地费用;
则规划区的充电站建设总成本为:
Cg_t=∑Cz_t(i) (13)
步骤2-3、建立用户便捷性指标:
用户便捷性以充电平均距离来表示;在规划区内平均取点,每个点的电动汽车容量密度为其所属区域的电动汽车容量密度,以每个点的电动汽车容量密度为权重对充电距离进行加权平均得到充电平均距离;则规划区的充电平均距离的算法为:
式中,εd(i)和Ld(i)分别为规划区第i个点的电动汽车容量密度和其到所属充电站的直角距离;
步骤2-4、建立上层综合目标:
充电站的上层规划以充电站建设成本、用户便捷性和充电体验指数为目标,对充电站进行选址定容优化;因此,上层综合目标的算法为:
式中,α1、α2和α3为权重系数,其作用为各指标的数量级相同,同时用于消掉各指标的单位,使得上层综合目标为无量纲,B1为上层综合目标;
步骤2-5、建立上层规划约束条件:
(1)相邻站间距离约束
为保证用户的充电需求,充电站的服务范围不能太大,同时为了避免资源浪费,应避免充电站分布过于集中;因此相邻充电站间的直角距离Lcc要满足约束条件:
Lmin<Lcc<Lmax (16)
式中,Lmin和Lmax分别为相邻站间距离约束的上下限;
(2)充电机数约束
为了满足用户的充电需求,同时出于对社会资源的充分利用,保证充电机有较高的利用率,充电站的充电机配置数量应满足充电机数约束:
Jmin≤JCD≤Jmax (17)
式中,Jmin和Jmax分别为最小充电机配置数和最大充电机配置数,JCD为充电站充电机的配置数量,算法为:
式中,ZCD是充电站服务范围内的最大充电需求;KCD为充电站服务范围内的电动汽车容量。
4.根据权利要求1所述的考虑用户充电体验和配电网运行风险的充电站规划方法,其特征在于,步骤3提出计及空间权重的元-点加权Voronoi图进行服务范围分区包括以下步骤:
步骤3-1、提出元-点加权Voronoi图:
针对传统的加权Voronoi图不能反映不同区域权重差异的缺陷,提出元-点加权Voronoi图来进行服务范围划分,以空间权重描述不同区域权重的大小;生成元生成Voronoi图的速度不仅与生成元的权重有关,还与各区域的空间权重有关;由于考虑空间权重,不同的行驶路径两点之间的加权距离不同,因此定义点A到点B的有向路径为A向B先纵向行驶,再横向行驶的路径,并定义加权有向距离为有向路径的加权距离;
基于此,元-点加权Voronoi图定义为:
V(Ri)={P∈V(Ri)|D(Ri,P)≤D(Rj,P),j=1,2,…,m1;j≠i} (19)
式中,D(Ri,P)为点Ri到点P的加权有向距离,m1为生成元的个数,D(Rj,P)为点Rj到点P的加权有向距离,若Ri的权重为vi,Ri到点P的有向路径经过区域Q1,Q2,…,Qn,n为Ri到点P的有向路径经过的区域个数,相应的空间权重为θ1,θ2,…,θn,Ri到点P的有向路径经过的边界点分别为P1,P2,…,Pn-1,d(Pn-1,P)为点Pn-1到点P的有向距离,则Ri到点P的加权有向距离为:
步骤3-2、建立充电站权重:
在地价低的地方建设服务范围更大、容量更大的充电站,能减小充电站建设总成本;因此在服务范围划分时,设充电站权重的平方与其占用土地单价的倒数成反比,则充电站i的权重为:
式中,Di为充电站i的占用土地单价;
步骤3-3、建立空间权重:
综合不同类型区域的充电路程用时要求和电动汽车容量密度来建立区域的空间权重,则区域类型i的空间权重的算法为:
式中,θq(i)、εq(i)、Tq(i)分别为区域类型i的空间权重、电动汽车容量密度和充电路程用时限值,β1和β2为权重系数;
步骤3-4、充电站服务范围划分:
为便于计算,在利用元-点加权Voronoi图进行充电站规划时,规划区内均匀生成点,邻近点的移动步长为Δx=Δy=Δl,用矩阵分别记录每个点的坐标、电动汽车容量密度、土地单价、充电路程用时限值信息,建立点数据库;在画Voronoi图时,每个点都用一个方块来代表,方块的边长为移动步长,点位于方块的中心;方块的电动汽车容量为方块的面积与点的电动汽车容量密度的乘积;
若权重为vi的充电站Ri到点O的有向路径依次经过O1,O2,…,On-1,各点相应的空间权重为θq1,θq2,…,θqn,则Ri到点O的加权有向距离为:
将各生成元到某点的加权有向距离进行比较,加权有向距离最小对应的生成元将该点划进其服务范围;
充电站的电动汽车容量为其服务范围内所有点的电动汽车容量之和,算法为:
式中,εd(i)为该充电站服务范围内第i个点的电动汽车容量密度,ΔS为单个方块的面积。
5.根据权利要求1所述的考虑用户充电体验和配电网运行风险的充电站规划方法,其特征在于,步骤4基于上层优化目标采用遗传禁忌混合算法建立上层规划流程包括:
采用遗传禁忌混合算法(GATS)进行充电站上层规划,主要步骤如下:
Ⅰ.采集点数据库信息,作为优化的基础数据;
Ⅱ.利用点数据库的电动汽车容量密度信息,求出规划区内的电动汽车容量,即Kz,进而计算规划区的最大充电需求Zz,算法为:
Ⅲ.根据规划区内的最大充电需求Zz,以充电站的最小充电机配置数Jmin和最大充电机配置数Jmax为约束,计算充电站的最小座数Nmin和最大座数Nmax,即:
Ⅳ.基于上层综合目标采用GATS和元-点加权Voronoi图进行选址定容优化,包括以下步骤:
(1)记充电站数为g,设置GA的最大迭代数Ngen、禁忌表长度Tlength1、群体规模Npop、重组概率pch和变异概率pm;TSM中禁忌表长度Tlength2、邻域解个数Nh和最大迭代数Nt参数;
(2)tr表示迭代次数,令tr=0,将规划区平均分为g个区,分别从各区中随机抽取一个地址,得到染色体的地址序列[x1,y1;…;xi,yi;…;xg,yg],(xi,yi)为染色体的地址序列中的第i个地址,地址之间不满足站间距离约束的重新生成;
(3)计算当代群体的上层综合目标;
(4)按滚轮盘方式选出Npop个染色体;
(5)生成0和1之间的随机数ri,i=1,2,…,Npop,若ri<pch,则第i个染色体为父代;对每对染色体进行交叉操作,随机互换两个基因,判断子代是否满足站间距离约束,不满足则选最优父代进下一代;若满足,则将子代的上层综合目标同渴望水平相比,若更好,则破禁进下一代;若更差,但不属于禁忌,则接受这个子代,若属于禁忌,则选最优父代进下一代;
(6)生成0和1之间的随机数qi,i=1,2,…,Npop;若qi<pm,则对第i个染色体进行变异操作:
a.以当前染色体作为最优染色体,将其上层综合目标作为最优上层综合目标,置TSM的禁忌表为空;
b.以当前地址基因为中心生成邻域地址基因(x+Δx,y)、(x-Δx,y)、(x,y+Δy)和(x,y-Δy),加上当前地址基因构成等位基因序列,进行全排列得到5^g个染色体并从中随机抽出Nh个;
c.将这Nh个染色体中不满足站间距离约束的剔除;剩余的不在禁忌表中的染色体为候选染色体,计算其分区后满足充电机数约束的上层综合目标;
d.以上层综合目标最小选出最优候选染色体,并用它替代当前染色体和最早进入不满足的染色体;
e.判断最优候选染色体的上层综合目标是否小于最优上层综合目标,是则用它替换最优染色体,并将其上层综合目标替换最优上层综合目标,然后转步骤f;否则直接转到步骤f;
f.判断是否达到最大迭代数Nt,或最优上层综合目标已连续30次保持不变,是则输出最优染色体;否则返回步骤b;
(7)tr=tr+1,若tr<Ngen,转步骤(3);否则输出最优解。
6.根据权利要求1所述的考虑用户充电体验和配电网运行风险的充电站规划方法,其特征在于,步骤5以新建馈线成本和配电网的运行风险为充电站负荷接入节点优化目标建立下层综合目标,并提出下层约束条件包括以下步骤:
步骤5-1、计算新建馈线成本:
充电站负荷要接入配电网,要通过建设新的馈线来接入配电网的节点;假定新建馈线是沿直线铺设的,则充电站的新建馈线成本算法为:
Ck=Dk∑Lk(i) (27)
式中,Dk为馈线的单位长度造价;Lk(i)为第i个充电站的新建馈线长度,即第i个充电站到接入节点的直角距离;
步骤5-2、计算配电网运行风险:
(1)元件风险指标
基于加权熵原理,以加权电压分布熵表征电压越限及电压分布不确定性引起的电压风险,以加权潮流分布熵表征潮流越限及潮流分布不确定性引起的潮流风险,算法为:
式中,Hv_i(t)和Hl_i(t)分别为t时刻节点i的加权电压分布熵和支路i的加权潮流分布熵;nv_i(t)和nl_i(t)分别为t时刻节点i的电压状态数和支路i的潮流状态数;Sv_j(t)和Pv_j(t)分别为t时刻节点i第j个电压状态的电压损失严重度和概率;Sl_j(t)和Pl_j(t)分别为t时刻支路i第j个潮流状态的潮流损失严重度和概率;
损失严重度S的算法为:
式中,μ为损失严重度的放大系数,以调节损失严重度的灵敏性;w为损失量;
(2)网络风险指标
综合电压风险和潮流风险对含大规模电动汽车的配电网进行运行风险评估;为考察配电网运行风险一天的变化情况,将t时刻网络运行风险定义为:
Hw(t)=γ1Hv(t)+γ2Hl(t) (30)
式中,Hw(t)为t时刻的网络运行风险;γ1和γ2为权重系数;Hv(t)和Hl(t)分别为t时刻的网络电压风险和网络潮流风险,算法为:
式中,z1和z2分别为配电网的节点数和支路数;
为研究配电网安全运行面对的最严重的风险情况,定义网络最大运行风险为网络运行风险一天的最大值,即:
Hwmax=max(Hw(t)) t=1,2,…,24 (32)
步骤5-3、提出下层综合目标:
充电站下层规划综合考虑新建馈线成本和配电网的运行风险,对充电站负荷接入节点进行优化;定义下层综合目标为:
B2=ρ1Ck+ρ2Hwmax (33)
式中,Ck和Hwmax分别为新建馈线成本和网络最大运行风险;ρ1和ρ2为权重系数;
步骤5-4、提出下层规划约束条件:
(1)馈线长度约束
当新建馈线太长时,不仅馈线成本过大,同时馈线的阻抗过大,会增大功率损耗和电压损耗,因此新建的馈线应满足馈线长度约束:
Lk(i)≤Lkmax (34)
式中,Lkmax为新建馈线的长度限值;
(2)充电站负荷约束
由于充电站的充电负荷不可能大于其安装容量,因此要满足约束:
PCD(i)≤αcJ(i)Pc (35)
式中,PCD(i)为充电站i的充电负荷。
7.根据权利要求1所述的考虑用户充电体验和配电网运行风险的充电站规划方法,其特征在于,步骤6基于下层综合目标采用循环替换法建立下层规划流程包括以下步骤:
(1)输入充电站地址及各分区电动汽车容量、充电过程参数、配电网络参数数据;
(2)通过蒙特卡罗抽样计算各个充电站分区的Nm次日充电负荷数据;
(3)将充电站负荷就近接入配电网节点作为当前方案,计算下层综合目标;
(4)记充电站的个数为g,确定各分区的优化顺序,bm=0;
(5)bm=bm+1,计算充电站bm到各节点的距离,将满足约束的节点作为候选节点;
(6)保持其他接入节点不变,将充电站bm的负荷分别接入各候选节点,计算下层综合目标,以下层综合目标最小选最优候选方案;
(7)判断最优候选方案与当前方案的优劣,若最优候选方案比当前方案更优,则以最优候选方案替换当前方案,否则保留当前方案不变;
(8)判断是否有bm<g,是则返回步骤(5),否则转到步骤(9);
(9)不断循环优化,直至充电站负荷接入方案不再改变为止。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910820177.9A CN110826841A (zh) | 2019-08-31 | 2019-08-31 | 考虑用户充电体验和配电网运行风险的充电站规划方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910820177.9A CN110826841A (zh) | 2019-08-31 | 2019-08-31 | 考虑用户充电体验和配电网运行风险的充电站规划方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110826841A true CN110826841A (zh) | 2020-02-21 |
Family
ID=69548095
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910820177.9A Pending CN110826841A (zh) | 2019-08-31 | 2019-08-31 | 考虑用户充电体验和配电网运行风险的充电站规划方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110826841A (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111609867A (zh) * | 2020-06-19 | 2020-09-01 | 北京交通大学 | 一种电动汽车路径规划方法 |
CN111628496A (zh) * | 2020-05-19 | 2020-09-04 | 南京工程学院 | 一种电动汽车充电站选址和最大充电负荷确定方法 |
CN111641595A (zh) * | 2020-05-11 | 2020-09-08 | 湖南大学 | 一种电力网络安全风险评估方法及系统 |
CN111882123A (zh) * | 2020-07-16 | 2020-11-03 | 成都市高博汇科信息科技有限公司 | 一种大数据智能决策的智慧社区物流终端选址方法 |
CN112468984A (zh) * | 2020-11-04 | 2021-03-09 | 国网上海市电力公司 | 一种电力无线专网基站选址方法及相关设备 |
CN112757954A (zh) * | 2020-12-30 | 2021-05-07 | 华南理工大学 | 联合专变共享模式下的电动汽车有序充电联合调节方法 |
CN115130810A (zh) * | 2022-04-07 | 2022-09-30 | 兰州理工大学 | 电动汽车充电站扩容规划方法、系统、介质、设备及终端 |
CN115848196A (zh) * | 2022-12-07 | 2023-03-28 | 南通国轩新能源科技有限公司 | 基于动态需求和新能源消纳的电动汽车有序充电引导方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102880921A (zh) * | 2012-10-16 | 2013-01-16 | 山东电力集团公司电力科学研究院 | 一种电动汽车充电站选址优化方法 |
US20160137076A1 (en) * | 2014-11-14 | 2016-05-19 | General Electric Company | System, method and apparatus for remote opportunity charging |
CN109034465A (zh) * | 2018-07-16 | 2018-12-18 | 河海大学 | 考虑充电站选址与出行路径耦合的充电站两层规划方法 |
CN109117486A (zh) * | 2017-06-23 | 2019-01-01 | 南京理工大学 | 一种电动汽车充电站最优规划方法 |
CN109840635A (zh) * | 2019-01-29 | 2019-06-04 | 三峡大学 | 基于电压稳定性和充电服务质量的电动汽车充电站规划方法 |
CN110033122A (zh) * | 2019-03-08 | 2019-07-19 | 中国电力科学研究院有限公司 | 一种电动汽车充电站选址方法及装置 |
-
2019
- 2019-08-31 CN CN201910820177.9A patent/CN110826841A/zh active Pending
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102880921A (zh) * | 2012-10-16 | 2013-01-16 | 山东电力集团公司电力科学研究院 | 一种电动汽车充电站选址优化方法 |
US20160137076A1 (en) * | 2014-11-14 | 2016-05-19 | General Electric Company | System, method and apparatus for remote opportunity charging |
CN109117486A (zh) * | 2017-06-23 | 2019-01-01 | 南京理工大学 | 一种电动汽车充电站最优规划方法 |
CN109034465A (zh) * | 2018-07-16 | 2018-12-18 | 河海大学 | 考虑充电站选址与出行路径耦合的充电站两层规划方法 |
CN109840635A (zh) * | 2019-01-29 | 2019-06-04 | 三峡大学 | 基于电压稳定性和充电服务质量的电动汽车充电站规划方法 |
CN110033122A (zh) * | 2019-03-08 | 2019-07-19 | 中国电力科学研究院有限公司 | 一种电动汽车充电站选址方法及装置 |
Non-Patent Citations (1)
Title |
---|
方柳佳: ""电动公交车充换电站布局规划研究"", 《中国优秀硕士学位论文全文数据库 工程科技II辑》, no. 12, pages 034 - 196 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111641595A (zh) * | 2020-05-11 | 2020-09-08 | 湖南大学 | 一种电力网络安全风险评估方法及系统 |
CN111628496A (zh) * | 2020-05-19 | 2020-09-04 | 南京工程学院 | 一种电动汽车充电站选址和最大充电负荷确定方法 |
CN111609867A (zh) * | 2020-06-19 | 2020-09-01 | 北京交通大学 | 一种电动汽车路径规划方法 |
CN111882123A (zh) * | 2020-07-16 | 2020-11-03 | 成都市高博汇科信息科技有限公司 | 一种大数据智能决策的智慧社区物流终端选址方法 |
CN111882123B (zh) * | 2020-07-16 | 2023-06-23 | 成都市高博汇科信息科技有限公司 | 一种大数据智能决策的智慧社区物流终端选址方法 |
CN112468984A (zh) * | 2020-11-04 | 2021-03-09 | 国网上海市电力公司 | 一种电力无线专网基站选址方法及相关设备 |
CN112468984B (zh) * | 2020-11-04 | 2023-02-10 | 国网上海市电力公司 | 一种电力无线专网基站选址方法及相关设备 |
CN112757954A (zh) * | 2020-12-30 | 2021-05-07 | 华南理工大学 | 联合专变共享模式下的电动汽车有序充电联合调节方法 |
CN115130810A (zh) * | 2022-04-07 | 2022-09-30 | 兰州理工大学 | 电动汽车充电站扩容规划方法、系统、介质、设备及终端 |
CN115848196A (zh) * | 2022-12-07 | 2023-03-28 | 南通国轩新能源科技有限公司 | 基于动态需求和新能源消纳的电动汽车有序充电引导方法 |
CN115848196B (zh) * | 2022-12-07 | 2024-01-05 | 南通国轩新能源科技有限公司 | 基于动态需求和新能源消纳的电动汽车有序充电引导方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110826841A (zh) | 考虑用户充电体验和配电网运行风险的充电站规划方法 | |
Li et al. | Public charging station location determination for electric ride-hailing vehicles based on an improved genetic algorithm | |
Asna et al. | Analysis of an optimal planning model for electric vehicle fast-charging stations in Al Ain City, United Arab Emirates | |
Pal et al. | Placement of public fast-charging station and solar distributed generation with battery energy storage in distribution network considering uncertainties and traffic congestion | |
Bayram et al. | Electric power allocation in a network of fast charging stations | |
Xiang et al. | Electric vehicles in smart grid: a survey on charging load modelling | |
Pan et al. | Deploying public charging stations for electric taxis: A charging demand simulation embedded approach | |
Meng et al. | Sequential construction planning of electric taxi charging stations considering the development of charging demand | |
CN110929925A (zh) | 基于环境温度与用户行为的电动汽车充电站规划方法 | |
CN106875075A (zh) | 一种基于出行行为学的电动汽车充电站布点方法 | |
CN113222241B (zh) | 考虑充电服务引导与客户需求的出租车快充电站规划方法 | |
CN113553769B (zh) | 一种考虑充电效能的电动汽车充电设施优化配置方法 | |
Liu et al. | Solving urban electric transit network problem by integrating Pareto artificial fish swarm algorithm and genetic algorithm | |
CN115907227B (zh) | 一种高速公路固定与移动充电设施双层协同优化方法 | |
Wang et al. | Mobile electric vehicles | |
Bayram et al. | Could petrol stations play a key role in transportation electrification? a gis-based coverage maximization of fast ev chargers in urban environment | |
Azin et al. | Infrastructure enabled and electrified automation: Charging facility planning for cleaner smart mobility | |
Zhang et al. | An immune algorithm for solving the optimization problem of locating the battery swapping stations | |
Krallmann et al. | Multi-objective optimization of charging infrastructure to improve suitability of commercial drivers for electric vehicles using real travel data | |
He et al. | Expansion planning of electric vehicle charging stations considering the benefits of peak‐regulation frequency modulation | |
Wang et al. | A real-time strategy for vehicle-to-station recommendation in battery swapping mode | |
Zhang et al. | Multi-period planning of locations and capacities of public charging stations | |
da Paixão et al. | EV fast charging microgrid on highways: A hierarchical analysis for choosing the installation site | |
Sadhukhan et al. | Optimal placement of electric vehicle charging stations in a distribution network | |
Sun et al. | Location and Capacity Determination Method of Electric Vehicle Charging Station Based on Simulated Annealing Immune Particle Swarm Optimization |
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 |