CN103822646A - 光纤光栅解调***中长距离光延迟效应的解调方法 - Google Patents

Abstract

本发明涉及一种光纤光栅解调***中长距离光延迟效应的解调方法，用于消除长距离光缆传输带来的波长解调误差。本发明方法包括扫频速率切换解调、检测通道速率选择、FBG波长顺序矫正、FBG延迟校准参数计算和校准解调五个步骤，消除光延迟影响，得到FBG真实波长值。本发明方法能够保证基于扫频光源的光纤光栅解调***数据的完整性、正确性，适用于长距离且分布特征不明的光纤光栅传感网络的远距离、大容量、高速的动态波长解调。

Description

光纤光栅解调***中长距离光延迟效应的解调方法

技术领域

本发明涉及光纤光栅解调技术领域，具体地指一种基于扫频光源的光纤光栅解调***中长距离光延迟效应的解调方法。

背景技术

随着光纤传感技术的发展，光纤光栅传感器以其抗电磁干扰、耐腐蚀、高绝缘性、便于分布式测量等优点，在机械装备、土木工程、石油化工、航空航天、医疗、船舶工业等领域取得了广泛应用，成为传感领域发展较快的技术之一。FBG（Fiber Bragg Grating，光纤布喇格光栅）可以封装成温度、应变、压力、振动等传感器，应用在桥梁、大坝、隧道、港口翻车机和大型吊车以及石油化工等安全监测与健康诊断领域。在实际应用中，例如港口岸桥起重机、翻车机等在线监测***中，***要求需要部署数十个甚至成百上千个FBG传感器构成传感器阵列完成远距离、分布式、实时、多通路的在线监测。同时，必须能够对复杂的机械装备结构实现多参数动态信息的状态监测，例如高转速设备振动信号测量等，要求解调***具有高解调速率和高分辨率。可见，光纤光栅解调***必须满足远距离测量、高分辨率、高解调速率、大容量等性能要求。

通常，光纤光栅解调***与FBG传感器之间会通过长距离的光缆连接，达到远程监控的目的。然而，基于扫频光源的光纤光栅解调***，对FBG中心波长反射光接收时刻因长距离光缆传输耗时与发射时刻出现延迟滞后现象，使得FBG波长位置发生变化，给FBG中心波长解调带来不可忽略的误差影响。这种由长距离光缆引起的光延迟效应是使用扫频光源在光纤光栅解调过程中所不可避免的。因此，亟待研究一种技术方案以消除光延迟效应给波长解调带来的误差，能够实现校准计算出真实中心波长值，此外该技术方案能够自适应测量环境的变化，如光缆长度变化、环境温度与压力的变化、通道配置改变、FBG传感器部署位置更改等，保证传感***数据的正确性、完整性。

发明内容

本发明目的在于克服上述现有技术的不足而提供一种光纤光栅解调***中长距离光延迟效应的解调方法，该方法针对使用扫频光源在光纤光栅解调过程中出现的光延迟效应，利用每个FBG传感器对应的延迟校准参数对传感器波长进行实时同步校准解调得到真实中心波长值，以消除光延迟效应的影响。

实现本发明目的采用的技术方案是一种光纤光栅解调***中长距离光延迟效应的调解方法，该方法包括：

（1）切换扫频速率获得各速率下参考通道和多路检测通道的峰值计数值数据，并根据标定算法完成各速率下波峰波长值的标定，以及“去噪平均”处理后得到算法计算所需的各波峰计数值与波长值相应解调数据；

（2）选出各路检测通道所对应的两个扫频速率RateⅠ、RateⅡ，在所述两个扫频速率下，检测通道不能发生“周期旋转”；

（3）在工作速率RateW下，对各检测通道是否发生“距离乱序”现象进行判断和顺序矫正，得到所述各检测通道中检测波峰排序与真实波长排序的对应关系；

（4）利用所述各检测通道对应扫频速率RateⅠ、RateⅡ下的解调数据完成全部FBG对应延迟校准参数ΔN^fbg的计算；

（5）解调***至工作速率RateW，利用以上步骤计算得到的各ΔN^fbg参数和对各检测通道FBG波长顺序矫正后得到的检测波峰排序与真实波长排序的对应关系对解调数据进行实时校准计算，得到真实中心波长值λ^fbg输出，消除光延迟效应对波长解调的影响。

本发明具有以下优点：

本发明适用于分布式FBG传感网络的远距离、大容量、高速的动态波长解调，各检测通道对应光缆距离长度可互不相同。解调***能够随工程实际测量场景和环境因素的变化实时更新FBG延迟校准参数，从而实现快速的自适应校准解调。延迟校准参数是测量环境中能够影响光信号传输延迟的众多因素的综合表述，包括了光缆距离长度，环境温度、压力等等。与传统通过光时域反射计（OTDR）测量距离来消除误差的方法相比，本发明操作简单，计算精度高，受环境影响小，且不会增加工程施工的难度和复杂度。

附图说明

图1为基于扫频光源的光纤光栅解调***原理框图。

图2为长距离光延迟效应示意图。

图3为“周期旋转”现象示意图。

