CN114676573A - 三电平光伏逆变器开路故障统一数学建模方法 - Google Patents
三电平光伏逆变器开路故障统一数学建模方法 Download PDFInfo
- Publication number
- CN114676573A CN114676573A CN202210313222.3A CN202210313222A CN114676573A CN 114676573 A CN114676573 A CN 114676573A CN 202210313222 A CN202210313222 A CN 202210313222A CN 114676573 A CN114676573 A CN 114676573A
- Authority
- CN
- China
- Prior art keywords
- photovoltaic inverter
- level photovoltaic
- open
- circuit fault
- representing
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 23
- 238000013178 mathematical model Methods 0.000 claims abstract description 31
- 239000011159 matrix material Substances 0.000 claims description 27
- 239000003990 capacitor Substances 0.000 claims description 14
- 230000007935 neutral effect Effects 0.000 claims description 6
- 230000027311 M phase Effects 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- 238000004088 simulation Methods 0.000 abstract description 15
- 238000010586 diagram Methods 0.000 description 7
- 229910052760 oxygen Inorganic materials 0.000 description 5
- 229910052757 nitrogen Inorganic materials 0.000 description 3
- 229910052698 phosphorus Inorganic materials 0.000 description 3
- 238000010248 power generation Methods 0.000 description 3
- 238000003745 diagnosis Methods 0.000 description 2
- 238000001338 self-assembly Methods 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/18—Network design, e.g. design based on topological or interconnect aspects of utility systems, piping, heating ventilation air conditioning [HVAC] or cabling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- Business, Economics & Management (AREA)
- Economics (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Health & Medical Sciences (AREA)
- Strategic Management (AREA)
- Primary Health Care (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Water Supply & Treatment (AREA)
- Computer Networks & Wireless Communication (AREA)
- Computational Mathematics (AREA)
- Public Health (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Marketing (AREA)
- Human Resources & Organizations (AREA)
- General Health & Medical Sciences (AREA)
- Inverter Devices (AREA)
Abstract
本发明公开了一种三电平光伏逆变器开路故障的统一数学建模方法,其步骤包括:1根据三电平光伏逆变器处于最大功率点跟踪运行时直流侧电压波动的特点,得出三电平光伏逆变器直流侧电压的数学表达式;2根据三电平光伏逆变器的主电路功率开关器件不同工作状况和电流路径,得出对应故障状况下的开关函数表达式,最后得出统一的开路故障开关函数表达式和三电平光伏逆变器输出相电压的表达式。本发明能够通过开关函数建模方式快速地建立正确的三电平光伏逆变器开路故障统一数学模型,从而解决现有三电平光伏逆变器故障仿真电路建模方法存在的模块数量多、计算量大、不同类故障切换仿真灵活性差等问题,减少了建模的工作量及其仿真工作时间。
Description
技术领域
本发明涉及一种三电平光伏逆变器开路故障的统一数学建模方法,属于光伏运维故障诊断技术领域。
背景技术
光伏逆变器是分布式光伏发电系统中的一个关键电能变换设备,分布式光伏发电地域分布分散、运行环境恶劣、运行工况多变,易导致光伏逆变器故障,对光伏发电可靠运行及系统安全造成危害。光伏逆变器故障一般分为结构性故障和参数性故障。结构性故障为功率器件的开路或短路故障,引起系统结构的变化。参数性故障为元器件的参数漂移,如电容、电感、变压器等元件特性的改变。参数性故障降低了系统的工作特性和安全裕度,如果得不到及时诊断与处理,往往会引发结构性故障,导致光伏逆变器退网。据统计,逆变器开关功率器件开路故障发生概率比较高。
目前文献研究的开路故障大多数是基于两电平光伏逆变器,然而如今分布式光伏电站中的光伏组串逆变器多数采用的是三电平光伏逆变器具有纹波和谐波含量少、开关损耗降低等优点。由于三电平光伏逆变器其内部开关器件数量增多,使其发生故障时出现九种情况且多达78种开路故障类型,导致三电平光伏逆变器故障的复杂性增加。现有的三电平光伏逆变器开关故障模型多数采用MATLAB/Simulink仿真软件进行搭建,但Simulink自带的仿真模块不适用于三电平光伏逆变器开关故障仿真。为了解决这问题,多数采用自搭模块,但自搭模块存在花费时间长且不同类故障切换仿真灵活性差等问题。因此可见,建立实际三电平光伏逆变器功率开关管开路故障的数学模型是有必要的。
发明内容
本发明是为了解决上述现有技术存在的不足之处,提出一种三电平光伏逆变器开路故障的统一数学建模方法,以期能通过故障机理建模方式建立较为简便的三电平光伏逆变器开关故障统一数学模型,从而能够更好地识别判断三电平光伏逆变器开路故障,并提升开路故障诊断的快速性,以解决现有三电平光伏逆变器故障仿真电路建模方法存在的模块数量多、计算量大、不同类故障切换仿真灵活性差等问题。
本发明为达到上述发明目的,采用如下技术方案:
本发明一种三电平光伏逆变器开路故障的统一数学建模方法的特点在于,包括以下步骤:
步骤1、根据三电平光伏逆变器处于最大功率点跟踪运行时直流侧电压波动的特点,得出三电平光伏逆变器直流侧电压的数学表达式;
步骤2、根据三电平光伏逆变器的主电路功率开关器件正常工作状态下的换流情况,得出正常时的交流侧输出相电压的数学表达式;
步骤3、根据故障物理模型,得到发生开路故障时的拓扑结构以及电流路径;
步骤4、构建不同开路故障下的开关函数所表达的开路故障统一的开关函数数学模型表达式,从而得到三电平光伏逆变器输出相电压的表达式。
本发明所述的三电平光伏逆变器开路故障的统一数学建模方法的特点也在于,所述步骤1中,是利用式(1)得到三电平光伏逆变器直流侧电压的数学表达式:
式(1)中:Udc代表三电平光伏逆变器的直流侧电压;Umpp代表三电平光伏逆变器处于最大功率点跟踪运行时的直流侧稳态电压;Δu代表三电平光伏逆变器处于最大功率点跟踪运行时的直流侧电压波动;Tm代表MPPT的控制采样周期;Cdc代表直流侧电容C1或C2的容值,且C1=C2=Cdc;Iin代表三电平光伏逆变器的直流侧输入电流;Sk′代表k相的三电平光伏逆变器统一开关函数;Ik代表k相的三电平光伏逆变器输出电流,k=A,B,C。
所述步骤2包括:
步骤2.1、利用式(2)得到正常时的三电平光伏逆变器开关函数的数学表达式:
S′k(0)=Sk=Sk1-Sk2 (2)
式(2)中:S′k(0)代表k相正常工作状态下的开关函数;Sk代表三电平光伏逆变器理论等效的开关函数;Sk1和Sk2是用于表示Sk的两个自定义开关量,且Sk=Sk1-Sk2;k=A,B,C;
步骤2.2、利用式(3)得到三电平光伏逆变器的输出相电压的数学表达式:
式(3)中,UkO代表三电平光伏逆变器k相的输出相电压;UkN代表三电平光伏逆变器k相输出端对直流侧电容中点N的电压;UNO代表直流侧电容中点N对电网中点O的电压;S′k(0)代表k相正常工作状态下的开关函数;S′m(0)代表m相正常工作状态下的开关函数;S′n(0)代表n相正常工作状态下的开关函数,k=A,B,C;m=A,B,C且m≠k;n=A,B,C且n≠k≠m。
所述步骤3中的开路故障包括:单个开关管发生开路故障以及两个开关管同时发生开路故障。
所述步骤4中是利用式(3)构建开路故障统一的开关函数数学模型表达式,从而利用式(4)得到三电平光伏逆变器输出相电压UkO的表达式:
式(3)和(4)中:δk表示三电平光伏逆变器输出电流的方向为流入或流出的逻辑量;δk=1,代表电流方向是从三电平光伏逆变器流向电网;δk=0,代表电流方向是从电网流向三电平光伏逆变器;D是代表三电平光伏逆变器工作状态的3×3单位矩阵;E是代表Qk1发生开路故障的3×3矩阵;F是代表Qk2发生开路故障的3×3矩阵;G是代表Qk3发生开路故障的3×3矩阵;H是代表Qk4发生开路故障的3×3矩阵;矩阵E、F、G、H中,当发生开路故障时,其对应的矩阵对角元素为1;未发生开路故障时,其对应的矩阵对角元素为0。
与现有技术相比,本发明的有益效果在于:
1、本发明以LCL型并网三电平光伏逆变器作为对象,以三电平光伏逆变器正常工作状态下的数学模型为基础,根据三电平光伏逆变器处于最大功率点跟踪运行时直流侧电压波动的特点,得出三电平光伏逆变器直流侧电压的数学表达式,从故障机理出发,根据三电平光伏逆变器开关管开关故障时的拓扑图以及实际电流的电流路径,推导出不同开路故障的三电平光伏逆变器开关故障数学模型以及统一的数学模型,使得该模型能够快速地描述正确的开关故障时的三电平光伏逆变器输出电流特性。
2、相比于现有的几种三电平光伏逆变器开关函数建模方法,本发明为三电平光伏逆变器开路故障统一数学建模,该模型直观地表述出发生开关管开路故障时对三电平光伏逆变器输出参数的影响且通过一个统一的数学模型得出不同开路故障下的三电平光伏逆变器输出参数,减小了计算量,同时能保证建模精度;本发明仅依靠三电平光伏逆变器理论等效的开关函数Sk、三电平光伏逆变器实际输出电流的方向δk、直流侧电压Udc三种逆变器参数来建立三电平光伏逆变器开关故障统一模型,减少了建模的工作量与工作时间、不同类故障切换仿真灵活性强。
附图说明
图1为本发明中三电平光伏逆变器开关故障统一建模方法的流程图;
图2为本发明中LCL型并网三电平光伏逆变器物理模型图(省略了输入级boost电路);
图3为本发明中三电平光伏逆变器正常工作状态下k相的电流路径(P状态)图;
图4为本发明中三电平光伏逆变器正常工作状态下k相的电流路径(O状态)图;
图5为本发明中三电平光伏逆变器正常工作状态下k相的电流路径(N状态)图;
图6a为本发明中Qk1发生开路故障时k相的电流路径(P状态)图;
图6b为本发明中QA1发生开路故障物理模型仿真结果与数学模型结果对比图;
图7a为本发明中Qk2发生开路故障时k相的电流路径(P和O状态)图;
图7b为本发明中QA2发生开路故障物理模型仿真结果与数学模型结果对比图;
图8为本发明中Qk3发生开路故障时k相的电流路径(O和N状态)图;
图9为本发明中Qk4发生开路故障时k相的电流路径(N状态)图;
图10为本发明中QA1和QB1发生开路故障物理模型仿真结果与数学模型结果对比图。
具体实施方式
本实施例中,一种三电平光伏逆变器开路故障统一建模方法是以LCL型并网三电平光伏逆变器作为对象,以三电平光伏逆变器正常工作状态下的数学模型为基础,根据三电平光伏逆变器处于最大功率点跟踪运行时直流侧电压波动的特点,得出三电平光伏逆变器直流侧电压的数学表达式,从故障机理出发,根据三电平光伏逆变器开关管开关故障时的拓扑图以及实际电流的电流路径,推导出不同开路故障的三电平光伏逆变器开关故障数学模型以及统一的数学模型,使得该模型能够快速地描述正确的开关故障时的三电平光伏逆变器输出电流特性。具体的说,如图1所示,建模方法包括以下步骤:
步骤1、根据三电平光伏逆变器处于最大功率点跟踪运行时直流侧电压波动的特点,得出三电平光伏逆变器直流侧电压的数学表达式:
LCL型并网三电平光伏逆变器物理模型图如图2所示,根据基尔霍夫电流定理:
式(1)和式(2)中,Iin代表三电平光伏逆变器直流侧输入电流;Idc1代表三电平光伏逆变器正极侧输入电流(Iin经过电容C1分流后);Idc2代表三电平光伏逆变器负极侧输入电流(Iin经过电容C2分流后);UC1代表输入侧的电容C1两端电压;UC2代表输入侧的电容C2两端电压;IC1代表流经输入侧的电容C1的电流;IC2代表流经输入侧的电容C2的电流。
设三电平光伏逆变器直流侧电容C1、C2的容值为Cdc,则C1=C2=Cdc,则根据式(1)和式(2)可以得出:
根据Idc1-Idc2=S′AIA+S′BIB+S′CIC可以将式(3)变换成:
式(4)中:S′k代表k相的三电平光伏逆变器统一开关函数;Ik代表k相的三电平光伏逆变器输出电流,k=A,B,C。
三电平光伏逆变器在处于最大功率点跟踪运行时,会导致三电平光伏逆变器直流侧电压存在波动,最后根据式(4)可以得知,电平光伏逆变器处于最大功率点跟踪运行时的直流侧电压波动Δu表达式可以写成:
式(5)中,Δu代表三电平光伏逆变器处于最大功率点跟踪运行时的直流侧电压波动;Tm代表MPPT的控制采样周期。
故此三电平光伏逆变器直流侧电压表达式:
式(6)中:Udc代表三电平光伏逆变器直流侧电压;Umpp代表三电平光伏逆变器处于最大功率点跟踪运行时的直流侧稳态电压。
步骤2、根据三电平光伏逆变器的主电路功率开关器件正常工作状态下的换流情况,得出正常时的交流侧输出相电压的数学表达式:
定义开关量Sk1和Sk2来表示理论等效的开关函数Sk,即Sk=Sk1-Sk2。三电平光伏逆变器工作状态有三个状态,分别为P、O和N。开关函数Sk可以用于表示三电平光伏逆变器的三种工作状态,描述如下:正常的工作状态下,P状态(Sk=1)下,此时所定义的开关量满足Sk1=1、Sk2=0;O状态(Sk=0)下,此时所定义的开关量满足Sk1=0、Sk2=0;N状态(Sk=-1)下,此时所定义的开关量满足Sk1=0、Sk2=1。P O N三个状态下的电流路径如图3、图4和图5所示。三电平光伏逆变器k相输出端对直流侧电容中点N的电压UkN可以用表达式来表述。
定义三电平光伏逆变器实际等效的开路故障统一的开关函数为Sk′,则当三电平光伏逆变器正常工作状态下,有S′k=Sk。规定三电平光伏逆变器实际输出电流方向从三电平光伏逆变器流向电网为正,从电网流向三电平光伏逆变器为负。定义其逻辑变量为δk(k=A,B,C下同)。δk代表的是k相电流的电流方向。δk=1,代表电流方向为正(从三电平光伏逆变器流向电网);δk=0,代表电流方向为负(从电网流向三电平光伏逆变器)。
当三电平光伏逆变器处于正常工作状态下,三电平光伏逆变器输出端对直流侧电容中点的电压UkN的真值表如表1所示:
表1三电平光伏逆变器正常工作状态下,电压UkN真值表
S<sub>k</sub> | S<sub>k</sub>′ | S<sub>k1</sub> | S<sub>k2</sub> | δ<sub>k</sub> | U<sub>kN</sub> |
1 | 1 | 1 | 0 | 1,0 | U<sub>dc</sub>/2 |
0 | 0 | 0 | 0 | 1,0 | 0 |
-1 | -1 | 0 | 1 | 1,0 | -U<sub>dc</sub>/2 |
根据表1可以得出,三电平光伏逆变器k点对N点电压UkN的数学描述为:
根据式(7)得知开关函数S′k可以写成:
S′k(0)=Sk=Sk1-Sk2 (8)
此时,三电平光伏逆变器输出相电压数学表达式:
其中,UkO代表三电平光伏逆变器k相的输出相电压;UkN代表三电平光伏逆变器k相输出端对直流侧电容中点N的电压;UNO代表直流侧电容中点N对电网中点O的电压;S′k(0)代表k相正常工作状态下的开关函数;S′m(0)代表m相正常工作状态下的开关函数;S′n(0)代表n相正常工作状态下的开关函数,k=A,B,C;m=A,B,C且m≠k;n=A,B,C且n≠k≠m。
步骤3、根据故障物理模型,得到发生开路故障时的拓扑结构以及电流路径:
当开关管发生开路故障时,使得三电平光伏逆变器工作的拓扑以及输出电流路径不同,此时三电平光伏逆变器输出相电压不仅仅与开关管开关状态有关,也与输出电流的电流方向有关。一般在逆变器实际运行工作中,很少出现3只或3只以上开关管同时故障的情况,所以下面只考虑单管以及双管发生开路故障。
当开关管Qk1发生开路故障时,P状态受到影响,而O和N状态均未受到影响。在P状态下,其电流路径如图6所示:
根据图6a和图6b可知,此时电流路径与正常工作状态下不同,由于Qk1发生了开路故障,使得当电流Ik>0时,电流Ik通过二极管Dk5以及开关管Qk2流向电网,不再是正常工作状态下Qk1和Qk2的电流路径。此时三电平光伏逆变器工作状态从P状态变成O状态。此时三电平光伏逆变器k点对N点的电压UkN与电流方向和工作状态之间的关系如表2所示:
表2 Qk1发生开路故障情况下,电压UkN真值表
当Qk2发生开路故障,则对三电平光伏逆变器状态P和O产生影响,对N状态未造成影响。
根据图7a和图7b可知Qk2发生了开路故障,在P状态下,当电流Ik>0时,电流Ik通过二极管Dk3以及二极管Dk4流向电网,不再是正常工作状态下Qk1和Qk2的电流路径。此时三电平光伏逆变器工作状态从P状态变成N状态。在O状态下,当电流Ik>0时,电流Ik通过二极管Dk3以及二极管Dk4流向电网,不再是正常工作状态下Dk5和Qk2的电流路径。此时三电平光伏逆变器工作状态从O状态变成N状态。分析得知Qk2发生开路故障时,三电平光伏逆变器k点对N点的电压UkN与电流方向和工作状态之间的关系如表3所示:
表3 Qk2发生开路故障情况下,电压UkN真值表
当Qk3发生开路故障,则对三电平光伏逆变器状态O和N产生影响,对P状态未造成影响。
根据图8可知Qk3发生了开路故障,在O状态下,当电流Ik<0时,电流Ik通过二极管Dk1以及二极管Dk2流向逆变器,不再是正常工作状态下Qk3和Dk6的电流路径。此时三电平光伏逆变器工作状态从O状态变成P状态。在N状态下,当电流Ik<0时,电流Ik通过二极管Dk1以及二极管Dk2流向逆变器,不再是正常工作状态下Qk3和Qk4的电流路径。此时三电平光伏逆变器工作状态从N状态变成P状态。分析得知出Qk3发生开路故障时,三电平光伏逆变器k点对N点的电压UkN与电流方向和工作状态之间的关系如表4所示:
表4 Qk3发生开路故障情况下,电压UkN真值表
当Qk4发生开路故障,则对三电平光伏逆变器状态N产生影响,对P和O状态均未造成影响。
根据图9可知,由于Qk4发生了开路故障,使得当电流Ik<0时,电流Ik通过二极管Qk3以及开关管Dk6流向三电平光伏逆变器,不再是正常工作状态下Qk3和Qk4的电流路径。分析得知当Qk4发生开路故障时,三电平光伏逆变器k点对N点的电压UkN与电流方向和工作状态之间的关系如表5所示:
表5 Qk4发生开路故障情况下端电压UkN真值表
发生双管开路故障时,可以看作为分别两个单管开路故障的组合,则开关函数Sk′的变化也可以看作两个单管开关管分别开路故障的开关函数之和。
步骤4、构建不同开路故障下的开关函数所表达的开路故障统一的开关函数数学模型表达式,从而得到三电平光伏逆变器输出相电压的表达式:
根据表2可以得知,当Qk1发生开路故障时,会导致其开关函数发生变化,开路故障后的开关函数Sk′满足式(10):
根据表3可以得知,当Qk2发生开路故障时,会导致其开关函数发生变化,开路故障后的开关函数Sk′满足式(11):
根据表4可以得知,当Qk3发生开路故障时,会导致其开关函数发生变化,开路故障后的开关函数Sk′满足式(12):
根据表5可以得出,当Qk4发生开路故障时,会导致其开关函数发生变化,开路故障后的开关函数Sk′满足式(13):
则根据式(8)、式(10)、式(11)、式(12)和式(13)可以得出,开路故障统一的开关函数Sk′的表达式为:
式(14)中,Sk=Sk1-Sk2;D是代表三电平光伏逆变器工作状态的3×3单位矩阵;E是代表Qk1发生开路故障的3×3矩阵;F是代表Qk2发生开路故障的3×3矩阵;G是代表Qk3发生开路故障的3×3矩阵;H是代表Qk4发生开路故障的3×3矩阵;矩阵EFGH中,当发生开路故障时,其对应的矩阵对角元素为1;未发生开路故障时,其对应的矩阵对角元素为0。
综上所述:三电平光伏逆变器k相输出端对直流侧电容中点N的电压UkN表达式可写成:
根据三相系统对称可知:
根据式(15)和式(16)可以得知三电平光伏逆变器输出相电压UkO表达式:
最后根据基尔霍夫定律,得出LCL型并网三电平光伏逆变器对应的电气参数表达式:
其中,Igk代表k相并网电流;Uck为k相LCL滤波器中的滤波电容两端电压;Ek为电网电压;UkO为三电平光伏逆变器k相输出相电压;Cf为LCL滤波器中的滤波电容容值。
为验证数学模型准确性,首先在Simulink中搭建三电平光伏逆变器开关故障物理模型,再使用上述所得出来的三电平光伏逆变器开关故障统一数学模型搭建数学模型。最后根据物理模型和数学模型的输出电流对比,便能得知开路故障数学模型的正确性。具体如下:
(1)设QA1发生开路故障,分别在物理模型和数学模型设置在0.2s时发生单管开路故障;三电平光伏逆变器在0.1s之前还处于最大功率跟踪的阶段,且QA1发生开路故障物理模型仿真结果与数学模型结果对比图如图6b所示,两者故障前与故障后的电流基本一致。
(2)设QA2发生开路故障,分别在物理模型和数学模型设置在0.2s时发生单管开路故障;三电平光伏逆变器在0.1s之前还处于最大功率跟踪的阶段,且QA2发生开路故障物理模型仿真结果与数学模型结果对比图如图7b所示,两者故障前与故障后的电流基本一致。
(3)设QA1和QB1发生开路故障,分别在物理模型和数学模型设置在0.2s时发生双管开路故障;三电平光伏逆变器在0.1s之前还处于最大功率跟踪的阶段,且QA1和QB1发生开路故障物理模型仿真结果与数学模型结果对比图如图10所示,两者故障前与故障后的电流基本一致。
图6b、图7b和图10的结果证明,本发明提出的三电平光伏逆变器开关故障统一数学模型搭建方法得出的三电平光伏逆变器开路故障数学模型具有较高的准确性等特点。
Claims (5)
1.一种三电平光伏逆变器开路故障的统一数学建模方法,其特征在于,包括以下步骤:
步骤1、根据三电平光伏逆变器处于最大功率点跟踪运行时直流侧电压波动的特点,得出三电平光伏逆变器直流侧电压的数学表达式;
步骤2、根据三电平光伏逆变器的主电路功率开关器件正常工作状态下的换流情况,得出正常时的交流侧输出相电压的数学表达式;
步骤3、根据故障物理模型,得到发生开路故障时的拓扑结构以及电流路径;
步骤4、构建不同开路故障下的开关函数所表达的开路故障统一的开关函数数学模型表达式,从而得到三电平光伏逆变器输出相电压的表达式。
3.根据权利要求2的三电平光伏逆变器开路故障的统一数学建模方法,其特征是,所述步骤2包括:
步骤2.1、利用式(2)得到正常时的三电平光伏逆变器开关函数的数学表达式:
S′k(0)=Sk=Sk1-Sk2 (2)
式(2)中:S′k(0)代表k相正常工作状态下的开关函数;Sk代表三电平光伏逆变器理论等效的开关函数;Sk1和Sk2是用于表示Sk的两个自定义开关量,且Sk=Sk1-Sk2;k=A,B,C;
步骤2.2、利用式(3)得到三电平光伏逆变器的输出相电压的数学表达式:
式(3)中,UkO代表三电平光伏逆变器k相的输出相电压;UkN代表三电平光伏逆变器k相输出端对直流侧电容中点N的电压;UNO代表直流侧电容中点N对电网中点O的电压;S′k(0)代表k相正常工作状态下的开关函数;S′m(0)代表m相正常工作状态下的开关函数;S′n(0)代表n相正常工作状态下的开关函数,k=A,B,C;m=A,B,C且m≠k;n=A,B,C且n≠k≠m。
4.根据权利要求3的三电平光伏逆变器开路故障的统一数学建模方法,其特征是,所述步骤3中的开路故障包括:单个开关管发生开路故障以及两个开关管同时发生开路故障。
5.根据权利要求4的三电平光伏逆变器开路故障统一数学建模方法,其特征是,所述步骤4中是利用式(3)构建开路故障统一的开关函数数学模型表达式,从而利用式(4)得到三电平光伏逆变器输出相电压UkO的表达式:
式(3)和(4)中:δk表示三电平光伏逆变器输出电流的方向为流入或流出的逻辑量;δk=1,代表电流方向是从三电平光伏逆变器流向电网;δk=0,代表电流方向是从电网流向三电平光伏逆变器;D是代表三电平光伏逆变器工作状态的3×3单位矩阵;E是代表Qk1发生开路故障的3×3矩阵;F是代表Qk2发生开路故障的3×3矩阵;G是代表Qk3发生开路故障的3×3矩阵;H是代表Qk4发生开路故障的3×3矩阵;矩阵E、F、G、H中,当发生开路故障时,其对应的矩阵对角元素为1;未发生开路故障时,其对应的矩阵对角元素为0。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210313222.3A CN114676573A (zh) | 2022-03-28 | 2022-03-28 | 三电平光伏逆变器开路故障统一数学建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210313222.3A CN114676573A (zh) | 2022-03-28 | 2022-03-28 | 三电平光伏逆变器开路故障统一数学建模方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN114676573A true CN114676573A (zh) | 2022-06-28 |
Family
ID=82076987
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210313222.3A Pending CN114676573A (zh) | 2022-03-28 | 2022-03-28 | 三电平光伏逆变器开路故障统一数学建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114676573A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115267606A (zh) * | 2022-07-26 | 2022-11-01 | 合肥工业大学 | Anpc型三电平并网逆变器中igbt的开路故障诊断方法 |
-
2022
- 2022-03-28 CN CN202210313222.3A patent/CN114676573A/zh active Pending
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115267606A (zh) * | 2022-07-26 | 2022-11-01 | 合肥工业大学 | Anpc型三电平并网逆变器中igbt的开路故障诊断方法 |
CN115267606B (zh) * | 2022-07-26 | 2024-04-05 | 合肥工业大学 | Anpc型三电平并网逆变器中igbt的开路故障诊断方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103825267B (zh) | 一种mmc-mtdc直流侧短路电流的计算方法 | |
CN110581565A (zh) | 光伏发电并网系统中的控制方法和装置 | |
CN110556871B (zh) | 基于结构保持方法的大规模光伏发电系统聚合等值方法 | |
CN102969888A (zh) | 基于rtds的mmc自定义子模块的设计方法 | |
CN113258607B (zh) | 基于不同渗透下的光伏发电系统动态离散等值模型建立方法 | |
CN109145452A (zh) | 基于mmc离散建模的交流短路电流计算方法 | |
Debdouche et al. | Robust integral backstepping control microgrid connected photovoltaic System with battery energy storage through multi-functional voltage source inverter using direct power control SVM strategies | |
CN114676573A (zh) | 三电平光伏逆变器开路故障统一数学建模方法 | |
CN110429839B (zh) | 一种三相电压型pwm整流器的分数阶建模方法 | |
CN109617079B (zh) | 一种直流配电系统潮流解的存在性与稳定性分析方法 | |
Ge et al. | An effective PV power generation control system using quasi-Z source inverter with battery | |
CN113991715A (zh) | 非理想电网下中压直挂不对称混合储能系统控制方法 | |
Kaushik et al. | Performance Analysis of Boost Converters in a PV System with P and O based MPPT Controller connected to a Battery Backup and Grid | |
CN111756066B (zh) | 光伏直流变换器的运行控制与孤岛检测方法和系统 | |
Zhang et al. | Design of T-type three-level energy storage inverter and grid-connected control strategy | |
CN109861287B (zh) | 一种光伏逆变器的统一控制方法和系统 | |
Jia et al. | Transient switching performance of VSC and the DC fault partitions | |
CN105914736A (zh) | 一种配电网短路下的逆变器电源建模方法 | |
CN111525567B (zh) | 一种光伏并网逆变器故障电流的计算方法和装置 | |
CN107026475A (zh) | 一种太阳能分布式发电系统入网控制方法 | |
Jokipii et al. | Dynamic characteristics of three-phase Z-source-based photovoltaic inverter with asymmetric impedance network | |
Pandey | Stability analysis of AC/DC system with multirate discrete-time HVDC converter model | |
CN111490541B (zh) | 一种适用于极端电能质量条件下的三相逆变器锁相技术 | |
CN112904743B (zh) | 一种用于计算刚性电力电子系统的离散状态事件驱动的方法 | |
CN117614256B (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 |