CN108920856B - 一种改进式模块化多电平换流器等效方法 - Google Patents

一种改进式模块化多电平换流器等效方法 Download PDF

Info

Publication number
CN108920856B
CN108920856B CN201810765310.0A CN201810765310A CN108920856B CN 108920856 B CN108920856 B CN 108920856B CN 201810765310 A CN201810765310 A CN 201810765310A CN 108920856 B CN108920856 B CN 108920856B
Authority
CN
China
Prior art keywords
bridge
equivalent
arm
modular multilevel
multilevel converter
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
CN201810765310.0A
Other languages
English (en)
Other versions
CN108920856A (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.)
State Grid Corp of China SGCC
State Grid Shanghai Electric Power Co Ltd
Northeast Electric Power University
Original Assignee
State Grid Corp of China SGCC
Northeast Dianli University
State Grid Shanghai Electric Power 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 State Grid Corp of China SGCC, Northeast Dianli University, State Grid Shanghai Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN201810765310.0A priority Critical patent/CN108920856B/zh
Publication of CN108920856A publication Critical patent/CN108920856A/zh
Application granted granted Critical
Publication of CN108920856B publication Critical patent/CN108920856B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Inverter Devices (AREA)

Abstract

本发明是一种改进式模块化多电平换流器等效方法,通过区分系统正常运行与启动或故障时的半桥型子模块工作状态,对系统处于正常与处于启动阶段或发生故障时的桥臂进行戴维南等效,并综合考虑闭锁情况和非闭锁情况,提出了改进式模块化多电平换流器等效方法,有别于只考虑系统正常工作状态的桥臂等效方法。针对半桥型子模块闭锁时,桥臂电流方向频繁改变,二极管将面临数值插值问题,从而导致桥臂电流出现畸变点,本发明在半桥型子模块闭锁时,利用PSCAD仿真软件中自带的二极管进行建模,使二极管数值插值问题得到了有效解决,提高了仿真精度,且提速效果明显。本发明具有科学合理,计算简单,简便易行,适用性强,仿真精度高,提速效果好的优点。

Description

