CN112489148B - 一种砂岩型铀矿平米铀量等值线图成图方法 - Google Patents

一种砂岩型铀矿平米铀量等值线图成图方法 Download PDF

Info

Publication number
CN112489148B
CN112489148B CN201910857956.6A CN201910857956A CN112489148B CN 112489148 B CN112489148 B CN 112489148B CN 201910857956 A CN201910857956 A CN 201910857956A CN 112489148 B CN112489148 B CN 112489148B
Authority
CN
China
Prior art keywords
uranium
contour
value
amount
minimum
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
CN201910857956.6A
Other languages
English (en)
Other versions
CN112489148A (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.)
Nuke Industry No216 Brigade
Original Assignee
Nuke Industry No216 Brigade
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 Nuke Industry No216 Brigade filed Critical Nuke Industry No216 Brigade
Priority to CN201910857956.6A priority Critical patent/CN112489148B/zh
Publication of CN112489148A publication Critical patent/CN112489148A/zh
Application granted granted Critical
Publication of CN112489148B publication Critical patent/CN112489148B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/002D [Two Dimensional] image generation
    • G06T11/001Texturing; Colouring; Generation of texture or colour
    • 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

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明属于数字勘察系统技术领域,具体涉及一种砂岩型铀矿平米铀量等值线图成图方法,包括五个步骤,第一步,确定纯数字十位数工作区代码,第二步,确定储量估算剖面图图幅图框大小,第三步,确定钻孔不同平米铀量等值线值得大小,第四步,平米铀量等值线区域的形成,第五步,对不同数值平米铀量等值线区域进行颜色填充形成平米铀量等值线图;本发明对基于数字煤矿勘查软件的煤矿地质勘查中砂岩型煤矿平米铀量等值线图编制需手工、繁琐、多次重复,实现平米铀量等值线图编制数字化、标准化、信息化。

Description

一种砂岩型铀矿平米铀量等值线图成图方法
技术领域
本发明属于数字勘察系统技术领域,具体涉及一种砂岩型铀矿平米铀量等值线图成图方法。
背景技术
以往的平米铀量等值线图多为人工手动完成,工作量大而且计算繁琐。随着工作区勘查程度的提高,积累了大量的钻孔数据,为完成平米铀量等值线图,必然要利用前人形成地质资料形成最新、最全的平米铀量等值线图,需要重新整理、制图,极大的增加了工作量,降低了工作效率。随着计算机技术的全面迅速发展,核工业系统研发了数字煤矿勘查软件,建立了煤矿勘查全过程的独立数据库,本专利利用计算机技术,由软件人机交互完成工作区确定、剖面图图幅范围、钻孔定位、计算参数选择、计算方法等复杂工序。实现了全流程自动化、智能化,有效的提高工作效率和成图精度。
发明内容
本发明的目的是针对现有技术存在的上述缺陷,提供一种砂岩型煤矿平米铀量等值线图成图方法,实现了全流程自动化、智能化,有效的提高工作效率和成图精度,为储量计算奠定了基础。
本发明的技术方案如下:
一种砂岩型铀矿平米铀量等值线图成图方法,包括五个步骤,第一步,确定纯数字十位数工作区代码,第二步,确定储量估算剖面图图幅图框大小,第三步,确定钻孔不同平米铀量等值线值得大小,第四步,平米铀量等值线区域的形成,第五步,对不同数值平米铀量等值线区域进行颜色填充形成平米铀量等值线图;
第一步,确定纯数字十位数工作区代码;
选定需要进行地质勘查的区域,绘制该区域矩形勘探线图,该矩形勘探线中包括若干钻孔;并通过所述钻孔确定该勘探线图中的基本勘探线。
一种砂岩型铀矿平米铀量等值线图成图方法,所述第二步,定储量估算剖面图图幅图框大小;
确定勘探线图幅范围,即长度和宽度:
以第一步中确定的矩形勘探线图的一个端点为原点,建立坐标系,则矩形勘探线图中1号钻孔至n号钻孔坐标依次(X1,Y1),(X2,Y2),(X3,Y3,),……(Xn,Yn);
用最小值、最大值公式来计算矩形综合成图中:
所有钻孔的最小XMinZkX=min(Zk1.x,Zk2.x,.......Zkn.x);
所有钻孔的最小YMinZkY=min(Zk1.y,Zk2.y.......Zkn.y);
所有钻孔的最大XMaxZkX=max(Zk1.x,Zk2.x,.......Zkn.x);
所有钻孔的最大YMaxZkY=max(Zk1.y,Zk2.y.......Zkn.y);
所有勘探线的最小XMinKtxX=min(Ktx1.x,Ktx2.x,.......Ktxn.x);
所有勘探线的最小YMinKtxY=min(Ktx1.y,Ktx2.y.......Ktxn.y);
所有勘探线的最大XMaxKtxX=max(Ktx1.x,Ktx2.x,.......Ktxn.x);
所有勘探线的最大YMaxKtxY=max(Ktx1.y,Ktx2.y.......Ktxn.y);
所有等值线的最小XMindzxX=min(dzx1.x,dzx2.x,dzx3.x,.....dzxn.x)
所有等值线的最小YMindzxY=min(dzx1.y,dzx2.y,dzx3.y,.....dzxn.y)
所有等值线的最大XMaxdzxX=max(dzx1.x,dzx2.x,dzx3.x,.....dzxn.x)
所有等值线的最大YMaxdzxY=max(dzx1.y,dzx2.y,dzx3.y,.....dzxn.y)
得到平米铀量等值线图:
图框的最小X=min(MinZkX,MinKtxXMindzxX);
图框的最小Y=min(MinZkY,MinKtxYMindzxY);
图框的最大X=max(MaxZkX,MaxKtxXMaxdzxX);
图框的最大Y=max(MaxZkY,MaxKtxYMaxdzxY);
其中:
zk1.x表示:1号钻孔孔口坐标的x值
ktx1.x表示:1号勘探线起点或终点的坐标的x值
dzx1.x.表示:第1条等值线坐标的x值
zk1.y表示:1号钻孔孔口坐标的y值
ktx1.y表示:1号勘探线起点或终点的坐标的y值
dzx1.y.表示:第1条等值线坐标的y值。
一种砂岩型铀矿平米铀量等值线图成图方法,所述矩形平米铀量等值线图图幅长度=图框的最大Y-图框的最小Y。
一种砂岩型铀矿平米铀量等值线图成图方法,所述矩形平米铀量等值线图图幅宽度=图框的最大X-图框的最小X。
一种砂岩型铀矿平米铀量等值线图成图方法,所述第三步,确定钻孔不同平米铀量等值线值得大小;
在平米铀量等值线图中确定钻孔不同平米铀量等值线值得大小;
分别测量1至n号钻孔的平米铀量,并定义
工业铀矿体:铀品位≥0.01%,平米铀量≥1千克/平方米,
铀矿化体:铀品位≥0.01%,平米铀量<1千克/平方米。
一种砂岩型铀矿平米铀量等值线图成图方法,所述第四步,平米铀量等值线区域的形成;
将zknUm表示n号钻孔平米铀量等值线值≥Um曲线值用光滑曲线连接并封闭成区形成;
将zkn表示n号钻孔平米铀量等值线值<1曲线值用光滑曲线连接并封闭成区形成。
一种砂岩型铀矿平米铀量等值线图成图方法,当m≥1时,zknUm表示n号钻孔平米铀量等值线值,zknUm≥m;
当m<1时,zknUm表示n号钻孔平米铀量等值线值,zknUm<1;
Um表示平米铀量等值线值,zkn表示钻孔号,m表示平米铀量值值数目,
所述第五步,对不同数值平米铀量等值线区域进行颜色填充形成平米铀量等值线图;
对不同数值平米铀量等值线值形成闭合区进行不同颜色填充。
本发明的有益效果在于:
本发明对基于数字煤矿勘查软件的煤矿地质勘查中砂岩型煤矿平米铀量等值线图编制需手工、繁琐、多次重复,实现平米铀量等值线图编制数字化、标准化、信息化。
附图说明
图1为平米铀量等值线图流程图。
具体实施方式
下面结合附图及具体实施例对本发明作进一步详细说明。
一种砂岩型铀矿平米铀量等值线图成图方法,包括五个步骤,第一步,确定纯数字十位数工作区代码,第二步,确定储量估算剖面图图幅图框大小,第三步,确定钻孔不同平米铀量等值线值得大小,第四步,平米铀量等值线区域的形成,第五步,对不同数值平米铀量等值线区域进行颜色填充形成平米铀量等值线图。
第一步、确定纯数字十位数工作区代码;
选定需要进行地质勘查的区域,绘制该区域矩形勘探线图,该矩形勘探线中包括若干钻孔;并通过所述钻孔确定该勘探线图中的基本勘探线;
第二步、定储量估算剖面图图幅图框大小;
确定勘探线图幅范围,即长度和宽度:
以第一步中确定的矩形勘探线图的一个端点为原点,建立坐标系,则矩形勘探线图中1号钻孔至n号钻孔坐标依次(X1,Y1),(X2,Y2),(X3,Y3,),……(Xn,Yn)。
用最小值、最大值公式来计算矩形综合成图中:
所有钻孔的最小X MinZkX=min(Zk1.x,Zk2.x,.......Zkn.x);
所有钻孔的最小Y MinZkY=min(Zk1.y,Zk2.y.......Zkn.y);
所有钻孔的最大X MaxZkX=max(Zk1.x,Zk2.x,.......Zkn.x);
所有钻孔的最大Y MaxZkY=max(Zk1.y,Zk2.y.......Zkn.y);
所有勘探线的最小X MinKtxX=min(Ktx1.x,Ktx2.x,.......Ktxn.x);
所有勘探线的最小Y MinKtxY=min(Ktx1.y,Ktx2.y.......Ktxn.y);
所有勘探线的最大X MaxKtxX=max(Ktx1.x,Ktx2.x,.......Ktxn.x);
所有勘探线的最大Y MaxKtxY=max(Ktx1.y,Ktx2.y.......Ktxn.y);
所有等值线的最小X MindzxX=min(dzx1.x,dzx2.x,dzx3.x,.....dzxn.x)
所有等值线的最小Y MindzxY=min(dzx1.y,dzx2.y,dzx3.y,.....dzxn.y)
所有等值线的最大X MaxdzxX=max(dzx1.x,dzx2.x,dzx3.x,.....dzxn.x)
所有等值线的最大Y MaxdzxY=max(dzx1.y,dzx2.y,dzx3.y,.....dzxn.y)
得到平米铀量等值线图:
图框的最小X=min(MinZkX,MinKtxX MindzxX);
图框的最小Y=min(MinZkY,MinKtxY MindzxY);
图框的最大X=max(MaxZkX,MaxKtxX MaxdzxX);
图框的最大Y=max(MaxZkY,MaxKtxY MaxdzxY);
矩形平米铀量等值线图图幅长度=图框的最大Y-图框的最小Y。
矩形平米铀量等值线图图幅宽度=图框的最大X-图框的最小X。
其中:
zk1.x表示:1号钻孔孔口坐标的x值
ktx1.x表示:1号勘探线起点或终点的坐标的x值
dzx1.x.表示:第1条等值线坐标的x值
zk1.y表示:1号钻孔孔口坐标的y值
ktx1.y表示:1号勘探线起点或终点的坐标的y值
dzx1.y.表示:第1条等值线坐标的y值。
第三步,确定钻孔不同平米铀量等值线值得大小;
在平米铀量等值线图中确定钻孔不同平米铀量等值线值得大小。
分别测量1至n号钻孔的平米铀量,并定义
工业铀矿体:铀品位≥0.01%,平米铀量≥1千克/平方米,
铀矿化体:铀品位≥0.01%,平米铀量<1千克/平方米。
第四步、平米铀量等值线区域的形成
将zknUm表示n号钻孔平米铀量等值线值≥Um曲线值用光滑曲线连接并封闭成区形成。
将zkn表示n号钻孔平米铀量等值线值<1曲线值用光滑曲线连接并封闭成区形成。
其中:当m≥1时zknUm表示n号钻孔平米铀量等值线值zknUm≥m;
当m<1时zknUm表示n号钻孔平米铀量等值线值zknUm<1;
Um表示平米铀量等值线值,zkn表示钻孔号,m表示平米铀量值值数目。
第五步、对不同数值平米铀量等值线区域进行颜色填充形成平米铀量等值线图
对不同数值平米铀量等值线值形成闭合区进行不同颜色填充。
本发明的效果在于提高工作效率、图件为数字化资料,可反复调用、修改并利用资料数据。

Claims (5)

1.一种砂岩型铀矿平米铀量等值线图成图方法,包括五个步骤,其特征在于:第一步,确定纯数字十位数工作区代码,第二步,确定储量估算剖面图图幅图框大小,第三步,确定钻孔不同平米铀量等值线值得大小,第四步,平米铀量等值线区域的形成,第五步,对不同数值平米铀量等值线区域进行颜色填充形成平米铀量等值线图;
所述第一步,确定纯数字十位数工作区代码包括选定需要进行地质勘查的区域,绘制该区域矩形勘探线图,该矩形勘探线中包括若干钻孔;并通过所述钻孔确定该勘探线图中的基本勘探线;
所述第二步,确定储量估算剖面图图幅图框大小;
确定勘探线图幅范围,即长度和宽度:
以第一步中确定的矩形勘探线图的一个端点为原点,建立坐标系,则矩形勘探线图中1号钻孔至n号钻孔坐标依次(X1,Y1),(X2,Y2),(X3,Y3,),……(Xn,Yn);
用最小值、最大值公式来计算矩形综合成图中:
所有钻孔的最小XMinZkX=min(Zk1.x,Zk2.x,.......Zkn.x);
所有钻孔的最小YMinZkY=min(Zk1.y,Zk2.y.......Zkn.y);
所有钻孔的最大XMaxZkX=max(Zk1.x,Zk2.x,.......Zkn.x);
所有钻孔的最大YMaxZkY=max(Zk1.y,Zk2.y.......Zkn.y);
所有勘探线的最小XMinKtxX=min(Ktx1.x,Ktx2.x,.......Ktxn.x);
所有勘探线的最小YMinKtxY=min(Ktx1.y,Ktx2.y.......Ktxn.y);
所有勘探线的最大XMaxKtxX=max(Ktx1.x,Ktx2.x,.......Ktxn.x);
所有勘探线的最大YMaxKtxY=max(Ktx1.y,Ktx2.y.......Ktxn.y);
所有等值线的最小XMindzxX=min(dzx1.x,dzx2.x,dzx3.x,.....dzxn.x)
所有等值线的最小YMindzxY=min(dzx1.y,dzx2.y,dzx3.y,.....dzxn.y)
所有等值线的最大XMaxdzxX=max(dzx1.x,dzx2.x,dzx3.x,.....dzxn.x)
所有等值线的最大YMaxdzxY=max(dzx1.y,dzx2.y,dzx3.y,.....dzxn.y)得到平米铀量等值线图:
图框的最小X=min(MinZkX,MinKtxXMindzxX);
图框的最小Y=min(MinZkY,MinKtxYMindzxY);
图框的最大X=max(MaxZkX,MaxKtxXMaxdzxX);
图框的最大Y=max(MaxZkY,MaxKtxYMaxdzxY);
其中:
zk1.x表示:1号钻孔孔口坐标的x值
ktx1.x表示:1号勘探线起点或终点的坐标的x值
dzx1.x.表示:第1条等值线坐标的x值
zk1.y表示:1号钻孔孔口坐标的y值
ktx1.y表示:1号勘探线起点或终点的坐标的y值
dzx1.y.表示:第1条等值线坐标的y值;
所述第三步,确定钻孔不同平米铀量等值线值得大小,具体为在平米铀量等值线图中确定钻孔不同平米铀量等值线值得大小;
分别测量1至n号钻孔的平米铀量,并定义
工业铀矿体:铀品位≥0.01%,平米铀量≥1千克/平方米,
铀矿化体:铀品位≥0.01%,平米铀量<1千克/平方米;
所述第四步,平米铀量等值线区域的形成;具体为将zknUm表示n号钻孔平米铀量等值线值≥Um曲线值用光滑曲线连接并封闭成区形成;
将zkn表示n号钻孔平米铀量等值线值<1曲线值用光滑曲线连接并封闭成区形成。
2.如权利要求1所述的一种砂岩型铀矿平米铀量等值线图成图方法,其特征在于:所述平米铀量等值线图图幅长度=图框的最大Y-图框的最小Y。
3.如权利要求1所述的一种砂岩型铀矿平米铀量等值线图成图方法,其特征在于:所述平米铀量等值线图图幅宽度=图框的最大X-图框的最小X。
4.如权利要求3所述的一种砂岩型铀矿平米铀量等值线图成图方法,其特征在于:
当m≥1时,zknUm表示n号钻孔平米铀量等值线值,zknUm≥m;
当m<1时,zknUm表示n号钻孔平米铀量等值线值,zknUm<1;
Um表示平米铀量等值线值,zkn表示钻孔号,m表示平米铀量值值数目。
5.如权利要求1所述的一种砂岩型铀矿平米铀量等值线图成图方法,其特征在于:所述第五步,对不同数值平米铀量等值线区域进行颜色填充形成平米铀量等值线图;
对不同数值平米铀量等值线值形成闭合区进行不同颜色填充。
CN201910857956.6A 2019-09-11 2019-09-11 一种砂岩型铀矿平米铀量等值线图成图方法 Active CN112489148B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910857956.6A CN112489148B (zh) 2019-09-11 2019-09-11 一种砂岩型铀矿平米铀量等值线图成图方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910857956.6A CN112489148B (zh) 2019-09-11 2019-09-11 一种砂岩型铀矿平米铀量等值线图成图方法

