CN106371139B - 一种应用层段振幅属性识别岩性油藏的方法及装置 - Google Patents

一种应用层段振幅属性识别岩性油藏的方法及装置 Download PDF

Info

Publication number
CN106371139B
CN106371139B CN201610676257.8A CN201610676257A CN106371139B CN 106371139 B CN106371139 B CN 106371139B CN 201610676257 A CN201610676257 A CN 201610676257A CN 106371139 B CN106371139 B CN 106371139B
Authority
CN
China
Prior art keywords
interval
reservoir
amplitude attribute
amplitude
attribute
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
CN201610676257.8A
Other languages
English (en)
Other versions
CN106371139A (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
Exploration and Development Research Institute of Sinopec Henan Oilfield Branch Co
Original Assignee
China Petroleum and Chemical Corp
Exploration and Development Research Institute of Sinopec Henan Oilfield Branch 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, Exploration and Development Research Institute of Sinopec Henan Oilfield Branch Co filed Critical China Petroleum and Chemical Corp
Priority to CN201610676257.8A priority Critical patent/CN106371139B/zh
Publication of CN106371139A publication Critical patent/CN106371139A/zh
Application granted granted Critical
Publication of CN106371139B publication Critical patent/CN106371139B/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. analysis, for interpretation, for correction
    • 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
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/62Physical property of subsurface

Abstract

本发明涉及一种应用层段振幅属性识别岩性油藏的方法及装置,属于油藏圈闭识别技术领域。本发明通过对已钻井进行合成记录制作并进行储层层位标定,建立目标区层系的地层格架;对不同层系分别进行储层岩石物理特征及油气水层的振幅规律统计,选取能够应用层段振幅属性刻画的储层段;优选参数提取所选取储层的层段振幅属性,并根据所提取层段振幅属性进行目标评价以确定岩性油藏的分布。本发明利用提取的层段振幅属性,通过一张平面振幅属性图就能够快速有效的发现多套层位的岩性圈闭,大大提高了工作效率并且解决了遗漏的问题,为整体的圈闭评价优选和勘探开发部署提供依据。

Description

一种应用层段振幅属性识别岩性油藏的方法及装置
技术领域
本发明涉及一种应用层段振幅属性识别岩性油藏的方法及装置,属于油藏圈闭识别技术领域。
背景技术
随着油气勘探从构造油气藏向岩性油气藏的转变,岩性油气藏已成为目前研究的主体,是许多油田勘探开发的重点油藏类型。大多数岩性油藏均具有复杂、隐蔽的特点,因此如何快速有效的识别岩性油藏是解决我们目前勘探过程中的一项非常重要的课题。目前所使用的岩性油藏识别方法通常是应用常规的沿层振幅属性,需要建立在多套目的层层位的精细地震资料解释基础上,大部分情况下受地震资料分辨率和信噪比等诸多因素的限制,合理、准确的解释受人为因素和解释方案的影响较大,存在一定的局限性。而且实际地质情况复杂,有些孤立成层性不好的砂体或者三角洲中发育的小型前积砂体在层位追踪中容易遗漏,这样的方法不仅耗时耗力,而且容易漏掉部分岩性圈闭,从而导致这类岩性油藏无法得以识别。
发明内容
本发明的目的是提供一种应用层段振幅属性识别岩性油藏的方法及装置,以解决现有岩性圈闭识别过程中应用常规沿层振幅属性耗时耗力并且容易遗漏许多岩性油藏的问题。
本发明为解决上述技术问题而提供一种应用层段振幅属性识别岩性油藏的方法,该方法包括以下步骤:
1)对目标区的地震资料进行频谱分析,选取子波对已钻井进行合成记录制作并进行储层层位标定;
2)根据实际沉积背景和地层接触关系建立目标区层系的地层格架;
3)对不同层系分别进行储层岩石物理特征油气水层的振幅规律统计,选取能够应用层段振幅属性刻画的储层层段;
4)提取层段振幅属性,并根据所提取层段振幅属性达到最大振幅值的3/5以上作为目标进行评价以确定岩性油藏的分布。
步骤1)中选取的子波为与地震数据主频相同频率的雷克子波或井旁道子波。
步骤2)建立地层格架所需要解释的层位为各层系及砂组之间的界面。
步骤3)能够应用层段振幅属性刻画的储层指的是砂泥岩具有较大速度差异的储层,即选取速度差异达到平均速度1/10以上的层段应用层段振幅属性刻画储层。
步骤4)提取层段振幅属性时需要消除地层之间由于不整合或者沉积背景不同本身存在的强阻抗差异,需要对解释层系及砂组界面通过层位时窗漂移完成,具体漂移时窗数据根据同相轴宽度确定。
步骤4)层段振幅属性的提取范围由实际储层标定位置确定。
步骤4)提取层段振幅属性的类型由储层与围岩的阻抗差异确定,当储层速度小于围岩速度时,提取层段最小振幅属性;当储层速度大于围岩速度时,提取层段最大振幅属性。
本发明还提供了一种应用层段振幅属性识别岩性油藏的装置,该装置包括层位标定模块、地层格架建立模块、储层选取模块和层段振幅属性提取模块。
所述层位标定模块用于对目标区的地震资料进行频谱分析,选取子波对已钻井进行合成记录制作及层位标定;
所述地层格架建立模块用于根据实际沉积背景和地层接触关系建立目标区层系的地层格架;
所述储层选取模块用于对不同层系分别进行储层岩石物理特征油气水层的振幅规律统计,即选取速度差异达到平均速度1/10以上的层段应用层段振幅属性刻画储层;
所述层段振幅属性提取模块用于提取层段振幅属性,并根据所提取层段振幅属性达到最大振幅值的3/5以上作为目标进行评价以确定岩性油藏的分布。
选取的子波为与地震数据主频相同频率的雷克子波或井旁道子波。
提取层段振幅属性的类型由储层与围岩的阻抗差异确定,当储层速度小于围岩速度时,提取层段最小振幅属性;当储层速度大于围岩速度时,提取层段最大振幅属性。
取层段最小振幅属性;当储层速度大于围岩速度时,提取层段最大振幅属性。
本发明的有益效果是:本发明在地层格架建立及重点目的层系岩石物理分析的基础上,针对实储层标定位置确定层段属性提取范围,根据储层与围岩的阻抗差异关系确定提取的振幅属性类型,从而实现在一张平面振幅属性图上发现多套目的层系岩性油藏的快速有效的方法,解决了应用常规的沿层振幅属性,需要针对每套储层进行精细地震资料解释,耗时耗力并且容易由于解释方案及人员等方面问题漏掉部分岩性圈闭,从而导致这类岩性油藏无法识别的问题。大大提高了工作效率并且解决了遗漏问题,为整体的圈闭评价优选和勘探开发部署提供依据。
附图说明
图1是本发明应用层段振幅属性识别岩性油藏的流程图;
图2是本发明实施例中探区已钻井的合成记录图;
图3是本发明实施例中探区地层格架连井剖面图;
图4是本发明实施例中探区目的层系储层及围岩速度对比直方图;
图5是本发明实施例中探区目的层系油水层振幅统计直方图;
图6是本发明实施例中探区目的层系具体砂层常规沿层振幅属性图;
图7是本发明实施例中探区目的层系优选层段振幅属性图。
具体实施方式
下面结合附图对本发明的具体实施方式做进一步的说明。
本发明的一种应用层段振幅属性识别岩性油藏的方法的实施例
层段振幅属性是计算某一层系顶界和底界之间的最大、最小或者均方根等类型振幅值,能够将所有强振幅区压缩在一张平面振幅属性图中,从而直观高效的反映“亮点“区所反映的岩性油藏。为此,本发明提供了一种应用层段振幅属性来识别岩性油藏的方法,该方法通过对目标区的地震资料进行频谱分析,选取子波对已钻井进行合成记录制作并进行储层层位标定;根据实际沉积背景和地层接触关系建立目标区层系的地层格架;对不同层系分别进行储层岩石物理特征油气水层的振幅规律统计,选取能够应用层段振幅属性刻画储的储层;提取所选取储层层段振幅属性,并根据所提取层段振幅属性进行目标评价以确定岩性油藏。该方法的流程如图1所示,具体实施过程如下。
1.合成记录制作及精细层位标定
对地区的地震资料进行频谱分析,选取合适的子波对已钻井进行合成记录制作。这其中应用到的声波和密度曲线首先需要进行曲线预处理,从而完成精细层位标定。
本实施例中合成记录如图2所示,对探区的地震资料进行频谱分析,地震数据有效频带为10-120Hz,主频65Hz;因此一般选取65hz的雷克子波或者井旁道子波进行合成记录制作。在合成记录道与地震道同相轴数目和能量匹配一致后,探区多口井的主要目的层储层标定在波谷中。
2.地层格架建立及地震资料解释
建立地层格架所需要解释的主要层位为各层系及主要砂组之间的界面。另外还要根据各层系实际反射特征,追踪解释具有类似反射特征及岩石物理基础层段的顶界和底界。
本实施例中在地层对比的基础上,井震结合,通过连井地震剖面确定了探区主要研究层系之间界面的层位精细解释,解释精度达到10*10,确定了地层之间的接触关系,建立了区块的地层格架(如图3所示)。应用此技术中关注目的层系顶界和底界的解释精度需达到20*20,解释范围需达到目的层发育的最大范围。
3.储层岩石物理特征及油气水层的振幅规律统计分析
储层岩石物理特征油气水层的振幅规律统计分析是指针对不同层系分别进行,选取砂泥岩具有较大速度差异,可以应用振幅属性区分储层的层系来应用层段振幅属性刻画储层,进而结合油气水层的振幅规律统计分析来最终识别油层。由于不同探区埋深及成岩作用的差异,速度基值不一样,同样的300ms米/秒的速度差异在整体速度较高地区应用的效果与整体速度值较低地区应用效果不一样,因此这里所说的较大速度差异是指通过统计及正演试验分析认为需达到工区平均速度的1/10以上就可以应用层段振幅属性来识别岩性油藏,本实施例中区块平均速度约为2500米/秒,所选取的砂岩速度小于围岩速度差异在300米/秒以上,大于平均速度的1/10。
本实施例中具体统计层段划分主要是依据相邻时期发育的具有相似沉积背景的层段根据已钻井的实际测井曲线数据来统计分析。通过储层岩石物理特征油气水层的振幅规律统计分析,发现目的层砂岩速度小于泥岩速度(如图4所示),这与合成记录中储层标定在波谷中是一致的。储层含油之后速度降低(如图5所示),导致发现的岩性油藏在地震剖面和振幅属性平面图上表现为强振幅的亮点特征(如图6所示)。因此,可以通过层段振幅属性的应用来识别发现未知的岩性油藏。
4.优选层段的顶界及底界时窗参数和提取属性类型,提取层段振幅属性,进行目标评价寻找出岩性油藏。
由于解释层系及砂组之间的界面一般都是在波峰中,而储层是标定在波谷中,提取层段属性时候需要消除地层之间由于不整合或者沉积背景不同本身存在的强阻抗差异,因此需要对解释层系及砂组界面通过层位时窗漂移完成,具体漂移时窗数据根据同相轴宽度确定。提取层段振幅属性的类型由储层与围岩的阻抗差异确定,当储层速度小于围岩速度时,提取层段最小振幅属性;当储层速度大于围岩速度时,提取层段最大振幅属性。
振幅绝对值较大的范围同储层与围岩速度差异原理类似,由于不同的地区以及地震资料的不同处理流程实际振幅值的大小差异非常大,通过统计及正演试验分析认为需达到最大振幅值的3/5以上认为是需要评价的岩性油藏发育区。储层速度小于围岩速度时,选取振幅负值的绝对值较大处;储层速度大于围岩速度时,选取振幅正值绝对值较大处。
本实施例中目的层底界超覆于底界不整合面上,为消除不整合面之上波谷反射强阻抗影响,选取层段底界向上漂移15ms(-15ms);顶界与目的层系本身存在沉积差异,为消除顶界下方波谷的强阻抗,选取层段顶界向下漂移7ms(+7ms)。从而提取出的层段最小振幅属性平面图(如图7所示)。图7中可以看出,之前花费大量时间和工作量发现负值振幅值绝对值大于30000的亮点岩性油藏,(最大振幅绝对值为50000,大于3/5),在这张层段振幅属性图中都可以找到,如井1、井2、井3等。另外还发现了许多之前解释漏掉的一些小砂层形成的岩性圈闭,如井4、井5所处的岩性圈闭。可见本发明的识别方法能够快速有效的识别圈闭,识别出的岩性油藏经评价优选部署的井位试油获得了高产工业油流,目前每天日产20多吨的稀油。另外还储备了一些有待下步井位部署的有利圈闭,如设计井6和井7所处圈闭。这些实例都证实了本发明的高效性和先进性,能够应用于具有同类型地质条件下岩性油藏的识别,减少工作量同时提高勘探开发效益。
本发明的一种应用层段振幅属性识别岩性油藏的装置的实施例:
本发明应用层段振幅属性识别岩性油藏的该装置包括层位标定模块、地层格架建立模块、储层选取模块和层段振幅属性提取模块,层位标定模块用于对目标区的地震资料进行频谱分析,选取子波对已钻井进行合成记录制作及层位标定;地层格架建立模块用于根据实际沉积背景和地层接触关系建立目标区层系的地层格架;储层选取模块用于对不同层系分别进行储层岩石物理特征油气水层的振幅规律统计,选取能够应用层段振幅属性刻画的储层;层段振幅属性提取模块用于提取层段振幅属性,并根据所提取层段振幅属性进行目标评价以确定岩性油藏。各模块的具体实现手段已在方法的实施例中进行了详细说明,这里不再赘述。
本发明通过采用应用层段属性技术,不需要解释每套砂体,直接在一张平面振幅属性图上就能够快速有效的发现大量岩性圈闭,直接筛选锁定目标解释,减少了每个小层一一解释的庞大工作量,大大提高了工作效率并且解决了遗漏问题,为整体的圈闭评价优选和勘探开发部署提供依据,为岩性油藏的发现提供了一种高效直观的方式。

