CN111366978A - 一种基于多次挤压小波变换的地震时频分析方法及系统 - Google Patents
一种基于多次挤压小波变换的地震时频分析方法及系统 Download PDFInfo
- Publication number
- CN111366978A CN111366978A CN202010354535.4A CN202010354535A CN111366978A CN 111366978 A CN111366978 A CN 111366978A CN 202010354535 A CN202010354535 A CN 202010354535A CN 111366978 A CN111366978 A CN 111366978A
- Authority
- CN
- China
- Prior art keywords
- frequency
- seismic
- time
- dimensional
- wavelet transform
- 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
Links
- 238000001125 extrusion Methods 0.000 title claims abstract description 42
- 238000004458 analytical method Methods 0.000 title claims abstract description 36
- 238000001228 spectrum Methods 0.000 claims abstract description 30
- 238000000034 method Methods 0.000 claims abstract description 12
- 230000009466 transformation Effects 0.000 claims abstract description 12
- 238000004590 computer program Methods 0.000 claims description 15
- 230000008569 process Effects 0.000 claims description 5
- 238000003860 storage Methods 0.000 claims description 4
- 230000001360 synchronised effect Effects 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000000605 extraction Methods 0.000 claims description 3
- 239000000126 substance Substances 0.000 claims description 3
- 230000008859 change Effects 0.000 description 7
- 238000005070 sampling Methods 0.000 description 3
- 230000002776 aggregation Effects 0.000 description 2
- 238000004220 aggregation Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000003825 pressing Methods 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 230000002441 reversible effect Effects 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 230000002194 synthesizing effect 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/30—Analysis
- G01V1/301—Analysis for determining seismic cross-sections or geostructures
- G01V1/302—Analysis for determining seismic cross-sections or geostructures in 3D data cubes
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/624—Reservoir parameters
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
本发明涉及一种基于多次挤压小波变换的地震时频分析方法及系统,其特征在于,包括以下内容:1)选取待测三维地震数据体中某一剖面的若干道地震信号;2)对选取的地震信号分别进行傅里叶变换;3)根据得到的平均谱图,确定选取的各道地震信号所对应的主频区间,并在该主频区间内选取参考频率数值;4)对选取的各道地震信号分别作多次挤压小波变换;5)根据选取的参考频率数值和得到的时频谱,提取待测三维地震数据体中该剖面的三维常频率数据体;6)对提取出的三维常频率数据体进行沿层切片;7)进入步骤1),直至实现待测三维地震数据体中所有剖面三维常频率数据体的沿层切片,本发明可以广泛应用于地震勘探技术领域中。
Description
技术领域
本发明是关于一种基于多次挤压小波变换的地震时频分析方法及系统,属于地震勘探技术领域。
背景技术
受地下复杂介质的影响,地面上所接收到的地震信号是一种典型的非平稳信号。傅里叶变换只能提供整个时间上的频率信息,不能描述局部频率的变化。时频分析能够将一维的时间信号映射到一个二维的时间-频率域上,用于表示信号的频率成分随时间的变化情况。时频分析也逐渐成为地震信号处理中的一个重要方式,可以用时频分析工具来描述地震信号的局部变化特征,进而利用这些特征分析所对应的地下结构和储层。
传统的时频分析方法,例如短时傅里叶变换、小波变换和S变换,已经广泛应用于地震时频分析中,其中,短时傅里叶变换的缺点是窗函数的选择问题,一旦选定窗函数,则短时傅里叶变换的时间分辨率和频率分辨率即固定。因此,短时傅里叶变换对处理强非平稳性信号不够灵活。为解决短时傅里叶变换的不足,研究人员提出并完善了小波变换,小波变换既有短时傅里叶变换的功能,实现局部频率变化的描述,又具有多分辨率特性。但是,小波变换在低频时的时间分辨率和在高频时的频率分辨率较低。在1996年,研究人员提出了S变换,S变换结合了短时傅里叶变换和小波变换的优势,且克服其缺陷,S变换是一个线性、多分辨率、无损可逆的时频分析方法。然而,S变换的窗函数随频率变化的趋势是固定不变的,不能根据处理的数据的不同而改变。
可以看出,传统的时频分析方法受限于测不准原理,时间分辨率和频率分辨率不能同时达到最优,这往往导致得到的时频谱的能量是扩散的。因此,在实际应用中,选择或设计一个能量更加聚集的时频分析工具是一个重要的问题。
发明内容
针对上述问题,本发明的目的是提供一种能量更加聚集的基于多次挤压小波变换的地震时频分析方法及系统。
为实现上述目的,本发明采取以下技术方案:一种基于多次挤压小波变换的地震时频分析方法,包括以下内容:1)选取待测三维地震数据体中某一剖面的若干道地震信号;2)对选取的地震信号分别进行傅里叶变换,得到对应的平均谱图;3)根据得到的平均谱图,确定选取的各道地震信号所对应的主频区间,并在该主频区间内选取参考频率数值;4)对选取的各道地震信号分别作多次挤压小波变换,得到各道地震信号多次挤压小波变换的时频谱;5)根据选取的参考频率数值和得到的时频谱,提取待测三维地震数据体中该剖面的三维常频率数据体;6)对提取的三维常频率数据体进行沿层切片;7)进入所述步骤1)重新选取待测三维地震数据体中某一剖面的若干道地震信号,直至完成待测三维地震数据体中所有剖面三维常频率数据体的沿层切片。
进一步地,所述步骤4)的具体过程为:4.1)对地震信号f(t)进行小波变换,得到地震信号的时频表示W(a,b):
…
其中,δ(·)为单位冲激函数。
进一步地,所述步骤5)的具体过程为:5.1)根据选取的参考频率数值和得到的时频谱,提取出待测三维地震数据体中该剖面某道地震信号的常频率切线T[M](ω0,b),其中,ω0为选取的参考频率数值;5.2)重复所述步骤5.1),直至提取出待测三维地震数据体中该剖面各道地震信号的常频率切线,各道地震信号的常频率切线集合为待测三维地震数据体中该剖面的三维常频率数据体。
一种基于多次挤压小波变换的地震时频分析系统,包括:地震信号选取模块,用于选取待测三维地震数据体中某一剖面的若干道地震信号;傅里叶变换模块,用于对选取的地震信号分别进行傅里叶变换,得到对应的平均谱图;参考频率数值选取模块,用于根据得到的平均谱图,确定选取的各道地震信号所对应的主频区间,并在该主频区间内选取参考频率数值;多次挤压小波变换模块,用于对选取的各道地震信号分别作多次挤压小波变换,得到各道地震信号多次挤压小波变换的时频谱;三维常频率数据体提取模块,用于根据选取的参考频率数值和得到的时频谱,提取待测三维地震数据体中该剖面的三维常频率数据体;沿层切片模块,用于对提取出的三维常频率数据体进行沿层切片。
进一步地,所述多次挤压小波变换模块包括:小波变换单元,用于对地震信号进行小波变换,得到地震信号的时频表示;频率估计值计算单元,用于根据地震信号的时频表示,计算地震信号的频率估计值;时频谱计算单元,用于根据地震信号的频率估计值,计算地震信号多次挤压小波变换的时频谱。
一种计算机程序,包括计算机程序指令,其中,所述计算机程序指令被处理器执行时用于实现上述基于多次挤压小波变换的地震时频分析方法对应的步骤。
一种计算机可读存储介质,所述计算机可读存储介质上存储有计算机程序指令,其中,所述计算机程序指令被处理器执行时用于实现上述基于多次挤压小波变换的地震时频分析方法对应的步骤。
本发明由于采取以上技术方案,其具有以下优点:本发明首次在小波变换的框架内,引入多次挤压的瞬时频率估计算子,该变换经过对小波变换的时频系数沿频率方向进行重排,使得时频表示的能量更加聚集。对比于二阶同步挤压变换,多次挤压小波变换更加适合分析较强的非平稳性地震信号。在实际地震数据的时频分析应用中,能够更聚集、更准确地表征地下结构的边界,更准确地描述强频率调制信号的局部频率信息,进而能够检测更多潜在性的地下储层,可以广泛应用于地震勘探技术领域中。
附图说明
图1是本发明实施例中合成信号不同方法的时频谱示意图,其中,图1(a)为采用二阶挤压小波变换的时频谱,图1(b)为采用多次挤压小波变换的时频谱;
图2是本发明实施例中某油田的二维地震剖面图;
图3是图2中部分地震信号的40Hz常频率数据体的沿层切片示意图,其中,图3(a)为采用二阶挤压小波变换的40Hz常频率数据体的沿层切片,图3(b)为采用多次挤压小波变换的40Hz常频率数据体的沿层切片。
具体实施方式
以下结合附图来对本发明进行详细的描绘。然而应当理解,附图的提供仅为了更好地理解本发明,它们不应该理解成对本发明的限制。
实施例一
本实施例提供一种基于多次挤压小波变换的地震时频分析方法,包括以下步骤:
1)选取待测三维地震数据体中某一剖面的若干道地震信号。
2)对选取的地震信号分别进行傅里叶变换,得到对应的平均谱图。
3)根据得到的平均谱图,确定选取的各道地震信号所对应的主频区间(例如为20-80Hz),并在该主频区间内选取一低频率、中频率和高频率作为下述提取三维常频率数据体的参考频率数值,例如选取30、50和70Hz,其中,低频、中频和高频可以根据实际情况进行选定,具体过程在此不多做赘述。
4)对选取的各道地震信号分别作多次挤压小波变换,得到各道地震信号多次挤压小波变换的时频谱,具体为:
4.1)对地震信号f(t)进行小波变换,得到地震信号的时频表示W(a,b):
其中,ψ(t)为母小波函数;a和b分别为尺度因子和时间变量。
…
其中,δ(·)为单位冲激函数。
5)根据选取的参考频率数值和得到的时频谱,提取待测三维地震数据体中该剖面的三维常频率数据体,具体为:
5.1)根据选取的参考频率数值和得到的时频谱,提取出待测三维地震数据体中该剖面某道地震信号的常频率切线T[M](ω0,b),其中,ω0为步骤3)中选取的参考频率数值,采用常频率切线T[M](ω0,b)替代地震信号在三维地震数据体中的位置。
5.2)重复步骤5.1),直至提取出待测三维地震数据体中该剖面各道地震信号的常频率切线,各道地震信号的常频率切线集合即为待测三维地震数据体中该剖面的三维常频率数据体。
6)对提取出的三维常频率数据体进行沿层切片,以供地质专业解释人员参考分析。
7)进入步骤1)重新选取待测三维地震数据体中某一剖面的若干道地震信号,直至完成待测三维地震数据体中所有剖面三维常频率数据体的沿层切片,完成待测三维地震数据体的地震时频分析。
下面通过具体实施例详细说明本发明的基于多次挤压小波变换的地震时频分析方法:
A)合成模型数据
为验证多次挤压小波变换提高时频谱能量聚集性的有效性,本实施例选取一个合成信号,其数学表达式为:
f(t)=f1(t)+f2(t)
其中:
如图1所示,为二阶同步挤压变换和多次挤压小波变换的时频结果,对比图1(a)和(b),可以发现在频率缓变的地方,两种不同的时频分析方法均能够很精确地描绘信号的频率变化规律,但是,在频率变化较快的地方,多次挤压小波变换的时频聚集性更好,能够更准确地刻画信号的时频特性。
B)实际地震资料
基于理论分析结果,进一步将多次挤压小波变换应用至地震时频分析中。如图2所示,为某油田的二维地震资料,该剖面包括361道,每道均有751个采样点,采样间隔为0.002s。采用本发明方法,选择201道、每道501个采样点的子区域进行分析,如图3所示。由于多次挤压小波变换能使得时频谱的能量更加聚集,所以能够更准确地观测储层的边界(箭头指示处),且能够检测到更多潜在可能性的结构信息。
实施例二
本实施例提供一种基于多次挤压小波变换的地震时频分析系统,包括:
地震信号选取模块,用于选取待测三维地震数据体中某一剖面的若干道地震信号;
傅里叶变换模块,用于对选取的地震信号分别进行傅里叶变换,得到对应的平均谱图;
参考频率数值选取模块,用于根据得到的平均谱图,确定选取的各道地震信号所对应的主频区间,并在该主频区间内选取参考频率数值;
多次挤压小波变换模块,用于对选取的各道地震信号分别作多次挤压小波变换,得到各道地震信号多次挤压小波变换的时频谱;
三维常频率数据体提取模块,用于根据选取的参考频率数值和得到的时频谱,提取待测三维地震数据体中该剖面的三维常频率数据体;
沿层切片模块,用于对提取出的三维常频率数据体进行沿层切片。
在一个优选的实施例中,多次挤压小波变换模块包括:小波变换单元,用于对地震信号进行小波变换,得到地震信号的时频表示;频率估计值计算单元,用于根据地震信号的时频表示,计算地震信号的频率估计值;时频谱计算单元,用于根据地震信号的频率估计值,计算地震信号多次挤压小波变换的时频谱。
实施例三
本实施例提供一种计算机程序,包括计算机程序指令,其中,计算机程序指令被处理器执行时用于实现上述基于多次挤压小波变换的地震时频分析方法对应的步骤。
实施例四
本实施例提供一种计算机可读存储介质,计算机可读存储介质上存储有计算机程序指令,其中,计算机程序指令被处理器执行时用于实现上述基于多次挤压小波变换的地震时频分析方法对应的步骤。
上述各实施例仅用于说明本发明,其中各部件的结构、连接方式和制作工艺等都是可以有所变化的,凡是在本发明技术方案的基础上进行的等同变换和改进,均不应排除在本发明的保护范围之外。
Claims (7)
1.一种基于多次挤压小波变换的地震时频分析方法,其特征在于,包括以下内容:
1)选取待测三维地震数据体中某一剖面的若干道地震信号;
2)对选取的地震信号分别进行傅里叶变换,得到对应的平均谱图;
3)根据得到的平均谱图,确定选取的各道地震信号所对应的主频区间,并在该主频区间内选取参考频率数值;
4)对选取的各道地震信号分别作多次挤压小波变换,得到各道地震信号多次挤压小波变换的时频谱;
5)根据选取的参考频率数值和得到的时频谱,提取待测三维地震数据体中该剖面的三维常频率数据体;
6)对提取的三维常频率数据体进行沿层切片;
7)进入所述步骤1)重新选取待测三维地震数据体中某一剖面的若干道地震信号,直至完成待测三维地震数据体中所有剖面三维常频率数据体的沿层切片。
3.如权利要求2所述的一种基于多次挤压小波变换的地震时频分析方法,其特征在于,所述步骤5)的具体过程为:
5.1)根据选取的参考频率数值和得到的时频谱,提取出待测三维地震数据体中该剖面某道地震信号的常频率切线T[M](ω0,b),其中,ω0为选取的参考频率数值;
5.2)重复所述步骤5.1),直至提取出待测三维地震数据体中该剖面各道地震信号的常频率切线,各道地震信号的常频率切线集合为待测三维地震数据体中该剖面的三维常频率数据体。
4.一种基于多次挤压小波变换的地震时频分析系统,其特征在于,包括:
地震信号选取模块,用于选取待测三维地震数据体中某一剖面的若干道地震信号;
傅里叶变换模块,用于对选取的地震信号分别进行傅里叶变换,得到对应的平均谱图;
参考频率数值选取模块,用于根据得到的平均谱图,确定选取的各道地震信号所对应的主频区间,并在该主频区间内选取参考频率数值;
多次挤压小波变换模块,用于对选取的各道地震信号分别作多次挤压小波变换,得到各道地震信号多次挤压小波变换的时频谱;
三维常频率数据体提取模块,用于根据选取的参考频率数值和得到的时频谱,提取待测三维地震数据体中该剖面的三维常频率数据体;
沿层切片模块,用于对提取出的三维常频率数据体进行沿层切片。
5.如权利要求4所述的一种基于多次挤压小波变换的地震时频分析系统,其特征在于,所述多次挤压小波变换模块包括:
小波变换单元,用于对地震信号进行小波变换,得到地震信号的时频表示;
频率估计值计算单元,用于根据地震信号的时频表示,计算地震信号的频率估计值;
时频谱计算单元,用于根据地震信号的频率估计值,计算地震信号多次挤压小波变换的时频谱。
6.一种计算机程序,其特征在于,包括计算机程序指令,其中,所述计算机程序指令被处理器执行时用于实现权利要求1-3中任一项所述的基于多次挤压小波变换的地震时频分析方法对应的步骤。
7.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质上存储有计算机程序指令,其中,所述计算机程序指令被处理器执行时用于实现权利要求1-3中任一项所述的基于多次挤压小波变换的地震时频分析方法对应的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010354535.4A CN111366978B (zh) | 2020-04-29 | 2020-04-29 | 一种基于多次挤压小波变换的地震时频分析方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010354535.4A CN111366978B (zh) | 2020-04-29 | 2020-04-29 | 一种基于多次挤压小波变换的地震时频分析方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111366978A true CN111366978A (zh) | 2020-07-03 |
CN111366978B CN111366978B (zh) | 2021-08-17 |
Family
ID=71211055
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010354535.4A Active CN111366978B (zh) | 2020-04-29 | 2020-04-29 | 一种基于多次挤压小波变换的地震时频分析方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111366978B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113723200A (zh) * | 2021-08-03 | 2021-11-30 | 同济大学 | 一种非平稳信号的时频谱结构特征提取方法 |
CN115390140A (zh) * | 2022-08-24 | 2022-11-25 | 中国科学院地质与地球物理研究所 | 一种地震属性体提取的方法及系统 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104880730A (zh) * | 2015-03-27 | 2015-09-02 | 西安交通大学 | 基于Synchrosqueezing变换的地震资料时频分析和衰减估计方法 |
CN107132576A (zh) * | 2017-07-05 | 2017-09-05 | 西安交通大学 | 基于二阶挤压小波变换的地震资料时频分析方法 |
CN109239780A (zh) * | 2017-07-10 | 2019-01-18 | 中国石油化工股份有限公司 | 基于同步挤压小波变换去除面波的方法 |
-
2020
- 2020-04-29 CN CN202010354535.4A patent/CN111366978B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104880730A (zh) * | 2015-03-27 | 2015-09-02 | 西安交通大学 | 基于Synchrosqueezing变换的地震资料时频分析和衰减估计方法 |
CN107132576A (zh) * | 2017-07-05 | 2017-09-05 | 西安交通大学 | 基于二阶挤压小波变换的地震资料时频分析方法 |
CN109239780A (zh) * | 2017-07-10 | 2019-01-18 | 中国石油化工股份有限公司 | 基于同步挤压小波变换去除面波的方法 |
Non-Patent Citations (2)
Title |
---|
GANG YU ET AL.: "Multisynchrosqueezing Transform", 《IEEE TRANSACTIONS ON INDUSTRIAL ELECTRONICS》 * |
XIANGXIANG ZHU ET AL.: "Multiple squeezes from adaptive chirplet transform", 《SIGNAL PROCESSING》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113723200A (zh) * | 2021-08-03 | 2021-11-30 | 同济大学 | 一种非平稳信号的时频谱结构特征提取方法 |
CN113723200B (zh) * | 2021-08-03 | 2024-01-12 | 同济大学 | 一种非平稳信号的时频谱结构特征提取方法 |
CN115390140A (zh) * | 2022-08-24 | 2022-11-25 | 中国科学院地质与地球物理研究所 | 一种地震属性体提取的方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN111366978B (zh) | 2021-08-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Chen | Dip-separated structural filtering using seislet transform and adaptive empirical mode decomposition based dip filter | |
US7764573B2 (en) | Method of processing seismic data and method of producing a mineral hydrocarbon fluid and a computer program product | |
CN111366978B (zh) | 一种基于多次挤压小波变换的地震时频分析方法及系统 | |
CN107132576B (zh) | 基于二阶同步挤压小波变换的地震资料时频分析方法 | |
CN106597532A (zh) | 一种结合井资料与层位资料的叠前地震数据频带拓展方法 | |
CN111505716B (zh) | 一种基于时间同步抽取广义Chirplet变换的地震时频分析方法 | |
CN107272062B (zh) | 一种数据驱动的地下介质q场估计方法 | |
CN102053273A (zh) | 一种对地震波信号进行反q滤波的方法 | |
CN109164489A (zh) | 一种基于vmd与tk能量算子的地震流体预测方法 | |
CN109031422A (zh) | 一种基于CEEMDAN与Savitzky-Golay滤波的地震信号噪声抑制方法 | |
CN108845357B (zh) | 一种基于同步挤压小波变换估计地层等效品质因子的方法 | |
CN107179550B (zh) | 一种数据驱动的地震信号零相位反褶积方法 | |
CN113608259B (zh) | 一种基于iceemdan约束广义s变换的地震薄层检测方法 | |
CN113777650B (zh) | 一种基于混合范数和小波变换的稀疏时频谱分解方法、装置、设备及存储介质 | |
CN104730576A (zh) | 基于Curvelet变换的地震信号去噪方法 | |
WO2005017564A1 (en) | Dip value in seismic images | |
Duchesne et al. | Did you smooth your well logs the right way for seismic interpretation? | |
CN104422956B (zh) | 一种基于稀疏脉冲反演的高精度地震谱分解方法 | |
CN102854530B (zh) | 基于对数时频域双曲平滑的动态反褶积方法 | |
CN107703548A (zh) | 基于沉积物品质因子和回波损失级曲线峰谷的浅地层层界划分方法 | |
CN104182944A (zh) | 基于曲波与小波变换相串联的光学图像去噪方法 | |
CN104932008B (zh) | 补偿j变换的复时‑频谱提高地震剖面分辨率的方法 | |
Guo et al. | Seismic data denoising under the morphological component analysis framework by dictionary learning | |
CN114200522B (zh) | 深度域地震子波提取方法和装置、存储介质及电子设备 | |
CA3203426A1 (en) | Calculation method and device for interval transit time, and storage medium |
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 |