基于有效衰减系数直接测量的人体光谱特征分析

刘 瑾， 付润娟， 韩同帅， 刘 蓉， 孙 迪

天津大学精密测试技术及仪器国家重点实验室， 天津 300072

引 言

近红外光谱具有无损、 高效、 实时等优点， 在人体成分的无创测量领域有广泛的应用[1]。 成功的应用有血氧饱和度[2]、 血红蛋白浓度[3]的测量， 均有成熟的仪器产品。 另外， 基于近红外光谱法的无创血糖测量是受到广泛关注、 被寄予厚望的一个重要应用领域[4]， 其难度较大。

在活体成分测量中， 人体组织的光谱间存在千差万别， 即使同一组织的光谱也时刻发生着变化， 这是影响活体组织成分测量精度的重要原因。 组织的漫反射光谱并不随光源-探测器距离(source detector separation, SDS)呈现线性变化， 在不同距离下的漫射光谱(或不同厚度下的漫透射光谱)， 不易相互转换和比较。 不同距离下的漫射光信息中所包含的吸收、 散射的比例均不同[5]， 在1 000～1 600 nm波段， 甚至出现对组织的散射差异不敏感、 只对吸收敏感的测量距离[5-7]。 因此， 若采用不同的测量距离进行测量， 其包含的信息量是不一致的。 因此， 非常有必要找到一个与测量距离无关、 有明确物理意义的光学参数来表征组织光谱的光学特性， 以此来考察活体组织的差异或者其变动情况。

利用实时测量的漫射光谱直接反推组织的基本光学参数有一定的难度。 常用的方法是反向倍-增法[8]， 它采用漫透射、 漫反射相结合的方式， 更适合离体组织的测量。 还有一些依靠智能算法， 如神经网络等， 进行反推的数学方法。 近年来， 随着多距离漫射光谱测量技术的发展， 学者们基于多个距离下丰富的漫射光信息， 直接获得了组织的一个综合光学参数， 即有效衰减系数[9-13]。 Chiarelli等对漫反射光强进行校正， 将其校正为随SDS线性变化的新变量， 并将其应用于脑组织活动状态的监测[9-10]； 刘瑾等采用双SDS差分测量， 详细推导了无限介质、 半无限介质下的差分测量公式(含有效衰减系数)[11]； 曹海青等利用双SDS测量了多种Intralipid溶液的有效衰减系数[12]； 韩同帅等利用双SDS差分法获得的有效衰减系数进行了血糖测量[13]。 上述研究均证实了差分测量获得的有效衰减系数是不依赖测量距离的光学参数。

本文将基于双/多SDSs的差分测量用于人体组织光谱有效衰减系数的测量， 考察人体组织的光谱特征。 设计仿体实验， 测试了葡萄糖、 血红蛋白、 颗粒密度、 温度四个因素影响下介质的有效衰减系数光谱， 比较了它们的光谱特征。 进一步， 进行了人体实验， 考察了人体组织光谱的差异。

1 实验部分

1.1 基于双光源-探测器距离差分的有效衰减系数估计

根据漫射方程的一阶近似解， 在半无限介质的组织表面、 SDS为ρ的距离下的漫反射光能流密度可表示为[14]

(1)

(2)

两个距离下(rA和rB)的漫反射光(IrA和IrB)的吸光度进行差分， 得到新的变量AD

(3)

将式(2)代入式(3)可得

AD=μeff(rA-rB)+Ksinf

(4)

(5)

上述推导为理想情况下， 且做了一定的简化， 这不影响式(5)的使用。 获得的有效衰减系数可能与真实值存在系统误差， 但其是一个与测量距离无关的、 稳定反映组织光学特性的量， 可看作有效衰减系数的估计值， 并可进一步通过数学方法进行系统误差的修正。

1.2 方法

采用双距离的漫反射测量系统进行仿体溶液和人体组织的测量。 图1是实验系统的示意图。 实验系统包括卤素灯、 光谱仪(爱万提斯， 型号： AvaSpec-NIR256/512-1.7(TEC)， 波长范围900～1 750 nm)， 光开关(中电科34所)、 光纤等。 测头包括定制的三环结构光纤， 中心环为光源入射， 另外两个为漫反射光接收环。 接收环的内外缘中心距入射光纤中心的距离即可视为光源-探测器距离， 第一个环(内环)对应的SDS约0.6 mm， 第二个环(外环)对应的SDS约2.0 mm。 利用光开关对两路光信号分时采集。

图1 双光源-探测器距离的漫反射光谱测量系统示意图

实验采用3%和3.5%的Intralipid溶液作为仿体溶液， 并分别对两种仿体溶液配置了葡萄糖浓度分别为0 mmol·L-1， 55.55 mmol·L-1(1 000 mg·dL-1)， 111.1 mmol·L-1(2 000 mg·dL-1)， 166.65 mmol·L-1(3 000 mg·dL-1)和222.2 mmol·L-1(4 000 mg·dL-1)的样品、 血红蛋白浓度分别为2和4 g·L-1的样品， 共14个样品溶液。 对上述14个样品在室温(约25 ℃)下测量漫反射光强。 采用恒温水浴箱对3%的Intralipid溶液样品进行加热和保温， 采用TES温度传感器(TES1310， 分辨率0.1 ℃)实时测试样品的温度， 分别测量溶液在27和29 ℃下的漫反射光强。

对每次测量获得的两个SDSs(0.6和2.0 mm)下的漫反射光强代入式(3)和式(5)， 计算得到有效衰减系数。

1.3 人体实验

对三名受试者(22岁女， 23岁女， 25岁男)的五个部位(手指、 手掌、 手背、 手臂外侧、 手臂内侧)进行了有效衰减系数光谱的测量。

2 结果与讨论

2.1 葡萄糖光谱测量结果

图2给出了3%intralipid溶液的有效衰减系数光谱和溶液中四种葡萄糖浓度引起的有效衰减系数变化量的光谱。 在1 000～1 300 nm波段葡萄糖吸收非常小， 因此该波段的光谱变化主要是由葡萄糖引起的介质折射率变化引起的。 1 mmol·L-1的葡萄糖浓度变化大约引起2.5×10-5的溶液折射率增加[15-16]。 溶液折射率的增加将导致介质散射系数的减小， 引起有效衰减系数的减小。

图2 3%Intralipid溶液中的不同浓度葡萄糖的测量结果

2.2 散射颗粒密度、温度、血红蛋白的光谱测量结果

3% intralipid溶液和3.5% intralipid溶液的差别可以看作是散射颗粒密度的差异。 图3(a)给出了几种溶液样品的有效衰减系数光谱， 以及在3%intralipid溶液中加入不同血红蛋白浓度、 以及不同温度下的有效衰减系数光谱变化量结果。 由图3(b)可以看出颗粒密度、 温度、 血红蛋白引起的有效衰减系数光谱， 其中血红蛋白的光谱与其吸收光谱相似， 温度光谱与温度引起的水吸收光谱形状相似， 而散射系数变化引起的光谱与有效衰减系数本身的形状一致， 因为散射系数本身随波长的变化没有明显的形状。 可见， 有效衰减系数光谱可以基本反映吸收光谱的影响， 但在血红蛋白、 温度的光谱中也包含了散射系数的光谱， 若想提取纯吸收光谱， 还需要继续开发散射校正的算法。 另外， 比较葡萄糖光谱、 散射系数变化的光谱可知， 两者都是由散射变化引起的， 不易区分。

图3 几种因素变化时的3%Intralipid溶液有效衰减系数光谱

2.3 人体测量结果

人体皮肤自身的随机变异性对测量造成影响。 对同一部位(某受试者手指和手臂内侧)进行了多次重复测量， 结果如图4所示。 在1 000～1 300 nm波段， 有效衰减系数的变化约为±0.5 cm-1， 产生的相对测量误差约为10%以内。 其中手指处微血管丰富， 脉搏显著， 其光谱稳定性稍差。

图4 某受试者手指和手臂内侧的有效衰减系数光谱

不同部位、 不同个体的测量结果如图5所示， 可以明显看出不同皮肤有效衰减系数的差异。 图5(a)中的五个部位， 即手掌、 手指、 手背、 手臂外侧、 手臂内侧之间均存在显著差异， 图5(b)给出了手掌、 手背、 手臂外侧、 手臂内侧与手指的光谱差异， 这些差异是组织吸收(水、 血红蛋白等)、 散射差异的共同作用结果。 图5(c)给出了三个人手指的有效衰减系数， 图5(d)给出了它们的光谱差异， 在血红蛋白吸收的1 000～1 200 nm波段和散射较明显的1 200～1 300 nm波段， 都能看到光谱差异。

图5 不同部位、 不同个体的人体有效衰减系数光谱测量结果

3 结 论

采用多光源-探测器距离的差分测量， 给出了与测量距离无关的组织的有效衰减系数光谱测量方法， 为组织光谱的测量提出了新的思路。

在1 000～1 300 nm波段评估了葡萄糖、 血红蛋白、 颗粒密度、 温度四个因素有效衰减系数光谱的特征。 其中葡萄糖、 散射颗粒密度的光谱特征接近。 可见， 在该波段对于血糖测量而言， 组织的颗粒密度构成了葡萄糖测量的干扰因素。 尤其是活体测量场合下， 如何保持组织密度的一致性给测量接口提出了严峻的考验。 另外， 血红蛋白浓度和温度变化均具有显著的特征， 与它们的吸收光谱形状相似， 是与葡萄糖光谱不同的， 可以通过多变量分析进行消除。 不同部位、 不同人体之间的光学参数差异， 也有望通过组织的有效衰减系数光谱进行区分。

基于有效衰减系数光谱还可以继续开发吸收系数、 散射系数的反演算法， 或者配合散射校正技术获取到组织的纯吸收光谱， 为精准的组织光谱分析提供有力的工具。

综上， 本文提出的有效衰减系数光谱测量方法有望成为一种活体生物组织的通用测量方法， 在生物组织光谱测量领域将有广泛的应用前景。