电力系统电磁暂态概述.docx

第1章 电力系统电磁暂态概述1.1 电力系统电磁暂态现象1.2 电力系统电磁暂态分析的目的1.3 电力系统电磁暂态研究的方法1.4 电力系统电磁暂态仿真的特点1.5 电力系统数字仿真思考及练习题1.1 电力系统电磁暂态现象21.2 电力系统电磁暂态分析的目的41.3 电力系统电磁暂态研究的方法51.4 电力系统电磁暂态的特点71.4.1 频率范围广71.4.2 元件模型因计算目的而异81.4.3 行波现象和分布参数101.4.4 非线性元件和开关操作161.4.5 元件参数的频率特性171.4.6 时间跨度的要求181.5 电力系统数字仿真181.5.1 电力系统数字仿真的分类181.5.2 电力系统数字仿真的优点201.5.3 电力系统数字仿真软件21第1章 电力系统电磁暂态概述1.1 电力系统电磁暂态现象 电力系统稳态运行时，发电厂发出的功率及用户所需要的功率及电网中损耗的功率相平衡，系统的电压与频率都是稳定的。但电力系统在运行过程中常常会发生故障或需要进展操作，常见的电力系统故障有：雷击电力设备等雷害故障，短路、接地故障与谐振等电气故障，断线等机械故障。常见的电力系统操作有：1断路器的投切操作，如合空载线路、合空载变压器、切空载线路、重合闸、甩负荷等。2隔离开关的投切操作，如母线投切等。电力系统发生故障或进展操作时，系统的运行参数发生急剧变化，系统的运行状态有可能急促地从一种运行状态过渡到另一种运行状态，也有可能使正常运行的电力系统局部甚至全部遭到破坏，其运行参数大大偏离正常值，如不采取特别措施，系统很难恢复正常运行，这将给国民经济生产与人民生活带来严重的后果。电力系统运行状态的改变，不是瞬时完成的，而要经历一个过渡状态，这种过渡状态称为暂态过程。电力系统的暂态过程通常可以分为电磁暂态过程与机电暂态过程。电磁暂态过程指电力系统各元件中电场与磁场以及相应的电压与电流的变化过程，机电暂态过程指由于发动机与电动机电磁转矩的变化所引起的电机转子机械运动的变化过程。虽然电磁暂态过程与机电暂态过程同时发生并且相互影响，但由于现代电力系统规模的不断扩大，构造愈益复杂，需要考虑的因素繁多，再加上这两个暂态过程的变化速度相差很大，要对它们统一分析是十分复杂的工作，因此在工程上通常近似地对它们分别进展分析。例如，在电磁暂态过程分析中，由于在刚开场的一段时间内，系统中的发电机与电动机等转动机械的转速由于惯性作用还来不及变化，暂态过程主要决定于系统各元件的电磁参数，故常不计发动机与电动机的转速变化，即忽略机电暂态过程。而在静态稳定性与暂态稳定性等机电暂态过程分析中，转动机械的转速已有了变化，暂态过程不仅及电磁参数有关，而且还及转动机械的机械参数转速、角位移有关，分析时往往近似考虑或甚至忽略电磁暂态过程。只在分析由发动机轴系引起的次同步谐振现象、计算大扰动后轴系的暂态扭矩等问题中，才不得不同时考虑电磁暂态过程与机电暂态过程。下面以一个简单开关接通RL电路的例子，以便获得对在电力系统暂态时起关键作用的物理过程的充分了解。一个正弦波电压接通到一个电感及电阻串联的电路上，如图1-1所示。这实际上是一个高压断路器闭合到短路的输电线路或短路的电缆的最简单单相表示法。电压源代表连接的同步发电机的电动势。电感L包括发电机的同步电感、电力变压器的漏电感及母线、电缆及输电线的电感，电阻R表示供电电路的电阻损耗。图1-1 正弦波电压源接通到RL串联电路假设时合闸，应用基尔霍夫电压定律，得到电路方程 1-1该方程为一阶常系数、线性、非齐次常微分方程，其解就是合闸电路的全电流，它由两局部组成：稳态分量与暂态分量，即其中稳态分量为 1-2暂态分量，也就是合闸电流的自由分量,记为 1-3式中，为特征方程的根；，为暂态分量电流衰减的时间常数；C为由初始条件决定的积分常数。在开关闭合之前，电感L中的磁通为0，根据磁通守恒定律，在闭合的瞬间，即由此得到 1-4从而得到合闸的全电流表达式为 1-5式1-5中的暂态分量含有衰减项，也称为直流分量，其系数为常数，数值大小取决于电流合闸瞬间，在为其中时，直流分量为0，电流立即进入稳态，换言之，不存在暂态振荡过程。但当开关闭合电路不在为其中时，合闸过程将引起电磁暂态过程，在为时，暂态过程将到达最大电流，如图1-2所示。图1-2 开关合闸的暂态过程电流波形1.2 电力系统电磁暂态分析的目的电磁暂态过程分析的主要目的在于分析与计算故障或进展操作后可能出现的暂态过电压与过电流，以便对电力设备进展合理设计。通常情况下，电力系统电磁暂态产生的过电压在确定设备绝缘水平中起决定作用，据此制定高电压试验电压标准，确定已有设备能否平安运行，并研究相应的限制与保护措施。此外，对于研究电力系统新型快速保护装置的动作原理及其工况分析，故障测距原理及定点方法以及电磁干扰等问题，也常需要进展电磁暂态过程分析。另外，调查事故原因，寻找对策；计算电力系统过电压发生概率，预测事故率；检查电气设备的动作责能，如断路器的暂态恢复电压与零点偏移；检查继电保护与平安自动装置的响应等，也离不开电磁暂态过程的计算与模拟。电磁暂态过程变化很快，一般需要分析与计算持续时间在ms级，甚至是s级以内的电压、电流瞬时值变化情况，因此，在分析中需要考虑元件的电磁耦合，计及输电线路分布参数所引起的波过程，有时甚至要考虑三相构造的不对称、线路参数的频率特性以及电晕等因素的影响。1.3 电力系统电磁暂态研究的方法为了保证电力系统运行的可靠性、平安性与经济性，在电力系统设计、运行、分析与研究中必须全面地了解实际系统的电磁暂态特性。目前，研究电力系统电磁暂态过程的手段有3种：1系统的现场实测方法。2应用暂态网络分析仪Transient Network Analyzer，简称TNA的物理模拟方法。3计算机的数字仿真或称数值计算方法。系统的现场实测方法是在实际的电力系统上直接进展试验与研究，六十年代之前经常要在实际电力系统进展短路、操作等试验，这种试验对电力系统的考验是真实与严格的，以确保电力系统运行的可靠性、平安性与经济性，但是系统的现场实测方法会对电力系统的正常运行与电气设备带来很大危害，短路点的电弧有可能烧坏电气设备，很大的短路电流通过设备会使发热增加，当持续一定时间后，可能使设备过热而损伤；很大的短路电流引起的电动力有可能使设备变形或遭到不同程度的破坏；操作试验过程中产生的过电压可能引起电气设备载流局部的绝缘损坏，加剧绝缘材料的老化。即便如此，实测对于研究电力系统电磁暂态过程仍是非常重要的，它一方面验证TNA及数字仿真的准确性，为系统平安运行提供依据；另一方面可以全面研究系统各类元件的参数特性，为TNA及数字仿真提供更准确的原始数据。系统的现场实测常常会遇到困难，有些困难甚至是不可能解决的，利用模型系统进展试验与分析就成为一种非常有效的途径。暂态网络分析仪就是一种用于研究电力系统动态特性的物理模型系统。TNA方法多用于模拟操作过电压与交流过电压的暂态现象，同时通过改变元件特性，TNA也可用来模拟更高频率下的暂态现象。它是在相似理论的指导下，把实际电力系统的各个局部，如同步发电机、变压器、输配电线路、电力负荷等按照相似条件设计、建造并组成一个电力系统模型，这样将一个高电压、大电流、体积庞大的电力系统，按照一定的比例转化为一个低电压、小电流、体积较小的模拟试验台，在模拟台中出现的电磁暂态现象，电压与电流的波形及它模拟的电力系统是一样的，用这种模型代替实际电力系统进展各种正常运行及故障状态的试验与分析。及系统的现场实测相比，TNA方法对电力系统的正常运行与电气设备不产生影响，为了缩小模拟装置的尺寸，节省电感元件与电容元件，减少模拟设备的昂贵费用，并考虑到现有的技术条件、模拟精度要求等，选择恰当的比例尺是非常重要的。TNA具有物理意义清晰，易于理解与使用的优点，可以屡次重复试验现象，便于观察与研究，北美不少大的电力公司都将TNA作为培训新员工的一种工具。随着现代电力系统的开展，电力系统的规模与复杂程度发生很大变化，采取物理模型的动态模拟方法受到很大限制。及此同时，数字计算机与数值计算技术飞速开展，数字计算机的性能价格比不断提高，出现了用数学模型代替物理模型的新型模型系统。电力系统数字仿真Digital Simulation of Power System就是将电力系统的电源、网络与负荷元件建立其数学模型，用数学模型在数字计算机上进展实验与分析的过程。电力系统数字仿真的主要步骤为建立各元件数学模型、建立数字仿真模型与进展仿真试验。建立数学模型是处理物理原型及数学模型之间的关系，建立数字仿真模型是处理数学模型及计算之间的关系。有些数学模型是利用数字计算机与模拟计算机的混合数学模型系统。电力系统数字仿真是一门新兴的技术科学，它的产生与开展是同现代科学技术开展分不开的，数字仿真及实际系统试验与动态物理模拟相比，不仅节省了大量的人力、物力与财力，而且不受外部条件的限制，几乎不受系统规模与时间跨度的约束，甚至不受各种暂态现象频率范围的限制理论上它可以对各类暂态过程进展计算，但是，它需要相关设备真实的频率特性，有时候，这种频率特性是很难得到的。具有无可比较的灵活性，能到达试验不可到达的广度与深度。譬如我国南北联网这样的课题，地理上相距数千公里，跨越了几个大电网，没有方法用试验来分析联网可能出现的问题，但通过数字仿真发现南北联网可能会出现低频震荡问题。今天实际系统的现场实测方法主要是为了建立数学模型，取得数学模型的参数。1.4 电力系统电磁暂态的特点 频率范围广电力系统中暂态现象的研究所涉及的频率范围广，从直流到大约50MHz的范围。高于系统频率的暂态现象通常涉及到电磁暂态，而低于系统频率的暂态现象主要涉及到机电暂态过程。表1-1给出了多种暂态现象的起因以及它们通常的频率范围。表1-1电力系统暂态的起因及频率范围起因频率范围投入变压器时的铁磁谐振DC0.1Hz 1kHz甩负荷0.1Hz 3 kHz故障去除50/60 Hz 3 kHz故障发生50/60 Hz 20 kHz线路充电50/60 Hz 20 kHz线路重合闸DC50/60 Hz 20 kHz断路器端部故障BTF50/60 Hz 20 kHz短路故障50/60 Hz 100 kHz断路器屡次重燃10 kHz 1 MHz雷击10 kHz 3 MHzGIS故障与隔离开关操作100 kHz 50 MHz网络中每个元件的模拟都要及所研究的特定暂态现象的频率范围相符合。当所研究现象的频率大于1 MHz时，如GIS中由于隔离开关操作所引起的快速暂态现象，那么不仅在母线上产生波的传播，而且施加在变压器、支柱绝缘子以及在某些情况下管形母线上的弯管处，它们非常小的电容与电感对模拟结果都将产生非常重要的影响。表1-1中所列电磁暂态现象的频率范围可以分成4组，对应于各暂态现象的频率范围之间存在着重叠，图1-3是国际大电网会议CIGRE对各种过电压的频率分类；各类的频率范围是及其所表示的过电压波形的实际陡度相关的。研究者必须清楚自己的研究对象所在的频率范围，确定被模拟设备的频率特性，只有这样，才能得到满意的电磁暂态分析计算结果。图1-3 各种过电压的频率范围 通常，频率越高，所考虑的现象(如过电压)在时间上空间上的衰减越快，因此所考虑的物理范围(模拟范围)越小，模拟时间越短。相对地，在工频或及此接近的频率领域，为了掌握现象的性质，需要大范围长时间的模拟。图1-4表示电力系统数字仿真中各种计算所考虑现象的时间幅度与计算涉及的系统规模。图1-4 电力系统各种现象的变化速度与计算范围 元件模型因计算目的而异电力系统由各种不同的元件所组成，元件的动态特性对于系统的暂态过程有直接的影响。为此，首先需要研究各元件的动态特性，建立它们的数学模型。在此根底上，根据系统的具体构造，即各元件之间的相互关系，组成全系统的数学模型，然后采用适当的数学方法进展求解，这便是电力系统暂态分析的一般方法。然而，由于各元件的动态响应有所不同，系统各种暂态过程的性质也不一样。因此，在不同目的暂态过程分析中，所考虑的元件种类与对它们数学模型的要求并不一样。例如，在电磁暂态过程分析中，所研究的暂态过程持续时间通常较短，在此情况下，一些动态响应比较缓慢的元件，如原动机及调速系统等的影响往往可以忽略不计，而发电机定子回路与电力网中的电磁暂态过程那么需加以考虑。相反，在电力系统稳定性分析中，那么通常忽略发电机定子回路与电力网中的电磁暂态过程，而将线路与变压器等元件用它们的等值阻抗来描述。另外，就同一种系统暂态过程来说，对于不同的分析精度与速度要求，元件所用数学模型的准确程度也不一样。一般地说，在进展系统规划与设计时，暂态分析的精度要求可以适当降低，这时各元件可以采用较粗略的数学模型，以便提高分析速度。因此，在建立元件数学模型时，不但需要研究它们的准确模型而且需要考虑各种简化模型，以适应不同的需要。在建立元件模型时还必须注意研究对象所处的频率范围。例如，在计算交流过电压时，变压器采用通常的以互感及绕组漏感与电阻表示的模型，如图1-5(a)所示；但在计算操作过电压时，除了上述要素外，还需要考虑绕组的对地电容与端子间电容及绕组间电容，如图1-5(b)所示；而在计算雷过电压时，变压器模型通常用冲击电容表示，无需考虑电感与电阻要素，如图1-5(c)所示。当研究现象的频率很高时，变压器与互感器的杂散电容、引线的微小电感对计算结果都有举足轻重的影响。图1-5 变压器的模型 行波现象与分布参数电力系统采用长线路将能源中心发出的电能输送给各电力用户，长线路的具体形式有架空输电线路与电缆线路两种，每微段的线路都呈现自感与对地电容，即线路是具有分布参数特性的电路元件。当电力系统中某一点突然发生雷电过电压或操作过电压时，这一变化并不能立即在系统其他各点出现，而是以电磁波的形式按一定的速度从电压或电流突变点向系统其他部位传播。例如，当架空输电线路遭受雷击时，雷击点导线将产生雷电过电压，该电压将沿着导线向两侧传播。这个沿线路传播的电压以及及其相伴而行的电流波称为行波。当行波到达变电站或其他节点时，由于电路参数的改变，将引起波的折射与反射。这种在分布参数电路中产生的暂态过程本质上是电磁波的传播过程，简称波过程。实际电力系统采用三相交流或双极直流输电，属于多导线线路，而且沿线路的电场、磁场与损耗情况也不尽一样，因此所谓的均匀无损单导线线路实际上是不存在的。但为了提醒线路波过程的物理本质与根本规律，可暂时不考虑线路的电阻与电导损耗，并假定沿线线路参数处处一样，即首先研究均匀无损单导线中的波过程。1. 波传播的物理概念假设有一无限长的均匀无损单导线，见图1-6(a)， =0时刻合闸直流电源，形成无限长直角波，单位长度线路的电容、电感分别为0、0，线路参数看成是由无数很小的长度单元构成，如图1-6(b)所示。图1-6 均匀无损的单导线a单根无损线首端合闸于E；b等效电路合闸后，电源向线路电容充电，在导线周围空间建立起电场，形成电压。靠近电源的电容立即充电，并向相邻的电容放电，由于线路电感的作用，较远处的电容要间隔一段时间才能充上一定数量的电荷，并向更远处的电容放电。这样电容依次充电，沿线路逐渐建立起电场，将电场能储存于线路对地电容中，也就是说电压波以一定的速度沿线路x方向传播。随着线路的充放电将有电流流过导线的电感，即在导线周围空间建立起磁场，因此与电压波相对应，还有电流波以同样的速度沿x方向流动。综上所述，电压波与电流波沿线路的传播过程实质上就是电磁波沿线路传播的过程，电压波与电流波是在线路中传播的伴随而行的统一体。2. 波动方程及解为了求出无损单导线线路行波的表达式，令x为线路首端到线路上任意一点的距离。线路每一单元长度dx具有电感0dx与电容0dx，如图1-7所示，线路上的电压与电流都是距离与时间的函数。图1-7 均匀无损单导线的单元等值电路根据节点电流方程可知，根据回路电压方程可知，整理得 1-6 1-7由方程(1-6)对再求导数，由方程(1-7)对再求导数，然后消去，并用类似的方法消去得。 1-8 1-9其中为单位长度电感与电容。通过拉普拉斯变换将(,)变换成(,)，(,)变换成(,)，并假定线路电压与电流初始条件为零，利用拉氏变换的时域导数性质，将式(1-8)、式(1-9)变换成 1-10 1-11其中。根据2阶齐次线性微分方程性质，令，那么式1-10、1-11解为 1-12 1-13将以上频域形式解变换到时域形式为 1-14 1-15式(1-14)、式(1-15)就是均匀无损单导线波动方程的解。3. 波速与波阻抗在波动方程中定义为波传播的速度。对于架空线路即沿架空线传播的电磁波波速等于空气中的光速度。而一般对于电缆，波速，其传播速度低于架空线，因此减小电缆介质的介电常数可提高电磁波在电缆中传播速度。定义波阻抗一般对单导线架空线而言，Z为500左右，考虑电晕影响时取400左右。由于分裂导线与电缆的较小而C0较大，故分裂导线架空线路与电缆的波阻抗都较小，电缆的波阻抗约为十几欧姆至几十欧姆不等。波阻抗Z表示了线路中同方向传播的电流涉及电压波的数值关系，但不同极性的行波向不同的方向传播，需要规定一个正方向。电压波的符号只取决于导线对地电容上相应电荷的符号，与运动方向无关。而电流波的符号不但及相应的电荷符号有关，而且及电荷运动方向有关，根据习惯规定：沿x正方向运动的正电荷相应的电流波为正方向。在规定行波电流正方向的前提下，前行涉及反行波总是同号，而反行电压涉及电流波总是异号，即必须指出，分布参数线路的波阻抗及集中参数电路的电阻虽然有一样的量纲，但物理意义上有着本质的不同：1波阻抗表示向同一方向传播的电压波与电流波之间比值的大小；电磁波通过波阻抗为Z的无损线路时，其能量以电磁能的形式储存于周围介质中，而不像通过电阻那样被消耗掉。2为了区别不同方向的行波，Z的前面应有正负号。3如果导线上有前行波，又有反行波，两波相遇时，总电压与总电流的比值不再等于波阻抗，即4波阻抗的数值Z只及导线单位长度的电感0与电容0有关，而及线路长度无关。4. 前行波与反行波下面用行波的概念来分析波动方程解的物理意义。首先讨论式(1-10)，电压u的第一个分量。设任意波形的电压波沿着线路x传播，如图1-8所示，假定当tt1时刻线路上任意位置点的电压值为ua，当时间t=t2时刻时(t2)，电压值为ua的点到达，那么应满足 即图1-8 行波运动由于v恒大于0，且由于(t2)，那么，由此可见表示前行波；同样的方法可以证明表示沿x反方向行进的电压波，称为反行电压波。，的证明过程类似。为方便将式(1-14)与式(1-15)可写成 1-16 1-17由式(1-16)与式(1-17)可知，线路中传播的任意波形的电压与电流传播的前行波与反方向传播的反行波，两个方向传播的波在线路中相遇时电压涉及电流波的值符合算术叠加定理，且前行电压涉及前行电流波的符号一样，反行电压涉及反行电流波的符号相反。5. 行波的折射与反射当波沿传输线传播，遇到线路参数发生突变，即波阻抗发生突变的节点时，都会在波阻抗发生突变的节点上产生折射与反射。如图1-9所示，当无穷长直角波沿线路1到达A点时后，在线路1上除、外又会产生新的行波、，因此线路上总的电压与电流为 1-18图1-9 波通过节点的折反射设线路2为无限长，或在线路2上未产生反射波前，线路2上只有前行波没有反行波，那么线路2上的电压与电流为 1-19然而节点A只能有一个电压电流，因此其左右两边的电压电流相等，即，因此有 1-20将，代入上式得， 1-21其中，分别为折射及反射系数。，计算如式1-22所示。 1-22另外，线路是用集中参数还是分布参数等值，除跟线路长度有关，还及暂态过程的频率有关。设线路长300m(约一档距) ，幅值为I的正弦波电流以光速(300m/s)传播，如图1-10所示，AB两点间的传播时间t为1s。如果是工频50Hz，两点的电流的差值最大为×10-4I ，这样可以看成同一值。但如果是100kHz，其差值最大可到达。因此在高频领域，即使距离很短，例如变电站的母线，也要考虑波的传播过程，即当成分布参数线路处理。变压器有时也要当作分布参数线路处理。图1-10 线路上的波过程 非线性元件与开关操作电力系统的暂态过程往往是因状态的变化而造成的。这种变化可以是断路器正常或故障操作而引起触头的闭合或开断；可以是雷电入侵波或操作过电压引起有间隙避雷器间隙击穿或电流过零时电弧的熄灭；也可以是系统发生故障造成相对地或相间突然短接等。在暂态计算中把电路中节点之间的闭合与开断用广义的开关操作来表示。因此，开关的计算模型以及正确处理开关操作所引起系统状态变化的程序方法，是电力系统电磁暂态计算的重要组成局部。电力系统中大局部元件属于线性元件，或可以近似地认为是线性元件，但也有一些元件具有明显的非线性特性，这些特性对暂态过程产生明显的影响。典型的非线性元件有避雷器的非线性电阻，如图1-11所示；变压器或电抗器等铁磁元件因铁心饱与而形成的非线性电感；以及断路器、保护间隙的电弧电阻等。因此，在暂态计算程序中应包括计及这些非线性元件特性的数学模型，并且含有一定的求解非线性电路的数值分析方法。对工程计算来说，还希望计算模型与分析方法尽可能实用，以便在尽可能短的计算时间里，得到具有一定准确度的结果。图1-11 避雷器的电压-电流特性在实际计算中，经常采用被称为分段线性化的方法来处理非线性元件，即把非线性元件的特性用几段具有不同斜率的直线线段来表示，把非线性元件局部等值为线性元件。 元件参数的频率特性在电力系统电磁暂态分析过程中，一个元件的特性模拟，不只是要作出正确的等值电路，还要模拟它的频率特性，因为这些频率特性有时对暂态现象有着决定性的影响。在暂态计算中，通常需要考虑频率特性的元件是架空线路与电缆。架空线路的正序电感1实际上是常数，在导线的趋肤效应不显著时，正序电阻根本上也是常数。零序电感与零序电阻那么因大地回路的趋肤效应而及频率密切相关。图1-12所示为架空线路电阻与电感的频率特性。变压器参数也有频率特性，但通常没有考虑。1图1-12 架空线路电阻与电感的频率特性 时间跨度的要求稳态计算的对象是一个时连续面，而暂态计算要模拟一个时间过程。数字计算机不可能连续地模拟暂态现象，只能在离散的时间点(步长t)求解，这将会导致累积误差。如何减少这类误差的积累是暂态仿真程序的重要课题。鉴于暂态计算的上述特点，暂态计算比稳态计算不管是程序编制还是应用难度都要大得多。1.5 电力系统数字仿真1.5.1 电力系统数字仿真的分类根据原型系统、数学模型与数字计算机三者的特征可以把电力系统数字仿真分成各种不同的类型。按照原型系统状态变化的时间过程，可分为连续系统仿真与离散事件系统仿真。连续系统仿真的系统状态量随时间连续变化，它的数学模型是一组方程式，包括连续模型(用微分方程描述)、离散时间模型(用差分方程描述)与连续及离散混合模型。离散事件系统仿真的系统状态量只在一些时间点上由某随机事件的驱动而发生变化，这类系统在两个事件之间其状态量保持不变，它的数学模型一般只用流程图或网络图描述。按照仿真目的，可分为以分析研究为目的的研究用系统仿真与以培训运行人员为目的的培训用系统仿真。研究用电力系统数字仿真，如电力系统电磁暂态计算程序EMTP，它可用于研究由开关操作、故障与雷击等引起的电磁暂态、电磁谐振与机电振荡，也可用于研究交直流换流器、控制系统与继电保护装置等的特性。除此以外，还有大量适合于专门功能的电力系统数字仿真程序，如电力系统综合程序BPA等。培训用电力系统数字仿真，如电力系统调度员培训仿真系统DTS、变电站培训仿真系统等，利用计算机及相关设备，将电力系统完整的模拟出来，并可以在上面进展正常操作训练及故障排除训练。培训用的仿真是为了训练系统运行、调度人员对系统环境的反响与判断能力，因此要求仿真的环境尽可能逼真，而对于仿真精度，只由培训的要求决定。研究用电力系统数字仿真又可分为系统稳态计算与暂态计算两大类。当需要研究电力系统处于相对平衡状态的运动特性时，采用系统稳态仿真；当研究系统处于受扰动状态的运动特性时，那么采用系统暂态仿真。两者的数学模型不同，仿真方法也不一样。稳态仿真中有潮流计算、故障计算，以及稳定计算与电压稳定计算中的静稳定计算。暂态仿真中有过电压计算、次同步振荡SSR计算、暂态恢复电压TRV计算、高次谐波计算，及稳定计算与电压稳定计算中的暂态稳定性计算。电力系统数字仿真按照研究电力系统运行状态的分类如图1-13所示。图1-13 电力系统数字仿真按照研究电力系统运行的状态的分类示意图按照计算方法，电力系统数字仿真又可以分成有效值计算与瞬时值计算两大类。有效值计算用于大规模系统的比较长时间的状态模拟，而瞬时值计算用于局部系统的短时间的暂态模拟。有效值计算通常只使用正序回路，而瞬时值计算需要使用三相回路。但是，在包含电力电子设备的情况，如高压直流系统，由于晶闸管等开关的频繁动作，即使是长时间模拟也需要采用瞬时值计算。另一方面，由于硬件设备与计算技术的进步，大规模系统的瞬时值在线(实时)计算也成为可能了。稳态计算与暂态计算中的动稳定计算属于有效值计算范畴，其它的暂态计算都属于瞬时值范畴。1.5.2 电力系统数字仿真的优点数字仿真在电力系统研究中得到快速的开展，除了计算机技术与软件技术的飞速开展外，电力系统数字仿真的独特优点是促使其快速开展的重要因素，这些优点有：1数字仿真不受被研究电力系统的规模与复杂性的限制。世界各国都在不断扩大电力系统的规模，大多数工业兴旺国家都建立了自己的全国统一电力系统，有些相邻国家间还建立了跨国联合电力系统。我国已实现跨区域电力系统的互联，依托三峡工程，实现以长江三峡为中心的全国统一电力系统。随着规模的增大，电力系统的构造也变得更加复杂。这些规模庞大与构造复杂的大系统，试验与研究的现场实测方法已很难进展。在电力系统动态模拟上做几十台发电机、几十条输电线路的电力系统暂态过程试验也是相当困难，而采用电力系统数字仿真就不存在这些困难，可以进展数百台发电机与上千条输电线路的大型电力系统数字仿真。2保证被研究系统的平安性。电力系统的故障试验、稳定性破坏试验、核电站控制试验等，直接在原型系统上做实验有很大的危险性，甚至是不允许的，这时，用电力系统数字仿真的方法是唯一可行的途径。3提高系统试验的经济性。在实际电力系统上做试验要暂停局部用户供电，需要配备各种测量设备、测量通道、通信工具，要求很多运行、调度人员与测试人员密切配合，花费大量人力、物力与财力，因此这种试验很难实现。如果用数字仿真做试验，所需费用要少得多。而且，数字仿真试验的设备一般都可重复使用，只需少数计算人员参加，试验时间很短。4增强对电力系统开展的预测性。需要对未来电力系统的特性做预测性的分析与研究，这些工作在实际电力系统中难以实现，而系统数字仿真可以对设计方案进展大量试验与计算，进展经济技术比较与优化，还可以对未来系统的假设条件的合理性进展验证。电力系统规划的方案是靠仿真得到的；新元件的接入、运行方式确实定是用仿真结果作为依据的；新方法研究、新装置设计、参数确定是用仿真来确认的。电力系统仿真软件试验已经成为电力系统设计、规划与运行阶段不可或缺的局部。1.5.3 电力系统数字仿真软件世界各国的电力系统数字仿真软件众多，目前国内外获得广泛应用的电力系统仿真软件主要有3大类：基于瞬时值计算的离线仿真软件，如EMTP、PSCAD/EMTDC、PSAPAC等；基于有效值计算的离线仿真软件，如BPA包括中国版BPA、PSASP、PSS/E、NETOMAC等；基于瞬时值计算的实时仿真软件，如RTDS、HYPERSIM。1. 电磁暂态分析程序EMTP电磁暂态分析程序EMTPElectro-Magnetic Transients Program是国际公认的电力系统电磁暂态分析的标准程序，其创始人是加拿大UBC大学任教的Dommel教授，目前EMTP有三个版本，即BPA-EMTP、ATP-EMTP与DCG-EMTP。BPA-EMTP是最早由美国邦纳维尔电力局Bonneville Power Adminstration，BPA无偿提供的EMTP版本，其用户现在大多已转用ATP-EMTP。ATP-EMTP是BPA的Scott-Meyer以自己的业余时间与资金开发的BPA-EMTP的替代程序，ATP-EMTP坚持无偿提供的原那么，在全世界拥有最多的用户，是目前国际上主流版本的EMTP程序。DCG-EMTP是1981年成立的DCGEMTP合作开发组织开发的EMTP版本，需有偿使用。EMTP具有分析功能多与元件模型全等优点，可以用于电力系统的稳态与暂态仿真分析，系统可由集中参数、分布参数元件、线性及非线性元件、具有频率相关参数的线路、各种类型开关、电力电子元件、变压器及电机、多种类型电源、控制电路的任意组合构成，只要是电路计算的范畴，对研究对象几乎没有限制。EMTP的计算精度经过了IEEE与CIGRE等国际权威组织的认定，因此计算结果的可信性很高。实际上，EMTP是开发其他的电磁暂态程序，如EMTDC、RTDS、HYPERSIM的根底。2直流电磁暂态程序PSCAD/EMTDCDennis Woodford博士于1976年在加拿大曼尼托巴直流研究中心开发完成了EMTDCElectro-Magnetic Transients in DC System的初版，它既可以研究交直流电力系统问题，又能完成电力电子仿真及其非线性控制，是一个离线仿真的电磁暂态计算程序，它有准确的直流元件模型、方便的数据输入方式以及强大的数据分析功能，是进展直流系统分析与工程研究的有力工具。事实上，EMTP程序所采用的电力系统模型与技术都可以应用于EMTDC中。PSCADPower System Computer Aided Design是其图形用户界面，PSCAD的开发成功，使得用户能更方便地使用EMTDC进展电力系统分析，使电力系统复杂局部可视化成为可能，而且软件可以作为实时数字仿真器的前置端。可模拟任意大小的交直流系统，在对直流系统电磁暂态仿真方面有绝对优势。操作环境为：UNIX OS, Windows95, 98,NT；Fortran 编辑器；浏览器与TCP/IP协议。主要功能如下： 1可以分析系统中断路器操作、故障及雷击时出现的过电压；2可对包含复杂非线性元件如直流输电设备的大型电力系统进展全三相的准确模拟，其输入、输出界面非常直观、方便；3进展电力系统时域或频域计算仿真；4电力系统谐波分析及电力电子领域的仿真计算；5实现高压直流输电、FACTS控制器的设计。3PSAPAC程序PSAPAC由美国EPRI开发，是一个全面分析电力系统静态与动态性能的软件工具。功能：DYNREDDynamic Reduction Program：网络化简及系统的动态等值，保存需要的节点。LOADSYNLoad Synthesis Program：模拟静态负荷模型与动态负荷模型。 IPFLOWInteractive Power Flow Program：采用快速分解法与牛顿拉夫逊法相结合的潮流分析方法，由电压稳态分析工具与不同负荷、事故及发电调度的潮流条件构成。TLIMTransfer Limit Program：快速计算电力潮流与各种负荷、事故及发电调度的输电线的传输极限。DIRECT：直接法稳定分析软件弥补了传统时域仿真工作量大、费时的缺陷，并且提供了计算稳定裕度的方法，增强了时域仿真的能力。LTSPLong Term Stability Program：LTSP是时域仿真程序，用来模拟大型电力系统受到扰动后的长期动态过程。为了保证仿真的准确性，提供了详细的模型与方法。VSTABVoltage Stability Program：该程序用来评价大型复杂电力系统的电压稳定性，给出接近于电压不稳定的信息与不稳定机理。为了估计电压不稳定状态，使用了一种增强的潮流程序，提供了一种接近不稳定的模式分析方法。ETMSPExtended Transient Midterm Stability Program： EPRI为分析大型电力系统暂态与中期稳定性而开发的一种时域仿真程序。为了满足大型电力系统的仿真，程序采用了稀疏技术，解网络方程时为得到最适宜的排 序采用了网络拓扑关系并采用了显式积分与隐式积分等数值积分法。SSSP(Smallsignal Stability Program)：该程序有助于局部电厂模式振荡与站间模式振荡的分析，由多区域小信号稳定程序(MASS)及大型系统特征值分析程序PEALS 两个子程序组成。MASS程序采用了QR变换法计算矩阵的所有特征值，由于系统的所有模式都计算，它对控制的设计与协调是理想的工具；PEALS使用了两 种技术：AESOPS算法与改良Arnoldi方法，这两种算法高效、可靠，而且在满足大型复杂电力系统的小信号稳定性分析的要求上互为补充。4BPA程序BPA程序是美国联邦政府能源部下属邦纳维尔电力局计算方法开发组自20世纪60年代初期开发的大型电力系统离线分析程序。该程序