用光谱数据分析进行无创的光学血糖检测的方法和系统.pdf

上传人:b*** 文档编号:6563112 上传时间:2019-09-02 格式:PDF 页数:17 大小:2.37MB
返回 下载 相关 举报
摘要
申请专利号:

CN201210419811.6

申请日:

20090417

公开号:

CN102988061B

公开日:

20150401

当前法律状态:

有效性:

有效

法律详情:

IPC分类号:

A61B5/1455

主分类号:

A61B5/1455

申请人:

密苏里大学董事会

发明人:

徐智

地址:

美国密苏里州

优先权:

61/055,303,61/089,152

专利代理机构:

中国国际贸易促进委员会专利商标事务所

代理人:

秦晨

PDF下载: PDF下载
内容摘要

公开了一种基于光谱数据无创测量生物样品中的血糖水平的系统和方法。公开了在光谱数据的采集中提高信噪比和计算由于样品中的血液引起的光的衰减率的各种技术。公开的技术包括:(1)与对数函数结合使用标准差运算,(2)使用归一化因数,(3)使用比例因数,(4)解决温度对例如电阻器的各种系统部件的影响,(5)通过执行校准来解决光检测器中的暗电流。

权利要求书

1.一种用于检测生物样品中的血糖的方法,包括:利用被配置为入射在样品的目标区域上的至少一个光源;利用被定位为接收来自所述至少一个光源的光并产生指示被检测的光的功率的输出信号的至少一个光检测器,该输出信号具有与时间相关的电流;用处理器从所述至少一个光检测器接收输出信号;并基于接收到的输出信号用处理器以归一化因数或比例因数计算由于存在于目标区域中的样品中的血液引起的衰减率;用处理器通过计算与时间相关的输出电流的对数的标准差来消除由所述至少一个光检测器的与温度相关的检测器响应所引起的不确定性的影响,其中所述与时间相关的输出电流是由来自所述生物样品的相同的目标区域的光功率产生的;以及用处理器基于计算的衰减率确定与存在于目标区域中的样品相关的血糖水平。 2.根据权利要求1所述的用于检测生物样品中的血糖的方法,其中,所述至少一个光检测器包括具有反馈电阻器的前置放大器。 3.根据权利要求2所述的用于检测生物样品中的血糖的方法,进一步包括用处理器计算归一化因数Q(C,T)和比例因数Y(C,T)中的至少一个:基于作为时间的函数的第i个前置放大器的输出电压V(t)和标准差σ根据公式计算归一化因数Q(C,T),或者基于作为时间的函数的第i个前置放大器的输出电压V(t)以及第j个前置放大器的输出电压V(t)和标准差σ根据公式计算比例因数Y(C,T),其中,T是生物样品的温度,C是生物样品中的血糖浓度。 4.根据权利要求2所述的用于检测生物样品中的血糖的方法,进一步包括使用具有数字化的电压输出的模数转换器。 5.根据权利要求4所述的用于检测生物样品中的血糖的方法,进一步包括用处理器计算归一化因数Q(C,T)和比例因数Y(C,T)中的至少一个:基于第i个模数转换器的电压输出Δ(ADC)根据公式计算归一化因数Q(C,T),其中σ是标准差,或者基于第i个模数转换器的电压输出Δ(ADC)以及第j个模数转换器的电压输出Δ(ADC)根据公式计算比例因数Y(C,T),其中σ是标准差其中,T是生物样品的温度,C是生物样品中的血糖浓度。 6.一种用于检测生物样品中的血糖的系统,包括:至少一个光源,所述至少一个光源被配置为入射在样品的目标区域上;至少一个光检测器,所述至少一个光检测器被定位为接收来自所述至少一个光源的光并产生指示被检测的光的功率的输出信号,该输出信号具有与时间相关的电流;处理器,该处理器从所述至少一个光检测器接收输出信号,基于接收到的输出信号以归一化因数或比例因数计算由于存在于目标区域中的样品中的血液引起的衰减率,通过计算与时间相关的输出电流的对数的标准差来消除由所述至少一个光检测器的与温度相关的检测器响应所引起的不确定性的影响,并且基于计算的衰减率确定与存在于目标区域中的样品相关的血糖水平,其中所述与时间相关的输出电流是由来自所述生物样品的相同的目标区域的光功率产生的。 7.根据权利要求6所述的用于检测生物样品中的血糖的系统,其中,所述至少一个光检测器包括具有反馈电阻器的前置放大器。 8.根据权利要求7所述的用于检测生物样品中的血糖的系统,其中,处理器计算归一化因数Q(C,T)和比例因数Y(C,T)中的至少一个:基于作为时间的函数的第i个前置放大器的输出电压V(t)和标准差σ根据公式计算归一化因数Q(C,T),或者基于作为时间的函数的第i个前置放大器的输出电压V(t)以及第j个前置放大器的输出电压V(t)和标准差σ根据公式计算比例因数Y(C,T),其中,T是生物样品的温度,C是生物样品中的血糖浓度。 9.根据权利要求6所述的用于检测生物样品中的血糖的系统,进一步包括具有数字化的电压输出的模数转换器。 10.根据权利要求9所述的用于检测生物样品中的血糖的系统,其中,处理器计算归一化因数Q(C,T)和比例因数Y(C,T)中的至少一个:基于第i个模数转换器的电压输出Δ(ADC)根据公式计算归一化因数Q(C,T),其中σ是标准差,或者基于第i个模数转换器的电压输出Δ(ADC)以及第j个模数转换器的电压输出Δ(ADC)根据公式计算比例因数Y(C,T),其中σ是标准差,其中,T是生物样品的温度,C是生物样品中的血糖浓度。

说明书

本申请是于2009年4月17日提交的、题为“用光谱数据分析进 行无创的光学血糖检测的方法和系统”的发明申请第200980126116.7 号的分案申请。

相关申请的交叉引用

本专利申请要求2008年5月22日提交的申请号为61/055,303的 美国临时专利申请的优先权,该美国临时专利申请的公开内容通过引 用被合并在此,本专利申请还要求2008年8月15日提交的申请号为 61/089,152的美国临时专利申请的优先权,该美国临时专利申请的公 开内容通过引用也被合并在此。

背景技术

糖尿病是一种慢性疾病,当不对其进行控制的时候,随着时间的 流逝会对许多身体系统(包括神经、血管、眼睛、肾脏和心脏)造成 严重的损害。美国国家糖尿病、消化系统病、肾疾病研究所 (NIDDK)估计2007年在美国有2360万人或者美国总人口的7.8% 患有糖尿病。在全世界,世界卫生组织(WHO)估计超过1.8亿人 患有糖尿病,到2030年他们预计这个数字将增加到3.66亿,美国占 3030万。根据WHO,在2005年估计有110万人死于糖尿病。他们 预测2006年到2015年,糖尿病致死率将总的增加超过50%且在中 上收入国家增加超过80%。

糖尿病对个人和整个社会造成的经济负担是重大的。根据美国糖 尿病协会,2007年美国糖尿病总年度经济成本估计达到1740亿美 元。自2002年以来这个数字增加了420亿美元。这32%的增长意味 着美元数量每年提升超过了80亿美元。

糖尿病管理一个很重要的要素就是糖尿病患者在家庭环境下自我 监测血糖(SMBG)浓度。通过经常测试血糖水平,糖尿病患者能更 好地管理用药、饮食和运动以持续控制和预防长期不良健康结果。实 际上,糖尿病控制与并发症实验(DCCT)若干年来密切关注了1441 名糖尿病患者,结果表明,与标准治疗组比较,那些每天按照强化控 制方案进行多种血糖测试的患者中只有1/4的人发展为糖尿病性眼疾 病,1/2的人发展为肾脏疾病,1/3的人发展为神经性疾病,且已经 有这三种并发症早期形态的患者中极少数人恶化。

然而,由于在分析之前通过皮肤抽取血液不方便并且带来痛苦, 这使得许多糖尿病患者做不到为了很好地控制血糖而应当做的那样勤 奋,因此目前的监测技术阻碍了正常的使用。结果,无创测量葡萄糖 浓度是糖尿病管理的一种可取和有利的发展。无创监测将使每天的多 次测试没有痛苦且更加受到儿童糖尿病患者的欢迎。根据2005年刊 登的一篇研究(J,Wagner,C.Malchoff,and G.Abbott,Diabetes  Technology&Therapeutics(糖尿病技术和疗法),7(4)2005,612- 619),患有糖尿病的人利用无创血糖监测装置更频繁地执行 SMBG,并且已经改善生活质量。

存在许多用于血糖确定的无创方法。一种无创的血液化学检测技 术涉及收集和分析光谱数据。

由于在正被感测的区域中除了血液之外还存在其它成分(比如, 皮肤、脂肪、肌肉、骨头、胞间液),所以从根据光谱学获得的光谱 数据或其它数据提取关于血液特性的信息(例如葡萄糖浓度)是一个 复杂的问题。这样的其它成分能够以改变读数(reading)这样的方 式影响这些信号。具体地讲,作为结果而获得的信号的幅度可以比信 号中与血液对应的部分的幅度大很多,因此,限制了精确提取血液特 性信息的能力。

本发明旨在克服上述的一个或多个问题。

发明内容

在本发明的一方面中,公开了一种用于检测生物样品中的葡萄糖 的系统。该系统包括:至少一个光源,所述至少一个光源被配置为入 射在样品的目标区域上;至少一个光检测器,所述至少一个光检测器 被定位为接收来自所述至少一个光源的光并产生指示被检测的光的功 率的输出信号,该输出信号具有与时间相关的电流(time dependent  current);和处理器,该处理器从所述至少一个光检测器接收输出 信号,基于接收到的输出信号,以至少20:1的信噪比计算由于存在 于目标区域中的样品中的血液引起的衰减率,并基于计算的衰减率确 定与存在于目标区域中的样品相关的血糖水平。

在本发明的另一方面中,公开了一种用于检测生物样品中的葡萄 糖的系统。该系统包括:至少一个光源,所述至少一个光源被配置为 入射在样品的目标区域上;至少一个光检测器,所述至少一个光检测 器被定位为接收来自所述至少一个光源的光并产生指示被检测的光的 功率的输出信号,该输出信号具有与时间相关的电流;和处理器,该 处理器从所述至少一个光检测器接收输出信号,基于接收到的输出信 号,以归一化因数计算由于存在于目标区域中的样品中的血液引起的 衰减率,并基于计算的衰减率确定与存在于目标区域中的样品相关的 血糖水平。

在本发明的另一方面中,公开了一种用于检测生物样品中的葡萄 糖的系统。该系统包括:至少一个光源,所述至少一个光源被配置为 入射在样品的目标区域上;至少一个光检测器,所述至少一个光检测 器被定位为接收来自所述至少一个光源的光并产生指示被检测的光的 功率的输出信号,该输出信号具有与时间相关的电流;和处理器,该 处理器从所述至少一个光检测器接收输出信号,基于接收到的输出信 号,以比例因数计算由于存在于目标区域中的样品中的血液引起的衰 减率,并基于计算的衰减率确定与存在于目标区域中的样品相关的血 糖水平。

在本发明的另一方面中,公开了一种用于检测生物样品中的葡萄 糖的系统。该系统包括:至少一个光源,所述至少一个光源被配置为 入射在样品的目标区域上;至少一个光检测器,所述至少一个光检测 器被定位为接收来自所述至少一个光源的光并产生指示被检测的光的 功率的输出信号,该输出信号具有与时间相关的电流;和处理器,该 处理器从所述至少一个光检测器接收输出信号,基于接收到的输出信 号计算由于存在于目标区域中的样品中的血液引起的衰减率并消除由 所述至少一个光检测器的与温度相关的检测器响应所引起的不确定性 的影响,并且,基于计算的衰减率确定与存在于目标区域中的样品相 关的血糖水平。

在本发明的另一方面中,公开了一种用于检测生物样品中的葡萄 糖的系统。该系统包括:至少一个光源,所述至少一个光源被配置为 入射在样品的目标区域上;至少一个光检测器,所述至少一个光检测 器包括具有反馈电阻器的前置放大器,被定位为接收来自所述至少一 个光源的光并产生指示被检测的光的功率的输出信号,该输出信号具 有与时间相关的电流;和处理器,该处理器从所述至少一个光检测器 接收输出信号,基于接收到的输出信号计算由于存在于目标区域中的 样品中的血液引起的衰减率并消除由所述至少一个光检测器的与温度 相关的检测器响应所引起的不确定性的影响,并且,基于计算的衰减 率确定与存在于目标区域中的样品相关的血糖水平。

在本发明的另一方面中,公开了一种用于检测生物样品中的葡萄 糖的系统。该系统包括具有前置放大器和反馈电阻器的至少一个光检 测器。

在本发明的另一方面中,公开了一种用于检测生物样品中的葡萄 糖的方法。该方法包括:利用被配置为入射在样品的目标区域上的至 少一个光源;利用被定位为接收来自所述至少一个光源的光并产生指 示被检测的光的功率的输出信号的至少一个光检测器,该输出信号具 有与时间相关的电流;用处理器从所述至少一个光检测器接收输出信 号,并基于接收到的输出信号,以至少20:1的信噪比计算由于存在 于目标区域中的样品中的血液引起的衰减率;以及用处理器基于计算 的衰减率确定与存在于目标区域中的样品相关的血糖水平。

在本发明的另一方面中,公开了一种用于检测生物样品中的葡萄 糖的方法。该方法包括:利用被配置为入射在样品的目标区域上的至 少一个光源;利用被定位为接收来自所述至少一个光源的光并产生指 示被检测的光的功率的输出信号的至少一个光检测器,该输出信号具 有与时间相关的电流;用处理器从所述至少一个光检测器接收输出信 号,并基于接收到的输出信号,用处理器以归一化因数计算由于存在 于目标区域中的样品中的血液引起的衰减率;以及用处理器基于计算 的衰减率确定与存在于目标区域中的样品相关的血糖水平。

在本发明的另一方面中,公开了一种用于检测生物样品中的葡萄 糖的方法。该方法包括:利用被配置为入射在样品的目标区域上的至 少一个光源;利用被定位为接收来自所述至少一个光源的光并产生指 示被检测的光的功率的输出信号的至少一个光检测器,该输出信号具 有与时间相关的电流;用处理器从所述至少一个光检测器接收输出信 号;用处理器基于接收到的输出信号,以比例因数计算由于存在于目 标区域中的样品中的血液引起的衰减率;以及用处理器基于计算的衰 减率确定与存在于目标区域中的样品相关的血糖水平。

在本发明的另一方面中,公开了一种用于检测生物样品中的葡萄 糖的方法。该方法包括:利用被配置为入射在样品的目标区域上的至 少一个光源;利用被定位为接收来自所述至少一个光源的光并产生指 示被检测的光的功率的输出信号的至少一个光检测器,该输出信号具 有与时间相关的电流;用处理器从所述至少一个光检测器接收输出信 号,并用处理器基于接收到的输出信号,以比例因数计算由于存在于 目标区域中的样品中的血液引起的衰减率;用处理器消除由所述至少 一个光检测器的与温度相关的检测器响应所引起的不确定性的影响; 以及用处理器基于计算的衰减率确定与存在于目标区域中的样品相关 的血糖水平。

这些只是本发明的众多方面中的一些,且不应被看作是关于本发 明的众多方面的全部清单。

附图说明

为了更好地理解本发明,可以参考附图,其中:

图1示出根据示例性实施例的对应于动脉血的光吸收的脉冲波 图;

图2示出用于获得光谱数据的示例性系统;

图3示出用图1中的数据根据公式(9)计算的A(t)的图;以及

图4是在数字化之前将光电流转换为电压的前置放大器电路的基 本示意图。

具体实施方式

在下面的详细描述中,为了提供对本发明的全面的理解,阐述了 大量的示例性的具体细节。然而,本领域技术人员将会理解,本发明 可以在没有这些具体细节或者对所述细节的各种修改的情况下实施。 在其他示例中,没有详细描述众所周知的方法、工艺和部件,以避免 使本发明模糊。

可用光谱学确定被生物样品(例如人的手指)所吸收和散射(即 衰减)的光量。通过测量被样品吸收的光量,可无创地确定被检者的 葡萄糖、胆固醇和血红素水平。由于指尖中的毛细血管的密集度大并 且在指尖发生从动脉血到静脉血的转换,所以通常优选指尖测量。然 而,本发明的方法并不限于用指尖。例如,生物样品可以是人的耳 垂。

当光透射穿过生物样品(例如人的手指)时,光被手指的各种成 分(包括皮肤、肌肉、骨头、脂肪、胞间液和血液)所衰减。然而, 观察到,通过人的手指的光衰减表现出与心跳对应的小范围的循环模 式。人们相信,这种循环模式将存在于许多其他的身体部位(耳垂是 众多例子中的一个)的测量中。

图1描绘了检测器光电流ID(t)的图102,该光电流对应于在光穿 过被检者的手指后被检测器接收的光功率。可以看出,检测器光电流 表现出循环模式。该循环模式是由于被检者的心跳而引起的,且这种 心跳循环地增加并减少被检者的毛细血管(或其他结构)中的血量。 尽管该循环模式的大小相比于检测器产生的总光电流而言小,但是可 以从图102的循环模式中提取出大量的信息。例如,假定人的心率是 60次每分钟,即,任何脉搏的开始和结束之间的时间是一秒钟。在 这一秒钟的期间内,光电流将具有最大值或峰值读数104和最小值或 谷值读数106。图的峰值读数104对应于当毛细血管中有最小量的血 时,谷值读数106对应于当毛细血管中有最大量的血时。通过使用所 述循环图中的峰值和谷值提供的信息,排除不在毛细血管中的主要的 手指成分(例如皮肤、脂肪、骨头、肌肉、胞间液)所造成的光吸收 和散射。排除这些不在毛细血管中的主要成分是由于其在一个心跳的 时间间隔期间内不可能变化。换句话说,可基于图102的峰值和谷值 检测被血液吸收和散射(即,衰减)的光。

假定光感测装置产生的循环光电流的峰值是IP,循环光电流的相 邻的谷值是Iv,在没有人手指的情况下光感测装置产生的光电流是 I0,对应于峰值和谷值光电流的透射率可定义为:

T V = I V I 0 - - - ( 1 ) ; ]]>

T P = I P I 0 - - - ( 2 ) . ]]>

对应的峰值和谷值吸光率为:

Av=-log(Tv)    (3);

AP=-log(TP)    (4)。

Av和AP之间的差表示手指中的血液(不包括非血液成分)导致的光吸 收和散射:

ΔA = A V - A P = log ( I P I V ) - - - ( 5 ) . ]]>

在公式(5)示出的算法中可以看出,ΔA不依赖于I0。所以,计 算ΔA不需要确定在无样品时由光感测装置产生的电流。监测对应于 透过样品的光功率的光电流足以计算ΔA。

图2描绘了用在示例性实施例中的示例性设备的简化框图。总体 上用数字200表示的光学测量系统用“脉动”的概念确定仅被样品 (在本示例性实施例中是人的手指)中的血液吸收和散射的光量。电 源201(例如电池)为光源202提供电力,该光源202产生多条向着 被检者的手指的顶部引导的光束204、206、208、210。在一个示例 性实施例中,光束204、206、208、210中的每一条具有相同的波长 或不同的波长范围,典型地在800nm到1600nm之间。尽管此处将 光学测量系统200描述为产生四(4)条光束,但是在其它实施例 中,设想可将光源202改装成产生更少的光束或者另外的光束。

第一光圈212确保光束204、206、208、210入射到手指的目标 区域上。第二光圈214确保透过手指的光束的一部分入射到透镜216 上。光束204、206、208、210被手指和光学测量系统200的各部件 衰减,所以衰减光束218、220、222、224从手指发射。衰减光束 218、220、222、224入射到透镜216上,透镜216收集衰减光束 218、220、222、224以使其更有效地入射到检测器块226上。

检测器块226位于透镜216的正下方,并且包括多个光感测装置 (LSD)228、230、232、234,例如光电二极管阵列。根据光学测量 系统200的一个方面,光感测装置228、230、232、234的每个分别 检测由对应的干涉滤光器(IF)236、238、240、242限定的具体波 长的光。干涉滤光器透射一个或多个光谱带或光线,并阻挡其它光谱 带或光线。

光感测装置228、230、232、234的每个产生与被特定光感测装 置接收的光的功率成比例的对应的光电流信号。可将光电二极管所产 生的光电流信号转换为另一形式的信号,例如模拟电压信号或数字信 号。处理器243与检测器块226耦合,并被构造为计算光电流信号 244、246、248、250的变化。

根据一方面,处理器243执行例如公式(5)中显示的算法以计 算仅由手指中的血液引起的光吸收变化(ΔA)。其后,可使用对血 液的光吸收的这个量化计算来确定血液的特性。例如,通过将计算的 光吸收值同与存储器(未显示)中存储的不同葡萄糖水平对应的预定 值进行比较,可确定被检者的血糖水平。

与基于手指的脉动检测方法相关的困难在于低的信噪(S/N) 比,因为循环模式的振幅(即,峰值和谷值之间的差)通常是由透过 手指的光功率产生的总光电流的1%-2%。为了在确定ΔA时获得 100:1的S/N比,在10Hz的带宽内,在吸光率(峰值到峰值)中用 于测量手指造成的光吸收的装置的基线噪声不应大于3.0×10-5。

然而,以一些电池供电的手持无创的血液化学测量装置使用的低 的光功率水平,难于获得在10Hz的带宽内3.0×10-5吸光率(峰值到峰 值)基线噪声水平。一种解决方案涉及数据取平均。为了提高S/N 比,下式定义的ΔA的平均值用于进一步的计算以提取血糖浓度:

ΔA ‾ = Σ j = 1 M ΔA j - - - ( 6 ) . ]]>

在公式(6)中,M是脉动测量的时间间隔期间的心跳的次数。 然而,由于心率在每秒钟一次的数量级上的事实而导致该方法需要较 长的数据采集时间。例如,为了使S/N比提高到五倍,将需要25秒 钟,为了使S/N比提高到十倍,将需要100秒钟。相较而言,现在的 商用的抽血血糖仪能在5秒钟内确定血糖水平。另外,较长的检测时 间会显著增加由于手指移动、光功率漂移、装置温度变化等引起的测 量误差。所以,需要新的技术来快速而准确地测量血糖水平。

通过标准差提高S/N比

图1中所示的与时间相关的检测器光电流输出ID(t)可表示为与心 跳相对应的小的与时间相关的循环光电流ΔI(t)、噪声电流n(t)和恒定 的基线光电流IB的和:

ID(t)=IB+ΔI(t)+n(t)    (7)。

上式可被重新整理为:

I D ( t ) I B = 1 + ΔI ( t ) + n ( T ) I B - - - ( 8 ) . ]]>

将常用的对数应用于公式(8)的两边,得到:

A ( t ) = log [ I D ( t ) I B ] = log ( 1 + ΔI ( t ) + n ( t ) I B ) - - - ( 9 ) . ]]>

总体上用数字300表示的图3示出了使用图1中的数据根据公式 (9)计算的典型的A(t)的图302。对于图3中示出的脉冲函数A(t), 在一个心跳的时间间隔期间存在下面的关键的关系:

σ[A(t)]=kΔA  (10),

其中σ[A(t)]是A(t)的标准差,k是比例常数。

考虑到IB是常数,且σ2(logIB)=0,得出:

σ[A(t)]=σ[logID(t)]  (11)。

所以,在一个心跳的时间间隔期间A(t)图的峰值到谷值的高度可 直接从ID(t)的对数的标准差中得到:

ΔA = σ [ A ( t ) ] k = σ [ log I D ( t ) ] k - - - ( 12 ) . ]]>

公式(12)的一个主要的好处是如下所解释的可在较短的数据采集时 间(约为一秒)内得到高的S/N比。

在图2所描绘的基于手指的脉动测量中,ΔI(t)+n(t)的和值通常小 于大的恒定的基线光电流IB的2%。所以,公式(9)可近似为:

A ( t ) = log [ I D ( t ) I B ] ≈ 1 ln 10 ΔI ( t ) + n ( t ) I B - - - ( 13 ) . ]]>

类似地,A(t)的标准差可近似为:

σ [ A ( t ) ] ≈ 1 ln 10 σ 2 [ ΔI ( t ) ] + σ 2 [ n ( t ) ] I B - - - ( 14 ) . ]]>

公式(14)说明了公式(12)的大的噪声降低功率。例如,对于 比值为(或10%)的相对高的基线噪声,基线噪声n(t) 对于σ[A(t)]的贡献估计少于0.005(或0.5%),其对应于在没有增加 检测时间的情况下S/N比增加到20倍。照此,可在不增加数据采集 时间的情况下得到惊人的噪声降低,可在一个心跳(近似为一秒)的 时间间隔内完成基于手指的脉动测量,且对于在确定ΔA时100:1的 S/N比的要求可以使用在10Hz的带宽内具有大约6.0×10-4吸光率(峰 值到峰值)的基线噪声的光学系统得到满足。应当指出的是,当由于 低的光照功率使得光学系统的基线噪声由散粒噪声主导时,噪声降低 到20倍等同于光照功率增加到202=400倍。

在很短的数据采集时间(例如少于一秒)内获得较高的S/N比的 这种能力将显著地减少由于测量期间的诸如手指移动、温度变化和光 功率漂移的因素引起的检测误差,所以惊人地提高脉动检测方法的准 确性和可再现性。

另外,因为所有波长的传输光由于心跳而具有相同的脉冲形状, 因此k的值不会随波长改变。结果是,常数k将在下面部分讨论的 数据归一化中被消除,且σ[logID(t)]将被用于进一步的回归分析以在光 学测量和血糖水平之间建立关联。因为σ[logID(t)]仅包含两个在大多数 普及的电子数据表程序(例如微软)中可用的标准数学函 数,所以这将大大简化了数据分析过程。是微软公司的联 邦注册商标,微软公司的办公地点在One Microsoft Way,Redmond, Washington 98052-6399。

归一化

在每个波长λi处,通过下面的公式将吸收ΔA(λi)与由于心跳导致 的指尖的光学感测区域中的血液量的增加(ΔB)相联系:

ΔA(λi)=ε(C,λi,T)ΔB    (15),

其中ε(C,λi,T)是血液在波长为λi、手指温度为T且血糖浓度为C时的 吸收/散射系数。众所周知,变量ΔB因人而异,且甚至对于同一个人 可能每天也会不同。

变量ΔB导致的不确定性可通过在每个波长λi处引入如下面的公式 所定义的归一化因数Qi(C,T)而被消除:

Q i ( C , T ) = ΔA ( λ i ) Σ i = 1 N ΔA ( λ i ) = ϵ ( C , λ i , T ) Σ i = 1 N ϵ ( C , λ i , T ) - - - ( 16 ) , ]]>

其中N是使用的波长的总数。优选地,N的范围通常为20到30。

基于公式(12)和(16),通过下面的公式将Qi(C,T)与每个波长 λi处的检测器光电流ID(λi,t)相联系:

Q i ( C , T ) = ΔA ( λ i ) Σ i = 1 N ΔA ( λ i ) = σ [ log I D ( λ i , t ) ] / k Σ i = 1 N σ [ log I D ( λ i , t ) ] / k = σ [ log I D ( λ i , t ) ] Σ i = 1 N σ [ log I D ( λ i , t ) ] - - - ( 17 ) , ]]>

如公式(17)所示,常数k被消除,且将σ[logID(t)]用于进一步的 回归分析以在光学测量和血糖水平之间建立关联。这可能是因为同时 从所有的检测通道(detection channel)获取数据。

光学测量和血糖浓度之间的关联可根据下式建立:

C optical = Σ i = 1 N a i ( T ) Q i ( C , T ) - - - ( 18 ) , ]]>

其中Coptical是通过光学测量预测的血糖浓度,Qi(C,T)由公式(16)和 (17)定义,且ai(T)是对应于波长λi的与温度相关的回归系数。ai(T) 的值可用适当的统计方法(例如,偏最小二乘(PLS)回归)提取。

公式(18)表示可在不同的手指温度进行大量的校准时的理想情 况。实际上,通常只能进行有限数量(例如,15到20)的校准,且 每次校准可能在不同的手指温度进行。在这种条件下,手指温度可被 视为独立的变量,且上面的公式可近似为:

C optical = Σ i = 1 N b i Q i ( C , T ) + ηT - - - ( 19 ) , ]]>

其中bi是对应于波长λi的与温度无关的回归系数,η是手指温度的回 归系数。bi和η的值可用适当的统计方法(例如,偏最小二乘 (PLS)回归)提取。

比例方法

可替代地,来自变量ΔB的不确定性可通过在波长λi处引入比例因 数Yij而被消除:

Y ij ( C , T ) = ΔA ( λ i ) ΔA ( λ j ) = ϵ ( C , λ i , T ) ϵ ( C , λ j , T ) = σ [ log I D ( λ i , t ) ] σ [ log I D ( λ j , t ) ] - - - ( 20 ) , ]]>

其中j可为从1到N的任意数,假定该装置在所有N个波长处收集 信号。

和前面讨论的归一化算法类似,光学测量和血糖水平之间的关联 可根据下式建立:

C optical = Σ i ≠ j N f i ( T ) Y ij ( C , T ) - - - ( 21 ) , ]]>

其中Coptical是通过光学测量预测的血糖浓度,Yij(C,T)由公式(20)定 义,且fi(T)是对应于波长λi的与温度相关的的回归系数。fi(T)的值可 用适当的统计方法(例如,偏最小二乘(PLS)回归)提取。

公式(21)表示可在不同的手指温度进行大量的校准时的理想情 况。实际上,通常只能进行有限数量(例如,15到20)的校准,且 每次校准可能在不同的手指温度进行。在这种条件下,手指温度可被 视为独立的变量,且上面的公式可近似为:

C optical = Σ i ≠ j N h i Y ij ( C , T ) + βT - - - ( 22 ) , ]]>

其中hi是对应于波长λi的与温度无关的回归系数,β是手指温度的回 归系数。hi和β的值可用适当的统计方法(例如,偏最小二乘 (PLS)回归)提取。

与温度相关的装置响应的影响的消除

众所周知,硅光电二极管检测器的检测器灵敏度是波长和温度的 函数。对于图2中所示的总体上由数字200表示的装置配置,由第i 个硅二极管检测器接收的对应于波长λi的光功率根据下式被转换为光 电流:

ID(λi,t)=P(λi,t)S0(λi)[1+γ(λi)(TDi(t)-25°C)]  (23)。

上面的公式(23)中,P(λi,t)是由检测器接收的光功率,S0(λi)是 在波长λi和25℃处的检测器的光灵敏度,γ(λi)是波长λi处的光灵敏度 的温度系数,TDi(t)是第i个光电二极管检测器的温度。温度系数γ(λi) 随波长变化。例如,对于Hamamatsu的S1337系列光电二极管检测 器,γ(λi)的范围从900nm处接近零到1100nm处超过1.0%/℃。这给 图2所示的装置配置施加了一个潜在的问题,因为很难在正常的家庭 /办公环境下在糖尿病人使用的手持设备中保持每个单独的二极管检 测器的温度恒定不变。

由于检测器温度TDi(t)引起的这一不确定性可用公式(11)和 (12)所示的算法来被消除。将常用的对数应用到公式(23)的两 边,得到:

考虑到S0(λi)是常数且检测器温度TDi(t)在约为一秒的很短的数据 采集时间间隔期间几乎保持恒定不变的事实,得到:

σ[logID(λi,t)]=σ[logP(λi,t)]    (25)。

如此,检测器温度TDi(t)所导致的不确定性可通过使用该标准差方法来 被消除。

电压检测模式

在图2所示的装置配置中,通常,在数字化之前用前置放大器将 第i个光电二极管检测器的光电流ID(λi,t)转换为电压。图4示出总体 上用数字400表示的典型的前置放大器的电路示意图。

与第i个光电二极管检测器408耦合的第i个前置放大器400的 输出电压412可表达为:

Vi(t)=RiID(λi,t)=R0i[1+χi(TRi(t)-25°C)]ID(λi,t)    (26)。

在上面的公式(26)中,R0i是在25℃处第i个前置放大器的反 馈电阻器402的电阻值,χi是电阻器的温度系数,TRi(t)是电阻器的温 度。将常用的对数应用到公式(26)的两边,得到:

logVi(t)=logRi0+log[1+χi(TRi(t)-25°C)]+logID(λi,t)    (27)。

考虑到R0i是常数且电阻器温度TRi(t)在约为一秒的很短的数据采集 时间间隔期间不变的事实,得到:

σ[logVi(t)]=σ[logID(λi,t)]    (28)。

将公式(25)代入公式(28),得到:

σ[logVi(t)]=σ[logP(λi,t)]    (29)。

如此,电阻器温度TR(t)所导致的不确定性被消除。

在电压检测模式下,公式(17)中的归一化因数可表达为:

Q i ( C , T ) = σ [ log V i ( t ) ] Σ i = 1 N σ [ log V i ( t ) ] - - - ( 30 ) . ]]>

然后,在相应的校准条件下,光学测量和血糖浓度之间的数学关 联可根据公式(18)或公式(19)建立。

类似地,由公式(20)定义的归一化因数可表达为:

Y ij ( C , T ) = σ [ log V i ( t ) ] [ log V j ( t ) ] - - - ( 31 ) . ]]>

然后,在相应的校准条件下,光学测量和血糖浓度之间的数学 关联可根据公式(21)或公式(22)建立。典型的前置放大器400的 电路示意图还包括反馈电容器404、运算放大器406和接地线 (ground connection)410。

数字化

来自前置放大器400的输出电压412通常使用模数转换器 (ADC)被数字化。然后,将数字化的信号发送至计算机以进行数 据分析。与和在波长λi处收集光功率的第i个光电二极管408耦合的 第i个前置放大器通信的第i个ADC的输出可表达为下式:

在上面的公式(32)中,(ADC)0i是第i个ADC的偏移量 (shift),Gi是检测期间使用的标称的ADC增益,ID(λi,t)是第i个光 电二极管检测器的光电流,IDark,i是第i个光电二极管检测器的暗电 流,是第i个前置放大器的反馈电阻器的电 阻,Aoi是第i个前置放大器的偏移量。

(ADC)0i、IDark,i和Aoi这三个因数的贡献可通过正好在相应的手指测 量之前或之后切断光源而实施暗测量来被去除。当切断电源时,上面 的公式(32)变为:

(ADC)Dark,i=(ADC)0i+Gi(IDark,iRi+A01)    (33)。

上面的两个公式(32)和(33)之间的差反映了对应于光电流的 ADC输出:

Δ(ADC)i=(ADC)i-(ADC)Dark,i=GiID(λi,t)Ri    (34)。

将常用的对数应用到公式(34)的两边,得到:

logΔ(ADC)i=logGi+logID(λi,t)+logRi    (35)。

只要手指测量和暗测量之间的时间间隔短,Gi和Ri就可以被视为 常数。如此,得到:

σ[logΔ(ADC)i]=σ[logID(λi,t)]    (36)。

将公式(25)代入公式(36),进一步得到:

σ[logΔ(ADC)i]=σ[logP(λi,t)]    (37)。

基于公式(36),由公式(17)定义的归一化因数可表达为:

Q i ( C,T ) = σ [ log Δ ( ADC ) i ] Σ i = 1 N σ [ log Δ ( ADC ) i ] - - - ( 38 ) . ]]>

然后,在相应的校准条件下,光学测量和血糖浓度之间的数学关 联可根据公式(18)或(19)建立。

和归一化类似,由公式(20)定义的比例因数可表达为:

Y ij ( C,T ) = σ [ log Δ ( ADC ) j ] σ [ log Δ ( ADC ) j ] - - - ( 39 ) . ]]>

然后,在相应的校准条件下,光学测量和血糖浓度之间的关联可 根据公式(21)或(22)建立。

因此,已经示出和描述了新颖的发明的一些实施例。从前面的描 述中可以明显看出,本发明的某些方面不局限于本文示出的例子的特 定细节,所以,本领域技术人员将会设想到其它的修改和应用或其等 同物。在前面的说明书中使用的术语“具有”、“包括”和“包含” 及其相似的术语取“可选的”或者“可包括”而不是“必需”的含 义。然而,在考虑到说明书和附图后,本发明的很多改变、修改、变 型以及其它的使用和应用对于本领域技术人员来说将是显而易见的。 不脱离本发明的范围和精神的所有这些改变、修改、变型以及其它的 使用和应用被视为由本发明所覆盖,本发明仅仅由下述的权利要求所 限定。应当理解,本文公开的实施例包括在任何从属权利要求中描述 的特征的任何组合和所有组合。

用光谱数据分析进行无创的光学血糖检测的方法和系统.pdf_第1页
第1页 / 共17页
用光谱数据分析进行无创的光学血糖检测的方法和系统.pdf_第2页
第2页 / 共17页
用光谱数据分析进行无创的光学血糖检测的方法和系统.pdf_第3页
第3页 / 共17页
点击查看更多>>
资源描述

《用光谱数据分析进行无创的光学血糖检测的方法和系统.pdf》由会员分享,可在线阅读,更多相关《用光谱数据分析进行无创的光学血糖检测的方法和系统.pdf(17页珍藏版)》请在专利查询网上搜索。

1、(10)授权公告号 (45)授权公告日 (21)申请号 201210419811.6 (22)申请日 2009.04.17 61/055,303 2008.05.22 US 61/089,152 2008.08.15 US 200980126116.7 2009.04.17 A61B 5/1455(2006.01) (73)专利权人 密苏里大学董事会 地址 美国密苏里州 (72)发明人 徐智 (74)专利代理机构 中国国际贸易促进委员会专 利商标事务所 11038 代理人 秦晨 (54) 发明名称 用光谱数据分析进行无创的光学血糖检测的 方法和系统 (57) 摘要 公开了一种基于光谱数据无创测。

2、量生物样品 中的血糖水平的系统和方法。公开了在光谱数据 的采集中提高信噪比和计算由于样品中的血液引 起的光的衰减率的各种技术。公开的技术包括 : (1) 与对数函数结合使用标准差运算,(2) 使用归 一化因数,(3) 使用比例因数,(4) 解决温度对例 如电阻器的各种系统部件的影响,(5) 通过执行校 准来解决光检测器中的暗电流。 (30)优先权数据 (62)分案原申请数据 (51)Int.Cl. 审查员 熊狮 (19)中华人民共和国国家知识产权局 (12)发明专利 权利要求书2页 说明书11页 附图3页 (10)授权公告号 CN 102988061 B (45)授权公告日 2015.04.0。

3、1 CN 102988061 B 1/2 页 2 1. 一种用于检测生物样品中的血糖的方法, 包括 : 利用被配置为入射在样品的目标区域上的至少一个光源 ; 利用被定位为接收来自所述至少一个光源的光并产生指示被检测的光的功率的输出 信号的至少一个光检测器, 该输出信号具有与时间相关的电流 ; 用处理器从所述至少一个光检测器接收输出信号 ; 并基于接收到的输出信号用处理器 以归一化因数或比例因数计算由于存在于目标区域中的样品中的血液引起的衰减率 ; 用处理器通过计算与时间相关的输出电流的对数的标准差来消除由所述至少一个光 检测器的与温度相关的检测器响应所引起的不确定性的影响, 其中所述与时间相关。

4、的输出 电流是由来自所述生物样品的相同的目标区域的光功率产生的 ; 以及 用处理器基于计算的衰减率确定与存在于目标区域中的样品相关的血糖水平。 2. 根据权利要求 1 所述的用于检测生物样品中的血糖的方法, 其中, 所述至少一个光 检测器包括具有反馈电阻器的前置放大器。 3. 根据权利要求 2 所述的用于检测生物样品中的血糖的方法, 进一步包括用处理器计 算归一化因数 Qi(C, T) 和比例因数 Yij(C, T) 中的至少一个 : 基于作为时间的函数的第 i 个 前置放大器的输出电压 Vi(t) 和标准差 根据公式计算归一化因数 Qi(C, T), 或者基于作为时间的函数的第 i 个前置放。

5、大器的输出电压 Vi(t) 以及第 j 个前置 放大器的输出电压Vj(t)和标准差根据公式计算比例因数Yij(C, T), 其中, T 是生物样品的温度, C 是生物样品中的血糖浓度。 4. 根据权利要求 2 所述的用于检测生物样品中的血糖的方法, 进一步包括使用具有数 字化的电压输出的模数转换器。 5. 根据权利要求 4 所述的用于检测生物样品中的血糖的方法, 进一步包括用处理器计 算归一化因数 Qi(C, T) 和比例因数 Yij(C, T) 中的至少一个 : 基于第 i 个模数转换器的电压 输出 (ADC)i根据公式计算归一化因数 Qi(C, T), 其中 是标 准差, 或者基于第 i 。

6、个模数转换器的电压输出 (ADC)i以及第 j 个模数转换器的电压输出 (ADC)i根据公式计算比例因数 Yij(C, T), 其中 是标准差 其中, T 是生物样品的温度, C 是生物样品中的血糖浓度。 6. 一种用于检测生物样品中的血糖的系统, 包括 : 至少一个光源, 所述至少一个光源被配置为入射在样品的目标区域上 ; 至少一个光检测器, 所述至少一个光检测器被定位为接收来自所述至少一个光源的光 并产生指示被检测的光的功率的输出信号, 该输出信号具有与时间相关的电流 ; 处理器, 该处理器从所述至少一个光检测器接收输出信号, 基于接收到的输出信号以 权 利 要 求 书 CN 102988。

7、061 B 2 2/2 页 3 归一化因数或比例因数计算由于存在于目标区域中的样品中的血液引起的衰减率, 通过计 算与时间相关的输出电流的对数的标准差来消除由所述至少一个光检测器的与温度相关 的检测器响应所引起的不确定性的影响, 并且基于计算的衰减率确定与存在于目标区域中 的样品相关的血糖水平, 其中所述与时间相关的输出电流是由来自所述生物样品的相同的 目标区域的光功率产生的。 7. 根据权利要求 6 所述的用于检测生物样品中的血糖的系统, 其中, 所述至少一个光 检测器包括具有反馈电阻器的前置放大器。 8. 根据权利要求 7 所述的用于检测生物样品中的血糖的系统, 其中, 处理器计算归一 化。

8、因数 Qi(C, T) 和比例因数 Yij(C, T) 中的至少一个 : 基于作为时间的函数的第 i 个前置放 大器的输出电压 Vi(t) 和标准差 根据公式计算归一化因数 Qi(C, T), 或者基于作为时间的函数的第i个前置放大器的输出电压Vi(t)以及第j个前置放大器 的输出电压 Vj(t) 和标准差 根据公式计算比例因数 Yij(C, T), 其中, T 是生物样品的温度, C 是生物样品中的血糖浓度。 9. 根据权利要求 6 所述的用于检测生物样品中的血糖的系统, 进一步包括具有数字化 的电压输出的模数转换器。 10. 根据权利要求 9 所述的用于检测生物样品中的血糖的系统, 其中,。

9、 处理器计算归 一化因数 Qi(C, T) 和比例因数 Yij(C, T) 中的至少一个 : 基于第 i 个模数转换器的电压输 出 (ADC)i根据公式计算归一化因数 Qi(C, T), 其中 是标准 差, 或者基于第 i 个模数转换器的电压输出 (ADC)i以及第 j 个模数转换器的电压输出 (ADC)i根据公式计算比例因数 Yij(C, T), 其中 是标准差, 其中, T 是生物样品的温度, C 是生物样品中的血糖浓度。 权 利 要 求 书 CN 102988061 B 3 1/11 页 4 用光谱数据分析进行无创的光学血糖检测的方法和系统 0001 本申请是于2009年4月17日提交的。

10、、 题为 “用光谱数据分析进行无创的光学血糖 检测的方法和系统” 的发明申请第 200980126116.7 号的分案申请。 0002 相关申请的交叉引用 0003 本专利申请要求 2008 年 5 月 22 日提交的申请号为 61/055,303 的美国临时专利 申请的优先权, 该美国临时专利申请的公开内容通过引用被合并在此, 本专利申请还要求 2008 年 8 月 15 日提交的申请号为 61/089,152 的美国临时专利申请的优先权, 该美国临时 专利申请的公开内容通过引用也被合并在此。 背景技术 0004 糖尿病是一种慢性疾病, 当不对其进行控制的时候, 随着时间的流逝会对许多身 体。

11、系统 (包括神经、 血管、 眼睛、 肾脏和心脏) 造成严重的损害。美国国家糖尿病、 消化系统 病、 肾疾病研究所 (NIDDK) 估计 2007 年在美国有 2360 万人或者美国总人口的 7.8% 患有糖 尿病。在全世界, 世界卫生组织 (WHO) 估计超过 1.8 亿人患有糖尿病, 到 2030 年他们预计 这个数字将增加到 3.66 亿, 美国占 3030 万。根据 WHO, 在 2005 年估计有 110 万人死于糖尿 病。他们预测 2006 年到 2015 年, 糖尿病致死率将总的增加超过 50% 且在中上收入国家增 加超过 80%。 0005 糖尿病对个人和整个社会造成的经济负担是。

