CN110455324B - 一种基于cp-φotdr的高重复率准分布式传感系统及其实现方法 - Google Patents

一种基于cp-φotdr的高重复率准分布式传感系统及其实现方法 Download PDF

Info

Publication number
CN110455324B
CN110455324B CN201910768983.6A CN201910768983A CN110455324B CN 110455324 B CN110455324 B CN 110455324B CN 201910768983 A CN201910768983 A CN 201910768983A CN 110455324 B CN110455324 B CN 110455324B
Authority
CN
China
Prior art keywords
signal
optical
representing
chirp
phase
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
CN201910768983.6A
Other languages
English (en)
Other versions
CN110455324A (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.)
University of Electronic Science and Technology of China
Original Assignee
University of Electronic Science and Technology of China
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 University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CN201910768983.6A priority Critical patent/CN110455324B/zh
Publication of CN110455324A publication Critical patent/CN110455324A/zh
Application granted granted Critical
Publication of CN110455324B publication Critical patent/CN110455324B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01D—MEASURING 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/00—Mechanical 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/26—Mechanical 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/32—Mechanical 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/34—Mechanical 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/353—Mechanical 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 influencing the transmission properties of an optical fibre
    • G01D5/35306—Mechanical 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 influencing the transmission properties of an optical fibre using an interferometer arrangement
    • G01D5/35309—Mechanical 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 influencing the transmission properties of an optical fibre using an interferometer arrangement using multiple waves interferometer
    • G01D5/35316—Mechanical 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 influencing the transmission properties of an optical fibre using an interferometer arrangement using multiple waves interferometer using a Bragg gratings

Landscapes

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

Abstract

本发明公开了一种基于CP‑ΦOTDR的高重复率准分布式传感系统及实现方法，属于光纤传感测量技术领域，包括激光器模块、耦合器、光信号调制模块、信号检测模块、光环形器、待测传感光纤，所述光信号调制模块将光信号调制为N个位于相同频带和不同频带的正负啁啾脉冲光信号，所述信号检测模块利用反射光信号和耦合器输入的本振光信号产生拍频信号，所述信号检测模块将所述拍频信号输入信号解调模块进行解调和输出，得到所述反射光信号的扰动信息；本发明在保证高空间分辨率的同时，使得测量重复率提升了N倍，不需要双通道探测，只需要单通道即可，降低了系统复杂度。

Description

