CN109594971B - 基于核磁共振测井增强扩散气层识别因子的流体性质识别方法 - Google Patents
基于核磁共振测井增强扩散气层识别因子的流体性质识别方法 Download PDFInfo
- Publication number
- CN109594971B CN109594971B CN201811574908.8A CN201811574908A CN109594971B CN 109594971 B CN109594971 B CN 109594971B CN 201811574908 A CN201811574908 A CN 201811574908A CN 109594971 B CN109594971 B CN 109594971B
- Authority
- CN
- China
- Prior art keywords
- gas layer
- magnetic resonance
- nuclear magnetic
- enhanced diffusion
- diffusion gas
- 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
Links
- 238000005481 NMR spectroscopy Methods 0.000 title claims abstract description 72
- 239000012530 fluid Substances 0.000 title claims abstract description 64
- 238000009792 diffusion process Methods 0.000 title claims abstract description 50
- 238000000034 method Methods 0.000 title claims abstract description 31
- 238000001228 spectrum Methods 0.000 claims abstract description 44
- 230000003595 spectral effect Effects 0.000 claims description 7
- 238000012360 testing method Methods 0.000 claims description 5
- 238000000354 decomposition reaction Methods 0.000 claims description 3
- 238000004611 spectroscopical analysis Methods 0.000 claims description 3
- 239000007789 gas Substances 0.000 description 40
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 8
- 239000003345 natural gas Substances 0.000 description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 4
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 239000001257 hydrogen Substances 0.000 description 2
- 229910052739 hydrogen Inorganic materials 0.000 description 2
- 239000011148 porous material Substances 0.000 description 2
- 238000000685 Carr-Purcell-Meiboom-Gill pulse sequence Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 238000002592 echocardiography Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 125000004435 hydrogen atom Chemical group [H]* 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000000655 nuclear magnetic resonance spectrum Methods 0.000 description 1
- 230000005311 nuclear magnetism Effects 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 230000010287 polarization Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK 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 OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明公开了一种基于核磁共振测井增强扩散气层识别因子的流体性质识别方法，包括以下步骤：1)采用双TE采集模式进行核磁共振测井，获取核磁共振测井数据；2)对步骤1)获取的核磁共振测井数据进行回波反演，得到长、短回波间隔的T2谱；3)提取步骤2)得到的长、短回波间隔T2谱的特征参数；4)根据提取的长、短回波间隔T2谱的特征参数构建增强扩散气层识别因子；5)根据步骤4)得到的增强扩散气层识别因子识别储层流体性质，完成基于核磁共振测井增强扩散气层识别因子的流体性质识别，该方法能够实现流体性质的定量识别。
Description
技术领域
本发明属于天然气储层评价技术领域，涉及一种基于核磁共振测井增强扩散气层识别因子的流体性质识别方法。
背景技术
非常规类型的低渗透致密天然气储层作为一种新的气藏类型已被人们所重视，这类气藏的主要特征是其储层岩性致密、物性差、孔隙度和渗透率都很低、孔隙结构复杂，一些现有的常规测井解释方法不能正确评价这类低渗透致密天然气储层。
核磁共振测井能够直接测量地层中流体氢原子信息，不受岩石骨架的影响，专门设计的采集模式可以帮助识别油气水层。核磁共振测井流体识别方法在主要包括差谱法和移谱法，差谱法是利用不同流体的极化率不同采用不同等待时间的观测模式来识别流体性质；移谱法是利用不同流体的扩散特性不同采用不同回波间隔的采集模式来识别流体性质。然而在低渗透致密天然气储层中，受孔隙空间小、含氢指数低等因素的影响，核磁采集信号较微弱，差谱法应用效果不理想；而移谱法主要依靠经验定性识别，无法定量化识别流体性质。
发明内容
本发明的目的在于克服上述现有技术的缺点，提供了一种基于核磁共振测井增强扩散气层识别因子的流体性质识别方法，该方法能够实现流体性质的定量识别。
为达到上述目的，本发明所述的基于核磁共振测井增强扩散气层识别因子的流体性质识别方法包括以下步骤：
1)采用双TE采集模式进行核磁共振测井，获取核磁共振测井数据；
2)对步骤1)获取的核磁共振测井数据进行回波反演，得到长、短回波间隔的T2谱；
3)提取步骤2)得到的长、短回波间隔T2谱的特征参数；
4)根据提取的长、短回波间隔T2谱的特征参数构建增强扩散气层识别因子；
5)根据步骤4)得到的增强扩散气层识别因子识别储层流体性质，完成基于核磁共振测井增强扩散气层识别因子的流体性质识别。
步骤1)中双TE采集模式为单TW/双TE采集模式或者双TW/双TE采集模式。
步骤2)的具体操作为：
所述核磁共振测井数据为由一系列符合指数衰减规律衰减曲线组成的回波串信号，需要从观测到的回波串信号中分解出各弛豫分量T2i对应的幅度φi，核磁共振测井中获得的回波串信号描述为以下形式：
其中，y(t)为在t时刻采样所得的回波串孔隙度，T2min为T2谱弛豫时间的初始值，T2max为T2谱弛豫时间的终止值，式(1)为第一类Fred Holm积分方程，可以通过奇异值分解法来求解式(1)，得到长、短回波间隔的T2谱。
步骤3)中长、短回波间隔T2谱的特征参数包括可动流体T2几何均值、T2谱主峰位置及T2谱可动峰宽度。
可动流体T2几何均值T2fgm的表达式为：
其中，φi为对应分量T2i的孔隙度分量，T2c为核磁共振T2截止值，T2nonzero为φi非零时对应的最大T2谱弛豫时间。
核磁共振T2截止值T2c的表达式为：
T2c＝aT2gm b (3)
其中，T2gm为核磁共振T2谱的几何平均值，a及b为模型参数，a及b通过核磁共振实验数据拟合获取。
T2谱主峰位置T2mian为孔隙度分量φi最大时对应的T2i。
T2谱可动峰宽度T2width的表达式为：
T2width＝log10(T2nonzero)-log10(T2c) (4)
其中，T2c为核磁共振T2截止值，T2nonzero为φi非零时对应的最大T2谱弛豫时间。
增强扩散气层识别因子FLAG的表达式为：
其中，T2fgmA、T2mianA及T2widthA分别表示短回波间隔A组T2谱的可动流体几何均值、主峰位置点及可动峰宽度，T2fgmD、T2mianD及T2widthD分别表示长回波间隔D组T2谱的可动流体几何均值、主峰位置点及可动峰宽度。
步骤5)的具体操作为：根据区域核磁共振测井增强扩散气层识别因子的计算结果与试气资料进行综合分析，建立区域增强扩散气层识别因子的流体性质识别标准，再根据建立的区域增强扩散气层识别因子的流体性质识别标准对储层进行流体性质识别。
本发明具有以下有益效果：
本发明所述的基于核磁共振测井增强扩散气层识别因子的流体性质识别方法在具体操作时，基于核磁共振传统移谱法的原理，深入挖掘核磁共振测井丰富信息，利用气、水扩散系数的差异引起核磁共振测井响应不同来识别流体性质，具体的，对核磁共振测井数据进行回波反演，得到长、短回波间隔的T2谱，再根据长、短回波间隔的T2谱的特征参数构建增强扩散气层识别因子，最后根据增强扩散气层识别因子识别储层流体性质，实现核磁共振测井的流体性质定量识别，克服传统核磁共振经典移谱法无法定量识别流体性质的缺陷，提高低渗透致密气层流体性质识别准确率。
附图说明
图1为本发明的流程图；
图2为本发明中核磁共振T2的谱特征参数示意图；
图3为核磁T2截止值与核磁T2几何均值的关系图；
图4为基于核磁共振测井增强扩散气层识别因子的流体性质识别实例图。
具体实施方式
下面结合附图对本发明做进一步详细描述：
参考图1，本发明所述的基于核磁共振测井增强扩散气层识别因子的流体性质识别方法包括以下步骤：
1)采用双TE采集模式进行核磁共振测井，获取核磁共振测井数据；
核磁共振测井是按照中华人民共和国石油天然气行业标准《核磁共振成像测井作业技术规范SY/T 67752-2009》规定的流程进行的，步骤1)中双TE采集模式为单TW/双TE采集模式或者双TW/双TE采集模式。
2)对步骤1)获取的核磁共振测井数据进行回波反演，得到长、短回波间隔的T2谱；
步骤2)的具体操作为：
核磁共振测井以地层中流体氢核为研究对象，通过重复使用CPMG脉冲得弛豫信号，所述核磁共振测井数据为由一系列符合指数衰减规律衰减曲线组成的回波串信号，但是在这个回波串信号中无法辨别出每一种衰减所对应的具体信息，需要从观测到的回波串信号中分解出各弛豫分量T2i对应的幅度φi，核磁共振测井中获得的回波串信号描述为以下形式：
其中，y(t)为在t时刻采样所得的回波串孔隙度，T2min为T2谱弛豫时间的初始值，T2max为T2谱弛豫时间的终止值，式(1)为第一类Fred Holm积分方程，可以通过奇异值分解法来求解式(1)，得到长、短回波间隔的T2谱。
3)提取步骤2)得到的长、短回波间隔T2谱的特征参数；
步骤3)中长、短回波间隔T2谱的特征参数包括可动流体T2几何均值、T2谱主峰位置及T2谱可动峰宽度，参考图2。
可动流体T2几何均值T2fgm的表达式为：
其中，φi为对应分量T2i的孔隙度分量，T2c为核磁共振T2截止值，T2nonzero为φi非零时对应的最大T2谱弛豫时间。
核磁共振T2截止值T2c的表达式为：
T2c＝aT2gm b (3)
其中，T2gm为核磁共振T2谱的几何平均值，a及b为模型参数，a及b通过核磁共振实验数据拟合获取，将核磁共振实验测量数据得到的核磁T2截止值与T2几何均值进行拟合，得到a＝1.415，b＝0.647，如图3所示。
T2谱主峰位置T2mian为孔隙度分量φi最大时对应的T2i。
T2谱可动峰宽度T2width的表达式为：
T2width＝log10(T2nonzero)-log10(T2c) (4)
其中，T2c为核磁共振T2截止值，T2nonzero为φi非零时对应的最大T2谱弛豫时间。
4)根据提取的长、短回波间隔T2谱的特征参数构建增强扩散气层识别因子；
增强扩散气层识别因子FLAG的表达式为：
其中，T2fgmA、T2mianA及T2widthA分别表示短回波间隔A组T2谱的可动流体几何均值、主峰位置点及可动峰宽度，T2fgmD、T2mianD及T2widthD分别表示长回波间隔D组T2谱的可动流体几何均值、主峰位置点及可动峰宽度。
5)根据步骤4)得到的增强扩散气层识别因子识别储层流体性质，完成基于核磁共振测井增强扩散气层识别因子的流体性质识别。
步骤5)的具体操作为：根据区域核磁共振测井增强扩散气层识别因子的计算结果与试气资料进行综合分析，建立区域增强扩散气层识别因子的流体性质识别标准(如表1所示)，再根据建立的区域增强扩散气层识别因子的流体性质识别标准对储层进行流体性质识别。
表1
| 储层流体性质 | 增强扩散气层识别因子 |
| 气层 | D≥0.1 |
| 差气层或气水同层 | 0.01≤D＜0.1 |
| 含气水层或气层 | D＜0.01 |
实际资料处理过程中，核磁共振测井增强扩散气层识别因子是通过编写程序实现。图4为基于核磁共振测井增强扩散气层识别因子的流体性质识别应用实例，第一道为岩性曲线，第二道为阵列侧向曲线，第三道为孔隙度曲线，第四道为深度道，第五道为深电阻率曲线，第六道为A组回波反演T2谱，第七道为D组回波反演T2谱，第八道曲线为气层识别因子FLAG。可以看出：12号层和14号层气层识别因子大于0.1，综合解释为气层，该井12层试气日产气22.49×104m3/d，产水1.2m3/d，试气为气层，验证了本发明的有效性。
上述本发明实施例序号仅仅为了描述，不代表实施例的优劣。
以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (10)
1.一种基于核磁共振测井增强扩散气层识别因子的流体性质识别方法，其特征在于，包括以下步骤：
1)采用双TE采集模式进行核磁共振测井，获取核磁共振测井数据；
2)对步骤1)获取的核磁共振测井数据进行回波反演，得到长、短回波间隔的T2谱；
3)提取步骤2)得到的长、短回波间隔T2谱的特征参数；
4)根据提取的长、短回波间隔T2谱的特征参数构建增强扩散气层识别因子；
5)根据步骤4)得到的增强扩散气层识别因子识别储层流体性质，完成基于核磁共振测井增强扩散气层识别因子的流体性质识别。
2.根据权利要求1所述的一种基于核磁共振测井增强扩散气层识别因子的流体性质识别方法，其特征在于，步骤1)中双TE采集模式为单TW/双TE采集模式或者双TW/双TE采集模式。
3.根据权利要求1所述的基于核磁共振测井增强扩散气层识别因子的流体性质识别方法，其特征在于，步骤2)的具体操作为：
所述核磁共振测井数据为由一系列符合指数衰减规律衰减曲线组成的回波串信号，需要从观测到的回波串信号中分解出各弛豫分量T2i对应的孔隙度分量φi，核磁共振测井中获得的回波串信号描述为以下形式：
其中，y(t)为在t时刻采样所得的回波串孔隙度，T2min为T2谱弛豫时间的初始值，T2max为T2谱弛豫时间的终止值，式(1)为第一类Fred Holm积分方程，可以通过奇异值分解法来求解式(1)，得到长、短回波间隔的T2谱。
4.根据权利要求1所述的基于核磁共振测井增强扩散气层识别因子的流体性质识别方法，其特征在于，步骤3)中长、短回波间隔T2谱的特征参数包括可动流体T2几何均值、T2谱主峰位置及T2谱可动峰宽度。
5.根据权利要求4所述的基于核磁共振测井增强扩散气层识别因子的流体性质识别方法，其特征在于，可动流体T2几何均值T2fgm的表达式为：
其中，φi为对应分量T2i的孔隙度分量，T2c为核磁共振T2截止值，T2nonzero为φi非零时对应的最大T2谱弛豫时间。
6.根据权利要求5所述的基于核磁共振测井增强扩散气层识别因子的流体性质识别方法，其特征在于，核磁共振T2截止值T2c的表达式为：
T2c＝aT2gm b (3)
其中，T2gm为核磁共振T2谱的几何平均值，a及b为模型参数，a及b通过核磁共振实验数据拟合获取。
7.根据权利要求4所述的基于核磁共振测井增强扩散气层识别因子的流体性质识别方法，其特征在于，T2谱主峰位置T2mian为孔隙度分量φi最大时对应的T2i。
8.根据权利要求4所述的基于核磁共振测井增强扩散气层识别因子的流体性质识别方法，其特征在于，T2谱可动峰宽度T2width的表达式为：
T2width＝log10(T2nonzero)-log10(T2c) (4)
其中，T2c为核磁共振T2截止值，T2nonzero为φi非零时对应的最大T2谱弛豫时间。
9.根据权利要求1所述的基于核磁共振测井增强扩散气层识别因子的流体性质识别方法，其特征在于，增强扩散气层识别因子FLAG的表达式为：
其中，T2fgmA、T2mianA及T2widthA分别表示短回波间隔A组T2谱的可动流体几何均值、主峰位置点及可动峰宽度，T2fgmD、T2mianD及T2widthD分别表示长回波间隔D组T2谱的可动流体几何均值、主峰位置点及可动峰宽度。
10.根据权利要求1所述的基于核磁共振测井增强扩散气层识别因子的流体性质识别方法，其特征在于，步骤5)的具体操作为：根据区域核磁共振测井增强扩散气层识别因子的计算结果与试气资料进行综合分析，建立区域增强扩散气层识别因子的流体性质识别标准，再根据建立的区域增强扩散气层识别因子的流体性质识别标准对储层进行流体性质识别。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811574908.8A CN109594971B (zh) | 2018-12-21 | 2018-12-21 | 基于核磁共振测井增强扩散气层识别因子的流体性质识别方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811574908.8A CN109594971B (zh) | 2018-12-21 | 2018-12-21 | 基于核磁共振测井增强扩散气层识别因子的流体性质识别方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109594971A CN109594971A (zh) | 2019-04-09 |
| CN109594971B true CN109594971B (zh) | 2022-07-05 |
Family
ID=65963358
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811574908.8A Active CN109594971B (zh) | 2018-12-21 | 2018-12-21 | 基于核磁共振测井增强扩散气层识别因子的流体性质识别方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109594971B (zh) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110244369B (zh) * | 2019-06-28 | 2020-10-13 | 中国石油大学(北京) | 储层束缚和可动流体分布确定方法、装置及系统 |
| CN110344825A (zh) * | 2019-06-28 | 2019-10-18 | 中国石油天然气股份有限公司 | 一种低孔低渗低电阻率砂岩油层综合判识方法 |
| CN111042810B (zh) * | 2019-12-27 | 2021-06-22 | 中国石油天然气集团有限公司 | 基于核磁量化表征的砂砾岩储层产能下限确定方法及系统 |
| CN111538099B (zh) * | 2020-05-06 | 2023-05-26 | 中国石油天然气集团有限公司 | 一种核磁共振测井速度确定方法 |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102141637B (zh) * | 2010-01-28 | 2012-10-17 | 中国石油天然气股份有限公司 | 一种利用核磁共振测井资料连续定量评价储集层孔隙结构的方法 |
| CN104215652B (zh) * | 2014-08-21 | 2016-08-31 | 中国石油天然气股份有限公司 | 确定油气饱和度的方法和装置 |
| US10739489B2 (en) * | 2016-01-15 | 2020-08-11 | Baker Hughes, A Ge Company, Llc | Low gradient magnetic resonance logging for measurement of light hydrocarbon reservoirs |
| CN108049866B (zh) * | 2017-12-11 | 2021-05-14 | 中石化石油工程技术服务有限公司 | 二维核磁共振测井致密气藏定量评价方法 |
-
2018
- 2018-12-21 CN CN201811574908.8A patent/CN109594971B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN109594971A (zh) | 2019-04-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109594971B (zh) | 基于核磁共振测井增强扩散气层识别因子的流体性质识别方法 | |
| CN101943669B (zh) | 低场核磁共振测量钻井液含油率的方法 | |
| CN109753755B (zh) | 一种确定储层含水饱和度的方法 | |
| CN110231272B (zh) | 致密砂岩孔径与核磁共振t2值转换关系的确定方法及系统 | |
| CN108049866B (zh) | 二维核磁共振测井致密气藏定量评价方法 | |
| CN109725016B (zh) | 一种用于含重油、沥青质岩心的核磁共振实验室测量方法 | |
| CN103528934A (zh) | 一种测量超低渗岩石渗透率应力敏感性的互相关技术 | |
| CN109932297B (zh) | 一种致密砂岩储层渗透率的计算方法 | |
| CN103628871A (zh) | 一种基于阿尔奇公式的电阻率侵入校正的新方法 | |
| CN109856688B (zh) | 基于核磁测井双tw极化增强法的流体性质识别方法 | |
| CN113503156B (zh) | 储层垂向非均质性定量测定方法、系统、介质、设备及终端 | |
| CN103630939A (zh) | 一种气层识别评价方法 | |
| CN106290103B (zh) | 一种页岩气储层中粘土微孔孔隙度的测定方法 | |
| CN110700820A (zh) | 松辽盆地北部致密油储层甜点分类方法 | |
| CN106897531A (zh) | 一种低渗透石灰岩储层渗透率的定量评价方法 | |
| CN110905493B (zh) | 一种测量井下地层流体污染率的方法 | |
| RU2009143585A (ru) | Способ разработки неоднородного массивного или многопластового газонефтяного или нефтегазоконденсатного месторождения | |
| CN110017136B (zh) | 一种基于视水层电阻率的水淹层识别与产水率预测方法 | |
| CN110057853B (zh) | 一种基于低场核磁共振响应的岩石杨氏模量计算方法 | |
| CN105590018A (zh) | 砂泥岩薄互层油藏的油水层判别方法 | |
| CN109655394B (zh) | 一种孔喉结构参数约束下的核磁共振t2谱渗透率计算方法 | |
| CN104122182A (zh) | 矿井储层有效厚度下限的获取方法 | |
| CN105089632A (zh) | 一种高温高压储层co2流体纵波时差骨架参数的获取方法 | |
| CN112012727B (zh) | 获得气相有效渗透率的方法及储层产能的预测方法 | |
| CN112031753B (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 |