图4为本发明方法流程图。

图5为本发明所涉及算法的流程图。

图6为单一FBG解调判断扫频速率是否发生“周期旋转”的示意图。

图7为多FBG解调判断扫频速率是否发生“周期旋转”的示意图。

图8为“乱序”现象示意图。

图9为FBG波长校准参数计算示意图。

具体实施方式

下面结合附图和具体实施例对本发明作进一步地详细说明。

本发明是针对在使用基于扫频光源的光纤光栅解调***的过程中，由于长距离光延迟效应产生调解误差而提出的一种解调方法。

本发明首先分析基于扫频光源的光纤光栅解调***产生解调误差的原因。

如图1所示，基于扫频光源的光纤光栅解调***包括：扫频光源模块、参考通道、多路检测通道、信号采集与处理模块和上位机软件***。

扫频光源模块中的扫频激光源在波长扫描驱动电压控制下，随电压的变化输出不同波长的光，扫频光进入参考通道（梳状滤波器）和多路检测通道。参考通道用来提高解调***的精度与稳定性，参考通道梳状滤波器透射光波峰与检测通道FBG反射光波峰由信号采集与处理模块进行检测和处理完成FBG传感器解调。扫频光源的波长扫描驱动电压信号频率决定了解调***的解调速率。

扫频激光源的出光波长与波长扫描驱动电压一一对应，同时，在光源使能控制信号作用下，只有前半周期有激光输出，而FBG反射波峰的检测则贯穿整个周期，如图2所示。参考通道透射波峰和检测通道FBG反射波峰信号被光电二极管检测采集后以峰值计数值N的形式表示其在一个驱动电压周期中的相对时间位置t，如图2中d所示。根据N对应的波长扫描驱动电压U标定出FBG反射中心波长值λ^fbg。

FBG中心波长对应光信号在t^Out时刻发射，光信号在光缆中传输耗时τ，FBG反射光在t^In（t^In=t^Out+τ）时刻被解调***接收。t^Out与t^In分别对应扫描电压U^Out、U^In，对应峰值计数值为N^Out、N^In，如图2中b、c、d所示。则解调***由N^In标定波长值λ^In与真实中心波长λ^fbg（即λ^Out）之间存在波长误差Δλ（Δλ=λ^In-λ^Out），波长解调发生错误。以4kHz解调频率为例，1km的距离会带来约4000的计数差值（ΔN=N^In-N^Out），造成约4nm的解调误差。由于在实际工程应用中，光缆距离可能会是几十甚至几百公里，因此，这种光延迟效应给基于扫频光源的解调***带来了不可忽略的解调误差。

随传输光缆距离的增加，延迟时间τ会进一步增长，过长的延时会导致本扫描周期发射的光信号被FBG反射后却延迟在下个周期或第n（n>1）周期被接收到而造成“周期旋转”延迟。如图3所示，扫频速率为4kHz，λ^out=λ^fbg=1290nm，光缆距离Len_optic cable=23Km，延迟导致ΔN约为92000、Δλ约为92nm，发射光在第二个周期被接收，而“周期旋转”延迟致使检测到Δλ=λ^In-λ^Out=-8nm。

本发明采用光源使能控制信号以控制扫频光源是否出光（H：有光输出；L：无光输出）；在一个工作周期中，波长扫描驱动电压上升阶段激光器输出扫频光，下降阶段无光输出；扩展峰值计数值N^min～N^max为0-100000，扩展波长解调范围λ^min～λ^max为100nm，以解决反射波峰出现在无光输出阶段而导致无法进行波长标定及校准计算的情况。图2中d图中虚线表示峰值计数值与波长解调范围的扩展部分。光延迟效应下的解调波长λ^In是带有误差的波长值，需要利用算法消除该误差得到真实中心波长值。如图4和图5所示，本发明基于扫频光源的光纤光栅解调***中长距离光延迟效应的调解方法包括以下步骤：

S100、解调***切换扫频速率进行数据记录与波长解调。解调***根据具体工程需要设置多个扫频速率档位，其中包含了高扫频速率档位（几千Hz）和低扫频速率档位（几百Hz或更低），工作速率RateW为最高速率。获得各速率下参考通道和多路检测通道解调的波峰峰值计数值数据。根据现有波长标定算法完成各速率下波峰波长值的标定，“去噪平均”等处理后得到算法计算所需的各波峰计数值与解调波长数据。多路检测通道可在光缆距离长度、FBG部署方案互不相同情况下进行同步解调。

多路检测通道可在光缆距离长度、FBG部署方案互不相同情况下进行同步解调。

S200、切换速率进行数据记录与波长解调完成后进行各检测通道速率选择。选出各检测通道对应扫频速率RateⅠ、RateⅡ，在所述两个扫频速率下，检测通道不能发生“周期旋转”，具体步骤为：

当检测通道只解调单一FBG传感器时，按速率从低到高的顺序，迭代对比各速率下峰值计数值，找到解调波长值（或峰值计数值）保持增长趋势的速率集合N（速率档位集合M表示，

），选择集合N中最大速率为RateⅡ，次高速率为RateⅠ，RateⅠ＜RateⅡ。

