CN111693425B - 一种基于压汞曲线的岩心薄膜束缚水饱和度的测定方法 - Google Patents

一种基于压汞曲线的岩心薄膜束缚水饱和度的测定方法 Download PDF

Info

Publication number
CN111693425B
CN111693425B CN201910192711.6A CN201910192711A CN111693425B CN 111693425 B CN111693425 B CN 111693425B CN 201910192711 A CN201910192711 A CN 201910192711A CN 111693425 B CN111693425 B CN 111693425B
Authority
CN
China
Prior art keywords
mercury
saturation
film
pore throat
formula
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.)
Expired - Fee Related
Application number
CN201910192711.6A
Other languages
English (en)
Other versions
CN111693425A (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 University of Petroleum Beijing
Original Assignee
China University of Petroleum Beijing
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 University of Petroleum Beijing filed Critical China University of Petroleum Beijing
Priority to CN201910192711.6A priority Critical patent/CN111693425B/zh
Publication of CN111693425A publication Critical patent/CN111693425A/zh
Application granted granted Critical
Publication of CN111693425B publication Critical patent/CN111693425B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/08Investigating permeability, pore-volume, or surface area of porous materials
    • G01N15/088Investigating volume, surface area, size or distribution of pores; Porosimetry

Landscapes

  • Chemical & Material Sciences (AREA)
  • Dispersion Chemistry (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

本发明公开了基于压汞曲线的岩心薄膜束缚水饱和度的测定方法。所述方法包括如下步骤:对岩心进行压汞实验,得到不同进汞压力下的进汞饱和度,进而得到岩心的孔喉半径和孔吼分布频率;根据进汞压力得到岩心的薄膜束缚水厚度;根据孔喉半径得到岩心的不同进汞区间孔喉平均半径;根据孔喉分布频率、薄膜束缚水厚度和不同进汞区间孔喉平均半径得到岩心的不同进汞区间薄膜束缚水饱和度,进而得到岩心的薄膜束缚水饱和度。本发明可以根据实验压汞曲线数据,快速得到岩心薄膜束缚水饱和度,具有简便实用、可操作性强的特点。本发明薄膜束缚水饱和度的测定结果可应用于储层评价、产能预测、储量计算以及油气藏出水机理研究等众多方面,具有广阔的应用前景。

Description

一种基于压汞曲线的岩心薄膜束缚水饱和度的测定方法
技术领域
本发明涉及一种基于压汞曲线的岩心薄膜束缚水饱和度的测定方法,属于油气田生产技术领域。
背景技术
束缚水是滞留在储层岩石微小孔隙或吸附于岩石颗粒表面的不可流动的水。束缚水饱和度为束缚水所占孔隙的体积与岩石总孔隙体积之比。束缚水饱和度是油藏单相与油水两相渗流的临界参数,也是进行储层评价、产能预测、油水渗流规律分析和储量计算的基础参数。
束缚水包括两种形式:一种是滞留在储层岩石微小孔隙中的毛细管束缚水(如图1a所示);另一种是吸附于岩石颗粒表面(或大孔隙壁上)的薄膜束缚水(如图1b所示)。因此,束缚水饱和度也可分为两种:即毛细管束缚水饱和度和薄膜束缚水饱和度。目前,束缚水饱和度主要通过压汞、离心毛管力、油水相渗、核磁共振等流体实验方法测定。由于这些实验手段不能对毛细管束缚水和薄膜束缚水进行有效区分,致使目前薄膜束缚水饱和度尚无法准确定量的测定,这在很大程度上制约了薄膜束缚水的理论研究和实践应用。因此,建立一种薄膜束缚水饱和度的测定方法具有重要的理论意义和实践价值。
发明内容
本发明的目的是提供一种基于压汞曲线的岩心薄膜束缚水饱和度的测定方法,具有简便实用、可操作性强的特点。
本发明所提供的岩心薄膜束缚水饱和度的测定方法,包括如下步骤:
1)对岩心进行压汞实验,得到不同进汞压力下的进汞饱和度,进而得到岩心的孔喉半径和孔吼分布频率;
2)根据所述进汞压力得到岩心的薄膜束缚水厚度;
3)根据所述孔喉半径得到岩心的不同进汞区间孔喉平均半径;
4)根据所述孔喉分布频率、所述薄膜束缚水厚度和所述不同进汞区间孔喉平均半径得到岩心的不同进汞区间薄膜束缚水饱和度,进而得到岩心的薄膜束缚水饱和度。
上述的测定方法中,步骤1)中,所述压汞实验过程中,通过测量不同进汞压力下的进汞饱和度,便可绘制压汞曲线,如图2所示,进而确定岩心孔喉分布,具体地,可按照下述公式得到:
根据所述进汞压力,结合式(1)得到所述孔喉半径;
Figure BDA0001994842690000021
式(1)中,ri表示第i个进汞压力对应的孔喉半径,μm;σ表示汞表面张力,为0.48N/m;θ表示汞润湿角,为140°;Pi表示第i个测试点进汞压力,MPa;n表示进汞压力测试总点数。
根据所述进汞饱和度,结合式(2)得到所述孔喉分布频率;
fi=ΔSi=Si+1-Si (i=1,2,……,n-1;) (2)
式(2)中,fi表示第i个区间孔喉分布频率,%;ΔSi表示第i个区间汞饱和度增量,%;Si+1表示第i+1个测试点进汞压力对应的汞饱和度,%;Si表示第i个测试点进汞压力对应的汞饱和度,%;n表示进汞压力测试总点数。
上述的测定方法中,步骤2)中,根据式(3)得到所述薄膜束缚水厚度;
Figure BDA0001994842690000022
式(3)中,Pn表示第n个测试点进汞压力,MPa;rn表示束缚水薄膜厚度,μm;σ表示汞表面张力,为0.48N/m;θ表示汞润湿角,为140°。
上述的测定方法中,步骤3)中,根据式(4)得到所述不同进汞区间孔喉平均半径;
Figure BDA0001994842690000023
式(4)中,rci表示第i个进汞区间孔喉平均半径,μm;ri表示第i个测试点进汞压力对应的孔喉半径,μm;ri+1表示第i+1个测试点进汞压力对应的孔喉半径,μm;n表示进汞压力测试总点数。
上述的测定方法中,步骤4)中,建立了薄模束缚水同心毛细管模型,如图3所示,进而根据式(5)得到所述不同进汞区间薄膜束缚水饱和度;
Figure BDA0001994842690000024
式(5)中,Swni表示第i个进汞区间薄膜束缚水饱和度,%;rci表示第i个进汞区间孔喉平均半径,μm;rn表示束缚水薄膜厚度,μm;ΔSi表示第i个进汞区间汞饱和度增量,即第i个区间孔喉分布频率,%;n表示进汞压力测试总点数。
上述的测定方法中,步骤4)中,根据式(6)得到所述薄膜束缚水饱和度;
Figure BDA0001994842690000031
式(6)中,Swni表示第i个进汞区间薄膜束缚水饱和度,%;Swn:岩心薄膜束缚水饱和度,%,n表示进汞压力测试总点数。
本发明由于采取以上技术方案,其具有以下优点:
1、本发明可以根据实验压汞曲线数据,快速得到岩心薄膜束缚水饱和度,具有简便实用、可操作性强的特点。
2、本发明有效弥补了目前实验方法不能对薄膜束缚水饱和度进行定量测定的缺点。
3、基于本发明方法,可以进一步对薄膜束缚水饱和度和毛细管束缚水饱和度构成进行定量评价。
4、本发明薄膜束缚水饱和度的测定结果可应用于储层评价、产能预测、储量计算以及油气藏出水机理研究等众多方面,具有广阔的应用前景。
附图说明
图1为束缚水的两种形式示意图,其中,图1a表示毛细管束缚水,图1b表示薄膜束缚水。
图2为由压汞曲线确定孔喉分布的示意图。
图3为薄模束缚水同心毛细管的模型图。
图4为本发明方法的流程图。
图5为本发明实施例中的压汞曲线确定孔喉分布图。
具体实施方式
下述实施例中所使用的实验方法如无特殊说明,均为常规方法。
下述实施例中所用的材料、试剂等,如无特殊说明,均可从商业途径得到。
如图4所示,本发明提供的基于压汞曲线的岩心薄膜束缚水饱和度的测定方法包括以下步骤:
1、根据压汞曲线数据,确定岩心孔喉分布
由压汞曲线数据,确定岩心孔喉分布计算公式如下:
1.1孔喉半径计算
设孔喉半径为ri,其计算公式为:
Figure BDA0001994842690000041
式中,Pi:第i个测试点进汞压力(以下简称第i个进汞压力),MPa;ri:第i个进汞压力对应的孔喉半径,μm;σ:汞表面张力,0.48mN/m;θ:汞润湿角,140°;n:进汞压力测试总点数。
1.2孔喉分布频率计算
压汞实验中两相邻测试点进汞压力(Pi,Pi+1)构成一个进汞区间i(简称区间)。其中,孔喉分布频率(fi)等于区间汞饱和度增量ΔSi,其计算公式为:
fi=ΔSi=Si+1-Si(i=1,2,……,n-1;)
式中,fi:第i个区间孔喉分布频率,%;ΔSi:第i个区间汞饱和度增量,%;Si+1:第i+1个进汞压力Pi+1对应的汞饱和度,%;Si:第i个进汞压力Pi对应的汞饱和度,%;n:进汞压力测试总点数。
2、确定薄膜束缚水厚度
设束缚水薄膜厚度为rn,其计算公式为:
Figure BDA0001994842690000042
式中,Pn:第n个进汞压力,MPa;rn:束缚水薄膜厚度,μm;σ:汞表面张力,0.48N/m;θ:汞润湿角,140°。
3、求取不同进汞区间薄膜束缚水饱和度
3.1计算不同进汞区间孔喉平均半径
设区间孔喉平均半径为rci,其计算公式为:
Figure BDA0001994842690000043
式中,rci:第i个进汞区间孔喉平均半径,μm;ri:第i个进汞压力对应的孔喉半径,μm;ri+1:第i+1个进汞压力对应的孔喉半径,μm;n:进汞压力测试总点数。
3.2计算不同进汞区间薄膜束缚水饱和度
建立了薄模束缚水同心毛细管模型,如图3所示,设区间薄膜束缚水饱和度为Swni,其计算公式为:
Figure BDA0001994842690000044
式中,rci:第i个进汞区间孔喉平均半径,μm;rn:束缚水薄膜厚度,μm;Swni:第i个进汞区间薄膜束缚水饱和度,%;ΔSi:第i个进汞区间汞饱和度增量,%;n:进汞压力测试总点数。
4、确定岩心薄膜束缚水饱和度
设岩心薄膜束缚水饱和度为Swn,其计算公式为:
Figure BDA0001994842690000051
式中,Swni:第i个进汞区间薄膜束缚水饱和度,%;Swn:岩心薄膜束缚水饱和度,%,n表示进汞压力测试总点数。
下面通过一个具体的实施例,进一步说明本发明的技术效果。
南梁油田位于鄂尔多斯盆地陕北斜坡西南部,地处甘肃省华池、庆阳县境内,面积约2600km2。主力层为三叠系上统延长组长6油组长63砂组,属于深湖-半深湖相重力流沉积,储层平均孔隙度9.09%、平均渗透率0.213×10-3μm2,属于特低孔-超低渗储层。以南梁油田长6油组山156井18#岩心(2060.1m)压汞曲线为实施例,该岩样压汞实验基础数据如表1所示,压汞实验测试数据如表2所示。
表1压汞实验基础数据表
Figure BDA0001994842690000052
表2压汞实验测试数据表
Figure BDA0001994842690000053
Figure BDA0001994842690000061
1)根据压汞曲线数据,确定岩心孔喉分布
由压汞曲线数据,确定岩心孔喉分布计算公式如下:
①孔喉半径计算
设孔喉半径为ri,其计算公式为:
Figure BDA0001994842690000062
式中,Pi:第i个测试点进汞压力(以下简称第i个进汞压力),MPa;ri:第i个进汞压力对应的孔喉半径,μm;σ:汞表面张力,0.48mN/m;θ:汞润湿角,140°;n:进汞压力测试总点数,计算结果如表3所示。
②孔喉分布频率计算
压汞实验中两相邻测试点进汞压力(Pi,Pi+1)构成一个进汞区间i(简称区间)。其中,孔喉分布频率(fi)等于区间汞饱和度增量ΔSi,其计算公式为:fi=ΔSi=Si+1-Si(i=1,2,……,n-1;)
式中,fi:第i个区间孔喉分布频率,%;ΔSi:第i个区间汞饱和度增量,%;Si+1:第i+1个进汞压力Pi+1对应的汞饱和度,%;Si:第i个进汞压力Pi对应的汞饱和度,%;n:进汞压力测试总点数,计算结果如表3所示,由18#岩心压汞曲线确定孔喉分布如图5。
2)确定薄膜束缚水厚度
设束缚水薄膜厚度为rn,其计算公式为:
Figure BDA0001994842690000063
式中,Pn:第n个进汞压力,MPa;rn:束缚水薄膜厚度,μm;σ:汞表面张力,0.48N/m;θ:汞润湿角,140°,取n=29,计算结果如表3所示。
3)求取不同进汞区间薄膜束缚水饱和度
①计算不同进汞区间孔喉平均半径
设区间孔喉平均半径为rci,其计算公式为:
Figure BDA0001994842690000071
式中,rci:第i个进汞区间孔喉平均半径,μm;ri:第i个进汞压力对应的孔喉半径,μm;ri+1:第i+1个进汞压力对应的孔喉半径,μm;n:进汞压力测试总点数,计算结果如表3所示。
②计算不同进汞区间薄膜束缚水饱和度
建立了薄模束缚水同心毛细管模型,如图3所示,设区间薄膜束缚水饱和度为Swni,其计算公式为:
Figure BDA0001994842690000072
式中,rci:第i个进汞区间孔喉平均半径,μm;rn:束缚水薄膜厚度,μm;Swni:第i个进汞区间薄膜束缚水饱和度,%;ΔSi:第i个进汞区间汞饱和度增量,%;n:进汞压力测试总点数,计算结果如表3所示。
4)确定岩心薄膜束缚水饱和度
设岩心薄膜束缚水饱和度为Swn,其计算公式为:
Figure BDA0001994842690000073
式中,Swni:第i个进汞区间薄膜束缚水饱和度,%;Swn:岩心薄膜束缚水饱和度,%,n:进汞压力测试总点数,计算结果如表3所示。
表3薄膜束缚水饱和度计算结果表
Figure BDA0001994842690000074
Figure BDA0001994842690000081
以上描述仅为本申请的实施例以及对所运用技术原理的说明。本领域技术人员应当理解,本申请中所涉及的发明范围,并不限于上述技术特征的特定组合而成的技术方案,同时也应涵盖在不脱离所述发明构思的情况下,由上述技术特征或其等同特征进行任意组合而形成的其它技术方案。例如上述特征与本申请中公开的(但不限于)具有类似功能的技术特征进行互相替换而形成的技术方案。

