CN107728227B - 一种井网成熟区快速判别地下河道砂体方法 - Google Patents

一种井网成熟区快速判别地下河道砂体方法 Download PDF

Info

Publication number
CN107728227B
CN107728227B CN201710973010.7A CN201710973010A CN107728227B CN 107728227 B CN107728227 B CN 107728227B CN 201710973010 A CN201710973010 A CN 201710973010A CN 107728227 B CN107728227 B CN 107728227B
Authority
CN
China
Prior art keywords
microfacies
well
single layer
river
sand
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
CN201710973010.7A
Other languages
English (en)
Other versions
CN107728227A (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 CN201710973010.7A priority Critical patent/CN107728227B/zh
Publication of CN107728227A publication Critical patent/CN107728227A/zh
Application granted granted Critical
Publication of CN107728227B publication Critical patent/CN107728227B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V9/00Prospecting or detecting by methods not provided for in groups G01V1/00 - G01V8/00
    • G01V9/007Prospecting or detecting by methods not provided for in groups G01V1/00 - G01V8/00 by detecting gases or particles representative of underground layers at or near the surface

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Revetment (AREA)

Abstract

本发明一种井网成熟区快速判别地下河道砂体方法,包括步骤1),地层划分和对比,建立等时地层对比格架;步骤2),绘制各个单层的砂体厚度等值线图;步骤3),在等时地层对比格架约束下,根据测井曲线形态和岩屑资料开展沉积微相解释;步骤4),以单层为基本单位,建立沉积微相数据库;步骤5),以任意一单层为研究对象,逐井开展河道砂体识别;步骤6)在作为研究对应的单层所有的井都参与对比后,则进行下一单层的河道砂体识别,直至所有单层全部完成河道砂体识别;步骤7),将任意一单层所有井的识别结果按照井位坐标投影到该单层的砂体厚度等值线图上去,以砂体厚度等值线图为约束,对识别出的地下河道砂体进行平面组合。

Description

