CN106226587A - 一种基于les‑‑hht的交流微电网电压暂降快速检测方法 - Google Patents

一种基于les‑‑hht的交流微电网电压暂降快速检测方法 Download PDF

Info

Publication number
CN106226587A
CN106226587A CN201610524871.2A CN201610524871A CN106226587A CN 106226587 A CN106226587 A CN 106226587A CN 201610524871 A CN201610524871 A CN 201610524871A CN 106226587 A CN106226587 A CN 106226587A
Authority
CN
China
Prior art keywords
signal
voltage
hht
les
voltage dip
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
Application number
CN201610524871.2A
Other languages
English (en)
Other versions
CN106226587B (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.)
Zhejiang University of Technology ZJUT
Original Assignee
Zhejiang University of Technology ZJUT
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 Zhejiang University of Technology ZJUT filed Critical Zhejiang University of Technology ZJUT
Priority to CN201610524871.2A priority Critical patent/CN106226587B/zh
Publication of CN106226587A publication Critical patent/CN106226587A/zh
Application granted granted Critical
Publication of CN106226587B publication Critical patent/CN106226587B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R19/00Arrangements for measuring currents or voltages or for indicating presence or sign thereof
    • G01R19/25Arrangements for measuring currents or voltages or for indicating presence or sign thereof using digital measurement techniques
    • G01R19/2503Arrangements for measuring currents or voltages or for indicating presence or sign thereof using digital measurement techniques for measuring voltage only, e.g. digital volt meters (DVM's)

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Measurement Of Resistance Or Impedance (AREA)

Abstract

一种基于LES‑‑HHT的交流微电网电压暂降快速检测方法,包括如下步骤:a.原始电压暂降信号采集;b.以步骤a中得到的信号为输入,通过HHT变换后得到希尔伯特频谱,进而得到瞬时频率;c.以步骤b变换后的信号为输入,结合连续电压信号N次非均匀采样的规律,得到新的采样频率,通过自适应采样得到采样数据;d.提取电压暂降信号的幅值和相位角:通过改进型LES算法提取电压暂降信号各基波分量和谐波分量的幅值和相位角。本发明通过自适应采样解决采样精度和检测速度矛盾的问题,通过改进型LES算法减少计算量,实现对电压跌落信号的快速检测。

Description

一种基于LES--HHT的交流微电网电压暂降快速检测方法
技术领域
本发明涉及交流微电网电压信号检测分析技术领域,尤其涉及电压暂降快速检测方法。
背景技术
微电网可将多种间歇式新能源和负荷联系起来,对电网为可控的单元,具有灵活便捷就地消纳等优点,因此微电网是分布式新能源的重要组织利用形式。但由于微电网存在大量间歇式电源,并采用电力电子装置接入主馈线,惯性几乎为零,抗扰动能力弱,同时由于微电网的容量有限,DG的波动、负荷的波动、离并网过程都会对馈线产生较大的影响。因此,相对于大电网,微电网交流母线的电压闪变和暂降更加频繁剧烈。电压暂降会产生瞬时电流突变,导致微网中敏感用电设备工作失常。因此,一种能快速检测电压暂降的算法是必要的,由此可以设计出相应的对策,减少电压暂降造成的危害。
传统的电压暂降检测算法主要有以下几种:快速傅里叶变换(Fast FourierTransform,FFT)法、小波变换法、瞬时电压d-q变换法、对称分量预估法、卡尔曼滤波法、希尔伯特黄变换(Hilbert-Huang Transform,HHT)法和最小方差(Least Error Squares,LES)滤波算法等。目前,瞬时电压d-q变换法在电压跌落检测中应用最广泛。但是,若电压波形存在谐波或畸变时,d-q变换后需要加入滤波器,d-q变换算法的响应速度就会受影响。FFT法存在一个工频周期的延时,因此检测误差仍较大。小波变换法和卡尔曼滤波法运算比较复杂,计算量过大,也不利于电压跌落的快速检测。HHT法利用经验模态分解(EmpiricalMode Decomposition,EMD)方法将信号分解为固有模态函数(Intrinsic Mode Function,IMF),再将其作Hilbert变换,结合瞬时频率得到电压的幅值和相位。但HHT算法存在端点效应、包络线和均值曲线拟合等问题,造成了较大的瞬时幅值计算偏差,降低检测精度。LES滤波法通过提取被测电压的幅值和相位,获取电压跌落的特征值。结合对称分量法可达到快速检测的目的,但其在电压跌落瞬间的计算波动较大,影响检测精度。可见,这些分析方法都不能满足电压暂降快速检测的要求。
发明内容
为了克服已有电压暂降检测方式的无法兼顾采样速度和检测精度的不足,本发明提供一种兼顾采样速度和检测精度的基于LES--HHT的交流微电网电压暂降快速检测方法,通过自适应采样解决采样精度和检测速度矛盾的问题,通过改进型LES算法减少计算量,实现对电压跌落信号的快速检测。
本发明解决其技术问题所采用的技术方案是:
一种基于LES--HHT的交流微电网电压暂降快速检测方法,所述检测方法包括如下步骤:
a.原始电压暂降信号采集
通过电压互感器采集交流测电压暂降信号作为输入,采样得到电压暂降信号;
b.提取瞬时频率
以步骤a中得到的信号为输入,通过HHT变换后得到希尔伯特频谱,进而得到瞬时频率;
c.自适应采样
以步骤b变换后的信号为输入,并结合连续电压信号N次非均匀采样的规律,得到新的采样频率,通过自适应采样得到采样数据。
d.提取电压暂降信号的幅值和相位角
以步骤c得到的采样数据为输入,通过改进型LES算法提取电压暂降信号各基波分量和谐波分量的幅值和相位角,过程如下:
通过改进型LES算法提取到的电压信号的幅值和相位角分别为:θa=arctan(X1/X2),其中,X1=kasin(θa),X2=kacos(θa),式中,ka为电压峰值,θa为初始角度;
设:
U=[u(t) u(t-Δt)…u(t-(N-1)Δt))]T
A = cos ( Nω 0 Δ t ) cos ( ( N - 1 ) ω 0 Δ t ) ... cos ( ω 0 Δ t ) sin ( Nω 0 Δ t ) sin ( ( N - 1 ) ω 0 Δ t ) ... sin ( ω 0 Δ t ) T
X=[kasin(θa) kacos(θa)]T=[X1 X2]T
其中,U为对电压信号进行N次采样,采样时间间隔为Δt,u(t)表示电压,A、X是为了便于计算自定义的矩阵,根据输入采样数据U,由U=A×X即可得到矩阵X的值,上述矩阵运算结果得出:其中,是[ATA]-1的分块矩阵表示,K表示采样的次序,对于每次新的采样来说,将计算过程中不变的量D11、D12、D21、D22、cos((N+i)ω0Δt)、cos(iω0Δt)、sin((N+i)ω0Δt)、sin(iω0Δt)分别存入DSP内存,即可实现电压暂降的快速检测。
进一步,在步骤b中,HHT变换的具体过程为,电压暂降信号通过经验模态分解法EMD分解,得到电压跌落信号的本征模态函数IMF,再经Hilbert变换后得到希尔伯特频谱。电压暂降信号经HHT变换最终得到的希尔伯特频谱为:式中,ai(t)为瞬时相位函数,fi(t)为每个IMF分量的瞬时频率。
再进一步,在步骤c中,自适应采样时间间隔Δt由步骤b得到的瞬时频率来决定,其关系可表达为:Δt=Γ-1((k+1)Ts)-Γ-1(kTs)≈Ts/fi(t),其中,对连续电压信号X(t)=f[Γ(t)]的N次非均匀采样的结果可表示为式中,Γ(t)为伽马函数,δ(t)为单位脉冲序列,Ts为一常数。
本发明的有益效果主要表现在:通过自适应采样,有效解决采样精度与检测速度矛盾的问题;通过使用改进型LES算法,可减少***计算时间,有效提高电压跌落的检测速度。
附图说明
图1是基于LES--HHT的交流微电网电压暂降快速检测方法的流程图。
具体实施方式
下面结合附图对本发明作进一步描述。
参照图1,一种基于LES--HHT的交流微电网电压暂降快速检测方法,所述检测方法包括如下步骤:
a.原始电压暂降信号采集
通过电压互感器采集交流测电压暂降信号作为输入,采样得到电压暂降信号;
b.提取瞬时频率
以步骤a中得到的信号为输入,通过HHT变换后得到希尔伯特频谱,进而得到瞬时频率;
c.自适应采样
以步骤b变换后的信号为输入,并结合连续电压信号N次非均匀采样的规律,得到新的采样频率,通过自适应采样得到采样数据。
d.提取电压暂降信号的幅值和相位角
以步骤c得到的采样数据为输入,通过改进型LES算法提取电压暂降信号各基波分量和谐波分量的幅值和相位角,过程如下:
通过改进型LES算法提取到的电压信号的幅值和相位角分别为:θa=arctan(X1/X2),其中,X1=kasin(θa),X2=kacos(θa),式中,ka为电压峰值,θa为初始角度;
设:
U=[u(t) u(t-Δt)…u(t-(N-1)Δt))]T
A = cos ( Nω 0 Δ t ) cos ( ( N - 1 ) ω 0 Δ t ) ... cos ( ω 0 Δ t ) sin ( Nω 0 Δ t ) sin ( ( N - 1 ) ω 0 Δ t ) ... sin ( ω 0 Δ t ) T
X=[kasin(θa) kacos(θa)]T=[X1 X2]T
其中,U为对电压信号进行N次采样,采样时间间隔为Δt,u(t)表示电压,A、X是为了便于计算自定义的矩阵,根据输入采样数据U,由U=A×X即可得到矩阵X的值,上述矩阵运算结果得出:其中,是[ATA]-1的分块矩阵表示,K表示采样的次序,对于每次新的采样来说,将计算过程中不变的量D11、D12、D21、D22、cos((N+i)ω0Δt)、cos(iω0Δt)、sin((N+i)ω0Δt)、sin(iω0Δt)分别存入DSP内存,即可实现电压暂降的快速检测。
进一步,在步骤b中,HHT变换的具体过程为,电压暂降信号通过经验模态分解法EMD分解,得到电压跌落信号的本征模态函数IMF,再经HHT变换后得到希尔伯特频谱。电压暂降信号经HHT变换最终得到的希尔伯特频谱为:式中,ai(t)为瞬时相位函数,fi(t)为每个IMF分量的瞬时频率。
再进一步,在步骤c中,自适应采样时间间隔Δt由步骤b得到的瞬时频率来决定,其关系可表达为:Δt=Γ-1((k+1)Ts)-Γ-1(kTs)≈Ts/fi(t),其中,对连续电压信号X(t)=f[Γ(t)]的N次非均匀采样的结果可表示为式中,Γ(t)为伽马函数,δ(t)为单位脉冲序列,Ts为一常数。
本实施例中,从微电网交流母线或其他关键支路上采集原始电压暂降信号作为待处理信号输入,通过自适应采样和改进型LES算法快速计算出暂降电压的幅值和相位,从而为电压暂降补偿提供可行的策略,保证微电网***的稳定运行。
本实例以20kVA的动态电压恢复器(DVR)测试平台为例分析。基于LES--HHT算法的交流微电网电压暂降快速检测方法包括如下步骤:
步骤1,原始电压暂降信号采集
本实例通过电压互感器在交流侧采集电压信号,互感器的输出信号经过运放调理后输入到数字信号处理器(DSP)的A/D采样口,初始的采样频率是10kHz。采样得到的信号X(t)通过经验模态分解法(EMD)分解,得到电压跌落信号的本征模态函数(IMF),可表示为经Hilbert变换后得到的希尔伯特频谱,每个IMF的瞬时频率为由公式Δt≈Ts/fi(t)来决定下一个周期的采样频率。
步骤2,通过改进型LES算法快速计算出电压暂降信号的幅值和相位。
以步骤1得到的新的采样频率对电压信号进行采样,采样结果用矩阵表示,为:U=[u(t)u(t-Δt)…u(t-(N-1)Δt)]T,定义如下矩阵:
A = cos ( Nω 0 Δ t ) cos ( ( N - 1 ) ω 0 Δ t ) ... cos ( ω 0 Δ t ) sin ( Nω 0 Δ t ) sin ( ( N - 1 ) ω 0 Δ t ) ... sin ( ω 0 Δ t ) T
X=[kasin(θa) kacos(θa)]T=[X1 X2]T
其关系为U=A×X。由矩阵运算可改写为X=[ATA]-1ATU,令D=[ATA]-1,并用分块矩阵表示 取中间变量矩阵X1=D11×E1+D12×E2,X2=D21×E1+D22×E2,计算E矩阵产生的变化量:
E1new=E1+u(t)cos((N+i)ω0Δt)-u(t-NΔt)cos(iω0Δt)
E2new=E2+u(t)sin((N+i)ω0Δt)-u(t-NΔt)sin(iω0Δt)
并将其中计算过程中不变的量D11、D12、D21、D22、cos((N+i)ω0Δt)、cos(iω0Δt)、sin((N+i)ω0Δt)、sin(iω0Δt)事先存入DSP内存。最后,暂降电压信号的幅值为瞬时相位为θa=arctan(X1/X2)。
最后,还需要注意的是,以上列举的仅是本发明的一个具体实施例。显然,本发明不限于以上实施例,还可以有许多变形。本领域的普通技术人员能从本发明公开的内容直接导出或联想到的所有变形,均应认为是本发明的保护范围。

