节点配电网的建模与仿真

《节点配电网的建模与仿真》由会员分享，可在线阅读，更多相关《节点配电网的建模与仿真（24页珍藏版）》请在装配图网上搜索。

1、大学生课外创新实验竞赛总结报告项目简介：在配电网的正常运行中， 随着用电负荷的变化和系统运行方式的改变， 网络中的损耗也将发生变化。 要严格保证所有的用户在任何时刻都有额定的电压是不可能的， 因此配电网的某个节点电压偏移额定值是无法避免的。 为保证电力系统的正常运行， 需要对各节点处的电气量进行监测。 在电力系统发生故障后， 各节点的电气量也会发生较大波动， 也需要对各节点处的电气量进行监测。 对配电网进行建模仿真， 可以了解配电网各个节点的电流， 电压特征，仿真结果对配电网的实际运行有着正要的实际意义。在PSCAD/EMTDC台上搭建13节点中性点不接地配电网的仿真模型，仿真了不同故障情况对

2、故障点的电压， 电流的幅值、 相位特征的影响， 并得到各节点处的节点电压和电流的波形。 将各节点电压和电流的数据导入 matlab 中对数据进行绘图。 并对数据进行频谱分析，在频域上观察电流电压数据的基本特征。为了提取各节点处的以及故障处电流电压的特征量，基于PSCAD/EMTDC仿真环境，搭建了 13 节点中性点不接地配电网的仿真模型。线路模型采用分布式参数模型进行仿真。对于故障分析， 改变故障的发生位置， 故障的发生时间以及故障的类型进行仿真， 观察各节点处电流、 电压的仿真曲线， 并对结果进行简单分析。 深入研究故障对整个配电网各节点处的电流电压特征量的影响，为故障的确定以及排除提供了理

3、论依据。关键词： PSCAD/EMTDC matlab 潮流分析 频谱分析 Ag 故障 配电网项目创新点 ：在PSCAD/EMTDC台上搭建13节点中性点不接地配电网的仿真模型，其中输电线路采取分布式参数模型进行建模。 改变故障的发生时间， 发生位置以及故障类型， 在 PSCAD/EMTDC平台上进行仿真。将在 PSCAD/EMTDC得到的各节点电压和电流的数据导入matlab中对数据进行绘图， 并对不同节点的电流电压波形进行对比， 简要分析造成差异的主要原因。 同时对各节点所采集的电流电压数据进行频谱分析。第一章绪论 1项目背景 1研究意义 1第二章PSCAD勺软件介绍与潮流计算简介2PSC

4、A瞅件的介绍 2潮流计算简介 3第三章13节点配电网的建模4仿真模块的设计 4配电网的建模 9第四章 仿真结果的分析 10不同故障下的仿真波形 10matlab 绘图 16频谱分析 18附件 20第一章 绪论项目背景配电网是由架空线路、电缆、杆塔、配电变压器、隔离开关、无功补偿电容以及一些附属设施等组成的，在电力网中起重要分配电能作用的网络。配电网按电压等级来分类可分为高压配电网(35110KV),中压配电网(610KV,苏州有 20KV 的) ，低压配电网( 220/380V ) ；在负载率较大的特大型城市， 220KV 电网也有配电功能。 按供电区的功能来分类可分为城市配电网， 农村配电网

5、和工厂配电网等。 在城市电网系统中，主网是指110KV及其以上电压等级的电网，主要起连接区域高压(220KV及以上)电网的作用配电网是指35KV 及其以下电压等级的电网，作用是给城市里各个配电站和各类用电负荷供给电源。配电网络的拓扑分析是根据配电电气元件的连接关系， 把整个配电网络看成线与点结合的拓扑图， 然后根据电源结点、 开关结点等进行整个网络的拓扑连线分析， 它是配电网络进行状态估计、潮流计算、故障定位、隔离及供电恢复、网络重构等其它分析的基础。配电网络的结构庞大且复杂， 网络结构由于故障或负荷转移操作中开关的开合， 经常发生变化。 作为配电网络分析的基础， 网络拓扑计算需要进一步提高，

6、 因此迫切需要一个好的网络拓扑算法。 好的网络拓扑算法应该有效且直观， 它不仅能满足配电网自动化中的不同高级功能的要求， 还应能实现配电网络连通性的快速跟踪和识别， 适应事件变化。 同时还应节省存储空间和其他高级计算功能的时间。 目前国内外在这方面现有的研究有关联表矩阵表示法、网基矩阵表示法、结点消去法、树搜索表示法、离散处理法等。研究意义电力系统的出现，使高效，无污染，使用方便，易于控制的电能得到广泛应用，推动了社会生产各个领域的发展， 开创了电力时代， 发生了第二次技术革命。 潮流计算是电力系统最基本最常用的计算。 根据系统给定的运行条件， 网络接线及元件参数， 通过潮流计算可以确定各母线

7、的电压幅值和相角， 各元件流过的功率， 整个系统的功率损耗。 潮流计算是实现电力系统安全经济发供电的必要手段和重要工作环节O潮流计算在电力系统的规划计算，生产运行， 调度管理及科学计算中都有着广泛的应用。对于正在规划的电力系统， 通过潮流计算， 可以为选择电网供电方案和电气设备提供依据。 潮流计算还可以为继电保护和自动装置定整计算、电力系统故障计算和稳定计算等提供原始数据。第二章PSCA欧件的介绍与潮流计算简介PSCAD 软件的介绍PSCAD 是一种有效的用户图形界面, 能够显著地提高电力系统电磁瞬时模拟研究的效率。利用PSCADC族的软件工具，使得电力系统工程师能够充分利用现代微机工作站的资

8、源 , 更为有效地使用马尼托巴高压直流研究中心的EMTD瞬时模拟软件。该族软件还可作为该中心的实时数字模拟器（RTDS的用户界面。PSCAD 由下述软件模块构成：档管理系统、建模（DRAFTS莫块、架空线（T-LINE）和电缆 （CABLE）模块、运行（RUNTIME）模块、单曲线绘图（UNIPLOT）和多曲线绘图（MULTIPLOT膜块。 PSCAD/EMTD面时间域描述和求解完整的电力系统及其控制的微分方程（包括电磁和机电两个系统 ） 。这一类的模拟工具不同于潮流和暂态视定的模拟工具。后者是用稳态解去描述电路（即电磁过程）。但是在解电机的机械动态 （即转动惯量）微分方程。PSCAD/EMT

9、DC结果是 作为时间的即时值被求解。但通过内置的转换器和测量功能（ 象实有效值表计或者快速育里叶变换频谱分析等） 。 这些结果能被转换为矢量的幅值和相角。 实际系统的测量能够通过很多途径来完成。由于潮流和稳定的程序是通过稳定方程来代表它们只能基频段幅值和相位。因此PSCAD的模拟结果能够产生电力系统所有频率的相应限制仅在于用户自己选择的时间 步长。这种时间步长可以在毫秒到秒之间变化。基于以上简单说明我们可以了解到PSCAD1电力专业十分有用的仿真软件。我们组的课题系统节点较多， 传统的手工计算显然不切实际。 于是要寻找一种简单科学的计算方法来替代传统的手工计算，从而提高计算效率。因此利用PSC

10、A昉真运行出结果就成了本次课程设 计最为关键的一个环节。利用PSCAD寸13节点的配电网进行的仿真可以快速准确得出各母线上的潮流分布及系统的其他运行状态输出显示快速、明了。1） PSCAD/EMTDCF 发背景最初EMTD玳表直流暂态。最初的EMTD跌码由丹尼斯.伍德福德在1975年开始编写。编写这个程序的原因是因为当时现存的研究工具不能够曼尼托巴电力局对尼尔逊河高压直流工程进行强有力和灵活的研究的要求。 自此之后程序被不断开发， 至今已被广泛地应用在电力系统许多类型的模拟研究，其中包括交流研究，雷电过电压和电力电子学研究。 EMTDC开始时在大型计算机上使用。然后在1986年被移植到Unix

11、系统和以后的PC机上。PSCAD代表电力系统计算机辅助设计，并且是EMTDC勺图形用户接口。PSCAD V1198孙首先在阿波罗工作站上使用，然后大约在1995 年 PSCADV2 开始应用。 （ 至今 PSCADV2 仍然广泛地应用于 Unix系统上。）PSCAD V3以PC Windows作为平台，在 1999年面世。PSCAD V3 的 Unix 版本尚未完成。2） PSCAD/EMTDC 应用范围PSCAD/EMTDC 在时间域描述和求解完整的电力系统及其控制的微分方程（ 包括电磁和机电两个系统） 。这一类的模拟工具不同于潮流和暂态视定的模拟工具。后者是用稳态解去描述电路（即电磁过程）

