在干涉型传感系统中运动补偿的方法及装置制造方法

在干涉型传感系统中运动补偿的方法及装置制造方法
【专利摘要】光学解调系统例如基于OFDR的系统测量经受时变干扰的传感光波导的折射率的局部变化。为一定长度的传感光波导检测的干涉测量信号被变换到谱域。从变换的干涉测量数据集确定时变信号。从时变信号确定补偿信号，该补偿信号用于为干涉测量数据集补偿时变干扰。通过使用求平均和应变补偿实现进一步的鲁棒性。该补偿技术可以沿光波导的长度被应用。
【专利说明】在干涉型传感系统中运动补偿的方法及装置
[0001] 本申请要求2011年12月5日提交的美国临时专利申请序列号61/566, 860和2012 年3月23日提交的美国临时专利申请序列号61/614, 662的优先权，这些申请的内容通过 引用并入本文。

【技术领域】
[0002] 本发明涉及干涉型传感应用。一种示例应用为光频域反射计（0FDR)传感应用。

【背景技术】
[0003] 光频域反射计（0FDR)已经证明是在光纤中以高空间分辨率来测量应变的有效系 统。例如，参见美国专利6545760、6566648、5798521和7538883。该高分辨率性能已经被 证明在形状传感应用中非常有用。例如，参见美国专利777254U7781724和美国专利申请 20110109898。在美国专利申请20110247427中描述了在单芯光纤中的简单的应变传感应 用。
[0004] 0FDR通过在激光调谐的宽频率范围内获取数据来实现高空间分辨率。在大多数应 用中，该调谐随着激光在一定频率范围内进行扫频而发生。在使用该技术时，假设在扫频的 时间期间测试光纤是静态的或者无变化的。然而，如果在激光扫频的时间期间正在被测量 的系统发生变化，那么所得到的测量结果可能质量降低。此类改变可以归因于测试光纤的 移动或者连接测试光纤与仪器的光纤引线的移动。在存在时变变化里实现更高质量的测量 将是可取的。


【发明内容】

[0005] 示例的实施例包括用于测量经受时变干扰的传感光波导的参数的方法及光学解 调/询问（interrogation)系统。示例的时变干扰是传感光波导的运动。光学解调系统包 括光学干涉解调器和光学检测电路，该光学检测电路耦合到光学干涉解调器，用于为一定 长度的传感光波导检测光学干涉测量信号。数据处理电路接收来自光学检测电路的干涉测 量信号，并且产生该一定长度的传感光波导的干涉测量数据集。该干涉测量数据集被变换 到谱域，并且从变换的干涉测量数据集确定时变信号。补偿信号从时变信号确定，并且该补 偿信号被用于为干涉测量数据集补偿时域干扰，以提高参数的测量。
[0006] 在一个示例实施方式中，光学解调系统是基于光频域反射的系统，该系统包括解 调光源，并且其中光学干涉测量信号表示沿传感光波导作为时间的函数的背散射幅度。
[0007] 在非限制实施例中，处理电路通过将干涉测量数据集与干涉参考数据集进行比较 来确定时变信号。为干涉测量数据集中的反射事件确定干涉测量数据集，并且接收到的干 涉数据在反射事件周围被加窗。时变信号可以是相位信号，在此情况下，数据处理电路被配 置成通过展开（unwrap)相位信号而从变换的干涉测量数据集中提取相位信号，并且通过 从展开的相位信号去除线性拟合确定描述对传感器的干扰的非线性信号。在该示例中，非 线性信号是补偿信号。然后，数据处理电路可以从干涉测量数据集减去非线性信号，以便为 干涉测量数据集补偿时变干扰。
[0008] 在另外的示例的实施例中，时变信号是相位信号，数据处理电路被配置成通过与 谱域中的基线数据比较而从谱域中的变换的干涉测量数据中提取相位信号。
[0009] 在其他的示例实施例中，数据处理电路被配置成：将干涉测量数据集分为多个测 量数据段；变换每个段；将变换的测量段的多个与对应的参考谱数据段结合；对结合的段 求平均；从平均的段中确定相位响应；以及基于相位响应确定时变信号。此外，数据处理电 路可以被配置成确定施加于测量数据段的应变，其中该应变表明在相邻的测量数据段的响 应之间的未对准的量，并且使用测量的应变来对准谱域中的测量数据段的响应。
[0010] 在另外的示例实施例中，数据处理电路被配置成比较谱域中的每个数据段与对应 的参考数据段；确定时间延迟，该时间延迟表明在时域中的参考数据段与数据段之间的未 对准的量；以及使用测量的时间延迟对准时域中的参考数据段和数据段。
[0011] 在传感光波导为光纤的示例中，其中数据处理电路可以被配置成为干涉测量数据 集补偿在光纤上的应变。
[0012] 在优选的但仍然为示例性的实施例中，数据处理电路被配置成确定沿着传感光波 导的运动、时间延迟以及应变，并且补偿沿着传感光波导的累积运动、时间不对准以及应 变。

【专利附图】

【附图说明】
[0013] 图1是基于光频域反射计（0FDR)的传感系统的非限制示例，其中传感器经历时变 干扰。
[0014] 图2使用0FDR系统的非限制的、示例的分布式测量过程的流程图。
[0015] 图3是使用0FDR测量的理想的单反射峰的理论幅度与光学延迟的图。
[0016] 图4是被展开（unwrapped)后的显示线性相位响应的图3所示的相位的图。
[0017] 图5是来自单个反射器的幅度与延迟的图，在单个反射器中在获取扫频期间相位 已因一些时变变化而产生失真。
[0018] 图6是图5所示的失真峰的相位与光频的图。
[0019] 图7是示出在点A处测量的校正能够为所有后续点校正在测量点A和仪器之间的 受干扰的非限制性示例图。
[0020] 图8是来自两个反射事件的幅度与延迟的图，在两个反射事件中两个峰在光频域 中都因一些时变相位而发生失真。
[0021] 图9是在图8中索引50处反射峰的周围加窗的图。
[0022] 图10示出在图9中仅描述第一反射峰的信息的窗口。
[0023] 图11描绘了从图10中的加窗数据的傅里叶变换计算的相位失真。
[0024] 图12是以灰色描绘在延迟域中的原始的失真的反射峰以及以黑色描绘数据已通 过使用相位校正而被校正的峰的图表。
[0025] 图13是在0FDR测量中用于提取和随后校正时变信号的非限制性示例过程的流程 图。
[0026] 图14是用于在参考0FDR数据集和测量0FDR数据集之间的谱域中提取时变相位 响应和对其求平均的非限制的示例过程的流程图。
[0027] 图15图示模拟散射幅度与延迟的示例的图。
[0028] 图16图示被施加到散射数据的示例的相位失真的图。
[0029] 图17是以黑色迹线示出原始的散射幅度与散射模式的延迟以及以灰色迹线示出 在频率域中已因图16的相位失真而产生失真之后的散射幅度的图。
[0030] 图18是从散射模式的单个延迟域段计算的谱相位失真的图。
[0031] 图19是与施加的相位失真（灰色）相比较通过对谱域中每个段的复数数据求平 均而计算出的相位失真与频率的图。
[0032] 图20将插值计算的相位失真与原始的施加的相位失真进行对比。
[0033] 图21示出了原始的散射模式和校正的散射模式的幅度。
[0034] 图22描绘了施加于光纤引线的4种不同振动的相位失真与光频。
[0035] 图23描绘了在具有和不具有振动校正的情况下以300HZ振动光纤之后沿传感光 纤所计算的应变。
[0036] 图24描绘了将光纤覆盖在电风扇上的情况下沿传感光纤所计算的应变。
[0037] 图25描绘了作为一定长度的无应变的光纤的波长的函数为多个部分所计算的失 真。
[0038] 图26示出与原始施加的相位失真相比较通过对图30中所示的数据求平均而计算 的相位失真。
[0039] 图27示出作为一定长度的应变的光纤的波长的函数为多个部分所计算的相位失 真。
[0040] 图28示出作为一定长度的应变的光纤的波长的函数为多个部分所计算的相位失 真，其中已为每段校正施加的应变。
[0041] 图29是示出去除在用于提取运动校正信号的段上的均匀轴向应变的作用的示例 步骤的流程图。
[0042] 图30是相位与延迟的图，其示出在延迟域中的参考与测量结果之间的复数相乘 的提取幅角（argument)如何提供由于施加到测量的应变而产生的光延迟变化的测量。
[0043] 图31是相位与延迟的图，其示出从参考和校正的测量结果之间的复数相乘的幅 角所提取的光学延迟响应。
[0044] 图32图示了示例性反馈过程，其可以沿着一定长度的光纤执行，其中为该一定长 度的光纤校正在参考和测量段之间的累积的振动、累积的应变以及延迟失配。

