CN109916943A - 岩心核磁共振t2谱图驰豫时间与压汞孔喉半径转换方法 - Google Patents
岩心核磁共振t2谱图驰豫时间与压汞孔喉半径转换方法 Download PDFInfo
- Publication number
- CN109916943A CN109916943A CN201910340453.1A CN201910340453A CN109916943A CN 109916943 A CN109916943 A CN 109916943A CN 201910340453 A CN201910340453 A CN 201910340453A CN 109916943 A CN109916943 A CN 109916943A
- Authority
- CN
- China
- Prior art keywords
- pore throat
- magnetic resonance
- throat radius
- nuclear magnetic
- rock core
- 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
- QSHDDOUJBYECFT-UHFFFAOYSA-N mercury Chemical compound [Hg] QSHDDOUJBYECFT-UHFFFAOYSA-N 0.000 title claims abstract description 80
- 229910052753 mercury Inorganic materials 0.000 title claims abstract description 80
- 239000011148 porous material Substances 0.000 title claims abstract description 70
- 238000005481 NMR spectroscopy Methods 0.000 title claims abstract description 49
- 239000011435 rock Substances 0.000 title claims abstract description 32
- 238000000034 method Methods 0.000 title claims abstract description 31
- 238000006243 chemical reaction Methods 0.000 title claims abstract description 29
- 230000001186 cumulative effect Effects 0.000 claims abstract description 13
- 238000002474 experimental method Methods 0.000 claims abstract description 13
- 238000005406 washing Methods 0.000 claims description 8
- 238000001035 drying Methods 0.000 claims description 5
- 229920006395 saturated elastomer Polymers 0.000 claims description 5
- 239000003350 kerosene Substances 0.000 claims description 4
- 238000004088 simulation Methods 0.000 claims description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 2
- 230000005311 nuclear magnetism Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000009466 transformation Effects 0.000 description 3
- 210000000867 larynx Anatomy 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明公开岩心核磁共振T2谱图驰豫时间与压汞孔喉半径转换方法，处理岩心，再对岩心进行核磁共振实验，测得核磁共振的横向驰豫时间T2与信号强度，并计算信号强度累积百分数；将岩心再处理以及进行高压压汞实验，测得孔喉半径r与进汞饱和度Shg，以及高压压汞曲线；在核磁共振横向驰豫时间T2的任意时刻T2(i)处的信号强度累积百分数为Sheci(i)，以Sheci(i)为标准对高压压汞曲线进行插值，得到与横向驰豫时间T2(i)对应的孔喉半径r(i)；分别选取插值得到的横向驰豫时间T2(i)与孔喉半径r(i)的两端值并对两端值进行分段拟合，求得两端处横向驰豫时间与孔喉半径的关系，得到的关系计算与剩余T2(i)对应的r(i)。该方法计算难度小、结果较准确，能够有效将核磁共振T2时间转换为孔喉半径。
Description
技术领域
本发明属于油气层物理技术领域，具体涉及岩心核磁共振T2谱图驰豫时间与压汞孔喉半径转换方法。
背景技术
核磁共振机理表明，当样品孔隙中充满流体时，H+在孔隙中作横向弛豫运动，H+与孔隙壁会产生碰撞作用，运动造成H+能量损失，使H+从高能态回到低能态，即H+的横向弛豫过程。碰撞越频繁，H+的能量损失越快，显然，孔隙的大小决定了H+与孔隙壁碰撞过程的多少，孔隙越小，H+在进行横向弛豫过程中与孔隙碰撞机率越大，因此可以认为孔隙的大小与H+的弛豫率呈反比关系，即孔隙越小，H+的横向弛豫率越高，对应的弛豫时间越短；孔隙越大，H+的横向弛豫率越低，对应的弛豫时间越长。同时，核磁共振弛豫谱分布范围越宽，幅度越低，孔喉分选性越差。这就是应用核磁共振资料研究岩石孔隙结构的理论基础。
目前，核磁共振横向驰豫时间与高压压汞孔喉半径的转换主要根据理论公式转换方法分为两种，一种是毛管压力图的转换，一种是孔喉半径频率图的转换，这两种方法的关键在于求取C值与n值。假定一个C值与n值，将核磁共振时间T2转换为孔喉半径r，将转换后的核磁共振曲线与高压压汞曲线放于同一个坐标图中。转换后的核磁共振曲线横坐标为r，纵坐标为核磁共振信号强度累积百分数；高压压汞曲线横轴为孔喉半径，纵轴为进汞饱和度。转换过程中不断变化C与n的值，使转换后的核磁曲线与高压压汞曲线最相似，即完成转换。这种转换方法计算复杂、步骤多，并且最终拟合结果误差大。
发明内容
为解决现有技术中存在的问题，本发明提供一种岩心核磁共振T2谱图驰豫时间与压汞孔喉半径转换方法，该方法计算难度小、结果较准确，能够有效将核磁共振T2时间转换为孔喉半径，有利于基于核磁共振技术对孔喉半径的研究。
本发明采用的技术方案如下：
岩心核磁共振T2谱图驰豫时间与压汞孔喉半径转换方法，包括如下步骤：
步骤1，将岩心洗油、烘干以及进行饱和，再对岩心进行核磁共振实验，测得核磁共振的横向驰豫时间T2与信号强度，将横向驰豫时间T2与信号强度按横向驰豫时间由大到小的顺序排列，并计算信号强度累积百分数；
步骤2，将岩心洗油、烘干以及进行高压压汞实验，测得孔喉半径r与进汞饱和度Shg，将进汞饱和度Shg转换为拟累积进汞饱和度S′hg，根据孔喉半径r与拟累积进汞饱和度S′hg获得高压压汞曲线；
所述拟累积进汞饱和度S′hg能够与核磁共振的信号强度累积百分数一致，达到100％；
步骤3，在核磁共振横向驰豫时间T2的任意时刻T2(i)处的信号强度累积百分数为Sheci(i)，以Sheci(i)为标准对高压压汞曲线进行插值，得到与横向驰豫时间T2(i)对应的孔喉半径r(i)；
步骤4，分别选取步骤3插值得到的横向驰豫时间T2(i)与孔喉半径r(i)的两端值并对两端值进行分段拟合，求得两端处横向驰豫时间与孔喉半径的关系，利用得到的关系计算与剩余T2(i)对应的r(i)。
步骤1中，取信号强度非零值的部分，将横向驰豫时间T2与信号强度按横向驰豫时间由大到小的顺序排列。
步骤2中，利用下式将进汞饱和度Shg转换为拟累积进汞饱和度S′hg；
式中，Shg(i)为某一孔喉半径r(i)处的进汞饱和度；S′hg(i)为某一孔喉半径r(i)处的拟累积进汞饱和度；∑Shg为高压压汞实验中最终进汞饱和度。
步骤4中，两端处横向驰豫时间与孔喉半径的关系如下：
r(i)＝CT2(i)n
式中，C和n均为预设常数。
岩心采用地层水或煤油进行饱和。
油为煤油或模拟油。
本发明具有如下有益效果：
本发明岩心核磁共振T2谱图驰豫时间与压汞孔喉半径转换方法以高压压汞数据为参考值，结合高压压汞数据，将核磁共振的横向驰豫时间T2值转换为相对应的孔喉半径值，从而做到对岩石孔喉半径的全面研究。由于每个岩心的核磁共振时间与孔喉半径的关系不是单纯的幂函数关系，在转换过程中，需要对插值结果进行分段拟合，其结果会更加准确。本发明的具体优点如下：第一：由核磁共振数据到高压压汞数据的方向插值，能得到更多的可供拟合的基础数据，使得拟合结果更为可靠；
第二：对插值结果进行分段拟合，即分别计算大孔喉与小孔喉的核磁共振的横向驰豫时间T2与孔喉半径r之间的关系，结果更可靠、准确；
第三：通过将核磁共振的横向驰豫时间T2转换成孔喉半径，能够更加详细、全面地研究孔喉半径分布，更有利于研究致密储层孔隙结构。
附图说明
图1为本发明实施例中的插值方法示意图。
具体实施方式
下面结合附图和实施例来对本发明做进一步的说明。
本发明的岩心核磁共振T2谱图驰豫时间与压汞孔喉半径转换方法，包括如下步骤：
步骤1，将岩心洗油、烘干以及进行饱和，再对岩心进行核磁共振实验，测得核磁共振的横向驰豫时间T2与信号强度，取信号强度非零值的部分，为了与高压压汞数据的孔喉半径对应，将横向驰豫时间T2与信号强度按横向驰豫时间由大到小的顺序排列，并计算信号强度累积百分数；
步骤2，将岩心洗油、烘干以及进行高压压汞实验，测得孔喉半径r与进汞饱和度Shg，为与核磁共振数据对照，利用公式(1)将进汞饱和度Shg转换为拟累积进汞饱和度S′hg，根据孔喉半径r与拟累积进汞饱和度S′hg获得高压压汞曲线；
式中，Shg(i)为某一孔喉半径r(i)处的进汞饱和度；S′hg(i)为某一孔喉半径r(i)处的拟累积进汞饱和度；∑Shg为高压压汞实验中最终进汞饱和度。
步骤3，在核磁共振横向驰豫时间T2的任意时刻T2(i)处的信号强度累积百分数为Sheci(i)，以Sheci(i)为标准对高压压汞曲线进行插值，得到与横向驰豫时间T2(i)对应的孔喉半径r(i)；
步骤4，由于步骤3插值只得到与部分T2(i)对应的r(i)，因此分别选取步骤3插值得到的横向驰豫时间T2(i)与孔喉半径r(i)的两端值并对两端值进行分段拟合，求得两端处横向驰豫时间与孔喉半径的关系r(i)＝CT2(i)n，利用得到的关系计算与剩余T2(i)对应的r(i)。
实施例
本实施例的岩心核磁共振T2谱图驰豫时间与压汞孔喉半径转换方法，包括如下步骤：
(1)选取一块砂岩岩心，对砂岩岩心依次进行洗油、烘干以及饱和煤油，然后进行核磁共振实验，测得核磁共振的横向驰豫时间T2与信号强度。取信号强度非零值的部分，为了与高压压汞数据的孔喉半径对应，将核磁共振的数据按横向驰豫时间T2由大到小的顺序排列，并计算信号强度累积百分数，具体如表1所示；
表1
(2)将岩心再次进行洗油和烘干，然后进行高压压汞实验，测得孔喉半径r与进汞饱和度Shg，为与核磁共振数据对照，将进汞饱和度Shg转换为拟累积进汞饱和度S′hg，具体的数据如表2所示。如图1所示，根据孔喉半径r与拟累积进汞饱和度S′hg获得高压压汞曲线。
表2
(3)如图1所示，利用高压压汞数据(拟累积进汞饱和度S′hg和孔喉半径r)，在核磁共振横向驰豫时间T2的任意时刻T2(i)处的信号强度累积百分数为Sheci(i)，以Sheci(i)为标准对压汞曲线进行插值，得到与横向驰豫时间T2(i)对应的孔喉半径r(i)，本步骤对应的参数如表3所示。
表3
(4)由于步骤3插值只得到与部分T2(i)对应的r(i)，因此分别选取步骤3插值得到的横向驰豫时间T2(i)与孔喉半径r(i)的两端值并对两端值进行分段拟合，求得两端处横向驰豫时间与孔喉半径的关系r(i)＝CT2(i)n，利用得到的关系计算与剩余T2(i)对应的r(i)，多的结果如表4所示。
表4
比较现有方法与本发明的转换结果，如表5所示：
现有转换方法是以压汞数据为标准对核磁共振数据进行插值，再求得横向驰豫时间与孔喉半径的关系，直接对所有数据进行拟合。
表5
由表5可以看出，对照高压压汞数据中拟累积进汞饱和度所对应的孔喉半径的大小，显然，本发明的转换结果更为准确；并且，在进行转换过程中，本发明直接进行插值、分段拟合、得出结果，计算简洁，但现有方法需要不断变换C与n的值以求最优解，计算复杂而误差大。
Claims (6)
1.岩心核磁共振T2谱图驰豫时间与压汞孔喉半径转换方法，其特征在于，包括如下步骤：
步骤1，将岩心洗油、烘干以及进行饱和，再对岩心进行核磁共振实验，测得核磁共振的横向驰豫时间T2与信号强度，将横向驰豫时间T2与信号强度按横向驰豫时间由大到小的顺序排列，并计算信号强度累积百分数；
步骤2，将岩心洗油、烘干以及进行高压压汞实验，测得孔喉半径r与进汞饱和度Shg，将进汞饱和度Shg转换为拟累积进汞饱和度S′hg，根据孔喉半径r与拟累积进汞饱和度S′hg获得高压压汞曲线；
所述拟累积进汞饱和度S′hg能够与核磁共振的信号强度累积百分数一致，达到100％；
步骤3，在核磁共振横向驰豫时间T2的任意时刻T2(i)处的信号强度累积百分数为Sheci(i)，以Sheci(i)为标准对高压压汞曲线进行插值，得到与横向驰豫时间T2(i)对应的孔喉半径r(i)；
步骤4，分别选取步骤3插值得到的横向驰豫时间T2(i)与孔喉半径r(i)的两端值并对两端值进行分段拟合，求得两端处横向驰豫时间与孔喉半径的关系，利用得到的关系计算与剩余T2(i)对应的r(i)。
2.根据权利要求1所述的岩心核磁共振T2谱图驰豫时间与压汞孔喉半径转换方法，其特征在于，步骤1中，取信号强度非零值的部分，将横向驰豫时间T2与信号强度按横向驰豫时间由大到小的顺序排列。
3.根据权利要求1所述的岩心核磁共振T2谱图驰豫时间与压汞孔喉半径转换方法，其特征在于，步骤2中，利用下式将进汞饱和度Shg转换为拟累积进汞饱和度S′hg；
式中，Shg(i)为某一孔喉半径r(i)处的进汞饱和度；S′hg(i)为某一孔喉半径r(i)处的拟累积进汞饱和度；∑Shg为高压压汞实验中最终进汞饱和度。
4.根据权利要求1所述的岩心核磁共振T2谱图驰豫时间与压汞孔喉半径转换方法，其特征在于，步骤4中，两端处横向驰豫时间与孔喉半径的关系如下：
r(i)＝CT2(i)n
式中，C和n均为预设常数。
5.根据权利要求1所述的岩心核磁共振T2谱图驰豫时间与压汞孔喉半径转换方法，其特征在于，岩心采用地层水或油进行饱和。
6.根据权利要求5所述的岩心核磁共振T2谱图驰豫时间与压汞孔喉半径转换方法，其特征在于，油为煤油或模拟油。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910340453.1A CN109916943B (zh) | 2019-04-25 | 2019-04-25 | 岩心核磁共振t2谱图驰豫时间与压汞孔喉半径转换方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910340453.1A CN109916943B (zh) | 2019-04-25 | 2019-04-25 | 岩心核磁共振t2谱图驰豫时间与压汞孔喉半径转换方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109916943A true CN109916943A (zh) | 2019-06-21 |
| CN109916943B CN109916943B (zh) | 2022-05-20 |
Family
ID=66978557
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910340453.1A Expired - Fee Related CN109916943B (zh) | 2019-04-25 | 2019-04-25 | 岩心核磁共振t2谱图驰豫时间与压汞孔喉半径转换方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109916943B (zh) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110231272A (zh) * | 2019-07-09 | 2019-09-13 | 中国地质大学（北京） | 致密砂岩孔径与核磁共振t2值转换关系的确定方法及系统 |
| CN110618158A (zh) * | 2019-10-28 | 2019-12-27 | 中石化石油工程技术服务有限公司 | 利用核磁共振信息构建岩心毛管压力曲线的方法 |
| CN111537544A (zh) * | 2020-06-08 | 2020-08-14 | 长安大学 | 一种提高核磁共振t2谱表征致密储层孔径分布精度的转换方法 |
| CN112067645A (zh) * | 2020-07-31 | 2020-12-11 | 中国石油大学(北京) | 核磁共振弛豫时间与孔喉半径的转换方法 |
| CN112284999A (zh) * | 2020-10-26 | 2021-01-29 | 中国石油大学(华东) | 砂岩孔径分布测定方法及其应用 |
| CN113008927A (zh) * | 2021-04-01 | 2021-06-22 | 中海石油(中国)有限公司 | 一种煤层核磁共振t2谱与孔隙分布的转换方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20120275658A1 (en) * | 2011-02-28 | 2012-11-01 | Hurley Neil F | Petrographic image analysis for determining capillary pressure in porous media |
| CN104634718A (zh) * | 2015-03-05 | 2015-05-20 | 中国石油大学(华东) | 应用核磁共振表征致密砂岩孔径分布的标定方法 |
| CN109060623A (zh) * | 2018-06-30 | 2018-12-21 | 西南石油大学 | 致密砂岩孔隙大小分布联合表征方法 |
-
2019
- 2019-04-25 CN CN201910340453.1A patent/CN109916943B/zh not_active Expired - Fee Related
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20120275658A1 (en) * | 2011-02-28 | 2012-11-01 | Hurley Neil F | Petrographic image analysis for determining capillary pressure in porous media |
| CN104634718A (zh) * | 2015-03-05 | 2015-05-20 | 中国石油大学(华东) | 应用核磁共振表征致密砂岩孔径分布的标定方法 |
| CN109060623A (zh) * | 2018-06-30 | 2018-12-21 | 西南石油大学 | 致密砂岩孔隙大小分布联合表征方法 |
Non-Patent Citations (2)
| Title |
|---|
| 刘春杨等: "用核磁共振和常规测井资料预测毛细管压力曲线方法", 《测井技术》 * |
| 李爱芬 等: "核磁共振研究致密砂岩孔隙结构的方法及应用", 《中国石油大学学报（自然科学版）》 * |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110231272A (zh) * | 2019-07-09 | 2019-09-13 | 中国地质大学（北京） | 致密砂岩孔径与核磁共振t2值转换关系的确定方法及系统 |
| CN110231272B (zh) * | 2019-07-09 | 2020-05-12 | 中国地质大学（北京） | 致密砂岩孔径与核磁共振t2值转换关系的确定方法及系统 |
| CN110618158A (zh) * | 2019-10-28 | 2019-12-27 | 中石化石油工程技术服务有限公司 | 利用核磁共振信息构建岩心毛管压力曲线的方法 |
| CN111537544A (zh) * | 2020-06-08 | 2020-08-14 | 长安大学 | 一种提高核磁共振t2谱表征致密储层孔径分布精度的转换方法 |
| CN111537544B (zh) * | 2020-06-08 | 2023-04-07 | 长安大学 | 一种提高核磁共振t2谱表征致密储层孔径分布精度的转换方法 |
| CN112067645A (zh) * | 2020-07-31 | 2020-12-11 | 中国石油大学(北京) | 核磁共振弛豫时间与孔喉半径的转换方法 |
| CN112067645B (zh) * | 2020-07-31 | 2022-07-05 | 中国石油大学(北京) | 核磁共振弛豫时间与孔喉半径的转换方法 |
| CN112284999A (zh) * | 2020-10-26 | 2021-01-29 | 中国石油大学(华东) | 砂岩孔径分布测定方法及其应用 |
| CN113008927A (zh) * | 2021-04-01 | 2021-06-22 | 中海石油(中国)有限公司 | 一种煤层核磁共振t2谱与孔隙分布的转换方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109916943B (zh) | 2022-05-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109916943A (zh) | 岩心核磁共振t2谱图驰豫时间与压汞孔喉半径转换方法 | |
| CN107894386B (zh) | 超临界二氧化碳注入对低渗透砂岩油藏孔喉结构影响的定量评价方法 | |
| CN103279647B (zh) | 一种基于孔隙特征参数计算致密砂岩储层渗透率的方法 | |
| CN104819923B (zh) | 基于核磁共振的低渗透砂岩储层孔隙结构定量反演方法 | |
| CN103654789B (zh) | 磁共振快速参数成像方法和系统 | |
| CN106093099B (zh) | 一种获得高分辨二维j分解谱的方法 | |
| CN103900942B (zh) | 基于高压压汞分析的储层微观孔喉参数的连续表征方法 | |
| CN105092629B (zh) | 一种测量氢‑氢j耦合常数的磁共振二维谱方法 | |
| CN110296931B (zh) | 一种致密砂岩油水相对渗透率信息的表征方法及系统 | |
| CN106050225A (zh) | 一种核磁共振测井t2谱100%纯水谱的确定方法 | |
| CN104278989A (zh) | 一种获取低孔低渗储层饱和度指数的方法 | |
| CN104268883A (zh) | 一种基于边缘检测的时频谱曲线提取方法 | |
| Ji et al. | Flow pattern identification based on EMD and LS-SVM for gas–liquid two-phase flow in a minichannel | |
| CN102425409B (zh) | 一种井斜和工具面测量装置及垂直钻井设备 | |
| CN109100812A (zh) | 基于核磁共振的岩石孔隙分形维数评价方法及装置 | |
| CN104405374A (zh) | 一种致密气藏储层应力敏感性的测量方法 | |
| CN107221334A (zh) | 一种音频带宽扩展的方法及扩展装置 | |
| CN111827996A (zh) | 基于力学性质的多参数综合定性致密砂岩储层分类方法 | |
| CN109856688A (zh) | 基于核磁测井双tw极化增强法的流体性质识别方法 | |
| CN104155621B (zh) | 一种准确测量静磁场b0分布的方法 | |
| CN103941205B (zh) | 一种不均匀场下改善核磁共振氢谱分辨率的方法 | |
| CN103334740A (zh) | 考虑启动压力梯度的确定泄油前缘的方法 | |
| CN107843531A (zh) | 低渗透低粘油油层孔隙结构的核磁共振参数表征方法 | |
| CN106845685B (zh) | 一种降低水平井试油压裂作业成本的参数优化方法 | |
| CN109617051A (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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20220520 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |