CN109594979A - 基于井轨迹的分流砂坝内部夹层参数计算方法 - Google Patents

基于井轨迹的分流砂坝内部夹层参数计算方法 Download PDF

Info

Publication number
CN109594979A
CN109594979A CN201811360634.2A CN201811360634A CN109594979A CN 109594979 A CN109594979 A CN 109594979A CN 201811360634 A CN201811360634 A CN 201811360634A CN 109594979 A CN109594979 A CN 109594979A
Authority
CN
China
Prior art keywords
well
sand dam
interlayer
shunting
dam
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
Application number
CN201811360634.2A
Other languages
English (en)
Other versions
CN109594979B (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 National Offshore Oil Corp CNOOC
CNOOC China Ltd Tianjin Branch
Original Assignee
China National Offshore Oil Corp CNOOC
CNOOC China Ltd Tianjin Branch
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 National Offshore Oil Corp CNOOC, CNOOC China Ltd Tianjin Branch filed Critical China National Offshore Oil Corp CNOOC
Priority to CN201811360634.2A priority Critical patent/CN109594979B/zh
Publication of CN109594979A publication Critical patent/CN109594979A/zh
Application granted granted Critical
Publication of CN109594979B publication Critical patent/CN109594979B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21BEARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B49/00Testing 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
    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21BEARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B47/00Survey of boreholes or wells

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Geology (AREA)
  • Mining & Mineral Resources (AREA)
  • Physics & Mathematics (AREA)
  • Environmental & Geological Engineering (AREA)
  • Fluid Mechanics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Geophysics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

一种基于井轨迹的分流砂坝内部夹层参数计算方法:是在分流砂坝型浅水三角洲沉积模式指导下,通过地震正演模拟提取对单一分流砂坝边界敏感的各项属性进行融合,根据融合属性定性刻画单一分流砂坝平面分布。同时以地震低角度前积为出发点,建立前积理论计算模型定量计算单一分流砂坝的边界。在单一分流砂坝内部,根据定向井和水平井井轨迹和夹层的接触关系建立分流砂坝内部夹层参数理论计算模型,并从钻遇夹层的井资料中提取相关地质参数,最终定量计算分流砂坝内部夹层的倾角、间距和厚度。本发明充分利用井轨迹信息,建立井轨迹信息与地质体定量关系,重现了单一分流砂坝内部夹层的展布,使刻画结果更加准确可靠。

Description

基于井轨迹的分流砂坝内部夹层参数计算方法
技术领域
本发明涉及一种分流砂坝内部夹层参数计算方法。特别是涉及一种基于井轨迹的分流砂坝内部夹层参数计算方法。
背景技术
区别于常规正常三角洲,浅水三角洲作为一种特殊的沉积类型,发育于水体较浅、构造相对平缓、整体沉降缓慢的坳陷湖盆中,根据沉积骨架的差异可以分为分流河道型和分流砂坝型。近年来,随着现代沉积观察及水槽沉积模拟被广泛应用于浅水三角洲储层研究,对分流砂坝型浅水三角洲的认识也进一步得到了加深。
目前针对分流砂坝型浅水三角洲单一分流砂坝内部夹层,最主要的研究手段是根据沉积模式利用定向井井资料进行井间预测。基本原理是在精细等时地层格架(小层或者单层)内部对单一期次的成因砂体进行识别,进而在单一成因砂体内部以模式为指导,根据实钻定向井资料及生态动态资料研究夹层的分布,该方法适用于密井网油田,研究难度大,多解性强,而且对研究人员的经验要求高。对于井网相对稀疏并且水平井比例较高的海上油田,该方法基本不具备可实施性,实现精确定量化研究难度大,不利于海上油田的高效开发。
现有技术中缺乏在井距相对稀疏的水平井-定向井交错井网条件下精细刻画分流砂坝型浅水三角洲单一分流砂坝内部夹层的方法的问题。因此,寻求一种适用于水平井-定向井交错井网条件下的分流砂坝内部夹层计算方法对于海上油田高效开发意义重大。
发明内容
本发明所要解决的技术问题是,提供一种适用于水平井-定向井交错井网条件下的分流砂坝内部夹层计算的基于井轨迹的分流砂坝内部夹层参数计算方法。
本发明所采用的技术方案是:一种基于井轨迹的分流砂坝内部夹层参数计算方法,包括如下步骤:
1)根据高分辨率层序地层学理论,针对目的层进行层序划分,选取稳定湖泛面为标志层,确定地层发育模式,通过井震结合、分级控制,建立研究区等时地层格架;
2)以取心井岩心观察为基础,识别取心井岩心内部单一分流砂坝,并获单一分流砂坝的规模、电性特征及接触关系;
3)根据岩心观察所获取的单一分流砂坝的规模及接触关系,构建地震正演模拟的原型模型;
4)根据地震正演模拟的原型模型,从地质研究一体化软件地震模块中选择表征单一分流砂坝边界的地震属性,并根据分流砂坝型浅水三角洲沉积模式定性刻画分流砂坝;
5)在定性刻画分流砂坝的基础上,以地震低角度前积为出发点,根据沉积过程理论,建立前积理论计算模型;
6)根据实际钻井资料提取前积理论计算模型相关计算参数,定量计算单一期次前积分流砂坝的边界,并依据该定量计算结果对步骤4)所刻画的分流砂坝进行修正,最终精确刻画出各单一分流砂坝平面分布;
7)在精确刻画出各单一分流砂坝平面分布的基础上,考虑到海上油田井距相对较大和水平井比例相对多的情况,综合利用水平井和定向井井轨迹及测井曲线信息建立适合海上油田的分流砂坝内部侧叠式夹层的定量计算模型。
步骤1)包括:综合对比岩心-测井-地震资料,结合高分辨率层序地层学理论,识别目的层之上稳定的泥岩段为最大湖泛面,识别三级层序界面,并在三级层序内部识别小于三级的层序单元,并最终建立目的层的单砂体等时地层格架。
步骤3)是根据单一分流砂坝的规模及接触关系,建立以下三种地震正演模拟模型,包括:非均质性单一分流砂坝侧向叠置,高程有差异的分流砂坝侧向叠置和分流砂坝孤立分布;
步骤4)所述的地质研究一体化软件地震模块,采用discovery或petrel地质研究一体化模块。
步骤5)所述的计算模型为:
Δs1=h2*Δs/Δh (1)
式中:Δs1为单一成因砂体井点外推尖灭距离;h2为井点钻遇储层的垂厚;Δs为钻遇单一成因砂体的两口井的间距;Δh为钻遇单一成因砂体的两口井的垂厚差。
步骤7)包括:
根据水平井-定向井交错井网,采用水平井确定夹层厚度、间距,定向井确定夹层倾角的定量刻画方法,通过小于150m井距的定向对子井实钻资料确定单一分流砂坝内部侧叠式夹层的倾角;然后根据水平井井轨迹不同井斜段实钻数据,抽象出斜列式夹层的厚度、间距定量计算公式如下:
β=arctan(H1/ΔL) (2)
式中:β为夹层倾角;H1为不同定向井点钻遇夹层的垂距;ΔL为不同定向井点钻遇夹层的水平间距;
S=(ΔD)sinα (3)
式中:S为夹层间距;ΔD为水平井钻遇相邻两个夹层间的轨迹长度;α为水平井钻遇点时的井斜角;
h=L3sin(arctan(H1/(L1sinα-L3) (4)
式中:h为夹层厚度;L1为井点钻遇夹层段的轨迹长度;H1为不同定向井点钻遇夹层的垂距;α为水平井钻遇点时的井斜角;L3为L1的水平投影和H1的水平投影之间的长度差。
本发明的基于井轨迹的分流砂坝内部夹层参数计算方法,充分挖掘了地震资料信息,从成因角度出发大幅提高了分流砂坝的定量刻画精度,将近似连片分布的砂体转换为平面上不同的单一分流砂坝,并充分利用海上油田水平井比例高,水平井资料侧向分辨率高的特点,建立适合海上油田的分流砂坝内部侧叠式夹层的定量计算模型。本发明提供一种针对分流砂坝内部夹层参数的获取方法,充分利用井轨迹信息,建立井轨迹信息与地质体定量关系,重现了单一分流砂坝内部夹层的展布,使刻画结果更加准确可靠。
附图说明
图1是本发明基于井轨迹的分流砂坝内部夹层参数计算方法的流程图;
图2是本发明中不同接触关系类型的分流砂坝地震正演模拟示意图;
图3是本发明中单一前积分流砂坝边界理论计算示意图;
图4是本发明中利用井轨迹定量求取单一分流砂坝内部夹层参数示意图。
具体实施方式
下面结合实施例和附图对本发明的基于井轨迹的分流砂坝内部夹层参数计算方法做出详细说明。
如图1所示,本发明的基于井轨迹的分流砂坝内部夹层参数计算方法,包括如下步骤:
1)根据高分辨率层序地层学理论,针对目的层进行层序划分,选取稳定湖泛面为标志层,确定地层发育模式,通过井震结合、分级控制,建立研究区等时地层格架;包括:综合对比岩心-测井-地震资料,结合高分辨率层序地层学理论,识别目的层之上稳定的泥岩段为最大湖泛面,识别三级层序界面,并在三级层序内部识别小于三级的层序单元,并最终建立目的层的单砂体等时地层格架。
2)以取心井岩心观察为基础,识别取心井岩心内部单一分流砂坝,并获单一分流砂坝的规模、电性特征及接触关系;
3)根据岩心观察所获取的单一分流砂坝的规模及接触关系,构建地震正演模拟的原型模型;是根据单一分流砂坝的规模及接触关系,建立以下三种地震正演模拟模型,包括:非均质性单一分流砂坝侧向叠置,高程有差异的分流砂坝侧向叠置和分流砂坝孤立分布,如图2所示;根据实际情况,可选用45Hz雷克子波激发,其中储层厚度为8m,砂岩速度为2600m/s,密度为2100kg/m3,泥岩速度为3100m/s,密度为2350kg/m3
4)根据地震正演模拟的原型模型,从地质研究一体化软件地震模块中选择表征单一分流砂坝边界的地震属性,所述的地质研究一体化软件地震模块,可以采用discovery或petrel等地质研究一体化模块。并根据分流砂坝型浅水三角洲沉积模式定性刻画分流砂坝;需要说明的是,从地震正演模拟的原型模型来看,相对于传统振幅类属性,振幅变化率对单一分流砂坝的边界具有更好的识别能力。同时考虑到最小振幅属性对砂体宏观分布具有良好的响应,故将两者属性融合。
5)考虑到应用地震沉积学方法刻画沉积体平面分布具有一定的主观性和多解性,在定性刻画分流砂坝的基础上,以地震低角度前积为出发点,根据沉积过程理论,建立前积理论计算模型,如图3所示;所述的前积理论计算模型为:
Δs1=h2*Δs/Δh (1)
式中:Δs1为单一成因砂体井点外推尖灭距离;h2为井点钻遇储层的垂厚;Δs为钻遇单一成因砂体的两口井的间距;Δh为钻遇单一成因砂体的两口井的垂厚差。
6)根据实际钻井资料提取前积理论计算模型相关计算参数,定量计算单一期次前积分流砂坝的边界,并依据该定量计算结果对步骤4)所刻画的分流砂坝进行修正,最终精确刻画出各单一分流砂坝平面分布;
研究区隔夹层存在泥质、钙质和物性三种类型,按照级次可划分为砂层组间隔夹层、小层间隔夹层和层内夹层,研究所述的隔夹层为层内夹层,其主要为物性夹层,部分区域为泥质夹层,层内夹层存在于单一分流砂坝内部不同增生体叠置的界面处,其主要向单一方向侧向斜列,即产状为侧叠式。
水平井在进行测井时,其井筒周围的环境与定向井有明显差异,但两者在夹层发育的位置具有相似的测井响应特征。分流砂坝内部夹层主要呈侧叠式,水平井水平分辨率高,在较小的间距范围内,井轨迹钻穿同一夹层。利用这一特性,结合钻遇同一分流砂坝内的定向井,可以有效提高分流砂坝内部夹层的刻画精度。
7)在精确刻画出各单一分流砂坝平面分布的基础上,考虑到海上油田井距相对较大(平均约250~350米)和水平井比例相对多的情况,综合利用水平井和定向井井轨迹及测井曲线信息建立适合海上油田的分流砂坝内部侧叠式夹层的定量计算模型。包括:
根据水平井-定向井交错井网,采用水平井确定夹层厚度、间距,定向井确定夹层倾角的定量刻画方法,通过小于150m井距的定向对子井实钻资料确定单一分流砂坝内部侧叠式夹层的倾角;然后根据水平井井轨迹不同井斜段实钻数据,抽象出斜列式夹层的厚度、间距定量计算公式如下(如图4所示):
β=arctan(H1/ΔL) (2)
式中:β为夹层倾角;H1为不同定向井点钻遇夹层的垂距;ΔL为不同定向井点钻遇夹层的水平间距;
S=(ΔD)sinα (3)
式中:S为夹层间距;ΔD为水平井钻遇相邻两个夹层间的轨迹长度;α为水平井钻遇点时的井斜角;
h=L3sin(arctan(H1/(L1sinα-L3) (4)
式中:h为夹层厚度;L1为井点钻遇夹层段的轨迹长度;H1为不同定向井点钻遇夹层的垂距;α为水平井钻遇点时的井斜角;L3为L1的水平投影和H1的水平投影之间的长度差。

Claims (6)

1.一种基于井轨迹的分流砂坝内部夹层参数计算方法,其特征在于,包括如下步骤:
1)根据高分辨率层序地层学理论,针对目的层进行层序划分,选取稳定湖泛面为标志层,确定地层发育模式,通过井震结合、分级控制,建立研究区等时地层格架;
2)以取心井岩心观察为基础,识别取心井岩心内部单一分流砂坝,并获单一分流砂坝的规模、电性特征及接触关系;
3)根据岩心观察所获取的单一分流砂坝的规模及接触关系,构建地震正演模拟的原型模型;
4)根据地震正演模拟的原型模型,从地质研究一体化软件地震模块中选择表征单一分流砂坝边界的地震属性,并根据分流砂坝型浅水三角洲沉积模式定性刻画分流砂坝;
5)在定性刻画分流砂坝的基础上,以地震低角度前积为出发点,根据沉积过程理论,建立前积理论计算模型;
6)根据实际钻井资料提取前积理论计算模型相关计算参数,定量计算单一期次前积分流砂坝的边界,并依据该定量计算结果对步骤4)所刻画的分流砂坝进行修正,最终精确刻画出各单一分流砂坝平面分布;
7)在精确刻画出各单一分流砂坝平面分布的基础上,考虑到海上油田井距相对较大和水平井比例相对多的情况,综合利用水平井和定向井井轨迹及测井曲线信息建立适合海上油田的分流砂坝内部侧叠式夹层的定量计算模型。
2.根据权利要求1所述的基于井轨迹的分流砂坝内部夹层参数计算方法,其特征在于,步骤1)包括:综合对比岩心-测井-地震资料,结合高分辨率层序地层学理论,识别目的层之上稳定的泥岩段为最大湖泛面,识别三级层序界面,并在三级层序内部识别小于三级的层序单元,并最终建立目的层的单砂体等时地层格架。
3.根据权利要求1所述的基于井轨迹的分流砂坝内部夹层参数计算方法,其特征在于,步骤3)是根据单一分流砂坝的规模及接触关系,建立以下三种地震正演模拟模型,包括:非均质性单一分流砂坝侧向叠置,高程有差异的分流砂坝侧向叠置和分流砂坝孤立分布。
4.根据权利要求1所述的基于井轨迹的分流砂坝内部夹层参数计算方法,其特征在于,步骤4)所述的地质研究一体化软件地震模块,采用discovery或petrel地质研究一体化模块。
5.根据权利要求1所述的基于井轨迹的分流砂坝内部夹层参数计算方法,其特征在于,步骤5)所述的计算模型为:
Δs1=h2*Δs/Δh(1)
式中:Δs1为单一成因砂体井点外推尖灭距离;h2为井点钻遇储层的垂厚;Δs为钻遇单一成因砂体的两口井的间距;Δh为钻遇单一成因砂体的两口井的垂厚差。
6.根据权利要求1所述的基于井轨迹的分流砂坝内部夹层参数计算方法,其特征在于,步骤7)包括:
根据水平井-定向井交错井网,采用水平井确定夹层厚度、间距,定向井确定夹层倾角的定量刻画方法,通过小于150m井距的定向对子井实钻资料确定单一分流砂坝内部侧叠式夹层的倾角;然后根据水平井井轨迹不同井斜段实钻数据,抽象出斜列式夹层的厚度、间距定量计算公式如下:
β=arctan(H1/ΔL)(2)
式中:β为夹层倾角;H1为不同定向井点钻遇夹层的垂距;ΔL为不同定向井点钻遇夹层的水平间距;
S=(ΔD)sinα (3)
式中:S为夹层间距;ΔD为水平井钻遇相邻两个夹层间的轨迹长度;α为水平井钻遇点时的井斜角;
h=L3sin(arctan(H1/(L1sinα-L3)(4)
式中:h为夹层厚度;L1为井点钻遇夹层段的轨迹长度;H1为不同定向井点钻遇夹层的垂距;α为水平井钻遇点时的井斜角;L3为L1的水平投影和H1的水平投影之间的长度差。
CN201811360634.2A 2018-11-15 2018-11-15 基于井轨迹的分流砂坝内部夹层参数计算方法 Active CN109594979B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811360634.2A CN109594979B (zh) 2018-11-15 2018-11-15 基于井轨迹的分流砂坝内部夹层参数计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811360634.2A CN109594979B (zh) 2018-11-15 2018-11-15 基于井轨迹的分流砂坝内部夹层参数计算方法

Publications (2)

Publication Number Publication Date
CN109594979A true CN109594979A (zh) 2019-04-09
CN109594979B CN109594979B (zh) 2022-04-22

Family

ID=65958147

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811360634.2A Active CN109594979B (zh) 2018-11-15 2018-11-15 基于井轨迹的分流砂坝内部夹层参数计算方法

Country Status (1)

Country Link
CN (1) CN109594979B (zh)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090204377A1 (en) * 2004-09-10 2009-08-13 Van Wagoner John C Method for Constructing Geologic Models of Subsurface Sedimentary Volumes
CN103218750A (zh) * 2013-04-03 2013-07-24 中国海洋石油总公司 一种海上沉积岩油田地质评估方法
CN103760617A (zh) * 2013-12-26 2014-04-30 中国石油天然气股份有限公司 浅水湖盆三角洲生长规律与沉积模式的获取方法及装置
CN104018827A (zh) * 2014-05-13 2014-09-03 长江大学 一种基于叠覆朵体的三角洲砂体内部结构分析方法
CN107490810A (zh) * 2017-08-17 2017-12-19 中国地质大学(北京) 基于风场、物源、盆地系统的储集砂体预测方法和装置
CN108763732A (zh) * 2018-05-24 2018-11-06 长江大学 一种分流砂坝型三角洲前缘训练图像建立方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090204377A1 (en) * 2004-09-10 2009-08-13 Van Wagoner John C Method for Constructing Geologic Models of Subsurface Sedimentary Volumes
CN103218750A (zh) * 2013-04-03 2013-07-24 中国海洋石油总公司 一种海上沉积岩油田地质评估方法
CN103760617A (zh) * 2013-12-26 2014-04-30 中国石油天然气股份有限公司 浅水湖盆三角洲生长规律与沉积模式的获取方法及装置
CN104018827A (zh) * 2014-05-13 2014-09-03 长江大学 一种基于叠覆朵体的三角洲砂体内部结构分析方法
CN107490810A (zh) * 2017-08-17 2017-12-19 中国地质大学(北京) 基于风场、物源、盆地系统的储集砂体预测方法和装置
CN108763732A (zh) * 2018-05-24 2018-11-06 长江大学 一种分流砂坝型三角洲前缘训练图像建立方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
赵汉卿等: "渤海湾盆地A油田新近系明下段沉积特征及演化模式", 《吉林大学学报(地球科学版)》 *

Also Published As

Publication number Publication date
CN109594979B (zh) 2022-04-22

Similar Documents

Publication Publication Date Title
CN104632200B (zh) 河流相多层及厚层油藏精细建模方法
CN109441422A (zh) 一种页岩气井间距优化开采方法
CN108301825B (zh) 一种高温高压储层孔隙压力分析方法
CN107783187B (zh) 一种将测井速度和地震速度结合建立三维速度场的方法
CN105182444A (zh) 一种高分辨率层序地层格架约束的地质统计学反演方法
CN106401574B (zh) 一种钻前高温地热井地层压力的预测方法
CN106837319B (zh) 一种密井网条件下的地下单一古河道识别方法
CN109143373B (zh) 页岩气储层地层压力计算方法及计算机可读存储介质
US20100305927A1 (en) Updating a reservoir model using oriented core measurements
Xiaoqing et al. An in-depth analysis of reservoir architecture of underwater distributary channel sand bodies in a river dominated delta: A case study of T51 Block, Fuyu Oilfield, Jilin
CN106501873B (zh) 一种基于岩心恢复三角洲前缘前积层沉积倾角的方法
CN109541685A (zh) 一种河道砂体识别方法
CN108678672A (zh) 深水浅部复杂岩性地层导管高效置入及井口稳定预测方法
CN106285642A (zh) 一种基于地震资料的地层孔隙压力预测方法
Zhifeng et al. Application of seismic architecture interpretation in enhancing oil recovery in late development Stage——Taking meandering river reservoir in Dongying depression as an example
CN110532507A (zh) 一种提高水平井致密油储层钻遇率的方法
Liu et al. The Control Theory and Application for Well Pattern Optimization of Heterogeneous Sandstone Reservoirs
CN109594979A (zh) 基于井轨迹的分流砂坝内部夹层参数计算方法
CN105239998A (zh) 一种基于海上油田储层细分的部井方法
Wei et al. Application of geological model dynamic updating method in Guanjiapu Oilfield, Dagang
Luo et al. A quantitative study of the scale and distribution of tight gas reservoirs in the Sulige gas field, Ordos Basin, northwest China
CN110244357A (zh) 一种构造破坏型古油藏的判识方法
Zhao et al. Research on meandering river reservoir deposition architecture and 3D modeling of the Gudao Oil field in the Bohai Bay Basin
CN113514879B (zh) “耳朵”层识别方法
CN107065034B (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