WO2010125712A1 - 温度測定システム及び温度測定方法 - Google Patents

温度測定システム及び温度測定方法 Download PDF

Info

Publication number
WO2010125712A1
WO2010125712A1 PCT/JP2009/071302 JP2009071302W WO2010125712A1 WO 2010125712 A1 WO2010125712 A1 WO 2010125712A1 JP 2009071302 W JP2009071302 W JP 2009071302W WO 2010125712 A1 WO2010125712 A1 WO 2010125712A1
Authority
WO
WIPO (PCT)
Prior art keywords
temperature
optical fiber
correction
measurement
temperature distribution
Prior art date
Application number
PCT/JP2009/071302
Other languages
English (en)
French (fr)
Japanese (ja)
Inventor
丈夫 笠嶋
和史 宇野
文雄 武井
Original Assignee
富士通株式会社
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 富士通株式会社 filed Critical 富士通株式会社
Priority to CN200980158413.XA priority Critical patent/CN102365536B/zh
Priority to JP2011511264A priority patent/JP5218648B2/ja
Priority to EP09844055.5A priority patent/EP2431719B1/de
Publication of WO2010125712A1 publication Critical patent/WO2010125712A1/ja
Priority to US13/229,933 priority patent/US8636408B2/en

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01K—MEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
    • G01K11/00—Measuring temperature based upon physical or chemical changes not covered by groups G01K3/00, G01K5/00, G01K7/00 or G01K9/00
    • G01K11/32—Measuring 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
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01K—MEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
    • G01K15/00—Testing or calibrating of thermometers
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01K—MEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
    • G01K11/00—Measuring temperature based upon physical or chemical changes not covered by groups G01K3/00, G01K5/00, G01K7/00 or G01K9/00
    • G01K11/32—Measuring 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/324—Measuring 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

