CN104113052A - 有源电力滤波器的保护方法 - Google Patents
有源电力滤波器的保护方法 Download PDFInfo
- Publication number
- CN104113052A CN104113052A CN201410305638.6A CN201410305638A CN104113052A CN 104113052 A CN104113052 A CN 104113052A CN 201410305638 A CN201410305638 A CN 201410305638A CN 104113052 A CN104113052 A CN 104113052A
- Authority
- CN
- China
- Prior art keywords
- current
- mtd
- active power
- mtr
- power 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 37
- 230000000670 limiting effect Effects 0.000 claims abstract description 54
- 230000003750 conditioning effect Effects 0.000 claims description 15
- 230000000903 blocking effect Effects 0.000 claims description 13
- 238000004422 calculation algorithm Methods 0.000 claims description 12
- 238000004364 calculation method Methods 0.000 claims description 9
- 230000001143 conditioned effect Effects 0.000 claims description 2
- 239000000284 extract Substances 0.000 claims description 2
- 230000009977 dual effect Effects 0.000 abstract 1
- 238000005070 sampling Methods 0.000 description 40
- 230000009466 transformation Effects 0.000 description 8
- 239000003990 capacitor Substances 0.000 description 7
- 238000001514 detection method Methods 0.000 description 7
- 230000000694 effects Effects 0.000 description 7
- 230000008569 process Effects 0.000 description 7
- 238000013461 design Methods 0.000 description 5
- 230000007547 defect Effects 0.000 description 4
- 230000001105 regulatory effect Effects 0.000 description 4
- 230000001360 synchronised effect Effects 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 230000007935 neutral effect Effects 0.000 description 3
- 238000012795 verification Methods 0.000 description 3
- 238000012935 Averaging Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 230000000737 periodic effect Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- BRUQQQPBMZOVGD-XFKAJCMBSA-N Oxycodone Chemical compound O=C([C@@H]1O2)CC[C@@]3(O)[C@H]4CC5=CC=C(OC)C2=C5[C@@]13CCN4C BRUQQQPBMZOVGD-XFKAJCMBSA-N 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 238000011217 control strategy Methods 0.000 description 1
- 230000001276 controlling effect Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000010183 spectrum analysis Methods 0.000 description 1
- 230000000087 stabilizing effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
Classifications
-
- 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/20—Active power filtering [APF]
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种有源电力滤波器的保护方法,本发明提出的有源电力滤波器的保护方法中通过对指令电流进行限流,使得有源电力滤波器实际输出的补偿电流限制在合理范围内,防止由于指令电流过大造成APF过流而损坏有源电力滤波器,且本发明中对指令电流依次进行比例限流和截断限流,通过两次限流实现双重保护,提高了***的可靠性,能够更好的保护有源电力滤波器,延长使用寿命。
Description
技术领域
本发明属于电力电子技术领域,具体涉及一种有源电力滤波器的保护方法。
背景技术
电力电子技术的日益发展成熟和电力自动化设备的广泛应用在推进工业生产力发展和电气自动化水平提高的同时,也引起了严重的电网谐波污染问题,有源电力滤波器(APF)应运而生。
有源电力滤波器(APF)通过包括主电路(变流器)、驱动电路和控制***。其补偿过程为先通过电压互感器(PT)、电流互感器(CT)获取***电压和非线性负载的电流信号,控制***对获取的电压和电流信号进行处理,计算出需要补偿的各次谐波对应的指令电流,再根据该指令电流,生成用于驱动变流器的驱动脉冲,实时快速的控制主电路中开关器件(功率器件)的通断,输出实际需要补偿电流,该补偿电流与负载电流中需要补偿的各次谐波进行抵消,最终得到期望的正弦化的电网电流,从而改善了电网的电能质量。
作为一个电网谐波和无功的动态补偿装置,APF不仅需要有良好的稳态补偿性能和动态响应速度,最主要的是要保证能稳定可靠的运行。APF稳定运行的前提是其补偿容量不能超过自身的额定容量,因为其变流器中功率单元对瞬时过流和过压非常敏感,短时的过流和过压也可能会导致开关管的损坏,所以要尽量避免APF的过容运行。
正常情况下,可能造成APF过容运行的情况有两种:(1)接入公用电网的非线性负载突增;(2)接入公用电网的负载发生短路或断路等故障,造成***不对称分量突增。这两种情况如果造成APF过容运行,其根本原因是负载电流谐波分量的突增导致指令电流的突增。所以要限制APF的过容运行,最关键的是要限制指令电流在自身补偿容量范围之内。
杨振宇等人在题为并联型有源电力滤波器限流补偿策略研究(电力自动化设备.2006,26(3):21-25)一文中提出一种软件限流控制策略,只对参考电流进行处理,控制实现简单,并取得了较好的限流效果。但它存在一个最关键的问题尚未解决:指令电流最大绝对值irefmax和***最大允许补偿电流值imax的定义。传统的截断限流和比例限流策略对于指令电流最大绝对值irefmax的定义就是指令电流峰值的绝对值,而***最大允许补偿电流值imax则根据***容量被定义为一个固定的值。
实际上,对于一个能够补偿0~2.5kHz以内任意次谐波的APF***来说,指令电流最大绝对值irefmax绝不能简单的定义为指令电流峰值的绝对值。除了要考虑补偿全频段谐波时指令电流的有效值达到***额定容量时的峰值,还要考虑补偿低次或单次谐波时指令电流的有效值达到***额定容量时的峰值。通常APF的额定容量是以每相能够补偿的额定谐波电流有效值来定义的。结合实际情况,当APF需要补偿全频段谐波时,若***的补偿电流有效值达到额定值,此时指令电流的峰值往往要比APF***只补偿低次或单次谐波而补偿电流有效值达到额定值时的指令电流峰值高得多。所以,若简单的将指令电流最大值irefmax定义为峰值,那么可能出现下面两种情况:
(1)若按APF补偿全频段谐波时有效值达到额定值时指令电流的峰值定义为***最大允许补偿电流峰值ipeak,那么当APF补偿低次或单次谐波时,即使补偿电流的有效值超过了额定值,由于此时指令电流的峰值未达到***设定的最大允许补偿电流峰值ipeak,指令电流会继续增大,直到峰值达到设定的ipeak,而此时指令电流的有效值已远远超过额定有效值,从而造成APF过流。
(2)若按APF补偿低次或单次谐波时有效值达到额定值时指令电流的峰值定义为***最大允许补偿电流峰值ipeak(即最大允许电流),那么当APF补偿全频段谐波时,由于当指令电流的峰值达到设定的ipeak时,指令电流的有效值还远远未达到***额定值,从而造成APF利用率不高。
同时,获得用于比例限流环节的指令电流有效值及实现APF的选择性谐波补偿需要采用一种有效的谐波检测方法。张树全等人在题为多同步旋转坐标系下指定次谐波电流控制(中国电机工程学报.2010(03):55-62)提出一种基于同步旋转坐标变换的选择性谐波电流检测算法,但是该方法需要经过多次旋转坐标变换,计算量大,并且需要设计一个低通或带通滤波器,谐波分量的检测精度取决于滤波器的设计,滤波器的参数设计不当会直接影响检测精度。
此外,大容量非线性负载突变时会引起APF直流侧电压波动,若直流侧出现过压,则并联在直流母线电容两端的开关器件有可能因为瞬态高压击穿而导致损坏。
发明内容
针对现有技术中存在的不足,本发明提供了一种有源电力滤波器的保护方法。
本发明的有源电力滤波器的保护方法,包括如下步骤:
(1)提取负载电流中需要补偿的各次谐波分量,并将提取的各次谐波分量叠加得到指令电流;
(2)根据有源电力滤波器的额定补偿电流有效值和需要补偿的最高谐波次数依次设定第一阈值和第二阈值;
(3)将第一阈值作为最大允许电流对所述的指令电流进行比例限流,得到第一指令电流;
(4)将第二阈值作为截断峰值电流对所述的第一指令电流进行截断限流,得到第二指令电流;
(5)根据第二指令电流生成用于驱动有源电力滤波器中变流器的驱动脉冲,使变流器输出相应的补偿电流。
本发明提出的有源电力滤波器的保护方法中通过对指令电流进行限流,使得有源电力滤波器实际输出的补偿电流限制在合理范围内,防止由于指令电流过大造成APF过流(指实际输出的补偿电流超过开关器件所能承受的最大允许值)而损坏有源电力滤波器,且本发明中对指令电流依次进行比例限流和截断限流,通过两次限流实现双重保护,提高了***的可靠性,能够更好的保护有源电力滤波器,延长使用寿命。
在进行比例限流和截断限流处理前,先根据用户设定的补偿频次设定比例限流和截断限流的相关参数。这样根据实际情况设定相关参数,与现有技术中任何工况下都固定使用单一参数比较,更为灵活,在APF仅补偿低次谐波的情况下,可以进一步提高APF的输出能力,改善了APF的利用效率,即使在全频段补偿的情况下,也能有效地将补偿电流限制在额定输出能力范围之内。该种参数可调型的比例限流和截断限流方法有助于实现不同工况下APF输出能力最大化和高可靠性的统一。
作为优选,所述步骤(1)基于滑窗迭代DFT算法提取需要补偿的各次谐波分量。
基于滑窗迭代DFT算法克服了基于同步旋转坐标系的谐波电流检测算法中坐标变换复杂、参数设计困难的缺陷,并通过取消周期平均滤波器对传统滑窗DFT算法进行改进,大大减小了对程序空间资源的占用,实现了任意次谐波的单次或者叠加提取,具有很高的精度。
基于滑窗迭代DFT算法提取n次谐波分类的具体过程如下:
a初始化一个滑窗,所述滑窗的长度为一个基波周期内的采样点的总数,记为N,假设滑窗采样的数据存储在一个大小为N的数组X[N]中,初始条件下滑窗中所有的元素的值为0,即X[0]~X[N-1]=0;
b采集一个新的数据x(kτ)(x(kτ)为某一时刻的采样数据,其中k表示当前的采样点编号,k的取值范围为0~N-1,τ为采样时间间隔),更新滑窗一个基波周期前该采样点对应的采样数据,;
c.计算需要提取的n次谐波在当前采样点的相位,具体为:
一个周期是2π,一个周期采样N个采样点,故每个采样点之间的相位差为本发明中假设第一个采样点处的相位为0,则第k个采样点处的相位就是由于是n次谐波,则第k个采样点处n次谐波的相位为
d.判断当前采样点的编号k是否满足条件k=N-1,N为一个基波周期内的采样点的总数:
若满足,则令:
当前采样点处用于谐波分量瞬时值计算的n次谐波分量的实部An和虚部Bn分别为:An=an,Bn=bn;
当前采样点处n次谐波分量瞬时值xn(k)为:
xn(k)=ancos(2πnk/N)+bnsin(2πnk/N),
其中,
an=an+x(kτ)*cos(2πnk/N),
bn=bn+x(kτ)*sin(2πnk/N)。
否则,令:
当前采样点处用于计算的n次谐波分量的瞬时值的实部An和虚部Bn分别为:An=A'n,Bn=B'n,其中A'n和B'n分别为上一基波周期n次谐波分量的实部An和虚部Bn,即计算当前采样点处的谐波分量所用的实部和虚部不进行更新,继续沿用上一周期的值;
进一步,根据实部An和虚部Bn得到当前采样点处n次谐波分量的瞬时值xn(k)为:
xn(k)=Ancos(2πnk/N)+Bnsin(2πnk/N),
e.采样点数k加1后,返回继续执行步骤b,直至k=N时结束。
所述步骤(2)根据如下公式计算第一阈值imax:
imax=Ki*ip,
其中,ip为额定补偿电流有效值;
第一比例系数为Ki:
其中,h为需要补偿的最高谐波次数。
所述步骤(2)根据如下公式计算第二阈值ipeak:
ipeak=Kp*imax,
其中,imax为第一阈值;
Kp为第二比例系数:
其中,h为需要补偿的最高谐波次数。
根据需要补偿的最高谐波次数设定第一比例系数和第二比例系数使APF运行效率最大化原理,保证APF在各种不同工况下既能安全可靠地稳定运行,又能实现变流器的开关器件(功率单元)的利用效率达到最高。具体选取过程如下:先通过Matlab的Simulink仿真功能搭建仿真相关APF模型,根据设定的需要补偿的谐波次数以及对应补偿谐波次数下的谐波电流波形特性(主要是谐波电流有效值和谐波电流峰值的关系)初始化设置该两个比例系数(即第一比例系数和第二比例系数,还包括频段如何划分和各个频段内比例系数值的初值选取),根据仿真的APF模型在初始化比例系数值的作用下实际峰值电流是否超过最大允许补偿电流峰值以及实际补偿电流有效值是否超过最大允许补偿电流有效值,若未超过,可适当加大比例系数再进行效果验证;反之,若超过,则适当减小进行效果验证。经过多次反馈调节和验证后可以得到最优的两个比例系数,本发明提供的梯度式设定值可以作为优选。
所述步骤(3)过程如下:
(3-1)计算指令电流的有效值,并用第一阈值除以指令电流的有效值和第一阈值中的较大者得到的商作为比例限流因子;
(3-2)根据公式:
i'ref=k*iref
计算第一指令电流i'ref,iref为指令电流,k为比例限流因子。
取指令电流的有效值irms和第一阈值imax中的较大者,记为iMax,再用第一阈值imax除以该较大者iMax得到的商作为比例限流因子k,即k=imax/iMAX。当指令电流的有效值大于第一阈值时时,iMAX=irms此时比例限流因子k小于1,进而得到的第一指令电流i'ref小于指令电流iref(原指令电流),即起到限流的作用。当指令电流的有效值小于第一阈值时,iMAX=imax,此时比例效率因子k为1,进而得到的第一指令电流等于指令电流,相当于将原先计算得到的指令电流原样进行输出。
本发明中根据:
计算指令电流的有效值irms,N为一个基波周期内的采样点数总数,h为最高补偿谐波次数,n表示补偿谐波次数。
所述步骤(4)中根据公式:
对第一指令电流进行截断限流,。
以比例限流环节的输出作为截断限流环节的输入,从而构成先比例限流后截断限流的双重限流,经过截断限流之后可以得到第二指令电流i″ref,该第二指令电流用于产生驱动变流器的驱动脉冲通过截断限流将第一指令电流的幅值限定在[-ipeak,ipeak]之间,进一步防止补偿电流过大。
所述的保护方法还包括实时检测有源电力滤波器的直流侧电压(即直流母线电容两端的电压),根据直流侧电压对有源电力滤波器进行硬件封锁限流,具体过程如下:
(S1)实时检测有源电力滤波器的直流侧电压;
(S2)将检测到的直流侧电压调理到设定的电压范围内得到相应的等效直流侧电压;
(S3)将调理得到的等效直流侧电压与第三阈值进行比较,并根据比较结果进行如下操作:
若大于或等于第三阈值,封锁有源电力滤波器中变流器的驱动脉冲,直至等效直流侧电压下降至第四阈值时解除封锁,使变流器输出的相应的补偿电流;
否则,则不操作。
由于负载突变很容易出现直流侧过压的情况,为了调节直流侧电压,使其恢复到正常水平,控制器会在原有的根据负载谐波电流检测得出的指令电流中叠加一个基波有功分量,两个分量一叠加,有可能超过APF的最大补偿电流允许值,从而造成APF过流,但是该种过流情况可以通过本发明提供的指令限流加以有效限制。但是在负载波动时,由于指令电流的变化及跟踪通常会存在一定延时,所以依靠指令电流中叠加一个基波有功分量来来完成直流侧电压的控制,需要一定的延时,很有可能指令电流还来不及调节时就由于直流侧过压,已经造成APF损坏了。本发明中通过实时检测直流侧电压,当检测到直流侧过压时直接封锁APF中功率模块的驱动脉冲,防止过压造成开关器件的损坏,同时直流侧由于二极管的续流作用也会逐渐下降,从而完成直流侧电压的调节,能够有效消除延时,实时性强,对冲击性电流具有良好的抑制作用,且控制灵活,实现简单。
硬件封锁限流的原理如下:直流侧电压瞬时值的等效电压值超过最大设定值(第三阈值)时,将驱动脉冲强制拉低(即封锁有源电力滤波器中变流器的驱动脉冲),封锁变流器正在导通的开关器件(功率单元),使电路进入续流状态,通过开关器件两端的二极管的续流作用使得直流侧电压逐渐降低,当直流侧电压瞬时值的等效电压值下降到最小设定值(第四阈值)时,驱动脉冲重新开放,功率单元恢复正常工作,进而起到直流侧过压情况下保护有源电力滤波器的作用。
由于APF正常运行情况下采集到的直流侧电压通常为较高的电压(大约700±50V左右),无法进行直接的比较保护,为便于实现,本发明中将采集到的有源电力滤波器的直流侧电压调理到设定的电压范围内得到相应的等效直流侧电压;然后利用等效直流侧电压完成比较以确定是否需要进行硬件封锁限流。本发明中根据实际应用需求设定合适的等效调理变比KT,通常情况下KT的取值范围为0.004~0.006。利用设定的等效调理变比KT根据以下公式将检测到的直流侧电压Vdc调理到设定的电压范围内得到等效直流侧电压V′dc:
V′dc=KT*Vdc。
实际上,若技术条件允许,也可以不进行调理,直接以采集到的直流侧电压进行后续的比较。
所述的第三阈值V1根据公式:
V1=KT*(Vdcp+ΔVh)
计算,其中,Vdcp表示直流侧电压额定值,ΔVh表示直流侧过压保护的过压裕量,KT表示等效调理变比。
所述的第四阈值V2根据公式:
V2=KT*(Vdcp-ΔVl)
计算,Vdcp表示直流侧电压额定值,ΔVl表示直流侧过压保护解除的欠压裕量,KT表示等效调理变比。
直流侧过压保护的过压裕量和直流侧过压保护解除的欠压裕量的选取会直接影响到保护的灵敏度,选取的原则一方面要考虑保护的可靠性,故取值不能太大,否则敏感度太低;另一方面也不能取的太小,否则有可能会发生误动作,影响APF工作效率。作为优选,直流侧过压保护的过压裕量ΔVh可取50~100V,直流侧过压保护解除的欠压裕量ΔVl可取100~150V。
与现有技术相比,本发明存在如下优点:
(1)综合了比例限流和截断限流各自的优势,采用先比例限流后截断限流的双重保护保证了PWM驱动脉冲的可靠性。
(2)综合考虑了APF补偿全频段谐波和只补偿低次或者单次谐波时的电流波形特性,根据实际应用情况灵活设定比例限流和截断限流的阈值,既能有效的限制指令电流的大小,进而起到对APF输出的补偿电流的限流作用,又能保证各类工况下APF的利用效率。
(3)在对指令电流依次进行比例限流和截断限流的同时,根据APF的直流侧电压进行基于硬件封锁的保护方法,有效消除了利用指令电流中叠加一个基波有功分量调节直流侧电压时存在一定延时的缺陷,能迅速地保护变流器的开关器件,以防止直流侧过压损开关器件,对冲击性负载突变具有良好的抑制作用,且控制灵活,实现简单。
(4)使用一种改进的滑窗迭代DFT算法克服了基于同步旋转坐标系的谐波电流检测算法坐标变换复杂、参数设计困难的缺陷,并通过取消周期平均滤波器对传统滑窗DFT算法进行改进,大大减小了对程序空间资源的占用,实现了任意次谐波的单次或者叠加提取,具有很高的精度。
附图说明
图1本实施例的有源电力滤波器的保护方法的流程图;
图2为本实施例的基准电压电路;
图3为本实施例的迟滞比较电路;
图4为本实施例的SAPF***结构示意图;
图5为本实施例的SAPF在额定条件下的全频段稳态补偿波形;
图6为本实施例的SAPF进行限流补偿的稳态补偿波形;
图7为本实施例的SAPF进行限流补偿的动态补偿波形。
具体实施方式
下面结合附图和具体实例对本发明作进一步的说明。
本实施例的有源电力滤波器的保护方法,如图1所示,包括如下步骤:
(1)提取负载电流中需要补偿的各次谐波分量,并将提取的各次谐波分量叠加得到指令电流;
本实施例中基于滑窗迭代DFT算法提取需要补偿的各次谐波分量,基于滑窗迭代DFT算法提取n次谐波分类的具体过程如下:
a初始化一个滑窗,滑窗的长度为一个基波周期内的采样点的总数,记为N(本实施例中N=200),假设滑窗采样的数据存储在一个大小为N的数组X[N]中,初始条件下滑窗中所有的元素的值为0,即X[0]~X[N-1]=0;
b采集一个新的数据x(kτ),k表示一个基波周期内的当前采样点对应的编号,τ为采样时间间隔(为1/N倍的基波周期),更新滑窗一个基波周期前该采样点对应的采样数据,;
c.计算需要提取的n次谐波在当前采样点的相位,具体为:
一个周期是2π,一个周期采样N个点,故每个采样点之间的相位差为本发明中假设第一个采样点处的相位为0,则第k个采样点处的相位就是由于是n次谐波,所以第k个采样点处n次谐波的相位为
d.判断当前采样点的编号k是否满足条件k=N-1,N为一个基波周期内的采样点的总数:
若满足,则令:
当前采样点处用于谐波分量瞬时值计算的n次谐波分量的实部An和虚部Bn分别为:An=an,Bn=bn;
当前采样点处n次谐波分量瞬时值为:
xn(k)=ancos(2πnk/N)+bnsin(2πnk/N),
其中,
an=an+x(kτ)*cos(2πnk/N),
bn=bn+x(kτ)*sin(2πnk/N),
否则,令:
当前采样点处用于谐波分量瞬时值计算的n次谐波分量的实部An和虚部Bn分别为:An=A'n,Bn=B'n,其中A'n和B'n分别为上一基波周期n次谐波分量的实部An和虚部Bn,即计算当前采样点处的谐波分量所用的实部和虚部不进行更新,继续沿用上一周期的值;
当前采样点处n次谐波分量的瞬时值为:
xn(k)=Ancos(2πnk/N)+Bnsin(2πnk/N),
e.采样点数k加1后,返回继续执行步骤b,直至k=N时结束。
每次循环结束,可以得到当前采样点处n次谐波分量(即n次谐波分量的瞬时值),以及对应的An和Bn。
(2)根据有源电力滤波器额定补偿电流有效值和需要补偿的最高谐波次数依次设定第一阈值和第二阈值;
根据如下公式计算第一阈值imax:
imax=Ki*ip,
其中,ip为额定补偿电流有效值;
Ki为第一比例系数:
其中,h为需要补偿的最高谐波次数。
根据如下公式计算第二阈值ipeak:
ipeak=Kp*imax,
其中,imax为第一阈值;
Kp为第二比例系数:
其中,h为需要补偿的最高谐波谐波次数。
(3)将第一阈值作为最大允许电流对指令电流进行比例限流,得到第一指令电流;
(3-1)根据:
计算指令电流的有效值irms,N为一个基波周期内的采样点数总数,h为需要补偿的最高谐波谐波次数,n表示需要补偿的谐波谐波次数(只补偿奇数次,即n取奇数);
并用第一阈值除以指令电流的有效值和第一阈值中的较大者得到的商作为比例限流因子;
(3-2)根据公式:
i'ref=k*iref
计算第一指令电流i'ref,iref为指令电流,k为比例限流因子。
(4)以第二阈值作为截断峰值电流对第一指令电流进行截断限流,得到第二指令电流;具体如下根据公式:
对第一指令电流进行截断限流,即以比例限流环节的输出作为截断限流环节的输入,从而构成先比例限流后截断限流的双重限流,经过截断限流之后可以得到第二指令电流i″ref,该i″ref即是最终的参考指令电流,通过截断限流将第一指令电流的幅值限定在[-ipeak,ipeak]之间,进一步防止补偿电流过大。
(5)根据第二指令电流生成用于驱动有源电力滤波器的变流器的驱动脉冲,使变流器输出相应的补偿电流。
本实施例的保护方法还包括实时检测有源电力滤波器的直流侧电压(直流母线电容两端的电压),根据直流侧电压对有源电力滤波器进行硬件封锁限流,具体过程如下:
(S1)实时采集有源电力滤波器的直流侧电压;
(S2)将检测到的直流侧电压调理到设定的电压范围内得到相应的等效直流侧电压,具体如下:
本实施例中将检测到的直流侧电压调理到5V以下,设定等效调理变比KT=0.00517,根据以下公式将检测到的直流侧电压Vdc调理到设定的电压范围内得到等效直流侧电压V′dc:
V′dc=KT*Vdc。
(S3)将调理得到的直流侧等效电压与第三阈值进行比较,并根据比较结果进行如下操作:
若大于或等于第三阈值,封锁有源电力滤波器中变流器的驱动脉冲,直至等效直流侧电压下降至第四阈值时解除封锁,进而完成对变流器输出的补偿电流的调节;
否则,则不操作。
本实施例中直流侧电压额定值为750V,直流侧过压保护的过压裕量ΔVdc取100V,直流侧过压保护解除的欠压裕量Δ'Vdc取150V,故实际的直流侧过压保护值为850V,过压保护解除值为600V。
根据公式:
V1=KT*(Vdcp+ΔVh)
计算第三阈值V1,其中,Vdcp表示直流侧电压额定值,ΔVh表示直流侧过压保护的过压裕量,KT表示等效调理变比,得到第三阈值V1=4.4V。
根据公式:
V2=KT*(Vdcp-ΔVl)
计算第四阈值V2,Vdcp表示直流侧电压额定值,ΔVl表示直流侧过压保护解除的欠压裕量,KT表示等效调理变比,得到第四阈值V1=3.1V。
硬件封锁限流基于以下***完成,该***包括电压调理电路、参考电压电路、迟滞比较电路、驱动保护电路,其中:
电压调理电路,用于将采集到的直流侧电压调理到设定的电压范围内输出等效直流侧电压;
参考电压电路,用于根据设定电压范围生成相应的参考电压;
迟滞比较电路,以参考电压和等效直流侧电压为输入,输出相应的电平信号且迟滞比较电路的对应的阈值分别为第三阈值和第四阈值,即分别为4.4V和3.1V;
驱动保护电路,接收迟滞比较电路的输出,当迟滞比较电路的输出电平为高电平时输出驱动保护信号,该驱动保护信号用于封锁APF中变流器的驱动脉冲。
参考电压电路如图2所示,包括稳压二极管U1,稳压二极管U1的第一引脚1通过限流电阻R1与电源DVCC相连,稳压二极管U1的第二引脚2接地GND,第一引脚1和第二引脚2之间连接有旁路电容C3,稳压二极管U1的第三引脚3通过第一分压电阻R2与第一引脚连接,通过第二分压电阻R3与第二引脚连接,为实现稳定输出基准电压,电源DVCC通过并联的第一滤波电容C1和第二滤波电容C2接地DGND。
该参考电压电路以稳压二极管U1的第一引脚作为输出端输出等效直流侧电压的参考电压Vref,输出的参考电压Vref的大小取决于R2、R3以及稳压二极管的基准电压,满足以下关系:
式中,VT为所用稳压管的基准电压。
根据设定的电压范围可以通过调节相关参数得到不同的参考电压。
本实施例中参考电压电路中稳压二极管U1为TL431型,基准电压为2.5V,本实施例R2=1.2kΩ,R3=3kΩ,其余参数根据常规选择设定,分别为C1=0.1uF,C2=0.01uF,C3=1nF,R1=120Ω,相应的输出的参考电压为:
图3为迟滞比较电路,包括比较器U2,等效直流侧电压Vdc通过电阻R4连接到比较器U2的正向输入端5(V+),基准电压Vref通过电阻R5连接到比较器U2的反向输入端6(V-)。正向输入端5还通过电容C4接地DGND,起滤波作用,比较器U2的输出端7通过电阻R7连接到电源DVCC(5V),起上拉作用,且正向输入端与输出端之间连接有反馈电阻R6,比较器的输出端连接到下一级驱动保护电路。
本实施例中设定的电压范围为5V以下,相应的第三阈值V1和第四阈值V2分别为4.4V和3.1V,需要合理设置迟滞比较电路中的参数,设定原理如下:
其中,VpH和VpL分别为滞环电路输出电压的高低电平。
本实施例中迟滞比较电路中VpH=0V,VpL=5V,相应的设定R4=6.2kΩ,R6=24kΩ。其余参数根据常规选择设定,分别如下:
C4=1nF,R5=5.1kΩ,R7=39kΩ。
正常情况下,输入的等效直流侧电压在3.1V~4.4V之间,此时输出的电平信号Vp为0V。当输入的等效直流侧电压大于或等于4.4V(第四阈值V2)时,输出的电平信号Vp为高电平,直至等效直流侧电压小于或等于3.1V(第三阈值V1)时,输出的电平信号Vp重新为低电平。
为了进一步验证本发明的有益效果,在搭建的一台14kVA三相四线制SAPF样机上进行了相应的实验验证,样机***结构如图4所示。图中,uSa、uSb、uSc为电网三相相电压,iSa、iSb、iSc、iSN分别为电网三相相电流和中线电流,Lga、Lgb、Lgc为电网各相漏感,iLa、iLb、iLc、iLN分别为负载三相相电流和中线电流,iCa、iCb、iCc、iCN分别为补偿三相相电流和中线电流。非线性负载采用三相不控整流桥带阻感负载,Lc、Ls、Rd和Cf组成交流侧输出滤波器,滤除输出电流中的开关纹波及不需要的高频分量。三相四线制SAPF样机***的主要参数如表1所示。
表1
符号 | 参数说明 | 数值 |
us | 电网相电压 | 220V |
fs | 电网频率 | 50Hz |
Vdc | 直流侧电压 | 750V |
Cdc | 直流侧总电容 | 11.75mF |
fsm | 采样频率 | 10kHz |
fsw | 开关频率 | 10kHz |
图5是该14kVA三相四线制SAPF样机在额定满载条件下的全频段稳态补偿波形,图中IL、IC、IS分别为负载电流、补偿电流和网侧电流。额定补偿电流有效值ip=20A。可以看出,额定工作条件下,基于指令电流限流的复合限流方法不起作用,指令电流值即为检测出的负载谐波电流值,补偿后的电网侧电流波形呈正弦波,总谐波失真THD(Total Harmonic Distortion)=3.6%,说明该样机具有良好的谐波补偿能力。
图6是在该样机上进行限流补偿的稳态实验结果,图中IL、IC、IS分别为负载电流、补偿电流和网侧电流。设定负载谐波电流有效值为60A,且选择性谐波补偿频次为3、5、7、11次。则需要补偿的最高谐波次数h=11,相应的指令电流限流中定义***所允许的最大补偿电流有效值(即第一阈值)imax=30A,峰值(即第二阈值)ipeak=69A。由图可见,尽管负载电流中待补偿的谐波电流有效值是该样机额定补偿能力的3倍,远远超过了它的补偿范围,但***补偿电流的有效值始终限制在imax,不会出现严重过流运行的情况。另外,利用wavestar对补偿电流进行频谱分析,可知其有效值为29.4A,很接近第一阈值imax=30A,误差在2%以内,具有很高的控制精度。同时也可以看出,相比于定义***最大补偿电流有效值始终取额定值的传统方法,在保证SAPF安全可靠运行的前提下多输出补偿了接近10A的谐波电流,补偿低次谐波的情况下装置利用效率提高了50%。
图7是在该样机上进行限流补偿的动态实验结果,图中IL、IC、IS分别为负载电流、补偿电流和网侧电流。设定空载-40A谐波负载动态切换,SAPF的补偿频次仍旧为3、5、7、11次。可以看出,突加负载时,虽然待补偿的负载谐波电流远远超过了该SAPF的补偿能力,但是样机输出的实际补偿电流有效值被很好的限制在了30A以内,没有出现负载突变时刻补偿电流超出***最大允许值而使得SAPF过流运行的情况,说明负载动态切换时本限流策略也能具有很好的控制效果。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。
Claims (8)
1.一种有源电力滤波器的保护方法,其特征在于,包括如下步骤:
(1)提取负载电流中需要补偿的各次谐波分量,并将提取的各次谐波分量叠加得到指令电流;
(2)根据有源电力滤波器的额定补偿电流有效值和需要补偿的最高谐波次数依次设定第一阈值和第二阈值;
(3)将第一阈值作为最大允许电流对所述的指令电流进行比例限流,得到第一指令电流;
(4)将第二阈值作为截断峰值电流对所述的第一指令电流进行截断限流,得到第二指令电流;
(5)根据第二指令电流生成用于驱动有源电力滤波器中变流器的驱动脉冲,使变流器输出相应的补偿电流。
2.如权利要求1所述的有源电力滤波器的保护方法,其特征在于,所述步骤(1)基于滑窗迭代DFT算法提取需要补偿的各次谐波分量。
3.如权利要求1所述的有源电力滤波器的保护方法,其特征在于,所述步骤(2)根据如下公式计算第一阈值imax:
imax=Ki*ip,
其中,ip为额定补偿电流有效值;
Ki为第一比例系数:
其中,h为需要补偿的最高谐波次数。
4.如权利要求1所述的有源电力滤波器的保护方法,其特征在于,所述步骤(2)根据如下公式计算第二阈值ipeak:
ipeak=Kp*imax,
其中,imax为第一阈值;
Kp为第二比例系数:
其中,h为需要补偿的最高谐波次数。
5.如权利要求1所述的有源电力滤波器的保护方法,其特征在于,所述步骤(3)过程如下:
(3-1)计算指令电流的有效值,并用第一阈值除以指令电流的有效值和第一阈值中的较大者得到的商作为比例限流因子;
(3-2)根据公式:
i'ref=k*iref
计算第一指令电流i'ref,iref为指令电流,k为比例限流因子。
6.如权利要求1所述的有源电力滤波器的保护方法,其特征在于,所述的保护方法还包括实时检测有源电力滤波器的直流侧电压,根据直流侧电压对有源电力滤波器进行硬件封锁限流,具体过程如下:
(S1)实时检测有源电力滤波器的直流侧电压;
(S2)将检测到的直流侧电压调理到设定的电压范围内得到相应的等效直流侧电压;
(S3)将调理得到的等效直流侧电压与第三阈值进行比较,并根据比较结果进行如下操作:
若大于或等于第三阈值,则封锁有源电力滤波器中变流器的驱动脉冲,直至等效直流侧电压下降至第四阈值时解除封锁,使变流器输出相应的补偿电流;
否则,则不操作。
7.如权利要求6所述的有源电力滤波器的保护方法,其特征在于,所述的第三阈值V1根据公式:
V1=KT*(Vdcp+ΔVh)
计算,Vdcp表示直流侧电压额定值,ΔVh表示直流侧过压保护的过压裕量,KT表示等效调理变比。
8.如权利要求7所述的有源电力滤波器的保护方法,其特征在于,所述的第四阈值V2根据公式:
V2=KT*(Vdcp-ΔVl)
计算,Vdcp表示直流侧电压额定值,ΔVl表示直流侧过压保护解除的欠压裕量,KT表示等效调理变比。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410305638.6A CN104113052B (zh) | 2014-06-30 | 2014-06-30 | 有源电力滤波器的保护方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410305638.6A CN104113052B (zh) | 2014-06-30 | 2014-06-30 | 有源电力滤波器的保护方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104113052A true CN104113052A (zh) | 2014-10-22 |
CN104113052B CN104113052B (zh) | 2017-05-17 |
Family
ID=51709739
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410305638.6A Active CN104113052B (zh) | 2014-06-30 | 2014-06-30 | 有源电力滤波器的保护方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104113052B (zh) |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104600682A (zh) * | 2015-02-05 | 2015-05-06 | 湖北工业大学 | 有源电力滤波器电流保护电路及方法 |
CN104749435A (zh) * | 2015-04-01 | 2015-07-01 | 四川厚天科技股份有限公司 | 一种无累积误差滑窗的dft谐波检测方法 |
CN105048466A (zh) * | 2015-09-10 | 2015-11-11 | 苏州市欧米伽智能控制技术有限公司 | 一种并联型混合有源电力滤波器的复合控制方法 |
CN105207217A (zh) * | 2015-09-09 | 2015-12-30 | 泰州学院 | 有源电力滤波器瞬时值截断比例限流控制方法 |
CN106253279A (zh) * | 2016-08-23 | 2016-12-21 | 南京亚派科技股份有限公司 | 一种应用于有源滤波器的防谐振控制算法 |
CN106684874A (zh) * | 2017-03-15 | 2017-05-17 | 河南理工大学 | 一种有源电力滤波器特定次谐波优化补偿方法 |
CN109038580A (zh) * | 2018-07-20 | 2018-12-18 | 乐山拉得电网自动化有限公司 | 一种解决负载突变引起滑动dft异常而导致APF故障的方法 |
CN109066630A (zh) * | 2018-09-21 | 2018-12-21 | 上海沫安能源科技有限公司 | 一种有源电力滤波器的限流装置及其限流方法 |
CN109546866A (zh) * | 2018-11-30 | 2019-03-29 | 广东美的制冷设备有限公司 | 空调***、变频控制器及其交流电流估算方法和装置 |
CN109592550A (zh) * | 2019-01-24 | 2019-04-09 | 通力电梯有限公司 | 控制柜及包括其的自动扶梯 |
CN110829432A (zh) * | 2019-10-14 | 2020-02-21 | 北京精密机电控制设备研究所 | 一种apf的限幅方法 |
CN111525584A (zh) * | 2020-07-06 | 2020-08-11 | 深圳华工能源技术有限公司 | 基于全方位电能质量优化治理的配用电***智能节电方法 |
CN111815880A (zh) * | 2020-08-07 | 2020-10-23 | 成都惠之慧安全科技有限公司 | 一种电气火灾及燃气智能监控报警*** |
CN112018731A (zh) * | 2020-07-08 | 2020-12-01 | 中车株洲电力机车研究所有限公司 | 整流滤波设备的监控方法、设备和*** |
CN113991672A (zh) * | 2021-10-25 | 2022-01-28 | 国网山东省电力公司青岛供电公司 | 一种基于三电平双向变流的分布式控制***及方法 |
CN114400642A (zh) * | 2022-01-25 | 2022-04-26 | 帝森克罗德集团有限公司 | 一种有源电力滤波器的谐波提取及限流方法 |
CN117595266A (zh) * | 2024-01-18 | 2024-02-23 | 深圳市海思瑞科电气技术有限公司 | 基于apf装置的谐波治理方法、装置、设备及存储介质 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050024131A1 (en) * | 2003-07-28 | 2005-02-03 | James Lazar | Current control method and apparatus for active power filters |
CN102208803A (zh) * | 2011-06-03 | 2011-10-05 | 上海理工大学 | 并联型多目标短路电流故障限流器 |
-
2014
- 2014-06-30 CN CN201410305638.6A patent/CN104113052B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050024131A1 (en) * | 2003-07-28 | 2005-02-03 | James Lazar | Current control method and apparatus for active power filters |
CN102208803A (zh) * | 2011-06-03 | 2011-10-05 | 上海理工大学 | 并联型多目标短路电流故障限流器 |
Non-Patent Citations (1)
Title |
---|
张国荣等: "有源电力滤波器的并联运行及其控制策略", 《低压电器》 * |
Cited By (24)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104600682B (zh) * | 2015-02-05 | 2017-05-24 | 湖北工业大学 | 有源电力滤波器电流保护电路及方法 |
CN104600682A (zh) * | 2015-02-05 | 2015-05-06 | 湖北工业大学 | 有源电力滤波器电流保护电路及方法 |
CN104749435A (zh) * | 2015-04-01 | 2015-07-01 | 四川厚天科技股份有限公司 | 一种无累积误差滑窗的dft谐波检测方法 |
CN105207217A (zh) * | 2015-09-09 | 2015-12-30 | 泰州学院 | 有源电力滤波器瞬时值截断比例限流控制方法 |
CN105048466A (zh) * | 2015-09-10 | 2015-11-11 | 苏州市欧米伽智能控制技术有限公司 | 一种并联型混合有源电力滤波器的复合控制方法 |
CN106253279A (zh) * | 2016-08-23 | 2016-12-21 | 南京亚派科技股份有限公司 | 一种应用于有源滤波器的防谐振控制算法 |
CN106684874A (zh) * | 2017-03-15 | 2017-05-17 | 河南理工大学 | 一种有源电力滤波器特定次谐波优化补偿方法 |
CN106684874B (zh) * | 2017-03-15 | 2019-12-27 | 河南理工大学 | 一种有源电力滤波器特定次谐波优化补偿方法 |
CN109038580B (zh) * | 2018-07-20 | 2021-11-26 | 乐山一拉得电网自动化有限公司 | 一种解决负载突变引起滑动dft异常而导致APF故障的方法 |
CN109038580A (zh) * | 2018-07-20 | 2018-12-18 | 乐山拉得电网自动化有限公司 | 一种解决负载突变引起滑动dft异常而导致APF故障的方法 |
CN109066630A (zh) * | 2018-09-21 | 2018-12-21 | 上海沫安能源科技有限公司 | 一种有源电力滤波器的限流装置及其限流方法 |
CN109546866A (zh) * | 2018-11-30 | 2019-03-29 | 广东美的制冷设备有限公司 | 空调***、变频控制器及其交流电流估算方法和装置 |
CN109592550A (zh) * | 2019-01-24 | 2019-04-09 | 通力电梯有限公司 | 控制柜及包括其的自动扶梯 |
CN110829432A (zh) * | 2019-10-14 | 2020-02-21 | 北京精密机电控制设备研究所 | 一种apf的限幅方法 |
CN111525584A (zh) * | 2020-07-06 | 2020-08-11 | 深圳华工能源技术有限公司 | 基于全方位电能质量优化治理的配用电***智能节电方法 |
CN112018731A (zh) * | 2020-07-08 | 2020-12-01 | 中车株洲电力机车研究所有限公司 | 整流滤波设备的监控方法、设备和*** |
CN112018731B (zh) * | 2020-07-08 | 2022-08-16 | 中车株洲电力机车研究所有限公司 | 整流滤波设备的监控方法、设备和*** |
CN111815880A (zh) * | 2020-08-07 | 2020-10-23 | 成都惠之慧安全科技有限公司 | 一种电气火灾及燃气智能监控报警*** |
CN113991672A (zh) * | 2021-10-25 | 2022-01-28 | 国网山东省电力公司青岛供电公司 | 一种基于三电平双向变流的分布式控制***及方法 |
CN113991672B (zh) * | 2021-10-25 | 2024-05-31 | 国网山东省电力公司青岛供电公司 | 一种基于三电平双向变流的分布式控制***及方法 |
CN114400642A (zh) * | 2022-01-25 | 2022-04-26 | 帝森克罗德集团有限公司 | 一种有源电力滤波器的谐波提取及限流方法 |
CN114400642B (zh) * | 2022-01-25 | 2023-11-10 | 帝森克罗德集团有限公司 | 一种有源电力滤波器的谐波提取及限流方法 |
CN117595266A (zh) * | 2024-01-18 | 2024-02-23 | 深圳市海思瑞科电气技术有限公司 | 基于apf装置的谐波治理方法、装置、设备及存储介质 |
CN117595266B (zh) * | 2024-01-18 | 2024-03-26 | 深圳市海思瑞科电气技术有限公司 | 基于apf装置的谐波治理方法、装置、设备及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN104113052B (zh) | 2017-05-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104113052B (zh) | 有源电力滤波器的保护方法 | |
CN103368181B (zh) | Lcl型滤波器的三相并网逆变器的电流谐波抑制方法 | |
CN105207186B (zh) | 一种用于含统一潮流控制器的输电线路的距离保护方法 | |
CN101793918A (zh) | 一种电压暂降检测方法 | |
CN110954743B (zh) | 一种分布式录波装置及小电流接地选线方法 | |
CN106950512B (zh) | 一种储能变流器并离网特性一体化检测***及方法 | |
CN103701303A (zh) | 具有可测性功能的开关电源及测试方法 | |
CN108627731A (zh) | 一种单相断电的快速检测方法 | |
CN109100600B (zh) | 一种磁控式并联电抗器故障确定方法及*** | |
Zhang et al. | Application of partial direct-pole-placement and differential evolution algorithm to optimize controller and LCL filter design for grid-tied inverter | |
Hosseini | The operation and model of UPQC in voltage sag mitigation using EMTP by direct method | |
CN111650422B (zh) | 高压直流***的同步触发方法、***及存储介质 | |
Tanti et al. | An ANN based approach for optimal placement of DSTATCOM for voltage sag mitigation | |
Bhasker et al. | Wavelet transform based discrimination between inrush and internal fault of indirect symmetrical phase shift transformer | |
CN106410769A (zh) | 一种基于参考滤波电感电流的下垂控制微电网的限流方法 | |
CN109066630B (zh) | 一种有源电力滤波器的限流装置及其限流方法 | |
Tomy et al. | Power quality improvement strategy for non-linear load in single phase system | |
Abdul-Hameed et al. | An investigation of the harmonic effects of nonlinear loads on power distributon network | |
Hosseinpour et al. | RDFT based control of a 3-phase 4-leg shunt active power filter under unbalanced load and non-ideal mains voltage conditions | |
Raut et al. | An overview and design of Dynamic Voltage Restorer to improve power quality in microgrid | |
Sahril et al. | A single phase dynamic voltage restorer (DVR) with direct AC-AC converter using dq transform to mitigate voltage sag | |
Ismail et al. | Modelling voltage sag mitigation using dynamic voltage restorer and analyzing power quality issue | |
CN108767821B (zh) | 一种选相合闸投切方法和*** | |
Chan et al. | Boundary controller for dynamic voltage restorers to achieve fast dynamic response | |
Cuma et al. | Implementation of a non-linear adaptive filter based sag detection method for dynamic voltage restorers under unbalanced fault conditions |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |