CN103186366A - 基于cuda并行计算实现电力***电磁暂态实时仿真测试方法 - Google Patents
基于cuda并行计算实现电力***电磁暂态实时仿真测试方法 Download PDFInfo
- Publication number
- CN103186366A CN103186366A CN2012102414106A CN201210241410A CN103186366A CN 103186366 A CN103186366 A CN 103186366A CN 2012102414106 A CN2012102414106 A CN 2012102414106A CN 201210241410 A CN201210241410 A CN 201210241410A CN 103186366 A CN103186366 A CN 103186366A
- Authority
- CN
- China
- Prior art keywords
- cuda
- thread
- simulation
- electromagnetic transient
- grid
- 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
Links
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种基于CUDA并行计算实现电力***电磁暂态实时仿真测试方法,电磁暂态仿真测试主机采用CPU+GPU方式,GPU应用CUDA平台并行计算架构,仿真测试主机通过千兆以太网和各路信号采集设备连接;主要针对数字化变电站的仿真而不是全网或全***的仿真,模型可以大大简化,能够将传统仿真速度提高几十至几百倍,实现变电站电磁暂态实时仿真计算要求;由于CUDA并行计算平台的成本相对于大型并行计算机或专业并行计算机的成本大大降低,本发明提出的方法在电力***电磁暂态实时仿真上性能大幅提高、成本巨幅降低,具有广泛推广的价值。
Description
技术领域
本发明涉及电力控制领域,具体涉及基于CUDA并行计算实现电力***电磁暂态实时仿真测试方法。
背景技术
电力***电磁暂态过程和机电暂态过程,是两个用不同数学模型表征、具有不同时间常数的物理过程。在传统的电力***分析工具中,通常对这两个过程分别进行数字仿真。其中,相对于机电暂态仿真,电磁暂态模型的非线性和参量变化的瞬时性,大大增加了仿真的计算量和计算时间,实现一定规模***的实时数字仿真要困难得多。随着区域电网的互联和直流输电的发展以及FACTS、SVC等电力电子装备在电力***中的大量应用,使得电网在规模上不断扩展,元件构成上日趋复杂,不同物理特性的动态过程-电磁暂态和机电暂态过程交织并行,成为现代电力***的主要特征,使得电力***分析和仿真技术面临新的挑战。如何将大规模复杂电力***的机电暂态仿真和对于局部***的电磁暂态仿真集成在一个进程中,并进一步使数字模拟和实际发生的物理过程同步,实现大规模电力***的电磁暂态和机电暂态混合实时仿真是近年来电力***数字仿真研究领域的前端课题。这一问题的解决,对分析研究电力***动态特性,改善电力***控制和保护水平,提高电力***运行的稳定性和可靠性,有着重要的理论价值和现实意义。
目前电力***的电磁暂态实时仿真平台均由大型计算机组成的并行计算仿真平台,动辄数百万甚至上千万人民币。比如,世界上应用最广的电力***实时数字仿真仪是由加拿大曼尼托巴RTDS公司开发制造的RTDS,是一种专门 设计用于研究电力***中电磁暂态现象的装置。但由于其昂贵的售价,让用户对其望而却步。
发明内容
本发明提供了一种处理速度快,造价低廉的实现电力***电磁暂态实时仿真方法。
为解决上述技术问题,本发明通过以下技术方案来实现:基于数字化变电站电磁暂态实时仿真测试***,其实现包括以下步骤:
a、电磁暂态仿真测试主机采用CPU+GPU方式,GPU应用CUDA平台并行计算架构,仿真测试主机通过千兆以太网和各路信号采集设备连接;
b、所述CUDA把计算分解成一个一个的线程(kernel):kernel_1,kernel_2……kernel_M,每一个线程(kernel)就交给每一个网格(Grid)来完成,网格(Grid)又把任务划分为一个个的块(Block),这里每一个网格(Grid)管理的块(Block)是有限的,块(Block)这个层直接管理每一个线程(Thread);
c、在进行电力***电磁暂态实时仿真的时候,就是将2k次的计算分配到不同的线程(Thread)中去,由于CUDA计算核的数量非常多,当2k小于这个数量时,在每一个2k的计算中可以将计算分解成更多的并行计算,通过开发适合CUDA程序执行模型的算法,并合理设置块(Block)和线程(Thread)数目,可充分利用GPU强大的并行处理能力,解决电磁暂态海量数据的处理问题。
d、具体计算方法采用线性常微分方程组,经过简单变换,可以写成如下形式:
dX/dt=F(X)+G(U) (1)
X为***状态变量矩阵
采用两点欧拉法(其它方法原理相同),即
Xi n+1=Xi n-1+(2h)(dXi n/dt) (2)
将(1)代入(2),得
Xi n+1=Xi n-1+(2h)(Fi(Xn)+Gi(Un)) (3)
为保证计算方法的数值稳定性,用预测-校正的方法解方程组(3),得到(4)式(这里,Xip n+1预测值,Xic n+1为校正值)
Xip n+1=Xi n-1+2h(Fi(Xn)+Gi(Un)) (4)
Xic n+1=Xi n+(h/2)(Fi(Xn)+Gi(Un)+Fi(Xp n+1)+Gi(Un+1)) (5)
其中,i=1,2,…,k
(4)(5)包含2k个等式计算,但(4)式k个等式必须先计算,(5)式因需用到Xp n+1,Un+1,所以必须等(4)式全部计算完才能计算;
为实现2k个等式同时计算,将Xp n+1,Un+1用Xp n,Un代替,得到下式:
Xic n+1=Xi n+(h/2)(Fi(Xn)+Gi(Un)+Fi(Xp n)+Gi(Un)) (6)
联合(4)(6)即可实现2k个等式的并行计算。
本发明具有以下优点:1、主要针对数字化变电站的仿真而不是全网或全***的仿真,模型可以大大简化,能够将传统仿真速度提高几十至几百倍,实现变电站电磁暂态实时仿真计算要求;2、由于CUDA并行计算平台的成本相对于大型并行计算机或专业并行计算机的成本大大降低,本专利提出的方法在电力***电磁暂态实时仿真上性能大幅提高、成本巨幅降低,具有广泛推广的价值。
附图说明
一下结合附图对本发明做详细的举例说明;
图1为本发明CUDA并行计算计算机实现框图;
图2为本发明CUDA并行计算实现流程框图。
具体实施方式
实施例一
图1所示,CUDA把计算分解成一个一个的线程kernel:kernel_1,kernel_2……kernel_M,每一个线程kernel就交给每一个网格Grid来完成,网格Grid又把任务划分为一个个的块Block,这里每一个网格Grid管理的块Block是有限的,Block这个层直接管理每一个线程T hread,在进行电力***电磁暂态实时仿真的时候,就是将2k次的计算分配到不同的线程Thread中去,由于CUDA计算核的数量非常多,当2k小于这个数量时,在每一个2k的计算中可以将计算分解成更多的并行计算,通过开发适合CUDA程序执行模型的算法,并合理设置Block和Thread数目,可充分利用GPU强大的并行处理能力,解决电磁暂态海量数据的处理问题。
实施例二
图2所示,电力***电磁暂态仿真要考虑发电机定子绕组及线路的电磁暂态过程,需全部用微分方程描述。这是些线性常微分方程组,经过简单变换,可以写成如下形式:
dX/dt=F(X)+G(U) (1)
X为***状态变量矩阵
采用两点欧拉法(其它方法原理相同),即
Xi n+1=Xi n-1+(2h)(dXi n/dt) (2)
将(1)代入(2),得
Xi n+1=Xi n-1+(2h)(Fi(Xn)+Gi(Un)) (3)
为保证计算方法的数值稳定性,用预测-校正的方法解方程组(3),得到(4)式。这里,Xip n+1预测值,Xic n+1为校正值
Xip n+1=Xi n-1+2h(Fi(Xn)+Gi(Un)) (4)
Xic n+1=Xi n+(h/2)(Fi(Xn)+Gi(Un)+Fi(Xp n+1)+Gi(Un+1)) (5)
其中,i=1,2,…,k
(4)(5)包含2k个等式计算,但(4)式k个等式必须先计算,(5)式因需用到Xp n+1,Un+1,所以必须等(4)式全部计算完才能计算。
为实现2k个等式同时计算,将Xp n+1,Un+1用Xp n,Un代替,得到下式:
Xic n+1=Xi n+(h/2)(Fi(Xn)+Gi(Un)+Fi(Xp n)+Gi(Un)) (6)
联合(4)(6)即可实现2k个等式的并行计算。
Claims (1)
1.基于CUDA并行计算实现电力***电磁暂态实时仿真测试方法,基于数字化变电站电磁暂态实时仿真测试***,其实现包括以下步骤:
a、电磁暂态仿真测试主机采用CPU+GPU方式,GPU应用CUDA平台并行计算架构,仿真测试主机通过千兆以太网和各路信号采集设备连接;
b、所述CUDA把计算分解成一个一个的线程(kernel):kernel_1,kernel_2……kernel_M,每一个线程(kernel)就交给每一个网格(Grid)来完成,网格(Grid)又把任务划分为一个个的块(Block),这里每一个网格(Grid)管理的块(Block)是有限的,块(Block)这个层直接管理每一个线程(Thread);
c、在进行电力***电磁暂态实时仿真的时候,就是将2k次的计算分配到不同的线程(Thread)中去,由于CUDA计算核的数量非常多,当2k小于这个数量时,在每一个2k的计算中可以将计算分解成更多的并行计算,通过开发适合CUDA程序执行模型的算法,并合理设置块(Block)和线程(Thread)数目,可充分利用GPU强大的并行处理能力,解决电磁暂态海量数据的处理问题。
d、具体计算方法采用线性常微分方程组,经过简单变换,可以写成如下形式:
dX/dt=F(X)+G(U) (1)
X为***状态变量矩阵
采用两点欧拉法(其它方法原理相同),即
Xi n+1=Xi n-1+(2h)(dXi n/dt) (2)
将(1)代入(2),得
Xi n+1=Xi n-1+(2h)(Fi(Xn)+Gi(Un)) (3)
为保证计算方法的数值稳定性,用预测-校正的方法解方程组(3),得到(4)式(这里,Xip n+1预测值,Xic n+1为校正值)
Xip n+1=Xi n-1+2h(Fi(Xn)+Gi(Un)) (4)
Xic n+1=Xi n+(h/2)(Fi(Xn)+Gi(Un)+Fi(Xp n+1)+Gi(Un+1)) (5)
其中,i=1,2,…,k
(4)(5)包含2k个等式计算,但(4)式k个等式必须先计算,(5)式因需用到Xp n+1,Un+1,所以必须等(4)式全部计算完才能计算;
为实现2k个等式同时计算,将Xp n+1,Un+1用Xp n,Un代替,得到下式:
Xic n+1=Xi n+(h/2)(Fi(Xn)+Gi(Un)+Fi(Xp n)+Gi(Un)) (6)
联合(4)(6)即可实现2k个等式的并行计算。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210241410.6A CN103186366B (zh) | 2012-07-12 | 2012-07-12 | 基于cuda并行计算实现电力***电磁暂态实时仿真测试方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210241410.6A CN103186366B (zh) | 2012-07-12 | 2012-07-12 | 基于cuda并行计算实现电力***电磁暂态实时仿真测试方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103186366A true CN103186366A (zh) | 2013-07-03 |
CN103186366B CN103186366B (zh) | 2017-05-03 |
Family
ID=48677546
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210241410.6A Expired - Fee Related CN103186366B (zh) | 2012-07-12 | 2012-07-12 | 基于cuda并行计算实现电力***电磁暂态实时仿真测试方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103186366B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105574809A (zh) * | 2015-12-16 | 2016-05-11 | 天津大学 | 基于矩阵指数的电磁暂态仿真图形处理器并行计算方法 |
CN105760599A (zh) * | 2016-02-16 | 2016-07-13 | 东南大学 | 一种用于大型电力***暂态稳定研究的小型化超级仿真器 |
CN106200413A (zh) * | 2015-05-05 | 2016-12-07 | 中国电力科学研究院 | 基于Aurora协议的电磁暂态实时仿真通讯方法和装置 |
CN109558661A (zh) * | 2018-11-22 | 2019-04-02 | 北京未尔锐创科技有限公司 | 一种电磁环境计算方法及装置 |
CN114117872A (zh) * | 2022-01-24 | 2022-03-01 | 广州中望龙腾软件股份有限公司 | 一种多gpu并行的时域有限差分电磁仿真方法、设备及介质 |
CN116521339A (zh) * | 2023-04-26 | 2023-08-01 | 中国人民解放军92942部队 | 一种面向多用户并发访问的信号级仿真实现方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101382969A (zh) * | 2008-10-31 | 2009-03-11 | 中国电力科学研究院 | 一种多步变步长电磁暂态仿真方法 |
US20100135418A1 (en) * | 2008-11-28 | 2010-06-03 | Thomson Licensing | Method for video decoding supported by graphics processing unit |
CN102403720A (zh) * | 2011-11-23 | 2012-04-04 | 昆明理工大学 | 一种基于暂态电压安全裕度的超实时重合时序整定方法 |
-
2012
- 2012-07-12 CN CN201210241410.6A patent/CN103186366B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101382969A (zh) * | 2008-10-31 | 2009-03-11 | 中国电力科学研究院 | 一种多步变步长电磁暂态仿真方法 |
US20100135418A1 (en) * | 2008-11-28 | 2010-06-03 | Thomson Licensing | Method for video decoding supported by graphics processing unit |
CN102403720A (zh) * | 2011-11-23 | 2012-04-04 | 昆明理工大学 | 一种基于暂态电压安全裕度的超实时重合时序整定方法 |
Non-Patent Citations (2)
Title |
---|
XUE-XIN LIU等: "Transient Analysis of Large Linear Dynamic Networks on Hybid GPU-Multicore Platforms", 《NEW CIRCUITS AND SYSTEMS CONFERENCE (NEWCAS), 2012 IEEE 10TH INTERNATIONAL》 * |
王蕾等: "电力***暂态稳定的两种仿真算法分析", 《计算技术与自动化》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106200413A (zh) * | 2015-05-05 | 2016-12-07 | 中国电力科学研究院 | 基于Aurora协议的电磁暂态实时仿真通讯方法和装置 |
CN106200413B (zh) * | 2015-05-05 | 2019-08-06 | 中国电力科学研究院 | 基于Aurora协议的电磁暂态实时仿真通讯方法和装置 |
CN105574809A (zh) * | 2015-12-16 | 2016-05-11 | 天津大学 | 基于矩阵指数的电磁暂态仿真图形处理器并行计算方法 |
CN105574809B (zh) * | 2015-12-16 | 2018-07-20 | 天津大学 | 基于矩阵指数的电磁暂态仿真图形处理器并行计算方法 |
CN105760599A (zh) * | 2016-02-16 | 2016-07-13 | 东南大学 | 一种用于大型电力***暂态稳定研究的小型化超级仿真器 |
CN105760599B (zh) * | 2016-02-16 | 2019-02-22 | 东南大学 | 一种用于大型电力***暂态稳定研究的小型化超级仿真器 |
CN109558661A (zh) * | 2018-11-22 | 2019-04-02 | 北京未尔锐创科技有限公司 | 一种电磁环境计算方法及装置 |
CN114117872A (zh) * | 2022-01-24 | 2022-03-01 | 广州中望龙腾软件股份有限公司 | 一种多gpu并行的时域有限差分电磁仿真方法、设备及介质 |
CN114117872B (zh) * | 2022-01-24 | 2022-06-14 | 广州中望龙腾软件股份有限公司 | 一种多gpu并行的时域有限差分电磁仿真方法、设备及介质 |
CN116521339A (zh) * | 2023-04-26 | 2023-08-01 | 中国人民解放军92942部队 | 一种面向多用户并发访问的信号级仿真实现方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103186366B (zh) | 2017-05-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Li et al. | Simulation of the smart grid communications: Challenges, techniques, and future trends | |
CN103186366A (zh) | 基于cuda并行计算实现电力***电磁暂态实时仿真测试方法 | |
CN103793562B (zh) | 基于fpga的有源配电网暂态实时仿真***设计方法 | |
Wang et al. | Reliability analysis of phasor measurement unit considering data uncertainty | |
CN106063070B (zh) | 用于在电力网络中注入或提取功率的方法和*** | |
CN103106314B (zh) | 太阳能光伏电源输出功率的时序概率建模方法 | |
Vanfretti et al. | SmarTS Lab—A laboratory for developing applications for WAMPAC Systems | |
CN103746370A (zh) | 一种风电场可靠性建模方法 | |
Jalili-Marandi et al. | A real-time dynamic simulation tool for transmission and distribution power systems | |
CN104104081B (zh) | 一种基于优化方法的非迭代不确定潮流分析方法 | |
CN107423476B (zh) | 基于多fpga的有源配电网实时仿真器并行通讯方法 | |
CN107862103A (zh) | 一种基于adpss的仿真模型生成方法及*** | |
CN105005652B (zh) | 一种基于fpga的小步长实时仿真*** | |
CN103473602A (zh) | 一种电网理论线损数据预测***及预测方法 | |
CN105207197A (zh) | 包含风电场的电力***可靠性评估方法 | |
CN104252550A (zh) | 一种基于fpga的高性能实时仿真方法 | |
CN103605293B (zh) | 一种变电站继电保护仿真测试***及其仿真测试方法 | |
CN106651144A (zh) | 一种带有同杆并架线路信息的杆塔与线路关系建模方法 | |
Cao et al. | Faster-than-real-time dynamic simulation of AC/DC grids on reconfigurable hardware | |
CN104865842A (zh) | 适用于稳定控制装置接口的混合仿真***及混合仿真方法 | |
CN106056305A (zh) | 一种基于状态聚类的发电***可靠性快速评估方法 | |
CN107122562B (zh) | 基于多fpga的有源配电网实时仿真器串行通讯方法 | |
Bai et al. | Automatic modeling and optimization for the digital twin of a regional multi-energy system | |
CN104779613B (zh) | 基于试验的含变流器电力元件等效建模方法 | |
CN206892220U (zh) | 一种电能质量监测装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170503 Termination date: 20210712 |
|
CF01 | Termination of patent right due to non-payment of annual fee |