一种改进式模块化多电平换流器等效方法
技术领域
本发明涉及模块化多电平换流器技术领域,是一种改进式模块化多电平换流器等效方法。
背景技术
随着模块化多电平换流器高压直流输电(ModularMultilevelConverterbasedHighVoltageDirectCurrent,MMC-HVDC)工程向更高电压等级、更大传输功率发展,所需的半桥型子模块数目也不断增加,在大规模MMC-HVDC系统仿真时,需要对大量超高阶矩阵求逆,使得求解过程非常缓慢、耗费大量时间,不利于在实际工程中的应用,因此有必要对详细模型进行等效。
现有的等效方法都是建立在MMC稳态状态下,未考虑半桥型子模块闭锁的情况,而在电磁暂态仿真中,当模块化多电平换流器处于启动阶段或发生故障时,会有许多子模块处于闭锁状态,用现有的等效方法处理处于闭锁状态下的半桥型子模块时,由于二极管插值模拟不当而导致桥臂电流出现畸变点,致使仿真精度低,效果差。
发明内容
本发明所要解决的技术问题是,克服现有等效方法的不足,提供一种科学合理,综合考虑半桥型子模块闭锁情况和非闭锁情况并解决二极管数值插值问题,适用性强,仿真精度高,提速效果好的改进式模块化多电平换流器等效方法。
1.一种改进式模块化多电平换流器等效方法,其特征是,它包括如下步骤:
1)根据半桥型子模块中IGBT器件K1和IGBT器件K2的开通情况进行区分,得到半桥型子模块当前工作状态;
半桥型子模块的三种工作状态:第一种为投入状态,半桥型子模块中IGBT器件K1导通、IGBT器件K2闭锁;第二种为切出状态,半桥型子模块中IGBT器件K1闭锁、IGBT器件K2导通;第三种为闭锁状态,半桥型子模块中IGBT器件K1和IGBT器件K2均闭锁;当系统中MMC正常工作时,半桥型子模块的工作状态包括投入状态和切出状态;当系统中MMC启动或故障时,半桥型子模块中IGBT处于闭锁状态;
2)当确定系统中模块化多电平换流器处于正常工作状态,含有N个半桥型子模块和一个电感串联组成的桥臂用戴维南进行等值,等效步骤为:
(1)半桥型子模块中IGBT器件K1和续流二极管用可变电阻R1等效代替,半桥型子模块中IGBT器件K2和续流二极管用可变电阻R2来等效代替,当SPWM调制输出的信号为高电平时,触发半桥型子模块中IGBT开通,用一个数值较小的电阻,取值0.001Ω来代替;当SPWM调制输出的信号为低电平时,触发半桥型子模块中IGBT关断,用一个数值较大的电阻,取值108Ω来代替;
可变电阻R1和R2由式(1)给出:
R=KiRT+(1-Ki)Ro (1)
其中,RT是通态电阻,Ro是阻断电阻,Ki是开关i的开关信号,为1时开关导通,为0时开关关断;
(2)对半桥型子模块中电容C进行Dommel电磁暂态等值计算,过程如下:
其暂态过程用电磁感应定律来描述,即:
Figure BDA0001728890160000021
其中,ic(t)为t时刻电容电流,C为电容值,uc(t)为t时刻电容电压;
式(2)积分形式为:
Figure BDA0001728890160000022
其中,ΔT为仿真步长,uc(t-ΔT)为t-ΔT时刻的电容电压;
应用梯形积分法则,将式(3)表示成式(4)的形式;
Figure BDA0001728890160000023
其中,ic(t-ΔT)为t-ΔT时刻电容电流;
将式(4)改写为式(5)、(6)的形式;
Figure BDA0001728890160000024
Figure BDA0001728890160000025
其中,Rc=ΔT/(2C)为电容C暂态等值电阻,Ic(t-ΔT)为t-ΔT时刻电容的等值电流源;
由式(6)可知,Ic(t-ΔT)与t-ΔT时刻的电流和电压值有关;
t-ΔT时刻的电容电流为:
Figure BDA0001728890160000031
将式(7)代入式(6),得到t-ΔT时刻电容的等值电流源Ic(t-ΔT)为:
Figure BDA0001728890160000032
对比式(6)和式(8),可知式(8)中不包含ic(t-ΔT),表明Ic(t-ΔT)可直接利用t-2ΔT时刻电容的等值电流源Ic(t-2ΔT)进行递推计算,使计算过程简单化,计算速度加快;
(3)通过戴维南定理计算获得半桥型子模块的等效电压Ueqsm和等效电阻Reqsm为:
Figure BDA0001728890160000033
(4)对桥臂中电感L进行Dommel电磁暂态等值计算,过程如下:
t-ΔT时刻的电感电流和电感的电流源递推公式为:
Figure BDA0001728890160000034
Figure BDA0001728890160000035
其中,iL(t)为t时刻电感电流,RL=(2L)/ΔT为电感L暂态等值电阻,IL(t-ΔT)为t-ΔT时刻电感的等值电流源,iL(t-ΔT)为t-ΔT时刻电感电流,uL(t-ΔT)为t-ΔT时刻电感电压;
与上述求解电容的等值电流源方法相同,可得电感的等值电流源递推公式为:
Figure BDA0001728890160000036
其中,IL(t-2ΔT)为t-2ΔT时刻电感的等值电流源;
(5)获得由N个半桥型子模块和一个电感串联组成的桥臂的戴维南等效电路,等效电路中等效电压∑Ueq(t)和等效电阻∑Req为:
Figure BDA0001728890160000037
Figure BDA0001728890160000041
其中,Ueqi(t)为桥臂中第i个半桥型子模块的戴维南等效电压,Reqi为桥臂中第i个半桥型子模块的戴维南等效电阻,i=1,2,3,···,N,N为桥臂中半桥型子模块数目,其值可由式(9)获得;
3)当确定系统中模块化多电平换流器处于故障或启动状态,由N个半桥型子模块和一个电感串联组成的桥臂将等效为一个具有等效参数的半桥不控子模块;
(1)半桥型子模块中的IGBT闭锁时,即半桥型子模块中IGBT处于断开的状态,由于每个桥臂的所有半桥型子模块是串联关系,因此N个半桥型子模块流过相同的电流,均等于桥臂电流iarm,当iarm>0,所有半桥型子模块中的IGBT器件K1的反并联二极管D1导通;当iarm<0,所有半桥型子模块的IGBT器件K2的反并联二极管D2导通;
(2)闭锁后的桥臂电压uarm(t)为:
Figure BDA0001728890160000042
其中,ucj(t)为桥臂中第j个半桥型子模块的电容电压,j=1,2,3,···,N,N为桥臂中半桥型子模块数目;iarm(t)为桥臂电流,Ron为二极管的通态电阻;
定义:
Figure BDA0001728890160000043
Ron∑=NRon (17)
其中,Rc=ΔT/(2C)为电容C暂态等值电阻,ucjeq(t-ΔT)为t-ΔT时刻桥臂中第j个半桥型子模块的电容等值电压源,Ron∑为桥臂中二极管的通态电阻之和;
式(14)可改写为:
Figure BDA0001728890160000044
在闭锁状态下,桥臂戴维南等效电阻Rarm和戴维南等效电压uarmeq(t)根据桥臂电流方向获得:
当iarm(t)≥0时,Rarm=NRc+Ron∑,
Figure BDA0001728890160000051
当iarm(t)≤0时,Rarm=Ron∑,uarmeq(t)=0;
(3)由于二极管是单向导通元件,当半桥型子模块中IGBT闭锁时,桥臂电流方向会频繁改变,采用PSCAD仿真软件中自带的二极管进行建模,利用自带的插值算法避免了二极管数值插值的问题,提高了仿真波形的精度;
4)综合考虑闭锁情况和非闭锁情况,模块化多电平换流器等效电压∑Ueq为:
Figure BDA0001728890160000052
模块化多电平换流器等效电阻∑Req为:
Figure BDA0001728890160000053
本发明是一种改进式模块化多电平换流器等效方法,通过区分系统正常运行与启动或故障时的半桥型子模块工作状态,分别对系统处于正常工作状态与处于启动阶段或发生故障时的桥臂进行戴维南等效,并综合考虑闭锁情况和非闭锁情况,提出了改进式模块化多电平换流器等效方法,有别于只考虑系统正常工作状态的传统桥臂等效方法。针对半桥型子模块闭锁时,桥臂电流方向会频繁改变,采用传统等效方法二极管将会面临数值插值,从而导致桥臂电流出现畸变点的问题,本发明在半桥型子模块处于闭锁情况时,利用PSCAD仿真软件中自带的二极管进行建模,使二极管数值插值的问题得到了有效的解决,从而提高了仿真波形的精度,且提速效果明显。本发明具有科学合理,计算简单,简便易行,适用性强,仿真精度高,提速效果好的优点。
附图说明
图1为模块化多电平换流器(MMC)的通用拓扑结构图;
图2为模块化多电平换流器(MMC)中半桥型子模块的拓扑结构图;
图3为工程中常用的半桥型子模块的工作状态图;半桥型子模块工作状态包括:投入、切出和闭锁;
图4为非闭锁情况下的半桥型子模块的戴维南等效电路图;半桥型子模块的电容和电感采用Dommel电磁暂态数值计算;
图5为电容元件进行等值计算前的电路图;
图6为电容元件的等值计算模型图;
图7为电感元件进行等值计算前的电路图;
图8为电感元件的等值计算模型图;
图9为模块化多电平换流器非闭锁情况下的桥臂戴维南等值图;将由N个半桥型子模块和一个电感串联组成的桥臂用戴维南进行等值;
图10为模块化多电平换流器闭锁且iarm>0情况下的桥臂电路结构图,所有半桥型子模块的D1导通;
图11为模块化多电平换流器闭锁且iarm<0情况下的桥臂电路结构图,所有半桥型子模块的D2导通;
图12为闭锁情况下的半桥型子模块的戴维南等效电路图;其中使用PSCAD仿真软件中自带的二极管进行搭建,避免桥臂电流出现畸变点;
图13为改进式模块化多电平换流器等效方法的桥臂戴维南等值图,考虑闭锁情况和非闭锁情况;
图14为模块化多电平换流器启动仿真的详细模型的直流电压波形图;
图15为模块化多电平换流器启动仿真的详细模型的桥臂电流波形图;
图16为模块化多电平换流器启动仿真的传统等效模型的直流电压波形图;
图17为模块化多电平换流器启动仿真的传统等效模型的桥臂电流波形图;
图18为模块化多电平换流器启动仿真的改进等效模型的直流电压波形图;
图19为模块化多电平换流器启动仿真的改进等效模型的桥臂电流波形图;
图20为仿真精度对比的模块化多电平换流器A相上桥臂稳态电流波形图;
图21为仿真精度对比的模块化多电平换流器A相电压稳态波形波形图;
图22为仿真精度对比的模块化多电平换流器A相上桥臂暂态电流波形图;
图23为仿真精度对比的模块化多电平换流器A相电压暂态波形波形图;
图24为改进式模块化多电平换流器等效方法的模型与详细模型仿真运行时间对比图;
图25为改进式模块化多电平换流器等效方法的模型的提速倍数图。
具体实施方式
下面将结合附图和实例对本发明的改进式模块化多电平换流器等效方法进行清楚、完整地描述。
本发明基于半桥型子模块的模块化多电平换流器(ModularMultilevelConverter,MMC),其通用拓扑结构如图1所示,由三相六桥臂构成,每个桥臂由N个半桥型子模块(Half-BridgeSub-Module,HBSM)及电抗器L相串联构成,半桥型子模块拓扑结构如图2所示,包括1个电容器、两个反向二极管和2个绝缘栅双极晶闸管(Insulated-gateBipolarTransistor,IGBT或IGT)。
S1:根据半桥型子模块工作状态确定模块化多电平换流器工作状态;
步骤S1定义半桥型子模块的工作状态,确定模块化多电平换流器在正常工作状态和处于故障或启动时的半桥型子模块工作状态,具体包括:
S11:利用半桥型子模块中IGBT器件K1和IGBT器件K2的开通情况定义半桥型子模块三种工作状态,如图3所示,半桥型子模块的三种工作状态:第一种为投入状态,半桥型子模块中IGBT器件K1导通,IGBT器件K2闭锁;第二种为切出状态,半桥型子模块中IGBT器件K1闭锁,IGBT器件K2导通;第三种为闭锁状态,半桥型子模块中IGBT器件K1和IGBT器件K2均闭锁。
S12:确定系统中模块化多电平换流器不同状态下的半桥型子模块工作状态:第一种为系统中模块化多电平换流器处于正常工作状态时,半桥型子模块工作状态包括投入状态和切出状态;第二种为系统中模块化多电平换流器处于故障或启动状态时,半桥型子模块中的IGBT处于闭锁状态;
S2:当确定系统中模块化多电平换流器处于正常工作状态,将由N个半桥型子模块和一个电感串联组成的桥臂采用戴维南定理进行等效,等效过程如图4所示。
步骤S2包括对半桥型子模块中的两组IGBT和续流二极管进行等效,对半桥型子模块中电容C的Dommel电磁暂态等值计算,通过戴维南定理获得半桥型子模块的等效电压Ueqsm和等效电阻Reqsm,对桥臂中电感L进行Dommel电磁暂态等值计算,得到系统中模块化多电平换流器处于正常工作状态下,即非闭锁情况下的桥臂的戴维南等效电路,具体包括:
S21:对半桥型子模块中的两组IGBT和续流二极管进行等效;
(6)本发明将半桥型子模块中的两组IGBT和续流二极管用具有两种状态的可变电阻R1、R2来等效代替,当SPWM调制输出的信号为高电平时,触发半桥型子模块中IGBT开通,用一个数值较小的电阻,取值0.001Ω来代替;当SPWM调制输出的信号为低电平时,触发半桥型子模块中IGBT关断,用一个数值较大的电阻,取值108Ω来代替;
可变电阻R1和R2由式(1)给出:
R=KiRT+(1-Ki)Ro (1)
其中,RT是通态电阻,Ro是阻断电阻,Ki是开关i的开关信号,为1时开关导通,为0时开关关断。
S22:对半桥型子模块中电容C进行Dommel电磁暂态等值计算,电容元件进行等值计算前的电路如图5所示;
过程如下:
其暂态过程用电磁感应定律来描述,即
Figure BDA0001728890160000081
其中,ic(t)为t时刻电容电流,C为电容值,uc(t)为t时刻电容电压;
式(2)积分形式为
Figure BDA0001728890160000082
其中,ΔT为仿真步长,uc(t-ΔT)为t-ΔT时刻的电容电压;
应用梯形积分法则,将式(3)表示成式(4)的形式;
Figure BDA0001728890160000083
其中,ic(t-ΔT)为t-ΔT时刻电容电流;
将式(4)改写为式(5)、(6)的形式;
Figure BDA0001728890160000084
Figure BDA0001728890160000085
其中,Rc=ΔT/(2C)为电容C暂态等值电阻,Ic(t-ΔT)为t-ΔT时刻电容的等值电流源;所得电容元件的等值电路如图6所示;
由式(6)可知,Ic(t-ΔT)与t-ΔT时刻的电流和电压值有关;
t-ΔT时刻的电容电流为
Figure BDA0001728890160000091
将式(7)代入式(6),得到t-ΔT时刻电容的等值电流源Ic(t-ΔT)为:
Figure BDA0001728890160000092
对比式(6)和式(8),可知式(8)中不包含ic(t-ΔT),表明Ic(t-ΔT)可直接利用t-2ΔT时刻电容的等值电流源Ic(t-2ΔT)进行递推计算,使计算过程简单化,计算速度加快。
S23:通过戴维南定理计算获得半桥型子模块的等效电压Ueqsm和等效电阻Reqsm为:
Figure BDA0001728890160000093
S24:对桥臂中电感L进行Dommel电磁暂态等值计算,电感元件进行等值计算前的电路如图7所示;
t-ΔT时刻的电感电流和电感的电流源递推公式为:
Figure BDA0001728890160000094
Figure BDA0001728890160000095
其中,iL(t)为t时刻电感电流,RL=(2L)/ΔT为电感L暂态等值电阻,IL(t-ΔT)为t-ΔT时刻电感的等值电流源,iL(t-ΔT)为t-ΔT时刻电感电流,uL(t-ΔT)为t-ΔT时刻电感电压,所得电感元件的等值电路如图8所示;
与上述求解电容的等值电流源方法相同,可得电感的等值电流源递推公式为:
Figure BDA0001728890160000096
其中,IL(t-2ΔT)为t-2ΔT时刻电感的等值电流源。
S25:得到系统中模块化多电平换流器处于正常工作状态,即非闭锁情况下桥臂的戴维南等效电路,如图9所示;
获得由N个半桥型子模块和一个电感串联组成的桥臂的戴维南等效电路,等效电路中等效电压∑Ueq(t)和等效电阻∑Req为:
Figure BDA0001728890160000101
Figure BDA0001728890160000102
其中,Ueqi(t)为桥臂中第i个半桥型子模块的戴维南等效电压,Reqi为桥臂中第i个半桥型子模块的戴维南等效电阻,i=1,2,3,···,N,N为桥臂中半桥型子模块数目,其值可由式(9)获得。
S3:当确定系统中模块化多电平换流器处于故障或启动状态,由N个半桥型子模块和一个电感串联组成的桥臂将等效为一个具有等效参数的半桥不控子模块,其桥臂等值电路如图12所示,并提出二极管数值插值问题的解决方法;
半桥型子模块闭锁时,IGBT处于断开的状态,由于每个桥臂的所有半桥型子模块是串联关系,因此N个半桥型子模块流过相同的电流,均等于桥臂电流iarm,当iarm>0,所有半桥型子模块中的IGBT器件K1的反并联二极管D1导通;当iarm<0,所有半桥型子模块中的IGBT器件K2的反并联二极管D2导通;
步骤S3包括计算闭锁状态下桥臂戴维南等效电阻Rarm和戴维南等效电压uarmeq(t),提出二极管数值插值问题的解决方法,具体包括:
S31:计算闭锁状态下桥臂戴维南等效电阻Rarm和戴维南等效电压uarmeq(t);
闭锁后的桥臂电压uarm(t)为:
Figure BDA0001728890160000103
其中,ucj(t)为桥臂中第j个半桥型子模块的电容电压,j=1,2,3,···,N,N为桥臂中半桥型子模块数目;iarm(t)为桥臂电流,Ron为二极管的通态电阻;
定义:
Figure BDA0001728890160000111
Ron∑=NRon (17)
其中,Rc=ΔT/(2C)为电容C暂态等值电阻,ucjeq(t-ΔT)为t-ΔT时刻桥臂中第j个半桥型子模块的电容等值电压源,Ron∑为桥臂中二极管的通态电阻之和;
式(14)可以改写为:
Figure BDA0001728890160000112
在闭锁状态下,桥臂戴维南等效电阻Rarm和戴维南等效电压uarmeq(t)根据桥臂电流方向获得:
当iarm(t)≥0时,Rarm=NRc+Ron∑,
Figure BDA0001728890160000113
当iarm(t)≤0时,Rarm=Ron∑,uarmeq(t)=0。
S32:提出二极管数值插值的问题的解决方法;
由于二极管是单向导通元件,当半桥型子模块中IGBT闭锁时,桥臂电流方向会频繁改变,采用传统的等效方法,二极管将会面临数值插值的问题,导致桥臂电流出现畸变点,因此采用PSCAD仿真软件中自带的二极管进行建模,利用自带的插值算法以避免数值插值的问题,从而提高了仿真波形的精度。
S4:综合考虑闭锁情况和非闭锁情况,对模块化多电平换流器桥臂进行戴维南等效,获得适用于正常工作状态和故障或启动工作状态的等效模型,即基于改进式模块化多电平换流器等效方法的模块化多电平换流器模型,并采用PSCAD仿真软件中自带的二极管进行建模以解决二极管数值插值的问题,其戴维南等效电路如图13所示;
基于改进式模块化多电平换流器等效方法的模块化多电平换流器等效电压∑Ueq和等效电阻∑Req为:
Figure BDA0001728890160000121
Figure BDA0001728890160000122
S5:在PSCAD仿真软件中分别搭建了5电平的模块化多电平换流器详细模型、模块化多电平换流器传统等效模型和本发明所提出的基于改进式模块化多电平换流器等效方法的模块化多电平换流器模型,进行仿真验证;
具体参数如下:半桥型子模块数为4,交流侧电压为5.5kV,半桥型子模块电容C=1900μH,桥臂电感L=2mH,开关频率为2kHz。
S51:模块化多电平换流器启动仿真验证;
PSCAD仿真软件中搭建模块化多电平换流器详细模型、模块化多电平换流器传统等效模型和本发明所提出的基于改进式模块化多电平换流器等效方法的模块化多电平换流器模型,分别进行换流器启动仿真,仿真步长为20μs,仿真时间为0.2s;
模块化多电平换流器详细模型、模块化多电平换流器传统等效模型和本发明所提出的基于改进式模块化多电平换流器等效方法的模块化多电平换流器模型的直流电压波形分别如图14、图16和图18所示,由图16可以看出,模块化多电平换流器传统等效模型的直流电压出现畸变点,而由图14和图16中可以看出基于改进式模块化多电平换流器等效方法的模块化多电平换流器模型中的直流电压波形与模块化多电平换流器详细模型中的直流电压波形吻合度很高,未出现明显的畸变点,验证了本发明提出的改进式模块化多电平换流器等效方法有效地解决了半桥型子模块发生闭锁后的插值问题;
模块化多电平换流器详细模型、模块化多电平换流器传统等效模型和本发明所提出的基于改进式模块化多电平换流器等效方法的模块化多电平换流器模型的桥臂电流波形分别如图15、图16和图19所示,由图15和图19的桥臂电流可以看出,模块化多电平换流启动时,基于改进式模块化多电平换流器等效方法的模块化多电平换流器模型和模块化多电平换流器详细模型充电时桥臂电流波形具有很高的吻合度,验证了该改进式模块化多电平换流器等效方法具有较高的仿真精度。
S52:仿真精度验证;
图20为稳态情况下模块化多电平换流器A相上桥臂的稳态电流,图21为稳态情况下模块化多电平换流器A相输出电压波形,由稳态仿真波形图可以看出,本发明提出的基于改进式模块化多电平换流器等效方法的模块化多电平换流器模型与搭建的模块化多电平换流器详细模型仿真图符合完好,两条曲线接近重合,验证了基于改进式模块化多电平换流器等效方法的模块化多电平换流器模型在稳态情况下具有较高的仿真精度;
在1.3s时直流线路发生故障,故障持续0.6s恢复正常,图22为故障期间模块化多电平换流器A相上桥臂电流,图23为故障期间模块化多电平换流器A相输出电压波形。由暂态仿真波形图可以看出,本发明所提出的基于改进式模块化多电平换流器等效方法的模块化多电平换流器模型与搭建的模块化多电平换流器详细模型仿真图也符合完好,验证了该模型在暂态情况下同样具有较高的精度。
S53:仿真速度验证;
本发明的仿真验证在Windows7操作系统下运行,处理器为3.4GHz双核CPU,安装内存为4GB,PSCAD/EMTDC的版本为专业版4.5,在PSCAD中分别搭建了5至31电平模块化多电平换流器详细模型和本发明所提出的基于改进式模块化多电平换流器等效方法的模块化多电平换流器模型,以验证改进式模块化多电平换流器等效方法的提速效果,仿真步长取20us,仿真时间2s。
表1为基于改进式模块化多电平换流器等效方法的模块化多电平换流器模型和模块化多电平换流器详细模型的运行时间;
表1基于改进式模块化多电平换流器等效方法的模块化多电平换流器模型和模块化多电平换流器详细模型的运行时间对比表
Figure BDA0001728890160000131
图24为模块化多电平换流器详细模型和基于改进式模块化多电平换流器等效方法的模块化多电平换流器模型运行时间的对比图,由图中看出,随着模块化多电平换流器半桥型子模块数的增加,模块化多电平换流器详细模型的仿真时间以指数的形式增加,而基于改进式模块化多电平换流器等效方法的模块化多电平换流器模型以线性形式增长。
图25为基于改进式模块化多电平换流器等效方法的模块化多电平换流器模型的提速倍数,从图中可以看出随着半桥型子模块数的增多,基于改进式模块化多电平换流器等效方法的模块化多电平换流器模型的提速倍数越来越高,因此本发明提出的改进式模块化多电平换流器等效方法在保证仿真精度较高的基础上,能显著提高仿真速度,验证了改进式模块化多电平换流器等效方法的可行性与有效性。
显然,上述实施例仅仅是为清楚地说明所作的举例,而并非对实施方式的限定。对于所属领域的普通技术人员来说,在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。而由此所引伸出的显而易见的变化或变动仍处于本发明创造的保护范围之中。

