电力系统的MATLAB_SIMULINK仿真与应用6

《电力系统的MATLAB_SIMULINK仿真与应用6》由会员分享，可在线阅读，更多相关《电力系统的MATLAB_SIMULINK仿真与应用6（162页珍藏版）》请在金锄头文库上搜索。

1、第6章 电力系统稳态与暂态仿真 第6章 电力系统稳态与暂态仿真 6.1 Powergui模块模块 6.2 电力系统稳态仿真电力系统稳态仿真 6.3 电力系统电磁暂态仿真电力系统电磁暂态仿真 6.4 电力系统机电暂态仿真电力系统机电暂态仿真 习题习题 第6章 电力系统稳态与暂态仿真 6.1 Powergui模块模块Powergui模块为电力系统稳态与暂态仿真提供了有用的图形用户分析界面。通过Powergui模块，可以对系统进行可变步长连续系统仿真、定步长离散系统仿真和相量法仿真，并实现以下功能：(1) 显示测量电压、测量电流和所有状态变量的稳态值；(2) 改变仿真初始状态；(3) 进行潮流计算并

2、对包含三相电机的电路进行初始化设置；(4) 显示阻抗的依频特性图；第6章 电力系统稳态与暂态仿真 (5) 显示FFT分析结果；(6) 生成状态空间模型并打开“线性时不变系统”(LTI)时域和频域的视窗界面；(7) 生成报表，该报表中包含测量模块、电源、非线性模块和电路状态变量的稳态值，并以后缀名.rep保存；(8) 设计饱和变压器模块的磁滞特性。6.1.1 主窗口功能简介主窗口功能简介MATLAB提供的Powergui模块在SimPowerSystems库中，图标如图6-1所示。第6章 电力系统稳态与暂态仿真 图6-1 Powergui模块图标第6章 电力系统稳态与暂态仿真 双击Powergu

3、i模块图标将弹出该模块的主窗口，如图6-2所示。该主窗口包含“仿真类型”(Simulation Type)和“分析工具”(Analysis Tools)两块内容，简介如下。1. 仿真类型仿真类型1) “相量法仿真”(Phasor simulation)单选框点击该单选框后，在该单选框下方的“频率”(Frequency)文本框中输入指定的频率，进行相量法分析。若未选中该单选框，“频率”文本框显示为灰色。第6章 电力系统稳态与暂态仿真 图6-2 Powergui模块主窗口第6章 电力系统稳态与暂态仿真 2) “离散系统仿真”(Discretize electrical model)单选框点击该单选

4、框后，在“采样时间”(Sample time)文本框中输入指定的采样时间(Ts0)，按指定的步长对离散化系统进行分析。若采样时间等于0，表示不对数据进行离散化处理，采用连续算法分析系统。若未选中该单选框，“采样时间”文本框显示为灰色。3) “连续系统仿真”(Continuous)单选框点击该单选框后，采用连续算法分析系统。4) “显示分析信息”(Show message during analysis)复选框选中该复选框后，命令窗口中将显示系统仿真过程中的相关信息。第6章 电力系统稳态与暂态仿真 2. 分析工具分析工具1) “稳态电压电流分析”(Steady-State Voltages an

5、d Currents)按键打开稳态电压电流分析窗口，显示模型文件的稳态电压和电流。2) “初始状态设置”(Initial States Setting)按键打开初始状态设置窗口，显示初始状态，并允许对模型的初始电压和电流进行更改。3) “潮流计算和电机初始化”(Load Flow and Machine Initialization)按键打开潮流计算和电机初始化窗口。第6章 电力系统稳态与暂态仿真 4) “LTI视窗”(Use LTI Viewer)按键打开窗口，使用“控制系统工具箱”(Control System Toolbox)的LTI视窗。5) “阻抗依频特性测量”(Impedance

6、vs. Frequency Measurement)按键打开窗口，如果模型文件中含阻抗测量模块，该窗口中将显示阻抗依频特性图。6) “FFT分析”(FFT Analysis)按键打开FFT分析窗口。7) “报表生成”(Generate Report)按键打开窗口，产生稳态计算的报表。第6章 电力系统稳态与暂态仿真 8) “磁滞特性设计工具”(Hysteresis Design Tool)按键打开窗口，对饱和变压器模块和三相变压器模块的铁芯进行磁滞特性设计。9) “计算RLC线路参数”(Compute RLC Line Parameters)按键打开窗口，通过导线型号和杆塔结构计算架空输电线的R

7、LC参数。6.1.2 稳态电压电流分析窗口稳态电压电流分析窗口打开“稳态电压电流分析”窗口如图6-3所示。该窗口中含有以下内容：第6章 电力系统稳态与暂态仿真 图6-3 “稳态电压电流分析”窗口第6章 电力系统稳态与暂态仿真 (1) “稳态值”(Steady state value)列表框：显示模型文件中指定的电压、电流稳态值。(2) “单位”(Units)下拉框：选择将显示的电压、电流值是“峰值”(Peak)还是“有效值”(RMS)。(3) “频率”(Frequency)下拉框：选择将显示的电压、电流相量的频率。该下拉框中列出模型文件中电源的所有频率。(4) “状态”(States)复选框：

8、显示稳态下电容电压和电感电流的相量值。默认状态为不选。(5) “测量”(Measurements)复选框：显示稳态下测量模块测量到的电压、电流相量值。默认状态为选中。第6章 电力系统稳态与暂态仿真 (6) “电源”(Sources)复选框：显示稳态下电源的电压、电流相量值。默认状态为不选。(7) “非线性元件”(Nonlinear elements)复选框：显示稳态下非线性元件的电压、电流相量值。默认状态为不选。(8) “格式”(Format)下拉框：在下拉列表框中选择要观测的电压和电流的格式。“浮点格式”(floating point)以科学计数法显示5位有效数字；“最优格式”(best o

9、f)显示4位有效数字并且在数值大于9999时以科学计数法表示；最后一个格式直接显示数值大小，小数点后保留2位数字。默认格式为“浮点格式”。(9) “更新稳态值”(Update Steady State Values)按键：重新计算并显示稳态电压、电流值。第6章 电力系统稳态与暂态仿真 6.1.3 初始状态设置窗口初始状态设置窗口仿真时，常常希望仿真开始时系统处于稳态，或者仿真开始时系统处于某种初始状态，这时，就可以使用“初始状态设置”按键。打开“初始状态设置”窗口如图6-4所示。该窗口中含有以下内容：(1) “初始状态”(Initial state values for simulation)

10、列表框：显示模型文件中状态变量的名称和初始值。(2) “设置到指定状态”(Set selected state)文本框：对“初始状态”列表框中选中的状态变量进行初始值设置。(3) “设置所有状态量”(Reset all States)：选择从“稳态”(To Steady State)或者“零初始状态”(To Zero) 开始仿真。第6章 电力系统稳态与暂态仿真 (4) “加载状态”(Reload States)：选择从“指定的文件”(From File)中加载初始状态或直接以“当前值”(From Diagram)作为初始状态开始仿真。第6章 电力系统稳态与暂态仿真 图6-4 “初始状态设置”窗

11、口第6章 电力系统稳态与暂态仿真 (5) “应用”(Apply)按键：用设置好的参数进行仿真。(6) “返回”(Revert)按键：返回到“初始状态设置”窗口打开时的原始状态。(7) “保存初始状态”(Save Initial States)按键：将初始状态保存到指定的文件中。(8) “格式”(Format)下拉框：选择观测的电压和电流的格式。格式类型见6.1.2节。默认格式为“浮点格式”。(9) “分类”(Sort values by)下拉框：选择初始状态值的显示顺序。“默认顺序”(Default order)是按模块在电路中的顺序显示初始值；“状态序号”(State number)是按状态

