JPH0626940A - Temperature distribution measuring method and apparatus - Google Patents

Temperature distribution measuring method and apparatus

Info

Publication number
JPH0626940A
JPH0626940A JP4182353A JP18235392A JPH0626940A JP H0626940 A JPH0626940 A JP H0626940A JP 4182353 A JP4182353 A JP 4182353A JP 18235392 A JP18235392 A JP 18235392A JP H0626940 A JPH0626940 A JP H0626940A
Authority
JP
Japan
Prior art keywords
light
optical fiber
stokes light
temperature 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.)
Granted
Application number
JP4182353A
Other languages
Japanese (ja)
Other versions
JP3106443B2 (en
Inventor
Toshinori Wakami
俊則 若見
Junichi Ota
順一 太田
Kazuhiro Okamoto
和弘 岡本
Kenji Tamano
研治 玉野
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.)
Sumitomo Electric Industries Ltd
Original Assignee
Sumitomo Electric Industries Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Sumitomo Electric Industries Ltd filed Critical Sumitomo Electric Industries Ltd
Priority to JP04182353A priority Critical patent/JP3106443B2/en
Publication of JPH0626940A publication Critical patent/JPH0626940A/en
Application granted granted Critical
Publication of JP3106443B2 publication Critical patent/JP3106443B2/en
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Landscapes

  • Measuring Temperature Or Quantity Of Heat (AREA)
  • Photometry And Measurement Of Optical Pulse Characteristics (AREA)

Abstract

(57)【要約】 【目的】 正確な温度分布の測定を行う。 【構成】 このシステムでは、ストークス光及び反スト
ークス光の検出強度を光ファイバの分光屈折率に応じた
サンプリング間隔でサンプリングし、前述したOTDR
の原理にもとづいて散乱点までの距離をもとめて温度分
布の測定が行われる。このシステムは、測定装置102
と、演算表示装置103で構成される。測定装置102
は、被測定光ファイバ101と光結合され、被測定光フ
ァイバ101に光PL を入射し、後方散乱光PB にて被
測定光ファイバ101の温度分布データを演算表示装置
103に出力するものである。演算表示装置103に
は、温度分布の表示,時間変化や異常高温の検出が行わ
れる。
(57) [Summary] [Purpose] Perform accurate temperature distribution measurement. [Configuration] In this system, the detected intensities of the Stokes light and the anti-Stokes light are sampled at sampling intervals according to the spectral refractive index of the optical fiber, and the OTDR described above is sampled.
Based on the principle of, the temperature distribution is measured based on the distance to the scattering point. This system comprises a measuring device 102
And the calculation display device 103. Measuring device 102
Is optically coupled to the optical fiber 101 to be measured, the light P L is incident on the optical fiber 101 to be measured, and the temperature distribution data of the optical fiber 101 to be measured is output to the arithmetic and display unit 103 as backscattered light P B. Is. The calculation display device 103 displays the temperature distribution and detects a temporal change and an abnormal high temperature.

Description

【発明の詳細な説明】Detailed Description of the Invention

【0001】[0001]

【産業上の利用分野】本発明は光ファイバ内で生じるラ
マン散乱光から光ファイバ上の各点の温度を測定するこ
とによって温度分布を測定する方法及び装置に関する。
BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to a method and apparatus for measuring temperature distribution by measuring the temperature of each point on an optical fiber from Raman scattered light generated in the optical fiber.

【0002】[0002]

【従来の技術】光ファイバを用いた従来技術として、O
TDR(Optical Time Domain Reflectometry )による
距離測定の原理と、ラマン散乱光の検出による温度測定
の原理と組み合わせた装置が知られ、例えば、文献「J.
P.Dakin et al:“Distrbuted optical fibre Raman tem
perature sensor using asemiconductor light sourcea
nd detection ”」に開示されている。
2. Description of the Related Art As a conventional technique using an optical fiber,
There is known an apparatus that combines the principle of distance measurement by TDR (Optical Time Domain Reflectometry) and the principle of temperature measurement by detection of Raman scattered light, for example, the document “J.
P. Dakin et al: “Distrbuted optical fiber Raman tem
perature sensor using asemiconductor light sourcea
nd detection "".

【0003】これを原理的に説明すると、図6(a)の
ようになる。まず、被測定光ファイバ101にパルス光
を入射すると、伝播の過程で後方散乱光が現れ、入射端
に戻ってくる。ここで、ある散乱点までの距離をL、パ
ルスの入射時点から後方散乱光の検出時点までの時間を
t、被測定光ファイバ101の屈折率をn、被測定光フ
ァイバ101中の光速をCとすると、 L=C・t/2n (1) となる。したがって、散乱点の位置が定量的に求まる。
This will be explained in principle as shown in FIG. 6 (a). First, when pulsed light is incident on the optical fiber 101 to be measured, backscattered light appears during the propagation process and returns to the incident end. Here, the distance to a certain scattering point is L, the time from the pulse incident time to the backscattered light detection time is t, the refractive index of the measured optical fiber 101 is n, and the speed of light in the measured optical fiber 101 is C. Then, L = C · t / 2n (1) Therefore, the position of the scattering point can be quantitatively obtained.

【0004】一方、後方散乱光にはレーリ光とストーク
ス光、反ストークス光が含まれ、入射パルス光の波長を
λo とするとレーリ光の波長はλo となり、ストークス
光の波長λs と反ストークス光の波長λasは、 1/λs =1/λo +ν (2a) 1/λas=1/λo −ν (2b) となる(図6(b)参照)。ここで、hはプランク定数
(J・S)、νはラマンシフト量(m-1)、kはボルツ
マン定数(J/K)である。
On the other hand, the backscattered light includes Rayleigh light, Stokes light, and anti-Stokes light. If the wavelength of the incident pulsed light is λo, the wavelength of Rayleigh light is λo, and the wavelength of Stokes light λs and that of anti-Stokes light are The wavelength λas is 1 / λs = 1 / λo + ν (2a) 1 / λas = 1 / λo−ν (2b) (see FIG. 6B). Here, h is Planck's constant (J · S), ν is Raman shift amount (m −1 ), and k is Boltzmann's constant (J / K).

【0005】ストークス光強度Is と反ストークス光強
度Iasの比は、被測定光ファイバ101中の散乱点の絶
対温度Tに依存し、Ias/Is がexp(−h・C・
ν/kT)に比例する関係となる。これにより、散乱点
の温度が定量的に求められる。そして、被測定光ファイ
バ101中の各散乱点からは絶対温度Tに依存したスト
ークス光及び反ストークス光が発生し、ストークス光強
度Is 及び反ストークス光強度Iasは温度分布に応じた
時間変化を示し、時間は各散乱点の距離に対応したもの
になっている。ストークス光強度Is 及び反ストークス
光強度Iasをサンプリングし、各サンプリング点で求め
た温度から被測定光ファイバ101中の温度分布が得ら
れることになる。
The ratio of the Stokes light intensity Is to the anti-Stokes light intensity Ias depends on the absolute temperature T of the scattering point in the optical fiber 101 to be measured, and Ias / Is is exp (-h.C.
The relationship is proportional to ν / kT). Thereby, the temperature of the scattering point is quantitatively obtained. Then, Stokes light and anti-Stokes light depending on the absolute temperature T are generated from each scattering point in the optical fiber 101 to be measured, and the Stokes light intensity Is and the anti-Stokes light intensity Ias show a time change according to the temperature distribution. , Time corresponds to the distance of each scattering point. The Stokes light intensity Is and the anti-Stokes light intensity Ias are sampled, and the temperature distribution in the optical fiber 101 to be measured can be obtained from the temperature obtained at each sampling point.

【0006】[0006]

【発明が解決しようとする課題】前述の従来技術におい
ては、ストークス光強度Is 及び反ストークス光強度I
asを同じ時間間隔でサンプリングし、同じ時間に対応す
るサンプリングデータの強度比から温度の計算を行って
いた。かつては光検出系の精度が悪く、長距離の測定で
はストークス光強度Is 及び反ストークス光強度Iasが
微弱なものとっていたので、それほど長い距離の測定に
は用いられなかった。現在、光検出系の精度が改善さ
れ、より微弱な光の測定し得るようになり、長い距離の
測定が可能になりつつある。
In the above-mentioned conventional technique, the Stokes light intensity Is and the anti-Stokes light intensity I are used.
As was sampled at the same time interval, and the temperature was calculated from the intensity ratio of the sampling data corresponding to the same time. In the past, the accuracy of the photodetection system was poor, and the Stokes light intensity Is and the anti-Stokes light intensity Ias were weak in long-distance measurement, so they were not used for measuring such a long distance. At present, the accuracy of the light detection system is improved, weaker light can be measured, and long distance measurement is becoming possible.

【0007】しかし、前述の従来技術では、測定する距
離がより長いものになると、正確な測定が行えなくなる
ことが明らかになってきた。この原因は、本件発明者の
鋭意検討の結果、ストークス光及び反ストークス光に対
する光ファイバの分光屈折率の違いに起因し、この屈折
率の違いによってストークス光と反ストークス光の散乱
点からの到達時間の違いが原因となっていることが分か
った。
However, in the above-mentioned conventional technique, it has become clear that accurate measurement cannot be performed if the distance to be measured becomes longer. The cause of this is that, as a result of diligent studies by the inventors of the present invention, it is caused by the difference in the spectral refractive index of the optical fiber with respect to Stokes light and anti-Stokes light, and due to this difference in refractive index, arrival from the scattering point of Stokes light and anti-Stokes light It turned out that the difference in time was the cause.

【0008】測定する距離がより長いものになると、例
えば、光ファイバの長さ方向にステップ状に温度が変化
している場合、ストークス光及び反ストークス光のサン
プリングデータ間に位置ずれがあるので、温度変化の前
後のデータの強度比を取ることがある。そのため、正確
な温度測定ができなくなるのである。また、正しく測定
されている場合では、ストークス光及び反ストークス光
の伝送損失は同じ量になるので、ストークス光強度Is
及び反ストークス光強度Iasの比をとることによって伝
送損失をキャンセルすることが可能である。しかし、前
述の従来技術では、同一の散乱点からのストークス光強
度Is 及び反ストークス光強度Iasでないので、伝送損
失をキャンセルすることができずに逆に誤差の原因にな
ることがある。
When the distance to be measured becomes longer, for example, when the temperature changes stepwise in the length direction of the optical fiber, there is a positional deviation between the sampling data of the Stokes light and the anti-Stokes light. The intensity ratio of the data before and after the temperature change may be taken. Therefore, accurate temperature measurement cannot be performed. Further, when measured correctly, the transmission loss of the Stokes light and the anti-Stokes light become the same amount, so the Stokes light intensity Is
It is possible to cancel the transmission loss by taking the ratio of the anti-Stokes light intensity Ias. However, in the above-described conventional technique, since the Stokes light intensity Is and the anti-Stokes light intensity Ias from the same scattering point are not used, the transmission loss cannot be canceled, which may cause an error.

【0009】[0009]

【課題を解決するための手段】上記課題を解決するため
に、本発明の温度分布測定方法は、光ファイバに入射し
たパルス光にて生じる後方散乱光に含まれるラマン散乱
光のうちストークス光と反ストークス光とを検出すると
共に、これらの検出強度のサンプリングをして検出強度
のデータをもとめ、データの比から後方散乱光の生じる
散乱点の温度を順次算出することによって光ファイバの
温度分布を測定する温度分布測定方法であって、ストー
クス光と反ストークス光に対する光ファイバの分光屈折
率に応じてストークス光及び反ストークス光の検出強度
のデータをもとめることを特徴とする。
In order to solve the above problems, the temperature distribution measuring method of the present invention uses Stokes light among Raman scattered lights included in backscattered light generated by pulsed light incident on an optical fiber. The temperature distribution of the optical fiber is calculated by detecting the anti-Stokes light, sampling the detected intensities, obtaining the detected intensity data, and sequentially calculating the temperature of the scattering point where the backscattered light is generated from the ratio of the data. A temperature distribution measuring method for measuring, which is characterized in that detection intensity data of Stokes light and anti-Stokes light is obtained according to the spectral refractive index of the optical fiber with respect to Stokes light and anti-Stokes light.

