CN109900382A - Fiber temperature signal demodulation method, device and fiber optic temperature (FBG) demodulator - Google Patents
Fiber temperature signal demodulation method, device and fiber optic temperature (FBG) demodulator Download PDFInfo
- Publication number
- CN109900382A CN109900382A CN201910307202.3A CN201910307202A CN109900382A CN 109900382 A CN109900382 A CN 109900382A CN 201910307202 A CN201910307202 A CN 201910307202A CN 109900382 A CN109900382 A CN 109900382A
- Authority
- CN
- China
- Prior art keywords
- signal
- stokes
- current signal
- temperature
- luminous flux
- 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.)
- Granted
Links
- 239000000835 fiber Substances 0.000 title claims abstract description 71
- 238000000034 method Methods 0.000 title claims abstract description 39
- 238000001069 Raman spectroscopy Methods 0.000 claims abstract description 53
- 230000004907 flux Effects 0.000 claims abstract description 53
- 239000013307 optical fiber Substances 0.000 claims description 25
- 238000004590 computer program Methods 0.000 claims description 8
- 238000006243 chemical reaction Methods 0.000 claims description 6
- 230000009466 transformation Effects 0.000 claims description 5
- 230000005611 electricity Effects 0.000 claims description 3
- 238000001228 spectrum Methods 0.000 claims description 3
- 230000003595 spectral effect Effects 0.000 description 19
- 238000010586 diagram Methods 0.000 description 11
- 230000006870 function Effects 0.000 description 9
- 230000005540 biological transmission Effects 0.000 description 4
- 239000000243 solution Substances 0.000 description 4
- 230000000694 effects Effects 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000003062 neural network model Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000000644 propagated effect Effects 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000000151 deposition Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 239000004615 ingredient Substances 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 239000000382 optic material Substances 0.000 description 1
Abstract
The invention discloses a kind of fiber temperature signal demodulation method, device and fiber optic temperature (FBG) demodulators, this method comprises: obtaining Raman diffused light subsignal；The Stokes signal is converted into the first current signal, and the Anti-Stokes signal is converted into the second current signal；First current signal and second current signal are debugged, so that the waveform of first current signal reaches default coincidence threshold value be overlapped degree of the waveform of second current signal in predetermined time section；The luminous flux of the luminous flux and the second current signal that reach default the first current signal for being overlapped threshold value according to coincidence degree calculates temperature value.Technical solution of the present invention eliminates influence of the difference of Stokes signal and Anti-Stokes signal to temperature value, improves temperature value computational accuracy by debugging Stokes signal and the current waveform of Anti-Stokes signal.
Description
Technical field
The present invention relates to fiber optic temperature (FBG) demodulator technical fields, in particular to a kind of fiber temperature signal demodulation side
Method, device and fiber optic temperature (FBG) demodulator.
Background technique
It is existing to be measured in optical fiber based on the thermometric instruments of Raman scattering principle using the Raman scattering effect of light
Temperature value.Thermometric instruments transmitting laser signal is propagated in a fiber, and meeting back reflection returns thermometric instruments, by anti-
The Stokes signal in signal and Anti-Stokes signal that are emitted back towards can determine the temperature value of optical fiber.
However, only only accounting in the existing thermometric instruments based on Raman scattering effect due to fiber optic materials pair
The difference that the decaying of light, the decaying of micro-bend introducing, the variation of injection light pulse power and thermometric instruments respond optical signal
The problems such as caused by temperature value measurement error, do not consider the difference of Stokes signal and Anti-Stokes signal itself to temperature
It is influenced caused by value.
Summary of the invention
In view of the above problems, the embodiment of the present invention is designed to provide a kind of fiber temperature signal demodulation method, device
And fiber optic temperature (FBG) demodulator, so as to solve the deficiencies in the prior art.
According to embodiment of the present invention, a kind of fiber temperature signal demodulation method is provided, this method comprises:
Obtain Raman diffused light subsignal, wherein the Raman diffused light subsignal include Stokes signal and it is anti-this
Lentor signal；
The Stokes signal is converted into the first current signal, and the Anti-Stokes signal is converted to second
Current signal；
First current signal and second current signal are debugged, so that the wave of first current signal
Shape reaches default coincidence threshold value be overlapped degree of the waveform of second current signal in predetermined time section；
Reach the luminous flux of default the first current signal for being overlapped threshold value and the light of the second current signal according to coincidence degree
Flux calculates temperature value.
In above-mentioned fiber temperature signal demodulation method, " the obtaining Raman diffused light subsignal " includes:
Obtain the reflection signal of laser in a fiber；
Reflection signal progress Fourier transformation is obtained into the reflection signal of frequency domain, according to Raman diffused light subsignal
Spectral range obtains Raman diffused light subsignal in the reflection signal of the frequency domain.
In above-mentioned fiber temperature signal demodulation method, after described " obtaining Raman diffused light subsignal " further include:
Denoising is carried out according to the Raman diffused light subsignal of the noise model pre-established to acquisition, wherein
The noise model is learnt to obtain by the noise in the Raman diffused light subsignal to preset predetermined number.
It is described " to reach default according to coincidence degree and be overlapped the of threshold value in above-mentioned fiber temperature signal demodulation method
The luminous flux of the luminous flux of one current signal and the second current signal calculates temperature value " include:
Reach default be overlapped in threshold portion waveform in the coincidence degree and presets reference position；
According to the temperature value of the reference position, the corresponding Stokes signal luminous flux in the reference position and this anti-support
Gram corresponding Stokes signal luminous flux of this signal luminous flux and position to be measured and Anti-Stokes signal luminous flux calculate institute
State the temperature value of position to be measured.
In above-mentioned fiber temperature signal demodulation method, it is calculated by the following formula the temperature value of the position to be measured:
Wherein, T is the temperature value of position to be measured, T0For the temperature value of reference position, h is Planck's constant, k be Bohr hereby
Graceful constant, Δ v are the Raman Phonon frequency of optical fiber, are a constant, ΦSIt (T) is the corresponding Stokes letter in the position to be measured
Number luminous flux, ΦASIt (T) is the corresponding Anti-Stokes signal luminous flux in the position to be measured, ΦS(T0) it is the reference position
Corresponding Stokes signal luminous flux, ΦAS(T0) it is the corresponding Anti-Stokes signal luminous flux in the reference position.
In above-mentioned fiber temperature signal demodulation method, further includes:
The position to be measured is determined according to the time that reaches of the launch time of the laser and the reflection signal.
It is described " according to the launch time of the laser and the reflection in above-mentioned fiber temperature signal demodulation method
The time that reaches of signal determines the position to be measured " include:
Determine the laser signal described according to the time that reaches of the launch time of the laser and the reflection signal
Propagation distance in optical fiber；
The position to be measured is determined according to the transmitting position of the laser and the propagation distance.
According to another implementation of the invention, a kind of fiber temperature signal demodulating equipment is provided, which includes:
Module is obtained, for obtaining Raman diffused light subsignal, wherein the Raman diffused light subsignal includes stoke
This signal and Anti-Stokes signal；
Conversion module, for the Stokes signal to be converted to the first current signal, and by the anti-Stokes
Signal is converted to the second current signal；
Debugging module, for being debugged to first current signal and second current signal, so that described
The waveform of one current signal reaches default weight be overlapped degree of the waveform of second current signal in predetermined time section
Close threshold value；
Computing module, for reaching the luminous flux and second of default the first current signal for being overlapped threshold value according to coincidence degree
The luminous flux of current signal calculates temperature value.
Another embodiment according to the present invention provides a kind of fiber optic temperature (FBG) demodulator, the fiber optic temperature (FBG) demodulator
Including memory and processor, the memory runs the computer journey for storing computer program, the processor
Sequence is so that the fiber optic temperature (FBG) demodulator executes above-mentioned fiber temperature signal demodulation method.
Yet another embodiment according to the present invention provides a kind of computer readable storage medium, is stored with above-mentioned
The computer program used in fiber optic temperature (FBG) demodulator.
The technical scheme provided by this disclosed embodiment may include it is following the utility model has the advantages that
A kind of fiber temperature signal demodulation method, device and fiber optic temperature (FBG) demodulator in the present invention, to Stokes signal
And the waveform of the corresponding current signal of Anti-Stokes signal is debugged, to eliminate Stokes signal and anti-Stokes letter
Number influence of the difference of itself to temperature；According to the lumen meter of Stokes signal and Anti-Stokes signal after debugging
Temperature value is calculated, temperature value computational accuracy is improved.
To enable the above objects, features and advantages of the present invention to be clearer and more comprehensible, preferred embodiment is cited below particularly, and cooperate
Appended attached drawing, is described in detail below.
Detailed description of the invention
In order to illustrate more clearly of technical solution of the present invention, letter will be made to attached drawing needed in the embodiment below
It singly introduces, it should be understood that the following drawings illustrates only certain embodiments of the present invention, therefore is not construed as to the present invention
The restriction of protection scope for those of ordinary skill in the art without creative efforts, can be with root
Other relevant attached drawings are obtained according to these attached drawings.
Fig. 1 shows a kind of flow diagram of fiber temperature signal demodulation method of first embodiment of the invention offer.
Fig. 2 shows a kind of waveform diagrams for reflection signal that first embodiment of the invention provides.
Fig. 3 shows the first current signal and the second current signal after a kind of debugging of first embodiment of the invention offer
Waveform diagram.
Fig. 4 shows a kind of flow diagram of fiber temperature signal demodulation method of second embodiment of the invention offer.
Fig. 5 shows a kind of flow diagram of fiber temperature signal demodulation method of third embodiment of the invention offer.
Fig. 6 shows a kind of structural schematic diagram of fiber temperature signal demodulating equipment of fourth embodiment of the invention offer.
Main element symbol description:
400- fiber temperature signal demodulating equipment；410- obtains module；420- conversion module；430- debugging module；440-
Computing module.
Specific embodiment
Below in conjunction with attached drawing in the embodiment of the present invention, technical solution in the embodiment of the present invention carries out clear, complete
Ground description, it is clear that described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.Usually exist
The component of the embodiment of the present invention described and illustrated in attached drawing can be arranged and be designed with a variety of different configurations herein.Cause
This, is not intended to limit claimed invention to the detailed description of the embodiment of the present invention provided in the accompanying drawings below
Range, but it is merely representative of selected embodiment of the invention.Based on the embodiment of the present invention, those skilled in the art are not doing
Every other embodiment obtained under the premise of creative work out, shall fall within the protection scope of the present invention.
Embodiment 1
Fig. 1 shows a kind of flow diagram of fiber temperature signal demodulation method of first embodiment of the invention offer.
The fiber temperature signal demodulation method includes the following steps:
In step s 110, Raman diffused light subsignal is obtained.
It is generated due to Raman scattering
Further, obtaining Raman diffused light subsignal includes:
Obtain the reflection signal of laser in a fiber；Reflection signal progress Fourier transformation is obtained into the reflection of frequency domain
Signal obtains Raman scattering photon letter according to the spectral range of Raman diffused light subsignal in the reflection signal of the frequency domain
Number.
Specifically, when measuring the temperature value of optical fiber, transmitting laser signal is to which in the optical fiber of testing temperature, laser is believed first
It number is propagated in the optical fiber, since the amorphous material of optical fiber is in the uneven texture of micro-space, has sub-fraction light that can send out
Raw scattering.Scattering envelope in optical fiber includes Rayleigh scattering, Brillouin scattering and Raman scattering, wherein in Raman scattering photon
Carried in Stokes photon and non-Stokes photon be optical fiber temperature value, be influence fiber optic temperature resolution ratio it is main
Factor.
Due to including Rayleigh scattering photon signal, Brillouin scattering in the reflection signal that is reflected back in the slave optical fiber that receives
Photon signal and Raman diffused light subsignal.Due to the spectral range and Brillouin-scattered photons signal of Rayleigh scattering photon signal
Spectral range, the spectral range of Raman diffused light subsignal it is different, then, the reflection signal that can will receive time domain carries out Fu
In leaf transformation, obtain the reflection signal of frequency domain.
For example, Fourier transformation can be carried out by reflection signal of the following formula to time domain:
Wherein, F (ω) is the reflection signal of frequency domain, and f (t) is the reflection signal of time domain.
It is illustrated in figure 2 the spectrum curve of the reflection signal of frequency domain, according to the spectral range of Rayleigh scattering photon signal, cloth
In the spectral range of deep scattered photon signal and the spectral range of Raman diffused light subsignal it is found that with incident laser signal frequency
The identical spectral line A1 of rate is the spectral line of Rayleigh scattering photon signal, neighbouring two spectrums with the spectral line A1 of Rayleigh scattering photon signal
Line A2 is the spectral line of Brillouin-scattered photons signal, spectral line A2 and the Rayleigh scattering photon signal of Brillouin-scattered photons signal
The frequency difference of spectral line A1 is generally 10-1~100cm magnitude.A little spectral line is to draw farther out by spectral line A1 apart from Rayleigh scattering photon signal
The spectral line A31 and A32 of graceful scattered photon signal.
In frequency domain, Raman diffused light subsignal is divided into Stokes signal and Anti-Stokes signal, Stokes letter
Number frequency are as follows:
vs=v0-Δv
Wherein, vsFor the frequency of Stokes signal, v0For the frequency of laser signal incident in optical fiber, Δ v is optical fiber sound
The vibration frequency of son, wherein Δ v=1.32 × 1013Hz。
The frequency of Anti-Stokes signal are as follows:
va=v0+Δv
Wherein, vaFor the frequency of Anti-Stokes signal.
Stokes signal is determined according to the frequency range of the frequency range of Stokes signal and Anti-Stokes signal
Spectral line A31 and Anti-Stokes signal spectral line A32.
In the step s 120, Stokes signal is converted into the first current signal, and Anti-Stokes signal is converted
For the second current signal.
Specifically, due in time domain Stokes signal and Anti-Stokes signal be equal optical signal, for the ease of root
Calculate temperature value according to the parameter of the Stokes signal and Anti-Stokes signal, can by photoelectric conversion mode will it is described this
Lentor signal energy transmitting electron makes its movement form the first current signal, by Anti-Stokes signal energy transmission to electricity
Son makes its movement form the second current signal.
In step s 130, the first current signal and the second current signal are debugged, so that the first current signal
Waveform reaches default coincidence threshold value be overlapped degree of the waveform of the second current signal in predetermined time section.
Specifically, include thermally sensitive frequency content and temperature-resistant frequency in the first current signal at
Point；It equally include thermally sensitive frequency content and temperature-resistant frequency content in second current signal.
As shown in figure 3, the corresponding waveform of thermally sensitive frequency content in the first current signal is S31, to temperature
The corresponding waveform of insensitive frequency content is N31；The corresponding waveform of thermally sensitive frequency content in second current signal
For S32, the corresponding waveform of temperature-resistant frequency content is N32.
It is worth noting that, the predetermined time is interregional to be thermally sensitive frequency in first current signal
The corresponding time interval of rate ingredient is also possible to the corresponding time zone of thermally sensitive frequency content in the second current signal
Between, due to, Stokes signal and Anti-Stokes signal are the heterogeneity in Raman signal, so, first electric current letter
Thermally sensitive frequency content pair in the corresponding time interval of thermally sensitive frequency content and the second current signal in number
The time interval answered is consistent, such as the section Δ t1 in Fig. 3.In addition, the predetermined time section can also be to belong to the section Δ t1
Some interior subinterval.
Specifically, as shown in figure 3, in predetermined time interval Δ t1, judge the waveform and the second electricity of the first current signal
Flow whether coincidence degree of the waveform of signal in predetermined time interval Δ t1 reaches default coincidence threshold value.
If the waveform of the first current signal being overlapped in predetermined time interval Δ t1 with the waveform of the second current signal
Degree does not reach default and is overlapped threshold value, can be by constantly regulate the photoelectric parameter of the first current signal and the second current signal
(such as bias voltage, gain etc.), the waveform of the waveform and the second current signal that finally make first current signal is in pre- timing
Between coincidence degree in interval Δ t1 reach default and be overlapped threshold value.
Wherein, in predetermined time interval Δ t1, the waveform and second of the coincidence degree and first current signal
The distance between waveform of current signal correlation, between the waveform of first current signal and the waveform of the second current signal
Distance is bigger, and the coincidence degree of the waveform of the waveform of the first current signal and the second current signal is smaller；The first electric current letter
Number waveform and the distance between the waveform of the second current signal it is smaller, the waveform of the first current signal and the second current signal
The coincidence degree of waveform is bigger.
In step S140, the luminous flux and second of default the first current signal for being overlapped threshold value is reached according to coincidence degree
The luminous flux of current signal calculates temperature value.
Further, described " to reach the luminous flux and the of default the first current signal for being overlapped threshold value according to coincidence degree
The luminous flux of two current signals calculates temperature value " include:
Reach default be overlapped in threshold portion waveform in the coincidence degree and presets reference position；According to the reference
The temperature value of position, the corresponding Stokes signal luminous flux in the reference position and Anti-Stokes signal luminous flux and to be measured
The corresponding Stokes signal luminous flux in position and Anti-Stokes signal luminous flux calculate the temperature value of the position to be measured.
Further, it is calculated by the following formula the temperature value of position to be measured:
Wherein, T is the temperature value of position to be measured, T0For the temperature value of reference position, h is Planck's constant, k be Bohr hereby
Graceful constant, Δ v are the Raman Phonon frequency of optical fiber, are a constant, ΦSIt (T) is the corresponding Stokes letter in the position to be measured
Number luminous flux, ΦASIt (T) is the corresponding Anti-Stokes signal luminous flux in the position to be measured, ΦS(T0) it is the reference position
Corresponding Stokes signal luminous flux, ΦAS(T0) it is the corresponding Anti-Stokes signal luminous flux in the reference position.
Specifically, the waveform S31 and the second current signal of the thermally sensitive frequency content in the first current signal
In the waveform S32 of thermally sensitive frequency content be overlapped degree in predetermined time interval Δ t1 and reach default and be overlapped threshold value
When, it calculates the first current signal and is overlapped when degree reaches default coincidence threshold value in corresponding predetermined time interval Δ t1 in waveform
The luminous flux phi of reference positionS(T0):
ΦS(T0)=KS·S·vs 4·Φe·exp[-(α0+αs)·L]·RS(T0)
Wherein, KSFor the related coefficient in stokes scattering section, S is the backscattering factor of optical fiber, vsFor Stokes
The frequency of signal, ΦeIn the photon flux of the laser signal of optical fiber incidence end, α0For the fiber transmission attenuation of incident laser signal,
αsFor the fiber transmission attenuation of Stokes signal, L is wave reference position, RS(T0) be and the population in optical fiber molecule low-lying level
The related coefficient of number, it is related with the temperature at optical fiber L.
RS(T0)=[1-exp (- h Δ v/kT0)]-1
Wherein, h is Planck's constant, and Δ v is the vibration frequency of optical fiber phonon, and k is Boltzmann constant, T0For reference bit
The temperature value set.
It calculates the second current signal and is overlapped corresponding predetermined time interval Δ when degree reaches default coincidence threshold value in waveform
The luminous flux phi of reference position in t1AS(T0):
ΦAS(T0)=KAS·S·va 4·Φe·exp[-(α0+αa)·L]·RAS(T0)
Wherein, KASFor coefficient related with anti-Stokes scattering section, vaFor the frequency of Anti-Stokes signal, αaFor
The fiber transmission attenuation of Anti-Stokes signal, RAS(T0) it is coefficient related with the i on population on optical fiber molecule high level, with
Temperature at optical fiber L is related.
RAS(T0)=[exp (h Δ v/kT0)-1]-1
Equally, the first current signal can be calculated in the above manner in the corresponding luminous flux phi in position to be measuredS(T) and second
Current signal is in the corresponding luminous flux phi in position to be measuredAS(T), in the present embodiment, the luminous flux is the width of the current signal
Value.
Embodiment 2
Fig. 4 shows a kind of flow diagram of fiber temperature signal demodulation method of second embodiment of the invention offer.
The fiber temperature signal demodulation method includes the following steps:
In step S210, Raman diffused light subsignal is obtained.
The step is identical as step S110, and details are not described herein.
In step S220, carried out at denoising according to Raman diffused light subsignal of the noise model pre-established to acquisition
Reason.
Specifically, due to the Stokes signal and Anti-Stokes signal of carrying temperature value in Raman diffused light subsignal
It is very faint, wherein containing a large amount of noise, therefore, how to realize and faint Raman diffused light subsignal is carried out effectively
Processing becomes a very important part during measurement fiber optic temperature.
In the present embodiment, can the noise in advance in the Raman diffused light subsignal to predetermined number learn, will be multiple
Noise obtains the corresponding noise model of Raman diffused light subsignal by way of fitting.
In some other embodiments, noise model can also be pre-established, which can be neural network
Model, by the Raman diffused light subsignal and the corresponding noise of Raman diffused light subsignal of predetermined number to the neural network
Model is trained, and obtains trained noise model.
Specifically, after obtaining noise model, i.e., the regularity of distribution of the noise in Raman diffused light subsignal.It is made an uproar by this
Raman diffused light subsignal of the acoustic model by the noise filtering in Raman diffused light subsignal, after being denoised.
In step S230, Stokes signal is converted into the first current signal, and Anti-Stokes signal is converted
For the second current signal.
This step is identical as step S120, and details are not described herein.
In step S240, the first current signal and the second current signal are debugged, so that the first current signal
Waveform reaches default coincidence threshold value be overlapped degree of the waveform of the second current signal in predetermined time section.
This step is identical as step S130, and details are not described herein.
In step s 250, the luminous flux and second of default the first current signal for being overlapped threshold value is reached according to coincidence degree
The luminous flux of current signal calculates temperature value.
This step is identical as step S140, and details are not described herein.
Embodiment 3
Fig. 5 shows a kind of flow diagram of fiber temperature signal demodulation method of third embodiment of the invention offer.
The fiber temperature signal demodulation method includes the following steps:
In step s310, Raman diffused light subsignal is obtained.
This step is identical as step S110, and details are not described herein.
In step s 320, Stokes signal is converted into the first current signal, and Anti-Stokes signal is converted
For the second current signal.
This step is identical as step S120, and details are not described herein.
In step S330, the first current signal and the second current signal are debugged, so that the first current signal
Waveform reaches default coincidence threshold value be overlapped degree of the waveform of the second current signal in predetermined time section.
This step is identical as step S130, and details are not described herein.
In step S340, the luminous flux and second of default the first current signal for being overlapped threshold value is reached according to coincidence degree
The luminous flux of current signal calculates temperature value.
This step is identical as step S140, and details are not described herein.
In step S350, position to be measured is determined according to the launch time of laser and the arrival time for reflecting signal.
Further, " determining position to be measured according to the launch time of laser and the arrival time for reflecting signal " includes:
Determine the laser signal described according to the time that reaches of the launch time of the laser and the reflection signal
Propagation distance in optical fiber；The position to be measured is determined according to the transmitting position of the laser and the propagation distance.
It specifically, can also be according to launch time of laser signal, reflection interval after temperature value is obtained by calculation
The spread speed of arrival time and laser signal in a vacuum calculates the corresponding location information of the temperature value.
For example, the propagation distance Len of laser signal in a fiber can be calculate by the following formula:
Wherein, C is the spread speed of laser signal in a vacuum, and t is that laser signal is emitted between reflection signal reception
Time interval, i.e., the absolute value of the receiving time time difference of the launch time of laser signal and reflection signal, n are optical fiber
Refractive index.
Position at the Len apart from laser signal transmitting position (emitting the position where the device of the laser signal)
It sets as position to be measured.
Embodiment 4
Fig. 6 shows a kind of structural schematic diagram of fiber temperature signal demodulating equipment of fourth embodiment of the invention offer.
The fiber temperature signal demodulating equipment 400 includes obtaining module 410, conversion module 420, debugging module 430 and meter
Calculate module 440.
Module 410 is obtained, for obtaining Raman diffused light subsignal, wherein the Raman diffused light subsignal includes this
Lentor signal and Anti-Stokes signal.
Conversion module 420, for the Stokes signal to be converted to the first current signal, and by the anti-stoke
This signal is converted to the second current signal.
Debugging module 430, for being debugged to first current signal and second current signal, so that described
The waveform of first current signal reaches default be overlapped degree of the waveform of second current signal in predetermined time section
It is overlapped threshold value.
Computing module 440, for reached according to coincidence degree default the first current signal for being overlapped threshold value luminous flux and
The luminous flux of second current signal calculates temperature value.
The present invention also provides a kind of fiber optic temperature (FBG) demodulator, the fiber optic temperature (FBG) demodulator includes memory and processing
Device, memory can mainly include storing program area and storage data area, wherein storing program area can storage program area, at least
Application program needed for one function etc.；Storage data area, which can be stored, uses created data etc. according to mobile phone.In addition, depositing
Reservoir may include high-speed random access memory, can also include nonvolatile memory, for example, at least a disk storage
Device, flush memory device or other volatile solid-state parts.
The processor is for running the computer program stored in the memory so that the fiber optic temperature demodulates
Instrument executes the function of each module in fiber temperature signal demodulation method or fiber temperature signal demodulating equipment in the above embodiments
Energy.
Optionally, processor may include one or more processing units；Preferably, processor can integrate application processor,
The main processing operation system of application processor, user interface and application program etc..Processor can integrate modem processor,
Modem processor can not also be integrated into the processor.
Fiber optic temperature is demodulated it will be understood by those skilled in the art that above-mentioned fiber optic temperature (FBG) demodulator structure is not constituted
The restriction of instrument may include perhaps combining certain components or different component layouts than illustrating more or fewer components.
The present embodiment additionally provides a kind of computer storage medium, for storing used in above-mentioned fiber optic temperature (FBG) demodulator
The computer program.
In several embodiments provided herein, it should be understood that disclosed device and method can also pass through
Other modes are realized.The apparatus embodiments described above are merely exemplary, for example, flow chart and structure in attached drawing
Figure shows the system frame in the cards of the device of multiple embodiments according to the present invention, method and computer program product
Structure, function and operation.In this regard, each box in flowchart or block diagram can represent a module, section or code
A part, a part of the module, section or code includes one or more for implementing the specified logical function
Executable instruction.It should also be noted that function marked in the box can also be to be different from the implementation as replacement
The sequence marked in attached drawing occurs.For example, two continuous boxes can actually be basically executed in parallel, they are sometimes
It can execute in the opposite order, this depends on the function involved.It is also noted that in structure chart and/or flow chart
The combination of each box and the box in structure chart and/or flow chart, can function or movement as defined in executing it is dedicated
Hardware based system realize, or can realize using a combination of dedicated hardware and computer instructions.
In addition, each functional module or unit in each embodiment of the present invention can integrate one independence of formation together
Part, be also possible to modules individualism, an independent part can also be integrated to form with two or more modules.
If the function is realized and when sold or used as an independent product in the form of software function module, can store one
In a computer-readable storage medium.Based on this understanding, technical solution of the present invention is substantially in other words to existing skill
The part of part or the technical solution that art contributes can be embodied in the form of software products, the computer software
Product is stored in a storage medium, including some instructions are used so that computer equipment (can be smart phone, a
People's computer, server or network equipment etc.) it performs all or part of the steps of the method described in the various embodiments of the present invention.
And storage medium above-mentioned includes: that USB flash disk, mobile hard disk, read-only memory (ROM, Read-Only Memory), arbitrary access are deposited
The various media that can store program code such as reservoir (RAM, Random Access Memory), magnetic or disk.
The above description is merely a specific embodiment, but scope of protection of the present invention is not limited thereto, any
Those familiar with the art in the technical scope disclosed by the present invention, can easily think of the change or the replacement, and should all contain
Lid is within protection scope of the present invention.
Claims (10)
1. a kind of fiber temperature signal demodulation method, which is characterized in that this method comprises:
Obtain Raman diffused light subsignal, wherein the Raman diffused light subsignal includes Stokes signal and anti-stoke
This signal；
The Stokes signal is converted into the first current signal, and the Anti-Stokes signal is converted into the second electric current
Signal；
First current signal and second current signal are debugged so that the waveform of first current signal with
Coincidence degree of the waveform of second current signal in predetermined time section reaches default and is overlapped threshold value；
Reach the luminous flux of default the first current signal for being overlapped threshold value and the luminous flux of the second current signal according to coincidence degree
Calculate temperature value.
2. fiber temperature signal demodulation method according to claim 1, which is characterized in that described " to obtain Raman diffused light
Subsignal " includes:
Obtain the reflection signal of laser in a fiber；
Reflection signal progress Fourier transformation is obtained into the reflection signal of frequency domain, according to the frequency spectrum of Raman diffused light subsignal
Range obtains Raman diffused light subsignal in the reflection signal of the frequency domain.
3. fiber temperature signal demodulation method according to claim 1, which is characterized in that described " to obtain Raman diffused light
After subsignal " further include:
Denoising is carried out according to the Raman diffused light subsignal of the noise model pre-established to acquisition, wherein described
Noise model is learnt to obtain by the noise in the Raman diffused light subsignal to preset predetermined number.
4. fiber temperature signal demodulation method according to claim 1, which is characterized in that described " to be reached according to coincidence degree
The luminous flux of luminous flux and the second current signal to default the first current signal for being overlapped threshold value calculates temperature value " include:
Reach default be overlapped in threshold portion waveform in the coincidence degree and presets reference position；
According to the temperature value of the reference position, the corresponding Stokes signal luminous flux in the reference position and anti-Stokes
The corresponding Stokes signal luminous flux of signal luminous flux and position to be measured and Anti-Stokes signal luminous flux calculate it is described to
The temperature value that location is set.
5. fiber temperature signal demodulation method according to claim 4, which is characterized in that be calculated by the following formula described
The temperature value of position to be measured:
Wherein, T is the temperature value of position to be measured, T0For the temperature value of reference position, h is Planck's constant, and k is that Boltzmann is normal
Number, Δ v are the Raman Phonon frequency of optical fiber, are a constant, ΦSIt (T) is the corresponding Stokes signal light in the position to be measured
Flux, ΦASIt (T) is the corresponding Anti-Stokes signal luminous flux in the position to be measured, ΦS(T0) it is that the reference position is corresponding
Stokes signal luminous flux, ΦAS(T0) it is the corresponding Anti-Stokes signal luminous flux in the reference position.
6. fiber temperature signal demodulation method according to claim 4, which is characterized in that further include:
The position to be measured is determined according to the time that reaches of the launch time of the laser and the reflection signal.
7. fiber temperature signal demodulation method according to claim 6, which is characterized in that described " according to the laser
Launch time and the time that reaches of the reflection signal determine the position to be measured " include:
Determine the laser signal in the optical fiber according to the time that reaches of the launch time of the laser and the reflection signal
In propagation distance；
The position to be measured is determined according to the transmitting position of the laser and the propagation distance.
8. a kind of fiber temperature signal demodulating equipment, which is characterized in that the device includes:
Module is obtained, for obtaining Raman diffused light subsignal, wherein the Raman diffused light subsignal includes Stokes letter
Number and Anti-Stokes signal；
Conversion module, for the Stokes signal to be converted to the first current signal, and by the Anti-Stokes signal
Be converted to the second current signal；
Debugging module, for being debugged to first current signal and second current signal, so that first electricity
The waveform for flowing signal reaches default coincidence threshold be overlapped degree of the waveform of second current signal in predetermined time section
Value；
Computing module, for reaching the luminous flux and the second electric current of default the first current signal for being overlapped threshold value according to coincidence degree
The luminous flux of signal calculates temperature value.
9. a kind of fiber optic temperature (FBG) demodulator, which is characterized in that the fiber optic temperature (FBG) demodulator includes memory and processor, institute
It states memory and runs the computer program so that the fiber optic temperature (FBG) demodulator for storing computer program, the processor
Perform claim requires 1 to 7 described in any item fiber temperature signal demodulation methods.
10. a kind of computer storage medium, which is characterized in that it stores institute in fiber optic temperature (FBG) demodulator described in claim 9
The computer program.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910307202.3A CN109900382B (en) | 2019-04-17 | 2019-04-17 | Optical fiber temperature signal demodulation method and device and optical fiber temperature demodulator |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910307202.3A CN109900382B (en) | 2019-04-17 | 2019-04-17 | Optical fiber temperature signal demodulation method and device and optical fiber temperature demodulator |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109900382A true CN109900382A (en) | 2019-06-18 |
| CN109900382B CN109900382B (en) | 2020-12-04 |
Family
ID=66955775
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910307202.3A Active CN109900382B (en) | 2019-04-17 | 2019-04-17 | Optical fiber temperature signal demodulation method and device and optical fiber temperature demodulator |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109900382B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112964387A (en) * | 2021-03-02 | 2021-06-15 | 南京大学 | Demodulation method for temperature along optical fiber in Raman optical time domain reflectometer |
| CN113008407A (en) * | 2021-02-26 | 2021-06-22 | 中煤科工集团重庆研究院有限公司 | Multi-parameter distributed optical fiber signal demodulator |
| CN117290669A (en) * | 2023-11-24 | 2023-12-26 | 之江实验室 | Optical fiber temperature sensing signal noise reduction method, device and medium based on deep learning |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103256999A (en) * | 2013-05-10 | 2013-08-21 | 深圳市英唐电气技术有限公司 | Distributed type optical fiber temperature measuring method |
| CN104748887A (en) * | 2015-04-15 | 2015-07-01 | 中国电子科技集团公司第八研究所 | Distributed fiber temperature sensor and temperature triple-demodulation algorithm thereof |
| CN105953941A (en) * | 2016-04-29 | 2016-09-21 | 深圳艾瑞斯通技术有限公司 | Distributed fiber temperature measurement method and device based on Raman scattering |
| CN205909944U (en) * | 2016-06-02 | 2017-01-25 | 国家电网公司 | Hollow parallel reactor distributed optical fiber temperature measurement device of portable dry -type |
| CN108332877A (en) * | 2018-01-31 | 2018-07-27 | 北京航天控制仪器研究所 | A kind of the fiber grating measuring device and method of light pulse and optical time domain reflection combination |
| US20190086244A1 (en) * | 2017-09-19 | 2019-03-21 | Analog Devices, Inc. | Fiber bragg grating interrogation and sensing system and methods |
-
2019
- 2019-04-17 CN CN201910307202.3A patent/CN109900382B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103256999A (en) * | 2013-05-10 | 2013-08-21 | 深圳市英唐电气技术有限公司 | Distributed type optical fiber temperature measuring method |
| CN104748887A (en) * | 2015-04-15 | 2015-07-01 | 中国电子科技集团公司第八研究所 | Distributed fiber temperature sensor and temperature triple-demodulation algorithm thereof |
| CN105953941A (en) * | 2016-04-29 | 2016-09-21 | 深圳艾瑞斯通技术有限公司 | Distributed fiber temperature measurement method and device based on Raman scattering |
| CN205909944U (en) * | 2016-06-02 | 2017-01-25 | 国家电网公司 | Hollow parallel reactor distributed optical fiber temperature measurement device of portable dry -type |
| US20190086244A1 (en) * | 2017-09-19 | 2019-03-21 | Analog Devices, Inc. | Fiber bragg grating interrogation and sensing system and methods |
| CN108332877A (en) * | 2018-01-31 | 2018-07-27 | 北京航天控制仪器研究所 | A kind of the fiber grating measuring device and method of light pulse and optical time domain reflection combination |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113008407A (en) * | 2021-02-26 | 2021-06-22 | 中煤科工集团重庆研究院有限公司 | Multi-parameter distributed optical fiber signal demodulator |
| CN113008407B (en) * | 2021-02-26 | 2024-04-09 | 中煤科工集团重庆研究院有限公司 | Multi-parameter distributed optical fiber signal demodulator |
| CN112964387A (en) * | 2021-03-02 | 2021-06-15 | 南京大学 | Demodulation method for temperature along optical fiber in Raman optical time domain reflectometer |
| CN112964387B (en) * | 2021-03-02 | 2022-03-15 | 南京大学 | Demodulation method for temperature along optical fiber in Raman optical time domain reflectometer |
| CN117290669A (en) * | 2023-11-24 | 2023-12-26 | 之江实验室 | Optical fiber temperature sensing signal noise reduction method, device and medium based on deep learning |
| CN117290669B (en) * | 2023-11-24 | 2024-02-06 | 之江实验室 | Optical fiber temperature sensing signal noise reduction method, device and medium based on deep learning |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109900382B (en) | 2020-12-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109900382A (en) | Fiber temperature signal demodulation method, device and fiber optic temperature (FBG) demodulator | |
| US10444391B2 (en) | Noise removal for distributed acoustic sensing data | |
| US10401519B2 (en) | Noise removal for distributed acoustic sensing data | |
| CN104864979B (en) | A kind of modification method of distributed Raman fiber temp measuring system measurement error | |
| CN103115693B (en) | A kind of distributed fiber Raman temp measuring system | |
| CN108225541A (en) | The distributed fiberoptic sensor and foreign body intrusion signal for identifying foreign body intrusion perceive processing method | |
| CN109029769A (en) | High-precision temperature demodulation method based on distributed fiber Raman sensing technology | |
| CN102279062A (en) | Temperature calibration and demodulation method for distributed optical fiber sensing temperature measurement system | |
| CN104034409B (en) | Distributed optical fiber vibration sensing method and system based on pulse code external modulation | |
| CN109540207A (en) | A kind of calculation type distributing optical fiber sensing method and system | |
| Olsen III et al. | Observed one‐dimensional return stroke propagation speeds in the bottom 170 m of a rocket‐triggered lightning channel | |
| CN104006948B (en) | Based on the method that multimodal division cycle demodulates polarization maintaining optical fibre polarization coupled point position | |
| CN103940501B (en) | A kind of BOTDA distributed vibration sensing system based on dynamic phasing demodulation | |
| CN109084883A (en) | Based on phase-BOTDR optical fiber distributed type Brillouin's vibrating sensing measurement method | |
| CN106643842A (en) | Distributed sensor and sensing method | |
| CN107478250B (en) | A kind of Φ-OTDR signal demodulating method based on three port phase demodulations | |
| Dean et al. | Distributed vibration sensing for seismic acquisition | |
| CN209605977U (en) | A kind of temperature-measuring system of distributed fibers containing auto-adaptive filter circuit | |
| CN111289142B (en) | Signal adaptive interference cancellation method for distributed optical fiber Raman temperature measurement system | |
| CN109444895B (en) | Vibration information positioning method for eliminating interference fading of distributed vibration sensor | |
| CN103257114B (en) | Underground gas detecting method based on optical fiber Bragg grating sensor network | |
| CN108692830A (en) | A kind of temperature-measuring system of distributed fibers | |
| Wang et al. | Luminosity characteristics of multiple dart leader/return stroke sequences measured with a high‐speed digital image system | |
| CN113834509A (en) | Weak reflection fiber grating string demodulator based on deep learning | |
| Hou et al. | Signal processing of single-mode fiber sensor system based on Raman scattering |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20210111 Address after: 510700 room 108, building S2, No.3, Fukang West Street, Yushu Industrial Park, Science City, Huangpu District, Guangzhou City, Guangdong Province Patentee after: Guangzhou Weitong Technology Co.,Ltd. Address before: 510000 402, building C1, 182 Kexue Avenue, Luogang District, Guangzhou City, Guangdong Province Patentee before: GUANGZHOU SHENG'AN MEASUREMENT AND CONTROL SCIENCE & TECHNOLOGY Co.,Ltd. |
|
| CP03 | Change of name, title or address | ||
| CP03 | Change of name, title or address |
Address after: Room 404, Building D, No. 7 Jingye Third Street, Yushu Industrial Park, Science City, Huangpu District, Guangzhou City, Guangdong Province, 510000 Patentee after: Guangzhou Weitong Shengan Technology Co.,Ltd. Address before: 510700 room 108, building S2, No.3, Fukang West Street, Yushu Industrial Park, Science City, Huangpu District, Guangzhou City, Guangdong Province Patentee before: Guangzhou Weitong Technology Co.,Ltd. |