Claims (6)

1.一种应用层段振幅属性识别岩性油藏的方法,其特征在于,该方法包括以下步骤:
1)对目标区的地震资料进行频谱分析,选取子波对已钻井进行合成记录制作并进行储层层位标定;
2)根据实际沉积背景和地层接触关系建立目标区层系的地层格架;
3)对不同层系分别进行储层岩石物理特征油气水层的振幅规律统计,选取能够应用层段振幅属性刻画的储层层段;
4)提取层段振幅属性,并根据所提取层段振幅属性达到最大振幅值的3/5以上作为目标进行评价以确定岩性油藏的分布;储层速度小于围岩速度时,选取振幅负值的绝对值较大处;储层速度大于围岩速度时,选取振幅正值绝对值较大处;
步骤2)建立地层格架所需要解释的层位为各层系及砂组之间的界面;还根据各层系实际反射特征,追踪解释具有类似反射特征及岩石物理基础层段的顶界和底界;
步骤3)能够应用层段振幅属性刻画的储层指的是砂泥岩具有较大速度差异的储层,即选取速度差异达到平均速度1/10以上的层段应用层段振幅属性刻画储层;
步骤4)提取层段振幅属性的类型由储层与围岩的阻抗差异确定,实现在一张平面振幅属性图上发现多套目的层系岩性油藏;当储层速度小于围岩速度时,提取层段最小振幅属性;当储层速度大于围岩速度时,提取层段最大振幅属性;
所述的储层层段是指包含纵向上具有类似岩石物理特征的多套目的层位的砂组段;所述的岩性油藏是指以砂岩作为储层,泥岩作为围岩的油藏。
2.根据权利要求1所述的应用层段振幅属性识别岩性油藏的方法,其特征在于,步骤1)中选取的子波为与地震数据主频相同频率的雷克子波或井旁道子波。
3.根据权利要求1所述的应用层段振幅属性识别岩性油藏的方法,其特征在于,步骤4)提取层段振幅属性时需要消除地层之间由于不整合或者沉积背景不同本身存在的强阻抗差异,需要对解释层系及砂组界面通过层位时窗漂移完成,具体漂移时窗数据根据同相轴宽度确定。
4.根据权利要求1所述的应用层段振幅属性识别岩性油藏的方法,其特征在于,步骤4)层段振幅属性的提取范围由实际储层标定位置确定。
5.一种应用层段振幅属性识别岩性油藏的装置,其特征在于,该装置包括层位标定模块、地层格架建立模块、储层选取模块和层段振幅属性提取模块:
所述层位标定模块用于对目标区的地震资料进行频谱分析,选取子波对已钻井进行合成记录制作及层位标定;
所述地层格架建立模块用于根据实际沉积背景和地层接触关系建立目标区层系的地层格架;
所述储层选取模块用于对不同层系分别进行储层岩石物理特征油气水层的振幅规律统计,即选取速度差异达到平均速度1/10以上的层段应用层段振幅属性刻画储层;
所述层段振幅属性提取模块用于提取层段振幅属性,并根据所提取层段振幅属性达到最大振幅值的3/5以上作为目标进行评价以确定岩性油藏的分布;储层速度小于围岩速度时,选取振幅负值的绝对值较大处;储层速度大于围岩速度时,选取振幅正值绝对值较大处;
还根据各层系实际反射特征,追踪解释具有类似反射特征及岩石物理基础层段的顶界和底界;
提取层段振幅属性的类型由储层与围岩的阻抗差异确定,实现在一张平面振幅属性图上发现多套目的层系岩性油藏;当储层速度小于围岩速度时,提取层段最小振幅属性;当储层速度大于围岩速度时,提取层段最大振幅属性。
6.根据权利要求5所述的应用层段振幅属性识别岩性油藏的装置,其特征在于,选取的子波为与地震数据主频相同频率的雷克子波或井旁道子波。
CN201610676257.8A 2016-08-16 2016-08-16 一种应用层段振幅属性识别岩性油藏的方法及装置 Active CN106371139B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610676257.8A CN106371139B (zh) 2016-08-16 2016-08-16 一种应用层段振幅属性识别岩性油藏的方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610676257.8A CN106371139B (zh) 2016-08-16 2016-08-16 一种应用层段振幅属性识别岩性油藏的方法及装置

