CN111061892B - 一种山前地区沉积相带分布范围勘测方法 - Google Patents

一种山前地区沉积相带分布范围勘测方法 Download PDF

Info

Publication number
CN111061892B
CN111061892B CN201911286459.1A CN201911286459A CN111061892B CN 111061892 B CN111061892 B CN 111061892B CN 201911286459 A CN201911286459 A CN 201911286459A CN 111061892 B CN111061892 B CN 111061892B
Authority
CN
China
Prior art keywords
curve
standard
logging curve
preset
logging
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
CN201911286459.1A
Other languages
English (en)
Other versions
CN111061892A (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.)
Chengdu Univeristy of Technology
Original Assignee
Chengdu Univeristy of Technology
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 Chengdu Univeristy of Technology filed Critical Chengdu Univeristy of Technology
Priority to CN201911286459.1A priority Critical patent/CN111061892B/zh
Publication of CN111061892A publication Critical patent/CN111061892A/zh
Application granted granted Critical
Publication of CN111061892B publication Critical patent/CN111061892B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/50Information retrieval; Database structures therefor; File system structures therefor of still image data
    • G06F16/53Querying
    • G06F16/532Query formulation, e.g. graphical querying
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/002D [Two Dimensional] image generation
    • G06T11/20Drawing from basic elements, e.g. lines or circles
    • G06T11/203Drawing of straight lines or curves
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/30Assessment of water resources

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Data Mining & Analysis (AREA)
  • Databases & Information Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)

Abstract

本发明提供一种山前地区沉积相带分布范围勘测方法,通过将所采集得到的测井曲线中的曲线形态与预设的标准测井曲线的标准曲线形态进行图形比对,根据预设的标准测井曲线的标准曲线形态确定采集到的测井曲线中的曲线形态的具体类型,避免了由于人为因素而导致测井曲线中曲线形态的误判,本发明基于简单的图像比对的方法,实现了准确确定测井曲线的曲线形态的类型,从而确保后续准确确定环境中沉积相带分布情况。

Description

