CN104866708B - 一种botdr光纤应变分布曲线拼接方法 - Google Patents
一种botdr光纤应变分布曲线拼接方法 Download PDFInfo
- Publication number
- CN104866708B CN104866708B CN201510178879.3A CN201510178879A CN104866708B CN 104866708 B CN104866708 B CN 104866708B CN 201510178879 A CN201510178879 A CN 201510178879A CN 104866708 B CN104866708 B CN 104866708B
- Authority
- CN
- China
- Prior art keywords
- data
- strain
- distribution curve
- botdr
- sddc
- 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
Links
- 239000000835 fiber Substances 0.000 title claims abstract description 22
- 238000000034 method Methods 0.000 title claims abstract description 20
- VMSRVIHUFHQIAL-UHFFFAOYSA-N sodium;dimethylcarbamodithioic acid Chemical compound [Na+].CN(C)C(S)=S VMSRVIHUFHQIAL-UHFFFAOYSA-N 0.000 claims description 27
- 239000013307 optical fiber Substances 0.000 abstract description 20
- 238000012360 testing method Methods 0.000 abstract description 13
- 238000000205 computational method Methods 0.000 abstract description 2
- 230000010354 integration Effects 0.000 abstract description 2
- 238000012544 monitoring process Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 239000004744 fabric Substances 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000010276 construction Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000005303 weighing Methods 0.000 description 1
Landscapes
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明公开了一种BOTDR光纤应变分布曲线拼接方法，该方法依据光纤两端测试得到的应变数据的重合相似特性，以应变数据关联权重计算方法为核心，利用滑动求商权重积分寻找两条光纤应变曲线的重合区域，实现了光纤应变曲线的自动拼接，有效的减少了光纤末端信噪比较低的应变数据对曲线拼接的影响，极大的提高了光纤应变分布数据的拼接效率以及拼接精度，降低了人工成本，提升了BOTDR产品的测试距离，实现了对超长距离光纤应变分布的测试精度，扩大了BOTDR产品的应用领域。
Description
技术领域
本发明涉及一种BOTDR光纤应变分布曲线拼接方法。
背景技术
布里渊光时域反射计(BOTDR)依靠测量光纤中后向布里渊散射光的布里渊频移分布计算光纤的应变分布，脉冲光以一定的频率自光纤一端入射，入射的脉冲光与光纤中的声子发生相互作用后产生布里渊散射，后向布里渊散射光沿光纤原路返回到入射端。由于光纤中布里渊散射光频移与光纤轴向应变和温度间存在线性关系，因此测量出光纤的布里渊散射频移分布即可计算出光纤中的应变分布。布里渊光时域反射计具有低能源依赖性、高环境耐受性、抗电磁干扰、抗腐蚀、防水、抗潮湿、温度适应性强等特性，并且由于可以单端测量，施工难度较低而广受关注。布里渊光时域反射计可用于岩土工程健康监测、地质灾害预警监测、电缆及管道的健康监测等领域，是工程领域用于取代传统点式传感器的最有力的产品之一。
BOTDR产品在测试光纤应变分布数据时，空间分辨率指标直接受探测光脉冲宽度影响，脉冲宽度越宽，空间分辨率越差，测试距离越长，在测量长距离光纤应变分布数据时，为保证布里渊后向散射信号的强度，只能使用大脉宽进行测量，测试数据的空间分辨率很差且难以提升。同时，现有BOTDR产品的最大测试量程为80km，难以满足距离超过80km的超长距离光纤的应变分布测试需求。根据BOTDR产品的测试原理，在光纤一端测试的应变曲线与在光纤另一端测试得到的应变曲线互为镜像，因此将光纤两端各自测试的应变曲线进行互补拼接，可以有效的改善单端测试光纤应变曲线的末端应变数据信噪比，既可以提升BOTDR产品的测试距离，也可以使用更窄的脉冲实现相同长度的光纤应变分布测试，从而提升BOTDR产品的空间分辨率。
但是，目前对光纤两端的BOTDR应变曲线进行拼接只能采用人工手段，通过人工寻找两端应变曲线的重合区域，然后手动对两组应变曲线进行拼接，由于测试距离长，两条测量应变曲线的末端数据信噪比很低，噪声很高，人工拼接耗费时间长，成本高，而且误差较大，极大的影响了BOTDR产品对超长距离光纤应变分布测试能力，以及在超长距离光纤应变分布测试领域的推广及应用。
发明内容
针对现有技术中存在的上述技术问题，本发明提出了一种BOTDR光纤应变分布曲线拼接方法，其采用如下技术方案：
一种BOTDR光纤应变分布曲线拼接方法，包括如下步骤：
步骤101、获取应变分布曲线的数据点数N；
步骤102、获取应变分布曲线A的数据SDA[0]～SDA[N-1]；
步骤103、获取应变分布曲线B的数据SDB[0]～SDB[N-1]；
步骤104、计算应变分布曲线B的镜像翻转数据SDBI[0]～SDBI[N-1]；
步骤105、计算应变数据关联权重数据SDC[0]～SDC[N-1]；
步骤106、寻找应变数据关联权重数据SDC[0]～SDC[N-1]的最大值所在位置SDMAXI；
步骤107、建立应变分布拼接数据SDDC[0]～SDDC[2N-1]，并全部初始化为0；
步骤108、将应变分布曲线A的数据SDA[0]～SDA[N-SDMAXI/2-1]依次赋值给SDDC[0]～SDDC[N-SDMAXI/2-1]；
步骤109、将应变分布曲线B的数据SDB[SDMAXI/2]～SDB[N-1]依次赋值给SDDC[N-SDMAXI/2]～SDDC[2N-SDMAXI/2-1]；
步骤110、输出数据SDDC[0]～SDDC[2N-1-SDMAXI/2]。
进一步，上述步骤104中，计算应变分布曲线B的镜像翻转数据的方法为：
步骤1041、建立应变分布曲线B的水平镜像数据SDBI[0]～SDBI[N-1]，全部初始化为0；
步骤1042、初始化I，赋值为0；
步骤1043、将SDBI[I]赋值为SDB[N-1-I]；
步骤1044、I的值增加1；
步骤1045、如果I<N，转到步骤1043，否则转步骤1046；
步骤1046、结束。
进一步，上述步骤105中，计算应变数据关联权重数据的方法为：
步骤1051、初始化应变数据关联权重数据SDC[0]～SDC[N-1]的值为0；
步骤1052、初始化T为0；
步骤1053、初始化S为0；
步骤1054、初始化I为T；
步骤1055、将X赋值为SDA[I]/SDB[I-T]；
步骤1056、计算关联权重CE，CE＝A×exp[-0.5*[(X-X0)/B]2]，其中，A、X0与B均为高斯函数系数；
步骤1057、将S的值增加CE；
步骤1058、如果I<N-T，转步骤1055，否则转步骤1059；
步骤1059、将SDC[T]赋值为S；
步骤10510、如果T<N-1，转步骤1053，否则转步骤10511；
步骤10511、结束。
本发明具有如下优点：
本发明方法依据光纤两端测试得到的应变数据的重合相似特性，以应变数据关联权重计算方法为核心，利用滑动求商权重积分寻找两条光纤应变曲线的重合区域，实现了光纤应变曲线的自动拼接，有效的减少了光纤末端信噪比较低的应变数据对曲线拼接的影响，极大的提高了光纤应变分布数据的拼接效率以及拼接精度，降低了人工成本，提升了BOTDR产品的测试距离，实现了对超长距离光纤应变分布的测试精度，扩大了BOTDR产品的应用领域。
附图说明
图1为本发明中一种BOTDR光纤应变分布曲线拼接方法的流程示意图；
图2为本发明中应变分布曲线B的镜像翻转数据计算方法示意图；
图3为本发明中应变数据关联权重计算方法示意图。
具体实施方式
下面结合附图以及具体实施方式对本发明作进一步详细说明：
结合图1所示，一种BOTDR光纤应变分布曲线拼接方法，包括如下步骤：
步骤101、获取应变分布曲线的数据点数N。
步骤102、获取应变分布曲线A的数据SDA[0]～SDA[N-1]。
步骤103、获取应变分布曲线B的数据SDB[0]～SDB[N-1]。
步骤104、计算应变分布曲线B的镜像翻转数据SDBI[0]～SDBI[N-1]。
如图2所示，应变分布曲线B的镜像翻转数据计算方法如下：
步骤1041、建立应变分布曲线B的水平镜像数据SDBI[0]～SDBI[N-1]，全部初始化为0；
步骤1042、初始化I，赋值为0；
步骤1043、将SDBI[I]赋值为SDB[N-1-I]；
步骤1044、I的值增加1；
步骤1045、如果I<N，转到步骤1043，否则转步骤1046；
步骤1046、结束。
步骤105、计算应变数据关联权重数据SDC[0]～SDC[N-1]。
如图3所示，计算应变数据关联权重数据的方法为：
步骤1051、初始化应变数据关联权重数据SDC[0]～SDC[N-1]的值为0；
步骤1052、初始化T为0；
步骤1053、初始化S为0；
步骤1054、初始化I为T；
步骤1055、将X赋值为SDA[I]/SDB[I-T]；
步骤1056、计算关联权重CE，CE＝A×exp[-0.5*[(X-X0)/B]2]；
其中，A、X0与B均为高斯函数系数，其典型值可分别设为1、1、0.5，A、X0与B的值不限于典型值，可根据实际情况确定；
步骤1057、将S的值增加CE；
步骤1058、如果I<N-T，转步骤1055，否则转步骤1059；
步骤1059、将SDC[T]赋值为S；
步骤10510、如果T<N-1，转步骤1053，否则转步骤10511；
步骤10511、结束。
步骤106、寻找应变数据关联权重数据SDC[0]～SDC[N-1]的最大值所在位置SDMAXI。
步骤107、建立应变分布拼接数据SDDC[0]～SDDC[2N-1]，并全部初始化为0。
步骤108、将应变分布曲线A的数据SDA[0]～SDA[N-SDMAXI/2-1]依次赋值给SDDC[0]～SDDC[N-SDMAXI/2-1]。
步骤109、将应变分布曲线B的数据SDB[SDMAXI/2]～SDB[N-1]依次赋值给SDDC[N-SDMAXI/2]～SDDC[2N-SDMAXI/2-1]。
步骤110、输出数据SDDC[0]～SDDC[2N-1-SDMAXI/2]。
当然，以上说明仅仅为本发明的较佳实施例，本发明并不限于列举上述实施例，应当说明的是，任何熟悉本领域的技术人员在本说明书的教导下，所做出的所有等同替代、明显变形形式，均落在本说明书的实质范围之内，理应受到本发明的保护。
Claims (2)
1.一种BOTDR光纤应变分布曲线拼接方法，其特征在于，包括如下步骤：
步骤101、获取应变分布曲线的数据点数N；
步骤102、获取应变分布曲线A的数据SDA[0]～SDA[N-1]；
步骤103、获取应变分布曲线B的数据SDB[0]～SDB[N-1]；
步骤104、计算应变分布曲线B的镜像翻转数据SDBI[0]～SDBI[N-1]；
步骤105、计算应变数据关联权重数据SDC[0]～SDC[N-1]；
步骤106、寻找应变数据关联权重数据SDC[0]～SDC[N-1]的最大值所在位置SDMAXI；
步骤107、建立应变分布拼接数据SDDC[0]～SDDC[2N-1]，并全部初始化为0；
步骤108、将应变分布曲线A的数据SDA[0]～SDA[N-SDMAXI/2-1]依次赋值给SDDC[0]～SDDC[N-SDMAXI/2-1]；
步骤109、将应变分布曲线B的数据SDB[SDMAXI/2]～SDB[N-1]依次赋值给SDDC[N-SDMAXI/2]～SDDC[2N-SDMAXI/2-1]；
步骤110、输出数据SDDC[0]～SDDC[2N-1-SDMAXI/2]；
所述步骤105中，计算应变数据关联权重数据的方法为：
步骤1051、初始化应变数据关联权重数据SDC[0]～SDC[N-1]的值为0；
步骤1052、初始化T为0；
步骤1053、初始化S为0；
步骤1054、初始化I为T；
步骤1055、将X赋值为SDA[I]/SDB[I-T]；
步骤1056、计算关联权重CE，CE＝A×exp[-0.5*[(X-X0)/B]2]，其中，A、X0与B均为高斯函数系数；
步骤1057、将S的值增加CE；
步骤1058、如果I<N-T，转步骤1055，否则转步骤1059；
步骤1059、将SDC[T]赋值为S；
步骤10510、如果T<N-1，转步骤1053，否则转步骤10511；
步骤10511、结束。
2.根据权利要求1所述的一种BOTDR光纤应变分布曲线拼接方法，其特征在于，所述步骤104中，计算应变分布曲线B的镜像翻转数据的方法为：
步骤1041、建立应变分布曲线B的水平镜像数据SDBI[0]～SDBI[N-1]，全部初始化为0；
步骤1042、初始化I，赋值为0；
步骤1043、将SDBI[I]赋值为SDB[N-1-I]；
步骤1044、I的值增加1；
步骤1045、如果I<N，转到步骤1043，否则转步骤1046；
步骤1046、结束。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510178879.3A CN104866708B (zh) | 2015-04-16 | 2015-04-16 | 一种botdr光纤应变分布曲线拼接方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510178879.3A CN104866708B (zh) | 2015-04-16 | 2015-04-16 | 一种botdr光纤应变分布曲线拼接方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104866708A CN104866708A (zh) | 2015-08-26 |
| CN104866708B true CN104866708B (zh) | 2017-08-29 |
Family
ID=53912532
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510178879.3A Active CN104866708B (zh) | 2015-04-16 | 2015-04-16 | 一种botdr光纤应变分布曲线拼接方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104866708B (zh) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110635840B (zh) * | 2018-12-19 | 2022-09-27 | 国网湖北省电力有限公司超高压公司 | 一种双向otdr测试曲线的拼接方法 |
| CN111490818B (zh) * | 2020-04-10 | 2022-03-11 | 武汉光迅科技股份有限公司 | 一种长距离光纤检测方法及装置、设备、系统、存储介质 |
| CN111609996B (zh) * | 2020-05-18 | 2022-08-02 | 昂纳信息技术（深圳）有限公司 | 一种otdr曲线拼接方法及多条otdr曲线拼接方法 |
| CN115001571B (zh) * | 2022-05-11 | 2023-09-12 | 青岛诺克通信技术有限公司 | 一种光时域反射计的自适应数字滤波方法及光时域反射计 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101226078A (zh) * | 2008-01-30 | 2008-07-23 | 广厦建设集团有限责任公司 | 一种基于分布式光纤传感器的长距离线性结构异常振动的检测方法 |
| CN102539113A (zh) * | 2011-11-21 | 2012-07-04 | 河南省电力通信自动化公司 | 一种基于布里渊光时域分析仪botda的曲线分析方法 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8922875B2 (en) * | 2011-10-03 | 2014-12-30 | Daisy Williams | Method for optical phase modulation, optical phase modulator and optical network using the same |
| US8982340B2 (en) * | 2012-07-20 | 2015-03-17 | Ramot At Tel-Aviv University Ltd. | Fast brillouin optical time domain analysis for dynamic sensing |
-
2015
- 2015-04-16 CN CN201510178879.3A patent/CN104866708B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101226078A (zh) * | 2008-01-30 | 2008-07-23 | 广厦建设集团有限责任公司 | 一种基于分布式光纤传感器的长距离线性结构异常振动的检测方法 |
| CN102539113A (zh) * | 2011-11-21 | 2012-07-04 | 河南省电力通信自动化公司 | 一种基于布里渊光时域分析仪botda的曲线分析方法 |
Non-Patent Citations (2)
| Title |
|---|
| 一种提高BOTDR系统信噪比的方法研究;李勇涛 等;《光通信研究》;20101031;第2010年卷(第5期);图1、图2、第0节倒数第一段、第2节第1段、第3.1节倒数第1段 * |
| 基于BOTDR的光纤复合海底电缆应变/温度监测;吕安强 等;《高电压技术》;20140228;第40卷(第2期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104866708A (zh) | 2015-08-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108519175B (zh) | 基于布拉格光纤光栅的可变量程的土体压力测量方法 | |
| CN108760109B (zh) | 基于布拉格光纤光栅的可变量程的土体压力测量装置和方法 | |
| Xu et al. | Full scale strain monitoring of a suspension bridge using high performance distributed fiber optic sensors | |
| CN104866708B (zh) | 一种botdr光纤应变分布曲线拼接方法 | |
| Zhao et al. | On-line monitoring system of 110 kV submarine cable based on BOTDR | |
| Zhang et al. | Performance evaluation of BOTDR-based distributed fiber optic sensors for crack monitoring | |
| Madjdabadi et al. | Experimental evaluation of a distributed Brillouin sensing system for measuring extensional and shear deformation in rock | |
| Kuang | Distributed damage detection of offshore steel structures using plastic optical fibre sensors | |
| Zheng et al. | Theoretical and experimental study on fiber-optic displacement sensor with bowknot bending modulation | |
| CN105371785A (zh) | 一种曲率测量方法 | |
| Li et al. | Design of an enhanced sensitivity FBG strain sensor and application in highway bridge engineering | |
| Zhou et al. | Fiber-reinforced polymer-packaged optical fiber sensors based on Brillouin optical time-domain analysis | |
| Marković et al. | Application of fiber-optic curvature sensor in deformation measurement process | |
| Bernini et al. | Identification of defects and strain error estimation for bending steel beams using time domain Brillouin distributed optical fiber sensors | |
| JP2005147900A (ja) | Ｏｆｄｒ方式の歪連続分布計測装置 | |
| CN110375780B (zh) | Ofdr断纤续接测量的方法 | |
| Sah et al. | Temperature independent FBG based displacement sensor for crack detection in civil structures | |
| Zhang et al. | Spatial resolution of DOFS and its calibration methods | |
| Nawrot et al. | Mechanical strain-amplifying transducer for fiber Bragg grating sensors with applications in structural health monitoring | |
| Li et al. | Distributed optical fiber bi-directional strain sensor for gas trunk pipelines | |
| Wang et al. | Research and application of an improved internal thrust force measurement system for rock and soil mass based on OFDR | |
| Zeni et al. | Distributed fibre optic sensing techniques for soil slope monitoring | |
| Cheng et al. | Quasidistributed optical fiber sensing system to monitor the displacement of joints on the face slab of concrete-faced rockfill dams | |
| Bernini et al. | Damage detection in bending beams through Brillouin distributed optic-fibre sensor | |
| CN110987946A (zh) | 一种用分光器改进裂缝传感器精度的装置及方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| EXSB | Decision made by sipo to initiate substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| CB03 | Change of inventor or designer information |
Inventor after: Yuan Ming Inventor after: Li Ligong Inventor after: Yan Jisong Inventor after: Zhang Zhihui Inventor after: Chen Xiaohui Inventor before: Yuan Ming Inventor before: Yan Jisong Inventor before: Zhang Zhihui Inventor before: Chen Xiaofeng |
|
| COR | Change of bibliographic data | ||
| CB03 | Change of inventor or designer information | ||
| CB03 | Change of inventor or designer information |
Inventor after: Yuan Ming Inventor after: Li Ligong Inventor after: Yan Jisong Inventor after: Zhang Zhihui Inventor after: Chen Xiaofeng Inventor before: Yuan Ming Inventor before: Li Ligong Inventor before: Yan Jisong Inventor before: Zhang Zhihui Inventor before: Chen Xiaohui |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20190226 Address after: 266555 Xiangjiang 98, Huangdao District, Qingdao City, Shandong Province Patentee after: China Electronics Technology Instrument and Meter Co., Ltd. Address before: 266555 No. 98 Xiangjiang Road, Qingdao economic and Technological Development Zone, Shandong Patentee before: The 41st Institute of CETC |
|
| TR01 | Transfer of patent right | ||
| CP01 | Change in the name or title of a patent holder |
Address after: 266555 Xiangjiang 98, Huangdao District, Qingdao City, Shandong Province Patentee after: CLP kesiyi Technology Co.,Ltd. Address before: 266555 Xiangjiang 98, Huangdao District, Qingdao City, Shandong Province Patentee before: CHINA ELECTRONIC TECHNOLOGY INSTRUMENTS Co.,Ltd. |
|
| CP01 | Change in the name or title of a patent holder |