CN113670350B - 高精度光频域反射方法、装置及系统 - Google Patents

高精度光频域反射方法、装置及系统 Download PDF

Info

Publication number
CN113670350B
CN113670350B CN202010415101.0A CN202010415101A CN113670350B CN 113670350 B CN113670350 B CN 113670350B CN 202010415101 A CN202010415101 A CN 202010415101A CN 113670350 B CN113670350 B CN 113670350B
Authority
CN
China
Prior art keywords
signal
wavelength
laser
main path
output
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
CN202010415101.0A
Other languages
English (en)
Other versions
CN113670350A (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.)
Dalian University of Technology
Original Assignee
Dalian University of Technology
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 Dalian University of Technology filed Critical Dalian University of Technology
Priority to CN202010415101.0A priority Critical patent/CN113670350B/zh
Publication of CN113670350A publication Critical patent/CN113670350A/zh
Application granted granted Critical
Publication of CN113670350B publication Critical patent/CN113670350B/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
    • 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/35325—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 interferometer with two arms in reflection, e.g. Mickelson interferometer
    • 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/35329—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 interferometer with two arms in transmission, e.g. Mach-Zender interferometer

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Instruments For Measurement Of Length By Optical Means (AREA)
  • Length Measuring Devices By Optical Means (AREA)

Abstract

本申请公开了一种高精度光频域反射方法，用于通过耦合于待测主体的传感光纤对所述待测主体的物理量变化进行测量，其将垂直腔面反射激光器应用于基于光频域反射技术的分布式物理量测量装置中，实现了大范围无跳模波长调谐范围，提高了分布式测量方法和装置的空间分辨力和测量量程。本申请还公开了与该方法相关联的装置和系统。

Description

