Temperature measurement device and temperature measurement method

Abstract

To perform temperature measurement with high accuracy over a long period of time.SOLUTION: A temperature measurement device 1 includes a control unit 31 that acquires a Brillouin frequency shift from a measurement position of an optical fiber 10 to be measured, and, when an exposure time (t), which is a period in which the measurement position of the optical fiber 10 is exposed to a prescribed temperature environment, has reached a prescribed time, corrects the Brillouin frequency shift in accordance with a drift to a high frequency side based on the exposure time (t), and calculates a measurement temperature of the measurement position on the basis of the corrected Brillouin frequency shift.SELECTED DRAWING: Figure 2

Description

[Technical field]

[0001] The present invention relates to a temperature measuring device and a temperature measuring method. [Background technology]

[0002] In recent years, there has been an increasing demand for temperature measurement in high-temperature environments in industries such as the plant industry, heavy electrical equipment industry, and the upstream oil and natural gas industry. Optical fibers made of silica glass have a high melting point (over 1,000°C) and can be used in high-temperature environments, and are therefore attracting attention as high-temperature sensors.

[0003] As a distributed temperature sensor (DTS) that uses an optical fiber as a sensor for measuring temperature, Raman Optical Time Domain Reflectometry (ROTDR) that uses Raman scattered light is generally known and has already been put to practical use.

[0004] However, due to problems such as changes (increases) in the transmission loss of the optical fiber that acts as the sensor in high-temperature environments, the temperature measurement range of ROTDR is currently limited to 300°C or less (Non-Patent Document 1). On the other hand, DTS using Brillouin scattering calculates temperature from the frequency difference between the incident light and the Brillouin scattered light (Brillouin Frequency Shift: BFS), so it is less susceptible to changes in the transmission loss of the optical fiber and is expected to realize high-temperature measurements of 300°C or more.

[0005] In recent years, many useful experimental results have been reported from many research institutes, and the relationship between BFS and temperature has been considered linear, but recent research has reported that it becomes nonlinear above 500°C (Non-Patent Document 2). The reason for this is thought to be that the Young's modulus of optical fiber is nonlinear with respect to temperature.

[0006] Brillouin scattering is scattering caused by acoustic waves in optical fibers. There are two types of Brillouin scattering: spontaneous Brillouin scattering and stimulated Brillouin scattering. Spontaneous Brillouin scattering is scattering caused by acoustic waves that exist naturally in optical fibers. Stimulated Brillouin scattering is a scattering phenomenon caused by the interaction of two lights: light that enters one end of the optical fiber (called pump light) and light that enters in the opposite direction from the other end of the optical fiber and has a frequency about 11 GHz lower than that of the pump light (called probe light).

[0007] Representative DTS using spontaneous Brillouin scattering include BOTDR (Brillouin Optical Time Domain Reflectometry) and BOCDR (Brillouin Optical Correlation Domain Reflectometry), while representative DTS using stimulated Brillouin scattering include BOTDA (Brillouin Optical Time Domain Analysis) and BOCDA (Brillouin Optical Correlation Domain Analysis).

[0008] Fig. 5 is an explanatory diagram for explaining an outline of a conventional temperature measuring device. As shown in Fig. 5, temperature (high temperature) measurement in DTS is performed using an optical fiber characteristic measuring device 100 and an optical fiber sensor 101 (measured fiber). Here, in the optical fiber sensor 101, position 102 is a position where the temperature is relatively high, and position 103 is a position where the temperature is relatively low.

[0009] Fig. 6 is an explanatory diagram explaining the relationship between temperature and Brillouin frequency shift. As shown in Fig. 6, the optical fiber characteristic measuring device 100 converts the BFS measured during temperature measurement into temperature using the following BFS-temperature conversion formula (1) calculated from the relationship between the BFS and temperature (T) of the optical fiber sensor 101 measured in advance.

[0010]

number

[0011] Here, A and B in the BFS-temperature conversion formula (1) are coefficients, and C is a constant. The coefficients A and B and the constant C use values ​​calculated from the relationship between the BFS and temperature (T) in the optical fiber sensor 101 that has been measured in advance.

[0012] Fig. 7 is a diagram showing an example of the Brillouin frequency shift with respect to the distance in the longitudinal direction of the optical fiber sensor 101. When a position 102 where the temperature is relatively high and a position 103 where the temperature is relatively low exist in the longitudinal direction of the optical fiber sensor 101, a BFS distribution according to the distance is measured as shown in Fig. 7.

[0013] Fig. 8 is a diagram showing an example of temperature versus distance in the longitudinal direction of the optical fiber sensor. The optical fiber characteristic measuring device 100 can measure the temperature distribution versus distance in the longitudinal direction as shown in Fig. 8 by applying the BFS-temperature conversion formula (1) to the BFS distribution versus distance in the longitudinal direction of the optical fiber sensor 101 measured in this way.

[0014] Furthermore, in the optical fiber characteristic measuring apparatus 100, by using BOCDR and BOCDA that can measure the BFS for each position of the optical fiber sensor 101 at high speed, it is also possible to measure the change in temperature over time at any position.

[0015] Fig. 9 is an explanatory diagram explaining the change over time of the Brillouin frequency shift for each position of the optical fiber sensor. As shown in Fig. 9, the optical fiber characteristic measuring device 100 can measure the change over time of the Brillouin frequency shift at positions a1, a2, and a3 from the optical fiber characteristic measuring device 100 side in the longitudinal direction of the optical fiber sensor 101 by using BOCDR and BOCDA.

[0016] Fig. 10 is an explanatory diagram illustrating the change in temperature over time at each position of the optical fiber sensor. As shown in Fig. 10, the optical fiber characteristic measuring device 100 can measure the temperature change at the positions a1 to a3 by applying the BFS-temperature conversion formula (1) to the change in Brillouin frequency shift over time at the measured positions a1 to a3.

[0017] However, since the BFS is sensitive to strain and temperature, it is necessary to devise a method to prevent strain from being applied to the optical fiber sensor 101 during temperature measurement. As one example, a commonly used method is to insert the optical fiber sensor 101 into a metal tube and then fix the metal tube to an external object to be measured. The optical fiber sensor 101 is not affected by strain because it is free inside the metal tube. There is also a temperature measurement method that eliminates the effects of strain by using a conventional technology (Non-Patent Document 3) that separates strain and temperature. In the following explanation, it is assumed that no strain is applied to the optical fiber sensor 101. [Prior art documents] [Patent documents]

[0018] [Patent Document 1] Patent No. 5012032 [Patent Document 2] Patent No. 5831189 [Non-patent literature]

[0019] [Non-Patent Document 1] I. Laarossi, MA Quintela, JM Lopez-Higuera, "Comparative Experimental Study of a High-Temperature Raman-Based Distributed Optical Fiber Sensor with Different Special Fibers," Sensors 19(3), 574-587 (2019) [Non-Patent Document 2] Y. Li, F. Zhang, and T. Yoshino, "Wide-range temperature dependence of Brillouin shift in a dispersion-shifted fiber and its annealing effect," J. Lightwave Technol. 21(7), 1663-1667 (2003) [Non-Patent Document 3] Weiwen Zou, Zuyuan He, and Kazuo Hotate, “Complete discrimination of strain and temperature using Brillouin frequency shift and birefringence in a polarization-maintaining fiber,” Optics Express, Vol. 17, No. 3, 2009, pp. 1248-1255 [Non-Patent Document 4] Kazuya Saito, Akira Ikushima, "Structural Relaxation and Optical Properties of Silica Glass", NEW GLASS Vol.14 No.4 1999 [Non-Patent Document 5] Kazuya Saito, "What is fictive temperature?" NEW GLASS Vol.27 No.104 2012 [Non-Patent Document 6] Akio Koike, "Spectroscopic evaluation method for fictive temperature of glass", NEW GLASS Vol.27 No.105 2012 [Non-Patent Document 7] Hirokazu MASAI, et al., “Relationship between the first sharp diffraction peak and physical properties of silicon dioxide (SiO2) glasses possessing different fictive temperatures,” Journal of the Ceramic Society of Japan 128

[12] 1038-1044 2020” Summary of the Invention [Problem to be solved by the invention]

[0020] In recent years, there has been an increasing demand for temperature measurement in high-temperature environments, and temperature measurement for long periods of time (tens of thousands of hours) is expected in high-temperature environments of 300°C or higher. Therefore, the inventors conducted an experiment assuming such long-term temperature measurement (tens of thousands of hours).

[0021] Fig. 11 shows the experimental results of the time change of the Brillouin frequency shift when the optical fiber sensor is exposed to a constant temperature environment. Experimental result R1 in Fig. 11 shows the time change of the Brillouin frequency shift when the optical fiber sensor is exposed to a constant temperature environment of 800°C. Experimental result R2 shows the time change of the Brillouin frequency shift when the optical fiber sensor is exposed to a constant temperature environment of 850°C.

[0022] From the experimental results R1 and R2 shown in Fig. 11, the inventors learned that the BFS of the optical fiber sensor measured in a constant temperature high-temperature environment changes (drifts) with exposure time. Here, the drift of the BFS to the low-frequency side with exposure time is called low-frequency side drift, and the drift to the high-frequency side is called high-frequency side drift.

[0023] It was newly discovered that the BFS, which changes with exposure time, not only exhibits a low-frequency drift, but also exhibits a high-frequency drift at a certain exposure time. For example, in the experimental result R1 at 800°C, a high-frequency drift of 26 MHz occurs after about 12,000 hours, and in the experimental result R2 at 850°C, a high-frequency drift of 50 MHz occurs after about 12,000 hours. For this reason, it was found that it is difficult to perform high-precision temperature measurements over long periods of time using conventional temperature measurements based on the BFS.

