CN105973451A - Optical fiber vibration model determination method and device - Google Patents

Optical fiber vibration model determination method and device Download PDF

Info

Publication number
CN105973451A
CN105973451A CN201610301609.1A CN201610301609A CN105973451A CN 105973451 A CN105973451 A CN 105973451A CN 201610301609 A CN201610301609 A CN 201610301609A CN 105973451 A CN105973451 A CN 105973451A
Authority
CN
China
Prior art keywords
vibration
signal
spectral characteristic
fiber
model
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.)
Pending
Application number
CN201610301609.1A
Other languages
Chinese (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.)
Shenzhen Ai Rui Stone Technology Co Ltd
Original Assignee
Shenzhen Ai Rui Stone Technology Co Ltd
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 Shenzhen Ai Rui Stone Technology Co Ltd filed Critical Shenzhen Ai Rui Stone Technology Co Ltd
Priority to CN201610301609.1A priority Critical patent/CN105973451A/en
Publication of CN105973451A publication Critical patent/CN105973451A/en
Pending legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
    • G01H9/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves by using radiation-sensitive means, e.g. optical means
    • G01H9/004—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves by using radiation-sensitive means, e.g. optical means using fibre optic sensors

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

The invention discloses an optical fiber vibration model determination method and device. The method includes obtaining a vibration signal of a target optical fiber position; extracting frequency spectrum features of the vibration signal; matching the frequency spectrum features with the frequency spectrum features of a preset vibration model; determining that the vibration signal belongs to the preset vibration model if the frequency spectrum features are matched with the frequency spectrum features of the preset vibration model. By adopting the above method, the vibration type of the target optical fiber position can be judged accurately and the optical fiber vibration detection efficiency is improved.

Description

A kind of fiber-optic vibration model determines method and device
Technical field
The present invention relates to fiber-optic vibration detection technique field, particularly relate to a kind of fiber-optic vibration model Determine method and device.
Background technology
Using optical fiber can realize distance detection for no reason at all as sensor, small disturbance is the most permissible Making optical fiber deform due to vibration, transmission locus also can respective change.
Phase sensitive optical time domain reflectometerIt it is a kind of novel distributing optical fiber sensing Device, in terms of the vibration signal in the range of fibre circuit is carried out long-range and real-time monitoring, has biography The irreplaceable advantage of sensor of system.With the conventional optical time domain reflection detected based on scattering strength Instrument (OTDR) is compared,It is possible not only to the intensity of detection fiber Rayleigh scattering echo, But also the phase fluctuation of echo can be detected, there is higher sensitivity and detect Dynamic Signal Ability.In recent years,Fiber fence, in terms of circumference security protection, receives more and more extensive Concern.FromThe large amount of complex signal detected judging, how optical fiber there occurs along the line Disturbance, the how invasion of character, this is the target of circumference security protection application actual concern.Therefore, The disturbing signal detected is carried out preliminary identification extremely important, in order to make quickly and accurately Response.
Summary of the invention
The technical problem that present invention mainly solves be to provide a kind of fiber-optic vibration model determine method and Device, it is possible to accurately judge the oscillatory type of target fiber position, and improve fiber-optic vibration detection Efficiency.
For solving above-mentioned technical problem, the technical scheme that the present invention uses is: provide a kind of light Fine model of vibration determines that method, the method include: obtain the vibration signal of target fiber position；Carry Take the spectral characteristic of vibration signal；The spectral characteristic of spectral characteristic Yu default model of vibration is carried out Join；When the spectral characteristic of spectral characteristic with default model of vibration is mated, determine that vibration signal belongs to Preset model of vibration.
Wherein, the spectral characteristic of spectral characteristic with default model of vibration is mated, including: meter Calculate the Euclidean distance between spectrum signature and the spectral characteristic of default model of vibration of vibration signal；Sentence Whether disconnected Euclidean distance meets predeterminable range threshold value；If meeting, it is determined that spectral characteristic is shaken with presetting The spectral characteristic coupling of movable model.
Wherein, extract the spectral characteristic of vibration signal, including: vibration signal is normalized place Reason；Vibration signal after normalized is carried out Fourier's change, to obtain the frequency of vibration signal Spectrum signature.
Wherein, obtain the vibration signal of target fiber position, including: target fiber position is reflected Fiber-optic signal in time domain, be divided into multiple signal frame；Judge whether vibrate overall strength or shake Dynamic intensity of variation meets pre-conditioned signal frame；If existing, then continuous print is met pre-conditioned Signal frame be defined as vibration signal.
Wherein, it may be judged whether there is vibration overall strength or vibration intensity of variation meets pre-conditioned letter Number frame, including: judge whether short-time energy be more than energy threshold and/or in short-term translation difference big Signal frame in differential threshold.
Wherein, in multiple signal frames, two adjacent signal frames partly overlap in time.
Wherein, fiber-optic signal is beat signal, by target fiber position reflect fiber-optic signal time Before being divided into multiple signal frame on territory, also include: optical signal is divided into first via optical signal and Second road optical signal；After being amplified by first via optical signal modulation, it is passed through optical fiber；First via light is believed Number reflect the reflected light signal obtained in a fiber and the second road optical signal couples；To couple To optical signal carry out opto-electronic conversion, to obtain the beat signal of target fiber position.
For solving above-mentioned technical problem, another technical solution used in the present invention is: provide one Fiber-optic vibration model determines device, and this device includes: signal acquisition module, is used for obtaining target light The vibration signal of fine position；Frequency spectrum extraction module, for extracting the spectral characteristic of vibration signal；? Join module, for the spectral characteristic of spectral characteristic with default model of vibration being mated；Vibration mould Type determines module, for when the spectral characteristic of spectral characteristic with default model of vibration is mated, determines Vibration signal belongs to default model of vibration.
Wherein, matching module includes: computing unit, for calculate the spectrum signature of vibration signal with Preset the Euclidean distance between the spectral characteristic of model of vibration；Judging unit, be used for judging European away from From whether meeting predeterminable range threshold value；Determine unit, for meeting predeterminable range threshold at Euclidean distance During value, determine that spectral characteristic is mated with the spectral characteristic of default model of vibration.
Wherein, frequency spectrum extraction module specifically for: vibration signal is normalized；To returning Vibration signal after one change processes carries out Fourier's change, to obtain the spectrum signature of vibration signal.
The invention has the beneficial effects as follows: be different from the situation of prior art, the fiber-optic vibration of the present invention Model determines that method includes: obtain the vibration signal of target fiber position；Extract the frequency of vibration signal Spectral property；The spectral characteristic of spectral characteristic with default model of vibration is mated；In spectral characteristic When mating with the spectral characteristic of default model of vibration, determine that vibration signal belongs to default model of vibration. By the way, the signal of time domain is converted into the signal of frequency domain, by calculating Euclidean distance pair Ratio, it is possible to accurately judge the oscillatory type of target fiber position, and improve the effect of fiber-optic vibration detection Rate.
Accompanying drawing explanation
Fig. 1 is the schematic flow sheet that fiber-optic vibration model of the present invention determines method the first embodiment；
Fig. 2 is the schematic flow sheet that fiber-optic vibration model of the present invention determines method the second embodiment；
Fig. 3 is the signal flow diagram that fiber-optic vibration model of the present invention determines method the second embodiment；
Fig. 4 is that fiber-optic vibration model of the present invention is when determining beat signal in method the second embodiment Territory signal schematic representation；
Fig. 5 is the schematic flow sheet that fiber-optic vibration model of the present invention determines method the 3rd embodiment；
Fig. 6 is distance-vibrational energy that fiber-optic vibration model of the present invention determines method one object lesson Schematic diagram；
Fig. 7 is short-time energy and the translation that fiber-optic vibration model of the present invention determines method one object lesson The schematic diagram of difference；
Fig. 8 is that fiber-optic vibration model of the present invention determines time domain letter after the normalization of method one object lesson Number schematic diagram；
Fig. 9 is the spectrum diagram that fiber-optic vibration model of the present invention determines method one object lesson；
Figure 10 is that fiber-optic vibration model of the present invention determines method one object lesson test result and time The contrast schematic diagram of the result that dynamic programming calculates；
Figure 11 is the structural representation that fiber-optic vibration model of the present invention determines device one embodiment；
Figure 12 is the structural representation that fiber-optic vibration model of the present invention determines equipment one embodiment.
Detailed description of the invention
Refering to Fig. 1, fiber-optic vibration model of the present invention determines the flow process signal of method the first embodiment Figure, the method includes:
S11: obtain the vibration signal of target fiber position.
Wherein, vibration signal can be the signal of time domain, and it represents that target fiber position is when difference Between Oscillation Amplitude.Optionally, the amplitude of vibration signal can be converted into the signal of telecommunication with optical signal Intensity size, it is also possible to analog electrical signal is converted into the Digital size of digital signal.
S12: extract the spectral characteristic of vibration signal.
Optionally, the vibration signal of time domain can be converted into the vibration signal of frequency domain.Wherein, shake The spectral characteristic of dynamic signal can be the frequency spectrum function of vibration signal, it is also possible to signal is carried out frequency spectrum Analyze, the amplitude of signal, phase place or energy conversion are represented with frequency coordinate axle, and then analyzes it Frequency characteristic.Signal is carried out spectrum analysis and can obtain more useful information, as tried to achieve dynamically letter Each frequency content in number and frequency distribution scope, obtain each frequency content amplitude distribution and Energy distribution, thus obtain the frequency values of main amplitude and Energy distribution.
S13: the spectral characteristic of spectral characteristic with default model of vibration is mated.
Optionally, the model of vibration preset can be in the amplitude in frequency spectrum, phase place or energy conversion The characteristic of some parameter, it is also possible to be characteristics of two or more combinations.Such as, may be used To arrange multiple threshold value to amplitude, phase place or energy conversion, when meeting the threshold value preset, i.e. frequency Spectral property mates with the spectral characteristic preset.
Optionally, model of vibration can also is that the curve chart drawn according to spectral characteristic, when vibration letter Number spectrum curve similar, i.e. it can be assumed that shake with some spectrum curve of model of vibration preset The spectral characteristic of dynamic signal is mated with the spectral characteristic of default model of vibration.
S14: when the spectral characteristic of spectral characteristic with default model of vibration is mated, determines vibration letter Number belong to default model of vibration.
Optionally, model of vibration can be according to the power of vibration according to grade separation, and such as one-level is shaken Dynamic, secondary vibration, three grades of vibrations etc..Can also classify according to the really scene causing vibration, Such as pedestrian walks, automobile is exercised, construction or natural wind etc..
Being different from prior art, the fiber-optic vibration model of present embodiment determines that method includes: obtain The vibration signal of target fiber position；Extract the spectral characteristic of vibration signal；By spectral characteristic with pre- If the spectral characteristic of model of vibration is mated；Special with the frequency spectrum of default model of vibration in spectral characteristic Property coupling time, determine that vibration signal belongs to default model of vibration.By the way, by time domain Signal is converted into the signal of frequency domain, contrasts by calculating Euclidean distance, it is possible to accurately judge target light The oscillatory type of fine position, and improve the efficiency of fiber-optic vibration detection.
It is that fiber-optic vibration model of the present invention determines method the second embodiment party refering to Fig. 2 and Fig. 3, Fig. 2 The flow chart of formula, Fig. 3 is the signal that fiber-optic vibration model of the present invention determines method the second embodiment Flow graph, wherein, connecting line represents optical fiber, two ellipse representation bonders, circular expression optical fiber end Mouthful, for being passed through optical signal or receiving the optical signal of reflection.
Present embodiment discloses the idiographic flow of S11 in above-mentioned first embodiment, and this S11 has Body includes:
S111: optical signal is divided into first via optical signal and the second road optical signal.
Optionally, optical fiber laser 31 can use distributed feedback fibre laser (Distributed Feedback Laser, DFB).
S112: after being amplified by first via optical signal modulation, be passed through optical fiber.
It is alternatively possible to by first via optical signal by acousto-optic modulator (AOM) 32 frequency modulation, It is modulated to the light pulse of certain pulse width, then after fiber amplifier (EDFA0) 33 amplifies Injection fibre.
S113: first via optical signal is reflected in a fiber the reflected light signal obtained and the second road light Signal couples.
Optionally, first via optical signal is reflected in a fiber the reflected light signal obtained and the second tunnel Optical signal interferes superposition on the bonder of 3dB.
S114: optical signal coupling obtained carries out opto-electronic conversion, to obtain target fiber position Beat signal.
Optionally, the signal after interfering superposition can be passed through photodetector 34, carry out photoelectricity and turn Change.Optionally, photodetector 34 can be double flat weighing apparatus photodetector (DB-PD), is used for Receive the output of bonder, obtain beat signal.
Optionally, after obtaining beat signal, it is also possible to by data collecting card 35, data are passed Deliver to processor 36, in order to the follow-up Treatment Analysis to data.Wherein, data collecting card 35 He Processor 36 can also be a part for a terminal unit, such as computer.
S115: beat signal is divided in time domain multiple signal frame.
Refering to Fig. 4, fiber-optic vibration model of the present invention determines beat signal in method the second embodiment Time-domain signal schematic diagram.Wherein, abscissa n express time, vertical coordinate x (n) represents vibration Energy.
Optionally, can start to divide from 0 moment, every m time span is a signal frame, I.e. 0-m is a signal frame, and m-2m is a signal frame, and 2m-3m is a signal frame, 3m-4m It it is a signal frame.In the diagram, m is 5, the most every 5 signaling points, one signal frame of composition. In other embodiments, it is also possible to any time from beat signal starts to choose signal frame.
Optionally, in multiple signal frames, two adjacent signal frames partly overlap in time.Weight Folded length can be the 1/3 of signal frame length, 1/4,1/5 etc. chosen, and not can exceed that 1/2, such as overlap length is 1/5 signal frame, then first signal frame is 0-4 totally 5, Second signal frame is then 4-8 totally 5, and wherein the 4th signaling point is lap, so has It is beneficial to ensure the seriality of signal frame.
S116: judge whether to vibrate overall strength or vibration intensity of variation meets pre-conditioned letter Number frame.
Wherein, it is judged that the vibration overall strength of signal frame can directly reflect that the vibration of this signal frame is big Little, it is judged that the vibration intensity of variation of signal frame can reflect the fluctuating quantity of amplitude.Implementing In, can be using one of them as criterion, it is also possible to both are combined as criterion.
Optionally, in a specific embodiment, it can be determined that whether there is short-time energy and be more than Energy threshold and/or in short-term translation difference are more than the signal frame of differential threshold.
Specifically, as a example by signal frame m-2m, short-time energy can be to be obtained by equation below:
E = Σ n = m 2 m - 1 x 2 ( n )
Wherein, this formula represents the quadratic sum of each energy value in signal frame m-2m.
Specifically, as a example by signal frame m-2m, translation difference can be to pass through equation below in short-term Obtain:
D = Σ n = m 2 m - 1 | x ( n ) - x ( n - 1 ) |
Wherein, this formula represents each energy value and the difference of previous energy value in signal frame m-2m The sum of value.
S117: if existing, then continuous print is met pre-conditioned signal frame and is defined as vibration signal.
The signal frame that continuous print meets above-mentioned condition is defined as vibration signal, and carries out first real Execute S12 and follow-up step in mode.
Refering to Fig. 5, fiber-optic vibration model of the present invention determines the flow process signal of method the 3rd embodiment Figure, present embodiment is mainly striven another embodiment of S12, S13 in the first embodiment, Wherein, S12 specifically may include that
S121: vibration signal is normalized.
Specifically, in order to when follow-up and default model of vibration spectral characteristic contrast, have phase Same amplitude criteria, needs exist for being normalized vibration signal, the signal frame i.e. chosen In each data divided by the maximum in these data, these data are all normalized to [-1,1] Interval in.
S122: the vibration signal after normalized is carried out Fourier's change, to obtain vibration letter Number spectrum signature.
Specifically, as a example by two signal frame m-3m altogether 2m signal value, Fourier transformation Expression formula is as follows:
F ( k ) = Σ n = 1 2 m x ( n ) e - j 2 π k n 2 m
Wherein, two signal frame length of seasonal effect in time series length, i.e. m-2m after 2m is extraction. Spectral data sequence k=[1 ..., 2m].Signal amplitude after Fourier transformation is as entering Invade the characteristic signal of signal.Through experiment test, the spectrum signature extracted is the most stable, mutually of the same race The spectrum signature of class signal is basically identical, and the impact that disturbance individual variation and noise introduce is relatively Little.Therefore, as characteristic signal, it is possible not only to distinguish different types of signal, Er Qieli In eliminating environment noise, the impact of individual invasion difference.
Wherein, S13 specifically may include that
S131: calculate between spectrum signature and the spectral characteristic of default model of vibration of vibration signal Euclidean distance.
Specifically, as a example by two signal frame m-3m altogether 2m signal value, Euclidean distance is calculated Formula as follows:
d = { Σ k = 1 2 m [ F s i g n a l ( k ) - F mod e l ( k ) ] 2 } 1 / 2
Wherein, FsignalAnd FmodelIt is respectively the frequency spectrum of vibration signal to be measured and default model of vibration Frequency spectrum, 2m is the length of spectrum sequence.
S132: judge whether Euclidean distance meets predeterminable range threshold value.
S133: if meeting, it is determined that spectral characteristic is mated with the spectral characteristic of default model of vibration.
The frequency spectrum of vibration signal to be measured and the frequency spectrum of different default models of vibration are brought into respectively Above-mentioned formula calculates, it will obtains different Euclidean distance numerical value, it is entered with predeterminable range threshold value Row compares, and can realize the quick identification of disturbing signal.
It should be understood that different Euclidean distance threshold values can be set for different models of vibration.
Below, with a specific example, embodiments of the present invention are described in detail:
Three kinds of disturbances are artificially set in target fiber position, including: A, play wall (toeing), B, Kick wall (heeling), C, running on spot (running).Target fiber length is 5m, distance optical fiber The distance of end is 4.3km.
The pulsewidth assuming light pulse is 100ns, by obtaining reflected signal and being calculated by the time Distance, can obtain the vibrational energy of different fiber position, shown in Fig. 6, wherein, the horizontal stroke of Fig. 6 Coordinate representation fiber distance, vertical coordinate represents vibrational energy.
Extract the time-domain signal of different fiber position the most respectively, this time domain of target fiber position is believed Number it is divided into multiple continuous print signal frame, calculates short-time energy and the translation difference of each signal frame, The curve as shown in Figure 7 obtained.Wherein abscissa represents signal frame number, and vertical coordinate represents signal Energy, curve above represents short-time energy, and curve below represents translation difference, EthRepresent short Time energy threshold, DthRepresent translation differential threshold.It can be seen from figure 7 that only A, B, Short-time energy at C tri-is more than Eth, and this translation difference at three is also greater than Dth, it is believed that This is exactly playing wall, kicking wall and three kinds of disturbances of running on spot of above-mentioned artificial setting at three.
Extract at above-mentioned three kinds of disturbances respectively and the time series of natural wind (wind) disturbance, and Data are normalized, as shown in Figure 8, wherein, abscissa express time, vertical coordinate table Indicating calibration signal intensity, from top to bottom four curves represent successively play wall, running on spot, kick wall and from So wind.Again aforementioned four time-domain signal is converted into frequency domain, shown in Fig. 9, the abscissa of Fig. 9 Representing frequency, vertical coordinate represents frequency spectrum, four curves from top to down of the enlarged fragmentary portion in Fig. 9 Represent natural wind successively, play wall, kick wall and running on spot.
As shown in Figure 10, calculate the Euclidean distance between frequency spectrum and the frequency spectrum of model obtained, and Contrast.Wherein, (a1), (b1), (c1) are for utilizing time dynamic programming (DTW) to calculate The result of method method；(a2), (b2), (c2) for utilize frequency spectrum Euclidean distance method (EDFS, i.e. Present embodiment).Abscissa represents result sequence, and vertical coordinate represents Euclidean distance.
As a example by running in place class, the threshold value of EDFS is 1000, from (c2) figure it can be seen that The data of all running in place results the most correctly identify；See again (c1), and in DTW algorithm two Point occurs in that erroneous judgement, and this European more than other disturbance types of Euclidean distance at erroneous judgement at two Distance, this erroneous judgement is soluble not by adjusting threshold value.
Optical fiber intrusion detection system is highly prone to the impact of natural wind in the application.To this end, have studied The feature of wind, as shown in Figure 8,9, 10.It can be seen that the frequency ratio of normal conditions leeward enters The frequency invading signal is low, utilizes spectrum analysis to be easier to eliminate the impact of environment noise.EDFS is Based on short time spectrum analysis, it is easier to eliminating the interference of low-frequency noise, the impact of white noise is also Can reduce.
Additionally, EDFS algorithm for pattern recognition is not only effective, and efficiently.Result shows, in phase Process identical data EDFS under the conditions of with and have only to 0.11s, and DTW algorithm needs 5.05 s。
What deserves to be explained is, Fig. 6-Figure 10 is only the test process of the present invention one object lesson, The principle of the present invention first, second, third embodiment is described in detail by main users, its In data, figure be not limiting as protection scope of the present invention.
Refering to Figure 11, fiber-optic vibration model of the present invention determines the structural representation of device one embodiment, This device includes:
Signal acquisition module 1101, for obtaining the vibration signal of target fiber position.
Frequency spectrum extraction module 1102, for extracting the spectral characteristic of vibration signal.
Matching module 1103, for carrying out the spectral characteristic of spectral characteristic Yu default model of vibration Join.
Model of vibration determines module 1104, for special with the frequency spectrum of default model of vibration in spectral characteristic Property coupling time, determine that vibration signal belongs to default model of vibration.
Optionally, in other embodiments, matching module 1103 can also include: calculates single Unit, for calculating the Europe between the spectrum signature of vibration signal and the spectral characteristic of default model of vibration Formula distance；Judging unit, is used for judging whether Euclidean distance meets predeterminable range threshold value；Determine list Unit, for when Euclidean distance meets predeterminable range threshold value, determines spectral characteristic and presets vibration mould The spectral characteristic coupling of type.
Optionally, in other embodiments, frequency spectrum extraction module specifically for: to vibration signal It is normalized；Vibration signal after normalized is carried out Fourier's change, to obtain The spectrum signature of vibration signal.
Refering to Figure 12, fiber-optic vibration model of the present invention determines the structural representation of equipment one embodiment, This equipment 120 includes processor 1201, memorizer 1202, receptor 1203 and transmitter 1204.
Optionally, this processor 1201, memorizer 1202, receptor 1203 and transmitter 1204 Connected by a bus.
Memorizer 1202 is used for storing system file, application software and various algorithm, threshold value letter Breath and history detection record etc..
Receptor 1203 is used for receiving vibration signal, and optionally, receptor 1203 can be to gather Card, is used for gathering vibration signal and being sent to processor 1201；Transmitter 1204 is used for sending place The result of reason device 1201 is to the outut device such as display, siren.
Processor 1201 is used for performing following steps:
Obtain the vibration signal of target fiber position；Extract the spectral characteristic of vibration signal；By frequency spectrum Characteristic is mated with the spectral characteristic of default model of vibration；At spectral characteristic and default model of vibration Spectral characteristic coupling time, determine that vibration signal belongs to default model of vibration.
Optionally, processor 1201 is additionally operable to perform: calculate the spectrum signature of vibration signal with pre- If the Euclidean distance between the spectral characteristic of model of vibration；Judge Euclidean distance whether meet preset away from From threshold value；If meeting, it is determined that spectral characteristic is mated with the spectral characteristic of default model of vibration.
Optionally, processor 1201 is additionally operable to perform: be normalized vibration signal； Vibration signal after normalized is carried out Fourier's change, special to obtain the frequency spectrum of vibration signal Levy.
Optionally, processor 1201 is additionally operable to perform: optical fiber letter target fiber position reflected Number in time domain, it is divided into multiple signal frame；Judge whether to vibrate overall strength or vibration change journey Degree meets pre-conditioned signal frame；If existing, then continuous print is met pre-conditioned signal frame It is defined as vibration signal.
Optionally, processor 1201 is additionally operable to perform: judge whether that short-time energy is more than energy Amount threshold value and/or in short-term translation difference are more than the signal frame of differential threshold.
In several embodiments provided by the present invention, it should be understood that disclosed method with And equipment, can realize by another way.Such as, equipment embodiment described above It is only schematically, such as, described module or the division of unit, it is only a kind of logic function Dividing, actual can have other dividing mode, the most multiple unit or assembly to tie when realizing Close or be desirably integrated into another system, or some features can be ignored, or not performing.
The described unit illustrated as separating component can be or may not be physically separate , the parts shown as unit can be or may not be physical location, i.e. may be located at One place, or can also be distributed on multiple NE.Can select according to the actual needs Select some or all of unit therein to realize the purpose of present embodiment scheme.
It addition, each functional unit in each embodiment of the present invention can be integrated in a process In unit, it is also possible to be that unit is individually physically present, it is also possible to two or more unit Integrated in a unit.Above-mentioned integrated unit both can realize to use the form of hardware, it is possible to To use the form of SFU software functional unit to realize.
If the integrated unit in other embodiments above-mentioned is real with the form of SFU software functional unit Now and as independent production marketing or when using, the storage of embodied on computer readable can be stored in In medium.Based on such understanding, technical scheme is the most in other words to prior art The part contributed or this technical scheme completely or partially can be with the form of software product Embodying, this computer software product is stored in a storage medium, uses including some instructions So that a computer equipment (can be personal computer, server, or the network equipment etc.) Or processor (processor) performs all or part of of method described in each embodiment of the present invention Step.And aforesaid storage medium includes: USB flash disk, portable hard drive, read only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), The various medium that can store program code such as magnetic disc or CD.
The foregoing is only embodiments of the present invention, not thereby limit the scope of the claims of the present invention, Every equivalent structure utilizing description of the invention and accompanying drawing content to be made or equivalence flow process conversion, or Directly or indirectly being used in other relevant technical fields, the patent being the most in like manner included in the present invention is protected In the range of protecting.

Claims (10)

1. a fiber-optic vibration model determines method, it is characterised in that including:
Obtain the vibration signal of target fiber position；
Extract the spectral characteristic of described vibration signal；
The spectral characteristic of described spectral characteristic with default model of vibration is mated；
When the spectral characteristic of described spectral characteristic with default model of vibration is mated, determine described vibration Signal belongs to described default model of vibration.
Method the most according to claim 1, it is characterised in that described by described spectral characteristic Mate with the spectral characteristic of default model of vibration, including:
Calculate between the spectrum signature of described vibration signal and the spectral characteristic of default model of vibration Euclidean distance；
Judge whether described Euclidean distance meets predeterminable range threshold value；
If meeting, it is determined that described spectral characteristic is mated with the spectral characteristic of default model of vibration.
Method the most according to claim 2, it is characterised in that described extraction described vibration letter Number spectral characteristic, including:
Described vibration signal is normalized；
Vibration signal after normalized is carried out Fourier's change, to obtain described vibration signal Spectrum signature.
Method the most according to claim 1, it is characterised in that described acquisition target optical fiber position The vibration signal put, including:
The fiber-optic signal that target fiber position reflects is divided into multiple signal frame in time domain；
Judge whether to vibrate overall strength or vibration intensity of variation meets pre-conditioned described letter Number frame；
If existing, then continuous print is met pre-conditioned signal frame and is defined as described vibration signal.
Method the most according to claim 4, it is characterised in that described in judge whether to shake Dynamic overall strength or vibration intensity of variation meet pre-conditioned described signal frame, including:
Judge whether that short-time energy is more than energy threshold and/or translates difference in short-term more than difference The described signal frame of threshold value.
Method the most according to claim 4, it is characterised in that in the plurality of signal frame, Two adjacent signal frames partly overlap in time.
Method the most according to claim 4, it is characterised in that described fiber-optic signal is beat frequency Signal, the described fiber-optic signal by the reflection of target fiber position is divided into multiple signal frame in time domain Before, also include:
Optical signal is divided into first via optical signal and the second road optical signal；
After being amplified by described first via optical signal modulation, it is passed through optical fiber；
Described first via optical signal is reflected in described optical fiber the reflected light signal obtained with described Second road optical signal couples；
Optical signal coupling obtained carries out opto-electronic conversion, to obtain the described bat of target fiber position Frequently signal.
8. a fiber-optic vibration model determines device, it is characterised in that including:
Signal acquisition module, for obtaining the vibration signal of target fiber position；
Frequency spectrum extraction module, for extracting the spectral characteristic of described vibration signal；
Matching module, for carrying out the spectral characteristic of described spectral characteristic Yu default model of vibration Join；
Model of vibration determines module, for special with the frequency spectrum of default model of vibration in described spectral characteristic Property coupling time, determine that described vibration signal belongs to described default model of vibration.
Device the most according to claim 8, it is characterised in that described matching module includes:
Computing unit, for calculating the spectrum signature of described vibration signal and the frequency of default model of vibration Euclidean distance between spectral property；
Judging unit, is used for judging whether described Euclidean distance meets predeterminable range threshold value；
Determine unit, for when described Euclidean distance meets predeterminable range threshold value, determine described frequency Spectral property mates with the spectral characteristic of default model of vibration.
Device the most according to claim 9, it is characterised in that described frequency spectrum extraction module has Body is used for:
Described vibration signal is normalized；
Vibration signal after normalized is carried out Fourier's change, to obtain described vibration signal Spectrum signature.
CN201610301609.1A 2016-05-09 2016-05-09 Optical fiber vibration model determination method and device Pending CN105973451A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610301609.1A CN105973451A (en) 2016-05-09 2016-05-09 Optical fiber vibration model determination method and device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610301609.1A CN105973451A (en) 2016-05-09 2016-05-09 Optical fiber vibration model determination method and device

Publications (1)

Publication Number Publication Date
CN105973451A true CN105973451A (en) 2016-09-28

Family

ID=56992785

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610301609.1A Pending CN105973451A (en) 2016-05-09 2016-05-09 Optical fiber vibration model determination method and device

Country Status (1)

Country Link
CN (1) CN105973451A (en)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106503642A (en) * 2016-10-18 2017-03-15 长园长通新材料股份有限公司 A kind of model of vibration method for building up for being applied to optical fiber sensing system
CN107727227A (en) * 2017-09-30 2018-02-23 南京大学 Ultra-high-tension power transmission line icing based on Φ OTDR waves monitoring method
CN109974836A (en) * 2019-04-09 2019-07-05 苏州珈全智能科技有限公司 A kind of device and method improving φ-OTDR frequency response
CN110135283A (en) * 2019-04-25 2019-08-16 上海大学 The signal recognition method of optical fiber perimeter defence system based on FastDTW algorithm
CN110440901A (en) * 2019-08-13 2019-11-12 郑州信大先进技术研究院 A kind of distributed optical fiber vibration sensing localization method and device based on pulse accumulation
CN111537056A (en) * 2020-07-08 2020-08-14 浙江浙能天然气运行有限公司 Pipeline along-line third-party construction dynamic early warning method based on SVM and time-frequency domain characteristics
CN112504429A (en) * 2020-11-24 2021-03-16 姚峰 High-precision demodulation algorithm for strong interference DVS
CN112539772A (en) * 2020-11-02 2021-03-23 上海大学 Positioning method of Sagnac distributed optical fiber sensing system based on convolutional neural network integrated learning

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101441092A (en) * 2008-10-16 2009-05-27 北京邮电大学 Perimeter protection sensing positioning system based on coherent light time domain reflection
CN101566497A (en) * 2009-04-29 2009-10-28 上海华魏光纤传感技术有限公司 Distribution-type fiber vibration sensor system based on phase detection and optical time domain reflection
CN101603856A (en) * 2009-07-16 2009-12-16 上海华魏光纤传感技术有限公司 A kind of long-distance distributed optical fiber vibration sensing system and method
WO2012063066A2 (en) * 2010-11-11 2012-05-18 Fotech Solutions Limited Distributed optical fibre sensor
US20140090474A1 (en) * 2009-09-28 2014-04-03 At&T Intellectual Property I, L.P. Long Distance Optical Fiber Sensing System and Method
CN105389917A (en) * 2015-09-18 2016-03-09 南京派光信息技术有限公司 Rapid early-warning method based on phase-sensitive optical time-domain reflectometer

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101441092A (en) * 2008-10-16 2009-05-27 北京邮电大学 Perimeter protection sensing positioning system based on coherent light time domain reflection
CN101566497A (en) * 2009-04-29 2009-10-28 上海华魏光纤传感技术有限公司 Distribution-type fiber vibration sensor system based on phase detection and optical time domain reflection
CN101603856A (en) * 2009-07-16 2009-12-16 上海华魏光纤传感技术有限公司 A kind of long-distance distributed optical fiber vibration sensing system and method
US20140090474A1 (en) * 2009-09-28 2014-04-03 At&T Intellectual Property I, L.P. Long Distance Optical Fiber Sensing System and Method
WO2012063066A2 (en) * 2010-11-11 2012-05-18 Fotech Solutions Limited Distributed optical fibre sensor
CN105389917A (en) * 2015-09-18 2016-03-09 南京派光信息技术有限公司 Rapid early-warning method based on phase-sensitive optical time-domain reflectometer

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
刘建霞: "Φ-OTDR分布式光纤传感监测技术的研究进展", 《激光与光电子学进展》 *

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106503642A (en) * 2016-10-18 2017-03-15 长园长通新材料股份有限公司 A kind of model of vibration method for building up for being applied to optical fiber sensing system
CN106503642B (en) * 2016-10-18 2019-09-20 长园长通新材料股份有限公司 A kind of model of vibration method for building up applied to optical fiber sensing system
CN107727227A (en) * 2017-09-30 2018-02-23 南京大学 Ultra-high-tension power transmission line icing based on Φ OTDR waves monitoring method
CN107727227B (en) * 2017-09-30 2019-05-21 南京大学 Ultra-high-tension power transmission line icing based on Φ-OTDR waves monitoring method
CN109974836A (en) * 2019-04-09 2019-07-05 苏州珈全智能科技有限公司 A kind of device and method improving φ-OTDR frequency response
CN110135283A (en) * 2019-04-25 2019-08-16 上海大学 The signal recognition method of optical fiber perimeter defence system based on FastDTW algorithm
CN110440901A (en) * 2019-08-13 2019-11-12 郑州信大先进技术研究院 A kind of distributed optical fiber vibration sensing localization method and device based on pulse accumulation
CN110440901B (en) * 2019-08-13 2021-08-17 郑州信大先进技术研究院 Distributed optical fiber vibration sensing positioning method and device based on pulse accumulation
CN111537056A (en) * 2020-07-08 2020-08-14 浙江浙能天然气运行有限公司 Pipeline along-line third-party construction dynamic early warning method based on SVM and time-frequency domain characteristics
CN112539772A (en) * 2020-11-02 2021-03-23 上海大学 Positioning method of Sagnac distributed optical fiber sensing system based on convolutional neural network integrated learning
CN112504429A (en) * 2020-11-24 2021-03-16 姚峰 High-precision demodulation algorithm for strong interference DVS

Similar Documents

Publication Publication Date Title
CN105973451A (en) Optical fiber vibration model determination method and device
CN106225907B (en) It is a kind of based on Φ-OTDR technique fiber-optic vibration identifying system and method
Tejedor et al. Toward prevention of pipeline integrity threats using a smart fiber-optic surveillance system
CN111854921A (en) Distributed optical fiber deceleration strip vibration early warning system and method
CN103183263B (en) The safety and protection system of works and possess its elevator
JP6974747B2 (en) Manhole position identification method and manhole position identification system
CN111373252A (en) Bridge damage rapid detection method and related device
CN103116957B (en) A kind of method of optical fiber perimeter security and protection system shielding climate impact
CN104538041A (en) Method and system for detecting abnormal sounds
CN102629470B (en) Consonant-segment detection apparatus and consonant-segment detection method
CN107782438A (en) Optical fiber raster vibration measuring system based on pulse code and edge filter demodulation method
CN104504911B (en) A kind of speed measurer for motor vehicle and speed-measuring method
CN107525849A (en) A kind of single-input single-output test modal analysis system and method based on fiber grating
CN107026731A (en) A kind of method and device of subscriber authentication
CN115100816A (en) Wearable near-electricity sensing equipment and monitoring system
CN109974836A (en) A kind of device and method improving φ-OTDR frequency response
CN110595599B (en) Method for reducing polarization fading of optical fiber vibration system and detection system applying same
CN106197646A (en) The detection of a kind of fiber-optic vibration reduces the method for error and fine vibration detection device
CN106706109B (en) Vibration source identification method and system based on time domain two-dimensional characteristics
CN105651373A (en) Two-point synchronous vibration measuring method based on polarization optical time-domain reflection technology
JP2024514865A (en) Locating hazardous events on utility poles
CN106501843A (en) A kind of identification determination methods to magnitude
CN109612568A (en) A kind of mobile method for interference source identification of vibration source
CN104964736B (en) Optical fiber invasion vibration source identification method based on time-frequency characteristic maximum expected classification
CN103886868A (en) Impact sound detection method and detection 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
WD01 Invention patent application deemed withdrawn after publication
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20160928