Definitions

  • the present invention relates to a temperature measurement system and a temperature measurement method.
  • a method of cooling an electronic device for example, there are a method of controlling an air flow in a room in order to prevent an excessive accumulation of heat, and a method of cooling a room by air conditioning.
  • the indoor cooling state In order to always keep the indoor cooling state optimal by these methods, it is effective to constantly monitor the indoor thermometer side point. In particular, if the three-dimensional temperature distribution in the room can be grasped, the temperature distribution can be used to control air conditioning and the like, and the indoor cooling state can be kept optimal.
  • the purpose is to measure temperature with high accuracy in the temperature measurement system and temperature measurement method.
  • a laser light source that outputs laser light, an optical fiber that is laid in a temperature measurement area and that receives the laser light, and backscattered light of the laser light that has entered the optical fiber
  • a temperature measurement unit that detects the temperature distribution of the temperature in the temperature measurement area along the optical fiber laying path and calculates the corrected temperature distribution by correcting the measurement temperature distribution,
  • the temperature measurement unit convolves the transfer function of the optical fiber along the laying path with the correction temperature distribution, and a square error between the measurement temperature distribution and the measurement temperature distribution is reduced with respect to the correction
  • laser light is incident on an optical fiber laid in a temperature measurement area, backscattered light of the laser light in the optical fiber is detected, and the optical fiber
  • the correction of the measured temperature distribution is sequentially performed so that the convolution of the transfer function of the optical fiber along the laying path and the corrected temperature distribution and the square error between the measured temperature distribution are reduced every time correction is performed.
  • There is provided a temperature measurement method that is performed a plurality of times and replaces the temperature after correction at a specific point of the laying path with the estimated temperature at the specific point for each correction.
  • the measurement temperature distribution is corrected sequentially, and the corrected temperature at a specific point of the optical fiber installation path is corrected at each specific point for each correction. It became clear that the corrected temperature distribution approaches the actual temperature distribution by substituting the estimated temperature.
  • FIG. 1 is a schematic diagram showing the structure of a computer room to which this embodiment is applied.
  • FIG. 2 is a cross-sectional view showing an example of laying optical fibers.
  • FIG. 3 is a schematic diagram showing the configuration of the temperature measurement system according to the present embodiment.
  • FIG. 4 is a diagram illustrating a spectrum of light backscattered by the optical fiber.
  • FIG. 5 is a diagram illustrating an example of a time series distribution of the intensity of Raman scattered light.
  • FIG. 6 is a diagram obtained by calculating the I 1 / I 2 ratio for each time based on the time-series distribution of the intensity of Raman scattered light in FIG.
  • FIG. 7 is a diagram illustrating a measured temperature distribution obtained using the thermometer side system according to the present embodiment.
  • FIG. 1 is a schematic diagram showing the structure of a computer room to which this embodiment is applied.
  • FIG. 2 is a cross-sectional view showing an example of laying optical fibers.
  • FIG. 3 is a schematic diagram
  • FIG. 8 is a transfer function of the thermometer side system obtained from the step-type actual temperature distribution.
  • FIG. 9 is a diagram illustrating a function obtained by performing a Fourier transform on the transfer function.
  • FIG. 10 is a diagram showing the measured temperature distribution obtained by measuring the actual temperature distribution in which the temperature changes at a relatively high spatial frequency with the temperature measurement system according to the present embodiment.
  • FIG. 11 is a diagram illustrating an example of an inverse filter used for correcting the measured temperature distribution.
  • FIG. 12 is a diagram showing a corrected temperature distribution obtained by correcting the measured temperature distribution by applying an inverse filter to the measured temperature distribution.
  • FIG. 13 is a cross-sectional view showing an example of laying an optical fiber useful for correcting the measured temperature distribution in the present embodiment.
  • FIG. 14 shows an embodiment in which the section of ⁇ 1 m to 1 m of the optical fiber is heated to 55 ° C. with the heating center as a reference, and the temperature of the other section is maintained at room temperature (about 23 ° C.). It is a figure which shows the measurement temperature distribution obtained by the temperature measurement system which concerns on.
  • FIG. 15 is a schematic diagram showing the laying length of the optical fiber in consideration of the influence from the heat source in the present embodiment.
  • FIG. 16 is a flowchart illustrating the temperature measurement method according to the present embodiment.
  • FIG. 17 is a diagram illustrating a measured temperature distribution obtained by the temperature measurement system according to the present embodiment.
  • FIG. 18 is a diagram illustrating a corrected temperature distribution when the measured temperature distribution is corrected once by the temperature measuring system according to the present embodiment.
  • FIG. 19 is a diagram illustrating a corrected temperature distribution when the temperature at a specific point is replaced after the correction is performed once by the temperature measurement system according to the present embodiment.
  • FIG. 20 is a diagram illustrating a corrected temperature distribution when the measured temperature distribution is corrected 100 times by the temperature measurement system according to the present embodiment.
  • FIG. 21 is a diagram illustrating a correction result according to the comparative example.
  • FIG. 22 is a schematic diagram illustrating an installed state of the optical fiber according to the first embodiment.
  • FIG. 23 is a diagram illustrating a temperature distribution measured by a thermocouple, a temperature distribution measured by an optical fiber temperature measurement system (before correction), and a temperature distribution after correction according to the present embodiment in Example 1.
  • FIG. 24 is a schematic diagram illustrating an installed state of the optical fiber in the second embodiment.
  • FIG. 25 is a diagram illustrating a temperature distribution measured by a thermocouple, a temperature distribution measured by an optical fiber temperature measurement system (before correction), and a temperature distribution after correction according to the present embodiment in
  • FIG. 1 is a schematic diagram showing the structure of a computer room to which this embodiment is applied.
  • the computer room is roughly divided into an equipment installation area 10 and a free access floor 15.
  • each rack 11 stores a plurality of computers such as blade servers. Moreover, between each server rack 11, the management space used at the time of management of the channel
  • the free access floor 15 is provided under the floor of the equipment installation area 10.
  • various cables 16 such as power cables and communication cables connected to the racks 11 are arranged in a state of being housed in a cable duct 17.
  • thermometer-side points are provided in the equipment installation area 10 and the free access floor 15, the temperature at each point is constantly monitored, and the air volume of the air conditioner 19 is controlled in real time according to the obtained temperature distribution. It is preferable to do this.
  • thermometer side method for measuring the temperature in real time there is a method of measuring the temperature of the optical fiber at multiple points from the intensity change of Raman scattered light generated when an infrared laser pulse is incident on the optical fiber. is there.
  • FIG. 2 is a sectional view showing an example of laying the optical fiber.
  • the optical fiber 24 is introduced into the server rack 11 and in the vicinity thereof from the free access floor 15 that is constantly maintained at a substantially constant temperature by cold air.
  • the optical fiber 24 is laid so as to pass through the vicinity of the intake port and the vicinity of the exhaust port of the rack 11 which is a point at which the temperature is to be measured.
  • FIG. 3 is a schematic diagram showing the configuration of the temperature measurement system 20 using the optical fiber 24.
  • FIG. 4 is a diagram showing a spectrum of light backscattered by the optical fiber 24.
  • thermometer side system 20 includes a laser light source 21, lenses 22a and 22b, a beam splitter 23, an optical fiber 24, a wavelength separator 25, a photodetector 26, and a temperature measurement. Part 27.
  • Laser light with a predetermined pulse width is output from the laser light source 21 at a constant cycle. This laser light enters the optical fiber 24 from the light source side end of the optical fiber 24 through the lens 22a, the beam splitter 23, and the lens 22b.
  • reference numeral 24 a indicates the cladding of the optical fiber 24, and reference numeral 24 b indicates the core of the optical fiber 24.
  • the backscattered light includes Rayleigh scattered light, Brillouin scattered light, and Raman scattered light.
  • Rayleigh scattered light is light having the same wavelength as incident light
  • Brillouin scattered light and Raman scattered light are light having wavelengths shifted from the incident wavelength.
  • Raman scattered light includes Stokes light shifted to a longer wavelength side than incident light and anti-Stokes light shifted to a shorter wavelength side than incident light.
  • the shift amount of the Stokes light and the anti-Stokes light depends on the wavelength of the laser light and the material of the optical fiber 24, but is usually about 50 nm.
  • the amount of change with temperature of the Stokes light is small, and the amount of change with temperature of the anti-Stokes light is large. That is, it can be said that the Stokes light has a small temperature dependency, and the anti-Stokes light has a large temperature dependency.
  • the wavelength separation unit 25 includes beam splitters 31a, 31b, and 31c that transmit or reflect light according to the wavelength, and optical filters 33a, 33b, and 33c that transmit only light of a specific wavelength. Further, the wavelength separation unit 25 includes condensing lenses 34a, 34b, and 34c that condense the light transmitted through the optical filters 33a, 33b, and 33c onto the light receiving units 26a, 26b, and 26c of the photodetector 26, respectively.
  • the light incident on the wavelength separator 25 is separated into Rayleigh scattered light, Stokes light, and anti-Stokes light by the beam splitters 31a, 31b, 31c and the optical filters 33a, 33b, 33c. , 26c.
  • signals corresponding to the intensity of Rayleigh scattered light, Stokes light, and anti-Stokes light are output from the light receiving sections 26a, 26b, and 26c.
  • the temperature measuring unit 27 such as a personal computer acquires a measured temperature distribution along the laying path of the optical fiber 24 based on the signal output from the photodetector 26, and the measured temperature distribution is described later. Correct as follows.
  • the pulse width of the backscattered light input to the photodetector 26 is related to the length of the optical fiber 24. Therefore, the interval between the laser pulses output from the laser light source 21 is set so that the backscattered light from each laser pulse does not overlap. In addition, if the power of the laser beam is too high, there is a possibility that a stimulated Raman scattering state occurs and correct measurement cannot be performed. Therefore, it is preferable to control the power of the laser light source 21 so that the stimulated Raman scattering state does not occur.
  • the temperature at the position where the backscattering can be evaluated by the ratio between the two.
  • the intensity ratio of Stokes light and anti-Stokes light is as follows when the angular frequency of incident light is ⁇ 0 , the angular frequency of optical phonons in the optical fiber is ⁇ k , the Planck constant is h, the Boltzmann constant is k, and the temperature is T. It is expressed by the following equation (1).
  • the temperature at the position where backscattering can be calculated from the equation (1).
  • the backscattered light generated in the optical fiber 24 is attenuated while returning through the optical fiber 24. Therefore, in order to correctly evaluate the temperature at the position where backscattering occurs, it is necessary to consider the attenuation of light.
  • FIG. 5 is a diagram showing an example of a time-series distribution of the intensity of Raman scattered light, with time on the horizontal axis and the signal intensity output from the light receiving unit of the photodetector on the vertical axis.
  • Stokes light and anti-Stokes light are detected by the photodetector 26 for a certain period immediately after the laser pulse is incident on the optical fiber 24.
  • the signal intensity decreases with the passage of time when the laser pulse is incident on the optical fiber 24 as a reference.
  • the time on the horizontal axis indicates the distance from the light source side end of the optical fiber 24 to the position where the backscattering occurs, and the decrease in signal intensity with time indicates the attenuation of light by the optical fiber 24. .
  • the signal intensity of Stokes light and anti-Stokes light is not attenuated uniformly. As shown in FIG. 5, peaks and valleys appear on the curve indicating the change in signal intensity with time.
  • the intensity of anti-Stokes light at a certain time t is I 1
  • the intensity of Stokes light is I 2 .
  • the temperature distribution in the length direction of the optical fiber 24 can be measured by calculating the intensity ratio (I 1 / I 2 ) between the anti-Stokes light and the Stokes light.
  • thermometer side system 20 using such an optical fiber will be examined.
  • FIG. 7 is a diagram showing a measured temperature distribution obtained using the thermometer side system 20.
  • the stepped actual temperature distribution rising from room temperature to 55 ° C. was given by immersing the predetermined portion 24a of the optical fiber 24 in hot water of 55 ° C.
  • the predetermined portion 24a has three lengths of 0.5 m, 1.0 m, and 2.0 m.
  • each measured temperature distribution has a rounded shape obtained by applying a weighted moving average to the actual temperature distribution. From this, it can be understood that the thermometer side system 20 has only a low spatial frequency response, that is, the position resolution is poor.
  • FIG. 8 is a transfer function h of the thermometer-side system 20 obtained from the step-type actual temperature distribution of FIG. 7, where the horizontal axis indicates the distance from the heating center, and the vertical axis indicates the relative intensity of the temperature.
  • the power spectrum of the function g has a very small value in a region where the spatial frequency is about 0.6 m ⁇ 1 or more. From this, it can be seen that the thermometer side system 20 functions as a low-pass filter that cuts off a region where the spatial frequency is about 0.6 m ⁇ 1 or more, and loses most of the frequency information in this region.
  • the actual temperature changes in a long cycle along the optical fiber laying path, so that a highly accurate position resolution is not required for the temperature measurement system.
  • FIG. 10 is a diagram showing the measured temperature distribution obtained by measuring the actual temperature distribution in which the temperature changes at a relatively high spatial frequency with the temperature measurement system 20.
  • the actual temperature distribution in FIG. 10 is a measured value distribution of the temperature by the thermocouple.
  • the measured temperature distribution is shaped like a weighted moving average of the actual temperature distribution through a low-pass filter.
  • the measured temperature distribution of the temperature measurement system 20 cannot be used as it is, and the measured temperature distribution is appropriately corrected to approach the actual temperature distribution. I understand that it is necessary.
  • FIG. 11 is a diagram showing the characteristics of such an inverse filter.
  • this inverse filter is designed on the assumption that the noise amplification is reduced when there is noise in the measured temperature distribution in the region where the spatial frequency is 0.6 m ⁇ 1 or higher.
  • FIG. 12 is a diagram showing a corrected temperature distribution obtained by correcting the measured temperature distribution by applying this inverse filter to the measured temperature distribution of FIG.
  • the frequency component in the spatial frequency domain lost by the action of the low-pass filter is smaller than the power spectrum of the frequency component in the spatial frequency domain included in the noise during measurement, the frequency component is good for the inverse filter Cannot recover.
  • the measured temperature distribution is corrected as follows and brought closer to the actual temperature distribution.
  • FIG. 13 is a cross-sectional view showing an example of laying the optical fiber 24 useful for correcting the measured temperature distribution.
  • the first winding portion 24x of the optical fiber 24 and the second winding are provided on the free access floor 15 where the temperature is kept constant by the cold air under the floor 12.
  • a portion 24y is provided.
  • a third winding portion 24z of the optical fiber 24 is provided in the vicinity of the exhaust port 11a of the server rack 11, a third winding portion 24z of the optical fiber 24 is provided.
  • the exhaust port 11a is a part where the air heated by the computer that has generated heat is exhausted, and the third winding part 24z is heated to a temperature higher than room temperature by the hot air.
  • the optical fiber 24 is laid so that the round trip path to the third winding portion 24z is the same path.
  • each winding part 24x, 24y, 24z is not particularly limited, but the lower limit is preferably twice the minimum bending radius (about 15 mm) allowed for the optical fiber 24.
  • the upper limit of the diameter of each of the winding parts 24x, 24y, and 24z is preferably set to a diameter that allows the winding parts to be accommodated in a region that can be regarded as the same spatial temperature, for example, 45 mm. If it does in this way, it can be considered that the area wound by each winding part 24x, 24y, 24z among the installation paths
  • the optical fiber 24 in the section wound around the first winding part 24x and the second winding part 24y is held at the temperature of the free access floor.
  • the optical fiber 24 of the area wound by the 3rd winding part 24z is hold
  • the length of the section wound around the first winding portion 24x and the second winding portion 24y in the laying path of the optical fiber 24 is not particularly limited, but in the present embodiment, the length is as follows. Determine the length of the section.
  • FIG. 14 shows the temperature measurement when the section of ⁇ 1 m to 1 m of the optical fiber 24 is heated to 55 ° C. with the heating center as the reference, and the temperature of the other section is maintained at room temperature (about 23 ° C.).
  • FIG. 3 is a diagram showing a measured temperature distribution obtained by the system 20.
  • the measured temperature distribution has a base portion outside the heated section of ⁇ 1 m to 1 m, and the measured temperature of the portion is not equal to the actual room temperature. This is because if there is a temperature difference between two points on the optical fiber installation path, the temperature measurement result at each point is affected by the temperature difference.
  • the winding unit 24x if more than the absolute value of the zero point X 3 the length of the section to be wound on each 24y, a heat source outside the installation path of the section Even if there is, there is a temperature measurement point corresponding to the actual temperature at the winding parts 24x and 24y that is not affected by the heat source.
  • FIG. 15 is a schematic diagram showing the laying length of the optical fiber 24 in consideration of the influence from such a heat source.
  • the length of the optical fiber 24 between the adjacent server racks 11 is D 1
  • the optical fiber 24 extends from the winding portions 24 x and 24 y to the region above the floor 12 that is a higher temperature portion than the free access floor 15.
  • a length of up to enter are the D 2.
  • the heat source is a computer in the server rack 11. Moreover, it can be considered that the area G in the free access floor 15 in the laying path of the optical fiber 24 is kept at a constant temperature by the cold air.
  • the section G is assigned to each winding section 24x, 24y, the starting point of the section G is the floor 12, and the end point is the middle point P of the adjacent rack 11.
  • the length L of the optical fiber 24 in that section G is each turn 24x, the length of the optical fiber 24 of the wound portion 24y When D 3, the D 1/2 + D 2 + D 3. If this length L is set to be equal to or greater than the absolute value of the zero point X 3 of the transfer function h described above, a temperature measurement point that is not affected by the heat of the computer in the server rack 11, for example, the midpoint P exists in the section G. To do.
  • the section G when there are a plurality of temperature measurement points in the section G, it can be estimated that one of the measurement points, for example, the temperature at the midpoint P and the temperature of the remaining measurement points in the section G are the same. .
  • the measurement temperature distribution is corrected by utilizing the identity of the measurement temperature in such a section G.
  • Each length D 1 , D 2 , D 3 is particularly suitable if the length L of the optical fiber 24 in the section G is equal to or greater than the absolute value (3.3 m) of the zero point X 3 of the transfer function h. It is not limited. In this example, by setting D 1 to 1.0 m, D 2 to 0.5 m, and D 3 to 2.3 m, the length L is set to 3.3 m, and the length L is set to the zero X 3 of the transfer function h. The absolute value of (3.3 m) or more.
  • the optical fiber 24 is laid so that the reciprocating path to the third winding portion 24z is the same path.
  • any point of the optical fiber 24 in the third winding portion 24z can be regarded as the overlapping point K i having substantially the same temperature, and the correction temperature at each overlapping point K i is A condition that the temperature is the same can be added.
  • the temperature measurement area by the temperature measurement system 20 is the equipment installation area 10 and the free access floor 15, and the emphasis is particularly on temperature measurement in the server rack 11 and in the vicinity thereof.
  • FIG. 16 is a flowchart showing the temperature measurement method according to this embodiment. Each step in the flowchart is performed in the temperature measuring unit 27 described above.
  • the temperature was measured at a plurality of temperature measurement points of the optical fiber 24, and the distance between each temperature measurement point was set to 0.1 m.
  • some of the temperature measurement points are provided with thermocouples for measuring the actual temperature.
  • the measured temperature distribution obtained by the temperature measuring system 20 deviates from the actual temperature distribution by the thermocouple.
  • the measured temperature distribution is corrected as follows to bring it closer to the actual distribution.
  • the subscript k in the component y k represents a measurement point along the optical fiber laying path
  • the component y k is a temperature measurement value in a region where there is no temperature change from the temperature measurement value at the measurement point k (see FIG. 15).
  • a value obtained by subtracting the temperature value T AB ) at the midpoint P is shown.
  • the subscript i in the component x i represents the measurement point
  • the component x i is the temperature measurement value in the region where the temperature does not change from the actual temperature at the measurement point i (in the example shown in FIG. 15).
  • a value obtained by subtracting the temperature value T AB ) at the midpoint P is shown.
  • the measured temperature distribution y can be expressed by the following equation (4) as a convolution of the actual temperature distribution x and the transfer function h.
  • the range of i is a range satisfying that the subscript k-i is 0 or more.
  • each component h ij of the transfer function can be calculated using the least square method or the like, with the equation (5) as a simultaneous equation for h j .
  • the actual temperature distribution x and the measured temperature distribution y when obtaining each component h ij of the transfer function for example, a step-type actual temperature distribution as shown in FIG. 14 and a measured temperature distribution corresponding thereto are used. Can do.
  • the transfer function h changes according to the distance from the light source because the optical fiber 24 has a group delay characteristic. Therefore, the transfer function h cannot be uniquely defined over the entire length of the optical fiber 24. However, in the short section of the optical fiber 24, the loss and delay of the optical signal in the optical fiber 24 are considered to be uniform, and the transfer function h can be uniquely defined in the section.
  • the transfer function h differs depending not only on the distance from the light source but also on the material of the optical fiber 24, the pulse waveform of the incident laser, and the pulse response characteristics of the photodetector 26. Therefore, when obtaining each component h ij of the transfer function, it is preferable to obtain each component h ij according to the same conditions as when actually measuring the temperature.
  • each component 0 in the region where there is no temperature change is meaningless which is not necessary for the calculation, and thus is 0 before and after the attention region where the temperature change exists from each component in Equation (2).
  • a column vector obtained by collecting only the components excluding all components is expressed as the following equation (7).
  • the numbers of column vector components in equations (6) and (7) are n + 1 and m + 1, respectively, but for m and n, m is larger than n. This is because, as shown in FIG. 14, the measured temperature distribution spreads more in the horizontal direction than the actual temperature distribution, and therefore the number of non-zero components is larger in the measured temperature distribution.
  • Equation (8) when the actual temperature distribution x and the measured temperature distribution y are finite-dimensional column vectors, and equation (5) is expressed in the following equation (8), [H ] Is configured based on the transfer function h and has (m + 1) ⁇ (n + 1) finite number of components. [H] configured in this way is referred to as a transfer function matrix display.
  • each component y i of y is m + 1 values obtained by temperature measurement, and [H] can be regarded as a (m + 1) ⁇ (n + 1) counting matrix of simultaneous equations. Since there is a relationship of m> n as described above, this simultaneous equation cannot be uniquely solved for x.
  • column vector X in the equation (9) is an n-dimensional vector having a component as in the following equation (10), like the actual temperature distribution.
  • Equation (9) The distribution X for reducing e in Equation (9) approximately satisfies Equation (8). Then, as e in Equation (9) becomes smaller, the accuracy of approximation increases, and the distribution X approaches the actual temperature distribution x.
  • the distribution X is referred to as a corrected temperature distribution of the measured temperature distribution y.
  • equation (9) can be said to be an equation for calculating the square error e between the convolution of the transfer function h of the optical fiber 24 along the laying path and the corrected temperature distribution X and the measured temperature distribution y. .
  • the gradient vector ⁇ e / ⁇ X of the square error e is calculated from the following equation (11) from the equation (9).
  • the least square method is to determine each component X i of X so that the gradient vector ⁇ e / ⁇ X becomes zero.
  • Equation (11) If the diagonal component of [H] t [H] in Equation (11) is slightly increased in consideration of the noise at the time of measurement, amplification of the high frequency component of noise can be suppressed and margin tolerance can be increased.
  • the above-described correction by the inverse filter corresponds to the correction calculated by the least square method.
  • the square error e decreases if the direction of the reverse sign - ⁇ e / ⁇ X is advanced.
  • k indicates the number of correction iterations
  • X (k) is a corrected temperature distribution when the correction is performed k times.
  • the component of X (k) can be expressed as in the following equation (13).
  • ⁇ is a positive correction count such that the expression (12) converges and can be selected in the range of 0.5 to 1 empirically. In the following calculation, ⁇ is set to 0.5.
  • Equation (11) is used.
  • the iterative calculation is performed using the equation (12), whereby the correction temperature distribution X (k + 1) where the square error e becomes smaller than X (k) is sequentially calculated a plurality of times. .
  • the temperatures of the plurality of measurement points i in the section G in the laying path of the optical fiber 24 are the same as the temperatures at the midpoint P. .
  • the component X i corresponding to the plurality of measurement points i included in the first and second winding portions 24x and 24y in the section G is obtained each time the correction calculation is performed according to Expression (12).
  • each component of the column vectors x, y, and X is obtained by subtracting the measured temperature T AB at the midpoint P from the actual value.
  • T AB is the value of the component at the measurement point corresponding to the middle point P among the plurality of components y i of the measured temperature distribution y. It has significance as a common estimated temperature in G.
  • an average value of a plurality of components y i of a measurement temperature distribution y corresponding to a plurality of measurement points i in a region where the actual temperature included in each winding part 24x, 24y is measured is an estimated temperature common to the measurement points It is good.
  • each turn 24x, winding portion of the length D 3 by longer than previously described 2.3m of optical fibers 24 of the wound portion 24y 24x, actual temperature comprised 24y is The number of measurement points i corresponding to the region to be measured is increased, and the temperature estimation accuracy is improved. Furthermore, by increasing the length D 3 in this way, it is possible to make it less susceptible to the temperature of the adjacent server rack 11.
  • the corrected temperatures X i1 (k) and X i2 (k) are the measurement points i 1 and x 2 corresponding to the overlapping points H 1 and H 2 among the plurality of components X i (k) of the corrected temperature distribution X (k) . It is the value of the component at i2, and the average value X avg1 has a significance as an estimated temperature common to the overlapping points H 1 and H 2 .
  • the component X i of the measured temperature distribution at each overlapping point K i is obtained every time correction calculation according to the equation (12) is performed. (k) is replaced with the average value X avg2 of the correction temperature X i (k) at each of these overlapping points K i .
  • the correction temperature X i (k) is the value of the component at the measurement point i corresponding to each overlapping point K i among the plurality of components X i (k) of the correction temperature distribution X (k). It is.
  • the average value X avg1 has significance as an estimated temperature common to the overlapping points K i .
  • the interval between the temperature measurement points of the optical fiber 24 is set to 0.1 m, so that the length of the optical fiber 24 in the portion wound around the third winding portion 24z is 0.
  • the distance is 0.5 m
  • each winding 24x as described above, a plurality of measurement points i in each 24y can be regarded as the temperature at the same temperature T AB.
  • Equation (6) has a value as mentioned, the components of the column vector x, y, X is obtained by subtracting the T AB from the value of the actual temperature when defining the (7). Therefore, in order to obtain the final corrected temperature distribution T IOMP _ i is again the temperature T AB as the following equation (14) after the calculation of the required number of iterations (n times) has been completed for formula (12) Should be added.
  • the first and second winding portions 24x and 24y are prevented from rising in temperature by the cold air, while the other portions are floors that are hotter than the free access floor 15. Therefore, the temperature does not fall below the winding portions 24x and 24y.
  • the component X i (k) is set to 0, and then the k + 1-th calculation is performed.
  • the corrected temperature of the portion using Equation (12) is lower than the predetermined temperature. In some cases, it is preferable to simplify the calculation by replacing the corrected temperature in the portion with the predetermined temperature.
  • the corrected temperature of the part using Expression (12) is higher than the predetermined temperature. Even when the temperature becomes higher, it is preferable to replace the corrected temperature in this portion with the predetermined temperature.
  • step S2 correction calculation is repeatedly performed using the equation (12) as described above, and an index for a reduction amount of the square error e, for example, e (n) -e (n-1) is less than a predetermined value.
  • the final corrected temperature distribution T iomp — i is obtained from X i (n) .
  • e (n) is a square error obtained from the equation (9) using X (n) obtained by performing the correction by the equation (12) n times.
  • the corrected temperature distribution X i (k) at a specific point on the laying path of the optical fiber 24 is replaced with a predetermined value for each correction using the equation (12).
  • the plurality of overlapping points K i in the third winding part 24z are also replaced with the average value X avg2 .
  • 18 to 20 are diagrams showing the measured temperature distribution after correction.
  • the horizontal axis indicates the distance from the end of the optical fiber 24, and the vertical axis indicates the temperature.
  • FIG. 18 is a diagram when the correction by the equation (12) is performed once. As shown in FIG. 18, the deviation between the actual temperature distribution and the corrected temperature distribution has not been eliminated by one correction.
  • the corrected temperature distribution of the replaced portion has a temperature value that is substantially the same as the actual distribution.
  • FIG. 20 is a diagram when such correction calculation is repeated 100 times. As shown in FIG. 20, when the correction calculation is performed 100 times, the corrected temperature distribution substantially matches the actual temperature distribution.
  • FIG. 21 is a diagram illustrating a correction result according to the comparative example.
  • the correction by the equation (12) is performed 100 times without replacing the corrected temperature distribution X i (k) at a specific point on the laying path of the optical fiber 24 with a predetermined value. .
  • FIG. 21 it can be seen that the deviation from the actual temperature distribution is not eliminated in the comparative example.
  • the correction temperature distribution X i (k) at a specific point on the laying path of the optical fiber 24 is replaced with a predetermined value every time correction calculation is performed according to Expression (12). A close correction temperature distribution can be obtained.
  • Example 1 As shown in FIG. 22, an optical fiber 24 was laid on the door on the exhaust side of the server rack 11. That is, the winding portions 24x and 24y of the optical fiber 24 are arranged at the same location in the free access floor, and the optical fiber 24 between the winding portions 24x and 24y is drawn into the rack 11 and along the exhaust door from the bottom to the top. Laying back and forth.
  • reference numeral 241 indicates an optical fiber on the forward path side
  • reference numeral 242 indicates an optical fiber on the return path side.
  • a plurality of thermocouples were arranged along the optical fiber 24, and the temperature distribution of the rack 11 and the free access floor was also measured using these thermocouples.
  • the horizontal axis indicates the position in the length direction of the optical fiber 24, and the vertical axis indicates the temperature distribution measured by the thermocouple, and the temperature distribution measured by the optical fiber temperature measurement system 20 (before correction). And a temperature distribution after correction according to the present embodiment.
  • the correction using the equation (12) is repeated 100 times, and the temperature of the measurement point in the region corresponding to the winding portions 24x and 24y is measured corresponding to the midpoint P in the section G of FIG. Replaced with point temperature.
  • the corrected temperature distribution almost coincides with the temperature distribution measured by the thermocouple.
  • the corrected temperature distribution can capture the difference in temperature distribution between the forward path 241 and the return path 242 of the optical fiber.
  • Example 2 As shown in FIG. 24, an optical fiber 24 was laid on the door on the exhaust side of the server rack. That is, the winding portions 24x and 24y of the optical fiber 24 are arranged at the same location in the free access floor, and the optical fiber 24 between the winding portions 24x and 24y is drawn into the rack 11 and along the exhaust door from the bottom to the top. Laying back and forth. However, in the forward path, optical fiber winding portions A, B, C, and D are arranged at positions corresponding to exhaust ports of four servers (not shown) mounted in the rack. The return optical fiber is laid so as to pass through the same position as the forward optical fiber except for the winding portions A, B, C, and D and the vicinity thereof. Further, for comparison, a plurality of thermocouples are arranged along the optical fiber 24, and the temperature distribution of the rack 11 and the free access floor is also measured by these thermocouples.
  • FIG. 25 shows the temperature distribution measured by the thermocouple and the temperature distribution measured by the optical fiber temperature measurement system 20 (before correction), with the horizontal axis representing the position in the length direction of the optical fiber 24 and the vertical axis representing the temperature. And a temperature distribution after correction according to the present embodiment.
  • the correction using the equation (12) is repeated 100 times, and the temperature of the measurement point in the region corresponding to the winding portions 24x and 24y is measured corresponding to the midpoint P in the section G of FIG. Replaced with point temperature.
  • the corrected temperature distribution almost coincides with the temperature distribution measured by the thermocouple. Further, a wide range of temperature distributions including the winding parts A, B, C, D and their peripheral parts can be captured.

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Measuring Temperature Or Quantity Of Heat (AREA)
PCT/JP2009/071302 2009-05-01 2009-12-22 温度測定システム及び温度測定方法 WO2010125712A1 (ja)

