(57)【特許請求の範囲】
【請求項１】
遠端に反射部を有する光線路の前記反射部から反射された第１試験光及び第２試験光を検出し、光線路の特性を解析する光線路特性解析装置であって、
周波数帯域幅が前記光線路のブリルアン利得帯域幅より広く、前記第１試験光と前記第２試験光の周波数帯域幅の差が前記光線路のブリルアン周波数シフトの距離方向の周波数変化量より大きく、中心周波数差が前記ブリルアン周波数シフトと略同一であり、互いに周波数の異なる前記第１試験光及び前記第２試験光をそれぞれ出力する第１の光源及び第２の光源と、
前記第１試験光及び前記第２試験光をそれぞれパルス化して第１試験光パルス及び第２試験光パルスを生成する第１のパルス化部及び第２のパルス化部と、
前記第１のパルス化部及び前記第２のパルス化部に対し、前記第１試験光パルス及び前記第２試験光パルスを発生させるタイミングの時間差を有するように制御する入射時間制御部と、
前記第１試験光パルス及び前記第２試験光パルスを合波し、合波した前記第１試験光パルス及び前記第２試験光パルスを前記光線路に入射する合波素子と、
前記反射部で反射され、前記光線路の入射端から出射された前記第１試験光パルス及び前記第２試験光パルスを戻り光として受光し、それぞれ電気信号へ変換する光受信部と、
前記電気信号をデジタル信号へ変換する変換器と、
前記デジタル信号を用いて前記第１試験光パルスと前記第２試験光パルスとの前記光線路へ入射させる時間差を検出し、
前記光線路で反射された戻り光のうち前記第１試験光パルスの光強度を用いてブリルアン利得を算出し、
前記時間差を繰り返し変化させ、変化させた時間差毎に前記ブリルアン利得を算出し、
前記時間差を繰り返し変化させて算出した前記ブリルアン利得を用いて前記光線路の損失分布を取得する演算処理装置と、
を備える光線路特性解析装置。
【請求項２】
前記光線路は、基幹光線路の一方端を光分岐器によって複数に分岐し、遠端に前記反射部を有する分岐光線路の近端を前記光分岐器の分岐端部それぞれに光結合してなる光線路であり、
前記第１のパルス化部は、
各分岐光線路間の長さの差の最小値に応じて決定したパルス幅でパルス化して前記第１試験光パルスを生成し、
前記演算処理装置は、
前記第１試験光パルスが前記分岐光線路の前記反射部で反射されて戻ってくる時間に応じて、複数の前記分岐光線路のうち前記第１試験光パルスを反射した分岐光線路を特定し、
反射された前記第１試験光パルスの光強度を用いて前記ブリルアン利得を算出し、
前記時間差を繰り返し変化させ、変化させた時間差毎に前記ブリルアン利得を算出し、
前記時間差を繰り返し変化させて算出した前記ブリルアン利得を用いて複数の前記分岐光線路のそれぞれの損失分布を取得する
ことを特徴とする請求項１に記載の光線路特性解析装置。
【請求項３】
前記第１のパルス化部から出力される前記第１試験光パルスのパルス幅τは、前記分岐光線路の前記反射部からの戻り光の時間差２ｎΔＬ／ｃより狭い
ことを特徴とする請求項２に記載の光線路特性解析装置。
ただし、ｎ：前記分岐光線路の屈折率、ΔＬ：前記各分岐光線路間の長さの差の最小値、ｃ：真空中の光速
【請求項４】
前記光受信部及び前記変換器は、
１／τより広い受光可能帯域を有する
ことを特徴とする請求項１から３のいずれかに記載の光線路特性解析装置。
ただし、τ：前記第１試験光パルスのパルス幅
【請求項５】
前記第１の光源及び前記第２の光源は、レーザダイオードであり、
前記第１のパルス化部及び前記第２のパルス化部は、前記レーザダイオードの注入電流を強度変調する手段である
ことを特徴とする請求項１から４のいずれかに記載の光線路特性解析装置。
【請求項６】
遠端に反射部を有する光線路の前記反射部から反射された第１試験光及び第２試験光を検出し、光線路の特性を解析する光線路特性解析方法であって、
周波数帯域幅が前記光線路のブリルアン利得帯域幅より広く、前記第１試験光と前記第２試験光の周波数帯域幅の差が前記光線路のブリルアン周波数シフトの距離方向の周波数変化量より大きく、中心周波数差が前記ブリルアン周波数シフトと略同一であり、互いに周波数の異なる前記第１試験光及び前記第２試験光をそれぞれ出力する第１の光源及び第２の光源による試験光出力手順と、
前記第１試験光及び前記第２試験光をそれぞれパルス化して第１試験光パルス及び第２試験光パルスを生成する第１のパルス化手順及び第２のパルス化手順と、
前記第１のパルス化手順及び前記第２のパルス化手順に対し、前記第１試験光パルス及び前記第２試験光パルスを発生させるタイミングの時間差を有するように制御する入射時間制御手順と、
前記第１試験光パルス及び前記第２試験光パルスを合波し、前記合波した前記第１試験光パルス及び前記第２試験光パルスを前記光線路に入射する合波素子による合波手順と、
前記反射部で反射され、前記光線路の入射端から出射された前記第１試験光パルス及び前記第２試験光パルスを戻り光として受光し、それぞれ電気信号へ変換する光受信手順と、
前記電気信号をデジタル信号へ変換する変換手順と、
前記デジタル信号を用いて前記第１試験光パルスと前記第２試験光パルスとの前記光線路へ入射させる時間差を検出し、
前記光線路で反射された戻り光のうち前記第１試験光パルスの光強度を用いてブリルアン利得を算出し、
前記時間差を繰り返し変化させ、変化させた時間差毎に前記ブリルアン利得を算出し、
前記時間差を繰り返し変化させて算出した前記ブリルアン利得を用いて前記光線路の損失分布を取得する演算処理手順と、
を順に行う光線路特性解析方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、光線路の損失分布を測定するための光線路特性解析装置とその解析方法に関する。
【背景技術】
【０００２】
光ファイバなどの光線路を使用する光通信システムでは、光線路の正常／異常の確認や光線路の故障を検出し、または故障位置を特定する技術が重要である。正常／異常の確認には、接続位置での接続損失値が規格内であるかの評価が行われる。
【０００３】
光線路の接続損失値を評価する方法として、遠端反射ブリルアン利得解析法が提案されている（例えば、非特許文献１参照）。非特許文献１では、パルス化したポンプ光とプローブ光を光線路に入射し、遠端で反射した一方の光と他方の光が光ファイバ中で衝突した際に生じるブリルアン散乱を測定することで接続損失の評価が可能である。この時、ポンプ光とプローブ光の光周波数差は、光ファイバ中のブリルアン周波数シフトと一致させる必要がある。
【０００４】
一方、分岐下部光ファイバの損失情報を個別に測定可能な技術が提案されている（例えば、特許文献１、非特許文献２参照。）。特許文献１では、ポンプ光パルス及びプローブ光パルスの二つの試験光パルスを入射し、両試験光の衝突位置でのブリルアン利得を解析することにより、分岐下部光ファイバ個別の損失分布を測定する。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】国際公開第２０１２／１６５５８７号公報
【非特許文献】
【０００６】
【非特許文献１】Ｈ．  Ｔａｋａｈａｓｈｉ  ｅｔ  ａｌ．，  “Ｃｅｎｔｒａｌｉｚｅｄ  Ｍｅａｓｕｒｅｍｅｎｔ  ｏｆ  Ａｃｔｕａｌ  Ｓｐｌｉｃｅ  Ｌｏｓｓ  ｆｏｒ  Ｉｎｓｔａｌｌｅｄ  Ｏｐｔｉｃａｌ  Ｆｉｂｅｒ  Ｃａｂｌｅ  Ｎｅｔｗｏｒｋｓ  Ｕｓｉｎｇ  Ｅｎｄ−ｒｅｆｌｅｃｔｉｏｎ−ａｓｓｉｓｔｅｄ  Ｂｒｉｌｌｏｕｉｎ  Ａｎａｌｙｓｉｓ，”  Ｐｒｏｃｅｅｄｉｎｇｓ  ｏｆ  ｔｈｅ  ６３ｒｄ  ＩＷＣＳ，  ｐａｐｅｒ  ３−６  （２０１４）．
【非特許文献２】Ｈ．  Ｔａｋａｈａｓｈｉ，  Ｆ．  Ｉｔｏ，  Ｃ．  Ｋｉｔｏ， ａｎｄ  Ｋ．  Ｔｏｇｅ，   “Ｉｎｄｉｖｉｄｕａｌ  ｌｏｓｓ  ｄｉｓｔｒｉｂｕｔｉｏｎ  ｍｅａｓｕｒｅｍｅｎｔ  ｉｎ  ３２−ｂｒａｎｃｈｅｄ  ＰＯＮ  ｕｓｉｎｇ  ｐｕｌｓｅｄ  ｐｕｍｐ−ｐｒｏｂｅ  Ｂｒｉｌｌｏｕｉｎ  Ａｎａｌｙｓｉｓ，”  Ｏｐｔｉｃｓ  Ｅｘｐｒｅｓｓ，  Ｖｏｌ．２１，  Ｎｏ．６，  ６７３９，（２０１３）．
【非特許文献３】Ｈ．  Ｔａｋａｈａｓｈｉ  ｅｔ  ａｌ．，“Ｉｎｄｉｖｉｄｕａｌ  ＰＯＮ  Ｍｏｎｉｔｏｒｉｎｇ  Ｕｓｉｎｇ  Ｍａｉｎｔｅｎａｎｃｅ  Ｂａｎｄ  Ｐｕｌｓｅｄ  Ｐｕｍｐ−Ｐｒｏｂｅ  Ｂｒｉｌｌｏｕｉｎ  Ａｎａｌｙｓｉｓ，”  ＯｐｔｏＥｌｅｃｔｒｏｎｉｃｓ  ａｎｄ  Ｃｏｍｍｕｎｉｃａｔｉｏｎｓ  Ｃｏｎｆｅｒｅｎｃｅ，ＴｈＰ１−４，２０１３．
【非特許文献４】高橋央， “広帯域プローブ光によるブリルアン損失測定法，” 信学会ソサイエティ大会２０１４．
【発明の概要】
【発明が解決しようとする課題】
【０００７】
この手法においても非特許文献１と同様に、光ファイバ中のブリルアン周波数シフトと一致させる必要がある。ブリルアン散乱は光ファイバの温度や歪、またはファイバパラメータの違いにより、ブリルアン周波数シフトは異なる。そのため、既存の光線路では、ブリルアン周波数シフトが長手方向で異なる。
【０００８】
これに対して、ポンプ光パルスとプローブ光パルスとの光周波数差を変化させて測定を繰り返し行うことでブリルアン散乱の強度を取得する方法が提案されている（例えば、非特許文献３参照）。しかしながら、非特許文献３の方法では、光周波数差を変化させて測定を繰り返し行う必要があるため、ポンプ光パルスとプローブ光パルスとの光周波数差を精密に制御しなければならない。このことは、試験装置構成及び測定方法の複雑化を招いていた。
【０００９】
上記を解決するために、非特許文献４では、プローブ光パルスの光周波数帯域幅を広帯域にすることで、装置の複雑化を緩和する方法が提案されているが、プローブ光を広帯域にすることにより、ブリルアン散乱に寄与するプローブ光エネルギーが減少する。
【００１０】
これは、ポンプ光とプローブ光の周波数差がブリルアン周波数シフトと一致した一部分のプローブ光エネルギーしかブリルアン散乱を起こさず、広帯域化したプローブ光の大部分がブリルアン散乱に寄与しなくなるからである。
【００１１】
この広帯域化したプローブ光を非特許文献４では、光受信器で受信した後、Ａ／Ｄ変換器によりデジタル電気信号に変換し、ブリルアン増幅を受けた広帯域プローブ光とブリルアン増幅を受けていない広帯域プローブ光のデジタル電気信号を比較することでブリルアン散乱強度（光ファイバの損失）を求めている。
【００１２】
この場合、ブリルアン散乱に寄与していない部分のプローブ光エネルギーのデジタル電気信号があることで、ブリルアン散乱のデジタル電気信号を検出するためにはより高分解能なＡ／Ｄ変換器が必要になる。また、実現可能な分解能には制限があるため、非特許文献２のように狭い帯域のポンプ光とプローブ光を利用した場合より線路損失測定の感度が劣化してしまい、３２分岐や６４分岐スプリッタ下部または損失の大きい線路を測定できなくなる。
【００１３】
また、非特許文献３では、ポンプ光とプローブ光の周波数差を変化させ、ブリルアン周波数シフトと一致していない測定において入射したプローブ光エネルギーは無駄になり、非特許文献４では、ブリルアン周波数シフトと一致していない周波数のプローブ光エネルギーは無駄になる。
【００１４】
そのため、ブリルアン利得帯域幅（３０−４０ＭＨｚ）に対して、補償すべきブリルアン周波数シフト変化量が大きくなると、（ブリルアン利得帯域幅：Ａ）／（補償すべきブリルアン周波数シフト変化量：Ｂ）の信号しか得られない。
【００１５】
これを補償するためには、測定回数をＢ／Ａ倍にする必要があるが、非特許文献２では、３２分岐スプリッタの下部において損失分布を精度良く測定するためにすでに数万回のオーダで測定を繰り返して平均値を求めるようにしているため、さらにプローブ光エネルギーを補償するためには、必要な測定時間はますます長くなるという課題もある。
【００１６】
前記課題を解決するために、本発明は、光分岐線路において、光スプリッタからユーザ装置側の分岐下部光ファイバ、および装置（光スプリッタなどの光デバイス）を監視するにあたり、新たに光デバイスや光線路構成を変更することなく（既設設備を変更することによりコストをかけることなく）、所外既設設備（同一特性の分岐光ファイバと同一波長光反射フィルタ）を使用するのみで光線路特性を個別に測定することを可能とし、分岐光ファイバの長手方向でブリルアン周波数シフトが異なる場合でも高速・高感度かつ装置構成が簡易な光試験システムを提供することにある。
【課題を解決するための手段】
【００１７】
上記目的を達成するため、本発明では、遠端反射ブリルアン利得解析法において、プローブ光パルスおよびポンプ光パルスの周波数帯域幅をそれぞれブリルアン利得帯域幅より広くし、プローブ光パルスとポンプ光パルスとの周波数帯域幅の差を被測定光線路のブリルアン周波数シフトの距離方向の周波数変化量より大きくし、プローブ光パルスとポンプ光パルスとの中心周波数差をブリルアン周波数シフト程度とする。
【００１８】
具体的には、本発明に係る光線路特性解析装置は、
遠端に反射部を有する光線路の前記反射部から反射された第１試験光及び第２試験光を検出し、光線路の特性を解析する光線路特性解析装置であって、
周波数帯域幅が前記光線路のブリルアン利得帯域幅より広く、前記第１試験光と前記第２試験光の周波数帯域幅の差が前記光線路のブリルアン周波数シフトの距離方向の周波数変化量より大きく、中心周波数差が前記ブリルアン周波数シフトと略同一であり、互いに周波数の異なる前記第１試験光及び前記第２試験光をそれぞれ出力する第１の光源及び第２の光源と、
前記第１試験光及び前記第２試験光をそれぞれパルス化して第１試験光パルス及び第２試験光パルスを生成する第１のパルス化部及び第２のパルス化部と、
前記第１のパルス化部及び前記第２のパルス化部に対し、前記第１試験光パルス及び前記第２試験光パルスを発生させるタイミングの時間差を有するように制御する入射時間制御部と、
前記第１試験光パルス及び前記第２試験光パルスを合波し、合波した前記第１試験光パルス及び前記第２試験光パルスを前記光線路に入射する合波素子と、
前記反射部で反射され、前記光線路の入射端から出射された前記第１試験光パルス及び前記第２試験光パルスを戻り光として受光し、それぞれ電気信号へ変換する光受信部と、
前記電気信号をデジタル信号へ変換する変換器と、
前記デジタル信号を用いて前記第１試験光パルスと前記第２試験光パルスとの前記光線路へ入射させる時間差を検出し、
前記光線路で反射された戻り光のうち前記第１試験光パルスの光強度を用いてブリルアン利得を算出し、
前記時間差を繰り返し変化させ、変化させた時間差毎に前記ブリルアン利得を算出し、
前記時間差を繰り返し変化させて算出した前記ブリルアン利得を用いて前記光線路の損失分布を取得する演算処理装置と、を備える。
【００１９】
本発明に係る光線路特性解析装置では、
前記光線路は、基幹光線路の一方端を光分岐器によって複数に分岐し、遠端に前記反射部を有する分岐光線路の近端を前記光分岐器の分岐端部それぞれに光結合してなる光線路であり、
前記第１のパルス化部は、
各分岐光線路間の長さの差の最小値に応じて決定したパルス幅でパルス化して前記第１試験光パルスを生成し、
前記演算処理装置は、
前記第１試験光パルスが前記分岐光線路の前記反射部で反射されて戻ってくる時間に応じて、複数の前記分岐光線路のうち前記第１試験光パルスを反射した分岐光線路を特定し、
反射された前記第１試験光パルスの光強度を用いて前記ブリルアン利得を算出し、
前記時間差を繰り返し変化させ、変化させた時間差毎に前記ブリルアン利得を算出し、
前記時間差を繰り返し変化させて算出した前記ブリルアン利得を用いて複数の前記分岐光線路のそれぞれの損失分布を取得してもよい。
【００２０】
本発明に係る光線路特性解析装置では、
前記第１のパルス化部から出力される前記第１試験光パルスのパルス幅τは、前記分岐光線路の前記反射部からの戻り光の時間差２ｎΔＬ／ｃより狭いことが望ましい。
ただし、ｎ：前記分岐光線路の屈折率、ΔＬ：前記各分岐光線路間の長さの差の最小値、ｃ：真空中の光速である。
【００２１】
本発明に係る光線路特性解析装置では、
前記光受信部及び前記変換器は、
１／τより広い受光可能帯域を有することが望ましい。
ただし、τ：前記第１試験光パルスのパルス幅である。
【００２２】
本発明に係る光線路特性解析装置では、
前記第１の光源及び前記第２の光源は、レーザダイオードであり、
前記第１のパルス化部及び前記第２のパルス化部は、前記レーザダイオードの注入電流を強度変調する手段であってもよい。
【００２３】
具体的には、本発明に係る光線路特性解析方法は、
遠端に反射部を有する光線路の前記反射部から反射された第１試験光及び第２試験光を検出し、光線路の特性を解析する光線路特性解析方法であって、
周波数帯域幅が前記光線路のブリルアン利得帯域幅より広く、前記第１試験光と前記第２試験光の周波数帯域幅の差が前記光線路のブリルアン周波数シフトの距離方向の周波数変化量より大きく、中心周波数差が前記ブリルアン周波数シフトと略同一であり、互いに周波数の異なる前記第１試験光及び前記第２試験光をそれぞれ出力する第１の光源及び第２の光源による試験光出力手順と、
前記第１試験光及び前記第２試験光をそれぞれパルス化して第１試験光パルス及び第２試験光パルスを生成する第１のパルス化手順及び第２のパルス化手順と、
前記第１のパルス化手順及び前記第２のパルス化手順に対し、前記第１試験光パルス及び前記第２試験光パルスを発生させるタイミングの時間差を有するように制御する入射時間制御手順と、
前記第１試験光パルス及び前記第２試験光パルスを合波し、前記合波した前記第１試験光パルス及び前記第２試験光パルスを前記光線路に入射する合波素子による合波手順と、
前記反射部で反射され、前記光線路の入射端から出射された前記第１試験光パルス及び前記第２試験光パルスを戻り光として受光し、それぞれ電気信号へ変換する光受信手順と、
前記電気信号をデジタル信号へ変換する変換手順と、
前記デジタル信号を用いて前記第１試験光パルスと前記第２試験光パルスとの前記光線路へ入射させる時間差を検出し、
前記光線路で反射された戻り光のうち前記第１試験光パルスの光強度を用いてブリルアン利得を算出し、
前記時間差を繰り返し変化させ、変化させた時間差毎に前記ブリルアン利得を算出し、
前記時間差を繰り返し変化させて算出した前記ブリルアン利得を用いて前記光線路の損失分布を取得する演算処理手順と、を順に行う。
【００２４】
なお、上記各発明は、可能な限り組み合わせることができる。
【発明の効果】
【００２５】
本発明によれば、光分岐線路において、光スプリッタからユーザ装置側の分岐下部光ファイバ、および装置（光スプリッタなどの光デバイス）を監視するにあたり、新たに光デバイスや光線路構成を変更することなく（既設設備を変更することによりコストをかけることなく）、所外既設設備（同一特性の分岐光ファイバと同一波長光反射フィルタ）を使用するのみで光線路特性を個別に測定することを可能とし、分岐光ファイバの長手方向でブリルアン周波数シフトが異なる場合でも高速・高感度かつ装置構成が簡易な光試験システムを提供することができる。
【００２６】
以上のように、本発明は、光スプリッタの分岐点で複数系統に分岐された各分岐光ファイバの長さの差ΔＬを利用し、波長の異なる二種の試験光を用意し、パルス化した後に二種のパルス試験光に入射時間差を与えて被測定光線路に入射することで、先に入射したパルス試験光（第１試験光）の反射光と、後に入射したパルス試験光（第２試験光）が対向伝搬することにより発生した誘導ブリルアン後方散乱光を光受信器で受信し、光受信器の出力電流を２ｎΔＬ／ｃ（ｃは光速）よりも高い時間分解能で解析する。
【００２７】
このことで、第１〜第Ｎのどの分岐光ファイバからの誘導ブリルアン散乱かを特定することが可能となり、第１試験光と第２試験光の周波数帯域幅を光線路全域のブリルアン周波数シフトのばらつき以上に設定することで、第１試験光と第２試験光の周波数差を制御することなくブリルアン利得を取得できる。
【００２８】
さらに、第１試験光と第２試験光の周波数帯域幅をほぼ同じ値に設定することで、第１試験光と第２試験光のエネルギーを無駄にすることなくブリルアン相互作用を発生させることができ、高感度な測定が可能となる。
【００２９】
したがって、本発明によれば、既存の所外設備を取り替えることなく、ブリルアン周波数シフトの異なる光線路の損失分布を高感度かつ装置構成を簡易に測定することのできる光線路特性測定装置及びその測定方法を提供することができる。
【図面の簡単な説明】
【００３０】
【図１】本実施形態に係る分岐光線路特性解析装置の構成の一例を示す。
【図２】本実施形態に係る光源の装置構成の一例を示す。
【図３】本実施形態に係る光源の装置構成の一例を示す。
【図４】第１試験光と第２試験光の周波数帯域幅の差が被測定光ファイバのブリルアン周波数シフト変化量より大きい場合の一例を示す。
【図５】本実施形態に係る測定手順の一例を示す。
【発明を実施するための形態】
【００３１】
以下、本発明の実施形態について、図面を参照しながら詳細に説明する。なお、本発明は、以下に示す実施形態に限定されるものではない。これらの実施の例は例示に過ぎず、本発明は当業者の知識に基づいて種々の変更、改良を施した形態で実施することができる。なお、本明細書及び図面において符号が同じ構成要素は、相互に同一のものを示すものとする。
【００３２】
（実施形態１）
以下、図面を参照してこの発明に係わる実施形態を説明する。図１は、本発明の一実施形態に係る分岐光線路特性解析装置の構成を示すブロック図である。図１に示される光線路特性解析装置として機能する分岐光線路特性解析装置は、第１試験光出力部１１及び第２試験光出力部１２、光帯域制御部１３、入射時間制御器１５、パルス化部１４、合波素子１６、サーキュレータ１７、光受信部２１、Ａ／Ｄ変換器２２及び演算処理装置２３を具備する。
【００３３】
試験光出力部として機能する第１試験光出力部１１および第２試験光出力部１２は、例えば、図２及び３に示す構成により実現され、連続光を出力する。第１試験光出力部１１および第２試験光出力部１２から出力される連続光を、以下ではそれぞれプローブ光、ポンプ光とする。プローブ光とポンプ光は、ブリルアン周波数シフト程度の異なる光周波数を出力する。
【００３４】
プローブ光及びポンプ光は、光帯域制御部１３により、距離方向の周波数変化量より広い帯域幅に設定される。図２に示されるように、試験光出力部は、レーザダイオード２５に注入する電流に変調を加えることで広帯域な光を出力する。
【００３５】
また、図３で示されるように、試験光出力部は、広帯域光源２６から出力される自然放出光を光フィルタ２７で帯域制限することで広帯域な第１及び第２試験光を出力する。なお、広帯域光源２６からの出力光は、光増幅器（不図示）からの自然放出光で代用しても構わない。
【００３６】
パルス化部１４は、音響光学素子をパルス駆動する音響光学スイッチ（不図示）を備える。なお、パルス化部１４は、電気光学素子をＬｉＮｂＯ３を用いてパルス駆動する導波路スイッチを備えるようにしてもよい。以下では、音響光学スイッチから成る音響光学変調器、又は、電気光学スイッチから成るＬｉＮｂＯ３変調器を光デバイスと称する。
【００３７】
光デバイスは、電気パルスにより駆動される時間で、パルス化部＃１及びパルス化部＃２により、それぞれ第１試験光及び第２試験光をパルス化する。また、パルス化部１４は、図２のレーザダイオード２５を用いる場合、レーザダイオード２５に注入する電流をパルス化することで代用しても構わない。
【００３８】
このとき、パルス化部＃１は、第１試験光パルスのパルス幅を２ｎΔＬ／ｃ以下とする。ここで、ΔＬは、分岐光線路として機能する被測定光ファイバ２０に含まれる各分岐光ファイバの長さの差の最小値である。ｃは、真空中の光速である。ｎは、光ファイバの屈折率である。
【００３９】
入射時間制御部は、パルス化部１４において、駆動用の電気パルスの変調時間を変化させることができる機能を有する。入射時間制御部１５は、パルス化部１４がパルスを変調するタイミングを変化させることで、光パルスが生成されるタイミングを制御する。
【００４０】
合波素子１６は、第１試験光パルスと第２試験光パルスを合波する。合波素子１６により合波された第１試験光パルスと第２試験光パルスとは、サーキュレータ１７を通過して被測定光ファイバ２０へ入射される。これにより、第１試験光パルスを被測定光ファイバ２０へ入射する時間と、第２試験光パルスを被測定光ファイバ２０へ入射する時間とに時間差を与えることが可能となる。
【００４１】
なお、第１試験光のパルス化部＃１における電気パルス、または、第２試験光のパルス化部＃２における電気パルスのうち一方の電気パルスの電圧を常にゼロにすることにより、第１試験光パルス又は第２試験光パルスの一方のみを被測定光ファイバ２０へ入射することも可能となる。
【００４２】
図１に示すように、被測定光ファイバ２０は、基幹光ファイバの一方端に光結合される光分岐器として機能する光スプリッタ１８、及び分岐光線路として機能する分岐光ファイバを備える。ここで、各分岐光ファイバの終端には反射部として機能する反射型光フィルタ１９が配置される。これらの反射型光フィルタ１９は、第１試験光パルス及び第２試験光パルスそれぞれの波長の光を反射し、それ以外の波長の光を透過する特性を有する。なお、前記反射型光フィルタ１９は光線路の開放端のフレネル反射でも構わない。
【００４３】
光スプリッタ１８は、供給される第１試験光パルス及び第２試験光パルスをＮ（例えば、Ｎ＝８）個にそれぞれ分岐する。光スプリッタ１８で分岐された第１試験光パルスと第２試験光パルスとは、分岐光ファイバ中でインタラクションする。第１試験光パルスと第２試験光パルスとのインタラクションにより、誘導ブリルアン散乱の後方散乱光が発生する。誘導ブリルアン後方散乱光、第１試験光パルス、及び第２試験光パルスは、サーキュレータ１７に戻り、サーキュレータ１７を通過して光受信部２１へ導出される。
【００４４】
光受信部２１は、被測定光ファイバ２０からの戻り光であるサーキュレータ１７の導出光を受信する。光受信部２１は、受信した導出光を電流信号へ変換する。Ａ／Ｄ変換器２２は、光受信部２１から出力される電流信号をデジタル信号に変換する。Ａ／Ｄ変換器２２は、デジタル信号を演算処理装置２３へ出力する。
【００４５】
演算処理装置２３は、入力された電流値に対して下記に説明するような演算処理を行い、分岐光ファイバについての距離に対する損失分布を求める。まず、演算処理装置２３は、入射時間差ｔ１で第１試験光パルスと第２試験光パルスとの両方を被測定光ファイバ２０へ入射した場合の第１試験光強度から、第１試験光パルスのみを被測定光ファイバ２０へ入射した場合の第１試験光強度を引くことで、誘導ブリルアン利得を求める。
【００４６】
次に、演算処理装置２３は、ｔ１を変化させて誘導ブリルアン利得の測定を繰り返し、距離に対する誘導ブリルアン利得を算出する。演算処理装置２３は、距離に対する誘導ブリルアン利得から、分岐光ファイバについての距離に対する損失分布を求める。
【００４７】
次に、上述したように構成される本実施形態の分岐光線路特性解析装置の動作について説明する。まず、第１試験光および第２試験光の出力光周波数及び周波数帯域幅、パルス化部１４、光受信部２１、及びＡ／Ｄ変換器２２は次の条件を満足する必要がある。
（条件１）  第１試験光と第２試験光の周波数帯域幅の差は、被測定光ファイバ２０に想定されるブリルアン周波数シフトの距離方向の周波数変化量より大きいこと。
（条件２）  第１試験光の出力光周波数の最小値は、第２試験光の出力光周波数の最小値と比べて、設定した想定されるブリルアン周波数シフトの最大値よりも低い周波数であること。
（条件３）  第１試験光の出力光周波数の最大値は、第２試験光の出力光周波数の最大値と比べて、設定した想定されるブリルアン周波数シフトの最小値よりも高い周波数であること。
（条件４）  パルス化部１４から出力される第１試験光パルスのパルス幅τは、分岐光ファイバの終端の反射型光フィルタ１９からの戻り光の時間差２ｎΔＬ／ｃより狭いこと。
（条件５）  光受信部２１の帯域及びＡ／Ｄ変換器２２の帯域は、パルス幅τを受光可能な帯域であること。
条件１〜５は次のような意味を持つ。
【００４８】
条件１〜３は、第１試験光パルスと第２試験光パルスとが誘導ブリルアン散乱を、被測定光ファイバ２０中の全ての距離で起こすために必要となる条件である。図４に、第１試験光及び第２試験光の周波数帯域幅の差が被測定光ファイバ２０のブリルアン周波数シフトの変化量より大きい場合の概念図を示す。
【００４９】
第１試験光の周波数帯域幅がブリルアン周波数シフトの変化量より大きい場合、分岐光線路特性解析装置は、第１試験光の光周波数を変更することなく、ファイバＡ中で衝突した第２試験光パルスによる誘導ブリルアン散乱と、ファイバＢ中で衝突した第２試験光パルスによる誘導ブリルアン散乱とを測定可能である。
【００５０】
条件４は、分岐光ファイバ毎の誘導ブリルアン散乱光が重ならないようにするための条件である。光パルス化部１４のパルス幅τが各分岐光ファイバの終端の反射型光フィルタ１９からの戻り光の時間差の最小値２ｎΔＬ／ｃより広いとき、分岐光ファイバ毎の誘導ブリルアン散乱光が重なる。これにより、分岐光ファイバ毎の誘導ブリルアン散乱光を時間的に切り分けることができなくなる。
【００５１】
条件５は、パルス幅τの光パルスを正確に測定するための条件である。すなわち、光受信部２１の帯域、及びＡ／Ｄ変換器２２の帯域は、１／τより広い必要があることを意味する。上記の条件を満足する場合の本発明を用いた分岐光線路特性解析方法を説明する。
【００５２】
第１試験光出力部１１及び第２試験光出力部１２としてそれぞれ機能する第１の光源及び第２の光源は、波長の異なる二つの試験光（第１試験光、第２試験光）を出力する。第１試験光はプローブ光であり、第２試験光はポンプ光である。ここで、第１試験光及び第２試験光の光周波数をｆ１±Δｆ１／２、ｆ２±Δｆ２／２とする。Δｆ１、Δｆ２は、第１試験光及び第２試験光の周波数帯域幅である。ｆＢ＿ｍｉｎ及びｆＢ＿ｍａｘは、被測定光ファイバ２０に想定される誘導ブリルアン後方散乱による光周波数シフト量の最小値及び最大値である。
【００５３】
まず、分岐光線路特性解析装置は、プローブ光パルスを被測定光ファイバ２０へ入射させる。そして、分岐光線路特性解析装置は、プローブ光パルスを入射してｔ１秒後に、ポンプ光パルスを被測定光ファイバ２０へ入射する。被測定光ファイバ２０に入射されたプローブ光パルスとポンプ光パルスとは、光スプリッタ１８によりＮ（例えば、Ｎ＝８）個に分岐される。
【００５４】
（ａ）広帯域試験光によるブリルアン利得
プローブ光とポンプ光との周波数差がｆＢである場合、プローブ光パルスとポンプ光パルスとがインタラクションすると、ブリルアン回折格子が発生し、この回折格子によりポンプ光が反射されることで、プローブ光パルスは増幅される。
【００５５】
ここで、ブリルアン格子強度Ｑ、プローブ光の電界振幅をＡｓ、ポンプ光の電界振幅をＡＰ、ブリルアン利得係数をｇＢとすると、ｚ地点でのプローブ光電界は式（１）で表される増幅を受ける。また、ブリルアン格子の減衰定数をΓとすると、ブリルアン格子強度は式（２）で表せる。
【数１】
0006325490.tif 000002
【数２】
0006325490.tif 000003
【００５６】
ここで、ブリルアンフォノンの格子は十分強められた後を考えると、式（３）で表せる。また、式（１）、（３）より、ｚ地点でのプローブ光が受ける増幅は、式（４）で表せる。
【数３】
0006325490.tif 000004
【数４】
0006325490.tif 000005
【００５７】
ここで、同じ帯域幅を持つ広帯域プローブ光と広帯域ポンプ光において、電界振幅の周波数密度をＡｓ＿ｐａｔｉａｌおよびＡｐ＿ｐａｔｉａｌ、全電界振幅をＡｓ＿ｔｏｔａｌおよびＡｐ＿ｔｏｔａｌとすると、式（３）はプローブ光とポンプ光の全電界振幅により増強されるので、式（５）で表せる。
【数５】
0006325490.tif 000006
【００５８】
また、ｚ地点でのある周波数におけるプローブ光電界の変化量は、式（６）で表される。プローブ光の全周波数の変化量の総和は、式（７）で表せ、ｚ地点におけるプローブ光強度の変化量（ブリルアン利得）は、プローブ光強度をＩｐｒｏｂｅ＿ｔｏｔａｌ，ポンプ光強度をＩｐｕｍｐ＿ｔｏｔａｌとすると、式（８）で表せる。
【数６】
0006325490.tif 000007
【数７】
0006325490.tif 000008
【数８】
0006325490.tif 000009
【００５９】
以上のように、本実施形態では、帯域幅が同じ場合の広帯域プローブ光及び広帯域ポンプ光によるブリルアン利得は、入射の全エネルギーをブリルアン増幅に用いることができる。
【００６０】
上述した広帯域のポンプ光とプローブ光のブリルアン増幅において、プローブ光とポンプ光の帯域幅が違う場合は、ポンプ光とプローブ光のブリルアン増幅に作用する帯域幅は、プローブ光とポンプ光の帯域幅の小さい方に制限されることとなる。
【００６１】
上記条件２，３より、ポンプ光帯域幅 Δｆ２が小さい場合を考えると、帯域幅の差｜Δｆ１−Δｆ２｜なので、ブリルアン増幅に寄与しないエネルギーの割合｜Δｆ１−Δｆ２｜／Δｆ２となる。ここで、帯域幅の差は想定されるブリルアン周波数シフトの変化量により設定し、上述した補償すべきブリルアン周波数シフト変化量Ｂと同じである。
【００６２】
つまり、ポンプ光帯域幅をブリルアン利得帯域幅より広く設定することにより、狭い帯域のポンプ光を使った場合よりも、平均測定化後において大きなブリルアン利得を得ることが可能である。
【００６３】
（ｂ）被測定光ファイバの線路損失の測定
プローブ光とポンプ光との周波数差がｆＢである場合、プローブ光パルスとポンプ光パルスとがインタラクションすると、ブリルアン回折格子が発生し、この回折格子によりポンプ光が反射されることで、プローブ光パルスは式（９）で表される増幅を受ける。
【数９】
0006325490.tif 000010
【００６４】
ここで、ｚは、分岐光ファイバの入射端から、プローブ光パルスとポンプ光パルスとがインタラクションした位置までの距離である。αＢ（ｚ，ｆ）は、入射端からの位置ｚ、かつ、ブリルアン周波数シフトｆのときの誘導ブリルアン利得である。
【００６５】
ｇＢ（ｆ）は、ブリルアン周波数シフトｆの場合の誘導ブリルアン散乱係数である。Ｉｐｕｍｐ（ｚ）は、分岐光ファイバの入射端から距離ｚだけ離れた位置におけるポンプ光パルスの強度である。τｐｕｍｐはポンプ光パルスのパルス幅である。
【００６６】
分岐光ファイバ（＃ｉ）の損失係数をαｉ、分岐光ファイバ（＃ｉ）を往復する場合の全損失を２Ｌｉ、ブリルアン周波数シフトをｆ、入射プローブパワーをＩｐｒｏｂｅ（０）とすると、終端の反射型光フィルタ（＃ｉ）で反射された後、入射端から距離ｚ１の位置でポンプ光パルスと衝突したプローブ光パルスの、分岐光ファイバの入射端での強度Ｉｐｒｏｂｅ（２Ｌｉ，ｚ１，ｆ）は、式（１０）で表される。
【数１０】
0006325490.tif 000011
【００６７】
式（１０）より、分岐光ファイバの入射端での第１試験光の強度Ｉｐｒｏｂｅ（２Ｌｉ，ｚ，ｆ）は、ｇＢ（ｆ）とＩｐｕｍｐ（ｚ）との関数となる。ここで、Ｉｐｕｍｐ（ｚ）は、入射ポンプパワーをＩｐｕｍｐ（０）とすると、式（１１）で表される。また、プローブ光のみを入射した場合に分岐光ファイバの入射端へ戻ってくる反射プローブ光強度Ｉｒｅｆ（２Ｌｉ）は、式（１２）で表される。
【数１１】
0006325490.tif 000012
【数１２】
0006325490.tif 000013
【００６８】
よって、式（１０）は、式（１１）及び式（１２）を用いると式（１３）として表される。下記（１３）式より、誘導ブリルアン散乱光の利得は、インタラクションした場所までの損失と誘導ブリルアン散乱係数との積をポンプ光パルス幅で積分した値となる。ここで、誘導ブリルアン利得係数ｇＢ（ｆ）は有限の周波数帯域をもち、周波数方向に積分した値は定数である。
【００６９】
そのため、広帯域プローブ光及び広帯域ポンプ光において、条件２、３を満たす場合、常に誘導ブリルアン利得係数の周波数方向の積分値を取得できる。つまり、上記式（１３−１及び１３−２）は、のみの関数となる。
【００７０】
そこで、以下式（１４）の演算を行うことで、ある地点ｚ０からの損失分布を取得することができる。よって、誘導ブリルアン散乱光の特性を解析すれば、ある地点ｚ０を基準にした被測定光ファイバ２０の線路損失を測定することができる。
【数１３】
0006325490.tif 000014
【数１４】
0006325490.tif 000015
【００７１】
（ｃ）分岐光ファイバの距離に対するブリルアン散乱光分布の測定
被測定光ファイバ２０の入射端から分岐光ファイバ（＃ａ）（１≦ａ≦Ｎの整数、ここではＮ＝８）の終端までの長さをＬａとする。第１試験光パルスは、分岐光ファイバ（＃ａ）の終端に設置された反射型光フィルタ（＃ａ）により反射される。
【００７２】
ここで、分岐光ファイバの終端からの距離をｌ、被測定光ファイバ２０の屈折率をｎ、真空中の光速をｃとすると、反射された第１試験光パルスはｔ１／２秒後にｌ＝ｃ／ｎ×ｔ１／２だけ進むので、被測定光ファイバ２０の入射端からの距離をｌｘ１とすると、その距離ｌｘ１は、式（１５）となる。
【数１５】
0006325490.tif 000016
【００７３】
また、第１試験光パルスが被測定光ファイバ２０に入射されてから光スプリッタ１８により分岐され、分岐光ファイバの終端の反射型光フィルタ（＃ａ）で反射されて、被測定光ファイバ２０の入射端からの距離ｌｘ１に到達する時間ｔは、式（１６）である。
【数１６】
0006325490.tif 000017
【００７４】
第２試験光を被測定光ファイバ２０に入射する時刻は、第１試験光を入射してからｔ１秒後とする。第２試験光がｔ秒後に到達する被測定光ファイバ２０の入射端からの距離をｌｘ２とすると、その距離ｌｘ２は式（１７）で表される。
【数１７】
0006325490.tif 000018
【００７５】
式（１５）、式（１７）より、ｃ・ｔ１／２ｎの位置でプローブ光パルスとポンプ光パルスとはインタラクションする。また、インタラクションする時間は、分岐光ファイバの終端の反射型光フィルタ（＃ａ）で反射された時間からｔ１／２秒後である。
【００７６】
つまり、プローブ光パルスとポンプ光パルスとを被測定光ファイバ２０に入射する時間差ｔ１を変化させることにより、第１試験光と第２試験光がインタラクションする位置を制御できる。このため、距離に対する誘導ブリルアン散乱の特性分布を求めることができる。
【００７７】
（ｄ）分岐光ファイバ（＃ａ）で反射された第１試験光パルスが光受信部２１に到達する時間の測定
第１試験光パルスが光受信部２１に到達する時間をｔｄａとする。第１試験光パルスは、分岐光ファイバの終端の反射型光フィルタ（＃ａ）により反射され、光受信部２１へ戻ってくる。そのため、到達時間は、式（１８）で表される。
【数１８】
0006325490.tif 000019
【００７８】
ここで、他の分岐光ファイバ（＃ｂ）（１≦ｂ≦Ｎの整数、ここではＮ＝８）から戻ってきた第１試験光パルスが光受信部２１に到達する時間ｔｄｂは、式（１９）で表される。よって、光受信部２１に戻る時間差は、式（２０）で表される。
【数１９】
0006325490.tif 000020
【数２０】
0006325490.tif 000021
【００７９】
Ｌａ≠Ｌｂのとき、光受信部２１に到達する時間が異なる。ここで、第１試験光パルスのパルス幅をτとすると、式（２１）のとき（条件４）、分岐光ファイバ（＃１）〜（＃８）から戻った第１試験光パルスは光スプリッタ１８で重ならない。
【数２１】
0006325490.tif 000022
【００８０】
そのため、光受信部２１の到達時間を測定することで、分岐光ファイバ（＃１）〜（＃８）のうちどの分岐光ファイバから出力された第１試験光パルスであるかを時間的に切り分けることができる。
【００８１】
上記（ｂ）〜（ｄ）により、分岐光線路特性解析装置は、測定を繰り返す間で光周波数制御することなく、分岐光ファイバ個別の損失分布を測定可能である。また、上記（ａ）により、広帯域なポンプ光とプローブ光を用いることにより、高速・高感度な測定が可能である。
【００８２】
以上より、本実施形態の分岐光線路特性解析装置は、被測定光ファイバ２０の距離方向において、温度及び歪、並びに、ファイバパラメータの違いによりブリルアン周波数シフトが変化した場合においても、第１試験光と第２試験光との周波数差を変更させずに、高速・高感度な分岐光ファイバそれぞれの損失分布を取得可能である。
【００８３】
図５は、本実施形態の分岐光線路特性解析装置による測定手順を示すフローチャートである。まず、第１試験光パルスと第２試験光パルスとの入射時間差ｔ１を設定する（ステップＳ１０１）。次に、分岐光線路特性解析装置は、第１試験光パルスのみを被測定光ファイバ２０へ出力する（ステップＳ１０２）。
【００８４】
演算処理装置２３は、戻り光の到達時間からどの分岐光ファイバの反射型光フィルタ１９で反射された第１試験光パルスのデジタル信号であるかを特定する（ステップＳ１０３）。また、演算処理装置２３は、反射プローブ光強度を取得する（ステップＳ１０４）。
【００８５】
次に、分岐光線路特性解析装置は、第１試験光パルス及び第２試験光パルスを被測定光ファイバ２０へ出力する（ステップＳ１０５）。演算処理装置２３は、分岐光ファイバの反射型光ファイバにより反射された第１試験光パルスから誘導ブリルアン散乱光を解析する（ステップＳ１０６）。演算処理装置２３は、解析結果を後段へ出力する（ステップＳ１０７）。
【００８６】
続いて、演算処理装置２３は、最長の分岐光ファイバの反射型光ファイバにより反射された第２試験光パルスが到達したか否かを判断する（ステップＳ１０８）。演算処理装置２３は、第２試験光パルスが到達した場合（ステップＳ１０８のＹｅｓ）、入力時間差が２ｎＬ／ｃと等しいか否かを判断する（ステップＳ１０９）。
【００８７】
入力時間差が２ｎＬ／ｃと等しくない場合（ステップＳ１０９のＮｏ）、演算処理装置２３は、入力時間差ｔをｔ＝ｔ１＋Δｔと設定して（ステップＳ１１１）、ステップＳ１０１の処理から解析処理を繰り返し行う。ステップＳ１０９で入力時間差が２ｎＬ／ｃと等しい場合（ステップＳ１０９のＹｅｓ）、演算処理装置２３は、一連の測定作業を終了する。
【００８８】
以上の処理をまとめると、本実施形態に係る解析方法は、以下のような処理手順となる。
手順１：第１試験光パルスと第２試験光パルスとの入射時間差ｔ１を設定。
手順２：第１試験光パルスの戻り時間によりどの分岐光ファイバで反射した第１試験光パルスであるかを特定。
手順３：第１試験光パルスの光強度により、誘導ブリルアン利得を出力。
手順４：第１試験光パルスと第２試験光パルスとの入射時間差ｔ１を変化させて上記手順１から手順３を繰り返し、ｔ＝２ｎＬ／ｃで終了。
【００８９】
このように、本実施形態では、分岐光線路特性解析装置は、周波数帯域幅の差が被測定光ファイバ２０における長手方向のブリルアン周波数シフトの変化量より大きく、周波数差がブリルアン周波数シフトと同程度の第１試験光と第２試験光とを用意する。
【００９０】
分岐光線路特性解析装置は、第１試験光を、第１試験光のパルス幅が分岐光ファイバの長さの差の最小値ΔＬを利用した条件を満たすようにパルス化すると共に、第２試験光をパルス化する。分岐光線路特性解析装置は、第１試験光パルス及び第２試験光パルスに入射時間差を与えて被測定光ファイバ２０へ入射する。
【００９１】
そして、分岐光線路特性解析装置は、入射時間差を変えながら、第１試験光パルス及び第２試験光パルスを被測定光ファイバ２０へ入射し、第１試験光パルス及び第２試験光パルスのインタラクションにより発生する誘導ブリルアン後方散乱光を解析する。
【００９２】
これにより、分岐光線路特性解析装置は、高速または高感度に、第１試験光と第２試験光との周波数差を変更させずに、どの分岐光ファイバからの誘導ブリルアン散乱かを特定することが可能となる。このため、分岐光線路特性解析装置は、第１試験光と第２試験光との周波数差を変更させることなく、分岐光ファイバ毎の損失分布を求めることができる。
【００９３】
したがって、上記実施形態の分岐光線路特性解析装置によれば、ＰＯＮ型の光分岐線路において、光スプリッタ１８からユーザ装置側の分岐光ファイバを監視するにあたり、所外既設設備を使用するのみで分岐光ファイバの光線路特性を個別に測定することが可能であり、かつ、分岐光ファイバの距離方向でブリルアン周波数シフトが異なる場合でも第１試験光パルス及び第２試験光パルスの光周波数差を変化させることなく精確に、高速かつ高感度に分岐光ファイバの損失特性を測定することができる。
【００９４】
なお、上記実施形態では、分岐光線路特性解析装置がサーキュレータ１７を備える場合を例に説明した。しかしながら、これに限定されない。分岐光線路特性解析装置は、例えば、サーキュレータ１７の代わりにカプラを備えていても構わない。
【００９５】
なお、この発明は、上記実施形態そのままに限定されるものではなく、実施段階ではその要旨を逸脱しない範囲で構成要素を変形して具体化できる。また、上記実施形態に開示されている複数の構成要素の適宜な組み合せにより種々の発明を形成できる。例えば、実施形態に示される全構成要素から幾つかの構成要素を削除してもよい。さらに、異なる実施形態に亘る構成要素を適宜組み合せてもよい。
【産業上の利用可能性】
【００９６】
本発明は情報通信産業に適用することができる。
【符号の説明】
【００９７】
１１：第１試験光出力部
１２：第２試験光出力部
１３：光帯域制御部
１４：パルス化部
１５：入射時間制御部
１６：合波素子
１７：サーキュレータ
１８：光スプリッタ
１９：反射型光フィルタ
２０：被測定光ファイバ
２１：光受信部
２２：Ａ／Ｄ変換器
２３：演算処理装置
２５：レーザダイオード
２６：広帯域光源
２７：光フィルタ
【図１】
0006325490.tif 000023
【図２】
0006325490.tif 000024
【図３】
0006325490.tif 000025
【図４】
0006325490.tif 000026
【図５】
0006325490.tif 000027