一种基于CP-ΦOTDR的高重复率准分布式传感系统及其实现 方法
技术领域
本发明涉及光纤传感测量技术领域，具体涉及一种基于CP-ΦOTDR的高重复率准分布式传感系统及其实现方法。
背景技术
传感器作为一种集成自动测量、记录等功能的检测装置，在我们的生活及生产活动中得到了广泛的应用。相比于传统的电学传感器，基于光纤的传感器具有轻便、抗腐蚀、抗电磁干扰、耐高温、探测灵敏度高等独特的优势，所以分布式光纤传感在近几十年来得到快速发展。
相位敏感形光时域反射仪是一种典型的分布式光纤传感系统，它通过解调光纤中瑞利散射光的干涉强度来进行对外界干扰的实时监测，由于其具有超高灵敏度的显著优势而备受关注，但是由于瑞利散射光的强度很微弱而且反射率不稳定，所以在信噪比和传感距离等方面有很多的限制。
由于FBG的反射强度很高而且很稳定，所以基于弱反射率FBG串的准分布式测量系统受到广泛应用。但是不管是分布式测量系统还是准分布式测量系统，测量重复率fscan和传感距离L是相互制约的,他们之间满足fscan≤c/2nL。
近年来，对于打破测量重复率和传感距离限制的研究，主要有以下两种：
1、基于频分复用的分布式光纤传感技术，通过向传感光纤中注入N个不同频率范围的啁啾脉冲，通过在频域加窗的方法分别解调出来各个频率范围的散射信号，这些散射信号可以分别恢复出来扰动信息，从而测量重复率提升了N倍，但是这种方法使得传感带宽增加了N倍，而且由于啁啾信号扫频范围的减小，使得系统的空间分辨率也降低了N倍。
2、基于正负频IQ解调的分布式光纤传感技术，通过向传感光纤中注入正负频率的啁啾脉冲，由于正频和负频在频域上是可以分开的，所以运用傅里叶变换技术，可以将携带反射光信息的正负频率的信号分开，它们可以分别恢复出来扰动信息，使得测量重复率提升了2倍。但是这种方法需要双通道探测，增加了系统的复杂度。
发明内容
本发明的目的在于：本发明提供了一种基于CP-ΦOTDR的高重复率准分布式传感系统及其实现方法，解决了目前的准分布式测量系统由于测量重复率与传感距离相互制约，无法保证在系统具有高空间分辨率的同时提高重复率的技术问题。
本发明采用的技术方案如下：
一种基于CP-ΦOTDR的高重复率准分布式传感系统，包括激光器模块，所述激光器模块产生光信号至耦合器，所述耦合器分别将光信号输入光信号调制模块和信号检测模块，
所述光信号调制模块将光信号调制为N个位于相同频带和不同频带的正负啁啾脉冲光信号，所述光信号调制模块的输出端连接光环形器的端口1，所述光环形器的端口2将啁啾脉冲光信号输入待测传感光纤并接收所述待测传感光纤返回的反射光信号，所述光环形器的端口3将所述反射光信号输入信号检测模块，
所述信号检测模块利用反射光信号和耦合器输入的本振光信号产生拍频信号，所述信号检测模块将所述拍频信号输入信号解调模块进行解调和输出，得到所述反射光信号的扰动信息。
进一步的，所述光信号调制模块包括波形发生器和电光调制器。
进一步的，所述待测传感光纤中的反射点为光纤布拉格光栅(FBG)。
进一步的，所述信号检测模块包括光混合器和信号采集器。
一种基于CP-ΦOTDR的高重复率准分布式传感系统的实现方法，包括以下步骤：
步骤1：激光器模块产生光信号，并将所述光信号发送至耦合器，所述耦合器将光信号分别输入光信号调制模块和信号检测模块；
步骤2：所述光信号调制模块将光信号调制为N个位于相同频带和不同频带的正负啁啾脉冲光信号，并将所述啁啾脉冲光信号发送至光环形器；
步骤3：所述光环形器将所述啁啾脉冲光信号发送至待测传感光纤，所述待测传感光纤向所述光环形器返回反射光信号；
步骤4：所述光环形器将所述反射光信号发送至信号检测模块；
步骤5：所述信号检测模块利用反射光信号和耦合器输入的本振光信号产生拍频信号，并将所述拍频信号输入信号解调模块；
步骤6：所述信号解调模块将所述拍频信号进行解调和输出，得到反射光信号中的扰动信息。
进一步的，所述步骤1中，激光器模块产生的光信号EL(t)为：
Figure BDA0002172939310000021
其中，
Figure BDA0002172939310000022
表示光信号的初始相位，t表示时间，θ(t)表示激光器相位噪声，EL表示光信号的幅度。
进一步的，所述步骤2中，啁啾脉冲光信号Ep(t)为：
Figure BDA0002172939310000031
td表示正负啁啾脉冲之间的时延，γ表示啁啾脉冲光信号的啁啾率；
正啁啾信号sdigital(t)为：
Figure BDA0002172939310000032
负啁啾信号snegdigital(t)为：
Figure BDA0002172939310000033
其中，
Figure BDA0002172939310000034
表示反射光的相位项，
Figure BDA0002172939310000035
表示啁啾脉冲初相，τp表示脉冲宽度，E0表示啁啾脉冲信号的幅度，fm表示啁啾脉冲信号的中心频率，i表示反射点的序号，ti表示第i个反射点处的时延。
进一步的，所述步骤3中，反射光信号中第i个反射点处反射光信号Eref(ti,t)为：
Figure BDA0002172939310000036
其中，ti表示第i个反射点处的时延，R(ti)表示第i个反射点的反射率，
Figure BDA0002172939310000037
表示第i个反射点处反射信号的相位。
进一步的，所述步骤5中，第i个反射点处反射光信号与本振光信号的拍频信号EB(ti,t)为
Figure BDA0002172939310000038
Figure BDA0002172939310000039
其中，A(ti)表示第i个反射点反射信号的幅度，
Figure BDA00021729393100000310
表示本振与信号之间的相位差。
进一步的，所述步骤6中，信号解调模块将所述拍频信号EB(ti,t)做希尔伯特变换之后得到Ebeat(ti,t)，再进行解调和输出，包括利用正啁啾信号sdigital(t)进行解调和利用负啁啾信号snegdigital(t)进行解调，具体为：
拍频信号经过希尔伯特变换后得到：
Figure BDA0002172939310000041
Figure BDA0002172939310000042
利用正啁啾信号sdigital(t)进行解调，即利用正啁啾信号sdigital(t)对Ebeat(ti,t)做匹配滤波，当两者完全重合时，获得互相关峰：
Figure BDA0002172939310000043
Figure BDA0002172939310000044
Eneg表示在只用正啁啾信号做匹配滤波时负啁啾信号分量的影响；S(ti)表示拍频信号和正啁啾信号的互相关峰值；
则
Figure BDA0002172939310000045
将公式12向基带方向移频fm之后，取相位可得：
Figure BDA0002172939310000046
Figure BDA0002172939310000047
其中，
Figure BDA0002172939310000048
表示反射光信号匹配滤波峰值处的相位，
Figure BDA0002172939310000049
表示
Figure BDA00021729393100000410
求和之后的相位，
Figure BDA00021729393100000411
表示第ti+1和第ti个反射点峰值处相位差；
利用负啁啾信号snegdigital(t)进行解调，即利用负啁啾信号snegdigital(t)对Ebeat(ti,t)做匹配滤波，当两者完全重合时，获得互相关峰：
Figure BDA0002172939310000051
Figure BDA0002172939310000052
Epos表示在只用负啁啾信号做匹配滤波时正啁啾信号分量的影响；Sneg(ti+td)表示拍频信号和负啁啾信号的互相关峰值；
Figure BDA0002172939310000053
将公式17向基带方向移频fm之后，取相位可得：
Figure BDA0002172939310000054
Figure BDA0002172939310000055
其中，
Figure BDA0002172939310000056
表示反射信号匹配滤波峰值处的相位，
Figure BDA0002172939310000057
表示
Figure BDA0002172939310000058
求和之后的相位，
Figure BDA0002172939310000059
表示第ti+1和第ti个反射点峰值处相位之差。
综上所述，由于采用了上述技术方案，本发明的有益效果是：
在保证系统有高空间分辨率的同时，使得测量重复率提高N倍，而且本发明使用单通道探测，降低了系统的复杂度。例如，在使用相隔5米的FBG串作为传感光纤时，设置啁啾信号的扫频范围为20MHz，即为5m的空间分辨率，向光纤中注入由正啁啾信号和负啁啾信号组成的两个脉冲的周期信号，正负啁啾脉冲间隔为半个周期，使得系统在保证5米空间分辨率的同时，将测量重复率提升了2倍。
本发明相较于基于频分复用技术的分布式光纤传感系统，由于相同频带宽度的正负啁啾脉冲不需要额外增加传感带宽，位于不同频带宽度的正负啁啾脉冲只需要增加少量传感带宽，而且在准分布式系统中，空间分辨率是由啁啾的扫频范围和反射点之间的间隔共同决定的，当反射点间隔大于啁啾的扫频范围对应的空间分辨率时，空间分辨率与啁啾信号的扫频范围无关，所以本发明在保证高空间分辨率的同时，使得测量重复率提升了N倍。
本发明相较于基于正负频IQ解调技术的分布式光纤传感系统，本发明不需要双通道探测，只需要单通道即可，降低了系统复杂度。
附图说明
为了更清楚地说明本发明实施例的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，应当理解，以下附图仅示出了本发明的某些实施例，因此不应被看作是对范围的限定，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他相关的附图。
图1是本发明实施例1提供的一种基于CP-ΦOTDR的高重复率准分布式传感系统的系统框图；
图2是本发明实施例1中的正负啁啾脉冲调制信号的时域和频域图；
图3是本发明实施例1提供的一种基于CP-ΦOTDR的高重复率准分布式传感实现方法流程图；
图4是本发明实施例1中的正负啁啾信号解调过程原理图；
图5是本发明实施例2提供的另一种基于CP-ΦOTDR的高重复率准分布式传感系统的系统框图；
图6是本发明实施例3提供的又一种基于CP-ΦOTDR的高重复率准分布式传感系统的系统框图；
附图标记：
1-激光器模块，2-耦合器，3-光信号调制模块，31-波形发生器，32-电光调制器，4-光环形器，5-待测传感光纤，6-信号检测模块，61-光混合器，62-信号采集器，7-信号解调模块，8-偏振控制模块。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明，即所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。通常在此处附图中描述和示出的本发明实施例的组件可以以各种不同的配置来布置和设计。
因此，以下对在附图中提供的本发明的实施例的详细描述并非旨在限制要求保护的本发明的范围，而是仅仅表示本发明的选定实施例。基于本发明的实施例，本领域技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例，都属于本发明保护的范围。
需要说明的是，术语“第一”和“第二”等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
以下结合实施例对本发明的特征和性能作进一步的详细描述。
实施例1
本发明实施例中的高重复率准分布式传感系统主要基于光时域反射仪(opticaltime-domain reflectometer，OTDR)实现，具体地，是基于啁啾脉冲相位敏感光时域反射仪(CP-ΦOTDR)实现的。
如图1所示，本实施例提供的一种基于CP-ΦOTDR的高重复率准分布式传感系统包括激光器模块1，所述激光器模块1产生光信号至耦合器2，耦合器2为光耦合器，所述耦合器2分别将光信号输入光信号调制模块3和信号检测模块6，所述光信号调制模块3将光信号调制为N个位于相同频带和不同频带的正负啁啾脉冲光信号，所述光信号调制模块3的输出端连接光环形器4的端口1，所述光环形器4的端口2将啁啾脉冲光信号输入待测传感光纤5并接收所述待测传感光纤5返回的反射光信号，所述光环形器4的端口3将所述反射光信号输入信号检测模块6，
所述信号检测模块6利用反射光信号和耦合器2输入的本振光信号产生拍频信号，所述信号检测模块6将所述拍频信号输入信号解调模块7进行解调和输出，得到所述反射光信号的扰动信息。
所述光信号调制模块3包括波形发生器31和电光调制器32。所述待测传感光纤5为沿轴向具有一连串反射点的光纤，所述反射点可以为光纤布拉格光栅(FBG)、FP腔等，本实施例的反射点为光纤布拉格光栅(FBG)。所述信号检测模块6包括光混合器61和信号采集器62。
其中，光环形器4(optical circulator)是一种多端口非互易光学器件，其具有光导向作用，其典型结构有N(N大于等于3)个端口，当光由其中任一个端口输入(一般是端口1)时，可以几乎无损失地按照数字顺序由下一个端口(端口2)输出，而其它端口(端口3)处几乎没有光输出；以此类推，当光由端口2输入时，也可以由端口3近乎无损失的输出，于此同时，端口1或其他端口上没有光输出。光环形器4类型可以是透射式或反射式的光环形器4。
本实施例还提供一种基于CP-ΦOTDR的高重复率准分布式传感系统的实现方法，包括以下步骤(如图2所示)：
步骤1(S101)：激光器模块1产生光信号，并将所述光信号发送至耦合器2，所述耦合器2将光信号分别输入光信号调制模块3和信号检测模块6；
激光器模块1产生的光信号EL(t)为：
Figure BDA0002172939310000081
其中，
Figure BDA0002172939310000082
表示光信号的初始相位，t表示时间，θ(t)表示激光器相位噪声，EL表示光信号的幅度。
步骤2(S102)：所述光信号调制模块3将光信号调制为N个位于相同频带和不同频带的正负啁啾脉冲光信号(如图3所示为正负啁啾脉冲调制信号的时域和频域图像)，并将所述啁啾脉冲光信号发送至光环形器4；
啁啾脉冲光信号Ep(t)为：
Figure BDA0002172939310000083
td表示正负啁啾脉冲之间的时延，γ表示啁啾脉冲光信号的啁啾率；
正啁啾信号sdigital(t)为：
Figure BDA0002172939310000084
负啁啾信号snegdigital(t)为：
Figure BDA0002172939310000085
其中，t表示时间，θ(t)表示激光器相位噪声，
Figure BDA0002172939310000086
表示反射光的相位项，
Figure BDA0002172939310000087
表示啁啾脉冲初相，τp表示脉冲宽度，E0表示啁啾脉冲信号的幅度，fm表示啁啾脉冲信号的中心频率，i表示反射点的序号，ti表示第i个反射点处的时延。
步骤3(S103)：所述光环形器4将所述啁啾脉冲光信号发送至待测传感光纤5，所述待测传感光纤5向所述光环形器4返回反射光信号；
第i个反射点处反射光信号Eref(ti,t)为：
Figure BDA0002172939310000091
其中，i表示反射点的序号，ti表示第i个反射点处的时延，R(ti)表示第i个反射点的反射率，
Figure BDA0002172939310000092
表示第i个反射点处反射信号的相位，τp表示脉冲宽度，
Figure BDA0002172939310000093
为啁啾脉冲初相，t表示时间。
步骤4(S104)：所述光环形器4将所述反射光信号发送至信号检测模块6；
步骤5(S105)：所述信号检测模块6利用反射光信号和耦合器2输入的本振光信号产生拍频信号，并将所述拍频信号输入信号解调模块7，其中本振光信号即激光器产生的光信号；
第i个反射点处反射光信号与本振光信号的拍频信号EB(ti,t)为
Figure BDA0002172939310000094
Figure BDA0002172939310000095
其中，A(ti)表示第i个反射点反射信号的幅度，
Figure BDA0002172939310000096
表示本振与信号之间的相位差。
步骤6(S106)：信号解调模块7将所述拍频信号EB(ti,t)做希尔伯特变换之后得到Ebeat(ti,t)，再进行解调和输出，得到反射光信号中的扰动信息。解调原理如图4所示，拍频信号在做互相关时，只有响应函数与原始信号相同时才会匹配出来峰，如果与原始信号不相同，信号就会被抑制，由此可以通过匹配滤波算法分离N个位于相同频带和不同频带的正负啁啾信号，从而可以分别通过所述正负啁啾信号解调出来扰动信息。
包括利用正啁啾信号sdigital(t)进行解调和利用负啁啾信号snegdigital(t)进行解调，具体为：
拍频信号经过希尔伯特变换后得到：
Figure BDA0002172939310000101
Figure BDA0002172939310000102
利用正啁啾信号sdigital(t)进行解调，即利用正啁啾信号sdigital(t)对Ebeat(ti,t)做匹配滤波，当两者完全重合时，获得互相关峰：
Figure BDA0002172939310000103
Figure BDA0002172939310000104
Eneg表示在只用正啁啾信号做匹配滤波时负啁啾信号分量的影响；S(ti)表示拍频信号和正啁啾信号的互相关峰值；
其中，激光器模块1的相位噪声项θ目前不考虑，之后对时间积分可得：
则
Figure BDA0002172939310000105
由于
Figure BDA0002172939310000106
这一项积分之后很小，即在只用正啁啾做匹配滤波时，负啁啾分量的影响很小，所以可以忽略不计；
将公式31向基带方向移频fm之后，取相位可得：
Figure BDA0002172939310000107
Figure BDA0002172939310000108
其中，
Figure BDA0002172939310000109
表示反射光匹配滤波峰值处的相位，
Figure BDA00021729393100001010
表示
Figure BDA00021729393100001011
求和之后的相位，
Figure BDA00021729393100001012
表示第ti+1和第ti个反射元峰值处相位之差；
利用负啁啾信号snegdigital(t)进行解调，即利用负啁啾信号snegdigital(t)对Ebeat(ti,t)做匹配滤波，当两者完全重合时，获得互相关峰：
Figure BDA0002172939310000111
Figure BDA0002172939310000112
Epos表示在只用负啁啾做匹配滤波时，正啁啾分量的影响；Sneg(ti+td)表示拍频信号和负啁啾信号的互相关峰值；
则
Figure BDA0002172939310000113
Figure BDA0002172939310000114
由于
Figure BDA0002172939310000115
这一项积分之后很小，即在只用负啁啾做匹配滤波时，正啁啾分量的影响很小，所以可以忽略不计；
将上式向基带方向移频fm之后，取相位可得：
Figure BDA0002172939310000116
Figure BDA0002172939310000117
其中，
Figure BDA0002172939310000118
表示反射光匹配滤波峰值处的相位，
Figure BDA0002172939310000119
表示
Figure BDA00021729393100001110
求和之后的相位，
Figure BDA00021729393100001111
表示第ti+1和第ti个反射元峰值处相位之差。
整体的反射信号和拍频信号，只需对所有的i按照时间顺序排列即可。
显然，两个峰值处相位之差与△t线性相关，所以根据拍频信号的相位信息，应用正负啁啾信号可以分别解调出来扰动信息，实现高重复率准分布式传感。
由此可见，本发明实施例中的基于CP-ΦOTDR的高重复率准分布式传感实现方法，与现有技术相比，有以下优越性：相较于基于频分复用技术的分布式光纤传感系统，由于相同频带宽度的正负啁啾脉冲不需要额外增加传感带宽，位于不同频带宽度的正负啁啾脉冲只需要增加少量传感带宽，而且在准分布式系统中，空间分辨率是由啁啾的扫频范围和反射点之间的间隔共同决定的，当反射点间隔大于啁啾的扫频范围对应的空间分辨率时，空间分辨率与啁啾信号的扫频范围无关，所以本发明在保证高空间分辨率的同时，使得测量重复率提升了N倍；相较于基于正负频IQ解调技术的分布式光纤传感系统，本发明不需要双通道探测，只需要单通道即可，降低了系统复杂度。
实施例2
如图5所示，本发明的实施例2提供的一种基于CP-ΦOTDR的高重复率准分布式传感系统。在本发明实施例2中，信号调制模块中包括波形发生器31和电光调制器32，其中，光耦合器2的分路输出1与电光调制器32连接，分路输出2与信号检测模块6连接，波形发生器31与电光调制器32连接，电光调制器32的输出端与光环形器4的1端口连接，光环形器4的2端口与待测传感光纤5连接，光环形器4的3端口与信号检测模块6,连接。
在本发明实施例2中，使用波形发生器31对预先写入的正负啁啾波形进行处理，将正负啁啾脉冲的数字信号转换成电信号，然后通过电光调制器32产生正负啁啾脉冲光信号，产生的信号被光环形器4导入待测传感光纤5中。信号检测模块6通过分别接收光耦合器2输出两个分路的本振信号和反射信号强度信息，再利用信号解调模块7进行解调，基于拍频数据得到扰动信息的计算结果。
本发明实施例的方案可以基于波形发生器31和电光调制器32产生预定频率的正负啁啾脉冲，利用激光器提供的本振信号实现相干检测，提高信号测试的信噪比。
实施例3
如图6所示为本发明实施例3提供的一种基于CP-ΦOTDR的高重复率准分布式传感系统。
该实施例与实施例2之间的区别在于：信号检测模块6中包括光混合器61和信号采集器62，其中光混合器61采用光耦合器，具体的检测信号是本振光和信号光同时进入光混合器61中，继而使信号采集器62采集输出的拍频信号，接着经过信号解调模块7进行希尔伯特变换得到光信号的相位信息，进而解调得到扰动信息；另外在光混合器61的分路2另设置了一个偏振控制模块8，控制本振光的偏振态，进一步提高测量信号的信噪比。
由此可见，本实施例3中通过光混合器61和信号采集器62能够实现拍频信号的强度探测，通过信号解调模块7进行希尔伯特变换得到相位信息，进而得到扰动信息，此外，该实施例通过在本振端引入了偏振控制模块8得到较高信噪比。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。