Priority Applications (4)

Application Number Priority Date Filing Date Title
CN200980158413.XA CN102365536B (zh) 2009-05-01 2009-12-22 温度测定系统以及温度测定方法
JP2011511264A JP5218648B2 (ja) 2009-05-01 2009-12-22 温度測定システム及び温度測定方法
EP09844055.5A EP2431719B1 (de) 2009-05-01 2009-12-22 System und verfahren zur temperaturmessung
US13/229,933 US8636408B2 (en) 2009-05-01 2011-09-12 Temperature measurement system and temperature measurement method

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2009-112156 2009-05-01
JP2009112156 2009-05-01

Related Child Applications (1)

Application Number Title Priority Date Filing Date
US13/229,933 Continuation US8636408B2 (en) 2009-05-01 2011-09-12 Temperature measurement system and temperature measurement method

Publications (1)

Publication Number Publication Date
WO2010125712A1 true WO2010125712A1 (ja) 2010-11-04

Family

ID=43031872

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2009/071302 WO2010125712A1 (ja) 2009-05-01 2009-12-22 温度測定システム及び温度測定方法

Country Status (5)

Country Link
US (1) US8636408B2 (de)
EP (1) EP2431719B1 (de)
JP (1) JP5218648B2 (de)
CN (1) CN102365536B (de)
WO (1) WO2010125712A1 (de)

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2012056567A1 (ja) * 2010-10-29 2012-05-03 富士通株式会社 温度測定システム及び温度測定方法
WO2012081092A1 (ja) * 2010-12-15 2012-06-21 富士通株式会社 空調制御システム及び風量調整装置
JP2013502006A (ja) * 2009-08-12 2013-01-17 インターナショナル・ビジネス・マシーンズ・コーポレーション データ・センタに関する知識ベースのモデル化のための方法、製品、および装置
WO2013024525A1 (ja) * 2011-08-15 2013-02-21 富士通株式会社 温度分布測定装置及び温度分布測定方法
CN103115696A (zh) * 2013-01-22 2013-05-22 山东大学 消除瑞利光串扰的分布式光纤拉曼测温系统及测温方法
WO2014064845A1 (ja) 2012-10-26 2014-05-01 富士通株式会社 温度測定システム及び異常検知方法
WO2014064771A1 (ja) 2012-10-23 2014-05-01 富士通株式会社 異常検知システム及び異常検知方法
WO2014064770A1 (ja) 2012-10-23 2014-05-01 富士通株式会社 異常検知システム及び異常検知方法
CN103782144A (zh) * 2011-08-31 2014-05-07 富士通株式会社 温度分布测定系统、温度分布测定装置以及温度分布测定方法
JP2015040643A (ja) * 2013-08-20 2015-03-02 富士通株式会社 氷床温度管理システム及び氷床温度管理方法
WO2015132899A1 (ja) * 2014-03-05 2015-09-11 富士通株式会社 異常検知システム及び異常検知方法
JP2019090696A (ja) * 2017-11-15 2019-06-13 横河電機株式会社 温度測定システムおよび温度測定方法
US10422751B2 (en) 2013-09-24 2019-09-24 Fujitsu Limited Optical fiber cord and abnormality detection system
US10704964B2 (en) 2015-05-13 2020-07-07 Fujitsu Limited Temperature measurement device, temperature measurement method, and computer-readable non-transitory medium

