CN104215271B - Positioning method for disturbance position in distributed optical fiber disturbance monitoring system - Google Patents

Positioning method for disturbance position in distributed optical fiber disturbance monitoring system Download PDF

Info

Publication number
CN104215271B
CN104215271B CN201410368811.7A CN201410368811A CN104215271B CN 104215271 B CN104215271 B CN 104215271B CN 201410368811 A CN201410368811 A CN 201410368811A CN 104215271 B CN104215271 B CN 104215271B
Authority
CN
China
Prior art keywords
optical fiber
disturbance
frequency spectrum
positioning
logarithm
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
CN201410368811.7A
Other languages
Chinese (zh)
Other versions
CN104215271A (en
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.)
Fudan University
Original Assignee
Fudan University
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 Fudan University filed Critical Fudan University
Priority to CN201410368811.7A priority Critical patent/CN104215271B/en
Publication of CN104215271A publication Critical patent/CN104215271A/en
Application granted granted Critical
Publication of CN104215271B publication Critical patent/CN104215271B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

The invention belongs to the technical field of distributed measurement, and particularly relates to a positioning method for a disturbance position in a distributed optical fiber disturbance monitoring system. The method includes firstly, performing Fourier transform on an acquired interferometric phase signal to acquire a frequency spectrum phi (F), taking the logarithm of the frequency spectrum to acquire a logarithm frequency spectrum, subjecting the logarithm frequency spectrum of the signal to filter subtraction or/and autocorrelation, and finally performing Fourier transform to acquire the disturbance position. The positioning method has the advantages of obvious peak value, centralized peak energy, accurate positioning and wide available positioning range. The positioning method is particularly applicable to distributed optical fiber monitoring and positioning systems.

Description