Claims (7)

1.一种岩心薄膜束缚水饱和度的测定方法,包括如下步骤:
1)对岩心进行压汞实验,得到不同进汞压力下的进汞饱和度,进而得到岩心的孔喉半径和孔喉分布频率;
2)根据所述进汞压力得到岩心的薄膜束缚水厚度;
3)根据所述孔喉半径得到岩心的不同进汞区间孔喉平均半径;
4)根据所述孔喉分布频率、所述薄膜束缚水厚度和所述不同进汞区间孔喉平均半径得到岩心的不同进汞区间薄膜束缚水饱和度,进而得到岩心的薄膜束缚水饱和度。
2.根据权利要求1所述的测定方法,其特征在于:步骤1)中,根据所述进汞压力,结合式(1)得到所述孔喉半径;
Figure FDA0002982246250000011
式(1)中,ri表示第i个进汞压力对应的孔喉半径,μm;σ表示汞表面张力,为0.48N/m;θ表示汞润湿角,为140°;Pi表示第i个测试点进汞压力,MPa;i=1,2,……,n-1,n表示进汞压力测试总点数。
3.根据权利要求1或2所述的测定方法,其特征在于:步骤1)中,根据所述进汞饱和度,结合式(2)得到所述孔喉分布频率;
fi=ΔSi=Si+1-Si (2)
式(2)中,fi表示第i个区间孔喉分布频率,%;ΔSi表示第i个区间汞饱和度增量,%;Si+1表示第i+1个测试点进汞压力对应的汞饱和度,%;Si表示第i个测试点进汞压力对应的汞饱和度,%;i=1,2,……,n-1,n表示进汞压力测试总点数。
4.根据权利要求3所述的测定方法,其特征在于:步骤2)中,根据式(3)得到所述薄膜束缚水厚度;
Figure FDA0002982246250000012
式(3)中,Pn表示第n个测试点进汞压力,MPa;rn表示束缚水薄膜厚度,μm;σ表示汞表面张力,为0.48N/m;θ表示汞润湿角,为140°。
5.根据权利要求4所述的测定方法,其特征在于:步骤3)中,根据式(4)得到所述不同进汞区间孔喉平均半径;
Figure FDA0002982246250000013
式(4)中,rci表示第i个进汞区间孔喉平均半径,μm;ri表示第i个测试点进汞压力对应的孔喉半径,μm;ri+1表示第i+1个测试点进汞压力对应的孔喉半径,μm;i=1,2,……,n-1,n表示进汞压力测试总点数。
6.根据权利要求5所述的测定方法,其特征在于:步骤4)中,根据式(5)得到所述不同进汞区间薄膜束缚水饱和度;
Figure FDA0002982246250000021
式(5)中,Swni表示第i个进汞区间薄膜束缚水饱和度,%;rci表示第i个进汞区间孔喉平均半径,μm;rn表示束缚水薄膜厚度,μm;ΔSi表示第i个进汞区间汞饱和度增量,即第i个区间孔喉分布频率,%;i=1,2,……,n-1,n表示进汞压力测试总点数。
7.根据权利要求6所述的测定方法,其特征在于:步骤4)中,根据式(6)得到所述薄膜束缚水饱和度;
Figure FDA0002982246250000022
式(6)中,Swn:岩心薄膜束缚水饱和度,%,Swni表示第i个进汞区间薄膜束缚水饱和度,%;n表示进汞压力测试总点数。
CN201910192711.6A 2019-03-14 2019-03-14 一种基于压汞曲线的岩心薄膜束缚水饱和度的测定方法 Expired - Fee Related CN111693425B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910192711.6A CN111693425B (zh) 2019-03-14 2019-03-14 一种基于压汞曲线的岩心薄膜束缚水饱和度的测定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910192711.6A CN111693425B (zh) 2019-03-14 2019-03-14 一种基于压汞曲线的岩心薄膜束缚水饱和度的测定方法

Publications (2)

Publication Number Publication Date
CN111693425A CN111693425A (zh) 2020-09-22
CN111693425B true CN111693425B (zh) 2021-05-11

Family

ID=72475089

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910192711.6A Expired - Fee Related CN111693425B (zh) 2019-03-14 2019-03-14 一种基于压汞曲线的岩心薄膜束缚水饱和度的测定方法

Country Status (1)

Country Link
CN (1) CN111693425B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114112841B (zh) * 2021-11-19 2023-06-09 中国石油大学(华东) 一种致密砂岩储层束缚水饱和度的计算方法

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2341372A1 (en) * 2009-12-16 2011-07-06 BP Exploration Operating Company Limited Method for measuring rock wettability
CN101799559B (zh) * 2010-02-09 2011-12-07 康志勇 基于砂岩油藏油层电阻率的地层数据处理方法
CN103278436B (zh) * 2013-02-01 2015-02-04 西安石油大学 特低渗透双重介质砂岩油藏微观孔隙结构的定量表征方法
CN104458525B (zh) * 2013-09-13 2017-04-05 中国石油天然气股份有限公司 一种常规压汞曲线表征微观孔隙喉道特征的采集处理方法
CN103900942B (zh) * 2014-01-10 2016-06-01 中国石油大学(华东) 基于高压压汞分析的储层微观孔喉参数的连续表征方法
US10495774B2 (en) * 2014-01-13 2019-12-03 Schlumberger Technology Corporation Method for estimating irreducible water saturation from mercury injection capillary pressure
CN104200105A (zh) * 2014-09-05 2014-12-10 中国石油大学(华东) 确定致密砂岩气充注物性下限的方法
CN107944671B (zh) * 2017-11-07 2021-11-02 中国石油天然气股份有限公司 一种确定地层束缚水饱和度的方法

