基于FPGA的控制系统实时仿真框架设计方法与流程
本发明属于电力系统中电气设备测试技术领域，特别是涉及一种基于fpga的控制系统实时仿真框架设计方法。
背景技术：
电力系统实时仿真是与现实时间完全同步的电磁暂态仿真，其主要功能在于通过硬件在环仿真(hardware-in-the-loop-simulation)测试出厂的保护和控制设备。实时仿真器可以与待测设备相连，并模拟较为真实的暂态条件，避免待测设备对真实系统产生影响。从功率等级看，硬件在环仿真可分为信号型硬件在环仿真以及功率型硬件在环仿真，在信号型硬件在环仿真中，待测设备通常为二次控制器、保护设备、智能终端等，与实时仿真器之间进行低功率的信号交互。相对来说，功率型硬件在环仿真中，待测设备一般为实际电气设备，比如发动机、电力电子设备等，需要与实时仿真器之间实现较大的功率流动，因此在实时仿真器与待测设备之间往往需要功率放大器进行接口。
在电力系统实时仿真中，一般将电力系统仿真对象划分为电气系统和控制系统两部分进行求解。在电气系统中，其元件特性是以元件的伏安关系描述的，如线路、变压器、开关元件、电机、阻抗等；在控制系统中，其元件特性是以输入输出关系描述的，如传递函数、延迟环节、限幅环节等，需要基于控制系统进行建模。在传统的电力系统电磁暂态仿真中，控制系统主要用于对发电机的励磁控制、原动部分进行建模，在含电力电子设备的电力系统暂态仿真中，控制系统不仅可以对由各种线性环节、动态环节以及非线性环节构成的二次系统进行建模，在现代电力系统中，控制系统还承担针对高压直流输电、柔性交流输电技术等大功率电力电子设备的二次系统进行仿真的任务，同时，随着可再生能源发电、分布式电源等清洁能源在配电网中广泛接入，具有复杂非线性特性的分布式电源同样可以在控制系统中进行建模和仿真，因此控制系统在电力系统实时仿真中具有十分重要的作用。
然而二次系统以及各类分布式电源本身就具有较强的非线性特征，另外，随着现代电力系统中分布式电源、高压直流输电技术、柔性交流输电技术的大量应用，控制系统的仿真规模在电力系统实时仿真中所占的比重越来越高，对于控制系统的计算往往成为整个实时仿真最大的计算负担，此时就亟需使用新型底层硬件提升控制系统实时仿真的计算能力。
传统电力系统实时仿真主要基于rtds、hypersim、emegasim等商业实时仿真平台，然而现代电力系统具有了新的特点：1)现代电力系统中引入了大量高频电力电子设备及其二次系统，这使得实时仿真器中控制系统部分的计算资源消耗越来越大，对实时仿真器的计算性能提出了较高的要求，以rtds等实时仿真平台研究相关问题时，如不借助系统等值等方法，研究的花费和代价过大，也不利于问题的扩展。另一方面，系统等值会使得电网的运行特性发生变化，在涉及电力系统特性研究时，会因电网的等值化简影响到仿真的准确性；2)高压直流、柔性交流输电技术应用大量电力电子装置，要准确反映系统的暂态过程，需要采用更小的仿真步长，这更进一步加剧了计算资源的矛盾，并对仿真的实时性提出了更高的要求。因此，电力系统实时仿真面临计算能力与精度和经济性的双重挑战。
为解决现代电力系统实时仿真中控制系统的计算需求，需要先进的底层计算硬件以及快速的仿真算法。在底层计算硬件层面，大多实时仿真器采用基于risc(精简指令级计算机)的工作站、多dsp、多cpu计算机、pc机群、多核技术等，通过并行技术达到实时计算能力。这类方法涉及的大部分数据处理工作仍是串行进行。另外，为了实现大规模系统仿真，需要设置大量数据处理单元，处理单元之间的数据通讯又会成为计算速度的主要瓶颈。相比之下，基于现场可编程逻辑门阵列(field-programmablegatearray，fpga)的全硬件计算为实时仿真提供了一种新思路。fpga本质上具有完全可配置的固有硬件并行结构，其逻辑资源可配置为很多并行处理单元并实现多层级高度并行计算；同时，fpga芯片上具有大量嵌入式块ram，可配置为大量分布式rom或ram，其数据和地址宽度、端口数量皆可配置，而传统实时仿真器中内存和总线大多是共享的，且端口有限，因此限制了数据的传输效率；fpga允许使用流水线技术，加强了数据处理效率，并且，fpga还拥有大量传输速度极快的内部连线，不会引入过大的通讯延迟；最后，fpga具有可重构特性，可以针对仿真算例量身定做地调整实时仿真器，以达到最快的计算速度。但目前尚缺少采用fpga对电力系统仿真中的控制系统部分进行实时仿真的方法。
技术实现要素：
为了解决上述问题，本发明的目的在于提供一种基于fpga的控制系统实时仿真框架设计方法。
为了达到上述目的，本发明提供的基于fpga的控制系统实时仿真框架设计方法包括按顺序进行的下列步骤：
第一步：在离线环境下，采用基于fpga的控制系统基本元件对拟研究的电力系统的二次系统、分布式电源分别进行建模；
第二步：在离线环境下，利用图的连通性以及深度优先搜索算法，对上述分布式电源、二次系统模型进行拓扑识别，划分为nc个控制子系统；
第三步：在离线环境下，确定上述nc个控制子系统需要由电气系统读入的接口变量及个数ni，以及需要向电气系统读出的接口变量及个数no；
第四步：在离线环境下，分别计算nc个控制子系统进行一个时步解算所需的时钟周期数，根据fpga的驱动时钟频率f以及各个控制子系统中的最长时钟周期数nc，计算控制系统所对应的每一时步仿真所需的计算时间tc，其中tc＝nc/f；
第五步：在离线环境下，设定控制系统实时仿真步长δt；
第六步：在离线环境下，为nc个控制子系统分配独立的fpga硬件资源，同时根据上述选定的控制系统实时仿真步长δt，计算nc个控制子系统中各元件的具体仿真参数，并根据nc个控制子系统中各元件的连接关系以及基于fpga的控制系统元件输入输出模型，自动生成基于元件连接顺序进行求解的nc个控制子系统求解模块；
第七步：在离线环境下，根据由电气系统读入的接口变量及个数ni，以及需要向电气系统读出的接口变量及个数no，自动生成控制系统读入模块和控制系统输出模块；
第八步：在在线环境下，设置仿真时刻t＝0；
第九步：开始下一时步的仿真，令t＝t+δt；
第十步：从控制系统读入模块的存储器rami中以串行形式读出ni个接口数据，分配给nc个控制子系统求解模块；同时控制系统每个时钟周期均检测电气系统是否将no个接口数据传递过来，将传递过来的接口数据存入存储器rami中；
第十一步：控制系统进行一个时步计算，对nc个控制子系统进行仿真计算，求解出需要向电气系统读出的no个接口变量；其中nc个控制子系统在基于fpga的仿真框架中是通过分配独立的计算资源且并行求解的，并且控制系统求解时间由nc个控制子系统中计算耗时最长者决定；
第十二步：将nc个控制子系统求解出的no个接口变量按顺序组成串行数据流，并以串行的形式存入控制系统输出模块的存储器ramo中，由存储器ramo统一读出并存入电气系统的电力电子开关模块和电源模块的对应存储器中；
第十三步：判断仿真时间是否达到仿真终了时刻，如达到仿真终了时刻，则仿真结束；否则返回第九步。
在第十步中，所述的“从控制系统读入模块的存储器rami中以串行形式读出ni个接口数据，分配给nc个控制子系统求解模块”和“每个时钟周期均检测电气系统是否将接口数据传递过来，将传递过来的接口数据存入存储器rami中”这两个过程是完全独立的，并且能够并行进行。
本发明提供的基于fpga的控制系统实时仿真框架设计方法充分利用了fpga高度并行的硬件架构，具有较好的可行性与适用性，可有效提高控制系统的仿真效率，实现了控制系统的高速计算，为实现含高压直流、柔性交流输电、分布式电源的现代电力系统实时仿真提供了一种很好的解决思路，缓解了实时仿真中控制系统部分的计算压力。
附图说明
图1是本发明提供的基于fpga的控制系统实时仿真框架设计方法中控制系统求解框架示意图。
图2是本发明提供的基于fpga的控制系统实时仿真框架设计方法流程图。
图3是实施例中选取的作为待研究电力系统的典型光蓄混合发电系统示意图。
图4是实施例中控制子系统耗时示意图。
图5是实施例中逆变器输出的a相电流图。
图6是实施例中逆变器输出功率图。
图7是实施例中蓄电池输出功率图。
图8是光伏阵列输出功率图。
图9是蓄电池输出电压图。
图10是直流母线电压图
具体实施方式
现以图3所示的典型光蓄发电系统为例对本发明提供的基于fpga的控制系统实时仿真框架设计方法进行详细说明。
在该系统中，蓄电池组通过dc/dc换流器与光伏阵列并于直流母线a。其中，光伏阵列采用mppt控制；蓄电池组放电和充电时分别采用boost升压电路和buck降压电路模式，用于维持直流母线a的电压恒定；逆变器采用pq控制，维持整个光蓄发电系统的输出有功和无功功率恒定。直流母线a的电压控制在750v，无功参考值qref设为0var，保证单位功率因数运行，温度设置为298.15k。设系统最初的光照强度为1000w/m2，逆变器有功功率指令为10kw，系统达到稳态后，光照强度由1000w/m2下降为800w/m2，随后1s后有功功率指令降低为4kw。
如图1、图2所示，本实施例提供的基于fpga的控制系统实时仿真框架设计方法包括按顺序进行的下列步骤：
第一步：在离线环境下，采用基于fpga的控制系统基本元件对拟研究的电力系统的二次系统、分布式电源分别进行建模；在本实施例中，主要采用基于fpga的控制系统基本元件对光伏阵列、蓄电池组及其二次系统分别进行建模；
第二步：在离线环境下，利用图的连通性以及深度优先搜索算法，对上述分布式电源、二次系统模型进行拓扑识别，划分为nc个控制子系统；本实施例中将光蓄发电系统分为5个控制子系统，即nc＝5，分别为光伏阵列、蓄电池组、boost电路的mppt控制系统、boost/buck电路的控制系统以及逆变器控制系统，控制系统的解算时间由耗时最长的逆变器控制系统决定，如图4所示；
第三步：在离线环境下，确定上述nc个控制子系统需要由电气系统读入的接口变量及个数ni，以及需要向电气系统读出的接口变量及个数no；在本实施例中，由电气系统读入的接口变量的个数ni＝11，向电气系统读出的接口变量的个数no＝13；
第四步：在离线环境下，分别计算nc个控制子系统进行一个时步解算所需的时钟周期数，根据fpga的驱动时钟频率f以及各个控制子系统中的最长时钟周期数nc，计算控制系统所对应的每一时步仿真所需的计算时间tc，其中tc＝nc/f；在本实施例中，逆变器控制系统的解算时间耗时最长，为222个时钟周期，tc＝1.644μs；
第五步：在离线环境下，设定控制系统实时仿真步长δt，因为控制系统的解算首先要保证实时性，即tc≤δt，因此应在保证仿真精度的前提下选择实时仿真步长；在本实施例中，δt取为1.696μs；
第六步：在离线环境下，为nc个控制子系统分配独立的fpga硬件资源，同时根据上述选定的控制系统实时仿真步长δt，计算nc个控制子系统中各元件的具体仿真参数，并根据nc个控制子系统中各元件的连接关系以及基于fpga的控制系统元件输入输出模型，自动生成基于元件连接顺序进行求解的nc个控制子系统求解模块；
第七步：在离线环境下，根据由电气系统读入的接口变量及个数ni，以及需要向电气系统读出的接口变量及个数no，自动生成控制系统读入模块和控制系统输出模块；
第八步：在在线环境下，设置仿真时刻t＝0；
第九步：开始下一时步的仿真，令t＝t+δt；
第十步：从控制系统读入模块的存储器rami中以串行形式读出ni个接口数据，分配给nc个控制子系统求解模块；同时控制系统每个时钟周期均检测电气系统是否将no个接口数据传递过来，将传递过来的接口数据存入存储器rami中；
第十一步：控制系统进行一个时步计算，对nc个控制子系统进行仿真计算，求解出需要向电气系统读出的no个接口变量。其中nc个控制子系统在基于fpga的仿真框架中是通过分配独立的计算资源且并行求解的，并且控制系统求解时间由nc个控制子系统中计算耗时最长者决定；
第十二步：将nc个控制子系统求解出的no个接口变量按顺序组成串行数据流，并以串行的形式存入控制系统输出模块的存储器ramo中，由存储器ramo统一读出并存入电气系统的电力电子开关模块和电源模块的对应存储器中；
第十三步：判断仿真时间是否达到仿真终了时刻，如达到仿真终了时刻，则仿真结束；否则返回第九步。
本实施例的执行环境为altera公司的
在仿真精度方面，附图5～10比较了采用本发明方法与商业仿真软件pscad/emtdc的仿真结果。从图中可以看出，pscad/emtdc仿真结果与本发明方法的仿真结果在稳态与暂态过程中都能够完全吻合，二者的动态响应特性保持了高度一致，体现出了良好的仿真精度，充分验证了本发明方法的可行性。
以上算例测试结果证明，本发明提供的基于fpga的控制系统实时仿真框架设计方法充分利用了fpga高度并行的硬件架构，具有较好的可行性与适用性，可有效提高控制系统的仿真效率，实现了控制系统的高速计算，为实现含高压直流、柔性交流输电、分布式电源的现代电力系统实时仿真提供了一种很好的解决思路，缓解了实时仿真中控制系统部分的计算压力。
- 还没有人留言评论。精彩留言会获得点赞！