Claims (8)

1.一种基于CP-ΦOTDR的高重复率准分布式传感系统，包括激光器模块(1)，所述激光器模块(1)产生光信号至耦合器(2)，所述耦合器(2)分别将光信号输入光信号调制模块(3)和信号检测模块(6)，其特征在于，
所述光信号调制模块(3)将光信号调制为N个位于相同频带和不同频带的正负啁啾脉冲光信号，所述光信号调制模块(3)的输出端连接光环形器(4)的端口1，所述光环形器(4)的端口2将啁啾脉冲光信号输入待测传感光纤(5)并接收所述待测传感光纤(5)返回的反射光信号，所述光环形器(4)的端口3将所述反射光信号输入信号检测模块(6)，
所述信号检测模块(6)利用反射光信号和耦合器(2)输入的本振光信号产生拍频信号，所述信号检测模块(6)将所述拍频信号输入信号解调模块(7)进行解调和输出，得到所述反射光信号的扰动信息；
激光器模块(1)产生的光信号EL(t)为：
Figure FDA0003007503480000011
其中，
Figure FDA0003007503480000012
表示光信号的初始相位，t表示时间，θ(t)表示激光器相位噪声，EL表示光信号的幅度；
啁啾脉冲光信号Ep(t)为：
Figure FDA0003007503480000013
td表示正负啁啾脉冲之间的时延，γ表示啁啾脉冲光信号的啁啾率；
正啁啾信号sdigital(t)为：
Figure FDA0003007503480000014
负啁啾信号snegdigital(t)为：
Figure FDA0003007503480000015
其中，
Figure FDA0003007503480000016
表示反射光的相位项，
Figure FDA0003007503480000017
表示啁啾脉冲初相，τp表示脉冲宽度，E0表示啁啾脉冲信号的幅度，fm表示啁啾脉冲信号的中心频率，i表示反射点的序号，ti表示第i个反射点处的时延。
2.根据权利要求1所述的一种基于CP-ΦOTDR的高重复率准分布式传感系统，其特征在于：所述光信号调制模块(3)包括波形发生器(31)和电光调制器(32)。
3.根据权利要求1所述的一种基于CP-ΦOTDR的高重复率准分布式传感系统，其特征在于：所述待测传感光纤(5)中的反射点为光纤布拉格光栅(FBG)。
4.根据权利要求1所述的一种基于CP-ΦOTDR的高重复率准分布式传感系统，其特征在于：所述信号检测模块(6)包括光混合器(61)和信号采集器(62)。
5.一种基于CP-ΦOTDR的高重复率准分布式传感系统的实现方法，其特征在于：包括以下步骤：
步骤1：激光器模块(1)产生光信号，并将所述光信号发送至耦合器(2)，所述耦合器(2)将光信号分别输入光信号调制模块(3)和信号检测模块(6)；
步骤2：所述光信号调制模块(3)将光信号调制为N个位于相同频带和不同频带的正负啁啾脉冲光信号，并将所述啁啾脉冲光信号发送至光环形器(4)；
步骤3：所述光环形器(4)将所述啁啾脉冲光信号发送至待测传感光纤(5)，所述待测传感光纤(5)向所述光环形器(4)返回反射光信号；
步骤4：所述光环形器(4)将所述反射光信号发送至信号检测模块(6)；
步骤5：所述信号检测模块(6)利用反射光信号和耦合器(2)输入的本振光信号产生拍频信号，并将所述拍频信号输入信号解调模块(7)；
步骤6：所述信号解调模块(7)将所述拍频信号进行解调和输出，得到反射光信号中的扰动信息；
所述步骤1中，激光器模块(1)产生的光信号EL(t)为：
Figure FDA0003007503480000021
其中，
Figure FDA0003007503480000022
表示光信号的初始相位，t表示时间，θ(t)表示激光器相位噪声，EL表示光信号的幅度；
所述步骤2中，啁啾脉冲光信号Ep(t)为：
Figure FDA0003007503480000023
td表示正负啁啾脉冲之间的时延，γ表示啁啾脉冲光信号的啁啾率；
正啁啾信号sdigital(t)为：
Figure FDA0003007503480000024
负啁啾信号snegdigital(t)为：
Figure FDA0003007503480000031
其中，
Figure FDA0003007503480000032
表示反射光的相位项，
Figure FDA0003007503480000033
表示啁啾脉冲初相，τp表示脉冲宽度，E0表示啁啾脉冲信号的幅度，fm表示啁啾脉冲信号的中心频率，i表示反射点的序号，ti表示第i个反射点处的时延。
6.根据权利要求5所述的一种基于CP-ΦOTDR的高重复率准分布式传感系统的实现方法，其特征在于：所述步骤3中，反射光信号中第i个反射点处反射光信号Eref(ti,t)为：
Figure FDA0003007503480000034
其中，ti表示第i个反射点处的时延，R(ti)表示第i个反射点的反射率，
Figure FDA0003007503480000035
表示第i个反射点处反射信号的相位。
7.根据权利要求6所述的一种基于CP-ΦOTDR的高重复率准分布式传感系统的实现方法，其特征在于：所述步骤5中，第i个反射点处反射光信号与本振光信号的拍频信号EB(ti,t)为
Figure FDA0003007503480000036
Figure FDA0003007503480000037
其中，A(ti)表示第i个反射点反射信号的幅度，
Figure FDA0003007503480000038
表示本振与信号之间的相位差。
8.根据权利要求7所述的一种基于CP-ΦOTDR的高重复率准分布式传感系统的实现方法，其特征在于：所述步骤6中，信号解调模块(7)将所述拍频信号EB(ti,t)做希尔伯特变换之后得到Ebeat(ti,t)，再进行解调和输出，包括利用正啁啾信号sdigital(t)进行解调和利用负啁啾信号snegdigital(t)进行解调，具体为：
拍频信号经过希尔伯特变换后得到：
Figure FDA0003007503480000041
Figure FDA0003007503480000042
利用正啁啾信号sdigital(t)进行解调，即利用正啁啾信号sdigital(t)对Ebeat(ti,t)做匹配滤波，当两者完全重合时，获得互相关峰：
Figure FDA0003007503480000043
Figure FDA00030075034800000412
Eneg表示在只用正啁啾信号做匹配滤波时负啁啾信号分量的影响；S(ti)表示拍频信号和正啁啾信号的互相关峰值；
Figure FDA0003007503480000044
将公式12向基带方向移频fm之后，取相位可得：
Figure FDA0003007503480000045
Figure FDA0003007503480000046
其中，
Figure FDA0003007503480000047
表示反射光信号匹配滤波峰值处的相位，
Figure FDA0003007503480000048
表示
Figure FDA0003007503480000049
求和之后的相位，
Figure FDA00030075034800000410
表示第ti+1和第ti个反射点峰值处相位差；
利用负啁啾信号snegdigital(t)进行解调，即利用负啁啾信号snegdigital(t)对Ebeat(ti,t)做匹配滤波，当两者完全重合时，获得互相关峰：
Figure FDA00030075034800000411
Figure FDA0003007503480000051
Epos表示在只用负啁啾信号做匹配滤波时正啁啾信号分量的影响；Sneg(ti+td)表示拍频信号和负啁啾信号的互相关峰值；
则
Figure FDA0003007503480000052
Figure FDA0003007503480000053
将公式17向基带方向移频fm之后，取相位可得：
Figure FDA0003007503480000054
Figure FDA0003007503480000055
其中，
Figure FDA0003007503480000056
表示反射信号匹配滤波峰值处的相位，
Figure FDA0003007503480000057
表示
Figure FDA0003007503480000058
求和之后的相位，
Figure FDA0003007503480000059
表示第ti+1和第ti个反射点峰值处相位之差。
CN201910768983.6A 2019-08-20 2019-08-20 一种基于cp-φotdr的高重复率准分布式传感系统及其实现方法 Active CN110455324B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910768983.6A CN110455324B (zh) 2019-08-20 2019-08-20 一种基于cp-φotdr的高重复率准分布式传感系统及其实现方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910768983.6A CN110455324B (zh) 2019-08-20 2019-08-20 一种基于cp-φotdr的高重复率准分布式传感系统及其实现方法

Publications (2)

Publication Number Publication Date
CN110455324A CN110455324A (zh) 2019-11-15
CN110455324B true CN110455324B (zh) 2021-07-16

Family

ID=68487943

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910768983.6A Active CN110455324B (zh) 2019-08-20 2019-08-20 一种基于cp-φotdr的高重复率准分布式传感系统及其实现方法

Country Status (1)

Country Link
CN (1) CN110455324B (zh)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111412935B (zh) * 2020-01-17 2021-08-10 电子科技大学 基于时分复用的高重复率准分布式传感系统
CN111609919B (zh) * 2020-06-09 2021-06-01 重庆大学 光纤分布式振动和损耗同时检测系统
CN111609875B (zh) * 2020-06-10 2021-12-28 电子科技大学 基于啁啾连续光的数字域可调分布式光纤传感系统及方法
CN111854815B (zh) * 2020-08-05 2022-04-01 国网四川省电力公司经济技术研究院 基于分数阶傅里叶变换的高速分布式光纤传感系统及方法
CN112484837B (zh) * 2020-11-24 2021-12-28 电子科技大学 光纤空间定位系统及其实现方法
CN113776567B (zh) * 2021-08-27 2024-06-14 武汉普赛斯电子技术有限公司 一种基于光频域反射计的偏振衰落消除装置及方法
CN115014407B (zh) * 2022-05-26 2024-03-29 中国人民解放军国防科技大学 一种光纤调制不稳定性演化特性分布式测量方法及装置
CN116907627B (zh) * 2023-09-13 2023-12-19 之江实验室 基于光程差辅助的大动态范围分布式相位传感方法和装置

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103115633B (zh) * 2013-01-29 2016-02-24 复旦大学 利用相位生成载波降低干涉路径散（反）射光干扰的方法
CN105157812A (zh) * 2015-09-18 2015-12-16 南京派光信息技术有限公司 基于数字增强干涉的高灵敏度准分布式光纤光栅振动传感器
CN106248119B (zh) * 2016-07-14 2018-11-13 南京大学 一种分布式超高速扰动定量检测方法及装置
CN109459072B (zh) * 2018-12-06 2021-03-26 电子科技大学 一种基于本振光正交相位调制的相干光时域反射计
CN110132330B (zh) * 2019-05-22 2021-08-10 电子科技大学 一种基于cp-φotdr的双折射分布式测量系统及其方法
CN110132331B (zh) * 2019-05-22 2021-05-07 电子科技大学 一种基于子啁啾脉冲提取的cotdr测量方法及其系统

