本发明涉及一种光器件时延测量方法，尤其涉及一种基于相推法的光器件时延测量方法及装置。
背景技术：
常用的光器件时延测量方法主要有脉冲法、频率扫描干涉法和相推法三种。脉冲法通过观测发射光脉冲与接收光脉冲的时间间隔计算出被测光器件的时延，由于光器件色散会对光脉冲进行展宽，恶化测量精度，因此脉冲法不适合对长光纤等色散量较大的光器件进行精确测量。此外，窄脉冲的频谱范围较宽，在测量通带范围较小的光器件时，无法全部通过，进而恶化测量精度(例如：超密集波分复用器、光控时延芯片等)。脉冲法存在着许多不可避免的误差，如仪器分辨力误差、光器件色散误差等。因此脉冲法的测量精度只是米量级，且随着光器件时延的增加，测量误差也随着增大。频率扫描干涉法需要使用连续扫频激光器，价格昂贵，而且受限于这种激光器的线宽跟扫频线性度，其测量范围较小，一般为公里(10微秒)量级，而且测量精度随着光器件时延的增大而明显减小。此外，频率扫描干涉法本质上是以大带宽换取高测量精度，所以在测量通带范围小的光器件时，测量精度会下降。相推法由于使用相位变化来推算光器件时延，精度较高，且可以规避大时延量恶化精度的问题。但是现有相推法在测量大时延的时候，需要精细的频率扫描，扫描频点数激增，测量时间大大加长，容易引入环境误差。
传统相推法使用光矢量分析仪测量光器件的相位响应，进而算出光器件时延响应。然而较高的时延测量精度需要极宽的扫频范围，对测量仪器要求较高，价格昂贵，且容易受到色散的影响。为克服这一问题，2019年李树鹏等人(s.p.li,x.c.wang,t.qing,s.f.liu,j.b.fu,m.xue,s.l.pan,"opticalfibertransferdelaymeasurementbasedonphase-derivedranging,"ieeephotonicstechnologyletters,vol.31,no.16,pp.1351-1354,aug.2019.)提出的基于相推法的高精度光纤时延测量系统，通过对窄线宽光源进行外调制，调制后的光信号通过环形器进入被测光纤，反射光经过光电转换之后通过鉴相器得到调制信号在待测光纤中经历的相位变化，再以固定的频率间隔在较小的一段频率范围内对调制信号进行线性扫频，得到一系列的相位变化，并由此推算出光纤时延。由于扫频范围小，不仅降低了对器件的要求，避免了色散的影响而且适用于测量通带范围小的光器件，可测光器件的种类范围较大。但是由于扫频的频率间隔决定了能测量的光器件时延，所以测量大时延的时候，需要较小的频率间隔，而频率范围不变，那么扫描的频点数就会激增，所需要的测量时间随之增长，不仅测量效率低，还会引入环境误差。
技术实现要素：
本发明所要解决的技术问题在于克服现有相推法光器件时延测量技术的不足，提供一种基于相推法的光器件时延测量方法，可大幅度减少扫描频点数，从而提高测量效率并减少环境误差对测量的影响。
本发明具体采用以下技术方案解决上述技术问题：
一种基于相推法的光器件时延测量方法，预先根据测量需求确定扫频频率范围并在其中选取多个扫频频点；在每一个扫频频点，用该频率的微波调制信号对光载波进行调制，并通过鉴相器测量出调制光信号经过待测光器件前后微波调制信号的相位变化；对所测得的一系列相位变化进行相位展开，并利用相位展开所得到的各扫频频点的展开相位计算出最大扫频频点的整周模糊度，最后根据所述最大扫频频点的整周模糊度计算出待测光器件的时延；所选取的最小扫频频点频率ωa、最大扫频频点频率ωb、扫频频点数量m以及各扫频频点的频率ωi具体如下：
其中，δθ为所述鉴相器的相位精度，δτ为时延测量目标精度，τ0为测量系统的时延，τmax为最大可测时延，λ为取值范围为(0,1]的修正系数。
优选地，所述相位展开的方法具体如下：
扫频频点ω2的展开相位φ(ω2)通过经典相位展开算法得到，其余扫频频点的展开相位通过下式得到：
其中，[…]为取整符号，表示对其中的数据进行“四舍五入”取整。
优选地，最大扫频频点的整周模糊度nb按照下式得到：
其中，[…]为取整符号，表示对其中的数据进行“四舍五入”取整。
优选地，所述修正系数λ的取值范围为[0.80,0.99]。
优选地，所述根据最大扫频频点的整周模糊度计算出待测光器件时延，具体根据以下公式：
其中，τd为待测光器件的时延，nb为最大扫频频点的整周模糊度。
根据相同的发明思路还可以得到以下技术方案：
一种基于相推法的光器件时延测量装置，包括：
频点确定单元，用于预先根据测量需求确定扫频频率范围并在其中选取多个扫频频点；相位测量单元，用于在每一个扫频频点，用该频率的微波调制信号对光载波进行调制，并通过鉴相器测量出调制光信号经过待测光器件前后微波调制信号的相位变化；
解算单元，用于对所测得的一系列相位变化进行相位展开，并利用相位展开所得到的各扫频频点的展开相位计算出最大扫频频点的整周模糊度，最后根据所述最大扫频频点的整周模糊度计算出待测光器件的时延；
频点确定单元所选取的最小扫频频点频率ωa、最大扫频频点频率ωb、扫频频点数量m以及各扫频频点的频率ωi具体如下：
其中，δθ为所述鉴相器的相位精度，δτ为时延测量目标精度，τ0为测量系统的时延，τmax为最大可测时延，λ为取值范围为(0,1]的修正系数。
优选地，所述相位展开的方法具体如下：
扫频频点ω2的展开相位φ(ω2)通过经典相位展开算法得到，其余扫频频点的展开相位通过下式得到：
其中，[…]为取整符号，表示对其中的数据进行“四舍五入”取整。
优选地，最大扫频频点的整周模糊度nb按照下式得到：
其中，[…]为取整符号，表示对其中的数据进行“四舍五入”取整。
优选地，所述修正系数λ的取值范围为[0.80,0.99]。
优选地，所述根据最大扫频频点的整周模糊度计算出待测光器件的时延，具体根据以下公式：
其中，τd为待测光器件的时延，nb为最大扫频频点的整周模糊度。
相比现有技术，本发明技术方案具有以下有益效果：
本发明对现有基于相推法的高精度光器件时延测量技术进行改进，采用非线性扫频，频率间隔按照指数级递增，在相同的扫描范围内，相比于现有的线性扫频方式，本方法需要扫描的频点数较少，扫描时间大大减少；此外，由于扫描频点大幅减少，对微波扫频源的要求也大幅降低，同时减少了环境误差的引入，提高了测量精度。
具体实施方式
针对现有基于相推法的高精度光器件时延测量技术由于扫描频点数过多所导致的测量效率低且易引入环境噪声的不足，本发明的解决思路是摈弃传统的线性扫描方式，转而采用非线性扫频，频率间隔按照指数级递增。这样，在相同的扫描范围内，需要扫描的频点数和扫描时间大大减少；此外，由于扫描频点大幅减少，对微波扫频源的要求也大幅降低，同时减少了环境误差的引入，提高了测量精度。
为便于公众理解，在对本发明技术方案进行详细说明之前，先对现有基于相推法的高精度光器件时延测量技术的基本原理进行介绍。
以(s.p.li,x.c.wang,t.qing,s.f.liu,j.b.fu,m.xue,s.l.pan,"opticalfibertransferdelaymeasurementbasedonphase-derivedranging,"ieeephotonicstechnologyletters,vol.31,no.16,pp.1351-1354,aug.2019.)中的测量系统为例，激光源发出一束光载波到mzm调制器，偏置点控制器控制mzm调制器的偏置点位于线性点，再将微波源输出的微波信号加载到mzm调制器的rf输入口，所产生的探测光信号均可表示为：
eo(t)＝a(1+mcos(ωet))exp(jωct)(1)
其中，a是光场幅度，ωe和ωc分别是微波信号和光载波的角频率，m是调幅系数。探测光经过光环形器到待测光器件传输后到达高灵敏度光电探测模块，此时的光场可表示为：
er(t)＝a(1+mcos(ωe(t-τ0-τd)))exp(jωc(t-τ0-τd))(2)
其中，τ0是探测光在测量系统中的总传输时间，τd是探测光在待测光器件中的总传输时间。返回的探测光经过光电转换后，通过滤波器提取出一倍频信号，其电场可表示为：
i(t)＝2ηma2cos(ωe(t-τ0-τd))(3)
其中，η为光电转换系数。从公式(3)可以看出，频率为ωe的微波信号的相位变化可表示为：
由于鉴相器测得的相位范围一般为-π到π，因此鉴相器得到的频率为ωe的微波信号的相位变化可表示为：
其中，n为整数,一般也叫做频率ωe的整周模糊度。
现有的相推法通过线性扫频来解算整周模糊度。首先确定扫频区间[ωa,ωb]，起始频率ωa和终止频率ωb可由如下公式确定：
其中，δθ为鉴相器的相位精度，δτ为所需要的光器件时延的测量精度。利用这种方法确定的扫频区间较小，所以不仅可以用来测量光纤也可以用于测量光芯片等通带范围较小的光器件。不过现有的相推法需要使用经典相位展开算法，所以它扫频的频率步进δω需要满足
其中，[…]为取整符号，此处为“进一法”取整，即：去掉数值的小数部分再加一，如果数值本身是整数则不变。由此，可以得到需要扫描的各个频率点：
依次扫描各个频点，可以得到相应的相位响应θ(ω1)，θ(ω2)…θ(ωm)。再通过经典相位展开算法得到展开的相位响应，可以用matlab的unwrap函数来实现。这种经典相位展开算法的描述如下：判断相邻两个频点测得的相位差的绝对值是否大于等于π，如果满足，就将后一个频率点的相位加上或者减去一个2π，使得相邻频点测得的相位差的绝对值小于π。展开的相位响应随频率呈线性变化，可利用最小二乘法线性拟合得到展开的相位响应随频率变化的线性函数φ(ωi)＝k·ωi+b，再将拟合得到的斜率k用于解算ωb的整周模糊度：
其中，[…]为取整运算符，此处为“四舍五入”取整，如果数值是负数，则先取绝对值，再“四舍五入”后加上负号。最后，结合公式(4)、(5)可得待测光器件的时延为：
在测量光纤这种传输介质时，其长度可以通过下式计算：
其中，c为真空中光速，n0为光纤折射率(可由光纤生产厂商提供的手册中查询得到)。如果是反射式测量，还需要除以二。
这种测量方法受限于经典相位展开算法，其扫频点数随着最大可测时延的增大呈线性增长，在测量大时延的时候，测量速度较慢。例如：测量一段公里级别长度的光纤，测量精度为0.1毫米的话，需要扫描成千上万个频点。
与该方法不同，本发明使用非线性频率扫描，后续频率与起始频率的频差呈指数级上升，具体的扫频频率确定方法如下。首先确定起始频率ωa和终止频率ωb：
其中，修正系数λ∈(0,1]，这是一个经验值，视系统稳定性而定，其优选的取值范围为[0.80,0.99]。接着确定扫频点数：
其中，[…]为取整符号，此处为“进一法”取整。最后确定起始频率ω1＝ωa与终止频率ωm＝ωb之间的频率：
至此，本发明可以确定好各个扫描频率，依次扫描各个频点，可以得到相应的相位响应θ(ω1)，θ(ω2)…θ(ωm)。由于采用非线性扫频，扫频步进呈指数级上升，所以需要的频点数量大大减少，测量速度可获得较大的提升，其相应的时延解算算法如下：
首先，利用经典相位展开算法计算出第二个频点的展开相位φ(ω2)，再按照如下公式计算出剩余频点的展开相位：
其中，[…]为取整符号，此处为“四舍五入”取整。由此可解算出终止频率的整周模糊度：
其中，[…]为取整运算符，此处为“四舍五入”取整。最后，可由公式(10)得到待测光器件的时延。测量光纤时，可由公式(11)计算出长度。
本发明可直接采用现有基于相推法的高精度光器件时延测量系统的硬件部分，只需要对软件部分进行简单修改即可实现。