CN114370889A - Ofdr系统中不同测量长度的自适应补偿测量方法 - Google Patents
Ofdr系统中不同测量长度的自适应补偿测量方法 Download PDFInfo
- Publication number
- CN114370889A CN114370889A CN202210015114.8A CN202210015114A CN114370889A CN 114370889 A CN114370889 A CN 114370889A CN 202210015114 A CN202210015114 A CN 202210015114A CN 114370889 A CN114370889 A CN 114370889A
- Authority
- CN
- China
- Prior art keywords
- coupler
- optical fiber
- frequency
- interferometer
- beat frequency
- 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.)
- Granted
Links
- 238000005259 measurement Methods 0.000 title claims abstract description 35
- 238000002168 optical frequency-domain reflectometry Methods 0.000 title claims abstract description 16
- 238000000691 measurement method Methods 0.000 title claims description 16
- 239000013307 optical fiber Substances 0.000 claims abstract description 66
- 238000000034 method Methods 0.000 claims abstract description 49
- 230000003044 adaptive effect Effects 0.000 claims abstract description 12
- 238000012937 correction Methods 0.000 claims abstract description 8
- 230000035559 beat frequency Effects 0.000 claims description 66
- 230000003287 optical effect Effects 0.000 claims description 17
- 239000000835 fiber Substances 0.000 claims description 11
- 238000005070 sampling Methods 0.000 claims description 11
- 238000001514 detection method Methods 0.000 claims description 4
- 238000012952 Resampling Methods 0.000 claims description 3
- 238000010408 sweeping Methods 0.000 claims description 2
- 230000000694 effects Effects 0.000 abstract description 8
- 239000000523 sample Substances 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 238000007796 conventional method Methods 0.000 description 3
- 230000001934 delay Effects 0.000 description 2
- 238000009795 derivation Methods 0.000 description 2
- 238000013519 translation Methods 0.000 description 2
- 230000014616 translation Effects 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 230000005672 electromagnetic field Effects 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
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/35338—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 other arrangements than interferometer arrangements
- G01D5/35354—Sensor working in reflection
- G01D5/35358—Sensor working in reflection using backscattering to detect the measured quantity
- G01D5/35361—Sensor working in reflection using backscattering to detect the measured quantity using elastic backscattering to detect the measured quantity, e.g. using Rayleigh backscattering
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Testing Of Optical Devices Or Fibers (AREA)
Abstract
本发明公开了一种OFDR系统中不同测量长度的自适应补偿测量方法，可以自适应校正激光器非线性调谐效应。自适应校正过程如下：首先测量出含有非线性信息的辅助干涉仪信号和主干涉仪信号，然后从辅助干涉仪信号提取出光源的非线性扫频曲线，根据需要测量的光纤反射点的不同位置，通过对扫频曲线进行相应的时移并与原扫频曲线作差，可以获取与不同位置相匹配的时延下的非线性扫频校正信号，从而可以对待测光纤不同位置的主干涉仪信号均进行最佳补偿，实现长距离、高空间分辨率的分布式光纤传感测量。
Description
技术领域
本发明涉及一种提高光频域反射计系统中测量范围和空间分辨率的方法，尤其涉及一种 OFDR系统中不同测量长度的自适应补偿测量方法，属于分布式光纤传感领域。
背景技术
光频域反射计(Optical Frequency Domain Reflectometer,OFDR)是分布式光纤传感中一项重要技术，最早由Eickhoff等人提出并证明。随着OFDR技术的成熟，它逐渐被应用到温度、应变、振动以及电磁场等参数的测量。在光频域反射系统中，可调谐光源(TLS)输出光被分为两束，一束为参考光，另一束为探测光。来自待测光纤菲涅尔反射光和后向散射光与参考光发生拍频干涉，产生拍频信号。待测点位置和拍频频率存在一一对应关系，通过傅里叶变换对拍频信号进行处理，就可以得到待测光纤上待测光纤的沿线反射点位置信息。
由于激光器的非线性扫频，导致拍频信号经过傅里叶变换后得到的频谱展宽，无法准确恢复出待测信息。因此对可调谐光源的非线性补偿是至关重要的。目前常见的非线性补偿方法有两种：硬件补偿法和软件补偿法。两种方法都是在原有的光路基础上添加一个相似结构的干涉仪，即辅助干涉仪，原来的测量光路为主干涉仪。硬件补偿法是辅助干涉仪输出的拍频信号作为采集卡的触发信号，实现对主干涉仪信号的等光频间隔采样，从而实现非线性补偿，但是最大测量距离只有辅助干涉仪两臂光路差的一半。而软件补偿法通过辅助干涉仪输出的拍频信号提取出可调谐光源的非线性信息，从而实现对主干涉仪拍频信号的补偿。传统的软件补偿方法有一维插值法、NUFFT、基于时间尺度因子补偿法、去斜滤波器法以及相位噪声补偿法等。一维插值法、NUFFT、基于时间尺度因子补偿法虽然降低了系统对辅助干涉仪两臂光路差的要求，但随着测量距离的增大，仍对辅助干涉仪两臂光路差有一定要求，测量距离受限。去斜滤波器法适用于长距离的测量，但是当非线性较大时，补偿效果不佳。
目前传统方法存在着测量距离受辅助干涉仪延迟光纤长度制约、空间分辨率差、运算速度慢等问题，难以实现在不同测量距离下仍保持高空间分辨率、高运算速度的测量。
发明内容
发明目的：为了解决上述问题，本发明提供一种OFDR系统中不同测量长度的自适应补偿测量方法，并且提高了系统的空间分辨率。
技术方案：本发明提出一种OFDR系统，由可调谐激光器、第一耦合器、主干涉仪、辅助干涉仪、第一和第二探测器、第一和第二采集卡以及计算机构成；所述主干涉仪由第二耦合器、第三耦合器、环形器及待测光纤组成，所述辅助干涉仪由第四耦合器、第五耦合器、短延迟光纤、声光移频器AOM组成。
所述第一耦合器的输入端连接可调谐激光器，所述第一耦合器的第一输出端连接第二耦合器的输入端，第一耦合器的第二输出端连接第四耦合器的输入端；所述第二耦合器的第一输出端连接环形器的第一端，第二耦合器的第二输出端连接第三耦合器的输入端，所述环形器的第二端连接第三耦合器的输入端，所述环形器的第三端连接待测光纤，所述第三耦合器的输出端连接第一探测器的输入端，所述第一探测器的输出端连接第一采集卡的输入端，第一采集卡的输出端连接计算机；所述第四耦合器的第一输出端连接第五耦合器的输入端，第四耦合器的第二输出端通过短延迟光纤连接AOM的输入端，所述AOM的输出端连接第五耦合器的输入端，所述第五耦合器的输出端连接第二探测器的输入端，所述第二探测器的输出端连接第二采集卡，所述第二采集卡的输出端连接计算机。
所述第一耦合器将可调谐激光器的输出分为两路，第一路进入主干涉仪后，经过第二耦合器分为探测光和参考光，探测光经过环形器注入待测光纤，待测光纤产生的后向瑞利散射光经过环形器进入第三耦合器，参考光直接进入第三耦合器与后向瑞利散射光发生拍频并产生第一拍频信号；第二路进入辅助干涉仪后，经过第四耦合器分为两束光，第一束光经过短延迟光纤和AOM后注入第五耦合器，第二束光直接注入第五耦合器与AOM的输出发生拍频并产生第二拍频信号；第一和第二拍频信分别经过第一和第二探测器转化为电信号，第一和第二采集卡分别对第一和第二探测器的输出进行采集，采集到的第一和第二信号输入计算机进行处理，以得到待测光纤的沿线反射点位置信息。
上述OFDR系统中不同测量长度的自适应补偿测量方法，包括以下步骤：
步骤一，测量在可调谐激光源非线性扫频影响下的辅助干涉仪拍频信号C(t)和主干涉仪拍频信号S(t)，其中辅助干涉仪短延迟光纤长度为Lc，主干涉仪待测光纤长度为Ls，分别对应两臂间时间延迟为τc、τs；
步骤二，求出辅助干涉仪拍频信号C(t)两两过零点之间的时间间隔，进而计算出拍频频率fc(t)；
步骤三，用步骤二中得到的拍频频率fc(t)，恢复出激光器的扫频速率γ(t)；
步骤四，将步骤三中得到的可调谐激光源的扫频速率γ(t)在时间上积分，得到激光器的扫频曲线ν(t)；
步骤五，将待测光纤分为N段，分别对每段按照以下方法进行补偿：；
步骤1，对待测光纤第i段，将扫频曲线ν(t)以该段中间位置Li产生的时延τi进行平移得到ν(t-τi)；其中，i＝1,2,…,N，N为待测光纤的段数；
步骤2，将ν(t)和ν(t-τi)作差值，得到新的拍频频率fci(t)；
步骤3，将新的拍频频率fci(t)代入辅助干涉仪信号表达式，得到新的辅助干涉仪信号Ci’(t)；
步骤4，采用新的辅助干涉仪信号Ci’(t)来补偿主干涉仪信号S(t)，得到待测光纤第i段沿线反射点的位置信息。
具体地，所述步骤二具体包括：采集卡的采样率为fs，辅助干涉仪拍频信号相邻两个过零点之间点数间隔结果记成数组d，点数间隔乘以采样间隔1/fs得到相邻过零点之间的时间间隔d/fs。将时间间隔d/fs乘以2再求倒数得到辅助干涉仪的拍频频率，减去AOM产生的移频频率fAOM得到辅助干涉仪真正的瞬时拍频频率fci(t)＝fs/2d-fAOM。
具体地，所述步骤三中具体包括：辅助干涉仪拍频频率是通过对实际扫频曲线ν(t)进行时移τc并与原扫频曲线作差得到，即fc(t)＝ν(t)-ν(t-τc)。为了获取扫频曲线，拍频频率公式两边都除以τc，得到fc(t)/τc＝[ν(t)-ν(t-τc)]/τc,因为辅助干涉仪两臂时延τc较小，表达式右边近似看作ν’(t)即扫频速率γ(t)，因此扫频速率γ(t)＝fc(t)/τc。
具体地，所述步骤六中：时延τi等于来回光纤长度2Li除以真空中光速c和光纤折射率n 即：τi＝2nLi/c。
具体地，所述步骤4中采用新的辅助干涉仪信号Ci’(t)来补偿主干涉仪信号S(t)，具体包括：该过程非线性补偿方法是基于时间尺度因子补偿法。具体步骤：求出新的辅助干涉仪信号Ci’(t)两两过零点，并根据过零点序号差得到校正系数，采用校正系数并结合重采样将主干涉仪采样点校正为等光频间隔点，从而实现非线性补偿。
与现有技术相比，本发明的技术方案具有以下有益技术效果：
(1)实现不同测量距离、高空间分辨率的测量
当辅助干涉仪两臂时延等于主干涉仪待测光纤产生的时延时，补偿效果达到最佳，本方法采用一个辅助干涉仪结构，可以产生两臂间不同时延的辅助干涉仪拍频信号，进而对不同测量距离的主干涉仪信号进行最佳补偿，实现不同测量距离、高空间分辨率的测量；
(2)恢复出准确的光源非线性扫频曲线
辅助干涉仪延迟臂采用短延迟光纤加AOM的结构，采用短延迟光纤可以提取出更准确的光源瞬时光频，AOM的引入可以增加提取出的光源瞬时光频点的数量，两者结合可以恢复出准确的光源非线性扫频曲线，从而有利于后期对主干涉仪信号的非线性补偿，提高测量空间分辨率；
(3)运算速度快
本发明所提出的非线性补偿方法中的算法简单，且所需要的辅助干涉仪延迟光纤短，算法数据量较小，所以本发明运算速度较于其他方法运算速度快。
附图说明
图1是一个实施例的光频域反射计结构示意图；
图2是一个实施例的光频域反射计中不同测量长度的自适应补偿测量方法的流程图；
图3是一个实施例的不同测量长度的自适应补偿测量方法的示意图，其中(a)是含有非线性信息的辅助干涉仪拍频信号和拍频频率时域图，(b)是拍频频率fc(t)恢复出激光器的扫频速率γ(t)过程,(c)是恢复出的可调谐激光源的扫频速率γ(t)和扫频曲线ν(t),(d)是求出新的辅助干涉仪拍频频率的过程；
图4是一个实施例的不同测量长度的自适应补偿测量方法补偿前后的结果，其中(a)是补偿前的结果，(b)是不同测量长度的自适应补偿测量方法补偿后的结果；
图5是一个实施例的采用传统方法加AOM和不加AOM补偿后的结果，其中(a)是加了20MHz的AOM结合传统补偿方法的补偿结果，(b)是无AOM的传统方法补偿结果。
具体实施方式
为了使本申请的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本申请进行进一步详细说明。应当理解，此处描述的具体实施例仅用以解释本申请，并不用于限定本申请。
在本文中提及“实施例”意味着，结合实施例描述的特定特征、结构或特性可以包含在本申请的至少一个实施例中。在说明书中的各个位置出现该短语并不一定均是指相同的实施例，也不是与其它实施例互斥的独立的或备选的实施例。本领域技术人员显式地和隐式地理解的是，本文所描述的实施例可以与其它实施例相结合。
图1为一个实施例的光频域反射计结构示意图。光频域反射计中不同测量长度的自适应补偿测量方法，所述光频域反射计系统由可调谐激光器、第一耦合器、主干涉仪、辅助干涉仪、第一和第二探测器、第一和第二采集卡以及计算机构成。主干涉仪由第二耦合器、第三耦合器、环形器及待测光纤组成，进入主干涉仪的扫频光经过第二耦合器分为探测光和参考光，探测光和参考光分别注入探测臂和参考臂，探测光经过环形器注入待测光纤，待测光纤产生的后向瑞利散射光经过环形器进入第三耦合器，参考光经过参考臂进入第三耦合器与后向瑞利散射光发生拍频并产生拍频信号S(t)。辅助干涉仪由第四耦合器、第五耦合器、短延迟光纤、声光移频器(AOM)组成，进入辅助干涉仪的扫频光经过第四耦合器分为两束，分别进入延迟臂和参考臂。注入延迟臂的光经过短延迟光纤和AOM移频到达第五耦合器，注入参考臂的光直接到达第五耦合器，两束光发生拍频干涉并产生拍频信号C(t)。拍频信号S(t) 和C(t)经过探测器转化为电信号，再通过采集卡进行数据采集，最后输入电脑进行信号处理。
图2为一个实施例的OFDR系统中不同测量长度的自适应补偿测量方法的流程图，具体包括如下步骤：
步骤一，测量在可调谐激光源非线性扫频影响下的辅助干涉仪拍频信号C(t)和主干涉仪拍频信号S(t)，其中辅助干涉仪短延迟光纤长度为Lc，主干涉仪待测光纤长度为Ls，分别对应两臂间时间延迟为τc、τs。
在一个实施例中，辅助干涉仪拍频信号和主干涉仪拍频信号理论推导过程如下：
由于可调谐激光源的扫频非线性的影响，扫频速率不再是定值，而是随时间随机波动的扫频速率γ(t)。扫频速率进行积分得到激光器的扫频曲线ν(t):
主干涉仪和辅助干涉仪随时间变化的拍频频率通过扫频曲线fs(t)和fc(t)经过平移作差得到：
fs(t)＝v(t)-v(t-τs) (2)
fc(t)＝v(t)-v(t-τc) (3)
其中，τs和τc分别为1.5*10-6s和1.5*10-7s，对应的主干涉仪中待测光纤长度为150m，辅助干涉仪延迟光纤长度为30m。通过对拍频频率进行积分得到拍频信号的相位，从而得到主干涉仪和辅助干涉仪拍频信号如下：
其中S0、C0分别为主干涉仪和辅助干涉仪拍频信号幅度。采集卡同时采集到辅助干涉仪拍频信号C(t)和主干涉仪拍频信号S(t)后，接下来采用不同测量长度的自适应补偿测量方法实现沿待测光纤的非线性补偿。
步骤二，求出辅助干涉仪拍频信号C(t)的两两过零点之间的间隔，进而计算出拍频频率 fc(t)；
该过程是根据辅助干涉仪拍频信号的过零点得到拍频频率，如图3中的(a)所示。辅助干涉仪的一个干涉臂是由短延迟光纤和声光移频器(AOM)构成。AOM用于提高输出拍频信号的过零点数量,为步骤四恢复可调谐激光源扫频曲线提供更多参考点，此时理论辅助干涉仪的拍频信号为：
其中C0为辅助干涉仪拍频信号的幅度，fAOM为移频频率，设置为20MHz。接着通过插值求出过零点的横坐标，并结合采样率求出两两相邻过零点之间的时间间隔tj(j＝1,2…k，k为总过零点间隔数)，将该时间间隔补充为一个完整的周期来求两两过零点之间的频率。减去AOM 的影响后，得到的辅助干涉仪拍频频率为：
fc(t)＝1/(2tj)-fAOM (7)
步骤三，用步骤二中拍频频率fc(t)恢复出激光器的扫频速率γ(t)；
在一个实施例中，如图3中的(b)所示，辅助干涉仪的拍频频率为经过时移τc后的实际扫频曲线与原实际扫频曲线之差。而对扫频曲线ν(t)的求导就是扫频速率γ(t)，扫频速率γ(t) 表达式如下：
其中Δt趋近于0。由于辅助干涉仪两臂间时延τc较小，式(9)可改为：
结合参考图3中的(c)以及式(10)可得，因为辅助干涉仪两臂间时延τc较小，即延迟臂上的延迟光纤较短，式(10)恢复出的激光器的扫频速率γ(t)较准确。
步骤四，将步骤三中的激光器的扫频速率γ(t)在时间上积分得到激光器的扫频曲线ν(t)，表达式为：
实际的采样点是离散的，因此上述的积分过程可改：
v(M)＝v(M+1)+γ(M-1)ΔT(M-1) (11)
其中M表示第M个过零点，△T是第M-1到第M个过零点之间的时间间隔。
步骤五，将待测光纤Ls分为N段进行补偿，每段中间位置分别为L1，L2，…Li…LN。其中i代表第i段光纤中间位置，i初始化的值为1；
步骤六，选取需要补偿的光纤位置Li并将扫频曲线ν(t)以该位置产生的时延τi进行平移得到ν(t-τi)；
进行非线性补偿的时候，当待测反射点位置在最佳补偿位置附近的时候，补偿效果最佳。最佳补偿位置即待测光纤某处产生的时延与辅助干涉仪两臂产生的时延相等。因此我们将待测光纤分为若干段，每段中间位置分别为L1，L2，…Li…Ls。选取第i段待测光纤段产生的时延为τi，接着对扫频曲线ν(t)以该位置产生的时延τi进行平移得到ν(t-τi)，与原扫频曲线作差便可得到以τi为时延的辅助干涉仪拍频频率，该过程如图3中的(d)所示。此时新的辅助干涉仪两臂时延与主干涉仪待测位置Li产生的时延相等，可以实现最佳补偿。因此我们通过选取不同段待测光纤，并求出对应的时延下的辅助干涉仪拍频频率，便可实现对待测光纤任意位置的最佳补偿。
步骤七，将步骤六中两条平移前后的扫频曲线作差值得到新的拍频频率fci(t)；
步骤八，将新求出的拍频频率fci(t)带入辅助干涉仪信号表达式，得到新的辅助干涉仪信号Ci’(t)；
将以时延τi时移后的扫频曲线与原扫频曲线作差得到新的辅助干涉仪拍频频率。以τi为时延的辅助干涉仪拍频频率表达式为：
fci(t)＝v(t)-v(t-τi) (12)
将新求出的拍频频率带入表达式(6)，得到新的辅助干涉仪拍频信号。新辅助干涉仪两臂时延等于所选取主干涉仪待测光纤第i段中间位置产生的时延，从而实现该段最佳补偿。
步骤九，基于时间尺度因子补偿法，采用新的辅助干涉仪信号Ci’(t)来补偿主干涉仪信号 S(t)，得到待测光纤第i段沿线反射点的位置信息。接着i＝i+1，重复步骤六到步骤九直至补偿整个待测光纤。
基于时间尺度因子补偿法(J Xing,Y Zhang,F Wang,et al.A Method Based onTime-Scale Factor for Correcting the Nonlinear Frequency Sweeping in an OFDRSystem[J].Photonics Journal, IEEE,2019,11(3):1-8)，具体步骤：求出新的辅助干涉仪信号Ci’(t)两两过零点，并根据过零点序号差得到校正系数，采用校正系数并结合重采样将主干涉仪采样点校正为等光频间隔点，从而实现非线性补偿。
下面进行实验验证。主干涉仪的待测光纤长度为150m，对应的两臂时延τs＝1.5*10-6s，辅助干涉仪中延迟光纤长度为30m，对应的两臂时延为τc＝1.5*10-7s。激光器的扫频速率为 70nm/s，非线性扰动项幅值为5。采样率为312.5MSa/s，采样点数为40M。AOM的频率为 20MHz,系统的理论空间分辨率为0.09mm。
将主干涉仪待测光纤均匀分为5段，每段中间位置依次为15m、45m、75m、105m以及135m，接着采用本发明提出的不同测量长度的自适应补偿测量方法进行补偿。参考图4是采用本发明所提出的方法补偿前后的结果，其中辅助干涉仪延迟光纤长度为30m，远小于主干涉仪中待测光纤的长度。相比于图4中的(a)，经过该方法补偿后，如图4中的(b)所示，空间分辨率达到0.14mm，接近理论空间分辨率的值。
图5是在相同参数设置下，采用传统非线性补偿方法的对比结果。图5中的(a)是加了20 MHz的AOM结合传统补偿方法的补偿结果，其中AOM的作用是为补偿提供更多参考点。图5中的(b)是无AOM的传统方法补偿结果。可以看到由于辅助干涉仪延迟光纤较短，即使加了AOM，补偿效果仍不佳，空间分辨率为60mm，远大于理论空间分辨率的值。
以上的实验结果说明：本发明所提出的可以实现较好的补偿效果，空间分辨率接近理论值。同时和传统方法相比，本发明所提出的方法可以实现当辅助干涉仪延迟光纤较短时，通过准确恢复出扫频曲线后，对扫频曲线进行时移并与原扫频曲线作差可以得到任意延迟光纤长度下的辅助干涉仪拍频信号，从而实现任意长度的主干涉仪信号的补偿，实现不同测量距离、高空间分辨率的测量。而传统方法测量长度受限于辅助干涉仪延迟光纤长度，在相同参数设置下，传统方法的补偿效果要远差于本发明提出的非线性补偿方法。
Claims (7)
1.一种OFDR系统，其特征在于，所述OFDR系统由可调谐激光器、第一耦合器、主干涉仪、辅助干涉仪、第一和第二探测器、第一和第二采集卡以及计算机构成；所述主干涉仪由第二耦合器、第三耦合器、环形器及待测光纤组成，所述辅助干涉仪由第四耦合器、第五耦合器、短延迟光纤、声光移频器AOM组成；
所述第一耦合器的输入端连接可调谐激光器，所述第一耦合器的第一输出端连接第二耦合器的输入端，第一耦合器的第二输出端连接第四耦合器的输入端；所述第二耦合器的第一输出端连接环形器的第一端，第二耦合器的第二输出端连接第三耦合器的输入端，所述环形器的第二端连接第三耦合器的输入端，所述环形器的第三端连接待测光纤，所述第三耦合器的输出端连接第一探测器的输入端，所述第一探测器的输出端连接第一采集卡的输入端，第一采集卡的输出端连接计算机；所述第四耦合器的第一输出端连接第五耦合器的输入端，第四耦合器的第二输出端通过短延迟光纤连接AOM的输入端，所述AOM的输出端连接第五耦合器的输入端，所述第五耦合器的输出端连接第二探测器的输入端，所述第二探测器的输出端连接第二采集卡，所述第二采集卡的输出端连接计算机；
所述第一耦合器将可调谐激光器的输出分为两路，第一路进入主干涉仪后，经过第二耦合器分为探测光和参考光，探测光经过环形器注入待测光纤，待测光纤产生的后向瑞利散射光经过环形器进入第三耦合器，参考光直接进入第三耦合器与后向瑞利散射光发生拍频并产生第一拍频信号；第二路进入辅助干涉仪后，经过第四耦合器分为两束光，第一束光经过短延迟光纤和AOM后注入第五耦合器，第二束光直接注入第五耦合器与AOM的输出发生拍频并产生第二拍频信号；第一和第二拍频信分别经过第一和第二探测器转化为电信号，第一和第二采集卡分别对第一和第二探测器的输出进行采集，采集到的第一和第二信号输入计算机进行处理，以得到待测光纤的沿线反射点位置信息。
2.根据权利要求1所述的一种OFDR系统，其特征在于，所述短延迟光纤的长度小于待测光纤长度的一半。
3.根据权利要求1或2所述的系统实现的不同测量长度的自适应补偿测量方法，其特征在于，包括以下步骤：
步骤一，测量在可调谐激光源非线性扫频影响下的第二拍频信号C(t)和第一拍频信号S(t)；
步骤二，求出第二拍频信号C(t)两两过零点之间的时间间隔，进而计算出拍频频率fc(t)；
步骤三，用步骤二中得到的拍频频率fc(t)，恢复出可调谐激光源的扫频速率γ(t)；
步骤四，将步骤三中得到的扫频速率γ(t)在时间上积分，得到激光器的扫频曲线ν(t)；
步骤五，将待测光纤分为N段，分别对每段按照以下方法进行补偿：
步骤1，对待测光纤第i段，将扫频曲线ν(t)以该段中间位置Li产生的时延τi进行平移得到ν(t-τi)；其中，i＝1,2,…,N，N为待测光纤的段数；
步骤2，将ν(t)和ν(t-τi)作差值，得到新的拍频频率fci(t)；
步骤3，根据新的拍频频率fci(t)，得到新的辅助干涉仪信号Ci’(t)；
步骤4，采用新的辅助干涉仪信号Ci’(t)来补偿主干涉仪信号S(t)，得到待测光纤第i段沿线反射点的位置信息。
4.根据权利要求1所述的不同测量长度的自适应补偿测量方法，其特征在于，所述步骤二中fc(t)＝fs/2d-fAOM，其中fs为第二采集卡的采样率，d为第二拍频信号C(t)相邻两个过零点之间时间间隔构成的数组，fAOM为AOM产生的移频频率。
5.根据权利要求1所述的不同测量长度的自适应补偿测量方法，其特征在于，所述步骤三‘中γ(t)＝fc(t)/τc，其中τc为第四耦合器输出的两束光进入辅助干涉仪后发生的时移。
6.根据权利要求1所述的不同测量长度的自适应补偿测量方法，其特征在于，所述步骤1中τi＝2nLi/c，其中n为光纤折射率，c为真空中光速。
7.根据权利要求1所述的不同测量长度的自适应补偿测量方法，其特征在于，所述步骤4中采用基于时间尺度因子补偿法进行补偿，具体步骤是：求出新的辅助干涉仪信号Ci’(t)两两过零点，并根据过零点序号差得到校正系数，采用校正系数并结合重采样将主干涉仪采样点校正为等光频间隔点，从而实现非线性补偿。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210015114.8A CN114370889B (zh) | 2022-01-07 | 2022-01-07 | Ofdr系统中不同测量长度的自适应补偿测量方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210015114.8A CN114370889B (zh) | 2022-01-07 | 2022-01-07 | Ofdr系统中不同测量长度的自适应补偿测量方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114370889A true CN114370889A (zh) | 2022-04-19 |
| CN114370889B CN114370889B (zh) | 2023-03-14 |
Family
ID=81141479
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210015114.8A Active CN114370889B (zh) | 2022-01-07 | 2022-01-07 | Ofdr系统中不同测量长度的自适应补偿测量方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114370889B (zh) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115127778A (zh) * | 2022-05-25 | 2022-09-30 | 深圳大学 | 一种简化ofdr系统的方法和装置以及ofdr简化系统、控制装置 |
| CN115327515A (zh) * | 2022-08-10 | 2022-11-11 | 哈尔滨工业大学 | 一种基于相位传递的双扫频干涉动态测量系统及测量方法 |
| CN116659395A (zh) * | 2023-06-07 | 2023-08-29 | 哈尔滨工业大学 | 一种基于相位补偿的扫频干涉测量非线性同步误差校正方法 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102997937A (zh) * | 2012-12-12 | 2013-03-27 | 天津大学 | 一种可抑制光源相位噪声的光频域反射装置和解调方法 |
| US20140140691A1 (en) * | 2011-06-28 | 2014-05-22 | Luna Innovations Incorporated | Fiber optic network interrogation tool for combined swept-heterodyne optical spectrum analysis and optical frequency-domain reflectometry |
| US20150346053A1 (en) * | 2012-12-24 | 2015-12-03 | Intuitive Surgical Operations, Inc. | Dispersion correction in optical frequency-domain reflectometry |
| CN111397644A (zh) * | 2020-03-26 | 2020-07-10 | 南京大学 | 一种用于光频域反射计的激光器非线性调谐效应补偿系统及补偿方法 |
| CN111578971A (zh) * | 2020-05-13 | 2020-08-25 | 武汉昊衡科技有限公司 | Ofdr分段采集实现长距离测量的装置和方法 |
| CN112923960A (zh) * | 2021-02-01 | 2021-06-08 | 南京大学 | 用于校正非线性调谐效应的光纤参数测量装置 |
-
2022
- 2022-01-07 CN CN202210015114.8A patent/CN114370889B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20140140691A1 (en) * | 2011-06-28 | 2014-05-22 | Luna Innovations Incorporated | Fiber optic network interrogation tool for combined swept-heterodyne optical spectrum analysis and optical frequency-domain reflectometry |
| CN102997937A (zh) * | 2012-12-12 | 2013-03-27 | 天津大学 | 一种可抑制光源相位噪声的光频域反射装置和解调方法 |
| US20150346053A1 (en) * | 2012-12-24 | 2015-12-03 | Intuitive Surgical Operations, Inc. | Dispersion correction in optical frequency-domain reflectometry |
| CN111397644A (zh) * | 2020-03-26 | 2020-07-10 | 南京大学 | 一种用于光频域反射计的激光器非线性调谐效应补偿系统及补偿方法 |
| CN111578971A (zh) * | 2020-05-13 | 2020-08-25 | 武汉昊衡科技有限公司 | Ofdr分段采集实现长距离测量的装置和方法 |
| CN112923960A (zh) * | 2021-02-01 | 2021-06-08 | 南京大学 | 用于校正非线性调谐效应的光纤参数测量装置 |
Non-Patent Citations (1)
| Title |
|---|
| FENG WANG，ET AL: "Enhancing the Effect of Nonlinear Frequency Sweep Correction in OFDR With Improved Reference Frequency", 《JOURNAL OF LIGHTWAVE TECHNOLOGY》 * |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115127778A (zh) * | 2022-05-25 | 2022-09-30 | 深圳大学 | 一种简化ofdr系统的方法和装置以及ofdr简化系统、控制装置 |
| CN115127778B (zh) * | 2022-05-25 | 2023-04-18 | 深圳大学 | 一种简化ofdr系统的方法和装置以及ofdr简化系统、控制装置 |
| CN115327515A (zh) * | 2022-08-10 | 2022-11-11 | 哈尔滨工业大学 | 一种基于相位传递的双扫频干涉动态测量系统及测量方法 |
| CN115327515B (zh) * | 2022-08-10 | 2023-04-18 | 哈尔滨工业大学 | 一种基于相位传递的双扫频干涉动态测量系统及测量方法 |
| CN116659395A (zh) * | 2023-06-07 | 2023-08-29 | 哈尔滨工业大学 | 一种基于相位补偿的扫频干涉测量非线性同步误差校正方法 |
| CN116659395B (zh) * | 2023-06-07 | 2024-02-23 | 哈尔滨工业大学 | 一种基于相位补偿的扫频干涉测量非线性同步误差校正方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114370889B (zh) | 2023-03-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN114370889B (zh) | Ofdr系统中不同测量长度的自适应补偿测量方法 | |
| CN102420650B (zh) | 一种抑制光频域反射仪的激光器非线性扫描装置和方法 | |
| CN112923960B (zh) | 用于校正非线性调谐效应的光纤参数测量装置 | |
| JP6236369B2 (ja) | 光ファイバの温度・歪み分布測定方法および装置 | |
| JP6647420B2 (ja) | ブリルアン散乱測定方法およびブリルアン散乱測定装置 | |
| CN111397644B (zh) | 一种用于光频域反射计的激光器非线性调谐效应的补偿方法 | |
| CN113804404B (zh) | 一种用于光频域偏振串音测量的光源扫频非线性校正方法 | |
| US11920974B2 (en) | Vibration distribution measuring system, vibration waveform analysis method, vibration waveform analyzing device, and analyzing program | |
| CN103090808A (zh) | 一种基于光谱位相的高精度大量程间距测量方法及系统 | |
| Xing et al. | A method based on time-scale factor for correcting the nonlinear frequency sweeping in an OFDR system | |
| CN103944644B (zh) | 基于匹配延迟采样光相干系统 | |
| Xiao et al. | Frequency response enhancement of Φ-OTDR using interval-sweeping pulse equivalent sampling based on compressed sensing | |
| CN113804301A (zh) | 一种基于光频域移频干涉的分布式偏振串音快速测量装置 | |
| CN104038281B (zh) | 非线性相位估计的长距离高分辨率光频域反射解调方法 | |
| CN110749420B (zh) | 一种ofdr检测装置 | |
| CN103763021B (zh) | 一种相干光时域反射测量方法及反射仪装置 | |
| Wang et al. | Enhancing the effect of nonlinear frequency sweep correction in OFDR with improved reference frequency | |
| Li et al. | Long-distance OFDR with high spatial resolution based on time division multiplexing | |
| Li et al. | Phase extraction of optical carrier-based microwave interferometry with all-phase fast Fourier transform for distance measurement | |
| CN106705863B (zh) | 一种提高光频域反射仪的最大测试距离的方法 | |
| JPH03120437A (ja) | 光導波路の歪または温度の測定方法および装置 | |
| CN115711633A (zh) | 一种环路结构参考干涉仪的相位噪声精确校正光频域反射计 | |
| CN117480364A (zh) | 振动测量器及振动测量方法 | |
| JP3759798B2 (ja) | 落雷点標定方法 | |
| CN112798025A (zh) | 提高ofdr测量空间分辨率的方法及ofdr系统 |
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 |