CN111352158B - 地震信号增强方法及装置 - Google Patents
地震信号增强方法及装置 Download PDFInfo
- Publication number
- CN111352158B CN111352158B CN201811580504.XA CN201811580504A CN111352158B CN 111352158 B CN111352158 B CN 111352158B CN 201811580504 A CN201811580504 A CN 201811580504A CN 111352158 B CN111352158 B CN 111352158B
- Authority
- CN
- China
- Prior art keywords
- seismic data
- ovt
- domain
- ovt domain
- gather
- 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 33
- 238000012545 processing Methods 0.000 claims abstract description 40
- 230000003044 adaptive effect Effects 0.000 claims abstract description 13
- 238000004590 computer program Methods 0.000 claims description 16
- 238000005070 sampling Methods 0.000 claims description 7
- 230000001427 coherent effect Effects 0.000 claims description 4
- 238000005516 engineering process Methods 0.000 abstract description 3
- 230000002708 enhancing effect Effects 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 19
- 230000005012 migration Effects 0.000 description 6
- 238000013508 migration Methods 0.000 description 6
- 230000006870 function Effects 0.000 description 5
- 238000013459 approach Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000008030 elimination Effects 0.000 description 2
- 238000003379 elimination reaction Methods 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000001629 suppression Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 230000009466 transformation Effects 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
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/362—Effecting static or dynamic corrections; Stacking
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)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供一种地震信号增强方法及装置,包括:获取CMP道集地震数据;将CMP道集地震数据分选至OVT域,得到OVT域道集地震数据;对OVT域道集地震数据中存在的空道地震数据进行插值运算,获得插值后的OVT域道集地震数据;对插值后的OVT域道集地震数据中的每一个插值后的OVT域道集地震数据分别进行开方处理、自适应三维倾斜叠加处理和幂次运算处理,获得去噪后的OVT域道集地震数据;将去噪后的OVT域道集地震数据分选至CMP域,得到去噪后的CMP域道集地震数据。该方案对地震数据开方,实现低信噪比地震信号的增强。将自适应三维倾斜叠加去噪技术应用于OVT域道集,进一步提高地震数据的信噪比。
Description
技术领域
本发明涉及石油地震勘探技术领域,特别涉及一种地震信号增强方法及装置。
背景技术
随着勘探目标构造更加复杂,剧烈起伏地表和逆冲断裂复杂构造等是目前地震勘探的重点区域。这类地区地震激发和接收条件非常差,采集到的地震反射信号几乎被强噪音淹没,从而对地震数据信噪比有更高的要求。然而低信噪比地震资料,即使采用多域多维的叠前噪音压制技术,其地震资料信噪比仍难以满足地质的要求。因此,提高低信噪比地震资料信噪比对油气勘探非常重要,实现弱信号增强成为地震信号处理的关键。
针对弱信号增强,去除强噪声干扰,国内外的学者进行了研究。地震勘探噪声消除是根据采集到地震信号的特点以及噪声特征,采用不同的方法实现噪音的消除。比如一维频率域去噪方法和f-k域滤波方法就是利用噪声与有效信号在频率、视速度方面的差异,实现噪音的消除,而f-k域滤波法去噪的同时伴随假频现象的产生,且去噪后的剖面保真度也受到破坏。因此,在去噪过程中尽量使得假频降到最低。另外,目前,应用小波变换去噪得到比较广泛的应用,即小波变换是基于多尺度、多分辨率的时频分解去噪方法能够有效进行去噪处理、小波变换与多重自相关相结合的强白噪声(或有色噪声)背景下微弱信号的提取。但是,小波变换对一维数据具有比较好的局部时频分析能力,在二维(或更高维)的信号处理方面,则存在一定的局限性;而且常用的二维小波也仅仅是一维小波的张量积,它对二维信号中的点状奇异性具有很好的刻画能力,而对于直线(或曲线)等存在的边缘特征,则难以用小波变换精确刻画,对于弱地震信号效果欠佳。
发明内容
本发明实施例提供了一种地震信号增强方法及装置,可以提高地震数据信噪比。
本发明实施例提供了一种地震信号增强方法,该方法包括:
获取CMP道集地震数据;
将所述CMP道集地震数据分选至OVT域,得到OVT域道集地震数据;
对所述OVT域道集地震数据中存在的空道地震数据进行插值运算,获得插值后的OVT域道集地震数据;
对所述插值后的OVT域道集地震数据中的每一个插值后的OVT域道集地震数据分别进行开方处理、自适应三维倾斜叠加处理和幂次运算处理,获得去噪后的OVT域道集地震数据;
将所述去噪后的OVT域道集地震数据分选至CMP域,得到去噪后的CMP域道集地震数据。
本发明实施例还提供了一种地震信号增强装置,该装置包括:
地震数据获取模块,用于获取CMP道集地震数据;
分选模块,用于将所述CMP道集地震数据分选至OVT域,得到OVT域道集地震数据;
插值运算模块,用于对所述OVT域道集地震数据中存在的空道地震数据进行插值运算,获得插值后的OVT域道集地震数据;
去噪模块,用于对所述插值后的OVT域道集地震数据中的每一个插值后的OVT域道集地震数据分别进行开方处理、自适应三维倾斜叠加处理和幂次运算处理,获得去噪后的OVT域道集地震数据;
所述分选模块还用于:
将所述去噪后的OVT域道集地震数据分选至CMP域,得到去噪后的CMP域道集地震数据。
本发明实施例还提供了一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现上述所述地震信号增强方法。
本发明实施例还提供了一种计算机可读存储介质,所述计算机可读存储介质存储有执行上述所述地震信号增强方法的计算机程序。
在一个实施例中,利用任意正实数开无穷次根号其结果趋近于零的特点,对地震数据进行开方运算,实现低信噪比地震信号的增强。根据插值后的OVT域数据与叠后数据的相似性,把倾斜叠加去噪应用于OVT域的三维叠前数据,进一步改善叠前数据的信噪比,从而提高偏移精度和成像质量。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本发明实施例提供的一种地震信号增强方法流程图;
图2是本发明实施例提供的一种去噪前作为输入的CMP道集示意图;
图3是本发明实施例提供的一种分选至OVT域的道集示意图;
图4是本发明实施例提供的一种对图2的OVT域道集进行空道插值示意图;
图5是本发明实施例提供的一种去噪前的OVT片示意图;
图6是本发明实施例提供的一种去噪后的OVT片示意图;
图7是本发明实施例提供的一种去噪后输出的CMP道集示意图;
图8是本发明实施例提供的一种去噪前的S/N示意图;
图9是本发明实施例提供的一种去噪后的S/N示意图;
图10是本发明实施例提供的一种去噪前的叠加剖面示意图;
图11是本发明实施例提供的一种去噪后的叠加剖面示意图;
图12是本发明实施例提供的一种去噪前的叠前时间偏移剖面示意图;
图13是本发明实施例提供的一种去噪后的叠前时间偏移剖面示意图;
图14是本发明实施例提供的一种地震信号增强装置结构示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
在本发明实施例中,为实现本发明的目的,提供了一种地震信号增强方法,该方法基于能量均衡原理和自适应三维倾斜叠加的地震信号增强技术,利用任意正实数开无穷次根号其结果趋近于零的特点,对地震数据进行开方运算,实现低信噪比地震信号的增强。然后根据插值后的OVT域数据与叠后数据的相似性,把倾斜叠加去噪应用于OVT域的三维叠前数据,进一步改善叠前数据的信噪比,从而提高偏移精度和成像质量,其实现过程包括以下步骤,如图1所示:
步骤101:获取CMP道集(Commom middle point,共中心点道集)地震数据;
步骤102:将所述CMP道集地震数据分选至OVT(OFFset Vector Tiles,炮检距矢量片)域,得到OVT域道集地震数据;
步骤103:对所述OVT域道集地震数据中存在的空道地震数据进行插值运算,获得插值后的OVT域道集地震数据;
步骤104:对所述插值后的OVT域道集地震数据中的每一个插值后的OVT域道集地震数据分别进行开方处理、自适应三维倾斜叠加处理和幂次运算处理,获得去噪后的OVT域道集地震数据;
步骤105:将所述去噪后的OVT域道集地震数据分选至CMP域,得到去噪后的CMP域道集地震数据。
在本发明实施例中,步骤104具体执行如下:
(1)进行开方运算:
保持地震数据的符号不变,对振幅进行开方处理,开方的计算公式为:
其中,x(t)表示插值后的OVT域道集地震数据;N表示开方的次数。
2)进行自适应三维倾斜叠加处理:
利用地震信号的相干性和随机噪声的不相干性,实现噪声的压制。首先进行倾角扫描,设置时窗、时窗滑动增量、Inline和Xline方向上的倾角范围等参数,进行相干性分析,计算相干值,相干值计算公式为:
其中,NL为L处的相干值,Aij为第i道采样点j的地震数据x(t)的振幅,t为道长,n为OVT域道集地震数据总道数,Kj为采样点j的道数。
然后,通过相干值确定地下真实反射界面的倾角等信息,再根据倾角确定的面与反射界面重合时相关性最大的特点,采用最大相关准则自适应地确定反射界面,对该反射界面内数据进行加权叠加获得纯信号道SP(表示某道的地震数据),纯信号道(地震数据)与输入数据加权,获得输出数据:
其中,SI为输入道,SP为纯信号道,W为相干加权系数,So为输出道。
3)进行幂次运算处理:
保持地震数据的符号不变,对去噪后的OVT域道集地震数据的振幅进行幂次运算,得到去噪后的CMP域道集地震数据。幂次运算计算公式为:
其中,x(t)表示插值后的OVT域道集地震数据,此时的x(t)表示的是So中的地震数据,也就是指的去噪后的CMP域道集地震数据;N的数值与公式(1)中相同。
实施例
按照如图1所示的方法进行地震信号增强,对输入的CMP道集(图2)分选至OVT域,得到图3的OVT域道集;由于图3的OVT域道集存在较多的空道,因此需要对其插值,得到图4所示的插值后的OVT域道集。对每一个OVT片(图5)分别进行开方、自适应三维倾斜叠加去噪、幂次运算等处理,获得去噪后的OVT片(图6)。组合去噪后的OVT片并重新分选到CMP域,得到去噪后的CMP道集(图7)。图8是去噪前的信噪比,图9是去噪后的信噪比,可以看出,去噪后信噪比得到较大的提高。图10是去噪前的叠加剖面,图11是去噪后的叠加剖面,可以看出,去噪后同相轴的连续性更好,特别深层的弱信号得到增强。图12是去噪前的叠前时间偏移剖面,图13是去噪后的叠前时间偏移剖面,可以看出信噪比得到大幅提高,断裂构造更加清晰,深层的弱信号得到增强,说明本发明能够有效实现地震信号的增强。其中,图2至图13中的横坐标表示inline线号,纵坐标表示深度。
基于同一发明构思,本发明实施例中还提供了一种地震信号增强装置,如下面的实施例所述。由于地震信号增强装置解决问题的原理与地震信号增强方法相似,因此地震信号增强装置的实施可以参见地震信号增强方法的实施,重复之处不再赘述。以下所使用的,术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现,但是硬件,或者软件和硬件的组合的实现也是可能并被构想的。
图14是本发明实施例的地震信号增强装置的一种结构框图,如图14所示,包括:
地震数据获取模块1401,用于获取CMP道集地震数据;
分选模块1402,用于将所述CMP道集地震数据分选至OVT域,得到OVT域道集地震数据;
插值运算模块1403,用于对所述OVT域道集地震数据中存在的空道地震数据进行插值运算,获得插值后的OVT域道集地震数据;
去噪模块1404,用于对所述插值后的OVT域道集地震数据中的每一个插值后的OVT域道集地震数据分别进行开方处理、自适应三维倾斜叠加处理和幂次运算处理,获得去噪后的OVT域道集地震数据;
所述分选模块1402还用于:
将所述去噪后的OVT域道集地震数据分选至CMP域,得到去噪后的CMP域道集地震数据。
在本发明实施例中,所述去噪模块1404具体用于:按照如下公式对所述插值后的OVT域道集地震数据中的每一个插值后的OVT域道集地震数据进行开方处理:
其中,x(t)表示插值后的OVT域道集地震数据;N表示开方的次数。
在本发明实施例中,所述去噪模块1404具体用于:按照如下公式对所述插值后的OVT域道集地震数据中的每一个插值后的OVT域道集地震数据进行自适应三维倾斜叠加处理:
其中,NL为L处的相干值,Aij为第i道采样点j的振幅,t为道长,n为OVT域道集地震数据总道数,Kj为采样点j的道数;
其中,SI为输入道,SP为纯信号道,W为相干加权系数,So为输出道。
在本发明实施例中,所述去噪模块1404具体用于:按照如下公式对所述插值后的OVT域道集地震数据中的每一个插值后的OVT域道集地震数据进行幂次运算处理:
其中,x(t)表示插值后的OVT域道集地震数据;N表示开方的次数。
本发明实施例还提供了一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现上述所述地震信号增强方法。
本发明实施例还提供了一种计算机可读存储介质,所述计算机可读存储介质存储有执行上述所述地震信号增强方法的计算机程序。
综上所述,本发明提出的地震信号增强方法及装置具有如下优点:
利用任意正实数开无穷次根号其结果趋近于零的特点,对地震数据进行开方运算,实现低信噪比地震信号的增强。然后根据插值后的OVT域数据与叠后数据的相似性,把自适应三维倾斜叠加去噪应用于OVT域的三维叠前数据,进一步改善叠前数据的信噪比。经过两次信号增强,能够大幅提高地震数据的信噪比,使弱地震信号得到增强。
本领域内的技术人员应明白,本发明的实施例可提供为方法、系统、或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明实施例可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (6)
1.一种地震信号增强方法,其特征在于,包括:
获取CMP道集地震数据;
将所述CMP道集地震数据分选至OVT域,得到OVT域道集地震数据;
对所述OVT域道集地震数据中存在的空道地震数据进行插值运算,获得插值后的OVT域道集地震数据;
对所述插值后的OVT域道集地震数据中的每一个插值后的OVT域道集地震数据分别进行开方处理、自适应三维倾斜叠加处理和幂次运算处理,获得去噪后的OVT域道集地震数据;
将所述去噪后的OVT域道集地震数据分选至CMP域,得到去噪后的CMP域道集地震数据;
按照如下公式对所述插值后的OVT域道集地震数据中的每一个插值后的OVT域道集地震数据进行开方处理:
按照如下公式对所述插值后的OVT域道集地震数据中的每一个插值后的OVT域道集地震数据进行幂次运算处理:
其中,x(t)表示插值后的OVT域道集地震数据;N表示开方的次数。
3.一种地震信号增强装置,其特征在于,包括:
地震数据获取模块,用于获取CMP道集地震数据;
分选模块,用于将所述CMP道集地震数据分选至OVT域,得到OVT域道集地震数据;
插值运算模块,用于对所述OVT域道集地震数据中存在的空道地震数据进行插值运算,获得插值后的OVT域道集地震数据;
去噪模块,用于对所述插值后的OVT域道集地震数据中的每一个插值后的OVT域道集地震数据分别进行开方处理、自适应三维倾斜叠加处理和幂次运算处理,获得去噪后的OVT域道集地震数据;
所述分选模块还用于:
将所述去噪后的OVT域道集地震数据分选至CMP域,得到去噪后的CMP域道集地震数据;
所述去噪模块具体用于:按照如下公式对所述插值后的OVT域道集地震数据中的每一个插值后的OVT域道集地震数据进行开方处理:
所述去噪模块具体用于:按照如下公式对所述插值后的OVT域道集地震数据中的每一个插值后的OVT域道集地震数据进行幂次运算处理:
其中,x(t)表示插值后的OVT域道集地震数据;N表示开方的次数。
5.一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现权利要求1至2任一所述地震信号增强方法。
6.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质存储有执行权利要求1至2任一所述地震信号增强方法的计算机程序。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811580504.XA CN111352158B (zh) | 2018-12-24 | 2018-12-24 | 地震信号增强方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811580504.XA CN111352158B (zh) | 2018-12-24 | 2018-12-24 | 地震信号增强方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111352158A CN111352158A (zh) | 2020-06-30 |
CN111352158B true CN111352158B (zh) | 2023-02-10 |
Family
ID=71195215
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811580504.XA Active CN111352158B (zh) | 2018-12-24 | 2018-12-24 | 地震信号增强方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111352158B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114076985B (zh) * | 2020-08-18 | 2024-04-09 | 中国石油化工股份有限公司 | 可控震源的噪音压制方法、装置、设备和存储介质 |
-
2018
- 2018-12-24 CN CN201811580504.XA patent/CN111352158B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN111352158A (zh) | 2020-06-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Huang et al. | Signal extraction using randomized-order multichannel singular spectrum analysis | |
Huang et al. | Mathematical morphological filtering for linear noise attenuation of seismic data | |
Liu et al. | A 1D time-varying median filter for seismic random, spike-like noise elimination | |
US10436924B2 (en) | Denoising seismic data | |
Klokov et al. | Selecting an optimal aperture in Kirchhoff migration using dip-angle images | |
CN109143363B (zh) | 海洋拖缆双检采集鬼波压制方法及系统 | |
Walda et al. | Determination of wavefront attributes by differential evolution in the presence of conflicting dips | |
Wang et al. | Separation and imaging of seismic diffractions using a localized rank-reduction method with adaptively selected ranks | |
Al-Dossary et al. | Lineament-preserving filtering | |
Huang et al. | Random noise attenuation by planar mathematical morphological filtering | |
CN111736224B (zh) | 一种压制叠前地震资料线性干扰方法、存储介质及设备 | |
Wang et al. | Robust singular value decomposition filtering for low signal-to-noise ratio seismic data | |
CN111781635B (zh) | 海底四分量弹性波高斯束深度偏移方法和装置 | |
CN111352158B (zh) | 地震信号增强方法及装置 | |
Mousa et al. | Enhancement of first arrivals using the τ-p transform on energy-ratio seismic shot records | |
Pilkington et al. | Noise reduction procedures for gravity-gradiometer data | |
CN111257931B (zh) | 一种去除海洋地震勘探过船干扰噪音的方法 | |
Eichkitz et al. | Enhanced coherence attribute imaging by structurally oriented filtering | |
Wu et al. | Applications of the trace transform in surface wave attenuation on seismic records | |
CN109782344B (zh) | 沉积层序边界识别方法及装置 | |
CN112684498A (zh) | 一种基于宽方位地震数据的储层裂缝预测方法及系统 | |
Tian et al. | Joint sparsity recovery for noise attenuation | |
CN112083490B (zh) | 地震数据噪声衰减方法及装置 | |
Deng et al. | Automatic multihorizons recognition for seismic data based on Kalman filter tracker | |
CN110596756B (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 |