CN114508331B - 一种基于霍尔曲线导数确定聚合物注入能力新方法 - Google Patents

一种基于霍尔曲线导数确定聚合物注入能力新方法 Download PDF

Info

Publication number
CN114508331B
CN114508331B CN202210063243.4A CN202210063243A CN114508331B CN 114508331 B CN114508331 B CN 114508331B CN 202210063243 A CN202210063243 A CN 202210063243A CN 114508331 B CN114508331 B CN 114508331B
Authority
CN
China
Prior art keywords
polymer
hall
injection
derivative
curve
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
Application number
CN202210063243.4A
Other languages
English (en)
Other versions
CN114508331A (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.)
Northeast Petroleum University
Original Assignee
Northeast Petroleum University
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 Northeast Petroleum University filed Critical Northeast Petroleum University
Priority to CN202210063243.4A priority Critical patent/CN114508331B/zh
Publication of CN114508331A publication Critical patent/CN114508331A/zh
Application granted granted Critical
Publication of CN114508331B publication Critical patent/CN114508331B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • 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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • 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—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/02—Agriculture; Fishing; Forestry; Mining
    • 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
    • Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
    • Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Mining & Mineral Resources (AREA)
  • Geology (AREA)
  • Theoretical Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Environmental & Geological Engineering (AREA)
  • Fluid Mechanics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • General Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Tourism & Hospitality (AREA)
  • Animal Husbandry (AREA)
  • Economics (AREA)
  • General Health & Medical Sciences (AREA)
  • Human Resources & Organizations (AREA)
  • Marketing (AREA)
  • Primary Health Care (AREA)
  • Strategic Management (AREA)
  • Marine Sciences & Fisheries (AREA)
  • General Business, Economics & Management (AREA)
  • Agronomy & Crop Science (AREA)
  • Geophysics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
  • Injection Moulding Of Plastics Or The Like (AREA)

Abstract

一种基于霍尔曲线导数确定聚合物注入能力新方法，该方法包括下列步骤：1、筛选能够反映研究区块储层物性参数的代表性取芯岩心，开展岩心不同渗透率级别、不同聚合物浓度条件下注入能力室内实验；2、基于步骤1聚合物注入能力实验数据，绘制不同渗透率级别、不同聚合物浓度条件下霍尔曲线和霍尔导数曲线；3、绘制不同聚合物浓度条件下霍尔导数随储层岩心渗透率变化关系曲线，根据聚驱阶段霍尔导数发生突变位置，确定临界点。该方法具有普适性，能够评价给出不同类型油藏、不同类型聚合物的注入能力，解决了传统经验方法的不足，为聚合物注入能力评价提供了新思路，提供了理论依据。

Description

