引言

入射激光脉冲在光纤中传播时产生拉曼散射，拉曼散射包含斯托克斯散射和反斯托克斯散射.因为反斯托克斯散射光对温度敏感，斯托克斯散射光对温度不敏感，可以用斯托克斯散射光解调反斯托克斯散射光，利用二者的比值得到温度信息，这是传统的基于拉曼散射的分布式光纤温度传感器的基本原理^[1-4].但是这种方法存在一些不足之处:不同的光纤类型斯托克斯散射光和反斯托克斯散射光的波长有100~200 nm的差值，由于波长不同，两者在光纤中传播时的损耗不同，在光纤不同距离下两者的比值是关于距离的一个函数，而不是一个定值，这就给测量带来了误差.

针对传统方法存在的问题^[5]，文献[6]提出了双光源的方案，该结构采用了两个不同波长的光源.主激光器的入射激光波长和副激光器的斯托克斯散射光的波长相匹配，这样斯托克斯散射光的波长和反斯托克斯散射光的波长相同，消除了两者的损耗差.但是由于采用双光源增加了整个系统的成本，而且特殊波长不容易相匹配.文献[7-8]提出了双端结构，只利用反斯托克斯散射光，消除了波长差引起的损耗偏差.这种结构采用一个光源的激光脉冲从光纤两端入射的办法，增加了系统的复杂性.文献[9]对该结构进行了优化，提出了光纤末端增加反射镜的方法，同样解决了斯托克斯散射光和反斯托克斯散射光波长差引起的损耗差问题.这种结构相对简单，但是依赖反射镜的反射率，而且探测到的信号中除了反斯托克斯散射光的有效成分还包含了很多的其他成分，增加了后期信号处理难度.本文提出了一种基于拉曼散射的折叠光纤结构，在两根相同光纤的末端将两根光纤相互熔接即可实现折叠结构.两根光纤所处环境相同，激光脉冲从入射端进入光纤之后，在温度变化处可以有两次拉曼散射，接收到两次反斯托克斯散射光.由于只利用反斯托克斯散射光，所以消除了斯托克斯散射光和反斯托克斯散射光波长差引起的损耗偏差.

1 基本原理

如图 1所示，折叠后的两根光纤紧密排布在一起置于测量场内，l处为温度变化区域，由于光纤的折叠结构，2L-l处也处于温度变化场内(2L为整个传感光纤长度).激光器发射的激光脉冲通过光旋转镜进入传感光纤，在l处产生第一次反斯托克斯背向拉曼散射，在入射光经过2L-l处产生第二次反斯托克斯背向拉曼散射.由于两次反斯托克斯背向拉曼散射发生在同一温度变化区域，所以两次反斯托克斯背向拉曼散射光强相同，最终两次反斯托克斯背向拉曼散射光被探测器接收.

l和2L-l处温度场相同，两位置处的后向反斯托克斯拉曼散射光强分别为

$ I\left( l \right) = {K_{{\rm{as}}}}S\nu _{{\rm{as}}}^4{I_0}R\left( T \right)\exp \left[{-\left( {{\alpha _0} + {\alpha _{{\rm{as}}}}} \right)l} \right], $

(1)

$ I\left( {2L- l} \right) = {K_{{\rm{as}}}}S\nu _{{\rm{as}}}^4{I_0}R\left( T \right)\exp \left[{-\left( {{\alpha _0} + {\alpha _{{\rm{as}}}}} \right)\left( {2L-l} \right)} \right], $

(2)

式中：R(T)是光纤l和2L-l处温度信息场函数，$R\left( T \right) = {\left[{\exp \left( {\hbar \Delta \nu /\kappa T} \right)-1} \right]^{ -1}}$，$\hbar $为普朗克常量，Δν为拉曼散射的频移量，κ为玻尔兹曼常量；K_as为反斯托克斯散射截面；S为背向散射因子；ν_as为反斯托克斯散射的光子频率；I₀为入射光强；α₀和α_as为光纤中入射光和反斯托克斯散射光的衰减系数.

将(1) 式和(2) 式相乘，可得

$ I = \sqrt {I\left( l \right)I\left( {2L- l} \right)} = {K_{{\rm{as}}}}S\nu _{{\rm{as}}}^4{I_0}R\left( T \right)\exp \left[{-\left( {{\alpha _0} + {\alpha _{{\rm{as}}}}} \right)L} \right], $

其中：exp[-(α₀+α_as)L]是与探测距离无关的一个常量，最后得到的任何一点的光强损耗都为入射光和反斯托克斯散射光在整根探测光纤中的损耗.选择距离起始处为参考点，参考温度为T₀，可得解调温度公式为

$ T = {\left( {\frac{\kappa }{{\hbar \Delta \nu }}\ln \left( {\frac{{I\left( {{T_0}} \right)}}{{T\left( T \right)}}\exp \left( {\frac{{\hbar \Delta \nu }}{{\kappa T}}-1} \right) + 1} \right)} \right)^{-1}}. $

(3)

该温度解调方法只利用了反斯托克斯散射光，所以不存在斯托克斯散射光和反斯托克斯散射光由于波长差引起的损耗差问题，而且利用光纤折叠结构使得温度变化区域处的反斯托克斯散射光在光纤中的损耗为在整根传感光纤中传播的损耗，且不随距离变化而改变，有效减少了误差.

2 结果与分析

由于反斯托克斯散射强度相对较小，探测器探测到的信号中含有大量的噪声.噪声的主要成分为探测器的热噪声，即白噪声，所以需要对探测到的信号进行降噪处理^[9-12].由于信噪比较小，有效信号被噪声所淹没，本文中采用两步降噪法对测量结果进行噪声处理，先累加平均再小波变换.首先对探测到的信号进行时域累加平均预处理.激光器发出N个激光脉冲，探测器探测到N个光强分布曲线，对N个光强分布曲线进行累加平均得到预处理后的信号，采用的公式为$I = \left( {\sum\limits_{i = 1}^N {{I_i}\left( {l, T} \right)} } \right)/N$.

光纤总长(2L)为8 km, 环境温度为20 ℃，光纤长度2.5 km和5.5 km处探测场的温度为60 ℃，对探测到的光强分布累加15 000次平均处理后得到的曲线如图 2所示，图 2中光强数值为归一化处理的结果.可以看出，光强变化的两处关于4 km对称，传感光纤上光强的分布随着距离的增加而减小.把探测到的光强信号带入到式(3) 中得到温度分布曲线，结果如图 3所示.可以看出，经过预处理的信号信噪比还是非常低，需要进一步的降噪处理.

由于温度分布曲线上的探测点以单一脉冲方式随机出现，单独的频域或者时域滤波效果都不好，所以本文采用时频滤波的方法小波变换去噪，用db5小波对预处理后的信号6层小波分解，得到1至6层小波系数细节如图 4所示.

小波变换模极大值法能有效保留温度突变点的信息，所以用此方法对不同细节分量d_i小波系数进行处理，阈值的确定采用软阈值方法，重构得到信号如图 5所示.可以看出，在2.5 km处的温度为61.5 ℃，在5.5 km处的温度为62.7 ℃，误差分别为1.5 ℃和2.7 ℃，没有随着距离的增加而出现较大的误差.光纤的前端和后端环境温度基本保持在20 ℃左右，系统测量误差±1.7 ℃，空间分辨率为1 m，不存在波长差引起的损耗差使得光纤上的温度分布随着距离的增加而减小.因此，使用折叠光纤结构使整条光纤上的温度分布得到了很好的矫正.

3 结论

对基于拉曼散射的分布式光纤温度传感器的温度解调方法进行分析，提出了一种光纤折叠的结构.传统温度解调的方法使用反斯托克斯散射和斯托克斯散射的比值来得到温度，但是这两种散射光波长不同，在光纤中传播会有损耗差.使用折叠光纤的方法使激光脉冲在温度场内有两次拉曼散射，从而可以只探测对温度敏感的反斯托克斯散射光，消除了损耗差，使温度的测量精度得到提升.由于探测信号存在大量噪声，提出了先用累加平均方法对信号预处理，然后再用小波变换的方法降噪，信噪比由1 dB提升到11 dB.