CN116451553A - 一种改进的变分模态分解与BiGRU融合的水质预测方法 - Google Patents
一种改进的变分模态分解与BiGRU融合的水质预测方法 Download PDFInfo
- Publication number
- CN116451553A CN116451553A CN202211398569.9A CN202211398569A CN116451553A CN 116451553 A CN116451553 A CN 116451553A CN 202211398569 A CN202211398569 A CN 202211398569A CN 116451553 A CN116451553 A CN 116451553A
- Authority
- CN
- China
- Prior art keywords
- water quality
- modal
- data
- decomposition
- representing
- 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
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 50
- 238000000354 decomposition reaction Methods 0.000 title claims abstract description 43
- 238000000034 method Methods 0.000 title claims abstract description 25
- 230000004927 fusion Effects 0.000 title claims abstract description 13
- 241000287127 Passeridae Species 0.000 claims abstract description 18
- 238000010845 search algorithm Methods 0.000 claims abstract description 11
- 230000002159 abnormal effect Effects 0.000 claims abstract description 9
- 238000010606 normalization Methods 0.000 claims abstract description 8
- 238000007781 pre-processing Methods 0.000 claims abstract description 8
- 238000012545 processing Methods 0.000 claims abstract description 7
- 230000000694 effects Effects 0.000 claims description 8
- 238000005457 optimization Methods 0.000 claims description 4
- 238000010276 construction Methods 0.000 claims description 3
- 229940050561 matrix product Drugs 0.000 claims description 3
- 238000001228 spectrum Methods 0.000 claims description 3
- 238000012549 training Methods 0.000 abstract description 4
- 230000006870 function Effects 0.000 description 4
- 230000002457 bidirectional effect Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 229910052760 oxygen Inorganic materials 0.000 description 2
- 239000001301 oxygen Substances 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000001133 acceleration Effects 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 description 1
- 239000010931 gold Substances 0.000 description 1
- 229910052737 gold Inorganic materials 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000003062 neural network model Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000004083 survival effect Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000003911 water pollution Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/27—Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
-
- 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
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A20/00—Water conservation; Efficient water supply; Efficient water use
- Y02A20/152—Water filtration
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Artificial Intelligence (AREA)
- Computer Hardware Design (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种改进的变分模态分解与BiGRU融合的水质预测方法,包括:收集水质历史数据,对原数据进行数据预处理,填补缺失值,剔除异常值,最后数据进行归一化处理,构建水质数据集;利用麻雀搜索算法对变分模态分解的参数进行寻优;改进后的变分模态分解算法对水质数据进行分解,分解为若干个模态分量,以降低数据复杂度;将每个模态分量输出到BiGRU预测模型中训练,并将每个模态预测结果进行叠加,得到最终预测结果。本发明对水质数据进行分解,能够很好解决水质数据波动性大、非线性的问题,提高预测准确性,为水质预测提供一种新的方法。
Description
技术领域
本发明涉及水质预测技术领域,具体涉及一种改进的变分模态分解与BiGRU融合的水质预测方法。
背景技术
水不仅是人类生存不可缺少的自然资源,也是国家经济发展的重要战略资源。近几年,随着工业化和城市进程化速度的加快,水污染事件经常发生,水资源短缺问题日益严重,水质预测作为水资源管理的前置环节,精确预测水质指标能及时掌握水质变化趋势,避免污染事件发生。水质数据是波动性大、非线性的时间序列,并且含有大量噪声信号,给水质预测带来难度,将原始数据分解成几个简单序列是一种分析数据的重要手段,常用分解方法有小波变换、经验模态分解(EMD)、集成经验模式分解(EEMD)等。VMD作为一种新的信号分解方法,能够克服EMD和EEMD存在问题,但是VMD分解层数和惩罚因子需要人为选取,选择的大小会影响VMD的分解效果。目前,水质预测应用比较多的模型有灰色系统理论模型、时间序列模型以及神经网络模型等,GRU网络作为LSTM的变体,同样能解决梯度下降问题,相比较LSTM结构更加简单,参数更少,但是单向GRU在训练过程中只能从前向后提取数据的特征,无法充分体现时间序列双向的相关性。
发明内容
本发明的目的是为了解决上述存在问题,利用麻雀搜索算法对VMD的两个参数进行寻优,再利用VMD对水质数据分解,有效克服水质数据非线性、非稳定性的问题。另外采用双向GRU作为预测模型,从正反两个反向提取水质数据特征,提高预测精度。
为实现以上目的,本发明通过以下技术方案予以实现:一种改进的变分模态分解与BiGRU融合的水质预测方法,包括以下步骤:
(1)收集水质历史数据,对原数据进行数据预处理,填补缺失值,剔除异常值,最后对数据进行归一化处理,构建水质数据;
(2)利用麻雀搜索算法对变分模态分解的参数进行寻优,使VMD分解算法达到最佳分解效果;
(3)改进后的变分模态分解对水质数据进行分解,分解为若干个模态分量,以降低数据复杂度;
(4)将每个模态分量输出到BiGRU预测模型中,并将每个模态预测结果进行叠加,得到最终预测结果。
所述的数据预处理具体包括:采用向前插值方式填充缺失值,采用孤立森林算法识别异常值,最后对数据归一化处理。
所述的麻雀搜索算法对变分模态分解的两个参数分解个数k和惩罚因子α进行寻优,具体步骤如下:
步骤1:参数初始化,根据时间序列大小,设置k和α合适范围;
步骤2:麻雀种群初始化,设置优化维数、迭代次数、预警值、种群大小等参数;
步骤3:选取模糊熵作为适应度函数,计算各模态分量模糊熵值,经过多次迭代,当模糊熵值最小时,此时麻雀种群位置为最优位置;
步骤4:确定k和α最优值。
所述的VMD将原数据f分解为k个中心频率为ωk的有限带宽的模态分量uk,具体步骤如下:
(1)构建变分问题:首先利用希尔伯特变换得到uk的解析信号,计算得到单边谱,再将每个模态解析信号与中心频率混合,调制到相应的基带。然后计算上述信号梯度L平方的范数,得到变分约束模型。
式中:uk(k=1,2,…)为模态分量;ωk表示各个模态分量中心频率;表示时间的偏导数;
(2)求解变分问题:首先引入拉格朗日乘法算子λ和惩罚因子α,将约束模型变为无约束的变分模型,然后采用交替方向乘子迭代算法,得到各模态分量和中心频率,uk和ωk表达式如下:
式中:f(ω)、λn(t)分别为f(t),/>λn(t)的傅里叶变换;n表示迭代次数。
如图3所示,BiGRU模型由一个正向和一个反向GRU构成,BiGRU预测模型算法表示为:
式中:zt表示重置门;rt表示更新门;Wr,Wz,Ur,Uz分别表示权重矩阵;ht-1表示上一级隐藏层状态;xt表示当前时刻输入向量;ht表示当前时刻候选的隐藏状态;ht表示当前时刻隐藏状态;“*”表示矩阵乘积;σ表示sigmoid函数。
所述的预测结果是将各个模态分量输入到建立BiGRU预测模型中,然后将每个模态分量预测值叠加,得到最终的预测结果。
本发明具有以下有益效果:
(1)本发明针对水质数据非稳定性和非线性等特点,利用VMD分解降低水质复杂程度,VMD相比EMD和EEMD相比有效防止模态混叠和虚假分量的现象,采用麻雀搜索算法对VMD参数寻优,使VMD分解效果达到最佳。
(2)相比于LSTM网络,GRU结构更加简单,参数更少,保留LSTM预测效果,并且收敛效果好。而采用双向BiGRU能够从正反两个方向充分挖掘水质深层特征,输出水质数据完整的特征,使得预测精度更高。
附图说明
图1为本发明VMD与BiGRU融合水质预测方法的流程图
图2为本发明麻雀搜索算法优化VMD的流程图
图3为本发明BiGRU网络结构图
图4为本发明SSA优化VMD的迭代图
图5为本发明VMD分解结果图
图6为本发明VMD-BiGRU预测结果图
具体实施方式
下面结合附图及实施例对本发明做进一步详细说明,本发明提出一种改进的变分模态分解与BiGRU融合的水质预测方法,如图1所示,具体步骤如下:
(1)收集水质历史数据,对原数据进行数据预处理,填补缺失值,剔除异常值,最后对数据进行归一化处理,构建水质数据;
(2)利用麻雀搜索算法对变分模态分解的参数进行寻优,使VMD分解算法达到最佳分解效果;
(3)改进后的变分模态分解对水质数据进行分解,分解为若干个模态分量,以降低数据复杂度;
(4)将每个模态分量输出到BiGRU预测模型中,并将每个模态预测结果进行叠加,得到最终预测结果;
所述的数据预处理具体包括:采用向前插值方式填充缺失值,采用孤立森林算法识别异常值,最后对数据归一化处理。
如图2所示,麻雀搜索算法对变分模态分解的两个参数分解个数k和惩罚因子α进行寻优,具体步骤如下:
步骤1:参数初始化,根据时间序列大小,设置k和α合适范围;
步骤2:麻雀种群初始化,设置优化维数、迭代次数、预警值、种群大小等参数;
步骤3:选取模糊熵作为适应度函数,计算各模态分量模糊熵值,经过多次迭代,当模糊熵值最小时,此时麻雀种群位置为最优位置;
步骤4:确定k和α最优值。
所述的VMD将原数据f分解为k个中心频率为ωk的有限带宽的模态分量uk,具体步骤如下:
(1)构建变分问题:首先利用希尔伯特变换得到uk的解析信号,计算得到单边谱,再将每个模态解析信号与中心频率混合,调制到相应的基带。然后计算上述信号梯度L平方的范数,得到变分约束模型。
式中:uk(k=1,2,…)为模态分量;ωk表示各个模态分量中心频率;表示时间的偏导数;
(2)求解变分问题:首先引入拉格朗日乘法算子λ和惩罚因子α,将约束模型变为无约束的变分模型,然后采用交替方向乘子迭代算法,得到各模态分量和中心频率,uk和ωk表达式如下:
式中:f(ω)、λn(t)分别为f(t),/>λn(t)的傅里叶变换;n表示迭代次数。
如图3所示,BiGRU模型由一个正向和一个反向GRU构成,BiGRU预测模型算法表示为:
式中:zt表示重置门;rt表示更新门;Wr,Wz,Ur,Uz分别表示权重矩阵;ht-1表示上一级隐藏层状态;xt表示当前时刻输入向量;ht表示当前时刻候选的隐藏状态;ht表示当前时刻隐藏状态;“*”表示矩阵乘积;σ表示sigmoid函数。
实例分析1
步骤1:为了验证本发明的有效性,收集水质历史数据,选取浙闽片河流域小金山断面2021年7至2021年12月溶解氧数据为例,构建VMD-BiGRU预测模型。
步骤2:数据预处理,首先采用向前插值方式填补缺失值,再采用孤立森林算法对异常值进行初步识别,去除异常数据后采用插值法修正,最后采用最大值最小值归一化方法。
步骤3:利用麻雀搜索算法对VMD参数寻优在寻优过程中,当适应度值达最小或者达到最大迭代次数时,输出的值为参数最优组合,如图4所示为麻雀算法迭代过程,迭代第6次时,适应度值达到最小,值为0.017782。
步骤4:麻雀搜索算法优化VMD后,得到VMD的分解个数为9,惩罚因子为2425,利用VMD对处理后数据进行分解,如图5所示为溶解氧分解结果。
步骤5:按照7:3比例将每个模态分量划分为测试集和训练集,分别输入到BiGRU预测模型中,BiGRU网络由一层BiGRU层、Dropout层、全连接层组成,设置epochs值为100,batch_size值为20,学习率为0.001。
步骤6:最后将每个模态预测结果相加得到最终预测结果,如图6所示为VMD与BiGRU融合预测结果。
步骤7:选取均方根误差(RMSE)、平均绝对百分(MAPE)、平均绝对误差(MAE)、决定系数(R2)作为评价指标,来分析融合模型预测效果,如表1所示,相比较单一GRU模型,融合模型溶解氧预测误差更小,预测精度更高。
表1模型评估表
模型 | RMSE | MAPE/% | MAE | R2 |
GRU | 0.1497 | 6.7572 | 0.0909 | 0.8876 |
VMD-BiGRU | 0.0764 | 0.07071 | 0.0586 | 0.9748 |
综上所述,本实施例提出一种改进的变分模态分解与BiGRU融合的水质预测方法,该方法在BiGRU预测模型的基础上融合VMD分解,有效解决水质非线性、非稳定性的问题,使预测精确度得到提升,预测误差减小,能够很好预测水质变化趋势,有利于管理人员掌握水质变化情况。
以上所述仅为本发明较佳的实施方式,并不用于限制本发明。本领域普通技术人员应当理解,在不脱离本发明的精神和原则的情况下,可以对实施例进行修改、替换、变型,本发明的范围由所附权利要求及其等同物限。
Claims (6)
1.一种改进的变分模态分解与BiGRU融合的水质预测方法,其特征在于,包括以下步骤:
(1)收集水质历史数据,对原数据进行数据预处理,填补缺失值,剔除异常值,最后数据进行归一化处理,构建水质数据;
(2)利用麻雀搜索算法对变分模态分解的参数进行寻优,使VMD分解算法达到最佳分解效果;
(3)改进后的变分模态分解对水质数据进行分解,分解为若干个模态分量,以降低数据复杂度;
(4)将每个模态分量输出到BiGRU预测模型中,并将每个模态预测结果进行叠加,得到最终预测结果。
2.根据权利要求1所述的一种改进的变分模态分解与BiGRU融合的水质预测方法,其特征在于,数据预处理具体包括:采用向前插值方式填充缺失值,采用孤立森林算法识别异常值,最后对数据归一化处理。
3.根据权利要求1所述的一种改进的变分模态分解与BiGRU融合的水质预测方法,其特征在于,利用麻雀搜索算法对变分模态分解的两个参数分解个数k和惩罚因子α进行寻优具体步骤如下:
步骤1:参数初始化,根据时间序列大小,设置k和α合适范围;
步骤2:麻雀种群初始化,设置优化维数、迭代次数、预警值、种群大小等参数;
步骤3:选取模糊熵作为适应度函数,计算各模态分量模糊熵值,经过多次迭代,当模糊熵值最小时,此时麻雀种群位置为最优位置;
步骤4:确定k和α最优值。
4.根据权利要求1所述的一种改进的变分模态分解与BiGRU融合的水质预测方法,其特征在于,VMD将原数据f分解为k个中心频率为ωk的有限带宽的模态分量uk,具体步骤如下:
(1)构建变分问题:首先利用希尔伯特变换得到uk的解析信号,计算得到单边谱,再将每个模态解析信号与中心频率混合,调制到相应的基带,然后计算上述信号梯度L平方的范数,得到变分约束模型;
式中:uk(k=1,2,…)为模态分量;ωk表示各个模态分量中心频率;表示时间的偏导数;
(2)求解变分问题:首先引入拉格朗日乘法算子λ和惩罚因子α,将约束模型变为无约束的变分模型,然后采用交替方向乘子迭代算法,得到各模态分量和中心频率,uk和ωk表达式如下:
式中:f(ω)、λn(t)分别为f(t),/>λn(t)的傅里叶变换;n表示迭代次数。
5.根据权利要求1所述的一种改进的变分模态分解与BiGRU融合的水质预测方法,其特征在于,将每个模态分量输出到BiGRU预测模型中,BiGRU由一个正向和一个反向GRU构成,BiGRU预测模型算法表示为:
式中:zt表示重置门;rt表示更新门;Wr,Wz,Ur,Uz分别表示权重矩阵;ht-1表示上一级隐藏层状态;xt表示当前时刻输入向量;ht表示当前时刻候选的隐藏状态;ht表示当前时刻隐藏状态;“*”表示矩阵乘积;σ表示sigmoid函数。
6.根据权利要求5所述的一种改进的变分模态分解与BiGRU融合的水质预测方法,其特征在于,各个模态分量输入到建立BiGRU预测模型中,最后将每个模态分量预测值叠加得到最终预测结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211398569.9A CN116451553A (zh) | 2022-11-09 | 2022-11-09 | 一种改进的变分模态分解与BiGRU融合的水质预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211398569.9A CN116451553A (zh) | 2022-11-09 | 2022-11-09 | 一种改进的变分模态分解与BiGRU融合的水质预测方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN116451553A true CN116451553A (zh) | 2023-07-18 |
Family
ID=87126104
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211398569.9A Pending CN116451553A (zh) | 2022-11-09 | 2022-11-09 | 一种改进的变分模态分解与BiGRU融合的水质预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116451553A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116701822A (zh) * | 2023-08-07 | 2023-09-05 | 云南电网有限责任公司 | 基于麻雀搜索算法与变分模态分解的负荷电能计量方法 |
CN116992256A (zh) * | 2023-09-26 | 2023-11-03 | 中国海洋大学 | 海洋潮汐预测精度提升方法 |
CN117670147A (zh) * | 2024-02-01 | 2024-03-08 | 江西省科学院微生物研究所(江西省流域生态研究所) | 一种湖泊水质预测方法及系统 |
-
2022
- 2022-11-09 CN CN202211398569.9A patent/CN116451553A/zh active Pending
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116701822A (zh) * | 2023-08-07 | 2023-09-05 | 云南电网有限责任公司 | 基于麻雀搜索算法与变分模态分解的负荷电能计量方法 |
CN116701822B (zh) * | 2023-08-07 | 2023-11-28 | 云南电网有限责任公司 | 基于麻雀搜索算法与变分模态分解的负荷电能计量方法 |
CN116992256A (zh) * | 2023-09-26 | 2023-11-03 | 中国海洋大学 | 海洋潮汐预测精度提升方法 |
CN116992256B (zh) * | 2023-09-26 | 2023-12-19 | 中国海洋大学 | 海洋潮汐预测精度提升方法 |
CN117670147A (zh) * | 2024-02-01 | 2024-03-08 | 江西省科学院微生物研究所(江西省流域生态研究所) | 一种湖泊水质预测方法及系统 |
CN117670147B (zh) * | 2024-02-01 | 2024-04-19 | 江西省科学院微生物研究所(江西省流域生态研究所) | 一种湖泊水质预测方法及系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN116451553A (zh) | 一种改进的变分模态分解与BiGRU融合的水质预测方法 | |
Li et al. | Applying various algorithms for species distribution modelling | |
CN109508360B (zh) | 一种基于元胞自动机的地理多元流数据时空自相关分析方法 | |
CN110334875B (zh) | 计及评估指标冲突的风电功率组合概率预测方法 | |
CN109886464B (zh) | 基于优化奇异值分解生成特征集的低信息损失短期风速预测方法 | |
CN110621026B (zh) | 一种基站流量多时刻预测方法 | |
CN112367130B (zh) | 一种无线网络信号传输强度计算方法及计算机存储介质 | |
Goh et al. | A multimodal approach to chaotic renewable energy prediction using meteorological and historical information | |
CN115906954A (zh) | 一种基于图神经网络的多变量时间序列预测方法及装置 | |
CN113592144A (zh) | 一种中长期径流概率预报方法及系统 | |
CN113379107A (zh) | 基于lstm和gcn的区域电离层tec预报方法 | |
CN114692947A (zh) | 一种基于VMD-Lasso-CNN-LSTM的用户用电能耗预测方法 | |
CN115796351B (zh) | 基于变分模态分解和微波衰减的降雨短临预测方法及装置 | |
CN112131731A (zh) | 一种基于空间特征向量滤波的城市生长元胞模拟方法 | |
CN115587666A (zh) | 基于季节趋势分解和混合神经网络的负荷预测方法及系统 | |
CN111476285A (zh) | 一种图像分类模型的训练方法及图像分类方法、存储介质 | |
CN113836823A (zh) | 一种基于负荷分解和优化双向长短期记忆网络的负荷组合预测方法 | |
Bustillo et al. | Predictive Hybridization Model integrating Modified Genetic Algorithm (MGA) and C4. 5 | |
Lu et al. | Quality-relevant feature extraction method based on teacher-student uncertainty autoencoder and its application to soft sensors | |
CN114219096A (zh) | 一种机器学习算法模型的训练方法、装置及存储介质 | |
CN116680548A (zh) | 一种针对多源观测数据的时间序列干旱因果分析方法 | |
CN115936196A (zh) | 基于时序卷积网络的月降水模型预测方法 | |
CN113538239B (zh) | 一种基于时空自回归神经网络模型的插值方法 | |
CN114070262B (zh) | 附加扰动的集成混合集合Kalman滤波天气预报同化方法及装置 | |
CN113761777B (zh) | 一种基于hp-ovmd的超短期光伏功率预测方法 |
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 |