CN110297269B - 一种基于改进型bp神经网络的地震数据双向预测插值方法 - Google Patents
一种基于改进型bp神经网络的地震数据双向预测插值方法 Download PDFInfo
- Publication number
- CN110297269B CN110297269B CN201810242781.3A CN201810242781A CN110297269B CN 110297269 B CN110297269 B CN 110297269B CN 201810242781 A CN201810242781 A CN 201810242781A CN 110297269 B CN110297269 B CN 110297269B
- Authority
- CN
- China
- Prior art keywords
- interpolation
- seismic data
- neural network
- training
- prediction
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 56
- 238000013528 artificial neural network Methods 0.000 title claims abstract description 35
- 230000002457 bidirectional effect Effects 0.000 title claims abstract description 27
- 238000012549 training Methods 0.000 claims description 36
- 238000004422 calculation algorithm Methods 0.000 claims description 11
- 238000012937 correction Methods 0.000 claims description 7
- 230000008569 process Effects 0.000 claims description 6
- 230000003044 adaptive effect Effects 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000004458 analytical method Methods 0.000 claims description 2
- 238000003384 imaging method Methods 0.000 claims description 2
- 230000000694 effects Effects 0.000 abstract description 6
- 230000008859 change Effects 0.000 description 9
- 230000006870 function Effects 0.000 description 7
- 238000011478 gradient descent method Methods 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 210000002569 neuron Anatomy 0.000 description 3
- 230000009471 action Effects 0.000 description 2
- 230000004913 activation Effects 0.000 description 2
- 239000000654 additive Substances 0.000 description 2
- 230000000996 additive effect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 230000001537 neural effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000002945 steepest descent method Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Feedback Control In General (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供了一种基于改进型BP神经网络的地震数据双向预测插值方法,包括:步骤I、获取地震数据信息,根据地震数据信息缺失的实际情况,确定是横向插值或纵向插值,并建立相对应的双向预测模式;步骤II、通过引入附加动量法和自适应学习速率法对BP神经网络改进,建立改进型BP神经网络;步骤III、基于步骤I中建立的双向预测模式和步骤II中建立的改进型BP神经网络,分别从相对的两个方向进行地震数据插值预测,并将预测结果输出;步骤IV、将步骤III中得到的相对的两个方向的地震数据插值预测结果进行融合,得到最终的地震数据插值结果。本发明的方法充分利用多道地震信息,提高了插值精度和效果。
Description
技术领域
本发明属于地震数据处理技术领域,涉及一种基于改进型BP神经网络的地震数据双向预测插值方法。
背景技术
在地震数据采集过程中,地震道缺失和空间采样不足是常见的,因此,在地震数据处理过程中,地震数据插值技术是常用的方法。
目前,常用的地震数据插值技术是基于BP(Back-Propagation)神经网络算法。BP神经网络算法是1985年Rumelhart在误差反向传播理论的基础上提出反向传播(Back-Propagation)学习算法。BP神经网络构建以多层前馈网络为基础,由输入、输出、隐含层组成。输入信号在各层神经元间传递,通过各神经元激活函数的作用,完成正向传播,同时根据正向传播的误差信号,BP神经网络利用梯度下降法进行有导师式反向传播学习,通过反复的正向、反向学习,使误差信号达到最小。BP神经网络具有强大的数据识别和模拟能力,在解决非线性系统问题时,优势明显,是一种可以被推广应用的前沿理论与技术。
然而,由于BP神经网络采用非线性规划中的最速下降方法,按误差函数的负梯度方向修改权值,因而通常存在以下问题:学习效率低,收敛速度慢,易陷入局部极小状态。因而,如果采用BP神经算法直接进行地震插值预测,其插值结果必然与实际存在一定的偏差。
另外,目前常用的插值方法主要是多项式插值和样条插值,多项式插值简单可逼近,然而,其仅考虑的整体性而忽略了局部独立性;样条插值是一种分段多项式的插值算法,各相邻段上的多项式又具有某种连接性,既保持了多项式的简单和逼近可行性,又在各段之间保持了局部的独立性,但是其无法充分利用多道地震信息,限制了其在地震资料中处理中的应用。
发明内容
为了解决上述问题,本发明基于附加动量法和自适应学习速率法的改进型BP神经网络算法,并且提出使用双向预测插值的方法,兼顾了整体性和局部独立性,在地震数据处理过程中,充分利用了多道地震信息,提高了插值精度和效果。
本发明的通过如下技术方案实现:
一种基于改进型BP神经网络的地震数据双向预测插值方法,包括:
步骤I、获取地震信息,根据实际插值的需要,确定插值方式,建立相应的双向预测模式;
步骤II、通过引入附加动量法和自适应学习速率法对BP神经网络改进,建立改进型BP神经网络;
步骤III、基于步骤I中建立的双向预测模式和步骤II中建立的改进型BP神经网络,分别从相对的两个方向进行地震数据插值预测,并将预测结果输出;
步骤IV、将步骤III中得到的相对的两个方向的地震数据插值预测结果进行融合,得到最终的地震数据插值结果。
根据本发明,所述双向预测模式包括左右双向预测模式或上下双向预测模式。
根据本发明,所述地震信息包括地震速度场、地震叠前炮记录和地震成像数据体。根据地震信息中缺失的数据,确定插值方式,插值方式包括横向插值和纵向插值。当横向插值时对应的为左右双向预测模式,当纵向插值时对应的为上下双向预测模式。
根据本发明,在步骤II建立改进型BP神经网络之后还包括确定预测参数,并输入已获取的地震数据信息进行学习训练;所述预测参数包括输入层、隐含层、输出层、预测误差和迭代次数。一般地,预测误差和迭代次数人为给定,一般预测误差控制在10-6以内,迭代次数在1000以内。
在进行学习训练的过程中,每训练一次,样本的特征矢量就要修正一次权值。如果在最大训练次数以内,网络的误差降到期望误差之下,训练停止,而且认为网络己经收敛;如果训练次数超过最大训练次数,而网络误差还没有降到期望误差之下,则训练停止,且认为网络不收敛。通常情况下只要训练次数足够大,一般网络都会收敛,但有时收敛时间很长。因此应当在达到误差要求的前提下选取适当的训练次数。
根据本发明的优选实施方式,在步骤IV之后还包括步骤V:对地震数据插值结果进行误差分析。
在步骤III中,分别从相对的两个方向进行地震数据插值预测是指,当横向插值时,包括从左向右和从右向左两个方向进行地震数据插值预测;当纵向插值时,包括从上到下和从下到上两个方向进行地震数据插值预测。
根据本发明的优选实施方式,在步骤IV中,采用反距离加权的方式将两个方向的地震数据插值预测结果进行融合。
基本BP算法包括两个方面:信号的正向传播和误差的反向传播,即计算实际输出时按从输入到输出的方向进行,而权值和阈值的修正从输出到输入的方向进行。图6是BP神经网络的网络结构图,其中:
xj表示输入层第j个节点的输入,j=1,…,M;
wij表示隐含层第i个节点到输入层第j个节点之间的权值,i=1,…,q;j=1,…,M;
θi表示隐含层第i个节点的阈值;
φ(x)表示隐含层的激活函数;
wki表示输出层第k个节点到隐含层第i个节点之间的权值,k=1,…,L;i=1,…,q;
ak表示输出层第k个节点的阈值,k=1,…,L;
ψ(x)表示输出层的激励函数;
ok表示输出层第k个节点的输出,k=1,…,L。
以上字母所表示的含义在以下式(4-1)-(4-20)中具有相同含义。
(1)信号的前向传播过程
隐含层第i个节点的输入neti:
隐含层第i个节点的输出yi:
输出层第k个节点的输入netk:
输出层第k个节点的输出ok:
(2)误差的反向传播过程
误差的反向传播,即首先由输出层开始逐层计算各层神经元的输出误差,然后根据误差梯度下降法来调节各层的权值和阈值,使修改后的网络的最终输出能接近期望值。
对于每一个样本p的二次型误差准则函数为Ep:
其中,Tk是第k个节点的期望误差。
系统对P个训练样本的总误差准则函数E为:
根据误差梯度下降法依次修正输出层权值的修正量Δwki,输出层阈值的修正量Δak,隐含层权值的修正量Δwij,隐含层阈值的修正量Δθi。
输出层权值调整公式:
输出层阈值调整公式:
隐含层权值调整公式:
隐含层阈值调整公式:
又因为:
其中,φ'(x)是φ(x)的一阶导数,ψ'(x)是ψ(x)的一阶导数。
所以最后得到以下公式:
根据本发明,所述改进型BP神经网络是对现有的BP神经网络中的误差的反向传播过程中的权值和阈值的修正步骤进行改进,附加动量法使网络在修正其权值时,不仅考虑误差在梯度上的作用,而且考虑在误差曲面上变化趋势的影响。在没有附加动量的作用下,网络可能陷入浅的局部极小值,利用附加动量的作用有可能滑过这些极小值。该方法是在反向传播法的基础上在每一个权值(或阈值)的变化上加上一项正比于前次权值(或阈值)变化量的值,并根据反向传播法来产生新的权值(或阈值)变化。引入动量因子和学习速率后的权值和阈值的调节公式为:
ΔW(d+1)=(1-mc)ηδjpk+mcΔW(d) (4-21)
Δb(d+1)=(1-mc)ηδj+mcΔb(d) (4-22)
其中,d为训练次数,mc为动量因子,δj为输出节点j的计算误差,pk为输入节点k的计算输入,ΔW(d)为第d次训练时权值增量,Δb(d)为第d次训练时阈值增量,η为学习速率。
附加动量法的实质是将最后一次权值(或阈值)变化的影响,通过一个动量因子来传递。当动量因子取值为零时,权值(或阈值)的变化仅是根据梯度下降法产生;当动量因子取值为1时,新的权值(或阈值)变化则是设置为最后一次权值(或阈值)的变化,而依梯度法产生的变化部分则被忽略掉了。以此方式,当增加了动量项后,促使权值的调节向着误差曲面底部的平均方向变化,当网络权值进入误差曲面底部的平坦区时,δi将变得很小,于是ΔW(d+1)=ΔW(d),从而防止了ΔWjk(d)=0的出现,有助于使网络从误差曲面的局部极小值中跳出。
根据附加动量法的设计原则,当修正的权值在误差中导致太大的增长结果时,新的权值应被取消而不被采用,并使动量作用停止下来,以使网络不进入较大误差曲面;当新的误差变化率对其旧值超过一个事先设定的最大误差变化率时,也得取消所计算的权值变化。其最大误差变化率可以是任何大于或等于1的值。典型的取值取1.04。所以,在进行附加动量法的训练程序设计时,必须加进条件判断以正确使用其权值修正公式。
根据本发明的优选实施方式,采用附加动量法的判断条件是:
其中,mc为动量因子,E(d)为第d次训练时误差平方和。
对于一个特定的问题,要选择适当的学习速率不是一件容易的事情。通常是凭经验或实验获取,但即使这样,对训练开始初期功效较好的学习速率,不见得对后来的训练合适。为了解决这个问题,人们自然想到在训练过程中,自动调节学习速率。通常调节学习速率的准则是:检查权值是否真正降低了误差函数,如果确实如此,则说明所选学习速率小了,可以适当增加一个量;若不是这样,而产生了过调,则应该减少学习速率的值。根据本发明的优选实施方式,自适应学习速率的调整公式为:
其中,η(d)为第d次训练时的学习速率,E(d)为第d次训练时误差平方和。
本发明的有益效果:
本发明地震数据预测插值方式使基于改进型BP神经网络,通过引入自适应学习速率法和附加动量法对传统BP神经网络进行改进,解决了传统BP神经网络算法中存在的学习效率低和易陷入局部极值的问题。
另外,本发明在地震数据预测插值方法中,提出了双向预测插值技术,采用反距离加权的方法将相对两个方向的地震预测结果进行融合,得到最终的地震数据插值结果,在地震数据预测插值技术中使用该方法兼顾了整体性和局部性,且充分利用多道地震信息,提高了插值精度和效果。使用本发明提供的方法进行地震数据插值,预测插值与实际结果的相对误差均在1.5%以内,精确度高。
附图说明
图1是本发明实施例1中真实速度场与预测插值前后速度场的比较;
其中,A是真实速度场,B是预测插值前速度场,C是预测插值后速度场。
图2是205道真实速度和预测速度对比。
图3是205道预测速度相对误差。
图4是215道真实速度和预测速度对比。
图5是215道预测速度相对误差。
图6是传统的BP神经网络的网络结构图。
具体实施方式
以下结合实施例对本发明进行详细说明,但本发明并不受下述实施例限定。
实施例1
一种基于改进型BP神经网络的地震数据双向预测插值方法,首先获地震速度场信息,根据图2中的地震数据信息可以确定插值方式为横向插值,进而可以确定双向预测模式为左右双向预测。基于引入附加动量法和自适应学习速率法后获得的改进型BP神经网络进行地震数据预测插值,选择8个训练样本数,用这8个点参与训练,每次训练用2点预测1点,训练结束形成训练模式,按照形成的模式进行缺失地震数据的预测;
首先从左向右预测地震数据插值,将所预测的地震数据信息输出;然后从右向左预测地震数据插值,将所预测的地震信息输出。将上述两个相对的方向预测得到的插值结果按照反距离加权方式进行融合,得到最终的预测插值结果。
选取205道和215道的真实速度与预测速度进行对比,真实速度和预测速度的比较如图2和4所示,预测的相对误差如图3和图5所示,使用本发明提供的方法预测插值与实际结果的相对误差都在1.5%以内,说明使用本发明的方法预测插值的速度场精确高。
以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (9)
1.一种基于改进型BP神经网络的地震数据双向预测插值方法,包括:
步骤I、获取地震数据信息,根据地震数据信息缺失的实际情况,确定插值方式,并建立相对应的双向预测模式,所述双向预测模式包括左右双向预测模式或上下双向预测模式,所述插值方式包括横向插值和竖向插值,当横向插值时,相对应的预测模式为左右双向预测,当纵向插值时,相对应的预测模式为上下双向预测;
步骤II、通过引入附加动量法和自适应学习速率法对BP神经网络改进,建立改进型BP神经网络;
步骤III、基于步骤I中建立的双向预测模式和步骤II中建立的改进型BP神经网络,分别从相对的两个方向进行地震数据插值预测,并将预测结果输出;
步骤IV、将步骤III中得到的相对的两个方向的地震数据插值预测结果进行融合,得到最终的地震数据插值结果。
2.根据权利要求1所述的方法,其特征在于,所述步骤II中,建立改进型BP神经网络后,还包括确定预测参数,并输入已获取的地震数据信息进行学习训练。
3.根据权利要求2所述的方法,其特征在于,所述步骤II中,所述预测参数包括输入层、隐含层、输出层、预测误差和迭代次数。
4.根据权利要求1-3中任一项所述的方法,其特征在于,所述地震信息包括地震速度场、地震叠前炮记录和地震成像数据体。
5.根据权利要求1-3中任一项所述的方法,其特征在于,在步骤IV之后还包括步骤V:对地震数据插值结果进行误差分析。
6.根据权利要求1-3中任一项所述的方法,其特征在于,在步骤IV中,采用反距离加权的方式将两个方向的地震数据插值预测结果进行融合。
7.根据权利要求1-3中任一项所述的方法,其特征在于,所述改进型BP神经网络是对现有的BP神经网络中的误差的反向传播过程中的权值和阈值的修正项进行改进,包括引入动量因子和学习速率,改进型BP算法中的权值和阈值的调节公式为:
ΔW(d+1)=(1-mc)ηδjpk+mcΔW(d) (4-21)
Δb(d+1)=(1-mc)ηδj+mcΔb(d) (4-22)
其中,d为训练次数,mc为动量因子,δj为输出节点j的计算误差,pk为输入节点k的计算输入,ΔW(d)为第d次训练时权值增量,Δb(d)为第d次训练时阈值增量,η为学习速率。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810242781.3A CN110297269B (zh) | 2018-03-23 | 2018-03-23 | 一种基于改进型bp神经网络的地震数据双向预测插值方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810242781.3A CN110297269B (zh) | 2018-03-23 | 2018-03-23 | 一种基于改进型bp神经网络的地震数据双向预测插值方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110297269A CN110297269A (zh) | 2019-10-01 |
CN110297269B true CN110297269B (zh) | 2020-12-04 |
Family
ID=68025849
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810242781.3A Active CN110297269B (zh) | 2018-03-23 | 2018-03-23 | 一种基于改进型bp神经网络的地震数据双向预测插值方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110297269B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110766234B (zh) * | 2019-10-29 | 2023-06-20 | 济南大学 | 基于信息融合的水泥冷却过程篦下压力预测方法 |
CN113624998A (zh) * | 2021-09-17 | 2021-11-09 | 辽宁旭能科技有限公司 | 基于电力大数据的电锅炉补热储热成本优化方法和装置 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2000002769A (ja) * | 1998-06-17 | 2000-01-07 | Nippon Koei Co Ltd | 地質構造の空間分布を予測して地質図を作成する方法及び装置 |
WO2010151354A1 (en) * | 2009-06-26 | 2010-12-29 | Exxonmobil Upstream Research Company | Constructing resistivity models from stochastic inversion |
CN105353412A (zh) * | 2015-12-14 | 2016-02-24 | 中国石油大学(华东) | 一种井震联合平均速度场的计算方法及系统 |
CN107167843A (zh) * | 2017-05-25 | 2017-09-15 | 中国石油天然气集团公司 | 多波时间域匹配方法及装置 |
CN107607992A (zh) * | 2017-08-24 | 2018-01-19 | 电子科技大学 | 基于卷积神经网络的多波匹配方法 |
CN107664773A (zh) * | 2017-09-26 | 2018-02-06 | 武汉大学 | 一种基于时间移位和熵约束的地震信号编码方法 |
-
2018
- 2018-03-23 CN CN201810242781.3A patent/CN110297269B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2000002769A (ja) * | 1998-06-17 | 2000-01-07 | Nippon Koei Co Ltd | 地質構造の空間分布を予測して地質図を作成する方法及び装置 |
WO2010151354A1 (en) * | 2009-06-26 | 2010-12-29 | Exxonmobil Upstream Research Company | Constructing resistivity models from stochastic inversion |
CN105353412A (zh) * | 2015-12-14 | 2016-02-24 | 中国石油大学(华东) | 一种井震联合平均速度场的计算方法及系统 |
CN107167843A (zh) * | 2017-05-25 | 2017-09-15 | 中国石油天然气集团公司 | 多波时间域匹配方法及装置 |
CN107607992A (zh) * | 2017-08-24 | 2018-01-19 | 电子科技大学 | 基于卷积神经网络的多波匹配方法 |
CN107664773A (zh) * | 2017-09-26 | 2018-02-06 | 武汉大学 | 一种基于时间移位和熵约束的地震信号编码方法 |
Non-Patent Citations (3)
Title |
---|
"3D interpolation of irregular data with POCS algorithm";Ray A.等;《Geophysics》;20061231;第71卷(第6期);第E91-E97页 * |
"基于粒子群优化的BP网络在地震属性融合技术中的应用";曹琳昱 等;《石油与天然气地质》;20101031;第31卷(第5期);第685-688页 * |
"自适应与附加动量BP神经网络的ECT流型辨识";王莉莉 等;《哈尔滨理工大学学报》;20180228;第23卷(第1期);第105-110页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110297269A (zh) | 2019-10-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108051999B (zh) | 基于深度强化学习的加速器束流轨道控制方法及系统 | |
CN103105246A (zh) | 一种基于遗传算法改进的bp神经网络的温室环境预测反馈方法 | |
CN106815782A (zh) | 一种基于神经网络统计学模型的房地产估值方法及系统 | |
CN110297269B (zh) | 一种基于改进型bp神经网络的地震数据双向预测插值方法 | |
CN109886343B (zh) | 图像分类方法及装置、设备、存储介质 | |
CN110427690A (zh) | 一种基于全局粒子群算法生成ato速度曲线的方法及装置 | |
CN110276441A (zh) | 一种基于深度学习的梯形重叠核脉冲估计方法 | |
CN104122796A (zh) | 一种智能装配序列规划方法 | |
CN112270442A (zh) | 一种基于ivmd-acmpso-cslstm组合电力负荷预测方法 | |
CN112163671A (zh) | 一种新能源场景生成方法及系统 | |
CN112348155A (zh) | 一种模糊神经网络模型的优化方法及系统 | |
CN111290282A (zh) | 火电机组协调系统的预见式预测控制方法 | |
CN112257847A (zh) | 一种基于CNN和LSTM预测地磁Kp指数的方法 | |
CN113762464B (zh) | 一种基于学习的列车运行参考曲线动态生成方法 | |
CN113420508B (zh) | 一种基于lstm的机组组合计算方法 | |
CN110782181A (zh) | 一种低压台区线损率的计算方法及可读存储介质 | |
CN109886405A (zh) | 一种抑制噪声的基于人工神经网络结构优化方法 | |
CN107273971B (zh) | 基于神经元显著性的前馈神经网络结构自组织方法 | |
CN105550774A (zh) | 设施投料预测方法和系统 | |
CN111091234A (zh) | 一种基于遗传和神经网络算法的电力市场电价预测方法 | |
CN110779526B (zh) | 一种路径规划方法、装置及存储介质 | |
CN112270270A (zh) | 一种基于abc优化神经网络的接触网故障类型识别方法 | |
CN111562740A (zh) | 基于利用梯度的多目标强化学习算法的自动控制方法 | |
CN108446506B (zh) | 一种基于区间反馈神经网络的不确定系统建模方法 | |
CN110991604A (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |