特許 6796043 光反射測定装置及びその方法

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】6796043

(24)【登録日】2020年11月17日

(45)【発行日】2020年12月2日

(54)【発明の名称】光反射測定装置及びその方法

(51)【国際特許分類】

   G01M 11/00 20060101AFI20201119BHJP

   H04B 10/071 20130101ALI20201119BHJP

【ＦＩ】

   G01M11/00 R

   H04B10/071

【請求項の数】6

【全頁数】15

(21)【出願番号】特願2017-177149(P2017-177149)

(22)【出願日】2017年9月14日

(65)【公開番号】特開2019-52938(P2019-52938A)

(43)【公開日】2019年4月4日

【審査請求日】2019年9月3日

(73)【特許権者】

【識別番号】000004226

【氏名又は名称】日本電信電話株式会社

(74)【代理人】

【識別番号】110001243

【氏名又は名称】特許業務法人 谷・阿部特許事務所

(72)【発明者】

【氏名】大野 槙悟

(72)【発明者】

【氏名】飯田 大輔

(72)【発明者】

【氏名】戸毛 邦弘

(72)【発明者】

【氏名】真鍋 哲也

【審査官】 藤田 都志行

(56)【参考文献】

【文献】 特表２００８−５２３３９６（ＪＰ，Ａ）

【文献】 特表２０１５−５００４８３（ＪＰ，Ａ）

【文献】 国際公開第２００８／１０５３２２（ＷＯ，Ａ１）

【文献】 特開２０１０−２７１１３７（ＪＰ，Ａ）

【文献】 米国特許出願公開第２０１７／００５２０９１（ＵＳ，Ａ１）

(58)【調査した分野】（Int.Cl.，ＤＢ名）

Ｇ０１Ｍ １１／００−１１／０２

Ｈ０４Ｂ １０／０７１

ＪＳＴＰｌｕｓ／ＪＳＴ７５８０（ＪＤｒｅａｍＩＩＩ）

(57)【特許請求の範囲】

【請求項1】

周波数掃引した連続光を出射する光源と、
前記連続光を第１の光経路と第２の光経路とに分波して前記第１の光経路伝搬後の連続光と前記第２の光経路伝搬後の連続光とを合波してビート信号を出力する干渉計であって、前記第１の光経路は遅延時間τ_Loopを与える周回光路を含み、前記第１の光経路伝搬後の連続光は前記周回光路伝搬前の連続光および前記周回光路伝搬後の連続光を被測定光ファイバに入射して前記被測定光ファイバで生じた反射光又は散乱光である、前記干渉計と、
前記ビート信号を電気信号に変換する光検出手段と、
前記電気信号を用いて反射光又は散乱光の光振幅分布波形を算出する演算処理手段と、
を備え、前記演算処理手段は、
前記ビート信号を時間領域でτ_Loop単位に分割してフーリエ変換し、前記τ_Loop単位毎の複素スペクトルを複数算出し、
複数の前記複素スペクトルのそれぞれの最も高周波側に現れる反射光又は散乱光振幅分布波形の任意地点の位相成分を抽出して前記反射光又は散乱光の位相の時間変化を算出し、
前記位相の時間変化に基づき算出した時間列で前記ビート信号をリサンプリングし、
前記リサンプリングされたビート信号をフーリエ変換し、前記被測定光ファイバの位相揺らぎが補償された反射光又は散乱光の光振幅分布波形を算出することを特徴とする光反射測定装置。

【請求項2】

周波数掃引した連続光を出射する光源と、
前記連続光を第１の光経路と第２の光経路とに分波して前記第１の光経路伝搬後の連続光と前記第２の光経路伝搬後の連続光とを合波してビート信号を出力する干渉計であって、前記第１の光経路は遅延時間τ_Loopを与える周回光路を含み、前記第１の光経路伝搬後の連続光は前記周回光路伝搬前の連続光および前記周回光路伝搬後の連続光を被測定光ファイバに入射して前記被測定光ファイバで生じた反射光又は散乱光である、前記干渉計と、
前記ビート信号を電気信号に変換する光検出手段と、
前記電気信号を用いて反射光又は散乱光の光振幅分布波形を算出する演算処理手段と、
を備え、前記演算処理手段は、
前記ビート信号を時間領域でτ_Loop単位に分割してフーリエ変換し、前記τ_Loop単位毎の複素スペクトルを複数算出し、
複数の前記複素スペクトルのそれぞれの最も高周波側に現れる反射光又は散乱光振幅分布波形の任意地点の位相成分を抽出して前記反射光又は散乱光の位相の時間変化を算出し、
前記位相の時間変化から前記被測定光ファイバに生じた振動の周波数を解析することを特徴とする光反射測定装置。

【請求項3】

前記第１の光経路は、前記周回光路と前記被測定光ファイバと間に位相参照用ダミー光ファイバを含み、
前記演算処理手段は、前記位相の時間変化として、前記被測定光ファイバの反射光又は散乱光と前記位相参照用ダミー光ファイバの反射光又は散乱光との位相差の時間変化を用いることを特徴とする請求項１又は２に記載の光反射測定装置。

【請求項4】

光源から周波数掃引した連続光を出射する第１ステップと、
前記連続光を第１の光経路と第２の光経路とに分波して前記第１の光経路伝搬後の連続光と前記第２の光経路伝搬後の連続光とを合波してビート信号を出力する第２ステップであって、前記第１の光経路は遅延時間τ_Loopを与える周回光路を含み、前記第１の光経路伝搬後の連続光は前記周回光路伝搬前の連続光および前記周回光路伝搬後の連続光を被測定光ファイバに入射して前記被測定光ファイバで生じた反射光又は散乱光である、前記ビート信号を出力する前記第２ステップと、
前記ビート信号を電気信号に変換する第３ステップと、
前記電気信号を用いて反射光又は散乱光の光振幅分布波形を算出する第４ステップと、
を有し、前記第４ステップは、
前記ビート信号を時間領域でτ_Loop単位に分割してフーリエ変換し、前記τ_Loop単位毎の複素スペクトルを複数算出する第４１ステップと、
複数の前記複素スペクトルのそれぞれの最も高周波側に現れる反射光又は散乱光振幅分布波形の任意地点の位相成分を抽出して前記反射光又は散乱光の位相の時間変化を算出する第４２ステップと、
前記位相の時間変化に基づき算出した時間列で前記ビート信号をリサンプリングする第４３ステップと、
前記リサンプリングされたビート信号をフーリエ変換し、前記被測定光ファイバの位相揺らぎが補償された反射光又は散乱光の光振幅分布波形を算出する第４４ステップと、
を有することを特徴とする光反射測定方法。

【請求項5】

光源から周波数掃引した連続光を出射する第１ステップと、
前記連続光を第１の光経路と第２の光経路とに分波して前記第１の光経路伝搬後の連続光と前記第２の光経路伝搬後の連続光とを合波してビート信号を出力する第２ステップであって、前記第１の光経路は遅延時間τ_Loopを与える周回光路を含み、前記第１の光経路伝搬後の連続光は前記周回光路伝搬前の連続光および前記周回光路伝搬後の連続光を被測定光ファイバに入射して前記被測定光ファイバで生じた反射光又は散乱光である、前記ビート信号を出力する前記第２ステップと、
前記ビート信号を電気信号に変換する第３ステップと、
前記電気信号を用いて反射光又は散乱光の光振幅分布波形を算出する第４ステップと、
を有し、前記第４ステップは、
前記ビート信号を時間領域でτ_Loop単位に分割してフーリエ変換し、前記τ_Loop単位毎の複素スペクトルを複数算出する第４１ステップと、
複数の前記複素スペクトルのそれぞれの最も高周波側に現れる反射光又は散乱光振幅分布波形の任意地点の位相成分を抽出して前記反射光又は散乱光の位相の時間変化を算出する第４２ステップと、
前記位相の時間変化から前記被測定光ファイバに生じた振動の周波数を解析する第４５ステップと、
を有することを特徴とする光反射測定方法。

【請求項6】

前記第１の光経路は、前記周回光路と前記被測定光ファイバと間に位相参照用ダミー光ファイバを含み、
前記第４２ステップでは、前記位相の時間変化として、前記被測定光ファイバの反射光又は散乱光と前記位相参照用ダミー光ファイバの反射光又は散乱光との位相差の時間変化を用いることを特徴とする請求項４または５記載の光反射測定方法。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、光ファイバ中の反射光もしくは後方レイリー散乱光振幅分布を測定する光反射測定技術に関する。

【背景技術】

【0002】