12、重大的。根据美国糖尿病协会, 2007 年美国糖尿病总年度经济成本估计达到 1740 亿美元。自 2002 年以来这个数字增加了 420 亿美元。这 32% 的增长意味着美元数量每年提升超过了 80 亿美元。 0006 糖尿病管理一个很重要的要素就是糖尿病患者在家庭环境下自我监测血糖 (SMBG) 浓度。 通过经常测试血糖水平, 糖尿病患者能更好地管理用药、 饮食和运动以持续控 制和预防长期不良健康结果。实际上, 糖尿病控制与并发症实验 (DCCT) 若干年来密切关注 了 1441 名糖尿病患者, 结果表明, 与标准治疗组比较, 那些每天按照强化控制方案进行多 种血糖测试的患者中只有 1/4 。

13、的人发展为糖尿病性眼疾病, 1/2 的人发展为肾脏疾病, 1/3 的人发展为神经性疾病, 且已经有这三种并发症早期形态的患者中极少数人恶化。 0007 然而, 由于在分析之前通过皮肤抽取血液不方便并且带来痛苦, 这使得许多糖尿 病患者做不到为了很好地控制血糖而应当做的那样勤奋, 因此目前的监测技术阻碍了正常 的使用。结果, 无创测量葡萄糖浓度是糖尿病管理的一种可取和有利的发展。无创监测将 使每天的多次测试没有痛苦且更加受到儿童糖尿病患者的欢迎。根据 2005 年刊登的一篇 研究 (J,Wagner,C.Malchoff, and G.Abbott,Diabetes Technology&The。

14、rapeutics (糖尿病 技术和疗法) ,7(4)2005,612-619) , 患有糖尿病的人利用无创血糖监测装置更频繁地执行 SMBG, 并且已经改善生活质量。 0008 存在许多用于血糖确定的无创方法。 一种无创的血液化学检测技术涉及收集和分 析光谱数据。 0009 由于在正被感测的区域中除了血液之外还存在其它成分 (比如, 皮肤、 脂肪、 肌肉、 骨头、 胞间液) , 所以从根据光谱学获得的光谱数据或其它数据提取关于血液特性的信息 说 明 书 CN 102988061 B 4 2/11 页 5 (例如葡萄糖浓度) 是一个复杂的问题。这样的其它成分能够以改变读数 (reading) 。

15、这样的 方式影响这些信号。具体地讲, 作为结果而获得的信号的幅度可以比信号中与血液对应的 部分的幅度大很多, 因此, 限制了精确提取血液特性信息的能力。 0010 本发明旨在克服上述的一个或多个问题。 发明内容 0011 在本发明的一方面中, 公开了一种用于检测生物样品中的葡萄糖的系统。该系统 包括 : 至少一个光源, 所述至少一个光源被配置为入射在样品的目标区域上 ; 至少一个光 检测器, 所述至少一个光检测器被定位为接收来自所述至少一个光源的光并产生指示被检 测的光的功率的输出信号, 该输出信号具有与时间相关的电流 (time dependent current) ; 和处理器, 该处理器。

16、从所述至少一个光检测器接收输出信号, 基于接收到的输出信号, 以至 少 20 : 1 的信噪比计算由于存在于目标区域中的样品中的血液引起的衰减率, 并基于计算 的衰减率确定与存在于目标区域中的样品相关的血糖水平。 0012 在本发明的另一方面中, 公开了一种用于检测生物样品中的葡萄糖的系统。该系 统包括 : 至少一个光源, 所述至少一个光源被配置为入射在样品的目标区域上 ; 至少一个 光检测器, 所述至少一个光检测器被定位为接收来自所述至少一个光源的光并产生指示被 检测的光的功率的输出信号, 该输出信号具有与时间相关的电流 ; 和处理器, 该处理器从所 述至少一个光检测器接收输出信号, 基于接。

17、收到的输出信号, 以归一化因数计算由于存在 于目标区域中的样品中的血液引起的衰减率, 并基于计算的衰减率确定与存在于目标区域 中的样品相关的血糖水平。 0013 在本发明的另一方面中, 公开了一种用于检测生物样品中的葡萄糖的系统。该系 统包括 : 至少一个光源, 所述至少一个光源被配置为入射在样品的目标区域上 ; 至少一个 光检测器, 所述至少一个光检测器被定位为接收来自所述至少一个光源的光并产生指示被 检测的光的功率的输出信号, 该输出信号具有与时间相关的电流 ; 和处理器, 该处理器从所 述至少一个光检测器接收输出信号, 基于接收到的输出信号, 以比例因数计算由于存在于 目标区域中的样品中。

18、的血液引起的衰减率, 并基于计算的衰减率确定与存在于目标区域中 的样品相关的血糖水平。 0014 在本发明的另一方面中, 公开了一种用于检测生物样品中的葡萄糖的系统。该系 统包括 : 至少一个光源, 所述至少一个光源被配置为入射在样品的目标区域上 ; 至少一个 光检测器, 所述至少一个光检测器被定位为接收来自所述至少一个光源的光并产生指示被 检测的光的功率的输出信号, 该输出信号具有与时间相关的电流 ; 和处理器, 该处理器从所 述至少一个光检测器接收输出信号, 基于接收到的输出信号计算由于存在于目标区域中的 样品中的血液引起的衰减率并消除由所述至少一个光检测器的与温度相关的检测器响应 所引起。

19、的不确定性的影响, 并且, 基于计算的衰减率确定与存在于目标区域中的样品相关 的血糖水平。 0015 在本发明的另一方面中, 公开了一种用于检测生物样品中的葡萄糖的系统。该系 统包括 : 至少一个光源, 所述至少一个光源被配置为入射在样品的目标区域上 ; 至少一个 光检测器, 所述至少一个光检测器包括具有反馈电阻器的前置放大器, 被定位为接收来自 所述至少一个光源的光并产生指示被检测的光的功率的输出信号, 该输出信号具有与时间 说 明 书 CN 102988061 B 5 3/11 页 6 相关的电流 ; 和处理器, 该处理器从所述至少一个光检测器接收输出信号, 基于接收到的输 出信号计算由于。

20、存在于目标区域中的样品中的血液引起的衰减率并消除由所述至少一个 光检测器的与温度相关的检测器响应所引起的不确定性的影响, 并且, 基于计算的衰减率 确定与存在于目标区域中的样品相关的血糖水平。 0016 在本发明的另一方面中, 公开了一种用于检测生物样品中的葡萄糖的系统。该系 统包括具有前置放大器和反馈电阻器的至少一个光检测器。 0017 在本发明的另一方面中, 公开了一种用于检测生物样品中的葡萄糖的方法。该方 法包括 : 利用被配置为入射在样品的目标区域上的至少一个光源 ; 利用被定位为接收来自 所述至少一个光源的光并产生指示被检测的光的功率的输出信号的至少一个光检测器, 该 输出信号具有与。

21、时间相关的电流 ; 用处理器从所述至少一个光检测器接收输出信号, 并基 于接收到的输出信号, 以至少 20 : 1 的信噪比计算由于存在于目标区域中的样品中的血液 引起的衰减率 ; 以及用处理器基于计算的衰减率确定与存在于目标区域中的样品相关的血 糖水平。 0018 在本发明的另一方面中, 公开了一种用于检测生物样品中的葡萄糖的方法。该方 法包括 : 利用被配置为入射在样品的目标区域上的至少一个光源 ; 利用被定位为接收来自 所述至少一个光源的光并产生指示被检测的光的功率的输出信号的至少一个光检测器, 该 输出信号具有与时间相关的电流 ; 用处理器从所述至少一个光检测器接收输出信号, 并基 于。

22、接收到的输出信号, 用处理器以归一化因数计算由于存在于目标区域中的样品中的血液 引起的衰减率 ; 以及用处理器基于计算的衰减率确定与存在于目标区域中的样品相关的血 糖水平。 0019 在本发明的另一方面中, 公开了一种用于检测生物样品中的葡萄糖的方法。该方 法包括 : 利用被配置为入射在样品的目标区域上的至少一个光源 ; 利用被定位为接收来自 所述至少一个光源的光并产生指示被检测的光的功率的输出信号的至少一个光检测器, 该 输出信号具有与时间相关的电流 ; 用处理器从所述至少一个光检测器接收输出信号 ; 用处 理器基于接收到的输出信号, 以比例因数计算由于存在于目标区域中的样品中的血液引起 的。

23、衰减率 ; 以及用处理器基于计算的衰减率确定与存在于目标区域中的样品相关的血糖水 平。 0020 在本发明的另一方面中, 公开了一种用于检测生物样品中的葡萄糖的方法。该方 法包括 : 利用被配置为入射在样品的目标区域上的至少一个光源 ; 利用被定位为接收来自 所述至少一个光源的光并产生指示被检测的光的功率的输出信号的至少一个光检测器, 该 输出信号具有与时间相关的电流 ; 用处理器从所述至少一个光检测器接收输出信号, 并用 处理器基于接收到的输出信号, 以比例因数计算由于存在于目标区域中的样品中的血液引 起的衰减率 ; 用处理器消除由所述至少一个光检测器的与温度相关的检测器响应所引起的 不确定。

