CN104993715A - 模块化多电平换流器阀损快速评估算法 - Google Patents
模块化多电平换流器阀损快速评估算法 Download PDFInfo
- Publication number
- CN104993715A CN104993715A CN201510365972.5A CN201510365972A CN104993715A CN 104993715 A CN104993715 A CN 104993715A CN 201510365972 A CN201510365972 A CN 201510365972A CN 104993715 A CN104993715 A CN 104993715A
- Authority
- CN
- China
- Prior art keywords
- igbt
- diode
- loss
- centerdot
- submodule
- 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.)
- Granted
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
- Y02E60/60—Arrangements for transfer of electric power between AC networks or generators via a high voltage DC link [HVCD]
Landscapes
- Power Conversion In General (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
本发明公开了一种模块化多电平换流器阀损快速评估算法,步骤如下:一、根据系统参数和控制方法,建立MMC的快速仿真模型;二、将仿真结果中的子模块电容电压、开关器件触发脉冲和桥臂电流保存为可读取数据文件;三、根据MMC子模块中开关器件的特性参数建立开关器件损耗计算和PN结温度估计的数学模型;四、将数据文件导入损耗计算程序,计算每个数据步长下的单个IGBT和二极管的阀损耗;五、根据IGBT和二极管的阀损耗功率与结温的关系,通过迭代计算进而叠加得到MMC的阀损耗功率。本发明解决了MMC子模块电压、电流以及换流阀工作温度无法通过解析方式精确计算的关键问题,可以广泛应用于柔性直流输电技术领域。
Description
技术领域
本发明涉及电力系统柔性输配电技术领域,特别是涉及一种模块化多电平换流器阀损快速评估算法。
背景技术
目前,模块化多电平换流器(MMC,Modular Multilevel Converter)在高压柔性直流输电领域占据主导地位。相比于采用PWM(pulsewidth modulation,脉冲宽度调制)调制的两电平或三电平换流器,MMC避免了IGBT(Insulated Gate Bipolar Transistor,绝缘栅双极晶体管)直接串联引起的一致触发等问题,适用于各种高压直流输电场合。MMC随着子模块串联数即电平数的增加,能够显著提升电压等级和功率容量,而且交流侧电压正弦度更高、损耗更低。损耗不仅是直流系统的一个重要评价指标,对于换流阀散热器设计和器件选型也起着关键的作用。因此需要寻找一种快捷且有效的方法来完成MMC系统阀损耗的评估。
对于采用PWM调制的两电平换流器,稳态运行下的直流电压波动很小,其开关动作时机和桥臂电流都可以用解析的方法准确表达。而MMC采用子模块串联构成桥臂的结构,在阶梯波调制(通常为最近电平调制NLC)下,子模块电容电压会出现较大范围的波动,需要设计额外的子模块均压控制方法。均压控制会带来额外的开关动作,而且这些动作时机难以用解析方法表达。故业界一般选用基于电磁暂态仿真的方法精确计算MMC阀损耗。IEC正在制定IEC 62751-1“高压直流输电系统的电压源型换流器阀损耗计算”标准。即推荐采用MMC电磁暂态仿真结果进行准确的阀损耗计算。现有技术一般采用详细电磁暂态仿真结果计算MMC的阀损。但MMC包含的子模块数量众多,采用详细电磁暂态仿真结果计算MMC的阀损耗存在着仿真速度慢从而导致阀损评估结果慢的缺陷。
发明内容
本发明的目的是为了克服上述背景技术的不足,提供一种模块化多电平换流器阀损快速评估算法,解决了MMC子模块电压、电流以及换流阀工作温度无法通过解析方式精确计算的关键问题,相较于原有的纯解析计算方法,具有更高的准确性和更强的适应性。
本发明提供的一种模块化多电平换流器阀损快速评估算法,包括以下步骤:步骤一、根据系统参数和控制方法,在电磁暂态仿真软件中建立MMC的快速仿真模型;步骤二、该快速仿真模型将仿真结果中的子模块电容电压、开关器件触发脉冲和桥臂电流保存为可读取数据文件;步骤三、根据MMC子模块中开关器件的特性参数,所述MMC子模块中开关器件即为IGBT及其反并联二极管,建立开关器件损耗计算和PN结温度估计的数学模型;步骤四、将仿真结果的数据文件导入损耗计算程序,计算出每个数据步长下的单个IGBT和单个二极管的阀损耗;步骤五、根据IGBT和二极管的阀损耗功率与结温的关系,通过迭代计算,得到单个IGBT和单个二极管修正的阀损耗功率,进而叠加得到MMC的阀损耗功率。
在上述技术方案中,所述步骤一中,MMC快速仿真模型的子模块电容电压动态简化公式表示为:
在上述技术方案中,所述步骤一中,MMC快速仿真模型的桥臂为等效受控电压源结构,其控制电压值瞬时值Uup(t)为:
其中,USMi(t)为第i个子模块的瞬时端口电压,Uconi(t)为开关管的导通压降,N为单个桥臂上串联子模块的数量,UCi(t)是子模块电容电压,SCi(t)是开关函数,取值0或1分别表示子模块切除或投入。
在上述技术方案中,所述步骤三中的损耗计算数学模型如下:
PTcon(t)=Uce(t)·IT(t)=(RT(TjT(t))·IT(t)+Uce0(TjT(t)))·IT(t)
PDcon(t)=UD(t)·ID(t)=(RD(TjD(t))×ID(t)+UD0(TjD(t)))·ID(t)
式中,IT(t)、ID(t)分别代表流经IGBT和二极管的瞬时电流,由仿真结果经处理得到;PTcon(t)、PTsw(t)和PToff(t)分别为IGBT导通损耗、开关损耗和关断损耗;PDcon(t)、PDsw(t)和PDoff(t)分别为二极管导通损耗、开关损耗和关断损耗;Uce0(t)、UD(t)分别为t时刻IGBT和二极管的导通压降;TjT(t)、TjD(t)分别为t时刻IGBT和二极管的结温,其初值选取为器件底板温度Tcase;RT是IGBT导通电阻,Uce0是IGBT擎住电压;Roff_T,Roff_D分别为IGBT和二极管的截止电阻,T0为采样间隔,N0为采样间隔时间内器件开关次数;Eon(j)、Eoff(j)和Erec(j)分别为IGBT单次开能量、IGBT单次关能量和二极管单次关断能量,j为采样间隔内开关次数的序号。
在上述技术方案中,所述步骤三中,所述的IGBT单次开能量Eon(j)、单次关能量Eoff(j)和二极管单次关断能量Erec(j)分别为:
式中,a1、a2、a3、b1、b2、b3、c1、c2、c3是根据器件表开关能量曲线经二次拟合得到的系数,UceN、UDN为IGBT和二极管的额定电压,ion_T(t)、ioff_T(t)为IGBT开关电流瞬时值,ioff_D(t)为二极管关断电流瞬时值,USM(t)为子模块瞬时电压,ρTon、ρToff、ρDrec为结温修正系数。
在上述技术方案中,所述步骤三中,所述的结温修正系数ρTon、ρToff、ρDrec计算方法为:
式中,Eon1、Eon2分别为IGBT在125°和25°时的单次开通能量;Eoff1、Eoff2分别为IGBT在125°和25°时的单次关断能量;Erec1、Erec2分别为二极管在125°和25°时的单次关断能量,TjT、TjD分别为IGBT和二极管的结温。
在上述技术方案中,所述步骤五中,IGBT和二极管结温通过下式求得:
TjT(t+Δt)=PT(TjT(t))·(Zth(JC_T)+Zth(CS))+TS
TjD(t+Δt)=PD(TjD(t))·(Zth(JC_D)+Zth(CS))+TS,
式中,PT、PD分别为IGBT和二极管总的阀损,Zth(JC_T)、Zth(JC_D)为IGBT和底板、二极管和底板间的温阻;Zth(CS)为底板和散热器间的温阻,TS为散热器温度。
本发明模块化多电平换流器阀损快速评估算法,具有以下有益效果:本发明搭建的适用于全面评估MMC阀损耗的快速仿真模型,能够保留损耗计算所需的子模块电压电流以及触发脉冲,在与原模型基本等效的条件下,简化原有模型以达到提高仿真速度的目的,使得本发明具有较强的实用性。
MMC阀损耗按如下计算方案获得:
一、建立能大幅提高仿真速度、适用于损耗计算的快速模型。快速模型编写自定义模块,采用梯形积分法简化电容充放电过程,以桥臂受控电压源替代串联子模块结构以提高仿真速度。仿真获得换流站每相的桥臂电流Iarm、子模块电压USM以及触发脉冲S。
二、将PSCAD(全称Power Systems Computer Aided Design,是世界上广泛使用的电磁暂态仿真软件)得到的仿真数据输出到MATLAB(矩阵工厂或矩阵实验室,是美国MathWorks公司出品的商业数学软件)中进行处理,得到流经各器件的电流以及开关时刻对应的电流瞬时值,结合仿真数据共同作为阀损耗计算的输入。
三、根据子模块中IGBT及其反并联二极管的器件参数表,建立IGBT、二极管的损耗模型,利用插值法获得对应温度下的参数值。
四、根据实际器件工作环境,建立器件-底板-散热器的热电路。
五、根据损耗功率与器件温度的关系迭代计算,得到各站阀损耗、器件温度及损耗率。
本发明采用录入仿真数据计算阀损耗的方式,主要有以下技术优势:
1、本发明与仿真实例贴合度高,误差较解析算法要小得多;
2、本发明无需修改计算方法,自适应各种控制策略、调制方式和拓扑结构下的阀损耗计算;
3、本发明得出了每个子模块上器件的阀损耗及工作温度,方便器件选型及损耗分布评估。
附图说明
图1为本发明模块化多电平换流器阀损快速评估算法的流程示意图;
图2为双极双端MMC-MTDC仿真电路图;
图3为采用实际开关器件布置如图1的电路结构图;
图4为本发明模块化多电平换流器阀损快速评估算法的步骤一中采用受控电压源等效的快速模型示意图;
图5为本发明中的快速模型与现有技术的详细模型仿真结果对比图;
图6为本发明模块化多电平换流器阀损快速评估算法的步骤二中子模块结构示意图;
图7为本发明模块化多电平换流器阀损快速评估算法的步骤二中子模块桥臂电流及触发脉冲波形示意图;
图8为本发明模块化多电平换流器阀损快速评估算法的步骤二中子模块上臂电流及标记的开关时间点示意图;
图9为本发明模块化多电平换流器阀损快速评估算法的步骤三中UCE-IC伏安特性曲线拟合示意图;
图10为本发明模块化多电平换流器阀损快速评估算法的步骤三中开关能量曲线拟合示意图;
图11为本发明模块化多电平换流器阀损快速评估算法的步骤五中开关器件结温、开关器件底板和散热器间的热路图;
图12为本发明送端正极换流站实施例中不同工况下损耗分布示意图;
图13为本发明送端正极换流站实施例中不同工况下结温变化曲线示意图。
具体实施方式
下面结合附图及实施例对本发明作进一步的详细描述,但该实施例不应理解为对本发明的限制。
参见图1,本发明模块化多电平换流器阀损快速评估算法,包括以下步骤:
步骤一、根据系统参数和控制方法,在电磁暂态仿真软件中建立MMC的快速仿真模型:
所述步骤一中,MMC快速仿真模型的子模块电容电压动态简化公式表示为:
其中:△t是积分步长,i为子模块序号,UCi(t)是子模块电容电压,UCi(t-△t)是前一积分步长的电压值,SCi(t)是开关函数,取值0或1分别表示子模块切除或投入,Iarm代表每一相的上桥臂或下桥臂电流,C为电容值;
另外,MMC快速仿真模型的桥臂为等效受控电压源结构,其控制电压值瞬时值Uup(t)为:
其中,USMi(t)为第i个子模块的瞬时端口电压,Uconi(t)为开关管的导通压降,N为单个桥臂上串联子模块的数量,UCi(t)是子模块电容电压,SCi(t)是开关函数,取值0或1分别表示子模块切除或投入;
步骤二、该快速仿真模型将仿真结果中的子模块电容电压、开关器件触发脉冲和桥臂电流保存为可读取数据文件;
步骤三、根据MMC子模块中开关器件的特性参数,所述MMC子模块中开关器件即为IGBT及其反并联二极管,建立开关器件损耗计算和PN结温度估计的数学模型:
所述损耗计算数学模型如下:
PTcon(t)=Uce(t)·IT(t)=(RT(TjT(t))·IT(t)+Uce0(TjT(t)))·IT(t)
PDcon(t)=UD(t)·ID(t)=(RD(TjD(t))×ID(t)+UD0(TjD(t)))·ID(t)
式中,IT(t)、ID(t)分别代表流经IGBT和二极管的瞬时电流,由仿真结果经处理得到;PTcon(t)、PTsw(t)和PToff(t)分别为IGBT导通损耗、开关损耗和关断损耗;PDcon(t)、PDsw(t)和PDoff(t)分别为二极管导通损耗、开关损耗和关断损耗;Uce0(t)、UD(t)分别为t时刻IGBT和二极管的导通压降;TjT(t)、TjD(t)分别为t时刻IGBT和二极管的结温,其初值选取为器件底板温度Tcase;RT是IGBT导通电阻,Uce0是IGBT擎住电压;Roff_T,Roff_D分别为IGBT和二极管的截止电阻,T0为采样间隔,N0为采样间隔时间内器件开关次数;Eon(j)、Eoff(j)和Erec(j)分别为IGBT单次开能量、IGBT单次关能量和二极管单次关断能量,j为采样间隔内开关次数的序号;
上式中,所述的IGBT单次开能量Eon(j)、单次关能量Eoff(j)和二极管单次关断能量Erec(j)分别为:
式中,a1、a2、a3、b1、b2、b3、c1、c2、c3是根据器件表开关能量曲线经二次拟合得到的系数,UceN、UDN为IGBT和二极管的额定电压,ion_T(t)、ioff_T(t)为IGBT开关电流瞬时值,ioff_D(t)为二极管关断电流瞬时值,USM(t)为子模块瞬时电压,ρTon、ρToff、ρDrec为结温修正系数;
进一步地,所述的结温修正系数ρTon、ρToff、ρDrec计算方法为:
式中,Eon1、Eon2分别为IGBT在125°和25°时的单次开通能量;Eoff1、Eoff2分别为IGBT在125°和25°时的单次关断能量;Erec1、Erec2分别为二极管在125°和25°时的单次关断能量,TjT、TjD分别为IGBT和二极管的结温;
步骤四、将仿真结果的数据文件导入损耗计算程序,计算出每个数据步长下的单个IGBT和单个二极管的阀损耗;
步骤五、根据IGBT和二极管的阀损耗功率与结温的关系,通过迭代计算,得到单个IGBT和单个二极管修正的阀损耗功率,进而叠加得到MMC的阀损耗功率:
IGBT和二极管结温通过下式求得:
TjT(t+Δt)=PT(TjT(t))·(Zth(JC_T)+Zth(CS))+TS
TjD(t+Δt)=PD(TjD(t))·(Zth(JC_D)+Zth(CS))+TS
式中,PT、PD分别为IGBT和二极管总的阀损,Zth(JC_T)、Zth(JC_D)为IGBT和底板、二极管和底板间的温阻;Zth(CS)为底板和散热器间的温阻,TS为散热器温度。
实施例
以下结合附图及实施例,对本发明进行进一步详细说明。
本发明提出了适用于MMC-MTDC(Modular Multilevel Converterbased multi-terminal high voltage direct current,MMC-MTDC,基于模块化多电平变流器的多端直流输电系统)系统快速阀损计算方法,步骤如下:
步骤1:根据实际MMC-MTDC系统的拓扑结构,在PSCAD/EMTDC(Power Systems Computer Aided Design/Electromagnetic Transients including DC)平台上搭建与之等效的快速仿真模型。选取双端双极性MMC-MTDC系统计算其阀损耗分布,系统拓扑结构如图2所示。
双端双极系统共四个MMC,每个MMC采用三相六桥臂结构(如图3所示)。对于含实际开关器件的详细模型中各桥臂,将其中的串联子模块部分用受控电压源代替,为保留子模块电容电压及电流开关动态,搭建自定义模块模拟调制过程如图4。
首先计算电容电压动态,采用梯形积分法等效电容动态过程:
其中,Δt是积分步长,UCi(t)是子模块电容电压,UCi(t-Δt)是前一积分步长的电压值,N为桥臂子模块个数。SCi(t)是开关函数,取值0或1表示子模块切除或投入,ISMi(t)表示流入子模块的电流,Iarm(t)代表每相上桥臂或下桥臂电流。
由子模块电流计算开关管压降Uconi(t):
Uconi(t)=Rconi·ISMi(t)+U0i,(i=1,2,…N) (2)
式中,Rconi为开关管导通电阻,Uoi是开关管的擎住电压。
根据电流方向和投入切除状态,Uconi(t)共有四种状态:
式中,U0i是开关管的擎住电压,Rconi代表开关管的导通电阻,其下标T1、T2分别代表位于子模块上、下臂的IGBT;D1、D2分别代表位于子模块上、下臂的二极管。据此求解子模块交流侧输出电压USM:
USMi(t)=Uconi(t)+SCi(t)·UCi(t),(i=1,2,…N) (4)
将桥臂上串联子模块的输出电压叠加得到桥臂电压,以每相上(下)桥臂电压Uup(t)为例:
编写自定义模块实现上述过程,将该电压值作为桥臂受控电压源的指令值,即可构造一个等效数字桥臂,进而得到MMC的快速计算模型,与详细模型对比如图5。仿真完成后,以数据流的形式保存自定义模块输出子模块电压、桥臂电流及触发脉冲的实时数据。
步骤2:将PSCAD得到的仿真数据流输出到MATLAB中处理。
时域仿真结果通过.out文件保存在对应算例的.emt文件夹中。仿真结果得到双端双极系统4个换流站每相的桥臂电流Iarm(t)、子模块电压USM(t)以及触发脉冲SCi(t)。子模块结构如图6,其触发脉冲表示为:
S1与S2分别为子模块中IGBT和二极管的触发脉冲,在开通时呈现互补的状态,桥臂电流Iarm(t)流过IGBT和二极管,计其正负分别为结合触发状态和桥臂电流正负计算各元件电流,如图7所示:
式中iT1(t)、iD1(t)、iT2(t)、iD2(t)分别为留过子模块上、下臂的IGBT和二极管对应瞬时电流。由仿真结果得到开关动作发生的时间及对应时刻的电流瞬时值(以IGBT为例,设IGBT开通及关断时刻为ton_T、toff_T,对应的电流瞬时值为ion_T(t)、ioff_T(t)(如图8所示):
步骤3:设置器件初始温度计算各部分损耗:
1)计算IGBT和二极管的导通损耗PTcon(t)、PDcon(t):
PTcon(t)=Uce(t)·IT(t)=(RT(TjT)·IT(t)+Uce0(TjT))·IT(t) (9)
PDcon(t)=UD(t)·ID(t)=(RD(TjD)×ID(t)+UD0(TjD))·ID(t) (10)
式中,结温TjT、TjD初值选取为器件底板温度Tcase。RT是IGBT导通电阻,Uce0是IGBT擎住电压,参数值由厂家提供的器件特性曲线(图9)拟合得到,都与器件PN结温度有关。由插值法得到对应结温下的参数值:
式中Uce0、UD0分别为IGBT、二极管擎住电压,RT、RD分别为IGBT、二极管导通电阻,角标1、2代表125°和25°时的器件参数值,由器件参数表给出。
2)计算截止损耗:IGBT和二极管在关断状态下由于漏电流引起的损耗,由子模块电压和截止电阻按如下式求取:
式中,Roff_T,Roff_D分别为IGBT和二极管的截止电阻,其取值在器件手册中给定,典型值为4×104~2×105欧姆。
3)计算开关损耗:计算器件单次开关能量,根据步骤2中得到的开关时刻电流值ion_T、ioff_T、ioff_D及子模块电压USM得到:
式中,a、b、c是根据器件表开关能量曲线经二次拟合得到的系数(如图10),UceN、UDN为IGBT和二极管的额定电压,ρTon、ρToff、ρDrec为结温修正系数,利用125°和25°时的器件参数值(角标为1、2)插值计算:
步骤4:叠加单位时间内的开关能量计算开关损耗功率:
式中N0为单位时间内开关次数。
计算IGBT、二极管器件总损耗:
PT=PTsw+PTcon+PToff
(16)
PD=PDsw+PDcon+PDoff
式中PTcon(t)、PTsw(t)和PToff(t)分别为IGBT导通损耗、开关损耗和关断损耗;PDcon(t)、PDsw(t)和PDoff(t)为二极管导通损耗、开关损耗和关断损耗;
步骤5:修正器件温度值:
根据计算出的器件损耗,结合图11所示的等效热电路结构,修正各器件工作温度:
TjT(t+Δt)=PT(TjT(t))·(Zth(JC_T)+Zth(CS))+TS
(17)
TjD(t+Δt)=PD(TjD(t))·(Zth(JC_D)+Zth(CS))+TS
式中Zth(JC_T)、Zth(JC_D)为IGBT和底板、二极管和底板间的温阻;Zth(CS)为底板和散热器间的温阻,TS为散热器温度。将修正温度TjT(t+Δt)、TjD(t+Δt)替代原始温度TjT(t)、TjD(t),重复步骤3和步骤4,直到两者相等,完成迭代过程确定PT(t)、PD(t)终值。
计算双端双极MMC-MTDC各换流站损耗及损耗率:
计算每个子模块的损耗功率PSM(t):
式中j=1,2代表半桥型子模块的上下臂。
每相桥臂上子模块个数为N,整个换流器阀损耗功率PMMC_loss(t)为:
式中Pup_loss,k、Pdn_loss,k为每相上、下桥臂总损耗,k取A、B、C三相。由此计算换流站损耗率η:
η=PMMC_loss(t)/PMMC(t) (21)
式中PMMC(t)为换流站输入有功功率。
为验证此方法的实效性,基于双极双端MMC-MTDC系统运行数据,在MATLAB中编写代码实现上述步骤,计算流程如图1所示。表1中给出了系统运行参数,表2列出了换流阀器件参数及开关能量拟合参数。
表1 换流站系统参数
系统参数 | 送端彭厝换流站 | 受端湖边换流站 |
交流系统额定电压 | 220kV | 220kV |
交流系统标称频率 | 50Hz | 50Hz |
换流变压器 | 230±8*1.25%/167kV | 230±8*1.25%/167kV |
直流额定电压 | ±320kV | ±320kV |
额定直流功率 | 500MW×2 | 500MW×2 |
桥臂额定直流电流 | 594A | 594A |
桥臂额定交流电流 | 985A | 985A |
单个桥臂子模块数 | 200 | 200 |
子模块电容值 | 10mF | 10mF |
桥臂电抗器 | 60mH | 90mH |
直流电抗器 | 50mH | 100mH |
调制比工作范围 | 0.7-0.95 | 0.7-0.95 |
表2 IGBT模块主要参数
IGBT型号 | 5SNA 1200E330100 |
关断电阻(Ohm) | 1000000 |
IGBT开通电阻(Ohm) | 0.0019 |
二极管开通电阻(Ohm) | 0.0011 |
IGBT导通压降(kV) | 0.00159 |
二极管导通压降(kV) | 0.00138 |
IGBT-BASE温阻(K/W) | 0.0085 |
Diode-BASE温阻(K/W) | 0.017 |
IGBT模块额定电压(kV) | 3.3 |
IGBT单次开关能量(125度)(mJ) | 3840 |
IGBT单次开关能量(25度)(mJ) | 2760 |
Diode单次开关能量(125度)(mJ) | 1530 |
Diode单次开关能量(25度)(mJ) | 840 |
换流站损耗计算结果:以送端换流站向受端换流站传输功率为正方向,
工况为:额定电压下双端传输功率均为P=1.0pu,Q=0.2pu。
1)送端正极换流站A相17号子模块:
2)送端正极A相上桥臂:
3)换流站整体损耗分布及损耗率:
更进一步的,研究各种工况下阀损耗的分布情况:
1)双端传输功率P=+500MW,Q=+200MVar
2)双端传输功率P=+250MW,Q=+200MVar
3)双端传输功率P=+25MW,Q=+200MVar
4)双端传输功率P=-250MW,Q=+200MVar
5)双端传输功率P=-500MW,Q=+200MVar
工况的改变会引起通态损耗的大幅变化,也会引起各个器件温度不同程度的波动。以送端正极换流站为例,其分布如图12所示。以送端正极换流站A相17号子模块为例,其在上述工况下的温度变化折线见图13,其中无功功率传输值均为+200MVar。
本实施方法可以得到每个器件的损耗分布和工作温度、单个MMC换流器损耗以及系统整体阀损耗和对应的损耗率,有助于换流站设计阶段的器件选型。得益于快速模型的建立,利用电磁暂态仿真结果的阀损耗计算方法更适用于工程应用,且具有较高的精确性。
显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也意图包含这些改动和变型在内。
本说明书中未作详细描述的内容属于本领域专业技术人员公知的现有技术。
Claims (7)
1.一种模块化多电平换流器阀损快速评估算法,其特征在于:包括以下步骤:
步骤一、根据系统参数和控制方法,在电磁暂态仿真软件中建立MMC的快速仿真模型;
步骤二、该快速仿真模型将仿真结果中的子模块电容电压、开关器件触发脉冲和桥臂电流保存为可读取数据文件;
步骤三、根据MMC子模块中开关器件的特性参数,所述MMC子模块中开关器件即为IGBT及其反并联二极管,建立开关器件损耗计算和PN结温度估计的数学模型;
步骤四、将仿真结果的数据文件导入损耗计算程序,计算出每个数据步长下的单个IGBT和单个二极管的阀损耗;
步骤五、根据IGBT和二极管的阀损耗功率与结温的关系,通过迭代计算,得到单个IGBT和单个二极管修正的阀损耗功率,进而叠加得到MMC的阀损耗功率。
2.根据权利要求1所述的模块化多电平换流器阀损快速评估算法,其特征在于:所述步骤一中,MMC快速仿真模型的子模块电容电压动态简化公式表示为:
其中:△t是积分步长,i为子模块序号,UCi(t)是子模块电容电压,UCi(t-△t)是前一积分步长的电压值,SCi(t)是开关函数,取值0或1分别表示子模块切除或投入,Iarm代表每一相的上桥臂或下桥臂电流,C为电容值。
3.根据权利要求2所述的模块化多电平换流器阀损快速评估算法,其特征在于:所述步骤一中,MMC快速仿真模型的桥臂为等效受控电压源结构,其控制电压值瞬时值Uup(t)为:
其中,USMi(t)为第i个子模块的瞬时端口电压,Uconi(t)为开关管的导通压降,N为单个桥臂上串联子模块的数量,UCi(t)是子模块电容电压,SCi(t)是开关函数,取值0或1分别表示子模块切除或投入。
4.根据权利要求1至3中任一项所述的模块化多电平换流器阀损快速评估算法,其特征在于:所述步骤三中的损耗计算数学模型如下:
PTcon(t)=Uce(t)·IT(t)=(RT(TjT(t))·IT(t)+Uce0(TjT(t)))·IT(t)
PDcon(t)=UD(t)·ID(t)=(RD(TjD(t))×ID(t)+UD0(TjD(t)))·ID(t)
式中,IT(t)、ID(t)分别代表流经IGBT和二极管的瞬时电流,由仿真结果经处理得到;PTcon(t)、PTsw(t)和PToff(t)分别为IGBT导通损耗、开关损耗和关断损耗;PDcon(t)、PDsw(t)和PDoff(t)分别为二极管导通损耗、开关损耗和关断损耗;Uce0(t)、UD(t)分别为t时刻IGBT和二极管的导通压降;TjT(t)、TjD(t)分别为t时刻IGBT和二极管的结温,其初值选取为器件底板温度Tcase;RT是IGBT导通电阻,Uce0是IGBT擎住电压;Roff_T,Roff_D分别为IGBT和二极管的截止电阻,T0为采样间隔,N0为采样间隔时间内器件开关次数;Eon(j)、Eoff(j)和Erec(j)分别为IGBT单次开能量、IGBT单次关能量和二极管单次关断能量,j为采样间隔内开关次数的序号。
5.根据权利要求4所述的模块化多电平换流器阀损快速评估算法,其特征在于:所述步骤三中,所述的IGBT单次开能量Eon(j)、单次关能量Eoff(j)和二极管单次关断能量Erec(j)分别为:
式中,a1、a2、a3、b1、b2、b3、c1、c2、c3是根据器件表开关能量曲线经二次拟合得到的系数,UceN、UDN为IGBT和二极管的额定电压,ion_T(t)、ioff_T(t)为IGBT开关电流瞬时值,ioff_D(t)为二极管关断电流瞬时值,USM(t)为子模块瞬时电压,ρTon、ρToff、ρDrec为结温修正系数。
6.根据权利要求5所述的模块化多电平换流器阀损快速评估算法,其特征在于:所述步骤三中,所述的结温修正系数ρTon、ρToff、ρDrec计算方法为:
式中,Eon1、Eon2分别为IGBT在125°和25°时的单次开通能量;Eoff1、Eoff2分别为IGBT在125°和25°时的单次关断能量;Erec1、Erec2分别为二极管在125°和25°时的单次关断能量,TjT、TjD分别为IGBT和二极管的结温。
7.根据权利要求6所述的模块化多电平换流器阀损快速评估算法,其特征在于:所述步骤五中,IGBT和二极管结温通过下式求得:
TjT(t+Δt)=PT(TjT(t))·(Zth(JC_T)+Zth(CS))+TS
TjD(t+Δt)=PD(TjD(t))·(Zth(JC_D)+Zth(CS))+TS
式中,PT、PD分别为IGBT和二极管总的阀损,Zth(JC_T)、Zth(JC_D)为IGBT和底板、二极管和底板间的温阻;Zth(CS)为底板和散热器间的温阻,TS为散热器温度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510365972.5A CN104993715B (zh) | 2015-06-29 | 2015-06-29 | 模块化多电平换流器阀损快速评估算法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510365972.5A CN104993715B (zh) | 2015-06-29 | 2015-06-29 | 模块化多电平换流器阀损快速评估算法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104993715A true CN104993715A (zh) | 2015-10-21 |
CN104993715B CN104993715B (zh) | 2017-11-10 |
Family
ID=54305482
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510365972.5A Active CN104993715B (zh) | 2015-06-29 | 2015-06-29 | 模块化多电平换流器阀损快速评估算法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104993715B (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2018019376A1 (en) * | 2016-07-28 | 2018-02-01 | Abb Schweiz Ag | Heat balancing in a power converter |
CN107863781A (zh) * | 2017-11-15 | 2018-03-30 | 华北电力大学 | 一种模块化多电平换流阀损耗确定方法及简化方法 |
CN108647447A (zh) * | 2018-05-11 | 2018-10-12 | 中电普瑞电力工程有限公司 | Mmc换流阀可靠性分析方法及装置 |
CN108829982A (zh) * | 2018-06-21 | 2018-11-16 | 福州大学 | 模块化多电平换流器能量等效建模方法 |
CN109001541A (zh) * | 2018-09-10 | 2018-12-14 | 广东电网有限责任公司 | 一种mmc子模块的电容量在线检测方法及装置 |
CN109991872A (zh) * | 2017-12-29 | 2019-07-09 | 上海科梁信息工程股份有限公司 | 一种模块化多电平变换器仿真方法 |
CN110826170A (zh) * | 2019-09-18 | 2020-02-21 | 中国东方电气集团有限公司 | 一种电力电子变流器功率器件损耗实时计算系统 |
CN111464061A (zh) * | 2020-03-26 | 2020-07-28 | 长沙理工大学 | 模块化中压大功率逆变电源单载波调制损耗建模方法 |
CN113489042A (zh) * | 2021-07-01 | 2021-10-08 | 中广核风电有限公司 | 柔性直流输电系统功率子模块损耗的计算方法及装置 |
CN113644834A (zh) * | 2021-10-14 | 2021-11-12 | 广东电网有限责任公司珠海供电局 | 一种模块化多电平换流器损耗评估方法及系统 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2811629A2 (en) * | 2013-06-05 | 2014-12-10 | LSIS Co., Ltd. | Multilevel inverter |
-
2015
- 2015-06-29 CN CN201510365972.5A patent/CN104993715B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2811629A2 (en) * | 2013-06-05 | 2014-12-10 | LSIS Co., Ltd. | Multilevel inverter |
Non-Patent Citations (3)
Title |
---|
周诗嘉 等: ""两电平VSC与MMC通用型平均值仿真模型"", 《电力系统自动化》 * |
张哲任 等: ""基于分段解析公式的MMC-HVDC阀损耗计算方法"", 《电力系统自动化》 * |
薛英林 等: ""采用不同子模块的MMC-HVDC阀损耗通用计算方法"", 《电力自动化设备》 * |
Cited By (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2018019376A1 (en) * | 2016-07-28 | 2018-02-01 | Abb Schweiz Ag | Heat balancing in a power converter |
CN107863781A (zh) * | 2017-11-15 | 2018-03-30 | 华北电力大学 | 一种模块化多电平换流阀损耗确定方法及简化方法 |
CN107863781B (zh) * | 2017-11-15 | 2019-10-01 | 华北电力大学 | 一种模块化多电平换流阀损耗确定方法及简化方法 |
CN109991872B (zh) * | 2017-12-29 | 2022-06-07 | 上海科梁信息科技股份有限公司 | 一种模块化多电平变换器仿真方法 |
CN109991872A (zh) * | 2017-12-29 | 2019-07-09 | 上海科梁信息工程股份有限公司 | 一种模块化多电平变换器仿真方法 |
CN108647447A (zh) * | 2018-05-11 | 2018-10-12 | 中电普瑞电力工程有限公司 | Mmc换流阀可靠性分析方法及装置 |
CN108647447B (zh) * | 2018-05-11 | 2022-07-01 | 中电普瑞电力工程有限公司 | Mmc换流阀可靠性分析方法及装置 |
CN108829982A (zh) * | 2018-06-21 | 2018-11-16 | 福州大学 | 模块化多电平换流器能量等效建模方法 |
CN108829982B (zh) * | 2018-06-21 | 2021-07-13 | 福州大学 | 模块化多电平换流器能量等效建模方法 |
CN109001541A (zh) * | 2018-09-10 | 2018-12-14 | 广东电网有限责任公司 | 一种mmc子模块的电容量在线检测方法及装置 |
CN110826170A (zh) * | 2019-09-18 | 2020-02-21 | 中国东方电气集团有限公司 | 一种电力电子变流器功率器件损耗实时计算系统 |
CN110826170B (zh) * | 2019-09-18 | 2023-07-11 | 中国东方电气集团有限公司 | 一种电力电子变流器功率器件损耗实时计算系统 |
CN111464061A (zh) * | 2020-03-26 | 2020-07-28 | 长沙理工大学 | 模块化中压大功率逆变电源单载波调制损耗建模方法 |
CN111464061B (zh) * | 2020-03-26 | 2023-01-03 | 长沙理工大学 | 模块化中压大功率逆变电源单载波调制损耗建模方法 |
CN113489042A (zh) * | 2021-07-01 | 2021-10-08 | 中广核风电有限公司 | 柔性直流输电系统功率子模块损耗的计算方法及装置 |
CN113489042B (zh) * | 2021-07-01 | 2023-07-18 | 中广核风电有限公司 | 柔性直流输电系统功率子模块损耗的计算方法及装置 |
CN113644834B (zh) * | 2021-10-14 | 2022-02-15 | 广东电网有限责任公司珠海供电局 | 一种模块化多电平换流器损耗评估方法及系统 |
CN113644834A (zh) * | 2021-10-14 | 2021-11-12 | 广东电网有限责任公司珠海供电局 | 一种模块化多电平换流器损耗评估方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN104993715B (zh) | 2017-11-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104993715A (zh) | 模块化多电平换流器阀损快速评估算法 | |
Shen et al. | Real-time device-level transient electrothermal model for modular multilevel converter on FPGA | |
Liu et al. | Modeling and analysis of HVDC converter by three-phase dynamic phasor | |
CN103324843A (zh) | 一种适用于不同子模块类型的mmc阀损耗计算方法 | |
CN102969888B (zh) | 基于rtds的mmc自定义子模块的设计方法 | |
CN106649927B (zh) | 一种基于fpga的电力电子元件实时仿真组合建模方法 | |
CN102403916B (zh) | 一种仿真提速电路的设计方法 | |
CN103914599B (zh) | 一种模块化多电平换流器的戴维南等效整体建模方法 | |
CN112464443B (zh) | 一种用于电力电子变流器igbt结温波动的计算方法 | |
CN102654565A (zh) | 一种电池模拟方法及模拟器 | |
CN102708225B (zh) | 一种交直流大电网电磁暂态仿真的分片调试方法 | |
CN101572408A (zh) | 一种提高高压直流输电系统仿真能力的方法 | |
CN103018583B (zh) | 基于mmc柔性直流输电系统电平数选择验证方法 | |
CN103605850A (zh) | 一种带子模块闭锁功能的mmc等效建模方法 | |
CN108536949A (zh) | 基于电感/电容开关模型的lcc-hvdc仿真模型和参数优化方法 | |
CN108631632A (zh) | 一种基于虚拟桥臂数学模型的mmc瞬时功率损耗精细计算方法 | |
CN104199997A (zh) | 一种mmc多子模块自定义集成元件的设计方法 | |
CN109241678A (zh) | 模块化多电平换流器多拓扑通用仿真模型 | |
CN103605828B (zh) | 含变流器电力元件快速仿真建模方法 | |
Gao et al. | Fast simulation model of voltage source converters with arbitrary topology using switch-state prediction | |
CN102593876A (zh) | 一种大规模交直流互联电力系统的连续潮流算法 | |
CN103580054B (zh) | 一种换流阀交流系统1s故障耐受特性分析方法 | |
CN106941260B (zh) | 一种模块化多电平换流器的性能评价方法及装置 | |
CN110427635A (zh) | Lcc-hvdc优化电磁暂态仿真方法 | |
Arifujjaman et al. | Dynamic simulation and control of a battery energy storage system |
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 |