6METALLURGICALPOWER冶金动力2014年第11期总第177期电力系统短路电流计算的算法设计与研究章燎（酒泉钢铁（集团）公司能源中心，甘肃嘉峪关735100）【摘要】根据酒钢集团公司发电厂和电网的实际需要，开发了一套短路电流计算软件。详细介绍了各序节点导纳矩阵的形成、考虑各种运行方式及短路信息下节点导纳矩阵的修正，应用连续回代法求取各序节点阻抗矩阵的方法，以及变结构电网的处理方法。最后给出了实际算例，证明了本论文提出方法的正确性。【关键词】短路电流计算；节点导纳矩阵；节点阻抗矩阵；变结构支路【中图分类号】TM76【文献标识码】A【文章编号】1006-6764(2014)11-0006-05DesignandStudyontheArithmeticofShortCircuitCurrentCalculationinPowerSystemsZHANGLiao(TheEnergyCenterofJiuquanIronandSteelCo.,Group,Jiayuguan,Gansu735100,China)【Abstract】TomeetthepracticalneedofJiuquanSteel’sself-supplypowerplant,acalculationsoftwareforshortcircuitcurrentwasdeveloped.Thispapermainlydemonstratesthecalculationofshortcircuitcurrentsystem,introducesindetailtheformationofnodeadmit-tancematrixofeachsequencewithconsiderationtothechangesunderdifferentoperationmodesandshortcircuitinformation,themethodofobtainingnodeimpedancematrixforeachsequencebyapplyingcontinuoussubstitutionandthetreatmentmethodsforvariablestructurepowernetworks.Finally,anexampleofcalculationwasprovidedtoshowthevalidityofthemethod.【Keywords】shortcircuitcalculation;nodeadmittancematrix;nodeimpedancematrix;variablestructurebranch1引言短路电流计算在电力系统地设计、设备选型和继电保护定值计算等方面应用广泛。随着国民经济的发展，电力系统结构越来越复杂，电力系统安全稳定问题也越来越重要。如何依据现代电力系统安全稳定运行的要求，开发一套适合实际工程需要的短路电流计算软件具有重要的现实意义。短路电流计算大致经历了三个阶段：人工计算阶段，早期程序计算阶段，新兴程序计算阶段。20世纪70年代及以前，各地电网结构还非常简单，短路电流计算通过人工实现。随着计算机技术的发展，在80年代末期，逐渐出现了基于DOS操作系统的短路电流计算软件。进入90年代以来，随着计算机硬件和软件技术的发展，涌现了大量基于Windows图形界面和面向对象技术的故障计算软件。故障计算的算法，涵盖了各种等级和形式的电网及简单/复杂故障，节点优化编号技术和稀疏技术被广泛应用，在这种情况下，许多电力科学工作者相继提出了一些适合配电网短路电流计算的方法。2短路电流计算软件的设计2.1软件设计目标及开发和运行环境根据酒泉钢铁（集团）公司自备发电厂和电网的实际需要，开发短路电流计算软件的具体目标为：（1）可视化，全图形操作；（2）图形与数据库结合；（3）通用性：要适于各类不同的电网；（4）实用性：结合实际电网，计及互感线路对零序电流的影响，且能计算任意复杂故障；（5）实效性：应尽可能简化程序，提高计算速率，满足实际工程的需要。软件采用C#作为开发工具，在WindowXp操作系统下运行。C#具有面向对象语言所应有的基本特性：封装、继承和多态性。同时，它与Web紧密结合，可以用统一的方法来处理本地的和远程的C#对2014年第11期总第177期METALLURGICALPOWER冶金动力7象，且C#组件能够方便地转为Web服务，并被其他平台上的各种编程语言调用。数据库采用SQL数据库。2.2软件总体设计依照系统要实现的目标，系统整体框架如图1所示。图1系统总框图整个系统按功能划分成四个模块，即数据库系统的建立、网络拓扑图处理、短路电流计算和文档管理。每一功能模块可继续进行划分，直到最小的单元模块为止，这样可保证系统具有良好的可扩充性和通用性。其中数据库系统位于服务器端，用于系统数据的建模，存储和管理员级的维护，并支持数据的多用户和网络共享，其他3个模块集成为客户端应用程序，网络拓扑图处理模块用于处理交互式绘制电网接线图；短路电流计算模块面向一个故障计算过程储方式。具体方法如下：完成核心任务；文档管理模块将以上各模块有机结合在一起，响应用户请求的总调度和管理功能。客户零元素的导纳值；端应用程序通过局域网向数据库服务器请求数据并返回计算处理结果。其中网络拓扑图处理目前已是元素的行标；一个独立的研究领域，应用较广。短路电流计算软件总体结构见图2。元素的列标；短路电流计算共包括如下几部分：节点编号优化、各序节点导纳矩阵的形成、按运行方式和短路信息修正节点导纳矩阵、按保留节点和非保留节点对中的位置，各序节点导纳阵进行行列变换，并用高斯消元法消去非保留节点，简化网络、连续回代法求各序节点阻个非零元素在抗阵、变结构支路的处理、最后根据不同的短路情况3.1.1计算对应的短路电流与短路电压。其中，节点编号优化、正负序、零序节点导纳矩相同，因此正序节点导纳阵的求法与负序节点导纳阵的形成三个子模块属于原始数据处理部分，为后阵也完全相同。本文使用支路扫描法形成正序和负面的计算提供原始电网数据。用户选择运行方式和序节点导纳阵。经过拓扑分析后的所有支路只有两短路信息后，各序节点导纳矩阵要做相应修正。个端点或一个端点。对所有支路扫描一次就可以形3节点导纳矩阵的形成及修改成节点导纳阵。而不必形成支路阻抗阵，支路导纳3.1节点导纳矩阵的形成阵，图2短路电流计算软件总体结构图为了提高计算速度，导纳阵的存储采用链表存Posi_Value———按行存储正序节点导纳阵中非Posi_Row———按行存储正序节点导纳阵中非零Posi_Col———按行存储正序节点导纳阵中非零Posi_NA———每行非零元素的个数；Posi_LINK———下一个非零元素在Posi_Value对每行最后一个非零元素，该值置为0；Posi_IA———存储正序节点导纳阵中每行第一Posi_Value中的位置。正负序节点导纳矩阵的形成正序支路阻抗阵的结构与负序支路阻抗阵完全可大大提高计算速度。也不必形成关联矩阵，8METALLURGICALPOWER冶金动力2014年第11期总第177期3.1.2零序节点导纳矩阵的形成当发生接地短路时，零序网中会出现一些空节点，即不连接任何支路的节点。编写程序时，应增加判别空节点的功能，避免运算时出错。为了确保发生接地短路时能快速求取零序电流，本软件采用正常运行时求初始零序节点导纳阵，接地短路时对该节点导纳阵进行修正的处理方式。本节主要介绍初始零序节点导纳阵的形成方法。（1）变压器支路的处理在读取两卷变压器和三卷变压器的时候，要考虑是否是零序终止支路。如两卷变压器为YN，d接线方式时，角接端节点即为零序中止节点。对于此种接线方式，驻端节点的自导纳用一非常小的实数代替，互导纳为0；Y端节点的自导纳为该变压器支路导纳，与驻端节点间的互导纳为0。读取三卷变压器的数据时，处理方法类似。（2）输电线路的处理为了处理方便，对于输电线路的处理，总体上分为考虑互感和不考虑互感两种情况，而开始均按照无互感处理。对于考虑互感的情况，只需用互感信息对已形成的节点导矩阵进行修正即可。（3）其他支路的处理对于其他支路，如发电机和外系统支路，当它在零序网中时，与在正序节点导纳阵中的处理方法类似；当该支路确定不在零序网中时，如中性点不接地的发电机支路，将该节点的自导纳用一非常小的实数代替，互导纳为0。当扫描完所有支路后，应对已存储的零序节点导纳阵中非零元素进行重新排序。重新排序的原则与正序节点导纳阵相同。当发生接地短路时，只需从短路点开始搜索对应的零序网并对最初形成的零序节点导纳阵进行修正。此外，由于本套短路电流计算软件最终是为满足酒泉电厂继电保护整定需要，为使计算出的短路电流最大，不计负荷影响。3.2节点导纳矩阵的修改实际电网运行中，需要考虑新建线路、检修线路以及故障等情况下电网运行状态的调整，故需要按运行方式的变化修改节点导纳阵。由于实际电网中，运行方式的改变有多种形式，为了方便表示，在运行方式变化表中设置了一标志位flag，表示不同的运行信息。flag＝1~10分别代表十种情况：（1）追加非互感树支，且为非接地支路；（2）追加非互感接地支路；（3）追加非互感连支；（4）切除非接地互感树支或非互感连支；（5）切除非互感接地支路；（6）追加互感树支；（7）追加互感连支；（8）完全切除互感树支；（9）单端切除互感连支；（10）停电检修互感连支。不同的操作信息对应各自的节点导纳矩阵的修正方式，以正序网flag=2为例，其具体处理方法如下：当追加的非互感树支为双端支路，即flag=1时，若其非零端点的节点号为i、j，数组Posi_Value、Posi_Row、Posi_Col的长度t均增加3。令t=t+3，则相关数组处理如下：Posi_Value(t)=yl；Posi_Row(t)＝j；Posi_Col(t)=j；Posi_Value(t－1)=－yl；Posi_Row(t－1)＝j；Posi_Col(t－1)=i；Posi_Value(t－2)=－yl；Posi_Row(t－2)＝i；Posi_Col(t－2)=j；Posi_Value(i，i)=Posi_Value(i，i)+yl；Posi_NA(i)=Posi_NA(i)+1；Posi_NA(j)=Posi_NA(j)+2。当扫描完运行方式变化表后，对此过程中涉及到的所有节点，应将各序节点导纳阵中这些节点对应的非零元素按照的行、列坐标的大小重新排序。重新排序的原则与前文相同。3.3短路信息的处理短路位置性质的不同会对各序节点导纳矩阵产生影响。当为母线短路时，节点导纳矩阵不必修正；当为支路短路时，需在故障位置增设节点。无互感支路发生短路和零序网的处理方法较为简单，以下以互感支路的处理为例讨论支路短路时节点导纳阵的修正方法。对于互感线路，由于它与其他线路间存在耦合关系，当发生支路故障时，新增节点与属于该互感组的其他线路的节点之间存在耦合导纳。所以，处理互感线路故障时，既需要象无互感线路一样将故障线路的零序阻抗按比例分解到故障新增节点两侧的线路上去，也需要将互感参数按比例分解，并且要保持正确的互感范围关系，以保证对多重复杂故障的正确处理。考虑任意重互感组和任意复杂部分互感结构的情况，设互感线路的重数为k，线路L（两端节点为i的bi%~ei%段与线路L（两端节点为的Pi、Qi）Pj、Qj）j100]，bj%~ej%段存在互感zmij，其中，bi，ei，bj，ej∈[0，设线路Li上的fi%位置发生短ei>bi，ej>bj，i，j∈[1，k]。路故障，新增节点为ri，则zmij的分解规则为：（1）如果bi