一种井网成熟区快速判别地下河道砂体方法
技术领域
本发明涉及油气田储层地质研究领域,具体为一种井网成熟区快速判别地下河道砂体方法。
背景技术
油气田开发后期,井网比较成熟,基于剩余油(气)分布规律研究的提高采收率技术是油气地质学家面临的最关键挑战,而剩余油(气)分布规律研究的基础是以地下河道识别为核心储层精细描述。
油气田开发后期,储层精细描述可以利用的资料丰富,包括钻井、测井、录井、地震等,在井网加密区块描述精度可以达到单砂体级别,丰富的描述成果为剩余油(气)分布预测及提高采收率措施调整奠定了坚实的基础。
然而需要注意的是,传统的储层精细描述工作研究过程必须先对研究工区建立剖面网格,在单层划分与对比基础上,开展沉积微相解释,然后逐一对各个剖面各个单层进行河道识别,用砂体厚度等值线图进行平面控制,两者结合预测河道空间分布。虽然该过程可以借助计算机相关软件进行,然而还是存在一些效率低下的问题,主要表现在:
(1)剖面网格建立
剖面网格建立因人而异,不同的地质研究人员建立的剖面网格可能会有所不同,这会对研究结果产生一定的影响。
(2)河道识别
河道砂体的识别必须先在剖面中进行识别,识别过程复杂,主要体现在①必须逐层进行,每完成一个单层,才能进行下一个单层,且在识别过程中必须要对单层进行顶拉平;②在针对每一个单层进行河道识别过程中,必须对该单层内所建立的所有的剖面进行识别,过程十分繁琐;③河道识别通常需要将所有河道完整的识别出来,包括河道主体、河道边部,然后再根据砂体厚度等值线图进行平面约束,预测河道空间分布,这意味着在剖面、平面两个角度都需要进行河道预测,显然这产生了一定量重复性工作。
(3)研究区域井数
一般的油气田在开发后期,井数众多,有的多达上万口井,这进一步加剧了研究工作的繁琐性,工作效率低下。因此,亟需开发一种快速、高效的油气田开发后期快速判别地下河道砂体方法。
发明内容
针对现有技术中存在的问题,本发明提供一种井网成熟区快速判别地下河道砂体方法,工作效率高,判断速度快,操作简便,判断准确。
本发明是通过以下技术方案来实现:
一种井网成熟区快速判别地下河道砂体方法,包括如下步骤;
步骤1),地层划分和对比,对研究区所有井地层划分至单层级别,建立等时地层对比格架;
步骤2),绘制各个单层的砂体厚度等值线图;
步骤3),在等时地层对比格架约束下,根据测井曲线形态和岩屑资料开展沉积微相解释,识别出河道微相、泛滥平原微相和溢岸沉积微相三种类型;
步骤4),以单层为基本单位,建立沉积微相数据库;每一口井的一个单层数据包括单层名称、井名、所有微相类型和各微相的顶部深度及底部深度;
步骤5),以任意一单层为研究对象,逐井开展河道砂体识别;其对每一口井的判断标准如下,
若仅识别出河道微相,则判断该井处于地下河道砂体内部;
若仅识别出溢岸微相,则判断该井处于地下河道砂体外部;
若识别出河道微相和泛滥平原微相两种微相,则河道微相顶部深度等于泛滥平原微相底部深度,判断该井处于地下河道砂体内部;
若识别出河道微相和溢岸微相两种微相,则溢岸微相顶部深度等于河道微相底部深度,判断该井处于地下河道砂体边部;
若识别出泛滥平原微相和溢岸微相两种微相,则溢岸微相顶部深度等于泛滥平原微相底部深度,判断该井处于地下河道砂体外部;
若识别出河道微相、泛滥平原微相和溢岸微相三种微相,则溢岸微相顶部深度等于河道微相底部深度,河道微相顶部深度等于泛滥平原微相底部深度,判断该井处于地下河道砂体边部;
步骤6)在作为研究对应的单层所有的井都参与对比后,则进行下一单层的河道砂体识别,直至所有单层全部完成河道砂体识别;
步骤7),将任意一单层所有井的识别结果按照井位坐标投影到该单层的砂体厚度等值线图上去,以砂体厚度等值线图为约束,对识别出的地下河道砂体进行平面组合,完成井网成熟区地下河道砂体的快速判断。
优选的,步骤1)中单层为垂向由单一沉积单元构成的地层单位;将不同井的相同时间内的沉积单元连起来即建立等时地层对比格架。
优选的,步骤2)中绘制各个单层的砂体厚度等值线图时,砂体厚度数据来源于测井资料的解释结果。
优选的,步骤4)中建立沉积微相数据库时;
首先对沉积微相进行编号,河道微相编号为Ⅰ,泛滥平原微相编号为Ⅱ,溢岸微相编号为Ⅲ;
其次以单层为单位,对所有井识别出的沉积微相进行数据入库。
优选的,步骤7)中进行平面组合时,对单层中的井进行如下判断;
若井的位置属于河道砂体内部,则其必然落在砂体厚度等值线图中厚度大的区域,厚度大的区域是指砂体厚度大于该单层等值线图中最大厚度的75%,小于等于该单层等值线图中的最大厚度;
若井的位置属于河道砂体边部,则其必然落在砂体厚度等值线图中厚度较大区域的左边或右边;厚度较大区域是指砂体厚度大于该单层等值线图中最大厚度的25%,小于等于该单层等值线图中最大厚度的75%;
若井的位置属于河道砂体边部,则其必然落在砂体厚度等值线图中厚度极小的区域;厚度极小的区域是指砂体厚度大于等于该单层等值线图中最大厚度的0%,小于等于该单层等值线图中最大厚度的25%。
进一步,若井的位置属于河道砂体边部,如果其在厚度等值线图中位置落在厚度较大区域的左边,则属于河道砂体的左边界。如果其在厚度等值线图中位置落在厚度较大区域的右边,则属于河道砂体的右边界。
与现有技术相比,本发明具有以下有益的技术效果:
本发明进行河道砂***置判断时,不需要依靠剖面网格,仅仅依靠单井即可完成;以单层中的单个井为研究对象,仅仅通过对单个井的沉积微相叠置关系进行分析,便可所处的河道砂体进行快速判断,不需要对单层进行顶拉平处理;同时可以灵活进行识别顺序调整,既可以对一口井进行所有单层识别后,再识别下一口井,也可以对某一个单层内所有井进行识别后,再对下一个单层进行识别,而不必像传统方法那样,必须先要对某一个单层内所有井进行识别后,才能进行下一个单层进行识别;油气田开发后期,井数众多,传统河道识别工作繁琐,效率低下,本发明方法,能够高效的识别地下河道砂体,同时操作简便。
附图说明
图1是本发明一种井网成熟区快速判别地下河道砂体方法流程图。
图2是本发明实例中所述利用井A的沉积微相进行河道砂体识别的示意图。
图3是本发明实例中所述利用井B和井B'的沉积微相进行河道砂体识别的示意图。
图4是本发明实例中所述利用井C的沉积微相进行河道砂体识别的示意图。
图5是本发明实例中所述利用井D和井D'的沉积微相进行河道砂体识别的示意图。
图6是本发明实例中所述利用井E的沉积微相进行河道砂体识别的示意图。
图7是本发明实例中所述利用井F的沉积微相进行河道砂体识别的示意图。
具体实施方式
下面结合具体的实施例对本发明做进一步的详细说明,所述是对本发明的解释而不是限定。
本发明一种井网成熟区快速判别地下河道砂体方法,用于油气田开发后期开发井数很多的区块,对河道砂体边界进行快速判断。如图1所示,其主要包括如下步骤。
1)地层划分和对比,对研究区所有井地层划分至单层级别,建立等时地层对比格架。
步骤1)中单层是指垂向由单一沉积单元构成的地层单位;将不同井的相同时间内的沉积单元连起来即建立等时地层对比格架。
2)绘制各个单层的砂体厚度等值线图。
步骤2)砂体厚度数据来源于测井资料解释结果。
3)依据测井曲线形态,对每一口井开展沉积微相解释。
步骤3)中沉积微相即沉积环境与沉积物的综合;沉积微相解释是指根据测井曲线形态、岩屑资料识别在等时地层对比格架约束下,识别河道、泛滥平原、溢岸沉积微相。每一口井某一个单层内微相的种类可能有3种情况:①只存在1种微相;②存在2种微相;③存在3种微相。
4)以单层为基本单位,建立沉积微相数据库。
步骤4)中首先对沉积微相进行编号,河道微相编号为Ⅰ,泛滥平原微相编号为Ⅱ,溢岸微相编号为Ⅲ;其次以单层为单位,对所有井识别出的沉积微相进行数据入库。如A井a单层,识别出Ⅰ、Ⅱ、Ⅲ三种微相,则分别记录单层名称a、井名A、微相Ⅰ、微相Ⅰ顶部深度、微相Ⅰ底部深度,微相Ⅱ、微相Ⅱ顶部深度、微相Ⅱ底部深度,微相Ⅲ、微相Ⅲ顶部深度、微相Ⅲ底部深度。
5)以某一单层为研究对象,逐井开展河道砂体识别;
步骤5)判断标准是:假设在a单层井A;
若仅识别出微相Ⅰ,则认为该井处于地下河道砂体内部;
若仅识别出微相Ⅲ,则认为该井处于地下河道砂体外部;
若识别出Ⅰ和Ⅱ两种微相,则微相Ⅰ顶部深度=微相Ⅱ底部深度,认为该井处于地下河道砂体内部;
若识别出Ⅰ和Ⅲ两种微相,则微相Ⅲ顶部深度=微相Ⅰ底部深度,则认为该井处于地下河道砂体边部;
若识别出Ⅱ和Ⅲ两种微相,则微相Ⅲ顶部深度=微相Ⅱ底部深度,则认为该井处于地下河道砂体外部;
若识别出Ⅰ、Ⅱ和Ⅲ三种微相,则微相Ⅲ顶部深度=微相Ⅰ底部深度,微相Ⅰ顶部深度=微相Ⅱ底部深度,则认为该井处于地下河道砂体边部;
具体应用后分别对如下井进行判断。
如图2所示,井A只在该单层a只钻遇了河道微相Ⅰ,表明该井处于河道砂体内部;
如图3所示,井B和井B'在该单层a垂向微相发育特征是:河道微相叠置在溢岸微相之上,表明该井处于河道砂体边部;
如图4所示,井C在该单层a垂向微相发育特征是:泛滥平原微相叠置在河道微相之上,表明该井处于河道砂体内部;
如图5所示,井D和井D'在该单层垂向微相发育特征是:泛滥平原微相叠置在河道微相之上,同时河道微相叠置在溢岸微相之上,表明该井处于河道砂体边部;
如图6所示,井E在该单层只钻遇了溢岸微相,表明该井处于河道砂体外部;
如图7所示,井F在该单层垂向微相发育特征是:泛滥平原微相叠置在溢岸微相之上,表明该井处于河道砂体外部。
6)在该单层所有的井都参与对比后,则进行下一单层的河道砂体识别,直至所有单层全部进行河道砂体识别;
7)将某单层识别结果投影到该单层的砂体厚度等值线图上去,参考砂体厚度等值线图,对识别出的地下河道砂体进行平面组合。
步骤7)中平面组合方法是指,假设在a单层井A;
若属于河道砂体内部,则其必然落在砂体厚度等值线图中厚度大的区域,厚度大的区域是指砂体厚度大于该单层等值线图中最大厚度的75%,小于等于该单层等值线图中的最大厚度;
若属于河道砂体边部,如果其在厚度等值线图中位置落在厚度较大区域的左边,则属于河道砂体的左边界。如果其在厚度等值线图中位置落在厚度较大区域的右边,则属于河道砂体的右边界。厚度较大区域是指砂体厚度大于该单层等值线图中最大厚度的25%,小于等于该单层等值线图中最大厚度的75%;
若属于河道砂体边部,则其必然落在砂体厚度等值线图中厚度极小(甚至为零)的区域;厚度极小的区域是指砂体厚度大于等于该单层等值线图中最大厚度的0%,小于等于该单层等值线图中最大厚度的25%。
本发明所述的方法在苏里格气田某区块(700口井)的应用下,工作效率提高了30%,识别结果应用于三维地质建模中作为约束条件,使得模型符合率提高了8%,效果明显。
以上内容仅为说明本发明的技术思想,不能以此限定本发明的保护范围,凡是按照本发明提出的技术思想,在技术方案基础上所做的任何改动,均落入本发明权利要求书的保护范围之内。

