JP3502329B2 - 光ファイバひずみ計測方法およびその装置 - Google Patents
光ファイバひずみ計測方法およびその装置Info
- Publication number
- JP3502329B2 JP3502329B2 JP2000123290A JP2000123290A JP3502329B2 JP 3502329 B2 JP3502329 B2 JP 3502329B2 JP 2000123290 A JP2000123290 A JP 2000123290A JP 2000123290 A JP2000123290 A JP 2000123290A JP 3502329 B2 JP3502329 B2 JP 3502329B2
- Authority
- JP
- Japan
- Prior art keywords
- light
- optical fiber
- power
- measurement
- frequency
- 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.)
- Expired - Lifetime
Links
Landscapes
- Length Measuring Devices By Optical Means (AREA)
- Testing Of Optical Devices Or Fibers (AREA)
Description
【０００１】
【発明の属する技術分野】本発明は、光ファイバをセン
サとして用い、その長さ方向についてひずみを連続的に
計測する光ファイバひずみ計測方法およびその装置に関
し、より詳細には、コンクリートや鉄鋼構造物、地盤な
どの計測対象にセンシング用の光ファイバを固定して、
計測対象に生じているひずみを固定されているセンシン
グ用光ファイバで計測するのに利用される光ファイバひ
ずみ計測方法およびその装置に関する。
サとして用い、その長さ方向についてひずみを連続的に
計測する光ファイバひずみ計測方法およびその装置に関
し、より詳細には、コンクリートや鉄鋼構造物、地盤な
どの計測対象にセンシング用の光ファイバを固定して、
計測対象に生じているひずみを固定されているセンシン
グ用光ファイバで計測するのに利用される光ファイバひ
ずみ計測方法およびその装置に関する。
【０００２】
【従来の技術】現在、ひずみの計測に一般的に用いられ
ているものとしてひずみゲージがある。このひずみゲー
ジは、空間分解能が数ｍｍから数ｃｍ程度のゲージ長と
同程度と高く、簡単に局所的なひずみを高精度に計測す
ることができるという利点がある。しかしながら、計測
はゲージが取り付けられた離散点であり、連続した計測
には多くのゲージを密に取り付ける必要がある。また、
信号線を兼ねた電力線を介してゲージに給電する必要が
あり、多点計測ではその電力線の取り扱いが煩雑であ
る。さらに落雷などの電磁雑音の影響を受けるなどの欠
点がある。それに対して光ファイバひずみ計測方法で
は、光ファイバに沿って連続的に数ｋｍから数１０ｋｍ
にわたる長距離計測が可能であるとともに、上述したよ
うな電磁雑音の影響を受けないという利点がある。
ているものとしてひずみゲージがある。このひずみゲー
ジは、空間分解能が数ｍｍから数ｃｍ程度のゲージ長と
同程度と高く、簡単に局所的なひずみを高精度に計測す
ることができるという利点がある。しかしながら、計測
はゲージが取り付けられた離散点であり、連続した計測
には多くのゲージを密に取り付ける必要がある。また、
信号線を兼ねた電力線を介してゲージに給電する必要が
あり、多点計測ではその電力線の取り扱いが煩雑であ
る。さらに落雷などの電磁雑音の影響を受けるなどの欠
点がある。それに対して光ファイバひずみ計測方法で
は、光ファイバに沿って連続的に数ｋｍから数１０ｋｍ
にわたる長距離計測が可能であるとともに、上述したよ
うな電磁雑音の影響を受けないという利点がある。
【０００３】このように、ひずみゲージは注目すべきひ
ずみの発生位置が既知あるいは予想可能であり、かつ計
測範囲が狭い場合に適し、光ファイバひずみ計測技術は
ひずみ発生位置が未知で広範囲の計測が要求される場合
に適している。このような特長から、光ファイバひずみ
計測技術はトンネルなどの大型構造物や堤防などの土木
構造物への適用が期待されている。
ずみの発生位置が既知あるいは予想可能であり、かつ計
測範囲が狭い場合に適し、光ファイバひずみ計測技術は
ひずみ発生位置が未知で広範囲の計測が要求される場合
に適している。このような特長から、光ファイバひずみ
計測技術はトンネルなどの大型構造物や堤防などの土木
構造物への適用が期待されている。
【０００４】図１は、従来の光ファイバひずみ計測装置
の一例として、光ファイバ歪測定装置（特願平８−２４
３７６０号「光ファイバ歪測定装置」）を示す図であ
る。この光ファイバひずみ計測装置において、１は光
源、２は光分岐器、３は光周波数シフタ、４はパルス化
装置、５は光方向性結合器、６は光合波器、７は光検出
器、８は信号処理部、１０〜１７は光ファイバ、１８は
信号線である。
の一例として、光ファイバ歪測定装置（特願平８−２４
３７６０号「光ファイバ歪測定装置」）を示す図であ
る。この光ファイバひずみ計測装置において、１は光
源、２は光分岐器、３は光周波数シフタ、４はパルス化
装置、５は光方向性結合器、６は光合波器、７は光検出
器、８は信号処理部、１０〜１７は光ファイバ、１８は
信号線である。
【０００５】光源１は、一定の光周波数の連続光を発生
するもので、単一波長の連続光を出射するものである。
光分岐器２は、発生された連続光を信号光と参照光とに
分岐するもので、入射端が光源１と光ファイバ１１によ
って接続され、光源１から出射された連続光を２つの出
射端に強度比１対１で出射するものである。光分岐器２
から出射される一方の分岐光を信号光と呼び、もう一方
の分岐光を参照光と呼ぶことにする。
するもので、単一波長の連続光を出射するものである。
光分岐器２は、発生された連続光を信号光と参照光とに
分岐するもので、入射端が光源１と光ファイバ１１によ
って接続され、光源１から出射された連続光を２つの出
射端に強度比１対１で出射するものである。光分岐器２
から出射される一方の分岐光を信号光と呼び、もう一方
の分岐光を参照光と呼ぶことにする。
【０００６】光周波数シフタ３は、光ファイバ内で発生
する微弱なブリルアン散乱光のパワーを高感度に検出す
るために、後述する光検出器７でコヒーレント検波を行
うために挿入されているものであり、その入射端は光分
岐器２の一方の出射端と光ファイバ１２で接続されてい
る。この光周波数シフタ３は、入射される信号光の光周
波数を後述のブリルアン周波数シフト分だけシフトさせ
る。パルス化装置４は、連続した信号光をパルス光に変
換するもので、その入射端と光周波数シフタ３の出射端
とは光ファイバ１３で接続されている。このパルス化装
置４は、入射光を時間領域において電界が矩形波になる
振幅変調を行い、時間幅１０ｎｓ〜１μｓ程度のパルス
光に変換する。
する微弱なブリルアン散乱光のパワーを高感度に検出す
るために、後述する光検出器７でコヒーレント検波を行
うために挿入されているものであり、その入射端は光分
岐器２の一方の出射端と光ファイバ１２で接続されてい
る。この光周波数シフタ３は、入射される信号光の光周
波数を後述のブリルアン周波数シフト分だけシフトさせ
る。パルス化装置４は、連続した信号光をパルス光に変
換するもので、その入射端と光周波数シフタ３の出射端
とは光ファイバ１３で接続されている。このパルス化装
置４は、入射光を時間領域において電界が矩形波になる
振幅変調を行い、時間幅１０ｎｓ〜１μｓ程度のパルス
光に変換する。
【０００７】光方向性結合器５は、パルス光をセンシン
グ用光ファイバ１０に入射するとともに光ファイバ１０
で発生したブリルアン散乱光を出射するものであり、入
射端、入出射端、出射端を有し、入射端から入射された
光を入出射端から出射し、入出射端から入射された光を
出射端から出射する。光方向性結合器５の入射端は、パ
ルス化装置４と光ファイバ１４で接続され、パルス光が
入射される。光方向性結合器５から出射された光は、セ
ンシング用光ファイバ１０に入射される。この出射され
た光により、センシング用光ファイバ１０内で生じた後
方散乱光がこの入出射端に入射する。
グ用光ファイバ１０に入射するとともに光ファイバ１０
で発生したブリルアン散乱光を出射するものであり、入
射端、入出射端、出射端を有し、入射端から入射された
光を入出射端から出射し、入出射端から入射された光を
出射端から出射する。光方向性結合器５の入射端は、パ
ルス化装置４と光ファイバ１４で接続され、パルス光が
入射される。光方向性結合器５から出射された光は、セ
ンシング用光ファイバ１０に入射される。この出射され
た光により、センシング用光ファイバ１０内で生じた後
方散乱光がこの入出射端に入射する。
【０００８】光合波器６は、出射されたブリルアン散乱
光と信号光とを合波するもので、２つの入射端と１つの
出射端を有する。光合波器６の一方の入射端と前述した
光分岐器２の出射端とは光ファイバ１５によって接続さ
れ、もう一方の入射端と光方向性結合器５の出射端とが
光ファイバ１６で接続されている。これにより光合波器
６にはセンシング用光ファイバ１０から出射される後方
散乱光と、光分岐器２から出射される参照光とが入射さ
れ、ここでこれらは合波される。光検出器７は、合波光
を検出して電気信号に変換しブリルアン散乱光のパワー
を計測してそのスペクトルを得る手段であり、入射光に
対してコヒーレント検波を行い、この入射光パワーを電
力に変換して出力する光検出器である。光検出器７の入
射端は光合波器６の出射端と光ファイバ１７で接続され
ている。
光と信号光とを合波するもので、２つの入射端と１つの
出射端を有する。光合波器６の一方の入射端と前述した
光分岐器２の出射端とは光ファイバ１５によって接続さ
れ、もう一方の入射端と光方向性結合器５の出射端とが
光ファイバ１６で接続されている。これにより光合波器
６にはセンシング用光ファイバ１０から出射される後方
散乱光と、光分岐器２から出射される参照光とが入射さ
れ、ここでこれらは合波される。光検出器７は、合波光
を検出して電気信号に変換しブリルアン散乱光のパワー
を計測してそのスペクトルを得る手段であり、入射光に
対してコヒーレント検波を行い、この入射光パワーを電
力に変換して出力する光検出器である。光検出器７の入
射端は光合波器６の出射端と光ファイバ１７で接続され
ている。
【０００９】信号処理部８は、得られた散乱光パワース
ペクトルに対して演算処理を行うためのもので、検出さ
れた散乱光パワーに演算処理を行い、センシング用光フ
ァイバ１０に生じたひずみを求める。
ペクトルに対して演算処理を行うためのもので、検出さ
れた散乱光パワーに演算処理を行い、センシング用光フ
ァイバ１０に生じたひずみを求める。
【００１０】上述した構成において、光源１から出射さ
れた連続光は光分岐器２に入射され、信号光と参照光と
に分岐される。この信号光は光周波数シフタ３に入射さ
れ、光周波数がシフトされる。光周波数シフタ３より出
射された光はパルス化装置４に入射され、時間幅が１０
ｎｓ〜１μｓ程度のパルス光に変換される。パルス化装
置４から出射されたパルス光は、光方向性結合器５を通
過してセンシング用光ファイバ１０に入射される。
れた連続光は光分岐器２に入射され、信号光と参照光と
に分岐される。この信号光は光周波数シフタ３に入射さ
れ、光周波数がシフトされる。光周波数シフタ３より出
射された光はパルス化装置４に入射され、時間幅が１０
ｎｓ〜１μｓ程度のパルス光に変換される。パルス化装
置４から出射されたパルス光は、光方向性結合器５を通
過してセンシング用光ファイバ１０に入射される。
【００１１】パルス光がセンシング用光ファイバ１０に
入射されると、このセンシング用光ファイバ１０中でレ
イリー散乱やブリルアン散乱を受け、後方散乱光が生じ
る。この後方散乱光は光方向性結合器５と光ファイバ１
６を介して光合波器６の一方の入射端に入射される。
入射されると、このセンシング用光ファイバ１０中でレ
イリー散乱やブリルアン散乱を受け、後方散乱光が生じ
る。この後方散乱光は光方向性結合器５と光ファイバ１
６を介して光合波器６の一方の入射端に入射される。
【００１２】前述した光分岐器２から出射された参照光
は、光ファイバ１５を介して光合波器６へ入射され、前
述した後方散乱光と合波される。光合波器６から出射さ
れる合波光は、光検出器７に入射されてコヒーレント検
波され、受信されたパワーに対応した検出信号が光検出
器７から信号処理部８に出力される。
は、光ファイバ１５を介して光合波器６へ入射され、前
述した後方散乱光と合波される。光合波器６から出射さ
れる合波光は、光検出器７に入射されてコヒーレント検
波され、受信されたパワーに対応した検出信号が光検出
器７から信号処理部８に出力される。
【００１３】次に、ブリルアン散乱光の特性とひずみと
の関係について説明し、その後で従来の演算方法につい
て述べる。ブリルアン散乱光は、物質に入射した光が物
質中を伝搬する際に、物質に屈折率の周期的変化を起こ
しながら散乱されて入射端に戻る光であり、上述した装
置を用いると自然ブリルアン後方散乱光のパワーＰ
B（ｚ，ν）は次のように与えられる。
の関係について説明し、その後で従来の演算方法につい
て述べる。ブリルアン散乱光は、物質に入射した光が物
質中を伝搬する際に、物質に屈折率の周期的変化を起こ
しながら散乱されて入射端に戻る光であり、上述した装
置を用いると自然ブリルアン後方散乱光のパワーＰ
B（ｚ，ν）は次のように与えられる。
【００１４】
【数１】
【００１５】
【数２】
【００１６】
【数３】
【００１７】
【数４】
【００１８】ここで、ｚは光ファイバに沿ったパルス光
入射端からの距離、νは自然ブリルアン後方散乱光の光
周波数、ｃは真空中での光速、ｎ0 は光ファイバの屈折
率、Ｐは入射パルス光の全パワー、αZ は光ファイバの
減衰係数である。ｇ（ν，ν B ）はブリルアン利得スペ
クトルであり、式（２）で表わされるローレンツ関数で
与えられる。νB はｇ（ν，νB ）がピークパワーｈと
なるときの光周波数であり、ｗはｇ（ν，νB ）の半値
全幅である。ｐ12 ,λ，ρ，ｖA はそれぞれ光ファイバ
の光弾性係数、入射光の波長、光ファイバの密度、光フ
ァイバ中での音速である。ｔはパルス光を入射してから
その散乱光を検出するまでの時間である。式（１）の最
後の項は、光ファイバ内での減衰のために散乱光パワー
が散乱位置ｚに依存して指数関数的に減衰していくこと
を示している。式（２）において、ｇ（ν，νB ）の形
状はｚに依存しないと仮定している。
入射端からの距離、νは自然ブリルアン後方散乱光の光
周波数、ｃは真空中での光速、ｎ0 は光ファイバの屈折
率、Ｐは入射パルス光の全パワー、αZ は光ファイバの
減衰係数である。ｇ（ν，ν B ）はブリルアン利得スペ
クトルであり、式（２）で表わされるローレンツ関数で
与えられる。νB はｇ（ν，νB ）がピークパワーｈと
なるときの光周波数であり、ｗはｇ（ν，νB ）の半値
全幅である。ｐ12 ,λ，ρ，ｖA はそれぞれ光ファイバ
の光弾性係数、入射光の波長、光ファイバの密度、光フ
ァイバ中での音速である。ｔはパルス光を入射してから
その散乱光を検出するまでの時間である。式（１）の最
後の項は、光ファイバ内での減衰のために散乱光パワー
が散乱位置ｚに依存して指数関数的に減衰していくこと
を示している。式（２）において、ｇ（ν，νB ）の形
状はｚに依存しないと仮定している。
【００１９】光周波数がｆ0 のパルス光を入射し、それ
によってピーク周波数がνB のブリルアン散乱光が生じ
たとする。この差ｓB （＝ｆ0 −νB ）はブリルアル周
波数シフトと呼ばれ、次式で与えられる。
によってピーク周波数がνB のブリルアン散乱光が生じ
たとする。この差ｓB （＝ｆ0 −νB ）はブリルアル周
波数シフトと呼ばれ、次式で与えられる。
【００２０】
【数５】
【００２１】ただしｖA は
【００２２】
【数６】
【００２３】で与えられる。ここで、Ｅはヤング率、κ
はポアソン比である。光ファイバにひずみが発生する
と、式（６）の関係にしたがってｖA が変化し、その結
果式（５）のｓB も変化する。したがって、ｆ0 を一定
に保つと、ひずみの変化に応じてνB が変化することに
なる。このνB の変化は、光ファイバに作用している応
力によって生じるひずみの大きさと比例関係があること
が見出されている。（T.Horiguchi,T.Kurashima,and M.
Tateda "Tensile strain dependence of Brillouinfreq
uency shift in silica optical fibers," IEEE Photon
ics Technol.Lett.,Vol.1, No.5,pp.107-108,1989.5）
そこで、予め、ひずみεの変化Δεと、ブリルアン散乱
光パワースペクトルのピーク周波数νB の変化ΔνB と
の関係を求めておくことにより、得られたνB の値から
ひずみεを求めることができる。ひずみの計測点、すな
わちセンシング用光ファイバ１０の入射端からの距離ｚ
は、パルス光を光ファイバに入射してからブリルアン散
乱光を検出するまでの時間ｔからｚ＝ｃｔ／（２ｎ0 ）
で求められる。また空間分解能Δｚは、光ファイバに入
射するパルス光時間をτとすると、Δｚ＝ｃτ／（２ｎ
0 ）で与えられる。
はポアソン比である。光ファイバにひずみが発生する
と、式（６）の関係にしたがってｖA が変化し、その結
果式（５）のｓB も変化する。したがって、ｆ0 を一定
に保つと、ひずみの変化に応じてνB が変化することに
なる。このνB の変化は、光ファイバに作用している応
力によって生じるひずみの大きさと比例関係があること
が見出されている。（T.Horiguchi,T.Kurashima,and M.
Tateda "Tensile strain dependence of Brillouinfreq
uency shift in silica optical fibers," IEEE Photon
ics Technol.Lett.,Vol.1, No.5,pp.107-108,1989.5）
そこで、予め、ひずみεの変化Δεと、ブリルアン散乱
光パワースペクトルのピーク周波数νB の変化ΔνB と
の関係を求めておくことにより、得られたνB の値から
ひずみεを求めることができる。ひずみの計測点、すな
わちセンシング用光ファイバ１０の入射端からの距離ｚ
は、パルス光を光ファイバに入射してからブリルアン散
乱光を検出するまでの時間ｔからｚ＝ｃｔ／（２ｎ0 ）
で求められる。また空間分解能Δｚは、光ファイバに入
射するパルス光時間をτとすると、Δｚ＝ｃτ／（２ｎ
0 ）で与えられる。
【００２４】このようにブリルアン散乱光を利用したひ
ずみ計測技術では、ひずみ散乱光のパワースペクトルの
ピーク周波数νB から間接的に求められるので、ひずみ
計測精度を向上させるためには、ピーク周波数νB を高
精度に決定しなければならない。
ずみ計測技術では、ひずみ散乱光のパワースペクトルの
ピーク周波数νB から間接的に求められるので、ひずみ
計測精度を向上させるためには、ピーク周波数νB を高
精度に決定しなければならない。
【００２５】以下、図２を参照しながら、具体的にこれ
までに提案されたνB の算出方法について説明する。
までに提案されたνB の算出方法について説明する。
【００２６】従来の方法の第１ステップは、センシング
用光ファイバに沿った各計測点について、光周波数を変
数として解析を行うための準備をするステップであり、
各計測点について、光周波数を変数として、計測された
ブリルアン散乱光のパワーをソートするステップであ
る。前述したように従来の装置では、光周波数を固定し
てパルス光を入射し、そのブリルアン散乱光のパワーを
計測する。そして次に光周波数をある所定量だけ変化さ
せて再び光周波数を固定し、散乱光パワーを計測するこ
とを繰り返す。この処理によって、各光周波数について
計測点を変数としたパワーが得られるので、これをもと
に各計測点毎に、光周波数を変数としたパワーデータを
得る。以下ではこのステップをソートとよぶ。より具体
的には、ｉ番目の計測点までの距離をｚi （ｉ＝１〜
Ｉ，Ｉは全計測点数）、ｊ番目の光周波数をνj （ｊ＝
１〜Ｊ，Ｊは全計測周波数点数）、その時のパワーの計
測値を
用光ファイバに沿った各計測点について、光周波数を変
数として解析を行うための準備をするステップであり、
各計測点について、光周波数を変数として、計測された
ブリルアン散乱光のパワーをソートするステップであ
る。前述したように従来の装置では、光周波数を固定し
てパルス光を入射し、そのブリルアン散乱光のパワーを
計測する。そして次に光周波数をある所定量だけ変化さ
せて再び光周波数を固定し、散乱光パワーを計測するこ
とを繰り返す。この処理によって、各光周波数について
計測点を変数としたパワーが得られるので、これをもと
に各計測点毎に、光周波数を変数としたパワーデータを
得る。以下ではこのステップをソートとよぶ。より具体
的には、ｉ番目の計測点までの距離をｚi （ｉ＝１〜
Ｉ，Ｉは全計測点数）、ｊ番目の光周波数をνj （ｊ＝
１〜Ｊ，Ｊは全計測周波数点数）、その時のパワーの計
測値を
【００２７】
【外１】
【００２８】とすると、まずνj を所定の光周波数値に
セットし、ｚi を変化させて
セットし、ｚi を変化させて
【００２９】
【外２】
【００３０】を計測し、ｚ1 までの計測が終わったとこ
ろでνj を次の光周波数へと変化させνJ まで計測を行
っていく。本ステップでは、あるｚi を選択・固定し、
νj について
ろでνj を次の光周波数へと変化させνJ まで計測を行
っていく。本ステップでは、あるｚi を選択・固定し、
νj について
【００３１】
【外３】
【００３２】をソートするステップである。その結果得
られた
られた
【００３３】
【外４】
【００３４】は、同一のｚi を有しνj を変数としてソ
ースされたものであること、また各ｚ i について独立に
処理がなされることから、表記を簡単化するために以下
ではこれを
ースされたものであること、また各ｚ i について独立に
処理がなされることから、表記を簡単化するために以下
ではこれを
【００３５】
【外５】
【００３６】と書くこととする。
【００３７】第２ステップは、第１ステップにおいて得
られたパワー
られたパワー
【００３８】
【外６】
【００３９】に対して、前述のブリルアン散乱光のパワ
ースペクトルを与えるピーク周波数ν B を算出するステ
ップである。式（２）はνについての非線形関数である
ため、解析解は得られない。そこで、以下のように線形
化して解く方法が提案されている。（C N Pannell,J Dh
liwayo,D J Webb,"How to estimate the accuracy of a
Brillouin distributed temperature sensor,"Proc.OF
S'97,pp.524-527,1997）この方法では、式（２）の逆数
をとり、評価関数Ｅr として
ースペクトルを与えるピーク周波数ν B を算出するステ
ップである。式（２）はνについての非線形関数である
ため、解析解は得られない。そこで、以下のように線形
化して解く方法が提案されている。（C N Pannell,J Dh
liwayo,D J Webb,"How to estimate the accuracy of a
Brillouin distributed temperature sensor,"Proc.OF
S'97,pp.524-527,1997）この方法では、式（２）の逆数
をとり、評価関数Ｅr として
【００４０】
【数７】
【００４１】を考える。ここで、係数ａ1 ，ａ2 ，ａ3
はそれぞれ、
はそれぞれ、
【００４２】
【数８】
【００４３】
【数９】
【００４４】
【数１０】
【００４５】であり、Ｗj は逆数をとったことによる変
化を避けるための重み、σj はパワーの計測値に含まれ
るノイズの標準偏差である。また計測値の全体に対して
あてはめ計算を行うのではなく、通常はあるしきい値を
設定しそれ以上の値を有するピーク周波数近傍に対して
あてはめ計算を行うので、その計算に用いる部分をあら
ためてｊ（＝１〜Ｊ）で表わす。簡単のためにσj が一
定とすると、（９），（１０）よりνB は、
化を避けるための重み、σj はパワーの計測値に含まれ
るノイズの標準偏差である。また計測値の全体に対して
あてはめ計算を行うのではなく、通常はあるしきい値を
設定しそれ以上の値を有するピーク周波数近傍に対して
あてはめ計算を行うので、その計算に用いる部分をあら
ためてｊ（＝１〜Ｊ）で表わす。簡単のためにσj が一
定とすると、（９），（１０）よりνB は、
【００４６】
【数１１】
と求められる。式（７）を最小化する条件から、最小２
乗法を用いてａ2 ，ａ3を求め、この結果を式（１１）
に代入することによりνB を算出する。
乗法を用いてａ2 ，ａ3を求め、この結果を式（１１）
に代入することによりνB を算出する。
【００４７】第３ステップは光ファイバに沿った求める
べき全ての計測点について第２ステップを繰り返すステ
ップである。すなわち、ｚi 変化させながら
べき全ての計測点について第２ステップを繰り返すステ
ップである。すなわち、ｚi 変化させながら
【００４８】
【外７】
【００４９】から
【００５０】
【外８】
【００５１】を得て、ｚ1 まで第２ステップを繰り返す
ステップである。
ステップである。
【００５２】
【発明が解決しようとする課題】しかしながら、以上説
明した方法では、入射パルス光の時間幅が十分に長く、
散乱光のパワースペクトルが式（２）で表されるローレ
ンツ関数で与えられる場合が想定されている。しかし、
空間分解能を上げるためにパルス光の時間幅を短くして
いくと、散乱光のパワースペクトルはローレンツ関数と
は異なり、パワーが広範囲に分布するようになることが
実験的に確認されている。（A.Fellay,L.Thevenaz,M.Fa
cchini,M.Nikles,and P.A.Robert,“Distributed sensi
ng using stimulated Brillouin scattering:towards u
ltimate resolution,”in Proc.OFS'97,pp324-327,199
7）また、このパワースペクトルの広がりは理論的にも
示されており、その結果としてパルス幅が１０ｎｓ程度
以下になると、ピーク周波数νB の決定精度すなわちひ
ずみ計測精度がパルス幅に逆比例して急激に劣化するこ
とも明らかにされている。（H.Naruse,and M.Tateda,
“Trade-off between the spatial and frequency reso
lution in measuring the power spectrum of the Bril
louin backscattered light in an optical fiber,”Ap
pl.Opt.,vol.38,no.31,pp6516-6521,1999.）
明した方法では、入射パルス光の時間幅が十分に長く、
散乱光のパワースペクトルが式（２）で表されるローレ
ンツ関数で与えられる場合が想定されている。しかし、
空間分解能を上げるためにパルス光の時間幅を短くして
いくと、散乱光のパワースペクトルはローレンツ関数と
は異なり、パワーが広範囲に分布するようになることが
実験的に確認されている。（A.Fellay,L.Thevenaz,M.Fa
cchini,M.Nikles,and P.A.Robert,“Distributed sensi
ng using stimulated Brillouin scattering:towards u
ltimate resolution,”in Proc.OFS'97,pp324-327,199
7）また、このパワースペクトルの広がりは理論的にも
示されており、その結果としてパルス幅が１０ｎｓ程度
以下になると、ピーク周波数νB の決定精度すなわちひ
ずみ計測精度がパルス幅に逆比例して急激に劣化するこ
とも明らかにされている。（H.Naruse,and M.Tateda,
“Trade-off between the spatial and frequency reso
lution in measuring the power spectrum of the Bril
louin backscattered light in an optical fiber,”Ap
pl.Opt.,vol.38,no.31,pp6516-6521,1999.）
【００５３】本発明は、このような問題に鑑みてなされ
たもので、その目的とするところは、ひずみ計測に最適
な入射パルス光形状を提供するとともに、その入射パル
ス光を用いて計測されたブリルアン散乱光のパワースペ
クトルに対して、そのピークパワー周波数を精度良く決
定することのできる光ファイバひずみ計測方法およびそ
の装置を提供することにある。
たもので、その目的とするところは、ひずみ計測に最適
な入射パルス光形状を提供するとともに、その入射パル
ス光を用いて計測されたブリルアン散乱光のパワースペ
クトルに対して、そのピークパワー周波数を精度良く決
定することのできる光ファイバひずみ計測方法およびそ
の装置を提供することにある。
【００５４】
【課題を解決するための手段】本発明は、このような目
的を達成するために、請求項１に記載の発明は、一定の
光周波数の連続光を発生する発生ステップと、発生され
た連続光を信号光と参照光に分岐する分岐ステップと、
連続した信号光をパルス光に変換する変換ステップと、
パルス光を光ファイバに入射するとともに光ファイバで
発生したブリルアン散乱光を出射する出射ステップと、
出射されたブリルアン散乱光と信号光とを合波する合波
ステップと、合波光を検出して電気信号に変換し、ブリ
ルアン散乱光のパワーを計測してそのスペクトルを得る
取得ステップと、得られたパワースペクトルからそのピ
ーク周波数を演算処理して求める演算ステップとからな
り、光ファイバに発生している長さ方向のひずみを求め
る光ファイバひずみ計測方法において、前記演算ステッ
プが、光ファイバに沿った長さ方向の各計測点につい
て、光周波数を変数としてブリルアン散乱光のパワーの
計測値をソートする第１のステップと、計測されたパワ
ーに含まれる計測誤差を、決定すべきパワースペクトル
のピーク周波数、ピーク値、半値全幅を変数に含んだ形
で記述し、全計測周波数についての計測誤差の２乗和が
最小になる条件からピーク周波数を決定する第２のステ
ップと、求めるべき長さ方向の全ての計測点について前
記ピーク周波数を繰り返し決定する第３のステップとを
有することを特徴とするものである。
的を達成するために、請求項１に記載の発明は、一定の
光周波数の連続光を発生する発生ステップと、発生され
た連続光を信号光と参照光に分岐する分岐ステップと、
連続した信号光をパルス光に変換する変換ステップと、
パルス光を光ファイバに入射するとともに光ファイバで
発生したブリルアン散乱光を出射する出射ステップと、
出射されたブリルアン散乱光と信号光とを合波する合波
ステップと、合波光を検出して電気信号に変換し、ブリ
ルアン散乱光のパワーを計測してそのスペクトルを得る
取得ステップと、得られたパワースペクトルからそのピ
ーク周波数を演算処理して求める演算ステップとからな
り、光ファイバに発生している長さ方向のひずみを求め
る光ファイバひずみ計測方法において、前記演算ステッ
プが、光ファイバに沿った長さ方向の各計測点につい
て、光周波数を変数としてブリルアン散乱光のパワーの
計測値をソートする第１のステップと、計測されたパワ
ーに含まれる計測誤差を、決定すべきパワースペクトル
のピーク周波数、ピーク値、半値全幅を変数に含んだ形
で記述し、全計測周波数についての計測誤差の２乗和が
最小になる条件からピーク周波数を決定する第２のステ
ップと、求めるべき長さ方向の全ての計測点について前
記ピーク周波数を繰り返し決定する第３のステップとを
有することを特徴とするものである。
【００５５】また、請求項２に記載の発明は、一定の光
周波数の連続光を発生する発生手段と、発生された連続
光を信号光と参照光に分岐する分岐手段と、連続した信
号光をパルス光に変換する変換手段と、パルス光を光フ
ァイバに入射するとともに光ファイバで発生したブリル
アン散乱光を出射する光結合手段と、出射されたブリル
アン散乱光と信号光とを合波する合波手段と、合波光を
検出して電気信号に変換しブリルアン散乱光のパワーを
計測してそのスペクトルを得る取得手段と、得られたパ
ワースペクトルからそのピーク周波数を演算処理して求
める演算手段とからなり、光ファイバに発生している長
さ方向のひずみを求める光ファイバひずみ計測装置にお
いて、前記演算手段が、光ファイバに沿った長さ方向の
各計測点について、光周波数を変数としてブリルアン散
乱光のパワーの計測値をソートする第１の手段と、計測
されたパワーに含まれる計測誤差を、決定すべきパワー
スペクトルのピーク周波数、ピーク値、半値全幅を変数
に含んだ形で記述し、全計測周波数についての計測誤差
の２乗和が最小になる条件からピーク周波数を決定する
第２の手段と、求めるべき長さ方向の全ての計測点につ
いて前記ピーク周波数を繰り返し決定する第３手段とを
有することを特徴とするものである。
周波数の連続光を発生する発生手段と、発生された連続
光を信号光と参照光に分岐する分岐手段と、連続した信
号光をパルス光に変換する変換手段と、パルス光を光フ
ァイバに入射するとともに光ファイバで発生したブリル
アン散乱光を出射する光結合手段と、出射されたブリル
アン散乱光と信号光とを合波する合波手段と、合波光を
検出して電気信号に変換しブリルアン散乱光のパワーを
計測してそのスペクトルを得る取得手段と、得られたパ
ワースペクトルからそのピーク周波数を演算処理して求
める演算手段とからなり、光ファイバに発生している長
さ方向のひずみを求める光ファイバひずみ計測装置にお
いて、前記演算手段が、光ファイバに沿った長さ方向の
各計測点について、光周波数を変数としてブリルアン散
乱光のパワーの計測値をソートする第１の手段と、計測
されたパワーに含まれる計測誤差を、決定すべきパワー
スペクトルのピーク周波数、ピーク値、半値全幅を変数
に含んだ形で記述し、全計測周波数についての計測誤差
の２乗和が最小になる条件からピーク周波数を決定する
第２の手段と、求めるべき長さ方向の全ての計測点につ
いて前記ピーク周波数を繰り返し決定する第３手段とを
有することを特徴とするものである。
【００５６】
【発明の実施の形態】まず、ひずみ計測に最適な入射パ
ルス光形状を理論的に求め、次に計測されたパワーを、
その入射パルス光によって生じる自然ブリルアン後方散
乱光のパワースペクトルにあてはめる方法について説明
する。
ルス光形状を理論的に求め、次に計測されたパワーを、
その入射パルス光によって生じる自然ブリルアン後方散
乱光のパワースペクトルにあてはめる方法について説明
する。
【００５７】時間領域において、図３に示すように電界
包絡線ＥL （ｔ）が変化する入射パルス光を考える。こ
の場合、ＥL ( ｔ）は以下のようにモデル化される。
包絡線ＥL （ｔ）が変化する入射パルス光を考える。こ
の場合、ＥL ( ｔ）は以下のようにモデル化される。
【００５８】（１）ＥL (ｔ）はｔ＝０に関して対称で
ある。すなわち、ＥL (ｔ）＝ＥL（−ｔ）である。 （２）パルスの電界をその最大値で正規化して表す。す
なわち、ＥL (０）＝１、ＥL （−∞）＝ＥL （∞）＝
０とする。 （３）パルス幅をその電界の半値全幅τで定義する。す
なわち、ＥL (−τ／２）＝ＥL (τ／２）＝１／２で
ある。 （４）Ｅ（ｔ）が０から１に立上がったり、１から０に
立下がるのに要する時間、すなわち立上がり／下がり幅
をΔτとする。 （５）立上がりの場合には、０≦ＥL （ｔ）≦１／２を
開始区間、１／２≦Ｅ L （ｔ）≦１を終了区間と呼び、
逆に立下がりの場合には、１≧ＥL （ｔ）≧１／２を開
始区間、１／２≧ＥL （ｔ）≧０を終了区間と呼ぶ。両
区間での電界は(−τ／２，１／２）または (τ／２，
１／２）に関して点対称とする。 （６）開始と終了区間の急岐さをｍ次関数（曲率パラメ
ータｍ）で与える。 （７）立上がり開始／下がり終了区間では上に凸でな
い、立上がり終了／下がり開始区間では下に凸でない、
すなわちｍ≧１とする。ｍ＝１は電界が直線的に立上が
る／下がる波形、ｍ＝∞はステップ的に立上がる／下が
る矩形波を与える。 （８）入射光は図３に示した振幅変調のみを受け、周波
数変調は受けていない。 （９）パルスの持続時間中、パルスは一定の光周波数ｆ
0 である。
ある。すなわち、ＥL (ｔ）＝ＥL（−ｔ）である。 （２）パルスの電界をその最大値で正規化して表す。す
なわち、ＥL (０）＝１、ＥL （−∞）＝ＥL （∞）＝
０とする。 （３）パルス幅をその電界の半値全幅τで定義する。す
なわち、ＥL (−τ／２）＝ＥL (τ／２）＝１／２で
ある。 （４）Ｅ（ｔ）が０から１に立上がったり、１から０に
立下がるのに要する時間、すなわち立上がり／下がり幅
をΔτとする。 （５）立上がりの場合には、０≦ＥL （ｔ）≦１／２を
開始区間、１／２≦Ｅ L （ｔ）≦１を終了区間と呼び、
逆に立下がりの場合には、１≧ＥL （ｔ）≧１／２を開
始区間、１／２≧ＥL （ｔ）≧０を終了区間と呼ぶ。両
区間での電界は(−τ／２，１／２）または (τ／２，
１／２）に関して点対称とする。 （６）開始と終了区間の急岐さをｍ次関数（曲率パラメ
ータｍ）で与える。 （７）立上がり開始／下がり終了区間では上に凸でな
い、立上がり終了／下がり開始区間では下に凸でない、
すなわちｍ≧１とする。ｍ＝１は電界が直線的に立上が
る／下がる波形、ｍ＝∞はステップ的に立上がる／下が
る矩形波を与える。 （８）入射光は図３に示した振幅変調のみを受け、周波
数変調は受けていない。 （９）パルスの持続時間中、パルスは一定の光周波数ｆ
0 である。
【００５９】この場合、入射パルス光の電界Ｅ（ｔ）
は、
は、
【００６０】
【数１２】
【００６１】で与えられる。ここで、ｉは虚数単位であ
る。Ｔ1 ，Ｔ2 はそれぞれ立下がり開始、終了時刻であ
り、
る。Ｔ1 ，Ｔ2 はそれぞれ立下がり開始、終了時刻であ
り、
【００６２】
【数１３】
【００６３】
【数１４】
【００６４】である。このパルス光のパワースペクトル
ＰP （ｆ，ｆ0 ）は、
ＰP （ｆ，ｆ0 ）は、
【００６５】
【数１５】
【００６６】で与えられる。ここで、ｆは入射光につい
ての光周波数を表わす変数である。一例として、τが
１，１０，１００ｎｓの場合について、曲率パラメータ
ｍを変化させてＰP （ｆ，ｆ0 ）を算出した結果を図４
に示す。図４の横軸βとパラメータｒは、それぞれ以下
の式（１６），（１７）で定義される正規化入射光周波
数と立上がり／下がり幅である。
ての光周波数を表わす変数である。一例として、τが
１，１０，１００ｎｓの場合について、曲率パラメータ
ｍを変化させてＰP （ｆ，ｆ0 ）を算出した結果を図４
に示す。図４の横軸βとパラメータｒは、それぞれ以下
の式（１６），（１７）で定義される正規化入射光周波
数と立上がり／下がり幅である。
【００６７】
【数１６】
【００６８】
【数１７】
【００６９】ここで、０≦ｒ≦１であり、ｒ＝０が矩形
波である。ＰP （ｆ，ｆ0 ）の半値全幅はτに逆比例す
るため、ＰP （ｆ，ｆ0 ）はτが十分大きい場合にはｆ
0 のまわりの狭いスペクトルとなり、逆にτが小さくな
ると広い周波数に分布するようになる。またその分布
は、τの値が同じ場合には、ｍが大きくなるほど、すな
わち立上がり／下がり幅が急峻になるほど広がり、矩形
波のパワースペクトルで最も広がっている。図４より、
τが１０ｎｓ程度以上になると、急峻さのパワースペク
トルに与える影響は無視できるようになることがわか
る。
波である。ＰP （ｆ，ｆ0 ）の半値全幅はτに逆比例す
るため、ＰP （ｆ，ｆ0 ）はτが十分大きい場合にはｆ
0 のまわりの狭いスペクトルとなり、逆にτが小さくな
ると広い周波数に分布するようになる。またその分布
は、τの値が同じ場合には、ｍが大きくなるほど、すな
わち立上がり／下がり幅が急峻になるほど広がり、矩形
波のパワースペクトルで最も広がっている。図４より、
τが１０ｎｓ程度以上になると、急峻さのパワースペク
トルに与える影響は無視できるようになることがわか
る。
【００７０】ここで以下の２つのことを仮定して、自然
ブリルアン後方散乱光のパワースペクトルを求める。
ブリルアン後方散乱光のパワースペクトルを求める。
【００７１】（１）ポンプ光の消耗がない、すなわち入
射パルス光から発生した自然ブリルアン後方散乱光への
エネルギーの伝達量が無視できる。 （２）光ファイバ中の多くの散乱点で発生したブリルア
ン後方散乱光の電界には互いに位相相関がなく、このた
め後方散乱光のパワーに関して重ね合せの原理が成り立
つ。別な言い方をすると、他の場所で生じた後方散乱光
によるブリルアン利得への影響は無視でき、このブリル
アン利得によってスペクトルが狭くなる効果が無視でき
る。
射パルス光から発生した自然ブリルアン後方散乱光への
エネルギーの伝達量が無視できる。 （２）光ファイバ中の多くの散乱点で発生したブリルア
ン後方散乱光の電界には互いに位相相関がなく、このた
め後方散乱光のパワーに関して重ね合せの原理が成り立
つ。別な言い方をすると、他の場所で生じた後方散乱光
によるブリルアン利得への影響は無視でき、このブリル
アン利得によってスペクトルが狭くなる効果が無視でき
る。
【００７２】入射パルス光とブリルアン後方散乱光のパ
ワースペクトルの関係を図５に示す。入射光周波数ｆに
対応するブリルアン後方散乱光は、中心周波数ｆ−ｓB
にもつローレンツ関数になる。ｓB は前述したブリルア
ン周波数シフトであり、ｆ0とνB の差である。上述の
仮定により、自然ブリルアン後方散乱光の周波数に依存
する因子Ｈ（ν）は、
ワースペクトルの関係を図５に示す。入射光周波数ｆに
対応するブリルアン後方散乱光は、中心周波数ｆ−ｓB
にもつローレンツ関数になる。ｓB は前述したブリルア
ン周波数シフトであり、ｆ0とνB の差である。上述の
仮定により、自然ブリルアン後方散乱光の周波数に依存
する因子Ｈ（ν）は、
【００７３】
【数１８】
【００７４】で与えられる。ｒが０．５と１の場合につ
いて、τを１，１０，１００ｎｓとして、ｍを変化させ
ながら式（１８）より算出した自然ブリルアン後方散乱
光のパワースペクトルを図６に示す。いずれのτとｒに
ついても、ｍ＝１が最も狭いスペクトルを与え、ｍが大
きくなるにつれて広がったスペクトルを与える。この計
算において、式（２）のローレンツ関数の半値全幅ｗの
値には、実験で得られた８１．４ＭＨｚを用いた。図６
の横軸αと縦軸は、それぞれ、式（１９）で定義される
正規化周波数と、パワースペクトルのピーク値Ｈ（０）
で正規化した値である。
いて、τを１，１０，１００ｎｓとして、ｍを変化させ
ながら式（１８）より算出した自然ブリルアン後方散乱
光のパワースペクトルを図６に示す。いずれのτとｒに
ついても、ｍ＝１が最も狭いスペクトルを与え、ｍが大
きくなるにつれて広がったスペクトルを与える。この計
算において、式（２）のローレンツ関数の半値全幅ｗの
値には、実験で得られた８１．４ＭＨｚを用いた。図６
の横軸αと縦軸は、それぞれ、式（１９）で定義される
正規化周波数と、パワースペクトルのピーク値Ｈ（０）
で正規化した値である。
【００７５】
【数１９】
【００７６】図７は、ｒが０．５と１の場合について、
τとｍを変化させて得られたブリルアン散乱光のパワー
スペクトルの広がりを求めたものである。図７の縦軸
は、図６より数値的に算出した半値全幅Ｗを上述のｗで
正規化した値である。図７では、τ≦１０ｎｓのときは
１ｎｓごとに、τ≧１０ｎｓのときは１０ｎｓごとにＷ
／ｗの値を算出し、それらを直線で結んで示している。
この結果より、Ｗ／ｗの値は、τ＞１００ｎｓではｍの
値にかかわらずほぼ１に収束しており、τが略１０ｎｓ
までは１．５程度まで緩やかに変化している。それに対
し、τ＜１０ｎｓではτの減少とともに急激に大きくな
っていく。また、ｍが１に近づくにつれて、同一空間分
解能を与えるパルス光であってもその散乱光パワースペ
クトルの半値全幅は小さくなり、この効果はｒが１に近
づくほど大きくなることがわかる。
τとｍを変化させて得られたブリルアン散乱光のパワー
スペクトルの広がりを求めたものである。図７の縦軸
は、図６より数値的に算出した半値全幅Ｗを上述のｗで
正規化した値である。図７では、τ≦１０ｎｓのときは
１ｎｓごとに、τ≧１０ｎｓのときは１０ｎｓごとにＷ
／ｗの値を算出し、それらを直線で結んで示している。
この結果より、Ｗ／ｗの値は、τ＞１００ｎｓではｍの
値にかかわらずほぼ１に収束しており、τが略１０ｎｓ
までは１．５程度まで緩やかに変化している。それに対
し、τ＜１０ｎｓではτの減少とともに急激に大きくな
っていく。また、ｍが１に近づくにつれて、同一空間分
解能を与えるパルス光であってもその散乱光パワースペ
クトルの半値全幅は小さくなり、この効果はｒが１に近
づくほど大きくなることがわかる。
【００７７】上述の散乱光パワースペクトルの解析に基
づいて、観測されたブリルアン後方散乱光に重畳してい
るノイズによって生じるピークパワー周波数計測誤差に
ついて見積もる。光電変換によって得られる電気信号パ
ワーは、受信光パワーの２乗に比例した値であるので、
受信光のパワースペクトルの最大値Ｈ（０）を信号Ｈ
S 、電気信号パワーのノイズの２乗平均値をＨN とする
と、信号対雑音比ＳＮＲは（ＨS ／ＨN ）2 で与えられ
る。
づいて、観測されたブリルアン後方散乱光に重畳してい
るノイズによって生じるピークパワー周波数計測誤差に
ついて見積もる。光電変換によって得られる電気信号パ
ワーは、受信光パワーの２乗に比例した値であるので、
受信光のパワースペクトルの最大値Ｈ（０）を信号Ｈ
S 、電気信号パワーのノイズの２乗平均値をＨN とする
と、信号対雑音比ＳＮＲは（ＨS ／ＨN ）2 で与えられ
る。
【００７８】ピークパワー周波数近傍において、自然ブ
リルアン後方散乱光のパワースペクトルＨ（α）を２次
関数で近似し、Ｈ（α）＋ＨN の最大値がα＝Δα（≠
０）で観測されたとする。この場合、ピークパワー周波
数計測誤差Δα，ＨS ，ＨNの関係は、
リルアン後方散乱光のパワースペクトルＨ（α）を２次
関数で近似し、Ｈ（α）＋ＨN の最大値がα＝Δα（≠
０）で観測されたとする。この場合、ピークパワー周波
数計測誤差Δα，ＨS ，ＨNの関係は、
【００７９】
【数２０】
【００８０】で与えられるため、Δαは、
【００８１】
【数２１】
【００８２】と求められる。ここでｐは、パルス幅τ、
立上がり／下がり幅Δτ、急峻さｍによって決まる定数
である。式（１５）と（１８）に基づいて解析的にｐを
求めるかわりに、図６に示したようにτとｍを変化させ
てパワースペクトルを求め、求められた各曲線について
パワーがピークパワーから１０％減少する周波数を読み
取りｐを数値的に算出した。このｐの値と、ＳＮＲの値
として１０ｄＢを式（２１）に代入してΔαを求めた。
その結果を表１に示す。
立上がり／下がり幅Δτ、急峻さｍによって決まる定数
である。式（１５）と（１８）に基づいて解析的にｐを
求めるかわりに、図６に示したようにτとｍを変化させ
てパワースペクトルを求め、求められた各曲線について
パワーがピークパワーから１０％減少する周波数を読み
取りｐを数値的に算出した。このｐの値と、ＳＮＲの値
として１０ｄＢを式（２１）に代入してΔαを求めた。
その結果を表１に示す。
【００８３】
【表１】
【００８４】矩形波の場合の解析結果（H.Naruse,and
M.Tateda,“Trade-off between thespatial and freque
ncy resolution in measuring the power spectrum of
theBrillouin backscattered light in an optical fib
er,”Appl.Opt.,vol.38,no.31,pp.6516-6521,1999）に
基づいて、τが１，１０，１００ｎｓのときのΔαを算
出すると、それらはそれぞれ１０．５１，１．４５，
０．７８であった。表１より、これらの算出結果と今回
の数値計算結果とは、ほぼ一致していることが確認され
る。τの値にかかわらず、正規化立上がり／下がり幅ｒ
が１に近づくにつれて、また、ｍが１に近づくにつれて
計測誤差が小さくなることがわかる。この結果、ｒとｍ
が１である三角波の場合に計測誤差が最小になることが
予想される。すなわち、三角波が最適な入射パルス光の
形状であると考えられる。
M.Tateda,“Trade-off between thespatial and freque
ncy resolution in measuring the power spectrum of
theBrillouin backscattered light in an optical fib
er,”Appl.Opt.,vol.38,no.31,pp.6516-6521,1999）に
基づいて、τが１，１０，１００ｎｓのときのΔαを算
出すると、それらはそれぞれ１０．５１，１．４５，
０．７８であった。表１より、これらの算出結果と今回
の数値計算結果とは、ほぼ一致していることが確認され
る。τの値にかかわらず、正規化立上がり／下がり幅ｒ
が１に近づくにつれて、また、ｍが１に近づくにつれて
計測誤差が小さくなることがわかる。この結果、ｒとｍ
が１である三角波の場合に計測誤差が最小になることが
予想される。すなわち、三角波が最適な入射パルス光の
形状であると考えられる。
【００８５】これまではｒが０．５と１の場合について
の解析だったので、ｍを１としてｒを変化させた場合に
ついて解析する。この場合、式（１５）を考慮して、式
（１８）の積分を実行すると、パワースペクトルＨ
（α）として
の解析だったので、ｍを１としてｒを変化させた場合に
ついて解析する。この場合、式（１５）を考慮して、式
（１８）の積分を実行すると、パワースペクトルＨ
（α）として
【００８６】
【数２２】
【００８７】が得られるので、これより解析的に正規化
ピークパワー周波数計測誤差Δαは次のように求められ
る。
ピークパワー周波数計測誤差Δαは次のように求められ
る。
【００８８】
【数２３】
【００８９】ただし、Ｘ，Ｙは、
【００９０】
【数２４】
【００９１】
【数２５】
【００９２】
【数２６】
【００９３】である。Ａが大きく連続光とみなせる場合
には、Δαはｒの値によらず、
には、Δαはｒの値によらず、
【００９４】
【数２７】
【００９５】となり、立上がり／下がり幅や急峻さの影
響を受けない。一方、Ａが小さい場合には、
響を受けない。一方、Ａが小さい場合には、
【００９６】
【数２８】
【００９７】と近似でき、ｒ＝０を代入すると矩形波の
解析結果に一致することが確認される。この結果は、Δ
αはτと√（１＋ｒ2 ）に逆比例し、ｒ＝１である三角
波のときピーク周波数計測誤差は最小になり、矩形波の
場合の１／√２になるとことを示している。
解析結果に一致することが確認される。この結果は、Δ
αはτと√（１＋ｒ2 ）に逆比例し、ｒ＝１である三角
波のときピーク周波数計測誤差は最小になり、矩形波の
場合の１／√２になるとことを示している。
【００９８】したがって、空間分解能を向上させるため
にパルス幅を短くすることは、ブリルアン散乱光のパワ
ースペクトルを広がらせることになり、その結果とし
て、ピークパワー周波数決定精度を低下させることにな
る。その低下の程度は、入射光のパルス波形に依存し、
矩形波の場合に最も著しく、パルス幅が１ｎｓの誤差は
１００ｎｓの誤差の１３．３倍になる。一方、三角波の
場合にはパルス幅によらずに劣化の程度は最も小さく、
パルス幅が短い場合にはその誤差は矩形波の誤差の１／
√２倍になる。このことから、三角波による計測が最も
高い精度を与えると結論される。
にパルス幅を短くすることは、ブリルアン散乱光のパワ
ースペクトルを広がらせることになり、その結果とし
て、ピークパワー周波数決定精度を低下させることにな
る。その低下の程度は、入射光のパルス波形に依存し、
矩形波の場合に最も著しく、パルス幅が１ｎｓの誤差は
１００ｎｓの誤差の１３．３倍になる。一方、三角波の
場合にはパルス幅によらずに劣化の程度は最も小さく、
パルス幅が短い場合にはその誤差は矩形波の誤差の１／
√２倍になる。このことから、三角波による計測が最も
高い精度を与えると結論される。
【００９９】以上の解析により、最適な入射パルス光の
形状と、それによって生じるブリルアン後方散乱光のパ
ワースペクトル
形状と、それによって生じるブリルアン後方散乱光のパ
ワースペクトル
【０１００】
【数２９】
【０１０１】が得られた。以下では、本解析で得られた
パワースペクトルに、計測されたパワーをあてはめる方
法について説明する。
パワースペクトルに、計測されたパワーをあてはめる方
法について説明する。
【０１０２】あてはめ方法として、ここでは共役勾配法
を例にとって説明する。共役勾配法は制約なしの代表的
な非線形最適化手法の一つであり、関数ｅ＝ｑ（ｙ）を
最小にするｙを逐次近似法の形で求める手法である。
（茨木俊秀、福島雅夫，“最適化の手法，”共立出版，
１９９６）この方法では第ｋ近似解をｙk とした場合、
これにΔｙk の修正を加えて次の第ｋ＋１近似解ｙk+1
を
を例にとって説明する。共役勾配法は制約なしの代表的
な非線形最適化手法の一つであり、関数ｅ＝ｑ（ｙ）を
最小にするｙを逐次近似法の形で求める手法である。
（茨木俊秀、福島雅夫，“最適化の手法，”共立出版，
１９９６）この方法では第ｋ近似解をｙk とした場合、
これにΔｙk の修正を加えて次の第ｋ＋１近似解ｙk+1
を
【０１０３】
【数３０】
【０１０４】と更新していく。この更新の際、Δｙk を
修正量μk と修正方向ベクトルｄk との積として
修正量μk と修正方向ベクトルｄk との積として
【０１０５】
【数３１】
【０１０６】と与え、これらのμk ，ｄk を次のように
して決定する。ここでは、ｄk として共役方向ベクトル
を用いることとする。この場合、ｅ＝ｑ（ｙ）が２次関
数であるとすると、共役方向ベクトルｄk は、勾配方向
ベクトルｓk すなわち
して決定する。ここでは、ｄk として共役方向ベクトル
を用いることとする。この場合、ｅ＝ｑ（ｙ）が２次関
数であるとすると、共役方向ベクトルｄk は、勾配方向
ベクトルｓk すなわち
【０１０７】
【数３２】
【０１０８】を用いて、
【０１０９】
【数３３】
【０１１０】と表せる。ただし、上式においてｄ1 ＝ｓ
1 である。一方、修正量μk についてはｄk を用いて関
数
1 である。一方、修正量μk についてはｄk を用いて関
数
【０１１１】
【数３４】
【０１１２】と考え、μk ＞０の条件下でφ（μk ）を
最小化する値として求める。このようにして決定された
μk ，ｄk からΔｙk を算出し、さらにこの結果を式
（３０）に代入する一連の繰り返し計算により、近似解
の精度を次々に上げていく。関数ｅの極小点においては
最小化する値として求める。このようにして決定された
μk ，ｄk からΔｙk を算出し、さらにこの結果を式
（３０）に代入する一連の繰り返し計算により、近似解
の精度を次々に上げていく。関数ｅの極小点においては
【０１１３】
【数３５】
【０１１４】となるので、式（３５）の条件を満たすま
で繰り返し計算を行ない、ｅの極小点を数値的に算出す
る。
で繰り返し計算を行ない、ｅの極小点を数値的に算出す
る。
【０１１５】上記の共役勾配法を用いた場合のあてはめ
計算方法を説明する。今、ｊ番目の計測光周波数をνj
で、そこで計測されたパワーを
計算方法を説明する。今、ｊ番目の計測光周波数をνj
で、そこで計測されたパワーを
【外９】
【０１１６】で表わし、評価関数Ｆr として
【０１１７】
【数３６】
【０１１８】を考える。ここでＪは計測に用いた光周波
数の数であるため予め与えられる。上述のｙk ，ｓk は
それぞれ
数の数であるため予め与えられる。上述のｙk ，ｓk は
それぞれ
【０１１９】
【数３７】
【０１２０】
【数３８】
【０１２１】となるので、実際にｒk の成分であるｖ
B ，ｈ，ｗの勾配方向のベクトル∂Ｅr／∂ｖB ，∂Ｅr
／∂ｈ，∂Ｅr ／∂ｗを計算すると、以下のようにな
る。
B ，ｈ，ｗの勾配方向のベクトル∂Ｅr／∂ｖB ，∂Ｅr
／∂ｈ，∂Ｅr ／∂ｗを計算すると、以下のようにな
る。
【０１２２】
【数３９】
【０１２３】
【数４０】
【０１２４】
【数４１】
【０１２５】ただし、
【０１２６】
【数４２】
【０１２７】
【数４３】
【０１２８】
【数４４】
【０１２９】である。
【０１３０】以下、図８を用いて本発明の実施例を具体
的に説明する。装置構成については、従来技術のパルス
化装置４において、連続した信号光を電界が三角波形状
を有するパルス光に変換すること以外は従来のものと同
様である。
的に説明する。装置構成については、従来技術のパルス
化装置４において、連続した信号光を電界が三角波形状
を有するパルス光に変換すること以外は従来のものと同
様である。
【０１３１】本発明の第１ステップは従来の方法の第１
ステップと同じで、センシング用光ファイバに沿った各
計測点について、光周波数を変換として計測されたブリ
ルアン散乱光のパワー
ステップと同じで、センシング用光ファイバに沿った各
計測点について、光周波数を変換として計測されたブリ
ルアン散乱光のパワー
【０１３２】
【外１０】
【０１３３】をソートするステップであり、本ステップ
では、各計測点において、光周波数についてソートした
パワー
では、各計測点において、光周波数についてソートした
パワー
【０１３４】
【外１１】
【０１３５】を得るステップである。
【０１３６】第２ステップは、計測されたパワーに含ま
れる計測誤差を、ピーク周波数νB、ピーク値ｈ、半値
全幅ｗとを変数に含んだ形で、電界が三角波形状を有す
るパルス光を光ファイバに入射した際に生じる自然ブリ
ルアン散乱光のパワースペクトルを記述し、計測値
れる計測誤差を、ピーク周波数νB、ピーク値ｈ、半値
全幅ｗとを変数に含んだ形で、電界が三角波形状を有す
るパルス光を光ファイバに入射した際に生じる自然ブリ
ルアン散乱光のパワースペクトルを記述し、計測値
【０１３７】
【外１２】
【０１３８】に対し、その全計測周波数についての計測
誤差が最小になる条件からνB を得るステップである。
なお、
誤差が最小になる条件からνB を得るステップである。
なお、
【０１３９】
【外１３】
【０１４０】は従来の技術で説明したように、ｊのうち
しきい値以上のパワーを有する部分である。決定すべき
変数νB ，ｈ，ｗを含んだ形で記述したブリルアン散乱
光の周波数分布Ｈ（ν）を表す式（１９），（２６），
（２９）と、全計測周波数について計測されたパワーに
含まれる計測誤差の合計を表す式（３６）に対し、共役
勾配法等を用いてその計測誤差を最小化するνB ，ｈ，
ｗを求める。そして、このνB からひずみを算出する。
しきい値以上のパワーを有する部分である。決定すべき
変数νB ，ｈ，ｗを含んだ形で記述したブリルアン散乱
光の周波数分布Ｈ（ν）を表す式（１９），（２６），
（２９）と、全計測周波数について計測されたパワーに
含まれる計測誤差の合計を表す式（３６）に対し、共役
勾配法等を用いてその計測誤差を最小化するνB ，ｈ，
ｗを求める。そして、このνB からひずみを算出する。
【０１４１】第３ステップは、求めるべき全ての計測点
のデータについて第２ステップを繰り返すステップであ
る。
のデータについて第２ステップを繰り返すステップであ
る。
【０１４２】
【発明の効果】以上説明したように本発明によれば、、
一定の光周波数の連続光を発生する発生ステップと、発
生された連続光を信号光と参照光に分岐する分岐ステッ
プと、連続した信号光をパルス光に変換する変換ステッ
プと、パルス光を光ファイバに入射するとともに光ファ
イバで発生したブリルアン散乱光を出射する出射ステッ
プと、出射されたブリルアン散乱光と信号光とを合波す
る合波ステップと、合波光を検出して電気信号に変換
し、ブリルアン散乱光のパワーを計測してそのスペクト
ルを得る取得ステップと、得られたパワースペクトルか
らそのピーク周波数を演算処理して求める演算ステップ
とからなり、光ファイバに発生している長さ方向のひず
みを求める光ファイバひずみ計測方法において、演算ス
テップが、光ファイバに沿った長さ方向の各計測点につ
いて、光周波数を変数としてブリルアン散乱光のパワー
の計測値をソートする第１のステップと、計測されたパ
ワーに含まれる計測誤差を、決定すべきパワースペクト
ルのピーク周波数、ピーク値、半値全幅を変数に含んだ
形で記述し、全計測周波数についての計測誤差の２乗和
が最小になる条件からピーク周波数を決定する第２のス
テップと、求めるべき長さ方向の全ての計測点について
ピーク周波数を繰り返し決定する第３のステップとを有
するので、入射パルス光の形状とそれによって生じる自
然ブリルアン後方散乱光パワースペクトルとの間の関係
を解析し、同一の空間分解能を得るための最適な入射パ
ルス光形状を提供することができるとともに、ひずみ計
測精度を決定することができ、自然ブリルアン後方散乱
光パワーが最大になるピーク周波数を算出するあてはめ
方法を提供することができる。これにより、高空間分解
能計測の際に、従来に比べ、ひずみ計測誤差を１／√２
に低減できる。
一定の光周波数の連続光を発生する発生ステップと、発
生された連続光を信号光と参照光に分岐する分岐ステッ
プと、連続した信号光をパルス光に変換する変換ステッ
プと、パルス光を光ファイバに入射するとともに光ファ
イバで発生したブリルアン散乱光を出射する出射ステッ
プと、出射されたブリルアン散乱光と信号光とを合波す
る合波ステップと、合波光を検出して電気信号に変換
し、ブリルアン散乱光のパワーを計測してそのスペクト
ルを得る取得ステップと、得られたパワースペクトルか
らそのピーク周波数を演算処理して求める演算ステップ
とからなり、光ファイバに発生している長さ方向のひず
みを求める光ファイバひずみ計測方法において、演算ス
テップが、光ファイバに沿った長さ方向の各計測点につ
いて、光周波数を変数としてブリルアン散乱光のパワー
の計測値をソートする第１のステップと、計測されたパ
ワーに含まれる計測誤差を、決定すべきパワースペクト
ルのピーク周波数、ピーク値、半値全幅を変数に含んだ
形で記述し、全計測周波数についての計測誤差の２乗和
が最小になる条件からピーク周波数を決定する第２のス
テップと、求めるべき長さ方向の全ての計測点について
ピーク周波数を繰り返し決定する第３のステップとを有
するので、入射パルス光の形状とそれによって生じる自
然ブリルアン後方散乱光パワースペクトルとの間の関係
を解析し、同一の空間分解能を得るための最適な入射パ
ルス光形状を提供することができるとともに、ひずみ計
測精度を決定することができ、自然ブリルアン後方散乱
光パワーが最大になるピーク周波数を算出するあてはめ
方法を提供することができる。これにより、高空間分解
能計測の際に、従来に比べ、ひずみ計測誤差を１／√２
に低減できる。
【図１】従来の技術による光ファイバを用いたひずみ測
定装置の一例を説明した図である。
定装置の一例を説明した図である。
【図２】計測されたブリルアン散乱光のピーク周波数を
算出するための従来方法を説明するための図である。
算出するための従来方法を説明するための図である。
【図３】入射パルス光の電界包絡線の形状を説明するた
めの図である。
めの図である。
【図４】入射パルス光のパワースペクトルを説明するた
めの図である。
めの図である。
【図５】入射パルス光とそれによって生じるブリルアン
散乱光のパワースペクトルの関係を説明するための図で
ある。
散乱光のパワースペクトルの関係を説明するための図で
ある。
【図６】入射パルス光形状に対するブリルアン後方散乱
光のパワースペクトルを説明するための図である。
光のパワースペクトルを説明するための図である。
【図７】入射パルス光の形状とブリルアン後方散乱光の
正規化半値全幅との関係を説明するための図である。
正規化半値全幅との関係を説明するための図である。
【図８】本発明によってピーク周波数を算出する方法を
説明するための図である。
説明するための図である。
１ 光源
２ 光分岐器
３ 光周波数シフタ
４ パルス化装置
５ 光方向性結合器
６ 光合波器
７ 光検出器
８ 信号処理部
１０〜１７ 光ファイバ
１８ 信号線
─────────────────────────────────────────────────────
フロントページの続き
(58)調査した分野(Int.Cl.7，ＤＢ名)
G01B 11/16
G01L 1/24
G01M 11/00
G01K 11/12
G01D 5/36
Claims (2)
- 【請求項１】 一定の光周波数の連続光を発生する発生
ステップと、発生された連続光を信号光と参照光に分岐
する分岐ステップと、連続した信号光をパルス光に変換
する変換ステップと、パルス光を光ファイバに入射する
とともに光ファイバで発生したブリルアン散乱光を出射
する出射ステップと、出射されたブリルアン散乱光と信
号光とを合波する合波ステップと、合波光を検出して電
気信号に変換し、ブリルアン散乱光のパワーを計測して
そのスペクトルを得る取得ステップと、得られたパワー
スペクトルからそのピーク周波数を演算処理して求める
演算ステップとからなり、光ファイバに発生している長
さ方向のひずみを求める光ファイバひずみ計測方法にお
いて、 前記変換ステップが、時間領域において電界包絡線が直
線的に立上がった後に立下がる三角波振幅変調を行って
連続した信号光をパルス光に変換する変換ステップであ
って、 前記演算ステップが、 光ファイバに沿った長さ方向の各計測点について、光周
波数を変数としてブリルアン散乱光のパワーの計測値を
ソートする第１のステップと、 計測されたパワーに含まれる計測誤差を、決定すべきパ
ワースペクトルのピーク周波数、ピーク値、半値全幅を
変数に含んだ形で記述し、全計測周波数についての計測
誤差の２乗和が最小になる条件からピーク周波数を決定
する第２のステップと、 求めるべき長さ方向の全ての計測点について前記ピーク
周波数を繰り返し決定する第３のステップとを有するこ
とを特徴とする光ファイバひずみ計測方法。 - 【請求項２】 一定の光周波数の連続光を発生する発生
手段と、発生された連続光を信号光と参照光に分岐する
分岐手段と、連続した信号光をパルス光に変換する変換
手段と、パルス光を光ファイバに入射するとともに光フ
ァイバで発生したブリルアン散乱光を出射する光結合手
段と、出射されたブリルアン散乱光と信号光とを合波す
る合波手段と、合波光を検出して電気信号に変換しブリ
ルアン散乱光のパワーを計測してそのスペクトルを得る
取得手段と、得られたパワースペクトルからそのピーク
周波数を演算処理して求める演算手段とからなり、光フ
ァイバに発生している長さ方向のひずみを求める光ファ
イバひずみ計測装置において、 前記変換手段が、時間領域において電界包絡線が直線的
に立上がった後に立下がる三角波振幅変調を行って連続
した信号光をパルス光に変換する変換手段であって、 前記演算手段が、 光ファイバに沿った長さ方向の各計測点について、光周
波数を変数としてブリルアン散乱光のパワーの計測値を
ソートする第１の手段と、 計測されたパワーに含まれる計測誤差を、決定すべきパ
ワースペクトルのピーク周波数、ピーク値、半値全幅を
変数に含んだ形で記述し、全計測周波数についての計測
誤差の２乗和が最小になる条件からピーク周波数を決定
する第２の手段と、 求めるべき長さ方向の全ての計測点について前記ピーク
周波数を繰り返し決定する第３手段とを有することを特
徴とする光ファイバひずみ計測装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2000123290A JP3502329B2 (ja) | 2000-04-24 | 2000-04-24 | 光ファイバひずみ計測方法およびその装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2000123290A JP3502329B2 (ja) | 2000-04-24 | 2000-04-24 | 光ファイバひずみ計測方法およびその装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2001304823A JP2001304823A (ja) | 2001-10-31 |
| JP3502329B2 true JP3502329B2 (ja) | 2004-03-02 |
Family
ID=18633618
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2000123290A Expired - Lifetime JP3502329B2 (ja) | 2000-04-24 | 2000-04-24 | 光ファイバひずみ計測方法およびその装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3502329B2 (ja) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP3930023B2 (ja) * | 2002-11-01 | 2007-06-13 | 欣増 岸田 | 分布型光ファイバセンサシステム |
| JP5002199B2 (ja) * | 2006-06-14 | 2012-08-15 | 富士重工業株式会社 | 接着部はく離形状同定装置 |
| JP4896814B2 (ja) * | 2007-05-18 | 2012-03-14 | ニューブレクス株式会社 | 分布型光ファイバセンサ |
| CN104568249B (zh) * | 2014-12-22 | 2017-02-22 | 天津大学 | 一种基于太赫兹时域光谱系统的应力测量方法 |
| CN105115446B (zh) * | 2015-05-11 | 2018-07-17 | 南昌航空大学 | 基于三角波条纹离焦的条纹反射三维测量方法 |
-
2000
- 2000-04-24 JP JP2000123290A patent/JP3502329B2/ja not_active Expired - Lifetime
Also Published As
| Publication number | Publication date |
|---|---|
| JP2001304823A (ja) | 2001-10-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US9702244B2 (en) | Distributed acoustic sensing gauge length effect mitigation | |
| Niklès | Fibre optic distributed scattering sensing system: Perspectives and challenges for high performance applications | |
| JP4441624B2 (ja) | 光ファイバを用いた歪・温度の分布測定方法及び測定装置 | |
| CN108603773A (zh) | 基于布里渊散射的光电分布式测量装置 | |
| JP2004069685A (ja) | 光ファイバー分散型センサーを使用した測定値の測定及び較正方法 | |
| KR100625807B1 (ko) | 브릴루앙 광섬유 센서를 이용하는 물리량 측정방법 | |
| Lu et al. | Numerical modeling of Fcy OTDR sensing using a refractive index perturbation approach | |
| Kechavarzi et al. | Distributed fibre optic sensing for monitoring reinforced concrete piles | |
| Lu et al. | Spectral properties of the signal in phase-sensitive optical time-domain reflectometry with direct detection | |
| US20230125375A1 (en) | Few-mode rayleigh-based distributed fiber sensor for simultaneous temperature and strain sensing | |
| Handerek et al. | Improved optical power budget in distributed acoustic sensing using enhanced scattering optical fibre | |
| Ravet et al. | Signal processing technique for distributed Brillouin sensing at centimeter spatial resolution | |
| JP3502329B2 (ja) | 光ファイバひずみ計測方法およびその装置 | |
| WO2022044174A1 (ja) | 振動分布測定装置及びその方法 | |
| Liang et al. | Novel fiber Bragg grating sensing method based on the sidelobe modulation for ultrasound detection | |
| RU2428682C1 (ru) | Способ теплового неразрушающего контроля теплотехнического состояния протяженных, сложнопрофильных и труднодоступных объектов | |
| JP3493158B2 (ja) | 光ファイバひずみ計測方法及びその方法を実現するための記録媒体 | |
| Chen et al. | Real-time monitoring of strain processes with large range and high spatial resolution using the method of weak reflection FBG measurement based on OFDR | |
| JP3686588B2 (ja) | 光ファイバひずみ計測方法及びその装置 | |
| CN110823530B (zh) | 一种获取微谐振腔品质因子的方法 | |
| Liu et al. | Optical fiber sensors for landslide monitoring | |
| Iten et al. | Study of a progressive failure in soil using BEDS | |
| JP3154450B2 (ja) | 分布型磁気光ファイバセンサ | |
| RU2248540C1 (ru) | Волоконно-оптический датчик температуры и деформации | |
| Lebang et al. | Enhanced sensitivity of glass optical fiber with various configuration for displacement and force sensor |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20031128 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20031204 |
|
| R151 | Written notification of patent or utility model registration |
Ref document number: 3502329 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R151 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20071212 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20081212 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20091212 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20101212 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20101212 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20111212 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20111212 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20121212 Year of fee payment: 9 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20121212 Year of fee payment: 9 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20131212 Year of fee payment: 10 |
|
| S531 | Written request for registration of change of domicile |
Free format text: JAPANESE INTERMEDIATE CODE: R313531 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| EXPY | Cancellation because of completion of term |