Also Published As

Publication number Publication date
CN111693425A (zh) 2020-09-22

Similar Documents

Publication Publication Date Title
CN110346258B (zh) 一种致密岩石油相相对渗透率测定的方法
CN107063967B (zh) 一种基于微观孔隙通道的饱和/非饱和渗透系数预测方法
CN107525746B (zh) 一种表征致密油藏低速非达西渗流特征的方法及装置
CN110598167B (zh) 低渗透油藏油水相对渗透率实验数据的处理方法
Muskat et al. Flow of gas through porous materials
Wan et al. Chloride content and pH value in the pore solution of concrete under carbonation
CN114720331B (zh) 岩土体迂曲度幂律浆液时变性的柱半球渗透半径确定方法
CN110516322B (zh) 一种不同水力梯度下黏土饱和非线性渗透系数的预测方法
CN111693425B (zh) 一种基于压汞曲线的岩心薄膜束缚水饱和度的测定方法
CN108444890B (zh) 测试中、高渗岩心液测渗透率的非稳态滴定装置及方法
Lu et al. Constant flow method for concurrently measuring soil-water characteristic curve and hydraulic conductivity function
CN115618643B (zh) 宾汉浆液时水效应岩土体迂曲度柱半球渗透半径确定方法
CN106777515B (zh) 一种基于岩心实验资料分析致密气井产能的方法
CN108133086A (zh) 一种应力敏感储层中产水压裂气井裂缝半长反演方法
Yanjie et al. Lower limit of tight oil flowing porosity: Application of high-pressure mercury intrusion in the fourth Member of Cretaceous Quantou Formation in southern Songliao Basin, NE China
CN109520891A (zh) 过滤材料接触角的测定方法
CN116383573B (zh) 一种基于多区相变传质渗流耦合的凝析气产能评价方法
CN109932296B (zh) 一种定量表征贾敏效应动态变化的方法
CN111208053A (zh) 一种基于恒压压汞实验的储层微毛细管孔隙度确定方法
CN114112841B (zh) 一种致密砂岩储层束缚水饱和度的计算方法
CN115083543A (zh) 一种超大掺量矿物掺合料水泥石毛细吸水性能的预测方法
CN115559715A (zh) 一种超高压低渗透气藏储层产水评价方法
CN207457017U (zh) 一种用于测试多孔工程材料孔隙分布特征的试验装置
CN109490162A (zh) 考虑应力敏感条件下的致密油藏毛细管压力计算方法
CN113945497B (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: 20210511

CF01 Termination of patent right due to non-payment of annual fee