CN105896565B - 基于比重变异粒子群算法的配电网无功优化方法 - Google Patents
基于比重变异粒子群算法的配电网无功优化方法 Download PDFInfo
- Publication number
- CN105896565B CN105896565B CN201610394699.3A CN201610394699A CN105896565B CN 105896565 B CN105896565 B CN 105896565B CN 201610394699 A CN201610394699 A CN 201610394699A CN 105896565 B CN105896565 B CN 105896565B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- particle
- value
- max
- 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.)
- Expired - Fee Related
Links
- 239000002245 particle Substances 0.000 title claims abstract description 122
- 238000005457 optimization Methods 0.000 title claims abstract description 37
- 238000000034 method Methods 0.000 title claims abstract description 26
- 230000035772 mutation Effects 0.000 title claims abstract description 12
- 238000004364 calculation method Methods 0.000 claims description 21
- 230000008602 contraction Effects 0.000 claims description 6
- 230000001174 ascending effect Effects 0.000 claims description 3
- 230000008859 change Effects 0.000 claims description 2
- 238000004422 calculation algorithm Methods 0.000 abstract description 24
- 230000002068 genetic effect Effects 0.000 description 4
- 235000019606 astringent taste Nutrition 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000003990 capacitor Substances 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000001627 detrimental effect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000006855 networking Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000001172 regenerating effect Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000009897 systematic effect Effects 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
- H02J3/18—Arrangements for adjusting, eliminating or compensating reactive power in networks
-
- 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
-
- 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
- Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
- Y02E40/30—Reactive power compensation
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
- Feedback Control In General (AREA)
Abstract
本发明基于传统的粒子群算法加入了比重变异的方法,当找寻最优解达到一定条件时,可根据每个粒子的情况对粒子的位置进行变异,使粒子在搜索过程中避免陷入局部最优;另外,本发明还对对传统线性变化的惯性权重有所改进,根据反馈值的大小动态的改变,当反馈值大于平均值时,w值将变大,粒子搜索范围扩大,反之,将在较小范围搜索;这样可以更加精确的对整个粒子空间进行搜索,收敛速度也得到提高。经验证,基于比重变异粒子群的配电网无功优化方法,可结合配电网现有可控设备对配电网进行高效管理,在运行时能保证更高的收敛性和优化程度。
Description
技术领域
本发明涉及一种配电网无功优化方法,具体涉及一种基于粒子群算法的配电网无功优化方法。
背景技术
近些年我国大力提倡可再生能源的使用与开发,分布式电源的接入和电动汽车充电桩建设逐步增多。这些新变化对电网系统的电压、损耗、可靠性都带来不利的影响。为了降低这些影响,确保电网的供电质量和系统经济安全运行,需要对我国电网中已经安装的有载调压的变压器和无功补偿装置进行调节和控制。根据算法要求,将可以控制的有载可调变压器,无功补偿装置和发电机端电压量化后作为控制变量,将网络损耗或者节点电压作为优化目标,通过相应的算法找出最优设置,从而达到优化配电网系统的目的。
传统的无功优化方法包括:粒子群算法算法(particle swarm optimization,PSO)、遗传算法(Genetic Algorithm,GA)等,但是传统的算法存在收敛性差,容易陷入局部最优等不足之处。
发明内容
为了克服上述算法容易陷入局部最优、收敛性差的和优化程度不够等缺点。本发明采用了基于比重变异粒子群算法,该算法能够有效的避免陷入局部最优,同时进一步的提高了收敛速度和优化深度。
本发明提供的一种基于比重变异粒子群算法的配电网无功优化方法,具体如下:
建立配电网的无功优化数学模型,即选取系统的有功网损作为适应函数,公式如下:
上式中:PL为网络有功损耗,Vi、Vj分别表示节点i和j的电压,δij表示节点i和节点j之间的相角差,Gij表示节点i和节点j之间的电导;控制变量包括发电机机端电压、变压器档位、无功补偿等控制变量等。
然后再进行如下步骤:
S1:对配电网无功优化系统进行初始化,获取配电网系统的节点信息和支路信息;设置控制变量的个数,并对各控制变量的取值范围、第一学习因子c1、第二学习因子c2、最大迭代次数Nmax以及初始种群的规模T进行设置;
S2:在设置完成的控制变量取值范围内,随机初始化后得到一个种群数目为T的粒子群,并获得初始粒子群每个粒子的位置X(i)和速度V(i);
S3:对初始粒子群进行潮流计算,计算所述初始粒子群中每个粒子的网损PL(i),取每个粒子的当前网损PL(i)作为每个个体的初始最优解Pi,并找出最小网损作为全局最优解Pg;
S4:判断当前迭代次数是否达到设定最大值Nmax,若达到最大值跳转到步骤S11,否则转到步骤S5;
S5:根据迭代更新公式对当前种群每个粒子的速度和位置进行更新,生成一个新粒子群;
S6:对新粒子群的每个粒子速度和位置进行评价,判断每个粒子的位置和速度是否超过设定值;如果超过设定值,则对该粒子进行修正,并取其设定值作为该粒子的位置和速度;
S7:对新粒子群进行潮流计算,获得每个粒子的网损;比较每个粒子更新前后的网损,取较小的为该个体最优解,同时找出新的网损最小值,并和上次的全局最优解比较,保留较小值作为新全局最优解;
S8:判断全局最优解是否保持n代没有变化;若是,转到步骤S9,否则转到步骤S4;n值可选择取6、7、8任一个数值;
S9:对粒子群粒子的位置进行比重变异;
S10:对变异后的每个粒子速度和位置按S6所述的方法进行评价和修正,再对变异粒子群按步骤S7所述方法进行潮流计算并找出最新的个体最优值和全局最优值,然后跳转到步骤S4;
S11:停止迭代计算,保存个体最优值和全局最优值,同时保存全局最优时对应各控制变量的值;此时,全局最优即为优化理想结果,控制变量值即为最优设置值。
对于上述步骤:
进一步,所述步骤S9中比重变异的具体步骤如下:
S91:将需要变异粒子对应的适应函数值用MATLAB自带sort函数进行升序排列,由此得到对应的粒子顺序;
S92:将新序列的粒子带入适应函数求出每个粒子的函数反馈值PL(i),即有功网损值;
S93:计算出每个粒子的适应度F(i):
S94:将每个粒子按照其适应度比重进行相应变异:X(i)*=X(i)+F(i)·randn;
上式中:PL(i)为对应粒子i的网损值;randn是MATLAB自带的正态分布随机函数。
进一步,所述步骤S5中,粒子的速度和位置的迭代更新公式如下:
vi,d,k+1=L(wvi,d,k+c1rand1(Pi,k-xi,d,k)+c2rand2(Pg,k-xi,d,k)),
xi,d,k+1=xi,d,k+vi,d,k+1
X(i)=[VG,i,QC,i,TN,i]=[VG1,i...VGn,i,Qc1,i...Qcn2,i,TN1,i...TNn3,i]
上式中:k和k+1表示迭代次数,i表示粒子编号,d表示第d维数分量;x和v分别表示粒子位置和速度;c1和c2为学习因子;rand1和rand2为随机变量;Pi,k和Pg,k分别为粒子个体最优解和全局最优解;w为惯性权重,L为收缩因子;其中,第一学习因子c1他第二学习因子c2影响粒子向个体最优和整体最优粒子的飞行速度;
其中,惯性权重w计算公式如下:
式中:wmax、wmin分别为惯性权重w的最大值和最小值,一般分别取0.9和0.4;Nmax和N分别为最大更新次数和当前更新次数,PL(i)为有功网损,Pave为整个群体平均有功网损,rand为随机函数。
其中,收缩因子L计算公式如下:
第一学习因子c1、第二学习因子c2计算公式如下:
上式中:c1和c2为学习因子,Nmax和N分别为最大更新次数和当前更新次数。
进一步,所述步骤S2中,配电系统中控制参数对应于粒子的位置如下式所示:
X=[VG,QC,TN]=[VG1...VGn1,Qc1...Qcn2,TN1...TNn3]
其中,第i代粒子位置是:
X(i)=[VG,i,QC,i,TN,i]=[VG1,i...VGn,i,Qc1,i...Qcn2,i,TN1,i...TNn3,i]
上式中:VG表示发电机机端电压,QC表示无功补偿量,TN表示有载可调变压器档位。
进一步,所述步骤S6中,对粒子群的每个粒子速度和位置进行评价,即保持:
其中,Vmin、Vmax,Qmin、Qmax、Tmin、Tmax分别表示电压、补偿电容和变压器调节的上下限;由于在实际优化过程中无功补偿量和可调变压器的档位都不是连续变化的,就需要对其评价并进行离散化,即,通过自带函数将原来的非整数四舍五入变为整数。
有益效果:
本发明通过加入比重变异,可以根据每个粒子的情况对粒子进行变异,使粒子在搜索过程中避免陷入局部最优,使优化结果更好改进粒子群算法具有更好的收敛性和优化程度;并且,该方法对配电网的无功优化也是实际可行的。
本发明还通过对传统线性变化惯性权重公式的改进,根据反馈值的大小动态的改变,当反馈值大于平均值时,w值将变大,这样粒子搜索范围就扩大;相反当反馈值小于平均值时w值就变小,粒子将在较小范围搜索;这样可以更加精确的对整个粒子空间进行搜索,收敛速度也得到提高。
附图说明
图1基于比重变异粒子群算法的流程图;
图2将比重变异粒子群算法的应用于IEEE30系统网络结构图;
图3在IEEE30节点上的比重变异粒子群算法和其他两种算法仿真结果对比图。
具体实施方式
基于比重变异粒子群算法的配电网无功优化方法,具体实施方式如下:
建立配电网的无功优化数学模型,即选取系统的有功网损作为适应函数,公式如下:
上式中:PL为网络有功损耗,Vi、Vj分别表示节点i和j的电压,δij表示节点i和节点j之间的相角差,Gij表示节点i和节点j之间的电导;控制变量包括发电机机端电压、变压器档位、无功补偿等控制变量等。
然后再进行如下步骤:
S1:对配电网无功优化系统进行初始化,获取配电网系统的节点信息和支路信息;设置控制变量的个数,并对各控制变量的取值范围、第一学习因子c1、第二学习因子c2、最大迭代次数Nmax以及初始种群的规模T进行设置;
S2:在设置完成的控制变量取值范围内,随机初始化后得到一个种群数目为T的粒子群,并获得初始粒子群每个粒子的位置X(i)和速度V(i);
配电系统中控制参数对应于粒子的位置如下式所示:
X=[VG,QC,TN]=[VG1...VGn1,Qc1...Qcn2,TN1...TNn3]
其中,第i代粒子位置是:
X(i)=[VG,i,QC,i,TN,i]=[VG1,i...VGn,i,Qc1,i...Qcn2,i,TN1,i...TNn3,i]
上式中:VG表示发电机机端电压,QC表示无功补偿量,TN表示有载可调变压器档位。
S3:对初始粒子群进行潮流计算,计算所述初始粒子群中每个粒子的网损PL(i),取每个粒子的当前网损PL(i)作为每个个体的初始最优解Pi,并找出最小网损作为全局最优解Pg;
S4:判断当前迭代次数是否达到设定最大值Nmax,若达到最大值跳转到步骤S11,否则转到步骤S5;
S5:根据迭代更新公式对当前种群每个粒子的速度和位置进行更新,生成一个新粒子群,
粒子的速度和位置的迭代更新公式如下:
vi,d,k+1=L(wvi,d,k+c1rand1(Pi,k-xi,d,k)+c2rand2(Pg,k-xi,d,k)),
xi,d,k+1=xi,d,k+vi,d,k+1
X(i)=[VG,i,QC,i,TN,i]=[VG1,i...VGn,i,Qc1,i...Qcn2,i,TN1,i...TNn3,i]
上式中:k和k+1表示迭代次数,i表示粒子编号,d表示第d维数分量;x和v分别表示粒子位置和速度;c1和c2为学习因子;rand1和rand2为随机变量;Pi,k和Pg,k分别为粒子个体最优解和全局最优解;w为惯性权重,L为收缩因子;其中,第一学习因子c1他第二学习因子c2影响粒子向个体最优和整体最优粒子的飞行速度;
其中,惯性权重w计算公式如下:
式中:wmax、wmin分别为惯性权重w的最大值和最小值,一般分别取0.9和0.4;Nmax和N分别为最大更新次数和当前更新次数,PL(i)为有功网损,Pave为整个群体平均有功网损,rand为随机函数。
其中,收缩因子L计算公式如下:
第一学习因子c1、第二学习因子c2计算公式如下:
上式中:c1和c2为学习因子,Nmax和N分别为最大更新次数和当前更新次数。
S6:对新粒子群的每个粒子速度和位置进行评价,判断每个粒子的位置和速度是否超过设定值;如果超过设定值,则对该粒子进行修正,并取其设定值作为该粒子的位置和速度,
对粒子群的每个粒子速度和位置进行评价,即保持:
其中,Vmin、Vmax,Qmin、Qmax、Tmin、Tmax分别表示电压、补偿电容和变压器调节的上下限;由于在实际优化过程中无功补偿量和可调变压器的档位都不是连续变化的,就需要对其评价并进行离散化,即,通过自带函数将原来的非整数四舍五入变为整数。
S7:对新粒子群进行潮流计算,获得每个粒子的网损;比较每个粒子更新前后的网损,取较小的为该个体最优解,同时找出新的网损最小值,并和上次的全局最优解比较,保留较小值作为新全局最优解;
S8:判断全局最优解是否保持6代没有变化;若是,转到步骤S9,否则转到步骤S4;n值还可以根据具体情况选择取6、7、8任一个数值;
S9:对粒子群各粒子的位置进行比重变异,具体步骤如下:
S91:将需要变异粒子对应的适应函数值用MATLAB自带sort函数进行升序排列,由此得到对应的粒子顺序;
S92:将新序列的粒子带入适应函数求出每个粒子的函数反馈值PL(i),即有功网损值;
S93:计算出每个粒子的适应度F(i):
S94:将每个粒子按照其适应度比重进行相应变异:X(i)*=X(i)+F(i)·randn;
上式中:PL(i)为对应粒子i的网损值;randn是MATLAB自带的正态分布随机函数。
S10:对变异后的每个粒子速度和位置按S6所述的方法进行评价和修正,再对变异粒子群按步骤S7所述方法进行潮流计算并找出最新的个体最优值和全局最优值,然后跳转到步骤S4;
S11:停止迭代计算,保存个体最优值和全局最优值,同时保存全局最优时对应各控制变量的值;此时,全局最优即为优化理想结果,控制变量值即为最优设置值。
进一步的,在本发明实施案例在如图2所示IEEE30节点系统进行仿真。IEEE30节点为标准网络系统。IEEE30节点系统有41条支路、21个负荷节点、6台发电机、4台可调变压器以及2个电容补偿点。节点1,2,5,8,11,13为发电机节点,在发电机中,节点1为平衡节点;节点2,5,8,11,13为PV节点;其余为PQ节点。发电机端电压取值范围为0.95-1.05;有载调压变压器调节范围为0.9-1.1,分为8个档位;节点10、24为补偿电容器,补偿范围为0-0.5。
本发明取粒子维度d为10,最大迭代次数Nmax为100,种群数目T为60。将未加入权重变异的粒子群算法算法(particle swarm optimization,PSO)、遗传算法(GeneticAlgorithm,GA)和加入比重变异的PSO算法进行对比,对比过程中保证它们的基本参数相同,即利用3种不同的算法对IEEE30节点进行仿真。经过30次仿真,统计结果如下表1所示。
表1 3种算法有功网损结果统计表
如图3所示,三种算法都有效的降低了系统网损,但基于本发明的加入比重变异和对惯性权利改进的的PSO算法收敛性和优化程度明显优于GA算法和PSO算法。
可见,通过引入比重变异,在对IEEE30节点案例仿真中,基于比重变异粒子群算法具有更好的收敛性和优化程度,同时也说明基于比重变异粒子群算法对配电网的无功优化是实际可行的。
上述虽然结合附图对本发明的具体实施方式进行了描述,但并非对本发明保护范围的限制,所属领域技术人员应该明白,在本发明的技术方案的基础上,本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本发明的保护范围以内。
Claims (7)
1.一种基于比重变异粒子群算法的配电网无功优化方法,其特征在于:建立配电网的无功优化数学模型,选取系统的有功网损作为适应函数,即:
<mrow>
<mi>min</mi>
<mi> </mi>
<msub>
<mi>P</mi>
<mi>L</mi>
</msub>
<mo>=</mo>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</munder>
<mrow>
<mo>(</mo>
<msubsup>
<mi>V</mi>
<mi>i</mi>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>V</mi>
<mi>j</mi>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<mn>2</mn>
<msub>
<mi>V</mi>
<mi>i</mi>
</msub>
<msub>
<mi>V</mi>
<mi>j</mi>
</msub>
<msub>
<mi>cos&delta;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>G</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
上式中:PL为网络有功损耗,Vi、Vj分别表示节点i和j的电压,δij表示节点i和节点j之间的相角差,Gij表示节点i和节点j之间的电导;
然后进行如下步骤:
S1:对配电网无功优化系统进行初始化,获取配电网系统的节点信息和支路信息;设置控制变量的个数,并对各控制变量的取值范围、第一学习因子c1、第二学习因子c2、最大迭代次数Nmax以及初始种群的规模T进行设置;
S2:在设置完成的控制变量取值范围内,随机初始化后得到一个种群数目为T的粒子群,并获得初始粒子群每个粒子的位置X(i)和速度V(i);
S3:对初始粒子群进行潮流计算,计算所述初始粒子群中每个粒子的网损PL(i),取每个粒子的当前网损PL(i)作为每个个体的初始最优解Pi,并找出最小网损作为全局最优解Pg;
S4:判断当前迭代次数是否达到设定最大值Nmax,若达到最大值跳转到步骤S11,否则转到步骤S5;
S5:根据迭代更新公式对当前种群每个粒子的速度和位置进行更新,生成一个新粒子群;粒子的速度和位置的迭代更新公式如下:
vi,d,k+1=L(wvi,d,k+c1rand1(Pi,k-xi,d,k)+c2rand2(Pg,k-xi,d,k)),
xi,d,k+1=xi,d,k+vi,d,k+1,
X(i)=[VG,i,QC,i,TN,i]=[VG1,i...VGn,i,Qc1,i...Qcn2,i,TN1,i...TNn3,i],
上式中:k和k+1表示迭代次数,i表示粒子编号,d表示第d维数分量;x和v分别表示粒子位置和速度;c1和c2为学习因子;rand1和rand2为随机变量;Pi,k和Pg,k分别为粒子个体最优解和全局最优解;w为惯性权重,L为收缩因子;
S6:对新粒子群的每个粒子速度和位置进行评价,判断每个粒子的位置和速度是否超过设定值;如果超过设定值,则对该粒子进行修正,并取其设定值作为该粒子的位置和速度;
S7:对新粒子群进行潮流计算,获得每个粒子的网损;比较每个粒子更新前后的网损,取较小的为该个体最优解,同时找出新的网损最小值,并和上次的全局最优解比较,保留较小值作为新全局最优解;
S8:判断全局最优解是否保持n代没有变化;若是,转到步骤S9,否则转到步骤S4;
S9:对粒子群粒子的位置进行比重变异;具体步骤如下:
S91:将需要变异粒子对应的适应函数值用MATLAB自带sort函数进行升序排列,由此得到对应的粒子顺序;
S92:将新序列的粒子带入适应函数求出每个粒子的函数反馈值PL(i),即有功网损值;
S93:计算出每个粒子的适应度F(i):
S94:将每个粒子按照其适应度比重进行相应变异:X(i)*=X(i)+F(i)·randn;
上式中:PL(i)为对应粒子i的网损值;randn是MATLAB自带的正态分布随机函数;
S10:对变异后的每个粒子速度和位置按S6所述的方法进行评价和修正,再对变异粒子群按步骤S7所述方法进行潮流计算并找出最新的个体最优值和全局最优值,然后跳转到步骤S4;
S11:停止迭代计算,保存个体最优值和全局最优值,同时保存全局最优时对应各控制变量的值。
2.根据权利要求1所述的配电网无功优化方法,其特征在于:惯性权重w计算公式如下:
<mrow>
<mi>w</mi>
<mo>=</mo>
<mrow>
<mo>{</mo>
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>w</mi>
<mi>max</mi>
</msub>
<mo>-</mo>
<mfrac>
<mrow>
<msub>
<mi>w</mi>
<mi>max</mi>
</msub>
<mo>-</mo>
<msub>
<mi>w</mi>
<mi>min</mi>
</msub>
</mrow>
<msub>
<mi>N</mi>
<mi>max</mi>
</msub>
</mfrac>
<mo>,</mo>
<mi>N</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>w</mi>
<mi>max</mi>
</msub>
<mo>-</mo>
<mfrac>
<mrow>
<msub>
<mi>w</mi>
<mi>max</mi>
</msub>
<mo>-</mo>
<msub>
<mi>w</mi>
<mi>min</mi>
</msub>
</mrow>
<msub>
<mi>N</mi>
<mi>max</mi>
</msub>
</mfrac>
<mo>+</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>L</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>a</mi>
<mi>v</mi>
<mi>e</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mi>r</mi>
<mi>a</mi>
<mi>n</mi>
<mi>d</mi>
<mo>,</mo>
<mi>N</mi>
<mo>></mo>
<mn>1</mn>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>,</mo>
</mrow>
</mrow>
</mrow>
上式中:wmax、wmin分别为惯性权重w的最大值和最小值;Nmax和N分别为最大更新次数和当前更新次数,PL(i)为有功网损,Pave为整个群体平均有功网损,rand为随机函数。
3.根据权利要求1所述的配电网无功优化方法,其特征在于:收缩因子L计算公式如下:
第一学习因子c1、第二学习因子c2计算公式如下:
<mrow>
<msub>
<mi>c</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<mi>N</mi>
<msub>
<mi>N</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
</mfrac>
<mo>,</mo>
<msub>
<mi>c</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mn>1</mn>
<mo>+</mo>
<mfrac>
<mi>N</mi>
<msub>
<mi>N</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
</mfrac>
<mo>,</mo>
</mrow>
上式中:c1和c2为学习因子,Nmax和N分别为最大更新次数和当前更新次数。
4.根据权利要求1所述的配电网无功优化方法,其特征在于:所述步骤S2中,配电系统中控制参数对应于粒子的位置如下式所示:
X=[VG,QC,TN]=[VG1...VGn1,Qc1...Qcn2,TN1...TNn3]
其中,第i代粒子位置是:
X(i)=[VG,i,QC,i,TN,i]=[VG1,i...VGn,i,Qc1,i...Qcn2,i,TN1,i...TNn3,i]
上式中:VG表示发电机机端电压,QC表示无功补偿量,TN表示有载可调变压器档位。
5.根据权利要求1所述的配电网无功优化方法,其特征在于:所述步骤S6中,对粒子群的每个粒子速度和位置进行评价,即保持:
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>V</mi>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>&le;</mo>
<msub>
<mi>V</mi>
<mrow>
<mi>G</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>&le;</mo>
<msub>
<mi>V</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Q</mi>
<mi>min</mi>
</msub>
<mo>&le;</mo>
<msub>
<mi>Q</mi>
<mrow>
<mi>C</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>&le;</mo>
<msub>
<mi>Q</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>T</mi>
<mi>min</mi>
</msub>
<mo>&le;</mo>
<msub>
<mi>T</mi>
<mrow>
<mi>N</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>&le;</mo>
<msub>
<mi>T</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中,Vmin、Vmax,Qmin、Qmax、Tmin、Tmax分别表示电压、补偿电容和变压器调节的上下限。
6.根据权利要求1所述的配电网无功优化方法,其特征在于:步骤S8中的n值可选择取6、7、8任一个数值。
7.根据权利要求2所述的配电网无功优化方法,其特征在于:惯性权重w的最大值wmax和最小值wmin分别取0.9和0.4。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610394699.3A CN105896565B (zh) | 2016-06-06 | 2016-06-06 | 基于比重变异粒子群算法的配电网无功优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610394699.3A CN105896565B (zh) | 2016-06-06 | 2016-06-06 | 基于比重变异粒子群算法的配电网无功优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105896565A CN105896565A (zh) | 2016-08-24 |
CN105896565B true CN105896565B (zh) | 2018-02-13 |
Family
ID=56710751
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610394699.3A Expired - Fee Related CN105896565B (zh) | 2016-06-06 | 2016-06-06 | 基于比重变异粒子群算法的配电网无功优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105896565B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106777449A (zh) * | 2016-10-26 | 2017-05-31 | 南京工程学院 | 基于二进制粒子群算法的配电网络重构方法 |
CN107272419A (zh) * | 2017-08-01 | 2017-10-20 | 成都雅骏新能源汽车科技股份有限公司 | 一种基于改进pso的驾驶员自适应方向控制方法 |
CN112297936B (zh) * | 2020-10-28 | 2022-03-29 | 杭州市电力设计院有限公司余杭分公司 | 一种电动汽车的充放电控制方法、装置、设备及存储介质 |
CN114389277B (zh) * | 2022-01-12 | 2024-09-06 | 福州大学 | 基于蜉蝣算法的配电网无功优化方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102723721A (zh) * | 2012-05-31 | 2012-10-10 | 西南交通大学 | 基于个体最优位置自适应变异扰动粒子群算法的电力系统无功优化方法 |
US9645558B2 (en) * | 2012-09-29 | 2017-05-09 | Operation Technology, Inc. | Dynamic parameter tuning using particle swarm optimization |
CN103972908A (zh) * | 2014-05-23 | 2014-08-06 | 国家电网公司 | 一种基于自适应混沌粒子群算法的多目标无功优化方法 |
CN104682405B (zh) * | 2015-03-31 | 2018-01-12 | 福州大学 | 一种基于禁忌粒子群算法的配电网无功优化方法 |
-
2016
- 2016-06-06 CN CN201610394699.3A patent/CN105896565B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN105896565A (zh) | 2016-08-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Wong et al. | Optimal placement and sizing of battery energy storage system for losses reduction using whale optimization algorithm | |
CN105896565B (zh) | 基于比重变异粒子群算法的配电网无功优化方法 | |
CN104037776B (zh) | 随机惯性因子粒子群优化算法的电网无功容量配置方法 | |
CN103972908A (zh) | 一种基于自适应混沌粒子群算法的多目标无功优化方法 | |
CN105449675A (zh) | 优化分布式能源接入点和接入比例的电力网络重构方法 | |
CN103701134B (zh) | 基于mcr与电容混合补偿的风电场并网点电压控制方法 | |
CN102856918A (zh) | 一种基于小生境粒子群算法的配电网的无功优化方法 | |
CN105186500B (zh) | 一种基于重加权加速Lagrangian的配电网能量分散协调优化方法 | |
Tang et al. | Study on day-ahead optimal economic operation of active distribution networks based on Kriging model assisted particle swarm optimization with constraint handling techniques | |
CN105243432B (zh) | 基于nsga-ⅲ的upfc选址和容量多目标配置方法 | |
CN114362267B (zh) | 考虑多目标优化的交直流混合配电网分散式协调优化方法 | |
CN105406485B (zh) | 基于改进蝙蝠算法的变电站动态无功优化方法与系统 | |
CN113378100B (zh) | 考虑碳排放的配电网源网荷储协同优化调度模型和方法 | |
CN116388302B (zh) | 一种协调网侧资源的配电网有功-无功联合优化方法 | |
CN110929931A (zh) | 考虑分布式电源和负荷时序特性的配电网协调规划方法 | |
CN103346573B (zh) | 基于黄金分割云粒子群优化算法的风电系统无功规划方法 | |
CN103368186A (zh) | 一种风电系统无功优化方法 | |
CN111277004A (zh) | 一种配电网“源-网-荷”两阶段多目标控制方法及系统 | |
CN105337290B (zh) | 一种适用于电力系统低频振荡辅助决策的无功调整方法 | |
CN104836240B (zh) | 一种基于电网固有结构特性的配电网最优无功补偿方法 | |
CN105720591A (zh) | 一种电力系统无功优化方法及系统 | |
CN109888817B (zh) | 对光伏电站和数据中心进行位置部署和容量规划方法 | |
CN105896547B (zh) | 一种风电接入下的大电网分级电压控制方法 | |
CN106950831A (zh) | 一种离线优化/在线切换的无功补偿方法 | |
CN117808151A (zh) | 一种基于粒子群-遗传融合算法的变电站无功优化方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180213 |