CN113250685A - 一种页岩油藏全尺度升级方法 - Google Patents

一种页岩油藏全尺度升级方法 Download PDF

Info

Publication number
CN113250685A
CN113250685A CN202110562214.8A CN202110562214A CN113250685A CN 113250685 A CN113250685 A CN 113250685A CN 202110562214 A CN202110562214 A CN 202110562214A CN 113250685 A CN113250685 A CN 113250685A
Authority
CN
China
Prior art keywords
scale
oil reservoir
shale oil
flow
shale
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
CN202110562214.8A
Other languages
English (en)
Other versions
CN113250685B (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 University of Petroleum East China
Original Assignee
China University of Petroleum East China
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 University of Petroleum East China filed Critical China University of Petroleum East China
Priority to CN202110562214.8A priority Critical patent/CN113250685B/zh
Publication of CN113250685A publication Critical patent/CN113250685A/zh
Application granted granted Critical
Publication of CN113250685B publication Critical patent/CN113250685B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • EFIXED CONSTRUCTIONS
    • E21EARTH DRILLING; MINING
    • E21BEARTH DRILLING, e.g. DEEP 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
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/28Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]

Abstract

本发明公开了一种页岩油藏全尺度升级方法,包括:S1.构建页岩油藏多尺度数字岩心;S2.采用滑移边界的N‑S方程描述页岩油藏孔隙尺度上的流动;S3.建立页岩油藏达西尺度上的有机质和无机质流动方程;S4.建立页岩油藏宏观等效流动模型,及等效参数,完成页岩油藏全尺度升级。本发明通过实验、分子模拟和两次基于均化理论的尺度升级,解决了页岩油藏微观流动机理与宏观流动模型间的尺度关联难题,实现由传统渗流力学向非常规油气渗流力学的跨越发展,为提高增产效果、保持页岩油藏长期稳定和高效生产提供理论依据。

Description

