CN110044397A - A kind of quantitatively measuring device and its method of Real-Time Optical fiber disturbance sensing - Google Patents
A kind of quantitatively measuring device and its method of Real-Time Optical fiber disturbance sensing Download PDFInfo
- Publication number
- CN110044397A CN110044397A CN201910191993.8A CN201910191993A CN110044397A CN 110044397 A CN110044397 A CN 110044397A CN 201910191993 A CN201910191993 A CN 201910191993A CN 110044397 A CN110044397 A CN 110044397A
- Authority
- CN
- China
- Prior art keywords
- data
- fpga module
- ddr3
- disturbance
- optical fiber
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 42
- 239000013307 optical fiber Substances 0.000 title claims abstract description 35
- 230000003287 optical effect Effects 0.000 claims abstract description 33
- 238000001514 detection method Methods 0.000 claims abstract description 19
- 238000003909 pattern recognition Methods 0.000 claims abstract description 11
- 230000008569 process Effects 0.000 claims abstract description 10
- 239000000835 fiber Substances 0.000 claims description 23
- 239000011159 matrix material Substances 0.000 claims description 17
- 238000001914 filtration Methods 0.000 claims description 16
- 230000010355 oscillation Effects 0.000 claims description 12
- 238000004422 calculation algorithm Methods 0.000 claims description 10
- 238000012545 processing Methods 0.000 claims description 7
- 238000005070 sampling Methods 0.000 claims description 5
- 230000017105 transposition Effects 0.000 claims description 3
- 235000013399 edible fruits Nutrition 0.000 claims 1
- 230000007704 transition Effects 0.000 abstract description 4
- 238000013527 convolutional neural network Methods 0.000 abstract description 3
- 238000003672 processing method Methods 0.000 abstract description 2
- 230000035559 beat frequency Effects 0.000 description 8
- 230000006872 improvement Effects 0.000 description 8
- 238000004458 analytical method Methods 0.000 description 6
- 230000008859 change Effects 0.000 description 5
- 102220492046 Phospholipid scramblase 1_H12A_mutation Human genes 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 3
- 230000005622 photoelectricity Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 230000011664 signaling Effects 0.000 description 3
- 238000001228 spectrum Methods 0.000 description 3
- 230000003068 static effect Effects 0.000 description 3
- 241001269238 Data Species 0.000 description 2
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 2
- 238000013528 artificial neural network Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 229940085805 fiberall Drugs 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 238000012546 transfer Methods 0.000 description 2
- 230000003321 amplification Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 210000004556 brain Anatomy 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000003628 erosive effect Effects 0.000 description 1
- 230000001939 inductive effect Effects 0.000 description 1
- 238000009413 insulation Methods 0.000 description 1
- 229910052742 iron Inorganic materials 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000003058 natural language processing Methods 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 238000000253 optical time-domain reflectometry Methods 0.000 description 1
- 230000008439 repair process Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 230000001953 sensory effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01D—MEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
- G01D5/00—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable
- G01D5/26—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light
- G01D5/32—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light
- G01D5/34—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light the beams of light being detected by photocells
- G01D5/353—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light the beams of light being detected by photocells influencing the transmission properties of an optical fibre
- G01D5/35338—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light the beams of light being detected by photocells influencing the transmission properties of an optical fibre using other arrangements than interferometer arrangements
- G01D5/35354—Sensor working in reflection
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
The invention discloses the quantitatively measuring devices and its method of a kind of Real-Time Optical fiber disturbance sensing, including phase sensitive optical time domain reflectometer；Double-channel analog/digital converter；FPGA module；DDR3 random access memory；Industrial control computer and GPU.Compared with prior art, the present invention demodulates the phase of phase sensitive optical time domain reflectometer using the mode of heterodyne detection, realizes the quantitative measurment to strength of turbulence；Its data processing method uses sliding window mode, and the transition that can be realized disturbance is shown, improves the readability of result；The characteristics of using novel fpga chip large capacity, high concurrent speed and DDR3 random access memory high bandwidth, high capacity, realize the multi-level buffer of data, multi-channel parallel calculates and pipeline processes, solving the problems, such as optical fiber disturbance sensor-based system over long distances and at high speed cannot realize simultaneously, convolutional neural networks pattern recognition and classification is carried out to disturbance using GPU simultaneously, there is preferable real-time and accuracy.
Description
Technical field
The invention belongs to sensory field of optic fibre and digital processing field, are related to a kind of determining for Real-Time Optical fiber disturbance sensing
Measuring device and its method.
Background technique
Distributed Optical Fiber Sensing Techniques not only have the advantage of ordinary optic fibre sensing technology, such as electromagnetism interference, concealment
Good, erosion resisting insulation etc., and the continuous of all kinds of physical quantitys (temperature, vibration, strain etc.) on whole optical fiber link may be implemented
Detection.Phase sensitive optical time domain reflectometer (Φ-OTDR) has as a kind of Distributed Optical Fiber Sensing Techniques of high sensitivity with it
There is dynamic monitoring performance and is widely studied.It can be realized vehicle location and tracking, iron using phase sensitive optical time domain reflectometer
Road safety detection, circumference security protection etc., with good application prospect.
Currently, the research of the Distributed Optical Fiber Sensing Techniques based on phase sensitive optical time domain reflectometer is mostly with direct detection side
Based on formula and intensity demodulation.Time difference usually is carried out to scattered light intensity curve, to determine the position of disturbance point.According to theory
Analysis, it is not linear relationship, even that the relationship between scattered light intensity (or photoelectricity field amplitude) and external disturbance is complex
Without monotonicity.Therefore, this mode can only carry out disturbance positioning and qualitative detection, it is difficult to realize quantitative measurment.
In addition, in the prior art, the output signal of phase sensitive optical time domain reflectometer generally all uses Data Acquisition Card
It is digitized, Data Acquisition Card is sent to host computer and is calculated and analyzed after acquiring data.It is being using this scheme
The design phase unite with preferable flexibility, change data processing algorithm can be facilitated at any time, but in practical applications, due to
The data volume that phase sensitive optical time domain reflectometer generates quicklys increase, this scheme with the increase of distance sensing and sensing speed
It is limited to the software processing speed of buffer memory capacity on the plate of capture card, capture card interface data transmission speed and host computer, nothing
Method is realized and is sensed over long distances with the optical fiber disturbance of high speed.
Summary of the invention
To solve the above problems, the purpose of the present invention is to provide a kind of Real-Time Optical fiber disturbance sensings of highly-parallel
Quantitatively measuring device and its method, and carry out hardware realization.
To achieve the above object, the technical scheme is that a kind of method for quantitative measuring of Real-Time Optical fiber disturbance sensing and
It realizes that hardware composition includes phase sensitive optical time domain reflectometer H1；Double-channel analog/digital converter H2；FPGA module H3；
DDR3 random access memory H4；Industrial control computer H5 and GPU H6, it is characterised in that: the phase sensitive optical time domain reflectometer H1's
Two-way analog signal connects double-channel analog/digital converter H2, the double-channel analog/digital converter H2 connection FPGA module H3, described
DDR3 random access memory H4 connection FPGA module H3, the FPGA module H3 connection industrial control computer H5, the GPUH6 connection work
Control computer H5.
Wherein, the phase sensitive optical time domain reflectometer is connected with double-channel analog/digital converter, for acquiring I, Q two-way
Analog signal；The FPGA module is connected with double-channel analog/digital converter, for carrying out data receiver and calculating short-time energy
Data Concurrent gives industrial control computer；The DDR3 random access memory is connected with FPGA module, for data cached；The work
Control computer is connected with FPGA module, is sent to GPU on the upper computer interface and simultaneously for showing the data received；
The GPU is connected with industrial control computer, for carrying out pattern-recognition and sending result to industrial control computer；
The phase sensitive optical time domain reflectometer H1, including the following contents:
Laser light source H10: laser light source is provided for the sensor-based system.
Sequence controller H11: for controlling the operation of the phase sensitive optical time domain reflectometer.
Photo-coupler H12: the coupling for multi-path laser.
Acousto-optic modulator H13: for continuous laser to be converted to pulsed light.
Erbium-doped fiber amplifier H14: for improving laser emitting power.
Circulator H15: the direction controlling for laser.
Sensor fibre H16: for incuding disturbance.
Double balanced detector H17: for converting optical signals to electric signal.
Quadrature demodulator H18: for converting electrical signals to orthogonal signalling.
Local local oscillation signal H19: local oscillation signal is provided for quadrature demodulation.
Its operation logic are as follows: laser light source H10 exports continuous light, connects photo-coupler H12A by optical fiber and divides laser
For two-way: being used as reference light all the way；Acousto-optic modulator H13 is connected by optical fiber all the way, the copped wave under sequence controller H11 control
As pulse laser, and generate upshift.Acousto-optic modulator connects erbium-doped fiber amplifier H14 by optical fiber, and pulsed light is put
Greatly, erbium-doped fiber amplifier output end connects the A mouth of circulator H15 by optical fiber, and the B mouth of circulator H15 connects sensor fibre
H16, laser enter sensor fibre H16 by the B mouth of circulator H15.The backward Rayleigh scattering of direct impulse is returned along sensor fibre
It returns, into the B mouth of circulator, is emitted through C mouthfuls.Two input terminal of photo-coupler H12B is separately connected reference light and scattering light, scattering
Light and reference light respectively enter two input terminals of photo-coupler H12B, carry out beat frequency.Photo-coupler H12B output end connects double flat
Weigh detector H17, and beat frequency output is converted to electric signal by double balanced detector H17.The output end of double balanced detectors and local
Electric signal and local local oscillation signal H19 are inputted quadrature demodulator H18, output by local oscillation signal H19 connection quadrature demodulator H18
Two-way I, Q analog signal.Phase sensitive optical time domain reflectometer is controlled it and run and export commutator pulse by sequence controller H11 to be believed
Number for follow-up signal processing commutator pulse is provided.
It is found on the basis of theoretical and practical studies according to applicant, the space differential phase for scattering light is disturbed with the external world
Move linear, external disturbance will result directly in the linear change of space differential phase, therefore scatter the phase information of light
It can be well reflected disturbing signal.The present invention uses heterodyne detection mode, and optocoupler is arranged in phase sensitive optical time domain reflectometer
Clutch H12B, double balanced detector H17, quadrature demodulator H18 and local local oscillation signal H19, scattering light and local light are carried out
The optical signal received is converted to electric signal by double balanced detectors (BPD), and electric signal is used quadrature demodulator by beat frequency
It is demodulated into I, Q orthogonal signalling.Orthogonal I/Q demodulation is a kind of common method in the demodulation of radio frequency (RF) signal, I, Q letter of output
Number be respectively scattering light EbackThe sin cos functions of phase.
Δ φ=arctan (I/Q)+2m π, m=..., -1,0,1 ...
Therefore the amplitude and phase information that scattering light can be demodulated in numeric field by orthogonal I/Q demodulation method, thus real
Now to the quantitative detection of strength of turbulence.
Scheme as a further improvement, the FPGA module H3, inside be divided into three levels: hardware alternation of bed
L1, data DMA layers of L2 and algorithm process layer L3.Wherein the hardware alternation of bed L1 includes hardware control module: LVDS receiver
H20, DDR3 Memory Controller Hub H21, PCI Express controller H22；DMA layers of L2 of data include data DMA transfer logic electricity
Road: AXI interconnected H23, AXI reading-writing port H24, PCIe DMA drives H25；Algorithm process layer L3 includes data process method
Circuit: adc data unit H26, data buffer storage unit H27, multichannel FFT computing unit H28, data filtering unit H29, in short-term
Energy calculation unit H30.
Scheme as a further improvement, FPGA module execute following steps:
Step S1: after the system brings into operation, the sequence controller H11 in phase sensitive optical time domain reflectometer H1 starts
The sampling pulse of certain frequency is generated, and exports Quadrature analog signal I, Q.
Step S2: raw data acquisition and pretreatment, comprising the following steps:
Step S21: two-way Quadrature analog signal I, Q are converted to digital signal by the double-channel analog/digital converter H2, raw
It embarks on journey data.Its sample rate is related to sensing resolution, and data volume is related to distance sensing.
Step S22: the FPGA module H3 acquires adc data, and calculates the modulus value of the quadrature digital signalAnd each row of data is deposited by first of DDR3 Memory Controller Hub H21 write-in DDR3 random access memory H4
Storage area domain, multiple data lines form a data frame.
Step S3: data transposition and calculating, comprising the following steps:
Step S31: the FPGA module H3 arranges multiple column to read data frame by DDR3 Memory Controller Hub H21, and
It is stored in the data buffer storage unit H27 of FPGA module H3.
Step S32: the FPGA module H3 sends multichannel FFT to data for the column in its data buffer storage unit H27
Computing unit H28 does FFT calculating.
Step S33: repeating step S32, calculates to form one until the FPGA module H3 carries out FFT to a data frame
FFT result matrix, that is, the disturbing spectrum information of each sensing point.
Step S4: data filtering and summation, comprising the following steps:
Step S41: the data filtering unit H29 in the FPGA module H3 carries out the matrix of consequence in the step S3
Laterally summation, and the maximum value in summed result is found, while matrix of consequence being stored to the second of DDR3 random access memory H4
A storage region.
Step S42: the maximum value place that the data filtering unit H29 in the FPGA module H3 is obtained according to step S41
Row, is filtered the matrix of consequence in step S3.
After step S43: short-time energy computing unit H30 in the FPGA module H3 filtering obtained according to step S42
FFT result matrix carry out short-time energy calculating, obtain the strength of turbulence data of each sensing point, and short-time energy data are deposited
Store up the third storage region of DDR3 random access memory H4.
Step S5: step S4 is obtained short-time energy data and drives H25 to be sent to by PCIe DMA by the FPGA module
Industrial control computer H5 host computer.
Step S6: the industrial control computer H5 host computer shows the data received on the upper computer interface, simultaneously will
Data, which are sent on GPU H6, carries out pattern-recognition, and pattern recognition result returns to host computer interface and shows.
In the above-mentioned technical solutions, short-time energy data are obtained by calculating due to FPGA module H3 of the present invention, in short-term can
Meter algorithm carries out frequency domain analysis to each sensing location: being filtered first in frequency domain to sensing data by filtering algorithm
Then wave adds up to spectrum energy.Compared to traditional time difference method, this method is capable of inherently making an uproar for filtering appts
Sound, while realizing the quantitative calculating to strength of turbulence.
Scheme as a further improvement, the realization step are run using Fully-pipelinedization.
Scheme as a further improvement, the FPGA module use heterodyne detection mode and the short-time energy meter to disturbance
Calculate the quantitative demodulation realized to strength of turbulence.
Scheme as a further improvement, the data frame in the FPGA module use sliding window mode, and data are
The combination of freshly harvested data and the partial data of previous frame.
Scheme as a further improvement, the pattern-recognition disturbed using GPU using CNN disaggregated model with point
Class.
Scheme as a further improvement, the address space of the DDR3 random access memory are divided into three regions, respectively
Store the different data for calculating step.
Compared with prior art, the present invention at least has the following beneficial effects:
(1) currently, the research of the Distributed Optical Fiber Sensing Techniques based on phase sensitive optical time domain reflectometer is mostly directly to visit
Based on survey mode and intensity demodulation.But the relationship between scattered light intensity (or photoelectricity field amplitude) and external disturbance is more multiple
Miscellaneous, this mode can only carry out disturbance positioning and qualitative detection, it is difficult to realize quantitative measurment.It is quick that the present invention has redesigned phase
Photosensitive domain reflectometer is led to using heterodyne detection mode using the phase information of scattering light to be well reflected disturbing signal
The amplitude and phase information that orthogonal I/Q demodulation method demodulates scattering light in numeric field is crossed, strength of turbulence is determined to realize
Amount detection.
(2) in the prior art generally all data are acquired using Data Acquisition Card after be sent to host computer carry out calculate with
Analysis.Cache size and the software processing speed of host computer on the limited plate with capture card of this scheme, cannot achieve it is long away from
It is sensed from the optical fiber disturbance of high speed.The present invention realizes hardware concurrent and the pipelining processing of algorithm using FPGA, from
And it realizes and is sensed over long distances with the optical fiber disturbance of high speed.
(3) its data frame is typically independent in the prior art, and front and back does not have relevance, will will lead to so in short-term
Between disturbance can only show a frame on a monitor, high refresh rate display when be unfavorable for observing.The present invention uses sliding window
Mode sampling, data frame not only includes newest sampled data, further includes a part of history samples data, can be realized disturbance
Transition show, improve the readability of result；
(4) present invention carries out pattern-recognition and classification, convolution to disturbance by convolutional neural networks after calculating data
Neural network is in speech recognition, recognition of face, generic object identification, motion analysis, natural language processing even brain wave analysis
Aspect has excellent performance, can carry out accurate modality identification and classification to disturbance.
Detailed description of the invention
Fig. 1 is the system construction drawing of the quantitatively measuring device of the Real-Time Optical fiber disturbance sensing of the embodiment of the present invention；
Fig. 2 is the structure chart of FPGA module in the embodiment of the present invention；
Fig. 3 is the software steps flow chart of the method for quantitative measuring of the Real-Time Optical fiber disturbance sensing of the embodiment of the present invention；
Specific embodiment
In order to make the objectives, technical solutions, and advantages of the present invention clearer, with reference to the accompanying drawings and embodiments, right
The present invention is further elaborated.It should be appreciated that the specific embodiments described herein are merely illustrative of the present invention, and
It is not used in the restriction present invention.
On the contrary, the present invention covers any substitution done on the essence and scope of the present invention being defined by the claims, repairs
Change, equivalent method and scheme.Further, in order to make the public have a better understanding the present invention, below to of the invention thin
It is detailed to describe some specific detail sections in section description.Part without these details for a person skilled in the art
The present invention can also be understood completely in description.
In the prior art, the research of the Distributed Optical Fiber Sensing Techniques based on phase sensitive optical time domain reflectometer is mostly with direct
Based on detection mode and intensity demodulation.Time difference usually is carried out to scattered light intensity curve, to determine the position of disturbance point.Root
According to theory analysis, it is not linearly to close that the relationship between scattered light intensity (or photoelectricity field amplitude) and external disturbance is complex
System, or even do not have monotonicity.Therefore, this mode can only carry out disturbance positioning and qualitative detection, it is difficult to realize quantitative measurment.
In view of the drawbacks of the prior art, applicant has found on the basis of theoretical and practical studies, scatters the space parallax of light
Split-phase position be it is linear with external disturbance, external disturbance will result directly in the linear change of space differential phase.Therefore
The phase information of scattering light can be well reflected disturbing signal, it is necessary to demodulate to the phase information of scattering light.It is true
The phase information for protecting scattering light can remain, and generally use heterodyne detection mode, i.e. scattering light and local light carries out beat frequency,
It is received by double balanced detectors (BPD).Orthogonal I/Q demodulation is a kind of common method in the demodulation of radio frequency (RF) signal, I, Q difference
To scatter light EbackThe sin cos functions of phase.
Δ φ=arctan (I/Q)+2m π, m=..., -1,0,1 ...
The amplitude and phase information of scattering light can be demodulated by orthogonal I/Q demodulation method, to realize to strength of turbulence
Quantitative detection.
Using above-mentioned principle, the present invention proposes a kind of quantitatively measuring device of Real-Time Optical fiber disturbance sensing, hardware composition
Including phase sensitive optical time domain reflectometer H1；Double-channel analog/digital converter H2；FPGA module H3；DDR3 random access memory H4；Work
Control computer H5 and GPU H6, it is characterised in that: the two-way analog signal connection of the phase sensitive optical time domain reflectometer H1 is double
Channel modulus converter H2, the double-channel analog/digital converter H2 connection FPGA module H3, the DDR3 random access memory H4 company
Meet FPGA module H3, the FPGA module H3 connection industrial control computer H5, the GPU H6 connection industrial control computer H5.
The phase sensitive optical time domain reflectometer H1, including the following contents:
Laser light source H10: laser light source is provided for the sensor-based system.
Sequence controller H11: for controlling the operation of the phase sensitive optical time domain reflectometer and FPGA module.
Photo-coupler H12: the coupling for multi-path laser.
Acousto-optic modulator H13: for continuous laser to be converted to pulsed light.
Erbium-doped fiber amplifier H14: for improving laser emitting power.
Circulator H15: the direction controlling for laser.
Sensor fibre H16: for incuding disturbance.
Double balanced detector H17: for converting optical signals to electric signal.
Quadrature demodulator H18: for converting electrical signals to orthogonal signalling.
Local local oscillation signal H19: local oscillation signal is provided for quadrature demodulation.
Its operation logic are as follows: laser light source H10 exports continuous light, connects photo-coupler H12A by optical fiber and divides laser
For two-way: being used as reference light all the way；Acousto-optic modulator H13 is connected by optical fiber all the way, the copped wave under sequence controller H11 control
As pulse laser, and generate upshift.Acousto-optic modulator connects erbium-doped fiber amplifier H14 by optical fiber, and pulsed light is put
Greatly, erbium-doped fiber amplifier output end connects the A mouth of circulator H15 by optical fiber, and the B mouth of circulator H15 connects sensor fibre
H16, laser enter sensor fibre H16 by the B mouth of circulator H15.The backward Rayleigh scattering of direct impulse is returned along sensor fibre
It returns, into the B mouth of circulator, is emitted through C mouthfuls.Two input terminal of photo-coupler H12B is separately connected reference light and scattering light, scattering
Light and reference light respectively enter two input terminals of photo-coupler H12B, carry out beat frequency.Photo-coupler H12B output end connects double flat
Weigh detector H17, and beat frequency output is converted to electric signal by double balanced detector H17.The output end of double balanced detectors and local
Electric signal and local local oscillation signal H19 are inputted quadrature demodulator H18, output by local oscillation signal H19 connection quadrature demodulator H18
Two-way I Q analog signal.Phase sensitive optical time domain reflectometer is controlled it and run and export commutator pulse by sequence controller H11 to be believed
Number for follow-up signal processing commutator pulse is provided.
Quadrature analog signal I, Q of phase sensitive optical time domain reflectometer H1 output are respectively the sine and cosine letter for scattering light phase
Number.The amplitude and phase information of scattering light can be demodulated in numeric field by orthogonal I/Q demodulation method, to realize strong to disturbing
The quantitative detection of degree.In the present embodiment, the double-channel analog/digital converter H2 realizes the number to Quadrature analog signal I, Q
Change, the FPGA module realizes that the magnitude calculation to I, Q digital signal, FFT calculate, short-time energy calculates, and realizes
To the quantitative detection of strength of turbulence.
Referring to fig. 2, show the theory structure block diagram of FPGA module H3, ask inside it and be divided into three levels: hardware is handed over
Alternating layers L1, data DMA layers of L2 and algorithm process layer L3.Wherein the hardware alternation of bed L1 includes hardware control module: LVDS connects
Receive device H20, DDR3 Memory Controller Hub H21, PCI Express controller H22；DMA layers of L2 of data include data DMA transfer logic
Circuit: AXI interconnected H23, AXI reading-writing port H24, PCIe DMA drives H25；Algorithm process layer L3 is patrolled comprising data processing
Volume circuit: adc data unit H26, data buffer storage unit H27, multichannel FFT computing unit H28, data filtering unit H29, short
When energy calculation unit H30.
The industrial control computer has multiple PCIe connectors, is equipped with FPGA module H3 and GPU H6, is also equipped with general
Processor essential to computer, display equipment, storage equipment and input-output equipment etc..
The address space of the DDR3 random access memory is divided into three regions, stores the different numbers for calculating step respectively
According to.In the present embodiment the capacity of DDR3 random access memory be 8Gbit, three storage regions store respectively IQ modulus value data,
FFT result data and short-time energy data.
By adopting the above technical scheme, the present invention uses the mode of heterodyne detection to the phase of phase sensitive optical time domain reflectometer
It is demodulated, realizes the quantitative measurment to strength of turbulence；Its data processing method uses sliding window mode, can be realized and disturbs
Dynamic transition is shown, improves the readability of result；It is deposited at random using novel fpga chip large capacity, high concurrent speed and DDR3
The characteristics of reservoir high bandwidth, high capacity, realizes the multi-level buffer of data, and multi-channel parallel calculates and pipeline processes, solves
The problem of optical fiber disturbance sensor-based system cannot be realized simultaneously over long distances and at high speed, while using GPU to disturbance progress convolution
Neural network filter and classification, have preferable real-time and accuracy.
Referring to Fig. 3, the invention also provides a kind of method for quantitative measuring of Real-Time Optical fiber disturbance sensing, specific at this
Its operational process in embodiment the following steps are included:
Step S1: after the system brings into operation, the laser light source H10 in phase sensitive optical time domain reflectometer H1, which is exported, to be connected
Laser is divided into two-way by continuous light, photo-coupler H12A: all the way as reference light；Become all the way by acousto-optic modulator H13 copped wave
Pulse laser, and generate upshift.Erbium-doped fiber amplifier H14 will input the A mouth of circulator H15, laser after pulse light amplification
Enter sensor fibre H16 by the B mouth of circulator H15.The backward Rayleigh scattering light of direct impulse is returned along sensor fibre, is entered
The B mouth of circulator is emitted through C mouthfuls.Scattering light and reference light respectively enter two input terminals of photo-coupler H12B, carry out beat frequency.
Beat frequency output light is converted to electric signal by double balanced detector H17.Electric signal and local local oscillation signal H19 are inputted into quadrature demodulation
Device H18, demodulation output two-way I, Q analog signal.In the present embodiment, the width of pulsed light is taken as 10ns, corresponding space point
Resolution is 1m；
Sequence controller H11 starts to generate the pulse f of certain frequency.The pulse frequency f of sequence controller H11, which is determined, is
The maximum perturbation frequency that system can detect, such as in the present embodiment, f=2kHz, according to nyquist sampling law, system
The maximum perturbation frequency that can be detected is 1kHz.
Step S2: raw data acquisition and pretreatment, comprising the following steps:
Step S21: phase sensitive optical time domain reflectometer H1 sampling pulse signal controls the double-channel analog/digital and turns
Parallel operation H2 acquires two-way I, Q analog signal, M data of each pulse collection, as data line by particular sample speed S.
ADC sample rate is related to sensing resolution, and each row of data amount is related to distance sensing, in the present embodiment, ADC
Sample rate S=100Msps, the sensing space resolution ratio of corresponding 1m, M=40000 data of each pulse collection, to inductive sensing away from
From L=40000m.
Step S22: the FPGA module H3 acquires adc data, and calculates the modulus value of the quadrature digital signalAnd each row of data is deposited by first of DDR3 Memory Controller Hub H21 write-in DDR3 random access memory H4
Storage area domain, multiple data lines form a data frame.
In the present embodiment, ADC resolution ratio R=14bit, DDR3 random access memory are divided into three storage regions, and first
Storage region stores ADC and acquires data, and sets fifo mode for this storage region, and FIFO depth is 2000 rows.
In the present embodiment, data frame is made of 500 row adc datas, and component content is divided to two situations:
Situation 1: in the case that system is sky in DDR3 memory at the beginning of powering on, data frame is by 500 initial row data
Composition.
Situation 2: system preserves preamble frame data, data frame is by 400 row of preamble in even running in DDR3 memory
Data and newest 100 row data composition, thereby realize sliding window calculating, and calculated result can show data transition variation
Effect.
Step S3: data transposition and calculating, comprising the following steps:
Step S31: the FPGA module H3 arranges multiple column to read data frame by DDR3 Memory Controller Hub H21, and
It is stored in the data buffer storage unit H27 of FPGA module H3.In the present embodiment, DDR3 random access memory is compared to static state
The read-write delay of SRAM is big, therefore in order to improve reading efficiency, this embodiment introduces a FPGA on piece static state RAM conducts
The caching of DDR3 memory once reads more data buffer storage to FPGA on piece static state RAM, in this way may be used when reading data
To save the more read latch time, reading efficiency is improved.
Step S32: the FPGA module H3 sends multichannel FFT to data for the column in its data buffer storage unit H27
Computing unit H28 does FFT calculating.In the present embodiment, 512 point FFT computing modules of 8 full flowing water are shared while carrying out FFT
It calculates, Quan Liushui FFT computing module can not stop persistently input data, and continue output and calculate the data completed.
Step S33: repeating step S32, calculates to form one until the FPGA module H3 carries out FFT to a data frame
FFT result matrix.In the present embodiment, the size A*B=40000*512 of FFT result matrix indicates at 40000 groups 512 points
FFT data, physical significance are the disturbing spectrum data of the measurement point every 1m of 40KM optical fiber.
Step S4: data filtering and summation, comprising the following steps:
Step S41: the data filtering unit H29 in the FPGA module H3 carries out the matrix of consequence in the step S3
Laterally summation, and the maximum value in summed result is found, while matrix of consequence being stored to the second of DDR3 random access memory H4
A storage region.
In the present embodiment, it arranges according to every row X 1 to the column of X 2 and sums, obtain 512 values, and find out the 52nd and arrive
Maximum value in 250th number (X1, X2 needs are modified according to application scenarios, such as X1=1000, X2=20000).
Step S42: the maximum value place that the data filtering unit H30 in the FPGA module H3 is obtained according to step S41
Row, is filtered the matrix of consequence in step S3.
In the present embodiment, according to the row ZB where the maximum value acquired in step S41, by matrix of consequence in step S3
(ZB-10)~(ZB+10) value all clear 0 (if ZB-10 or ZB+10 reach boundary 52 or 250, then can only get boundary
Value).
After step S43: short-time energy computing unit H31 in the FPGA module H3 filtering obtained according to step S42
FFT result matrix carry out short-time energy calculating, short-time energy obtains from FFT result rectangular array to summation, and by result data
Store the third storage region of DDR3 random access memory H4.In the present embodiment, there are 40000 result datas, physics
Meaning is the strength of turbulence of each sensing point, distance sensing 40KM, sensing resolution 1m.
Step S5: step S4 is obtained short-time energy data and drives H25 to be sent to by PCIe DMA by the FPGA module
Industrial control computer H5 host computer.
Step S6: the industrial control computer H5 host computer shows the data received on the upper computer interface, simultaneously will
Data are sent to GPU H6, and GPU H6 carries out pattern-recognition according to strength of turbulence.
Repeat above-mentioned steps, until artificial command program end of run.
The foregoing is merely illustrative of the preferred embodiments of the present invention, is not intended to limit the invention, all in essence of the invention
Made any modifications, equivalent replacements, and improvements etc., should all be included in the protection scope of the present invention within mind and principle.
Claims (9)
1. a kind of quantitatively measuring device of Real-Time Optical fiber disturbance sensing, which is characterized in that including phase sensitive optical time domain reflectometer,
Double-channel analog/digital converter, FPGA module, DDR3 random access memory, industrial control computer and GPU, wherein
The phase sensitive optical time domain reflectometer is connected with double-channel analog/digital converter, for acquiring I, Q two-way analog signal；
The FPGA module is connected with double-channel analog/digital converter, send for carrying out data receiver and calculating short-time energy Data Concurrent
To industrial control computer；The DDR3 random access memory is connected with FPGA module, for data cached；The industrial control computer with
FPGA module is connected, and is sent to GPU on the upper computer interface and simultaneously for showing the data received；The GPU with
Industrial control computer is connected, for carrying out pattern-recognition and sending result to industrial control computer；
The phase sensitive optical time domain reflectometer further comprises: laser light source, sequence controller, photo-coupler, acousto-optic modulation
Device, erbium-doped fiber amplifier, circulator, sensor fibre, double balanced detectors, quadrature demodulator, local local oscillation signal, wherein
The laser light source connects photo-coupler by optical fiber, and laser is divided into two-way: all the way as reference light, another way passes through optical fiber
It is sequentially connected acousto-optic modulator, erbium-doped fiber amplifier, A mouthfuls of circulator and sensor fibre, scattering light returns along B mouthfuls, goes out through C mouthfuls
It penetrates；Two input terminal of photo-coupler is separately connected reference light and scattering light, and output end connects double balanced detectors；Double balanced detectors
Output end connect quadrature demodulator, quadrature demodulator output two-way I, Q analog signal with local local oscillation signal；Phase sensitive light
Domain reflectometer by sequence controller control its run and for follow-up signal processing commutator pulse is provided.
2. the quantitatively measuring device of Real-Time Optical fiber disturbance sensing according to claim 1, which is characterized in that the FPGA mould
Block, inside be divided into three levels: hardware alternation of bed, DMA layers of data and algorithm process layer；The hardware alternation of bed includes
LVDS receiver, DDR3 Memory Controller Hub, PCI Express controller；Described data DMA layers include AXI interconnected, AXI
Reading-writing port, PCIe DMA driving；The algorithm process layer includes adc data unit, data buffer storage unit, multichannel FFT meter
Calculate unit, data filtering unit, short-time energy computing unit.
3. a kind of method for quantitative measuring of Real-Time Optical fiber disturbance sensing, which comprises the following steps:
Step S1: the sequence controller in phase sensitive optical time domain reflectometer starts to generate the sampling pulse of certain frequency, and defeated
Quadrature analog signal I, Q out；
Step S2: raw data acquisition and pretreatment, comprising the following steps:
Step S21: two-way Quadrature analog signal I, Q are converted to digital signal and are sent to FPGA mould by double-channel analog/digital converter
Block；
The modulus value of step S22:FPGA module calculating quadrature digital signalAnd each row of data is passed through in DDR3
First storage region of DDR3 random access memory is written in memory controller, and multiple data lines form a data frame；
Step S3: data transposition and calculating, comprising the following steps:
DDR3 Memory Controller Hub in step S31:FPGA module arranges multiple column to read data frame, and is stored to described
In the data buffer storage unit of FPGA module；
Step S32:FPGA module sends multichannel FFT computing unit to data for the column in its data buffer storage unit and is FFT
It calculates；
Step S33: repeating step S32, calculates to form a FFT knot until the FPGA module carries out FFT to a data frame
Fruit matrix；
Step S4: data filtering and summation, comprising the following steps:
Data filtering unit in step S41:FPGA module carries out the specific region of the matrix of consequence in the step S3 horizontal
To summation, and the maximum value in summed result specific region is found, while matrix of consequence being stored to DDR3 random access memory
Second storage region；
Step S42: the data filtering unit in the FPGA module is expert at according to the maximum value that step S41 is obtained, to step
Matrix of consequence in S3 is filtered；
The filtered FFT result matrix that short-time energy computing unit in step S43:FPGA module is obtained according to step S42
Short-time energy calculating is carried out, obtains the strength of turbulence data of each sensing point, and the storage of short-time energy data is random to DDR3
The third storage region of memory；
Step S4 is obtained short-time energy data and is sent in industrial control computer by PCIe DMA driving by step S5:FPGA module
Position machine；
Step S6: industrial control computer host computer shows the data received on the upper computer interface, transmits data to simultaneously
Pattern-recognition is carried out on GPU, pattern recognition result returns to host computer interface and shows.
4. implementation method according to claim 3, which is characterized in that in the step S6, disturbed using GPU
Pattern recognition and classification further includes steps of
Step S61: the data set of input disturbance, training CNN disaggregated model；
Step S62: the data of disturbance to be identified are inputted into trained model and are identified.
5. implementation method according to claim 3, which is characterized in that above-mentioned steps are run using Fully-pipelinedization.
6. implementation method according to claim 3, which is characterized in that FPGA module is using heterodyne detection mode and to disturbance
Short-time energy calculate and realize quantitative demodulation to strength of turbulence.
7. implementation method according to claim 3, which is characterized in that the data frame in the FPGA module uses sliding window
Mouth mold formula, data frame are the combination of the partial data of freshly harvested data and previous frame.
8. implementation method according to claim 3, which is characterized in that the address space of the DDR3 random access memory divides
For three regions, the different data for calculating step are stored respectively.
9. implementation method according to claim 3, which is characterized in that the FPGA module receives modulus by LVDS interface
ADC data.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910191993.8A CN110044397A (en) | 2019-03-14 | 2019-03-14 | A kind of quantitatively measuring device and its method of Real-Time Optical fiber disturbance sensing |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910191993.8A CN110044397A (en) | 2019-03-14 | 2019-03-14 | A kind of quantitatively measuring device and its method of Real-Time Optical fiber disturbance sensing |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN110044397A true CN110044397A (en) | 2019-07-23 |
Family
ID=67273694
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910191993.8A Pending CN110044397A (en) | 2019-03-14 | 2019-03-14 | A kind of quantitatively measuring device and its method of Real-Time Optical fiber disturbance sensing |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110044397A (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112710379A (en) * | 2019-10-25 | 2021-04-27 | 中兴通讯股份有限公司 | Phase acquisition method and device for optical fiber disturbance event and storage medium |
| CN115355931A (en) * | 2022-08-25 | 2022-11-18 | 中国科学院半导体研究所 | Distributed disturbance measurement method and system based on sensing optical fiber |
| CN115355931B (en) * | 2022-08-25 | 2024-09-27 | 中国科学院半导体研究所 | Distributed disturbance measurement method and disturbance measurement system based on sensing optical fiber |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102645268A (en) * | 2012-04-26 | 2012-08-22 | 中国科学院上海光学精密机械研究所 | Optical frequency division multiplexing phase-sensitive optical time domain reflectometer |
| CN103900623A (en) * | 2014-04-16 | 2014-07-02 | 武汉理工光科股份有限公司 | Optical time domain reflectometer based on double acoustic-optical modulators and common-mode rejection method of optical time domain reflectometer |
| CN105466548A (en) * | 2015-12-16 | 2016-04-06 | 上海大学 | Phase sensitive optical time domain reflection fiber sensing system positioning method |
| WO2018076551A1 (en) * | 2016-10-25 | 2018-05-03 | 北京大学 | Distributed optical fibre sensing system capable of realizing all-phase demodulation, and measuring method therefor |
| CN108180978A (en) * | 2018-01-30 | 2018-06-19 | 广州天长通信技术有限公司 | A kind of combination PGC technologies and the method and device of Φ-OTDR technique detection optical fiber vibration |
| CN108225541A (en) * | 2017-12-29 | 2018-06-29 | 鞍山睿科光电技术有限公司 | The distributed fiberoptic sensor and foreign body intrusion signal for identifying foreign body intrusion perceive processing method |
| CN108489598A (en) * | 2018-05-21 | 2018-09-04 | 浙江杰昆科技有限公司 | Phase sensitive optical time domain reflectometer phase simulates coherent demodulation system and method |
-
2019
- 2019-03-14 CN CN201910191993.8A patent/CN110044397A/en active Pending
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102645268A (en) * | 2012-04-26 | 2012-08-22 | 中国科学院上海光学精密机械研究所 | Optical frequency division multiplexing phase-sensitive optical time domain reflectometer |
| CN103900623A (en) * | 2014-04-16 | 2014-07-02 | 武汉理工光科股份有限公司 | Optical time domain reflectometer based on double acoustic-optical modulators and common-mode rejection method of optical time domain reflectometer |
| CN105466548A (en) * | 2015-12-16 | 2016-04-06 | 上海大学 | Phase sensitive optical time domain reflection fiber sensing system positioning method |
| WO2018076551A1 (en) * | 2016-10-25 | 2018-05-03 | 北京大学 | Distributed optical fibre sensing system capable of realizing all-phase demodulation, and measuring method therefor |
| CN108225541A (en) * | 2017-12-29 | 2018-06-29 | 鞍山睿科光电技术有限公司 | The distributed fiberoptic sensor and foreign body intrusion signal for identifying foreign body intrusion perceive processing method |
| CN108180978A (en) * | 2018-01-30 | 2018-06-19 | 广州天长通信技术有限公司 | A kind of combination PGC technologies and the method and device of Φ-OTDR technique detection optical fiber vibration |
| CN108489598A (en) * | 2018-05-21 | 2018-09-04 | 浙江杰昆科技有限公司 | Phase sensitive optical time domain reflectometer phase simulates coherent demodulation system and method |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112710379A (en) * | 2019-10-25 | 2021-04-27 | 中兴通讯股份有限公司 | Phase acquisition method and device for optical fiber disturbance event and storage medium |
| CN112710379B (en) * | 2019-10-25 | 2023-08-01 | 中兴通讯股份有限公司 | Method, device and storage medium for acquiring phase of optical fiber disturbance event |
| CN115355931A (en) * | 2022-08-25 | 2022-11-18 | 中国科学院半导体研究所 | Distributed disturbance measurement method and system based on sensing optical fiber |
| CN115355931B (en) * | 2022-08-25 | 2024-09-27 | 中国科学院半导体研究所 | Distributed disturbance measurement method and disturbance measurement system based on sensing optical fiber |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102262246B (en) | Method and apparatus for marine seismic data acquisition | |
| CN103196465B (en) | Method for separating noise of sensing signal and extracting signal of phase-sensitive optical time domain reflectometer | |
| CN104808208B (en) | Measurement system and method for detecting orientation and size of underwater target on basis of laser sound source | |
| Johnston et al. | Internal tidal beams and mixing near Monterey Bay | |
| CN112883564B (en) | Water body temperature prediction method and prediction system based on random forest | |
| CN103226028B (en) | Method for identifying and detecting disturbance signals of phase-sensitive optical time domain reflectometer | |
| CN107179097A (en) | The optical fiber interference type sensing modulation-demo-demodulation method and device modulated based on phase shift signal | |
| CN107782438A (en) | Optical fiber raster vibration measuring system based on pulse code and edge filter demodulation method | |
| CN108415067A (en) | A kind of earthquake wave measuring system based on microstructured optical fibers distribution sound wave sensing | |
| CN103499383A (en) | Self-correlation positioning method for improving positioning accuracy of optical fiber vibration sensor | |
| CN103438980B (en) | Method and device for liquid surface wave detection based on linear array CCD and linear infrared laser | |
| CN101603817B (en) | Glass thickness detection device and method | |
| Bai et al. | Noise reduction method of Φ-OTDR system based on EMD-TFPF algorithm | |
| CN110044397A (en) | A kind of quantitatively measuring device and its method of Real-Time Optical fiber disturbance sensing | |
| CN109752720A (en) | Based on distributed vibration sensor vibration information localization method and system | |
| Lin et al. | Monitoring ocean currents during the passage of Typhoon Muifa using optical-fiber distributed acoustic sensing | |
| CN204228605U (en) | A kind of measurement mechanism of water turbidity and measuring system | |
| Cornuelle | Acoustic tomography | |
| CN107014409A (en) | A kind of long range optical frequency domain reflection-based optical fiber Distributed Multi destabilization sensing method | |
| CN116026449A (en) | Vibration positioning monitoring system based on single-core optical fiber sensing | |
| CN104833657A (en) | Laser wireless sand meter having side-direction compensation | |
| CN115752696A (en) | Distributed optical fiber vibration sensing system based on forward transmission and signal processing method | |
| CN102193082A (en) | Device for positioning leak source of three-sensor multi-scale constrained pipe network | |
| Kozick et al. | Algorithms for tracking with an array of magnetic sensors | |
| CN105572742A (en) | Method and device for determining depth of seawater |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190723 |
|
| RJ01 | Rejection of invention patent application after publication |