一种基于霍尔曲线导数确定聚合物注入能力新方法
技术领域：
本发明涉及石油工程中化学驱提高采收率领域，具体涉及一种基于霍尔曲线导数确定聚合物注入能力新方法。
背景技术：
化学剂注入能力是保证化学驱提高原油采收率效果的前提，因此对化学剂注入能力的评价工作十分重要。目前通常采用注聚压力、粘度保留率、封堵能力、阻力系数和残余阻力系数等参数来评价化学驱过程中聚合物的注入能力。
目前油田开发现场工程技术人员在进行聚合物驱油作业时，通常利用残余阻力系数经验值(5和10)来评价聚合物注入能力。当残余阻力系数超过10时，认为聚合物注入压力高、注入能力差(专利：一种低渗透油藏聚合物注入能力评价新方法)。根据实际现场聚合物注入难易程度的经验确定残余阻力系数临界点，缺少理论依据，经验方法具有偶然性和误差大的缺点，同时不能保证该经验值适用于所有储层。
发明内容：
为了解决背景技术中存在的技术问题，本发明提供一种基于霍尔曲线导数确定聚合物注入能力新方法。该方法具有普适性，能够评价给出不同类型油藏、不同类型聚合物的注入能力，解决了传统经验方法的不足，为聚合物注入能力评价提供了新思路，提供了理论依据。
本发明的技术方案是：一种基于霍尔曲线导数确定聚合物注入能力新方法，该方法包括下列步骤：
1、筛选能够反映研究区块储层物性参数的代表性取芯岩心，开展岩心不同渗透率级别、不同聚合物浓度条件下注入能力室内实验，为了确保实验准确性，取芯岩心数量≥4块，每个渗透率级别下各1块岩心，聚合物浓度≥4种；
具体地，聚合物注入能力评价室内实验流程如下：
(1)利用覆压孔渗仪测定取芯岩心孔隙度和渗透率参数，筛选实验用岩心并对岩心进行编号；
(2)将取芯岩心装入岩心夹持器抽真空后、饱和地层水，并测定其孔隙度与渗透率；
(3)根据聚驱方案要求配制聚合物溶液并测定其粘度，剪切至聚合物溶液粘度降半为止并装入中间容器备用；
(4)以恒定速度水驱至压力平稳、注入步骤(3)配制的聚合物体系至压力平稳、后续水驱至压力平稳，每隔30min记录注入量、累积注入量和注入压力；
(5)整理实验数据并绘制阻力系数/残余阻力系数随注入PV数变化关系曲线；
2、基于步骤1聚合物注入能力实验数据，绘制不同渗透率级别、不同聚合物浓度条件下霍尔曲线和霍尔导数曲线。
具体地，霍尔曲线和霍尔导数曲线绘制过程如下：
(1)基于聚合物注入能力实验数据，利用公式(1)计算第i个实验数据点对应的霍尔积分：
Figure BDA0003479043980000021
其中，Ih,i—第i个实验数据点对应的霍尔积分(i≥2)，MPa·min；△ti—第i-1个与第i个数据记录点时间间隔，min；Pi-1—第i-1个数据点对应的注入压力，MPa；Pi—第i个数据点对应的注入压力，MPa；n—实验数据总记录点；
(2)基于步骤(1)计算结果，绘制霍尔积分随累积注入量变化关系曲线，即霍尔曲线；
(3)利用公式(2)，计算霍尔导数：
Figure BDA0003479043980000022
其中，f(Ih,i)—第i个实验数据点对应的霍尔导数(i≥2)，MPa·min/mL；Ih,i-1—第i-1个实验数据点对应的霍尔积分，MPa·min；Qi-1—第i-1个数据点对应的累积注入量，mL；Qi—第i个数据点对应的累积注入量，mL；
(4)基于步骤(3)计算结果，绘制霍尔导数随累积注入量变化关系曲线，即霍尔导数曲线。
3、绘制不同聚合物浓度条件下霍尔导数随储层岩心渗透率变化关系曲线，根据聚驱阶段霍尔导数发生突变位置，确定拐点，由拐点处低一级别渗透率对应的残余阻力系数值作为该浓度聚合物注入能力突变点，取不同浓度聚合物注入能力突变点对应的残余阻力系数平均值，作为评价聚合物注入能力的临界点。若聚合物残余阻力系数值大于临界点值，则表明该聚合物注入困难；若聚合物残余阻力系数值小于临界点值，则表明该聚合物能够顺利注入。
本发明具有如下有益效果：通过本申请的方法建立如下注入能力评价标准：残余阻力系数＞注入能力临界点，代表注入困难；残余阻力系数＜注入能力临界点，代表顺利注入。该方法具有普适性，能够评价给出不同类型油藏、不同类型聚合物的注入能力，解决了传统经验方法取临界点为10的不足，为聚合物注入能力评价提供了新思路，提供了理论依据。
附图说明：
图1是抗盐聚合物粘度-浓度曲线
图2是阻力系数/残余阻力系数随注入PV数变化关系(聚合物浓度3000mg/L)；
图3是阻力系数/残余阻力系数随注入PV数变化关系(聚合物浓度2300mg/L)；
图4是阻力系数/残余阻力系数随注入PV数变化关系(聚合物浓度1500mg/L)；
图5是阻力系数/残余阻力系数随注入PV数变化关系(聚合物浓度800mg/L)；
图6是不同渗透率岩心条件下霍尔曲线(聚合物浓度3000mg/L)；
图7是不同渗透率岩心条件下霍尔曲线(聚合物浓度2300mg/L)；
图8是不同渗透率岩心条件下霍尔曲线(聚合物浓度1500mg/L)；
图9是不同渗透率岩心条件下霍尔曲线(聚合物浓度800mg/L)；
图10是霍尔导数随渗透率变化关系(聚合物浓度3000mg/L)；
图11是霍尔导数随渗透率变化关系(聚合物浓度2300mg/L)；
图12是霍尔导数随渗透率变化关系(聚合物浓度1500mg/L)；
图13是霍尔导数随渗透率变化关系(聚合物浓度800mg/L)；
具体实施方式：
下面结合具体实验对本发明作进一步说明：
一、聚合物注入能力评价室内实验
1、筛选与研究区块储层物性参数相匹配的天然露头方岩心(尺寸：4.5×4.5×30cm左右)，确保方岩心气测渗透率在50-100、100-150、150-200、200-250、250-300、300-350mD级别(基于研究区块储层物性渗透率变化范围确定)。实验用材料及相关参数如下：实验用模拟地层水(矿化度19534mg/L，钙镁离子含量840mg/L)；实验用抗盐聚合物浓度分别为800mg/L、1500mg/L、2300mg/L和3000mg/L；驱替速度0.3mL/min；实验温度85℃。
2、实验方案设计：
开展岩心渗透率为50-100、100-150、150-200、200-250、250-300、300-350mD级别条件下4种浓度(800mg/L、1500mg/L、2300mg/L、3000mg/L)聚合物溶液注入能力评价实验，具体实验方案设计参见表1。
表1抗盐聚合物注入能力实验方案设计
Figure BDA0003479043980000041
3、实验步骤：
(1)测定天然露头方岩心孔渗参数，筛选实验用岩心并对岩心进行编号；
(2)将天然露头方岩心装入岩心夹持器抽真空，饱和水后测定其孔隙度与渗透率，见表2；
(3)配制实验用抗盐聚合物溶液(800mg/L、1500mg/L、2300mg/L、3000mg/L)并测定粘度-浓度曲线，粘度-浓度曲线见图1，剪切至聚合物溶液粘度降半为止并装入中间容器备用；
(4)以0.3mL/min注入速度水驱至压力平稳，之后以0.3mL/min注入聚合物体系(浓度800mg/L)至压力平稳，再之后以0.3mL/min注入速度开展后续水驱至压力平稳，每隔30min记录注入量、累积注入量和注入压力。浓度为1500mg/L、2300mg/L、3000mg/L聚合物注入程序与800mg/L聚合物相同；
(5)整理实验数据并绘制4种聚合物浓度在6个不同渗透率级别条件下阻力系数/残余阻力系数随注入PV数变化关系曲线，见图2、图3、图4及图5。
表2实验用方岩心基本物性参数
Figure BDA0003479043980000051
Figure BDA0003479043980000061
4、实验结果：
基于24组抗盐聚合物注入能力实验记录数据，统计给出了不同聚合物浓度在不同渗透率级别条件下的阻力系数/残余阻力系数，见表3。
表3阻力系数/残余阻力系数统计结果
Figure BDA0003479043980000062
二、霍尔曲线绘制和聚驱阶段霍尔导数确定
霍尔曲线代表累积注入量与霍尔积分之间的关系。当注入流体不同时，霍尔曲线上会出现不同的直线段，利用不同注入阶段稳定时期曲线斜率即可评价注入流体介质的有效性。
基于聚合物注入能力实验数据和霍尔积分计算结果，绘制了抗盐聚合物不同渗透率、不同浓度条件下的霍尔曲线，见图6-图9。基于霍尔曲线和霍尔导数计算结果，给出了不同浓度聚合物在不同渗透率条件下聚驱阶段的霍尔导数，阻力系数/残余阻力系数/聚驱阶段霍尔导数对比结果见表4。
表4阻力系数/残余阻力系数/聚驱阶段霍尔导数对比结果
Figure BDA0003479043980000071
三、霍尔导数随储层岩心渗透率变化关系曲线绘制
基于聚驱阶段霍尔导数计算结果，绘制了不同聚合物浓度条件下霍尔导数随储层岩心渗透率变化关系曲线，见图10-图13。根据抗盐聚合物聚驱阶段霍尔导数发生突变位置，确定不同聚合物浓度条件下对应的霍尔导数拐点。根据拐点对应的储层渗透率和聚合物浓度，确定拐点位置对应的残余阻力系数值。基于上述原则，该抗盐聚合物在浓度800mg/L、1500mg/L、2300mg/L和3000mg/L条件下拐点位置处对应的渗透率分别为145mD、203mD、208mD和209mD。由拐点处低一级别渗透率对应的残余阻力系数值作为该浓度聚合物注入能力突变点，取不同浓度聚合物注入能力突变点对应的残余阻力系数平均值，作为评价聚合物注入能力的临界点。该抗盐聚合物在浓度800mg/L、1500mg/L、2300mg/L和3000mg/L条件下突变点对应的残余阻力系数值分别为8.14、6.58、7.52和8.12，则该抗盐聚合物注入能力的临界点为7.59。由此可知：若聚合物残余阻力系数值大于临界点值7.59，则表明该聚合物注入困难；若聚合物残余阻力系数值小于临界点值7.59，则表明该聚合物能够顺利注入。