【0010】また、光ファイバの分光屈折率に応じたサ
ンプリング間隔でストークス光及び反ストークス光の検
出強度をサンプリングすることを特徴としても良い。
Further, the detection intensity of the Stokes light and the anti-Stokes light may be sampled at a sampling interval according to the spectral refractive index of the optical fiber.

【0011】さらに、ストークス光及び反ストークス光
の検出強度から光ファイバの分光屈折率に応じたサンプ
リング間隔に応じた補完データを求め、これらの補完デ
ータの比から後方散乱光の生じる散乱点の温度を順次算
出することによって被測定光ファイバの温度分布を測定
するを特徴としても良い。
Further, complementary data corresponding to the sampling interval corresponding to the spectral refractive index of the optical fiber is obtained from the detected intensities of the Stokes light and the anti-Stokes light, and the temperature of the scattering point where the backscattered light is generated is calculated from the ratio of these complementary data. The temperature distribution of the optical fiber to be measured may be measured by sequentially calculating

【0012】そして、本発明の温度分布測定装置は、光
ファイバに入射したパルス光にて生じる後方散乱光に含
まれるラマン散乱光のうちストークス光と反ストークス
光とを検出する第1の手段と、光ファイバの分光屈折率
に応じたサンプリング間隔でストークス光と反ストーク
ス光の検出強度のサンプリングをして検出強度のデータ
を出力する第2の手段と、データの比から後方散乱光の
生じる散乱点の温度を順次算出する第3の手段とを備え
る。
The temperature distribution measuring device of the present invention comprises first means for detecting Stokes light and anti-Stokes light among the Raman scattered light contained in the backscattered light generated by the pulsed light incident on the optical fiber. , Second means for sampling the detection intensity of the Stokes light and anti-Stokes light at a sampling interval according to the spectral refractive index of the optical fiber and outputting the detection intensity data, and the scattering of backscattered light from the ratio of the data And a third means for sequentially calculating the temperature of the point.

【0013】また、本発明の温度分布測定装置は、光フ
ァイバに入射したパルス光にて生じる後方散乱光に含ま
れるラマン散乱光のうちストークス光と反ストークス光
とを検出する第1の手段と、光ファイバの分光屈折率に
応じたサンプリング間隔でストークス光と反ストークス
光の検出強度のサンプリングをして検出強度のデータを
出力する第2の手段と、ストークス光及び反ストークス
光の検出強度から光ファイバの分光屈折率に応じたサン
プリング間隔に応じた補完データを求め、これらの補完
データの比から後方散乱光の生じる散乱点の温度を順次
算出する第3の手段とを備える。
Further, the temperature distribution measuring apparatus of the present invention comprises first means for detecting Stokes light and anti-Stokes light among Raman scattered light contained in backscattered light generated by pulsed light incident on the optical fiber. From the detection intensity of the Stokes light and the anti-Stokes light, the second means for sampling the detection intensity of the Stokes light and the anti-Stokes light at a sampling interval according to the spectral refractive index of the optical fiber and outputting the detection intensity data, Third means for obtaining complementary data according to the sampling interval according to the spectral refractive index of the optical fiber and sequentially calculating the temperature of the scattering point where the backscattered light is generated from the ratio of these complementary data.

【0014】[0014]

【作用】本発明の温度分布測定方法では、パルス光が光
ファイバ内部の各所で散乱し、各散乱点で生じたストー
クス光と反ストークス光とが光ファイバ内部を検出され
る地点まで伝搬して検出される。そして、これらの検出
強度のサンプリングをして検出強度のデータをもとめら
れる。このとき、ストークス光及び反ストークス光に対
して光ファイバの分光屈折率は異なった値を持つので、
散乱点から検出される地点までの伝搬時間は異なった値
になっているのだが、ストークス光と反ストークス光に
対する光ファイバの分光屈折率に応じたストークス光及
び反ストークス光の検出強度のデータがもとめられる。
即ち、光ファイバの分光屈折率に応じたサンプリング間
隔でサンプリングを行うか、或いは、補完データを求め
ることで、温度を算出するデータは各散乱点に対応した
ものになる。これによって光ファイバの分光屈折率に影
響されないより正確な温度分布がもとめられる。
In the temperature distribution measuring method of the present invention, pulsed light is scattered at various points inside the optical fiber, and Stokes light and anti-Stokes light generated at each scattering point propagate to the point where it is detected inside the optical fiber. To be detected. Then, the detection intensity data can be obtained by sampling these detection intensities. At this time, since the spectral refractive index of the optical fiber has different values for Stokes light and anti-Stokes light,
Although the propagation time from the scattering point to the detected point has different values, the detection intensity data of Stokes light and anti-Stokes light according to the spectral refractive index of the optical fiber for Stokes light and anti-Stokes light are available. It is sought after.
That is, by performing sampling at sampling intervals according to the spectral refractive index of the optical fiber or obtaining complementary data, the data for calculating the temperature corresponds to each scattering point. As a result, a more accurate temperature distribution that is not affected by the spectral refractive index of the optical fiber can be obtained.

【0015】本発明の温度分布測定装置では、上述の温
度分布測定方法で光ファイバ内の各散乱点の温度が順次
算出され、光ファイバの分光屈折率に影響されないより
正確な温度分布がもとめられる。
In the temperature distribution measuring device of the present invention, the temperature of each scattering point in the optical fiber is sequentially calculated by the above-mentioned temperature distribution measuring method, and a more accurate temperature distribution which is not affected by the spectral refractive index of the optical fiber can be obtained. .

【0016】[0016]

【実施例】本発明の実施例を図面を参照して説明する。
図1は、本発明の第1の方法を実現するための一実施例
の構成を示したものである。この光ファイバの温度分布
測定システムでは、ストークス光及び反ストークス光の
検出強度を光ファイバの分光屈折率に応じたサンプリン
グ間隔でサンプリングし、前述したOTDRの原理にも
とづいて散乱点までの距離をもとめて温度分布の測定が
行われる。このシステムは、測定装置102と、演算表
示装置103で構成される。
Embodiments of the present invention will be described with reference to the drawings.
FIG. 1 shows the configuration of an embodiment for implementing the first method of the present invention. In this optical fiber temperature distribution measurement system, the detection intensities of Stokes light and anti-Stokes light are sampled at sampling intervals according to the spectral refractive index of the optical fiber, and the distance to the scattering point is determined based on the above-mentioned OTDR principle. Then, the temperature distribution is measured. This system includes a measuring device 102 and a calculation display device 103.

【0017】測定装置102は、被測定光ファイバ10
1と光結合され、被測定光ファイバ101に光PL を入
射し、後方散乱光PB にて被測定光ファイバ101の温
度分布データを演算表示装置103に出力するものであ
る。演算表示装置103には、いわゆるパーソナルコン
ピュータが用いられ、RS232C又はGP−IBなど
の規格化されたインターフェイスにより測定装置102
と温度や距離のデータなどをやり取りし、被測定光ファ
イバ101の温度分布の表示,時間変化や異常高温の検
出が行われる。測定装置102の各部の構成はつぎのよ
うになっている。
The measuring device 102 includes an optical fiber 10 to be measured.
The light P L is incident on the optical fiber 101 to be measured, and the backscattered light P B outputs the temperature distribution data of the optical fiber 101 to the arithmetic and display unit 103. A so-called personal computer is used as the arithmetic display device 103, and the measuring device 102 is provided by a standardized interface such as RS232C or GP-IB.
The temperature distribution and the temperature data are exchanged to display the temperature distribution of the optical fiber 101 to be measured and to detect a temporal change and an abnormal high temperature. The structure of each part of the measuring device 102 is as follows.

【0018】タイミング発生器132は、水晶発振機と
分周器とにより構成され、レーザパルス光PL のタイミ
ングトリガ(レーザダイオード121が発光する時間を
制御するタイミングパルス)をレーザドライバ128に
順次与える。また、レーザダイオード121の1回の発
光ごとに、A/D変換器131s,a,平均化器127
s,a,信号処理器126へ測定開始のタイミングパル
スを与える。レーザドライバ128は、タイミングトリ
ガにしたがってレーザダイオード121を駆動し、パル
ス光PL を発生させている。ビームスプリッタ122
は、パルス光PLを透過して被測定光ファイバ101に
与え、被測定光ファイバ101からの後方散乱光PB
波長分離器123に与える。
The timing generator 132 is composed of a crystal oscillator and a frequency divider, and sequentially supplies the laser driver 128 with a timing trigger of the laser pulse light P L (a timing pulse for controlling the time when the laser diode 121 emits light). . In addition, each time the laser diode 121 emits light, the A / D converters 131 s, a and the averaging device 127
A timing pulse for starting measurement is given to s, a and the signal processor 126. The laser driver 128 drives the laser diode 121 according to the timing trigger to generate the pulsed light P L. Beam splitter 122
Transmits the pulsed light P L to the measured optical fiber 101 and supplies the backscattered light P B from the measured optical fiber 101 to the wavelength separator 123.

【0019】被測定光ファイバ101の入力部にはマル
チモード光ファイバをコイル状に巻回した屈曲部111
が設けられ、高次光が取り除かれるようになっている。
波長分離器123は、被測定光ファイバ101からの後
方散乱光PB からストークス光Ps,反ストークス光P
aを分離する。アバランシェフォトダイオードによる光
検出器APDs,APDaは、ストークス光Ps,反ス
トークス光Paの検出強度を電気信号にして出力し、増
幅器125s,aはこれらの検出強度を増幅する。
At the input portion of the optical fiber 101 to be measured, a bent portion 111 in which a multimode optical fiber is wound in a coil shape.
Is provided so that higher-order light can be removed.
The wavelength demultiplexer 123 converts the backscattered light P B from the measured optical fiber 101 into Stokes light Ps and anti-Stokes light P B.
Separate a. The photodetectors APDs and APDa with avalanche photodiodes output the detected intensities of the Stokes light Ps and the anti-Stokes light Pa as electric signals, and the amplifiers 125s and a amplify the detected intensities.

【0020】A/D変換器131s,aは、ストークス
光Ps,反ストークス光Paの検出強度をディジタルデ
ータに変換するもので、高速タイプのもので構成され
る。A/D変換器131s,aのサンプリングクロック
は外部(この場合、演算表示装置103)からの入力に
よって可変であり、設定可能になっている。A/D変換
器131s,aのサンプリング間隔Ts,Taは、それ
ぞれストークス光Ps,反ストークス光Paに対する光
ファイバの分光屈折率に応じてつぎの式(3)で示され
る関係の値に設定される。
The A / D converters 131s and a are for converting the detection intensities of the Stokes light Ps and the anti-Stokes light Pa into digital data, and are of high speed type. The sampling clock of the A / D converters 131s, a is variable and settable by an input from the outside (in this case, the arithmetic display device 103). The sampling intervals Ts and Ta of the A / D converters 131s and a are set to the values of the relationship shown by the following equation (3) according to the spectral refractive index of the optical fiber with respect to the Stokes light Ps and the anti-Stokes light Pa, respectively. It

【0021】[0021]

【数1】 [Equation 1]

【0022】平均化器127s,aは、A/D変換器1
31s,aからのディジタルデータについてそれぞれ数
万回の加算平均を行い、ランダムなノイズ成分を相殺し
S/N比の改善を行う。信号処理器126は、平均化器
127s,aでノイズを除去されたストークス光Ps,
反ストークス光Paの検出強度のディジタルデータを演
算処理して温度分布の波形のデータに変換する。このと
きの演算は、A/D変換器131s,aそれぞれに対応
する各サンプリングクロックついて行われ、その結果を
信号処理器126は演算表示装置103にインターフェ
イス仕様に合わせて出力する。
The averaging unit 127s, a is an A / D converter 1
The digital data from 31s and a are averaged tens of thousands of times to cancel random noise components and improve the S / N ratio. The signal processor 126 includes the Stokes light Ps, from which the noise is removed by the averaging unit 127s, a.
The digital data of the detected intensity of the anti-Stokes light Pa is arithmetically processed to be converted into temperature distribution waveform data. The calculation at this time is performed for each sampling clock corresponding to each of the A / D converters 131s and a, and the signal processor 126 outputs the result to the calculation display device 103 in accordance with the interface specifications.

【0023】つぎにこのシステムの動作について説明す
る。
Next, the operation of this system will be described.

【0024】半導体レーザ121からのパルス光PL
ビームスプリッタ122を透過して被測定光ファイバ1
01に入射される。被測定光ファイバ101からの後方
散乱光は、高次光が取り除かれたのちビームスプリッタ
122で光路変換され、波長分離部123で波長λs の
ストークス光と波長λasの反ストークス光に分離され
る。光検出器APDs,APDaでそれぞれ検出強度が
電気信号に変換され増幅される。
The pulsed light P L from the semiconductor laser 121 passes through the beam splitter 122 and the measured optical fiber 1
It is incident on 01. The backscattered light from the optical fiber 101 to be measured is subjected to optical path conversion by the beam splitter 122 after the higher-order light is removed, and is separated into the Stokes light having the wavelength λs and the anti-Stokes light having the wavelength λas by the wavelength separation unit 123. The detection intensity is converted into an electric signal and amplified by each of the photodetectors APDs and APDa.

【0025】ここで、ストークス光,反ストークス光に
ついて光ファイバの分光屈折率が異なったものである
と、光ファイバ内での伝搬速度が異なったものになる。
パルス光PL がトリガされてからの時間tにおいて観測
されるストークス光,反ストークス光の発生地点(散乱
点)までの距離ls,laは理論的に式(4)で示され
るものになる。これらの検出強度の波形は散乱点までの
距離に応じて時間的にずれたものになっている。これら
の波形の時間的ずれΔtはつぎの式(5)で示される。
Here, if the spectral refractive index of the optical fiber is different for Stokes light and anti-Stokes light, the propagation velocities in the optical fiber will be different.
The distances ls, la to the generation points (scattering points) of the Stokes light and anti-Stokes light observed at time t after the pulsed light P L is triggered are theoretically expressed by the equation (4). The waveforms of these detected intensities are temporally displaced according to the distance to the scattering point. The time shift Δt of these waveforms is expressed by the following equation (5).

【0026】[0026]

【数2】 [Equation 2]

【0027】[0027]

【数3】 [Equation 3]

【0028】このシステムでは、A/D変換器131
s,aのサンプリングタイムに前述の式(3)の関係を
持たせてサンプリングを行い、対応するストークス検出
強度Is と反ストークス検出強度Iasのデータから温度
分布を求めている。図2は、パルス光PL の波長λを
0.91μm、ストークス光の波長λsを0.94μ
m、反ストークス光の波長λasを0.88μm、分光屈
折率na ,ns をそれぞれ1.4800,1.4786
とした場合の一例を示したものである。
In this system, the A / D converter 131
Sampling times of s and a are related by the above equation (3) to perform sampling, and the temperature distribution is obtained from the corresponding data of the Stokes detection intensity Is and the anti-Stokes detection intensity Ias. In FIG. 2, the wavelength λ of the pulsed light P L is 0.91 μm, and the wavelength λs of the Stokes light is 0.94 μm.
m, the wavelength λ as of the anti-Stokes light is 0.88 μm, and the spectral refractive indices n a and n s are 1.4800 and 1.4786, respectively.
This is an example of the case.

【0029】被測定光ファイバ101の温度分布が図2
(a)に示すものであると、ストークス検出強度Is と
反ストークス検出強度Iasは、図2(b),(c)に示
すような互いに時間的ずれを持つ波形になっていて、図
3は、このずれを理論(実線)と実験(プロット)を比
較して示したものであり、前述の式(5)と良い一致を
しめしている。A/D変換器131s,aのサンプリン
グタイムは前述の式(3)の関係を持ち、A/D変換器
131s,aでは、ストークス検出強度Is と反ストー
クス検出強度Iasのサンプリングがそれぞれ行われる。
そして、A/D変換されて、同じ散乱点に対応するスト
ークス光と反ストークスの検出強度Is,Iasのディジ
タルデータが得られる。
The temperature distribution of the optical fiber 101 to be measured is shown in FIG.
As shown in FIG. 3A, the Stokes detection intensity Is and the anti-Stokes detection intensity Ias are waveforms having a time lag with respect to each other as shown in FIGS. 2B and 2C, and FIG. This deviation is shown by comparing the theory (solid line) and the experiment (plot), and shows good agreement with the above-mentioned equation (5). The sampling times of the A / D converters 131s and a have the relationship of the above-mentioned formula (3), and the A / D converters 131s and a sample the Stokes detection intensity Is and the anti-Stokes detection intensity Ias, respectively.
Then, A / D conversion is performed to obtain digital data of detected intensities Is and Ias of Stokes light and anti-Stokes corresponding to the same scattering point.

【0030】平均化器127s,aでノイズ成分を相殺
してS/N比が改善される。そして、信号処理器126
で、図2(b),(c)の横軸において時間t=…,n
−1,n,n+1,…に対応する検出強度Is ,Iasの
ディジタルデータの比を求めて温度が順次算出される。
即ち、同じ散乱点に対応する検出強度Is ,Iasのディ
ジタルデータの比から温度が順次求められる。時間tに
対応する検出強度比Is /Iasは被測定光ファイバ10
1上の位置に対応した温度となっており、検出強度Is
,Iasは同じ散乱点に対応するものとなっている。そ
のため、順次求められた温度は、前述の式(5)で示さ
れた時間ずれを含まないものになっており、これから正
確な温度分布が得られる。この正確な温度分布が、演算
表示装置103に出力されて表示される(図2
(d))。図2(d)は図2(a)と非常に良い一致を
示しており、非常に正確に温度分布が測定されている。
The noise component is canceled by the averaging unit 127s, a to improve the S / N ratio. Then, the signal processor 126
Then, on the horizontal axis of FIGS. 2B and 2C, time t = ..., N
The temperatures are sequentially calculated by obtaining the ratio of the digital data of the detected intensities Is and Ias corresponding to -1, n, n + 1, ....
That is, the temperature is sequentially obtained from the ratio of the digital data of the detection intensities Is and Ias corresponding to the same scattering point. The detected intensity ratio Is / Ias corresponding to the time t is the measured optical fiber 10
The temperature corresponds to the position above 1, and the detection intensity Is
, Ias correspond to the same scattering point. Therefore, the temperatures sequentially obtained do not include the time lag indicated by the above equation (5), and an accurate temperature distribution can be obtained from this. This accurate temperature distribution is output to and displayed on the arithmetic and display unit 103 (FIG. 2).
(D)). FIG. 2D shows a very good agreement with FIG. 2A, and the temperature distribution is measured very accurately.

【0031】検出強度Is ,Iasの波形のずれは、被測
定光ファイバ101が長いほど伝搬時間が掛かるので大
きいものになる。そのため、大きな距離についての温度
分布を求める際に、長距離のものほど検出強度Is ,I
asのずれが大きくなるのだが、このシステムにおいて
は、被測定光ファイバ101上の位置に対応した検出強
度Is ,Iasのディジタルデータの比で温度を求めるの
で、長距離のものでも正確でその効果は大きい。また、
光ファイバは波長λが1.3μm近傍において波長分散
は小さいのであるが、本発明では、これとは離れた波長
のパルス光PL を用いても正確に測定でき、温度測定に
最も良好な波長のものを用いることができる。
The deviation between the waveforms of the detection intensities Is and Ias becomes large because the longer the optical fiber 101 to be measured, the longer the propagation time. Therefore, when obtaining the temperature distribution for a large distance, the longer the distance is, the higher the detection intensities Is and I are.
Although the deviation of as becomes large, in this system, the temperature is obtained by the ratio of the digital data of the detected intensities Is and Ias corresponding to the position on the optical fiber 101 to be measured, so that the effect is accurate even for long distances. Is big. Also,
The optical fiber has a small chromatic dispersion in the vicinity of the wavelength λ of 1.3 μm, but in the present invention, it is possible to accurately measure even the pulsed light P L having a wavelength apart from this wavelength, and it is the best wavelength for temperature measurement. Can be used.

【0032】つぎに、本発明の第2の方法を実現するた
めの実施例について説明する。
Next, an embodiment for realizing the second method of the present invention will be described.

【0033】第2の方法による温度分布測定システムで
は、サンプリング間隔でサンプリングされたストークス
光及び反ストークス光の検出強度から光ファイバの分光
屈折率に応じて同じ散乱点に対応する検出強度を求める
ものである。このシステムは、前述の図1と同等の構成
を持つが、測定装置102のA/D変換器131s,a
のサンプリングクロックを同じものとしている点と、信
号処理器126での処理とが異なっている。
In the temperature distribution measuring system according to the second method, the detection intensity corresponding to the same scattering point is obtained from the detection intensities of the Stokes light and the anti-Stokes light sampled at the sampling intervals according to the spectral refractive index of the optical fiber. Is. This system has the same configuration as that of FIG. 1 described above, but the A / D converters 131s, a of the measuring apparatus 102 are used.
The same sampling clock is used, and the processing by the signal processor 126 is different.

【0034】この信号処理器126の演算は、ソフトウ
ェア的になされるもので、平均化器127s,aでノイ
ズを除去された反ストークス光Paの検出強度のデータ
から比例配分若しくは近傍3点の2次曲線近似によって
ストークス光Psと同じ散乱点に対応する検出強度を求
め、これと対応するストークス光Psの検出強度とから
温度分布を求めることで為される。この演算は、信号処
理器126ではなく、演算表示装置103で行っても等
価である。これによって前述の式(5)の波形のずれが
補正され、正確な温度分布が信号処理器126に出力さ
れる。
The calculation of the signal processor 126 is performed by software, and is proportionally distributed from the detected intensity data of the anti-Stokes light Pa from which the noise is removed by the averaging units 127s and a, or two points of three neighboring points. This is performed by obtaining the detected intensity corresponding to the same scattering point as the Stokes light Ps by the following curve approximation, and obtaining the temperature distribution from this and the detected intensity of the corresponding Stokes light Ps. This calculation is equivalent even if the calculation display device 103 is used instead of the signal processor 126. As a result, the deviation of the waveform of the above equation (5) is corrected, and the accurate temperature distribution is output to the signal processor 126.

【0035】図4は、この実施例について前述の図2と
同様のものを示したものである。被測定光ファイバ10
1の温度分布が図4(a)に示すものであると、ストー
クス検出強度Is と反ストークス検出強度Iasは、図4
(b),(c)に示すような互いに時間的ずれを持つ波
形になる。この図の横軸において時間t=…,n−1,
n,n+1,…でサンプリングされA/D変換されて、
検出強度Is ,Iasのディジタルデータが得られる。そ
して、検出強度Iasのデータの補完データが信号処理器
126で求められる。比例配分によるものとすれば、時
間t=nに対応するデータは、時間t=n−1,nのデ
ータと時間ずれΔtから求められ((c)の点C)、2
次曲線近似であれば、時間t=n−1,n,n+1から
求められる。
FIG. 4 shows the same as FIG. 2 described above for this embodiment. Optical fiber to be measured 10
When the temperature distribution of No. 1 is as shown in FIG. 4A, the Stokes detection intensity Is and the anti-Stokes detection intensity Ias are as shown in FIG.
The waveforms have time lags as shown in (b) and (c). On the horizontal axis of this figure, time t = ..., N-1,
is sampled by n, n + 1, ..., A / D converted,
Digital data of the detection intensities Is and Ias can be obtained. Then, the complementary data of the data of the detection intensity Ias is obtained by the signal processor 126. If proportional distribution is used, the data corresponding to time t = n is obtained from the data at time t = n−1, n and the time difference Δt (point C in (c)), 2
If the curve is approximated, the time t = n-1, n, n + 1 can be obtained.

【0036】このようにして得た補完データとストーク
ス光の時間t=nの検出強度Is とから該当する散乱点
の温度が求められる((c)の点D)。(d)の点
C’,D’は、点C,Dの値による温度を示したもの
で、近似曲線の次数の高いものほど正確さを増してい
る。補完データは、検出強度Iasの波形に含まれる時間
ずれΔtを補正し、前述の実施例のものとほぼ等価なも
のになっている。そのため、これを用いて比を求めるこ
とで温度分布の誤差が抑えられ正確な温度分布の測定を
可能にしている(図4(d))。
The temperature of the corresponding scattering point can be obtained from the complementary data thus obtained and the detected intensity Is of the Stokes light at time t = n (point D in (c)). Points C ′ and D ′ in (d) represent temperatures according to the values of points C and D, and the accuracy increases as the order of the approximate curve increases. The complementary data corrects the time shift Δt included in the waveform of the detected intensity Ias, and is substantially equivalent to that of the above-described embodiment. Therefore, by calculating the ratio using this, an error in the temperature distribution can be suppressed and an accurate temperature distribution can be measured (FIG. 4 (d)).

【0037】図5は、補完データを求めないで検出強度
比Is /Iasを取ったものについて示したもので、図2
または図4との比較からも分かるように正確な温度分布
が得られないものになっている。
FIG. 5 shows the detection intensity ratio Is / Ias obtained without obtaining the complementary data.
Alternatively, as can be seen from the comparison with FIG. 4, an accurate temperature distribution cannot be obtained.

【0038】この実施例では、前述の実施例の効果に加
えて、ソフトウェア的な手段によって補完データを得て
いるので、ハードウェアに変更を加えることがない。そ
のため、市販品を使用することができ、特別な費用を要
しない。また、分光屈折率のことなるファイバを用いて
も、その屈折率の補正を簡単に行うことができる。
In this embodiment, in addition to the effects of the above-mentioned embodiment, since complementary data is obtained by software means, no change is made to hardware. Therefore, a commercially available product can be used without any special cost. Further, even if a fiber having a different spectral refractive index is used, the refractive index can be easily corrected.

【0039】[0039]

【発明の効果】以上の通り本発明の温度分布測定方法及
び装置によれば、温度を算出するデータは各散乱点に対
応したものになるので、長い距離の測定であっても、ス
トークス光及び反ストークス光の伝搬時間の差が補償さ
れて光ファイバの分光屈折率に影響されないより正確な
温度分布がもとめられる。
As described above, according to the temperature distribution measuring method and apparatus of the present invention, since the data for calculating the temperature corresponds to each scattering point, the Stokes light and the Stokes light can be measured even at a long distance. A more accurate temperature distribution that is not affected by the spectral refractive index of the optical fiber is obtained by compensating for the difference in the propagation time of the anti-Stokes light.

【図面の簡単な説明】[Brief description of drawings]

【図1】実施例の構成図。FIG. 1 is a configuration diagram of an embodiment.

【図2】温度分布の測定の説明図。FIG. 2 is an explanatory diagram of measurement of temperature distribution.

【図3】波形のずれを示す図。FIG. 3 is a diagram showing a waveform shift.

【図4】温度分布の測定の説明図。FIG. 4 is an explanatory diagram of measurement of temperature distribution.

【図5】温度分布の測定の説明図。FIG. 5 is an explanatory diagram of measurement of temperature distribution.

【図6】従来の測定を示す図。FIG. 6 is a diagram showing a conventional measurement.

【符号の説明】[Explanation of symbols]

101…光ファイバ、102…測定装置、103…演算
表示装置、123…波長分離器、APDs,APDa…
光検出器、126…信号処理器、131s,a…A/D
変換器、132…タイミング発生器。
101 ... Optical fiber, 102 ... Measuring device, 103 ... Arithmetic display device, 123 ... Wavelength separator, APDs, APDa ...
Photodetector, 126 ... Signal processor, 131s, a ... A / D
Converter, 132 ... Timing generator.

───────────────────────────────────────────────────── フロントページの続き (72)発明者 玉野 研治 神奈川県横浜市栄区田谷町1番地 住友電 気工業株式会社横浜製作所内 ─────────────────────────────────────────────────── ─── Continuation of the front page (72) Inventor Kenji Tamano 1 Taya-cho, Sakae-ku, Yokohama-shi, Kanagawa Sumitomo Electric Industries, Ltd. Yokohama Works

Claims (5)

