[go: up one dir, main page]
More Web Proxy on the site http://driver.im/

JP2010223831A - Temperature measuring device - Google Patents

Temperature measuring device Download PDF

Info

Publication number
JP2010223831A
JP2010223831A JP2009072637A JP2009072637A JP2010223831A JP 2010223831 A JP2010223831 A JP 2010223831A JP 2009072637 A JP2009072637 A JP 2009072637A JP 2009072637 A JP2009072637 A JP 2009072637A JP 2010223831 A JP2010223831 A JP 2010223831A
Authority
JP
Japan
Prior art keywords
light
temperature
stokes light
distribution
stokes
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
JP2009072637A
Other languages
Japanese (ja)
Inventor
Nobunari Takeuchi
伸成 竹内
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Yokogawa Electric Corp
Original Assignee
Yokogawa Electric Corp
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 Yokogawa Electric Corp filed Critical Yokogawa Electric Corp
Priority to JP2009072637A priority Critical patent/JP2010223831A/en
Publication of JP2010223831A publication Critical patent/JP2010223831A/en
Pending legal-status Critical Current

Links

Images

Landscapes

  • Measuring Temperature Or Quantity Of Heat (AREA)

Abstract

<P>PROBLEM TO BE SOLVED: To provide a temperature measuring device capable of measuring temperature precisely considering distribution of Raman shift wavenumbers. <P>SOLUTION: The temperature measuring device 1 includes: photoelectric conversion circuits 21a, 21b for outputting photoelectric conversion signals by receiving Raman scattering light (Stokes light P<SB>s</SB>and anti-Stokes light P<SB>a</SB>) obtained by applying laser beams R<SB>L</SB>to an optical fiber 13; a memory 26 for storing a correction table considering distribution of the Raman shift wavenumbers of the Raman scattering light; and an arithmetic processing section 25 for performing prescribed operation for obtaining temperature of the optical fiber 13 by photoelectric conversion signals outputted from the photoelectric conversion circuits 21a, 21b assuming that the Raman scattering light does not have any distribution of Raman shift wavenumbers and correcting operation results obtained by the prescribed operation by the correction table stored in the memory 26. <P>COPYRIGHT: (C)2011,JPO&INPIT

Description

本発明は、被測定対象にレーザ光を入射させて得られるラマン散乱光を用いて被測定対象の温度を測定する温度測定装置に関する。   The present invention relates to a temperature measuring apparatus that measures the temperature of a measurement target using Raman scattered light obtained by making laser light incident on the measurement target.

近年、ラマン散乱光を用いて温度を測定する温度測定装置の研究開発が盛んに行われている。この温度測定装置の1つに、光ファイバ内で生ずるラマン散乱光(ストークス光及び反ストークス光)を検出して光ファイバの長さ方向における温度分布を測定するR−OTDR(Raman Optical Time Domain Reflectmetry)がある。具体的に、この温度測定装置は、光ファイバに入射させるパルス状のレーザ光を射出する光源部と、光ファイバにレーザ光を入射させて得られる後方散乱光に含まれるストークス光と反ストークス光とを個別に受光する受光部と、受光部から出力される受光信号を用いて光ファイバの長さ方向における温度分布を求める演算処理部とを備える。   In recent years, research and development of temperature measuring devices that measure temperature using Raman scattered light have been actively conducted. One of the temperature measuring devices is an R-OTDR (Raman Optical Time Domain Reflectmetry) that detects Raman scattered light (Stokes light and anti-Stokes light) generated in an optical fiber and measures the temperature distribution in the length direction of the optical fiber. ) Specifically, the temperature measuring device includes a light source unit that emits pulsed laser light that is incident on an optical fiber, and Stokes light and anti-Stokes light included in backscattered light that is obtained by allowing the laser light to enter the optical fiber. And a calculation processing unit for obtaining a temperature distribution in the length direction of the optical fiber using a light reception signal output from the light receiving unit.

ここで、光ファイバから射出されるストークス光及び反ストークス光は、その強度比が温度に比例して変化する。このため、演算処理部は、レーザ光を光ファイバに入射させた後のある時点で得られるストークス光の受光信号と反ストークス光の受光信号との強度比を求めることにより、光ファイバの長さ方向におけるある点の温度を測定している。また、レーザ光を光ファイバに入射させた時点以降のストークス光の受光信号と反ストークス光の受光信号との強度比の時間変化を求めることで、光ファイバの長さ方向における温度分布を測定している。   Here, the intensity ratio of the Stokes light and the anti-Stokes light emitted from the optical fiber changes in proportion to the temperature. For this reason, the arithmetic processing unit obtains the intensity ratio of the Stokes light reception signal and the anti-Stokes light reception signal obtained at a certain point after the laser light is incident on the optical fiber, thereby obtaining the length of the optical fiber. The temperature at a point in the direction is measured. In addition, the temperature distribution in the length direction of the optical fiber is measured by determining the temporal change in the intensity ratio between the received light signal of Stokes light and the received light signal of anti-Stokes light after the point where the laser light is incident on the optical fiber. ing.

光ファイバの長さ方向におけるある点の温度の具体的な求め方は以下の通りである。レーザ光を光ファイバに入射させた後のある時点で得られる反ストークス光の強度A及びストークス光の強度Sは、以下の(1)式で表される。

Figure 2010223831
The specific method for obtaining the temperature at a certain point in the length direction of the optical fiber is as follows. The intensity of anti-Stokes light obtained at some time after the laser beam is incident on the optical fiber A s and Stokes light intensities S t is expressed by the following equation (1).
Figure 2010223831

上記(1)式中におけるA,Sは温度測定装置の性能及び光ファイバの特性に依存する固有の定数であり、νは光源部から射出されるレーザ光(光ファイバに入射するレーザ光)の波数であり、νは光ファイバに入射するレーザ光の波数とストークス光又は反ストークス光の波数との差であるラマンシフト波数である。また、hはプランク定数であり、cは光速であり、kはボルツマン定数であって、Tは温度(絶対温度)である。上記(1)式を参照すると、反ストークス光の強度A及びストークス光の強度Sは何れも温度Tの関数で表されるのが分かる。 In the equation (1), A 0 and S 0 are specific constants depending on the performance of the temperature measuring device and the characteristics of the optical fiber, and ν 0 is a laser beam emitted from the light source unit (a laser incident on the optical fiber). ) 1 is a Raman shift wave number that is the difference between the wave number of laser light incident on the optical fiber and the wave number of Stokes light or anti-Stokes light. Further, h is a Planck constant, c is the speed of light, k is a Boltzmann constant, and T is a temperature (absolute temperature). Referring to equation (1) it is seen that both the intensity A s and Stokes light intensities S t of the anti-Stokes light is expressed by a function of the temperature T.

ストークス光の波長をλ、反ストークス光の波長をλとし、A=Sとすると、反ストークス光とストークス光との強度比R(T)=A/Sは以下の(2)式で表される。

Figure 2010223831
Assuming that the wavelength of the Stokes light is λ s , the wavelength of the anti-Stokes light is λ a, and A 0 = S 0 , the intensity ratio R (T) = A s / S t between the anti-Stokes light and the Stokes light is as follows ( 2) It is expressed by the formula.
Figure 2010223831

ここで、既知の温度Tにおける反ストークス光とストークス光との強度比をR(T)とすると、以下の(3)式を用いて温度Tを求めることができる。

Figure 2010223831
Here, when the intensity ratio between the anti-Stokes light and the Stokes light at a known temperature T r is R (T r ), the temperature T can be obtained using the following equation (3).
Figure 2010223831

つまり、光ファイバの既知の温度Tにおける反ストークス光とストークス光との強度比を予め求めておき、その光ファイバの未知の温度における反ストークス光とストークス光とを測定して強度比を求めれば、上記(3)式から未知の温度を求めることができる。尚、ラマンシフト波数νは光ファイバ固有の定数であり、他の手段によって予め測定しておく。レーザ光を光ファイバに入射させてからストークス光と反ストークス光との強度比を順次求めて上記(3)式に順次代入すれば、光ファイバの長さ方向における温度分布を測定することができる。尚、従来の温度測定装置の詳細については、例えば以下の特許文献1,2及び非特許文献1を参照されたい。 That is, the intensity ratio between the anti-Stokes light and the Stokes light at a known temperature Tr of the optical fiber is obtained in advance, and the intensity ratio can be obtained by measuring the anti-Stokes light and the Stokes light at the unknown temperature of the optical fiber. For example, an unknown temperature can be obtained from the above equation (3). The Raman shift wavenumber ν 1 is a constant inherent to the optical fiber and is measured in advance by other means. If the intensity ratio between Stokes light and anti-Stokes light is sequentially obtained after the laser light is incident on the optical fiber and is sequentially substituted into the above equation (3), the temperature distribution in the length direction of the optical fiber can be measured. . For details of the conventional temperature measuring device, refer to, for example, Patent Documents 1 and 2 and Non-Patent Document 1 below.

特開平6−26940号公報JP-A-6-26940 英国特許出願公開第2140554号明細書UK Patent Application No. 2140554

ジェイ・ピー・ダイキン、他3名(J. P. Dakin et al.),「ディストリビューティド・アンチストークス・レシオ・サーモメトリー(Distributed antistokes ratio thermometry)」,カンファレンスペーパー(Conference Paper),オプティカル・ファイバー・センサーズ(Optical Fiber Sensors),サンディエゴ・カリフォルニア(San Diego, CA),1985年1月1日,p.PDS3−1Jay Pea Daikin, 3 others (JP Dakin et al.), "Distributed antistokes ratio thermometry", Conference Paper, Optical Fiber Sensors (Optical Fiber Sensors), San Diego, CA, January 1, 1985, p. PDS3-1

ところで、従来の温度測定装置では、前述した(1)〜(3)式中におけるラマンシフト波数νを定数として取り扱っている。しかしながら、光ファイバ内で実際に生ずるストークス光及び反ストークス光の波数はある広がりを有しており、ラマンシフト波数も定数(線スペクトル)ではなくある分布を有している。このため、ラマンシフト波数を定数として取り扱う従来の温度測定装置では測定される温度に誤差が生ずるという問題があった。 By the way, in the conventional temperature measuring apparatus, the Raman shift wave number ν 1 in the above-described equations (1) to (3) is handled as a constant. However, the wave numbers of Stokes light and anti-Stokes light actually generated in the optical fiber have a certain spread, and the Raman shift wave number also has a distribution rather than a constant (line spectrum). For this reason, the conventional temperature measuring apparatus that handles the Raman shift wave number as a constant has a problem that an error occurs in the measured temperature.

ここで、前述した(2)式を変形すると、以下の(4)式で表される。

Figure 2010223831
上記(4)式に示される通り、ラマンシフト波数νは指数関数の指数内に現れるばかりではなく、4乗項の分母及び分子に現れる。このため、ラマンシフト波数νの分布によって反ストークス光とストークス光との強度比R(T)は大きく影響され、この結果として測定温度に誤差が生ずることが容易に分かる。 Here, when the above-described equation (2) is modified, it is expressed by the following equation (4).
Figure 2010223831
As shown in the above equation (4), the Raman shift wavenumber ν 1 appears not only in the exponent of the exponential function but also in the denominator and numerator of the fourth power term. Therefore, the intensity ratio between anti-Stokes light and the Stokes light by the distribution of Raman shift wavenumber ν 1 R (T) is greatly affected, it is easily seen that the error occurs in the result as a measurement temperature.

図8は、光ファイバ内で実際に生ずるストークス光の波数分布を示す図の一例である。尚、図8においては、横軸に波数をとり、縦軸に相対強度をとっている。尚、図8中の波数ν100は光ファイバに入射するレーザ光の波数である。図8に示す通り、光ファイバ内で実際に生ずるストークス光は、波数ν101〜ν102の範囲で図示の通りの分布を有する。ラマンシフト波数はストークス光の波数ν101〜ν102とレーザ光の波数ν100との差であるため、その分布は図8に示すストークス光の分布に応じた分布となる。尚、図8に示す例では、波数が44,000[1/m]のときにストークス光の強度がピークになり、ストークス光の分布の半値全幅は26,000[1/m]である。 FIG. 8 is an example of a diagram showing the wave number distribution of Stokes light actually generated in the optical fiber. In FIG. 8, the horizontal axis represents the wave number and the vertical axis represents the relative intensity. Note that the wave number ν 100 in FIG. 8 is the wave number of the laser light incident on the optical fiber. As shown in FIG. 8, the Stokes light actually generated in the optical fiber has a distribution as shown in the range of wave numbers ν 101 to ν 102 . Since the Raman shift wave number is the difference between the Stokes light wave numbers ν 101 to ν 102 and the laser light wave number ν 100 , the distribution is a distribution corresponding to the Stokes light distribution shown in FIG. 8. In the example shown in FIG. 8, the intensity of Stokes light peaks when the wave number is 44,000 [1 / m], and the full width at half maximum of the Stokes light distribution is 26,000 [1 / m].

図9は、従来の温度測定装置における温度測定誤差の一例を示す図である。尚、図8においては、横軸に温度をとり、縦軸に温度測定誤差をとっている。図9において、符号D101を付した曲線は、前述した(1)〜(3)式中のラマンシフト波数νに図8中のラマンシフト波数ν201を設定した場合の温度測定誤誤差を示すグラフであり、符号D102を付した曲線はラマンシフト波数νに図8中のラマンシフト波数ν202を設定した場合の温度測定誤誤差を示すグラフである。 FIG. 9 is a diagram illustrating an example of a temperature measurement error in a conventional temperature measurement device. In FIG. 8, the horizontal axis represents temperature, and the vertical axis represents temperature measurement error. In FIG. 9, a curve denoted by reference numeral D101 indicates a temperature measurement error error when the Raman shift wave number ν 201 in FIG. 8 is set to the Raman shift wave number ν 1 in the above-described equations (1) to (3). The curve with the symbol D102 is a graph showing a temperature measurement error when the Raman shift wave number ν 202 in FIG. 8 is set to the Raman shift wave number ν 1 .

図9を参照すると、光ファイバの温度が前述した(3)式中の既知の温度Tである場合には、グラフD101とグラフD102とが交差する。しかしながら、既知の温度Tから外れた温度においては、温度測定誤差がグラフD101とグラフD102とで異なる(ラマンシフト波数νの値によって異なる)ことが分かる。例えば、既知の温度Tが300[K]である場合に、光ファイバの実際の温度が350[K]であるとすると、約3[K]の測定誤差が生ずる。 Referring to FIG. 9, when the temperature of the optical fiber is the known temperature Tr in the above-described equation (3), the graph D101 and the graph D102 intersect. However, it can be seen that at a temperature deviating from the known temperature Tr , the temperature measurement error differs between the graph D101 and the graph D102 (depending on the value of the Raman shift wavenumber ν 1 ). For example, if the known temperature T r is 300 [K] and the actual temperature of the optical fiber is 350 [K], a measurement error of about 3 [K] occurs.

従来は、温度測定装置の性能がさほど高くなく、要求される測定精度もさほど高くなかったため、以上の測定誤差は許容される範囲であった。しかしながら、近年においては温度測定装置の性能が飛躍的に向上してストークス光及び反ストークス光を精度良く測定することができるようになったため、以上のラマンシフト波数の分布に起因する温度誤差が無視できなくなってきた。   Conventionally, the performance of the temperature measuring device is not so high, and the required measurement accuracy is not so high, so the above measurement error is in an allowable range. However, in recent years, the performance of temperature measuring devices has improved dramatically, and it has become possible to measure Stokes light and anti-Stokes light with high accuracy, so the temperature error due to the above-mentioned Raman shift wavenumber distribution is ignored. I can't do it.

本発明は上記事情に鑑みてなされたものであり、ラマンシフト波数の分布を考慮した高精度の温度測定を行うことができる温度測定装置を提供することを目的とする。   The present invention has been made in view of the above circumstances, and an object of the present invention is to provide a temperature measurement apparatus that can perform highly accurate temperature measurement in consideration of the distribution of Raman shift wavenumbers.

上記課題を解決するために、本発明の温度測定装置は、被測定対象(13)にレーザ光(P)を入射させて得られるラマン散乱光(P、P)を用いて前記被測定対象の温度を測定する温度測定装置(1)において、前記ラマン散乱光を受光して受光信号を出力する受光部(21a、21b)と、前記ラマン散乱光のラマンシフト波数の分布を考慮した補正テーブルを記憶する記憶部(26)と、前記ラマン散乱光が前記ラマンシフト波数の分布を有さないものとして前記受光部から出力される受光信号を用いて前記被測定対象の温度を求めるための所定の演算を行い、当該所定の演算によって得られる演算結果を前記記憶部に記憶された前記補正テーブルを用いて補正する演算処理部(25)とを備える。
この発明によると、被測定対象にレーザ光を入射させて得られるラマン散乱光が受光部で受光され、受光部から出力される受光信号を用いて被測定対象の温度を求めるための所定の演算が行われ、この所定の演算によって得られる演算結果が記憶部に記憶されたラマン散乱光のラマンシフト波数の分布を考慮した補正テーブルを用いて補正される。
また、本発明の温度測定装置は、前記補正テーブルが、前記ラマン散乱光に含まれるストークス光及び反ストークス光のラマンシフト波数の分布を考慮して求められた温度毎の前記ストークス光と前記反ストークス光との強度比を示すテーブルであり、前記演算処理部は、前記所定の演算として前記ストークス光と前記反ストークス光との強度比を求める演算を行い、当該演算により得られた強度比から求められる温度を前記補正テーブルを用いて補正することを特徴としている。
或いは、本発明の温度測定装置は、前記補正テーブルが、前記ラマン散乱光に含まれるストークス光及び反ストークス光が前記ラマンシフト波数の分布を有さないものとした場合に求められる温度と、前記ストークス光及び前記反ストークス光のラマンシフト波数の分布を考慮した場合に求められる温度との差分を温度毎に示すテーブルであり、前記演算処理部は、前記所定の演算として前記ストークス光と前記反ストークス光との強度比から前記被測定対象の温度を求める演算を行い、当該演算により得られた温度を前記補正テーブルを用いて補正することを特徴としている。
また、本発明の温度測定装置は、前記演算処理部が、前記記憶部に記憶される前記補正テーブルを、前記ラマン散乱光のラマンシフト波数の分布を示すデータを用いて作成することを特徴としている。
また、本発明の温度測定装置は、前記ラマン散乱光に含まれるストークス光の所定の成分及び反ストークス光の所定の成分をそれぞれ個別に通過させるフィルタ部を備えることを特徴としている。
また、本発明の温度測定装置は、前記被測定対象が光ファイバであり、前記受光部は、前記光ファイバにパルス状のレーザ光を入射させて得られる前記ラマン散乱光に含まれるストークス光及び反ストークス光を個別に受光し、前記演算処理部は、前記受光部から順次出力される受光信号を用いて前記所定の演算を順次行うことにより、前記光ファイバの長さ方向における温度分布を求めることを特徴としている。
In order to solve the above-described problem, the temperature measuring device of the present invention uses the Raman scattered light (P a , P s ) obtained by making the laser beam (P L ) incident on the measurement target (13). In the temperature measuring device (1) for measuring the temperature of the measurement object, the light receiving unit (21a, 21b) that receives the Raman scattered light and outputs a received light signal, and the distribution of the Raman shift wave number of the Raman scattered light are considered. In order to obtain the temperature of the object to be measured using a storage unit (26) for storing a correction table and a light reception signal output from the light receiving unit on the assumption that the Raman scattered light does not have a distribution of the Raman shift wave number. And a calculation processing unit (25) that corrects a calculation result obtained by the predetermined calculation using the correction table stored in the storage unit.
According to this invention, the Raman scattered light obtained by making the laser beam incident on the measurement target is received by the light receiving unit, and the predetermined calculation for obtaining the temperature of the measurement target using the received light signal output from the light receiving unit The calculation result obtained by this predetermined calculation is corrected using a correction table that takes into account the distribution of the Raman shift wave number of the Raman scattered light stored in the storage unit.
In the temperature measuring device according to the present invention, the correction table may include the Stokes light and the antireflection for each temperature obtained in consideration of a distribution of Raman shift wave numbers of Stokes light and anti-Stokes light included in the Raman scattered light. It is a table showing an intensity ratio with the Stokes light, and the arithmetic processing unit performs an operation for obtaining an intensity ratio between the Stokes light and the anti-Stokes light as the predetermined operation, and from the intensity ratio obtained by the operation The required temperature is corrected using the correction table.
Alternatively, in the temperature measuring device of the present invention, the correction table includes a temperature required when the Stokes light and the anti-Stokes light included in the Raman scattered light do not have a distribution of the Raman shift wavenumber, and the temperature It is a table showing the difference from the temperature obtained when considering the distribution of Raman shift wavenumbers of the Stokes light and the anti-Stokes light for each temperature, and the arithmetic processing unit performs the predetermined calculation as the Stokes light and the anti-Stokes light. A calculation for obtaining the temperature of the object to be measured from the intensity ratio with the Stokes light is performed, and the temperature obtained by the calculation is corrected using the correction table.
In the temperature measuring device of the present invention, the arithmetic processing unit creates the correction table stored in the storage unit using data indicating a distribution of Raman shift wavenumbers of the Raman scattered light. Yes.
In addition, the temperature measuring device of the present invention includes a filter unit that individually passes a predetermined component of Stokes light and a predetermined component of anti-Stokes light included in the Raman scattered light.
In the temperature measuring device of the present invention, the object to be measured is an optical fiber, and the light receiving unit includes Stokes light included in the Raman scattered light obtained by making pulsed laser light incident on the optical fiber, and The anti-Stokes light is individually received, and the calculation processing unit obtains a temperature distribution in the length direction of the optical fiber by sequentially performing the predetermined calculation using light reception signals sequentially output from the light receiving unit. It is characterized by that.

本発明によれば、被測定対象にレーザ光を入射させて得られるラマン散乱光を受光部で受光し、受光部から出力される受光信号を用いて被測定対象の温度を求めるための所定の演算を行い、この所定の演算によって得られる演算結果を記憶部に記憶されたラマン散乱光のラマンシフト波数の分布を考慮した補正テーブルを用いて補正しているため、ラマンシフト波数の分布を考慮した高精度の温度測定を行うことができるという効果がある。   According to the present invention, the Raman scattered light obtained by making the laser beam incident on the measurement target is received by the light receiving unit, and the predetermined temperature for obtaining the temperature of the measurement target using the light reception signal output from the light receiving unit. The calculation results obtained by this predetermined calculation are corrected using a correction table that takes into account the distribution of the Raman shift wave number of the Raman scattered light stored in the storage unit, so the distribution of the Raman shift wave number is taken into account. It is possible to perform temperature measurement with high accuracy.

本発明の第1実施形態による温度測定装置の要部構成を示すブロック図である。It is a block diagram which shows the principal part structure of the temperature measuring device by 1st Embodiment of this invention. 本発明の第1実施形態による温度測定装置の動作を示すフローチャートである。It is a flowchart which shows operation | movement of the temperature measurement apparatus by 1st Embodiment of this invention. 光ファイバ13内で生ずるストークス光及び反ストークス光のラマンシフト波数の分布を示す分布データの一例を示す図である。It is a figure which shows an example of distribution data which shows distribution of the Raman shift wave number of Stokes light and anti-Stokes light which arise in the optical fiber. 本発明の第1実施形態における初期化ステップS1で作成される補正テーブルの一例を示す図である。It is a figure which shows an example of the correction table produced by initialization step S1 in 1st Embodiment of this invention. 本発明の第2実施形態による温度測定装置の動作を示すフローチャートである。It is a flowchart which shows operation | movement of the temperature measurement apparatus by 2nd Embodiment of this invention. 本発明の第2実施形態における初期化ステップS3で作成される補正テーブルの一例を示す図である。It is a figure which shows an example of the correction table produced by initialization step S3 in 2nd Embodiment of this invention. 本発明の第3実施形態による温度測定装置が備える光フィルタの透過特性の一例を示す図である。It is a figure which shows an example of the transmission characteristic of the optical filter with which the temperature measuring device by 3rd Embodiment of this invention is provided. 光ファイバ内で実際に生ずるストークス光の波数分布を示す図の一例である。It is an example of the figure which shows the wave number distribution of the Stokes light actually produced within an optical fiber. 従来の温度測定装置における温度測定誤差の一例を示す図である。It is a figure which shows an example of the temperature measurement error in the conventional temperature measuring apparatus.

以下、図面を参照して本発明の実施形態による温度測定装置について詳細に説明する。   Hereinafter, a temperature measuring apparatus according to an embodiment of the present invention will be described in detail with reference to the drawings.

〔第1実施形態〕
図1は、本発明の第1実施形態による温度測定装置の要部構成を示すブロック図である。図1に示す通り、本実施形態の温度測定装置1は、パルス光源回路11、フィルタ部12、光ファイバ13(被測定対象)、処理部14、タイミング発生回路15、及び表示操作装置16を備える。この温度測定装置1は、所謂R−OTDR(Raman Optical Time Domain Reflectmetry)といわれるものであり、光ファイバ13内で生ずるラマン散乱光(ストークス光及び反ストークス光)を検出して光ファイバ13の長さ方向における温度分布を測定する。
[First Embodiment]
FIG. 1 is a block diagram showing a main configuration of the temperature measuring device according to the first embodiment of the present invention. As shown in FIG. 1, the temperature measurement device 1 of this embodiment includes a pulse light source circuit 11, a filter unit 12, an optical fiber 13 (measurement target), a processing unit 14, a timing generation circuit 15, and a display operation device 16. . This temperature measuring apparatus 1 is so-called R-OTDR (Raman Optical Time Domain Reflectmetry), and detects the Raman scattered light (Stokes light and anti-Stokes light) generated in the optical fiber 13 to detect the length of the optical fiber 13. Measure the temperature distribution in the vertical direction.