如图6所示，检测通道解调一个FBG传感器（λ^fbg=1310nm），光缆距离Len_optic cable=45Km，4kHz、2kHz、1kHz、500Hz和100Hz构成速率集合M（4kHz为工作速率）。从100Hz开始迭代对比解调波长值（或峰值计数值），100Hz、500Hz和1kHz下保持增长趋势（λ_0.1λ_0.5λ_1k，N_0.1N_0.5N_1k，），则100Hz、500Hz和1kHz构成集合N。可以判断集合N中速率未发生“周期旋转”，其余速率发生。选择1kHz作为RateⅡ，选择500Hz为RateⅠ。

当检测通道解调多FBG传感器时，首先判断检测通道在各速率下是否发生“周期旋转”。方法是在某速率下按检测通道波峰计数值顺序依次计算两个相邻FBG波峰峰值计数值之差ΔN_fbg的绝对值|ΔN_fbg|，取最大值|ΔN_fbg|_max与峰值计数值中值N_mid比较，若|ΔN_fbg|_max≥N_mig，则发生“周期旋转”；反之则无。在各速率中找出发生“周期旋转”的最小速率，选择比之更小的两个速率档位依次作为RateⅠ、RateⅡ，RateⅠ＜RateⅡ。

如图7所示，检测通道中接有五个FBG传感器（λ₁=1290nm、λ₂=1300nm、λ₃=1310nm、λ₄=1320nm、λ₅=1325nm），光缆距离Len_optic cable=17Km，4kHz、2kHz、1kHz、500Hz和100Hz构成速率集合M（500Hz和100Hz解调结果不会发生“周期旋转”，图中未列出，4kHz为工作速率）。按照上述检测通道是否发生“周期旋转”的方法，可以判断：4kHz时发生“周期旋转”，2kHz、1kHz未发生。选择2kHz作为RateⅡ，选择1kHz为RateⅠ。

各检测通道在选出的扫频速率RateⅠ、扫频速率RateⅡ后，从S100结果中获得RateⅠ、RateⅡ和工作速率RateW（本实施例设为4kHz）下参考通道和检测通道峰值计数值数据与解调波长值数据，分别表示如下：

分别为RateⅠ、RateⅡ扫频速率下参考通道梳状滤波器通道各波峰计数值（SZ代表梳状滤波器通道，m代表波峰编号）；

梳状滤波器通道各波峰中心波长值，常量；

RateⅠ扫频速率下某检测通道中FBG传感器解调波长值；

RateⅠ扫频速率下某检测通道中FBG传感器波峰计数值；

RateⅡ扫频速率下某检测通道中FBG传感器解调波长值；

RateⅡ扫频速率下某检测通道中FBG传感器波峰计数值；

工作速率下某检测通道中FBG传感器解调波长值；

工作速率下某检测通道中FBG传感器波峰计数值。

S300、检测通道FBG波长顺序矫正。

在“周期旋转”的影响下，检测通道在进行多FBG传感器同步解调时会出现两种FBG检测波峰排列“乱序”现象（检测通道在进行单一FBG传感器解调时无“乱序”现象）。第一种是在单一扫频速率下，检测波峰排序与波峰真实波长排序不一致带来的“距离乱序”；第二种是在相同距离情况下，不同扫频速率下检测波峰排序互不相同而造成“速率乱序”。

在五个FBG传感器同图7设置情况下，如图8所示，光缆距离为Len_optic cable=42Km，（4kHz为工作速率）检测波峰无延迟时排序应与真实波长顺序一致，对应关系为{λ₁、λ₂、λ₃、λ₄、λ₅}→{N₁、N₂、N₃、N₄、N₅}，而4kHz扫描速率下“距离乱序”导致检测波峰排序与真实波长排序对应关系为{λ₁、λ₂、λ₃、λ₄、λ₅}→{N₄、N₅、N₁、N₂、N₃}，2kHz扫描速率下“距离乱序”导致检测波峰排序与真实波长排序对应关系为{λ₁、λ₂、λ₃、λ₄、λ₅}→{N₂、N₃、N₄、N₅、N₁}；4kHz与2kHz之间由于“距离乱序”现象不同而存在“速率乱序”。

可见，在工作速率下进行波长校准解调之前必须对各检测通道中多FBG完成波长顺序矫正，找到检测波峰排序与其真实波长排序的对应关系。

本实施例中，由于检测通道在RateⅠ、RateⅡ速率下无“周期旋转”与“速率乱序”现象发生，所以只需对工作速率下各检测通道进行“周期旋转”现象的判断，若发生“周期旋转”则找出“距离乱序”下检测波峰排序与真实波长排序的对应关系，完成FBG波长顺序矫正。

首先，按下面关系式判断检测通道在工作速率下是否发生“周期旋转”：

若发生“周期旋转”，则工作速率存在“距离乱序”。记录下峰值差值|ΔN_fbg|=|ΔN_fbg|_max的两个波峰N_m与N_m+1（N_m+1>N_m）。将峰值计数值比N_m小的若干波峰分为A组，峰值计数值比N_m+1大的若干波峰分为B组。A组波峰为真实波长较大的FBG传感器对应延迟波峰，B组为真实波长值较小的传感器对应延迟波峰。将两者按{B、A}的顺序拼接即为真实波长排序，从而找出检测波峰排序与真实波长排序的对应关系。