光ファイバの損失分布測定や光ファイバセンシングに用いられる光反射測定法の１つに光周波数領域反射測定法（ＯＦＤＲ）がある（例えば非特許文献１参照）。ＯＦＤＲは、周波数掃引した連続光を２分岐し、そのうち一方を試験光として被測定ファイバに入射し、被測定ファイバからの反射光もしくは後方散乱光と分岐した他方の連続光とを合波して得られるビート信号の周波数領域から反射（散乱）地点を解析する方法である。ＯＦＤＲにおける試験光をＥ（ｔ）、ビート信号をＩ（ｔ）とするとそれぞれ以下のように表される。

【0003】

【数1】

【0004】

【数2】

【0005】

ここでＥ₀、ν₀、γはそれぞれ試験光の光振幅、初期光周波数、周波数掃引速度、ｒ_FUTとτ_FUTは、被測定ファイバの反射（散乱）地点の反射率と光往復伝搬時間、θ_FUT（ｔ）は被測定ファイバに加わる振動等の外乱による位相雑音である。なお、ここでは簡単のため光源由来の位相雑音は無視している。また、τ_FUTが十分小さいことからτ_FUTの２乗項は無視している。

【0006】

式（２）からわかるように、位相雑音θ_FUT（ｔ）を無視するとＩ（ｔ）は周波数γτ_FUTの正弦波であるため、Ｉ（ｔ）のフーリエ変換により反射（散乱）光の反射位置と振幅が解析される。

【先行技術文献】

【非特許文献】

【0007】

【非特許文献1】W.Eickhohh et al, “Optical frequency domain reflectometry in single-mode fiber,” Applied Physics Letters 39(9), pp.693-695 (1981).

【非特許文献2】S.Kreger et al, “Distributed Rayleigh scatter dynamic strain sensing above the scan rate with optical frequency domain reflectometry,” Proc. SPIE 9480, Fiber Optic Sensors and Applications XII, 948006 (2015).

【発明の概要】

【発明が解決しようとする課題】

【0008】

しかしながら、前述では振動等の外乱による位相雑音θ_FUT（ｔ）を無視できることとしているが、実際の測定場面では外乱によりビート信号のスペクトルが広がるため、ＯＦＤＲの空間分解能が劣化するという課題がある。

【0009】

本発明は上記事情を鑑みてなされたものであり、その目的は、被測定ファイバに振動等の外乱が加わる環境下において、外乱を解析する、または補償して高空間分解能を実現する光反射測定装置及びその方法を提供することにある。

【課題を解決するための手段】

【0010】

上記の課題を解決するために、本発明は、光反射測定装置であって、周波数掃引した連続光を出射する光源と、前記連続光を第１の光経路と第２の光経路とに分波して前記第１の光経路伝搬後の連続光と前記第２の光経路伝搬後の連続光とを合波してビート信号を出力する干渉計であって、前記第１の光経路は遅延時間τ_Loopを与える周回光路を含み、前記第１の光経路伝搬後の連続光は前記周回光路伝搬前の連続光および前記周回光路伝搬後の連続光を被測定光ファイバに入射して前記被測定光ファイバで生じた反射光又は散乱光である、前記干渉計と、前記ビート信号を電気信号に変換する光検出手段と、前記電気信号を用いて反射光又は散乱光の光振幅分布波形を算出する演算処理手段と、を備え、前記演算処理手段は、前記ビート信号を時間領域でτ_Loop単位に分割してフーリエ変換し、前記τ_Loop単位毎の複素スペクトルを複数算出し、複数の前記複素スペクトルのそれぞれの最も高周波側に現れる反射光又は散乱光振幅分布波形の任意地点の位相成分を抽出して前記反射光又は散乱光の位相の時間変化を算出し、前記位相の時間変化に基づき算出した時間列で前記ビート信号をリサンプリングし、前記リサンプリングされたビート信号をフーリエ変換し、前記被測定光ファイバの位相揺らぎが補償された反射光又は散乱光の光振幅分布波形を算出することを特徴とする光反射測定装置。

【0011】

請求項２に記載の発明は、光反射測定装置であって、周波数掃引した連続光を出射する光源と、前記連続光を第１の光経路と第２の光経路とに分波して前記第１の光経路伝搬後の連続光と前記第２の光経路伝搬後の連続光とを合波してビート信号を出力する干渉計であって、前記第１の光経路は遅延時間τ_Loopを与える周回光路を含み、前記第１の光経路伝搬後の連続光は前記周回光路伝搬前の連続光および前記周回光路伝搬後の連続光を被測定光ファイバに入射して前記被測定光ファイバで生じた反射光又は散乱光である、前記干渉計と、前記ビート信号を電気信号に変換する光検出手段と、前記電気信号を用いて反射光又は散乱光の光振幅分布波形を算出する演算処理手段と、を備え、前記演算処理手段は、前記ビート信号を時間領域でτ_Loop単位に分割してフーリエ変換し、前記τ_Loop単位毎の複素スペクトルを複数算出し、複数の前記複素スペクトルのそれぞれの最も高周波側に現れる反射光又は散乱光振幅分布波形の任意地点の位相成分を抽出して前記反射光又は散乱光の位相の時間変化を算出し、前記位相の時間変化から前記被測定光ファイバに生じた振動の周波数を解析することを特徴とする。

【0012】

請求項３に記載の発明は、請求項１又は２に記載の光反射測定装置であって、前記第１の光経路は、前記周回光路と前記被測定光ファイバと間に位相参照用ダミー光ファイバを含み、前記演算処理手段は、前記位相の時間変化として、前記被測定光ファイバの反射光又は散乱光と前記位相参照用ダミー光ファイバの反射光又は散乱光との位相差の時間変化を用いることを特徴とする。

【0013】

請求項４に記載の発明は、光反射測定方法であって、光源から周波数掃引した連続光を出射する第１ステップと、前記連続光を第１の光経路と第２の光経路とに分波して前記第１の光経路伝搬後の連続光と前記第２の光経路伝搬後の連続光とを合波してビート信号を出力する第２ステップであって、前記第１の光経路は遅延時間τ_Loopを与える周回光路を含み、前記第１の光経路伝搬後の連続光は前記周回光路伝搬前の連続光および前記周回光路伝搬後の連続光を被測定光ファイバに入射して前記被測定光ファイバで生じた反射光又は散乱光である、前記ビート信号を出力する前記第２ステップと、前記ビート信号を電気信号に変換する第３ステップと、前記電気信号を用いて反射光又は散乱光の光振幅分布波形を算出する第４ステップと、を有し、前記第４ステップは、前記ビート信号を時間領域でτ_Loop単位に分割してフーリエ変換し、前記τ_Loop単位毎の複素スペクトルを複数算出する第４１ステップと、複数の前記複素スペクトルのそれぞれの最も高周波側に現れる反射光又は散乱光振幅分布波形の任意地点の位相成分を抽出して前記反射光又は散乱光の位相の時間変化を算出する第４２ステップと、前記位相の時間変化に基づき算出した時間列で前記ビート信号をリサンプリングする第４３ステップと、前記リサンプリングされたビート信号をフーリエ変換し、前記被測定光ファイバの位相揺らぎが補償された反射光又は散乱光の光振幅分布波形を算出する第４４ステップと、を有することを特徴とする。

【0014】

請求項５に記載の発明は、光反射測定方法であって、光源から周波数掃引した連続光を出射する第１ステップと、前記連続光を第１の光経路と第２の光経路とに分波して前記第１の光経路伝搬後の連続光と前記第２の光経路伝搬後の連続光とを合波してビート信号を出力する第２ステップであって、前記第１の光経路は遅延時間τ_Loopを与える周回光路を含み、前記第１の光経路伝搬後の連続光は前記周回光路伝搬前の連続光および前記周回光路伝搬後の連続光を被測定光ファイバに入射して前記被測定光ファイバで生じた反射光又は散乱光である、前記ビート信号を出力する前記第２ステップと、前記ビート信号を電気信号に変換する第３ステップと、前記電気信号を用いて反射光又は散乱光の光振幅分布波形を算出する第４ステップと、を有し、前記第４ステップは、前記ビート信号を時間領域でτ_Loop単位に分割してフーリエ変換し、前記τ_Loop単位毎の複素スペクトルを複数算出する第４１ステップと、複数の前記複素スペクトルのそれぞれの最も高周波側に現れる反射光又は散乱光振幅分布波形の任意地点の位相成分を抽出して前記反射光又は散乱光の位相の時間変化を算出する第４２ステップと、前記位相の時間変化から前記被測定光ファイバに生じた振動の周波数を解析する第４５ステップと、を有することを特徴とする。

【0015】

請求項６に記載の発明は、請求項４または５記載の光反射測定方法において、前記第１の光経路は、前記周回光路と前記被測定光ファイバと間に位相参照用ダミー光ファイバを含み、前記第４２ステップでは、前記位相の時間変化として、前記被測定光ファイバの反射光又は散乱光と前記位相参照用ダミー光ファイバの反射光又は散乱光との位相差の時間変化を用いることを特徴とする。