Claims (2)

1.一种基于霍尔曲线导数确定聚合物注入能力新方法，该方法包括下列步骤：
筛选能够反映研究区块储层物性参数的代表性取芯岩心，开展岩心不同渗透率级别、不同聚合物浓度条件下注入能力室内实验，为了确保实验准确性，取芯岩心数量≥4块，每个渗透率级别下各1块岩心，聚合物浓度≥4种；
具体地，聚合物注入能力评价室内实验流程如下：
(1)利用覆压孔渗仪测定取芯岩心孔隙度和渗透率参数，筛选实验用岩心并对岩心进行编号；
(2)将取芯岩心装入岩心夹持器抽真空后、饱和地层水，并测定其孔隙度与渗透率；
(3)根据聚驱方案要求配制聚合物溶液并测定其粘度，剪切至聚合物溶液粘度降半为止并装入中间容器备用；
(4)以恒定速度水驱至压力平稳、注入步骤(3)配制的聚合物体系至压力平稳、后续水驱至压力平稳，每隔30min记录注入量、累积注入量和注入压力；
(5)整理实验数据并绘制阻力系数/残余阻力系数随注入PV数变化关系曲线；
基于步骤(1)聚合物注入能力实验数据，绘制不同渗透率级别、不同聚合物浓度条件下霍尔曲线和霍尔导数曲线，
具体地，霍尔曲线和霍尔导数曲线绘制过程如下：
(1)基于聚合物注入能力实验数据，利用公式(1)计算第i个实验数据点对应的霍尔积分：
Figure FDA0003722701710000011
其中，Ih,i—第i个实验数据点对应的霍尔积分i≥2，MPa·min；△ti—第i-1个与第i个数据记录点时间间隔，min；Pi-1—第i-1个数据点对应的注入压力，MPa；Pi—第i个数据点对应的注入压力，MPa；n—实验数据总记录点；
(2)基于步骤(1)计算结果，绘制霍尔积分随累积注入量变化关系曲线，即霍尔曲线；
(3)利用公式(2)，计算霍尔导数：
Figure FDA0003722701710000012
其中，f(Ih,i)—第i个实验数据点对应的霍尔导数i≥2，MPa·min/mL；Ih,i-1—第i-1个实验数据点对应的霍尔积分，MPa·min；Qi-1—第i-1个数据点对应的累积注入量，mL；Qi—第i个数据点对应的累积注入量，mL；
(4)基于步骤(3)计算结果，绘制霍尔导数随累积注入量变化关系曲线，即霍尔导数曲线。
2.根据权利要求1所述的基于霍尔曲线导数确定聚合物注入能力新方法，其特征在于：绘制不同聚合物浓度条件下霍尔导数随储层岩心渗透率变化关系曲线，根据聚驱阶段霍尔导数发生突变位置，确定拐点，由拐点处低一级别渗透率对应的残余阻力系数值作为该浓度聚合物注入能力突变点，取不同浓度聚合物注入能力突变点对应的残余阻力系数平均值，作为评价聚合物注入能力的临界点，若聚合物残余阻力系数值大于临界点值，则表明该聚合物注入困难；若聚合物残余阻力系数值小于临界点值，则表明该聚合物能够顺利注入。
CN202210063243.4A 2022-01-20 2022-01-20 一种基于霍尔曲线导数确定聚合物注入能力新方法 Active CN114508331B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210063243.4A CN114508331B (zh) 2022-01-20 2022-01-20 一种基于霍尔曲线导数确定聚合物注入能力新方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210063243.4A CN114508331B (zh) 2022-01-20 2022-01-20 一种基于霍尔曲线导数确定聚合物注入能力新方法

