1. 引言
食品安全是当今社会的一个重大问题,食用油作为人们不可缺少的食品,关系着人们的健康。尤其是目前出现了很多“地沟油”出现在食用油市场,严重的危害了人们的健康和生活品质。针对这一问题,研究学者们提出了多种方法来鉴别食用油品质,有光学,电化学分析、分子生物学等技术分析方法 [1] [2] [3] 。光学方法中使用较多的有太赫兹波段光谱 [4] 、红外光谱 [5] 和荧光光谱 [6] 分析等。这些光谱分析法主要通过对比观察正常油品与劣质油品的吸收光谱或发射光谱之间的差别来进行油品的品质鉴别,但实际的精炼地沟油或同一品种的复热油分子结构近似,与正常油品之间的差异很难被区分出来 [7] ,需要处理样品及专业的分析。
表征使用油品质的指标有很多,折射率是其中一个重要指标 [8] 。一般情况下,油品的混合浓度、油品是否掺假可以通过其折射率之间的差别来进行有效的甄别 [9] 。通常油品的折射率可以通过阿贝折射率仪、分光仪和干涉法来进行测量 [10] [11] [12] 。阿贝折射仪通过全反射临界角的测定来测量液体折射率,不需要处理样品,但测量精度只有0.005,且清洗过程复杂。利用分光计最小偏向角法需待测液加工成三棱镜状,测量时间长。采用干涉法利用折射率引起的干涉条纹的变化来检测液体折射率,快捷,较高分辨率,但是光路需要仔细的调节,对操作人的技术水平有一定要求。
光纤迈克尔逊干涉仪将光干涉法在空间传输的光路集成在光纤内,不需要进行精细的调整,即可形成干涉 [13] 。本文提出了一种通过观测谐振峰的中心波长漂移来快速测量食用油折射率的方法:在常见光纤迈克尔逊干涉臂后连接一段玻璃毛细管,待测油品装入毛细管内,油品折射率的不同,形成的干涉光谱的谐振峰中心波长会发生改变。模拟和测试表明,这种方法可快速得到不同样品的折射率。
2. 测量装置及原理
迈克尔逊干涉仪是光学实验中的一种经典测量方法,其原理简单、测量精度高,但是干涉条纹不容易分辨,读数困难。光纤迈克尔逊干涉仪将光路集成在光纤中,通过检测光谱的变化来感测外界变化。但是由于油质折射率要高于光纤纤芯折射率,因而不好测量,我们对普通光纤迈克尔逊干涉仪的结构进行了改造,如图1所示(圆框内把结构局部放大,以方便观察)。在传感臂光纤的尾端连接一段毛细玻璃管作为测量样品区,并在毛细玻璃的尾部紧贴一个平面镜来保证每次测量样品的长度,以及提高反射率。
光源的光经过一个2 × 2的光纤耦合器,在另外一端被分成两束光。其中一个光纤端面切平,光在光纤端面部分光反射回耦合器,另外一路在光纤切平的端面后连接了毛细玻璃管,这样会形成两个反射面,一个是光纤和待测油形成的界面,另外一个是待测油后的平面镜,这两个反射面都会反射部分光返回耦合器,经过耦合器合束后,由光谱仪接收检测。实验时只需要将毛细玻璃管内注入需要检测的油样品即可,每次更换样品前将样玻璃毛细管清理干净,如图2所示。
光源采用一个宽带连续光I,经过2 × 2耦合器后分为两束光,其中一束进入参考臂,在参考臂末端反射后,有部分光I1又返回耦合器;另外一束进入传感臂,在光纤和油的端面,有部分光I2又返回耦合器,还有一部分光会进入玻璃毛细管中,在平面镜处反射返回耦合器,为I3。三束光可分别表示为
(1)
(2)
(3)
式中,q为耦合器的分光比,n1为光纤的纤芯折射率,n2为玻璃管内待测油的折射率,α为光在玻璃管内的传输损耗,l1、l2和l3分别为光纤1和光纤2和玻璃管的长度,l1,l2,l3分别为三束光走过的距离,r1,r2,r3分别为三个端面的反射率,根据菲涅尔公式,反射率取决于端面两端的折射率,可表示为:
(4)
(5)
图1. 提出的光纤迈克尔干涉仪结构示意图
图2. 提出的迈克尔逊干涉仪干涉原理示意图
假设平面镜的反射率r3 = 0.99。三束光经过耦合器的耦合同时到达探测器,探测器接收到的相对光强可表示为:
(6)
光纤纤芯的折射率n1为1.4682,而合格的油的折射率约为1.472,因此r2的值大约为10−6,非常的小,因此,式(6)的第5项起主要作用。可以看出,探测器最后接收的光强I4关于入射波波长λ是a为周期的余弦函数,其中a为
(7)
当入射波波长满足下式的条件时,光强达到最小值
(8)
式中,λN称为第N个谐振峰的中心波长,可以看出λN与待测油的折射率成正比,可用于检测待测油折射率的变化。
3. 模拟分析
图3(a)是利用波传输法模拟光在装有待测油的玻璃毛细管内的传输情况。模拟设置玻璃毛细管内径130 μm,外径250 μm,待测油折射率1.4726,长度2 mm。从图3(a)中可以看出,光在光纤内传输损耗几乎为0,但在一进入待测油区域,光强会迅速分散降低,在1 mm后由于毛细管管壁的反射作用会有略微的起伏变化,传输距离与最后剩余的光强的关系如图3(b)所示。从图3(a)中可以看出,剩余的光强在传输距离大约为1.5 mm时最大,因此后面的模拟设置待测油区域长度为1.5 mm,传输损耗
。
式(6)的第5项的系数决定了接收到余弦峰的对比度,对比度越大越容易观察。根据式(6)可以看出,当
时,光谱中谐振峰的对比度最大,代入数据可得,约在q = 0.13时,谐振峰对比度最大,根据常见的耦合器比例,我们选取最接近的q = 0.1。利用matlab对探测接收到的光谱进行
(a) (b)
Figure 3. (a) Beam propagation in the glass capillary filled with oil; (b) optical power change with the length of the glass capillary
图3. (a) 光在装满油的玻璃毛细管内的传输情况,(b) 能量随玻璃毛细管的长度变化
(a) (b)
Figure 4. (a) Beam propagation in the glass capillary filled with oil; (b) optical power change with the length of the glass capillary
图4. (a) 光在装满油的玻璃毛细管内的传输情况,(b) 能量随玻璃毛细管的长度变化
模拟分析。模拟结构参数设置:光源为1250 nm~1650 nm的宽带光源,n1 = 1.4682,r3 = 0.99,l1 = 201.5 mm,l2 = 200.0 mm,l3 = 1.5 mm,
。图4(a)为模拟得到的待测油的折射率分别为1.4726、1.47262、1.47264、1.47266、1.47268、1.4727时接收到的光谱。可以看出在1250 nm到1650 nm范围内,大约有3个峰,每个峰的对比度约9 dB,且峰随着待测油折射率的升高向长波方向漂移。图4(b)为1550 n附近峰的中心波长随待测油折射率的变化。可以看出,峰的中心波长随油的折射率呈线性变化,具有较好的线性度(0.9998),且具有较高的灵敏度35418.8 nm/RIU。探测器为光谱仪,分辨率为0.02 nm,那么用提出的光纤迈克尔逊干涉仪结构测量折射率最小分辨率0.02/35418.8 = 5.6 × 10−7,即分辨率小于0.000001 RIU。
4. 结论
对于同一种油品,例如油中含杂质量不同,或重复加热,其折射率会相应的改变,用光纤迈克尔逊干涉仪进行检测时表现为光谱中峰的中心波长的变化。因此,以已知的合格产品作为标准,就可以快速知道某种食用油是否有掺假、不达标的情况。
本文检测方法原理简单,精度高,为食用油品质的检测或地沟油的快速鉴别提供了一种简便、实用的方法。
基金项目
大学生创新创业训练计划项目,201610060045。