24、性的影响 ; 以及用处理器基于计算的衰减率确定与存在于目标区域中的样品相关的 血糖水平。 0021 这些只是本发明的众多方面中的一些, 且不应被看作是关于本发明的众多方面的 全部清单。 附图说明 说 明 书 CN 102988061 B 6 4/11 页 7 0022 为了更好地理解本发明, 可以参考附图, 其中 : 0023 图 1 示出根据示例性实施例的对应于动脉血的光吸收的脉冲波图 ; 0024 图 2 示出用于获得光谱数据的示例性系统 ; 0025 图 3 示出用图 1 中的数据根据公式 (9) 计算的 A(t) 的图 ; 以及 0026 图 4 是在数字化之前将光电流转换为电压的前置。

25、放大器电路的基本示意图。 具体实施方式 0027 在下面的详细描述中, 为了提供对本发明的全面的理解, 阐述了大量的示例性的 具体细节。 然而, 本领域技术人员将会理解, 本发明可以在没有这些具体细节或者对所述细 节的各种修改的情况下实施。在其他示例中, 没有详细描述众所周知的方法、 工艺和部件, 以避免使本发明模糊。 0028 可用光谱学确定被生物样品 (例如人的手指) 所吸收和散射 (即衰减) 的光量。通 过测量被样品吸收的光量, 可无创地确定被检者的葡萄糖、 胆固醇和血红素水平。 由于指尖 中的毛细血管的密集度大并且在指尖发生从动脉血到静脉血的转换, 所以通常优选指尖测 量。然而, 本发。

26、明的方法并不限于用指尖。例如, 生物样品可以是人的耳垂。 0029 当光透射穿过生物样品 (例如人的手指) 时, 光被手指的各种成分 (包括皮肤、 肌 肉、 骨头、 脂肪、 胞间液和血液) 所衰减。然而, 观察到, 通过人的手指的光衰减表现出与心跳 对应的小范围的循环模式。 人们相信, 这种循环模式将存在于许多其他的身体部位 (耳垂是 众多例子中的一个) 的测量中。 0030 图1描绘了检测器光电流ID(t)的图102, 该光电流对应于在光穿过被检者的手指 后被检测器接收的光功率。可以看出, 检测器光电流表现出循环模式。该循环模式是由于 被检者的心跳而引起的, 且这种心跳循环地增加并减少被检者。

27、的毛细血管 (或其他结构) 中 的血量。尽管该循环模式的大小相比于检测器产生的总光电流而言小, 但是可以从图 102 的循环模式中提取出大量的信息。例如, 假定人的心率是 60 次每分钟, 即, 任何脉搏的开始 和结束之间的时间是一秒钟。在这一秒钟的期间内, 光电流将具有最大值或峰值读数 104 和最小值或谷值读数106。 图的峰值读数104对应于当毛细血管中有最小量的血时, 谷值读 数 106 对应于当毛细血管中有最大量的血时。通过使用所述循环图中的峰值和谷值提供的 信息, 排除不在毛细血管中的主要的手指成分 (例如皮肤、 脂肪、 骨头、 肌肉、 胞间液) 所造成 的光吸收和散射。 排除这些。

28、不在毛细血管中的主要成分是由于其在一个心跳的时间间隔期 间内不可能变化。换句话说, 可基于图 102 的峰值和谷值检测被血液吸收和散射 (即, 衰减) 的光。 0031 假定光感测装置产生的循环光电流的峰值是 IP, 循环光电流的相邻的谷值是 Iv, 在没有人手指的情况下光感测装置产生的光电流是 I0, 对应于峰值和谷值光电流的透射率 可定义为 : 0032 0033 和 说 明 书 CN 102988061 B 7 5/11 页 8 0034 0035 对应的峰值和谷值吸光率为 : 0036 Av -log(Tv) (3) ; 0037 和 0038 AP -log(TP) (4)。 003。

29、9 Av和 AP之间的差表示手指中的血液 (不包括非血液成分) 导致的光吸收和散射 : 0040 0041 在公式 (5) 示出的算法中可以看出, A 不依赖于 I0。所以, 计算 A 不需要确定 在无样品时由光感测装置产生的电流。监测对应于透过样品的光功率的光电流足以计算 A。 0042 图 2 描绘了用在示例性实施例中的示例性设备的简化框图。总体上用数字 200 表 示的光学测量系统用 “脉动” 的概念确定仅被样品 (在本示例性实施例中是人的手指) 中的 血液吸收和散射的光量。电源 201 (例如电池) 为光源 202 提供电力, 该光源 202 产生多条向 着被检者的手指的顶部引导的光束。

30、 204、 206、 208、 210。在一个示例性实施例中, 光束 204、 206、 208、 210 中的每一条具有相同的波长或不同的波长范围, 典型地在 800nm 到 1600nm 之 间。尽管此处将光学测量系统 200 描述为产生四 (4) 条光束, 但是在其它实施例中, 设想可 将光源 202 改装成产生更少的光束或者另外的光束。 0043 第一光圈212确保光束204、 206、 208、 210入射到手指的目标区域上。 第二光圈214 确保透过手指的光束的一部分入射到透镜 216 上。光束 204、 206、 208、 210 被手指和光学 测量系统 200 的各部件衰减, 。

31、所以衰减光束 218、 220、 222、 224 从手指发射。衰减光束 218、 220、 222、 224入射到透镜216上, 透镜216收集衰减光束218、 220、 222、 224以使其更有效地 入射到检测器块 226 上。 0044 检测器块 226 位于透镜 216 的正下方, 并且包括多个光感测装置 (LSD) 228、 230、 232、 234, 例如光电二极管阵列。根据光学测量系统 200 的一个方面, 光感测装置 228、 230、 232、 234 的每个分别检测由对应的干涉滤光器 (IF) 236、 238、 240、 242 限定的具体波长的 光。干涉滤光器透射一。

32、个或多个光谱带或光线, 并阻挡其它光谱带或光线。 0045 光感测装置 228、 230、 232、 234 的每个产生与被特定光感测装置接收的光的功率 成比例的对应的光电流信号。可将光电二极管所产生的光电流信号转换为另一形式的信 号, 例如模拟电压信号或数字信号。处理器 243 与检测器块 226 耦合, 并被构造为计算光电 流信号 244、 246、 248、 250 的变化。 0046 根据一方面, 处理器 243 执行例如公式 (5) 中显示的算法以计算仅由手指中的血 液引起的光吸收变化 (A) 。 其后, 可使用对血液的光吸收的这个量化计算来确定血液的特 性。例如, 通过将计算的光吸。

33、收值同与存储器 (未显示) 中存储的不同葡萄糖水平对应的预 定值进行比较, 可确定被检者的血糖水平。 0047 与基于手指的脉动检测方法相关的困难在于低的信噪 (S/N) 比, 因为循环模式的 振幅 (即, 峰值和谷值之间的差) 通常是由透过手指的光功率产生的总光电流的 1%-2%。为 说 明 书 CN 102988061 B 8 6/11 页 9 了在确定 A 时获得 100:1 的 S/N 比, 在 10Hz 的带宽内, 在吸光率 (峰值到峰值) 中用于测 量手指造成的光吸收的装置的基线噪声不应大于 3.010-5。 0048 然而, 以一些电池供电的手持无创的血液化学测量装置使用的低的光。

34、功率水平, 难于获得在 10Hz 的带宽内 3.010-5吸光率 (峰值到峰值) 基线噪声水平。一种解决方案涉 及数据取平均。为了提高 S/N 比, 下式定义的 A 的平均值用于进一步的计算以提取血糖 浓度 : 0049 0050 在公式 (6) 中, M 是脉动测量的时间间隔期间的心跳的次数。然而, 由于心率在每 秒钟一次的数量级上的事实而导致该方法需要较长的数据采集时间。例如, 为了使 S/N 比 提高到五倍, 将需要 25 秒钟, 为了使 S/N 比提高到十倍, 将需要 100 秒钟。相较而言, 现在 的商用的抽血血糖仪能在 5 秒钟内确定血糖水平。另外, 较长的检测时间会显著增加由于 。

35、手指移动、 光功率漂移、 装置温度变化等引起的测量误差。所以, 需要新的技术来快速而准 确地测量血糖水平。 0051 通过标准差提高 S/N 比 0052 图 1 中所示的与时间相关的检测器光电流输出 ID(t) 可表示为与心跳相对应的小 的与时间相关的循环光电流 I(t)、 噪声电流 n(t) 和恒定的基线光电流 IB的和 : 0053 ID(t) IB+I(t)+n(t) (7)。 0054 上式可被重新整理为 : 0055 0056 将常用的对数应用于公式 (8) 的两边, 得到 : 0057 0058 总体上用数字 300 表示的图 3 示出了使用图 1 中的数据根据公式 (9) 计算。

36、的典型 的 A(t) 的图 302。对于图 3 中示出的脉冲函数 A(t), 在一个心跳的时间间隔期间存在下面 的关键的关系 : 0059 A(t) kA (10), 0060 其中 A(t) 是 A(t) 的标准差, k 是比例常数。 0061 考虑到 IB是常数, 且 2(logIB) 0, 得出 : 0062 A(t) logID(t) (11)。 0063 所以, 在一个心跳的时间间隔期间 A(t) 图的峰值到谷值的高度可直接从 ID(t) 的 对数的标准差中得到 : 0064 0065 公式 (12) 的一个主要的好处是如下所解释的可在较短的数据采集时间 (约为一 秒) 内得到高的 。