一种山前地区沉积相带分布范围勘测方法
技术领域
本发明属于石油天然气勘探与开发技术领域,尤其涉及一种山前地区沉积相带分布范围勘测方法。
背景技术
沉积相带是指在一定沉积环境中形成的所有沉积岩的综合,有利的沉积相带是控制油气大规模聚集的重要条件,沉积相带的分布规律是制约油气勘探进程的重要因素,沉积相带的延伸距离能够很好的反映沉积相带的平面展布情况。
目前,挤压盆地的山前地区是我国西部重要的油气勘探领域,厘定山前地区沉积相带的延伸距离能够有效的预测沉积相带的展布范围,进而明确油气勘探的有利地区。现有技术中,勘测沉积相带分布范围的方法包括基于地震数据的分析方法、基于砂砾岩百分含量平面分布法、基于测井曲线分析的方法等。其中,地震数据的提取需要高品质的三维地震资料,砂砾岩百分含量的提取需要高密度的钻井资料,而受山前地区复杂的构造特征及复杂的自然地理条件的限制,山前地区的三维地震资料提取非常困难,同时钻井提取砂砾岩百分含量的成本又较高。相对来说,使用测井曲线勘测沉积相带分布范围更加便捷和高效。
测井曲线是指在测井时形成的曲线反应出不同岩性、层位特征,进而根据所得曲线判断出具体岩性、层位等,根据测井曲线的形态、幅度、光滑程度、组合特征等方面进行测井相的分析,既可识别不同的沉积相分布情况。但是,在对测井曲线的形态进行分析时,当由于测量数据或地质环境的不同,实际测量得出的各自形态下的测井曲线的结构可能非常接近(见附图1,例如对称齿形和反形齿形),对于工作经验较少的勘测人员来说极容易出现错别的判断,即容易基于测井曲线的形态结构得出错误的沉积相带分布结论。因此,如何准确识别出测井曲线的形态,从而得出准确的沉积相带分布结论是当前有待解决问题。
发明内容
针对上述现有技术存在的缺陷,本发明提供一种山前地区沉积相带分布范围勘测方法,该方法基于准确、高效识别出测井曲线中曲线形态的基础上,以测井曲线分析可以准确得出山前地区沉积相带分布范围。
为达到上述目的,本发明采用如下技术方案:
一种山前地区沉积相带分布范围勘测方法,具体包括以下步骤:
S1、采集获取当前待测钻井的测井曲线;
S2、提取所述当前待测钻井的测井曲线的曲线幅度及幅度差、曲线形态、曲线旋回形态、旋回幅度、曲线光滑程度、包络线;
S3、将所提取的测井曲线中的曲线形态与预设的标准测井曲线的标准曲线形态进行图形比对,当提取的测井曲线中的曲线形态与预设的标准测井曲线的标准曲线形态之间的重合度在预设阈值范围时,则根据预设的标准测井曲线的标准曲线形态确定所提取的测井曲线中的曲线形态的具体类型;
S4、基于所确定的测井曲线中的曲线形态的具体类型以及所提取的测井曲线的曲线幅度及幅度差、曲线旋回形态、旋回幅度、曲线光滑程度、包络线确定当前环境中沉积相带分布情况。
进一步的,所述待测钻井的测井曲线的曲线形态类型包括钟形、漏斗形、箱形、对称齿形、反向齿形、正向齿形;所述预设的标准测井曲线的标准曲线形态类型包括标准钟形、标准漏斗形、标准箱形、标准对称齿形、标准反向齿形、标准正向齿形。
进一步的,由于实际测量中可能存在环境因素等干扰,使得所提取的测井曲线中的钟形、漏斗形、箱形、对称齿形、反向齿形、正向齿形形态曲线中没有与预设的标准测井曲线中的标准钟形、标准漏斗形、标准箱形、标准对称齿形、标准反向齿形、标准正向齿形标准形态曲线相匹配的形态图像,为了提出由于干扰导致所提取的测井曲线中的形态不客观、不准确的情况,所述步骤S3中进一步包括,当提取的测井曲线中的曲线形态与预设的标准测井曲线的标准曲线形态之间的重合度不在预设阈值范围时,则返回步骤S1。
进一步的,所述提取的测井曲线中的曲线形态与预设的标准测井曲线的标准曲线形态的图像尺寸和像素相同。
进一步的,所述提取的测井曲线中的曲线形态与预设的标准测井曲线的标准曲线形态之间的重合度,具体指具有所提取的测井曲线中的曲线形态的图像与具有预设的标准测井曲线的标准曲线形态的图像两者之间的重合比例。
进一步的,所述预设阈值范围为95%-100%。
进一步的,所述测井曲线为自然电位测井曲线、自然伽马曲线、微球电阻率曲线或侧向电阻率曲线。
本发明与现有技术相比,有以下突出的实质性的特点和显著的优势:
本发明通过将所采集得到的测井曲线中的曲线形态与预设的标准测井曲线的标准曲线形态进行图形比对,根据预设的标准测井曲线的标准曲线形态确定采集到的测井曲线中的曲线形态的具体类型,避免了由于人为因素而导致测井曲线中曲线形态的误判,本发明通过简单的图像比对的方法,实现了准确确定测井曲线的曲线形态的类型,从而确保后续准确确定环境中沉积相带分布情况。
附图说明
为了更清楚地说明本发明具体实施方式或相关技术中的技术方案,下面将对具体实施方式或相关技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施方式,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是现有技术中测井曲线识别沉积相带分布情况的曲线形态及特征表。
图2是本发明涉及的一种山前地区沉积相带分布范围勘测方法流程图。
具体实施方式
下面将结合附图对本发明的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
实施例一:
一种山前地区沉积相带分布范围勘测方法,具体包括以下步骤:
S1、采集获取当前待测钻井的测井曲线;
S2、提取所述当前待测钻井的测井曲线的曲线幅度及幅度差、曲线形态、曲线旋回形态、旋回幅度、曲线光滑程度、包络线;
S3、将所提取的测井曲线中的曲线形态与预设的标准测井曲线的标准曲线形态进行图形比对,当提取的测井曲线中的曲线形态与预设的标准测井曲线的标准曲线形态之间的重合度在预设阈值范围时,则根据预设的标准测井曲线的标准曲线形态确定所提取的测井曲线中的曲线形态的具体类型,否则返回步骤S1;
S4、基于所确定的测井曲线中的曲线形态的具体类型以及所提取的测井曲线的曲线幅度及幅度差、曲线旋回形态、旋回幅度、曲线光滑程度、包络线确定当前环境中沉积相带分布情况。
在本实施例中,所述待测钻井的测井曲线的形态包括钟形、漏斗形、箱形、对称齿形、反向齿形、正向齿形;所述预设的标准测井曲线的标准形态包括标准钟形、标准漏斗形、标准箱形、标准对称齿形、标准反向齿形、标准正向齿形。
在本实施例中,所述步骤S3具体包括步骤:
S31、工作人员根据所提取的测井曲线中的形态初步判断出该形态所代表的类型;
S32、然后将初步得出的测井曲线的曲线形态类型与预设的标准测井曲线中对应同类型的标准曲线形态进行比对;
S33、当初步得出的测井曲线的形态类型与预设的标准测井曲线中对应同类型的标准形态之间的重合度在预设阈值范围时,则根据预设的标准测井曲线的标准形态确定所提取的测井曲线中的形态的具体类型;
S34、当初步得出的测井曲线的形态类型与预设的标准测井曲线中对应同类型的标准形态之间的重合度不在预设阈值范围,则将初步得出的测井曲线的形态类型与预设的标准测井曲线中其它类型的标准形态进行比对,当初步得出的测井曲线的形态类型与预设的标准测井曲线中其它类型的标准形态之间的重合度在预设阈值范围时,则根据预设的标准测井曲线的标准形态确定所提取的测井曲线中的形态的具体类型;
S35、若初步得出的测井曲线的形态类型与预设的标准测井曲线中对所有类型的标准形态之间的重合度都不在预设阈值范围,则返回步骤S1。
在本实施例中,所述预设阈值范围为98%-100%。
实施例二:
一种山前地区沉积相带分布范围勘测方法,具体包括以下步骤:
S1、采集获取当前待测钻井的测井曲线;
S2、提取所述当前待测钻井的测井曲线的曲线幅度及幅度差、曲线形态、曲线旋回形态、旋回幅度、曲线光滑程度、包络线;
S3、将所提取的测井曲线中的曲线形态与预设的标准测井曲线的标准曲线形态进行图形比对,当提取的测井曲线中的曲线形态与预设的标准测井曲线的标准曲线形态之间的重合度在预设阈值范围时,则根据预设的标准测井曲线的标准曲线形态确定所提取的测井曲线中的曲线形态的具体类型,否则返回步骤S1;
S4、基于所确定的测井曲线中的曲线形态的具体类型以及所提取的测井曲线的曲线幅度及幅度差、曲线旋回形态、旋回幅度、曲线光滑程度、包络线确定当前环境中沉积相带分布情况。
在本实施例中,所述待测钻井的测井曲线的形态包括钟形、漏斗形、箱形、对称齿形、反向齿形、正向齿形;所述预设的标准测井曲线的标准形态包括标准钟形、标准漏斗形、标准箱形、标准对称齿形、标准反向齿形、标准正向齿形。
在本实施例中,所述步骤S3具体为:
S31、将所提取的测井曲线中的曲线形态与预设的标准测井曲线中所有的标准曲线形态进行一一比对;
S32、当所提取的测井曲线中的曲线形态与预设的标准测井曲线中某种类型的标准曲线形态之间的重合度在预设阈值范围时,则该类型为所提取的测井曲线中的形态的具体类型;
S33、所提取的测井曲线中的曲线形态与预设的标准测井曲线中所有的标准曲线形态之间的重合度都不在预设阈值范围,则返回步骤S1。
在本实施例中,所述预设阈值范围为95%-100%。
最后应说明的是:以上各实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述各实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分或者全部技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。