Claims (3)

1.一种基于LES--HHT的交流微电网电压暂降快速检测方法,其特征在于:所述检测方法包括如下步骤:
a.原始电压暂降信号采集
通过电压互感器采集交流测电压暂降信号作为输入,采样得到电压暂降信号;
b.提取瞬时频率
以步骤a中得到的信号为输入,通过HHT变换后得到希尔伯特频谱,进而得到瞬时频率;
c.自适应采样
以步骤b变换后的信号为输入,并结合连续电压信号N次非均匀采样的规律,得到新的采样频率,通过自适应采样得到采样数据。
d.提取电压暂降信号的幅值和相位角
以步骤c得到的采样数据为输入,通过改进型LES算法提取电压暂降信号各基波分量和谐波分量的幅值和相位角,过程如下:
通过改进型LES算法提取到的电压信号的幅值和相位角分别为:θa=arctan(X1/X2),其中,X1=kasin(θa),X2=kacos(θa),式中,ka为电压峰值,θa为初始角度;
设:
U=[u(t)u(t-Δt)…u(t-(N-1)Δt))]T
A = cos ( Nω 0 Δ t ) cos ( ( N - 1 ) ω 0 Δ t ) ... cos ( ω 0 Δ t ) sin ( Nω 0 Δ t ) sin ( ( N - 1 ) ω 0 Δ t ) ... sin ( ω 0 Δ t ) T
X=[kasin(θa) kacos(θa)]T=[X1 X2]T
其中,U为对电压信号进行N次采样,采样时间间隔为Δt,u(t)表示电压,A、X是为了便于计算自定义的矩阵,根据输入采样数据U,由U=A×X即可得到矩阵X的值,上述矩阵运算结果得出:
其中,是[ATA]-1的分块矩阵表示,K表示采样的次序,对于每次新的采样来说,将计算过程中不变的量D11、D12、D21、D22、cos((N+i)ω0Δt)、cos(iω0Δt)、sin((N+i)ω0Δt)、sin(iω0Δt)分别存入DSP内存,即可实现电压暂降的快速检测。
2.如权利要求1所述的一种基于LES--HHT的交流微电网电压暂降快速检测方法,其特征在于,在步骤b中,HHT变换的具体过程为,电压暂降信号通过经验模态分解法EMD分解,得到电压跌落信号的本征模态函数IMF,再经Hilbert变换后得到希尔伯特频谱。电压暂降信号经HHT变换最终得到的希尔伯特频谱为:式中,ai(t)为瞬时相位函数,fi(t)为每个IMF分量的瞬时频率。
3.如权利要求1或2所述的一种基于LES--HHT的交流微电网电压暂降快速检测方法,其特征在于,在步骤c中,自适应采样时间间隔Δt由步骤b得到的瞬时频率来决定,其关系可表达为:Δt=Γ-1((k+1)Ts)-Γ-1(kTs)≈Ts/fi(t),其中,对连续电压信号X(t)=f[Γ(t)]的N次非均匀采样的结果可表示为式中,Γ(t)为伽马函数,δ(t)为单位脉冲序列,Ts为一常数。
CN201610524871.2A 2016-07-01 2016-07-01 一种基于les--hht的交流微电网电压暂降快速检测方法 Active CN106226587B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610524871.2A CN106226587B (zh) 2016-07-01 2016-07-01 一种基于les--hht的交流微电网电压暂降快速检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610524871.2A CN106226587B (zh) 2016-07-01 2016-07-01 一种基于les--hht的交流微电网电压暂降快速检测方法

