郭伟+基于pq分解模型的静态安全分析 - 图文

更新时间:2024-06-29 21:39:01 阅读量: 综合文库 文档下载

说明:文章内容仅供预览,部分内容可能不全。下载后的文档,内容与下面显示的完全一致。下载之前请确认下面内容是否您想要的,是否完整无缺。

攀枝花学院本科毕业设计(论文)

基于P-Q分解模型的静态安全分析

学生姓名: 郭 伟 学生学号: 201010502032 院(系): 电气信息工程学院 年级专业: 2010级电气工程与自动化 指导教师: 户 秀 琼 讲 师 助理指导教师:

二〇一四年五月

攀枝花学院本科毕业设计(论文) 摘要

摘 要

在电力系统中,线路传输功率是否过载、节点电压是否越限直接影响电力系统的静态安全,因此,有必要对电力系统进行静态安全分析。

论文主要针对基于P-Q分解模型的静态安全分析方法及其在MATLAB中的实现进行了深入的研究。考虑到基于P-Q分解模型的静态安全分析方法来源于潮流计算中的P-Q分解模型,论文首先研究了基于P-Q分解模型的潮流计算。而后,针对支路断线的预想故障,对基于P-Q分解模型的静态安全分析方法进行了研究。最后,根据给定的5节点环网系统数据,采用MATLAB实现了基于P-Q分解模型的潮流仿真计算;在此基础上,实现了基于P-Q分解模型的静态安全分析仿真研究,并对断线后系统的运行状态进行了分析,给出了系统断线后的安全状态。

关键词 电力系统,静态安全分析,P-Q分解模型,MATLAB,5节点环网

攀枝花学院本科毕业设计(论文) ABSTRACT

ABSTRACT

In the power system, whether the transmission power is overloaded , or node voltage is over the limits will absolutely take a effect on the static security of the power system. Therefore, it is necessary to take a static security analysis on the system. This paper mainly conducts in deeply research on static security analysis method based on P-Q decomposition model and its implementation in the MATLAB environment. first of the paper, power flow calculation using P-Q decomposition is the primary task, considering the method, static security analysis method, simplified from P-Q decomposition model. Then, the static security analysis method based on P-Q model will be studied, assuming the one branch breakdown. In the end, according to the given 5_node loop networks data, achieved this mission, flow calculating was successful using P-Q decomposition in MATLAB environment. Further more, the simulation of static security analysis based on P-Q decomposition model was achieved. While comparing the system operation state after one branch breakdown with the stable state, no one branch breakdown. Naturally the final result includes secure state of one branch-outage system.

Key words Power system, Static security analysis, P-Q decomposition method,

MATLAB, 5-bus loop network

攀枝花学院本科毕业设计(论文) 目录

目 录

摘 要 ................................................................................................................................. I ABSTRACT ..................................................................................................................... II

1 绪论 .............................................................................................................................. 1

1.1课题的研究背景和研究意义 ................................................................................... 1 1.2 课题的国内外研究现状 ......................................................................................... 1

1.2.1 电力系统静态安全分析的概念及分类 ··································································· 1 1.2.2 静态安全分析的发展历史 ······················································································ 2 1.3 论文的主要研究内容 .............................................................................................. 3

2基于P-Q分解模型的静态安全分析原理 .................................................................. 4

2.1 基于P-Q分解模型的潮流计算原理 ........................................................................ 4 2.2 基于P-Q分解模型的静态安全分析 ....................................................................... 11

3 仿真分析 ..................................................................................................................... 15

3.1 MATLAB语言简介 ............................................................................................... 15

3.1.1 MATLAB简介 ······································································································· 15 3.1.2 MATLAB语言简介 ······························································································· 15 3.2基于P-Q分解模型的静态安全分析在MATLAB中的实现 ..................................... 16

3.2.1 总体流程图 ··········································································································· 16 3.2.2 原始数据表及分析 ································································································ 16 3.2.3 基于P-Q分解模型的潮流计算仿真分析 ····························································· 20 3.2.4 基于P-Q分解模型的静态安全分析的仿真研究 ·················································· 21

结论 ................................................................................................................................. 25 参考文献 ......................................................................................................................... 28 附录A:基础数据生成子程序 ..................................................................................... 29 附录B:潮流计算子程序 .............................................................................................. 31 附录C:支路传输功率计算子程序 ............................................................................. 33 附录D:支路开断时首末节点注入功率增量计算子函数 ......................................... 35 附录E:修改节点注入功率子程序 .............................................................................. 35 附录F:静态安全分析主程序 ...................................................................................... 37 致 谢 ............................................................................................................................... 42

攀枝花学院本科毕业设计(论文) 1绪论

1 绪论

1.1课题的研究背景和研究意义

电力系统作为当今规模最大,结构和技术含量最高的人工实时自动化系统,其发展的自动化程度可以作为一个国家综合实力的重要标志。目前,我国电网的规模已居世界第一位,无论是在电网规模电压等级还是自动化水平上都可以说世界领先,但同时其安全性和稳定性也正经受着严峻的考验[1-3]。尤其是在见证了国内外电力系统事故造成的损失后,电力系统的稳定性和安全性问题日趋凸显,比如:1978年法国电网因负荷的迅速攀升,导致电网电压迫降而失稳,造成大面积停电事故;2003年由于输电线路过负荷引发的美加大停电,造成了巨大的经济损失,对社会造成了严重的恶劣的影响。因此,为了保证系统的安全稳定运行,并且保证系统供电的持续性,迫切需要对电力系统进行实时运行状态评估,以提高其应对事故的能力。电力系统静态安全分析作为电力系统实时运行状态评估的方法之一,在电力系统的研究中占据了重要的位置[2]。

利用静态安全分析可以对系统进行事故预想,对于一个输电系统规划方案而言,可以校验其承受事故的能力;对运行中的系统而言,可以校验其运行方式及接线方式的安全性,进而给出事故前后应采用的防范措施或校正措施。迄今为止,国内外开发研究了许多专门用于静态安全分析的方法。其中基于P-Q分解模型的静态安全分析方法由于具有较好的精度,可以提供电力系统所有状态变量的变化情况,在电力系统的静态安全分析中得到了广泛使用[2]。

随着计算机技术的发展以及电网规模的增大,数字仿真技术对于电力系统的分析计算变得至关重要。MATLAB由于其强大的计算功能,为电力系统的仿真分析提供了技术支持[4,5]。

因此,研究基于P-Q分解模型的静态安全分析方法及其在MATLAB中的仿真实现,对于电网的安全稳定运行意义重大。

1.2 课题的国内外研究现状

1.2.1 电力系统静态安全分析的概念及分类

安全性指电力系统在运行中承受故障扰动(例如突然失去电力系统的元件或短路故障等)的能力[6-9]。而安全分析分为静态安全分析和动态安全分析。静态安全分析假设电力系统从事故前的静态直接转移到事故后得另一个静态,不考虑中间暂态过程,其主要是指应用N-1原则,逐个无故障断开线路,变压器等元件,

1

攀枝花学院本科毕业设计(论文) 1绪论

检查其他元件是否因此过负荷以及电网电压是否越限,用以检验电网结构强度和运行方式是否满足安全运行要求。而且,电力静态安全分析仅考虑事故后稳态运行情况的安全性,研究电力系统中元件断开引起支路有功潮流及母线电压越限,如果出现越限就要采取相应的校正控制策略以保证系统的正常运行[2]。

从不同角度可将静态安全分析分成若干类,下面逐一说明。

①依据生成导纳矩阵,即得到事故后的网络拓扑分类中有三种模式[10,11]。 1)第一种:首先假设具体某条线路断开,然后针对事故后的网络生成导纳矩阵,再进行潮流计算,得到状态参数,最后和规定的章程来对比判断是否越限;

2)第二种:首先计算潮流,然后假设某条线路断开,此时并不重新计算故障后的潮流,而是从用原来的潮流的参数结合特定潮流算法的特点导出故障后的状态参数,再查规定的章程来判断是否越限。很显然,第一种最贴近实际,但计算量会随着剧增,第二种就不存在这种问题。

3)第三种:首先假设每条线路开断,由于单条线路开断只影响导纳矩阵中的四个参数,依据规律得出断线后的导纳矩阵,在进行潮流计算,得到状态参数,最后和规定的章程来对比是否越限。