Families Citing this family (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3035021A1 (de) * 2012-12-11 2016-06-22 Fujitsu Limited Temperaturmesssystem und temperaturmessverfahren
JP5742861B2 (ja) * 2013-02-28 2015-07-01 横河電機株式会社 光ファイバ温度分布測定装置
CN103364658A (zh) * 2013-06-28 2013-10-23 国网电力科学研究院武汉南瑞有限责任公司 基于光纤光栅测温系统进行变压器寿命预测的方法
WO2015047267A1 (en) 2013-09-26 2015-04-02 Intel Corporation Methods of forming dislocation enhanced strain in nmos structures
CN103606240B (zh) * 2013-11-27 2016-06-29 宁波振东光电有限公司 采用分布式光纤温度传感器系统进行火灾报警的方法
CN104089652B (zh) * 2014-07-14 2016-09-21 国家电网公司 一种光纤光栅变压器在线监测系统及其监测方法
US10861682B2 (en) 2014-07-31 2020-12-08 iSenseCloud, Inc. Test wafer with optical fiber with Bragg Grating sensors
JP2016053443A (ja) * 2014-09-03 2016-04-14 富士通株式会社 温度分布予測方法および空調管理システム
JP6428336B2 (ja) * 2015-02-12 2018-11-28 富士通株式会社 温度測定システム、温度測定方法及びプログラム
JP6428350B2 (ja) * 2015-02-18 2018-11-28 富士通株式会社 温度測定システム、温度測定方法及びプログラム
CN104792436B (zh) * 2015-04-10 2017-09-08 青岛创立科技开发有限公司 一种分布式光纤测温方法
CN105466594B (zh) * 2015-11-20 2018-01-30 国网上海市电力公司 一种变压器热点温度最小二乘法修正方法
CN107462899B (zh) * 2017-06-30 2019-08-16 华中科技大学鄂州工业技术研究院 一种大气环境参数的测量方法及系统
US11899080B2 (en) 2017-09-06 2024-02-13 United States Department Of Energy Optical sensing of magnetic fields using magnetically sensitive materials
US11067458B2 (en) * 2017-09-06 2021-07-20 Carnegie Mellon University Distributed sensing of electromagnetic components using optical fiber based methods
CN110081998A (zh) * 2018-01-25 2019-08-02 博立信（洛阳）科技有限公司 温度采集系统
CN109280726B (zh) * 2018-09-18 2021-08-27 安徽工业大学 一种基于多元线性回归算法预测炼铁高炉炉芯死料柱温度的方法
CN110579249B (zh) * 2019-09-17 2021-06-29 西北大学 一种基于掺钴多模光纤光栅的热线式流量传感器及制造方法
CN111811694B (zh) * 2020-07-13 2021-11-30 广东博智林机器人有限公司 一种温度校准方法、装置、设备及存储介质

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH05172656A (ja) * 1991-12-24 1993-07-09 Asahi Glass Co Ltd 分布型光ファイバーセンサー
JPH06109557A (ja) 1992-09-30 1994-04-19 Toshiba Corp 測定装置およびその装置に用いるシース形光ファイバの製造法
JPH1011681A (ja) * 1996-06-21 1998-01-16 Yamato Protec Co 設備異常警報方法及び設備異常警報システム
JP3065832B2 (ja) * 1993-01-06 2000-07-17 株式会社東芝 温度分布検出装置
JP2004028748A (ja) 2002-06-25 2004-01-29 Hitachi Cable Ltd 温度測定方法及び温度測定装置

Family Cites Families (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4576486A (en) * 1983-08-23 1986-03-18 The United States Of America As Represented By The Secretary Of Commerce Optical fiber thermometer
JPS61133826A (ja) * 1984-12-04 1986-06-21 Mitsubishi Electric Corp 光計測装置
US5178465A (en) * 1990-07-11 1993-01-12 Fujikura Ltd. Optical fiber laying structure for electric power cable line trouble occurrence location detecting system
US5308162A (en) * 1992-02-13 1994-05-03 Fujikura Ltd. Temperature abnormality detecting structure for fluid pipe
KR0133488B1 (en) * 1993-01-06 1998-04-23 Toshiba Kk Temperature distribution detector using optical fiber
US5638396A (en) * 1994-09-19 1997-06-10 Textron Systems Corporation Laser ultrasonics-based material analysis system and method
FR2764376B1 (fr) * 1997-06-04 1999-07-16 Commissariat Energie Atomique Dispositif integre de lecture des raies spectrales contenues dans un spectre optique
US20050149264A1 (en) * 2003-12-30 2005-07-07 Schlumberger Technology Corporation System and Method to Interpret Distributed Temperature Sensor Data and to Determine a Flow Rate in a Well
CN100439860C (zh) * 2004-06-25 2008-12-03 光纳株式会社 分布式光纤传感器
US7480460B2 (en) * 2005-03-29 2009-01-20 University Of New Brunswick Dynamic strain distributed fiber optic sensor
US7628531B2 (en) * 2006-03-13 2009-12-08 SensorTran, Inc Methods and apparatus for dual source calibration for distributed temperature systems
WO2007108239A1 (ja) 2006-03-20 2007-09-27 Sharp Kabushiki Kaisha 表示装置
US7478552B2 (en) * 2006-03-21 2009-01-20 Veeco Instruments Inc. Optical detection alignment/tracking method and apparatus
JP4775173B2 (ja) * 2006-08-24 2011-09-21 住友電気工業株式会社 光ファイバ温度センサ

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH05172656A (ja) * 1991-12-24 1993-07-09 Asahi Glass Co Ltd 分布型光ファイバーセンサー
JPH06109557A (ja) 1992-09-30 1994-04-19 Toshiba Corp 測定装置およびその装置に用いるシース形光ファイバの製造法
JP3065832B2 (ja) * 1993-01-06 2000-07-17 株式会社東芝 温度分布検出装置
JPH1011681A (ja) * 1996-06-21 1998-01-16 Yamato Protec Co 設備異常警報方法及び設備異常警報システム
JP2004028748A (ja) 2002-06-25 2004-01-29 Hitachi Cable Ltd 温度測定方法及び温度測定装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
"Fujitsu Develops Technology Enabling Real-time Multiple-Point Temperature Measurement", 4 April 2008, FUJITSU LABORATORIES LTD.

Cited By (28)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2013502006A (ja) * 2009-08-12 2013-01-17 インターナショナル・ビジネス・マシーンズ・コーポレーション データ・センタに関する知識ベースのモデル化のための方法、製品、および装置
JPWO2012056567A1 (ja) * 2010-10-29 2014-03-20 富士通株式会社 温度測定システム及び温度測定方法
WO2012056567A1 (ja) * 2010-10-29 2012-05-03 富士通株式会社 温度測定システム及び温度測定方法
JP5664658B2 (ja) * 2010-10-29 2015-02-04 富士通株式会社 温度測定システム及び温度測定方法
JP5601379B2 (ja) * 2010-12-15 2014-10-08 富士通株式会社 空調制御システム
WO2012081092A1 (ja) * 2010-12-15 2012-06-21 富士通株式会社 空調制御システム及び風量調整装置
JPWO2012081092A1 (ja) * 2010-12-15 2014-05-22 富士通株式会社 空調制御システム及び風量調整装置
WO2013024525A1 (ja) * 2011-08-15 2013-02-21 富士通株式会社 温度分布測定装置及び温度分布測定方法
US9797782B2 (en) 2011-08-15 2017-10-24 Fujitsu Limited Temperature distribution measurement apparatus and temperature distribution measurement method
JP5673830B2 (ja) * 2011-08-15 2015-02-18 富士通株式会社 温度分布測定装置及び温度分布測定方法
JPWO2013024525A1 (ja) * 2011-08-15 2015-03-05 富士通株式会社 温度分布測定装置及び温度分布測定方法
US9488532B2 (en) 2011-08-31 2016-11-08 Fujitsu Limited Temperature distribution measurement system, temperature distribution measurement apparatus, and temperature distribution measurement method
CN103782144A (zh) * 2011-08-31 2014-05-07 富士通株式会社 温度分布测定系统、温度分布测定装置以及温度分布测定方法
WO2014064771A1 (ja) 2012-10-23 2014-05-01 富士通株式会社 異常検知システム及び異常検知方法
US9915620B2 (en) 2012-10-23 2018-03-13 Fujitsu Limited Abnormality detection system and abnormality detection method
WO2014064770A1 (ja) 2012-10-23 2014-05-01 富士通株式会社 異常検知システム及び異常検知方法
US9528860B2 (en) 2012-10-23 2016-12-27 Fujitsu Limited Abnormality detection system and abnormality detection method
US9347803B2 (en) 2012-10-23 2016-05-24 Fujitsu Limited Abnormality detection system and abnormality detection method
WO2014064845A1 (ja) 2012-10-26 2014-05-01 富士通株式会社 温度測定システム及び異常検知方法
US9816878B2 (en) 2012-10-26 2017-11-14 Fujitsu Limited Temperature measurement system and abnormality detection method
CN103115696A (zh) * 2013-01-22 2013-05-22 山东大学 消除瑞利光串扰的分布式光纤拉曼测温系统及测温方法
JP2015040643A (ja) * 2013-08-20 2015-03-02 富士通株式会社 氷床温度管理システム及び氷床温度管理方法
US10422751B2 (en) 2013-09-24 2019-09-24 Fujitsu Limited Optical fiber cord and abnormality detection system
WO2015132899A1 (ja) * 2014-03-05 2015-09-11 富士通株式会社 異常検知システム及び異常検知方法
JPWO2015132899A1 (ja) * 2014-03-05 2017-03-30 富士通株式会社 異常検知システム及び異常検知方法
US10088342B2 (en) 2014-03-05 2018-10-02 Fujitsu Limited Abnormality detection system and abnormality detection method
US10704964B2 (en) 2015-05-13 2020-07-07 Fujitsu Limited Temperature measurement device, temperature measurement method, and computer-readable non-transitory medium
JP2019090696A (ja) * 2017-11-15 2019-06-13 横河電機株式会社 温度測定システムおよび温度測定方法

Also Published As

Publication number Publication date
EP2431719A1 (de) 2012-03-21
US20120033709A1 (en) 2012-02-09
EP2431719B1 (de) 2018-07-04
EP2431719A4 (de) 2016-01-20
US8636408B2 (en) 2014-01-28
JPWO2010125712A1 (ja) 2012-10-25
JP5218648B2 (ja) 2013-06-26
CN102365536A (zh) 2012-02-29
CN102365536B (zh) 2014-06-11

Similar Documents

Publication Publication Date Title
JP5218648B2 (ja) 温度測定システム及び温度測定方法
JP5664658B2 (ja) 温度測定システム及び温度測定方法
JP5825352B2 (ja) 温度分布測定システム、温度分布測定装置及び温度分布測定方法
JP5975108B2 (ja) 温度測定システム及び異常検知方法
JP5673830B2 (ja) 温度分布測定装置及び温度分布測定方法
JP5613974B2 (ja) 温度測定方法
AU2019202787B2 (en) Temperature measurement device, temperature measurement method, and temperature measurement program
JP5228798B2 (ja) 温度計測システム及び温度計測方法
AU2015394726A1 (en) Temperature measurement device, temperature measurement method, and temperature measurement program
AU2015394728B2 (en) Temperature measurement device, temperature measurement method, and temperature measurement program
JP5136429B2 (ja) 温度測定システム
JP6791374B2 (ja) 温度測定装置、温度測定方法および温度測定プログラム
JPH0242179B2 (de)
JP2017110921A (ja) ケーブル診断システムおよびセンシング・ケーブル
JP6631175B2 (ja) 温度測定装置、温度測定方法および温度測定プログラム
CN114279594A (zh) 一种分布式光纤传感器高温动态标定系统及方法
JP6003296B2 (ja) 湿度検出装置及び空調システム
JPH04208877A (ja) 布設電力ケーブルジョイントの導体温度監視方法
Hayashi et al. Infrared thermometry for detecting flaws in large-core polymer optical fibers
CN114577367A (zh) 一种光纤温度传感器标定的方法、装置和计算机设备

Legal Events

Date Code Title Description
WWE Wipo information: entry into national phase

Ref document number: 200980158413.X

Country of ref document: CN

121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 09844055

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2011511264

Country of ref document: JP

Kind code of ref document: A

WWE Wipo information: entry into national phase

Ref document number: 2009844055

Country of ref document: EP

NENP Non-entry into the national phase

Ref country code: DE