Publications (2)

Publication Number Publication Date
CN106371139A CN106371139A (zh) 2017-02-01
CN106371139B true CN106371139B (zh) 2019-04-26

Family

ID=57878853

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610676257.8A Active CN106371139B (zh) 2016-08-16 2016-08-16 一种应用层段振幅属性识别岩性油藏的方法及装置

Country Status (1)

Country Link
CN (1) CN106371139B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107576984B (zh) * 2017-07-27 2019-05-07 中国石油天然气集团公司 地震振幅属性提取过程中时窗的确定方法和装置
CN112392469A (zh) * 2019-08-12 2021-02-23 中国石油天然气集团有限公司 储层特征分类确定方法及装置
CN110703354B (zh) * 2019-10-23 2021-11-16 中国石油化工股份有限公司 一种基于井震结合的岩性边界识别方法
CN113514883B (zh) * 2021-06-18 2023-03-17 中国石油化工股份有限公司 一种断层-岩性油藏刻画方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102147478A (zh) * 2010-12-29 2011-08-10 中国海洋大学 复杂油藏的叠前低频信号识别方法
CN102305943A (zh) * 2011-07-22 2012-01-04 中国石油天然气股份有限公司 基于地震子波衰减谱的油气检测方法及装置
CN103149589A (zh) * 2013-02-22 2013-06-12 中国石油天然气股份有限公司 火成岩油气勘探方法及装置

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8185316B2 (en) * 2007-05-25 2012-05-22 Prime Geoscience Corporation Time-space varying spectra for seismic processing

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102147478A (zh) * 2010-12-29 2011-08-10 中国海洋大学 复杂油藏的叠前低频信号识别方法
CN102305943A (zh) * 2011-07-22 2012-01-04 中国石油天然气股份有限公司 基于地震子波衰减谱的油气检测方法及装置
CN103149589A (zh) * 2013-02-22 2013-06-12 中国石油天然气股份有限公司 火成岩油气勘探方法及装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
四扣洼陷沙四段灰岩地震识别与预测;朱定蓉,等;《油气地球物理》;20090430;第7卷(第2期);34-38

