CN111155974A - 缝洞型油藏注水替油的油水置换状态预测方法 - Google Patents
缝洞型油藏注水替油的油水置换状态预测方法 Download PDFInfo
- Publication number
- CN111155974A CN111155974A CN202010090575.2A CN202010090575A CN111155974A CN 111155974 A CN111155974 A CN 111155974A CN 202010090575 A CN202010090575 A CN 202010090575A CN 111155974 A CN111155974 A CN 111155974A
- Authority
- CN
- China
- Prior art keywords
- oil
- water
- grid
- saturation
- well
- 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
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 144
- 238000006073 displacement reaction Methods 0.000 title claims abstract description 54
- 238000000034 method Methods 0.000 title claims abstract description 38
- 238000002347 injection Methods 0.000 title claims abstract description 24
- 239000007924 injection Substances 0.000 title claims abstract description 24
- 230000005484 gravity Effects 0.000 claims abstract description 35
- 230000008859 change Effects 0.000 claims abstract description 15
- 230000035699 permeability Effects 0.000 claims abstract description 11
- 238000012360 testing method Methods 0.000 claims abstract description 5
- 239000003921 oil Substances 0.000 claims description 54
- 238000004364 calculation method Methods 0.000 claims description 14
- 238000004891 communication Methods 0.000 claims description 11
- 239000010779 crude oil Substances 0.000 claims description 9
- 239000007788 liquid Substances 0.000 claims description 9
- 239000000463 material Substances 0.000 claims description 5
- 230000001133 acceleration Effects 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 238000012935 Averaging Methods 0.000 claims description 2
- 230000008569 process Effects 0.000 abstract description 8
- 208000010392 Bone Fractures Diseases 0.000 description 9
- 206010017076 Fracture Diseases 0.000 description 9
- 239000012071 phase Substances 0.000 description 8
- 239000012530 fluid Substances 0.000 description 3
- 238000013178 mathematical model Methods 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 239000008346 aqueous phase Substances 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000003129 oil well Substances 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 235000020681 well water Nutrition 0.000 description 1
- 239000002349 well water 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
- E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
- E21B43/16—Enhanced recovery methods for obtaining hydrocarbons
- E21B43/20—Displacing by water
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Mining & Mineral Resources (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Computational Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Fluid Mechanics (AREA)
- Pure & Applied Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Operations Research (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提出了缝洞型油藏注水替油的油水置换状态预测方法,油水的置换首先发生在水体的边缘,再逐步发展到井点处,基于重力差建立油水置换流动模型,将井区划分为多个网格,利用试井解释的溶洞渗透率,计算井区的油水饱和度分布随闷井时间的变化,从而判断油水置换程度,优化闷井期的长短;本发明能更好的描述油水置换基本过程,对含水饱和度在分布形态和变化模式上的预测也更加准确。
Description
技术领域
本发明涉及缝洞型油藏注水开发领域,尤其涉及缝洞型油藏注水替油的油水置换状态预测方法。
背景技术
缝洞型油藏的溶洞孔隙尺度大,油水在溶洞内重力分异作用显著,利用该特点形成了缝洞型油藏特有的注水替油工艺,即对单井进行周期注水、闷井置换、吞吐采油。大排量注水补充地层能量后,在井周形成类似径向的水体,人工水体在闷井期依靠重力作用下沉,实现油水的空间置换,再开井采出上部的原油。
注水替油工艺实施的关键参数之一是闷井时间;闷井时间过短,油水置换不足,注入水回吐量高,替油效率低;闷井时间过长,替油效率虽高,但油井的生产时效又低。因此,如何预测油水的置换状态、优化闷井时间成为注水替油工艺设计的关键技术。
发明内容
本发明的目的在于,针对上述问题,提出缝洞型油藏注水替油的油水置换状态预测方法,可有效判断油水置换程度、优化闷井时间。
缝洞型油藏注水替油的油水置换状态预测方法,油水的置换首先发生在水体的边缘,再逐步发展到井点处,基于重力差建立油水置换流动模型,将井区划分为多个网格,利用试井解释的溶洞渗透率,计算井区的含水饱和度分布随闷井时间的变化,从而判断油水置换程度、优化闷井时间。
进一步的,所述油水置换流动模型的建立包括以下步骤:
S1:计算网格i与网格i+1间重力压差产生的横向流动;
S2:建立水体径向流动置换模型,并根据置换模型建立水相的横向流动方程;
S3:将重力流速vg及其梯度代入水相的横向流动方程,得到饱和度方程。
进一步的,所述步骤S1包括如下子步骤:
S11:计算网格i与网格i+1间液柱重力压差的最大值并平均化;
S12:取平均重力压差为网格间的重力压差;
S13:计算重力压差形成的网格压力梯度和横向压力梯度;
S14:计算重力流速vg以表示重力压差产生的横向流动。
进一步的,所述网格间液柱重力压差的最大值ΔPg,max的计算方法为:
ΔPg,max=h(ρw-ρo)g(Sw,i+1-Sw,i);
其中,含水饱和度Sw,i用于表示网格i中的油水分布比例,ρo为原油密度,ρw为注入水密度,h为网格液柱高度,g为重力加速度,Δx为网格步长,x为沿流动方向的距离坐标,Pg为重力压强。
进一步的,所述重力流速vg的计算方法如下:
其中,k为溶洞等效渗透率,μo为原油粘度。
进一步的,所述水体径向流动置换模型是以井点为原点的径向坐标系,在半径r至r+Δr的圆环微元体中,根据物质平衡关系,建立水相的横向流动方程如下:
其中含水饱和度Sw用于表示圆环微元体中油水分布比例,φ为溶洞的孔隙度,t为闷井时间。
进一步的,所述饱和度方程为:
进一步的,计算井区的含水饱和度分布随闷井时间的变化的方法如下:
S41:将饱和度方程离散化并进行第n+1时步的饱和度计算;
S42:对边界条件进行处理。
进一步的,所述步骤S41通过差分法数值求解,将饱和度方程进行离散化,得到:
第n+1时步的饱和度计算方法为:
进一步的,所述对边界条件进行处理包括以下子步骤:
本发明的有益效果:本方法基于重力差建立油水置换的一维流动模型,能更好的描述油水置换基本过程,对含水饱和度在分布形态和变化模式上的预测更加准确,为判断油水置换程度、优化闷井时间提供了一种更加简洁的技术方案。
附图说明
图1是本发明网格间油水组成不同时的重力差示意图;
图2是本发明微元体Δr环中的物质平衡关系图;
图3是本发明闷井期水饱和度分布随时间的变化图;
图4是本发明井点的含水饱和度随闷井时间的变化图。
具体实施方式
为了对本发明的技术特征、目的和效果有更加清楚的理解,现对照附图说明本发明的具体实施方式。
缝洞型油藏注水替油的油水置换状态预测方法,油水的置换首先发生在水体的边缘,再逐步发展到井点处,基于重力差建立油水置换流动模型,将井区划分为多个网格,利用试井解释的溶洞渗透率,计算井区的含水饱和度分布随闷井时间的变化,从而判断油水置换程度、优化闷井时间。
所述油水置换流动模型的目的在于建立数学模型以描述油水置换基本过程,其通过基本油水置换模型和水体径向流动置换模型得到油水置换的水饱和度方程,如图1所示的基本油水置换模型,假设注水期近井区的水已经沉降到溶洞底部,以相邻2个网格来分析油水的置换流动,以含水饱和度Sw,i形式表示网格i中的油水分布比例,原油的密度为ρo,注入水的密度为ρw,网格液柱高度为h,g为重力加速度,则网格间液柱重力压差的最大值为:
ΔPg,max=h(ρw-ρo)g(Sw,i+1-Sw,i);
取平均重力压差为网格间压差:
由重力压差形成的网格压力梯度为:
其中,Δx为网格步长。
重力压差形成的横向压力梯度可表示为:
其中,x为沿流动方向的距离坐标。
利用达西定律近似描述重力压差产生的横向流动,称为重力流速vg:
其中,k为溶洞等效渗透率,μo为原油粘度。由于网格净流出的水等于净流入的油,油相的流动性相对较差,油相的流动能力控制了流体的交换速度,因此,这里用油的粘度μo来计算重力流速。
如图2所示的微元体Δr环中的物质平衡关系,是以井点为原点的径向坐标系中,在半径r至r+Δr的圆环微元体中,建立的水体径向流动置换模型,该模型中根据物质平衡关系:流入质量-流出质量=内部的质量变化率,建立水相的横向流动方程为:
其中,φ为溶洞的孔隙度,t为闷井时间。
展开水相的横向流动方程,忽略Δr的二阶、三阶项,将其简化为:
其中重力流速vg及其梯度表示为:
将重力流速vg及其梯度代入水相的横向流动方程得:
整理得闷井期人工水体的分布方程,即饱和度方程为:
计算井区的含水饱和度分布随闷井时间的变化的方法如下:
S1:将饱和度方程离散化并进行第n+1时步的饱和度计算;
S2:对边界条件进行处理。
饱和度方程难以解析求解,可通过差分法数值求解。将饱和度方程离散化得到:
第n+1时步的饱和度计算方法为:
所述边界条件处理包括:
当取实用单位制时,渗透率k为μm2;油藏厚度h为m;流体密度为kg/m3;流体粘度为mPa.s;时间t为天。系数C的计算方式为:
在一个具体实施例中,取溶洞等效渗透率为5μm2,液柱高度20m,溶洞孔隙度45%,原油粘度5mPa.s,原油密度850kg/m3,注入水密度1150kg/m3,时间步长0.02天,网格步长2m,注水期水淹的网格数量10,网格总数量15,对应的注水量11310m3、溶洞体积25447m3。
根据第n+1时步的饱和度计算方法计算出闷井期0天~90天的含水饱和度分布曲线如图3所示,井点的含水饱和度随闷井时间的变化曲线如图4所示。
图3反映出置换流动首先发生在水体边缘,逐步向井筒端延伸;初期油水高度差异大、置换速度高,后期随水体的下降,油水高度差变小、置换速度下降;井筒端即井点处的水体置换最晚,井点的含水饱和度相对其它位置最高;长期闷井后含水饱和度分布将变成水平线。
含水饱和度分布的计算结果表明:本案所建立的数学模型即饱和度方程描述了油水置换基本过程,计算的饱和度分布形态、变化模式正常。
利用预测的井点含水饱和度及饱和度下降率随闷井时间变化曲线,可以辅助优选闷井时间,如图4所示,闷井9天含水饱和度下降率达到最大值,此后含水饱和度下降率逐步变小,3倍的下降率最大值时间以后,井点的含饱和度变化缓慢,饱和度下降率低于最大值的1/2,该规律在不同参数下的计算情况类似。因此,从时效上考虑推荐最优闷井时间为3倍下降率最大值时间。
本方案提出一种预测缝洞型油藏注水替油的井区油水置换状态的计算方法,基于重力差建立油水置换流动模型,将井区划分为多个网格,利用试井解释的溶洞渗透率,计算井区的油水饱和度分布随闷井时间的变化,从而判断油水置换程度,优化闷井期的长短;能更好的描述油水置换基本过程,对含水饱和度在分布形态和变化模式上的预测也更加准确。
以上显示和描述了本发明的基本原理和主要特征和本发明的优点。本行业的技术人员应该了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的只是说明本发明的原理,在不脱离本发明精神和范围的前提下,本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书界定。
Claims (10)
1.缝洞型油藏注水替油的油水置换状态预测方法,其特征在于,基于重力差建立油水置换流动模型,将井区划分为多个网格,利用试井解释的溶洞渗透率,计算井区的含水饱和度分布随闷井时间的变化,从而判断油水置换程度、优化闷井时间。
2.根据权利要求1所述的缝洞型油藏注水替油的油水置换状态预测方法,其特征在于,所述油水置换流动模型的建立包括以下步骤:
S1:计算网格i与网格i+1间重力压差产生的横向流动;
S2:建立水体径向流动置换模型,并根据置换模型建立水相的横向流动方程;
S3:将重力流速vg及其梯度代入水相的横向流动方程,得到饱和度方程。
3.根据权利要求2所述的缝洞型油藏注水替油的油水置换状态预测方法,其特征在于,所述步骤S1包括如下子步骤:
S11:计算网格i与网格i+1间液柱重力压差的最大值并平均化;
S12:取平均重力压差为网格间的重力压差;
S13:计算重力压差形成的网格压力梯度和横向压力梯度;
S14:计算重力流速vg以表示重力压差产生的横向流动。
8.根据权利要求1或7所述的缝洞型油藏注水替油的油水置换状态预测方法,其特征在于,计算井区的含水饱和度分布随闷井时间的变化的方法如下:
S41:将饱和度方程离散化并进行第n+1时步的饱和度计算;
S42:对边界条件进行处理。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010090575.2A CN111155974B (zh) | 2020-02-13 | 2020-02-13 | 缝洞型油藏注水替油的油水置换状态预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010090575.2A CN111155974B (zh) | 2020-02-13 | 2020-02-13 | 缝洞型油藏注水替油的油水置换状态预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111155974A true CN111155974A (zh) | 2020-05-15 |
CN111155974B CN111155974B (zh) | 2022-03-11 |
Family
ID=70565600
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010090575.2A Expired - Fee Related CN111155974B (zh) | 2020-02-13 | 2020-02-13 | 缝洞型油藏注水替油的油水置换状态预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111155974B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113586023A (zh) * | 2021-07-26 | 2021-11-02 | 中国石油大学(北京) | 页岩油藏压后闷井时间的确定方法及设备 |
CN117077574A (zh) * | 2023-10-16 | 2023-11-17 | 西安石油大学 | 一种缝洞油藏模型驱替机理定量表征方法及装置 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE3300685A1 (de) * | 1983-01-11 | 1984-07-12 | Mobil Oil Corp., New York, N.Y. | Verfahren zur bestimmung der fluessigkeitssaettigung in einer untergrundformation |
SU1571229A1 (ru) * | 1988-07-05 | 1990-06-15 | Украинский Научно-Исследовательский Геологоразведочный Институт | Способ определени водонасыщенности горных пород |
CN104747180A (zh) * | 2014-02-26 | 2015-07-01 | 中国石油化工股份有限公司 | 一种用于注水开发的缝洞型油藏分析方法及其应用 |
CN105587298A (zh) * | 2015-12-17 | 2016-05-18 | 西南石油大学 | 多流管模式的水驱油井含水率反演方法 |
CN105676309A (zh) * | 2016-01-28 | 2016-06-15 | 中国石油大学(华东) | 裂缝性油藏逆向渗吸采收率预测方法 |
CN207177792U (zh) * | 2017-06-28 | 2018-04-03 | 中国石油天然气股份有限公司 | 堵塞器用投捞工具 |
CN109060863A (zh) * | 2018-08-20 | 2018-12-21 | 西安石油大学 | 一种评价页岩气储层自吸前缘运移能力的方法 |
CN109555514A (zh) * | 2018-11-07 | 2019-04-02 | 西南石油大学 | 缝洞型油藏的示踪剂解释模型及其构建方法与应用 |
CN110334868A (zh) * | 2019-07-08 | 2019-10-15 | 西南石油大学 | 一种耦合流体流动与地质应力预测最优焖井时间的方法 |
-
2020
- 2020-02-13 CN CN202010090575.2A patent/CN111155974B/zh not_active Expired - Fee Related
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE3300685A1 (de) * | 1983-01-11 | 1984-07-12 | Mobil Oil Corp., New York, N.Y. | Verfahren zur bestimmung der fluessigkeitssaettigung in einer untergrundformation |
SU1571229A1 (ru) * | 1988-07-05 | 1990-06-15 | Украинский Научно-Исследовательский Геологоразведочный Институт | Способ определени водонасыщенности горных пород |
CN104747180A (zh) * | 2014-02-26 | 2015-07-01 | 中国石油化工股份有限公司 | 一种用于注水开发的缝洞型油藏分析方法及其应用 |
CN105587298A (zh) * | 2015-12-17 | 2016-05-18 | 西南石油大学 | 多流管模式的水驱油井含水率反演方法 |
CN105676309A (zh) * | 2016-01-28 | 2016-06-15 | 中国石油大学(华东) | 裂缝性油藏逆向渗吸采收率预测方法 |
CN207177792U (zh) * | 2017-06-28 | 2018-04-03 | 中国石油天然气股份有限公司 | 堵塞器用投捞工具 |
CN109060863A (zh) * | 2018-08-20 | 2018-12-21 | 西安石油大学 | 一种评价页岩气储层自吸前缘运移能力的方法 |
CN109555514A (zh) * | 2018-11-07 | 2019-04-02 | 西南石油大学 | 缝洞型油藏的示踪剂解释模型及其构建方法与应用 |
CN110334868A (zh) * | 2019-07-08 | 2019-10-15 | 西南石油大学 | 一种耦合流体流动与地质应力预测最优焖井时间的方法 |
Non-Patent Citations (2)
Title |
---|
廖书宇: ""缝洞型碳酸盐岩油藏高效吞吐开发研究"", 《中国优秀硕士学位论文全文数据库 工程科技I辑》 * |
杨帆 等: ""碳酸盐岩油藏注水吞吐闷井时间初探"", 《内蒙古石油化工》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113586023A (zh) * | 2021-07-26 | 2021-11-02 | 中国石油大学(北京) | 页岩油藏压后闷井时间的确定方法及设备 |
CN117077574A (zh) * | 2023-10-16 | 2023-11-17 | 西安石油大学 | 一种缝洞油藏模型驱替机理定量表征方法及装置 |
CN117077574B (zh) * | 2023-10-16 | 2024-02-23 | 西安石油大学 | 一种缝洞油藏模型驱替机理定量表征方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN111155974B (zh) | 2022-03-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104504230B (zh) | 低渗气井采收率与极限动用半径估计方法 | |
CN111155974B (zh) | 缝洞型油藏注水替油的油水置换状态预测方法 | |
CN105243182A (zh) | 致密油压裂水平井的动态储量计算方法 | |
CN107066769B (zh) | 适用于超深层裂缝型碳酸盐岩储层的高效酸化设计方法 | |
CN107832481B (zh) | 一种用于稠油油藏组合蒸汽吞吐的分区方法 | |
CN102507412A (zh) | 碳酸盐岩油藏等效模型的裂缝-基质渗透率级差判别方法 | |
CN104727790A (zh) | 水驱油藏剩余油再富集成藏周期预测方法 | |
CN115270662A (zh) | 一种欠饱和储层煤层气井生产边界定量预测的计算方法 | |
CN111502639A (zh) | 一种救援井最小压井排量的确定方法 | |
CN113360984B (zh) | 一种考虑壁面捕获效应的支撑剂输送数值模拟方法 | |
CN114528780A (zh) | 一种考虑稠油非牛顿特性的数值模拟方法 | |
CN110924935A (zh) | 致密油藏井底流压调控方案确定方法、装置和设备 | |
CN109403957B (zh) | 一种高压地层压力获取方法 | |
CN112377182B (zh) | 大型孔洞裂缝型碳酸盐岩储集体参数确定方法及装置 | |
CN111608647B (zh) | 一种注水注聚井区地层压力预测方法 | |
CN104407397A (zh) | 一种致密油气储层物性下限的确定方法 | |
CN110242261B (zh) | 垂向注气孔隙油气水微观渗流规律的预测方法及系统 | |
CN108843384B (zh) | 海洋超浅层油气的长水平井砾石充填方法 | |
CN108386170B (zh) | 一种油藏开发过程中地下能耗表征方法 | |
Sagen et al. | A coupled dynamic reservoir and pipeline model–development and initial experience | |
CN106227903B (zh) | 底水层钻开厚度的确定方法及装置 | |
CN111636454B (zh) | 承压水地层控制地下结构沉降的完整井水面设计方法 | |
CN110608031B (zh) | 一种井下节流器选井方法 | |
CN108229713A (zh) | 断块油藏多层合采方案优化设计方法 | |
Taghizadeh et al. | Geomechanical and thermal reservoir simulation during steam flooding |
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: 20220311 |
|
CF01 | Termination of patent right due to non-payment of annual fee |