12、。但是在解电机的机械动态 （即转动惯量）微分方程。PSCAD/EMTDC结果是作为时间的即时值被求解。但通过内置的转换器和测量功能（象实有效值表计，或者快速育里叶变换频谱分析等） 。这些结果能被转换为矢量的幅值和相角。实际系统的测量能够通过很多途径来完成。 由于潮流和稳定的程序是通过稳定方程来代表，它们只能基频段幅值和相位。因此PSCAD勺模拟结果能够产生电力系统所有频率的相应，限制仅在于用户自己选择的时间步长。 这种时间步长可以在毫秒到秒之间变化。 典型的研究包括：研究电力系统中由于故障或开关操作引起的过电压。它也能模拟变压器的非线性 （ 即饱和） 这一决定性因素。多运行工具（Mnltipl

13、e mnfacilities） 经常用来进行数以百计的模拟从而在下 列不同情况下发生故障时最坏的情况。故障发生在波形的不同位置，故障的类型不同，故障点不同。潮流计算简介潮流计算是电力系统分析中的一种最基本的计算，指在给定电力系统网络拓扑、 元件参 数和发电、 负荷参量条件下，计算有功功率、 无功功率及电压在电力网中的分布。 潮流计算是根据给定的电网结构、 参数和发电机、 负荷等元件的运行条件， 确定电力系统各部分稳态运行状态参数的计算。通常给定的运行条件有系统中各电源和负荷点的功率、枢纽点电压、衡点的电压和相位角。 待求的运行状态参量包括电网各母线节点的电压幅值和相角， 以及各支路的功率分布、

14、 网络的功率损耗等。 在电力系统运行和规划中， 都需要研究电力系统稳定运行情况， 确定电力系统的稳态运行状态。 给定电力系统的网络结构、 参数和决定电力系统运行状况的边界条件，确定电力系统运行的方法之一是朝流计算。从数学上说： 朝流计算是要求解一组有潮流方程描述的非线性方程组。 电力系统潮流计算是电力系统分析中最重要最基本的计算， 是电力运行、 规划以及安全性、 可靠性分析和优化的基础， 也是各种电磁暂态和机电暂态分析的基础和出发点。 电力系统潮流计算是电力系统最基本的计算， 也是最重要的计算。 所谓潮流计算， 就是已知电网的接线方式与参数及运行条件， 计算电力系统稳态运行各母线电压、 个支路

15、电流与功率及网损。 对于正在运行的电力系统，通过潮流计算可以判断电网母线电压、支路电流和功率是否越限， 如果有越限，就应采取措施，调整运行方式。对于正在规划的电力系统， 通过潮流计算，可以为选择电网供电方案和电气设备提供依据。 潮流计算还可以为继电保护和自动装置定整计算、 电力系统故障计算和稳定计算等提供原始数据。 此外， 在进行电力系统稳态及暂态稳定计算时， 要利用潮流计算的结果作为其计算的基础； 一些故障分析以及优化计算也需要有相应的潮流计算作配合； 潮流计算往往成为上述计算程序的一个重要组成部分。 以上这些， 主要是在系统规划设计及运行方式安排中的应用， 属于离线计算范畴。 随着现代化的

16、调度中心的建立， 为了对电力系统进行实时安全监控，需要根据实时数据库所提供的信息。第三章 13 节点配电网的建模仿真模块的设计1）设置仿真时间和步长新建的仿真工程，先应对“工程”的仿真时间、步长进行设置（也可在建好模型仿真开始前完成)。在“工程”模型窗口空白处鼠标右击，选择 Project Setting ,出现设置窗口，如图3-1所示，在这里可对本“工程”的仿真时间、计算步长、PSCA侬图步长等进行设定。仿真时间Duration of run ” 设为 1s,计算步长EMTDC time step ( us )”设为 5,绘图步长 “ PSCAD plot step ( us )” 设为 5

17、0。图1仿真时间设置2)电源的设置：电源模型如下图所示:R-0-0-图2电源模型在下面的界面中设置电源的名字，电源内部阻抗，以及电源是是否接地。source_3 Tfiree Ph3s-e Voltage- Source M_图3电源定义在下面的界面中设置电源的基本参数，具体设置如下图所示3）主变压器的设定主变压器模型如下图所示:scxjy 刃 Three Phase Voltage，凸urre M.图4电源参数设置图5变压器模型具体设置如下图所示:图6变压器设置4）开关/断路器的设置开关/断路器的模型如下图所示:图7断路器模型开关（断路器）从 master库到Breakers单元中选取，库中

18、有：单相、三相开关,相开关有单线和三线连接方式，开关必须与一个Logic模块配合使用，即：需将 “BRK和“BRK- Time Breaker Logic 都复制到工程中去。具体的设置如下图所示:tbrsakn TimFdi Breaker Logic图8断路器逻辑设置5）架空线（分布参数线路）模型 架空线的模型如下图所示：Tlme4Tline4图9架空线模型具体的参数设置如下：根据每一段线路的具体长度进行设置：下面以 Tline3段的线路参数设置为例:Tri-isni&sian LieCan HelpWK T阻把5KT*aegeRr汕 1t二- 5k1 . :图10架空线参数设置5）故障的设

19、置：故障模块如下图所示:图11故障模型故障发生在线路中间，因此将线路分为两段进行建模。故障的设置以及逻辑控制如下图所示:图12故障设置6）负荷的设置eiu 6-11eub 川qF平也啤创一宦图13负载种类对于无变压器容量的节点，我们以一个固定负载进行仿真分析；对于标有变压器容量的 节点，我们以仪的变压器 +固定负载进行仿真分析，变压器的一次侧和二次侧电压根据前后 负载节点的额定电压进行确定。配电网的建模综上所述，在 PSCAD/EMTDC对如下图所示的13节点的配电网进行建模。拓扑结构如下图所示:220kV650010kV646-1.5km架空线6452km架空线6323km架空线1.2km

20、架空线6330.8km架空线634T6116846710.5km架空线1km架空线B 16926754km架空线2.5km电缆线0.6km架空线3.5km架空线65268014拓扑结构在PSCA砰建立的仿真模型如下图所示:图15整体仿真模型第四章仿真结果的分析不同故障下的仿真波形1）故障类型为ABC三相接地，故障发生时间为，持续时间为。在ABC三相接地时，配电 网仍然是对称的，故我们只观察单相电压电流的仿真波形:图16节点632的单相电流电压图17节点633的单相电流电压*!事! 罩图18节点634的单相电流电压图19节点645的单相电流电压!一-图20节点646的单相电流电压图21节点671

21、的单相电流电压图22节点692的单相电流电压图23节点675的单相电流电压图24节点684的单相电流电压图25节点611的单相电流电压图26节点652的单相电流电压2）改变故障类型发生不对称故障，例如 AB相发生接地故障。此时节点的三相电流电压将发生变化，在此我们只给出632节点，684节点，692节点的仿真波形。仿真结果如下图所示：图27节点632三相的电流与电压图28节点684的三相电流与电压fifgI jg |祠， * * *餐qR图29节点692的三相电流与电压3）改变故障发生的时间故障类型ABC三相接地故障。故障时间改为 ,持续时间。选取632节点，684节点，692节点的仿真波形：

22、图30节点632的单相电流与电压图31节点684的单相电流与电压图32节点692的单相电流与电压4）改变故障发生的地点:线路684-671在、距671处发生Ag故障。图33节点点632的单相电流与电压图34节点684的单相电流与电压图35节点692的单相电流与电压matlab绘图在将所有的仿真结果导入到matlab中，对各个电气量在 matlab中进行绘图:以ABC三相接地故障为例，故障发生时间为，故障发生位置为线路 684-671在、距671处。将电流电压的数据导入matlab,将生成一个20001*77的矩阵，第一列代表时间，后面每列数据代表一个电气量。我们将第一个电气量在matlab绘图

23、如下所示:信号p(t)pp图37节点632的单相电流图36节点632的单相电压信号p(t)所有的电气量在 matlab中的绘图如下所示:加0 0101 20 赫。邮 10川科唧u.5 10. 0.05，也刑NM 0Mi:t).510 0.5 10 0.5 10 01 20 200.5 100.5 100.5 100.05 0.05 1Q-0.5 1恻。.：电:1喷2。加n rfc.1.10 010).5 1 h 0.5 1-。0.5 1-。0.5 10 050 05 1200.5 18 0.5 1 30 020 0 05_).5 1 1。 0.5 10 05比0 M腓。：血也1d fl n n

24、kIr恻.0 0.5 10 0.5 10 0.5 10 0.5 10 0.5 10 0.5 10 0.5 10 0.5 10.0510 20-11。-10.05-0.05Z111-0.010抑耐。IB0|【I帅地则0.5土BF! M-0.0i 0n nrz-101 0 0111 0 0随1 0 01,- -20.5 10,10-J *i rLr n-0.05-10 01 n n|d_00.50.5T -101 01。0.5- -0.0510-I 0 50.5 0.5 10n j0.5t -0.0110200.5-0.0110 0.51。之01m11414.5 1 - 0, 0 1_1州(I10

25、 吊1 -0 0.5 1 - .0 mif 1 010_脚/-0.000。1-0.1 000.50.5-101 00 10.5-0.5-100 05-0.50.50.5。10 0.51 0.1U .5 10 0 1_叫呼力赚愀】。 c nt-ih)10 010州也 .5 1 0 . .5 1 0IL1 0 01_mi0MHi . 0.50.5-0.11 0200.5-0.051020_0.50.50.5-0.110 0.5r 0 05r.收/，MC0.5 J;1用计.修W州hi 1 0 |晒1 0 01史。的0.50-0.1100.5 100.5-201 00.50.50.5-0.1-100.

26、5-0.0510 0.5111X牌喇 00.51：0嬲即由心 0 0.5 10 0.5 11 0.021 0.021 0.0210.50Ml小。啊w册时赭010.021 -0.021 -0.0210 0.5 10 0.5 10 0.5 10 0.5 110.50.02 0.02 0.02 0.50删阳,0卿如。.删制。脏融-0.5-0.020.02-1-0.02 -0.500.510 0.5 10 0.5 10 0.5 10 0.5 1图38单相电流或者电压其中：每一个代表某一节点的单相电流或者电压。频谱分析利用傅里叶变换的方法对振动的信号进行分解，并按频率顺序展开， 使其成为频率的函数，进而

27、在频率域中对信号进行研究和处理的一种过程，称为频谱分析。由时间函数求频谱函数的傅里叶变换公式就是将该时间函数乘以以频率为系数的指数函数之后，在从负无限大到正无限大的整个区间内，对时间进行积分， 这样就得到了与这个时间函数对应的，以频率为自变量的频谱函数。 频谱函数是信号的频域表示方式。根据上述傅里叶变换公式， 可以求出常数（直流信号）的频谱函数为频域中位于零频率处的一个冲激函数，表示直流信号就是一个频率等于零的信号。与此相反，冲激函数的频谱函数等于常数，表示冲激函数含有无限多个、频率无限密集的正弦成分。同样的，单个正弦波的频谱函数就是频域中位于该正弦波 频率处的一对冲激函数。在matlab中编

28、程做出节点632单相电压波形的频谱图，如下图所示:0.10.20.30.40.50.60.70.80.91100150200250300350400450500550Frequency(Hz)OL.SrrrrrI100150200250300350400450500550Frequency(Hz)图39频谱分析图附件附件1 PSCAD仿真模型工三之二z二附件二matlab绘图程序a1=load();a2=load();a3=load();a4=load();a5=load();a6=load();a7=load();a=a1,a2,a3,a4,a5,a6,a7;琳个电气量绘图t=a(：,1)；

29、p = a(:,26)figure(1)plot(t,p)；grid ontitle( 信号 p(t);xlabel(t)ylabel(p)%所有电气量绘图t=a(:,1);for i=2:77;x=a(:,i);subplot(9,9,i)plot(t,x)end附件 3 matlab 频谱分析Fs=1000;t=a(:,1);y=a(:,2);subplot(3,1,1);plot(t,y)N =10000;Y = fft(y,N)/(N/2); % /(N/2)归一化到 1 以内k = (0:N/2-1)*Fs/N; % kmax = (N/2)*Fs/N = Fs/2(Nyquist), Fs/N 为频率分辨率subplot(3,1,2);plot(k+50),abs(Y(1:N/2);xlabel(Frequency(Hz);ylabel(|Y|);subplot(3,1,3);plot(k+50),angle(Y(1:N/2);xlabel(Frequency(Hz);ylabel(sita);