Claims (6)

1.一种井网成熟区快速判别地下河道砂体方法,其特征在于,包括如下步骤;
步骤1),地层划分和对比,对研究区所有井地层划分至单层级别,建立等时地层对比格架;
步骤2),绘制各个单层的砂体厚度等值线图;
步骤3),在等时地层对比格架约束下,根据测井曲线形态和岩屑资料开展沉积微相解释,识别出河道微相、泛滥平原微相和溢岸沉积微相三种类型;
步骤4),以单层为基本单位,建立沉积微相数据库;每一口井的一个单层数据包括单层名称、井名、步骤3)中的三种微相类型和各微相的顶部深度及底部深度;
步骤5),以任意一单层为研究对象,逐井开展河道砂体识别;其对每一口井的判断标准如下,
若仅识别出河道微相,则判断该井处于地下河道砂体内部;
若仅识别出溢岸微相,则判断该井处于地下河道砂体外部;
若识别出河道微相和泛滥平原微相两种微相,则河道微相顶部深度等于泛滥平原微相底部深度,判断该井处于地下河道砂体内部;
若识别出河道微相和溢岸微相两种微相,则溢岸微相顶部深度等于河道微相底部深度,判断该井处于地下河道砂体边部;
若识别出泛滥平原微相和溢岸微相两种微相,则溢岸微相顶部深度等于泛滥平原微相底部深度,判断该井处于地下河道砂体外部;
若识别出河道微相、泛滥平原微相和溢岸微相三种微相,则溢岸微相顶部深度等于河道微相底部深度,河道微相顶部深度等于泛滥平原微相底部深度,判断该井处于地下河道砂体边部;
步骤6)在作为研究对应的单层所有的井都参与对比后,则进行下一单层的河道砂体识别,直至所有单层全部完成河道砂体识别;
步骤7),将任意一单层所有井的识别结果按照井位坐标投影到该单层的砂体厚度等值线图上去,以砂体厚度等值线图为约束,对识别出的地下河道砂体进行平面组合,完成井网成熟区地下河道砂体的快速判断。
2.根据权利要求1所述的一种井网成熟区快速判别地下河道砂体方法,其特征在于,步骤1)中单层为垂向由单一沉积单元构成的地层单位;将不同井的相同时间内的沉积单元连起来即建立等时地层对比格架。
3.根据权利要求1所述的一种井网成熟区快速判别地下河道砂体方法,其特征在于,步骤2)中绘制各个单层的砂体厚度等值线图时,砂体厚度数据来源于测井资料的解释结果。
4.根据权利要求1所述的一种井网成熟区快速判别地下河道砂体方法,其特征在于,步骤4)中建立沉积微相数据库时;
首先对沉积微相进行编号,河道微相编号为Ⅰ,泛滥平原微相编号为Ⅱ,溢岸微相编号为Ⅲ;
其次以单层为单位,对所有井识别出的沉积微相进行数据入库。
5.根据权利要求1所述的一种井网成熟区快速判别地下河道砂体方法,其特征在于,步骤7)中进行平面组合时,对单层中的井进行如下判断;
若井的位置属于河道砂体内部,则其必然落在砂体厚度等值线图中厚度大的区域,厚度大的区域是指砂体厚度大于该单层等值线图中最大厚度的75%,小于等于该单层等值线图中的最大厚度;
若井的位置属于河道砂体边部,则其必然落在砂体厚度等值线图中厚度较大区域的左边或右边;厚度较大区域是指砂体厚度大于该单层等值线图中最大厚度的25%,小于等于该单层等值线图中最大厚度的75%;
若井的位置属于河道砂体边部,则其必然落在砂体厚度等值线图中厚度极小的区域;厚度极小的区域是指砂体厚度大于等于该单层等值线图中最大厚度的0%,小于等于该单层等值线图中最大厚度的25%。
6.根据权利要求5所述的一种井网成熟区快速判别地下河道砂体方法,其特征在于,若井的位置属于河道砂体边部,如果其在厚度等值线图中位置落在厚度较大区域的左边,则属于河道砂体的左边界;如果其在厚度等值线图中位置落在厚度较大区域的右边,则属于河道砂体的右边界。
CN201710973010.7A 2017-10-18 2017-10-18 一种井网成熟区快速判别地下河道砂体方法 Active CN107728227B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710973010.7A CN107728227B (zh) 2017-10-18 2017-10-18 一种井网成熟区快速判别地下河道砂体方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710973010.7A CN107728227B (zh) 2017-10-18 2017-10-18 一种井网成熟区快速判别地下河道砂体方法

