CN107300717B - 一种三维空间的方位偏移距道集匀化的方法 - Google Patents
一种三维空间的方位偏移距道集匀化的方法 Download PDFInfo
- Publication number
- CN107300717B CN107300717B CN201610236671.7A CN201610236671A CN107300717B CN 107300717 B CN107300717 B CN 107300717B CN 201610236671 A CN201610236671 A CN 201610236671A CN 107300717 B CN107300717 B CN 107300717B
- Authority
- CN
- China
- Prior art keywords
- trace gather
- point
- gather
- offset distance
- cdp
- 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 35
- 238000007619 statistical method Methods 0.000 claims abstract description 4
- 230000003068 static effect Effects 0.000 claims description 6
- 238000001914 filtration Methods 0.000 claims description 3
- 238000011084 recovery Methods 0.000 claims description 3
- 238000000265 homogenisation Methods 0.000 abstract description 6
- 238000010586 diagram Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000007781 pre-processing Methods 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
-
- 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/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
-
- 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/50—Corrections or adjustments related to wave propagation
- G01V2210/53—Statics correction, e.g. weathering layer or transformation to a datum
-
- 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/63—Seismic attributes, e.g. amplitude, polarity, instant phase
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道集进行统计分析,寻找最大覆盖的CDP点;对最大覆盖次数的单个CDP点的道集进行分方位分偏移距统计;对最大覆盖次数的单个CDP点的道集进行宏面元匀化;对所有叠前地震道集依次进行分方位分偏移距面元匀化;获得分方位分偏移距三维匀化数据的道集。本发明基于三维空间的方位角度道集匀化技术,以此得到一种相对方位角度比较均匀的叠前方位角度(偏移距)道集,可以为叠前AVAZ分析以及叠前各向异性反演提供较好的道集质量。
Description
技术领域
本发明涉及地震勘探技术,是一种对叠前方位偏移距道集进行自动匀化,获得较为均一化的方位偏移距道集的技术。
背景技术
目前的叠前AVAZ各向异性分析以及叠前AVAZ各向异性道集反演中,都需要利用高质量的叠前方位角度(偏移距)道集,道集质量直接决定了反演成果的好坏;理想的AVAZ反演道集应该是方位-角度(偏移距)道集在各个方位和角度(偏移距)都是均匀的,但是由于采集观测系统的限制,实际地震勘探中得到的叠前方位角度道集都存在较大的不均一性,因此,本领域亟需开发一种基于三维空间的方位偏移距度道集匀化技术,以此得到一种相对方位角度比较均匀的叠前方位偏移距度道集。
发明内容
为了解决现有技术中叠前方位角度(偏移距)道集存在较大的不均一性的问题,本发明提供一种基于三维空间的方位偏移距道集匀化技术,以此得到一种相对方位角度比较均匀的叠前方位偏移距道集,可以为叠前AVAZ分析以及叠前各向异性反演提供较好的道集质量。
根据本发明的一个方面,提供一种三维空间的方位偏移距道集匀化的方法,该方法包括如下步骤:抽取宏面元道集;对叠前CMP道集进行统计分析,寻找最大覆盖的CDP点;对最大覆盖次数的单个CDP点的道集进行分方位分偏移距统计;对最大覆盖次数的单个CDP点的道集进行宏面元匀化;对所有叠前地震道集依次进行分方位分偏移距面元匀化;获得分方位分偏移距三维匀化数据的道集。
进一步地,该方法还包括在抽取宏面元道集之前对叠前地震道集进行预处理。
进一步地,对叠前地震道集进行预处理是对前期的地震资料作保幅处理,包括道编辑、带通滤波、真振幅恢复、静校正、速度分析、剩余静校正、地表振幅一致性补偿、叠前反褶积及动校正。
进一步地,抽取宏面元道集包括选取宏面元尺度,得到当前CDP点的宏面元道集地震数据的偏移距-方位角分布图,将分布图中当前CDP点周围宏面元尺度倍数个CDP点的地震数据放置到当前的CDP的道集上。
进一步地,寻找最大覆盖的CDP点包括对工区内所有CDP点的CMP道集进行统计,选择CDP点满覆盖的区域,统计出的最大覆盖次数。
进一步地,对最大覆盖次数的单个CDP点的道集进行分方位分偏移距统计包括,针对最大覆盖区域内的CDP点找出最大满偏移距,把单个CMP道集分成多个偏移距和多个方位角个体的分方位分偏移距道集数据体,然后进行分方位角分偏移距统计。
进一步地,对最大覆盖次数的单个CDP点的道集进行宏面元匀化包括,在分方位分偏移距道集数据体中寻找道集数量最少的数据体,把超过该最少道集数量的地震道集切掉。
进一步地,对所有叠前地震道集依次进行分方位分偏移距面元匀化包括,针对每个CDP点重复进行分方位分偏移距统计以及宏面元匀化,得到整个满覆盖区域内的每个方位角和每个偏移距都相等的叠前地震道集数据体。
利用本发明的技术系列,可以部分消除地震方位偏移距道集偏移以及反演过程中的不一致性,消除采集因素的影响,可以为后续的AVAZ反演获取一个较高质量的叠前方位角度(偏移距)道集。
附图说明
通过结合附图对本公开示例性实施方式进行更详细的描述,本公开的上述以及其它目的、特征和优势将变得更加明显,其中,在本公开示例性实施方式中,相同的参考标号通常代表相同部件。
图1显示了某工区某个CDP点的方位角-偏移距分布图。
图2显示了某工区某个CDP点宏面元扩大后的方位角-偏移距分布图。
图3显示了某工区某个CDP点宏面元扩大后的满覆盖区域的方位角-偏移距分布图(由网格划分出来的18个区域即为分方位分偏移距道集的数据体分布情况)。
图4显示了原始方位角-偏移距道集。
图5显示了匀化后的方位角-偏移距道集。
图6显示了根据本发明实施方式的方法流程图。
具体实施方式
下面将参照附图更详细地描述本公开的优选实施方式。虽然附图中显示了本公开的优选实施方式,然而应该理解,可以以各种形式实现本公开而不应被这里阐述的实施方式所限制。相反,提供这些实施方式是为了使本公开更加透彻和完整,并且能够将本公开的范围完整地传达给本领域的技术人员。
本发明涉及地震勘探技术,是一种对叠前道集进行分方位角、偏移距部分匀化,获得规则化方位角度(偏移距)道集的技术。方位角度道集和方位偏移距道集在实际应用中常常是一种道集数据的两种不同的表现形式,两者之间可以互相转换:实际野外地震勘探中采集的第一手地震资料为方位偏移距道集,但是在后期反演应用中常常要用到方位角度道集,方位偏移距道集可以通过地震勘探中的速度转换为方位角度道集。
在叠前AVAZ各向异性分析以及叠前AVAZ各向异性道集反演中,都需要利用高质量的叠前方位角度(偏移距)道集,道集质量直接决定了反演成果的好坏;理想的AVAZ反演道集应该是方位-角度(偏移距)道集在各个方位和角度(偏移距)都是均匀的,但是由于采集观测系统的限制,实际地震勘探中得到的叠前方位角度(偏移距)道集都存在较大的不均一性,因此,本发明提出一种基于三维空间的方位角度道集匀化技术,以此得到一种相对方位角度比较均匀的叠前方位角度(偏移距)道集。
参照图6描述根据本发明实施方式的三维空间的方位偏移距道集匀化的方法,该方法包括如下步骤:
抽取宏面元道集;
对叠前CMP道集进行统计分析,寻找最大覆盖的CDP点;
对最大覆盖次数的单个CDP点的道集进行分方位分偏移距统计;
对最大覆盖次数的单个CDP点的道集进行宏面元匀化;
对所有叠前地震道集依次进行分方位分偏移距面元匀化;
获得分方位分偏移距三维匀化数据的道集。
优选地,该方法还包括在抽取宏面元道集之前对叠前地震道集的预处理。
在本发明中叠前地震资料的预处理是指针对前期的地震资料作保幅处理,包括道编辑、带通滤波、真振幅恢复、静校正、速度分析、剩余静校正、地表振幅一致性补偿、叠前反褶积及动校正。上述预处理的方法都是本领域已知的技术,在此不再赘述。
优选地,抽取宏面元道集包括选取宏面元尺度,得到当前CDP点的宏面元道集地震数据的偏移距-方位角分布图,将分布图中当前CDP点周围宏面元尺度倍数个CDP点的地震数据放置到当前的CDP的道集上。例如,假设当前CDP点的叠前地震道集的道数为Num,选取的宏面元尺度为3×3,则该CDP点的宏面元道集道数为9×Num道,该9×Num道地震数据均对应一个偏移距和一个方位角,因此得到一个9×Num道的地震数据的偏移距-方位角分布图,我们把周围9个CDP点得到的地震数据放置到当前的CDP的道集上,相当于把当前CDP点的道集数量扩大了9倍。优选地,可以根据需要选择不同的宏面元尺度,把当前CDP点的道集数量扩大相应的倍数。
优选地,寻找最大覆盖的CDP点包括对工区内所有CDP点的CMP道集进行统计,选择CDP点满覆盖的区域,统计出的最大覆盖次数。例如,对整个工区的所有单点的CMP道集进行统计,比如对于一个M条线N个道的地震工区而言,一共M×N(如果M=300,N=200,则M×N=60000)个CDP点,每个点对应一个CMP道集。对工区内所有CDP点的CMP道集进行统计,单个CDP点的CMP道集的方位偏移距分布如图1所示。统计完所有的CDP点,选择整个地震数据满覆盖的那个区域(一般这个区域都小于M×N,我们假定这个区域为M-100,N-100),这个区域内的CDP点基本都是满覆盖的,统计出的最大覆盖次数为P(举例为1800),常规的覆盖次数经常在几十次到两三百次不等,但是经过之前的抽取宏面元道集的步骤,相当于把覆盖次数提高了9倍次。
接下来,对最大覆盖次数的单个CDP点的道集进行分方位分偏移距统计包括,针对最大覆盖区域内的CDP点找出最大满偏移距,把单个CMP道集分成多个偏移距和多个方位角个体的分方位分偏移距道集数据体,然后进行分方位角分偏移距统计。例如,针对最大覆盖区域内的CDP点(Max_M,Max_N),找出它最大满偏移距(即所有方位角都有数据的最大偏移距),如图1所示的例子大概为4200米。然后对其进行分方位角分偏移距统计,一般为了叠前AVAZ反演的需要,最少需要分4个方位,在覆盖技术较好的区域,一般可以分为6个方位,偏移距根据需要一般分为3个,当满偏移距为4200米时,一般分0米到1400米为一个偏移距,1400米到2800米为一个偏移距,2800米到4200米为一个偏移距(其他情况可以按照这样进行类推)。依次类推,可以把单个CMP道集分成三个偏移距六个方位角一共18个体的分方位分偏移距道集数据体。
优选地,对最大覆盖次数的单个CDP点的道集进行宏面元匀化包括,在分方位分偏移距道集数据体中寻找道集数量最少的数据体,把超过该最少道集数量的地震道集切掉。例如,在对单个CDP点的叠前地震道进行匀化的步骤中,假设经过宏面元扩大以后的当前最大覆盖次数为1800次,对于每个CDP点而言,18个体的分方位分偏移距道集的数据的道集数量大约为100个。但是由于采集观测系统的限制,每个方位体内部的道集数量必然是不均匀的,因此需要对18个数据体的道集进行匀一化,由于采用了宏面元扩大的办法,当前18个数据体内的道集数量是足够多的,因此我们直接采用最简单的办法,在18个数据体中寻找道集数量最少的那个数据体,比如当前18个数据体中道集数量最少的为60个,数量最多的为110个;直接把超过60个地震道的地震道集直接切掉,就形成了18个数据体道集数量均为60个的地震道集数据体。
最后,对工区内满覆盖区域的叠前地震道集依次进行分方位分偏移距面元匀化。对整个三维地震工区内的满覆盖区域(假设为M-100,N-100),每个CDP点重复进行分方位分偏移距统计以及宏面元匀化,得到整个满覆盖区域内的每个方位角和每个偏移距都相等的叠前地震道集数据体。
为便于理解本发明实施例的方案及其效果,以下给出一个具体应用示例。本领域技术人员应理解,该示例仅为了便于理解本发明,其任何具体细节并非意在以任何方式限制本发明。
如图1所示的某个地震工区经过预处理的某个CDP点的方位角-偏移距道集分布示意图。图2为该CDP点经过9点宏面元扩大后的方位角-偏移距道集分布示意图。图3为满覆盖(最大偏移距为4200米)的方位角-偏移距道集分布示意图,从图3可以看到分成了18个部分,6个不同的方位,3个不同的偏移距。图4为该CDP点原始的道集,图5为该CDP点经过匀化之后的道集,可以看到道集能量有较好的改善。
以上已经描述了本公开的各实施例,上述说明是示例性的,并非穷尽性的,并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下,对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。本文中所用术语的选择,旨在最好地解释各实施例的原理、实际应用或对市场中的技术改进,或者使本技术领域的其它普通技术人员能理解本文披露的各实施例。
Claims (8)
1.一种三维空间的方位偏移距道集匀化的方法,其特征在于,该方法包括如下步骤:
抽取宏面元道集;
对叠前CMP道集进行统计分析,寻找最大覆盖的CDP点;
对最大覆盖次数的单个CDP点的道集进行分方位分偏移距统计;
对最大覆盖次数的单个CDP点的道集进行宏面元匀化;
对所有叠前地震道集依次进行分方位分偏移距面元匀化;
获得分方位分偏移距三维匀化数据的道集。
2.根据权利要求1所述的方法,其特征在于,该方法还包括在抽取宏面元道集之前对叠前地震道集进行预处理。
3.根据权利要求2所述的方法,其特征在于,对叠前地震道集进行预处理是对前期的地震资料作保幅处理,包括道编辑、带通滤波、真振幅恢复、静校正、速度分析、剩余静校正、地表振幅一致性补偿、叠前反褶积及动校正。
4.根据权利要求1所述的方法,其特征在于,抽取宏面元道集包括选取宏面元尺度,得到当前CDP点的宏面元道集地震数据的偏移距-方位角分布图,将分布图中当前CDP点周围宏面元尺度倍数个CDP点的地震数据放置到当前的CDP的道集上。
5.根据权利要求1所述的方法,其特征在于,寻找最大覆盖的CDP点包括对工区内所有CDP点的CMP道集进行统计,选择CDP点满覆盖的区域,统计出的最大覆盖次数。
6.根据权利要求1所述的方法,其特征在于,对最大覆盖次数的单个CDP点的道集进行分方位分偏移距统计包括,针对最大覆盖区域内的CDP点找出最大满偏移距,把单个CMP道集分成多个偏移距和多个方位角个体的分方位分偏移距道集数据体,然后进行分方位角分偏移距统计。
7.根据权利要求6所述的方法,其特征在于,对最大覆盖次数的单个CDP点的道集进行宏面元匀化包括,在分方位分偏移距道集数据体中寻找道集数量最少的数据体,把超过该最少道集数量的地震道集切掉。
8.根据权利要求1所述的方法,其特征在于,对所有叠前地震道集依次进行分方位分偏移距面元匀化包括,针对每个CDP点重复进行分方位分偏移距统计以及宏面元匀化,得到整个满覆盖区域内的每个方位角和每个偏移距都相等的叠前地震道集数据体。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610236671.7A CN107300717B (zh) | 2016-04-15 | 2016-04-15 | 一种三维空间的方位偏移距道集匀化的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610236671.7A CN107300717B (zh) | 2016-04-15 | 2016-04-15 | 一种三维空间的方位偏移距道集匀化的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107300717A CN107300717A (zh) | 2017-10-27 |
CN107300717B true CN107300717B (zh) | 2019-06-04 |
Family
ID=60137307
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610236671.7A Active CN107300717B (zh) | 2016-04-15 | 2016-04-15 | 一种三维空间的方位偏移距道集匀化的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107300717B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109471163B (zh) * | 2018-10-29 | 2020-07-21 | 中国海洋石油集团有限公司 | 基于地质体连续性建模的高精度井拓反演方法 |
CN110068863A (zh) * | 2019-04-04 | 2019-07-30 | 中国科学院计算机网络信息中心 | 一种共深度点cdp的确定方法及装置 |
CN112305599A (zh) * | 2019-07-30 | 2021-02-02 | 中国石油天然气集团有限公司 | 高密度宽方位的地震数据分类方法及装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102033242A (zh) * | 2010-10-22 | 2011-04-27 | 中国石油化工股份有限公司 | 一种深层倾斜裂缝储层地震振幅预测方法 |
CN103076623A (zh) * | 2011-10-25 | 2013-05-01 | 中国石油化工股份有限公司 | 一种基于叠前相干的裂缝检测方法 |
EP2728382A2 (en) * | 2012-11-06 | 2014-05-07 | CGG Services SA | Spatial expansion seismic data processing method and apparatus. |
CN105445787A (zh) * | 2014-08-11 | 2016-03-30 | 中国石油化工股份有限公司 | 一种最优方位子体相干的裂缝预测方法 |
-
2016
- 2016-04-15 CN CN201610236671.7A patent/CN107300717B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102033242A (zh) * | 2010-10-22 | 2011-04-27 | 中国石油化工股份有限公司 | 一种深层倾斜裂缝储层地震振幅预测方法 |
CN103076623A (zh) * | 2011-10-25 | 2013-05-01 | 中国石油化工股份有限公司 | 一种基于叠前相干的裂缝检测方法 |
EP2728382A2 (en) * | 2012-11-06 | 2014-05-07 | CGG Services SA | Spatial expansion seismic data processing method and apparatus. |
CN105445787A (zh) * | 2014-08-11 | 2016-03-30 | 中国石油化工股份有限公司 | 一种最优方位子体相干的裂缝预测方法 |
Non-Patent Citations (2)
Title |
---|
塔河油田南部缝孔型碳酸盐岩储层叠前地震预测;管路平;《石油与天然气地质》;20081031;第29卷(第5期);第654-667、675页 |
碳酸盐岩裂缝型储层叠前地震预测技术应用;邓道静;《勘探地球物理进展》;20090430;第32卷(第2期);第133-137页 |
Also Published As
Publication number | Publication date |
---|---|
CN107300717A (zh) | 2017-10-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101031821B (zh) | 用于迁移和速度分析的偏移-方位建仓 | |
CN107300717B (zh) | 一种三维空间的方位偏移距道集匀化的方法 | |
Jones et al. | P-wave tomography of potential convective downwellings and their source regions, Sierra Nevada, California | |
Allam et al. | Ten kilometer vertical Moho offset and shallow velocity contrast along the Denali fault zone from double-difference tomography, receiver functions, and fault zone head waves | |
CN104280775B (zh) | 一种基于全波形矢量偏移叠加的微地震监测定位方法 | |
CN104483705B (zh) | 一种三维剩余静校正方法 | |
CN103217710B (zh) | 获取共偏移距共方位角域道集的方法及系统 | |
CN103744115B (zh) | 在共偏移距共方位角域进行数据规则化处理的方法及系统 | |
CN104199103B (zh) | 转换波静校正方法及装置 | |
Yang et al. | Stress distribution near the seismic gap between Wenchuan and Lushan earthquakes | |
CN105806208A (zh) | 一种基于gnss网形变化的形变异常检测方法 | |
CN103149592A (zh) | 一种变偏移距vsp波场分离方法 | |
Wei et al. | Imaging the Topography of Crust–Mantle Boundary from a High‐Density Seismic Array beneath the Middle‐Lower Yangtze River, Eastern China | |
Petit et al. | A maximum entropy approach to detect close-in giant planets around active stars | |
CN104914469B (zh) | 一种转换横波中的静校正方法及装置 | |
Bulut et al. | Spatiotemporal earthquake clusters along the North Anatolian fault zone offshore Istanbul | |
CN107918157A (zh) | 三分量p波初动震源机制反演方法及装置 | |
CN106908836A (zh) | 采集脚印压制方法及系统 | |
Douglas et al. | Compilation and critical review of GMPEs for the GEM-PEER global GMPEs project | |
Haggag et al. | Seismicity and 3D velocity structure in the Aswan Reservoir Lake area, Egypt | |
Okada et al. | Source processes of 15 September 1998 M 5.0 Sendai, northeastern Japan, earthquake and its M 3.8 foreshock by waveform inversion | |
Strobbia et al. | Double-array refraction microtremors | |
CN104133239A (zh) | 一种地震单炮记录目的层时窗的定义方法 | |
Salathe et al. | Energy reconstruction of an n-type segmented inverted coaxial point-contact HPGe detector | |
Hejrani et al. | Teleseismic tomography when stations follow profiles: pitfalls and remedies |
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 |