红外偏振光治疗仪偏振度测量方法研究

牛礼军1,王安意1,唐征宇2a,迟玉刚2b

1.中国电子科技集团公司第四十一研究所 研发一部,安徽 蚌埠 233000;2.蚌埠市第一人民医院 a.骨科;b.影像科,安徽 蚌埠 233000

[摘 要] 目前空间光的偏振度测量方法不适用于红外偏振光治疗仪偏振度的测量。本文阐述了偏振度的测量方法,利用现有的设备搭建测量光路,利用激光窄带滤光片模拟单色光,对红外偏振光治疗仪的偏振度进行了测量,测得红外偏振光治疗仪各个治疗头输出光的偏振度都在95%以上。分析了影响偏振度测量的各种因素,在今后的工作中可以改进测量方法,提高红外偏振光治疗仪偏振度的测量精度。

[关键词] 红外偏振光治疗仪;单色光;偏振度;测量光路

引言

红外偏振光治疗仪利用对人体透射效果最好的0.6~1.6 μm光谱范围的线偏振光照射人体。线偏振光输出功率达到2200 mW,对人体组织的有效作用深度超过5 cm,可以进行神经节、神经干和神经根的照射,从而赋予光疗以全新的治疗概念——神经照射法。红外偏振光治疗仪作用于人体相应的疼痛部位,有抑制神经兴奋、松弛肌肉、促进活性物质产生、加速致痛物质的代谢与清理和调节内环境稳定等作用[1]。红外偏振光治疗仪主要依靠光的偏振特性对人体产生影响,从而可以治疗不同类型的疼痛[2-6]。目前,红外偏振光治疗仪已经广泛应用于神经内科、疼痛科、康复理疗科、外科、骨伤科、皮肤科等科室。红外偏振光治疗仪有B、C、D和SG四种类型的治疗头,分别用于治疗不同的症状,要求四种治疗头输出光的偏振度都大于等于95%。偏振度是衡量红外偏振光治疗仪治疗性能的重要指标。本文阐述了红外偏振光治疗仪输出光偏振度的测量方法,并测得了各个治疗头的偏振度。

1 测量方法

偏振度测量方法有直接测量法和调制法两大类[7-15]。直接测量法采用检偏器,在不同的透射方向(通常为0°、60°、120°)测量入射光强,计算出目标光的偏振度。该方法原理简单,但受检偏器的透过轴方向定位精度、测量偶然性误差等因素的影响,测量精度一般在2%左右[16]。调制法利用光波的振幅、频率和强度等信息,使他们按调制信号的规律变化。强度调制方法采用在光路中引入起偏、相位延迟或者机械旋转以改变偏振方向的光学元件,通过测量调制信号求得斯托克斯参量。引入相位延迟光学元件的偏振度测量调制法的优点是快速同步获取目标光谱和偏振信息,相对于直接测量法,测量精度有所提升;缺点是引入了偏振度测量的不确定度贡献[17-19]。本文利用THORLABS公司的基于强度调制原理的偏振测量仪进行偏振度测量。

光的偏振态可用斯托克斯矢量表示[20]。斯托克斯矢量包含4个元素,既S=[S0+S1+S2+S3]T,各元素均可用光强表示,能够直接测量。利用斯托克斯矢量定义的偏振度表达为:

在大多数情形下,自然景物的反射光是部分偏振光,可看作是完全线偏振光和非偏振光的混合。根据研究数据,这类偏振光的斯托克斯矢量中,表征左旋和右旋偏振光强度差异的S3分量接近于0。偏振度定义可简化为线偏振度的定义:

斯托克斯参量的测量方法一般采用光强调制法,在待测光路中引入起偏器和相位延迟器,并对它进行调制,通过测量调制光强求得斯托克斯参量。,旋转波片法测量斯托克斯矢量的原理图如图1所示。让待测光依次通过一个匀速旋转的1/4波片和一个线偏振片,波片的角速度为ω

图1 旋转波片法原理图

入射光斯托克斯矢量为S=[S0+S1+S2+S3]T,旋转角度为θ时,1/4波片的穆勒矩阵为:

水平方向的线偏振片的穆勒矩阵为:

出射光的斯托克斯矢量为

代入θ=ωt,上式可改写为

其中

式(7)的形式为离散傅里叶级数的展开式,包括一个直流分量(A),一个2ω项(B),和两个4ω项(C和D)。对I傅里叶分析可以得到系数计算式如下:

其中θ=ωt,于是可以通过式(8)和(9)求解斯托克斯矢量。

