CN109157224B - 一种增加参考光源校准的脉搏血氧监控系统及方法 - Google Patents

一种增加参考光源校准的脉搏血氧监控系统及方法 Download PDF

Info

Publication number
CN109157224B
CN109157224B CN201810587483.8A CN201810587483A CN109157224B CN 109157224 B CN109157224 B CN 109157224B CN 201810587483 A CN201810587483 A CN 201810587483A CN 109157224 B CN109157224 B CN 109157224B
Authority
CN
China
Prior art keywords
blood oxygen
light source
pulse
reference light
calibration
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
CN201810587483.8A
Other languages
English (en)
Other versions
CN109157224A (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.)
University of Electronic Science and Technology of China
Original Assignee
University of Electronic Science and Technology of 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 University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CN201810587483.8A priority Critical patent/CN109157224B/zh
Publication of CN109157224A publication Critical patent/CN109157224A/zh
Application granted granted Critical
Publication of CN109157224B publication Critical patent/CN109157224B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/145Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue
    • A61B5/1455Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue using optical sensors, e.g. spectral photometrical oximeters
    • A61B5/14551Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue using optical sensors, e.g. spectral photometrical oximeters for measuring blood gases
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/145Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue
    • A61B5/1455Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue using optical sensors, e.g. spectral photometrical oximeters
    • A61B5/14551Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue using optical sensors, e.g. spectral photometrical oximeters for measuring blood gases
    • A61B5/14552Details of sensors specially adapted therefor
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B2560/00Constructional details of operational features of apparatus; Accessories for medical measuring apparatus
    • A61B2560/02Operational features
    • A61B2560/0223Operational features of calibration, e.g. protocols for calibrating sensors
    • A61B2560/0228Operational features of calibration, e.g. protocols for calibrating sensors using calibration standards
    • A61B2560/0233Optical standards

Abstract

本发明公开了一种通过增加参考光源实现更精确校准的脉搏血氧监控系统及方法,属于脉搏血氧监控技术领域;可以解决现有计算过程中近似处理造成精度较低的问题。采用的方法:在现有血氧探头采用660nm和940nm光作为测量光源的基础上,加入波长为800nm的参考光源用于系统校准,在系统校准时信号获取部分采用改进后的三光源探头,避免信息处理部分计算过程公式化简中出现近似处理,使得到的血氧饱和度计算公式更加精确。在进行校准时分别使用800nm和660nm、800nm和940nm、660和940nm光源,在实际测量时仍使用660nm和940nm光源。本发明的优点是:获得更加精确的系统参数,使血氧含量的测定值更加准确。

Description