一种页岩油藏全尺度升级方法
技术领域
本发明属于油气田开发工程油藏数值模拟技术领域,特别是涉及一种页岩油藏全尺度升级方法。
背景技术
页岩油藏是重要的非常规油气资源,具有分布广、储量大的特点,与常规油气藏相比,页岩油藏赋存方式多样,孔隙结构复杂且具有多尺度特征,既具有纳米级的基质孔隙、又具有微米-厘米尺度的天然裂缝系统,以及大规模压裂所产生的米到百米级的人工裂缝。不同的多孔介质特征尺度具有不同的流动机制、控制方程、渗流模型以及流动模拟方法,传统的油气渗流理论以及单一尺度上的油藏渗流理论和流动模拟方法无法应用于页岩油藏,因此有必要揭示页岩油藏多尺度孔隙结构特征以及建立从孔隙尺度-达西尺度-宏观尺度的页岩油藏全尺度升级的渗流模拟理论与方法,将页岩油藏中的微观分布特征有效地表征到页岩油藏宏观开发动态中,为页岩油藏的高效开发提供理论指导。
发明内容
本发明提出一种页岩油藏全尺度升级方法,通过实验、分子模拟与尺度升级相结合,建立了一种从孔隙尺度-达西尺度-宏观尺度的页岩油藏全尺度渗流模拟理论与方法,揭示页岩油藏多尺度孔隙结构特征,解决页岩油藏微观流动机理与宏观流动模型间的尺度关联难题。
为实现上述目的,本发明提供了如下方案:
一种页岩油藏全尺度升级方法,包括如下步骤:
S1.通过不同方式对页岩油藏岩样开展多级扫描,构建页岩油藏多尺度数字岩心;
S2.基于所述页岩油藏多尺度数字岩心和页岩油藏微观孔隙结构,建立不同结构纳米级规则孔隙内流体流动模型,得到不同结构纳米孔隙内流体的滑移长度和有效粘度参数,采用滑移边界的N-S方程描述页岩油藏孔隙尺度上的流动;
S3.升级所述页岩油藏孔隙尺度上的流动,建立页岩油藏达西尺度上的有机质和无机质流动方程;
S4.升级所述有机质和无机质流动方程,建立页岩油藏宏观流动模型,完成页岩油藏全尺度升级。
优选的,所述S1中构建页岩油藏多尺度数字岩心包括如下步骤:
S1.1.分别建立微米尺度数字岩心、亚微米尺度数字岩心、纳米尺度数字岩心;
S1.2.基于所述微米尺度数字岩心、所述亚微米尺度数字岩心、所述纳米尺度数字岩心,构建所述页岩油藏多尺度数字岩心。
优选的,所述S1.1中:
通过获取所述页岩油藏岩样的微米级岩心裂缝形态和基质分布特征的投影,建立所述微米尺度数字岩心;
通过获取所述页岩油藏岩样的亚微米级孔隙的几何和拓扑结构特征,建立所述亚微米尺度数字岩心;
通过获取所述页岩油藏岩样的有机质孔隙特征,建立所述纳米尺度数字岩心。
优选的,所述S1.1中:
通过对所述页岩油藏岩样进行CT扫描,获取所述微米级岩心裂缝形态和基质分布特征的投影;
通过对所述页岩油藏岩样进行SEM-MAP扫描实验,获取所述亚微米级孔隙的几何和拓扑结构特征;
通过对所述页岩油藏岩样进行SEM纳米扫描,获取所述有机质孔隙特征。
优选的,所述S2中,建立不同结构纳米级规则孔隙内流体流动模型的方法为:基于所述页岩油藏多尺度数字岩心和页岩油藏微观孔隙结构,使用分子模拟方法,建立不同结构的所述纳米级规则孔隙内流体流动模型。
优选的,所述S3中将所述页岩油藏孔隙尺度上的流动从孔隙尺度升级到达西尺度的方法为:基于所述页岩油藏达西尺度上的有机质和无机质流动方程,采用均化理论,建立所述页岩油藏在达西尺度上的有机质和无机质流动方程。
优选的,所述S4中,所述页岩油藏宏观流动模型包括微观流动机制和孔隙结构特征。
优选的,所述S4中升级所述有机质和无机质流动方程的方法为:基于有机质和无机质的分布特征和流动机理,采用均化理论将所述有机质和无机质流动方程从达西尺度升级到宏观尺度,建立所述页岩油藏宏观等效流动模型,及等效参数,完成页岩油藏多级尺度升级。
本发明的有益效果为:
本发明公开了一种页岩油藏全尺度升级方法,通过实验、分子模拟和两次基于均化理论的尺度升级,解决页岩油藏微观流动机理与宏观流动模型间的尺度关联难题,揭示了页岩油藏达西尺度上的流动规律,建立考虑微观流动机制及孔隙结构特征的页岩油藏宏观流动模型,本发明突破了传统油气渗流理论的局限,实现由传统渗流力学向非常规油气渗流力学的跨越发展,为提高增产效果、保持页岩油藏长期稳定和高效生产提供理论依据。
附图说明
为了更清楚地说明本发明的技术方案,下面对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例的全尺度升级方法流程示意图;
图2为本发明实施例页岩油藏尺度升级流程图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
如图1所示,为本发明提供的页岩油藏全尺度升级方法流程示意图,包括如下步骤:
S1.通过不同方式对页岩油藏岩样开展多级扫描,构建页岩油藏多尺度数字岩心;具体步骤包括:
对页岩油藏岩样进行不同分辨率下的CT扫描,获取大尺度的微米级的岩心裂缝形态和基质分布特征的投影,建立微米尺度数字岩心;
对页岩油藏岩样进行SEM-MAP扫描实验,获取得到页岩油藏岩样亚微米级孔隙的几何和拓扑结构特征,基于有机质和无机质孔隙,构建亚微米尺度数字岩心;
对页岩油藏岩样进行SEM高精度纳米扫描,基于有机质孔隙的特征,建立纳米尺度数字岩心;
基于上述微米尺度数字岩心、亚微米尺度数字岩心、纳米尺度数字岩心,针对页岩油藏储层孔隙结构多尺度特点,精确的刻画页岩油藏不同尺度下的孔隙结构,构建的页岩油藏多尺度数字岩心。
S2.基于页岩油藏多尺度数字岩心和页岩油藏微观孔隙结构,采用分子模拟方法研究不同结构的纳米级规则孔隙内流体流动,建立不同结构纳米级规则孔隙内流体流动模型,得到不同结构纳米孔隙内流体的滑移长度和有效粘度等参数,采用滑移边界的N-S方程描述页岩油藏孔隙尺度上的流动。
S3.针对采用滑移边界的N-S方程描述页岩油藏孔隙尺度上的流动,采用均化理论,将页岩油藏从孔隙尺度升级到达西尺度,建立页岩油藏达西尺度上的有机质和无机质流动方程,同时也揭示了页岩油藏达西尺度上的流动规律。
S4.针对S3得到的页岩油藏达西尺度上的有机质及无机质流动方程,基于有机质和无机质的分布特征及不同的流动机理,再次基于均化理论从达西尺度升级到宏观尺度,建立包含微观流动机制及孔隙结构特征的页岩油藏宏观等效流动模型,及相关等效参数,完成页岩油藏全尺度升级。
如图2所示,为本发明实施例两次尺度升级概述,即从孔隙尺度升级到达西尺度,再升级到宏观尺度。
本发明通过实验、分子模拟和两次基于均化理论的尺度升级,解决页岩油藏微观流动机理与宏观流动模型间的尺度关联难题,揭示了页岩油藏达西尺度上的流动规律,建立考虑微观流动机制及孔隙结构特征的页岩油藏宏观流动模型,实现由传统渗流力学向非常规油气渗流力学的跨越发展,为提高增产效果、保持页岩油藏长期稳定和高效生产提供理论依据。
以上所述的实施例仅是对本发明优选方式进行的描述,并非对本发明的范围进行限定,在不脱离本发明设计精神的前提下,本领域普通技术人员对本发明的技术方案做出的各种变形和改进,均应落入本发明权利要求书确定的保护范围内。

Claims (8)

1.一种页岩油藏全尺度升级方法,其特征在于,包括如下步骤:
S1.通过对页岩油藏岩样开展多级分辨率扫描,构建页岩油藏多尺度数字岩心;
S2.基于所述页岩油藏多尺度数字岩心和页岩油藏微观孔隙结构,建立不同结构纳米级规则孔隙内流体流动模型,得到不同结构纳米孔隙内流体的滑移长度和有效粘度参数,采用滑移边界的N-S方程描述页岩油藏孔隙尺度上的流动;
S3.升级所述页岩油藏孔隙尺度上的流动,建立页岩油藏达西尺度上的有机质和无机质流动方程;
S4.升级所述有机质和无机质流动方程,建立页岩油藏宏观流动模型,完成页岩油藏全尺度升级。
2.根据权利要求1所述的页岩油藏全尺度升级方法,其特征在于,所述S1中构建页岩油藏多尺度数字岩心包括如下步骤:
S1.1.分别建立微米尺度数字岩心、亚微米尺度数字岩心、纳米尺度数字岩心;
S1.2.基于所述微米尺度数字岩心、所述亚微米尺度数字岩心、所述纳米尺度数字岩心,构建所述页岩油藏多尺度数字岩心。
3.根据权利要求2所述的页岩油藏全尺度升级方法,其特征在于,所述S1.1中:
通过获取所述页岩油藏岩样的微米级岩心裂缝形态和基质分布特征的投影,建立所述微米尺度数字岩心;
通过获取所述页岩油藏岩样的亚微米级孔隙的几何和拓扑结构特征,建立所述亚微米尺度数字岩心;
通过获取所述页岩油藏岩样的有机质孔隙特征,建立所述纳米尺度数字岩心。
4.根据权利要求3所述的页岩油藏全尺度升级方法,其特征在于,所述S1.1中:
通过对所述页岩油藏岩样进行CT扫描,获取所述微米级岩心裂缝形态和基质分布特征的投影;
通过对所述页岩油藏岩样进行SEM-MAP扫描实验,获取所述亚微米级孔隙的几何和拓扑结构特征;
通过对所述页岩油藏岩样进行SEM纳米扫描,获取所述有机质孔隙特征。
5.根据权利要求1所述的页岩油藏全尺度升级方法,其特征在于,所述S2中,建立不同结构纳米级规则孔隙内流体流动模型的方法为:基于所述页岩油藏多尺度数字岩心和页岩油藏微观孔隙结构,采用分子模拟方法,建立不同结构的所述纳米级规则孔隙内流体流动模型。
6.根据权利要求5所述的页岩油藏全尺度升级方法,其特征在于,所述S3中将所述页岩油藏孔隙尺度上的流动从孔隙尺度升级到达西尺度的方法为:基于所述页岩油藏达西尺度上的有机质和无机质流动方程,采用均化理论,建立所述页岩油藏在达西尺度上的有机质和无机质流动方程。
7.根据权利要求1所述的页岩油藏全尺度升级方法,其特征在于,所述S4中,所述页岩油藏宏观流动模型包括微观流动机制和孔隙结构特征。
8.根据权利要求7所述的页岩油藏全尺度升级方法,其特征在于,所述S4中升级所述有机质和无机质流动方程的方法为:基于有机质和无机质的分布特征和流动机理,采用均化理论将所述有机质和无机质流动方程从达西尺度升级到宏观尺度,建立所述页岩油藏宏观等效流动模型,及等效参数,完成页岩油藏全尺度升级。
CN202110562214.8A 2021-05-24 2021-05-24 一种页岩油藏全尺度升级方法 Active CN113250685B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110562214.8A CN113250685B (zh) 2021-05-24 2021-05-24 一种页岩油藏全尺度升级方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110562214.8A CN113250685B (zh) 2021-05-24 2021-05-24 一种页岩油藏全尺度升级方法