高精度光频域反射方法、装置及系统
技术领域
本发明属于光纤传感技术领域，尤其涉及一种高精度光频域反射方法、装置及系统。
背景技术
光频域反射技术原理实现的分布式物理量测量是一种可以高空间分辨力实现物理量分布式测量的一种技术手段，相关的较早的文献包括：
Distributed measurement of static strain in an optical fiber withmultiple Bragg gratings at nominally equal wavelengths[J].Applied Optics,1998,37(10):1741-1746.
High-spatial-resolution distributed strain measurement in opticalfiber with Rayleigh scatter[J].Applied Optics,1998,37(10):1735-1740.
基于光频域反射原理的分布式传感系统中使用可调谐激光器作为系统光源，空间分辨力及物理量测量量程受限于可调谐激光器输出光信号的调制范围。系统空间分辨力在数值上反比于可调谐激光器输出光信号的调制范围且系统空间分辨力的极限取决于可调谐激光器输出光信号的调制范围。调谐范围越大，系统空间分辨力越高。同时，调谐范围越大，被测物理量如应变温度等扰动的测量量程也越大。可见，要想提高光纤传感的空间分辨力和物理量测量量程，需要提高可调谐激光器的扫频或调谐范围。此外，系统测量时间为调谐光源完成一次全波段调谐的时间，因此系统测量速度直接取决于可调谐激光器波长调谐速度。
现有技术或装置中，通常采用外腔式调谐激光器或半导体激光器作为光源，但是外腔式调谐激光器成本昂贵且容易出现跳模。同时，外腔激光器通常采用littrow或者littman结构，为机械方式实现的波长调谐，因此调谐速度慢，通常为几十nm/s，影响了在一些要求高速测量的场合下的测量速度。
基于垂直腔面反射激光器(VCSEL)是半导体激光器的一种，垂直腔面反射激光器的波长可随电流或温度连续调谐,且温度和电流调谐的灵敏度都很高。同时相较于分布式反馈激光器(DFB)调谐范围更大。以韩国Raycan公司的垂直腔面反射激光器RC32xxx1-PFAmt为例，温度每变化1摄氏度，波长变换0.1nm左右，调谐范围一般低于5nm。另一方面，激光器驱动电流调谐实现的波长调谐，通常在标称安全电流以下每mA对应波长变换量为0.5nm左右(Altabas J A,Izquierdo D,Lazaro J A,et al.Chirp-based direct phasemodulation of VCSELs for cost-effective transceivers[J].Optics letters,2017,42(3):583-586.)。且电流调谐下其调制速度可以上kHz甚至几十kHz。
本发明将垂直腔面反射激光器应用于基于光频域反射技术的分布式传感系统中，研究一种高精度光频域反射方法及装置。
发明内容
本发明的一些实施例提供了一种高精度光频域反射方法，用于通过耦合于待测主体的传感光纤对所述待测主体的物理量变化进行测量，该方法包括以下步骤：对至少两个垂直腔面反射激光器中的每一个调整至不同的工作温度以使得所述各垂直腔面反射激光器的激光输出具有不同的输出波长；对每个所述垂直腔面反射激光器进行电流调谐以提供调谐激光输出，其中所述调谐激光输出包括所述每个垂直腔面反射激光器提供的调谐激光输出波段，其中波长相邻的所述调谐激光输出波段有重叠区域；在包括所述物理量变化的测量态下，同步地采集包括所述调谐激光输出的绝对波长的测量态波长监测信号；以及响应所述调谐激光输出于所述传感光纤对调谐激光输出的干涉的测量态主路干涉信号；根据所述测量态波长监测信号确定各个垂直腔面反射激光器输出的测量态信号拼接位置，根据该位置对所述测量态主路干涉信号拼接得到拼接后测量态主路干涉信号；以及
基于所述拼接后测量态干涉光信号和在无所述物理量变化的参考态获得的拼接后参考态干涉光信号解算所述物理量变化。
本发明的另一些实施例提供了高精度光频域反射装置，用于通过耦合于待测主体的传感光纤对所述待测主体的物理量变化进行测量，该装置包括：至少两个垂直腔面反射激光器，所述每一个垂直腔面反射激光器被配置为作以下调谐：调整至不同的工作温度以使得所述各垂直腔面反射激光器的激光输出具有不同的输出波长；以及对每个所述垂直腔面反射激光器进行电流调谐以提供调谐激光输出，其中所述调谐激光输出包括所述每个垂直腔面反射激光器提供的调谐激光输出波段，其中波长相邻的所述调谐激光输出波段有重叠区域；电流驱动单元，用于为所述垂直腔面反射激光器提供驱动电流信号，进而对垂直腔面反射激光器进行波长调谐；温度控制单元，用于为所述垂直腔面反射激光器中的垂直腔面反射激光器提供的温度控制信号；波长监测单元，用于根据波长监测信号中的光谱重叠部分来确定相邻波段的波长拼接位置；主路干涉仪单元，用于测量传感光纤上的散射或者反射信号；采集单元，用于采集激光器波长监测单元输出的信号或者其读数，用于采集辅助干涉仪信号，用于采集主路干涉仪信号；控制单元，用于上述光源调谐电流驱动单元、温度控制单元、采集单元的时序控制；数据处理单元，配置为对采集的数据进行运算，包括确定各路信号的拼接位置；对主路干涉仪信号进行截取和拼接得到拼接后主路干涉仪信号；对所述主路干涉仪信号、波长监测单元输出的信号进行非线性校正；以及对分布式物理量进行解算。。
本申请还提供了基于上述装置和/或方法的分布式测量系统。
本发明的有益效果：本发明将垂直腔面反射激光器应用于基于光频域反射技术的分布式物理量测量装置中，实现了大范围无跳模波长调谐范围，提高了分布式测量方法和装置的空间分辨力和测量量程。同时，测量速度可以达到kHz以上，可以将原本受限于系统光源调谐速度只能工作在静态或准静态测量的分布式物理量测量系统可以应用于动态测量场合。该方法和装置具有控制方法和装置简单，成本低，且体积小，利于系统集成的优点。
附图说明
图1为根据本申请的实施例的一种垂直腔面反射激光器结构示意图；
图2为根据本申请的实施例的高精度光频域反射装置；
图3为根据本申请的实施例的半导体光放大器及其闭环光功率控制示意图；
图4为根据本申请的实施例的氰化氢气室吸收谱；
图5为根据本申请的实施例的拼接点确定示意图；
图6为根据本申请的实施例的光纤干涉仪信号；
图1中：44为热电冷却器，45为热敏电阻，43为垂直腔面反射激光器，46为激光二极管。
图2中：28为第一激光器，29为第二激光器，30为第五光纤耦合器，24为第一电流驱动单元，25为第一温度控制单元，26为第二电流驱动单元，27为第二温度控制单元，3为第一光纤耦合器，4为第二光纤耦合器，12为第三光纤耦合器，7为第四光纤耦合器，23为延时光纤，20为波长监测单元，21为辅助干涉仪单元，22为主路干涉仪单元，19为采集单元，18为第一光电探测器，15为第二光电探测器，10为第三光电探测器，31为存储单元，32为数据处理单元，13为第一法拉第旋转镜，14为第二法拉第旋转镜，6为传感光纤，39为外界作用拉伸应变，2为闭环功率控制单元，5为光纤环形器，1为控制单元。
图3中，91为半导体光放大器，92为第十耦合器，93为第四光电探测器，94为比较器。95为滤波器，96为运算放大器，97为电流驱动器。
图5中，70为第一主路干涉仪信号拼接位置，71为第二主路干涉仪信号拼接位置，72为第一主路干涉仪信号，73为第二主路干涉仪信号，74为第一波长监测单元信号，75为第二波长监测单元信号，77为拼接后主路干涉仪输出光信号。
具体实施方式
本发明将垂直腔面反射激光器应用于基于光频域反射技术的分布式传感系统中，研究一种高精度光频域反射方法及装置。本专利所提出的高精度光频域反射方法及装置将以垂直腔面反射激光器作为系统光源。在本发明中将采用驱动电流调谐来控制激光器的输出波长的调谐，通过将多个垂直腔面反射激光器配置在不同温度来使得各个激光器本身具有一定波长间隔，然后在特定驱动电流下依次对各个垂直腔面反射激光器进行波长调谐。对各个激光器输出激光在波长上进行复用和拼接，相当于扩大了分布式物理量测量装置光源的扫频或波长调谐范围。同时由于驱动电流的大小同时影响激光器输出的光功率，因为为了获得稳定功率的激光输出，在激光器后面增加以半导体光放大器为核心的闭环功率控制单元。此外，利用波长监测单元来确定相邻波长激光器调谐过程的拼接位置，使得各个垂直腔面反射激光器的波长可以实现无间隔覆盖。在获取了全波段的激光输出后，将该激光注入到分布式物理量测量系统中，通过考察参考态和测量态的相对光谱移动进而实现高速高分辨力的分布式物理量测量。下面将就本发明的具体细节进行阐述。
图1为一种典型垂直腔面反射激光器结构示意图。通常情况下，垂直腔面反射激光器43上具有可以用电流控制的用于加热或者制冷的热电冷却器44以及阻值随温度变化而变化的热敏电阻45。以韩国Raycan公司的RC32xxx1-PFAmt型号垂直腔面反射激光器来说，垂直腔面反射激光器输出波长对温度和电流均有所响应。其阈值电流为2mA，最大工作电流为14mA,在电流从2mA变化的14mA下，其输出波长变化约为6nm。若为同一生产批次，则激光器在某一驱动电流和温度，如7mA，15度下的输出中心波长应相等或偏差小于0.5nm。在固定驱动电流下，温度值从10度改变到50度，其输出波长变化约为5nm。不失一般性，下面以该垂直腔面反射激光器及其参数作为一种高精度光频域反射装置的光源来说明如何实现传感。
测量装置中需要有两个或两个以上的垂直腔面反射激光器，以两个垂直腔面反射激光器复用作为测量装置的光源为例进行阐述。假设这两个垂直腔面反射激光器在2mA，10度下具有相同的输出波长，如1530nm。首先通过外围控制电路包括电流驱动和温度驱动来改变激光器驱动电流和温度，并通过激光器波长监测单元如光谱仪来确定激光器在10度下，驱动电流从2mA变化到14mA下输出波长变化量，一种可能的情况为：驱动电流从2mA变化到14mA下输出波长从1530nm变化到1536nm。然后，在驱动电流为2mA下，从10度逐渐增加激光器温度，同时观察光谱仪读数，记录激光器输出波长值。在波长值为1535nm时停止温度增长，记录此时温度。一种可能的情况为：此时温度为50度。重新配置两个激光器，两个激光器的驱动电流均设置为2mA，温度为10度的激光器被记为1号激光器，温度为50度的激光器被记为2号激光器。
图2中控制单元1控制着四路共两组信号，分别是针对第一激光器28的第一电流驱动单元24和第一温度控制单元25，以及针对第二激光器29的第二电流驱动单元26和第二温度控制单元27。控制单元1控制第一温度控制单元25使其控制第一激光器28温度为10度，控制单元1控制第二温度控制单元26使其控制第二激光器温度为50度。其后，控制单元1控制第二电流驱动单元24为第二激光器施加阈值电流之下的恒定电流，如0或者1mA，此时第二激光器无激光输出。然后控制单元1控制第一电流驱动单元26为第一激光器施加阈值电流之上的起始驱动电流2mA。然后控制单元1控制第一电流驱动单元26为第一激光器施加起始驱动电流2mA到终止驱动电流14mA的电流信号，该电流信号可以是斜波信号，也就是线性增加的电流信号，完成一次驱动电流从2mA到14mA的时间可以为10微秒。此时第一激光器28输出激光的波长从1530到1536nm的波长调谐激光。该激光经过第五光纤耦合器30进入闭环功率控制单元2，其后进入到第一光纤耦合器3，激光在第一光纤耦合器3被分成三束输出光分别进入到波长监测单元20，辅助干涉仪单元21与主路干涉仪单元22。波长监测单元20包括可以输出特征信号的氰化氢分子气室，其吸收谱如图4所示，在特定可溯源的波长位置被吸收，经过氰化氢分子气室的光被第一光电探测器18探测并光电转换被采集单元19的一个采集通道采集并传入到存储单元31。采集单元19可以为多通道示波器或采集卡。同时，一部分光经过第一光纤耦合器3的另一个输出端口c输出到辅助干涉仪单元21，辅助干涉仪单元21可以为图2所示的迈克尔孙干涉仪结构组成，由第三光纤耦合器12，延时光纤23，第一法拉第旋转镜13，第二法拉第旋转镜14，第二光电探测器15构成。对于该图所示的迈克尔孙结构干涉仪而言，在激光器调谐过程下干涉仪输出的信号为正弦信号，正弦的周期与扫过的波长范围有关且与延时光纤23长度有关。延时光纤23长度越长，正弦的周期越小，每个正弦对应扫过的波长范围越小。同时该正弦信号与光源输出信号的相位直接对应，因此该干涉仪可以对调谐激光器输出的激光进行波长或者相位追迹，可以用于后续的非线性校正。辅助干涉仪单元21输出的正弦信号被采集单元19采集。第一光纤耦合器3的另一个输出端口d输出光到主路干涉仪单元22，主路干涉仪也是完成测量或者传感的光路结构，包括具有马赫曾德结构或者迈克尔孙结构的光纤干涉仪，用于测量传感光纤上的散射或者反射信号。图2所示的为一个马赫曾德结构的光纤干涉仪，第二光纤耦合器4分出一束参考臂直接连接到第四光纤耦合器7，另一路为测量臂，从第二光纤耦合器4输出到光纤环形器5的a端口，光纤环形器5特性为a进c出，c进b出。因此，从光纤环形器5的a端口进去的光进入到测量传感光纤6，传感光纤6上的散射或者反射信号返回到光纤环形器5并进入到光纤环形器5的b端口进而进入到第四光纤耦合器7。两束光在光纤耦合器7中合束后被第三光电探测器10探测并被采集单元19采集。三路信号被采集单元19采集后传到存储单元31进而传入数据处理单元32。上述三路信号的原始数据分别记为第一波长监测数据，第一辅助干涉仪数据，第一主路干涉仪数据，第一指的是第一激光器。
接下来，其后，控制单元1控制第一电流驱动单元24为第一激光器施加阈值电流之下的恒定电流，如0或者1mA，此时第一激光器28无激光输出。然后控制单元1控制第二电流驱动单元26为第二激光器28施加阈值电流之上的起始驱动电流2mA。然后控制单元1控制第二电流驱动单元26为第二激光器29施加起始驱动电流2mA到终止驱动电流14mA的电流信号，此时第二激光器29输出激光的波长从1535到1541nm的波长调谐激光。，其他采集过程同上段所述。三路信号被采集单元19采集后传到存储单元31进而传入数据处理单元32。上述三路信号的原始数据分别记为第二波长监测数据，第二辅助干涉仪数据，第二主路干涉仪数据，第二指的是第二激光器。
下面将两个垂直腔面反射激光器调谐下的原始数据进行拼接。
为了在波长重叠区域内确定一个拼接位置以便对辅助干涉仪和主路干涉仪信号进行截取与拼接来得到无重叠的连续的输出信号，需要对垂直腔面反射激光器调谐过程的波长进行监测或者追迹。在下面介绍利用波长监测单元20来确定辅助干涉仪信号和主路干涉仪信号拼接位置的过程。
波长监测单元20其本身可以是直接测量波长的装置如光谱仪或者波长计，其读数即为激光波长，此时只需要选取重合区域中任意一点即可，优选地，可以选择重合区域正中的一点作为拼接位置。
波长监测单元20也可以是可以表征波长特征或变化的气体分子气室如氰化氢分子气室或者乙炔分子气室。下面用图6说明以氰化氢气室作为绝对波长监测单元20来确定拼接点位置，并对相邻两段直接光强探测单元信号进行截取与拼接的过程。如图5为氰化氢分子气室特征谱线，在输入光信号进行波长调谐中，氰化氢分子气室的透射光具有图5的吸收谱线。在第一激光器28调谐范围1530-1536nm和第二激光器29调谐范围1535-1541nm，两段波长重合的部分为1535-1536nm，因此可以用落在重叠区域的任意一个氰化氢分子气室吸收峰作为拼接位置，在1535-1536nm存在一个吸收峰，R9:1535.53981nm。如图6中72和73为第一和第二直接光强探测单元信号，74和75是同步采集的第一和第二绝对波长监测单元信号(这里是氰化氢气室的透射信号)。对氰化氢气室透射信号而言，在两段波长重合的部分为1535-1536nm，区域，第一绝对波长监测单元信号74在采样点70处存在吸收峰，第二绝对波长监测单元信号75在采样点71处存在吸收峰，这两个吸收峰都为R9:1535.53981nm。因此可以用该位置作为拼接波长位置。对于第一直接光强探测单元信号72舍弃采样点70后边的数据，对于第二直接光强探测单元信号73舍弃采样点71前边的数据。新得到的第一和第二直接光强探测单元信号拼接得到新的直接光强探测单元信号77记为拼接后直接光强探测单元信号。同理通过绝对波长监测单元信号可以得到拼接后相对波长监测单元信号。
所述的波长监测单元包括输出特征信号的气室或中心波长已知的光纤光栅，或可以直接得到波长大小的光谱仪或波长计，或光纤干涉仪或FP标准器或光学谐振腔，或以上几种的组合。
辅助干涉仪单元21可以是为典型的马赫增德干涉仪或者迈克尔孙干涉仪。图6为调谐光信号经过光纤干涉仪的信号，正弦信号周期和干涉仪两个臂的光程差有关。正弦信号周期决定了该干涉仪的自由光谱范围，也就是每个正弦代表的光频间距。将该信号希尔伯特展开后可以得到光信号相位变化情况，因此该输出信号亦可作为波长追迹的信号(Ahn,T.and D.Y.Kim,Analysis of nonlinear frequency sweep in high-speedtunable laser sources using a self-homodyne measurement and Hilberttransformation.2007.46(13):p.2394.)。
如图3为半导体光放大器及其闭环功率稳定单元，其构成包括：半导体光放大器91，可以通过改变其上的驱动电流来对注入至半导体光放大器的激光功率施加不同增益；光电探测器93，用于将光信号转化为与光功率成比例的电学信号；比较器94，用于将测得的光强信号与设定值比较得到误差信号；比较器94输出的误差信号经过滤波及放大后的电压信号输入到电流驱动器97成为电流驱动信号。第十耦合器分出10％的光进入到闭环控制环路，剩下的90％作为输出光注入到后边的测量光路作为测量装置信号光。可以看出闭环功率控制单元40将输入的光稳定在设定的某一功率上。半导体光放大器的相关介绍可参见[[1]https://www.thorlabs.de/newgrouppage9.cfm？objectgroup_id＝3901
[2]https://courses.cit.cornell.edu/ece533/Lectures/handout9.pdf]
由于电流调谐带来的激光器波长调谐具有非线性，也就是输出光频随时间不是线性增长，此时若用固定采样率对各个输出信号进行采样，则采样点并非等光频间隔，这一效应恶化对传感或测量系统的空间分辨力。上述所述的利用辅助干涉仪信号对测量态主路干涉仪信号及测量态波长监测信号进行非线性校正方法，具有若干种不同的实现方式：可以为用采集单元19将该辅助干涉仪信号与其他路信号同步用固定采样率采集，然后在数据处理单元中进行对主路干涉仪信号或波长监测信号进行非线性校正。方法为将辅助干涉仪信号进行希尔伯特展开，进行相位解卷，然后对相位进行等分，如按照π弧度等分，得到对应的采样点，然后利用这些采样点对主路干涉仪信号与波长监测信号进行重采样，重采样后的主路干涉仪信号与波长监测信号为校正了非线性的信号。此外，还有非均匀傅里叶变换，去斜滤波器，PNC相位补偿等用后期软件处理方式校正非线性的实现形式。此外，还可以将辅助干涉仪输出的正弦信号作为采集单元19的时钟，用该时钟作为主路干涉仪信号和波长监测信号的采集时钟对这两路进行采集。此外，还可以采用光电锁相环技术进行非线性校正。考虑到存在的现有技术，这一部分不再进行赘述。相关文献可见(1.丁振扬,几种改进OFDR性能方法的提出及验证,2013,天津大学.2.Fan,X.,Y.Koshikiya and F.Ito,Phase-noise-compensated optical frequency domain reflectometry with measurementrange beyond laser coherence length realized using concatenative referencemethod.Optics letters,2007.32(22):p.3227 3.Swept-wavelength InterferometricInterrogation of Fiber Rayleigh Scatter for Distributed Sensing Applications4.Song,J.,et al.,Long-Range High Spatial Resolution Distributed Temperatureand Strain Sensing Based on Optical Frequency-Domain Reflectometry.IEEEPhotonics Journal,2014.6(3):p.1-8.5.张浩,电流调谐半导体激光器绝对距离测量技术研究,2016,哈尔滨工业大学.第15页6.Satyan N,Vasilyev A,Rakuljic G,et al.Precisecontrol of broadband frequency chirps using optoelectronic feedback.Opticsexpress,2009,17(18):15991-15999.)
基于光频域反射原理的分布式物理量测量或传感，其原理是利用传感或测量光纤上不同位置对应的拍频量大小不一样来进行微米量级的空间定位的。理论上，光纤上两点空间分辨力为Δz＝c/2nΔF,其中c为真空中光速，n为光纤中折射率，ΔF为可调谐激光器扫过的光频范围。所以为了降低光纤上两点空间分辨力的数值，提高空间分辨力应该扩大调谐范围。同时基于光频域反射原理的分布式物理量测量或传感通过光谱的偏移量解算物理量的变化量，因此调谐范围越大，分布式物理量的量程也越大。
鉴于基于光频域反射原理的分布式物理量测量或传感是相对测量，所以首先需要确定参考态信号，这一参考态在下文中是在外界第一个状态下采集得到的并将其存于计算机内存中，而下面所述的外界第二个状态表示测量态，相对于第一个参考态，传感光纤可能发生了扰动。高精度光频域反射过程如下：
第1步、外界第一个状态下得到参考态第一波长监测单元信号，参考态第二波长监测单元信号，参考态第一辅助干涉仪信号，参考态第二辅助干涉仪信号，参考态第一主路干涉仪信号，参考态第二主路干涉仪信号；
第2步、外界第二个状态下得到测量态第一波长监测单元信号，测量态第二波长监测单元信号，测量态第一辅助干涉仪信号，测量态第二辅助干涉仪信号，测量态第一主路干涉仪信号，测量态第二主路干涉仪信号；
第3步、确定拼接位置，得到拼接后参考态主路干涉仪信号，拼接后参考态辅助干涉仪信号，拼接后测量态主路干涉仪信号，和拼接后测量态主路干涉仪信号；
第4步、分别利用拼接后参考态辅助干涉仪信号和拼接后测量态辅助干涉仪信号对拼接后参考态主路干涉仪信号和拼接后测量态主路干涉仪信号进行非线性校正，得到最终参考态主路干涉仪信号和最终测量态主路干涉仪信号；
第5步、分布式物理量解算：对所述最终参考态干涉信号和所述最终测量态干涉信号分别进行快速傅里叶变换得到所述最终参考态干涉信号和所述最终测量态干涉信号的距离域信号，分别对所述最终参考态干涉信号和所述最终测量态干涉信号的距离域信号利用移动窗在距离域上的相同位置选取空间传感单元，将所述移动窗选定的空间传感单元信号进行傅里叶逆变换，得到所述移动窗对应的空间传感单元所对应的参考态瑞利散射光谱信号和测量态瑞利散射光谱信号；对两者进行互相关运算，求取互相关运算结果的最大值所在位置，该最大值所在位置对应所述的该位置的空间传感单元上的被测物理量变化；通过在距离域信号上滑动所述移动窗来选取距离域上不同位置的所述空间传感单元进而得到光纤上不同位置的物理量变化。
对未进行非线性校正的情形，则基于所述拼接后参考态干涉信号和所述拼接后测量态干涉信号分别进行快速傅里叶变换得到所述拼接后参考态干涉信号和所述拼接后测量态干涉信号的距离域信号，并进行后续解算。
传感光纤6可以为普通单模光纤，或者为刻写了中心波长相等的弱反射光纤光栅阵列的光纤，例如(Use of 3000Bragg grating strain sensors distributed on foureight-meter optical fibers during static load tests of a compositestructure.)公开的传感光纤，或者是瑞利散射增强的光纤，例如(Loranger,S.,et al.,Rayleigh scatter based order of magnitude increase in distributed temperatureand strain sensing by simple UV exposure of optical fibre.Scientific Reports,2015.5:p.11177.)公开的光纤等。
若传感光纤为具有等中心波长的弱反射光纤光栅阵列构成，则第5步的分布式物理量解算过程为：包括对所述最终参考态干涉信号和所述最终测量态干涉信号分别进行快速傅里叶变换得到所述最终参考态干涉信号和所述最终测量态干涉信号的距离域信号，利用窗函数选取所述距离域信号中各个光纤光栅对应的部分，分别利用傅里叶逆变换将所述选取的部分转换到光频域，分别得到各个光纤光栅的参考态信号和测量态信号下的光栅光谱信号；对所述参考态信号和测量态信号的光栅光谱信号进行包络探测并找到峰值所在位置，所述各个光纤光栅的参考态信号和测量态信号的光栅光谱信号峰值位置的差值代表该位置光栅上被测物理量大小。(李政颖,孙文丰与王洪海,基于光频域反射技术的超弱反射光纤光栅传感技术研究.光学学报,2015(08):第64-71页.)；若传感光纤为具有瑞利散射的光纤，则由于其没有光频域为随机分布的信号，没有单一的峰值，则其光谱偏移量可以用互相关求取，互相关的峰值位置对应光谱的偏移量(可参考Cui J,Zhao S,Yang D,etal.Investigation of the interpolation method to improve the distributedstrain measurement accuracy in optical frequency domain reflectometry systems[J].Applied optics,2018,57(6):1424-1431.)。此外，可以使用其他已知的成熟的方法对存在峰值的光谱信号进行偏移量解算，如最大值法，能量重心法等,相关文献见(Tosi,D.,Review and Analysis of Peak Tracking Techniques for Fiber Bragg GratingSensors.Sensors,2017.17(10):p.2368.)。
对未进行非线性校正的情形，则基于所述拼接后参考态干涉信号和所述拼接后测量态干涉信号分别进行快速傅里叶变换得到所述拼接后参考态干涉信号和所述拼接后测量态干涉信号的距离域信号，并进行后续解算。
分布式物理量解算中无论是用互相关还是峰值检测，直接得到的结果是光谱的偏移量，而光谱的偏移量是分布式物理量的响应函数。图2仅仅用拉伸39代表了被测的物理量，也就是分布式应变。但是分布式物理量可以为应变，或温度，或可以导致光纤发生应变或温度变化的其他物理量。这些物理量和光谱的偏移量之间一般相差一个系数(灵敏度)或者为了更加精确，符合多项式函数的关系，或者可以通过标定实验得到系数值(可参考CuiJ,Zhao S,Yang D,et al.Investigation of the interpolation method to improvethe distributed strain measurement accuracy in optical frequency domainreflectometry systems[J].Applied optics,2018,57(6):1424-1431.)。
应当理解，在本发明的一些实施例中，非线性校正并非必须的，本发明的测量方法可以以在不进行非线性校正的情形下实现。在本发明的一些实施例中的参考态信息，例如参考态的光谱可以预先存储作为基准信号，或者可以如上述实施例所描述的那样通过实时测量获得。
应当理解，在本发明的一些实施例中，非线性校正并非必须的，本发明的测量方法可以以在不进行非线性校正的情形下实现，即无需提供辅助干涉仪单元以及其所提供的辅助干涉光信号。在本发明的一些实施例中的参考态信息，例如参考态的包括拼接后参考态干涉光信号的光谱可以预先存储作为基准信号从而在物理量解算时直接调用，或者可以如上述实施例所描述的那样通过实时测量获得。
在本申请中，如无相反定义，度指摄氏度。
在本申请中，存储单元可以存储数据处理单元得到信号拼接位置，存储分布式物理量解算结果或直接存储原始采集信号以便后期离线处理。
本申请中，所述的闭环功率控制单元可以包括半导体光放大器，配置为通过改变其上的驱动电流来对注入至半导体光放大器的激光功率施加不同增益的功能；光电探测器，配置为将光信号转化为与光功率成比例的电学信号；比较器，配置为将测得的光强信号与设定值比较得到误差信号；以及执行器，配置为将该误差信号转换为电流驱动信号。
在本申请中，主路干涉光信号指由主路干涉仪单元或其他具有相同或实质相同功能的单元提供的信号；辅助干涉光信号指由辅路干涉仪单元或其他具有相同或实质相同功能的单元提供的信号；激光输出波长监测信号或简称波长监测信号是指由激光器波长监测单元或其他具有相同或实质相同功能单元提供的信号。
以上所述，仅为本发明较佳的具体实施方式，这些具体实施方式都是基于本发明整体构思下的不同实现方式，而且本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。本发明还可通过以下公开的示例实现：
1.一种高精度光频域反射方法，用于通过耦合于待测主体的传感光纤对所述待测主体的物理量变化进行测量，其特征在于该方法包括以下步骤：
对至少两个垂直腔面反射激光器中的每一个调整至不同的工作温度以使得所述各垂直腔面反射激光器的激光输出具有不同的输出波长；
对每个所述垂直腔面反射激光器进行电流调谐以提供调谐激光输出，其中所述调谐激光输出包括所述每个垂直腔面反射激光器提供的调谐激光输出波段，其中波长相邻的所述调谐激光输出波段有重叠区域；
在包括所述物理量变化的测量态下，同步地采集包括所述调谐激光输出的绝对波长的测量态波长监测信号；以及响应所述调谐激光输出于所述传感光纤对调谐激光输出的干涉的测量态主路干涉信号；
根据所述测量态波长监测信号确定各个垂直腔面反射激光器输出的测量态信号拼接位置，根据该位置对所述测量态主路干涉信号拼接得到拼接后测量态主路干涉信号；以及
基于所述拼接后测量态干涉光信号和在无所述物理量变化的参考态获得的拼接后参考态干涉光信号解算所述物理量变化。
2.根据示例1的方法，其特征在于：
所述提供调谐激光输出包括：
配置两个或两个以上垂直腔面反射激光器及其外围驱动单元，使得每个激光器的驱动电流可调，温度可控；
记录每个激光器在某一温度下，从起始驱动电流到终止驱动电流的输出激光的波长变化量；
固定驱动电流为所述起始驱动电流到所述终止驱动电流之间的某一电流值，将各个激光器配置为不同的温度，使得此时各个激光器的输出波长从小到大依次排列且依次存在小于所述波长变化量的波长间隔；各个温度下的各个激光器按照输出波长从小到大依次排列，其顺序被固定并保持各个激光器维持上述波长状态下的温度不变；
控制所述排列顺序的各个垂直腔面反射激光器依次改变其驱动电流得到波长调谐激光输出；其中所述驱动电流配置为所述起始驱动电流到所述终止驱动电流。
3.根据示例2的方法，其特征在于：对所述垂直腔面反射激光器的激光输出进行闭环功率控制来对垂直腔面反射激光器的输出激光的光功率进行实时调节，以输出功率稳定激光输出。
4.根据上述示例中任意一项所述的方法，其特征在于：还包括在所述测量态下同步地采集包括所述调谐激光输出的绝对波长的测量态波长监测信号以及响应所述调谐激光输出于所述传感光纤对调谐激光输出的干涉的测量态主路干涉信号的步骤中同步地采集响应所述调谐激光输出的相对干涉的测量态辅助干涉信号；根据所述测量态波长监测信号确定的各个垂直腔面反射激光器输出的测量态信号拼接位置对所述测量态辅助干涉信号进行信号拼接得到拼接后测量态辅助干涉信号；利用拼接后测量态辅助干涉仪信号对测量态主路干涉仪信号进行非线性校正，得到最终测量态主路干涉仪信号；并对所述最终测量态主路干涉仪信号进行快速傅里叶变换，将光频域信息转换到传感光纤中各个位置的距离域信息，利用移动窗扫描传感光纤中的各个位置的信号，利用傅里叶逆变换转换到光频域，对应各个位置测量态信号的瑞利散射光谱信号；求取所述测量态信号相对参考态信号的光谱偏移量，进而得到所述物理量变化。
5.根据上述示例中任意一项所述的方法，其特征在于：所述同步地采集包括所述调谐激光输出的绝对波长的测量态波长监测信号以及响应所述调谐激光输出于所述传感光纤对调谐激光输出的干涉的测量态主路干涉信号包括采集响应所述调谐激光输出的相对干涉的测量态辅助干涉信号，以及以所述相对干涉的测量态辅助干涉信号作为外时钟同步地采集所述测量态波长监测信号以及所述测量态主路干涉信号。
6.根据上述示例中任意一项所述的方法，其特征在于：
所述参考态信号预先保存或通过以下方式获取：
在不包括所述物理量变化的参考态下，同步地采集包括所述调谐激光输出的绝对波长的参考态波长监测信号；响应所述调谐激光输出的相对干涉的参考态辅助干涉信号；
以及响应所述调谐激光输出于所述传感光纤对调谐激光输出的干涉的参考态主路干涉信号；
根据所述参考态波长监测信号确定各个垂直腔面反射激光器输出的参考态信号拼接位置，根据该位置分别对所述参考态主路干涉信号和参考态辅助干涉信号进行信号拼接得到拼接后参考态辅助干涉信号和拼接后参考态主路干涉信号；
利用拼接后参考态辅助干涉信号对参考态主路干涉信号进行非线性校正，得到最终参考态主路干涉信号；以及
对最终参考态主路干涉信号进行快速傅里叶变换，将光频域信息转换到传感光纤中各个位置的距离域信息，利用移动窗扫描传感光纤中的各个位置的信号，利用傅里叶逆变换转换到光频域，对应各个位置参考态信号的瑞利散射光谱信号；求取所述参考态信号相对参考态信号的光谱偏移量，进而得到所述物理量变化。
7.根据上述示例中任意一项所述的方法，其特征在于：所述垂直腔面反射激光器的输出波长随电流和温度的增加而增加。
8.根据上述示例中任意一项所述的方法，其特征在于：所述的非线性校正过程包括：分别用参考态辅助干涉信号以及所述测量态辅助干涉信号作为所述同步采集的外时钟，触发主路干涉仪信号的同步采集；或者，用固定采样率的内时钟同步采集辅助干涉信号和主路干涉信号，根据辅助干涉信号估算激光器输出信号的相位，并由此对主路干涉信号进行非线性补偿，如重采样，使得主路干涉信号具有等光频间隔；或者采用光电锁相环实现非线性校正。
9.根据上述示例中任意一项所述的方法，其特征在于：所述的分布式物理量为应变，或温度，或可以导致光纤发生应变或温度变化的其他物理量。
10.根据以上示例中任意一项所述的方法，其特征在于：所述解算包括对所述拼接后参考态干涉光信号和所述拼接后测量态干涉光信号分别进行快速傅里叶变换得到所述拼接后参考态干涉光信号和所述拼接后测量态干涉光信号的距离域信号，分别对所述拼接后参考态干涉光信号和所述拼接后测量态干涉光信号的距离域信号利用移动窗在距离域上的相同位置选取空间传感单元，将所述移动窗选定的空间传感单元信号进行傅里叶逆变换，得到所述移动窗对应的空间传感单元所对应的参考态瑞利散射光谱信号和测量态瑞利散射光谱信号；对两者进行互相关运算，求取互相关运算结果的最大值所在位置，该最大值所在位置对应所述的该位置的空间传感单元上的被测物理量变化；
通过在距离域信号上滑动所述移动窗来选取距离域上不同位置的所述空间传感单元进而得到光纤上不同位置的物理量变化。
11.根据以上示例中任意一项所述的方法，其特征在于：对于具有等中心波长的弱反射光纤光栅阵列作为所述光纤传感器的情形；所述解算包括：包括对所述拼接后参考态干涉光信号和所述拼接后测量态干涉光信号分别进行快速傅里叶变换得到所述拼接后参考态干涉光信号和所述拼接后测量态干涉光信号的距离域信号，利用窗函数选取所述距离域信号中各个光纤光栅对应的部分，分别利用傅里叶逆变换将所述选取的部分转换到光频域，分别得到各个光纤光栅的参考态信号和测量态信号下的光栅光谱信号；
对所述参考态信号和测量态信号的光栅光谱信号进行包络探测并找到峰值所在位置，所述各个光纤光栅的参考态信号和测量态信号的光栅光谱信号峰值位置的差值代表该位置光栅上被测物理量大小。
12.一种高精度光频域反射装置，用于通过耦合于待测主体的传感光纤对所述待测主体的物理量变化进行测量，其特征在于该装置包括：
至少两个垂直腔面反射激光器，所述每一个垂直腔面反射激光器被配置为作以下调谐：
调整至不同的工作温度以使得所述各垂直腔面反射激光器的激光输出具有不同的输出波长；以及对每个所述垂直腔面反射激光器进行电流调谐以提供调谐激光输出，其中所述调谐激光输出包括所述每个垂直腔面反射激光器提供的调谐激光输出波段，其中波长相邻的所述调谐激光输出波段有重叠区域；
电流驱动单元，用于为所述垂直腔面反射激光器提供驱动电流信号，进而对垂直腔面反射激光器进行波长调谐；
温度控制单元，用于为所述垂直腔面反射激光器中的垂直腔面反射激光器提供的温度控制信号；
波长监测单元，用于对输出激光的波长进行测量得到激光输出波长值，以及用于根据波长监测信号中的光谱重叠部分来确定相邻波段的波长拼接位置；
主路干涉仪单元，用于测量传感光纤上的散射或者反射信号；
采集单元，用于采集激光器波长监测单元输出的信号或者其读数，用于采集辅助干涉仪信号，用于采集主路干涉仪信号；
控制单元，用于上述光源调谐电流驱动单元、温度控制单元、采集单元的时序控制；
数据处理单元，配置为对采集的数据进行运算，包括确定各路信号的拼接位置；对主路干涉仪信号进行截取和拼接得到拼接后主路干涉仪信号；对所述主路干涉仪信号、波长监测单元输出的信号进行非线性校正；以及对分布式物理量进行解算。
13.根据示例12所述的装置，其特征在于：还包括辅助干涉仪单元，用于校正主路干涉仪单元输出信号的光源波长调谐存在的非线性；
14.根据示例12所述的装置，其特征在于：所述辅助干涉仪单元包括具有马赫曾德结构或者迈克尔孙结构的光纤干涉仪。
15.根据上述示例中任意一项所述的装置，其特征在于：所述主路干涉仪单元包括具有马赫曾德结构或者迈克尔孙结构的光纤干涉仪。
16.根据上述示例中任意一项所述的装置，其特征在于：还包括闭环功率控制单元，用于调节垂直腔面反射激光器输出激光的光功率，以提供具有稳定光功率的激光输出。
17.根据示例16所述的装置，其特征在于：其特征在于所述的闭环功率控制单元包括：
半导体光放大器，可以通过改变其上的驱动电流来对注入至半导体光放大器的激光功率施加不同增益的功能；
光电探测器，用于将光信号转化为与光功率成比例的电学信号；
比较器，用于将测得的光强信号与设定值比较得到误差信号；
执行器，用于将该误差信号转换为电流驱动信号，该电流驱动信号用于驱动半导体光放大器，进而实时改变半导体光放大器的增益；
18.根据上述示例中任意一项所述的装置，其特征在于：还包括存储单元用于存储数据处理单元得到信号拼接位置，存储分布式物理量解算结果或直接存储原始采集信号以便后期离线处理。
19.根据上述示例中任意一项所述的装置，其特征在于：所述传感光纤位于主路干涉仪单元的测量臂上，可以为普通单模光纤，或者为刻写了中心波长相等的弱反射光纤光栅阵列的光纤，或者是瑞利散射增强的光纤。
20.根据上述示例中任意一项所述的装置，其特征在于：所述的波长监测单元包括输出特征信号的气室或中心波长已知的光纤光栅，或可以直接得到波长大小的光谱仪或波长计，或光纤干涉仪或FP标准器或光学谐振腔，或以上几种的组合。
21.一种高精度光频域反射系统，用于对所述待测主体的物理量变化进行测量，其特征在于该系统包括：
耦合于所述待测主体的传感光纤，
至少两个垂直腔面反射激光器，所述每一个垂直腔面反射激光器被配置为作以下调谐：
调整至不同的工作温度以使得所述各垂直腔面反射激光器的激光输出具有不同的输出波长；以及对每个所述垂直腔面反射激光器进行电流调谐以提供调谐激光输出，其中所述调谐激光输出包括所述每个垂直腔面反射激光器提供的调谐激光输出波段，其中波长相邻的所述调谐激光输出波段有重叠区域；
电流驱动单元，用于为所述垂直腔面反射激光器提供驱动电流信号，进而对垂直腔面反射激光器进行波长调谐；
温度控制单元，用于为所述垂直腔面反射激光器中的垂直腔面反射激光器提供的温度控制信号；
波长监测单元，用于对输出激光的波长进行测量得到激光输出波长值，以及用于根据波长监测信号中的光谱重叠部分来确定相邻波段的信号拼接位置；
主路干涉仪单元，用于测量传感光纤上的散射或者反射信号；
采集单元，用于采集激光器波长监测单元输出的信号或者其读数，用于采集辅助干涉仪信号，用于采集主路干涉仪信号；
控制单元，用于上述光源调谐电流驱动单元、温度控制单元、采集单元的时序控制；
数据处理单元，配置为对采集的数据进行运算，包括确定各路信号的拼接位置；对主路干涉仪信号进行截取和拼接得到拼接后主路干涉仪信号；对所述主路干涉仪信号、波长监测单元输出的信号进行非线性校正；以及对分布式物理量进行解算。
22.根据示例21所述的系统，其特征在于：还包括辅助干涉仪单元，用于校正主路干涉仪单元输出信号的光源波长调谐存在的非线性；
23.根据示例21所述的系统，其特征在于：所述辅助干涉仪单元包括具有马赫曾德结构或者迈克尔孙结构的光纤干涉仪。
24.根据上述示例中任意一项所述的系统，其特征在于：所述主路干涉仪单元包括具有马赫曾德结构或者迈克尔孙结构的光纤干涉仪。
25.根据上述示例中任意一项所述的系统，其特征在于：还包括闭环功率控制单元，用于调节垂直腔面反射激光器输出激光的光功率，以提供具有稳定光功率的激光输出。
26.根据示例25所述的系统，其特征在于：其特征在于所述的闭环功率控制单元包括：
27.半导体光放大器，可以通过改变其上的驱动电流来对注入至半导体光放大器的激光功率施加不同增益的功能；
光电探测器，用于将光信号转化为与光功率成比例的电学信号；
比较器，用于将测得的光强信号与设定值比较得到误差信号；
执行器，用于将该误差信号转换为电流驱动信号，该电流驱动信号用于驱动半导体光放大器，进而实时改变半导体光放大器的增益；
28.根据上述示例中任意一项所述的系统，其特征在于：还包括存储单元用于存储数据处理单元得到信号拼接位置，存储分布式物理量解算结果或直接存储原始采集信号以便后期离线处理。
29.根据上述示例中任意一项所述的系统，其特征在于：所述传感光纤位于主路干涉仪单元的测量臂上，可以为普通单模光纤，或者为刻写了中心波长相等的弱反射光纤光栅阵列的光纤，或者是瑞利散射增强的光纤。
30.根据上述示例中任意一项所述的系统，其特征在于：所述的波长监测单元包括输出特征信号的气室或中心波长已知的光纤光栅，或可以直接得到波长大小的光谱仪或波长计，或光纤干涉仪或FP标准器或光学谐振腔，或以上几种的组合。

Claims (29)

1.一种高精度光频域反射方法，用于通过耦合于待测主体的传感光纤对所述待测主体的物理量变化进行测量，其特征在于该方法包括以下步骤：
对至少两个垂直腔面反射激光器中的每一个调整至不同的工作温度以使得各所述垂直腔面反射激光器的激光输出具有不同的输出波长；
对每个所述垂直腔面反射激光器进行电流调谐以提供调谐激光输出，其中所述调谐激光输出包括每个所述垂直腔面反射激光器提供的调谐激光输出波段，其中波长相邻的所述调谐激光输出波段有重叠区域；
在包括所述物理量变化的测量态下，同步地采集包括所述调谐激光输出的绝对波长的测量态波长监测信号；以及响应所述调谐激光输出于所述传感光纤对调谐激光输出的干涉的测量态主路干涉信号；
根据所述测量态波长监测信号确定各个垂直腔面反射激光器输出的测量态主路干涉信号拼接位置，根据该位置对所述测量态主路干涉信号拼接得到拼接后测量态主路干涉信号；以及
基于所述拼接后测量态主路干涉信号和在无所述物理量变化的参考态获得的拼接后参考态主路干涉信号解算所述物理量变化。
2.根据权利要求1所述的方法，其特征在于：
所述提供调谐激光输出包括：
配置两个或两个以上垂直腔面反射激光器及其外围驱动单元，使得每个激光器的驱动电流可调，温度可控；
记录每个激光器在某一温度下，从起始驱动电流到终止驱动电流的输出激光的波长变化量；
固定驱动电流为所述起始驱动电流到所述终止驱动电流之间的某一电流值，将各个激光器配置为不同的温度，使得此时各个激光器的输出波长从小到大依次排列且依次存在小于所述波长变化量的波长间隔；各个温度下的各个激光器按照输出波长从小到大依次排列，其顺序被固定并保持各个激光器维持上述波长状态下的温度不变；
控制所述顺序的各个垂直腔面反射激光器依次改变其驱动电流得到波长调谐激光输出；其中所述驱动电流配置为所述起始驱动电流到所述终止驱动电流。
3.根据权利要求2所述的方法，其特征在于：对所述垂直腔面反射激光器的激光输出进行闭环功率控制来对垂直腔面反射激光器的输出激光的光功率进行实时调节，以输出功率稳定激光输出。
4.根据权利要求1-3中任意一项所述的方法，其特征在于：还包括在所述测量态下同步地采集包括所述调谐激光输出的绝对波长的测量态波长监测信号以及响应所述调谐激光输出于所述传感光纤对调谐激光输出的干涉的测量态主路干涉信号的步骤中同步地采集响应所述调谐激光输出的相对干涉的测量态辅助干涉信号；根据所述测量态波长监测信号确定的各个垂直腔面反射激光器输出的测量态信号拼接位置对所述测量态辅助干涉信号进行信号拼接得到拼接后测量态辅助干涉信号；利用拼接后测量态辅助干涉信号对测量态主路干涉信号进行非线性校正，得到最终测量态主路干涉信号；并对所述最终测量态主路干涉信号进行快速傅里叶变换，将光频域信息转换到传感光纤中各个位置的距离域信息，利用移动窗扫描传感光纤中的各个位置的信号，利用傅里叶逆变换转换到光频域，对应各个位置测量态信号的瑞利散射光谱信号；求取所述测量态信号相对参考态信号的光谱偏移量，进而得到所述物理量变化。
5.根据权利要求4所述的方法，其特征在于：所述同步地采集包括所述调谐激光输出的绝对波长的测量态波长监测信号以及响应所述调谐激光输出于所述传感光纤对调谐激光输出的干涉的测量态主路干涉信号包括采集响应所述调谐激光输出的相对干涉的测量态辅助干涉信号，以及以所述相对干涉的测量态辅助干涉信号作为外时钟同步地采集所述测量态波长监测信号以及所述测量态主路干涉信号。
6.根据权利要求5所述的方法，其特征在于：
所述参考态信号预先保存或通过以下方式获取：
在不包括所述物理量变化的参考态下，同步地采集包括所述调谐激光输出的绝对波长的参考态波长监测信号；响应所述调谐激光输出的相对干涉的参考态辅助干涉信号；以及响应所述调谐激光输出于所述传感光纤对调谐激光输出的干涉的参考态主路干涉信号；
根据所述参考态波长监测信号确定各个垂直腔面反射激光器输出的参考态信号拼接位置，根据该位置分别对所述参考态主路干涉信号和参考态辅助干涉信号进行信号拼接得到拼接后参考态辅助干涉信号和拼接后参考态主路干涉信号；
利用拼接后参考态辅助干涉信号对参考态主路干涉信号进行非线性校正，得到最终参考态主路干涉信号；以及
对最终参考态主路干涉信号进行快速傅里叶变换，将光频域信息转换到传感光纤中各个位置的距离域信息，利用移动窗扫描传感光纤中的各个位置的信号，利用傅里叶逆变换转换到光频域，对应各个位置参考态信号的瑞利散射光谱信号；求取所述参考态信号相对参考态信号的光谱偏移量，进而得到所述物理量变化。
7.根据权利要求6所述的方法，其特征在于：所述垂直腔面反射激光器的输出波长随电流和温度的增加而增加。
8.根据权利要求6所述的方法，其特征在于：所述非线性校正过程包括：分别用参考态辅助干涉信号以及所述测量态辅助干涉信号作为所述同步地采集的外时钟，触发主路干涉信号的同步地采集；或者，用固定采样率的内时钟同步地采集辅助干涉信号和主路干涉信号，根据辅助干涉信号估算激光器输出信号的相位，并由此对主路干涉信号进行非线性补偿，如重采样，使得主路干涉信号具有等光频间隔；或者采用光电锁相环实现非线性校正。
9.根据权利要求1所述的方法，其特征在于：所述的物理量为应变，或温度，或导致光纤发生应变或温度变化的其他物理量。
10.根据权利要求1所述的方法，其特征在于：所述解算包括对拼接后参考态干涉光信号和拼接后测量态干涉光信号分别进行快速傅里叶变换得到所述拼接后参考态干涉光信号和所述拼接后测量态干涉光信号的距离域信号，分别对所述拼接后参考态干涉光信号和所述拼接后测量态干涉光信号的距离域信号利用移动窗在距离域上的相同位置选取空间传感单元，将所述移动窗选定的空间传感单元信号进行傅里叶逆变换，得到所述移动窗对应的空间传感单元所对应的参考态瑞利散射光谱信号和测量态瑞利散射光谱信号；对两者进行互相关运算，求取互相关运算结果的最大值所在位置，该最大值所在位置对应该位置的空间传感单元上的被测物理量变化；通过在距离域信号上滑动所述移动窗来选取距离域上不同位置的所述空间传感单元进而得到光纤上不同位置的物理量变化。
11.根据权利要求1所述的方法，其特征在于：对于具有等中心波长的弱反射光纤光栅阵列作为所述传感光纤的情形；所述解算包括：包括对拼接后参考态干涉光信号和拼接后测量态干涉光信号分别进行快速傅里叶变换得到所述拼接后参考态干涉光信号和所述拼接后测量态干涉光信号的距离域信号，利用窗函数选取所述距离域信号中各个光纤光栅对应的部分，分别利用傅里叶逆变换将选取的部分转换到光频域，分别得到各个光纤光栅的参考态信号和测量态信号下的光栅光谱信号；对所述参考态信号和测量态信号的光栅光谱信号进行包络探测并找到峰值所在位置，各个所述光纤光栅的参考态信号和测量态信号的光栅光谱信号峰值位置的差值代表该位置光栅上被测物理量大小。
12.一种高精度光频域反射装置，用于通过耦合于待测主体的传感光纤对所述待测主体的物理量变化进行测量，其特征在于该装置包括：
至少两个垂直腔面反射激光器，每一个所述垂直腔面反射激光器被配置为作以下调谐：调整至不同的工作温度以使得各所述垂直腔面反射激光器的激光输出具有不同的输出波长；以及对每个所述垂直腔面反射激光器进行电流调谐以提供调谐激光输出，其中所述调谐激光输出包括每个所述垂直腔面反射激光器提供的调谐激光输出波段，其中波长相邻的所述调谐激光输出波段有重叠区域；
电流驱动单元，用于为所述垂直腔面反射激光器提供驱动电流信号，进而对垂直腔面反射激光器进行波长调谐；
温度控制单元，用于为所述垂直腔面反射激光器中的垂直腔面反射激光器提供的温度控制信号；
波长监测单元，用于对输出激光的波长进行测量得到激光输出波长值，以及用于根据波长监测信号中的光谱重叠部分来确定相邻波段的波长拼接位置；
主路干涉仪单元，用于测量传感光纤上的散射或者反射信号；
采集单元，用于采集激光器波长监测单元输出的信号或者其读数，用于采集辅助干涉仪信号，用于采集主路干涉仪信号；
控制单元，用于光源调谐电流驱动单元、温度控制单元、采集单元的时序控制；
数据处理单元，配置为对采集的数据进行运算，包括确定各路信号的拼接位置；对主路干涉仪信号进行截取和拼接得到拼接后主路干涉仪信号；对所述主路干涉仪信号、波长监测单元输出的信号进行非线性校正；以及对分布式物理量进行解算。
13.根据权利要求12所述的装置，其特征在于：还包括辅助干涉仪单元，用于校正主路干涉仪单元输出信号的光源波长调谐存在的非线性。
14.根据权利要求13所述的装置，其特征在于：所述辅助干涉仪单元包括具有马赫曾德结构或者迈克尔孙结构的光纤干涉仪。
15.根据权利要求12-14中任意一项所述的装置，其特征在于：所述主路干涉仪单元包括具有马赫曾德结构或者迈克尔孙结构的光纤干涉仪。
16.根据权利要求12中所述的装置，其特征在于：还包括闭环功率控制单元，用于调节垂直腔面反射激光器输出激光的光功率，以提供具有稳定光功率的激光输出。
17.根据权利要求16所述的装置，其特征在于：其特征在于所述的闭环功率控制单元包括：半导体光放大器，通过改变其上的驱动电流来对注入至半导体光放大器的激光功率施加不同增益的功能；
光电探测器，用于将光信号转化为与光功率成比例的电学信号；
比较器，用于将测得的光强信号与设定值比较得到误差信号；
执行器，用于将该误差信号转换为电流驱动信号，该电流驱动信号用于驱动半导体光放大器，进而实时改变半导体光放大器的增益。
18.根据权利要求12所述的装置，其特征在于：还包括存储单元用于存储数据处理单元得到信号拼接位置，存储分布式物理量解算结果或直接存储原始采集信号以便后期离线处理。
19.根据权利要求12所述的装置，其特征在于：所述传感光纤位于主路干涉仪单元的测量臂上，为普通单模光纤，或者为刻写了中心波长相等的弱反射光纤光栅阵列的光纤，或者是瑞利散射增强的光纤。
20.根据权利要求12所述的装置，其特征在于：所述的波长监测单元包括输出特征信号的气室或中心波长已知的光纤光栅，或直接得到波长大小的光谱仪或波长计，或光纤干涉仪或FP标准器或光学谐振腔，或以上几种的组合。
21.一种高精度光频域反射系统，用于对待测主体的物理量变化进行测量，其特征在于该系统包括：
耦合于所述待测主体的传感光纤，
至少两个垂直腔面反射激光器，每一个所述垂直腔面反射激光器被配置为作以下调谐：调整至不同的工作温度以使得各所述垂直腔面反射激光器的激光输出具有不同的输出波长；以及对每个所述垂直腔面反射激光器进行电流调谐以提供调谐激光输出，其中所述调谐激光输出包括每个所述垂直腔面反射激光器提供的调谐激光输出波段，其中波长相邻的所述调谐激光输出波段有重叠区域；
电流驱动单元，用于为所述垂直腔面反射激光器提供驱动电流信号，进而对垂直腔面反射激光器进行波长调谐；
温度控制单元，用于为所述垂直腔面反射激光器中的垂直腔面反射激光器提供的温度控制信号；
波长监测单元，用于对输出激光的波长进行测量得到激光输出波长值，以及用于根据波长监测信号中的光谱重叠部分来确定相邻波段的信号拼接位置；
主路干涉仪单元，用于测量传感光纤上的散射或者反射信号；
采集单元，用于采集激光器波长监测单元输出的信号或者其读数，用于采集辅助干涉仪信号，用于采集主路干涉仪信号；
控制单元，用于光源调谐电流驱动单元、温度控制单元、采集单元的时序控制；
数据处理单元，配置为对采集的数据进行运算，包括确定各路信号的拼接位置；对主路干涉仪信号进行截取和拼接得到拼接后主路干涉仪信号；对所述主路干涉仪信号、波长监测单元输出的信号进行非线性校正；以及对分布式物理量进行解算。
22.根据权利要求21所述的系统，其特征在于：还包括辅助干涉仪单元，用于校正主路干涉仪单元输出信号的光源波长调谐存在的非线性。
23.根据权利要求22所述的系统，其特征在于：所述辅助干涉仪单元包括具有马赫曾德结构或者迈克尔孙结构的光纤干涉仪。
24.根据权利要求21-23中任意一项所述的系统，其特征在于：所述主路干涉仪单元包括具有马赫曾德结构或者迈克尔孙结构的光纤干涉仪。
25.根据权利要求24所述的系统，其特征在于：还包括闭环功率控制单元，用于调节垂直腔面反射激光器输出激光的光功率，以提供具有稳定光功率的激光输出。
26.根据权利要求25所述的系统，其特征在于：其特征在于所述的闭环功率控制单元包括：半导体光放大器，通过改变其上的驱动电流来对注入至半导体光放大器的激光功率施加不同增益的功能；
光电探测器，用于将光信号转化为与光功率成比例的电学信号；
比较器，用于将测得的光强信号与设定值比较得到误差信号；
执行器，用于将该误差信号转换为电流驱动信号，该电流驱动信号用于驱动半导体光放大器，进而实时改变半导体光放大器的增益。
27.根据权利要求21所述的系统，其特征在于：还包括存储单元用于存储数据处理单元得到信号拼接位置，存储分布式物理量解算结果或直接存储原始采集信号以便后期离线处理。
28.根据权利要求21所述的系统，其特征在于：所述传感光纤位于主路干涉仪单元的测量臂上，为普通单模光纤，或者为刻写了中心波长相等的弱反射光纤光栅阵列的光纤，或者是瑞利散射增强的光纤。
29.根据权利要求21所述的系统，其特征在于：所述的波长监测单元包括输出特征信号的气室或中心波长已知的光纤光栅，或直接得到波长大小的光谱仪或波长计，或光纤干涉仪或FP标准器或光学谐振腔，或以上几种的组合。
CN202010415101.0A 2020-05-15 2020-05-15 高精度光频域反射方法、装置及系统 Active CN113670350B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010415101.0A CN113670350B (zh) 2020-05-15 2020-05-15 高精度光频域反射方法、装置及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010415101.0A CN113670350B (zh) 2020-05-15 2020-05-15 高精度光频域反射方法、装置及系统

Publications (2)

Publication Number Publication Date
CN113670350A CN113670350A (zh) 2021-11-19
CN113670350B true CN113670350B (zh) 2022-06-21

Family

ID=78537836

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010415101.0A Active CN113670350B (zh) 2020-05-15 2020-05-15 高精度光频域反射方法、装置及系统

Country Status (1)

Country Link
CN (1) CN113670350B (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107990926A (zh) * 2017-11-23 2018-05-04 武汉理工大学 基于波长偏移的光纤光栅逐个锁定解调系统及方法
CN108562237A (zh) * 2018-01-04 2018-09-21 大连理工大学 一种采用hcn气室在光频域反射传感系统中进行光谱校准的装置和方法
CN109186644A (zh) * 2018-08-29 2019-01-11 北京航天控制仪器研究所 光谱拼接光频域反射型分布式光纤传感器及信号解调方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
IES20030281A2 (en) * 2003-04-14 2004-10-20 Intune Technologies Ltd Method and apparatus for continuous sweeping of a tunable laser
KR100488221B1 (ko) * 2003-09-08 2005-05-10 주식회사 파이버프로 광섬유 브래그 격자 센서 시스템

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107990926A (zh) * 2017-11-23 2018-05-04 武汉理工大学 基于波长偏移的光纤光栅逐个锁定解调系统及方法
CN108562237A (zh) * 2018-01-04 2018-09-21 大连理工大学 一种采用hcn气室在光频域反射传感系统中进行光谱校准的装置和方法
CN109186644A (zh) * 2018-08-29 2019-01-11 北京航天控制仪器研究所 光谱拼接光频域反射型分布式光纤传感器及信号解调方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
幅度取样光纤光栅拼接制作技术的分析与研究;张丽萍等;《半导体光电》;20110215;第32卷(第01期);64-68 *

Also Published As

Publication number Publication date
CN113670350A (zh) 2021-11-19

Similar Documents

Publication Publication Date Title
US11721948B2 (en) Swept light source and drive data generation method and optical deflector for swept light source
CN103828146B (zh) 从激光器动态自适应地生成波长连续的且规定的波长对时间的扫描的系统和方法
US7333680B2 (en) Fiber Bragg grating sensor system
US8392138B2 (en) System and method for correcting sampling errors associated with radiation source tuning rate fluctuations in swept-wavelength interferometry
EP2980537A1 (en) Multi-peak reference grating
CN112082581B (zh) 基于分布式反馈激光器阵列的准分布式物理量测量方法、装置及系统
Zhu et al. Fabry–Perot sensor using cascaded chirped fiber Bragg gratings with opposite chirp directions
CN112086857B (zh) 连续可调谐激光拼接位置的获取方法及装置以及对应的连续可调谐激光的输出方法及装置
CN113670350B (zh) 高精度光频域反射方法、装置及系统
CN112082585B (zh) 基于干涉光信号拼接的驱动电流调谐分布式测量方法、装置及系统
CN113670352B (zh) 低成本准分布式物理量测量方法，装置及系统
CN112082584B (zh) 基于激光器调谐控制的光纤分布式物理量测量方法、装置及系统
CN113670347B (zh) 一种高分辨力准分布式物理量测量方法、装置及系统
CN112082586B (zh) 基于分布式反馈激光器阵列的光纤光栅阵列传感方法、装置及系统
CN113670348B (zh) 高精度光纤分布式物理量测量方法、装置和系统
Liu et al. An integrated design of ultra-high-speed FBG interrogation system-based on FDML laser
CN113670349B (zh) 基于光频域反射技术的高速光纤分布式物理量测量方法、装置及系统
CN113670351B (zh) 基于vcsel复用的光纤光栅阵列测量方法、装置及系统
Chen et al. A low bandwidth DFB laser-based interrogator for terahertz-range fiber Bragg grating sensors
Zhao et al. Nonlinear correction in the frequency scanning interferometry system by a fiber resonator
CN109724633B (zh) 一种低反射率外腔光纤f-p干涉型传感器的解调系统
US6816264B1 (en) Systems and methods for amplified optical metrology
CN114061637B (zh) 一种基于光纤瑞利散射谱相关性的激光器探测信号拼接及其拓展方法
Li et al. Instant laser wavelength measurement based on dual Fabry–Pérot Etalons
Choi et al. Wavelength stabilization and intensity correction of fiber Fabry–Pérot tunable filter-based wavelength-swept laser

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