Claims (8)

1.一种山前地区沉积相带分布范围勘测方法,具体包括以下步骤:
S1、采集获取当前待测钻井的测井曲线;
S2、提取所述当前待测钻井的测井曲线的曲线幅度及幅度差、曲线形态、曲线旋回形态、旋回幅度、曲线光滑程度、包络线;
S3、将所提取的测井曲线中的曲线形态与预设的标准测井曲线的标准曲线形态进行图形比对,当提取的测井曲线中的曲线形态与预设的标准测井曲线的标准曲线形态之间的重合度在预设阈值范围时,则根据预设的标准测井曲线的标准曲线形态确定所提取的测井曲线中的曲线形态的具体类型,当提取的测井曲线中的曲线形态与预设的标准测井曲线的标准曲线形态之间的重合度不在预设阈值范围时,则返回步骤S1;
S4、基于所确定的测井曲线中的曲线形态的具体类型以及所提取的测井曲线的曲线幅度及幅度差、曲线旋回形态、旋回幅度、曲线光滑程度、包络线确定当前环境中沉积相带分布情况。
2.根据权利要求1所述的山前地区沉积相带分布范围勘测方法,其特征在于:所述待测钻井的测井曲线的曲线形态类型包括钟形、漏斗形、箱形、对称齿形、反向齿形、正向齿形;所述预设的标准测井曲线的标准曲线形态类型包括标准钟形、标准漏斗形、标准箱形、标准对称齿形、标准反向齿形、标准正向齿形。
3.根据权利要求1所述的山前地区沉积相带分布范围勘测方法,其特征在于:所述预设阈值范围为95%-100%。
4.根据权利要求1所述的山前地区沉积相带分布范围勘测方法,其特征在于:所述测井曲线为自然电位测井曲线、自然伽马曲线、微球电阻率曲线或侧向电阻率曲线。
5.根据权利要求1所述的山前地区沉积相带分布范围勘测方法,其特征在于:所述步骤S3具体包括步骤:
S31、工作人员根据所提取的测井曲线中的形态初步判断出该形态所代表的类型;
S32、然后将初步得出的测井曲线的曲线形态类型与预设的标准测井曲线中对应同类型的标准曲线形态进行比对;
S33、当初步得出的测井曲线的形态类型与预设的标准测井曲线中对应同类型的标准形态之间的重合度在预设阈值范围时,则根据预设的标准测井曲线的标准形态确定所提取的测井曲线中的形态的具体类型;
S34、当初步得出的测井曲线的形态类型与预设的标准测井曲线中对应同类型的标准形态之间的重合度不在预设阈值范围,则将初步得出的测井曲线的形态类型与预设的标准测井曲线中其它类型的标准形态进行比对,当初步得出的测井曲线的形态类型与预设的标准测井曲线中其它类型的标准形态之间的重合度在预设阈值范围时,则根据预设的标准测井曲线的标准形态确定所提取的测井曲线中的形态的具体类型;
S35、若初步得出的测井曲线的曲线形态类型与预设的标准测井曲线中对所有类型的标准曲线形态之间的重合度都不在预设阈值范围,则返回步骤S1。
6.根据权利要求1所述的山前地区沉积相带分布范围勘测方法,其特征在于:所述步骤S3具体包括步骤:
S31、将所提取的测井曲线中的曲线形态与预设的标准测井曲线中所有的标准曲线形态进行一一比对;
S32、当所提取的测井曲线中的曲线形态与预设的标准测井曲线中某种类型的标准曲线形态之间的重合度在预设阈值范围时,则该类型为所提取的测井曲线中的形态的具体类型;
S33、所提取的测井曲线中的曲线形态与预设的标准测井曲线中所有的标准曲线形态之间的重合度都不在预设阈值范围,则返回步骤S1。
7.根据权利要求5所述的山前地区沉积相带分布范围勘测方法,其特征在于:所述预设阈值范围为98%-100%。
8.根据权利要求6所述的山前地区沉积相带分布范围勘测方法,其特征在于:所述预设阈值范围为95%-100%。
CN201911286459.1A 2019-12-13 2019-12-13 一种山前地区沉积相带分布范围勘测方法 Active CN111061892B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911286459.1A CN111061892B (zh) 2019-12-13 2019-12-13 一种山前地区沉积相带分布范围勘测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911286459.1A CN111061892B (zh) 2019-12-13 2019-12-13 一种山前地区沉积相带分布范围勘测方法

