CN109002565B - 压裂液临界返排速度计算方法 - Google Patents

压裂液临界返排速度计算方法 Download PDF

Info

Publication number
CN109002565B
CN109002565B CN201710416025.3A CN201710416025A CN109002565B CN 109002565 B CN109002565 B CN 109002565B CN 201710416025 A CN201710416025 A CN 201710416025A CN 109002565 B CN109002565 B CN 109002565B
Authority
CN
China
Prior art keywords
proppant
settling
flowback
fracturing fluid
fracture
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
CN201710416025.3A
Other languages
English (en)
Other versions
CN109002565A (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 Petroleum and Chemical Corp
Sinopec Research Institute of Petroleum Engineering Shengli Co
Original Assignee
China Petroleum and Chemical Corp
Sinopec Research Institute of Petroleum Engineering Shengli Co
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 Petroleum and Chemical Corp, Sinopec Research Institute of Petroleum Engineering Shengli Co filed Critical China Petroleum and Chemical Corp
Priority to CN201710416025.3A priority Critical patent/CN109002565B/zh
Publication of CN109002565A publication Critical patent/CN109002565A/zh
Application granted granted Critical
Publication of CN109002565B publication Critical patent/CN109002565B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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
    • 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/25—Methods for stimulating production
    • E21B43/26—Methods for stimulating production by forming crevices or fractures
    • E21B43/267—Methods for stimulating production by forming crevices or fractures reinforcing fractures by propping

Landscapes

  • Engineering & Computer Science (AREA)
  • Mining & Mineral Resources (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Business, Economics & Management (AREA)
  • Physics & Mathematics (AREA)
  • Geology (AREA)
  • Economics (AREA)
  • General Health & Medical Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Agronomy & Crop Science (AREA)
  • Animal Husbandry (AREA)
  • Marine Sciences & Fisheries (AREA)
  • Environmental & Geological Engineering (AREA)
  • Health & Medical Sciences (AREA)
  • Fluid Mechanics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Human Resources & Organizations (AREA)
  • Marketing (AREA)
  • Primary Health Care (AREA)
  • Strategic Management (AREA)
  • Tourism & Hospitality (AREA)
  • General Business, Economics & Management (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明提供一种压裂液临界返排速度计算方法，该压裂液临界返排速度计算方法包括：步骤1，支撑剂在破胶液中进行沉降速度实验，以确定支撑剂沉降速度；步骤2，将支撑剂在裂缝中的沉降模型简化；步骤3，将支撑剂在孔眼中的沉降模型简化；步骤4，推导临界返排速度计算公式。该压裂液临界返排速度计算方法排除了目前现场经验的不确定性对返排速度的影响，确定了返排速度的临界值，解决了由于返排速度过快导致的地层吐砂等问题。

Description

压裂液临界返排速度计算方法
技术领域
本发明涉及油田开采领域，特别是涉及到一种压裂液临界返排速度计算方法。
背景技术
压裂作为开采低渗透油气藏的有效手段，应用越来越广泛。压裂后的压裂液返排，目前主要有两种方法：
一是快速返排：为了减少压裂液对地层的污染，压裂之后，采取措施立刻强制返排。但弊端非常明显，常常会出现压裂液没有完全破胶就携带支撑剂返排出来，导致堵塞油嘴。
二是根据压力控制返排：待压裂液破胶后，根据压力大小，更换不同尺寸的油嘴，控制放喷返排。
上述第二种方法，目前是压裂液返排的主流措施。但其依据主要是现场经验，因此常会出现裂缝吐砂，导致砂埋管柱、砂卡封隔器等重大井下事故。为此我们发明了一种新的压裂液临界返排速度计算方法，解决了以上技术问题。
发明内容
本发明的目的是提供一种根据支撑剂沉降实验和模型合理简化，给出了压裂液临界返排速度的计算方法和现场操作方式的压裂液临界返排速度计算方法。
本发明的目的可通过如下技术措施来实现：压裂液临界返排速度计算方法，该压裂液临界返排速度计算方法包括：步骤1，支撑剂在破胶液中进行沉降速度实验，以确定支撑剂沉降速度；步骤2，将支撑剂在裂缝中的沉降模型简化；步骤3，将支撑剂在孔眼中的沉降模型简化；步骤4，推导临界返排速度计算公式。
本发明的目的还可通过如下技术措施来实现：
在步骤1中，裂缝内支撑剂的运动，在垂直方向做自由沉降运动，在水平方向由于受压裂液的滞带作用而运动，因此可视为沉降作用与压裂液的滞带的合成。
在步骤1中，不考虑砂浓度和裂缝壁面的影响，单粒支撑剂在水平方向的速度与压裂液的水平流速一样。
在步骤1中，在计算垂直方向的沉降速度时，针对油田目前施工使用的支撑剂，进行压裂液不同粘度下的支撑剂沉降实验，得到支撑剂沉降实验的曲线图，根据该支撑剂沉降实验的曲线图，获得不同返排液粘度下的支撑剂沉降速度。
在步骤1中，油田目前施工使用的支撑剂为20/40目支撑剂，粒径规格为0.425-0.85mm，闭合压力为69MPa，体积密度为1.88g/cm3，球度≥0.8，圆度≥0.8，粒径规格内支撑剂质量所占比例≥90％。
在步骤1中，对油田目前施工使用的支撑剂，进行压裂液不同粘度下的支撑剂沉降实验，获得如下数据：返排液粘度为10mp.s时，支撑剂沉降速度为0.4cm/s；返排液粘度为8mp.s时，支撑剂沉降速度为1cm/s；返排液粘度为6mp.s时，支撑剂沉降速度为4cm/s；返排液粘度为3mp.s时，支撑剂沉降速度为12cm/s；返排液粘度为1mp.s时，支撑剂沉降速度为17cm/s；根据上述数据绘制该支撑剂沉降实验的曲线图。
在步骤2中，支撑剂回流涉及到裂缝形态和支撑剂在裂缝中的沉降规律，裂缝形态简化为：裂缝为垂直缝，KGD或PKN模型，支撑剂在裂缝内均匀分布；支撑剂由于粒径小，受裂缝壁面和支撑剂之间相互影响，在裂缝中的沉降视为小雷诺数运动。
在步骤2中，对裂缝模型和支撑剂沉降模型进一步简化，简化的裂缝模型认为，压裂液在进入射孔孔眼时，会冲蚀扩大孔眼，裂缝在射孔孔眼端部启裂；压裂液返排期间的支撑剂的运动，简化为在孔眼内的运动。
在步骤3中，返排液可能从每一个孔眼流出，假设孔眼的作用都一样，把所有的孔眼进一步简化成一个大的孔眼。
在步骤3中，假设所有的返排液、支撑剂都从该大的简化孔眼中流出，忽略孔壁和支撑剂相互之间的影响，那么支撑剂在此简化孔眼中的流动简化为返排液的拖曳运动和自身的沉降运动。
在步骤4中，根据步骤2和步骤3中对模型的简化，考虑一个最极端情况，一个支撑剂颗粒随返排液从射孔孔眼的端部流出，只要此支撑剂颗粒在到达井筒之前落在最下端的一个射孔孔眼的距离内，那么就不会发生吐砂造成封隔器卡封事故，据此可以计算出返排速度。
在步骤4中，推导出的临界返排速度计算公式为：
Figure BDA0001312874480000021
其中：
V返排：返排速度，m3/h；
V沉：支撑剂沉降速度，m/s；
n：射孔有效孔数；
d：孔眼直径，m；
L：射孔穿深，m。
本发明中的压裂液临界返排速度计算方法，以支撑剂在压裂液破胶液中的沉降速度实验为基础，对孔眼和裂缝端部进行了合理的模型简化，由此推导出了压裂液临界返排速度计算公式。排除了目前现场经验的不确定性对返排速度的影响，确定了返排速度的临界值，解决了由于返排速度过快导致的地层吐砂等问题。
附图说明
图1为本发明的压裂液临界返排速度计算方法的一具体实施例的流程图；
图2为本发明的一具体实施例中支撑剂沉降实验的曲线图；
图3为本发明的一具体实施例中裂缝模型的示意图；
图4为本发明的一具体实施例中支撑剂运动模型的示意图；
图5为本发明的一具体实施例中简化裂缝模型的示意图；
图6为本发明的一具体实施例中简化孔眼的示意图；
图7为本发明的一具体实施例中支撑剂运动的示意图；
图8为本发明的一具体实施例中KGD模型的示意图；
图9为本发明的一具体实施例中PKN模型的示意图。
具体实施方式
为使本发明的上述和其他目的、特征和优点能更明显易懂，下文特举出较佳实施例，并配合附图所示，作详细说明如下。
如图1所示，图1为本发明的压裂液临界返排速度计算方法的流程图。
步骤101，支撑剂沉降实验。支撑剂在破胶液中的沉降速度实验，以确定支撑剂沉降速度.
裂缝内支撑剂的运动，在垂直方向做自由沉降运动，在水平方向由于受压裂液的滞带作用而运动，因此可视为沉降作用与压裂液的滞带的合成。不考虑砂浓度和裂缝壁面的影响，单粒支撑剂在水平方向的速度最终应和压裂液的水平流速一样。垂直方向的沉降速度可以根据公式计算。但公式涉及的参数较多，计算复杂。针对油田目前施工使用的支撑剂，主要是20/40目支撑剂，据此做了压裂液不同粘度下的支撑剂沉降实验。
表1支撑剂指标
Figure BDA0001312874480000031
Figure BDA0001312874480000041
表2支撑剂沉降实验
Figure BDA0001312874480000042
图2为本发明的一具体实施例中支撑剂沉降实验的曲线图，由图2看出，压裂液破胶越彻底，支撑剂沉降就越快，返排时受的拖曳力就越小，地层吐砂的可能性就越小。目前压后关井2小时，返排液粘度一般分布在3-8mp.s。
步骤102，支撑剂在裂缝中沉降模型简化。
支撑剂回流涉及到裂缝形态和支撑剂在裂缝中的沉降规律。裂缝形态一般简化为：裂缝为垂直缝，GDK或PKN模型，支撑剂在裂缝内均匀分布等；支撑剂由于粒径小，受裂缝壁面和支撑剂之间相互影响，在裂缝中的沉降视为小雷诺数运动，可以通过公式计算出来。
人工裂缝在延伸的过程中，涉及到流体力学、固体力学、断裂力学和热力学等。KGD和PKN模型是对人工裂缝进行简化的二维模型。如图8所示，KGD模型特点：假设每一水平截面独立作用，即假设裂缝面任一点处裂缝宽度沿垂向变化远比水平方向的变化慢。在缝高远大于缝长或者储积层边界产生完全滑移的条件下成立缝端区域起着很重要的作用，而缝内压力可以估算。如图9所示，PKN模型特点：假设每一垂向截面独立作用，即假设截面的压力是由高度控制的而非由缝长控制的。在缝长远大于缝高的条件下成立没有考虑断裂力学和缝端的影响，而主要考虑了缝内流体的流动以及相应的压力梯度的影响。
图3和图4为裂缝模型和支撑剂运动模型的示意图。图中，H：缝高；L(t):缝长；W(x,z,t):缝宽。上述模型在计算支撑剂返排沉降过程中，公式较复杂。为简化计算，对裂缝模型和支撑剂沉降模型进一步简化如下，如图5所示。
简化的裂缝模型认为，压裂液在进入射孔孔眼时，会冲蚀扩大孔眼，裂缝在射孔孔眼端部启裂。因此，压裂液返排期间的支撑剂的运动，可以简化为在孔眼内的运动。
步骤103，支撑剂在孔眼中沉降模型简化。
返排液可能从每一个孔眼流出，假设孔眼的作用都一样，这样可以把所有的孔眼进一步简化成一个大的孔眼，如图6所示。
假设所有的返排液、支撑剂都从此孔眼中流出，忽略孔壁和支撑剂相互之间的影响，那么支撑剂在此简化孔眼中的流动可简化为返排液的拖曳运动和自身的沉降运动，如图7所示。
步骤104，临界返排速度计算公式推导。
根据上述简化，考虑一个最极端的情况，就是裂缝的高度和射孔高度完全一致，也就是裂缝的缝高完全贯通了各射孔孔眼。那么一粒支撑剂从最上端的一个射孔孔眼的端部(相当于图7的A点)随返排液流出，在到达井筒之前，只要落在最下端的一个射孔孔眼的距离内(图7中BC相当于最下端一个射孔孔眼的长度)，那么就不会发生“吐砂”造成封隔器卡封事故。据此可以计算出返排速度：
Figure BDA0001312874480000051
其中：
V返排：返排速度，m3/h；
V沉：支撑剂沉降速度，m/s；
n：射孔有效孔数；
d：孔眼直径，m；
L：射孔穿深，m。
在应用本发明的一具体实施例中，包括了以下步骤:
1、压裂之后，按规定的时间关井；
2、关井结束，开始放喷返排，待井筒的液体返排干净之后，取返排液检测粘度；
3、根据公式(1)，计算返排速度；
4、实例：
某井射孔段长10米，有效孔数136孔，孔径13mm，射孔穿深0.7m，胍胶压裂液、陶粒支撑剂压裂，压后关井2小时放喷，检测返排液粘度为3.5mp.s，从图2可查出此时支撑剂沉降速度为0.011m/s，根据公式(1)可以计算出合理的临界返排速度为：
V_返排<3.3m3/h。

Claims (6)

1.压裂液临界返排速度计算方法，其特征在于，该压裂液临界返排速度计算方法包括：
步骤1，支撑剂在破胶液中进行沉降速度实验，以确定支撑剂沉降速度；
步骤2，将支撑剂在裂缝中的沉降模型简化；
步骤3，将支撑剂在孔眼中的沉降模型简化；
步骤4，推导临界返排速度计算公式；
在步骤2中，支撑剂回流涉及到裂缝形态和支撑剂在裂缝中的沉降规律，裂缝形态简化为：裂缝为垂直缝，KGD或PKN模型，支撑剂在裂缝内均匀分布；支撑剂由于粒径小，受裂缝壁面和支撑剂之间相互影响，在裂缝中的沉降视为小雷诺数运动；对裂缝模型和支撑剂沉降模型简化，简化的裂缝模型认为，压裂液在进入射孔孔眼时，会冲蚀扩大孔眼，裂缝在射孔孔眼端部启裂；压裂液返排期间的支撑剂的运动，简化为在孔眼内的运动；
在步骤3中，假设孔眼的作用都一样，把所有的孔眼简化成一个大的孔眼；假设所有的返排液、支撑剂都从该大的简化孔眼中流出，忽略孔壁和支撑剂相互之间的影响，那么支撑剂在此简化孔眼中的流动简化为返排液的拖曳运动和自身的沉降运动；
在步骤4中，根据步骤2和步骤3中对模型的简化，考虑一个最极端情况，一个支撑剂颗粒随返排液从射孔孔眼的端部流出，只要此支撑剂颗粒在到达井筒之前落在最下端的一个射孔孔眼的距离内，那么就不会发生吐砂造成封隔器卡封事故，据此计算出临界返排速度；临界返排速度计算公式为：
Figure DEST_PATH_IMAGE002
（1）
其中：
Figure DEST_PATH_IMAGE004
m3/h；
Figure DEST_PATH_IMAGE006
：支撑剂沉降速度，m/s；
n：射孔有效孔数；
d：孔眼直径，m；
L ：射孔穿深，m。
2.根据权利要求1所述的压裂液临界返排速度计算方法，其特征在于，在步骤1中，裂缝内支撑剂的运动，在垂直方向做自由沉降运动，在水平方向由于受压裂液的滞带作用而运动，因此视为沉降作用与压裂液的滞带的合成。
3.根据权利要求2所述的压裂液临界返排速度计算方法，其特征在于，在步骤1中，不考虑砂浓度和裂缝壁面的影响，单粒支撑剂在水平方向的速度与压裂液的水平流速一样。
4.根据权利要求2所述的压裂液临界返排速度计算方法，其特征在于，在步骤1中，在计算垂直方向的沉降速度时，针对油田目前施工使用的支撑剂，进行压裂液不同粘度下的支撑剂沉降实验，得到支撑剂沉降实验的曲线图，根据该支撑剂沉降实验的曲线图，获得不同返排液粘度下的支撑剂沉降速度。
5.根据权利要求4所述的压裂液临界返排速度计算方法，其特征在于，在步骤1中，油田目前施工使用的支撑剂为20/40目支撑剂，粒径规格为0.425-0.85mm，闭合压力为69MPa，体积密度为1.88g/cm3，球度≥0.8，圆度≥0.8，粒径规格内支撑剂质量所占比例≥90%。
6.根据权利要求4所述的压裂液临界返排速度计算方法，其特征在于，在步骤1中，对油田目前施工使用的支撑剂，进行压裂液不同粘度下的支撑剂沉降实验，获得如下数据：返排液粘度为10mp.s时，支撑剂沉降速度为0.4cm/s；返排液粘度为8mp.s时，支撑剂沉降速度为1cm/s；返排液粘度为6mp.s时，支撑剂沉降速度为4cm/s；返排液粘度为3mp.s时，支撑剂沉降速度为12cm/s；返排液粘度为1mp.s时，支撑剂沉降速度为17cm/s；根据上述数据绘制该支撑剂沉降实验的曲线图。
CN201710416025.3A 2017-06-05 2017-06-05 压裂液临界返排速度计算方法 Active CN109002565B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710416025.3A CN109002565B (zh) 2017-06-05 2017-06-05 压裂液临界返排速度计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710416025.3A CN109002565B (zh) 2017-06-05 2017-06-05 压裂液临界返排速度计算方法

Publications (2)

Publication Number Publication Date
CN109002565A CN109002565A (zh) 2018-12-14
CN109002565B true CN109002565B (zh) 2021-08-27

Family

ID=64573457

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710416025.3A Active CN109002565B (zh) 2017-06-05 2017-06-05 压裂液临界返排速度计算方法

Country Status (1)

Country Link
CN (1) CN109002565B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110219630A (zh) * 2019-06-04 2019-09-10 西南石油大学 一种裂缝性砂岩气藏压裂水平井的压裂液返排量计算方法
CN111879668B (zh) * 2020-04-30 2023-09-26 中国石油天然气股份有限公司 暂堵球性能确定方法和装置
CN112112609B (zh) * 2020-09-14 2022-04-05 西南石油大学 一种实现气藏压后返排油嘴尺寸调控的方法
CN113109221B (zh) * 2021-04-09 2022-07-19 中国石油大学(北京) 一种支撑剂沉降速率测量装置

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6732800B2 (en) * 2002-06-12 2004-05-11 Schlumberger Technology Corporation Method of completing a well in an unconsolidated formation
CN104329072A (zh) * 2014-10-30 2015-02-04 中国石油天然气股份有限公司 一种压裂液返排控制装置及方法
CN105574283A (zh) * 2015-12-24 2016-05-11 中国石油化工股份有限公司 一种压裂后放喷油嘴尺寸的确定方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6732800B2 (en) * 2002-06-12 2004-05-11 Schlumberger Technology Corporation Method of completing a well in an unconsolidated formation
CN104329072A (zh) * 2014-10-30 2015-02-04 中国石油天然气股份有限公司 一种压裂液返排控制装置及方法
CN105574283A (zh) * 2015-12-24 2016-05-11 中国石油化工股份有限公司 一种压裂后放喷油嘴尺寸的确定方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
支撑剂返排控制优化;吴亚红 等;《断块油气田》;20120925;第19卷(第5期);第662-665页 *
水力压裂返排优化设计研究;赵英杰;《中国优秀硕士学位论文全文数据库工程科技Ⅰ辑（月刊）》;20150615;正文第13，39页 *

Also Published As

Publication number Publication date
CN109002565A (zh) 2018-12-14

Similar Documents

Publication Publication Date Title
CN109002565B (zh) 压裂液临界返排速度计算方法
CN110134984B (zh) 一种页岩压裂过程中复杂裂缝扩展影响因素的分析方法
CN104508079A (zh) 改进水力裂缝网络的方法
CN112112609B (zh) 一种实现气藏压后返排油嘴尺寸调控的方法
US11613992B2 (en) Methods and systems for characterizing and/or monitoring wormhole regimes in matrix acidizing
CA2953618A1 (en) Compound cluster placement in fractures
CN105574283A (zh) 一种压裂后放喷油嘴尺寸的确定方法
Barree Potential Issues with Extreme Limited Entry in Horizontal Wells
WO2012143666A2 (en) System and method for improved propped fracture geometry for high permeability reservoirs
CN110173225A (zh) 页岩气水平井压裂后控砂返排工艺
AU2015256622B2 (en) Real time tool erosion prediction monitoring
CN109854194A (zh) 钻井液循环系统、减小钻井起钻抽吸压力的方法及设备
CN111810108B (zh) 一种页岩气水平井压后返排油嘴动态调整系统及方法
Zhang et al. Use of far-field diverters to mitigate parent-and infill-well-fracture interactions in shale formations
Brown et al. Factors influencing optimum ball sealer performance
WO2013090049A1 (en) Method and apparatus for modeling high solids content fluid fracturing
RU2547530C1 (ru) Способ разработки газонефтяных залежей
Barree Processes of Screenout Development and Avoidance
Tran et al. Correlating flowing time and condition for perforation plugging by suspended particles
Hwang et al. Filtration in frac packs and its impact on injector performance
US11215034B2 (en) Controlling redistribution of suspended particles in non-Newtonian fluids during stimulation treatments
CN116498287B (zh) 一种加砂压裂人工裂缝高度控制方法
Hardesty Improved Hydraulic Fracturing Perforation Efficiency Observed with Constant Entry Hole and Constant Penetration Perforating System
CN114607335A (zh) 确定体积压裂对应的压裂液流量方法
US20230383639A1 (en) Automatic real time screen-out mitigation

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