CN111307322A - Distributed optical fiber temperature sensing system based on annular light path - Google Patents

Distributed optical fiber temperature sensing system based on annular light path Download PDF

Info

Publication number
CN111307322A
CN111307322A CN201911248571.6A CN201911248571A CN111307322A CN 111307322 A CN111307322 A CN 111307322A CN 201911248571 A CN201911248571 A CN 201911248571A CN 111307322 A CN111307322 A CN 111307322A
Authority
CN
China
Prior art keywords
optical fiber
stokes
scattering
scattered light
rayleigh
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
Application number
CN201911248571.6A
Other languages
Chinese (zh)
Other versions
CN111307322B (en
Inventor
潘亮
杨潇君
杜瑞
朱晓非
李娜
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Aerospace Wanda Hi Tech Ltd
Original Assignee
Beijing Aerospace Wanda Hi Tech Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beijing Aerospace Wanda Hi Tech Ltd filed Critical Beijing Aerospace Wanda Hi Tech Ltd
Priority to CN201911248571.6A priority Critical patent/CN111307322B/en
Publication of CN111307322A publication Critical patent/CN111307322A/en
Application granted granted Critical
Publication of CN111307322B publication Critical patent/CN111307322B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01KMEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
    • G01K11/00Measuring temperature based upon physical or chemical changes not covered by groups G01K3/00, G01K5/00, G01K7/00 or G01K9/00
    • G01K11/32Measuring temperature based upon physical or chemical changes not covered by groups G01K3/00, G01K5/00, G01K7/00 or G01K9/00 using changes in transmittance, scattering or luminescence in optical fibres
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01KMEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
    • G01K11/00Measuring temperature based upon physical or chemical changes not covered by groups G01K3/00, G01K5/00, G01K7/00 or G01K9/00
    • G01K11/32Measuring temperature based upon physical or chemical changes not covered by groups G01K3/00, G01K5/00, G01K7/00 or G01K9/00 using changes in transmittance, scattering or luminescence in optical fibres
    • G01K11/324Measuring temperature based upon physical or chemical changes not covered by groups G01K3/00, G01K5/00, G01K7/00 or G01K9/00 using changes in transmittance, scattering or luminescence in optical fibres using Raman scattering

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Measuring Temperature Or Quantity Of Heat (AREA)

Abstract

The invention discloses a distributed optical fiber temperature sensing system based on an annular light path, which comprises: the device comprises a pulse laser (1), a wavelength division multiplexer (2), an optical switch (3), a first optical fiber ring (5), a second optical fiber ring (6), a temperature measurement optical fiber (7), a photoelectric detector (8), a data acquisition card (9) and an upper computer (10); the distributed optical fiber temperature sensing system of the annular optical path can perform laser power feedback adjustment and Rayleigh scattering noise compensation for inhibiting nonlinear effects. The invention improves the long-term stability of the distributed optical fiber temperature sensor in the practical application environment; the influence of the Rayleigh scattering light on the temperature demodulation precision is eliminated, so that the direct demodulation of the temperature field information of the sensing optical fiber is realized, the complexity of system debugging is reduced, and the temperature measurement precision and the long-term measurement stability are improved.

Description

Distributed optical fiber temperature sensing system based on annular light path
Technical Field
The invention belongs to the technical field of optical fiber sensing, and particularly relates to a distributed optical fiber temperature sensing system based on an annular light path.
Background
A distributed optical fiber temperature sensor based on Raman scattering utilizes the principle that Raman scattering light intensity and temperature in a sensing optical fiber are sensitive, and collects backward Raman scattering light of each point on the sensing optical fiber, wherein each point is modulated by temperature, by injecting detection light pulses into the sensing optical fiber. The Raman scattering light comprises an up-shifted anti-Stokes scattering light and a down-shifted Stokes scattering light, and because the light intensity is very weak, the signal to noise ratio of a signal is improved by a high-sensitivity avalanche photodetector and accumulation average processing. The anti-Stokes scattered light is sensitive to temperature and is used for temperature demodulation, and the Stokes scattered light is used for correcting the anti-Stokes scattered light and weakening the influence of light source power fluctuation and sensing optical fiber loss on temperature demodulation performance.
When the Raman is subjected to weak light detection, the signal-to-noise ratio of Raman scattering light can be improved by improving the output power of the laser, but the output power of the laser cannot be infinitely increased in order to avoid the occurrence of stimulated Raman scattering. In practical engineering use, the laser power changes with time, different losses of connecting devices and other factors. Therefore, to maximize the raman scattering signal-to-noise ratio, the output power of the laser needs to be adjusted in real time.
In addition, in practical engineering applications, the performance of the sensor is affected by environmental factors, such as performance changes of the light source and the photodetector due to changes in operating temperature, and nonlinear loss caused by bending or stretching of the sensing optical fiber or erosion of a strong corrosive substance, which all degrade the performance of the sensor, increase temperature demodulation errors, and fail in severe cases.
The external factors acting on the sensing fiber have uncertainty, and the loss degree of the sensing fiber is related to the wavelength, so that the external factors cannot be eliminated or corrected. The annular sensing light path is adopted to collect Raman scattered light of the sensing optical fiber in the positive direction and the negative direction, and optical fiber loss related to the position can be effectively eliminated through geometric average processing, so that the environmental adaptability of the sensor is greatly improved.
The scattered light generated in the optical fiber includes not only raman scattered light but also rayleigh scattered light and brillouin scattered light, and a common method is to select the raman scattered light by using a wavelength division multiplexer. However, since the rayleigh scattered light intensity is high and is about 30dB higher than the raman scattered light intensity, and the isolation of a commercially available wavelength division multiplexer is usually 35dB to 40dB, a large amount of rayleigh scattered light noise exists in the measured raman scattered light. In a small dynamic range (< 100 ℃), rayleigh scattered light has less influence on temperature demodulation errors, but temperature errors gradually deteriorate as the dynamic range increases. Therefore, it is necessary to compensate for rayleigh scattered light noise without changing the device performance.
Disclosure of Invention
The technical problem solved by the invention is as follows: the distributed optical fiber temperature sensing system based on the annular light path is provided, through the annular light path structure, nonlinear loss caused by factors such as bending and corrosion of strong corrosive substances in the optical fiber to Raman scattering light intensity is eliminated, long-term stability of the distributed optical fiber temperature sensor in an actual application environment is improved, and the application field of the distributed optical fiber temperature sensing system is expanded. And (4) passing. Meanwhile, Rayleigh scattered light in the back Raman scattered light is calculated, and the influence of the Rayleigh scattered light on temperature demodulation precision is eliminated, so that the temperature field information of the sensing optical fiber is directly demodulated, the complexity of system debugging is reduced, and the temperature measurement precision and the long-term measurement stability are improved.
The purpose of the invention is realized by the following technical scheme: a distributed optical fiber temperature sensing system based on a ring-shaped optical path comprises: the device comprises a pulse laser, a wavelength division multiplexer, an optical switch, a first optical fiber ring, a second optical fiber ring, a temperature measuring optical fiber, a photoelectric detector, a data acquisition card and an upper computer; the output ends of the pulse laser and the pulse laser are connected with the input end of the wavelength division multiplexer; the common end of the wavelength division multiplexer is connected with the input end of the optical switch; two output ends of the optical switch are respectively connected with the first optical fiber ring and the second optical fiber ring; the output ends of the first optical fiber ring and the second optical fiber ring are respectively connected with two ends of the temperature measuring optical fiber; two output ends of the wavelength division multiplexer are respectively connected with two input ends of the photoelectric detector; two output ends of the photoelectric detector are respectively connected with two input ends of a data acquisition card, and the output end of the data acquisition card is connected with the input end of an upper computer; the first optical fiber ring and the second optical fiber ring are both placed in a temperature measuring device, and the output end of the temperature measuring device is connected with the input end of an upper computer; the distributed optical fiber temperature sensing system of the annular optical path can perform laser power feedback adjustment and Rayleigh scattering noise compensation for inhibiting nonlinear effects.
In the above distributed optical fiber temperature sensing system based on the annular optical path, the laser power feedback adjustment for suppressing the occurrence of the nonlinear effect includes the following steps: the method comprises the following steps: two groups of anti-Stokes scattered light V are acquired by a data acquisition cardas1And Vas2And Stokes scattered light Vs1And Vs2Respectively carrying out geometric mean calculation to obtain VasmAnd Vsm(ii) a Step two: calculating the ratio of the anti-Stokes scattered light to the Stokes scattered light after geometric averaging, VR=Vasm/VsmAnd calculating the ratio difference DeltaV of the positions of the first optical fiber ring and the second optical fiber ringR=VR(L0)-VR(L1) Wherein L is0And L1The distances between the positions of the first optical fiber ring and the second optical fiber ring are respectively; step three: starting from 0mA, the pumping current of the pulse laser is gradually increased, and the delta V is observed in real timeRWhen Δ VR>10%×MIN[VR(L0),VR(L1)]Reducing the pumping current of the pulse laser to delta V<1%×MIN[VR(L0),VR(L1)]。
In the above distributed optical fiber temperature sensing system based on the annular optical path, the rayleigh scattering noise compensation includes the following steps: (1) placing the temperature measuring optical fiber, the first optical fiber ring and the second optical fiber ring in the same temperature environment, and keeping the environment temperature stable; acquiring anti-Stokes scattering signals and Stokes scattering signals from forward and backward directions respectively to obtain forward acquired anti-Stokes signals Vas_forWith Stokes scattering signal Vs_forReversely collected anti-stokes signal Vas_backWith Stokes scattering signal Vs_back(ii) a (2) Anti-stokes based on forward acquisitionObtaining optical fiber L from the scattering signals of stokes and stokesnForward anti-stokes scattering signal loss at
Figure BDA0002308382310000031
Loss of signal associated with forward stokes scattering
Figure BDA0002308382310000032
And obtaining compensated forward anti-Stokes scattering signal loss
Figure BDA0002308382310000033
With compensated forward stokes scattering signal loss
Figure BDA0002308382310000034
(3) Obtaining the optical fiber L according to the reversely collected anti-Stokes and Stokes scattering signalsnLoss of reverse anti-stokes scattering signal
Figure BDA0002308382310000035
Loss of signal associated with reverse Stokes scattering
Figure BDA0002308382310000036
And obtaining compensated reverse anti-Stokes scattering signal loss
Figure BDA0002308382310000037
With compensated reverse Stokes scattering signal loss
Figure BDA0002308382310000038
(4) Based on compensated forward anti-Stokes scattering signal loss
Figure BDA0002308382310000039
With compensated forward stokes scattering signal loss
Figure BDA00023083823100000310
Obtaining compensated forward anti-Stokes scattering signal loss containing Rayleigh scattering noise
Figure BDA00023083823100000311
With compensated forward stokes scattering signal loss
Figure BDA00023083823100000312
Based on compensated reverse anti-Stokes scattering signal loss
Figure BDA00023083823100000313
With compensated reverse Stokes scattering signal loss
Figure BDA0002308382310000041
Obtaining compensated reverse anti-stokes scattering signal loss containing Rayleigh scattering noise
Figure BDA0002308382310000042
With compensated reverse Stokes scattering signal loss
Figure BDA0002308382310000043
In the distributed optical fiber temperature sensing system based on the annular light path, in the step (2), the optical fiber LnOptical fiber L is obtained from anti-Stokes and Stokes scattering signals of forward collectionnForward anti-stokes scattering signal loss at
Figure BDA0002308382310000044
Loss of signal associated with forward stokes scattering
Figure BDA0002308382310000045
Respectively expressed as:
Figure BDA0002308382310000046
Figure BDA0002308382310000047
wherein m is the number of translation steps of the loss point of the optical fiber.
In the distributed optical fiber temperature sensing system based on the annular optical path, in the step (2), the compensated forward anti-stokes scattering signal loss
Figure BDA0002308382310000048
With compensated forward stokes scattering signal loss
Figure BDA0002308382310000049
Respectively expressed as:
Figure BDA00023083823100000410
Figure BDA00023083823100000411
in the above distributed optical fiber temperature sensing system based on the annular optical path, in step (3), the optical fiber LnLoss of reverse anti-stokes scattering signal
Figure BDA00023083823100000412
Loss of signal associated with reverse Stokes scattering
Figure BDA00023083823100000413
Respectively expressed as:
Figure BDA00023083823100000414
Figure BDA00023083823100000415
in the distributed optical fiber temperature sensing system based on the annular optical path, in the step (3), the compensated reverse anti-stokes scattering signal loss
Figure BDA00023083823100000416
And after compensationLoss of reverse stokes scattered signal
Figure BDA00023083823100000417
Respectively expressed as:
Figure BDA0002308382310000051
Figure BDA0002308382310000052
in the above distributed optical fiber temperature sensing system based on the annular optical path, in step (4), the compensated forward anti-stokes scattering signal loss containing rayleigh scattering noise
Figure BDA0002308382310000053
Comprises the following steps:
Figure BDA0002308382310000054
wherein l is the position distance of a certain point of the optical fiber, h is Planck constant, Deltav is Raman scattering frequency shift in the quartz optical fiber, k is Boltzmann constant, T is the environment temperature of the optical fiber, αasFor the attenuation coefficient of anti-Stokes scattered light in silica fiber, αsThe attenuation coefficient of Stokes' scattered light in a silica optical fiber, αrAttenuation coefficient of Rayleigh scattered light in quartz fiber, AforScattering coefficient of anti-Stokes scattered light obtained for forward detection, BforDetecting Rayleigh scattering coefficient, C, in channel for anti-Stokes scattered lightforScattering coefficient of anti-Stokes scattered light, D, obtained for forward detectionforRayleigh scattering coefficient, A, in a channel for Stokes scattering light detectionbackScattering coefficient of anti-Stokes scattered light obtained for reverse detection, BbackDetecting Rayleigh scattering coefficient, C, in channel for anti-Stokes scattered lightbackStokes scattering coefficient, D, for reverse detectionbackIs reversedRayleigh scattering coefficient in the detection channel of Stokes scattered light to be detected, αr0The attenuation coefficient of Rayleigh scattered light in a forward anti-Stokes detection signal in a quartz optical fiber is detected.
In the above distributed optical fiber temperature sensing system based on the ring optical path, in step (4), the compensated forward stokes scattering signal loss containing rayleigh scattering noise
Figure BDA0002308382310000055
Comprises the following steps:
Figure BDA0002308382310000056
wherein l is the position distance of a certain point of the optical fiber, h is Planck constant, Deltav is Raman scattering frequency shift in the quartz optical fiber, k is Boltzmann constant, T is the environment temperature of the optical fiber, αasFor the attenuation coefficient of anti-Stokes scattered light in silica fiber, αsThe attenuation coefficient of Stokes' scattered light in a silica optical fiber, αrAttenuation coefficient of Rayleigh scattered light in quartz fiber, CforScattering coefficient of anti-Stokes scattered light, D, obtained for forward detectionforRayleigh scattering coefficient in the detection channel for Stokes scattered light, αr1The attenuation coefficient of Rayleigh scattered light in a forward Stokes detection signal in a quartz optical fiber is detected.
In the above distributed optical fiber temperature sensing system based on the ring optical path, in step (4), the compensated reverse anti-stokes scattering signal loss containing rayleigh scattering noise
Figure BDA0002308382310000061
Comprises the following steps:
Figure BDA0002308382310000062
wherein l is the position distance of a certain point of the optical fiber, h is Planck constant, Deltav is Raman scattering frequency shift in the quartz optical fiber, k is Boltzmann constant,t is the ambient temperature of the fiber, αasFor the attenuation coefficient of anti-Stokes scattered light in silica fiber, αsThe attenuation coefficient of Stokes' scattered light in a silica optical fiber, αrAttenuation coefficient of Rayleigh scattered light in quartz fiber, AbackScattering coefficient of anti-Stokes scattered light obtained for reverse detection, BbackDetection of Rayleigh scattered light scattering coefficients in a channel for anti-Stokes scattered light, αr2The attenuation coefficient of Rayleigh scattered light in a signal in a quartz fiber is detected for reverse anti-Stokes.
In the above distributed optical fiber temperature sensing system based on the ring optical path, in step (4), the compensated reverse stokes scattering signal loss containing rayleigh scattering noise
Figure BDA0002308382310000063
Comprises the following steps:
Figure BDA0002308382310000064
wherein l is the position distance of a certain point of the optical fiber, h is Planck constant, Deltav is Raman scattering frequency shift in the quartz optical fiber, k is Boltzmann constant, T is the environment temperature of the optical fiber, αasFor the attenuation coefficient of anti-Stokes scattered light in silica fiber, αsThe attenuation coefficient of Stokes' scattered light in a silica optical fiber, αrAttenuation coefficient of Rayleigh scattered light in quartz fiber, CbackStokes scattering coefficient, D, for reverse detectionbackRayleigh scattering coefficient in the channel is detected for Stokes scattering light for reverse detection, αr3The attenuation coefficient of Rayleigh scattered light in a signal in a quartz fiber is detected for the reverse Stokes.
Compared with the prior art, the invention has the following beneficial effects:
(1) compared with fixed laser parameters, the pulsed light regulated by feedback control can ensure that the fiber-entering power of the laser is improved as much as possible on the premise of not generating stimulated Raman scattering, thereby ensuring that the signal-to-noise ratio of Raman scattering light intensity is maximized.
(2) Compared with the temperature measurement result without a Rayleigh scattering noise filtering device or a Rayleigh scattering noise compensation algorithm, the Rayleigh scattering noise compensation algorithm can effectively make up the Rayleigh scattering noise influence caused by insufficient isolation of the wavelength division multiplexer. By changing the temperature of the fiber ring, the rayleigh scattering intensity distribution in the sensing fiber can be quickly calculated and subtracted from the demodulation formula.
Drawings
Various other advantages and benefits will become apparent to those of ordinary skill in the art upon reading the following detailed description of the preferred embodiments. The drawings are only for purposes of illustrating the preferred embodiments and are not to be construed as limiting the invention. Also, like reference numerals are used to refer to like parts throughout the drawings. In the drawings:
fig. 1 is a block diagram of a distributed optical fiber temperature sensing system of a ring optical path according to an embodiment of the present invention.
Detailed Description
Exemplary embodiments of the present disclosure will be described in more detail below with reference to the accompanying drawings. While exemplary embodiments of the present disclosure are shown in the drawings, it should be understood that the present disclosure may be embodied in various forms and should not be limited to the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the scope of the disclosure to those skilled in the art. It should be noted that the embodiments and features of the embodiments may be combined with each other without conflict. The present invention will be described in detail below with reference to the embodiments with reference to the attached drawings.
Fig. 1 is a block diagram of a distributed optical fiber temperature sensing system of a ring optical path according to an embodiment of the present invention. As shown in fig. 1, the distributed optical fiber temperature sensing system of the ring-shaped optical path includes: the device comprises a pulse laser 1, a wavelength division multiplexer 2, an optical switch 3, a first optical fiber ring 5, a second optical fiber ring 5, a temperature measuring optical fiber 7, a photoelectric detector 8, a data acquisition card 9 and an upper computer 10; wherein,
the output ends of the pulse laser 1 and the pulse laser 1 are connected with the input end of the wavelength division multiplexer 2; the common end of the wavelength division multiplexer 2 is connected with the input end of the optical switch 3; two output ends of the optical switch 3 are respectively connected with a first optical fiber ring 5 and a second optical fiber ring 6; the output ends of the first optical fiber ring 5 and the second optical fiber ring 6 are respectively connected with two ends of a temperature measuring optical fiber 7; two output ends of the wavelength division multiplexer 2 are respectively connected with two input ends of the photoelectric detector 8; two output ends of the photoelectric detector 8 are respectively connected with two input ends of a data acquisition card 9, and the output end of the data acquisition card is connected with the input end of an upper computer 10; the first optical fiber ring 5 and the second optical fiber ring 6 are both placed in the temperature measuring device 4, and the output end of the temperature measuring device 4 is connected with the input end of the upper computer 10;
the distributed optical fiber temperature sensing system of the annular optical path can perform laser power feedback adjustment and Rayleigh scattering noise compensation for inhibiting nonlinear effects.
Distributed optical fiber temperature sensing system based on annular light path includes: the device comprises a pulse laser with the center wavelength of 1550nm, a band-pass wavelength division multiplexer with the center wavelengths of 1450nm and 1661nm, a 1 multiplied by 2 optical switch, an optical fiber ring, a multimode temperature measuring optical fiber with the length of 2km, two optical fiber rings with the length of 100m, a two-channel synchronous acquisition photoelectric detector with the bandwidth of 100MHz, a two-channel data acquisition card of 100MSPS, an upper computer and a PT1000 platinum resistance measuring device.
The laser power feedback regulation for inhibiting the nonlinear effect comprises the following steps:
the method comprises the following steps: two groups of anti-Stokes scattered light (V) are acquired by a data acquisition cardas1And Vas2) And Stokes scattered light (V)s1And Vs2) Respectively carrying out geometric mean calculation to obtain VasmAnd Vsm
Step two: calculating the ratio of the anti-Stokes scattered light to the Stokes scattered light after geometric averaging, VR=Vasm/VsmAnd calculating the ratio difference DeltaV of the positions of the first optical fiber ring and the second optical fiber ringR=VR(L0)-VR(L1) Wherein L is0=50m,L1=2150m;
Step three: starting from 20mA, the pumping current of the laser is gradually increased, and the delta V is observed in real timeRWhen Δ VR>10%×MIN[VR(L0),VR(L1)]Reducing the pumping current of the laser to delta V<1%×MIN[VR(L0),VR(L1)]。
This embodiment passes through two optical fiber rings at the annular light path head and end, realizes the real-time supervision to stimulated raman scattering threshold value, and the pumping current of feedback regulation laser instrument realizes injecting pulse optical power maximize under spontaneous raman scattering state, and then promotes sensing system's wholeness ability.
In the invention: the temperature sensing system has a Rayleigh scattering noise compensation function, and the compensation algorithm comprises the following steps:
the method comprises the following steps: and placing the temperature measuring optical fiber, the first optical fiber ring and the second optical fiber ring in an environment at 25 ℃, and keeping the environment temperature stable. Respectively collecting anti-Stokes scattering signals and Stokes scattering signals from positive and negative directions, Vas_for,Vs_for,Vas_backAnd Vs_back
Step two: and respectively calculating the loss of each point including the optical fiber connector and the optical fiber fusion point for the anti-Stokes and Stokes scattering signals collected in the forward direction. Optical fiber LnThe anti-stokes scattering signal loss and the stokes scattering signal loss at (a) can be expressed as:
Figure BDA0002308382310000091
Figure BDA0002308382310000092
wherein m is 5 m. The compensated raman scattering curve is:
Figure BDA0002308382310000093
Figure BDA0002308382310000094
wherein L is 2200 m.
Step three: and (3) for the reversely collected anti-Stokes and Stokes scattering signals, obtaining the loss of each point of the optical fiber calculated by the direction data by adopting the same processing method as the step two:
Figure BDA0002308382310000095
Figure BDA0002308382310000096
the compensated raman scattering curve should be:
Figure BDA0002308382310000097
Figure BDA0002308382310000098
step three: the two-way anti-stokes scattering signal and the stokes scattering signal after loss compensation respectively contain Rayleigh scattering noise with different degrees, and can be represented by the following formulas:
Figure BDA0002308382310000101
Figure BDA0002308382310000102
Figure BDA0002308382310000103
Figure BDA0002308382310000104
wherein l is the position distance of a certain point of the optical fiber, h is Planck constant, Δ v ═ 13.2THz is Raman scattering frequency shift in the quartz optical fiber, k is Boltzmann constant, T is the ambient temperature of the optical fiber, αasFor the attenuation coefficient of anti-Stokes scattered light in silica fiber, αsThe attenuation coefficient of Stokes' scattered light in a silica optical fiber, αrAttenuation coefficient of Rayleigh scattered light in quartz fiber, AforScattering coefficient of anti-Stokes scattered light obtained for forward detection, BforDetecting Rayleigh scattering coefficient, C, in channel for anti-Stokes scattered lightforScattering coefficient of anti-Stokes scattered light, D, obtained for forward detectionforRayleigh scattering coefficient, A, in a channel for Stokes scattering light detectionbackScattering coefficient of anti-Stokes scattered light obtained for reverse detection, BbackDetecting Rayleigh scattering coefficient, C, in channel for anti-Stokes scattered lightbackStokes scattering coefficient, D, for reverse detectionbackThe rayleigh scattering coefficient in the channel is detected for the stokes scattered light for reverse detection.
According to formula (1) and formula (2), the time when the first fiber loop (l) is recorded0) And a second optical fiber ring (l)1) Are respectively at a temperature T030 ℃ and T1Intensity of bidirectional raman scattered light at 80 ℃:
Figure BDA0002308382310000105
Figure BDA0002308382310000106
Figure BDA0002308382310000107
Figure BDA0002308382310000108
Figure BDA0002308382310000111
Figure BDA0002308382310000112
Figure BDA0002308382310000113
Figure BDA0002308382310000114
where L2100 m is the fiber length between the first fiber loop and the second fiber loop. A is calculated by the formula (5) and the formula (6)forAnd BforAnd the data are substituted into the formula (7) and the formula (8), and the formula (8) is subtracted from the formula (7) to calculate αr0
C can be calculated from the formula (9), the formula (10), the formula (11) and the formula (12) by the same method as described abovefor、DforAnd αr1
Similarly, A for reverse detection of anti-Stokes scattered light and Stokes scattered light can be calculatedfor、Bfor、Cfor、Dfor、αr2And αr3
Wherein, αrGet αr0、αr1、αr2And αr3Average value of (a).
The signal after forward and backward raman scattering with rayleigh scattering noise removed can be expressed as:
Figure BDA0002308382310000115
Figure BDA0002308382310000116
Figure BDA0002308382310000117
Figure BDA0002308382310000118
geometric mean processing is performed on the formulas (13) and (15), and the formulas (14) and (16), so that:
Figure BDA0002308382310000121
Figure BDA0002308382310000122
in this case, equations (17) and (18) are related only to the temperature T at each position of the optical fiber, and not to the distance.
Dividing equation (17) by equation (18) yields:
Figure BDA0002308382310000123
selecting the intensity of the scattering signal of the first optical fiber ring or the second optical fiber ring and the real-time environment temperature where the scattering signal is located as a reference temperature, and obtaining:
Figure BDA0002308382310000124
dividing equation (19) by equation (20) yields:
Figure BDA0002308382310000125
Figure BDA0002308382310000126
the nonlinear loss caused by factors such as bending and corrosion of strong corrosive substances to the Raman scattering light intensity in the optical fiber is eliminated through the annular light path structure, the long-term stability of the distributed optical fiber temperature sensor in the practical application environment is improved, and the application field of the distributed optical fiber temperature sensor is expanded. And (4) passing. Meanwhile, Rayleigh scattered light in the back Raman scattered light is calculated, and the influence of the Rayleigh scattered light on temperature demodulation precision is eliminated, so that the temperature field information of the sensing optical fiber is directly demodulated, the complexity of system debugging is reduced, and the temperature measurement precision and the long-term measurement stability are improved.
Although the present invention has been described with reference to the preferred embodiments, it is not intended to limit the present invention, and those skilled in the art can make variations and modifications of the present invention without departing from the spirit and scope of the present invention by using the methods and technical contents disclosed above.

Claims (11)

1. A distributed optical fiber temperature sensing system based on an annular light path is characterized by comprising: the device comprises a pulse laser (1), a wavelength division multiplexer (2), an optical switch (3), a first optical fiber ring (5), a second optical fiber ring (6), a temperature measurement optical fiber (7), a photoelectric detector (8), a data acquisition card (9) and an upper computer (10); wherein,
the output ends of the pulse laser (1) and the pulse laser (1) are connected with the input end of the wavelength division multiplexer (2); the common end of the wavelength division multiplexer (2) is connected with the input end of the optical switch (3); two output ends of the optical switch (3) are respectively connected with the first optical fiber ring (5) and the second optical fiber ring (6); the output ends of the first optical fiber ring (5) and the second optical fiber ring (6) are respectively connected with the two ends of the temperature measuring optical fiber (7); two output ends of the wavelength division multiplexer (2) are respectively connected with two input ends of the photoelectric detector (8); two output ends of the photoelectric detector (8) are respectively connected with two input ends of a data acquisition card (9), and the output end of the data acquisition card is connected with the input end of an upper computer (10); the first optical fiber ring (5) and the second optical fiber ring (6) are both placed in the temperature measuring device 4, and the output end of the temperature measuring device 4 is connected with the input end of the upper computer (10);
the distributed optical fiber temperature sensing system of the annular optical path can perform laser power feedback adjustment and Rayleigh scattering noise compensation for inhibiting nonlinear effects.
2. The distributed optical fiber temperature sensing system based on the annular optical path as claimed in claim 1, wherein: the laser power feedback regulation for inhibiting the nonlinear effect comprises the following steps:
the method comprises the following steps: two groups of anti-Stokes scattered light V are acquired by a data acquisition card (9)as1And Vas2And Stokes scattered light Vs1And Vs2Respectively carrying out geometric mean calculation to obtain VasmAnd Vsm
Step two: calculating the ratio of the anti-Stokes scattered light to the Stokes scattered light after geometric averaging, VR=Vasm/VsmAnd calculating the difference DeltaV of the ratio of the positions of the first optical fiber ring (5) and the second optical fiber ring (6)R=VR(L0)-VR(L1) Wherein L is0And L1The positions of the first optical fiber ring (5) and the second optical fiber ring (6) are respectively distant;
step three: starting from 0mA, the pumping current of the pulse laser (1) is gradually increased, and the delta V is observed in real timeRWhen Δ VR>10%×MIN[VR(L0),VR(L1)]The pumping current of the pulse laser (1) is reduced to delta V<1%×MIN[VR(L0),VR(L1)]。
3. The distributed optical fiber temperature sensing system based on the annular optical path as claimed in claim 1, wherein: the rayleigh scattering noise compensation comprises the steps of:
(1) placing the temperature measuring optical fiber (7), the first optical fiber ring (5) and the second optical fiber ring (6) in the same temperature environment, and keeping the environment temperature stable; collecting anti-Stokes scattering signals from positive and negative directions respectivelyObtaining a forward-collected anti-Stokes signal V from a Stokes scattering signalas_forWith Stokes scattering signal Vs_forReversely collected anti-stokes signal Vas_backWith Stokes scattering signal Vs_back
(2) Obtaining optical fiber L according to forward collected anti-Stokes and Stokes scattering signalsnForward anti-stokes scattering signal loss at
Figure RE-FDA0002479182620000021
Loss of signal associated with forward stokes scattering
Figure RE-FDA0002479182620000022
And obtaining compensated forward anti-Stokes scattering signal loss
Figure RE-FDA0002479182620000023
With compensated forward stokes scattering signal loss
Figure RE-FDA0002479182620000024
(3) Obtaining the optical fiber L according to the reversely collected anti-Stokes and Stokes scattering signalsnLoss of reverse anti-stokes scattering signal
Figure RE-FDA0002479182620000025
Loss of signal associated with reverse Stokes scattering
Figure RE-FDA0002479182620000026
And obtaining compensated reverse anti-Stokes scattering signal loss
Figure RE-FDA0002479182620000027
With compensated reverse Stokes scattering signal loss
Figure RE-FDA0002479182620000028
(4) Based on compensated forward anti-Stokes scattering signal loss
Figure RE-FDA0002479182620000029
With compensated forward stokes scattering signal loss
Figure RE-FDA00024791826200000210
Obtaining compensated forward anti-Stokes scattering signal loss containing Rayleigh scattering noise
Figure RE-FDA00024791826200000211
With compensated forward stokes scattering signal loss
Figure RE-FDA00024791826200000212
Based on compensated reverse anti-Stokes scattering signal loss
Figure RE-FDA00024791826200000213
With compensated reverse Stokes scattering signal loss
Figure RE-FDA00024791826200000214
Obtaining compensated reverse anti-stokes scattering signal loss containing Rayleigh scattering noise
Figure RE-FDA00024791826200000215
With compensated reverse Stokes scattering signal loss
Figure RE-FDA00024791826200000216
4. The distributed optical fiber temperature sensing system based on the annular optical path as claimed in claim 3, wherein: in step (2), the optical fiber LnForward acquired anti-stokes and stokes scattered signalsTo obtain an optical fiber LnForward anti-stokes scattering signal loss at
Figure RE-FDA00024791826200000217
Loss of signal associated with forward stokes scattering
Figure RE-FDA00024791826200000218
Respectively expressed as:
Figure RE-FDA0002479182620000031
Figure RE-FDA0002479182620000032
wherein m is the number of translation steps of the loss point of the optical fiber.
5. The distributed optical fiber temperature sensing system based on the annular optical path as claimed in claim 3, wherein: in step (2), the compensated forward anti-Stokes scattering signal loss
Figure RE-FDA0002479182620000033
With compensated forward stokes scattering signal loss
Figure RE-FDA0002479182620000034
Respectively expressed as:
Figure RE-FDA0002479182620000035
Figure RE-FDA0002479182620000036
6. the distributed optical fiber temperature sensing system based on the annular optical path as claimed in claim 5, wherein: in thatIn step (3), the optical fiber LnLoss of reverse anti-stokes scattering signal
Figure RE-FDA0002479182620000037
Loss of signal associated with reverse Stokes scattering
Figure RE-FDA0002479182620000038
Respectively expressed as:
Figure RE-FDA0002479182620000039
Figure RE-FDA00024791826200000310
7. the distributed optical fiber temperature sensing system based on the annular optical path as claimed in claim 6, wherein: in step (3), the compensated reverse anti-Stokes scattering signal loss
Figure RE-FDA00024791826200000311
With compensated reverse Stokes scattering signal loss
Figure RE-FDA00024791826200000312
Respectively expressed as:
Figure RE-FDA00024791826200000313
Figure RE-FDA00024791826200000314
8. the distributed optical fiber temperature sensing system based on the annular optical path as claimed in claim 7, wherein: in step (4), the compensated forward inverse gaussian noise including the rayleigh scattering noise is includedLoss of Thauss scattering signal
Figure RE-FDA0002479182620000041
Comprises the following steps:
Figure RE-FDA0002479182620000042
wherein l is the position distance of a certain point of the optical fiber, h is Planck constant, Deltav is Raman scattering frequency shift in the quartz optical fiber, k is Boltzmann constant, T is the environment temperature of the optical fiber, αasFor the attenuation coefficient of anti-Stokes scattered light in silica fiber, αsThe attenuation coefficient of Stokes' scattered light in a silica optical fiber, αrAttenuation coefficient of Rayleigh scattered light in quartz fiber, AforScattering coefficient of anti-Stokes scattered light obtained for forward detection, BforDetecting Rayleigh scattering coefficient, C, in channel for anti-Stokes scattered lightforScattering coefficient of anti-Stokes scattered light, D, obtained for forward detectionforRayleigh scattering coefficient, A, in a channel for Stokes scattering light detectionbackScattering coefficient of anti-Stokes scattered light obtained for reverse detection, BbackDetecting Rayleigh scattering coefficient, C, in channel for anti-Stokes scattered lightbackStokes scattering coefficient, D, for reverse detectionbackRayleigh scattering coefficient in the channel is detected for Stokes scattering light for reverse detection, αr0The attenuation coefficient of Rayleigh scattered light in a forward anti-Stokes detection signal in a quartz optical fiber is detected.
9. The distributed optical fiber temperature sensing system based on the annular optical path as claimed in claim 7, wherein: in step (4), the compensated forward stokes scattering signal loss containing rayleigh scattering noise
Figure RE-FDA0002479182620000043
Comprises the following steps:
Figure RE-FDA0002479182620000044
wherein l is the position distance of a certain point of the optical fiber, h is Planck constant, Deltav is Raman scattering frequency shift in the quartz optical fiber, k is Boltzmann constant, T is the environment temperature of the optical fiber, αasFor the attenuation coefficient of anti-Stokes scattered light in silica fiber, αsThe attenuation coefficient of Stokes' scattered light in a silica optical fiber, αrAttenuation coefficient of Rayleigh scattered light in quartz fiber, CforScattering coefficient of anti-Stokes scattered light, D, obtained for forward detectionforRayleigh scattering coefficient in the detection channel for Stokes scattered light, αr1The attenuation coefficient of Rayleigh scattered light in a forward Stokes detection signal in a quartz optical fiber is detected.
10. The distributed optical fiber temperature sensing system based on the annular optical path as claimed in claim 7, wherein: in step (4), the compensated anti-Stokes scattering signal loss containing Rayleigh scattering noise
Figure RE-FDA0002479182620000051
Comprises the following steps:
Figure RE-FDA0002479182620000052
wherein l is the position distance of a certain point of the optical fiber, h is Planck constant, Deltav is Raman scattering frequency shift in the quartz optical fiber, k is Boltzmann constant, T is the environment temperature of the optical fiber, αasFor the attenuation coefficient of anti-Stokes scattered light in silica fiber, αsThe attenuation coefficient of Stokes' scattered light in a silica optical fiber, αrAttenuation coefficient of Rayleigh scattered light in quartz fiber, AbackScattering coefficient of anti-Stokes scattered light obtained for reverse detection, BbackDetection of Rayleigh scattered light scattering coefficients in a channel for anti-Stokes scattered light, αr2The attenuation coefficient of Rayleigh scattered light in a signal in a quartz fiber is detected for reverse anti-Stokes.
11. The distributed optical fiber temperature sensing system based on the annular optical path as claimed in claim 7, wherein: in step (4), the compensated reverse Stokes scattering signal loss containing Rayleigh scattering noise
Figure RE-FDA0002479182620000053
Comprises the following steps:
Figure RE-FDA0002479182620000054
wherein l is the position distance of a certain point of the optical fiber, h is Planck constant, Deltav is Raman scattering frequency shift in the quartz optical fiber, k is Boltzmann constant, T is the environment temperature of the optical fiber, αasFor the attenuation coefficient of anti-Stokes scattered light in silica fiber, αsThe attenuation coefficient of Stokes' scattered light in a silica optical fiber, αrAttenuation coefficient of Rayleigh scattered light in quartz fiber, CbackStokes scattering coefficient, D, for reverse detectionbackRayleigh scattering coefficient in the channel is detected for Stokes scattering light for reverse detection, αr3The attenuation coefficient of Rayleigh scattered light in a signal in a quartz fiber is detected for the reverse Stokes.
CN201911248571.6A 2019-12-09 2019-12-09 Distributed optical fiber temperature sensing system based on annular light path Active CN111307322B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911248571.6A CN111307322B (en) 2019-12-09 2019-12-09 Distributed optical fiber temperature sensing system based on annular light path

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911248571.6A CN111307322B (en) 2019-12-09 2019-12-09 Distributed optical fiber temperature sensing system based on annular light path

