CN113554259B - 一种渗流阻力与驱替压力双场匹配优化评价方法 - Google Patents

一种渗流阻力与驱替压力双场匹配优化评价方法 Download PDF

Info

Publication number
CN113554259B
CN113554259B CN202010341083.6A CN202010341083A CN113554259B CN 113554259 B CN113554259 B CN 113554259B CN 202010341083 A CN202010341083 A CN 202010341083A CN 113554259 B CN113554259 B CN 113554259B
Authority
CN
China
Prior art keywords
ith
decimal
seepage resistance
seepage
parameters
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
CN202010341083.6A
Other languages
English (en)
Other versions
CN113554259A (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
Exploration and Development Research Institute of Sinopec Shengli Oilfield Co
Original Assignee
China Petroleum and Chemical Corp
Exploration and Development Research Institute of Sinopec Shengli Oilfield 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, Exploration and Development Research Institute of Sinopec Shengli Oilfield Co filed Critical China Petroleum and Chemical Corp
Priority to CN202010341083.6A priority Critical patent/CN113554259B/zh
Publication of CN113554259A publication Critical patent/CN113554259A/zh
Application granted granted Critical
Publication of CN113554259B publication Critical patent/CN113554259B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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
    • G06Q10/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0639Performance analysis of employees; Performance analysis of enterprise or organisation operations
    • G06Q10/06393Score-carding, benchmarking or key performance indicator [KPI] analysis
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/02Agriculture; Fishing; Forestry; Mining

Landscapes

  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Engineering & Computer Science (AREA)
  • Strategic Management (AREA)
  • Economics (AREA)
  • Theoretical Computer Science (AREA)
  • General Business, Economics & Management (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Educational Administration (AREA)
  • Development Economics (AREA)
  • General Physics & Mathematics (AREA)
  • Marketing (AREA)
  • Physics & Mathematics (AREA)
  • Tourism & Hospitality (AREA)
  • Agronomy & Crop Science (AREA)
  • Primary Health Care (AREA)
  • Mining & Mineral Resources (AREA)
  • Animal Husbandry (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Marine Sciences & Fisheries (AREA)
  • Health & Medical Sciences (AREA)
  • Game Theory and Decision Science (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明涉及油田开发技术领域,具体地涉及一种渗流阻力与驱替压力双场匹配优化评价方法,该方法的建立主要针对水淹水窜严重的高含水油藏的无效水循环区域,含油饱和度场与压力场不匹配的问题,根据流线分布状况判识无效循环带,并依据静态渗流屏障(断层、夹层、侧积层、岩性边界等)及无效水循环动态边界(水窜通道)划分局部流动体系,建立潜力参数、渗流阻力参数、驱动力参数及驱动力与渗流阻力适配参数,根据各种措施手段改变流动阻力的分布状况,驱替高饱和油,控制高饱和水,根据渗流阻力场与驱替压力场的适配性参数评价优选最佳方案。

Description

一种渗流阻力与驱替压力双场匹配优化评价方法
技术领域
本发明涉及油田开发技术领域,具体地涉及一种渗流阻力与驱替压力双场匹配优化评价方法。
背景技术
储层非均质、注采井距差异性以及边界的不规则性等原因导致局部区域含水差异大、无效水循环严重等问题,主要是由于局部渗流阻力与驱替压力不匹配,水窜通道渗流阻力小,但驱替压力大,导致无效水循环更为严重;而难动用区域渗流阻力大,却没有建立有效驱替,导致剩余油富集。
目前现场针对此问题主要开展定性方面的注采调配,没有对无效水循环等现象进行更深层次的剖析,也没有对措施后渗流阻力场及驱替压力场匹配关系进行定量化评价。因此,有必要发明一种渗流阻力与驱替压力双场匹配优化评价方法,对目前的渗流阻力场与驱替压力场的匹配性进行评价,并根据措施后双场变化情况优选最佳方案。
发明内容
本发明的目的在于提供一种渗流阻力与驱替压力双场匹配优化评价方法,对目前的渗流阻力场与驱替压力场的匹配性进行评价,并根据措施后双场变化情况优选最佳方案。
本发明解决其技术问题所采取的技术方案是:
一种渗流阻力与驱替压力双场匹配优化评价方法,包括
(a)无效循环动态边界判识,确定无效循环带;
(b)划分局部流动体系;
(c)确定潜力参数、渗流阻力参数和驱动力参数;
(d)建立驱动力与渗流阻力适配参数,对渗流阻力场与驱替压力场匹配程度评价;
(e)改变渗流压力和阻力分布场,提高双场匹配程度。
进一步的,所述步骤(a)的具体过程为:根据区块的流线分布状况可判识单相流水循环带,在出现窜流通道时,流线分布密集,密集流线经过的网格Sw≥0.95(1-Sor)作为流线切割的无效循环带,式中,Sw—含水饱和度,%;Sor—残余油饱和度,%。
进一步的,所述步骤(b)的具体过程为:将步骤(a)确定的流线切割的无效循环带作为动态边界,结合断层、夹层、侧积层、岩性边界等静态渗流屏障,将流动体系划分为若干个流动单元。
进一步的,所述步骤(c)中,
定义潜力参数:系统平均含油饱和度
式中,n—网格数,整数;—第i个网格的含油饱和度,小数;Vi—第i个网格的体积,m3;V—模型网格总体积,m3
定义渗流阻力参数:平均渗流阻力
式中,n—网格数,整数;Ki—第i个网格的渗透率,mD;—第i个网格的含水饱和度,小数;/>—含水饱和度为/>时的油相相对渗透率,小数;/>—含水饱和度为时的水相相对渗透率,小数;μo—地下原油粘度,mPa·s;μw—地层水粘度,mPa·s。
定义驱动力参数:驱动潜力
式中,n—网格数,整数;—Hamilton算子,用来表示梯度和散度;Pi—第i个网格的压力,MPa;/>—第i个网格的含油饱和度,小数;Vi—第i个网格的体积,m3,V—模型网格总体积,m3
进一步的,所述步骤(d)中,建立驱动力与渗流阻力适配参数:
其中,
式中,n—网格数,整数;Ki—第i个网格的渗透率,mD;—第i个网格的含水饱和度,小数;/>—含水饱和度为/>时的油相相对渗透率,小数;/>—含水饱和度为时的水相相对渗透率,小数;μo—地下原油粘度,mPa·s;μw—地层水粘度,mPa·s;Ri—第i个渗流阻力,MPa,/>—Hamilton算子,用来表示梯度和散度;Pi—第i个网格的压力,MPa。
进一步的,所述步骤(e)中,通过封堵、射孔、转注等各种措施手段改变流动阻力的分布状况,驱替高饱和油,控制高饱和水,根据渗流阻力场与驱替压力场的适配性参数评价优选最佳方案,适配参数越大,开发效果越好。
本发明的技术效果:
与现有技术相比,本发明的一种渗流阻力与驱替压力双场匹配优化评价方法,针对水淹水窜严重的高含水油藏的无效水循环区域,含油饱和度场与压力场不匹配的问题,根据流线分布状况判识无效循环带,并依据静态渗流屏障(断层、夹层、侧积层、岩性边界等)及无效水循环动态边界划分局部流动体系,建立潜力参数、渗流阻力参数、驱动力参数及驱动力与渗流阻力适配参数,根据各种措施手段改变流动阻力的分布状况,驱替高饱和油,控制高饱和水,根据渗流阻力场与驱替压力场的适配性参数评价优选最佳方案。
附图说明
图1为本发明的流线切割的无效循环带示意图;
图2为本发明的流动体系划分示意图;
图3为实施例中的井组含油饱和度图(调整前);
图4为实施例中的渗流阻力场与驱替压力场分布匹配程度(调整前);
图5为实施例中的井组含油饱和度图(方案一);
图6为实施例中的渗流阻力场与驱替压力场分布匹配程度(方案一);
图7为实施例中的井组含油饱和度图(方案二);
图8为实施例中的渗流阻力场与驱替压力场分布匹配程度(方案二)。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。通常在此处附图中描述和示出的本发明实施例的组件可以以各种不同的配置来布置和设计。因此,以下对在附图中提供的本发明的实施例的详细描述并非旨在限制要求保护的本发明的范围,而是仅仅表示本发明的选定实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
应注意到:相似的标号和字母在下面的附图中表示类似项,因此,一旦某一项在一个附图中被定义,则在随后的附图中不需要对其进行进一步定义和解释。
下面通过具体实施例并结合附图对本发明作进一步说明。
实施例1:
本实施例涉及的渗流阻力与驱替压力双场匹配优化评价方法,包括:
(1)根据区块的流线分布状况可判识单相流水循环带,在出现窜流通道时,流线分布密集,密集流线经过的网格(Sw≥0.95(1-Sor))作为流线切割的无效循环带,如图1所示;
(2)将确定的流线切割的无效循环带作为动态边界,结合断层、夹层、侧积层、岩性边界等静态渗流屏障,将流动体系划分为若干个流动单元,如图2所示;
(3)定义评价驱动力与渗流阻力适配性的参数。
定义潜力参数:系统平均含油饱和度
式中,n—网格数,234260;—第i个网格的含油饱和度,小数;Vi—第i个网格的体积,m3;V—模型网格总体积,187408000m3
定义渗流阻力参数:平均渗流阻力
式中,n—网格数,234260;Ki—第i个网格的渗透率,mD;—第i个网格的含水饱和度,小数;/>—含水饱和度为/>时的油相相对渗透率,小数;/>—含水饱和度为/>时的水相相对渗透率,小数;μo—地下原油粘度,1.5mPa·s;μw—地层水粘度,0.35mPa·s。
定义驱动力参数:驱动潜力
式中,n—网格数,234260;—Hamilton算子,用来表示梯度和散度;Pi—第i个网格的压力,MPa;/>—第i个网格的含油饱和度,小数;Vi—第i个网格的体积,m3,V—模型网格总体积,187408000m3
(4)定义评价驱动力与渗流阻力适配性的参数。
建立驱动力与渗流阻力适配参数:
其中,
式中,n—网格数,234260;Ki—第i个网格的渗透率,mD;—第i个网格的含水饱和度,小数;/>—含水饱和度为/>时的油相相对渗透率,小数;/>—含水饱和度为/>时的水相相对渗透率,小数;μo—地下原油粘度,1.5mPa·s;μw—地层水粘度,0.35mPa·s;Ri—第i个渗流阻力,MPa,/>—Hamilton算子,用来表示梯度和散度;Pi—第i个网格的压力,MPa。
本实施例井间局部低效水循环如图3所示,驱动力与渗流阻力适配参数为0.25,渗流阻力场与驱替压力场分布如图4所示,匹配程度差。
本实施例根据井组开发状况,制定如下两种措施方案,以改变压力和流动阻力场,并分别评价两种方案的驱动力与渗流阻力适配参数,适配参数越大,开发效果越佳。
方案一:低效水循环井间调堵
结果如图5所示,原低效水循环区域有所改善,但其他方向水循环效率有所降低,渗流阻力场与驱替压力场分布匹配程度如图6所示,驱动力与渗流阻力适配参数为0.45,仍有待提高。
方案二:低效水循环油井转注,水井关停
结果如图7所示,原低效水循环区域改善良好,其他方向水循环效率也相对均匀,渗流阻力场与驱替压力场分布匹配程度如图8所示,驱动力与渗流阻力适配参数为0.75,效果良好。
因此,方案二,即低效水循环油井转注,水井关停为最佳方案。
上述具体实施方式仅是本发明的具体个案,本发明的专利保护范围包括但不限于上述具体实施方式的产品形态和式样,任何符合本发明权利要求书且任何所属技术领域的普通技术人员对其所做的适当变化或修饰,皆应落入本发明的专利保护范围。

Claims (2)

1.一种渗流阻力与驱替压力双场匹配优化评价方法,其特征在于:包括
步骤(a)无效循环动态边界判识,确定无效循环带;
步骤(b)划分局部流动体系;
步骤(c)确定潜力参数、渗流阻力参数和驱动力参数;
步骤(d)建立驱动力与渗流阻力适配参数,对渗流阻力场与驱替压力场匹配程度评价;
步骤(e)改变驱替压力场和渗流阻力场,提高双场匹配程度;
所述步骤(a)的具体过程为:根据区块的流线分布状况判识单相流水循环带,在出现窜流通道时,流线分布密集,密集流线经过的网格Sw≥0.95(1-Sor)作为流线切割的无效循环带,式中,Sw—含水饱和度,%;Sor—残余油饱和度,%;
所述步骤(b)的具体过程为:将步骤(a)确定的流线切割的无效循环带作为动态边界,结合断层、夹层、侧积层、岩性边界静态渗流屏障,将流动体系划分为若干个流动单元;
所述步骤(c)中,定义潜力参数为系统平均含油饱和度:
式中,n—网格数,整数;—第i个网格的含油饱和度,小数;Vi—第i个网格的体积,m3;V—模型网格总体积,m3
定义渗流阻力参数为平均渗流阻力:
式中,n—网格数,整数;Ki—第i个网格的渗透率,mD;—第i个网格的含水饱和度,小数;/>—含水饱和度为/>时的油相相对渗透率,小数;/>—含水饱和度为/>时的水相相对渗透率,小数;μo—地下原油粘度,mPa·s;μw—地层水粘度,mPa·s;
定义驱动力参数为驱动潜力:
式中,n—网格数,整数;—Hamilton算子,用来表示梯度和散度;Pi—第i个网格的压力,MPa;/>—第i个网格的含油饱和度,小数;Vi—第i个网格的体积,m3,V—模型网格总体积,m3
所述步骤(d)中,建立驱动力与渗流阻力适配参数:
其中,/>
式中,n—网格数,整数;Ki—第i个网格的渗透率,mD;—第i个网格的含水饱和度,小数;/>—含水饱和度为/>时的油相相对渗透率,小数;/>—含水饱和度为/>时的水相相对渗透率,小数;μo—地下原油粘度,mPa·s;μw—地层水粘度,mPa·s;Ri—第i个渗流阻力,MPa,/>—Hamilton算子,用来表示梯度和散度;Pi—第i个网格的压力,MPa。
2.根据权利要求1所述的渗流阻力与驱替压力双场匹配优化评价方法,其特征在于:通过封堵、射孔、转注的方式,以改变驱替压力场和流动阻力场,并分别评价驱动力与渗流阻力适配参数,驱动力与渗流阻力适配参数越大,开发效果越好。
CN202010341083.6A 2020-04-26 2020-04-26 一种渗流阻力与驱替压力双场匹配优化评价方法 Active CN113554259B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010341083.6A CN113554259B (zh) 2020-04-26 2020-04-26 一种渗流阻力与驱替压力双场匹配优化评价方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010341083.6A CN113554259B (zh) 2020-04-26 2020-04-26 一种渗流阻力与驱替压力双场匹配优化评价方法

Publications (2)

Publication Number Publication Date
CN113554259A CN113554259A (zh) 2021-10-26
CN113554259B true CN113554259B (zh) 2024-04-16

Family

ID=78129930

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010341083.6A Active CN113554259B (zh) 2020-04-26 2020-04-26 一种渗流阻力与驱替压力双场匹配优化评价方法

Country Status (1)

Country Link
CN (1) CN113554259B (zh)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2061221C1 (ru) * 1994-03-01 1996-05-27 Внедренческий научно-исследовательский инженерный центр "Нефтегазтехнология" Способ определения модифицированных функций относительных фазовых проницаемостей для слоисто-неоднородных пористых сред
RU2148169C1 (ru) * 1998-09-17 2000-04-27 Внедренческий научно-исследовательский инженерный центр "Нефтегазтехнология" Способ контроля за разработкой нефтяного месторождения со слоисто-неоднородными пластами
CN105089612A (zh) * 2014-05-04 2015-11-25 中国石油化工股份有限公司 低渗透油藏人工裂缝压裂缝长与井排距确定方法
CN105545267A (zh) * 2015-12-09 2016-05-04 东北石油大学 一种实现变渗流阻力驱油的方法
CN106126829A (zh) * 2016-06-27 2016-11-16 西南石油大学 一种评价驱油缔合聚合物油藏适应性的方法
CN107346518A (zh) * 2016-05-06 2017-11-14 中国石油化工股份有限公司 致密低渗透油藏油水两相流最大渗流阻力梯度的获取方法
CN109753701A (zh) * 2018-12-21 2019-05-14 昆明理工大学 一种随机渗流场作用下的土质边坡可靠度分析下限法
CN109858177A (zh) * 2019-02-23 2019-06-07 中国海洋石油集团有限公司 一种基于拟流线法的底水油藏水平井水驱数值模拟方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2969206B1 (fr) * 2010-12-15 2013-11-01 IFP Energies Nouvelles Procede d'exploitation d'un gisement petrolier a partir d'un modele d'ecoulement cale au moyen d'une mise a l'echelle de cartes de pression et de saturation

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2061221C1 (ru) * 1994-03-01 1996-05-27 Внедренческий научно-исследовательский инженерный центр "Нефтегазтехнология" Способ определения модифицированных функций относительных фазовых проницаемостей для слоисто-неоднородных пористых сред
RU2148169C1 (ru) * 1998-09-17 2000-04-27 Внедренческий научно-исследовательский инженерный центр "Нефтегазтехнология" Способ контроля за разработкой нефтяного месторождения со слоисто-неоднородными пластами
CN105089612A (zh) * 2014-05-04 2015-11-25 中国石油化工股份有限公司 低渗透油藏人工裂缝压裂缝长与井排距确定方法
CN105545267A (zh) * 2015-12-09 2016-05-04 东北石油大学 一种实现变渗流阻力驱油的方法
CN107346518A (zh) * 2016-05-06 2017-11-14 中国石油化工股份有限公司 致密低渗透油藏油水两相流最大渗流阻力梯度的获取方法
CN106126829A (zh) * 2016-06-27 2016-11-16 西南石油大学 一种评价驱油缔合聚合物油藏适应性的方法
CN109753701A (zh) * 2018-12-21 2019-05-14 昆明理工大学 一种随机渗流场作用下的土质边坡可靠度分析下限法
CN109858177A (zh) * 2019-02-23 2019-06-07 中国海洋石油集团有限公司 一种基于拟流线法的底水油藏水平井水驱数值模拟方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
低渗透储层水驱油渗流阻力特征;苏海波;于金彪;张同伍;张明安;刘洋;;油气地质与采收率;20171207(01);全文 *
低渗透油藏毛细作用评价方法研究;耿向飞;罗健辉;丁彬;王平美;彭宝亮;贺丽鹏;;油田化学;20180118(04);全文 *
水驱砂岩油藏优势渗流通道识别;刘薇薇;龚丽荣;罗福全;温玉焕;王力那;;复杂油气藏;20200325(01);全文 *
特低渗透油藏水驱规律及最佳驱替模式;王文环;《石油天然气工业》;20190228;第40卷(第1期);全文 *

Also Published As

Publication number Publication date
CN113554259A (zh) 2021-10-26

Similar Documents

Publication Publication Date Title
CN105626010B (zh) 一种分段注水井中注水层段合理划分方法
CN110439515B (zh) 注采参数优化方法及装置
CN103939066B (zh) 一种一注多采井组定注水量确定油井产液量的方法
CN103422849B (zh) 注水开发老油田井网重建方法及装置
CN109209312B (zh) 一种适用聚合物驱的资源开发方法
CN111648751B (zh) 一种超高轮次吞吐后期多层系稠油油藏开发方法
CN113554259B (zh) 一种渗流阻力与驱替压力双场匹配优化评价方法
US11702933B2 (en) Method for determining vertical well stratified fracturing crack lengths based on interlayer equilibrium displacement
DE102008007627B3 (de) Erdwärmeanlage
CN116108572A (zh) 页岩凝析气井体积压裂外区产能贡献分析方法
CN107035348A (zh) 一种油田调剖用多因素分析选井方法及装置
Denbina et al. Modelling cold production for heavy oil reservoirs
CN111749688B (zh) 一种优势渗流通道发育层位和方向的预测方法
CN109989747B (zh) 一种油田注采井组开发非均匀程度定量表征方法
RU2290501C1 (ru) Способ разработки нефтяной залежи
CN113090232B (zh) 提高复杂断块高饱和砂岩油藏采收率的方法
CN112502677A (zh) 一种基于多元线性回归的注水开发效果评价方法
RU2563463C1 (ru) Способ разработки слоистой нефтяной залежи с высоковязкой нефтью
CN110939425A (zh) 一种正对水平井井网的视反七点法交替开采方法及装置
CN115081694B (zh) 海上薄互层砂岩油田加密潜力量化评价方法及装置
CN113187472B (zh) 一种层状砂岩油藏水驱开发渗流优势通道的识别方法
CN114198072B (zh) 多层合采油藏水驱开发井网优化部署方法
CN116291346B (zh) 一种用于优化纵向非均质稠油油藏泡沫调驱体系的图版确定方法
CN114592838B (zh) 一种层状超稠油油藏蒸汽驱潜力评价方法及系统
CN114718529B (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