【発明の効果】

【0016】

本発明により、振動等の外乱が加わる環境下でも高分解能で損失分布測定や光ファイバセンシング等のＯＦＤＲを用いた応用を実現できる。さらに、本発明で得られる反射（散乱）光の位相変化をフーリエ変換して周波数解析することで、被測定ファイバに生じた振動の周波数を解析できるため、光ファイバ振動センシングとしての応用も可能である。ＯＦＤＲを用いた従来の振動センシング手法（例えば非特許文献２）では事前に無振動状態で取得した参照データを必要とするが、本発明による手法では参照データが不要であるため、従来法に比べて容易に振動センシングを実施できる。

【図面の簡単な説明】

【0017】

【図1】本発明における反射（散乱）光の位相解析方法の概念図でり、（ａ）はビート信号を示す図であり、（ｂ）はビート信号をフーリエ変換して得られる複素スペクトルを示す図であり、（ｃ）は反射（散乱）光の位相の時間変化を示す図である。

【図2】本発明で計算されるビート信号の位相の時間変化の一例である。

【図3】本発明で得られるビート信号のパワースペクトル波形の一例である。

【図4】本発明の第１及び第２の実施形態に係る光反射測定装置の構成の一例を示すブロック図である。

【図5】本発明の第１の実施形態に係る光反射測定装置における測定の流れを示すフローチャートである。

【図6】本発明の第２の実施形態に係る光反射測定装置における測定の流れを示すフローチャートである。

【発明を実施するための形態】

【0018】

本発明ではＯＦＤＲ測定系に周回光路を設け、試験光に用いられる周波数掃引連続光を周回光路に通すことで周回時間ごとに複数の試験光を被測定ファイバに入射し、これら複数の試験光による反射（散乱）光の位相を個別解析することで周回時間ごとの位相変化を測定し、得られた位相変化情報に基づいて位相雑音の影響を補償することで、上記課題を解決する。

【0019】

抑々ＯＦＤＲにおいて外乱由来の位相雑音の影響を補償するためには、ビート信号取得時間内における反射（散乱）光の位相の時間変化情報を得ることが必要となるが、従来のＯＦＤＲ測定の試験光は１回のＯＦＤＲ測定につき１つの周波数掃引連続光であるため、１回の測定中に発生する位相変化情報を得ることが不可能であった。

【0020】

一方、本発明では周回光路を通すことにより複数の周波数掃引連続光を作り出すことができ、またそれらは被測定ファイバ到達時刻がそれぞれ異なるため、各周波数掃引連続光による反射（散乱）光信号を個別解析することで反射（散乱）光の位相変化情報を得ることができる。

【0021】

さらに周回光路伝搬中も光源で周波数掃引を続ける（掃引を繰り返すのではなく、１回の掃引を長く続ける）ことで、周回光路伝搬光による信号と同時に従来のＯＦＤＲビート信号（非周回光による信号）を同時に得ることができ、周回光路伝搬光から得られた位相変化情報を用いて非周回光による信号に雑音補償を施すことができる。本発明における周回光路伝搬後の試験光Ｅ’（ｔ）は次式で表される。

【0022】

【数3】

【0023】

【数4】

【0024】

【数5】

【0025】

ここでｎは周回光路における周回数、τ_LoopとΔφは１周回あたりの伝搬遅延時間と位相変化である。なお、τ_Loopが十分小さいことからτ_Loopの２乗項は無視している。このとき、光源の周波数掃引に係る時間幅ＴはＴ＞２τ_Loopを満たし、Ｎ（ｔ）は１以上かつＴ／τ_Loop以下の整数とする。Ｅ’（ｔ）を試験光に用いると、得られるビート信号Ｉ’（ｔ）は次式のようになる。

【0026】

【数6】

【0027】

なお、ここでは遅延時間τ_FUT、τ_Loopの２乗項は無視している。

【0028】

図１に、ビート信号Ｉ’（ｔ）を用いた反射（散乱）光の位相解析の概念図を示す。反射（散乱）光の位相変化情報は、ビート信号Ｉ’（ｔ）をｉτ_Loop＜ｔ＜（ｉ＋１）τ_Loop（ｉは０以上の整数）の範囲で短時間フーリエ変換して得られる複素スペクトルＳ_i（ｆ）を周回光路における周回数毎に計算し、算出した各複素スペクトルＳ_i（ｆ）の位相成分から得る。正の周波数成分（ｆ＞０）のみを考慮すると、Ｉ’（ｔ）の短時間フーリエ変換は次式のようになる。

