杨汉瑞 李勇勇 徐士博 张经纬 栾 宁

随着我国供热面积不断扩大以及供热管道运行时限逐年增加,管道老化变形等问题日益加剧,泄漏问题时有发生,这对城市交通以及居民的生命财产安全造成了极大的威胁,大量人力物力的投入更加剧了资源的浪费.其中热力管道中膨胀节作为一种能根据外力或温差而发生弹性伸缩的补偿元件,极易因形变过大而致使管道产生泄漏甚至炸裂等安全事故.然而目前对于膨胀节健康状态检测的相关研究也仅局限于实验检验.而对于膨胀节疲劳实验前后的波距变化率检测更是处于起步甚至空白状态,膨胀节膨胀量等状态参数的准确测量问题亟待解决.

膨胀节结构复杂,其检测较为困难,相关研究成果尚少.文献[1]从理论上系统地分析了波纹管补偿器在供热管网中的选择方法,并给出选择依据取得了良好的效果.Chugunov 等在文献[2]中利用SolidWorks 软件依据相关标准进行波纹管优化设计,提高了设计效率.Zhang 等在文献[3−5]中分别对波纹管的稳定性进行了评估,分别对管道不同部件进行力学分析,为膨胀量大小进行了限定.以上文献中主要是就波纹管设计和材料力学性能等方面展开的研究.文献[6]中崔亚飞等提到的膨胀节波距测量方法采用的是人工在波纹管上划线标记,完成疲劳测试后利用游标卡尺对膨胀节内波纹管的波距变化进行测量,进而得出膨胀量.这种方法效率低下,检测误差较大,对测量人员的专业技能要求较高.另外文献[7]中王晓嘉提到利用三坐标测量机进行波距测量,但是针对大型膨胀节的测量,由于三坐标测量机价格昂贵就需要付出更高的成本.此外,还有人提到利用激光器对波距进行检测的方法,这种方法虽然在一定意义上弥补了之前膨胀节检测方法的不足,但是复杂环境下设备的安装以及造价等问题也是需要面临的现实问题.为此,亟待提出新的检测方法.

分布式光纤温度传感技术是利用传感光纤中光强大小与温度变化之间的关系进行温度测量的一种先进技术[8−11].光纤既作为光传输介质又作为温度传感介质,具有体积小,质量轻,抗干扰能力强等优点,另外光纤还能够适应不同的复杂环境[12−13].目前依据研究机理主要有三种系统,分别为基于瑞利散射、基于布里渊散射以及基于拉曼散射的分布式光纤温度传感系统.相较而言,基于拉曼散射的分布式光纤传感系统测量精度较高,造价较低并且有长距离测量的相关研究成果.近年来,分布式光纤传感技术在桥梁、管道、高架线等领域的应用在不断增加[14−15].而其在供热管道上的应用还较少,因此利用其分布式连续测温的优点将该系统应用于供热管网来检测膨胀节及各个管段的健康运行状态具有很高的潜在价值.同时,温度信息对于电厂热力管道的运行以及参数调节有指导性的意义.因此,为推进分布式光纤测温技术在热网健康状态中的应用,本文提出了一种基于分布式光纤传感技术的膨胀节膨胀量检测的新方法,通过实验验证了方法的可行性,并对实验数据进行了分析给出了利用分布式光纤温度传感系统检测膨胀节膨胀量准确性的量化指标.

1 理论基础

1.1 基于拉曼散射的分布式光纤测温原理

拉曼分布式测温系统基本组成如图1 所示,主要包括:电源模块,光源模块,分光模块,传感模块,光电转换模块,数据采集模块[16].光源在控制单元发出指令后发出光脉冲,通过耦合器进入传感光纤.光在光纤中传播时会发生散射,携带有温度信息的后向散射光经波分复用器(Wavelength division multiplexer,WDM)后,产生两束波长不同的光,即斯托克斯光和反斯托克斯光,这两束光通过雪崩光电二极管(Avalanche photodiode,APD)进行光电转换,进而被高速采集卡所采集,将采集到的数据经过处理后得到不同位置的温度信息[17].

1.2 温度解调理论

本部分主要分析系统的温度解调方法以及位置信息的获取理论,为后续推导并建立膨胀节膨胀量检测模型提供必要的依据.

图1 系统组成图Fig.1 System composition diagram

设注入光纤中的光脉冲峰值功率为P0,光脉冲沿光纤传输到距离L处产生Raman 散射(包括斯托克斯光和反斯托克斯光),在光脉冲注入端可以探测到后向斯托克斯光与反斯托克斯光,其光功率可表示为:

式中:νs、νas分别为斯托克斯光和反斯托克斯光在光纤中传播的频率,P0、Ps(L)、Pas(L)分别为入射光、斯托克斯和反斯托克斯信号光功率,Ks、Kas为与两种散射有关的系数,L为传感光纤的轴向距离,α0、αs、αas分别为光纤中入射光、斯托克斯光、反斯托克斯光的衰减系数,Rs(T)、Ras(T)分别为斯托克斯和反斯托克斯的玻尔兹曼因子,表达式如下[14−15]:

其中:∆ν为光纤分子振动频率,h是planck常数,k为boltzman常数,T(L)为距离L处的温度值.斯托克斯光信号与反斯托克斯光信号的功率之比为:

为消除光源振动、微弯损耗、光纤接头损耗等引起的光功率波动对温度测量精度的影响,需要利用斯托克斯光对反斯托克斯光进行解调,因此传感系统通常使用双光路探测.计算方法如下.

将光纤移至参考温度为T0的情况再次测定,此时得到的比值是:

当参考温度变为T时,得到的比值为:

从而将式(6)和式(7)作比之后导出温度解调表达式:

由光时域反射原理可知,测量出光纤轴向距离Z处的温度值后,利用Z与时间t以及激光在光纤中的传播速度v的关系,得出测量点的位置信息:

2 膨胀节膨胀量检测方法

2.1 膨胀节组成和工作原理

图2 为膨胀节半剖图,由图可以看出膨胀节主要由波纹管、接管、内衬管、垫环、耳板、双头螺柱以及螺母组成.其作用是缓冲热力管道因膨胀或收缩而发生的形变[4],其膨胀量计算的理论公式是:

其中,α为线膨胀系数(钢铁一般为(12∼18)×(10−6m)/(m·◦C)),L为所需补偿管道两个支架之间的距离,单位是m.∆T为注水时的水温与环境温度的差值,单位是◦C.

与之对应的,管道在受热时产生的热应力为:

式中,σ的单位是MPa,E100为100◦C 时钢的弹性模量.管道因膨胀产生的应力为:

式中,F为流体所在管道的横截面积,单位是cm2.

图2 膨胀节半剖图Fig.2 Semi caesarean section of expansion joint

由上述理论公式可以看出,由于管道内温度的变化会使得在热应力与膨胀应力的作用下产生一个膨胀量,同时对支架也产生一个轴向应力,严重时会使管道畸形[5].通常情况下膨胀量的大小受水温与环境温度差值变化影响较大,因此有必要对膨胀节的膨胀量进行检测,提前预判防止事故的发生.

2.2 检测方法

结合前面的分析,论文提出依据分布式光纤温度传感系统测量温度突变点与光纤线缆位置的对应关系来建立膨胀量检测模型的方法,测量方法示意图如图3 所示.对光纤线缆上某两点加热,通过分布式光纤温度传感系统对信号采集解调处理后,便可以在上位机显示界面温度分布曲线上看出这两点的温度突变点间的实际距离.而膨胀节发生伸长或者收缩时,会带动膨胀节环板10 前后移动,环板的温度与管道6 中水温相当.基于这个温度信息,焊接金属导热体4 至感温光纤线缆2.一方面,由于膨胀节本身伸缩量较小,另一方面受到分布式传感系统空间分辨率的影响,需要对光纤线缆的长度进行“放大”— 即利用柱形绝热体3 绕制光纤线缆2,以此来对微小的膨胀进行几十倍甚至上百倍的放大.这样,给定膨胀节环板所导出温度的初始位置,当膨胀节产生伸长量或者收缩量时,便可以从上位机1 上读出末端位置.这样膨胀量的变化距离就转化为光纤线缆的温度变化点之间距离,只要知道温度变化点之间的距离就能计算得到膨胀量值[18],两者具体关系见第2.3 节模型建立.

2.3 模型建立

图4 所示为膨胀量检测原理图,根据图中所标定的物理量可以对分布式光纤温度传感系统解调后的温度变化信息及位置信息建立膨胀量检测数学模型.通过1.2 温度解调理论可知,温度变化信息反映于斯托克斯光强与反斯托克斯光强之中,利用式(8)即可实现对温度变化点的温度进行解调.位置信息能够利用光时域反射原理进行测定,由于光纤线缆长度Z与温度变化点间的距离L相对应,因此利用式(9)即可得到温度变化点间的距离L.

设定膨胀量为K,柱形绝热体的直径为D,光纤线缆的宽度为Q,环绕在柱形绝热体上的感温光缆圈数为N,由图4 可知,温度突变点间的距离L与柱形绝热体直径D之间存在线性关系:

当膨胀节的膨胀量发生变化时,金属导热体触头便会产生位移,如图3 所示.由图4 可知,膨胀量K与绕制于柱形绝热体上的光纤线缆的宽度Q及圈数N存在着如下关系:

从式(14)可知,当光纤线缆固定即光纤线缆的宽度Q固定,只需知道温度突变点间的距离L以及柱形绝热体的直径D便可得到膨胀量K的值.该方法的可行性以及各参数之间的联系详见下方实验分析.

图3 测试系统示意图Fig.3 The diagram of test system

3 实验与结果分析

3.1 实验

采用分布式光纤温度传感系统作为测量媒介,温度突变采用抽取恒温水浴箱中温度恒定的水进行模拟,将抽水泵一端系入恒温水浴箱,出水端口压制为2 mm 宽的水缝,延长至柱形光纤信号采集器(绕有光纤线缆的柱形绝热体)表面.将恒温水浴箱中水加热至一定温度,利用抽水泵沿柱形光纤信号采集器自左至右缓慢移动,观察上位机界面温度曲线显示情况,并选取特定点进行计算验证.分布式光纤温度传感系统,型号为DTS-10,空间分辨率为1 m,温度分辨率为1◦C,温度范围为−20◦C∼100◦C.具体实验器材如表1 所示.

图4 检测原理图Fig.4 The detection schematic diagram

表1 实验器材表Table 1 Experimental equipment table

将分布式光纤温度传感系统标定校准后,选定柱形光纤信号采集器一端光纤线缆作为起点,利用抽水泵对选定的不同刻度依次加热,观察温度分布曲线显示结果.而后,依据本文第3 节所给出的模型进行理论值计算,得出的结果与实际膨胀量作对比,验证方法的可行性与准确性.

3.2 实验数据分析

根据第4.1 节中实验思路进行实验,由柱形光纤信号采集器一端前1 m 绕制成环作为0 起始点,每隔5 圈选取一个位置进行模拟膨胀节移动,直至第75 圈处的温度变化情况,实验结果如表2 所示.

表中圈数代表绕至于柱形缠绕体的光纤线缆圈数;光纤线缆长度表示实验中所用光纤线缆的实际长度;显示距离表示光纤线缆绕制的温度变化点于上位机界面的显示距离;实际膨胀量表示抽水泵在柱形光纤采集器上移动的实际距离;计算膨胀量表示利用第3.2 节所推出数学模型计算出的理论值.

3.3 误差分析

实验过程中膨胀量的计算值与实际值会存着一定的误差.对于不同等大小的被测量,需采用相对误差形式来表示测量的精确度.针对本实验,实际相对误差指的是膨胀量的绝对误差∆与被测量的约定真值(膨胀量实际值)A0之比,即

表2 实验数据表Table 2 Experimental data table

其中,x为膨胀量计算值.根据表2 中实验数据以及上述误差分析式(15),计算得到相应的相对误差列于表2 第6 列.