Claims (1)

1.一种改进式模块化多电平换流器等效方法,其特征是,它包括如下步骤:
1)根据半桥型子模块中IGBT器件K1和IGBT器件K2的开通情况进行区分,得到半桥型子模块当前工作状态;
半桥型子模块的三种工作状态:第一种为投入状态,半桥型子模块中IGBT器件K1导通、IGBT器件K2闭锁;第二种为切出状态,半桥型子模块中IGBT器件K1闭锁、IGBT器件K2导通;第三种为闭锁状态,半桥型子模块中IGBT器件K1和IGBT器件K2均闭锁;当系统中MMC正常工作时,半桥型子模块的工作状态包括投入状态和切出状态;当系统中MMC启动或故障时,半桥型子模块中IGBT处于闭锁状态;
2)当确定系统中模块化多电平换流器处于正常工作状态,含有N个半桥型子模块和一个电感串联组成的桥臂用戴维南进行等值,等效步骤为:
(1)半桥型子模块中IGBT器件K1和续流二极管用可变电阻R1等效代替,半桥型子模块中IGBT器件K2和续流二极管用可变电阻R2来等效代替,当SPWM调制输出的信号为高电平时,触发半桥型子模块中IGBT开通,用一个数值较小的电阻,取值0.001Ω来代替;当SPWM调制输出的信号为低电平时,触发半桥型子模块中IGBT关断,用一个数值较大的电阻,取值108Ω来代替;
可变电阻R1和R2由式(1)给出:
R=KiRT+(1-Ki)Ro (1)
其中,RT是通态电阻,Ro是阻断电阻,Ki是开关i的开关信号,为1时开关导通,为0时开关关断;
(2)对半桥型子模块中电容C进行Dommel电磁暂态等值计算,过程如下:
其暂态过程用电磁感应定律来描述,即:
Figure FDA0001728890150000011
其中,ic(t)为t时刻电容电流,C为电容值,uc(t)为t时刻电容电压;
式(2)积分形式为:
Figure FDA0001728890150000012
其中,ΔT为仿真步长,uc(t-ΔT)为t-ΔT时刻的电容电压;
应用梯形积分法则,将式(3)表示成式(4)的形式;
Figure FDA0001728890150000021
其中,ic(t-ΔT)为t-ΔT时刻电容电流;
将式(4)改写为式(5)、(6)的形式;
Figure FDA0001728890150000022
Figure FDA0001728890150000023
其中,Rc=ΔT/(2C)为电容C暂态等值电阻,Ic(t-ΔT)为t-ΔT时刻电容的等值电流源;
由式(6)可知,Ic(t-ΔT)与t-ΔT时刻的电流和电压值有关;
t-ΔT时刻的电容电流为:
Figure FDA0001728890150000024
将式(7)代入式(6),得到t-ΔT时刻电容的等值电流源Ic(t-ΔT)为:
Figure FDA0001728890150000025
对比式(6)和式(8),可知式(8)中不包含ic(t-ΔT),表明Ic(t-ΔT)可直接利用t-2ΔT时刻电容的等值电流源Ic(t-2ΔT)进行递推计算,使计算过程简单化,计算速度加快;
(3)通过戴维南定理计算获得半桥型子模块的等效电压Ueqsm和等效电阻Reqsm为:
Figure FDA0001728890150000026
(4)对桥臂中电感L进行Dommel电磁暂态等值计算,过程如下:
t-ΔT时刻的电感电流和电感的电流源递推公式为:
Figure FDA0001728890150000027
Figure FDA0001728890150000028
其中,iL(t)为t时刻电感电流,RL=(2L)/ΔT为电感L暂态等值电阻,IL(t-ΔT)为t-ΔT时刻电感的等值电流源,iL(t-ΔT)为t-ΔT时刻电感电流,uL(t-ΔT)为t-ΔT时刻电感电压;
与上述求解电容的等值电流源方法相同,可得电感的等值电流源递推公式为:
Figure FDA0001728890150000031
其中,IL(t-2ΔT)为t-2ΔT时刻电感的等值电流源;
(5)获得由N个半桥型子模块和一个电感串联组成的桥臂的戴维南等效电路,等效电路中等效电压∑Ueq(t)和等效电阻∑Req为:
Figure FDA0001728890150000032
Figure FDA0001728890150000033
其中,Ueqi(t)为桥臂中第i个半桥型子模块的戴维南等效电压,Reqi为桥臂中第i个半桥型子模块的戴维南等效电阻,i=1,2,3,···,N,N为桥臂中半桥型子模块数目,其值可由式(9)获得;
3)当确定系统中模块化多电平换流器处于故障或启动状态,由N个半桥型子模块和一个电感串联组成的桥臂将等效为一个具有等效参数的半桥不控子模块;
(1)半桥型子模块中的IGBT闭锁时,即半桥型子模块中IGBT处于断开的状态,由于每个桥臂的所有半桥型子模块是串联关系,因此N个半桥型子模块流过相同的电流,均等于桥臂电流iarm,当iarm>0,所有半桥型子模块中的IGBT器件K1的反并联二极管D1导通;当iarm<0,所有半桥型子模块的IGBT器件K2的反并联二极管D2导通;
(2)闭锁后的桥臂电压uarm(t)为:
Figure FDA0001728890150000034
其中,ucj(t)为桥臂中第j个半桥型子模块的电容电压,j=1,2,3,···,N,N为桥臂中半桥型子模块数目;iarm(t)为桥臂电流,Ron为二极管的通态电阻;
定义:
Figure FDA0001728890150000041
Ron∑=NRon (17)
其中,Rc=ΔT/(2C)为电容C暂态等值电阻,ucjeq(t-ΔT)为t-ΔT时刻桥臂中第j个半桥型子模块的电容等值电压源,Ron∑为桥臂中二极管的通态电阻之和;
式(14)可改写为:
Figure FDA0001728890150000042
在闭锁状态下,桥臂戴维南等效电阻Rarm和戴维南等效电压uarmeq(t)根据桥臂电流方向获得:
当iarm(t)≥0时,Rarm=NRc+Ron∑
Figure FDA0001728890150000045
当iarm(t)≤0时,Rarm=Ron∑,uarmeq(t)=0;
(3)由于二极管是单向导通元件,当半桥型子模块中IGBT闭锁时,桥臂电流方向会频繁改变,采用PSCAD仿真软件中自带的二极管进行建模,利用自带的插值算法避免了二极管数值插值的问题,提高了仿真波形的精度;
4)综合考虑闭锁情况和非闭锁情况,模块化多电平换流器等效电压∑Ueq为:
Figure FDA0001728890150000043
模块化多电平换流器等效电阻∑Req为:
Figure FDA0001728890150000044
CN201810765310.0A 2018-07-12 2018-07-12 一种改进式模块化多电平换流器等效方法 Active CN108920856B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810765310.0A CN108920856B (zh) 2018-07-12 2018-07-12 一种改进式模块化多电平换流器等效方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810765310.0A CN108920856B (zh) 2018-07-12 2018-07-12 一种改进式模块化多电平换流器等效方法

Publications (2)

Publication Number Publication Date
CN108920856A CN108920856A (zh) 2018-11-30
CN108920856B true CN108920856B (zh) 2022-11-08

Family

ID=64411432

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810765310.0A Active CN108920856B (zh) 2018-07-12 2018-07-12 一种改进式模块化多电平换流器等效方法

Country Status (1)

Country Link
CN (1) CN108920856B (zh)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110108986A (zh) * 2019-05-29 2019-08-09 华北电力大学 直流电网故障电流复频域计算方法
CN111177932B (zh) * 2019-12-31 2021-04-27 天津大学 一种用于提升柔性变电站电磁暂态仿真速度的建模方法
CN111382550B (zh) * 2020-03-10 2023-04-18 南方电网科学研究院有限责任公司 模块化多电平换流器的动态组合实时仿真方法及使用方法
CN111682523B (zh) * 2020-05-21 2022-02-01 中国南方电网有限责任公司超高压输电公司检修试验中心 考虑桥臂闭锁时二极管续流作用的mmc快速仿真模型和方法
CN111753413B (zh) * 2020-06-15 2024-05-14 国网经济技术研究院有限公司 混合式制动电阻变换器的电磁暂态等效建模方法和系统
CN112417667B (zh) * 2020-11-16 2024-02-02 国网浙江省电力有限公司电力科学研究院 一种基于mmc高效电磁暂态桥臂等效模型的仿真方法
CN113281678B (zh) * 2021-04-12 2022-11-08 西安理工大学 模块化多电平变换器半桥子模块上管开路故障的定位方法
CN113705144B (zh) * 2021-07-28 2024-06-11 国网冀北电力有限公司经济技术研究院 一种考虑闭锁的半桥型mmc电磁暂态仿真方法
CN115808646B (zh) * 2023-01-17 2023-05-02 国网山西省电力公司电力科学研究院 基于改进半桥型mmc的柔性直流系统线路故障性质识别方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103593520A (zh) * 2013-10-31 2014-02-19 南方电网科学研究院有限责任公司 一种模块化多电平换流器的等值仿真计算建模方法
CN106055752A (zh) * 2016-05-23 2016-10-26 华北电力大学 一种提高mmc高效电磁暂态模型仿真精度的方法
CN106787890A (zh) * 2017-01-13 2017-05-31 中国电力科学研究院 一种cdsm‑mmc桥臂电磁暂态建模方法和装置
CN108229021A (zh) * 2018-01-03 2018-06-29 华北电力大学 基于实时数字仿真器的模块化多电平换流器闭锁建模方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103593520A (zh) * 2013-10-31 2014-02-19 南方电网科学研究院有限责任公司 一种模块化多电平换流器的等值仿真计算建模方法
CN106055752A (zh) * 2016-05-23 2016-10-26 华北电力大学 一种提高mmc高效电磁暂态模型仿真精度的方法
CN106787890A (zh) * 2017-01-13 2017-05-31 中国电力科学研究院 一种cdsm‑mmc桥臂电磁暂态建模方法和装置
CN108229021A (zh) * 2018-01-03 2018-06-29 华北电力大学 基于实时数字仿真器的模块化多电平换流器闭锁建模方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
"Modeling, control, and reduced-order representation of modular multilevel";Andres E. Leona等;《Electric Power Systems Research》;20180627;第163卷;196-210页 *
"一种二极管箝位型MMC的高效等值建模方法";周月宾等;《中国电机工程学报》;20160405;第36卷(第7期);1925-1932页 *
"模块化多电平柔性直流输电数字_模拟混合实时仿真技术";刘栋等;《电力自动化设备》;20130228;第33卷(第2期);68-80页 *

