CN112200410A - 一种基于离散猴群算法的配变群集负载均衡优化方法 - Google Patents
一种基于离散猴群算法的配变群集负载均衡优化方法 Download PDFInfo
- Publication number
- CN112200410A CN112200410A CN202010942314.9A CN202010942314A CN112200410A CN 112200410 A CN112200410 A CN 112200410A CN 202010942314 A CN202010942314 A CN 202010942314A CN 112200410 A CN112200410 A CN 112200410A
- Authority
- CN
- China
- Prior art keywords
- cluster
- feeder line
- distribution transformer
- monkey
- load
- 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
- 238000000034 method Methods 0.000 title claims abstract description 69
- 241000282693 Cercopithecidae Species 0.000 title claims abstract description 46
- 238000005457 optimization Methods 0.000 title claims abstract description 39
- 230000008569 process Effects 0.000 claims description 44
- 230000006870 function Effects 0.000 claims description 13
- 230000009194 climbing Effects 0.000 claims description 10
- 230000009193 crawling Effects 0.000 claims description 9
- 230000008859 change Effects 0.000 claims description 7
- 230000009466 transformation Effects 0.000 claims description 7
- 238000012546 transfer Methods 0.000 claims description 6
- 230000008878 coupling Effects 0.000 claims description 3
- 238000010168 coupling process Methods 0.000 claims description 3
- 238000005859 coupling reaction Methods 0.000 claims description 3
- 230000009471 action Effects 0.000 description 6
- 239000011159 matrix material Substances 0.000 description 6
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 4
- RYGMFSIKBFXOCR-UHFFFAOYSA-N Copper Chemical compound [Cu] RYGMFSIKBFXOCR-UHFFFAOYSA-N 0.000 description 3
- 229910052802 copper Inorganic materials 0.000 description 3
- 239000010949 copper Substances 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 239000000243 solution Substances 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000004134 energy conservation Methods 0.000 description 2
- 229910052742 iron Inorganic materials 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 239000007983 Tris buffer Substances 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000009191 jumping Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000002028 premature Effects 0.000 description 1
- 230000005855 radiation Effects 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/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
- G06Q10/06313—Resource planning in a project environment
-
- 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
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/27—Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/004—Artificial life, i.e. computing arrangements simulating life
- G06N3/006—Artificial 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]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- 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]
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Evolutionary Computation (AREA)
- Strategic Management (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Health & Medical Sciences (AREA)
- Entrepreneurship & Innovation (AREA)
- Marketing (AREA)
- General Health & Medical Sciences (AREA)
- Software Systems (AREA)
- Tourism & Hospitality (AREA)
- Computer Hardware Design (AREA)
- General Business, Economics & Management (AREA)
- Life Sciences & Earth Sciences (AREA)
- Artificial Intelligence (AREA)
- Biophysics (AREA)
- Primary Health Care (AREA)
- Computational Linguistics (AREA)
- Data Mining & Analysis (AREA)
- Quality & Reliability (AREA)
- Operations Research (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- Game Theory and Decision Science (AREA)
- Mathematical Physics (AREA)
- Educational Administration (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- Biomedical Technology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Medical Informatics (AREA)
- Development Economics (AREA)
- Biodiversity & Conservation Biology (AREA)
- Computer Networks & Wireless Communication (AREA)
Abstract
本发明公开了一种基于离散猴群算法的配变群集负载均衡优化方法,包括如下步骤:S1、根据配电网网络参数利用联络开关以及分段开关的开合状态,求解变电站群集中每台配变自身馈线上所供负荷;S2、求解通过联络线路转移到每台配变馈线上的负荷;S3、求解配变群集的有功损耗;S4、建立以开关设备的开合状态为决策变量的配变群集负载时空优化模型;S5、利用离散猴群算法求解优化模型,该方法在避免开关设备频繁操作的前提下,可以有效降低变电站群集中配变的有功损耗。
Description
技术领域
本发明适用于含有多台变电站的配电网节能降损领域,尤其涉及一种基于离散猴群算法的配变群集负载均衡优化方法。
背景技术
目前配电网的调控还不够灵活,配变的各条出线向固定的负荷点供电,没有很好地实现配变群集的负载时空均衡优化。通过优化配电网中分段开关以及联络开关的闭合状态,可以有效降低变电站群集的总损耗,能够使配变群集在满足供电可靠性的前提下负载率更均衡,进而实现节能减排。
配变之间往往通过联络线路等设备存在联络关系,通过优化配电网中分段开关以及联络开关的闭合状态,对配电网进行正常态的网络重构,可以使得变电站的负载率更均衡,进而有效降低变电站群集的总损耗。有关于配变群集的负载均衡研究还存在以下两点不足:(1) 开关设备的开合操作会影响供电可靠性和设备的使用寿命,已有研究没能有效计及操作频率约束;(2)将变压器等效为一段含阻抗的线路,关注配变下游的馈线网络重构,对变电站群集的负载均衡优化关注不够。此外,配网重构是一个非线性组合优化问题,开关状态组合数量多,求解难度大,一般需要借助智能算法求解,遗传算法及其改进算法因操作的简单性而得到了广泛应用,但仍存在着早熟收敛和收敛性能差的缺点。
发明内容
本发明将提出一种基于离散猴群算法的配变群集负载均衡优化方法,该方法以一段时间内而不是某一时间断面上的网损成本最小为目标函数,在满足配电系统“N-1”安全约束的前提下,完成负载时空均衡优化。其难点在于如何建立以开关设备的开合状态为决策变量的配变群集负载时空优化模型,并快速求解多变量优化模型。为此,本发明建立了以网损和开关设备动作成本最小为目标函数的优化模型,计及了非孤岛运行、开关操作频率、N-1安全等约束条件,并对传统的离散猴群算法加以改进,从而实现对优化模型的快速求解。
本发明通过如下技术方案实现:
一种基于离散猴群算法的配变群集负载均衡优化方法,包括如下步骤:
S1、根据配电网网络参数,以联络开关以及分段开关的开合状态为变量,求解变电站群集中每台配变自身馈线上所供负荷;
S2、求解通过联络线路转移到每台配变馈线上的负荷;
S3、求解配变群集的有功损耗;
S4、建立以开关设备的开合状态为决策变量的配变群集负载时空优化模型;
S5、利用离散猴群算法求解优化模型。
所述步骤S2中求解通过联络线路转移到每台配变馈线上的负荷过程:
2.1、根据联络开关状态更新每条馈线指示变量SIi,j,l(t);
2.2、根据配网集群的联络情况,通过下式计算Pl,t,tr;
其中:Pl,t,tr表示其他馈线经第l个联络开关向第j条馈线转移的负荷量。
2.3、根据更新每条馈线指示变量SIi,j,l(t)和Pl,t,tr通过如下公式计算Pi,t,tr;
其中,Pi,t,tr表示通过联络线路转供到第i台配变Ni,f条馈线的负载,第i台配变第j条馈线的联络开关数量记为Ni,j,tie,联络开关编号为1,2,…,Ni,j,tie;第j条馈线Ni,j,oth段馈线中有且仅有一段馈线与第l个联络开关相连,与第l个联络开关相连的馈线段编号记为g(l),该对应关系由网络拓扑决定。
所述S5步骤中利用离散猴群算法求解优化模型过程:
5.1、初始化:考虑到决策变量存在时间维度上的耦合关系,令每条馈线上联络开关处于打开状态,其他开关均处于关闭状态,然后,对某一条馈线上的分段开关进行随机变异,以此为第一只猴子的初始位置X1。以此类推,直至产生种群内所有猴子的初始位置;
5.2、爬过程
爬过程具体步骤如下:
对于第i只猴子而言,爬过程具体步骤如下:
步骤1:判断该猴子馈线中是否存在负荷转供,如果有,转至下一步骤,如果没有,爬过程结束。
步骤2:对第j条馈线中发生开闭状态变化的常闭开关的打开时刻进行寻优,令:
如果满足:
则把Xi更新为Xi'。如果满足:
则把Xi更新为Xi”。重复此步骤,直到目标函数值不发生变化。
步骤3.对第j条馈线发生开闭状态变化的常闭开关位置进行寻优,对Xi中非零列的x1,j,在[0,N1,j]中逐个遍历,并重复步骤2,修正Xi,爬过程结束。
对于第i只猴子,该过程步骤如下:
步骤2:如果G(Bi)<G(Xi),则令Xi=Bi,调用爬过程;
步骤3:重复步骤1和步骤2,直至满足望-跳过程的终止条件(比如达到预先设定的望跳次数)。
5.4、翻过程
翻过程具体步骤如下:
步骤1.计算猴群所处的平均位置:
步骤2.在[-10Ni,j,oth,10Ni,j,oth]内产生实数β,对Xi如下变化:
步骤3.判断Xi是否满足约束,如果满足,转至步骤6,如果不满足,转至步骤4;
步骤4:创建与Xi相同维度的零矩阵Z,Z的列元素依次替换为Xi中的元素,每次替换均判断是否满足式(7)。对于第k次(k≤Nto)替换,如果满足约束,则继续第k+1次替换;如果不满足,表明前k-1次替换中,有且仅有一次非零列的替换与第k次替换存在冲突。存在冲突的两列一定存在联络关系,根据联络矩阵T确定该非零列,这两个非零列随机选取一个置零,并继续上述替换过程,直至完成所有元素的替换。
步骤5:重复步骤2-步骤4,直至完成对所有猴子的修正;
步骤6:对新生成的猴群,调用爬过程和望-跳过程,直至达到收敛终止条件。
有益效果
本发明基于配电网中分段开关和联络开关时空维度的约束,该方法避免了开关频繁操作影响部分负荷用电的同时,可以有效降低变电站群集中配变的有功损耗,并且通过改进的猴群算法对非线性整数规划模型进行求解,算例分析结果表明本发明所提出的方法能够满足运行约束,具有较好的寻优能力。
附图说明
图1是本发明基于改进猴群算法的变电站群集负载时空均衡优化方法制定流程图;
图2是网络的拓扑图;
图3是不同迭代次数下目标函数的变化过程。
具体实施方式:
以下结合实施案例和附图,对本发明的技术和方法进行详细描述,以下案例和附图用于说明本发明的构成,但不是用来限定本发明的范围。
本发明的目的在于提出一种基于离散猴群算法的配变群集负载均衡优化方法,通过优化配电网中分段开关以及联络开关的闭合状态,对配电网进行正常态的网络重构,可以使得变电站的负载率更均衡,进而有效降低变电站群集的总损耗。本发明:①以一段时间内而不是某一时间断面上的网损成本最小为目标函数,在满足配电系统“N-1”安全约束的前提下,完成负载均衡优化;②在时间维度上计及了开关操作频率约束,在空间维度上计及了潮流约束、电压约束、传输容量约束、配电网辐射运行约束;③由于优化模型中决策变量维度高,模型求解存在困难,对离散猴群算法加以改进,用于配变群集的负载均衡研究。具体而言,本发明将通过以下五个步骤通过网络重构降低变电站群集的总损耗,通过如下步骤实现:
如图1所示,1.求解变电站群集中每台配变自身馈线上所供负荷
辐射状运行的配电网中每段馈线的负荷有且仅有一条供电路径,当且仅当该路径上每个开关均闭合时,这段馈线上的负荷由第i台配变供电,否则该段馈线上的负荷通过联络开关转供到了其他配变上(正常态不考虑缺供电)。第i台配变第j条馈线的常闭开关(包括断路器、分段开关等)的数量记为Ni,j,oth,这些开关编号为1,2,…,Ni,j,oth;第j条馈线上存在 Ni,j,oth段馈线,每段馈线可以向负荷供电,令这段馈线上的负荷编号与这段馈线首端的常闭开关编号一致;第k段馈线供电路径上的开关集合记为Ωi,j,k,k∈[1,Ni,j,oth],该集合是由馈线拓扑结构决定的,集合中的开关开合状态记为SWm(t),m∈Ωi,j,k,1表示闭合,0表示打开。则存在:
其中,Pi,t,or表示第i台配变母线Ni,f条出线的负载,PLi,j,k,t表示第i台配变第j条馈线第k段馈线t时刻的视在功率。
2.求解通过联络线路转移到每台配变馈线上的负荷
步骤一:分析第i台配变第j条馈线是否具备转供条件。
配电系统在正常态一般不进行二次转供,即当且仅当第i台配变第j条馈线向本馈线所有负荷供电时,相邻馈线才可能会向该馈线转移负荷,定义指示变量SIi,j,l(t):
其中,sign(·)是符号函数,当函数变量大于0时函数值取1,当函数变量等于0时函数值取0,否则函数值取-1。联络开关的开合状态记为SWTi,j,l(t),1表示闭合,0表示打开,SIi,j,l(t)等于 1时,标志着第i台配变第j条馈线通过第l个联络开关向相邻馈线的负荷供电。
步骤二:根据配网集群的联络情况计算Pl,t,tr
Pl,t,tr表示其他馈线经第l个联络开关向第j条馈线转移的负荷量。配电系统馈线总数量为则存在阶方阵T,T中的每个要素表示行号和列号所对应的馈线的联络情况,存在联络记为1,否则记为0,对角线元素取0,该矩阵由配电系统的拓扑结构唯一确定。当 SIi,j,l(t)=0时,Pl,t,tr为零;当SIi,j,l(t)=1时,根据T确定与第i台配变第j条馈线通过第l个联络开关存在联络关系的线路为第i′台配变的第j′条馈线,此时可得到:
步骤三:根据步骤一和步骤二的所得结果计算Pi,t,tr
Pi,t,tr表示通过联络线路转供到第i台配变Ni,f条馈线的负载,第i台配变第j条馈线的联络开关数量记为Ni,j,tie,联络开关编号为1,2,…,Ni,j,tie。第j条馈线Ni,j,oth段馈线中有且仅有一段馈线与第l个联络开关相连,与第l个联络开关相连的馈线段编号记为g(l),该对应关系由网络拓扑决定。因此存在:
3.求解配变群集的有功损耗
变压器有功损耗包括铁耗和铜耗,其中铁耗同电压水平相关,铜耗同负载率相关,本专利主要关注配变的铜耗。含有Ns台变压器的配变群集,其有功损耗如下:
其中,ΔPi,t表示第i台配变t时刻的有功损耗;表示第i台配变的功率因数,一般由无功补偿设备调整到0.95以上;Ri表示第i台配变的电阻;Ui,t表示第i台配变t时刻的电压;Pi,t表示第i台配变t时刻的视在功率,包括该配变母线Ni,f条出线的负载Pi,t,or和通过联络线路转供到Ni,f条馈线的负载Pi,t,tr。
4.建立以开关设备的开合状态为决策变量的配变群集负载时空优化模型
以优化周期内配变群集的有功网损和开关设备的动作成本最小为目标函数:
其中,CkWh表示单位电能的成本,Nch表示优化周期T内所有开关设备的动作总次数,Cch表示开关设备单次动作成本。为了保证配电系统的安全运行,需要满足以下空间和时间维度的约束条件:
1)网络拓扑无孤岛约束
在改变开关设备开闭状态的过程中,不允许出现某些馈线段孤岛(即不与任何配变相连) 的情况,该约束可通过如下等式表示:
2)馈线不过载约束
开关设备改变开合状态后,不允许出现负荷转供导致的馈线过载情况:
其中,Si,j,rate表示第i台配变第j条馈线的额定容量。
3)辐射状运行约束
配电网正常态的网络重构,需要满足重构后的系统仍然保持辐射状的结构。在满足式(7) 的条件下,在任意时刻需要满足:
该式表示在网络重构过程中,由打开到闭合的开关数量与由闭合到打开的开关数量相等。
4)开关设备动作频率约束
开关设备开合状态的改变对设备寿命和系统供电可靠性都存在影响,工程上一般要求其动作频率不能太高:
式中⊕是异或数学运算符,Ni,j,l和NT分别表示优化周期T内一个联络开关的最大允许动作次数和所有联络开关的总动作次数。
5)配变群集“N-1”安全约束
第i台配变的负载率为:
其中,Si,rate表示第i台配变的额定容量。根据电力系统“N-1”安全约束,配变群集中的每台变压器存在负载率上限,文献[1]考虑主变互联关系,提出了配变群集中每台主变的最大负载率计算方法,故存在如下约束:
其中,Ti(Ns-1)表示第i台配变满足“N-1”校验的负载率上限,具体计算方法参考文献[1],此处不再展开。
5.利用离散猴群算法求解优化模型
本发明所建立模型是0-1整数非线性规划问题,时间维度的约束条件式(10)使得配变群集的运行呈现为非马尔科夫过程,求解难度增大,因此本专利对离散猴群算法进行改进,实现模型的快速求解。
在算法执行过程中随机产生的开关组合往往难以满足所优化模型的约束条件,因此对目标函数进行改进:
其中,Z表示随机产生的开关变量不满足约束条件时取一个极大值。
尽管猴群算法具有较好的寻优能力和收敛速度,但是由于爬过程用“伪梯度”法寻优,而配变群集的优化模型中决策变量为离散变量,不能直接使用,需加以改进。设猴群的规模为M,对猴群中的第i只猴子,算法执行过程中,第i只猴子其所在位置可定义为:
其中,n为变电站群集中的馈线总数,N1表示变电站群集中馈线上分段开关的数目,x1,j表示需要打开第j条馈线上第x1,j个分段开关,x2,n表示该分段开关的变化时刻,当馈线中某个分段开关发生改变的时候,线路中的联络开关状态也可以确定,这样每只猴子的当前位置即为所优化模型一个优化周期内的一个解。
(1)初始化
考虑到决策变量存在时间维度上的耦合关系,令每条馈线上联络开关处于打开状态,其他开关均处于关闭状态,然后,对某一条馈线上的分段开关进行随机变异,以此为第一只猴子的初始位置X1,以此类推,直至产生种群内所有猴子的初始位置。
(2)爬过程
对于第i只猴子而言,爬过程具体步骤如下:
步骤1.判断该猴子馈线中是否存在负荷转供,如果有,转至下一步骤,如果没有,爬过程结束。
步骤2.对第j条馈线中发生开闭状态变化的常闭开关的打开时刻进行寻优,令:
如果满足:
则把Xi更新为Xi'。如果满足:
则把Xi更新为Xi”。重复此步骤,直到目标函数值不发生变化。
步骤3.对第j条馈线发生开闭状态变化的常闭开关位置进行寻优,对Xi中非零列的x1,j,在[1,N1,j]中逐个遍历,并重复步骤2,修正Xi,爬过程结束。
(3)望-跳过程
对于第i只猴子,该过程步骤如下:
步骤2:如果G(Bi)<G(Xi),则令Xi=Bi,调用爬过程;
步骤3:重复步骤1和步骤2,直至满足望-跳过程的终止条件(比如达到预先设定的望- 跳次数)。
(4)翻过程
翻过程具体步骤如下:
步骤1.计算猴群所处的平均位置:
步骤2.在[-10Ni,j,oth,10Ni,j,oth]内产生实数β,对Xi如下变化:
步骤3.判断Xi是否满足式(7)约束,如果满足,转至步骤6,如果不满足,转至步骤4;
步骤4:创建与Xi相同维度的零矩阵Z,Z的列元素依次替换为Xi中的元素,每次替换均判断是否满足约束。对于第k次(k≤Nto)替换,如果满足约束,则继续第k+1次替换;如果不满足,表明前k-1次替换中,有且仅有一次非零列的替换与第k次替换存在冲突。存在冲突的两列一定存在联络关系,根据联络矩阵T确定该非零列,这两个非零列随机选取一个置零,并继续上述替换过程,直至完成所有元素的替换。
步骤5:重复步骤2-步骤4,直至完成对所有猴子的修正;
步骤6:对新生成的猴群,调用爬过程和望-跳过程,直至达到收敛终止条件。
(5)扰动过程
在翻过程的步骤2中,当所有猴子的某列元素均相同时,翻过程将失去效用。为此,需要在翻过程中引入随机扰动:随机选取Xk,令Xk的该列元素为零,之后再进行翻过程。随机扰动机制的引入增加了猴群的多样性,降低了搜索过程陷入局部最优的风险。
假设某一配电网中存在3座变电站,每座变电站含有2台110/10kV的主变,容量为50MW, 根据主变参数计算可得RT=0.219Ω,XT=3.111Ω,每条馈线线路传输容量为20MW,其拓扑图如图2所示,算法的收敛效果如图3所示,算法经过5次迭代之后趋于收敛,最优解目标函数值为2964.4元,相比不采取负载均衡时的网损成本3581.2元,网损降低了17.2%。
本发明并不限于上文描述的实施方式。以上对具体实施方式的描述旨在描述和说明本发明的技术方案,上述的具体实施方式仅仅是示意性的,并不是限制性的。在不脱离本发明宗旨和权利要求所保护的范围情况下,本领域的普通技术人员在本发明的启示下还可做出很多形式的具体变换,这些均属于本发明的保护范围之内。
Claims (3)
1.一种基于离散猴群算法的配变群集负载均衡优化方法,包括如下步骤:
S1、根据配电网网络参数,以联络开关以及分段开关的开合状态为变量,求解变电站群集中每台配变自身馈线上所供负荷;
S2、求解通过联络线路转移到每台配变馈线上的负荷;
S3、求解配变群集的有功损耗;
S4、建立以开关设备的开合状态为决策变量的配变群集负载时空优化模型;
S5、利用离散猴群算法求解优化模型。
2.根据权利要求1所述的一种面向变电站群集的负载时空均衡优化方法,其特征在于,所述步骤S2中求解通过联络线路转移到每台配变馈线上的负荷:
2.1、根据联络开关状态更新每条馈线指示变量SIi,j,l(t);
2.2、根据配网集群的联络情况,通过下式计算Pl,t,tr;
其中:Pl,t,tr表示其他馈线经第l个联络开关向第j条馈线转移的负荷量。
2.3、根据更新后每条馈线指示变量SIi,j,l(t)和Pl,t,tr通过如下公式计算Pi,t,tr;
其中,Pi,t,tr表示通过联络线路转供到第i台配变Ni,f条馈线的负载,第i台配变第j条馈线的联络开关数量记为Ni,j,tie,联络开关编号为1,2,…,Ni,j,tie;第j条馈线Ni,j,oth段馈线中有且仅有一段馈线与第l个联络开关相连,与第l个联络开关相连的馈线段编号记为g(l),该对应关系由网络拓扑决定。
3.根据权利要求1所述的一种基于离散猴群算法的配变群集负载均衡优化方法,其特征在于,所述S5步骤中利用离散猴群算法求解优化模型过程:
5.1、初始化:考虑到决策变量存在时间维度上的耦合关系,令每条馈线上联络开关处于打开状态,其他开关均处于关闭状态,然后,对某一条馈线上的分段开关进行随机变异,以此为第一只猴子的初始位置X1。以此类推,直至产生种群内所有猴子的初始位置。
5.2、爬过程
对于第i只猴子而言,爬过程具体步骤如下:
步骤1:判断该猴子馈线中是否存在负荷转供,如果有,转至下一步骤,如果没有,爬过程结束;
步骤2:对第j条馈线中发生开闭状态变化的常闭开关的打开时刻进行寻优,令:
如果满足:
则把Xi更新为Xi',如果满足:
则把Xi更新为Xi”,重复此步骤,直到目标函数值不发生变化;
步骤3:对第j条馈线发生开闭状态变化的常闭开关位置进行寻优,对Xi中非零列的x1,j,在[0,N1,j]中逐个遍历,并重复步骤2,修正Xi,爬过程结束。
5.3、望-跳过程
对于第i只猴子,该过程步骤如下:
步骤2:如果G(Bi)<G(Xi),则令Xi=Bi,调用爬过程;
步骤3:重复步骤1和步骤2,直至满足望-跳过程的终止条件(比如达到预先设定的望跳次数);
5.4、翻过程
翻过程具体步骤如下:
步骤1:计算猴群所处的平均位置:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010942314.9A CN112200410A (zh) | 2020-09-09 | 2020-09-09 | 一种基于离散猴群算法的配变群集负载均衡优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010942314.9A CN112200410A (zh) | 2020-09-09 | 2020-09-09 | 一种基于离散猴群算法的配变群集负载均衡优化方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN112200410A true CN112200410A (zh) | 2021-01-08 |
Family
ID=74014462
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010942314.9A Pending CN112200410A (zh) | 2020-09-09 | 2020-09-09 | 一种基于离散猴群算法的配变群集负载均衡优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112200410A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112798895A (zh) * | 2021-02-09 | 2021-05-14 | 南京软核科技有限公司 | 一种联络开关位置识别方法 |
CN117913818A (zh) * | 2024-01-24 | 2024-04-19 | 广东电网有限责任公司广州供电局 | 一种配电网馈线组互济运行优化方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100084211A1 (en) * | 2008-10-08 | 2010-04-08 | Richard Bruce Seidel | Hybrid-powered slow-speed vehicle, vehicle travel-range extension method, and method of hybrid-vehicle manufacture |
CN101814735A (zh) * | 2010-05-19 | 2010-08-25 | 天津大学 | 一种输电网扩展规划方法 |
CN103489047A (zh) * | 2013-09-30 | 2014-01-01 | 国家电网公司 | 一种中压配电网络分段开关的优化配置方法 |
CN105205740A (zh) * | 2015-10-23 | 2015-12-30 | 国家电网公司 | 考虑主变互联关系有效扩展的配电网供电能力评估方法 |
CN105207213A (zh) * | 2015-10-23 | 2015-12-30 | 国家电网公司 | 一种考虑分段负荷转移的配电网n-1校验方法 |
CN106779197A (zh) * | 2016-12-05 | 2017-05-31 | 河海大学 | 基于memetic算法的配电网开关优化配置研究方法 |
CN110061500A (zh) * | 2019-06-04 | 2019-07-26 | 广东电网有限责任公司 | 一种配电网供电能力的快速评估方法及系统 |
CN111564845A (zh) * | 2020-04-28 | 2020-08-21 | 国网福建省电力有限公司 | 一种考虑负荷特性和馈线分段的配电网供电能力评估方法 |
-
2020
- 2020-09-09 CN CN202010942314.9A patent/CN112200410A/zh active Pending
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100084211A1 (en) * | 2008-10-08 | 2010-04-08 | Richard Bruce Seidel | Hybrid-powered slow-speed vehicle, vehicle travel-range extension method, and method of hybrid-vehicle manufacture |
CN101814735A (zh) * | 2010-05-19 | 2010-08-25 | 天津大学 | 一种输电网扩展规划方法 |
CN103489047A (zh) * | 2013-09-30 | 2014-01-01 | 国家电网公司 | 一种中压配电网络分段开关的优化配置方法 |
CN105205740A (zh) * | 2015-10-23 | 2015-12-30 | 国家电网公司 | 考虑主变互联关系有效扩展的配电网供电能力评估方法 |
CN105207213A (zh) * | 2015-10-23 | 2015-12-30 | 国家电网公司 | 一种考虑分段负荷转移的配电网n-1校验方法 |
CN106779197A (zh) * | 2016-12-05 | 2017-05-31 | 河海大学 | 基于memetic算法的配电网开关优化配置研究方法 |
CN110061500A (zh) * | 2019-06-04 | 2019-07-26 | 广东电网有限责任公司 | 一种配电网供电能力的快速评估方法及系统 |
CN111564845A (zh) * | 2020-04-28 | 2020-08-21 | 国网福建省电力有限公司 | 一种考虑负荷特性和馈线分段的配电网供电能力评估方法 |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112798895A (zh) * | 2021-02-09 | 2021-05-14 | 南京软核科技有限公司 | 一种联络开关位置识别方法 |
CN117913818A (zh) * | 2024-01-24 | 2024-04-19 | 广东电网有限责任公司广州供电局 | 一种配电网馈线组互济运行优化方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
WO2022016622A1 (zh) | 真双极柔性直流输电系统故障下的自适应优化调控方法 | |
Niknam | An efficient hybrid evolutionary algorithm based on PSO and HBMO algorithms for multi-objective distribution feeder reconfiguration | |
CN103117545B (zh) | 一种智能变电站负荷自动分配方法 | |
CN107979092A (zh) | 一种考虑分布式电源与软开关接入的配电网动态重构方法 | |
Sedighizadeh et al. | Reconfiguration of distribution systems to improve re-liability and reduce power losses using imperialist com-petitive algorithm | |
CN113516278B (zh) | 有源配电网多时间尺度有功无功协调优化调度方法及系统 | |
CN112200410A (zh) | 一种基于离散猴群算法的配变群集负载均衡优化方法 | |
CN111490542B (zh) | 一种多端柔性多状态开关的选址定容方法 | |
CN108808738B (zh) | 一种考虑约束优先级的电网安全校正控制方法 | |
CN114362267B (zh) | 考虑多目标优化的交直流混合配电网分散式协调优化方法 | |
CN112217214B (zh) | 一种多级高压配电网协调转供优化方法 | |
CN109768546A (zh) | 基于多智能软开关协调的有源配电网供电恢复方法 | |
Hosseinpour et al. | Optimal placement of on-load tap changers in distribution networks using SA-TLBO method | |
CN114665479A (zh) | 一种考虑网络重构的配电网供电恢复方法及系统 | |
CN107528321B (zh) | 互联电网交直流并联输电通道功率随机优化分配方法 | |
CN112467748A (zh) | 三相不平衡主动配电网双时标分布式电压控制方法及系统 | |
CN116345466A (zh) | 一种考虑配网重构的主动配电网两阶段潮流优化方法 | |
CN116247737A (zh) | 一种基于有源配电网规划动态负荷裕度的分析方法 | |
CN114069664B (zh) | 一种规模化储能系统的配电网电压分布式控制方法 | |
CN115423242A (zh) | 一种基于配电网重构的供电能力挖掘方法 | |
CN111242389A (zh) | 一种智能储能软开关规划方法、系统、设备及介质 | |
Abd El-Hamed et al. | Self-healing restoration of a distribution system using hybrid Fuzzy Control/Ant-Colony Optimization Algorithm | |
CN105207226A (zh) | 一种高压电网无功补偿优化配置方法 | |
CN111091307A (zh) | 配电网资源处理方法、存储介质以及处理器 | |
CN115765035A (zh) | 一种适应于全时段动态重构的柔性配电网扰动恢复方法 |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20210108 |
|
RJ01 | Rejection of invention patent application after publication |