CN105098788B - 一种时序递进的无功电压优化调度控制方法 - Google Patents
一种时序递进的无功电压优化调度控制方法 Download PDFInfo
- Publication number
- CN105098788B CN105098788B CN201510530930.2A CN201510530930A CN105098788B CN 105098788 B CN105098788 B CN 105098788B CN 201510530930 A CN201510530930 A CN 201510530930A CN 105098788 B CN105098788 B CN 105098788B
- Authority
- CN
- China
- Prior art keywords
- reactive
- load
- voltage
- discrete
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 47
- 230000000750 progressive effect Effects 0.000 title claims abstract description 15
- 238000005457 optimization Methods 0.000 claims abstract description 42
- 230000006870 function Effects 0.000 claims abstract description 25
- 230000008859 change Effects 0.000 claims abstract description 13
- 230000009977 dual effect Effects 0.000 claims abstract description 12
- 239000000243 solution Substances 0.000 claims description 29
- 230000000630 rising effect Effects 0.000 claims description 11
- 239000003990 capacitor Substances 0.000 claims description 10
- 230000009471 action Effects 0.000 claims description 9
- 230000005611 electricity Effects 0.000 claims description 8
- 238000002347 injection Methods 0.000 claims description 8
- 239000007924 injection Substances 0.000 claims description 8
- 238000013178 mathematical model Methods 0.000 claims description 8
- 230000005540 biological transmission Effects 0.000 claims description 6
- 238000010079 rubber tapping Methods 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 5
- 238000012545 processing Methods 0.000 claims description 5
- 239000005367 kimax Substances 0.000 claims description 4
- 230000009466 transformation Effects 0.000 claims description 4
- 230000033228 biological regulation Effects 0.000 claims description 3
- 230000008569 process Effects 0.000 abstract description 6
- 230000011218 segmentation Effects 0.000 abstract description 6
- 238000005096 rolling process Methods 0.000 abstract description 4
- 230000008878 coupling Effects 0.000 abstract description 3
- 238000010168 coupling process Methods 0.000 abstract description 3
- 238000005859 coupling reaction Methods 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 4
- 238000009826 distribution Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 238000000354 decomposition reaction Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 206010003830 Automatism Diseases 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000004888 barrier function Effects 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000005312 nonlinear dynamic Methods 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012772 sequence design Methods 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 238000012800 visualization Methods 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
- Control Of Electrical Variables (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明提供一种时序递进的无功电压优化调度控制方法;其从日前优化、到日内修正、再到实时控制的时序递进;优化过程中考虑系统的动态特性,负荷按变化趋势即单调性分段,不同分段优化目标函数不同,兼顾了系统的经济性和安全性。同时采用基于非线性互补的原对偶内点法进行多时段耦合的求解。本发明的方法可根据电网动态变化趋势进行滚动无功优化控制,解决AVC只针对当前电网状态进行优化控制的问题,使得电网以更优的状态运行;同时使得无功优化控制更加合理,保障电网运行的经济性、安全性及可靠性;并保证了连续变量和离散变量共存的大规模非线性优化问题的计算收敛性和求解的全局最优性,对电网进行精细化无功电压控制具有重要的指导意义。
Description
技术领域
本发明涉及自动化技术领域,具体涉及一种时序递进的无功电压优化调度控制方法。
背景技术
电力系统是由发电厂、送变电线路、供配电所和用电等环节组成的电能生产与消费系统。它的功能是将自然界的一次能源通过发电动力装置转化成电能,再经输电、变电和配电将电能供应到各用户。为实现这一功能,电力系统在各个环节和不同层次还具有相应的信息与控制系统,对电能的生产过程进行测量、调节、控制、保护、通信和调度,以保证用户获得安全、优质的电能。而无功优化是在电力系统安全稳定运行、降低网损、提高电压质量、减少无功设备调节次数等多目标之间进行分解协调控制;在空间维上,在高压、中压、低压三个电压等级之间进行分解协调;在时间维上,经济性体现在长过程,电压安全性体现在短过程,需要进行不同时间尺度的调节。
目前,对电网进行无功优化通常使用AVC自动电压控制方法,该方法只针对当前电网状态进行优化控制,而无法根据电网动态变化趋势进行滚动的无功优化调度控制,从而导致无功设备频繁反复调节、完整调度周期内优化控制效果不明显等问题。
因此,如何设计使得无功优化控制更加合理、同时保障电网运行的经济性、安全性以及可靠性的无功电压优化调度控制方法,是本领域的技术人员亟待解决的问题。
发明内容
有鉴于此,本发明提供的一种时序递进的无功电压优化调度控制方法,该方法可根据电网动态变化趋势进行滚动无功优化控制,解决AVC只针对当前电网状态进行优化控制的问题,使得电网以趋于更优的状态运行;同时使得无功优化控制更加合理,保障电网运行的经济性、安全性以及可靠性;并保证了连续变量和离散变量共存的大规模非线性优化问题的计算收敛性和求解的全局最优性。
本发明的目的是通过以下技术方案实现的:
一种时序递进的无功电压优化调度控制方法,所述方法用于对无功电压系统中的无功电压控制设备进行优化调度;所述无功电压控制设备包括连续无功调节设备和离散无功调节设备,所述连续无功调节设备包括发电机、SVC和SVG设备,所述离散无功调节设备包括有载调压变压器抽头、电容器组和电抗器组;所述控制方法包括如下步骤:
步骤1.日前基于负荷单调性将其未来24小时的预测负荷分段;并根据所述负荷的变化趋势不同对所述无功电压系统运行的稳定性影响不同,对其进行分时段多目标优化,得到各时段的连续和离散无功调节设备在每一负荷分段内的动作值;
步骤2.日内进行超短期负荷预测,根据其预测结果对日前的短期负荷预测数据进行修正,并重新按负荷分段进行分时段多目标优化;得到最优目标值以及各时段各调节设备的动作值,进而对日前的无功优化结果进行修正;
步骤3.在实时控制阶段,在日前、日内无功电压控制方案基础上,根据所述无功电压系统实时状态决定离散无功调节设备的投切时刻,协调所述离散与连续无功调节设备进行优化控制,实现无功电压精细化调节;
步骤4.根据非线性互补原对偶内点法求解所述无功电压系统中的离散变量的优化值,完成对所述无功电压系统的无功电压的优化控制。
优选的,所述步骤1,包括:
1-1.当所述负荷预测曲线的变化趋势为负荷低谷期时,以所述无功电压系统网损最小为目标进行优化;
1-2.当所述负荷预测曲线的变化趋势为负荷上升期和高峰期时,以所述无功电压系统网损最小和电压稳定性最高为目标函数进行优化;
1-3.当所述负荷预测曲线的变化趋势为负荷下降期时,以所述无功电压系统网损最小和无功备用最多为目标进行优化;
1-4.优化数学模型的等式和不等式约束条件,包含离散无功调节设备和连续无功调节设备的约束条件;
1-5.优化求解得到日前的最优目标值以及连续和离散无功调节设备在每一负荷分段内的动作值。
优选的,所述步骤1-1,包括:
当所述负荷预测曲线的变化趋势为负荷低谷期时,以所述无功电压系统网损最小为目标函数即F1进行优化:
其中:
式(1)中Ploss为所述无功电压系统网损,Ploss.0为网损初始值,t1为负荷低谷期的时间;式(2)中n为系统的节点数,Vi和Vj为节点i和j的电压值,Gij和Bij分别为节点导纳矩阵的电导和电纳,θij为节点i和j的相角差。
优选的,所述步骤1-2,包括:
当所述负荷预测曲线的变化趋势为负荷上升期和高峰期时,以所述无功电压系统网损最小和电压稳定性最高为目标函数即F2进行优化:
其中:
式(3)中Dv为所述无功电压系统电压稳定指标,Dv.0为电压稳定指标初始值,α、β分别为网损和电压稳定指标的权重系数,tu、th分别为负荷上升期和高峰期的时间;式(4)中n为所述无功电压系统的支路数,中Vi和Vj为支路两端节点i和j的电压值,R和X为相应线路上的电阻和电抗,P′ij和Q′ij为线路传输的有功功率和无功功率。
优选的,所述步骤1-3,包括:
当所述负荷预测曲线的变化趋势为负荷下降期时,以所述无功电压系统网损最小和无功备用最多为目标函数即F3进行优化:
其中:
式(5)中Qback为所述无功电压系统无功备用,即为各个无功调节设备的最大出力与当前出力之差的总和,Qback.0为所述无功电压系统无功备用初始值,α、γ分别为网损和无功备用的权重系数;td为负荷下降期的时间;
式(6)中Qg.back为发电机的无功备用,Qc.back为无功补偿设备的无功备用,m为所述无功电压系统发电机台数,Qgi.max为发电机i所能发的最大无功,Qgi为其当前所发无功,Qcj.max为无功补偿设备j所能提供的最大无功补偿量,Qcj为其当前提供的无功补偿量。
优选的,所述步骤1-4,包括:
优化数学模型的等式和不等式约束条件为:
式(7)中Pi为节点i的注入有功功率,Qi为节点i的注入无功功率,ki、kimin和kimax为变压器分接头档位以及其上下限,δij为节点j与节点i的相角差;Qg、Qgmin和Qgmax为发电机无功及其上下限,Bc、Bcmin和Bcmax为无功补偿电纳及其上下限,Ui、Uimin和Uimax为节点电压及其上下限,Ncimax为无功补偿设备的日调节次数Nci的上限值,Timax为变压器档位的日调节次数Ti的上限值。
优选的,所述步骤3,包括:
3-1.根据系统实时状态进行优化决定离散无功调节设备的投切时刻;
3-2.确定该周期所有离散设备动作方案后进行调节,然后在该周期内将离散设备作为常量,针对所述连续无功调节设备进行优化调节;
3-3.当离散调节设备达到在该负荷分段内的目标值后,在该负荷分段内将其当作常量不再进行调节,利用连续调节设备进行无功电压优化控制。
优选的,所述步骤4,包括:
4-1.非线性互补原对偶内点法需要对离散变量以及连续变量进行求解,先将变压器抽头和电容器组数等离散变量按连续变量处理,进行无功优化计算,对其处理如式(7)所示;ki、Bc为离散变量,将其当作连续量给出其上下限约束进行优化;
4-2.按4-1方法优化计算可得离散变量的两界,然后用所得结果作为初始解,以离散变量的两界构造其非线性互补约束条件,并进行松弛化处理以保证计算收敛性,进行优化求解,如式(8)所示;
式(8)中kin、ki(n+1)为变压器抽头当作连续变量优化得到的两个边界,Qcin、Qci(n+1)为电容器当作连续变量优化得到的两个边界,Qci为无功补偿设备i的无功补偿量;ε为松弛变量,随着优化迭代次数的增大,ε也在逐渐减小;按式(8)处理可求得离散变量的最终动作值。
从上述的技术方案可以看出,本发明提供了一种从日前优化、到日内修正、再到实时控制的时序递进无功电压优化调度控制方法;其优化过程中考虑系统的动态特性,负荷按变化趋势即单调性分段,不同分段优化目标函数不同,即负荷上升期和高峰期以降低网损和提高电压稳定裕度为目标,负荷下降期以降低网损和提高无功储备为目标,负荷低谷期以降低网损为目标;同时采用基于非线性互补的原对偶内点法进行多时段耦合的求解,可根据电网动态变化趋势进行滚动无功优化控制,解决AVC只针对当前电网状态进行优化控制的问题,使得电网以趋于更优的状态运行;同时使得无功优化控制更加合理,保障电网运行的经济性、安全性以及可靠性;并保证了连续变量和离散变量共存的大规模非线性优化问题的计算收敛性和求解的全局最优性,对电网进行精细化无功电压控制具有重要的指导意义。
与最接近的现有技术比,本发明提供的技术方案具有以下优异效果:
1、本发明所提供的技术方案提出的日前优化、到日内修正、再到实时控制的时序递进的无功电压优化控制方法,可根据电网动态变化趋势进行滚动无功优化控制,解决AVC只针对当前电网状态进行优化控制的问题,使得电网以趋于更优的状态运行。
2、本发明所提供的技术方案,考虑负荷不同变化趋势给电网带来的影响,分时段进行多目标无功优化,使得无功优化控制更加合理,保障电网运行的经济性、安全性以及可靠性。
3、本发明所提供的技术方案,利用非线性互补原对偶内点法进行多时段耦合的求解,保证了连续变量和离散变量共存的大规模非线性优化问题的计算收敛性和求解的全局最优性。
4、本发明所提供的技术方案,本发明丰富了自动电压控制的理论体系,对电网进行精细化无功电压控制具有重要的指导意义。
5、本发明提供的技术方案,应用广泛,具有显著的社会效益和经济效益。
附图说明
图1是本发明的一种时序递进的无功电压优化调度控制方法的流程图;
图2是本发明的方法的步骤1的流程示意图;
图3是本发明的方法的步骤4的流程示意图;
图4是本发明的方法的步骤5的流程示意图;
图5是本发明的一种时序递进的无功电压优化调度控制方法的应用例示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
如图1所示,本发明提供一种时序递进的无功电压优化调度控制方法,方法用于对无功电压系统中的无功电压控制设备进行优化调度;无功电压控制设备包括连续无功调节设备和离散无功调节设备,连续无功调节设备包括发电机、SVC和SVG设备,离散无功调节设备包括有载调压变压器抽头、电容器组和电抗器组;
包括如下步骤:
步骤1.日前基于负荷单调性将其未来24小时的预测负荷分段;并根据负荷的变化趋势不同对无功电压系统运行的稳定性影响不同,对其进行分时段多目标优化,得到各时段的连续和离散无功调节设备在每一负荷分段内的动作值;
步骤2.日内进行超短期负荷预测,根据其预测结果对日前的短期负荷预测数据进行修正,并重新按负荷分段进行分时段多目标优化;得到最优目标值以及各时段各调节设备的动作值,进而对日前的无功优化结果进行修正;
步骤3.在实时控制阶段,在日前、日内无功电压控制方案基础上,根据无功电压系统实时状态决定离散无功调节设备的投切时刻,协调离散与连续无功调节设备进行优化控制,实现无功电压精细化调节;
步骤4.根据非线性互补原对偶内点法求解无功电压系统中的离散变量的优化值,完成对无功电压系统的无功电压的优化控制。
如图2所示,步骤1,包括:
1-1.当负荷预测曲线的变化趋势为负荷低谷期时,以无功电压系统网损最小为目标进行优化;
1-2.当负荷预测曲线的变化趋势为负荷上升期和高峰期时,以无功电压系统网损最小和电压稳定性最高为目标函数进行优化;
1-3.当负荷预测曲线的变化趋势为负荷下降期时,以无功电压系统网损最小和无功备用最多为目标进行优化;
1-4.优化数学模型的等式和不等式约束条件,包含离散无功调节设备和连续无功调节设备的约束条件;
1-5.优化求解得到日前的最优目标值以及连续和离散无功调节设备在每一负荷分段内的动作值。
其中,步骤1-1,包括:
当所述负荷预测曲线的变化趋势为负荷低谷期时,以所述无功电压系统网损最小为目标函数即F1进行优化:
其中:
式(1)中Ploss为所述无功电压系统网损,Ploss.0为网损初始值,t1为负荷低谷期的时间;式(2)中n为系统的节点数,Vi和Vj为节点i和j的电压值,Gij和Bij分别为节点导纳矩阵的电导和电纳,θij为节点i和j的相角差。
其中,所述步骤1-2,包括:
当所述负荷预测曲线的变化趋势为负荷上升期和高峰期时,以所述无功电压系统网损最小和电压稳定性最高为目标函数即F2进行优化:
其中:
式(3)中Dv为所述无功电压系统电压稳定指标,Dv.0为电压稳定指标初始值,α、β分别为网损和电压稳定指标的权重系数,tu、th分别为负荷上升期和高峰期的时间;式(4)中n为所述无功电压系统的支路数,中Vi和Vj为支路两端节点i和j的电压值,R和X为相应线路上的电阻和电抗,P′ij和Q′ij为线路传输的有功功率和无功功率。
其中,所述步骤1-3,包括:
当所述负荷预测曲线的变化趋势为负荷下降期时,以所述无功电压系统网损最小和无功备用最多为目标函数即F3进行优化:
其中:
式(5)中Qback为所述无功电压系统无功备用,即为各个无功调节设备的最大出力与当前出力之差的总和,Qback.0为所述无功电压系统无功备用初始值,α、γ分别为网损和无功备用的权重系数;td为负荷下降期的时间;
式(6)中Qg.back为发电机的无功备用,Qc.back为无功补偿设备的无功备用,m为所述无功电压系统发电机台数,Qgi.max为发电机i所能发的最大无功,Qgi为其当前所发无功,Qcj.max为无功补偿设备j所能提供的最大无功补偿量,Qcj为其当前提供的无功补偿量。
其中,所述步骤1-4,包括:
优化数学模型的等式和不等式约束条件为:
式(7)中Pi为节点i的注入有功功率,Qi为节点i的注入无功功率,ki、kimin和kimax为变压器分接头档位以及其上下限,δij为节点j与节点i的相角差;Qg、Qgmin和Qgmax为发电机无功及其上下限,Bc、Bcmin和Bcmax为无功补偿电纳及其上下限,Ui、Uimin和Uimax为节点电压及其上下限,Ncimax为无功补偿设备的日调节次数Nci的上限值,Timax为变压器档位的日调节次数Ti的上限值。
如图3所示,步骤3,包括:
3-1.根据系统实时状态进行优化决定离散无功调节设备的投切时刻;
3-2.确定该周期所有离散设备动作方案后进行调节,然后在该周期内将离散设备作为常量,针对连续无功调节设备进行优化调节;
3-3.当离散调节设备达到在该负荷分段内的目标值后,在该负荷分段内将其当作常量不再进行调节,利用连续调节设备进行无功电压优化控制。
如图4所示,步骤4,包括:
4-1.非线性互补原对偶内点法需要对离散变量以及连续变量进行求解,先将变压器抽头和电容器组数等离散变量按连续变量处理,进行无功优化计算,对其处理如式(7)所示;ki、Bc为离散变量,将其当作连续量给出其上下限约束进行优化;
4-2.按4-1方法优化计算可得离散变量的两界,然后用所得结果作为初始解,以离散变量的两界构造其非线性互补约束条件,并进行松弛化处理以保证计算收敛性,进行优化求解,如式(8)所示;
式(8)中kin、ki(n+1)为变压器抽头当作连续变量优化得到的两个边界,Qcin、Qci(n+1)为电容器当作连续变量优化得到的两个边界,Qci为无功补偿设备i的无功补偿量;ε为松弛变量,随着优化迭代次数的增大,ε也在逐渐减小;按式(8)处理可求得离散变量的最终动作值。
如图5所示,本发明提供一种时序递进的无功电压优化调度控制方法的应用例,本应用例中,优化算法是基于MATLAB编程实现,最终得到优化目标值以及连续和离散无功调节设备的动作值。所涉及到的软件介绍如下:
MATLAB是由美国mathworks公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中,为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案,并在很大程度上摆脱了传统非交互式程序设计语言(如C、Fortran)的编辑模式,代表了当今国际科学计算软件的先进水平。
本应用例的详细过程如下:
步骤1:日前将未来24小时的负荷预测曲线根据其变化趋势按单调性分段,负荷变化趋势不同对系统的影响也不同,因此本发明进行分时段多目标优化。具体如下:
1)在负荷低谷期系统的稳定性较好,此时以传统的优化目标即网损最小为目标函数进行优化,如式(1)所示。
式(1)中n为系统的节点数,Vi和Vj为节点i和j的电压值,Gij和Bij分别为节点导纳矩阵的电导和电纳,θij为节点i和j的相角差。
2)负荷上升期和高峰期除了需要考虑运行的经济性外,系统的电压稳定性也在接受考验,在这两个时段内以系统经济性和稳定性为目标进行优化,即以网损更小和电压稳定性更好为目标函数。由于本专利用到的电压稳定指标也是越小稳定性越高,所以以两个目标之和最小为总目标函数。二者的量纲和数量级不同,本发明以两者各自的初始值为基准值进行标幺化处理统一量纲,至于两目标权重系数的分配不是本发明研究的内容,只是进行了简单处理,两者以1:1比例相加。目标函数为:
其中
式(2)中Ploss为系统网损,Ploss.0为网损初始值,Dv为电压稳定指标,Dv.0为电压稳定指标初始值。式(3)中Vi和Vj为节点i和j的电压值,R和X为相应线路上的电阻和电抗,P′ij和Q′ij为线路传输的有功功率和无功功率。
3)负荷下降期线路传输功率不断减小,此时系统稳定性较好,此时除了考虑经济性,将无功储备设为目标,使得当电网运行出现状况时系统可调控的无功资源越多,更有利于进行无功电压控制。此时无功备用是越多越有利于调控,因此总目标函数为网损与无功备用之差最小,目标函数处理与第2)部分相同。总目标函数为
其中
式(4)中Ploss为系统网损,Ploss.0为网损初始值,Qback为系统无功备用,Qback.0为系统无功备用初始值。式(5)中Qg.back为发电机的无功备用,Qc.back为无功补偿设备的无功备用,m为系统发电机台数,Qgi.max为发电机i所能发的最大无功,Qgi为其当前所发无功,Qcj.max为无功补偿设备j所能提供的最大无功补偿量,Qcj为其当前提供的无功补偿量。
各负荷分段优化目标函数确定后,优化数学模型的等式约束条件以及不等式约束条件如下:
式(6)中Pi为节点i的注入有功功率,Qi为节点i的注入无功功率,kimin和kimax为变压器分接头档位的上下限,Qgmin和Qgmax为发电机无功的上下限,Bcmin和Bcmax为无功补偿电纳的上下限,Uimin和Uimax为节点电压的上下限,Nc.imax为无功补偿设备的日调节次数限值,Timax为变压器档位的日调节次数限值。
优化求解得到最优目标值以及连续和离散无功调节设备在每一负荷分段内的动作值。
步骤2:日内每y分钟进行一次超短期负荷预测,对日前的负荷预测曲线进行不断修正,并重新按负荷分段进行无功优化计算,优化数学模型同上,不再赘述,优化后修正日前各时段连续和离散调节设备的动作值。
步骤3:实时控制层面目的是对连续与离散无功调节设备进行协调控制,对无功电压进行更精细化的调控。本发明假设实时优化控制周期为x分钟(如5分钟)。
由于离散无功调节设备比如电容器或电抗器需分组投切,有载调压变压器档位不能连续调节,因此在实时控制阶段,根据系统实时状态决定离散无功调节设备的投切时刻,直至相应设备达到该分段内的优化动作值。假设根据日内优化结果在某分段内某离散无功调节设备需要动作m次。每x分钟根据系统当前运行状态进行一次优化,若连续两个周期该离散设备优化差值不小于1,则该离散设备动作一次;若其优化差值小于1,则该离散设备不动作。在确定该周期所有离散设备动作方案后进行调节,然后在该周期x分钟内固定离散设备动作值,只针对发电机等连续无功调节设备进行优化求解,得到该周期内连续调节设备优化值进行调节。当离散调节设备达到该负荷分段内的优化动作值后,该负荷分段内离散调节设备当作常量,对连续调节设备进行优化求解进行无功电压控制,这样既避免了离散设备的频繁调节,又保证了当离散设备不能调节时系统的优化控制。
步骤4:基于MATLAB软件编写无功优化算法—非线性互补原对偶内点法。内点法是制约函数法的一种,其整个计算过程均在可行域内进行,即保证每次无约束非线性优化的解点都在可行域之内,直至得出极小点为止。原-对偶内点法是基于对数障碍函数的内点法,它在保持解的原始可行性和对偶可行性的同时,沿着原-对偶路径找到目标函数的最优解。本发明所用算法先将变压器抽头和电容器组数等离散变量按连续化处理,进行无功优化计算,快速寻求离散变量的两界,用所得结果作为初始解,以离散变量的两界构造其非线性互补约束条件进行优化求解,最终得到离散变量的优化值。IEEE14节点算例仿真验证了本发明提出的时序递进的无功电压优化控制方法的有效性。
以上实施例仅用以说明本发明的技术方案而非对其限制,尽管参照上述实施例对本发明进行了详细的说明,所属领域的普通技术人员依然可以对本发明的具体实施方式进行修改或者等同替换,而这些未脱离本发明精神和范围的任何修改或者等同替换,其均在申请待批的本发明的权利要求保护范围之内。
Claims (8)
1.一种时序递进的无功电压优化调度控制方法,其特征在于,所述方法用于对无功电压系统中的无功电压控制设备进行优化调度;所述无功电压控制设备包括连续无功调节设备和离散无功调节设备,所述连续无功调节设备包括发电机、SVC和SVG设备,所述离散无功调节设备包括有载调压变压器抽头、电容器组和电抗器组;
所述控制方法包括如下步骤:
步骤1.日前基于负荷单调性将其未来24小时的预测负荷分段;并根据所述负荷的变化趋势不同对所述无功电压系统运行的稳定性影响不同,对其进行分时段多目标优化,得到各时段的连续和离散无功调节设备在每一负荷分段内的动作值;
步骤2.日内进行超短期负荷预测,根据其预测结果对日前的短期负荷预测数据进行修正,并重新按负荷分段进行分时段多目标优化;得到最优目标值以及各时段各调节设备的动作值,进而对日前的无功优化结果进行修正;
步骤3.在实时控制阶段,在日前、日内无功电压控制方案基础上,根据所述无功电压系统实时状态决定离散无功调节设备的投切时刻,协调所述离散与连续无功调节设备进行优化控制,实现无功电压精细化调节;
步骤4.根据非线性互补原对偶内点法求解所述无功电压系统中的离散变量的优化值,完成对所述无功电压系统的无功电压的优化控制。
2.如权利要求1所述的方法,其特征在于,所述步骤1,包括:
1-1.当所述负荷预测曲线的变化趋势为负荷低谷期时,以所述无功电压系统网损最小为目标进行优化;
1-2.当所述负荷预测曲线的变化趋势为负荷上升期和高峰期时,以所述无功电压系统网损最小和电压稳定性最高为目标函数进行优化;
1-3.当所述负荷预测曲线的变化趋势为负荷下降期时,以所述无功电压系统网损最小和无功备用最多为目标进行优化;
1-4.优化数学模型的等式和不等式约束条件,包含离散无功调节设备和连续无功调节设备的约束条件;
1-5.优化求解得到日前的最优目标值以及连续和离散无功调节设备在每一负荷分段内的动作值。
3.如权利要求2所述的方法,其特征在于,步骤1-1,包括:
当所述负荷预测曲线的变化趋势为负荷低谷期时,以所述无功电压系统网损最小为目标函数即F1进行优化:
其中:
式(1)中Ploss为所述无功电压系统网损,Ploss.0为网损初始值,t1为负荷低谷期的时间;式(2)中n为系统的节点数,Vi和Vj为节点i和j的电压值,Gij和Bij分别为节点导纳矩阵的电导和电纳,θij为节点i和j的相角差。
4.如权利要求2所述的方法,其特征在于,步骤1-2,包括:
当所述负荷预测曲线的变化趋势为负荷上升期和高峰期时,以所述无功电压系统网损最小和电压稳定性最高为目标函数即F2进行优化:
其中:
式(3)中Dv为所述无功电压系统电压稳定指标,Dv.0为电压稳定指标初始值,α、β分别为网损和电压稳定指标的权重系数,tu、th分别为负荷上升期和高峰期的时间;式(4)中n为所述无功电压系统的支路数,式(4)中Vi和Vj为支路两端节点i和j的电压值,R和X为相应线路上的电阻和电抗,Pij'和Qi'j为线路传输的有功功率和无功功率。
5.如权利要求2所述的方法,其特征在于,步骤1-3,包括:
当所述负荷预测曲线的变化趋势为负荷下降期时,以所述无功电压系统网损最小和无功备用最多为目标函数即F3进行优化:
其中:
式(5)中Qback为所述无功电压系统无功备用,即为各个无功调节设备的最大出力与当前出力之差的总和,Qback.0为所述无功电压系统无功备用初始值,α、γ分别为网损和无功备用的权重系数;td为负荷下降期的时间;
式(6)中Qg.back为发电机的无功备用,Qc.back为无功补偿设备的无功备用,m为所述无功电压系统发电机台数,Qgi.max为发电机i所能发的最大无功,Qgi为其当前所发无功,Qcj.max为无功补偿设备j所能提供的最大无功补偿量,Qcj为其当前提供的无功补偿量。
6.如权利要求2所述的方法,其特征在于,步骤1-4,包括:
优化数学模型的等式和不等式约束条件为:
式(7)中Pi为节点i的注入有功功率,Qi为节点i的注入无功功率,ki、kimin和kimax为变压器分接头档位以及其上下限,δij为节点j与节点i的相角差;Qg、Qgmin和Qgmax为发电机无功及其上下限,Bc、Bcmin和Bcmax为无功补偿电纳及其上下限,Ui、Uimin和Uimax为节点电压及其上下限,Ncimax为无功补偿设备的日调节次数Nci的上限值,Timax为变压器档位的日调节次数Ti的上限值。
7.如权利要求1所述的方法,其特征在于,所述步骤3,包括:
3-1.根据系统实时状态进行优化决定离散无功调节设备的投切时刻;
3-2.确定周期所有离散设备动作方案后进行调节,然后在周期内将离散设备作为常量,针对所述连续无功调节设备进行优化调节;
3-3.当离散调节设备达到在该负荷分段内的目标值后,在该负荷分段内将其当作常量不再进行调节,利用连续调节设备进行无功电压优化控制。
8.如权利要求1所述的方法,其特征在于,所述步骤4,包括:
4-1.非线性互补原对偶内点法需要对离散变量以及连续变量进行求解,先将变压器抽头和电容器组数等离散变量按连续变量处理,进行无功优化计算,对其处理如式(7)所示;ki、Bc为离散变量,将其当作连续量给出其上下限约束进行优化;
4-2.按4-1方法优化计算可得离散变量的两界,然后用所得结果作为初始解,以离散变量的两界构造其非线性互补约束条件,并进行松弛化处理以保证计算收敛性,进行优化求解,如式(8)所示;
式(8)中kin、ki(n+1)为变压器抽头当作连续变量优化得到的两个边界,Qcin、Qci(n+1)为电容器当作连续变量优化得到的两个边界,Qci为无功补偿设备i的无功补偿量;ε为松弛变量,随着优化迭代次数的增大,ε也在逐渐减小;按式(8)处理可求得离散变量的最终动作值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510530930.2A CN105098788B (zh) | 2015-08-26 | 2015-08-26 | 一种时序递进的无功电压优化调度控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510530930.2A CN105098788B (zh) | 2015-08-26 | 2015-08-26 | 一种时序递进的无功电压优化调度控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105098788A CN105098788A (zh) | 2015-11-25 |
CN105098788B true CN105098788B (zh) | 2018-07-20 |
Family
ID=54578655
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510530930.2A Active CN105098788B (zh) | 2015-08-26 | 2015-08-26 | 一种时序递进的无功电压优化调度控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105098788B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109755943A (zh) * | 2019-03-22 | 2019-05-14 | 中国电建集团青海省电力设计院有限公司 | 一种基于新能源接入的无功电压时序递进优化控制方法 |
CN109861307A (zh) * | 2019-03-26 | 2019-06-07 | 国网北京市电力公司 | 电力系统的发电机调节方法及装置 |
CN110768265A (zh) * | 2019-11-26 | 2020-02-07 | 华北电力大学 | 一种考虑时序的配电网调度方法 |
CN111130150B (zh) * | 2019-12-17 | 2021-03-19 | 国网冀北电力有限公司 | 基于IPM-intPSO的两阶段动态无功优化方法 |
CN112086983B (zh) * | 2020-07-28 | 2022-02-15 | 国网湖南省电力有限公司 | 一种特高压直流输电系统的运行电压优化方法 |
CN112734590B (zh) * | 2020-11-18 | 2023-12-19 | 国网青海省电力公司 | 基于内点法和变领域搜索的动态规划法的二阶段区域电网无功优化方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102013697A (zh) * | 2010-10-13 | 2011-04-13 | 中电普瑞科技有限公司 | 一种风电场的电压/无功综合智能控制方法 |
CN104065079A (zh) * | 2014-07-14 | 2014-09-24 | 东方电子股份有限公司 | 一种风电离散和连续无功补偿的协调控制方法 |
-
2015
- 2015-08-26 CN CN201510530930.2A patent/CN105098788B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102013697A (zh) * | 2010-10-13 | 2011-04-13 | 中电普瑞科技有限公司 | 一种风电场的电压/无功综合智能控制方法 |
CN104065079A (zh) * | 2014-07-14 | 2014-09-24 | 东方电子股份有限公司 | 一种风电离散和连续无功补偿的协调控制方法 |
Non-Patent Citations (2)
Title |
---|
基于负荷预测的地区电网电压无功控制方案;于汀,郭瑞鹏;《电力系统保护与控制》;20120616;第40卷(第12期);第121-124,155页 * |
电网AGC与AVC协调控制方法;于汀 等;《电力系统保护与控制》;20150801;第43卷(第15期);第42-47页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105098788A (zh) | 2015-11-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105098788B (zh) | 一种时序递进的无功电压优化调度控制方法 | |
CN108683179B (zh) | 基于混合整数线性规划的主动配电网优化调度方法及系统 | |
CN105226664B (zh) | 一种主动配电网无功电压分层分布协调控制方法 | |
CN101752903B (zh) | 一种时序递进式电力调度方法 | |
Ma et al. | Reactive power optimization in power system based on improved niche genetic algorithm | |
CN108023364B (zh) | 基于凸差规划的配电网分布式电源最大接入能力计算方法 | |
CN103248059A (zh) | 配电网无功电压优化方法与系统 | |
CN110690732A (zh) | 一种光伏无功分区计价的配电网无功优化方法 | |
CN108321810A (zh) | 抑制光伏并网点电压波动的配网多时间尺度无功控制方法 | |
CN113516278B (zh) | 有源配电网多时间尺度有功无功协调优化调度方法及系统 | |
CN106992526A (zh) | 一种分层分布式的光伏接入配电网电压控制系统及方法 | |
CN105633949A (zh) | 一种经济运行实时发电调度和控制方法 | |
CN109390953A (zh) | 含分布式电源和电动汽车的低压配电网无功电压协调控制方法和系统 | |
CN114362267B (zh) | 考虑多目标优化的交直流混合配电网分散式协调优化方法 | |
CN104701867A (zh) | 基于分支定界法以及原对偶内点法的日前无功优化方法 | |
CN103595061A (zh) | 基于综合效益分析的企业电网无功优化方法及系统 | |
CN107947185A (zh) | 一种基于多目标优化的电网自动电压控制方法及系统 | |
CN115313403A (zh) | 一种基于深度强化学习算法的实时电压调控方法 | |
CN103490433A (zh) | 一种配电网无功优化方法 | |
CN113067344A (zh) | 一种基于模型预测控制的主动配电网无功优化方法 | |
CN108365637A (zh) | 一种送电计划与抽水储能发电计划优化方法及系统 | |
CN115481856A (zh) | 考虑综合需求响应的综合能源系统多尺度调度方法及系统 | |
CN106169760B (zh) | 主站avc系统与子站svg系统的协调控制方法 | |
CN109066769B (zh) | 风电全消纳下虚拟电厂内部资源调度控制方法 | |
CN108365610A (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 |