(57)【要約】【課題】分岐光線路の個別の特性分布を高速で測定可能な分岐光線路特性解析方法および分岐光線路特性解析装置を提供すること、並びに当該分岐光線路特性解析装置で特性分布を高速で測定可能な分岐光線路の構造を特定できる分岐光線路設計方法を提供することを目的とする。【解決手段】本発明に係る分岐光線路特性解析方法は、光スプリッタの上流側から光周波数の等しいＸ個（Ｘは２以上の自然数）のプローブ光パルスを入射し、その後、プローブ光パルスと異なる光周波数の１つのポンプ光パルスを入射し、発生したブリルアン散乱光を受光して解析することとした。また、分岐光線路測定対象の総延長が与えられた際に、測定性能を最大化するための、分岐損失が大きすぎず、ポンプディプレッションが生じにくい最適な分岐数を所定式を最大化する値とすることとした。【選択図】図１
2017026512.tif 000002
【特許請求の範囲】
【請求項１】
光スプリッタに接続される複数の分岐光線路に、前記分岐光線路でブリルアン（Ｂｒｉｌｌｏｕｉｎ）散乱を発生させるポンプ光パルスと、前記分岐光線路の前記光スプリッタと反対側の遠端で反射させて前記ポンプ光パルスと相互作用させるプローブ光パルスと、を伝搬させ、前記分岐光線路の特性を解析する分岐光線路特性解析方法であって、
前記プローブ光パルスと前記ポンプ光パルスとの光周波数差である設定ブリルアン周波数シフトｆＢを設定する光周波数差設定ステップと、
パルス間隔が等しく且つ光周波数の等しいＸ個（Ｘは２以上の自然数）の前記プローブ光パルスを発生し、Ｘ個の前記プローブ光パルスを前記光スプリッタを経由して前記分岐光線路に入射するプローブ光パルス入射ステップと、
前記プローブ光パルス入射ステップ後に、前記ポンプ光パルスを発生し、Ｘ個のプローブ光パルスが入射された後に前記ポンプ光パルスを前記光スプリッタを経由して前記分岐光線路に入射するポンプ光パルス入射ステップと、
前記ポンプ光パルス入射ステップの後、前記分岐光線路からの戻り光パルスを前記光スプリッタを介して受光し、電気信号に変換する戻り光受光ステップと、
前記戻り光受光ステップで変換した前記電気信号から設定ブリルアン周波数シフトｆＢにおける前記戻り光パルスのブリルアン利得スペクトル分布を取得する特性解析ステップと、
を行うことを特徴とする分岐光線路特性解析方法。
【請求項２】
前記光周波数差設定ステップの前に、最長の前記分岐光線路と最短の前記分岐光線路との長さの差ΔＬｍａｘを検出する光路長差検出ステップを行うことを特徴とする請求項１に記載の分岐光線路特性解析方法。
【請求項３】
前記プローブ光パルス入射ステップは、前記分岐光線路それぞれからの前記戻り光パルスの重畳を防止するように前記プローブ光パルスのパルス間隔Ｔｓを設定することを特徴とする請求項１又は２に記載の分岐光線路特性解析方法。
【請求項４】
前記分岐光線路の線長差が等間隔ΔＬである場合、前記プローブ光パルス入射ステップは、τｐｒｏｂｅを前記プローブ光パルスのパルス幅、νを前記分岐光線路中の光速、Ｎを前記分岐光線路の分岐数とすれば、
Ｔｓ＝τｐｒｏｂｅ＋２ΔＬｍａｘ／ν
でパルス間隔Ｔｓを設定することを特徴とする請求項３に記載の分岐光線路特性解析方法。
ただし、２ΔＬ＝τｐｒｏｂｅνを満たし、ΔＬｍａｘ＝ΔＬ（Ｎ−１）を満たす。
【請求項５】
所定の光周波数範囲内で前記設定ブリルアン周波数シフトｆＢを変化させ、前記光周波数差設定ステップから前記特性解析ステップを繰り返すことを特徴とする請求項１から３のいずれかに記載の分岐光線路特性解析方法。
【請求項６】
請求項１から５のいずれかに記載の分岐光線路特性解析方法で特性を解析する分岐光線路の分岐数Ｎを、ブリルアン利得が最大となるように決定する分岐光線路設計方法。
ただし、前記ブリルアン利得は数（Ｃ１）で定義される。
【数Ｃ１】
2017026512.tif 000022
ここで、ｇＢはブリルアン利得係数、Ａｅｆｆは光ファイバのコアの実効断面積、ｌｃは前記プローブ光パルスと前記ポンプ光パルスの相関長、Ｌは分岐下部心線の長さ、αは光ファイバの光損失係数、Ｐｒ（Ｌ）は前記分岐光線路の遠端での前記プローブ光パルスのパワー、Ｐｕ（Ｌ）は前記分岐光線路の遠端での前記ポンプ光パルスのパワーである。
【請求項７】
光スプリッタに接続される複数の分岐光線路に、前記分岐光線路でブリルアン（Ｂｒｉｌｌｏｕｉｎ）散乱を発生させるポンプ光パルスと、前記分岐光線路の前記光スプリッタと反対側の遠端で反射させて前記ポンプ光パルスと相互作用させるプローブ光パルスと、を伝搬させ、前記分岐光線路の特性を解析する分岐光線路特性解析装置であって、
パルス間隔が等しく且つ光周波数の等しいＸ個（Ｘは２以上の自然数）の前記プローブ光パルスを発生し、Ｘ個の前記プローブ光パルスを前記光スプリッタを経由して前記分岐光線路に入射するプローブ光パルス入射部と、
前記ポンプ光パルスを発生し、Ｘ個のプローブ光パルスが入射された後に前記ポンプ光パルスを前記光スプリッタを経由して前記分岐光線路に入射するポンプ光パルス入射部と、
前記分岐光線路からの戻り光パルスを前記光スプリッタを介して受光し、電気信号に変換する戻り光受光部と、
前記戻り光受光部が変換した前記電気信号から設定ブリルアン周波数シフトｆＢにおける前記戻り光パルスのブリルアン利得スペクトル分布を取得する特性解析部と、
前記プローブ光パルスと前記ポンプ光パルスとの光周波数差である前記設定ブリルアン周波数シフトｆＢを設定し、
前記プローブ光パルス入射部に、前記分岐光線路それぞれからの前記戻り光パルスの重畳を防止するように前記プローブ光パルスのパルス間隔Ｔｓを設定させ、
前記プローブ光パルス入射部に前記プローブ光パルスを前記分岐光線路に入射させた後に前記ポンプ光パルス入射部に前記ポンプ光パルスを前記分岐光線路に入射させる制御部と、
を備えることを特徴とする分岐光線路特性解析装置。
【請求項８】
前記制御部は、所定の光周波数範囲内で前記設定ブリルアン周波数シフトｆＢを変化させ、前記プローブ光パルス入射部と前記ポンプ光パルス入射部に前記プローブ光パルスと前記ポンプ光パルスを繰り返し前記分岐光線路に入射させることを特徴とする請求項７に記載の分岐光線路特性解析装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、光ファイバを分岐する光スプリッタの下流側の光ファイバ（以下では分岐光線路と称する）の特性を、光スプリッタの上流側から解析する技術に関する。
【背景技術】
【０００２】
光センシングは構造物の周辺の歪みや温度変化をモニタリングする技術である。非特許文献１に記載のブリルアンＯＴＤＲ（Ｏｐｔｉｃａｌ  Ｔｉｍｅ  Ｄｏｍａｉｎ  Ｒｅｆｌｅｃｔｏｍｅｔｅｒ）や、非特許文献２に記載のＢＯＣＤＡ（Ｂｒｉｌｌｏｕｉｎ  Ｏｐｔｉｃａｌ  Ｃｏｒｒｅｌａｔｉｏｎ  −  ｄｏｍａｉｎ  ａｎａｌｙｓｉｓ）などが知られている。非特許文献３に記載の、ＦＢＧ（Ｆｉｂｅｒ  Ｂｒａｇｇ  Ｇｒａｔｉｎｇ）を配置する方法も提案されている。
【０００３】
ＦＴＴＨ（Ｆｉｂｅｒ  ｔｏ  ｔｈｅ  Ｈｏｍｅ）を志向する光通信分野では、ＰＯＮ（Ｐａｓｓｉｖｅ  Ｏｐｔｉｃａｌ  Ｎｅｔｗｏｒｋ）型ネットワークにおける分岐光線路の特性（損失など）を個別にモニタリングしたいというニーズがある。関連技術が非特許文献４に開示される。非特許文献４に示される技術は、光スプリッタの上流側の光ファイバにプローブ光パルスとポンプ光パルスとを入射し、両光パルスの衝突位置でのブリルアン利得を解析して分岐光線路の特性分布を個別に測定するという技術である。さらに、関連する技術が非特許文献５にも開示される。
【先行技術文献】
【非特許文献】
【０００４】
【非特許文献１】Ｈ．  Ｏｈｎｏ，  Ｈ．  Ｎａｒｕｓｅ，  Ｍ．  Ｋｉｈａｒａ，  ａｎｄ  Ａ．  Ｓｈｉｍａｄａ，“Ｉｎｄｕｓｔｒｉａｌ  Ａｐｐｌｉｃａｔｉｏｎｓ  ｏｆ  ｔｈｅ  ＢＯＴＤＲ  Ｏｐｔｉｃａｌ  Ｆｉｂｅｒ  Ｓｔｒａｉｎ  Ｓｅｎｓｏｒ”，  Ｏｐｔｉｃａｌ  Ｆｉｂｅｒ  Ｔｅｃｈｎｏｌｏｇｙ  ７，  ４５−６４，（２００１）．
【非特許文献２】Ｋ．  Ｈｏｔａｔｅ  ａｎｄ  Ｔ．  Ｈａｓｅｇａｗａ，  ＩＥＩＣＥ  Ｔｒａｎｓ．  Ｅｌｅｃｔｒｏｎ．  Ｅ８３−Ｃ，  ４０５，  （２０００）．
【非特許文献３】Ｙ．Ｊ．  Ｒａｏ，  “Ｒｅｃｅｎｔ  ｐｒｏｇｒｅｓｓ  ｉｎ  ａｐｐｌｉｃａｔｉｏｎｓ  ｏｆ  ｉｎ−ｆｉｂｒｅ  Ｂｒａｇｇ  ｇｒａｔｉｎｇ  ｓｅｎｓｏｒｓ，”  Ｏｐｔｉｃｓ  ａｎｄ  Ｌａｓｅｒｓ  ｉｎ  Ｅｎｇｉｎｅｅｒｉｎｇ  ３１，  ２９７−３２４，  （１９９９）．
【非特許文献４】Ｈ．  Ｔａｋａｈａｓｈｉ，  Ｆ．  Ｉｔｏ，  Ｃ．  Ｋｉｔｏ，  ａｎｄ  Ｋ．  Ｔｏｇｅ，  “Ｉｎｄｉｖｉｄｕａｌ  ｌｏｓｓ  ｄｉｓｔｒｉｂｕｔｉｏｎ  ｍｅａｓｕｒｅｍｅｎｔ  ｉｎ  ３２−ｂｒａｎｃｈｅｄ  ＰＯＮ  ｕｓｉｎｇ  ｐｕｌｓｅｄ  ｐｕｍｐ−ｐｒｏｂｅ  Ｂｒｉｌｌｏｕｉｎ  Ａｎａｌｙｓｉｓ”，  Ｏｐｔｉｃｓ  Ｅｘｐｒｅｓｓ，  Ｖｏｌ．２１，  Ｎｏ．６，  ６７３９，  （２０１３）．
【非特許文献５】Ｈ．  Ｔａｋａｈａｓｈｉ，  Ｋ．  Ｔｏｇｅ，  Ｃ．  Ｋｉｔｏ，  ａｎｄ  Ｆ．  Ｉｔｏ，“Ｉｎｄｉｖｉｄｕａｌ  ＰＯＮ  Ｍｏｎｉｔｏｒｉｎｇ  Ｕｓｉｎｇ  ｍａｉｎｔｅｎａｎｃｅ  Ｂａｎｄ  Ｐｕｌｓｅｄ  Ｐｕｍｐ−Ｐｒｏｂｅ  Ｂｒｉｌｌｏｕｉｎ  Ａｎａｌｙｓｉｓ”，  ２０１３  １８ｔｈ  ＯＥＣＣ，  ＴｈＰ１−４，  （２０１３）．
【発明の概要】
【発明が解決しようとする課題】
【０００５】
非特許文献４および５に示される技術では、被測定光ファイバに一対のプローブ光パルスとポンプ光パルス（以下、パルス対と称する）が入射された後は、このパルス対が被測定光ファイバから排出されるのを待たねばならない。つまり第１波のパルス対を被測定光ファイバに入射後このパルス対が被測定光ファイバから消滅するまで、第２波以降のパルス対を被測定光ファイバに入射することができない。これは、プローブ光とポンプ光とが光ファイバ内で多重に衝突するとブリルアン利得ピークから正確な特性分布情報を取得できなくなるからで、よって測定に時間のかかることが課題である。
【０００６】
しかも非特許文献４および５に示される技術では、ＳＮ比（信号対雑音比）を改善するために数万回のオーダで計測を繰り返して平均値を求めるようにしているので、必要な時間はますます長くなる。まして被測定線路が長いケースや分解能向上のため測定点を多点化したケースにおいては結果を得られるまでの時間はさらに長くなるので、抜本的な対処を求められている。
【０００７】
そこで、本発明は、上記課題を解決するために、分岐光線路の個別の特性分布を高速で測定可能な分岐光線路特性解析方法および分岐光線路特性解析装置を提供することを目的とする。さらに、当該分岐光線路特性解析装置で特性分布を高速で測定可能な分岐光線路の構造を特定できる分岐光線路設計方法を提供することも目的とする。
【課題を解決するための手段】
【０００８】
本発明に係る分岐光線路特性解析方法及び分岐光線路特性解析装置は、光スプリッタの上流側から光周波数の等しいＸ個（Ｘは２以上の自然数）のプローブ光パルスを入射し、その後、プローブ光パルスと異なる光周波数の１つのポンプ光パルスを入射し、発生したブリルアン散乱光を受光して解析することとした。
【０００９】
具体的には、本発明に係る分岐光線路特性解析方法は、光スプリッタに接続される複数の分岐光線路に、前記分岐光線路でブリルアン（Ｂｒｉｌｌｏｕｉｎ）散乱を発生させるポンプ光パルスと、前記分岐光線路の前記光スプリッタと反対側の遠端で反射させて前記ポンプ光パルスと相互作用させるプローブ光パルスと、を伝搬させ、前記分岐光線路の特性を解析する分岐光線路特性解析方法であって、
前記プローブ光パルスと前記ポンプ光パルスとの光周波数差である設定ブリルアン周波数シフトｆＢを設定する光周波数差設定ステップと、
パルス間隔が等しく且つ光周波数の等しいＸ個（Ｘは２以上の自然数）の前記プローブ光パルスを発生し、Ｘ個の前記プローブ光パルスを前記光スプリッタを経由して前記分岐光線路に入射するプローブ光パルス入射ステップと、
前記プローブ光パルス入射ステップ後に、前記ポンプ光パルスを発生し、Ｘ個のプローブ光パルスが入射された後に前記ポンプ光パルスを前記光スプリッタを経由して前記分岐光線路に入射するポンプ光パルス入射ステップと、
前記ポンプ光パルス入射ステップの後、前記分岐光線路からの戻り光パルスを前記光スプリッタを介して受光し、電気信号に変換する戻り光受光ステップと、
前記戻り光受光ステップで変換した前記電気信号から設定ブリルアン周波数シフトｆＢにおける前記戻り光パルスのブリルアン利得スペクトル分布を取得する特性解析ステップと、
を行うことを特徴とする。
【００１０】
本分岐光線路特性解析方法は、ポンプ光パルス１つに対してＸ個のプローブ光パルスを入射するため、１回の作業で分岐光線路の長さ方向にＸ箇所の測定点におけるブリルアン散乱光を受光することができる。そして、本分岐光線路特性解析方法は、プローブ光パルス列とポンプ光パルスとの入射時間差を変化させることで測定点をずらし、この作業を繰り返すことで１回目の作業の測定点間のブリルアン散乱光を受光することができる。すなわち、本分岐光線路特性解析方法は、既存の手法と比較して、単位時間当たりＸ倍の速度で測定点の情報を得ることができ、測定時間をＸ分の１に短縮することが可能となる。
【００１１】
従って、本発明は分岐光線路の個別の特性分布を高速で測定可能な分岐光線路特性解析方法を提供することができる。
【００１２】
本発明に係る分岐光線路特性解析方法は、前記光周波数差設定ステップの前に、最長の前記分岐光線路と最短の前記分岐光線路との長さの差ΔＬｍａｘを検出する光路長差検出ステップを行うことを特徴とする。それぞれの分岐光線路の長さは異なる。そこで、本分岐光線路特性解析方法は、予め最長の分岐光線路の長さと最短の分岐光線路の長さとの差を取得しておく。
【００１３】
本発明に係る分岐光線路特性解析方法は、前記プローブ光パルス入射ステップは、前記分岐光線路それぞれからの前記戻り光パルスの重畳を防止するように前記プローブ光パルスのパルス間隔Ｔｓを設定することを特徴とする。それぞれの分岐光線路からのブリルアン散乱光が重畳すると分岐光線路の特性を測定できなくなるので、プローブ光パルスのパルス間隔でブリルアン散乱光（戻り光）が重畳しないように調整する。
【００１４】
なお、本発明に係る分岐光線路特性解析方法は、前記分岐光線路の線長差が等間隔ΔＬである場合、前記プローブ光パルス入射ステップは、τｐｒｏｂｅを前記プローブ光パルスのパルス幅、νを前記分岐光線路中の光速、Ｎを前記分岐光線路の分岐数とすれば、
Ｔｓ＝τｐｒｏｂｅ＋２ΔＬｍａｘ／ν
でパルス間隔Ｔｓを設定することを特徴とする。ただし、２ΔＬ＝τｐｒｏｂｅνを満たし、ΔＬｍａｘ＝ΔＬ（Ｎ−１）を満たす。
【００１５】
なお、本発明に係る分岐光線路特性解析方法は、所定の光周波数範囲内で前記設定ブリルアン周波数シフトｆＢを変化させ、前記光周波数差設定ステップから前記特性解析ステップを繰り返すことを特徴とする。ｆＢを変化させることで測定対象のブリルアン散乱光の光周波数を変えることができる。詳細には、光ファイバの特性（温度や曲げ等）は、固有のブリルアン散乱光の光周波数を持つ。このため、ｆＢを変化させることで、測定対象の光ファイバ特性を変えることができる。
【００１６】
本発明に係る分岐光線路特性解析方法を実現できる分岐光線路特性解析装置の構造は、次の通りである。
【００１７】
本発明に係る分岐光線路特性解析装置は、光スプリッタに接続される複数の分岐光線路に、前記分岐光線路でブリルアン（Ｂｒｉｌｌｏｕｉｎ）散乱を発生させるポンプ光パルスと、前記分岐光線路の前記光スプリッタと反対側の遠端で反射させて前記ポンプ光パルスと相互作用させるプローブ光パルスと、を伝搬させ、前記分岐光線路の特性を解析する分岐光線路特性解析装置であって、
パルス間隔が等しく且つ光周波数の等しいＸ個（Ｘは２以上の自然数）の前記プローブ光パルスを発生し、Ｘ個の前記プローブ光パルスを前記光スプリッタを経由して前記分岐光線路に入射するプローブ光パルス入射部と、
前記ポンプ光パルスを発生し、Ｘ個のプローブ光パルスが入射された後に前記ポンプ光パルスを前記光スプリッタを経由して前記分岐光線路に入射するポンプ光パルス入射部と、
前記分岐光線路からの戻り光パルスを前記光スプリッタを介して受光し、電気信号に変換する戻り光受光部と、
前記戻り光受光部が変換した前記電気信号から設定ブリルアン周波数シフトｆＢにおける前記戻り光パルスのブリルアン利得スペクトル分布を取得する特性解析部と、
前記プローブ光パルスと前記ポンプ光パルスとの光周波数差である前記設定ブリルアン周波数シフトｆＢを設定し、
前記プローブ光パルス入射部に、前記分岐光線路それぞれからの前記戻り光パルスの重畳を防止するように前記プローブ光パルスのパルス間隔Ｔｓを設定させ、
前記プローブ光パルス入射部に前記プローブ光パルスを前記分岐光線路に入射させた後に前記ポンプ光パルス入射部に前記ポンプ光パルスを前記分岐光線路に入射させる制御部と、
を備えることを特徴とする。
【００１８】
さらに、本発明に係る分岐光線路特性解析装置の前記制御部は、所定の光周波数範囲内で前記設定ブリルアン周波数シフトｆＢを変化させ、前記プローブ光パルス入射部と前記ポンプ光パルス入射部に前記プローブ光パルスと前記ポンプ光パルスを繰り返し前記分岐光線路に入射させることを特徴とする。
【００１９】
従って、本発明は分岐光線路の個別の特性分布を高速で測定可能な分岐光線路特性解析装置を提供することができる。
【００２０】
続いて、本発明に係る分岐光線路特性解析装置で特性分布を高速で測定可能な分岐光線路の構造を特定できる分岐光線路設計方法を説明する。
【００２１】
本発明に係る分岐光線路特性解析装置および分岐光線路特性解析方法の測定性能（所望の感度を達成するための測定時間、あるいは一定時間に達成する測定感度）は、被測定線路が有する分岐数によって変化する。分岐数が多い場合、分岐損失により分岐下部での試験光パワーが小さくなるため、信号対雑音比（ＳＮＲ）は劣化する。一方で、分岐数が少ない場合、分岐損失が小さいためポンプディプレッションを相対的に大きく生じやすく、ポンプディプレッションを抑えるために低い試験光パワーを使用することとなり、やはりＳＮＲが劣化する。すなわち、分岐光線路測定対象の総延長が与えられた際に、測定性能を最大化するための、分岐損失が大きすぎず、ポンプディプレッションが生じにくい最適な分岐数が必ず存在する。
【００２２】
具体的には、本発明に係る分岐光線路設計方法は、  請求項１から５のいずれかに記載の分岐光線路特性解析方法で特性を解析する分岐光線路の分岐数Ｎを、ブリルアン利得が最大となるように決定する分岐光線路設計方法である。
ただし、前記ブリルアン利得は数（Ｃ１）で定義される。
【数Ｃ１】
2017026512.tif 000003
ここで、ｇＢはブリルアン利得係数、Ａｅｆｆは光ファイバのコアの実効断面積、ｌｃは前記プローブ光パルスと前記ポンプ光パルスの相関長、Ｌは分岐下部心線の長さ、αは光ファイバの光損失係数、Ｐｒ（Ｌ）は前記分岐光線路の遠端での前記プローブ光パルスのパワー、Ｐｕ（Ｌ）は前記分岐光線路の遠端での前記ポンプ光パルスのパワーである。。
【００２３】
従って、本発明は、前記分岐光線路特性解析装置で特性分布を高速で測定可能な分岐光線路の構造を特定できる分岐光線路設計方法を提供することができる。
【発明の効果】
【００２４】
本発明は、分岐光線路の個別の特性分布を高速で測定可能な分岐光線路特性解析方法および分岐光線路特性解析装置を提供することができる。また、本発明は、前記分岐光線路特性解析装置で特性分布を高速で測定可能な分岐光線路の構造を特定できる分岐光線路設計方法を提供することもできる。
【図面の簡単な説明】
【００２５】
【図１】本発明に係る分岐光線路特性解析装置を説明するブロック図である。
【図２】本発明に係る分岐光線路特性解析方法を説明するフローチャートである。
【図３】本発明に係る分岐光線路特性解析方法の原理を説明する模式図である。
【図４】本発明に係る分岐光線路特性解析方法で測定した分岐光線路毎のブリルアン利得スペクトル分布を示す図である。
【図５】本発明に係る分岐光線路特性解析装置において、最適な入射プローブ光パワーの計算結果を示す図である。
【図６】本発明に係る分岐光線路設計方法において、信号対雑音比を最大化する分岐数の計算結果の例である。
【発明を実施するための形態】
【００２６】
添付の図面を参照して本発明の実施の形態を説明する。以下に示される実施の形態は一例であり、本発明は以下の実施の形態に制限されるものではない。なお、本明細書及び図面において符号が同じ構成要素は、相互に同一のものを示すものとする。
【００２７】
［分岐光線路特性解析装置］
図１は、本実施形態に係る分岐光線路特性解析装置の一例を示すブロック図である。
本分岐光線路特性解析装置は、光スプリッタに接続される複数の分岐光線路に、前記分岐光線路でブリルアン（Ｂｒｉｌｌｏｕｉｎ）散乱を発生させるポンプ光パルスと、前記分岐光線路の前記光スプリッタと反対側の遠端で反射させて前記ポンプ光パルスと相互作用させるプローブ光パルスと、を伝搬させ、前記分岐光線路の特性を解析する分岐光線路特性解析装置であって、
パルス間隔が等しく且つ光周波数の等しいＸ個（Ｘは２以上の自然数）の前記プローブ光パルスを発生し、Ｘ個の前記プローブ光パルスを前記光スプリッタを経由して前記分岐光線路に入射するプローブ光パルス入射部と、
前記ポンプ光パルスを発生し、Ｘ個のプローブ光パルスが入射された後に前記ポンプ光パルスを前記光スプリッタを経由して前記分岐光線路に入射するポンプ光パルス入射部と、
前記分岐光線路からの戻り光パルスを前記光スプリッタを介して受光し、電気信号に変換する戻り光受光部と、
前記戻り光受光部が変換した前記電気信号から設定ブリルアン周波数シフトｆＢにおける前記戻り光パルスのブリルアン利得スペクトル分布を取得する特性解析部と、
前記プローブ光パルスと前記ポンプ光パルスとの光周波数差である前記設定ブリルアン周波数シフトｆＢを設定し、
前記プローブ光パルス入射部に、前記分岐光線路それぞれからの前記戻り光パルスの重畳を防止するように前記プローブ光パルスのパルス間隔Ｔｓを設定させ、
前記プローブ光パルス入射部に前記プローブ光パルスを前記分岐光線路に入射させた後に前記ポンプ光パルス入射部に前記ポンプ光パルスを前記分岐光線路に入射させる制御部と、
を備えることを特徴とする。
【００２８】
図１において、プローブ光パルス入射部は、光源１０、分岐素子１１、光周波数変更手段１２、光パルス化手段１４、光増幅器１８、合波素子２０、及びサーキュレータ２１が相当する。ポンプ光パルス入射部は、光源１０、分岐素子１１、光パルス化手段１３、光増幅器１９、合波素子２０、及びサーキュレータ２１が相当する。光受光部は、サーキュレータ２１、光フィルタ２５、及び光受信器２６が相当する。特性解析部は、Ａ／Ｄ変換器２７及び演算処理装置２８が相当する。制御部は、正弦波発生器１７及び入射時間制御手段（１５、１６）が相当する。
【００２９】
本分岐光線路特性解析装置は、光スプリッタ２２を介して被測定光ファイバに試験光を入射し、その戻り光を解析して被測定光ファイバの光線路特性を解析することが可能である。光線路特性は、例えば、距離に対する光減衰量、曲げ障害の位置、曲げの程度、断線障害の位置、距離に対する温度変化量などである。試験光には２種類あり、一方はプローブ光と称される。他方はポンプ光と称され、分岐光線路２３の遠端に設置された反射フィルタ２４で反射したプローブ光と対向伝搬してプローブ光にブリルアン利得を与える。
【００３０】
プローブ光をパルス化してプローブ光パルスが得られる。ポンプ光をパルス化してポンプ光パルスが得られる。分岐光線路の遠端で反射され、ポンプ光パルスによるブリルアン増幅を受けて戻るプローブ光パルスを、戻り光パルスと称する。
【００３１】
図１において、レーザ光源１０から出力される光周波数ｆ０の連続光は分岐素子１１で２分岐され、一方は光周波数変更手段（光周波数変更手段）１２および光パルス化手段１４を経てプローブ光パルス（第１試験光）となる。他方は光パルス化手段１３を経てポンプ光パルス（第２試験光）となる。
【００３２】
光周波数変更手段１２は、ＬｉＮｂＯ３を用いた位相変調器、あるいはＳＳＢ（Ｓｉｎｇｌｅ  Ｓｉｄｅ  Ｂａｎｄ）変調器である。光周波数変更手段１２にオシレータ（正弦波発生器）１７から正弦波を与え、正弦波の周波数に応じて変調側波帯の周波数を変化させてプローブ光を外部変調し、周波数シフトを与えることができる。光周波数変更手段１２は、プローブ光の光周波数を測定しようとするブリルアン散乱光の光周波数へ変化させる。光周波数変更手段１２が変化させる当該光周波数（光周波数差Δｆ）を設定ブリルアン周波数シフトｆＢと称する。
【００３３】
より詳細に説明すれば、分岐光線路内で発生するブリルアン散乱光は、曲げや温度変化等でその光周波数が異なる。そしてプローブ光の光周波数に一致する光周波数のブリルアン散乱光が増幅され、測定に利用される。このため、光周波数変更手段１２に、曲げに対応する設定ブリルアン周波数シフトｆＢを設定してプローブ光の光周波数をｆ０から変化させｆ０−ｆＢとすれば曲げの発生を捉えることができ、温度変化に対応する設定ブリルアン周波数シフトｆＢを設定してプローブ光の光周波数をｆ０から変化させｆ０−ｆＢとすれば温度変化を捉えることができる。
【００３４】
そして、所定の光周波数範囲内で前記設定ブリルアン周波数シフトｆＢを変化させ、前記プローブ光パルスと前記ポンプ光パルスを繰り返し前記分岐光線路に入射させることで、分岐光線路の全体的な特性（曲げや温度変化等）を解析することができる。
【００３５】
周波数シフトを受けたプローブ光は、入射時間制御手段１６の制御に基づいて光パルス化手段１４でパルス化される。実施形態では特に、光周波数の等しい複数のプローブ光パルスを発生させる。一方、ポンプ光は、入射時間制御手段１５の制御に基づいて光パルス化手段１３でパルス化される。
【００３６】
光パルス化手段（１３、１４）は、例えば音響光学素子をパルス駆動する音響光学スイッチ、あるいはＬｉＮｂＯ３を用いた電気光学素子をパルス駆動する導波路スイッチなどを用いることができる。つまり光デバイス（音響光学変調器やＬｉＮｂＯ３変調器）を電気パルスで変調すれば光パルスを得られる。電気パルスで駆動された時間に応じて連続光がパルス化される。
【００３７】
入射時間制御手段１６は、複数の分岐光線路からの戻り光パルス同士が干渉することを防止するために、複数のプローブ光パルスのパルス間隔ＴＳを制御する。また、入射時間制御手段１５は、プローブ光パルスが被測定光ファイバに入射される時間（タイミング）とポンプ光パルスが被測定光ファイバに入射される時間（タイミング）とに、時間差を与える。
【００３８】
入射時間制御手段（１５、１６）は、例えば光パルス化手段を駆動する電気パルスの変調時間を変化させることで実施可能である。上記光デバイスを変調するタイミングを変化させることで光パルスの生成されるタイミングを制御することができる。光パルス化手段（１３、１４）の各出力光を合波素子２０で合波すれば、入射時間の異なるプローブ光パルス列（複数のプローブ光パルス）およびポンプ光パルスを生成することが可能となる。なお、光パルス化手段（１３、１４）に与えられる電気パルスのうち一方の電圧をゼロにすれば、プローブ光パルスまたはポンプ光パルスの一方だけを入射することも可能となる。
【００３９】
光パルス化手段１４からのプローブ光パルスは光増幅器１８で測定に必要なレベルにまで増幅され、合波素子２０に入射される。また、光パルス化手段１３からのポンプ光パルスは光増幅器１９で測定に必要なレベルにまで増幅され、合波素子２０に入射される。プローブ光およびポンプ光は合波素子２０から光サーキュレータ２１を介して基幹光ファイバに入射され、さらに被測定光ファイバに入射される。
【００４０】
被測定光ファイバは、光スプリッタ２２と、光スプリッタ２２の分岐下流側の分岐光線路２３とを備える。各分岐光線路２３はその終端に光反射フィルタ２４を有する。なお終端面が十分な反射率を有していれば光反射フィルタを省略してもよい。
【００４１】
光反射フィルタ２４で反射された試験光（プローブ光、ポンプ光）、あるいはブリルアン後方散乱光などの戻り光は光サーキュレータ２１経由で光フィルタ２５に入射される。光フィルタ２５は解析に必要なプローブ光だけを通過させる特性を持ち、プローブ光だけが光受信器２６に達する。光受信器２６は受光したプローブ光を光／電変換し、プローブに基づく電気信号（出力電流）をＡ／Ｄ変換器２７に入力する。Ａ／Ｄ変換器２７は電気信号をＡ／Ｄ（アナログ／ディジタル）変換し、ディジタルデータを生成する。このディジタルデータはパーソナルコンピュータ（ＰＣ）などの演算処理装置２８に入力される。
【００４２】
続いて、本分岐光線路特性解析装置の作用を詳しく説明する。
本分岐光線路特性解析装置では、以下に示される（条件１）〜（条件４）が満たされることが必要である。
【００４３】
（条件１）
光周波数変更手段１２による周波数シフトのレンジは、被測定光ファイバで生じるブリルアン周波数シフトの距離方向の周波数変化量よりも広いこと（光源１０の光周波数から測定しようとする特性に対応するブリルアン散乱光の光周波数へシフト可能であること）。条件１は、被測定光ファイバ中の全ての距離（位置）でプローブ光パルスとポンプ光パルスとがインタラクション（相互作用）して誘導ブリルアン散乱を生じるために必要な条件である。
【００４４】
（条件２）
プローブ光パルスのパルス幅τｐｒｏｂｅが、分岐光線路遠端（終端）からの戻り光の時間差２ｎΔＬ／ｃの最小値より狭いこと。なおΔＬは分岐光線路の長さの差の最小値である。また真空中の光速をｃとし、ｎは光ファイバの屈折率である。すなわちｃ／ｎは光ファイバ中の光速νとなる。
【００４５】
条件２は、分岐光線路毎の戻り光（誘導ブリルアン散乱光）が光受信器おいて重なること（干渉）を防ぐために必要な条件である。条件２が満たされなければ分岐光線路毎の戻り光（誘導ブリルアン散乱光）が光受信器において干渉し、各分岐光線路を時間的に切り分けることができなくなる。つまり各分岐光線路を区別できなくなる。
【００４６】
（条件３）
光受信器およびＡ／Ｄ変換器の帯域が、パルス幅τｐｒｏｂｅを受光可能な程度に広いこと。一般に、パルス幅τの光パルスを精確に測定するためには、光受信器およびＡ／Ｄ変換器の帯域が１／τより広い必要がある。
【００４７】
（条件４）
最長の分岐光線路と最短の分岐光線路との長さの差をΔＬｍａｘとし、光ファイバ中の光速をνとし、隣り合うプローブ光パルスの入射時間差をＴｓとしたとき、
Ｔｓ≧τｐｒｏｂｅ＋２ΔＬｍａｘ／ν
が満たされること。条件４も条件２と同様に、分岐光線路毎の戻り光パルスが光受信器おいて重なることを防ぐために必要な条件である。
【００４８】
なお、分岐光線路の長さの差のばらつきによっては、条件４の数式が測定速度に対して最適とならない場合もある。光受信器への到達時に戻り光パルス同士が干渉しない範囲で、Ｔｓを自由に設定して良い。ただし、分岐線路の長さの差が一様（線長差が等間隔ΔＬ）である場合、すなわちΔＬｍａｘ＝（Ｎ−１）ΔＬとなる場合に、プローブ光送出数を最大化でき、測定時間を短縮できる。この時、プローブ光送出数は、分岐下部心線長をＬとすると、Ｌ／（ＮΔＬ）となる。
【００４９】
［分岐光線路特性解析方法］
図２は、図１の分岐光線路特性解析装置が行う分岐光線路特性解析方法の一例を示すフローチャートである。本分岐光線路特性解析方法は、光スプリッタに接続される複数の分岐光線路に、前記分岐光線路でブリルアン（Ｂｒｉｌｌｏｕｉｎ）散乱を発生させるポンプ光パルスと、前記分岐光線路の前記光スプリッタと反対側の遠端で反射させて前記ポンプ光パルスと相互作用させるプローブ光パルスと、を伝搬させ、前記分岐光線路の特性を解析する分岐光線路特性解析方法であって、
前記プローブ光パルスと前記ポンプ光パルスとの光周波数差である設定ブリルアン周波数シフトｆＢを設定する光周波数差設定ステップと、
パルス間隔が等しく且つ光周波数の等しいＸ個（Ｘは２以上の自然数）の前記プローブ光パルスを発生し、Ｘ個の前記プローブ光パルスを前記光スプリッタを経由して前記分岐光線路に入射するプローブ光パルス入射ステップと、
前記プローブ光パルス入射ステップ後に、前記ポンプ光パルスを発生し、Ｘ個のプローブ光パルスが入射された後に前記ポンプ光パルスを前記光スプリッタを経由して前記分岐光線路に入射するポンプ光パルス入射ステップと、
前記ポンプ光パルス入射ステップの後、前記分岐光線路からの戻り光パルスを前記光スプリッタを介して受光し、電気信号に変換する戻り光受光ステップと、
前記戻り光受光ステップで変換した前記電気信号から設定ブリルアン周波数シフトｆＢにおける前記戻り光パルスのブリルアン利得スペクトル分布を取得する特性解析ステップと、
を行うことを特徴とする。
なお、図２において、ステップＳ１とＳ９が光周波数差設定ステップ、ステップＳ２とＳ３がプローブ光パルス入射ステップ、ステップＳ４がポンプ光パルス入射ステップ、ステップＳ５からＳ８が戻り光受光ステップと特性解析ステップに相当する。
【００５０】
本実施形態では、分岐光線路の特性の解析に先立ちΔＬｍａｘを測定する。つまり、本分岐光線路特性解析方法は、前記光周波数差設定ステップの前に、最長の前記分岐光線路と最短の前記分岐光線路との長さの差ΔＬｍａｘを検出する光路長差検出ステップを行うことを特徴とする。ΔＬｍａｘは例えば、プローブ光パルスを１波だけ入射して各分岐光線路からの戻り時間を測定し、最短の戻り時間と最長の戻り時間（いずれも往復時間）との差にνを乗算して２で割ることで求めることができる。
【００５１】
演算処理装置２８は、先ず、設定ブリルアン周波数シフトｆＢを設定する（ステップＳ１）。ｆＢは、プローブ光パルスとポンプ光パルスの光周波数差Δｆである。つまりプローブ光パルスとポンプ光パルスとの光周波数差Δｆ＝ｆＢであり、各プローブ光パルスの光周波数は光周波数ｆ０−ｆＢとなる。なお、ｆ０はポンプ光パルスの光周波数、ｆＢはブリルアン後方散乱光のうち測定対象とする光ファイバの特性に対応する光周波数である。ｆ０からｆＢだけ変化させるので光周波数シフト量（Ｂｒｉｌｌｏｕｉｎ  Ｆｒｅｑｕｅｎｃｙ  Ｓｈｉｆｔ：ＢＦＳ）と称する。
【００５２】
次に演算処理装置２８は、プローブ光パルス列に含まれる各プローブ光パルスとポンプ光パルスとの入射時間差ｔ１，ｔ２，ｔ３，・・・，ｔＸを設定する（ステップＳ２）。ステップＳ２において、前記分岐光線路それぞれからの前記戻り光パルスの重畳を防止するように前記プローブ光パルスのパルス間隔Ｔｓを設定することが必要である。
【００５３】
ここで、前記分岐光線路の線長差が等間隔ΔＬである場合、ステップＳ２では、τｐｒｏｂｅを前記プローブ光パルスのパルス幅、νを前記分岐光線路中の光速、Ｎを前記分岐光線路の分岐数とすれば、
Ｔｓ＝τｐｒｏｂｅ＋２ΔＬｍａｘ／ν
でパルス間隔Ｔｓを設定することを特徴とする。ただし、２ΔＬ＝τｐｒｏｂｅνを満たし、ΔＬｍａｘ＝ΔＬ（Ｎ−１）を満たす。
【００５４】
図３は、本実施形態における解析の手順を示す模式図である。図３において、白抜き矢印Ａが進行するプローブ光パルス列を示し、黒塗り矢印Ｂは進行するポンプ光パルスを示す。図３の（Ｉ）に示されるように、パルス幅τｐｒｏｂｅで入射されるプローブ光パルスは、互いに条件４に基づく時間差Ｔｓで入射される。すなわちプローブ光パルス列は、パルス間隔Ｔｓで並ぶＸ個のプローブ光パルスを含む。以下の説明では、Ｘ個（Ｘ波）を含むプローブ光パルス列のｉ番目のプローブ光パルスとポンプ光パルスの入射時間差をｔｉとする。この入射時間差ｔｉによりプローブ光パルスとポンプ光パルスとの衝突位置が変化し、各分岐光線路における長手方向の測定位置が決定される。つまり、ｔｉを変化させることで分岐光線路の距離方向全体を測定することができる。
【００５５】
図２に戻り、演算処理装置２８は、制御部（正弦波発生器１７、入射時間制御手段（１５、１６））に指示し、プローブ光パルス列を被測定光ファイバに入射させ、Ｘ個のプローブ光パルスの入射が完了した後にポンプ光パルスを被測定光ファイバに入射させる（ステップＳ３，Ｓ４）。上記したようにプローブ光パルス列に含まれるプローブ光パルスとポンプ光パルスはそれぞれｔ１，ｔ２，ｔ３，・・・，ｔＸの入射時間差を持つ。図３にはＸ＝３のケースが示される。
【００５６】
次に、演算処理装置２８は戻り光パルスの光受信器２６への到達時間（戻り時間）に基づいて、各戻り光パルスを反射した分岐光線路を個別に識別する（ステップＳ５）。つまりＡ／Ｄ変換器２７から出力されるディジタルデータがどの分岐光線路の情報を反映するのかが判断される。そして、得られたディジタルデータに基づいて、演算処理装置２８は誘導ブリルアン散乱光を解析する（ステップＳ６）。
【００５７】
図３の（ＩＩ）に示されるように、光スプリッタ２２でＮ分岐されたプローブ光パルス列とポンプ光パルスは分岐光線路中でインタラクションし、誘導ブリルアン散乱による後方散乱光が発生する。これによりプローブ光パルス列はポンプ光パルスからブリルアン利得を得る。プローブ光パルス列はＸ個のプローブ光パルスを含むので、ポンプ光パルスは異なる位置でＸ回のインタラクションを生じることになる。図３においては分岐光線路毎に３か所の衝突位置がある。
【００５８】
ブリルアン増幅を受けたＸ×Ｎ個のプローブ光パルスは、ポンプ光パルスを伴って光スプリッタ２２で合波されたのち光サーキュレータ２１に戻り、光フィルタ２５に入射される。光フィルタ２５はポンプ光パルスを除去し、これによりブリルアン増幅を受けたＸ×Ｎ個のプローブ光パルスだけが光受信器２６で受光される。（条件４）が満たされているので、図３の（ＩＩＩ）に示されるように、光受信器２６において戻り光パルスが重なることはない。
【００５９】
光受信器２６からの電気信号はＡ／Ｄ変換器２７でディジタルデータに変換され、演算処理装置２８に入力される。演算処理装置２８は図２のステップＳ６の解析を行い、解析結果を出力する（ステップＳ７）。
【００６０】
最長の分岐光線路からポンプ光パルスが受光されると（ステップＳ８でＹｅｓ）１つの小シーケンスが終了する（ステップＳ２〜ステップＳ７）。この小シーケンスは入射時間差ｔ１，ｔ２，ｔ３，・・・，ｔＸをδｔだけ増加させつつ（ステップＳ１０）、ｔｉがｔｉ＋１の初期値と等しくなるまで繰り返される。
【００６１】
小シーケンスが完了すると、所定の光周波数範囲内で前記設定ブリルアン周波数シフトｆＢを変化させ、前記光周波数差設定ステップから前記特性解析ステップを繰り返す。
具体的には、演算処理装置はプローブ光パルスとポンプ光パルスとの光周波数差Δｆ（＝ｆＢ）をδｆだけ変化させつつ（ステップＳ１２）、ステップＳ１〜ステップＳ９の大シーケンスを繰り返す。この大シーケンスはｆＢ＝ＦＢになるまで繰り返される（ステップＳ１１でＹｅｓ）。ここでＦＢは設定ブリルアン周波数の全帯域である。ステップＳ１１までの手順が完了すると、既に得られたブリルアン利得スペクトル分布に基づいて距離ごとのブリルアン利得スペクトラムのピーク分布を解析することができる。
【００６２】
続いて、数式を用いて上記作用を詳しく説明する。
【００６３】
（ｉ）プローブ光パルス列とポンプ光パルスによる誘導ブリルアン散乱の測定について
プローブ光パルス列に含まれるプローブ光パルスとポンプ光パルスとがインタラクションすると誘導ブリルアン散乱が発生する。プローブ光パルス列とポンプ光パルスとの光周波数差ΔｆがｆＢのとき、プローブ光パルスは誘導ブリルアン散乱により式（１）に示される増幅を受ける。
【数１】
2017026512.tif 000004
【００６４】
式（１）のαＢ（ｚ，ｆ）は入射端からｚの位置でインタラクションし、ブリルアン周波数差がｆのときの誘導ブリルアンによる利得を示す。Ｗｐｕｍｐはポンプ光のパルス幅を距離の次元で表したものである。ｇＢ（ｆ）はプローブ光パルスとポンプ光パルスの周波数差がｆであるときの誘導ブリルアン散乱係数を示す。ｚは分岐光線路入射端からプローブ光パルスとポンプ光パルスがインタラクションした位置までの距離を示す。Ｉｐｕｍｐ（ｚ）は分岐光線路入射端から距離ｚだけ離れた位置におけるポンプ光パルスの強度を示す。
【００６５】
分岐光線路＃ｉ（図３ではｉは１〜４の自然数）の損失係数をαｉ、分岐光線路＃ｉを往復する光の全損失をＬｉとする。そうすると分岐光線路に入射されたのち終端で反射され、入射端から距離ｚの位置でポンプ光パルスと衝突したプローブ光パルスの当該入射端における強度Ｉｐｒｏｂｅ（２Ｌｉ，ｚ）は、式（２）で表される。
【数２】
2017026512.tif 000005
【００６６】
式（２）に示されるように、分岐光線路入射端でのプローブ光パルスの強度Ｉｐｒｏｂｅ（２Ｌｉ，ｚ，ｆ）は、ｇＢ（ｆ）とＩｐｕｍｐ（ｚ）の関数になる。Ｉｐｕｍｐ（ｚ）は式（３）で表される。
【数３】
2017026512.tif 000006
また、プローブ光パルスだけを入射した場合の入射端に戻る反射プローブ光パルス強度Ｉｒｅｆ（２Ｌｉ）は、式（４）に示される。
【数４】
2017026512.tif 000007
従って式（２）を式（３）および式（４）を用いて変換すると式（５）が得られる。
【数５】
2017026512.tif 000008
【００６７】
ブリルアン散乱光の利得を解析することにより、式（５）に基づいてインタラクションした場所までの損失をポンプ光パルス幅で積分した値とブリルアン散乱係数との積を算出できる。従って、被測定光ファイバにおける損失を加味したブリルアン利得強度を得ることができる。
【００６８】
なお本実施形態では、ポンプ光パルスが複数のプローブ光パルスと衝突するので複数回のインタラクションによるポンプディプレッションを生じ、ブリルアン利得強度が式（５）よりも抑圧されることが考えられる。しかし実施形態の手法によれば、後述の通り、ブリルアン利得強度の絶対値情報は必要無いので、ポンプディプレッションへの特段の配慮は必要ない。必要な情報は、ブリルアン利得強度がピーク値を取るブリルアン周波数シフト量であり、ブリルアン利得強度のピーク値をモニタ（ピークサーチ）できれば良いからである。測定距離ごとのブリルアン利得スペクトルは、プローブ光パルスとポンプ光パルスとの周波数差ｆを変化させて測定を繰り返すことで取得される。
【００６９】
（ｉｉ）分岐光線路の距離に対するブリルアン散乱光分布の測定について
簡単のため、プローブ光パルス列のうち先頭のプローブ光パルスが、ポンプ光パルスに対して時間ｔ１だけ先に入射する場合だけを説明する。プローブ光パルス列に含まれる他のプローブ光パルスについても、時間ｔ１をｔ２，ｔ３，・・・，ｔＸとすれば同様の説明が成り立つ。
【００７０】
被測定光ファイバの入射端から分岐光線路＃ａ（ａは１≦ａ≦Ｎの自然数）の終端までの長さをＬａとする。プローブ光パルスは光反射フィルタにより反射される。分岐光線路終端からの距離をｌとし、光ファイバの屈折率をｎとし、真空中の光速をｃとすると、反射されたプローブ光パルスはｔ１／２秒後にｌ＝ｃ／ｎ×ｔ１／２だけ進むので、入射端からの距離をｌｘ１とすると式（６）が得られる。
【数６】
2017026512.tif 000009
【００７１】
プローブ光パルスが入射されてから分岐光線路に進入し、光反射フィルタで反射してｌｘ１に達するまでの時間ｔは、式（７）で表される。
【数７】
2017026512.tif 000010
【００７２】
プローブ光パルスが入射されてからｔ１秒後にポンプ光パルスを入射するとする。ポンプ光パルスがｔ秒後に到達する入射端からの距離ｌｘ２は、式（８）で表される。
【数８】
2017026512.tif 000011
【００７３】
ｌｘ１＝ｌｘ２の位置、つまり式（６）＝式（８）となる位置でポンプ光パルスはプローブ光パルスとインタラクションする。インタラクションのタイミングは光反射フィルタでプローブ光パルスが反射されてからｔ１／２秒後である。つまり、プローブ光パルスとポンプ光パルスとを被測定光ファイバに入射する時間差を変化させることで、プローブ光パルスとポンプ光パルスとがインタラクションする位置を制御できる。これにより、距離に対する誘導ブリルアン散乱の特性分布を求めることができる。
【００７４】
本実施形態において、プローブ光パルス列とポンプ光パルスが被測定光ファイバを一往復する間にインタラクションする位置は、プローブ光パルス列に含まれるプローブ光パルスの数Ｘだけ存在する。つまり多点（Ｘ点）の測定点のブリルアン利得情報を一括して取得することができるので、プローブ光パルスとポンプ光パルスの入射時間差を変化させて測定する回数をＸ分の１に削減でき、測定時間をＸ分の１に削減することができる。
【００７５】
（ｉｉｉ）分岐光線路＃ａの誘導ブリルアン散乱光が光受信器に到達する時間の測定について
プローブ光パルスが光受信器に到達する時間をｔｄａとする。プローブ光パルスは、分岐光線路終端の光フィルタで反射されて光受信器に入射される。プローブ光パルスの光受信器への到達時間は式（９）で表される。
【数９】
2017026512.tif 000012
【００７６】
他の分岐光線路＃ｂ（１≦ｂ≦Ｎの整数）から戻るプローブ光パルスが光受信器に到達する時間ｔｄｂは、式（１０）で表される。
【数１０】
2017026512.tif 000013
よって、光受信器に戻る時間差は、式（１１）で表される。
【数１１】
2017026512.tif 000014
【００７７】
Ｌａ≠Ｌｂであれば、光受信器に到達する時間が異なる。
（条件２）は、式（１２）で表される。
【数１２】
2017026512.tif 000015
【００７８】
式（１２）および（条件４）が満たされていれば、各分岐光線路から戻ったプローブ光パルスは光スプリッタで合波された後も重ならない。よって光受信器への到達時間を測定することで、どの分岐光線路から戻ってきたプローブ光パルスであるかを特定できる。つまり各分岐光線路を時間的に切り分けることができる。
【００７９】
上記（ｉ）〜（ｉｉｉ）により、分岐光線路個別のブリルアン利得スペクトル分布を、既存の技術（特に非特許文献４）のＸ倍の速度で高速に測定することが可能になる。測定されたブリルアン利得ピークの分布から温度、歪み、ファイバ構造パラメータの違いの特性分布を算出することももちろん可能である。
【００８０】
以上述べたように実施形態では、演算処理装置は、或る設定ブリルアン周波数シフトｆＢでの測定が完了した後、設定ブリルアン周波数シフトを設定ブリルアン周波数シフト間隔変化させて測定を繰り返し、設定ブリルアン周波数全帯域幅ＦＢになると測定を終了する。
【００８１】
得られた距離に対するブリルアン利得分布を周波数毎に並べることで、距離に対するブリルアン利得スペクトル分布を求める。このブリルアン利得スペクトル分布において、演算処理により、距離ごとにブリルアン利得スペクトルの強度ピークをとる周波数がブリルアン周波数シフト量であり、被測定光ファイバ長手方向のブリルアン周波数シフト量の測定結果を出力することで、距離に対する特性分布を求める。つまり温度や歪、ファイバ構造パラメータの違いなどの、ブリルアン周波数シフト量の変化をもたらす物理量の分岐光線路毎の距離方向の特性分布を算出することができる。
【００８２】
以上の測定において、実施形態によれば、既存の手法に比べて単位時間当たりＸ倍の速度で測定点の情報を得ることができる。従って測定時間をＸ分の１にすることが可能になる。
【００８３】
図４は、分岐光線路毎のブリルアン利得スペクトル分布を示す図である。図４は＃１〜＃８の８つの分岐光線路を測定した結果である。各グラフは直行する３つの軸に距離、光周波数シフト量（Ｂｒｉｌｌｏｕｉｎ  Ｆｒｅｑｕｅｎｃｙ  Ｓｈｉｆｔ：ＢＦＳ）および損失をそれぞれプロットしたグラフである。なお分岐光線路＃１，＃７，＃８に、基幹光ファイバとファイバ構造パラメータの異なる光ファイバを用いた。
【００８４】
グラフ上の段差から一見して分かるのは、９．５ｋｍの距離に光スプリッタがあり、９．５ｋｍ以遠が分岐光線路であることである。また各グラフによれば、分岐光線路＃１、＃７、＃８のブリルアン周波数シフト量が基幹光ファイバと異なることが分かる。ブリルアン周波数シフト量は歪みや温度変化に対して線形な応答変化を示すので、実施形態によれば各分岐光線路を識別可能となり、かつ、歪みや温度変化を分岐光線路毎に測定することができる。
【００８５】
以上のように実施形態では、光スプリッタで複数に分岐された各分岐光線路の長さの差ΔＬを利用し、波長の異なる二種のパルス試験光を用意する。予め設定された入射時間差を持つＸ個の試験光パルス列（プローブ光パルス列）を先に入射する。この入射時間差は分岐光線路の長さの最大値に基づいて最適化されるので、光受信器でプローブ光パルスどうしが重なることはない。全てのプローブ光パルスの入射が完了したのち、ポンプ光パルスが入射される。
【００８６】
プローブ光パルス列は分岐光線路遠端、あるいは遠端に設置される反射フィルタで反射される。これによりポンプ光パルスが分岐光線路内の各位置（Ｘ箇所）で対向衝突し、誘導ブリルアン後方散乱光が生じる。この誘導ブリルアン後方散乱光を光受信器で受信し、光受信器の出力電流を２ｎΔＬ／ｃ（ｃは光速）よりも高い時間分解能で解析することで、分岐光線路のどの心線から戻った誘導ブリルアン後方散乱光であるかを特定することが可能となる。
【００８７】
また、プローブ光パルスとポンプ光パルスとの光周波数差Δｆ（＝ｆＢ）を変化させながら測定を繰り返すことで分岐光線路毎のブリルアン利得スペクトル分布を取得する。そして、取得されたブリルアン利得スペクトル分布において、距離ごとにブリルアン利得スペクトルがピーク値を取るブリルアン周波数シフト量を取得（ピークサーチ）する。ブリルアン周波数シフト量は温度変化や歪み変化に対して線形な変化応答を示すので、ブリルアン周波数シフト量に寄与する物理量変化を分岐光線路毎に測定することができる。
【００８８】
要するに、或る設定ブリルアン周波数シフトｆＢでの距離方向のブリルアン利得スペクトルを得るのに、既存の技術ではＸ回（分岐光線路長とサンプリング距離間隔に依存）の測定を要していたのに対し、実施形態によれば１回の測定で足りる。すなわち、予めプローブ光パルスを１波だけ入射して（Ｎ分岐の被測定光ファイバからの戻り時間から）分岐光線路の最長と最短の差ΔＬｍａｘを測定したうえで、周波数の等しいＸ個のプローブ光パルス（周波数ｆ０−ｆＢ）をτ＋２ΔＬｍａｘ／ν以上の時間差をもって入射するプローブ光パルス列と、プローブ光パルス列とは周波数の異なる１個のポンプ光パルス（周波数ｆ０）を更に時間差をもって入射することで、プローブ光パルス列とポンプ光パルスとの対向伝搬によりブリルアン利得スペクトルを求める。
【００８９】
そして、或る設定ブリルアン周波数シフト（ｆＢ）での測定が完了した後、プローブ光パルス列とポンプ光パルスの周波数差を設定ブリルアンシフト間隔（δｆ）ずつ変化させて設定ブリルアン周波数全帯域（ＦＢ）まで測定を繰り返し、距離方向に対するブリルアン利得分布を得ることで、光線路の特性分布を求めることができる。
【００９０】
以上から実施形態によれば、分岐光線路の個別の特性を測定することができる。しかも分岐光線路毎にＸ個のデータを一気に取得できるので、試験光入射待機時間を大幅に削減でき、従って測定時間を大幅に短縮することができる。
【００９１】
これらのことから本実施形態によれば、分岐光線路の個別の特性分布を高速で測定することの可能な分岐光線路の特性解析装置および分岐光線路の特性解析方法を提供することが可能となる。
【００９２】
さらに、本発明に関わる分岐線路の特性解析装置および特性解析方法の測定時間を最も短くするためには、測定対象の分岐数Ｎを全センシング長（≒ＮＬ）に対して最適化する必要がある。また、入射するポンプディプレッションを抑え、なおかつＳＮＲが最大化できる入射光パワーを選定して使用する必要がある。以下、ＳＮＲを最大化する入射光パワーの選定方法と最適分岐数の設計について述べる。なお、本発明装置及び方法においてＳＮＲが最も悪くなる測定地点は、常に分岐下部心線遠端であるため、分岐下部心線遠端でのＳＮＲが最大となる入射光パワーに焦点を当てる。
【００９３】
光ファイバ中で生じるプローブ光とポンプ光のブリルアン利得は式（１３），（１４）の結合方程式にて表現できる。
【数１３】
2017026512.tif 000016
【数１４】
2017026512.tif 000017
【００９４】
式（１３），（１４）において、Ｐｒはポンプ光と衝突時のプローブ光パワー、Ｐｕはプローブ光と衝突時のポンプ光パワー、ｇＢはブリルアン利得係数、Ａｅｆｆはコアの実効断面積、αは光ファイバの光損失係数である。式（１３），（１４）の結合方程式を数値解析的に解くことで、プローブ光およびポンプ光のパワーの距離に対する変化、および遠端でのプローブ光およびポンプ光パワーを明らかにすることができる。詳細を以下に述べる。
【００９５】
まず、プローブ光パワーの距離に対する変化について、分岐数Ｎに応じた分岐損失と片端測定であることを考慮し、ｚ地点でポンプ光と衝突するプローブパワーは、
【数１５】
2017026512.tif 000018
Ｐｒｉは分岐点直前におけるプローブ光パワー、αは光ファイバ損失係数、Ｌは分岐下部心線長である。
【００９６】
式（１５）より分岐下部心線遠端Ｌ地点でのプローブ光パワーは、
【数１６】
2017026512.tif 000019
となる。
【００９７】
微小区間ｄｚ＝ｚｎ−ｚｎ−１におけるポンプ光減衰量は、分岐点直前におけるポンプパワーＰｕｉを初期値とし、ポンプ光とプローブ光の相関長をｌｃとして、下記漸化式で表せる。
【数１７】
2017026512.tif 000020
【００９８】
遠端Ｌ地点で生じるブリルアン利得受光パワーは、
【数１８】
2017026512.tif 000021
【００９９】
式（１８）が最大となるＰｒｉおよびＰｕｉを数値解析で求めることで、ＳＮＲを最大化することができる。ただし、一般に実験系の制約から各々の入射パワー上限が存在する。入射ポンプ光パワーは不要な非線形光学効果（誘導ブリルアン散乱など）が十分抑えられる範囲で大きいほどＳＮＲは向上する。入射ポンプ光パワーはポンプディプレッションを十分に抑えられる範囲で大きいほどＳＮＲは向上する。例として、１５  ｄＢｍを上限とする入射プローブ光のパワーの最適値の計算結果の例（分岐数Ｎ＝８）を図５に示す。１５ｄＢｍの入射ポンプ光パワーでは誘導ブリルアン散乱は十分に小さい。
【０１００】
図５より、ポンプディプレッションを抑えつつ、最大のブリルアン利得が生じるような入射プローブ光パワーの最適値が存在することがわかる。最適な入射プローブ光パワーは被測定線路の分岐数、長さ、分岐心線長差、線路損失に依存して変化するため、入射プローブ光パワーの最適値を計算する前に予め、被測定線路の分岐数、長さ、分岐心線長さ、線路損失を把握しておく必要がある。被測定線路の分岐数、長さ、分岐心線長差、線路損失は、例えば、単一パルス光を線路に入射することにより、その戻りのパルス数、戻り時間、戻り時間差、受光パワーからそれぞれ容易に把握することができる。
【０１０１】
また、式（１８）について、分岐数をパラメータとして数値解析することで、ＳＮＲを最大化する分岐数もわかる。最適分岐数の計算結果の例を図６に示す。図は、分岐下部心線遠端で生じたブリルアン利得の受光パワーについて、分岐数が１の場合に対する各分岐構成時のＳＮＲのＦＯＭ（Ｆｉｇｕｒｅ  ｏｆ  ｍｅｒｉｔ）を示している。センシングすべき総センシング長が与えられた場合に、分岐数の変化に対して、ＦＯＭは極値を持つことがわかる。すなわち、極値をとるＳＮＲを最大化できる分岐数が存在する。ＦＯＭが最大となる分岐数を選択して分岐線路を構成することで、最大のＳＮＲにて本発明の分岐光線路の特性解析装置及び方法を実施可能である。以上の数値解析を実施することで、本発明は、総センシング長が与えられた場合に（分岐数に制約がない場合に）ＳＮＲを最大化する最適分岐数の設計手法を提供するものである。
【０１０２】
なお、図６の横軸に示される任意の分岐数Ｎは、入射パワーを等分配する２分岐、３分岐、５分岐の光合分波器の組み合わせにて構成することができる。例えば、３分岐の光合分波器は１対２の光合分波器で分波したのち、高出力ポートの出力をさらに１対１の光合分波器で分けることで実現できる。５分岐の光合分波器は２対３の光合分波器で分波したのち、低出力ポートの出力をさらに１対１の光合分波器で分け、高出力ポートに上記３分岐の合分波器を接続することで実現できる。ただし、一般に光合分波器には、理論的な分岐損失（Ｎ分岐に対して光強度がＮ分の１になる）の他、過剰損失と呼ばれる光損失が生じるため、合分波器の挿入損失が各分岐下部心線に接続する光合分波器の出力ポート毎に均一になるように構成することが望ましい。
【０１０３】
なお、この発明は上記実施形態そのままに限定されるものではなく、実施段階ではその要旨を逸脱しない範囲で構成要素を変形して具体化できる。また、上記各実施形態に開示されている複数の構成要素の適宜な組み合せにより種々の発明を形成できる。例えば、各実施形態に示される全構成要素から幾つかの構成要素を削除してもよい。さらに、異なる実施形態に亘る構成要素を適宜組み合せてもよい。
【０１０４】
［付記］
以下は、本実施形態の分岐光線路特性解析装置及びその方法を説明したものである。
【０１０５】
（１）：
光スプリッタに接続される複数の分岐光線路に入射される試験光パルスを発生する試験光パルス発生部と、
前記入射された試験光パルスに由来する戻り光パルスを受光して電気信号を発生する受光部と、
前記電気信号を処理する信号処理部とを具備し、
前記試験光パルス発生部は、
光周波数の等しいＸ個（Ｘは２以上の自然数）のプローブ光パルスを発生するプローブ光パルス発生部と、
前記複数の分岐光線路の遠端でそれぞれ反射された前記Ｘ個のプローブ光パルスに前記分岐光線路内で対向伝搬し衝突したプローブ光パルスとインタラクションするポンプ光パルスを発生するポンプ光パルス発生部と、
前記受光された戻り光パルス同士の干渉を防止すべく前記Ｘ個のプローブ光パルスのパルス間隔を制御する制御部とを備え、
前記信号処理部は、
前記電気信号を解析して得られる前記戻り光パルスのブリルアン（Ｂｒｉｌｌｏｕｉｎ）利得スペクトル分布に基づいて前記複数の分岐光線路の特性を個別に解析する特性解析装置において、受光するブリルアン利得が最大となる分岐数を予め計算により求めることを特徴とする、分岐光線路の特性解析装置。
【０１０６】
（２）：
前記ブリルアン利得が最大となる分岐数は、分岐数Ｎ、ブリルアン利得係数ｇＢ、光ファイバのコアの実効断面積Ａｅｆｆ、プローブ光とポンプ光の相関長をｌｃ、分岐下部心線の長さＬ、光ファイバの光損失係数αとし、被測定線路の分岐下部線路遠端でのプローブ光パワーＰｒ（Ｌ）とポンプ光パワーＰｕ（Ｌ）について結合方程式を解くことにより求め、式（Ｃ１）が最大となる分岐数Ｎを求めることで決定することを特徴とする、上記（１）に記載の特性解析装置。
【０１０７】
（３）：
分岐光線路の線長差はΔＬで等間隔とし、前記プローブ光パルスのパルス幅をτｐｒｏｂｅとし、光ファイバ中の光速をνとしたとき、２ΔＬ＝τｐｒｏｂｅνを満たし、分岐数Ｎとすると、最長と最短の分岐線路の長さの差ΔＬｍａｘは、ΔＬｍａｘ＝ΔＬ（Ｎ−１）を満たす構成の分岐光線路に接続し、
前記制御部は、Ｔｓ＝τｐｒｏｂｅ＋２ΔＬｍａｘ／νを満たすべく前記Ｘ個のプローブ光パルスのパルス間隔を制御する、上記（１）に記載の分岐光線路の特性解析装置。
【０１０８】
（４）：
前記信号処理部は、
前記プローブ光パルスと前記ポンプ光パルスとの光周波数差を変化させ、
前記戻り光パルスの光強度から得られる距離毎のブリルアン利得スペクトルの最大値を取るブリルアン周波数シフト量に基づいて、当該ブリルアン周波数シフト量の変化に寄与する物理量の特性分布を前記複数の分岐光線路毎に取得する、上記（１）に記載の分岐光線路の特性解析装置。
【０１０９】
（５）：
光スプリッタに接続される複数の分岐光線路に試験光パルスを入射し前記試験光パルスに由来する戻り光パルスを受光する、分岐光線路の特性解析方法において、
光周波数の等しいＸ個（Ｘは２以上の自然数）のプローブ光パルスを発生することと、
前記受光された戻り光パルス同士の干渉を防止すべく前記Ｘ個のプローブ光パルスのパルス間隔を制御することと、
前記Ｘ個のプローブ光パルスを入射することと、
前記複数の分岐光線路の遠端でそれぞれ反射された前記Ｘ個のプローブ光パルスに前記分岐光線路内で対向伝搬し衝突したプローブ光パルスとインタラクションするポンプ光パルスを発生することと、
前記Ｘ個のプローブ光パルスが入射された後に前記ポンプ光パルスを入射することと、
前記受光された戻り光パルスを電気信号に変換することと、
前記電気信号を解析して得られる前記戻り光パルスのブリルアン（Ｂｒｉｌｌｏｕｉｎ）利得スペクトル分布に基づいて前記複数の分岐光線路の特性を個別に解析することとを具備する特性解析方法において、受光するブリルアン利得が最大となる分岐数を予め計算により求めることを特徴とする、分岐光線路の特性解析方法。
【０１１０】
（６）：
前記ブリルアン利得が最大となる分岐数は、分岐数Ｎ、ブリルアン利得係数ｇＢ、光ファイバのコアの実効断面積Ａｅｆｆ、プローブ光とポンプ光の相関長をｌｃ、分岐下部心線の長さＬ、光ファイバの光損失係数αとし、被測定線路の分岐下部線路遠端でのプローブ光パワーＰｒ（Ｌ）とポンプ光パワーＰｕ（Ｌ）について結合方程式を解くことにより求め、式（Ｃ１）が最大となる分岐数Ｎを求めることで決定することを特徴とする、上記（５）に記載の特性解析方法。
【０１１１】
（７）：
分岐光線路の線長差はΔＬで等間隔とし、前記プローブ光パルスのパルス幅をτｐｒｏｂｅとし、光ファイバ中の光速をνとしたとき、２ΔＬ＝τｐｒｏｂｅνを満たし、分岐数Ｎとすると、最長と最短の分岐線路の長さの差ΔＬｍａｘは、ΔＬｍａｘ＝ΔＬ（Ｎ−１）を満たす構成の分岐光線路に接続し、
前記制御部は、Ｔｓ＝τｐｒｏｂｅ＋２ΔＬｍａｘ／νを満たすべく前記Ｘ個のプローブ光パルスのパルス間隔を制御する、上記（５）に記載の分岐光線路の特性解析方法。
【０１１２】
（８）：
前記解析することは、
前記プローブ光パルスと前記ポンプ光パルスとの光周波数差を変化させ、
前記戻り光パルスの光強度から得られる距離毎のブリルアン利得スペクトルの最大値を取るブリルアン周波数シフト量に基づいて、当該ブリルアン周波数シフト量の変化に寄与する物理量の特性分布を前記複数の分岐光線路毎に取得する、上記（５）に記載の分岐光線路の特性解析方法。
【０１１３】
（効果）
本発明によれば、既存の手法と比較して、単位時間当たりＸ倍の速度で測定点の情報を得ることができ、測定時間をＸ分の１に短縮することが可能な分岐光線路の特性解析装置および分岐光線路の特性解析方法を提供できる。さらに、当該分岐光線路の特性解析装置および分岐光線路の特性解析方法において測定性能を最大化する分岐光線路の設計方法を提供することが可能になる。
【符号の説明】
【０１１４】
１０：光源
１１：分岐素子
１２：光周波数変更手段
１３：光パルス化手段
１４：光パルス化手段
１５、１６：入射時間制御手段
１７：正弦波発生器
１８、１９：光増幅器
２０：合波素子
２１：サーキュレータ
２２：光スプリッタ
２３：分岐光線路
２４：光反射フィルタ
２５：光フィルタ
２６：光受信器
２７：Ａ／Ｄ変換器
２８：演算処理装置
【図１】
2017026512.tif 000023
【図２】
2017026512.tif 000024
【図３】
2017026512.tif 000025
【図５】
2017026512.tif 000026
【図６】
2017026512.tif 000027
【図４】
2017026512.tif 000028