12、空间模型中状态变量的序号来显示初始值；“类型”(Type)是按电容和电感来分类显示初始值。默认格式为“默认顺序”。第6章 电力系统稳态与暂态仿真 6.1.4 潮流计算和电机初始化窗口潮流计算和电机初始化窗口打开“潮流计算和电机初始化”窗口如图6-5所示。该窗口中含有以下内容：(1) “电机潮流分布”(Machines load flow)列表框：显示“电机”(Machines)列表框中选中电机的潮流分布。(2) “电机”(Machines)列表框：显示简化同步电机、同步电机、非同步电机和三相动态负荷模块的名称。选中该列表框中的电机或负荷后，才能进行参数设置。第6章 电力系统稳态与暂态仿真 图6

13、-5 潮流计算和电机初始化窗口第6章 电力系统稳态与暂态仿真 (3) “节点类型”(Bus type)下拉框：选择节点类型。对于“PV节点”(P&V Generator)，可以设置电机的端口电压和有功功率；对于“PQ节点”(P&Q Generator)，可以设置电机的有功和无功功率；对于“平衡节点”(Swing Bus)，可以设置终端电压UAN的有效值和相角，同时需要对有功功率进行预估。如果选择了非同步电机模块，则仅需要输入电机的机械功率；如果选择了三相动态负荷模块，则需要设置该负荷消耗的有功和无功功率。(4) “终端电压UAB”(Terminal voltage UAB)文本框：对选中电机的

14、输出线电压进行设置(单位：V)。第6章 电力系统稳态与暂态仿真 (5) “有功功率”(Active power)文本框：设置选中的电机或负荷的有功功率(单位：W)。(6) “预估的有功功率”(Active power guess)文本框：如果电机的节点类型为平衡节点的话，设置迭代起始时刻电机的有功功率。(7) “无功功率”(Reactive power)文本框：设置选中的电机或负荷的无功功率(单位：var)。(8) “电压UAN的相角”(Phase of UAN voltage)文本框：当电机的节点类型为平衡节点时，该文本框被激活。指定选中电机a相相电压的相角。第6章 电力系统稳态与暂态仿真

15、(9) “负荷频率”(Load flow frequency)下拉框：对潮流计算的频率进行设置，通常为60 Hz或者50 Hz。(10) “负荷潮流初始状态”(Load flow initial condition)下拉框：常常选择默认设置“自动”(Auto)，使得迭代前系统自动调节负荷潮流初始状态。如果选择“从前一个结果开始” (Start from previous solution)，则负荷潮流的初始值为上次仿真结果。如果改变电路参数、电机的功率分布和电压后负荷潮流不收敛，就可以选择这个选项。第6章 电力系统稳态与暂态仿真 (11) “更新电路和测量结果”(Update Circuit

16、& Measurements)按键：更新电机列表，更新电压相量和电流相量，更新“电机潮流分布”列表框中的功率分布。其中的电机电流是最近一次潮流计算的结果。该电流值储存在电机模块的“初始状态参数”(Initial conditions)文本框中。(12) “更新潮流分布”(Update Load Flow)按键：根据给定的参数进行潮流计算。第6章 电力系统稳态与暂态仿真 6.1.5 LTI视窗视窗打开“LTI视窗”窗口如图6-6所示。该窗口中含有以下内容：(1) “系统输入”(System inputs)列表框：列出电路状态空间模型中的输入变量，选择需要用到LTI视窗的输入变量。(2) “系统输

17、出”(System outputs)列表框：列出电路状态空间模型中的输出变量，选择需要用到LTI视窗的输出变量。(3) “打开新的LTI视窗”(Open New LTI Viewer)按键：产生状态空间模型并打开选中的输入和输出变量的LTI视窗。(4) “打开当前LTI视窗”(Open in current LTI Viewer)按键：产生状态空间模型并将选中的输入和输出变量叠加到当前LTI视窗。第6章 电力系统稳态与暂态仿真 图6-6 “LTI视窗”窗口第6章 电力系统稳态与暂态仿真 图6-7 “阻抗依频特性测量”窗口6.1.6 阻抗依频特性测量窗口阻抗依频特性测量窗口打开“阻抗依频特性测量

18、”窗口如图6-7所示。第6章 电力系统稳态与暂态仿真 该窗口中含有以下内容： (1) 图表：窗口左上侧的坐标系表示阻抗频率特性，左下侧的坐标系表示相角频率特性。(2) “测量模块”(Measurement)列表框：列出模型文件中的阻抗测量模块，选择需要显示依频特性的阻抗测量模块。使用“Ctrl”键可选择多个阻抗显示在同一个坐标中。(3) “范围”(Range)文本框：指定频率范围(单位：Hz)。该文本框中可以输入任意有效的MATLAB表达式。(4) “对数阻抗”(Logarithmic Impedance)单选框：坐标系纵坐标的阻抗以对数值形式表示。第6章 电力系统稳态与暂态仿真 (5) “线

19、性阻抗”(Linear Impedance)单选框：坐标系纵坐标的阻抗以线性形式表示。(6) “对数频率”(Logarithmic Frequency)单选框：坐标系横坐标的频率以对数值形式表示。(7) “线性频率”(Linear Frequency)单选框：坐标系横坐标的频率以线性形式表示。(8) “网格”(Grid)复选框：选中该复选框，阻抗频率特性图和相角频率特性图上将出现网格。默认设置为无网格。第6章 电力系统稳态与暂态仿真 (9) “更新后保存数据”(Save data when updated)复选框：选中该复选框后，该复选框下面的“工作间变量名” (Workspace varia

20、ble name)文本框被激活，数据以该文本框中显示的变量名被保存在工作间中。复数阻抗和对应的频率保存在一起。其中频率保存在第1列，阻抗保存在第2列。默认设置为不保存。(10) “显示/保存”(Display/Save)按键：开始阻抗依频特性测量并显示结果，如果选择了“更新后保存数据”(Save data when updated)复选框，数据将保存到指定位置。6.1.7 FFT分析窗口分析窗口打开“FFT分析”窗口如图6-8所示。第6章 电力系统稳态与暂态仿真 图6-8 “FFT分析”窗口第6章 电力系统稳态与暂态仿真 该窗口中含有以下内容：(1) 图表：窗口左上侧的图形表示被分析信号的波形

21、，窗口左下侧的图形表示该信号的FFT分析结果。(2) “结构”(Structure)下拉框：列出工作间中带时间的结构变量的名称。使用下拉菜单选择要分析的结构变量。这些结构变量名可以由“示波器”(Scope)模块产生。打开示波器模块参数对话框，选中“数据历史”(Data history)标签页，如图6-9所示，在“变量名”(Variable name)文本框中输入该结构变量的名称，在“存储格式”(Format)下拉框中选择“带时间的结构变量”(Structure with time)。第6章 电力系统稳态与暂态仿真 图6-9 由“示波器”模块产生结构变量第6章 电力系统稳态与暂态仿真 这些结构变

22、量名也可以由“到工作间”(To Workspace)模块产生，如图6-10所示。打开该模块的对话窗口，并在“变量名”(Variable name)文本框中输入该结构变量的名称，“存储格式”(Save format)为“结构变量”(Structure)。第6章 电力系统稳态与暂态仿真 图6-10 由“到工作间”模块产生结构变量(a) “到工作间”模块图标；(b) “到工作间”模块参数对话框第6章 电力系统稳态与暂态仿真 (3) “输入变量”(Input)下拉框：列出被选中的结构变量中包含的输入变量名称，选择需要分析的输入变量。(4) “信号路数”(Signal number)下拉框：列出被选中的