红外偏振光治疗仪采用0.6~1.6 μm的连续光谱照射人体。目前,连续光谱的偏振度测量还是一个难点,还没有商业化的产品用于连续光谱的偏振度测量。本文所采用的红外偏振光治疗仪的偏振度测量方法如图2所示。由于连续光谱的偏振度还不能测量,本文采用不同波长的激光窄带滤光片模拟单色光,分别测量红外偏振光治疗仪在600、700、800 nm等波长处的偏振度。偏振态测试仪的测量头可以接受Ø3.0 mm的自由空间准直光束,要求入射光束为准直的平行光。本文设计了自由空间准直光路用于偏振度的测量。将红外偏振光治疗仪治疗头固定在光学平台上,通过一个可调光阑把红外偏振光治疗仪发光面模拟成为一个点光源,光束通过光阑后再通过一个准直透镜将光变成平行光,平行光通过激光窄带滤光片模拟成为单色光,最后光线进入偏振态测试仪测得红外偏振光治疗仪输出光的偏振度。

图2 偏振度测量示意图

2 偏振度测量实验及结果

偏振态测试仪选用THORLABS公司的PAX5710系列。该系列偏振测试仪采用旋转波片的方法,测量准直入射光的偏振态。该系列有三种型号可选,涉及可见光到近红外波段。红外偏振光治疗仪的发光光谱恰好位于可见光及近红外波段。偏振测量仪要求入射光必须为单色相干光,且光会在模块内终止。偏振态测量仪具有70 dB的高动态范围,邦加球上方位角和椭偏度的准确度为±0.25°。可变光阑选用THORLABS公司的拨杆控制的连续可变光阑,用来模拟点光源。准直透镜选用THORLABS公司的双胶合消色差透镜,用来模拟平行光。滤光片选用THORLABS公司的激光窄带滤光片。

为了验证窄带滤光片出射光的光谱宽度,本文利用Gooch&Housego的双单色仪OL750-M-D对滤光片的透过率进行了测量,600、700和800 nm滤光片的透过率如图3所示。由图3可知,滤光片在通带区域内具有比较高的透过率。虽然利用滤光片模拟单色光的光谱宽度比激光器输出光的光谱宽度要宽,但是利用窄带滤光片产生的光也可认为是单色光。

图3 不同波长滤光片的透过率测量曲线

注:a.600 nm;b.700 nm;c.800 nm。

按照红外偏振光治疗仪输出光的偏振度测量原理,搭建的测量系统如图4所示。按照此方法对红外偏振光治疗仪进行偏振度测量,以C型治疗头为例测量治疗头在600、700和800 nm的偏振度,测得结果如图5所示。红外偏振光治疗仪的C型治疗头在600、700和800 nm处的偏振度分别是95.6%、99.3%和96.6%。

图4 偏振度测量系统

图5 C型治疗头在不同波长处的偏振度

注:a.600 nm;b.700 nm;c.800 nm。

按照同样的方法对红外偏振光治疗仪的其他治疗头的输出光的偏振度进行测量,测得结果如表1所示。

表1 各治疗头偏振度测量结果(%)

SG C D 600 96.5 95.1 95.6 96.8 700 96.8 95.8 99.3 98.9 800 96.0 95.0 96.6 96.3波长 (nm) 治疗头B

在偏振度的测量过程中,有诸多因素影响偏振度测量。偏振态测试仪要求入射光是单色平行光,入射光束的平行度和单色性都会影响偏振度的测量。而本文采用模拟点光源,最后由准直透镜形成平行光的方法进行偏振度测量。点光源不是理想的点光源,由于装调误差准直透镜的出射光也不是严格的平行光,这些因素都会对偏振度测量造成影响。运用激光窄带滤光片模拟平行光,而滤光片的中心波长有±2 nm的误差,半带宽(10±2)nm,所以出射光不是严格的单色光,这样也会对偏振度的测量造成一定的影响。偏振态测试仪本身也会带来一定的测量误差。偏振态测试仪主要由1/4波片、检偏器和探测器组成,波片的相位延迟量偏差,快轴初始定位角度偏差与检偏器透光轴角度偏差是影响基于旋转波片法的斯托克斯参量检测的主要影响因子。在偏振度测量过程中应尽量减少测量误差提高偏振度测量的准确性。

3 结论

红外偏振光治疗仪要求每个治疗头输出光的偏振度都大于等于95%,其偏振度测量是检验红外偏光治疗仪性能优劣的关键步骤。本文通过对偏振度测量原理的阐述,模拟单色平行光,搭建测量光路,测得了红外偏振光治疗仪各个治疗头输出光的偏振度,为红外偏振光治疗仪的偏振度测量提供参考。但是在测量过程中还存在很多影响测量结果的因素,在以后的工作中可以改进测量方法,提高偏振度测量精度。

[参考文献]