图5 为膨胀量测量相对误差分布图,给出的是以40◦C 时为例的光纤线缆绕制圈数与膨胀量测量相对误差之间的关系.在柱形绝热体直径一定的条件下,由于所绕制的圈数与光纤长度有对应线性关系,且读数便捷,因此给出了圈数与膨胀量测量误差之间的关系,即代表所用光纤线缆长度与膨胀量测量误差之间的关系.从图中可以看出,除第1 圈外,其他圈数点对应的膨胀节测量误差均在3% 以内.由于受系统空间分辨率的影响即本系统分辨率为1米,上位机界面显示的只有1 的整数倍.故本实验测得的如表2 实验数据中第1 圈实际光纤线缆长度785 毫米,与上位机界面显示距离1 000 毫米之差为215 毫米,造成了如图5 膨胀量测量相对误差分布图中第1 圈误差变化达27%.为此,在实际分析过程中,我们不考虑对第1 圈进行误差分析.

随后对表2 中相对误差值进行了数据拟合,得到图5 中膨胀量测量相对误差趋势曲线.由拟合曲线可以看出,相对误差的整体趋势随着圈数的增加而减小.此外,论文还对膨胀节膨胀量大小与光纤线缆圈数之间的关系进行了分析,综合表2 中数据以及图5 中膨胀量测量相对误差趋势曲线,可以得到以下分析结果:当柱形绝热体的直径为250 mm 时,膨胀量长度大于3.3 mm 时,利用本文中提到的方法检测膨胀节的膨胀量,测量结果均可以控制在3%以内;实验结果表明,随着光纤线缆圈数的增加即绕制光纤线缆长度的增加,膨胀量测量值与实际值之间的相对误差呈递减的趋势.

图5 膨胀量测量相对误差分布图Fig.5 Distribution of relative error in measurement of expansion

为进一步证实结果的规律性并增强实验的说服力,分别采用了直径为300 mm 与350 mm 两种型号的柱形绝缘体进行光纤线缆绕制.在同一温度下,依照第4.1 节实验中的步骤进行对比实验.得到以直径分别为250 mm、300 mm 以及350 mm 的柱形绝缘体作为载体绕制光纤线缆的情况下,如图6 所示的不同载体膨胀量测量相对误差分布图.

图6 不同载体膨胀量测量相对误差分布图Fig.6 Distribution diagram of the expansion relative error in different carriers

从图6 中可以得出以下分析结果:实验中采用三种不同型号的柱形绝缘体测得的膨胀量相对误差,均有随圈数增加而减小的趋势,再次印证了前述规律;三种不同直径的柱形绝缘体在绕制圈数大于10圈时,相同圈数下得到的膨胀量相对误差大小与柱形绝缘体的直径长短近似成反比;光纤线缆绕制圈数大于10 圈时,其测量误差均在3% 以内.这为膨胀量的高精度测量提供了一种新思路,即按要求的精度事先绕制一定圈数光纤线缆并以绕制结束点定为触头起点,而后运用本文提出的方法进行膨胀量测量从而达到更高的测量精度.

4 结论

论文根据拉曼分布式光纤测温系统的温度及位置的解调原理,综合热力管网中膨胀节膨胀量检测难点,提出了一种利用分布式光纤温度传感技术的膨胀量测量方法,很好地解决了热网状态实时检测中的关键问题,并通过膨胀量与上位机显示之间的关系推导出了膨胀量数学模型;实验结果表明,实际应用中,当柱形绝缘体的直径为确定值250 mm 且膨胀量≥3.3 mm 时,该模型的测量误差可控制在3% 以内,随着膨胀量的增大其测量误差呈近似线性减小趋势;不同柱形绝缘体直径下,膨胀量测量相对误差随其直径的增大而减小,且总能找到能够接受范围内的误差值.进一步证明了利用此方法的可行性与准确性.

基于拉曼分布式光纤温度监测系统的膨胀节膨胀量检测方法具有检测方便,读数便捷,以及实时检测的优点.对于热力管网中膨胀节部分的健康状态检测有很大的利用价值,能够有效准确地给出膨胀量随温度变化位置的关系,对危险情况起到提前预判的作用.