パルス光源回路11は、例えば半導体レーザ等の光源を備えており、タイミング発生回路15で規定されるタイミングでパルス状のレーザ光Pを射出する。尚、以下の説明では、パルス光源回路11から射出されるレーザ光Pの波数をνとする。フィルタ部12は、方向性結合器12aと光フィルタ12bとを備える。方向性結合器12aは、パルス光源回路11から射出されたレーザ光Pが光ファイバ13に導かれ、且つ、光ファイバ13で生じた後方散乱光Pが光フィルタ12に導かれるよう、パルス光源回路11、光ファイバ13、及び光フィルタ12bを光学的に結合する。 Pulse light source circuit 11 includes, for example, a light source such as a semiconductor laser emits a pulsed laser beam P L at a timing defined by the timing generation circuit 15. In the following description, the wave number of the laser light P L emitted from the pulse light source circuit 11 and [nu 0. The filter unit 12 includes a directional coupler 12a and an optical filter 12b. The directional coupler 12a, the laser beam P L emitted from the pulse light source circuit 11 is guided to the optical fiber 13, and so that backscattered light P B generated in the optical fiber 13 is guided to the optical filter 12, a pulse The light source circuit 11, the optical fiber 13, and the optical filter 12b are optically coupled.

光フィルタ12bは、方向性結合器12aからの後方散乱光Pに含まれるラマン散乱光(ストークス光P及び反ストークス光P)を抽出するとともに、ストークス光Pと反ストークス光Pとを分離して出力するフィルタである。尚、光ファイバ13で生ずるラマンシフト波数をνとすると、ストークス光Pの波数はν−νで表され、反ストークス光Pの波数はν+νで表される。光ファイバ13は、例えば数km〜数十km程度の長さを有する石英系マルチモード光ファイバを用いることができる。尚、シングルモード光ファイバを用いてもよい。 Optical filter 12b is configured to extract the Raman scattered light (Stokes light P s and the anti-Stokes light P a) contained in the backscattered light P B from the directional coupler 12a, the Stokes beam P s and the anti-Stokes light P a Is a filter that separates and outputs. Incidentally, when the Raman shift wave number occurring in the optical fiber 13 and [nu 1, wave number of the Stokes beam P s is represented by [nu 0 -v 1, wave number of the anti-Stokes light P a is represented by ν 0 + ν 1. As the optical fiber 13, for example, a silica-based multimode optical fiber having a length of about several kilometers to several tens of kilometers can be used. A single mode optical fiber may be used.

処理部14は、光電変換回路21a,21b(受光部)、増幅回路22a,22b、A/D変換回路23a,23b、平均化処理回路24a,24b、演算処理回路25(演算処理部)、及びメモリ26(記憶部)を備える。光電変換回路21a,21bは、例えばアバランシェ・フォトダイオード等の受光素子を備えており、フィルタ部12の光フィルタ12bから出力されるストークス光P及び反ストークス光Pを受光して光電変換信号(受光信号)をそれぞれ出力する。増幅回路22a,22bは光電変換回路21a,21bから出力される光電変換信号をそれぞれ所定の増幅率で増幅する。 The processing unit 14 includes photoelectric conversion circuits 21a and 21b (light receiving units), amplification circuits 22a and 22b, A / D conversion circuits 23a and 23b, averaging processing circuits 24a and 24b, an arithmetic processing circuit 25 (an arithmetic processing unit), and A memory 26 (storage unit) is provided. The photoelectric conversion circuit 21a, 21b has, for example, a light receiving element such as an avalanche photodiode, the Stokes beam P s and the photoelectric conversion signal by receiving anti-Stokes light P a is output from the optical filter 12b of the filter unit 12 (Light reception signal) is output. The amplifier circuits 22a and 22b amplify the photoelectric conversion signals output from the photoelectric conversion circuits 21a and 21b, respectively, with a predetermined amplification factor.

A/D変換回路23a,23bは、増幅回路22a,22bで増幅された光電変換信号をタイミング発生回路15で規定されるタイミングでサンプリングし、ディジタル化されたサンプルデータを出力する。ここで、図1に示す通り、タイミング発生回路15から出力されるタイミング信号TG2はA/D変換回路23a,23bの双方に入力されているため、A/D変換回路23a,23bは増幅回路22a,22bで増幅された光電変換信号を同一のタイミングでそれぞれサンプリングする。   The A / D conversion circuits 23a and 23b sample the photoelectric conversion signals amplified by the amplification circuits 22a and 22b at a timing defined by the timing generation circuit 15, and output digitized sample data. Here, as shown in FIG. 1, since the timing signal TG2 output from the timing generation circuit 15 is input to both of the A / D conversion circuits 23a and 23b, the A / D conversion circuits 23a and 23b are connected to the amplification circuit 22a. , 22b are sampled at the same timing.

平均化処理回路24aは、パルス光源回路11から射出されるレーザ光Pが複数回に亘って光ファイバ13に入射される度にA/D変換回路23aから出力されるサンプルデータを平均化する。同様に、平均化処理回路24bは、パルス光源回路11から射出されるレーザ光Pが複数回に亘って光ファイバ13に入射される度にA/D変換回路23bから出力されるサンプルデータを平均化する。光ファイバ13で生ずるラマン散乱光(ストークス光P及び反ストークス光P)は微弱であるため、光ファイバ13に対して複数回に亘ってレーザ光Pを入射させて得られるサンプルデータを平均化することにより、所望の信号対雑音比(S/N比)を得ている。 Averaging circuit 24a, the laser beam P L emitted from the pulse light source circuit 11 averages the sample data output from the A / D converter circuit 23a each time it is injected into the optical fiber 13 over a plurality of times . Similarly, averaging processing circuit 24b is a sample data laser beam P L emitted from the pulse light source circuit 11 is output from the A / D converter 23b each time it is injected into the optical fiber 13 over a plurality of times Average. Since Raman scattered light generated in the optical fiber 13 (Stokes light P s and the anti-Stokes light P a) is weak, the sample data obtained by a plurality of times with respect to the optical fiber 13 is incident laser light P L A desired signal-to-noise ratio (S / N ratio) is obtained by averaging.

演算処理回路25は、平均化処理回路24aで平均化処理が行われたストークス光Pに関するサンプルデータと、平均化処理回路24aで平均化処理が行われた反ストークス光Pに関するサンプルデータとを用いて、サンプルポイント(光ファイバ13内における位置と同義)毎の強度比を求め、この強度比からサンプルポイント毎の温度を求める演算を行う。 Arithmetic processing circuit 25, the sample data related to the Stokes beam P s that averaging processing is performed by the averaging processing circuit 24a, and examples of the anti-Stokes light P a of averaging process by the averaging processing circuit 24a is performed data Is used to calculate the intensity ratio for each sample point (synonymous with the position in the optical fiber 13) and calculate the temperature for each sample point from this intensity ratio.

ここで、光ファイバ13で生じた後方散乱光Pに含まれるストークス光P及び反ストークス光Pの波数はある広がりを有しており、ラマンシフト波数も定数(線スペクトル)ではなくある分布を有している。演算処理回路15は、上記のサンプルポイント毎の強度比を、ストークス光P及び反ストークス光Pが共にラマンシフト波数の分布を有さないものとして求める。 Here, the wave number of the Stokes beam P s and the anti-Stokes light P a is included in the backscattered light P B generated in the optical fiber 13 has a certain extent, it is not the Raman shift wavenumber be constant (line spectrum) Have a distribution. Arithmetic processing circuit 15, the intensity ratio for each sample point of the obtained as Stokes light P s and the anti-Stokes light P a is not both have a distribution of the Raman shift wavenumber.

また、演算処理回路25は、上記の演算によって求められるサンプルポイント毎の強度比を、メモリ26に記憶される補正テーブルを用いて補正する。この補正テーブルは、ストークス光P及び反ストークス光Pのラマンシフト波数の分布を考慮して求められた温度毎のストークス光と反ストークス光との強度比を示すテーブルであり、演算処理回路15において求められる強度比(ストークス光P及び反ストークス光Pが共にラマンシフト波数の分布を有さないものとして求められる強度比)を、ストークス光P及び反ストークス光Pのラマンシフト波数の分布を考慮した強度比に補正するものである。 In addition, the arithmetic processing circuit 25 corrects the intensity ratio for each sample point obtained by the above calculation using a correction table stored in the memory 26. This correction table is a table showing the intensity ratio of the Stokes beam P s and anti-Stokes light P Stokes light for each temperature distribution determined in consideration of the Raman shift wavenumber of a and anti-Stokes light, the arithmetic processing circuit intensity ratio obtained in 15 (intensity ratio Stokes light P s and the anti-Stokes light P a is both determined as having no distribution of the Raman shift wavenumber), Raman shift of the Stokes beam P s and the anti-Stokes light P a The intensity ratio is corrected in consideration of the wave number distribution.

補正テーブルを用いて補正された強度比からサンプルポイント毎の温度を求めることにより、ストークス光P及び反ストークス光Pのラマンシフト波数の分布を考慮した光ファイバ13の長さ方向における温度分布が高精度に求められる。更に、演算処理回路25は、ストークス光P及び反ストークス光Pのラマンシフト波数の分布を示すデータを用いてメモリ26に記憶される補正テーブルを作成することもできる。尚、補正テーブルを用いた補正方法及び補正テーブルの作成方法の詳細については後述する。 By determining the temperature of each sample point from the corrected intensity ratio using the correction table, the temperature distribution in the longitudinal direction of the optical fiber 13 in consideration of the distribution of the Raman shift wavenumber of the Stokes beam P s and the anti-Stokes light P a Is required with high accuracy. Furthermore, the arithmetic processing circuit 25 can also create a correction table stored in the memory 26 by using the data showing the distribution of Raman shift wavenumber of the Stokes beam P s and the anti-Stokes light P a. Details of the correction method using the correction table and the method of creating the correction table will be described later.

メモリ26は、例えば揮発性又は不揮発性のRAM(Random Access Memory)等であり、演算処理回路25で用いられる各種のデータを記憶する。例えば、ストークス光P及び反ストークス光Pのラマンシフト波数の分布を示すデータ、演算処理回路25によって作成される補正テーブル、演算処理回路25で演算を行う際に一時的に用いられるデータ等を記憶する。尚、このメモリ26は、演算処理回路25とは別に設けられていてもよく、演算処理回路25内に内蔵されていても良い。 The memory 26 is, for example, a volatile or nonvolatile RAM (Random Access Memory) or the like, and stores various data used in the arithmetic processing circuit 25. For example, data indicating the distribution of the Raman shift wavenumber of the Stokes beam P s and the anti-Stokes light P a, the correction table, data and the like which are temporarily used when performing arithmetic operation processing circuit 25 that is created by the processing circuit 25 Remember. The memory 26 may be provided separately from the arithmetic processing circuit 25 or may be built in the arithmetic processing circuit 25.

タイミング発生回路15は、パルス光源回路11に対して、パルス状のレーザ光Pを射出させるタイミングを規定するタイミング信号TG1を出力する。また、タイミング発生回路15は、A/D変換回路23a,23bに対して、これらが増幅回路22a,22bから出力される光電変換信号をそれぞれサンプリングするタイミングを規定するタイミング信号TG2を出力する。 Timing generating circuit 15, the pulse light source circuit 11, and outputs a timing signal TG1 that defines the timing to emit the pulsed laser light P L. Further, the timing generation circuit 15 outputs to the A / D conversion circuits 23a and 23b a timing signal TG2 that defines the timing at which the photoelectric conversion signals output from the amplification circuits 22a and 22b are sampled.

表示操作装置16は、例えば、液晶表示装置やCRT(Cathode Ray Tube)等の表示装置とユーザにより操作されるキーボードやマウス等の入力装置とを備えたコンピュータにより実現される。この表示操作装置16は、演算処理回路25から出力される演算結果(光ファイバ13の長さ方向における温度分布)を表示装置に表示するとともに、ユーザの操作に応じて、処理部14及びタイミング発生回路15等を制御する。   The display operation device 16 is realized by, for example, a computer including a display device such as a liquid crystal display device or a CRT (Cathode Ray Tube) and an input device such as a keyboard or a mouse operated by a user. The display operation device 16 displays the calculation result (temperature distribution in the length direction of the optical fiber 13) output from the calculation processing circuit 25 on the display device, and the processing unit 14 and timing generation according to a user operation. The circuit 15 and the like are controlled.

次に、上記構成による温度測定装置1の動作について説明する。図2は、本発明の第1実施形態による温度測定装置の動作を示すフローチャートである。図2に示す通り、温度測定装置1の動作は、補正テーブルを作成する初期化ステップS1と、光ファイバ13の長手方向における温度分布を測定する測定ステップS2とに大別される。以下、これらのステップで行われる処理の詳細について順に説明する。   Next, the operation of the temperature measuring apparatus 1 having the above configuration will be described. FIG. 2 is a flowchart showing the operation of the temperature measuring device according to the first embodiment of the present invention. As shown in FIG. 2, the operation of the temperature measuring device 1 is roughly divided into an initialization step S1 for creating a correction table and a measuring step S2 for measuring a temperature distribution in the longitudinal direction of the optical fiber 13. Hereinafter, details of processing performed in these steps will be described in order.

初期化ステップS1が開始されると、まずユーザの操作によって、光ファイバ13内で生ずるストークス光及び反ストークス光のラマンシフト波数の分布を示すデータ(分布データ)が温度測定装置1のメモリ26に記憶される(ステップS11)。この分布データは、予め光スペクトラムアナライザ等を用いて光ファイバ13を測定することにより得られるデータであり、光ファイバ13をなす材質(分子構造)によって決まるものであるため1度だけ測定すればよい。   When the initialization step S1 is started, data (distribution data) indicating the distribution of Raman shift wavenumbers of Stokes light and anti-Stokes light generated in the optical fiber 13 is first stored in the memory 26 of the temperature measurement device 1 by a user operation. Stored (step S11). This distribution data is data obtained by measuring the optical fiber 13 in advance using an optical spectrum analyzer or the like, and is determined by the material (molecular structure) forming the optical fiber 13, so that it needs to be measured only once. .

図3は、光ファイバ13内で生ずるストークス光及び反ストークス光のラマンシフト波数の分布を示す分布データの一例を示す図である。図3に示す通り、分布データは、ラマンシフト波数毎のストークス光及び反ストークス光の相対強度を示すデータである。尚、図3においては、ラマンシフト波数毎のストークス光の相対強度を示す分布データのみを図示しているが、ラマンシフト波数毎の反ストークス光の相対強度を示す分布データも同様のデータである。   FIG. 3 is a diagram illustrating an example of distribution data indicating the distribution of Raman shift wave numbers of Stokes light and anti-Stokes light generated in the optical fiber 13. As shown in FIG. 3, the distribution data is data indicating the relative intensity of Stokes light and anti-Stokes light for each Raman shift wavenumber. In FIG. 3, only the distribution data indicating the relative intensity of the Stokes light for each Raman shift wave number is illustrated, but the distribution data indicating the relative intensity of the anti-Stokes light for each Raman shift wave number is similar data. .

尚、上記の分布データは、ストークス光及び反ストークス光のラマンシフト波数の分布を正確に反映したデータであることが望ましいが、その分布の近似値を示すデータであっても良い。以下では、ストークス光のラマンシフト波数の分布をS(ν)とし、反ストークス光のラマンシフト波数の分布をA(ν)とする。   The distribution data is preferably data that accurately reflects the distribution of Raman shift wavenumbers of Stokes light and anti-Stokes light, but may be data indicating an approximate value of the distribution. Hereinafter, the Raman shift wave number distribution of Stokes light is S (ν), and the Raman shift wave number distribution of anti-Stokes light is A (ν).

次に、演算処理回路25によって、メモリ26に記憶された分布データが読み出され、ストークス光及び反ストークス光のラマンシフト波数の分布を考慮したストークス光と反ストークス光との強度比が算出される(ステップS12)。この処理では、まずメモリ26から読み出された分布データを正規化し、次いで正規化した分布データを考慮した既知の温度Tにおけるストークス光及び反ストークス光の強度をそれぞれ求め、最後にストークス光と反ストークス光との強度比を求める処理が順に行われる。 Next, the distribution data stored in the memory 26 is read out by the arithmetic processing circuit 25, and the intensity ratio between the Stokes light and the anti-Stokes light in consideration of the Raman shift wave number distribution of the Stokes light and the anti-Stokes light is calculated. (Step S12). In this process, the distribution data read from the memory 26 is first normalized, then the intensities of the Stokes light and the anti-Stokes light at a known temperature Tr taking into account the normalized distribution data are obtained, and finally the Stokes light and Processing for obtaining the intensity ratio with the anti-Stokes light is sequentially performed.

具体的には、以下の処理が順に行われる。分布S(ν),A(ν)を正規化した分布(正規化分布)をそれぞれS(ν),A(ν)とすると、メモリ26から読み出された分布データに対して例えば以下の(5)式に示す演算が行われて正規化分布A(ν),S(ν)が求められる。

Figure 2010223831
Specifically, the following processing is performed in order. Assuming that distributions (normalized distributions) obtained by normalizing the distributions S (ν) and A (ν) are respectively S * (ν) and A * (ν), for example, the distribution data read from the memory 26 is as follows. The normalized distributions A * (ν) and S * (ν) are obtained by performing the calculation shown in equation (5).
Figure 2010223831

次に、既知の温度Tにおけるストークス光及び反ストークス光の強度をラマンシフト波数毎に求める処理が行われ、以下の(6)式に示す通り、求められたストークス光及び反ストークス光の強度に対して上記(5)式に示す正規化分布S(ν),A(ν)をそれぞれ乗算して波数νで積分する処理が行われる。この処理によってラマンシフト波数の分布を考慮したストークス光及び反ストークス光の強度がそれぞれ求められる。

Figure 2010223831
Next, a process for obtaining the intensity of Stokes light and anti-Stokes light at a known temperature Tr for each Raman shift wavenumber is performed, and the intensity of the obtained Stokes light and anti-Stokes light is obtained as shown in the following equation (6). Is multiplied by the normalized distributions S * (ν) and A * (ν) shown in the above equation (5) and integrated at the wave number ν. By this processing, the intensities of Stokes light and anti-Stokes light in consideration of the Raman shift wave number distribution are obtained.
Figure 2010223831

最後に、以下の(7)式に示す通り、上記(6)式に示される反ストークス光の強度を示す式をストークス光の強度を示す式で除算する処理が行われる。これにより、ストークス光及び反ストークス光のラマンシフト波数の分布を考慮したストークス光と反ストークス光との強度比が算出される。

Figure 2010223831
Finally, as shown in the following equation (7), a process for dividing the equation indicating the anti-Stokes light intensity shown in the above equation (6) by the equation indicating the intensity of the Stokes light is performed. Thereby, the intensity ratio between the Stokes light and the anti-Stokes light in consideration of the Raman shift wave number distribution of the Stokes light and the anti-Stokes light is calculated.
Figure 2010223831

次いで、既知の温度Tにおけるストークス光及び反ストークス光の強度が従来の方法により測定され、その強度比が従来の方法で算出される(ステップS13)。つまり、ストークス光及び反ストークス光が共にラマンシフト波数の分布を有さないものとしてストークス光及び反ストークス光の強度が測定されるとともに強度比が算出される。ここで求められる強度比R(T)は、以下の(8)式で表される。

Figure 2010223831
Next, the intensities of Stokes light and anti-Stokes light at a known temperature Tr are measured by a conventional method, and the intensity ratio is calculated by a conventional method (step S13). That is, the intensity of the Stokes light and the anti-Stokes light is measured and the intensity ratio is calculated on the assumption that both the Stokes light and the anti-Stokes light do not have a Raman shift wave number distribution. The intensity ratio R (T r ) obtained here is expressed by the following equation (8).
Figure 2010223831

そして、以下の(9)式に示す通り、ステップS12で得られた強度比をステップS13で得られた強度比で除算する処理を行う。かかる処理は、メモリ26に記憶される分布データを得た光スペクトラムアナライザ等と温度測定装置1との測定感度の違いを解消するために行われる処理である。

Figure 2010223831
Then, as shown in the following equation (9), a process of dividing the intensity ratio obtained in step S12 by the intensity ratio obtained in step S13 is performed. Such processing is performed in order to eliminate the difference in measurement sensitivity between the optical spectrum analyzer or the like that obtained the distribution data stored in the memory 26 and the temperature measurement device 1.
Figure 2010223831

以上の処理によって求められた強度比R(T)は、ある既知の温度Tにおけるものである。このため、測定すべき温度範囲内において温度Tを変化させながら以上の処理を繰り返すことによって、温度補正テーブルが作成される(ステップS14)。図4は、本発明の第1実施形態における初期化ステップS1で作成される補正テーブルの一例を示す図である。図4に示す通り、補正テーブルは、ストークス光及び反ストークス光のラマンシフト波数の分布を考慮して求められた温度毎のストークス光と反ストークス光との強度比を示すテーブルである。 The intensity ratio R (T r ) obtained by the above processing is at a known temperature T r . Therefore, a temperature correction table is created by repeating the above processing while changing the temperature Tr within the temperature range to be measured (step S14). FIG. 4 is a diagram showing an example of the correction table created in the initialization step S1 in the first embodiment of the present invention. As shown in FIG. 4, the correction table is a table showing the intensity ratio of Stokes light and anti-Stokes light at each temperature obtained in consideration of the Raman shift wave number distribution of Stokes light and anti-Stokes light.

尚、上記のステップS13では、ストークス光と反ストークス光の強度比を従来の方法で求めるために、既知の温度Tにおけるストークス光及び反ストークス光の強度を従来の方法で測定している。しかしながら、既知の温度Tを変えながら従来の方法による測定を予め行ってストークス光及び反ストークス光の強度をそれぞれ求めておけば、ステップS13で測定を行うことなく演算処理回路25の演算のみで補正テーブルを作成することも可能である。尚、作成された補正テーブルは例えばメモリ26に記憶される。 In step S13, the intensity of Stokes light and anti-Stokes light at a known temperature Tr is measured by a conventional method in order to obtain the intensity ratio of Stokes light and anti-Stokes light by a conventional method. However, if the intensity of the Stokes light and the anti-Stokes light is obtained in advance by performing the measurement by the conventional method while changing the known temperature Tr , only the calculation of the arithmetic processing circuit 25 is performed without performing the measurement in step S13. It is also possible to create a correction table. The created correction table is stored in the memory 26, for example.

以上の初期化ステップS1が終了すると測定ステップS2が行われる。測定ステップS2が開始されると、まず温度が不明である光ファイバ13の測定が行われる(ステップS21)。光ファイバ13の測定が開始されると、図1に示すタイミング発生回路15からタイミング信号TG1が出力され、このタイミング信号TG1に基づいてパルス光源回路11からパルス状のレーザ光Pが射出される。このレーザ光Pは、光フィルタ部12を介して光ファイバ13に入射し、光ファイバ13内を伝播する。これにより、光ファイバ13内でラマン散乱光(ストークス光P及び反ストークス光P)を含む後方散乱光Pが発生する。 When the initialization step S1 is completed, a measurement step S2 is performed. When the measurement step S2 is started, first, the optical fiber 13 whose temperature is unknown is measured (step S21). When the measurement of the optical fiber 13 is started, the output timing signal TG1 from the timing generating circuit 15 shown in FIG. 1, the pulsed laser light P L is emitted from the pulse light source circuit 11 based on the timing signal TG1 . The laser beam P L enters the optical fiber 13 via the optical filter portion 12, propagates through the optical fiber 13. As a result, backscattered light P B including Raman scattered light (Stokes light P s and anti-Stokes light P a ) is generated in the optical fiber 13.

この後方散乱光Pは、光ファイバ13内をレーザ光Pの進行方向とは逆方向に進み、フィルタ部12の光フィルタ12bに入射する。そして、ストークス光Pと反ストークス光Pとが抽出されて分離される。ストークス光P及び反ストークス光Pは、光電変換回路21a,21bでそれぞれ光電変換されて、それらの光電変換信号が増幅回路22a,22bでそれぞれ増幅される。増幅回路22a,22bで増幅された光電変換信号は、A/D変換回路23a,23bにおいて、タイミング発生回路15から出力されるタイミング信号TG2で規定されるタイミングでそれぞれサンプリングされる。 The backscattered light P B is the traveling direction of the laser beam P L through the optical fiber 13 proceeds in the reverse direction, is incident to the optical filter 12b of the filter unit 12. Then, the Stokes beam P s and the anti-Stokes light P a is separated is extracted. Stokes light P s and the anti-Stokes light P a photoelectric conversion circuit 21a, are respectively converted photoelectrically 21b, their photoelectric conversion signals are amplified respectively in the amplifier circuits 22a, 22b. The photoelectric conversion signals amplified by the amplifier circuits 22a and 22b are sampled at the timings defined by the timing signal TG2 output from the timing generation circuit 15 in the A / D conversion circuits 23a and 23b, respectively.

光ファイバ13にパルス状のレーザ光Pが入射される度に、以上の処理が繰り返し行われ、A/D変換回路23aから出力されるストークス光Pに関するサンプルデータが平均化処理回路24aでサンプルポイント毎に平均化されるとともに、A/D変換回路23bから出力される反ストークス光Pに関する複数のサンプルデータが平均化処理回路24bでサンプルポイント毎に平均化される。 Each time the optical fiber 13 to the pulsed laser light P L is incident, the above process is repeated, A / D conversion circuit 23a Stokes light P s for the sample data averaging processing circuit 24a which is output from the Averaging is performed for each sample point, and a plurality of sample data relating to the anti-Stokes light Pa output from the A / D conversion circuit 23b is averaged for each sample point by the averaging processing circuit 24b.

A/D変換回路23a,23bの各々でサンプルポイント毎に平均化されたサンプルデータは演算処理回路25に入力され、ストークス光Pに関するサンプルデータと反ストークス光Pに関するサンプルデータとの比がサンプルポイント毎に求められる(ステップS22)。つまり、ストークス光及び反ストークス光が共にラマンシフト波数の分布を有さないものとして、サンプルポイント毎のストークス光と反ストークス光との強度比が求められる。 A / D conversion circuit 23a, a sample data averaged for each sample point in each of 23b are input to the arithmetic processing circuit 25, the ratio of the sample data and the anti-Stokes light P a for the sample data related to the Stokes beam P s It is obtained for each sample point (step S22). That is, the intensity ratio between the Stokes light and the anti-Stokes light for each sample point is obtained on the assumption that both the Stokes light and the anti-Stokes light do not have a Raman shift wave number distribution.

次に、ステップS22で求められた強度比を、初期化ステップS1で作成された補正テーブルを用いて補正して温度を算出する処理が行われる(ステップS23)。例えば、あるサンプルポイントの強度比が「1.00000」であった場合には、図4に示す補正テーブルを用いて、そのサンプルポイントにおける温度が300[K]であると求められる。   Next, a process of calculating the temperature by correcting the intensity ratio obtained in step S22 using the correction table created in initialization step S1 is performed (step S23). For example, when the intensity ratio of a certain sample point is “1.000000”, the temperature at that sample point is determined to be 300 [K] using the correction table shown in FIG.

ここで、前述した(3)式中の強度比R(T)に上記の強度比「1.00000」を代入すれば従来の方法による温度が求められる。本実施形態では、図4に示す補正テーブルを用いて温度を求めることにより、(3)式中から求められる温度(ストークス光P及び反ストークス光Pが共にラマンシフト波数の分布を有さないものとして求められる温度)を、ストークス光P及び反ストークス光Pのラマンシフト波数の分布を考慮して求めた温度に補正している。 Here, if the above-described intensity ratio “1.00000” is substituted into the intensity ratio R (T) in the above-described equation (3), the temperature by the conventional method can be obtained. In the present embodiment, by determining the temperature using the correction table shown in FIG. 4, (3) have a both distribution of Raman shift wavenumber sought temperature (Stokes light P s and the anti-Stokes light P a from formula the temperature) determined as not, is corrected to a temperature obtained by considering the distribution of Raman shift wavenumber of the Stokes beam P s and the anti-Stokes light P a.

次いで、演算処理回路25によって、全てのサンプルポイントの温度を算出したか否かが判断される(ステップS24)。温度の算出を行っていないサンプルポイントが存在すると判断した場合(判断結果が「NO」の場合)には、ステップS23の処理が繰り返され、全てのサンプルポイントについて温度の算出を行ったと判断した場合(判断結果が「YES」の場合)には、一連の測定ステップS2が終了する。このようにして、ストークス光P及び反ストークス光Pのラマンシフト波数の分布を考慮した光ファイバ13の長手方向における温度分布が高精度に求められる。 Next, the arithmetic processing circuit 25 determines whether or not the temperatures of all the sample points have been calculated (step S24). When it is determined that there is a sample point for which the temperature is not calculated (when the determination result is “NO”), the process of step S23 is repeated, and it is determined that the temperature has been calculated for all the sample points. When the determination result is “YES”, the series of measurement steps S2 is completed. In this way, the temperature distribution in the longitudinal direction of the optical fiber 13 in consideration of the distribution of the Raman shift wavenumber of the Stokes beam P s and the anti-Stokes light P a is determined with high accuracy.

尚、図4に示す補正テーブルの温度刻みは5[K]であるが、より細かい温度刻みの補正テーブルを作成すれば温度分布の測定精度を向上させることができる。また、図4に示す補正テーブルを用いて温度を測定する場合において、ステップS22で算出された強度比に応じて補正テーブルを補完することにより、より高い精度で温度分布を測定することができる。また、既知の温度Tが一定ならば、一度補正テーブルを作成してしまえば繰り返し同じ補正テーブルを用いることができる。このため、温度分布の測定を行った光ファイバ13を再び測定する場合には、図2に示す初期化処理S1を省略することができる。更に、既知の温度Tが変動しないように既知の温度Tを求める光ファイバ13の区間を恒温槽内に配置する等の対策を行っても良い。 Although the temperature increment of the correction table shown in FIG. 4 is 5 [K], the measurement accuracy of the temperature distribution can be improved by creating a finer temperature increment correction table. In the case where the temperature is measured using the correction table shown in FIG. 4, the temperature distribution can be measured with higher accuracy by complementing the correction table according to the intensity ratio calculated in step S22. If the known temperature Tr is constant, once the correction table is created, the same correction table can be used repeatedly. Therefore, when the optical fiber 13 for which the temperature distribution has been measured is measured again, the initialization process S1 shown in FIG. 2 can be omitted. Furthermore, measures may be performed such as to place the section of optical fiber 13 to obtain the known temperature T r as known temperature T r does not vary in a thermostatic chamber.

〔第2実施形態〕
次に、本発明の第2実施形態による温度測定装置について詳細に説明する。本実施形態の温度測定装置の構成は以上説明した第1実施形態による温度測定装置の構成とほぼ同じである。但し、演算処理回路25における補正テーブルの作成方法及び補正テーブルを用いた温度分布の補正方法が異なる。
[Second Embodiment]
Next, the temperature measuring device according to the second embodiment of the present invention will be described in detail. The configuration of the temperature measurement device of this embodiment is substantially the same as the configuration of the temperature measurement device according to the first embodiment described above. However, the correction table creation method and the temperature distribution correction method using the correction table in the arithmetic processing circuit 25 are different.

具体的に、演算処理回路25は、光ファイバ13内で生ずるストークス光及び反ストークス光が共にラマンシフト波数の分布を有さないものとした場合に求められる温度と、ストークス光及び反ストークス光のラマンシフト波数の分布を考慮して求められる温度との差分を温度毎に示す補正テーブルを作成する。また、演算処理回路25は、光ファイバ13内で生ずるストークス光及び反ストークス光が共にラマンシフト波数の分布を有さないものとして従来の方法により各サンプリングポイントの温度を求め、この温度の各々を補正テーブルを用いて補正する処理を行う。   Specifically, the arithmetic processing circuit 25 calculates the temperature required when both Stokes light and anti-Stokes light generated in the optical fiber 13 do not have a Raman shift wave number distribution, and the Stokes light and anti-Stokes light. A correction table is created that shows the difference from the temperature determined in consideration of the Raman shift wave number distribution for each temperature. Further, the arithmetic processing circuit 25 obtains the temperature of each sampling point by a conventional method assuming that both the Stokes light and the anti-Stokes light generated in the optical fiber 13 do not have a Raman shift wave number distribution, and calculates each of these temperatures. A correction process is performed using the correction table.

以下、本実施形態の温度測定装置の動作について説明する。図5は、本発明の第2実施形態による温度測定装置の動作を示すフローチャートである。図5に示す通り、温度測定装置の動作は、補正テーブルを作成する初期化ステップS3と、光ファイバ13の長手方向における温度分布を測定する測定ステップS3とに大別される。以下、これらのステップで行われる処理の詳細について順に説明する。   Hereinafter, the operation of the temperature measuring apparatus of the present embodiment will be described. FIG. 5 is a flowchart showing the operation of the temperature measuring apparatus according to the second embodiment of the present invention. As shown in FIG. 5, the operation of the temperature measurement device is roughly divided into an initialization step S3 for creating a correction table and a measurement step S3 for measuring a temperature distribution in the longitudinal direction of the optical fiber 13. Hereinafter, details of processing performed in these steps will be described in order.

初期化ステップS3が開始されると、まずストークス光及び反ストークス光の強度が従来の方法により測定され、その強度比が従来の方法で算出される(ステップS31)。つまり、ストークス光及び反ストークス光が共にラマンシフト波数の分布を有さないものとしてストークス光及び反ストークス光の強度が測定されて強度比が算出される。次いで、ユーザの操作によって、光ファイバ13内で生ずるストークス光及び反ストークス光のラマンシフト波数の分布を示す分布データが温度測定装置のメモリ26に記憶される(ステップS32)。   When the initialization step S3 is started, first, the intensities of Stokes light and anti-Stokes light are measured by the conventional method, and the intensity ratio is calculated by the conventional method (step S31). That is, the intensity ratio is calculated by measuring the intensity of the Stokes light and the anti-Stokes light, assuming that both the Stokes light and the anti-Stokes light do not have a Raman shift wave number distribution. Next, distribution data indicating the distribution of Raman shift wavenumbers of Stokes light and anti-Stokes light generated in the optical fiber 13 is stored in the memory 26 of the temperature measurement device (step S32).

次に、演算処理回路25によって、メモリ26に記憶された分布データが読み出され、この分布データを乗じたストークス光及び反ストークス光の温度毎の強度を示すテーブルが作成される(ステップS33)。具体的には、(1)式に示す式に対して温度Tをパラメータとして指定して温度毎の反ストークス光の強度A及びストークス光の強度Sをそれぞれ求め、これら温度毎の反ストークス光の強度A及びストークス光の強度Sに対して反ストークス光及びストークス光の分布データをそれぞれ乗算する処理が行われる。尚、ラマンシフト波数νとしては、実際の光ファイバ13の有するラマンシフト波数の分布に含まれる所定の波数を用いればよい。 Next, the distribution data stored in the memory 26 is read out by the arithmetic processing circuit 25, and a table indicating the intensity for each temperature of the Stokes light and the anti-Stokes light multiplied by this distribution data is created (step S33). . Specifically, (1) respectively obtained the anti-Stokes light intensity A s and Stokes light intensities S t for each temperature by specifying a temperature T as parameters for the expressions expression, anti-Stokes These per temperature processing of multiplying each distribution data of the anti-Stokes light and the Stokes light to the intensity of the light a s and Stokes light intensities S t is performed. As the Raman shift wave number ν, a predetermined wave number included in the distribution of Raman shift wave numbers of the actual optical fiber 13 may be used.

以上の処理が終了すると、演算処理回路25によって、ストークス光及び反ストークス光が共にラマンシフト波数の分布を有さないものとしてストークス光と反ストークス光との温度毎の強度比を示すテーブルが作成される(ステップS34)。具体的には、既知の温度Tにおける強度比がステップS31で算出された強度比と等しくなるように、ステップS33で作成されたテーブルを補正(線形補正)する処理が行われる。 When the above processing is completed, a table indicating the intensity ratio of Stokes light and anti-Stokes light for each temperature is created by the arithmetic processing circuit 25 assuming that both the Stokes light and anti-Stokes light do not have a Raman shift wave number distribution. (Step S34). Specifically, a process of correcting (linear correction) the table created in step S33 is performed so that the intensity ratio at the known temperature Tr is equal to the intensity ratio calculated in step S31.

ここで、ステップS31で測定された既知の温度Tにおけるストークス光及び反ストークス光の強度をそれぞれA(T),S(T)とすると、既知の温度Tにおける反ストークス光とストークス光との強度比Pは、以下の(10)式で表される。

Figure 2010223831
Here, assuming that the intensities of the Stokes light and the anti-Stokes light at the known temperature T r measured in step S31 are A s (T r ) and S t (T r ), respectively, the anti-Stokes light at the known temperature T r . the intensity ratio P 1 of the Stokes light is expressed by the following equation (10).
Figure 2010223831

よって、線形補正されたテーブルは、温度Tをパラメータにした以下の(11)式に示す強度比Pで表される。

Figure 2010223831
Thus, linear correction tables is represented by the intensity ratio P 2 shown in the following (11) in which the temperature T as a parameter.
Figure 2010223831

次に、演算処理回路25によって、ラマンシフト波数分布を考慮したストークス光と反ストークス光との温度毎の強度比を示すテーブルが作成される(ステップS35)。この処理では、まず正規化した分布データを考慮したストークス光及び反ストークス光の温度毎の強度を示すテーブルを作成し、次にストークス光と反ストークス光との温度毎の強度比を示すテーブルを作成し、最後に既知の温度Tにおける強度比がステップS31で算出された強度比と等しくなるようにテーブルを補正(線形補正)する処理が順に行われる。 Next, the arithmetic processing circuit 25 creates a table indicating the intensity ratio of Stokes light and anti-Stokes light for each temperature in consideration of the Raman shift wave number distribution (step S35). In this process, first, a table indicating the intensity for each temperature of the Stokes light and the anti-Stokes light considering the normalized distribution data is created, and then a table indicating the intensity ratio for each temperature of the Stokes light and the anti-Stokes light. Finally, a process of correcting (linear correction) the table so that the intensity ratio at the known temperature Tr is equal to the intensity ratio calculated in step S31 is sequentially performed.

具体的には、以下に示す処理が順に行われる。まず、以下の(12)式に示す通り、ストークス光及び反ストークス光の強度に対して正規化分布S(ν),A(ν)をそれぞれ乗算して波数νで積分し、温度Tをパラメータとしてストークス光及び反ストークス光の強度を示すテーブルを作成する処理が行われる。尚、以下の(12)式中におけるA,Sは、温度測定装置の性能及び光ファイバの特性に依存する固有の定数である。

Figure 2010223831
Specifically, the following processing is performed in order. First, as shown in the following equation (12), the intensity of Stokes light and anti-Stokes light is multiplied by normalized distributions S * (ν) and A * (ν), respectively, and integrated by wave number ν, and temperature T Is used as a parameter to create a table indicating the intensity of Stokes light and anti-Stokes light. In the following equation (12), A 1 and S 1 are inherent constants depending on the performance of the temperature measuring device and the characteristics of the optical fiber.
Figure 2010223831

次に、以下の(13)式に示す通り、上記(12)式に示される反ストークス光の強度を示す式をストークス光の強度を示す式で除算する処理が行われる。これにより、ストークス光及び反ストークス光のラマンシフト波数の分布を考慮したストークス光と反ストークス光との強度比が算出される。

Figure 2010223831
Next, as shown in the following expression (13), a process for dividing the expression indicating the intensity of anti-Stokes light shown in the above expression (12) by the expression indicating the intensity of Stokes light is performed. Thereby, the intensity ratio between the Stokes light and the anti-Stokes light in consideration of the Raman shift wave number distribution of the Stokes light and the anti-Stokes light is calculated.
Figure 2010223831

最後に、既知の温度Tにおける強度比がステップS31で算出された強度比と等しくなるようにテーブルを補正(線形補正)する処理が順に行われる。ここで、ステップS31で測定された既知の温度Tにおけるストークス光及び反ストークス光の強度をそれぞれA(T),S(T)とすると、既知の温度Tにおける反ストークス光とストークス光との強度比Pは、以下の(14)式で表される。

Figure 2010223831
Finally, the process of correcting the table (linear correction) is sequentially performed so that the intensity ratio at the known temperature Tr is equal to the intensity ratio calculated in step S31. Here, assuming that the intensities of the Stokes light and the anti-Stokes light at the known temperature T r measured in step S31 are A s (T r ) and S t (T r ), respectively, the anti-Stokes light at the known temperature T r . the intensity ratio P 3 of the Stokes light is expressed by the following equation (14).
Figure 2010223831

よって、線形補正されたテーブルは、温度Tをパラメータにした以下の(15)式に示す強度比Pで表される。

Figure 2010223831
Thus, linear correction tables is represented by the intensity ratio P 4 shown in the following equation (15) in which the temperature T as a parameter.
Figure 2010223831

そして、演算処理回路25によって、ステップ34で作成されたテーブルとステップS35で作成されたテーブルとの温度毎の差分を算出する処理が行われる。以上の処理によって、補正テーブルが作成される。図6は、本発明の第2実施形態における初期化ステップS3で作成される補正テーブルの一例を示す図である。図6に示す通り、補正テーブルは、光ファイバ13内で生ずるストークス光及び反ストークス光が共にラマンシフト波数の分布を有さないものとした場合に求められる温度と、ストークス光及び反ストークス光のラマンシフト波数の分布を考慮して求められる温度との差分を温度毎に示すテーブルである。尚、作成された補正テーブルは例えばメモリ26に記憶される。   Then, the arithmetic processing circuit 25 performs a process of calculating a difference for each temperature between the table created in step 34 and the table created in step S35. The correction table is created by the above processing. FIG. 6 is a diagram showing an example of the correction table created in the initialization step S3 in the second embodiment of the present invention. As shown in FIG. 6, the correction table includes the temperature required when both the Stokes light and the anti-Stokes light generated in the optical fiber 13 do not have a Raman shift wave number distribution, and the Stokes light and the anti-Stokes light. It is a table which shows the difference with temperature calculated | required in consideration of distribution of a Raman shift wave number for every temperature. The created correction table is stored in the memory 26, for example.

以上の初期化ステップS3が終了すると測定ステップS4が行われる。測定ステップS4が開始されると、まず温度が不明である光ファイバ13の測定が行われる(ステップS41)。尚、光ファイバ13の測定時における温度測定装置の動作は第1実施形態と同様であるため、ここでは説明を省略する。   When the initialization step S3 is completed, a measurement step S4 is performed. When the measurement step S4 is started, first, the optical fiber 13 whose temperature is unknown is measured (step S41). The operation of the temperature measuring device at the time of measuring the optical fiber 13 is the same as that of the first embodiment, and thus the description thereof is omitted here.

ステップS41の測定を行って得られたサンプルポイント毎のサンプルデータ(A/D変換回路23a,23bの各々でサンプルポイント毎に平均化されたサンプルデータ)は演算処理回路25に入力され、A/D変換回路23aから出力されるサンプルデータとA/D変換回路23bから出力されるサンプルデータとの比がサンプルポイント毎に求められる。そして、この比からサンプルポイント毎の温度が算出される(ステップS42)。つまり、ストークス光及び反ストークス光が共にラマンシフト波数の分布を有さないものとして、サンプルポイント毎のストークス光と反ストークス光との強度比が求められて温度が算出される。   Sample data for each sample point (sample data averaged for each sample point in each of the A / D conversion circuits 23a and 23b) obtained by performing the measurement in step S41 is input to the arithmetic processing circuit 25, and A / A ratio of sample data output from the D conversion circuit 23a and sample data output from the A / D conversion circuit 23b is obtained for each sample point. And the temperature for every sample point is calculated from this ratio (step S42). That is, assuming that both the Stokes light and the anti-Stokes light do not have a Raman shift wave number distribution, the intensity ratio between the Stokes light and the anti-Stokes light for each sample point is obtained and the temperature is calculated.

次に、ステップS42で求められた温度を、初期化ステップS3で作成された補正テーブルを用いて補正してする処理が行われる(ステップS43)。例えば、あるサンプルポイントの温度が320[K]であった場合には、図6に示す補正テーブルを用いて、そのサンプルポイントにおける温度が320[K]から−0.68[K]だけ減じられた319.32[K]に補正される。   Next, a process of correcting the temperature obtained in step S42 using the correction table created in initialization step S3 is performed (step S43). For example, when the temperature of a certain sample point is 320 [K], the temperature at that sample point is reduced by −0.68 [K] from 320 [K] using the correction table shown in FIG. To 319.32 [K].

次いで、演算処理回路25によって、全てのサンプルポイントの温度を算出したか否かが判断される(ステップS44)。温度の算出を行っていないサンプルポイントが存在すると判断した場合(判断結果が「NO」の場合)には、ステップS43の処理が繰り返され、全てのサンプルポイントについて温度の算出を行ったと判断した場合(判断結果が「YES」の場合)には、一連の測定ステップS4が終了する。このようにして、ストークス光P及び反ストークス光Pのラマンシフト波数の分布を考慮した光ファイバ13の長手方向に温度分布が高精度に求められる。 Next, the arithmetic processing circuit 25 determines whether or not the temperatures of all sample points have been calculated (step S44). When it is determined that there is a sample point for which the temperature is not calculated (when the determination result is “NO”), the process of step S43 is repeated, and it is determined that the temperature has been calculated for all the sample points. When the determination result is “YES”, a series of measurement steps S4 is completed. In this way, the longitudinal temperature distribution of the Stokes beam P s and the anti-Stokes light P a optical fiber 13 in consideration of the distribution of the Raman shift wavenumber of is required with high accuracy.

尚、本実施形態においても、第1実施形態と同様に、より高い精度で温度分布を測定するために、より細かい温度刻みの補正テーブルを作成しても良く、ステップS42で算出された温度に応じて補正テーブルを補完しても良い。また、第1実施形態と同様に、既知の温度Tが変動しないように既知の温度Tを求める光ファイバ13の区間を恒温槽内に配置する等の対策を行っても良い。 In the present embodiment, as in the first embodiment, in order to measure the temperature distribution with higher accuracy, a correction table with finer temperature increments may be created, and the temperature calculated in step S42 may be set. The correction table may be complemented accordingly. As in the first embodiment, measures may be carried out, such as known temperature T r to place a section of the optical fiber 13 to obtain the known temperature T r so as not to vary in a thermostatic chamber.

〔第3実施形態〕
次に、本発明の第3実施形態による温度測定装置について詳細に説明する。前述した第1,第2実施形態では、ストークス光及び反ストークス光のラマンシフト波数の分布を考慮した補正テーブルを作成し、ストークス光及び反ストークス光が共にラマンシフト波数の分布を有さないものとして求められた温度分布を補正テーブルを用いて補正することにより、光ファイバ13の長さ方向における温度分布を高精度に測定するものであった。
[Third Embodiment]
Next, a temperature measuring device according to a third embodiment of the present invention will be described in detail. In the first and second embodiments described above, a correction table that takes into account the distribution of Raman shift wave numbers of Stokes light and anti-Stokes light is created, and both the Stokes light and anti-Stokes light do not have a distribution of Raman shift wave numbers. The temperature distribution in the length direction of the optical fiber 13 is measured with high accuracy by correcting the temperature distribution obtained as follows using a correction table.

これに対し、本実施形態の温度測定装置は、ラマンシフト波数の分布を有するストークス光P及び反ストークス光Pを個別に線スペクトル状にすることで、ラマンシフト波数の分布によって生ずる温度誤差の低減を図り、もって温度分布の測定精度を高めるものである。本実施形態の温度測定装置の構成は前述した第1,第2実施形態による温度測定装置の構成とほぼ同じである。但し、光フィルタ12bが、ストークス光P及び反ストークス光Pを線スペクトル状にする透過特性を有する点が相違する。 In contrast, the temperature measuring device of the present embodiment, by the Stokes beam P s and the anti-Stokes light P a have a distribution of Raman shift wavenumber individually line spectrum shape, the temperature error caused by the distribution of Raman shift wavenumber Therefore, the measurement accuracy of the temperature distribution is improved. The configuration of the temperature measurement device of this embodiment is substantially the same as the configuration of the temperature measurement device according to the first and second embodiments described above. However, the optical filter 12b is different in that it has a transmission characteristic that makes the Stokes light P s and the anti-Stokes light Pa a linear spectrum.

図7は、本発明の第3実施形態による温度測定装置が備える光フィルタの透過特性の一例を示す図であって、(a)はストークス光に係る光フィルタの透過特性を示す図であり、(b)は反ストークス光に係る光フィルタの透過特性を示す図である。尚、図7中に示す波数νはパルス光源回路11から射出されるレーザ光P(光ファイバ13に入射するレーザ光)の波数である FIG. 7 is a diagram illustrating an example of the transmission characteristics of an optical filter provided in the temperature measurement device according to the third embodiment of the present invention, and FIG. 7A is a diagram illustrating the transmission characteristics of an optical filter related to Stokes light. (B) is a figure which shows the permeation | transmission characteristic of the optical filter based on anti-Stokes light. The wave number ν 0 shown in FIG. 7 is the wave number of the laser beam P L (laser beam incident on the optical fiber 13) emitted from the pulse light source circuit 11.

図7(a)に示す通り、ストークス光に係る光フィルタは、ラマンシフト波数の分布を有するストークス光Pのうちの波数ν11〜ν12を有する波数成分のみを透過させるバンドパスフィルタの透過特性C1を有する。また、図7(b)に示す通り、反ストークス光に係る光フィルタは、ラマンシフト波数の分布を有する反ストークス光Pのうちの波数ν21〜ν22を有する波数成分のみを透過させるバンドパスフィルタの透過特性C2を有する。 As shown in FIG. 7A, the optical filter related to Stokes light transmits a band-pass filter that transmits only wave number components having wave numbers ν 11 to ν 12 of the Stokes light P s having a Raman shift wave number distribution. It has characteristic C1. Further, as shown in FIG. 7 (b), the optical filter according to the anti-Stokes light, the band which transmits only the wave number component having a wave number ν 2122 of the anti-Stokes light P a having a distribution of Raman shift wavenumber It has a transmission characteristic C2 of a pass filter.

尚、これらのフィルタの透過特性は、複数の波数成分を透過させる透過特性である必要は必ずしも無く、反ストークス光及びストークス光の各々について少なくとも1つの波数成分を透過させる透過特性であればよい。また、これらのフィルタは、波数成分の透過帯域が可変であるのが望ましい。ここで、反ストークス光の強度は、ラマンシフト波数と温度とによって変動するため、必ずしも強度が最大となる波数が最大の温度依存性を示すとは限らない。透過帯域が可変のフィルタを用いることで、反ストークスの温度依存性が最大となる波数を選択するといったことが可能となる。   Note that the transmission characteristics of these filters are not necessarily transmission characteristics that transmit a plurality of wave number components, but may be transmission characteristics that transmit at least one wave number component for each of anti-Stokes light and Stokes light. Further, it is desirable that these filters have a variable transmission band of wave number components. Here, since the intensity of the anti-Stokes light varies depending on the Raman shift wave number and the temperature, the wave number with the maximum intensity does not necessarily exhibit the maximum temperature dependence. By using a filter with a variable transmission band, it is possible to select a wave number that maximizes the temperature dependence of anti-Stokes.

ストークス光及び反ストークス光を線スペクトル状にする透過特性を有する光フィルタを光フィルタ12bとして用いれば、ストークス光及び反ストークス光が共にラマンシフト波数の分布を有さないものとして温度分布を測定する従来の方法を用いても、測定誤差を低減することができる。勿論、上記の透過特性を有する光フィルタを用いた場合であっても、前述した第1,第2実施形態のように測定された温度分布を補正テーブルを用いて補正して、測定精度をより高めることも可能である。   If an optical filter having a transmission characteristic for making Stokes light and anti-Stokes light into a line spectrum is used as the optical filter 12b, the temperature distribution is measured assuming that neither the Stokes light nor the anti-Stokes light has a distribution of Raman shift wavenumbers. Even if a conventional method is used, measurement errors can be reduced. Of course, even when an optical filter having the above transmission characteristics is used, the measured temperature distribution is corrected using the correction table as in the first and second embodiments described above, so that the measurement accuracy is further improved. It can also be increased.

以上説明した通り、本発明の第1,第2実施形態では、ストークス光及び反ストークス光のラマンシフト波数の分布を考慮した補正テーブルを作成し、ストークス光及び反ストークス光が共にラマンシフト波数の分布を有さないものとして求められた温度分布を補正テーブルを用いて補正しているため、光ファイバ13の長さ方向における温度分布を高精度に測定することができる。また、本発明の第3実施形態のように、ラマンシフト波数の分布を有するストークス光及び反ストークス光を線スペクトル状にする光フィルタを用いれば、測定精度をより高めることができる。   As described above, in the first and second embodiments of the present invention, a correction table that takes into account the distribution of Raman shift wave numbers of Stokes light and anti-Stokes light is created, and both the Stokes light and anti-Stokes light have Raman shift wave numbers. Since the temperature distribution obtained as having no distribution is corrected using the correction table, the temperature distribution in the length direction of the optical fiber 13 can be measured with high accuracy. In addition, as in the third embodiment of the present invention, if an optical filter that makes Stokes light and anti-Stokes light having a distribution of Raman shift wavenumbers into a line spectrum is used, measurement accuracy can be further improved.

以上、本発明の実施形態による温度測定装置について説明したが、本発明は上記実施形態に制限されることなく、本発明の範囲内で自由に変更が可能である。例えば、上記実施形態では、パルス状のレーザ光Pを光ファイバ13に入射させて光ファイバ13の長さ方向における温度分布を測定する温度測定装置について説明した。しかしながら、本発明は、光ファイバ以外の被測定対象の温度を測定する温度測定装置、連続光を用いたラマン分光による温度測定装置にも適用することができる。 The temperature measuring device according to the embodiment of the present invention has been described above. However, the present invention is not limited to the above-described embodiment, and can be freely changed within the scope of the present invention. For example, in the above-described embodiment, the temperature measurement device that measures the temperature distribution in the length direction of the optical fiber 13 by causing the pulsed laser light PL to enter the optical fiber 13 has been described. However, the present invention can also be applied to a temperature measuring device that measures the temperature of an object to be measured other than an optical fiber, and a temperature measuring device using Raman spectroscopy that uses continuous light.

また、被測定対象に入射させるレーザ光の周波数を変化させながら被測定対象から得られる光の周波数特性・位相特性を求めることにより温度を測定するR−OFDR(Raman Optical Frequency Domain Reflectometry:光周波数領域反射測定装置)と言われる温度測定装置にも適用することができる。また、被測定対象に入射させるパルス状のレーザ光としてゴーレイ符号、バーカー符号、シンプレックス符号等を用いて符号変調した光パルス列を用いる温度測定装置にも適用することができる。   Further, an R-OFDR (Raman Optical Frequency Domain Reflectometry: optical frequency domain) that measures temperature by obtaining frequency characteristics and phase characteristics of light obtained from the measurement target while changing the frequency of the laser light incident on the measurement target. The present invention can also be applied to a temperature measuring device called a reflection measuring device. Further, the present invention can also be applied to a temperature measurement apparatus that uses an optical pulse train that is code-modulated using a Golay code, a Barker code, a simplex code, or the like as a pulsed laser beam incident on a measurement target.

更に、上記実施形態のように、ストークス光Pと反ストークス光Pとを測定して光ファイバ13の長手方向における温度分布を測定する温度測定装置のみならず、ストークス光のみ、反ストークス光のみ、或いはレーリー散乱光と反ストークス光とから温度を測定する温度測定装置にも本発明を適用することができる。つまり、本発明は、ラマン散乱光を用いて温度を測定するあらゆる温度測定装置に適用可能である。 Further, as described in the above embodiment, not only the temperature measuring device for measuring the temperature distribution in the longitudinal direction of the optical fiber 13 to measure the Stokes beam P s and the anti-Stokes light P a, Stokes light only, the anti-Stokes light The present invention can also be applied to a temperature measuring device that measures the temperature only from Rayleigh scattered light and anti-Stokes light. That is, the present invention can be applied to any temperature measuring apparatus that measures temperature using Raman scattered light.

1 温度測定装置
13 光ファイバ
21a,21b 光電変換回路
25 演算処理回路
26 メモリ
反ストークス光
レーザ光
ストークス光
1 a temperature measuring device 13 optical fiber 21a, 21b photoelectric conversion circuit 25 arithmetic processing circuit 26 memory P a anti-Stokes light P L laser beam P s Stokes light

Claims (6)

被測定対象にレーザ光を入射させて得られるラマン散乱光を用いて前記被測定対象の温度を測定する温度測定装置において、
前記ラマン散乱光を受光して受光信号を出力する受光部と、
前記ラマン散乱光のラマンシフト波数の分布を考慮した補正テーブルを記憶する記憶部と、
前記ラマン散乱光が前記ラマンシフト波数の分布を有さないものとして前記受光部から出力される受光信号を用いて前記被測定対象の温度を求めるための所定の演算を行い、当該所定の演算によって得られる演算結果を前記記憶部に記憶された前記補正テーブルを用いて補正する演算処理部と
を備えることを特徴とする温度測定装置。
In the temperature measuring device for measuring the temperature of the measurement object using Raman scattered light obtained by making laser light incident on the measurement object,
A light receiving unit that receives the Raman scattered light and outputs a light reception signal;
A storage unit for storing a correction table in consideration of the distribution of the Raman shift wave number of the Raman scattered light;
The Raman scattered light does not have the distribution of the Raman shift wave number, and performs a predetermined calculation for obtaining the temperature of the measurement target using a light reception signal output from the light receiving unit, and by the predetermined calculation A temperature measurement device comprising: an arithmetic processing unit that corrects an obtained calculation result using the correction table stored in the storage unit.
前記補正テーブルは、前記ラマン散乱光に含まれるストークス光及び反ストークス光のラマンシフト波数の分布を考慮して求められた温度毎の前記ストークス光と前記反ストークス光との強度比を示すテーブルであり、
前記演算処理部は、前記所定の演算として前記ストークス光と前記反ストークス光との強度比を求める演算を行い、当該演算により得られた強度比から求められる温度を前記補正テーブルを用いて補正する
ことを特徴とする請求項1記載の温度測定装置。
The correction table is a table showing an intensity ratio between the Stokes light and the anti-Stokes light for each temperature obtained in consideration of the distribution of Raman shift wave numbers of the Stokes light and the anti-Stokes light included in the Raman scattered light. Yes,
The arithmetic processing unit performs an operation for obtaining an intensity ratio between the Stokes light and the anti-Stokes light as the predetermined operation, and corrects a temperature obtained from the intensity ratio obtained by the operation using the correction table. The temperature measuring device according to claim 1.
前記補正テーブルは、前記ラマン散乱光に含まれるストークス光及び反ストークス光が前記ラマンシフト波数の分布を有さないものとした場合に求められる温度と、前記ストークス光及び前記反ストークス光のラマンシフト波数の分布を考慮した場合に求められる温度との差分を温度毎に示すテーブルであり、
前記演算処理部は、前記所定の演算として前記ストークス光と前記反ストークス光との強度比から前記被測定対象の温度を求める演算を行い、当該演算により得られた温度を前記補正テーブルを用いて補正する
ことを特徴とする請求項1記載の温度測定装置。
The correction table includes a temperature required when Stokes light and anti-Stokes light included in the Raman scattered light do not have a distribution of the Raman shift wavenumber, and a Raman shift of the Stokes light and the anti-Stokes light. It is a table showing the difference from the temperature obtained when considering the wave number distribution for each temperature,
The calculation processing unit performs a calculation for obtaining the temperature of the measurement target from the intensity ratio between the Stokes light and the anti-Stokes light as the predetermined calculation, and uses the correction table to calculate the temperature obtained by the calculation. The temperature measuring device according to claim 1, wherein the temperature measuring device is corrected.
前記演算処理部は、前記記憶部に記憶される前記補正テーブルを、前記ラマン散乱光のラマンシフト波数の分布を示すデータを用いて作成することを特徴とする請求項1から請求項3の何れか一項に記載の温度測定装置。   The said arithmetic processing part produces the said correction table memorize | stored in the said memory | storage part using the data which show distribution of the Raman shift wave number of the said Raman scattered light, Any one of Claim 1 to 3 characterized by the above-mentioned. The temperature measuring device according to claim 1. 前記ラマン散乱光に含まれるストークス光の所定の成分及び反ストークス光の所定の成分をそれぞれ個別に通過させるフィルタ部を備えることを特徴とする請求項1から請求項4の何れか一項に記載の温度測定装置。   5. The filter unit according to claim 1, further comprising a filter unit that individually passes a predetermined component of Stokes light and a predetermined component of anti-Stokes light included in the Raman scattered light. Temperature measuring device. 前記被測定対象は光ファイバであり、
前記受光部は、前記光ファイバにパルス状のレーザ光を入射させて得られる前記ラマン散乱光に含まれるストークス光及び反ストークス光を個別に受光し、
前記演算処理部は、前記受光部から順次出力される受光信号を用いて前記所定の演算を順次行うことにより、前記光ファイバの長さ方向における温度分布を求める
ことを特徴とする請求項1から請求項5の何れか一項に記載の温度測定装置。
The object to be measured is an optical fiber,
The light receiving unit individually receives Stokes light and anti-Stokes light included in the Raman scattered light obtained by allowing pulsed laser light to enter the optical fiber,
The said arithmetic processing part calculates | requires the temperature distribution in the length direction of the said optical fiber by performing the said predetermined calculation sequentially using the light reception signal sequentially output from the said light-receiving part. The temperature measuring device according to claim 5.
JP2009072637A 2009-03-24 2009-03-24 Temperature measuring device Pending JP2010223831A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2009072637A JP2010223831A (en) 2009-03-24 2009-03-24 Temperature measuring device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2009072637A JP2010223831A (en) 2009-03-24 2009-03-24 Temperature measuring device

Publications (1)

Publication Number Publication Date
JP2010223831A true JP2010223831A (en) 2010-10-07

Family

ID=43041167

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2009072637A Pending JP2010223831A (en) 2009-03-24 2009-03-24 Temperature measuring device

Country Status (1)

Country Link
JP (1) JP2010223831A (en)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2012027001A (en) * 2010-06-22 2012-02-09 Yokogawa Electric Corp Optical fiber temperature distribution measuring apparatus
CN103335742A (en) * 2013-07-12 2013-10-02 杭州欧忆光电科技有限公司 High-precision distributed optical fiber temperature sensing system
WO2014064770A1 (en) * 2012-10-23 2014-05-01 富士通株式会社 Abnormality detecting system and abnormality detecting method
CN105910728A (en) * 2016-06-16 2016-08-31 威海北洋光电信息技术股份公司 High-spatial-resolution Raman temperature measurement sensor and temperature measurement method
CN105973501A (en) * 2016-07-26 2016-09-28 威海北洋光电信息技术股份公司 Long-distance high-spatial-resolution Raman temperature measurement sensor and realization method thereof
CN109405995A (en) * 2018-09-17 2019-03-01 中国电子科技集团公司第五十五研究所 A kind of analysis method improving the junction temperature of chip measuring accuracy based on Raman spectroscopy

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2012027001A (en) * 2010-06-22 2012-02-09 Yokogawa Electric Corp Optical fiber temperature distribution measuring apparatus
WO2014064770A1 (en) * 2012-10-23 2014-05-01 富士通株式会社 Abnormality detecting system and abnormality detecting method
JP5867618B2 (en) * 2012-10-23 2016-02-24 富士通株式会社 Anomaly detection system and anomaly detection method
US9915620B2 (en) 2012-10-23 2018-03-13 Fujitsu Limited Abnormality detection system and abnormality detection method
CN103335742A (en) * 2013-07-12 2013-10-02 杭州欧忆光电科技有限公司 High-precision distributed optical fiber temperature sensing system
CN105910728A (en) * 2016-06-16 2016-08-31 威海北洋光电信息技术股份公司 High-spatial-resolution Raman temperature measurement sensor and temperature measurement method
CN105973501A (en) * 2016-07-26 2016-09-28 威海北洋光电信息技术股份公司 Long-distance high-spatial-resolution Raman temperature measurement sensor and realization method thereof
CN109405995A (en) * 2018-09-17 2019-03-01 中国电子科技集团公司第五十五研究所 A kind of analysis method improving the junction temperature of chip measuring accuracy based on Raman spectroscopy

Similar Documents

Publication Publication Date Title
EP2966426B1 (en) Optical fiber temperature distribution measuring device
US8858069B2 (en) Optical fiber temperature distribution measuring device
US9322721B2 (en) Optic fiber distributed temperature sensor system with self-correction function and temperature measuring method using thereof
EP2587238B1 (en) Optical fibre temperature distribution measurement apparatus
US7744275B2 (en) Optical fiber temperature sensing device
US9689751B2 (en) Optical fiber temperature distribution measurement device and method of measuring optical fiber temperature distribution
JP2010223831A (en) Temperature measuring device
JP2008051643A (en) Optical fiber temperature sensor
JP5743200B2 (en) Optical fiber temperature distribution measuring device
EP2944933A1 (en) Optical fiber temperature distribution measuring device
EP4097422A1 (en) Distributed optical fibre sensor
US9952104B2 (en) Optical fiber temperature distribution measuring device
KR101480101B1 (en) optical fiber bragg grating sensor system
JP5382364B2 (en) Optical fiber temperature distribution measuring device
US11994376B2 (en) Distributed vibration measuring device and method
JP6314678B2 (en) Optical fiber temperature distribution measuring device
JPH07243920A (en) Optical fiber temperature distribution measuring system
KR102590280B1 (en) apparatus for measuring distribution temperature by using optical fiber and method of correction temperature thereof
JP2020169889A (en) Temperature measuring device, temperature measuring method, and temperature measuring program
JP2977373B2 (en) Optical fiber temperature distribution sensor
JP4803083B2 (en) Optical fiber measuring apparatus and sampling method
JP3106443B2 (en) Temperature distribution measuring method and device therefor
US20160146642A1 (en) Optical fiber sensing optical system and optical fiber sensing system
JP6775024B2 (en) Optical wavelength measuring device
JP2011038871A (en) Temperature distribution measuring apparatus