The localization method of disturbance location in a kind of distributed optical fiber disturbance monitoring system
Technical field
The invention belongs to distributed measurement technical field, specifically, it is related to a kind of distributed optical fiber disturbance monitoring system The localization method of middle disturbance location.
Background technology
Distributed optical fiber disturbance monitoring and positioning system is widely used in measurement, simultaneously measurement signal be also used for right The specificity analysis of monitoring objective.For example, in oil pipeline, replace manually patrolling using distributed optical fiber disturbance monitoring and positioning system Patrol, find to steal oil, oil leak, the behavior such as destruction pipeline in time, and be accurately positioned the position of behavior generation, intelligent management pipeline, Guarantee pipeline safety.In safety-security area, due to optical fiber disguise with anti-electromagnetically, distributed optical fiber disturbance monitoring positioning system System is very suitable for all kinds of accident of real-time monitoring, and it is high to monitor the more common electro-optical system of success rate.
Location algorithm in traditional distributed optical fiber disturbance monitoring and positioning system mainly has trap wave point method and time delay Method.Trap wave point method chooses the single frequency depression points of interferometric phase signal spectrum or several frequency depression points, by these frequencies Relation between the position of point and disturbance point position is positioned.This method mainly has two shortcomings.One is these Frequency points Position easily affected to shift by other factors, lead to positioning precision not high;Two is that multiple peak-seeking algorithm effect is bad, and Ratio is relatively time-consuming, usually chooses the depression frequency of mistake, leads to finally position failure.Time-delay method is using the two paths of signals of construction Between delay inequality positioned with the relation of disturbance point distance, this time delay is extracted by Appropriate method, and utilizes corresponding relation Change into orientation distance.The major defect of the method is to adopt iterative algorithm, leads to calculate slowly, processing equipment is required relatively Height, and light channel structure is complicated.
Content of the invention
The present invention is directed to distributed optical fiber disturbance detecting system it is proposed that disturbing in a kind of distributed optical fiber disturbance monitoring system The localization method of dynamic position;Its setting accuracy is high.
The present invention is directed to distributed optical fiber disturbance detecting system shown in Fig. 5 it is assumed that there is disturbance, interference signal at 5 in Fig. 5 Amplitude conditioning through amplifying or decay, obtains phase signal (t) by phase place retrieving algorithm, and through a FFT Obtain the frequency spectrum of phase signal, form can be expressed as:
Wherein orientation distance information is included inXiang Zhong, byDetermine, the relation with orientation distance is,
According to theory, directly frequency spectrum is carried out a Fourier transform again, obtain, you can obtain disturbance location.But It is due to the frequency spectrum of disturbance itselfAnd other factors, lead to containing time delayFeature inconspicuous, impact positioning Effect.It is necessary to nurse one's health in second Fourier transform advance traveling wave shape.
Carry out disturbance location positioning waveform conditioning using two kinds of technology, the first is filtering subtraction technique, that is, in the present invention By log spectrum through LPF, obtain the envelope of log spectrum, deduct this envelope with former frequency spectrum, make log spectrum " drawing high ", near 0 value, reduces again the low frequency background of FFT(The i.e. strong background in nearly orientation distance region)So that it is near Distance positioning accuracy greatly improves.Second is autocorrelation technique, and the frequency spectrum after will processing carries out auto-correlation so that frequency spectrum Modulation trend become apparent from, beneficial to the effective energy concentration degree carrying out during FFT again, thus obtaining in full orientation range Obtain and more accurately position.
In the present invention, filtering subtraction technique and autocorrelation technique can be used alone or be superimposed use.Concrete scheme is situated between Continue as follows.
The present invention provides a kind of localization method of disturbance location in distributed optical fiber disturbance monitoring system, and it is individually using filter Ripple subtraction technique, comprises the following steps that:
(1)First two-way interference signal PIN I and PIN II on the optical fiber collecting is obtained interfering phase through phase place reduction Potentiometer signal, then it is carried out Fourier transform and obtains frequency spectrum, then to frequency spectrumTake the logarithm, obtain logarithm Frequency spectrum
(2)Using filtering subtraction technique and/or autocorrelation technique, carry out disturbance location positioning waveform conditioning;Wherein, described Filtering subtraction technique is by log spectrumBy low pass filter, obtain log spectrum envelope;Again will Log spectrumDeduct log spectrum envelopeObtain;Described autocorrelation technique is by logarithm frequency SpectrumOrCarry out autocorrelation calculation;
(3)Finally, to step(2)Result, carry out Fourier transformation and obtain orientation distance.
The beneficial effects of the present invention is:
(1)The present invention is the place of the disturbance location information of the distributed optical fiber disturbance positioning system with structure shown in Fig. 5 Reason provides an effective method.
(2)The present invention can extract out the positional information of disturbing source, and can get the spectral characteristic of disturbing source, is to disturbing source Specificity analysis, judge provide solution.
(3)Filtering subtracted signal in pilot process of the present invention, autocorrelation signal can be used for judging perturbation features so that Positioning result is accurate.
Brief description
Fig. 1 is this invention difference corresponding design sketch of implementation steps.
Fig. 2 is this bright position error figure at 100km in an embodiment.
Fig. 3 is a series of auditory localization cues figures in the embodiment of the present invention.
Fig. 4 is the algorithm flow chart of one embodiment of the present invention.
Fig. 5 is the light path that embodiment is adopted.
In figure label:1 is 3*3 optical fiber splitter;2 is 2*2 optical fiber splitter;3 is time delay optical fiber, and one end is divided with optical fiber The 1b1 port of road device 1 is connected, and the other end is connected with the 2b1 port of optical fiber splitter 2;4 is sensing optic cable, from optical fiber splitter 2 2a port draw, end is connected with 6;5 is disturbance point active position, amount as to be detected;6 is perflectometer.
Specific embodiment
In the present embodiment, optical fiber distributed perturbation alignment system light path concrete structure is as shown in Figure 5.Centre wavelength exists The wide spectrum light source of 1550nm and the port 1a1 input of optical fiber splitter 1, are divided into two bundles, respectively from another after optical fiber splitter 1 Port 1b1,1b2 output of side.Optical interference circuit is accordingly:
A1:1b1→3→2b1→2a→4→5→6→5→4→2a→2b2→1b2
A2:1b2→2b2→2a→4→5→6→5→4→2a→2b1→3→1b1
In the present embodiment, the light source of 1550nm is all the super radiation light emitting tube of research institute of parent company of electronics group 44 production (SLD).Optical fiber splitter, wavelength division multiplexer are that Wuhan post and telecommunications research institute produces.1 is the fiber coupler that 3*3 divides equally, and 2 is 2*2 The fiber coupler divided equally.Feedback device 6 is made up of optical fiber connector plating reflectance coating.Monitoring optical cable(Optical fiber 4)It is layed in needs The vicinity of monitoring pipeline, fiber optic interferometric module need to be placed in noise insulation equipment to shield external interference.From optical fiber splitter 1a2,1a3 The interference signal of output is switched to electric signal with the InGaAs photodetector of 44 models GT322C500 being produced.? In the range of 10km to 130km, every 10km produces 20 disturbances.By National Instruments company data capture card Electrical signal is gathered and carries out signal transacting into computer by PCI-6122, data processing software Labview software programming, low Bandpass filter is Butterworth filter, is worth for 0.75, exponent number 2 rank.Observe the disturbance positioning distribution of each perturbation distance.
Fig. 1 is this invention difference corresponding design sketch of implementation steps, lists no corrective measure respectively, only filtering is subtracted each other, Auto-correlation, the positioning result schematic diagram of four kinds of implementations of auto-correlation is subtracted each other in filtering, it can be found that superposition applies filtering subtracting each other Best with autocorrelative effect.
Fig. 2 is this bright position error figure at 100km in an embodiment.Positioning result is distributed in adjacent four and takes Value point, each data point represent interval be 25m about it can be seen that locating effect is very good.
Fig. 3 is a series of auditory localization cues figures in the embodiment of the present invention, covers the orientation range of 10km ~ 130km, and peak value is bright Aobvious, in peak energy quantity set.