【具体实施方式】
[0045] 下列描述出于解释而非限制性目的提出具体细节，诸如特定的实施例。但是本领 域技术人员将理解的是，偏离这些具体细节的其他的实施例可以被使用。在一些实例中， 众所周知的方法、接口、电路和设备的详细描述被省略，以便不会因不必要的细节而模糊本 说明书。在一些附图中，示出了独立的块。本领域的技术人员将理解这些块的功能可以通 过使用独立的硬件电路、使用连同适合的可编程的数字微处理器或者通用计算机的软件程 序和数据、和/或使用专用集成电路（ASIC)和/或使用一个或者更多个的数字信号处理器 (DSP)而实现。软件程序指令和数据可以被存储在非暂时性的计算机可读存储介质上，并且 当指令被计算机或者其他的适合的处理器控制执行时，计算机或者处理器执行与这些指令 有关的功能。
[0046] 因此，例如，本领域技术人员将理解本文中的附图能够表现说明性的电路或者其 他功能单元的概念上的视图。同样地，将被理解的是，任何流程图、状态转换图、伪代码等表 现了各种过程，这些过程可以大致被体现在计算机可读介质并且由此被计算机或者处理器 执行，不管这些计算机或者处理器是否被明确的示出。
[0047] 各种说明的块可以通过使用硬件诸如电路硬件和/或能够执行以存储在计算机 可读介质中的编码指令的方式的软件的硬件而被提供。因此，这些功能和说明的功能块将 被理解为硬件实现的和/或计算机实现的，以及机器实现的。
[0048] 关于硬件实现，功能块可以包括或者包含但不限于数字信号处理器（DSP)硬件、 精简指令集处理器、包括但不限于（多个）专用集成电路（ASIC)和/或（多个）现场可编 程门阵列（FPGA)的硬件（例如，数字或者模拟）电路以及（在适合的地方）能构执行这些 功能的状态机。
[0049] 关于计算机实现，计算机一般被理解为包括一个或更多个的处理器或者一个或更 多个的控制器，并且术语计算机、处理器和控制器可以交换使用。当由计算机、处理器或者 控制器提供功能时，功能可以通过单个的专用计算机或者处理器或者控制器、通过单个的 共享的计算机或者处理器或者控制器，或者通过多个独立的计算机或者处理器或者控制器 (它们中一些可以是共享的或者分布式的）来提供。此外，术语"处理器"或者"控制器"也 指其他能够执行这些功能和/或执行软件的硬件，诸如上面列举的示例的硬件。
[0050] 在本应用中的技术提供了用于补偿传感器运动的基于0FDR的方法和装置。术语 传感器包括基于0FDR的测量技术可以被应用到其中的任何波导。使得0FDR测量发生失真 的时变信号包括传感器的任何形式的运动，其发生在时间量程上，可与0FDR测量相当。振 动是一种非限制的示例，而突然的移动是另一个非限制性示例。在0FDR测量中，恒定的频 率和幅度的振动具有明显的轮廓，而在0FDR测量期间试验台的震动具有不明确的轮廓。在 这两种情况下，运动可能对于预期的测量的准确度是不利的。因此，发明人意识到必须开发 一种方法来处理广范范围的运动，在配置的传感环境中其能够被给定的传感器感受到。术 语"运动"被用于包含0FDR测量期间改变传感光纤的位置的任何类型的误差源。
[0051] 0FDR在执行沿波导长度的散射分布的高分辨率分布式测量时是高效的。在给定位 置处，沿波导的光的散射与局部的折射率有关。可以比较两种连续的测量以通过检测散射 分布的变化来检测沿波导长度的折射率的局部变化。
[0052] 图1是0FDR系统20的非限制的示例的设置，该0FDR系统20被用于监测沿光纤 传感器10的长度的折射率的局部变化，该光纤传感器10在一个或更多个测量和/或传感 应用中是有用的。在一些应用中，光纤传感器起传感器的作用，并且在其他的应用中它可以 是待测设备（DUT)或者其他实体。所示的时变干扰或者运动12影响光纤10 (术语光纤出 于方便的目的被使用，但该技术适用于任何适合的波导）的特定的位置。可谐调光源30在 光谱的范围内扫频。该光通过光耦合器的使用被分开，并且发送到两个单独的干涉仪。第 一干涉仪用作干涉解调器（interferometric interrogator) 34,其被连接到传感光纤10。 光通过干涉解调器34的测量臂进入传感光纤10。然后，沿光纤10长度的散射光被沿干涉 解调器34的参考臂传播的光干涉。激光监测网络32内的第二干涉仪在光源30扫描整个 频率范围时测量调谐率的波动。该激光监测网络32也包括氰化氢（HCN)气体单元，其在整 个测量扫描中被用于提供绝对波长参考。一系列的光学检测器40将来自激光监测网络32、 气体单元HCN的检测的光信号，以及来自传感光纤35的干涉模式转换成用于数据获取单元 36的电信号。数据处理器28使用来自数据获取单元36的获取的电信号以提取沿传感器 10长度的散射分布，如结合图2更详细解释的。
[0053] 图2是使用0FDR系统的非限制性的示例分布式测量过程的流程图。在步骤S1中， 可谐调的光源在光频的范围内扫频，并且经由干涉解调器的测量臂导向到传感器10中（步 骤S2)。沿传感器10长度的散射光与通过干涉解调器的参考路径传播的光干涉。为测量扫 频提供绝对波长参考（步骤S3)，并且测量调谐率变化（步骤S4)。光学检测器将检测的光 信号转换成电信号（步骤S5)以便由数据处理器28处理。传感光纤的干涉模式通过使用 激光监控信号被优选地重新采样，以保证检测的信号以光频的恒定增量被采样（步骤S6)。 一旦被重新采样，就执行傅里叶变换以产生时域上的传感器10散射信号。在时域中，散射 信号将散射事件的幅度描述成沿传感器10长度的延迟的函数（步骤S8)。使用光以给定的 时间增量传播的距离，该延迟被转换成沿传感器10长度的信号测量。也就是说，该信号将 每个散射事件描述成沿传感器10的距离的函数。采样周期被称作空间分辨率，并且采样周 期与测量期间可谐调的光源扫过的频率范围成反比。
[0054] 由于传感光纤10应变，局部的散射随着光纤在物理长度上的变化而移位。可以证 明的是这些失真是高度地可重复的。因此，0FDR测量结果可以被保留在内存中，用作在无 应变状态下的光纤的参考模式。随后的测量结果可以与该参考模式进行比较，以获得沿传 感光纤长度的局部散射的延迟移位的测量。当与参考散射模式比较时，这种延迟的移位表 示为连续的、缓慢变化的光学相位信号。该光学相位信号的导数与传感芯的物理长度的变 化成正比。物理长度的变化可以与应变成比例，产生沿传感光纤10的应变的连续测量。
[0055] 如上详述的，在0FDR测量期间，随着激光扫过光频范围，在参考路径和测量路径 之间的光学干涉被记录。在系统的光学检测器处所得的干涉模式包含关于从传感光纤反射 的光的幅度和相位的信息。光的相位和幅度在激光扫频期间作为时间的函数被记录。系统 的基本假设是在解调/询问（interrogation)下的干涉系统（包括传感光纤10)在扫频期 间不变，以及当激光扫频时作为光频的函数的系统的相位和幅度响应在时间上被编码。因 此，如果激光频率已知是时间的函数，那么系统的相位和幅度响应已知为激光波长的函数。
[0056] 但是如果系统在测量期间不是静态的，那么测量会失真。例如，考虑当激光随时间 线性地扫频时单个反射的响应。单个反射事件导致检测电子器件处的正弦干涉条纹。来自 单个反射器的干涉信号的傅里叶变换提供作为光延迟的函数的单个干净峰。图3描绘了使 用0FDR测量的理想的单个反射峰的理论上的幅度与光学延迟。
[0057] 延迟中的反射的位置与谱域中相位斜率的大小成正比。进一步，按照定义，在时 域中来自单个事件的理想反射在谱域中具有线性相位斜率。例如，图4描绘了"展开的 (unwrapped)"谱域相位。展开（unwrapping)可以被认为是随着光沿着光纤移动，光学相位 增加。在复平面上，由于相位围绕单位圆增加和移动，相位旋转一周并且接下来继续。展开 相位通过跟踪完整旋转的总数而有效地使旋转的相位信号线性化，以获得总的光学相位变 化的测量。
[0058] 现考虑这种情况，其中在干涉仪中在事件的反射和光束重组之间的光学延迟改 变，即在激光扫频时相位已被一些时变变化扭曲。如图5所示，这导致在测量的干涉信号上 的额外的相位调制。在存在这种随时间变化的相位调制的情况下，如图5描述的，所得的干 涉信号的傅里叶变换失真。峰不再是"干净"的，在峰的底部周围有噪声出现。
[0059] 图6描绘了图5所示的失真的峰的相位与光频。该失真的峰的谱域相位表现出从 线性相位的小偏离。从线性的这些偏差表明在0FDR系统执行测量时传感光纤的位置不是 恒定的。从线性相位的这些偏差在谱域中响应，并且所得的理想的峰的失真对分布式的应 变测量来说是有问题的。如上所述，将分布式的0FDR应变测量与传感光纤的基线参考扫描 进行比较，例如，通过执行互相关。由于传感器的运动，测量的失真大大降低这种比较的准 确度。发明人认识到这个问题并且开发了用于补偿测量时的运动的影响的技术，以便提高 0FDR测量的准确度。
[0060] 作为解决方案的开始，发明人认识到0FDR测量的基本特性允许从0FDR测量高效 地并且有效地去除运动的影响。在0FDR获取期间，来自延迟域中的理想单个反射事件的在 谱域中的相位响应的大小与事件的延迟成比例。因此，相位调制必须是沿传感器的长度添 加的。也就是说，在光纤中任何点能看到的相位调制是在感兴趣的点之前光纤中的所有相 位调制的和。这意味着在光纤中的任一点处测量相位调制提供了校正沿光纤的剩余部分的 这些调制的能力。
[0061] 图7示出在点A处测量的校正能够为所有后续点校正测量点A处与0FDR仪器之 间的干扰的图。也就是说，通过测量传感光纤10中的点A处的运动，可以为光纤10的剩余 长度校正在点A和仪器之间的所有的运动影响。因此，点B处的测量将不再受点A之前的 运动失真的影响。这个概念在图8中进一步被说明，即通过观察由于在点A之前的时变失 真而引起的在延迟域中的点A和点B二者之间的反射失真。在光学频域中峰已被一些时变 相位扭曲。
[0062] 如果已知在索引50处或其周围仅有单个反射，那么与该反射有关的复数数据集， 即从0到100之间的点可以被用于表征引起该反射所发生的振动（相位调制）。从数据集 中去除这些振动影响去除了沿光纤10的点A处与点B处的两个反射事件中观察到的失真。
[0063] 在这种情况下，第一峰的振动信息可以通过在该峰的周围加窗获得，如图9中的 以索引50处的反射峰为中心的矩形窗口所示。相乘之后，没有其他峰的信息存留。回想在 这些前100个点中所包含的信息是包含了相位和幅度两种信息的复数值。在加窗操作之后 对应的幅度数据在图10中描绘，其中该窗与从两种反射所获取的数据相乘，产生仅来自第 一反射峰的复数数据集。对于加窗后的复数数据集进行傅里叶变换，以便提取在谱域中的 相位失真。该相位失真描述了与来自点A处的原始的理想的单个反射的相位响应相关的非 线性，并且在图11中进行描绘。
[0064] 然后，该相位失真被用作相位失真的校正，其可以通过从全部获得的数据集的测 量的相位中减去该相位失真而被应用到全部获得的数据集，从而在延迟域中恢复两个位置 处的干净的峰。这在图12中被示出，其中灰色线示出在延迟域中的原始的失真的反射峰， 黑色线示出通过使用相位校正对数据进行校正后的峰，该相位校正从在索引50处的峰计 算而来。现在两个峰值都是锐利的。
[0065] 在图13中示出了略述非限制性示例步骤以提取由于运动带来的时变信号和随后 的数据的校正的流程图。以步骤S10开始，对与反射事件有关的时间延迟域中的复数值 0FDR数据集例如在数据集的开始处执行加窗操作。应用傅里叶逆变换将数据集变换到频 域（步骤S12)。谱域中的相位响应通过展开相位信号而被提取（步骤S14)。全部的累积 相位与延迟域中的反射事件的位置成比例。在延迟域中的每个延迟索引累积谱域中的两个 Pi ( π )相位变化，例如在索引500处的反射事件将具有在展开后谱域中1000 Pi ( π )的全 部的累积的相位响应。从该累积的相位响应中减去线性拟合去除了与在延迟域中的反射事 件的位置有关的相位，提供了描述传感器运动的非线性时变信号的测量（步骤S16)。非线 性的时变信号是补偿信号，例如，将其从谱域中的原始测量OFDR数据集的相位中减掉（步 骤S18)。然后，运动补偿的谱响应通过使用傅里叶变换操作被变换到时间延迟域（步骤 S20)，这可以被用于识别和/或展现反射的干净的谱峰。
[0066] 该示例过程通过提取时变信号直至传感器光纤10中的一点并且随后使用该信号 以补偿沿着传感器10的长度所产生的数据，以此来去除或者至少大致减小对0FDR测量的 运动的影响。该技术补偿沿传感器光纤10长度的一个或者多个运动，并且保证在不同的 0FDR应用中配置传感器的广泛的环境范围中的鲁棒性。
[0067] 在0FDR传感应用中的常见类型的光纤是具有瑞利散射分布的工业标准、低弯曲 损耗的光纤。给定长度的传感器光纤的瑞利散射分布表现为随机的宽带信号，并且由从光 纤自身的玻璃分子散射的光造成。该散射模式是高度地可重复的，并且可以被认为是给定 传感器的物理特性。因此，这种反射"指纹"被用于获得光纤的基线或者参考测量。
[0068] 在单个反射事件的示例的情况下，从谱域中的相位响应中去除线性拟合，以隔离 非线性分量，该非线性分量捕获测量中的运动的影响。在实际的系统中，理想的单个反射事 件一般不发生，并且相邻的反射事件的谱响应将掩饰期望的非线性分量。这能够通过比较 光纤的测量和相同光纤的基线测量来克服。光纤的散射分布是高度地可重复的，并且谱域 中的测量段与谱域中的对应段的基线扫描的比较将允许隔离测量期间的运动影响。也就是 说，如果在测量期间没有运动存在，则谱域中的测量与基线之间的比较将产生零值相位差。 如果存在运动，将看到非线性的响应。
[0069] 三个问题必须被克服。第一，瑞利散射是微弱信号，并且因此，对运动的补偿易受 到噪声的影响。发明人开发了一种平均方法以增加运动补偿的鲁棒性。第二，配置的传感 器易受到许多不同的环境因素的影响。鲁棒的运动补偿应该保证环境的变化（诸如温度变 化）不降低补偿的效率。发明人开发了一种从运动提取区域去除轴向应变的影响的方法。 第三，不能保证配置的传感器仅在传感光纤中感兴趣区域之前经受时变干扰。因此，补偿技 术应该沿传感器的长度操作，以便系统地去除沿传感器长度的各种运动影响。
[0070] 瑞利散射的低信号水平问题通过对光纤中的多个邻近段的散射使用平均操作来 解决。在一般的应用中，可以假设对于短长度的传感光纤来说，运动的影响是相似的。在延 迟域中选择小段的传感光纤允许这种假设，即相邻的段将观察到相似的时变失真。首先，时 变的谱失真从每个段中提取。将一段的测量0FDR数据和相同段的参考/基线0FDR数据两 者傅里叶变换至谱域，可以比较复数值谱数据以提取由于在测量段中的运动引起的相位失 真的测量。可以为每个相邻的段重复该过程，并且谱失真的平均值被计算。这种过程的示 例步骤被描述在图14中的流程图中。
[0071] 对基线环境中（例如，没有影响光纤的运动、控制温度等）的传感光纤10执行传 感光纤的0FDR扫描（步骤S33)，从而为光纤的多个段的每段产生0FDR参考数据（步骤 S34)，并且对应用环境中（例如，可能存在影响光纤的运动、温度变化等）的传感光纤10执 行传感光纤的OFDR扫描（步骤S30)，从而为光纤的多个段的每段产生OFDR测量数据（步 骤S31)。段可以通过使用例如上述的加窗技术来被分区。对每个测量段（步骤S32)和每 个参考段（步骤S35)执行傅里叶变换操作，以将段数据变换至谱域。执行傅里叶变换的一 种示例方法是使用快速傅里叶变换（FFT)。然后，每个段的测量谱数据与每个相应段的参考 谱数据的复数共轭进行复数相乘（数据是复数）（步骤S36)。在数学上，复数值信号与第二 信号的共轭相乘等同于计算复数信号的相位值的差。由于存在使该相位测量失真的噪声， 所以期望对相邻段求平均。为计算相邻段的平均值，来自相邻段的复数乘积的实部和虚部 分别被求平均（S37)。通过计算平均的复数乘积的实部和虚部的反正切来提取平均的相位 响应。该相位差如上述被随后展开。然后，该相位响应能够被用于补偿由于运动引起的失 真。例如，相位响应能够被线性插值以匹配测量结果的大小，并且从原始的测量数据的谱域 被减掉以去除运动的影响。
[0072] 虽然在步骤S30-S32和S33-S35中示出了两个平行的轨迹，但是这些步骤不需要 在同一时间并行执行或者即使一对一地执行。例如，参考扫描步骤S33-S35可以被提前执 行，并且参考段数据被存储在内存以随后用于测量的段数据。一种测量扫描被示出，但应当 理解的是步骤S30-S32和S36-S39可以对一个或更多个其他测量扫描重复进行。
[0073] 现从数学角度描述图14中所示的运动补偿的计算的处理程流。κ p为在一定长度 的均匀应变的光纤中的复数散射测量的数组，其中p表示在延迟中沿数组的位置。通过将 κ ρ的连续的段变换回到谱域来产生子数组集。Κκη代表测量段的谱，η和κ是谱索引：
[0074] KKn = FFT{KnS··· κ (n+1)s}，其中S是段的点数的长度。使用参考测量％重复该过 程，
[0075] N^FFTU，. n(n+1)s}
[0076] 然后，捕获运动影响的信号由下式给出
[0077]

