CN110196997A - 一种电磁暂态异步并行计算方法及装置 - Google Patents

一种电磁暂态异步并行计算方法及装置 Download PDF

Info

Publication number
CN110196997A
CN110196997A CN201910361932.1A CN201910361932A CN110196997A CN 110196997 A CN110196997 A CN 110196997A CN 201910361932 A CN201910361932 A CN 201910361932A CN 110196997 A CN110196997 A CN 110196997A
Authority
CN
China
Prior art keywords
subnet
parallel computation
time
calculating
constraint condition
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201910361932.1A
Other languages
English (en)
Other versions
CN110196997B (zh
Inventor
刘文焯
李霞
穆世霞
叶小晖
许克
许冰
宋新立
吴国旸
苏志达
刘涛
戴汉扬
左克江
肖雄
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Electric Power Research Institute Co Ltd CEPRI
State Grid Shanghai Electric Power Co Ltd
Original Assignee
China Electric Power Research Institute Co Ltd CEPRI
State Grid Shanghai Electric Power Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China Electric Power Research Institute Co Ltd CEPRI, State Grid Shanghai Electric Power Co Ltd filed Critical China Electric Power Research Institute Co Ltd CEPRI
Priority to CN201910361932.1A priority Critical patent/CN110196997B/zh
Publication of CN110196997A publication Critical patent/CN110196997A/zh
Application granted granted Critical
Publication of CN110196997B publication Critical patent/CN110196997B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Energy or water supply

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Health & Medical Sciences (AREA)
  • Economics (AREA)
  • General Physics & Mathematics (AREA)
  • Water Supply & Treatment (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Strategic Management (AREA)
  • Tourism & Hospitality (AREA)
  • Human Resources & Organizations (AREA)
  • General Business, Economics & Management (AREA)
  • Marketing (AREA)
  • Public Health (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Data Exchanges In Wide-Area Networks (AREA)

Abstract

本发明公开了一种电磁暂态异步并行计算方法，包括：对大电网进行网络分割，确定与每个子网相连的其他所有子网；并根据传输线理论计算任意两个相连子网间的最小延时；判断任意两个相连子网之间的计算时间差是否在最小延时范围内，并将所述判断条件作为任意子网异步并行计算的时间约束条件；当任意子网计算完毕，根据所述约束条件，判断所述任意子网是否满足异步并行计算的时间约束条件，若所述任意子网不满足异步并行计算条件，则需要挂起等待，当收到解锁通知后进行所述约束条件的再次判断；当所述任意子网满足异步并行计算的时间约束条件，则所述任意子网进入下一时刻的计算。本发明解决了现有电磁暂态仿真程序计算效率低的问题。

Description

一种电磁暂态异步并行计算方法及装置
技术领域
本申请涉及电力系统的仿真领域，具体涉及一种基于大电网网络分割的电磁暂态异步并行计算方法，同时涉及一种基于大电网网络分割的电磁暂态异步并行计算装置。
背景技术
电力系统动态仿真是大规模电力系统动态特性机理研究、严重事故特征分析以及稳定措施研究的重要方法。随着电网规模的不断地扩大，电力系统仿真计算任务也越来越重。且随着电力系统大量高压直流设备、新能源、柔性交直流输电系统的接入，传统的机电暂态仿真已经无法满足大规模系统对仿真精度的要求，电磁暂态仿真逐渐成为精确模拟当今及未来电网的有效手段。但是传统电磁暂态存在步长小，仿真速度慢的缺点。
目前，电力系统仿真方法可以分为串行计算方法和并行计算方法两种。其中，并行计算主要包括空间并行算法、时间并行算法和时间空间组合并行算法。空间并行算法的核心思想是将电力系统按拓扑结构划分为多个子区域，而每个子区域分配给一个进(线)程进行单步求解。时间并行算法是将微分方程组在多个时步上的差分方程组联立。时间空间组合算法是将空间并行和时间并行相结合的一种并行算法。上述并行算法均属于同步并行范畴，虽然异步并行在电力系统动态仿真方面也有一定的研究，但是目前在国内外实用程序中都不具备异步并行计算功能。
发明内容
本申请提供一种电磁暂态异步并行计算方法，解决了现有电磁暂态仿真程序计算效率低的问题。
本申请提供一种电磁暂态异步并行计算方法，其特征在于，包括：
对大电网进行网络分割，确定与每个子网相连的其他所有子网；并根据传输线理论计算任意两个相连子网间的最小延时；
判断任意两个相连子网之间的计算时间差是否在最小延时范围内，并将所述判断条件作为任意子网异步并行计算的时间约束条件；
当任意子网计算完毕，并记录计算完毕的时刻，根据任意子网异步并行计算的时间约束条件，判断所述任意子网是否满足异步并行计算的时间约束条件，若所述任意子网不满足异步并行计算的时间约束条件，则需要挂起等待，当收到解锁通知后进行所述约束条件的再次判断；当所述任意子网满足异步并行计算的时间约束条件时，则所述任意子网进入下一时刻的计算。
优选的，所述若所述任意子网不满足异步并行计算的时间约束条件，则需要挂起等待，当收到解锁通知后进行所述约束条件的再次判断，包括：
若所述任意子网不满足异步并行计算的时间约束条件，则所述任意子网挂起等待，与所述任意子网相连的其他子网在设置的时刻若满足条件则继续进行计算，当与所述任意子网相连的其他任一子网计算完毕，则主动发出解锁通知，并更新任一相连子网的计算时刻；当所述处于挂起等待的任意子网接收到相连子网的解锁通知后，进行所述约束条件的再次判断。
优选的，在对大电网进行网络分割，确定与每个子网相连的其他所有子网的步骤后，还包括：将所述相连的子网组成对应的相连子网集。
优选的，所述根据传输线理论计算任意两个相连子网间的最小延时，包括：
根据传输线理论，两个相连子网间包括多个传输线延时τ1，τ2，…，τm，其中，m为正整数，且满足m≥2；
在所述多个传输线延时中，求取最小延时τ的公式为，
τ＝min{τ1,τ2,...,τm}；
根据求取最小延时τ的公式，分别获得子网i与相连子网i1，i2，…，ik，…，in之间的最小延时分别为τi1、τi2、…、τik、…、τin。
优选的，所述判断任意两个相连子网之间的计算时间差是否在最小延时范围内，包括：
根据传输线理论，对于任意一个计算子网络，必须使用与其相连的其它子网络的历史数据，对于子网络i的计算，必须要使用相连子网i1，i2，…，ik，…， in的历史数据，假设子网络i，i1，i2，…，ik，…，in的计算时刻分别为Ti，Ti1，Ti2，…， Tin，则需要满足以下条件，
从上面公式可以得出，任意两个相连子网之间的计算时间差，必须要在两者之间的最小延时范围内，并将所述判断条件作为任意子网异步并行计算的时间约束条件。
优选的，当任意子网在时刻计算完毕的步骤之后，还包括：
所述任意子网主动通知各相连的子网，完成通知解锁，并更新自身的计算时刻。
优选的，若所述任意子网满足异步并行计算的时间约束条件，则所述任意子网进入下一时该的计算，包括：
当子网i在ti时刻计算完毕后，若子网i满足异步并行计算条件，则进入ti+ ΔT时刻计算。
优选的，当与任意子网相连的其他任一子网计算完毕，则所述任意子网进入下一时刻计算的步骤之后，还包括：
所述任意子网发出计算完毕的解锁通知信号，通知与其相连的子网。
本申请同时提供一种电磁暂态异步并行计算装置，其特征在于，包括：
系统初始化计算单元，对大电网进行网络分割，确定与每个子网相连的其他所有子网；并根据传输线理论计算任意两个相连子网间的最小延时；
约束条件确定单元，判断任意两个相连子网之间的计算时间差是否在最小延时范围内，并将所述判断条件作为任意子网异步并行计算的时间约束条件；
并行计算执行单元，当任意子网计算完毕，并记录计算完毕的时刻，根据任意子网异步并行计算的时间约束条件，判断所述任意子网是否满足异步并行计算的时间约束条件，若所述任意子网不满足异步并行计算的时间约束条件，则需要挂起等待，当收到解锁通知后进行所述约束条件的再次判断；当所述任意子网满足异步并行计算的时间约束条件时，则所述任意子网进入下一时刻的计算。
优选的，所述并行计算执行单元，包括：
挂起等待子单元，若所述任意子网不满足异步并行计算的时间约束条件，则所述任意子网挂起等待，与所述任意子网相连的其他子网在设置的时刻若满足条件则继续进行计算，当与所述任意子网相连的其他任一子网计算完毕，则主动发出解锁通知，并更新任一相连子网的计算时刻；当所述处于挂起等待的任意子网接收到相连子网的解锁通知后，进行所述约束条件的再次判断。
优选的，还包括：
通知解锁单元，所述任意子网计算完毕后会主动通知各相连的子网，完成通知解锁；
时间更新单元，当所述任意子网计算完毕，则更新所述任意子网的计算时刻。
本申请提供一种电磁暂态异步并行计算方法，通过对大电网的网络进行分割，采用挂起等待/通知解锁模式实现了各子网之间的异步并行计算，提高了仿真程序计算效率，满足现有大电网对仿真效率的需求。
附图说明
图1是本申请实施例提供的一种电磁暂态异步并行计算方法；
图2是本申请实施例涉及的大电网的网络分割示意图；
图3是本申请实施例涉及的大电网的各子网之间的时序图；
图4是本申请实施例涉及的一种电磁暂态异步并行计算方法的步骤示意图；
图5是本申请实施例提供的一种电磁暂态异步并行计算装置。
具体实施方式
在下面的描述中阐述了很多具体细节以便于充分理解本申请。但是本申请能够以很多不同于在此描述的其它方式来实施，本领域技术人员可以在不违背本申请内涵的情况下做类似推广，因此本申请不受下面公开的具体实施的限制。
请参看图1，图1是本申请实施例提供的一种电磁暂态异步并行计算方法，下面结合图1对申请提供的方法进行详细说明。
步骤S101，对大电网进行网络分割，确定与每个子网相连的其他所有子网；并根据传输线理论计算任意两个相连子网间的最小延时。
如图2所示，对大电网进行网络分割，分为子网络i，i1，i2，…，ik，…， in等多个子网，n为正整数。
确定与每个子网相连的其他所有子网，将所述相连的子网组成对应的相连子网集。由图2可以，与子网i相连的各子网分别为i1，i2，…，ik，…，in；同理，子网i1的相连子网分别为子网i，i2；依次确认各子网的相连子网，并形成相应的相连子网集。对于子网i的相连子网集为：{子网i1，i2，…，ik，…，in}。然后根据传输线理论计算任意两个相连子网间的最小延时。还是根据图2所示，根据传输线理论，相连两个子网间(子网i与相连子网i1)包括多个传输线延时，τ1，τ2，…，τm，其中，m为正整数，且满足m≥2；在所述多个传输线延时中，求取最小延时τ的公式为，τ＝min{τ1,τ2,...,τm}；根据求取最小延时τ的公式，获得子网i与相连子网i1，i2，…，ik，…，in之间的最小延时分别为τi1、τi2、…、τik、…、τin。
步骤S102，判断任意两个相连子网之间的计算时间差是否在最小延时范围内，并将所述判断条件作为任意子网异步并行计算的时间约束条件。
根据传输线理论，对于任意一个计算子网络，必须使用与其相连的其它子网络的历史数据，如图2所示的网络，对于子网络i的计算，必须要使用相连子网i1，i2，…，ik，…，in的历史数据，假设子网络i，i1，i2，…，ik，…，in的计算时刻分别为Ti，Ti1，Ti2，…，Tin，则需要满足以下条件，
从上面公式可以得出，任意两个相连子网之间的计算时间差，必须要在两者之间的最小延时范围内。
然后将所述判断条件作为任意子网异步并行计算的时间约束条件。
步骤S103，当任意子网计算完毕，并记录计算完毕的时刻，根据任意子网异步并行计算的时间约束条件，判断所述任意子网是否满足异步并行计算的时间约束条件，若所述任意子网不满足异步并行计算的时间约束条件，则需要挂起等待，当收到解锁通知后进行所述约束条件的再次判断；当所述任意子网满足异步并行计算的时间约束条件时，则所述任意子网进入下一时刻的计算。
当任意子网计算完毕后，并记录计算完毕的时刻，所述任意子网主动通知各相连的子网，完成通知解锁。因此，在Ti时刻计算完成后，计算子网i会主动发出计算完毕的通知信号，通知其他所有相连子网当前计算已完成。如图3 所示，当子网i在Ti时刻计算完毕后，需要判断所述子网i是否满足异步并行计算的时间约束条件，以便进入Ti+ΔT时刻计算。判断所述任意子网是否满足异步并行计算的时间约束条件，当子网i在Ti时刻计算完毕后，若子网i满足异步并行计算条件，则进入Ti+ΔT时刻计算。也就是判断子网i与其相连的各子网i1，i2，…，ik，…，in的时间差值是否在最小延时τi1、τi2、…，τik、…、τin范围内，满足步骤S102中约束条件的公式。
若所述任意子网不满足异步并行计算的时间约束条件，则所述任意子网挂起等待，与所述任意子网相连的其他子网继续进行计算。当子网i不满足异步并行计算的约束条件时，需要进入挂起等待模式。在本实施例中，假如计算子网i与任意相连子网ik之间的计算时间差值大于两者之间的最小延时，即不满足-τik≤Ti-Tik≤τik，则子网i需要循环等待，进入挂起等待模式。典型地在本实施例中，可以通过调用C++11新标准库中的std::condition_variable的成员函数wait()实现等待。
如果有子网计算完毕，会发出计算完毕的通知信号，主动通知与其相连的其他所有子网，并更新自身的计算时刻值。在本实施例中，当上述任意子网ik计算完毕后，会主动通知与其相连的子网，即子网ik计算完毕后自动发出通知，并更新计算时刻值Tik，完成通知解锁。当子网i接收到相连子网ik的通知信号后，进入解锁模式。典型地可以采用C++11新标准库std::condition_variable的成员函数notify_all()完成通知解锁。
对于接收到相连子网通知解锁信号的子网i，需要进行下一步条件判断，即计算子网i进行异步并行计算约束条件的重新判断。
判断所述任意子网是否满足约束条件，当满足异步并行计算的约束条件时，即当计算子网与相连子网之间的计算时间差值在最小延时范围内，便可进入下一时刻的计算。对于任意计算子网i在Ti时刻计算结束后，判断当满足异步并行计算条件后，则计算子网i自动进入Ti+Δt时刻的计算。
图4是本发明实施例涉及的一种电磁暂态异步并行计算方法的步骤示意图，从图中可以看出，将本发明的方法的步骤分为三个部分，包括初始化计算、异步并行原则、异步并行计算实现方法。
初始化计算分为三个步骤，步骤101：大电网的网络分割，即分网。根据实际计算系统的特性，可以采用相应的分网方法，典型的可采用传输线分网技术。
步骤102：确定与子网i相连的所有子网i1，i2，…，ik，…，in。对分网后的各子网络进行统计分析，确定各计算子网相连的所有其他子网络，形成相连子网集。其余为不相连子网络。不相连的两个网络之间的计算独立进行，互不约束。
步骤103：根据传输线理论计算最小延时。
对于相连的两个子网，子网之间通过传输线相连，根据传输线理论，可知两子网之间包括多个传输线延时τ1，τ2，…，τm。其中，m为正整数，且满足m≥2。
在上述计算得到的所有延时中，求取最小延时τ，即：
τ＝min{τ1,τ2,...,τm} (1)
根据公式(1)得到任意计算子网i与其各相连子网之间的最小延时，即计算子网i与相连子网i1，i2，…，ik，…，in之间的最小延时分别为τi1、τi2、…、τik、…、τin。
优选的，所述异步并行原则，包括：
步骤201：得到子网i异步并行计算的时间约束条件。
各子网之间通过传输线相连，根据传输线计算理论可知，对于任意一个计算子网络，必须使用与其相连的其它子网络的历史数据。即对于子网络i的计算，必须要使用相连子网i1，i2，…，ik，…，in的历史数据，假设子网络i，i1， i2，…，ik，…，in的计算时刻分别为Ti，Ti1，Ti2，…，Tin，则需要满足以下条件：
即任意两个相连子网络之间的计算时间差，必须要在两者之间的最小延时范围内，这就是该异步并行计算的时间约束条件，也是该异步并行计算的基本原则。
优选的，所述异步并行计算实现方法，包括：
步骤301：判断计算子网i的Ti时刻是否计算完毕。
当子网计算完毕后，均会主动通知各相连子网，完成通知解锁。因此在上一时刻计算完成后，计算子网i会发出计算完毕的通知信号，通知其他所有相连子网Ti时刻计算完成。
步骤302：判断Ti+Δt时刻是否满足约束条件，即是否满足异步并行计算的原则。
当计算子网i计算完毕后，需要判断是否满足异步并行计算条件，以便进入下一时刻计算。异步并行的约束条件如所述步骤201，即当计算子网i在Ti时刻计算完毕后，判断计算子网i与其相连的各子网i1，i2，…，ik，…，in的时间差值是否在最小延时τi1、τi2、…，τik、…、τin范围内，满足公式(2)。
步骤303：子网i等待通知。
对不满足异步并行计算约束条件的子网i进行挂起等待。所述步骤302中，需要判断是否满足异步并行计算的约束条件，对于任意不满足异步并行条件的子网i，需要进行等待，即该计算子网i进入了挂起等待模式。特别地在本发明中，优先推荐采用C++11新标准库，调用std::condition_variable的成员函数wait() 实现等待。
步骤304：其他相连子网的解锁通知。
所述步骤303中，对于不满足异步并行条件的计算子网进行挂起等待，并判断与其相连的其他子网是否计算完毕。如果有子网计算完毕，，会发出计算完毕的通知信号，通知与其相连的子网，并更新自身的计算时刻。
特别地在本发明中，优先推荐采用C++11新标准库std::condition_variable 的成员函数notify_all()，当任意子网计算完毕，可采用上述函数notify_all()，主动通知等待的子网的线程，即完成通知解锁。
对于接收到相连子网通知解锁信号的子网i，需要进行约束条件的再次判断，即计算子网i进入上述步骤302，进行异步并行计算约束条件的重新判断。
步骤305：对满足约束条件的子网络，开始Ti+Δt时刻的计算，即进行下一时刻计算。
所述步骤303中，判断计算子网是否满足约束条件，当满足异步并行计算的约束条件时，即当计算子网与相连子网之间的计算时间差值在最小延时范围内，便可进入下一时刻的计算。对于任意计算子网i在Ti时刻计算结束后，判断当满足异步并行计算条件后，则计算子网i自动进入Ti+Δt时刻的计算。
优选的，异步并行的实现方法，典型的计算子网i的时序图如附图3所示，首先判断计算子网i的Ti时刻计算是否结束；当Ti时刻计算结束后，判断Ti+ ΔT时刻是否满足相连子网络之间的时间约束条件，即计算子网i与相连子网之间的计算时间差是否在两相连网络的最小延时范围内；如果计算子网络i不满足以上条件，则调用wait()函数进入挂起等待模式。其他相连子网i1，i2，…， ik，…，in计算结束后，更新计算时刻，并主动向子网i发送解锁通知信号。当子网i接收到相连子网的解锁通知信号后，并根据更新后的计算时刻值，子网 i重新判断异步并行计算条件是否满足，当满足条件后，子网i进入Ti+ΔT时刻计算。
与本实施例提供的一种电磁暂态异步并行计算方法相对应的，本申请同时提供一种电磁暂态异步并行计算装置500，其特征在于，包括：
系统初始化计算单元510，对大电网进行网络分割，确定与每个子网相连的子网；并根据传输线理论计算任意两个相连子网间的最小延时；
约束条件确定单元520，判断任意两个相连子网之间的计算时间差是否在最小延时范围内，并将所述判断条件作为任意子网异步并行计算的时间约束条件；
并行计算执行单元530，当任意子网计算完毕，并记录计算完毕的时刻，根据任意子网异步并行计算的时间约束条件，判断所述任意子网是否满足异步并行计算的时间约束条件，若所述任意子网不满足异步并行计算的时间约束条件，则需要挂起等待，当收到解锁通知后进行所述约束条件的再次判断；当所述任意子网满足异步并行计算的时间约束条件时，则所述任意子网进入下一时刻的计算。
优选的，所述并行计算执行单元，包括：挂起等待子单元，若所述任意子网不满足异步并行计算的时间约束条件，则所述任意子网挂起等待，与所述任意子网相连的其他子网在设置的时刻若满足条件则继续进行计算，当与所述任意子网相连的其他任一子网计算完毕，则主动发出解锁通知，并更新任一相连子网的计算时刻。当所述处于挂起等待的任意子网接收到相连子网的解锁通知后，进行所述约束条件的再次判断。
优选的，还包括：通知解锁单元，所述任意子网计算完毕后会主动通知各相连的子网，完成通知解锁。
时间更新单元，当所述任意子网计算完毕，则更新所述任意子网的计算时刻。
本申请提供的一种电磁暂态异步并行计算方法，通过对大电网的网络进行分割，采用挂起等待/通知解锁模式实现了各子网之间的异步并行计算，提高了仿真程序计算效率，满足现有大电网对仿真的需求。与现有技术相比，可缩短电力系统电磁暂态仿真的计算时间，有效提高电磁暂态仿真的计算效率。本发明有效地利用C++11新标准库，不需要戴维南等值，采用挂起等待和通知解锁模式，提高了电磁暂态仿真精度并简化了建模过程，编程简单，易于实现。
以上实施例仅用以说明本发明的技术方案而非对其限制，尽管参照上述实施例对本发明进行了详细的说明，所属领域的普通技术人员依然可以对本发明的具体实施方式进行修改或者等同替换，而这些未脱离本发明精神和范围的任何修改或者等同替换，其均在申请待批的本发明的权利要求保护范围之内。

Claims (11)

1.一种电磁暂态异步并行计算方法，其特征在于，包括：
对大电网进行网络分割，确定与每个子网相连的其他所有子网；并根据传输线理论计算任意两个相连子网间的最小延时；
判断任意两个相连子网之间的计算时间差是否在最小延时范围内，并将所述判断条件作为任意子网异步并行计算的时间约束条件；
当任意子网计算完毕，并记录计算完毕的时刻，根据任意子网异步并行计算的时间约束条件，判断所述任意子网是否满足异步并行计算的时间约束条件，若所述任意子网不满足异步并行计算的时间约束条件，则需要挂起等待，当收到解锁通知后进行所述约束条件的再次判断；当所述任意子网满足异步并行计算的时间约束条件时，则所述任意子网进入下一时刻的计算。
2.根据权利要求1所述的方法，其特征在于，所述若所述任意子网不满足异步并行计算的时间约束条件，则需要挂起等待，当收到解锁通知后进行所述约束条件的再次判断，包括：
若所述任意子网不满足异步并行计算的时间约束条件，则所述任意子网挂起等待，与所述任意子网相连的其他子网在设置的时刻若满足条件则继续进行计算，当与所述任意子网相连的其他任一子网计算完毕，则主动发出解锁通知，并更新任一相连子网的计算时刻；当所述处于挂起等待的任意子网接收到相连子网的解锁通知后，进行所述约束条件的再次判断。
3.根据权利要求1所述的方法，其特征在于，在对大电网进行网络分割，确定与每个子网相连的其他所有子网的步骤后，还包括：将所述相连的子网组成对应的相连子网集。
4.根据权利要求1所述的方法，其特征在于，所述根据传输线理论计算任意两个相连子网间的最小延时，包括：
根据传输线理论，两个相连子网间包括多个传输线延时τ1，τ2，…，τm，其中，m为正整数，且满足m≥2；
在所述多个传输线延时中，求取最小延时τ的公式为，
τ＝min{τ1,τ2,...,τm}；
根据求取最小延时τ的公式，分别获得子网i与相连子网i1，i2，…，ik，…，in之间的最小延时分别为τi1、τi2、…、τik、…、τin。
5.根据权利要求1所述的方法，其特征在于，所述判断任意两个相连子网之间的计算时间差是否在最小延时范围内，包括：
根据传输线理论，对于任意一个计算子网络，必须使用与其相连的其它子网络的历史数据，对于子网络i的计算，必须要使用相连子网i1，i2，…，ik，…，in的历史数据，假设子网络i，i1，i2，…，ik，…，in的计算时刻分别为Ti，Ti1，Ti2，…，Tin，则需要满足以下条件，
从上面公式可以得出，任意两个相连子网之间的计算时间差，必须要在两者之间的最小延时范围内，并将所述判断条件作为任意子网异步并行计算的时间约束条件。
6.根据权利要求1所述的方法，其特征在于，当任意子网计算完毕的步骤之后，还包括：
所述任意子网主动通知各相连的子网，完成通知解锁，并更新自身的计算时刻。
7.根据权利要求1所述的方法，其特征在于，若所述任意子网满足异步并行计算的时间约束条件，则所述任意子网进入下一时刻的计算，包括：
当子网i在ti时刻计算完毕后，若子网i满足并行计算条件，则进入ti+ΔT时刻计算。
8.根据权利要求1所述的方法，其特征在于，当与任意子网相连的其他任一子网计算完毕，则所述任意子网进入下一时刻计算的步骤之后，还包括：
所述任意子网发出计算完毕的解锁通知信号，通知与其相连的子网。
9.一种电磁暂态异步并行计算装置，其特征在于，包括：
系统初始化计算单元，对大电网进行网络分割，确定与每个子网相连的其他所有子网；并根据传输线理论计算任意两个相连子网间的最小延时；
约束条件确定单元，判断任意两个相连子网之间的计算时间差是否在最小延时范围内，并将所述判断条件作为任意子网异步并行计算的时间约束条件；
并行计算执行单元，当任意子网计算完毕，并记录计算完毕的时刻，根据任意子网异步并行计算的时间约束条件，判断所述任意子网是否满足异步并行计算的时间约束条件，若所述任意子网不满足异步并行计算的时间约束条件，则需要挂起等待，当收到解锁通知后进行所述约束条件的再次判断；当所述任意子网满足异步并行计算的时间约束条件时，则所述任意子网进入下一时刻的计算。
10.根据权利要求9所述的装置，其特征在于，所述并行计算执行单元，包括：
挂起等待子单元，若所述任意子网不满足异步并行计算的时间约束条件，则所述任意子网挂起等待，与所述任意子网相连的其他子网在设置的时刻若满足条件则继续进行计算，当与所述任意子网相连的其他任一子网计算完毕，则主动发出解锁通知，并更新任一相连子网的计算时刻；当所述处于挂起等待的任意子网接收到相连子网的解锁通知后，进行所述约束条件的再次判断。
11.根据权利要求9所述的装置，其特征在于，还包括：
通知解锁单元，所述任意子网计算完毕后会主动通知各相连的子网，完成通知解锁；
时间更新单元，当所述任意子网计算完毕，则更新所述任意子网的计算时刻。
CN201910361932.1A 2019-04-30 2019-04-30 一种电磁暂态异步并行计算方法及装置 Active CN110196997B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910361932.1A CN110196997B (zh) 2019-04-30 2019-04-30 一种电磁暂态异步并行计算方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910361932.1A CN110196997B (zh) 2019-04-30 2019-04-30 一种电磁暂态异步并行计算方法及装置

Publications (2)

Publication Number Publication Date
CN110196997A true CN110196997A (zh) 2019-09-03
CN110196997B CN110196997B (zh) 2024-04-26

Family

ID=67752293

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910361932.1A Active CN110196997B (zh) 2019-04-30 2019-04-30 一种电磁暂态异步并行计算方法及装置

Country Status (1)

Country Link
CN (1) CN110196997B (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6819685B1 (en) * 1998-12-03 2004-11-16 Alcatel Method of and system for controlling a frequency via an asynchronous transmission network and mobile telephone network including the system
CN102750416A (zh) * 2012-06-19 2012-10-24 中国电力科学研究院 一种含开关特性电路的电磁暂态仿真的拓扑分网方法
CN106886616A (zh) * 2015-12-15 2017-06-23 中国电力科学研究院 一种大规模电磁暂态电网仿真的自动分网方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6819685B1 (en) * 1998-12-03 2004-11-16 Alcatel Method of and system for controlling a frequency via an asynchronous transmission network and mobile telephone network including the system
CN102750416A (zh) * 2012-06-19 2012-10-24 中国电力科学研究院 一种含开关特性电路的电磁暂态仿真的拓扑分网方法
CN106886616A (zh) * 2015-12-15 2017-06-23 中国电力科学研究院 一种大规模电磁暂态电网仿真的自动分网方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
岳程燕, 周孝信, 李若梅: "电力系统电磁暂态实时仿真中并行算法的研究", 中国电机工程学报, no. 12, 17 December 2004 (2004-12-17) *

Also Published As

Publication number Publication date
CN110196997B (zh) 2024-04-26

Similar Documents

Publication Publication Date Title
Chenine et al. Modeling and simulation of wide-area communication for centralized PMU-based applications
Dekel et al. Optimal distributed online prediction
Cope Regret and convergence bounds for a class of continuum-armed bandit problems
CN109449970B (zh) 一种适用于高比例分布式光伏配电网的分区方法
CN108122032A (zh) 一种神经网络模型训练方法、装置、芯片和系统
CN110083179A (zh) 一种预定时间多智能体系统一致性跟踪控制方法
CN107204874A (zh) 保证时延最小的sdn网络多控制器部署方法
CN105227601A (zh) 流处理系统中的数据处理方法、装置和系统
CN107069801B (zh) 一种基于最小用户停电损失的配电网孤岛划分方法
CN110198280A (zh) 一种基于bp神经网络的sdn链路分配方法
CN106356840B (zh) 基于同步相量量测的地区电力系统状态估计方法及系统
CN104063281B (zh) 一种调控虚拟机物理内存的方法和装置
Yuan et al. Graph computing based distributed fast decoupled power flow analysis
CN104462657A (zh) 基于Petri网复杂系统多设计方案的选择方法
CN106941683A (zh) 一种lte频间互操作自优化的方法及装置
Li et al. Model-distributed dnn training for memory-constrained edge computing devices
CN110059094A (zh) 主动配电网分布式状态估计的分区方法、介质及设备
Cai et al. A network model for the real-time communications of a smart grid prototype
CN110196997A (zh) 一种电磁暂态异步并行计算方法及装置
CN109857562A (zh) 一种众核处理器上访存距离优化的方法
CN113094887A (zh) 移频电磁暂态仿真的优化方法、装置和电子设备
Xu et al. OMCO: Online multiple coflow scheduling in optical circuit switch
Papadopoulos et al. Measurement based method for online characterization of generator dynamic behaviour in systems with renewable generation
CN110376919A (zh) 一种配电网信息物理系统联合时序仿真方法和系统
CN105656682B (zh) 一种sdn控制器选址方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant