用于微创手术探针的光纤力传感器设计
2022-10-28于昌新何彦霖何超江祝航威祝连庆
于昌新,何彦霖,何超江,祝航威,祝连庆
(北京信息科技大学 光电测试技术与仪器教育部重点实验室 北京 100192)
1 引言
近年来临床手术逐渐朝着微创化的方向发展,其中穿刺手术是微创外科手术中最为常见的手术之一[1-2]。穿刺手术通常在核磁共振成像(Magnetic Resonance Imaging,MRI)和计算机断层扫描(Computed Tomography,CT)成像的辅助下进行药物放置、组织活检和放射治疗[3]。目前的穿刺手术机器人多数只能起到导航和定位的作用,缺乏力反馈技术[4-5],而穿刺力的准确检测可以提高手术的成功率、灵活性和安全性[6]。
微创穿刺针的细长结构决定了传感器的体积和质量必须足够小,体内的应用方式使其必须抗电磁干扰,这种需求使得光纤传感器成为测量探针末端力的首选传感器。光纤传感技术自诞生以来,被广泛应用于医疗、工程、国防等诸多领域[7-9],可以用来检测外力、温度、应变、折射率、湿度等环境因素的变化[10-13],其中力作为一个重要的物理量成为测量的热点。
目前光纤力传感按传感原理可分为光强式、相位式和光栅式。2017年,Liu等人[14]提出一种基于悬臂锥形插入光纤微腔的应变力传感器,该传感器由一种特殊制造的法布里-珀罗(Fabry-Perot,F-P)腔构成,虽然灵敏度高,但其制备过程复杂、不易加工。2018年,Lai等人[15]提出一种用于测量肌腱张力的光纤力传感器,该传感器的灵敏度为24.28 pm/N,但其在解耦力时忽略了温度的影响。2021年,Sun等人[16]提出一种用于腹腔镜手术中握爪力检测的光纤力传感器,该传感器的测力范围为0~10 N,轴向力灵敏度为47.5 pm/N。2021年,Abdulfatah等 人[17]提出 一种基于长、短波长的光纤布拉格光栅(Fiber Bragg Grating,FBG)光纤传感器,研究了传感器对横向力和温度的变化特性,其长、短波长的力灵敏度分别为71.2 pm/N、57.3 pm/N。2022年,孙世政等人[18]提出一种双层十字梁结构FBG传感器,该传感器可以实现三维力传感并且具有较好的非线性解耦能力,但是其尺寸较大,不适合应用于微创手术中。上述这些力传感器虽然在某些方面具有不错的效果,但由于未能解决力与温度的交叉问题,难以满足微创外科穿刺手术探针力的测量需求。
本文针对以上问题,设计了一种基于F-P腔级联FBG的光纤力传感器,该传感器由F-P腔和FBG在单根光纤上级联形成,直径为0.25 mm,可集成于手术探针上实现探针尖端的力传感。通过研究传感器的力响应和温度响应对传感器进行标定,通过引入灵敏度矩阵的方式解决了力与温度的交叉影响问题。
2 光纤F-P腔级联FBG传感机理
本文采用的F-P腔级联FBG结构如图1所示。F-P结构传感的原理是双光束干涉,当一束光传播至F-P腔时,由于腔体与纤芯的折射率不同,使得在F-P腔的两个端面分别发生反射,产生光程差,形成双光束干涉效果。反射光的光程差为:
图1 F-P腔级联FBG结构Fig.1 Structure of F-P cavity cascaded FBG
其中:n为光纤纤芯的折射率,d为F-P腔的长度。第m级干涉条纹的波谷处干涉波长满足:
其中:m为整数,φ0为初始相位,λm为第m级干涉条纹波谷处的反射波长。由式(2)可得:
则光程差L与波谷波长λm成正比,对波长λm求导可得:
其中:ΔL为F-P腔的腔长变化量,Δλm为m级波谷波长变化量。
当对F-P腔施加轴向外力F时,会导致其腔结构上产生应变,由广义胡克定律有:
其中:ε为所受应变,F为所施加外力,A为端面的横截面积,E为光纤的弹性模量。
当F-P腔上发生应变时,应变效应和光弹效应会使其光程差发生变化,有:
其中:μ为光纤的应变系数,ρ为光纤的光弹系数。
空白试验是指除用纯水代替样品外,其他所加试剂和操作步骤与样品测定完全相同的操作过程。空白值的大小及其分散程度,对分析结果的精密度和分析方法的检出限度都有很大影响,并在一定程度上反映一个实验室及其分析人员的水平,如纯水质量、试剂纯度、量器及容器的洁净度、仪器的性能、滴定终点误差等对试验结果的影响。所以空白实验值应符合质控要求,否则就要从上述各方面查找原因。
根据式(4)~(6)有:
当F-P腔上温度发生变化时,热光效应和热膨胀效应会使其光程差发生变化,即有:
其中:α为光纤的热光系数,ζ为光纤的热膨胀系数。
根据式(4)和(8)有:
对于FBG,当应变ε和温度T同时作用时,有:
其中:λ0为初始布拉格波长,ΔλFBG为FBG的波长漂移量,pe为光纤的光弹系数,αΛ为光纤的热膨胀系数,αn为光纤的热光系数。
根据式(5)和(10)有:
当外力F和温度T同时作用于传感器上时,光纤F-P腔和FBG的波长漂移对二者的响应可表示为矩阵:
其中:KF1、KT1分 别 为F-P腔 的 力灵 敏 度 系数 和温度灵敏度系数,KF2、KT2分别为FBG的力灵敏度系数和温度灵敏度系数。
对式(12)求逆可得:
式(13)表明,当测得光纤传感器对力和温度的灵敏度系数后,通过测量F-P腔的m级波谷波长漂移量和FBG的特征波长漂移量,即可实现探针尖端受力的测量,同时也可对探针所处环境的温度进行测量,并且避免了二者之间的交叉影响。
3 实验系统建立
本文采用的实验系统包括力控制系统(如图2)和光纤传感系统(如图3)两部分。力控制系统主要包括手术探针、位移平台、六维力传感器、数据采集卡。手术探针为镍钛合金材料,直径为1 mm,长度为400 mm,内孔直径为0.3 mm,用于植入传感器件。位移平台采用微纳光科公司的WN103TM13H型手动位移平台,该平台的行程为13 mm,精度为10 μm,六维力传感器型号为NaNo17(SI-25-0.25,ATI Industrial Automation,NC,USA),数据采集卡型号 为USB-6210(National Instruments Corporation,USA)。光纤传感部分主要包括光谱仪、光源、上位机、光纤。光源为ASE光源,光平坦度小于2 dB。光谱仪为YOKOGAWA AQ6375,最小分辨率0.05 nm,波长范围1 200~2 400 nm,功率为-70 dBm~+20 dBm,快速测量时间为0.2 s,跨度为100 nm。
图2 力控制系统Fig.2 Force control system
图3 光纤传感系统Fig.3 Fiber sensing system
对于传感器的F-P腔结构,本文采用化学腐蚀的方式进行刻蚀,首先将SM-28单模光纤的一端剥去涂覆层后擦拭干净,将端面切平后浸入浓度为40%的氢氟酸之中,在通风厨中腐蚀20分钟。由于纤芯与包层中的硅、锗含量存在差异,使纤芯的腐蚀速度大于包层,从而使得在光纤纤芯上呈现一个凹槽结构,将腐蚀完成的光纤用清水冲洗干净,去除端面残留的氢氟酸。将腐蚀好的光纤凹槽与一根切平的单模光纤在熔接机(80S,Fujikura)中 熔 接,制 备 的F-P腔 腔 长 约33.4 μm,结构如图4所示。本文采用紫外光刻写掩膜版法制备传感器中的FBG结构,利用波长为248 nm的紫外光激光器与相位掩膜版相结合,在经载氢后的SM-28单模光纤上刻写,FBG的反射波长为1 549.7 nm、反射率为80%、3 dB带宽为0.3 nm,栅区长度为10 mm。
图4 F-P微腔结构Fig.4 Structure of F-P cavity
制备完成的F-P微腔与FBG相距3.5 cm熔接成一根光纤,形成光纤力传感器。为了将该传感器集成于探针之中,首先采用3M公司的DP2216环氧胶,将A、B胶按照等比例搅拌均匀,静置使气泡完全消失以提高力的传递效率,然后将F-P腔级联FBG传感器结构置于探针之中,使用制备好的胶将传感器密封于探针之内。
在实验室温度下测得无外力时传感器的反射光谱如图5所示。传感器的光谱范围为1 520 nm~1 610 nm,FBG的中心波长为1549.7 nm,F-P腔的反射光谱选为1 584.4 nm,整个光谱的最大对比度为-32.3 dB。
图5 传感器反射光谱Fig.5 Reflection spectrum of the sensor
4 传感器标定实验测试
4.1 轴向力标定
为了研究F-P腔与FBG对轴向力的敏感程度,对传感器进行了轴向力标定实验。实验过程中先将ATI置零,然后扭动微分头,每次旋转0.1 mm对应拉力约为0.55 N,分别记录ATI测得的探针尖端受力值与光谱仪上对应的光谱。施加拉力范围为0 N~5.5 N,测力步长为0.55 N。F-P腔的特征光谱和力与波长关系分别如图6和图7所示;FBG的特征光谱和力与波长关系分别如图8和图9所示。
图6 力变化时光纤F-P腔特征光谱Fig.6 Reflective spectrum of F-P cavity during force variation
图7 光纤F-P腔力与波长关系Fig.7 Relationship between force and wavelength for fiber F-P cavity
图8 力变化时FBG特征光谱Fig.8 Reflective spectrum of FBG during force variation
图9 FBG力与波长关系Fig.9 Relationship between force and wavelength for FBG
图6和图8表明,光纤F-P腔和FBG对力的响应趋势相同。在受力由0 N增加至5.5 N再减少到0 N的过程中,光纤F-P腔和FBG的特征波长均出现先红移再蓝移的现象。其中F-P腔的特征波长漂移1.9 nm,力灵敏度为331.8 pm/N,线性度大于0.98;FBG特征波长漂移0.84 nm,力灵敏度为159.9 pm/N,线性度大于0.99。
4.2 温度标定
为了研究F-P腔与FBG对温度的敏感程度,对传感器进行了温度标定实验。将嵌有光纤传感器的探针固定在手动位移平台上,然后放入高精度温度控制箱中,在20℃~50℃范围内,温度每变化5℃采集一次传感器的光谱。F-P腔的特征光谱和温度与波长关系分别如图10和图11所示;FBG的特征光谱和温度与波长关系分别如图12和 图13所 示。
图11 光纤F-P腔温度与波长关系Fig.11 Relationship between temperature and wavelength for fiber F-P cavity
图13 FBG温度与波长关系Fig.13 Relationship between temperature and wavelength for FBG
图10和图12表明,该光纤传感器中的F-P腔与FBG对温度的响应具有相同的趋势。在温度由20℃升高至50℃再降至20℃的过程中,光纤F-P腔和FBG的特征波长均出现先红移再蓝移的现象。其中,光纤F-P腔的波长漂移量为1.06 nm,温度灵敏度为33.7 pm/℃,线性度大于0.99;FBG的特征波长漂移为0.49 nm,温度灵敏度为16.6 pm/℃,线性度大于0.99。
图10 温度变化光纤F-P腔特征光谱Fig.10 Reflective spectrum of F-P cavity during temperature variation
图12 温度变化FBG特征光谱Fig.12 Reflective spectrum of FBG during temperature variation
根据测得的温度与力的灵敏度系数,代入式(13)可得到该光纤传感器对力与温度的响应矩阵方程为:
在实际测量过程中,当外界力与温度变化时,通过光谱仪测得光纤F-P腔与FBG的波长漂移量,代入式(14)即可求得不受温度影响时探针所受轴向力。
4.3 稳定性测试
为了测试传感器的稳定性,将集成传感器的探针固定于位移平台上,位移平台置于光学平台之上,实验在实验室室温下进行。每10分钟监测一次传感器的反射光谱,共监测10组,通过分别对比FBG和F-P腔的反射波长变化来说明传感器的稳定性,采集的光谱稳定性如图14所示。
图14 光谱稳定性测试Fig.14 Stability test of reflection spectrum
由图14可以看出FBG的最大波长漂移为0.01 nm,F-P腔的最大波长漂移为0.04 nm,平均波长漂移量为25 pm,结果表明该力传感器的稳定性较好。
5 结论
针对穿刺手术探针末端力测量的实际需求,提出并设计了基于F-P腔级联FBG的光纤力传感器,同时对其性能进行了实验测试,结果表明:
(1)该光纤力传感器在0 N—5.5 N—0 N加载卸载过程中,光纤F-P腔和FBG的特征光谱均出现先红移再蓝移的现象。F-P腔力灵敏度为331.8 pm/N,FBG力灵敏度为159.9 pm/N。
(2)该光纤力传感器在20℃—50℃—20℃升温降温过程中,光纤F-P腔和FBG的特征光谱均出现先红移再蓝移的现象。F-P腔温度灵敏度为33.7 pm/N,FBG温度灵敏度为16.6 pm/N。
(3)该光纤传感器在90分钟内具有良好的稳定性。其中FBG最大波长漂移为0.01 nm,光纤F-P腔最大波长漂移为0.04 nm,平均波长漂移量为25 pm。
本文提出的F-P腔级联FBG光纤力传感器不仅可以集成于探针之中实现其末端力测量,还可避免温度对探针受力的交叉影响,实现高灵敏度和高线性度,在微创外科穿刺手术机器人穿刺针力反馈测量领域有广阔的前景。在后续工作中将进一步研究在多场并存时的力传感方法和传感器性能。