Publications (2)

Publication Number Publication Date
CN111061892A CN111061892A (zh) 2020-04-24
CN111061892B true CN111061892B (zh) 2023-05-23

Family

ID=70301630

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911286459.1A Active CN111061892B (zh) 2019-12-13 2019-12-13 一种山前地区沉积相带分布范围勘测方法

Country Status (1)

Country Link
CN (1) CN111061892B (zh)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103510947A (zh) * 2012-06-21 2014-01-15 中国石油化工股份有限公司 滩坝砂岩微相识别模式的建立及其在微相识别中的应用
CN103969682A (zh) * 2013-01-28 2014-08-06 中国石油集团东方地球物理勘探有限责任公司 一种钻测井资料深度匹配方法及系统
US9097096B1 (en) * 2010-09-10 2015-08-04 Selman and Associates, Ltd. Multi dimensional model for directional drilling
CN105484739A (zh) * 2015-11-26 2016-04-13 中国科学院武汉岩土力学研究所 碳酸盐岩地层孔隙压力测试方法和装置
CN106199755A (zh) * 2016-06-24 2016-12-07 北京金双狐油气技术有限公司 点线面三维度综合沉积相的研究方法
CN109725348A (zh) * 2017-10-30 2019-05-07 中国石油化工股份有限公司 一种基于地震资料识别沉积相的方法
CN110208872A (zh) * 2019-06-05 2019-09-06 中国石油大港油田勘探开发研究院 一种冲积扇隔夹层综合识别方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11814931B2 (en) * 2018-05-31 2023-11-14 Schlumberger Technology Corporation Machine-learning for sedimentary facies prediction

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9097096B1 (en) * 2010-09-10 2015-08-04 Selman and Associates, Ltd. Multi dimensional model for directional drilling
CN103510947A (zh) * 2012-06-21 2014-01-15 中国石油化工股份有限公司 滩坝砂岩微相识别模式的建立及其在微相识别中的应用
CN103969682A (zh) * 2013-01-28 2014-08-06 中国石油集团东方地球物理勘探有限责任公司 一种钻测井资料深度匹配方法及系统
CN105484739A (zh) * 2015-11-26 2016-04-13 中国科学院武汉岩土力学研究所 碳酸盐岩地层孔隙压力测试方法和装置
CN106199755A (zh) * 2016-06-24 2016-12-07 北京金双狐油气技术有限公司 点线面三维度综合沉积相的研究方法
CN109725348A (zh) * 2017-10-30 2019-05-07 中国石油化工股份有限公司 一种基于地震资料识别沉积相的方法
CN110208872A (zh) * 2019-06-05 2019-09-06 中国石油大港油田勘探开发研究院 一种冲积扇隔夹层综合识别方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
《测井曲线的沉积相自动识别》;高海焦;《中国优秀硕士学位论文全文数据库-工程科技I辑》;20150315;全文 *

