CN113608023A - 采样频率自适应的谐波检测方法及系统 - Google Patents
采样频率自适应的谐波检测方法及系统 Download PDFInfo
- Publication number
- CN113608023A CN113608023A CN202110949594.0A CN202110949594A CN113608023A CN 113608023 A CN113608023 A CN 113608023A CN 202110949594 A CN202110949594 A CN 202110949594A CN 113608023 A CN113608023 A CN 113608023A
- Authority
- CN
- China
- Prior art keywords
- points
- data
- frequency
- fft
- harmonic detection
- 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
- 238000001514 detection method Methods 0.000 title claims abstract description 28
- 238000005070 sampling Methods 0.000 title claims abstract description 22
- 238000012952 Resampling Methods 0.000 claims abstract description 23
- 238000000034 method Methods 0.000 claims abstract description 17
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 16
- 238000004364 calculation method Methods 0.000 claims description 9
- 230000003044 adaptive effect Effects 0.000 claims description 5
- 239000000872 buffer Substances 0.000 claims description 4
- 238000004590 computer program Methods 0.000 claims description 4
- 238000000605 extraction Methods 0.000 claims description 3
- 230000006870 function Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 238000012544 monitoring process Methods 0.000 description 5
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000000630 rising effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000003139 buffering effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000010891 electric arc Methods 0.000 description 1
- 238000012806 monitoring device Methods 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 238000003908 quality control method Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
- G01R23/16—Spectrum analysis; Fourier analysis
Landscapes
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Current Or Voltage (AREA)
Abstract
本发明涉及智能电网技术领域,公开一种采样频率自适应的谐波检测方法及系统,以匹配IEC规定并提升整体的自适应性。本发明方法包括:利用DSP的AD转换器,采集电压和电流信号,根据IEC谐波检测标准,分别缓存10周期数据;针对缓存的数据点数N0,依据2m或2m的整数倍,确定要重采样的目标点数N;根据N0和N的关系,确定相应的重采样算法将N0点数据重采至N点;将重采样后的N个数据进行FFT运算,利用DSP的FFT函数快速求得,得到频率间隔为5Hz的各频率分量的实部和虚部结果,并求得各频率分量的幅值信息;根据IEC标准中的子组算法计算各次谐波含量,50Hz的整数倍谐波检测的相位作为该次谐波的相位。
Description
技术领域
本发明涉及智能电网技术领域,尤其涉及一种采样频率自适应的谐波检测方法及系统。
背景技术
电能质量领域包括电能质量治理和电能质量监测,对于电能质量治理装置而言,具备电能质量监测功能也必不可少。由于台区中存在开关器件和相关的设备,比如公变台区有电视机、空调等设备,专变台区有电弧炉等,因此线路中避免不了存在谐波电流。对于电能质量监测装置而言,设备主要具备监测功能,可以使用高速采样频率和高采样位数,保证检测精度。对于电能质量治理装置而言,在产品设计时主要以电能质量的治理功能为主,一些采样频率或者控制频率要考虑到治理功能。但是电能质量监测功能同样重要,在电能质量监测的参数中,谐波由于其频率高、变化复杂以及检测运算复杂度高等因素,成为比较难实现的检测参数。
现有的行业内电能质量治理设备对于谐波检测的方法以SDFT和SOGI等算法为主,存在与IEC规定的方法不兼容的技术问题。
发明内容
本发明目的在于公开一种采样频率自适应的谐波检测方法及系统,以匹配IEC规定并提升整体的自适应性。
为达上述目的,本发明公开一种采样频率自适应的谐波检测方法,包括以下步骤:
步骤1:利用DSP的AD转换器,采集电压和电流信号,并缓存至数组里,根据IEC谐波检测标准,分别缓存10周期数据;
步骤2:针对缓存的数据点数N0,依据2m或2m的整数倍,确定要重采样的目标点数N;
步骤3:根据N0和N的关系,确定相应的重采样算法将N0点数据重采至N点;
步骤4:将重采样后的N个数据进行FFT运算,利用DSP的FFT函数快速求得,得到频率间隔为5Hz的各频率分量的实部和虚部结果,并求得各频率分量的幅值信息;
步骤5:根据IEC标准中的子组算法计算各次谐波含量,50Hz的整数倍谐波检测的相位作为该次谐波的相位。
优选地,当N0小于N时,所述步骤4具体包括:通过均匀抽点的方式,做至少两个低点数的FFT,再根据所需检测的频率,提取对应频率的低点数FFT计算结果,根据抽取的顺序依次乘以对应的系数,再将这几个乘系数后的结果相加得到等效为高点数FFT得到的该频率的结果。
优选地,当重采样点数不为2m时,对缓存数据进行分包处理,使得每一小包数据量为2m,然后对每一小包进行低点数FFT计算。
可选地,所述步骤1具体包括:对电压电流信号每个周期采样102.4个点,使得每10个周期采样1024个点。进一步地,所述步骤1以等间隔对缓存数据时进行缓存。
为达上述目的,本发明还公开一种采样频率自适应的谐波检测系统,包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现上述方法的步骤。
本发明具有以下有益效果:
能结合内存使用量以及谐波检测需求的矛盾因素,对于一些谐波检测次数不是很高的工况下,缓存数据时可以隔几个点进行缓存,这样既能满足检测需求,又能提高内存效率。
考虑到采样频率不确定性,因此N选择最接近于N0的2m的整数倍,而N与N0之间的关系不确定,因此重采样功能包含升点数重采样和降点数重采样。
为保证FFT运算速度、FFT的运算精度以及内存使用量的平衡,或在重采样点数不为2m时,可以对缓存数据进行分包处理,每一小包数据量为2m,因此每一小包进行FFT运算,对运算后的数据进行处理,可以得到高点数重采样数据的FFT的运算结果。
下面将参照附图,对本发明作进一步详细的说明。
附图说明
构成本申请的一部分的附图用来提供对本发明的进一步理解,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。在附图中:
图1是发明实施例中一种采样频率自适应的谐波检测方法流程图。
图2是发明实施例中从高点数到低点数的重采样示意图。
图3是发明实施例中从低点数到高点数的重采样示意图。
图4是发明实施例中从低点数FFT转换至高点数FFT的算法流程图。
图5是发明实施例中谐波检测的子组算法示意图。
具体实施方式
以下结合附图对本发明的实施例进行详细说明,但是本发明可以由权利要求限定和覆盖的多种不同方式实施。
实施例1
本实施例公开一种采样频率自适应的谐波检测方法,如图1所示,包括以下步骤:
步骤1:利用DSP的AD转换器,采集电压和电流信号,并缓存至数组里,根据IEC谐波检测标准,分别缓存10周期数据。
步骤2:针对缓存的数据点数N0,依据2m或2m的整数倍,确定要重采样的目标点数N。
步骤3:根据N0和N的关系,对信号进行不同方案的重采样算法,将N0点数据重采至N点,且保持误差较小。
步骤4:将重采样后的N个数据进行FFT运算,利用DSP的FFT函数快速求得,得到频率间隔为5Hz的各频率分量的实部和虚部结果,并求得各频率分量的幅值信息。
步骤5:根据IEC标准中的子组算法计算各次谐波含量,50Hz的整数倍谐波检测的相位作为该次谐波的相位。
重采样算法示意图见图2和图3,图2是从高点数到低点数的降点数重采样,图3是从低点数到高点数的升点数重采样,二者得到重采样数据的方式都是通过插值的方法,考虑到算法实现的复杂性以及信号的随机性,因此插值采用简单的一阶线性插值方法,降点数重采样方法和升点数重采样方法的区别在于,如图2所示,对于降点数重采样方法而言,n≤k,重采样信号的个数逐次递增,但是所用到的Xk会有不连续的情况;如图3所示,对于升点数重采样方法而言,n≥k,这就使得会出现Yn和Yn+1都在Xk和Xk+1之间的情况,因此两种重采样方案要分别进行实现。
考虑到DSP的实现以及重采样算法对于精度的影响,DSP芯片对于FFT运算所需的缓存和点数的局限性,并且谐波检测的最高次数一般为50次(即2.5kHz左右),因此根据奈奎斯特采样定理,对电压电流信号每个周期采样102.4个点,即每10个周期采样1024个点就能保证FFT分析至50次谐波。因此,如图4所示,对于一些采样点数为2m的整数倍或者较高的2m个数的信号,可以通过均匀抽点的方式,做多个低点数的FFT,再根据所需检测的频率,提取对应频率的低点数FFT计算结果,根据抽取的顺序依次乘以对应的系数,再将这几个乘系数后的结果相加,即等效为高点数FFT得到的该频率的结果。
图5是子组算法示意图,以基波检测为例,10个基波周期信号的FFT检测得到的是以5Hz为频率间隔的频谱信息,因此不仅仅可以得到50Hz的结果,也可以得到45Hz和55Hz的结果,在计算基波成分时,联合45Hz、50Hz和55Hz检测结果共同计算有效值,即
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (6)
1.一种采样频率自适应的谐波检测方法,其特征在于,包括以下步骤:
步骤1:利用DSP的AD转换器,采集电压和电流信号,并缓存至数组里,根据IEC谐波检测标准,分别缓存10周期数据;
步骤2:针对缓存的数据点数N0,依据2m或2m的整数倍,确定要重采样的目标点数N;
步骤3:根据N0和N的关系,确定相应的重采样算法将N0点数据重采至N点;
步骤4:将重采样后的N个数据进行FFT运算,利用DSP的FFT函数快速求得,得到频率间隔为5Hz的各频率分量的实部和虚部结果,并求得各频率分量的幅值信息;
步骤5:根据IEC标准中的子组算法计算各次谐波含量,50Hz的整数倍谐波检测的相位作为该次谐波的相位。
2.根据权利要求1所述的方法,其特征在于,当N0小于N时,所述步骤4具体包括:
通过均匀抽点的方式,做至少两个低点数的FFT,再根据所需检测的频率,提取对应频率的低点数FFT计算结果,根据抽取的顺序依次乘以对应的系数,再将这几个乘系数后的结果相加得到等效为高点数FFT得到的该频率的结果。
3.根据权利要求2所述的方法,其特征在于,还包括:
当重采样点数不为2m时,对缓存数据进行分包处理,使得每一小包数据量为2m,然后对每一小包进行低点数FFT计算。
4.根据权利要求1至3任一所述的方法,其特征在于,所述步骤1具体包括:
对电压电流信号每个周期采样102.4个点,使得每10个周期采样1024个点。
5.根据权利要求4所述的方法,其特征在于,所述步骤1还包括:
以等间隔对缓存数据时进行缓存。
6.一种采样频率自适应的谐波检测系统,包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现上述权利要求1至5任一所述方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110949594.0A CN113608023A (zh) | 2021-08-18 | 2021-08-18 | 采样频率自适应的谐波检测方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110949594.0A CN113608023A (zh) | 2021-08-18 | 2021-08-18 | 采样频率自适应的谐波检测方法及系统 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113608023A true CN113608023A (zh) | 2021-11-05 |
Family
ID=78308921
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110949594.0A Pending CN113608023A (zh) | 2021-08-18 | 2021-08-18 | 采样频率自适应的谐波检测方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113608023A (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114624513A (zh) * | 2022-01-27 | 2022-06-14 | 清华大学 | 周期信号的抗谐波干扰的相位检测方法及装置 |
CN117169590A (zh) * | 2023-08-15 | 2023-12-05 | 嘉兴市科讯电子有限公司 | 一种基于软件变采样率的电力谐波分析的方法和装置 |
CN117574030A (zh) * | 2023-11-17 | 2024-02-20 | 青岛艾诺仪器有限公司 | 动态采样的谐波分量计算系统及其控制方法 |
CN118091252A (zh) * | 2024-04-22 | 2024-05-28 | 青岛艾诺仪器有限公司 | 一种基于fpga平台的谐波计算方法及计算系统 |
Citations (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH095362A (ja) * | 1995-06-16 | 1997-01-10 | Koyo Denki Kk | 波形検出方法およびその装置 |
JP2006209892A (ja) * | 2005-01-28 | 2006-08-10 | Sanyo Electric Co Ltd | Pll回路およびディスク再生装置 |
CN101915874A (zh) * | 2010-07-20 | 2010-12-15 | 北海市深蓝科技发展有限责任公司 | 一种基于傅立叶变换的谐波检测方法 |
CN103353558A (zh) * | 2013-05-31 | 2013-10-16 | 深圳市康必达控制技术有限公司 | 一种电能质量监测方法 |
CN103439590A (zh) * | 2013-08-08 | 2013-12-11 | 深圳市康必达中创科技有限公司 | 一种电能监测dsp控制系统 |
JP2014081352A (ja) * | 2012-09-27 | 2014-05-08 | Daihen Corp | 周波数解析装置、当該周波数解析装置を用いた信号処理装置、および、当該信号処理装置を用いた高周波測定装置 |
CN103995181A (zh) * | 2014-05-13 | 2014-08-20 | 武汉中元华电科技股份有限公司 | 一种数字化变电站电能质量谐波分析方法 |
CN104076203A (zh) * | 2014-07-07 | 2014-10-01 | 国家电网公司 | 一种考虑负频率影响的超低频间谐波检测方法 |
CN104639085A (zh) * | 2015-02-06 | 2015-05-20 | 国家电网公司 | 一种微弱电信号低功耗滤波电路 |
CN105718423A (zh) * | 2016-01-19 | 2016-06-29 | 清华大学 | 一种流水可重构的单精度浮点fft/ifft协处理器 |
CN106405229A (zh) * | 2016-08-30 | 2017-02-15 | 威胜集团有限公司 | 一种基波和谐波电能计量方法 |
CN106896360A (zh) * | 2017-04-21 | 2017-06-27 | 南京航空航天大学 | 一种sar信号处理算法的fpga实现方法 |
CN107085144A (zh) * | 2017-04-28 | 2017-08-22 | 珠海泰芯半导体有限公司 | 一种快速测量电力系统谐波的方法 |
CN109088617A (zh) * | 2018-09-20 | 2018-12-25 | 电子科技大学 | 比率可变数字重采样滤波器 |
CN109507480A (zh) * | 2018-12-12 | 2019-03-22 | 中国电力科学研究院有限公司 | 一种邻近基波/谐波的间谐波检测方法和装置 |
CN109633262A (zh) * | 2019-01-29 | 2019-04-16 | 国网湖南省电力有限公司 | 基于组合窗多谱线fft的三相谐波电能计量方法、装置 |
CN109829132A (zh) * | 2019-01-21 | 2019-05-31 | 东南大学 | 一种嵌入式环境下长数据序列快速谱分析方法 |
CN110596635A (zh) * | 2019-07-26 | 2019-12-20 | 广西电网有限责任公司 | 一种电能表后故障检测的方法及电能表 |
CN112180161A (zh) * | 2020-08-11 | 2021-01-05 | 国网山西省电力公司电力科学研究院 | 一种非同步高采样率采样条件下谐波间谐波群测量方法 |
US20210141854A1 (en) * | 2019-04-28 | 2021-05-13 | Harbin Institute Of Technology | Fourier Analysis Method with Variable Sampling Frequency |
CN113157637A (zh) * | 2021-04-27 | 2021-07-23 | 电子科技大学 | 一种基于fpga的大容量可重构的fft运算ip核 |
-
2021
- 2021-08-18 CN CN202110949594.0A patent/CN113608023A/zh active Pending
Patent Citations (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH095362A (ja) * | 1995-06-16 | 1997-01-10 | Koyo Denki Kk | 波形検出方法およびその装置 |
JP2006209892A (ja) * | 2005-01-28 | 2006-08-10 | Sanyo Electric Co Ltd | Pll回路およびディスク再生装置 |
CN101915874A (zh) * | 2010-07-20 | 2010-12-15 | 北海市深蓝科技发展有限责任公司 | 一种基于傅立叶变换的谐波检测方法 |
JP2014081352A (ja) * | 2012-09-27 | 2014-05-08 | Daihen Corp | 周波数解析装置、当該周波数解析装置を用いた信号処理装置、および、当該信号処理装置を用いた高周波測定装置 |
CN103353558A (zh) * | 2013-05-31 | 2013-10-16 | 深圳市康必达控制技术有限公司 | 一种电能质量监测方法 |
CN103439590A (zh) * | 2013-08-08 | 2013-12-11 | 深圳市康必达中创科技有限公司 | 一种电能监测dsp控制系统 |
CN103995181A (zh) * | 2014-05-13 | 2014-08-20 | 武汉中元华电科技股份有限公司 | 一种数字化变电站电能质量谐波分析方法 |
CN104076203A (zh) * | 2014-07-07 | 2014-10-01 | 国家电网公司 | 一种考虑负频率影响的超低频间谐波检测方法 |
CN104639085A (zh) * | 2015-02-06 | 2015-05-20 | 国家电网公司 | 一种微弱电信号低功耗滤波电路 |
CN105718423A (zh) * | 2016-01-19 | 2016-06-29 | 清华大学 | 一种流水可重构的单精度浮点fft/ifft协处理器 |
CN106405229A (zh) * | 2016-08-30 | 2017-02-15 | 威胜集团有限公司 | 一种基波和谐波电能计量方法 |
CN106896360A (zh) * | 2017-04-21 | 2017-06-27 | 南京航空航天大学 | 一种sar信号处理算法的fpga实现方法 |
CN107085144A (zh) * | 2017-04-28 | 2017-08-22 | 珠海泰芯半导体有限公司 | 一种快速测量电力系统谐波的方法 |
CN109088617A (zh) * | 2018-09-20 | 2018-12-25 | 电子科技大学 | 比率可变数字重采样滤波器 |
CN109507480A (zh) * | 2018-12-12 | 2019-03-22 | 中国电力科学研究院有限公司 | 一种邻近基波/谐波的间谐波检测方法和装置 |
CN109829132A (zh) * | 2019-01-21 | 2019-05-31 | 东南大学 | 一种嵌入式环境下长数据序列快速谱分析方法 |
CN109633262A (zh) * | 2019-01-29 | 2019-04-16 | 国网湖南省电力有限公司 | 基于组合窗多谱线fft的三相谐波电能计量方法、装置 |
US20210141854A1 (en) * | 2019-04-28 | 2021-05-13 | Harbin Institute Of Technology | Fourier Analysis Method with Variable Sampling Frequency |
CN110596635A (zh) * | 2019-07-26 | 2019-12-20 | 广西电网有限责任公司 | 一种电能表后故障检测的方法及电能表 |
CN112180161A (zh) * | 2020-08-11 | 2021-01-05 | 国网山西省电力公司电力科学研究院 | 一种非同步高采样率采样条件下谐波间谐波群测量方法 |
CN113157637A (zh) * | 2021-04-27 | 2021-07-23 | 电子科技大学 | 一种基于fpga的大容量可重构的fft运算ip核 |
Non-Patent Citations (1)
Title |
---|
朱明星: "一种多窗口宽度的主导间谐波频谱分布的算法", 电力电容器与无功补偿, 5 February 2018 (2018-02-05), pages 1 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114624513A (zh) * | 2022-01-27 | 2022-06-14 | 清华大学 | 周期信号的抗谐波干扰的相位检测方法及装置 |
CN117169590A (zh) * | 2023-08-15 | 2023-12-05 | 嘉兴市科讯电子有限公司 | 一种基于软件变采样率的电力谐波分析的方法和装置 |
CN117169590B (zh) * | 2023-08-15 | 2024-03-29 | 嘉兴市科讯电子有限公司 | 一种基于软件变采样率的电力谐波分析的方法和装置 |
CN117574030A (zh) * | 2023-11-17 | 2024-02-20 | 青岛艾诺仪器有限公司 | 动态采样的谐波分量计算系统及其控制方法 |
CN118091252A (zh) * | 2024-04-22 | 2024-05-28 | 青岛艾诺仪器有限公司 | 一种基于fpga平台的谐波计算方法及计算系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113608023A (zh) | 采样频率自适应的谐波检测方法及系统 | |
CN107085144B (zh) | 一种快速测量电力系统谐波的方法 | |
CN109462404B (zh) | 基于相似度分段的自适应波形数据压缩方法 | |
CN111521914B (zh) | 一种用于确定高压输电直流线路起晕场强的方法及系统 | |
CN104459354B (zh) | 一种三相交流电网相序检测方法及装置 | |
CN109655644B (zh) | 一种降低任意波信号输出抖动的方法和装置 | |
CN115372698A (zh) | 一种用于电力系统高次谐波抑制的测控装置及方法 | |
CN110618338A (zh) | 一种电能质量瞬态测量方法、装置和设备 | |
CN104111373A (zh) | 一种用于智能变电站数字化电能的计量方法 | |
CN113484596A (zh) | 一种电能质量监测方法、装置、设备及可读存储介质 | |
CN115864449B (zh) | 新能源接入系统振荡监测方法、装置、系统、设备及介质 | |
CN206132884U (zh) | 一种电能质量监测装置 | |
CN113406453B (zh) | 一种基于mcu的prpd/prps图谱数据处理方法及检测装置 | |
CN109782064B (zh) | 一种风电场输出阻抗频率特性测试分析方法 | |
Qiu et al. | Lossless compression of synchro-waveform measurements for smart grid monitoring | |
CN111400649A (zh) | 谐波责任量化方法、装置、计算机设备和存储介质 | |
An et al. | Supraharmonics Measurement Based on Hybrid Online Measurement and Offline Analysis | |
CN114184838A (zh) | 基于sn互卷积窗的电力系统谐波检测方法、系统及介质 | |
CN110244144A (zh) | 一种避雷器状态监测方法及交流数据采集方法 | |
CN116520022B (zh) | 电力谐波动态检测方法、装置、电子设备及介质 | |
CN103604990A (zh) | 一种电力谐波检测分析设备 | |
CN108318761A (zh) | 基于压缩感知的风力发电机组电能质量检测方法 | |
Qiu et al. | Design of A Hybrid Compression Algorithm for High-fidelity Synchro-waveform Measurements | |
CN115575706B (zh) | 适用于电能表校验仪的频率连续跟踪测量系统和方法 | |
CN112836390B (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 |