CN104834791B - 一种光伏发电***的并行实时仿真方法 - Google Patents

一种光伏发电***的并行实时仿真方法 Download PDF

Info

Publication number
CN104834791B
CN104834791B CN201510268404.3A CN201510268404A CN104834791B CN 104834791 B CN104834791 B CN 104834791B CN 201510268404 A CN201510268404 A CN 201510268404A CN 104834791 B CN104834791 B CN 104834791B
Authority
CN
China
Prior art keywords
simulation
emulation
generating system
network
length
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
Application number
CN201510268404.3A
Other languages
English (en)
Other versions
CN104834791A (zh
Inventor
汪新星
刘旭
刘浩
杜航
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai Yuan Kuan Energy Science Co Ltd
Original Assignee
Shanghai Yuan Kuan Energy Science Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shanghai Yuan Kuan Energy Science Co Ltd filed Critical Shanghai Yuan Kuan Energy Science Co Ltd
Priority to CN201510268404.3A priority Critical patent/CN104834791B/zh
Publication of CN104834791A publication Critical patent/CN104834791A/zh
Application granted granted Critical
Publication of CN104834791B publication Critical patent/CN104834791B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Inverter Devices (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明涉及一种光伏发电***的并行实时仿真方法，光伏发电***的电力电子部分一般都通过含有电感和电容的滤波器并网，利用这个拓扑特点、以及电感电流和电容电压在几十微秒的时间内不会突变的原理，以滤波器电感和滤波器电容为界将光伏发电***分割，设计了受控电压源和受控电流源的接口方式，分别建立电力电子仿真部分和交流电网仿真部分，将电力电子仿真部分在FPGA上进行几微秒步长的实时仿真，将交流电网仿真部分放在CPU上进行几十微秒步长的实时仿真。与现有技术相比，本发明分别建立电力电子仿真部分和交流电网仿真部分，充分利用了***的拓扑结构和元件的电气特性实现***分割，可以对两个仿真部分实现不同步长的联合仿真。

Description

一种光伏发电***的并行实时仿真方法
技术领域
本发明涉及电力***仿真技术领域，尤其是涉及一种光伏发电***的并行实时仿真方法。
背景技术
出于人们对环境保护与可持续发展的关注，太阳能作为一种清洁、取之不尽用之不竭的能源受到了越来越多的重视。光伏发电是太阳能发电的一种主要形式，光伏发电***一般都包含了光伏电池，电力电子变换电路(DC-DC变换电路，DC-AC逆变电路)，LC滤波器，交流电网等部分。光伏发电***的工作依赖于控制器对DC-DC变换电路、DC-AC逆变电路等电力电子变换电路的精确控制来实现最大功率追踪、并网发电等正常运行功能；同时随着越来越多的光伏接入电网，光伏发电***也需要能够确保电网的安全稳定运行，因此光伏发电***的控制希望往往也被要求实现低电压穿越、孤岛检测等应对电网故障的功能。
为了确保光伏控制器的功能，无论在研发和还是在量产的阶段都需要对控制器做许多的测试。传统的测试方式是将控制器直接通过实物***上进行测试，但即使是通过一些小功率的微缩实物***来测试，也存在实验成本高、较危险(尤其是电网故障实验)、难于实现测试自动化等挑战。实时仿真器是在一个实时的硬件平台上利用数学模型来模拟实际***行为的装置，通过实时仿真器来测试控制器，可以对控制设备进行非常接近真实情况的测试与验证。这个方式具有安全，容易实现测试自动化等优点。
含光伏的电力***的实时仿真测试的挑战在于电力电子电路由于还有快速开关的器件，一般需要几微秒级的仿真步长才能准确仿真。目前工业界基本是利用FPGA(FieldProgrammable Gate Array，即可编程门阵列)的硬件并行性来实现这样小步长的实时仿真。但是FPGA上编程较复杂，乘法器等资源也有限，如将其于交流电网的实时仿真难度较大。交流电网一般只需几十微秒大步长即可准确仿真，且基于CPU的交流电网大步长实时仿真技术较为成熟。理想的话，最好是可以实时电力电子部分和交流电网部分的分割，然后分别利用FPGA与CPU并行仿真。
用于分割电气***的经典理论有戴维南定理和诺顿定理，这些理论的应用有限制，它们一般只适用于线性电路。但是光伏发电***一般不能满足这个限制，比如电力电子部分的电力电子开关，交流电网部分的断路器等都是非线性元件。
发明内容
本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种光伏发电***的并行实时仿真方法，利用逆变器的LC滤波器分割光伏发电***，将电力电子部分在FPGA上进行几微秒步长的仿真，将交流电网部分在CPU上进行几十微秒步长的仿真，从而实现联合仿真。
本发明的目的可以通过以下技术方案来实现：
一种光伏发电***的并行实时仿真方法，包括步骤：
1)根据光伏发电***的架构，分别建立电力电子仿真部分和交流电网仿真部分；
2)在电力电子仿真部分进行几微秒步长的小步长仿真，在交流电网仿真部分进行几十微秒步长的大步长仿真，以进行两个仿真部分的联合仿真。
所述步骤1)具体包括步骤：
101)利用光伏发电***中滤波器电感和滤波器电容将光伏发电***分割为电力电子子***与交流电网子***；
102)将交流电网子***等效为受控电压源，并与电力电子子***连接构成电力电子仿真部分；
103)将电力电子子***等效为受控电流源，并与交流电网子***连接构成交流电网仿真部分。
所述步骤101)具体为：以光伏发电***中滤波器电感和滤波器电容为界对光伏发电***进行分割，将滤波器电感一侧作为电力电子子***，将滤波器电容一侧作为交流电网子***。
所述电力电子仿真部分建立在FPGA上，所述交流电网仿真部分建立在CPU上。
所述光伏发电***包括光伏电池，该光伏电池依次通过电力电子电路、滤波器电感和滤波器电容与交流电网连接，所述光伏电池在电力电子仿真部分中被等效为受控电流源，且光伏电池的模型运行于CPU上。
所述电力电子仿真部分和交流电网仿真部分采用同一时钟源。
所述时钟源为FPGA上的高速时钟。
所述步骤2)中，大步长仿真的步长为小步长仿真步长的整数倍。
所述大步长仿真的每个仿真周期的起始时刻，两个仿真部分交换电容电压、光伏电流、以及滑动平均后的电感电流和光伏端电压。
所述滑动平均的窗口长度为大步长仿真的步长。
与现有技术相比，本发明具有以下优点：
1)根据光伏发电***的架构，分别建立电力电子仿真部分和交流电网仿真部分，充分利用了***的拓扑结构和元件的电气特性实现***分割，可以对两个仿真部分实现不同步长的联合仿真。
2)交流电网子***在电力电子子***中等效为受控电压源，利用了接口处有电容，且电容电压在几十微秒这样一个短的时间窗内不会突变的特性；电力电子子***在交流电网子***中用受控电流源代替，利用了接口处有电感，且电感电流在几十微秒这样一个短的时间窗内不会突变的特性，分割过程基于并网逆变器后的滤波器实现，尤其是通过为电力电子和交流电网选择各自合适的实时计算平台与仿真步长，降低光伏发电***实时仿真实现难度。
3)电力电子仿真部分建立在FPGA上，交流电网仿真部分建立在CPU上，可以实现对于电力电子仿真部分的小步长仿真，而对交流电网仿真部分实现大步长仿真。
4)电力电子仿真部分和交流电网仿真部分采用同一时钟源，确保了两个仿真部分的时间上的同步。
5)时钟源为FPGA上的高速时钟则可以实现更小的时间分割，提高仿真效果。
6)大步长仿真的步长为小步长仿真步长的整数倍，便于两个并行仿真部分交换数据。
附图说明
图1为本发明提出的光伏发电***并行实时仿真框架的示意图；
图2为本发明中电力电子仿真部分与交流电网仿真部分的仿真循环由同一时钟定时确保同步的示意图；
图3为本发明实施实例的双极光伏发电***的拓扑图；
图4为本发明实施实例的电力电子仿真部分具体建模框图；
图5为本发明实施实例的交流电网仿真部分具体建模框图；
图6为本发明实施实例的仿真结果波形图。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。本实施例以本发明技术方案为前提进行实施，给出了详细的实施方式和具体的操作过程，但本发明的保护范围不限于下述的实施例。
一种光伏发电***的并行实时仿真方法，其仿真框架如图1所示，包括步骤：
1)根据光伏发电***的架构，分别建立电力电子仿真部分和交流电网仿真部分，具体包括步骤：
101)利用光伏发电***中滤波器电感和滤波器电容将光伏发电***分割为电力电子子***与交流电网子***，具体为：以光伏发电***中滤波器电感和滤波器电容为界对光伏发电***进行分割，将滤波器电感一侧作为电力电子子***，将滤波器电容一侧作为交流电网子***；
102)将交流电网子***等效为受控电压源，并与电力电子子***连接构成电力电子仿真部分，利用了接口处有电容，且电容电压在几十微秒这样一个短的时间窗内不会突变的特性；
103)将电力电子子***等效为受控电流源，并与交流电网子***连接构成交流电网仿真部分，利用了接口处有电感，且电感电流在几十微秒这样一个短的时间窗内不会突变的特性。
经过步骤102)和步骤103)的等效过程，受控电流源的电流数值受滤波器的电感电流控制，受控电压源的电压数值受滤波器的电容电压控制。这样一个电压源和电流源的接口将光伏发电***的电力电子部分和交流电网部分的电气联系转为通过值传递互相影响，完成仿真模型的建模，从而实现***的拆分和并行仿真，特别指出，电力电子仿真部分和交流电网仿真部分实质上为电力电子仿真模型和交流电网仿真模型，它们之间通过交换电容电压，以及滑动平均后的电感电流等信息实现并行仿真，本实施例中，电力电子仿真模型和电力电子仿真部分，以及交流电网仿真模型和交流电网仿真部分都应该被理解为同一含义。
2)并行仿真：在电力电子仿真模型进行几微秒步长的小步长仿真，在交流电网仿真模型进行几十微秒步长的大步长仿真，以进行两个仿真模型的联合仿真，其中电力电子仿真模型建立在FPGA上，交流电网仿真模型建立在CPU上。
光伏发电***包括光伏电池，该光伏电池依次通过电力电子电路、滤波器电感和滤波器电容与交流电网连接，且在电力电子仿真模型中被等效为受控电流源，光伏电池本身的模型(光伏端电压和光伏输出电流的关系)实时运行在CPU部分，光伏电池和电力电子电路共同构成电力电子部分。
如图2所示，为了确保两个仿真模型的时间上的同步，FPGA仿真循环和CPU实时仿真循环由同一个时钟源来定时，即FPGA上的高速时钟。同时为了便于两个并行仿真平台交换数据，大步长的CPU实时仿真的仿真步长需要是小步长的FPGA仿真步长的整数倍。在每个CPU大步长仿真周期开始时，两个仿真循环交换电容电压、光伏电流，以及滑动平均后的电感电流和光伏端电压等信息，具体的，FPGA上滤波电感电流、光伏端电压经过滑动平均后传递给CPU的仿真循环，其滑动平均的窗长为CPU的仿真步长，而电容电压和光伏电流则由交流电网仿真模型传递给电力电子仿真模型。
下面具体以一个典型的双极光伏发电***为具体的实施实例，其拓扑具体如图3所示。该双极光伏发电***包含依次连接的光伏电池、平波电容Cpv、Boost电路、直流电容Cdc、DC-AC逆变电路、逆变器的RLC滤波器、变压器、负荷、输电线路和无穷大电源。
为了实现对此双极光伏发电***的实时仿真，利用逆变器的的并网滤波器的电感和电容将***分割为电力电子子***与交流电网子***，这两个子***采用受控电压源与受控电流源的方式来接口。即交流电网子***等效为受控电压源，并与电力电子子***连接构成电力电子仿真模型，电力电子子***等效为受控电流源，并与交流电网子***连接构成交流电网仿真模型，其中电力电子仿真模型如图4所示，交流电网仿真模型在CPU上的具体建模框图如图5所示。
将***分割后，电力电子仿真模型在FPGA上用几微秒的小步长实时仿真。为了在FPGA上实现的便利性，在本实例中Boost电路和DC-AC逆变电路采用受控源的方式等效建模。定义电力电子电路中受控开关的开关函数为S,当开关接收到PWM脉冲为高时，S为1，PWM脉冲为低时，S为0。
Boost电路中的受控电压源和受控电流源与其开关函数的关系为：
Vboost＝(1-Sboost)VCdc
Iboost＝(1-Sboost)ILpv
DC-AC全桥逆变电路等效模型中的受控电压源和受控电流源与其开关函数的关系为：
Vinv_a＝VCdcS1
Vinv_b＝VCdcS3
Vinv_c＝VCdcS5
Iinv_dc＝-(ILaS1+ILbS3+ILcS5)
其中：VCdc为Cdc的电压，ILpv为Lpv的电流，ILa,ILb,ILc分别为滤波器电感电流，Sboost、S1、S3、S5为图3中对应开关的开关函数。
在本实例中的光伏电池是等效为受控电流源，光伏电池的具体模型实时运行在CPU上，光伏电池的具体的数学模型如下，设在参考条件下，Isc为短路电流，Voc为开路电压，Im、Vm为最大功率点电流和电压，则当光伏阵列电压为V，其对应点电流为I：
C2＝(Vm/Voc-1)/ln(1-Im/Isc)
如上是光伏电池在参考的辐照强度和温度，即Gref、Tref时的数学模型，一般分别为1000W/m2和25℃，当太阳辐射变化和温度影响时，其数学模型修正如下：
ΔI＝α·G/Gref·ΔT+(G/Gref-1)·Isc
ΔV＝-β·ΔT-Rs·ΔI
ΔT＝Tc-Tref
其中：G为光伏电池当前太阳辐射强度(W/m2)，Tc为当前温度(℃)；α、β分别为参考太阳辐射强度下的电流温度系数(A/℃)和电压系数(V/℃)，对于单晶硅及多晶硅太阳电池其通常取得值为α＝0.0012Isc(A/℃)和β＝0.005Voc(V/℃)；Rs为光伏电池内阻。
取电力电子仿真模型的仿真步长为2.5微秒，交流电网仿真模型的仿真步长为50微秒。在1.5秒的时刻，太阳辐射强度从800W/m2阶跃到1200W/m2，图6中为并行实时仿真时的逆变器电感A相电流的波形。

Claims (4)

1.一种光伏发电***的并行实时仿真方法，其特征在于，包括步骤：
1)根据光伏发电***的架构，分别建立电力电子仿真部分和交流电网仿真部分，
所述光伏发电***包括光伏电池，该光伏电池依次通过电力电子电路、滤波器电感和滤波器电容与交流电网连接，
2)在电力电子仿真部分进行几微秒步长的小步长仿真，在交流电网仿真部分进行几十微秒步长的大步长仿真，以进行两个仿真部分的联合仿真；
所述步骤1)具体包括步骤：
101)利用光伏发电***中滤波器电感和滤波器电容将光伏发电***分割为电力电子子***与交流电网子***，
102)将交流电网子***等效为受控电压源，并与电力电子子***连接构成电力电子仿真部分，
所述光伏电池在电力电子仿真部分中被等效为受控电流源，且光伏电池的模型运行于CPU上；
103)将电力电子子***等效为受控电流源，并与交流电网子***连接构成交流电网仿真部分；
所述步骤101)具体为：以光伏发电***中滤波器电感和滤波器电容为界对光伏发电***进行分割，将滤波器电感一侧作为电力电子子***，将滤波器电容一侧作为交流电网子***；
所述电力电子仿真部分和交流电网仿真部分采用同一时钟源；
所述时钟源为FPGA上的高速时钟；
所述步骤2)中，大步长仿真的步长为小步长仿真步长的整数倍。
2.根据权利要求1所述的一种光伏发电***的并行实时仿真方法，其特征在于，所述电力电子仿真部分建立在FPGA上，所述交流电网仿真部分建立在CPU上。
3.根据权利要求1所述的一种光伏发电***的并行实时仿真方法，其特征在于，所述大步长仿真的每个仿真周期的起始时刻，两个仿真部分交换电容电压、光伏电流、以及滑动平均后的电感电流和光伏端电压。
4.根据权利要求3所述的一种光伏发电***的并行实时仿真方法，其特征在于，所述滑动平均的窗口长度为大步长仿真的步长。
CN201510268404.3A 2015-05-22 2015-05-22 一种光伏发电***的并行实时仿真方法 Active CN104834791B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510268404.3A CN104834791B (zh) 2015-05-22 2015-05-22 一种光伏发电***的并行实时仿真方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510268404.3A CN104834791B (zh) 2015-05-22 2015-05-22 一种光伏发电***的并行实时仿真方法