Also Published As

Publication number Publication date
CN111061892A (zh) 2020-04-24

Similar Documents

Publication Publication Date Title
US10324229B2 (en) System and method of pore type classification for petrophysical rock typing
CN103454678B (zh) 一种地震切片等时性的确定方法及系统
CN107402176B (zh) 裂缝孔隙度的确定方法和装置
CN103993871A (zh) 针对薄互层地层的测井资料标准化处理方法及装置
CN107908919B (zh) 一种致密油储层层理缝有效性的评价方法
CN105388525A (zh) 储层预测方法及装置
CN105114067A (zh) 一种岩性测井相的方法
CN102562051B (zh) 处理气测图版原始数据的方法
US10152783B2 (en) Method and system for two dimensional sediment particle shape characterization based on borehole image
CN104459800A (zh) 一种预测砂体尖灭的方法和装置
US20220120933A1 (en) Method of detection of hydrocarbon horizontal slippage passages
CN103628866A (zh) 获取地层的地层压力系数的方法
CN111061892B (zh) 一种山前地区沉积相带分布范围勘测方法
AU2017279838B1 (en) Method for classifying deep rock geofacies based on data mining
Holland et al. Evolution of fractures in a highly dynamic thermal, hydraulic, and mechanical system-(II) Remote sensing fracture analysis, Jabal Shams, Oman Mountains
CN108181655B (zh) 暗河垮塌体系形态的确定方法和装置
CN106842289B (zh) 一种适用于测井约束反演的波阻抗曲线去压实处理方法
CN109738955B (zh) 一种基于成分-结构分类下的变质岩岩性综合判别方法
Ross et al. An algorithm for automated identification of fault zone trapped waves
Koch et al. Does fluvial channel-belt clustering predict net sand to gross rock volume? Architectural metrics and point-pattern analysis of a digital outcrop model
CN111596351B (zh) 碳酸盐岩输导体系定量评价方法、系统、装置及存储介质
CN109598049B (zh) 钻孔岩石裂隙发育程度及区域岩石裂隙发育规律的方法
CN107942400B (zh) 一种预测钻井风险的方法及装置
US11320550B2 (en) Processing a 4D seismic signal based on noise model
CN107831537B (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