基于脉冲模板函数的光纤传感阵列的脉冲延时测量方法
【专利摘要】本发明公开了一种基于脉冲模板函数的光纤传感阵列的脉冲延时测量方法。本发明读取信号数据，确定第一个定位脉冲位置，结合脉冲周期数划分信号数据得到数据矩阵A，计算得到方差值向量c，根据脉冲的时间参数构建脉冲模板函数h，计算方差值向量c与脉冲模板函数h的相关值得到相关系数向量d，找出d中最大的K个极大值点得到相关系数向量极大值索引值向量p，对p求差分得到脉冲延时值；本发明的方法极大提高了处理速度；并且，对噪声的抵抗能力显著增强，可以在其他方法无法进行正常判断的情况下仍能准确的计算出脉冲延时值；同时指定输入的参数少，仅需要事先知道脉冲的宽度、上升时间、下降时间和阵列中的传感器数目，因此非常的简便、可靠。
【专利说明】
基于脉冲模板函数的光纤传感阵列的脉冲延时测量方法
技术领域
[0001] 本发明涉及光纤传感器，具体涉及一种基于脉冲模板函数的光纤传感阵列的脉冲 延时测量方法。
【背景技术】
[0002] 光纤传感器是一类的新型传感器，因其抗电磁干扰，耐腐蚀，传输距离远，可大规 模复用等优势在水声探测、油田开发、大型结构物检测等领域有着日益广泛的应用。为了进 一步发挥光纤传感器体积小、重量轻等优势，目前主流的应用趋势是尽可能将多个光纤传 感器用各种不同的手段进行复用，形成光纤传感阵列或者网络来感知外界信息，如在水声 探测领域形成的海底阵缆(0BC)，在石油勘探领域采用的VSP检波器阵列等，都是目前已经 商用的典型光纤传感阵列。主要的复用方法包括时分复用，空分复用，波分复用，频分复用 等，这其中又以时分复用效率最高，成本最低。在上述提到的这些应用中，都是以时分复用 为基础，结合波分复用和空分复用形成大规模的传感阵列或者网络。
[0003] 时分复用的方法是通过向光纤传感阵列中注入周期性的脉冲光，接收返回的一系 列脉冲光再进行数据处理。由于阵列上的不同传感器返回的光脉冲在时域上相互不重叠， 因此就可以从返回的一串光脉冲中将不同传感器对应的脉冲分离开，进而解调出其分别对 应的传感信息。在数据处理之前，一个很重要的工作就是需要测定阵列上每一个传感器对 应的脉冲延时值，只有获知准确的脉冲延时值，才能精确地分离出每个传感器对应的脉冲。 由于在每一个光纤传感阵列中，光纤的长度均有所差异，因此，对于每一个时分复用的光纤 传感阵列，其脉冲延时值都有所不同，都需要在使用之前仔细测定其脉冲延时值。
[0004] 在以往的应用中，时分复用的光纤传感阵列的脉冲延时值基本上位手动测量，即 通过人为采集、观察、计算返回的光信号数据进行脉冲延时判定。即使存在自动测量方法， 也十分依赖信号的质量，对信号的幅度时分敏感，不够稳健。总之现有的方案存在如下问 题:效率低下，过程繁琐；同时很容易受到噪声的干扰导致误判，即对信号的信噪比要求较 高。因此，在大规模的光纤传感阵列中，需要开发简便易行，可靠度高的脉冲延时自动测量 方法。
【发明内容】
[0005] 针对以上现有技术中存在的问题，本发明提出了一种基于脉冲模板函数的光纤传 感阵列的脉冲延时测量方法，能准确的计算出脉冲延时值;同时简便、可靠。
[0006] 对于时分复用的光纤传感阵列，结合某阵列返回的光脉冲信号，首先给出以下基 本术语：
[0007] 定位脉冲：用于确定脉冲延时值的基准脉冲，不包含在光纤传感阵列返回的脉冲 序列中，是系统主动产生并添加进接收信号中的，该脉冲的幅值大于正常的脉冲序列中所 有的数值。
[0008] 脉冲周期:相邻两个定位脉冲之前的时间是一个脉冲周期，在一个传感阵列中，这 个脉冲周期是事先确定且不变的。
[0009] 脉冲序列:光纤传感阵列返回的一系列脉冲，会周期性出现，出现的周期值等同于 脉冲周期，即在一个脉冲周期内包含一个脉冲序列。其中不发生干涉现象的脉冲为非干涉 脉冲，发生干涉现象的为干涉脉冲。一般非干涉脉冲的幅度值在不同脉冲周期内维持不变， 干涉脉冲的幅度值会发生明显的变化，因为其携带了传感信息。脉冲延时值测量只针对干 涉脉冲。
[0010] 脉冲延时：各干涉脉冲中之间的时间差值，其中第一个脉冲延时值为第一个干涉 脉冲距离定位脉冲的距离。
[0011] 本发明的基于脉冲模板函数的光纤传感阵列的脉冲延时测量方法，包括以下步 骤：
[0012] (1)读取信号数据:从已经保存的文本文件中，或者从上位机接收的数据缓冲区中 光纤传感阵列的信号数据，包括多个脉冲周期，信号数据中包括添加进入的定位脉冲数据 和光纤传感阵列返回的脉冲序列数据；
[0013] (2)确定第一个定位脉冲位置:在信号数据中找出第一个数值最大且固定不变的 脉冲，即为第一个定位脉冲；
[0014] (3)得到数据矩阵:根据第一个定位脉冲的位置和已知的脉冲周期，将读取的信号 数据划分到一个m*n的矩阵中，记为数据矩阵A，其中m为两个定位脉冲之间的数据点数，η为 读取的信号数据中所包含的完整的脉冲周期的个数，则每一个行向量代表距离定位脉冲特 定距离数据点在各个脉冲周期内的变化值，记为行向量a 1;每一个列向量代表一个脉冲周 期内的数据点，记为列向量bj，其中，? = 0，1，2···ηι-1，j = 0，1，2···η-1;
[0015] (4)得到方差值向量:计算每一个行向量&1的方差值，得到由方差值构成的方差值 向量c，其中c的元素Ci = var(ai);
[0016] (5)构建脉冲模板函数：已知光纤传感阵列的脉冲的时间长度τ，脉冲的上升时间 ^，脉冲的下降时间1￡，在数字系统中将参数1，^和^按照数字系统的采样率离散化成1'，1'1· 和T f:
[0020]其中，round代表四舍五入操作，匕为系统的采样率，然后构建脉冲模板函数h，h中 的元素In的表达式为：
[0022] 其中，1为整数且0彡1ST;
[0023] (6)计算得到相关系数向量:从方差值向量c的第一元素开始，依次计算c中的元素 与脉冲模板函数h的相关值，得到相关系数向量d，相关系数向量d中的元素di表达式为：
[0025] (7)得到相关系数向量极大值索引值向量:根据时分复用光纤传感阵列中所包含 的传感器的数目K，找出相关系数向量中最大的K个极大值点，将这些极大值点在相关系数 向量d中索引值按升序排列成一个新的1XK维的向量，从而得到相关系数向量极大值索引 值向量P;
[0026] (8)得到脉冲延时值向量:对p求差分，可得到K-1个脉冲延时值，第一个值为p的第 一个元素加上T/2，即脉冲延时值构成的向量q中的元素表达式为：
[0028] 其中，round代表四舍五入操作，pk为相关系数向量极大值索引值向量p中的元素， 1? = 0"···，Κ-1〇
[0029] 其中，在步骤(1)中，一般要求读取的信号数据中至少包含20个以上的连续脉冲周 期。
[0030] 在步骤(2)中，信号数据中会包含定位脉冲，此定位脉冲并不是光纤传感阵列实际 返回的信号，而是人为添加进去的，作用在于辅助确定脉冲延时值的数值;一般定位脉冲的 幅度值会远远大于正常的信号数据的数值，为最大值且固定不变，因此第一个定位脉冲的 位置非常容易确定。
[0031] 在步骤(5)中，脉冲的上升时间为从最大值的5%~20%上升至80%~90%，脉冲 的下降时间为从最高值的80 %~90 %下降至5 %~20 %。
[0032]本发明的优点：
[0033]本发明读取信号数据，确定第一个定位脉冲位置，结合脉冲周期数划分信号数据 得到数据矩阵Α，计算得到方差值向量c，根据脉冲的时间参数构建脉冲模板函数h，计算方 差值向量c与脉冲模板函数h的相关值得到相关系数向量d，找出d中最大的K个极大值点得 到相关系数向量极大值索引值向量P，对P求插分得到脉冲延时值;与手动测试方法相比，本 发明的方法处理速度大为提高；同其他的脉冲延时值测量方法相比，本发明的方法对噪声 的抵抗能力显著增强，可以在其他方法无法进行正常判断的情况下仍能准确的计算出脉冲 延时值；同时需要指定输入的参数少，仅需要事先知道脉冲的宽度，上升时间，下降时间以 及阵列中的传感器数目，因此非常的简便、可靠。
【附图说明】
[0034] 图1为根据本发明的基于脉冲模板函数的光纤传感阵列的脉冲延时测量方法读取 信号数据的示意图；
[0035] 图2为本发明的基于脉冲模板函数的光纤传感阵列的脉冲延时测量方法的流程 图。
【具体实施方式】
[0036]下面结合附图，通过具体实施例，进一步阐述本发明。
[0037]在本实施例中，光纤传感阵列包括8个传感器;读取65组脉冲序列，在一组脉冲序 列中首尾两个脉冲是非干涉脉冲，其余为干涉脉冲。
[0038] 如图2所示，本实施例的基于脉冲模板函数的光纤传感阵列的脉冲延时测量方法， 包括以下步骤：
[0039] (1)读取信号数据:从已经保存的文本文件中读取信号数据，包括8个脉冲周期，如 图1所示；
[0040] (2)确定第一个定位脉冲位置:在信号数据中找出第一个数值最大且固定不变的 脉冲，即为第一个定位脉冲；
[0041] (3)得到数据矩阵:根据第一个定位脉冲的位置和已知的脉冲周期，将读取信号数 据划分到一个m*n的矩阵中，记为数据矩阵A，其中m为两个定位脉冲之间的数据点数，η为读 取的信号数据中所包含的完整的脉冲周期的个数，ri = 8,则每一个行向量代表距离定位脉 冲特定距离数据点在各个脉冲周期内的变化值，记为行向量a1;每一个列向量代表一个脉 冲周期内的数据点，记为列向量bj，其中，? = 0，1，2···πι-1，」_ = 0，1，2···7;
[0042] (4)得到方差值向量:计算每一个行向量&1的方差值，得到由方差值构成的方差值 向量c，其中c的元素 Ci = var(ai);
[0043] (5)构建脉冲模板函数：已知光纤传感阵列的脉冲的时间长度τ，脉冲的上升时间 为从最高值的10%上升至90% 脉冲的下降时间为从最高值的90%下降至10%τ?，τ = 140ns，Tr = 30ns，Tf = 30ns，在数字系统中将参数τ，1^和1^按照数字系统的采样率离散化成 T，T#Tf:
[0047]其中，round代表四舍五入操作，匕为系统的采样率，构建脉冲模板函数h，h中的元 素 hi的表达式为：
[0049] 其中，1为整数且0彡1ST;
[0050] (6)计算得到相关系数向量:从方差值向量c的第一元素开始，依次计算c中的元素 与脉冲模板函数h的相关值，得到相关系数向量d，相关系数向量d中的元素 di表达式为：
[0052] (7)得到相关系数向量极大值索引值向量:根据时分复用光纤传感阵列中所包含 的传感器的数目K，K = 8,找出相关系数向量中最大的8个极大值点，将这些极大值点在相关 系数向量d中索引值按升序排列成一个新的1X8维的向量，从而得到相关系数向量极大值 索引值向量P;
[0053] (8)得到脉冲延时值向量:对p求差分，可得到7个脉冲延时值，第一个值为p的第一 个元素加上T/2，即脉冲延时值构成的向量q中的元素表达式为：
[0055]其中，round代表四舍五入操作，pk为相关系数向量极大值索引值向量p中的元素， k = 0,1···，7〇
[0056] 本实施例在Matlab程序运行结果，win7 64位系统，Intel i5处理器，最高主频 2.7GHz，计算出8个延时参数，时间不超过0.5s。如果是采用该平台下的Labview程序，贝1J可 以做到根据返回的数据每隔l〇ms对参数进行更新和调整。
[0057]最后需要注意的是，公布实施例的目的在于帮助进一步理解本发明，但是本领域 的技术人员可以理解:在不脱离本发明及所附的权利要求的精神和范围内，各种替换和修 改都是可能的。因此，本发明不应局限于实施例所公开的内容，本发明要求保护的范围以权 利要求书界定的范围为准。
【主权项】
1. 一种基于脉冲模板函数的光纤传感阵列的脉冲延时测量方法，其特征在于，所述测 量方法包括W下步骤： (1) 读取信号数据:从已经保存的文本文件中，或者从上位机接收的数据缓冲区中光纤 传感阵列的信号数据，包括多个脉冲周期，信号数据中包括添加进入的定位脉冲数据和光 纤传感阵列返回的脉冲序列数据； (2) 确定第一个定位脉冲位置:在信号数据中找出第一个数值最大且固定不变的脉冲， 即为第一个定位脉冲； (3) 得到数据矩阵:根据第一个定位脉冲的位置和已知的脉冲周期，将读取的信号数据 划分到一个m*n的矩阵中，记为数据矩阵A，其中m为两个定位脉冲之间的数据点数，η为读取 的信号数据中所包含的完整的脉冲周期的个数，则每一个行向量代表距离定位脉冲特定距 离数据点在各个脉冲周期内的变化值，记为行向量ai;每一个列向量代表一个脉冲周期内 的数据点，记为列向量bj，其中，i = 0，l，2…m-1，j = 0，1，2…n-1; (4) 得到方差值向量:计算每一个行向量ai的方差值，得到由方差值构成的方差值向量 C，其中C的元素 Ci = var(ai); (5) 构建脉冲模板函数：已知光纤传感阵列的脉冲的时间长度τ，脉冲的上升时间Tr，脉 冲的下降时间Tf，在数字系统中将参数T，Tr和Tf按照数字系统的采样率离散化成T，Tr和Tf:其中，round代表四舍五入操作，fs为系统的采样率，然后构建脉冲模板函数h，h中的元 素 hi的表达式为：其中，1为整数且〇《1《Τ; (6) 计算得到相关系数向量:从方差值向量C的第一元素开始，依次计算C中的元素与脉 冲模板函数h的相关值，得到相关系数向量d，相关系数向量d中的元素 di表达式为：(7) 得到相关系数向量极大值索引值向量:根据时分复用光纤传感阵列中所包含的传 感器的数目K，找出相关系数向量中最大的K个极大值点，将运些极大值点在相关系数向量d 中索引值按升序排列成一个新的1XK维的向量，从而得到相关系数向量极大值索引值向量 P; (8)得到脉冲延时值向量:对p求差分，可得到Κ-1个脉冲延时值，第一个值为p的第一个 元素加上T/2，即脉冲延时值构成的向量q中的元素表达式为：其中，round代表四舍五入操作，pk为相关系数向量极大值索引值向量P中的元素 ，k = 0，2. 如权利要求1所述的测量方法，其特征在于，在步骤(1)中，读取的信号数据中至少包 含20个W上的连续脉冲周期。3. 如权利要求1所述的测量方法，其特征在于，在步骤（2)中，信号数据中包含定位脉 冲，此定位脉冲并不是光纤传感阵列实际返回的信号，而是人为添加进去的，作用在于辅助 确定脉冲延时值的数值;定位脉冲的幅度值会远远大于正常的信号数据的数值，为最大值 且固定不变。4. 如权利要求1所述的测量方法，其特征在于，在步骤(5)中，脉冲的上升时间为从最大 值的5%~20%上升至80%~90%，脉冲的下降时间为从最高值的80%~90%下降至5%~ 20%。
【文档编号】G01J11/00GK106092340SQ201610564506
【公开日】2016年11月9日
【申请日】2016年7月18日
【发明人】张敏, 刘飞, 邱晓康
【申请人】北京大学, 清华大学