电网潮流计算

潮流计算就是根据给定电力系统的网络结构、参数和决定电力系统运行状况的边界条件，确定电力系统稳态运行状态的方法。电力系统潮流计算是研究电力系统稳态运行状况的一种计算，它根据给定的运行条件及系统接线情况确定整个电力系统的运行状态:各母线的电压，各元件中流过的功率，系统的功率损耗等等。从数学上讲，潮流计算是要求解一组由潮流方程描述的非线性代数方程组。电力系统潮流计算是电力系统分析中最基本的最重要的计算，是电力系统无功优化的前提与基础，是无功优化最基础的计算工具。电力系统潮流计算分为离线计算和在线计算两种，前者主要用于系统规划设计和安排系统的运行方式，后者则用于正在运行系统的正常监视及实时控制。本文所研究基于前一种离线分析的计算方法。潮流方程对于N个节点的电力网络(地作为参考节点不包括在内)，如果网络结构和网络元件参数己知，则流入节点i可用1nIijjjIYV表示示。式中Y是节点导纳矩阵。以极坐标的形式写出1niijjijjjIYV电力系统计算中，给定的运行变量是节点注入功率，计算用的方程如下:*iiiiPjQVINi,,2,1用极坐标表示，则有:1niiiijjijjijPjQVYV分离出实部虚部可得1cosnijiijijijjPVVY1sinnijiijijijjQVVY上式为用极坐标表示的最基本的潮流计算方程。对于N个节点的电力系统，每个节点有四个运行变量，分别为有功注入P,无功注入Q、电压模值V、电压相角。一般说来，每个节点的4个变量中给定两个，求解另两个。哪两个作为给定的变量由该节点的类型决定。对于节点类型的给定应遵循一定的原则，并不是任意指定，当节点类型选择不当时，会出现潮流不收敛，或者潮流计算结果明显偏离实际系统的运行情况，在潮流运算初始化的时候，可以按照如下情况来指定。（1）PQ节点（负荷节点）事先给定的是节点功率(P、Q），待求的是节点电压向量(V、)。通常变电所母线也就是负荷节点通常都指定为PQ节点，当某些发电机的出力P、Q给定时，也可作为PQ节点。显然，在潮流计算中，系统大部分节点属于PQ节点。（2）PV节点（电压控制节点）给出的参数是节点的有功功率P及电压幅值V，待求量就是该节点的无功功率Q及电压向量的相角。通常选择有一定无功功率贮备的发电机母线或者有无功补偿设备的变电所母线作PV节点。PV节点上的发电机称之为PV机（或PV给定型发电机）。（3）平衡节点给定的运行参数是V和，，而待求量是该节点的P、Q，因此又称为V节点。在潮流计算中，这类节点一般只设一个。关于平衡节点的选择，一般选择系统中担任调频调压的某一发电厂（或发电机），或者将外网作为等值机来处理，作为功率平衡的调节节点，必须具备较大的有功和无功调节能力。、对每个负荷节点我们得到两个等式即1.1和1.2电压控制节点有一个等式，即1.2。于是我们得到两个非线性方程组或7-=0iiiijjjiPjQUYU。和方程组11cossinnijiijijijjnijiijijijjPVVYQVVY抽象成下列非线性方程组11221212,,,0,,,0,,,0nnnnfxxxfxxxfxxx（1.1）其中1,2,,ifin为给定在n维欧氏空间的nR的中的区域D上的实值函数，引进向量记号，令112233()0()0()(),,0,00()nnfxxfxxfxFxxxxfx则1.1可写成()0Fx（1.2）若存在*xD，使得*()0Fx，则称*x是方程1.2的解。用迭代法求解1.2，先将2.1化为等价的方程：xGx1.3这里映像:nnGDRR，例如可取GxxBFx，其中nBLR为非奇异矩阵。显然1.3和1。2等价。高斯赛德而迭代法：如果kx是变量x的厨师估计值，于是迭代格式变为1kkxGx柯西收敛定理：数列{xn}有极限的充要条件是：对任意给定的ε0，有一正整数N，当m,nN时，有|xn−xm|ε成立。又定理柯西收敛定理有，当连续迭代结果的差小于某一特定值时，就得到方程的解。+1kkxx这里是要求的精度。定理1.1设*x是方程1.2的解，:nnGDRR，若存在一个开球**,,0SSxxxxD和常数0,1a，使得对一切xS，有**-GxGxaxx，则对任意0xS，*x是迭代序列1,0,1,2,,kkxGxk的吸引点。看潮流方程7-=0iiiijjjiPjQUYU。P，Q已知，给定初值iU，jU进行迭代就能求出所要的值。高斯迭代法需要多次迭代才能获得所要的精度，并且不能保证收敛。牛顿法通过逐步线性化可以构造牛顿程序考虑方程()0Fx，这里我们假定映像nnFDRR：与开凸集D中二次G-可导，且''F于D连续。设*xD是方程组的解。0x为*x的初始近似，0xD，利用泰勒公式，我们有120000000'FxFxFxxxFxtxxxxdt一般0x充分接近*x，略去高阶无穷小量，因而可用线性方程组000'0FxFxxx近似代替方程。设方程的解为1x,则1'1000xxFxFx一般的我们有1'1,1,2,3.kkkkxxFxFxk实际计算可采用下列形式：1'0kkkkkkxxxFxxFx从上式可以看出牛顿法每步都要解一个n阶线性方程组。这个方程组的解kx可看成对前次近似kx的修正量，即kx加上修正量x就是新的近似1kx。牛顿法的局部收敛性定理设:nnFDRR在*x地开邻域0SD上F可导且*'Fx非奇异，*x为方程0Fx的解。那么存在闭球*0,SSxS，使映像1',kkkGxxFxFx在S上有定义，且牛顿迭代法产生的序列kx超线性收敛到*x。又若假定'**'FxFxaxx其中0a为常数，则牛顿迭代法至少二阶收敛。将1.1和1.2在初始估计值处进行泰勒级展开并忽略高阶项，可以得到下列一组线性方程套用1'0kkkkkkxxxFxxFx有2222222222222222222kKKKnKKKKKnnnnKnnnKKKKKnnKnKKKKnnnnnnPPPPVVPPPPPVVPQQQQQVVQQQQQVV22nKKnVV\简写为1234JJPVJJQ由节点导纳矩阵可算出节点的有功功率和无功功率，在分别对i和iV求偏导数渴求的雅克比矩阵各元素，利用公式和设定的电压相角初值计算,PQ，,PQ为计划值于计算值的的差叫做功率余额，kishkiiikishkiiiPPPQQQ解得新的节点电压和相角新的节点电压新的估计值为11kkkiiikkkiiiVVV，代入新的电压相角再迭代，直到满足收敛条件。计算步骤如下，首先要输入网络的原始数据以及各节点的给定值并形成节点导纳矩阵。输入节点电压初值(0)iV和(0)i置迭代计数k=0。然后开始进入牛顿法的迭代过程。在进入第k+1次迭代时其计算步骤如下（1）按上一次迭代算出的个节点电压，利用式(2.41)计算各类节点的不平衡量()()2(),,KKKiiiPQV；（2）按照条件()()2()1max,,KKKiiiPQV校验收敛；（3）计算雅克比矩阵的各元素；（4）解修正方程式求各节点的修正量()(),KKiiV；（5）修正各节点电压(1)()()(1)()(),KKKKKKiiiiiiVVV；（6）判断是否满足收敛条件。若满足，则停止计算；反之，令k=k+1，进行新一轮迭代。牛顿法是局部收敛的，也就是说初值选择必须接近真解*x附近，选择初值就成了牛顿法的弱点，估计牛顿法的初值可用康托罗维奇定理定理：设:nnFDRR及初始近似值0x满足下列条件，：11'0Fx存在且1'0Fx，2在0x地邻域0,Sx内，'0Fx存在冰满足李普希兹条件：'**'FxFxxx,0,,xySx。并且12112则方程在0,Sx内有解*x，且牛顿程序产生的序列nx收敛于*x，并有估计式121*2120kkjjxx其中112112当牛顿法用于潮流计算时，选择初值变得很容易，在正常运行方式下，其电压标幺值的幅值范围在1附近，节点电压相角在0附近，在运用上述定理，就可以得到更加准确地初值，使牛顿法收敛的更快。牛顿法的改进1简化牛顿法牛顿法虽然收敛速度快但是每次迭代都要计算一次矩阵的逆，同时计算2n个偏导数，计算量还是很客观的如果将'kFx改为固定的'0Fx，则只需要计算一次'kFx，迭代公式为1'10kkkxxFxFx2修正牛顿法吧m步简化牛顿法组成一次牛顿步，于是得到如下迭代程序0111'1kkiiikkkkmkkxxxxFxFxxx这种方法收敛阶为M+1阶，当M=1时为普通牛顿法，当2M时这种方法收敛阶更高，收敛速度更快但是计算量更大。带参数的牛顿法潮流方程的雅克比矩阵病态或者奇异的时候无法用普通牛顿法进行迭代，这时我们就需要带入一个参数使得雅克比矩阵变成绝对对角占优矩阵，绝对对角占优矩阵式必然有解的。此时得到迭代程序1'1kkkkkxxFxIFx，0,1,2,,k称k为阻尼因子，且,k，其中2minRe02Reiii2minRe02Reiii而i是*'Fx的特征值，Rei为i的实部。阻尼牛顿法xxap其中x为0Fx的某一近似解，并设从x出发下一次迭代近似为x，a为阻尼因子0,1a，选择参数a的原则是使之满足1FxFxapaFx这里是一个固定的数，0,1。由下式确定搜索方向P,、''FxFxpFx其中0'1下面给出阻尼牛顿法的步骤1：选初值0xD并计算0Fx，置1,1ka，选精度要求2：确定P3：确定阻尼因子a，如上一步用的阻尼因子比较小那么下一步也不太可能用到比较大的阻尼因子可用下式得到min1,2aaaa如果在前一步是减小的其他情形4计算xxap，若xx则停止计算5若xD或者1FxaFx，则2aa，并转第四步，否则1kk，转第二步。拟牛顿法类似于简单牛顿法，如果不矩阵的导数或者少计算导数，自然会大大减少计算量，减少编程难度，但是简单牛顿法也就是割线法对单个方程比较成功，对方程组较困难。若能仿照割线法的思路构造出收敛阶高而且计算量小的新算法就有重大意义。拟牛顿算法如下111111,0,1,2,,.1kkkkkkkkkkkkkxxAFxkAxxFxFxAAAran