基于脉冲模板函数的光纤传感阵列的脉冲延时测量方法技术

本发明专利技术公开了一种基于脉冲模板函数的光纤传感阵列的脉冲延时测量方法。本发明专利技术读取信号数据，确定第一个定位脉冲位置，结合脉冲周期数划分信号数据得到数据矩阵A，计算得到方差值向量c，根据脉冲的时间参数构建脉冲模板函数h，计算方差值向量c与脉冲模板函数h的相关值得到相关系数向量d，找出d中最大的K个极大值点得到相关系数向量极大值索引值向量p，对p求差分得到脉冲延时值；本发明专利技术的方法极大提高了处理速度；并且，对噪声的抵抗能力显著增强，可以在其他方法无法进行正常判断的情况下仍能准确的计算出脉冲延时值；同时指定输入的参数少，仅需要事先知道脉冲的宽度、上升时间、下降时间和阵列中的传感器数目，因此非常的简便、可靠。

【技术实现步骤摘要】

本专利技术涉及光纤传感器，具体涉及一种基于脉冲模板函数的光纤传感阵列的脉冲延时测量方法。
技术介绍
光纤传感器是一类的新型传感器，因其抗电磁干扰，耐腐蚀，传输距离远，可大规模复用等优势在水声探测、油田开发、大型结构物检测等领域有着日益广泛的应用。为了进一步发挥光纤传感器体积小、重量轻等优势，目前主流的应用趋势是尽可能将多个光纤传感器用各种不同的手段进行复用，形成光纤传感阵列或者网络来感知外界信息，如在水声探测领域形成的海底阵缆(OBC)，在石油勘探领域采用的VSP检波器阵列等，都是目前已经商用的典型光纤传感阵列。主要的复用方法包括时分复用，空分复用，波分复用，频分复用等，这其中又以时分复用效率最高，成本最低。在上述提到的这些应用中，都是以时分复用为基础，结合波分复用和空分复用形成大规模的传感阵列或者网络。时分复用的方法是通过向光纤传感阵列中注入周期性的脉冲光，接收返回的一系列脉冲光再进行数据处理。由于阵列上的不同传感器返回的光脉冲在时域上相互不重叠，因此就可以从返回的一串光脉冲中将不同传感器对应的脉冲分离开，进而解调出其分别对应的传感信息。在数据处理之前，一个很重要的工作就是需要测定阵列上每一个传感器对应的脉冲延时值，只有获知准确的脉冲延时值，才能精确地分离出每个传感器对应的脉冲。由于在每一个光纤传感阵列中，光纤的长度均有所差异，因此，对于每一个时分复用的光纤传感阵列，其脉冲延时值都有所不同，都需要在使用之前仔细测定其脉冲延时值。在以往的应用中，时分复用的光纤传感阵列的脉冲延时值基本上位手动测量，即通过人为采集、观察、计算返回的光信号数据进行脉冲延时判定。即使存在自动测量方法，也十分依赖信号的质量，对信号的幅度时分敏感，不够稳健。总之现有的方案存在如下问题：效率低下，过程繁琐；同时很容易受到噪声的干扰导致误判，即对信号的信噪比要求较高。因此，在大规模的光纤传感阵列中，需要开发简便易行，可靠度高的脉冲延时自动测量方法。
技术实现思路
针对以上现有技术中存在的问题，本专利技术提出了一种基于脉冲模板函数的光纤传感阵列的脉冲延时测量方法，能准确的计算出脉冲延时值；同时简便、可靠。对于时分复用的光纤传感阵列，结合某阵列返回的光脉冲信号，首先给出以下基本术语：定位脉冲：用于确定脉冲延时值的基准脉冲，不包含在光纤传感阵列返回的脉冲序列中，是系统主动产生并添加进接收信号中的，该脉冲的幅值大于正常的脉冲序列中所有的数值。脉冲周期：相邻两个定位脉冲之前的时间是一个脉冲周期，在一个传感阵列中，这个脉冲周期是事先确定且不变的。脉冲序列：光纤传感阵列返回的一系列脉冲，会周期性出现，出现的周期值等同于脉冲周期，即在一个脉冲周期内包含一个脉冲序列。其中不发生干涉现象的脉冲为非干涉脉冲，发生干涉现象的为干涉脉冲。一般非干涉脉冲的幅度值在不同脉冲周期内维持不变，干涉脉冲的幅度值会发生明显的变化，因为其携带了传感信息。脉冲延时值测量只针对干涉脉冲。脉冲延时：各干涉脉冲中之间的时间差值，其中第一个脉冲延时值为第一个干涉脉冲距离定位脉冲的距离。本专利技术的基于脉冲模板函数的光纤传感阵列的脉冲延时测量方法，包括以下步骤：(1)读取信号数据：从已经保存的文本文件中，或者从上位机接收的数据缓冲区中光纤传感阵列的信号数据，包括多个脉冲周期，信号数据中包括添加进入的定位脉冲数据和光纤传感阵列返回的脉冲序列数据；(2)确定第一个定位脉冲位置：在信号数据中找出第一个数值最大且固定不变的脉冲，即为第一个定位脉冲；(3)得到数据矩阵：根据第一个定位脉冲的位置和已知的脉冲周期，将读取的信号数据划分到一个m*n的矩阵中，记为数据矩阵A，其中m为两个定位脉冲之间的数据点数，n为读取的信号数据中所包含的完整的脉冲周期的个数，则每一个行向量代表距离定位脉冲特定距离数据点在各个脉冲周期内的变化值，记为行向量ai；每一个列向量代表一个脉冲周期内的数据点，记为列向量bj，其中，i＝0,1,2…m-1，j＝0,1,2…n-1；(4)得到方差值向量：计算每一个行向量ai的方差值，得到由方差值构成的方差值向量c，其中c的元素ci＝var(ai)；(5)构建脉冲模板函数：已知光纤传感阵列的脉冲的时间长度τ，脉冲的上升时间τr，脉冲的下降时间τf，在数字系统中将参数τ，τr和τf按照数字系统的采样率离散化成T，Tr和Tf： T = r o u n d ( τ f s ) ]]> T r = r o u n d ( τ r f s ) ]]> T f = r o u n d ( τ f f s ) ]]>其中，round代表四舍五入操作，fs为系统的采样率，然后构建脉冲模板函数h，h中的元素hl的表达式为： h l = l T r ( 0 ≤ l < T r ) 1 ( T r ≤ l < T - T f ) 本文档来自技高网...

