电力系统自动化 第6版 [李先彬 主编] 2014年_201-250.docxVIP

194电力系统自动化(第六版)

由于考虑了平衡节点，所以式中B′为n-1阶方阵，考虑有r个PV节点，B为n-r-1阶方阵，而△P与△Q为估计功率与每次交叉迭代功率计算值的差值。

现在利用式(6-10a)来讨论两种预想事故下的安全分析。

(1)节点i、j间的连线断开，即B;;=0。从第五章式(5-20a)的推导过程可知B′;=B?1+Bi2+…+Bin

B′;;=-B;

式中B;;——各条线路的串联导纳。

式(6-11)中右端

为n-1个节点电压相角增量列向量，而

为n-r-1个节点电压的幅值增量列向量，及

(6-12a)

现以图5-13的电力系统为例，母线4为平衡节点，引用表5-8及表5-9中的数

值得

以结点1为PV节点，则

现假定预想事故为3-4连线断开，即B??=0。由于母线4是平衡节点，所以断开3—4连线，只需计算B?3的变化，即

B3—(B′23+B′13)=+(0.4651+0.1331)=0.5982所以，在B?=0的安全分析中

第六章电力系统的安全调度与运行动态检测195

(6-12b)

可见在预想断线事故的安全分析中，P-Q分解法的方程组为

(6-12c)

在式(6-12c)中，右端B′?与B1已在式(6-12b)中表示了，按照P-Q分解法的定义，左端△P与△Q各元素的表示式为

(6-13a)

式中Pi、Qi——状态估计提供的节点i的实时有功功率与无功功率，在预想断线事故的安全分析中始终保持不变；

P、Q:—每次迭代所得的节点i的有功功率与无功功率的计算值，式(5-20a)

就是P;与Q;的计算公式。

为简化三角函数的计算手续，可令

(6-13b)

至此就可以按电力网络中求正常运行状态时的P-Q分解法的步骤，求解在预想断线事故后各节点的状态变量U;与θi,并由此算出各条连线上的潮流，具体步骤与计算数据不再赘述。

式(6-13a)的存在说明安全分析中的P-Q分解法也是一个交叉迭代求解的过程，这是与直流潮流法根本不同的地方。运用式(6-10b)不但可解出在预想事故下各联络线的潮流分布，用以估计是否有过负荷的情况，而且还能求出各节点的电压幅值，用以估计是否有过电压的情况。因此安全分析中的P-Q分解法可以提供电力系统所有状态变量的变化情况，从而使调度人员对某种预想事故的后果有一个全面、清晰的了解。如有必要，当然可以采取相应的措施，以达到安全调度的目的。

上述例子是较为简单的，仅说明了问题的有关方面。在一般的情况中，如果i-j连线断开，即B;;=0,这只能影响式(6-12a)右端两方阵中B′;、B′;;及B′;、B′;的数值。令M为i、j两元素分别为1、-1,而其余元素为0的行向量，则

并有

(6-14)

196电力系统自动化(第六版)

其中

X′=B′-1MT

X=B-1MT

C′、C”是数量，X、X”是列向量。由此得

(6-15)

其中B′?-1与B1-1可以运用式(6-14),在原有B′?-1与B1-1的基础上求到。

按

按B,=0启动

输入BB′-1

BB-1

计算C′,C,X′,X”得出BB{-1B”B”-1

输入时实时数据P?θ:U?

KP=KQ=1

按式(6-15)第一式计算△θ?θk+1=θk+△θk

否|△0maxEo?

KP=1

是KQ=1?否

按式(6-15)第二式计算△θkUk+1=U+△U

否

△0-=e

KQ=1

否

KP=1?

是

KP=0

是

KQ=0

输出

是

2

图6-8P-Q分解法安全分析程序框图

△P与△Q;按式(6-13a)计算，式(5-20a)是P、Q;的计算公式。Pis、Qis为状态估计提供的实时数据。因此按一般P-Q分解法的步骤，将式(6-15)两式进行交叉迭代，就可得出一般断线预想事故的安全分析结果。其程序框图如图6-8所示，程序框图中的符号可参看图5-15。

如果在调度计算机中，存有线性方程的三角因子分解程序(如状态估计中使用此类程序时),则安全分析可按式(