基于Matpower的电力系统潮流计算原理及仿真设计（详细）

---

第一部分 前言

  由于之前有关电力系统潮流计算的文章分成了两部分，不方便大家理解潮流计算，因此我将潮流计算的原理和仿真合在一起完成本篇文章。
  电力系统分析中, 潮流计算的意义十分重大,它是在电网正常或故障情况下的稳定运行状态的计算，为了更加深入地理解潮流计算，本文基于牛顿-拉夫逊潮流计算原理，设计出一个8机28节点的电力系统网络结构，然后基于Matlab进行仿真运行分析。
  各类配电网潮流算法性能通常从以下几个方面进行分析: 算法的收敛速度、稳定性、算法的复杂程度。潮流计算的稳定性对于维护电力系统稳定潮流计算。它根据给定的初始值和网络结构确定整个系统的运行状态, 确定各个母线上的电压幅值及相角、网络中的功率分布以及网络损耗等。潮流计算的基本方程是由电力系统的网络方程得到的。在电力系统中都存在一些静态装置, 例如变压器、输电线路、并联电容器和电抗器等, 它们可以由电阻、电感和电容等基础元件构成的等值电路来模拟。

---

第二部分 牛顿-拉夫逊算法潮流计算的基本原理

  牛顿-拉夫逊法实质上就是切线法，是逐步线性化的方法。牛顿-拉夫逊计算法不仅用于求解单变量方程，还可以求解多变量非线性代数方程。

1.牛顿-拉夫逊计算法潮流计算原理

  设有单变量非线性方程

  在求解此方程的时候，需要先给出解的近似值$x^{(0)}$,它的与真解的误差将满足下列方程式

  将上述式子左边的函数在$x^{(0)}$附近展开成泰勒级数得到：

  如果差值 的值很小，差值的二次及以上阶次的各项均可以省去，于是可以简化成：

  这是对于变量的修正量线性方程式（修正方程式），据此可以得到修正量

  但是修正后的近似解$x^{(1)}$与真解仍然存在误差，为了进一步逼近真解，可以一直迭代下去，最终可以得到迭代计算通式：

  而迭代过程中的收敛判据为：

  由此可见，牛顿-拉夫逊法实质上就是切线法，是逐步线性化的方法。牛顿-拉夫逊计算法不仅用于求解单变量方程，还可以求解多变量非线性代数方程。
设有n个联立的非线性代数方程
  假定已经给出各变量的初值$x{1}^{(0)}$,$x{2}^{(0)}$,…,$x{n}^{(0)}$,令
  分别为各个变量的修正量，使其满足方程组，即
  将n个多元函数在初始值附近分别展开成泰勒级数，并略去二次及以上阶次的各项，可以得到牛顿-拉夫逊法的修正方程式
  利用高斯消去法或者三角分解法可以解出修正量,然后对初始值近似解进行修正
  如此反复迭代，在进行第k+1次迭代时，从而求出修正方程式：
  得到修正量
  并对各变量进行修正，最终可以表示为
  其中，**$X$**是由n个变量组成的n维矩阵列向量，$F（X）$,是由n个多元函数组成的n维向量:$J$是n*n阶方阵，称为雅可比矩阵。迭代过程一直进行到满足收敛判据

2.牛顿-拉夫逊计算法潮流直角坐标计算原理

  由于节点电压可以采用不同的坐标系表示，牛顿-拉夫逊潮流计算也将相应地采用不同的计算公式。
  采用直角坐标时，节点电压可以表示为
  导纳矩阵元素可以表示为
  假定系统中的第1,2，…，m号节点为PQ节点，第i个节点的给定功率设为 $Pis$和$Qis$，对于该节点可以写出如下方程
  假定系统中的第m+1，m+2，…，n-1号节点为PU节点，对于其中每一个节点可以写出如下方程
  第n号节点为平衡节点，其中电压
  是给定的，所以不参加迭代，对于PQ节点方程式和PU节点方程式来说，总共包含2(n-1)个方程式，带求的变量也是2（n-1）个，因此，我们可以得到如下的修正方程式
  其中修正方程式还可以写成分块矩阵的形式
  其中，
  对于PQ节点

  对于PU节点
  雅可比矩阵各元素都是节点电压的函数，它们的数值将在迭代的过程中不断地改变，因此在用牛顿-拉夫逊法计算潮流首先要输入网络的原始数据以及各节点的给定值并且形成节点导纳矩阵，输入节点电压初值，设置迭代计数k=0,然后开始进入牛顿-拉夫逊算法的迭代过程。

3.牛顿-拉夫逊计算法潮流计算流程

第三部分 Matpower仿真软件

1.Matpower软件简介

   MATPOWER是基于 Matlab 语言的m文件组建包，作为电力系统潮流计算仿真软件。潮流计算仿真中，按照MATPOWER中潮流计算程序的格式修改或重新编写一个新的程序“caseX.m”，在MATLAB软件中的命令窗口输入runpf（‘程序名’）就可以通过MATPOWER已经编好的程序进行潮流计算，而通过调用函数runpf的参数计算潮流。

2.Matpower矩阵参数

（1） baseMVA
  baseMVA是一个标量，用来设置基准容量。对于计算中采用有名值，可以根据实际情况设置，在本文设计中设置为100MVA。

---

（2） bus data矩阵
  bus矩阵是关于电网母线节点参数设计的矩阵，通过设置母线的参数来达到计算的目的，如表3.1所示，bus data有如下参数

Bus_i	type	Pd	Qd	Gs	Bs	Area	Vm	Va	baseKV	zone	Vmax	Vmin
1	1	97.6	44.2	0	0	1	1.0393896	-13.5366	311	1	1.06	0.94

  其中，各项参数的含义为：

• bus_i用来设置母线编号，范围为1~29997。
• type用来设置母线类型，1为PQ节点，2为PV节点，3为平衡节点；
• Pd用来设置母线注入负荷的有功功率
• Qd用来设置母线注入负荷的无功功率
• Gs用来设置与母线并联的电导
• Bs用来设置与母线并联的电纳
• area 用来设置电网断面号，可设置范围为1~100，一般设置为1
• Vm用来设置母线电压的幅值初值
• Va用来设置母线电压的相角初值
• baseKV用来设置该母线的基准电压
• zone用来设置省耗分区号，可设置范围为1~999，一般设置为1
• Vmax用来设置工作时母线电压最高幅值
• Vmin用来设置工作时母线电压最低幅值

（3） generator data矩阵
  generator data矩阵是用来设置电网中发电机的参数，各项参数名称如表3.2所示：

Bus	Pg	Qg	Qmax	Qmin	Vg	mBase	status	Pmax	Pmin	Pc1	Pc2	Qc1min	Qc1max	Qc2min	Qc2max	amp_agc	Ramp_10	Ramp_30	Ramp_q
17	650	207	300	150	1.05256	100	1	725	0	0