P-Q分解法分析算法原理.docxVIP

四、P-Q分解法潮流计算 P-Q分解法潮流计算派生于极坐标表示时的牛顿-拉夫逊法。两者的区别在修正方程和计算步骤。以下仅着重讨论这两方面。 4.1潮流计算时的修正方程 PQ分解法潮流计算时的修正方程是计及电力系统的特点后对牛顿-拉夫逊法修正方程的简化。为说明这个简化，先将式(4-44)重新排列如下: （n-1） （n-1） （4-43） (m-1） （n-1） （m-1） 丨 (4-44) 重新排列时不再留空行、空列。显然，这种重新排列并不影响修正方程的内容。 对修正方程的第一个简化是，计及电力网络中各元件的电抗-般远大于电阻，?以致各节点电压相位角的改变主要影响各元件中的有功功率潮流从而各节点的注入有功功率，各节点电压大小的改变主要影响各元件中的无功功率潮流从而各节点的注人无功功率，可将式(4-44)中的子阵N、J略去，而将修正方程式简化为 (4-45) 电力网络的上述特点，不仅可从第3章中对纵、横向附加电动势的讨论，即式(3-48)中得到证实，也还可用图4-9说明。图4-9a所示为两节点电压大小相等、相位不同的情况。图4-9b所示则为两节点电压大小不等、相位相同的情况。由图可见，由于两节点间连接元件的电抗远大于电阻，其中流过的电流滞后两节点间电压降落相量dU基本达90°，以致对第一种情况，它基本为有功电流,对第二种情况，它基本为无功电流，而某节点的注人电流或注人功率则是与该节点有直接联系的所有元件电流或功率的总和。 对修正方程的第二个简化是，基于对状态变量δ,的约束条件|δi-δj|＜|δi-δj|max,，即|δi-δj|=|δij|不宜过大，计及这一条件,再计及GijBij，可以认为 cosδij≈1，Gijsinδij＜＜Bij 图4-9元件电流与节点电压的关系 a)?U1=U2,δ1≠δ2?b)?U1≠U2,δ1=δ2 于是，雅可比矩阵的各个元索可简化为 Hij=-UiUjBij Lij=-UiUjBij 再按自导纳的定义，上两式中的Ui2Bii项应为各元件电抗远大f电阻的前提下，除节点i外其他节点都接地时，由节点i注入的无功功率。这功率必远大于正常运行时节点i的注入无功功率，也即Ui2BiiQi，上两式又可化简为 Hii=-Ui2Bii Lii=-Ui2Bii 这样，雅可比矩阵中两个子阵H、L的元素将具有相同的表示式，但是它们的阶数不同，前者为(n-1)?阶、后者为(m-1)阶。 这两个子阵都可展开如下式所示: (4-46) 将式(4-46)代人式(4-45),展开，可得 (4-47a) （4-47b） 将式(4-47a)、?式(4-47b)等号左右都前乘以 可得 （4-48a） （4-48b） 它们可简写为 （4-49a） （4-49b） 这就是P-Q分解法的修正方程。式(4-49)中，等号左侧列向量中的有功、无功功率不平衡量?Pi、?Qi:?仍如式(4-39);?但等号右侧的系数矩阵B、B却并不总直接由导纳矩阵的虚数部分组成。往往为了加速收敛，在B中去除那些与有功功率、电压相位关系较小的因素，在B中去除那些与无功功率、电压相位关系较小的因素，以致这两个矩阵不仅阶数不同，而且它们相应元素的数值也不完全相等。与牛顿拉夫逊法相比，P-Q分解法的修正方程有如下的特点:1)以一个(n-1)阶和一个(m-1)阶系数矩阵B、B”替代原有的(n+m-2)阶系数矩阵J,提高计算速度，对存储容量的要求。2)以迭代过程中保持不变的系数矩阵B、B”替代起变化的系数矩阵J,显著地提高了计算速度。3)以对称的系数矩阵B、B”替代不对称的系数矩阵J,使求逆等运算量和所需的存储容量大为减少。但应强调指出，导出这修正方程时所作的种种简化毫不影响用这种方法计算的准确度。因采用这种方法时，迭代收敛的判据仍是△Pi≤ε、△Qi≤ε，而其中的△Pi、△Qi;已如上述，仍