【0029】

【数7】

【0030】

なお、ここでは周回光路の１周回あたりの伝搬遅延時間τ_Loopが被測定ファイバ２００の位相雑音θ_FUT（ｔ）の変動周期よりも十分短く、ｉτ_Loop＜ｔ＜（ｉ＋１）τ_Loopの時間間隔ではθ_FUT（ｔ）はほぼ一定であり、θ_FUT（ｉτ_Loop）とみなせることとした。

【0031】

式（６）に示されるように、Ｓ₀（ｆ）では１個の反射（散乱）光振幅分布波形がビート周波数領域で現れ、Ｓ_i（ｆ）（ｉが１〜Ｎ−１の整数）ではＮ個の反射（散乱）光振幅分布波形がビート周波数領域で現れる（図１（ｂ））。すなわちｉが１〜Ｎ−１の整数の場合にはＳ_i（ｆ）は周回光路の周回数の異なる光とのビート信号に対応する複数の反射（散乱）光振幅分布波形を含む。Ｓ_i（ｆ）のｆ＝γ（τ_FUT＋ｉτ_Loop）（ｉは０〜Ｎ−１の整数）での値は被測定ファイバ２００の伝搬遅延時間τ_FUTに対応する地点におけるｉτ_Loop＜ｔ＜（ｉ＋１）τ_Loop時点の反射（散乱）光の電界複素振幅を示す。

【0032】

Ｓ_i（ｆ）を異なるｉについて計算すると、各Ｓ_i（ｆ）において最も高周波側の反射（散乱）光振幅分布波形（ｆ＝γ（τ_FUT＋ｉτ_Loop））、すなわち周回光路の周回数が最も多い光とのビート信号に対応する反射（散乱）光振幅分布波形は、それぞれ異なる時刻に発生した同一光周波数の反射（散乱）光振幅分布波形を示す。したがって、Ｓ_i（γ（τ_FUT＋ｉτ_Loop））の位相成分を各ｉの値について求めることで、図１（ｃ）に示すように反射（散乱）光の位相の時間変化を解析することができる。Ｓ_i（γ（τ_FUT＋ｉτ_Loop））の位相成分は次式のようになる。

【0033】

【数8】

【0034】

なお、ここではＳ_i（ｆ）において現れる複数の反射（散乱）光振幅分布波形のビート周波数差（γτ_FUT）がスペクトル分解能（〜１／τ_Loop）よりも十分大きいことから、下記のように近似した。

【0035】

【数9】

【0036】

ｉτ_Loop＝ｔ’と置き換えると、式（７）は次式のように書き換えられる。

【0037】

【数10】

【0038】

したがって、位相雑音θ_FUT（ｔ’）は次式により求められる。

【0039】

【数11】

【0040】

次に、図２に示すように位相雑音θ_FUT（ｔ）が重畳されたＯＦＤＲビート信号の位相の時間変化を算出し、単位位相変化に対応する時間列を算出し、各時間列におけるＩ’（ｔ）を時間軸上で等間隔に並べる（リサンプリングする）ことにより、Ｉ’（ｔ）における位相雑音の影響を補償する。すなわち、次式を満たす時間列ｔ_m（ｍは自然数）でＩ’（ｔ）をリサンプリングする。

【0041】

【数12】

【0042】

ここでωは任意の角周波数、δは任意の位相定数である。

【0043】

【数13】

【0044】

と置き換えると、リサンプリング後のビート信号

【0045】

【数14】

【0046】

は次式のようになる。

【0047】

【数15】

【0048】

【数16】

【0049】

のフーリエ変換の一例を図３に示す。ビート周波数領域において、低周波側から順にｎ＝０，１，２，・・・の反射（散乱）光振幅分布波形が現れる。最も低周波側（ｎ＝０）の波形において、遅延τ_FUT＝ω／２πγに対応する地点で位相雑音θ_FUT（ｔ）項が０となり、位相雑音が補償される。他の地点についても、上記の位相解析からリサンプリングまでの一連の処理を異なるτ_FUTについて実施することで位相雑音補償が可能である。

【0050】

添付の図面を参照して本発明の実施形態を説明する。以下の実施形態は本発明の構成の一例であり、本発明は以下の実施形態に制限されるものではない。

【0051】

（第１の実施形態）
以下に説明する第１の実施形態は、被測定ファイバに振動等の外乱が加わる環境下において、外乱による位相雑音を補償し、高空間分解能光反射測定を行う目的で実施される。

【0052】

図４に、本発明の第１の実施形態に係る光反射測定装置１００の構成を示す。光源には周波数掃引手段を有する周波数掃引光源１０１を用い、時間に対して線形に周波数掃引された連続光が出射される。出射された連続光を光分波器１１１で２分岐し、一方は光合分波器１１２から周回光路１０２に入射し、もう一方は反射（散乱）光をコヒーレント検波する際のローカル光とする。

【0053】

周回光路１０２に入射した連続光は、被測定ファイバ長の２倍以上の長さを有する遅延ファイバ１０４を伝搬後、光増幅器１０３において遅延ファイバ伝搬による光損失が補償される。なお、周回光路１０２伝搬による光損失が生じても反射（散乱）光の位相解析に十分な信号強度が得られる場合は、ここでの光増幅器１０３は必ずしも用いなくてもよい。

【0054】

周回光路１０２伝搬後の連続光は、光合分波器１１２、光サーキュレータ１０８および光位相参照用ダミーファイバ１０９を介して被測定ファイバ２００に入射され、位相参照用ダミーファイバ１０９及び被測定ファイバ２００中で反射、レイリー散乱される。なお、ここでの位相参照用ダミーファイバ１０９は振動等の外乱が加わらない環境に設置する。

【0055】

位相参照用ダミーファイバ１０９及び被測定ファイバ２００からの反射（散乱）光は、光サーキュレータ１０８を介して光合波器１１３に入射されてローカル光と合波され、光合波器１１３からビート信号が出力される。このように光反射測定装置１００において、光分波器１１１、光合波器１１３、およびそれらの間の光経路により光干渉計を構成する。

【0056】

光合波器１１３から出力されるビート信号を受光器１０５で電気信号に変換し、Ａ／Ｄ変換器１０６でデジタル信号に変換する。Ａ／Ｄ変換器１０６では時間幅Ｎτ_Loop（Ｎは自然数、τ_Loopは周回光路１０２の１周伝搬に係る遅延時間）のビート信号を収録する。

【0057】

次に演算処理装置１０７において、被測定ファイバ２００中で生じた外乱による位相雑音を補償する信号処理を施し、被測定ファイバ２００長手方向の反射（散乱）光振幅分布を得る。

【0058】

図５に、本発明の第１の実施形態に係る光反射測定装置１００の演算処理装置１０７で行われる演算処理の流れを示すフローチャートを示す。

【0059】

初めにステップＳ５０１において、図１に示されるようにビート信号をｉτ_Loop＜ｔ＜（ｉ＋１）τ_Loop（ｉ＝０〜Ｎ−１）の範囲で短時間フーリエ変換し、各分割信号をフーリエ変換してＮ個の複素スペクトルＳ_i（ｆ）を得る。

【0060】

次にステップＳ５０２において、Ｓ_i（γ（τ_Ref＋ｉτ_Loop））の位相成分を各ｉについて求めることにより、位相参照用ダミーファイバ１０９の任意地点における散乱光の位相変化Θ_Ref（ｉτ_Loop）を解析する。ここでのτ_Refは、位相参照用ダミーファイバ１０９のファイバ長をＬ_Refとして２Ｌ_Ref／ｖ（ｖはファイバ中の光速）以下の任意の値を用いる。このとき、Θ_Ref（ｉτ_Loop）は次式で与えられる。

【0061】

【数17】

【0062】

なお、ここでは位相参照用ダミーファイバ１０９が外乱の加わらない環境に設置されていることから、位相参照用ダミーファイバ１０９を伝搬中の位相雑音は無視できることする。

【0063】

次にステップＳ５０３において、Ｓ_i（γ（τ_FUT＋ｉτ_Loop））の位相成分を各ｉについて求めることにより、被測定ファイバ２００の任意地点における散乱光の位相変化Θ_FUT（ｉτ_Loop）を解析する。ここでのτ_FUTは、被測定ファイバ長をＬ_FUTとして２Ｌ_Ref／ｖ＜τ_FUT＜２Ｌ_FUT／ｖを満たす任意の値とする。Θ_FUT（ｉτ_Loop）は次式で与えられる。

