(57)【要約】【課題】ＤＦＢレーザを用いた場合であっても、ファイバレーザや外部共振型レーザを用いたときと同等の空間分解能を維持できる光周波数多重型コヒーレントＯＴＤＲ、試験方法、信号処理装置、及びプログラムを提供することを目的とする。【解決手段】本発明に係るＯＴＤＲは、  光源からの光の光周波数を所定時間毎に所定周波数間隔で変化させて試験光パルスを生成し、前記試験光パルスを順次、被試験光ファイバに入射する光入射手段と、前記光源からの光を局発光として前記被試験光ファイバからの後方散乱光をコヒーレント検波して受信信号を取得する光受信手段と、前記受信信号を前記所定周波数毎に分離し、周波数分離された信号の振幅をそれぞれ自乗して自乗値を生成し、前記自乗値に対してウィナーフィルタ処理を施した後に前記試験光パルスを前記被試験光ファイバに入射した時の遅延時間をそれぞれ補償して加算平均する演算手段と、を備える。【選択図】図１
2019219298.tif 000002
【特許請求の範囲】
【請求項１】
光周波数多重型コヒーレントＯＴＤＲ（Ｏｐｔｉｃａｌ  Ｔｉｍｅ  Ｄｏｍａｉｎ  Ｒｅｆｌｅｃｔｏｍｅｔｒｙ）であって、
光源からの光の光周波数を所定時間毎に所定周波数間隔で変化させて試験光パルスを生成し、前記試験光パルスを順次、被試験光ファイバに入射する光入射手段と、
前記光源からの光を局発光として前記被試験光ファイバからの後方散乱光をコヒーレント検波して受信信号を取得する光受信手段と、
前記受信信号を前記所定周波数毎に分離し、周波数分離された信号の振幅をそれぞれ自乗して自乗値を生成し、前記自乗値に対してウィナーフィルタ処理を施した後に前記試験光パルスを前記被試験光ファイバに入射した時の遅延時間をそれぞれ補償して加算平均する演算手段と、
を備えることを特徴とする光周波数多重型コヒーレントＯＴＤＲ。
【請求項２】
前記ウィナーフィルタ処理は、フーリエ変換した前記自乗値に、フーリエ変換した前記光源の周波数スペクトルの二乗に任意値を加算した値でフーリエ変換した前記光源の周波数スペクトルの複素共役を除した値を乗算し、逆フーリエ変換する処理であることを特徴とする請求項１に記載の光周波数多重型コヒーレントＯＴＤＲ。
【請求項３】
前記光入射手段は、前記光源の光の波長と異なる波長のダミー光を前記試験光パルスに重畳することを特徴とする請求項１又は２に記載の光周波数多重型コヒーレントＯＴＤＲ。
【請求項４】
光周波数多重型コヒーレントＯＴＤＲが行う試験方法であって、
光源からの光の光周波数を所定時間毎に所定周波数間隔で変化させて試験光パルスを生成し、前記試験光パルスを順次、被試験光ファイバに入射する光入射手順と、
前記光源からの光を局発光として前記被試験光ファイバからの後方散乱光をコヒーレント検波して受信信号を取得する光受信手順と、
前記受信信号を前記所定周波数毎に分離し、周波数分離された信号の振幅をそれぞれ自乗して自乗値を生成し、前記自乗値に対してウィナーフィルタ処理を施した後に前記試験光パルスを前記被試験光ファイバに入射した時の遅延時間をそれぞれ補償して加算平均する演算手順と、
を行うことを特徴とする試験方法。
【請求項５】
前記ウィナーフィルタ処理は、フーリエ変換した前記自乗値に、フーリエ変換した前記光源の周波数スペクトルの二乗に任意値を加算した値でフーリエ変換した前記光源の周波数スペクトルの複素共役を除した値を乗算し、逆フーリエ変換する処理であることを特徴とする請求項４に記載の試験方法。
【請求項６】
光源からの光の光周波数を所定時間毎に所定周波数間隔で変化させて試験光パルスを生成し、前記試験光パルスを順次、被試験光ファイバに入射する光入射手段と、
前記光源からの光を局発光として前記被試験光ファイバからの後方散乱光をコヒーレント検波して受信信号を取得する光受信手段と、
を備える光周波数多重型コヒーレントＯＴＤＲの信号処理装置であって、
前記受信信号を前記所定周波数毎に分離し、周波数分離された信号の振幅をそれぞれ自乗して自乗値を生成し、前記自乗値に対してウィナーフィルタ処理を施した後に前記試験光パルスを前記被試験光ファイバに入射した時の遅延時間をそれぞれ補償して加算平均することを特徴とする信号処理装置。
【請求項７】
前記ウィナーフィルタ処理は、フーリエ変換した前記自乗値に、フーリエ変換した前記光源の周波数スペクトルの二乗に任意値を加算した値でフーリエ変換した前記光源の周波数スペクトルの複素共役を除した値を乗算し、逆フーリエ変換する処理であることを特徴とする請求項６に記載の信号処理装置。
【請求項８】
請求項６又は７に記載の信号処理装置としてコンピュータを機能させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本開示は、光線路の光損失分布や断線位置等を測定するための光パルス試験方法に関する。
【背景技術】
【０００２】
従来より、被試験光ファイバ（Ｆｉｂｅｒ  ｕｎｄｅｒ  ｔｅｓｔ、以後ＦＵＴと称する）の距離や光損失分布、断線位置等を測定する技術として、光パルス試験器（Ｏｐｔｉｃａｌ  Ｔｉｍｅ  Ｄｏｍａｉｎ  Ｒｅｆｌｅｃｔｏｍｅｔｒｙ、以後ＯＴＤＲと称する）がある。このＯＴＤＲは、ＦＵＴに試験光パルスを送出し、試験光パルスによってＦＵＴ中で生じる反射光やレイリー後方散乱光（以後、単に後方散乱光と称する）のパワーを時間領域で測定することでＦＵＴの各地点における光の反射率分布（以後、ＯＴＤＲ波形と称する）を測定する装置、方法である。
【０００３】
海底光ケーブル試験用のＯＴＤＲとして特許文献１に光周波数多重型コヒーレントＯＴＤＲ（以後、ＦＤＭ−ＯＴＤＲと称する）が開示されている。ＦＤＭ−ＯＴＤＲは、それぞれ光周波数の異なる複数の光パルスを時間軸に沿って並べた試験光パルス列をＦＵＴに入射し、ＦＵＴ中で生じた複数の光パルス分の後方散乱光を一度に受信処理することで、複数回分の測定を一度に測定できる特徴を持つ。
【０００４】
以下に、特許文献１に記載されるＦＤＭ−ＯＴＤＲについて簡単に説明する。図４は、特許文献１に記載されるＦＤＭ−ＯＴＤＲの構成を説明する図である。当該ＦＤＭ−ＯＴＤＲは、光源１１からの出力光を分岐素子１２で２つに分岐し、この分岐された光の一方を局発光とし、他方を試験光とし、光周波数制御器１３で所定時間間隔毎に試験光の光周波数を所定間隔だけ変化させ、試験光を光パルス化し、被試験光ファイバ１８に繰り返し入射する（図５参照）。そして、当該ＦＤＭ−ＯＴＤＲは、被試験光ファイバ１８の各地点で反射または散乱により発生する後方散乱光と前記局発光を結合素子１９で結合し、その結合光を光受信して電流に変換し、この電流を周波数毎に分離し、試験光の複数の周波数成分による被試験光ファイバ１８からの反射光および後方散乱光の反射率分布を演算処理装置２５で求める。
【０００５】
ここで、光源１１からの出力光について、光周波数制御器１３、演算処理装置２５によって周波数を分離する方法は、次の条件を満足する必要がある。
（条件１）光源１１からの出力光の線幅は、光周波数制御器１３により所定周波数を持続させる時間の逆数よりも小さいこと。
（条件２）光周波数制御器１３による周波数シフト間隔は、周波数を持続させる所定時間間隔の逆数の自然数倍であること。
（条件３）光周波数制御器１３による所定周波数シフト総量は、数値化処理器２４のサンプリングレートの１／２以下であること。
（条件4）演算処理によるフーリエ変換の周波数分解能の自然数倍は、光周波数制御器１３による周波数シフト間隔であること。
【０００６】
通常、光パルス試験器は一測定につき１つのＯＴＤＲ波形しか取得できない（通常光パルス試験器）が、上記条件を満たすように、試験光の周波数をＹ個、所定時間間隔毎に所定周波数間隔だけ変化させ、周波数分離して測定することで、Ｙ個分のＯＴＤＲ波形を同時に取得できる。これらの波形を加算平均処理することで、一回の測定当たりの加算数をＹ倍多くすることができる。このため、当該ＦＤＭ−ＯＴＤＲは、通常光パルス試験器と同じ測定時間で、５ｌｏｇ√Ｙ（ｄＢ）分のＳＷＤＲを向上させることができる（図６参照）。また、このことは、通常光パルス試験器と同じダイナミックレンジを得るのに、当該ＦＤＭ−ＯＴＤＲは測定時間を１／Ｙ短縮できることも意味する。
【先行技術文献】
【特許文献】
【０００７】
【特許文献１】特開２０１１−１６４０７５号公報
【発明の概要】
【発明が解決しようとする課題】
【０００８】
しかしながら、特許文献１記載のＦＤＭ−ＯＴＤＲでは試験光パルスを発生させるレーザ光源のスペクトル線幅が試験光パルス幅の逆数となる信号帯域より十分小さくない場合、デッドゾーンが大きくなり、正しいファイバの損失分布および故障点での反射イベントを示した波形が精度良く得られないといった課題がある。
【０００９】
一般に、レーザ光源の線幅が試験光パルスの帯域に比べ十分小さい条件を満たさない場合、それぞれの光周波数の後方散乱光信号スペクトルサードローブが高くなり、異なる周波数を持った試験光パルス同士のスペクトル成分の重なりが増加し、十分な性能を有した周波数応答フィルタを用いたとしても周波数分離が実行できなくなる。上記の理由により異なる周波数を持った試験光パルスを時間軸に沿って並べたＦＤＭ−ＯＴＤＲは、周波数分離が十分でない場合、異なる時間に入射した試験パルス光の影響により測定で得られる波形が時間的に広がる傾向を示し、反射や中継機の利得といった急峻な反射率変動が起きる地点において大きなデッドゾーンが生じてしまうためである。
【００１０】
上記の理由により、ＦＤＭ−ＯＴＤＲでは、試験光パルス幅の逆数となる信号帯域より十分小さいスペクトル線幅を持つ試験光パルス発生用のレーザ光源を用いる必要がある。例えば、海底光ケーブル測定でＯＴＤＲ波形に要求される空間分解能１ｋｍでは、試験パルス幅は１０μｓとなり、信号帯域は１００ｋＨｚとなる。よって、空間分解能１ｋｍの条件下では、信号帯域１００ｋＨｚより十分小さい１０ｋＨｚ以下の線幅のレーザ光源を用いることが必要となる。
【００１１】
しかしながら、通常、光通信で用いられている分布帰環型（Ｄｉｓｔｒｉｂｕｔｅｄ  Ｆｅｅｄｂａｃｋ：以下、ＤＦＢと称する）レーザは線幅１００ｋＨｚ程度となり、線幅１０ｋＨｚ以下の光学性能は満たさない。
【００１２】
一般にレーザの線幅は、レーザの共振器長に依存するため、共振器長を長く取ることができる構造を持ったレーザ、例えばファイバレーザや外部共振型レーザといったレーザであれば線幅１０ｋＨｚ以下の光学特性を満足する。このようなレーザは、構成部品数も多くＤＦＢレーザに比べれば一般に高価となる。
【００１３】
以上より、ＦＤＭ−ＯＴＤＲには、空間分解能を高めて正しいファイバの損失分布および故障点での反射イベントを示した波形を測定するためには、高価なファイバレーザや外部共振型レーザを試験光パルス生成用の光源として使用しなければならならず、装置価格を低減することが困難という課題がある。
【００１４】
そこで、本発明は上記課題を解決するために、ＤＦＢレーザを用いた場合であっても、ファイバレーザや外部共振型レーザを用いたときと同等の空間分解能を維持できる光周波数多重型コヒーレントＯＴＤＲ、試験方法、信号処理装置、及びプログラムを提供することを目的とする。
【課題を解決するための手段】
【００１５】
上記目的を達成するために、本発明に係るＦＤＭ−ＯＴＤＲでは、試験光の被試験光ファイバからの後方散乱光の反射率分布を求める演算処理において、フーリエ変換後の信号に対してウィナーフィルタの処理を施すこととした。
【００１６】
具体的には、本発明に係る光周波数多重型コヒーレントＯＴＤＲは、
光源からの光の光周波数を所定時間毎に所定周波数間隔で変化させて試験光パルスを生成し、前記試験光パルスを順次、被試験光ファイバに入射する光入射手段と、
前記光源からの光を局発光として前記被試験光ファイバからの後方散乱光をコヒーレント検波して受信信号を取得する光受信手段と、
前記受信信号を前記所定周波数毎に分離し、周波数分離された信号の振幅をそれぞれ自乗して自乗値を生成し、前記自乗値に対してウィナーフィルタ処理を施した後に前記試験光パルスを前記被試験光ファイバに入射した時の遅延時間をそれぞれ補償して加算平均する演算手段と、
を備える。
【００１７】
また、本発明に係る試験方法は、
光源からの光の光周波数を所定時間毎に所定周波数間隔で変化させて試験光パルスを生成し、前記試験光パルスを順次、被試験光ファイバに入射する光入射手順と、
前記光源からの光を局発光として前記被試験光ファイバからの後方散乱光をコヒーレント検波して受信信号を取得する光受信手順と、
前記受信信号を前記所定周波数毎に分離し、周波数分離された信号の振幅をそれぞれ自乗して自乗値を生成し、前記自乗値に対してウィナーフィルタ処理を施した後に前記試験光パルスを前記被試験光ファイバに入射した時の遅延時間をそれぞれ補償して加算平均する演算手順と、
を行う。
【００１８】
さらに、本発明に係る信号処理装置は、
光源からの光の光周波数を所定時間毎に所定周波数間隔で変化させて試験光パルスを生成し、前記試験光パルスを順次、被試験光ファイバに入射する光入射手段と、
前記光源からの光を局発光として前記被試験光ファイバからの後方散乱光をコヒーレント検波して受信信号を取得する光受信手段と、
を備える光周波数多重型コヒーレントＯＴＤＲの信号処理装置であって、
前記受信信号を前記所定周波数毎に分離し、周波数分離された信号の振幅をそれぞれ自乗して自乗値を生成し、前記自乗値に対してウィナーフィルタ処理を施した後に前記試験光パルスを前記被試験光ファイバに入射した時の遅延時間をそれぞれ補償して加算平均することを特徴とする。
【００１９】
ウィナーフィルタを用いることで、ＦＵＴからの戻り光を受信した受信信号に含まれるレーザ光源のスペクトル線幅の影響を補償することができる。このため、安価なＤＦＢレーザを光源に使用してもデッドゾーンの広がりを低減することができる。従って、本発明は、ＤＦＢレーザを用いた場合であっても、ファイバレーザや外部共振型レーザを用いたときと同等の空間分解能を維持できる光周波数多重型コヒーレントＯＴＤＲ、試験方法、及び信号処理装置を提供することができる。
【００２０】
ここで、ウィナーフィルタ処理は、フーリエ変換した前記自乗値に、フーリエ変換した前記光源の周波数スペクトルの二乗に任意値を加算した値でフーリエ変換した前記光源の周波数スペクトルの複素共役を除した値を乗算し、逆フーリエ変換する処理であることが好ましい。
【００２１】
また、前記光入射手段は、前記光源の光の波長と異なる波長のダミー光を前記試験光パルスに重畳することが好ましい。ダミー光を試験光パルスに重畳させて試験光全体の強度変動を抑えることで、その強度を通信用の信号光強度とほぼ同程度に調整し、光サージの影響を抑制することができる。
【００２２】
また、本発明に係るプログラムは、前記信号処理装置としてコンピュータを機能させるためのプログラムである。前記信号処理装置はコンピュータとプログラムによっても実現でき、プログラムを記録媒体に記録することも、ネットワークを通して提供することも可能である。
【発明の効果】
【００２３】
本発明は、ＤＦＢレーザを用いた場合であっても、ファイバレーザや外部共振型レーザを用いたときと同等の空間分解能を維持できる光周波数多重型コヒーレントＯＴＤＲ、試験方法、信号処理装置、及びプログラムを提供することができる。
【図面の簡単な説明】
【００２４】
【図１】本発明に係る光周波数多重型コヒーレントＯＴＤＲの構成を説明する図である。
【図２】本発明に係る光周波数多重型コヒーレントＯＴＤＲの信号処理方法を説明する図である。
【図３】本発明に係る光周波数多重型コヒーレントＯＴＤＲの効果を説明する図である。
【図４】特許文献１に記載される光周波数多重型コヒーレントＯＴＤＲの構成を説明する図である。
【図５】光周波数多重型コヒーレントＯＴＤＲの試験光を説明する図である。
【図６】光周波数多重型コヒーレントＯＴＤＲの効果を説明する図である。
【発明を実施するための形態】
【００２５】
添付の図面を参照して本発明の実施形態を説明する。以下に説明する実施形態は本発明の実施例であり、本発明は、以下の実施形態に制限されるものではない。なお、本明細書及び図面において符号が同じ構成要素は、相互に同一のものを示すものとする。
【００２６】
図１は、本実施形態の光周波数多重型コヒーレントＯＴＤＲ３０１（以下、光パルス試験装置３０１と記載することがある。）の構成を示す図である。光パルス試験装置３０１は、試験光の各周波数成分によるＦＵＴからの反射光および後方散乱光の反射率分布を求めることができる。
【００２７】
光パルス試験装置３０１は、光周波数多重型コヒーレントＯＴＤＲであって、
光源からの光の光周波数を所定時間毎に所定周波数間隔で変化させて試験光パルスを生成し、前記試験光パルスを順次、被試験光ファイバに入射する光入射手段と、
前記光源からの光を局発光として前記被試験光ファイバからの後方散乱光をコヒーレント検波して受信信号を取得する光受信手段と、
前記受信信号を前記所定周波数毎に分離し、周波数分離された信号の振幅をそれぞれ自乗して自乗値を生成し、前記自乗値に対してウィナーフィルタ処理を施した後に前記試験光パルスを前記被試験光ファイバに入射した時の遅延時間をそれぞれ補償して加算平均する演算手段と、
を備える。
【００２８】
［光入射手段］
スペクトル線幅Δνの光波を発生する第一の光源１１からの出力光は合分波器１２で二系統に分岐され、分岐された光の一方を局発光とし、他方は試験光として光周波数制御器１３に入射される。ここで、合分波器１２は具体的には光カプラ等で構成される。
【００２９】
光周波数制御器１３に入射された試験光は、光周波数制御器１３によって所定時間間隔Ｔ毎に所定周波数ｆｋ’（ｋ＝１、２、…、Ｎ、Ｎは周波数多重数）だけ周波数シフトされる。本実施例では、Ｔ＝１０μｓ、Ｎ＝４０、ｆｋ’＝１０８．４＋（ｋ−１）×０．８ＭＨｚとする。ここで光周波数制御器１３は、搬送波や高次変調側波帯を抑圧し、バイアス電圧調整により、＋１次もしくは−１次の変調側波帯のみ出力できる搬送波抑圧光単側波帯変調器（ＳＳＢ−ＳＣ変調器）を用いる。本実施例では、＋１次の変調側波帯が出力されるようバイアス調整を行う。
【００３０】
上記によって周波数制御を受けた試験光は、光パルス化処理器１６に入力され、パルス発生器２８で制御されたタイミングおよびパルス幅で光パルス化される。本実施例では、光パルスの時間波形は矩形波とする。この光パルス化処理器１６は、具体的には音響光学変調器をパルス駆動した音響光学スイッチである。ここで、音響光学スイッチからの出力光は予め音響光学スイッチ製造時に設定されている固定の周波数シフト（ｆＡＯＭとする）を受けるため、各周波数の試験光パルスは局発光に対して
｜ｆｋ’＋ｆＡＯＭ｜＝ｆｋ
の周波数シフト量を持つことになる。
本実施例では、ｆＡＯＭ＝−１００［ＭＨｚ］とするため、ｆｋ＝８．４＋ｋ×０．８［ＭＨｚ］となる。
【００３１】
なお、光周波数制御器１３と光パルス化処理器１６は、信号タイミング制御器１７によって同期された正弦波発生器１４とパルス発生器２８によってそれぞれ駆動されており、光周波数制御器１３で周波数制御を受けた時間の試験光のみが光パルス化して出力されるようにタイミング調整されている。
【００３２】
第二の光源３１は、第一の光源１１とは波長の異なる光源である。例えば、ＦＵＴが海底光増幅中継伝送システムの光ファイバの場合に、第二の光源３１の光をダミー光として試験光パルスに重畳させて試験光全体の強度変動を抑え、その強度を通信用の信号光強度とほぼ同程度に調整することで光サージの影響を抑制することができる。第二の光源３１からのダミー光は、２つの光パルス化処理器（１６、３６）によって消光比を高くした状態で試験光パルスに重畳される。
【００３３】
光パルス化処理器１６より出力される試験光パルスとダミーパルスは光増幅器１５により増幅された後、サーキュレータ１７を通過し、ＦＵＴに入射される。
【００３４】
［光受信手段］
試験光パルスによってＦＵＴ中で生じた後方散乱光は、サーキュレータ１７を通過した後、偏波によるコヒーレント検波効率の変動を抑えるため偏波制御器２９によって測定ごとに偏波状態を変えられた局発光と合分波器１９で合波し、バランス型光受信器２０で受信される。バランス型光受信器２０から出力される後方散乱光と局発光のビート信号は帯域ろ過フィルタ２３によって、不要な高周波成分がカットされた後、数値化処理器２４でサンプリングされる。
【００３５】
［演算手段］
サンプリングされた後の各周波数成分のビート信号は数値演算処理器２５によって周波数分離され、全て加算平均処理のために足し合わされる。最後に、一連の測定および演算処理を繰り返し行い、得られた結果について加算平均処理し、処理された数値列を対数表示し、最終的にＯＴＤＲ波形を得ることができる。
【００３６】
図２（ｉ）は、ＯＴＤＲ波形を得るために、特許文献１の演算処理器２５にて行われる、サンプリングされた信号に対する周波数分離演算処理の方法を説明する図である。ここで、反射光と局発光の偏波面は常に一致しているものと仮定して記述する。まず任意の位置τｒにて反射したフレネル反射光の受信信号ｉ（ｔ）は、以下のように記述される。
【数１】
2019219298.tif 000003
【数２】
2019219298.tif 000004
ここで、ｗｐ（ｔ）は時刻ｔにおける入射試験パルスの強度、θ（ｔ）は時刻ｔにおける位相雑音を表す。式（１）は特許文献１の式（１８）に相当する。
【００３７】
受信信号ｉ（ｔ）を周波数分離処理するために、以下の短時間フーリエ変換を行う。
【数３】
2019219298.tif 000005
ここで、ｗｒは窓関数を表す。式（３）は窓関数ｗｒを時間軸上でτだけ移動させながら信号ｉ（ｔ）との積をとった信号ｗｒ（ｔ−τ）ｉ（ｔ）に対して、フーリエ変換を行っていることを示す。式（３）は特許文献１の式（１９）に相当する。
【００３８】
上記の演算処理を行うことで、受信信号ｉ（ｔ）について、周波数分離処理された中心周波数ｆｋの試験光パルスによる反射光の振幅Ｉｗ（ｆｋ，τ）を得ることができる。
【００３９】
周波数分離された反射光の振幅Ｉｗ（ｆｋ，τ）を自乗し、各周波数の信号毎にパルス入射時の遅延時間（ｋ−１）Ｔをそれぞれ時間シフトさせ、Ｎ波分の周波数信号を加算平均処理する。即ち、式（３−１）よりＦＤＭ−ＯＴＤＲにおけるフレネル反射波形を得られる。
【数３−１】
2019219298.tif 000006
【００４０】
以上の方法が、図２（ｉ）に示した通常の周波数分離の信号処理方法となる。
【００４１】
ＯＴＤＲ波形においても同じ信号処理方法にて波形を表示する。ＯＴＤＲ波形は、フェーディング雑音と偏波の影響を無視すれば、下記の１次元インパルス応答からシミュレートすることができる。
【数３−２】
2019219298.tif 000007
【００４２】
以下、光源１１のレーザ線幅が信号帯域に比べ十分小さくない場合の周波数分離信号処理について説明する。
まず、式（３）は、式（６）のように変形できる。
【数６】
2019219298.tif 000008
ここで、Ｆ［・］はフーリエ変換を示し、Ｗｒ（ｆ）とＩ（ｆ）はそれぞれ
【数７】
2019219298.tif 000009
【数８】
2019219298.tif 000010
を示す。
【００４３】
ここで、受信信号のパワースペクトルＳｄ（ｆ）＝｜Ｉ（ｆ）｜２は、試験光パルスのパワースペクトルＳｐ（ｆ）とレーザ光源の周波数スペクトルＳＬ（ｆ）を用いて
【数９】
2019219298.tif 000011
なる畳み込み演算で表すことができる。ここで、Ｓｐ（ｆ）はパルス形状のみで一意に決まり、レーザ光源の位相雑音特性には依存しない。レーザ光源のＦＭ雑音スペクトルが白色雑音であると仮定すると、レーザ光源の周波数スペクトルＳＬ（ｆ）は以下のローレンツ関数で記述できる。
【数１０】
2019219298.tif 000012
ここで、Δνは光源１１のスペクトル線幅を示している。
【００４４】
これより式（３）の短時間フーリエ変換で求めた受信信号のパワースペクトログラムは式（１３）で記述できる。
【数１３】
2019219298.tif 000013
【００４５】
式（１３）においてＳＬ（ｆ）がレーザ光源の線幅の影響を表す。つまり、式（１３）を解き、
【数１３−１】
2019219298.tif 000014
を求めることでＯＴＤＲ波形においてレーザ光源の線幅の影響を補償することができる。
【００４６】
以下、式（１３）の解き方を示す。簡単のため、
【数１３−２】
2019219298.tif 000015
と書くと、
【数１４】
2019219298.tif 000016
となり、両辺フーリエ変換すると
【数１５】
2019219298.tif 000017
となる。ここで、Ｙ（ｔ）、Ｈ（ｔ）、Ｘ（ｔ）はそれぞれ、
Ｙ（ｔ）＝Ｆ［ｙ（ｆ）］、
Ｈ（ｔ）＝Ｆ［ｈ（ｆ）］、
Ｘ（ｔ）＝Ｆ［ｘ（ｆ）］
を示す。式（１５）を変形し逆フーリエ変換すると、
【数１６】
2019219298.tif 000018
となり、式（１６）が式（１３）の解となる。
【００４７】
この解は、受信信号｜Ｉｗ（ｆ、τ）｜２に雑音が付加されていない条件下では厳密に正しく、式（１３）で表される線形応答に対してＹ（ｔ）／Ｈ（ｔ）は逆フィルタと呼ばれる。｜Ｉｗ（ｆ、τ）｜２に雑音ｎ（ｆ）が付加されている場合は、式（１７）のウィナーフィルタにより解く事が可能となる。
【数１７】
2019219298.tif 000019
ここで＊は複素共役を意味し、Ｎ（ｔ）＝Ｆ［ｎ（ｆ）］で定義する。一般に
【数１８】
2019219298.tif 000020
は未知の関数となるため、定数Γで置き換える。
【００４８】
上述した通り、式（１７）で表されるウィナーフィルタを用いることで式（３）にて示した受信信号｜Ｉｗ（ｆ、τ）｜２に含まれるレーザ光源のスペクトル線幅の影響を補償する事が可能となる。
【００４９】
すなわち、図２（ｉｉ）に示されるように、本実施形態の光パルス試験装置３０１の数値演算処理器２５は、図２（ｉ）の光パルス試験装置の数値演算処理器２５の処理に対して、周波数分離された反射光の振幅Ｉｗ（ｆｋ，τ）の自乗後にウィナーフィルタ処理を施した後にＮ波分の周波数信号を加算平均処理する。
【００５０】
当該ウィナーフィルタ処理は、式（１７）であって、フーリエ変換した前記自乗値（Ｙ（ｔ））に、フーリエ変換した前記光源の周波数スペクトルの二乗（｜Ｈ（ｔ）｜２）に任意値（定数Γ）を加算した値でフーリエ変換した前記光源の周波数スペクトルの複素共役（Ｈ＊（ｔ））を除した値を乗算し、逆フーリエ変換（Ｆ−１［・］）する処理である。
【００５１】
（実施例）
図３は、全長４００ｋｍの光増幅中継線路に対して（ｉ）特許文献１の信号処理方法と（ｉｉ）光パルス試験装置３０１を用いたＯＴＤＲ波形のシミュレーション結果である。
図３（ａ）−１及び（ｂ）−１は光増幅中継線路全体に対するＯＴＤＲ波形、図３（ａ）−２及び（ｂ）−２は距離１５０ｋｍにおける増幅器利得４０ｄＢの反射率変動点近傍の拡大図である。
【００５２】
特許文献１の信号処理方法では、レーザ線幅が４ｋＨｚと信号帯域（Ｂ＝１／Ｔ）１００ｋＨｚに比べ十分小さい場合、デッドゾーンは３．１ｋｍであるのに対し、レーザ線幅が３５ｋＨｚと信号帯域（Ｂ＝１／Ｔ）１００ｋＨｚに比べ十分小さくない場合、デッドゾーンは１５．０ｋｍと広くなる。
【００５３】
一方、光パルス試験装置３０１を用いることでレーザ線幅が３５ｋＨｚの場合においてもデッドゾーンは３．１ｋｍとなり、これはレーザ線幅が４ｋＨｚの結果と一致する。よって、数値演算処理器２５でウィナーフィルタを用いることによって、光源１１のレーザ線幅によるデッドゾーン広がりを補償することができる。
【００５４】
（発明の効果）
本発明によれば、受信後の信号処理においてにウィナーフィルタを用いることで、高価なスペクトル線幅の小さい挟線幅レーザ光源を用いること無く、通常の光通信で用いられる安価なＤＦＢレーザによる光周波数多重型コヒーレントＯＴＤＲを提供することができる。
【符号の説明】
【００５５】
１１：第一の光源
１２：合分波器
１３：光周波数制御器
１４：正弦波発生器
１５：光増幅器
１６：光パルス化処理器
１７：サーキュレータ
１９：合分波器
２０：バランス型光受信器
２１：ミキサー
２２：正弦波発生器
２３：帯域濾過フィルタ
２４：数値化処理器
２５：数値演算処理器
２７：信号タイミング制御器
２８：パルス発生器
２９：偏波制御器
３１：第二の光源
３２：合分波器
３６：光パルス化処理器
【図１】
2019219298.tif 000021
【図２】
2019219298.tif 000022
【図３】
2019219298.tif 000023
【図４】
2019219298.tif 000024
【図５】
2019219298.tif 000025
【図６】
2019219298.tif 000026