CN103308083B - 一种基于能量分布的布里渊谱识别方法 - Google Patents

一种基于能量分布的布里渊谱识别方法 Download PDF

Info

Publication number
CN103308083B
CN103308083B CN201310264845.7A CN201310264845A CN103308083B CN 103308083 B CN103308083 B CN 103308083B CN 201310264845 A CN201310264845 A CN 201310264845A CN 103308083 B CN103308083 B CN 103308083B
Authority
CN
China
Prior art keywords
brillouin spectrum
energy
data
energy density
bandwidth
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
CN201310264845.7A
Other languages
English (en)
Other versions
CN103308083A (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.)
CLP Kesiyi Technology Co Ltd
Original Assignee
CETC 41 Institute
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 CETC 41 Institute filed Critical CETC 41 Institute
Priority to CN201310264845.7A priority Critical patent/CN103308083B/zh
Publication of CN103308083A publication Critical patent/CN103308083A/zh
Application granted granted Critical
Publication of CN103308083B publication Critical patent/CN103308083B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Length Measuring Devices By Optical Means (AREA)

Abstract

本发明提供一种基于能量分布的布里渊谱识别方法,包括以下步骤:确定布里渊谱测试数据为L;根据探测脉冲宽度PW选择相应的能量积分窗口宽度W,得到能量密度数据为P;在能量密度数据P中搜索能量密度最高的频率范围SW;在能量密度数据P中,在所述能量密度最高的频率范围SW内,计算其峰值是否超过能量阈值PT;在能量密度数据P中,搜索所述最高的频率范围内峰值对应的频率点。采用上述方案,一方面根据判断结果只对含有布里渊谱的数据进行曲线拟合和谱分析,避免了很多无谓的布里渊谱分析过程,节省整体数据输出时间,提高了整体测试速度,另一方面提高了各种信噪比情况下布里渊谱分析的准确性,提高了应变测试精度。

Description

一种基于能量分布的布里渊谱识别方法
技术领域
本发明属于布里渊谱识别技术领域,尤其涉及的是一种基于能量分布的布里渊谱识别方法。
背景技术
在光纤传感系统中,光纤应变分布测试仪(简称BOTDR)是不可或缺的测试仪器,广泛应用于桥梁、堤坝、建筑等结构的健康状态监测以及石油、煤矿、管道等的区域安全监测等领域。其基本原理是从光纤一端注入脉冲光并接收光纤沿线各点返回的自发布里渊散射光,经过时域采样及分析,即可得到沿线各点的布里渊散射光的频谱,进而经分析得到频移,然后再通过计算,即可得到光纤沿线各点的应变分布。
在BOTDR的内部数据处理中,通过扫频方式获得光纤沿线各点的布里渊谱分布数据后,需首先通过分析确定各点的布里渊谱中心位置,然后经复杂的洛仑兹曲线拟合及布里渊散射谱分析才得到光纤沿线各点的应变分布。在确定各点的布里渊谱中心位置过程中,通常的做法是将测试量程内采集的全部数据直接通过寻找最大值的方法,该方法有两方面的问题:一是由于电路噪声的干扰,对采集的数据通过寻找最大值的方法确定的布里渊谱中心位置有时会与真实位置偏差很大,尤其是在信噪比较差的情况下,该现象比较明显,这将导致后续的洛仑兹曲线拟合和布里渊谱分析结果误差很大,进而造成错误的应变计算结果;另一方面,由于在测试过程中,被测光纤长度小于测试量程,在当前测试量程内进行数据采集时,超过光纤长度后的采样数据都是高斯噪声,无任何有用信息,没有必要再进行复杂的洛仑兹曲线拟合及布里渊散射谱分析,只会占用很多无谓的数据处理时间,从而影响了整体测试速度。
因此,现有技术存在缺陷,需要改进。
发明内容
本发明所要解决的技术问题是针对现有技术的不足,提供一种布里渊谱自适应分析方法。
本发明的技术方案如下:
一种基于能量分布的布里渊谱识别方法,其中,包括以下步骤:
步骤100,确定布里渊谱测试数据为L;
步骤101,根据探测脉冲宽度PW选择相应的能量积分窗口宽度W,并根据L及W计算,得到能量密度数据为P;
步骤102,在能量密度数据P中搜索预定带宽DW条件下能量密度最高的频率范围SW,并将此频率范围提供给步骤103;
步骤103,在能量密度数据P中,在所述能量密度最高的频率范围SW内,计算其峰值是否超过能量阈值PT,如超过能量阈值PT则进入步骤104,如不超过则进入步骤108;
步骤104,在能量密度数据P中,搜索所述最高的频率范围内峰值对应的频率点;
步骤105,在布里渊谱数据L中,以所述峰值为中心频率计算3dB带宽,并与3dB带宽阈值比较,如超过阈值则进入步骤109,如未超过则进入步骤106;
步骤106,在布里渊谱数据L中,以所述峰值为中心频率计算6dB带宽,并与6dB带宽阈值比较,如超过阈值则进入步骤109,如未超过则进入步骤107;
步骤107,在布里渊谱数据L中,以所述峰值为中心频率计算9dB带宽,并与9dB带宽阈值值比较,如超过阈值则进入步骤109,如未超过则进入步骤120;
步骤108,未超过则认为此布里渊谱数据L中都为噪声,无布里渊谱存在;
步骤109,识别过程结束,输出中心频率及其在数据中的位置;
步骤120,结束,无布里渊谱存在。
所述的基于能量分布的布里渊谱识别方法,其中,所述步骤101中,所述脉冲宽度PW越窄,能量积分窗口宽度越宽。
所述的基于能量分布的布里渊谱识别方法,其中,所述步骤101中,所述布里渊谱测试数据为L根据积分窗口宽度W进行积分运算。
所述的基于能量分布的布里渊谱识别方法,其中,所述步骤102中,所述预定带宽DW的选择范围应小于能量积分窗口W。
所述的基于能量分布的布里渊谱识别方法,其中,所述步骤102中,所述预定带宽DW为能量积分窗口宽度W的1/2。
所述的基于能量分布的布里渊谱识别方法,其中,所述步骤103中,所述能量阈值PT是根据量程最远的100个点的布里渊谱测试数据L进行平均计算获取。
所述的基于能量分布的布里渊谱识别方法,其中,所述步骤105、步骤106及步骤107中,所述带宽阈值是根据布里渊谱数据L的探测脉冲宽度确定。
采用上述方案,能够对采集数据是否含有布里渊谱及其位置进行准确判断和定位,一方面根据判断结果只对含有布里渊谱的数据进行曲线拟合和谱分析,避免了很多无谓的布里渊谱分析过程,节省整体数据输出时间,提高了整体测试速度,另一方面提高了各种信噪比情况下布里渊谱分析的准确性,提高了应变测试精度。
附图说明
图1为本发明布里渊谱自适应分析方法流程图。
具体实施方式
以下结合附图和具体实施例,对本发明进行详细说明。
实施例1
如图1所示,本发明提供一种布里渊谱自适应分析方法,该方法包括以下实施步骤:
步骤100,确定布里渊谱测试数据为L(包括布里渊谱频率及其对应的功率值),其探测脉冲宽度PW为10ns;
步骤101,根据探测脉冲宽度PW选择相应的能量积分窗口宽度W,探测脉冲宽度越窄,能量积分窗口宽度越宽,其选择方案可如表1所示,当脉冲宽度PW为10ns时,可选择能量积分窗口宽度W为100Mhz,并根据L及W对布里渊谱测试数据L进行积分运算(即以能量积分窗口宽度W为积分宽度,对L内所有数据进行积分),得到能量密度数据为P;
表1:不同PW对应的能量积分窗口宽度W
脉冲宽度PW 能量积分窗口宽度W
5ns 200MHz
10ns 100Mhz
20ns 50MHz
50ns 25MHz
100ns 20MHz
200ns 15MHz
步骤102,在能量密度数据P中搜索预定带宽DW条件下能量密度最高的频率范围SW,并将此频率范围提供给步骤103;预定带宽DW的选择范围应小于能量积分窗口W,可以为能量积分窗口宽度W的1/2;
步骤103,在能量密度数据P中,在所述能量密度最高的频率范围SW内,计算其峰值是否超过能量阈值PT,如超过能量阈值PT则进入步骤104,如不超过则进入步骤108;能量阈值PT可以根据量程最远的100个点的布里渊谱测试数据L进行平均计算得到;
步骤104,在能量密度数据P中,搜索所述最高的频率范围内峰值对应的频率点;
步骤105,在布里渊谱数据L中,以所述峰值为中心频率计算3dB带宽,并与3dB带宽阈值比较,如超过阈值则进入步骤109,如未超过则进入步骤106;带宽阈值可以根据布里渊谱数据L的探测脉冲宽度确定;
步骤106,在布里渊谱数据L中,以所述峰值为中心频率计算6dB带宽,并与6dB带宽阈值比较,如超过阈值则进入步骤109,如未超过则进入步骤107;带宽阈值可以根据布里渊谱数据L的探测脉冲宽度确定;
步骤107,在布里渊谱数据L中,以所述峰值为中心频率计算9dB带宽,并与9dB带宽阈值值比较,如超过阈值则进入步骤109,如未超过则进入步骤120;带宽阈值可以根据布里渊谱数据L的探测脉冲宽度确定;
步骤108,未超过则认为此布里渊谱数据L中都为噪声,无布里渊谱存在;
步骤109,识别过程结束,输出中心频率及其在数据中的位置;
步骤120,结束,无布里渊谱存在。采用上述方案,能够对采集数据是否含有布里渊谱及其位置进行准确判断和定位,一方面根据判断结果只对含有布里渊谱的数据进行曲线拟合和谱分析,避免了很多无谓的布里渊谱分析过程,节省整体数据输出时间,提高了整体测试速度,另一方面提高了各种信噪比情况下布里渊谱分析的准确性,提高了应变测试精度。
应当理解的是,对本领域普通技术人员来说,可以根据上述说明加以改进或变换,而所有这些改进和变换都应属于本发明所附权利要求的保护范围。

