CN106294897B - 一种适用于电磁暂态多时间尺度实时仿真接口的实现方法 - Google Patents
一种适用于电磁暂态多时间尺度实时仿真接口的实现方法 Download PDFInfo
- Publication number
- CN106294897B CN106294897B CN201510264005.XA CN201510264005A CN106294897B CN 106294897 B CN106294897 B CN 106294897B CN 201510264005 A CN201510264005 A CN 201510264005A CN 106294897 B CN106294897 B CN 106294897B
- Authority
- CN
- China
- Prior art keywords
- time
- fpga
- time scale
- platform
- real
- 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.)
- Active
Links
- 238000004088 simulation Methods 0.000 title claims abstract description 58
- 238000000034 method Methods 0.000 title claims abstract description 24
- 230000001052 transient effect Effects 0.000 title claims abstract description 18
- 238000004364 calculation method Methods 0.000 claims abstract description 25
- 239000011159 matrix material Substances 0.000 claims abstract description 16
- 230000005540 biological transmission Effects 0.000 claims abstract description 5
- 238000007781 pre-processing Methods 0.000 claims abstract description 4
- 238000013213 extrapolation Methods 0.000 claims description 3
- 230000008859 change Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 2
- 101000579490 Solanum lycopersicum Suberization-associated anionic peroxidase 1 Proteins 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000006855 networking Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
Abstract
本发明涉及一种适用于电磁暂态多时间尺度实时仿真接口的实现方法,所述方法包括下述步骤:根据电磁暂态多时间尺度,将仿真系统划分为多个通过解耦元件传输线连接的子网络,并确定所述子网络的时间尺度;通过戴维南等值电路将解耦元件分解为两个内阻固定的受控源,并入所述子网络中,进行网络矩阵初始化;进行仿真计算,接收FPGA平台的数据;预处理FPGA平台的数据,获得中间变量,继续仿真计算。
Description
技术领域
本发明涉及一种实现方法,具体涉及一种适用于电磁暂态多时间尺度实时仿真接口的实现方法。
背景技术
电力系统中,为了验证新的电力自动化设备的功能和性能,在设备投入实际系统运行前需要进行大量测试。电力系统实时数字仿真系统可以实时模拟电力系统各种运行工况,具有体积小、功耗低、通用性好、可重复性强、价格较动态模拟和数模混合式仿真装置低廉等优点,因而在电力自动化设备的测试中获得了广泛应用。
小步长仿真系统是一个完整的电磁暂态仿真系统的重要组成部分。随着仿真步长变小,系统的仿真负担也随之增加,从而需要使用FPGA硬件加速平台用以实现小步长仿真系统的实时的双精度浮点运算。
但是由于仿真系统的仿真能力的限制,很多大型的仿真算例并不能完全在一个平台上运行。因此必须把大型的仿真算例划分为规模较小的网络,并通过接口互联,使得仿真结果与一个完整的大型仿真算例完全一致,而FPGA上只进行小步长仿真的局部网络。于是产生了在FPGA平台上进行多时间尺度分网并行仿真的需求。
FPGA平台的设计必须保证三个原则:
1.FPGA的程序应该实现保证并行,减少数据间的依赖关系。
2.FPGA为了提升仿真能力,应该缩短单一计算流水线的长度。
3.FPGA的LUT既可以用作存储又可以用于计算;为了提升FPGA的计算能力,必须减少FPGA的存储空间的使用。
现有的多时间尺度分网算法的设计方法利用了传输线自然延迟特性。其中,等值计算电路如图2所示,
电流源递推公式为:
由此递推公式可知,每一个子网上都需要保留2τmi/dt的变量,同时,由于τmi是一个变量,导致变量缓存区的大小不能固定;同时,上述递推公式中变量是模态量,交流的模态量转换为仿真系统需要的瞬时量,该转换是一个长流水线计算,会占用大量时序;最后,为了提升FPGA平台的仿真能力,必须尽量优化FPGA侧的计算,对于复杂的接口计算必须进行简化。
发明内容
为了实现上述目的,本发明提供一种适用于电磁暂态多时间尺度实时仿真接口的实现方法,有效的解决电力系统的FPGA小步长仿真系统与服务器系统之间的互联问题,解决了仿真规模限制的问题,极大提升了仿真效率。
本发明的目的是采用下述技术方案实现的:
一种适用于电磁暂态多时间尺度实时仿真接口的实现方法,所述方法包括下述步骤:
(1)根据电磁暂态多时间尺度,将仿真系统划分为多个通过解耦元件传输线连接的子网络,并确定所述子网络的时间尺度;
(2)通过戴维南等值电路将解耦元件分解为两个内阻固定的受控源,并入所述子网络中,进行网络矩阵初始化;
(3)进行仿真计算,接收FPGA平台的数据;
(4)预处理FPGA平台的数据,获得中间变量,继续仿真计算。
优选的,所述步骤(1)中,根据系统的动态时间常数确定所述子网络的时间尺度包括最小时间尺度和非最小时间尺度。
进一步地,所述最小时间尺度的网络在实时FPGA仿真平台上运行,非最小时间尺度的网络在实时服务器平台上运行。
优选的,所述步骤(2)中网络矩阵初始化包括:采用外插值法对电磁暂态仿真系统的初始电压和电流进行反推,估算当前时刻之前的系统电压和电流,并记录该时间与解耦元件相关的电压和电流;
完成网络矩阵初始化后,将该网络矩阵存储至FPGA平台的内存中。
优选的,所述步骤(3)中,通过FPGA平台进行时序控制,同步启动所述FPGA平台与实时服务器平台,每隔一个非最小时间尺度向所述实时服务器平台传输在FPGA系统上运行的网络矩阵的数据。
优选的,所述步骤(4)中,当实时服务器平台接收到FPGA传输的数据后,与实时服务器平台上的历史接口数据进行比较,获得FPGA平台下一个非最小时间尺度的中间变量,并在一个非最小时间尺度内传送至FPGA平台。
进一步地,在一个非最小时间尺度内,如果FPGA接收到实时服务器平台的数据,则重复步骤(3);如未收到,则利用所述下一个非最小时间尺度的中间变量进行仿真计算。
与最接近的现有技术比,本发明达到的有益效果是:
1.本发明提供了电磁暂态多时间尺度实时仿真接口的实现方法,在具体的仿真计算前,通过初始化把解耦元件分解为两个内阻固定的受控源,通过在服务器侧的预计算,将受控源在仿真开始前并入了子网络,防止由于受控源计算的变化使核心网络方程造成影响。由此固化了计算过程,有利于FPGA的实现。
2.提出了解耦元件预计算的解决方案,结合了FPGA运算的特点,把可能发生的情况都在预计算过程中进行处理,即将变量变化导致计算和存储变化的计算区域提取出来,放入了独立的服务器芯片中进行计算,把必须保留的固定的计算部分保留在FPGA中实现。大大节省了FPGA的资源,提升了FPGA的网络计算能力。
3.提出了对接口仿真计算的拆分和设计,解决了FPGA在多尺度并行分网计算中接口计算负荷的不平衡、占用过多计算资源的问题,提高了计算速度和仿真规模。
附图说明
图1为本发明提供的电磁暂态多时间尺度实时仿真接口的实现方法流程图;
图2为背景技术提供的等值计算电路结构示意图。
具体实施方式
如图1所示,一种适用于电磁暂态多时间尺度实时仿真接口的实现方法,所述方法包括下述步骤:
(1)根据电磁暂态多时间尺度,将仿真系统划分为多个通过解耦元件传输线连接的子网络,并确定所述子网络的时间尺度;
所述步骤(1)中,根据系统的动态时间常数确定所述子网络的时间尺度包括最小时间尺度和非最小时间尺度。
根据不同的仿真精度,可选择以下动态时间常数(微秒):1,10,50,100,1000。
所述最小时间尺度的网络在实时FPGA仿真平台上运行,非最小时间尺度的网络在实时服务器平台上运行。
(2)通过戴维南等值电路将解耦元件分解为两个内阻固定的受控源,并入所述子网络中,进行网络矩阵初始化;
所述步骤(2)中,网络矩阵初始化包括:采用外插值法对电磁暂态仿真系统的初始电压和电流进行反推,估算当前时刻之前的系统电压和电流,并记录该时间与解耦元件相关的电压和电流;
完成网络矩阵初始化后,将该网络矩阵存储至FPGA平台的内存(SRAM)中。
(3)进行仿真计算,接收FPGA平台的数据;
所述步骤(3)中,通过FPGA平台进行时序控制,同步启动所述FPGA平台与实时服务器平台,每隔一个非最小时间尺度向所述实时服务器平台传输在FPGA系统上运行的网络矩阵的数据。
(4)预处理FPGA平台的数据,获得中间变量(SUB4_LC_TMP1和SUB4_LC_TMP5),继续仿真计算。
所述步骤(4)中,当实时服务器平台接收到FPGA传输的数据后,与实时服务器平台上的历史接口数据进行比较,获得FPGA平台下一个非最小时间尺度的中间变量,并在一个非最小时间尺度内传送至FPGA平台。
在一个非最小时间尺度内,如果FPGA接收到实时服务器平台的数据,则重复步骤(3);如未收到,则利用所述下一个非最小时间尺度的中间变量进行仿真计算。
最后应当说明的是:以上实施例仅用以说明本发明的技术方案而非对其限制,尽管参照上述实施例对本发明进行了详细的说明,所属领域的普通技术人员应当理解:依然可以对本发明的具体实施方式进行修改或者等同替换,而未脱离本发明精神和范围的任何修改或者等同替换,其均应涵盖在本发明的权利要求范围当中。
Claims (5)
1.一种适用于电磁暂态多时间尺度实时仿真接口的实现方法,所述方法包括下述步骤:
(1)根据电磁暂态多时间尺度,将仿真系统划分为多个通过解耦元件传输线连接的子网络,并确定所述子网络的时间尺度;
(2)通过戴维南等值电路将解耦元件分解为两个内阻固定的受控源,并入所述子网络中,进行网络矩阵初始化;
(3)进行仿真计算,接收FPGA平台的数据;
(4)预处理FPGA平台的数据,获得中间变量,继续仿真计算;
所述步骤(2)中网络矩阵初始化包括:采用外插值法对电磁暂态仿真系统的初始电压和电流进行反推,估算当前时刻之前的系统电压和电流,并记录该时间与解耦元件相关的电压和电流;
完成网络矩阵初始化后,将该网络矩阵存储至FPGA平台的内存中;
所述步骤(4)中,当实时服务器平台接收到FPGA传输的数据后,与实时服务器平台上的历史接口数据进行比较,获得FPGA平台下一个非最小时间尺度的中间变量,并在一个非最小时间尺度内传送至FPGA平台。
2.如权利要求1所述的电磁暂态多时间尺度实时仿真接口的实现方法,其特征在于,所述步骤(1)中,根据系统的动态时间常数确定所述子网络的时间尺度包括最小时间尺度和非最小时间尺度。
3.如权利要求2所述的电磁暂态多时间尺度实时仿真接口的实现方法,其特征在于,所述最小时间尺度的网络在实时FPGA仿真平台上运行,非最小时间尺度的网络在实时服务器平台上运行。
4.如权利要求1所述的电磁暂态多时间尺度实时仿真接口的实现方法,其特征在于,所述步骤(3)中,通过FPGA平台进行时序控制,同步启动所述FPGA平台与实时服务器平台,每隔一个非最小时间尺度向所述实时服务器平台传输在FPGA系统上运行的网络矩阵的数据。
5.如权利要求1所述的电磁暂态多时间尺度实时仿真接口的实现方法,其特征在于,在一个非最小时间尺度内,如果FPGA接收到实时服务器平台的数据,则重复步骤(3);如未收到,则利用所述下一个非最小时间尺度的中间变量进行仿真计算。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510264005.XA CN106294897B (zh) | 2015-05-22 | 2015-05-22 | 一种适用于电磁暂态多时间尺度实时仿真接口的实现方法 |
PCT/CN2016/089823 WO2016188503A2 (zh) | 2015-05-22 | 2016-07-12 | 一种适用于电磁暂态多时间尺度实时仿真接口的实现方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510264005.XA CN106294897B (zh) | 2015-05-22 | 2015-05-22 | 一种适用于电磁暂态多时间尺度实时仿真接口的实现方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106294897A CN106294897A (zh) | 2017-01-04 |
CN106294897B true CN106294897B (zh) | 2020-03-17 |
Family
ID=57393517
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510264005.XA Active CN106294897B (zh) | 2015-05-22 | 2015-05-22 | 一种适用于电磁暂态多时间尺度实时仿真接口的实现方法 |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN106294897B (zh) |
WO (1) | WO2016188503A2 (zh) |
Families Citing this family (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106777827B (zh) * | 2017-01-24 | 2023-08-15 | 中国电力科学研究院 | 一种机电-电磁混合仿真方法及系统 |
CN108628182A (zh) * | 2017-03-24 | 2018-10-09 | 中国电力科学研究院 | 一种电磁暂态小步长仿真方法及系统 |
CN108984847B (zh) * | 2018-06-21 | 2022-12-06 | 武汉大学 | 一种基于分频阻抗补偿的实时数字混合仿真接口方法 |
CN109002660B (zh) * | 2018-09-07 | 2022-12-06 | 天津大学 | 基于fpga的有源配电网实时仿真解算器通用化设计方法 |
CN109508479B (zh) * | 2018-10-19 | 2022-12-16 | 天津大学 | 基于fpga的有源配电网实时仿真器参数配置通用化方法 |
CN109190320A (zh) * | 2018-11-02 | 2019-01-11 | 贵州电网有限责任公司 | 适用于adpss双馈风机闭环试验的并行异构仿真方法 |
CN109783845B (zh) * | 2018-12-05 | 2023-07-18 | 中国电力科学研究院有限公司 | 边界点分群解耦的机电-电磁混合仿真分网方法及装置 |
CN109711001B (zh) * | 2018-12-10 | 2023-03-24 | 三峡大学 | 基于柔性化子系统划分和场景属性聚类的水电站多维全数字建模方法 |
CN109768542B (zh) * | 2018-12-14 | 2022-10-11 | 中国电力科学研究院有限公司 | 一种大规模电网电磁暂态自动建模的继承方法及系统 |
CN111695258B (zh) * | 2020-06-11 | 2023-05-26 | 哈尔滨工业大学 | 一种电磁继电器动态特性仿真动能注入仿真方法 |
CN112084624A (zh) * | 2020-07-31 | 2020-12-15 | 清华大学 | 电磁-机电混合仿真电磁暂态侧接口功率计算方法与装置 |
CN113033024B (zh) * | 2021-04-25 | 2022-09-02 | 海南省电力学校(海南省电力技工学校) | 输电网细粒度并行电磁暂态仿真方法、系统、终端及介质 |
CN113887052A (zh) * | 2021-10-09 | 2022-01-04 | 清鸾科技(成都)有限公司 | 综合能源系统及其仿真计算方法、装置和计算机可读存储介质 |
CN116505522B (zh) * | 2023-06-28 | 2023-09-26 | 国网浙江省电力有限公司宁波供电公司 | 一种电力系统运行仿真模拟方法、仿真平台及设备 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102184297A (zh) * | 2011-05-12 | 2011-09-14 | 天津大学 | 适于微网暂态并行仿真的电气与控制系统解耦预测方法 |
CN103793562A (zh) * | 2014-01-05 | 2014-05-14 | 天津大学 | 基于fpga的有源配电网暂态实时仿真系统设计方法 |
CN103942372A (zh) * | 2014-04-04 | 2014-07-23 | 天津大学 | 基于fpga的有源配电网暂态实时仿真多速率接口方法 |
CN104361159A (zh) * | 2014-10-31 | 2015-02-18 | 国家电网公司 | 一种大规模电力系统暂态稳定时间空间并行仿真方法 |
CN104376158A (zh) * | 2014-11-05 | 2015-02-25 | 天津大学 | 一种用于矩阵指数的暂态仿真多时间尺度输出方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH06261453A (ja) * | 1993-03-08 | 1994-09-16 | Hitachi Ltd | 電力系統シミュレーション装置 |
CN101719182B (zh) * | 2009-12-11 | 2012-01-18 | 中国电力科学研究院 | 一种交直流电力系统分割并行电磁暂态数字仿真方法 |
US8397054B2 (en) * | 2009-12-23 | 2013-03-12 | L-3 Communications Integrated Systems L.P. | Multi-phased computational reconfiguration |
CN102592003B (zh) * | 2011-12-15 | 2013-11-13 | 天津大学 | 一种用于机电与电磁暂态混合仿真的数据交换方法 |
CN104573164B (zh) * | 2014-07-07 | 2017-10-31 | 贵州电力试验研究院 | 小步长系统等值历史电流源信息更新的硬件实现方法 |
CN104218600B (zh) * | 2014-09-22 | 2016-04-20 | 国家电网公司 | 一种多时间尺度动态仿真中电池储能系统模型的构建方法 |
-
2015
- 2015-05-22 CN CN201510264005.XA patent/CN106294897B/zh active Active
-
2016
- 2016-07-12 WO PCT/CN2016/089823 patent/WO2016188503A2/zh active Application Filing
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102184297A (zh) * | 2011-05-12 | 2011-09-14 | 天津大学 | 适于微网暂态并行仿真的电气与控制系统解耦预测方法 |
CN103793562A (zh) * | 2014-01-05 | 2014-05-14 | 天津大学 | 基于fpga的有源配电网暂态实时仿真系统设计方法 |
CN103942372A (zh) * | 2014-04-04 | 2014-07-23 | 天津大学 | 基于fpga的有源配电网暂态实时仿真多速率接口方法 |
CN104361159A (zh) * | 2014-10-31 | 2015-02-18 | 国家电网公司 | 一种大规模电力系统暂态稳定时间空间并行仿真方法 |
CN104376158A (zh) * | 2014-11-05 | 2015-02-25 | 天津大学 | 一种用于矩阵指数的暂态仿真多时间尺度输出方法 |
Also Published As
Publication number | Publication date |
---|---|
WO2016188503A3 (zh) | 2016-12-29 |
WO2016188503A2 (zh) | 2016-12-01 |
CN106294897A (zh) | 2017-01-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106294897B (zh) | 一种适用于电磁暂态多时间尺度实时仿真接口的实现方法 | |
CN111709107B (zh) | 一种含大规模电力电子设备的电力网络仿真方法及系统 | |
CN110058936B (zh) | 用于确定专用处理资源的资源量的方法、设备和计算机程序产品 | |
WO2022063247A1 (zh) | 神经网络结构搜索方法及装置 | |
CN109066669B (zh) | 电力系统仿真方法、电力系统仿真装置及电子设备 | |
CN109558248B (zh) | 一种用于确定面向海洋模式计算的资源分配参数的方法及系统 | |
CN111431996B (zh) | 用于资源配置的方法、装置、设备和介质 | |
CN112818535A (zh) | 建立电热仿真模型及获得电热仿真值的方法和装置 | |
CN115203126A (zh) | 一种算子融合处理方法、装置、设备及存储介质 | |
CN107301266B (zh) | 一种磷酸铁锂电池loc估算方法和系统 | |
Zhao et al. | High efficient parallel numerical surface wave model based on an irregular quasi-rectangular domain decomposition scheme | |
CN112557905B (zh) | 电池组及其数据处理方法、计算机设备、介质和车辆 | |
CN117009283A (zh) | 一种多核多芯片数据处理方法、装置、芯片及存储介质 | |
CN116227599A (zh) | 一种推理模型的优化方法、装置、电子设备及存储介质 | |
CN111092755B (zh) | 一种基于资源占用的边缘服务迁移仿真方法 | |
CN114020454A (zh) | 一种内存管理方法、装置、设备及介质 | |
CN115293083B (zh) | 集成电路时序预测方法、装置、电子设备及存储介质 | |
CN109193660A (zh) | 一种多端直流配电网实时仿真装置及其仿真方法 | |
Wei et al. | Segmented transmission delay based decoupling for parallel simulation of a distribution network | |
CN117709129A (zh) | 多尺度仿真方法、装置、设备及存储介质 | |
CN116204500A (zh) | 一种电网设备数据的处理方法、装置、设备及介质 | |
KR20230134307A (ko) | 심층 신경망 기반 분할 학습을 위한 모델 재분할 방법 및 장치 | |
CN116069506A (zh) | 容器资源配置方法、装置和电子设备 | |
Orlandi | Announcing the best transactions paper for the year 2015:" field-to-wire coupling model for the common mode in random bundles of twisted-wire pairs" | |
CN116862064A (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 |