CN112505760B - 一种基于虚拟井的断溶体储层刻画方法 - Google Patents

一种基于虚拟井的断溶体储层刻画方法 Download PDF

Info

Publication number
CN112505760B
CN112505760B CN202011102576.0A CN202011102576A CN112505760B CN 112505760 B CN112505760 B CN 112505760B CN 202011102576 A CN202011102576 A CN 202011102576A CN 112505760 B CN112505760 B CN 112505760B
Authority
CN
China
Prior art keywords
reservoir
fracture
well
virtual
broken
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
CN202011102576.0A
Other languages
English (en)
Other versions
CN112505760A (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.)
Petrochina Co Ltd
Original Assignee
Petrochina Co Ltd
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 Petrochina Co Ltd filed Critical Petrochina Co Ltd
Priority to CN202011102576.0A priority Critical patent/CN112505760B/zh
Publication of CN112505760A publication Critical patent/CN112505760A/zh
Application granted granted Critical
Publication of CN112505760B publication Critical patent/CN112505760B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30Analysis
    • G01V1/306Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/282Application of seismic models, synthetic seismograms
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30Analysis
    • G01V1/301Analysis for determining seismic cross-sections or geostructures
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/40Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/61Analysis by combining or comparing a seismic data set with other data
    • G01V2210/616Data from specific type of measurement
    • G01V2210/6161Seismic or acoustic, e.g. land or sea measurements
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/61Analysis by combining or comparing a seismic data set with other data
    • G01V2210/616Data from specific type of measurement
    • G01V2210/6169Data from specific type of measurement using well-logging
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/62Physical property of subsurface
    • G01V2210/622Velocity, density or impedance
    • G01V2210/6226Impedance
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/62Physical property of subsurface
    • G01V2210/624Reservoir parameters
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/64Geostructures, e.g. in 3D data cubes
    • G01V2210/642Faults

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

本发明公开了一种基于虚拟井的断溶体储层刻画方法。该断溶体储层刻画方法包括以下步骤:S100、选取基准井,在断裂两侧的非储层位置建立虚拟井;S200、将所述基准井的声波曲线赋予所述虚拟井;利用所述虚拟井进行地震反演,以此建立反映地层格架的波阻抗模型;S300、以目标区内实钻井参与进行地震反演;S400、将S300得到的地震反演结果与所述反映地层格架的波阻抗模型做差,刻画所述断溶体储层。本发明可有效刻画断裂破碎带储层,解决了断溶体储层内部反射弱、地层界面反射影响大、储层预测难度大等问题,对指导深层断控储层的勘探开发具有重要意义。

Description

一种基于虚拟井的断溶体储层刻画方法
技术领域
本发明涉及油气勘探技术领域,具体涉及一种基于虚拟井的断溶体储层刻画方法。
背景技术
碳酸盐岩断溶体储层主要受控于走滑断裂,断裂破碎带为主要储集空间,在地震上主要表现为杂乱反射特征,其反射强度较弱,同时受地层界面影响,常规地震属性和地震反演难以对断溶体内幕储层进行有效刻画,现有的结构张量等地震结构类属性可刻画断溶体边界,但存在分辨率不足等问题,断溶体内幕储层结构难以有效刻画;常规的地震反演对以杂乱弱反射为主的断溶体储层预测效果不理想。
本发明通过在断裂破碎带两侧致密地层位置建立若干虚拟井,借助未钻遇储层井信息建立反映地层骨架的波阻抗模型,最后将实钻井参与的地震反演结果和地层格架阻抗求差,该残差突出断裂破碎带低阻抗信息,可有效刻画断裂破碎带储层,对指导深层断控储层的勘探开发具有重要意义。
发明内容
本发明的目的在于提供一种基于虚拟井的断溶体储层刻画方法,以解决断溶体储层内部反射弱、地层界面反射影响大、储层预测难度大等问题,为断控缝洞油藏的勘探开发提供重要依据。
为了实现以上目的,本发明采用以下技术方案:
本发明提供一种基于虚拟井的断溶体储层刻画方法,该断溶体储层刻画方法包括以下步骤:
S100、选取基准井,在断裂两侧的非储层位置建立虚拟井;
S200、将所述基准井的声波曲线赋予所述虚拟井;利用所述虚拟井进行地震反演,以此建立反映地层格架的波阻抗模型;
S300、以目标区内所有已钻井参与进行地震反演;
S400、求取S300得到的地震反演结果和所述反映地层格架的波阻抗模型的残差,刻画所述断溶体储层。
本发明通过在断裂破碎带两侧致密地层位置建立若干虚拟井,借助未钻遇储层井信息建立反映地层骨架的波阻抗模型,最后将实际井参与的地震反演结果和地层格架阻抗求残差,突出断裂破碎带低阻抗信息,可有效刻画断裂破碎带储层,对指导深层断控储层的勘探开发具有重要意义。
以下针对每一步骤进行详细说明:
S100、选取基准井,在断裂两侧的非储层位置建立虚拟井。
目标区内至少包括钻遇缝洞储层和未钻遇储层两种类型已钻井,根据本发明的断溶体储层刻画方法,优选的,以未钻遇所述断溶体储层的已钻井为基准井。若存在多口未钻遇储层的已钻井,选取反映地层格架的典型井,例如电测曲线特征清楚、完整,井震关系优良的。
根据本发明的断溶体储层刻画方法,优选的,根据断裂规模与长度,在断裂两侧非储层位置每10km各建立5~8口虚拟井。
根据本发明的断溶体储层刻画方法,优选的,所述虚拟井均匀分布在断裂破碎带两侧,且垂直断裂破碎带距离大于200米,根据相干和不连续性属性所描述断裂及破碎带位置。
S200、将所述基准井的声波曲线赋予所述虚拟井;利用所述虚拟井进行地震反演,以此建立反映地层格架的波阻抗模型。
根据本发明的断溶体储层刻画方法,优选的,S200中利用所述虚拟井的声波曲线进行地层格架地震反演,虚拟井含有非储层信息,由于受未钻遇储层虚拟井约束,反演的结果主要包含地层格架的波阻抗数据L。
S300、以目标区内所有已钻井参与进行地震反演。
此处参与进行地震反演的已钻井包括钻遇储层和未钻遇储层的已钻井。
根据本发明的断溶体储层刻画方法,优选的,S S300中所述利用目标区内所有已钻井参与进行地震反演计算,获取包含地层格架和储层信息的波阻抗数据M。
S400、将S300得到的地震反演结果与所述反映地层格架的波阻抗模型做差,以残差刻画所述断溶体储层。
S300中反演的结果阻抗包含断溶体破碎带的阻抗信息和地层格架的波阻信息,断裂破碎带总体阻抗较低,但破碎带低阻抗信息淹没在地层格架阻抗中,很难进行有效刻画;与所述反映地层格架的波阻抗模型做差,消除致密地层格架波阻抗的影响,仅保留断裂破碎带低阻抗残余信息,突出断溶体储层。
根据本发明的断溶体储层刻画方法,优选的,S400具体包括:将所述波阻抗数据M与波阻抗数据L做差,刻画所述断溶体储层。
根据本发明的断溶体储层刻画方法,优选的,所述断溶体储层刻画方法包括以下步骤:
S100、选取目标区内未钻遇所述断溶体储层的实钻井作为基准井,在断裂两侧的非储层位置分别建立若干虚拟井;
S200、将所述基准井的声波曲线赋予所述虚拟井;利用所述虚拟井的声波曲线进行地层格架地震反演,建立反映地层格架的波阻抗模型;
S300、以目标区内所有已钻井参与进行地震反演;反演的结果阻抗包含断溶体破碎带的阻抗信息和地层格架的波阻信息,断裂破碎带总体阻抗较低,但破碎带低阻抗信息淹没在地层格架阻抗中,很难进行有效刻画;
S400、将S300得到的地震反演结果与所述反映地层格架的波阻抗模型做差,消除致密地层格架波阻抗的影响,仅保留断裂破碎带低阻抗残余信息,突出断溶体储层,完成刻画所述断溶体储层。
优选的,S100中根据断裂规模与长度,在断裂两侧非储层位置每10km各建立5~8口虚拟井;所述虚拟井均匀分布在断裂破碎带两侧,与断裂破碎带垂直距离大于200米。
本发明提供的断溶体储层刻画方法包括:选取基准井在断裂两侧建立若干虚拟井作为控制点、利用虚拟井进行反演以建立地层格架波阻抗模型、实际井参与进行地震反演、地震反演与地层格架反演模型做差进行断溶体储层刻画四个关键步骤,解决了断溶体储层内部反射弱、地层界面反射影响大、储层预测难度大等问题,为断控缝洞油藏的勘探开发提供重要依据。
附图说明
图1为本发明应用例中布置的虚拟井平面示意图。
图2为本发明应用例中布置的虚拟井剖面示意图。
图3为本发明应用例中建立的反映断裂破碎带致密地层格架的波阻抗模型。
图4为本发明应用例中利用实钻井进行波阻抗反演结果图。
图5为本发明应用例中断溶体储层的刻画结果。
图6a为对比案例中过S102井地震剖面图。
图6b为对比案例中过S102井常规反演剖面。
图6c为对比案例中过S102井采用本发明的基于虚拟井反演剖面。
具体实施方式
为了更清楚地说明本发明,下面结合优选实施例对本发明做进一步的说明。本领域技术人员应当理解,下面所具体描述的内容是说明性的而非限制性的,不应以此限制本发明的保护范围。
应用实例:
采用本发明的刻画方法对塔里木盆地某区断控缝洞储层进行预测与刻画,具体包括以下步骤:
S100、选取基准井,在断裂两侧的非储层位置建立虚拟井。
断控缝洞储层具有沿走滑断裂局限发育的特征,根据地震属性可获得断溶体储层发育的平剖基本形态,实钻井C未钻遇缝洞储层,钻遇地层均为致密灰岩,测井曲线完整,典型曲线反映非储层地层结构,实钻井D在断裂破碎带钻遇缝洞储层。
选取实钻井C为基准井,如图1和图2所示,在断控储层两侧致密灰岩位置建立虚拟井A1-A7、B1-B6,利用实钻井D和虚拟井进行控制,为后续断控储层刻画做准备。
S200、将所述基准井的声波曲线赋予所述虚拟井;利用所述虚拟井进行地震反演,以此建立反映地层格架的波阻抗模型。
将未钻遇缝洞储层实钻井C声波曲线赋予虚拟井,利用虚拟井A1、A2、A3、A4、A5、A5、A6、A7以及B1、B2、B3、B4、B5、B6的声波曲线进行地层格架地震反演,建立反映断裂破碎带致密地层格架的波阻抗模型;如图3所示。
S300、以目标区内实钻井参与进行地震反演。
实钻井C和实钻井D参与进行地震反演,利用分形分维算法,建立适合碳酸盐岩缝洞储层非均质特征的地质模型,在此基础上进行波阻抗反演,结果如图4所示,反演的结果阻抗包含断溶体破碎带的阻抗信息和地层格架的波阻信息,断裂破碎带总体阻抗较低,但破碎带低阻抗信息淹没在地层格架阻抗中,很难进行有效刻画。
S400、将S300得到的地震反演结果与所述反映地层格架的波阻抗模型做差,刻画所述断溶体储层。
实钻井参与地震反演与地层格架反演模型做差,消除致密地层格架波阻抗的影响,仅保留断裂破碎带低阻抗残余信息,突出断溶体储层,完成断溶体储层的刻画,结果如图5所示,可见断溶体储层主要受控于断裂破碎带,边界轮廓比较清楚,且断溶体内幕储层强弱关系符合地质认识,主断裂储层相对较发育。
对比案例:
如图6a-图6c所示,其中图6a为过S106井地震剖面图,图6b为常规反演剖面,图6c为采用本发明的基于虚拟井反演剖面;断溶体内幕反射较弱,常规反演方法只能对强反射储层进行刻画,同时反演结果受地层界面影响,而本发明基于虚拟井反演方法可同时刻画断溶体边界及内幕储层特征。
显然,本发明的上述实施例仅仅是为清楚地说明本发明所作的举例,而并非是对本发明的实施方式的限定,对于所属领域的普通技术人员来说,在上述说明的基础上还可以做出其它不同形式的变化或变动,这里无法对所有的实施方式予以穷举,凡是属于本发明的技术方案所引伸出的显而易见的变化或变动仍处于本发明的保护范围之列。