Claims (1)

1. in a kind of distributed optical fiber disturbance monitoring system the localization method of disturbance location it is characterised in that comprising the following steps that:
(1)First two-way interference signal PIN I and PIN II on the optical fiber collecting is obtained interference phase difference through phase place reduction Signal, then it is carried out Fourier transform and obtains frequency spectrum, then to frequency spectrumTake the logarithm, obtain log spectrum
(2)Using filtering subtraction technique and/or autocorrelation technique, carry out disturbance location positioning waveform conditioning;Wherein, described filtering Subtraction technique is by log spectrumBy low pass filter, obtain log spectrum envelope;Again by logarithm frequency SpectrumDeduct log spectrum envelopeObtain;Described autocorrelation technique is by log spectrumOrCarry out autocorrelation calculation;
(3)Finally, to step(2)Result, carry out Fourier transformation and obtain orientation distance.
CN201410368811.7A 2014-07-30 2014-07-30 Positioning method for disturbance position in distributed optical fiber disturbance monitoring system Active CN104215271B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410368811.7A CN104215271B (en) 2014-07-30 2014-07-30 Positioning method for disturbance position in distributed optical fiber disturbance monitoring system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410368811.7A CN104215271B (en) 2014-07-30 2014-07-30 Positioning method for disturbance position in distributed optical fiber disturbance monitoring system

Publications (2)

Publication Number Publication Date
CN104215271A CN104215271A (en) 2014-12-17
CN104215271B true CN104215271B (en) 2017-02-15

Family

ID=52097004

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410368811.7A Active CN104215271B (en) 2014-07-30 2014-07-30 Positioning method for disturbance position in distributed optical fiber disturbance monitoring system

Country Status (1)