【权利要求】
1. 一种用于测量经受时变干扰的传感光波导的参数的光学解调系统，包括 光学干涉解调器； 光学检测电路，其被耦合到所述光学干涉解调器，用于检测一定长度的所述传感光波 导的光学干涉测量信号；和 数据处理电路，其被配置成： 从所述光学检测电路接收干涉测量信号， 产生所述一定长度的传感光波导的干涉测量数据集； 将所述干涉测量数据集变换到谱域； 根据变换的干涉测量数据集确定时变信号；以及 根据所述时变信号确定补偿信号； 使用所述补偿信号为所述干涉测量数据集补偿所述时变干扰，以增强所述参数的测 量。
2. 根据权利要求1所述的光学解调系统，其中所述光学解调系统是包括解调光源的基 于光频域反射即OFDR的系统，并且其中所述光学干涉测量信号指示沿所述传感光波导的 作为时间的函数的背散射幅度。
3. 根据权利要求1所述的光学解调系统，其中所述处理电路被配置成通过比较所述干 涉测量数据集和干涉参考数据集来确定时变信号。
4. 根据权利要求1所述的光学解调系统，其中所述时变干扰包括所述传感光波导的运 动。
5. 根据权利要求1所述的光学解调系统，其中为所述干涉测量数据集中的反射事件确 定所述干涉测量数据集，以及 其中所述数据处理电路被配置成在所述反射事件的周围为接收的干涉数据加窗。
6. 根据权利要求1所述的光学解调系统，其中所述时变信号是相位信号，并且其中所 述数据处理电路被配置成： 通过展开所述相位信号从所述变换的干涉测量数据集提取所述相位信号；以及 通过从展开的相位信号去除线性拟合来确定描述对所述传感器的所述时变干扰的非 线性信号， 其中所述非线性信号是所述补偿信号。
7. 根据权利要求6所述的光学解调系统，其中所述数据处理电路被配置成从所述干涉 测量数据集减去所述非线性信号，以便为所述干涉测量数据集补偿所述时变干扰。
8. 根据权利要求1所述的光学解调系统，其中所述时变信号是相位信号，并且其中所 述数据处理电路被配置成： 通过比较在所述谱域中的基线数据集，从所述谱域中的变换的干涉测量数据中提取所 述相位信号。
9. 根据权利要求1所述的光学解调系统，其中所述数据处理电路被配置成： 将所述干涉测量数据集分成多个测量数据段； 变换每个所述段； 将变换的测量段中的多个与对应的参考谱数据段结合； 对结合的段求平均； 根据平均的段确定相位响应；以及 基于所述相位响应确定所述时变信号。
10. 根据权利要求1所述的光学解调系统，其中所述数据处理电路被配置成： 将所述干涉测量数据集分成多个测量数据段，以及 确定应用于测量数据段的应变，其中所述应变指示相邻的测量数据段的响应之间的未 对准的量，以及 使用测量的应变以对准所述谱域中的测量数据段的响应。
11. 根据权利要求1所述的光学解调系统，其中所述数据处理电路被配置成： 将所述干涉测量数据集分成多个测量数据段； 将每个数据段与对应的参考数据段比较； 确定时间延迟，所述时间延迟指示时域中的参考数据段与数据段之间的未对准的量； 以及 使用测量的时间延迟以对准所述时域中的参考数据段与数据段。
12. 根据权利要求1所述的光学解调系统，其中所述传感光波导是光纤，并且其中所述 数据处理电路被配置成： 为所述干涉测量数据集补偿光纤上的应变。
13. 根据权利要求1所述的光学解调系统，其中所述数据处理电路被配置成确定沿所 述传感光波导的运动、时间延迟和应变，并且补偿沿所述传感光波导的累积的运动、时间未 对准和应变。
14. 一种用于测量经受时变干扰的传感光波导的参数的方法，其包括： 检测一定长度的所述传感光波导的光学干涉测量信号； 根据检测的光学干涉测量信号产生所述一定长度的所述传感光波导的干涉测量数据 集； 将所述干涉测量数据集变换至所述谱域； 根据所述变换的干涉测量数据集确定时变信号； 根据所述时变信号确定补偿信号；以及 使用所述补偿信号为所述干涉测量数据集补偿所述时变干扰，以增强所述参数的测 量。
15. 根据权利要求14所述的方法，其中所述方法使用包括解调光源的基于光频域反射 即OFDR的系统，并且其中所述光学干涉测量信号指示沿所述传感光波导的作为时间的函 数的背散射幅度。
16. 根据权利要求14所述的方法，进一步包括通过比较所述干涉测量数据集与干涉参 考数据集来确定所述时变信号。
17. 根据权利要求14所述的方法，其中所述时变干扰包括所述传感光波导的运动。
18. 根据权利要求14所述的方法，其中为所述干涉测量数据集中的反射事件确定所述 干涉测量数据集，所述方法进一步包括在所述反射事件的周围为接收的干涉数据加窗。
19. 根据权利要求14所述的方法，其中所述时变信号是相位信号，并且所述方法进一 步包括： 通过展开所述相位信号从所述变换的干涉测量数据集提取所述相位信号；以及 通过从展开的相位信号去除线性拟合来确定描述对所述传感器的所述时变干扰的非 线性信号， 其中所述非线性信号是所述补偿信号。
20. 根据权利要求19所述的方法，进一步包括从所述干涉测量数据集中减去所述非线 性信号，以便为所述干涉测量数据集补偿所述时变干扰。
21. 根据权利要求14所述的方法，其中所述时变信号是相位信号，并且所述方法进一 步包括： 通过比较所述谱域中的基线数据集，从所述谱域中的所述变换的干涉测量数据中提取 所述相位信号。
22. 根据权利要求14所述的方法，进一步包括： 将所述干涉测量数据集分成多个测量数据段； 变换每个所述段； 将变换的测量段中的多个与对应的参考谱数据段结合； 对结合的段求平均； 从平均的段确定相位响应；以及 基于所述相位响应确定所述时变信号。
23. 根据权利要求14所述的方法，进一步包括： 将所述干涉测量数据集分成多个测量数据段，以及 确定应用于测量数据段的应变，其中所述应变指示相邻的测量数据段的响应之间的未 对准的量，以及 使用测量的应变对准所述谱域中的测量数据段的响应。
24. 根据权利要求14所述的方法，进一步包括： 将所述干涉测量数据集分成多个测量数据段； 将每个数据段与对应的参考数据段比较； 确定时间延迟，所述时间延迟指示时域中的参考数据段与数据段之间的未对准的量； 以及 使用测量的时间延迟对准所述时域中的参考数据段与数据段。
25. 根据权利要求14所述的方法，其中所述传感光波导是光纤，并且所述方法进一步 包括为所述干涉测量数据集补偿在所述光纤上的应变。
26. 根据权利要求14所述的方法，进一步包括： 确定沿所述传感光波导的运动、时间延迟和应变，并且 补偿沿所述传感光波导的累积的运动、时间未对准以及应变。
【文档编号】G01N21/45GK104126103SQ201280068990
【公开日】2014年10月29日 申请日期:2012年12月3日 优先权日:2011年12月5日
【发明者】M·E·佛罗凯特, A·K·桑, D·K·吉福德, J·W·克莱因 申请人:直观外科手术操作公司