Publications (2)

Publication Number Publication Date
CN104834791A CN104834791A (zh) 2015-08-12
CN104834791B true CN104834791B (zh) 2018-09-28

Family

ID=53812675

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510268404.3A Active CN104834791B (zh) 2015-05-22 2015-05-22 一种光伏发电***的并行实时仿真方法

Country Status (1)

Country Link
CN (1) CN104834791B (zh)

Families Citing this family (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105608245A (zh) * 2015-10-29 2016-05-25 贵州黔燎原电气科技有限公司 大规模配电网多cpu独立解算的建模方法
CN106849152A (zh) * 2015-12-03 2017-06-13 甘肃省电力公司风电技术中心 一种光伏电站的综合仿真平台建立方法
CN107656143A (zh) * 2016-12-14 2018-02-02 中国电力科学研究院有限公司 新能源并网变流器的正负序阻抗测量模块、***及方法
CN106951646A (zh) * 2017-03-23 2017-07-14 南方电网科学研究院有限责任公司 一种电力***实时仿真中的多速率接口方法和装置
CN108063442A (zh) * 2017-12-30 2018-05-22 贵州大学 一种电力***交流电网实时仿真装置及其仿真方法
CN108416132B (zh) * 2018-02-28 2021-11-12 东南大学 一种分布式光伏集群的自动变步长仿真加速方法
CN112989545B (zh) * 2019-12-13 2023-02-03 中车唐山机车车辆有限公司 一种列车供电***的仿真建模方法及仿真建模装置
CN111596567B (zh) * 2020-04-27 2023-04-21 南方电网科学研究院有限责任公司 一种交直流电力***电磁暂态仿真装置
CN111858026B (zh) * 2020-06-10 2021-08-31 中国人民解放军海军航空大学航空作战勤务学院 高效的大规模多粒度仿真模型并行调度方法
CN112131812B (zh) * 2020-08-12 2023-10-27 北京华电天仁电力控制技术有限公司 一种多步长并行电力实时仿真***和方法
CN112928754B (zh) * 2021-02-22 2022-05-17 清华大学 基于状态离散事件驱动的光伏微网***仿真方法
CN115598994A (zh) * 2022-09-07 2023-01-13 贵州电网有限责任公司（Cn） 一种永磁直驱式风力发电***的多速率实时仿真方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102521473A (zh) * 2012-01-06 2012-06-27 江苏省电力试验研究院有限公司 基于DIgSILENT的并网型光伏仿真发电***
CN103942372A (zh) * 2014-04-04 2014-07-23 天津大学 基于fpga的有源配电网暂态实时仿真多速率接口方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3984117B2 (ja) * 2002-07-15 2007-10-03 シャープ株式会社 太陽光発電装置設計方法および太陽光発電装置設計支援装置

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102521473A (zh) * 2012-01-06 2012-06-27 江苏省电力试验研究院有限公司 基于DIgSILENT的并网型光伏仿真发电***
CN103942372A (zh) * 2014-04-04 2014-07-23 天津大学 基于fpga的有源配电网暂态实时仿真多速率接口方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
基于RTDS的光伏并网***实时仿真平台研究;苏丽萍等;《电力***保护与控制》;20120801;第110-115页 *
大功率光伏并网发电***RTDS的建模与仿真;黄鑫等;《电力***自动化》;20141125;第32-37页 *

Also Published As

Publication number Publication date
CN104834791A (zh) 2015-08-12

Similar Documents

Publication Publication Date Title
CN104834791B (zh) 一种光伏发电***的并行实时仿真方法
CN103997059B (zh) 三相单级光伏并网的电磁与机电暂态模型的建模方法
Can Model of a photovoltaic panel emulator in MATLAB-Simulink
CN105608244B (zh) 基于pscad/emtdc的分布式光伏发电***
CN104865847A (zh) 功率在环型数字与物理混合实时仿真的试验***及其试验方法
CN103020385A (zh) 基于RTDS电网500kV主网建模仿真***
CN108336758A (zh) 一种基于纹波关联法的光伏组件mppt算法
CN103124080B (zh) 一种光伏发电***模型的建模方法
CN103106328A (zh) 基于整体矢量拟合法的频率相关网络等值的生成方法
CN105978373A (zh) 实现微电网稳定的三相逆变器逆推滑模控制方法及***
CN103605850A (zh) 一种带子模块闭锁功能的mmc等效建模方法
CN103474989A (zh) 一种基于灵敏度分析的网络重构方法
CN204332198U (zh) 一种光伏发电并网模拟装置
CN103809650A (zh) 一种光伏发电***的等效建模方法
CN102521473B (zh) 基于DIgSILENT的并网型光伏仿真发电***
Saxena et al. Performance Analysis of P&O and Incremental Conductance MPPT Algorithms Under Rapidly Changing Weather Conditions.
CN108205595A (zh) 一种适用于谐波分析的pscad双馈风电机组仿真模型***
CN103279626B (zh) 一种仿真模型建立方法及***
Shen et al. Parameter identification of photovoltaic discrete-time equivalent model using the bat algorithm
CN204856067U (zh) 功率在环型数字与物理混合实时仿真的试验***
CN106253341B (zh) 并网逆变器损耗在线建模仿真方法
CN107563062A (zh) 一种基于动态相量法的交‑直‑交变频器建模方法
Qingsheng et al. Multi-rate Real-time Simulation Method for Photovoltaic Power Generation System Based on FPGA
Iero et al. SPICE modelling and experiments on a complete photovoltaic system including cells, storage elements, inverter and load
CN106849152A (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