CN105226981B - 一种具有中点电位平衡控制的空间矢量调制方法 - Google Patents

一种具有中点电位平衡控制的空间矢量调制方法 Download PDF

Info

Publication number
CN105226981B
CN105226981B CN201510701913.0A CN201510701913A CN105226981B CN 105226981 B CN105226981 B CN 105226981B CN 201510701913 A CN201510701913 A CN 201510701913A CN 105226981 B CN105226981 B CN 105226981B
Authority
CN
China
Prior art keywords
mrow
msub
mtd
mtr
mfrac
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.)
Expired - Fee Related
Application number
CN201510701913.0A
Other languages
English (en)
Other versions
CN105226981A (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.)
Suzhou Power Supply Co of State Grid Jiangsu Electric Power Co Ltd
Original Assignee
China University of Mining and Technology CUMT
Suzhou Power Supply Co of State Grid Jiangsu 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 China University of Mining and Technology CUMT, Suzhou Power Supply Co of State Grid Jiangsu Electric Power Co Ltd filed Critical China University of Mining and Technology CUMT
Priority to CN201510701913.0A priority Critical patent/CN105226981B/zh
Publication of CN105226981A publication Critical patent/CN105226981A/zh
Application granted granted Critical
Publication of CN105226981B publication Critical patent/CN105226981B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
  • Inverter Devices (AREA)

Abstract

一种具有中点电位平衡控制的空间矢量调制方法,属于空间矢量调制方法。该方法首先对传统虚拟空间矢量进行改进在虚拟小矢量和虚拟中矢量中分别引入调节因子,然后针对虚拟空间矢量在三电平变流器中点电位出现直流不平衡时无法自行恢复的缺点,在调制策略中加入一种基于模型预测的多目标优化直流不平衡控制策略,以实现虚拟空间矢量对三电平变流器中点电位的完全控制。改进型虚拟空间矢量调制方法,克服了传统虚拟空间矢量无法抑制中点电位直流不平衡的缺点,所加入的多目标优化控制策略很好的解决了直流不平衡恢复过程中恢复速度和品质之间的矛盾。同时改进型虚拟空间矢量在稳态下与传统虚拟空间矢量一样可以消除中点电位低频振荡。

Description

