CN109884698B - 基于目的层的地震勘探观测系统定量评价方法 - Google Patents

基于目的层的地震勘探观测系统定量评价方法 Download PDF

Info

Publication number
CN109884698B
CN109884698B CN201910216727.6A CN201910216727A CN109884698B CN 109884698 B CN109884698 B CN 109884698B CN 201910216727 A CN201910216727 A CN 201910216727A CN 109884698 B CN109884698 B CN 109884698B
Authority
CN
China
Prior art keywords
observation system
illumination energy
target layer
shot
designed
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
Application number
CN201910216727.6A
Other languages
English (en)
Other versions
CN109884698A (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
Geophysical Research Institute of Sinopec Shengli Oilfield Co
Original Assignee
China Petroleum and Chemical Corp
Geophysical Research Institute of Sinopec Shengli Oilfield Co
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, Geophysical Research Institute of Sinopec Shengli Oilfield Co filed Critical China Petroleum and Chemical Corp
Priority to CN201910216727.6A priority Critical patent/CN109884698B/zh
Publication of CN109884698A publication Critical patent/CN109884698A/zh
Application granted granted Critical
Publication of CN109884698B publication Critical patent/CN109884698B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明提供一种基于目的层的地震勘探观测系统定量评价方法，包括：根据工区地质任务和技术要求设计观测系统参数；根据设计的观测系统参数和满覆盖次数边界建立炮点集、检波点集及炮检关系集；根据工区地震资料解释的层位数据建立所有目的层的模型；选择一个目的层,设置照明计算参数；根据设计的观测系统参数计算所选目的层每个面元的照明能量；统计目的层每个面元最大照明能量比；计算目的层面元照明能量均方差；计算设计的观测系统综合质量因子；根据需要选择其他目的层进行计算；根据计算的所有目的层观测系统综合质量因子，对观测系统做出评价。该方法对照明度能量进行细分，按照偏移距、方位角等多种进行综合评价。

Description

基于目的层的地震勘探观测系统定量评价方法
技术领域
本发明涉及油气地震勘探数据采集设计领域，特别是涉及到一种基于目的层的地震勘探观测系统定量评价方法。
背景技术
观测系统评价工作主要分为两大类：第一类是基于共中心点面元属性分析方法来开展的。这种方法在地震资料采集施工中占有重要的地位。其主要包含覆盖次数分析、方位角均匀性分析、偏移距均匀性分析以及玫瑰图均匀性分析等手段，在很多地区的观测系统设计中证明是切实有效的。这类方法在地表极其复杂、地下构造极其复杂的地区暴露出较大的问题。第二类是基于共反射点面元的分析方法。这类方法也包含了覆盖次数、方位角、偏移距等分析手段，但与第一类方法不同的是，这类方法是基于真实反射点进行的分析，大量的理论和实际资料证明，在这些复杂地区，进行面向目标地质体的共反射点面元分析才是观测系统评价有效途径。为了刚直观地进行分析，通常采用目标地质体的照明度来进行综合分析。目前生产中广泛认可第二类分析方法，这类方法又被称为面向目标地质体的观测系统评价方法。一般包含了模型建立、正演照明分析、观测系统评价三个步骤。
照明度是一个能量的概念，作为面向目标的观测系统设计的基础数据，照明度的计算尤为重要。照明度可以分为很多种不同的统计方式，最核心的数据为炮点-面元-检波点形成的炮检能量对，这些基础数据对观测系统评价具有重要的意义。目前通过照明度进行观测系统评价主要是计算照明度的均匀性，均匀性好的则认为观测系统好，反之，观测系统质量较差。这种分析方式虽然可以直观地获得不同观测系统的评价值，但是也存在较大的问题。照明能量是某一面元内不同方位角、不同偏移距能量的叠加，只评价面元整体能量值的分布，观测系统中包含的方位角、偏移距、覆盖次数信息，地震记录中包含的频率、能量信息等无法综合的表现出来。因此对照明度能量进行细分，按照偏移距、方位角等多种进行综合评价，可以更有助于对观测系统做出评价。为此我们发明了一种新的基于目的层的地震勘探观测系统定量评价方法，解决了以上技术问题。
发明内容
本发明的目的是提供一种对照明度能量进行细分，按照偏移距、方位角等多种进行综合评价的基于目的层的地震勘探观测系统定量评价方法。
本发明的目的可通过如下技术措施来实现：基于目的层的地震勘探观测系统定量评价方法，该基于目的层的地震勘探观测系统定量评价方法包括：步骤1：根据工区地质任务和技术要求设计观测系统参数；步骤2：根据设计的观测系统参数和满覆盖次数边界建立炮点集、检波点集及炮检关系集；步骤3：根据工区地震资料解释的层位数据建立所有目的层的模型；步骤4：选择一个目的层,设置照明计算参数；步骤5:根据设计的观测系统参数计算所选目的层每个面元的照明能量；步骤6:根据步骤5计算结果，统计目的层每个面元最大照明能量比；步骤7：根据步骤5计算结果，计算目的层面元照明能量均方差；步骤8：根据步骤5、6、7计算结果，计算设计的观测系统综合质量因子；步骤9：返回步骤4，根据需要选择其他目的层进行计算；步骤10：根据计算的所有目的层观测系统综合质量因子，对观测系统做出评价。
本发明的目的还可通过如下技术措施来实现：
在步骤5中，根据设计的观测系统参数计算所选目的层每个面元的照明能量，所使用的正演方法包括射线追踪正演、波动方程正演和高斯射线束正演方法。
在步骤5中，根据设计的观测系统参数计算所选目的层每个面元的照明能量，计算面元内每个炮检对的照明能量的同时还要记录面元内每个炮检对的偏移距方位角信息。
在步骤5中，根据设计的观测系统参数计算所选目的层每个面元的照明能量，对于不同的震源，计算与该震源关联的炮检对照明能量时需要乘上一个震源系数，对于井炮激发，该震源系数定义为实际药量与最大药量的比值；对于可控震源，该系数根据可控震源与实际井炮最大药量激发的单炮对比关系设置。
在步骤6中，根据步骤5计算结果，统计目的层每个面元最大照明能量比，首先对单个反射面元的炮检对照明能量按照方位角和偏移距进行统计，具体角度间隔和偏移距间隔根据实际资料情况确定，方位角按照10 度划分为36个区间，偏移距根据最大偏移距等分为3个区间，i表示偏移距区间号，j表示方位角区间号，Ei,j表示位于第i个偏移距区间同时位于第j个方位角区间的所有炮检对能量之和；最后，计算单个面元内炮检对能量总和与Ei,j中的最大值之比作为该面元最大照明能量比。
在步骤7中，根据步骤5计算结果，计算目的层面元照明能量均方差，是将一个面元内的炮检能量对进行均方差计算，从而获得每个面元内照明能量的均匀性分布情况。
在步骤8中，根据步骤5、6、7计算结果，计算设计的观测系统综合质量因子，是按照以下公式进行计算的：
Figure BDA0002001142220000031
其中，Q为观测系统综合质量因子，Iave为所有面元照明能量的平均值， Emax为所有面元照明能量均方差的最大值，Evar为所有面元照明能量均方差的最大值与最小值之差，Rvar为所有面元最大照明能量比的最大值与最小值之差。
在步骤9中，对工区所有目的层进行计算，综合各目的层的计算结果对设计的观测系统进行评价。
在步骤10中，根据计算的所有目的层观测系统综合质量因子，对观测系统进行评价时，综合质量因子越大，观测系统质量越好。
在步骤10中，对观测系统进行评价时是在步骤1时设计若干个备选的观测系统方案，然后对备选的观测系统进行评价，从中选择综合质量因子最高的观测系统。
本发明中的基于目的层的地震勘探观测系统定量评价方法，可以快速的量化不同观测系统的质量。相对于其他的观测系统评价方法，该方法有不可比拟的优势，其具体优势和特点表现在：(1)原理可靠。常用的覆盖次数分布均匀性分析、方位角和偏移距分布玫瑰图分析、照明能量均方差分析都在该发明中有所体现，通过照明综合质量因子将他们结合在一起，较单一评价方法更为可靠。(2)分析简便。通过输入面元的照明能量，自动统计能量分布特征，最终每个观测系统只有一个数值，便于分析对比。
附图说明
图1为本发明的基于目的层的地震勘探观测系统定量评价方法的一具体实施例的流程图；
图2是实施例1中同观测系统参数图；
图3是实施例1中观测系统1常规覆盖次数分布图；
图4是实施例1中观测系统2常规覆盖次数分布图；
图5是实施例1中观测系统1针对某一目的层的照明结果图；
图6是实施例1中观测系统2针对某一目的层的照明结果图；
图7是实施例2中同观测系统参数图；
图8是实施例2中观测系统1常规覆盖次数分布图；
图9是实施例2中观测系统2常规覆盖次数分布图；
图10是实施例2中观测系统1针对某一目的层的照明结果图；
图11是实施例2中观测系统2针对某一目的层的照明结果图。
具体实施方式
为使本发明的上述和其他目的、特征和优点能更明显易懂，下文特举出较佳实施例，并配合附图所示，作详细说明如下。
如图1所示，图1为本发明的基于目的层的地震勘探观测系统定量评价方法的流程图。
步骤101：根据工区地质任务和技术要求设计观测系统参数。
步骤102：根据设计的观测系统参数和满覆盖次数边界建立炮点集、检波点集及炮检关系集。
步骤103：根据工区地震资料解释的层位数据建立所有目的层的模型。
步骤104：选择一个目的层,设置照明计算参数。
步骤105:根据设计的观测系统参数计算所选目的层每个面元的照明能量。根据设计的观测系统参数计算所选目的层每个面元的照明能量，所使用的正演方法包括射线追踪正演、波动方程正演和高斯射线束正演方法。
根据设计的观测系统参数计算所选目的层每个面元的照明能量，计算面元内每个炮检对的照明能量的同时还要记录面元内每个炮检对的偏移距方位角信息。
根据设计的观测系统参数计算所选目的层每个面元的照明能量，对于不同的震源，计算与该震源关联的炮检对照明能量时需要乘上一个震源系数，对于井炮激发，该震源系数定义为实际药量与最大药量的比值；对于可控震源，该系数根据可控震源与实际井炮最大药量激发的单炮对比关系设置。
平滑处理算法采用正方形邻域移动平均算法以节省搜索时间。
步骤106:根据步骤105计算结果，统计目的层每个面元最大照明能量比。根据步骤105计算结果，统计目的层每个面元最大照明能量比，首先对单个反射面元的炮检对照明能量按照方位角和偏移距进行统计，具体角度间隔和偏移距间隔根据实际资料情况确定，通常方位角按照10度划分为36个区间，偏移距根据最大偏移距等分为3个区间，假设i表示偏移距区间号，j表示方位角区间号，Ei,j表示位于第i个偏移距区间同时位于第 j个方位角区间的所有炮检对能量之和。最后，计算单个面元内炮检对能量总和与Ei,j中的最大值之比作为该面元最大照明能量比。
步骤107：根据步骤105计算结果，计算目的层面元照明能量均方差。根据步骤105计算结果，计算目的层面元照明能量均方差，是将一个面元内的炮检能量对进行均方差计算，从而获得每个面元内照明能量的均匀性分布情况。
步骤108：根据步骤105、106、107计算结果，计算设计的观测系统综合质量因子。根据步骤105、106、107计算结果，计算设计的观测系统综合质量因子，是按照以下公式进行计算的：
Figure BDA0002001142220000061
其中，Q为观测系统综合质量因子，Iave为所有面元照明能量的平均值， Emax为所有面元照明能量均方差的最大值，Evar为所有面元照明能量均方差的最大值与最小值之差，Rvar为所有面元最大照明能量比的最大值与最小值之差。
步骤109：返回步骤104，根据需要选择其他目的层进行计算。一般需要对工区所有目的层进行计算，综合各目的层的计算结果对设计的观测系统进行评价。
步骤110：根据计算的所有目的层观测系统综合质量因子，对观测系统做出评价。
根据计算的所有目的层观测系统综合质量因子，对观测系统进行评价时，综合质量因子越大，观测系统质量越好。
对观测系统进行评价时一般是在步骤110时设计若干个备选的观测系统方案，然后对备选的观测系统进行评价，从中选择综合质量因子最高的观测系统。
以下为应用本发明的几个具体实施例：
实施例1，来源于中石化胜利油田一个三维地震采集观测系统设计项目，该工区位于山东省东营市内，具体实施方式为：
(1)根据工区地质任务和技术指标设计了两套观测系统，观测系统参数如图2；
(2)根据设计的观测系统参数和满覆盖次数边界建立炮点集、检波点集及炮检关系集；
(3)计算常规观测系统属性，很难直观地选择一个相对更合适的观测系统(如图3和图4)；
(3)根据工区地震资料解释的层位数据建立所有目的层的模型；
(4)选择最深的目的层,设置照明计算参数。
(5)根据设计的观测系统参数计算所选目的层每个面元的照明能量，从目的层照明来看，两套观测系统具有明显差异(如图5和图6)，但仍不能量化评价哪个更好。
(6)根据步骤5计算结果，统计目的层每个面元最大照明能量比。
(7)根据步骤5计算结果，计算目的层面元照明能量均方差。
(8)根据步骤5、6、7计算结果，计算设计的观测系统综合质量因子，方案1为6.7，方案2为5.1，因此方案1质量更好。
实施例2，来源于中石化胜利油田新疆探区一个三维地震采集观测系统设计项目，该工区位于新疆维吾尔自治区内，具体实施方式为：
(1)根据工区地质任务和技术指标设计了两套观测系统，观测系统参数如图7；
(2)根据设计的观测系统参数和满覆盖次数边界建立炮点集、检波点集及炮检关系集；
(3)计算常规观测系统属性，很难直观地选择一个相对更合适的观测系统(如图8和图9)；
(3)根据工区地震资料解释的层位数据建立所有目的层的模型；
(4)选择最深的目的层,设置照明计算参数。
(5)根据设计的观测系统参数计算所选目的层每个面元的照明能量，从目的层照明来看，两套观测系统具有明显差异(如图10和图11)，但仍不能量化评价哪个更好。
(6)根据步骤5计算结果，统计目的层每个面元最大照明能量比。
(7)根据步骤5计算结果，计算目的层面元照明能量均方差。
(8)根据步骤5、6、7计算结果，计算设计的观测系统综合质量因子，方案1为7.5，方案2为4.8，因此方案1质量更好。

Claims (9)

1.基于目的层的地震勘探观测系统定量评价方法，其特征在于，该基于目的层的地震勘探观测系统定量评价方法包括：
步骤1：根据工区地质任务和技术要求设计观测系统参数；
步骤2：根据设计的观测系统参数和满覆盖次数边界建立炮点集、检波点集及炮检关系集；
步骤3：根据工区地震资料解释的层位数据建立所有目的层的模型；
步骤4：选择一个目的层,设置照明计算参数；
步骤5:根据设计的观测系统参数计算所选目的层每个面元的照明能量；
步骤6:根据步骤5计算结果，统计目的层每个面元最大照明能量比；
步骤7：根据步骤5计算结果，计算目的层面元照明能量均方差；
步骤8：根据步骤5、6、7计算结果，计算设计的观测系统综合质量因子；
步骤9：返回步骤4，根据需要选择其他目的层进行计算；
步骤10：根据计算的所有目的层观测系统综合质量因子，对观测系统做出评价；
在步骤6中，根据步骤5计算结果，统计目的层每个面元最大照明能量比，首先对单个反射面元的炮检对照明能量按照方位角和偏移距进行统计，具体角度间隔和偏移距间隔根据实际资料情况确定，方位角按照10度划分为36个区间，偏移距根据最大偏移距等分为3个区间，i表示偏移距区间号，j表示方位角区间号，Ei,j表示位于第i个偏移距区间同时位于第j个方位角区间的所有炮检对能量之和；最后，计算单个面元内炮检对能量总和与Ei,j中的最大值之比作为该面元最大照明能量比。
2.根据权利要求1所述的基于目的层的地震勘探观测系统定量评价方法，其特征在于，在步骤5中，根据设计的观测系统参数计算所选目的层每个面元的照明能量，所使用的正演方法包括射线追踪正演、波动方程正演和高斯射线束正演方法。
3.根据权利要求1所述的基于目的层的地震勘探观测系统定量评价方法，其特征在于，在步骤5中，根据设计的观测系统参数计算所选目的层每个面元的照明能量，计算面元内每个炮检对的照明能量的同时还要记录面元内每个炮检对的偏移距方位角信息。
4.根据权利要求1所述的基于目的层的地震勘探观测系统定量评价方法，其特征在于，在步骤5中，根据设计的观测系统参数计算所选目的层每个面元的照明能量，对于不同的震源，计算与该震源关联的炮检对照明能量时需要乘上一个震源系数，对于井炮激发，该震源系数定义为实际药量与最大药量的比值；对于可控震源，该系数根据可控震源与实际井炮最大药量激发的单炮对比关系设置。
5.根据权利要求1所述的基于目的层的地震勘探观测系统定量评价方法，其特征在于，在步骤7中，根据步骤5计算结果，计算目的层面元照明能量均方差，是将一个面元内的炮检能量对进行均方差计算，从而获得每个面元内照明能量的均匀性分布情况。
6.根据权利要求1所述的基于目的层的地震勘探观测系统定量评价方法，其特征在于，在步骤8中，根据步骤5、6、7计算结果，计算设计的观测系统综合质量因子，是按照以下公式进行计算的：
Figure FDA0002823777680000021
其中，Q为观测系统综合质量因子，Iave为所有面元照明能量的平均值，Emax为所有面元照明能量均方差的最大值，Evar为所有面元照明能量均方差的最大值与最小值之差，Rvar为所有面元最大照明能量比的最大值与最小值之差。
7.根据权利要求1所述的基于目的层的地震勘探观测系统定量评价方法，其特征在于，在步骤9中，对工区所有目的层进行计算，综合各目的层的计算结果对设计的观测系统进行评价。
8.根据权利要求1所述的基于目的层的地震勘探观测系统定量评价方法，其特征在于，在步骤10中，根据计算的所有目的层观测系统综合质量因子，对观测系统进行评价时，综合质量因子越大，观测系统质量越好。
9.根据权利要求1所述的基于目的层的地震勘探观测系统定量评价方法，其特征在于，在步骤10中，对观测系统进行评价时是在步骤1时设计若干个备选的观测系统方案，然后对备选的观测系统进行评价，从中选择综合质量因子最高的观测系统。
CN201910216727.6A 2019-03-20 2019-03-20 基于目的层的地震勘探观测系统定量评价方法 Active CN109884698B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910216727.6A CN109884698B (zh) 2019-03-20 2019-03-20 基于目的层的地震勘探观测系统定量评价方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910216727.6A CN109884698B (zh) 2019-03-20 2019-03-20 基于目的层的地震勘探观测系统定量评价方法

Publications (2)

Publication Number Publication Date
CN109884698A CN109884698A (zh) 2019-06-14
CN109884698B true CN109884698B (zh) 2021-04-27

Family

ID=66933431

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910216727.6A Active CN109884698B (zh) 2019-03-20 2019-03-20 基于目的层的地震勘探观测系统定量评价方法

Country Status (1)

Country Link
CN (1) CN109884698B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113156494B (zh) * 2020-01-07 2023-10-31 中国石油天然气集团有限公司 面元等效覆盖次数的计算方法及装置
CN113126158B (zh) * 2021-04-20 2022-06-07 北京众合华创生物环保科技有限公司 一种面向全波场的超高密度三维地震采集方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103543465A (zh) * 2012-07-12 2014-01-29 中国石油天然气集团公司 基于目的层照明能量确定有效激发点的方法
CN104345336A (zh) * 2013-07-25 2015-02-11 中国石油天然气集团公司 一种基于目标区域照明度的观测系统优化方法
CN104502962A (zh) * 2014-12-16 2015-04-08 中国石油天然气集团公司 一种设计炮点的方法及装置
CN105824042A (zh) * 2015-01-08 2016-08-03 中石化石油工程地球物理有限公司胜利分公司 基于照明能量最优的最大纵向距设计方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10551517B2 (en) * 2014-09-04 2020-02-04 Westerngeco L.L.C. Multiples mitigation and imaging with incomplete seismic data

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103543465A (zh) * 2012-07-12 2014-01-29 中国石油天然气集团公司 基于目的层照明能量确定有效激发点的方法
CN104345336A (zh) * 2013-07-25 2015-02-11 中国石油天然气集团公司 一种基于目标区域照明度的观测系统优化方法
CN104502962A (zh) * 2014-12-16 2015-04-08 中国石油天然气集团公司 一种设计炮点的方法及装置
CN105824042A (zh) * 2015-01-08 2016-08-03 中石化石油工程地球物理有限公司胜利分公司 基于照明能量最优的最大纵向距设计方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Geometry Matching Technology of Non-Repeating Acquired Time-Lapse Seismic Data Processing;CUI Qinghui;《Advances in Petroleum Exploration and Development》;20160331;第11卷(第1期);第55-60页 *
提高目的层阴影区成像质量的观测系统优化设计方法;许银坡等;《石油地球物理勘探》;20151215;第50卷(第06期);第1048-1053页 *

Also Published As

Publication number Publication date
CN109884698A (zh) 2019-06-14

Similar Documents

Publication Publication Date Title
CN100514087C (zh) 地震采集三维观测系统定量分析方法
US6791900B2 (en) Method of calculating a throw volume for quantitative fault analysis
US9915742B2 (en) Method and system for geophysical modeling of subsurface volumes based on label propagation
US10073190B2 (en) Method and system for geophysical modeling of subsurface volumes based on computed vectors
CN105044799B (zh) 确定三维地震观测系统面元属性均匀程度及均匀化的方法
CN109884698B (zh) 基于目的层的地震勘探观测系统定量评价方法
WO2014099201A1 (en) Geophysical modeling of subsurface volumes based on horizon extraction
CN104834006B (zh) 结构成因沉积参数图的处理方法和装置
CN109425900A (zh) 一种地震储层预测方法
CN103777248A (zh) 一种适用于不规则发射回线的tem一维正演方法
CN109855593A (zh) 一种钛合金车削加工三维表面形貌分析方法
US20220221604A1 (en) Method of analysing seismic data
CN104345336A (zh) 一种基于目标区域照明度的观测系统优化方法
CN109163674B (zh) 一种面结构光自动化三维测量中传感器测量视点规划方法
CN112946751B (zh) 一种预测不同储层平面分布的方法及系统
CN110469327B (zh) 一种五步逼近式碳酸盐岩储层测井立体评价方法
CN112114358B (zh) 一种基于三维地震资料表征的地下火山通道识别方法
CN104570062B (zh) 一种以激发为中心的vsp观测系统设计方法
CN103777244A (zh) 一种地震裂缝属性体的定量分析方法
CN116027454A (zh) 一种砂岩厚度等值线图的成图方法及智能绘制装置
CN107015289B (zh) 圈闭地质资源量确定方法及装置
CN107942373A (zh) 基于裂缝性油气储层断裂系统检测的相干算法
Qi* et al. Segmentation of salt domes, mass transport complexes on 3D seismic data volumes using Kuwahara windows and multiattribute cluster analysis
CN106896402B (zh) 基于地质单元体的地震正演方法及装置
CN111965705A (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