Publications (2)

Publication Number Publication Date
CN111307322A true CN111307322A (en) 2020-06-19
CN111307322B CN111307322B (en) 2021-12-07

Family

ID=71156207

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911248571.6A Active CN111307322B (en) 2019-12-09 2019-12-09 Distributed optical fiber temperature sensing system based on annular light path

Country Status (1)

Country Link
CN (1) CN111307322B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20220166579A (en) * 2021-06-10 2022-12-19 한국광기술원 apparatus for measuring distribution temperature by using optical fiber and method of correction temperature thereof
CN115931168A (en) * 2023-03-13 2023-04-07 之江实验室 Attenuation difference compensation method for distributed optical fiber temperature sensing system

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB9307660D0 (en) * 1993-04-14 1993-06-02 York Ltd Distributed sensing
US20130100984A1 (en) * 2011-10-24 2013-04-25 Yokogawa Electric Corporation Opticalfiber temperature distribution measurement apparatus
CN103575424A (en) * 2013-11-27 2014-02-12 威海北洋电气集团股份有限公司 Temperature fluctuation controllable distributed fiber sensor as well as temperature measuring system and method
CN106768278A (en) * 2017-01-06 2017-05-31 天津大学 A kind of distributed optical fiber vibration and the double physical quantity sensing positioning systems of temperature
CN107091698A (en) * 2017-06-16 2017-08-25 苏州光格设备有限公司 Brillouin optical time domain analysis system and method
CN110307920A (en) * 2019-06-12 2019-10-08 太原理工大学 Based on noise-modulated fiber optic temperature, stress sensing system and measurement method

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB9307660D0 (en) * 1993-04-14 1993-06-02 York Ltd Distributed sensing
US20130100984A1 (en) * 2011-10-24 2013-04-25 Yokogawa Electric Corporation Opticalfiber temperature distribution measurement apparatus
CN103575424A (en) * 2013-11-27 2014-02-12 威海北洋电气集团股份有限公司 Temperature fluctuation controllable distributed fiber sensor as well as temperature measuring system and method
CN106768278A (en) * 2017-01-06 2017-05-31 天津大学 A kind of distributed optical fiber vibration and the double physical quantity sensing positioning systems of temperature
CN107091698A (en) * 2017-06-16 2017-08-25 苏州光格设备有限公司 Brillouin optical time domain analysis system and method
CN110307920A (en) * 2019-06-12 2019-10-08 太原理工大学 Based on noise-modulated fiber optic temperature, stress sensing system and measurement method

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20220166579A (en) * 2021-06-10 2022-12-19 한국광기술원 apparatus for measuring distribution temperature by using optical fiber and method of correction temperature thereof
KR102590280B1 (en) * 2021-06-10 2023-10-17 한국광기술원 apparatus for measuring distribution temperature by using optical fiber and method of correction temperature thereof
CN115931168A (en) * 2023-03-13 2023-04-07 之江实验室 Attenuation difference compensation method for distributed optical fiber temperature sensing system
CN115931168B (en) * 2023-03-13 2023-06-16 之江实验室 Attenuation difference compensation method of distributed optical fiber temperature sensing system

Also Published As

Publication number Publication date
CN111307322B (en) 2021-12-07

Similar Documents

Publication Publication Date Title
WO2020034558A1 (en) High-precision temperature demodulation method oriented toward distributed fiber raman sensor
CN108458814B (en) Self calibration detection device and temperature demodulation method towards fiber Raman temperature-sensing system
US10018517B2 (en) Optical fiber temperature distribution measuring device
CN104596670B (en) A kind of method for solving distributed fiber Raman temperature-sensing system temperature drift
CN104864979B (en) A kind of modification method of distributed Raman fiber temp measuring system measurement error
CN111307322B (en) Distributed optical fiber temperature sensing system based on annular light path
CN111006788B (en) High-precision optical fiber Raman temperature detection method based on anti-Stokes light self-demodulation
CN111006786B (en) Double-path high-precision temperature demodulation method based on distributed optical fiber Raman sensing system
US20140254629A1 (en) Optical fiber temperature distribution measurement device and method of measuring optical fiber temperature distribution
CN104101447A (en) Distributed optical fiber temperature sensor and method for removing nonlinear error of same
CN109029770B (en) Distributed optical fiber Raman temperature and strain demodulation method based on loop demodulation
CN110307920B (en) Optical fiber temperature and stress sensing system based on noise modulation and measuring method
CN106404215A (en) Design of distributed fiber sensing system based on Brillouin scattering
CN102279062A (en) Temperature calibration and demodulation method for distributed optical fiber sensing temperature measurement system
CN106225816B (en) A kind of grating sensing apparatus and method based on Brillouin&#39;s wave filter
CN111006787B (en) Distributed optical fiber Raman double-end temperature demodulation method based on differential temperature compensation
CN114942087B (en) Raman Stokes light intensity dynamic self-calibration distributed temperature detection system and method
CN111795759A (en) Method and system for dispersion correction of distributed Raman fiber temperature measurement system
Li et al. R-DTS with heat transfer functional model for perceiving the surrounding temperature
CN113654683A (en) Calibration method and device for distributed optical fiber temperature measurement system
CN111141414B (en) Temperature and strain simultaneous measurement device and method based on chaos BOCDA
CN112880865B (en) Ultra-long-distance high-spatial-resolution Raman optical fiber dual-parameter sensing system and method
CN108120461B (en) Fiber Bragg grating demodulation method and device based on chaotic light source
CN105092087A (en) Photoelectric conversion module, temperature compensation method for photoelectric conversion module, and distributed light sensing system
CN110057386A (en) The demodulation method and its demodulating equipment of optical fiber FP sensor based on full spectrum

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