Loading [MathJax]/jax/output/SVG/jax.js
  • 匿名盲审
  • 学术期刊非营利性
  • 全球免费开放获取全文
  • 最新科研成果提供绿色通道

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

MMC子模块电容电压波动及谐波研究

邹凯凯 李钢 邹强 董云龙 卢宇

邹凯凯, 李钢, 邹强, 董云龙, 卢宇. MMC子模块电容电压波动及谐波研究[J]. 南方能源建设, 2023, 10(5): 129-138. doi: 10.16516/j.gedi.issn2095-8676.2023.05.016
引用本文: 邹凯凯, 李钢, 邹强, 董云龙, 卢宇. MMC子模块电容电压波动及谐波研究[J]. 南方能源建设, 2023, 10(5): 129-138. doi: 10.16516/j.gedi.issn2095-8676.2023.05.016
ZOU Kaikai, LI Gang, ZOU Qiang, DONG Yunlong, LU Yu. Capacitor Voltage Ripple and Harmonics Analysis in MMC[J]. SOUTHERN ENERGY CONSTRUCTION, 2023, 10(5): 129-138. doi: 10.16516/j.gedi.issn2095-8676.2023.05.016
Citation: ZOU Kaikai, LI Gang, ZOU Qiang, DONG Yunlong, LU Yu. Capacitor Voltage Ripple and Harmonics Analysis in MMC[J]. SOUTHERN ENERGY CONSTRUCTION, 2023, 10(5): 129-138. doi: 10.16516/j.gedi.issn2095-8676.2023.05.016

MMC子模块电容电压波动及谐波研究

DOI: 10.16516/j.gedi.issn2095-8676.2023.05.016
CSTR: 32391.14.j.gedi.issn2095-8676.2023.05.016
基金项目: 国家电网有限公司总部管理科技项目资助“新能源经柔性直流送出系统在多类型谐波扰动下的穿越能力及控制策略研究”(5100-202256371A-2-0-ZN)
详细信息
    作者简介:

    邹凯凯,1991-,男,工程师,华北电力大学电气工程及其自动化硕士,主要研究方向为直流输电仿真建模、控制保护策略开发(e-mail)zoukk@nrec.com

    李钢,1983-,男,高级工程师,电力系统及其自动化,主要从事柔性直流控制保护策略开发相关工作(e-mail)lig@nrec.com

    邹强,1981-,男,高级工程师,电力系统及其自动化,主要从事柔性直流控制保护策略开发相关工作(e-mail)zouq@nrec.com

    董云龙,1977-,男,教授级高级工程师,电力系统及其自动化,主要从事电力电子设备及控制保护研究相关工作(e-mail)dongyl@nrec.com

    卢宇,1979-,男,教授级高级工程师,电力系统及其自动化,主要从事电力电子设备及控制保护研究相关工作(e-mail)luy@nrec.com

    通讯作者:

    邹凯凯,1991-,男,工程师,华北电力大学电气工程及其自动化硕士,主要研究方向为直流输电仿真建模、控制保护策略开发(e-mail)zoukk@nrec.com

  • 中图分类号: TM7;TM712

Capacitor Voltage Ripple and Harmonics Analysis in MMC

  • 摘要:   目的  柔性直流的快速发展为MMC(Modular Multi-Level Converter,模块化多电平变流器)拓扑确定了实践支撑。文章以半桥MMC换流器为例,研究MMC稳态电容电压波动、谐波交互规律及桥臂开关函数的关系,以对MMC系统设计、控制策略给予理论支持。  方法  通过建立子模块电容电压、桥臂开关函数的解析表达式,得到MMC换流器上下桥臂能量波动及电容电压波动的数学近似表达式,进而研究了上下桥臂电流、交流阀侧、直流侧中各次谐波的分布规律及相间环流的正负序关系,基于谐波公式得到二次环流的估算方法;进一步,采用变量对照法研究了桥臂开关函数、子模块电容电压随阀侧功率因数、子模块电容大小的变化关系;最后基于实际工程参数搭建了详细的PSCAD/EMTDC模型,开展了理论解析与离线仿真的一致性对比。  结果  研究表明:MMC相间环流只含偶数次谐波,且其中6k+2次环流呈负序、6k+4次环流呈正序、6k次环流呈零序特性,同时稳态下MMC阀侧出口交流电压电流只含奇数次谐波,MMC直流侧电压电流只含偶数次谐波,且在功率因数角φ<0的某些范围内桥臂开关函数峰值会大于1,在低功率因数下子模块电容电压波动峰值大于高功率因数下的波动。  结论  可由此开展换流器电容容值的选型设计及暂稳态策略研究,并通过理论解析掌握MMC拓扑的运行本征特性。
  • 随着电力电子器件的发展,直流输电进入到了第三代新型直流输电,其基于全控型器件IGBT(Insulated-Gate Bipolar Transistor,绝缘栅双极晶体管)。采用全控型器件可以使换流器在四象限运行,实现有功无功解耦控制,因此特别适用于海上风电联网、无源网络供电、异步交流电网互联等。

    传统的两电平或三电平电压源换流器(Voltage Source Converter,VSC)由于谐波含量大、容量低、电压等级低、器件串联动态均压并联均流等问题,发展受到了制约。2001年德国学者提出了模块化多电平换流器MMC(Modular Multilevel Converter,模块化多电平变流器)概念[1-3],西门子公司于2010年将其用于实际工程。近年来,基于MMC的工程越来越多,据不完全统计,目前已投运30项工程。国外工程集中在欧洲,主要将其用于海上风电的送出和电网互联。近年来国内已有南汇、厦门、舟山、鲁西、渝鄂、张北柔性直流电网、乌东德混合直流、三峡如东海风柔性直流、白鹤滩混合直流等工程。

    文献[4-6]给出了MMC的基本数学关系和静态特性;文献[7]研究了采用PWM( Pulse Width Modulation,脉宽调制)方式的谐波特性;文献[8-10]对MMC的动态特性进行了建模研究;文献[11-13]研究了MMC联接弱交流电网的稳定运行区间;文献[14-16]开展了MMC孤岛送新能源场景下的控制研究;文献[17-18]就MMC处于直流电压控制下的稳定性开展了研究分析;文献[19-21]就MMC从电容电压控制出发开展的控制手段进行了研究。

    近年来,对于MMC的研究主要集中于MMC与传统LCC(Line-Commuted Converter,换相换流器)直流的混合特性、主动故障隔离及故障清除策略、接入弱电网时的稳定性分析、新能源经柔性直流送出的稳定性分析、自适应强弱电网等方向。

    其中对于柔性直流稳定性的分析成为重中之重,当前的分析手段是针对柔性直流系统的稳态工作点建立小信号模型,基于小信号模型得到柔性直流系统的阻抗特性,再根据柔性直流及外部电网的阻抗特性进行稳定性分析[22-24]。因此,对MMC稳态建模是稳定性分析的基石。

    此外,高压大容量MMC换流器设备造价主要取决于IGBT器件及并联电容器的造价。当系统容量一定时,IGBT器件选型较为固定,并联电容器的大小选型的重要参考是子模块电容电压的波动。因此,对MMC的稳态工作点的精细化建模特别是模块电容电压分析具有重要意义。已发表的相关文献侧重点在于MMC内部环流产生原因及环流抑制的相关策略,文章则侧重电容能量变化规律,从子模块电容能量变化角度出发给出一种精细化数学建模,在推导出桥臂电容能量的波动公式前提下,深入研究了桥臂电容电压随时间、功率因数角、子模块电容大小的波动规律,根据电容电压的波动得到MMC谐波分布特性,并得出二次谐波的估算公式;在电容电压波动前提下,研究了桥臂子模块投入个数和模块电容电压随时间、功率因数的关系;最后对得到的规律进行了详细的仿真验证。

    典型的半桥MMC换流器拓扑如图1所示,各相单元含上下两个桥臂,各桥臂单元由N个MMC子模块及桥臂电抗L组成,图中R为子模块等效内阻。子模块采用半桥结构,内含两个具有反并联二极管的IGBT和电容组成,实际工程中每个子模块还含一个保护晶闸管和快速开关。在直流侧通过钳位电容或钳位大电阻接地提供接地钳位点时[25],此时正负极母线上电压分别为Vd/2及−Vd/2。

    图  1  MMC换流器拓扑图
    Fig.  1  Topological structure of MMC

    对相单元与交直流侧形成回路的利用基尔霍夫电压定律( Kirchhoff's Voltage Law,KVL)有:

    {Vd2vuRiuLdiudt=vaVd2+vl+Ril+Ldildt=va 1

    式中:

    L ——桥臂电抗值(H);

    R ——子模块的等效内阻(Ω);

    Vd ——直流电压(kV);

    vuvl ——上、下桥臂端间电压(kV);

    iuil ——上、下桥臂电流(kA),电流方向如图1所示;

    va ——阀侧A相相电压(kV)。

    对式(1)进行数学变换得:

    {L2disdt=vu+vl2vsvaR2isLdicdt=vu+vl2vc+Vd2Ric 2

    式中:

    is ——阀侧相电流(kA);

    ic ——流过桥臂的直流电流分量(kA);

    vs ——上下桥臂的差模电压(kV);

    vc ——上下桥臂的共模电压(kV)。

    式(2)中用到了如下关系[26]

    {iu=is2+icil=is2+ic 3

    根据每个子模块投入切除状态定义单个子模块jj=1,2,···,N)投入切除的开关状态数nju,l(为1表示投入,为0表示切除),在单桥臂中所有子模块电压均衡时,可得到单个桥臂总投入子模块等效电压为:

    vcu,l=Nj=1nju,lvjcu,l=vΣcu,lNNj=1nju,l=nu,lvΣcu,l 4

    式中:

    vΣcu,l ——上/下桥臂所有子模块电压之和(kV);

    N ——上/下桥臂子模块个数;

    nu,l ——上/下桥臂的总开关函数,nu为上桥臂的总开关函数,nl为下桥臂的总开关函数。

    考虑每个子模块电容电压波动后将其叠加得到桥臂上所有MMC子模块两端电压波动,以上桥臂为例,有:

    Cdvjcudt=njuiu 5
    CNj=1dvjcudtdvΣcu/dt=Nj=1njuiu=iuNj=1njuNnu 6

    式中:

    C ——MMC子模块电容容值(F)。

    替换下标u为表示下桥臂的l,公式依然成立。将式(6)重写得到上下桥臂的总子模块电容电压波动方程如下:

    {CNdvΣcudt=nu(is2+ic)CNdvΣcldt=nl(is2+ic) 7

    将式(4)代入到vsvc的公式中,可得到桥臂的总开关函数为:

    nu=vcvsvΣcunl=vc+vsvΣcl 8

    结合式(7)和式(8),可以得到上下桥臂电容能量的波动方程:

    {vΣcuCNdvΣcudt=C2Nd(vΣcu)2dtdWu/dt=(vcvs)(is2+ic)vΣclCNdvΣcldt=C2Nd(vΣcl)2dtdWl/dt=(vc+vs)(is2+ic) 9

    从式(9)进一步得到某相单元中上下桥臂总电容能量的波动及上下桥臂电容能量之差的波动方程为:

    {WΣ=Wu+WlWΔ=WuWl{dWΣdt=2vcicvsisdWΔdt=vcis2vsic 10

    表达式中vc=Vd/2;vs=Vscos(w1t);is=Iscos(w1tφ)。

    将上述关系代入式(10)中,可得:

    {dWΣdt=VdicVsIs2cosφVsIs2cos(2w1tφ)dWΔdt=VdIs2cos(w1tφ)2Vsiccos(w1t) 11

    从式(11)的第一个公式可以看出,相单元的总能量变化主要由三部分构成:第一部分为直流分量对应到单相中的有功功率;第二部分为交流系统注入MMC换流器的有功功率;第三部分为表示相间交换的有功,其以二倍频为主,表示相间环流以二倍频为主。

    将式(11)两边积分得:

    {WΣ=W0VsIs4w1sin(2w1tφ)ΔWΣWΔ=WΔ0+VdIs2w1sin(w1tφ)2Vsicw1sin(w1t)ΔWΔ 12

    式中:

    ΔWΣ ——相单元总能量的波动量(MJ);

    ΔWΔ ——上下桥臂能量差的波动量(MJ)。

    由于能量储存于子模块电容中,假设不考虑实际工程中的子模块冗余个数,则有:

    {WΣ0=2N×12C(VdN)2=CV2dNWΔ0=0 13

    将式(13)代入式(12)可得到上下桥臂总子模块电容能量为:

    {Wu=WΣ+WΔ2=WΣ0+ΔWΣ+ΔWΔ2Wl=WΣWΔ2=WΣ0+ΔWΣΔWΔ2 14

    根据桥臂能量与电容电压之间的关系,结合式(12)~式(14)可最终得到上下桥臂MMC子模块电容电压的详细波动关系如下:

    {vΣcu=2NCWu=Vd1+NCV2d(ΔWΣ+ΔWΔ)Vd+N2CVd(ΔWΣ+ΔWΔ)ΔvΣcuvΣcl=2NCWl=Vd1+NCV2d(ΔWΣΔWΔ)Vd+N2CVd(ΔWΣΔWΔ)ΔvΣcl 15

    式中:

    ΔvΣcuΔvΣcl——上、下桥臂子模块电容电压之和的波动量(kV)。

    为定性初步分析基波分量与二倍频分量的规律,在泰勒展开时忽略了更高频分量,从式(15)可初步得到如下结论:

    1)桥臂电容电压和的波动与电容C成反比。

    2)上下桥臂的电容电压和的波动均含两部分,一部分为ΔWΣ二倍频分量,一部分为ΔWΔ基频分量。

    3)二倍频分量在上下桥臂总大小相等,方向相同;基波分量在上下桥臂上大小相等,方向相反。

    将式(15)泰勒展开后取其任一项有:

    {(ΔWΣ+ΔWΔ)n=nk=0CknΔWΣkΔWΔnk(ΔWΣΔWΔ)n=nk=0CknΔWΣk(ΔWΔ)nk 16

    对典型项CknΔWΣkΔWΔnkCknΔWΣk(ΔWΔ)nk进行分析,由于ΔWΣ为二倍频分量,​ΔWΔ为基波分量,可以看出ΔWΣk项仍为偶数倍频,谐波电压的奇偶关系与ΔWΔnk(ΔWΔ)nkn−k的奇偶性相同。因此,可以得到上下桥臂的电容电压波动量的奇数次谐波呈谐波大小相等,方向相反的关系;偶数次谐波大小相等,方向相同的关系。由于电容电压的波动会带来桥臂上的谐波电流,因此,我们可以得出桥臂电流上谐波电流分布规律:

    1)上下桥臂的奇数次谐波大小相等,方向相反。

    2)上下桥臂的偶数次谐波大小相等,方向相同。

    3)相间环流为偶数次环流,且以二倍频为主。

    4)第1)条等同于MMC阀侧出口交流电压电流只含奇数次谐波。

    5)第2)条等同于MMC直流侧电压电流只含偶数次谐波。

    6)根据ΔWΣ和​ΔWΔ的计算公式,可以得到相间环流只含偶数次谐波,其中6k+2​次环流呈负序;6k+4​次环流呈正序;6k​次环流呈零序[12]

    从式(17)、式(18)的电容电压波动泰勒展开式出发,考虑取2w1t的相关部分后可以达到二次谐波电流的估算大小。由于高频分量在进行三角函数计算时也有可能产生二倍频分量,这里忽略这些高频分量因中间计算过程带来的二倍频分量的大小,忽略的原因在于随着展开的幂次数的增加,其分量的幅值越来越小,图2给出了泰勒展开式的系数函数随展开阶数的关系。

    图  2  泰勒展开式系数函数
    Fig.  2  Coefficient function of Talyor Series
    vΣcu=vd1+NCv2d(ΔWΣ+ΔWΔ)=vd+N2Cvd(ΔWΣ+ΔWΔ)++12×(121)×(122)×(12k+1)vdk!×[NCv2d(ΔWΣ + ΔWΔ)]k+ 17
    vΣcl=vd1+NCv2d(ΔWΣΔWΔ)=vd+N2Cvd(ΔWΣΔWΔ)++12×(121)×(122)××(12k+1)k!×vd[NCv2d(ΔWΣΔWΔ)]k+ 18

    得到二次谐波电流的估算大小为:

    i2harm=CNdΔvcudt=Asin(2w1t)+Bcos(2w1t)=A2+B2cos(2w1t+φ) 19
    A=VsIs4VdsinφNV2si2c2w1CV3dNV2dI2s32w1CV3dcos(2φ)+NVdicVsIs8w1CV3dcosφ 20
    B=VsIs4VdcosφNV2dI2s32w1CV3dsin(2φ)+NVdicVsIs8w1CV3dsinφ 21

    式(22)、式(23)给出了上、下桥臂子模块电容电压的波动规律,从式子可以知道在忽略高频谐波的情况下桥臂电容电压和的波动与电容C成反比,在直流电压、直流电流、阀侧相电压相电流幅值一定的情况下,子模块电容电压的波动主要取决于阀侧功率因数角φ和时间t。为深入研究,本章节以我国首个海风柔性直流江苏如东工程参数为算列开展,其主参数如表1所示。

    表  1  MMC参数
    Tab.  1  Parameters for MMC
    参数名
    额定直流电压/kV ±400
    额定直流电流/kA 1.375
    额定有功/GW 1.100
    额定视在功率/MVA 1230
    子模块个数不考虑冗余/个 400
    子模块电容容值/mF 4, 6, 8, 9, 10, 12
    桥臂电抗器/H 0.133
    阀侧相电压峰值/kV 340
    阀侧相电流峰值/kA 2.411
    下载: 导出CSV 
    | 显示表格
    vΣcu=vd1+NCv2d(ΔWΣ+ΔWΔ)=vd1+NCv2d(VsIs4w1sin(2w1tφ)+vdIs2w1sin(w1tφ)2Vsicw1sin(w1t)) 22
    vΣcl=vd1+NCv2d(ΔWΣΔWΔ)=vd1+NCv2d(VsIs4w1sin(2w1tφ)vdIs2w1sin(w1tφ)+2Vsicw1sin(w1t)) 23

    根据式(22)、式(23),假设阀侧电压与阀侧电流的相位差在[−π, π]之间变动,考虑4个工频周期下的波动,得到如下三维关系图,其中x轴坐标为时间t(单位为s),y轴坐标为相角差φ(单位为rad),z轴坐标为桥臂子模块电容电压的波动量ΔvΣcu(单位为kV),根据前面式子,得到上下桥臂所有投入的子模块电容电压的波动随时间及阀侧电压初始相位差的三维关系图,如图3图4所示。

    图  3  上桥臂子模块电容电压的波动
    Fig.  3  Capacitor voltage ripple of upper bridge
    图  4  下桥臂子模块电容电压的波动
    Fig.  4  Capacitor voltage ripple of lower bridge

    从上图可知,在φ=0时,MMC换流阀与交流侧没有无功交换,此时桥臂子模块电容电压波动的峰值最小。在φ减小或增大的过程中,波动量均在一定程度上增大。

    根据式(8)同样可得到上下桥臂的总开关函数nunl随时间及功率因数角φ的曲线图。其中x轴坐标为时间t(单位为s),y轴坐标为相角差φ(单位为rad),z轴坐标为桥臂总子模块开关函数,图5为上桥臂总开关函数图,图6为下桥臂总开关函数图。

    图  5  上桥臂总开关函数
    Fig.  5  Switching function of upper bridge
    图  6  下桥臂总开关函数
    Fig.  6  Switching function of lower bridge

    图5图6可以明显看出,功率因数角φ一定时桥臂总开关函数是时间t的三角工频周期函数,在时间t一定时,开关函数是功率因数角的周期函数,且在φ<0的一定范围内,开关函数的最大值大于1,这是由于φ<0时,换流器向交流系统注入容性无功,在[−90°, 0]范围内,φ越小注入容性无功越大,此时同等有功条件下需要投入的子模块个数越多。将式(15)代入式(8)以后,可得到上下桥臂的开关函数的表达式,这里仅列出下桥臂开关函数:

    nl=Vd/2+Vscosw1tVd+ΔvΣcl 24
    ΔvΣcl=N2Cvd[VsIs4w1sin(2w1tφ)vdIs2w1sin(w1tφ)+2Vsicw1sin(w1t)] 25

    t=0时刻,有桥臂最大的开关函数为:

    nl|max=Vd/2+VsVd+NIs4w1C(Vs2Vd+1)sinφ 26

    在开关函数最大值为1的情况下,此时得出:

    sinφ=4w1CNIs×VsVd2Vs2Vd+1 27

    由于阀侧相电压最大值小于直流电压的一半,可以看出在开关函数的最大值为1时,功率因数角φ∈[−π/2, 0],这与图中结论一致。

    为进一步分析开关函数与模块电容电压随功率因数角的变换关系,给出了不同功率因数角下的开关函数及桥臂子模块电容电压波动的关系图(图7)。从图中及式(24)可得到如下结论:

    图  7  不同功率因数角下桥臂开关函数及模块电压波动
    Fig.  7  Switching function and capacitor voltage ripple for different phi

    1)桥臂开关函数的最大值在φ∈[−π/2, π/2]范围内单调递减,且在[−π/2, 0]内某点φ0处其最大值为1。

    2)在额定视在功率一定的情况下,桥臂开关函数的最大值在φ<0的某些范围内会呈现大于1的情况,此时功率因数较小,理论交换无功较大,φ=−90°时换流器工作在STATCOM工况下,实际工程中采用限幅环节对开关函数进行限幅。

    3)子模块电容电压纹波存在明显的二倍频和基波分量。

    4)在低功率因数φ=±90°下,子模块电容电压纹波比高功率因数下的纹波更大,且在低功率因数下子模块电容电压波动呈现正负半周期不对称现象。

    考虑φ=0时,子模块电容电压及桥臂开关函数随电容容值的变化关系。图8给出了表1参数下的桥臂开关函数及电容电压波动随时间和电容容值的变化。

    图  8  不同子模块电容容值下桥臂开关函数及模块电压波动
    Fig.  8  Switching function and capacitor voltage ripple for different C

    图8中可知,在φ=0时,桥臂开关函数的最大值不随电容电压的变化而变化,进一步从式(20)中可知,在φ≠0时,开关函数的最大值一定程度上随电容容值的变大而变小。子模块电容电压的波动量随着电容容值的变大而减少。

    为验证上述结论,在PSCAD/EMTDC仿真软件中搭建了上述表1参数下的仿真模型。

    为验证子模块电容电压波动随子模块电容容值的变化,图9给出了计算与仿真的对比关系图,其横坐标为子模块电容容值,纵坐标为电容电压的波动幅值相对基准值的占比,可看出子模块电容电压波动幅值随容值的变大而减小,理论计算值与仿真值较为吻合。

    图  9  电容电压波动幅值随电容的变化关系
    Fig.  9  Capacitor voltage ripple for different capacitance

    为验证子模块电容电压波动实时曲线及桥臂开关函数的实时曲线,在仿真模型中任取一个PQ功率点(1100 MW, 123 Mvar)下对两者进行了对比。从图10对比中可看出在任取的功率点下,上下桥臂开关函数理论计算值与仿真值吻合度较高。图11给出了上下桥臂子模块电容电压实时波动的对比,可看出计算与仿真总体较为吻合,在模块电容电压波动过零点处的不一致主要是由于计算中未考虑二次以上的高倍频谐波波动导致。

    图  10  实时桥臂开关函数计算与仿真对比
    Fig.  10  Switching function comparison between calculation and simulation
    图  11  模块总电容电压实时波动对比
    Fig.  11  Capacitor voltage ripple comparison between calculation and simulation

    文章从能量角度出发得到了桥臂电容电压、桥臂开关的理论计算值,在此基础上研究了桥臂谐波的特性,并深入研究了子模块电容电压的波动和桥臂开关函数的情况,主要得出的结论与仿真结果一致,为实际工程中MMC换流器设计提供了参考。

    但是,文章对换流器建模尚未结合换流器自身控制系统,后续可开展结合控制系统的换流器阻抗建模研究及换流器自身调制环节、控制链路等MMC稳定性影响因子分析;另一方面,通过注入倍频环流或先进控制方式可降低模块电容电压的波动,以优化系统参数的选型,这也是后续待研究之一。

  • 图  1  MMC换流器拓扑图

    Fig.  1  Topological structure of MMC

    图  2  泰勒展开式系数函数

    Fig.  2  Coefficient function of Talyor Series

    图  3  上桥臂子模块电容电压的波动

    Fig.  3  Capacitor voltage ripple of upper bridge

    图  4  下桥臂子模块电容电压的波动

    Fig.  4  Capacitor voltage ripple of lower bridge

    图  5  上桥臂总开关函数

    Fig.  5  Switching function of upper bridge

    图  6  下桥臂总开关函数

    Fig.  6  Switching function of lower bridge

    图  7  不同功率因数角下桥臂开关函数及模块电压波动

    Fig.  7  Switching function and capacitor voltage ripple for different phi

    图  8  不同子模块电容容值下桥臂开关函数及模块电压波动

    Fig.  8  Switching function and capacitor voltage ripple for different C

    图  9  电容电压波动幅值随电容的变化关系

    Fig.  9  Capacitor voltage ripple for different capacitance

    图  10  实时桥臂开关函数计算与仿真对比

    Fig.  10  Switching function comparison between calculation and simulation

    图  11  模块总电容电压实时波动对比

    Fig.  11  Capacitor voltage ripple comparison between calculation and simulation

    表  1  MMC参数

    Tab.  1.   Parameters for MMC

    参数名
    额定直流电压/kV ±400
    额定直流电流/kA 1.375
    额定有功/GW 1.100
    额定视在功率/MVA 1230
    子模块个数不考虑冗余/个 400
    子模块电容容值/mF 4, 6, 8, 9, 10, 12
    桥臂电抗器/H 0.133
    阀侧相电压峰值/kV 340
    阀侧相电流峰值/kA 2.411
    下载: 导出CSV
  • [1] RODRIGUEZ J, LAI J S, PENG F Z. Multilevel inverters: a survey of topologies, controls, and applications [J]. IEEE transactions on industrial electronics, 2002, 49(4): 724-738. DOI:  10.1109/TIE.2002.801052.
    [2] LESNICAR A, MARQUARDT R. An innovative modular multilevel converter topology suitable for a wide power range [C]//2003 IEEE Bologna Power Tech Conference Proceedings, Bologna, Italy, June 23-26, 2003. New York: IEEE, 2003: 6. DOI: 10.1109/PTC.2003.1304403.
    [3] LESNICAR A, MARQUARDT R. A new modular voltage source inverter topology [C]//Proceedings of 10th European Conference on Power Electronics and Applications (EPE), Toulouse, France, September 2-4, 2003. Oxfordshire, UK: Taylor & Francis, 2015.
    [4] NORRGA S, ÄNGQUIST L, ILVES K, et al. Decoupled steady-state model of the modular multilevel converter with half-bridge cells [C]//6th IET International Conference on Power Electronics, Machines and Drives, Bristol, UK, March 27-29, 2012. Hertfordshire, UK: IET, 2012: 1-6. DOI: 10.1049/cp.2012.0232.
    [5] 王姗姗, 周孝信, 汤广福, 等. 模块化多电平电压源换流器的数学模型 [J]. 中国电机工程学报, 2011, 31(24): 1-8. DOI:  10.13334/j.0258-8013.pcsee.2011.24.005.

    WANG S S, ZHOU X X, TANG G F, et al. Modeling of modular multi-level voltage source converter [J]. Proceedings of the CSEE, 2011, 31(24): 1-8. DOI:  10.13334/j.0258-8013.pcsee.2011.24.005.
    [6] SONG Q, LIU W H, LI X Q, et al. A steady-state analysis method for a modular multilevel converter [J]. IEEE transactions on power electronics, 2013, 28(8): 3702-3713. DOI:  10.1109/TPEL.2012.2227818.
    [7] ILVES K, ANTONOPOULOS A, NORRGA S, et al. Steady-state analysis of interaction between harmonic components of arm and line quantities of modular multilevel converters [J]. IEEE transactions on power electronics, 2012, 27(1): 57-68. DOI:  10.1109/TPEL.2011.2159809.
    [8] SAEEDIFARD M, IRAVANI R. Dynamic performance of a modular multilevel back-to-back HVDC system [J]. IEEE transactions on power delivery, 2010, 25(4): 2903-2912. DOI:  10.1109/TPWRD.2010.2050787.
    [9] HARNEFORS L, ANTONOPOULOS A, NORRGA S, et al. Dynamic analysis of modular multilevel converters [J]. IEEE transactions on industrial electronics, 2013, 60(7): 2526-2537. DOI:  10.1109/TIE.2012.2194974.
    [10] HARNEFORS L, ANTONOPOULOS A, ILVES K, et al. Global asymptotic stability of current-controlled modular multilevel converters [J]. IEEE transactions on power electronics, 2015, 30(1): 249-258. DOI:  10.1109/TPEL.2014.2298560.
    [11] 刘昇, 徐政. 联于弱交流系统的VSC-HVDC稳定运行区域研究 [J]. 中国电机工程学报, 2016, 36(1): 133-144. DOI:  10.13334/j.0258-8013.pcsee.2016.01.014.

    LIU S, XU Z. Study on stable operating region of VSC-HVDC connected to weak AC systems [J]. Proceedings of the CSEE, 2016, 36(1): 133-144. DOI:  10.13334/j.0258-8013.pcsee.2016.01.014.
    [12] KRIGE E. Harmonic interaction between weak AC systems and VSC-based HVDC schemes [D]. Stellenbosch: Stellenbosch University, 2012.
    [13] ZHANG L D, HARNEFORS L, NEE H P. Modeling and control of VSC-HVDC links connected to island systems [J]. IEEE transactions on power systems, 2011, 26(2): 783-793. DOI:  10.1109/TPWRS.2010.2070085.
    [14] GIVAKI K, XU L. Stability analysis of large wind farms connected to weak AC networks incorporating PLL dynamics [C]//International Conference on Renewable Power Generation, Beijing, China, October 17-18, 2015. Hertfordshire, UK: IET, 2016: 1-6. DOI: 10.1049/cp.2015.0440.
    [15] ZHANG L D, HARNEFORS L, NEE H P. Interconnection of two very weak AC systems by VSC-HVDC links using power-synchronization control [J]. IEEE transactions on power systems, 2011, 26(1): 344-355. DOI:  10.1109/TPWRS.2010.2047875.
    [16] ADIB A, MIRAFZAL B, WANG X F, et al. On stability of voltage source inverters in weak grids [J]. IEEE access, 2018, 6: 4427-4439. DOI:  10.1109/ACCESS.2017.2788818.
    [17] HUANG Y H, YUAN X M, HU J B, et al. Modeling of VSC connected to weak grid for stability analysis of DC-link voltage control [J]. IEEE journal of emerging and selected topics in power electronics, 2015, 3(4): 1193-1204. DOI:  10.1109/JESTPE.2015.2423494.
    [18] HUANG Y H, YUAN X M, HU J B, et al. DC-bus voltage control stability affected by AC-bus voltage control in VSCs connected to weak AC grids [J]. IEEE journal of emerging and selected topics in power electronics, 2016, 4(2): 445-458. DOI:  10.1109/JESTPE.2015.2480859.
    [19] 黄守道, 彭也伦, 廖武. 模块化多电平型变流器电容电压波动及其抑制策略研究 [J]. 电工技术学报, 2015, 30(7): 62-71. DOI:  10.3969/j.issn.1000-6753.2015.07.008.
    [20] HUANG S D, PENG Y L, LIAO W. Study of capacitor voltage fluctuation and its suppression for modular multilevel converter [J]. Transactions of China electrotechnical society, 2015, 30(7): 62-71. DOI:  10.3969/j.issn.1000-6753.2015.07.008.
    [21] ANGQUIST L, ANTONOPOULOS A, SIEMASZKO D, et al. Open-loop control of modular multilevel converters using estimation of stored energy [J]. IEEE transactions on industry applications, 2011, 47(6): 2516-2524. DOI:  10.1109/TIA.2011.2168593.
    [22] 周月宾, 江道灼, 郭捷, 等. 模块化多电平换流器子模块电容电压波动与内部环流分析 [J]. 中国电机工程学报, 2012, 32(24): 8-14. DOI:  10.13334/j.0258-8013.pcsee.2012.24.002.

    ZHOU Y B, JIANG D Z, GUO J, et al. Analysis of sub-module capacitor voltage ripples and circulating currents in modular multilevel converters [J]. Proceedings of the CSEE, 2012, 32(24): 8-14. DOI:  10.13334/j.0258-8013.pcsee.2012.24.002.
    [23] YANG L G, XU Z H, FENG L, et al. Analysis on harmonic resonance of offshore wind farm transmitted by MMC-HVDC system [C]//2019 IEEE Innovative Smart Grid Technologies - Asia, Chengdu, China, May 21-24, 2019. New York, USA: IEEE, 2019: 2296-2301. DOI: 10.1109/ISGT-Asia.2019.8881373.
    [24] ZOU C Y, RAO H, XU S K, et al. Analysis of resonance between a VSC-HVDC converter and the AC grid [J]. IEEE transactions on power electronics, 2018, 33(12): 10157-10168. DOI:  10.1109/TPEL.2018.2809705.
    [25] ILVES K, ANTONOPOULOS A, HARNEFORS L, et al. Capacitor voltage ripple shaping in modular multilevel converters allowing for operating region extension [C]//IECON 2011 - 37th Annual Conference of the IEEE Industrial Electronics Society, Melbourne, Australia, November 7-10, 2011. New York, USA: IEEE, 2011: 4403-4408. DOI: 10.1109/IECON.2011.6120033.
    [26] WANG H T, TANG G F, HE Z Y, et al. Efficient grounding for modular multilevel HVDC converters (MMC) on the AC side [J]. IEEE transactions on power delivery, 2014, 29(3): 1262-1272. DOI:  10.1109/TPWRD.2014.2311796.
  • [1] 陈娅莉, 卢波, 白兴宇, 马悦心, 许凌飞, 梁军, 左览.  TH2103A型速调管阳极电流测量和保护系统 . 南方能源建设, 2024, 11(3): 146-151. doi: 10.16516/j.ceec.2024.3.16
    [2] 王腾.  托卡马克超导磁体失超电压探测技术 . 南方能源建设, 2024, 11(3): 65-74. doi: 10.16516/j.ceec.2024.3.07
    [3] 林聪, 李湘峰, 郭芳.  直流微电网的故障电流控制器研究 . 南方能源建设, 2023, 10(5): 50-56. doi: 10.16516/j.gedi.issn2095-8676.2023.05.007
    [4] 劳志烜, 郑炳耀, 郭芳, 梅洪灯, 朱文峰, 王瑞阳, 方俊杰.  基于故障电流控制的直流微网过电流保护研究 . 南方能源建设, 2023, 10(5): 65-71. doi: 10.16516/j.gedi.issn2095-8676.2023.05.009
    [5] 李学斌, 刘剑, 马麟, 韩云阳, 赵洁琼, 李富春.  储能参与配电网削峰填谷对网络损耗及电压偏差的影响分析 . 南方能源建设, 2023, 10(6): 51-63. doi: 10.16516/j.gedi.issn2095-8676.2023.06.006
    [6] 刘军伟, 梁展弘, 刘展志, 钟杰峰, 张紫凡.  统一潮流控制器注入电压及串联换流器容量优化计算研究 . 南方能源建设, 2023, 10(5): 157-165. doi: 10.16516/j.gedi.issn2095-8676.2023.05.019
    [7] 吴小刚, 何锐.  全桥MMC型直流融冰装置关键参数及装置容量解耦分析研究 . 南方能源建设, 2023, 10(5): 72-79. doi: 10.16516/j.gedi.issn2095-8676.2023.05.010
    [8] 黄华, 谢静媛, 汪正玲, 申葳, 张伟鹏.  基于时频分析法的开关柜避雷器泄漏电流检测 . 南方能源建设, 2022, 9(4): 150-158. doi: 10.16516/j.gedi.issn2095-8676.2022.04.019
    [9] 劳永钊, 吴任博, 肖健, 徐全, 陈吕鹏.  基于同步相量测量的配电网电压暂降溯源系统设计研究 . 南方能源建设, 2022, 9(1): 115-121. doi: 10.16516/j.gedi.issn2095-8676.2022.01.017
    [10] 王健声, 茆华风, 茆智伟, 许留伟.  整流变压器偏磁对托卡马克电源系统谐波不稳定的分析 . 南方能源建设, 2022, 9(2): 70-76. doi: 10.16516/j.gedi.issn2095-8676.2022.02.009
    [11] 余欣梅, 陈豪君, 王星华.  基于深度神经网络的DFIG低电压穿越技术研究 . 南方能源建设, 2021, 8(3): 122-130. doi: 10.16516/j.gedi.issn2095-8676.2021.03.018
    [12] 陈亮, 周伟, 杨源, 谭任深.  某越南电厂自动电压控制的实现及安全控制策略 . 南方能源建设, 2020, 7(S2): 96-100. doi: 10.16516/j.gedi.issn2095-8676.2020.S2.015
    [13] 姜绍艳, 李京平, 蔡志平.  智能分布式加电压-时间型双策略就地故障判别方法 . 南方能源建设, 2020, 7(S1): 48-52. doi: 10.16516/j.gedi.issn2095-8676.2020.S1.009
    [14] 刘子俊, 晋龙兴, 余锐, 陈永稳, 曹飞.  基于多维虚拟电流差动母线保护研究与应用 . 南方能源建设, 2019, 6(S1): 54-57. doi: 10.16516/j.gedi.issn2095-8676.2019.S1.011
    [15] 袁康龙.  运行方式对同塔线路的感应电压与电流影响分析 . 南方能源建设, 2018, 5(S1): 65-70. doi: 10.16516/j.gedi.issn2095-8676.2018.S1.012
    [16] 戴仲覆, 张英杰, 钟杰峰, 蒋泽甫, 钟以林.  贵州电网“十三五”静态电压稳定裕度的计算方法 . 南方能源建设, 2017, 4(4): 15-18. doi: 10.16516/j.gedi.issn2095-8676.2017.04.003
    [17] 邱国斌, 汪晶毅, 梁盼望, 王衍东.  同走廊内多回交流线路对在建线路的感应电压电流影响研究 . 南方能源建设, 2017, 4(4): 76-81. doi: 10.16516/j.gedi.issn2095-8676.2017.04.015
    [18] 徐伟, 廖民传, 郑志源, 汪晶毅.  500 kV海底电缆雷电过电压研究 . 南方能源建设, 2016, 3(2): 62-66. doi: 10.16516/j.gedi.issn2095-8676.2016.02.012
    [19] 郝为瀚, 简翔浩.  混合直流输电柔性直流换流站直流场接线研究 . 南方能源建设, 2016, 3(S1): 99-102. doi: 10.16516/j.gedi.issn2095-8676.2016.S1.021
    [20] 张青立.  Matlab/SPS与EMTP的操作过电压仿真分析及比较 . 南方能源建设, 2015, 2(S1): 35-37,92. doi: 10.16516/j.gedi.issn2095-8676.2015.S1.008
  • 期刊类型引用(2)

    1. 史小军,于铄航,王伟,公铮. 基于直挂储能控制的煤矿交直流混合配电网电压波动抑制研究. 工矿自动化. 2025(01): 78-84+94 . 百度学术
    2. 杨德林,谭令其,欧阳文艳,林曙光,宋啓能. 半桥型柔性直流换流器功率模块典型故障分析. 内蒙古电力技术. 2024(01): 84-90 . 百度学术

    其他类型引用(0)

  • 加载中