【技术保护点】
一种基于脉冲模板函数的光纤传感阵列的脉冲延时测量方法，其特征在于，所述测量方法包括以下步骤：(1)读取信号数据：从已经保存的文本文件中，或者从上位机接收的数据缓冲区中光纤传感阵列的信号数据，包括多个脉冲周期，信号数据中包括添加进入的定位脉冲数据和光纤传感阵列返回的脉冲序列数据；(2)确定第一个定位脉冲位置：在信号数据中找出第一个数值最大且固定不变的脉冲，即为第一个定位脉冲；(3)得到数据矩阵：根据第一个定位脉冲的位置和已知的脉冲周期，将读取的信号数据划分到一个m*n的矩阵中，记为数据矩阵A，其中m为两个定位脉冲之间的数据点数，n为读取的信号数据中所包含的完整的脉冲周期的个数，则每一个行向量代表距离定位脉冲特定距离数据点在各个脉冲周期内的变化值，记为行向量ai；每一个列向量代表一个脉冲周期内的数据点，记为列向量bj，其中，i＝0,1,2…m‑1，j＝0,1,2…n‑1；(4)得到方差值向量：计算每一个行向量ai的方差值，得到由方差值构成的方差值向量c，其中c的元素ci＝var(ai)；(5)构建脉冲模板函数：已知光纤传感阵列的脉冲的时间长度τ，脉冲的上升时间τr，脉冲的下降时间τf，在数字系统中将参数τ，τr和τf按照数字系统的采样率离散化成T，Tr和Tf：T=round(τfs)]]>Tr=round(τrfs)]]>Tf=round(τffs)]]>其中，round代表四舍五入操作，fs为系统的采样率，然后构建脉冲模板函数h，h中的元素hl的表达式为：hl=lTr(0≤l<Tr)1(Tr≤l<T-Tf)T-lTf(T-Tf≤l≤T)]]>其中，l为整数且0≤l≤T；(6)计算得到相关系数向量：从方差值向量c的第一元素开始，依次计算c中的元素与脉冲模板函数h的相关值，得到相关系数向量d，相关系数向量d中的元素di表达式为：di=Σl=0Tc(i+l)hl;]]>(7)得到相关系数向量极大值索引值向量：根据时分复用光纤传感阵列中所包含的传感器的数目K，找出相关系数向量中最大的K个极大值点，将这些极大值点在相关系数向量d中索引值按升序排列成一个新的1×K维的向量，从而得到相关系数向量极大值索引值向量p；(8)得到脉冲延时值向量：对p求差分，可得到K‑1个脉冲延时值，第一个值为p的第一个元素加上T/2，即脉冲延时值构成的向量q中的元素表达式为：qk=p0+round(T/2)(k=0)pk-pk-1(1≤k≤K-1)]]>其中，round代表四舍五入操作，pk为相关系数向量极大值索引值向量p中的元素，k＝0,1…,K‑1。...

【技术特征摘要】
1.一种基于脉冲模板函数的光纤传感阵列的脉冲延时测量方法，其特征在于，所述测量方法包括以下步骤：(1)读取信号数据：从已经保存的文本文件中，或者从上位机接收的数据缓冲区中光纤传感阵列的信号数据，包括多个脉冲周期，信号数据中包括添加进入的定位脉冲数据和光纤传感阵列返回的脉冲序列数据；(2)确定第一个定位脉冲位置：在信号数据中找出第一个数值最大且固定不变的脉冲，即为第一个定位脉冲；(3)得到数据矩阵：根据第一个定位脉冲的位置和已知的脉冲周期，将读取的信号数据划分到一个m*n的矩阵中，记为数据矩阵A，其中m为两个定位脉冲之间的数据点数，n为读取的信号数据中所包含的完整的脉冲周期的个数，则每一个行向量代表距离定位脉冲特定距离数据点在各个脉冲周期内的变化值，记为行向量ai；每一个列向量代表一个脉冲周期内的数据点，记为列向量bj，其中，i＝0,1,2…m-1，j＝0,1,2…n-1；(4)得到方差值向量：计算每一个行向量ai的方差值，得到由方差值构成的方差值向量c，其中c的元素ci＝var(ai)；(5)构建脉冲模板函数：已知光纤传感阵列的脉冲的时间长度τ，脉冲的上升时间τr，脉冲的下降时间τf，在数字系统中将参数τ，τr和τf按照数字系统的采样率离散化成T，Tr和Tf： T = r o u n d ( τ f s ) ]]> T r = r o u n d ( τ r f s ) ]]> T f = r o u n d ( τ f f s ) ]]>其中，round代表四舍五入操作，fs为系统的采样率，然后构建脉冲模板函数h，h中的元素hl的表达式为： h l = l T r ( 0 ≤ l < T r ) 1 ( T r ≤ l < T - T f ) T - l T f ...

【专利技术属性】
技术研发人员：张敏，刘飞，邱晓康，
申请(专利权)人：北京大学，清华大学，
类型：发明
国别省市：北京;11