Disclosure of Invention
Aiming at the problems, the invention provides a method for demodulating the temperature along the optical fiber in a Raman optical time domain reflectometer, which is characterized in that a standard reference temperature and a reference optical fiber are set in the Raman optical time domain reflectometer, a reference temperature conversion coefficient is obtained by utilizing a reference section temperature-Raman intensity ratio curve and the position temperature, and the temperature distribution along the optical fiber is demodulated by combining a temperature calibration curve, so that the problem of deviation of a demodulation temperature result under different reference temperatures caused by light attenuation, unstable devices and the like is solved, and the accuracy and the stability of an optical time domain reflectometer system are improved.
In order to achieve the purpose of the invention, the invention provides a method for demodulating the temperature along the optical fiber in the Raman optical time domain reflectometer, which comprises the following steps:
s10, setting a certain section of the optical fibers as a reference optical fiber, installing a thermometer on the section, and setting the rest optical fibers as sections to be tested;
s20, placing the reference optical fiber in a thermostat with a specific reference temperature, placing the optical fiber to be measured in another thermostat, and adjusting the temperature of the thermostat of the section to be measured to obtain the anti-Stokes and Stokes light intensity ratios of the two sections at different temperatures to be measured;
s30, adjusting the temperature of the thermostat at the reference fiber to obtain the anti-Stokes and Stokes light intensity ratio of the reference fiber at different reference temperatures;
s50, fitting according to the anti-Stokes and Stokes light intensity ratios of the reference optical fibers at different reference temperatures to obtain anti-Stokes and Stokes light intensity ratio curves of the reference optical fibers at different reference temperatures;
s70, determining a proportionality coefficient of the conversion reference temperature according to the thermometer result, the optical signal intensity ratio and the anti-Stokes-to-Stokes light intensity ratio at the reference optical fiber;
and S80, substituting the intensity ratio and the proportionality coefficient of the section to be measured into a specific demodulation formula by using the specific reference temperature and the Raman intensity ratio thereof, and calculating the temperature of the optical fiber to be measured by combining a calibration formula.
In one embodiment, after step S30, the method further includes:
and S40, obtaining a calibration formula at the specific reference temperature by using the calibration formula and the known temperature to be measured.
Specifically, the original formula includes:
the calibration formula comprises:
T=1.0917T1-10.792,
wherein, T0Indicating a specific reference temperature, T indicating the true temperature of the section to be measured, T1The temperature of the section to be measured which is calculated by using a demodulation formula below and is not calibrated in the last step is represented, h represents a Planck constant, c represents an optical speed, k represents a Boltzmann constant, v represents a Raman frequency shift, R (T) represents the anti-Stokes-to-Stokes light intensity ratio of the optical fiber of the section to be measured, and R (T)0) Representing the anti-stokes to stokes light intensity ratio of the reference segment fiber.
In one embodiment, the specific demodulation formula includes:
wherein, T0Denotes a specific reference temperature, h denotes a Planckian constant, c denotes an optical speed, k denotes a Boltzmann constant, v denotes a Raman frequency shift, R (T)1) Representing the anti-Stokes to Stokes light intensity ratio, T, of the section of fibre to be measured1Representing the temperature of the section to be measured calculated by using a demodulation formula, which needs to be further calibrated to obtain the final real temperature of the section to be measured, x representing the proportionality coefficient of the conversion reference temperature, R (T)0) Representing the anti-stokes to stokes light intensity ratio of the reference segment fiber.
In one embodiment, after step S50, the method further includes:
s60, performing field measurement to obtain the anti-stokes optical signal intensity ratio and the temperature at the reference fiber.
The demodulation method for the temperature along the optical fiber in the Raman optical time domain reflectometer has the following technical effects:
1. and the accuracy and stability of the measurement result are improved by using the data calculation of the reference optical fiber.
2. The reference optical fiber does not need to be constant in temperature, so that the system cost is reduced, and the application range of the system is widened.
Detailed Description
In order to make the objects, technical solutions and advantages of the present application more apparent, the present application is described in further detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the present application and are not intended to limit the present application.
Reference herein to "an embodiment" means that a particular feature, structure, or characteristic described in connection with the embodiment can be included in at least one embodiment of the application. The appearances of the phrase in various places in the specification are not necessarily all referring to the same embodiment, nor are separate or alternative embodiments mutually exclusive of other embodiments. It is explicitly and implicitly understood by one skilled in the art that the embodiments described herein can be combined with other embodiments.
Referring to fig. 1, fig. 1 is a flowchart of a demodulation method of temperature along a fiber in a raman optical time domain reflectometer according to an embodiment, including the following steps:
s10, a certain section of the optical fibers is set as a reference optical fiber, a thermometer is installed in the section, and the remaining optical fibers are set as a section to be measured.
S20, placing the reference optical fiber in a thermostat with a specific reference temperature, placing the optical fiber to be measured in another thermostat, and adjusting the temperature of the thermostat of the section to be measured to obtain the anti-Stokes and Stokes light intensity ratio of the two sections at different temperatures to be measured.
The steps can realize calibration, and a calibration curve is made according to the anti-Stokes and Stokes light intensity ratios of the two sections at different temperatures to be measured for correcting results in subsequent measurement.
And S30, adjusting the temperature of the thermostat at the reference fiber to obtain the anti-Stokes and Stokes light intensity ratio of the reference fiber at different reference temperatures.
And S50, fitting according to the anti-Stokes and Stokes light intensity ratios of the reference optical fibers at different reference temperatures to obtain anti-Stokes and Stokes light intensity ratio curves of the reference optical fibers at different reference temperatures.
And S70, determining a proportionality coefficient of the conversion reference temperature according to the thermometer result, the optical signal intensity ratio and the anti-Stokes-to-Stokes light intensity ratio at the reference optical fiber.
Specifically, the anti-stokes and stokes light intensity ratios of the reference fiber at different reference temperatures of the reference segment (reference fiber) are obtained in step S30, and a temperature-light intensity ratio curve is obtained after fitting. The light intensity ratio of the corresponding temperature in step S30 is divided by the actually measured light intensity ratio, so that the proportionality coefficient of the conversion reference temperature can be obtained.
And S80, substituting the intensity ratio and the proportionality coefficient of the section to be measured into a specific demodulation formula by using the specific reference temperature and the Raman intensity ratio thereof, and calculating the temperature of the optical fiber to be measured by combining a calibration formula.
This step sets the specific reference temperature T0And the anti-Stokes to Stokes light intensity ratio R (T) of the reference section fiber0) And the calculated proportionality coefficient x and the Stokes light and anti-Stokes light intensity ratio R (T) of the segment to be measured1) Substituting into a specific demodulation formula to obtain the calculated temperature T1And finally, obtaining the measured temperature T by using a calibration formula.
In this embodiment, the raman optical time domain reflectometer includes a laser, a photodetector, a wavelength division multiplexer, an optical fiber, and a collection card, where the wavelength division multiplexer is configured to separate stokes light and anti-stokes light and transmit the separated stokes light and anti-stokes light to the photodetector in two paths, and the separated stokes light and anti-stokes light are collected by the collection card.
The demodulation method of the temperature along the optical fiber in the Raman optical time domain reflectometer comprises the steps of setting a certain section of the optical fiber as a reference optical fiber, installing a thermometer on the section, setting the rest optical fibers as sections to be detected, placing the reference optical fiber in a constant temperature box with a specific reference temperature, placing the optical fiber to be detected in another constant temperature box, adjusting the temperature of the constant temperature box at the section to be detected to obtain the anti-Stokes and Stokes light intensity ratios of the two sections at different temperatures to be detected, adjusting the temperature of the constant temperature box at the reference optical fiber to obtain the anti-Stokes and Stokes light intensity ratios of the reference optical fiber at different reference temperatures, fitting the anti-Stokes and Stokes light intensity ratios of the reference optical fiber at different reference temperatures to obtain anti-Stokes and Stokes light intensity ratio curves of the reference optical fiber at different reference temperatures, and fitting the anti-Stokes and Stokes light intensity ratios according to the thermometer result, the optical signal intensity ratio and the anti-Stokes and Stokes light intensity ratios at the reference optical fiber, determining a proportional coefficient of the conversion reference temperature, substituting the specific reference temperature and the Raman intensity ratio thereof into a specific demodulation formula, and calculating the temperature of the optical fiber to be measured by combining a calibration formula, wherein the obtained temperature of the optical fiber to be measured has higher accuracy; the data calculation of the reference optical fiber is utilized, the accuracy and the stability of the measuring result are improved, the reference optical fiber does not need to be constant in temperature, the system cost is reduced, and the application range of the system is widened.
In one embodiment, after step S30, the method further includes:
and S40, obtaining a calibration formula at the specific reference temperature by using the calibration formula and the known temperature to be measured.
Specifically, the calculation process mentioned in step S70 is to first obtain the proportionality coefficient x of the conversion reference temperature by using the reference segment real temperature and the anti-stokes and stokes light intensity ratio curves at different reference temperatures calculated in step S40.
Specifically, the original formula includes:
the calibration formula comprises:
T=1.0917T1-10.792,
wherein, T0Indicating a specific reference temperature, T indicating the true temperature of the section to be measured, T1The temperature of the section to be measured which is calculated by using a demodulation formula and is not calibrated is represented, h represents a Planck constant, c represents an optical speed, k represents a Boltzmann constant, v represents a Raman frequency shift, R (T) represents the anti-Stokes-to-Stokes light intensity ratio of the section to be measured, and R (T)0) Representing the anti-stokes to stokes light intensity ratio of the reference segment fiber.
In one embodiment, the specific demodulation formula includes:
wherein, T0Denotes a specific reference temperature, h denotes a Planckian constant, c denotes an optical speed, k denotes a Boltzmann constant, v denotes a Raman frequency shift, R (T)1) Representing the anti-Stokes to Stokes light intensity ratio, T, of the segment to be measured1Expressing the temperature of the section to be measured obtained by using a demodulation formula, x expressing the proportionality coefficient of the conversion reference temperature, R (T)0) Representing the anti-stokes to stokes light intensity ratio of the reference segment fiber.
In one embodiment, after step S50, the method further includes:
s60, performing field measurement to obtain the anti-stokes optical signal intensity ratio and the temperature at the reference fiber.
In an embodiment, the ROTDR structure diagram can be shown in fig. 2, and the ROTDR system is a distributed optical fiber temperature sensing system based on raman effect and optical time domain reflectometry, and includes optoelectronic devices such as a laser, a wavelength division multiplexer, an optical-to-electrical converter, an acquisition card, and an optical fiber. The working principle is that a laser emits pulse laser, a Raman scattering effect occurs when the pulse laser is transmitted in an optical fiber, Stokes light and anti-Stokes light are generated, the Stokes light and the anti-Stokes light are separated by a wavelength division multiplexer, light intensity data are collected through a photoelectric converter and a collection card, the Raman intensity ratio of the Stokes light divided by the anti-Stokes light is calculated, and the temperature along the optical fiber can be obtained through demodulation by combining a temperature calibration curve under a reference temperature T0. The demodulation formula is:
wherein, T0Representing reference temperature, T representing temperature to be measured, h, c, k and v representing Planck constant, light speed, Boltzmann constant and Raman frequency shift, R (T) representing anti-Stokes and Stokes light intensity ratio of optical fiber to be measured, and R (T)0) And represents the anti-stokes and stokes light intensity ratio of the reference section fiber to the reference section fiber.
In practical application, due to attenuation and loss of light propagating in the optical fiber, parameter offset of devices such as a laser and an acquisition card during long-time operation, and the like, the temperature result obtained by directly demodulating the optical signal by using the formula is not accurate. Some processing of the temperature calculation results is usually required using calibration equations, however, the calibration equations are not the same at different reference temperatures. Therefore, once the reference fiber of the calibration formula is determined, the reference fiber must be placed in an incubator, otherwise, the demodulated result has deviation, which limits the use of the raman optical time domain reflectometer and reduces the accuracy of measurement. The demodulation method for the temperature along the optical fiber in the raman optical time domain reflectometer provided in this embodiment is to solve this problem, and reference may be made to fig. 3 for a corresponding working flow.
After the ROTDR system is constructed, the reference fiber and the fiber to be measured are respectively placed in two thermostats, if the temperature of the reference fiber is 16 ℃ and the temperature of the fiber to be measured is 60 ℃, after the collected optical signal is measured, T0 is 16 ℃, and R (T) and R (T0) are substituted into the original formula, and the result is shown in fig. 4(a), wherein about the first 300 meters in the figure are the reference section, and 300 and 500 meters are the sections to be measured, and the following figures are the same. It can be seen that the accuracy of the measurement result is poorer than the actual temperature of 60 ℃. And changing the temperature of the optical fiber to be measured, for example, measuring the temperature once every 5 ℃ or 10 ℃ within the range of 20 ℃ to 900 ℃, obtaining the original formula calculation value and the real temperature value under different temperatures, and obtaining a linear calibration formula, as shown in fig. 5. Wherein T is the true temperature, T1Temperature calculated for direct demodulation: 1.0917T1-10.792。
The temperature of the thermostat of the section to be measured is changed again, for example, 70 ℃, the calculation is carried out by using the original formula and the calibration formula, and the obtained result is shown in fig. 4(b), so that the calibration formula is effective at the reference temperature of 16 ℃.
If the temperature of the reference optical fiber is changed to 25 ℃ and the temperature of the optical fiber to be measured is set to 60 ℃, the temperature of the two thermostats is measured again, the original calibration formula is still used for demodulation to obtain a result as shown in fig. 6, the accuracy of the result is reduced, and the experimental formula obtained under the condition of 16 ℃ is proved to be not suitable for the condition that the reference temperature is 25 ℃.
And changing the temperature of the constant temperature box of the reference section to obtain the anti-Stokes and Stokes light intensity ratio of the reference section at different reference temperatures, for example, measuring every 5 ℃ within the range of 15 ℃ to 40 ℃, and fitting the anti-Stokes and Stokes light intensity ratio curves of the points, as shown in figure 7.
Performing actual measurement, for example, measuring the known reference temperature of 25 deg.C and its anti-Stokes-to-Stokes light intensity ratio, combining the reference segment anti-Stokes-to-Stokes light intensity ratio curve, and calculating the light intensity ratio of the corresponding temperature in the curve at 25 deg.C to the actual temperatureThe ratio of the measured light intensity ratio, i.e. the value of the corresponding point in the curve is divided by the value obtained by measurement, so as to obtain the proportionality coefficient x of the conversion reference temperature which is 1.032:
and then the obtained x and the collected anti-Stokes light intensity ratio R (T) of the segment to be measured to the Stokes light intensity ratio R1) A specific reference temperature T016 ℃ and R (T)0) Calculating the temperature T1The result of the section to be measured can be obtained by combining a calibration formula as shown in FIG. 6, and the accuracy can be seen to be higher.
The demodulation method for the temperature along the optical fiber in the Raman optical time domain reflectometer provided by the embodiment has the following technical effects:
1. and the accuracy and stability of the measurement result are improved by using the data calculation of the reference optical fiber.
2. The reference optical fiber does not need to be constant in temperature, so that the system cost is reduced, and the application range of the system is widened.
The technical features of the above embodiments can be arbitrarily combined, and for the sake of brevity, all possible combinations of the technical features in the above embodiments are not described, but should be considered as the scope of the present specification as long as there is no contradiction between the combinations of the technical features.
It should be noted that the terms "first \ second \ third" referred to in the embodiments of the present application merely distinguish similar objects, and do not represent a specific ordering for the objects, and it should be understood that "first \ second \ third" may exchange a specific order or sequence when allowed. It should be understood that "first \ second \ third" distinct objects may be interchanged under appropriate circumstances such that the embodiments of the application described herein may be implemented in an order other than those illustrated or described herein.
The terms "comprising" and "having" and any variations thereof in the embodiments of the present application are intended to cover non-exclusive inclusions. For example, a process, method, apparatus, product, or device that comprises a list of steps or modules is not limited to the listed steps or modules but may alternatively include other steps or modules not listed or inherent to such process, method, product, or device.
The above-mentioned embodiments only express several embodiments of the present application, and the description thereof is more specific and detailed, but not construed as limiting the scope of the invention. It should be noted that, for a person skilled in the art, several variations and modifications can be made without departing from the concept of the present application, which falls within the scope of protection of the present application. Therefore, the protection scope of the present patent shall be subject to the appended claims.