华中科技大学现代电力系统分析潮流计算作业

  • 格式:docx
  • 大小:143.03 KB
  • 文档页数:13

现代电力系统分析作业 ------基于Matpower的电力系统潮流计算 专 业: 班 级: 姓 名: 学 号: 目录 基于Matpower的电力系统潮流计算 1.本次潮流计算的目的及意义 本次潮流计算的目的及意义主要是了解电力系统潮流计算及其意义,在此基础上,了解电力系统潮流计算的模型以及常规的潮流计算的方法,掌握并熟练使用电力系统潮流计算软件Matpower。最后,利用Matpower计算美国西部电网WSCC三机九节点系统的静态潮流及最优潮流并给出分析。 2.电力系统潮流计算及其意义 电力系统潮流计算是研究电力系统稳态运行情况的基本电气计算,电力系统潮流计算的任务是根据给定的网络结构及运行条件,求出电网的运行状态,其中包括各母线的电压、各支路的功率分布以及功率损耗等。 潮流计算分为离线计算和在线计算两大类。离线计算可以用于电力系统调度,并确定系统的运行方式;离线计算的结果还可以用于电力系统规划方案的分析以及优化系统的运行状态;此外离线潮流计算可以作为初值,用于配合系统的故障分析以及稳定性分析。在线计算可以用于能量管理系统EMS中,主要用作静态分析;在线计算结果还可以帮助调度员合理安排系统潮流以及进行系统的安全性分析;此外在线潮流计算还可以帮助优化系统的潮流及运行状态。 3.电力系统潮流计算常规方法 牛顿-拉夫逊法

电力系统潮流计算本质上是多元非线性方程组的求解问题,而求解的方法有很多种。最常见的方法则是牛顿—拉夫逊法、高斯-塞德尔法、快速解耦法及直流潮流法。 牛顿—拉夫逊法能有效地解决非线性方程组,有较好的收敛性;它利用了导纳矩阵的对称性、稀疏性以及节点标号顺序优化的技巧,已经成为电力系统潮流计算的最常用方法。因此,本次潮流计算将重点了解牛顿—拉夫逊法。 牛顿—拉夫逊法实质上就是切线法,是一种逐步线性化的方法。设有非线性方程组 





12112212,,,0,,,0,,,0n

n

nn

fxxx

fxxx

fxxx







(1.1)

其近似解为(0)1x,(0)2x,...,(0)nx。设近似解与精确解分别相差(0)1x,(0)2x,…,(0)nx,则可得到 (0)(0)(0)(0)(0)(0)11122(0)(0)(0)(0)(0)(0)21122(0)(0)(0)(0)(0)(0)1122(,,...,)0(,,...,)0 (,,...,)0nn

nn

nnn

fxxxxxxfxxxxxxfxxxxxx



(1.2)

将上式(1.2)中的多元函数在初值附近分别展开成泰勒级数。则可得到 (0)(0)(0)(0)(0)(0)111

1121212000(0)(0)(0)0(0)(0)2222121212000(0)(0)(0)012110(,,...,)...0(,,...,)...0 (,,...,)nn

n

nnn

nnnnffffxxxxxxxxxffffxxxxxxxxxfffxxxxxx()()(0)(0)2200...0nn

n

fxxx



(1.3)

忽略上式(1.3)中(0)1x,(0)2x,…,(0)nx的二次及以上阶次的各项,整理便可以得到 11112000(0)(0)(0)(0)1121222(0)(0)(0)(0)212212000(0)(0)(0)(1212000(,,...,)(,,...,)(,,...,)nn

nn

nnnnnn

n

fffxxxfxxxxffffxxxxxxxfxxxxfffxxx

0)



(1.4)

上式(1.4)是对于修正量(0)1x,(0)2x,…,(0)nx的线性方程组,称为牛顿法的修正方程式。利用高斯消去法或三角分解法可以求出修正量(0)1x,(0)2x,…,(0)nx。然后对初始解近似解进行修正 (1)()()(1,2,...,)kkkiiixxxin, (1.5)

经过k次迭代即可将上式(1.4)写成如下形式 ()()()(1)()()()kkkkkkFXJX

XXX (1.6)

如此反复的迭代,直至满足收敛判据 ()()()121max(,,...,)kkk

infxxx (1.7)

或 ()2maxk