Publications (2)

Publication Number Publication Date
CN113250685A true CN113250685A (zh) 2021-08-13
CN113250685B CN113250685B (zh) 2022-04-29

Family

ID=77183820

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110562214.8A Active CN113250685B (zh) 2021-05-24 2021-05-24 一种页岩油藏全尺度升级方法

Country Status (1)

Country Link
CN (1) CN113250685B (zh)

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140305649A1 (en) * 2012-04-12 2014-10-16 Hongxin Tang Submicron particles surfactant method for improved oil recovery from subterranean reservoirs
CN106716114A (zh) * 2015-09-11 2017-05-24 数岩科技(厦门)股份有限公司 多孔介质分析系统和方法
CN107060746A (zh) * 2017-04-27 2017-08-18 中国石油大学(华东) 一种复杂裂缝性油藏流动模拟的方法
US20170298714A1 (en) * 2016-04-18 2017-10-19 International Business Machines Corporation Method, apparatus and computer program product providing simulator for enhanced oil recovery based on micron and submicron scale fluid-solid interactions
WO2018204463A1 (en) * 2017-05-03 2018-11-08 Schlumberger Technology Corporation Fractured reservoir simulation
CN110702581A (zh) * 2019-10-23 2020-01-17 山东省科学院海洋仪器仪表研究所 一种强非均质多孔介质的多尺度渗透率计算方法
CN111428321A (zh) * 2020-04-03 2020-07-17 中国石油天然气股份有限公司 一种基于简化数字岩心的砾岩储层孔隙网络模型建模方法
CN111553108A (zh) * 2020-05-20 2020-08-18 中国石油大学(华东) 一种页岩气藏流固耦合多尺度数值模拟方法
CN111597721A (zh) * 2020-05-20 2020-08-28 中国石油大学(华东) 一种基于均匀化理论的页岩基质流固耦合尺度升级方法
CN112199903A (zh) * 2020-09-25 2021-01-08 中国地质大学(武汉) 基于多参数的离散元纳米颗粒封堵页岩孔隙数值模拟方法

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140305649A1 (en) * 2012-04-12 2014-10-16 Hongxin Tang Submicron particles surfactant method for improved oil recovery from subterranean reservoirs
CN106716114A (zh) * 2015-09-11 2017-05-24 数岩科技(厦门)股份有限公司 多孔介质分析系统和方法
US20170298714A1 (en) * 2016-04-18 2017-10-19 International Business Machines Corporation Method, apparatus and computer program product providing simulator for enhanced oil recovery based on micron and submicron scale fluid-solid interactions
CN107060746A (zh) * 2017-04-27 2017-08-18 中国石油大学(华东) 一种复杂裂缝性油藏流动模拟的方法
WO2018204463A1 (en) * 2017-05-03 2018-11-08 Schlumberger Technology Corporation Fractured reservoir simulation
CN110702581A (zh) * 2019-10-23 2020-01-17 山东省科学院海洋仪器仪表研究所 一种强非均质多孔介质的多尺度渗透率计算方法
CN111428321A (zh) * 2020-04-03 2020-07-17 中国石油天然气股份有限公司 一种基于简化数字岩心的砾岩储层孔隙网络模型建模方法
CN111553108A (zh) * 2020-05-20 2020-08-18 中国石油大学(华东) 一种页岩气藏流固耦合多尺度数值模拟方法
CN111597721A (zh) * 2020-05-20 2020-08-28 中国石油大学(华东) 一种基于均匀化理论的页岩基质流固耦合尺度升级方法
CN112199903A (zh) * 2020-09-25 2021-01-08 中国地质大学(武汉) 基于多参数的离散元纳米颗粒封堵页岩孔隙数值模拟方法

Non-Patent Citations (16)

