レーザ光測定方法及びその測定装置
【課題】 スペクトル線幅が狭い（コヒーレンス時間が長い）レーザ光のコヒーレンス時間を正確に測定する。
【解決手段】 遅延時間τの関数である被測定レーザ光のコヒーレンス関数γ(τ)を測定する場合に、被測定レーザ光（１１）の光周波数を線形に掃引し（１２）、掃引されたレーザ光を分岐し（１３１）、分岐された一方のレーザ光をサーキュレータ（１４）により光ファイバ（１５）に入射してレイリー散乱光を生じさせ、光ファイバ（１５）で生じたレイリー散乱光と分岐されたもう一方のレーザ光とを合波し（１３２）、合波された光を受光して光電流を検出し（１６）、検出された光電流を数値化してフーリエ変換し（１７，１８）、フーリエ変換された光電流値より、遅延０における振幅の標準偏差と遅延τにおける振幅の標準偏差とを算出し、両者の比により遅延τにおけるコヒーレンス関数の絶対値を求める（１９）。
【解決手段】 遅延時間τの関数である被測定レーザ光のコヒーレンス関数γ(τ)を測定する場合に、被測定レーザ光（１１）の光周波数を線形に掃引し（１２）、掃引されたレーザ光を分岐し（１３１）、分岐された一方のレーザ光をサーキュレータ（１４）により光ファイバ（１５）に入射してレイリー散乱光を生じさせ、光ファイバ（１５）で生じたレイリー散乱光と分岐されたもう一方のレーザ光とを合波し（１３２）、合波された光を受光して光電流を検出し（１６）、検出された光電流を数値化してフーリエ変換し（１７，１８）、フーリエ変換された光電流値より、遅延０における振幅の標準偏差と遅延τにおける振幅の標準偏差とを算出し、両者の比により遅延τにおけるコヒーレンス関数の絶対値を求める（１９）。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、レーザ光のコヒーレンス時間を正確に測定するレーザ光測定方法及びその測定装置に関する。
【背景技術】
【０００２】
従来、比較的スペクトル線幅が狭い（コヒーレンス時間が長い）レーザ光のコヒーレンス時間を正確に測定するために、レーザ光のコヒーレンス時間を関数として測定する方法が提案されている。
このレーザ光のコヒーレンス関数の測定に関わる方法としては、非特許文献１に記載のレーザ光のスペクトル測定方法が知られている。この測定方法は、被測定レーザ光を２分岐し、一方を十分な長さの光ファイバに入射し、もう一方を音響光学素子などによりその光周波数を一定量だけシフトする。その後、両者を合波し、受光素子によりビートスペクトルを測定する。
【０００３】
上記ビートスペクトル幅の１／√２が、レーザ光の発振スペクトル幅であることが知られていることから、上記ビートスペクトル幅からレーザ光のスペクトル線幅を測定することができる。また、このレーザ光のスペクトル線幅をΔνとすると、レーザ光のコヒーレンス時間は１／Δνであるので、これよりレーザ光のコヒーレンス時間を求めることができる。また、これに光速を乗算すればコヒーレンス長が得られる。
【０００４】
但し、上記の測定方法には以下のような問題点が知られている。すなわち、上記説明の中で、分岐されたレーザ光の一方は十分に長い光ファイバにより、十分な時間遅延を与えてからもう一方と合波する必要があり、被測定レーザ光のコヒーレンス長よりも十分に長い遅延が要求される。即ち、従来の測定方法は、まさに被測定量であり、よって現段階では不明のレーザ光のコヒーレンス長よりも長い遅延が得られていることを前提に有効となる。したがって、従来の測定方法による測定結果は、レーザ光スペクトル幅の一定の推測を可能にするものの、不確実性を含んだものにならざるを得なかった。
【０００５】
上記の懸念は、特にスペクトル線幅が狭い（コヒーレンス時間が長い）レーザ光を測定する場合に顕著であり、近年においてファイバレーザなどの数ｋＨｚ程度のスペクトル線幅を有するレーザが商用化されている状況では、その測定のあいまいさが問題となっていた。
【先行技術文献】
【非特許文献】
【０００６】
【非特許文献１】T. Okoshi, K. Kikuchi, and A. Nakayama, “Novel method for high resolution measurement of laser output spectrum,” Electronics Letters, Vol. 16, No. 16, pp. 630-631, 1980
【非特許文献２】Y. Koshikiya, X. Fan, and F. Ito, “Long range and cm-level spatial resolution measurement using coherent optical frequency domain reflectmetry with SSB-SC modulator and narrow linewidth fiber laser”, IEEE/OSA J. Lightwave Technol. Vol. 26, No. 18, pp. 3287-3294, 2008
【発明の概要】
【発明が解決しようとする課題】
【０００７】
しかしながら、上記のような従来の測定法では、分岐されたレーザ光の一方を十分に長い光ファイバを伝送させ、被測定レーザ光のコヒーレンス長よりも十分に長い遅延を生じさせた後、もう一方と合波し測定を行うが、ここでの測定量は、遅延を生じさせるために用意された光ファイバによって生じるコヒーレンス長であり、測定時には不明であるレーザ光のコヒーレンス長よりも長い遅延が得られていることが前提条件となっている。したがって、測定結果は、レーザ光スペクトル幅の推測は可能であるものの、その測定結果は不確実性を含んでいるという問題点があった。
【０００８】
本発明は、上記の事情に着目してなされたもので、その目的とするところは、スペクトル線幅が狭い（コヒーレンス時間が長い）レーザ光のコヒーレンス時間を正確に測定することのできるレーザ光測定方法及びその測定装置を提供することにある。
【課題を解決するための手段】
【０００９】
本発明に係るレーザ光測定方法は以下のような態様の構成とする。
（１）遅延時間τの関数である被測定レーザ光のコヒーレンス関数γ(τ)を測定する方法であって、前記被測定レーザ光の光周波数を線形に掃引し、前記光周波数が掃引されたレーザ光を分岐し、前記分岐された一方のレーザ光を光ファイバに入射してレイリー散乱光を生じさせ、前記光ファイバで生じたレイリー散乱光と前記分岐されたもう一方のレーザ光とを合波し、前記合波された光を受光して光電流を検出し、前記検出された光電流を数値化してフーリエ変換し、フーリエ変換された光電流値より、遅延０における振幅の標準偏差と遅延τにおける振幅の標準偏差とを算出し、両者の比により遅延τにおけるコヒーレンス関数の絶対値を求める態様とする。
【００１０】
（２）（１）の構成において、さらに、前記コヒーレンス関数の絶対値が1/eになるτをもってコヒーレンス時間を求める態様とする。
本発明に係るレーザ光測定装置は以下のような態様の構成とする。
（３）遅延時間τの関数であるレーザ光のコヒーレンス関数γ(τ)を測定する装置であって、被測定レーザ光の光周波数を線形に掃引する光周波数掃引手段と、前記光周波数掃引がされたレーザ光を分岐する分岐手段と、前記分岐された一方のレーザ光を入射してレイリー散乱光を生じさせる光ファイバと、前記光ファイバで生じたレイリー散乱光と前記分岐されたもう一方のレーザ光とを合波する合波手段と、前記合波された光を受光して光電流を検出する検出手段と、前記検出された光電流を数値化してフーリエ変換し、フーリエ変換された光電流値より、遅延０における振幅の標準偏差と遅延τにおける振幅の標準偏差とを算出し、両者の比により遅延τにおけるコヒーレンス関数の絶対値を求める演算手段とを具備する態様とする。
【００１１】
（４）（３）の構成において、さらに、前記コヒーレンス関数の絶対値が1/eになるτをもってコヒーレンス時間を求める態様とする。
【発明の効果】
【００１２】
以上のように、本発明は、測定値に対し統計処理を行うことにより、従来よりも信頼性を向上させることを特徴としている。具体的には、周波数掃引された被測定レーザ光をC-OFDRを応用しレイリー散乱強度を測定し得られたレイリー散乱強度の絶対値の標準偏差を算出し、τ（コヒーレンス時間）におけるコヒーレンス関数の絶対値を求めるものである。
【００１３】
本発明によれば、用意した光ファイバの全長の光の往復時間よりもτが小さい場合であっても、コヒーレンス関数の絶対値を求めることが可能となる。さらに、その値が1/eに減ずるτを求めることにより、レーザ光のコヒーレンス時間を求めることも可能となる。
【００１４】
要するに本発明の構成によれば、スペクトル線幅が狭い（コヒーレンス時間が長い）レーザ光のコヒーレンス時間を正確に測定することのできるレーザ光測定方法及びその測定装置を提供することができる。
【図面の簡単な説明】
【００１５】
【図１】本発明に係るレーザ光測定装置の一実施形態の構成を示すブロック図。
【図２】図１に示す実施形態の周波数掃引によって与えられる光周波数の変化を示す特性図。
【発明を実施するための形態】
【００１６】
以下、図面を参照して本発明の実施形態を説明する。
図１は本発明に係るレーザ光測定装置の一実施形態の構成を示すブロック図である。図１において、１１は被測定レーザ光源、１２は被測定レーザ光源１１で発生されるレーザ光を入射して光周波数を掃引する機能を有する光周波数掃引装置である。この光周波数掃引装置１２で光周波数が掃引されたレーザ光は、光伝送路に送出され、光分岐器１３１により２分岐される。
【００１７】
ここで分岐された一方のレーザ光は光サーキュレータ１４を介して光ファイバ１５に入射される。この光ファイバ１５内ではレイリー散乱と呼ばれる光散乱が生じる。この散乱光は光ファイバ１５を逆行し、光サーキュレータ１４によって光合波器１３２に向けて送出される。光分岐器１３１で分岐された他方のレーザ光はそのまま光合波器１３２に送られ、光サーキュレータ１４からのレーザ光と合波された後、再度２系統に分配されてバランス型受光素子１４の受光面に照射される。
【００１８】
上記バランス型受光素子１４では、２つの光合波出力それぞれを個別に受光し、それぞれの受光量に応じた光電流を発生する。このバランス型受光素子１４で得られた２つの光電流出力はデータ取得部１７に送られる。
このデータ取得部１７は、例えばＡ／Ｄコンバータ等による数値化処理部と数値化されたデータを格納するメモリによって構成され、バランス型受光素子１６に生じた光電流を数値化し、メモリに記録しておく。
【００１９】
一方、フーリエ変換部１８及びコヒーレンス関数解析部１９は、パーソナルコンピュータ等の数値演算部で実現されるものであり、上記光電流をフーリエ変換してコヒーレンス関数解析を行う。
上記構成において、以下に処理動作を説明する。
まず、被測定レーザ光源１１が発するレーザ光の電界振幅Ｅ(t)を（１）式のように表すこととする。
【００２０】
【数１】
ここでθ(t)は、レーザ光の位相雑音を表す確率変数である。
レーザ光のコヒーレンス関数γ(τ)は、（２）式で与えられる。
【００２１】
【数２】
ここで、かっこ記号の＜＞は統計平均を意味するが、通常のレーザ光においてはこれを時間平均で置き換えることができる。その場合、コヒーレンス関数は（３）式により定義される。
【００２２】
【数３】
【００２３】
【数４】
【００２４】
【数５】
で表される。また、＊は位相共役を意味する。
被測定レーザ光源１１が発した被測定レーザ光は、光周波数掃引装置１２によって、図２に示すようにその光周波数が時間に対してＴ秒間線形に掃引される。光周波数掃引装置１２としては、例えば非特許文献２に記載の単一側波帯変調器を用いることにより実現できる。
【００２５】
光周波数掃引装置１２により光周波数が掃引された光波の電界振幅ＥSW(t)は、（４）式で表される。
【００２６】
【数６】
ここでｇは光周波数の掃引速度（Hz/s）である。
周波数掃引されたレーザ光は、光分岐器１３１により２分岐され、一方は光サーキュレータ１４を介して光ファイバ１５に入射される。光ファイバ１５内ではレイリー散乱と呼ばれる光散乱が生じ、その散乱光は光ファイバ１５を逆方向に伝搬して光サーキュレータ１４に戻り、このサーキュレータ１４の機能によって光合波器１５に送り出される。分岐されたもう一方のレーザ光は、そのまま光合波器１３２に向かって進行する。光合波器１３２ではこれらのレーザ光が合波される。
【００２７】
本構成は、同じく非特許文献２に説明されているコヒーレント光周波数領域反射計(Coherent Optical Frequency Domain Reflectometry：C-OFDR)と呼ばれる反射分布測定装置の構成を応用したものであり、光ファイバ１５内の１つの反射点Ｚmまでの光の往復時間をτm とすると、この反射点からの散乱光によってバランス受光素子１６に生じる光電流ｉm(t)は、（５）式により表される。
【００２８】
【数７】
ここで複素数ｒm は反射点ｚm のレイリー反射係数であり、ランダムな値をとる確率変数である。
（５）式のフーリエ変換を（６）式に変換すると、
【００２９】
【数８】
（６）式に（５）式を代入してｆ＝ｇτm と置くことにより（７）式を得る。
【００３０】
【数９】
ここでＴが十分に大きく、
【００３１】
【数１０】
が成立すると、（７）式は（８）式に示すようになる。
【００３２】
【数１１】
一般に、光ファイバにコヒーレントなレーザ光を入射すると、そのレイリー散乱光はランダムな強度揺らぎを持つことが知られているが、上記の装置構成においてもそれは同様である。このことから、（８）式は光ファイバの微小体積からのコヒーレントな散乱量を表し、コヒーレントな散乱の全散乱強度Ω(gτm )はそれらの無数の和で決まる。すなわち（９）式に示すようになる。
【００３３】
【数１２】
尚、バランス型受光素子１６に生じる全電流（あらゆる反射点からの電流の和）をｉ(t)とすれば、
【００３４】
【数１３】
であるので、Ω(gτm )は全電流ｉ(t)のフーリエ変換によって以下のように求められることは明らかである。
【００３５】
【数１４】
（８）式によれば、遅延時間がτm の場合には、微小体積からのコヒーレントな散乱の大きさがγ(τm )倍に小さくなる（コヒーレンス関数γ(τ)は複素数であるが、その絶対値は０以上１以下である）ので、その近傍でのレイリー散乱強度の揺らぎの大きさもγ(τm )倍に小さくなることになる。
【００３６】
これを利用して、上記観測からコヒーレンス関数を算出するためには、観測可能量であるレイリー散乱強度の絶対値｜Ω(gτm )｜の統計量を計算する。まずτ＝０においては、γ(τ)＝γ(0)＝１であり、τm 近傍ではγ(τm )はほぼ一定と考えられるので、散乱強度Ω(0)の絶対値の標準偏差は、
【００３７】
【数１５】
の標準偏差σ|γ|に等しい。すなわち（１０）式に示すようになる。
【００３８】
【数１６】
（９）式において、ｒm とγ(τm )はいずれも複素数であり、かつｒm はランダムな確率複素変数であって、その偏角は０〜２πまで均一に分布する。よって、γ(τm )の偏角をｒm に組み入れて新しい変数ｒ’m を定義し、γ(τm )はその絶対値のみを残す形に表記しても問題ない。すなわち
【００３９】
【数１７】
が成立する。ｒ’m はｒm と全く同じ性質の確率変数である。遅延τm におけるΩ(gτm )の絶対値の標準偏差は、（１１）式より直ちに
【００４０】
【数１８】
と求めることができる。（１０）式、（１２）式の比をとってσ|γ|を消去すれば、以下のようにコヒーレンス関数を求めることができる。
【００４１】
【数１９】
標準偏差を求めるには多数の標本を必要とする。C-OFDRの理論によれば、観測されるレイリー散乱強度は、光周波数１／Ｔごとに独立である。したがって、例えばΩ(0)の標準偏差を求めるに当たっては、τ=０近傍において、Δｆ＝１／Ｔごとに標本
【００４２】
【数２０】
を収集し、これらの標準偏差を求めればよい。Ω(gτm )の標準偏差を求める場合も同様である。
以上の知見にもとづき、図１の構成において次の手順でコヒーレンス関数を求める。
データ取得部１７は、例えばＡ／Ｄコンバータによる数値化装置と数値化された値を格納するメモリによって構成され、バランス型受光素子１６に生じた光電流ｉ(t)を数値化し、メモリに記録する。フーリエ変換部１８は、パーソナルコンピュータ等の数値演算処理装置で構成され、上記光電流ｉ(t)をフーリエ変換し、
【００４３】
【数２１】
を得る。
最後に、コヒーレンス関数解析部１９は、Ω(0)近傍及びΩ(gτ)近傍でのそれぞれの標準偏差を算出し、
【００４４】
【数２２】
によりτにおけるコヒーレンス関数の絶対値を求める。
この作業は、使用する光ファイバ１５の全長の光の往復時間よりも小さいτの場合であっても可能であり、その範囲でのコヒーレンス関数の絶対値を求めることができる。その値が1/eに減ずるτを求めれば、光源のコヒーレンス時間を求めることもできる。
【００４５】
但し、標準偏差の代わりに変動係数（標準偏差／標準偏差の算出に使用したポイント数）を用いてもコヒーレンス関数を求めることができる。
上記実施形態によるコヒーレンス関数の測定方法並びにコヒーレンス時間の測定方法は、従来技術に対して以下の優位性を持つ。
【００４６】
第１に、コヒーレンス時間と同程度の光ファイバ長により測定が可能な上、仮にコヒーレンス時間が光ファイバ長よりもはるかに長かった場合には、そのこと自体が測定により明らかにされるので、より長い光ファイバに置き換えて測定をやり直すなどの処置を取ることができる。よって、本実施形態は、従来技術よりも信頼性に優れた測定結果を提供するといえる。
【００４７】
第２に、単にコヒーレンス時間の測定が可能となるだけでなく、そのコヒーレンス関数の絶対値の全容が掌握される。これにより、レーザ光の性質をより詳しく測定できると考えられる。
以上、実施形態に記載の構成によれば、スペクトル線幅が狭い（コヒーレンス時間が長い）レーザ光のコヒーレンス時間を正確に測定することのできるレーザ光測定方法及びその測定装置が実現する。
【００４８】
尚、本発明は上記実施形態そのままに限定されるものではなく、実施段階ではその要旨を逸脱しない範囲で構成要素を変形して具体化できる。また、上記実施形態に開示されている複数の構成要素の適宜な組み合わせにより、種々の発明を形成できる。例えば、実施形態に示される全構成要素から幾つかの構成を削除してもよい。さらに、異なる実施形態例に亘る構成要素を適宜組み合わせてもよい。
【符号の説明】
【００４９】
１１…被測定レーザ光源、１２…光周波数掃引装置、１３１…光分岐器、１３２…光合波器、１４…光サーキュレータ、１５…光ファイバ、１４…バランス型受光素子、１７…データ取得部、１８…フーリエ変換部、１９…コヒーレンス関数解析部。
【技術分野】
【０００１】
本発明は、レーザ光のコヒーレンス時間を正確に測定するレーザ光測定方法及びその測定装置に関する。
【背景技術】
【０００２】
従来、比較的スペクトル線幅が狭い（コヒーレンス時間が長い）レーザ光のコヒーレンス時間を正確に測定するために、レーザ光のコヒーレンス時間を関数として測定する方法が提案されている。
このレーザ光のコヒーレンス関数の測定に関わる方法としては、非特許文献１に記載のレーザ光のスペクトル測定方法が知られている。この測定方法は、被測定レーザ光を２分岐し、一方を十分な長さの光ファイバに入射し、もう一方を音響光学素子などによりその光周波数を一定量だけシフトする。その後、両者を合波し、受光素子によりビートスペクトルを測定する。
【０００３】
上記ビートスペクトル幅の１／√２が、レーザ光の発振スペクトル幅であることが知られていることから、上記ビートスペクトル幅からレーザ光のスペクトル線幅を測定することができる。また、このレーザ光のスペクトル線幅をΔνとすると、レーザ光のコヒーレンス時間は１／Δνであるので、これよりレーザ光のコヒーレンス時間を求めることができる。また、これに光速を乗算すればコヒーレンス長が得られる。
【０００４】
但し、上記の測定方法には以下のような問題点が知られている。すなわち、上記説明の中で、分岐されたレーザ光の一方は十分に長い光ファイバにより、十分な時間遅延を与えてからもう一方と合波する必要があり、被測定レーザ光のコヒーレンス長よりも十分に長い遅延が要求される。即ち、従来の測定方法は、まさに被測定量であり、よって現段階では不明のレーザ光のコヒーレンス長よりも長い遅延が得られていることを前提に有効となる。したがって、従来の測定方法による測定結果は、レーザ光スペクトル幅の一定の推測を可能にするものの、不確実性を含んだものにならざるを得なかった。
【０００５】
上記の懸念は、特にスペクトル線幅が狭い（コヒーレンス時間が長い）レーザ光を測定する場合に顕著であり、近年においてファイバレーザなどの数ｋＨｚ程度のスペクトル線幅を有するレーザが商用化されている状況では、その測定のあいまいさが問題となっていた。
【先行技術文献】
【非特許文献】
【０００６】
【非特許文献１】T. Okoshi, K. Kikuchi, and A. Nakayama, “Novel method for high resolution measurement of laser output spectrum,” Electronics Letters, Vol. 16, No. 16, pp. 630-631, 1980
【非特許文献２】Y. Koshikiya, X. Fan, and F. Ito, “Long range and cm-level spatial resolution measurement using coherent optical frequency domain reflectmetry with SSB-SC modulator and narrow linewidth fiber laser”, IEEE/OSA J. Lightwave Technol. Vol. 26, No. 18, pp. 3287-3294, 2008
【発明の概要】
【発明が解決しようとする課題】
【０００７】
しかしながら、上記のような従来の測定法では、分岐されたレーザ光の一方を十分に長い光ファイバを伝送させ、被測定レーザ光のコヒーレンス長よりも十分に長い遅延を生じさせた後、もう一方と合波し測定を行うが、ここでの測定量は、遅延を生じさせるために用意された光ファイバによって生じるコヒーレンス長であり、測定時には不明であるレーザ光のコヒーレンス長よりも長い遅延が得られていることが前提条件となっている。したがって、測定結果は、レーザ光スペクトル幅の推測は可能であるものの、その測定結果は不確実性を含んでいるという問題点があった。
【０００８】
本発明は、上記の事情に着目してなされたもので、その目的とするところは、スペクトル線幅が狭い（コヒーレンス時間が長い）レーザ光のコヒーレンス時間を正確に測定することのできるレーザ光測定方法及びその測定装置を提供することにある。
【課題を解決するための手段】
【０００９】
本発明に係るレーザ光測定方法は以下のような態様の構成とする。
（１）遅延時間τの関数である被測定レーザ光のコヒーレンス関数γ(τ)を測定する方法であって、前記被測定レーザ光の光周波数を線形に掃引し、前記光周波数が掃引されたレーザ光を分岐し、前記分岐された一方のレーザ光を光ファイバに入射してレイリー散乱光を生じさせ、前記光ファイバで生じたレイリー散乱光と前記分岐されたもう一方のレーザ光とを合波し、前記合波された光を受光して光電流を検出し、前記検出された光電流を数値化してフーリエ変換し、フーリエ変換された光電流値より、遅延０における振幅の標準偏差と遅延τにおける振幅の標準偏差とを算出し、両者の比により遅延τにおけるコヒーレンス関数の絶対値を求める態様とする。
【００１０】
（２）（１）の構成において、さらに、前記コヒーレンス関数の絶対値が1/eになるτをもってコヒーレンス時間を求める態様とする。
本発明に係るレーザ光測定装置は以下のような態様の構成とする。
（３）遅延時間τの関数であるレーザ光のコヒーレンス関数γ(τ)を測定する装置であって、被測定レーザ光の光周波数を線形に掃引する光周波数掃引手段と、前記光周波数掃引がされたレーザ光を分岐する分岐手段と、前記分岐された一方のレーザ光を入射してレイリー散乱光を生じさせる光ファイバと、前記光ファイバで生じたレイリー散乱光と前記分岐されたもう一方のレーザ光とを合波する合波手段と、前記合波された光を受光して光電流を検出する検出手段と、前記検出された光電流を数値化してフーリエ変換し、フーリエ変換された光電流値より、遅延０における振幅の標準偏差と遅延τにおける振幅の標準偏差とを算出し、両者の比により遅延τにおけるコヒーレンス関数の絶対値を求める演算手段とを具備する態様とする。
【００１１】
（４）（３）の構成において、さらに、前記コヒーレンス関数の絶対値が1/eになるτをもってコヒーレンス時間を求める態様とする。
【発明の効果】
【００１２】
以上のように、本発明は、測定値に対し統計処理を行うことにより、従来よりも信頼性を向上させることを特徴としている。具体的には、周波数掃引された被測定レーザ光をC-OFDRを応用しレイリー散乱強度を測定し得られたレイリー散乱強度の絶対値の標準偏差を算出し、τ（コヒーレンス時間）におけるコヒーレンス関数の絶対値を求めるものである。
【００１３】
本発明によれば、用意した光ファイバの全長の光の往復時間よりもτが小さい場合であっても、コヒーレンス関数の絶対値を求めることが可能となる。さらに、その値が1/eに減ずるτを求めることにより、レーザ光のコヒーレンス時間を求めることも可能となる。
【００１４】
要するに本発明の構成によれば、スペクトル線幅が狭い（コヒーレンス時間が長い）レーザ光のコヒーレンス時間を正確に測定することのできるレーザ光測定方法及びその測定装置を提供することができる。
【図面の簡単な説明】
【００１５】
【図１】本発明に係るレーザ光測定装置の一実施形態の構成を示すブロック図。
【図２】図１に示す実施形態の周波数掃引によって与えられる光周波数の変化を示す特性図。
【発明を実施するための形態】
【００１６】
以下、図面を参照して本発明の実施形態を説明する。
図１は本発明に係るレーザ光測定装置の一実施形態の構成を示すブロック図である。図１において、１１は被測定レーザ光源、１２は被測定レーザ光源１１で発生されるレーザ光を入射して光周波数を掃引する機能を有する光周波数掃引装置である。この光周波数掃引装置１２で光周波数が掃引されたレーザ光は、光伝送路に送出され、光分岐器１３１により２分岐される。
【００１７】
ここで分岐された一方のレーザ光は光サーキュレータ１４を介して光ファイバ１５に入射される。この光ファイバ１５内ではレイリー散乱と呼ばれる光散乱が生じる。この散乱光は光ファイバ１５を逆行し、光サーキュレータ１４によって光合波器１３２に向けて送出される。光分岐器１３１で分岐された他方のレーザ光はそのまま光合波器１３２に送られ、光サーキュレータ１４からのレーザ光と合波された後、再度２系統に分配されてバランス型受光素子１４の受光面に照射される。
【００１８】
上記バランス型受光素子１４では、２つの光合波出力それぞれを個別に受光し、それぞれの受光量に応じた光電流を発生する。このバランス型受光素子１４で得られた２つの光電流出力はデータ取得部１７に送られる。
このデータ取得部１７は、例えばＡ／Ｄコンバータ等による数値化処理部と数値化されたデータを格納するメモリによって構成され、バランス型受光素子１６に生じた光電流を数値化し、メモリに記録しておく。
【００１９】
一方、フーリエ変換部１８及びコヒーレンス関数解析部１９は、パーソナルコンピュータ等の数値演算部で実現されるものであり、上記光電流をフーリエ変換してコヒーレンス関数解析を行う。
上記構成において、以下に処理動作を説明する。
まず、被測定レーザ光源１１が発するレーザ光の電界振幅Ｅ(t)を（１）式のように表すこととする。
【００２０】
【数１】
ここでθ(t)は、レーザ光の位相雑音を表す確率変数である。
レーザ光のコヒーレンス関数γ(τ)は、（２）式で与えられる。
【００２１】
【数２】
ここで、かっこ記号の＜＞は統計平均を意味するが、通常のレーザ光においてはこれを時間平均で置き換えることができる。その場合、コヒーレンス関数は（３）式により定義される。
【００２２】
【数３】
【００２３】
【数４】
【００２４】
【数５】
で表される。また、＊は位相共役を意味する。
被測定レーザ光源１１が発した被測定レーザ光は、光周波数掃引装置１２によって、図２に示すようにその光周波数が時間に対してＴ秒間線形に掃引される。光周波数掃引装置１２としては、例えば非特許文献２に記載の単一側波帯変調器を用いることにより実現できる。
【００２５】
光周波数掃引装置１２により光周波数が掃引された光波の電界振幅ＥSW(t)は、（４）式で表される。
【００２６】
【数６】
ここでｇは光周波数の掃引速度（Hz/s）である。
周波数掃引されたレーザ光は、光分岐器１３１により２分岐され、一方は光サーキュレータ１４を介して光ファイバ１５に入射される。光ファイバ１５内ではレイリー散乱と呼ばれる光散乱が生じ、その散乱光は光ファイバ１５を逆方向に伝搬して光サーキュレータ１４に戻り、このサーキュレータ１４の機能によって光合波器１５に送り出される。分岐されたもう一方のレーザ光は、そのまま光合波器１３２に向かって進行する。光合波器１３２ではこれらのレーザ光が合波される。
【００２７】
本構成は、同じく非特許文献２に説明されているコヒーレント光周波数領域反射計(Coherent Optical Frequency Domain Reflectometry：C-OFDR)と呼ばれる反射分布測定装置の構成を応用したものであり、光ファイバ１５内の１つの反射点Ｚmまでの光の往復時間をτm とすると、この反射点からの散乱光によってバランス受光素子１６に生じる光電流ｉm(t)は、（５）式により表される。
【００２８】
【数７】
ここで複素数ｒm は反射点ｚm のレイリー反射係数であり、ランダムな値をとる確率変数である。
（５）式のフーリエ変換を（６）式に変換すると、
【００２９】
【数８】
（６）式に（５）式を代入してｆ＝ｇτm と置くことにより（７）式を得る。
【００３０】
【数９】
ここでＴが十分に大きく、
【００３１】
【数１０】
が成立すると、（７）式は（８）式に示すようになる。
【００３２】
【数１１】
一般に、光ファイバにコヒーレントなレーザ光を入射すると、そのレイリー散乱光はランダムな強度揺らぎを持つことが知られているが、上記の装置構成においてもそれは同様である。このことから、（８）式は光ファイバの微小体積からのコヒーレントな散乱量を表し、コヒーレントな散乱の全散乱強度Ω(gτm )はそれらの無数の和で決まる。すなわち（９）式に示すようになる。
【００３３】
【数１２】
尚、バランス型受光素子１６に生じる全電流（あらゆる反射点からの電流の和）をｉ(t)とすれば、
【００３４】
【数１３】
であるので、Ω(gτm )は全電流ｉ(t)のフーリエ変換によって以下のように求められることは明らかである。
【００３５】
【数１４】
（８）式によれば、遅延時間がτm の場合には、微小体積からのコヒーレントな散乱の大きさがγ(τm )倍に小さくなる（コヒーレンス関数γ(τ)は複素数であるが、その絶対値は０以上１以下である）ので、その近傍でのレイリー散乱強度の揺らぎの大きさもγ(τm )倍に小さくなることになる。
【００３６】
これを利用して、上記観測からコヒーレンス関数を算出するためには、観測可能量であるレイリー散乱強度の絶対値｜Ω(gτm )｜の統計量を計算する。まずτ＝０においては、γ(τ)＝γ(0)＝１であり、τm 近傍ではγ(τm )はほぼ一定と考えられるので、散乱強度Ω(0)の絶対値の標準偏差は、
【００３７】
【数１５】
の標準偏差σ|γ|に等しい。すなわち（１０）式に示すようになる。
【００３８】
【数１６】
（９）式において、ｒm とγ(τm )はいずれも複素数であり、かつｒm はランダムな確率複素変数であって、その偏角は０〜２πまで均一に分布する。よって、γ(τm )の偏角をｒm に組み入れて新しい変数ｒ’m を定義し、γ(τm )はその絶対値のみを残す形に表記しても問題ない。すなわち
【００３９】
【数１７】
が成立する。ｒ’m はｒm と全く同じ性質の確率変数である。遅延τm におけるΩ(gτm )の絶対値の標準偏差は、（１１）式より直ちに
【００４０】
【数１８】
と求めることができる。（１０）式、（１２）式の比をとってσ|γ|を消去すれば、以下のようにコヒーレンス関数を求めることができる。
【００４１】
【数１９】
標準偏差を求めるには多数の標本を必要とする。C-OFDRの理論によれば、観測されるレイリー散乱強度は、光周波数１／Ｔごとに独立である。したがって、例えばΩ(0)の標準偏差を求めるに当たっては、τ=０近傍において、Δｆ＝１／Ｔごとに標本
【００４２】
【数２０】
を収集し、これらの標準偏差を求めればよい。Ω(gτm )の標準偏差を求める場合も同様である。
以上の知見にもとづき、図１の構成において次の手順でコヒーレンス関数を求める。
データ取得部１７は、例えばＡ／Ｄコンバータによる数値化装置と数値化された値を格納するメモリによって構成され、バランス型受光素子１６に生じた光電流ｉ(t)を数値化し、メモリに記録する。フーリエ変換部１８は、パーソナルコンピュータ等の数値演算処理装置で構成され、上記光電流ｉ(t)をフーリエ変換し、
【００４３】
【数２１】
を得る。
最後に、コヒーレンス関数解析部１９は、Ω(0)近傍及びΩ(gτ)近傍でのそれぞれの標準偏差を算出し、
【００４４】
【数２２】
によりτにおけるコヒーレンス関数の絶対値を求める。
この作業は、使用する光ファイバ１５の全長の光の往復時間よりも小さいτの場合であっても可能であり、その範囲でのコヒーレンス関数の絶対値を求めることができる。その値が1/eに減ずるτを求めれば、光源のコヒーレンス時間を求めることもできる。
【００４５】
但し、標準偏差の代わりに変動係数（標準偏差／標準偏差の算出に使用したポイント数）を用いてもコヒーレンス関数を求めることができる。
上記実施形態によるコヒーレンス関数の測定方法並びにコヒーレンス時間の測定方法は、従来技術に対して以下の優位性を持つ。
【００４６】
第１に、コヒーレンス時間と同程度の光ファイバ長により測定が可能な上、仮にコヒーレンス時間が光ファイバ長よりもはるかに長かった場合には、そのこと自体が測定により明らかにされるので、より長い光ファイバに置き換えて測定をやり直すなどの処置を取ることができる。よって、本実施形態は、従来技術よりも信頼性に優れた測定結果を提供するといえる。
【００４７】
第２に、単にコヒーレンス時間の測定が可能となるだけでなく、そのコヒーレンス関数の絶対値の全容が掌握される。これにより、レーザ光の性質をより詳しく測定できると考えられる。
以上、実施形態に記載の構成によれば、スペクトル線幅が狭い（コヒーレンス時間が長い）レーザ光のコヒーレンス時間を正確に測定することのできるレーザ光測定方法及びその測定装置が実現する。
【００４８】
尚、本発明は上記実施形態そのままに限定されるものではなく、実施段階ではその要旨を逸脱しない範囲で構成要素を変形して具体化できる。また、上記実施形態に開示されている複数の構成要素の適宜な組み合わせにより、種々の発明を形成できる。例えば、実施形態に示される全構成要素から幾つかの構成を削除してもよい。さらに、異なる実施形態例に亘る構成要素を適宜組み合わせてもよい。
【符号の説明】
【００４９】
１１…被測定レーザ光源、１２…光周波数掃引装置、１３１…光分岐器、１３２…光合波器、１４…光サーキュレータ、１５…光ファイバ、１４…バランス型受光素子、１７…データ取得部、１８…フーリエ変換部、１９…コヒーレンス関数解析部。
【特許請求の範囲】
【請求項１】
遅延時間τの関数である被測定レーザ光のコヒーレンス関数γ(τ)を測定する方法であって、
前記被測定レーザ光の光周波数を線形に掃引し、
前記光周波数が掃引されたレーザ光を分岐し、
前記分岐された一方のレーザ光を光ファイバに入射してレイリー散乱光を生じさせ、
前記光ファイバで生じたレイリー散乱光と前記分岐されたもう一方のレーザ光とを合波し、
前記合波された光を受光して光電流を検出し、
前記検出された光電流を数値化してフーリエ変換し、フーリエ変換された光電流値より、遅延０における振幅の標準偏差と遅延τにおける振幅の標準偏差とを算出し、両者の比により遅延τにおけるコヒーレンス関数の絶対値を求める
ことを特徴とするレーザ光測定方法。
【請求項２】
さらに、前記コヒーレンス関数の絶対値が1/eになるτをもってコヒーレンス時間を求めることを特徴とする請求項１に記載のレーザ光測定方法。
【請求項３】
遅延時間τの関数であるレーザ光のコヒーレンス関数γ(τ)を測定する装置であって、
被測定レーザ光の光周波数を線形に掃引する光周波数掃引手段と、
前記光周波数掃引がされたレーザ光を分岐する分岐手段と、
前記分岐された一方のレーザ光を入射してレイリー散乱光を生じさせる光ファイバと、
前記光ファイバで生じたレイリー散乱光と前記分岐されたもう一方のレーザ光とを合波する合波手段と、
前記合波された光を受光して光電流を検出する検出手段と、
前記検出された光電流を数値化してフーリエ変換し、フーリエ変換された光電流値より、遅延０における振幅の標準偏差と遅延τにおける振幅の標準偏差とを算出し、両者の比により遅延τにおけるコヒーレンス関数の絶対値を求める演算手段と
を具備することを特徴とするレーザ光測定装置。
【請求項４】
前記演算手段は、さらに前記コヒーレンス関数の絶対値が1/eになるτをもってコヒーレンス時間を求めることを特徴とする請求項３に記載のレーザ光測定装置。
【請求項１】
遅延時間τの関数である被測定レーザ光のコヒーレンス関数γ(τ)を測定する方法であって、
前記被測定レーザ光の光周波数を線形に掃引し、
前記光周波数が掃引されたレーザ光を分岐し、
前記分岐された一方のレーザ光を光ファイバに入射してレイリー散乱光を生じさせ、
前記光ファイバで生じたレイリー散乱光と前記分岐されたもう一方のレーザ光とを合波し、
前記合波された光を受光して光電流を検出し、
前記検出された光電流を数値化してフーリエ変換し、フーリエ変換された光電流値より、遅延０における振幅の標準偏差と遅延τにおける振幅の標準偏差とを算出し、両者の比により遅延τにおけるコヒーレンス関数の絶対値を求める
ことを特徴とするレーザ光測定方法。
【請求項２】
さらに、前記コヒーレンス関数の絶対値が1/eになるτをもってコヒーレンス時間を求めることを特徴とする請求項１に記載のレーザ光測定方法。
【請求項３】
遅延時間τの関数であるレーザ光のコヒーレンス関数γ(τ)を測定する装置であって、
被測定レーザ光の光周波数を線形に掃引する光周波数掃引手段と、
前記光周波数掃引がされたレーザ光を分岐する分岐手段と、
前記分岐された一方のレーザ光を入射してレイリー散乱光を生じさせる光ファイバと、
前記光ファイバで生じたレイリー散乱光と前記分岐されたもう一方のレーザ光とを合波する合波手段と、
前記合波された光を受光して光電流を検出する検出手段と、
前記検出された光電流を数値化してフーリエ変換し、フーリエ変換された光電流値より、遅延０における振幅の標準偏差と遅延τにおける振幅の標準偏差とを算出し、両者の比により遅延τにおけるコヒーレンス関数の絶対値を求める演算手段と
を具備することを特徴とするレーザ光測定装置。
【請求項４】
前記演算手段は、さらに前記コヒーレンス関数の絶対値が1/eになるτをもってコヒーレンス時間を求めることを特徴とする請求項３に記載のレーザ光測定装置。
【図１】
【図２】
【図２】
【公開番号】特開２０１２−１１２６５７（Ｐ２０１２−１１２６５７Ａ）
【公開日】平成２４年６月１４日（２０１２．６．１４）
【国際特許分類】
【出願番号】特願２０１０−２５９１５６（Ｐ２０１０−２５９１５６）
【出願日】平成２２年１１月１９日（２０１０．１１．１９）
【出願人】（０００００４２２６）日本電信電話株式会社 (13,992)
【Ｆターム（参考）】
【公開日】平成２４年６月１４日（２０１２．６．１４）
【国際特許分類】
【出願日】平成２２年１１月１９日（２０１０．１１．１９）
【出願人】（０００００４２２６）日本電信電話株式会社 (13,992)
【Ｆターム（参考）】
[ Back to top ]