ix (1.8) 将牛顿-拉夫逊法用于潮流计算,要求潮流方程写成形如方程式(1.1)的形式,由于节点电压可以采用不同的坐标表示,牛顿—拉夫逊潮流计算也将相应地采用不同的计算公式。 节点电压用直角坐标表示时的牛顿—拉夫逊潮流计算

采用直角坐标系时,节点电压可表示为

iiiVejf (1.9) 导纳矩阵元素表示为 ijijijYGjB (1.10)

则可得到节点相应的有功功率和无功功率为 

1111nniiijjijjiijjijjjjnniiijjijjiijjijjjjPeGeBffGfBeQfGeBfeGfBe











 (1.11)

假定系统中的第1,2,…,m号节点为PQ节点;假定系统中的第m+1,m+2,…,n-1号节点为PV节点,第i个节点的给定功率设为isP和isQ,对节点可列写如下方程



1111222201,2,...,101,2,...,=0 nniisiisiijjijjiijjijjjjnniisiisiijjijjiijjijjjjiisiiPPPPeGeBffGfBeinQQQQfGeBfeGfBeimVVef



,

, 1,...,1imn, (1.12)

由于第n号节点为平衡节点,其电压nnnVejf是给定的,故不参加迭代。 式(1.12)总共包含了2(n-1)个变量,而由(1.9)可知,待求的变量也是2(n-1)个。还注意到(1.12)已经具备了方程组(1.1)的形式。因此即可得到如下的修正方程式 WJV (1.13)

式(1.13)中 而J则为雅克比矩阵。 有了这些,即可按照上面的牛顿—拉夫逊法计算电力系统的潮流。由于极坐标下的牛顿—拉夫逊潮流计算的推导过程与上类似,这里就不在详尽介绍。 牛顿—拉夫逊法潮流计算程序框图 牛顿—拉夫逊法潮流计算程序框图如下图所示。 图 3-1牛顿—拉夫逊法潮流计算程序框图 4.选用的潮流计算的系统 本次选用的潮流计算系统为美国西部电网WSCC三机九节点系统。该系统比较经典,经常被业内人士用来作为分析的典型例子。其接线图如下图所示 图 4-1三机九节点系统连接图 统的基准值和频率取为: 系统中设备和线路的参数如表1~4所示: 表4-1 发电机参数标幺值 参数名称 VN/kV xd x'd xq x'q xl T'd0 T'q0 H D G1 0s 0 G2 18 6s 0 G3 0 表2 变压器参数的标幺值 变压器名称 T1 T2 T3 电抗X(pu) 变比K 1 1 1 表3 线路参数的标幺值 线路名称 line1 line2 line3 line4 line5 line6 阻抗 Z(pu) R L 导纳B(pu) 表4 负荷参数 负荷名称 LumpA LumpB LumpC 有功/MW 125 90 100 无功/Mvar 50 30 35 5.利用软件matpower计算潮流 matpower简介 Matpower更准确地讲,应该是matlab中的一个工具箱,而不是一个独立的软件。它是由卡奈尔大学电气学院电力系统工程研究中心的RAY ,CARLOS 和甘德强在ROBERT THOMAS的指导下开发出来的。最初的基于MATLAB的电力潮流和最优潮流代码是为POWERWEB项目的需要而编写的。 Matpower是一个基于matlab m文件的组建包,用来解决电力潮流和优化潮流的问题。它致力于为研究人员和教育从业者提供一种易于使用和可更新的仿真工具。Matpower的设计理念是用尽可能简单、易懂,可更新的代码来实现最优秀的功能。 基态潮流计算 基态潮流计算条件说明及数据输入 本次潮流计算时选取Bus1为平衡节点,Bus2,Bus3为PV节点,其余的节点为PQ节点;取发电机G2、G3的输出有功功率为其额定功率,则依据matpower用户手册即可输入以下数据(数据说明见注释): function [baseMVA,bus,gen,branch,areas,gencost] = case9hbyq %定义case9hbyq为一个含变压器的三机九节点电力系统的潮流算例 version ='2'; %定义系统的基值为100MVA baseMVA = 100; %bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vmin %%母线参数,定义母线的序号、节点类型(PV、PQ、Swing节点)、注入的有功/无功、分区、 %%母线电压标幺值、基值、相角、电压标幺值的上下限等 bus = [ 1 3 0 0 0 0 2 1 0 1 ; 2 2 0 0 0 0 1 1 0 18 1 ;