[0024] The present invention has been made in view of the above, and an object of the present invention is to provide a temperature measuring device and a temperature measuring method that enable temperature measurement with high measurement accuracy over a long period of time. [Means for solving the problem]

[0025] The present invention provides a temperature measuring device comprising a control unit that acquires a Brillouin frequency shift from a measurement position of a measured optical fiber, and when a predetermined exposure time, which is the time during which the measurement position of the measured optical fiber is exposed to a predetermined temperature environment, has elapsed, corrects the Brillouin frequency shift in response to a drift to the high frequency side based on the exposure time, and calculates a measured temperature at the measurement position based on the corrected Brillouin frequency shift. Effect of the Invention

[0026] According to the present invention, there is an advantage that temperature measurement can be performed with high measurement accuracy over a long period of time. [Brief description of the drawings]

[0027] [Figure 1] FIG. 2 is an explanatory diagram illustrating the time dependence of the Brillouin frequency shift. [Diagram 2] 1 is a block diagram illustrating an example of a functional configuration of a temperature measuring device according to an embodiment. [Diagram 3] FIG. 4 is an explanatory diagram illustrating an example of recording data. [Figure 4A] 4 is a flowchart illustrating an example of the operation of the temperature measuring device according to the embodiment. [Figure 4B] 4 is a flowchart illustrating an example of the operation of the temperature measuring device according to the embodiment. [Diagram 5] FIG. 1 is an explanatory diagram illustrating an overview of a conventional temperature measuring device. [Figure 6] FIG. 2 is an explanatory diagram illustrating the relationship between temperature and Brillouin frequency shift. [Figure 7] FIG. 2 is a diagram showing an example of the Brillouin frequency shift over the longitudinal distance of an optical fiber sensor. [Figure 8] FIG. 13 is a diagram showing an example of temperature versus distance in the longitudinal direction of an optical fiber sensor. [Figure 9]FIG. 2 is an explanatory diagram illustrating the change over time of the Brillouin frequency shift for each position of the optical fiber sensor. [Figure 10] FIG. 4 is an explanatory diagram illustrating the change in temperature over time at each position of the optical fiber sensor. [Figure 11] FIG. 13 is a diagram showing an experimental result of the change in Brillouin frequency shift over time when an optical fiber sensor is exposed to a constant temperature environment. DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS

[0028] Hereinafter, a temperature measuring device and a temperature measuring method according to an embodiment will be described with reference to the drawings. Configurations having the same functions in the embodiments are given the same reference numerals, and duplicated descriptions will be omitted. Note that the temperature measuring device and the temperature measuring method described in the following embodiments are merely examples, and do not limit the embodiments. In addition, the following embodiments may be appropriately combined within a range that does not cause contradiction.

[0029] The temperature measuring device according to the embodiment obtains the BFS at each measurement position in the measured optical fiber of the optical fiber sensor. Then, the temperature measuring device according to the embodiment calculates the temperature at each measurement position by applying the BFS-temperature conversion formula (1) to the obtained BFS.

[0030] Here, the time dependency of the Brillouin frequency shift, more specifically, the high-frequency drift of the BFS occurring at a certain exposure time is experimentally obtained in advance for each annealing temperature for the optical fiber under test. The annealing temperature refers to the heat treatment temperature of the optical fiber under test, i.e., the temperature to which the optical fiber under test is exposed. The annealing time refers to the time during which the optical fiber under test is treated at the annealing temperature, i.e., the time during which the optical fiber under test is exposed to the annealing temperature.

[0031] FIG. 1 is an explanatory diagram illustrating the time dependency of the Brillouin frequency shift. As shown in FIG. 1, case C1 is an experimental result in which the annealing time when the annealing temperature Ta is XX degrees Celsius is plotted on the horizontal axis, and the Brillouin frequency shift (BFS) is plotted on the vertical axis. Case C2 is an experimental result in which the annealing time when the annealing temperature Ta is △△ degrees Celsius (〇〇<△△) is plotted on the horizontal axis, and the Brillouin frequency shift (BFS) is plotted on the vertical axis. In both cases C1 and C2, a high-frequency drift in the BFS occurs even though the annealing temperature is constant beyond a certain annealing time. In the following explanation, the annealing time at which such a high-frequency drift starts is referred to as the high-frequency drift start time t1.

[0032] As shown in cases C1 and C2, the high-frequency drift of the BFS increases with increasing annealing time after the high-frequency drift start time t1. Therefore, by using a regression calculation or the like that fits the experimental results (cases C1 and C2), formula (2) for calculating the BFS correction amount (ΔBFS) to reproduce the high-frequency drift according to the exposure time of the BFS after the high-frequency drift start time t1 is obtained for each annealing temperature.

[0033]

number

[0034] Here, in the BFS correction amount calculation formula (2), t is the exposure time, t1 is the high frequency drift start time, a and b are coefficients, and c is a constant. The high frequency drift start time t1, the coefficients a and b, and the constant c are different values ​​for each annealing temperature. Note that while the above calculation formula (2) is expressed as a quadratic polynomial, it goes without saying that a cubic, quartic, or other polynomial may also be used.