Also Published As

Publication number Publication date
CN108920856A (zh) 2018-11-30

Similar Documents

Publication Publication Date Title
CN108920856B (zh) 一种改进式模块化多电平换流器等效方法
Yaragatti Design and implementation of active neutral-point-clamped nine-level reduced device count inverter: an application to grid integrated renewable energy sources.
Oates et al. A comparison of two methods of estimating losses in the modular multi-level converter
Nishida et al. A predictive instantaneous-current PWM controlled rectifier with AC-side harmonic current reduction
CN112464443B (zh) 一种用于电力电子变流器igbt结温波动的计算方法
Li et al. Loss calculation method and loss characteristic analysis of MMC based VSC-HVDC system
CN108229021B (zh) 基于实时数字仿真器的模块化多电平换流器闭锁建模方法
CN105356780B (zh) 子模块混合型模块化多电平变流器的调制方法及系统
CN108631632B (zh) 基于虚拟桥臂数学模型的mmc瞬时功率损耗计算方法
CN112134472A (zh) 基于mmc换流器的双端系统直流侧谐振控制方法及系统
CN109428340B (zh) 一种柔性直流输电装置的仿真方法及系统
CN110429839B (zh) 一种三相电压型pwm整流器的分数阶建模方法
CN108540042A (zh) 空间矢量脉冲宽度幅值调制方法
CN115296554A (zh) 一种高调制比混合式mmc及其控制方法
CN108376992B (zh) 一种混合子模块mmc等效仿真方法
Lee et al. High performance current controller for sparse matrix converter based on model predictive control
CN112016043B (zh) 模块化多电平换流器的稳态故障电流的计算方法
CN107863781B (zh) 一种模块化多电平换流阀损耗确定方法及简化方法
CN109256972B (zh) 一种基于五段式五电平变换器svpwm调制方法
CN112001142A (zh) 一种半桥型模块化多电平换流器的实时仿真方法
Bendre et al. Modeling and design of a neutral point regulator for a three level diode clamped rectifier
CN115425861A (zh) 多电平逆变电路的控制方法和并网逆变器
CN112886840B (zh) 一种模块化多电平变换器损耗优化控制方法
CN112052638B (zh) 一种全桥型模块化多电平换流器的实时仿真方法
CN114447943A (zh) 静止无功发生器限流值的确定方法及装置

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
GR01 Patent grant
GR01 Patent grant