(57)【特許請求の範囲】
【請求項１】
プローブ光を生成する光源部と、
前記プローブ光により測定対象となる光ファイバで発生する後方ブリルアン散乱光を、第１光路及び第２光路に２分岐する分岐部と、
前記第１光路に設けられた、第１周波数の周波数シフトを与える第１光周波数シフタ部と、
前記第２光路に設けられた、第２周波数の周波数シフトを与える第２光周波数シフタ部と、
前記第１光路及び前記第２光路のいずれか一方に設けられた遅延部と、
前記第１光路及び前記第２光路を伝播する光を合波して合波光を生成する合波部と、
前記合波光をヘテロダイン検波して差周波を第１電気信号として出力するコヒーレント検波部と、
前記第１電気信号と同じ周波数を持つ第２電気信号を生成する電気信号生成部と、
前記第１電気信号が２分岐された一方と前記第２電気信号とをホモダイン検波して、周波数シフト量δνＳＤＨを取得する周波数シフト量取得部と、
前記第１電気信号から後方ブリルアン散乱光の強度情報δＰＢ／ＰＢを取得する信号強度取得部と、
前記周波数シフト量δνＳＤＨ及び前記強度情報δＰＢ／ＰＢから歪みδε及び温度変化δＴを分離して取得する信号処理部と
を備えることを特徴とする光ファイバ歪み及び温度測定装置。
【請求項２】
前記信号処理部は、
前記周波数シフト量δνＳＤＨ及び前記強度情報δＰＢ／ＰＢと、予め求めておいた、光ファイバ中の後方ブリルアン散乱の周波数シフトの歪み依存係数Ｃνε、温度依存係数ＣνＴ、及び後方ブリルアン散乱の散乱係数の歪み依存係数ＣＰε、温度依存係数ＣＰＴとから、下記の２元連立方程式（１）を解くことにより、光ファイバ中の歪みδεと温度変化δＴを取得する
ことを特徴とする請求項１に記載の光ファイバ歪み及び温度測定装置。
【数１】
0006358277.tif 000012
【請求項３】
プローブ光を生成する過程と、
前記プローブ光により測定対象となる光ファイバで発生する後方ブリルアン散乱光を、第１光路及び第２光路に２分岐する過程と、
前記第１光路を伝播する光に、第１周波数の周波数シフトを与える過程と、
前記第２光路を伝播する光に、第２周波数の周波数シフトを与える過程と、
前記第１光路及び前記第２光路のいずれか一方を伝搬する光に遅延を与える過程と、
前記第１光路及び前記第２光路を伝播する光を合波して合波光を生成する過程と、
前記合波光をヘテロダイン検波して差周波を第１電気信号として出力する過程と、
前記第１電気信号と同じ周波数を持つ第２電気信号を生成する過程と、
前記第１電気信号が２分岐された一方と前記第２電気信号とをホモダイン検波して、周波数シフト量δνＳＤＨを取得する過程と、
前記第１電気信号から後方ブリルアン散乱光の強度情報δＰＢ／ＰＢを取得する過程と、
前記周波数シフト量δνＳＤＨ及び前記強度情報δＰＢ／ＰＢから歪みδε及び温度変化δＴを分離して取得する過程と
を備えることを特徴とする光ファイバ歪み及び温度測定方法。
【請求項４】
前記歪みδε及び温度変化δＴを、前記周波数シフト量δνＳＤＨ及び前記強度情報δＰＢ／ＰＢと、予め求めておいた、光ファイバ中の後方ブリルアン散乱の周波数シフトの歪み依存係数Ｃνε、温度依存係数ＣνＴ、及び後方ブリルアン散乱の散乱係数の歪み依存係数ＣＰε、温度依存係数ＣＰＴとから、下記の２元連立方程式（１）を解くことにより取得する
ことを特徴とする請求項３に記載の光ファイバ歪み及び温度測定方法。
【数２】
0006358277.tif 000013
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、ブリルアン散乱光を用いた、光ファイバ歪み及び温度測定装置並びに光ファイバ歪み及び温度測定方法に関する。
【背景技術】
【０００２】
光ファイバ通信の発展とともに、光ファイバ自体をセンシング媒体とする分布型光ファイバセンシングが盛んに研究されている。分布型光ファイバセンシングでは、光ファイバの片端から光パルスを入射し、光ファイバ中で後方散乱された光を時間に対して測定する時間領域リフレクトメトリ（ＯＴＤＲ：Ｏｐｔｉｃａｌ  Ｔｉｍｅ  Ｄｏｍａｉｎ  Ｒｅｆｌｅｃｔｏｍｅｔｒｙ）が代表的である。光ファイバ中の後方散乱には、レイリー散乱、ブリルアン散乱及びラマン散乱がある。この中で自然ブリルアン散乱を測定するものはＢＯＴＤＲ（Ｂｒｉｌｌｏｕｉｎ  ＯＴＤＲ）と呼ばれる（例えば、非特許文献１参照）。
【０００３】
ブリルアン散乱は、光ファイバに入射される光パルスの中心周波数に対して、ストークス側及び反ストークス側にＧＨｚ程度周波数シフトした位置に観測され、そのスペクトルはブリルアン利得スペクトル（ＢＧＳ：Ｂｒｉｌｌｏｕｉｎ  Ｇａｉｎ  Ｓｐｅｃｔｒｕｍ）と呼ばれる。ＢＧＳの周波数シフト及びスペクトル線幅は、それぞれブリルアン周波数シフト（ＢＦＳ：Ｂｒｉｌｌｏｕｉｎ  Ｆｒｅｑｕｅｎｃｙ  Ｓｈｉｆｔ）及びブリルアン線幅と呼ばれる。ＢＦＳ及びブリルアン線幅は、光ファイバの材質および入射光波長によって異なる。例えば、石英系のシングルモード光ファイバの場合、波長１．５５μｍにおけるＢＦＳの大きさ及びブリルアン線幅は、それぞれ約１１ＧＨｚ及び約３０ＭＨｚとなることが報告されている。また、非特許文献１からシングルモードファイバ中の歪み、温度の変化に伴うＢＦＳの大きさは波長１．５５μｍにおいて、それぞれ０．０４９ＭＨｚ／με、１．０ＭＨｚ／℃である。
【０００４】
ここで、ＢＦＳは歪みと温度に対して依存性を持つため、ＢＯＴＤＲは橋梁やトンネルなどに代表される大型建造物や、地滑りが発生する恐れのある箇所などの監視目的に注目されている。  
【０００５】
ＢＯＴＤＲは、光ファイバ中で発生する自然ブリルアン散乱光のスペクトル波形を測定するため、別途用意した参照光とのヘテロダイン検波を行うのが一般的である。自然ブリルアン散乱光の強度はレイリー散乱光の強度に比べて２〜３桁小さい。このため、ヘテロダイン検波は最小受光感度を向上させる上でも有用となる。
【０００６】
図１１を参照して、従来のＢＯＴＤＲについて説明する（例えば、特許文献１参照）。図１１は、従来の光ファイバ歪み測定装置の模式的なブロック図である。
【０００７】
光源１１２から出射された連続光は、光カプラ１４２によって、２分岐される。２分岐された一方は、参照光として用いられ、他方は、光周波数シフタ１４３によってブリルアン周波数に相当する周波数シフトを受けた後、光パルス発生器１１４により、パルス状のプローブ光となる。
【０００８】
このプローブ光は、光カプラ１２０を経て測定対象となる光ファイバ（被測定光ファイバ）１００に入射される。被測定光ファイバ１００からの後方ブリルアン散乱光は、光カプラ１５０において参照光と合波された後、バランス型フォトダイオード（ＰＤ）１６２及びＦＥＴ増幅器１６４からなるレシーバ１６０によってヘテロダイン検波される。
【０００９】
ここで、プローブ光は、光周波数シフタ１４３によってブリルアン周波数程度の周波数シフトが施されているため、ヘテロダイン検波されて生成されるビート信号の周波数は低くなる。これによって、受信側では低域なＰＤ１６２及びＦＥＴ増幅器１６４を使用することが可能になる。ビート信号をミキサー１７０、電気フィルタ１７８により周波数をダウンシフトさせた後、検波回路１７２により２乗検波もしくは包絡線検波することにより得られるＩＦ（Ｉｎｔｅｒｍｅｄｉａｔｅ  Ｆｒｅｑｕｅｃｙ）信号のパワーや振幅を測定する。この結果は、信号処理装置１７４に送られる。
【００１０】
なお、ＢＯＴＤＲは、光ファイバの長手方向に対する周波数スペクトル分布の情報を扱うため、時間、振幅及び周波数の３次元の情報を取得する必要がある。図１２を参照してＢＯＴＤＲにおいて、時間、振幅及び周波数の３次元の情報の取得方法について説明する。図１２は、従来の光ファイバ歪み測定装置における時間、振幅及び周波数の３次元の情報の取得方法を説明するための模式図である。上述の特許文献１に開示の技術では、ブリルアン周波数スペクトル全体を測定するには、時間ｔ及び振幅Ｉの２次元情報を、局発電気信号源１８３の周波数ｆを掃引して取得する。  
【００１１】
ここで、ＢＯＴＤＲに限らずブリルアン散乱を用いた分布型光ファイバセンシングでは、上述のように、歪み及び温度の両方に対してＢＦＳが生じる。従って、歪みと温度を区別することは必須の課題である。この課題に対して、光ファイバ中の後方ブリルアン散乱係数の、歪み依存係数及び温度依存係数を利用した方法が提案されている（例えば、非特許文献２又は３参照）。
【００１２】
後方ブリルアン散乱では、周波数シフトだけでなく、その散乱係数も温度、歪み依存性を有することが報告されている。ＢＦＳの歪み依存係数及び温度依存係数をそれぞれＣνε及びＣνＴとし、ブリルアン散乱係数の歪み依存係数及び温度依存係数をそれぞれＣＰε及びＣＰＴとして、これらの係数をあらかじめ測定しておけば、次式で示す２元連立方程式（ａ）を解くことによって歪と温度の分離が可能となる。
【００１３】
【数１】
0006358277.tif 000002
【００１４】
ここで、δνＢはＢＦＳの大きさであり、δＰＢ／ＰＢは、ブリルアン散乱強度の相対変化である。これら、δνＢ及びδＰＢ／ＰＢは、ＢＯＴＤＲで測定される値である。また、δε及びδＴは、それぞれ、歪み及び温度の変化量である。
【先行技術文献】
【特許文献】
【００１５】
【特許文献１】特開２００１−１６５８０８号公報
【非特許文献】
【００１６】
【非特許文献１】Ｔ．Ｋｕｒａｓｈｉｍａ  ｅｔ  ａｌ．，“Ｂｒｉｌｌｏｕｉｎ  Ｏｐｔｉｃａｌ−ｆｉｂｅｒ  ｔｉｍｅ  ｄｏｍａｉｎ  ｒｅｆｌｅｃｔｏｍｅｔｒｙ”，ＩＥＩＣＥ  Ｔｒａｎｓ．  Ｃｏｍｍｕｎ．，  ｖｏｌ．Ｅ７６−Ｂ，  ｎｏ．４，  ｐｐ．３８２−３９０  （１９９３）
【非特許文献２】Ｔ．Ｒ．Ｐａｒｋｅｒ  ｅｔ  ａｌ．，“Ｓｉｍｕｌｔａｎｅｏｕｓ  ｄｉｓｔｒｉｂｕｔｅｄ  ｍｅａｓｕｒｅｍｅｎｔ  ｏｆ  ｓｔｒａｉｎ  ａｎｄ  ｔｅｍｐｅｒａｔｕｒｅ  ｆｒｏｍ  ｎｏｉｓｅ−ｉｎｉｔｉａｔｅｄ  Ｂｒｉｌｌｏｕｉｎ  ｓｃａｔｔｅｒｉｎｇ  ｉｎ  ｏｐｔｉｃａｌ  ｆｉｂｅｒｓ”，ＩＥＥＥ  Ｊ．Ｑｕａｎｔｕｍ  Ｅｌｅｃｔｒｏｎ．，  ｖｏｌ．３４，  Ｎｏ．４，  ｐｐ．６４５−６５９  （１９９８）
【非特許文献３】Ｙ．Ｓａｋａｉｒｉ  ｅｔ  ａｌ．，“Ａｓｙｓｔｅｍ  ｆｏｒ  ｍｅａｓｕｒｉｎｇ  ｔｅｍｐｅｒａｔｕｒｅ  ａｎｄ  ｓｔｒａｉｎ  ｓｅｐａｒａｔｅｌｙ  ｂｙ  ＢＯＴＤＲ  ａｎｄ  ＯＴＤＲ”，  Ｐｒｏｃｅｅｄｉｎｇ  ｏｆ  ＳＰＩＥ，  ｖｏｌ．４９２０，  ｐｐ．２７４−２８４  （２００２）
【発明の概要】
【発明が解決しようとする課題】
【００１７】
ここで、自然ブリルアン散乱光は非常に微弱なため、ヘテロダイン検波を適用しても十分な信号雑音比（Ｓ／Ｎ）を確保できない。その結果、Ｓ／Ｎ改善のための平均化処理が必要となる。この平均化処理と上述の３次元情報の取得のため、従来の光ファイバ歪み測定装置では、測定時間の短縮が難しい。
【００１８】
この発明の発明者は、上述の問題点に鑑みて、光の周波数変化をコヒーレント検波により与えられるビート信号の位相差として測定することにより、時間及び位相の２次元の情報を取得する、自然ブリルアン散乱光を用いた、光ファイバ歪み測定装置及び光ファイバ歪み測定方法を検討し、その検討結果の一部を特許出願（特願２０１５−０７２５４６号：以後、「先の出願」という。）している。
【００１９】
この先の出願の光ファイバ歪み測定装置及び光ファイバ歪み測定方法によれば、自己遅延ヘテロダイン型のＢＯＴＤＲ（ＳＤＨ−ＢＯＴＤＲ：Ｓｅｌｆ−Ｄｅｌａｙｅｄ  Ｈｅｔｅｒｏｄｙｎｅ  ＢＯＴＤＲ）の技術を用いて、光の周波数変化をコヒーレント検波により与えられるビート信号の位相差として測定することにより、時間及び位相の２次元の情報を取得する。このＳＤＨ−ＢＯＴＤＲでは、周波数掃引を必要としないため、３次元の情報の取得が必要な従来技術に比べて、測定時間が短縮される。
【００２０】
しかしながら、ＳＤＨ−ＢＯＴＤＲでは、ＢＦＳに起因する位相差だけでなく、ブリルアン散乱係数の変化に伴う強度変化も位相差として捉えてしまう。すなわち、ＳＤＨ−ＢＯＴＤＲで観測される周波数シフトには、ＢＦＳの歪み依存係数Ｃνε及び温度依存係数ＣνＴだけでなく、ブリルアン散乱係数の歪み及び温度の変化に起因する変化分が重畳されている。このため、従来のＢＯＴＤＲと同様の手法で、歪みと温度変化を分離して取得することができない。
【００２１】
この発明は、上述の問題点に鑑みてなされたものであり、この発明の目的は、ＳＤＨ−ＢＯＴＤＲにおいて、歪みと温度変化を分離して取得することが可能な、光ファイバ歪み及び温度測定装置並びに光ファイバ歪み及び温度測定方法を提供することにある。
【課題を解決するための手段】
【００２２】
上述した目的を達成するために、この発明の光ファイバ歪み及び温度測定装置は、光源部と、分岐部と、第１光周波数シフタ部と、第２光周波数シフタ部と、遅延部と、合波部と、コヒーレント検波部と、電気信号生成部と、周波数シフト量取得部と、信号強度取得部と、信号処理部とを備えて構成される。
【００２３】
光源部は、プローブ光を生成する。プローブ光は、測定対象となる光ファイバ（被測定光ファイバ）に入射される。分岐部は、プローブ光により被測定光ファイバで発生する後方ブリルアン散乱光を、第１光路及び第２光路に２分岐する。第１光周波数シフタ部は、第１光路に設けられており、第１周波数の周波数シフトを与える。第２光周波数シフタ部は、第２光路に設けられており、第２周波数の周波数シフトを与える。遅延部は、第１光路及び第２光路のいずれか一方に設けられており、第１光路及び第２光路を伝播する光の間に遅延時間差を与える。合波部は、第１光路及び第２光路を伝播する光を合波して合波光を生成する。コヒーレント検波部は、合波光をヘテロダイン検波して差周波を第１電気信号として出力する。電気信号生成部は、第１周波数と第２周波数の差周波として生成される、いわゆるビート信号である、第２電気信号を生成する。周波数シフト量取得部は、第１電気信号が２分岐された一方と第２電気信号とをホモダイン検波して、周波数シフト量δνＳＤＨを取得する。この第１電気信号は、いわゆるビート信号である。信号強度取得部は、第１電気信号から後方ブリルアン散乱光の強度情報δＰＢ／ＰＢを取得する。信号処理部は、周波数シフト量δνＳＤＨ及び強度情報δＰＢ／ＰＢから歪みδε及び温度変化δＴを分離して取得する。
【００２６】
この光ファイバ歪み及び温度測定装置の好適実施形態によれば、信号処理部は、周波数シフト量δνＳＤＨ及び強度情報δＰＢ／ＰＢと、予め求めておいた、光ファイバ中の後方ブリルアン散乱の周波数シフトの歪み依存係数Ｃνε、温度依存係数ＣνＴ、及び後方ブリルアン散乱の散乱係数の歪み依存係数ＣＰε及び温度依存係数ＣＰＴとから、以下の２元連立方程式（１）を解くことにより、光ファイバ中の歪みδεと温度変化δＴを取得する。ただし、被測定ファイバ中の伝送損失は無視するものとする。
【００２７】
【数２】
0006358277.tif 000003
【００２８】
また、この発明の光ファイバ歪み及び温度測定方法は、以下の過程を備えて構成される。
【００２９】
先ず、プローブ光を生成する。プローブ光は、被測定光ファイバに入射される。次に、プローブ光により被測定光ファイバで発生する後方ブリルアン散乱光を、第１光路及び第２光路に２分岐する。次に、第１光路を伝播する光に対して、第１周波数の周波数シフトを与え、第２光路を伝播する光に対して、第２周波数の周波数シフトを与える。次に、第１光路及び第２光路を伝播する光の間に遅延時間差を与える。次に、第１光路及び第２光路を伝播する光を合波して合波光を生成する。次に、合波光をヘテロダイン検波して第１電気信号を生成する。次に、第１周波数と第２周波数の差周波として生成される、いわゆるビート信号である、第２電気信号を生成する。次に、第１電気信号が２分岐された一方と第２電気信号を、ホモダイン検波して、周波数シフト量δνＳＤＨを取得する。また、第１電気信号から後方ブリルアン散乱光の強度情報δＰＢ／ＰＢを取得する。次に、周波数シフト量δνＳＤＨ及び強度情報δＰＢ／ＰＢから歪みδε及び温度変化δＴを分離して取得する。
【００３１】
この光ファイバ歪み及び温度測定方法の実施に当たり、好適には、歪みδε及び温度変化δＴを、周波数シフト量δνＳＤＨ及び強度情報δＰＢ／ＰＢと、予め求めておいた、光ファイバ中の後方ブリルアン散乱の周波数シフトの歪み依存係数Ｃνε、温度依存係数ＣνＴ、及び後方ブリルアン散乱の散乱係数の歪み依存係数ＣＰε、温度依存係数ＣＰＴとから、上記２元連立方程式（１）を解くことにより取得する。
【発明の効果】
【００３２】
この発明の光ファイバ歪み及び温度測定装置及び光ファイバ歪み及び温度測定方法によれば、光の周波数変化をコヒーレント検波により与えられるビート信号の位相差として測定することにより、時間及び位相の２次元の情報を取得する。このため、３次元の情報の取得が必要な従来技術に比べて、測定時間が短縮される。
【００３３】
また、ビート信号の位相差と強度の両方の情報を用いることで、歪みδε及び温度変化δＴを取得することができる。
【図面の簡単な説明】
【００３４】
【図１】光ファイバ歪み測定装置の基本構成の模式的なブロック図である。
【図２】ブリルアンシフトとビート信号の位相変化を示す模式図である。
【図３】光ファイバに与えた歪み変化及び温度変化を示す図である。
【図４】基本構成の光ファイバ歪み測定装置で得られるビート信号波形を示す図である。
【図５】ビート信号に対して、ホモダイン検波及びフィルタ処理を施すことにより得られた、周波数シフト量を示す図である。
【図６】第１光ファイバ歪み及び温度測定装置の模式的なブロック図である。
【図７】δＰＢ／ＰＢを示す図である。
【図８】分離後のδεを示す図である。
【図９】分離後のδＴを示す図である。
【図１０】第２光ファイバ歪み及び温度測定装置の模式的なブロック図である。
【図１１】従来の光ファイバ歪み測定装置の模式的なブロック図である。
【図１２】従来の光ファイバ歪み測定装置における時間、振幅及び周波数の３次元の情報の取得方法を説明するための模式図である。
【発明を実施するための形態】
【００３５】
以下、図を参照して、この発明の実施の形態について説明するが、各図は、この発明が理解できる程度に概略的に示したものに過ぎない。また、以下、この発明の好適な構成例につき説明するが、単なる好適例にすぎない。従って、この発明は以下の実施の形態に限定されるものではなく、この発明の構成の範囲を逸脱せずにこの発明の効果を達成できる多くの変更又は変形を行うことができる。
【００３６】
（基本構成）
この発明の理解に資するため、図１を参照して、光の周波数変化をコヒーレント検波により与えられるビート信号の位相差として測定する、光ファイバ歪み測定装置の基本構成について説明する。図１は、光ファイバ歪み測定装置の基本構成の模式的なブロック図である。
【００３７】
光ファイバ歪み測定装置は、光源部１０、サーキュレータ２０、光増幅器３０、光バンドパスフィルタ３２、自己遅延ヘテロダイン干渉計４１及びタイミング制御器９０を備えて構成される。
【００３８】
光源部１０は、プローブ光を生成する。光源部１０は、連続光を生成する光源１２と、連続光から光パルスを生成する光パルス発生器１４を備えて構成される。
【００３９】
ここで、基本構成の光ファイバ歪み測定装置は、周波数変化に応じた位相差を測定する。このため、光源１２の周波数揺らぎ及び周波数スペクトル線幅（以下、単に線幅とも称する。）は、ブリルアンシフトよりも十分に小さくなければならない。そこで、光源１２として周波数安定化狭線幅光源が用いられる。例えば、測定対象となる光ファイバ（以下、被測定光ファイバとも称する。）１００の歪みを０．００８％としたとき、ブリルアンシフトは４ＭＨｚに相当する。このため、０．００８％程度の歪みを測定するには、光源１２の周波数揺らぎ及び線幅は４ＭＨｚより十分に小さく、数１０ｋＨｚ以下であることが望ましい。なお、周波数揺らぎ及び線幅が１０ｋＨｚ程度若しくはそれ以下の狭線幅レーザが、既製品として一般に入手可能である。
【００４０】
光パルス発生器１４は、任意好適な従来周知の、音響光学（ＡＯ：Ａｃｏｕｓｔ  Ｏｐｔｉｃａｌ）変調器又は電気光学（ＥＯ：Ｅｌｅｃｔｒｉｃ  Ｏｐｔｉｃａｌ）変調器を用いて構成される。光パルス発生器１４は、タイミング制御器９０で生成された電気パルスに応じて、連続光から光パルスを生成する。この光パルスの繰り返し周期は、被測定光ファイバ１００を光パルスが往復するのに要する時間よりも長く設定される。この光パルスが、プローブ光として、光源部１０から出力される。
【００４１】
この光源部１０から出力されたプローブ光は、サーキュレータ２０を経て、被測定光ファイバ１００に入射される。なお、サーキュレータ２０に換えて、光カプラを用いても良い。
【００４２】
被測定光ファイバ１００からの後方散乱光は、サーキュレータ２０を経て、例えば、エルビウム添加光ファイバ増幅器（ＥＤＦＡ）などで構成される光増幅器３０に送られる。光増幅器３０で増幅された後方散乱光は、光バンドパスフィルタ３２に送られる。光バンドパスフィルタ３２は、１０ＧＨｚ程度の透過帯域を有しており、自然ブリルアン散乱光のみを透過する。この自然ブリルアン散乱光は、自己遅延ヘテロダイン干渉計４１に送られる。この光バンドパスフィルタ３２から出射される自然ブリルアン散乱光の時刻ｔにおける信号Ｅ０（ｔ）は、以下の式（２）で表される。
【００４３】
【数３】
0006358277.tif 000004
【００４４】
ここで、Ａ０は振幅、ηＢ（ｔ）はブリルアン散乱係数、ｆＢ（ｔ）はブリルアン散乱光の光周波数、φ０は初期位相を示している。なお、ブリルアン散乱係数ηＢ（ｔ）及びブリルアン散乱光の光周波数ｆＢ（ｔ）は、光ファイバ中での局所的な歪みや温度変化により変化するため、時間ｔの関数としている。
【００４５】
自己遅延ヘテロダイン干渉計４１は、分岐部４２、光周波数シフタ部４３、遅延部４８、合波部５０、コヒーレント検波部６０、電気信号生成部８０及び信号処理装置７４を備えて構成される。
【００４６】
電気信号生成部８０の局発電気信号源８３は、周波数ｆＡＯＭの電気信号を生成する。
【００４７】
分岐部４２は、プローブ光により被測定光ファイバ１００で発生する後方ブリルアン散乱光を、光バンドパスフィルタ３２を経て受け取り、第１光路及び第２光路に２分岐する。
【００４８】
光周波数シフタ部４３は、第１光路に設けられている。光周波数シフタ部４３は、局発電気信号源８３で生成された周波数ｆＡＯＭの電気信号を用いて、第１光路を伝播する光に対して、周波数ｆＡＯＭの周波数シフトを与える。
【００４９】
従来の例えば特許文献１に開示されている測定装置では、ＢＦＳに対応する数十ＧＨｚ程度の周波数シフトを与える。これに対し、この光ファイバ歪み測定装置では、周波数ｆＡＯＭは、数十ＭＨｚ程度である。このため、従来の測定装置に比べて周波数シフタとして小型でかつ安価なものを用いることができる。
【００５０】
また、この構成例では、第２光路に遅延部４８が設けられている。遅延部４８は、第２光路を伝播する光に時間τの遅延を与える。
【００５１】
合波部５０は、第１光路及び第２光路を伝播する光を合波して合波光を生成する。合波部５０に入射される、第１光路を伝播する光信号Ｅ１（ｔ）及び第２光路を伝播する光信号Ｅ２（ｔ−τ）は、それぞれ、以下の式（３）及び（４）で表される。
【００５２】
【数４】
0006358277.tif 000005
【００５３】
ここで、Ａ１及びＡ２は、それぞれＥ１（ｔ）及びＥ２（ｔ−τ）の振幅であり、φ１及びφ２は、それぞれＥ１（ｔ）及びＥ２（ｔ−τ）の初期位相である。
【００５４】
コヒーレント検波部６０は、合波光をヘテロダイン検波してビート信号を生成する。コヒーレント検波部６０は、例えば、バランス型フォトダイオード（ＰＤ）６２とＦＥＴ増幅器６４を備えて構成される。ヘテロダイン検波により与えられるビート信号Ｉ１２は、以下の式（５）で表される。
【００５５】
【数５】
0006358277.tif 000006
【００５６】
コヒーレント検波部６０で生成されたビート信号Ｉ１２は第１電気信号としてミキサー部７０に送られる。また、局発電気信号源８３で生成された電気信号は第２電気信号としてミキサー部７０に送られる。
【００５７】
ミキサー部７０は、第１電気信号と、第２電気信号とをホモダイン検波して、ホモダイン信号を生成する。局発電気信号源８３で生成された電気信号ＩＡＯＭを以下の式（６）で表す。
【００５８】
【数６】
0006358277.tif 000007
【００５９】
ミキサー部７０で生成されたホモダイン信号は、上記式（５）及び式（６）を乗算して得られる、以下の式（７）で表される。
【００６０】
【数７】
0006358277.tif 000008
【００６１】
上記式（７）中の和周波成分をローパスフィルタ（ＬＰＦ）７２で除去すると、以下の式（８）で表される信号が得られる。
【００６２】
【数８】
0006358277.tif 000009
【００６３】
上記式（８）のφ1−φ2−φＡＯＭと、遅延時間τは一定であるため、ブリルアン周波数の変化ｆＢ（ｔ）のみが出力強度の差として出力される。
【００６４】
ブリルアン周波数ｆＢ（ｔ）は、光源１２の発振周波数の揺らぎと被測定光ファイバ１００の歪みの２つの要因によって変化する。しかし、光源１２として周波数安定化狭線幅光源を用いることで、被測定光ファイバ１００の歪みによる影響が支配的となる。
【００６５】
図２（Ａ）及び（Ｂ）は、ブリルアンシフトとビート信号の位相変化を示す模式図である。図２（Ａ）は、横軸に時間ｔを取って示し、縦軸に周波数を取って示している。また、図２（Ｂ）は、横軸に時間ｔを取って示し、縦軸に、電圧を取って示している。
【００６６】
この横軸の時間は、ブリルアン散乱が起こった場所を示している。すなわち、プローブ光が出射された時間に対して、時間ｔ経過後に後方ブリルアン散乱光が入射された場合、被測定光ファイバ内の光の伝播速度をｖとすると、被測定光ファイバの入射端からｖｔ／２の位置で後方ブリルアン散乱が生じたことになる。
【００６７】
図２（Ａ）及び（Ｂ）では、時刻ｔ１からｔ２までの時間Ｔに対応する区間において、周波数シフトが生じた例を示している。このとき、自己遅延ヘテロダイン干渉計で遅延時間τが与えられているため、位相変化はｔ１からｔ１＋τまでの間に変化し、時刻ｔ２からｔ２＋τまでの間に元の状態に戻る。すなわち、光ファイバ歪み測定装置で位相差を測定するには、Ｔ≧τの関係を満たす必要があり、測定可能な時間分解能（すなわち、空間分解能）がτによって定まる。さらに測定可能な周波数変化もτの大きさで定まる。すなわち、τが大きくなると、測定可能な周波数範囲が小さくなるが、空間分解能は大きくなる。一方、τが小さくなると、空間分解能は小さくなるが、測定可能な周波数範囲が大きくなる。このように、遅延時間と測定可能な周波数の間にトレードオフの関係がある。
【００６８】
図３〜５を参照して、光ファイバに歪み変化及び温度変化を与えたときの、ＳＤＨ−ＢＯＴＤＲで得られる周波数シフトについて説明する。
【００６９】
図３は、光ファイバに与えた歪み変化及び温度変化を示す図である。図３では、横軸に光ファイバにおける位置［単位：ｍ］を取って示し、縦軸に歪み変化［単位：με］及び温度変化［単位：℃］を取って示している。ここでは、光ファイバの長さを１ｋｍとし、光ファイバ歪み測定装置側の端部から３００〜３２０ｍの区間と９４０〜９６０ｍの区間に、２００μεの歪みを付与し、６２０〜６４０ｍの区間と９４０〜９６０ｍの区間に、２０℃の温度変化を付与した場合について考える。
【００７０】
図４は、上述の場合に、光ファイバ歪み測定装置で得られるビート信号波形を示す図である。図４（Ａ）及び図４（Ｂ）では、横軸に光ファイバにおける位置［単位：ｍ］を取って示し、縦軸に信号強度［単位：Ｖ］で取って示している。図４（Ｂ）は、図４（Ａ）のＩで示す部分を拡大して示している。
【００７１】
ここでは、ＢＦＳの歪み依存係数Ｃνεを０．０４９ＭＨｚ／μεとし、温度依存係数ＣνＴを１．０ＭＨｚ／℃としている。また、ブリルアン散乱係数の歪み依存係数ＣＰεを、−７．７×１０−４％／℃とし、温度依存係数ＣＰＴを０．３６％／℃としている。
【００７２】
図４（Ｂ）に示すように、ビート信号中には歪み、温度変化を付与した箇所においてＢＦＳに起因する位相シフト（図４（Ｂ）中、ＩＩで示す部分）とブリルアン散乱係数に起因する強度変化（図４（Ｂ）中、ＩＩＩで示す部分）が生じる。
【００７３】
図５は、このビート信号に対して、ホモダイン検波及びフィルタ処理を施すことにより得られた、周波数シフト量δνＳＤＨを示す図である。図５では、横軸に光ファイバにおける位置［単位：ｍ］を取って示し、縦軸に周波数シフト量δνＳＤＨ［単位：ＭＨｚ］を取って示している。
【００７４】
周波数シフト量δνＳＤＨは上記式（１）で与えられる周波数シフト量δνＢにブリルアン散乱係数の変化に起因する強度変化分が重畳した値で以下の式（９）を用いて計算される。
【００７５】
【数９】
0006358277.tif 000010
【００７６】
基本構成の光ファイバ歪み測定装置は、ブリルアン散乱の周波数変化を測定する構成である。このため、基本構成の光ファイバ歪み測定装置では、得られた情報から歪みと温度変化を分離するのは困難である。
【００７７】
（第１実施形態）
図６を参照して、第１実施形態の光ファイバ歪み及び温度測定装置（以下、第１測定装置とも称する。）について説明する。図６は、第１測定装置の模式的なブロック図である。なお、以下の説明において、光ファイバ歪み測定装置の基本構成と重複する説明を省略することがある。
【００７８】
第１測定装置では、コヒーレント検波部６０の出力が２分岐され、一方が、ミキサー部７０、ＬＰＦ７２を経て信号処理装置７５に送られ、他方が、２乗回路９２、ローパスフィルタ（ＬＰＦ）９４、及び１／２乗回路９６を通って信号処理装置７５に送られる点が基本構成の光ファイバ歪み測定装置と異なっている。
【００７９】
２乗回路９２、ＬＰＦ９４、及び１／２乗回路９６は、信号強度取得部９１を構成し、コヒーレント検波部６０の出力であるビート信号の包絡線検波の機能を実現する。この結果、ビート信号が、２乗回路９２、ＬＰＦ９４、及び１／２乗回路９６を順に通過すると、ビート信号の強度情報のみが得られる。すなわち、２乗回路９２、ＬＰＦ９４、及び１／２乗回路９６から、上記式（１）におけるδＰＢ／ＰＢが得られる。
【００８０】
図７は、このビート信号に対して、包絡線検波により得られた、δＰＢ／ＰＢを示す図である。図７では、横軸に光ファイバにおける位置［単位：ｍ］を取って示し、縦軸にδＰＢ／ＰＢ［単位：％］を取って示している。
【００８１】
また、ミキサー部７０、ＬＰＦ７２は、周波数シフト量取得部７１を構成する。周波数シフト量取得部７１から、図５に示した、上記式（９）におけるδνＳＤＨが得られる。
【００８２】
信号処理装置では、ビート信号から得られたδνＳＤＨとδＰＢ／ＰＢと、予め取得しておいた各係数を用いて、上記式（１）の２元連立方程式を解くことにより、δε及びδＴを取得する。
【００８３】
上記式（１）の２元連立方程式を解くと、δε及びδＴは、それぞれ以下の式（１０）で与えられる。
【００８４】
【数１０】
0006358277.tif 000011
【００８５】
図５に示すδνＳＤＨと、図７に示すδＰＢ／ＰＢから、上記式（１０）により、歪みδε及び温度変化δＴが分離されて求められる。
【００８６】
図８は、分離後のδεを示す図である。図８では、横軸に光ファイバにおける位置［単位：ｍ］を取って示し、縦軸にδε［単位：με］を取って示している。また、図９は、分離後のδＴを示す図である。図９では、横軸に光ファイバにおける位置［単位：ｍ］を取って示し、縦軸にδＴ［単位：℃］を取って示している。
【００８７】
このように、第１測定装置によれば、δε及びδＴをそれぞれ分離して取得することができる。
【００８８】
（第２実施形態）
図１０を参照して、第２実施形態の光ファイバ歪み及び温度測定装置（以下、第２測定装置とも称する。）について説明する。第２測定装置では、第１光路に第１光周波数シフタ部４４が設けられている。また、第２光路に第２光周波数シフタ部４６と遅延部４８が設けられている。
【００８９】
電気信号生成部８１は第１局発電気信号源８２、第２局発電気信号源８４、ミキサー部８６及びローパスフィルタ（ＬＰＦ）８８を備えて構成される。なお、第１局発電気信号源８２、第２局発電気信号源８４は電気信号生成部８１の外部に在っても良い。第１局発電気信号源８２は、第１周波数ｆ１の電気信号を生成する。第２局発電気信号源８４は、第２周波数ｆ２の電気信号を生成する。ミキサー部８６は、第１周波数ｆ１の電気信号と、第２周波数ｆ２の電気信号から、第１周波数ｆ１及び第２周波数ｆ２の和周波数成分と差周波数成分を生成する。ＬＰＦ８８はミキサー部８６で生成される信号から差周波数成分ｆＡＯＭ（＝ｆ１−ｆ２）のビート信号を出力する。
【００９０】
第１光周波数シフタ部４４は、第１光路に設けられている。第１光周波数シフタ部４４は、第１局発電気信号源８２で生成された第１周波数ｆ１の電気信号を用いて、第１光路を伝播する光に対して、第１周波数ｆ１の周波数シフトを与える。
【００９１】
第２光周波数シフタ部４６は、第２光路に設けられている。第２光周波数シフタ部４６は、第２局発電気信号源８４で生成された第２周波数ｆ２の電気信号を用いて、第２光路を伝播する光に対して、第２周波数ｆ２の周波数シフトを与える。
【００９２】
第２測定装置は、自己遅延ヘテロダイン干渉計４０において、第１光路及び第２光路の両者に光周波数シフタ部を備える点と、電気信号生成部８１の構成が、第１測定装置と異なっている。その他の構成は、第１測定装置と同様なので、重複する説明を省略する。
【００９３】
第１測定装置は、光周波数シフタ部や局発電気信号源が１つであるため、第２測定装置に比べて、製造コストの面で有利である。一方、第２測定装置は、合波部で合波される２つの光の周波数が近い値であるため、ホモダイン検波を行うという観点では、より高精度の測定を行うことができる。
【符号の説明】
【００９４】
１０    光源部
２０    サーキュレータ
３０    光増幅器
３２    光バンドパスフィルタ
４０、４１    自己遅延ヘテロダイン干渉計
４２    分岐部
４３    光周波数シフタ部
４４    第１光周波数シフタ部
４６    第２光周波数シフタ部
４８    遅延部
５０    合波部
６０    コヒーレント検波部
６２    バランス型ＰＤ
６４    ＦＥＴ増幅器
７０、８６    ミキサー部
７１    周波数シフト量取得部
７２、８８、９４    ローパスフィルタ（ＬＰＦ）
７４、７５    信号処理装置
８０、８１    電気信号生成部
８２    第１局発電気信号源
８３    局発電気信号源
８４    第２局発電気信号源
９０    タイミング制御器
９１    信号強度取得部
９２    ２乗回路
９６    １／２乗回路
【図１】
0006358277.tif 000014
【図２】
0006358277.tif 000015
【図３】
0006358277.tif 000016
【図４】
0006358277.tif 000017
【図５】
0006358277.tif 000018
【図６】
0006358277.tif 000019
【図７】
0006358277.tif 000020
【図８】
0006358277.tif 000021
【図９】
0006358277.tif 000022
【図１０】
0006358277.tif 000023
【図１１】
0006358277.tif 000024
【図１２】
0006358277.tif 000025