Publications (2)

Publication Number Publication Date
CN107728227A CN107728227A (zh) 2018-02-23
CN107728227B true CN107728227B (zh) 2019-03-12

Family

ID=61211918

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710973010.7A Active CN107728227B (zh) 2017-10-18 2017-10-18 一种井网成熟区快速判别地下河道砂体方法

Country Status (1)

Country Link
CN (1) CN107728227B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109441438A (zh) * 2018-10-18 2019-03-08 中国石油化工股份有限公司 一种构造复杂区岩性边界表征方法
CN111624677B (zh) * 2019-02-27 2022-07-26 中国石油化工股份有限公司 一种致密砂岩气藏单期河道砂体边界刻画方法
CN110376661B (zh) * 2019-07-19 2020-10-30 中国地质大学(北京) 一种地下井间单河道砂体的对比方法和装置
CN110656934B (zh) * 2019-10-08 2024-04-30 中国石油天然气股份有限公司 一种致密砂岩储层去压实地层对比方法
CN112684512B (zh) * 2021-01-05 2023-09-26 中国石油天然气股份有限公司 一种河道识别方法

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9869799B2 (en) * 2013-09-09 2018-01-16 Schlumberger Technology Corporation Object-based well correlation
CN104632200B (zh) * 2013-11-14 2017-06-16 中国石油化工股份有限公司 河流相多层及厚层油藏精细建模方法
CN106569287B (zh) * 2016-09-14 2019-02-26 中国石油化工股份有限公司 一种基于复合砂体构型模式的扇三角洲前缘储层预测方法
CN106837319B (zh) * 2016-12-27 2019-10-11 中国石油天然气股份有限公司 一种密井网条件下的地下单一古河道识别方法
CN107203651A (zh) * 2016-12-30 2017-09-26 中国石油天然气股份有限公司 一种辫状河致密砂岩气藏储层构型解剖方法
CN108119133B (zh) * 2017-11-13 2021-09-28 中国石油天然气股份有限公司 储层沉积过程的确定方法和装置

