潮流计算的基本算法及使用方法一、 潮流计算的基本算法1. 牛顿－拉夫逊法1．1 概述牛顿－拉夫逊法是目前求解非线性方程最好的一种方法;这种方法的特点就是把对非线性方程的求解过程变成反复对相应的线性方程求解的过程,通常称为逐次线性化过程,就是牛顿－拉夫逊法的核心;牛顿-拉夫逊法的基本原理是在解的某一邻域内的某一初始点出发,沿着该点的一阶偏导数——雅可比矩阵,朝减小方程的残差的方向前进一步,在新的点上再计算残差和雅可矩阵继续前进,重复这一过程直到残差达到收敛标准,即得到了非线性方程组的解;因为越靠近解,偏导数的方向越准,收敛速度也越快,所以牛顿法具有二阶收敛特性;而所谓“某一邻域”是指雅可比方向均指向解的范围,否则可能走向非线性函数的其它极值点,一般来说潮流由平电压即各母线电压相角为0,幅值为1启动即在此邻域内; 1．2 一般概念对于非线性代数方程组即 ()0,,,21=n i x x x f ()n i ,2,1= 1－1在待求量x 的某一个初始计算值()0x 附件,将上式展开泰勒级数并略去二阶及以上的高阶项,得到如下的线性化的方程组()()()()()0000=∆'+x x f x f 1－2上式称之为牛顿法的修正方程式;由此可以求得第一次迭代的修正量()()()[]()()0100x f x f x -'-=∆ 1－3将()0x ∆和()0x 相加,得到变量的第一次改进值()1x ;接着再从()1x 出发,重复上述计算过程;因此从一定的初值()0x 出发,应用牛顿法求解的迭代格式为()()()()()k k k x f x x f -=∆' 1－4()()()k k k x x x ∆+=+1 1－5上两式中：()x f '是函数()x f 对于变量x 的一阶偏导数矩阵,即雅可比矩阵J ；k 为迭代次数;由式1－4和式子1－5可见,牛顿法的核心便是反复形成求解修正方程式;牛顿法当初始估计值()0x 和方程的精确解足够接近时,收敛速度非常快,具有平方收敛特性;1．3 潮流计算的修正方程运用牛顿－拉夫逊法计算潮流分布时,首先要找出描述电力系统的非线性方程;这里仍从节点电压方程入手,设电力系统导纳矩阵已知,则系统中某节点i 节点电压方程为从而得∑=**••=nj j ij i i U Y U S 1进而有()01=-+*=*•∑j n j ij i i i U Y U jQ P1－6式1－6中,左边第一项为给定的节点注入功率,第二项为由节点电压求得的节点注入功率;他们二者之差就是节点功率的不平衡量;现在有待解决的问题就是各节点功率的不平衡量都趋近于零时,各节点电压应具有的价值;由此可见,如将式1－6作为牛顿－拉夫逊中的非线性函数()0=X F ,其中节点电压就相当于变量X ;建立了这种对应关系,就可列出修正方程式,并迭代求解;但由于节点电压可有两种表示方式——以直角做表或者极坐标表示,因而列出的迭代方程相应地也有两种,下面分别讨论;1．3．1 直角坐标表示的修正方程节点电压以直角坐标表示时,令i i i jf e U +=•、j j j jf e U +=•,且将导纳矩阵中元素表示为ij ij ij jB G Y +=,则式1－7改变为 ()()()()01=--+-+∑=nj j j ij ij i i i i jf e jB G jf e jQ P1－7再将实部和虚部分开,可得()()[]()()[]⎪⎪⎭⎪⎪⎬⎫=+---=++--∑∑==0011nj j ij j ij i j ij j ij i i nj j ij j ij i j ij j ij i i e B f G e f B e G f Q e B f G f f B e G e P 1－8这就是直角坐标下的功率方程;可见,一个节点列出了有功和无功两个方程;对于PQ 节点1,,21-=m i ，,给定量为节点注入功率,记为i P '、i Q ',则由式2－8可得功率的不平衡量,作为非线性方程()()[]()()[]⎪⎪⎭⎪⎪⎬⎫+---'=∆++--'=∆∑∑==nj j ij j ij i j ij j ij i i i nj j ij j ij i j ij j ij i i i e B f G e f B e G f Q Q e B f G f f B e G e P P 11 1－9式中i P ∆、i Q ∆——分别表示第i 节点的有功功率的不平衡量和无功功率的不平衡量;对于PV 节点n m m i ,,2,1 ++=,给定量为节点注入有功功率及电压数值,记为i P '、i U ',因此,可以利用有功功率的不平衡量和电压的不平衡量表示出非线性方程,即有()()[]()⎪⎭⎪⎬⎫+-'=∆++--'=∆∑=22221i i i i nj j ij j ij i j ij j ij i i i f e U U e B f G f f B e G e P P1－10式中i U ∆为电压的不平衡量;对于平衡节点m i =,因为电压数值及相位角给定,所以S s S jf e U +=•也确定,不需要参加迭代求节点电压;因此,对于n 个节点的系统只能列出()12-n 个方程,其中有功功率方程()1-n 个,无功功率方程()1-m 个,电压方程()m n -个;将式1－9、式1－10 非线性方程联立,称为n 个节点系统的非线性方程组,且按泰勒级数在()0i f 、()0i e m i n i ≠=,,,2,1 展开,并略去高次项,得到以矩阵形式表示的修正方程如下⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡∆∆∆∆∆∆∆∆⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡∆∆∆∆∆∆∆∆n n p p nn nnnpnpn n n n nn nn np np n n n n pn pn pp pp p p p p pn pn pp pp p p p p n n p p n n p p n n p p n n p p n n pp e f e f e f e f S R S R S R S R N H N H N H N H S R S R S R S R N H N H N H N H L J L J L J L J N H N H N H N H L J L J L J L J N H N H N H N H U P U P Q P Q P 22112211221122112211222222222121222222222121111112121111111112121111222211 1－11 上式中雅可比矩阵的各个元素则分别为 将1－11写成缩写形式[]⎥⎦⎤⎢⎣⎡∆∆=⎥⎦⎤⎢⎣⎡∆∆⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡∆∆∆e f e f S RL J N HU Q P J 2 1－12 对雅可比矩阵各元素可做如下讨论：当i j ≠时,对于特定的j ,只有该特定点的i f 和i e 是变量,于是雅可比矩阵中各非对角元素表示为当i j =时,雅可比矩阵中各对角元素的表示式为由上述表达式可知,直角坐标的雅可比矩阵有以下特点：1) 雅可比矩阵是()12-n 阶方阵,由于ji ij H H ≠、ji ij N N ≠等等,所以它是一个不对称的方阵;2) 雅可比矩阵中诸元素是节点电压的函数,在迭代过程中随电压的变化而不断地改变;3) 雅可比矩阵的非对角元素与节点导纳矩阵B Y 中对应的非对角元素有关,当B Y 中的ij Y 为零时,雅可比矩阵中相应的ij H 、ij N 、ij J 、ij L 也都为零,因此,雅可比矩阵也是一个稀疏矩阵;1．3．2 极坐标表示的修正方程在牛顿－拉夫逊计算中,选择功率方程∑=**•=-+nj j ij i i i U Y U jQ P 10作为非线性函数方程,把式中电压向量表示为极坐标形式 则节点功率方程变为 将上式分解成实部和虚部这就是功率方程的极坐标形式,由此可得到描述电力系统的非线性方程;对于PQ 节点,给定了()()⎪⎪⎭⎪⎪⎬⎫--'=+-'=∆∑∑==nj ij ij ij ij j i i i nj ij ij ij ij j i i i B G U U Q Q B G U U P P 11cos sin sin cos δδδδ ()121-=m i 、、 1－13对于PV 节点,给定了i P '、i U ',而i Q '未知,式1－13中i Q ∆将失去作用,于是PV 节点仅保留i P ∆方程,以求得电压的相位角;1－14对于平衡节点,同样因为s U 、s δ已知,不参加迭代计算;将式1－13、式1－14联立,且按泰勒级数展开,并略去高次项后,得出矩阵形式的修正方程⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡∆∆∆∆∆∆⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡∆∆∆∆∆∆n p nnnpn n n n pn pp p p p p n p n p n p n p n p U UU U H H N H N H H H N H N H L J L J L J N H N H N H L J L J L J H H N H N H P P Q P Q P δδδδ 22211122112211222121212122212121211112121111111212111122111－15雅可比矩阵终,对PV 节点,仍可写出两个方程的形式,但其中的元素以零元素代替,从而显示了雅可比矩阵的高度稀疏性;式中电压幅值的修正量采用U U ∆的形式,并没有什么特殊意义,仅是为了雅可比矩阵中各元素具有相似的表达式; 雅可比矩阵的各元素如下 将式1－15写成缩写形式⎥⎥⎦⎤⎢⎢⎣⎡∆∆⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡∆∆U U L J N HQ P δ1－16以上得到了两种坐标系下的修正方程,这是牛顿－拉夫逊潮流计算中需要反复迭代求解的基本方程式;2. 快速分解法2．1 概述快速分解法的基本思想是：把节点功率表示为电压向量的极坐标方程式,抓主要矛盾,以有功功率误差作为修正电压向量角度的依据,以无功功率误差作为修正电压幅值的依据,把有功功率和无功功率的迭代分开来进行;快速分解法根据电力系统实际运行状态的物理特点,对牛顿-拉夫逊法潮流计算的数学模型进行合理的简化;2．2 基本公式在交流高压电网中,输电线路的电抗要比电阻大得多,系统中母线有功功率的变化主要受电压相位的影响,无功功率的变化主要受母线电压幅值变化的影响;在修正方程式的系数矩阵中,偏导数δ∂∆∂Q 和V P ∂∆∂的数值相对于偏导数V Q ∂∆∂和δ∂∆∂P是相当小的,作为简化的第一步,可以将方程式2－1中的子块N 和K 略去不计,即认为它们的元素都等于零;这样,m n +-1阶的方程式便分解为一个1-n 阶和一个m 阶的方程式,即将式2－1简化为式2－2和式2－3;⎥⎦⎤⎢⎣⎡∆∆⎥⎦⎤⎢⎣⎡-=⎥⎦⎤⎢⎣⎡∆∆-V V L K N H Q P 1D δ 2－1 δ∆-=∆H P 2－2V LV Q 1D ∆-=∆- 2－3上述的简化大大地节省了计算机的内存和解题时间,但是矩阵H 和L 的元素都是节点电压幅值和相角差的函数,其数值在迭代过程中是不断变化的;因此,快速分解法潮流计算的第二个简化,也是最关键的一步简化就在于把系数矩阵H 和L 简化成在迭代过程中不变的常数对称矩阵;在一般情况下,线路两端电压的相角差是不大的通常不超过 10~ 20因此可以认为1cos ≈ij δ , ij ij ij B G <<δsin 2－4此外,与系统各节点无功功率相适应的导纳LDi B 必远小于该节点自导纳的虚部,即ii iiLDi B V Q B <<=2 或 ii i i B V Q 2<< 考虑到上面的关系,矩阵H 和L 的元素的表达式便被简化为ij j i ij B V V H = i,j=1,2,…,n -1 2－5 ij j i ij B V V L = i,j=1,2,…,m 2－6⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=------------11,1122,1111,1111,222222121211,1121211111n n n n n n n n n n n n V B V V B V V B V V B V V B V VB V V B V V B V V B VH 2－7⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=m mm m m m m m m m m m V B V V B V V B V V B V V B V V B V V B V V B V V B V221122222212121121211111L 2－8将式2－7和式2－8分别代入式2－2和2－3,便得到：用11-D V 和12-D V 分别左乘以上两式便得简化了的修正方程式,可展开写成：⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡∆∆∆⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡-=⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡∆∆∆--------=-1122111,12,11,11,222211,11211112211δδδn n n n n n n n n n V V V B B B B B B B B B V P V P V P2－9⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡∆∆∆⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡-=⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡∆∆∆m mm m m m m m m V V V B B B B B B B B B V Q V Q V Q212122221112112211 2－10 式2－9和式2－10就是快速分解法潮流计算的修正方程式,其中系数矩阵都是由节点导纳矩阵的虚部构成,只是阶次不同,矩阵B '为1-n 阶,不含平衡节点对应的行和列,矩阵B ''为m 阶,不含平衡节点和PV 节点对应的行和列;()∑=+-=-=∆nj ij ij ij ij j i is i is i B G V V P P P P 1δsin δcos 2－11∑=+-=-=∆nj ij ij ij ij j i is i is i B G V V Q Q Q Q 1)δcos δsin ( 2－12修正方程式2－9和2－10与功率误差方程式2－11和2－12构成了快速分解法迭代的基本计算公式; 2．3 快速分解法的特点快速分解法与牛顿法潮流计算的主要差别表现在它们的修正方程上;快速分解法通过对电力系统具体特点的分析,对牛顿法修正方程式的雅克比矩阵进行了有效的简化和改进,得到式2－9、式2－10所示的修正方程式;这两组方程式和牛顿法的修正方程相比主要有三个特点：a 快速分解法的修正方程式用两个n 阶线性方程组代替了一个n 2阶线方程组;b 快速分解法的修正方程式中系数矩阵的所有元素在迭代过程中维持常数不变;c 快速分解法的修正方程式中系数矩阵是对称矩阵;这些特点在提高计算速度和减少内存方面的作用是很明显的：首先,因为修正方程式的系数矩阵是导纳矩阵的虚部,因此在迭代过程中不必像牛顿法那样每次都要重新计算雅克比矩阵,这样不仅减少了运算量,而且也大大简化了程序；其次,由于系数矩阵在迭代过程中维持不变,因此在求解修正方程式时,不必每次都对系数矩阵进行消去运算,只需要在进入迭代过程以前,将系数矩阵用三角分解形成因子表,然后反复利用因子表对不同的常数项V P /∆或V Q /∆进行消去和回代运算,就可以迅速求得修正量,从而显着提高了迭代速度；第三,由于对称矩阵三角分解后,其上三角矩阵和下三角矩阵有非常简单的关系,所以在计算机中可以只存储上三角矩阵或下三角矩阵,从而也进一步节约了内存;快速分解法所采用的一系列简化假定只影响了修正方程的结构,也就是说只影响了迭代过程,但未影响最终结果;因为快速分解法和牛顿法都采用同样的数学模型,最后计算功率误差和判断收敛条件都是严格按照精确公式进行的,所以快速分解法和牛顿法一样都可以达到很高的精确度;为了改善快速分解法的收敛特性,修正方程的系数矩阵B '与B ''一般并不简单的是电力系统导纳矩阵的虚部,下面讨论一下B '与B ''的构成;B '与B ''的阶数是不同的,B '为 1-n 阶,B ''低于1-n 阶;因为式2－10不包含于PV 节点有关的项,所以,如果系统有r 个PV 节点,则B ''应为1--r n 阶;式2－9以有功功率误差为依据修正电压向量的角度,式2－10以无功功率误差依据修正电压幅值;为了加速收敛,使它们能够更有效地进行修正,可以考虑在B '中尽量去掉那些与有功功率及电压向量角度无关或影响较小的因素,而在B ''中尽量去掉与无功功率及电压幅值影响较小的因素;所以,我们以电力系统导纳矩阵的虚部作为B '和B ''时,可以在B '去掉充电电容和变压器变比的影响,在B ''中去掉输电线路电阻对B ''的影响;B '和B ''的非对角元素和对角元素可分别按式2－13和2－14计算：22ijijij ijxr x B +-=' ∑∈+='ij ijijij iixr x B 22ij x B 1-='' ∑∈-=''i j i ijb x B 01式2－13中ij r 和ij x 分别为支路ij 的电阻和感抗,式2－14中0i b 为节点i 接地支路的电纳;快速分解法改变了牛顿法迭代公式的结构,因此就改变了迭代过程的收敛性;牛顿法在迭代开始时收敛得较慢,当收敛到一定程度后,它的收敛速度非常之快,而快速分解法几乎是按同一速度收敛的,快速分解法每次迭代的计算量很小,因此快速分解法的计算速度比牛顿法有明显的提高;二、 潮流计算的使用方法1. 初始方式准备对任何潮流模拟操作计算,总是在某一个初始的运行方式上进行;这种初始方式可以是状态估计提供的实时运行方式,也可以是以往保存的历史运行方式;2. 调度操作模拟在准备好的初始潮流断面上,可以继续修改方式,模拟预想的潮流运行方式,再进行详细的潮流分析;模拟操作包括：1)开关刀闸变位模拟2)发电机功率调整3)负荷功率设置4)发电机分接头设置5)线路停运、投入6)变压器停运、投入7)母线停运、投入8)厂站停运、投入3.运行参数维护潮流计算参数画面上可以设置算法、收敛判据、迭代次数、单/多平衡机等运行参数;平衡发电机是电气岛内的电压相角参考点,当采用“单平衡机”模式时,电网的不平衡功率包括发电、负荷和网损都将由设定平衡机吸收;当采用“多平衡机”模式时,电网的不平衡功率将由多台发电机负责平衡,多台发电机之间的不平衡功率分配方式包括容量、系数和平均三种方式;选择容量时将根据发电机的可调容量分配,选择系数时根据人工设置的系数按比例分配,选择平均时则平均分配不平衡功率;在分配过程中,确保发电机的出力在最大出力和最小出力范围内;发电机参数中可以设置发电机的调节特性,包括节点类型平衡节点、PQ 节点、PV节点等,对于PV节点可以设定控制机端电压还是高压侧母线电压以及控制的目标电压值,对于按指定系数参与有功调节的机组可以设置比例系数;4.计算结果分析潮流计算结束后,计算结果分析包括：1)潮流计算状态2)电气岛、迭代信息3)潮流计算结果4)设备越限和重载监视5)运行信息5.误差统计在潮流模拟计算完成后,如果现场很快发生了模拟的动作,可以从统计每个测点模拟计算值和实际量测值相比的误差,并统计出全网平均误差,统计方法遵循实用化考核细则;在表格中全部列出所有测点的SCADA量测值、潮流模拟计算值、考核基准值以及测点误差等内容;如果只关心部分厂站的误差情况,局部误差统计中选择好需要关心的厂站,局部误差统计程序会过滤掉没有选中的厂站,只列出选中厂站的所有设备的误差统计情况;如果想要保存当前误差统计的运行断面,可以保存误差断面,将当前的断面以及误差统计结果一起保存起来;误差统计历史记录包含所有保存的误差断面,内容包括统计时间、平均误差、执行用户、值班主机、断面名称以及操作信息等内容;。