Claims (5)

1.一种基于虚拟井的断溶体储层刻画方法,其特征在于,所述断溶体储层刻画方法包括以下步骤:
S100、选取基准井,在断裂两侧的非储层位置建立虚拟井;断控缝洞储层具有沿走滑断裂局限发育的特征,根据地震属性获得断溶体储层发育的平剖基本形态,实钻井C未钻遇缝洞储层,钻遇地层均为致密灰岩,测井曲线完整,典型曲线反映非储层地层结构,实钻井D在断裂破碎带钻遇缝洞储层;
选取目标区内未钻遇所述断溶体储层的实钻井C作为基准井,在断控储层两侧致密灰岩位置分别建立若干虚拟井;
利用实钻井D和虚拟井进行控制,为后续断控储层刻画做准备;
S200、将所述基准井的声波曲线赋予所述虚拟井;利用所述虚拟井的声波曲线进行地层格架地震反演,建立反映断裂破碎带致密地层格架的波阻抗模型;
S300、以目标区内所有已钻井参与进行地震反演:实钻井C和实钻井D参与进行地震反演,利用分形分维算法,建立适合碳酸盐岩缝洞储层非均质特征的地质模型,在此基础上进行波阻抗反演,反演的结果阻抗包含断溶体破碎带的阻抗信息和地层格架的波阻信息,断裂破碎带总体阻抗较低;
S400、将步骤S300得到的地震反演结果与所述反映断裂破碎带致密地层格架的波阻抗模型做差,消除致密地层格架波阻抗的影响,仅保留断裂破碎带低阻抗残余信息,突出断溶体储层,完成断溶体储层的刻画。
2. 根据权利要求1所述的断溶体储层刻画方法,其特征在于,步骤S100中根据断裂规模与长度,在断裂两侧非储层位置每10 km各建立5~8口虚拟井;所述虚拟井均匀分布在断裂破碎带两侧,与断裂破碎带垂直距离大于200米。
3.根据权利要求1所述的断溶体储层刻画方法,其特征在于,步骤S200中利用所述虚拟井的声波曲线进行地层格架地震反演,由于受未钻遇储层虚拟井约束,反演的结果主要包含地层格架的波阻抗数据L。
4.根据权利要求3所述的断溶体储层刻画方法,其特征在于,步骤S300中利用目标区内所有已钻井参与进行地震反演计算,获取包含地层格架和储层信息的波阻抗数据M。
5.根据权利要求4所述的断溶体储层刻画方法,其特征在于,步骤S400具体包括:将所述波阻抗数据M与波阻抗数据L做差,刻画所述断溶体储层。
CN202011102576.0A 2020-10-15 2020-10-15 一种基于虚拟井的断溶体储层刻画方法 Active CN112505760B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011102576.0A CN112505760B (zh) 2020-10-15 2020-10-15 一种基于虚拟井的断溶体储层刻画方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011102576.0A CN112505760B (zh) 2020-10-15 2020-10-15 一种基于虚拟井的断溶体储层刻画方法

Publications (2)

Publication Number Publication Date
CN112505760A CN112505760A (zh) 2021-03-16
CN112505760B true CN112505760B (zh) 2024-05-28

Family

ID=74953869

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011102576.0A Active CN112505760B (zh) 2020-10-15 2020-10-15 一种基于虚拟井的断溶体储层刻画方法

Country Status (1)

Country Link
CN (1) CN112505760B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115822563B (zh) * 2022-12-30 2023-06-09 中国石油化工股份有限公司 一种超深层断控缝洞型油气藏高产井设计方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5798982A (en) * 1996-04-29 1998-08-25 The Trustees Of Columbia University In The City Of New York Method for inverting reflection trace data from 3-D and 4-D seismic surveys and identifying subsurface fluid and pathways in and among hydrocarbon reservoirs based on impedance models
CN101980053A (zh) * 2010-10-18 2011-02-23 中国石油化工股份有限公司 一种复杂礁滩储层预测方法
CN105527653A (zh) * 2015-12-01 2016-04-27 中国地质大学(北京) 一种基于地质信息的虚拟井构建方法
CN106354895A (zh) * 2015-07-16 2017-01-25 中国石油化工股份有限公司 断块油藏多层砂岩储层层间穿层建模方法
CN109387867A (zh) * 2017-08-10 2019-02-26 中国石油化工股份有限公司 一种致密砂岩储层建模方法
CN110794476A (zh) * 2018-08-01 2020-02-14 中国石油化工股份有限公司 一种基于断溶体相控的反演方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5798982A (en) * 1996-04-29 1998-08-25 The Trustees Of Columbia University In The City Of New York Method for inverting reflection trace data from 3-D and 4-D seismic surveys and identifying subsurface fluid and pathways in and among hydrocarbon reservoirs based on impedance models
CN101980053A (zh) * 2010-10-18 2011-02-23 中国石油化工股份有限公司 一种复杂礁滩储层预测方法
CN106354895A (zh) * 2015-07-16 2017-01-25 中国石油化工股份有限公司 断块油藏多层砂岩储层层间穿层建模方法
CN105527653A (zh) * 2015-12-01 2016-04-27 中国地质大学(北京) 一种基于地质信息的虚拟井构建方法
CN109387867A (zh) * 2017-08-10 2019-02-26 中国石油化工股份有限公司 一种致密砂岩储层建模方法
CN110794476A (zh) * 2018-08-01 2020-02-14 中国石油化工股份有限公司 一种基于断溶体相控的反演方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
虚拟井技术在琼东南深水区储层预测中的应用;何丽娟 等;物探与化探;第40卷(第2期);390-397 *

