1、 (1)处理得到简化的修方程式 (2)在潮流计算中,(2)式中两个修正方程式交替迭代,Stott把在此基础发展起来的潮流算法称为快速分解法。快速分解法是目前电力系统进行潮流计算的主要方法。由于该方法将节点电压的相位和幅值的迭代过程分开进行,使修正方程的系数矩阵维数降低很多,系数矩阵为常数且对称,因而快速分解法具有计算速度快,占用内存少等特点。2 理论基础以定雅克比牛顿-拉夫逊迭代方程为出发点,具体过程如下:通过高斯消去法,把牛顿-拉夫逊法的每一次迭代等价地细分为三步计算;对每一步运算作详细分析,证明了在连续的两次牛顿-拉夫逊迭代中,上一次迭代的第三步和下一次迭代的第一步可以合并,从而导出等效的
2、两步式分解算法;论证了该两步式分解算法的系数矩阵与快速分解算法的系数矩阵是一致的。推导过程并未因用任何解耦的假设。将极坐标型定雅克比法的修正公式中和用和代替,用代替。 (3)其中2.1将原问题分解为,子问题用高斯消去法消去子块 (4)记因,得到于是定义解为 (5)在给定的电压幅值和相角初值附近,保持电压相角不变,考虑只有电压幅值的变化时,有功功率的偏差量为 (6)综合上述结果,如果当前的迭代点为,则第次迭代可以分为三步 (7) (8) (9)2.2简化无功迭代步骤按式(7)、式(8)和式(9)完成第次迭代后,下面考虑第+1次迭代,有 (10)利用式(9),上式中的无功功率偏差为 (11)代入式
3、(10)中,经整理得 (12)式(12)说明,如果将第次迭代计算出的和用于计算第迭代的无功功率偏差量,则所求得的第次迭代的电压修正量将自动包含第次迭代计算出的。所以,的计算可以省略,相当于将式(9)与式(10)合并,因此,第次迭代可以两步完成 (13) (14)整理式(13)和式(14)得到简化的修正方程式 (15)因和即为极坐标型定雅克比法的修正公式中和,所以式(15)与式(2)格式相同。与XB型快速分解法的修正公式相比,系数矩阵是导纳矩阵的虚部,与相同。由的定义,有 (16)对于一般的电网,可能有较复杂的结构。为了对有直观的认识,假定网络中无接点,则式(16)中各矩阵的维数相等,并且接点导
4、纳矩阵可用节点支路关联矩阵和支路导纳对角矩阵(分别用的和表示电纳和电导)表示。下面将证明,对于树形电网或所有支路的比值都相同的环形网络,与相等。如果网络是树状的,其关联矩阵是方阵且非奇异,此时对式(16)有 (17)式中,为以为支路电纳组成的对角线矩阵;为以为支路电纳建立的节点电纳矩阵。这说明对树形电网,就是XB型快速分解法中的阵。对于环形网络,如果电网是均一网,即对任一条支路有,则得并有所以,故有 (18)如果电网不是均一网,上述结论不再严格成立。但和相比,在的零元素处,相应的元素近似等于零;在的非零元素处,相应的元素近似和的非零元素相等。这可以用下面的例子来说明。图1 四节点电力系统以图1
5、所示的四节点系统为例,图中给出了支路阻抗。该例中,和分别为可见,比更接近于,而用代替即得到XB型快速分解法。以上推导过程中,只在对有功功率偏差量和无功功率偏差量的计算处做了线性化近似,既没有r1);%找出非PQ节点的编号 B22=B; %BB矩阵为中间变量 B22(:,npq)=;%非PQ节点编号对应行置空 B22(npq,: B22 B2=B22; B2=inv(B2);%求得B2矩阵 %-计算各节点有功功率不平衡量dPi k=0; %k为迭代次数 ep=0; %计算P不平衡量dPi的收敛标志(0表示不收敛,1表示收敛) eq=0; %计算U不平衡量dQi的收敛标志(0表示不收敛,1表示收敛) np=find(TYPE(:,1)3);%找出非平衡节点编号 dPi=zeros(x-1,1);%dPi为x*1阶矩阵,x即为节点数 pq=find(TYPE(:,1)=1);%找出PQ节点编号 pqn=size(B2); pqn=pqn(1);%求PQ节点的个数(因B1矩阵的维数等于PQ节点数) dQi=zeros(pqn,1);%dQi为pqn*1阶矩阵while(eq|ep)&(k100) k=k+1;
copyright@ 2008-2022 冰豆网网站版权所有
经营许可证编号:鄂ICP备2022015515号-1