一种具有中点电位平衡控制的空间矢量调制方法
技术领域
本发明涉及一种空间矢量调制方法,特别是一种具有中点电位平衡控制的空间矢量调制方法。
背景技术
脉冲宽度调制技术(PWM)是变频器控制领域中的核心技术,调制策略的好坏将直接影响变频器的各项性能指标。现有的脉冲宽度调制技术大致可以分为两类:一种是基于载波比较的调制策略(CB-PWM),另一类是基于空间矢量合成的调制策略(SVM)。两类调制策略各自具有鲜明的优缺点,基于载波的调制技术算法简单硬件资源开销较少,但是其直流母线电压利用率较低且交流输出端电压畸变率较高;而基于空间矢量的调制技术却拥有相对较高的直流母线电压利用率和较低的交流输出端电压畸变率,但是其算法复杂较高硬件开销大。随着微控制芯片制造技术的发展,微控制器的运算速度与运算能力得到了大幅度的提升同时制造成本大幅降低,制约空间矢量调制技术发展的高算法复杂度得到了很好地解决,因此具有良好性能的空间矢量技术在变频器控制领域得到了广泛的应用。
在中压大功率应用场合中点电位钳位型三电平拓扑结构是一种成熟的解决方案,但是三电平拓扑也有其自身的缺陷,即中点电位在变频器运行过程中会出现偏差。这种中点电位偏差主要包括两个部分:第一,稳态情况下的中点电位直流不平衡;第二,暂态情况下的中点电位低频振荡(三倍于基波频率)。传统三电平空间矢量技术能够很好的实现变频器输出电压控制,但是当基于传统三电平空间矢量调制技术控制的三电平逆变器工作在高调制度和低功率因数的工况下,三电平变频器的中点电位会出现严重的中点电位低频振荡。而且随着调制度的增加和功率因数的降低这种中点电位低频震荡会变得更加严重,不仅会影响变频器的控制输出效果严重的甚至还会威胁设备安全。
虚拟空间矢量是一种从传统空间矢量衍生出来的新型矢量调制技术,虚拟空间矢量很好的解决了中点电位低频振荡问题,在全范围内实现了中点电位可控。然而虚拟空间矢量在每一个矢量作用时间内都保证平均中点电位为零,当中点电位因某种原因产生直流不平衡时虚拟空间矢量调制技术将无法自行恢复这种中点电位直流不平衡。
发明内容
本发明的目的是要提供一种具有中点电位平衡控制的空间矢量调制方法,实现对三电平变频器的输出脉冲宽度控制和中点电位两个方面问题的完全控制。
本发明的目的是这样实现的:该方法在传统虚拟空间矢量的基础上分别针对虚拟小矢量和虚拟中矢量的分配系数引入动态调整因子,同时建立中点电位预测模型;根据中点电位的反馈信号,判断当前变频器运行状态下是否存在中点电位直流不平衡,然后将反馈信号和预测模型的输出结果同时送入价值判断函数,计算价值函数最小的动态调节因子;当价值函数取值最小时获得最优的虚拟矢量动态调整因子,修正后的调整因子用来生成新的输出脉冲,实现中点电位的完全控制。
在传统虚拟空间矢量的基础上分别针对虚拟小矢量和虚拟中矢量的分配系数引入动态调整因子;具体步骤如下:
步骤a.在传统虚拟空间矢量按照矢量的长度分为四类:虚拟零矢量VVZ,虚拟小矢量VVSx,虚拟中矢量VVMx,虚拟大矢量VVLx,其中x=1,2,3,4,5,6;在虚拟中矢量和虚拟小矢量中引入分配系数动态调整因子,式(1)为传统虚拟空间矢量,式(2)为改进型虚拟空间矢量;
式中,KMx,KSx∈[-1,1],VMx,VLx分别是三电平空间矢量平面中的基本中矢量、基本大矢量、负组基本小矢量和正组基本小矢量,VNNN,VOOO,VPPP三者为三电平空间矢量平面中的三个零矢量;改进虚拟中矢量将传统虚拟矢量分割成两个部分:第一,基本中矢量部分;第二,相邻正负基本小矢量合成部分;在改进型虚拟矢量中无论引入的调节因子如何变化,都不影响最终合成的虚拟空间矢量的相位和幅值;
步骤b.根据边界条件确定参考电压矢量所在区域,所述的区域包括空间矢量平面大扇区和扇区内的三角形区域;然后选择构成三角形三个顶点所对应的虚拟空间矢量来合成参考电压矢量;根据伏秒平衡原理通过式(3)计算出被选中的三个虚拟空间矢量的作用时间;
步骤c.依据基于模型预测的多目标优化中点电位直流不平衡控制策略,实时计算出动态调节因子K,送入步骤b重新调整各矢量的时间分配,消除中点电位直流不平衡。
所述的计算价值函数最小的动态调节因子;具体步骤如下:
步骤a.根据直流侧电路简化模型建立中点电位离散数学模型,根据基尔霍夫电流定律流入三电平变频器中点的电流满足如下关系:
将式(4)中第1,2式带入第3式中得到其中vNP=uC1-uC2,vNP0为中点电位初始时刻偏差值;通常采样频率远高于工频时,用离散求和的方式来代替积分项,中点电位的离散模型表示为
其中,im(k)为第k个采样时刻m相所对应的相电流值,tm0(k)为在第k个采样周期内m相电流对应流入变频器中点的时间;获得中点电位的变化率,以式(6)表示;
步骤b.建立目标函数如下
其中,λ为权重因子,为中点电位理想参考值通常为零;目标函数由两部分构成:a)为第一部分,当第一部分满足取值最小时可以保证在下一个采样时刻中点电位以最快的速度趋近于给定参考值b)λ·vNP(k)·ΔvNP(k)2为第二部分,当第二部分取最小值是可以保证直流不平衡恢复过程最为平稳;通过选取不同的调节因子可以在速度和品质之间获得平衡;
步骤c.将当前时刻的中点电位值中点电位变化率送入目标函数中,解得使目标函数获得最小值时所对应的K值的最优解为
其中,tm0_max(k)为K取边界值时所对应的作用时间,sgn(*)为符号函数。
有益效果,由于采用了上述方案,三电平钳位型变频器中点电位完全控制,就是消除暂态不平衡中的低频振荡和抑制稳态不平衡中的直流不平衡;在传统虚拟空间矢量的基础上对于组成虚拟空间的固定分配比例做出修改,在虚拟小矢量和虚拟中矢量中引入可以动态调节分配比例的调节因子,使得改进的虚拟空间矢量同时具备了抑制中点电位低频振荡和直流不平衡的能力。在此基础上同时加入了一种基于模型预测的多目标优化控制环节,当中点电位出现直流不平衡的时候,多目标优化策略很好的解决了直流不平衡恢复过程中的速度和品质的矛盾。因此,是一种具备中点电位完全控制能力的空间矢量调制方法。该虚拟空间矢量调制方法需要对变频器直流侧上下电容的电压同步采样。
优点:该改进型虚拟空间矢量调制方法能确使得钳位型三电平变频器中点电位得到完全控制,即同时消除暂态不平衡中的低频振荡和抑制稳态不平衡中的直流不平衡,适用于各类钳位型三电平拓扑结构的变频器、变流器、逆变器等电力电子电能变化装置的脉冲宽度调制控制。
附图说明
图1为本发明三电平基本空间示意图。
图2为本发明三电平虚拟空间矢量示意图。
图2(a)为本发明全部虚拟空间矢量在整个空间矢量平面内的分布图。
图2(b)为本发明部分虚拟空间矢量在空间矢量平面第一扇区内的分布图以及判断参考电压矢量落点的区域划分(第一扇区内)示意图。
图3基于多目标优化的改进型虚拟空间矢量的三电平变频器系统框图。
图4(a)为采用本发明的改进型虚拟空间矢量控制的三电平变频器中点电位直流不平衡恢复过程。
图4(b)为采用本发明的改进型虚拟空间矢量控制的三电平变频器中点电位直流不平衡恢复过程细节。
具体实施方式
实施例1:该方法在传统虚拟空间矢量的基础上分别针对虚拟小矢量和虚拟中矢量的分配系数引入动态调整因子,同时建立中点电位预测模型;根据中点电位的反馈信号,判断当前变频器运行状态下是否存在中点电位直流不平衡,然后将反馈信号和预测模型的输出结果同时送入价值判断函数,计算价值函数最小的动态调节因子;当价值函数取值最小时获得最优的虚拟矢量动态调整因子,修正后的调整因子用来生成新的输出脉冲,实现中点电位的完全控制。
在传统虚拟空间矢量的基础上分别针对虚拟小矢量和虚拟中矢量的分配系数引入动态调整因子;具体步骤如下:
步骤a.在传统虚拟空间矢量按照矢量的长度分为四类:虚拟零矢量VVZ,虚拟小矢量VVSx,虚拟中矢量VVMx,虚拟大矢量VVLx,其中x=1,2,3,4,5,6;在虚拟中矢量和虚拟小矢量中引入分配系数动态调整因子,式(1)为传统虚拟空间矢量,式(2)为改进型虚拟空间矢量;
式中,KMx,KSx∈[-1,1],VMx,VLx分别是三电平空间矢量平面中的基本中矢量、基本大矢量、负组基本小矢量和正组基本小矢量,VNNN,VOOO,VPPP三者为三电平空间矢量平面中的三个零矢量;改进虚拟中矢量将传统虚拟矢量分割成两个部分:第一,基本中矢量部分;第二,相邻正负基本小矢量合成部分。在改进型虚拟矢量中无论引入的调节因子如何变化,都不影响最终合成的虚拟空间矢量的相位和幅值;
步骤b.根据边界条件确定参考电压矢量所在区域,所述的区域包括空间矢量平面大扇区和扇区内的三角形区域;然后选择构成三角形三个顶点所对应的虚拟空间矢量来合成参考电压矢量;根据伏秒平衡原理通过式(3)计算出被选中的三个虚拟空间矢量的作用时间;
步骤c.依据基于模型预测的多目标优化中点电位直流不平衡控制策略,实时计算出动态调节因子K,送入步骤b重新调整各矢量的时间分配,消除中点电位直流不平衡。
所述的计算价值函数最小的动态调节因子;具体步骤如下:
步骤a.根据直流侧电路简化模型建立中点电位离散数学模型,根据基尔霍夫电流定律流入三电平变频器中点的电流满足如下关系:
将式(4)中第1,2式带入第3式中得到其中vNP=uC1-uC2,vNP0为中点电位初始时刻偏差值;通常采样频率远高于工频时,用离散求和的方式来代替积分项,中点电位的离散模型表示为
其中,im(k)为第k个采样时刻m相所对应的相电流值,tm0(k)为在第k个采样周期内m相电流对应流入变频器中点的时间;获得中点电位的变化率,以式(6)表示;
步骤b.建立目标函数如下
其中,λ为权重因子,为中点电位理想参考值通常为零。目标函数由两部分构成:a)为第一部分,当第一部分满足取值最小时可以保证在下一个采样时刻中点电位以最快的速度趋近于给定参考值b)λ·vNP(k)·ΔvNP(k)2为第二部分,当第二部分取最小值是可以保证直流不平衡恢复过程最为平稳;通过选取不同的调节因子可以在速度和品质之间获得平衡;
步骤c.将当前时刻的中点电位值中点电位变化率送入目标函数中,解得使目标函数获得最小值时所对应的K值的最优解为
其中,tm0_max(k)为K取边界值时所对应的作用时间,sgn(*)为符号函数。

Claims (2)

1.一种具有中点电位平衡控制的空间矢量调制方法,其特征是:该方法在传统虚拟空间矢量的基础上分别针对虚拟小矢量和虚拟中矢量的分配系数引入动态调整因子,同时建立中点电位预测模型;根据中点电位的反馈信号,判断当前变频器运行状态下是否存在中点电位直流不平衡,然后将反馈信号和预测模型的输出结果同时送入价值判断函数,计算价值函数最小的动态调节因子;当价值函数取值最小时获得最优的虚拟矢量动态调整因子,修正后的调整因子用来生成新的输出脉冲,实现中点电位的完全控制;
在传统虚拟空间矢量的基础上分别针对虚拟小矢量和虚拟中矢量的分配系数引入动态调整因子;具体步骤如下:
步骤a.在传统虚拟空间矢量按照矢量的长度分为四类:虚拟零矢量VVZ,虚拟小矢量VVSx,虚拟中矢量VVMx,虚拟大矢量VVLx,其中x=1,2,3,4,5,6;在虚拟中矢量和虚拟小矢量中引入分配系数动态调整因子,式(1)为传统虚拟空间矢量,式(2)为改进型虚拟空间矢量;
<mrow> <mtable> <mtr> <mtd> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mi>V</mi> <mrow> <mi>V</mi> <mi>M</mi> <mi>x</mi> </mrow> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mn>3</mn> </mfrac> <mrow> <mo>(</mo> <msub> <mi>V</mi> <mrow> <mi>M</mi> <mi>x</mi> </mrow> </msub> <mo>+</mo> <msubsup> <mi>V</mi> <mrow> <mi>S</mi> <mi>x</mi> </mrow> <mo>-</mo> </msubsup> <mo>+</mo> <msubsup> <mi>V</mi> <mrow> <mi>S</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> <mo>+</mo> </msubsup> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>V</mi> <mrow> <mi>V</mi> <mi>S</mi> <mi>x</mi> </mrow> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <mrow> <mo>(</mo> <msubsup> <mi>V</mi> <mrow> <mi>S</mi> <mi>x</mi> </mrow> <mo>+</mo> </msubsup> <mo>+</mo> <msubsup> <mi>V</mi> <mrow> <mi>S</mi> <mi>x</mi> </mrow> <mo>-</mo> </msubsup> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>V</mi> <mrow> <mi>V</mi> <mi>L</mi> <mi>x</mi> </mrow> </msub> <mo>=</mo> <msub> <mi>V</mi> <mrow> <mi>L</mi> <mi>x</mi> </mrow> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>V</mi> <mrow> <mi>V</mi> <mi>Z</mi> </mrow> </msub> <mo>=</mo> <msub> <mi>V</mi> <mrow> <mi>N</mi> <mi>N</mi> <mi>N</mi> </mrow> </msub> <msub> <mi>orV</mi> <mrow> <mi>O</mi> <mi>O</mi> <mi>O</mi> </mrow> </msub> <msub> <mi>orV</mi> <mrow> <mi>P</mi> <mi>P</mi> <mi>P</mi> </mrow> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> </mtd> <mtd> <mrow> <mi>x</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>3</mn> <mo>,</mo> <mn>5</mn> </mrow> </mtd> <mtd> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mi>V</mi> <mrow> <mi>V</mi> <mi>M</mi> <mi>x</mi> </mrow> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mn>3</mn> </mfrac> <mrow> <mo>(</mo> <msub> <mi>V</mi> <mrow> <mi>M</mi> <mi>x</mi> </mrow> </msub> <mo>+</mo> <msubsup> <mi>V</mi> <mrow> <mi>S</mi> <mi>x</mi> </mrow> <mo>+</mo> </msubsup> <mo>+</mo> <msubsup> <mi>V</mi> <mrow> <mi>S</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> <mo>-</mo> </msubsup> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>V</mi> <mrow> <mi>V</mi> <mi>S</mi> <mi>x</mi> </mrow> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <mrow> <mo>(</mo> <msubsup> <mi>V</mi> <mrow> <mi>S</mi> <mi>x</mi> </mrow> <mo>+</mo> </msubsup> <mo>+</mo> <msubsup> <mi>V</mi> <mrow> <mi>S</mi> <mi>x</mi> </mrow> <mo>-</mo> </msubsup> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>V</mi> <mrow> <mi>V</mi> <mi>L</mi> <mi>x</mi> </mrow> </msub> <mo>=</mo> <msub> <mi>V</mi> <mrow> <mi>L</mi> <mi>x</mi> </mrow> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>V</mi> <mrow> <mi>V</mi> <mi>Z</mi> </mrow> </msub> <mo>=</mo> <msub> <mi>V</mi> <mrow> <mi>N</mi> <mi>N</mi> <mi>N</mi> </mrow> </msub> <msub> <mi>orV</mi> <mrow> <mi>O</mi> <mi>O</mi> <mi>O</mi> </mrow> </msub> <msub> <mi>orV</mi> <mrow> <mi>P</mi> <mi>P</mi> <mi>P</mi> </mrow> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> </mtd> <mtd> <mrow> <mi>x</mi> <mo>=</mo> <mn>2</mn> <mo>,</mo> <mn>4</mn> <mo>,</mo> <mn>6</mn> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <mtable> <mtr> <mtd> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mi>V</mi> <mrow> <mi>V</mi> <mi>M</mi> <mi>x</mi> </mrow> </msub> <mo>=</mo> <mfrac> <mrow> <mn>1</mn> <mo>+</mo> <msub> <mi>K</mi> <mrow> <mi>M</mi> <mi>x</mi> </mrow> </msub> </mrow> <mn>3</mn> </mfrac> <msub> <mi>V</mi> <mrow> <mi>M</mi> <mi>x</mi> </mrow> </msub> <mo>+</mo> <mfrac> <mrow> <mn>1</mn> <mo>-</mo> <msub> <mi>K</mi> <mrow> <mi>M</mi> <mi>x</mi> </mrow> </msub> </mrow> <mn>3</mn> </mfrac> <mrow> <mo>(</mo> <msubsup> <mi>V</mi> <mrow> <mi>S</mi> <mi>x</mi> </mrow> <mo>-</mo> </msubsup> <mo>+</mo> <msubsup> <mi>V</mi> <mrow> <mi>S</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> <mo>+</mo> </msubsup> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>V</mi> <mrow> <mi>V</mi> <mi>S</mi> <mi>x</mi> </mrow> </msub> <mo>=</mo> <mfrac> <mrow> <mn>1</mn> <mo>-</mo> <msub> <mi>K</mi> <mrow> <mi>S</mi> <mi>x</mi> </mrow> </msub> </mrow> <mn>2</mn> </mfrac> <msubsup> <mi>V</mi> <mrow> <mi>S</mi> <mi>x</mi> </mrow> <mo>+</mo> </msubsup> <mo>+</mo> <mfrac> <mrow> <mn>1</mn> <mo>+</mo> <msub> <mi>K</mi> <mrow> <mi>S</mi> <mi>x</mi> </mrow> </msub> </mrow> <mn>2</mn> </mfrac> <msubsup> <mi>V</mi> <mrow> <mi>S</mi> <mi>x</mi> </mrow> <mo>-</mo> </msubsup> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>V</mi> <mrow> <mi>V</mi> <mi>L</mi> <mi>x</mi> </mrow> </msub> <mo>=</mo> <msub> <mi>V</mi> <mrow> <mi>L</mi> <mi>x</mi> </mrow> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>V</mi> <mrow> <mi>V</mi> <mi>Z</mi> </mrow> </msub> <mo>=</mo> <msub> <mi>V</mi> <mrow> <mi>N</mi> <mi>N</mi> <mi>N</mi> </mrow> </msub> <msub> <mi>orV</mi> <mrow> <mi>O</mi> <mi>O</mi> <mi>O</mi> </mrow> </msub> <msub> <mi>orV</mi> <mrow> <mi>P</mi> <mi>P</mi> <mi>P</mi> </mrow> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> </mtd> <mtd> <mrow> <mi>x</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>3</mn> <mo>,</mo> <mn>5</mn> </mrow> </mtd> <mtd> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mi>V</mi> <mrow> <mi>V</mi> <mi>M</mi> <mi>x</mi> </mrow> </msub> <mo>=</mo> <mfrac> <mrow> <mn>1</mn> <mo>+</mo> <msub> <mi>K</mi> <mrow> <mi>M</mi> <mi>x</mi> </mrow> </msub> </mrow> <mn>3</mn> </mfrac> <msub> <mi>V</mi> <mrow> <mi>M</mi> <mi>x</mi> </mrow> </msub> <mo>+</mo> <mfrac> <mrow> <mn>1</mn> <mo>-</mo> <msub> <mi>K</mi> <mrow> <mi>M</mi> <mi>x</mi> </mrow> </msub> </mrow> <mn>3</mn> </mfrac> <mrow> <mo>(</mo> <msubsup> <mi>V</mi> <mrow> <mi>S</mi> <mi>x</mi> </mrow> <mo>+</mo> </msubsup> <mo>+</mo> <msubsup> <mi>V</mi> <mrow> <mi>S</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> <mo>-</mo> </msubsup> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>V</mi> <mrow> <mi>V</mi> <mi>S</mi> <mi>x</mi> </mrow> </msub> <mo>=</mo> <mfrac> <mrow> <mn>1</mn> <mo>-</mo> <msub> <mi>K</mi> <mrow> <mi>S</mi> <mi>x</mi> </mrow> </msub> </mrow> <mn>2</mn> </mfrac> <msubsup> <mi>V</mi> <mrow> <mi>S</mi> <mi>x</mi> </mrow> <mo>+</mo> </msubsup> <mo>+</mo> <mfrac> <mrow> <mn>1</mn> <mo>+</mo> <msub> <mi>K</mi> <mrow> <mi>S</mi> <mi>x</mi> </mrow> </msub> </mrow> <mn>2</mn> </mfrac> <msubsup> <mi>V</mi> <mrow> <mi>S</mi> <mi>x</mi> </mrow> <mo>-</mo> </msubsup> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>V</mi> <mrow> <mi>V</mi> <mi>L</mi> <mi>x</mi> </mrow> </msub> <mo>=</mo> <msub> <mi>V</mi> <mrow> <mi>L</mi> <mi>x</mi> </mrow> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>V</mi> <mrow> <mi>V</mi> <mi>Z</mi> </mrow> </msub> <mo>=</mo> <msub> <mi>V</mi> <mrow> <mi>N</mi> <mi>N</mi> </mrow> </msub> <msub> <mi>orV</mi> <mrow> <mi>O</mi> <mi>O</mi> <mi>O</mi> </mrow> </msub> <msub> <mi>orV</mi> <mrow> <mi>P</mi> <mi>P</mi> <mi>P</mi> </mrow> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> </mtd> <mtd> <mrow> <mi>x</mi> <mo>=</mo> <mn>2</mn> <mo>,</mo> <mn>4</mn> <mo>,</mo> <mn>6</mn> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow>
式中,KMx,KSx∈[-1,1],VMx,VLx分别是三电平空间矢量平面中的基本中矢量、基本大矢量、负组基本小矢量和正组基本小矢量,VNNN,VOOO,VPPP三者为三电平空间矢量平面中的三个零矢量;改进虚拟中矢量将传统虚拟矢量分割成两个部分:第一,基本中矢量部分;第二,相邻正负基本小矢量合成部分;在改进型虚拟矢量中无论引入的调节因子如何变化,都不影响最终合成的虚拟空间矢量的相位和幅值;
步骤b.根据边界条件确定参考电压矢量所在区域,所述的区域包括空间矢量平面大扇区和扇区内的三角形区域;然后选择构成三角形三个顶点所对应的虚拟空间矢量来合成参考电压矢量;根据伏秒平衡原理通过式(3)计算出被选中的三个虚拟空间矢量的作用时间;
<mrow> <msub> <mi>V</mi> <mrow> <mi>r</mi> <mi>e</mi> <mi>f</mi> </mrow> </msub> <mo>=</mo> <mfrac> <msub> <mi>t</mi> <mrow> <mi>V</mi> <mn>1</mn> </mrow> </msub> <msub> <mi>T</mi> <mi>s</mi> </msub> </mfrac> <msub> <mi>V</mi> <mrow> <mi>V</mi> <mn>1</mn> </mrow> </msub> <mo>+</mo> <mfrac> <msub> <mi>t</mi> <mrow> <mi>V</mi> <mn>2</mn> </mrow> </msub> <msub> <mi>T</mi> <mi>s</mi> </msub> </mfrac> <msub> <mi>V</mi> <mrow> <mi>V</mi> <mn>2</mn> </mrow> </msub> <mo>+</mo> <mfrac> <msub> <mi>t</mi> <mrow> <mi>V</mi> <mn>0</mn> </mrow> </msub> <msub> <mi>T</mi> <mi>s</mi> </msub> </mfrac> <msub> <mi>V</mi> <mrow> <mi>V</mi> <mn>0</mn> </mrow> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow>
步骤c.依据基于模型预测的多目标优化中点电位直流不平衡控制策略,实时计算出动态调节因子K,送入步骤b重新调整各矢量的时间分配,消除中点电位直流不平衡。
2.根据权利要求1所述的一种具有中点电位平衡控制的空间矢量调制方法,其特征是:所述的计算价值函数最小的动态调节因子;具体步骤如下:
步骤a.根据直流侧电路简化模型建立中点电位离散数学模型,根据基尔霍夫电流定律流入三电平变频器中点的电流满足如下关系:
<mrow> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mi>i</mi> <mn>1</mn> </msub> <mo>=</mo> <msub> <mi>C</mi> <mn>1</mn> </msub> <mfrac> <mrow> <msub> <mi>du</mi> <mrow> <mi>C</mi> <mn>1</mn> </mrow> </msub> </mrow> <mrow> <mi>d</mi> <mi>t</mi> </mrow> </mfrac> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>i</mi> <mn>2</mn> </msub> <mo>=</mo> <msub> <mi>C</mi> <mn>2</mn> </msub> <mfrac> <mrow> <msub> <mi>du</mi> <mrow> <mi>C</mi> <mn>2</mn> </mrow> </msub> </mrow> <mrow> <mi>d</mi> <mi>t</mi> </mrow> </mfrac> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>i</mi> <mrow> <mi>N</mi> <mi>P</mi> </mrow> </msub> <mo>=</mo> <msub> <mi>i</mi> <mn>1</mn> </msub> <mo>+</mo> <msub> <mi>i</mi> <mn>2</mn> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow>
将式(4)中第1,2式带入第3式中得到其中vNP=uC1-uC2,vNP0为中点电位初始时刻偏差值;通常采样频率远高于工频时,用离散求和的方式来代替积分项,中点电位的离散模型表示为
<mrow> <mtable> <mtr> <mtd> <mrow> <msub> <mi>v</mi> <mrow> <mi>N</mi> <mi>P</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>1</mn> <mi>C</mi> </mfrac> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <msub> <mi>i</mi> <mrow> <mi>N</mi> <mi>P</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <msub> <mi>T</mi> <mi>S</mi> </msub> <mo>+</mo> <msub> <mi>v</mi> <mrow> <mi>N</mi> <mi>P</mi> <mn>0</mn> </mrow> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>=</mo> <mfrac> <mn>1</mn> <mi>C</mi> </mfrac> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <munder> <mo>&amp;Sigma;</mo> <mrow> <mi>m</mi> <mo>=</mo> <mi>a</mi> <mo>,</mo> <mi>b</mi> <mo>,</mo> <mi>c</mi> </mrow> </munder> <msub> <mi>i</mi> <mi>m</mi> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>&amp;CenterDot;</mo> <msub> <mi>t</mi> <mrow> <mi>m</mi> <mn>0</mn> </mrow> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>v</mi> <mrow> <mi>N</mi> <mi>P</mi> <mn>0</mn> </mrow> </msub> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> </mrow>
其中,im(k)为第k个采样时刻m相所对应的相电流值,tm0(k)为在第k个采样周期内m相电流对应流入变频器中点的时间;获得中点电位的变化率,以式(6)表示;
<mrow> <mtable> <mtr> <mtd> <mrow> <msub> <mi>&amp;Delta;v</mi> <mrow> <mi>N</mi> <mi>P</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>v</mi> <mrow> <mi>N</mi> <mi>P</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>v</mi> <mrow> <mi>N</mi> <mi>P</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>=</mo> <mfrac> <mn>1</mn> <mi>C</mi> </mfrac> <munder> <mo>&amp;Sigma;</mo> <mrow> <mi>m</mi> <mo>=</mo> <mi>a</mi> <mo>,</mo> <mi>b</mi> <mo>,</mo> <mi>c</mi> </mrow> </munder> <msub> <mi>i</mi> <mi>m</mi> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>&amp;CenterDot;</mo> <msub> <mi>t</mi> <mrow> <mi>m</mi> <mn>0</mn> </mrow> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow>
步骤b.建立目标函数如下
<mrow> <mi>J</mi> <mo>=</mo> <msup> <mrow> <mo>&amp;lsqb;</mo> <msubsup> <mi>v</mi> <mrow> <mi>N</mi> <mi>P</mi> </mrow> <mo>*</mo> </msubsup> <mo>-</mo> <msub> <mi>v</mi> <mrow> <mi>N</mi> <mi>P</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>&amp;rsqb;</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <mi>&amp;lambda;</mi> <mo>&amp;CenterDot;</mo> <mo>|</mo> <msub> <mi>v</mi> <mrow> <mi>N</mi> <mi>P</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>|</mo> <mo>&amp;CenterDot;</mo> <msub> <mi>&amp;Delta;v</mi> <mrow> <mi>N</mi> <mi>P</mi> </mrow> </msub> <msup> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow>
其中,λ为权重因子,为中点电位理想参考值通常为零;目标函数由两部分构成:a)为第一部分,当第一部分满足取值最小时可以保证在下一个采样时刻中点电位以最快的速度趋近于给定参考值b)λ·vNP(k)·ΔvNP(k)2为第二部分,当第二部分取最小值是可以保证直流不平衡恢复过程最为平稳;通过选取不同的调节因子可以在速度和品质之间获得平衡;
步骤c.将当前时刻的中点电位值中点电位变化率送入目标函数中,解得使目标函数获得最小值时所对应的K值的最优解为
<mrow> <mi>K</mi> <mo>=</mo> <mo>-</mo> <mfrac> <mn>1</mn> <mrow> <mo>(</mo> <mi>&amp;lambda;</mi> <mo>+</mo> <mn>1</mn> <mo>/</mo> <mo>|</mo> <mrow> <msub> <mi>v</mi> <mrow> <mi>N</mi> <mi>P</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> </mrow> <mo>|</mo> <mo>)</mo> <munder> <mo>&amp;Sigma;</mo> <mrow> <mi>m</mi> <mo>=</mo> <mi>a</mi> <mo>,</mo> <mi>b</mi> <mo>,</mo> <mi>c</mi> </mrow> </munder> <msub> <mi>i</mi> <mi>m</mi> </msub> <mo>(</mo> <mi>k</mi> <mo>)</mo> <mo>&amp;CenterDot;</mo> <mfrac> <mrow> <msub> <mi>t</mi> <mrow> <mi>m</mi> <mn>0</mn> <mo>_</mo> <mi>max</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> </mrow> <msub> <mi>T</mi> <mi>S</mi> </msub> </mfrac> </mrow> </mfrac> <mo>&amp;CenterDot;</mo> <mi>sgn</mi> <mrow> <mo>(</mo> <msub> <mi>v</mi> <mrow> <mi>N</mi> <mi>P</mi> </mrow> </msub> <mo>(</mo> <mi>k</mi> <mo>)</mo> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>8</mn> <mo>)</mo> </mrow> </mrow>
其中,tm0_max(k)为K取边界值时所对应的作用时间,sgn(*)为符号函数。
CN201510701913.0A 2015-10-26 2015-10-26 一种具有中点电位平衡控制的空间矢量调制方法 Expired - Fee Related CN105226981B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510701913.0A CN105226981B (zh) 2015-10-26 2015-10-26 一种具有中点电位平衡控制的空间矢量调制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510701913.0A CN105226981B (zh) 2015-10-26 2015-10-26 一种具有中点电位平衡控制的空间矢量调制方法

Publications (2)

Publication Number Publication Date
CN105226981A CN105226981A (zh) 2016-01-06
CN105226981B true CN105226981B (zh) 2017-11-24

Family

ID=54995747

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510701913.0A Expired - Fee Related CN105226981B (zh) 2015-10-26 2015-10-26 一种具有中点电位平衡控制的空间矢量调制方法

Country Status (1)

Country Link
CN (1) CN105226981B (zh)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106385192B (zh) * 2016-10-31 2020-03-31 东北大学秦皇岛分校 三电平四桥臂逆变器控制方法
CN106655854B (zh) * 2016-11-12 2019-02-26 国网江西省电力公司电力科学研究院 一种分区域调制的npc逆变器中点电位平衡控制方法
CN107104604B (zh) * 2017-04-14 2019-05-14 华南理工大学 一种三电平并网逆变器模型预测直接功率控制方法
CN107872167B (zh) * 2017-11-21 2020-06-09 国家电网公司 一种电压型三电平中点钳位变流器虚拟空间矢量调制方法
CN108390582B (zh) * 2018-03-21 2019-02-22 四川大学 一种单相三电平变流器空间矢量脉宽调制优化方法
CN108683351B (zh) * 2018-06-08 2019-10-18 山东大学 一种z源三电平逆变器的混合调制方法、控制器及系统
CN109600065B (zh) * 2018-12-14 2020-01-10 合肥工业大学 一种三电平变流器不连续脉宽调制方法
CN113098353B (zh) * 2021-05-24 2022-09-06 广东美芝制冷设备有限公司 基于变频器的脉宽调制方法、装置、变频器及存储介质
CN113162451B (zh) * 2021-05-27 2022-11-11 华北电力大学(保定) 一种多电平逆变器的中点电位平衡的控制方法及装置

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103138591B (zh) * 2011-11-24 2015-09-23 深圳古瑞瓦特新能源股份有限公司 一种三电平变换中点平衡控制方法、装置及三电平变换器
JP6225418B2 (ja) * 2012-04-19 2017-11-08 株式会社明電舎 3レベルインバータの制御装置
US9214874B2 (en) * 2012-07-31 2015-12-15 Yashomani Y. Kolhatkar Intelligent level transition systems and methods for transformerless uninterruptible power supply

Also Published As

Publication number Publication date
CN105226981A (zh) 2016-01-06

Similar Documents

Publication Publication Date Title
CN105226981B (zh) 一种具有中点电位平衡控制的空间矢量调制方法
CN103701350B (zh) 低频工况下模块化多电平变流器电容电压波动抑制方法
CN109446653B (zh) 混合型mmc状态空间模型的混合型mmc小信号模型的建模方法
CN103746585B (zh) 基于混合调制的多电平逆变器中点电压平衡控制方法
CN103746584B (zh) 基于载波偏置的多电平逆变器中点电压平衡控制方法
CN104836464B (zh) 一种vienna整流器直流侧中点电位平衡控制装置及方法
CN105610339B (zh) 一种虚拟空间矢量的三电平调制方法
CN111371337B (zh) 二极管箝位型三电平逆变器的中性点电位平衡控制方法
CN107317500A (zh) 一种中点箝位四电平变换器的母线电容电压平衡控制方法
CN105978374A (zh) 三电平逆变器中点电压平衡和共模电压抑制的方法
CN107196536A (zh) 一种具有中点平衡和共模电压抑制能力的三电平svpwm方法
CN104393779B (zh) 一种基于载波层叠调制的模块化多电平变流器控制方法
CN106357135A (zh) 一种vienna整流器中点电位交直流分量平衡控制方法
CN104601028B (zh) 参数在线整定的中点电压控制系统及方法
CN107465359A (zh) 模块化多电平换流器的环流抑制方法及装置
CN106533230A (zh) 基于简化平衡因子的三电平虚拟空间矢量均压调制方法
CN105226982A (zh) 一种基于中点电流的三电平npc逆变器中点电位均衡控制方法
CN103441502B (zh) 并联单相h桥级联型有源电力滤波器控制装置及其方法
CN106169879A (zh) 修正注入零序分量的vienna整流器调制方法、控制器及系统
CN105720587A (zh) 一种级联型星形连接statcom直流母线相间均压控制方法
CN108471248B (zh) 基于α’-β’坐标的星形级联逆变器空间矢量调制方法
CN105591548A (zh) 基于多端口高频变压器的自平衡式电力电子变压器
CN104810854B (zh) 串联型微电网及其微源之间的功率协调控制方法
CN204290329U (zh) 一种svg直流侧电容中点电位平衡装置
CN102694385A (zh) D-statcom的线电流不对称补偿的相电流平衡限幅方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
CB02 Change of applicant information

Address after: 221116 Xuzhou Copper Mountain University Road, China University of Mining and Technology Research Institute,

Applicant after: China University of Mining & Technology

Address before: 221116 Research Institute, China University of Mining and Technology, Xuzhou University, Jiangsu, China,

Applicant before: China University of Mining & Technology

CB02 Change of applicant information
CB03 Change of inventor or designer information
CB03 Change of inventor or designer information

Inventor after: Tang Die

Inventor after: Zhang Chuanjin

Inventor after: Gong Lifeng

Inventor after: Wang Ke

Inventor after: Shen Jia

Inventor after: Li Peng

Inventor before: Tang Die

Inventor before: Zhang Chuanjin

Inventor before: Wang Ke

Inventor before: Han Dong

Inventor before: Yu Bo

TA01 Transfer of patent application right

Effective date of registration: 20171023

Address after: 221116 Research Institute of China University of Mining and Technology,, Jiangsu

Applicant after: China University of Mining & Technology

Applicant after: SUZHOU POWER SUPPLY COMPANY OF STATE GRID JIANGSU ELECTRIC POWER COMPANY

Address before: 221116 Xuzhou Copper Mountain University Road, China University of Mining and Technology Research Institute,

Applicant before: China University of Mining & Technology

TA01 Transfer of patent application right
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20171124

Termination date: 20201026