CN115423242A - 一种基于配电网重构的供电能力挖掘方法 - Google Patents

一种基于配电网重构的供电能力挖掘方法 Download PDF

Info

Publication number
CN115423242A
CN115423242A CN202210840243.0A CN202210840243A CN115423242A CN 115423242 A CN115423242 A CN 115423242A CN 202210840243 A CN202210840243 A CN 202210840243A CN 115423242 A CN115423242 A CN 115423242A
Authority
CN
China
Prior art keywords
power
mayflies
power supply
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.)
Pending
Application number
CN202210840243.0A
Other languages
English (en)
Inventor
燕跃豪
鲍薇
吴松梅
刘晓辉
王彬楠
陈启龙
王莹
刘雪珂
石云松
郑琰
王建新
蒋建东
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
State Grid Henan Electric Power Co Zhengzhou Power Supply Co
Zhengzhou University
State Grid Corp of China SGCC
Original Assignee
State Grid Henan Electric Power Co Zhengzhou Power Supply Co
Zhengzhou University
State Grid Corp of China SGCC
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by State Grid Henan Electric Power Co Zhengzhou Power Supply Co, Zhengzhou University, State Grid Corp of China SGCC filed Critical State Grid Henan Electric Power Co Zhengzhou Power Supply Co
Priority to CN202210840243.0A priority Critical patent/CN115423242A/zh
Publication of CN115423242A publication Critical patent/CN115423242A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0631Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • G06Q10/06313Resource planning in a project environment
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/004Artificial life, i.e. computing arrangements simulating life
    • G06N3/006Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06Energy or water supply
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Physics & Mathematics (AREA)
  • Economics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Strategic Management (AREA)
  • General Business, Economics & Management (AREA)
  • General Health & Medical Sciences (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Tourism & Hospitality (AREA)
  • Marketing (AREA)
  • Biomedical Technology (AREA)
  • General Engineering & Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Computational Linguistics (AREA)
  • Biophysics (AREA)
  • Public Health (AREA)
  • Water Supply & Treatment (AREA)
  • Mathematical Physics (AREA)
  • Evolutionary Computation (AREA)
  • Primary Health Care (AREA)
  • Artificial Intelligence (AREA)
  • Software Systems (AREA)
  • Computing Systems (AREA)
  • Biodiversity & Conservation Biology (AREA)
  • Development Economics (AREA)
  • Educational Administration (AREA)
  • Molecular Biology (AREA)
  • Game Theory and Decision Science (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Power Engineering (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明公开了一种基于配电网重构的供电能力挖掘方法,通过计算供电系统的最大供电能力来指导电网的规划,改变配电网中联络开关和分段开关的状态,引入分布式电源的方法,将配电网络的拓扑结构改变,馈线和变电站之间的平均负荷分配实现,提高负荷供应能力;同时,采用蜉蝣算法计算最大供电能力,拓宽问题的求解方法。本发明以网络潮流为约束条件,通过计算系统的最大供电能力,提升配电网的负荷供应能力,将最新的蜉蝣算法应用到最大供电能力的计算中,拓宽问题的求解方法;本发明将蜉蝣算法应用于配电网供电能力挖掘,通过算例模型的计算,得到蜉蝣算法,具有更好的收敛精度和收敛速度,在保证经济性的同时,有效的提升了配电网的供电能力。

Description

一种基于配电网重构的供电能力挖掘方法
技术领域:
本发明属于电力工程技术领域,具体涉及一种基于配电网重构的供电能力挖掘方法。
背景技术:
配电系统是联系用户和输电系统的纽带,是保证可靠供电的重要环节。配电系统的安全、经济、可靠运行,直接影响到人民生活水平质量的提高、城市经济的发展以及投资环境的优化。配电网直接面向用户,其最直接的目标是最大限度地满足供电区域内用户的用电需求。目前深入研究配电系统的需要已日趋显现。用户对可靠性要求越来越高,但大部分影响可靠性的停电都在配电系统。配电系统是城市基础设施,与发电系统、输电系统同样对供电安全具有重大影响。配电投资增幅最快,网络损耗最大的环节也在配电网。智能电网下分布式能源和微电网的接入在配电系统,智能电网相对传统电网变革最大的环节也在配电网。因此无论从安全、可靠、经济还是发展角度,配电地位都上升到与发、输电相同的高度。特别是智能电网下很多运行条件将发生重大变化,配电系统将实现信息化,配电自动化将普及。而且随着我国国民经济持续健康的发展,整个社会对电能的需求也越来越大,为了满足广大电力用户对电能的需求,配电系统的规模必须不断扩大。随着城乡电网的大规模建设改造,电力系统面临着日益繁重的规划任务,这对配电网提出了更高的要求,供电能力逐渐成为评价配电网建设水平的一个新指标。
城市电网是包括了220kV电压等级的输电网、110kV的高压配电网及以下电压等级配电网的电力系统。由于服务范围广,负荷分布的空间差异较大,负荷高峰时期局部220kV系统容易出线阻塞现象。因此工程上调度员常常通过改变110kV供电路径的主备关系,改变潮流分布,缓解输电阻塞。故分析城市片区电网供电能力时,有必要计及高压配电网重构能力对于腾挪供电空间、挖掘供电潜能的积极作用。
发明内容:
本发明所要解决的技术问题是:提供一种基于配电网重构的供电能力挖掘方法,通过计算供电系统的最大供电能力来指导电网的规划,改变配电网中联络开关和分段开关的状态,引入分布式电源的方法,将配电网络的拓扑结构改变,馈线和变电站之间的平均负荷分配实现,提高负荷供应能力;同时,采用蜉蝣算法计算最大供电能力,拓宽问题的求解方法,在配电网发生事故时,基于N-1安全准则,故障点将迅速消除,供电恢复,减少停电区域,减少停电引起的生产损失。
本发明为解决技术问题所采取的技术方案是:
一种基于配电网重构的供电能力挖掘方法,通过计算供电系统的最大供电能力来指导电网的规划,改变配电网中联络开关和分段开关的状态,引入分布式电源的方法,将配电网络的拓扑结构改变,馈线和变电站之间的平均负荷分配实现,提高负荷供应能力;同时,采用蜉蝣算法计算最大供电能力,拓宽问题的求解方法,具体包括以下几个步骤:
A、收集配电网的信息;
B、根据收集的配电网信息,构建具体的配电网网络;
C、根据配电网网络构建基于配电网重构的最大供电能力计算的数学模型:
第一步:建立目标函数:
Figure BDA0003750660340000031
式中:F为目标函数,Lp,i为线路i上的有功功率,ΦN代表所有的总线集和,N表示配电网络有N条线路;
目标是使系统中所有负荷的有功功率最大化,为了简化,设定负载的功率因数为常数(0.9),当确定了同一母线上的有功需求时,该线路自动获得所需的无功功率;
第二步:提供网络的约束条件:
Figure BDA0003750660340000032
式中:f代表变压器出线故障的场景,xi,j表示分支i-j的电抗,Nmode表示不包含根总线的所有总线数量,B表示与故障变压器直接连接的一组支路;
将分支的状态描述为Boolean值,对于径向拓扑,所有连接的线路数等于节点数减去线路条数;第二行的公式表示故障变压器支路i-j由于偶然事故而被断开;在分析网状网络的过程中,将
Figure BDA0003750660340000033
表示为
Figure BDA0003750660340000034
其中K代表独立环路的个数,i-j∈f表示分支i-j与变压器fth相连;
dij=-dij.dij∈{-1,1} (3)
式中dij描述了分支的正方向,其在计算模型之前自定义,方向参数在优化前是恒定且已知的;
Figure BDA0003750660340000041
式中:LP,i和LQ,i为总线i处的有功负荷和无功负荷,Pij,i和Qij,i为支路i-j的终端i的有功电源和无功电源,
Figure BDA00037506603400000410
为功率因数;
上式描述了功率平衡的约束条件:节点注入的有功功率等于连接该节点的各支路的功率之和,无功功率由有功功率和相应的功率因数求得;
分布式发电机DG母线的功率输出约束:
Figure BDA0003750660340000042
式中:
Figure BDA0003750660340000043
Figure BDA0003750660340000044
分别为连接到总线i上的DG的最大有功和无功功率输出,N(i)为通过分支连接到总线i的总线,ΦDG为连接有分布式电源的总线集;
各支路的热额定值约束:
Figure BDA0003750660340000045
式中:
Figure BDA0003750660340000046
为分支i-j的功率限制;
在计算的过程中采用一次变量代替二次变量,如下式:
Figure BDA0003750660340000047
式中:
Figure BDA0003750660340000048
为母线i的电压的平方,
Figure BDA0003750660340000049
为分支i-j的电流的平方,连接支路两端节点间的电压关系表示为:
Figure BDA0003750660340000051
式中:Pij,i、Qij,i和Xij,i分别为分支i-j的终端i的有功功率、无功功率和电抗,Xij为分支i-j的电抗;
将上述式子中的节点电压限定在上下边界内:
Figure BDA0003750660340000052
式中:
Figure BDA0003750660340000053
Figure BDA0003750660340000054
分别为母线i的电压大小的平方的上边界和下边界,支路电流的表达式为:
Figure BDA0003750660340000055
无功损耗和有功损耗的表达式为:
Figure BDA0003750660340000056
变压器的额定功率约束为:
Figure BDA0003750660340000057
式中:
Figure BDA0003750660340000058
为变压器的额定值,如果变压器发生故障且故障在考虑的场景之一被隔离,则该变压器的额定值为0,式中
Figure BDA0003750660340000059
为电压的平方根;
上述式子中f的取值如式(13)所示,表示系统共有的几种情况,包含几种单个变压器故障的情况和一种正常情况,用f=0表示正常情况,即所有的变压器均无故障:
f=0,1,2,K Nt (13)
由于未断开支路的两端不应有式(8)中所描述的电压关系,所以式(8)不能用于未断开支路,采用下式代替:
Figure BDA0003750660340000061
当分支i-j连接时
Figure BDA0003750660340000062
等于1,
Figure BDA0003750660340000063
等于0,因此,式(14)等于式(8);当分支i-j断开时,
Figure BDA0003750660340000064
等于0,
Figure BDA0003750660340000065
等于M0;因此,
Figure BDA0003750660340000066
的取值范围为[-M0,M0];
为了解决式(10)中不同变量的乘积模型不可解的问题,通过SOCP松弛使用下式代替式(10):
Figure BDA0003750660340000067
Figure BDA0003750660340000068
其中式(15)为网络损耗最小化,SOCP松弛用来凸化TSC编程模型,修正后的模型表示为:
Figure BDA0003750660340000069
s,t,(2)-(7),(9),(11)-(14),(16)
根据相应的网络配置,利用上述约束条件求解功率流,对任何变压器故障而言,系统中的所有负载都在满足约束的情况下为系统提供功率;系统中的变量和约束条件的总集和具有相同的目标,即考虑网络所能提供总负荷的最大化,该模型考虑了潮流求解约束,通过网络重构来挖掘整个网络对系统负荷的供电能力。
所述的式(17)给出的模型采用蜉蝣算法求解实际配电网考虑重构的最大供电能力,具体的步骤如下:
第一步:随机生成雄性蜉蝣和雌性蜉蝣种群,初始化参数;
第二步:根据步骤C中的目标函数计算各个蜉蝣的适应度值并排序,得到个体最优和群体最优;
第三步:雄性蜉蝣的运动:
雄性蜉蝣一般成群结队活动,这就意味着每个雄虫的位置是根据自己和它的邻居的经验来调整的,假设
Figure BDA0003750660340000071
是第i个蜉蝣在t时刻搜索空间中的位置,蜉蝣的位置是通过在当前位置上增加一个速度
Figure BDA0003750660340000072
来改变的,因此,蜉蝣的第t+1时刻的位置表达式为:
Figure BDA0003750660340000073
考虑到雄性蜉蝣总是在离水面几米高的地方进行舞蹈,因此他们的速度并不会很高,所以,一只雄性蜉蝣i的速度计算为:
Figure BDA0003750660340000074
式中:
Figure BDA0003750660340000075
为第i个雄性蜉蝣在搜索空间j=1,…,n中t时刻的速度,
Figure BDA0003750660340000076
表示t时刻的位置,a1和a2分别为衡量规模认知和社会贡献的积极吸引常数,pbesti表示雄性蜉蝣i在搜索空间中的最好位置,其表达式为:
Figure BDA0003750660340000077
式中:F()为目标函数,同步骤C中的式子(1)所示;
对于算法性能来说,群体中的最优雄性蜉蝣继续执行其上下舞蹈,因此,最优的雄性蜉蝣不断更新他们的速度,其表达式为:
Figure BDA0003750660340000081
其中d为舞蹈系数,r为[-1,1]范围内的随机值;
第四步:雌性蜉蝣的运动:
与雄性蜉蝣不同,雌雄蜉蝣不会有聚集性活动,它们会飞向雄性进行繁殖,假设
Figure BDA0003750660340000082
是t时刻雌性蜉蝣在搜索空间中的位置,同时通过添加速度
Figure BDA0003750660340000083
来改变其位置,即:
Figure BDA0003750660340000084
考虑雄性蜉蝣和雌性蜉蝣相互吸引,因此,它们的速度计算公式为:
Figure BDA0003750660340000085
式中:
Figure BDA0003750660340000086
为第i个雌性蜉蝣在搜索空间j=1,…,n中t时刻的速度,
Figure BDA0003750660340000087
表示t时刻的位置;a2为相互之间的吸引常数,β是一个固定的可见性系数;fl是一个随机游走系数,当雌性不被雄性吸引时使用;r为[-1,1]之间的随机数;rmf为雌性和雄性蜉蝣之间的笛卡尔距离,其表达式为:
Figure BDA0003750660340000088
第五步:蜉蝣的交配:
每对蜉蝣交配后产生两个后代,表达式为:
offspring1=L*male+(1-L)*female (25)
offspring2=L*female+(1-L)*male (26)
式中male为父系,female为母系;L为特定范围内的随机数,后代的初始速度设定为0;
第六步:计算适应度值并更新个体最优和群体最优;
第七步:判断是否满足停止田间,若满足,则停止迭代,输出最优结果,否则转到第三步继续执行。
随着城乡电网的大规模建设改造,电力系统面临着日益繁重的规划任务,这对配电网提出了更高的要求,对配电网供电能力合理挖掘能够最大程度利用配电网现有资源,有助于对配电网进行合理投资。
以往的配电网优化规划基本只考虑了网损、无功优化等问题,而没有考虑供电能力的提升,经济性也得不到有效保证。由于配电网复杂的电路设计和负荷较高,故障的发生可能性更大。
本发明的一种基于配电网重构的供电能力挖掘方法,通过计算供电系统的最大供电能力来指导电网的规划,改变配电网中联络开关和分段开关的状态,引入分布式电源的方法,将配电网络的拓扑结构改变,馈线和变电站之间的平均负荷分配实现,提高负荷供应能力;同时,采用蜉蝣算法计算最大供电能力,拓宽问题的求解方法,在配电网发生事故时,基于N-1安全准则,故障点将迅速消除,供电恢复,减少停电区域,减少停电引起的生产损失。
本发明的积极有益效果如下:
1、本发明考虑配电网重构的供电能力挖掘能够最大程度的利用配电网现有资源,相比不考虑配电网重构,能够为用户提供更大的供电负荷。在保证经济性的同时,有效的提升了配电网的供电能力。
2、本发明将蜉蝣算法应用于配电网供电能力挖掘,通过算例模型的计算,得到蜉蝣算法,相比粒子群算发和遗传算法拥有更好的收敛精度和收敛速度。
四、附图说明:
图1为本发明实施例中的算例网络图;
图2是本发明实施例中采用的标砖IEEE-33测试配电网;
图3是本发明实施例中采用蜉蝣算法的收敛图;
图4是本发明实施例中配网最大供电能力计算的算法流程图。
五、具体实施方式:
下面结合附图和具体实施例对本发明做进一步的解释和说明:
一种基于配电网重构的供电能力挖掘方法,通过计算供电系统的最大供电能力来指导电网的规划,改变配电网中联络开关和分段开关的状态,引入分布式电源的方法,将配电网络的拓扑结构改变,馈线和变电站之间的平均负荷分配实现,提高负荷供应能力;同时,采用蜉蝣算法计算最大供电能力,拓宽问题的求解方法,具体包括以下几个步骤:
A、收集配电网的信息;
B、根据收集的配电网信息,构建具体的配电网网络;
C、根据配电网网络构建基于配电网重构的最大供电能力计算的数学模型:
第一步:建立目标函数:
Figure BDA0003750660340000101
式中:F为目标函数,Lp,i为线路i上的有功功率,ΦN代表所有的总线集和,N表示配电网络有N条线路;
目标是使系统中所有负荷的有功功率最大化,为了简化,设定负载的功率因数为常数(0.9),当确定了同一母线上的有功需求时,该线路自动获得所需的无功功率;
第二步:提供网络的约束条件:
Figure BDA0003750660340000111
式中:f代表变压器出线故障的场景,xi,j表示分支i-j的电抗,Nmode表示不包含根总线的所有总线数量,B表示与故障变压器直接连接的一组支路;
将分支的状态描述为Boolean值,对于径向拓扑,所有连接的线路数等于节点数减去线路条数;第二行的公式表示故障变压器支路i-j由于偶然事故而被断开;在分析网状网络的过程中,将
Figure BDA0003750660340000112
表示为
Figure BDA0003750660340000113
其中K代表独立环路的个数,i-j∈f表示分支i-j与变压器fth相连;
dij=-dij.dij∈{-1,1} (3)
式中dij描述了分支的正方向,其在计算模型之前自定义,方向参数在优化前是恒定且已知的;
Figure BDA0003750660340000114
式中:LP,i和LQ,i为总线i处的有功负荷和无功负荷,Pij,i和Qij,i为支路i-j的终端i的有功电源和无功电源,
Figure BDA0003750660340000115
为功率因数;
上式描述了功率平衡的约束条件:节点注入的有功功率等于连接该节点的各支路的功率之和,无功功率由有功功率和相应的功率因数求得;
分布式发电机DG母线的功率输出约束:
Figure BDA0003750660340000121
式中:
Figure BDA0003750660340000122
Figure BDA0003750660340000123
分别为连接到总线i上的DG的最大有功和无功功率输出,N(i)为通过分支连接到总线i的总线,ΦDG为连接有分布式电源的总线集;
各支路的热额定值约束:
Figure BDA0003750660340000124
式中:
Figure BDA0003750660340000125
为分支i-j的功率限制;
在计算的过程中采用一次变量代替二次变量,如下式:
Figure BDA0003750660340000126
式中:
Figure BDA0003750660340000127
为母线i的电压的平方,
Figure BDA0003750660340000128
为分支i-j的电流的平方,连接支路两端节点间的电压关系表示为:
Figure BDA0003750660340000129
式中:Pij,i、Qij,i和Xij,i分别为分支i-j的终端i的有功功率、无功功率和电抗,Xij为分支i-j的电抗;
将上述式子中的节点电压限定在上下边界内:
Figure BDA00037506603400001210
式中:
Figure BDA00037506603400001211
Figure BDA00037506603400001212
分别为母线i的电压大小的平方的上边界和下边界,支路电流的表达式为:
Figure BDA0003750660340000131
无功损耗和有功损耗的表达式为:
Figure BDA0003750660340000132
变压器的额定功率约束为:
Figure BDA0003750660340000133
式中:
Figure BDA0003750660340000134
为变压器的额定值,如果变压器发生故障且故障在考虑的场景之一被隔离,则该变压器的额定值为0,式中
Figure BDA0003750660340000135
为电压的平方根;
上述式子中f的取值如式(13)所示,表示系统共有的几种情况,包含几种单个变压器故障的情况和一种正常情况,用f=0表示正常情况,即所有的变压器均无故障:
f=0,1,2,K Nt (13)
由于未断开支路的两端不应有式(8)中所描述的电压关系,所以式(8)不能用于未断开支路,采用下式代替:
Figure BDA0003750660340000136
当分支i-j连接时
Figure BDA0003750660340000137
等于1,
Figure BDA0003750660340000138
等于0,因此,式(14)等于式(8);当分支i-j断开时,
Figure BDA0003750660340000141
等于0,
Figure BDA0003750660340000142
等于M0;因此,
Figure BDA0003750660340000143
的取值范围为[-M0,M0];
为了解决式(10)中不同变量的乘积模型不可解的问题,通过SOCP松弛使用下式代替式(10):
Figure BDA0003750660340000144
Figure BDA0003750660340000145
其中式(15)为网络损耗最小化,SOCP松弛用来凸化TSC编程模型,修正后的模型表示为:
Figure BDA0003750660340000146
s,t,(2)-(7),(9),(11)-(14),(16)
根据相应的网络配置,利用上述约束条件求解功率流,对任何变压器故障而言,系统中的所有负载都在满足约束的情况下为系统提供功率;系统中的变量和约束条件的总集和具有相同的目标,即考虑网络所能提供总负荷的最大化,该模型考虑了潮流求解约束,通过网络重构来挖掘整个网络对系统负荷的供电能力。
所述的式(17)给出的模型采用蜉蝣算法求解实际配电网考虑重构的最大供电能力,具体的步骤如下:
第一步:随机生成雄性蜉蝣和雌性蜉蝣种群,初始化参数;
第二步:根据步骤C中的目标函数计算各个蜉蝣的适应度值并排序,得到个体最优和群体最优;
第三步:雄性蜉蝣的运动:
雄性蜉蝣一般成群结队活动,这就意味着每个雄虫的位置是根据自己和它的邻居的经验来调整的,假设
Figure BDA0003750660340000151
是第i个蜉蝣在t时刻搜索空间中的位置,蜉蝣的位置是通过在当前位置上增加一个速度
Figure BDA0003750660340000152
来改变的,因此,蜉蝣的第t+1时刻的位置表达式为:
Figure BDA0003750660340000153
考虑到雄性蜉蝣总是在离水面几米高的地方进行舞蹈,因此他们的速度并不会很高,所以,一只雄性蜉蝣i的速度计算为:
Figure BDA0003750660340000154
式中:
Figure BDA0003750660340000155
为第i个雄性蜉蝣在搜索空间j=1,…,n中t时刻的速度,
Figure BDA0003750660340000156
表示t时刻的位置,a1和a2分别为衡量规模认知和社会贡献的积极吸引常数,pbesti表示雄性蜉蝣i在搜索空间中的最好位置,其表达式为:
Figure BDA0003750660340000157
式中:F()为目标函数,同步骤C中的式子(1)所示;
对于算法性能来说,群体中的最优雄性蜉蝣继续执行其上下舞蹈,因此,最优的雄性蜉蝣不断更新他们的速度,其表达式为:
Figure BDA0003750660340000158
其中d为舞蹈系数,r为[-1,1]范围内的随机值;
第四步:雌性蜉蝣的运动:
与雄性蜉蝣不同,雌雄蜉蝣不会有聚集性活动,它们会飞向雄性进行繁殖,假设
Figure BDA0003750660340000159
是t时刻雌性蜉蝣在搜索空间中的位置,同时通过添加速度
Figure BDA00037506603400001510
来改变其位置,即:
Figure BDA0003750660340000161
考虑雄性蜉蝣和雌性蜉蝣相互吸引,因此,它们的速度计算公式为:
Figure BDA0003750660340000162
式中:
Figure BDA0003750660340000163
为第i个雌性蜉蝣在搜索空间j=1,…,n中t时刻的速度,
Figure BDA0003750660340000164
表示t时刻的位置;a2为相互之间的吸引常数,β是一个固定的可见性系数;fl是一个随机游走系数,当雌性不被雄性吸引时使用;r为[-1,1]之间的随机数;rmf为雌性和雄性蜉蝣之间的笛卡尔距离,其表达式为:
Figure BDA0003750660340000165
第五步:蜉蝣的交配:
每对蜉蝣交配后产生两个后代,表达式为:
offspring1=L*male+(1-L)*female (25)
offspring2=L*female+(1-L)*male (26)
式中male为父系,female为母系;L为特定范围内的随机数,后代的初始速度设定为0;
第六步:计算适应度值并更新个体最优和群体最优;
第七步:判断是否满足停止田间,若满足,则停止迭代,输出最优结果,否则转到第三步继续执行。
实施例:
测试系统如图1所示,变压器母线出线处包含一个断路器,其中每个馈线包含图2所示的标准的IEEE-33总线网络,假设各个节点均装有无功就地补偿装置且所有馈线的功率因数为0.9;节点电压的上下限分别为1.1p.u和0.9p.u,支路允许最大功率为9MVA。所有的负荷均为恒定的,因此在重构的过程中电压变化不会影响这些母线的有功和无功的需求。其中变压器A和C的容量均为31.5MW,C的变电容量为50MW。
根据步骤B构建图1所示的网络,其中包含两个变电站,变电站I中有一个变压器A,变电站II中有两个变压器B和C,其中变电站I的馈线AI、AII和AⅢ分别连接变电站II的馈线BI、BII和CI,且变电站II的馈线BII和CI相互连接。
根据步骤C建立测试网络的考虑配电网重构的最大供电能力计算的数学模型。
当变压器C发生故障时,考虑实际运行约束时,采用步骤D中的蜉蝣算法计算得到考虑配电网重构的最大供电能力计算值为32.92MVA,各条馈线所承担的负荷如表1所示。此时,线路CI面临失电的风险,可由线路AⅢ和BII转供其负荷。而不考虑配电网重构的供电能力计算值为29.34MW,明显低于本发明,配电网的供电能力并没有得到很好的挖掘。因此,本发明给出的考虑配电网重构的供电能力挖掘能够很好的挖掘配电网的供电能力。
表1:考虑配电网重构的不同馈线的最大供电能力
馈线 AI AII AⅢ BI BII CI
供电能力(MVA) 6.31 6.34 6.57 4.13 6.19 3.38

Claims (2)

1.一种基于配电网重构的供电能力挖掘方法,通过计算供电系统的最大供电能力来指导电网的规划,改变配电网中联络开关和分段开关的状态,引入分布式电源的方法,将配电网络的拓扑结构改变,馈线和变电站之间的平均负荷分配实现,提高负荷供应能力;同时,采用蜉蝣算法计算最大供电能力,拓宽问题的求解方法,具体包括以下几个步骤:
A、收集配电网的信息;
B、根据收集的配电网信息,构建具体的配电网网络;
C、根据配电网网络构建基于配电网重构的最大供电能力计算的数学模型:
第一步:建立目标函数:
Figure FDA0003750660330000011
式中:F为目标函数,Lp,i为线路i上的有功功率,ΦN代表所有的总线集和,N表示配电网络有N条线路;
目标是使系统中所有负荷的有功功率最大化,为了简化,设定负载的功率因数为常数(0.9),当确定了同一母线上的有功需求时,该线路自动获得所需的无功功率;
第二步:提供网络的约束条件:
Figure FDA0003750660330000012
式中:f代表变压器出线故障的场景,xi,j表示分支i-j的电抗,Nmode表示不包含根总线的所有总线数量,B表示与故障变压器直接连接的一组支路;
将分支的状态描述为Boolean值,对于径向拓扑,所有连接的线路数等于节点数减去线路条数;第二行的公式表示故障变压器支路i-j由于偶然事故而被断开;在分析网状网络的过程中,将
Figure FDA0003750660330000021
表示为
Figure FDA0003750660330000022
其中K代表独立环路的个数,i-j∈f表示分支i-j与变压器fth相连;
dij=-dij.dij∈{-1,1} (3)
式中dij描述了分支的正方向,其在计算模型之前自定义,方向参数在优化前是恒定且已知的;
Figure FDA0003750660330000023
式中:LP,i和LQ,i为总线i处的有功负荷和无功负荷,Pij,i和Qij,i为支路i-j的终端i的有功电源和无功电源,
Figure FDA0003750660330000027
为功率因数;
上式描述了功率平衡的约束条件:节点注入的有功功率等于连接该节点的各支路的功率之和,无功功率由有功功率和相应的功率因数求得;
分布式发电机DG母线的功率输出约束:
Figure FDA0003750660330000024
式中:
Figure FDA0003750660330000025
Figure FDA0003750660330000026
分别为连接到总线i上的DG的最大有功和无功功率输出,N(i)为通过分支连接到总线i的总线,ΦDG为连接有分布式电源的总线集;
各支路的热额定值约束:
Figure FDA0003750660330000031
式中:
Figure FDA0003750660330000032
为分支i-j的功率限制;
在计算的过程中采用一次变量代替二次变量,如下式:
Figure FDA0003750660330000033
式中:
Figure FDA0003750660330000034
为母线i的电压的平方,
Figure FDA0003750660330000035
为分支i-j的电流的平方,连接支路两端节点间的电压关系表示为:
Figure FDA0003750660330000036
式中:Pij,i、Qij,i和Xij,i分别为分支i-j的终端i的有功功率、无功功率和电抗,Xij为分支i-j的电抗;
将上述式子中的节点电压限定在上下边界内:
Figure FDA0003750660330000037
式中:
Figure FDA0003750660330000038
Figure FDA0003750660330000039
分别为母线i的电压大小的平方的上边界和下边界,支路电流的表达式为:
Figure FDA00037506603300000310
无功损耗和有功损耗的表达式为:
Figure FDA00037506603300000311
变压器的额定功率约束为:
Figure FDA0003750660330000041
式中:
Figure FDA0003750660330000042
为变压器的额定值,如果变压器发生故障且故障在考虑的场景之一被隔离,则该变压器的额定值为0,式中
Figure FDA0003750660330000043
为电压的平方根;
上述式子中f的取值如式(13)所示,表示系统共有的几种情况,包含几种单个变压器故障的情况和一种正常情况,用f=0表示正常情况,即所有的变压器均无故障:
f=0,1,2,K Nt (13)
由于未断开支路的两端不应有式(8)中所描述的电压关系,所以式(8)不能用于未断开支路,采用下式代替:
Figure FDA0003750660330000044
当分支i-j连接时
Figure FDA0003750660330000045
等于1,
Figure FDA0003750660330000046
等于0,因此,式(14)等于式(8);当分支i-j断开时,
Figure FDA0003750660330000047
等于0,
Figure FDA0003750660330000048
等于M0;因此,
Figure FDA0003750660330000049
的取值范围为[-M0,M0];
为了解决式(10)中不同变量的乘积模型不可解的问题,通过SOCP松弛使用下式代替式(10):
Figure FDA00037506603300000410
Figure FDA0003750660330000051
其中式(15)为网络损耗最小化,SOCP松弛用来凸化TSC编程模型,修正后的模型表示为:
Figure FDA0003750660330000052
根据相应的网络配置,利用上述约束条件求解功率流,对任何变压器故障而言,系统中的所有负载都在满足约束的情况下为系统提供功率;系统中的变量和约束条件的总集和具有相同的目标,即考虑网络所能提供总负荷的最大化,该模型考虑了潮流求解约束,通过网络重构来挖掘整个网络对系统负荷的供电能力。
2.根据权利要求1所述的基于配电网重构的供电能力挖掘方法,其特征在于:所述的式(17)给出的模型采用蜉蝣算法求解实际配电网考虑重构的最大供电能力,具体的步骤如下:
第一步:随机生成雄性蜉蝣和雌性蜉蝣种群,初始化参数;
第二步:根据步骤C中的目标函数计算各个蜉蝣的适应度值并排序,得到个体最优和群体最优;
第三步:雄性蜉蝣的运动:
雄性蜉蝣一般成群结队活动,这就意味着每个雄虫的位置是根据自己和它的邻居的经验来调整的,假设
Figure FDA0003750660330000053
是第i个蜉蝣在t时刻搜索空间中的位置,蜉蝣的位置是通过在当前位置上增加一个速度
Figure FDA0003750660330000054
来改变的,因此,蜉蝣的第t+1时刻的位置表达式为:
Figure FDA0003750660330000061
考虑到雄性蜉蝣总是在离水面几米高的地方进行舞蹈,因此他们的速度并不会很高,所以,一只雄性蜉蝣i的速度计算为:
Figure FDA0003750660330000062
式中:
Figure FDA0003750660330000063
为第i个雄性蜉蝣在搜索空间j=1,…,n中t时刻的速度,
Figure FDA0003750660330000064
表示t时刻的位置,a1和a2分别为衡量规模认知和社会贡献的积极吸引常数,pbesti表示雄性蜉蝣i在搜索空间中的最好位置,其表达式为:
Figure FDA0003750660330000065
式中:F()为目标函数,同步骤C中的式子(1)所示;
对于算法性能来说,群体中的最优雄性蜉蝣继续执行其上下舞蹈,因此,最优的雄性蜉蝣不断更新他们的速度,其表达式为:
Figure FDA0003750660330000066
其中d为舞蹈系数,r为[-1,1]范围内的随机值;
第四步:雌性蜉蝣的运动:
与雄性蜉蝣不同,雌雄蜉蝣不会有聚集性活动,它们会飞向雄性进行繁殖,假设
Figure FDA0003750660330000067
是t时刻雌性蜉蝣在搜索空间中的位置,同时通过添加速度
Figure FDA0003750660330000068
来改变其位置,即:
Figure FDA0003750660330000069
考虑雄性蜉蝣和雌性蜉蝣相互吸引,因此,它们的速度计算公式为:
Figure FDA0003750660330000071
式中:
Figure FDA0003750660330000072
为第i个雌性蜉蝣在搜索空间j=1,…,n中t时刻的速度,
Figure FDA0003750660330000073
表示t时刻的位置;a2为相互之间的吸引常数,β是一个固定的可见性系数;fl是一个随机游走系数,当雌性不被雄性吸引时使用;r为[-1,1]之间的随机数;rmf为雌性和雄性蜉蝣之间的笛卡尔距离,其表达式为:
Figure FDA0003750660330000074
第五步:蜉蝣的交配:
每对蜉蝣交配后产生两个后代,表达式为:
offspring1=L*male+(1-L)*female (25)
offspring2=L*female+(1-L)*male (26)
式中male为父系,female为母系;L为特定范围内的随机数,后代的初始速度设定为0;
第六步:计算适应度值并更新个体最优和群体最优;
第七步:判断是否满足停止田间,若满足,则停止迭代,输出最优结果,否则转到第三步继续执行。
CN202210840243.0A 2022-07-18 2022-07-18 一种基于配电网重构的供电能力挖掘方法 Pending CN115423242A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210840243.0A CN115423242A (zh) 2022-07-18 2022-07-18 一种基于配电网重构的供电能力挖掘方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210840243.0A CN115423242A (zh) 2022-07-18 2022-07-18 一种基于配电网重构的供电能力挖掘方法

Publications (1)

Publication Number Publication Date
CN115423242A true CN115423242A (zh) 2022-12-02

Family

ID=84196965

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210840243.0A Pending CN115423242A (zh) 2022-07-18 2022-07-18 一种基于配电网重构的供电能力挖掘方法

Country Status (1)

Country Link
CN (1) CN115423242A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116094169A (zh) * 2023-01-28 2023-05-09 国网江苏省电力有限公司连云港供电分公司 一种配电网拓扑模型生成方法及终端设备

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116094169A (zh) * 2023-01-28 2023-05-09 国网江苏省电力有限公司连云港供电分公司 一种配电网拓扑模型生成方法及终端设备
CN116094169B (zh) * 2023-01-28 2024-04-12 国网江苏省电力有限公司连云港供电分公司 一种配电网拓扑模型生成方法及终端设备

Similar Documents

Publication Publication Date Title
Syahputra Fuzzy multi-objective approach for the improvement of distribution network efficiency by considering DG
Aghaei et al. Distribution expansion planning considering reliability and security of energy using modified PSO (Particle Swarm Optimization) algorithm
Syahputra et al. Reconfiguration of distribution network with DG using fuzzy multi-objective method
Syahputra et al. Distribution network efficiency improvement based on fuzzy multi-objective method
CN104362623A (zh) 一种适用于主动配电网的多目标网络重构方法
Fu et al. A short-term preventive maintenance scheduling method for distribution networks with distributed generators and batteries
CN108462210B (zh) 基于数据挖掘的光伏可开放容量计算方法
CN104037765A (zh) 基于改进遗传算法选取有源配电网供电恢复方案的方法
Parizad et al. Optimal distribution systems reconfiguration for short circuit level reduction using PSO algorithm
CN112434905A (zh) 考虑多次转供对可靠性影响的配电系统供电能力评估方法
CN115423242A (zh) 一种基于配电网重构的供电能力挖掘方法
CN109004639B (zh) 基于完全分布式算法的配电网分区供电恢复策略优化方法
CN114285085A (zh) 一种基于网络重构的电网弹性资源广域调度方法
He et al. Multi-objective operation mode optimization of medium voltage distribution networks based on improved binary particle swarm optimization
Arias Planning models for distribution grid: A brief review
Ping et al. Research on comprehensive optimization of power intensive uhv/ehv sending power grid
Wang et al. A Novel Method for Islanding in Active Distribution Network Considering Distributed Generation.
Zheng et al. Collaborative Planning of Primary and Secondary Equipment in Distribution Network
Wang et al. Distribution Network Fault Location and Recovery Considering Load Importance
Wang et al. Fault Recovery Strategy of Distribution Network with Distributed Generation Based on Island Division
Wu Multi-objective distribution network reconstruction based on decimal coding
Liu et al. Non-cooperative game theory based active distribution network recovery strategy
Bai Reconfiguration performance of the urban power distribution system based on the genetic-ant colony fusion algorithm
Lu et al. Research on Power Grid Optimization Partition Based on Improved MOPSO Algorithm
Jianlin et al. Research on the optimization strategy of 220kV power grid based on the partitioning connecting channel

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