含大量继保和安自装置的连锁故障并行仿真软件研发.pdf

第４２卷第２ｌ期２０１４年１１月１日电力系统保护与控制ＰｏｗｅｒＳｙｓｔｅｍＰｒｏｔｅｃｔｉｏｎａｎｄＣｏｎｔｒｏｌＶｌ０１．４２ＮＯ．２１ＮＯＶ．１，２０１４含大量继保和安自装置的连锁故障并行仿真软件研发徐得超１，刘巍，朱旭凯，夏天，朱全胜，徐树文，陈德扬（１．中国电力科学研究院，北京１００１９２；２．国网河南省电力公司电力科学研究院，河南郑州４５００５２）摘要：基于ＭＰＩ实现了一种基于用户自定义建模思想的含有大量继电保护和安全自动装置的连锁故障并行仿真软件。软件首先采用面向对象的用户自定义建模方法，完成大量继电保护和安全自动装置建模；然后采用分网并行和分组并行计算技术，将大规模电网进行网络分割，分配给子网进程；将大量继电保护和安全自动装置模型进行分组划分，分配给外接进程。任务分配后的数据被提交到并行集群上，执行连锁故障的并行仿真。由于子网进程和外接进程由不同的ＣＰＵ并行计算，程序效率得以提高，能在用户可接受的时间内完成大电网的连锁故障仿真。关键词：用户自定义建模；连锁故障；并行计算；安全自动装置；继电保护；面向对象设计Ｒｅｓｅａｒｃｈａｎｄｄｅｖｅｌｏｐｍｅｎｔｏｆｃａｓｃａｄｉｎｇｆａｉｌｕｒｅｓｉｍｕｌａｔｉｎｇｓｏｆｔｗａｒｅｃｏｎｔａｉｎｉｎｇｍａｎｙｒｅｌａｙｉｎｇｐｒｏｔｅｃｔｉｏｎｓａｎｄａｕｔｏｍａｔｉｃｓａｆｅｔｙｄｅｖｉｃｅｓ—ｘｕＤｅｃｈａｏ，ＬＩＵＷｅｉ，ＺＨＵＸｕ．ｋａｉ，ＸＩＡＴｉａｎ，ＺＨＵＱｕａｎ．ｓｈｅｎｇ，ＸＵＳｈｕ．ｗｅｎ，ＣＨＥＮＤｅ．ｙａｎｇ（１．ＣｈｉｎａＥｌｅｃｔｒｉｃＰｏｗｅｒＲｅｓｅａｒｃｈＩｎｓｔｉｔｕｔｅ，Ｂｅｉｊｉｎｇ１００１９２，Ｃｈｉｎａ；２．ＳｔａｔｅＧｒｉｄＨｅｎａｎＥｌｅｃ￣ｉｃＰｏｗｅｒＲｅｓｅａｒｃｈＩｎｓｔｉｔｕｔｅ，Ｚｈｅｎｇｚｈｏｕ４５００５２，Ｃｈｉｎａ）Ａｂｓｔｒａｃｔ：Ａｃａｓｃａｄｉｎｇｆａｉｌｕｒｅｐａｒａｌｌｅｌｓｉｍｕｌｍｉｎｇｓｏｆｔｗａｒｅｗｈｉｃｈｃａｎｓｉｍｕｌａｔｅｗｏｒｋｉｎｇｏｆｍａｎｙｒｅｌａｙｐｒｏｔｅｃｔｉｏｎｓａｎｄａｕｔｏｍａｔｉｃ—ｓａｆｅｔｙｄｅｖｉｃｅｓｉｎｐｏｗｅｒｓｙｓｔｅｍｉｓｉｍｐｌｅｍｅｎｔｅｄｂａｓｅｄｏｎＭＰＩａｎｄｕｓｅｒｄｅｆｉｎｅｄｍｏｄｅｌｌｉｎｇｔｅｃｈｎｉｑｕｅ．Ｆｉｒｓｔ，ｂｙｄｅｓｉｇｎｉｎｇｕｓｅｒ－ｄｅｆｉｎｅｄ—ｍｏｄｅｌｉｎｇａｐｐｒｏａｃｈｕｓｉｎｇｏｂｊｅｃｔｏｒｉｅｎｔｅｄｍｅｔｈｏｄ，ｍｏｄｅｌｉｎｇｏｆｍａｎｙｒｅｌａｙｐｒｏｔｅｃｔｉｏｎｓａｎｄａｕｔｏｍａｔｉｃｓａｆｅｔｙｄｅｖｉｃｅｓｉｓｃｏｍｐｌｅｔｅｄ．Ｔｈｅｎｕｓｉｎｇｐａｒａｌｌｅｌｃｏｍｐｕｔｉｎｇｔｅｃｈｎｉｑｕｅｓｏｆｓｐｌｉａｉｎｇｐｏｗｅｒｇｒｉｄａｎｄｇｒｏｕｐｉｎｇｍｏｄｅｌｓ，ｔｈｅｌａｒｇｅ－ｓｃａｌｅｐｏｗｅｒｇｒｉｄｎｅｔｗｏｒｋｉｓｐａｒｔｉｔｉｏｎｅｄｉｎｔｏｓｅｖｅｒａｌｓｕｂｎｅｔｓａｓｓｉｇｎｅｄｔｏｔｈｅｓｅｐｒｏｃｅｓｓｅｓｎａｍｅｄｓｕｂｎｅｔ－ｐｒｏｃｅｓｓ；ａｎｄｔｈｅｇｒｏｕｐｅｄｍｏｄｅｌｓｆｏｒｒｅｌａｙｐｒｏｔｅｃｔｉｏｎｓａｎｄａｕｔｏｍａｔｉｃｓａｆｅｔｙｄｅｖｉｃｅｓａｒｅａｓｓｉｇｎｅｄｔｏｔｈｅｓｅｐｒｏｃｅｓｓｎａｍｅｄｅｘｔｅｒｎａｌｐｒｏｃｅｓｓｂｙｇｒｏｕｐ．Ｔｈｅｄａｔａａｆｔｅｒｓｐｌｉｔｔｉｎｇｇｒｉｄａｎｄｇｒｏｕｐｉｎｇｍｏｄｅｌｓａｒｅｓｕｂｍｉｔｔｅｄｔｏｔｈｅｐａｒａｌｌｅｌｃｌｕｓｔｅｒａｎｄｔｈｅｐａｒａｌｌｅｌｓｉｍｕｌａｔｉｎｇｆ—ｏｒｃａｓｃａｄｉｎｇｆａｉｌｕｒｅｉｓｓｔａｒｔｉｎｇ．ＢｅｃａｕｓｅｔｈｅｓｕｂｎｅｔｐｒｏｃｅｓｓａｎｄｅｘｔｅｒｎａｌｐｒｏｃｅｓｓａｒｅｅｘｅｃｕｔｅｄｉｎｐａｒａｌｌｅｌｂｙｄｉｆｆｅｒｅｎｔＣＰＵｓ，ｓｉｍｕｌａｔｉｎｇｅｆｆｉｃｉｅｎｃｙｃａｎｂｅｉｍｐｒｏｖｅｄ，ａｎｄｃａｓｃａｄｉｎｇｆａｉｌｕｒｅｓｉｍｕｌａｔｉｎｇＣａｌｌｂｅｃｏｍｐｌｅｔｅｄｗｉｔｈｉｎａｎａｃｃｅｐｔａｂｌｅｔｉｍｅｏｆｕｓｅｒｓ．ＴｈｉｓｗｏｒｋｉｓｓｕｐｐｏｒｔｅｄｂｙＮａｔｉｏｎａｌＮａｔｕｒａｌＳｃｉｅｎｃｅＦｏｕｎｄａｔｉｏｎｏｆＣｈｉｎａ（８６３ｐｒｏｇｒａｍ）（Ｎｏ．２０１１ＡＡ０５Ａ１０５）．Ｋｅｙｗｏｒｄｓ：ｕｓｅｒ－ｄｅｆｉｎｅｄｍｏｄｅｌｉｎｇ；ｃａｓｃａｄｉｎｇｆａｉｌｕｒｅ；ｐａｒａｌｌｅｌｃｏｍｐｕｔｉｎｇ；ｓｅｃｕｒｉｔｙａｕｔｏｍａｔｉｃｅｑｕｉｐｍｅｎｔ；ｒｅｌａｙｐｒｏｔｅｃｔｉｏｎ；ｏｂｊｅｃｔ－ｏｒｉｅｎｔｅｄｄｅｓｉｇｎｉｎｇ中图分类号：ＴＭ７１文献标识码：Ａ—文章编号：１６７４－３４１５（２０１４）２１－０１３２００７０引言”正常情况下的电力系统是一个动态平衡的系。但如果电网发生故障导致电压电流发生大幅度波动并越限，电网内继电保护及安全自动装置将相继动作，可能导致系统状态的进一步恶化，引起连锁故障。机电暂态仿真主要研究电力系统受到大扰动基金项目：国家高技术研究发展技术（８６３计划）课题（２０１１ＡＡ０５Ａ１０５）；国家电网公司科技项目（ＸＴ７１－１２－０１１）后的暂态稳定性，可用于分析受到诸如短路故障，切除线路、发电机、负荷，发电机失去励磁或者冲击性负荷等大扰动作用下，电力系统的动态行为和保持同步稳定运行的能力，是进行连锁故障仿真的手段之－－［－。通过机电暂态仿真研究大电网连锁故障，主要过程是建立符合实际的电网模型和元件模型，计算分析连锁故障的演变过程及其后果，并给出预防连锁故障扩展蔓延的措施。其难点之一在于电网中大量自动装置的建模，主要包括继电保护装置和安全徐得超，等含大量继保和安自装置的连锁故障并行仿真软件研发．１３３．自动装置【５】。继电保护装置种类多，既有主保护又有后备保护，原理复杂，而且随着电网发展，新型的继电保护装置层出不穷。安全自动装置涉及的控制结构更为复杂、在电网中经常跨区域控制。保护种类多和原理的纷繁复杂给继电保护和安自装置建模带来了软件编码和模型维护的难度和工作量，是大电网连锁故障仿真的挑战之一【７－１４］。连锁故障仿真的另一个难点在于，在大规模电网中，继电保护装置和安全自动装置数量多，规模庞大，导致计算量激增。以继电保护装置为例，如果采用用户自定义方式建模，一个继保装置由几十到上百个基本功能框组成，一个万节点的电力系统安装的继电保护装置将有几万个。如果采用交替求解法进行机电暂态仿真，一个时步内还需要进行多次迭代。如此下来，计算量将十分惊人。针对机电暂态仿真研究大电网连锁故障在建模和模型维护方面的难点，本文设计面向对象的用户自定义建模方法，充分利用面向对象的封装、继承、多态特性解决建模过程中遇到的可扩展性和可维护性差的问题。针对连锁故障仿真过程中计算量激增，本文结合分网并行和分组并行技术，利用子网进程￣ＨＪ，ｔ－接进程，构建连锁故障并行仿真系统，共同承担电网和模型的计算任务。应用表明，本文构建的并行仿真系统能有效减少仿真计算时间，达到用户可接受的时间范围内。本文的主要贡献在于将原有研发的ｊｊ取立功能集成，构建连锁故障并行仿真系统，使得大电网背景下大量继电保护和安自装置建模和连锁故障仿真实用化。１连锁故障并行仿真系统流程本文连锁故障并行仿真系统仿真流程如图１所示。图中进程分为两类：一类是子网进程，承担分网并行计算任务，负责联络系统或者某一个子网的连锁故障仿真。对于输入量和输出量位于同一子网的继保和自动装置模型，也由该子网进程负责计算，以减少进程间通信量。在子网进程中，进程号为０的进程称为主控进程，负责数据与模型分发、联络系统计算以及本子网的仿真。其他子网进程负责各自子网与模型的仿真。子网进程对应着图１中中间栏与右侧栏。另一类是外接进程，承担分组并行计算任务，负责用户指定的继电保护和安自装置模型（以下简称模型）以及输入输出量跨子网边界的模型仿真计算。根据用户的指定情况，可以有多个外接进程，每个外接进程负责～组模型的计算。外接进程对应图１中左侧栏。连锁故障并行仿真过程共分８个大步骤，分别为装置建模、网络分割与模型分发、通信接口初始化、输入量计算、模型初始化计算与输出量处理、电网与模型仿真和结果输出。每个大步中又分为若干小步骤，如在第２步网络分割与模型分发中，主控进程０读入全网数据，网络分割方案，模型拓扑②结构描述以及外接模型配置文件等信息（如），然③后确定设备和模型所属子网进程号并分发（如），④非主控子网进程则接收本子网计算数据（如）以⑤及本进程负责计算的模型信息（如），外接进程则⑥读入模型拓扑结构描述信息（如）。第３步通信接口初始化则分析外接模型的输入输出量所关联的设备信息，确定它们所属的子，然后发送给负责计⑧算该模型的外接进程（如），由外接进程接收处理⑦⑨⑩并排序后发送给对应子网进程（如图，，，⑩）。在第４步输入计算中，子网进程除了计算本子网仿真模型的输入量外，还需要计算本子网内部外接模型的输入量，然后发送给对应外接进程。：０厂丽：广；一：外接进程枷个）’于进程０（末摔仪１十）。于网进程个）图１连锁故障并行仿真流程Ｆｉｇ。１Ｆｌｏｗｃｈａｒｔｏｆｃａｓｃａｄｉｎｇｆａｉｌｕｒｅｐａｒａｌｌｅｌｓｉｍｕｌａｔｉｎｇ２装置建模方法本文采用用户自定义建模来完成大规模继电保护装置和安全自动装置的建模，并设计统一的仿真方法进行模型仿真。这样既解决了新增装置的建模仿真问题，又能保证建模过程的可扩展性和良好的维护性。装置建模过程分以下三步。１）定义基本功能框。基本功能框是用户自定义电力系统保护与控制建模的具有基本运算功能的最小功能单元，包含微分积分框、代数运算框、逻辑函数框、线性非线性控制等。２）定义基本功能框之间的拓扑。为了表示基本功能框之间的拓扑，定义基本功能框之间的连接公共点为端子。基本功能框之间的端子既表示前一个功能框的输出量，又代表后一个功能框的输入量。基本功能框和端子的关系可用图２（ａ１表示。图中，功能框１，２，３之间有两个端子，分别是１和２。２号端子是１号框的输出，同时是２，３号框的输入。３）定义基本功能框构成图的不相交割集为组合功能框。组合功能框可以看成一个抽象的基本功能框。和基本功能框一样，组合功能框具有输入量和输出量。其定义过程如图２（ｂ）和图２（ｃ）。（ｃ）图２端子与基本功能框及其组合功能框Ｆｉｇ．２Ｄｉａｇｒａｍｏｆｔｅｒｍｉｎａｌｓ，ｂａｓｉｃｆｕｎｃｔｉｏｎａｎｄｔｈｅｃｏｍｂｉｎａｔｉｏｎｏｆｔｈｅｍ图２（ｂ）中２、３、１０、１１四个基本功能框组成图中一个割集，形成组合功能框Ａ，８、９两个基本功能框组成图中与Ａ不相交的一个割集，形成组合功能框Ｂ。割集操作后形成新的用户自定义模型图如图２（ｃ）所示。图２（ｂ）和图２（ｃ）在功能上是等价的。提出组合功能框的意义在于：一些基本功能框可以组合成常用的组合功能框，并加入用户自定义模型库，从而便于模型的复用与扩展，方便用户自定义建模。结合上述用户自定义建模过程，本文采用面向对象的软件开发方法实现装置的建模过程，如图３所示。图３（ａ）描述了参数类、端子类和基本功能框类的关系。其中参数类描述模型中如时间常数、定值等属性。端子类描述模型中各功能模块之间的拓扑。基本功能框可以看作是装置中的各种功能元件的抽象，包含有与其他功能框的拓扑描述以及功能元件的一些公共的功能，如初始化、参数设定、计算等，这些功能被设计成虚化的函数，所有其他功能框均可继承该基本功能框的属性与功能，以实现面向对象的多态。如图３（ｂ）中加法框类、积分框类、限值框类等。图３（ｃ）中引入模型类。模型类从基本功能框类继承，具备基本功能框的一切属性，兼具有管理功能框、端子以及参数等的能力。（ａ）功能框基类依赖关系图（ｂ）各功能框继承关系图ｆｃ）模型类依赖和继承关系图图３面向对象的继保与安自装置建模Ｆｉｇ．３Ｏｂｊｅｃｔ－ｏｒｉｅｎｔｅｄｒｅｌａｙｐｒｏｔｅｃｔｉｏｎａｎｄｓａｆｅｔｙａｕｔｏｍａｔｉｃｄｅｖｉｃｅｍｏｄｅｌｉｎｇ３并行仿真方法本文连锁故障并行仿真方法由分网并行方法和分组并行方法，以及它们之间的接口组成。分网并行方法求解电网主系统，包括电网方程组和动态元件有关的微分方程组，由子网进程负责计算。分组并行方法求解用户指定的用户自定义模型，由外接进程负责。当有用户白定义模型输入输出量跨子网时，软件能自动指定其由外接进程负责。３．１分网并行方法Ｌ１＂】电力系统机电暂态仿真中网络方程的求解可归结为大型稀疏线性方程组的求解问题：Ａｘ＝其中：Ａ为对称稀疏阵，即导纳矩阵：为待求的未知数向量，即电压修正向量；６为已知的右端项向量。稀疏线性方程组的通用解法是三角分解法，解向量可由前代回代步骤得出。本文采用一种适于分网并行的方法一端口逆矩西匪一徐得超，等含大量继保和安自装置的连锁故障并行仿真软件研发阵法［１５－１６］来求解线性方程组。该方法的基本原理是根据阵逆阵在边界点上的对应元素组成的边界点系数阵，以及矩阵去掉边界点元素后求得的相应边界点初步解，求取边界点的准确解。然后求出各子网边界点右端项修正量，从而实现各子网的并行求解。该方法能够保持原有串行程序计算过程完整性，减少并行化开发工作量，提高软件可靠性。３．２分组仿真方法继保和安自装置等电网二次控制系统在连锁故障仿真过程中由于数量大而占用大量计算时间。由于网络分割，部分继保和安自装置的模型会随着分网而被分配到其关联设备子网中，计算自然并行化。但是对于跨子网的模型，以及那些分配过多模型的子网中模型，可以设置多个分组，设计分网并行和分组并行之间的接口，在多个外接进程中并行仿真，实现分组并行计算。分组并行计算仿真算法包括模型的初始化计算和仿真两部分，二者均基于用户白定义模型的传递函数框图，关键在于确定每个基本功能框的计算顺序，本文采用了文献［１２】中算法。３．３分网并行和分组并行接口本文连锁故障并行仿真系统中，分网并行计算由子网进程负责，而分组并行计算主要由外接进程负责。它们之间的关系是多对多的关系，也就是说，一个子网进程可能与多个外接进程通信，一个外接进程也可能与多个子网进程通信。图４描述了子网进程和外接进程之间并行仿真接口关系，虚线左侧是子网进程的计算过程，右侧是外接进程的计算过程。图４中代表仿真时步的计数，ｋ代表一个步长内网络方程组与微分方程组交替迭代的计数，初始化为０。左侧框１到９是一个时步内子网进程计算过程，框１０到１２是外接进程计算过程。当ｋ＝０时，需要保存上一时步电压量和状态量仿真结果，如框１和框１１所示。框２中子网进程根据电网元件存储位置和外接进程输入变量类型，按照第ｎ时步电压量和状态量求取外接进程输入变量值，由ＭＰＩ通信接口发送到外接进程。外接进程接收到输入量后，采用３．２节用户自定义仿真算法计算输出量。与此同时，框３可以并行进行微分方程组数值积分，求网络方程组右端项注入电流。框４通过ＭＰＩ接口从外接进程取得外接模型的输出量，并修正网络方程组右端项注入电流。框５求解网络方程组。框７进行迭代收敛性检查，如果收敛，则结束：步长计算，否则进行下一次迭代，转到框２。⑨图４子网进程和外接进程并行仿真方法Ｆｉｇ．４Ｐａｒａｌｌｅｌｓｉｍｕｌａｔｉｎｇｍｅｔｈｏｄｂｅｔｗｅｅｎｓｕｂｎｅｔｐｒｏｃｅｓｓａｎｄｅｘｔｅｍｐｒｏｃｅｓｓ图４中，框２到框１０，以及框ｌ２到框４是子网进程和外接进程的ＭＰＩ通信接口部分，它们是一对多的关系。为提高仿真效率，本文在仿真初始化阶段，根据外接模型的进程分配、输入输出量相关子网进程分布情况，对ＭＰＩ通信接口部分进行分析，确定数据交换的顺序与结构。具体来讲，在子网进程侧，按照外接进程排序，分别形成输入量的元件位置、变量类型信息数组，并申请对应顺序的变量值缓存；而在外接进程侧，则按照外接模型的排序，同样分别形成输入输出量元件位置、变量类型信息、和变量值缓存数组，并形成这两个缓存数组的位置映射关系。这样在仿真计算阶段，可以根据形成的顺序和索引快速检索，大大提高程序整体效率。４仿真算例本文采用河南电网２０１３年冬季大方式为例，验证本文实现方法的计算效率与效果。计算机环境为Ｉｎｔｅｌ至强处理器Ｘ５６５０，主频２．６７ＧＨｚ共１台，内存１６Ｇ，缓存１２Ｍ，同时支持１６个进程并行计算。电网和采用用户白定义模型的继保和安自装置仿真规模如表１和表２所示。４．１连锁故障仿真效果测试以河南电网周口地区为例测试连锁故障仿真效果，其与河南主网接线示意图如图５所示。测试案电力系统保护与控制例运行方式安排如下：周口地区５００ｋＶ电压等级层面通过邵陵．周口５００ｋＶ双回线与漯河地区相连，其中邵周Ｉ回因检修停运；２２０ｋＶ电压等级层面通过邵陵至川汇、邵陵至淮阳、薛坡至桐丘三回线与河南主网相连。此运行方式下，周口地区整体供电能力不足，其发电出力仅包括隆达电厂单台１３５ＭＷ机组，地区负荷主要由与河南主网相连的联络线承担。表１电网测试规模有效厂站４２８注：有效厂站指含厂站主接线的厂站，虚拟厂站未计算在内。表２用户自定义模型测试规模Ｔａｂｌｅ２Ｔｅｓｔｉｎｇｓｃａｌｅｏｆｕｓｅｒ－ｄｅｆｉｎｅｄｍｏｄｅｌｓ图５周口地区电网与河南主网接线示意图Ｆｉｇ．５ＣｏｎｎｅｃｔｉｏｎｄｉａｇｒａｍｏｆＺｈｏｕｋｏｕｇｒｉｄａｎｄＨｅｎａｎｍａｉｎｇｒｉｄ周口地区继电保护和安自装置的白定义模型主要安装在９个变电站内，其中５００ｋＶ站１个，２２０ｋＶ站８个，涉及线路１８条。自定义模型均依据实际的继电保护和安自装置建模，具体涉及到许继电气的ＷＢＨ．８０１变压器保护，ＷＸＨ．８０３和ＷＸＨ．８０２线路保护装置，国电南自ＰＳＬ．６０２Ｇ和ＰＳＬ．６０３Ｇ线路保护装置，国电南瑞ＲＣＳ一９３１ＡＭ，ＲＣＳ．９２５Ａ，ＲＣＳ．９３１Ｂ线路保护装置和ＲＣＳ．９９０Ａ稳控装置，以及四方公司的ＷＧＱ．８７１线路保护装置等。周口地区ＲＣＳ．９９０Ａ稳控系统由５００ｌ（Ｖ周口变、２２０ｋＶ川汇变、淮阳变、水寨变４个厂站组成。其中，周口变为稳控主站，川汇变、淮阳变为稳控子站，水寨变为远切负荷执行站。稳控系统主要功能如下：监视５００ｋＶ邵周双线的运行状态，接收川汇变发送的邵川线电流、功率；接收淮阳变发送的邵淮线、薛淮线电流、功率。当系统判断出５００ｋＶ邵周双线跳闸（检测２条线路功率代数和的绝对值“”小于低功率门槛定值）时，若同时检测到邵川线、邵淮线、薛淮线中的任一条发生过载时，装置进入过载动作逻辑，经延时后分轮次向切负荷执行站（川汇、淮阳、水寨）发远方切负荷命令。本测试案例中连锁故障的触发故障为５００ｋＶ邵周ＩＩ线首端三相短路故障，引发继电保护装置动作，切除故障线路邵周ＩＩ线。在此故障情况下，连锁故障软件在邵川线、邵淮线和薛桐线处观测到的仿真结果（功率，电压和电流）如图６所示。由图６中曲线可以看出，当５００ｋＶ邵周ＩＩ线因故障断开隋况下，原５００ｋＶ线路承担功率转移至三个２２０ｋＶ断面，其中邵淮线、邵川线电流过载越限，超出周口安稳装置整定值，稳控装置动作，７．２８Ｓ第一轮切负荷（川汇变）后，邵川线过载解除，１０．７１ｓ第二轮切负荷（淮阳变）后，邵淮线过载解除。从算例仿真结果可以看出，周口地区继保和安自装置动作与实际情况吻合，说明了连锁故障仿真软件的正确性。４．２连锁故障仿真效率测试采用４．１节所述的算例，测试连锁故障仿真软件的效率，测试规模见表１和表２，仿真时间２０Ｓ。文中分别采用串行和分网并行计算进行比较，其中有外接ＵＤ并行方案中电网被划分为８个子，河南电网被划分为豫西、豫北电网，剩余部分电网被划分到子网２中。自定义模型共有１２１８个继保与安自装置，其中因为网络分割共有８７个模型的输入输出来源于不同的子，因此被划分为３组，放入３个外接进程中。仿真效率对比结果如表３所示。表３连锁故障仿真效率测试结果Ｔａｂｌｅ３Ｒｅｓｕｌｔｓｆｏｒｓｉｍｕｌａｔｉｎｇｅｆｆｉｃｉｅｎｃｙｔｅｓｔｉｎｇ从表３的结果可以看出，将用户白定义模型置入外接进程中仿真，并通过ＭＰＩ接口与仿真主程序接口，通信带来的效率减低影响并不大，只有１．３９Ｓ，占仿真总时间的０．８％。但在此案例中，采用分网并行计算加速比并不高，这是因为所有自定义模型均集中安装在河南电，因此计算的效率实际上被河南电网所在子网计算效率所限制。徐得超，等含大量继保和安自装置的连锁故障并行仿真软件研发ｌ一自Ｐｌ一一无功ＤＩ－．’Ｉ０ｔ／ｓｆａＨｌｌＪｑ线－－－＂ｈｖ矧ＪＪＰ．１电流，ＩＥ－ｔ』『、、、（ｂ）邵淮线●＿一ＰＰ．一无功口口一－电Ｋ湖』．＼ｔ｜ｓ（ｃ）薛桐线ｌ一第１轮切鲍荷’｛一第２轮切负荷ｆ第３轮切荷ｔ／ｓ（ｄ）切负荷图６周口地区连锁故障仿真结果Ｆｉｇ．６ＣａｓｃａｄｉｎｇｆａｉｌｕｒｅｓｉｍｕｌａｔｉｎｇｒｅｓｕｌｔｓｆｏｒＺｈｏｕｋｏｕ尽管如此，采用外接进程仿真自定义模型可以解决输入输出跨子网的模型的仿真问题，并能实现不同进程间数值积分的部分并行计算，对提升程序效率有积极作用。５结论本文针对含有大量继电保护和安全自动装置的大规模电网连锁故障仿真，提出一种基于ＭＰＩ的并行仿真软件构建方法。在该方法中，大量继电保护和安全自动装置采用面向对象设计的用户自定义建模方法实现，以便于继保和安自装置模型库的维护与扩展；大规模电网和大量采用用户自定义建模的二次控制装置按照分网并行和分组并行思想，划分为子网进程和外接进程，一起并行仿真计算。案例分析表明，该并行仿真方法能够有效模拟大规模电网的连锁故障过程，并能有效提高仿真效率，使用户在可接受的时间内完成连锁故障仿真。限于实际电网继保和安自装置整定值的收资与建模工作繁琐以及工作量巨大，本文只在河南电网主网和部分地区２２０ｋＶ电网部署了继保和安自装置模型，对于大规模电网全网部署模型的效果尚待实际案例检验。参考文献［１］ＫＵＮＤＵＲＰＰｏｗｅｒｓｙｓｔｅｍｓｔａｂｉｌｉｔｙａｎｄｃｏｎｔｒｏｌ［Ｍ］．ＮｅｗＹｏｒｋ，ＵＳＡ：ＭｅＧｒａｗ－Ｈｉｌｌ，１９９４．［２］于洋，黄民翔，辛焕海．基于动态仿真的连锁故障分析【Ｊ】．电力系统自动化，２００８，３２（２０）：１５－２１．—ＹＵＹａｎｇ，ＨＵＡＮＧＭｉｎ・ｘｉａｎｇ，ＸＩＮＨｕａｎｈａｉ．Ｃａｓｃａｄｉｎｇｏｕｔａｇｅａｎａｌｙｓｉｓｂａｓｅｄｏｎｄｙｎａｍｉｃｓｉｍｕｌａｔｉｏｎ［Ｊ］．ＡｕｔｏｍａｔｉｏｎｏｆＥｌｅｃｔｒｉｃＰｏｗｅｒＳｙｓｔｅｍｓ，２００８，３２（２０）：１５－２１．［３］汤涌．电力系统数字仿真技术的现状与发展［Ｊ］．电力—系统自动化，２００２，２６（１７）：６６７０．ＴＡＮＧＹｏｎｇ．Ｐｒｅｓｅｎｔｓｉｔｕａｔｉｏｎａｎｄｄｅｖｅｌｏｐｍｅｎｔｏｆｐｏｗｅｒｓｙｓｔｅｍｓｉｍｕｌａｔｉｏｎｔｅｃｈｎｏｌｏｇｉｅｓ［Ｊ］．ＡｕｔｏｍａｔｉｏｎｏｆＥｌｅｃｔｒｉｃＰｏｗｅｒＳｙｓｔｅｍｓ，２００２，２６（１７）：６６－７０．［４］朱雨辰，赵冬梅，刘世良．基于机电暂态．电磁暂态混合仿真的电网合环分析计算系统［Ｊ］．电力系统保护与控制，２０１２，４０（２３）：７３．７９．—ＺＨＵＹｕ－ｃｈｅｎ，ＺＨＡＯＤｏｎｇｍｅｉ，ＬＩＵＳｈｉ－ｌｉａｎｇ．Ｌｏｏｐｃｌｏｓｉｎｇａｎａｌｙｔｉｃａｌｃａｌｃｕｌａｔｉｏｎｓｙｓｔｅｍｂａｓｅｄｏｎｅｌｅｃｔｒｏｍａｇｎｅｔｉｃ－ｅｌｅｃｔｒｏｍｅｃｈａｎｉｃａｌｔｒａｎｓｉｅｎｔｓｉｍｕｌａｔｉｏｎｆｏｒｐｏｗｅｒｎｅｔｗｏｒｋ［Ｊ］．￣．ｏｗｅｒＳｙｓｔｅｍＰｒｏｔｅｃｔｉｏｎａｎｄＣｏｎｔｒｏｌ，２０１２，４０（２３）：７３－７９．［５］孙元章，杨军，张晓东．电力系统动态仿真中的继电保护模型［Ｊ］．电力系统自动化，２００９，３３（２０）：４７－５１．ＳＵＮＹｕａｎ－ｚｈａｎｇ，ＹＡＮＧＪｕｎ，ＺＨＡＮＧＸｉａｏ－ｄｏｎｇ．Ｒｅｌａｙｐｒｏｔｅｃｔｉｏｎｍｏｄｅｌｆｏｒｄｙｎａｍｉｃｓｉｍｕｌａｔｉｏｎｏｆｐｏｗｅｒｓｙｓｔｅｍ［Ｊ］．ＡｕｔｏｍａｔｉｏｎｏｆＥｌｅｃｔｒｉｃＰｏｗｅｒＳｙｓｔｅｍｓ，２００９，３３（２０）：４７－５１．Ｏ５Ｏ５Ｏ５Ｏ５Ｏ５Ｏ５４４３３２２●ＯＯ．１３８．电力系统保护与控制［６］汤必强，李峰，於益军．计及自动装置模拟的大电网在线静态安全分析［Ｊ】．电力系统自动化，２０１１，３５（１４）：１Ｏ７一ｌ１１．ＴＡＮＧＢｉ－ｑｉａｎｇ，ＬＩＦｅｎｇ，ＹＵＹｉ－ｊｕｎ．Ｏｎｌｉｎｅｓｔａｔｅｓｅｃｕｒｉｔｙａｎａｌｙｓｉｓｃｏｎｓｉｄｅｒｉｎｇ［Ｊ］．ＡｕｔｏｍａｔｉｏｎｏｆＥｌｅｃｔｒｉｃＰｏｗｅｒＳｙｓｔｅｍｓ，２０１１，３５（１４）：１０７－１１１．［７］王增平，姜宪国，张执超，等．智能电网环境下的继电保护［Ｊ］．电力系统保护与控制，２０１３，４１（２）：１３．１８．—ＷＡＮＧＺｅｎｇ・ｐｉｎｇ，ＪＩＡＮＧＸｉａｎ－ｇｕｏ，ＺＨＡＮＧＺｈｉｃｈａｏ，ｅｔａ１．Ｒｅｌａｙｐｒｏｔｅｃｔｉｏｎｆｏｒｓｍａｒｔｇｒｉｄ［Ｊ］．ＰｏｗｅｒＳｙｓｔｅｍ—ＰｒｏｔｅｃｔｉｏｎａｎｄＣｏｎｔｒｏｌ，２０１３，４ｌ（２）：１３１８．［８］唐艳茹．区域系统安全自动装置专题研究报告［Ｊ］．继电器，１９９９，２７（４）：３０－３４．ＴＡＮＧＹａｎ・ｒｕ．Ｓｔｕｄｙｒｅｐｏｒｔｏｎｔｈｅｒｅｇｉｏｎａｌｓｙｓｔｅｍｓｅｃｕｒｉｔｙａｕｔｏｍａｔｉｃｄｅｖｉｃｅ［Ｊ］．Ｒｅｌａｙ，１９９９，２７（４）：３０－３４．［９］吴国呖，宋新立，汤涌．电力系统动态仿真中的安全稳定装置系统建模［Ｊ］．电力系统自动化，２０１２，３６（３）：７１．７５．—ＷＵＧｕｏ－ｙａｎｇ，ＳＯＮＧＸｉｎｌｉ，ＴＡＮＧＹｏｎｇ．Ｍｏｄｅｌｉｎｇｏｆｓｅｃｕｒｉｔｙａｎｄｓｔａｂｉｌｉｔｙｃｏｎｔｒｏｌｓｙｓｔｅｍｆｏｒｐｏｗｅｒｓｙｓｔｅｍｄｙｎａｍｉｃｓｉｍｕｌａｔｉｏｎ［Ｊ］．ＡｕｔｏｍａｔｉｏｎｏｆＥｌｅｃｔｒｉｃＰｏｗｅｒ—Ｓｙｓｔｅｍｓ，２０１２，３６（３）：７１７５．［１Ｏ］王治国，李兴建，王言国．基于统一建模的继电保护测试装置开发研究【Ｊ】．电力系统保护与控制，２０１０，—３８（１９）：１８０１８５．—ＷＡＮＧＺｈｉ－ｇｕｏ，ＬＩＸｉｎｇ－ｊｉａｎ，ＷＡＮＧＹａｎｇｕｏ．Ｒｅｓｅａｒｃｈｏｎｔｈｅｄｅｖｅｌｏｐｍｅｎｔｏｆｒｅｌａｙｐｒｏｔｅｃｔｉｏｎｔｅｓｔｅｒｂａｓｅｄｏｎｕｎｉｆｉｅｄｂｕｉｌｄｉｎｇｍｏｄｅ［Ｊ］．ＰｏｗｅｒＳｙｓｔｅｍＰｒｏｔｅｃｔｉｏｎａｎｄＣｏｎｔｒｏｌ，２０１０，３８（１９）：１８０－１８５．［１１］吴中习，周孝信，包忠明．暂态稳定计算中的用户自定义模型［Ｊ】．电网技术，１９９１，１５（３）：４６．５０．—ＷＵＺｈｏｎｇ－ｘｉ，ＺＨＯＵＸｉａｏｘｉｎ，ＢＡＯＺｈｏｎｇ－ｍｉｎｇ．Ｕｓｅｒ－ｄｅｆｉｎｅｄｍｏｄｅｌｉｎｇｉｎｐｏｗｅｒｓｙｓｔｅｍｔｒａｎｓｉｅｎｔｓｉｍｕｌａｔｉｏｎ［Ｊ］．ＰｏｗｅｒＳｙｓｔｅｍＴｅｃｈｎｏｌｏｇｙ，１９９１，１５（３）：４６．５Ｏ．［１２］孟新军，徐得超，李亚楼．ＰＳＡＳＰ用户自定义建模优化—算法研究［Ｊ］．电网技术，２０１２，３６（７）：２１４２１８．ＭＥＮＧＸｉｎ－ｊｕｎ，ＸＵＤｅ・ｃｈａｏ，ＬＩＹａ－ｌｏｕ．Ｒｅｓｅａｒｃｈｏｎｏｐｔｉｍｉｚａｔｉｏｎａｌｇｏｒｉｔｈｍｆｏｒｕｓｅｒ－ｄｅｆｉｎｅｄｍｏｄｅｌｉｎｇｏｆｐｏｗｅｒｓｙｓｔｅｍａｎａｌｙｓｉｓｓｏｆｔｗａｒｅｐａｃｋａｇｅ［Ｊ］．Ｐｏｗｅｒ—ＳｙｓｔｅｍＴｅｃｈｎｏｌｏｇｙ，２０１２，３６（７）：２１４２１８．［１３］黄英怀．电力系统分析综合程序的用户自定义模型［Ｊ］．广东电力，２００２，１５（５）：４６－４８．ＨＵＡＮＧＹｉｎｇ－ｈｕａｉ．Ｕｓｅｒ－ｄｅｆｉｎｅｄｍｏｄｅｌｏｆｐｏｗｅｒｓｙｓｔｅｍａｎａｌｙｓｉｓｓｏｆｔｗａｒｅｐａｃｋａｇｅ［Ｊ］．ＧｕａｎｇｄｏｎｇＥｌｅｃｔｒｉｃＰｏｗｅｒ，—２００２，１５（５）：４６４８．［１４］杨晓东，黄东平，陈世勇．基于全图形建模方式的电力系统仿真系统继电保护库的设计与应用［Ｊ］．电力系统保护与控制，２０１３，４１（１５）：１１６．１１９．———ＹＡＮＧＸｉａｏｄｏｎｇ，ＨＵＡＮＧＤｏｎｇｐｉｎｇ，ＣＨＥＮＳｈｉｙｏｎｇ．Ｄｅｓｉｇｎａｎｄａｐｐｌｉｃａｔｉｏｎｏｆｔｈｅｒｅｌａｙｐｒｏｔｅｃｔｉｏｎｍｏｄｕｌｅｌｉｂｒａｒｙｂａｓｅｄｏｎｔｈｅｆｕｌｌｇｒａｐｈｉｃｓｍｏｄｅｌｉｎｇｉｎｐｏｗｅｒｓｉｍｕｌａｔｉｏｎｓｙｓｔｅｍ［Ｊ］．ＰｏｗｅｒＳｙｓｔｅｍＰｒｏｔｅｃｔｉｏｎａｎｄ—Ｃｏｎｔｒｏｌ，２０１３，４１（１５）：１１６１１９．［１５］李芳，郭剑，吴中习．基于ＰＣ机群的电力系统小干扰稳定分布式并行算法【Ｊ】．中国电机工程学报，２００７，２７（３１、：７－１３．ＬＩＦａｎｇ，ＧＵＯＪｉａｎ，ＷＵＺｈｏｎｇ－ｘｉ．ＤｉｓｔｒｉｂｕｔｅｄｐａｒａｌｌｅｌｃｏｍｐｕｔｉｎｇａｌｇｏｒｉｔｈｍｓｆｏｒｐｏｗｅｒｓｙｓｔｅｍｓｍａｌｌｓｉｇｎａｌｓｔａｂｉｌｉｔｙｂａｓｅｄｏｎＰＣｃｌｕｓｔｅｒ［Ｊ］．ＰｒｏｃｅｅｄｉｎｇｓｏｆｔｈｅＣＳＥＥ，２００７，２７（３１）：７－１３．［１６］李芳．电力系统小干扰稳定分布式并行算法研究［Ｄ】．北京：中国电力科学研究院，２００６．［１７］李亚楼．大规模电力系统机电暂态实时仿真算法及软件的研究【Ｄ】．北京：中国电力科学研究院，２００３．——收稿日期：２０１４－０２２６；修回日期：２０１４－０５０５作者简介：徐得超（１９７５一），男，硕士，高工，主要从事电力系统数字仿真技术，配电自动化应用软件研究等；Ｅ．ｍａｉｌ：ｘｄｃ＠ｅｐｒｉ．ｓｇｃｃ．ｃｏｍ．ｃｎ刘巍（１９７７－），女，硕士，高工，研究方向为电力系统计算与仿真分析等；朱旭凯（１９８２－），男，博士后，高工，研究方向为电力系统连锁故障、混合仿真技术。