CN115065266B - 一种基于分数阶lcl滤波器的单相并网逆变器建模方法 - Google Patents
一种基于分数阶lcl滤波器的单相并网逆变器建模方法 Download PDFInfo
- Publication number
- CN115065266B CN115065266B CN202210654764.7A CN202210654764A CN115065266B CN 115065266 B CN115065266 B CN 115065266B CN 202210654764 A CN202210654764 A CN 202210654764A CN 115065266 B CN115065266 B CN 115065266B
- Authority
- CN
- China
- Prior art keywords
- fractional
- order
- grid
- lcl
- lcl filter
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 26
- 238000010586 diagram Methods 0.000 claims abstract description 26
- 238000012546 transfer Methods 0.000 claims abstract description 25
- 238000004088 simulation Methods 0.000 claims abstract description 20
- 238000013178 mathematical model Methods 0.000 claims abstract description 15
- 238000013016 damping Methods 0.000 claims description 16
- 230000008569 process Effects 0.000 claims description 12
- 239000003990 capacitor Substances 0.000 claims description 10
- 238000011217 control strategy Methods 0.000 claims description 4
- 230000007423 decrease Effects 0.000 claims description 3
- 230000000694 effects Effects 0.000 description 5
- 238000011160 research Methods 0.000 description 5
- 230000009471 action Effects 0.000 description 3
- 230000002401 inhibitory effect Effects 0.000 description 3
- 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 2
- 238000013461 design Methods 0.000 description 2
- 230000005764 inhibitory process Effects 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02M—APPARATUS FOR CONVERSION BETWEEN AC AND AC, BETWEEN AC AND DC, OR BETWEEN DC AND DC, AND FOR USE WITH MAINS OR SIMILAR POWER SUPPLY SYSTEMS; CONVERSION OF DC OR AC INPUT POWER INTO SURGE OUTPUT POWER; CONTROL OR REGULATION THEREOF
- H02M7/00—Conversion of ac power input into dc power output; Conversion of dc power input into ac power output
- H02M7/42—Conversion of dc power input into ac power output without possibility of reversal
- H02M7/44—Conversion of dc power input into ac power output without possibility of reversal by static converters
- H02M7/48—Conversion of dc power input into ac power output without possibility of reversal by static converters using discharge tubes with control electrode or semiconductor devices with control electrode
- H02M7/53—Conversion of dc power input into ac power output without possibility of reversal by static converters using discharge tubes with control electrode or semiconductor devices with control electrode using devices of a triode or transistor type requiring continuous application of a control signal
- H02M7/537—Conversion of dc power input into ac power output without possibility of reversal by static converters using discharge tubes with control electrode or semiconductor devices with control electrode using devices of a triode or transistor type requiring continuous application of a control signal using semiconductor devices only, e.g. single switched pulse inverters
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02M—APPARATUS FOR CONVERSION BETWEEN AC AND AC, BETWEEN AC AND DC, OR BETWEEN DC AND DC, AND FOR USE WITH MAINS OR SIMILAR POWER SUPPLY SYSTEMS; CONVERSION OF DC OR AC INPUT POWER INTO SURGE OUTPUT POWER; CONTROL OR REGULATION THEREOF
- H02M1/00—Details of apparatus for conversion
- H02M1/0003—Details of control, feedback or regulation circuits
- H02M1/0038—Circuits or arrangements for suppressing, e.g. by masking incorrect turn-on or turn-off signals, e.g. due to current spikes in current mode control
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02M—APPARATUS FOR CONVERSION BETWEEN AC AND AC, BETWEEN AC AND DC, OR BETWEEN DC AND DC, AND FOR USE WITH MAINS OR SIMILAR POWER SUPPLY SYSTEMS; CONVERSION OF DC OR AC INPUT POWER INTO SURGE OUTPUT POWER; CONTROL OR REGULATION THEREOF
- H02M1/00—Details of apparatus for conversion
- H02M1/08—Circuits specially adapted for the generation of control voltages for semiconductor devices incorporated in static converters
- H02M1/088—Circuits specially adapted for the generation of control voltages for semiconductor devices incorporated in static converters for the simultaneous control of series or parallel connected semiconductor devices
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02M—APPARATUS FOR CONVERSION BETWEEN AC AND AC, BETWEEN AC AND DC, OR BETWEEN DC AND DC, AND FOR USE WITH MAINS OR SIMILAR POWER SUPPLY SYSTEMS; CONVERSION OF DC OR AC INPUT POWER INTO SURGE OUTPUT POWER; CONTROL OR REGULATION THEREOF
- H02M1/00—Details of apparatus for conversion
- H02M1/12—Arrangements for reducing harmonics from ac input or output
- H02M1/126—Arrangements for reducing harmonics from ac input or output using passive filters
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02M—APPARATUS FOR CONVERSION BETWEEN AC AND AC, BETWEEN AC AND DC, OR BETWEEN DC AND DC, AND FOR USE WITH MAINS OR SIMILAR POWER SUPPLY SYSTEMS; CONVERSION OF DC OR AC INPUT POWER INTO SURGE OUTPUT POWER; CONTROL OR REGULATION THEREOF
- H02M1/00—Details of apparatus for conversion
- H02M1/32—Means for protecting converters other than automatic disconnection
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Power Conversion In General (AREA)
Abstract
本发明公开了一种基于分数阶LCL滤波器的单相并网逆变器建模方法,其包括如下步骤:S1、搭建并网逆变器的分数阶电感和分数阶电容的数学模型;S2、在频域建立其数学模型并进行谐振特性分析;S3、推导分数阶LCL滤波器主电路的频域表达式,并分析谐振特性,S4、将整数阶LCL滤波器与分数阶LCL滤波器进行比较分析;S5、比较分析六种分数阶LCL滤波器的伯德图和整数阶LCL滤波器的伯德图;S6、对分数阶微分算子进行离散化近似得到拟合传递函数;S7、搭建分数阶电感和分数阶电容的分抗链电路模型;S8、分别搭建整数阶单相LCL并网逆变器电路模型和分数阶单相LCL并网逆变器电路模型,记录其仿真结果。
Description
技术领域
本发明涉及技术领域,具体涉及一种基于分数阶LCL滤波器的单相并网逆变器建模方法。
背景技术
并网逆变器是可再生能源发电单元、公共电网以及本地负载之间的能量转换接口装置,其性能好坏直接决定着并网电流质量。精确的逆变器模型实际系统的研究有重要意义。
并网逆变器通常采用脉冲宽度调制(PWM)策略,其输出的PWM电压中存在丰富的开关谐波,这会使得并网电流中存在开关谐波。针对开关谐波,通常采用体积更小,成本更低,抑制效果更好的整数阶LCL滤波器进行抑制。而整数阶LCL滤波器存在谐振尖峰,谐振尖峰将会引起系统不稳定,针对这一问题现有的策略大致分为:有源阻尼和无源阻尼两种控制方式。前者会增加成本和控制复杂度,后者会增加成本和系统损耗。
分数阶微积分是整数阶微积分的微积分阶次从整数扩展到非整数后得到的。研究表明自然界本质上是分数维的,采用分数阶微积分可以建立更为精准的系统数学模型。且国内外的研究都已经表明,电容和电感的电学性质的本质上均为分数阶维度。而针对单相LCL并网逆变器的建模大多都将电感电容视为整数阶次带入计算,未充分考虑到其系统的分数阶特性,因此对分数阶LCL滤波器的谐振特性进行分析,并将其代替传统单相并网逆变器中的整数阶LCL滤波器进行研究分析将变得极有意义。
发明内容
本发明的目的在于提供一种基于分数阶LCL滤波器的单相并网逆变器建模方法,以解决整数阶LCL滤波器存在谐振尖峰,需要采用有源阻尼或无源阻尼进行抑制这一问题。
为实现上述目的,本发明提供如下技术方案:
一种基于分数阶LCL滤波器的单相并网逆变器建模方法,包括以下操作步骤:
S1、基于分数阶微积分理论搭建并网逆变器的分数阶电感和分数阶电容的数学模型,过程如下:
式1:
式2:
其中,μ,γ分别为电感阶数和电容阶数,且0<μ,γ<2,iL为电感电流,uC为电容电压,UL为电感电压,IC为电容电流。
S2、推导LCL滤波器的主电路传递函数,在频域建立其数学模型并进行谐振特性分析,过程如下:
式3:
特别地,当μ、γ都为1时,式(3)为整数阶LCL滤波器主电路传递函数,如下:
式4:
S3、推导分数阶LCL滤波器主电路的频域表达式,并分析该表达式的谐振特性,过程如下:
令s=jω,带入式(3)化简后如下:
式5:
为简化分析,令μ2=μ1=μ,将(jω)μ=ejμπ/2=ωμcos(μπ/2)+jωμsin(μπ/2)带入式(4),化简后可以得到分数阶LCL滤波器的频域数学模型及其幅频特性和相频特性表达式:
令D=(L1+L2)/(L1L2C),μ+γ∈(0,4);
式6:
幅频表达式:
式7:
相频表达式:
式8:
式9:D=(L1+L2)/(L1L2C),μ+γ∈(0,4)。
由上式分析可得,当μ+γ∈(0,1]∪[3,4)时,即cos[(μ+γ)π/2≥0],
随着角频率从0逐渐增大,幅频特性|Ggi(jω)|的分母增大,|Ggi(jω)|减小,分数阶LCL滤波器的幅频特性不存在谐振现象。
当μ+γ∈(1,3)时,即cos[(μ+γ)π/2<0]。
假设频率
则幅频表达式化简有式10:
当μ+γ=2时,|Ggi(jω)|=∞,此时分数阶LCL滤波器的幅频特性曲线出现谐振尖峰,将μ+γ=2带入式(10),可求得谐振频率由此可看出谐振频率只有L1,L2,C的值决定,与分数阶电容电感的阶次无关。
综上,μ+γ=2为分数阶LCL滤波器存在谐振尖峰的充要条件,而整数阶LCL滤波器的元器件阶次μ=1,γ=1,μ+γ=2正好也印证了这一点。
故针对传统整数阶LCL滤波器产生谐振尖峰这一问题,可以通过引入分数阶微积分这一工具对LCL滤波器进行分数阶建模,从本质上消除谐振尖峰,使得单相并网逆变器可省略有源阻尼或无源阻尼控制,简化控制策略。
比较图3、图4可以看出整数阶单相LCL并网逆变器需要采用电容电流反馈有源阻尼来抑制谐振尖峰,从而保证系统稳定,而分数阶单相LCL并网逆变器在通过合理选取电感和电容的分数阶阶次,达到省去有缘阻尼控制,且也能使系统稳定。
S4、从幅频特性的角度将整数阶LCL滤波器与分数阶LCL滤波器进行比较分析,过程如下:
通过对电感阶次取μ=0.8,电容取γ=0.8/1.0/1.2和电容阶次取γ=0.8,电感阶次取μ=0.8/1.0/1.2六种组合构造6种不同的分数阶LCL滤波器,并推导对应分数阶LCL滤波器组合的主电路传递函数,然后在MATLAB中将其与整数阶LCL滤波器的主电路传递函数分别进行伯德图仿真。
S5、记录S4中六种分数阶LCL滤波器的伯德图和整数阶LCL滤波器的伯德图,即为图5,并从幅值裕度和相位裕度两方面对两者进行比较分析,可以看出,当μ=0.8,γ=1.2和γ=0.8,μ=1.2时,会出现谐振尖峰现象,这会导致单相并网逆变器系统不稳定,而其他4种组合,μ+γ≠2,即无谐振尖峰,并网逆变器系统稳定。
S6、由于分数阶微分算子sμ,sγ是无理函数,在数值仿真以及实际应用中不能直接实现,故先在MATALB中应用Oustaloup拟合算法对分数阶微分算子sμ,sγ进行离散化近似,得到其拟合传递函数。
S7、根据分数阶微分算子sμ,sγ的拟合传递函数和电感电容的取值在Sumilink仿真平台下搭建分数阶电感和分数阶电容的分抗链电路模型。
S8、在Sumilink仿真平台下分别搭建整数阶单相LCL并网逆变器电路模型和分数阶单相LCL并网逆变器电路模型,分别记录其仿真结果,输入电压Vdc取值为360V,电网电压Vg取值为220V,基波频率fo为50Hz,开关频率fsw取值为10Khz,采用单极性倍频SPWM脉冲调制,逆变器侧电感值L1为0.6,滤波电容C取值为10,网侧电感取值为0.15,载波幅值为3.05V,电流采样系数为Hi2取值为0.15,PI控制器参数选取Kp为0.45,Ki为2100。
优选的,所述交流侧采用由分数阶电感和分数阶电容构成的分数阶LCL滤波器,并忽略分数阶电感和分数阶电容的内阻。
优选的,所述L1为变流器侧分数阶滤波电感,其阶次为α1,单位为 L2为电网侧分数阶滤波电感,其阶次为μ1,单位为/>
优选的,所述C为交流侧分数阶滤波电容,其阶次为γ,单位为F/sec(1-γ)。阶次μ1,μ2,的取值范围均为(0,2),UDC为逆变器直流侧电压,Ui为逆变器交流侧电压,Uc分数阶滤波电容的端电压,Ug为并网点电网电压。
优选的,所述idc为直流侧电流,i1为逆变器输出电流,ic为分数阶电容电流,i2网侧电流,T1-T4为功率开关元件,D1-D4为续流二极管。
图7为整数阶单相LCL并网逆变器模型的输出电流频谱图,采用整数阶LCL滤波器抑制开关谐波,用电容电流反馈有源阻尼抑制LCL滤波器谐振尖峰,在该单相并网逆变器电路模型下,并网电流为27.31A,总谐波畸变率(THD)为3.78%,图9为分数阶单相LCL并网逆变器模型的输出电流频谱图,采用无谐振尖峰的分数阶LCL滤波器抑制开关谐波的情况下,单相并网逆变器电路模型的仿真并网电流为27.44A,THD为0.62%,采用分数阶单相LCL并网逆变器的并网电流THD为0.62%低于整数阶单相LCL并网逆变器的并网电流THD3.78%,相比单相整数阶LCL并网逆变器,采用分数阶LCL滤波器的单相并网逆变器简化系统控制的同时可以获得更好的控制效果。
与现有技术相比,本发明提供的基于分数阶LCL滤波器的单相并网逆变器建模方法,具备以下有益效果:
1、该基于分数阶LCL滤波器的单相并网逆变器建模方法,相比传统整数阶LCL滤波器,本发明借用分数阶微积分这一工具引入了分数阶LCL滤波器,并对其谐振特性进行分析,得出LCL滤波器谐振尖峰产生的本质,即电感阶数μ加电容阶数γ等于2,(μ+γ=2)。
2、该基于分数阶LCL滤波器的单相并网逆变器建模方法,通过电感阶数和电容阶数的合理选取,避开LCL滤波器谐振尖峰的产生,可省去传统抑制谐振尖峰的手段,简化系统结构设计,降低系统成本,且能取得更好的控制效果。
附图说明
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图:
图1为本发明实施例分数阶单相LCL并网逆变器主电路拓扑示意图;
图2为本发明实施例分数阶LCL滤波器数学结构示意图;
图3为本发明实施例整数阶单相LCL并网逆变器控制结构示意图;
图4为本发明实施例分数阶单相LCL并网逆变器控制结构示意图;
图5为本发明实施例整数阶LCL滤波器与分数阶LCL滤波器仿真结果对比图;
图6为本发明实施例整数阶单相LCL并网逆变器模型的输出电流仿真波形图;
图7为本发明实施例整数阶单相LCL并网逆变器模型的输出电流频频图;
图8为本发明实施例分数阶单相LCL并网逆变器模型的输出电流仿真波形图;
图9为本发明实施例分数阶单相LCL并网逆变器模型的输出电流频谱图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
在本发明中,除非另有明确的规定和限定,术语“安装”、“相连”、“连接”、“固定”等术语应做广义理解,例如,可以是固定连接,也可以是可拆卸连接,或成一体。可以是机械连接,也可以是电连接。可以是直接相连,也可以通过中间媒介间接相连,可以是两个元件内部的连通或两个元件的相互作用关系。对于本领域的普通技术人员而言,可以根据具体情况理解上述术语在本发明中的具体含义。
实施例
请参阅图1-9,本发明实施例提供的基于分数阶LCL滤波器的单相并网逆变器建模方法,包括以下操作步骤:
S1、基于分数阶微积分理论搭建并网逆变器的分数阶电感和分数阶电容的数学模型,过程如下:
式1:
式2:
其中,μ,γ分别为电感阶数和电容阶数,且0<μ,γ<2,iL为电感电流,uC为电容电压,UL为电感电压,IC为电容电流。
S2、推导LCL滤波器的主电路传递函数,在频域建立其数学模型并进行谐振特性分析,过程如下:
式3:
特别地,当μ、γ都为1时,式(3)为整数阶LCL滤波器主电路传递函数,如下:
式4:
S3、推导分数阶LCL滤波器主电路的频域表达式,并分析该表达式的谐振特性,过程如下:
令s=jω,带入式(3)化简后如下:
式5:
为简化分析,令μ2=μ1=μ,将(jω)μ=ejμπ/2=ωμcos(μπ/2)+jωμsin(μπ/2)带入式(4),化简后可以得到分数阶LCL滤波器的频域数学模型及其幅频特性和相频特性表达式:
令D=(L1+L2)/(L1L2C),μ+γ∈(0,4);
式6:
幅频表达式:
式7:
相频表达式:
式8:
式9:D=(L1+L2)/(L1L2C),μ+γ∈(0,4)。
由上式分析可得,当μ+γ∈(0,1]∪[3,4)时,即cos[(μ+γ)π/2≥0],
随着角频率从0逐渐增大,幅频特性|Ggi(jω)|的分母增大,|Ggi(jω)|减小。分数阶LCL滤波器的幅频特性不存在谐振现象。
当μ+γ∈(1,3)时,即cos[(μ+γ)π/2<0]。
假设频率
则幅频表达式化简有式10:
当μ+γ=2时,|Ggi(jω)|=∞,此时分数阶LCL滤波器的幅频特性曲线出现谐振尖峰,将μ+γ=2带入式(10),可求得谐振频率由此可看出谐振频率只有L1,L2,C的值决定,与分数阶电容电感的阶次无关。
综上,μ+γ=2为分数阶LCL滤波器存在谐振尖峰的充要条件,而整数阶LCL滤波器的元器件阶次μ=1,γ=1,μ+γ=2正好也印证了这一点。
故针对传统整数阶LCL滤波器产生谐振尖峰这一问题,可以通过引入分数阶微积分这一工具对LCL滤波器进行分数阶建模,从本质上消除谐振尖峰,使得单相并网逆变器可省略有源阻尼或无源阻尼控制,简化控制策略。
比较图3、图4可以看出整数阶单相LCL并网逆变器需要采用电容电流反馈有源阻尼来抑制谐振尖峰,从而保证系统稳定,而分数阶单相LCL并网逆变器在通过合理选取电感和电容的分数阶阶次,达到省去有缘阻尼控制,且也能使系统稳定。
S4、从幅频特性的角度将整数阶LCL滤波器与分数阶LCL滤波器进行比较分析,过程如下:
通过对电感阶次取μ=0.8,电容取γ=0.8/1.0/1.2和电容阶次取γ=0.8,电感阶次取μ=0.8/1.0/1.2六种组合构造6种不同的分数阶LCL滤波器,并推导对应分数阶LCL滤波器组合的主电路传递函数,然后在MATLAB中将其与整数阶LCL滤波器的主电路传递函数分别进行伯德图仿真。
S5、记录S4中六种分数阶LCL滤波器的伯德图和整数阶LCL滤波器的伯德图,即为图5,并从幅值裕度和相位裕度两方面对两者进行比较分析,可以看出,当μ=0.8,γ=1.2和γ=0.8,μ=1.2时,会出现谐振尖峰现象,这会导致单相并网逆变器系统不稳定,而其他4种组合,μ+γ≠2,即无谐振尖峰,并网逆变器系统稳定。
S6、由于分数阶微分算子sμ,sγ是无理函数,在数值仿真以及实际应用中不能直接实现,故先在MATALB中应用Oustaloup拟合算法对分数阶微分算子sμ,sγ进行离散化近似,得到其拟合传递函数。
S7、根据分数阶微分算子sμ,sγ的拟合传递函数和电感电容的取值在Sumilink仿真平台下搭建分数阶电感和分数阶电容的分抗链电路模型。
S8、在Sumilink仿真平台下分别搭建整数阶单相LCL并网逆变器电路模型和分数阶单相LCL并网逆变器电路模型,分别记录其仿真结果,输入电压Vdc取值为360V,电网电压Vg取值为220V,基波频率fo为50Hz,开关频率fsw取值为10Khz,采用单极性倍频SPWM脉冲调制,逆变器侧电感值L1为0.6,滤波电容C取值为10,网侧电感取值为0.15,载波幅值为3.05V,电流采样系数为Hi2取值为0.15,PI控制器参数选取Kp为0.45,Ki为2100。
交流侧采用由分数阶电感和分数阶电容构成的分数阶LCL滤波器,并忽略分数阶电感和分数阶电容的内阻。
L1为变流器侧分数阶滤波电感,其阶次为α1,单位为L2为电网侧分数阶滤波电感,其阶次为μ1,单位为/>
C为交流侧分数阶滤波电容,其阶次为γ,单位为F/sec(1-γ),阶次μ1,μ2,的取值范围均为(0,2),UDC为逆变器直流侧电压,Ui为逆变器交流侧电压,Uc分数阶滤波电容的端电压,Ug为并网点电网电压。
idc为直流侧电流,i1为逆变器输出电流,ic为分数阶电容电流,i2网侧电流,T1-T4为功率开关元件,D1-D4为续流二极管。
图7为整数阶单相LCL并网逆变器模型的输出电流频谱图,采用整数阶LCL滤波器抑制开关谐波,用电容电流反馈有源阻尼抑制LCL滤波器谐振尖峰,在该单相并网逆变器电路模型下,并网电流为27.31A,总谐波畸变率(THD)为3.78%,图9为分数阶单相LCL并网逆变器模型的输出电流频谱图,采用无谐振尖峰的分数阶LCL滤波器抑制开关谐波的情况下,单相并网逆变器电路模型的仿真并网电流为27.44A,THD为0.62%,采用分数阶单相LCL并网逆变器的并网电流THD为0.62%低于整数阶单相LCL并网逆变器的并网电流THD3.78%,相比单相整数阶LCL并网逆变器,采用分数阶LCL滤波器的单相并网逆变器简化系统控制的同时可以获得更好的控制效果。
相比传统整数阶LCL滤波器,本发明借用分数阶微积分这一工具引入了分数阶LCL滤波器,并对其谐振特性进行分析,得出LCL滤波器谐振尖峰产生的本质,即电感阶数μ加电容阶数γ等于2,(μ+γ=2),通过电感阶数和电容阶数的合理选取,避开LCL滤波器谐振尖峰的产生,可省去传统抑制谐振尖峰的手段,简化系统结构设计,降低系统成本,且能取得更好的控制效果。
本发明上述实施例,充分考虑了自然界中电感和电容的分数阶本质,借用分数阶微积分这一研究工具,对分数阶LCL滤波器的谐振特性进行了研究,相比传统整数阶LCL滤波器产生谐振尖峰这一问题,通过对分数阶电感和电容阶次的合理选取可以有效避开谐振尖峰的出现,针对LCL滤波器在单相并网逆变器中的应用,分别建立了分数阶单相LCL并网逆变器模型和整数阶单相LCL并网逆变器模型,通过理论分析和仿真结果证明:采用分数阶LCL滤波器的情况下,省略传统抑制谐振尖峰的控制策略,可以使单相并网逆变器系统稳定,且并网电流的谐波畸变率更低。
需要说明的是,在本文中,诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个......”限定的要素,并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
尽管已经示出和描述了本发明的实施例,对于本领域的普通技术人员而言,可以理解在不脱离本发明的原理和精神的情况下可以对这些实施例进行多种变化、修改、替换和变型,本发明的范围由所附权利要求及其等同物限定。
Claims (1)
1.一种基于分数阶LCL滤波器的单相并网逆变器建模方法,其特征在于,包括以下操作步骤:
S1、基于分数阶微积分理论搭建并网逆变器的分数阶电感和分数阶电容的数学模型,过程如下:
式1:
式2:
其中,μ,γ分别为电感阶数和电容阶数,且0<μ,γ<2,iL为电感电流,uC为电容电压,UL为电感电压,IC为电容电流;
S2、推导LCL滤波器的主电路传递函数,在频域建立其数学模型并进行谐振特性分析,过程如下:
式3:
当μ、γ都为1时,式3为整数阶LCL滤波器主电路传递函数,如下:
式4:
S3、推导分数阶LCL滤波器主电路的频域表达式,并分析该表达式的谐振特性,过程如下:
令s=jω,带入式3化简后如下:
式5:
为简化分析,令μ2=μ1=μ,将(jω)μ=ejμπ/2=ωμcos(μπ/2)+jωμsin(μπ/2)带入式4,化简后可以得到分数阶LCL滤波器的频域数学模型及其幅频特性和相频特性表达式:
式6:
式7:
式8:
式9:D=(L1+L2)/(L1L2C),μ+γ∈(0,4);
由上式分析可得,当μ+γ∈(1,1]∪[3,4)时,即cos[(μ+γ)π/2≥0];
随着角频率从0逐渐增大,幅频特性|Ggi(jω)|的分母增大,|Ggi(jω)|减小,分数阶LCL滤波器的幅频特性不存在谐振现象;
当μ+γ∈(1,3)时,即cos[(μ+γ)π/2<0];
频率
则幅频表达式化简有式10:
当μ+γ=2时,|Ggi(jω)|=∞,此时分数阶LCL滤波器的幅频特性曲线出现谐振尖峰,将μ+γ=2带入式10,可求得谐振频率可看出谐振频率只有L1,L2,C的值决定,与分数阶电容电感的阶次无关;
μ+γ=2为分数阶LCL滤波器存在谐振尖峰的充要条件,而整数阶LCL滤波器的元器件阶次μ=1,γ=1,μ+γ=2正好也印证了这一点;
故针对传统整数阶LCL滤波器产生谐振尖峰这一问题,通过引入分数阶微积分这一工具对LCL滤波器进行分数阶建模,从本质上消除谐振尖峰,使得单相并网逆变器可省略有源阻尼或无源阻尼控制,简化控制策略;
S4、从幅频特性的角度将整数阶LCL滤波器与分数阶LCL滤波器进行比较分析,过程如下:
通过对电感阶次取μ=0.8,电容取γ=0.8/1.0/1.2和电容阶次取γ=0.8,电感阶次取μ=0.8/1.0/1.2六种组合构造6种不同的分数阶LCL滤波器,并推导对应分数阶LCL滤波器组合的主电路传递函数,然后在MATLAB中将其与整数阶LCL滤波器的主电路传递函数分别进行伯德图仿真;
S5、记录S4中六种分数阶LCL滤波器的伯德图和整数阶LCL滤波器的伯德图,并从幅值裕度和相位裕度两方面对两者进行比较分析;
S6、由于分数阶微分算子sμ,sγ是无理函数,在数值仿真以及实际应用中不能直接实现,故先在MATALB中应用Oustaloup拟合算法对分数阶微分算子sμ,sγ进行离散化近似,得到其拟合传递函数;
S7、根据分数阶微分算子sμ,sγ的拟合传递函数和电感电容的取值在Sumilink仿真平台下搭建分数阶电感和分数阶电容的分抗链电路模型;
S8、在Sumilink仿真平台下分别搭建整数阶单相LCL并网逆变器电路模型和分数阶单相LCL并网逆变器电路模型,分别记录其仿真结果;
交流侧采用由分数阶电感和分数阶电容构成的分数阶LCL滤波器,并忽略分数阶电感和分数阶电容的内阻;
所述L1为变流器侧分数阶滤波电感,其阶次为α1,单位为L2为电网侧分数阶滤波电感,其阶次为μ1,单位为/>
所述C为交流侧分数阶滤波电容,其阶次为γ,单位为F/sec(1-γ),阶次μ1,μ2,的取值范围均为(0,2),UDC为逆变器直流侧电压,Ui为逆变器交流侧电压,Uc分数阶滤波电容的端电压,Ug为并网点电网电压;
所述idc为直流侧电流,i1为逆变器输出电流,ic为分数阶电容电流,i2网侧电流,T1-T4为功率开关元件,D1-D4为续流二极管。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210654764.7A CN115065266B (zh) | 2022-06-10 | 2022-06-10 | 一种基于分数阶lcl滤波器的单相并网逆变器建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210654764.7A CN115065266B (zh) | 2022-06-10 | 2022-06-10 | 一种基于分数阶lcl滤波器的单相并网逆变器建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115065266A CN115065266A (zh) | 2022-09-16 |
CN115065266B true CN115065266B (zh) | 2024-05-14 |
Family
ID=83201384
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210654764.7A Active CN115065266B (zh) | 2022-06-10 | 2022-06-10 | 一种基于分数阶lcl滤波器的单相并网逆变器建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115065266B (zh) |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2017099084A (ja) * | 2015-11-20 | 2017-06-01 | 株式会社明電舎 | 多慣性共振システムにおける共振抑制制御装置 |
CN206226316U (zh) * | 2016-11-25 | 2017-06-06 | 华南理工大学 | 一种阶数和容值可调的分数阶电容 |
CN107069723A (zh) * | 2016-12-08 | 2017-08-18 | 河海大学常州校区 | 有源电力滤波器分数阶PIλ及自抗扰混合控制方法 |
CN107147117A (zh) * | 2017-05-12 | 2017-09-08 | 大连理工大学 | 一种新型lcl滤波器及其设计方法 |
CN108521235A (zh) * | 2018-06-15 | 2018-09-11 | 平顶山学院 | 一种分数阶d类并联谐振逆变器及参数优化设计方法 |
CN110165920A (zh) * | 2019-05-24 | 2019-08-23 | 南京理工大学 | 基于状态空间平均法的分数阶单相逆变器建模方法 |
CN110492774A (zh) * | 2019-09-24 | 2019-11-22 | 华南理工大学 | 一种大功率分数阶阻抗元件实现电路及其控制方法 |
CN113541196A (zh) * | 2021-08-10 | 2021-10-22 | 太原理工大学 | 一种单相lc型并网逆变器分数阶控制方法 |
-
2022
- 2022-06-10 CN CN202210654764.7A patent/CN115065266B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2017099084A (ja) * | 2015-11-20 | 2017-06-01 | 株式会社明電舎 | 多慣性共振システムにおける共振抑制制御装置 |
CN206226316U (zh) * | 2016-11-25 | 2017-06-06 | 华南理工大学 | 一种阶数和容值可调的分数阶电容 |
CN107069723A (zh) * | 2016-12-08 | 2017-08-18 | 河海大学常州校区 | 有源电力滤波器分数阶PIλ及自抗扰混合控制方法 |
CN107147117A (zh) * | 2017-05-12 | 2017-09-08 | 大连理工大学 | 一种新型lcl滤波器及其设计方法 |
CN108521235A (zh) * | 2018-06-15 | 2018-09-11 | 平顶山学院 | 一种分数阶d类并联谐振逆变器及参数优化设计方法 |
CN110165920A (zh) * | 2019-05-24 | 2019-08-23 | 南京理工大学 | 基于状态空间平均法的分数阶单相逆变器建模方法 |
CN110492774A (zh) * | 2019-09-24 | 2019-11-22 | 华南理工大学 | 一种大功率分数阶阻抗元件实现电路及其控制方法 |
CN113541196A (zh) * | 2021-08-10 | 2021-10-22 | 太原理工大学 | 一种单相lc型并网逆变器分数阶控制方法 |
Non-Patent Citations (2)
Title |
---|
Grid-connected inverter based on a resonance-free fractional-order LCL filter;Xiaogang Wang, et al;《FRACTAL AND FRACTIONAL》;20220701;第6卷(第7期);第1-19页 * |
分数阶LCL并网逆变器建模与控制研究;蔡俊辉;《中国知网学位论文全文数据库》;20240215;正文第21-39页 * |
Also Published As
Publication number | Publication date |
---|---|
CN115065266A (zh) | 2022-09-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Li et al. | Active DC-link power filter for single phase PWM rectifiers | |
Liu et al. | An adaptive Discontinuous Pulse Width Modulation (DPWM) method for three phase inverter | |
Cossoli et al. | Proportional-resonant controller and LCL filter design for single-phase grid-connected PV micro-inverters | |
CN115313346A (zh) | 单相电压型逆变器整流器负载输入电流低频谐波抑制方法 | |
CN114142751B (zh) | 一种不平衡电网电压下三相csr比例积分谐振控制方法 | |
Aly et al. | A model predictive control method for common grounded photovoltaic multilevel inverter | |
CN113541196A (zh) | 一种单相lc型并网逆变器分数阶控制方法 | |
Liu et al. | Model predictive control of single phase grid-connected inverter with LC filter | |
Delpino et al. | Line harmonics on systems using reduced DC-link capacitors | |
Kraemer et al. | Robust design of control structure for three-phase grid-tied inverters | |
CN112217381B (zh) | 一种单相变换器直流母线电压二次纹波抑制电路和方法 | |
Farhadi et al. | An Improved Model Predictive Control Method to Drive an Induction Motor Fed by Three-Level Diode-Clamped Indirect Matrix Converter. | |
CN115065266B (zh) | 一种基于分数阶lcl滤波器的单相并网逆变器建模方法 | |
CN111355257A (zh) | 一种弱电网下光伏并网逆变器无源准pr控制方法 | |
CN106532703B (zh) | 储能变流系统与变流器及双向变流器的控制方法和控制器 | |
CN210167979U (zh) | 抑制共模电流的非隔离单相光伏并网逆变器 | |
Shen et al. | Design and modeling of a two-stage pv inverter for single phase and three phase applications | |
Khan et al. | Electrolytic capacitor-less dual buck inverter with CM and DM active resonance damping control for Non-isolated grid-connected PV applications | |
Abusara et al. | Design of a robust digital current controller for a grid connected interleaved inverter | |
CN116915033A (zh) | 一种并网逆变lcl滤波器优化控制方法 | |
Mohamadian et al. | LCL filter design for single-phase grid-connected PV inverters with double-frequency unipolar PWM switching | |
Alsolami et al. | DC ripple current reduction on multilevel, multiport, single-phase DC/AC converter for renewable energy applications | |
Tyagi et al. | Third harmonic compensation in a single-phase H-bridge PFC | |
Xu et al. | Current control based on zero-placement strategy for grid-connected LCL-filtered inverters | |
Widagdo et al. | Harmonic Mitigation Using Passive Filters in 3-Phase Inverters to Improve Power Quality on Microgrid |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |