第31卷,第7期 2 0 1 1年7月 光谱学与光谱分析 Vo1.31,No.7,pp1862—1867 Spectroscopy and Spectral Analysis July,2011 光纤分布式扰动传感器中瑞利和受激布里渊 散射光谱特性及其影响研究 张春熹,梁 生,冯秀娟,林文台,李 琛,李 勤,钟 翔,李立京 北京航空航天大学仪器科学与光电工程学院,北京 100191 摘要针对光纤分布式扰动传感器长距离应用定位误差增大的问题,研究了传感器长距离光纤中的瑞利 散射和受激布里渊散射的光谱特性及其对系统定位精度的影响。建立了长距离光纤中瑞利和受激布里渊散 射的数学模型,通过仿真研究发现散射光波的杂波干涉和强度噪声会严重降低干涉仪的光学信噪比,从而 导致定位误差。通过实验验证了仿真结果的正确性。研究结论可以为提高长距离应用定位精度提供理论指 导。 关键词光纤分布式扰动传感器;马赫一泽德干涉仪;瑞利散射;受激布里渊散射;定位误差;光谱特性 文献标识码:A EOI:10.3964 ̄.issn.1000-0593(2011)07-1862—06 (C-MZI)和逆时针方向MZI(CC-MZI)。 中图分类号:TH744.1;TN253 引 言 光纤分布式扰动传感器(fiber-optic distributed disturb— ance sensor,FDDS) ̄.够对传感光纤上任意一点的扰动进行 在工程应用中,传感光纤n和b、传导光纤c和d通常布 置于同一光缆中。扰动信号_厂(£)作用于P点时,探测器PDl 和Plk接收的来自C-MZI和CC—MZI的干涉光强L和Jc 分 别为 探测和定位,可以应用于周界安防、油气管线预警、通信线 路监测以及结构监测[1 ]。长距离光纤中的散射是影响 FDDS长距离应用定位精度的主要因素。一些文献已经报道 了瑞利或受激布里渊散射对光纤放大器【3“]、光纤通信系 统_5]、光纤激光器[ 、波长转换器[ 、光纤射频信号系统[8]、 光纤气体传感器[9]、光纤陀螺[10,11]、光纤温度传感器 ]以及 其他干涉型光纤传感器【】 等光纤光学系统性能的影响,但 对于长距离FDDS中散射影响的研究工作还属空白。 本文针对FDDS长距离应用问题,研究了F【)I)S长距离 J (£)一百1 J {1+cos[Afo+B(1一K)f(t—Z'I)])(1) k(£)一÷J {1+cos[Ag ̄+B(1一K)f(t—r2)])(2) 其中, 为光源输出光强,△伽为MZI臂长差导致的初始相 位差,B(1一K) ( — )和B(1一K),(£一r2)为扰动导致的 相位调制,B为扰动相位调制灵敏度,K描述光纤n和6之 间扰动灵敏度的差异,c为真空中的光速, 为光纤纤芯折射 率,L 和L 分别是扰动点P距离C2和 的长度, 为 传感和传导光纤中的瑞利和受激布里渊散射对系统定位精度 的影响,从而为改善长距离定位精度提供理论指导。 c。至PDl间光纤的长度,时延 一 {生 和 一L.A t 。由 式(1)和(2)所示,干涉信号I 和J 为时域相关信号,存在 与扰动位置有关的时延r 和r2,则J 和Jcc之间的时延差 为 -L1 △r—r】一r2一—Lz-kLa-—1光纤分布式扰动传感器原理 FDDS的原理结构如图1所示,系统光路包括窄线宽激 光光源,2×2光纤耦合器(Coupler)C ,C2和C3,光电探测 一 (3) C C 其中,L=L +L2+L。=2L ,L I为监测距离。计算J 和 k的互相关函数R(J ,k),根据R(J ,I )的极值对应的时 延可以确定△r,通过下式可以实现对扰动的定位 L】一—L.器(Photo Detector,PD)PEh和PD2,传感光纤n,b以及传 导光纤C,d。系统光路是一个Mach-Zehnder干涉仪(MZI), 激光光波从两个方向同时入射,分别构成了顺时针方向MZI 收稿日期:2010—07—08。修订日期:2010一l1一O3 cAr/n (4) 基金项目:北京航空航天大学博士研究生创新基金项目资助 作者简介:梁生,1981年生,北京航空航天大学博士研究生 e-mail:zchunxi@hotmail.conl 第7期 光谱学与光谱分析 1863 …………………一一 Fig.1 Schematic illustration of FDDS and distributions of RB and DRB 传感光纤“和b巾,顺时针方向与PW f涉的RB的光强分 2瑞利散射影响 光纤中折射率分布的微观不均匀性会导致瑞利散射 j。 在FDDS长距离光纤中,需要考虑信号光波(primary wave, PW)、一阶瑞利后向散射光波(rayleigh backscattering,RB) 和二阶瑞利后向散射光波(double rayleigh backscattering, DRB)。 别为 J【 Rl{( )一,∞(£)d SE1一exp(——d- )]/2口 j ( Rl{(£)=f∞(z)口 sE1一exp(--0L )3/2a (7) (8) 其中,8—0.25 dB或0.058 Np・km‘ 为单模光纤衰减系 数l1 。, 为散射导致的光纤衰减系数,在单模光纤中,有:a ≈a 。S一1.5×10。为再新捕捉因子。 (f)为C。处的入射 信号光强,则 J∞(f)一 0(f)exp( … 1)/4 (9) 在顺时针和逆时针方向,分别有PW,RB和DRB到达 探测器PD1和P 。以Pr)l为例,到达PD 的光波电场复振 幅为 E。 一点 w+E +El籼+E Pw+E Rn+E r (5) 其中.凡(f)为激光光源输出光强。经过c{后,通过传导光 纤 传输L… 距离,最终到达PDl的RJ3光强为 JcRn(£)一生 口 Sexp(一2也 I)[1一exp(一aL )]/ (i0) 其中,Epw, 表示为 ,ERB,E R1{,Er,ea{和 IⅫ 分别为MZI两臂 光纤a和b的Pw,RB和DRB的复振幅,则1t涉光强可以 j 一E IE 一js+ (6) J 【、RH( )一 — a Sexp(一2 , 1)[1一exp(一口L )]/4d (11) 其中,,s和 分别表示=f:涉光强中的信号项和噪声项。信 号项为Pw之问的干涉光强,RB和DRB之问及其与PW之 间的干涉导致r噪声项J 。 2.1顺时针方向MZI中的瑞利散射 与RB不同的是,整个光路中DRB与相同传播方向的 Pw在光秤差小于L 的范围内均能发生干涉。为r计算 DRB j PW的于涉光强,将光纤分为N段,每段光纤长度 为L 一L/2,认为每段光纤中的DRB可以与Pw完全干涉, 小在同一段中的DRB与PW完全不干涉。将总长度为2L 的光纤分为N=2L 1/(L /2)=4L /L 段。第,2段光纤中 RB的光强为 在C-MZI中,只需要考虑PW和RB以及PW和DRB 之间的干涉光强。忽略PW,RB和DRB在单模光纤中偏振 态变化的影响_ 。如图1所示,则在光纤长度L大于激光光 1源相干长度L 的情况下,只有光程差为0的空间伯霞(C。) 前后(~L /2,L /2)范周内的RB与Pw之问发q三干涉。在 J K】 (")一Lo(”)d sE1——exp(——2。 )]/2口 (12) 其【1l, —J。exp[--a(n一1)/L ]/4为第 段光纤的入射信 1864 光谱学与光谱分析 第 段光纤a和b中RB的光强分别为 第3l卷 号光强。则第 段光纤中DRB的光强为 ImB(n)一 (n)口 s [1一exp(一2aL )] /(2a) (13) 到达探测器PDl处的DRB的光强分别为 f。RB(£)一 1o口 S2exp(一2aL,o ̄)E1一exp(一2aL )] /(2a)。 JKB_ ( )一 ( )口 s[1一exp(2aL )]/2a ItRB ̄t( )一J ( )a sE1一exp(2aL )]/2a 为 (18) (19) 第 段光纤a和b中最终经过 到达PD2的光强分别 (14) I'mB(t)一粤a:Szexp(-2aL ta1)[1一exp(-2aL )] /(2口) (15) add( )=o.5鲁口 sexp(一2anL )[1一exp(-2aL )]/2口 (20) 其中, (n)和 ̄OtDRB( )为随机变量,描述DRB光波相位 J RB_ ( )一0.5鲁口 sexp(-2onL )[1-exp(一2也 )]/2a (21) 的随机变化。最终,DRB与PW的干涉项可以表示为 2 N cos(卿一 )一 第n段光纤a和b中的RB之间发生干涉,最终导致的 噪声项为 2∑ ̄/—IcPwIC—DRB(n)cos[ 一 DRB( )](16a) n=I ecI 一2 ̄/ i N cos[仇x: s 一 2. f cos(卿一 c-砌)一 cc a ]一2∑ Io 3 dd(n)f CcR1{ (,z)・ cosEqa:c ̄_ dd(以)一 c【=RB_ ( )] (22) 2∑ cos[卿一 一(n)](16b) 一 一)一 2/7 『7 cos( 2E ̄/—I'cewI'—cmB(n)cos[ c 一 2 ̄/J∞RHJ c c0s(卿N ( )](16c) 3受激布里渊散射影响 光纤中入射光功率超过受激布里渊散射阈值功率P 时 会发生受激布里渊散射(stimulated brillouin scattering, SBS)_1 。 。在实际应用中,光路输出至PD的信号光功率不 能低于某~数值,即最小探测光功率 小值 可以表示为 Io一2J exp(2缸,。 ) (23) 一 ( )一 2 :;j cos[ ( )一 ](16d) 2.2逆时针方向MZI中的瑞利散射 。I…由电路系统的 在CC-MZI中,与PW发生干涉的RB的分布与C-MZI 相似,只是RB与PW发生干涉的有效区域长度变为L (如 图1)。CC-MZI中与PW发生干涉的DRB分布与C-MZI相 同,这里不再赘述。 信噪比决定,通常约为1~5 w。则激光光源输出光功率最 相应地,传导光纤C和传感光纤a和b中入射光功率的 最小值可以分别表示为 Jf in—Io/2一J exp(2aL ta1) Io~一j :lo/4=Iminexp(2aL。。 /2) (24a) (24b) 与C-MZI不同的是,在CC-MZI中,传感光纤a和b中 RB之间的干涉光强不能忽略。图1中C2至L tal—L /2处, 光纤a和b中的RB之间发生干涉。将传感光纤n和b从C2 至L l—L /2长度分为~段,每段长度为L 一L ,则N一 (L 一L /2)/L 。第 段光纤中的输入光强为 ( )一Ioexp[一a(n一1) I/4 CC.PW CC—Stokes 假设 一2“W,Iamin,h rrifn和Lmln随监测距离增加增 大,而P 则随监测距离增加减小,当监测距离分别超过约 (17) 68和65 krn时,将会在传感光纤n和b以及传导光纤c中发 Fig.2 Distribution of stokeswaveinFDDS 第7期 光谱学与光谱分析 1865 生SBS。传导光纤f巾首先发生SBS时,CC-MZI的PW由 开始出现定位误差;同时考虑RB、DRB和SBS时,在监测 距离为65 km时,开始出现定位误差,最大定位误差约为2O 于SBS导致严重衰减,此时为了保证探测器接收PW的最小 光功率,需要提高入射光功率以补偿SBS导致的衰减。补偿 后的入射功率为:f 。一10/(1一 )。其中,a 是SBS导致的 近似线性衰减系数。 km。为了描述MZI干涉信号中信号与噪声幅值的关系,定 义光学信噪[:t(optical signal—noise ratio,OSNR)为 ()SNR一101g(js/iN) (28) 补偿入射光功率后,导致监测距离为65 km时传感光纤 a和b巾发生SBS。如图2所示,传感光纤a和b发生SBS时 通过计算C—MZI和CC—MZI的OSNR,可以发现,CC- MZI中光纤n和b中RB之问的1二涉和SBS发生时Stokes光 产生的Stokes光波沿逆时针方向经过 进入PI)!。PDl和 PD2接收到的信号光波时戒平均光强分别为 r, <J )一 i 0exp(--2aL一口B) (25a) ‘ r, <IpD2)一 i o exp(--2aL一口B)(25b) 6 考虑到实际入射光功率比较接近P 小于SBS饱和光 功率,SBS产生的Stokes光波光强的时间平均值可以通过近 似的线性模型得到 < )一(f。一Plh) I{ (26) 其中, 为发生SBS光纤的入射光强,k 为近似线性斜率, 可以通过实验标定得到。传感光纤a和b中的Stokes光波在 发牛F涉。受激布 渊散射中的Stokes光波的线宽为2O ~100 MHz 其相干性远低于PW,在光纤n和方中,忽略 Stokes光波与PW之问的光学拍作用,只考虑光强叠加。 SBS发生时,进入PI=)2的光强』z(£)可以表示为 I2(£)一j (£)+I(r Rt (£)+Icy】)Klj(t)+flnI( )+jl}A(,) 斗 (f)+2yI (r”)v/fHA(£)f B(£)cos[△m (f)](27) 其中, P、v(f), .RB(f),J x、啪j(f)和,ln (£)分别为逆时针 方向的Pw光强,RB光强,DRB光强及以上各项之问的于 涉光强。“ (£)和h ( )分别为传感光纤n和b中的Stokes 光强,M ( )为Stokes光波的相干函数,rn为光纤n和6中 Stokes光波之间的时延差,由△L决定,△ (f)为Stokes光 波之间的相位差。Stokes光波的光强 (f)和I (f)的时间 平均值叮通过式(26)得到,其瞬时值需要考虑Stokes光波的 强度噪声]1 7,18 。自发布里渊散射的随机性会导致Stokes光 波光强随机起伏。因此,SBS对系统的影响体现在以下三个 方面:(1)发生SBS时,会导致PW严重衰减;(2)在传感光 纤“和b rI1 SBS产q二的Stokes光波之问会发生杂波干涉; (3)Stokes光波光强随时间起伏,即强度噪声。 4数值仿真与讨论 通过Matlab对F【)DS长距离光纤中的RB、DRB以及 SBS导致定位误差的机理进行数值仿真研究。数值仿真涉及 的参数如下:激光光源线宽Av=25 kHz,I~一2 w,△L一 5 m,Stokes光波线宽△ {一35 MHz。 由于噪声作用的统计特性,在一定条件下,定位误差发 生随机波动。通常取同一条件下多个定位误差的最大值,即 最大定位误差作为评价系统定位精度的技术指标。监测距离 在5O~100 km范围内变化时最大定位误差的数值仿真结果 如图3所示,其中,对应每一个监测距离重复计算100次。 在图3 If】,仅考虑RB和DRB时,在监测距离为95 km时, 波之问的_f涉及其强度噪声是导致定位误差的关键因素。 RB之 Ⅱ的f涉以及Stokes光波之问的干涉及其强度噪声导 致CC-MZI的OSNR降低,当OSNR小于约一10 dB时,定 位算法中的互相关函数出现多个随机峰值,导致最大定位误 差约为2O km。 l 20 ;15 = l 10 5 主05 O 6O 7O 8O 90 l00 Monitored length/kin Fig.3 Simulated relationship between maximal location errors and monitored length 5实验与讨论 为了验征数值仿真结论的正确性,通过FDDS样机进行 实验验证。分别对30 ̄75 km监测距离情况进行定位误差 测试,相同监测距离测试100次,最大定位误差与监测距离 的关系见图4。其中,监测距离小于65 km时,最大定位误 差小于100 111;监测距离大于65 km时,最大定位误差约为 30 km。图5给出_『监测距离分别为6O km和75 km时02 MZI的散射光波光谱,可以发现,监测距离为75 km时,同 时发生受激布里渊散射和瑞利散射;监测距离为60 km时, 仅有瑞利散射发生。实验与理论分析及数值仿真基本一致。 其中,监测距离小于65 km时产生的定位误差可认为主要来 自电路噪声和偏振衰落的影响。 3O 2O lO 6结论 本文对FDDS长距离光纤中瑞利散射和受激布里渊散射 们 0毫。1866 光谱学与光谱分析 O 第31卷 日p ≥od日。 一目日oZ 一 一 一一 一 一 一 —lO 0加 如∞∞∞加 -20 30 --40 50 -1 550 1 550.2 1 550.4 1 550.6 l 550.8 1 55l Wavelength/nm 1 550 1 550.2 I 550.4 1 550.6 l 550.8 l 55I Wavelength/nm Fig.5 Spectrums of scattering 对系统定位精度的影响进行了理论和实验研究。仿真研究发 现传感光纤中的瑞利散射和受激布里渊散射导致CC_MZI的 通过实验验证了理论分析和数值仿真结论的正确性。本文的 研究结论可以为提高FDDS长距离应用的定位精度提供理论 指导。 OSNR降低,当OSNR小于约--10 dB时,会导致定位误差。 References m— /】u≥o口勺 N 时昌l0Z [1] SUN Qi—zhen,LIU De-ruing,WANG Jian(孙琪真,刘德明,王健).Acta Physica Sinica(物理学报),2007,56(10):5903 . ang Sheng,Zhang Chunxi,Lin Wentai,et a1.Optics Letters,2009,34(12):1858[2] Li[3] Zhang W,Peng J,Liu X,et a1.Optics Communications,2001,199:231. [4] Willems F W,van der Plaats J C,DiGiovanni D J.Electronics Letters,1994,30(8):645 Guo-qiang,Guo Pan,et a1.Proc.of 2009 International Forum on Information Technology and Applications.209 [5] Wang Yan-hong,Nioche M,Gilles H,Girard S Optics Communications,2008,281:2243. [6] Lar Hasegawa,Hiromu Toba.IEEE Photonics Technology Letters,1995,7(3):327. [7] Kyo Inoue,Takasi[8] Zhang S I ,O’Reilly J J.IEEE Photonics Technology Letters,1997,9(3):395. [9] Jin W,Demokan M S,Stewart G,et a1.Proc.of SPIE,1997,3105:266. [1o] Bums W K,Moeller R P Journal of Lightwave Technology,l983,LT_1(2):381. ntosh,Brian Culshaw.Journal of Lightwave Technology,1989,7(9):1323. [11] James M Macki Science and Technology,2006,17:1065. [12] K De Souza.Measurementidge A,Kirkendall C K.IEEE Photonics Technology Letters,2003,15(11):1582. [13] Cranch G A,Dandfey.IEEE Transactions on Co[14] Peter Healmmunications,1987,Com-35(2):210. [15] M Oskar van Deventer.Journal of Lightwave Technology,1993,11(12):1895. exander L Gaeta,Robert W Boy&[16] AlPhysical Review A,1991,44(5):3205. tz N I,Chraplyvy A R,Tkach R W,et a1.IEEE Photonics Technology Letters,1997,9(1):124. [17] Horowi[18] Oiupsjobacka,Jacobsen G,Tromborg B Journal of I ightwave Technology,2000,18(3):416. Spectra and Influences of Rayleigh and Stimulated Brillouin Scattering in Fiber—Optic Distributed Disturbance Sensor ZHANG Chun-xi,LIANG Sheng,FENG Xiu-juan,LIN Wen-tai,LI Chen,LI Qin,ZHONG Xiang,LI Li-jing School of Instrument ciSence and Opto-Electronics Engineering,Beihang University,Beijing 100191,China Abstract Spectra and influences of Rayleigh backscattering(RB)and stimulated Brillouin scattering(SBS)in the fiber-optic distributed disturbance sensor(FDDS)were investigated.Models of RB,double RB(DRB)and SBS in long fibers were estab— lished.By numerical simulation,it was found that optical signal-to-noise ratio is extraordinarily reduced due to SBS and RB, which results in location errors.Numerical results were confirmed by experiments and helpful to improving the location precision for applications with long monitored length. 第7期 光谱学与光谱分析 1867 国际光谱会议预告 January 2011 17—21 IFPAC 2010 25th International Forum Process Analytical Technology。Baltimore,MD Contact:E-mail:info@ifpacnet.org Wehsite http f f .ifpacpat.org/ March 2011 13一l8 Pittcon 2011.Atlanta.GA Contact:E—mail:info@pittcon.org Websitet http |f 7、 .pittcon.org/ 27 3l 241st ACS National Meeting&Exposition.Anaheim,CA ontCact:E-mail:help@acs.org Website:http://porta1.acs.org/ May 2011 7—13 International Society for Magnetic Resonance in Medicine:19th Scientific Meeting& Exhibition,Montreal,Quebec, Canada Contact:ISMRM/SMRT Central Office, 2030 Addison Street,7th Floor,Berkeley,CA 94704,Ph:510 841 1899 Fax:510 84l 2340 Website http |{ .ismrm.org/meetings/index.html;E-mail info@ismrm.org 13—2O 15th International Conference on Near Infrared Spectroscopy(ICNRS).Cape Town.South Africa Contact:Conference Chair Marena Manley,Stellenbosch University,at:marena@nir2011.org; Conference Organizer Deidre Cloete,Conferences et al,at deidre@nir2011.org Ph+27(0)2l 963 22l5 Fax:+27(0)21 963 2215 Wehsite:http.//wwW.nit2011.org/