Country Link
CN (1) CN104215271B (en)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104568122B (en) * 2015-01-21 2017-10-31 东南大学 A kind of distributed optical fiber vibration sensing system disturbance event identification and localization method
CN105277272A (en) * 2015-10-25 2016-01-27 复旦大学 Distributed optical fiber vibration sensing multi-point disturbance localization algorithm
CN107063584B (en) * 2017-04-22 2020-05-05 中南大学 Boiler tube leakage judging and positioning method
CN107270952B (en) * 2017-07-27 2020-03-31 天津求实飞博科技有限公司 Long-distance optical fiber distributed disturbance sensing signal processing method based on optical frequency domain reflection
CN111337116A (en) * 2018-12-18 2020-06-26 中兴通讯股份有限公司 Method and device for determining position of disturbance source, storage medium and electronic device
CN109668620B (en) * 2019-02-27 2021-12-14 深圳大学 Positioning method based on interference type distributed optical fiber vibration sensor
CN110057387B (en) * 2019-05-10 2021-01-19 南昌航空大学 Positioning method of linear Sagnac distributed optical fiber sensing system
CN110146116B (en) * 2019-06-19 2021-04-27 南昌航空大学 Positioning method for Sagnac optical fiber sensing under multipoint disturbance
CN112161778B (en) * 2020-08-17 2022-08-02 南昌航空大学 Distributed optical fiber disturbance positioning method based on regression probability distribution

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101393268A (en) * 2008-11-06 2009-03-25 复旦大学 Earthquake early-warning system based on communication optical cable
CN101393269A (en) * 2008-11-06 2009-03-25 复旦大学 Method for monitoring geology by utilizing communication optical cable
CN102360519B (en) * 2011-07-29 2013-04-03 电子科技大学 Intrusion detecting and positioning method for distributed optical fiber fence
US9417103B2 (en) * 2011-09-20 2016-08-16 Schlumberger Technology Corporation Multiple spectrum channel, multiple sensor fiber optic monitoring system
CN102840909B (en) * 2012-08-21 2014-04-30 天津大学 Optical frequency domain reflection-distributed vibration frequency sensing and locating device and demodulation method
CN102879081B (en) * 2012-09-17 2016-06-29 北京航天时代光电科技有限公司 A kind of data processing method in distributed optical fiber vibration system
CN103090961B (en) * 2013-02-01 2014-11-26 华中科技大学 Disturbance source positioning method of distributed type optical fiber sensing system
CN103499383B (en) * 2013-09-06 2016-01-13 北京航天时代光电科技有限公司 A kind of auto-correlation localization method improving optical fibre vibration sensor positioning precision
CN103743461A (en) * 2013-12-24 2014-04-23 南京邮电大学 Distributed optical fiber vibration sensor vibration-source positioning method

Also Published As

Publication number Publication date
CN104215271A (en) 2014-12-17

Similar Documents

Publication Publication Date Title
CN104215271B (en) Positioning method for disturbance position in distributed optical fiber disturbance monitoring system
CN103900623B (en) Optical time domain reflectometer and its common mode inhibition method based on alliteration optical modulator
CN103090961B (en) Disturbance source positioning method of distributed type optical fiber sensing system
CN107270952B (en) Long-distance optical fiber distributed disturbance sensing signal processing method based on optical frequency domain reflection
CN104077267A (en) Location method for double Fourier transform and disturbance
CN111222743B (en) Method for judging vertical offset distance and threat level of optical fiber sensing event
CN102566533B (en) On-line monitoring device and method for preparing tobacco essence perfume
CN112525201B (en) Underwater target tracking method based on electromagnetic field characteristic multi-information fusion
CN102920464A (en) Instrument for measuring haemoglobin concentration and blood oxygen saturation and measuring method
CN203519025U (en) High-speed fiber grating demodulation system
CN103412137B (en) With speed-measuring method and device in twiddle factor
CN106353003A (en) Distributed optical fiber temperature measuring method and system
CN111413596A (en) Cable partial discharge positioning method and device, electronic equipment and positioning system
CN102564476B (en) Multipoint disturbance positioning method
CN105181152A (en) Calculation method for frequency shift of distributed Brillouin scattered spectrum
CN106323878A (en) Laser remote sensing detector for natural gas concentration and distance
CN111811549A (en) Demodulation method and positioning device of asymmetric double Mach-Zehnder sensing system
CN106452570A (en) Optical fiber fault detection system and detection method based on optical fiber phase mediation principle
CN102657521B (en) Full optical fiber heart rate measuring equipment
CN103823175B (en) A kind of photodetection circuit frequency response function test method based on OTDR
CN109696409A (en) A kind of spectral scan probe
CN111157102B (en) Positioning method for eliminating frequency interference in distributed optical fiber sensing system
CN101726405A (en) Frequency response parameter test system for front-mounted amplifier circuit of photoelectric detector and test method thereof
CN109798448A (en) Concrete duct leakage experiment device and method based on anti-Stokes light filtering
CN105319479B (en) Two ends of electric transmission line fault localization system

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