[1]张琼晖,彭子明.直线偏振光近红外治疗仪[J].中国医疗器械杂志,2002,26(4):307-308.

[2]金海龙,田宏毅,洪文学,等.复合光物理治疗系统的光谱特性研究[J].中国生物医学工程学报,2008,27(1):152-156.

[3]Iordanou P,Baltopoulos G,Giannakopoulou M,et al.Effect of polarized light in the healing process of pressure ulcers[J].Int J Nurs Pract,2002,8(1):49-55.

[4]唐勇,何劼,刘瑞芳,等.红外偏振光治疗空军飞行人员轴性颈痛的疗效观察[J].中国疗养医学,2016,25(4):337-338.

[5]周沫.红外偏振光治疗仪在产科临床应用总结[J].现代医学仪器与应用,2007,19(4):42-44.

[6]郭继祥.红外偏振光照射星状神经节结合牵引治疗椎动脉型颈椎病疗效观察[J].河北联合大学学报(医学版),2013(2):206-207.

[7]Gordon R,Brolo AG,McKinnon A,et al.Strong polarization in the optical transmission through elliptical nanohole arrays[J].Phys Rev Lett,2004,92(3):037401.

[8]宋志平,洪津,乔延利,et al.一种目标光辐射线偏振度测量新方法[J].半导体光电,2008,29(5):759-761.

[9]Chen L,Meng F,Yuan Y,et al.High-precision variable polarization light source[J].Procedia Eng,2012,29:1835-1839.

[10]Berrry HG,Gabrielse G,Livingston AE.Measurement of the stokes parameters of light[J].Appl Opt,1977,16(12):3200-3205.

[11]Collett E.Measurement of the four Stokes polarization parameters with a single circular polarizer[J].Opt Commun,1984,52(15):77-80.

[12]支丹丹,李健军,高冬阳,等.基于旋转波片的斯托克斯参量检测与精度分析[J].光谱学与光谱分析,2016,36(8):2655-2659.

[13]覃兆宇,程兆谷,张志平,等.偏振光斯托克斯参量的高速实时测量技术[J].光学学报,2007,27(4):659-662.

[14]阴亚芳,方强,刘毓.偏振度测试方法研究[J].光通信研究,2005,(4):68-70.

[15]罗英达,陈哲.偏振光斯托克斯参量原理与测量方法设计[J].激光与光电子学进展,2009,46(6):56-58.

[16]袁银麟,郑小兵,吴浩宇,等.基于傅里叶变换的强度调制型光谱偏振分析仪[J].光学学报,2014,34(12):271-278.

[17]宋志平,洪津,乔延利,等.一种目标光辐射线偏振度测量新方法[J].半导体光电,2008,20(5):759-762.

[18]王新全,相里斌,黄旻,等.静态成像光谱偏振仪[J].光电子激光,2011,22(5):689-692.

[19]陈友华,王召巴,王志斌,等.弹光调制型成像光谱偏振仪中的高精度偏振信息探测研究[J].物理学报,2013,62(6):73-80.

[20]廖延彪.偏振光学[M].北京:科学出版社,2003.

Study on the Polarization Measurement Method of Infrared Polarized Light Therapy Apparatus

NIU Lijun1, WANG Anyi1, TANG Zhengyu2a, CHI Yugang2b
1.Research and Development Division, The 41st Institute of China Electronics Technology Group Corporation, Bengbu Anhui 233000, China; 2.a.Department of Orthopedics; b.Department of Imaging, The First People's Hospital of Bengbu, Bengbu Anhui 233000, China

Abstract: The current polarization measurement method is not suitable for the polarization measurement of the infrared light therapy apparatus.In the present study, the measurement method of polarization degree was described.The polarization of the infrared light therapy apparatus was measured using the existing equipment to build the measurement light path and using the laser narrow-band filter to simulate the monochromatic light.For each treatment head of the infrared light therapy apparatus,the polarization degree of the output light was above 95%.Meanwhile, various factors affecting the polarization measurement were analyzed, which would contribute to the improvement of the measurement method and measurement accuracy of polarization degree of infrared polarization light therapy apparatus.

Key words: infrared polarized light therapy apparatus; monochromatic light; polarization degree; measurement optical path

收稿日期:2018-09-30

修回日期:2018-11-05

基金项目:安徽省重点研究与开发项目(1704a0802144)。

通讯作者:王安意,高级工程师,主要研究方向为康复医疗光疗。

通讯作者邮箱:114958528@qq.com

[中图分类号] O4363.3;TH773

[文献标识码] A

doi:10.3969/j.issn.1674-1633.2019.05.009

[文章编号] 1674-1633(2019)05-0035-04

本文编辑 尹娟