Publications (2)

Publication Number Publication Date
CN114508331A CN114508331A (zh) 2022-05-17
CN114508331B true CN114508331B (zh) 2022-08-12

Family

ID=81550323

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210063243.4A Active CN114508331B (zh) 2022-01-20 2022-01-20 一种基于霍尔曲线导数确定聚合物注入能力新方法

Country Status (1)

Country Link
CN (1) CN114508331B (zh)

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9902894B2 (en) * 2014-10-31 2018-02-27 Chevron U.S.A. Inc. Polymer compositions
US10619087B2 (en) * 2015-12-08 2020-04-14 Chevron U.S.A. Inc. Methods for hydrocarbon recovery
CN105715241B (zh) * 2016-01-15 2017-11-17 中国石油大学(华东) 一种聚合物驱相对渗透率曲线的测量方法
CN109184644B (zh) * 2018-09-28 2021-06-15 中国海洋石油集团有限公司 一种考虑聚合物非牛顿性和渗流附加阻力的早期注聚效果评价方法
CN112145139B (zh) * 2020-10-26 2022-02-01 西南石油大学 一种原油黏度合理映射聚合物驱使用浓度的方法

Also Published As

Publication number Publication date
CN114508331A (zh) 2022-05-17

Similar Documents

Publication Publication Date Title
CN108362614B (zh) 测量页岩油二氧化碳吞吐过程中扩散系数的装置及方法
CN109883917B (zh) 一种模拟基质-裂缝静动态渗吸排油的实验装置及方法
CN106596371B (zh) 饱和凝析气藏衰竭式开发近井区反凝析伤害实验评价方法
WO2016061854A1 (zh) 一种可渗透边界层天然气水合物开采模拟实验装置
CN113027418B (zh) 油田注水井污染解堵实验测试装置及测试效果评价方法
CN106884635A (zh) 一种低、特低渗透油藏co2驱最小混相压力的测定方法
CN106837269A (zh) 一种低、特低渗透油藏co2驱近混相压力区域确定方法
CN104865165B (zh) 全过程通量守恒的现场土壤溶质运移弥散系数测定的方法
CN107402286B (zh) 一种聚合物微球岩心内缓膨性能测试实验方法
CN106988711A (zh) 一种提高强纵向非均质性油藏驱油效果的方法
CN108344853A (zh) 平面非均质性定容干气气藏直井绝对无阻流量测试方法
CN104879102A (zh) 一种薄互层碳酸盐岩底水油藏co2吞吐实验测试方法
Liu et al. Experimental investigation into formation/dissociation characteristics of methane hydrate in consolidated sediments with resistance measurement
CN112081560A (zh) 一种海上深层高温超压气藏开发方法
CN105651665A (zh) 一种钻完井液对岩心油水渗透率影响评价方法
CN113075108B (zh) 一种考虑束缚水饱和度的岩心多次应力敏感测试方法
Liwei et al. Experimental study on gelling property and plugging effect of inorganic gel system (OMGL)
CN111878075A (zh) 一种倾斜油藏气水协同注入分区采出程度的测试方法
CN111472736A (zh) 海上油田组合调驱优化设计方法及装置
CN107402290A (zh) 研究co2驱过程中原油中沥青质沉淀对储层物性及开发效果影响的实验方法
CN110879196B (zh) 富含油凝析气藏油水相渗测试方法
CN114508331B (zh) 一种基于霍尔曲线导数确定聚合物注入能力新方法
CN114278266B (zh) 一种确定就地条件下酸蚀填砂裂缝有效缝长的方法
CN110715879A (zh) 基于气水分布的高演化页岩储层微孔隙吸附气量评价方法
CN116201538B (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