Also Published As

Publication number Publication date
CN106371139A (zh) 2017-02-01

Similar Documents

Publication Publication Date Title
CN106526670B (zh) 一种碎屑岩储层中地震属性砂体空间分布描述及评价的方法
CN101158724B (zh) 基于偶极小波的储层厚度预测方法
CN102109611B (zh) 一种利用低频地震属性预测油藏优质储层的方法
CN106371139B (zh) 一种应用层段振幅属性识别岩性油藏的方法及装置
US20030200030A1 (en) Method of processing acquired seismic data
CN105445800A (zh) 一种厚层砂体顶部分异岩性油藏的识别方法
CN109725348A (zh) 一种基于地震资料识别沉积相的方法
Naseer et al. Characterization of shallow-marine reservoirs of Lower Eocene carbonates, Pakistan: Continuous wavelet transforms-based spectral decomposition
CN105116449B (zh) 一种弱反射储层的识别方法
CN109541685A (zh) 一种河道砂体识别方法
CN105842733A (zh) 页岩储层地震识别方法
Karstens et al. Deep-seated focused fluid migration as indicator for hydrocarbon leads in the East Shetland Platform, North Sea Province
CN114910964B (zh) 一种断陷湖盆陡坡带砂砾岩体甜点区的预测方法
CN104965225A (zh) 一种基于子波旁瓣信息的油藏识别方法
CN102721979A (zh) 一种基于地震资料的薄层自动解释及厚度预测方法和装置
CN107817519A (zh) 基于重磁震综合解释技术的基底岩性识别方法及系统
CN106199725A (zh) 一种基于正振幅求和属性的煤岩厚度预测方法及装置
CN112505754B (zh) 基于高精度层序格架模型的井震协同划分沉积微相的方法
CN114488305B (zh) 一种无井新探区地震资料地质层位精细标定方法
Lyons et al. Seismic Exploration for Stratigraphic Traps: Geophysical Exploration Methods
Jarvis et al. Aquifer heterogeneity from SH-wave seismic impedance inversion
CN110795513B (zh) 河流相源储异位型致密油气甜点区分布的预测方法
RU2145101C1 (ru) Способ оценки эксплуатационных свойств нефтегазовой залежи
RU2263935C2 (ru) Способ поисков месторождений и залежей нефти и газа
Martinez et al. Utility of ground-penetrating radar in near-surface, high-resolution imaging of Lansing-Kansas City (Pennsylvanian) limestone reservoir analogs

Legal Events

Date Code Title Description
C06 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