一种增加参考光源校准的脉搏血氧监控系统及方法
技术领域
本发明涉及一种增加参考光源校准的脉搏血氧监控系统及方法,用于脉搏血氧监控,属于脉搏血氧监控技术领域。
背景技术
脉搏血氧监控系统,一般指指夹式血氧仪,是一种通过测定手指的光传导强度获得动脉血的光吸收参数,进行血氧含量计算的测量系统。因其无创、测量方便、可实时监控的特点,被广泛应用与临床和运动监控领域。监控系统的传感器,即指夹式血氧探头,测量时固定在人体手指端,以手指作为血红蛋白的容器,将透过手指的光信号传递给信号处理系统。处理器计算吸光比率,并将结果与存储器中的饱和度对照表进行对应,得到相应的血氧饱和度。
血氧仪在投入市场前必须进行校准,获得正常的饱和度对照表,保证仪器达到目标精度。而经过一段时间的使用后,由于电磁干扰、精密零件移动、水气和尘埃的影响,血氧仪的饱和度对照表与真值会发生一定程度的偏移,当偏移超过精度,为保证测试准确需要再次进行仪器校准,更新饱和度对照表。因此,血氧仪的校准是保证测量精度的关键。
目前的血氧仪使用波长660nm的红光和940nm的近红外光作测试光源,运用血液中含氧血红蛋白(HbO2)和还原血红蛋白(Hb)对两束光的吸收特性不同的原理,通过分析手指分别对两束光的透光曲线,计算血氧饱和度(SpO2)的值。而在系统校准过程中,在计算时不可避免地对其中一束光的吸光特性进行了取0的近似操作,造成参数误差。
申请公布号为CN1107714050A的专利申请,三波长血氧饱和度检测方法与装置以及可穿戴设备,用于消除活动引起的噪声导致的检测误差。具体通过增加等吸收点波长光源,利用等吸收点波长光源通过光电容积描记法分别检测用户静止状态下和活动状态下的归一化交流光信号,将检测到的活动状态下的归一化交流光信号减去静止状态下的归一化交流光信号得到活动噪声信号。在实际测量中先减去活动噪声信号再进行血氧饱和度计算以减小误差。但本方法得到的血氧饱和度的测量仍然不精确。
发明内容
本发明的目的在于:提出了一种增加参考光源校准的脉搏血氧监控系统及方法,解决现有三波长血氧饱和度检测方法对血氧饱和度的测量仍然不精确的问题。
本发明采用的技术方案如下:
一种增加参考光源校准的脉搏血氧监控方法,其特征在于:
步骤1、获取三束交替工作的脉搏血氧信号,三束脉搏血氧信号为两束测试光源和一束参考光源;
步骤2、基于血氧饱和度原始计算公式,将获取的脉搏血氧信号做精确处理;
步骤3、输出精确处理的结果。
进一步,所述步骤1中,两束测试光源分别为660nm和940nm,参数光源为800nm。
进一步,所述步骤2中,精确处理后得到的血氧饱和度计算公式为:
Figure GDA0001873284370000021
其中,
Figure GDA0001873284370000022
式中A1,A2,B1,B2,B,C分别为:
Figure GDA0001873284370000023
Figure GDA0001873284370000024
进一步,一种增加参考光源校准的脉搏血氧监控系统包括,
血氧探头:血氧探头包括三个沿手指宽度方向并列设置的发光二极管和一个正对发光二极管的光敏二极管。信号获取部分主要获取三束交替工作的脉搏血氧信号,三束脉搏血氧信号为两束测试光源和一束参考光源;
信号处理器:基于血氧饱和度原始计算公式,将获取的脉搏血氧信号作精确处理;
显示模块:输出精确处理的结果。
综上所述,由于采用了上述技术方案,本发明的有益效果是:
1.本发明中,一种增加参考光源校准的脉搏血氧监控方法,在两束原始测量光源基础上增加一束参考光源,进行校准时三束光源交替工作。血氧探头获取三束脉搏血氧光信号后转化为电信号并传递给信号处理部分,信号处理部分在处理器中进行计算时,血氧饱和度原始计算公式化简过程中将避免
Figure GDA0001873284370000025
这一近似处理,通过精确的化简步骤得到精确的系统参数,使得在实际测量时血氧饱和度计算公式更加精确。
2.本发明中,血氧探头中三个发光二极管沿手指宽度方向紧密并列放置,总面积不超过接收透过手指的光信号的光敏二极管的1/2,减小因手指形状和发光二极管位置不同造成的误差。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,应当理解,以下附图仅示出了本发明的某些实施例,因此不应被看作是对范围的限定,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他相关的附图。
图1是脉搏血氧监控系统框图;
图2是增加参考光源后的血氧探头侧视图;
图3是发光二极管电路图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明,即所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。通常在此处附图中描述和示出的本发明实施例的组件可以以各种不同的配置来布置和设计。
因此,以下对在附图中提供的本发明的实施例的详细描述并非旨在限制要求保护的本发明的范围,而是仅仅表示本发明的选定实施例。基于本发明的实施例,本领域技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例,都属于本发明保护的范围。
需要说明的是,术语“第一”和“第二”等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
一种增加参考光源校准的脉搏血氧监控系统及方法,解决现有三波长血氧饱和度检测方法,对血氧饱和度检测方法对血氧饱和度的测量仍然不精确的问题。
采用的方法如下:
一种增加参考光源校准的脉搏血氧监控方法,其特征在于:
步骤1、获取三束交替工作的脉搏血氧信号,三束脉搏血氧信号为两束测试光源和一束参考光源;
步骤2、基于血氧饱和度计算公式,将获取的脉搏血氧信号做精确处理;
步骤3、输出精确处理的结果。
进一步,所述步骤1中,两束测试光源分别为660nm和940nm,参数光源为800nm。
进一步,所述步骤2中,精确处理后得到的血氧饱和度计算公式为:
Figure GDA0001873284370000041
其中,
Figure GDA0001873284370000042
进一步,一种增加参考光源校准的脉搏血氧监控包括,
血氧探头:获取三束交替工作的脉搏血氧信号,三束脉搏血氧信号为两束测试光源和一束参考光源,血氧探头上的发光二极管沿手指宽度方向并列设置;
信号处理器:基于血氧饱和度计算公式,将获取的脉搏血氧信号作精确处理;
显示模块:输出精确处理的结果。
以下结合实施例对本发明的特征和性能作进一步的详细描述,在校准时分别使用800nm和660nm、800nm和940nm、660和940nm光源,在实际测量时仍使用660nm和940nm光源。
如图1所示,一种增加参考光源校准的脉搏血氧监控主要分为人体信号获取和信号处理两个部分。信号获取部分主要由血氧探头组成,通过发光二极管发出光信号,处于发光二极管正对面的光敏二极管接收透过手指的光强度变化信号,并将光信号转化为电信号传递给信号处理部分。信号处理分别由滤波放大、采样和处理器分析计算几个步骤组成。最后由显示模块输出得到的血氧饱和度值。
如图2所示,本发明的指夹式血氧探头中包含三个发光二极管,在其正对面放置用于接收透过手指的光强度变化信号并将光信号转化为电信号的光敏二极管。为减小因手指形状和发光二极管位置不同造成的误差,三个发光二极管并列紧密放置且总面积不超过光敏二极管的1/2。
图3为三个发光二极管的电路图,发光波长分别为660nm,800nm和940nm。其中,660nm和940nm发光二极管两者反向并联,为实际使用时的测量管,800nm发光二极管仅用于系统校准。
在本实施案例中,三束光源分别采用λ1=660nm和λ2=940nm的原始测量光源和λ3=800nm的参考光源。
步骤1、获取脉搏血氧信号:血氧探头上的三个沿手指宽度方向发光二极管交替发光,光敏二极管接收透过手指的光强度变化信号并转化为电信号,然后传递到信号处理部分。
步骤2、在信号处理部分光敏二极管传递过来的电信号经过滤波放大和采样后通过处理器进行分析计算,计算过程如下:
双波长测量时,两束测量光源波长分别为λ1=660nm和λ2=940nm,两束光交替发光,血氧饱和度计算原始公式为:
Figure GDA0001873284370000051
公式化简:
做近似处理
Figure GDA0001873284370000052
得到:
Figure GDA0001873284370000053
参数化简:
Figure GDA0001873284370000054
最终化简结果为:
SpO2=A+B·R
因在化简过程中做了近似处理,参数A与B的值不精确。
为使系统参数更加精确,在使用两束波长分别为λ1=660nm和λ2=940nm作为测量光源的基础上增加参考光源,并且参考光源波长为λ3=800nm,三束光交替发光,
血氧饱和度计算原始公式为:
Figure GDA0001873284370000055
Figure GDA0001873284370000056
式中:
Figure GDA0001873284370000057
可直接消去,化简得到:
Figure GDA0001873284370000058
Figure GDA0001873284370000059
参数化简:
Figure GDA00018732843700000510
Figure GDA0001873284370000061
Figure GDA0001873284370000062
Figure GDA0001873284370000063
式中C为常数;
将上述所得参数带入双波长测量时原始计算公式得:
Figure GDA0001873284370000064
式中,a,b,c均为精确值,
其中:
Figure GDA0001873284370000065
c=C。
由上述计算过程可知,为将血氧饱和度原始计算公式化简,需要条件:
Figure GDA0001873284370000066
在双波长模型中,采用的两束测量光源λ1=660nm和λ2=940nm都无法满足上述条件,此时利用现有方法做近似处理
Figure GDA0001873284370000067
得到的系统参数A与B的值不精确。因此,我们增加λ3=800nm的参考光,含氧血红蛋白和还原血红蛋白对参考光吸光系数相同,满足
Figure GDA0001873284370000068
这一条件,原始计算公式通过精确化简,可以得到精确的系统参数。将所得精确参数带入双波长测量时原始计算公式得到精确的血氧饱和度计算公式。
得到上述精确的血氧饱和度计算公式后,在实际测量之前进行系统校准时,分别使用800nm和660nm、800nm和940nm、660和940nm光源:
1、交替导通660nm二极管和800nm二极管,对获得的一系列R值和实际的血氧饱和度(SpO2)值依据公式SpO2=A+B·R进行线性拟合,确定参数A1,B1
2、交替导通940nm二极管和800nm二极管,对获得的一系列R值和实际的血氧饱和度(SpO2)值依据公式SpO2=A+B·R进行线性拟合,确定参数A2,B2
3、交替导通660nm二极管和940nm二极管,对获得的一系列R值和对应的血氧饱和度(SpO2)值依据公式SpO2=A+B·R进行线性拟合,确定参数A,。
4、确定以上参数后,令常数
Figure GDA0001873284370000071
代入上述精确的血氧饱和度计算公式
Figure GDA0001873284370000072
通过此式,脉搏血氧饱和度监控系统在实际测量(采用660nm和940nm光源)时可得到更精确的血氧饱和度值。
步骤3、通过显示模块输出得到的血氧饱和度值。
本实施案例通过增加参考光源,避免处理器中计算过程的近似处理,优化血氧饱和度计算公式的化简过程,使得系统参数更加精确,最后得到精确的血氧饱和度计算公式。由此精确公式得到校准后仪器的饱和度对照表,此对照表更加精确,使系统实际测得的血氧饱和度准确度更高。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神原
则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。

