CN102623988A - 一种电网中日前计划潮流的自动生成方法 - Google Patents
一种电网中日前计划潮流的自动生成方法 Download PDFInfo
- Publication number
- CN102623988A CN102623988A CN2012100847517A CN201210084751A CN102623988A CN 102623988 A CN102623988 A CN 102623988A CN 2012100847517 A CN2012100847517 A CN 2012100847517A CN 201210084751 A CN201210084751 A CN 201210084751A CN 102623988 A CN102623988 A CN 102623988A
- Authority
- CN
- China
- Prior art keywords
- electrical network
- plan
- active power
- power
- generating set
- 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.)
- Granted
Links
Images
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明涉及一种电网中日前计划潮流的自动生成方法,属于电力***调度自动化技术领域。本发明将日前计划潮流生成问题分解为有功子问题和无功子问题分步进行求解,并将子问题转化为带约束非线性优化问题,采用内点法进行求解。通过有功子问题初步消除***发电机组发电计划与母线负荷预测数据之间的失配量,协调发电机组发电计划与联络断面功率传输计划之间的不一致,通过无功子问题以及有功再调整步骤进一步消除***发电功率与负荷功率之间的失配量,按照计划值控制联络断面传输的有功功率,给出合理、优化的日前无功电压计划。本发明能有效地为日前计划安全校核模块提供可信的日前计划交流潮流,具有良好的可靠性和实用性。
Description
技术领域
本发明涉及一种电网中日前计划潮流的自动生成方法,属于电力***调度自动化技术领域。
背景技术
电力***是高维、强非线性的复杂动态***,其运行稳定性是电网安全生产的首要问题。长期以来,电网日前运行计划安排是电网调度部门确定电网未来运行方式、指导电力生产的主要依据,并且在向运行计划与静态、暂态安全校核联合的方向发展。合理的日前调度运行计划直接关系到电网日后的安全经济运行。现有的日前计划潮流生成主要根据有功发电计划生成直流潮流解,然后配上典型的发电机机端电压,由此获得日前交流潮流解,这种做法存在以下问题:
(1)收敛性问题:对大规模电网,典型的发电机机端电压与有功计划匹配不理想,可能导致交流潮流迭代计算发散,发散后极难调整,进而导致N-1静态安全校核和其他动态安全校核均无法进行。据一些电网统计,半数以上的日前计划潮流无法收敛,只能转为直流潮流计算。
(2)合理性问题:由于计划数据的多源性和交流潮流的非线性,容易导致有功计划的失配,通常可表现为:联络断面功率交换计划与发电计划的不兼容、计划发电与母线负荷预测的功率不平衡、考虑电网有功损耗后引起的发电负荷不平衡、交流潮流引起的联络断面功率传输计划非线性误差等。与此同时,即使能获得收敛的交流潮流解,由于可能存在不合理的无功环流和电压分布,可能导致过高或过低的发电机内电势,将显著降低动态安全校核结果的可信度。
(3)可行性问题:关键母线电压限制(包括:稳定约束和绝缘过压约束)、分区域无功裕度约束、关口功率因素约束、重要机组进迟相约束,交织一起,无功电压的调整难度极大,难以获得可行的交流潮流解。
因此,如何获得收敛、合理和可行的有功无功联合的计划潮流解,是日前安全校核和计划调整顺利进行的基础和难点问题。
发明内容
本发明的目的是提出一种电网日前计划潮流的自动生成方法,解决现有的基于直流潮流生成日前计划潮流这一方法的不足。依据由电力控制中心的发电计划/检修计划管理子***或者由其他方式生成的管辖范围内电网计划日有功发电计划、电网计划日设备检修计划、电网计划日联络断面有功功率传输计划以及电网各母线负荷预测有功无功值,通过优化方法给出合理的计划日交流潮流解,为后续的安全校核和日前运行状态评估模块提供收敛、合理的计划日计划潮流。
本发明提出的电网中日前计划潮流的自动生成方法,包括以下步骤:
(1)根据计划日之前的电网拓扑结构以及计划日当天发电机组的开停机计划、设备检修计划以及负荷投切计划,生成一个计划日电网拓扑结构;
(2)设电网中按发电计划发电的发电机组为计划发电机组,无发电计划的发电机组为非计划发电机组,并设非计划发电机组的计划发电值为计划日之前的与计划日相同测量时刻点的发电值,建立一个调整计划日发电机有功功率的优化模型:
s.t.PG-PD-B′θ=0
上式中,PGi为电网中各发电机组的有功功率;为计划发电机组安排的计划发电值,为非计划发电机组设定的计划发电值,N1为计划发电机组数目,N2为所有发电机组数目,PTie,k为电网各联络断面传输的有功功率,为电网各联络断面计划传输的有功功率,wi、wi’和wk分别为上述优化模型目标函数的三个部分的加权系数,PG为上述优化模型直流潮流方程组中电网发电机的有功功率,PG=[PG1,PG2,…,PGN2],PD为上述优化模型直流潮流方程组中电网母线负荷预测的有功功率,PD=[Pd1,Pd2,…,PdNd],Nd为电网负荷个数,B’为上述优化模型中的直流潮流方程组导纳矩阵,θm和θn分别为电网中的支路(m,n)首端节点m和末端节点n的电压相角,xmn为电网中的支路(m,n)的电抗;Pmn为电网中的支路(m,n)的有功功率,和分别为电网中的支路(m,n)有功功率的上限值和下限值;PTie,k、和分别为电网联络断面k传输的有功功率、有功功率上限值和有功功率下限值;S为电网的联络断面集;
上述优化模型中的目标函数为三个部分的加权和,三个部分分别为:计划发电机组的有功功率与该发电机组的计划发电值的偏差的平方、非计划发电机组的有功功率与该发电机组的计划发电值的偏差的平方,以及电网中各联络断面的有功功率与联络断面有功功率传输计划值的偏差的平方,前后三个加权系数的比例为:
wi∶wi’∶wk=(10~20)∶1∶(25~100);
上述优化模型的控制变量为电网中发电机组的有功功率;
上述优化模型的约束条件为:
电网的直流潮流方程组,即发电机有功总功率=负荷有功总功率+电网有功功率损耗,
电网中各发电机的有功功率上限值和有功功率下限值,
电网中各输电线路的有功功率极限值以及电网中各变压器绕组的有功功率极限值,
电网中各联络断面的有功功率极限值;
求解上述优化模型,根据计划日的发电机有功功率计划值,得到计划日电网中发电机的有功功率,且计划日发电机的有功功率与计划日的发电机计划发电值之间的偏差最小;
(3)建立一个用于求解计划日电网无功电压分布的优化模型,以得到与上述步骤(2)的计划日电网中发电机的有功功率相匹配的计划日电网中发电机的无功功率和发电机机端电压:
min PGSlack
s.t.f(V,θ)=0
上式中,PGSlack为电网中容量最大的非计划发电机组的有功功率,f(V,θ)=0为电网中节点的交流潮流方程组;Vj为电网内节点j的电压幅值,和分别为电网内节点j的电压幅值的上限值和下限值;PGi为电网内发电机组i的有功功率,和分别为电网内发电机组i的有功功率上限值和下限值;QGi为电网内发电机无功功率,和分别为发电机组i的无功功率的上限值和下限值;Pmn为电网中节点m和节点n之间的支路(m,n)传输的有功功率,及分别为电网中的支路(m,n)传输的有功功率上限值和下限值;Qzonel为电网中各设定管理区域的无功备用,和分别为电网中各设定管理区域无功备用的上限值和下限值;
上述优化模型的目标函数是电网有功损耗最小;
上述优化模型的控制变量为所有发电机组的无功功率、变压器分接头档位、并联电容器及并联电抗器的投退;
上述模型的状态变量为全网节点电压幅值Vi及相角θi;
上述优化模型的约束条件有:
电网中各节点的交流潮流方程,
电网中各支路的有功功率约束,
电网中各节点电压幅值约束,
电网中各发电机的无功功率约束,
电网中各设定管理区域的发电机无功备用约束;
求解上述优化模型,得到电网中各发电机机端电压和发电机无功功率,且电网有功损耗最小;
(4)根据步骤(2)得到的电网发电机的有功功率以及步骤(3)得到的计划日无功电压分布,得到电网中日前计划潮流结果。
上述电网中日前计划潮流的自动生成方法中,在其中的步骤(3)以后,还可以包括以下步骤:
建立一个重新调整电网联络断面的有功功率优化模型:
s.t.dPG-B′dθ=0
上式中,dPGi为再次进行机组有功功率调整时的调整量,dPTie,k为电网各联络断面有功功率的调整量,大于零代表上调,小于零代表下调,为步骤(3)中得到的各联络断面有功功率与有功传输计划值的偏差,dPG为基于直流潮流的发电机组有功功率调整量的向量,dPG=[dPG1,dPG2,…,dPGN2],dθ为基于直流潮流的电网节点相角调整量的向量,dθ=[dθ1,dθ2,…,dθn];
上述优化模型为基于直流潮流的最优增量调节模型,目标函数包括三项:计划机组的有功功率再调整量,非计划机组有功功率再调整量,电网联络断面传输有功的再调整量;
控制变量为电网中各发电机组的有功功率再调整量;
等式约束为基于直流潮流的功率平衡方程;
求解上述优化模型,得到电网中各发电机组的有功功率再调整量;
根据上述步骤(2)得到的电网发电机的有功功率、上述步骤(3)得到的计划日无功电压分布以及上述最后得到的电网中各发电机组的有功功率再调整量,得到电网中日前计划潮流结果。
本发明提出的电网中日前计划潮流的自动生成方法,与已有的在直流潮流基础上配上典型机端电压的方法,具有以下优点:
1、现有的其他日前计划潮流生成方法需要一个收敛的初始基态潮流结果,本发明提出的日前计划潮流自动生成方法,可在不收敛的初始基态潮流上启动,自动给出各约束可行域中的收敛的交流潮流优化解;
2、利用本发明提出的电网中日前计划潮流的自动生成方法,可以实现联络线计划与发电计划的兼容,发电与负荷及电网有功损耗匹配,电网联络断面传输功率与计划值匹配,能够消除不合理的无功环流和电压分布,提高了安全校核结果的可信度;
3、本发明提出的电网中日前计划潮流的自动生成方法,采用了非线性带约束优化模型求解,可计及电网的各种复杂运行约束,如:联络断面计划、管理区域无功备用、各节点电压约束、支路功率约束等,传统的方法无法灵活处理这些约束。
附图说明
图1是使用本发明方法进行的电网日前计划安全校核的流程框图。
具体实施方式
本发明提出的电网中日前计划潮流的自动生成方法,其流程如图1所示,包括以下步骤:
(1)根据计划日之前的电网拓扑结构以及计划日当天发电机组的开停机计划、设备检修计划以及负荷投切计划,生成一个计划日电网拓扑结构;
本发明方法中,计划日之前的电网拓扑结构包括计划日之前电网内所有设备的参数、发电机的运行状态以及开关刀闸的投退状态。对于计划日开机发电的机组以及日前检修计划中检修竣工重新投入电网的设备,需要进行拓扑搜索,从发电机的机端以及设备的两端出发,搜索设备与当前带电电网的连接关系,通过合上搜索到的通电路径上的开关刀闸,将开机发电的机组及检修竣工的设备重新投入到电网中。
(2)鉴于潮流计算中有功和无功存在弱耦合关系,本发明将计划潮流生成问题分解为有功子问题和无功子问题,分为两步进行求解;有功子问题的目标是调整计划日电网中发电机组的有功功率,在初步消除交流潮流计算中的电网有功失配量的同时,使之与发电计划之间的偏差最小,并匹配联络断面功率计划值,该失配量是指发电机组有功发电计划与母线负荷预测值及交流潮流的有功损耗之间的失配;实际应用中,发电机组的有功发电计划、电网母线负荷预测以及电网联络断面功率计划数据往往出自不同的来源,存在不一致的情况,因此需要在有功子问题中进行调整,协调各部分数据的不一致;无功子问题的目标是在有功子问题结果的基础上,给出计划日发电机组无功功率、发电机的机端电压以及电网中其他各节点的电压幅值,从而得到完整的交流潮流解;
设电网中按发电计划发电的发电机组为计划发电机组,无发电计划的发电机组为非计划发电机组,并设非计划发电机组的计划发电值为计划日之前的与计划日相同测量时刻点的发电值,本步骤解决有功子问题,建立一个调整计划日发电机有功功率的优化模型,即对计划日有功发电计划做出合理调整:
s.t.PG-PD-B′θ=0
上式中,PGi为电网中各发电机组的有功功率;为计划发电机组安排的计划发电值,为非计划发电机组设定的计划发电值,N1为计划发电机组数目,N2为所有发电机组数目,PTie,k为电网各联络断面传输的有功功率,为电网各联络断面计划传输的有功功率,wi、wi’和wk分别为上述优化模型目标函数的三个部分的加权系数,PG为上述优化模型直流潮流方程组中电网发电机的有功功率,PG=[PG1,PG2,…,PGN2],PD为上述优化模型直流潮流方程组中电网母线负荷预测的有功功率,PD=[Pd1,Pd2,…,PdNd],Nd为电网负荷个数,B’为上述优化模型中的直流潮流方程组导纳矩阵,θm和θn分别为电网中的支路(m,n)首端节点m和末端节点n的电压相角,xmn为电网中的支路(m,n)的电抗;Pmn为电网中的支路(m,n)的有功功率,和分别为电网中的支路(m,n)有功功率的上限值和下限值;PTie,k、和分别为电网联络断面k传输的有功功率、有功功率上限值和有功功率下限值;S为电网的联络断面集;
上述优化模型中的目标函数为三个部分的加权和,三个部分分别为:计划发电机组的有功功率与该发电机组的计划发电值的偏差的平方、非计划发电机组的有功功率与该发电机组的计划发电值的偏差的平方,以及电网中各联络断面的有功功率与联络断面有功功率传输计划值的偏差的平方,前后三个加权系数的比例为:
wi∶wi∶wk=(10~20)∶1∶(25~100);
一般情况下联络断面的输电计划是首要关心的项目,需要按计划执行,所以联络断面有功功率偏差的权重值最大,其次为计划机组有功功率与计划值偏差的权重值,非计划机组的有功功率与计划值偏差的权重值最小,负责消纳大部分发电计划与母线负荷预测的偏差。
上述优化模型的控制变量为电网中发电机组的有功功率;
上述优化模型的约束条件为:
电网的直流潮流方程组,即发电机有功总功率=负荷有功总功率+电网有功功率损耗,
电网中各发电机的有功功率上限值和有功功率下限值,
电网中各输电线路的有功功率极限值以及电网中各变压器绕组的有功功率极限值,
电网中各联络断面的有功功率极限值;
求解上述优化模型,根据计划日的发电机有功功率计划值,对全网发电机功率进行调整,协调联络断面计划与发电机功率计划之间的不一致,得到计划日电网中发电机的有功功率,且计划日发电机的有功功率与计划日的发电机计划发电值之间的偏差最小;
(3)在有功子问题的基础上给出与有功子问题结果相匹配的计划日电网中发电机的无功功率和发电机机端电压,电力***无功是个区域性的量,不能用处理有功问题的方法来处理无功问题,需要通过分配各分区内无功功率、减小区域间远距离无功输送来优化***有功损耗,提高区域电压水平,该步骤所使用的数学模型如下:
min PGSlack
s.t.f(V,θ)=0
上式中,PGSlack为电网中容量最大的非计划发电机组的有功功率,f(V,θ)=0为电网中节点的交流潮流方程组;Vj为电网内节点j的电压幅值,和分别为电网内节点j的电压幅值的上限值和下限值;PGi为电网内发电机组i的有功功率,和分别为电网内发电机组i的有功功率上限值和下限值;QGi为电网内发电机无功功率,和分别为发电机组i的无功功率的上限值和下限值;Pmn为电网中节点m和节点n之间的支路(m,n)传输的有功功率,及分别为电网中的支路(m,n)传输的有功功率上限值和下限值;Qzonel为电网中各设定管理区域的无功备用,和分别为电网中各设定管理区域无功备用的上限值和下限值;
上述优化模型的目标函数是电网有功损耗最小;
上述优化模型的控制变量为所有发电机组的无功功率、变压器分接头档位、并联电容器及并联电抗器的投退;
上述模型的状态变量为全网节点电压幅值Vi及相角θi;
上述优化模型的约束条件有:
电网中各节点的交流潮流方程,
电网中各支路的有功功率约束,
电网中各节点电压幅值约束,
电网中各发电机的无功功率约束,
电网中各设定管理区域的发电机无功备用约束;
求解上述优化模型,得到电网中各发电机机端电压和发电机无功功率,且电网有功损耗最小;
(4)本发明在求解有功子问题的时候使用直流潮流模型,将电网联络断面传输的有功功率控制在计划交换值,而在求解无功子问题时,采用交流潮流模型,考虑了电网有功损耗的影响,并对PGSlack进行了调整,这将使得联络断面有功功率产生变化,偏离计划值,实际试验表明,会产生6%~8%的偏离,为了更精确地将联络断面传输的有功功率控制在计划值,本步骤设计一个联络断面有功功率再调整模块,在无功子问题结果的基础上重新校准联络断面功率,该步骤的数学模型下:
s.t.dPG-B′dθ=0
上式中,dPGi为再次进行机组有功功率调整时的调整量,dPTie,k为电网各联络断面有功功率的调整量,大于零代表上调,反之则代表下调,为无功子问题求解结束后各联络断面有功功率与有功传输计划值的偏差,dPG为基于直流潮流的发电机组有功功率调整量的向量,dPG=[dPG1,dPG2,…,dPGN2],dθ为基于直流潮流的电网节点相角调整量的向量,dθ=[dθ1,dθ2,…,dθn],其他符号意义与前文相同;
上述优化模型为基于直流潮流的最优增量调节模型,等式约束中不再考虑电网负荷;本步骤的目的是在求解无功子问题得到的交流潮流结果的基础上,通过对发电机有功功率进行二次调整来达到控制联络断面交换的有功功率的目的;由于无功子问题求解结果联络断面交换功率与计划值偏差大约为6%~8%,已经不大,进行本步骤再调整时,各发电机有功功率调整量是很小的;实例试验表明,在发电机有功功率二次调整后执行一次PQ分解法潮流迭代,最终结果联络断面交换功率较计划偏差可控制在1.5%以内。
(5)根据上述步骤(2)得到的电网发电机的有功功率、上述步骤(3)得到的计划日无功电压分布以及上述步骤(4)得到的电网中各发电机组的有功功率再调整量,得到电网中日前计划潮流结果。
本发明方法基于调整后的发电机有功计划以及日前无功电压计划,得到完整的交流潮流解,为日前计划安全校核模块提供所需的日前计划潮流;安全校核结果以及本方法计算过程中对机组功率计划的调整结果将返回到发电计划管理模块,对机组计划进行调整,从而形成闭环的计划校核调整过程,如附图1所示。
本发明将分解的有功子问题和无功子问题转化为优化问题进行求解,采用内点法求解式1)以及式2)的非线性优化问题,式1)以及式2)均可以写为如下非线性规划标准形式:
min f(x)
s.t.h(x)=0 6)
其中:x∈R(n),h(x)≡[h1(x),…,hm(x)]T,g(x)≡[g1(x),…,gr(x)]T。
引入松弛变量(l,u)∈R(r),将6)式变换成带等式约束和简单变量约束的非线性规划形式如下:
minimize f(x)
subjectto h(x)=0
g(x)-l-g=0 7)
(l,u)≥0
定义7)式的拉格朗日函数为:
根据Karush-Kuhn-Tucker定理,最优解必须满足如下KKT***:
(l,u)≥0, y≠0 17)
其中,(L,U,Z,W)∈R(r×r)依据数学规划中的惯例,表示元素为l,u,z,w的对角线矩阵,e≡[1,…,1]T∈R(r)。
将13)、14)式和互补条件15)、16)式合并,再根据扰动KKT理论,引入扰动参数μ>0,松弛合并后的互补条件,将13)-16)替换为:
直接应用牛顿法求解上述扰动KKT***9)-12),18)-19),得到如下修正***:
化简(20-25),消去Δl,Δu,Δz,Δw,得到最终的简约KKT***如下:
其中:
及消去后,Δl,Δu,Δz,Δw的表达式:
27)式中,H(·)由两部分组成:Hh是f(x),h(x)和g(x)海森矩阵的线性组合;Hg称作障碍矩阵(Barrier Matrix);Hg项防止不等式约束g(x)违反它们的双边极限。
由此产生的路径跟踪法内点迭代格式如下:
(0)初始化。置k=0,Kmax,中心参数(centering parameter)σ∈(0,1],收敛精度ε=10-6;选择[l,u]T>0,[z>0,w<0,y=0]T,这里,k,Kmax分别是迭代计数和最大允许迭代次数;
(1)如果k<Kmax,继续;否则,打印不收敛;
(2)计算互补间隙(complementary gap)CGap:
如果CGap<ε,停止,输出最优解;
(3)求解修正***26)-33),获取修正量[Δx,Δy],[Δl,Δu]和[Δz,Δw];
(4)比值检验,确定原始和对偶空间中的最大步长:
i=1,2,…,r
(5)更新原始和对偶变量:
(6)k=k+1,回(1);
以下结合附图2,介绍本发明方法的一个实施例。
选取IEEE39节点测试***对本文提出的算法进行测试,附图2给出IEEE39***接线图。
各发电机的有功发电计划如表1所示,这里以发电机所挂接的母线号作为发电机名;***的母线负荷预测数据如表2所示。
表1发电机有功/无功功率限值及有功发电计划
Tab.1 Active/Reactive output limits and active output schedule of generators
表2电网母线负荷预测有功/无功
Tab.2 Bus bad forecasting
设定三个电网的联络断面,如附图2及表3所示。
表3IEEE39***联络断面定义
Tab.3 Definition of flow gates in IEEE 39 buses system
表3中,LN-39-9代表节点39和节点9之间的输电线路。
对本实例进行计算,各步骤的时间性能如下:
表4IEEE39测试***计算耗时
Tab.4 CPU Time elapsed of IEEE 39 buses system
有功子问题与有功再调整问题的数学模型中等式约束均采用直流潮流方程,为线性约束,模型收敛性好,只需数步即可迭代收敛。无功子问题为非线性约束优化问题,使用预测-校正内点法在可行域内部搜索优化解,具有良好的收敛性能,并应用稀疏技术大大降低计算耗时。
依次求解有功子问题、无功子问题、有功再调整子问题,最终计算结果各发电机组功率情况见表5,联络断面交换功率的变化见表6。
表5IEEE 39***发电机功率计算结果
Tab.5 Generation results of IEEE 39 buses system
表6IEEE 39***联络断面功率结果
Tab.6 Power flow through the flowgates in IEEE 39 buses system
将表5和表1、2作对比,可看到大多数发电机功率结果都较计划值进行了约10.0MW的调整,但是机组bus35、bus36功率反而较计划值上调约50.0MW,机组bus38功率较计划值下调了约40.0MW(表5阴影数字);这是因为,由附图2***单线图可以得知,断面②和断面③为割集断面,调整机组bus38以及机组bus35、bus36的功率将直接影响到两断面传输的有功功率,由于制定的联络断面功率传输计划与发电机有功发电计划存在矛盾;因此计算时算法优先考虑了满足联络断面功率传输计划,对机组bus38以及机组bus35、bus36的功率计划进行了较大的调整;由表6可见,无功子问题后,联络断面传输的有功功率与计划值存在明显偏差,通过有功再调整,这一偏差大大减小了。
Claims (2)
1.一种电网中日前计划潮流的自动生成方法,其特征在于该方法包括以下步骤:
(1)根据计划日之前的电网拓扑结构以及计划日当天发电机组的开停机计划、设备检修计划以及负荷投切计划,生成一个计划日电网拓扑结构;
(2)设电网中按发电计划发电的发电机组为计划发电机组,无发电计划的发电机组为非计划发电机组,并设非计划发电机组的计划发电值为计划日之前的与计划日相同测量时刻点的发电值,建立一个调整计划日发电机有功功率的优化模型:
s.t.PG-PD-B′θ=0
上式中,PGi为电网中各发电机组的有功功率;为计划发电机组安排的计划发电值,为非计划发电机组设定的计划发电值,N1为计划发电机组数目,N2为所有发电机组数目,PTie,k为电网各联络断面传输的有功功率,为电网各联络断面计划传输的有功功率,wi、wi’和wk分别为上述优化模型目标函数的三个部分的加权系数,PG为上述优化模型直流潮流方程组中电网发电机的有功功率,PG=[PG1,PG2,…,PGN2],PD为上述优化模型直流潮流方程组中电网母线负荷预测的有功功率,PD=[Pd1,Pd2,…,PdNd],Nd为电网负荷个数,B’为上述优化模型中的直流潮流方程组导纳矩阵,θm和θn分别为电网中的支路(m,n)首端节点m和末端节点n的电压相角,xmn为电网中的支路(m,n)的电抗;Pmn为电网中的支路(m,n)的有功功率,和分别为电网中的支路(m,n)有功功率的上限值和下限值;PTie,k、和分别为电网联络断面k传输的有功功率、有功功率上限值和有功功率下限值;S为电网的联络断面集;
上述优化模型中的目标函数为三个部分的加权和,三个部分分别为:计划发电机组的有功功率与该发电机组的计划发电值的偏差的平方、非计划发电机组的有功功率与该发电机组的计划发电值的偏差的平方,以及电网中各联络断面的有功功率与联络断面有功功率传输计划值的偏差的平方,前后三个加权系数的比例为:
wi∶wi’∶wk=(10~20)∶1∶(25~100);
上述优化模型的控制变量为电网中发电机组的有功功率;
上述优化模型的约束条件为:
电网的直流潮流方程组,即发电机有功总功率=负荷有功总功率+电网有功功率损耗,
电网中各发电机的有功功率上限值和有功功率下限值,
电网中各输电线路的有功功率极限值以及电网中各变压器绕组的有功功率极限值,
电网中各联络断面的有功功率极限值;
求解上述优化模型,根据计划日的发电机有功功率计划值,得到计划日电网中发电机的有功功率,且计划日发电机的有功功率与计划日的发电机计划发电值之间的偏差最小;
(3)建立一个用于求解计划日电网无功电压分布的优化模型,以得到与上述步骤(2)的计划日电网中发电机的有功功率相匹配的计划日电网中发电机的无功功率和发电机机端电压:
min PGSlack
s.t.f(V,θ)=0
上式中,PGSlack为电网中容量最大的非计划发电机组的有功功率,f(V,θ)=0为电网中节点的交流潮流方程组;Vj为电网内节点j的电压幅值,和分别为电网内节点j的电压幅值的上限值和下限值;PGi为电网内发电机组i的有功功率,和分别为电网内发电机组i的有功功率上限值和下限值;QGi为电网内发电机无功功率,和分别为发电机组i的无功功率的上限值和下限值;Pmn为电网中节点m和节点n之间的支路(m,n)传输的有功功率,及分别为电网中的支路(m,n)传输的有功功率上限值和下限值;Qzonel为电网中各设定管理区域的无功备用,和分别为电网中各设定管理区域无功备用的上限值和下限值;
上述优化模型的目标函数是电网有功损耗最小;
上述优化模型的控制变量为所有发电机组的无功功率、变压器分接头档位、并联电容器及并联电抗器的投退;
上述模型的状态变量为全网节点电压幅值Vi及相角θi;
上述优化模型的约束条件有:
电网中各节点的交流潮流方程,
电网中各支路的有功功率约束,
电网中各节点电压幅值约束,
电网中各发电机的无功功率约束,
电网中各设定管理区域的发电机无功备用约束;
求解上述优化模型,得到电网中各发电机机端电压和发电机无功功率,且电网有功损耗最小;
(4)根据步骤(2)得到的电网发电机的有功功率以及步骤(3)得到的计划日无功电压分布,得到电网中日前计划潮流结果。
2.一种电网中日前计划潮流的自动生成方法,其特征在于该方法包括以下步骤:
(1)根据计划日之前的电网拓扑结构以及计划日当天发电机组的开停机计划、设备检修计划以及负荷投切计划,生成一个计划日电网拓扑结构;
(2)设电网中按发电计划发电的发电机组为计划发电机组,无发电计划的发电机组为非计划发电机组,并设非计划发电机组的计划发电值为计划日之前的与计划日相同测量时刻点的发电值,建立一个调整计划日发电机有功功率的优化模型:
s.t.PG-PD-B′θ=0
上式中,PGi为电网中各发电机组的有功功率;为计划发电机组安排的计划发电值,为非计划发电机组设定的计划发电值,N1为计划发电机组数目,N2为所有发电机组数目,PTie,k为电网各联络断面传输的有功功率,为电网各联络断面计划传输的有功功率,wi、wi’和wk分别为上述优化模型目标函数的三个部分的加权系数,PG为上述优化模型直流潮流方程组中电网发电机的有功功率,PG=[PG1,PG2,…,PGN2],PD为上述优化模型直流潮流方程组中电网母线负荷预测的有功功率,PD=[Pd1,Pd2,…,PdNd],Nd为电网负荷个数,B’为上述优化模型中的直流潮流方程组导纳矩阵,θm和θn分别为电网中的支路(m,n)首端节点m和末端节点n的电压相角,xmn为电网中的支路(m,n)的电抗;Pmn为电网中的支路(m,n)的有功功率,和分别为电网中的支路(m,n)有功功率的上限值和下限值;PTie,k、和分别为电网联络断面k传输的有功功率、有功功率上限值和有功功率下限值;S为电网的联络断面集;
上述优化模型中的目标函数为三个部分的加权和,三个部分分别为:计划发电机组的有功功率与该发电机组的计划发电值的偏差的平方、非计划发电机组的有功功率与该发电机组的计划发电值的偏差的平方,以及电网中各联络断面的有功功率与联络断面有功功率传输计划值的偏差的平方,前后三个加权系数的比例为:
wi∶wi’∶wk=(10~20)∶1∶(25~100);
上述优化模型的控制变量为电网中发电机组的有功功率;
上述优化模型的约束条件为:
电网的直流潮流方程组,即发电机有功总功率=负荷有功总功率+电网有功功率损耗,
电网中各发电机的有功功率上限值和有功功率下限值,
电网中各输电线路的有功功率极限值以及电网中各变压器绕组的有功功率极限值,
电网中各联络断面的有功功率极限值;
求解上述优化模型,根据计划日的发电机有功功率计划值,得到计划日电网中发电机的有功功率,且计划日发电机的有功功率与计划日的发电机计划发电值之间的偏差最小;
(3)建立一个用于求解计划日电网无功电压分布的优化模型,以得到与上述步骤(2)的计划日电网中发电机的有功功率相匹配的计划日电网中发电机的无功功率和发电机机端电压:
min PGSlack
s.t.f(V,θ)=0
上式中,PGSlack为电网中容量最大的非计划发电机组的有功功率,f(V,θ)=0为电网中节点的交流潮流方程组;Vj为电网内节点j的电压幅值,和分别为电网内节点j的电压幅值的上限值和下限值;PGi为电网内发电机组i的有功功率,和分别为电网内发电机组i的有功功率上限值和下限值;QGi为电网内发电机无功功率,和分别为发电机组i的无功功率的上限值和下限值;Pmn为电网中节点m和节点n之间的支路(m,n)传输的有功功率,及分别为电网中的支路(m,n)传输的有功功率上限值和下限值;Qzonel为电网中各设定管理区域的无功备用,和分别为电网中各设定管理区域无功备用的上限值和下限值;
上述优化模型的目标函数是电网有功损耗最小;
上述优化模型的控制变量为所有发电机组的无功功率、变压器分接头档位、并联电容器及并联电抗器的投退;
上述模型的状态变量为全网节点电压幅值Vi及相角θi;
上述优化模型的约束条件有:
电网中各节点的交流潮流方程,
电网中各支路的有功功率约束,
电网中各节点电压幅值约束,
电网中各发电机的无功功率约束,
电网中各设定管理区域的发电机无功备用约束;
求解上述优化模型,得到电网中各发电机机端电压和发电机无功功率,且电网有功损耗最小;
(4)建立一个重新调整电网联络断面的有功功率优化模型:
s.t.dPG-B′dθ=0
上式中,dPGi为再次进行机组有功功率调整时的调整量,dPTie,k为电网各联络断面有功功率的调整量,大于零代表上调,小于零代表下调,为步骤(3)中得到的各联络断面有功功率与有功传输计划值的偏差,dPG为基于直流潮流的发电机组有功功率调整量的向量,dPG=[dPG1,dPG2,…,dPGN2],dθ为基于直流潮流的电网节点相角调整量的向量,dθ=[dθ1,dθ2,…,dθn];
上述优化模型为基于直流潮流的最优增量调节模型,目标函数包括三项:计划机组的有功功率再调整量,非计划机组有功功率再调整量,电网联络断面传输有功的再调整量;
控制变量为电网中各发电机组的有功功率再调整量;
等式约束为基于直流潮流的功率平衡方程;
求解上述优化模型,得到电网中各发电机组的有功功率再调整量;
(5)根据步骤(2)得到的电网发电机的有功功率、步骤(3)得到的计划日无功电压分布以及步骤(4)得到的电网中各发电机组的有功功率再调整量,得到电网中日前计划潮流结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210084751.7A CN102623988B (zh) | 2012-03-27 | 2012-03-27 | 一种电网中日前计划潮流的自动生成方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210084751.7A CN102623988B (zh) | 2012-03-27 | 2012-03-27 | 一种电网中日前计划潮流的自动生成方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102623988A true CN102623988A (zh) | 2012-08-01 |
CN102623988B CN102623988B (zh) | 2014-01-29 |
Family
ID=46563731
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210084751.7A Active CN102623988B (zh) | 2012-03-27 | 2012-03-27 | 一种电网中日前计划潮流的自动生成方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102623988B (zh) |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103077480A (zh) * | 2013-01-07 | 2013-05-01 | 中国电力科学研究院 | 一种电力***的安全校核方法 |
CN103178551A (zh) * | 2013-01-09 | 2013-06-26 | 中国海洋石油总公司 | 海上油田群电网的功率优化控制方法 |
CN103280797A (zh) * | 2013-04-27 | 2013-09-04 | 国家电网公司 | 一种日前静态安全校核方法 |
CN104218569A (zh) * | 2014-08-15 | 2014-12-17 | 国家电网公司 | 一种大规模电网静态安全校核的评价分析方法 |
CN105024379A (zh) * | 2015-08-07 | 2015-11-04 | 内蒙古东润能源科技有限公司 | 一种风机检修管理调度*** |
CN105044517A (zh) * | 2015-08-05 | 2015-11-11 | 河北省电力建设调整试验所 | 基于潮流特征值的智能变电站保护向量诊断方法 |
CN105071385A (zh) * | 2015-08-10 | 2015-11-18 | 国家电网公司 | 一种电网运行数据实时分析*** |
CN105356447A (zh) * | 2015-10-27 | 2016-02-24 | 河海大学 | 一种电-气互联综合能源***静态安全域分析方法 |
CN105518963A (zh) * | 2014-11-18 | 2016-04-20 | 深圳大学 | 交流电力网的发电和用电有功实时调度的对称方法 |
CN107196343A (zh) * | 2017-07-19 | 2017-09-22 | 国家电网公司 | 一种多端柔性直流孤岛电网送端的日前电压计划生成方法 |
CN108599168A (zh) * | 2018-03-30 | 2018-09-28 | 中国电力科学研究院有限公司 | 一种用于对大电网计划潮流进行合理性评估的方法及*** |
CN109980639A (zh) * | 2019-04-04 | 2019-07-05 | 重庆大学 | 互联电网联络线功率可行域的快速确定方法 |
CN110994587A (zh) * | 2019-10-28 | 2020-04-10 | 国网新疆电力有限公司 | 一种电网日前发电计划的安全校核方法 |
CN111799802A (zh) * | 2020-07-22 | 2020-10-20 | 国网重庆市电力公司 | 基于线性组合的潮流方程线性化方法 |
CN108695905B (zh) * | 2018-04-17 | 2021-04-06 | 华南理工大学 | 一种含110kV等级电网安全约束优化调度方法 |
CN112865107A (zh) * | 2021-01-11 | 2021-05-28 | 国网山西省电力公司忻州供电公司 | 一种基于负荷预测的地区电网电压无功日前预报方法 |
CN113872238A (zh) * | 2021-09-26 | 2021-12-31 | 国网江苏省电力有限公司 | 一种电力***自动电压控制方法、装置、电子设备及存储介质 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030040846A1 (en) * | 2001-05-21 | 2003-02-27 | Christian Rehtanz | Stability prediction for an electric power network |
CN102082434A (zh) * | 2011-03-02 | 2011-06-01 | 国电南瑞科技股份有限公司 | 多目标断面潮流控制的最优决策及性能评估方法 |
-
2012
- 2012-03-27 CN CN201210084751.7A patent/CN102623988B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030040846A1 (en) * | 2001-05-21 | 2003-02-27 | Christian Rehtanz | Stability prediction for an electric power network |
CN102082434A (zh) * | 2011-03-02 | 2011-06-01 | 国电南瑞科技股份有限公司 | 多目标断面潮流控制的最优决策及性能评估方法 |
Non-Patent Citations (1)
Title |
---|
郭烨等: "潮流预报在电网安全预警***中的应用", 《电力***自动化》 * |
Cited By (32)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103077480B (zh) * | 2013-01-07 | 2014-03-12 | 中国电力科学研究院 | 一种电力***的安全校核方法 |
CN103077480A (zh) * | 2013-01-07 | 2013-05-01 | 中国电力科学研究院 | 一种电力***的安全校核方法 |
CN103178551A (zh) * | 2013-01-09 | 2013-06-26 | 中国海洋石油总公司 | 海上油田群电网的功率优化控制方法 |
CN103280797A (zh) * | 2013-04-27 | 2013-09-04 | 国家电网公司 | 一种日前静态安全校核方法 |
CN103280797B (zh) * | 2013-04-27 | 2015-05-13 | 国家电网公司 | 一种日前静态安全校核方法 |
CN104218569A (zh) * | 2014-08-15 | 2014-12-17 | 国家电网公司 | 一种大规模电网静态安全校核的评价分析方法 |
CN104218569B (zh) * | 2014-08-15 | 2016-02-24 | 国家电网公司 | 一种大规模电网静态安全校核的评价分析方法 |
CN105518963B (zh) * | 2014-11-18 | 2018-02-02 | 深圳大学 | 交流电力网的发电和用电有功实时调度的对称方法 |
CN105518963A (zh) * | 2014-11-18 | 2016-04-20 | 深圳大学 | 交流电力网的发电和用电有功实时调度的对称方法 |
CN105044517B (zh) * | 2015-08-05 | 2018-11-30 | 河北省电力建设调整试验所 | 基于潮流特征值的智能变电站保护向量诊断方法 |
CN105044517A (zh) * | 2015-08-05 | 2015-11-11 | 河北省电力建设调整试验所 | 基于潮流特征值的智能变电站保护向量诊断方法 |
CN105024379A (zh) * | 2015-08-07 | 2015-11-04 | 内蒙古东润能源科技有限公司 | 一种风机检修管理调度*** |
CN105024379B (zh) * | 2015-08-07 | 2018-11-06 | 内蒙古东润能源科技有限公司 | 一种风机检修管理调度*** |
CN105071385A (zh) * | 2015-08-10 | 2015-11-18 | 国家电网公司 | 一种电网运行数据实时分析*** |
CN105071385B (zh) * | 2015-08-10 | 2019-02-12 | 国家电网公司 | 一种电网运行数据实时分析*** |
CN105356447B (zh) * | 2015-10-27 | 2018-01-19 | 河海大学 | 一种电‑气互联综合能源***静态安全域分析方法 |
CN105356447A (zh) * | 2015-10-27 | 2016-02-24 | 河海大学 | 一种电-气互联综合能源***静态安全域分析方法 |
CN107196343A (zh) * | 2017-07-19 | 2017-09-22 | 国家电网公司 | 一种多端柔性直流孤岛电网送端的日前电压计划生成方法 |
CN107196343B (zh) * | 2017-07-19 | 2020-01-10 | 国家电网公司 | 一种多端柔性直流孤岛电网送端的日前电压计划生成方法 |
CN108599168A (zh) * | 2018-03-30 | 2018-09-28 | 中国电力科学研究院有限公司 | 一种用于对大电网计划潮流进行合理性评估的方法及*** |
CN108599168B (zh) * | 2018-03-30 | 2020-12-04 | 中国电力科学研究院有限公司 | 一种用于对大电网计划潮流进行合理性评估的方法及*** |
CN108695905B (zh) * | 2018-04-17 | 2021-04-06 | 华南理工大学 | 一种含110kV等级电网安全约束优化调度方法 |
CN109980639B (zh) * | 2019-04-04 | 2020-09-01 | 重庆大学 | 互联电网联络线功率可行域的快速确定方法 |
CN109980639A (zh) * | 2019-04-04 | 2019-07-05 | 重庆大学 | 互联电网联络线功率可行域的快速确定方法 |
CN110994587A (zh) * | 2019-10-28 | 2020-04-10 | 国网新疆电力有限公司 | 一种电网日前发电计划的安全校核方法 |
CN110994587B (zh) * | 2019-10-28 | 2023-02-28 | 国网新疆电力有限公司 | 一种电网日前发电计划的安全校核方法 |
CN111799802A (zh) * | 2020-07-22 | 2020-10-20 | 国网重庆市电力公司 | 基于线性组合的潮流方程线性化方法 |
CN111799802B (zh) * | 2020-07-22 | 2023-12-05 | 国网重庆市电力公司 | 基于线性组合的潮流方程线性化方法 |
CN112865107A (zh) * | 2021-01-11 | 2021-05-28 | 国网山西省电力公司忻州供电公司 | 一种基于负荷预测的地区电网电压无功日前预报方法 |
CN112865107B (zh) * | 2021-01-11 | 2022-08-30 | 国网山西省电力公司忻州供电公司 | 一种基于负荷预测的地区电网电压无功日前预报方法 |
CN113872238A (zh) * | 2021-09-26 | 2021-12-31 | 国网江苏省电力有限公司 | 一种电力***自动电压控制方法、装置、电子设备及存储介质 |
CN113872238B (zh) * | 2021-09-26 | 2024-01-30 | 国网江苏省电力有限公司 | 一种电力***自动电压控制方法、装置、电子设备及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN102623988B (zh) | 2014-01-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102623988B (zh) | 一种电网中日前计划潮流的自动生成方法 | |
CN102709918B (zh) | 基于avc***的电网关口无功功率控制方法 | |
CN103077480B (zh) | 一种电力***的安全校核方法 | |
Qin | Voltage control in the future power transmission systems | |
CN104158199B (zh) | 对电力***实时状态进行无功电压优化控制的***和方法 | |
CN104810853A (zh) | 基于有向图的含分布式发电配电网孤岛划分方法 | |
Nentwig et al. | Application of DC choppers in HVDC grids | |
CN105391094A (zh) | 孤立微电网分布式控制方法 | |
Li et al. | Research on power flow calculation method of true bipolar VSC-HVDC grids with different operation modes and control strategies | |
CN109494750B (zh) | 一种高中压配电网的分层分布式电压优化控制方法 | |
CN109888791A (zh) | 一种基于混合式配电变压器的主动配电网及其控制方法 | |
CN106208099A (zh) | 一种基于二层规划的电力***无功优化方法及其应用 | |
Wang et al. | LMP decomposition with three-phase DCOPF for distribution system | |
Nabatirad et al. | Decentralized secondary controller in islanded dc microgrids to enhance voltage regulation and load sharing accuracy | |
CN103730900A (zh) | 电力***多时间尺度的省地县一体化无功优化方法 | |
Gebbran et al. | Multiperiod DER coordination using ADMM-based three-block distributed AC optimal power flow considering inverter volt-var control | |
CN106655150A (zh) | 一种配电网运行方式优化方法 | |
Zolfaghari et al. | Adaptive gain-regulating-based control of parallel-connected γ-Z-source power converters in hybrid microgrids | |
Zhang et al. | Successive MISOCP algorithm for islanded distribution networks with soft open points | |
CN104659812A (zh) | 一种基于预测控制的多微网协调控制方法 | |
CN106877501A (zh) | 基于控保一体化概念的直流变电站层级运行***及方法 | |
Zhang et al. | Day-ahead stochastic optimal dispatch of LCC-HVDC interconnected power system considering flexibility improvement measures of sending system | |
CN106026158A (zh) | 基于最优潮流的二层规划模型的交直流混合电网优化方法 | |
CN110601212B (zh) | 一种基于无功优化配置确定电网主变变比的方法及*** | |
CN105096209A (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 |