②依据开断原件类型,静态安全分析可分为两种模式[2]。 1)第一种:针对发电机开断的静态安全分析。 2)第二种:针对输电线路开断的静态安全分析。

③依据导出静态安全分析的潮流算法不同,静态安全分析可分为三种模式[2]。 1)第一种:基于直流潮流法的静态安全分析。

2)第二种:与Newton-Raphson潮流算法相结合的直接法。

3)第三种:与P-Q分解潮流算法(快速解耦潮流算法)相结合的直接法,即是基于P-Q分解模型的静态安全分析方法。

④依据开断元件的个数,静态安全分析可分为两种模式[12]。

1)第一种:针对特定元件(重负荷线路或者病态电路)进行的静态安全分析。 2)第二种:对全网元件逐个进行开断分析的静态安全分析。

无论静态安全分析分类如何,其思想始终是研究某元件开断后进入新稳态的电网是否能在章程规定的参数范围内稳定运行[13]。具体来说,即是研究元件开断后,系统各个节点以及各条支路的状态参数是否越限。

而本论文主要针对系统的断线故障,使用与P-Q分解潮流算法相结合的直接法对线路开断后的电网进行静态安全状态的判断。

1.2.2静态安全分析的发展历史

自上世纪60年代以来,大面停电造成的巨大损失,使各国对电力系统的安全分析,给予了足够的重视,成为上世纪七、八十年代非常活跃的研究领域。近年

2

攀枝花学院本科毕业设计(论文) 1绪论

来,国外的一些控制中心也配备了具有实时功能的静态安全分析软件,表明文的研究已进入了实用化阶段。目前电力系统中的很多安全稳定分析功能都处于离线应用状态。离线运用在计算时间上要求不高,但是不能实时跟踪电力系统的变化。在线运用要求计算时间短,这对于那些规模较大的电力系统以及安全稳定分析方法计算量较大的情况而言,确实是一个较为困难的问题。因此,现在的电力系统针对安全稳定分析这方面都在强调在线运用,以提高其实用性[14]。

1.3 论文的主要研究内容

本文主要是研究基于P-Q潮流模型的静态安全分析方法及其在MALTAB编程环境下的仿真实现。其主要研究内容如下:

①研究基于P-Q分解模型的潮流计算原理

基于P-Q分解模型的静态安全分析其基本的思想来源于基于P-Q分解模型的潮流计算,因此,有必要先对基于P-Q分解模型的潮流计算原理进行研究,得到节点电压幅值与系统无功之间的关系、系统有功与节点电压相角之间的关系,从而为后续静态安全分析中判断系统的安全状态奠定基础。

②研究基于P-Q分解模型的静态安全分析

在研究基于P-Q分解模型的潮流计算原理的基础上,针对系统的断线故障,研究基于P-Q分解模型的静态安全分析方法。考虑到实际系统中,N-1故障发生的概率远远大于其他多重故障,因此,论文主要针对单条线路断开的故障进行研究。通过线路的断线模拟分析,得到单条线路断开后,系统节点电压幅值的表达式,以便于对系统断线后的安全状态进行判断。

③基于P-Q分解模型的静态安全分析在MATLAB中的仿真实现

针对给定的5节点环网系统,得到其等值电路之后,在学习MATLAB语言的基础上,首先实现该系统基于P-Q分解模型的潮流仿真计算,得到系统在正常运行状态下,各个节点的电压幅值以及各条支路上传输的有功功率。而后,针对单条线路断开的预想故障,实现基于P-Q分解模型的静态安全分析仿真研究,判断系统断线后的安全状态。

3

攀枝花学院本科毕业设计(论文) 2基于P-Q分解模型的静态安全分析原理

2 基于P-Q分解模型的静态安全分析原理

2.1 基于P-Q分解模型的潮流计算原理

要研究基于P-Q分解模型的潮流计算原理,首先从所有潮流算法共同体现的支路潮流计算开始[6-8]。图2.1是支路潮流的示意图。图中,i,j表示支路的首末节点,gij+jbij表示ij支路之间的导纳,bi0表示ij支路的对地电纳,Vi,Vj分别表示节点i,j的节点电压相量,Pij?jQij表示ij支路传输的功率,Iij表示ij支路流过的电流,

'表示除去对地支路上流过的电流后ij支路上的Ii0表示对地支路上流过的电流,Iij电流。

Pij?jQijgij?jbijVi?i??'IijIijjbi0?Ii0Vjjjbi0

图 2.1 支路潮流示意图

由图2.1可以得到如式(2.1)所示的支路传输功率的表达式。

*'**Pij+jQij=ViIij=Vi(Iij+Ii0)=Vi[(Vi-Vj)(gij+jbij)+jVbii0]

式(2.1)

由于有:

Vi?Viei 式(2.2) Vj?Vjej 式(2.3)

j?j?ej?ij?cos?ij?jsin?ij 式(2.4)

式(2.2)~式(2.4)中,?i,?j,?ij分别表示节点i,j的电压相角及其相角差。

将式(2.2)~式(2.4)代入式(2.1),则有:

2P?Vijigij?VVij(gijcos?ij?bijsin?ij) 式(2.5) 2Q?VV(bcos?gsin)??V(b) 式(2.6) ijijijijijijiijbi0??式(2.5)和式(2.6)是所有潮流计算方法中所共有的数学模型。 下面就针对潮流计算方法中的Newton-Raphson法进行阐述,该方法是基于P-Q分解模型的潮流计算方法的基础。

①Newton-Raphson潮流算法 1)节点功率平衡方程

由式(2.5)和式(2.6),可以得到节点注入功率,其为与该节点相连的所有支路的功率之和,其表达式如式(2.7)和式(2.8)所示。

4

攀枝花学院本科毕业设计(论文) 2基于P-Q分解模型的静态安全分析原理

Pi??Pij??[Vi2gij?VVij(gijcos?ij?bijsin?ij)] 式(2.7)

j?1j?1nn Qi??Q??[VV(bijijj?1j?1nnijcos?ij?gijsin?ij)?Vi2(bij?bi0)] 式(2.8)

假定节点i的给定注入功率分别为Pis与Qis,结合式(2.7)和式(2.8),可得到式(2.9)和式(2.10)所示的节点功率平衡方程。

nn?Pi?Pis??Pij?Pis?Vi?Vj(gijcos?ij?bijsin?ij)j?1j?1 nn?Qi?Qis??Qij?Qis?Vi?Vj(bijcos?ij?gijsin?ij) j?1j?12)Newton-Raphson法的数学基础 设有非线性方程如式(2.11)所示。

f1(x1)?0 其近似解为x(0)1,其精确解与近似解之间相差?x1,即是:

x(0)1=x1??x1 因此有:

f(0)1(x1??x1)?0 按Taylor级数展开并略去高次项,则有:

f1(x(0)1)?dfdx?x1?0 ?x(0)df11??f1(x1)/dx(0) 1令:

?x(0)1??f(0)1(x1)/df1dx(0) 1x(0)(01?x1??x)1

再令:

x(1)1?x(0)1??x(0)1 ?x(2)?f((1)df11?1x1)/dx(1) 1如此迭代,其通式如下:

5

式(2.9)

式(2.10)

式(2.11)

式(2.12)

式(2.13)

式(2.14) 式(2.15)

式(2.16) 式(2.17) 式(2.18)

式(2.19)

攀枝花学院本科毕业设计(论文) 2基于P-Q分解模型的静态安全分析原理

(t?1)(t)?x??f(x111)/df1 式(2.20) (t)dx1(t?1)()t(t) 式(2.21) x?x?x11?1(t) |xx??1?1|1 式(2.22)

式(2.22)中,?1为给定的允许误差。

将上面的过程扩展到如式(2.23)所示的多变量方程组:

f1(x1,x2,,xn)?0?f?2(x1,x2,,xn)?0?? ?fn(x1,x2,,xn)?0??按Taylor公式展开,并略去高次项,得:

f(0)(0)?f11(x1,x2,,x(0)n)??f1?x(0)?x(0)1?1?x(0)?x(0)?n?0n?

??? f(x(0)(0)?fn?fn(0)n1,x2,,x(0)n)??x(0)?x(0)1?1?x(0)?x?0?n?n??将其写成矩阵的形式有:

??f1?f1?f1??f(0)(0)(0)??x?x1(x1,x2,,xn)??1,2,?xn??,?(0)?????x1??????????? ?fn(x(0)1,x(0)(0)2,,xn)?????fn?fn?f(0)n???????x1,?x2,,?x??xn?n??进一步可以表示成式(2.26):

F?J?X 其中:

?(0)(0)(0)F??f1(x1,x2,,xn)???? ??fn(x(0)(0)(0)1,x2,,xn)????x(0)??X(0)=?1??? ???x(0)n?? 6

式(2.23) 式(2.24)

式(2.25) 式(2.26)

式(2.27a)

式(2.27b)

攀枝花学院本科毕业设计(论文) 2基于P-Q分解模型的静态安全分析原理

??f1?f1??x?x?1,2,?J????fn?fn??x1?x2?,,?进而有:

?f1??xn?,??? 式(2.27c) ?fn???xn?,??X?J?1F 式(2.28)

类似于一维方程,n维方程的迭代通式可用式(2.29)来表示。

X(t?1)=X(t)-?X(t) 式(2.29a)

?X?J?1F 式(2.29b)

?X?? 式(2.29c)

式(2.29c)中,?为给定的允许误差。 3)Newton-Raphson法求解潮流的原理 由前面已经知道,节点的功率平衡方程:

nn?Pi?Pis??Pij?Pis?Vi?Vj(gijcos?ij?bijsin?ij) 式(2.30)

j?1j?1?Qi?Qis??Qij?Qis?Vi?Vj(bijcos?ij?gijsin?ij) 式(2.31)

j?1j?1nn设方程中的变量为节点电压幅值Vi、Vj、?i、?j 。对系统每个节点的功率平衡方程采用多元函数Taylor公式展开,并略去二阶及其高次项,则有:

??P?H111??????????Q??J111??????N11L11???1??????????V? 式(2.32) ??1?V1????????P??H N ???????Q?=?J L???V/V? 式(2.33) ???????H N ?令Jacobin 矩阵 A=?,则有: ??J L???P???????Q??A??V/V? 式(2.34) ????

7

攀枝花学院本科毕业设计(论文) 2基于P-Q分解模型的静态安全分析原理

?????1??P??A??V/V???Q? 式(2.35) ????

而Jacobin 矩阵的非对角元素如式(2.36)~(2.39)所示。

Hij???Pi?-VVij(gijsin?ij-bijcos?ij)??j(i?j)

式(2.36)

Nij??PiVj??VVij(gijcos?ij?bijsin?ij)?Vj??Pi?VVij(gijcos?ij?bijsin?ij)??j(i?j) 式(2.37)

Jij?(i?j) 式(2.38)

Lij???QiVj??VVij(gijsin?ij?bijcos?ij)?Vj(i?j) 式(2.39)

Jacobin 矩阵的对角元素如式(2.40)~(2.43)所示。

n??PiHii???Vi?Vj(gijsin?ij?bijcos?ij)(i?j) 式(2.40)

??ij?1j?iNii???PiVi??Vi?Vj(gijcos?ij?bijsin?ij)?2Vi2gij?Vij?ij?i(i?j) 式(2.41)

Jii???Pi??VVij(gijcos?ij?bijsin?ij)??i(i?j) 式(2.42)

n?QiLii?Vi??Vi?Vj(gijsin?ij?bijcos?ij)?2Vi2bii?Vij?1j?i(i?j) 式(2.43)

由于Newton-Raphson 法采用Taylor 公式对方程进行近似线形化,而Taylor公式要求?X是很小的数,即要求在Newton-Raphson中,设定的初值要接近于真值,否则会带来收敛问题[8]。 ②基于P-Q分解模型的潮流算法

P-Q分解法的思想是结合电网实际情况,在Newton-Raphson潮流极坐标算法数学的基础上,针对高压电网的特点进行了如下简化。

in?0,cos?11)由于两节点间的相角差?ij?10~20,认为:s; ij?ij?2)线路参数bij??gij;

3)节点电压一般在额定值附近,即对于i, j两节点有Vi?Vj ;

4)由于输电元件(线路或变压器)中的电抗极大于电阻,所以一个节点传送到另一个节点的有功功率,主要取决于两节点间的电压相角差,而不是两个节

8

攀枝花学院本科毕业设计(论文) 2基于P-Q分解模型的静态安全分析原理

点的电压幅值。而一个节点传送到另一个节点的无功功率,主要取决于两节点间的电压幅值差,而不受到电压相位差的影响[8]。

基于以上几点假设,将式(2.33)进行简化,可以得到式(2.44)。

??P??H??=?0??????? 式(2.44) ???Q??0L???V/V?即有:

??P???H?????

??Q???L?[?V/V] 同时有:

H2ii?ViBii Hij?ViVjBij LV2ij?iBij Lij?ViVBjij 当网络拓扑中没有PV节点时,H=L。

?V1VBn1n?H?L??V21B11???? ?VVBn1n1V2nBnn??将式(2.45)和式(2.46)用矩阵乘积形式来表示:

???P1??V0??B1n????1????B11??V10????1????P??n????0V??n????Bn1B?????nn??????0Vn???????n????V10??B 1n??????B11??V??1???0V??n????Bn1B???nn????V???n?????Qn?0????V1????B11B1n?????????Qn????0V??n????Bn1B???V1??? V?nn?????n??进而可以得到节点电压幅值的修正量和节点电压相角的修正量:

9

式(2.45)

式(2.46)

式(2.47)

式(2.48)

式(2.49)

式(2.50)

式(2.51) 式(2.52)

式(2.53) 攀枝花学院本科毕业设计(论文) 2基于P-Q分解模型的静态安全分析原理

??V1??B11??????????Vn????Bn1B1n???Bnn???1?V1????00???Vn???1??Q1??B11??????????Qn????Bn1B1n???Bnn???1?1??Q1/V1??? ?????Qn/Vn??式(2.54)

???1??V1???????????n????0?V1?????00???Vn??0???Vn???1?B11????Bn1?B11????Bn1B1n???Bnn??B1n???Bnn???1?V1????00???Vn????P1????????Pn???1?1??P1/V1????????Pn/Vn?? 式(2.55)

得到电压幅值和电压相角的修正量之后,运用迭代法求方程组的思想,解节点

功率平衡方程组。具体步骤如下:

1) 设定变量电压和相角的初值;

2) 运用公式(2.31),求解有功功率误差量; 3) 运用公式(2.55),求解相角误差量;

(t?1)(t)(t)4) 修正相角的误差? ; ?????5) 运用公式(2.30),求解无功功率误差量; 6) 运用公式(2.54)求解电压误差量;

(t?1)(t)(t)7) 修正电压V ; ?V??V8) 依收敛判据?V(t)??1,??(t)??2不收敛,则从2)开始进入下一次循环。

(?1,?2为允许误差),检查是否收敛,如

针对有n个节点的网络(假定该系统中有一个平衡节点,m个PV节点,n-m-1个PQ节点),分别对Newton-Raphson 法和P-Q分解法的特点进行如下简要分析。

Newton-Raphson 法,用Taylor公式进行近似线形化,精度高,呈平方收敛,由于Jacobin为n-1+n-m-1阶矩阵,且每次运算都需要更新,单次迭代耗时较多。

P-Q分解法在Newton-Raphson 法基础上,结合高压电网特点,做出简化,由于基础模型相同(式(2.7)和式( 2.8)),其在相同收敛判据下最终结果和Newton-Raphson法结果有相同的精度。而P-Q分解法把n-1+n-m-1维的每次迭代都变化的矩阵,依高压电网能够实际条件拆解成n-1和n-m-1维的两个常数矩阵,在单次迭代时间上相对于Newton-Raphson法大大提高。虽然单次迭代P-Q分解法精度低于Newton-Raphson法,但在相同网络拓扑下,P-Q分解法在整个潮流计算上消耗机时要少得多,即P-Q分解法用单次计算时间的减少以提高单次运算速度,依相同的收敛判据通过额外的运算达到和Newton-Raphson法相同精度[8]。

10

攀枝花学院本科毕业设计(论文) 2基于P-Q分解模型的静态安全分析原理

2.2 基于P-Q分解模型的静态安全分析