Claims (4)

1.一种增加参考光源校准的脉搏血氧监控方法,其特征在于:
步骤1、获取三束交替工作的脉搏血氧信号,三束脉搏血氧信号为两束测试光源和一束参考光源;
步骤2、基于血氧饱和度计算原始公式,将获取的脉搏血氧信号做精确处理;
步骤3、输出精确处理的结果;
所述步骤2中,精确处理后得到的血氧饱和度计算公式为:
Figure FDA0002991203630000011
其中,
Figure FDA0002991203630000012
式中A1、A2、B1、B2分别为:
Figure FDA0002991203630000013
Figure FDA0002991203630000014
2.根据权利要求1所述的一种增加参考光源校准的脉搏血氧监控方法:其特征在于:所述步骤1中,两束测试光源波长分别为660nm和940nm,参数光源波长为800nm。
3.一种增加参考光源校准的脉搏血氧监控系统,其特征在于:包括,
血氧探头:获取三束交替工作的脉搏血氧信号,三束脉搏血氧信号为两束测试光源和一束参考光源;
信号处理器:基于血氧饱和度计算公式,将获取的脉搏血氧信号作精确处理;
显示模块:输出精确处理的结果;
精确处理后得到的血氧饱和度计算公式为:
Figure FDA0002991203630000015
其中,
Figure FDA0002991203630000016
式中A1、A2、B1、B2分别为:
Figure FDA0002991203630000021
Figure FDA0002991203630000022
4.根据权利要求3所述的一种增加参考光源校准的脉搏血氧监控系统,其特征在于:所述血氧探头包括一个光敏二极管和三个发光二极管,三个发光二极管沿手指宽度方向并列设置,发光二极管的总面积不超过光敏二极管总面积的1/2。
CN201810587483.8A 2018-06-06 2018-06-06 一种增加参考光源校准的脉搏血氧监控系统及方法 Active CN109157224B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810587483.8A CN109157224B (zh) 2018-06-06 2018-06-06 一种增加参考光源校准的脉搏血氧监控系统及方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810587483.8A CN109157224B (zh) 2018-06-06 2018-06-06 一种增加参考光源校准的脉搏血氧监控系统及方法