Also Published As

Publication number Publication date
CN110455324A (zh) 2019-11-15

Similar Documents

Publication Publication Date Title
CN110455324B (zh) 一种基于cp-φotdr的高重复率准分布式传感系统及其实现方法
US11193818B2 (en) Frequency modulation demodulator based on fiber grating sensor array
EP3483572B1 (en) Distributed fibre sensing system and vibration detection and positioning method therefor
US7336365B2 (en) Method and apparatus for suppression of crosstalk and noise in time-division multiplexed interferometric sensor systems
EP2435795B1 (en) Method and apparatus for optical sensing
CN111854815B (zh) 基于分数阶傅里叶变换的高速分布式光纤传感系统及方法
CN110132329B (zh) 应力、温度和振动复合检测光纤传感器及信号处理方法
US7515273B2 (en) Method for measuring the brillouin shift distribution along an optical fiber based on the optical demodulation of the signals, and relevant apparatus
CN110132331B (zh) 一种基于子啁啾脉冲提取的cotdr测量方法及其系统
CN110108346B (zh) 基于延迟调相啁啾脉冲对的光纤振动传感器
CN110470327A (zh) 一种光时域分析仪和分析方法
CN111609875B (zh) 基于啁啾连续光的数字域可调分布式光纤传感系统及方法
US11828631B2 (en) Method and apparatus for suppression of noise due to local oscillator instability in a coherent fiber optical sensor
CN114543973B (zh) 一种分布式超高频振动信号测量方法及光纤传感器
CN111412935B (zh) 基于时分复用的高重复率准分布式传感系统
CN116907627B (zh) 基于光程差辅助的大动态范围分布式相位传感方法和装置
CN115900787A (zh) 光谱域反射仪的实现方法及系统
CN115290176A (zh) 基于低相干光源的准分布式振动传感系统及方法
CN109084816B (zh) 一种光纤传感信号的抗偏振衰落报警算法
CN219935149U (zh) 基于锁相环的光纤传感解调系统
CN113639774B (zh) 一种基于双波长双脉冲光源的准分布式传感装置
US20230288272A1 (en) Multi-wavelength laser for synchronously monitoring temperature and pressure of ocean
CN115355931B (zh) 基于传感光纤的分布式扰动测量方法、扰动测量系统
US20240175727A1 (en) Fibre Optic Sensing
CN117824818A (zh) 一种基于双光频率梳的光纤分布式声波传感方法及系统

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