Publications (2)

Publication Number Publication Date
CN112489148A CN112489148A (zh) 2021-03-12
CN112489148B true CN112489148B (zh) 2024-05-14

Family

ID=74920138

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910857956.6A Active CN112489148B (zh) 2019-09-11 2019-09-11 一种砂岩型铀矿平米铀量等值线图成图方法

Country Status (1)

Country Link
CN (1) CN112489148B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112951069A (zh) * 2019-11-26 2021-06-11 核工业二一六大队 一种砂岩型铀矿矿砂比等值线图成图方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4059760A (en) * 1974-11-18 1977-11-22 Murphy John A Georadiological surveying for oil and gas and subsurface structure contour mapping
CN103837908A (zh) * 2014-03-05 2014-06-04 核工业北京地质研究院 一种适用于隐伏砂岩型铀矿快速找矿定位方法
CN104504047A (zh) * 2014-12-19 2015-04-08 张夏林 一种基于双层数据库的固体矿产资源储量估算系统
CN107505659A (zh) * 2017-07-05 2017-12-22 佛山诸广矿业有限公司 一种寻找花岗岩地区隐伏铀矿的方法
CN109270567A (zh) * 2018-09-11 2019-01-25 湖南中核环保科技有限公司 一种在铀矿勘查中制成探槽γ照射量率等值图的方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
PE20141748A1 (es) * 2011-05-09 2014-12-07 Hydrocarbon Imaging Services Inc Sistema y metodo para la deteccion de hidrocarburos
CN105738970B (zh) * 2016-02-29 2017-04-05 山东科技大学 一种共生共存四位一体矿产协同勘查方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4059760A (en) * 1974-11-18 1977-11-22 Murphy John A Georadiological surveying for oil and gas and subsurface structure contour mapping
CN103837908A (zh) * 2014-03-05 2014-06-04 核工业北京地质研究院 一种适用于隐伏砂岩型铀矿快速找矿定位方法
CN104504047A (zh) * 2014-12-19 2015-04-08 张夏林 一种基于双层数据库的固体矿产资源储量估算系统
CN107505659A (zh) * 2017-07-05 2017-12-22 佛山诸广矿业有限公司 一种寻找花岗岩地区隐伏铀矿的方法
CN109270567A (zh) * 2018-09-11 2019-01-25 湖南中核环保科技有限公司 一种在铀矿勘查中制成探槽γ照射量率等值图的方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
砂岩型铀矿找矿勘探中钻孔资料的管理与应用技术;蔡煜琦;黄树桃;赵永安;李胜祥;韩效忠;郑恩玖;韩绍阳;单明;;世界核地质科学(04);全文 *