Also Published As

Publication number Publication date
CN107728227A (zh) 2018-02-23

Similar Documents

Publication Publication Date Title
CN107728227B (zh) 一种井网成熟区快速判别地下河道砂体方法
RU2661489C1 (ru) Способ комплексирования исходных данных для уточнения фильтрационного строения неоднородных карбонатных коллекторов
CN101236257A (zh) 油井位置确定技术方法
CN105701319B (zh) 一种水平井规则开发井网下的沉积微相建模方法
CN104049275A (zh) 复杂盆缘超剥带隐蔽圈闭边界的识别方法及识别***
CN111624677B (zh) 一种致密砂岩气藏单期河道砂体边界刻画方法
CN106285642B (zh) 一种基于地震资料的地层孔隙压力预测方法
CN111257962A (zh) 一种利用三维地质建模技术定位预测矿体的方法
CN112505754B (zh) 基于高精度层序格架模型的井震协同划分沉积微相的方法
CN109239783B (zh) 一种井震结合中取舍地震信息的方法
KR101131517B1 (ko) 물리검층 자료의 통계학적 해석결과를 활용한 오일샌드 지역 지하정보의 처리시스템 및 이를 이용한 오일샌드 지역 지하정보의 처리방법
CN111768503B (zh) 一种基于三维地质模型的海砂资源量估算方法
CN104047596A (zh) 三角洲前缘沉积小层对比中小断点识别的方法
CN109117534B (zh) 一种基于地质图件的断裂和裂缝预测方法
CN109441438A (zh) 一种构造复杂区岩性边界表征方法
CN114114411B (zh) 三维断层疏导性的定量确定方法及装置
CN113376692B (zh) 致密砂岩气水平井压裂改造方案优化方法及装置
CN107060744A (zh) 一种测井地质评价***
CN112285774A (zh) 一种基于多地质参数快速评价走滑断层封闭能力的方法
Amanipoor Static modeling of the reservoir for estimate oil in place using the geostatistical method
CN112069682A (zh) 一种基于三维建模的成矿预测方法
CN114594518B (zh) 基于井震交替的开发后期复杂断块精细地层对比方法
CN114415237B (zh) 基于三维地震资料的砂岩型铀矿控矿断裂识别方法及***
CN113514883B (zh) 一种断层-岩性油藏刻画方法
CN107991702A (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