CN113917522B - 用于指导采集观测系统设计的地震正演方法 - Google Patents
用于指导采集观测系统设计的地震正演方法 Download PDFInfo
- Publication number
- CN113917522B CN113917522B CN202010661116.5A CN202010661116A CN113917522B CN 113917522 B CN113917522 B CN 113917522B CN 202010661116 A CN202010661116 A CN 202010661116A CN 113917522 B CN113917522 B CN 113917522B
- Authority
- CN
- China
- Prior art keywords
- wave number
- scheme
- observation system
- target area
- scattering wave
- 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
- 238000013461 design Methods 0.000 title claims abstract description 37
- 238000000034 method Methods 0.000 title claims abstract description 30
- 239000013598 vector Substances 0.000 claims abstract description 69
- 238000013508 migration Methods 0.000 claims abstract description 20
- 230000005012 migration Effects 0.000 claims abstract description 20
- 238000004364 calculation method Methods 0.000 claims abstract description 14
- 230000005284 excitation Effects 0.000 claims description 23
- 238000009826 distribution Methods 0.000 claims description 4
- 230000015572 biosynthetic process Effects 0.000 claims description 3
- 238000003384 imaging method Methods 0.000 abstract description 6
- 238000012545 processing Methods 0.000 abstract description 6
- 238000004519 manufacturing process Methods 0.000 abstract description 5
- 230000008569 process Effects 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000012938 design process Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
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/282—Application of seismic models, synthetic seismograms
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/16—Receiving elements for seismic signals; Arrangements or adaptations of receiving elements
- G01V1/20—Arrangements of receiving elements, e.g. geophone pattern
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/10—Aspects of acoustic signal generation or detection
- G01V2210/12—Signal generation
- G01V2210/129—Source location
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/10—Aspects of acoustic signal generation or detection
- G01V2210/14—Signal detection
- G01V2210/142—Receiver location
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (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，从步骤3计算得到的散射波数向量组中抽取相应的散射波数向量，即可得到该方案的散射波数向量组；步骤5，应用傅里叶变换，得到波数域目标模型数据；步骤6，计算该方案的目标区域的叠前深度偏移剖面。该用于指导采集观测系统设计的地震正演方法可以高效快速的获得叠前深度偏移剖面，避免了常规技术方法中先生成叠前炮记录，然后再处理偏移得到成像剖面的步骤，可以用于地震采集观测系统设计实际生产。
Description
技术领域
本发明涉及油田开发技术领域，特别是涉及到一种用于指导采集观测系统设计的地震正演方法。
背景技术
地震采集观测系统设计是地震勘探中的重要步骤。对地下目标区域，需要合理的设计观测系统，以实现对目标区域的高质量成像。因此，为了检验评价观测系统设计方案的合理性，一般采用地震正演技术进行数值模拟，进而得到叠前成像剖面。这种技术方法的思路一般为，对形成的观测系统设计方案，结合地质模型，进行(基于波动方程或者射线追踪的)地震正演模拟，得到(二维或三维)地震炮记录；然后对地震炮记录进行处理、偏移等操作，得到叠前偏移剖面；通过分析目标区域的偏移剖面，对观测系统设计方案进行确定或分析修改。
但是这种技术方法存在效率难以满足生产需求的问题。一方面，观测系统设计一般在短时间内需要设计多套方案，或者需要经过多次修改完善，因此需要地震正演的次数较多，并且频率较高；另一方面，为了更充分的获得地下复杂地质构造信息，目前的观测系统方案大多设计为宽方位、大偏移距，因此导致单次地震正演的计算量巨大，即时采用一些效率优化方案，其计算周期仍然较长。而且正演得到的炮记录，还需要经过常规处理、偏移等流程才能得到成像剖面，这一过程同样十分耗时。总之，该技术方案难以满足目前观测系统设计对地震正演的高频次使用需求，缺乏可行性。
现有的技术方法中正演模拟过程和炮记录数据处理偏移过程均十分耗时，周期较长，难以满足观测系统设计生产中的对地震正演的高频次使用需求。
为此我们发明了一种新的用于指导采集观测系统设计的地震正演方法，解决了以上技术问题。
发明内容
本发明的目的是提供一种可以高效快速的获得叠前深度偏移剖面的用于指导采集观测系统设计的地震正演方法。
本发明的目的可通过如下技术措施来实现：用于指导采集观测系统设计的地震正演方法，该用于指导采集观测系统设计的地震正演方法包括：步骤1，设计初始的观测系统方案；步骤2，建立覆盖全部观测系统方案区域的地质模型；步骤3，计算方案目标区域中全部网格点上的散射波数向量组；步骤4，从步骤3计算得到的散射波数向量组中抽取相应的散射波数向量，即可得到该方案的散射波数向量组；步骤5，应用傅里叶变换，得到波数域目标模型数据；步骤6，计算该方案的目标区域的叠前深度偏移剖面。
本发明的目的还可通过如下技术措施来实现：
在步骤1中，选择对多个观测系统方案进行分组，将某一方案和对该方案的激发或接收进行加密或粗化的方案分在一组中。
在步骤2中，建立覆盖全部观测系统方案区域的地质模型，并确定目标区域；对地质模型进行网格划分。
在步骤3中，对每一个观测系统方案分组，选取分组中面元最小的方案，计算该方案目标区域中全部网格点上的散射波数向量组。
在步骤3中，散射波数向量组的计算使用以下方式：通过求解程函方程计算地震波由激发点到地下目标区域网格点的走时和入射角度，得到入射波向量，向量的大小等于波数，方向为入射角度；同样计算由地下目标区域网格点到接收点的出射波向量；通过公式(1)计算得到该网格点对应于一组激发点和接收点的散射波数向量；计算该网格点对应的全部激发点和接收点的散射波数向量，得到散射波数向量组；对目标区域中全部网格点，应用以上步骤，即可完成计算；
公式(1)为：
其中ω为激发源子波频率，为入射走时梯度，/>为入射走时梯度，cn为网格点附近的地层速度，ks为沿着激发点到网格点的射线路径上的波数向量，kr为沿着网格点到接收点的射线路径上的波数向量，二者的向量差即为该点的散射波数向量。
在步骤4中，对于分组中的其他方案，根据该方案的激发点和接收点分布情况，从上一步计算得到的散射波数向量组中抽取相应的散射波数向量，即可得到该方案的散射波数向量组；由此可得到全部设计方案的散射波数向量组。
在步骤6中，对于每一个设计方案，将波数域目标模型数据与该方案的散射波数向量组相乘，并应用傅里叶反变换，即可得到该方案的目标区域的叠前深度偏移剖面。
该用于指导采集观测系统设计的地震正演方法还包括，在步骤6之后，通过对叠前深度偏移剖面进行分析，确定或改进观测系统设计方案。
本发明中的用于指导采集观测系统设计的地震正演方法，可以高效快速的获得叠前深度偏移剖面，避免了常规技术方法中先生成叠前炮记录，然后再处理偏移得到成像剖面的步骤，可以用于地震采集观测系统设计实际生产。本发明与现有技术相比，具有以下优点：
第一、本发明避免了传统技术方法中先生成叠前炮记录，然后在处理偏移得到成像剖面的步骤，效率高，可满足实际生产需求。
第二、本发明规避了观测系统设计过程中大量的重复正演计算，提高了应用效率。
附图说明
图1为本发明的用于指导采集观测系统设计的地震正演方法的一具体实施例的流程图；
图2为本发明的一具体实施例中使用的模型及在模型中选定目标区域；
图3为本发明的一具体实施例中选取的目标区域的放大图；
图4为本发明的一具体实施例中目标区域中某一点(记为A)的散射波向量组的极坐标显示；
图5为本发明的一具体实施例中目标区域中A点处的模型转化到波数域的波数域振幅谱；
图6为本发明的一具体实施例中根据(a)-(e)观测系统方案得到的目标区域的叠前深度偏移剖面。
具体实施方式
为使本发明的上述和其他目的、特征和优点能更明显易懂，下文特举出较佳实施例，并配合附图所示，作详细说明如下。
如图1所示，图1为本发明的用于指导采集观测系统设计的地震正演方法的流程图。
步骤101，设计形成初始的一个或多个观测系统方案。为了提高效率，可以选择对多个观测系统方案进行分组，将某一方案和对该方案的激发或接收进行加密或粗化的方案分在一组中。
步骤102，建立覆盖全部观测系统方案区域的地质模型，并确定目标区域。对地质模型进行网格划分。
步骤103，对每一个观测系统方案分组，选取分组中面元最小的方案，计算该方案目标区域中全部网格点上的散射波数向量组。散射波数向量组的计算可以使用(但不限于)以下方式：通过求解程函方程计算地震波由激发点到地下目标区域网格点的走时和入射角度，得到入射波向量，向量的大小等于波数，方向为入射角度；同样可计算由地下目标区域网格点到接收点的出射波向量；通过公式(1)计算得到该网格点对应于一组激发点和接收点的散射波数向量；计算该网格点对应的全部激发点和接收点的散射波数向量，得到散射波数向量组；对目标区域中全部网格点，应用以上步骤，即可完成计算。
其中ω为激发源子波频率，为入射走时梯度，/>为入射走时梯度，cn为网格点附近的地层速度，ks为沿着激发点到网格点的射线路径上的波数向量，kr为沿着网格点到接收点的射线路径上的波数向量，二者的向量差即为该点的散射波数向量。
步骤104，对于分组中的其他方案，根据该方案的激发点和接收点分布情况，从上一步计算得到的散射波数向量组中抽取相应的散射波数向量，即可得到该方案的散射波数向量组。由此可得到全部设计方案的散射波数向量组。采用分组和抽取的方式，省去大量重复的计算量。设计方案的分组至此步骤可以取消。
步骤105，对地质模型的目标区域应用傅里叶变换，得到波数域目标模型数据。
步骤106，对于每一个设计方案，将波数域目标模型数据与该方案的散射波数向量组相乘，并应用傅里叶反变换，即可得到该方案的目标区域的叠前深度偏移剖面。
步骤107，通过对叠前深度偏移剖面进行分析，确定或改进观测系统设计方案。
在应用本发明的一具体实施例中，以经典Marmousi模型为例，说明本发明的具体实施流程。
(1)初步设计4套观测系统方案，具体参数如下：每炮的接收道数为炮点右侧100道接收，道间距为10米，最小偏移距为10米，最大偏移距1000米；起始炮点位置为0米，4套方案的炮间距分别为(a)100米、(b)50米、(c)30米、(d)20米、(e)10米，相应的覆盖次数分别为(a)5次、(b)10次、(c)17次、(d)25次、(e)50次。
(2)使用经典的Marmousi模型，如图2所示
选定如图2中黑框中所示区域为目标区域。使用2米乘以2米的网格将地质模型进行网格划分，如图3所示。
(3)对于方案(e)，设置激发源子波为25赫兹的Ricker子波，应用基于求解程函方程的射线追踪算法，计算目标区域中全部网格点上的散射波数向量组，如图4所示。
(4)对于其他方案，根据该方案的激发点和接收点分布情况，从方案(e)的散射波数向量组中抽取相应的散射波数向量，即可得到该方案的散射波数向量组。
(5)将地质模型目标区域应用傅里叶变换，转换到波数域，如图5所示(以方案(e)为例)。
(6)将波数域目标模型数据与该方案的散射波数向量组相乘，并傅里叶反变换，即可得到目标区域的叠前深度偏移剖面，如图6所示。
(7)通过对叠前深度偏移剖面进行分析，综合成像质量和经济实用高效等因素，建议采用20米到30米之间的炮间距进行地震采集。
Claims (6)
1.用于指导采集观测系统设计的地震正演方法，其特征在于，该用于指导采集观测系统设计的地震正演方法包括：
步骤1，设计初始的观测系统方案；
步骤2，建立覆盖全部观测系统方案区域的地质模型；
步骤3，计算方案目标区域中全部网格点上的散射波数向量组；
步骤4，从步骤3计算得到的散射波数向量组中抽取相应的散射波数向量，即可得到该方案的散射波数向量组；
步骤5，应用傅里叶变换，得到波数域目标模型数据；
步骤6，计算该方案的目标区域的叠前深度偏移剖面；
在步骤3中，对每一个观测系统方案分组，选取分组中面元最小的方案，计算该方案目标区域中全部网格点上的散射波数向量组；
散射波数向量组的计算使用以下方式：通过求解程函方程计算地震波由激发点到地下目标区域网格点的走时和入射角度，得到入射波向量，向量的大小等于波数，方向为入射角度；同样计算由地下目标区域网格点到接收点的出射波向量；通过公式(1)计算得到该网格点对应于一组激发点和接收点的散射波数向量；计算该网格点对应的全部激发点和接收点的散射波数向量，得到散射波数向量组；对目标区域中全部网格点，应用以上步骤，即可完成计算；
公式(1)为：
其中ω为激发源子波频率，为入射走时梯度，/>为出射走时梯度，cn为网格点附近的地层速度，ks为沿着激发点到网格点的射线路径上的波数向量，kr为沿着网格点到接收点的射线路径上的波数向量，二者的向量差即为该点的散射波数向量。
2.根据权利要求1所述的用于指导采集观测系统设计的地震正演方法，其特征在于，在步骤1中，选择对多个观测系统方案进行分组，将某一方案和对该方案的激发或接收进行加密或粗化的方案分在一组中。
3.根据权利要求1所述的用于指导采集观测系统设计的地震正演方法，其特征在于，在步骤2中，建立覆盖全部观测系统方案区域的地质模型，并确定目标区域；对地质模型进行网格划分。
4.根据权利要求1所述的用于指导采集观测系统设计的地震正演方法，其特征在于，在步骤4中，对于分组中的其他方案，根据该方案的激发点和接收点分布情况，从上一步计算得到的散射波数向量组中抽取相应的散射波数向量，即可得到该方案的散射波数向量组；由此可得到全部设计方案的散射波数向量组。
5.根据权利要求1所述的用于指导采集观测系统设计的地震正演方法，其特征在于，在步骤6中，对于每一个设计方案，将波数域目标模型数据与该方案的散射波数向量组相乘，并应用傅里叶反变换，即可得到该方案的目标区域的叠前深度偏移剖面。
6.根据权利要求1所述的用于指导采集观测系统设计的地震正演方法，其特征在于，该用于指导采集观测系统设计的地震正演方法还包括，在步骤6之后，通过对叠前深度偏移剖面进行分析，确定或改进观测系统设计方案。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010661116.5A CN113917522B (zh) | 2020-07-10 | 2020-07-10 | 用于指导采集观测系统设计的地震正演方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010661116.5A CN113917522B (zh) | 2020-07-10 | 2020-07-10 | 用于指导采集观测系统设计的地震正演方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113917522A CN113917522A (zh) | 2022-01-11 |
| CN113917522B true CN113917522B (zh) | 2024-03-19 |
Family
ID=79232060
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010661116.5A Active CN113917522B (zh) | 2020-07-10 | 2020-07-10 | 用于指导采集观测系统设计的地震正演方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113917522B (zh) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101625417A (zh) * | 2008-07-08 | 2010-01-13 | 中国石油集团东方地球物理勘探有限责任公司 | 一种优化垂直地震剖面观测系统设计的方法 |
| CN102313900A (zh) * | 2011-09-01 | 2012-01-11 | 中国石油集团川庆钻探工程有限公司 | 三维地震采集观测系统的激发位置确定方法 |
| CN104050359A (zh) * | 2014-05-30 | 2014-09-17 | 中国石油大学(华东) | 一种基于三维观测系统排列片数据分割的正演模拟方法 |
| CN104090297A (zh) * | 2014-06-14 | 2014-10-08 | 吉林大学 | 一种优化地震采集观测系统的逆向照明方法 |
| CN107490808A (zh) * | 2017-06-14 | 2017-12-19 | 中国石油化工股份有限公司 | 一种高可靠性地震勘探观测系统的建立方法 |
| CN109212589A (zh) * | 2017-07-06 | 2019-01-15 | 中国石油化工股份有限公司 | 一种基于gpu/cpu协同并行的地震采集观测系统设计方法 |
| CN110888159A (zh) * | 2019-11-15 | 2020-03-17 | 西安理工大学 | 基于角度分解与波场分离的弹性波全波形反演方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US10295683B2 (en) * | 2016-01-05 | 2019-05-21 | Schlumberger Technology Corporation | Amplitude inversion on partitioned depth image gathers using point spread functions |
-
2020
- 2020-07-10 CN CN202010661116.5A patent/CN113917522B/zh active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101625417A (zh) * | 2008-07-08 | 2010-01-13 | 中国石油集团东方地球物理勘探有限责任公司 | 一种优化垂直地震剖面观测系统设计的方法 |
| CN102313900A (zh) * | 2011-09-01 | 2012-01-11 | 中国石油集团川庆钻探工程有限公司 | 三维地震采集观测系统的激发位置确定方法 |
| CN104050359A (zh) * | 2014-05-30 | 2014-09-17 | 中国石油大学(华东) | 一种基于三维观测系统排列片数据分割的正演模拟方法 |
| CN104090297A (zh) * | 2014-06-14 | 2014-10-08 | 吉林大学 | 一种优化地震采集观测系统的逆向照明方法 |
| CN107490808A (zh) * | 2017-06-14 | 2017-12-19 | 中国石油化工股份有限公司 | 一种高可靠性地震勘探观测系统的建立方法 |
| CN109212589A (zh) * | 2017-07-06 | 2019-01-15 | 中国石油化工股份有限公司 | 一种基于gpu/cpu协同并行的地震采集观测系统设计方法 |
| CN110888159A (zh) * | 2019-11-15 | 2020-03-17 | 西安理工大学 | 基于角度分解与波场分离的弹性波全波形反演方法 |
Non-Patent Citations (6)
| Title |
|---|
| Finite Difference Simulations of Seismic Scattering: Implications for the Propagation of Short-Period Seismic Waves in the Crust and Models of Crustal Heterogeneity;A Frankel 等;Journal of Geophysical Research: Solid Earth;第91卷(第B6期);第6465-6489页 * |
| 三维叠前傅里叶变换法模型并行算法;吴艳辉 等;石油地球物理勘探;第41卷(第01期);第8-12页 * |
| 地震散射波场特征的数值模拟研究;尹军杰;中国博士学位论文全文数据库 基础科学辑(第06期);全文 * |
| 基于GPU 集群的大规模三维有限差分正演模拟并行策略;廉西猛 等;物探与化探;第39卷(第03期);第615-620页 * |
| 基于波动方程正演的地震观测系统设计;李万万;;石油地球物理勘探;第43卷(第02期);第134-141页 * |
| 高密度三维地震观测系统设计技术与应用;崔庆辉 等;石油物探;第59卷(第01期);第12-22页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113917522A (zh) | 2022-01-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108037526B (zh) | 基于全波波场vsp/rvsp地震资料的逆时偏移方法 | |
| CN111158049B (zh) | 一种基于散射积分法的地震逆时偏移成像方法 | |
| CN103238158B (zh) | 利用互相关目标函数进行的海洋拖缆数据同时源反演 | |
| CN102707316B (zh) | 复杂地质目标的2.5次三维地震勘探方法 | |
| CN110058303B (zh) | 声波各向异性逆时偏移混合方法 | |
| CN111766628A (zh) | 一种预条件的时间域弹性介质多参数全波形反演方法 | |
| US11397273B2 (en) | Full waveform inversion in the midpoint-offset domain | |
| CN101545986A (zh) | 基于最大能量旅行时计算的三维积分叠前深度偏移方法 | |
| CN113534243B (zh) | 一种被动源Marchenko成像方法及系统 | |
| CN101021568A (zh) | 基于最大能量旅行时计算的三维积分叠前深度偏移方法 | |
| CN106154315A (zh) | 一种对称均匀的高密度地震采集方法 | |
| CN112305612A (zh) | 高分辨率复谱分解时频空间域振幅随偏移距变化校正方法 | |
| CN105573963A (zh) | 一种电离层水平不均匀结构重构方法 | |
| CN106353798A (zh) | 多分量联合高斯束叠前逆时偏移成像方法 | |
| CN110850469A (zh) | 一种基于克希霍夫积分解的地震槽波深度偏移的成像方法 | |
| CN108508481B (zh) | 一种纵波转换波地震数据时间匹配的方法、装置及系统 | |
| CN113534259A (zh) | 一种可控震源高效采集实时叠前时间偏移成像方法 | |
| CN113917522B (zh) | 用于指导采集观测系统设计的地震正演方法 | |
| CN111123361B (zh) | 垂直地震剖面地震数据规则化重建方法及装置、存储介质 | |
| CN106291676A (zh) | 一种基于匹配追踪算法的地震数据重构方法 | |
| CN108363739A (zh) | 一种基于稀疏采集的地震资料高低频拓展方法 | |
| CN107807392A (zh) | 一种自适应抗频散的分块时空双变逆时偏移方法 | |
| CN111665550A (zh) | 地下介质密度信息反演方法 | |
| CN117826248B (zh) | 基于多尺度观测背景噪声聚束的面波频散提取方法 | |
| Jiaotong et al. | Time-varient slip sweep excitation time-distance optimization technology |
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 |