Also Published As

Publication number Publication date
CN112489148A (zh) 2021-03-12

Similar Documents

Publication Publication Date Title
Zhang et al. Collaborative 3D geological modeling analysis based on multi-source data standard
CN101114276A (zh) 固体矿床三维可视化储量计算系统及计算方法
CN103902777A (zh) 一种开采沉陷变形预测预报方法
CN112489148B (zh) 一种砂岩型铀矿平米铀量等值线图成图方法
CN105205864A (zh) 基于多源数据的地质结构面三维模型自动建模方法和系统
CN112150582B (zh) 一种面向多模态数据的地质剖面图近似表达方法
CN104459800A (zh) 一种预测砂体尖灭的方法和装置
CN112907693B (zh) 一种砂岩型铀矿砂体顶低板等高线图成图方法
CN111239822B (zh) 一种消除井震闭合差的断层下速度建模方法及处理终端
CN112489163B (zh) 一种砂岩型铀矿综合成果图成图方法
CN112908145A (zh) 一种砂岩型铀矿砂体等厚度图成图方法
CN114626230B (zh) 一种地质矿产资源储量的模拟方法
CN116699715A (zh) 一种高光谱遥感信息参与的铀成矿有利区预测方法
CN115046515A (zh) 基于单滑面d型深孔测斜曲线的滑动面位置精准确定方法
CN105741329B (zh) 一种基于孔壁图像的相邻钻孔结构面连通性分析方法
CN104778296B (zh) 钻孔柱状图岩性描述布局数学规划方法
CN112415628B (zh) 一种砂岩型铀矿取样平面分布图成图方法
Amadu et al. Geostatistical Approach for The Estimation of Shear-Hosted Gold Deposit: A Case Study of The Obuasi Gold Deposit, Ghana
CN110990936A (zh) 一种岩体质量指标的确定方法及系统
Korus et al. Integrating AEM and borehole data for regional hydrogeologic synthesis: Tools and examples from Nebraska, USA
CN114187346A (zh) 一种砂岩型铀矿地质剖面图矿体倾角计算方法
CN112211621A (zh) 一种砂岩型煤矿储量估算剖面图成图方法
CN114742937B (zh) 一种三维地质解析方法和装置
CN117251739B (zh) 用于砂岩型铀资源定量预测的样本集的构建方法
CN115809815B (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