CN115021602A - 基于虚拟同步机的三电平并网逆变器的双环模型预测方法 - Google Patents
基于虚拟同步机的三电平并网逆变器的双环模型预测方法 Download PDFInfo
- Publication number
- CN115021602A CN115021602A CN202210573684.9A CN202210573684A CN115021602A CN 115021602 A CN115021602 A CN 115021602A CN 202210573684 A CN202210573684 A CN 202210573684A CN 115021602 A CN115021602 A CN 115021602A
- Authority
- CN
- China
- Prior art keywords
- voltage
- value
- current
- output
- 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.)
- Withdrawn
Links
- 238000000034 method Methods 0.000 title claims abstract description 53
- 230000001360 synchronised effect Effects 0.000 title claims abstract description 45
- 239000013598 vector Substances 0.000 claims abstract description 73
- 238000005070 sampling Methods 0.000 claims abstract description 37
- 238000012216 screening Methods 0.000 claims abstract description 9
- 238000004364 calculation method Methods 0.000 claims description 13
- 239000003990 capacitor Substances 0.000 claims description 7
- 230000000295 complement effect Effects 0.000 claims description 6
- 239000000126 substance Substances 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 abstract description 2
- 230000033228 biological regulation Effects 0.000 description 8
- 238000010586 diagram Methods 0.000 description 8
- 230000005284 excitation Effects 0.000 description 4
- 238000011217 control strategy Methods 0.000 description 3
- 230000001276 controlling effect Effects 0.000 description 3
- 102100029469 WD repeat and HMG-box DNA-binding protein 1 Human genes 0.000 description 2
- 101710097421 WD repeat and HMG-box DNA-binding protein 1 Proteins 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000011084 recovery Methods 0.000 description 2
- 230000001105 regulatory effect Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 238000013016 damping Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000002401 inhibitory effect Effects 0.000 description 1
- 230000009191 jumping Effects 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 230000007935 neutral effect Effects 0.000 description 1
- 230000002441 reversible effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02M—APPARATUS FOR CONVERSION BETWEEN AC AND AC, BETWEEN AC AND DC, OR BETWEEN DC AND DC, AND FOR USE WITH MAINS OR SIMILAR POWER SUPPLY SYSTEMS; CONVERSION OF DC OR AC INPUT POWER INTO SURGE OUTPUT POWER; CONTROL OR REGULATION THEREOF
- H02M7/00—Conversion of ac power input into dc power output; Conversion of dc power input into ac power output
- H02M7/42—Conversion of dc power input into ac power output without possibility of reversal
- H02M7/44—Conversion of dc power input into ac power output without possibility of reversal by static converters
- H02M7/48—Conversion of dc power input into ac power output without possibility of reversal by static converters using discharge tubes with control electrode or semiconductor devices with control electrode
- H02M7/483—Converters with outputs that each can have more than two voltages levels
- H02M7/487—Neutral point clamped inverters
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/381—Dispersed generators
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02M—APPARATUS FOR CONVERSION BETWEEN AC AND AC, BETWEEN AC AND DC, OR BETWEEN DC AND DC, AND FOR USE WITH MAINS OR SIMILAR POWER SUPPLY SYSTEMS; CONVERSION OF DC OR AC INPUT POWER INTO SURGE OUTPUT POWER; CONTROL OR REGULATION THEREOF
- H02M1/00—Details of apparatus for conversion
- H02M1/12—Arrangements for reducing harmonics from ac input or output
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02M—APPARATUS FOR CONVERSION BETWEEN AC AND AC, BETWEEN AC AND DC, OR BETWEEN DC AND DC, AND FOR USE WITH MAINS OR SIMILAR POWER SUPPLY SYSTEMS; CONVERSION OF DC OR AC INPUT POWER INTO SURGE OUTPUT POWER; CONTROL OR REGULATION THEREOF
- H02M7/00—Conversion of ac power input into dc power output; Conversion of dc power input into ac power output
- H02M7/42—Conversion of dc power input into ac power output without possibility of reversal
- H02M7/44—Conversion of dc power input into ac power output without possibility of reversal by static converters
- H02M7/48—Conversion of dc power input into ac power output without possibility of reversal by static converters using discharge tubes with control electrode or semiconductor devices with control electrode
- H02M7/53—Conversion of dc power input into ac power output without possibility of reversal by static converters using discharge tubes with control electrode or semiconductor devices with control electrode using devices of a triode or transistor type requiring continuous application of a control signal
- H02M7/537—Conversion of dc power input into ac power output without possibility of reversal by static converters using discharge tubes with control electrode or semiconductor devices with control electrode using devices of a triode or transistor type requiring continuous application of a control signal using semiconductor devices only, e.g. single switched pulse inverters
- H02M7/5387—Conversion of dc power input into ac power output without possibility of reversal by static converters using discharge tubes with control electrode or semiconductor devices with control electrode using devices of a triode or transistor type requiring continuous application of a control signal using semiconductor devices only, e.g. single switched pulse inverters in a bridge configuration
- H02M7/53871—Conversion of dc power input into ac power output without possibility of reversal by static converters using discharge tubes with control electrode or semiconductor devices with control electrode using devices of a triode or transistor type requiring continuous application of a control signal using semiconductor devices only, e.g. single switched pulse inverters in a bridge configuration with automatic control of output voltage or current
- H02M7/53875—Conversion of dc power input into ac power output without possibility of reversal by static converters using discharge tubes with control electrode or semiconductor devices with control electrode using devices of a triode or transistor type requiring continuous application of a control signal using semiconductor devices only, e.g. single switched pulse inverters in a bridge configuration with automatic control of output voltage or current with analogue control of three-phase output
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2300/00—Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
- H02J2300/20—The dispersed energy generation being of renewable origin
- H02J2300/22—The renewable source being solar energy
- H02J2300/24—The renewable source being solar energy of photovoltaic origin
-
- 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
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/50—Photovoltaic [PV] energy
- Y02E10/56—Power conversion systems, e.g. maximum power point trackers
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Inverter Devices (AREA)
Abstract
本发明涉及电力电子变换控制技术领域，尤其涉及一种基于虚拟同步机的三电平并网逆变器的双环模型预测方法，通过VSG给定电流值及电压值计算出k+1时刻电流参考值；通过电压给定值筛选出三电平变换器输出的27个电压矢量中合适的电压向量；根据并网电流采样值和电压矢量计算k+1时刻电流预测值；根据k+1时刻电流预测值计算k时刻参考电压值，通过中点电压变化量、电压向量和参考电压值通过价值函数计算出差值最小的价值函数对应的序列，并输出开关状态。本发明输出稳定，产生的更低的总谐波失真。
Description
技术领域
本发明涉及电力电子变换控制技术领域，尤其涉及一种基于虚拟同步机的三电平并网逆变器的双环模型预测方法。
背景技术
虚拟同步机(VSG)的一般控制方法采用PI电流内环、PQ功率外环的控制方案或者简单地PR控制方法等，对于更优的控制方案没有更多的研究。近年来，由于数字信号处理器(Digital Signal Processors，DSP)和现场可编程门阵列(Field Programmable GateArray，FPGA)的快速发展，一批新型的控制算法不断涌现，其中模型预测控制算法(ModelPredictive Control，MPC)拥有实现方便，稳态性能好，动态响应快，无超调的优点受到青睐。MPC包括有限控制集模型预测控制(Finite Control Set Model Predictive Control，FCS-MPC)和连续控制集模型预测控制(Continuous Control Set Model PredictiveControl，CCS-MPC)，其中FCS-MPC因为其动态快、无需调制环节，一直是国内外的热点。通过MPC控制和VSG的结合，可以进一步提升VSG提供电网支撑的能力。
为了进一步优化光伏并网的输出容量以及电流的谐波，逆变器拓扑方面可以选择使用中点箝位型(Neutral Point Clamped,NPC)三电平逆变器它的优点是功率管电压应力为直流母线电压的一半；相对于其他三电平逆变器拓扑，其等效电容少；理论分析统一性强，可直接扩展到更多电平逆变器拓扑等。它主要缺点是箝位二极管是快恢复二极管且需要承受全功率电流和反向恢复电压；存在直流母线电容电压不平衡问题，需要进行中点电压平衡控制。因此VSG下的MPC控制中点电压平衡将是研究的重点之一。
发明内容
本发明的目的是为了提供一种响应较快且输出更稳定的虚拟同步机(VSG)的模型预测控制(MPC)的控制方法。
为解决以上技术问题，本发明的技术方案为：基于虚拟同步机的三电平并网逆变器的双环模型预测方法，包括：
步骤1：虚拟同步机根据下垂控制原理输出k时刻的电压给定值v(k)，并通过定子电压方程得出电流给定值i*(k)，同时采集得出并网电流采样值i(k)和并网电压采样值e(k)；
步骤2：通过k时刻的电流给定值i*(k)预测出k+1时刻电流参考值i*(k+1)；
步骤3：获取三电平变换器的输出电压并分为P组不同状态的电压矢量u(k)，每个电压矢量都有不同的开关序列；
步骤4：根据电压给定值v(k)和P个电压矢量u(k)作差值比较，筛选出P组电压矢量u(k)中离电压给定值v(k)最接近的Q组，并将这Q组作为新的逆变器输出电压向量ui(k)，作为双环模型预测的电压环，实现分步骤控制中的电压控制；
步骤5：根据k时刻的并网电流采样值i(k)、并网采样电压值e(k)和电压矢量u(k)，计算k+1时刻采样电流预测值i(k+1)，作为模型预测控制下的电流环，实现分步骤控制中的电流控制；步骤6：根据k+1时刻电流参考值i*(k+1)、电压矢量u(k)和并网电流采样值i(k)，计算k时刻参考电压值u*(k)；
步骤7：根据逆变器的中点电流计算k+1时刻的中点电压变化量ΔVNP(k+1)；
步骤8：根据k+1时刻中点电压变化量ΔVNP(k+1)、逆变器输出电压向量ui(k)及k时刻参考电压值u*(k)，通过价值函数gmin计算出参考电压值和电压向量差值最小时的价值函数，选取最小价值函数时的逆变器输出电压向量ui(k)，并输出ui(k)的开关序列；
步骤9：依据步骤8的开关序列给予逆变器选通脉冲。
本发明具有如下有益效果：
与现有技术相比，本发明的显著优点在于：(1)采用MPC控制代替一般的PI控制，改进虚拟同步机并提高VSG的输出质量；(2)在保持原有MPC控制方法的基础上，不增加控制方法复杂度，通过电压电流分步控制的方法达到电压电流双环控制的效果；(3)能够减少输出电压、电流谐波，提高输出波形质量。
附图说明
图1为本发明流程示意图；
图2为本发明实施例中NPC型三电平的拓扑结构图；
图3为本发明实施例中三电平空间矢量分布图；
图4为基于VSG的MPC控制框图；
图5为图4中VSG控制部分框图；
图6为图4中MPC部分框图；
图7为实施例1中采用电压筛选前后THD值对比图；其中，
图7(a)为采用一般的模型预测控制方法的THD值示意图；
图7(b)为采用改进的双环模型预测方法的THD值示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，下面结合附图和具体实施例对本发明作进一步详细说明。
请参考图1至图7，本发明采用改进虚拟同步机的三电平并网逆变器的双环模型预测方法，通过模型预测控制代替PI控制改进了对VSG的控制策略，并实现了对电压电流两个变量的控制与调节。
结合图4，实现本发明目的的技术解决方案为：一种改进虚拟同步机的三电平并网逆变器的双环模型预测方法，通过模型预测控制代替传统PI控制改进了VSG的控制策略，体现了模型预测控制在VSG上的优越性。双环模型预测策略的具体实现是电压电流的分步骤控制；首先检测多电平逆变器直流母线的中点与正极之间的电容瞬时电压、并网采样电压值与并网采样电流值，并确定该两个电容瞬时电压的偏差。之后传感器采样的有关数据经过VSG生成当前时刻的参考电压及电流；通过电压给定值筛选电压向量作为双环模型预测的电压环，基于电压筛选后的电流模型预测控制作为电流环。以上的电压电流环可以认为是分别通过电压控制步骤和电流控制步骤实现。
本发明同时考虑了电压和电流对MPC的影响，提出的FSF-MPC方法通过分步骤控制电压电流的方法，达到了模型预测双环控制的功能，通过使用优化的电压向量筛选方法，从而获得了与传统VSG的PI电压电流双环控制相比更高的响应能力。同时，与经典的MPC方法相比，使用所提出的基于虚拟同步电机电压电流双闭环MPC(FSF-MPC)方法分步控制电压电流两个变量，首先对电压进行筛选，再以筛选过的电压向量作为电流模型预测控制的参考电压矢量，从而在输出电压电流中产生的更低的总谐波失真(THD)。
本发明基于的虚拟同步机的三电平并网逆变器的双环模型预测方法，步骤包括：
步骤1：设定合适的有功功率和无功功率的值输入虚拟同步机(VirtualSynchronous Generator，VSG)，虚拟同步机根据下垂控制原理输出k时刻的电压给定值v(k)，并通过定子电压方程得出电流给定值i*(k)，同时采集得出并网电流采样值i(k)和并网电压采样值e(k)；
其中，步骤1中测出电流给定值i*(k)，具体步骤如下：
步骤101，通过电压型虚拟同步机的有功功率控制部分计算出输出电压给定值v(k)的相角θ，通过电压型虚拟同步机的无功功率控制部分计算出输出电压给定值v(k)的电压幅值E，由此可得出v(k)在αβ坐标系下的数学表达式；
步骤102，分析虚拟同步机的电机特性，将输出的电压给定值v(k)通过定子电压方程转换为线电流即电流给定值i*(k)。
步骤2：通过步骤1中k时刻的电流给定值i*(k)预测出k+1时刻电流参考值i*(k+1)；
其中，预测计算方法在于：假定电流给定值和电压给定值是恒定的正弦波，便可在k时刻的值预测出k+1时刻的值。
步骤3：获取中点钳位型三电平变换器的输出电压并分为27个不同状态的电压矢量u(k)，27组电压矢量u(k)通过αβ轴坐标系表示；每个电压矢量都有不同的开关序列，27个输出的电压矢量u(k)总共有3种不同的值，分别为-0.5vdc、0、0.5vdc；
步骤4：在此步骤，主要目标是双环模型预测方法下达成对电压的控制，即分步骤控制中的电压控制步骤，作为双环模型预测的电压环。根据步骤1的电压给定值v(k)和步骤3的27个电压矢量u(k)作差值比较，筛选出27组电压矢量u(k)中离电压给定值v(k)最接近的12组，并根据分步控制电压电流的策略将这12组作为新的逆变器输出电压向量ui(k)，输出电压向量ui(k)的步骤如下：
步骤401：计算27组电压矢量u(k)与电压给定值v(k)的27个不同的差值Δui，其中i∈{1,2,3…27}：
Δui＝|v(k)-u(k)|
步骤402：取最小的12组Δui，并将取得这12组差值Δui的12个u(k)归为新的逆变器输出电压向量ui(k)：
…
其中，
为u1(k)与u(k)的补集，
为u1(k)与u2(k)并集后与u(k)的补集，同理，
为u1(k)与u2(k)…u12(k)并集后与u(k)的补集。
双环模型预测方法下要求分步骤控制输出电压和输出电流，这组ui(k)便作为其中输出电压控制步骤的输出结果。
步骤5：在此步骤，主要目标是双环模型预测方法下达成对电流的控制，可以看作是模型预测控制下的电流环。根据步骤1的k时刻的并网电流采样值i(k)、并网采样电压值e(k)和步骤3的电压矢量u(k)，计算k+1时刻采样电流预测值i(k+1)；
计算电流预测值i(k+1)的方法如下：
S501：根据基尔霍夫电压定律：
式中，L表示三电平逆变器交流侧滤波电感的感值；R表示滤波电感的等效电阻值；u表示三电平逆变器输出的电压矢量，i和e分别表示并网电流矢量和并网电压矢量；
S502：将方程离散化：
整理得：
式中，i(k+1)为预测得到的k+1时刻的电流预测值，i(k)表示步骤1的并网采样电流值，u(k)表示步骤3中逆变器输出的27组电压矢量，e(k)表示步骤1的并网采样电压值，Ts为离散周期，亦是采样周期。
步骤6：在此步骤，根据步骤2的k+1时刻电流参考值i*(k+1)、步骤3的电压矢量u(k)和步骤1的并网电流采样值i(k)，计算k时刻参考电压值u*(k)，即参考电压矢量；参考电压值u*(k(计算方法具体如下：
参考公式：
已知k+1时刻电流参考值为i*(k+1)，用i*(k+1)替换上式中的i(k+1)，则k时刻参考电压值u*(k)为：
步骤7：根据逆变器的中点电流计算k+1时刻的中点电压变化量ΔVNP(k+1)；
根据本发明使用的NPC三电平逆变器，需要考虑逆变器得中点电压，如图1中O点的电压；具体步骤为：
步骤701，根据图1，假设直流母线C1、C2的容值相等，即C1＝C2＝C；
由上式可得中点电压变化量ΔVNP的表达式为
式中Vc1、Vc2满足Vc1+Vc2＝Vdc，可得中点电压变化量ΔVNP与中点电流直接相关；
步骤702，令每相电压的输出电平状态由
表示，1代表P状态，0代表O状态，-1代表N状态；中点电流的表达式为：
io＝(1-|Sa|)·ia+(1-|Sb|)·ib+(1-|Sc|)·ic
＝ia+ib+ic-(|Sa|·ia+|Sb|·ib+|Sc|·ic)
步骤703，对于三相三线制逆变器系统，三相输出电流满足表达式：ia+ib+ic＝0。所以io＝-(|Sa|·ia+|Sb|·ib+|Sc|·ic)，代入上一步可得中点电压变化量ΔVNP的表达式为：
步骤704，假设各相电流在采样周期中保持不变，从而可以将步骤703的公式离散化为：
其中，Sa(k)表示k时刻a相电压的输出电平状态，Sb(k)表示k时刻b相电压的输出电平状态，c(k)表示k时刻c相电压的输出电平状态；ia(k)表示k时刻a相输出电流，ib(k)表示k时刻b相输出电流，ic(k)表示k时刻c相输出电流；C为直流母线电容容值，Ts表示采样周期。
步骤8：在此步骤，根据步骤7计算出的k+1时刻中点电压变化量ΔVNP)k+1)和步骤4计算出的逆变器输出电压向量ui(k)及步骤6计算出的k时刻参考电压值u*(k)，通过价值函数gmin算出最合适的开关序列，即选取价值函数gmin最小时的逆变器输出电压向量ui(k)，ui(k)的开关序列便为此时所采用的开关序列，参考电压值和电压向量差值最小时获得的价值函数gmin最小；具体步骤如下；
步骤801，根据步骤7的中点电压变化量ΔVNP(k+1)、步骤2中的k+1时刻电流参考值i*(k+1)、步骤5中的k+1时刻电流预测值i(k+1)，计算价值函数为：
gmin＝min(i(k+1)-i*(k+1))+λ|ΔVNP(k+1)|
其中，λ为ΔVNP(k+1)权重系数；
步骤802，根据步骤7，i(k+1)-i*(k+1)简化为u(k)-u*(k)，则价值函数为：
gmin＝min(u(k)-u*(k))+λ|ΔVNP(k+1)|
步骤803，根据步骤4，用ui(k)代替u(k)，根据k时刻参考电压值u*(k)，计算价值函数为：
gmin＝min(ui(k)-u*(k))+λ|ΔVNP(k+1)|
步骤804，选取gmin最小时的逆变器输出电压向量ui(k)，其开关序列便为此时所采用的开关序列。
步骤9：依据上一步骤给予逆变器选通脉冲。
下面以三电平逆变器为例，结合附图说明本发明的具体实施方式过程。
1、虚拟同步机实施过程
虚拟同步机VSG控制策略可以让并网逆变器模拟同步发电机组的特性，参与电网的一次调频与一次调压，并且提供虚拟的旋转惯量以帮助电网维持频率稳定。
图5显示，电压型VSG分为有功功率控制部分和无功功率控制部分。其中有功功率部分依据下垂控制的原则对输出电压的频率进行调节，又称为有功-频率控制器；其中无功功率部分对输出电压的幅值进行调节，又称无功-电压控制器。
VSG的有功-频率控制器主要模拟同步发电机的旋转惯量和频率调节。因此首先考虑模拟同步发电机的转子方程，即：
式中，Pe为输出有功功率，Pm为虚拟机械有功功率，ωn为并网逆变器额定角频率，J为VSG的虚拟转动惯量。另外，考虑对同步发电机有功－频率调节的模拟，通常采用P-ω下垂方程，如下式所示：
Pm＝Pset+kp(ωn-ω) (2)
式中，Pm为虚拟机械有功功率，Pset为有功环输出有功功率设定值，kp为P-ω下垂系数，ω为VSG的角速度。综上所述，考虑式(1)、(2)，虚拟同步发电机有功－频率环可以表示为：
式中，
将上式与同步发电机的转子机械运动方程对比可知，下垂系数在Dp微分方程的作用与阻尼系数相同，因此此处的下垂系数能够提供一定虚拟阻尼作用。另外Tset与Te还可以表示为：
在电力系统中，电压的变化通常由系统无功变化引起，因此通常通过对同步发电机励磁的调节，从而实现对电压的调整：—方面，无功负荷的变化会影响线路压降，从而改变系统各个节点处电压。另一方面，电压调节一般采用就近调节，同步机通过检测电压幅值的偏移量，从而对其励磁电压进行调节，最终调整同步发电机的无功输出值，起到对系统电压变化的抑制作用。
模拟上述同步发电机通过调节励磁从而实现调压的过程，VSG的无功-电压控制器可通过检测电压偏移量，从而取改变逆变器无功输出值，通常可通过Q-U下垂方程来实现，如式下：
由上式，Qset为无功设定值，Qm为VSG无功输出值，Dq为Q-U下垂系数。Un为额定输出电压的有效值，U0为输出电压的有效值。同样，为了抑制电压的突变，仿照同步发电机转子方程，引入一个类似结构的微分方程来模拟同步发电机励磁调节系统的电气惯性，如下式所示：
式中，Qm为虚拟无功功率输入量，Qe为输出无功功率。Em为逆变器输出内电势有效值，K为电压方程式的电气惯性参数。综上所述，考虑式(5)、(6)，虚拟同步发电机无功－电压环可以表示为：
经过电压型VSG有功功率控制部分和无功功率控制部分计算后得出双环模型预测控制的电压控制步骤所需的电压给定值v(k)。
根据图4，依据图4的定子电压方程计算出电流给定值i*(k)。因为电流给定值可看作正弦波，则预测k+1时刻的电流参考值i*(k+1)。
2、模型预测控制实施过程
表1详细表述了三电平逆变器开关状态与输出三相电压电平的关系，其拓扑结构如图2所示，P表示正母线电压
电平，N表示负母线电压
电平，直流母线中点0表示参考零电位。共有12个开关，每个桥臂有4个，分为三组，分别为SX1、SX2、SX3和SX4，代表每个桥臂的4个开关，下标x为a、b和c，代表被三个桥臂分为的3组。
表1三电平逆变器开关状态与输出三相电压电平的关系参考表
在理想条件下，当中点电压平衡(即相对于零电位中点电压变化量ΔVNP＝0)时，直流母线电容C1、C2的电压
当中点电压不平衡时即偏离零电位，直流母线电容C1、C2的电压不相等时，中点电压变化量ΔVNP的表达式为
如图2所示，中点电流io与流经上下直流母线电容C1、C2的电流i1、i2的关系如下：令直流母线C1、C2的容值相等，即C1＝C2＝C。
由上式可得中点电压变化量ΔVNP的表达式为
式中Vc1、Vc2满足Vc1+Vc2＝Vdc，可得中点电压变化量ΔVNP与中点电流直接相关。
令每相电压的输出电平状态由
表示，1代表P状态，0代表O状态，-1代表N状态。中点电流的表达式为：
ip＝(1-|Sa|)·ia+(1|Sb|)·ib+(1-|Sc|)·ic
＝ia+ib+ic-(|Sa|·ia+|Sb|·ib+|Sc|·ic) (10)
对于三相三线制逆变器系统，三相输出电流满足表达式：ia+ib+ic＝0
所以io＝-(|Sa|·ia+|Sb|·ib+|Sc|·ic)，代入式(9)可得中点电压变化量ΔVNP的表达式为：
相比于电网电压周期(毫秒级)，系统采样周期(微秒级)十分微小，因此可以合理假设各相电流在采样周期中保持不变，从而可以将式(11)离散化为：
由于每相桥臂有3种不同的相开关状态，所以总共有33＝27种不同的相开关状态组合。为了对应表示三电平逆变器的27个不同的相开关状态组合，定义基本开关状态Sn：
Sn＝[Sa n Sb n Sc n] (13)
其中n＝1,2…27；Si n(i＝a,b,c)为相开关状态。
以逆变器直流侧中点(O)作为参考零点，对应于相开关状态Si n的三电平逆变器输出相电压ui n可以表示为：
其中，Vdc为直流母线总电压值。
因此，在αβ坐标系下可以得到对应于基本开关状态Sn的基本电压矢量un：
图3给出αβ坐标系下的基本电压矢量un的空间分布，un的位置如图中的高亮红点所示，图中的u1,u2,…,u27(依次对应于S1,S2,…,S27)构成FCS-MPC入网电流预测计算所需的基本输入控制集合。根据VSG的电压给定值v(k)和图中的27个电压矢量作差值比较，筛选出27组电压矢量中离电压给定值v(k)最接近的12组向量，并将它们作为新的逆变器输出电压向量ui(k)。由此可见对基本开关状态Sn的评估实质上是对基本电压矢量un的评估。
根据基尔霍夫电压定律，可得：
其中，L(La＝Lb＝Lc＝L)表示三电平逆变器交流侧滤波电感的感值；R(Ra＝Rb＝Rc＝R)表示滤波电感的等效电阻值；u表示三电平逆变器输出的电压矢量，且u∈{u1,u2,…,u27}，i和e分别表示入网电流矢量和电网电压矢量。
因为FCS-MPC需要利用被控变量的离散数学模型进行预测计算，所以本发明采用一阶前向差分方程对式(16)进行离散化后可得：
其中，Ts为离散周期，亦是采样周期。
由式(17)可知，根据k时刻的并网电流采样值i(k)、并网采样电压值e(k)和当前采样周期三电平逆变器所输出的电压矢量u(k)，可以预测计算出k+1时刻的电流预测值i(k+1)。
对于三电平并网逆变器而言，一般的控制目标有：入网功率跟踪控制、入网电流跟踪控制、直流母线电容电压平衡控制、共模电压抑制等。本发明设计价值函数为：
gmin＝min(ui(k)-u*(k))+λ|ΔVNP(k+1)|
ui(k)∈(1,2,3…12)；λ为中点电压的权重系数。
其中u*(k)为参考电压矢量，其计算方式为：
其中i*(k+1)为参考电流矢量，即电流参考值，由VSG电压给定值v(k)经过定子方程后获得。
实施例1
本实施例利用MATLAB中的Simulink工具搭建了三电平逆变电路，直流电经直流母线电容后由三电平电路逆变输出三相电压，经基于三电平变换器的虚拟同步机的双环模型预测控制方法输出光滑的三相正弦电压。仿真过程中的电气参数设置如下表：
具体程序流程如图6所示，通过流程输出合适的12个参考电压。第一步采样k时刻的电压值；第二步，27组向量数组逆变器输出电压us(k)(s∈{1,2,…27})与电压给定值v(k)的27个不同的差值Δui，循环27次；第三步将27次中每次循环中最小的Δui所对应的27组之一的元素us(k)作为新数组ui(k)的元素，同时将最小的Δui赋值为max(Δui)；第四步，跳回第三步，将会得出新的最小的Δui，如此循环13次，将得出的新数组ui(k)(i∈{1,2,3,4…13})作为新的参考电压。
图7(a)为在上述电气参数设置下的输出电流的谐波失真(THD)，图7(b)为在上述电气参数设置下但没有进行电压筛选方法的输出电流的谐波失真(THD)。在未使用双环模型预测方法时，THD为0.83％，有较大的谐波，而加入双环模型预测方法时，THD为0.39％，说明本发明所用的方法可以有效的改善输出谐波分量，能够有效的稳定输出，可见本发明的控制方法控制效果显著。
以上内容是结合具体的实施方式对本发明所作的进一步详细说明，不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干简单推演或替换，都应当视为属于本发明的保护范围。
Claims (9)
1.基于虚拟同步机的三电平并网逆变器的双环模型预测方法，其特征在于：步骤包括：
步骤1：虚拟同步机根据下垂控制原理输出k时刻的电压给定值v(k)，并通过定子电压方程得出电流给定值i*(k)，同时采集得出并网电流采样值i(k)和并网电压采样值e(k)；
步骤2：通过k时刻的电流给定值i*(k)预测出k+1时刻电流参考值i*(k+1)；
步骤3：获取三电平变换器的输出电压并分为P组不同状态的电压矢量u(k)，每个电压矢量都有不同的开关序列；
步骤4：根据电压给定值v(k)和P个电压矢量u(k)作差值比较，筛选出P组电压矢量u(k)中离电压给定值v(k)最接近的Q组，并将这Q组作为新的逆变器输出电压向量ui(k)，作为双环模型预测的电压环，实现分步骤控制中的电压控制；
步骤5：根据k时刻的并网电流采样值i(k)、并网采样电压值e(k)和电压矢量u(k)，计算k+1时刻采样电流预测值i(k+1)，作为模型预测控制下的电流环，实现分步骤控制中的电流控制；
步骤6：根据k+1时刻电流参考值i*(k+1)、电压矢量u(k)和并网电流采样值i(k)，计算k时刻参考电压值u*(k)；
步骤7：根据逆变器的中点电流计算k+1时刻的中点电压变化量ΔVNP(k+1)；
步骤8：根据k+1时刻中点电压变化量ΔVNP(k+1)、逆变器输出电压向量ui(k)及k时刻参考电压值u*(k)，通过价值函数gmin计算出参考电压值和电压向量差值最小时的价值函数，选取最小价值函数时的逆变器输出电压向量ui(k)，并输出ui(k)的开关序列；
步骤9：依据步骤8的开关序列给予逆变器选通脉冲。
2.根据权利要求1所述的基于虚拟同步机的三电平并网逆变器的双环模型预测方法，其特征在于：步骤1中测出电流给定值i*(k)，具体步骤如下：
步骤101：通过电压型虚拟同步机的有功功率控制部分计算出输出的电压给定值v(k)的相角θ，通过电压型虚拟同步机的无功功率控制部分计算出输出的电压给定值v(k)的电压幅值E，得出v(k)在αβ轴坐标系下的数学表达式；
步骤102：将输出的电压给定值v(k)通过定子电压方程转换为线电流即电流给定值i*(k)。
3.根据权利要求1所述的基于虚拟同步机的三电平并网逆变器的双环模型预测方法，其特征在于：步骤2的预测计算方法为：假定电流给定值和电压给定值是恒定的正弦波，根据正弦波在k时刻的值预测出k+1时刻的值。
4.根据权利要求1所述的基于虚拟同步机的三电平并网逆变器的双环模型预测方法，其特征在于：步骤3中的电压矢量u(k)的P组为27组，27个输出的电压矢量u(k)总共有3种不同的值，分别为-0.5vdc、0、0.5vdc。
5.根据权利要求1所述的基于虚拟同步机的三电平并网逆变器的双环模型预测方法，其特征在于：步骤4中，Q组为12组；输出电压向量ui(k)的步骤如下：
步骤401：计算27组电压矢量u(k)与电压给定值v(k)的27个不同的差值Δui，其中i∈{1，2，3…27}：
Δui＝|v(k)-u(k)|
步骤402：取最小的12组Δui，并将取得这12组差值Δui的12个u(k)归为新的逆变器输出电压向量ui(k)：
其中，
为u1(k)与u(k)的补集，
为u1(k)与u2(k)并集后与u(k)的补集，同理，
为u1(k)与u2(k)...u12(k)并集后与u(k)的补集。
双环模型预测方法下要求分步骤控制输出电压和输出电流，这组ui(k)便作为其中输出电压控制步骤的输出结果。
6.根据权利要求1所述的基于虚拟同步机的三电平并网逆变器的双环模型预测方法，其特征在于：步骤5中，计算电流预测值i(k+1)的方法如下：
式中，i(k+1)为预测得到的k+1时刻的电流预测值，i(k)表示步骤1的并网采样电流值，u(k)表示步骤3中逆变器输出的电压矢量，e(k)表示步骤1的并网采样电压值，Ts为离散周期即采样周期，L为三电平逆变器交流侧滤波电感的感值，R表示滤波电感的等效电阻值。
7.根据权利要求1所述的基于虚拟同步机的三电平并网逆变器的双环模型预测方法，其特征在于：步骤6中，k时刻参考电压值u*(k)计算方法具体如下：
其中，i*(k+1)为步骤2的k+1时刻电流参考值、u(k)为步骤3的电压矢量、i(k)为步骤1的并网电流采样值。
8.根据权利要求1所述的基于虚拟同步机的三电平并网逆变器的双环模型预测方法，其特征在于：步骤7中，k+1时刻的中点电压变化量ΔVNP(k+1)计算方法为：
其中，三电平逆变器每相电压的输出电平状态由
表示；其中，Sa(k)表示k时刻a相电压的输出电平状态，Sb(k)表示k时刻b相电压的输出电平状态，Sc(k)表示k时刻c相电压的输出电平状态；ia(k)表示k时刻a相输出电流，ib(k)表示k时刻b相输出电流，ic(k)表示k时刻c相输出电流；C为直流母线电容容值，Ts表示采样周期。
9.根据权利要求1所述的基于虚拟同步机的三电平并网逆变器的双环模型预测方法，其特征在于：步骤8具体包括：
步骤801，根据步骤7的中点电压变化量ΔVNP(k+1)、步骤2中的k+1时刻电流参考值i*(k+1)、步骤5中的k+1时刻电流预测值i(k+1)，计算价值函数为：
gmin＝min(i(k+1)-i*(k+1))+λ|ΔVNP(k+1)|
其中，λ为ΔVNP(k+1)权重系数；
步骤802，根据步骤7，i(k+1)-i*(k+1)简化为u(k)-u*(k)，则价值函数为：
gmin＝min(u(k)-u*(k))+λ|ΔVNP(k+1)|
步骤803，根据步骤4，用ui(k)代替u(k)，根据k时刻参考电压值u*(k)，计算价值函数为：
gmin＝min(ui(k)-u*(k))+λ|ΔVNP(k+1)|
步骤804，选取gmin最小时的逆变器输出电压向量ui(k)，其开关序列便为此时所采用的开关序列。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210573684.9A CN115021602A (zh) | 2022-05-25 | 2022-05-25 | 基于虚拟同步机的三电平并网逆变器的双环模型预测方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210573684.9A CN115021602A (zh) | 2022-05-25 | 2022-05-25 | 基于虚拟同步机的三电平并网逆变器的双环模型预测方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN115021602A true CN115021602A (zh) | 2022-09-06 |
Family
ID=83069254
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210573684.9A Withdrawn CN115021602A (zh) | 2022-05-25 | 2022-05-25 | 基于虚拟同步机的三电平并网逆变器的双环模型预测方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115021602A (zh) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115579872B (zh) * | 2022-10-13 | 2023-07-25 | 深圳市禾望电气股份有限公司 | 一种虚拟同步机电流跟踪控制方法 |
| CN117239711A (zh) * | 2023-11-13 | 2023-12-15 | 四川大学 | 改善抽油机井群供电质量的储能控制方法和装置 |
| CN117614300A (zh) * | 2024-01-23 | 2024-02-27 | 苏州大学 | T型三电平三相四桥臂逆变器连续模型预测控制方法 |
| CN117833248A (zh) * | 2024-03-06 | 2024-04-05 | 电子科技大学 | 一种t型三电平并联有源电力滤波器无模型预测控制方法 |
-
2022
- 2022-05-25 CN CN202210573684.9A patent/CN115021602A/zh not_active Withdrawn
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115579872B (zh) * | 2022-10-13 | 2023-07-25 | 深圳市禾望电气股份有限公司 | 一种虚拟同步机电流跟踪控制方法 |
| CN117239711A (zh) * | 2023-11-13 | 2023-12-15 | 四川大学 | 改善抽油机井群供电质量的储能控制方法和装置 |
| CN117239711B (zh) * | 2023-11-13 | 2024-02-02 | 四川大学 | 改善抽油机井群供电质量的储能控制方法和装置 |
| CN117614300A (zh) * | 2024-01-23 | 2024-02-27 | 苏州大学 | T型三电平三相四桥臂逆变器连续模型预测控制方法 |
| CN117614300B (zh) * | 2024-01-23 | 2024-04-05 | 苏州大学 | T型三电平三相四桥臂逆变器连续模型预测控制方法 |
| CN117833248A (zh) * | 2024-03-06 | 2024-04-05 | 电子科技大学 | 一种t型三电平并联有源电力滤波器无模型预测控制方法 |
| CN117833248B (zh) * | 2024-03-06 | 2024-05-10 | 电子科技大学 | 一种t型三电平并联有源电力滤波器无模型预测控制方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN115021602A (zh) | 基于虚拟同步机的三电平并网逆变器的双环模型预测方法 | |
| Song et al. | Predictive duty cycle control of three-phase active-front-end rectifiers | |
| CN110867887B (zh) | Npc三电平逆变器最优序列模型预测控制装置及方法 | |
| Xia et al. | Multi-objective optimal model predictive control for three-level ANPC grid-connected inverter | |
| CN112701710B (zh) | 基于交流恒功率负荷动态性能的储能变换器稳定控制方法 | |
| Pranith et al. | Improved Gaussian filter based solar PV-BES microgrid with PLL based islanding detection and seamless transfer control | |
| Wu et al. | Model predictive direct power control for modular multilevel converter under unbalanced conditions with power compensation and circulating current reduction | |
| CN115864878A (zh) | 一种模块化多电平换流器控制方法 | |
| CN107634657B (zh) | 一种矩阵变换器的预测控制方法及装置 | |
| CN113437894B (zh) | 一种三相八开关逆变器的共模电压抑制方法 | |
| CN113162115B (zh) | 一种三相并网逆变器加权滑模模型预测电流控制方法 | |
| CN114567022A (zh) | 一种非理想电网条件下pwm整流器优化预测功率控制方法 | |
| CN114785166A (zh) | 基于滑模控制的t型整流器三矢量模型预测控制方法 | |
| CN114928261A (zh) | 三相五电平pwm整流器的模型预测及零序电压平衡控制方法 | |
| CN110707949B (zh) | 定频pwm整流器的控制方法 | |
| Zhang et al. | An improved robust model predictive and repetitive combined control for three-phase four-leg active power filters with fixed switching frequency | |
| CN113452069A (zh) | 一种向无源网络供电的逆变站多步模型预测控制方法 | |
| CN114256874B (zh) | 一种t型逆变器快速有限集模型预测控制方法 | |
| Gu et al. | Control Strategy of Three-Level NPC Inverter Based on Variable Coefficient Virtual Vector Model Predictive Control | |
| Lamterkati et al. | Direct Power Fuzzy Control of Three-Phase AC-DC Converter using a Space Vector Modulation | |
| He et al. | Control strategy of high-performance active rectifier using modular multilevel converter without DC-link voltage sensors | |
| CN113890073B (zh) | 不等容的多逆变器并联飞轮储能系统改进模型预测方法 | |
| CN109950927B (zh) | 一种模块化多电平变流器的功率补偿策略 | |
| CN113765118A (zh) | 一种基于预测控制的vsg延时补偿方法 | |
| Draou et al. | Performance analysis of advanced static VAr compensator using three-level IGBT inverter |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| WW01 | Invention patent application withdrawn after publication | ||
| WW01 | Invention patent application withdrawn after publication |
Application publication date: 20220906 |