CN110485993A - 一种地震资料套变点预测方法 - Google Patents
一种地震资料套变点预测方法 Download PDFInfo
- Publication number
- CN110485993A CN110485993A CN201910721467.8A CN201910721467A CN110485993A CN 110485993 A CN110485993 A CN 110485993A CN 201910721467 A CN201910721467 A CN 201910721467A CN 110485993 A CN110485993 A CN 110485993A
- Authority
- CN
- China
- Prior art keywords
- fold
- little
- curvature
- amount
- seismic
- 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
Links
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
- E21B47/10—Locating fluid leaks, intrusions or movements
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Mining
Abstract
本发明公开了一种地震资料套变点预测方法,涉及油气田勘探开发的地震资料应用技术领域,a、高分辨率处理;b、精细层位对比;c、褶皱自动识别:沿地震层位求取地层曲率,发现褶皱,同时判断曲率正负;d、褶皱性质判断:根据曲率正负变化判断是背斜还是向斜;e、提取振幅差异:以微幅褶皱中心计算目的层范围内各个方向振幅差异;f、计算微幅褶皱量:根据曲率正负符号、曲率大小和振幅差异值三个参数,综合计算微幅褶皱量;g、套变点预测:根据井轨迹坐标,从微幅褶皱量数据中拾取可能的套变点位置及其微幅褶皱程度,沿水平井井轨迹提取微幅褶皱量,即得到套变点预测结果。本方法根据微幅褶皱来预测水平井的套变点的准确性和效率显著提高。
Description
技术领域
本发明涉及油气田勘探开发的地震资料应用技术领域,确切地说涉及一种地震资料套变点预测方法。
背景技术
四川盆地页岩气水平井套管作业过程中,频繁发生套变,即下到井中的套管发生变形破坏不能作业,轻则影响生产时效,重则造成井下作业报废,严重影响了页岩气勘探开发工程进展和开发效果,造成巨大损失和浪费。套变发生的原因可能是以下的一种或几种:地质情况突变、井壁不稳、套管材料不满足地层变化条件、工艺措施不适应等因素。
导致套变发生的根本地质原因在于施加在套管上的应力不均衡使得套管破损,直接地质原因是地质因素发生了突变,就可能导致套变,如地层形态扭曲、岩性变化、裂缝出现等,一种间接地质原因则是在很大程度上由于没有预先观察、判断地质因素的变化,井下作业时仍按照正常情况进行,没有及时采取应变措施,导致套变发生。其中微幅褶皱导致的套变情况在很大程度上是由于没有事先识别出微幅褶皱。其主要原因在于地震资料的精度不够,难以识别出微幅褶皱。
地震资料精度不够主要是指地震资料的分辨率低。通常的地震资料(纵向)分辨率低主要表现在主频低(30Hz左右)、频带宽度窄(10-70Hz),分辨的地层厚度大(15m以上)。这样的地震资料撇开断层外,地震反射同相轴一般具有粗、光滑、连续性较强的特点,难以观察到微幅褶皱。由于微幅褶皱观察不到,设计的井轨迹就会直接穿越褶皱,套管作业时发生应力不均,造成变形。为了准确识别微幅褶皱,必须显著提高分辨率。
提高分辨率的方法很多,但一般方法仅能有限提高分辨率,如提高30%左右,不能根本解决微幅褶皱的识别问题。本发明人的发明专利(ZL 2013 1 0333410.3)实现了大幅度提高分辨率的效果,可提高100%以上,为微幅褶皱的识别奠定了基础,形成了基于高分辨率地震资料的微幅褶皱识别方法。根据这样的方法,可以有效计算识别出微幅褶皱。
作者为朱洪昌、朱莉、玄长虹、刘升余、余学兵等在刊名为《石油地球物理勘探》的期刊上发表了题名为“运用高分辨率地震资料处理技术识别薄储层及微幅构造”的期刊文献,发表日期为2010年11月第45卷增刊1。该文献主要公开了:为了识别薄储层、微幅构造、小断层等小尺度地质体, 本文研究了以保护和提高地震资料高频成分信噪比为基础的提高分辨率处理技术, 分别在叠前和叠后有效压缩子波、拓宽地震资料的频带, 提高了地震资料的分辨率。在胜利探区永新工区地震资料处理中的应用结果表明, 成果剖面的弱反射层次清楚, 微幅构造成像清晰,地质信息丰富, 勘探目的层的主频由原来的20~ 23H z提高到46 ~ 50Hz , 对提高地震数据的处理质量, 解决陆相薄互层油气藏的微幅构造成像问题非常有效。
由此可见,现有技术已经有利用提高地震资料分辨率的方法来识别微幅构造的文献报道。同时,现有技术中,褶皱识别还可以通过计算曲率方式来获得,这是一种通常的做法,但是,以上述专利和期刊文献为代表的现有技术。并未涉及到根据识别的微幅褶皱,对水平井套变点进行预测。同时,由于上述现有技术中对微幅构造的识别仅采用单一的计算曲率的方式获得,或者单一的采用提高地震分辨率来识别微幅构造。在实际应用过程中,即使用上述现有技术来对水平井套变点进行预测,由于微幅构造的识别较为单一,因此识别效果不理想,因而导致对套变点的预测准确性和效率也不理想。
发明内容
本发明旨在针对上述现有技术所存在的缺陷和不足,提供一种地震资料套变点预测方法,本方法通过利用微幅褶皱来预测水平井的套变点,微幅褶皱识别准确性得到显著提高,识别效率较高,由此导致根据微幅褶皱来预测水平井的套变点的准确性和效率显著提高。
本发明是通过采用下述技术方案实现的:
一种地震资料套变点预测方法,其特征在于:
a、高分辨率处理:将常规的较低分辨率的地震资料处理成高分辨率资料;
b、精细层位对比:按照通常的地震层位对比技术,精确拾取地震层位;
c、褶皱自动识别:沿地震层位求取地层曲率,发现褶皱,同时判断曲率正负;
d、褶皱性质判断:根据曲率正负变化判断是背斜还是向斜;
e、提取振幅差异:以微幅褶皱中心计算目的层范围内各个方向振幅差异;
f、计算微幅褶皱量:根据曲率正负符号、曲率大小和振幅差异值三个参数,综合计算微幅褶皱量;
g、套变点预测:根据井轨迹坐标,从微幅褶皱量数据中拾取可能的套变点位置及其微幅褶皱程度,沿水平井井轨迹提取微幅褶皱量,即得到套变点预测结果。
所述a步骤中,若地震资料分辨率已经较高,可忽略此步,高分辨率处理方法是任意的,但分辨率应当有显著的提高。
所述b步骤中,拾取的地震层位应当处于同一个反射特征点或地质界面上,所述的反射特征点为波峰或波谷。
所述d步骤中,是根据曲率正负得到褶皱方向,褶皱方向为上凸则代表是背斜,褶皱方向是下凹则代表是向斜或挠曲,向斜对应的褶皱性质为压性应力,背斜对应的褶皱性质为张性应力。
所述f步骤中,微幅褶皱量=Σ曲率正负符号×曲率大小×振幅差异值,Σ表示各个方向的微幅褶皱量相加,得到褶皱位置的微幅褶皱量,再对地层内各点进行计算,得到地层内各点的微幅褶皱量。
与现有技术相比,本发明所达到的有益效果如下:
一、本发明能够利用高分辨地震资料自动将微幅褶皱识别出来。与期刊文献“运用高分辨率地震资料处理技术识别薄储层及微幅构造”等现有技术相比,本发明识别的微幅褶皱包含了地层曲率信息、褶皱性质和振幅变化信息,明显有别于单一的曲率计算,也有别于其它单一参数的裂缝预测方法,提出了一种新的微幅褶皱识别方法,对储层的识别具有重要的指示意义,在实际运用过程中,其识别准确性得到有效提高,识别的效率提高。由此,在微幅褶皱被准确识别的基础上,根据井轨迹坐标,从微幅褶皱量数据中拾取可能的套变点位置及其微幅褶皱程度。这样的处理方式,导致根据微幅褶皱来预测水平井的套变点的准确性和效率显著提高。从一个重要方面解决套变点预测难题,减少套变带来的经济损失。
二,本发明中,曲率反映了褶皱强度,曲率正负反映了褶皱性质,间接反映了应力性质(张性或压性);振幅变化反映了储层的变化。这些都可反映褶皱的储层特性。本发明提出的微幅褶皱量=Σ曲率正负符号×曲率大小×振幅差异值这一计算公式,可更加准确、有效地判断储层综合特性。该方法根据油气勘探开发实际提出,具有科学性、创新性,可用于储层预测,指导水平井轨迹设计和地质导向,也可用于套变、井漏等工程异常预测,减小钻井及井下作业的工程复杂。
附图说明
下面将结合说明书附图和具体实施方式对本发明作进一步的详细说明,其中:
图1为本发明的工艺流程图。
图2为微幅褶皱计算原理中纵向计算剖面图;
图3为微幅褶皱计算原理中横向计算俯视图;
图4为四川盆地某页岩气地震剖面(上)及对应的微幅褶皱识别剖面(下)图——严重不清晰,不符合专利法要求,要求不能有彩色,只能是黑白的,同时,字体不清晰,请尽量补充,修改为字体清晰的。
图5为四川盆地某页岩气层微幅褶皱平面分布图——严重不清晰,不符合专利法要求,要求不能有彩色,只能是黑白的,同时,字体不清晰。请尽量补充,修改为字体清晰的。
具体实施方式
实施例1
作为本发明最佳实施方式,其公开了一种地震资料井漏预测方法,其步骤为:
a、高分辨率处理:将常规的较低分辨率的地震资料处理成高分辨率资料;
b、精细层位对比:按照通常的地震层位对比技术,精确拾取地震层位;
c、褶皱自动识别:沿地震层位求取地层曲率,发现褶皱,同时判断曲率正负;
d、褶皱性质判断:根据曲率正负变化判断是背斜还是向斜;
e、提取振幅差异:以微幅褶皱中心计算目的层范围内各个方向振幅差异;
f、计算微幅褶皱量:根据曲率正负符号、曲率大小和振幅差异值三个参数,综合计算微幅褶皱量;
g、套变点预测:根据井轨迹坐标,从微幅褶皱量数据中拾取可能的套变点位置及其微幅褶皱程度,沿水平井井轨迹提取微幅褶皱量,即得到套变点预测结果。
上述第a步提到的高分辨率是相对的,根据实际预测效果,若地震资料分辨率已经较高,不需先作高分辨率处理,否则,要先作高分辨率处理,提高分辨率的方法是任意的,任何高分辨率处理方法都包含在内;第f步的微幅褶皱量计算公式不是固定的,但都包含曲率大小、曲率正负符号及振幅差异值。第g步的水平井不仅限于页岩气,包含了各种储层类型的水平井。
术语解释
褶皱:是指组成地壳的岩层受构造应力的作用发生波状弯曲而未失去连续性的构造,一般包含背斜、向斜、挠曲等构造。
微幅褶皱:亦即微褶皱,是一个相对概念,是指幅度较小的褶皱,其在纵向或横向上的规模相对较小。微幅褶皱在分辨率较低的地震剖面上不易识别。本发明指的微幅褶皱即是指较低分辨率地震资料上难以识别,而在相对较高分辨率地震资料能够识别的褶皱,包括向斜、背斜以及挠曲。
实施例2
具体应用实例:参照2和3,可以明确具体计算方法原理。
图2中包含一个背斜和一个向斜,仅以背斜为例。确定地层中心,地层厚度(在地震剖面为时间域时为时窗长度)为地层顶、底界之间的范围。根据地层中心拾取到的层位数据,可计算曲率,发现褶皱,并得到曲率的正负(背、向斜性质),确定出褶皱核部及褶皱中心所在地震道——中心道,然后从上到下逐点计算中心道两侧振幅差。实际计算时,还根据图3计算中心道各个方向多道间的振幅差,然后按照公式:Σ曲率×振幅差,即得到该中心道微幅褶皱量,再对地层内各点进行计算,得到地层内各点的微幅褶皱量。对每个褶皱重复上述步骤,即可得到井轨迹上的所有微幅褶皱识别最终结果。
图4为某地区一段地震剖面实施本方法得到微幅褶皱的实例。图上部为原始地震剖面,下部为页岩气层微幅褶皱检测剖面。中间的黑色线为页岩气层层位线,层位线上下的彩色为微幅褶皱综合值——含振幅差异、曲率大小和方向。该图可清晰识别出各微幅褶皱及其差异,反映出褶皱强弱变化。图5为该地区页岩气层微幅褶皱识别平面分布及实际发生的所有套变点分布图。从图中可见套变点一般都分布在微幅褶皱上或与微幅褶皱密切相关。
Claims (5)
1.一种地震资料套变点预测方法,其特征在于:
a、高分辨率处理:将常规的较低分辨率的地震资料处理成高分辨率资料;
b、精细层位对比:按照通常的地震层位对比技术,精确拾取地震层位;
c、褶皱自动识别:沿地震层位求取地层曲率,发现褶皱,同时判断曲率正负;
d、褶皱性质判断:根据曲率正负变化判断是背斜还是向斜;
e、提取振幅差异:以微幅褶皱中心计算目的层范围内各个方向振幅差异;
f、计算微幅褶皱量:根据曲率正负符号、曲率大小和振幅差异值三个参数,综合计算微幅褶皱量;
g、套变点预测:根据井轨迹坐标,从微幅褶皱量数据中拾取可能的套变点位置及其微幅褶皱程度,沿水平井井轨迹提取微幅褶皱量,即得到套变点预测结果。
2.根据权利要求1所述的一种地震资料套变点预测方法,其特征在于:所述a步骤中,若地震资料分辨率已经较高,忽略此步,高分辨率处理方法是任意的。
3.根据权利要求1所述的一种地震资料套变点预测方法,其特征在于:所述b步骤中,拾取的地震层位应当处于同一个反射特征点或地质界面上,所述的反射特征点为波峰或波谷。
4.根据权利要求1所述的一种地震资料套变点预测方法,其特征在于:所述d步骤中,是根据曲率正负得到褶皱方向,褶皱方向为上凸则代表是背斜,褶皱方向是下凹则代表是向斜或挠曲,向斜对应的褶皱性质为压性应力,背斜对应的褶皱性质为张性应力。
5.根据权利要求1所述的一种地震资料套变点预测方法,其特征在于:所述f步骤中,微幅褶皱量=Σ曲率正负符号×曲率大小×振幅差异值,Σ表示各个方向的微幅褶皱量相加,得到褶皱位置的微幅褶皱量,再对地层内各点进行计算,得到地层内各点的微幅褶皱量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910721467.8A CN110485993B (zh) | 2019-08-06 | 2019-08-06 | 一种地震资料套变点预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910721467.8A CN110485993B (zh) | 2019-08-06 | 2019-08-06 | 一种地震资料套变点预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110485993A true CN110485993A (zh) | 2019-11-22 |
CN110485993B CN110485993B (zh) | 2022-05-06 |
Family
ID=68549907
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910721467.8A Active CN110485993B (zh) | 2019-08-06 | 2019-08-06 | 一种地震资料套变点预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110485993B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110456416A (zh) * | 2019-08-06 | 2019-11-15 | 中国石油集团川庆钻探工程有限公司 | 一种地震资料微幅褶皱识别方法 |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB1270590A (en) * | 1968-12-30 | 1972-04-12 | Chevron Res | Method of and system for attenuating multiple seismic signals in the determination of in-line and cross dips employing cross steered seismic data |
US20040064257A1 (en) * | 2002-09-27 | 2004-04-01 | Tobias Steven M. | Method of using underbalanced well data for seismic attribute analysis |
US20080033656A1 (en) * | 2005-01-07 | 2008-02-07 | Herwanger Jorg V | Processing a Seismic Monitor Survey |
CN101487898A (zh) * | 2009-02-27 | 2009-07-22 | 中国石油集团川庆钻探工程有限公司 | 采用纵波地震勘探叠后资料进行油气水识别的方法 |
CN102937720A (zh) * | 2011-08-15 | 2013-02-20 | 中国石油化工股份有限公司 | 井控提高地震资料分辨率的方法 |
CN103696780A (zh) * | 2014-01-07 | 2014-04-02 | 山东大学 | 安装于隧道掘进机刀盘中心的单孔地质雷达自动钻探装置 |
CN203835399U (zh) * | 2014-01-07 | 2014-09-17 | 山东大学 | 安装于隧道掘进机刀盘中心的单孔地质雷达自动钻探装置 |
CN106154326A (zh) * | 2016-08-18 | 2016-11-23 | 中国石油天然气集团公司 | 一种纵弯褶皱裂缝密度评价的方法及装置 |
CN107515430A (zh) * | 2017-06-06 | 2017-12-26 | 青岛大学 | 一种地震法探测盐湖卤水的方法 |
US20190080122A1 (en) * | 2017-09-14 | 2019-03-14 | Saudi Arabian Oil Company | Subsurface reservoir model with 3d natural fractures prediction |
CN109751037A (zh) * | 2017-11-01 | 2019-05-14 | 中国石油化工股份有限公司 | 一种常压页岩气藏高频变排量体积压裂方法 |
-
2019
- 2019-08-06 CN CN201910721467.8A patent/CN110485993B/zh active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB1270590A (en) * | 1968-12-30 | 1972-04-12 | Chevron Res | Method of and system for attenuating multiple seismic signals in the determination of in-line and cross dips employing cross steered seismic data |
US20040064257A1 (en) * | 2002-09-27 | 2004-04-01 | Tobias Steven M. | Method of using underbalanced well data for seismic attribute analysis |
US20080033656A1 (en) * | 2005-01-07 | 2008-02-07 | Herwanger Jorg V | Processing a Seismic Monitor Survey |
CN101487898A (zh) * | 2009-02-27 | 2009-07-22 | 中国石油集团川庆钻探工程有限公司 | 采用纵波地震勘探叠后资料进行油气水识别的方法 |
CN102937720A (zh) * | 2011-08-15 | 2013-02-20 | 中国石油化工股份有限公司 | 井控提高地震资料分辨率的方法 |
CN103696780A (zh) * | 2014-01-07 | 2014-04-02 | 山东大学 | 安装于隧道掘进机刀盘中心的单孔地质雷达自动钻探装置 |
CN203835399U (zh) * | 2014-01-07 | 2014-09-17 | 山东大学 | 安装于隧道掘进机刀盘中心的单孔地质雷达自动钻探装置 |
CN106154326A (zh) * | 2016-08-18 | 2016-11-23 | 中国石油天然气集团公司 | 一种纵弯褶皱裂缝密度评价的方法及装置 |
CN107515430A (zh) * | 2017-06-06 | 2017-12-26 | 青岛大学 | 一种地震法探测盐湖卤水的方法 |
US20190080122A1 (en) * | 2017-09-14 | 2019-03-14 | Saudi Arabian Oil Company | Subsurface reservoir model with 3d natural fractures prediction |
CN109751037A (zh) * | 2017-11-01 | 2019-05-14 | 中国石油化工股份有限公司 | 一种常压页岩气藏高频变排量体积压裂方法 |
Non-Patent Citations (2)
Title |
---|
梁顺军等: "地震解释经纬跟踪分析", 《石油地球物理勘探》 * |
武加鹤等: "四川盆地涪陵焦石坝地区五峰—龙马溪组低序级断层识别技术及应用效果", 《石油实验地质》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110456416A (zh) * | 2019-08-06 | 2019-11-15 | 中国石油集团川庆钻探工程有限公司 | 一种地震资料微幅褶皱识别方法 |
CN110456416B (zh) * | 2019-08-06 | 2021-05-18 | 中国石油集团川庆钻探工程有限公司 | 一种地震资料微幅褶皱识别方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110485993B (zh) | 2022-05-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101738639B (zh) | 提高岩石裂缝参数精度的方法 | |
Chopra et al. | Curvature attribute applications to 3D surface seismic data | |
Chopra et al. | Emerging and future trends in seismic attributes | |
US6498989B1 (en) | Method for predicting dynamic parameters of fluids in a subterranean reservoir | |
CN103454685B (zh) | 利用测井约束波阻抗反演预测砂体厚度的方法和装置 | |
CN103728659B (zh) | 一种提高地下岩溶探测精度的方法 | |
CA2823710C (en) | Methods and systems regarding models of underground formations | |
CN104678434B (zh) | 一种预测储层裂缝发育参数的方法 | |
CN104142519B (zh) | 一种泥岩裂缝油藏预测方法 | |
CN103777245B (zh) | 基于地震资料的油气成藏条件定量评价方法 | |
CN104636980A (zh) | 针对河道砂油藏类型油气汇集条件的地球物理表征方法 | |
WO2023000257A1 (zh) | 一种砂岩型铀矿成矿有利部位地质-地震三维预测方法 | |
Xia et al. | Application of 3D fine seismic interpretation technique in Dawangzhuang area, Bohai Bay Basin, Northeast China | |
CN107728205B (zh) | 一种地层压力预测方法 | |
CN106033127B (zh) | 基于横波速度变化率的地应力方位地震预测方法 | |
CN110485993A (zh) | 一种地震资料套变点预测方法 | |
Elebiju et al. | Investigation of links between Precambrian basement structure and Paleozoic strata in the Fort Worth Basin, Texas, USA, using high-resolution aeromagnetic (HRAM) data and seismic attributes | |
Mahlalela et al. | Structural characteristics and 3D seismic detection of gas migration pathways in the deep‐water Orange Basin, South Africa | |
CN104749618B (zh) | 泥页岩缓倾角裂缝叠后概率定量表征方法 | |
Harris et al. | Subseismic Fault Identification Using the Fault Likelihood Attribute: Application to Geosteering in the DJ Basin | |
Su et al. | Low relief structure interpretation and mapping of the Donghe Sandstone thin reservoir based on seismic data | |
CN110439534A (zh) | 一种地震资料井漏预测方法 | |
CN110456416A (zh) | 一种地震资料微幅褶皱识别方法 | |
CN105259576A (zh) | 一种利用地震统计特征的油气藏识别方法 | |
Zhang et al. | Identifying minor faults on top of coalfield Ordovician limestone stratum using seismic attributes derived from azimuthally stacked data |
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 |