(57)【要約】【課題】光ファイバ素線の切り出しや加工をすることなく音響モード伝搬速度を測定できる音響モード伝搬速度測定方法及び音響モード伝搬速度測定装置を提供することを目的とする。【解決手段】本発明に係る音響モード伝搬速度測定方法は、光ファイバ中で発生するブリルアン散乱光の周波数シフトスペクトルを取得し、前記周波数シフトスペクトルに対してガウス関数でフィッティングし、前記ガウス関数でフィッティングした曲線からスペクトル半値幅ｗを取得し、スペクトル半値幅ｗと音響モードの伝搬速度ＶＡとの一次関数に取得した前記スペクトル半値幅ｗを代入して前記光ファイバの音響モードの伝搬速度ＶＡを計算することを特徴とする。【選択図】図５
2020143935.tif 000002
【特許請求の範囲】
【請求項１】
光ファイバ中で発生するブリルアン散乱光の周波数シフトスペクトルを取得し、
前記周波数シフトスペクトルに対してガウス関数でフィッティングし、
前記ガウス関数でフィッティングした曲線からスペクトル半値幅ｗを取得し、
スペクトル半値幅ｗと音響モードの伝搬速度ＶＡとの一次関数に取得した前記スペクトル半値幅ｗを代入して前記光ファイバの音響モードの伝搬速度ＶＡを計算すること
を特徴とする音響モード伝送速度測定方法。
【請求項２】
前記一次関数が数Ｃ１であることを特徴とする請求項１に記載の音響モード伝送速度測定方法。
【数Ｃ１】
2020143935.tif 000008
である。
【請求項３】
光ファイバ中で発生するブリルアン散乱光の周波数シフトスペクトルを取得し、
前記周波数シフトスペクトルに対して数Ｃ２の関数でフィッティングし、
前記関数でフィッティングした曲線から音響モードの寿命τを取得し、
音響モードの寿命τと音響モードの伝搬速度ＶＡとの一次関数に取得した前記音響モードの寿命τを代入して前記光ファイバの音響モードの伝搬速度ＶＡを計算すること
を特徴とする音響モード伝送速度測定方法。
【数Ｃ２】
2020143935.tif 000009
ここで、Ａ、Ｇはフィッティングパラメータであり、Ｓ（ｆ）はスペクトル強度、ｆは周波数、Γは音響モードの減衰率を示し、音響モードの寿命τはτ＝１／Γである。
【請求項４】
前記一次関数が数Ｃ３であることを特徴とする請求項３に記載の音響モード伝送速度測定方法。
【数Ｃ３】
2020143935.tif 000010
である。
【請求項５】
光ファイバ中で発生するブリルアン散乱光の周波数シフトスペクトルを取得する受信部と、
前記周波数シフトスペクトルに対してガウス関数でフィッティングし、前記ガウス関数でフィッティングした曲線からスペクトル半値幅ｗを取得する解析部と、
スペクトル半値幅ｗと音響モードの伝搬速度ＶＡとの一次関数に取得した前記スペクトル半値幅ｗを代入して前記光ファイバの音響モードの伝搬速度ＶＡを計算する計算部と、
を備えることを特徴とする音響モード伝送速度測定装置。
【請求項６】
前記一次関数が数Ｃ１であることを特徴とする請求項５に記載の音響モード伝送速度測定装置。
【数Ｃ１】
2020143935.tif 000011
である。
【請求項７】
光ファイバ中で発生するブリルアン散乱光の周波数シフトスペクトルを取得する受信部と、
前記周波数シフトスペクトルに対して数Ｃ２の関数でフィッティングし、前記関数でフィッティングした曲線から音響モードの寿命τを取得する解析部と、
音響モードの寿命τと音響モードの伝搬速度ＶＡとの一次関数に取得した前記音響モードの寿命τを代入して前記光ファイバの音響モードの伝搬速度ＶＡを計算する計算部と、
を備えることを特徴とする音響モード伝送速度測定装置。
【数Ｃ２】
2020143935.tif 000012
ここで、Ａ、Ｇはフィッティングパラメータであり、Ｓ（ｆ）はスペクトル強度、ｆは周波数、Γは音響モードの減衰率を示し、音響モードの寿命τはτ＝１／Γである。
【請求項８】
前記一次関数が数Ｃ３であることを特徴とする請求項７に記載の音響モード伝送速度測定装置。
【数Ｃ３】
2020143935.tif 000013
である。
【発明の詳細な説明】
【技術分野】
【０００１】
本開示は、光ファイバ素線中および敷設済みケーブル内の光ファイバ中を伝搬する音響モード伝搬特性を測定する方法及び装置に関するものである。
【背景技術】
【０００２】
現在の光ネットワークでは伝送遅延時間の低減に対する要求が顕在化し、特に長距離通信線路では光ファイバによる伝搬遅延時間が支配的であるため、光ファイバの群遅延時間に対する注目が高まっている。従来の測定では光パルス送信機から入射された光を、光ファイバ通過後に受信機で検出し伝搬時間を計測することによって群遅延時間を求めるインパルス応答法が利用されていた。一方、最近では、後方ブリルアン散乱を用いた手法が提案されており、従来手法と異なり光ファイバの片端測定が可能であるため、敷設済みの光ファイバに対しても応用が期待されている（非特許文献１）。
【０００３】
上記の後方ブリルアン散乱を用いた測定手法では光の散乱に寄与する縦波の音響モードについて光ファイバ中を伝搬する速度を特定する必要がある。光の散乱に寄与する音響モードの伝搬速度を特定する手法には、
（１）非特許文献２の光ファイバのコア屈折率測定値から経験的に求める方法、
（２）非特許文献３で紹介されている超音波パルスの照射による方法
の２つがある。
【０００４】
手法（１）は光ファイバの屈折率分布を測定する事によって実施可能である。手法（２）は被測定光ファイバを超音波パルス発射装置に光ファイバ断面が向くように台座に設置し、被測定光ファイバを通過して台座で反射した超音波パルスと、光ファイバを通過せず台座で反射した超音波パルスとの受信時間差を計算することによって実施可能である。
【先行技術文献】
【非特許文献】
【０００５】
【非特許文献１】寒河江、松井、辻川、中島“誘導ブリルアン散乱による光ファイバ群遅延時間測定法に関する検討，”電子情報通信学会ソサイエティ大会，  Ｂ−１３−１３，  ９月，  ２０１８
【非特許文献２】Ｙ．  Ｋｏｙａｍａｄａ，  Ｓ．  Ｓａｔｏ，  Ｓ．  Ｎａｋａｍｕｒａ，  Ｈ．  Ｓｏｔｏｂａｙａｓｈｉ  ａｎｄ  Ｗ．  Ｃｈｕｊｏ，  “Ｓｉｍｕｌａｔｉｎｇ  ａｎｄ  Ｄｅｓｉｇｎｉｎｇ  Ｂｒｉｌｌｏｕｉｎ  Ｇａｉｎ  Ｓｐｅｃｔｒｕｍ  ｉｎ  Ｓｉｎｇｌｅ−Ｍｏｄｅ  Ｆｉｂｅｒｓ，”  Ｊ．  Ｌｉｇｈｔｗａｖｅ  Ｔｅｃｈｏｌ．，  Ｖｏｌ．  ２２，  Ｎｏ．  ２，  ｐｐ．  ６３１−６３９，  Ｆｅｂｒｕａｒｙ，  ２００４．
【非特許文献３】Ｎ．  Ｈａｙａｓｈｉ，  Ｙ．  Ｍｉｚｕｎｏ，  Ｄ．  Ｋｏｙａｍａ  ａｎｄ  Ｋ．  Ｎａｋａｍｕｒａ，  “Ｄｅｐｅｎｄｅｎｃｅ  ｏｆ  Ｂｒｉｌｌｏｕｉｎ  Ｆｒｅｑｕｅｎｃｙ  Ｓｈｉｆｔ  ｏｎ  Ｔｅｍｐｅｒａｔｕｒｅ  ａｎｄ  Ｓｔｒａｉｎ  ｉｎ  Ｐｏｌｙ（ｍｅｔｈｙｌ  ｍｅｔｈａｃｒｙｌａｔｅ）−Ｂａｓｅｄ  Ｐｏｌｙｍｅｒ  Ｏｐｔｉｃａｌ  Ｆｉｂｅｒｓ  Ｅｓｔｉｍａｔｅｄ  ｂｙ  Ａｃｏｕｓｔｉｃ  Ｖｅｌｏｃｉｔｙ  Ｍｅａｓｕｒｅｍｅｎｔ，”  Ａｐｐｌ．  Ｐｈｙｓ．  Ｅｘｐ．  ５，  ０３２５０２，  Ｍａｒｃｈ，  ２０１２．
【非特許文献４】Ａ．  Ｙｅｎｉａｙ，  Ｊ−Ｍ．  Ｄｅｌａｖａｕｘ  ａｎｄ  Ｊ．  Ｔｏｕｌｏｕｓｅ，“Ｓｐｏｎｔａｎｅｏｕｓ  ａｎｄ  Ｓｔｉｍｕｌａｔｅｄ  Ｂｒｉｌｌｏｕｉｎ  Ｓｃａｔｔｅｒｉｎｇ  Ｇａｉｎ  Ｓｐｅｃｔｒａ  ｉｎ  Ｏｐｔｉｃａｌ  Ｆｉｂｅｒｓ，”  Ｊ．  Ｌｉｇｈｔｗａｖｅ  Ｔｅｃｈｎｏｌ．，  Ｖｏｌ．  ２０，  Ｎｏ．  ８，  ｐｐ．  １４２５−１４３２，  Ａｕｇｕｓｔ，  ２００２．
【非特許文献５】Ｈ．  Ａ．  Ａｌ−Ａｓａｄｉ，  Ｍ．  Ｈ．  Ａｌ−Ｍａｎｓｏｏｒｉ，  Ｓ．  Ｈｉｔａｍ，  Ｍ．  Ｉ．  Ｓａｒｉｐａｎ  ａｎｄ  Ｍ．  Ａ．  Ｍａｈｄｉｌ  “Ｂｒｉｌｌｏｕｉｎ  Ｌｉｎｅｗｉｄｔｈ  Ｃｈａｒａｃｔｅｒｉｚａｔｉｏｎ  ｉｎ  Ｓｉｎｇｌｅ  Ｍｏｄｅ  Ｌａｒｇｅ  Ｅｆｆｅｃｔｉｖｅ  Ａｒｅａ  Ｆｉｂｅｒ  ｔｈｒｏｕｇｈ  Ｃｏ−Ｐｕｍｐｅｄ  Ｔｅｃｈｎｉｑｕｅ，”  ＩＪＥＣＣＴ  ２０１０，  Ｖｏｌ．  １，  １，  ｐｐ．  １６−２０，  ２０１０．
【発明の概要】
【発明が解決しようとする課題】
【０００６】
上記手法（１）および（２）は、どちらも後方ブリルアン散乱測定の前または後に、音響モードの伝搬速度の測定のために光ファイバ素線の切出しや加工が必要となる。そのため、上記手法（１）および（２）には、敷設済みケーブル内の光ファイバなど、光ファイバ素線の切り出しや加工ができない被測定光ファイバに対して音響モード伝搬速度を測定することが困難という課題があった。そこで、本発明は、上記課題を解決するために、光ファイバ素線の切り出しや加工をすることなく音響モード伝搬速度を測定できる音響モード伝搬速度測定方法及び音響モード伝搬速度測定装置を提供することを目的とする。
【課題を解決するための手段】
【０００７】
上記目的を達成するために、本発明に係る音響モード伝搬速度測定方法及び音響モード伝搬速度測定装置は、試験光入射によって光ファイバ中で発生する後方ブリルアン散乱スペクトルを解析することによって音響モードの伝搬速度を計算することとした。
【０００８】
非特許文献２によると縦波音響波の伝搬速度と音響波の寿命は共に伝搬媒体への不純物添加量に比例することが報告されている。そのため、光ファイバのコアを伝搬する、光の散乱に寄与する縦波の音響モードの伝搬速度は寿命に比例する相関があると考えられる。また、後方ブリルアン散乱スペクトルの半値幅が音響モードの寿命によって決まることを考えると、光の散乱に寄与する縦波の音響モードの伝搬速度は後方ブリルアン散乱スペクトルの半値幅に比例する相関があると考えられる。よって、後方ブリルアン散乱スペクトルから半値幅、もしくは音響モードの寿命を解析することによって、光の散乱に寄与する音響モードの伝搬速度が測定可能である。
【０００９】
具体的には、本発明に係る第１の音響モード伝搬速度測定方法は、
光ファイバ中で発生するブリルアン散乱光の周波数シフトスペクトルを取得し、
前記周波数シフトスペクトルに対してガウス関数でフィッティングし、
前記ガウス関数でフィッティングした曲線からスペクトル半値幅ｗを取得し、
スペクトル半値幅ｗと音響モードの伝搬速度ＶＡとの一次関数に取得した前記スペクトル半値幅ｗを代入して前記光ファイバの音響モードの伝搬速度ＶＡを計算すること
を特徴とする。
【００１０】
また、第１の音響モード伝搬速度測定方法を行う第１の音響モード伝搬速度測定装置は、
光ファイバ中で発生するブリルアン散乱光の周波数シフトスペクトルを取得する受信部と、
前記周波数シフトスペクトルに対してガウス関数でフィッティングし、前記ガウス関数でフィッティングした曲線からスペクトル半値幅ｗを取得する解析部と、
スペクトル半値幅ｗと音響モードの伝搬速度ＶＡとの一次関数に取得した前記スペクトル半値幅ｗを代入して前記光ファイバの音響モードの伝搬速度ＶＡを計算する計算部と、
を備える。
【００１１】
このとき、前記一次関数が数Ｃ１であることを特徴とする。
【数Ｃ１】
2020143935.tif 000003
である。
【００１２】
一方、本発明に係る第２の音響モード伝搬速度測定方法は、
光ファイバ中で発生するブリルアン散乱光の周波数シフトスペクトルを取得し、
前記周波数シフトスペクトルに対して数Ｃ２の関数でフィッティングし、
前記関数でフィッティングした曲線から音響モードの寿命τを取得し、
音響モードの寿命τと音響モードの伝搬速度ＶＡとの一次関数に取得した前記音響モードの寿命τを代入して前記光ファイバの音響モードの伝搬速度ＶＡを計算すること
を特徴とする。
【数Ｃ２】
2020143935.tif 000004
ここで、Ａ、Ｇはフィッティングパラメータであり、Ｓ（ｆ）はスペクトル強度、ｆは周波数、Γは音響モードの減衰率を示し、音響モードの寿命τはτ＝１／Γである。
【００１３】
また、第２の音響モード伝搬速度測定方法を行う第２の音響モード伝搬速度測定装置は、
光ファイバ中で発生するブリルアン散乱光の周波数シフトスペクトルを取得する受信部と、
前記周波数シフトスペクトルに対して数Ｃ２の関数でフィッティングし、前記関数でフィッティングした曲線から音響モードの寿命τを取得する解析部と、
音響モードの寿命τと音響モードの伝搬速度ＶＡとの一次関数に取得した前記音響モードの寿命τを代入して前記光ファイバの音響モードの伝搬速度ＶＡを計算する計算部と、
を備えることを特徴とする。
【００１４】
このとき、前記一次関数が数Ｃ３であることを特徴とする。
【数Ｃ３】
2020143935.tif 000005
である。
【００１５】
なお、上記各発明は、可能な限り組み合わせることができる。
【発明の効果】
【００１６】
本発明は、光ファイバ素線の切り出しや加工をすることなく音響モード伝搬速度を測定できる音響モード伝搬速度測定方法及び音響モード伝搬速度測定装置を提供することができる。
【図面の簡単な説明】
【００１７】
【図１】光ファイバ屈折率分布の測定結果を説明する図である。
【図２】光ファイバの音響波速度分布の計算結果と光の散乱に寄与する音響モードの強度分布の計算結果を説明する図である。
【図３】本発明に係る音響モード伝搬速度測定方法で取得される後方ブリルアン散乱スペクトルの取得例と解析例を説明する図である。
【図４】本発明に係る音響モード伝搬速度測定方法で計算された音響モードの伝搬速度と後方ブリルアン散乱スペクトル線幅との相関を説明する特性図である。
【図５】本発明に係る音響モード伝搬速度測定装置を説明する図である。
【図６】本発明に係る音響モード伝搬速度測定方法で取得される後方ブリルアン散乱スペクトルの取得例と解析例を説明する図である。
【図７】本発明に係る音響モード伝搬速度測定方法で計算された音響モードの伝搬速度と音響モードの寿命との相関を説明する特性図である。
【図８】本発明に係る音響モード伝搬速度測定装置を説明する図である。
【図９】本発明に係る音響モード伝搬速度測定方法で計算した音響モードの伝搬速度と公知の手法で計算した計算値とを比較する表である。
【図１０】本発明に係る音響モード伝搬速度測定方法で計算した音響モードの伝搬速度と公知の手法で計算した計算値とを比較する表である。
【発明を実施するための形態】
【００１８】
添付の図面を参照して本発明の実施形態を説明する。以下に説明する実施形態は本発明の実施例であり、本発明は、以下の実施形態に制限されるものではない。なお、本明細書及び図面において符号が同じ構成要素は、相互に同一のものを示すものとする。
【００１９】
（実施形態１）
本実施形態の音響モード伝送速度測定方法は、後方ブリルアン散乱スペクトルの半値幅解析による光ファイバ片端測定であって、
光ファイバ中で発生するブリルアン散乱光の周波数シフトスペクトルを取得し、
前記周波数シフトスペクトルに対してガウス関数でフィッティングし、
前記ガウス関数でフィッティングした曲線からスペクトル半値幅ｗを取得し、
スペクトル半値幅ｗと音響モードの伝搬速度ＶＡとの一次関数に取得した前記スペクトル半値幅ｗを代入して前記光ファイバの音響モードの伝搬速度ＶＡを計算すること
を特徴とする。
【００２０】
図１にＧｅＯ２添加コアを有する複数の光ファイバについて測定した光ファイバ径方向の屈折率分布を示す。光ファイバ（ａ）は汎用のシングルモード光ファイバ（ＳＭＦ）、（ｂ）は分散シフト光ファイバ（ＤＳＦ）、（ｃ）はフッ素添加によるトレンチ層を有する高ＧｅＯ２添加コア光ファイバ（ＨＧＦ）である。
【００２１】
図２に非特許文献２によって示されている屈折率と音響波速度に関する経験式に基づいて計算した光ファイバ（ａ）〜（ｃ）の光ファイバ径方向の音響波速度分布（左軸、実線）を示す。また、同時に本音響波速度分布に基づいて有限要素法で計算した音響モード強度分布（右軸、破線）で示した。以上の速度分布Ｖｌ（ｒ）と強度分布Ｉ（ｒ）を用いて式（１）の重なり積分により光ファイバ中の実効的な音響モードの伝搬速度ＶＡを計算することができる。各図に光ファイバの音響波速度分布をもとに式（１）で求めたＶＡを示した。
【数１】
2020143935.tif 000006
【００２２】
図３に上記光ファイバ（ａ）〜（ｃ）にて測定した後方ブリルアン散乱スペクトルの測定波形を示す。ここで入射光は波長１．５５μｍの連続光とした。非特許文献４によると試験光強度が十分高強度な場合、後方ブリルアン散乱スペクトルはガウシアン型となることが知られており、図１中の実線はガウス関数によるフィッティング曲線である。スペクトル強度がピーク強度に対して−３ｄＢとなる区間ではいずれの光ファイバの後方ブリルアン散乱スペクトルについても良くフィッティングできる。よってガウス関数のフィッティングパラメータからスペクトル半値幅ｗを解析することが可能である。
【００２３】
図４は、光ファイバ（ａ）〜（ｃ）にて図２にて解析した後方ブリルアン散乱スペクトルの半値幅ｗを横軸とし、式（１）で計算した音響モードの伝搬速度ＶＡを縦軸とした相関図である。半値幅ｗと実効的な音響モードの伝搬速度ＶＡは負の相関を持つことが示されている。エラーバーの範囲内で半値幅ｗは変動し、図中の実線はフィッティング直線の一例であり
［数２］
ＶＡ＝−１３８．９ｗ＋７１５６                                          （２）
の関係式で表される。半値幅ｗの変動を考慮すると、ＶＡとｗの関係式をＶＡ＝Ａ１ｗ＋Ｂ１とするとき
［数３］
−１５９．８≦Ａ１≦−１２４．６
６９９２≦Ｂ１≦７３９７                                                （３）
の範囲の値をとる。
【００２４】
図９の表１に、光ファイバ（ａ）〜（ｃ）にて式（１）で計算した音響モードの伝搬速度と式（２）によりスペクトル半値幅から求めた音響モードの伝搬速度を示した。両差の誤差は１．０％未満であり、式（１）で計算する音響モードの伝搬速度をスペクトル半値幅から式（２）で精度よく計算できることが分かる。非特許文献５によると後方ブリルアン散乱スペクトル半値幅は光ファイバ長に依って変化することが報告されているが、１０ｋｍ以上の光ファイバではその変化が飽和するため、１０ｋｍ以上の長尺の光ファイバでの実施が好ましい。以上の様に後方ブリルアン散乱スペクトル解析によってスペクトル半値幅を求め、光の散乱に寄与する実効的な音響モードの伝搬速度を求めることができる。
【００２５】
図５は、上記の音響モード伝搬速度測定方法を実現する音響モード伝搬速度測定装置３０１を説明する構成例である。音響モード伝搬速度測定装置３０１は、
光ファイバ中で発生するブリルアン散乱光の周波数シフトスペクトルを取得する受信部１１と、
前記周波数シフトスペクトルに対してガウス関数でフィッティングし、前記ガウス関数でフィッティングした曲線からスペクトル半値幅ｗを取得する解析部１２と、
スペクトル半値幅ｗと音響モードの伝搬速度ＶＡとの一次関数に取得した前記スペクトル半値幅ｗを代入して前記光ファイバの音響モードの伝搬速度ＶＡを計算する計算部１３と、
を備える。
【００２６】
後方ブリルアン散乱を発生させる強度を持ったレーザ光源からの試験光が分岐部２１で励起光と参照光に分岐される。励起光が送信部２２から被測定光ファイバＦＵＴに入射されることで後方散乱光が発生する。後方散乱光が光サーキュレータ２４を経由し、合波部２３で参照光と合波される。受信部１１はビート光を受光し、そのスペクトル（ブリルアン散乱光の周波数シフトスペクトル）を測定する。解析部１２はガウス関数を用いて当該周波数シフトスペクトルの半値幅ｗを解析する。計算部１３は式（３）の係数を用いたＶＡとｗの関係式ＶＡ＝Ａ１ｗ＋Ｂ１により音響モード伝搬速度を計算する。
【００２７】
本実施形態の測定手法により、音響モードの伝搬速度の測定のための光ファイバ切出しや加工が不要となり、光ファイバの非破壊検査を実現することができる。
【００２８】
（実施形態２）
本実施形態の音響モード伝送速度測定方法は、後方ブリルアン散乱スペクトルの音響モード寿命解析による光ファイバ片端測定であって、
光ファイバ中で発生するブリルアン散乱光の周波数シフトスペクトルを取得し、
前記周波数シフトスペクトルに対して式（４）の関数でフィッティングし、
前記関数でフィッティングした曲線から音響モードの寿命τを取得し、
音響モードの寿命τと音響モードの伝搬速度ＶＡとの一次関数に取得した前記音響モードの寿命τを代入して前記光ファイバの音響モードの伝搬速度ＶＡを計算すること
を特徴とする。
【００２９】
一般的にブリルアン散乱スペクトル線幅は光の散乱に寄与する音響モードの寿命に依存することが知られている。図６に図１の被測定光ファイバで測定した後方ブリルアン散乱スペクトルに対して、非特許文献５で報告されている音響モード寿命を解析可能な以下の関数でフィッティングした結果を示す。
【数４】
2020143935.tif 000007
ここで、Ｓ（ｆ）は周波数毎のスペクトル強度、Ａ、Ｇはフィッティングパラメータ、ｆは周波数である。Γは音響モードの減衰率を示し、音響モードの寿命τはτ＝１／Γとなる。図６より測定結果が式（４）によって表現できることがわかる。
【００３０】
図７は、図６で解析した音響モードの寿命τを横軸とし、式（１）で計算した音響モードの伝搬速度ＶＡを縦軸としたの相関図である。両パラメータは直線的な相関を持つことが示されており、非特許文献２の経験式から示唆される相関と整合する。エラーバーの範囲で音響モードの寿命τは変動し、図中の実線はフィッティング直線の一例であり
［数５］
ＶＡ＝３６５．４τ＋３８９０                                            （５）
の関係式で表される。音響モードの寿命τの変動を考慮すると、ＶＡとτの関係式をＶＡ＝Ａ２τ＋Ｂ２とするとき
［数６］
３５４．５≦Ａ２≦３７６．９
３８４０≦Ｂ２≦３９３９                                              （６）
の範囲の値をとる。
【００３１】
図１０の表２に、光ファイバ（ａ）〜（ｃ）にて式（１）で計算した音響モードの伝搬速度と式（５）により音響波モードの寿命から求めた音響モードの伝搬速度を示した。両者の誤差は１．０％未満であり、式（５）で計算する音響モードの伝搬速度を後方ブリルアン散乱スペクトル解析によって得られる音響モードの寿命から式（１）で精度よく計算できることが分かる。実施形態１と同様に本発明の測定手法は１０ｋｍ以上の長尺の光ファイバでの実施が好ましい。以上の様に後方ブリルアン散乱スペクトル解析によって音響モード寿命を求め、光の散乱に寄与する音響モードの伝搬速度を求めることができる。
【００３２】
図５の音響モード伝搬速度測定装置３０１の構成で上記の音響モード伝搬速度測定方法を実現することができる。この場合、音響モード伝搬速度測定装置３０１は、解析部１２が前記周波数シフトスペクトルに対して式（４）の関数でフィッティングし、前記関数でフィッティングした曲線から音響モードの寿命τを取得し、計算部１３が音響モードの寿命τと音響モードの伝搬速度ＶＡとの一次関数に取得した前記音響モードの寿命τを代入して前記光ファイバの音響モードの伝搬速度ＶＡを計算する。
【００３３】
解析部１２は式（４）を用いて受信部１１が受光したビート光のスペクトル（ブリルアン散乱光の周波数シフトスペクトル）から音響モードの寿命τを計算する。計算部１３は式（６）の係数を用いたＶＡとτの関係式ＶＡ＝Ａ２τ＋Ｂ２により音響モード伝搬速度を計算する。
【００３４】
本実施形態の測定手法により、音響モードの伝搬速度の測定のための光ファイバ切出しや加工が不要となり、光ファイバの非破壊検査を実現することができる。
【００３５】
（実施形態３）
図８は、本実施形態の音響モード伝搬速度測定方法を実現する音響モード伝搬速度測定装置３０１を説明する構成例である。本実施形態では、ＢＯＴＤＡを用いた音響モードの伝搬速度の分布測定を説明する。本実施形態の測定系と実施形態１及び２の測定系との違いは、本実施形態の測定系が励起光をパルス化するパルス化部２５を備えていることである。
【００３６】
後方ブリルアン散乱を発生させる強度を持ったレーザ光源からの試験光が分岐部２１で励起光と参照光に分岐される。励起光はパルス化部２５でパルス化され、送信部２２から被測定光ファイバＦＵＴに入射される。当該励起光による後方散乱光は、光サーキュレータ２４を経由し、合波部２３で参照光と合波される。音響モード伝搬速度測定装置３０１は、実施形態１または実施形態２で説明したスペクトル解析を行う。励起光がパルス化されているため、光の散乱に寄与する音響モードの伝搬速度を被測定光ファイバに対して分布的に求めることができる。
【産業上の利用可能性】
【００３７】
本発明は、光ファイバの特性評価また光ファイバによる計測技術に用いることができる。
また、本発明に係る音響モード伝搬速度測定装置は、コンピュータとプログラムによっても実現でき、プログラムを記録媒体に記録することも、ネットワークを通して提供することも可能である。
【符号の説明】
【００３８】
１１：受信部
１２：解析部
１３：計算部
２１：分岐部
２２：送信部
２３：合波部
２４：光サーキュレータ
２５：パルス化部
３０１：音響モード伝搬速度測定装置
【図１】
2020143935.tif 000014
【図２】
2020143935.tif 000015
【図３】
2020143935.tif 000016
【図４】
2020143935.tif 000017
【図５】
2020143935.tif 000018
【図６】
2020143935.tif 000019
【図７】
2020143935.tif 000020
【図８】
2020143935.tif 000021
【図９】
2020143935.tif 000022
【図１０】
2020143935.tif 000023