Claims (7)

1.一种基于能量分布的布里渊谱识别方法,其特征在于,包括以下步骤:
步骤100,确定布里渊谱测试数据为L;
步骤101,根据探测脉冲宽度PW选择相应的能量积分窗口宽度W,并根据L及W计算,得到能量密度数据为P;
步骤102,在能量密度数据P中搜索预定带宽DW条件下能量密度最高的频率范围SW,并将此频率范围提供给步骤103;
步骤103,在能量密度数据P中,在所述能量密度最高的频率范围SW内,计算其峰值是否超过能量阈值PT,如超过能量阈值PT则进入步骤104,如不超过则进入步骤108;
步骤104,在能量密度数据P中,搜索所述最高的频率范围内峰值对应的频率点;
步骤105,在布里渊谱测试数据L中,以所述峰值为中心频率计算3dB带宽,并与3dB带宽阈值比较,如超过阈值则进入步骤109,如未超过则进入步骤106;
步骤106,在布里渊谱测试数据L中,以所述峰值为中心频率计算6dB带宽,并与6dB带宽阈值比较,如超过阈值则进入步骤109,如未超过则进入步骤107;
步骤107,在布里渊谱测试数据L中,以所述峰值为中心频率计算9dB带宽,并与9dB带宽阈值比较,如超过阈值则进入步骤109,如未超过则进入步骤120;
步骤108,未超过则认为此布里渊谱测试数据L中都为噪声,无布里渊谱存在;
步骤109,识别过程结束,输出中心频率及其在数据中的位置;
步骤120,结束,无布里渊谱存在。
2.如权利要求1所述的基于能量分布的布里渊谱识别方法,其特征在于,所述步骤101中,所述脉冲宽度PW越窄,能量积分窗口宽度越宽。
3.如权利要求1所述的基于能量分布的布里渊谱识别方法,其特征在于,所述步骤101中,所述布里渊谱测试数据L根据积分窗口宽度W进行积分运算。
4.如权利要求1所述的基于能量分布的布里渊谱识别方法,其特征在于,所述步骤102中,所述预定带宽DW的选择范围应小于能量积分窗口宽度W。
5.如权利要求1所述的基于能量分布的布里渊谱识别方法,其特征在于,所述步骤102中,所述预定带宽DW为能量积分窗口宽度W的1/2。
6.如权利要求1所述的基于能量分布的布里渊谱识别方法,其特征在于,所述步骤103中,所述能量阈值PT是根据量程最远的100个点的布里渊谱测试数据L进行平均计算获取。
7.如权利要求1所述的基于能量分布的布里渊谱识别方法,其特征在于,所述步骤105、步骤106及步骤107中,所述带宽阈值是根据布里渊谱测试数据L的探测脉冲宽度确定。
CN201310264845.7A 2013-06-27 2013-06-27 一种基于能量分布的布里渊谱识别方法 Active CN103308083B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310264845.7A CN103308083B (zh) 2013-06-27 2013-06-27 一种基于能量分布的布里渊谱识别方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310264845.7A CN103308083B (zh) 2013-06-27 2013-06-27 一种基于能量分布的布里渊谱识别方法

Publications (2)

Publication Number Publication Date
CN103308083A CN103308083A (zh) 2013-09-18
CN103308083B true CN103308083B (zh) 2015-12-09

Family

ID=49133566

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310264845.7A Active CN103308083B (zh) 2013-06-27 2013-06-27 一种基于能量分布的布里渊谱识别方法

Country Status (1)