如图7中4kHz（工作速率）解调结果图，N_m对应为N₅，N_m+1对应为N₁。A组即为{N₄、N₅}，B组即为{N₁、N₂、N₃}，拼接为{N₄、N₅、N₁、N₂、N₃}，得到检测波峰排序与真实波长排序对应关系为{λ₁、λ₂、λ₃、λ₄、λ₅}→{N₄、N₅、N₁、N₂、N₃}。

S400、计算FBG波长校准参数。

以单个FBG传感器为单位依次对各检测通道中所接FBG传感器进行逐一分析计算，包括以下步骤：

S401、计算理想波长值λ_Theo。

理论波长值是在扫描驱动电压与激光源出光波长严格线性条件下使用理论波长公式推导得到的。

理论波长通过下式计算：

$\left\{\begin{array}{c}({\mathrm{\λ}}_{\mathrm{RateII}}^{\mathrm{fbg}}-{\mathrm{\λ}}_{\mathrm{Theo}})\×{\mathrm{Coeff}}^{\mathrm{RateII}}=\frac{2\×{\mathrm{Len}}_{\mathrm{opticcable}}\×n}{c}\\ ({\mathrm{\λ}}_{\mathrm{RateI}}^{\mathrm{rbg}}-{\mathrm{\λ}}_{\mathrm{Theo}})\×{\mathrm{Coeff}}^{\mathrm{RateI}}=\frac{2\×{\mathrm{Len}}_{\mathrm{poticcable}}\×n}{c}\end{array}\right.$

式中，

为RateⅡ解调速率下FBG传感器解调波长值，

为RateⅠ解调速率下FBG传感器解调波长值，c代表光在真空中传播速度；n代表光缆介质的光折射率，Coeff^RateⅡ、Coeff^RateⅠ表示与解调***有关的***常量参数。

将λ_Theo与进行对比，检测出与λTheo相邻的两个梳状滤波器波峰，编号分别为Theo-peak-mR、Theo-peak-mL，并记录下对应波峰计数值数据。

其中波长较大的一个表示为：

${\mathrm{\&lambda;}}_{\mathrm{Theo}-\mathrm{peak}-mR}^{\mathrm{SZ}},({\mathrm{\&lambda;}}_{\mathrm{Theo}}<{\mathrm{\&lambda;}}_{\mathrm{Theo}-\mathrm{peak}-\mathrm{mR}}^{\mathrm{SZ}})$

称为λ_Theo“右侧梳状峰”；

波长较小的一个表示为：

${\mathrm{\&lambda;}}_{\mathrm{Theo}-\mathrm{peak}-\mathrm{mL}}^{\mathrm{SZ}},({\mathrm{\&lambda;}}_{\mathrm{Theo}}<{\mathrm{\&lambda;}}_{\mathrm{Theo}-\mathrm{peak}-\mathrm{mL}}^{\mathrm{SZ}})$

称为λ_Theo“左侧梳状峰”。

分别对应两个速率下的波峰计数值为：

$\left\{\begin{array}{cc}{N}_{\mathrm{RateI}-\mathrm{Theo}-\mathrm{peak}-\mathrm{mL}}^{\mathrm{SZ}},& N_{\mathrm{RateI}-\mathrm{Theo}-\mathrm{peak}-\mathrm{mR}}^{\mathrm{SZ}}\\ N_{\mathrm{RateII}-\mathrm{Theo}-\mathrm{peak}-\mathrm{mL},}^{\mathrm{SZ}}& N_{\mathrm{RateII}-\mathrm{Theo}-\mathrm{peak}-\mathrm{mR}}^{\mathrm{SZ}}\end{array}\right.$

S402、进行校准波长值λ_Cali计算。

在理论波长值λ_Theo的基础上，利用参考通道的多个透射波峰对扫频激光源输出光谱进行分段细分线性矫正后得到λ_Cali。

校准波长通过下式计算：

${\mathrm{\&lambda;}}_{\mathrm{Cali}}={\mathrm{\&lambda;}}_{\mathrm{Theo}-\mathrm{peak}-\mathrm{mR}}^{\mathrm{SZ}}-\frac{{\mathrm{\&beta;}}_{\mathrm{\&lambda;NRateI}}\&times;{\mathrm{\&beta;}}_{\mathrm{\&lambda;NRateII}}\&times;(N_{\mathrm{RateII\&lambda;R}}\&times;T_{\mathrm{cntRateII}}-N_{\mathrm{RateI\&lambda;R}}\&times;T_{\mathrm{cntRateI}})}{({\mathrm{\&beta;}}_{\mathrm{\&lambda;NRateII}}\&times;T_{\mathrm{cntRateI}}-{\mathrm{\&beta;}}_{\mathrm{\&lambda;NRateI}}\&times;T_{\mathrm{cntRateII}})}$

式中：

$\left\{\begin{array}{c}{\mathrm{\&beta;}}_{\mathrm{\&lambda;NRateI}}=\frac{({\mathrm{\&lambda;}}_{\mathrm{Theo}-\mathrm{peak}-\mathrm{mR}}^{\mathrm{SZ}}-{\mathrm{\&lambda;}}_{\mathrm{Theo}-\mathrm{peak}-\mathrm{mL}}^{\mathrm{SZ}})}{(N_{\mathrm{RateI}-\mathrm{Theo}-\mathrm{peak}-\mathrm{mR}}^{\mathrm{SZ}}-N_{\mathrm{RateI}-\mathrm{Theo}-\mathrm{peak}-\mathrm{mL}}^{\mathrm{SZ}})}\\ {\mathrm{\&beta;}}_{\mathrm{\&lambda;NRateII}}=\frac{({\mathrm{\&lambda;}}_{\mathrm{Theo}-\mathrm{peak}-\mathrm{mR}}^{\mathrm{SZ}}-{\mathrm{\&lambda;}}_{\mathrm{Theo}-\mathrm{peak}-\mathrm{mL}}^{\mathrm{SZ}})}{(N_{\mathrm{RateII}-\mathrm{theo}-\mathrm{peak}-\mathrm{mR}}^{\mathrm{SZ}}-N_{\mathrm{RateII}-\mathrm{Theo}-\mathrm{peak}-\mathrm{mL}}^{\mathrm{SZ}})}\end{array}\right.$

$\left\{\begin{array}{c}{N}_{\mathrm{RateI\&lambda;R}}=N_{\mathrm{RateI}}^{\mathrm{fbg}}-N_{\mathrm{RateI}-\mathrm{Theo}-\mathrm{peak}-\mathrm{mR}}^{\mathrm{SZ}}\\ N_{\mathrm{RateII\&lambda;R}}=N_{\mathrm{RateII}}^{\mathrm{fbg}}-N_{\mathrm{RateII}-\mathrm{Theo}-\mathrm{peak}-\mathrm{mR}}^{\mathrm{SZ}}\end{array}\right.$

T_cntRateⅠ、T_cntRateⅡ分别表示RateⅠ、RateⅡ速率下峰值计数值的计数时钟周期。

将λ_Cali与

进行对比，同理检测出与λCali相邻的两个梳状滤波器波峰，编号分别为Cali-peak-mR、Cali-peak-mL，并记录下对应波峰计数值数据。

其中波长较大的一个表示为：

${\mathrm{\&lambda;}}_{\mathrm{Cali}-\mathrm{peak}-mR}^{\mathrm{SZ}},({\mathrm{\&lambda;}}_{\mathrm{Cali}}<{\mathrm{\&lambda;}}_{\mathrm{Cali}-\mathrm{peak}-\mathrm{mR}}^{\mathrm{SZ}})$

称为λ_Cali“右侧梳状峰”；

波长较小的一个表示为：

${\mathrm{\&lambda;}}_{\mathrm{Cali}-\mathrm{peak}-\mathrm{mL}}^{\mathrm{SZ}},({\mathrm{\&lambda;}}_{\mathrm{Cali}}<{\mathrm{\&lambda;}}_{\mathrm{Cali}-\mathrm{peak}-\mathrm{mL}}^{\mathrm{SZ}})$

称为λ_Cali“左侧梳状峰”。

分别对应两个速率下的波峰计数值为：

$\left\{\begin{array}{cc}{N}_{\mathrm{RateI}-\mathrm{Cali}-\mathrm{peak}-\mathrm{mL}}^{\mathrm{SZ}},& N_{\mathrm{RateI}-\mathrm{Cali}-\mathrm{peak}-\mathrm{mR}}^{\mathrm{SZ}}\\ N_{\mathrm{RateII}-\mathrm{Cali}-\mathrm{peak}-\mathrm{mL},}^{\mathrm{SZ}}& N_{\mathrm{RateII}-\mathrm{Cali}-\mathrm{peak}-\mathrm{mR}}^{\mathrm{SZ}}\end{array}\right.$

S403、进行条件比较与再计算，获得延迟校准参数。

位置比较条件是：λ_Cali与λ_Theo的“左、右梳状峰”是否对应相同（即Theo-peak-mR是否等于Cali-peak-mR，Theo-peak-mL是否等于Cali-peak-mL），如图9所示，若满足位置比较条件，则直接计算延迟校准参数；若不满足位置比较条件，则需要再计算和比较，直到满足比较条件为止。

再计算则是把λ_Cali的“左、右梳状峰”替换λ_Theo的“左、右梳状峰”，利用校准波长公式再次计算λ_Cali。后续步骤如上。最终得到满足位置比较条件的λ_Cali。

得到最终λ_Cali后，利用延迟校准参数计算公式计算各FBG对应的波长校准的参数ΔN^fbg。

计算公式：

$\mathrm{\&Delta;}{N}^{\mathrm{fbg}}=\frac{(N_{\mathrm{RateII\&lambda;R}}+\frac{({\mathrm{\&lambda;}}_{\mathrm{Cali}-\mathrm{peak}-\mathrm{mR}}^{\mathrm{SZ}}-{\mathrm{\&lambda;}}_{\mathrm{Cali}})}{{\mathrm{\&beta;}}_{\mathrm{\&lambda;NRateII}}})}{T_{\mathrm{cntRateII}}}$

延迟校准参数ΔN^fbg是在具体工程应用场合中能够影响光信号传输延迟时间长度因素的综合体现，包括了光缆距离长度（Len_optic cable）、环境温度、压力等。

S500、校准解调。

解调***至工作速率，利用以上步骤计算得到的各ΔN^fbg参数和FBG波长顺序矫正后得到的检测波峰排序与真实波长排序的对应关系对解调数据进行实时校准计算，得到真实中心波长值λ^fbg输出，消除光延迟效应对波长解调的影响。

首先，对采集到的FBG峰值计数值

进行实时校准计算，得到

峰值计数值校准计算公式如下：

$N_{\mathrm{RateW}-\mathrm{Cali}}^{\mathrm{fbg}}=N_{\mathrm{Ratew}}^{\mathrm{fbg}}-\mathrm{\&Delta;}{N}^{\mathrm{fbg}}$

然后，根据

标定得到FBG中心波长λ^fbg，即消除光信号传输延迟影响后的FBG中心波长值。

Claims (8)

1.一种光纤光栅解调***中长距离光延迟效应的解调方法，其特征在于：

（1）切换扫频速率获得各速率下参考通道和多路检测通道的峰值计数值数据，并根据标定算法完成各速率下波峰波长值的标定，以及“去噪平均”处理后得到算法计算所需的各波峰计数值与波长值相应解调数据；

（2）选出各路检测通道所对应的两个扫频速率RateⅠ、RateⅡ，在所述两个扫频速率下，检测通道不能发生“周期旋转”；

（3）在工作速率RateW下，对各检测通道是否发生“距离乱序”现象进行判断和顺序矫正，得到所述各检测通道中检测波峰排序与真实波长排序的对应关系；

（4）利用所述各检测通道对应扫频速率RateⅠ、RateⅡ下的解调数据完成全部FBG对应延迟校准参数ΔN^fbg的计算；

（5）解调***至工作速率RateW，利用以上步骤计算得到的各ΔN^fbg参数和对各检测通道FBG波长顺序矫正后得到的检测波峰排序与真实波长排序的对应关系对解调数据进行实时校准计算，得到真实中心波长值λ^fbg输出，消除光延迟效应对波长解调的影响。

2.根据权利要求1所述光纤光栅解调***中长距离光延迟效应的解调方法，其特征在于，步骤（1）包括：

设置多个扫频速率档位，其中工作速率RateW作为最高速率；

由高到低，解调***切换所述扫频速率档位进行解调获得各速率下参考通道和各检测通道解调得到的波峰峰值计数值数据。

3.根据权利要求1所述光纤光栅解调***中长距离光延迟效应的解调方法，其特征在于，步骤（2）包括：

分析检测通道在各速率下的解调数据，根据检测通道解调FBG传感器的情况，各检测通道选出不会发生“周期旋转”所对应的两个扫频速率RateⅠ、RateⅡ，检测通道在对应RateⅠ、RateⅡ、工作速率RateW下FBG传感器解调波长值分别为RateⅠ、RateⅡ、工作速率RateW下传感器波峰计数值分别为

RateⅠ、RateⅡ解调速率下参考通道各波峰计数值分别为

参考通道各波峰中心波长值为

4.根据权利要求3所述光纤光栅解调***中长距离光延迟效应的解调方法，其特征在于：

若检测通道只解调单一FBG传感器，按速率从低到高的顺序，迭代对比各速率下峰值计数值，找到解调波长值保持增长趋势的速率集合N，

M表示速率档位集合，选择集合N中最大速率为RateⅡ，次高速率为RateⅠ，RateⅠ<RateⅡ。

5.根据权利要求3所述光纤光栅解调***中长距离光延迟效应的解调方法，其特征在于：

若检测通道解调多FBG传感器，首先按下式判断检测通道在各速率下是否发生“周期旋转”，

上式中，|ΔN_fbg|_max为按检测通道波峰计数值顺序依次计算两个相邻FBG波峰峰值计数值之差ΔN_fbg的绝对值|ΔN_fbg|的最大值；N_mid为峰值计数值中值。

在各速率中找出发生“周期旋转”的最小速率，选择比之更小的两个速率档位依次作为RateⅠ、RateⅡ，RateⅠ<RateⅡ。

6.根据权利要求1所述光纤光栅解调***中长距离光延迟效应的解调方法，其特征在于，所述步骤（3）包括：

若检测通道解调单一FBG传感器则不存在“乱序”现象；若检测通道解调多FBG传感器且发生“周期旋转”，则工作速率RateW存在“距离乱序”，记录下峰值差值|ΔN_fbg|=|ΔN_fbg|_max的两个波峰N_m与N_m+1，将峰值计数值比N_m小的全部波峰分为A组，峰值计数值比N_m+1大的全部波峰分为B组，A组波峰为真实波长较大的FBG传感器对应延迟波峰，B组为真实波长值较小的传感器对应延迟波峰，将两者按{B、A}的顺序拼接即为真实波长排序，从而找出检测波峰排序与真实波长排序的对应关系。

7.根据权利要求1所述光纤光栅解调***中长距离光延迟效应的解调方法，其特征在于，步骤（4）包括：

S401根据下式进行理想波长值λ_Theo的计算，

$\left\{\begin{array}{c}({\mathrm{\&lambda;}}_{\mathrm{RateII}}^{\mathrm{fbg}}-{\mathrm{\&lambda;}}_{\mathrm{Theo}})\&times;{\mathrm{Coeff}}^{\mathrm{RateII}}=\frac{2\&times;{\mathrm{Len}}_{\mathrm{opticcable}}\&times;n}{c}\\ ({\mathrm{\&lambda;}}_{\mathrm{RateI}}^{\mathrm{rbg}}-{\mathrm{\&lambda;}}_{\mathrm{Theo}})\&times;{\mathrm{Coeff}}^{\mathrm{RateI}}=\frac{2\&times;{\mathrm{Len}}_{\mathrm{poticcable}}\&times;n}{c}\end{array}\right.$

式中，为RateⅡ扫频速率下FBG传感器解调波长值，

为RateⅠ扫频速率下FBG传感器解调波长值，c代表光在真空中传播速度；n代表光缆介质的光折射率，Coeff^RateⅡ、Coeff^RateⅠ表示与解调***的***常量参数；

将λ_Theo与

进行对比，检测出与λ_Theo相邻的两个梳状滤波器波峰，编号分别为Theo-peak-mR、Theo-peak-mL，并记录下对应波峰计数值数据，其中波长较大的一个表示为：

${\mathrm{\&lambda;}}_{\mathrm{Theo}-\mathrm{peak}-mR}^{\mathrm{SZ}},({\mathrm{\&lambda;}}_{\mathrm{Theo}}<{\mathrm{\&lambda;}}_{\mathrm{Theo}-\mathrm{peak}-\mathrm{mR}}^{\mathrm{SZ}})$

称为λ_Theo“右侧梳状峰”；

波长较小的一个表示为：

${\mathrm{\&lambda;}}_{\mathrm{Theo}-\mathrm{peak}-\mathrm{mL}}^{\mathrm{SZ}},({\mathrm{\&lambda;}}_{\mathrm{Theo}}<{\mathrm{\&lambda;}}_{\mathrm{Theo}-\mathrm{peak}-\mathrm{mL}}^{\mathrm{SZ}})$

称为λ_Theo“左侧梳状峰”。

分别对应两个扫频速率下的波峰计数值为：

$\left\{\begin{array}{cc}{N}_{\mathrm{RateI}-\mathrm{Theo}-\mathrm{peak}-\mathrm{mL}}^{\mathrm{SZ}},& N_{\mathrm{RateI}-\mathrm{Theo}-\mathrm{peak}-\mathrm{mR}}^{\mathrm{SZ}}\\ N_{\mathrm{RateII}-\mathrm{Theo}-\mathrm{peak}-\mathrm{mL},}^{\mathrm{SZ}}& N_{\mathrm{RateII}-\mathrm{Theo}-\mathrm{peak}-\mathrm{mR}}^{\mathrm{SZ}}\end{array}\right.$

S402、通过下式计算校准波长值λ_Cali

${\mathrm{\&lambda;}}_{\mathrm{Cali}}={\mathrm{\&lambda;}}_{\mathrm{Theo}-\mathrm{peak}-\mathrm{mR}}^{\mathrm{SZ}}-\frac{{\mathrm{\&beta;}}_{\mathrm{\&lambda;NRateI}}\&times;{\mathrm{\&beta;}}_{\mathrm{\&lambda;NRateII}}\&times;(N_{\mathrm{RateII\&lambda;R}}\&times;T_{\mathrm{cntRateII}}-N_{\mathrm{RateI\&lambda;R}}\&times;T_{\mathrm{cntRateI}})}{({\mathrm{\&beta;}}_{\mathrm{\&lambda;NRateII}}\&times;T_{\mathrm{cntRateI}}-{\mathrm{\&beta;}}_{\mathrm{\&lambda;NRateI}}\&times;T_{\mathrm{cntrateII}})}$

式中：

$\left\{\begin{array}{c}{\mathrm{\&beta;}}_{\mathrm{\&lambda;NRateI}}=\frac{({\mathrm{\&lambda;}}_{\mathrm{Theo}-\mathrm{peak}-\mathrm{mR}}^{\mathrm{SZ}}-{\mathrm{\&lambda;}}_{\mathrm{Theo}-\mathrm{peak}-\mathrm{mL}}^{\mathrm{SZ}})}{(N_{\mathrm{RateI}-\mathrm{Theo}-\mathrm{peak}-\mathrm{mR}}^{\mathrm{SZ}}-N_{\mathrm{RateI}-\mathrm{Theo}-\mathrm{peak}-\mathrm{mL}}^{\mathrm{SZ}})}\\ {\mathrm{\&beta;}}_{\mathrm{\&lambda;NRateII}}=\frac{({\mathrm{\&lambda;}}_{\mathrm{Theo}-\mathrm{peak}-\mathrm{mR}}^{\mathrm{SZ}}-{\mathrm{\&lambda;}}_{\mathrm{Theo}-\mathrm{peak}-\mathrm{mL}}^{\mathrm{SZ}})}{(N_{\mathrm{RateII}-\mathrm{theo}-\mathrm{peak}-\mathrm{mR}}^{\mathrm{SZ}}-N_{\mathrm{RateII}-\mathrm{Theo}-\mathrm{peak}-\mathrm{mL}}^{\mathrm{SZ}})}\end{array}\right.$

$\left\{\begin{array}{c}{N}_{\mathrm{RateI\&lambda;R}}=N_{\mathrm{RateI}}^{\mathrm{fbg}}-N_{\mathrm{RateI}-\mathrm{Theo}-\mathrm{peak}-\mathrm{mR}}^{\mathrm{SZ}}\\ N_{\mathrm{RateII\&lambda;R}}=N_{\mathrm{RateII}}^{\mathrm{fbg}}-N_{\mathrm{RateII}-\mathrm{Theo}-\mathrm{peak}-\mathrm{mR}}^{\mathrm{SZ}}\end{array}\right.$

T_cntRateⅠ、T_cntRateⅡ分别表示RateⅠ、RateⅡ速率下峰值计数值的计数时钟周期；

将λ_Cali与

进行对比，同理检测出与λCali相邻的两个梳状滤波器波峰，编号分别为Cali-peak-mR、Cali-peak-mL，并记录下对应波峰计数值数据，其中波长较大的一个表示为：

${\mathrm{\&lambda;}}_{\mathrm{Cali}-\mathrm{peak}-mR}^{\mathrm{SZ}},({\mathrm{\&lambda;}}_{\mathrm{Cali}}<{\mathrm{\&lambda;}}_{\mathrm{Cali}-\mathrm{peak}-\mathrm{mR}}^{\mathrm{SZ}})$

称为λ_Cali“右侧梳状峰”；

波长较小的一个表示为：

${\mathrm{\&lambda;}}_{\mathrm{Cali}-\mathrm{peak}-\mathrm{mL}}^{\mathrm{SZ}},({\mathrm{\&lambda;}}_{\mathrm{Cali}}<{\mathrm{\&lambda;}}_{\mathrm{Cali}-\mathrm{peak}-\mathrm{mL}}^{\mathrm{SZ}})$

称为λ_Cali“左侧梳状峰”；

分别对应两个速率下的波峰计数值为：

$\left\{\begin{array}{cc}{N}_{\mathrm{RateI}-\mathrm{Cali}-\mathrm{peak}-\mathrm{mL}}^{\mathrm{SZ}},& N_{\mathrm{RateI}-\mathrm{Cali}-\mathrm{peak}-\mathrm{mR}}^{\mathrm{SZ}}\\ N_{\mathrm{RateII}-\mathrm{Cali}-\mathrm{peak}-\mathrm{mL},}^{\mathrm{SZ}}& N_{\mathrm{RateII}-\mathrm{Cali}-\mathrm{peak}-\mathrm{mR}}^{\mathrm{SZ}}\end{array}\right.$

S403、λ_Cali与λ_Theo的“左、右梳状峰”是否对应相同，若相同，则直接计算FBG延迟校准参数；若不相同，则需要再计算和比较，直到满足位置比较条件为止；

将λ_Cali的“左、右梳状峰”替换λ_Theo的“左、右梳状峰”，利用校准波长公式再次计算λ_Cali，后续步骤如上，最终得到满足位置比较条件的λ_Cali。，得到最终λ_Cali后，利用校准参数计算公式计算各FBG对应的延迟校准参数ΔN^fbg

$\mathrm{\&Delta;}{N}^{\mathrm{fbg}}=\frac{(N_{\mathrm{RateII\&lambda;R}}+\frac{({\mathrm{\&lambda;}}_{\mathrm{Cali}-\mathrm{peak}-\mathrm{mR}}^{\mathrm{SZ}}-{\mathrm{\&lambda;}}_{\mathrm{Cali}})}{{\mathrm{\&beta;}}_{\mathrm{\&lambda;NRateII}}})}{T_{\mathrm{cntRateII}}}$

8.根据权利要求1所述光纤光栅解调***中长距离光延迟效应的解调方法，其特征在于，步骤（5）包括：

在工作速率RateW下利用ΔN^fbg以及FBG波长顺序矫正后得到的检测通道中检测波峰排序与真实波长排序的对应关系，对FBG计数值

进行对应校准得到校准后计数值通过下式进行计数值校准：

$N_{\mathrm{RateW}-\mathrm{Cali}}^{\mathrm{fbg}}=N_{\mathrm{Ratew}}^{\mathrm{fbg}}-\mathrm{\&Delta;}{N}^{\mathrm{fbg}}$

然后，根据

利用现有波长标定技术标定得到FBG中心波长λ^fbg，即消除光信号传输延迟影响后的FBG真实中心波长值。

Images