直角坐标系下牛顿法潮流计算

武汉理工大学《电力系统分析》课程设计说明书1直角坐标系下牛顿法潮流计算1电力系统潮流计算潮流计算是电力系统分析中的一种最基本的计算，它的任务是对给定的运行条件确定系统的运行状态，如母线上的电压(幅值及相角)、网络中的功率分布及功率损耗等。在电力系统规划设计和现有电力系统运行方式的研究中，都需要利用潮流计算来定量地分析比较供电方案或运行方式的合理性.可靠性和经济性。此外，电力系统潮流计算也是计算系统动态稳定和静态稳定的基础。2节点导纳矩阵的形成在图1（a）的简单电力系统中，若略去变压器的励磁功率和线路电容，负荷用阻抗表示，便可以得到一个有5个节点（包括零电位点）和7条支路的等值网络，如图1（b）所示。将接于节点1和4的电势源和阻抗的串联组合变换成等值的电流源和导纳的并联组合，变得到图1（c）的等值网络，其中1101IyE和4404IyE分别称为节点1和4的注入电流源。1243(a)Ė1Ė41234y10y12y24y20y23y34y40(b)İ1İ41234y12y24y23y34y40y20y´10(c)图1电力系统及其网络以零电位点作为计算节点电压的参考点，根据基尔霍夫定律，可以写出4个独立节点的电流平衡方程如下：武汉理工大学《电力系统分析》课程设计说明书21011212112212022323242423323434244234434044()()()()0()()0()()yUyUUIyUUyUyUUyUUyUUyUUyUUyUUyUI（2-1）上述方程组经过整理可以写成1111221211222233244322333344422433444400YUYUIYUYUYUYUYUYUYUYUYUYUI（2-2）式中，111012Yyy；2220232412Yyyyy；332334Yyy；44402434Yyyy；122112YYy；233223YYy；244224YYy；344334YYy。一般的，对于有n个独立节点的网络，可以列写n个节点方程11112211211222221122nnnnnnnnnnYUYUYUIYUYUYUIYUYUYUI（2-3）也可以用矩阵写成1111121212222212nnnnnnnnUIYYYYYYUIYYYUI（2-4）或缩写为YUI（2-5）矩阵Y称为节点导纳矩阵。它的对角线元素iiY称为节点i的自导纳，其值等于接于节点i的所有支路导纳之和。非对角线元素ijY称为节点i、j间的互导纳，它等于直接接于节点i、j间的支路导纳的负值。若节点i、j间不存在直接支路，则有0ijY。由此可知节点导纳矩阵是一个稀疏的对称矩阵。武汉理工大学《电力系统分析》课程设计说明书33牛顿-拉夫逊法潮流计算3.1牛顿-拉夫逊法的基本原理牛顿—拉夫逊法（Newton—Raphson法）是求解非线性方程代数方程组的有效迭代计算方法。在牛顿—拉夫逊法的每一次迭代过程中，对非线性方程通过线性化处理逐步近似。下面以单变量加以说明。设有单变量非线性方程()0fx(3-1)求解此方程时。先给出解的近似值(0)x它与真解的误差为(0)x，则(0)(0)xxx将满足方程，即(0)(0)()0fxx(3-2)将(3-8)式左边的函数在(0)x附近展成泰勒级数，于是便得2'''()(0)(0)(0)(0)(0)(0)(0)(0)(0)()()()()......()....2!!()()nnffnxxfffxxxxxxx（3-3）式中'(0)()fx,……()(0)()nfx分别为函数()fx在(0)x处的一阶导数，….，n阶导数。如果差值(0)x很小，3-9式右端(0)x的二次及以上阶次的各项均可略去。于是，3-9便简化为'(0)(0)(0)(0)(0)()()()fffxxxxx＝0(3-4)这是对于变量的修正量(0)x的现行方程式，亦称修正方程式。解此方程可得修正量(0)(0)'(0)()()fxxfx(3-5)用所求的(0)x去修正近似解，变得(0)(1)(0)(0)(0)'(0)()()fxxxxxfx(3-6)武汉理工大学《电力系统分析》课程设计说明书4由于3-10是略去高次项的简化式，因此所解出的修正量(0)x也只是近似值。修正后的近似解(1)x同真解仍然有误差。但是，这样的迭代计算可以反复进行下去，迭代计算的通式是()(1)()'()()()kkkkfxxxfx(3-7)迭代过程的收敛判据为()1()kfx(3-8)或()2kx(3-9)式中1，2为预先给定的小正数。这种解法的几何意义可以从图3－1得到说明。函数y＝f(x)为图中的曲线。f(x)＝0的解相当于曲线与x轴的交点。如果第k次迭代中得到()kx，则过()()(),()kkkfyxx点作一切线，此切线同x轴的交点便确定了下一个近似值(1)kx。由此可见，牛顿－拉夫逊法实质上就是切线法，是一种逐步线性化的方法。应用牛顿法求解多变量非线性方程组3-1时，假定已给出各变量的初值1(0)x，2(0)x….(0)nx，令1(0)x，2(0)x，…..(0)nx分别为各变量的修正量，使其满足方程3-2即11122211221122(0)(0)(0)(0)(0)(0)(,,....,)0(0)(0)(0)(0)(0)(0)(,,....,)0......(0)(0)(0)(0)(0)(0)(,,....,)0nnnnnnnfxxxxxxfxxxxxxfxxxxxx(3-10)将上式中的n个多元函数在初始值附近分别展成泰勒级数,并略去含有)0(1x，)0(2x，……，)0(nx二次及以上阶次的各项，便得武汉理工大学《电力系统分析》课程设计说明书511100011212121110002121212101211(0)(0)(0)(0)(0)(0)(,,...,)...0(0)(0)(0)(0)(0)(0)(,,...,)...0......(0)(0)(0)(0)(,,...,)|||||||nnnnnnnnffffxxxxxxxxxffffxxxxxxxxxffxxxxx110022(0)(0)...0||nnffxxxx.(3-11)方程式3-17也可以写成矩阵形式11100012112222212000121200012...(0)(0)(0)(,,...,)(0)(0)(0)(,,...,).....................(0)(0)(0)(,,...,)...|||||||||nnnnnnnnnnfffxxxfxxxffffxxxxxxfxxxfffxxx12(0)(0)...(0)nxxx(3-12)方程式3-18是对于修正量)0(1x，)0(2x，……，)0(nx的线性方程组,称为牛顿法的修正方程式.利用高斯消去法或三角分解法可以解出修正量)0(1x，)0(2x，……，)0(nx。然后对初始近似值进行修正(1)(0)(0)iiixxx(i=1,2,….,n)(3-13)如此反复迭代，在进行k＋1次迭代时，从求解修正方程式11112112222212121212...()()()(,,...,)()()()(,,...,).....................()()()(,,...,)...|||||||||kkknnnkkknnnnnnkkknkkkkkkkkkfffxxxfxxxffffxxxxxxfxxxfffxxx12()()...()nkkkxxx(3-14)武汉理工大学《电力系统分析》课程设计说明书6得到修正量1()kx，2()kx，()nkx，并对各变量进行修正(1)()()iiikkkxxx(i=1,2,…,n)(3-15)式3-20和3-21也可以缩写为)()()(kkkxJxF(3-16)和)()()1(kkkxxx(3-17)式中的X和X分别是由n个变量和修正量组成的n维列向量；F(X)是由n个多元函数组成的n维列项量；J是n阶方阵，称为雅可比矩阵，它的第i、j个元素iijifJx是第n个函数12(,,...,,)nifxxx对第j个变量jx的偏导数；上角标(k)表示J阵的每一个元素都在点,,,()()()(...,)12ikkknfxxx处取值。迭代过程一直到满足收敛判据112()()()max(,,...,)inkkkfxxx(3-18)或2()maxikx(3-19)为止。1和2为预先给定的小正数。将牛顿－拉夫逊法用于潮流计算，要求将潮流方程写成形如方程式3-1的形式。由于节点电压可以采用不同的坐标系表示，牛顿－拉夫逊法潮流计算也将相应的采用不同的计算公式。武汉理工大学《电力系统分析》课程设计说明书73.2节点电压用直角坐标表示是的牛顿-拉夫逊法潮流计算采用直角坐标时，节点电压可表示为iiijfeV导纳矩阵元素则表示为ijijijjBGY将上述表示式代入nijiiiiiijjiSPjQUIUYU的右端，展开并分出实部和虚部，便得11()()nniiijjijjiijjijjjjPeGeBffGfBe（11-45）假定系统中的第1，2，3···，m号节点为PQ节点，第i个节点的给定功率设为isP和isQ，对对该节点可列写方程0)()(0)()(1111njjijjijinjjijjijiisiisinjjijjijinjjijjijiisiisieBfGffBeGePPPPeBfGffBeGePPPP（i=1,2，···，m）（11-46）假定系统中的第m+1，m+2，···，n-1号节点为PV节点，则对其中每一个节点可以列写方程0)(0)()(22222211iiisiisinjnjjijjijijijjijiisiisifeVVVVeBfGffBeGePPPP（i=m+1，m+2，···，n-1）（11-47）第n号节点为平衡点，其电压nnnjfeV是给定的，故不参加迭代。式（11-46）和式（11-47）总共包含了2（n-1）个方程，待求的变量有1111,,...,,nnfefe也是2（n-1）个。我们还可看到，方程（11-46）和式（11-47）已经具备了方程组（11-34）的形式。因此，不难写出如下的修正方程式VJW（11-48）式中11()()nniiijjijjiijjijjjjQfGeBfeGfBe武汉理工大学《电力系统分析》课程设计说明书8TnnmmmmVPVPQPQPW21121111......Tnnmmmmfe