1. 引言
糖尿病是一种由于胰腺无法正常生产胰岛素,从而导致细胞无法从血液中有效吸收葡萄糖的代谢紊乱[1],作为全球性的主要健康问题,其治疗和管理在20世纪初迎来了革命性的变化。目前,糖尿病已成为全球性的健康问题,世界卫生组织估计全球有超过4.63亿人患有糖尿病,且这一数字还将继续上升。
传统的血糖监测方法例如自我血糖监测(Self-Monitoring of Blood Glucose, SMBG) [2]、连续血糖监测(Continuous Glucose Monitoring, CGM)系统[3] [4]等,都存在着需要侵入式采血或者成本较高的问题。为了克服现有方法的局限性,研究人员大力推动了无创血糖监测技术[5]的发展,探索了包括光学/光谱学[5] [6]、超声波和微波传感[7] [8]在内的多种技术。这些方法可在无需采血的情况下测量血糖水平,为患者提供了更为友好的监测方式。
射频/微波传感技术已被证明在非侵入式血糖监测方面更有前景,[7] [8]中探索了天线或者波导的多种反射和传输方法应用于血糖的检测,这些方法依赖于矢量网络分析仪(Vector Network Analyzer, VNA) [9],以表征不同设置下血样的电磁特性(即介电常数和损耗特性),值得注意的是,所测量的血液特性与其葡萄糖含量之间存在明显的相关性。然而,考虑到成本和使用的复杂性,大多数射频系统并不支持方便、快捷的血糖检测。研究人员还在探索毫米(mm)波段的传感结构,以实现更短的波长和相对可接受的穿透深度,达到检测身体薄组织区域葡萄糖浓度的目的。
本文介绍了一种改进型便携式平面微波传感器设计,该传感器利用圆形互补开口环谐振器(Complementary Split-Ring Resonator, CSRR)提高无创血糖检测的灵敏度。CSRR在感应共振区域集中电磁场,增强传感器对葡萄糖分子引起的介电特性变化的敏感性。传感器配置采用中心对称的菱形排列CSRR单元,蚀刻在FR4介质基板的接地金属层上,通过基板顶部的平面微带线耦合馈电。这种设计不仅提高了血糖监测的灵敏度,而且确保了便携性和易用性,对需要长期血糖监测的糖尿病患者非常重要。本文重点研究了谐振频率对检测灵敏度的影响,仿真和数据分析结果表明,谐振频率与检测灵敏度之间存在线性关系,为精确测量血糖水平提供了设计参考。
2. 传感器设计
2.1. 葡萄糖浓度检测原理
CSRR结构是本传感器检测葡萄糖浓度的核心所在,下面使用简单的集总元件模型来描述集成CSRRs的工作机制[10]-[12]。该模型由三部分组成[10],如图1所示。
Figure 1. Lumped element model
图1. 集总元件模型
第一部分(红框)是用于激励谐振器的传输线(MTL, Microstrip Line),通过电感L1、L2建模。第二部分(黄框)是介质基板,用于将MTL与接地金属板上的谐振器进行电耦合。第三部分(绿框)是四个圆形谐振器单元,由RCL并联谐振电路建模,四个圆形单元中的介电环(宽度)绕组充当电感L3,介电环开口间隙和双环间距产生并联电容C2,传导损耗和介电损耗由电阻R2建模。第四部分(蓝框)为待测样品。当葡萄糖溶液放置在CSRRs表面的顶部时,模型中通过添加一个并联RC电路来调整,C3与样品的相对介电常数直接相关,而R3主要取决于其损耗特性。此部分在检测葡萄糖样本浓度方面起关键作用,通过测量CSRRs与葡萄糖溶液相互作用下的电磁特性变化,实现高灵敏度的血糖检测。
2.2. 传感器建模及仿真
传感器模型如图2所示。
Figure 2. Proposed sensor model
图2. 传感器模型
Figure 3. CSRR structure
图3. CSRR结构
该传感器将四个相同的圆形互补开口环谐振器(CSRR)单元集成在接地金属板上,与输入输出端口、传输线以及介质基板共同构成主体结构。其中,介质基板为FR-4,介电常数ε = 4.3,正切损耗tanδ = 0.02,整体尺寸为66 × 24 × 1.5 mm3。无源CSRR的磁壁垂直于基板上表面上的MTL,确保了MTL与接地平面之间随时变化的电场可以谐振激励CSRR。将MTL尺寸优化为66 × 2.8 mm2,以获得与输入端口匹配的50欧姆特性阻抗。如图3所示,CSRRs结构位于接地金属板中心,四个单元呈中心对称分布的菱形排列。横向两个圆形单元间距CL = 14 mm,竖向两个圆形单元间距CW = 12.6 mm。每个CSRR的单元由两个相互嵌套的圆形同心环组成,外环直径a = 7.6 mm,其间的耦合间隙t = 0.4 mm。每个环的环路都有一个边宽s = 0.4 mm的介电槽,末端有一个宽度g = 0.4 mm的金属狭缝,两个环的开口间隙位于彼此直径相对的两侧。
考虑到实测环境中使用的VNA在反射系数(幅度和相位)读数方面具有更高的不确定性。因此,本文仅考虑传输系数S21的频移来评估和比较所提出的传感器的灵敏度范围。图4为模型空载时的S21曲线。从S21曲线可知,模型的谐振点为2.88 GHz,具有较深且陡峭的谐振曲线,意味着拥有较大的Q值,在涉及频偏的检测上存在较大的潜力,可以作为血糖检测模型进行进一步测试。
由于外环直径a与介电槽s的长度成正比,因此增加a将直接增加所得的电感,谐振点也会因此向低频移动,如图5所示。
将葡萄糖样品加载到CSRRs表面,如图6所示。为达到最佳测试效果,样品需完全覆盖CSRRs结构,以a = 7.6 mm为例,此时取样本溶液半径为11 mm。样品与接地金属层之间留有0.15 mm空隙,以防止二者接触导致的短路,同时也为后续加载玻璃容器做准备。
本实验采用0~500 mg/dL的葡萄糖溶液模型代替血浆模型,该浓度范围涵盖了低血糖(<70 mg/dL)和高血糖(>140 mg/dL)。许多其他研究[13]-[15]也表明了这种近似是有效的。式(1)给出了一阶Debye模型[10] [16],用于建立不同浓度葡萄糖溶液的介电特性数值模型[17]。目前,该模型被认为是近似血浆样本最合理的模型,尚未有其他针对实际血糖溶液开发的数学模型被提出。
Figure 4. Simulation result: S21 curve
图4. 仿真结果:S21曲线
Figure 5. Variation of resonance point with parameter a
图5. 谐振点随参数a的变化
Figure 6. A model for loading glucose solution samples
图6. 装在葡萄糖溶液样本的模型
(1)
其中,
是葡萄糖浓度
(单位为mg/dL)的水溶液在角频率w下的复介电常数,
、
、
是与浓度相关的德拜系数[18],
是静态电导率,
是自由空间介电常数。如图7所示,加载的葡萄糖溶液模型干扰了介电基板敏感区域中产生的电场,导致谐振器周围介质的有效介电常数改变,从而引起谐振频率向高频移动。
Figure 7. S21 curve when unloaded/loaded
图7. 空载/负载时的S21曲线
Figure 8. S21 curve of glucose samples with different concentrations (
(70~80))
图8. 不同浓度葡萄糖样品(
(70~80))的S21曲线
图8显示了模拟传感器对同体积不同浓度葡萄糖样品(
(70~80))的S21。在4.6~4.9 GHz内,传感器产生了明显的谐振频率,并且随葡萄糖浓度产生偏移。同时,共振幅度也产生了显著变化,这意味着样本的损耗特性与传感器的频率响应也有相关性。此外,图8还表明了传感器在整个预期范围内检测葡萄糖的能力。Debye模型显示[16],葡萄糖溶液在每相差100 mg/dl时,介电常数变化为
,以此为该模型的灵敏度参考,其值约为5.06 MHz/(100 mg/dl)。
3. 模型优化与数据处理
3.1. 模型优化
考虑到一般传感器在高频时更易受到外部环境的干扰与影响,因此,如果增加传感器的工作频率,将有可能提升葡萄糖浓度的检测敏感度。由2.2节中的图5可知,CSRRs的外环直径a对该模型的谐振点起主要的影响作用,通过修改a的取值,我们可以调整模型检测的频段。同时,介质基板与样本的尺寸也要做相应调整,以实现检测频点更好的谐振深度。
具体测试结果S21如图9所示。图9(a)~(d)分别表示外环直径a = 4.6、5.6、6.6、7.6 (mm)时,传感器负载葡萄糖溶液模型(
(70~80))后,谐振点随样本介电常数的变化。以2.2节中提到的灵敏度参考标准可以计算出对应的灵敏度分别为9.00、6.46、5.60、5.06 (MHz/(100 mg/dL))。随着传感器谐振点调整至更高频率,其对葡萄糖溶液样本的灵敏度也在不断提高,可见该方案对于改善此类型传感器的检测灵敏度具有一定的作用。
3.2. 数据处理
为了更好地理解数据的变化趋势,量化理论预测与实验数据之间的差异,本节将对仿真得到的数据进行拟合分析,以便得到该传感器在用于无创血糖检测方面的通用模型。具体的数据处理在MATLAB
Figure 9. S21 results of sensors on glucose samples of different concentrations (
(70~80)) at various detection frequencies
图9. 各检测频点下传感器对不同浓度葡萄糖样品(
(70~80))的S21结果
中进行。图10为外环直径a和频偏的多项式拟合结果,其中纵轴为传感器检测灵敏度的频偏,横轴为a的取值,图10(a)~(d)分别为一次至四次的拟合结果。从拟合结果来看,传感器的检测灵敏度与CSRR的尺寸成线性相关,随着整体尺寸的缩小,传感器对于葡萄糖溶液介电常数的变化更为灵敏。此外,将传感器空载时的谐振点与检测的频偏进行拟合得到图11,很好地体现了两者之间的线性关系,其中一次拟合的R2就超过了0.98。
4. 结论
本文提出了一种低成本、高灵敏度的微波传感器,用于实现无创血糖检测。传感器部分采用微带结构实现,在接地金属板上蚀刻CSRR单元。实验结果表明,所提出的传感器在0~20 GHz内表现出了不错
Figure 10. Fitting results of parameter a and frequency offset
图10. 参数a与频偏的拟合结果
Figure 11. Fitting results of sensor resonance point and frequency offset
图11. 传感器谐振点与频偏的拟合结果
的检测灵敏度,通过对模型尺寸的反复调整验证,检测灵敏度由最初的5.06 MHz/(100 mg/dL)提升至11.3 MHz/(100 mg/dL)。随着检测频点向高频移动,检测灵敏度呈现出线性的增长,数据的拟合结果也表明传感器的空载谐振点与检测灵敏度呈现高度的相关性。
基金项目
2023浙江省教育厅一般科研项目(项目编号:3100JYN4118001G-155)。