CN106340873A - 一种采用基于无向生成树的并行遗传算法的配网重构方法 - Google Patents
一种采用基于无向生成树的并行遗传算法的配网重构方法 Download PDFInfo
- Publication number
- CN106340873A CN106340873A CN201610887323.6A CN201610887323A CN106340873A CN 106340873 A CN106340873 A CN 106340873A CN 201610887323 A CN201610887323 A CN 201610887323A CN 106340873 A CN106340873 A CN 106340873A
- Authority
- CN
- China
- Prior art keywords
- node
- undirected
- population
- spanning tree
- distribution
- 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
- 238000004422 calculation algorithm Methods 0.000 title claims abstract description 64
- 238000000034 method Methods 0.000 title claims abstract description 36
- 230000002068 genetic effect Effects 0.000 title claims abstract description 34
- 238000005070 sampling Methods 0.000 claims abstract description 16
- 238000004364 calculation method Methods 0.000 claims abstract description 11
- 238000000342 Monte Carlo simulation Methods 0.000 claims abstract description 4
- 238000013507 mapping Methods 0.000 claims abstract description 3
- 230000035772 mutation Effects 0.000 claims description 12
- 210000000349 chromosome Anatomy 0.000 claims description 9
- 230000008859 change Effects 0.000 claims description 4
- 230000005855 radiation Effects 0.000 claims description 4
- 238000002347 injection Methods 0.000 claims description 3
- 239000007924 injection Substances 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 230000008901 benefit Effects 0.000 abstract description 3
- 230000010354 integration Effects 0.000 abstract 1
- 238000012360 testing method Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 238000005457 optimization Methods 0.000 description 5
- 238000010845 search algorithm Methods 0.000 description 4
- 239000000243 solution Substances 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 239000002245 particle Substances 0.000 description 3
- 238000011156 evaluation Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 241000196324 Embryophyta Species 0.000 description 1
- 235000006508 Nelumbo nucifera Nutrition 0.000 description 1
- 240000002853 Nelumbo nucifera Species 0.000 description 1
- 235000006510 Nelumbo pentapetala Nutrition 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 238000010353 genetic engineering Methods 0.000 description 1
- 244000144980 herd Species 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/18—Network design, e.g. design based on topological or interconnect aspects of utility systems, piping, heating ventilation air conditioning [HVAC] or cabling
-
- 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
-
- H02J3/383—
-
- H02J3/386—
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
-
- 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
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/70—Wind energy
- Y02E10/76—Power conversion electric or electronic aspects
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Biophysics (AREA)
- Health & Medical Sciences (AREA)
- Geometry (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Biology (AREA)
- Pure & Applied Mathematics (AREA)
- Artificial Intelligence (AREA)
- Computer Hardware Design (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physiology (AREA)
- Genetics & Genomics (AREA)
- Power Engineering (AREA)
- Biomedical Technology (AREA)
- Computational Linguistics (AREA)
- Data Mining & Analysis (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明涉及一种采用基于无向生成树的并行遗传算法的配网重构方法,包括步骤:获取参数;进行蒙特卡洛模拟采样;随机产生拓扑可行的初始种群,设迭代次数n的初值为1;进行潮流计算;计算目标函数值,判断是否满足约束条件,若否则返回步骤重新产生初始种群,若是则对现有种群分为多个子种群进行并行遗传操作;产生一个1~Nsub的随机排列P=[p1,p2,…,pNsub],建立目标子种群i与源子种群pi间的映射关系;将各个目标子种群的最差个体替换为相应的源子种群的最优个体;判断迭代次数n是否达到要求,若否,则迭代次数n增加1,返回进行潮流计算,若是,则输出配电网重构方案。与现有技术相比,本发明具有计算效率高、整合性强、紧密结合实际等优点。
Description
技术领域
本发明涉及配网重构方法,尤其是涉及一种采用基于无向生成树的并行遗传算法的配网重构方法。
背景技术
配电网重构是配电网优化运行的一种手段,在满足闭环结构、开环运行以及电能质量要求的前提下,通过改变联络开关或分段开关的状态,实现配电网安全、经济运行。
分布式电源(distributed generation,DG)因配置灵活,并且具有较好的环境效益,在国内外得到了广泛应用。随着配电网中分布式电源的渗透率的提高,配电网正由传统的简单辐射无源网络向多端有源网络发展,进入了主动配电网发展的新阶段。分布式电源作为一种重要元素,也应该在配电网重构建模中予以考虑。
经对现有文献进行检索发现,现有文献中,王少林、唐巍、白牧可等在《电力系统保护与控制》(2012,40(18):117-122)上发表的《考虑分布式电源出力调整的多目标配电网重构》以开关状态、DG出力为优化变量,建立了以网络损耗、负荷均衡化率最小为目标函数的多目标优化模型,将生成树、蚁群算法和遗传算法相结合,提出了求解上述模型的多目标混合优化方法,以实现配电网结构和DG出力的协同优化;丛鹏伟、唐巍、张璐等在《电网技术》(2013,37(9):2573-2579)上发表的《基于机会约束规划考虑DG与负荷多状态的配电网重构》考虑自然资源和负荷的随机性,依据概率密度函数分别对不可控DG和负荷建立了多状态模型,进而获得了配电系统的多状态模型,以开关状态为优化变量,高于置信水平的网络损耗最小为目标函数,节点电压和支路功率越限概率满足置信水平要求为约束条件,建立了含多种DG的配电网重构机会约束规划模型;向小蓉、刘涤尘、向农等在《电网技术》(2012,36(8):100-105)上发表的《基于并行禁忌搜索算法的配电网重构》将配电网简化为只含支路组的环网,并给出了配电网重构的必要条件,重构时首先选择断开的支路组,采用深度优先搜索算法形成父子链表,然后引入可降低复杂度的分治思想,以可断开支路组为解空间进行分组,应用多处理器进行并行禁忌搜索。以上文献对DG接入后的配电网重构进行了研究,但对于DG出力及负荷的不确定性考虑不多,算法在复杂电网中的计算效率问题值得深入研究。
发明内容
本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种计算效率高、整合性强、紧密结合实际的一种采用基于无向生成树的并行遗传算法的配网重构方法。
本发明的目的可以通过以下技术方案来实现:
一种采用基于无向生成树的并行遗传算法的配网重构方法,其特征在于,该方法包括下列步骤:
(1)获取配网参数、风机和光伏参数以及并行遗传算法参数;
(2)根据步骤(1)中获取的参数对风速、辐射强度和负荷进行基于拉丁超立方的蒙特卡洛模拟采样;
(3)随机产生包含Nind个拓扑可行的染色体的初始种群作为目标种群,并设置迭代次数n的初值为1;
(4)根据步骤(2)中得到的采样值,对目标种群中的每个染色体所对应的网络拓扑进行潮流计算,得到系统网损期望、节点电压及支路潮流越限概率;
(5)根据步骤(4)中潮流计算得到的结果计算目标函数值,判断是否满足约束条件,若是则进入步骤(6),若否则返回步骤(3);
(6)对现有种群分为多个子种群,进行基于无向生成树的并行遗传操作,得到多个目标子种群;
(7)产生一个1~Nsub的随机排列P=[p1,p2,…,pNsub],建立目标子种群i与源子种群pi间的映射关系;
(8)将各个目标子种群的最差个体替换为相应的源子种群的最优个体,合并各个目标子种群作为目标种群;
(9)判断迭代次数n是否达到要求,若是,则进入步骤(10),若否,则迭代次数n增加1,返回步骤(4);
(10)根据得到的目标种群输出配电网重构方案。
所述的配网参数包括网络拓扑结构、各节点负荷、网络支路阻抗及导纳、网络电压等级;所述的风机参数包括切入风速、额定风速、切出风速、额定功率、功率因数、Weibull分布参数;所述的输入模块中的光伏参数包括光伏面板面积、光电效率、额定功率、功率因数、Beta分布参数。
所述的Weibull分布用以描述风电场风速概率分布,概率密度函数可表示为:
其中,v为风速,c和k分别为Weilbull分布的尺度参数和形状参数;
风机出力与风速的关系可表示为:
其中,vci、vcr和vco分别为切入风速、额定风速和切出风速;Pr为风机额定输出有功功率。
所述的Beta分布用以描述光照辐射强度概率分布,其概率密度可表示为:
其中,α和β为Beta分布的形状参数;r和rmax分别为某一时段内的太阳辐照强度及其最大值;
光伏输出有功功率可表示为:
PM=rMAη
其中,PM为光伏阵列输出有功功率;M为光伏阵列中太阳能电池组件个数;A为每个太阳能电池组件的面积;η为光伏阵列的光电转换效率。
所述的基于拉丁超立方的蒙特卡洛模拟采样具体步骤为:
(21)对每个输入变量进行采样,确保随机分布区域能够被采样点完全覆盖;
(22)改变变量采样值的排列顺序,使相对独立的随机变量采样值的相关性趋于最小。
所述的目标函数为:
式中,PLoss为系统网损期望值,;E{·}表示计算随机变量的数学期望;NL为系统支路数;Pk,Loss为随机变量,对应于第k条支路的网损;NLHS为LHS采样次数;Ii,k为第i次采样中流过第k条支路的潮流,其中i=1,2,…,NLHS;rk为第k条支路电阻,其中k=1,2,…,NL。
所述的约束条件包括:
1)潮流方程约束
其中,Pis和Qis分别为节点i处注入有功功率和无功功率;Vi为节点i处电压幅值;Ωi表示与节点i相连的其他节点集合;Gij、Bij分别为支路ij的电导和电纳;θij为节点i与节点j间相角差,其中i,j=1,2,...,n,n为节点总数;
2)节点电压约束
Pr{Vi min≤Vi≤Vi max,i∈Θ}≥βV
其中,Pr{·}为事件成立概率;Vi min和Vi max分别为节点i处电压幅值的上下限;βV为电压约束的置信水平;Θ为系统节点集合;
3)支路潮流约束
其中,Pj max为第j条支路传输有功功率的上限;βI为支路潮流约束的置信水平;Λ为系统支路集合;
4)网络拓扑约束
g∈G
其中,g为当前网络拓扑结构;G为所有满足辐射状结构的网络拓扑集合。
所述的基于无向生成树的并行遗传操作包括选择操作、交叉操作、变异操作。
所述的选择操作采用轮盘赌方法;所述的交叉操作具体步骤为:
(61)随机选择节点N1与N2,使得N1≠N2;
(62)对于父代生成树T1与T2,通过深度优先搜索确定N1与N2之间的无向路径,其中T1中的路径为P1,T2中的路径为P2;
(63)分别将T1与P2、T2与P1进行交叉操作,得到新个体T1 *与T2 *;
所述的变异操作具体步骤为:
(64)获取父代拓扑T;
(65)随机选择节点N1与N2,使得N1≠N2;
(66)判断T中是否存在以N1和N2为端点的支路,若否,则进入步骤(67),若是,则返回步骤(65)。
(67)通过深度优先搜索找出N1和N2之间的无向路径P,随机选取P中的任意一条支路,将其删除,同时在T中增加支路N1N2,从而产生新个体T*。
与现有技术相比,本发明采用了基于无向生成树的并行遗传算法,具有以下有益效果:
1)有效解决了分布式电网接入后的配电网重构问题,降低电网系统的网络损耗,有效改善馈线末端节点电压分布,缓解节点电压越线,具有较强的实用性;
2)相较于STGA、PSO、IHSA、ACS等算法,具有不错的稳定性、较好的收敛速度以及高效的并行能力;
3)整个方法最终进行了模块化处理,便于程序运行,具有较强的整合性。
附图说明
图1为本发明的结构示意图;
图2为并行计算部分流程图;
图3为交叉操作流程图;
图4为交叉操作示意图;
图5为变异操作流程图;
图6为变异操作示意图;
图7为IEEE33节点系统接线图;
图8为不同情形下节点电压分布情况;
图9为重构前后节点18处电压概率分布情况(含DG);
图10为重构前后节点33处电压概率分布情况(含DG);
图11为PSTGA、STGA及PSO稳定性比较;
图12为PSTGA、STGA及PSO收敛速度比较。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。
实施例1
如图1所示,一种采用基于无向生成树的并行遗传算法的配网重构方法,由依次连接的输入模块、初始化模块、随机潮流计算模块、并行遗传操作模块和输出模块实现。
初始化模块从输入模块获取配网参数、风机和光伏参数、PSTGA参数,并随机产生拓扑可行的染色体初始种群;随机潮流计算模块对初始化模块产生的初始化种群中的每个染色体所对应的网络拓扑进行随机潮流计算,得到系统网损期望、节点电压和支路潮流越限概率,进而计算计及惩罚项的目标函数值;并行遗传操作模块对现有种群分为多个子种群进行并行遗传操作,产生下一代种群;输出模块输出配电网重构方案。本发明总体结构示意图如附图1所示。
输入模块中的配网参数包括网络拓扑结构、各节点负荷、网络支路阻抗及导纳、网络电压等级;输入模块中的风机参数包括切入风速、额定风速、切出风速、额定功率、功率因数、Weibull分布参数;输入模块中的光伏参数包括光伏面板面积、光电效率、额定功率、功率因数、Beta分布参数。
Weibull分布用以描述风电场风速概率分布,概率密度函数可表示为:
其中,v为风速,c和k分别为Weilbull分布的尺度参数和形状参数。
风机出力与风速的关系可表示为:
其中,vci、vcr和vco分别为切入风速、额定风速和切出风速,单位m/s;Pr为风机额定输出有功功率,单位kW。
Beta分布用以描述光照辐射强度概率分布,其概率密度可表示为:
其中,α和β为Beta分布的形状参数;r和rmax分别为某一时段内的太阳辐照强度及其最大值,单位kW/m2。
光伏输出有功功率可表示为:
PM=rMAη
其中,PM为光伏阵列输出有功功率,单位kW;M为光伏阵列中太阳能电池组件个数;A为每个太阳能电池组件的面积,单位m2;η为光伏阵列的光电转换效率。
初始化模块中的基于拉丁超立方的蒙特卡洛模拟方法用以对风速、光照强度及负荷进行随机抽样,以进行概率潮流计算。拉丁超立方采样属于分层采样,首先对每个输入变量进行采样,确保随机分布区域能够被采样点完全覆盖;随后改变变量采样值的排列顺序,使相对独立的随机变量采样值的相关性趋于最小。
随机潮流计算模块将各组采样值依次代入潮流方程,求出节点电压、支路潮流及系统网损,得到系统网损均值,节点电压和支路潮流越限概率等,进而计算目标函数值,并判断是否满足约束条件。
目标函数为:
其中,PLoss为系统网损期望值,单位MW;E{·}表示计算随机变量的数学期望;NL为系统支路数;Pk,Loss为随机变量,对应于第k条支路的网损,单位MW;NLHS为LHS采样次数;Ii,k为第i次采样中流过第k条支路的潮流,单位kA,其中i=1,2,…,NLHS;rk为第k条支路电阻,单位Ω,其中k=1,2,…,NL。
约束条件为:
1)潮流方程约束
其中,Pis和Qis分别为节点i处注入有功功率和无功功率,单位MW;Vi为节点i处电压幅值,单位kV;Ωi表示与节点i相连的其他节点集合;Gij、Bij分别为支路ij的电导和电纳,单位S;θij为节点i与节点j间相角差,其中i,j=1,2,...,n,n为节点总数。
2)节点电压约束
Pr{Vi min≤Vi≤Vi max,i∈Θ}≥βV
其中,Pr{·}为事件成立概率;Vi min和Vi max分别为节点i处电压幅值的上下限,单位kV;βV为电压约束的置信水平;Θ为系统节点集合。
3)支路潮流约束
其中,Pj max为第j条支路传输有功功率的上限,单位MW;βI为支路潮流约束的置信水平;Λ为系统支路集合。
4)网络拓扑约束
g∈G
其中,g为当前网络拓扑结构;G为所有满足辐射状结构的网络拓扑集合。
并行遗传操作模块中的遗传操作包括选择算子、交叉算子、变异算子;并行计算将所有染色体进行等分为若干子种群,并行地对各子种群进行选择、交叉和变异操作,接着产生一个1~Nsub的随机排列P=[p1,p2,…,pNsub],从而建立目标子种群i与源子种群pi间的映射关系;最后将各个目标子种群的最差个体替换为相应的源子种群的最优个体。并行计算部分的流程图如附图2所示。
选择操作采用轮盘赌方法,使本代适应值较高的个体有较高概率遗传到下一代;交叉算子首先随机选择节点N1与N2,使得N1≠N2。然后,对于父代生成树T1与T2,通过深度优先搜索(Depth-First Search,DFS)确定N1与N2之间的无向路径,不妨设T1中的路径为P1,T2中的路径为P2,然后分别将T1与P2、T2与P1进行交叉操作,得到新个体T1 *与T2 *,交叉操作的流程图和示意图可见附图3、4;变异算子对父代拓扑T,随机选择节点N1与N2,使得N1≠N2,如果T中不存在以N1和N2为端点的支路,那么通过深度优先搜索(DFS)找出N1和N2之间的无向路径P,随机选取P中的任意一条支路,将其删除,同时在T中增加支路N1N2,从而产生新个体T*;反之,则重新选择节点N1和N2,变异操作的流程图和示意图见附图5、6。
输出模块给出开断支路情况,输出网络重构结构。
实施例2
本实例将一种采用基于无向生成树的并行遗传算法的配网重构方法应用于IEEE33节点标准配电系统以验证该方法的有效性。IEEE 33节点系统的电压等级为12.66kV,总有功负荷、无功负荷分别为3715kW和2300kVar,网络拓扑如附图7所示。假定系统各节点的负荷(有功负荷、无功负荷)均服从方差为0.1μL正态分布,其中μL为负荷期望。采用丹麦Bonus 1MW/54型风机,技术参数如表1所示;光伏阵列采用Pilkington SFM144Hx250wp型太阳能电池组件,其技术参数如表2所示。假定分别在节点25和节点30处安装2台风机,在节点14处安装4个光伏阵列。节点电压允许范围为0.93~1.07p.u.,节点1与节点6间支路有功潮流上限为5.6MW,其余支路有功潮流上限为3MW。
表1风机技术参数
表2光伏阵列技术参数
PSTGA算法参数为:选择概率0.8、交叉概率0.9、变异概率0.4、染色体总数100条、子种群数4个、最大迭代次数为50次。LHS采样规模为100次,待采随机变量数共69个,包含2个风速变量、1个太阳辐照强度变量、33个节点有功负荷变量以及33个节点无功负荷变量。本文测试环境为2011b,所使用的计算机CPU为AMD AthlonTMII X4-635,主频2.91GHz,内存为2.00GB。
表3给出了未安装DG时IEEE 33节点系统重构结果,可见重构前后节点电压有了较明显的改善,系统网损也随之显著降低,表明在未安装任何分布式电源的情况下,重构对于改善电压、降低网损有一定积极意义。
表3IEEE 33节点系统重构结果(未安装DG)
针对含DG情况,本文取置信度βU和βI为相同数值,分别给出了置信度为0.90,0.95,0.99时IEEE 33节点系统重构结果,如表4所示。由表4可见,不同置信度下的重构方案相同,之所以会出现这种情况是因为在该重构方案下系统中各节点电压均不越限、各支路潮流也不存在越限,而且系统总网损期望值最小。
表4不同置信度下IEEE 33节点系统重构结果(安装DG)
附图8给出系统各节点在有无DG(DWG/PV)时重构前后的电压分布情况,可以看出:①未安装DG时,系统中较多节点存在电压越下限的情况;安装DG以后,系统节点电压越限情况稍有改善,但是未完全消除(如节点18);②未安装DG时,网络重构对于改善节点电压有一定作用,但是在负荷波动的情况下可能存在部分节点电压越限的问题;在安装DG的基础上进行网络重构能够更大程度地改善节点电压分布,尤其是馈线末端节点处的电压。
针对DG接入情况下馈线末端节点可能存在电压越限问题,本文选取了部分有代表性的节点,分析了网络重构对上述节点处电压概率分布的影响,如附图9和附图10所示。可以看出,网络重构前节点18和节点33均存在不同程度的电压越限可能;网络重构后节点18和节点33处节点电压均处于正常范围。这充分说明在DG接入情况下网络重构能有效改善馈线末端节点电压分布。
为验证本发明算法性能,从稳定性、收敛速度、计算耗时三个角度,通过与基于无向生成树的遗传算法(STGA)、粒子群算法(PSO)、改进和声搜索算法(IHSA)、蚁群搜索算法(ACS)等算法的比较,以验证所提出的PSTGA算法的有效性。其中,PSTGA算法参数同前;STGA算法参数:染色体总数100条、最大迭代次数50次、选择概率0.8、交叉概率0.9、变异概率0.4;PSO算法参数:粒子总数100个、最大迭代次数50、加速常数c1=c2=2、惯性权重w∈[0.8,1.2]、飞翔速度v∈[-4,4];IHSA算法参数:RHM=8,RHMmin=0.7,RHMmax=0.9,NI=300、ACS算法参数设置:迭代次数100,蚁群M=8,α=0.3,β=0.6,Δτij(0)=0,τmax=τij(0),τmin=τmax/1024。
(1)稳定性性能分析
为了测试PSTGA算法的稳定性,本文分别使用PSTGA算法与其他算法求解不含DG情况下IEEE 33节点系统重构问题(记为情形1,下同),并在一定次数的独立重复试验的基础上,对结果进行了统计处理,如表5、附图11所示。可见,PSTGA算法具有较高的寻优率,稳定性明显优于改进PSO算法、IHSA算法,但稍逊于STGA算法;PSTGA算法得到的最优解与改进PSO算法、ACS算法给出的结果相同,且优于IHSA算法的优化结果。IHSA算法计算结果见陈春、汪沨、刘蓓等在《电力系统自动化》(2014,38(6):55-60)上发表的《基于基本环矩阵与改进和声搜索算法的配电网重构》;ACS算法计算结果见Wu Y K、Lee C Y、Liu L C等在《IEEETransactions on Power Delivery》(2010,25(3):1678-1685)上发表的《Study ofreconfiguration for the distribution system with distributed generators》;改进PSO算法计算结果见卢志刚、杨国良、张晓辉等在《电力系统保护与控制》(2009,37(7):30-34)上发表的《改进二进制粒子群优化算法在配电网络重构中的应用》。
表5PSTGA、STGA以及其他算法稳定性比较
注:最优值(网损)统一采用相同的潮流计算方法计算;ACS算法文献并未给出稳定性测试结果,故略去;寻优率为收敛到全局最优解次数与试验次数的比值
(2)收敛速度性能分析
同样,在情形1下,本文分别对PSTGA、STGA和PSO的收敛速度进行测试,并且独立重复试验次数为50,结果如表6和附图12所示。改进PSO算法和STGA算法收敛速度较快,而PSTGA收敛速度稍慢。
表6PSTGA、STGA和PSO收敛速度比较
(3)并行效果性能分析
为了测试PSTGA算法的并行效果,本文针对情形1(同上)和情形2(DG接入下考虑机会约束的IEEE 33节点系统网络重构问题,置信度为0.90)分别进行了多次独立重复试验,统计结果见表7。对于情形1,和STGA算法相比PSTGA算法的耗时减少了约36%,并未达到预期目标(耗时减少75%),这主要是以下两方面因素综合的结果:①情形1计算规模较小;②并行计算本身存在一定的时间开销。对于情形2,鉴于采用了基于拉丁超立方采样的蒙特卡洛模拟法,计算规模较大,因而PSTGA算法的并行优势体现得比较明显,耗时约为STGA算法的28.6%。总之,PSTGA算法因采用了并行方式,计算用时较STGA算法有不同程度的减少,具体减少程度应视计算规模而定。
表7并行效果测试结果
通过本实施例的验证,可知一种采用基于无向生成树的并行遗传算法的配网重构方法有效解决了DG接入后的配电网重构问题;在DG接入情况下进行配网重构有助于降低系统网损,有效改善馈线末端节点电压分布,缓解节点电压越限;同时该方法相比以往的几种算法,如STGA、PSO、IHSA、ACS等具有不错的稳定性、较好的收敛速度和高效的并行能力。
Claims (10)
1.一种采用基于无向生成树的并行遗传算法的配网重构方法,其特征在于,该方法包括下列步骤:
(1)获取配网参数、风机和光伏参数以及并行遗传算法参数;
(2)根据步骤(1)中获取的参数对风速、辐射强度和负荷进行基于拉丁超立方的蒙特卡洛模拟采样;
(3)随机产生包含Nind个拓扑可行的染色体的初始种群作为目标种群,并设置迭代次数n的初值为1;
(4)根据步骤(2)中得到的采样值,对目标种群中的每个染色体所对应的网络拓扑进行潮流计算,得到系统网损期望、节点电压及支路潮流越限概率;
(5)根据步骤(4)中潮流计算得到的结果计算目标函数值,判断是否满足约束条件,若是则进入步骤(6),若否则返回步骤(3);
(6)对现有种群分为多个子种群,进行基于无向生成树的并行遗传操作,得到多个目标子种群;
(7)产生一个1~Nsub的随机排列P=[p1,p2,…,pNsub],建立目标子种群i与源子种群pi间的映射关系;
(8)将各个目标子种群的最差个体替换为相应的源子种群的最优个体,合并各个目标子种群作为目标种群;
(9)判断迭代次数n是否达到要求,若是,则进入步骤(10),若否,则迭代次数n增加1,返回步骤(4);
(10)根据得到的目标种群输出配电网重构方案。
2.根据权利要求1所述的一种采用基于无向生成树的并行遗传算法的配网重构方法,其特征在于,所述的配网参数包括网络拓扑结构、各节点负荷、网络支路阻抗及导纳、网络电压等级;所述的风机参数包括切入风速、额定风速、切出风速、额定功率、功率因数、Weibull分布参数;所述的输入模块中的光伏参数包括光伏面板面积、光电效率、额定功率、功率因数、Beta分布参数。
3.根据权利要求2所述的一种采用基于无向生成树的并行遗传算法的配网重构方法,其特征在于,所述的Weibull分布用以描述风电场风速概率分布,概率密度函数可表示为:
其中,v为风速,c和k分别为Weilbull分布的尺度参数和形状参数;
风机出力与风速的关系可表示为:
其中,vci、vcr和vco分别为切入风速、额定风速和切出风速;Pr为风机额定输出有功功率。
4.根据权利要求2所述的一种采用基于无向生成树的并行遗传算法的配网重构方法,其特征在于,所述的Beta分布用以描述光照辐射强度概率分布,其概率密度可表示为:
其中,α和β为Beta分布的形状参数;r和rmax分别为某一时段内的太阳辐照强度及其最大值;
光伏输出有功功率可表示为:
PM=rMAη
其中,PM为光伏阵列输出有功功率;M为光伏阵列中太阳能电池组件个数;A为每个太阳能电池组件的面积;η为光伏阵列的光电转换效率。
5.根据权利要求1所述的一种采用基于无向生成树的并行遗传算法的配网重构方法,其特征在于,所述的基于拉丁超立方的蒙特卡洛模拟采样具体步骤为:
(21)对每个输入变量进行采样,确保随机分布区域能够被采样点完全覆盖;
(22)改变变量采样值的排列顺序,使相对独立的随机变量采样值的相关性趋于最小。
6.根据权利要求1所述的一种采用基于无向生成树的并行遗传算法的配网重构方法,其特征在于,所述的目标函数为:
式中,PLoss为系统网损期望值,;E{·}表示计算随机变量的数学期望;NL为系统支路数;Pk,Loss为随机变量,对应于第k条支路的网损;NLHS为LHS采样次数;Ii,k为第i次采样中流过第k条支路的潮流,其中i=1,2,…,NLHS;rk为第k条支路电阻,其中k=1,2,…,NL。
7.根据权利要求1所述的一种采用基于无向生成树的并行遗传算法的配网重构方法,其特征在于,所述的约束条件包括:
1)潮流方程约束
其中,Pis和Qis分别为节点i处注入有功功率和无功功率;Vi为节点i处电压幅值;Ωi表示与节点i相连的其他节点集合;Gij、Bij分别为支路ij的电导和电纳;θij为节点i与节点j间相角差,其中i,j=1,2,...,n,n为节点总数;
2)节点电压约束
其中,Pr{·}为事件成立概率;和分别为节点i处电压幅值的上下限;βV为电压约束的置信水平;Θ为系统节点集合;
3)支路潮流约束
其中,为第j条支路传输有功功率的上限;βI为支路潮流约束的置信水平;Λ为系统支路集合;
4)网络拓扑约束
g∈G
其中,g为当前网络拓扑结构;G为所有满足辐射状结构的网络拓扑集合。
8.根据权利要求1所述的一种采用基于无向生成树的并行遗传算法的配网重构方法,其特征在于,所述的基于无向生成树的并行遗传操作包括选择操作、交叉操作、变异操作。
9.根据权利要求8所述的一种采用基于无向生成树的并行遗传算法的配网重构方法,其特征在于,所述的选择操作采用轮盘赌方法;所述的交叉操作具体步骤为:
(61)随机选择节点N1与N2,使得N1≠N2;
(62)对于父代生成树T1与T2,通过深度优先搜索确定N1与N2之间的无向路径,其中T1中的路径为P1,T2中的路径为P2;
(63)分别将T1与P2、T2与P1进行交叉操作,得到新个体T1 *与T2 *。
10.根据权利要求8所述的一种采用基于无向生成树的并行遗传算法的配网重构方法,其特征在于,所述的变异操作具体步骤为:
(64)获取父代拓扑T;
(65)随机选择节点N1与N2,使得N1≠N2;
(66)判断T中是否存在以N1和N2为端点的支路,若否,则进入步骤(67),若是,则返回步骤(65)。
(67)通过深度优先搜索找出N1和N2之间的无向路径P,随机选取P中的任意一条支路,将其删除,同时在T中增加支路N1N2,从而产生新个体T*。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610887323.6A CN106340873A (zh) | 2016-10-11 | 2016-10-11 | 一种采用基于无向生成树的并行遗传算法的配网重构方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610887323.6A CN106340873A (zh) | 2016-10-11 | 2016-10-11 | 一种采用基于无向生成树的并行遗传算法的配网重构方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106340873A true CN106340873A (zh) | 2017-01-18 |
Family
ID=57838720
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610887323.6A Pending CN106340873A (zh) | 2016-10-11 | 2016-10-11 | 一种采用基于无向生成树的并行遗传算法的配网重构方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106340873A (zh) |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107066709A (zh) * | 2017-03-29 | 2017-08-18 | 西安电子科技大学 | 基于多智能体遗传算法的电力运输网络拓扑结构设计方法 |
CN107332290A (zh) * | 2017-08-30 | 2017-11-07 | 国网江苏省电力公司南京供电公司 | 一种基于直流线路的区域负荷转供方法 |
CN107979092A (zh) * | 2017-12-18 | 2018-05-01 | 国网宁夏电力有限公司经济技术研究院 | 一种考虑分布式电源与软开关接入的配电网动态重构方法 |
CN108334696A (zh) * | 2018-02-02 | 2018-07-27 | 国网江西省电力有限公司经济技术研究院 | 一种考虑功率随机性的配电网日前网络重构方法 |
CN110033405A (zh) * | 2019-03-20 | 2019-07-19 | 上海交通大学 | 一种基于图像处理器加速的符号执行方法及系统 |
CN110705813A (zh) * | 2019-07-23 | 2020-01-17 | 电子科技大学 | 一种考虑可靠性风电场集电系统的混合线缆连接方法 |
CN111083577A (zh) * | 2019-10-18 | 2020-04-28 | 国网浙江省电力有限公司嘉兴供电公司 | 一种电力通信拓扑结构优化方法 |
CN111082470A (zh) * | 2020-01-15 | 2020-04-28 | 合肥工业大学 | 含低风速分散式风电的配电网多目标动态鲁棒重构方法 |
CN111628528A (zh) * | 2020-06-22 | 2020-09-04 | 华北电力大学(保定) | 风电参与系统恢复的潮流越限消除方法及装置 |
CN115828489A (zh) * | 2023-02-22 | 2023-03-21 | 中国电力科学研究院有限公司 | 基于关键量测布点位置搜索的感知设备部署方法及系统 |
CN116693163A (zh) * | 2023-07-21 | 2023-09-05 | 耀昶嵘相变材料科技(广东)有限公司 | 污泥干化系统的控制方法、终端及系统 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103903055A (zh) * | 2014-03-21 | 2014-07-02 | 国家电网公司 | 基于无向图所有生成树的网络重构遗传算法 |
CN104269886A (zh) * | 2014-10-16 | 2015-01-07 | 东北电力大学 | 用于风电并网系统容量效益裕度计算及分配的方法 |
CN104361405A (zh) * | 2014-10-28 | 2015-02-18 | 广东电网有限责任公司电力科学研究院 | 一种基于容量限值约束的微网储能装置设计方法 |
CN105741193A (zh) * | 2016-04-20 | 2016-07-06 | 河海大学 | 计及分布式电源和负荷不确定性的多目标配网重构方法 |
-
2016
- 2016-10-11 CN CN201610887323.6A patent/CN106340873A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103903055A (zh) * | 2014-03-21 | 2014-07-02 | 国家电网公司 | 基于无向图所有生成树的网络重构遗传算法 |
CN104269886A (zh) * | 2014-10-16 | 2015-01-07 | 东北电力大学 | 用于风电并网系统容量效益裕度计算及分配的方法 |
CN104361405A (zh) * | 2014-10-28 | 2015-02-18 | 广东电网有限责任公司电力科学研究院 | 一种基于容量限值约束的微网储能装置设计方法 |
CN105741193A (zh) * | 2016-04-20 | 2016-07-06 | 河海大学 | 计及分布式电源和负荷不确定性的多目标配网重构方法 |
Non-Patent Citations (2)
Title |
---|
于晗等: "采用拉丁超立方采样的电力系统概率潮流计算方法", 《电力系统自动化》 * |
黄红程等: "基于无向生成树的并行遗传算法在配电网重构中的应用", 《电力系统自动化》 * |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107066709B (zh) * | 2017-03-29 | 2019-05-21 | 西安电子科技大学 | 基于多智能体遗传算法的电力运输网络拓扑结构设计方法 |
CN107066709A (zh) * | 2017-03-29 | 2017-08-18 | 西安电子科技大学 | 基于多智能体遗传算法的电力运输网络拓扑结构设计方法 |
CN107332290A (zh) * | 2017-08-30 | 2017-11-07 | 国网江苏省电力公司南京供电公司 | 一种基于直流线路的区域负荷转供方法 |
CN107332290B (zh) * | 2017-08-30 | 2023-06-06 | 国网江苏省电力公司南京供电公司 | 一种基于直流线路的区域负荷转供方法 |
CN107979092A (zh) * | 2017-12-18 | 2018-05-01 | 国网宁夏电力有限公司经济技术研究院 | 一种考虑分布式电源与软开关接入的配电网动态重构方法 |
CN108334696B (zh) * | 2018-02-02 | 2021-06-08 | 国网江西省电力有限公司经济技术研究院 | 一种考虑功率随机性的配电网日前网络重构方法 |
CN108334696A (zh) * | 2018-02-02 | 2018-07-27 | 国网江西省电力有限公司经济技术研究院 | 一种考虑功率随机性的配电网日前网络重构方法 |
CN110033405A (zh) * | 2019-03-20 | 2019-07-19 | 上海交通大学 | 一种基于图像处理器加速的符号执行方法及系统 |
CN110705813B (zh) * | 2019-07-23 | 2022-11-22 | 电子科技大学 | 一种考虑可靠性风电场集电系统的混合线缆连接方法 |
CN110705813A (zh) * | 2019-07-23 | 2020-01-17 | 电子科技大学 | 一种考虑可靠性风电场集电系统的混合线缆连接方法 |
CN111083577A (zh) * | 2019-10-18 | 2020-04-28 | 国网浙江省电力有限公司嘉兴供电公司 | 一种电力通信拓扑结构优化方法 |
CN111082470A (zh) * | 2020-01-15 | 2020-04-28 | 合肥工业大学 | 含低风速分散式风电的配电网多目标动态鲁棒重构方法 |
CN111082470B (zh) * | 2020-01-15 | 2022-09-02 | 合肥工业大学 | 含低风速分散式风电的配电网多目标动态鲁棒重构方法 |
CN111628528A (zh) * | 2020-06-22 | 2020-09-04 | 华北电力大学(保定) | 风电参与系统恢复的潮流越限消除方法及装置 |
CN115828489A (zh) * | 2023-02-22 | 2023-03-21 | 中国电力科学研究院有限公司 | 基于关键量测布点位置搜索的感知设备部署方法及系统 |
CN116693163A (zh) * | 2023-07-21 | 2023-09-05 | 耀昶嵘相变材料科技(广东)有限公司 | 污泥干化系统的控制方法、终端及系统 |
CN116693163B (zh) * | 2023-07-21 | 2023-12-05 | 耀昶嵘相变材料科技(广东)有限公司 | 污泥干化系统的控制方法、终端及系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106340873A (zh) | 一种采用基于无向生成树的并行遗传算法的配网重构方法 | |
Ehsan et al. | Optimal integration and planning of renewable distributed generation in the power distribution networks: A review of analytical techniques | |
CN108306303B (zh) | 一种考虑负荷增长和新能源出力随机的电压稳定评估方法 | |
CN104578157B (zh) | 一种分布式电源接入电网的潮流计算方法 | |
CN103903073B (zh) | 一种含分布式电源及储能的微电网优化规划方法及系统 | |
Yang et al. | Modelling, applications, and evaluations of optimal sizing and placement of distributed generations: A critical state‐of‐the‐art survey | |
CN103904644B (zh) | 一种基于分布式电源接入的智能变电站负荷自动分配方法 | |
CN104734153A (zh) | 一种含分布式电源的配电网络重构方法 | |
CN106532778A (zh) | 一种计算分布式光伏并网最大准入容量的方法 | |
CN107069814A (zh) | 配网分布式电源容量布点的模糊机会约束规划方法与系统 | |
CN108023364A (zh) | 基于凸差规划的配电网分布式电源最大接入能力计算方法 | |
Saba et al. | Development of an ontology-based generic optimisation tool for the design of hybrid energy systems | |
CN105186499A (zh) | 一种配电网多目标概率最优潮流模糊建模与求解方法 | |
CN105354655A (zh) | 计及功率相关性的光伏电站群置信容量评估方法 | |
Askarzadeh | Optimisation of solar and wind energy systems: A survey | |
Liu et al. | Research on short-term optimization for integrated hydro-PV power system based on genetic algorithm | |
CN106532710A (zh) | 计及电压稳定约束的微电网潮流优化方法 | |
Javed et al. | A comparative study of maximum power point tracking techniques for solar systems | |
Valarmathi et al. | Simulation and analysis of wind energy and photo voltaic hybrid system | |
CN110380419A (zh) | 一种混合随机和区间变量的不确定潮流计算方法 | |
Ye et al. | Multi-scenario stochastic optimal scheduling for power systems with source-load matching based on pseudo-inverse Laguerre polynomials | |
CN108233373B (zh) | 考虑天气场景的分布式光伏接入配电网的概率谐波分析方法 | |
Sayedin et al. | Evolutionary optimization approaches for direct coupling photovoltaic-electrolyzer systems | |
CN105262107A (zh) | 一种风电并网的无功优化方法 | |
Mohamed | Optimal Type-2 Fuzzy Logic Control for Sizing and Control Grid Connected PV System with Battery |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170118 |