* Cited by examiner, † Cited by third party
Title
FAN WEIPENG ET AL.: "Homogenization approach for liquid flow within shale system considering slip effect", 《JOURNAL OF CLEANER PRODUCTION》 *
MINZHANG ET AL.: "Triple-continuum modeling of shale gas reservoirs considering the effect of kerogen", 《JOURNAL OF NATURAL GAS SCIENCE AND ENGINEERING》 *
SUN HAI ET AL.: "Characterization of gas transport behaviors in shale gas and tight gas reservoirs by digital rock analysis", 《INTERNATIONAL JOURNAL OF HEAT AND MASS TRANSFER》 *
VISWASANTHI CHANDRA ET AL.: "Improved Reservoir Characterization Through Rapid Visualization and Analysis of Multiscale Image Data Using A Digital Core Analysis Ecosystem", 《 THE ABU DHABI INTERNATIONAL PETROLEUM EXHIBITION & CONFERENCE》 *
WANG DONGYING ET AL.: "Image-based core-scale real gas apparent permeability from pore-scale experimental data in shale reservoirs", 《FUEL》 *
WEIPENGFAN ET AL.: "An upscaled transport model for shale gas considering multiple mechanisms and heterogeneity based on homogenization theory", 《JOURNAL OF PETROLEUM SCIENCE AND ENGINEERING》 *
YAN XIA ET AL.: "Hierarchical Modeling of Hydromechanical Coupling in Fractured Shale Gas Reservoirs with Multiple Porosity Scales", 《ENERGY & FUELS》 *
Z JIANG ET AL.: "Representation of multiscale heterogeneity via multiscale pore networks", 《WATER RESOURCES RESEARCH》 *
严侠等: "基于均匀化理论的页岩基质流固耦合尺度升级", 《中国石油大学学报(自然科学版)》 *
姚军等: "现代油气渗流力学体系及其发展趋势", 《科学通报》 *
姚军等: "页岩气藏开发中的关键力学问题", 《中国科学:物理学 力学 天文学》 *
姜黎明等: "基于数字岩心的天然气储层岩石声电特性数值模拟研究", 《中国优秀博士学位论文全文数据库基础科学辑》 *
孙海等: "基于均匀化理论的页岩基岩运移机制尺度升级研究", 《中国科学:物理学 力学 天文学》 *
张庆福等: "多尺度嵌入式离散裂缝模型模拟方法", 《计算力学学报》 *
符东宇等: "基于格子Boltzmann方法的页岩气藏微观渗流机理研究", 《中国优秀硕士学位论文全文数据库工程科技Ⅰ辑》 *
胡洋铭等: "基于数字岩芯的碳酸盐岩弹性性质模拟", 《中国优秀硕士学位论文全文数据库基础科学辑》 *

Also Published As

Publication number Publication date
CN113250685B (zh) 2022-04-29

Similar Documents

Publication Publication Date Title
Yang et al. Pore-scale simulation of shale oil flow based on pore network model
CN111553108B (zh) 一种页岩气藏流固耦合多尺度数值模拟方法
CN111428321B (zh) 一种基于简化数字岩心的砾岩储层孔隙网络模型建模方法
CN105160134B (zh) 致密储层多重介质中油气流动的混合介质模拟方法及装置
CN107622165A (zh) 一种页岩气水平井重复压裂产能计算方法
CN111236931A (zh) 一种气藏大斜度井非稳态试井图版生成方法和系统
CN111428426A (zh) 一种基于格子Boltzmann的页岩气多相流动模拟方法
CN111274528B (zh) 储层裂缝渗吸质量预测方法及系统
CN105869496A (zh) 一种可视化微观孔隙结构仿真物理模型及制作方法
CN111062129A (zh) 页岩油复杂缝网离散裂缝连续介质混合数值模拟方法
CN111980654A (zh) 一种非均匀页岩油藏分段压裂水平井产能计算方法
CN112360422B (zh) 一种页岩气藏压裂水平井产量预测方法及系统
CN109829217A (zh) 压裂性裂缝油藏产能模拟方法及装置
CN111980695A (zh) 基于不同占比的页岩气藏单相气体双孔双渗模型构建方法
CN113250685B (zh) 一种页岩油藏全尺度升级方法
Shen et al. Spontaneous imbibition in asymmetric branch-like throat structures in unconventional reservoirs
Svec et al. Flow simulation of fiber reinforced self compacting concrete using Lattice Boltzmann method
Zhao et al. Laws of gas and water flow and mechanism of reservoir drying in tight sandstone gas reservoirs
CN116858630B (zh) 一种真实页岩干酪根孔隙模型的构建方法及系统
CN113297777A (zh) 碳酸盐岩油气藏酸化反应流多尺度数值模拟方法及系统
Li et al. Similar structure of solution for triple media shale gas reservoir
CN109723420A (zh) 一种改造高灰质致密油气藏的组合缝网压裂方法
CN109025983B (zh) 一种模拟致密油藏微观模型制作方法
Xu et al. Numerical study of water–air distribution in unsaturated soil by using lattice Boltzmann method
CN117408191B (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