• 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 2019-05-13
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-12
  • 韩国再提4国联合申办世界杯 中国网友无视:我们自己来 2019-05-11
  • 中国人为什么一定要买房? 2019-05-11
  • 十九大精神进校园:风正扬帆当有为 勇做时代弄潮儿 2019-05-10
  • 粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 2019-05-09
  • 太原设禁鸣路段 设备在测试中 2019-05-09
  • 拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 2019-05-08
  • “港独”没出路!“梁天琦们”该醒醒了 2019-05-07
  • 陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 2019-05-06
  • 人民日报客户端辟谣:“合成军装照”产品请放心使用 2019-05-05
  • 【十九大·理论新视野】为什么要“建设现代化经济体系”?   2019-05-04
  • 聚焦2017年乌鲁木齐市老城区改造提升工程 2019-05-04
  • 【专家谈】上合组织——构建区域命运共同体的有力实践者 2019-05-03
  • 【华商侃车NO.192】 亲!楼市火爆,别忘了买车位啊! 2019-05-03
    • / 21
    • 下载费用:30 金币  

    重庆时时彩一码必中口诀: 基于FPGA的有源配电网暂态实时仿真系统设计方法.pdf

    摘要
    申请专利号:

    重庆时时彩单双窍门 www.4mum.com.cn CN201410028769.4

    申请日:

    2014.01.05

    公开号:

    CN103793562A

    公开日:

    2014.05.14

    当前法律状态:

    授权

    有效性:

    有权

    法律详情: 授权|||著录事项变更IPC(主分类):G06F 17/50变更事项:申请人变更前:天津大学变更后:广东电网公司变更事项:地址变更前:300072 天津市南开区卫津路92号变更后:510620 广东省广州市越秀区东风东路757号变更事项:申请人变更前:南方电网科学研究院有限责任公司 广东电网公司变更后:天津大学 南方电网科学研究院有限责任公司|||实质审查的生效IPC(主分类):G06F 17/50申请日:20140105|||公开
    IPC分类号: G06F17/50 主分类号: G06F17/50
    申请人: 天津大学; 南方电网科学研究院有限责任公司; 广东电网公司
    发明人: 王成山; 丁承第; 李鹏; 于浩; 董旭柱; 于力; 黄红远; 黄小耘
    地址: 300072 天津市南开区卫津路92号
    优先权:
    专利代理机构: 天津市北洋有限责任专利代理事务所 12201 代理人: 杜文茹
    PDF完整版下载: PDF下载
    法律状态
    申请(专利)号:

    CN201410028769.4

    授权公告号:

    |||||||||

    法律状态公告日:

    2016.09.14|||2016.02.17|||2014.06.11|||2014.05.14

    法律状态类型:

    授权|||著录事项变更|||实质审查的生效|||公开

    摘要

    一种基于FPGA的有源配电网暂态实时仿真系统设计方法,分为离线仿真环境和在线仿真环境,离线仿真环境负责计算每时步计算的总时钟周期数ntotal以及仿真用时ttotal,根据实时仿真用时ttotal设置仿真步长Δt,并将根据仿真步长Δt计算出的相关参数以及读取到的基本参数信息上传至基于FPGA的在线仿真环境;在线仿真环境完成仿真的实时计算,通过有限状态机控制仿真状态,每一时步包括每个步长中各类元件的历史项电流源列向量计算、总历史项电流源列向量形成、线性方程组求解及更新步骤,各类元件的历史项电流源列向量计算以及更新步骤是完全独立的,可并发处理。本发明保证了整个系统暂态仿真过程中的实时性,具有较好的可行性与适用性。

    权利要求书

    权利要求书
    1.  一种基于FPGA的有源配电网暂态实时仿真系统设计方法,其特征在于,有源配电网暂态实时仿真系统分为离线仿真环境和在线仿真环境,其中离线仿真环境负责实时仿真用时ttotal的计算,根据实时仿真用时ttotal设置仿真步长Δt,并将根据仿真步长Δt计算出的相关参数以及读取到的基本参数信息上传至基于FPGA的在线仿真环境;在线仿真环境完成仿真的实时计算,包括每个步长中各类元件的历史项电流源列向量计算、总历史项电流源列向量形成、线性方程组求解以及更新步骤,具体包括如下步骤:
    第一步:在离线环境下,采用基本元件对有源配电网进行建模,读取基本无源元件、线路元件、电源元件、断路器元件、电力电子开关元件的基本参数信息,统计得到所述各类元件的数量;
    第二步:在离线环境下,分别计算基本无源元件、线路元件、电源元件、断路器元件、电力电子开关元件的历史项电流源列向量求解的时钟周期数:nh,RLC,nh,LINE,nh,SOURCE,nh,BREAKER,nh,PE,以及更新步骤的时钟周期数:nu,RLC,nu,LINE,nu,SOURCE,nu,BREAKER,nu,PE,其中,nh表示历史项电流源列向量求解的时钟周期数,取大于等于0的整数,nu表示更新步骤的时钟周期数,取大于等于0的整数,RLC表示基本无源元件,LINE表示线路元件,SOURCE表示电源元件,BREAKER表示断路器元件,PE表示电力电子开关元件;
    第三步:在离线环境下,计算第一步所述的各类元件历史项电流源列向量形成总历史项电流源列向量所需的时钟周期数nhist,计算采用并行的矩阵向量乘法实现线性方程组求解的时钟周期数nmatrix;
    第四步:在离线环境下,计算每一步仿真所需的计算时钟周期总数为ntotal=max(nh,RLC,nh,Line,nh,SOURCE,nh,BREAKER,nh,PE)+nhist+nmatrix+max(nu,RLC,nu,Line,nu,SOURCE,nu,BREAKER,nu,PE)+nother,其中第一步所述的各类元件的历史项电流源列向量求解以及更新步骤的计算并发进行,max函数表示取最大值,nother表示一些零散操作的总时钟周期数;
    第五步:在离线环境下,根据FPGA的驱动时钟频率f以及时钟周期总数ntotal,得到第四步中每一步仿真所需的实际时间ttotal,ttotal=ntotal/f,设定仿真计算步长Δt,Δt需满足ttotal<Δt以保证仿真实时性;
    第六步:在离线环境下,根据仿真步长Δt计算第一步所述的各类元件的等效电导,形成节点电导矩阵,计算所述的各类元件模型中历史项电流源以及更新运算所需的计算参数,计算节点电导矩阵的逆矩阵;
    第七步:将第六步已得到的等效电导、各类元件模型中历史项电流源以及更新运算所需的计算参数、节点电导矩阵的逆矩阵、第一步所述的基本参数信息、第五步所述的仿真步长Δt上传至基于FPGA的在线仿真环境;
    第八步:在在线环境下,设置仿真时刻t=0,全局控制??橹杏邢拮刺姆抡孀刺障凶刺?IDLE),启动仿真;
    第九步:在在线环境下,仿真状态进入状态一(STEP_I),仿真计时器开始计时,计算第一步所述的各类元件的历史项电流源,生成历史项电流源列向量,其中所述的各类元件的 计算是并行的,各类元件的计算任务完成后将各自的结束信号置高电平;
    第十步:在在线环境下,对第九步所述的各类元件的结束信号进行逻辑与操作,当该信号为高电平,即所有元件完成历史项电流源计算后,将所有元件生成的历史项电流源列向量组合,得到总的历史项电流源列向量Ih,并进行存储,生成状态一结束信号(end_STEP_I);
    第十一步:在在线环境下,仿真状态进入状态二(STEP_II),根据断路器开关状态读取相应的节点电导矩阵的逆矩阵G-1,采用并行的矩阵向量乘法由逆矩阵G-1及总的历史项电流源列向量Ih计算节点电压列向量v,生成状态二结束信号(end_STEP_II);
    第十二步:在在线环境下,仿真状态进入状态三(STEP_III),对所述的各类元件分别进行更新运算,计算每个元件的端电压以及支路电流,并进行存储,其中所述的各类元件的更新计算是并行的;
    第十三步:在在线环境下,将用户指定的仿真结果传回离线环境,以便于用户查看;
    第十四步:在在线环境下,判断仿真计时器是否计时至Δt时,如满足条件,则生成步长结束信号(end_Δt),否则仿真器等待直至计时至Δt;
    第十五步:判断仿真时间是否达到仿真终了时刻,如达到仿真终了时刻,则仿真结束,否则返回第九步。

    2.  根据权利要求1所述的基于FPGA的有源配电网暂态实时仿真系统设计方法,其特征在于,第一步所述的基本元件包括:基本无源元件、线路元件、电源元件、断路器元件、电力电子开关元件。

    3.  根据权利要求1所述的基于FPGA的有源配电网暂态实时仿真系统设计方法,其特征在于,第一步所述的基本参数信息包括:基本无源元件的电阻、电感和电容参数,线路元件的电阻和电感矩阵,电源元件的幅值、频率和初始相位,断路器元件的开路电阻、闭合电阻、断路器动作时间,电力电子开关元件表示开路的电阻、电容以及表示闭合的电感,各类元件的节点编号,节点电导矩阵维数。

    4.  如权利要求1所述的一种基于FPGA的有源配电网暂态实时仿真系统设计方法,其特征在于,实时仿真中的仿真状态一(STEP_I)、仿真状态二(STEP_II)和仿真状态三(STEP_III)的控制调度,是由有限状态机实现的。

    说明书

    说明书基于FPGA的有源配电网暂态实时仿真系统设计方法
    技术领域
    本发明涉及一种电力系统仿真设计方法。特别是涉及一种基于FPGA的有源配电网暂态实时仿真系统设计方法。
    背景技术
    有源配电网暂态实时仿真是一种与现实时钟同步的数字仿真方式,通过将实时仿真器与实际物理设备相连可开展各种控制与?;ぷ爸玫目⒂氩馐怨ぷ?,并可模拟光照、风速变化、电压跌落、短路故障、甩负荷等多种运行场景下的有源配电网复杂暂态过程,可有效降低研发及试验成本,避免待测设备对实际系统的影响。有源配电网暂态实时仿真主要用于分布式发电(储能)单元控制器设计与试验、有源配电网能量管理系统测试、有源配电网?;ぷ爸玫魇缘?。
    在有源配电网中,为了精确刻画包括电力电子开关动作等在内的一系列快动态过程,暂态仿真的步长变得越来越小。从仿真角度看,电力电子设备的存在会造成计算矩阵时变、步长间开关动作、数值震荡等一系列问题,对于这些问题的精确求解需要耗费较长的计算时间与较多的计算资源。更重要的是,这些额外的计算时间的消耗往往是不可估计的。由于分布式电源种类多样,其控制系统的数学模型十分复杂,且逻辑判断多,具有强非线性特征,从而导致控制系统求解规模庞大?;诔9鍯PU处理器或DSP等串行硬件的传统实时仿真器,在较小的仿真步长下很难实现详细建模的暂态实时仿真。相比之下,基于现场可编程逻辑门阵列(field-programmable gate array,FPGA)的全硬件计算为实时仿真提供了一种新思路。FPGA本质上具有完全可配置的固有硬件并行结构,其逻辑资源可配置为很多并行处理单元并实现多层级高度并行计算;同时,FPGA芯片上具有大量嵌入式块RAM,可配置为大量分布式ROM或RAM,其数据和地址宽度、端口数量皆可配置,而传统实时仿真器中内存和总线大多是共享的,且端口有限,限制了数据的传输效率;FPGA允许使用流水线技术,加强了数据处理效率,并且,FPGA还拥有大量传输速度极快的内部连线,不会引入过大的通讯延迟。
    因此基于FPGA的有源配电网暂态实时仿真系统允许更快的计算速度和更小的计算步长,可以为有源配电网暂态实时仿真中对于步长有要求严格的部分提供速度和精度支持,体积小巧,在开发周期与成本的经济性上也更具优势;而且考虑到传统实时仿真器的模型与算法等内容基本上不对用户开放,进一步研究开发准确、高效的仿真算法不大可能。相比之下,基于FPGA的有源配电网暂态实时仿真器具有可开发性、可拓展性,能够为有源配电网控制、?;げ呗缘难芯?,新设备调试等提供测试平台。
    暂态实时仿真问题的基本求解方法可以分为节点分析法(nodal analysis)以及状态变量分析法(state space analysis)两类。相对于状态变量分析,节点分析法在算法实现难度、仿真计算速度等方面具有较大优势,因此在EMTP、PSCAD/EMTDC等暂态离线仿真工具以及 RTDS、HYPERSIM等暂态实时仿真工具中,都以节点分析法作为基本框架。
    暂态仿真节点分析法包含2个基本步骤:
    1)采用某种数值积分方法(如梯形积分法)对系统中动态元件的特性方程进行差分化,得到等效计算电导与历史项电流源并联形式的诺顿等效电路。以附图1所示的电感支路为例,其基本伏安关系方程如式(1)所示,应用梯形积分法后可得到式(2)和(3)形式的差分方程。
    uk(t)-um(t)=Ldikm(t)dt---(1)]]>
    ikm(t)=Δt2L[uk(t)-um(t)]+Ih[t-Δt]---(2)]]>
    Ih(t-Δt)=ikm(t-Δt)+Δt2L[uk(t-Δt)-um(t-Δt)]]]>
    差分方程式(2)可看作是对一个值为GeqΔt/(2L)的电导与历史项电流源并联的诺顿等效电路形式。
    2)联立整个电气系统的差分方程,可形成如式(4)中所示节点电导矩阵G,求解该方程可得到系统中各节点电压的瞬时值,进而得到支路电压和支路电流。该求解过程的不断推进即可完成整个系统的暂态过程求解。
    Gv=Ih   (4)
    式(4)所示的节点电导方程为线性方程组,可采用各种成熟的数值软件包进行求解,而对于电气系统中的各种非线性元件,如非线性阻抗、电机等可采用分段线性化、伪非线性、预测校正法、补偿法等方法进行局部处理,而系统整体上仍是对式(4)形式的线性方程组的求解。在基于FPGA的有源配电网暂态实时仿真器中,需要充分发挥FPGA硬件并行优势以及流水线技术的系统求解框架,才能满足基于小步长的有源配电网暂态实时仿真的需求。为此,本发明提出了一种基于FPGA的有源配电网暂态实时仿真系统设计方法。该设计方法提出的系统框架能够充分发挥FPGA硬件并行性以及流水线架构的技术优势,保证了整个系统暂态仿真过程中的实时性,并为采用较小步长精确求解有源配电网的详细动态过程提供了可靠的基础。
    发明内容
    本发明所要解决的技术问题是,提供一种基于FPGA的有源配电网暂态实时仿真系统设计方法。其能够充分发挥FPGA硬件并行性以及流水线架构技术优势,保证了整个系统暂态仿真过程中的实时性,并为采用较小步长精确求解有源配电网的详细动态过程提供了可靠基础。
    本发明所采用的技术方案是:一种基于FPGA的有源配电网暂态实时仿真系统设计方法,有源配电网暂态实时仿真系统分为离线仿真环境和在线仿真环境,其中离线仿真环境负责实时仿真用时ttotal的计算,根据实时仿真用时ttotal设置仿真步长Δt,并将根据仿真步长Δt计算出的相关参数以及读取到的基本参数信息上传至基于FPGA的在线仿真环境;在线仿真环境完成仿真的实时计算,包括每个步长中各类元件的历史项电流源列向量计算、总历史项电流 源列向量形成、线性方程组求解以及更新步骤,具体包括如下步骤:
    第一步:在离线环境下,采用基本元件对有源配电网进行建模,读取基本无源元件、线路元件、电源元件、断路器元件、电力电子开关元件的基本参数信息,统计得到所述各类元件的数量;
    第二步:在离线环境下,分别计算基本无源元件、线路元件、电源元件、断路器元件、电力电子开关元件的历史项电流源列向量求解的时钟周期数:nh,RLC,nh,LINE,nh,SOURCE,nh,BREAKER,nh,PE,以及更新步骤的时钟周期数:nu,RLC,nu,LINE,nu,SOURCE,nu,BREAKER,nu,PE,其中,nh表示历史项电流源列向量求解的时钟周期数,取大于等于0的整数,nu表示更新步骤的时钟周期数,取大于等于0的整数,RLC表示基本无源元件,LINE表示线路元件,SOURCE表示电源元件,BREAKER表示断路器元件,PE表示电力电子开关元件;
    第三步:在离线环境下,计算第一步所述的各类元件历史项电流源列向量形成总历史项电流源列向量所需的时钟周期数nhist,计算采用并行的矩阵向量乘法实现线性方程组求解的时钟周期数nmatrix;
    第四步:在离线环境下,计算每一步仿真所需的计算时钟周期总数为ntotal=max(nh,RLC,nh,Line,nh,SOURCE,nh,BREAKER,nh,PE)+nhist+nmatrix+max(nu,RLC,nu,Line,nu,SOURCE,nu,BREAKER,nu,PE)+nother,其中第一步所述的各类元件的历史项电流源列向量求解以及更新步骤的计算并发进行,max函数表示取最大值,nother表示一些零散操作的总时钟周期数;
    第五步:在离线环境下,根据FPGA的驱动时钟频率f以及时钟周期总数ntotal,得到第四步中每一步仿真所需的实际时间ttotal,ttotal=ntotal/f,设定仿真计算步长Δt,Δt需满足ttotal<Δt以保证仿真实时性;
    第六步:在离线环境下,根据仿真步长Δt计算第一步所述的各类元件的等效电导,形成节点电导矩阵,计算所述的各类元件模型中历史项电流源以及更新运算所需的计算参数,计算节点电导矩阵的逆矩阵;
    第七步:将第六步已得到的等效电导、各类元件模型中历史项电流源以及更新运算所需的计算参数、节点电导矩阵的逆矩阵、第一步所述的基本参数信息、第五步所述的仿真步长Δt上传至基于FPGA的在线仿真环境;
    第八步:在在线环境下,设置仿真时刻t=0,全局控制??橹杏邢拮刺姆抡孀刺障凶刺?,启动仿真;
    第九步:在在线环境下,仿真状态进入状态一,仿真计时器开始计时,计算第一步所述的各类元件的历史项电流源,生成历史项电流源列向量,其中所述的各类元件的计算是并行的,各类元件的计算任务完成后将各自的结束信号置高电平;
    第十步:在在线环境下,对第九步所述的各类元件的结束信号进行逻辑与操作,当该信号为高电平,即所有元件完成历史项电流源计算后,将所有元件生成的历史项电流源列向量组合,得到总的历史项电流源列向量Ih,并进行存储,生成状态一结束信号;
    第十一步:在在线环境下,仿真状态进入状态二,根据断路器开关状态读取相应的节点电导矩阵的逆矩阵G-1,采用并行的矩阵向量乘法由逆矩阵G-1及总的历史项电流源列向量Ih计算节点电压列向量v,生成状态二结束信号;
    第十二步:在在线环境下,仿真状态进入状态三,对所述的各类元件分别进行更新运算,计算每个元件的端电压以及支路电流,并进行存储,其中所述的各类元件的更新计算是并行的;
    第十三步:在在线环境下,将用户指定的仿真结果传回离线环境,以便于用户查看;
    第十四步:在在线环境下,判断仿真计时器是否计时至Δt时,如满足条件,则生成步长结束信号,否则仿真器等待直至计时至Δt;
    第十五步:判断仿真时间是否达到仿真终了时刻,如达到仿真终了时刻,则仿真结束,否则返回第九步。
    第一步所述的基本元件包括:基本无源元件、线路元件、电源元件、断路器元件、电力电子开关元件。
    第一步所述的基本参数信息包括:基本无源元件的电阻、电感和电容参数,线路元件的电阻和电感矩阵,电源元件的幅值、频率和初始相位,断路器元件的开路电阻、闭合电阻、断路器动作时间,电力电子开关元件表示开路的电阻、电容以及表示闭合的电感,各类元件的节点编号,节点电导矩阵维数。
    实时仿真中的仿真状态一、仿真状态二和仿真状态三的控制调度,是由有限状态机实现的。
    本发明的基于FPGA的有源配电网暂态实时仿真系统设计方法,能够充分发挥FPGA硬件并行性以及流水线架构的技术优势,保证了整个系统暂态仿真过程中的实时性,并为采用较小步长精确求解有源配电网的详细动态过程提供了可靠的基础。本发明的仿真结果与商业仿真软件的仿真结果在稳态与暂态过程中都能够完全吻合,二者的动态响应特性保持了高度一致,体现出了良好的仿真精度,充分验证了本发明提出的适于FPGA实现的有源配电网暂态实时仿真系统框架设计方法的可行性。本发明具有较好的可行性与适用性,为实现含分布式电源、储能设备的有源配电网暂态实时仿真中提供了一种很好的解决思路。
    附图说明
    附图1是电感支路示意图;
    附图2是基于FPGA的有源配电网暂态实时仿真系统求解框架;
    附图3是基于有限状态机的全局控制??槭疽馔?;
    附图4是基于FPGA的基本无源元件历史量求解框图;
    附图5是基于FPGA的线路元件历史量求解框图;
    附图6是基于FPGA的电源类元件历史量求解框图;
    附图7是基于FPGA的电力电子开关元件历史量求解框图;
    附图8是基于FPGA的线性方程组求解框图;
    附图9是本发明的基于FPGA的有源配电网暂态实时仿真系统设计方法的流程图;
    附图10是算例系统接线图;
    附图11是节点11流向节点12的三相电流;
    附图12是节点12的三相电压。
    具体实施方式
    下面结合实施例和附图对本发明的基于FPGA的有源配电网暂态实时仿真系统设计方法做出详细说明。
    本发明的基于FPGA的有源配电网暂态实时仿真系统设计方法,如附图2所示,有源配电网暂态实时仿真系统分为离线仿真环境和在线仿真环境,其中离线仿真环境负责实时仿真用时ttotal的计算,根据实时仿真用时ttotal设置仿真步长Δt,并将根据仿真步长Δt计算出的相关参数以及读取到的基本参数信息上传至基于FPGA的在线仿真环境;在线仿真环境完成仿真的实时计算,包括每个步长中各类元件的历史项电流源列向量计算、总历史项电流源列向量形成、线性方程组求解以及更新步骤,具体如图9所示,包括如下步骤:
    第一步:在离线环境下,采用基本元件对有源配电网进行建模,读取各类基本元件的基本参数信息,统计得到所述各类元件的数量,其中,所述的基本元件包括:基本无源元件、线路元件、电源元件、断路器元件、电力电子开关元件等;所述的基本参数信息包括:基本无源元件的电阻、电感和电容参数,线路元件的电阻和电感矩阵,电源元件的幅值、频率和初始相位,断路器元件的开路电阻、闭合电阻、断路器动作时间,电力电子开关元件表示开路的电阻、电容以及表示闭合的电感,各类元件的节点编号,节点电导矩阵维数。
    第二步:在离线环境下,分别计算基本无源元件、线路元件、电源元件、断路器元件、电力电子开关元件等各类元件的历史项电流源列向量求解的时钟周期数:nh,RLC,nh,LINE,nh,SOURCE,nh,BREAKER,nh,PE,以及更新步骤的时钟周期数:nu,RLC,nu,LINE,nu,SOURCE,nu,BREAKER,nu,PE,其中,nh表示历史项电流源列向量求解的时钟周期数,取大于等于0的整数,nu表示更新步骤的时钟周期数,取大于等于0的整数;
    有源配电网中包括电源、电机、变压器、线路与开关模型等多类基本元件??悸堑接性磁涞缤猩婕暗脑隙?,这里仅针对具有典型处理方式的基本无源元件、线路元件、电源类元件和电力电子开关元件进行说明,并给出离线环境下这几类元件的历史项电流源列向量求解和更新步骤的时钟周期总数的计算方法。
    (1)基本无源元件
    基本无源元件包括电阻R、电感L和电容C以及由它们组成的RL、RC复合元件等。对于此类元件,其历史量求解和更新步骤的一般形式如式(5)和式(6)所示。式中,根据R、L、C组合形式不同,A1、A2、Geq的具体表达式也有所不同,如表1所示。
    Ih(t-Δt)=A1v(t-Δt)+A2i(t-Δt)   (5)
    i(t)=Geqv(t)+Ih(t-Δt)   (6)
    表1基本无源元件计算参数表达式


    对于基本无源元件的仿真,可充分利用FPGA的高度并行性并实现深度流水线结构。以历史项求解为例,如附图4所示。根据历史项求解计算公式,A1、A2及上一时步的支路电压v(t-Δt)和支路电流i(t-Δt)可同时从4个RAM中读出,并通过2个浮点数乘法器同时进行乘法运算,最后通过浮点数加法器求出历史项Ih(t-Δt),存入RAMhist中,用于下一时步计算。在该??橹?,由RAM读取数据需要2个时钟周期,浮点数乘法和加法运算使用了集成编译环境自带的IP核,其计算时钟数nmultiply和nadd分别为5和7个时钟周期。因此,??槌跏佳邮蔽?4个时钟周期。在流水线架构下,第一个基本无源元件的计算结果将在14个时钟周期后得到,当NRLC个无源元件依次通过历史项求解??榻写硎?,其总时钟周期数nh,RLC=14+(NRLC-1)个时钟周期。同理,其更新步骤包含两个浮点数加法和一个浮点数乘法,其总时钟周期数nu,RLC=21+(NRLC-1)个时钟周期。
    (2)线路元件
    线路元件的历史项计算公式以及更新计算公示如式(7)和式(8)所示:
    Ih1(t-Δt)Ih2(t-Δt)Ih3(t-Δt)=A11A12A13A21A22A23A31A32A33v1(t-Δt)v2(t-Δt)v3(t-Δt)+B11B12B13B21B22B23B31B32B33i1(t-Δt)i2(t-Δt)i3(t-Δt)---(7)]]>
    i1(t)i2(t)i3(t)=A11A12A13A21A22A23A31A32A33v1(t)v2(t)v3(t)+Ih1(t-Δt)Ih2(t-Δt)Ih3(t-Δt)]]>
    以历史项求解为例,如附图5所示。[A]和[B]均为3x3矩阵,可分别使用RAMA1、RAMA2、RAMA3、RAMB1、RAMB2、RAMB3存储每一个线路元件[A]和[B]矩阵的一行。仿真开始时,[A]和[B]以及上一时步的支路电压向量和支路电流向量同时从RAM中读出,并依照附图5所示流程进行并行矩阵向量点乘运算以及并行向量加法运算。当每一个线路元件的历史项电流源Ih1(t-Δt)、Ih2(t-Δt)、Ih3(t-Δt)计算得到后,将其重新转化为浮点数,并通过延时、整合、提取操作后存入存储器RAMhist。附图5中的Delay_nclk??楸硎窘衝个时钟周期的延时,多为数据对齐而设置,Fix2Float和Float2Fix??橛糜诙ǖ闶透〉闶南嗷プ?。当NLINE个线路元件依次通过历史项求解??榻写硎?,求解过程的时钟周期数nh,LINE=22+(3*NLINE-1)。另外,其更新过程的时钟周期数nu,LINE=33+(3*NLINE-1)。
    (3)电源类元件
    电源??榭梢源戆ǖハ嗟缪乖?、单相电流源、三相电压源、三相电流源以及受控电压源和受控电流源等多种电源。以交流源为例,其输出量一般都如式(9)所示:
    f(t)=Asin(ωt+θ)   (9)
    其中,A表示幅值,ω表示角频率,θ表示初始相位。采用查表法实现这一函数的求解,其原理如附图6所示。
    首先使用小于仿真步长Δt的采样步长Ts提前计算出一个周期内所有的单位正弦值sin(ωt+θ),并存入RAMsin,同时每一个交流源的初始相位θ推算出对应于RAMsin的初始地址并依次存入RAMθ中,同时RAMθ中存储的每一个相位都随仿真推进而递增,其增量Inc为 仿真步长与采样步长的比值Δt/Ts。这样就无需通过复杂的算法计算非线性正弦函数,只需通过RAMθ的地址数据从RAMsin中读取单位正弦值并与幅值A进行浮点数乘法即可,每个交流源的幅值A同样提前存入存储器RAMA中。当NSOURCE个电源元件依次通过历史项求解??榻写硎?,求解过程的时钟周期数nh,SOURCE=9+(NSOURCE-1)。另外,其更新步骤用于求解的时钟周期数nu,SOURCE=7+(NSOURCE-1)。
    (4)电力电子开关元件
    本发明使用如下方法进行开关建模,具体方法为:开关闭合时使用小电感模拟,开关断开时使用小电容与电阻串联形式模拟,式(10)列写了使用梯形法差分得到的开关断开和闭合时的特性方程。
    i(t)=Δt2Lsu(t)+Δt2Lsu(t-Δt)+i(t-Δt)]]>
    i(t)1Δt2Cs+Rsu(t)-1Δt2Cs+Rsu(t-Δt)-Δt2Cs-RsΔt2Cs+Rsi(t-Δt)---(10)]]>
    式(10)中Ls表示开关闭合时的电感,Cs表示开关断开时的电容,Rs表示开关断开时的电阻,i(t)和i(t-Δt)表示本时步以及上一时步的开关电流,u(t)和u(t-Δt)表示本时步以及上一时步的开关两侧电压,Δt表示仿真步长。若式(10)中开关断开和闭合时的导纳满足式(11)则保证了开关状态切换时导纳矩阵不变,仅依靠历史量变化即可切换开关状态。
    Δt2Ls=1Δt2Cs+Rs---(11)]]>
    使用该方法时,不论开关如何变化,导纳矩阵始终保持不变,因此只需存储一个逆矩阵即可,极大地缓解了逆矩阵的存储压力。
    对于开关元件的仿真实现可充分利用FPGA的高度并行性,并实现深度流水线结构。以其历史量求解为例进行说明,其求解公式如式(12)所示,由式(10)中的两个表达式综合而成,具体FPGA实现框图如附图7所示。
    i(t)=Gu(t)+Ih(t-Δt)   (12)
    Ih(t-Δt)=A1u(t-Δt)+A2i(t-Δt)   
    式(12)中,A1和A2的取值由开关状态决定,RAMv、RAMi及RAMhist用于存储开关电压、开关电流以及历史量,RAMA1_open、RAMA1_closed、RAMA2_open和RAMA2_closed分别用于存储每一个开关断开或闭合时A1和A2的取值A1_open、A1_closed、A2_open和A2_closed,RAMstate用于存储每个IGBT的开关状态state。如附图8所示,每一时步计算开始时,A1_open、A1_closed、A2_open、A2_closed以及state同时从存储器中读出,其中A1_open、A1_closed、A2_open和A2_closed分别通过两个数据分配器,并由通过开关状态state判断该时步A1和A2的取值。上一时步的开关端电压v(t-Δt)和开关电流i(t-Δt)延迟一个时钟读出,与数据分配器的输出结果A1、A2保持数据对齐,并同时通过2个浮点数乘法器进行乘法运算,最后通过浮点数加法器求出历史量Ih(t-Δt),存入RAMhist中,用于下一时步计算。当NPE个电力电子开关元件依次通过历史项求解??榻?处理时,求解过程的时钟周期数nh,PE=15+(NPE-1)。另外,其更新过程的时钟周期数nu,PE=21+(NPE-1)。
    第三步:在离线环境下,计算第一步所述的各类元件历史项电流源列向量形成总历史项电流源列向量所需的时钟周期数nhist,计算采用并行的矩阵向量乘法实现线性方程组求解的时钟周期数nmatrix;
    附图2中包含总历史项电流源列向量生成??橐约盎诓⑿芯卣笙蛄砍朔ǖ南咝苑匠套榍蠼饽??,其所需的时钟周期数如下面所述。
    在总历史项电流源列向量生成??橹?,需要根据各类元件计算出的历史项电流源列向量,叠加求出总历史项电流源。在该过程中叠加运算是在定点数数制下完成的,需要考虑浮点数与定点数的相互转换所消耗的时钟周期nconvert,叠加运算的时钟周期数等于矩阵维数Nmatrix,其中浮点数与定点数的相互转换使用了集成编译环境自带的IP核,其计算时钟数nconvert为6,因此该步骤的时钟周期数nhist=2*nconvert+Nmatrix。
    线性方程组求解往往是暂态仿真中最为耗时的部分,为了达到较快的计算速度,保证仿真实时性,可以预存电导矩阵逆矩阵,其计算公式如式(10)所示。对于维数为Nmatrix的矩阵,使用Nmatrix个RAM进行存储,每个RAM存储矩阵的一行数据。此时,可将线性方程组求解转化为Nmatrix个可并行的向量点乘运算,如式(11)所示。该??樵诶奂哟砩嫌胂呗吩嗤?,这里不再赘述,其计算需要先后考虑浮点数乘法、浮点数与定点数的相互转换以及定点数累加,因此其时钟周期数nmatrix=nmultiply+2nconvert。??槭迪址绞饺绺酵?所示,通过该??榭杉扑愕玫矫扛鼋诘愕乃彩钡缪怪?。
    v=YIh,Y=G-1   (10)
    v1v2...vN=Y1,1Y1,2....Y1,NY2,1Y2,2....Y2,N................YN,1YN,2....YN,NIh,1Ih,2....Ih,N---(11)]]>
    考虑到断路器使用双电阻模型进行模拟,即断路器断开时,使用阻值较大的电阻表示,断路器闭合时,使用阻值很小的电阻表示,断路器状态变化造成了系统电导矩阵的变化。因此在实时仿真中,所有可能的矩阵结构都需要提前计算并存储至RAM中,需要存储的计算矩阵个数为2Nbreaker,其中Nbreaker为断路器个数,通过每个断路器的当前状态形成状态序列,从而读取相应的电导矩阵逆矩阵。
    第四步:在离线环境下,如附图2所示,第一步所述的各类元件的历史项电流源列向量求解以及更新步骤的计算并发进行,它们的求解时钟周期数由耗时最长者确定,因此计算每一步仿真所需的计算时钟周期总数为ntotal=max(nh,RLC,nh,Line,nh,SOURCE,nh,BREAKER,nh,PE,......)+nhist+nmatrix+max(nu,RLC,nu,Line,nu,SOURCE,nu,BREAKER,nu,PE,......)+nother,max函数表示取最大值,nother表示一些零散操作的总时钟周期数,如生成功能??榻崾藕诺牟僮鞯?;
    第五步:在离线环境下,根据FPGA的驱动时钟频率f以及时钟周期总数ntotal,得到第四步中每一步仿真所需的实际时间ttotal,ttotal=ntotal/f,设定仿真计算步长Δt,Δt需满足ttotal<Δt以保证仿真实时性;
    第六步:在离线环境下,根据仿真步长Δt计算第一步所述的各类元件的等效电导,形成节点电导矩阵,计算所述的各类元件模型中历史项电流源以及更新运算所需的计算参数,计算节点电导矩阵的逆矩阵;
    第七步:将第六步已得到的等效电导、各类元件模型中历史项电流源以及更新运算所需的计算参数、节点电导矩阵的逆矩阵、第一步所述的基本参数信息、第五步所述的仿真步长Δt等上传至基于FPGA的在线仿真环境;
    如上所述的第一步~第七步,是在离线环境下完成的,而在线环境中的仿真状态控制是通过有限状态机实现的,如附图3所示,设置了四种运行状态并为每种状态分配一个3位二进制数表示,分别为空闲状态(IDLE:000);仿真状态一(STEP_I:001);仿真状态二(STEP_II:010);仿真状态三(STEP_III:100)。
    该有限状态机的运行流程为:仿真器初始处于空闲状态IDLE,当仿真开始信号Start拉高后,系统进入仿真状态一STEP_I;当总历史项电流源列向量生成后,仿真状态一结束信号end_STEP_I置1,促使仿真进入仿真状态二STEP_II;当线性方程组求解完成后,系统拉高仿真状态二结束信号end_STEP_II,促使仿真进入仿真状态三STEP_III;当仿真计时器计时至Δt,拉高步长结束信号end_Δt,该时步计算结束,促使仿真重新回到仿真状态一STEP_I并开始下一时步的仿真。由此,实时仿真在仿真状态一STEP_I、状态二STEP_II、状态三STEP_III三个状态间循环推进,直至仿真结束。当复位信号Reset拉高时,系统可重新回到空闲状态。具体流程如下所述。
    第八步:在在线环境下,设置仿真时刻t=0,全局控制??橹杏邢拮刺姆抡孀刺障凶刺?IDLE),启动仿真;
    第九步:在在线环境下,仿真状态进入状态一(STEP_I),仿真计时器开始计时,计算第一步所述的各类元件的历史项电流源,生成历史项电流源列向量,其中所述的各类元件的计算是并行的,各类元件的计算任务完成后将各自的结束信号置高电平;
    第十步:在在线环境下,对第九步所述的各类元件的结束信号进行逻辑与操作,当该信号为高电平,即所有元件完成历史项电流源计算后,将所有元件生成的历史项电流源列向量组合,得到总的历史项电流源列向量Ih,并进行存储,生成状态一结束信号(end_STEP_I);
    第十一步:在在线环境下,仿真状态进入状态二(STEP_II),根据断路器开关状态读取相应的节点电导矩阵的逆矩阵G-1,采用并行的矩阵向量乘法由逆矩阵G-1及总的历史项电流源列向量Ih计算节点电压列向量v,生成状态二结束信号(end_STEP_II);
    第十二步:在在线环境下,仿真状态进入状态三(STEP_III),对所述的各类元件分别进行更新运算,计算每个元件的端电压以及支路电流,并进行存储,其中所述的各类元件的更新计算是并行的;
    第十三步:在在线环境下,将用户指定的仿真结果传回离线环境,以便于用户查看;
    第十四步:在在线环境下,判断仿真计时器是否计时至Δt时,如满足条件,则生成步长结束信号(end_Δt),否则仿真器等待直至计时至Δt;
    第十五步:判断仿真时间是否达到仿真终了时刻,如达到仿真终了时刻,则仿真结束,否则返回第九步。
    该设计方法提出的系统框架能够充分发挥FPGA硬件并行性以及流水线架构的技术优势,保证了整个系统暂态仿真过程中的实时性,并为采用较小步长精确求解有源配电网的详细动态过程提供了可靠的基础。
    本发明中所述的实时仿真中的仿真状态一(STEP_I)、仿真状态二(STEP_II)和仿真状态三(STEP_III)的控制调度,是由有限状态机实现的。
    本发明提出的基于FPGA的有源配电网暂态实时仿真系统设计方法,属于电力系统仿真领域,特别适用于含分布式电源、储能装置的有源配电网实时暂态仿真。本发明提出了基于FPGA的有源配电网暂态实时仿真系统RTDG(Real-Time Transient Simulator for Distributed Generation and Microgrid)。这里以欧盟低压有源配电网算例系统为例进行说明,如附图10所示,系统中含有多种线路与负荷类型。实施步骤详细说明如下:
    第一步:在离线环境下,采用基本元件对有源配电网进行建模,读取基本无源元件、线路元件、电源元件、断路器元件、电力电子开关元件等各类基本元件的基本参数信息,统计得到各类元件的数量,在该算例中,电源元件3个、基本无源元件42个、线路元件15个,断路器元件1个,电力电子开关元件0个;
    第二步:在离线环境下,计算基本无源元件、线路元件、电源元件、断路器元件、电力电子开关等元件的历史项电流源列向量求解的时钟周期总数(nh,RLC,nh,LINE,nh,SOURCE,nh,BREAKER,nh,PE,......)以及更新步骤的时钟周期总数(nu,RLC,nu,LINE,nu,SOURCE,nu,BREAKER,nu,PE,......)。在该算例中nh,RLC=55,nh,LINE=66,nh,SOURCE=11,nh,BREAKER=0,nh,PE=0,nu,RLC=62,nu,LINE=77,nu,SOURCE=13,nu,BREAKER=22,nu,PE=0;
    第三步:在离线环境下,计算各类元件历史项电流源列向量形成总历史项电流源列向量所需的处理时钟周期数nhist,计算采用并行矩阵向量乘法实现线性方程组求解的时钟周期数nmatrix;在该算例中Nmatrix=67,nhist=79,nmatrix=17;
    第四步:在离线环境下,计算每一步仿真所需的总计算时钟周期数为ntotal=max(nh,RLC,nh,Line,nh,SOURCE,nh,BREAKER,nh,PE,......)+nhist+nmatrix+max(nu,RLC,nu,Line,nu,SOURCE,nu,BREAKER,nu,PE,......)+nother,其中各类元件的历史项电流源列向量求解以及更新步骤的计算并发进行,max函数表示取最大值,nother=25,在该算例中,依照公式计算出的ntotal为264;
    第五步:在离线环境下,根据FPGA的驱动时钟频率f以及ntotal,得到每一步计算所需的实际时间ttotal,ttotal=ntotal/f,设定仿真计算步长Δt,Δt需满足ttotal<Δt以保证仿真实时性。在该算例中,时钟频率f为135MHz,因此ttotal为1.956μs,Δt设定为2μs;
    第六步:在离线环境下,根据仿真步长Δt计算第一步所述的各类元件的等效电导,形成节点电导矩阵,计算所述的各类元件模型中历史项电流源以及更新运算所需的计算参数,计算节点电导矩阵的逆矩阵;
    第七步:将第六步已得到的等效电导、各类元件模型中历史项电流源以及更新运算所需的计算参数、节点电导矩阵的逆矩阵、第一步所述的基本参数信息、第五步所述的仿真步长Δt上传至基于FPGA的在线仿真环境;
    第八步:在在线环境下,设置仿真时刻t=0,全局控制??橹杏邢拮刺姆抡孀刺障凶刺?IDLE),启动仿真;
    第九步:在在线环境下,仿真状态进入状态一(STEP_I),仿真计时器开始计时,计算第一步所述的各类元件的历史项电流源,生成历史项电流源列向量,其中所述的各类元件的计算是并行的,各类元件的计算任务完成后将各自的结束信号置高电平;
    第十步:在在线环境下,对第九步所述的各类元件的结束信号进行逻辑与操作,当该信号为高电平,即所有元件完成历史项电流源计算后,将所有元件生成的历史项电流源列向量组合,得到总的历史项电流源列向量Ih,并进行存储,生成状态一结束信号(end_STEP_I);
    第十一步:在在线环境下,仿真状态进入状态二(STEP_II),根据断路器开关状态读取相应的节点电导矩阵的逆矩阵G-1,采用并行的矩阵向量乘法由逆矩阵G-1及总的历史项电流源列向量Ih计算节点电压列向量v,生成状态二结束信号(end_STEP_II);
    第十二步:在在线环境下,仿真状态进入状态三(STEP_III),对所述的各类元件分别进行更新运算,计算每个元件的端电压以及支路电流,并进行存储,其中所述的各类元件的更新计算是并行的;
    第十三步:在在线环境下,将用户指定的仿真结果传回离线环境,以便于用户查看;
    第十四步:在在线环境下,判断仿真计时器是否计时至Δt时,如满足条件,则生成步长结束信号(end_Δt),否则仿真器等待直至计时至Δt;
    第十五步:判断仿真时间是否达到仿真终了时刻,如达到仿真终了时刻,则仿真结束,否则返回第九步。
    该算例的执行环境为Altera公司的IV GX FPGA530官方开发板??迮溆蠸tratix IV系列FPGA EP4SGX530KH40C2N,该芯片包含531200个逻辑单元,212480个自适应逻辑???,1280个M9K存储器,64个M144K存储器,1024个18x18专用乘法器,8个PLL以及744个I/O。除了EP4SGX530KH40C2N芯片,开发板还提供了多个频率的时钟电路,3个用户可配置按钮,大量外部存储器,PCI Express插槽,10/100/1000Ethernet接口等外围电路。
    在仿真速度方面,通过本发明提出的基于FPGA的有源配电网暂态实时仿真系统设计方法,能够充分发挥FPGA硬件并行性以及流水线架构,保证整个系统暂态仿真的实时性;附图11~12比较了采用本发明方法的基于FPGA的有源配电网暂态实时仿真系统与商业仿真软件Matlab/SimPowerSystems的仿真结果。从附图中可以看出,Matlab/SimPowerSystems的仿真结果与FPGA的仿真结果在稳态与暂态过程中都能够完全吻合,二者的动态响应特性保持了高度一致,体现出了良好的仿真精度,充分验证了本发明提出的适于FPGA实现的有源配电网暂态实时仿真系统框架设计方法的可行性。
    以上算例测试结果证明,本发明提出的一种基于FPGA的有源配电网暂态实时仿真系统设计方法具有较好的可行性与适用性,为实现含分布式电源、储能设备的有源配电网暂态实时仿真中提供了一种很好的解决思路。

    关 键 词:
    基于 FPGA 有源 配电网 实时 仿真 系统 设计 方法
      专利查询网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
    关于本文
    本文标题:基于FPGA的有源配电网暂态实时仿真系统设计方法.pdf
    链接地址://www.4mum.com.cn/p-6156888.html
    关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服客服 - 联系我们

    [email protected] 2017-2018 www.4mum.com.cn网站版权所有
    经营许可证编号:粤ICP备17046363号-1 
     


    收起
    展开
  • 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 2019-05-13
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-12
  • 韩国再提4国联合申办世界杯 中国网友无视:我们自己来 2019-05-11
  • 中国人为什么一定要买房? 2019-05-11
  • 十九大精神进校园:风正扬帆当有为 勇做时代弄潮儿 2019-05-10
  • 粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 2019-05-09
  • 太原设禁鸣路段 设备在测试中 2019-05-09
  • 拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 2019-05-08
  • “港独”没出路!“梁天琦们”该醒醒了 2019-05-07
  • 陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 2019-05-06
  • 人民日报客户端辟谣:“合成军装照”产品请放心使用 2019-05-05
  • 【十九大·理论新视野】为什么要“建设现代化经济体系”?   2019-05-04
  • 聚焦2017年乌鲁木齐市老城区改造提升工程 2019-05-04
  • 【专家谈】上合组织——构建区域命运共同体的有力实践者 2019-05-03
  • 【华商侃车NO.192】 亲!楼市火爆,别忘了买车位啊! 2019-05-03