CN111896990B - 一种基于Frechet距离的放射源活度监测方法 - Google Patents
一种基于Frechet距离的放射源活度监测方法 Download PDFInfo
- Publication number
- CN111896990B CN111896990B CN202010662140.0A CN202010662140A CN111896990B CN 111896990 B CN111896990 B CN 111896990B CN 202010662140 A CN202010662140 A CN 202010662140A CN 111896990 B CN111896990 B CN 111896990B
- Authority
- CN
- China
- Prior art keywords
- activity
- radioactive source
- decay
- monitoring
- theoretical
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
- G01T1/16—Measuring radiation intensity
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T7/00—Details of radiation-measuring instruments
- G01T7/12—Provision for actuation of an alarm
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Health & Medical Sciences (AREA)
- Mathematical Analysis (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Molecular Biology (AREA)
- High Energy & Nuclear Physics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Theoretical Computer Science (AREA)
- Algebra (AREA)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Operations Research (AREA)
- Measurement Of Radiation (AREA)
Abstract
本发明公开一种基于Frechet距离的放射源活度监测方法,通过对放射源的实测活度值进行采集,建立源数据,自适应产生理论活度衰变曲线,从而通过比较分析理论活度衰变曲线与实测转换后的活度衰变曲线的相似度,来实现对放射源活度监测,把控性更强,解决针对某种特定放射源检测并提高准确度的问题,从而投入实际的监测工作中,监测放射源是否被移动或者调换,十分的便捷。同时对监测系统的数据进行不断的调试,提高监测的精确度。并且需要监测不同的放射源时,只用自定义输入不同放射源的相关参数即可,不需要再重新构建监测网络,因此能具有灵活性和可移植性。
Description
技术领域
本发明属于放射性安全于检测、预警技术领域。确切的讲,是涉及一种基于Frechet距离的放射源活度监测方法。
背景技术
伴随着核技术在各个行业扮演着越来越重要的角色,放射性核素的应用也愈加广泛。而放射源是指由放射性核素制作而成的辐射源,主要应用于工业、医疗、军事、农业等方面。放射源的利用,对我国乃至全世界的经济、社会发展具有重要意义。但是,放射源的不当使用或丢失,会导致一系列的问题,严重的甚至会威胁生命。放射源放出的射线会产生生物学效应,进而能够破坏人体器官组织,诱发各种疾病,甚至造成永久性损伤。由于放射源释放的射线无法用人眼直接观测,因此,放射源丢失或被盗事故屡见不鲜。
目前,在实际工作中,具有放射源使用许可的单位或多或少存在一定的监测管理问题,以至于放射源的污染问题层出不穷,如放射源的泄露、丢失、不当转移、偷窃等,严重影响我国社会的发展,威胁人民健康安全。因此,快速、有效的放射源活度监测方法,放射源安全监测和管理的不可或缺的组成部分。针对以上基础,本文提出一种利用放射源活度曲线图进来监测放射源,目的是防止放射源的掉包、丢失。
发明内容
针对现有技术的不足,本发明提供了基于Frechet距离的放射源活度监测方法,本发明的目的在于打破原有思想的束缚,提供以基于Frechet距离的放射源活度监测方法代替传统监测方法,通过计算监测放射性活度曲线与自适应活度衰变曲线的相似度,以及各种调试,最终做到减小误差,以提高放射源活度监测的准确度。
为实现以上目的,本发明通过以下技术方案予以实现:基于Frechet距离的放射源活度监测方法方法,包括以下过程步骤:
步骤01、在放射源活度监测系统自定义设置放射源的衰变常数λ、半衰期T1/2及初始放射性活度A0。
步骤02、在放射源活度监测系统自定义设置放射源的活度衰变系数ε与阈值。
步骤03、放射源活度监测系统开始记录放射源的实测活度值,并根据放射源的理论活度衰变公式,计算出放射源的理论活度值。
步骤04、放射源活度监测系统根据实测活度值及理论活度值,分别绘制放射源的实测活度衰变曲线及理论活度衰变曲线。
步骤05、建立线性回归方程,根据Frechet距离表达式算出实测活度衰变曲线及理论活度衰变曲线的空间曲线相似度。
步骤06、算出实测活度衰变曲线及理论活度衰变曲线的Frechet距离,并与步骤02中的阈值相互比较。
步骤07、返回比较结果,大于阈值则报警,小于阈值则继续监测。
步骤08、将放射源的实际位移情况与放射源活度监测系统监测结果相比较,得出准确度。
步骤09、准确度较高,监测系统具有可行性,则建立统一的放射性活度衰变曲线和线性回归方程;准确度过低,则调整放射源活度衰变系数ε及阈值,并重复步骤01至步骤08。
步骤10:利用优化后监测系统监测放射源的实测活度衰变曲线及理论活度衰变曲线的Frechet距离,并进行报警。
进一步地,基于Frechet距离的放射源活度监测方法还包括对放射源种类进行分类,根据强度不同分为低放射性、中放射性及高放射性,便于调整放射源活度衰变系数ε。
进一步地,放射源的实测活度值为衰变时间t、初始放射性活度A0以及衰变常数λ共同作用得出,放射源的理论活度值为自定义的初始放射性活度A0、衰变常数λ、活度衰变系数ε以及监测系统工作时间共同作用得出。
与现有技术相比,本发明具有以下有益效果:
(1)本发明基于Frechet距离的放射源活度监测方法,通过对放射源的实测活度值进行采集,建立源数据,产生理论活度衰变曲线,从而通过比较分析理论活度衰变曲线与实测转换后的活度衰变曲线的相似度,来实现对放射源活度监测,把控性更强,解决针对某种特定放射源检测并提高准确度的问题,从而投入实际的监测工作中,监测放射源是否被移动或者调换,十分的便捷。
(2)本发明通过对放射源活度监测系统的数据进行不断的完善,提高监测的精确度,同时,需要监测不同的放射源时,只用自定义输入不同放射源的相关参数即可,不需要再重新构建监测网络,因此能具有灵活性和可移植性。
附图说明
图1是本发明所提供的基于Frechet距离的放射源活度监测流程图;
图2是放射源活度监测正常情况下曲线图;
图3是放射源活度监测正常情况下结果图;
图4是放射源活度监测异常情况下曲线图;
图5是放射源活度监测异常情况下结果图。
具体实施方式
下面将结合本发明的实施例,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图1-5及具体实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
如图1至图5所示,一种基于Frechet距离的放射源活度监测方法,包括以下步骤:
步骤01、根据放射源种类,查阅资料,找到放射源的衰变常数λ、半衰期T1/2、初始放射性活度A0作为源数据。
步骤02、设置活度衰变系数,用于记录自适应活度衰变值,设置阈值,用于后续作为理论活度衰变曲线与理实测活度衰变曲线的Frechet距离是否达到报警的临界值,若大于阈值则报警,小于阈值则为正常现象继续监测。
步骤03、放射源活度监测系统开始记录放射源的实测活度值,根据活度衰变公式(1-1)和半衰期T1/2与衰变常数λ直接的关系,公式(1-2):
A=A0e-λT (1-1)
改写活度的另一种关系式(1-3):
其中,如待测放射源确定后,-ln2/T1/2就为一常数,可用D来表示,因此,某种确定后的放射源活度衰变公式可以用公式(1-4)来表示:
步骤04、放射源活度监测系统放射源活度监测系统,计算出放射源的理论活度值。放射源的理论活度衰变满足公式(1-5):
其中,ε为理论放射源的活度衰变系数,它的值是根据探测总效率及其他误差影响所定。
步骤05、放射源活度监测系统根据实测活度值及理论活度值,分别绘制放射源的实测活度衰变曲线及理论活度衰变曲线,将放射源的实测活度衰变曲线与理论活度衰变曲线进行分别绘制,使曲线变化趋势更加清晰明显。
步骤06、建立回归方程,根据更加严格精确度的Frechet距离表达式(1-6)计算出放射源的实测活度衰变曲线及理论活度衰变曲线的空间曲线相似度。
其中,α、β是在0到1区间上的任意连续不递减函数,P,Q为给定的两条曲线,δF为P、Q求解的Frechet距离。
步骤07、算出实测活度衰变曲线及理论活度衰变曲线的Frechet距离,将结果与设定的阈值相互比较,如果大于设定的阈值则报警,小于阈值则继续监测。
步骤08、将放射源的实际位移情况与放射源活度监测系统监测结果相比较,得出准确度。例如,放射源活度监测系统监测的实测活度衰变曲线及理论活度衰变曲线的Frechet距离大于阈值,证明放射源被移动或者调换,但实验测试人员并未移动或者调换放射源,则证明放射源活度监测系统的准确度不高,反之,则可以认为放射源活动监测系统具有较高的准确率,具有可行性,可以运用到实际的放射源监控工作。
步骤09、若放射源活度监测系统具有较高的准确率,具有可行性,则建立放射性活度衰变曲线和线性回归方程的统一,即经过不断的调整活度衰变系数ε和阈值后,系统的准确率不断提高满足所需,此时的放射源的理论活度衰变曲线即作为同类放射源监测的标准曲线,放射源的实测活度衰变曲线通过线性回归方程直接算出与标准曲线的Frechet距离,便于快速准备的监测到放射源的异动情况。
如放射源活度监测系统准确度过低,则调整放射源活度衰减系数ε及阈值,并重复步骤01至步骤08。如此经过多次实验,并对参数进行调整,例如公式中的放射源活度衰变系数ε以及阈值,既要提高放射源活动监控系统的精确度也避免整个过拟合化。
步骤10、利用优化后的放射源活度监测系统监测放射源的实测活度衰变曲线及理论活度衰变曲线的Frechet距离,并进行报警。
影响因素具体是指放射源的衰变常数λ、半衰期T1/2、初始放射性活度A0以及与之相关的所有此方法可验证为可利用的数据。
基于Frechet距离的放射源活度监测方法代替传统监测方法,通过计算监测放射性的实测活度衰变曲线与理论活度衰变曲线的相似度,算出实测活度衰变曲线及理论活度衰变曲线的Frechet距离,通过返回值与设置的阈值比较,大于阈值则报警,小于阈值则继续监测,能够快速的监测出放射源的各种异动情况。同时通过监测结果与放射源的实际位移情况做比较,不断的调试参数,使得放射源活度监测系统的精确率满足监控所需,以提高放射源活度监测的准确度。并且,对放射源种类进行分类,根据强度不同分为低放射性、中放射性及高放射性,便于调整放射源活度衰变系数ε。
利用理论放射源的活度衰变曲线进行放射源的异动预测,根据实际需要监测的放射源的情况和要求重新定义放射源活度监测系统的参数即可,不需要再重新构建网络,因此能具有灵活性和可移植性。理论活度衰变曲线与实测活度衰变曲线的距离是不同的,Frechet距离就是寻找一对最小化理论活度值和实测活度值之间的最长距离的函数来进行分析比较,提高了监测系统的可靠性。
以上显示和描述了本发明的基本原理和主要特征和本发明的优点。本行业的技术人员应该了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的只是说明本发明的原理,在不脱离本发明精神和范围的前提下,本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。
Claims (3)
1.一种基于Frechet距离的放射源活度监测方法,其特征在于,包括以下步骤:
步骤01、在放射源活度监测系统自定义设置放射源的衰变常数λ、半衰期T1/2及初始放射性活度A0;
步骤02、在放射源活度监测系统自定义设置放射源的活度衰变系数ε与阈值;
步骤03、放射源活度监测系统开始记录放射源的实测活度值,并根据放射源的理论活度衰变公式,计算出放射源的理论活度值;
所述放射源的理论活度衰变公式为:
步骤04、放射源活度监测系统根据实测活度值及理论活度值,分别绘制放射源的实测活度衰变曲线及理论活度衰变曲线;
步骤05、建立线性回归方程,根据Frechet距离表达式算出实测活度衰变曲线及理论活度衰变曲线的空间曲线相似度;
所述Frechet距离表达式为:
其中,α、β是在0到1区间上的任意连续不递减函数,P,Q为给定的两条曲线,δF为P、Q求解的Frechet距离;
步骤06、算出实测活度衰变曲线及理论活度衰变曲线的Frechet距离,并与步骤02中的阈值相互比较;
步骤07、返回比较结果,大于阈值则报警,小于阈值则继续监测;
步骤08、将放射源的实际位移情况与监测系统监测结果相比较,得出准确度;
步骤09、准确度较高,监测系统具有可行性,则建立统一的放射性活度衰变曲线和线性回归方程;准确度过低,则调整放射源的活度衰减系数ε及阈值,并重复步骤01至步骤08;
步骤10、利用优化后的放射源活度监测系统监测放射源的实测活度衰变曲线及理论活度衰变曲线的Frechet距离,并进行报警。
2.根据权利要求1所述的方法,其特征在于:还包括对放射源种类进行分类,根据强度不同分为低放射性、中放射性及高放射性,调整放射源活度衰变系数ε。
3.根据权利要求1所述的方法,其特征在于:所述放射源的实测活度值为衰变时间t、初始放射性活度A0以及衰变常数λ共同作用得出;放射源的理论活度值为自定义的初始放射性活度A0、衰变常数λ、活度衰变系数ε以及监测系统工作时间共同作用得出。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010662140.0A CN111896990B (zh) | 2020-07-10 | 2020-07-10 | 一种基于Frechet距离的放射源活度监测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010662140.0A CN111896990B (zh) | 2020-07-10 | 2020-07-10 | 一种基于Frechet距离的放射源活度监测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111896990A CN111896990A (zh) | 2020-11-06 |
CN111896990B true CN111896990B (zh) | 2023-01-24 |
Family
ID=73192235
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010662140.0A Active CN111896990B (zh) | 2020-07-10 | 2020-07-10 | 一种基于Frechet距离的放射源活度监测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111896990B (zh) |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2000013042A1 (en) * | 1998-09-01 | 2000-03-09 | Lockheed Martin Idaho Technologies Company | Radioactive decay detection method and apparatus |
JP2005249747A (ja) * | 2004-03-08 | 2005-09-15 | Japan Electric Cable Technology Center Inc | 大面積放射線線量分布測定法 |
EP2887103A1 (en) * | 2013-12-23 | 2015-06-24 | Services Pétroliers Schlumberger | Artificially activated radioactive source for a downhole tool |
CN107064989A (zh) * | 2017-03-09 | 2017-08-18 | 上海新漫传感技术研究发展有限公司 | 液体闪烁谱仪、获取淬灭校正曲线和放射性活度的方法 |
CN107643537A (zh) * | 2016-07-21 | 2018-01-30 | 环境保护部核与辐射安全中心 | 航空辐射监测方法 |
CN108873047A (zh) * | 2018-06-30 | 2018-11-23 | 上海联影医疗科技有限公司 | 检测放射源活度的方法、系统、计算机设备和存储介质 |
CN109633729A (zh) * | 2018-12-29 | 2019-04-16 | 哈尔滨工业大学 | 一种移动机器人核环境下的放射源定位及强度估计方法 |
CN110082810A (zh) * | 2019-05-09 | 2019-08-02 | 成都理工大学 | 一种拟合闪烁室k值的刻度装置 |
CN110146915A (zh) * | 2019-06-12 | 2019-08-20 | 成都理工大学 | 一种低活度γ能谱多峰稳谱方法 |
RU2722203C1 (ru) * | 2019-07-15 | 2020-05-28 | Российская Федерация, от имени которой выступает Государственная корпорация по атомной энергии "Росатом" Госкорпорация "Росатом" | Способ контроля содержания радионуклидов в радиоактивных отходах |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
MXPA04005655A (es) * | 2001-12-12 | 2005-04-19 | Univ Alberta The University Of | Ion radioactivo. |
EP2077453A1 (de) * | 2008-01-01 | 2009-07-08 | SMA Solar Technology AG | Auswerteverfahren |
DE102010033609B3 (de) * | 2010-08-06 | 2011-12-22 | Siemens Aktiengesellschaft | Verfahren zur Schätzung einer Strahlendosis einer Röntgenstrahlung und zugehörige Röntgenvorrichtung |
US20140224964A1 (en) * | 2013-02-08 | 2014-08-14 | Image Insight Inc. | Apparatus and method for detection of radiation |
FR3085492B1 (fr) * | 2018-08-28 | 2021-06-04 | Icohup | Procede de spectrometrie et dispositif de detection de rayonnements ionisants pour sa mise en œuvre |
-
2020
- 2020-07-10 CN CN202010662140.0A patent/CN111896990B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2000013042A1 (en) * | 1998-09-01 | 2000-03-09 | Lockheed Martin Idaho Technologies Company | Radioactive decay detection method and apparatus |
JP2005249747A (ja) * | 2004-03-08 | 2005-09-15 | Japan Electric Cable Technology Center Inc | 大面積放射線線量分布測定法 |
EP2887103A1 (en) * | 2013-12-23 | 2015-06-24 | Services Pétroliers Schlumberger | Artificially activated radioactive source for a downhole tool |
CN107643537A (zh) * | 2016-07-21 | 2018-01-30 | 环境保护部核与辐射安全中心 | 航空辐射监测方法 |
CN107064989A (zh) * | 2017-03-09 | 2017-08-18 | 上海新漫传感技术研究发展有限公司 | 液体闪烁谱仪、获取淬灭校正曲线和放射性活度的方法 |
CN108873047A (zh) * | 2018-06-30 | 2018-11-23 | 上海联影医疗科技有限公司 | 检测放射源活度的方法、系统、计算机设备和存储介质 |
CN109633729A (zh) * | 2018-12-29 | 2019-04-16 | 哈尔滨工业大学 | 一种移动机器人核环境下的放射源定位及强度估计方法 |
CN110082810A (zh) * | 2019-05-09 | 2019-08-02 | 成都理工大学 | 一种拟合闪烁室k值的刻度装置 |
CN110146915A (zh) * | 2019-06-12 | 2019-08-20 | 成都理工大学 | 一种低活度γ能谱多峰稳谱方法 |
RU2722203C1 (ru) * | 2019-07-15 | 2020-05-28 | Российская Федерация, от имени которой выступает Государственная корпорация по атомной энергии "Росатом" Госкорпорация "Росатом" | Способ контроля содержания радионуклидов в радиоактивных отходах |
Non-Patent Citations (6)
Title |
---|
Developments in radiation, stimulation and observation facilities in luminescence measurements;L.Bøtter-Jensen 等;《Radiation Measurements》;20030408;第37卷(第4-5期);535-541 * |
γ射线测量过程控制及其测量数据处理研究;邱化冬;《中国优秀硕士学位论文全文数据库 (工程科技Ⅱ辑)》;20070615(第6期);C040-8 * |
一种折中的小波阈值方法对γ能谱平滑处理的研究;颜瑜成 等;《原子核物理评论》;20140920;第31卷(第3期);407-410 * |
激光诱导的光核中子和等离子体对撞聚变中子产生研究;齐伟;《中国博士学位论文全文数据库 (工程科技Ⅱ辑)》;20200115(第1期);C040-31 * |
煤样吸附/解吸滞后效应定量分析;党广兴 等;《煤炭科学技术》;20170515;第45卷(第5期);187-191、222 * |
重庆市国省干线公路交通特征与路面使用性能衰变模型;赵桂香;《中国优秀硕士学位论文全文数据库 (工程科技Ⅱ辑)》;20160615(第6期);C034-34 * |
Also Published As
Publication number | Publication date |
---|---|
CN111896990A (zh) | 2020-11-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7698690B2 (en) | Identifying code that wastes time performing redundant computation | |
KR101194746B1 (ko) | 침입코드 인식을 위한 코드 모니터링 방법 및 장치 | |
CN109871764A (zh) | 一种异常行为识别方法、装置及存储介质 | |
Park et al. | Anomaly intrusion detection by clustering transactional audit streams in a host computer | |
TW201935237A (zh) | 異常檢測的方法及裝置 | |
CN113992340B (zh) | 用户异常行为识别方法、装置、设备和存储介质 | |
CN111896990B (zh) | 一种基于Frechet距离的放射源活度监测方法 | |
CN116112292A (zh) | 基于网络流量大数据的异常行为检测方法、系统和介质 | |
CN109308225A (zh) | 一种虚拟机异常检测方法、装置、设备及存储介质 | |
US8898648B2 (en) | Methodology for fast detection of false sharing in threaded scientific codes | |
CN113141276A (zh) | 一种基于知识图谱的信息安全方法 | |
US9465678B2 (en) | Event correlation | |
You et al. | ZeroSpy: exploring software inefficiency with redundant zeros | |
CN105825130A (zh) | 一种信息安全预警方法及装置 | |
Luktarhan et al. | Multi-stage attack detection algorithm based on hidden markov model | |
CN117955863A (zh) | 基于人工智能的数据安全检测方法及系统 | |
Ince et al. | An empirical and theoretical analysis of an information flow-based system design metric | |
WO2023098410A1 (zh) | 一种核电厂纵深防御层次独立性量化评价方法及实施该方法的系统 | |
Talebi‐Ghane et al. | Joint frailty model for recurrent events and death in presence of cure fraction: Application to breast cancer data | |
CN115567241A (zh) | 一种多站点网络感知检测系统 | |
WO2022239030A1 (en) | Method and system for anomaly detection in the banking system with graph neural networks (gnns) | |
CN117521762A (zh) | 一种基于通道李普希兹审查的后门移除方法 | |
CN117540372B (zh) | 智能学习的数据库入侵检测与响应系统 | |
Li et al. | Denoising Network of Dynamic Features for Enhanced Malware Classification | |
CN102929774B (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 |