Also Published As

Publication number Publication date
CN112505760A (zh) 2021-03-16

Similar Documents

Publication Publication Date Title
CN109116428B (zh) 一种缝洞型碳酸盐岩油藏不确定性建模方法及其装置
CN109388817B (zh) 一种储层裂缝三维建模方法
CN107390264B (zh) 碳酸盐岩断溶体内部结构的表征方法
CN106014399B (zh) 一种非均质地层高精度三维地应力模型建立方法
CN110308487B (zh) 一种断溶体型油藏定量表征方法
CN110309518B (zh) 缝洞型碳酸盐岩油藏溶蚀孔洞分级分类建模方法
CN103403768B (zh) 有关地下岩层的模型的方法和系统
Che et al. Three-dimensional geological modeling of coal seams using weighted Kriging method and multi-source data
CN103472484A (zh) 基于rs三维敏感地震属性分析的水平井轨迹优化方法
CN112394404B (zh) 一种渐进式储层精细表征方法
CN107558999B (zh) 一种孤立缝洞体剩余油定量预测方法
CN111612899B (zh) 一种碳酸盐岩缝洞型油藏地质建模方法
CN106894761A (zh) 利用时间域地震体的地质导向模型的地质导向方法
CN115166853A (zh) 页岩气藏天然裂缝模型建立方法、装置、电子设备及介质
CN111695228A (zh) 一种缝洞型碳酸盐岩油藏的多尺度裂缝建模方法
CN110308497B (zh) 缝洞型油藏不同类型储集体模型融合方法
CN112505760B (zh) 一种基于虚拟井的断溶体储层刻画方法
CN105404972A (zh) 一种油藏开发不确定性研究及风险控制的方法
CN113779811B (zh) 一种构造变形对现今地应力方位扰动的定量分析方法
CN115877447A (zh) 一种直平联合井网条件下地震约束三维地质建模的储层预测方法
CN111190228B (zh) 一种海底扇水道储层构型识别方法
WO2023130074A1 (en) Geologic modeling framework
CN114332391A (zh) 一种三维地质体概率模型建模方法
CN114153002A (zh) 储层天然裂缝三维地质建模方法、装置、电子设备及介质
CN114114452A (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