23、输入变量中包含的各路信号的名称。例如，若要把a、b、 c三相电压绘制在同一个坐标中，可以通过把这三个电压信号同时送入示波器的一个通道来实现，这个通道就对应一个输入变量，该变量含有3路信号，分别为a相、b相和c相电压。(5) “开始时间”(Start time)文本框：指定FFT分析的起始时间。(6) “周期个数”(Number of cycles)文本框：指定需要进行FFT分析的波形的周期数。第6章 电力系统稳态与暂态仿真 (7) “显示FFT窗/显示完整信号”(Display FFT window/Display entire signal)下拉框：选择“显示完整信号”(Display en

24、tire signal)，将在左上侧插图中显示完整的波形；选择“显示FFT窗”(Display FFT window)将在左上侧插图中显示指定时间段内的波形。(8) “基频”(Fundamental frequency)文本框：指定FFT分析的基频(单位：Hz)。(9) “最大频率”(Max Frequency)文本框：指定FFT分析的最大频率(单位：Hz)。(10) “频率轴”(Frequency axis)下拉框：在下拉框中选择“赫兹”(Hertz)使频谱的频率轴单位为Hz，选择“谐波次数” (Harmonic order)使频谱的频率轴单位为基频的整数次倍数。第6章 电力系统稳态与暂态仿

