CN111947696B - 矢量光学啁啾链布里渊时域分析系统 - Google Patents

矢量光学啁啾链布里渊时域分析系统 Download PDF

Info

Publication number
CN111947696B
CN111947696B CN202010817451.XA CN202010817451A CN111947696B CN 111947696 B CN111947696 B CN 111947696B CN 202010817451 A CN202010817451 A CN 202010817451A CN 111947696 B CN111947696 B CN 111947696B
Authority
CN
China
Prior art keywords
brillouin
vector
optical
electro
polarization
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202010817451.XA
Other languages
English (en)
Other versions
CN111947696A (zh
Inventor
张敬栋
朱涛
郑华
吴昊庭
郭南
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Chongqing University
Original Assignee
Chongqing University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Chongqing University filed Critical Chongqing University
Priority to CN202010817451.XA priority Critical patent/CN111947696B/zh
Publication of CN111947696A publication Critical patent/CN111947696A/zh
Application granted granted Critical
Publication of CN111947696B publication Critical patent/CN111947696B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01DMEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
    • G01D5/00Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable
    • G01D5/26Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light
    • G01D5/32Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light
    • G01D5/34Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light the beams of light being detected by photocells
    • G01D5/344Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light the beams of light being detected by photocells using polarisation

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Optical Transform (AREA)

Abstract

本发明提供一种矢量光学啁啾链布里渊时域分析系统,包括激光器,激光器与耦合器连接,耦合器的两个输出端分别连接第一电光调制器、第二电光调制器的第一输入端,任意波形发生器的两个输出端分别连接第一电光调制器、第二电光调制器的第二输入端,第一电光调制器的输出端连接第一密集波分复用器的输入端,第一密集波分复用器的两个输出端分别通过第一偏振控制器、第二偏振控制器连接偏振合束器的输入端,偏振合束器的输出端连接传感光纤的第一端;第二电光调制器的输出端连接环形器的第一端,环形器的第二端连接传感光纤的第二端且第三端连接第二密集波分复用器的输入端,第二密集波分复用器的两个输出端分别连接第一光电探测器和第二光电探测器。

Description

矢量光学啁啾链布里渊时域分析系统
技术领域
本发明属于光纤传感领域,具体涉及一种矢量光学啁啾链布里渊时域分析系统。
背景技术
基于布里渊散射的分布式光纤传感技术可在长距离传感光纤上以较高的灵敏度和准确性对温度或应变进行分布式测量。这些系统具有诸如低成本,抗腐蚀,耐用性和抗电磁干扰性等优点,这使其成为用于电力线,桥梁以及石油和天然气管道等大型结构的在线监测工具。
其中布里渊光学时域分析系统是一项成熟的技术,其原理是探测光通过受激布里渊散射非线性效应与反向传播的泵浦脉冲光相互作用。通过逐步扫描布里渊频移附近的泵浦探测频率偏移,重建了被测光纤的布里渊频谱分布。快速的布里渊采集方法避免了费时的频率扫描和轨迹平均,因此也引起了广泛的研究兴趣。其中新型的基于光学啁啾链探测光的布里渊光学时域分析系统,通过光学啁啾链的瞬时布里渊响应对光纤的布里渊频移进行提取,具有动态探测性能好的优点。但是该系统仅利用了布里渊瞬态响应的强度部分,需要较多的数据平均才能进行信号解调,存在传感精度低,信号弱的缺点。
发明内容
本发明提供一种矢量光学啁啾链布里渊时域分析系统,以解决目前进行光学布里渊时域分析时采集数据次数高、传感精度低且信号弱的问题。
根据本发明实施例的第一方面,提供一种矢量光学啁啾链布里渊时域分析系统,包括激光器,该激光器的输出端与耦合器的输入端连接,该耦合器的第一输出端与第一电光调制器的第一输入端连接,该第一电光调制器的第二输入端与任意波形发生器的第一输出端连接且其输出端与第一密集波分复用器的输入端连接,该第一密集波分复用器的第一输出端和第二输出端分别通过第一偏振控制器、第二偏振控制器与偏振合束器的输入端连接,该偏振合束器的输出端连接传感光纤的第一端;
该耦合器的第二输出端与第二电光调制器的第一输入端连接,该第二电光调制器的第二输入端与该任意波形发生器的第二输出端连接且其输出端与环形器的第一端连接,该环形器的第二端连接该传感光纤的第二端且其第三端连接第二密集波分复用器的输入端,该第二密集波分复用器的第一输出端和第二输出端分别连接第一光电探测器和第二光电探测器;
所述激光器输出的激光被所述耦合器分成两路激光,其中一路激光在所述第一电光调制器的调制下产生两个边带,所述两个边带被该第一密集波分复用器分成上边带和下边带,所述第一偏振控制器和第二偏振控制器分别对上边带和下边带进行偏振控制,以使上边带和下边带的偏振态正交,偏振态正交的上边带和下边带经该偏振合束器合束后传输给该传感光纤的第一端,其中该上边带和下边带中均包括啁啾链和参考光;
另一路激光被所述第二电光调制器调制成脉冲光信号,该脉冲光信号通过该环形器被传输给该传感光纤的第二端,在该传感光纤中上边带和下边带中的啁啾链分别与该脉冲光信号发生受激布里渊散射,分别产生上边带的探测光和下边带的探测光,该上边带中的探测光和参考光、以及该下边带的探测光和参考光通过该环形器被传输给第二密集波分复用器,该第二密集波分复用器将该上边带中的探测光和参考光、以及该下边带的探测光和参考光分别输出,由该第一光电探测器对该上边带中的探测光和参考光进行拍频,获得第一拍频信号,由第二光电探测器对该下边带中的探测光和参考光进行拍频,获得第二拍频信号;对该第一拍频信号和第二拍频信号进行数字相干解调,分别获得上边带和下边带中探测光的幅相信息,处理器将上边带的幅相信息与下边带的幅相信息相减,消除啁啾链与脉冲光信号偏振变化引起的偏振衰落现象,获得无偏振衰减的探测光幅相信息,基于该无偏振衰减的探测光幅相信息,进行矢量光学布里渊时域分析。
在一种可选的实现方式中,在该偏振合束器的输出端与该传感器的第一端之间还设置有隔离器。
在另一种可选的实现方式中,在该第二电光调制器的输出端与该环形器的第一端之间还设置有光放大器。
在另一种可选的实现方式中,所述任意波形发生器输出的波形满足:
Figure BDA0002633236510000031
其中,V0是输出信号的幅值,T为一个啁啾脉冲的宽度,K啁啾脉冲的啁啾率,fb为光纤的布里渊频移,fR为参考光与探测光之间的频差,t表示时间。
在另一种可选的实现方式中,所述第一电光调制器输出的光场表示为:
Figure BDA0002633236510000032
其中,E0是电光调制器输入光场的复振幅,fc为光载波的频率,T为一个啁啾脉冲的宽度,K啁啾脉冲的啁啾率,fb为光纤的布里渊频移,fR为参考光与探测光之间的频差,t表示时间。
在另一种可选的实现方式中,基于该无偏振衰减的探测光幅相信息,进行矢量光学布里渊时域分析时,按照以下步骤进行:
步骤S110、将该无偏振衰减的探测光幅相信息转换为矢量布里渊谱
Figure BDA0002633236510000033
Figure BDA0002633236510000034
其中A(f)表示该传感光纤上每段光纤的布里渊幅度谱,
Figure BDA0002633236510000035
表示该传感光纤上每段光纤的布里渊相移谱;
步骤S120、根据该矢量布里渊谱
Figure BDA0002633236510000036
采用复数主成分分析方法对布里渊频移信息进行提取。
在另一种可选的实现方式中,所述步骤S120包括:
步骤S121、通过数值仿真方法改变光纤布里渊中心频移,获得不同的矢量布里渊谱,将不同的矢量布里渊谱组成M×N的参考矩阵XM×N,其中M为布里渊中心频移个数,N为一个矢量布里渊谱的数据点数;
步骤S122、对参考矩阵XM×N进行归一化形成零均值参考矩阵
Figure BDA0002633236510000041
步骤S123、计算零均值参考矩阵
Figure BDA0002633236510000042
的协方差矩阵:
Figure BDA0002633236510000043
其中
Figure BDA0002633236510000044
Figure BDA0002633236510000045
的共轭转置;
步骤S124、对CN×N进行特征分解:
Figure BDA0002633236510000046
其中ΛN×N=diag(λ1,...,λN)为实数对角阵,包含了由大到小的特征值λj;UN×N为复数方针,其列向量uj为特征值λj的特征向量,又叫做参考矩阵的主成分,λj代表了主成分uj的能量;
步骤S124、将前j个主成分的积累能量定义为:
Figure BDA0002633236510000047
选择远小于N的数字L,使得其积累能量gL大于设定阈值α;
步骤S125、将复数方阵UN×N的前L列组成复数转换矩阵WN×L
步骤S126、零均值参考矩阵
Figure BDA0002633236510000048
在复数转换矩阵WN×L的投影为:
Figure BDA0002633236510000049
其中SM×L的行向量si为零均值参考矩阵
Figure BDA00026332365100000410
相应行xi的分数向量,对于待测的矢量布里渊谱向量y1×N,其在复数转换矩阵WN×L的投影为分数向量z1×L=y1×NWN×L
步骤S127、为了获得待测的矢量布里渊谱向量y的布里渊频移,对分数向量z分别和向量si求复数向量空间欧氏距离:
Figure BDA00026332365100000411
其中||z-si||2为z-si的欧式范数;
步骤S128、将欧式距离最小的第m个分数向量sm对应的布里渊中心频移作为待测矢量布里渊谱向量的布里渊频移。
本发明的有益效果是:
1、本发明不仅可以解调得到布里渊增益谱,还可以得到布里渊相移谱,实现了矢量布里渊分布式传感;本发明使用相干探测的结构,有效调高了信噪比,提高了系统传感精度;本发明同时使用布里渊增益谱和相移谱,进一步提高了系统的传感精度;本发明采用了偏振分集技术,有效消除系统中的偏振衰落现象,提高了系统的传感速度;本发明采用大带宽的光学啁啾链,具有更高的光学非线性阈值,因此啁啾链和脉冲光信号的能量更高,提高了信噪比,同时系统具有更大的动态应变响应范围;
2、本发明提出了复数主成分分析方法对测得的矢量布里渊谱进行布里渊频移解调,可以大幅提高光学啁啾链布里渊光时域分析系统的传感精度。
附图说明
图1是本发明矢量光学啁啾链布里渊时域分析系统的一个实施例结构示意图;
图2是本发明工作原理图;
图3是幅度和相移迹线示意图,其中图3(a)中两条曲线分别表示上下边带的幅度解调迹线,图3(b)中两条曲线分别表示上下边带的相移解调迹线;图3(c)表示上下边带的幅度解调迹线相减后的幅度迹线,图3(d)表示上下边带的相移解调迹线相减后的相移迹线;
图4是矢量布里渊的幅度响应谱和相移响应谱;
图5是本发明矢量光学啁啾链布里渊时域分析系统的另一实施例结构示意图。
具体实施方式
为了使本技术领域的人员更好地理解本发明实施例中的技术方案,并使本发明实施例的上述目的、特征和优点能够更加明显易懂,下面结合附图对本发明实施例中技术方案作进一步详细的说明。
在本发明的描述中,除非另有规定和限定,需要说明的是,术语“连接”应做广义理解,例如,可以是机械连接或电连接,也可以是两个元件内部的连通,可以是直接相连,也可以通过中间媒介间接相连,对于本领域的普通技术人员而言,可以根据具体情况理解上述术语的具体含义。
参见图1所示,为本发明矢量光学啁啾链布里渊时域分析系统的一个实施例结构示意图。该矢量光学啁啾链布里渊时域分析系统可以包括激光器,该激光器的输出端与耦合器的输入端连接,该耦合器的第一输出端与第一电光调制器的第一输入端连接,该第一电光调制器的第二输入端与任意波形发生器的第一输出端连接且其输出端与第一密集波分复用器的输入端连接,该第一密集波分复用器的第一输出端和第二输出端分别通过第一偏振控制器、第二偏振控制器与偏振合束器的输入端连接,该偏振合束器的输出端连接传感光纤的第一端;该耦合器的第二输出端与第二电光调制器的第一输入端连接,该第二电光调制器的第二输入端与该任意波形发生器的第二输出端连接且其输出端与环形器的第一端连接,该环形器的第二端连接该传感光纤的第二端且其第三端连接第二密集波分复用器的输入端,该第二密集波分复用器的第一输出端和第二输出端分别连接第一光电探测器和第二光电探测器。
所述激光器输出的激光被所述耦合器分成两路,其中一路激光在所述第一电光调制器的调制下产生两个边带,所述两个边带被该第一密集波分复用器分成上边带和下边带,所述第一偏振控制器和第二偏振控制器分别对上边带和下边带进行偏振控制,以使上边带和下边带的偏振态正交,偏振态正交的上边带和下边带经该偏振合束器合束后传输给该传感光纤的第一端,其中该上边带和下边带中均包括啁啾链和参考光;另一路激光被所述第二电光调制器调制成脉冲光信号,该脉冲光信号通过该环形器被传输给该传感光纤的第二端,在该传感光纤中上边带和下边带中的啁啾链分别与该脉冲光信号发生受激布里渊散射,分别产生上边带的探测光和下边带的探测光,该上边带中的探测光和参考光、以及该下边带的探测光和参考光通过该环形器被传输给第二密集波分复用器,该第二密集波分复用器将该上边带中的探测光和参考光、以及该下边带的探测光和参考光分别输出,由该第一光电探测器对该上边带中的探测光和参考光进行拍频,获得第一拍频信号,由第二光电探测器对该下边带中的探测光和参考光进行拍频,获得第二拍频信号;对该第一拍频信号和第二拍频信号进行数字相干解调,分别获得上边带和下边带中探测光的幅相信息,处理器将上边带的幅相信息与下边带的幅相信息相减,消除啁啾链与脉冲光信号偏振变化引起的偏振衰落现象,获得无偏振衰减的探测光幅相信息,基于该无偏振衰减的探测光幅相信息,进行矢量光学布里渊时域分析。
本实施例中,该激光器可以为窄线宽激光器,例如分布反馈半导体激光器,为了形成基于相干探测的光学啁啾链布里渊光时域分析系统,所述任意波形发生器输出的波形满足:
Figure BDA0002633236510000071
其中,V0是输出信号的幅值,T为一个啁啾脉冲的宽度,K啁啾脉冲的啁啾率,fb为光纤的布里渊频移,fR为参考光与探测光之间的频差,t表示时间。将啁啾脉冲首尾相接即可组成啁啾脉冲链。由于本发明采用的是双边带调制,因此所述第一电光调制器输出的光场表示为:
Figure BDA0002633236510000072
其中,E0是电光调制器输入光场的复振幅,fc为光载波的频率,T为一个啁啾脉冲的宽度,K啁啾脉冲的啁啾率,fb为光纤的布里渊频移,fR为参考光与探测光之间的频差,t表示时间。假设调制电压很小,则只有一阶边带保留,高阶边带可以忽略,如图2所示,本发明只存在一阶边带。为了消除单模光纤中的偏振衰落现象,本发明采用偏振分集技术,通过第一密集波分复用器将上边带和下边带分别传输给第一偏振控制器和第二偏振控制器进行偏振控制,使得上边带和下边带的偏振态正交,偏振态正交的上边带和下边带被传输至传感光纤后,上边带中的啁啾链与脉冲光信号发生受激布里渊散射,在此过程中上边带的啁啾链产生损耗,下边带中的啁啾链与脉冲光信号也发生受激布里渊散射,但在此过程中下边带的啁啾链产生增益,结合图2所示。在该传感光纤中上边带和下边带中的啁啾链分别与该脉冲光信号发生受激布里渊散射,分别产生上边带的探测光和下边带的探测光,该第一光电探测器对该上边带中的探测光和参考光进行拍频,获得第一拍频信号,对该第一拍频信号和第二拍频信号进行数字相干解调,获得上边带和下边带的幅度解调迹线,如图3(a)所示,以及上边带和下边带的相移解调迹线,如图3(b)所示,将上边带和下边带的幅度解调迹线相减,得到无偏振衰减的幅度迹线,如图3(c)所示,将上边带和下边带的相移解调迹线相减,得到无偏振衰减的相移迹线,如图3(d)所示。
其中,本发明基于该无偏振衰减的探测光幅相信息,进行矢量光学布里渊时域分析时,可以按照以下步骤进行:
步骤S110、将该无偏振衰减的探测光幅相信息转换为矢量布里渊谱
Figure BDA0002633236510000081
Figure BDA0002633236510000082
其中A(f)表示该传感光纤上每段光纤的布里渊幅度谱,
Figure BDA0002633236510000083
表示该传感光纤上每段光纤的布里渊相移谱。其中,根据啁啾链时间和瞬时频率的对应关系,可以将无偏振衰减的探测光幅度和相移响应时域迹线变换为矢量布里渊时域分析的布里渊幅度和相移迹线,如图4所示,其中图4(a)是布里渊幅度响应谱,图4(b)是布里渊相移响应谱。
步骤S120、根据该矢量布里渊谱
Figure BDA0002633236510000084
采用复数主成分分析方法对布里渊频移信息进行提取。所述步骤S120包括:
步骤S121、通过数值仿真方法改变光纤布里渊中心频移,获得不同的矢量布里渊谱,将不同的矢量布里渊谱组成M×N的参考矩阵XM×N,其中M为布里渊中心频移个数,N为一个矢量布里渊谱的数据点数;
步骤S122、对参考矩阵XM×N进行归一化形成零均值参考矩阵
Figure BDA0002633236510000085
步骤S123、计算零均值参考矩阵
Figure BDA0002633236510000086
的协方差矩阵:
Figure BDA0002633236510000087
其中
Figure BDA0002633236510000088
Figure BDA0002633236510000089
的共轭转置;
步骤S124、对CN×N进行特征分解:
Figure BDA0002633236510000091
其中ΛN×N=diag(λ1,...,λN)为实数对角阵,包含了由大到小的特征值λj;UN×N为复数方针,其列向量uj为特征值λj的特征向量,又叫做参考矩阵的主成分,λj代表了主成分uj的能量;
步骤S124、将前j个主成分的积累能量定义为:
Figure BDA0002633236510000092
选择远小于N的数字L,使得其积累能量gL大于设定阈值α;
步骤S125、将复数方阵UN×N的前L列组成复数转换矩阵WN×L
步骤S126、零均值参考矩阵
Figure BDA0002633236510000093
在复数转换矩阵WN×L的投影为:
Figure BDA0002633236510000094
其中SM×L的行向量si为零均值参考矩阵
Figure BDA0002633236510000095
相应行xi的分数向量,对于待测的矢量布里渊谱向量y1×N,其在复数转换矩阵WN×L的投影为分数向量z1×L=y1×NWN×L
步骤S127、为了获得待测的矢量布里渊谱向量y的布里渊频移,对分数向量z分别和向量si求复数向量空间欧氏距离:
Figure BDA0002633236510000096
其中||z-si||2为z-si的欧式范数;
步骤S128、将欧式距离最小的第m个分数向量sm对应的布里渊中心频移作为待测矢量布里渊谱向量的布里渊频移。本发明采用复数主成分分析方法,可以获得待测传感光纤的布里渊频移随距离的分布曲线,完成矢量啁啾时域分析系统的解调。
由上述实施例可见,本发明不仅可以解调得到布里渊增益谱,还可以得到布里渊相移谱,实现了矢量布里渊分布式传感;本发明使用相干探测的结构,有效调高了信噪比,提高了系统传感精度;本发明同时使用布里渊增益谱和相移谱,进一步提高了系统的传感精度;本发明采用了偏振分集技术,有效消除系统中的偏振衰落现象,提高了系统的传感速度;本发明采用大带宽的光学啁啾链,具有更高的光学非线性阈值,因此啁啾链和脉冲光信号的能量更高,提高了信噪比。此外,本发明提出了复数主成分分析方法对测得的矢量布里渊谱进行布里渊频移解调,可以大幅提高光学啁啾链布里渊光时域分析系统的传感精度和动态响应能力。
参见图5,为本发明矢量光学啁啾链布里渊时域分析系统的另一实施例结构示意图。图5与图1所示矢量光学啁啾链布里渊时域分析系统的区别在于,在该偏振合束器的输出端与该传感器的第一端之间还设置有隔离器,在该第二电光调制器的输出端与该环形器的第一端之间还设置有光放大器,该光放大器对第二电光调制器输出的脉冲光信号进行放大。
本领域技术人员在考虑说明书及实践这里公开的发明后,将容易想到本发明的其它实施方案。本申请旨在涵盖本发明的任何变型、用途或者适应性变化,这些变型、用途或者适应性变化遵循本发明的一般性原理并包括本发明未公开的本技术领域中的公知常识或惯用技术手段。说明书和实施例仅被视为示例性的,本发明的真正范围和精神由下面的权利要求指出。
应当理解的是,本发明并不局限于上面已经描述并在附图中示出的精确结构,并且可以在不脱离其范围进行各种修改和改变。本发明的范围仅由所附的权利要求来管制。

Claims (7)

1.一种矢量光学啁啾链布里渊时域分析系统,其特征在于,包括激光器,该激光器的输出端与耦合器的输入端连接,该耦合器的第一输出端与第一电光调制器的第一输入端连接,该第一电光调制器的第二输入端与任意波形发生器的第一输出端连接且其输出端与第一密集波分复用器的输入端连接,该第一密集波分复用器的第一输出端和第二输出端分别通过第一偏振控制器、第二偏振控制器与偏振合束器的输入端连接,该偏振合束器的输出端连接传感光纤的第一端;
该耦合器的第二输出端与第二电光调制器的第一输入端连接,该第二电光调制器的第二输入端与该任意波形发生器的第二输出端连接且其输出端与环形器的第一端连接,该环形器的第二端连接该传感光纤的第二端且其第三端连接第二密集波分复用器的输入端,该第二密集波分复用器的第一输出端和第二输出端分别连接第一光电探测器和第二光电探测器;
所述激光器输出的激光被所述耦合器分成两路激光,其中一路激光在所述第一电光调制器的调制下产生两个边带,所述两个边带被该第一密集波分复用器分成上边带和下边带,所述第一偏振控制器和第二偏振控制器分别对上边带和下边带进行偏振控制,以使上边带和下边带的偏振态正交,偏振态正交的上边带和下边带经该偏振合束器合束后传输给该传感光纤的第一端,其中该上边带和下边带中均包括啁啾链和参考光;
另一路激光被所述第二电光调制器调制成脉冲光信号,该脉冲光信号通过该环形器被传输给该传感光纤的第二端,在该传感光纤中上边带和下边带中的啁啾链分别与该脉冲光信号发生受激布里渊散射,分别产生上边带的探测光和下边带的探测光,该上边带中的探测光和参考光、以及该下边带的探测光和参考光通过该环形器被传输给第二密集波分复用器,该第二密集波分复用器将该上边带中的探测光和参考光、以及该下边带的探测光和参考光分别输出,由该第一光电探测器对该上边带中的探测光和参考光进行拍频,获得第一拍频信号,由第二光电探测器对该下边带中的探测光和参考光进行拍频,获得第二拍频信号;对该第一拍频信号和第二拍频信号进行数字相干解调,分别获得上边带和下边带中探测光的幅相信息,处理器将上边带的幅相信息与下边带的幅相信息相减,消除啁啾链与脉冲光信号偏振变化引起的偏振衰落现象,获得无偏振衰减的探测光幅相信息,基于该无偏振衰减的探测光幅相信息,进行矢量光学布里渊时域分析。
2.根据权利要求1所述的矢量光学啁啾链布里渊时域分析系统,其特征在于,在该偏振合束器的输出端与该传感光纤的第一端之间还设置有隔离器。
3.根据权利要求1或2所述的矢量光学啁啾链布里渊时域分析系统,其特征在于,在该第二电光调制器的输出端与该环形器的第一端之间还设置有光放大器。
4.根据权利要求1所述的矢量光学啁啾链布里渊时域分析系统,其特征在于,所述任意波形发生器输出的波形满足:
Figure FDA0002633236500000021
其中,V0是输出信号的幅值,T为一个啁啾脉冲的宽度,K啁啾脉冲的啁啾率,fb为光纤的布里渊频移,fR为参考光与探测光之间的频差,t表示时间。
5.根据权利要求1所述的矢量光学啁啾链布里渊时域分析系统,其特征在于,所述第一电光调制器输出的光场表示为:
Figure FDA0002633236500000022
其中,E0是电光调制器输入光场的复振幅,fc为光载波的频率,T为一个啁啾脉冲的宽度,K啁啾脉冲的啁啾率,fb为光纤的布里渊频移,fR为参考光与探测光之间的频差,t表示时间。
6.根据权利要求1所述的矢量光学啁啾链布里渊时域分析系统,其特征在于,基于该无偏振衰减的探测光幅相信息,进行矢量光学布里渊时域分析时,按照以下步骤进行:
步骤S110、将该无偏振衰减的探测光幅相信息转换为矢量布里渊谱
Figure FDA0002633236500000031
Figure FDA0002633236500000032
其中A(f)表示该传感光纤上每段光纤的布里渊幅度谱,
Figure FDA0002633236500000033
表示该传感光纤上每段光纤的布里渊相移谱;
步骤S120、根据该矢量布里渊谱
Figure FDA0002633236500000034
采用复数主成分分析方法对布里渊频移信息进行提取。
7.根据权利要求6所述的矢量光学啁啾链布里渊时域分析系统,其特征在于,所述步骤S120包括:
步骤S121、通过数值仿真方法改变光纤布里渊中心频移,获得不同的矢量布里渊谱,将不同的矢量布里渊谱组成M×N的参考矩阵XM×N,其中M为布里渊中心频移个数,N为一个矢量布里渊谱的数据点数;
步骤S122、对参考矩阵XM×N进行归一化形成零均值参考矩阵
Figure FDA0002633236500000035
步骤S123、计算零均值参考矩阵
Figure FDA0002633236500000036
的协方差矩阵:
Figure FDA0002633236500000037
其中
Figure FDA0002633236500000038
Figure FDA0002633236500000039
的共轭转置;
步骤S124、对CN×N进行特征分解:
Figure FDA00026332365000000310
其中ΛN×N=diag(λ1,...,λN)为实数对角阵,包含了由大到小的特征值λj;UN×N为复数方阵 ,其列向量uj为特征值λj的特征向量,又叫做参考矩阵的主成分,λj代表了主成分uj的能量;
步骤S124、将前j个主成分的积累能量定义为:
Figure FDA00026332365000000311
选择远小于N的数字L,使得其积累能量gL大于设定阈值α;
步骤S125、将复数方阵UN×N的前L列组成复数转换矩阵WN×L
步骤S126、零均值参考矩阵
Figure FDA00026332365000000312
在复数转换矩阵WN×L的投影为:
Figure FDA00026332365000000313
其中SM×L的行向量si为零均值参考矩阵
Figure FDA0002633236500000041
相应行xi的分数向量,对于待测的矢量布里渊谱向量y1×N,其在复数转换矩阵WN×L的投影为分数向量z1×L=y1×NWN×L
步骤S127、为了获得待测的矢量布里渊谱向量y的布里渊频移,对分数向量z分别和向量si求复数向量空间欧氏距离:
Figure FDA0002633236500000042
其中||z-si||2为z-si的欧式范数;
步骤S128、将欧式距离最小的第m个分数向量sm对应的布里渊中心频移作为待测矢量布里渊谱向量的布里渊频移。
CN202010817451.XA 2020-08-14 2020-08-14 矢量光学啁啾链布里渊时域分析系统 Active CN111947696B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010817451.XA CN111947696B (zh) 2020-08-14 2020-08-14 矢量光学啁啾链布里渊时域分析系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010817451.XA CN111947696B (zh) 2020-08-14 2020-08-14 矢量光学啁啾链布里渊时域分析系统

Publications (2)

Publication Number Publication Date
CN111947696A CN111947696A (zh) 2020-11-17
CN111947696B true CN111947696B (zh) 2021-11-23

Family

ID=73343150

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010817451.XA Active CN111947696B (zh) 2020-08-14 2020-08-14 矢量光学啁啾链布里渊时域分析系统

Country Status (1)

Country Link
CN (1) CN111947696B (zh)

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100080564A1 (en) * 2008-09-26 2010-04-01 Oki Electric Industry Co., Ltd. Optical phase locked loop
CN102307061A (zh) * 2011-08-26 2012-01-04 上海交通大学 超短光纤中布里渊散射高精度测量系统
CN103743354A (zh) * 2014-01-06 2014-04-23 桂林电子科技大学 一种基于布里渊相移检测的动态应变测量方法及测量装置
CN104567960A (zh) * 2015-01-04 2015-04-29 西南交通大学 一种基于相位调制探测光的相干布里渊光时域分析传感系统
US20150377603A1 (en) * 2013-01-29 2015-12-31 Fudan University Method for reducing interference from scattered light/reflected light of interference path by generating carrier through phase
CN107764297A (zh) * 2017-10-19 2018-03-06 重庆大学 基于泵浦脉冲扫频的动态布里渊光时域分析系统
CN108827177A (zh) * 2018-07-05 2018-11-16 重庆大学 长距离分布式布里渊动态应变测量系统
CN108981768A (zh) * 2018-09-05 2018-12-11 哈尔滨工业大学 基于光学啁啾链的单端快速分布式布里渊光学时域反射计
CN109163748A (zh) * 2018-09-17 2019-01-08 哈尔滨工业大学 基于捷变频技术的单端动态分布式布里渊反射装置及方法
CN109724529A (zh) * 2019-01-04 2019-05-07 重庆大学 基于多斜坡辅助的大动态范围布里渊快速测量系统

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100080564A1 (en) * 2008-09-26 2010-04-01 Oki Electric Industry Co., Ltd. Optical phase locked loop
CN102307061A (zh) * 2011-08-26 2012-01-04 上海交通大学 超短光纤中布里渊散射高精度测量系统
US20150377603A1 (en) * 2013-01-29 2015-12-31 Fudan University Method for reducing interference from scattered light/reflected light of interference path by generating carrier through phase
CN103743354A (zh) * 2014-01-06 2014-04-23 桂林电子科技大学 一种基于布里渊相移检测的动态应变测量方法及测量装置
CN104567960A (zh) * 2015-01-04 2015-04-29 西南交通大学 一种基于相位调制探测光的相干布里渊光时域分析传感系统
CN107764297A (zh) * 2017-10-19 2018-03-06 重庆大学 基于泵浦脉冲扫频的动态布里渊光时域分析系统
CN108827177A (zh) * 2018-07-05 2018-11-16 重庆大学 长距离分布式布里渊动态应变测量系统
CN108981768A (zh) * 2018-09-05 2018-12-11 哈尔滨工业大学 基于光学啁啾链的单端快速分布式布里渊光学时域反射计
CN109163748A (zh) * 2018-09-17 2019-01-08 哈尔滨工业大学 基于捷变频技术的单端动态分布式布里渊反射装置及方法
CN109724529A (zh) * 2019-01-04 2019-05-07 重庆大学 基于多斜坡辅助的大动态范围布里渊快速测量系统

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
80 km Fading Free Phase-Sensitive Reflectometry;Jingdong Zhang 等;《JOURNAL OF LIGHTWAVE TECHNOLOGY》;20190915;第37卷(第18期);第4748-4754页 *
基于脉冲调制的分布式振动传感技术研究;张敬栋;《中国博士学位论文全文数据库 (基础科学辑)》;20190415;第1-142页 *
布里渊光时域分析动态应变传感技术研究进展;朱涛 等;《应用科学学报》;20200331;第38卷(第2期);第197-212页 *

Also Published As

Publication number Publication date
CN111947696A (zh) 2020-11-17

Similar Documents

Publication Publication Date Title
US10066973B2 (en) Brillouin scattering measurement method and brillouin scattering measurement system
US8988671B2 (en) BOTDA system that combined optical pulse coding techniques and coherent detection
CN110632388B (zh) 一种基于混频的光电探测器频响测量方法及装置
CN107036734B (zh) 一种全分布式光纤温度或应变的传感方法与传感器
US11725965B2 (en) Method and device for reconstructing a backscattered electromagnetic vector wave
CN109297581A (zh) 一种用于补偿相位敏感光时域反射计中频率漂移的二次相位差分测量方法
CN103414513B (zh) 一种具有高动态范围的脉冲光动态消光比测量装置及方法
CN113654679B (zh) 一种分布式光纤温度与应变同时传感系统
CN112033568B (zh) 一种双脉冲调制的温度与应变光纤传感系统
CN114543973B (zh) 一种分布式超高频振动信号测量方法及光纤传感器
Li et al. Pulse coding in distributed optical fiber vibration sensor: A review
CN109724529B (zh) 基于多斜坡辅助的大动态范围布里渊快速测量系统
CN113340571A (zh) 基于光矢量分析的光时延测量方法及装置
CN111637910B (zh) 时域差分高速混沌布里渊光相干域监测装置及方法
CN111947696B (zh) 矢量光学啁啾链布里渊时域分析系统
CN113686424A (zh) 基于波长分集技术的高信噪比声传感系统及多波长合并方法
CN114353685B (zh) 基于混沌布里渊相位谱的高频动态应变测量装置和方法
CN213842395U (zh) 动静态联合测量分布式光纤传感系统
CN111609919B (zh) 光纤分布式振动和损耗同时检测系统
CN111721438B (zh) 一种噪声调制线阵ccd采集的免扫频botda装置
Wang et al. Interference fading suppression with fault-tolerant Kalman filter in phase-sensitive OTDR
CN114485905A (zh) 一种基于不同脉宽光信号的分布式光纤振动与应变传感方法
Kito et al. Fast acquirable long-range measurement with frequency-swept probe BOTDA
Okamoto et al. Ultrafast measurement of optical DPSK signals using 1-symbol delayed dual-channel linear optical sampling
Ishimaru et al. Fast BOTDA acquiring method based on broadband light as a probe signal

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant