CN107300717A - 一种三维空间的方位偏移距道集匀化的方法 - Google Patents

一种三维空间的方位偏移距道集匀化的方法 Download PDF

Info

Publication number
CN107300717A
CN107300717A CN201610236671.7A CN201610236671A CN107300717A CN 107300717 A CN107300717 A CN 107300717A CN 201610236671 A CN201610236671 A CN 201610236671A CN 107300717 A CN107300717 A CN 107300717A
Authority
CN
China
Prior art keywords
trace gather
point
gather
offset distance
trace
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
Application number
CN201610236671.7A
Other languages
English (en)
Other versions
CN107300717B (zh
Inventor
梁志强
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Petroleum and Chemical Corp
Sinopec Geophysical Research Institute
Original Assignee
China Petroleum and Chemical Corp
Sinopec Geophysical Research Institute
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China Petroleum and Chemical Corp, Sinopec Geophysical Research Institute filed Critical China Petroleum and Chemical Corp
Priority to CN201610236671.7A priority Critical patent/CN107300717B/zh
Publication of CN107300717A publication Critical patent/CN107300717A/zh
Application granted granted Critical
Publication of CN107300717B publication Critical patent/CN107300717B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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点重复进行分方位分偏移距统计以及宏面元匀化，得到整个满覆盖区域内的每个方位角和每个偏移距都相等的叠前地震道集数据体。
CN201610236671.7A 2016-04-15 2016-04-15 一种三维空间的方位偏移距道集匀化的方法 Active CN107300717B (zh)

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 true CN107300717A (zh) 2017-10-27
CN107300717B 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)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109471163A (zh) * 2018-10-29 2019-03-15 中国海洋石油集团有限公司 基于地质体连续性建模的高精度井拓反演方法
CN110068863A (zh) * 2019-04-04 2019-07-30 中国科学院计算机网络信息中心 一种共深度点cdp的确定方法及装置
CN112305599A (zh) * 2019-07-30 2021-02-02 中国石油天然气集团有限公司 高密度宽方位的地震数据分类方法及装置

Citations (4)

* Cited by examiner, † Cited by third party
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 中国石油化工股份有限公司 一种最优方位子体相干的裂缝预测方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
Title
管路平: "塔河油田南部缝孔型碳酸盐岩储层叠前地震预测", 《石油与天然气地质》 *
邓道静: "碳酸盐岩裂缝型储层叠前地震预测技术应用", 《勘探地球物理进展》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109471163A (zh) * 2018-10-29 2019-03-15 中国海洋石油集团有限公司 基于地质体连续性建模的高精度井拓反演方法
CN110068863A (zh) * 2019-04-04 2019-07-30 中国科学院计算机网络信息中心 一种共深度点cdp的确定方法及装置
CN112305599A (zh) * 2019-07-30 2021-02-02 中国石油天然气集团有限公司 高密度宽方位的地震数据分类方法及装置

Also Published As

Publication number Publication date
CN107300717B (zh) 2019-06-04

Similar Documents

Publication Publication Date Title
CN101031821B (zh) 用于迁移和速度分析的偏移-方位建仓
CN102692645B (zh) 利用纵波、转换波数据联合反演储层纵横波速度比的方法
CN102053261B (zh) 一种地震数据处理方法
CN104199103B (zh) 转换波静校正方法及装置
CN104536034A (zh) 多震源并行激发采集与混合地震记录分离方法
CN104656141A (zh) 共姿态道集地表一致性振幅补偿方法及装置
CN104155691A (zh) 转换波各向异性速度分析方法及装置
CN107300717A (zh) 一种三维空间的方位偏移距道集匀化的方法
CN109765615A (zh) 一种地层品质因子反演方法及装置
CN105866833A (zh) Vsp-cdp叠加方法及三维vsp覆盖次数计算方法
CN105093319A (zh) 基于三维地震数据的地面微地震静校正方法
CN103217710A (zh) 获取共偏移距共方位角域道集的方法及系统
Scarponi et al. Joint seismic and gravity data inversion to image intra-crustal structures: the Ivrea Geophysical Body along the Val Sesia profile (Piedmont, Italy)
CN104267432B (zh) 一种基于规则化的转换波共转换点道集高精度抽取方法
CN102692651A (zh) 速度空变的初至波剩余静校正方法
CN106199704A (zh) 一种三维三分量海底电缆地震资料速度建模方法
CN104422955B (zh) 一种利用旅行时变化量进行各向异性参数提取的方法
CN106842302A (zh) 一种批量编辑初至的方法及装置
CN105866839A (zh) 基于共姿态道集的静校正方法及装置
CN106054252A (zh) 一种叠前时间偏移的方法及装置
CN106257309A (zh) 叠后地震数据体处理方法及装置
CN103472482B (zh) 基于基因排序体系的多波地震资料时间域匹配方法及系统
CN106908836A (zh) 采集脚印压制方法及系统
CN106199705A (zh) 地震勘探采集数据的检测方法
CN111077577B (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