CN102856917B - 一种配电网无功优化方法 - Google Patents
一种配电网无功优化方法 Download PDFInfo
- Publication number
- CN102856917B CN102856917B CN201210269539.8A CN201210269539A CN102856917B CN 102856917 B CN102856917 B CN 102856917B CN 201210269539 A CN201210269539 A CN 201210269539A CN 102856917 B CN102856917 B CN 102856917B
- Authority
- CN
- China
- Prior art keywords
- particle
- max
- population
- wind
- power
- 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
- 238000005457 optimization Methods 0.000 title claims abstract description 18
- 238000000034 method Methods 0.000 title claims abstract description 15
- 239000002245 particle Substances 0.000 claims abstract description 45
- 239000011159 matrix material Substances 0.000 claims description 3
- 230000001105 regulatory effect Effects 0.000 claims description 3
- 230000008859 change Effects 0.000 claims description 2
- 230000008569 process Effects 0.000 claims description 2
- 238000004364 calculation method Methods 0.000 abstract description 4
- 230000003044 adaptive effect Effects 0.000 abstract description 2
- 238000013178 mathematical model Methods 0.000 abstract 2
- 239000000243 solution Substances 0.000 description 13
- 230000005611 electricity Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 239000003990 capacitor Substances 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000001172 regenerating effect Effects 0.000 description 1
Classifications
-
- 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
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开一种配电网无功优化方法,包括以下几个步骤:(1)建立风机的数学模型;(2)初始化节点电压;(3)建立无功规划数学模型,转化成多目标优化问题;(4)输入原始数据;(5)形成初始粒子群,初始化粒子速度和位置;(6)潮流计算,得到粒子适应值和当前最优值;(7)用小生境权重飞行时间法修正粒子速度和位置;(8)计算整个种群适应度,得到适应度值,并更新当前最优解pBestid;(9)对个体粒子中越界点的控制变量进行调整,修正越界状态变量;;(10)判断是否满足终止条件,若满足则输出结果,否则返回(6)。
Description
技术领域
本发明属于电力信息技术领域,其涉及一种配电网无功优化方法,其基小生境权重飞行粒子群优化算法。
背景技术
风能是廉价的可再生能源,基于风力的发电技术越来越受到人们的关注。由于风存在随机性,间歇性和不可调控等性质,当风力发电并网的时候,必然会对电力***的运行造成不可估量的影响。
优化运行是智能电网建设的一项重要任务,风电的接入为地区电网运行的优化提供了新的调节手段。目前,我国含风电场的电网都存在着问题:电压质量不够理想,负荷峰谷时电压波动过大。而影响电压水平重要因素的是无功功率,无功功率的合理分布是保证电压质量的基本条件。一方面,无功过剩将会抬高***电压,过剩的无功在***中流动,不仅会损耗有功功率,占用输电线路和设备容量,而且会在线路上产生电压降落,影响电压质量和稳定性。另一方面,无功不足会导致***电压降低,用电设备不能充分利用,过低的电压水平甚至会导致电压崩溃等严重事故。所以控制无功的合理流动,保持无功平衡,不仅能保证电压质量,提高***运行的安全性和稳定性,而且能降低电能损耗,获得经济效益。因此,无功优化作为含风电场电网规划的一个重要组成部分,通过无功补偿,可以实现电网电压控制、改善电网稳定性、减少网络损耗及保证有较宽的运行裕度。
发明内容
为了克服现有技术的缺陷,本发明公开了一种配电网无功优化方法,其建立了无功规划数学模型,并能够防止粒子后期在最优解附近“振荡”,并提高了普通算法的速度,且结构简单、计算量小、可操作较高,并能有效的提高风电***的节点电压水平,减小网损。
本发明的发明构思是:建立无功补偿在风电***中优化配置的数学模型,以无功设备投资和***有功网损的综合费用最省作为目标函数,同时将节点电压越限和发电机无功出力越限以罚函数的方式进行处理。对建立的模型,采用基于小生境权重飞行粒子群算法(NWFPSO)进行补偿容量的选取优化,以不但能有效地提高算法的收敛精度和收敛速度,使***的电压安全性也得到提高。
本发明解决技术问题所采取的技术方案是:该方法包括以下几个步骤:
(1)给定风电场风速,并在此风速下确定风电机输出的功率;
(2)在含风电场的电力***中形成节点导纳矩阵,设定各节点电压初值;
(3)建立无功规划的数学模型,包括目标函数、功率约束方程;
(4)输入原始数据,获取***节点信息和支路信息,获取控制变量的个数及各自的取值范围;
(5)将上述数学模型的解看成一个粒子,形成初始种群,获取粒子群的群体规模等参数,设置最大迭代次数,随机产生全部粒子及其自身的初始位置和初始速度;
(6)对粒子群中的每个粒子带入进行潮流计算,计算群中每个粒子的网损,得到各粒子的适应度值和当前最优解pBestij;
(7)带入小生境权重飞行法(NWFPSO),调整粒子的速度和位置,计算整个种群适应度。得到各粒子的适应度值和当前最优解pBestij,并跟上次比较,更新最优;
(8)对个体粒子中越界点的控制变量进行调整,修正越界状态变量;
(9)若粒子群迭代次数已经达到最大迭代次数或粒子群迄今为止搜索到的最优位置满足预定的最小适应阈值,则执行步骤(10);若不满足,则回到步骤(6);
(10)停止迭代,输出最优解和无功配置。
与现有技术相比,本发明的有益效果如下:
1.建立了无功规划数学模型
2.能够防止粒子后期在最优解附近“振荡”
3.提高了普通算法的速度
4.结构简单、计算量小、可操作较高
5.有效的提高风电***的节点电压水平,减小网损
附图说明
图1为本发明具体实施例一种配电网无功优化方法的流程图;
图2为本发明具体实施例配电网***图拓扑结构示意图。
具体的实施方式
下方结合附图和具体实施例对本发明做进一步的描述:
实施例
参看图1,一种配电网无功优化方法,包括以下步骤:
步骤(1):给定风电场风速,并在此风速下确定风电机输出的功率;
风电场输出功率的变化源于风速和风向的波动。风速服从威布尔分布,如式(1)所示:
式中:φ(V)为风速V的概率密度;k为形状***;c为尺度系数。
风电场发电功率与风速的关系:
式中:Vco为切除风速;Vci为切入风速;Vr为风电机组额定风速;Pr为风发电机组的额定功率;Vw为风机轮毂高度的风速。
对于异步风力发电机,其输出到电网中的有功功率计算公式如下:
式中:Xk=X1+X2,Pe>0,S<0。指的是异步风力发电机定子电抗;指的是异步风力发电机转子电抗;指的是折算后的转子电阻;u为异步风力发电机端电压。
步骤(2):在含风电场的电力***中形成节点导纳矩阵,设定各节点电压初值;
步骤(3):建立无功规划数学模型
无功规划目标函数为有功网损最小,如式(4)所示:
式中:Ploss为全网的有功网损;Gij、δij分别为节点i、j之间的电导和电压相角差;Vi、Vj为节点i、j的电压幅值。
节点功率平衡方程式为等式约束,如式(5)和(6)所示:
式中Pi,Qi为节点i注入的有功、无功功率;Bij为节点i、j之间的电纳;N表示与节点i直接相连的节点。
步骤(4):获取***节点信息和支路信息,获取控制变量的个数及各自的取值范围,
步骤(5):将无功功率的数学模型的解看成一个粒子,形成初始种群x=[x1,x2,...,xb]T,获取粒子群的群体规模等参数,设置最大迭代次数,随机产生全部粒子及其自身的初始位置和初始速度;
步骤(6):对粒子群中的每个粒子进行潮流计算,计算群中每个粒子的网损,得到各粒子的适应度值fi和当前个体最优解pBestij;
其中:
fi=Ploss (7)
若fi>fi+1,则需要更新各粒子的适应度值和当前最优解pBestij,否则fi+1=fi,则当前最优解pBestij不变。
步骤(7):运用小生境权重飞行时间粒子群算法,调整粒子的速度和位置计算得到新的适应度值,并更新最优解pBestij。
算法如下所示:
w=w0-Δw·t/tmax (9)
式中:t为当前迭代次数,c1为学习因子,取值一般为2;r1是介于[0,1]的随机数;vij∈[-vmax,vmax],vmax为粒子进化过程中得最大速度;w0为权重系数的初始值;Δw权重系数变化率;t为当前迭代次数;tmax为最大迭代次数;Tt表示第t代粒子的飞行时间;T0表示粒子最长飞行时间;t0为比例系数,起调节作用。
步骤(8):个体粒子中越界点的控制变量进行调整,修正越界状态变量;
式(12)中:UGK,min、UGK,max为可调发电机的上、下限;Ci,min、Ci,max为补偿电容器投切组数的上、下限;Tj,min、Tj,max为可调变压器分接头档位的上、下限;NG为所有可调发电机节点总数;NC为所有无功补偿节点总数;NT为所有变压器支路总数。Ui,min、Ui,max为PQ节点电压的上,下限;QGI,min、QGI,max为发电机无功出力的上,下限;NL为PQ节点总数。
步骤(9):若粒子群迭代次数已经达到最大迭代次数或粒子群迄今为止搜索到的最优位置满足预定的最小适应阈值,则执行步骤(10);若不满足,则回到步骤(7);
步骤(10):停止迭代,输出最优解和无功配置。
以国内某配电网为例,配电网***图拓扑结构如图2所示,配电网共17个节点,节点1为发电机,(2,3)、(6,7)、(11,12)为变压器支路。将风电场接入***的节点12。并联电容器分别安装于节点2、节点3、节点7和节点11。设出线首端根节点的电压为恒定值1.0kV。用计算程序进行了检验计算来验证基于改进的粒子群优化算法在风电***无功规划中的应用。无功补偿前后电压结果如表1所示:
表1电压数据表
节点 | PSO | NWFPSO |
1 | 1.0000 | 1.0000 |
2 | 0.9660 | 0.9685 |
3 | 0.9790 | 0.9820 |
4 | 0.9583 | 0.9689 |
5 | 0.9452 | 0.9643 |
6 | 0.9488 | 0.9612 |
7 | 0.9567 | 0.9705 |
8 | 0.9654 | 0.9705 |
9 | 0.9648 | 0.9712 |
10 | 0.9532 | 0.9623 |
11 | 0.9467 | 0.9567 |
12 | 0.9356 | 0.9695 |
13 | 0.9467 | 0.9667 |
14 | 0.9654 | 0.9648 |
15 | 0.9564 | 0.9629 |
16 | 0.9581 | 0.9629 |
17 | 0.9648 | 0.9705 |
表1中编号为对应节点编号。如表1所示,通过小生境权重飞行粒子群优化群算法补偿后各节点电压的改善度大于普通粒子群算法后各节点电压的改善,并且其中电压较低的节点电压改善越明显,如节点11、12;电压较高的节点改善较小,如节点3。数据表示采用自适应权重飞行粒子群优化算法后***电压水平及网损情况都得到了一定的提高。
本发明优选实施例只是用于帮助阐述本发明。优选实施例并没有详尽叙述所有的细节,也不限制该发明仅为所述的具体实施方式。显然,根据本说明书的内容,可作很多的修改和变化。本说明书选取并具体描述这些实施例,是为了更好地解释本发明的原理和实际应用,从而使所属技术领域技术人员能很好地利用本发明。本发明仅受权利要求书及其全部范围和等效物的限制。
Claims (1)
1.一种配电网无功优化方法,其基于小生境权重飞行时间粒子群优化算法,其特征在于,包括以下步骤:
(1)建立风机的数学模型,给定风电场风速,并在此风速下确定风机输出的功率, 其中,V为风速,Vco为切除风速;Vci为切入风速;Vr为风机组额定风速;Pr为风机的额定功率;Vw为风机轮毂高度的风速;
(2)在含风电场的电力***中形成节点导纳矩阵,设定各节点电压初值;
(3)建立无功规划的数学模型,将此数学模型看作是一个多目标优化问题,其中包括目标函数、功率约束方程;
(4)获取***节点信息和支路信息,获取控制变量的个数及各自的取值范围;
(5)将无功功率的数学模型的解看成一个粒子,形成初始种群,获取粒子群的群体规模等参数,设置最大迭代次数,随机产生全部粒子及其自身的初始位置和初始速度;
(6)对粒子群中的每个粒子带入进行潮流计算,计算群中每个粒子的网损,得到各粒子的适应度值和当前最优解pBestid;
(7)运用小生境权重飞行时间粒子群优化算法,调整粒子的速度和位置计算得到新的适应度值,并更新最优解pBestij;
w=w0-Δw·t/tmax;
式中:t为当前迭代次数,c1为学习因子;r1是介于[0,1]的随机数;Vij∈[-Vmax,Vmax],Vmax为粒子进化过程中得最大速度;w0为权重系数的初始值;Δw权重系数变化率;tmax为最大迭代次数;Tt表示第t代粒子的飞行时间;T0表示粒子最长飞行时间;t0为比例系数,起调节作用;
(8)对个体粒子中越界点的控制变量进行调整,修正越界状态变量;
(9)若粒子群迭代次数已经达到最大迭代次数或粒子群迄今为止搜索到的最优位置满足预定的最小适应阈值,则执行步骤(10);若不满足,则回到步骤(6);
(10)停止迭代,输出最优解和无功配置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210269539.8A CN102856917B (zh) | 2012-07-31 | 2012-07-31 | 一种配电网无功优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210269539.8A CN102856917B (zh) | 2012-07-31 | 2012-07-31 | 一种配电网无功优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102856917A CN102856917A (zh) | 2013-01-02 |
CN102856917B true CN102856917B (zh) | 2014-11-19 |
Family
ID=47403226
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210269539.8A Expired - Fee Related CN102856917B (zh) | 2012-07-31 | 2012-07-31 | 一种配电网无功优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102856917B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106130861A (zh) * | 2016-06-15 | 2016-11-16 | 北京工业大学 | 一种基于小生境的多目标粒子群优化的虚拟网络映射方法 |
Families Citing this family (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103337861A (zh) * | 2013-05-31 | 2013-10-02 | 上海交通大学 | 基于黄金混沌小生境粒子群算法的配电网无功优化方法 |
CN103368186A (zh) * | 2013-06-14 | 2013-10-23 | 上海交通大学 | 一种风电***无功优化方法 |
CN103514570B (zh) * | 2013-08-14 | 2016-08-10 | 国家电网公司 | 一种含分布式电源的配电网扩展规划综合优化方法 |
CN104036331B (zh) * | 2014-05-30 | 2017-09-01 | 国家电网公司 | 一种基于改进粒子群算法的电力***动态经济调度方法 |
CA2955008A1 (en) * | 2014-07-17 | 2016-01-21 | 3M Innovative Properties Company | Systems and methods for maximizing expected utility of signal injection test patterns in utility grids |
CN104701849B (zh) * | 2015-03-02 | 2017-10-31 | 清华大学 | 一种主动配电网全分布式自律电压控制方法 |
CN105207274B (zh) * | 2015-10-29 | 2017-12-08 | 许继集团有限公司 | 一种自适应调节无功输出的分布式光伏发电控制方法 |
CN105842581B (zh) * | 2016-03-23 | 2019-05-03 | 国网山东省电力公司章丘市供电公司 | 一种配电网故障预警方法以及*** |
CN106684885B (zh) * | 2016-08-30 | 2020-07-28 | 上海交通大学 | 基于多场景分析的风电机组配电网无功优化方法 |
CN106292291B (zh) * | 2016-10-14 | 2019-08-02 | 南方电网科学研究院有限责任公司 | 一种电网自动发电控制***控制器参数优化方法 |
CN107465197B (zh) * | 2017-08-01 | 2020-07-21 | 国网江苏省电力公司南京供电公司 | 一种基于动态多种群粒子群算法的配电网无功优化方法 |
CN108448658A (zh) * | 2018-02-07 | 2018-08-24 | 上海电力学院 | 一种主动配电网最优运行调度方法 |
CN110518644B (zh) * | 2018-05-21 | 2021-05-04 | 北京天诚同创电气有限公司 | 风电场的无功优化方法、装置及风电场 |
CN109614688B (zh) * | 2018-12-07 | 2022-11-29 | 新奥数能科技有限公司 | 一种负荷管理模型的优化方法、装置、介质及电子设备 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2000116003A (ja) * | 1998-09-29 | 2000-04-21 | Kansai Electric Power Co Inc:The | 電圧無効電力制御方法 |
CN102280889B (zh) * | 2011-06-03 | 2013-10-09 | 上海电力学院 | 基于克隆-粒子群混杂算法的电力***无功优化方法 |
-
2012
- 2012-07-31 CN CN201210269539.8A patent/CN102856917B/zh not_active Expired - Fee Related
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106130861A (zh) * | 2016-06-15 | 2016-11-16 | 北京工业大学 | 一种基于小生境的多目标粒子群优化的虚拟网络映射方法 |
CN106130861B (zh) * | 2016-06-15 | 2019-09-27 | 北京工业大学 | 一种基于小生境的多目标粒子群优化的虚拟网络映射方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102856917A (zh) | 2013-01-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102856917B (zh) | 一种配电网无功优化方法 | |
CN103401248B (zh) | 一种含风电场配电网随机无功优化方法 | |
CN102684201B (zh) | 一种基于电压越限概率的含风电场电网无功优化方法 | |
Haque | Evaluation of power flow solutions with fixed speed wind turbine generating systems | |
CN105322535B (zh) | 含统一潮流控制器的电力***两阶段最优潮流计算方法 | |
CN103368186A (zh) | 一种风电***无功优化方法 | |
CN106026113A (zh) | 一种具有无功自动补偿的微电网***的监控方法 | |
CN102856918A (zh) | 一种基于小生境粒子群算法的配电网的无功优化方法 | |
CN103593711B (zh) | 一种分布式电源优化配置方法 | |
CN107666155A (zh) | 基于Markov模型的多能互补***随机稳定性分析方法 | |
CN103346573B (zh) | 基于黄金分割云粒子群优化算法的风电***无功规划方法 | |
CN103094920A (zh) | 一种直驱式风电机组风电场等值方法 | |
CN110206686A (zh) | 一种用于风力发电机组的自适应最大功率跟踪控制方法 | |
CN103986193B (zh) | 一种最大风电并网容量获取的方法 | |
CN104617578B (zh) | 一种含风电场电力***的可用输电能力的获取方法 | |
CN103455688A (zh) | 建立通用风力发电机组动态仿真模型的方法 | |
CN106712032A (zh) | 一种考虑风电机组有功电压调节能力的最优潮流模型构建方法 | |
CN106786766B (zh) | 一种基于p2g技术的提高风电最大并网容量的方法 | |
CN105958530A (zh) | 一种具有无功自动补偿的微电网*** | |
CN103956767B (zh) | 一种考虑尾流效应的风电场并网稳定性分析方法 | |
CN106655281A (zh) | 一种可再生能源分散接入配电网的监控装置 | |
CN104037806B (zh) | 一种基于风力发电机组基本模型的电力***潮流计算方法 | |
CN112487752B (zh) | 一种基于最优潮流的储能电站选址方法 | |
CN105576712B (zh) | 一种基于自适应极点配置的风电场有功功率在线控制方法 | |
CN109586300A (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 | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20141119 Termination date: 20170731 |
|
CF01 | Termination of patent right due to non-payment of annual fee |