CN110011312A - 基于内模的频率自适应锁相环建模方法 - Google Patents
基于内模的频率自适应锁相环建模方法 Download PDFInfo
- Publication number
- CN110011312A CN110011312A CN201910446122.6A CN201910446122A CN110011312A CN 110011312 A CN110011312 A CN 110011312A CN 201910446122 A CN201910446122 A CN 201910446122A CN 110011312 A CN110011312 A CN 110011312A
- Authority
- CN
- China
- Prior art keywords
- frequency
- formula
- internal model
- angular frequency
- sef
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/01—Arrangements for reducing harmonics or ripples
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
- Y02E40/40—Arrangements for reducing harmonics
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Stabilization Of Oscillater, Synchronisation, Frequency Synthesizers (AREA)
- Filters That Use Time-Delay Elements (AREA)
- Feedback Control In General (AREA)
Abstract
本发明公开了一种基于内模的频率自适应锁相环建模方法,涉及信号处理技术领域。所述方法包括如下步骤:构建同步旋转坐标系锁相环SRF‑PLL模型;采用频率自适应滤波器IMAF,结合正负零序分离计算模块PNZSC提取电力系统的正序基波电压,以此替代SRF‑PLL模型中的低通滤波器模块;同时采用源端解耦合归一化锁频环SUN‑FLL为所述自适应滤波器IMAF提供估算频率,最终得到IMAF‑PLL锁相环模型。所述方法具有较好的选频效果和稳定性,且能够提高锁相环的抗电压波动和频率波动的鲁棒性。
Description
技术领域
本发明涉及信号处理技术领域,尤其涉及一种基于内模的频率自适应锁相环建模方法。
背景技术
独立电力系统或独立电网的运行过程中,除了电压畸变外,还会产生由负载有功波动导致的系统频率波动或偏移。最为常用的基于同步旋转坐标系锁相环(SynchronousRotating Frame-PLL,SRF-PLL),能准确提取电压正序分量,具有良好的动态响应速度和较小的误差,但在电压畸变和频率波动同时存在的环境下,SRF-PLL的低通滤波器由于不具备频率自适应能力,在带宽的选取上存在矛盾,同步效果差;最近提出的双二阶广义积分锁相环(Double Second Order Generalized Integrator-PLL,DSOGI-PLL),其滤波器对三相不平衡、谐波电压有较好的抑制效果,且具有频率自适应特性,但仍受频率波动和电压幅值波动的影响,在独立电网频率波动的运行环境下具有较大相位误差。
采用自适应滤波器代替传统的低通滤波器,可有效提高PLL抗频率波动的鲁棒性,而锁频环FLL估算频率决定自适应滤波器的中心频率ω0是否能够准确、快速跟踪上实际角频率ω,但实际上,被测源端参数(如,系统电压幅值、系统频率)的动态变化都会影响FLL估算频率的快速性与准确性。
发明内容
本发明所要解决的技术问题是如何提供一种具有较好的选频效果和稳定性,且能够提高抗电压波动和频率波动鲁棒性的频率自适应锁相环建模方法。
为解决上述技术问题,本发明所采取的技术方案是:一种基于内模的频率自适应锁相环建模方法,其特征在于包括如下步骤:
构建同步旋转坐标系锁相环SRF-PLL模型;
采用频率自适应滤波器IMAF,通过正负零序分离计算模块PNZSC提取电力系统的正序基波电压,替代所述SRF-PLL模型中的低通滤波器模块;
同时采用源端解耦合归一化锁频环SUN-FLL为所述自适应滤波器IMAF提供估算频率,最终得到基于内模的频率自适应锁相环建IMAF-PLL模型。
进一步的技术方案在于,所述频率自适应滤波器IMAF的建模方法包括如下步骤:
设计基于内模原理的滤波器;
设计频率跟踪模块,用来实时计算或跟踪电力的系统实际角频率;
通过所述基于内模原理的滤波器以及频率跟踪模块构建频率自适应滤波器模型。
进一步的技术方案在于:已知输入信号u为交流信号,其特征数学算子为cos或sin,又要求输出信号uα和uβ相互正交,即uα⊥uβ,可知,输入信号u表达式为cos(ωt)或sin(ωt),其拉普拉斯变换为:
上式即为目标传递函数G0(s);
根据内模原理,设计传递函数如下:
其中,u(t)为输入信号,u1(t)、u2(t)为输出信号,ω为系统角频率,为微分算子,ξ1、ξ2为误差控制参数;
根据公式(2)得到基于内模原理的滤波器的数学表达式为:
其中:是指对u1(t)求一阶导,是指对u2(t)求一阶导,εu(t)是输入误差。
进一步的技术方案在于,所述频率跟踪模块的设计方法包括如下步骤:
为使估算角频率能跟随实际角频率ω的波动而波动,故令其表达式为
该式具有唯一局部平衡点,即,当时,则呈增大趋势;当时,则呈减小趋势;
根据滤波器模型现有参量输入信号u(t)、输出信号u1(t)和输出信号u2(t)来设计频率跟踪模块的f(t):采用估算角频率来做实际角频率ω,得
可知εu(t)很小,故有
根据式(3)可得
将式(6)代入式(7),可得
再将式(8)乘以输出信号u2(t),并加入一个控制系数Γ(Γ>0),并令其为由此构造出频率跟踪模型的数学表达式为
其中,
进一步的技术方案在于,所述通过所述基于内模原理的滤波器以及频率跟踪模块构建频率自适应滤波器模型的方法如下:
将所设计的滤波器中心频率以频率跟踪模块实时估算角频率来代替,即,式(3)中的角频率ω均由估算角频率来代替,最终得到基于内模原理的自适应滤波器的数学模型,如下:
其中,输入误差εu(t)和频率误差εω(t)为
估计频率微分表达式为
进一步的技术方案在于,所述源端解耦合归一化锁频环SUN-FLL的建模方法包括如下步骤:
平均归一化估算角频率,得到基于内模原理的自适应滤波器的平均估算角频率;
通过基于内模原理的自适应滤波器的平均估算角频率,得到归一化锁频环跟踪时间常数τ的数学模型;
通过将锁频环跟踪时间常数τ与源端运行状态解耦合,得到源端解耦合的动态归一化系数σall数学模型;
通过源端解耦合的动态归一化系数σall数学模型得到源端解耦合归一化锁频环数学模型。
进一步的技术方案在于,所述得到基于内模原理的自适应滤波器的平均估算角频率的方法如下:
自适应滤波器IMAF中,锁频环的估算角频率的表达式为:
其中为对估算角频率求一阶导数,Γ为控制系数(Γ>0),εu(t)为输入误差,u2(t)为滤波器输出信号,则得到α、β、γ轴上的自适应滤波器锁频环估算角频率分别为:
对αβγ轴上计算所得估算频率取平均值,并取归一化,得到基于内模原理的自适应滤波的平均估算角频率:
式中,σ为锁频环的归一化系数。
进一步的技术方案在于,计算源端解耦合归一化锁频环跟踪时间常数τ的方法如下:
将αβγ轴电压表示为正序分量、负序分量和零序分量之和:
其中,V+、V-、V0分别为电压正、负、零序分量的幅值,为负序分量的相位,为零序分量相位;ω为实际角频率;
已知自适应滤波器IMAF的传递函数如下:
其中u为输入信号,u1、u2为输出信号,且u1正交于u2,为估算角频率,ξ1、ξ2为误差控制参数;
令估算角频率与实际角频率的关系为得
其中,下标x=α、β、γ,代表αβγ三轴;
如果频率变化很小,即,m≈1,则有Gdx≈Gqx=G,可得
将上式代入式(15),得
在处对上式线性化,令avgT0[f(x)]为对函数f(x)在周期T0内取平均值,得其角频率导数的平均值为:
将式(18)代入,取m≈1,得到统一归一化锁频环数学模型:
从公式(22)中即可得到锁频环中估算角频率对实际角频率ω的跟踪时间常数τ的表达式为:
进一步的技术方案在于,获得源端解耦合的动态归一化系数σall数学模型的方法如下:
已知锁频环FLL的响应速度由跟踪时间常数τ决定,公式(23)表明,时间常数τ,除了与参数ξ1、ξ2、Γ有关,同时还与源端运行状态(被测源端频率以及被测源端电压正负零序幅值V+、V-、V0)相关;
常用归一化系数σ一般取常数,为了使锁频环FLL的响应速度与源端运行状态解耦,即,将τ与V+、V-、V0解耦,根据式(23),将锁频环归一系数σ设定为动态参数,将相关源端状态参数V+、V-、V0馈入归一化过程,得到源端解耦合的动态归一化系数σall的表达式为:
其中,被测源端频率可采用前馈,电压幅值V+、V-、V0可通过正负零序的分离计算模块获取;
将式(24)代入式(22),由此得到源端解耦合的归一化锁频环数学模型表达式为:
采用上述技术方案所产生的有益效果在于:所述方法中基于内模设计的IMAF滤波器有更好的选频效果和稳定性,对电压谐波有更好的抑制作用;同时,IMAF的频率自适应特性使滤波器在频率波动下,仍能保持良好的谐波抑制效果;且设计了基于源端解耦的改进归一化锁频环SUN-FLL,通过PNZSC得到电压幅值,以及频率前馈,计算归一化系数σall,使FLL的响应时间常数τ与被测源端参数(系统电压幅值、频率)的相关性解耦,从而提高了FLL抗电压波动和频率波动的鲁棒性,保证了IMAF滤波器中心频率对实际频率跟踪的抗干扰性。
附图说明
下面结合附图和具体实施方式对本发明作进一步详细的说明。
图1是本发明实施例中基于内模原理的滤波器结构图;
图2是本发明实施例中加入频率跟踪模块的内模原理滤波器结构;
图3是本发明实施例中基于内模原理的自适应滤波器结构图;
图4是本发明实施例中SOGI的结构图;
图5是本发明实施例中IMAF与SOGI的Gd(s)幅相频特性比较图;
图6是本发明实施例中IMAF与SOGI的Gq(s)幅相频特性比较图;
图7是本发明实施例中E(s)的幅相频特性图;
图8是基于SRF-PLL的三相四线IMAF-PLL结构图;
图9是正负零序分离计算PNZSC结构图;
图10是改进FLL的三相四线制IMAF-PLL结构图;
图11是本发明实施例所述方法的流程图。
具体实施方式
下面结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明的一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
在下面的描述中阐述了很多具体细节以便于充分理解本发明,但是本发明还可以采用其他不同于在此描述的其它方式来实施,本领域技术人员可以在不违背本发明内涵的情况下做类似推广,因此本发明不受下面公开的具体实施例的限制。
内模原理是把外部作用信号的动力学模型植入控制器来构成高精度反馈控制系统的一种设计原理。原理的核心在于:要求一个反馈控制系统能够良好抵消外部扰动,或跟踪参考输入信号,其反馈回路必须包含一个与外部输入信号相同的动力学模型,这个内部模型即称为“内模”。
基于内模原理的滤波器设计:
已知输入信号u为交流信号,其特征数学算子为cos或sin,又要求输出信号uα、uβ相互正交,即uα⊥uβ。可知,输入信号u表达式为cos(ωt)或sin(ωt),其Laplace(拉普拉斯)变换(传递函数)为:
上式即为目标传递函数G0(s)。
根据内模原理,只要设计出的控制系统传递函数G(s)与G0(s)很接近,则可抑制各种干扰对系统的影响,实现控制系统的无误差跟踪。
由于G0(s)难以实现,但可选用与其相近的控制结构,因此设计的滤波器结构如图1所示。
其中,u(t)为输入信号,u1(t)、u2(t)为输出信号,ω为系统角频率(后文将对其进行自适应处理),为微分算子,ξ1、ξ2为误差控制参数。该结构滤波器的传递函数为:
对比G(s)与G0(s),除了分母含有单次分项ξ2ωs外,其余部分都相似,同时,由于分子项与分母单次分项分别由不同的参数ξ1、ξ2控制,则可通过单独减小参数ξ2,来使G(s)接近G0(s)。
根据公式(2)得到该滤波器模型的数学表达式为:
系统角频率的自适应:
由模型传递函数,即式(2)的bode图可以看出,该滤波器的中心频率为ω,理论上ω应该为实际系统频率,但实际系统频率并不能真实获取,常规滤波器将该值选取为一个常数,即ω0=2π×50Hz=100πrad/s,但这样选取的问题在于,当系统频率发生偏移或波动时,即实际系统频率ω偏离ω0时,将导致滤波效果严重下降。因此,所述方法设计一个频率跟踪模块,用来实时计算或跟踪实际系统频率,该模块实时计算的角频率用来表示,则得到基于内模原理的自适应滤波器结构如图2所示,实时估算角频率使滤波器中心频率具有自适应能力,能自动跟随系统的波动频率,从而使滤波器在限制带宽的同时仍具有良好的响应速度。
频率跟踪模块数学模型的设计:
为使估算角频率能跟随实际频率ω的波动而波动,故令其表达式为:
该式具有唯一局部平衡点,即,当时,则呈增大趋势;当时,则呈减小趋势。可见,总是跟随实际角频率ω运动,直至收敛至
根据滤波器模型现有参量输入信号u(t)、输出信号u1(t)和输出信号u2(t)来设计频率跟踪模块的f(t):采用估算角频率来做实际角频率ω,根据图2可得:
可知εu(t)很小,故有:
根据式(3)可得:
将式(6)代入式(7),可得:
再将式(8)乘以输出信号u2(t),并加入一个控制系数Γ(Γ>0),并令其为由此构造出频率跟踪模型的数学表达式为:
满足公式(4)的要求。
由此得到内模原理的自适应滤波器结构如图3所示,其中,根据式(9)给出了频率跟踪模块的结构。
基于内模原理的自适应滤波器的数学模型:
由此,将所设计的滤波器中心频率以频率跟踪模块实时估算角频率来代替,即,式(3)中的角频率ω均由估算角频率来代替,最终得到基于内模原理的自适应滤波器(Internal Model Adaptive Filter,IMAF)的数学模型,如下:
其中,输入误差εu(t)和频率误差εω(t)为
其中,
估计频率微分表达式为:
其中,-Γ为控制参数,为了方便分析,频率跟踪部分采用单相锁频环(Single-phase Frequency Locked Loop,SFLL)表达,为输入信号u(t)的估算角频率。
令IMAF算法中两个状态变量为x1(t),x2(t),x(t)=[x1(t)x2(t)]T,写出其状态方程式:
该方程满足解的存在性与唯一性条件。采用消元法,将上式写为:
求解其状态变量x(t)=[x1(t)x2(t)]T,y=x1(t),得到微分方程的特征根为
设状态变量的初始值x(0)=[x10 x20]T,分析该微分方程的齐次解的三种情况:
(1)当ξ2=2时,有则其齐次解为:
(2)当ξ2>2时,令有:
其中,
(3)当ξ2<2时,令有:
其中,
从式(15)~(17)可以看出,当ξ2>0,该算法齐次解一致渐进稳定,可保证x(t)以指数规律衰减至特解,且增大参数ξ1和ξ2,可加快x(t)暂态过程的衰减速度。
已知输入u(t)为正弦信号,设其为:
得到特解为:
其中,故上式可写为:
可见稳态时,x1(t)与x2(t)正交,x2(t)超前x1(t)90相角;且当时,有:
滤波输出信号:
可见,输入信号u(t)经过IMAF输出u’(t)与qu’(t)相互正交,即有, 则有:
当ξ1=ξ2时,该滤波器对的正弦交流信号可实现无静态误差的跟踪,这是由于基于内模原理,该滤波器结构的Laplace变换与sin信号的Laplace变换结构相同,当ξ1>ξ2时,对该信号有增强通过的效果。
SFLL频率跟踪误差分析:
IMAF的中心频率ω0自动跟随单相锁频环SFLL的估算角频率下面对SFLL的频率跟踪误差进行分析。
当偏离实际频率ω时,将(14)代入(13),可得:
将(10)代入(11)得到εu(t),再代入(11)得到频率跟踪模块的稳态误差为:
再将(24)代入(25),得:
得到频率估算跟踪方程:
该式具有唯一局部平衡点,即,当时,则呈增大趋势;当时,则呈减小趋势。可见,总是跟随实际角频率ω运动,直至收敛至
可见,IMAF频率自适应滤波器结构稳定且收敛,其中心频率可自适应跟踪实际角频率ω的动态变化,由此,滤波器带宽可不受被测频率波动的影响,其性能与参数ξ1、ξ2有关。
与SOGI的性能比较:
双二阶广义积分锁相环DSOGI_PLL(Second Order Generalized Integrator-PLL)中的二阶广义积分器SOGI也是一个自适应滤波器,如图4所示,同样采用实时计算的频率作为中心频率。IMAF算法与SOGI算法结构相比,输入误差信号的处理不同,下面将对两者传递函数及性能做出比较。
根据图3得出IMAF滤波输出信号u’(t)、qu’(t)与输入信号u(t)的传递函数,以及SOGI滤波器输出信号x1(t)、x2(t)与u(t)的传递函数为:
两者误差εu(t)对于输入u(t)的传递函数分别为
与SOGI相比,IMAF的性能与2个参数ξ1、ξ2有关,而SOGI只有1个参数k,尤其当ξ1=ξ2=k时,IMAF与SOGI性能相同。下面通过幅相频特性图来比较两者的滤波性能,如下图5-图6所示。
可以看出:
1)IMAF算法与SOGI算法都表现出二阶带通滤波器和二阶低通滤波器的特性,且以为中心频率,都具有频率自适应能力;
2)两者随着参数值(IMAF:ξ1、ξ2;SOGI:k)的减小,选频性能表现得越好;
3)SOGI算法在处的增益总是为1,与参数k无关;IMAF算法当ξ1>ξ2时,在处比SOGI有更高的增益,在附近,比SOGI能更快速衰减,故具有更好的选频效果;
4)计算两者的幅值裕度Gm和相位裕度Pm:
IMAF(ξ1=0.05,ξ2=0.01):Gmd=Inf,Pmd=24.1°;Gmq=Inf,Pmq=113.5°
SOGI(k=0.05):Gmd=Inf,Pmd=-90°;Gmq=Inf,Pmq=-180°
可见,SOGI相位裕度不足,在状态切换时易出现振荡;而IMAF有更充足的相位裕度,稳定性更好。
将IMAF的输入误差传递函数E(s)与Gq(s)绘入同一幅幅相频特性图中,如图7所示。
可以看出,当时,εu与qu’同相,当时,εu与qu’反相,故式(11)定义εω=εu×qu’。可见,当时,当时,因而总是跟随实际角频率ω运动,直至收敛至这与上文中通过数学推导所得频率误差分析的结论一致。
正负零序分离计算:
针对三相四线制的PL-IPS,基于SRF-PLL结构设计的三相四线锁相环IMAF-PLL如下图8所示。
图中可以看出,Clark变换后的静止坐标系电压uαβγ,经过自适应滤波IMAF环节后得到基波电压分量,但该分量中仍同时包含有正序、负序和零序分量,因此,还需对其进行正负零序的分离计算(Positive Negative Zero Sequence Calculation,PNZSC),下面对PNZSC计算方法进行推导。
已知三相四线制电压uabc可分解为正序u+ abc、负序u- abc和零序电压u0 abc,可通过下式得到
其中,算子a=ej(2π/3)代表相位前移120°。三相abc坐标系转换至静止αβγ坐标系:
由此可以得到正负零序电压的分解计算式PNZSC为
可以看出,u’+ γ=0,u’- γ=0,u’0 α=0,u’0 β=0,零序分量不参与正负序分量的计算,而正负序则相互耦合。同时,根据公式(23)可得正负零序电压幅值分别为
PNZSC计算结构如图9所示,其中正序电压u’+ α、u’+ β经由Clark变换得到的d轴电压ud0即为基波正序电压。
从图8可知,单相锁频环SFLL估算频率决定IMAF的中心频率ω0是否能够准确、快速跟踪上实际角频率ω,但实际上,被测源端参数(如,PL_EPS的电压幅值、频率)的动态变化都会影响SFLL的性能。
所述方法对图8中的锁相环进行改进,设计基于源端解耦的改进归一化锁频环(Source Uncoupled Normalization FLL,SUN-FLL),通过加入归一化系数σ,使FLL的响应时间常数τ与被测源端参数的相关性解耦,从而提高FLL的快速性和抗干扰性。
自适应滤波器IMAF中,锁频环的估算角频率的表达式为:
其中为对估算角频率求一阶导数,Γ为控制系数(Γ>0),εu(t)为输入误差,u2(t)为滤波器输出信号,则得到αβγ轴上的自适应滤波器锁频环估算角频率分别为:
对αβγ轴上计算所得估算频率取平均值,并取归一化,得到基于内模原理的自适应滤波的平均估算角频率:
式中,σ为锁频环的归一化系数;
将三相四线电压表示为三相正序分量、负序分量和零序分量之和:
其中,V+、V-、V0分别为电压正、负、零序分量的幅值,为负序分量的相位,为零序分量相位;ω为实际角频率;
已知IMAF的传递函数如下:
其中u为输入信号,u1、u2为输出信号,且u1正交于u2,为估算角频率,ξ1、ξ2为误差控制参数。
令估算角频率得:
如果频率变化很小,m≈1,则有Gdα≈Gqα=G,可得:
将上式代入式(36),得:
在处对上式线性化,令avgT0[f(x)]为对函数f(x)在周期T0内取平均值,得其角频率导数的平均值为:
将式(39)代入,取m≈1,得到统一归一化锁频环数学模型:
分析公式(43),可以得到FLL中对实际角频率ω跟踪的时间常数为
可以看出,用于反映FLL对实际频率ω的跟踪速度的时间常数τ,除了与参数ξ1、ξ2、Γ有关,同时还与被测源端频率以及被测源端电压正负零序幅值V+、V-、V0相关。
归一化系数σ一般为常数,为了将FLL的响应速度与源端运行状态解耦,根据式(11),将锁频环归一系数σ设定为动态参数,将相关源端状态参数V+、V-、V0馈入归一化过程,得到:
其中,被测源端频率可采用前馈,电压幅值V+、V-、V0可通过正负零序的分离计算模块获取。
将式(45)代入式(43),由此得到源端解耦合的归一化锁频环表达式为:
式(45)表明,FLL对实际频率ω的跟踪速度与参数ξ1、ξ2、Γ有关,同时还与被测源端频率以及源端电压正负零序幅值V+、V-、V0相关。可见:
1)ξ1/ξ2虽然增大了ω处的选频增益,提高了谐波抑制效果,但同时也增大了响应时间常数τ,因此,需权衡参数ξ1/ξ2的选取;
2)在PL-EPS中源端频率产生波动,而τ随源端角频率ω的动态变化而变化;
3)在电压不平衡状态下,V+、V-、V0的动态变化也很大,也会影响频率响应速度;
综上,FLL的频率响应时间常数τ与被测源端频率、电压幅值紧密相关,可通过在计算式中引入源端动态参数,解除两者的耦合关系,达到提高FLL的抗干扰能力。
可见,所述方法的源载解耦归一化(SUN,Source Uncoupled Normalization)系数σall同时引入了被测源端两个动态参数,使FLL的频率跟踪性能与源端频率和电压幅值波动解耦,从而提高PLL的快速性和鲁棒性。
综上,如图11所示,本发明实施例公开了一种基于内模的频率自适应锁相环建模方法,包括如下步骤:
采用频率自适应滤波器IMAF,结合正负零序分离计算模块PNZSC提取电力系统的正序基波电压,同时采用源端解耦合归一化锁频环SUN-FLL为所述自适应滤波器IMAF提供估算频率,最终得到IMAF-PLL锁相环模型,其结构如图10所示。
在该锁相环中,αβγ轴自适应滤波器IMAF分别得到滤除谐波的基波电压u’αβγ、qu’αβγ;输出经PNZSC分别得到αβγ轴基波正序电压u’+ α、u’+ β,经Park变换后用于锁定相位同时,PNZSC得到的基波正负零序电压经式(3-38)得到正负零序的基波电压幅值V+、V-、V0,用以进一步计算SUN-FLL的归一化系数σ;另外,PNZSC模块还可以进一步得到d轴基波正序电压ud1。
综上,基于内模的频率自适应锁相环IMAF-PLL在SRF-PLL结构基础上,对滤波器、锁频环进行了设计和改进,具有以下优点:
(1)基于内模设计的IMAF滤波器较之SOGI滤波器有更好的选频效果和稳定性,对电压谐波有更好的抑制作用;同时,IMAF的频率自适应特性使滤波器在频率波动下,仍能保持良好的谐波抑制效果;
(2)设计了基于源端解耦的改进归一化锁频环SUN-FLL,通过PNZSC得到电压幅值,以及频率前馈,计算归一化系数σ,使FLL的响应时间常数τ与被测源端参数(系统电压幅值、频率)的相关性解耦,从而提高了FLL抗电压波动和频率波动的鲁棒性,保证了IMAF滤波器中心频率对实际频率跟踪的抗干扰性。
Claims (9)
1.一种基于内模的频率自适应锁相环建模方法,其特征在于包括如下步骤:
构建同步旋转坐标系锁相环SRF-PLL模型;
采用频率自适应滤波器IMAF,结合正负零序分离计算模块PNZSC提取电力系统的正序基波电压,替代所述SRF-PLL模型中的低通滤波器模块;
同时采用源端解耦合归一化锁频环SUN-FLL为所述自适应滤波器IMAF提供估算频率,最终得到基于内模的频率自适应锁相环建IMAF-PLL模型。
2.如权利要求1所述的基于内模的频率自适应锁相环建模方法,其特征在于,所述频率自适应滤波器IMAF的建模方法包括如下步骤:
设计基于内模原理的滤波器;
设计频率跟踪模块,用来实时计算或跟踪电力系统的实际角频率;
通过所述基于内模原理的滤波器以及频率跟踪模块构建频率自适应滤波器模型。
3.如权利要求2所述的基于内模的频率自适应锁相环建模方法,其特征在于:
已知输入信号u为交流信号,其特征数学算子为cos或sin,又要求输出信号uα和uβ相互正交,即uα⊥uβ,可知,输入信号u表达式为cos(ωt)或sin(ωt),其拉普拉斯变换为:
上式即为目标传递函数G0(s);
根据内模原理,设计传递函数如下:
其中,u(t)为输入信号,u1(t)、u2(t)为输出信号,ω为系统角频率,为微分算子,ξ1、ξ2为误差控制参数;
根据公式(2)得到基于内模原理的滤波器的数学表达式为:
其中:是指对u1(t)求一阶导,是指对u2(t)求一阶导,εu(t)是输入误差。
4.如权利要求3所述的基于内模的频率自适应锁相环建模方法,其特征在于,所述频率跟踪模块的设计方法包括如下步骤:
为使估算角频率能跟随实际角频率ω的波动而波动,故令其表达式为
该式具有唯一局部平衡点,即,当时,则呈增大趋势;当时,则呈减小趋势;
根据滤波器模型现有参量输入信号u(t)、输出信号u1(t)和输出信号u2(t)来设计频率跟踪模块的f(t):采用估算角频率来做实际角频率ω,得
可知εu(t)很小,故有
根据式(3)可得
将式(6)代入式(7),可得
再将式(8)乘以输出信号u2(t),并加入一个控制系数Γ(Γ>0),并令其为由此构造出频率跟踪模型的数学表达式为
其中,
5.如权利要求4所述的基于内模的频率自适应锁相环建模方法,其特征在于,所述通过所述基于内模原理的滤波器以及频率跟踪模块构建频率自适应滤波器模型的方法如下:
将所设计的滤波器中心频率以频率跟踪模块实时估算角频率来代替,即,式(3)中的角频率ω均由估算角频率来代替,最终得到基于内模原理的自适应滤波器的数学模型,如下:
其中,输入误差εu(t)和频率误差εω(t)为
估计频率微分表达式为
6.如权利要求1所述的基于内模的频率自适应锁相环建模方法,其特征在于,所述源端解耦合归一化锁频环SUN-FLL的建模方法包括如下步骤:
平均归一化估算角频率,得到基于内模原理的自适应滤波器的平均估算角频率;
通过基于内模原理的自适应滤波器的平均估算角频率,得到归一化锁频环跟踪时间常数τ的数学模型;
通过将锁频环跟踪时间常数τ与源端运行状态解耦合,得到源端解耦合的动态归一化系数σall数学模型;
通过源端解耦合的动态归一化系数σall数学模型得到源端解耦合归一化锁频环数学模型。
7.如权利要求6所述的基于内模的频率自适应锁相环建模方法,其特征在于所述得到基于内模原理的自适应滤波器的平均估算角频率的方法如下:
自适应滤波器IMAF中,锁频环的估算角频率的表达式为:
其中为对估算角频率求一阶导数,Γ为控制系数(Γ>0),εu(t)为输入误差,u2(t)为滤波器输出信号,则得到α、β、γ轴上的自适应滤波器锁频环估算角频率分别为:
对αβγ轴上计算所得估算频率取平均值,并取归一化,得到基于内模原理的自适应滤波的平均估算角频率:
式中,σ为锁频环的归一化系数。
8.如权利要求7所述的基于内模的频率自适应锁相环建模方法,其特征在于,计算源端解耦合归一化锁频环跟踪时间常数τ的方法如下:
将αβγ轴电压表示为正序分量、负序分量和零序分量之和:
其中,V+、V-、V0分别为电压正、负、零序分量的幅值,为负序分量的相位,为零序分量相位;ω为实际角频率;
已知自适应滤波器IMAF的传递函数如下:
其中u为输入信号,u1、u2为输出信号,且u1正交于u2,为估算角频率,ξ1、ξ2为误差控制参数;
令估算角频率与实际角频率的关系为得
其中,下标x=α、β、γ,代表αβγ三轴;
如果频率变化很小,即,m≈1,则有Gdx≈Gqx=G,可得
将上式代入式(15),得
在处对上式线性化,令avgT0[f(x)]为对函数f(x)在周期T0内取平均值,得其角频率导数的平均值为:
将式(18)代入,取m≈1,得到统一归一化锁频环数学模型:
从公式(22)中即可得到锁频环中估算角频率对实际角频率ω的跟踪时间常数τ的表达式为:
9.如权利要求8所述的基于内模的频率自适应锁相环建模方法,其特征在于,获得源端解耦合的动态归一化系数σall数学模型的方法如下:
根据式(23),将锁频环归一系数σ设定为动态参数,将相关源端状态参数V+、V-、V0馈入归一化过程,得到源端解耦合的动态归一化系数σall的表达式为:
其中,被测源端频率可采用前馈,电压幅值V+、V-、V0可通过正负零序的分离计算模块获取;
将式(24)代入式(22),由此得到源端解耦合的归一化锁频环SUN-FLL数学模型表达式为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910446122.6A CN110011312B (zh) | 2019-05-27 | 2019-05-27 | 基于内模的频率自适应锁相环建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910446122.6A CN110011312B (zh) | 2019-05-27 | 2019-05-27 | 基于内模的频率自适应锁相环建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110011312A true CN110011312A (zh) | 2019-07-12 |
CN110011312B CN110011312B (zh) | 2020-11-10 |
Family
ID=67177839
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910446122.6A Expired - Fee Related CN110011312B (zh) | 2019-05-27 | 2019-05-27 | 基于内模的频率自适应锁相环建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110011312B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111360252A (zh) * | 2020-03-05 | 2020-07-03 | 西安理工大学 | 基于3d打印机的模型辨识与内模控制系统 |
CN111537809A (zh) * | 2020-06-21 | 2020-08-14 | 陕西航空电气有限责任公司 | 一种适用于交流供电系统的数字式相序检测方法 |
CN111555752A (zh) * | 2020-05-14 | 2020-08-18 | 贵州理工学院 | 一种单相频率自适应锁相环 |
CN112260319A (zh) * | 2020-09-25 | 2021-01-22 | 国网湖南综合能源服务有限公司 | 一种基于开环的电网同步方法、系统及介质 |
CN112462140A (zh) * | 2021-02-02 | 2021-03-09 | 成都能通科技有限公司 | 一种提供电能参数分析的频率跟踪方法 |
CN116231720A (zh) * | 2023-03-28 | 2023-06-06 | 山东大学 | 新能源经柔性直流并网系统暂态稳定性提升的方法及系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150214963A1 (en) * | 2014-01-27 | 2015-07-30 | Nvidia Corporation | Phase lock loop (pll/fll) clock signal generation with frequency scaling to power supply voltage |
CN106953634A (zh) * | 2017-03-03 | 2017-07-14 | 燕山大学 | 一种基于双自调谐二阶广义积分器的锁频环方法 |
-
2019
- 2019-05-27 CN CN201910446122.6A patent/CN110011312B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150214963A1 (en) * | 2014-01-27 | 2015-07-30 | Nvidia Corporation | Phase lock loop (pll/fll) clock signal generation with frequency scaling to power supply voltage |
CN106953634A (zh) * | 2017-03-03 | 2017-07-14 | 燕山大学 | 一种基于双自调谐二阶广义积分器的锁频环方法 |
Non-Patent Citations (1)
Title |
---|
何宇等: "一种嵌入重复控制内模的三相锁相环的设计与实现", 《电工技术学报》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111360252A (zh) * | 2020-03-05 | 2020-07-03 | 西安理工大学 | 基于3d打印机的模型辨识与内模控制系统 |
CN111555752A (zh) * | 2020-05-14 | 2020-08-18 | 贵州理工学院 | 一种单相频率自适应锁相环 |
CN111555752B (zh) * | 2020-05-14 | 2023-06-16 | 贵州理工学院 | 一种单相频率自适应锁相环 |
CN111537809A (zh) * | 2020-06-21 | 2020-08-14 | 陕西航空电气有限责任公司 | 一种适用于交流供电系统的数字式相序检测方法 |
CN111537809B (zh) * | 2020-06-21 | 2022-06-07 | 陕西航空电气有限责任公司 | 一种适用于交流供电系统的数字式相序检测方法 |
CN112260319A (zh) * | 2020-09-25 | 2021-01-22 | 国网湖南综合能源服务有限公司 | 一种基于开环的电网同步方法、系统及介质 |
CN112462140A (zh) * | 2021-02-02 | 2021-03-09 | 成都能通科技有限公司 | 一种提供电能参数分析的频率跟踪方法 |
CN116231720A (zh) * | 2023-03-28 | 2023-06-06 | 山东大学 | 新能源经柔性直流并网系统暂态稳定性提升的方法及系统 |
CN116231720B (zh) * | 2023-03-28 | 2023-10-27 | 山东大学 | 新能源经柔性直流并网系统暂态稳定性提升的方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN110011312B (zh) | 2020-11-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110011312A (zh) | 基于内模的频率自适应锁相环建模方法 | |
CN107706929B (zh) | 基于最小方差滤波的自适应锁相环方法及系统 | |
CN106953634A (zh) | 一种基于双自调谐二阶广义积分器的锁频环方法 | |
CN110165709A (zh) | 考虑采样延时的虚拟同步机并网逆变系统稳定性提升方法 | |
CN109617550A (zh) | 基于二阶广义积分器的单相锁相环的控制方法 | |
CN106936427B (zh) | 一种具有谐波抑制功能的单相锁相环优化设计方法 | |
CN105425011A (zh) | 一种适用于单相畸变电网的非线性幅相检测方法 | |
CN109861296B (zh) | 一种孤岛虚拟同步发电机系统的惯性和阻尼辨识方法 | |
CN105785788A (zh) | 一种快速三相电压锁相环方法及其动态响应性能分析方法 | |
CN105680854A (zh) | 一种锁相环及锁相方法 | |
CN108155643A (zh) | 一种基于滑模观测器的单相电网电压参数的鲁棒估计方法 | |
CN108414838B (zh) | 一种逆变器并联系统线路阻抗测量方法 | |
Qi et al. | Adaptive optimal PMU placement based on empirical observability gramian | |
CN111327075B (zh) | 两级式光伏并网系统低电压穿越动态评估方法及系统 | |
CN111969597A (zh) | 带有svc的多机无穷大电力系统动态面积分滑模控制器 | |
CN106842953B (zh) | 一种无人直升机自适应低阶控制器 | |
CN108181617A (zh) | 一种基于张量积模型变换的非线性调频系统的滤波方法 | |
CN103425516B (zh) | 一种锁相环实现方法 | |
CN106227040A (zh) | 基于极值控制的伺服系统输入补偿方法 | |
CN111711204A (zh) | 一种正负序分量检测方法、控制装置及存储介质 | |
CN112511056A (zh) | 一种基于相量测量的鲁棒发电机动态状态估计方法 | |
CN110165675A (zh) | 基于内模的频率自适应滤波器建模方法 | |
CN110175401A (zh) | 源端解耦合归一化锁频环建模方法 | |
CN110266036A (zh) | 一种变流器多维频域阻抗的降维方法 | |
Zhang et al. | A distributed optimization algorithm for attack-resilient wide-area monitoring of power systems: Theoretical and experimental methods |
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 | ||
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: 20201110 Termination date: 20210527 |