CN104142810A - 一种形成节点导纳矩阵的并行方法 - Google Patents
一种形成节点导纳矩阵的并行方法 Download PDFInfo
- Publication number
- CN104142810A CN104142810A CN201410334960.1A CN201410334960A CN104142810A CN 104142810 A CN104142810 A CN 104142810A CN 201410334960 A CN201410334960 A CN 201410334960A CN 104142810 A CN104142810 A CN 104142810A
- Authority
- CN
- China
- Prior art keywords
- admittance matrix
- nodal admittance
- bus admittance
- network
- parallel method
- 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.)
- Pending
Links
Landscapes
- Traffic Control Systems (AREA)
Abstract
本发明是一种形成节点导纳矩阵的并行方法。包括有如下步骤:1)将形成节点导纳矩阵的任务映射到图形处理器中进行;2)枚举网络所有支路根据支路信息,判断支路类型,在节点导纳矩阵中相应的位置加上相应的值;3)将节点导纳矩阵结果返回到主机中。本发明主要是根据不同支路对节点导纳矩阵的贡献相互独立的原理枚举网络所有支路进行叠加来形成节点导纳矩阵。随着电网的发展和扩大联网,电网的规模越来越大,要形成节点导纳矩阵计算量相对较大,同时节点导纳矩阵在电网潮流计算、短路计算、机电暂态仿真等方面都是必不可少的。本发明利用并行计算来形成节点导纳矩阵可以大大加快计算速度。
Description
技术领域
本发明是一种形成节点导纳矩阵的并行方法,属于电力系统分析的创新领域。
背景技术
电力系统网络模型可用网络元件参数和网络元件连接关系确定。在实际电力系统网络计算中,我们希望有更为简单的网络模型的描述方法,即用一个既包含网络元件参数又包含了元件的连接关系的矩阵来描述电力系统网络模型。节点导纳矩阵和节点阻抗矩阵具有这样的特点,其中节点导纳矩阵由于具有稀疏性而被广泛使用。电网的节点导纳矩阵是电力系统的重要属性之一,不管对电力系统进行什么分析(比如潮流计算,暂态分析等)都离不开节点导纳矩阵的形成。随着电网的发展和扩大联网,电网的规模越来越大,形成节点导纳矩阵的计算量越来越大。近年来,随着对电网运行的要求越来越精细化,对电网的实时分析和在线校核就显得尤为重要,所以对分析电网的速度要求也越来越高。提高节点导纳矩阵的形成速度,必然会提高电力系统分析的速度。所以考虑采用并行方式实现提高形成节点导纳矩阵的速度。
并行计算技术从最早的分布式计算发展到后来的多核CPU,再到近几年发展迅速的GPU,相对比较成熟。GPU具有低成本,高计算效率的优点,所以其在计算领域的地位越来越显著。随着NVIDIA推出的通用并行计算构架CUDA(Compute Unified Device Architecture),使得GPU解决复杂的计算问题能够充分发挥其优势。综上所述,基于CUDA并行计算形成节点导纳矩阵,具有可执行性强、速度快的优势。
发明内容
本发明的目的在于考虑上述问题而提供一种执行简单、速度快的形成节点导纳矩阵的并行方法。本发明比常规的形成节点导纳矩阵的时间要短,具有速度快的优点。
本发明的技术方案是:本发明形成节点导纳矩阵的并行方法,包括有如下步骤:
1)将形成节点导纳矩阵的任务映射到图形处理器中进行;
2)枚举网络所有支路根据支路信息,判断支路类型,在节点导纳矩阵中相应的位置加上相应的值;
3)将节点导纳矩阵结果返回到主机中。
上述步骤1)中,是基于通用并行计算构架将形成节点导纳矩阵的任务映射到图形处理器中进行。
上述步骤2)中各支路的分析是并行的。
本发明形成节点导纳矩阵的并行方法,利用并行计算实现能在速度上更快地形成节点导纳矩阵。
本发明将形成节点导纳矩阵的任务映射到图形处理器中进行,比常规的形成节点导纳矩阵的时间要短,具有速度快的优点。本发明是一种具有高效率、易执行的形成节点导纳矩阵的并行方法。
附图说明
图1为算法流程图;
图2为接地支路示意图;
图3为非标准变比的变压器支路示意图;
图4为普通非接地支路示意图。
具体实施方式
本发明形成节点导纳矩阵的并行方法,包括有如下步骤:
1)基于CUDA(Computer Unified Device Architecture)将形成节点导纳矩阵的任务映射到GPU(Graphic Processing Unit图形处理器)中进行。
2)枚举网络所有支路根据支路信息,判断支路类型,在节点导纳矩阵中相应的位置加上相应的值(各支路的分析是并行的)。
3)将节点导纳矩阵结果返回到主机中。
上述步骤1)中是将形成节点导纳矩阵的任务基于CUDA(ComputerUnified Device Architecture)映射到GPU中进行。
节点导纳矩阵Y的对角线元素Yii称为节点i的自导纳,其值等于接于节点i的所有支路导纳矩阵之和。非对角元素Yij称为节点i、j之间的互导纳,它等于直接联接于节点i、j间的支路导纳的负值。若节点i、j间不存在直接支路,则有Yij=0。根据支路信息(包括开始节点编号、终止节点编号、支路等效导纳)计算节点导纳矩阵,各个支路之间互不干扰,故可用并行计算实现。下面分析网络中不同支路对节点导纳矩阵的贡献。
对于一个网络中有三种类型的支路:
1、接地支路即开始节点编号等于终止节点编号,假设节点编号为i,支路等效导纳值为yi(如图2所示)。则只需把yi加到行号和列号都是i的导纳矩阵元素Yii上,即Yii=Yii+yi。
2、非标准变比的变压器支路,假设其开始节点编号为i,终止节点编号为j,变压器的非标准变比为k,变压器的等效导纳为yt(如图3所示)。则需要在行号为i列号为j和行号为j列号为i的导纳矩阵元素Yij、Yji上减去kyt即、;在行号和列号都为j的导纳矩阵元素上加上k2yt即Yjj=Yjj+k2yt;在行号和列号都为i的导纳矩阵元素上加上yt,即Yii=Yii+yt。
3、普通非接地支路,假设其开始节点编号为i,终止节点编号为j,支路导纳为yij(如图4所示)。则需要在行号为i列号为j和行号为j列号为i的导纳矩阵元素Yij、Yji上减去yij即、;在行号和列号都为j的导纳矩阵元素上加上yij即Yjj=Yjj+yij;在行号和列号都为i的导纳矩阵元素上加上yij,即Yii=Yii+yij。
根据上述原理,把形成节点导纳矩阵的任务映射到GPU中执行需要以下工作:
假设该网络有n条支路,则需要有n个线程,每个线程处理一条支路。在GPU中,各个线程同步进行,先判断该线程所处理支路的类型,按照上述原理修改节点导纳矩阵中的对应元素。直到穷举所有支路,节点导纳矩阵形成。把节点导纳矩阵返回到主机中。
Claims (4)
1.一种形成节点导纳矩阵的并行方法,其特征在于包括有如下步骤:
1)将形成节点导纳矩阵的任务映射到图形处理器中进行;
2)枚举网络所有支路根据支路信息,判断支路类型,在节点导纳矩阵中相应的位置加上相应的值;
3)将节点导纳矩阵结果返回到主机中。
2.根据权利要求1所述的形成节点导纳矩阵的并行方法,其特征在于上述步骤1)中,是基于通用并行计算构架将形成节点导纳矩阵的任务映射到图形处理器中进行。
3.根据权利要求1所述的形成节点导纳矩阵的并行方法,其特征在于上述步骤2)中各支路的分析是并行的。
4.根据权利要求1所述的形成节点导纳矩阵的并行方法,其特征在于利用并行计算实现能在速度上更快地形成节点导纳矩阵。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410334960.1A CN104142810A (zh) | 2014-07-14 | 2014-07-14 | 一种形成节点导纳矩阵的并行方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410334960.1A CN104142810A (zh) | 2014-07-14 | 2014-07-14 | 一种形成节点导纳矩阵的并行方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN104142810A true CN104142810A (zh) | 2014-11-12 |
Family
ID=51851994
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410334960.1A Pending CN104142810A (zh) | 2014-07-14 | 2014-07-14 | 一种形成节点导纳矩阵的并行方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104142810A (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105224290A (zh) * | 2015-09-24 | 2016-01-06 | 清华大学深圳研究生院 | 一种图像负载均衡处理方法及装置 |
CN107741528A (zh) * | 2017-11-17 | 2018-02-27 | 南方电网科学研究院有限责任公司 | 电力系统断线阻抗扫描方法和装置 |
CN108898279A (zh) * | 2018-06-01 | 2018-11-27 | 南方电网科学研究院有限责任公司 | 一种基于私有云的节点导纳计算方法和系统 |
CN112084198A (zh) * | 2020-09-16 | 2020-12-15 | 云南电网有限责任公司 | 一种压缩存储形式的节点导纳矩阵直接列写方法及装置 |
CN112615377A (zh) * | 2020-12-29 | 2021-04-06 | 南方电网科学研究院有限责任公司 | 电压支撑强度指标批量计算方法、装置及存储介质 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1601848A (zh) * | 2003-09-28 | 2005-03-30 | 中国电力科学研究院 | 电力系统数字仿真方法 |
US20100305888A1 (en) * | 2007-05-07 | 2010-12-02 | Siemens Aktiengesellschaft | Method and Device for Determining Load Flow in an Electrical Power Supply System |
CN103617150A (zh) * | 2013-11-19 | 2014-03-05 | 国家电网公司 | 一种基于gpu的大规模电力系统潮流并行计算系统及其方法 |
-
2014
- 2014-07-14 CN CN201410334960.1A patent/CN104142810A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1601848A (zh) * | 2003-09-28 | 2005-03-30 | 中国电力科学研究院 | 电力系统数字仿真方法 |
US20100305888A1 (en) * | 2007-05-07 | 2010-12-02 | Siemens Aktiengesellschaft | Method and Device for Determining Load Flow in an Electrical Power Supply System |
CN103617150A (zh) * | 2013-11-19 | 2014-03-05 | 国家电网公司 | 一种基于gpu的大规模电力系统潮流并行计算系统及其方法 |
Non-Patent Citations (1)
Title |
---|
郭春辉: ""基于GPU的电力系统并行计算的研究"", 《中国优秀硕士学位论文全文数据库 工程科技II辑》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105224290A (zh) * | 2015-09-24 | 2016-01-06 | 清华大学深圳研究生院 | 一种图像负载均衡处理方法及装置 |
CN105224290B (zh) * | 2015-09-24 | 2017-11-07 | 清华大学深圳研究生院 | 一种图像负载均衡处理方法及装置 |
CN107741528A (zh) * | 2017-11-17 | 2018-02-27 | 南方电网科学研究院有限责任公司 | 电力系统断线阻抗扫描方法和装置 |
CN107741528B (zh) * | 2017-11-17 | 2019-12-27 | 南方电网科学研究院有限责任公司 | 电力系统断线阻抗扫描方法和装置 |
CN108898279A (zh) * | 2018-06-01 | 2018-11-27 | 南方电网科学研究院有限责任公司 | 一种基于私有云的节点导纳计算方法和系统 |
CN112084198A (zh) * | 2020-09-16 | 2020-12-15 | 云南电网有限责任公司 | 一种压缩存储形式的节点导纳矩阵直接列写方法及装置 |
CN112084198B (zh) * | 2020-09-16 | 2022-09-16 | 云南电网有限责任公司 | 一种压缩存储形式的节点导纳矩阵直接列写方法及装置 |
CN112615377A (zh) * | 2020-12-29 | 2021-04-06 | 南方电网科学研究院有限责任公司 | 电压支撑强度指标批量计算方法、装置及存储介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109740747B (zh) | 运算方法、装置及相关产品 | |
CN104142810A (zh) | 一种形成节点导纳矩阵的并行方法 | |
CN103617150A (zh) | 一种基于gpu的大规模电力系统潮流并行计算系统及其方法 | |
CN106250349A (zh) | 一种高能效异构计算系统 | |
CN103106183A (zh) | 基于mapreduce的大规模稀疏矩阵乘法运算的方法 | |
CN104570081A (zh) | 一种积分法叠前时间偏移地震资料处理方法及系统 | |
CN103679545A (zh) | 一种配电网可靠性理论评估系统 | |
CN103279330A (zh) | 一种基于虚拟机GPU计算下的MapReduce并行编程模型 | |
CN109066687A (zh) | 一种电力系统潮流计算方法、系统及电子设备 | |
CN108958852A (zh) | 一种基于fpga异构平台的系统优化方法 | |
CN104967121B (zh) | 一种大规模电力系统节点的潮流计算方法 | |
CN103294599A (zh) | 一种基于云的嵌入式软件交叉测试方法 | |
CN105389666A (zh) | 配网模型拓扑连通性校验方法 | |
CN106528490A (zh) | 一种fpga异构加速计算装置及系统 | |
CN103279824A (zh) | 一种继电保护整定计算系统的建模方法 | |
CN106168959A (zh) | 网页布局方法及装置 | |
Boghetti et al. | Building energy models with Morphological urban-scale parameters: A case study in Turin | |
CN104615437A (zh) | 一种基于gpu的软件体系架构及uml和adl结合描述的方法 | |
CN103885842B (zh) | 一种带加速节点的片上网络的优化整体的任务映射方法 | |
US20140143746A1 (en) | Direct current circuit analysis based clock network design | |
Antunes et al. | Partitioning and mapping on NoC-based MPSoC: an energy consumption saving approach | |
CN104636311A (zh) | 多处理器片上网络互连阵列尺寸的开发方法 | |
Zhao et al. | GPU based parallel matrix exponential algorithm for large scale power system electromagnetic transient simulation | |
CN109165426A (zh) | 电磁暂态仿真方法和电磁暂态仿真装置 | |
CN105205854A (zh) | 一种数据中心3d建模方法和装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20141112 |
|
RJ01 | Rejection of invention patent application after publication |