基于动态连续潮流的电力系统静态电压稳定的评估方法

1、一种基于动态连续潮流的电力系统静态电压稳定的评估方法，其特征在于该方法包括以下步骤： (1)构建电力系统的动态有功潮流方程为PGi-βiΔPΣ-PDi-Pi(U，θ)＝0， 其中，i＝1，2，...，N，PGi、PDi分别为电力系统节点i上当前有功发电输出功率和有功负荷，ΔPΣ为电力系统的功率不平衡量，βi为电力系统节点i上发电机分担的不平衡功率份 N额，βi满足 ∑βi＝1，若节点i没有接发电机，或该节点发电机输出功率不可调，则βi＝0， i＝1U为电力系统节点电压幅值，θ为节点电压的相角； (2)根据上述电力系统的动态有功潮流方程，构建一个扩展动态连续潮流方程： f′(U，θ，λ)＝F′(U，θ)+λb＝0 其中，F′(U，θ)为电力系统的上述动态有功潮流方程和常规无功潮流方程组成的潮流方程，λ为表征电力系统中各节点功率变化的参数，b为表征电力系统中各节点功率注入方向和大小的向量； (3)采用常规的连续潮流预测校正法，对上述扩展动态连续潮流方程进行求解，得到电力系统负荷变化时的稳定裕度。
2、如权利要求1所述的方法，其特征在于其中所述的构建动态有功潮流方程的过程包括以下步骤：(1) 设电力系统的有功潮流方程为：formula see original document page 2(2) 设节点i上发电机承担的不平衡功率的份额为A，根据电力系统的功率不平衡里：formula see original document page 2得到动态连续潮流方程：其中A。M (^， 0是电力系统总网损。

基于动态连续潮流的电力系统静态电压稳定的评估方法
技术领域
本发明涉及一种基于动态连续潮流的电力系统静态电压稳定的评估方法，属于电力系 统运行和控制技术领域。
背景技术
静态电压稳定分析是电网自动电压控制系统（AVC)和能量管理系统（EMS)中的重要 功能，对电网的安全运行意义重大。连续潮流法(CPF， Continuation Power Flow)是电力 系统静态电压稳定性分析中的一种基本分析方法。连续潮流法由于扩展的参数方程的引 入，使得其雅可比矩阵在原潮流方程鞍结型分岔点处非奇异，从而克服了普通潮流方程在 鞍结型分岔点附近不收敛的缺点。按照扩展的参数方程不同，连续潮流模型可以分为负荷 型连续潮流、支路型连续潮流、控制型连续潮流、三相不对称连续潮流。负荷型连续潮流 是连续潮流法的典型模型，连续潮流问题的最初提出就是为了研究和模拟系统中单个节 点、多个节点、某区域或全网负荷(和发电)变化对于系统的非线性影响。
一个典型的负荷型连续潮流方法，其具体实施步骤为：
1、构造扩展参数化潮流方程：引入负荷参数A，参数化后的负荷和发电可以表示为-尸C尸"0"《,' 'S〜
2"义）-2".,。 + "^.〃 'e〜 （1) 尸"(义）-L。 + "户"'/e〜 式中&"， ^/，'.。， />"'。分别为初始运行点的负荷''的有功、无功和发电机J'的有功， A々'，Ag々'，^/^分别为预定的负荷''的有功、无功增加量和发电机/的有功增加量，nrf
和G《分别为参与负荷集合和参与发电机集合。
常规潮流方程在引入参数A后的参数化潮流方程可表示成，
/(JC，A)-F(:c) +义"0 (2)
式中：引入负荷参数A, ^为系统静态状态向量(即节点电压幅值和相角）；6为表示
节点功率注入方向和大小的向量：^"对应常规潮流方程。
根据不同的参数化方法可以构造不同的扩展参数化潮流方程，现有的参数化方法有局 部参数化、弧长参数化、拟弧长参数化、正交参数化等，以局部参数化方法为例，引入扩 展方程
formula see original document page 4 (3)
式中：^为计算步长，在计算中是已知量：^为^的分量，下标^的取法为
formula see original document page 4 ④
式中：i一2，…'^为变量^&'''，、的梯度。 由此得到的扩展潮流方程为
formula see original document page 4 (5)
2、预测环节：根据当前点及其以往几点来给出解轨迹上下一个点的估计值，从而有 利于下一点求解的快速收敛。连续潮流中通常釆用的预测方法有一阶微分方法(如正切预 测法)和多项式外插方法(如二分法等）。以正切预测法为例，首先通过指定切向量中与控 制参数相对应的分量为1或-l来计算切向量，如下式所示。
formula see original document page 4 (6)
式中^只有与控制参数相对应的第k个分量为l，其余为0。下一个计算步的控制参 数取切向量中变化最剧烈的分量^^ (绝对值最大者)，上式中对应扩展方程右端选+1或 -l由其控制分量对应&4的符号决定。解向量的预测值可由下式计算
formula see original document page 4 (7)
式中^为步长。
3、校正环节：校正环节就是以预测值为初值计算扩展潮流方程，显然，预测值越靠 近解点，方程的收敛性就越好。校正解的计算可通过下式求解
formula see original document page 4 (8)
7是预测解中对应控制参数的第^个分量，以巧，^)为初值，可釆用牛顿法或拟牛顿 法求解上述方程。
4、 步长控制，.理想的步长控制方法应能够随曲线形状的变化而调整，当曲线在拐点 附近时，由于预测解可能超出了下一个计算点的收敛域，而使连续型计算方法无法正常通
过拐点进入p「曲线的下半支，所以，必须采取一定的措施使程序在拐点附近自适应的减
小步长，如采用试探-回退的方法，当校正环节判断计算发散时，退回上一计算点，将步 长减半，继续计算，这样持续沿控制方向进行试探和减小步长，直到校正计算收敛为止。
5、 回到第2步，直到求解完解曲线上所有希望求取的点，结束计算。 由于常规连续潮流在处理系统负荷增量时，要满足下式，
Zapg,2ap。,;i (9)
事实上，随着系统负荷的增长，系统中的网损也在增加，当系统重负荷时，祠损的增 量尤其显著，在常规连续潮流计算中，这部分有功不平衡量将完全由平衡节点承担。网损 增量完全由平衡节点承担的处理方法，会造成在连续潮流的计算过程中， 一方面，功率传 输方向随着网损的增大不断地发生改变，在重负荷情况下，由于网损增加显著，功率传输方向也将与初始的功率传输方向产生较大的偏差；另一方面，功率传输方向将因平衡节点 的不同而发生改变，当平衡节点位于负荷增长中心或平衡节点远离负荷增长中心时，这种 变化尤为明显。功率传输方向的变化， 一方面，使得计算出的负荷裕度并不是初始设定功 率传输方向上的负荷裕度，另一方面，使得选择不同平衡节点时，计算出的负荷裕度存在 较大的差异。因此不能对电力系统的电压稳定性进行精确评估。
发明内容
本发明的目的是提出一种基于动态连续潮流（Dynamic Continuation Power Flow, 以下简称DCPF)的电力系统静态电压稳定的评估方法，自动合理地在全网发电机上分配 电力系统中的不平衡功率，以解决常规连续潮流计算结果依赖于平衡节点选择的问题。
本发明提出的基于动态连续潮流的电力系统静态电压稳定的评估方法，包括以下步
骤：
(1)构建电力系统的动态有功潮流方程为尸c7,—AA尸s-尸d,—《("0 = 0 , 其中，z' = l，2，...,W，尸g,、 /^分别为电力系统节点t'上当前有功发电输出功率和有功负
荷，A为电力系统节点/上发电机分担的不平衡功率份额，若节点/没有接发电机，或该
节点发电机输出功率不可调，则A-o, ^满足^《=1， c/为电力系统节点电压幅值，e
为节点电压的相角；
(2) 根据上述电力系统的动态有功潮流方程，构建一个扩展动态连续潮流方程-
/'(f/，e，义）-r(f/，6o+^-o 其中，r((/，e〕为电力系统的上述动态有功潮流方程和常规无功潮流方程组成的潮流方
程，A为表征电力系统中各节点功率变化的参数，A为表征电力系统中各节点功率注入方 向和大小的向量；
(3) 釆用常规的连续潮流预测校正法，对上述扩展动态连续潮流方程进行求解，得 到电力系统负荷变化时的稳定裕度。
上述方法中，构建动态有功潮流方程的过程包括以下步骤： '
(1) 设电力系统的有功潮流方程为-
H，，e)-O, !'-l，2，…,W
(2) 设节点i上发电机承担的不平衡功率的份额为A,根据电力系统的功率不平衡
量：
得到动态连续潮流方程：AA^-尸d,-/^/，0-0 。
本发明提出的基于动态连续潮流的电力系统静态电压稳定的评估方法，是在已有的连 续潮流方法基础上的一种改进，在已有的电力系统连续潮流方程中考虑了网损不平衡量在 全网发电机上的分配，通过求解动态连续潮流方程，合理地分配不平衡功率，解决了常规 连续潮流计算结果依赖于平衡节点选择的问题，提高了电力系统静态电压稳定评估的准确 性。
附图说明
图1为本发明方法的一个实施例中所用的IEEE 39节点系统图。
图2为本发明的一个实施例与现有的连续潮流方法的计算结果对比图。
具体实施方式
本发明提出的基于动态连续潮流的电力系统静态电压稳定的评估方法，首先构建电力
系统的动态有功潮流方程为& - AA^ = 0 ，
其中，/ = l,2，...,iV， &、 A,分别为电力系统节点/上当前有功发电输出功率和有功负
荷，/?,为电力系统节点/上发电机分担的不平衡功率份额，若节点z'没有接发电机，或该 节点发电机输出功率不可调，则A-0，《满足|>,-1， C/为电力系统节点电压幅值，0
为节点电压的相角；
根据上述电力系统的动态有功潮流方程，构建一个扩展动态连续潮流方程：
《义）- r(t/, e)+- o 其中，F'(f/,。为电力系统的上述动态有功潮流方程和常规无功潮流方程组成的潮流方 程，义为表征电力系统中各节点功率变化的参数，6为表征电力系统中各节点功率注入方
向和大小的向量；
采用常规的连续潮流预测校正法，对上述扩展动态连续潮流方程进行求解，得到电力 系统负荷变化时的稳定裕度。
上述方法中，构建动态有功潮流方程的过程为设电力系统的有功潮流方程为-H，,。-0， "l，2，…，iV
设节点i上发电机承担的不平衡功率的份额为A，根据电力系统的功率不平衡量：
得到动态连续潮流方程：AA^-&-/Kt/^) = o 。
以下详细介绍本发明方法的具体过程。
形成动态潮流方程：设系统的有功潮流方程，
HC-o, "l，2，…，iV (10) 式（10)中包括平衡节点的方程。
如果系统中出现的功率不平衡量为，
/豕1 /«1
尸C"尸C是节点t'当前的有功发电输出功率和有功负荷，A。M是系统总网损。这一差额
应由所有发电机共同分担，式（10)变成
尸G,-AdC)-0， '、1，2,…,W (12)
其中，A为节点i'上发电机分担的不平衡功率份额，若节点/没有接发电机，或该节 点发电机输出功率不可调，则/?,-0。 A满足i]A-l。
由（12)式和原来的无功潮流方程一起，就得到了动态潮流方程-
r(,-o (13)
可见，常规潮流中的F^节点的发电机有功功率也需要事先给定。整个系统的功率差 额，包括潮流计算中事先无法知道潮流计算结束时才能确定的网损将由所有发电机来平衡。
构造参数化的扩展动态潮流方程：基于式（13)，引入参数;i，可以得到参数化的动
态潮流方程：
/'(c/，e,;i)-尸'(f/，e)+^-o (14) 按照传输功率计算的定义，网损的增量由功率输出区域的发电机承担，因此，A可以 由方向向量6完全描述。这样，功率不平衡量，
式中，A&。M("。为当前网损与基态网损的变化量，即
A/^"C/") - A似("")—尸丄。,""。"。） （16)
因此，发电机节点有功方程，改写为
尸c,一6A^ —S(f/刀)-0， ,'e发电机节点 （17)
以局部参数化方法为例，引入扩展方程
;?(C/,0，义)=- = 0 (18)
式中，AS为计算步长，在计算中是已知量。^为C/的分量，下标A的取法为：
-KJ-max^l'l^l,...,!^} (19)
式中，&1，...，^为变量《1，...，"„的梯度。
得到的扩展动态潮流方程，
^ 、 " （20)
可见，与基于常规潮流的扩展潮流方程相比，未知数个数与方程个数都一致，差别仅 在于式（17)的发电机节点有功方程，需要加入对网损变化量Ai^,(C/")的分配量。
网损尸^,("。的表达式可写为，
=Zf = i",Z"A ="…，W (21)
'=>1 /»1 乂e/
式中，ye/表示所有和；相联的节点y'，包括7 = /。
式(17)为本发明的关键内容，从计算求解的角度出发，可以在此基础上进行一定程度 的化简或扩展，例如，在形成雅可比矩阵时忽略^^,0/，。的引入对雅可比矩阵的修正量， 直接使用原雅可比矩阵进行迭代求解；引入网损状态变量，再扩展一个网损方程进行求解。
采用如同常规连续潮流法一样的预测-校正方法进行求解。
以下是本发明方法的一个实施例：以IEEE39节点系统进行仿真实验作实施例，进一 步说明如下： •
1EEE39节点系统结构如图l所示，图中粗实线表示母线（即节点模型中的节点），母 线旁数字表示节点编号，向下的箭头表示负荷，符号圆中带'G'字符表示发电机。其DCPF
方法具体包括以下步骤：
1) 形成动态潮流方程； 对于本实施例，根据公式列式。
2) 构造参数化的扩展动态潮流方程：
本实施例中，釆用扩展网损方程的方法，根据公式列出参数化的扩展动态潮流方程。
3) 采用如同常规连续潮流法一样的预测-校正方法进行求解：
在本实施例中，采用如同常规连续潮流法一样的预测-校正方法进行求解。求解结果 与已有技术的比较，如图2所示。图中，横坐标是不同的平衡节点，纵坐标是计算出的全 网稳定裕度，斜线柱条是采用常规连续潮流法的计算结果，点柱条是采用动态连续潮流法 的计算结果。可见，采用常规连续潮流时，选择不同的发电机作为平衡节点，计算出的稳 定裕度差别很大：而采用改进的基于动态潮流方程的连续潮流方法，选择不同的发电机节 点作为平衡节点，所计算出的稳定裕度是一致的。因此，采用动态连续潮流方法进行计算， 稳定裕度的结果不依赖于平衡节点的选择，提高了电力系统静态电压稳定评估的准确性。