[0035] In the temperature measuring device according to the embodiment, the BFS after the high-frequency drift start time t1 is corrected (subtracted) according to the drift to the high-frequency side based on the exposure time by using a calculation formula for the BFS correction amount for each annealing temperature. This makes it possible for the temperature measuring device according to the embodiment to perform temperature measurement with high measurement accuracy over a long period of time.

[0036] 2 is a block diagram showing an example of a functional configuration of a temperature measuring device according to an embodiment of the present invention, the temperature measuring device 1 includes an optical fiber sensor 2 and an optical fiber characteristic measuring device 3. As shown in FIG.

[0037] The optical fiber sensor 2 comprises a measured optical fiber 10 used as a distributed temperature sensor, a light-emitting unit 11 that injects light (pump light, probe light) into the measured optical fiber 10, and a light-receiving unit 12 that receives scattered light due to Brillouin scattering generated at each measurement position by the light incident on the measured optical fiber 10 by the light-emitting unit 11.

[0038] The optical fiber characteristic measuring device 3 acquires the BFS due to Brillouin scattering occurring at each measurement position of the measured optical fiber 10 based on the detection result of the light receiving unit 12, and calculates the temperature at each measurement position by applying the BFS-temperature conversion formula (1) to the acquired BFS.

[0039] The storage unit 30 is a non-volatile memory or the like, and stores setting information 30a, recorded data 30b, etc. The setting information 30a is data describing various setting contents including a BFS-temperature conversion formula (1), a calculation formula (2) for a BFS correction value for each temperature (annealing temperature), a high-frequency drift start time for each temperature (annealing temperature), etc.

[0040] The recorded data 30b is data in which the history of the measurement results (BFS, measured temperature) at each measurement position of the measured optical fiber 10 is listed in chronological order with, for example, a time stamp (date and time of measurement) added.

[0041] Fig. 3 is an explanatory diagram for explaining an example of the recorded data 30b. As shown in Fig. 3, the recorded data 30b records the measurement date and time [year / month / day / time], BFS [Hz], and measurement temperature [°C] at the measurement positions AAA.A [m], BBB.B [m]...XXX.X [m] of the measured optical fiber 10. That is, the recorded data 30b records the exposure time (total of measurement dates and times) exposed to the temperature environment of the measurement temperature [°C].

[0042] The control unit 31 is, for example, a computer realized by a processor, and includes a BFS measurement unit 31a, a BFS correction unit 31b, a temperature calculation unit 31c, and a recording unit 31d.

[0043] The BFS measurement unit 31a is a processing unit that measures (acquires) the Brillouin frequency shift (BFS) at each measurement position of the measured optical fiber 10. Specifically, the BFS measurement unit 31a inputs pump light and probe light for each measurement position from the light emitter 11 into the measured optical fiber 10. Next, the BFS measurement unit 31a acquires the Brillouin frequency shift (BFS) by receiving scattered light due to Brillouin scattering generated at each measurement position with the light receiver 12.

[0044] The BFS correction unit 31b is a processing unit that corrects the Brillouin frequency shift (BFS) at each measurement position acquired by the BFS measurement unit 31a based on the exposure time recorded in the recording data 30b. Specifically, when the exposure time recorded in the recording data 30b is shorter than the high-frequency drift start time t1 set in the setting information 30a, the BFS correction unit 31b performs correction according to the low-frequency drift that occurs according to the exposure time.

[0045] The low-frequency drift is a phenomenon caused by a change in the glass structure of the measured optical fiber 10, called structural relaxation under high temperature conditions. The change in the glass structure of the measured optical fiber 10 changes the physical properties of the glass, and as a result, a low-frequency drift of the BFS occurs over the course of exposure time even though the annealing temperature is constant. Since this low-frequency drift is known as a function of exposure time, the BFS correction unit 31b uses this function to perform correction according to the low-frequency drift.

[0046] The low-frequency drift is a phenomenon caused by a change in the glass structure of the measured optical fiber 10, called structural relaxation, in a high-temperature environment. The change in the glass structure of the measured optical fiber 10 changes the physical properties of the glass, and as a result, a low-frequency drift in the BFS occurs over exposure time even if the annealing temperature is constant. The BFS correction unit 31b predicts the amount of change in the BFS due to structural relaxation of the measured optical fiber 10 from the thermal history of the measured optical fiber 10, and corrects the BFS acquired by the BFS measurement unit 31a using the predicted amount of change.

[0047] The BFS correction unit 31b predicts the amount of change in BFS due to the structural relaxation of the optical fiber 10 under test described above, using a fictive temperature prediction formula that indicates the change over time in the fictive temperature of the optical fiber 10 under test. Here, the fictive temperature is a temperature that indicates the temperature at which the glass structure corresponds to the stable structure of a supercooled liquid. The annealing time is t, and the change over time in the fictive temperature is T f (t), the change in fictive temperature over time T f (t) can be expressed by the following equation (3).

[0048]

number

[0049] In addition, in equation (3), τ1 is the relaxation time, α and β are coefficients, and T f0is the fictive temperature in the initial state of the measured optical fiber 10. The initial state of the measured optical fiber 10 is the state before the measured optical fiber 10 is subjected to an annealing process.

[0050] Here, the change in the fictive temperature over time when the test optical fiber 10 is annealed is expressed as the fictive temperature T f0 and the annealing temperature Ta. Therefore, the change in fictive temperature over time T f (t) is the fictive temperature T f0 and the annealing temperature Ta.

[0051] Change in BFS over time BFS(t) and change in fictive temperature over time T f The relationship with (t) is expressed by the following equation (4): In the following equation (4), k is a coefficient, and n is a constant.

[0052]

number

[0053] The change amount ΔBFS(t) of the BFS when the annealing time t has elapsed since the start of the annealing process on the optical fiber is expressed by the following formula (5). f0 is the fictive temperature at the start of the annealing process.

[0054]

number

[0055] Here, if it is considered that the annealing treatment of the measured optical fiber 10 is performed by the heat of the object to be measured, the above formula (5) can be said to be a formula showing the amount of change ΔBFS in BFS when the exposure time t during which the measurement position of the measured optical fiber 10 is exposed to a predetermined temperature environment has elapsed. Therefore, the amount of change ΔBFS in BFS of the measured optical fiber 10 can be predicted from the above formula (5). In this embodiment, the BFS correction unit 31b predicts the amount of change ΔBFS in BFS using the above formula (5), and corrects the BFS acquired by the BFS measurement unit 31a using the predicted amount of change ΔBFS in BFS.

[0056] Furthermore, when the exposure time recorded in the recorded data 30b is equal to or longer than the high-frequency drift start time t1 set in the setting information 30a, the BFS correction unit 31b performs correction according to the high-frequency drift occurring according to the exposure time by using the calculation formula (2) for the BFS correction value for each temperature (annealing temperature) included in the setting information 30a. That is, when the exposure time becomes equal to or longer than the high-frequency drift start time t1, the BFS correction unit 31b switches from correction according to the low-frequency drift to correction according to the high-frequency drift.

[0057] The temperature calculation unit 31c is a processing unit that calculates the temperature at each measurement position by applying the BFS-temperature conversion formula (1) to the BFS acquired by the BFS measurement unit 31a and corrected by the BFS correction unit 31b.

[0058] The recording unit 31d is a processing unit that records, as the recorded data 30b, the history of the measurement results (BFS, measurement temperature, measurement date and time) at each measurement position of the measured optical fiber 10. Specifically, the recording unit 31d records the BFS at each measurement position of the measured optical fiber 10 acquired by the BFS measurement unit 31a and the measurement temperature calculated by the temperature calculation unit 31c, with a timestamp (measurement date and time) or the like, as the recorded data 30b. Note that the recording unit 31d may include the BFS correction value calculated by the BFS correction unit 31b using the calculation formula (2) or the like in the history of the measurement results at each measurement position of the measured optical fiber 10.

[0059] 4A and 4B are flowcharts showing an example of the operation of the temperature measuring device according to the embodiment. As shown in Fig. 4A and 4B, when the process starts, the recording unit 31d starts measuring the exposure time based on the timestamp (measurement date and time) of the recorded data 30b (S11).

[0060] Next, the BFS measurement unit 31a measures the BFS at each measurement position of the measured optical fiber 10 (S12). Next, the temperature calculation unit 31c calculates the temperature at each measurement position from the measured BFS by using the BFS-temperature conversion formula (1) (S13).

[0061] Next, the temperature calculation unit 31c selects a virtual temperature prediction formula from the calculated temperature (S14). This virtual temperature prediction formula is a formula for selecting a virtual temperature corresponding to the calculated temperature, and the virtual temperature prediction formula for each temperature is included in advance in the setting information 30a.

[0062] Next, the recording unit 31d records various data such as the BFS measured at each measurement position and the calculated temperature in the record data 30b with a time stamp (date and time of measurement) (S15).

[0063] Next, the BFS measuring unit 31a measures the BFS at each measurement position of the measured optical fiber 10 (S16). Next, the temperature calculating unit 31c calculates the fictive temperature using the selected fictive temperature prediction formula (S17).

[0064] Next, the BFS correction unit 31b calculates a BFS correction value according to the low-frequency drift described above from the calculated virtual temperature (S18). Next, the BFS correction unit 31b corrects the BFS measured in S16 using the calculated BFS correction value (S19).

[0065] Next, the temperature calculation unit 31c calculates the temperature from the corrected BFS by using the BFS-temperature conversion formula (1) (S20). Next, the recording unit 31d records various data such as the BFS measured at each measurement position, the BFS correction value, and the calculated temperature in the recorded data 30b with a timestamp (measurement date and time) (S21).

[0066] Next, the control unit 31 refers to the recorded data 30b and calculates the difference between the current temperature and the previous temperature at each measurement position (S22). Next, the control unit 31 determines whether or not there is a temperature change point (measurement position) based on the calculated difference (S23).

[0067] If there is a temperature change point (S23: Yes), the temperature calculation unit 31c selects a new virtual temperature prediction formula based on the current temperature (S24) and proceeds to S25. If there is no temperature change point (S23: No), the temperature calculation unit 31c skips S24 and proceeds to S25.

[0068] In S25, the control unit 31 determines whether or not to end the temperature measurement based on the presence or absence of an end instruction from the user, etc. (S25) If the temperature measurement is to be ended (S25: Yes), the control unit 31 ends the process.

[0069] If the temperature measurement is not to be ended (S25: No), the BFS correction unit 31b refers to the recorded data 30b and determines whether the exposure time t at each measurement position is equal to or greater than the (high frequency) drift start time t1 (S100).

[0070] Here, the BFS correction unit 31b may obtain the exposure time t at each measurement position as the total of the measurement times. The BFS correction unit 31b may also obtain the exposure time t at each measurement position as the total of the times during which a predetermined temperature (e.g., XXX°C, △△°C) was measured. The high-frequency drift start time t1 may be a unique value (e.g., a value for each temperature) included in the setting information 30a. For example, the BFS correction unit 31b may obtain the high-frequency drift start time t1 corresponding to the temperature used to obtain the exposure time t from the values ​​for each temperature included in the setting information 30a.

[0071] If the exposure time t is not equal to or greater than the high-frequency drift start time t1 (S100: No), the control unit 31 returns the process to S16.

[0072] If the exposure time t is equal to or greater than the high-frequency drift start time t1 (S100: Yes), the BFS measurement unit 31a measures the BFS at each measurement position of the measured optical fiber 10 (S101). Next, the BFS correction unit 31b calculates a BFS correction value for the high-frequency drift occurring according to the exposure time by using the calculation formula (2) for the BFS correction value based on the exposure time recorded in the recording data 30b. Here, the calculation formula (2) is set in advance for each temperature (annealing temperature) in the setting information 30a of the storage unit 30. The BFS correction unit 31b selects the calculation formula (2) corresponding to the current temperature from the calculation formulas (2) for the BFS correction value for each temperature (annealing temperature) included in the setting information 30a. Next, the BFS correction unit 31b corrects the BFS using the calculated BFS correction value (S102).

[0073] Next, the temperature calculation unit 31c calculates the temperature from the corrected BFS by using the BFS-temperature conversion formula (1) (S103). Next, the recording unit 31d records each piece of data, such as the BFS measured at each measurement position, the BFS correction value, and the calculated temperature, in the recorded data 30b with a timestamp (measurement date and time) (S104).

[0074] Next, the control unit 31 refers to the recorded data 30b and calculates the difference between the current temperature and the previous temperature at each measurement position (S105). Next, the control unit 31 determines whether or not there is a temperature change point (measurement position) that is equal to or greater than the design value in the calculation formula (2) based on the calculated difference (S106).

[0075] If there is a temperature change point (S106: Yes), the temperature calculation unit 31c ends the measurement of the temperature change position (S107) and proceeds to S108. If there is no temperature change point (S106: No), the temperature calculation unit 31c skips S107 and proceeds to S108.

[0076] In S108, the control unit 31 determines whether to continue the temperature measurement based on the presence or absence of an end instruction from the user or the like (S108). Note that the determination in S108 may be that the temperature measurement should not be continued if there is a temperature change position where the measurement has ended.

[0077] If the temperature measurement is to be continued (S108: Yes), the control unit 31 returns the process to S101. If the temperature measurement is not to be continued (S108: No), the control unit 31 ends the measurement.

[0078] As described above, the temperature measuring device 1 measures the exposure time during which the measurement position of the measured optical fiber 10 is exposed to a predetermined temperature environment. The temperature measuring device 1 acquires the Brillouin frequency shift (BFS) from the measurement position of the measured optical fiber 10. When the exposure time has elapsed for a predetermined time (time when the high frequency drift starts), the temperature measuring device 1 corrects the Brillouin frequency shift (BFS) in response to the drift to the high frequency side based on the exposure time. The temperature measuring device 1 calculates the temperature at the measurement position based on the Brillouin frequency shift (BFS). This allows the temperature measuring device 1 to respond to the high frequency drift of the BFS that occurs after the predetermined exposure time, and enables temperature measurement with high measurement accuracy over a long period of time.

[0079] Furthermore, the temperature measuring device 1 corrects the Brillouin frequency shift in response to the drift to the low frequency side based on the exposure time until the exposure time t has elapsed for a predetermined time (the high frequency drift start time). This allows the temperature measuring device 1 to respond to the low frequency drift of the BFS that occurs with the exposure time until the high frequency drift starts, and enables temperature measurement with high measurement accuracy.

[0080] Furthermore, the temperature measuring device 1 records in the record data 30b the exposure time and temperature at each of the multiple measurement positions of the measured optical fiber 10. This allows the temperature measuring device 1 to measure the history of the exposure time and temperature at each of the multiple measurement positions of the measured optical fiber 10.