【特許請求の範囲】[Claims] 【請求項1】 光ファイバに入射したパルス光にて生じ
る後方散乱光に含まれるラマン散乱光のうちストークス
光と反ストークス光とを検出すると共に、これらの検出
強度のサンプリングをして前記検出強度のデータをもと
め、前記データの比から前記後方散乱光の生じる散乱点
の温度を順次算出することによって前記光ファイバの温
度分布を測定する温度分布測定方法であって、 前記ストークス光と前記反ストークス光に対する前記光
ファイバの分光屈折率に応じて前記ストークス光及び前
記反ストークス光の前記検出強度のデータをもとめるこ
とを特徴とする温度分布測定方法。
1. Stokes light and anti-Stokes light among Raman scattered light included in backscattered light generated by pulsed light incident on an optical fiber are detected, and the detected intensities are sampled to detect the detected intensities. Is a temperature distribution measuring method for measuring the temperature distribution of the optical fiber by sequentially calculating the temperature of the scattering point where the backscattered light is generated from the ratio of the data, the Stokes light and the anti-Stokes A temperature distribution measuring method, wherein data of the detected intensities of the Stokes light and the anti-Stokes light is obtained according to a spectral refractive index of the optical fiber with respect to light.
【請求項2】 前記光ファイバの分光屈折率に応じたサ
ンプリング間隔で前記ストークス光及び前記反ストーク
ス光の検出強度をサンプリングすることを特徴とする請
求項1記載の温度分布測定方法。
2. The temperature distribution measuring method according to claim 1, wherein the detection intensities of the Stokes light and the anti-Stokes light are sampled at sampling intervals according to the spectral refractive index of the optical fiber.
【請求項3】 前記ストークス光及び前記反ストークス
光の検出強度から前記光ファイバの分光屈折率に応じた
サンプリング間隔に応じた補完データを求め、これらの
補完データの比から前記後方散乱光の生じる散乱点の温
度を順次算出することによって前記被測定光ファイバの
温度分布を測定することを特徴とする請求項1記載の温
度分布測定方法。
3. Complementary data corresponding to a sampling interval corresponding to the spectral refractive index of the optical fiber is obtained from the detected intensities of the Stokes light and the anti-Stokes light, and the backscattered light is generated from the ratio of these complementary data. The temperature distribution measuring method according to claim 1, wherein the temperature distribution of the optical fiber to be measured is measured by sequentially calculating the temperatures of the scattering points.
【請求項4】 光ファイバに入射したパルス光にて生じ
る後方散乱光に含まれるラマン散乱光のうちストークス
光と反ストークス光とを検出する第1の手段と、 前記光ファイバの分光屈折率に応じたサンプリング間隔
で前記ストークス光と反ストークス光の検出強度のサン
プリングをして前記検出強度のデータを出力する第2の
手段と、 前記データの比から前記後方散乱光の生じる散乱点の温
度を順次算出する第3の手段とを備える温度分布測定装
置。
4. A first means for detecting Stokes light and anti-Stokes light among Raman scattered light included in backscattered light generated by pulsed light incident on an optical fiber; and a spectral refractive index of the optical fiber. Second means for sampling the detection intensities of the Stokes light and the anti-Stokes light at a sampling interval according to the output and outputting the data of the detection intensities, and the temperature of the scattering point where the backscattered light is generated from the ratio of the data. A temperature distribution measuring device comprising a third means for sequentially calculating.
【請求項5】 光ファイバに入射したパルス光にて生じ
る後方散乱光に含まれるラマン散乱光のうちストークス
光と反ストークス光とを検出する第1の手段と、 前記光ファイバの分光屈折率に応じたサンプリング間隔
で前記ストークス光と反ストークス光の検出強度のサン
プリングをして前記検出強度のデータを出力する第2の
手段と、 前記ストークス光及び前記反ストークス光の検出強度か
ら前記光ファイバの分光屈折率に応じたサンプリング間
隔に応じた補完データを求め、これらの補完データの比
から前記後方散乱光の生じる散乱点の温度を順次算出す
る第3の手段とを備える温度分布測定装置。
5. A first means for detecting Stokes light and anti-Stokes light among Raman scattered light included in backscattered light generated by pulsed light incident on an optical fiber; and a spectral refractive index of the optical fiber. Second means for sampling the detection intensities of the Stokes light and the anti-Stokes light at a sampling interval according to the second means for outputting data of the detection intensity, and detecting the intensity of the Stokes light and the anti-Stokes light from the optical fiber of the optical fiber. A temperature distribution measuring device comprising: third means for obtaining complementary data corresponding to a sampling interval corresponding to the spectral refractive index, and sequentially calculating the temperature of the scattering point where the backscattered light is generated from the ratio of these complementary data.
JP04182353A 1992-07-09 1992-07-09 Temperature distribution measuring method and device therefor Expired - Lifetime JP3106443B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP04182353A JP3106443B2 (en) 1992-07-09 1992-07-09 Temperature distribution measuring method and device therefor

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP04182353A JP3106443B2 (en) 1992-07-09 1992-07-09 Temperature distribution measuring method and device therefor

Publications (2)

Publication Number Publication Date
JPH0626940A true JPH0626940A (en) 1994-02-04
JP3106443B2 JP3106443B2 (en) 2000-11-06

Family

ID=16116832

Family Applications (1)

Application Number Title Priority Date Filing Date
JP04182353A Expired - Lifetime JP3106443B2 (en) 1992-07-09 1992-07-09 Temperature distribution measuring method and device therefor

Country Status (1)

Country Link
JP (1) JP3106443B2 (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008232888A (en) * 2007-03-22 2008-10-02 Yokogawa Electric Corp Optical fiber measuring apparatus and sampling method
JP2017120210A (en) * 2015-12-28 2017-07-06 鹿島建設株式会社 Distortion distribution data processing device and distortion distribution data processing method

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6787029B2 (en) 2001-08-31 2004-09-07 Cabot Corporation Material for chromatography

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008232888A (en) * 2007-03-22 2008-10-02 Yokogawa Electric Corp Optical fiber measuring apparatus and sampling method
JP2017120210A (en) * 2015-12-28 2017-07-06 鹿島建設株式会社 Distortion distribution data processing device and distortion distribution data processing method

Also Published As

Publication number Publication date
JP3106443B2 (en) 2000-11-06

Similar Documents

Publication Publication Date Title
US5272334A (en) Device for detecting a temperature distribution having a laser source with temperature and wavelength adjusting means
CN111896137B (en) A distributed optical fiber Raman sensing device and method with centimeter-level spatial resolution
US20210116311A1 (en) Self-calibration detection device and temperature demodulation method oriented to fiber raman temperature sensing system
US6486961B1 (en) System and method for measuring group delay based on zero-crossings
JP2006267107A (en) Multimode optical fiber delay time difference measuring device
US8897608B2 (en) High sampling resolution DTS system and method
CN110207733B (en) Optical fiber interferometer arm length difference measuring device and method based on sweep frequency laser
JP3106443B2 (en) Temperature distribution measuring method and device therefor
US20080018907A1 (en) Optical apparatus and method for distance measuring
JP2584478B2 (en) Method for processing received signal of optical fiber backscattered light
JPH08159882A (en) Method and apparatus for measuring temperature distribution
JPH07218353A (en) Method and apparatus for measuring temperature distribution by OTDR
CN116499606A (en) An optical fiber Raman temperature measurement device based on detection power correction
JP3377067B2 (en) Brillouin frequency shift distribution measuring method and apparatus
JP2001124664A (en) Measuring device of optical wavelength dispersion and measuring method of optical wavelength dispersion
US5090812A (en) Ring laser cavity backscatter measurement
JPH04332835A (en) Corrective processing method of distributed temperature data
JP2009281813A (en) Optical fiber measuring device
JPH04248426A (en) Apparatus for measuring strain distribution of optical fiber
JP2007240294A (en) Optical fiber strain measurement device
JP3486223B2 (en) Distance measuring device
JP4803083B2 (en) Optical fiber measuring apparatus and sampling method
JP2885980B2 (en) Temperature distribution detector
JPH02306131A (en) Optical fiber temperature distribution measurement method
JPH04318432A (en) Distribution temperature measurement method using optical fiber sensor

Legal Events

Date Code Title Description
FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20070908

Year of fee payment: 7

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20080908

Year of fee payment: 8

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20090908

Year of fee payment: 9

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20090908

Year of fee payment: 9

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100908

Year of fee payment: 10

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100908

Year of fee payment: 10

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110908

Year of fee payment: 11

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110908

Year of fee payment: 11

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120908

Year of fee payment: 12

EXPY Cancellation because of completion of term