解决方案采用注入功率模拟线路开断,实际线路并未开断,只是用附加的节点注入功率改变了与节点相连的支路功率,使其和断线后的功率相同。被模拟开断的线路在附加注入功率和原始线路功率的叠加效果让线路传输功率为0,即支路开断了。附加注入功率的求取成为关键,因基于潮流算法的不同,其推导过程不同,但最基本的模拟支路开断的思想是相同的[2]。

模拟线路开断的思想如图2.2~2.4所示[2]。

Pi?jQiPj?jQjji Pi??jQi?Pij?jQij 图2.2 被开断线路断开前的功率分布

Pi?jQii'Pi??jQi'? Pj?jQjPji?jQjiPj??jQj?

j Pi?jQiPj'??jQ'j?

图2.3被开断线路断开后的功率分布

Pj?jQjji ?Pi?j?Qi?Pj?j?Qj'''Pi??jQi'?Pij?jQij图2.2~图2.4中,节点i, j表示被断开支路的首末节点;Pi?jQi,Pj?jQj则分别表示节点i, j实际的注入功率;Pij?jQij,Pji?jQji分别表示从节点i到节点j的传输功率以及从节点j到节点i的传输功率;Pi??jQi?,Pj??jQj?分别表示除ij支路外

'与节点i相连的支路、与节点j相连的支路功率之和;Pi??jQi'?,Pj'??jQ'j?分别表示

ij支路断开后,与与节点i相连的支路、与节点j相连的支路功率之和;?Pi?j?Qi,

'?Pj?j?Qj则分别表示模拟ij支路断开时,节点i、节点j的注入功率增量;Pij'?jQij,'Pji?jQ'ji则表示模拟断线时,支路ij由于节点注入功率的变化而得到的新的支路传

'Pji?jQ'jiPj'??jQ'j?

图2.4对被开断线路进行模拟开断

输功率。

11

攀枝花学院本科毕业设计(论文) 2基于P-Q分解模型的静态安全分析原理

线路开断前后节点注入功率的差值,即作为模拟线路开断的节点注入功率。以注入功率模拟线路开断,端口效果应相同,即增加注入功率后,与节点相连其他支路的功率应和真实断线后的情况一致。

由图2.2可以列写功率平衡方程式 :

式(2.56) P??jQ?jQ?P?jQiiPi?i?ijij 式(2.57) P??jQ?jQ?P?jQjjPj?j?jiji由图2.3列写功率平衡方程式:

'' 式(2.58) PjQP?jQj?j?j?j?'' 式(2.59) PjQP?jQi?i?i?i?由于用注入功率模拟线路开断,应满足除被开断支路外,支路两端节点注入

'功率和断开后的相同,即有Pi?jQi,Pj?jQj,Pi??jQi'?,Pj'??jQ'j?在图2.3和2.4中

的是相同的,因此可得:

''?PjQ??jQi?iPij?ij 式(2.60) '' 式(2.61) ?PjQ??jQj?jPji?ji而:

''P??jQ??jQP?j?Q 式(2.61) ijijPijij?ijij''P??jQ??jQP?j?Q 式(2.62) jijiPjiji?jiji其中,?Pij,?Qij和?Pji,?Qji分别表示由于节点注入功率的变化所得到的支路ij所传输功率的增量。

进一步有:

?Pij???Pi???Pij??? 式(2.63) ??P????Pj????P????ji??ji???QQ?Qij???ij?i?? 式(2.64) ???????Q??????Q?Q?ji??j??ji?以Pij, Pji为函数,Pi,Pj为变量,则写出?Pij,?Pji的全微分,有:

?P?P?PP???ijijij?ij??P?ij???P?Pij???????P?P?P?PP?P?P?????????ij??iji?j?i?i??????? 式(2.65) ?=??????????????P?P?P?P?P?P?P?P?PP???i?ijijiji?j??ji??j??j??j?ij??????P+?Pij???P?P?P?P????????ijij???? 对应无功有:

12

攀枝花学院本科毕业设计(论文) 2基于P-Q分解模型的静态安全分析原理

?QQ?ij?ij??Q?ij????????Q?QQ?Q?Q???????ij??i?j?i?i??????????? 式(2.66) ???????Q?Q?QQQ?Q?PPj??j??????ji?ji?ji?ji?i?j???????QQ??????i?j????其中,bij,bi0,bj0分别为支路ij的实际电纳、支路ij在节点i侧的对地电纳、支路ij在节点j侧的对地电纳。将公式(2.65)和公式(2.66)分别代入公式(2.63)和(2.64)求得为模拟线路开断而注入的节点功率。

?P?P?a?i??11a12?ij?? 式(2.67) ??????Pa21a22??ji??j???P?Q?Q?c?ij?i??11c12?式(2.68) ????????QccQ?j??2122???ji??

其中:

???bS??a10??????ijbij?iiSij?11a12式(2.69) ????????????????aa01b?bSS?????2122??ijij?ijj????j??

?1??S?b2bb??i?????c10ij?i0ij??????iiSj?11c12?? 式(2.70) ???????????????c01S?b?b?2b??????21c22??jiSjj????ijijj0??????? 式(2.71) S????B?1?1 S???B?? 式(2.72)

在得到式(2.67)和式(2.68)之后,就可以对节点i, j的节点注入功率进行修正,进而采用第2.1节中的式(2.54)即可对模拟断线后各个节点的电压幅值增量进行求解,并由此可判断断线后各个节点的电压幅值是否在运行范围内。

也就是说,有了模拟开断注入功率,修改稳态时节点注入功率,调用潮流计算函数,调用支路传输功率函数,计算状态量,比如:节点电压的幅值和相角,支路的传输功率。最后就可在主函数里比较这些状态量是否在容许范围内,输出各个线路开断后其它节点和支路是否安全的判断结果。

图2.5给出了论文所研究的基于P-Q分解法的静态安全分析流程图。 本章主要对基于P-Q分解模型的潮流计算原理以及基于P-Q分解模型的静态安全分析方法进行了详细的理论推导。有了这个理论基础,后面就可以针对具体给定的系统,进行MATLAB的仿真研究。

?1 13

攀枝花学院本科毕业设计(论文) 2基于P-Q分解模型的静态安全分析原理

稳态时潮流计算的基础数据调用稳态时支路传输功率求出模拟开断所需的节点注入功率修改稳态时节点注入功率潮流计算子函数支路传输功率子函数在主函数中与上下限比较分析结果:即静态安全分析结果 图2.5静态安全分析流程图

14

攀枝花学院本科毕业设计(论文) 3仿真分析

3 仿真分析

3.1 MATLAB语言简介

3.1.1MATLAB简介

MATLAB全称为matrix laboratory,矩阵实验室,最初功能仅为矩阵运算,由20世纪70年代,美国新墨西哥大学计算机科学系主任Cleve Moler为了减轻学生编程的负担,用FORTRAN编写了最早的MATLAB。1984年由Little、Moler、Steve Bangert合作成立了的MathWorks公司正式把MALAB推向市场。到20世纪90年代,MATLAB已成为国际控制界的标准计算软件[15,16]。

MATLAB是一种用于数值计算、可视化及编程的高级语言和交互式环境。使用 MATLAB,可以分 析数据,开发算法,创建模型和应用程序。借助其语言、工具和内置数学函数,可以探求多种方法,比电子表格或传统编程语言(如 C/C++ 或 Java?)更快地求取结果。MATLAB 应用广泛,其中包括信号处理和通信、图像和视频处理、控制系统、测试和测量、计算金融学及计算生物学等众多应用领域。在各行业和学术机构中,有一百多万工程师和科学家使用 MATLAB 这一技术计算语言[17]。

3.1.2MATLAB语言简介

MATLAB 提供了一种高级语言和开发工具,使用户可以迅速地开发并分析算法和应用程序[16]。

MATLAB 语言对向量运算和矩阵运算提供内在支持,这些运算是解决工程和科学问题的基础,能够实现快速开发和执行。使用 MATLAB 语言,编程和开发算法的速度较使用传统语言大幅提高,这是因为无须执行注入声明变量、指定数据类型以及分配内存等低级管理任务。在很多情况下,支持向量运算和矩阵运算就无需使用for 循环。因此,一行 MATLAB 代码通常等同于数行 C 代码或 C++ 代码。

MATLAB 提供了传统编程语言的多项功能,其中包括流控制、错误处理以及面向对象编程 (OOP)。用户既可以使用基本的数据类型或高级数据结构,也可以定义自定义数据类型。采用交互方式,一次仅执行一个命令,可以即时生成结果。这种方法可以快速试探多个选项,通过反复迭代,找出最佳的解决方案。通过捕获交互式步骤,生成可以重复使用的脚本和函数,并实现任务的自动化。

MATLAB 附加产品可针对信号处理和通信、图像和视频处理、控制系统以及许多其他领域提供各种内置算法。通过将这些算法与自己的算法结合使用,可以构建复杂的程序和应用程序[15]。

15

攀枝花学院本科毕业设计(论文) 3仿真分析

MATLAB 应用程序可以与其他语言编写的应用程序集成。在 MATLAB 中,可以直接调用以 C、C++、Java 和 .NET 编写的代码。使用 MATLAB 引擎库,可从 C、C++ 或 Fortran 应用程序调用MATLAB 代码。

MATLAB 采用处理器优化库,可以快速执行矩阵运算和向量运算。对于通用的标量计算,MATLAB 使用其即时 (JIT) 编译技术,提供了可与传统编程语言相媲美的执行速度。为了充分利用多核和多处理器计算机,MATLAB 提供了众多的多线程线性代数和数值函数。这些函数可在单个 MATLAB 会话中自动执行多个计算线程,从而得以在多核计算机上提高执行速度。

通过附加的并行计算产品,可以充分利用多核台式机和其他的高性能计算资源,如GPU 和群集。这些产品可提供高级构造,通过附加的并行计算产品,可以充分利用多核台式机和其他的高性能计算资源,如GPU 和群集。这些产品可提供高级构造,只需稍加改动 MATLAB 代码即可实现应用程序的并行化。

本章就将针对具体给定的系统,采用MATLAB对该系统进行基于P-Q分解模型的潮流计算仿真以及基于P-Q分解模型的静态安全分析方法仿真。

3.2基于P-Q分解模型的静态安全分析在MATLAB中的实现

3.2.1 总体流程图

在实施基于P-Q分解模型的静态安全分析方法的仿真研究时,其必须以基于P-Q分解模型的潮流计算为基础,因此,其总体框架如图3.1所示。

原始数据处理潮流计算仿真获取静态安全分析所用数据静态安全分析模块 图3.1 静态安全分析总体流程图

从图3.1可以知道,本章的仿真内容主要有三部分: ①原始数据处理

对原始数据进行处理,得到给定系统在标幺值下的等值电路,系统的节点信息,系统的支路信息,系统的发电机信息,系统的负荷信息。进而得到系统的节点导纳矩阵,这些数据都是潮流计算的基础,也是静态安全分析的基础。

②潮流计算仿真获取静态安全分析所用数据

这一模块即是针对给定系统,采用MATLAB对基于P-Q分解模型的潮流计算进行仿真分析。得到系统在正常运行情况下节点电压的幅值,以及支路传输的

16

攀枝花学院本科毕业设计(论文) 3仿真分析

功率。

③静态安全分析模块

该模块主要是在前面潮流计算的基础上,针对单条线路断开的预想故障,对给定系统进行静态安全分析。其详细流程可参见第2章的图2.5。

3.2.2原始数据表及分析

①网络接线及参数

给定系统的网络接线结构如图3.2所示。

42351 图3.2系统接线图

系统各项参数如表3.1~表3.3所示。需要说明的是,这些参数均是标幺值。

表3.1 节点数据

节点号 1 2 3 4 5 电压幅值 - - - 1.05 1.05 电压相角 - - - - 0 发电机有功 - - - 5 -

表3.2 发电机出力限制 节点号 4 5 有功上限 8 8 有功下限 1 1

发电机无功 - - - - - 负荷有功 1.6 2 3.7 0.0 0.0 负荷无功 0.8 1 1.3 0.0 0.0 无功上限 3 5 无功下限 -3 -2.1 17

攀枝花学院本科毕业设计(论文) 3仿真分析

表3.3支路参数

支路号 1 2 3 首末节点号 1-2 1-3 2-3 支路电阻 0.04 0.1 0.08

表3. 4变压器数据

首末节点号 2-4 3-5 电阻 0.0 0.0 电抗 0.015 0.03 变比 1.05 1.05 变比上限 1.1 1.1 变比下限 0.9 0.9 变比步长 0.05 0.05 支路电抗 0.25 0.35 0.3 对地导纳(b/2) 0.25 - 0.25 ②节点分类

系统的节点可分为三类 :PV节点,PQ节点,平衡节点(习惯用slack表示)。 1) PV节点:给定节点注入有功功率和电压幅值。 2) PQ节点:给定节点注入有功功率和无功功率。

3)平衡节点,又称为参考节点,即为全网电压幅值和相角的参考点,已知电压幅值,相角为0度。

依据以上标准,本文的测试系统中,节点1,2,3为PQ节点;节点4为PV节点;节点5为平衡节点。

③系统在标幺值下的等值电路

在对电力系统做仿真研究时,常常要将其参数表示成标幺值。而要对给定系统进行潮流仿真计算以及静态安全分析,则必须要将其标幺值下的等值电路给出。电力系统在标幺值下的等值电路实际上是由系统中各个元件的标幺值构成的,下面将对系统各元件的处理进行阐述。

1)发电机模型

由于静态安全分析研究的电路为稳态电路,则发电机模型为理想电源。 2)线路模型选取

线路模型依距离考虑三种线路模型。 a. ?一?型模型,适用于小于100km架空线。

b. 中距离的?π?型和?T?模型,适用于100~300km架空线,小于100km的电缆。 c. 长距离的?π?型,考虑波过程,适用于大于300km的架空线,大于100km的电缆。

因此,依据本文具体数据,线路模型处理如下:支路1-2有电纳、阻抗,用π型等值电路;支路1-3只有阻抗,用“一”型等值电路;支路2-3有阻抗、电纳,用π型等值电路。

18

攀枝花学院本科毕业设计(论文) 3仿真分析

3)变压器模型

变压器模型分为以下3种。 a. Τ型等值模型

Τ型等值模型为变压器基本模型,考虑磁路和电路原理,推导出的电路模型。 b. Γ型等值模型

在Τ型等值模型基础上,由于变压器正常运行情况下,励磁支路通过的电流远小于额定电流,故将变压器一次侧,二次测电抗合并,简化后的模型即为Γ型。

c. π型等值模型

在Γ型等值模型基础下,忽略励磁支路,依据功率守恒和基尔霍夫电流定律将Γ型等值模型等值为纯阻抗模型,或纯导纳模型,即将理想变压器等效了。

Γ型Τ型等值电路中,由于模型不能反映变压器实际具有的电压变换功能,需归算参数。而π型等值模型可体现电压变换,在多电压等级网络计算中采用该种变压器模型后,无须进行参数和变量的归算。所以变压器模型采用π型等值模型。在实际编程时发现,可以将π型等值模型采用理想变压器串阻抗模型,使得系统的等值电路更为简化[7]。

4)给定系统的等值电路图

按上述方法,对给定系统中的各个元件进行处理后,得到给定系统的等值电路如图3.3所示。

j0.015④1:1.05②0.08+j0.3j0.25j0.25③1.05:1j0.03⑤3.7+j1.32+j1j0.25j0.25①1.6+j0.8图3.3等值电路图

得到系统的等值电路之后,即可进行基于P-Q分解法的潮流计算仿真分析。

19

0.1+j0.250.040+j.25攀枝花学院本科毕业设计(论文) 3仿真分析

3.2.3 基于P-Q分解模型的潮流计算仿真分析

①潮流计算流程图

基于P-Q分解模型的潮流计算流程如图3.4所示。

原始数据?P?Pis?Pi?P'?1(B)VV??? 修正??(t?1)??(t)???(t)?Q?Qis?Qi?V??Q\?1?B?V 修正VV(t?1)=V(t)-?V(t)依收敛判据判断是否收敛是收敛否

图3.4基于P-Q分解模型的潮流计算流程

该流程图是根据第2.1节中基于P-Q分解模型的潮流计算步骤来绘制的。在本节,根据该流程,将采用MATLAB进行编程,进而得到基于P-Q分解模型的潮流计算结果。

②系统的电导矩阵和电纳矩阵

在原始数据中,将根据图3.3所示的等值电路,采用MATLAB语言计算可得到系统的电导矩阵以及电纳矩阵。其程序参见附录A。

对于本文的测试系统而言,其系统的电导矩阵为:

?1.3787?0.624?0.7547??0.6241.4539?0.8299?g???0.7547?0.82991.5846?000??000?

20

00?00??00?

?00?00??攀枝花学院本科毕业设计(论文) 3仿真分析

系统的电纳矩阵为:

2.641500??6.29173.9002??3.9002?66.98083.1120?63.49210??b0??2.64153.1120?35.7379031.7460?

??063.49210?66.66670???031.74600?33.3333??0?③系统正常运行状态下的节点电压

基于P-Q分解模型的潮流子程序参见附录B。通过程序的运行,发现P-Q分解法迭代13次即可得到本文测试系统的潮流计算。在正常运行状态下,系统各节点电压的幅值和相角如表3.5所示。

表3.5 系统正常运行状态下节点电压幅值和相角 节点编号 1 2 3 4 5 节点电压幅值 0.8622 1.0779 1.0364 1.0500 1.0500 节点电压相角 -0.0834 0.3116 -0.0747 0.3812 0 由表3.5可以看出,1号节点的电压低于0.95,是因为该点为纯负荷点,远离两个电源点,同时节点负荷阻抗较大,这两个因素导致1号节点电压偏低。由此也可以看出,如果要抬高1号节点的电压,则应在该节点处一应添加无功补偿装置,采用就地补偿的方式进行无功补偿,以保证系统的电压质量[13]。

3.2.4基于P-Q分解模型的静态安全分析的仿真研究

在得到本文测试系统在正常运行状态下的节点电压之后,即可采用如图2.5所示的静态安全分析的流程,针对单条线路断开的预想故障,对本文测试系统进行静态安全分析。此处再次将该流程展现出来,以便于对照程序结果进行说明,如下图3.5所示。

下面将根据该流程,基于给定的系统,详细说明基于P-Q分解模型的静态安全分析的仿真实现。

①系统正常运行状态下的数据

根据该流程,根据前一节所得到的系统正常运行状态下的潮流,再次把系统在正常运行状态下的节点信息展示出来,如表3.6所示。

表3.6的第2列和第3列给出了节点给定的注入有功功率和无功功率,在后面的静态安全分析中,模拟断线的时候,需要对节点的注入功率进行修改,此处

21

攀枝花学院本科毕业设计(论文) 3仿真分析

给出静态安全分析前的初始值。表3.6的第4列和第5列给出的是前面利用P-Q分解法潮流计算所得到的节点电压幅值和相角,该数据有利于后面静态安全分析时,与断线后系统的节点电压进行对比,以判断系统的安全状态。

稳态时潮流计算的基础数据调用稳态时支路传输功率求出模拟开断所需的节点注入功率修改稳态时节点注入功率潮流计算子函数支路传输功率子函数在主函数中与上下限比较分析结果:即静态安全分析结果 图3.5静态安全分析流程图

表3.11 系统正常运行状态下的节点信息

节点 1 2 3 4 5 节点注入有功 -1.6 -2 -3.7 5 0 节点注入无功 -0.8 -1 -1.3 0 0 电压幅值 电压相角 0.86215 1.077916 1.036411 1.05 1.05 -0.0834 0.311603 -0.07473 0.381238 0 ②系统正常运行状态支路传输功率

在得到系统的各个节点的电压之后,即可求取系统在正常运行状态下各支路传输的功率。其程序参见附录C。利用该程序,得到系统在正常运行状态下各支路传输的功率如表3.6所示。根据第2章式(2.67)和式(2.68),得到系统正常运行状态下支路的传输功率,有助于求解模拟断线时节点的注入功率增量。

表3.6 系统正常运行状态下支路传输功率 支路传输功率 支路1-2 支路1-3 支路2-3 支路2-4 支路3-5 -5 5 -2.57943 2.579427 支路首端传输的有功功率 0.707755 -0.54637 -1.9249 支路末端传输的有功功率 -1.58455 -0.15679 1.27736 支路首端传输的无功功率 0.037425 0.205098 -0.04614 1.686916 -0.5494

22

攀枝花学院本科毕业设计(论文) 3仿真分析

支路末端传输的无功功率 -0.96303 -0.73985 -0.74024 -2.08471 0.082142 ③ 模拟断线时节点的注入功率增量

模拟断线时,求取节点的注入功率增量的程序见附录D。其程序运行结果如表3.7所示。

表3.7 模拟断线时节点的注入功率增量 节点注入功率增量 断线支路首节点有功注入功率增量 断线支路末节点有功注入功率增量 断线支路首节点无功注入功率增量 断线支路末节点无功注入功率增量 支路1-2断开 0.679078 -1.55587 0.108261 -1.0409 支路1-3断开 支路2-3断开 -0.54634 -0.15682 0.317586 -0.85592 -1.92518 1.27764 0.026242 -0.81263 由于如图3.2所示的系统接线图中,有3条线路支路,有2条变压器支路。在本文中,对于单条支路的断开,仅针对于线路支路。因此,对于线路支路1-2,1-3,2-3来说,在模拟断线时,其节点注入功率是要发生变化的。而对于变压器支路2-4,3-5来说,其连着发电机,此时其分析方法与本文所研究的分析方法有很大区别,因此本文未考虑其开断情况。

④ 修改节点注入功率

在得到单条线路断线时相应节点的注入功率增量后,可对系统正常运行状态下节点的注入功率进行修改,其程序见附录E。其运行结果如表3.8和表3.9所示。

表3.8 模拟线路开断时修改后的节点注入有功功率 节点注入有功功率 支路1-2开断 支路1-3开断 支路2-3开断 节点1 节点2 节点3 节点4 节点5 -0.920922 -3.55587 -3.7 5 0 -2.146343 -2 -3.856817 5 0 -1.6 -3.925177 -2.42236 5 0 表3.9 模拟线路开断时修改后的节点注入无功功率 节点注入无功功率 支路1-2开断 支路1-3开断 支路2-3开断 节点1 节点2 节点3 节点4 节点5 -0.69174 -2.0409 -1.3 0 0 -0.48241 -1 -2.15592 0 0 -0.8 -0.973758 -2.112628 0 0 表3.8表示模拟线路断开时修改后的节点注入有功功率,表3.9表示模拟线路

23

攀枝花学院本科毕业设计(论文) 3仿真分析

断开时修改后的节点注入无功功率。由于只研究线路开断,所以与线路支路相连的节点注入功率会发生变化,而发电机节点的注入功率则不会发生变化。

⑤静态安全分析结果

在得到模拟断线时,各节点修改后的注入功率之后,就可以采用式(2.54)对模拟断线后各个节点的电压幅值增量进行求解,并由此可判断断线后各个节点的电压幅值是否在运行范围内。程序见附录F。并得到安全分析结果如表3.10所示。

表3.10 针对节点电压的静态安全分析及结果 节点号 1 2 3 4 5 电压 上限 0.95 0.95 0.95 0.95 0.95 电压下限 1-2开断 1.05 1.05 1.05 1.05 1.05 0.95761 1.0713 1.0532 1.05 1.05 支路断开后的节点电压幅值及其越限情况 结果 安全 越上限 越上限 安全 安全 1-3开断 结果 2-3开断 结果 0.858136 越下限 0.87159 越下限 1.075916 越上限 1.078678 越上限 1.015741 1.05 1.05 安全 安全 安全 1.03253 1.05 1.05 安全 安全 安全 由表3.10可以知道,当支路1-2断开的时候,节点2和节点3的电压幅值越限;当支路1-3断开的时候,节点1和节点2的电压幅值越限;当支路2-3断开的时候,节点1和节点2的电压幅值越限。这个分析结果说明,在这三条线路支路分别断开的情况下,都需要进一步给出安全稳定控制对策,以保证系统的安全稳定运行。

24

攀枝花学院本科毕业设计(论文) 结论

结 论

本文主要对基于P-Q分解模型的静态安全分析方法及其在MATLAB中的实现进行了深入研究。其所做的主要工作如下:

①研究了基于P-Q分解模型的潮流计算原理

基于P-Q分解模型的静态安全分析其基本的思想来源于基于P-Q分解模型的潮流计算,因此,本文首先对基于P-Q分解模型的潮流计算原理进行了研究,从而为后续静态安全分析中判断系统的安全状态奠定了分析基础。

②研究了基于P-Q分解模型的静态安全分析

在研究基于P-Q分解模型的潮流计算原理的基础上,针对系统的N-1断线故障,研究了基于P-Q分解模型的静态安全分析方法。明确了如何去判断系统断线后的安全状态。

③在MATLAB中的实现了基于P-Q分解模型的静态安全分析仿真研究 针对给定的5节点环网系统,得到其等值电路之后,首先实现了该系统基于P-Q分解模型的潮流仿真计算,得到了系统在正常运行状态下的参数。而后,针对单条线路断开的预想故障,实现了基于P-Q分解模型的静态安全分析仿真研究,对系统断线后的状态进行了判断。

通过本文的研究,本人体会如下: ①关于程序编写的总结

编程前,对问题都有一个可行性分析的过程,前期的理论推导,证明了理论上是可行。下一步就是模块化的解决,将具体问题分成一个一个小的问题。当逐个小的问题解决时,大的问题也就解决了。而前期的可行性分析,中期的用编程语言在头脑中尝试实现各个模块,将认为成熟的代码,写在纸上,在纸上进行运行分析。当两步都实现了,再在编程环境中检验,即编程环境的作用仅作为检验调试程序的工具。

②关于编写时考虑的问题

编写程序时考虑的因素有机时、内存、本文具体数据。同时需要考虑以下问题:

1)跨文件读取变量的三种方式

在考虑将基础数据模块化时,考虑用从文件读取,应用全局变量和函数调用三种方式。

从硬盘文件读取数据的形式适合于电网,即多节点,大数据量,这种方式让平时运行的程序在运行中实时使用的变量较少,则占用内存少。变量只会在需要

25

攀枝花学院本科毕业设计(论文) 结论

时从文件中调去,不用时就清空工作空间,节省运行内存。

由于全局变量的工作空间在整个程序的运行过程作中始终存在。文中节点虽少,如果采用全局变量方式读取数据则需要设p,q,voltage,ang,b0,b1,b2,g八个全局变量,且全局变量的存在易让程序结构混乱,所以未采用全局变量。

以子函数调用的形式读取变量时适合本文,因其节点少,只有5个,读取内存的速度远快于读取硬盘文件,但应明确此种读取数据模式只针对本文,如果运算设备有充足的物理内存,这种模式有利于系统的响应时间的提高,是应该被采用的。并且子函数调用得到的返回值只能对稳态数据的读取,不能对源文件中的原始数据进行修改,因此有以下好处:

a.由于静态稳定只用稳态时的数据,则相对3次静态安全分析,无需担心调用数据对整个初始稳态值的修改,全局变量参与运算则会修改线路开断前得到的稳态解,而需增加数据恢复的步骤。

b.考虑到每次调用函数耗费的机时远大于从内存读取数据,所以基础数据函数只宜调用一次,函数返回值则应保存在静态分析的主程序中,主程序不结束,其对应的工作空间不会消失并且基础数据的变量名可用于主程序中。

2)考虑电路模型

在考虑电路模型时,编写导纳矩阵程序的时候发现运用理想变压器模型和导纳的定义去求解导纳矩阵和运用π型等值模型效果相同,因π型等值电路模型源于对理想变压器模型的等值变换。而采用定义求解导纳矩阵更加直接。在程序算法总采用的是变压器理想模型,在画电路模型图时依然采用理想变压器模型。

3)数据结构的优化

对数据结构的优化可以提高程序的运行效率,本文只是一个5节点的小程序,对其进行数据结构的优化再电脑反应上暂时没有明显的差异,但对于真实的电网合理优化节点应该是提高运行效率的一种方法。

4)方案思考

P-Q分解模型在众多模型中是结合电网实际对Newton-Raphson法进行了简化的,模型本身是优秀的。由于每个程序员的个人思路的不同,程序的运行效率依然存在很大差异。静态安全分析中本文中采用的是用潮流程序对修改后的节点功率网络再次进行迭代,计算出相应断线情况下的节点状态参数:电压幅值和相角。三次静态分析所需的迭代次数分别为12次,12次,10次。另一种方式是利用公式通过功率变化值,直接求出角度和电压的变化值,最后对电压幅值和相角进行修正,这种运行方式无需调用潮流子函数,则无需迭代,3次支路开断模拟已经少了34次的迭代,这对程序消耗机时的减少是很有效的。所以,本文中运用了优秀的潮流分析方法,因学生思索深度不够,未能及时提出认识到新方案。

26

攀枝花学院本科毕业设计(论文) 结论

5)对用注入功率模拟线路开断的方法的思考

本文中用注入功率模拟线路开断,相对于变导纳矩阵的方法固然是中进步,但关于用注入功率模拟线路开断的方法是否能够在理论上做出真实模拟,提出疑问 :一条线路开断将引起全网功率的分布的改变,而用注入功率模拟线路的开断只是对被模拟开断线路两侧节点进行了节点注入功率的修改,理论上并未对线路开断进行真实的模拟。对于该问题,在后续的研究中,初步拟定具体测试步骤如下:首先,在线路参数中删除一条线路的参数,即该线路不存在,即断开了,此时运行潮流程序对各个节点的电压幅值相角求解;然后,用修改注入功率模拟线路开断,两种情况下比对即可得出结论。

6)程序的适用性

关于程序的适用性上,此程序只限于5节点数据,是由于其初始化导纳矩阵维数的定义上的限制,如果对程序中的变量维数进行优化,维数用一具体变量替换,则该程序的适用性将不受节点个数的限制。此外由于时间和基础原因,程序运行结果不是友好界面,只是输出最终的运行数据和判断结果,即应用中至少应该是图形界面的。这是可以提高的。最后一点思考到的不足是,没有预留实时读取数据的接口函数。基础数据为手动输入数据,在真实开断模拟中,应实时跟踪采集电网数据,周期性的进行静态安全分,实时显示可能故障的支路,已供调度人员参考。

7)静态安全分析方法的优点

在本文的静态安全分析方法中,潮流迭代算法只用在第一次求系统正常运行状态参数上。随后的静态安全分程序则应用公式求出电压和相角的修正量,直接获得开断后节点的电压和相角。如果针对本文测试系统,在每一个断线预想故障状态下,都进行一次潮流计算的话,则会增加34次的迭代运算。因此,本文所研究的静态安全分析方法少了34次的迭代运算,这是在程序算法上的提高,运行效率提升很明显,也充分体现了本文所研究的静态安全分析方法的优点。

27

攀枝花学院本科毕业设计(论文) 参考文献

参 考 文 献

[1]KUNDUR. P. Power System Stability and Control [M]. 北京:中国电力出版社,2002. [2]吴际舜.电力系统静态安全分析[M].上海交通大学出版社,1985. [3]王锡凡, 方万良, 杜正春.现代电力系统分析[M].科学出版社,2003.

[4]吴天鹏,谢小介,彭斌. MATLAB电力系统设计与分析[M]. 北京:国防工业出版社,2004. [5]王正林,刘明. 精通MATLAB7[M].电子工业出版社,2006. [6]何仰赞.电力系统分析[M].武汉:武汉华中科技大学出版社,1985. [7]于永源.电力系统分析[M].北京:中国电力出版社,2007.

[8]陈珩.电力系统稳态分析(第三版)[M]. 北京:中国电力出版社,2007. [9]李光琦. 电力系统暂态分析[M].北京:中国电力出版社,1995.

[10]李文沅. 电力系统风险评估模型-方法和应用[M].北京:科学出版社,2006. [11]李文沅. 电力系统经济安全运行模型与方法[M].重庆:重庆大学出版社,1989. [12]周荣光.电力系统故障分析[M].北京:北京清华大学出版社,1988. [13]中华人民共和国国家经济贸易委员会. 电力系统安全稳定导则DL755-2001. [14]许高达.电力系统之计算机应用[M].新兴图书公司,1979. [15]朱仁峰. 精通MATLAB7[M].清华大学出版社,2006.

[16]Math Works. 产品和服务[EB/OL]. http://www.mathworks.cn/products/matlab/,2014-5-13. [17]刘浏. 百度百科[EB/OL].http:// baike.http://www.wodefanwen.com//view/10598.htm?fr=aladdin,2014-5-13.

28

攀枝花学院本科毕业设计(论文) 附录A:基础数据生成子程序

附录A:基础数据生成子程序

程序_基础数据生成子程序

function[p,q,voltage,ang,g,b0,b1,b2,y0,node_f,node_e]=basic_data();

%首节点 末节点 支路阻抗 对地导纳的1/2 变比

Z=[

1 2 0.04+j*0.25 j*0.25 1 3 0.1+j*0.35 0 2 3 0.08+j*0.3 j*0.25 4 2 j*0.015 0 5 3 j*0.03 0 ];

admittance=zeros(5); %初始化导纳矩阵 for m=1:5

p=Z(m,1); %p=1,1,2,4,5 q=Z(m,2); %q=2,3,3,2,3

admittance(p,q)=admittance(p,q)-1./(Z(m,3)*Z(m,5));

admittance(q,p)=admittance(p,q); %互导纳相等 %计算自导纳

admittance(q,q)=admittance(q,q)+1./(Z(m,3)*Z(m,5)^2)+Z(m,4); %是否添加对地导纳参数%计算高压侧自导纳 admittance(p,p)=admittance(p,p)+1./Z(m,3)+Z(m,4); %是否添加对地导纳参数%计算低压侧自导纳 end

conductance=real(admittance); %电导英文:conductance susceptance=imag(admittance); %电纳英文:susceptance b1=susceptance(1:4,1:4); %4*4

b1=inv(b1);%P-Q分解法中对应p中的电纳矩阵//仅用作求逆 b2=susceptance(1:3,1:3);

b2=inv(b2);%P-Q分解法中对应q的电纳矩阵 b0=susceptance;

29

1 1.05 1.05 1

1

攀枝花学院本科毕业设计(论文) 附录A:基础数据生成子程序

g=conductance;

var=[%发电机功率 负荷功率 电压初值 节点类型

0 1.6-j*0.8 1 2 0 2-j 1 2 0 3.7-j*1.3 1 2 5 0 1.05 3 0 0 1.05 1 ];

%/1平衡节点,2P-Q节点,3pv节点 voltage=zeros(5,1); %创建电压列向量 v=zeros(4,1);

ang =zeros(5,1); %创建相角列向量 voltage=var(:,3); %对电压列向量初始化

for m=1:5

s(m)=var(m,1)-var(m,2); end s=s'; p=real(s); q=imag(s);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%计算支路潮流所需的参数

y0 =Z(:,4); %取对地导纳矩阵 y0 =imag(y0); %对地导纳的共轭矩阵 node_f=zeros(5,1); node_e=zeros(5,1); node_f=Z(:,1); node_e=Z(:,2);

[voltage,ang]=P-Q(voltage,ang,p,q,g,b0,b1,b2); Return

30

攀枝花学院本科毕业设计(论文) 附录B:潮流计算子程序

附录B:潮流计算子程序

%,输入voltage,ang,p,q,g,b0,b1,b2,输出voltage,ang function [voltage,ang]=P-Q(voltage,ang,p,q,g,b0,b1,b2) %volatge,ang,为迭代计算的初值; iteration=1;

tol=1e-6; %%收敛容许误差 while(iteration<20) v=voltage(1:4,1); p_input=zeros(4,1);

for m=1:4 for n=1:5

p_input(m)=voltage(n)*(g(m,n)*cos(ang(m)-ang(n))+b0(m,n)*sin(ang(m)-ang(n)))+p_input(m);

end

p_input(m)=p_input(m)*voltage(m); end

p_correction=zeros(4,1); p_correction=p(1:4,1)-p_input;

ang_correction=zeros(4,1);

ang_correction=b1*p_correction./v./v;

ang(1:4,1)=ang(1:4,1)-ang_correction;

q_input=zeros(3,1); for m=1:3 for n=1:5

q_input(m,1)=voltage(n,1)*(g(m,n)*sin(ang(m)-ang(n))-b0(m,n)*cos(ang(m)-ang(n)))

31

攀枝花学院本科毕业设计(论文) 附录B:潮流计算子程序

+q_input(m);

end

q_input(m,1)=voltage(m,1)*q_input(m,1); end

q_correction=zeros(3,1); q_correction=q(1:3,1)-q_input; %% check for convergence F = [ p_correction;

q_correction ]; normF = norm(F, inf);

if normF < tol

fprintf('\\nNewton''s method power iterations.\\n', iteration);

return;%跳出对电压相角的计算 end

v2=v(1:3,1);

v_correction=zeros(3,1);

v_correction=b2*q_correction./v2;

voltage(1:3,1)=voltage(1:3,1)-v_correction; iteration=iteration+1; end

32

flow converged in %d 攀枝花学院本科毕业设计(论文) 附录C:支路功率计算子程序

附录C:支路传输功率计算子程序

%功能求支路流过的功率

%注明:1)pij和pji在数值上是不同的。2)默认i为z矩阵中的第一列元素,为低压侧数据,即默认branch_p中第一列总是存放:从低压侧到高压侧的支路功率

%以上的默认是对数据的结构体化,对数据的结构体一样固定其格式,便于对程序的优化,从中受益已不止一次。

function [branch_p,branch_q]=branch_s(v,ang,g,b0,y0,node_f,node_e) for m=1:5 %pij i=node_f, j=node_f

%为程序的可读性,设置如下3个中间变量;为方便模拟短线注入功率的求解将pij和qij设置为2*5的格式

n1=node_f(m,1); n2=node_e(m,1);

ang_mid=ang(n1,1)-ang(n2,1); %pij

branch_p(1,m)= v(n1,1)^2*g(n1,n2)-v(n1,1)*v(n2,1)*(g(n1,n2)*sin(ang_mid)+b0 (n1,n2)*sin(ang_mid));

%此处由于网络拓扑数据的特殊,没有电纳,故此处公式为简化后的公式 %qij

branch_q(1,m)=-v(n1,1)^2*(y0(n1,1)+b0(n1,n2))-v(n1,1)*v(n2,1)*(g(n1,n2)*sin (ang_mid)-b0(n1,n2)*cos(ang_mid));

%pji

branch_p(2,m)=v(n2,1)^2*g(n2,n1)-v(n2,1)*v(n1,1)*(g(n1,n2)*cos(-ang_mid)+b0

(n2,n1)*sin(-ang_mid));

%qji

branch_q(2,m)=-v(n2,1)^2*(y0(n2,1)+b0(n2,n1))-v(n2,1)*v(n1,1)*(g(n2,n1)*sin(-

ang_mid)-b0(n2,n1)*cos(-ang_mid));

end

33

攀枝花学院本科毕业设计(论文) 附录C:支路功率计算子程序

end

34

攀枝花学院本科毕业设计(论文) 附录D:支路开断时首末节点注入功率增量计算子函数

附录D:支路开断时首末节点注入功率增量计算子函数

function

[p_inject,q_inject]=simulink_dynamic(p,q,branch_p,branch_q,b0,b1,b2,y0,node_f,node_e)

p_inject=zeros(2,5); q_inject=zeros(2,5); for m=1:3

n1=node_f(m,1); n2=node_e(m,1);

p_inject(:,m)=inv(eye(2)-[-1,1;1,-1]*[b1(n1,n2),b1(n1,n2);b1(n2,n1),b1(n2,n2)]*b1(n1,n2))*branch_p(:,m);

c=1+2*y0(n1,1)/b0(n1,n2);

q_inject(:,m)=inv(eye(2)-[-c,1;1,-c]*[b2(n1,n1),b2(n1,n2);b2(n2,n1),b2(n2,n2)*b0(n1,n2)])*branch_q(:,m);

end

end

35

攀枝花学院本科毕业设计(论文) 附录E:修改节点注入功率子程序

附录E:修改节点注入功率子程序

%输入稳态时,节点注入功率, %输出断线后,节点注入功率

function[p_revise,q_revise]=pq_revise(p,q,p_inject,q_inject,node_f,node_e) p0=p; q0=q; for m=1:3

n1=node_f(m); n2=node_e(m);

p(n1,1)=p(n1,1)+p_inject(1,m); p(n2,1)=p(n2,1)+p_inject(2,m); p_revise(:,m)=p;

q(n1,1)=q(n1,1)+q_inject(1,m); q(n2,1)=q(n2,1)+q_inject(2,m); q_revise(:,m)=q; p=p0; q=q0; end end

36

本文来源:https://www.bwwdw.com/article/l663.html

Top