基于matlab的电力系统潮流计算毕业论文.doc
- 文档编号:359633
- 上传时间:2023-04-29
- 格式:DOC
- 页数:50
- 大小:1.64MB
基于matlab的电力系统潮流计算毕业论文.doc
《基于matlab的电力系统潮流计算毕业论文.doc》由会员分享,可在线阅读,更多相关《基于matlab的电力系统潮流计算毕业论文.doc(50页珍藏版)》请在冰点文库上搜索。
1引言
1.1本课题的目的和意义
电力系统潮流计算是对复杂电力系统正常和故障条件下稳态运行状态的计算。
其目的是求取电力系统在给定运行方式下的节点电压和功率分布,用以检查系统各元件是否过负荷、各点电压是否满足要求、功率分布和分配是否合理以及功率损耗等,是电力系统计算分析中的一种最基本的计算[1]。
潮流计算是电力系统的各种计算的基础,同时它又是研究电力系统的一项重要分析功能,是进行故障计算,继电保护鉴定,安全分析的工具。
电力系统潮流计算是计算系统动态稳定和静态稳定的基础。
在电力系统规划设计和现有电力系统运行方式的研究中,都需要利用电力系统潮流计算来定量的比较供电方案或运行方式的合理性、可靠性和经济性[1]。
对于正在规划的电力系统,通过潮流计算,可以为选择电网供电方案和电气设备提供依据。
潮流计算还可以为继电保护和自动装置整定计算、电力系统故障计算和稳定计算等提供原始数据。
潮流计算的目的在于:
确定是电力系统的运行方式;检查系统中的各元件是否过压或过载;为电力系统继电保护的整定提供依据;为电力系统的稳定计算提供初值,为电力系统规划和经济运行提供分析的基础。
因此,电力系统潮流计算是电力系统中一项最基本的计算,既具有一定的独立性,又是研究其他问题的基础[1]。
1.2国内外发展现状
利用电子计算机进行潮流计算从20世纪50年代中期就已经开始。
此后,潮流计算曾采用了各种不同的方法,这些方法的发展主要是围绕着对潮流计算的一些基本要求进行的。
对潮流计算的要求可以归纳为下面几点:
(1)算法的可靠性或收敛性
(2)计算速度和内存占用量
(3)计算的方便性和灵活性
电力系统潮流计算属于稳态分析范畴,不涉及系统元件的动态特性和过渡过程。
因此其数学模型不包含微分方程,是一组高阶非线性方程。
非线性代数方程组的解法离不开迭代,因此,潮流计算方法首先要求它是能可靠的收敛,并给出正确答案。
随着电力系统规模的不断扩大,潮流问题的方程式阶数越来越高,目前已达到几千阶甚至上万阶,对这样规模的方程式并不是采用任何数学方法都能保证给出正确答案的。
这种情况促使电力系统的研究人员不断寻求新的更可靠的计算方法[2]。
1.2.1高斯-赛德尔迭代法
在用数字计算机求解电力系统潮流问题的开始阶段,人们普遍采用以节点导纳矩阵为基础的高斯-赛德尔迭代法(一下简称导纳法)。
这个方法的原理比较简单,要求的数字计算机的内存量也比较小,适应当时的电子数字计算机制作水平和电力系统理论水平,于是电力系统计算人员转向以阻抗矩阵为主的逐次代入法(以下简称阻抗法)。
20世纪60年代初,数字计算机已经发展到第二代,计算机的内存和计算速度发生了很大的飞跃,从而为阻抗法的采用创造了条件。
阻抗矩阵是满矩阵,阻抗法要求计算机储存表征系统接线和参数的阻抗矩阵。
这就需要较大的内存量。
而且阻抗法每迭代一次都要求顺次取阻抗矩阵中的每一个元素进行计算,因此,每次迭代的计算量很大[3]。
阻抗法改善了电力系统潮流计算问题的收敛性,解决了导纳法无法解决的一些系统的潮流计算,在当时获得了广泛的应用,曾为我国电力系统设计、运行和研究作出了很大的贡献。
但是,阻抗法的主要缺点就是占用计算机的内存很大,每次迭代的计算量很大。
当系统不断扩大时,这些缺点就更加突出。
为了克服阻抗法在内存和速度方面的缺点,后来发展了以阻抗矩阵为基础的分块阻抗法。
这个方法把一个大系统分割为几个小的地区系统,在计算机内只需存储各个地区系统的阻抗矩阵及它们之间的联络线的阻抗,这样不仅大幅度的节省了内存容量,同时也提高了计算速度[4]。
1.2.2牛顿-拉夫逊法和P-Q分解法
克服阻抗法缺点的另一途径是采用牛顿-拉夫逊法(以下简称牛拉法)。
牛拉法是数学中求解非线性方程式的典型方法,有较好的收敛性。
解决电力系统潮流计算问题是以导纳矩阵为基础的,因此,只要在迭代过程中尽可能保持方程式系数矩阵的稀疏性,就可以大大提高牛顿潮流程序的计算效率。
自从20世纪60年代中期采用了最佳顺序消去法以后,牛拉法在收敛性、内存要求、计算速度方面都超过了阻抗法,成为直到目前仍被广泛采用的方法。
在牛拉法的基础上,根据电力系统的特点,抓住主要矛盾,对纯数学的牛拉法进行了改造,得到了P-Q分解法。
P-Q分解法在计算速度方面有显著的提高,迅速得到了推广[5]。
牛拉法的特点是将非线性方程线性化。
20世纪70年代后期,有人提出采用更精确的模型,即将泰勒级数的高阶项也包括进来,希望以此提高算法的性能,这便产生了保留非线性的潮流算法。
另外,为了解决病态潮流计算,出现了将潮流计算表示为一个无约束非线性规划问题的模型,即非线性规划潮流算法[6]。
近20多年来,潮流算法的研究仍然非常活跃,但是大多数研究都是围绕改进牛拉法和P-Q分解法进行的。
此外,随着人工智能理论的发展,遗传算法、人工神经网络、模糊算法也逐渐被引入潮流计算。
但是,到目前为止这些新的模型和算法还不能取代牛拉法和P-Q分解法的地位。
由于电力系统规模的不断扩大,对计算速度的要求不断提高,计算机的并行计算技术也将在潮流计算中得到广泛的应用,成为重要的研究领域[7]。
通过几十年的发展,潮流算法日趋成熟。
近几年,对潮流算法的研究仍然是如何改善传统的潮流算法,即高斯-塞德尔法、牛拉法和快速解耦法。
牛拉法,由于其在求解非线性潮流方程时采用的是逐次线性化的方法,为了进一步提高算法的收敛性和计算速度,人们考虑采用将泰勒级数的高阶项或非线性项也考虑进来,于是产生了二阶潮流算法。
后来又提出了根据直角坐标形式的潮流方程是一个二次代数方程的特点,提出了采用直角坐标的保留非线性快速潮流算法[8]。
1.2.3基于MATLAB的电力系统潮流计算发展前景
MATLAB自1980年问世以来,以其学习简单、使用方便以及其它高级语言所无可比拟的强大的矩阵处理功能越来越受到世人的关注。
目前,它已成为国际控制界最流行、使用最广泛的语言了。
它的强大的矩阵处理功能给电力系统的分析、计算带来许多方便。
在处理潮流计算时,其计算机软件的速度已无法满足大电网模拟和实时控制的仿真要求,而高效的潮流问题相关软件的研究已成为大规模电力系统仿真计算的关键。
随着计算机技术的不断发展和成熟,对MATLAB潮流计算的研究为快速、详细地解决大电网的计算问题开辟了新思路。
MATLAB语言允许用户以数学形式的语言编写程序,其比BASIC语言和FORTRAN等更为接近书写的数学表达格式,且程序易于调试。
在计算要求相同的情况下,使用MATLAB编程,工作量将会大为减少[9]。
基于MATLAB的电力系统潮流计算使计算机在计算、分析、研究复杂的电力系统潮流分布问题上又前进了一步。
矩阵输入、输出格式简单,与数学书写格式相似;以双精度类型进行数据的存储和运算,数据精确度高,能进行潮流计算中的各种矩阵运算,包括求逆、求积和矩阵LR分解等,其程序的编写也因MATLAB提供了许多功能函数而变得简单易行。
另外,MATLAB稀疏矩阵技术的引入,使电力系统潮流计算由传统方法转变为优化算法成为可能[10]。
2简单电力系统潮流计算的手工方法
2.1简单辐射网络的潮流计算
大约半个多世纪以前,数字计算机还没有出现的时候,潮流计算都是采用手工的计算方法。
虽然潮流计算的本质是解电力系统的节点功率方程,然而手工的计算方法是不可能用解上述节点功率方程的方法来进行潮流计算的。
手工潮流计算是根据一个简单支路的电压和功率传输关系,将较为复杂的电力系统分解为若干个简单支路来进行潮流计算的。
因此任何复杂的潮流计算都是从一个简单支路的潮流分布和电压降落的计算开始的。
2.1.1简单支路的潮流分布和电压降落
如图1所示的简单支路,节点1和2之间的阻抗为已知;两端的电压分别为和,从节点1注入该支路的复功率为,从节点2流出的功率为,阻抗消耗的功率为。
根据电路理论,、和、这四个变量,任何两个变量已知都可以求出另外两个变量。
图2.1简单支路示意图
(1)已知一侧的电压和功率求另一侧的电压和功率
假设已知节点2的电压和流出的功率,可知道流过该支路的电流为:
式(2.1)
如果以作为参考相量,阻抗Z引起的电压降落和功率损耗分别为:
式(2.2)
式(2.3)
因此另一端节点1的电压为:
式(2.4)
流过节点1的复功率为:
式(2.5)
两端电压的关系还可以从如图2所示的相量图中得到(以为参考相量),为末端电压和电流的夹角,称为功率因数角。
从相量图中,不难得到阻抗Z引起的电压降落的横分量和纵分量分别为:
式(2.6)
可得到首端的电压幅值和相角分别为:
式(2.7)
式(2.8)
j
d
j
j
j
如果已知首端(节点1)的电压和功率,求末端的电压和功率,其基本原理同上,读者可以自行推导分析。
图2.2两端电压相量示意图
(2)已知一端的电压和流过另一端的复功率
假如已知首端电压和末端的功率,要求首端的功率和末端的电压,我们可以利用两端电压的关系以及两端功率的关系列出如下方程组(以为参考相量):
式(2.9)
式(2.10)
直接求解上面这个相量方程组是很麻烦的,可以通过迭代法来求解:
先给定一个末端电压的初值,这个初值可以设定为该节点的平均额定电压,然后将之代入2.9,得到,然后再利用根据2.10得到,重复上面的过程,直到误差满足要求为止。
由于潮流计算通常是在电力系统的稳态运行条件下,此时节点电压与平均额定电压差别不大,因此,在手工近似计算中,将上述的迭代过程只进行一次。
即先设定未知的电压为平均额定电压,利用2.3式,根据末端的功率计算支路的功率损耗,然后利用2.5式计算出首端的功率,再利用首端的功率和首端的电压计算系统的电压损耗,最后计算出末端的电压。
2.1.2辐射型网络的手工潮流计算方法
所谓辐射型网络就是单电源供电的非环形网络,系统中所有的负荷都由一个电源供电,辐射型网络是由若干个简单支路树枝状串级联接而成的。
对于辐射型网络中的接地支路可以做如下处理:
(1)将对电力系统中的接地支路等效为该支路消耗的功率,对地支路的电压用额定电压来替代,例如,对地支路的导纳为,那么这个对地支路的消耗的功率;
(2)将同一节点消耗的功率进行合并。
通过这样处理,辐射型网络就化减为若干简单支路的级联,可以利用简单支路的潮流和电压计算方法逐级进行潮流计算。
辐射型网络的手工潮流计算一般从系统末端开始,因为通常辐射型网络的末端的负荷为已知,首先计算潮流的近似分布,然后再从电源端开始根据潮流分布计算出各个节点的电压。
因此,辐射型网络的手动潮流估算仅包含三步:
第一步,根据电力系统各个元件的电机参数,建立电力系统的等值计算电路;然后将对地支路等效为支路消耗的功率,并将各个节点消耗的功率进行合并。
第二步,首先将系统中各个节点的未知电压设为系统平均额定电压,然后从辐射型网络的末端开始,依次计算各个支路的功率损耗,最后得到潮流在辐射型网络中的近似分布。
第三步,根据估算出的潮流分布,从电源端开始,根据前面简单支路的电压计算公式依次计算各个节点的电压。
通过一个实例来说明潮流计算的过程,如图3所示的辐射型单电源的简单电力系统,已知节点1(发电机节点)的电压和各个节点的负荷、、、,求该系统的功率和电压的分布。
图2.3单电源辐射型电力系统
已知电力系统的各个元件的参数如下所示:
变压器T1:
额定容量,额定变比,空载损耗,空载电流百分数,短路损耗,短路电压百分数;
输电线路L:
每公里长的正序阻抗,每公里长的对地电纳,线路长度;
变压器T2:
额定容量,额定变比,空载损耗,空载电流百分数,短路损耗,短路电压百分数。
第一步作出等效电路及其参数:
首先做电力系统的等值电路,根据上述各个元件的参数,我们可以得到各个元件的等效电路及其电路参数,等效电路如图2.4所示。
在计算等值电路中各个元件参数之前,先选择功率和电压的基准值,,,。
变压器T1(根据等值电路,变压器参数都归算到高压侧):
;;;
;;;式(2.11)
输电线路:
;式(2.12)
变压器T2(根据等值电路,变压器参数都归算到高压侧):
;;
;;;式(2.13)
图2.4等值电路I
第二步,将对地支路简化为对地功率损耗:
如果电压基准值的选取与变压器的实际变比相匹配,那么,如果不匹配,则需要将变压器的变比的标么值等效到电路中,把变压器的阻抗支路,变为PI型等效电路。
为了说明问题,我们假设电压基准值选取与变压器实际变比匹配,或者忽略非标准变比的影响。
对地支路假设为对地损耗功率,其对地支路的损耗用该点的额定电压来计算,等效电路变为如图2.5所示。
图2.5等值电路II
其中:
;;
;;
第三步,节点功率合并:
然后,将1、2、3、4各个节点上的所有功率合并,如图2.6所示:
图2.6等值电路III
其中:
;;;。
第四步,从末端开始,根据末端功率计算功率分布:
先用各个节点的额定电压以及流出支路的功率来计算各个支路损耗以及功率分布:
;;
;;
;;
这样,就求得了功率的分布和节点1的注入功率。
第五步,从首端开始,根据首端电压计算电压损耗和各个节点的电压:
;;
;
;
2.2简单环形网络的潮流计算
环形网可以等效成两端供电网,两端供电网也可以等效成环形网。
2.2.1两端电压相等
如图下图所示、可将(a)图等效成(b)图。
S2
Z12
Z13
Z13
S3
Sa
Sb
(a)
Sa
S2
S3
Sb
Z12
Z23
Z31
1
2
3
(b)
1
ˊ
U1
·
U1=
ˊ
U1
·
·
图2.7简单环形网等效两端供电网
(a)环形网(b)两端供电网
1
2
·3
·
·
·
·
·
·
·
·
式(2.14)
2.2.2两端电压不相等
两端电压不相等的网络,可以等效成回路电压不为零的单一环网。
←
Z34
S2
S3
Z12
Sa
U
·
du
Z23
Sc
U1
·
U4
·
(b)
·
·
·
·
Sa
S2
S3
Z12
Z23
Z31
1
2
3
4
U1
·
U4≠U1
·
·
(a)
·
·
·
图2.8两端电压不等的网与环网等值
(a)两端供电网(b)环形网
式(2.15)
其中
式(2.16)
称为循环功率。
对环形网的潮流分布,首先求出、,然后求各支路上的流动功率,即初步的潮流分布,没有计及网络各段的电压降落、功率损耗。
初步潮流分布的目的,在于找出功率分点,以便在功率分点把闭环网打开成两个辐射网。
然后,以功率分点为末端,对这两个辐射网分别用逐段推算法进行潮流分布计算。
从中要计及各段的电压降落和功率损耗,所运用的公式与计算辐射网时完全相同。
在两端供电网中,当两端电压相量不等,不论是模值还是相位不等都将产生循环功率。
在环网中,循环功率是由于环网中有多台变压器,而变压器的变比不匹配引起的。
所谓变比不匹配则是指环网中有两台及以上变压器时,由于变压器变比的不同使得网络空载且开环时开口两侧有电压差,即开口两侧感应电势不同,因而闭环后,即使空载也有环路电流,产生循环功率。
应该特别注意正确地确定环网中循环功率的方向。
循环功率的正方向取决于电压降落的正方向。
环网和两端供电网中的循环功率可改变网络中功率的分布。
2.3手工计算算例
2.3.1网络结构图
10kV配电网络的电网结构如图所示。
已知各节点的负荷功率及线路参数如下:
Z12=1.2+j2.4Ω,Z23=1.0+j2.0Ω,Z24=1.5+j3.0Ω。
S2=0.3+j0.2MVA,S3=0.5+j0.3MVA,S4=0.2+j0.15MVA。
设母线1的电压为10.5kV,线路始端功率容许误差为0.3%。
图2.910kv配电网络
2.3.2计算各支路的功率损耗和功率分布。
假设各节点电压均为额定电压,功率损耗计算的支路顺序为3-2、4-2、2-1,第一轮计算依上列支路顺序计算各支路的功率损耗和功率分布。
则 MVA
MVA
MVA
又MVA
2.2.3求出线路各点电压,计算中忽略电压降落横分量。
第二步用已知的线路始端电压U1=10.5kV及上述求得的线路始端功率S12,按上列相反的顺序求出线路各点电压,计算中忽略电压降落横分量。
2.2.4根据上述求得的线路各点电压,重新计算各线路的功率损耗和线路始端功率
故
MVA
则 MVA
又
从而可得线路始端功率
经过两轮迭代计算,结果与第一步所得的计算结果比较相差小于0.3%,计算到此结束。
最后一次迭代结果可作为最终计算结果。
3复杂电力系统潮流计算的计算机方法
3.1潮流计算的计算机算法简介
潮流计算的计算机算法是以电网络理论为基础的,应用数值计算方法求解一组描述电力系统稳态特性的方程。
从数学上讲是一组多元的非线性方程式的求解问题,这类方程的求解过程都离不开迭代。
由于电力系统结构及参数的一些特点,同时随着电力系统不断扩大,潮流问题的方程式的阶数也越来越高,这样的非线性方程式并不是任何数学方法都能保证给出正确答案的。
这种情况就成为促使电力系统计算人员不断寻求新的且更可靠方法的一个重要因素。
电网潮流计算的性能优劣一般依据的是能否可靠收敛,计算速度的快慢,内存占有多少,使用是否方便灵活,调整和修改是否容易,是否满足工程需要等来判别,其中以是否可靠收敛作为评价的主要标准。
常用的分析法包括高斯-塞德尔法、牛顿-拉夫逊潮流算法、快速解耦算法(PQ分解法)等。
3.2潮流计算的约束条件
电力系统运行必须满足一定技术和经济上的要求。
这些要求够成了潮流问题中某些变量的约束条件,常用的约束条件如下:
3.2.1节点电压应满足:
式(3.1)
从保证电能质量和供电安全的要求来看,电力系统的所有电气设备都必须运行在额定电压附近。
PU节点电压幅值必须按上述条件给定。
因此,这一约束条件对PQ节点而言。
3.2.2节点的有功功率和无功功率应满足:
式(3.2)
PQ节点的有功功率和无功功率,以及PU节点的有功功率,在给定是就必须满足上述条件,因此,对平衡节点的P和Q以及PU节点的Q应按上述条件进行检验。
3.2.3节点之间电压的相位差应满足:
式(3.3)
为了保证系统运行的稳定性,要求某些输电线路两端的电压相位不超过一定的数值。
这一约束的主要意义就在于此。
因此,潮流计算可以归结为求解一组非线性方程组,并使其解答满足一定的约束条件。
在计算过程中,或得出结果之后用约束条件进行检验。
如果不能满足要求,则应修改某些变量的给定值,甚至修改系统的运行方式,重新进行计算。
3.3节点导纳矩阵的形成与修改
3.3.1节点电压方程
(1)自、互导纳的物理意义
自导纳在数值上等于与该节点I直接连接的所有支路导纳的总和。
如。
互导纳在数值上等于连接节点、支路导纳的负值,即。
如。
(2)节点导纳矩阵YB为对称方阵。
(3)节点导纳矩阵YB为稀疏矩阵。
(4)节点导纳矩阵具有对角优势。
3.3.2节点导纳矩阵的形成
用直接形成法形成节点导纳矩阵YB。
节点导纳矩阵即可根据自导纳和互导纳的定义直接形成,也可用支路——节点关联矩阵计算。
3.3.3节点导纳矩阵的修改
(1)从原有网络引出一支路,同时增加一节点,节点导纳矩阵将增加一阶。
新增的对角元,;
新增的非对角元,;
原有矩阵中的对角元将增加,。
(2)在原有网络的节点、之间增加一支路。
,
(3)在原有网络的节点,之间切除一支路
,,
(4)原有网络的节点、之间的导纳由改变为:
,,
(5)原有网络节点i、j之间变压器的变比由改变为
;;
3.4高斯-赛德尔法
3.4.1高斯-赛德尔迭代法的基本原理
为了方便理解这个n维方程组的叠代求解方法,先从一元非线性方程的求解开始。
假设有一维方程,高斯法的基本原理是,先将方程转化为:
那么给定一个初值,代入就可以得到一个新值,第k次叠代的值为:
一直叠代到误差满足要求为止,即
其中为事先设定的允许误差。
其计算流程如图3.1所示。
图3.1高斯迭代法的计算流程
这个解方程的方法称为高斯叠代法。
这个叠代求解的过程可以这样来理解:
的解可以认为是两个曲线和的交点的横坐标,首先给定一个初值,与斜线的交点的横坐标即为叠代后的新解,与斜线的交点的横坐标即为叠代后的新解,如此围绕交点往复循环,不断地逼近方程的解,如图所示。
图3.2高斯迭代法的几何解释
高斯迭代法可以推广到n维非线性代数方程组,假设n为方程组为:
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 matlab 电力系统 潮流 计算 毕业论文
文档标签
- 基于matlab电力系统短路
- 基于Matlab电力系统故障
- 基于Matlab低压电力系统
- matlab电力系统潮流计算
- 基于MATLAB电力系统保护
- MATLAB毕业论文
- 基于PSASP电网潮流
- 基于MATLAB的潮流计算论文精品基于MATLAB
- 毕业论文基于matlab转换
- 基于MATLAB电力系统稳定
- Matlab计算潮流电力系统
- 基于Matlab电力系统无功
- 毕业设计matlab电力系统分析
- 基于浪潮K1系统
- Matlab毕业论文
- 基于MATLAB电力系统潮流
- 电力系统潮流计算MATLAB
- MATLAB潮流计算论文
- 电力系统课程短路电流
- 基于MATLAB整流电路
- 电力系统三相短路电流
- 基于matlab风力发电
- 基于MATLAB短路电流
- 电力系统对称短路计算