Publications (2)

Publication Number Publication Date
CN106226587A true CN106226587A (zh) 2016-12-14
CN106226587B CN106226587B (zh) 2019-07-05

Family

ID=57519602

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610524871.2A Active CN106226587B (zh) 2016-07-01 2016-07-01 一种基于les--hht的交流微电网电压暂降快速检测方法

Country Status (1)

Country Link
CN (1) CN106226587B (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109164294A (zh) * 2018-09-06 2019-01-08 上海电力学院 一种用于dvr中的自适应电压跌落检测方法
CN109307795A (zh) * 2018-11-30 2019-02-05 广东电网有限责任公司 波形点检测方法及装置
CN110363130A (zh) * 2019-07-08 2019-10-22 国网四川省电力公司电力科学研究院 基于变分模态分解的电压暂降源辨识方法及辨识装置
CN113311223A (zh) * 2021-06-03 2021-08-27 上海市计量测试技术研究院 一种电压跌落检测仪及检测方法
CN113514686A (zh) * 2021-07-13 2021-10-19 北京英博电气股份有限公司 电压基波幅值的检测方法、装置、设备和存储介质
CN114325076A (zh) * 2021-12-13 2022-04-12 国网北京市电力公司 电压暂降的检测方法、检测装置和处理器

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102831433A (zh) * 2012-06-06 2012-12-19 西南交通大学 一种基于时频域多特征量的电能质量混合扰动分类方法
CN103424600A (zh) * 2013-08-20 2013-12-04 昆明理工大学 一种基于Hilbert-Huang变换和小波包能量谱的电压暂降源识别方法
US20160011268A1 (en) * 2014-07-10 2016-01-14 National Cheng Kung University Motor fault detecting method and motor fault detecting system
CN105510687A (zh) * 2015-12-24 2016-04-20 合肥工业大学 一种基于经验模态分解的电压异常特征识别方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102831433A (zh) * 2012-06-06 2012-12-19 西南交通大学 一种基于时频域多特征量的电能质量混合扰动分类方法
CN103424600A (zh) * 2013-08-20 2013-12-04 昆明理工大学 一种基于Hilbert-Huang变换和小波包能量谱的电压暂降源识别方法
US20160011268A1 (en) * 2014-07-10 2016-01-14 National Cheng Kung University Motor fault detecting method and motor fault detecting system
CN105510687A (zh) * 2015-12-24 2016-04-20 合肥工业大学 一种基于经验模态分解的电压异常特征识别方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
MARYAM MOEINDARBARI等: "Design of an automated system for detection and classification of power quality disturbances", 《2015 5TH INTERNATIONAL CONFERENCE ON COMPUTER AND KNOWLEDGE ENGINEERING (ICCKE)》 *
于婷: "电压暂降的检测与定位方法研究", 《中国优秀硕士学位论文全文数据库工程科技II辑》 *
刘云潺等: "HHT在电压暂降检测中的应用", 《科技创新导报》 *

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109164294A (zh) * 2018-09-06 2019-01-08 上海电力学院 一种用于dvr中的自适应电压跌落检测方法
CN109307795A (zh) * 2018-11-30 2019-02-05 广东电网有限责任公司 波形点检测方法及装置
CN109307795B (zh) * 2018-11-30 2020-11-06 广东电网有限责任公司 波形点检测方法及装置
CN110363130A (zh) * 2019-07-08 2019-10-22 国网四川省电力公司电力科学研究院 基于变分模态分解的电压暂降源辨识方法及辨识装置
CN110363130B (zh) * 2019-07-08 2023-01-13 国网四川省电力公司电力科学研究院 基于变分模态分解的电压暂降源辨识方法及辨识装置
CN113311223A (zh) * 2021-06-03 2021-08-27 上海市计量测试技术研究院 一种电压跌落检测仪及检测方法
CN113514686A (zh) * 2021-07-13 2021-10-19 北京英博电气股份有限公司 电压基波幅值的检测方法、装置、设备和存储介质
CN113514686B (zh) * 2021-07-13 2024-02-06 北京英博电气股份有限公司 电压基波幅值的检测方法、装置、设备和存储介质
CN114325076A (zh) * 2021-12-13 2022-04-12 国网北京市电力公司 电压暂降的检测方法、检测装置和处理器
CN114325076B (zh) * 2021-12-13 2023-10-24 国网北京市电力公司 电压暂降的检测方法、检测装置和处理器

Also Published As

Publication number Publication date
CN106226587B (zh) 2019-07-05

Similar Documents

Publication Publication Date Title
CN106226587A (zh) 一种基于les‑‑hht的交流微电网电压暂降快速检测方法
CN101907437B (zh) 一种基于小波差分算法的电缆故障测距方法
CN106019076B (zh) 一种高压直流输电线路故障测距方法
CN103345162B (zh) 一种功率级数模混合仿真***
CN105319447B (zh) 一种介质损耗角正切值测试方法及测试仪
CN103050942B (zh) 一种基于hht的ct饱和检测方法
CN103675544B (zh) 基于优化算法的电力***故障信号检测与波形识别方法
CN102305891A (zh) 一种电力***低频振荡在线监测方法
CN104714075B (zh) 一种电网电压闪变包络参数提取方法
CN106980044B (zh) 一种适应风电接入的电力***谐波电流估计方法
CN107478896A (zh) 一种基于级联广义积分器的频率自适应谐波电流检测方法
CN105487034A (zh) 一种0.05级电子式互感器校验方法及***
CN103529294A (zh) 基于hht的光伏***并网逆变器谐波检测***及方法
CN102508031A (zh) 一种基于傅里叶级数的局部放电脉冲相角测量方法
CN103904693A (zh) 基于频率自适应虚拟磁链估测的电网同步方法
CN101329374A (zh) 一种差分滤波器加全周复小波幅值计算方法
CN106250904A (zh) 基于改进s变换的电能扰动分析仪及分类方法
CN102636693A (zh) 一种结合fft与非线性最小二乘的谐波分析算法
CN103901273A (zh) 电网谐波检测方法和电网谐波检测装置
CN103675447A (zh) 一种电气化铁路的高精度实时谐波分析方法
ElRefaie et al. A novel technique to eliminate the effect of decaying DC component on DFT based phasor estimation
Damala et al. A simple decision tree-based disturbance monitoring system for VSC-based HVDC transmission link integrating a DFIG wind farm
CN103091561A (zh) 从交直流叠加信号中提取直流信号的装置及方法
CN101609112B (zh) 基于广域测量噪声信号的电力***节点相位关系辨识方法
CN101393237A (zh) 一种三相在线谐波电流实时监测***及其工作方法

Legal Events

Date Code Title Description
C06 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