图(11) / 表 (1)
计量
  • 文章访问数:  450
  • HTML全文浏览量:  187
  • PDF下载量:  41
  • 被引次数: 2
出版历程
  • 收稿日期:  2023-03-20
  • 修回日期:  2023-06-26
  • 网络出版日期:  2023-09-06
  • 刊出日期:  2023-09-10

MMC子模块电容电压波动及谐波研究

DOI: 10.16516/j.gedi.issn2095-8676.2023.05.016
CSTR: 32391.14.j.gedi.issn2095-8676.2023.05.016
    基金项目:  国家电网有限公司总部管理科技项目资助“新能源经柔性直流送出系统在多类型谐波扰动下的穿越能力及控制策略研究”(5100-202256371A-2-0-ZN)
    作者简介:

    邹凯凯,1991-,男,工程师,华北电力大学电气工程及其自动化硕士,主要研究方向为直流输电仿真建模、控制保护策略开发(e-mail)zoukk@nrec.com

    李钢,1983-,男,高级工程师,电力系统及其自动化,主要从事柔性直流控制保护策略开发相关工作(e-mail)lig@nrec.com

    邹强,1981-,男,高级工程师,电力系统及其自动化,主要从事柔性直流控制保护策略开发相关工作(e-mail)zouq@nrec.com

    董云龙,1977-,男,教授级高级工程师,电力系统及其自动化,主要从事电力电子设备及控制保护研究相关工作(e-mail)dongyl@nrec.com

    卢宇,1979-,男,教授级高级工程师,电力系统及其自动化,主要从事电力电子设备及控制保护研究相关工作(e-mail)luy@nrec.com

    通讯作者: 邹凯凯,1991-,男,工程师,华北电力大学电气工程及其自动化硕士,主要研究方向为直流输电仿真建模、控制保护策略开发(e-mail)zoukk@nrec.com
  • 中图分类号: TM7;TM712

摘要:   目的  柔性直流的快速发展为MMC(Modular Multi-Level Converter,模块化多电平变流器)拓扑确定了实践支撑。文章以半桥MMC换流器为例,研究MMC稳态电容电压波动、谐波交互规律及桥臂开关函数的关系,以对MMC系统设计、控制策略给予理论支持。  方法  通过建立子模块电容电压、桥臂开关函数的解析表达式,得到MMC换流器上下桥臂能量波动及电容电压波动的数学近似表达式,进而研究了上下桥臂电流、交流阀侧、直流侧中各次谐波的分布规律及相间环流的正负序关系,基于谐波公式得到二次环流的估算方法;进一步,采用变量对照法研究了桥臂开关函数、子模块电容电压随阀侧功率因数、子模块电容大小的变化关系;最后基于实际工程参数搭建了详细的PSCAD/EMTDC模型,开展了理论解析与离线仿真的一致性对比。  结果  研究表明:MMC相间环流只含偶数次谐波,且其中6k+2次环流呈负序、6k+4次环流呈正序、6k次环流呈零序特性,同时稳态下MMC阀侧出口交流电压电流只含奇数次谐波,MMC直流侧电压电流只含偶数次谐波,且在功率因数角φ<0的某些范围内桥臂开关函数峰值会大于1,在低功率因数下子模块电容电压波动峰值大于高功率因数下的波动。  结论  可由此开展换流器电容容值的选型设计及暂稳态策略研究,并通过理论解析掌握MMC拓扑的运行本征特性。

English Abstract

邹凯凯, 李钢, 邹强, 董云龙, 卢宇. MMC子模块电容电压波动及谐波研究[J]. 南方能源建设, 2023, 10(5): 129-138. doi: 10.16516/j.gedi.issn2095-8676.2023.05.016
引用本文: 邹凯凯, 李钢, 邹强, 董云龙, 卢宇. MMC子模块电容电压波动及谐波研究[J]. 南方能源建设, 2023, 10(5): 129-138. doi: 10.16516/j.gedi.issn2095-8676.2023.05.016
ZOU Kaikai, LI Gang, ZOU Qiang, DONG Yunlong, LU Yu. Capacitor Voltage Ripple and Harmonics Analysis in MMC[J]. SOUTHERN ENERGY CONSTRUCTION, 2023, 10(5): 129-138. doi: 10.16516/j.gedi.issn2095-8676.2023.05.016
Citation: ZOU Kaikai, LI Gang, ZOU Qiang, DONG Yunlong, LU Yu. Capacitor Voltage Ripple and Harmonics Analysis in MMC[J]. SOUTHERN ENERGY CONSTRUCTION, 2023, 10(5): 129-138. doi: 10.16516/j.gedi.issn2095-8676.2023.05.016
    • 随着电力电子器件的发展,直流输电进入到了第三代新型直流输电,其基于全控型器件IGBT(Insulated-Gate Bipolar Transistor,绝缘栅双极晶体管)。采用全控型器件可以使换流器在四象限运行,实现有功无功解耦控制,因此特别适用于海上风电联网、无源网络供电、异步交流电网互联等。

      传统的两电平或三电平电压源换流器(Voltage Source Converter,VSC)由于谐波含量大、容量低、电压等级低、器件串联动态均压并联均流等问题,发展受到了制约。2001年德国学者提出了模块化多电平换流器MMC(Modular Multilevel Converter,模块化多电平变流器)概念[1-3],西门子公司于2010年将其用于实际工程。近年来,基于MMC的工程越来越多,据不完全统计,目前已投运30项工程。国外工程集中在欧洲,主要将其用于海上风电的送出和电网互联。近年来国内已有南汇、厦门、舟山、鲁西、渝鄂、张北柔性直流电网、乌东德混合直流、三峡如东海风柔性直流、白鹤滩混合直流等工程。

      文献[4-6]给出了MMC的基本数学关系和静态特性;文献[7]研究了采用PWM( Pulse Width Modulation,脉宽调制)方式的谐波特性;文献[8-10]对MMC的动态特性进行了建模研究;文献[11-13]研究了MMC联接弱交流电网的稳定运行区间;文献[14-16]开展了MMC孤岛送新能源场景下的控制研究;文献[17-18]就MMC处于直流电压控制下的稳定性开展了研究分析;文献[19-21]就MMC从电容电压控制出发开展的控制手段进行了研究。

      近年来,对于MMC的研究主要集中于MMC与传统LCC(Line-Commuted Converter,换相换流器)直流的混合特性、主动故障隔离及故障清除策略、接入弱电网时的稳定性分析、新能源经柔性直流送出的稳定性分析、自适应强弱电网等方向。

      其中对于柔性直流稳定性的分析成为重中之重,当前的分析手段是针对柔性直流系统的稳态工作点建立小信号模型,基于小信号模型得到柔性直流系统的阻抗特性,再根据柔性直流及外部电网的阻抗特性进行稳定性分析[22-24]。因此,对MMC稳态建模是稳定性分析的基石。

      此外,高压大容量MMC换流器设备造价主要取决于IGBT器件及并联电容器的造价。当系统容量一定时,IGBT器件选型较为固定,并联电容器的大小选型的重要参考是子模块电容电压的波动。因此,对MMC的稳态工作点的精细化建模特别是模块电容电压分析具有重要意义。已发表的相关文献侧重点在于MMC内部环流产生原因及环流抑制的相关策略,文章则侧重电容能量变化规律,从子模块电容能量变化角度出发给出一种精细化数学建模,在推导出桥臂电容能量的波动公式前提下,深入研究了桥臂电容电压随时间、功率因数角、子模块电容大小的波动规律,根据电容电压的波动得到MMC谐波分布特性,并得出二次谐波的估算公式;在电容电压波动前提下,研究了桥臂子模块投入个数和模块电容电压随时间、功率因数的关系;最后对得到的规律进行了详细的仿真验证。

    • 典型的半桥MMC换流器拓扑如图1所示,各相单元含上下两个桥臂,各桥臂单元由N个MMC子模块及桥臂电抗L组成,图中R为子模块等效内阻。子模块采用半桥结构,内含两个具有反并联二极管的IGBT和电容组成,实际工程中每个子模块还含一个保护晶闸管和快速开关。在直流侧通过钳位电容或钳位大电阻接地提供接地钳位点时[25],此时正负极母线上电压分别为Vd/2及−Vd/2。

      图  1  MMC换流器拓扑图

      Figure 1.  Topological structure of MMC

    • 对相单元与交直流侧形成回路的利用基尔霍夫电压定律( Kirchhoff's Voltage Law,KVL)有:

      {Vd2vuRiuLdiudt=vaVd2+vl+Ril+Ldildt=va (1)

      式中:

      L ——桥臂电抗值(H);

      R ——子模块的等效内阻(Ω);

      Vd ——直流电压(kV);

      vuvl ——上、下桥臂端间电压(kV);

      iuil ——上、下桥臂电流(kA),电流方向如图1所示;

      va ——阀侧A相相电压(kV)。

      对式(1)进行数学变换得:

      {L2disdt=vu+vl2vsvaR2isLdicdt=vu+vl2vc+Vd2Ric (2)

      式中:

      is ——阀侧相电流(kA);

      ic ——流过桥臂的直流电流分量(kA);

      vs ——上下桥臂的差模电压(kV);

      vc ——上下桥臂的共模电压(kV)。

      式(2)中用到了如下关系[26]

      {iu=is2+icil=is2+ic (3)

      根据每个子模块投入切除状态定义单个子模块jj=1,2,···,N)投入切除的开关状态数nju,l(为1表示投入,为0表示切除),在单桥臂中所有子模块电压均衡时,可得到单个桥臂总投入子模块等效电压为:

      vcu,l=Nj=1nju,lvjcu,l=vΣcu,lNNj=1nju,l=nu,lvΣcu,l (4)

      式中:

      vΣcu,l ——上/下桥臂所有子模块电压之和(kV);

      N ——上/下桥臂子模块个数;

      nu,l ——上/下桥臂的总开关函数,nu为上桥臂的总开关函数,nl为下桥臂的总开关函数。

    • 考虑每个子模块电容电压波动后将其叠加得到桥臂上所有MMC子模块两端电压波动,以上桥臂为例,有:

      Cdvjcudt=njuiu (5)
      CNj=1dvjcudtdvΣcu/dt=Nj=1njuiu=iuNj=1njuNnu (6)

      式中:

      C ——MMC子模块电容容值(F)。

      替换下标u为表示下桥臂的l,公式依然成立。将式(6)重写得到上下桥臂的总子模块电容电压波动方程如下:

      {CNdvΣcudt=nu(is2+ic)CNdvΣcldt=nl(is2+ic) (7)

      将式(4)代入到vsvc的公式中,可得到桥臂的总开关函数为:

      nu=vcvsvΣcunl=vc+vsvΣcl (8)

      结合式(7)和式(8),可以得到上下桥臂电容能量的波动方程:

      {vΣcuCNdvΣcudt=C2Nd(vΣcu)2dtdWu/dt=(vcvs)(is2+ic)vΣclCNdvΣcldt=C2Nd(vΣcl)2dtdWl/dt=(vc+vs)(is2+ic) (9)

      从式(9)进一步得到某相单元中上下桥臂总电容能量的波动及上下桥臂电容能量之差的波动方程为:

      {WΣ=Wu+WlWΔ=WuWl{dWΣdt=2vcicvsisdWΔdt=vcis2vsic (10)

      表达式中vc=Vd/2;vs=Vscos(w1t);is=Iscos(w1tφ)。

      将上述关系代入式(10)中,可得:

      {dWΣdt=VdicVsIs2cosφVsIs2cos(2w1tφ)dWΔdt=VdIs2cos(w1tφ)2Vsiccos(w1t) (11)

      从式(11)的第一个公式可以看出,相单元的总能量变化主要由三部分构成:第一部分为直流分量对应到单相中的有功功率;第二部分为交流系统注入MMC换流器的有功功率;第三部分为表示相间交换的有功,其以二倍频为主,表示相间环流以二倍频为主。

      将式(11)两边积分得:

      {WΣ=W0VsIs4w1sin(2w1tφ)ΔWΣWΔ=WΔ0+VdIs2w1sin(w1tφ)2Vsicw1sin(w1t)ΔWΔ (12)

      式中:

      ΔWΣ ——相单元总能量的波动量(MJ);

      ΔWΔ ——上下桥臂能量差的波动量(MJ)。

      由于能量储存于子模块电容中,假设不考虑实际工程中的子模块冗余个数,则有:

      {WΣ0=2N×12C(VdN)2=CV2dNWΔ0=0 (13)

      将式(13)代入式(12)可得到上下桥臂总子模块电容能量为:

      {Wu=WΣ+WΔ2=WΣ0+ΔWΣ+ΔWΔ2Wl=WΣWΔ2=WΣ0+ΔWΣΔWΔ2 (14)

      根据桥臂能量与电容电压之间的关系,结合式(12)~式(14)可最终得到上下桥臂MMC子模块电容电压的详细波动关系如下:

      {vΣcu=2NCWu=Vd1+NCV2d(ΔWΣ+ΔWΔ)Vd+N2CVd(ΔWΣ+ΔWΔ)ΔvΣcuvΣcl=2NCWl=Vd1+NCV2d(ΔWΣΔWΔ)Vd+N2CVd(ΔWΣΔWΔ)ΔvΣcl (15)

      式中:

      ΔvΣcuΔvΣcl——上、下桥臂子模块电容电压之和的波动量(kV)。

      为定性初步分析基波分量与二倍频分量的规律,在泰勒展开时忽略了更高频分量,从式(15)可初步得到如下结论:

      1)桥臂电容电压和的波动与电容C成反比。

      2)上下桥臂的电容电压和的波动均含两部分,一部分为ΔWΣ二倍频分量,一部分为ΔWΔ基频分量。

      3)二倍频分量在上下桥臂总大小相等,方向相同;基波分量在上下桥臂上大小相等,方向相反。

    • 将式(15)泰勒展开后取其任一项有:

      {(ΔWΣ+ΔWΔ)n=nk=0CknΔWΣkΔWΔnk(ΔWΣΔWΔ)n=nk=0CknΔWΣk(ΔWΔ)nk (16)

      对典型项CknΔWΣkΔWΔnkCknΔWΣk(ΔWΔ)nk进行分析,由于ΔWΣ为二倍频分量,​ΔWΔ为基波分量,可以看出ΔWΣk项仍为偶数倍频,谐波电压的奇偶关系与ΔWΔnk(ΔWΔ)nkn−k的奇偶性相同。因此,可以得到上下桥臂的电容电压波动量的奇数次谐波呈谐波大小相等,方向相反的关系;偶数次谐波大小相等,方向相同的关系。由于电容电压的波动会带来桥臂上的谐波电流,因此,我们可以得出桥臂电流上谐波电流分布规律:

      1)上下桥臂的奇数次谐波大小相等,方向相反。

      2)上下桥臂的偶数次谐波大小相等,方向相同。

      3)相间环流为偶数次环流,且以二倍频为主。

      4)第1)条等同于MMC阀侧出口交流电压电流只含奇数次谐波。

      5)第2)条等同于MMC直流侧电压电流只含偶数次谐波。

      6)根据ΔWΣ和​ΔWΔ的计算公式,可以得到相间环流只含偶数次谐波,其中6k+2​次环流呈负序;6k+4​次环流呈正序;6k​次环流呈零序[12]

      从式(17)、式(18)的电容电压波动泰勒展开式出发,考虑取2w1t的相关部分后可以达到二次谐波电流的估算大小。由于高频分量在进行三角函数计算时也有可能产生二倍频分量,这里忽略这些高频分量因中间计算过程带来的二倍频分量的大小,忽略的原因在于随着展开的幂次数的增加,其分量的幅值越来越小,图2给出了泰勒展开式的系数函数随展开阶数的关系。

      图  2  泰勒展开式系数函数

      Figure 2.  Coefficient function of Talyor Series

      vΣcu=vd1+NCv2d(ΔWΣ+ΔWΔ)=vd+N2Cvd(ΔWΣ+ΔWΔ)++12×(121)×(122)×(12k+1)vdk!×[NCv2d(ΔWΣ + ΔWΔ)]k+ (17)
      vΣcl=vd1+NCv2d(ΔWΣΔWΔ)=vd+N2Cvd(ΔWΣΔWΔ)++12×(121)×(122)××(12k+1)k!×vd[NCv2d(ΔWΣΔWΔ)]k+ (18)

      得到二次谐波电流的估算大小为:

      i2harm=CNdΔvcudt=Asin(2w1t)+Bcos(2w1t)=A2+B2cos(2w1t+φ) (19)
      A=VsIs4VdsinφNV2si2c2w1CV3dNV2dI2s32w1CV3dcos(2φ)+NVdicVsIs8w1CV3dcosφ (20)
      B=VsIs4VdcosφNV2dI2s32w1CV3dsin(2φ)+NVdicVsIs8w1CV3dsinφ (21)
    • 式(22)、式(23)给出了上、下桥臂子模块电容电压的波动规律,从式子可以知道在忽略高频谐波的情况下桥臂电容电压和的波动与电容C成反比,在直流电压、直流电流、阀侧相电压相电流幅值一定的情况下,子模块电容电压的波动主要取决于阀侧功率因数角φ和时间t。为深入研究,本章节以我国首个海风柔性直流江苏如东工程参数为算列开展,其主参数如表1所示。

      表 1  MMC参数

      Table 1.  Parameters for MMC

      参数名
      额定直流电压/kV ±400
      额定直流电流/kA 1.375
      额定有功/GW 1.100
      额定视在功率/MVA 1230
      子模块个数不考虑冗余/个 400
      子模块电容容值/mF 4, 6, 8, 9, 10, 12
      桥臂电抗器/H 0.133
      阀侧相电压峰值/kV 340
      阀侧相电流峰值/kA 2.411
      vΣcu=vd1+NCv2d(ΔWΣ+ΔWΔ)=vd1+NCv2d(VsIs4w1sin(2w1tφ)+vdIs2w1sin(w1tφ)2Vsicw1sin(w1t)) (22)
      vΣcl=vd1+NCv2d(ΔWΣΔWΔ)=vd1+NCv2d(VsIs4w1sin(2w1tφ)vdIs2w1sin(w1tφ)+2Vsicw1sin(w1t)) (23)
    • 根据式(22)、式(23),假设阀侧电压与阀侧电流的相位差在[−π, π]之间变动,考虑4个工频周期下的波动,得到如下三维关系图,其中x轴坐标为时间t(单位为s),y轴坐标为相角差φ(单位为rad),z轴坐标为桥臂子模块电容电压的波动量ΔvΣcu(单位为kV),根据前面式子,得到上下桥臂所有投入的子模块电容电压的波动随时间及阀侧电压初始相位差的三维关系图,如图3图4所示。

      图  3  上桥臂子模块电容电压的波动

      Figure 3.  Capacitor voltage ripple of upper bridge

      图  4  下桥臂子模块电容电压的波动

      Figure 4.  Capacitor voltage ripple of lower bridge

      从上图可知,在φ=0时,MMC换流阀与交流侧没有无功交换,此时桥臂子模块电容电压波动的峰值最小。在φ减小或增大的过程中,波动量均在一定程度上增大。

      根据式(8)同样可得到上下桥臂的总开关函数nunl随时间及功率因数角φ的曲线图。其中x轴坐标为时间t(单位为s),y轴坐标为相角差φ(单位为rad),z轴坐标为桥臂总子模块开关函数,图5为上桥臂总开关函数图,图6为下桥臂总开关函数图。

      图  5  上桥臂总开关函数

      Figure 5.  Switching function of upper bridge

      图  6  下桥臂总开关函数

      Figure 6.  Switching function of lower bridge

      图5图6可以明显看出,功率因数角φ一定时桥臂总开关函数是时间t的三角工频周期函数,在时间t一定时,开关函数是功率因数角的周期函数,且在φ<0的一定范围内,开关函数的最大值大于1,这是由于φ<0时,换流器向交流系统注入容性无功,在[−90°, 0]范围内,φ越小注入容性无功越大,此时同等有功条件下需要投入的子模块个数越多。将式(15)代入式(8)以后,可得到上下桥臂的开关函数的表达式,这里仅列出下桥臂开关函数:

      nl=Vd/2+Vscosw1tVd+ΔvΣcl (24)
      ΔvΣcl=N2Cvd[VsIs4w1sin(2w1tφ)vdIs2w1sin(w1tφ)+2Vsicw1sin(w1t)] (25)

      t=0时刻,有桥臂最大的开关函数为:

      nl|max=Vd/2+VsVd+NIs4w1C(Vs2Vd+1)sinφ (26)

      在开关函数最大值为1的情况下,此时得出:

      sinφ=4w1CNIs×VsVd2Vs2Vd+1 (27)

      由于阀侧相电压最大值小于直流电压的一半,可以看出在开关函数的最大值为1时,功率因数角φ∈[−π/2, 0],这与图中结论一致。

      为进一步分析开关函数与模块电容电压随功率因数角的变换关系,给出了不同功率因数角下的开关函数及桥臂子模块电容电压波动的关系图(图7)。从图中及式(24)可得到如下结论:

      图  7  不同功率因数角下桥臂开关函数及模块电压波动

      Figure 7.  Switching function and capacitor voltage ripple for different phi

      1)桥臂开关函数的最大值在φ∈[−π/2, π/2]范围内单调递减,且在[−π/2, 0]内某点φ0处其最大值为1。

      2)在额定视在功率一定的情况下,桥臂开关函数的最大值在φ<0的某些范围内会呈现大于1的情况,此时功率因数较小,理论交换无功较大,φ=−90°时换流器工作在STATCOM工况下,实际工程中采用限幅环节对开关函数进行限幅。

      3)子模块电容电压纹波存在明显的二倍频和基波分量。

      4)在低功率因数φ=±90°下,子模块电容电压纹波比高功率因数下的纹波更大,且在低功率因数下子模块电容电压波动呈现正负半周期不对称现象。

    • 考虑φ=0时,子模块电容电压及桥臂开关函数随电容容值的变化关系。图8给出了表1参数下的桥臂开关函数及电容电压波动随时间和电容容值的变化。

      图  8  不同子模块电容容值下桥臂开关函数及模块电压波动

      Figure 8.  Switching function and capacitor voltage ripple for different C

      图8中可知,在φ=0时,桥臂开关函数的最大值不随电容电压的变化而变化,进一步从式(20)中可知,在φ≠0时,开关函数的最大值一定程度上随电容容值的变大而变小。子模块电容电压的波动量随着电容容值的变大而减少。

    • 为验证上述结论,在PSCAD/EMTDC仿真软件中搭建了上述表1参数下的仿真模型。

      为验证子模块电容电压波动随子模块电容容值的变化,图9给出了计算与仿真的对比关系图,其横坐标为子模块电容容值,纵坐标为电容电压的波动幅值相对基准值的占比,可看出子模块电容电压波动幅值随容值的变大而减小,理论计算值与仿真值较为吻合。

      图  9  电容电压波动幅值随电容的变化关系

      Figure 9.  Capacitor voltage ripple for different capacitance

      为验证子模块电容电压波动实时曲线及桥臂开关函数的实时曲线,在仿真模型中任取一个PQ功率点(1100 MW, 123 Mvar)下对两者进行了对比。从图10对比中可看出在任取的功率点下,上下桥臂开关函数理论计算值与仿真值吻合度较高。图11给出了上下桥臂子模块电容电压实时波动的对比,可看出计算与仿真总体较为吻合,在模块电容电压波动过零点处的不一致主要是由于计算中未考虑二次以上的高倍频谐波波动导致。

      图  10  实时桥臂开关函数计算与仿真对比

      Figure 10.  Switching function comparison between calculation and simulation

      图  11  模块总电容电压实时波动对比

      Figure 11.  Capacitor voltage ripple comparison between calculation and simulation

    • 文章从能量角度出发得到了桥臂电容电压、桥臂开关的理论计算值,在此基础上研究了桥臂谐波的特性,并深入研究了子模块电容电压的波动和桥臂开关函数的情况,主要得出的结论与仿真结果一致,为实际工程中MMC换流器设计提供了参考。

      但是,文章对换流器建模尚未结合换流器自身控制系统,后续可开展结合控制系统的换流器阻抗建模研究及换流器自身调制环节、控制链路等MMC稳定性影响因子分析;另一方面,通过注入倍频环流或先进控制方式可降低模块电容电压的波动,以优化系统参数的选型,这也是后续待研究之一。

参考文献 (26)

目录

/

返回文章
返回