CN105241390B - Quick Brillouin optical time domain analysis type strain gauge means and data processing method - Google Patents
Quick Brillouin optical time domain analysis type strain gauge means and data processing method Download PDFInfo
- Publication number
- CN105241390B CN105241390B CN201510686227.0A CN201510686227A CN105241390B CN 105241390 B CN105241390 B CN 105241390B CN 201510686227 A CN201510686227 A CN 201510686227A CN 105241390 B CN105241390 B CN 105241390B
- Authority
- CN
- China
- Prior art keywords
- signal
- frequency
- modulator
- optical
- output
- 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.)
- Expired - Fee Related
Links
Abstract
A kind of high-performance fast B OTDA (Brillouin optical time domain analysis type) strain gauge means being combined based on the processing of excited Brillouin gain spectrum width, pulse code and wavelet transformation technique and data processing method, belong to Distributed Optical Fiber Sensing Techniques field.By laser 1, first photo-coupler 2, first modulator 3, first microwave signal source 20, first optical filter 4, second modulator 5, second microwave signal source 21, DC power supply 22, image intensifer 6, second photo-coupler 7, optoisolator 8, second modulator 9, frequency synthesizer 14, second optical filter 10, 4th modulator 11, pulse signal generator 15, scrambler 12, optical circulator 13, sensor fibre 19, photodetector 16, data collecting card 17 and computer 18 form, the present invention realizes that optical fiber stimulated Brillouin scattering gain spectrum widthization is handled by using three pump signals, improve the precision of BOTDA strain measurement systems;Raising signal-to-noise ratio is combined using pulse code and wavelet transformation technique, and then improves measurement accuracy, shorten time of measuring.
Description
Technical field
The invention belongs to Distributed Optical Fiber Sensing Techniques field, and in particular to one kind is based on excited Brillouin gain spectral line width
The high-performance fast B OTDA (Brillouin optical time domain analysis type) that processing, pulse code and wavelet transformation technique are combined that narrows is answered
Become measuring device and data processing method.
Background technology
Fibre optical sensor has small, light-weight, high sensitivity, high pressure resistant, corrosion-resistant, electrical insulating property is good, anti-electromagnetism
The advantages that interference, be widely used in the structure detections such as communications optical cable, river levee, concrete, pipeline, tunnel, bridge.Based on by
The distributed fiberoptic sensor for swashing Brillouin scattering effect is the linear relationship utilized between excited Brillouin frequency shift amount and strain,
Corresponding strain variation at the position is measured by measuring the Brillouin shift of optical fiber everywhere, so as to fulfill the distribution of strain
Measurement.
Excited Brillouin distributed fiberoptic sensor mainly has three types:Brillouin optical time domain analysis (BOTDA) type, cloth
In deep optical frequency domain analysis (BOFDA) type and Brillouin scattering optical time domain reflection (BOTDR) type.BOTDA systems are due to utilizing optical fiber
Stimulated Brillouin scattering effect measurement strains, and brillouin gain line width directly influences the measurement accuracy of strain, brillouin gain
Line width is wider, and the measurement accuracy of strain is lower.The time of measuring of BOTDA systems is also a very crucial factor, is used at present
The reduction time method be using coding techniques improve signal-to-noise ratio so that reduce pendulous frequency, shorten time of measuring, but compile
When the exponent number of code is excessive, the time of measuring of system is also unfavorable for reducing.
The content of the invention
The object of the present invention is to provide one kind to be based on the processing of excited Brillouin gain spectrum width, pulse code and small echo
The high performance BOTDA strain gauge means that converter technique is combined.
The structure of high performance BOTDA strain measurement systems of the present invention is as shown in Figure 1, by laser 1, the first light
Coupler 2, the first modulator 3, the first microwave signal source 20, the first optical filter 4, the second modulator 5, the second microwave signal source
21, DC power supply 22, image intensifer 6, the second photo-coupler 7, optoisolator 8, the second modulator 9, frequency synthesizer 14, second
Optical filter 10, the 4th modulator 11, pulse signal generator 15, scrambler 12, optical circulator 13, sensor fibre 19, photoelectricity
Detector 16, data collecting card 17 and computer 18 form, the first modulator 3, the second modulator 5, the 3rd modulator the 9, the 4th
Modulator 11 is intensity modulator.
1 output frequency of laser is fcOptical signal (Fig. 2 (1)) through 2 one dividing into three of the first photo-coupler, be denoted as respectively
One branch 101,201 and the 3rd branch 301 of the second branch, the light splitting of the first branch 101, the second branch 201, the 3rd branch 301
Than for 3:3:4, wherein, the optical signal of the first branch is sent in the first modulator 3 as light carrier, is then believed by the first microwave
The frequency that number source 20 exports is 2fB(fBFor the stimulated Brillouin scattering frequency shift amount of sensor fibre 19) microwave signal modulation, due to
It is microwave signal modulation by a small margin, only considers carrier wave and single order sideband (fc、fc+2fB、fc-2fB), the letter of the first modulator 3 output
Number it is input in the first optical filter 4, filters carrier wave and single order upper side band, only surplus frequency is fc-2fBSingle order lower sideband signal
(Fig. 2 (2)), the frequency that the first optical filter 4 is exported are fc-2fBSignal is re-fed into the second modulator 5, is believed by the second microwave
The frequency that number source 21 exports is fLSignal modulation, the DC voltage that adjustment DC power supply 22 exports changes second modulator 5
DC offset voltage, makes it be operated in the double-side band output state of carrier wave suppression, its output frequency is fc-2fB+fLAnd fc-2fB-
fLTwo signals (Fig. 2 (3)), the second modulator 5 output signal (fc-2fB+fLAnd fc-2fB-fL) intensity can pass through light
Amplifier 6 amplifies;Frequency is f in the second branch 201cOptical signal and the first branch 101 in the light that amplifies through image intensifer 6 believe
Number it is sent to again through optoisolator 8 in sensor fibre 19 after the second photo-coupler 7, as excited Brillouin in sensor fibre 19
The pump signal of scattering effect, the frequency of three pump signals is respectively fc、fc-2fB+fLAnd fc-2fB-fL(intensity can pass through
Image intensifer 6 controls), pump signal fcThe centre frequency of the gain spectral of generation is fc-fB, pump signal fc-2fB+fLAnd fc-
2fB-fLThe centre frequency for producing loss spectra is fc-fB+fLAnd fc-fB-fL(Fig. 2 (4)), the amplification factor for changing image intensifer 6 can
To change pump signal fc-2fB+fLAnd fc-2fB-fLIntensity, so as to change loss spectra fc-fB+fLAnd fc-fB-fLIntensity, in
Frequency of heart is fc-fBGain spectral and centre frequency be fc-fB+fLAnd fc-fB-fLLoss spectra interacts, and realizes excited Brillouin
The processing that narrows (Fig. 2 (5)) of gain spectral line width, so as to improve measurement accuracy.Optical signal in 3rd branch 301 is input to the 3rd
In modulator 9, the frequency exported by frequency synthesizer 14 is fT(fTFrequency be with fBCentered in the range of 300M press one
Determine frequency interval change) small amplitude signal modulation, the second optical filter 10 filters carrier wave and single order upper side band in modulated signal
Signal (due to being small signal modulation, only carrier wave and single order sideband), a reserve frequency is fc-fTSingle order lower sideband signal, frequency
Rate is fc-fTSignal be input in the 4th modulator 11, the pulse signal modulation (pulse exported by pulse signal generator 15
The high level width of signal is 10ns~100ns, and pulse frequency is 10kHz~100kHz), the signal of the 4th modulator 11 output
It is input in scrambler 12, scrambler 12 is to be randomized the polarization state of light, excited Brillouin is dissipated with eliminating polarization state
The influence of effect is penetrated, pulsed light of the optical signal handled through scrambler 12 as stimulated Brillouin scattering effect in sensor fibre 19
Inputted from I port of optical circulator 13, II port is output and then enter in sensor fibre 19.Narrow from what optoisolator 8 exported
The pulsed light that three pump signals handled and II port from circulator export interacts in sensor fibre 19, works as two-way
Meet stimulated Brillouin scattering phenomenon can occur during stimulated Brillouin scattering condition between the difference on the frequency of light, carry excited Brillouin
Port II of the pump light signals of scattered information through optical circulator inputs, and is visited after being exported from the port III of optical circulator 13 by photoelectricity
Survey device 16 to detect, then handled the data collected feeding computer 18 by data collecting card 17, computer 18 also needs to
The frequency values of the output of frequency synthesizer 14 and pulse signal generator 15 is controlled to export the initial time of pulse, finally by spectrum
Fitting can obtain excited Brillouin gain spectral crest frequency, and then determine the size of strain, time for being sent by pulse and
The time difference that gain spectrum peak occurs may determine that the position that strain occurs.
The strain measurement accuracy formula of BOTDA measuring systems isWherein, Δ νBBrillouin
Gain spectral line width, SNR are system signal noise ratio, it is seen then that reduce brillouin gain line width Δ νBIt can be improved with increase signal-to-noise ratio
The measurement accuracy of system.
The present invention is combined by using S-codes (Simplex codes) codings and wavelet transformation technique and is to improve
The to-noise ratio (SNR) of system and then improve measurement accuracy and shorten the time of measuring of system.It is logical using S-codes coding techniques
Cross and Hadamard (hadmard) conversion carried out to s-matrix (S-matrix), then reduce noise by being superimposed, improve signal-to-noise ratio,
S-matrix is a unipolarity matrix, is made of 1 and 0, is encoded for S-codes, when the exponent number increase of s-matrix, decoding
Time is elongated, makes the time of measuring of system increase, and using the exponent number of appropriate s-matrix, at data processing end, decoding collects
Recycle wavelet transformation technique to carry out denoising after data, improve the to-noise ratio of system and then improve measurement accuracy and contract
The time of measuring of short system.The key of the coding techniques is structure S-matrix matrixes, to measure different positions along optical fiber at the same time
The Brillouin intensity put.Illustrate to improve SNR using the coding techniques by taking a three rank S-matrix matrixes as an example, it is assumed that
For a pulse P1(t) it is launched into measured optical fiber, obtained response is ψ1(t).Define new pulse P at the same time2
(t)=P1(t- τ), P3(t)=P1(t-2 τ), their response are respectively:ψ2(t)=ψ1(t- τ), ψ3(t)=ψ1(t-2τ);It is logical
Crossing transmitting S-codes sequences can obtain:
Then S-matrix can be used to be write as the form of hadamard conversion:
E (t) is the noise of measurement in formula.From above formula, code length is the exponent number of S-matrix.In order to recover ψ1
(t) need to carry out Hadmard inverse transformations, obtain following expression:
Then
Then the result after can be processed is:
Last Mean Square Error is:
And for common pulse situation, the Mean Square Error of iteration is σ three times2/ 3, obtained when being 3 for code length
Pattern gain beTherefore for the S-codes that length is L, the pattern gain that can derive it is
It is thus possible to make system signal-to-noise ratio be promoted to accordingly it is originalTimes.
Encoded for S-codes, when the exponent number increase of s-matrix, decoding time is elongated, increases the time of measuring of system
Add, the present invention recycles wavelet transformation using the exponent number of appropriate s-matrix after decoding the data collected at data processing end
Technology carries out denoising, can greatly save the time of one-shot measurement.
Following three steps of wavelet threshold shrinkage method denoising point:
1) orthogonal wavelet transformation of signals and associated noises is calculated.Signals and associated noises X for length for N, N=2 might as well be setJ, utilize
The fast algorithm of orthogonal wavelet transformation obtains the scale coefficient { v under low resolution L (0≤L < J)L,k, k=1 ..., 2J, and
Wavelet coefficient { w under each resolution ratioj,k, j=L, L+1 ..., J-1, k=1 ..., 2j, wherein scale coefficient and wavelet coefficient
Common N.When handling border, frequently with periodic extension method.
2) non-linear threshold processing is carried out to wavelet coefficient.To keep the global shape of signal constant, retain all low
Frequency coefficient vL,k, k=1 ..., 2L.Take threshold valueTo each wavelet coefficient, at hard thresholding method
Reason:
Hard -threshold:
That is, the absolute value of the wavelet coefficients of signals and associated noises compared with selected threshold value λ, less than or equal to threshold value
Point vanishing, the point more than threshold value remain unchanged.
3) inverse wavelet transform is carried out.By all low frequency scale coefficients, and do via the wavelet coefficient after threshold process inverse
Wavelet transformation is reconstructed, the estimate for the original signal being restored.
Finally, be derived from be measured field occur mutation peak come judge damage appearance position.
When measurement starts, the output frequency of laser 1 is arranged to fc=193.41THz (corresponding wavelength 1550nm), meter
Calculation machine sends instruction initiation culture synthesizer 14 and starts frequency sweep, while 15 output pulse width of starting impulse signal generator is 10
~100ns, frequency press the pulse signal of s-matrix coding for 10kHz~100kHz, and the exponent number of coding can take 3~255 ranks, with 3
Exemplified by rank coding, input first shown in its form of the pulsed optical signals of first group of sequence such as Fig. 3 (a), 18 recording impulse of computer
Signal generator 15 sends the time of pulse signal, and this time is the measurement time started, and the excited Brillouin of sensor fibre 19 is frequently
Shifting amount fBFor 9~11GHz, the centre frequency of frequency synthesizer 14 is set to fB, swept frequency range is 100~300MHz, and frequency step is
0.1~0.3MHz, the original frequency of setpoint frequency synthesizer 14, the microwave signal pass through second again after the 3rd modulator 9
Optical filter 10 filters, and retains lower sideband signal, is then modulated by the pulse signal shown in Fig. 3 (a) by the 4th modulator 11
Sequentially entered afterwards as pulsed light in sensor fibre 19,19 other end input energy of sensor fibre makes excited Brillouin line width narrow
Three pump signals, the pulsed optical signals of coding propagate along sensor fibre 19 successively, the pulsed optical signals of first group of sequence
When being propagated in sensor fibre 19, the light for carrying stimulated Brillouin scattering information can be received in the receiving terminal of photodetector 16
Signal, is sent in computer 18 through data collecting card 17, the frequency values of keep frequency synthesizer 14, sequentially inputs 3 rank s-matrix
The pulsed optical signals of coded residual are shown in Fig. 3 (b) and 3 (c), these signals detected are decoded in computer 18 and are obtained along light
The light intensity of each point in fine length, then carries out wavelet transform process to the response that decoding obtains in computer 18, just completes one
Secondary sweep measurement.Then the frequency step value set by frequency synthesizer 14 is successively according to above-mentioned coding and wavelet transform process
Process completes the sweep measurement of each Frequency point, finally obtains sensor fibre 19 and goes up the light intensity of each point along its length with frequency values
Change curve, frequency when brillouin gain reaches maximum can be obtained by Lorentz fit, maximum frequency values are exactly
Optical fiber stimulated Brillouin frequency shift amount, since strain and frequency shift amount are in a linear relationship, and then determines the size of strain, meanwhile, calculate
Machine 18 sends the time of pulse by pulse signal generator 15 and the time difference of gain spectrum peak appearance may determine that strain hair
Raw position, can realize the positioning to strain.
The present invention realizes that optical fiber stimulated Brillouin scattering gain spectrum widthization is handled by using three pump signals, carries
The precision of high BOTDA strain measurement systems;Raising signal-to-noise ratio is combined using pulse code and wavelet transformation technique, and then is improved
Measurement accuracy, shorten time of measuring.
Brief description of the drawings
Fig. 1:High-performance BOTDA strain measurement system schematic diagrames;
Fig. 2:Stimulated Brillouin scattering composes the processing procedure schematic diagram that narrows;
Fig. 3:The pulsed optical signals of three ranks coding;
Fig. 4:Individually using different coding order corresponding response curve during S codings;
Fig. 5:S codings different coding order corresponding response curve when being combined with wavelet transformation.
Embodiment
Embodiment 1:
The TSL-510 tunable lasers of Santec companies are selected to make carrier wave light source, the wave-length coverage of laser 1 is 1510nm
~1630nm wavelength, sets wavelength as 1550nm (respective frequencies fc=193.41THz);First modulator 3, the second modulator
5th, the 3rd modulator 9 and the 4th modulator 11 are the MXAN-LN-40 of Photline companies, bandwidth 32GHz;First light is filtered
4 and second optical filter 10 of ripple device be Santec companies adjustable light wave-filter, model OTF-950, wavelength tuning range
For 1548nm to 1552nm, line width is less than 10GHz;Image intensifer 6 be Co., Ltd of ZTE Corporation erbium-doped fiber amplifier, ripple
Long scope is 1530~1565nm, and amplification factor is more than 35 times;First microwave signal source 20 and the second microwave signal source 21 are peace
The 8257D of Jie Lun companies;Frequency synthesizer 14 is the RJUFS020180-1K of Chengdu Ren Jian microwave technologies Co., Ltd, and output is frequently
Rate scope is 2-18GHz;Pulse signal generator 15 is agilent company 81131A, frequency range 1Hz-400MHz;Photoelectricity
Detector 16 is the SD-48, bandwidth 35GHz of Imtech;Sensor fibre 19 is the optical fiber of Yangtze Optical Fiber and Cable Company Ltd,
Excited Brillouin gain line width is ΓB=40MHz, Brillouin shift amount fB=10GHz, length are 500 meters, gain and loss peak
Value coefficient is 5;The isolation of optoisolator 8 is more than 40dB;Scrambler 12 is the PCD-104 of high Micron Technology Co., Ltd;Data
Capture card 17 be NI companies PCI-5112 data collecting cards, sample frequency 100MSPS.
Corresponding instrument and equipment is connected by Fig. 1, opens instrument power source, the wavelength of laser 1 is set as 1550nm, first
The output frequency of microwave signal source 20 is 20GHz, and the set of frequency of the second microwave signal source 21 is 20MHz, and image intensifer 6 is put
Big multiple is arranged to 25 times, by the effect of 101 branches and 201 branches of Fig. 1, the line width that excited Brillouin is composed can be made to narrow
To 4.14MHz.
When measurement starts, the output frequency of laser 1 is arranged to fc=193.41THz (corresponding wavelength 1550nm), meter
Calculation machine sends instruction initiation culture synthesizer 14 and starts frequency sweep, while 15 output pulse width of starting impulse signal generator is
10ns, frequency press the pulse signal of s-matrix coding for 10kHz, the exponent number of coding is respectively set to 3,7,15,61,63,127,
255 ranks, 18 recording pulse signal generator of computer send the time of pulse signal, and this time is to measure time started, optical fiber
Excited Brillouin frequency shift amount be 10GHz, the original frequency of setpoint frequency synthesizer 14 is 9.9GHz, and swept frequency range is
The microwave signal of 200MHz, frequency step 0.1MHz, this 9.9GHz are filtered by the second light again after the 3rd modulator 9
Device 10 filters, and retains lower sideband signal fc- 9.9GHz, the arteries and veins encoded by s-matrix then exported by pulse signal generator 15
Rush signal to sequentially enter in sensor fibre 19 by being used as pulsed light after the modulation of the 4th modulator 11,19 other end of sensor fibre
Input energy makes three pump signals that excited Brillouin line width narrows, by the pulsed optical signals of s-matrix coding successively along sensing
Optical fiber 19 is propagated, when the pulsed optical signals of first group of sequence are propagated in sensor fibre 19, in the receiving terminal of photodetector 16
The optical signal for carrying stimulated Brillouin scattering information can be received, is sent to through data collecting card 17 in computer 18, keeps frequency
The frequency values 9.9GHz of rate synthesizer 14, the exponent number encoded by s-matrix sequentially inputs remaining pulsed optical signals, by what is detected
These signals decode in computer 18 and obtain frequency when being 9.9GHz, the light intensity of each point, Ran Hou along in fiber length
Wavelet transform process is carried out to the response that decoding obtains in computer 18, just completes a sweep measurement.Then frequency synthesis is pressed
The frequency step value that device 14 is set completes sweeping in the range of 200MHz according to the process of above-mentioned coding and wavelet transform process successively
Frequency measures, and finally obtains the light intensity of each point in 19 length of sensor fibre with the change curve of frequency values, can by Lorentz fit
To obtain frequency when brillouin gain reaches maximum, maximum frequency values are exactly optical fiber stimulated Brillouin frequency shift amount, due to answering
Become and frequency shift amount is in a linear relationship, and then determine the size of strain, meanwhile, computer 18 is sent by pulse signal generator 15
The time difference that the time of pulse and gain spectrum peak occur may determine that the position that strain occurs, and can realize and strain is determined
Position.
At the 470m of sensor fibre 19 plus strain, when individually being encoded using 3,7,15,31,63,127,255 rank S, decoding
After meet with a response shown in curve such as Fig. 4 (b) -4 (h), correspond to response during 3,7,15,31,63,127,255 ranks coding respectively,
Fig. 4 (a) is decoded response curve when not using coding, as can be seen that the order of coding is 255 from Fig. 4 (a) -4 (h)
The strain at 470 meters can be told during rank, after decoding, processing time at this time is 4439.7244s.
Fig. 5 is that 3,7,15,31 rank S are encoded and wavelet transformation junction understands response after code, Fig. 5 (a)-Fig. 5 (d) points
Not Wei 3,7,15,31 rank S codings and response during wavelet transformation, from figure 5 it can be seen that when coding exponent number increases to 31 ranks solution
The strain at 470m can be told with reference to wavelet transformation after code, at this time, 31 ranks coding amounts to processing time with reference to wavelet analysis
For 45.1673s.Signal processing time can be obviously shortened by being combined using S codings and wavelet transformation, increase the quick of system
Property.
Claims (3)
- A kind of 1. quick Brillouin optical time domain analysis type strain gauge means, it is characterised in that:By laser (1), the first optocoupler It is clutch (2), the first modulator (3), the first microwave signal source (20), the first optical filter (4), the second modulator (5), second micro- Ripple signal source (21), DC power supply (22), image intensifer (6), the second photo-coupler (7), optoisolator (8), the second modulator (9), frequency synthesizer (14), the second optical filter (10), the 4th modulator (11), pulse signal generator (15), scrambler (12), optical circulator (13), sensor fibre (19), photodetector (16), data collecting card (17) and computer (18) composition;Laser (1) output frequency is fcOptical signal through the first photo-coupler (2) one dividing into three, be denoted as the first branch respectively (101), the second branch (201) and the 3rd branch (301), wherein, the optical signal of the first branch (101) is sent to as light carrier In first modulator (3), the frequency then exported by the first microwave signal source (20) is 2fBMicrowave signal modulation, due to being Microwave signal is modulated by a small margin, only considers carrier wave and single order sideband, and the signal of the first modulator (3) output is input to the filter of the first light In ripple device (4), carrier wave and single order upper side band are filtered, only surplus frequency is fc-2fBSingle order lower sideband signal, by the first optical filter (4) frequency of output is fc-2fBSignal is re-fed into the second modulator (5), the frequency exported by the second microwave signal source (21) For fLSignal modulation, adjustment DC power supply (22) output DC voltage change the second modulator (5) DC offset voltage, It is set to be operated in the double-side band output state of carrier wave suppression, its output frequency is fc-2fB+fLAnd fc-2fB-fLTwo signals, The intensity of second modulator (5) output signal can be amplified by image intensifer (6);Frequency is f in the second branch (201)cOptical signal and the first branch (101) in the optical signal that amplifies through image intensifer (6) It is sent to again through optoisolator (8) in sensor fibre (19) after the second photo-coupler (7), as being excited in sensor fibre (19) The pump signal of Brillouin scattering effect, the frequency of three pump signals is respectively fc、fc-2fB+fLAnd fc-2fB-fL, pumping letter Number fcThe centre frequency of the gain spectral of generation is fc-fB, pump signal fc-2fB+fLAnd fc-2fB-fLProduce the center frequency of loss spectra Rate is fc-fB+fLAnd fc-fB-fL, pump signal f can be changed by changing the amplification factor of image intensifer (6)c-2fB+fLAnd fc- 2fB-fLIntensity, so as to change loss spectra fc-fB+fLAnd fc-fB-fLIntensity, centre frequency fc-fBGain spectral and center Frequency is fc-fB+fLAnd fc-fB-fLLoss spectra interacts, and realizes the processing that narrows of excited Brillouin gain spectral line width, so that Improve measurement accuracy;Optical signal in 3rd branch (301) is input in the 3rd modulator (9), is by the frequency that frequency synthesizer (14) exports fTSmall amplitude signal modulation, the second optical filter (10) filters carrier wave and single order upper side band signal in modulated signal, only retains Frequency is fc-fTSingle order lower sideband signal, frequency fc-fTSignal be input in the 4th modulator (11), by pulse signal The pulse signal modulation of generator (15) output, the signal of the 4th modulator (11) output are input in scrambler (12), disturb partially Device (12) is in order to which the polarization state of light is randomized, to eliminate influence of the polarization state to stimulated Brillouin scattering effect, through disturbing partially The optical signal of device (12) processing as stimulated Brillouin scattering effect in sensor fibre (19) pulsed light from optical circulator (13) I port input, II port is output and then enter in sensor fibre (19);From the three of the processing that narrows of optoisolator (8) output A pump signal and pulsed light interaction in sensor fibre (19) from the output of II port of circulator, when the frequency of two-way light Meet stimulated Brillouin scattering phenomenon can occur during stimulated Brillouin scattering condition between rate difference, carry stimulated Brillouin scattering letter Port II of the pump light signals of breath through optical circulator inputs, by photodetector after being exported from the port III of optical circulator (13) (16) detect, the data collected then are sent into computer (18) by data collecting card (17) is handled, and computer (18) is also Need to control the frequency values of frequency synthesizer (14) output and the initial time of pulse signal generator (15) output pulse, finally Excited Brillouin gain spectral crest frequency can be obtained by spectrum simulation, and then determine the size of strain, sent by pulse Time and gain spectrum peak occur time difference may determine that strain occur position.
- A kind of 2. quick Brillouin optical time domain analysis type strain gauge means as claimed in claim 1, it is characterised in that:First Modulator (3), the second modulator (5), the 3rd modulator (9) and the 4th modulator (11) are intensity modulator.
- A kind of 3. quick Brillouin optical time domain analysis type strain gauge means as claimed in claim 1, it is characterised in that:First The splitting ratio of branch (101), the second branch (201) and the 3rd branch (301) is 3:3:4.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510686227.0A CN105241390B (en) | 2015-10-21 | 2015-10-21 | Quick Brillouin optical time domain analysis type strain gauge means and data processing method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510686227.0A CN105241390B (en) | 2015-10-21 | 2015-10-21 | Quick Brillouin optical time domain analysis type strain gauge means and data processing method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105241390A CN105241390A (en) | 2016-01-13 |
CN105241390B true CN105241390B (en) | 2018-05-15 |
Family
ID=55039122
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510686227.0A Expired - Fee Related CN105241390B (en) | 2015-10-21 | 2015-10-21 | Quick Brillouin optical time domain analysis type strain gauge means and data processing method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105241390B (en) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105674905B (en) * | 2016-01-15 | 2018-04-03 | 华北电力大学(保定) | The single-ended vector B OTDA dynamic strain measurement methods of the pre- pumping of pulse and device |
CN105571507B (en) * | 2016-01-15 | 2018-04-03 | 华北电力大学(保定) | A kind of method and its measurement apparatus of single-ended vector B OTDA dynamic strain measurements |
CN105783758B (en) * | 2016-04-07 | 2018-06-08 | 华北电力大学(保定) | A kind of single-ended vector B OTDA dynamic strain measurement methods of self-heterodyne and device |
CN105783762B (en) * | 2016-05-10 | 2018-04-06 | 太原理工大学 | The brillouin distributed optical fiber sensing device and method of chaos correlation method positioning |
CN105784195B (en) * | 2016-05-10 | 2018-04-06 | 太原理工大学 | The distribution type optical fiber sensing equipment and method of single-ended chaos Brillouin optical time domain analysis |
CN106525098B (en) * | 2016-12-01 | 2018-11-13 | 南京大学 | Eliminating light amplification leads to the system and method for coded pulse luminous power unevenness phenomenon |
CN107065390B (en) * | 2017-06-28 | 2019-06-21 | 吉林大学 | Microwave signal generating method and device based on stimulated Brillouin scattering effect and frequency comb |
CN112016201B (en) * | 2020-08-27 | 2022-12-02 | 安徽理工大学 | DFOS strain-based deep stope advanced support pressure evolution model reconstruction method |
CN113310563A (en) * | 2021-04-22 | 2021-08-27 | 成都海恒之星科技合伙企业(有限合伙) | Distributed optical fiber vibration sensing device and method for improving positioning accuracy |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1752835A (en) * | 2004-09-23 | 2006-03-29 | 朗迅科技公司 | Multiple pump parametric apparatus having no idler broadening |
CN202195827U (en) * | 2011-08-09 | 2012-04-18 | 中国计量学院 | Super remote pulse code distributed fiber Brillouin sensor with integration of fiber Brillouin frequency shifter |
CN202255424U (en) * | 2011-08-09 | 2012-05-30 | 中国计量学院 | Pulse coding optical fiber Brillouin optical time domain analyzer |
CN102680136A (en) * | 2012-05-31 | 2012-09-19 | 北京邮电大学 | Distributed stimulated Brillouin temperature strain sensing system based on double-sideband modulation |
CN203519150U (en) * | 2013-09-06 | 2014-04-02 | 山东大学 | Shake monitoring system based on Brillouin distributed optical fiber sensing |
CN103715480A (en) * | 2014-01-20 | 2014-04-09 | 吉林大学 | Single-band-pass tunable microwave photonic filter with ultrahigh quality factor |
CN104776871A (en) * | 2015-01-30 | 2015-07-15 | 佛山科学技术学院 | Optical fiber Brillouin distributed type measuring light path, device and method |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2010161646A (en) * | 2009-01-08 | 2010-07-22 | Sumitomo Electric Ind Ltd | Optical transmitter and optical communication system |
-
2015
- 2015-10-21 CN CN201510686227.0A patent/CN105241390B/en not_active Expired - Fee Related
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1752835A (en) * | 2004-09-23 | 2006-03-29 | 朗迅科技公司 | Multiple pump parametric apparatus having no idler broadening |
CN202195827U (en) * | 2011-08-09 | 2012-04-18 | 中国计量学院 | Super remote pulse code distributed fiber Brillouin sensor with integration of fiber Brillouin frequency shifter |
CN202255424U (en) * | 2011-08-09 | 2012-05-30 | 中国计量学院 | Pulse coding optical fiber Brillouin optical time domain analyzer |
CN102680136A (en) * | 2012-05-31 | 2012-09-19 | 北京邮电大学 | Distributed stimulated Brillouin temperature strain sensing system based on double-sideband modulation |
CN203519150U (en) * | 2013-09-06 | 2014-04-02 | 山东大学 | Shake monitoring system based on Brillouin distributed optical fiber sensing |
CN103715480A (en) * | 2014-01-20 | 2014-04-09 | 吉林大学 | Single-band-pass tunable microwave photonic filter with ultrahigh quality factor |
CN104776871A (en) * | 2015-01-30 | 2015-07-15 | 佛山科学技术学院 | Optical fiber Brillouin distributed type measuring light path, device and method |
Non-Patent Citations (2)
Title |
---|
The extraction of the Brillouin frequency shift;L.Q.Li et.al;《2014 International Conference on Artificial Intelligence and Industrial Application 》;20140921;第67卷;第881-887页 * |
宽调谐范围的单通带微波光子滤波器的研究;郭静 等;《光电子·激光》;20140730;第25卷(第7期);第1274-1278页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105241390A (en) | 2016-01-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105241390B (en) | Quick Brillouin optical time domain analysis type strain gauge means and data processing method | |
CN105675031B (en) | Brillouin optical time domain analysis instrument based on pre- pumping pulse and gray encoding | |
CN103604446A (en) | Multi-channel fiber bragg grating absolute wavelength demodulation system based on single detector and method thereof | |
CN109297425B (en) | Brillouin optical time domain reflectometer modulated by physical random number | |
CN104180833A (en) | Optical time domain reflectometer simultaneously sensing temperature and stress | |
CN110108346B (en) | Optical fiber vibration sensor based on delay phase modulation chirp pulse pair | |
CN104019836B (en) | Based on relevant dual-pulse time series technique Brillouin optical time domain analysis instrument and the method utilizing this analyser suppression common-mode noise | |
CN102589617B (en) | Full-fiber type multi-parameter monitoring system based on chirped fiber grating | |
CN105784195A (en) | Single-end chaotic Brillouin optical time-domain analysis distributed fiber sensing device and method | |
EP3207340A1 (en) | Distributed brillouin sensing using correlation | |
CN104048684A (en) | OTDR device and method based on coded pulse optical signals | |
CN108917804A (en) | Quick long-distance distributed Brillouin light fiber sensing equipment based on chirp chain | |
CN102829811A (en) | Method for increasing detection speed of BOTDA (Brillouin Optical Time Domain Analyzer) on basis of optical combing and coding technologies | |
CN103604450A (en) | Seed injection BOTDR distributed optical fiber sensing system | |
CN108955733A (en) | A method of it is encoded based on assembled pulse and improves BOTDA system signal noise ratio | |
CN106610303B (en) | A kind of BOTDR method for sensing and system based on FFT and coding | |
CN103674082B (en) | A kind of High-spatial-resolutoptical optical frequency domain reflectometer system based on four-wave mixing process | |
Wang et al. | Robust and fast temperature extraction for Brillouin optical time-domain analyzer by using denoising autoencoder-based deep neural networks | |
CN105403257A (en) | Distribution type Brillouin time domain analysis system | |
CN113049014A (en) | Time-frequency multiplexing BOTDA system based on pumping light frequency sweep and sensing method | |
CN102564481A (en) | Method and device for improving signal-to-noise ratio of distributed optical fiber Brillouin sensor | |
CN102706475B (en) | Distribution-type optical-fiber Raman temperature sensor adopting circulating pulse coding and decoding and Rayleigh demodulation | |
CN111623902B (en) | Distributed optical fiber Raman temperature sensor based on intensity modulation chirp pulse compression | |
CN112325911B (en) | Silicon-based micro-ring pulse coding time division multiplexing dynamic wavelength demodulation method | |
CN104296673A (en) | Brillouin spectrum signal quality improving method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180515 Termination date: 20201021 |
|
CF01 | Termination of patent right due to non-payment of annual fee |