[0081] Furthermore, the temperature measuring device 1 corrects the Brillouin frequency shift (BFS) in response to the drift to the high frequency side based on the exposure time and temperature when the recorded exposure time has elapsed a high frequency side drift start time based on the recorded temperature history for each of the multiple measurement positions of the measured optical fiber 10. This allows the temperature measuring device 1 to more accurately correct the high frequency side drift of the BFS based on the temperature history recorded for each of the multiple measurement positions of the measured optical fiber 10.

[0082] The temperature measuring device 1 also calculates ΔBFS, which is the correction value of the Brillouin frequency shift, based on formula (2) where t is the exposure time, t1 is the high-frequency drift start time based on temperature, and a, b, and c are coefficients according to temperature. This enables the temperature measuring device 1 to more accurately calculate the correction value of the Brillouin frequency shift.

[0083] In addition, the temperature measuring device 1 may use the measured optical fiber 10 with a slow high-frequency drift start time t1 as a sensor. As a result, if the life of the sensor due to the light-emitting unit 11 and the light-receiving unit 12 is equal to or shorter than the high-frequency drift start time t1, there is no need to correct the measured BFS, and temperature measurement with high measurement accuracy can be performed for a long period of time. Even if the life of the sensor is not equal to or shorter than the high-frequency drift start time t1, the period for correction can be shortened.

[0084] Furthermore, the temperature measuring device 1 may use the optical fiber 10 under test whose core is not doped with Ge. For example, in a high-temperature environment, it is considered that the diffusion of Ge doped in the core into the cladding causes a decrease in the Young's modulus of the core, which is a cause of high-frequency drift. Therefore, by using the optical fiber 10 under test whose core is not doped with Ge, the high-frequency drift may not occur in the temperature measuring device 1, and it is possible to perform temperature measurement with high measurement accuracy for a long period of time without the need to correct the measured BFS.

[0085] Furthermore, a prediction formula may be created in advance that predicts the BFS time change of the high frequency side drift of the measured optical fiber 10 for each annealing temperature from the relationship between the annealing temperature and the diffusion coefficient of Ge. By using this prediction formula, the temperature measuring device 1 can predict the BFS time change of the high frequency side drift when the annealing temperature is low, which is difficult to obtain by experiment, for example.

[0086] Some examples of combinations of the disclosed technical features are set out below.

[0087] (1) Obtaining the Brillouin frequency shift from a measurement position on the optical fiber under test; when a predetermined exposure time has elapsed, the Brillouin frequency shift is corrected in response to a drift toward a high frequency side based on the exposure time, the exposure time being a time during which the measurement position of the optical fiber to be measured is exposed to a predetermined temperature environment; calculating a measured temperature at the measurement position based on the corrected Brillouin frequency shift; A temperature measuring device comprising a control unit.

[0088] (2) The control unit recording the exposure time and the measured temperature for each of a plurality of measurement locations of the test optical fiber; The temperature measuring device according to (1) above.

[0089] (3) the predetermined time is a time determined by the predetermined temperature environment, The control unit is correcting the Brillouin frequency shift for each of a plurality of measurement positions having different predetermined temperature environments of the optical fiber to be measured, based on the predetermined temperature environment, the predetermined time, and the exposure time; The temperature measuring device according to (2) above.

[0090] (4) The control unit A correction value ΔBFS of the Brillouin frequency shift is calculated based on Equation (2) in which t is the exposure time, t1 is the predetermined time based on the predetermined temperature environment, and a, b and c are coefficients corresponding to the predetermined temperature environment. A temperature measuring device according to any one of (1) to (3).

[0091] (5) The control unit obtaining the formula (2) corresponding to the predetermined temperature environment from the formulas (2) for the predetermined temperature environments stored in a storage unit, and calculating the ΔBFS; The temperature measuring device according to (4) above.

[0092] (6) The control unit correcting the Brillouin frequency shift in response to a drift toward a lower frequency side based on the exposure time until the exposure time has elapsed for a predetermined time; A temperature measuring device according to any one of (1) to (3).

[0093] (7) The control unit When the exposure time has elapsed for a predetermined time, the correction of the Brillouin frequency shift corresponding to the drift toward the low frequency side is switched to the correction of the Brillouin frequency shift corresponding to the drift toward the high frequency side. The temperature measuring device according to (6) above.

[0094] (8) The control unit A correction value ΔBFS of the Brillouin frequency shift is calculated based on Equation (5) where t is the exposure time and k is a coefficient according to the predetermined temperature environment. The temperature measuring device according to (6) above.

[0095] (9) The computer: Obtaining the Brillouin frequency shift from the measurement position of the optical fiber under test; when a predetermined exposure time has elapsed, the Brillouin frequency shift is corrected in response to a drift toward a high frequency side based on the exposure time, the exposure time being a time during which the measurement position of the optical fiber to be measured is exposed to a predetermined temperature environment; calculating a measured temperature at the measurement position based on the corrected Brillouin frequency shift; A temperature measurement method comprising:

[0096] (10) The computer, recording the exposure time and the measured temperature for each of a plurality of measurement locations of the test optical fiber; The temperature measuring method according to (9) above.

[0097] (11) The predetermined time is a time determined by the predetermined temperature environment, The computer includes: correcting the Brillouin frequency shift for each of a plurality of measurement positions having different predetermined temperature environments of the optical fiber to be measured, based on the predetermined temperature environment, the predetermined time, and the exposure time; The temperature measuring method according to (10) above.

[0098] (12) The computer, A correction value ΔBFS of the Brillouin frequency shift is calculated based on Equation (2) in which t is the exposure time, t1 is the predetermined time based on the predetermined temperature environment, and a, b and c are coefficients corresponding to the predetermined temperature environment. A temperature measuring method according to any one of (9) to (11).

[0099] (13) The computer obtaining the formula (2) corresponding to the predetermined temperature environment from the formulas (2) for the predetermined temperature environments stored in a storage unit, and calculating the ΔBFS; The temperature measuring method according to (12) above.

[0100] (14) The computer, correcting the Brillouin frequency shift in response to a drift toward a lower frequency side based on the exposure time until the exposure time has elapsed for a predetermined time; A temperature measuring method according to any one of (9) to (11).

[0101] (15) The computer When the exposure time has elapsed for a predetermined time, the correction of the Brillouin frequency shift corresponding to the drift toward the low frequency side is switched to the correction of the Brillouin frequency shift corresponding to the drift toward the high frequency side. The temperature measuring method according to (14) above.

[0102] (16) The computer, A correction value ΔBFS of the Brillouin frequency shift is calculated based on Equation (5) where t is the exposure time and k is a coefficient according to the predetermined temperature environment. The temperature measuring method according to (14) above. [Explanation of symbols]

[0103] 1...Temperature measuring device 2, 101...Optical fiber sensor 3, 100...Optical fiber characteristic measuring device 10...Optical fiber to be measured 11…Light emitting part 12...Light receiving section 30...Storage section 30a…Setting information 30b…Recorded data 31...Control unit 31a...BFS measurement section 31b…BFS correction section 31c...Temperature calculation section 31d…Recording section 102, 103…position C1, C2…Case R1, R2…Experimental results

Claims

1. Obtaining the Brillouin frequency shift from the measurement position of the optical fiber under test; when a predetermined exposure time has elapsed, the Brillouin frequency shift is corrected in response to a drift toward a high frequency side based on the exposure time, the exposure time being a time during which the measurement position of the optical fiber to be measured is exposed to a predetermined temperature environment; calculating a measured temperature at the measurement position based on the corrected Brillouin frequency shift; A temperature measuring device comprising a control unit.

2. The control unit is recording the exposure time and the measured temperature for each of a plurality of measurement locations of the test optical fiber; 2. The temperature measuring device according to claim 1 .

3. the predetermined time is a time determined by the predetermined temperature environment, The control unit is correcting the Brillouin frequency shift for each of a plurality of measurement positions having different predetermined temperature environments of the optical fiber to be measured, based on the predetermined temperature environment, the predetermined time, and the exposure time; 3. The temperature measuring device according to claim 2.

4. The control unit is A correction value ΔBFS of the Brillouin frequency shift is calculated based on Equation (1) in which t is the exposure time, t1 is the predetermined time based on the predetermined temperature environment, and a, b and c are coefficients corresponding to the predetermined temperature environment. [0010] 4. The temperature measuring device according to claim 1,

5. The control unit is obtaining the formula (1) corresponding to the predetermined temperature environment from the formulas (1) for the predetermined temperature environments stored in a storage unit, and calculating the ΔBFS; 5. The temperature measuring device according to claim 4.

6. The control unit is correcting the Brillouin frequency shift in response to a drift toward a lower frequency side based on the exposure time until the exposure time has elapsed for a predetermined time; 4. The temperature measuring device according to claim 1,

7. The control unit is When the exposure time has elapsed for a predetermined time, the correction of the Brillouin frequency shift corresponding to the drift toward the low frequency side is switched to the correction of the Brillouin frequency shift corresponding to the drift toward the high frequency side.

7. The temperature measuring device according to claim 6.

8. The control unit is A correction value ΔBFS of the Brillouin frequency shift is calculated based on the formula (2) in which the exposure time is t and the coefficient according to the predetermined temperature environment is k. [0025] 7. The temperature measuring device according to claim 6.

9. The computer Obtaining the Brillouin frequency shift from the measurement position of the optical fiber under test; when a predetermined exposure time has elapsed, the Brillouin frequency shift is corrected in response to a drift toward a high frequency side based on the exposure time, the exposure time being a time during which the measurement position of the optical fiber to be measured is exposed to a predetermined temperature environment; calculating a measured temperature at the measurement position based on the corrected Brillouin frequency shift; A temperature measurement method comprising:

Images