【0064】

【数18】

【0065】

なお、ステップＳ５０２とステップＳ５０３は実施する順序が逆でもよい。

【0066】

次にステップＳ５０４において、Θ_FUT（ｉτ_Loop）とΘ_Ref（ｉτ_Loop）の位相差から位相雑音を算出する。Θ_FUT（ｉτ_Loop）とΘ_Ref（ｉτ_Loop）の位相差ΔΘ（ｉτ_Loop）は次式のようになる。

【0067】

【数19】

【0068】

ｉτ_Loop＝ｔ’と置き換えると、式（１５）は次式のように書き換えられる。

【0069】

【数20】

【0070】

次にステップＳ５０５において、ステップＳ５０６で実施されるリサンプリング処理に用いる時間列を算出する。時間列ｔ_mは次式により算出する。

【0071】

【数21】

【0072】

次にステップＳ５０６において、時間列ｔ_mに基づいてビート信号をリサンプリングする。リサンプリングの結果、ビート信号は式（１２）に示される形になる。

【0073】

最後にステップＳ５０７において、ステップＳ５０６でリサンプリングしたビート信号をフーリエ変換することにより、反射（散乱）光振幅分布波形を得る。フーリエ変換後に複数現れる振幅分布波形のうち、最も低周波側が位相雑音補償後の波形を示す。

【0074】

上記演算処理により、遅延τ_FUTに対応する地点の反射（散乱）光の位相雑音が補償される。異なる地点について位相雑音補償する場合、異なるτ_FUTについて上記ステップＳ５０３からステップＳ５０７までを実施する。異なるτ_FUTについて得た反射（散乱）光振幅分布波形の各τ_FUTを中心とする任意区間の波形データを切り出し、それらをつなぎ合わせることにより、複数地点で位相雑音補償された反射（散乱）光振幅分布波形を得ることができる。

【0075】

（第２の実施形態）
以下に説明する第２の実施形態は、被測定ファイバ２００に加わった振動の周波数を解析する目的で実施される。

【0076】

本実施形態で用いられる装置構成は第１の実施形態と同様であり、図４における演算処理装置１０７で行われる演算処理が第１の実施形態と異なる。

【0077】

初めに図４に示される装置を用いて被測定ファイバ２００の反射光とローカル光との合波によるビート信号を取得する。ビート信号は第１の実施形態と同様の方法で取得する。

【0078】

次に演算処理装置１０７において、上記ビート信号を用いて被測定ファイバ２００に加わった振動の周波数を解析する。

【0079】

図６に、本発明の第２の実施形態に係る光反射測定装置１００の演算処理装置１０７で行われる演算処理の流れを示すフローチャートを示す。

【0080】

初めにステップＳ６０１からステップＳ６０４の処理により、被測定ファイバの反射（散乱）光の位相の時間変化ΔΘ（ｔ’）を解析する。ステップＳ６０１からステップＳ６０４までは第１の実施形態と同じ処理のため、ここでは説明を省略する。

【0081】

次にステップＳ６０５において、ΔΘ（ｔ’）を用いて被測定ファイバに加わった振動の周波数を解析する。被測定ファイバ２００に振動が加わると、反射（散乱）光は振動周波数で位相変調される。すなわち、ΔΘ（ｔ’）は次式のように記述される。

【0082】

【数22】

【0083】

ここでθ_vは振動による位相変調振幅、ｆ_vは振動周波数である。式（１８）から明らかなように、ΔΘ（ｔ’）をフーリエ変換して周波数解析することにより、被測定ファイバ２００に加わった振動の周波数を解析することができる。

【0084】

尚、第１および第２の実施形態における演算処理装置１０７は、汎用のコンピュータと上記光反射測定方法を実行させるためのプログラムとによっても実現でき、そのプログラムを記録媒体に記録することも、通信ネットワークを通して提供することも可能である。

【符号の説明】

【0085】

１００ 光反射測定装置
１０１ 周波数掃引光源
１０２ 周回光路
１０３ 光増幅器
１０４ 遅延ファイバ
１０５ 受光器
１０６ Ａ／Ｄ変換器
１０７ 演算処理装置
１０８ 光サーキュレータ
１０９ 光位相参照用ダミーファイバ
１１１ 光分波器
１１２ 光合分波器
１１３ 光合波器
２００ 被測定ファイバ

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】