Publications (2)

Publication Number Publication Date
CN109157224A CN109157224A (zh) 2019-01-08
CN109157224B true CN109157224B (zh) 2021-06-08

Family

ID=64897123

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810587483.8A Active CN109157224B (zh) 2018-06-06 2018-06-06 一种增加参考光源校准的脉搏血氧监控系统及方法

Country Status (1)

Country Link
CN (1) CN109157224B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112858196A (zh) * 2020-10-20 2021-05-28 重庆大学 一种三波长静脉血氧浓度测量的方法

Family Cites Families (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2328279B (en) * 1997-08-12 2001-10-10 Abbott Lab Optical glucose detector
JP3763687B2 (ja) * 1998-12-25 2006-04-05 三井金属鉱業株式会社 血糖値測定装置
US7343186B2 (en) * 2004-07-07 2008-03-11 Masimo Laboratories, Inc. Multi-wavelength physiological monitor
US9861305B1 (en) * 2006-10-12 2018-01-09 Masimo Corporation Method and apparatus for calibration to reduce coupling between signals in a measurement system
CN101347334B (zh) * 2007-07-19 2012-09-05 深圳迈瑞生物医疗电子股份有限公司 血氧饱和度测量方法和装置
US10092226B2 (en) * 2011-12-23 2018-10-09 General Electric Company Method, arrangement, sensor, and computer program product for non-invasively measuring hemoglobin concentrations in blood
CN102755157B (zh) * 2012-07-17 2014-11-05 桂林电子科技大学 一种指端血流速度的测量方法
CN103876747B (zh) * 2014-01-10 2015-11-18 中国计量科学研究院 血氧饱和度模拟仪校准装置
CN104068865B (zh) * 2014-04-24 2017-05-31 辛勤 一种测量血氧饱和度的方法以及便携式设备
CN104257390B (zh) * 2014-09-04 2016-04-20 深圳市前海安测信息技术有限公司 无创血糖测定方法及系统
CN106031639A (zh) * 2015-03-07 2016-10-19 上海贝瑞电子科技有限公司 一种高精度血氧检测管理方法及系统
CN104757949A (zh) * 2015-03-31 2015-07-08 电子科技大学 可穿戴设备
US10595777B2 (en) * 2016-05-06 2020-03-24 Lvl Technologies, Inc. Method and device for hydration monitoring
CN105943056A (zh) * 2016-05-19 2016-09-21 苏州阿尔星科电子科技有限公司 新生儿经皮黄疸仪、客户端
CN107714050B (zh) * 2017-10-13 2021-01-15 广东乐心医疗电子股份有限公司 一种三波长血氧饱和度检测方法与装置以及可穿戴设备

Also Published As

Publication number Publication date
CN109157224A (zh) 2019-01-08

Similar Documents

Publication Publication Date Title
EP2624754B1 (en) System, method and computer program product for optical measurement of blood parameters
US4714080A (en) Method and apparatus for noninvasive monitoring of arterial blood oxygen saturation
KR100376649B1 (ko) 저포화도에대해최적화된산소포화도측정기및센서
US4796636A (en) Noninvasive reflectance oximeter
US5842979A (en) Method and apparatus for improved photoplethysmographic monitoring of oxyhemoglobin, deoxyhemoglobin, carboxyhemoglobin and methemoglobin
US9380969B2 (en) Systems and methods for varying a sampling rate of a signal
US6801797B2 (en) Pulse oximeter sensor with piece-wise function
US4854699A (en) Backscatter oximeter
KR102033914B1 (ko) 혈당 측정방법 및 이를 이용한 인체착용형 혈당 측정장치
US9642576B2 (en) Cancellation of light shunting
KR102166444B1 (ko) 근적외선을 이용한 비채혈식 혈당측정기 및 이를 이용한 혈당 측정 방법
JP6125821B2 (ja) 酸素飽和度測定装置及び酸素飽和度算出方法
US11350861B2 (en) System and method for a non-invasive medical sensor
JP6373511B2 (ja) 光学的分析システム及び方法
CN109157224B (zh) 一种增加参考光源校准的脉搏血氧监控系统及方法
US20180317825A1 (en) Device and method for measuring the concentration of a chemical compound in blood
US20140187884A1 (en) Systems and methods for ensemble averaging in pulse oximetry
US20090171172A1 (en) Method and system for pulse gating
US20240023844A1 (en) Infant blood oxygen saturation monitoring method and intelligent monitoring device
JP6741485B2 (ja) パルスフォトメータ、および血中吸光物質濃度の算出値の信頼性評価方法
WO2019005577A1 (en) PHYSIOLOGICAL MONITORING AND ASSOCIATED METHODS
Ajith et al. Integration of low cost SpO2 sensor in a wearable monitor
JP7091090B2 (ja) パルスオキシメータ及び血液特性測定装置
Gazivoda et al. Measurement firmware for pulse oximetry sensor

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