37、S/N 比。 说 明 书 CN 102988061 B 9 7/11 页 10 0066 在图 2 所描绘的基于手指的脉动测量中, I(t)+n(t) 的和值通常小于大的恒定 的基线光电流 IB的 2%。所以, 公式 (9) 可近似为 : 0067 0068 类似地, A(t) 的标准差可近似为 : 0069 0070 公 式 (14)说 明 了 公 式 (12)的 大 的 噪 声 降 低 功 率。 例 如, 对 于 比 值 为 (或 10%) 的相对高的基线噪声, 基线噪声 n(t) 对于 A(t) 的贡献估计 少于 0.005(或 0.5%) , 其对应于在没有增加检测时间的情况下 S/N。

38、 比增加到 20 倍。照此, 可在不增加数据采集时间的情况下得到惊人的噪声降低, 可在一个心跳 (近似为一秒) 的时 间间隔内完成基于手指的脉动测量, 且对于在确定 A 时 100:1 的 S/N 比的要求可以使用 在 10Hz 的带宽内具有大约 6.010-4吸光率 (峰值到峰值) 的基线噪声的光学系统得到满 足。 应当指出的是, 当由于低的光照功率使得光学系统的基线噪声由散粒噪声主导时, 噪声 降低到 20 倍等同于光照功率增加到 202=400 倍。 0071 在很短的数据采集时间 (例如少于一秒) 内获得较高的S/N比的这种能力将显著地 减少由于测量期间的诸如手指移动、 温度变化和光功。

39、率漂移的因素引起的检测误差, 所以 惊人地提高脉动检测方法的准确性和可再现性。 0072 另外, 因为所有波长的传输光由于心跳而具有相同的脉冲形状, 因此 k 的值不会 随波长改变。 结果是, 常数k将在下面部分讨论的数据归一化中被消除, 且logID(t)将 被用于进一步的回归分析以在光学测量和血糖水平之间建立关联。 因为logID(t)仅包 含两个在大多数普及的电子数据表程序 (例如微软) 中可用的标准数学函数, 所以这将大大简化了数据分析过程。是微软公司的联邦注册商标, 微软公司的 办公地点在 One Microsoft Way,Redmond,Washington 98052-6399。

40、。 0073 归一化 0074 在每个波长 i处, 通过下面的公式将吸收 A(i) 与由于心跳导致的指尖的光 学感测区域中的血液量的增加 (B) 相联系 : 0075 A(i) (C, i, T)B (15), 0076 其中 (C,i,T) 是血液在波长为 i、 手指温度为 T 且血糖浓度为 C 时的吸收 / 散射系数。众所周知, 变量 B 因人而异, 且甚至对于同一个人可能每天也会不同。 0077 变量 B 导致的不确定性可通过在每个波长 i处引入如下面的公式所定义的归 一化因数 Qi(C,T) 而被消除 : 0078 0079 其中 N 是使用的波长的总数。优选地, N 的范围通常为 2。

41、0 到 30。 0080 基于公式 (12) 和 (16) , 通过下面的公式将 Qi(C,T) 与每个波长 i处的检测器光 说 明 书 CN 102988061 B 10 8/11 页 11 电流 ID(i,t) 相联系 : 0081 0082 如公式 (17)所示, 常数 k 被消除, 且将 logID(t) 用于进一步的回归分析以 在光学测量和血糖水平之间建立关联。这可能是因为同时从所有的检测通道 (detection channel) 获取数据。 0083 光学测量和血糖浓度之间的关联可根据下式建立 : 0084 0085 其中 Coptical是通过光学测量预测的血糖浓度, Qi(C。

42、,T) 由公式 (16) 和 (17) 定义, 且 ai(T) 是对应于波长 i的与温度相关的回归系数。ai(T) 的值可用适当的统计方法 (例 如, 偏最小二乘 (PLS) 回归) 提取。 0086 公式 (18) 表示可在不同的手指温度进行大量的校准时的理想情况。实际上, 通常 只能进行有限数量 (例如, 15 到 20) 的校准, 且每次校准可能在不同的手指温度进行。在这 种条件下, 手指温度可被视为独立的变量, 且上面的公式可近似为 : 0087 0088 其中 bi是对应于波长 i的与温度无关的回归系数, 是手指温度的回归系数。 bi和 的值可用适当的统计方法 (例如, 偏最小二乘 。

43、(PLS) 回归) 提取。 0089 比例方法 0090 可替代地, 来自变量 B 的不确定性可通过在波长 i处引入比例因数 Yij而被消 除 : 0091 0092 其中 j 可为从 1 到 N 的任意数, 假定该装置在所有 N 个波长处收集信号。 0093 和前面讨论的归一化算法类似, 光学测量和血糖水平之间的关联可根据下式建 立 : 0094 0095 其中 Coptical是通过光学测量预测的血糖浓度, Yij(C,T) 由公式 (20) 定义, 且 fi(T) 是对应于波长 i的与温度相关的的回归系数。fi(T) 的值可用适当的统计方法 (例如, 偏 最小二乘 (PLS) 回归) 提。

44、取。 0096 公式 (21) 表示可在不同的手指温度进行大量的校准时的理想情况。实际上, 通常 只能进行有限数量 (例如, 15 到 20) 的校准, 且每次校准可能在不同的手指温度进行。在这 种条件下, 手指温度可被视为独立的变量, 且上面的公式可近似为 : 说 明 书 CN 102988061 B 11 9/11 页 12 0097 0098 其中 hi是对应于波长 i的与温度无关的回归系数, 是手指温度的回归系数。 hi和 的值可用适当的统计方法 (例如, 偏最小二乘 (PLS) 回归) 提取。 0099 与温度相关的装置响应的影响的消除 0100 众所周知, 硅光电二极管检测器的检测。

45、器灵敏度是波长和温度的函数。对于图 2 中所示的总体上由数字 200 表示的装置配置, 由第 i 个硅二极管检测器接收的对应于波长 i的光功率根据下式被转换为光电流 : 0101 ID(i, t) P(i, t)S0(i)1+(i)(TDi(t)-25 C) (23)。 0102 上面的公式 (23) 中, P(i, t) 是由检测器接收的光功率, S0(i) 是在波长 i和 25处的检测器的光灵敏度, (i)是波长i处的光灵敏度的温度系数, TDi(t)是第i个 光电二极管检测器的温度。温度系数 (i) 随波长变化。例如, 对于 Hamamatsu 的 S1337 系列光电二极管检测器, (。

46、i) 的范围从 900nm 处接近零到 1100nm 处超过 1.0%/。这 给图 2 所示的装置配置施加了一个潜在的问题, 因为很难在正常的家庭 / 办公环境下在糖 尿病人使用的手持设备中保持每个单独的二极管检测器的温度恒定不变。 0103 由于检测器温度 TDi(t) 引起的这一不确定性可用公式 (11) 和 (12) 所示的算法来 被消除。将常用的对数应用到公式 (23) 的两边, 得到 : 0104 0105 考虑到 S0(i) 是常数且检测器温度 TDi(t) 在约为一秒的很短的数据采集时间间 隔期间几乎保持恒定不变的事实, 得到 : 0106 logID(i, t) logP(i,。

47、 t) (25)。 0107 如此, 检测器温度 TDi(t) 所导致的不确定性可通过使用该标准差方法来被消除。 0108 电压检测模式 0109 在图2所示的装置配置中, 通常, 在数字化之前用前置放大器将第i个光电二极管 检测器的光电流 ID(i,t) 转换为电压。图 4 示出总体上用数字 400 表示的典型的前置放 大器的电路示意图。 0110 与第 i 个光电二极管检测器 408 耦合的第 i 个前置放大器 400 的输出电压 412 可 表达为 : 0111 Vi(t) RiID(i, t) R0i1+i(TRi(t)-25 C)ID(i, t) (26)。 0112 在上面的公式 。

48、(26) 中, R0i是在 25处第 i 个前置放大器的反馈电阻器 402 的电 阻值, i是电阻器的温度系数, TRi(t) 是电阻器的温度。将常用的对数应用到公式 (26) 的 两边, 得到 : 0113 logVi(t) logRi0+log1+i(TRi(t)-25 C)+logID(i, t) (27)。 0114 考虑到 R0i是常数且电阻器温度 TRi(t) 在约为一秒的很短的数据采集时间间隔期 间不变的事实, 得到 : 0115 logVi(t) logID(i, t) (28)。 0116 将公式 (25) 代入公式 (28) , 得到 : 说 明 书 CN 102988061 B 12 10/11 页 13 0117 logVi(t) logP(i, t) (29)。 0118 如此, 电阻器温度 TR(t) 所导致的不确定性被消除。 0119 在电压检测模式下, 公式 (17) 中的归一化因数可表达为 : 0120 0121 然后, 在相应的校准条件下, 光学测量和血糖浓度之间的数学关联可根据公式 。

展开阅读全文
相关资源
猜你喜欢
相关搜索

当前位置:首页 > 人类生活必需 > 医学或兽医学;卫生学


copyright@ 2017-2020 zhuanlichaxun.net网站版权所有
经营许可证编号:粤ICP备2021068784号-1