CN107462929A - 一种井中铜镍矿产品位测量装置及方法 - Google Patents
一种井中铜镍矿产品位测量装置及方法 Download PDFInfo
- Publication number
- CN107462929A CN107462929A CN201710615018.6A CN201710615018A CN107462929A CN 107462929 A CN107462929 A CN 107462929A CN 201710615018 A CN201710615018 A CN 201710615018A CN 107462929 A CN107462929 A CN 107462929A
- Authority
- CN
- China
- Prior art keywords
- cupro
- well
- nickel
- gamma
- mineral
- 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
Links
- 229910052500 inorganic mineral Inorganic materials 0.000 title claims abstract description 61
- 239000011707 mineral Substances 0.000 title claims abstract description 61
- 229910000570 Cupronickel Inorganic materials 0.000 title claims abstract description 35
- 238000000034 method Methods 0.000 title claims abstract description 22
- 229910052751 metal Inorganic materials 0.000 claims abstract description 50
- 239000002184 metal Substances 0.000 claims abstract description 41
- 238000005259 measurement Methods 0.000 claims abstract description 22
- 238000001228 spectrum Methods 0.000 claims abstract description 15
- 239000002245 particle Substances 0.000 claims abstract description 14
- 238000012545 processing Methods 0.000 claims abstract description 10
- 210000002445 nipple Anatomy 0.000 claims abstract description 8
- WFKWXMTUELFFGS-UHFFFAOYSA-N tungsten Chemical compound [W] WFKWXMTUELFFGS-UHFFFAOYSA-N 0.000 claims abstract description 5
- 229910052721 tungsten Inorganic materials 0.000 claims abstract description 5
- 239000010937 tungsten Substances 0.000 claims abstract description 5
- 229910001200 Ferrotitanium Inorganic materials 0.000 claims abstract description 4
- 238000001514 detection method Methods 0.000 claims description 17
- RYGMFSIKBFXOCR-UHFFFAOYSA-N Copper Chemical compound [Cu] RYGMFSIKBFXOCR-UHFFFAOYSA-N 0.000 claims description 13
- 229910052802 copper Inorganic materials 0.000 claims description 13
- 239000010949 copper Substances 0.000 claims description 13
- 238000012360 testing method Methods 0.000 claims description 3
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 claims 1
- 239000010931 gold Substances 0.000 claims 1
- 229910052737 gold Inorganic materials 0.000 claims 1
- 230000009286 beneficial effect Effects 0.000 abstract description 2
- 230000004907 flux Effects 0.000 abstract description 2
- 238000011002 quantification Methods 0.000 abstract description 2
- PXHVJJICTQNCMI-UHFFFAOYSA-N Nickel Chemical compound [Ni] PXHVJJICTQNCMI-UHFFFAOYSA-N 0.000 description 16
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 10
- 229910052759 nickel Inorganic materials 0.000 description 8
- 229910052742 iron Inorganic materials 0.000 description 5
- 239000011435 rock Substances 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 230000005251 gamma ray Effects 0.000 description 3
- 239000011159 matrix material Substances 0.000 description 3
- 238000013461 design Methods 0.000 description 2
- 241001637516 Polygonia c-album Species 0.000 description 1
- 238000007792 addition Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 229910052951 chalcopyrite Inorganic materials 0.000 description 1
- DVRDHUBQLOKMHZ-UHFFFAOYSA-N chalcopyrite Chemical compound [S-2].[S-2].[Fe+2].[Cu+2] DVRDHUBQLOKMHZ-UHFFFAOYSA-N 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- YOCUPQPZWBBYIX-UHFFFAOYSA-N copper nickel Chemical compound [Ni].[Cu] YOCUPQPZWBBYIX-UHFFFAOYSA-N 0.000 description 1
- 239000013078 crystal Substances 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000004992 fission Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 229910052954 pentlandite Inorganic materials 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V5/00—Prospecting or detecting by the use of ionising radiation, e.g. of natural or induced radioactivity
- G01V5/04—Prospecting or detecting by the use of ionising radiation, e.g. of natural or induced radioactivity specially adapted for well-logging
- G01V5/08—Prospecting or detecting by the use of ionising radiation, e.g. of natural or induced radioactivity specially adapted for well-logging using primary nuclear radiation sources or X-rays
- G01V5/10—Prospecting or detecting by the use of ionising radiation, e.g. of natural or induced radioactivity specially adapted for well-logging using primary nuclear radiation sources or X-rays using neutron sources
- G01V5/104—Prospecting or detecting by the use of ionising radiation, e.g. of natural or induced radioactivity specially adapted for well-logging using primary nuclear radiation sources or X-rays using neutron sources and detecting secondary Y-rays as well as reflected or back-scattered neutrons
Landscapes
- Physics & Mathematics (AREA)
- High Energy & Nuclear Physics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种井中铜镍矿产品位测量装置及方法。该装置包括井下设备和地面数据处理设备,井下设备和地面数据处理设备通过电缆相连,井下设备包括壳体,壳体为钛钢筒状壳体,壳体内自下而上依次设置有D‑T可控中子源、钨屏蔽体、CLYC双粒子探测器和遥测短节。该方法主要包括:首先根据测量得到的俘获伽马能谱,选取高能段伽马计数信息,计算不同金属元素的伽马产额,然后对伽马产额与热中子计数做比值,消除热中子通量不稳定对伽马产额的影响,最后根据伽马产额‑热中子计数比值与金属矿产品位的线性关系,计算金属矿产品位。该方法可在井中提供连续深度的金属矿产品位,有益于金属矿产勘探的定量化,提升了金属矿产勘探效率。
Description
技术领域
本发明涉及金属矿产勘探领域,尤其涉及一种井中铜镍矿产品位测量装置及方法。
背景技术
目前对金属矿产勘探的手段主要包括航空物探、地面物探及综合测井方法,利用这些方法可对矿体可能埋藏的位置、规模和产状做出判断,了解矿体在地下的分布和延伸情况。针对金属矿产品位确定方法,常采用对矿体和围岩进行系统的取样和分析,了解矿体和围岩的物质成分。
瞬发伽马中子活化分析技术(PGNAA)常应用于地面的在线金属矿产品位计算,利用PGNAA技术在井中测量过程中,中子源发射快中子进入井眼和地层,与矿产金属元素发生辐射俘获反应,产生俘获特征伽马射线,被探测器记录,通过在俘获伽马能谱上截取特征能窗的方法获取金属元素特征伽马计数,利用伽马能窗计数与金属品位的拟合关系,最终计算得到金属矿产品位。常见铜镍金属矿产中,铜镍铁元素常相互伴生聚集成矿,以铜镍硫化物矿床为例,其矿石中的金属矿物以磁黄铁矿、镍黄铁矿和黄铜矿为主。因此,利用能窗法获取不同元素特征俘获伽马计数的过程中,易受到相互之间特征伽马射线的影响,特别是,铁、铜、镍三种金属元素产生的俘获伽马射线均集中与7.0-9.5MeV处,会对伽马产额计算结果产生较大影响。
发明内容
本发明的目的是针对上述不足,提出了一种通过在井中采集不同深度的伽马能谱和热中子计数信息,根据特征伽马产额与热中子计数的比值确定金属矿产品位井中铜镍矿产品位测量装置及方法。
本发明具体采用如下技术方案:
一种井中铜镍矿产品位测量装置,包括井下设备和地面数据处理设备,所述井下设备包括壳体,壳体内自下而上依次设置有D-T可控中子源、钨屏蔽体、CLYC双粒子探测器和遥测短节。
优选地,所述壳体为钛钢筒状壳体。
一种井中铜镍矿产品位测量方法,采用如上所述的一种井中铜镍矿产品位测量装置,包括以下步骤:
步骤一:D-T可控中子源以脉冲形式向地层发射快中子,CLYC双粒子探测器同时探测俘获伽马能谱和热中子计数,遥测短节将测量得到的数据上传至地面数据处理设备;
步骤二:截取7.0-9.5MeV处的俘获伽马能谱,利用最小二乘方法,计算得到铜和镍元素的伽马产额;
步骤三:根据伽马产额/热中子计数与铜镍金属矿产品位的线性关系,计算得到铜镍金属矿产品位。
优选地,所述步骤一中的D-T可控中子源脉冲周期为400μs,一个脉冲周期内,0-40μs发射快中子,40-400μs快中子停止发射。
优选地,所述步骤一中的CLYC双粒子探测器探测周期和D-T可控中子源脉冲周期一致,探测周期为400μs,一个探测周期内,0-50μs停止探测,50-400μs探测热中子和俘获伽马能谱。
优选地,所述步骤三中的铜镍金属矿产品位计算时,伽马产额-热中子计数比值与铜镍金属矿产品位的线性关系在刻度井中获取。
本发明具有如下有益效果:本发明利用可控D-T可控中子源,一个CLYC双粒子探测器,同时测量俘获伽马能谱和热中子计数,避免了化学中子源对操作人员的,减少了热中子通量对伽马产额的影响,克服了铜、镍、铁等元素产生的俘获伽马计数之间的相互影响,可提供连续深度的铜和镍金属矿产品位,有益于金属矿产勘探的定量化,提升了金属矿产勘探效率。
附图说明
图1为井中铜镍金属矿产品位测量结构原理示意图;
图2为铜、镍、铁三种金属元素在7.0-9.5MeV处的俘获伽马能谱;
图3为模拟井条件下,伽马产额-热中子计数比与铜镍金属矿产品位的关系图;
图4为铜镍金属矿产品位井中测量效果图。
其中,1为D-T可控中子源,2为钨屏蔽体,3为CLYC双粒子探测器,4为井眼,5为遥测短节,6为金属矿藏,7为矿藏围岩,8为壳体,9为地面数据处理设备。
具体实施方式
下面结合附图和具体实施例对本发明的具体实施方式做进一步说明:
如图1所示,一种井中铜镍矿产品位测量装置,包括井下设备和地面数据处理设备9,井下设备和地面数据处理设备9通过电缆相连,井下设备包括壳体9,壳体9为钛钢筒状壳体,壳体9内自下而上依次设置有D-T可控中子源1、钨屏蔽体2、CLYC双粒子探测器3和遥测短节5。
上述测量装置在井眼4中进行连续深度的测量,井眼4周围为金属矿藏6,或者金属矿藏围岩7。
上述D-T可控中子源1发射中子的平均能量约为14MeV。
上述CLYC双粒子探测器3对伽马光子探测的能量分辨率约为4.1%(相对Cs-137的0.662MeV伽马射线)。
上述CLYC双粒子探测器3利用晶体中的6Li同位素与热中子作用产生带电粒子,实现热中子的测量。
一种井中铜镍矿产品位测量方法,采用如上所述的一种井中铜镍矿产品位测量装置,包括以下步骤:
步骤一:D-T可控中子源1以脉冲形式向地层发射快中子,CLYC双粒子探测器3同时探测俘获伽马能谱和热中子计数,遥测短节5将测量得到的数据上传至地面数据处理设备9;
步骤二:截取7.0-9.5MeV处的俘获伽马能谱,利用最小二乘方法,计算得到铜和镍元素的伽马产额;
步骤三:根据伽马产额/热中子计数与铜镍金属矿产品位的线性关系,计算得到铜镍金属矿产品位。
步骤一中的D-T可控中子源1脉冲周期为400μs,一个脉冲周期内,0-40μs发射快中子,40-400μs快中子停止发射。
步骤一中的CLYC双粒子探测器3探测周期和D-T可控中子源1脉冲周期一致,探测周期为400μs,一个探测周期内,0-50μs停止探测,50-400μs探测热中子和俘获伽马能谱。
步骤三中的铜镍金属矿产品位计算时,伽马产额-热中子计数比值与铜镍金属矿产品位的线性关系在刻度井中获取。
如图2所示,铜、镍、铁三种金属元素在7.0-9.5MeV能量范围伽马计数分布的不同,建立元素标准谱矩阵A,矩阵的每列数据代表某种金属元素在不同能量道的伽马计数,矩阵每行数据代表不同金属元素。将测量得到的7.0-9.5MeV能量范围伽马计数设为向量X,按照式(1)计算各金属元素伽马产额Y。
Y=(ATA)-1ATX (1)
计算得到各元素的伽马产额之后,可以利用式(2)计算同一深度点处的各金属矿品位。
其中,Gi为第i种金属矿产品位,yi为第i种金属元素伽马产额,Nn为热中子计数,K1和K2为线性拟合系数。
如图3所示,利用上述金属矿产测量方法与装置,在模拟井条件下,测量得到的伽马产额-热中子计数比值与不同铜镍金属矿产品位的线性关系图,可以看出,伽马-热中子计数比与金属矿产品位具有较高的线性关系,根据该线性关系,可计算得到金属矿产品位。
如图4所示,为本发明的井中测量效果图,第一道为深度;第二道为地层岩性剖面,第三道为计算得到铜和镍金属元素伽马产额和热中子计数,第四道为计算得到的镍金属矿产品位(虚线)和理论镍金属矿产品位(灰度充填),第五道为计算得到的铜金属矿产品位(虚线)和理论铜金属矿产品位(灰度充填)。
当然,上述说明并非是对本发明的限制,本发明也并不仅限于上述举例,本技术领域的技术人员在本发明的实质范围内所做出的变化、改型、添加或替换,也应属于本发明的保护范围。
Claims (6)
1.一种井中铜镍矿产品位测量装置,其特征在于,包括井下设备和地面数据处理设备,所述井下设备包括壳体,壳体内自下而上依次设置有D-T可控中子源、钨屏蔽体、CLYC双粒子探测器和遥测短节。
2.如权利要求1所述的一种井中铜镍矿产品位测量装置,其特征在于,所述壳体为钛钢筒状壳体。
3.一种井中铜镍矿产品位测量方法,采用如权利要求1或2所述的一种井中铜镍矿产品位测量装置,其特征在于,包括以下步骤:
步骤一:D-T可控中子源以脉冲形式向地层发射快中子,CLYC双粒子探测器同时探测俘获伽马能谱和热中子计数,遥测短节将测量得到的数据上传至地面数据处理设备;
步骤二:截取7.0-9.5MeV处的俘获伽马能谱,利用最小二乘方法,计算得到铜和镍元素的伽马产额;
步骤三:根据伽马产额/热中子计数与铜镍金属矿产品位的线性关系,计算得到铜镍金属矿产品位。
4.如权利要求3所述的一种井中铜镍矿产品位测量方法,其特征在于,所述步骤一中的D-T可控中子源脉冲周期为400μs,一个脉冲周期内,0-40μs发射快中子,40-400μs快中子停止发射。
5.如权利要求4所述的一种井中铜镍矿产品位测量方法,其特征在于,所述步骤一中的CLYC双粒子探测器探测周期和D-T可控中子源脉冲周期一致,探测周期为400μs,一个探测周期内,0-50μs停止探测,50-400μs探测热中子和俘获伽马能谱。
6.如权利要求3所述的一种井中铜镍矿产品位测量方法,其特征在于,所述步骤三中的铜镍金属矿产品位计算时,伽马产额-热中子计数比值与铜镍金属矿产品位的线性关系在刻度井中获取。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710615018.6A CN107462929B (zh) | 2017-07-25 | 2017-07-25 | 一种井中铜镍矿产品位测量装置及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710615018.6A CN107462929B (zh) | 2017-07-25 | 2017-07-25 | 一种井中铜镍矿产品位测量装置及方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107462929A true CN107462929A (zh) | 2017-12-12 |
CN107462929B CN107462929B (zh) | 2019-04-26 |
Family
ID=60547135
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710615018.6A Active CN107462929B (zh) | 2017-07-25 | 2017-07-25 | 一种井中铜镍矿产品位测量装置及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107462929B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114264681A (zh) * | 2021-12-29 | 2022-04-01 | 清华大学 | 分析金矿石品位的方法和系统 |
CN115291288A (zh) * | 2022-09-29 | 2022-11-04 | 中石化经纬有限公司 | 基于铁中子标记的随钻脉冲中子孔隙度智能处理方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3947683A (en) * | 1973-06-05 | 1976-03-30 | Texaco Inc. | Combination of epithermal and inelastic neutron scattering methods to locate coal and oil shale zones |
CN201835828U (zh) * | 2010-07-23 | 2011-05-18 | 中国石油化工集团公司 | 双脉冲中子伽马谱测井仪 |
CN102518431A (zh) * | 2011-12-26 | 2012-06-27 | 中国石油大学(华东) | 基于可控中子源的随钻多参数测井方法 |
WO2013059394A1 (en) * | 2011-10-21 | 2013-04-25 | Schlumberger Canada Limited | Elpasolite scintillator-based neutron detector for oilfield applications |
CN103513287A (zh) * | 2012-06-19 | 2014-01-15 | 王新光 | 一种利用直流可控中子源计算地层密度的测井方法 |
CN103696765A (zh) * | 2013-11-06 | 2014-04-02 | 中国石油大学(华东) | 基于可控中子源的双LaBr3探测器元素能谱测井仪及测井方法 |
US20160370493A1 (en) * | 2013-04-02 | 2016-12-22 | Morpho Detection, Llc | Cross-correlated gamma ray and neutron detector |
-
2017
- 2017-07-25 CN CN201710615018.6A patent/CN107462929B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3947683A (en) * | 1973-06-05 | 1976-03-30 | Texaco Inc. | Combination of epithermal and inelastic neutron scattering methods to locate coal and oil shale zones |
CN201835828U (zh) * | 2010-07-23 | 2011-05-18 | 中国石油化工集团公司 | 双脉冲中子伽马谱测井仪 |
WO2013059394A1 (en) * | 2011-10-21 | 2013-04-25 | Schlumberger Canada Limited | Elpasolite scintillator-based neutron detector for oilfield applications |
CN103890615A (zh) * | 2011-10-21 | 2014-06-25 | 普拉德研究及开发股份有限公司 | 用于油田应用的基于钾冰晶石闪烁体的中子探测器 |
CN102518431A (zh) * | 2011-12-26 | 2012-06-27 | 中国石油大学(华东) | 基于可控中子源的随钻多参数测井方法 |
CN103513287A (zh) * | 2012-06-19 | 2014-01-15 | 王新光 | 一种利用直流可控中子源计算地层密度的测井方法 |
US20160370493A1 (en) * | 2013-04-02 | 2016-12-22 | Morpho Detection, Llc | Cross-correlated gamma ray and neutron detector |
CN103696765A (zh) * | 2013-11-06 | 2014-04-02 | 中国石油大学(华东) | 基于可控中子源的双LaBr3探测器元素能谱测井仪及测井方法 |
Non-Patent Citations (2)
Title |
---|
张明荣: ""Ce3+激活的钾冰晶石类非氟卤化物闪烁晶体的研发现状"", 《中国稀土学会2017学术年会摘要集》 * |
朱力 等: ""MC法模拟地层元素测井仪优化中子屏蔽体"", 《核电子学与探测技术》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114264681A (zh) * | 2021-12-29 | 2022-04-01 | 清华大学 | 分析金矿石品位的方法和系统 |
CN115291288A (zh) * | 2022-09-29 | 2022-11-04 | 中石化经纬有限公司 | 基于铁中子标记的随钻脉冲中子孔隙度智能处理方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107462929B (zh) | 2019-04-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105093343B (zh) | 一种地层元素俘获伽马能谱的测量装置 | |
US9835757B2 (en) | Pulsed-neutron tool methods and systems for monitoring casing corrosion | |
CN108457640B (zh) | 融合瞬发中子时间谱修正自然γ总量的铀矿测井定量方法 | |
CN110908010B (zh) | 一种行之有效的找800米以内浅砂岩型铀矿地球物理方法 | |
CN108825220B (zh) | 融合自然γ能谱与中子时间谱的铀矿测井仪及铀定量方程 | |
CN106250619B (zh) | 一种确定地层矿物含量的方法和装置 | |
CN104963677B (zh) | 一种利用支撑剂探测确定压裂裂缝高度的方法 | |
CN107229080B (zh) | 一种元素测井中子俘获伽马能谱的获取方法 | |
CN103678778B (zh) | 一种放射性物化探信息综合方法 | |
CN103470252B (zh) | 基于超热中子时间谱的瞬发中子测井及铀矿定量方法 | |
CN104747179A (zh) | 基于氘-氚加速器中子源的地层密度随钻测量仪 | |
CN109915127A (zh) | 一种基于d-d可控中子源的密度测量方法 | |
CN107462929B (zh) | 一种井中铜镍矿产品位测量装置及方法 | |
CN109239106A (zh) | 一种井中地层水矿化度测量装置及方法 | |
US3792253A (en) | Method and apparatus for detection of copper | |
CN107829728A (zh) | 多源距随钻中子孔隙度测量装置及其测量方法 | |
CN110469324A (zh) | 一种基于脉冲中子测井的计算地层密度方法 | |
CN103711479A (zh) | 基于超热中子与热中子比值的铀裂变瞬发中子测井技术 | |
CN206035475U (zh) | 一种碳氢比脉冲中子多参数测井仪 | |
CN107075939A (zh) | 自然伽马射线工具的增益稳定 | |
CN112630849A (zh) | 一种基于能谱测井特征谱峰的铀矿定量剥离系数求法 | |
Steingrimsson | Geothermal well logging: Geological wireline logs and fracture imaging | |
CN112649888B (zh) | 一种基于能谱测井特征谱峰的铀矿定量刻度系数求法 | |
Stockmann et al. | Exploring short-term soil landscape formation in the Hunter Valley, NSW, using gamma ray spectrometry | |
Keys et al. | Well logging with californium-252 |
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 |