基于MATLAB的输电线路故障行波仿真平台.doc

基于MATLAB的输电线路故障行波仿真平台马永明1，陈 平1，刘万超2（1.山东理工大学 电气与电子工程学院，山东淄博255049；2.山东省电力工程咨询院有限公司 电控部，山东 济南 250013） 收稿日期：2009-11-19作者简介：马永明（1984-），男，硕士. E-mail: 摘要：为了便于对输电线路故障行波的学习和研究，提出了基于MATLAB/SIMULINK的输电线路故障行波仿真方法，详细介绍了SIMULINK提供的电力系统工具箱（PSB）中的部分模块的功能和使用方法及线路故障行波的提取方法, 在此基础上建立了基于分布参数模型的输电线路故障行波仿真平台。最后，以输电线路的单相短路故障为例，应用该平台和MATLAB编程实现了对输电线路故障的仿真以及故障行波的提取。关键词：行波；MATLAB/SIMULINK；仿真平台；分布参数中图分类号：TM711 文献识别码：ATransmission Line Traveling-wave Simulation Platform Based on MATLAB/SIMULINKMA yong-ming1, CHEN ping1, LIU wan-chao2(1. School of Electrical and Electronic Engineering, Shandong University of Technology, Zibo 255049, China;2. Shandong Electric Power Engineering Consulting Institute LTD., Jinan 250013, China)Abstract: In order to facilitate the learning and research of transmission line traveling-wave, a simulation method of transmission line traveling-wave based on MATLAB/SIMULINK is proposed, the function and use of some blocks of PSB (Power System Blockset) in SIMULINK and the calculation of traveling-wave are discussed, on this base, the transmission line traveling-wave simulation platform based on distributed parameter model is established. At last, take phase A fault of transmission line for example, it realized the simulation of transmission line fault and the calculation of traveling-wave with the help of the platform and MATLAB programming. Keywords: traveling-wave; MATLAB/SIMULINK；simulation platform; distributed parameter随着经济的快速发展，电力系统的规模越来越大，其内部结构也日益复杂。为了保证电网运行的安全性和经济性，对超高压电网，人们已要求能够精确测出故障(尤其是重复性故障和永久性故障)发生的确切位置，以便迅速排除故障，尽快恢复系统正常运行，最大限度地减少经济损失1。在故障测距方法中，由于行波测距与传统的阻抗法和故障分析法相比具有测距准确可靠且经济方便等优点，日益成为关注的焦点和热点。为了便于对输电线路故障行波的学习和研究，利用功能强大的MATLAB/SIMULINK建立一个输电线路故障行波仿真平台就变得十分必要而且富有意义。1 MATLAB/SIMULINK与故障行波仿真1.1 MATLAB/SIMULINK简介MATLAB仿真软件是美国Mathworks公司于1998年推出的,在公司以及许多专家的努力下，经过多次扩充和修改，现已发行到MATLAB7.0以上版本，已成为当今国际上最流行的控制系统辅助设计（CACSD）技术的语言和软件工具。MATLAB具有3大特点：功能强大；界面友好，语言自然；开放性强。正是由于MATLAB的这些优点，使它获得了对应用学科（特别是边缘学科和交叉学科）的极强适应力，并很快成为应用学科计算机辅助分析设计、仿真、教学乃至科技文字处理中不可缺少的基础软件，成为欧美高等院校、科研机构教学与科研必备的基本工具2。MATLAB提供的动态仿真工具SIMULINK是众多仿真软件中功能最强大、最优秀、最容易使用的一种，可用来对动态系统进行建模、仿真和分析。SIMULINK为用户提供了用方框图进行建模的模型接口，它与传统仿真软件包利用微分方程和差分方程建模的方法相比，具有更直观、方便、灵活的优点。用户只需用鼠标选中各模块，并放进新建的模块图，最后将各模块连接起来，即完成对特定系统模型的创建3。1.2 SIMULINK的自定义模块功能虽然SIMULINK的模块库提供了丰富的模块，但是，并非所有的数学模型都能利用SIMULINK 模块轻易搭建起来，有的模型可能更适合于MATLAB 或C 语言的编写方法。因此，用户仍然需要定义自己的模块，来实现特殊算法、特定功能的子系统以及仿真过程的图形动画等。S-函数正好弥补了这一不足，它提供了扩展SIMULINK 功能的有力工具，并且支持参数传递。用户可以建立一个能和SIMULINK 模块库中的模块一起使用的新模块，利用封装产生一个与S-函数文件相对应的对话框和模块图标，不断扩充SIMULINK 的仿真功能4。这使得SIMULINK获得了更为广阔的发展空间。1.3 故障行波仿真当电网发生故障时，根据叠加原理,在故障瞬间,相当于在故障点突然附加一个与故障前电压大小相等、相位相反的虚拟电源,虚拟电源会产生向线路两端运动的电压、电流行波。行波实际上就是输电线路发生故障时产生的非工频暂态分量。由于行波的传播速度接近光速,且不受故障点电阻、线路结构及互感器变换误差等因素的影响,因此在故障测距中与传统的阻抗法和故障分析方法相比具有较高的测量精度5，并且在电网故障测距中日益得到广泛应用。MATLAB/SIMULINK不仅具有强大的绘图和数据处理功能，而且提供了完备的电力系统工具箱（Power System Blockset以下简称PSB），PSB提供了一种类似电路建模的方式进行模型绘制，在仿真前自动将其变换成状态方程描述的系统形式，然后在SIMULINK下进行仿真分析6。因此，可以在MATLAB/SIMULINK环境下，建立输电线路故障行波仿真平台，对输电线路的各种故障进行仿真，并对线路发生故障后的电压、电流数据和故障时产生的行波进行分析和研究。2 输电线路故障行波仿真平台2.1仿真平台的构建打开MATLAB软件，单击File菜单选择NewModel命令就可以打开一个新的仿真平台编辑窗口； 在MATLAB运行界面中键入powerlib命令，就可以打开电力系统元件库Library；在电力系统元件库中双击所需模块库，就会弹出一个新的窗口，在这个窗口中有该模块库中所包含各种元件模型。选择你所需要的元件模型，并将其拖入已打开的仿真平台编辑窗口中；当所有需要的元件模型都拖入到编辑窗口中后，用线将不同的元件模型按一定的秩序连接起来，就建成了一个仿真平台。2.2 仿真平台的基本框架按照上述步骤，建立了一个基于MATLAB 7.1的输电线路故障行波仿真平台，其基本框架如图1所示。此平台是由3个三相电源和4段分布参数输电线构成的环形网络。其中Three-Phase V-I Measurement模块为处于检测点的三相电压电流测量模块； Scope V和Scope I用于观察检测点的三相电压、电流波形；Three-Phase Fault为处于故障点的线路故障设置元件；其中故障点的位置可通过改变输电线的长度来调整。2.3 仿真平台的应用首先，根据不同的仿真要求设置各元件模块的参数。参数设置是仿真过程中的重要一环，参数设置的好坏直接影响仿真结果的准确性。然后设置仿真的起止时间、仿真步长和仿真算法等；最后，点击仿真启动按钮，待仿真结束后就可以在示波器Scope V和Scope I看到检测点的电压、电流波形，同时仿真后得到的电压、电流数据可存储到MATLAB的Workspace中，以便于对仿真数据的绘图、分析和处理。3 仿真平台各模块的功能与实现3.1 三相电源三相电源（Three-Phase Source）是仿真平台中的重要元件，它的运行特性对整个系统的运行状态起着决定性的影响。双击Three-Phase Source元件，则会弹出三相电源参数设置对话框。在该对话框的参数区域中包含7个选项：线电压有效值（Phase-to-phase rms voltage）；A相相角（Phase angle of phase A）为电源A相的相位角；频率图1 输电线路故障行波仿真平台框图（Frequency）；内部连接方式（Internal connection）；短路阻抗值（Specify impedance using short-circuit level）；三相电源内阻（Source resistance）和三相电源电感（Source inductance）。3.2 分布参数输电线分布参数输电线（Distributed ParametersLine）的参数区域中有7各选项：相数（Number of phases）；线路频率（Freqency used for RLC specification）；每公里输电线的电阻（Resistance per unit length）；每公里输电线的电感（Reductance per unit length）；每公里输电线的电容（Capacitance per unit length）；输电线路的长度（Line length）和测量量选择（Measurements）。其中，每公里输电线的电阻、电感和电容值都分为正序和零序两部分：R1、L1和C1表示线路的正序电阻、电感和电容；R0、L0和C0表示线路的零序电阻、电感和电容。3.3 三相电压电流测量模块三相电压电流测量模块（Three-Phase V-I Measurement）用来测量检测点的三相电压、电流。其参数区域分为两个部分：测量电压选项（Voltage measurement）和测量电流选项（Current measurement）。Voltage measurement有3个下拉选项：no, phase-to-ground和phase-to-phase，用来选择是否测量检测点的电压以及测量电压的类型（相电压或者线电压）；同样Current measurement也有2个下拉选项：yes和no，来选择是否测量检测点的电流值；勾选Use a label选项，可以在其下出现输入框中输入字符对被测量的电压、电流信号进行标记；勾选Voltage in pu 和Current in pu选项，可以设置电压和电流的标幺值。3.4 三相线路故障模块三相线路故障模块（Three-Phase Fault）是仿真平台中最重要的模块，也是选项最多、参数设置最复杂的一个模块，可用来设置故障点的故障类型等参数。该模块参数区域中共有12个选项：通过勾选Phase A Fault、Phase B Fault和Phase C Fault中的一相或者几相来选择短路故障相；Fault resistance用来设置短路阻抗；Ground Fault选项用来选择短路故障是否为短路接地故障；Ground resistance用来设置接地阻抗的大小；选中External control of fault timing可以添加控制信号来控制该模块故障的启动和停止；Transition status 和Transition times用来设置故障开关的状态及相应的动作时间； Snubbers resistance和Snubbers capacitance用来设置过渡电阻和过渡电容；Measurements同样是用来选择测量量。所以通过三相线路故障模块可以设置故障的类型、起止时间和被测量参数等。4 基于仿真平台的正反向行波的提取如果把行波从故障点指向线路两端的传播方向规定为正方向，则行波就可分正向行波和反向行波。正向行波和反向行波的提取对行波测距的研究具有重要意义。利用该仿真平台对线路故障进行仿真后，在MATLAB的wokspce中会得到以变量形式存储的检测点的三相电压和三相电流数据，根据该三相电压和三相电流数据就可以提取故障发生时电压的正向行波和反向行波，具体提取方法如下：提取三相电压和三相电流的故障分量，用故障后一段时段内的三相电压、电流值减去故障前相应的一时段内的三相电压、电流值，就得到了三相电压、电流的暂态量和；将三相电压、电流的暂态量和进行克拉克模量变换，就会得到电压、电流的模量值和，其计算过程如下： （1） （2）其中克拉克变换矩阵计算电压1模正向行波和反向行波，计算过程如下： （3） （4）其中、分别表示和的1模分量，也就是由矩阵和的第2行构成的矩阵；，、分别为每公里输电线路的正序电感和正序电容7。可用MATLAB语言将上述算法编写成程序，仿真之后直接运行该程序就可以求出电压1模正向行波和反向行波并绘制出相应的波形图。5. 应用实例不妨应用该仿真平台对输电线路的单相短路故障进行仿真，可设为A相短路。在三相线路故障模块（Three-Phase Fault）的参数区域中选中Phase A Fault和Ground Fault，Transition status为1 0，Transition times为0.035 0.080；3个三相电压源（Three-Phase Source）的相电压有效值（Phase-to-phase rms voltage）均为500kV，A相相角（Phase angle of phase A）依次为0度、30度和60度；在示波器（Scope）的Data history参数区域中取消Limit data pionts to last选项，而勾选Save data to workspace选项，并将变量的存储格式设置为数组(Array)格式。由于篇幅原因，其它模块参数的设置不再赘述。设置好模块参数后，再设置仿真参数：仿真的起止时间分别为0.0s和0.10s；仿真采用固定步长（Fixed-step）1e-6；仿真算法为ode3。然后点击仿真启动按钮，仿真完成后就可以在示波器（Scope）中看到检测点的三相电压、电流波形，如图2和图3所示。图2 检测点A相、B相和C相的电压波形图3 检测点A相、B相和C相的电流波形 通过检测点的三相电压和电流波形可以看出，仿真后得到的电压、电压波形符合线路发生A相短路故障后的电压、电流特征，从而可以说明仿真平台的正确性。运行事先用MATLAB语言编写的电压1模正向行波和反向行波的提取程序，就可以得到从仿真后的三相电压、电流数据中提取到电压1模正向行波和反向行波，如图4所示。图 4 电压1模正向行波和反向行波（1）表示检测点1模电压正向行波（