Country Link
CN (1) CN103308083B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104296673A (zh) * 2014-10-22 2015-01-21 中国电子科技集团公司第四十一研究所 一种布里渊谱信号质量提升方法
CN104807568B (zh) * 2015-04-02 2017-11-03 南京大学 基于残缺谱拼接的布里渊谱寻峰方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0504834A2 (en) * 1991-03-20 1992-09-23 Nec Corporation Optical transmitter
JP3408789B2 (ja) * 2000-10-30 2003-05-19 日本電信電話株式会社 後方散乱光の測定方法およびその装置
JP4201995B2 (ja) * 2001-03-01 2008-12-24 日本電信電話株式会社 光ファイバひずみ計測方法およびその装置
CN103067087A (zh) * 2013-01-15 2013-04-24 杭州电子科技大学 一种布里渊增益谱的优化方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2012047699A (ja) * 2010-08-30 2012-03-08 Mie Univ 光ファイバ構造物変位計測装置及びその計測方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0504834A2 (en) * 1991-03-20 1992-09-23 Nec Corporation Optical transmitter
JP3408789B2 (ja) * 2000-10-30 2003-05-19 日本電信電話株式会社 後方散乱光の測定方法およびその装置
JP4201995B2 (ja) * 2001-03-01 2008-12-24 日本電信電話株式会社 光ファイバひずみ計測方法およびその装置
CN103067087A (zh) * 2013-01-15 2013-04-24 杭州电子科技大学 一种布里渊增益谱的优化方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Brillouin power spectrum analysis for partially uniformly strained optical fiber;Dan Zhang, etc;《Optics and Lasers in Engineering》;20090517;第976~981页 *
一种新的光纤布里渊传感散射谱拟合方法;肖尚辉,等;《光学技术》;20091130;第35卷(第6期);第897~900页 *

Also Published As

Publication number Publication date
CN103308083A (zh) 2013-09-18

Similar Documents

Publication Publication Date Title
US11808639B2 (en) High-precision temperature demodulation method oriented toward distributed fiber Raman sensor
CN106443203A (zh) 一种脉冲信号检测系统及方法
CN103344288A (zh) 一种基于零点分析的时差式超声波流量计测量方法
CN105572499A (zh) 一种基于矢量网络分析仪的眼图生成方法
CN105223481A (zh) 基于差值能量函数的局部放电特高频信号起始时刻确定方法
CN111413596A (zh) 电缆局部放电定位方法及装置、电子设备及定位系统
CN104237749A (zh) 超高频局放信号初始时刻判别方法
CN104777326A (zh) 具有流速自动监测实时校准功能的颗粒物计数仪
CN103308083B (zh) 一种基于能量分布的布里渊谱识别方法
CN108303396A (zh) 太赫兹时域脉冲信号采集方法、装置和系统
CN112363017A (zh) 一种基于小波变换的线路故障定位方法
CN115062071A (zh) 一种用于河道的水资源污染扩散分析方法及系统
CN103823175B (zh) 一种基于otdr的光电探测电路频响特性测试方法
CN114034327B (zh) 基于稀疏采样和人工神经网络的布里渊散射信号测量方法
CN102590598B (zh) 基于多阈值比较的超声波信号过零点预测方法
CN103490830B (zh) 基于物联网电力测温设备的去噪声射频频谱峰值获取方法
CN106646422B (zh) 增强相干测风雷达多普勒频移信号信噪比的预处理系统
CN104931785A (zh) 一种基于众数法的pcb特性阻抗的测试方法
CN105467270B (zh) 基于频谱相似性评价的单端行波测距反射波辨识算法
CN106771928A (zh) 一种局部放电脉冲初至时刻在线拾取方法
CN102944255B (zh) 数字化光纤光栅解调仪表及方法
CN104267338A (zh) 微波高速开关的开关时间测量系统及测量方法
CN110792564B (zh) 一种风机叶片引下线故障检测与定位方法及装置
RU2584266C1 (ru) Способ определения расстояния до места повреждения на линии электропередачи
CN103336063B (zh) 一种声发射信号初至点检测方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CP02 Change in the address of a patent holder
CP02 Change in the address of a patent holder

Address after: No. 726, Huaguang Avenue, Bengbu, Anhui Province, Anhui

Patentee after: The 41st Institute of CETC

Address before: 266000 Qingdao economic and Technological Development Zone, Shandong Xiangjiang Road, No. 98

Patentee before: The 41st Institute of CETC

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20190314

Address after: 266000 No. 98 Xiangjiang Road, Huangdao District, Qingdao City, Shandong Province

Patentee after: China Electronics Technology Instrument and Meter Co., Ltd.

Address before: 233000 726 Huaguang Road, Bengbu, Anhui

Patentee before: The 41st Institute of CETC

CP03 Change of name, title or address
CP03 Change of name, title or address

Address after: Huangdao Xiangjiang Road 266555 Shandong city of Qingdao Province, No. 98

Patentee after: CLP kesiyi Technology Co.,Ltd.

Address before: 266000 No. 98 Xiangjiang Road, Huangdao District, Shandong, Qingdao

Patentee before: CHINA ELECTRONIC TECHNOLOGY INSTRUMENTS Co.,Ltd.