25、真 (11) “显示类型”(Display style)下拉框：频谱的显示类型可以是“以基频或直流分量为基准的柱状图”(Bar(relative to Fund. or DC)、“以基频或直流分量为基准的列表”(list (relative to Fund. or DC)、“指定基准值下的柱状图”(Bar (relative to specified base)、“指定基准值下的列表”(List (relative to specified base)四种类型。(12) “基准值”(Base value)文本框：当“显示类型”下拉框中选择“指定基准值下的柱状图”或“指定基准值下的列表”时，该文

26、本框被激活，输入谐波分析的基准值。(13) “显示”(Display)按键：显示FFT分析结果。第6章 电力系统稳态与暂态仿真 6.1.8 报表生成窗口报表生成窗口打开“报表生成”窗口如图6-11所示。该窗口中含有以下内容：(1) “报表中包含的内容”(Items to include in the report)：包括“稳态”(Steady state)复选框、“初始状态”(Initial states)复选框和“电机负荷潮流”(Machine load flow)复选框，这三个复选框可以任意组合。(2) “报表中的频率”(Frequency to include in the report

27、)下拉框：选择报表中包含的频率。可以是60 Hz或者全部，默认为60 Hz。第6章 电力系统稳态与暂态仿真 (3) “单位”(Units)下拉框：选择以“峰值”(Peak)还是“有效值”(Units)显示数据。(4) “格式”(Format)下拉框：与6.1.2节相关内容相同。(5) “报表生成”(Create Report)按键：生成报表并保存。第6章 电力系统稳态与暂态仿真 图6-11 “报表生成”窗口第6章 电力系统稳态与暂态仿真 6.1.9 磁滞特性设计工具窗口磁滞特性设计工具窗口打开“磁滞特性设计工具”窗口如图6-12所示。该窗口中含有以下内容：(1) “磁滞曲线”(Hysteres

28、is curve for file)图表：显示设计的磁滞曲线。(2) “分段”(Segments)下拉框：将磁滞曲线做分段线性化处理，并设置磁滞回路第1象限和第4象限内曲线的分段数目。左侧曲线和右侧曲线关于原点对称。(3) “剩余磁通”(Remanent flux Fr)文本框：设置零电流对应的剩磁。(4) “饱和磁通”(Saturation flux Fs)文本框：设置饱和磁通。第6章 电力系统稳态与暂态仿真 图6-12 “磁滞特性设计工具”窗口第6章 电力系统稳态与暂态仿真 (5) “饱和电流”(Saturation current Is)文本框：设置饱和磁通对应的电流。(6) “矫顽电流

29、”(Coercive current Ic)文本框：设置零磁通对应的电流。(7) “矫顽电流处的斜率”(dF/dI at coercive current)文本框：指定矫顽电流点的斜率。(8) “饱和区域电流”(Saturation region currents)文本框：设置磁饱和后磁化曲线上各点所对应的电流值，仅需设置第1象限值。注意该电流向量的长度必须和“饱和区域磁通”的向量长度相同。第6章 电力系统稳态与暂态仿真 (9) “饱和区域磁通”(Saturation region fluxes)文本框：设置磁饱和后磁化曲线上各点所对应的磁通值，仅需要设置第1象限值。注意该向量的长度必须和“饱

30、和区域电流”的向量长度相同。(10) “变压器额定参数”(Transfo Nominal Parameters)文本框：指定额定功率(单位：VA)、一次绕组的额定电压值(单位：V)和额定频率(单位：Hz)。(11) “参数单位”(Parameter units)下拉框：将磁滞特性曲线中电流和磁通的单位由国际单位制(SI)转换到标幺制(p.u.)或者由标幺制转换到国际单位制。第6章 电力系统稳态与暂态仿真 (12) “放大磁滞区域”(Zoom around hysteresis)复选框：选中该复选框，可以对磁滞曲线进行放大显示。默认设置为“可放大显示”。6.1.10 计算计算RLC线路参数窗口线

31、路参数窗口打开“计算RLC线路参数”窗口如图6-13所示。该窗口可分为三个子窗口，左上窗口输入常用参数(单位、频率、大地电阻和文件注释)，右上窗口输入线路的几何结构，下方窗口输入导线的特性。第6章 电力系统稳态与暂态仿真 图6-13 “计算RLC线路参数”窗口第6章 电力系统稳态与暂态仿真 1. 常用参数子窗口常用参数子窗口(1) “单位”(Units)下拉框：在下拉菜单中，选择以“米制”(metric)为单位时，以厘米作为导线直径、几何平均半径GMR和分裂导线直径的单位，以米作为导线间距离的单位；选择以“英制”(english)为单位时，以英寸作为导线直径、几何平均半径GMR和分裂导线直径的

32、单位，以英尺作为导线间距离的单位。(2) “频率”(Frequency)文本框：指定RLC参数所用的频率(单位：Hz)。(3) “大地电阻”(Ground resistivity)文本框：指定大地电阻(单位：m)。输入0表示大地为理想导体。第6章 电力系统稳态与暂态仿真 (4) “注释”(Comments)多行文本框：输入关于电压等级、导线类型和特性等的注释。该注释将与线路参数一同被保存。2. 线路几何结构子窗口线路几何结构子窗口(1) “导线相数”(Number of phase conductors(bundle)文本框：设置线路的相数。(2) “地线数目”(Number of groun

33、d wires(bundle)文本框：设置大地导线的数目。(3) 导线结构参数表：输入导线的“相序”(Phase number)、“水平挡距”(X)、“垂直挡距”(Y tower)、“挡距中央的高度”(Y min)、“导线的类型”(Conductor(bundle)type)共五个参数。第6章 电力系统稳态与暂态仿真 3. 导线特性子窗口导线特性子窗口(1) “导线类型的个数”(Number of conductor types or bundle types)文本框：设置需要用到导线类型(单导线或分裂导线)的数量。假如需要用到架空导线和接地导线，该文本框中就要填“2”。(2) “导线内电感计

34、算方法”(Conductor internal inductance evaluated from)下拉框：选择用“直径/厚度”(T/D ratio)、“几何平均半径”(Geometric Mean Radius(GMR)或者“1英尺(米)间距的电抗”(Reactance Xa at 1-foot spacing)进行内电感计算。(3) “考虑导线集肤效应”(Include conductor skin effect)复选框：选中该复选框后，在计算导线交流电阻和电感时将考虑集肤效应的影响。若未选中，电阻和电感均为常数。第6章 电力系统稳态与暂态仿真 (4) 导线特性参数表：输入导线“外径”(C

35、onductor Outside diameter)、“T/D”(Conductor T/D ratio)、“GMR” (Conductor GMR)、“直流电阻”(Conductor DC resistance)、“相对磁导率”(Conductor relative permeability)、“分裂导线中的子导线数目”(Number of conductors per bundle)、“分裂导线的直径”(Bundle diameter)、“分裂导线中1号子导线与水平面的夹角”(Angle of conductor 1)共八个参数。(5) “计算RLC参数”(compute RLC para

36、meters)按键：点击该按键后，将弹出RLC参数的计算结果窗口。第6章 电力系统稳态与暂态仿真 (6) “保存”(Save)按键：点击该按键后，线路参数以及相关的GUI信息将以后缀名.mat被保存。(7) “加载”(Load)按键：点击该按键后，将弹出窗口，选择“典型线路参数”(Typical line data)或“用户定义的线路参数” (User defined line data)将线路参数信息加载到当前窗口。第6章 电力系统稳态与暂态仿真 6.2 电力系统稳态仿真电力系统稳态仿真6.2.1 连续系统仿真连续系统仿真以例4.4为例，说明Powergui模块在电力系统稳态仿真分析中的应用

37、。【例【例6.1】计算例4.4的潮流分布，并利用Powergui模块实现连续系统的稳态分析。解：解：例4.4看似一个极简单的电路，可是它的潮流分析却不简单。已知电源侧的电压幅值、相角和额定电压下的负荷大小，这意味着在求出负荷侧电压前，负荷未知，因此按潮流计算的理论进行分析，这种电路必须经过一些假设，并反复迭代计算，直到逼近真实解。但是，利用Powergui模块却很简单。第6章 电力系统稳态与暂态仿真 图6-14 例6.1的系统仿真图(1) 重新布置系统仿真图，如图6-14所示。相对例4.4，本例添加了新的模块，新增模块的名称及提取路径如表6-2所示。第6章 电力系统稳态与暂态仿真 表表6-2

38、例例6.1新增模块的名称及提取路径新增模块的名称及提取路径(2) Powergui仿真。打开Powergui模块窗口，选中“连续系统仿真”(Continuous)单选框后，点击“稳态电压电流分析”按键，出现稳态电压电流分析窗口，如图6-15所示。第6章 电力系统稳态与暂态仿真 图6-15 例6.1的稳态电压电流分析窗口第6章 电力系统稳态与暂态仿真 图中，状态变量用电流和电压的符号加上电感或电容的模块名表示，例如“I1_110Mvar”表示110Mvar负荷上的电流大小，“Uc_input：Pi Line”表示PI形线路左侧并联电容器上的电压大小，“Uc_output：Pi Line”表示PI

39、形线路右侧并联电容器上的电压大小，“I1_section _1：Pi Line”表示第一段PI形线路串联电感上的电流大小。电压源上电压的名称与系统电压源名称一致，如“Vs”表示电压源Vs上的电压大小。测量模块测得的电压值用测量模块的名称表示，如“V1”表示电压表模块V1测得的电压大小(电源侧电压)，“V2”表示电压表模块V2测得的电压大小，“I2”表示电流表模块I2测得的电流大小。第6章 电力系统稳态与暂态仿真 由图6-15可见，PI形电路左侧的电压相量为244.88 0.19 kV，PI形电路右侧的电压相量为166.413.66kV，PI形电路上的电流为529.21.0886.12 A，负荷

40、侧电流为610.1786.15 A。因此，负荷大小为(6-1) 第6章 电力系统稳态与暂态仿真 图6-16所示为直接通过测量模块得到的PI形电路两侧电压和实际负荷大小。图中波形从上到下依次为PI形电路左侧电压、PI形电路右侧电压、负荷侧有功功率、负荷侧无功功率。该结果与Powergui所得结论一致。利用测量模块还可以得到电源侧电流的大小，这样就可以很容易求到线路各处的功率分布、功率损耗和电压损耗了。读者可以动手试一试，并和电力系统分析课程中潮流计算的相关内容进行比较。仿真时，特别是在对含电机的电路或者节点较多的电力系统，往往希望仿真开始时系统处于某种指定的初始状态。下例将说明对系统初始状态的设

41、置方法。第6章 电力系统稳态与暂态仿真 图6-16 例6.1的仿真波形图第6章 电力系统稳态与暂态仿真 【例【例6.2】改变例6.1的初始状态，使仿真开始时，电路处于零初始状态。解：解：打开Powergui模块窗口，点击“初始状态设置”按键，出现初始状态设置窗口，如图6-17所示。点击按键“to zero”将全部状态变量的初始值设为0，回到模型文件主窗口，重新开始仿真，得到仿真波形如图6-18所示。第6章 电力系统稳态与暂态仿真 图6-17 例6.2的初始状态设置窗口第6章 电力系统稳态与暂态仿真 图6-18 例6.2的仿真波形图第6章 电力系统稳态与暂态仿真 由仿真波形图可见，当仿真从零初始

42、状态开始时，线路两侧的电压在仿真初期不是理想正弦波。简单地说，就是合闸引起了一个暂态过程。从“波”的角度来讲，就是合闸后在线路上引起了“波”过程，“波”在阻抗不连续点(线路两端)不断进行折射和反射，形成了图6-18所示的波形。如果系统中存在电机模块，电机模块的初始状态需要通过“潮流计算和电机初始化”按键进行设置。第6章 电力系统稳态与暂态仿真 【例【例6.3】利用Powergui模块设置例4.2同步发电机的初始状态。解：解：例4.2中直接在同步发电机参数对话框中输入了初始状态，这些初始状态实际上是由Powergui模块自动产生的。首先在同步发电机参数对话窗口(见图6-19)的初始状态文本框中输

43、入初始状态0 0 0 0 0 0 0 0 1，确认后关闭该对话框。第6章 电力系统稳态与暂态仿真 图6-19 例6.3的同步发电机参数对话窗口第6章 电力系统稳态与暂态仿真 点击“潮流计算和电机初始化”按键，进入潮流计算和电机初始化设置窗口，如图6-20所示。该窗口中只显示了一个电机模块，名为SM_p.u.，选中该模块，并点击“更新潮流分布”按键，更新负荷潮流，窗口左侧将显示更新后的负荷潮流分布情况。现在重新打开同步发电机参数对话窗口(见图6-21)，可以看到，初始状态已经自动更新了。同步发电机的初始状态设置好后，剩下的仿真工作就可交给读者了。第6章 电力系统稳态与暂态仿真 图6-20 例6.

44、3的潮流计算和电机初始化设置窗口第6章 电力系统稳态与暂态仿真 图6-21 更新后的同步发电机参数对话窗口第6章 电力系统稳态与暂态仿真 6.2.2 离散系统仿真离散系统仿真连续系统仿真通常采用变步长积分算法。对小系统而言，变步长算法通常比定步长算法快，但是对含大量状态变量或非线性模块(如电力电子开关)的系统而言，采用定步长离散算法的优越性更为明显。对系统进行离散化时，仿真的步长决定了仿真的精确度。步长太大可能导致仿真精度不足，步长太小又可能大大增加仿真运行时间。判断步长是否合适的唯一方法就是用不同的步长试探并找到最大时间步长。对于50 Hz或60 Hz的系统，或者带有整流电力电子设备的系统，

45、通常2050 s的时间步长都能得到较好的仿真结果。第6章 电力系统稳态与暂态仿真 对于含强迫换流电力电子开关器件的系统，由于这些器件通常都运行在高频下，因此需要适当地减小时间步长。例如，对运行在8 kHz左右的脉宽调制(PWM)逆变器的仿真，需要的时间步长为1 s。【例【例6.4】将例4.4中的PI形电路的段数改为10，对系统进行离散化仿真并比较离散系统和连续系统的仿真结果。解：解： (1) 重新布置系统仿真图，如图6-22所示。第6章 电力系统稳态与暂态仿真 图6-22 例6-4的系统仿真图第6章 电力系统稳态与暂态仿真 (2) 参数设置。双击例4-4模型文件中PI形电路模块，打开参数对话框

46、，将分段数改为10，如图6-23所示。打开Powergui模块，选择“离散系统仿真”单选框，设置采样时间为25e-6 s，如图6-24所示。仿真时该系统将以25 s的采样率进行离散化。第6章 电力系统稳态与暂态仿真 图6-23 例6-4的PI形电路参数对话窗 第6章 电力系统稳态与暂态仿真 图6-24 例6-4的Powergui模块参数对话窗 第6章 电力系统稳态与暂态仿真 图6-25 例6-4的仿真参数设置对话窗由于系统离散化了，因此在该系统中无连续的状态变量，所以不需要采用变步长的积分算法进行仿真。打开菜单SimulationConfiguration parameters对话框，按图6-

47、25设置仿真参数，选择“定步长”(Fixed-step)和“离散”(discrete(no continuous states)选项并设置步长为25 s。第6章 电力系统稳态与暂态仿真 图6-26 例6-4的仿真运行时间(3) 仿真运行时间比较。为了得到仿真运行时间，在MATLAB命令窗口输入如下命令：tic; sim(gcs); toc仿真结束后，仿真所用的时间将以秒为单位显示在MATLAB命令窗口中，如图6-26所示。第6章 电力系统稳态与暂态仿真 可见，离散化系统后，仿真运行时间为0.188 s。将离散系统的采样时间设为0并回到连续系统的仿真状态，仿真算法改为连续积分算法ode23tb，

48、可以得到连续系统仿真需要的运行时间为0.219 s。因此，离散积分算法比连续积分算法更快。(4) 仿真精度比较。为了比较两种方法的精确度，执行以下三种仿真： 连续系统仿真，Ts=0 s； 离散系统仿真，Ts=25 s； 离散系统仿真，Ts=50 s。第6章 电力系统稳态与暂态仿真 如图6-27所示，双击并打开V2示波器模块，选择“参数”(Parameters)项，在打开的窗口中选择“数据历史”(Data history)，去掉“仅保留最新的数据点”(Limit data points to last)复选框，这样可以观察到整个仿真过程中的波形变化。选中“将数据保存到工作空间”(Save dat

49、a points to workspace)复选框，将变量名指定为V2，格式为“列”(Array)。第6章 电力系统稳态与暂态仿真 图6-27 例6-4示波器V2的参数设置(a) 波形；(b) 参数标签页 第6章 电力系统稳态与暂态仿真 开始连续系统仿真，仿真结束时间选为0.02 s。仿真结束后，在MATLAB命令窗口中输入命令：V2C=V2;这样，电压V2被保存在变量V2C中。重新开始仿真，将系统离散化，设置仿真步长Ts=25 s，注意仿真参数中的步长设置也要改为25 s，仿真结束时间为0.02 s。仿真结束后，将电压V2保存在变量V2d25中。再次仿真，设置仿真步长为Ts=50 s。仿真结

50、束后，将电压V2保存在变量V2d50中。第6章 电力系统稳态与暂态仿真 在MATLAB命令窗口中输入如下语句，可画出三种情况下的电压波形，如图6-28所示。plot(V2C(:,1),V2C(:,2),V2d25(:,1),V2d25(:,2), V2d50(:,1),V2d50(:,2)第6章 电力系统稳态与暂态仿真 图6-28 三种仿真方法波形比较第6章 电力系统稳态与暂态仿真 使用图形窗口中的放大功能，将目标集中到0.0045 s附近观察三种仿真的差别。如图6-29所示，25 s下的仿真结果与50 s的仿真结果一致，连续系统的仿真结果除了步长不同，结果也相同。可见，本例中，选择50 s的

51、步长不但可以提高计算速度而且不影响仿真的精确度。第6章 电力系统稳态与暂态仿真 图6-29 放大后三种仿真波形比较第6章 电力系统稳态与暂态仿真 6.2.3 相量法仿真相量法仿真相量是代表特定频率下的正弦电压和电流的复数，可以用直角坐标或者极坐标表示。相量法是电力系统正弦稳态分析的主要手段。它只关心系统中电压电流的相角和幅值，不需要求解电力系统状态方程，不需要特殊的算法，因此计算速度快得多。必须清楚的是，相量法给出的解是在特定频率下的解。【例【例6.5】用相量法分析例6.4。解：解：(1) 参数设置。打开Powergui模块，选择“相量法分析”单选框，并在“频率”对话框中将频率改为50 Hz。

52、关闭Powergui模块，模型文件主窗口中的Powergui模块图标显示为“相量法”(Phasors)分析，如图6-30所示。第6章 电力系统稳态与暂态仿真 图6-30 例6.5的Powergui模块相量法分析图标第6章 电力系统稳态与暂态仿真 图6-31 例6.5的电压测量模块V1打开电压测量模块V1，选择“幅值相角”(Magnitude-Angle)模式，如图6-31所示。电压测量模块V2也选择幅值相角模式。第6章 电力系统稳态与暂态仿真 (2) 仿真。开始仿真，得到输电线路送端V1和受端V2的电压幅值和相角，如图6-32所示。可见，V1侧电压幅值为1 p.u.，相角为0.19；V2侧电压

53、幅值为0.67 p.u.，相角为3.66。这和图6-16稳态分析的结论一致。第6章 电力系统稳态与暂态仿真 图6-32 例6.5的仿真结果(a) V1；(b) V2第6章 电力系统稳态与暂态仿真 6.3 电力系统电磁暂态仿真电力系统电磁暂态仿真SIMULINK的电力系统暂态仿真过程通过机械开关设备，如“断路器”(circuit breakers)模块或者电力电子设备的开断实现。6.3.1 断路器模块断路器模块SimPowerSystems库提供的断路器模块可以对开关的投切进行仿真。断路器合闸后等效于电阻值为Ron的电阻元件。Ron是很小的值，相对外电路可以忽略。断路器断开时等效于无穷大电阻，熄

54、弧过程通过电流过零时断开断路器完成。开关的投切操作可以受外部或内部信号的控制。第6章 电力系统稳态与暂态仿真 外部控制方式时，断路器模块上出现一个输入端口，输入的控制信号必须为0或者1，其中0表示切断，1表示投合；内部控制方式时，切断时间由模块对话框中的参数指定。如果断路器初始设置为1(投合)，SimPowerSystems库自动将线性电路中的所有状态变量和断路器模块的电流进行初始化设置，这样仿真开始时电路处于稳定状态。断路器模块包含Rs-Cs缓冲电路。如果断路器模块和纯电感电路、电流源和空载电路串联，则必须使用缓冲电路。第6章 电力系统稳态与暂态仿真 图6-33 单相断路器模块图标带有断路器

55、模块的系统进行仿真时需要采用刚性积分算法，如ode23tb、odel5s，这样可以加快仿真速度。1. 单相断路器模块单相断路器模块外部控制方式、带缓冲电路和不带缓冲电路的单相断路器模块图标如图6-33所示。第6章 电力系统稳态与暂态仿真 双击断路器模块，弹出该模块的参数对话框如图6-34。该对话框中含有如下参数：(1) “断路器电阻”(Breaker resistance Ron)文本框：断路器投合时的内部电阻(单位：)。断路器电阻不能为0。(2) “初始状态”(Initial state)文本框：断路器初始状态。断路器为合闸状态，输入1，对应的图标显示投合状态；输入0，表示断路器为断开状态。

56、(3) “缓冲电阻”(Snubber resistance Rs)文本框：并联缓冲电路中的电阻值(单位：)。缓冲电阻值设为inf时，将取消缓冲电阻。 第6章 电力系统稳态与暂态仿真 图6-34 单相断路器模块参数对话框第6章 电力系统稳态与暂态仿真 (4) “缓冲电容”(Snubber capacitance Cs)文本框：并联缓冲电路中的电容值(单位：F)。缓冲电容值设为0时，将取消缓冲电容；缓冲电容值设为inf时，缓冲电路为纯电阻性电路。(5) “开关动作时间”(Switching times)文本框：采用内部控制方式时，输入一个时间向量以控制开关动作时间。从开关初始状态开始，断路器在每个

57、时间点动作一次。例如，初始状态为0，在时间向量的第一个时间点，开关投合，第二个时间点，开关打开。如果选中外部控制方式，该文本框不可见。第6章 电力系统稳态与暂态仿真 (6) “外部控制”(External control of switching times)复选框：选中该复选框，断路器模块上将出现一个外部控制信号输入端。开关时间由外部逻辑信号(0或1)控制。(7) “测量参数”(Measurements)下拉框：对以下变量进行测量。 “无”(None)：不测量任何参数。 “断路器电压”(Branch voltages)：测量断路器电压。 “断路器电流”(Branch currents)：测量

58、断路器电流，如果断路器带有缓冲电路，测量的电流仅为流过断路器器件的电流。第6章 电力系统稳态与暂态仿真 图6-35 三相断路器模块图标 “所有变量”(Branch voltages and currents)：测量断路器电压和电流。选中的测量变量需要通过万用表模块进行观测。 2. 三相断路器模块三相断路器模块外部控制方式、带缓冲电路和不带缓冲电路的三相断路器模块图标如图6-35所示。第6章 电力系统稳态与暂态仿真 双击三相断路器模块，弹出该模块的参数对话框如图6-36所示。该对话框中含有以下参数：(1) “断路器初始状态”(Initial status of breakers)下拉框：断路器三

59、相的初始状态相同，选择初始状态后，图标会显示相应的切断或者投合状态。(2) “A相开关”(Switching of phase A)复选框：选中该复选框后表示允许A相断路器动作，否则A相断路器将保持初始状态。(3) “B相开关”(Switching of phase B)复选框：选中该复选框后表示允许B相断路器动作，否则B相断路器将保持初始状态。第6章 电力系统稳态与暂态仿真 图6-36 三相断路器模块参数对话框 第6章 电力系统稳态与暂态仿真 (4) “C相开关”(Switching of phase C)复选框：选中该复选框后表示允许C相断路器动作，否则C相断路器将保持初始状态。(5) “

60、切换时间 (Transition times)文本框：采用内部控制方式时，输入一个时间向量以控制开关动作时间。如果选中外部控制方式，该文本框不可见。(6) “外部控制”(External control of switching times)复选框：选中该复选框，断路器模块上将出现一个外部控制信号输入口。开关时间由外部逻辑信号(0或1)控制。(7) “断路器电阻”(Breaker resistance Ron)文本框：断路器投合时内部电阻(单位：W)。断路器电阻不能为0。第6章 电力系统稳态与暂态仿真 (8) “缓冲电阻”(Snubber resistance Rp)文本框：并联的缓冲电路中的

61、电阻值(单位：)。缓冲电阻值设为inf时，将取消缓冲电阻。 (9) “缓冲电容”(Snubber capacitance Cp)文本框：并联的缓冲电路中的电容值(单位：F)。缓冲电容值设为0时，将取消缓冲电容；缓冲电容值设为inf时，缓冲电路为纯电阻性电路。(10)“测量参数”(Measurements)下拉框：对以下变量进行测量。 “无”(None)：不测量任何参数。第6章 电力系统稳态与暂态仿真 “断路器电压”(Branch voltages)：测量断路器的三相终端电压。 “断路器电流”(Branch currents)：测量流过断路器内部的三相电流，如果断路器带有缓冲电路，测量的电流仅为

62、流过断路器器件的电流。 “所有变量”(Branch voltages and currents)：测量断路器电压和电流。选中的测量变量需要通过万用表模块进行观察。测量变量用“标签”加“模块名”加“相序”构成，例如断路器模块名称为B1时，测量变量符号如表6-3所示。第6章 电力系统稳态与暂态仿真 表表6-3 三相断路器测量变量符号三相断路器测量变量符号第6章 电力系统稳态与暂态仿真 图6-37 三相故障模块等效电路 3. 三相故障模块三相故障模块三相故障模块是由三个独立的断路器组成的、能对相相故障和相地故障进行模拟的模块。该模块的等效电路如图6-37所示。外部控制方式和内部控制方式下的三相故障模

63、块图标如图6-38所示。 第6章 电力系统稳态与暂态仿真 图6-38 三相故障模块图标 第6章 电力系统稳态与暂态仿真 双击三相故障模块，弹出该模块的参数对话框如图6-39所示。在该对话框中含有以下参数：(1) “A相故障”(Phase A Fault)复选框：选中该复选框后表示允许A相断路器动作，否则A相断路器将保持初始状态。(2) “B相故障”(Phase B Fault)复选框：选中该复选框后表示允许B相断路器动作，否则B相断路器将保持初始状态。(3) “C相故障”(Phase C Fault)复选框：选中该复选框后表示允许C相断路器动作，否则C相断路器将保持初始状态。(4) “故障电阻

64、”(Fault resistances Ron)文本框：断路器投合时的内部电阻(单位：)。故障电阻不能为0。第6章 电力系统稳态与暂态仿真 图6-39 三相故障模块参数对话框第6章 电力系统稳态与暂态仿真 (5) “接地故障”(Ground Fault)复选框：选中该复选框后表示允许接地故障。通过和各个开关配合可以实现多种接地故障。未选中该复选框时，系统自动设置大地电阻为106 。(6) “大地电阻”(Ground resistance Rg)文本框：接地故障时的大地电阻(单位：)。大地电阻不能为0。选中接地故障复选框后，该文本框可见。 (7) “外部控制”(External control

65、of fault timing)复选框：选中该复选框，三相故障模块上将增加一个外部控制信号输入端。开关时间由外部逻辑信号(0或1)控制。 第6章 电力系统稳态与暂态仿真 (8) “切换状态”(Transition status)文本框：设置断路器的开关状态，断路器按照该文本框设置状态进行切换。采用内部控制方式时，该文本框可见。断路器的初始状态默认为与该文本框中第一个状态量相反的状态。(9) “切换时间”(Transition times)文本框：设置断路器的动作时间，断路器按照该文本框设置的时间进行切换。(10) “断路器初始状态”(Initial status of fault)文本框：设置

66、断路器的初始状态。采用外部控制方式时，该文本框可见。(11) “缓冲电阻”(Snubber resistance Rp)文本框：并联的缓冲电路中的电阻值(单位：)。缓冲电阻值设为inf时，将取消缓冲电阻。第6章 电力系统稳态与暂态仿真 (12) “缓冲电容”(Snubber capacitance Cp)文本框：并联的缓冲电路中的电容值(单位：F)。缓冲电容值设为0时，将取消缓冲电容；缓冲电容值设为inf时，缓冲电路为纯电阻性电路。(13) “测量参数”(Measurements)下拉框：对以下变量进行测量。 “无”(None)：不测量任何参数。 “故障电压”(Branch voltages)

67、：测量断路器的三相端口电压。 “故障电流”(Branch currents)：测量流过断路器的三相电流，如果断路器带有缓冲电路，测量的电流仅为流过断路器器件的电流。第6章 电力系统稳态与暂态仿真 “所有变量”(Branch voltages and currents)：测量断路器电压和电流。选中的测量变量需要通过万用表模块进行观察。测量变量用“标签”加“模块名”加“相序”构成，例如三相故障模块名称为F1时，测量变量符号如表6-4所示。第6章 电力系统稳态与暂态仿真 表表6-4 三相故障模块测量参数符号三相故障模块测量参数符号第6章 电力系统稳态与暂态仿真 6.3.2 暂态仿真分析暂态仿真分析【

68、例【例6.6】 线电压为300kV的电压源经过一个断路器和300km的输电线路向负荷供电。搭建电路对该系统的高频振荡进行仿真，观察不同输电线路模型和仿真类型的精度差别。解：解：(1) 按图6-40搭建仿真单相电路图，选用的各模块的名称及提取路径见表6-5。第6章 电力系统稳态与暂态仿真 图6-40 例6.6的仿真电路图第6章 电力系统稳态与暂态仿真 表表6-5 例例6.6仿真电路模块的名称及提取路径仿真电路模块的名称及提取路径第6章 电力系统稳态与暂态仿真 (2) 设置模块参数和仿真参数。并联RLC模块Z_eq的参数设置如图6-41所示。断路器模块Breaker的参数设置如图6-42所示。其余

69、元件参数与例4.4相同，仿真参数的设置也与例4.4相同。仿真结束时间取为0.02 s。第6章 电力系统稳态与暂态仿真 图6-41 例6.6的Z_eq参数设置 第6章 电力系统稳态与暂态仿真 图6-42 例6.6的Breaker参数设置 第6章 电力系统稳态与暂态仿真 (3) 不同输电线路模型下的仿真。按例6.4的方法，设置线路为1段PI形电路、10段PI形电路和分布参数线路，把仿真得到的V2处电压分别保存在变量V21、V210和V2d中，并画出对应的波形如图6-43所示。由图6-43可见，断路器在0.005 s合闸时，系统中产生了高频振荡。其中由1段PI形电路模块构成的系统未反映高于206 H

70、z的振荡(见例4.4)，由10段PI形电路模块构成的系统较好地反映了这种高频振荡，分布参数线路由于波传导过程在断路器合闸后存在1.03 ms的时间延迟。第6章 电力系统稳态与暂态仿真 图6-43 例6.6不同线路模型下的电压波形比较 第6章 电力系统稳态与暂态仿真 (4) 不同仿真类型下的仿真。用10段PI型输电线路按例6.4的方法，执行以下三种仿真： 连续系统仿真，Ts=0 s； 离散系统仿真，Ts=25 s； 离散系统仿真，Ts=50 s。把仿真得到的V2处电压分别保存在变量V2c、V2d25和V2d50中，并画出对应的波形如图6-44所示。由图6-44可见，25 s步长下的仿真结果(短虚

71、线)与连续系统的结果(实线)很接近，而50 ms步长下的仿真结果(长虚线)已经小有误差了。三种仿真的运算时间分别为0.25 s、0.17s和0.15 s。因此，本例中选择25 s的步长不但仿真精度满足要求，还可以提高运算速度。第6章 电力系统稳态与暂态仿真 图6-44 例6.6不同仿真类型下的电压波形比较第6章 电力系统稳态与暂态仿真 图6-45 例6.7的系统图【例【例6.7】供电系统如图6-45所示，其中线路L的参数为长50 km，r=0.17 /km，x=0.402 /km。变压器T的参数为Sn=10 MVA，Vs%=10.5，KT=110/11。假定供电点电压Vi为106.5 kV，保

72、持恒定，当空载运行时变压器低压母线发生三相短路。试构建系统进行仿真，并观察短路电流周期分量、冲击电流大小。第6章 电力系统稳态与暂态仿真 解：解：(1) 理论分析。将供电点等效为理想电压源，同时忽略线路和变压器中的并联导纳，可得到线路电阻RL、线路电抗XL和变压器电抗XT分别为RL=500.17=8.5 (6-2) XL=500.402=20.1 (6-3)(6-4) 第6章 电力系统稳态与暂态仿真 因此，变压器低压侧短路电流周期分量的幅值Ipm为(6-5) 第6章 电力系统稳态与暂态仿真 冲击电流iim为(2) 按图6-46搭建仿真电路图，选用的各模块的名称及提取路径见表6-6。(6-6)

73、第6章 电力系统稳态与暂态仿真 图6-46 例6.6的仿真电路图 第6章 电力系统稳态与暂态仿真 表表6-6 例例6.7仿真电路模块的名称及提取路径仿真电路模块的名称及提取路径第6章 电力系统稳态与暂态仿真 (3) 设置模块参数和仿真参数。简化同步电机参数设置如图6-47所示。三相故障模块Breaker在0.02 s时三相合闸，对应的参数设置如图6-48所示。第6章 电力系统稳态与暂态仿真 图6-47 例6.7的简化同步电机模块参数设置 第6章 电力系统稳态与暂态仿真 图6-48 例6.7的三相故障模块参数设置 第6章 电力系统稳态与暂态仿真 并联RLC负荷为有功功率负荷，负荷大小为5 MW，

74、其余元件参数按题目已知条件设置。选择ode23tb算法，仿真结束时间取为0.6 s。仿真开始前，利用Powergui模块对电机进行初始化设置，初始化后，与简化同步电机模块输入端口相连的两个常数模块Pm和VLLrms的参数被自动设置为4.70961e6和106665。(4) 仿真。开始仿真，得到变压器低压侧的a相电流如图6-49所示。第6章 电力系统稳态与暂态仿真 图6-49 例6.7的仿真波形图第6章 电力系统稳态与暂态仿真 由图6-49可见，短路电流周期分量的幅值为5.85 kA，冲击电流为10.5 kA，和理论计算分别存在0.85%和2.98%的误差。这是由于实际仿真中，供电点并不是理想电

75、压源，发生短路后，供电点电压将降低，因此计算得到的短路电流和冲击电流值偏大。利用SimPowerSystems/Extra Library/Measurements子库中的“FFT模块”(Fourier)和“三相序分量模块”(3-Phase Sequence Analyzer)还可以进行电流直流和倍频分量的分析或者正序、负序和零序的分析。限于篇幅，本节不再给出各序分量和各次谐波分量的电流波形图，读者可以动手试试，看看结果是否和理论一致。第6章 电力系统稳态与暂态仿真 6.4 电力系统机电暂态仿真电力系统机电暂态仿真当电力系统受到大的扰动时，表征系统运行状态的各种电磁参数都要发生急剧的变化。但是

76、，由于原动机调速器具有较大的惯性，它必须经过一定时间后才能改变原动机的功率。这样，发电机的电磁功率与原动机的机械功率之间便失去了平衡，于是产生了不平衡转矩。在不平衡转矩作用下，发电机开始改变转速，使各发电机转子间的相对位置发生变化(机械运动)。发电机转子相对位置，即相对角的变化，反过来又将影响到电力系统中电流、电压和发电机电磁功率的变化。第6章 电力系统稳态与暂态仿真 所以，由大扰动引起的电力系统暂态过程，是一个电磁暂态过程和发电机转子间机械运动暂态过程交织在一起的复杂过程。如果计及原动机调速器、发电机励磁调节器等调节设备的暂态过程，则过程将更加复杂。精确地确定所有电磁参数和机械运动参数在暂态

77、过程中的变化是困难的，对于解决一般的工程实际问题往往也是不必要的。通常，暂态稳定性分析计算的目的在于确定系统在给定的大扰动下发电机能否继续保持同步运行。因此，只需研究表征发电机是否同步的转子运动特性，即功角d随时间变化特性便可以了。这就是通常说的机电暂态过程，即稳定性问题。第6章 电力系统稳态与暂态仿真 本节将对一个含两台水轮发电机组的输电系统进行暂态稳定性的仿真演示。为提高系统的暂态稳定性和阻尼振荡的能力，该系统中配置了静止无功补偿器(SVC)以及电力系统稳定器(PSS)。打开SimPowerSystems库的demo子库中的模型文件power_svc_pss，可以直接得到如图6-50所示的

78、仿真系统，以文件名circuit_pss另存。这样，用户可对该原始模型进行进一步的调整。该系统的具体实现方法和参数设置可以参考有关文献10。对于初学者来说，本节有一定的难度。第6章 电力系统稳态与暂态仿真 图6-50 电力系统暂态稳定性分析的仿真系统图第6章 电力系统稳态与暂态仿真 6.4.1 输电系统的描述输电系统的描述图6-50是一个简单的500 kV输电系统图。图中，一个1000 MVA的水轮发电厂(M1)通过500 kV、700 km输电线路与5000 MW的负荷中心相连，另一容量为5000 MVA的本地发电厂(M2)也向该负荷供电。为了提高故障后系统的稳定性，在输电线路中点并联了一个

79、容量为200 Mvar的静止无功补偿器。两个水轮发电机组均配置水轮机调速器、励磁系统和电力系统稳定器。单击并进入“涡轮和调速器”(Turbine & Regulators)子系统，其结构如图6-51所示。第6章 电力系统稳态与暂态仿真 图6-51 “涡轮和调速器”子系统结构图第6章 电力系统稳态与暂态仿真 该子系统中，与励磁系统相连的稳定器模块有两种类型：一种是“普通PSS”(Generic Power System Stabilizer)模块，另一种是“多频段PSS”(Multi-Band Power System Stabilizer)模块。这两种稳定器模块都可以从SimPowerSyst

80、ems/Machines库中直接提取。通过手动设置图6-50左下方的“开关”模块可以选择不同的PSS，或者将系统设置为不含PSS的工作状态。第6章 电力系统稳态与暂态仿真 图6-50中的SVC模块是SimPowerSysterms/Phasor Elements库中的相量模块。打开SVC模块的参数对话框，在“显示”(Display)下拉框中选择“功率数据”(Power data)选项，将显示功率数据参数对话框(见图6-52(a)，确定SVC的额定容量是+/200Mvar；若在“显示”(Display)下拉框中选择“控制参数”(Control parameters)选项，将显示控制参数对话框(见

81、图6-52(b)，在该窗口中，可以选择SVC的运行模式为“电压调整”(Voltage regulation)或“无功控制”(Var control)，默认设置为“无功控制”模式。若不希望投入SVC，直接将电纳设置为Bref=0即可。第6章 电力系统稳态与暂态仿真 图6-52 SVC模块参数对话框(a) SVC功率数据；(b) SVC控制参数第6章 电力系统稳态与暂态仿真 图6-50中的母线B1上连接有一个三相故障模块。通过该故障模块设置不同类型的故障，可观测PSS和SVC对系统稳定性的影响。仿真开始前，打开Powergui模块参数对话框，选中“相量法分析”单选框以加快仿真速度。点击Powerg

82、ui模块的“潮流计算和电机初始化”按键进行初始化设置。将发电机M1定义为PV节点(V=13 800 V，P=950 MW)，发电机M2定义为平衡节点(V=13 800 V，a相电压相角为0，估计要送出的有功功率为4000 MW)。潮流计算和初始化工作完成后，两个发电机参数对话框中的初始条件、两个发电机输入端口的参考功率都被自动更新，其中Pref1=0.95 p.u.(950 MW)，Pref2=0.8091 p.u.(4046 MW)。第6章 电力系统稳态与暂态仿真 更新后发电机的初始状态如图6-53。进入“涡轮和调速器”子系统，可以看见两个励磁系统输入端口上的参考电压被自动更新为Vref=V

83、ref1=1.0 p.u.。第6章 电力系统稳态与暂态仿真 图6-53 更新后的发电机初始参数(a) 发电机M1；(b) 发电机M2第6章 电力系统稳态与暂态仿真 6.4.2 单相故障单相故障本节将对不使用SVC时的单相故障进行仿真，并观测系统的暂态稳定性。电力系统中发电机经输电线路并列运行时，在扰动下会发生发电机转子间的相对摇摆，并在缺乏阻尼时引起持续振荡。此时，输电线路上功率也会发生相应振荡。由于其振荡频率很低，一般为0.22.5 Hz，故称为低频振荡。电力系统低频振荡在国内外均有发生，常出现在长距离、重负荷输电线路上，在采用现代快速、高顶值倍数励磁系统的条件下更容易发生。这种低频振荡可以

84、通过电力系统稳定器得到有效抑制。第6章 电力系统稳态与暂态仿真 此外，从理论分析上可知，当转子间相角差为90时，发电机输出的电磁功率达到最大值。若系统长期在功角大于90的状况下运行，电机将失去同步，系统不稳定。设置SVC的参数Bref=0，即不使用SVC。设置三相故障模块在0.1 s时发生a相接地故障，0.2 s时清除故障。分别对投入普通PSS、投入多频段PSS、退出PSS三种情况进行暂态仿真。将这三种情况下的仿真结果叠加比较，如图6-54所示。图中波形从上到下依次为转子间相角差、电机转速和SVC端口上的正序电压幅值。第6章 电力系统稳态与暂态仿真 图6-54 单相接地故障时的暂态仿真波形第6

85、章 电力系统稳态与暂态仿真 从图中可见，在故障期间，由于电机M1的电磁功率小于机械功率，因此电机M1的转速增大。未安装PSS时，转子相角差在3.8 s时超过90，并且振荡失稳，因此系统是暂态不稳定的。普通PSS和多频段PSS下，最大转子相角差分别为62.5和57.8，5 s时，相角差在53左右重新达到平衡，因此系统具有暂态稳定性。本例中，普通PSS有效抑制了0.6 Hz的低频振荡，但对0.025 Hz的低频振荡作用不明显。若将仿真时间延长到50 s，则可以很清楚地观察到故障清除后电机发生了0.025 Hz的低频振荡，而多频段PSS有效抑制了0.6 Hz和0.025 Hz的低频振荡。第6章 电力

86、系统稳态与暂态仿真 可见，发生单相接地故障后，尽管未使用SVC，电机之间仍然能够重新恢复同步运行，因此具有暂态稳定性。故障清除后，0.6 Hz的低频振荡迅速衰减。6.4.3 三相故障三相故障本节对SVC和PSS均投入使用时的三相故障进行仿真，并观测系统的暂态稳定性。设置三相故障模块在0.1 s时发生三相接地故障，0.2 s时故障消失。打开SVC参数对话框，选择显示“控制参数”参数对话框，将SVC的运行模式改为“电压调整”，将“参考电压Vref”(Reference voltage Vref)文本框中的值改为1.009(1.009 p.u.为未投入SVC时SVC端口的电压稳态值)，其余参数不变。

87、开始仿真，仿真结果如图6-55所示。第6章 电力系统稳态与暂态仿真 图6-55 三相接地故障时的暂态仿真波形 第6章 电力系统稳态与暂态仿真 图中波形从上到下依次为转子间相角差、电机转速、SVC端口的正序电压、从SVC端口看入的等效电纳。为了方便比较，将投入普通PSS、未投入SVC(Bref=0)时的三相接地故障仿真波形叠加到图6-55中。由图6-55可见，未使用SVC时，两个电机在0.3 s时迅速单调失去同步，电机转速单调增大，SVC等效电纳为0，表示不向系统吸收无功，也不向系统发送无功。安装SVC后， SVC的等效电纳在正值和负值间波动，正的电纳表示向系统发送无功，负的电纳表示从系统吸收无

88、功。转子间相角差虽然有短暂时间超过了90，但最终以衰减振荡的形式稳定在53附近。因此，尽管发生了最为严重的三相接地故障，但系统仍然具有暂态稳定性。第6章 电力系统稳态与暂态仿真 题6-1图习习 题题6-1 按题6-1图设计交流电压源的叠加电路，分析线路首端电压的变化情况。两个单相交流电压源分别为v1=100 sin(120pt+p/6) V和v2=75 sin(100pt+p/3) V。第6章 电力系统稳态与暂态仿真 题6-2图6-2 按题6-2图设计一个5次谐波滤波器电路并进行稳态运行分析。第6章 电力系统稳态与暂态仿真 题6-3图6-3 按题6-3图建模，利用指令分析该电路的结构特征和状态

89、方程。第6章 电力系统稳态与暂态仿真 题6-4图6-4 按题6-4图建模，利用park变换模块将三相电压从abc坐标系转换为dq0坐标系。6-5 按题6-5图建模，利用三相序分量模块分析A相接地后的正序、负序、零序分量。第6章 电力系统稳态与暂态仿真 题6-5图第6章 电力系统稳态与暂态仿真 6-6 将题6-5中的“三相电压源模块”(Three-Phase Source)替换为“三相可编程电压源模块”(Three-Phase Programmable Voltage Source)。在电压源的基频分量中叠加谐波分量，其中3次谐波分量的大小为基频分量的13%，7次谐波分量的大小为基频分量的5%。利用傅里叶分析模块分析输电线路A相接地前后线路两端的电压谐波含量。6-7 利用Powergui模块中的“稳态电压电流分析”功能分析题6-7图电路。第6章 电力系统稳态与暂态仿真 题6-7图第6章 电力系统稳态与暂态仿真 题6-8图6-8 按题6-8图建模，利用Powergui模块中的“相量分析”功能分析输电线路两端的电压和线路中电流的变化。第6章 电力系统稳态与暂态仿真 结束语结束语若有不当之处，请指正，谢谢！若有不当之处，请指正，谢谢！