(57)【要約】【課題】短時間で距離の測定ができるテラヘルツレーダ装置を実現すること。【解決手段】モード同期周波数制御手段１２によってフェムト秒レーザ１０、１１のモード同期周波数を安定させ、モード同期周波数の差を一定とする。フェムト秒レーザ１０からのポンプ光を光伝導アンテナ１３に入射させ、テラヘルツコムを発生させる。一方、フェムト秒レーザ１１からのプローブ光は、光伝導アンテナ１４に入射させ、光伝導アンテナ１４に測定物Ｘにより反射されたテラヘルツコムを入射させることにより、ＲＦコムを発生させる。このＲＦコムのコムモードを位相差計１７により抽出し、抽出したコムモードの位相を測定することにより、測定物Ｘまでの距離を算出する。複数のコムモードを選択、抽出することで、高ダイナミックレンジな測距が可能となる。【選択図】図１
2011080860.tif 000002
【特許請求の範囲】
【請求項１】
モード同期周波数が異なる２台のフェムト秒レーザと、
２台の前記フェムト秒レーザのそれぞれのモード同期周波数が安定化され、かつモード同期周波数の差が一定値に保持されるよう制御するモード同期周波数制御手段と、
一方のフェムト秒レーザの出力光をポンプ光として用いてテラヘルツコムを発生させるテラヘルツ波発生手段と、
前記テラヘルツ波発生手段により発生させたテラヘルツコムを任意の方向に走査して測定物に照射するテラヘルツ波照射手段と、
他方のフェムト秒レーザの出力光をプローブ光として用い、前記測定物により反射されたテラヘルツコムと、プローブ光とを入射し、多周波ヘテロダイン光伝導検出により、プローブ光のモード同期周波数とテラヘルツコムのモード同期周波数との差周波数をコム間隔とする周波数コムである電気信号を生成して検出するテラヘルツ波検出手段と、
前記電気信号から単数または複数のコムモードを抽出し、各コムモードで位相を測定する位相測定手段と、
を備えていることを特徴とするテラヘルツレーダ装置。
【請求項２】
前記位相測定手段は、位相の変化量により前記測定物の移動量を測定する手段である、ことを特徴とする請求項１に記載のテラヘルツレーダ装置。
【請求項３】
前記テラヘルツ波検出手段から出力される前記電気信号の周波数スペクトルを測定する周波数スペクトル測定手段と、
測定した周波数スペクトルを、テラヘルツ吸収または散乱の指紋スペクトルのデータベースと比較して前記測定物を識別する測定物識別手段と、
をさらに備えていることを特徴とする請求項１または請求項２に記載のテラヘルツレーダ装置。
【請求項４】
前記位相測定手段は、ロックインアンプを有し、前記ロックインアンプを用いて参照信号に同期させることで前記電気信号からコムモードを抽出する、
ことを特徴とする請求項１ないし請求項３に記載のテラヘルツレーダ装置。
【請求項５】
前記モード同期周波数制御手段は、ポンプ光およびプローブ光からモード同期周波数の差周波数のｎ倍周波数の信号を生成して出力する手段を有し、この信号を前記ロックインアンプの前記参照信号とする、
ことを特徴とする請求項４に記載のテラヘルツレーダ装置。
【請求項６】
前記モード同期周波数制御手段は、前記参照信号の周波数を掃引可能であることを特徴とする請求項５に記載のテラヘルツレーダ装置。
【請求項７】
モード同期周波数が異なる２台のフェムト秒レーザと、
２台の前記フェムト秒レーザのそれぞれのモード同期周波数が安定化され、かつモード同期周波数の差が一定値に保持されるよう制御するモード同期周波数制御手段と、
一方のフェムト秒レーザの出力光をポンプ光として用いてテラヘルツパルス波を発生させるテラヘルツ波発生手段と、
前記テラヘルツ波発生手段により発生させたテラヘルツパルス波を任意の方向に走査して測定物に照射するテラヘルツ波放射手段と、
他方のフェムト秒レーザの出力光をプローブ光として用い、前記測定物により反射されたテラヘルツパルス波と、プローブ光とを入射し、非同期光サンプリング法により、テラヘルツパルス波とプローブ光が同じタイミングで入射したときのみ、テラヘルツパルス波の電場に比例した電気信号を生成して検出するテラヘルツ波検出手段と、
前記ポンプ光と前記プローブ光とから、時間原点信号であるＳＦＧ相互相関信号を生成する時間原点信号生成手段と、
前記電気信号を増幅して前記電気信号の時間波形を、前記時間原点信号をトリガーとして高速デジタイザーで測定し、前記時間原点信号を基準とした前記電気信号の遅れから、前記測定物までの距離を算出する距離算出手段と、
を備えていることを特徴とするテラヘルツレーダ装置。
【請求項８】
前記モード同期周波数制御手段は、モード同期周波数の差を所望の値に設定することができる、ことを特徴とする請求項７に記載のテラヘルツレーダ装置。
【請求項９】
前記モード同期周波数制御手段は、出力光をポンプ光として用いる方の前記フェムト秒レーザのモード同期周波数を掃引可能可能である、
ことを特徴とする請求項７または請求項８に記載のテラヘルツレーダ装置。
【請求項１０】
前記距離算出手段は、ポンプ光のモード同期周波数を変化させた場合の、前記電気信号の遅れの変化量から前記測定物までの距離を算出する手段である、
ことを特徴とする請求項９に記載のテラヘルツレーダ装置。
【請求項１１】
前記テラヘルツ波生成手段は、光伝導アンテナまたは非線形光学結晶であることを特徴とする請求項１ないし請求項１０のいずれか１項に記載のテラヘルツレーダ装置。
【請求項１２】
前記テラヘルツ波検出手段は、光伝導アンテナまたは電気光学結晶であることを特徴とする請求項１ないし請求項１１のいずれか１項に記載のテラヘルツレーダ装置。
【請求項１３】
前記モード同期周波数制御手段は、原子時計を基準として２台の前記フェムト秒レーザのそれぞれのモード同期周波数を制御する、ことを特徴とする請求項１ないし請求項１２のいずれか１項に記載のテラヘルツレーダ装置。
【請求項１４】
前記測定物の距離画像を生成し、その距離画像から前記測定物の識別を行う手段をさらに備える、ことを特徴とする請求項１ないし請求項１３のいずれか１項に記載のテラヘルツレーダ装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、テラヘルツ波を測定物に照射し、測定物によって反射されたテラヘルツ波から距離測定や測定物の識別を行うテラヘルツレーダ装置に関する。
【背景技術】
【０００２】
近年、光と電波の境界の帯域であるテラヘルツ波（周波数０．１〜１０ＴＨｚ）が注目されている。このテラヘルツ領域は、これまで良質なレーザ光源や高感度検出器の開発が遅れていたため、ほとんど研究が行われていなかったが、安定な超短パルスレーザの出現、超高速光技術の発達により、パルス状のテラヘルツ波が比較的容易に生成・検出できるようになり、センシングやイメージングなどの応用開発が加速している。特に、生体分子、ビタミン、糖、医薬品、農薬など、さまざまな物質がテラヘルツ領域において固有の吸収スペクトル（指紋スペクトル）を有することが明らかとなっており、テラヘルツ分光法が新しいセンシング技術として注目されている。
【０００３】
特許文献１には、テラヘルツパルス波を用いたＴＯＦ（タイム・オブ・フライト）法により、塗装膜の厚さを測定する方法が示されている。また、テラヘルツパルス波の時間波形は、プローブ光を機械的ステージを用いて時間遅延させることにより、テラヘルツパルス波とプローブ光の重なるタイミングを少しずつずらしていき、この重なるタイミングを測定することによりテラヘルツパルス波の時間波形をサンプリング測定して求めることが示されている。
【０００４】
また、特許文献２には、非同期光サンプリング法（ＡＯＳ法）と呼ばれる技術が記載されている。従来は、特許文献１に記載のように、機械的ステージによるプローブ光の時間遅延走査によりテラヘルツパルス波の時間波形を求めていたが、機械的ステージによる走査は時間がかかるため、テラヘルツパルス波の時間波形を高速に測定することはできなかった。しかし、ＡＯＳ法を用いると、機械的ステージを用いることなくテラヘルツパルス波の時間波形を高速に測定することができる。
【０００５】
このＡＯＳ法の原理について、図８を参照に説明する。モード同期周波数が僅かに異なる２台のフェムト秒レーザを用い、一方はポンプ光（モード同期周波数ｆ１）、他方はプローブ光（モード同期周波数ｆ２）として用いる。ポンプ光は光伝導アンテナなどに入射させてテラヘルツパルス波を発生させる。テラヘルツパルス波と、プローブ光の重なるタイミングは、モード同期周波数が異なるため、パルスごとにＳ（Ｓ＝Δｆ／（ｆ１＊ｆ２）、Δｆ＝ｆ１−ｆ２）ずつずれていく。このテラヘルツパルス波とプローブ光の重なるタイミングを測定することは、テラヘルツパルス波の時間波形を時間的に拡大した波形をサンプリング測定することと等価である。この時間的に拡大されたテラヘルツパルス波の時間波形は、オシロスコープで直接リアルタイムに観測することができる。
【０００６】
また、非特許文献１には、多周波ヘテロダイン光伝導検出法によりテラヘルツコムをＲＦ帯の周波数コムにダウンスケーリングする技術が示されている。
【０００７】
多周波ヘテロダイン光伝導検出法の原理について、図９を参照に説明する。モード同期周波数が僅かに異なる２台のフェムト秒レーザを用い、一方はポンプ光（モード同期周波数ｆ１）、他方はプローブ光（モード同期周波数ｆ２）として用いる。ポンプ光およびプローブ光は、周波数スペクトル領域で見ると、モード同期周波数の間隔で規則的に多数の周波数モード列が並んだ離散スペクトル構造を有しており、周波数コムと呼ばれている。ポンプ光は光伝導アンテナに入射させ、テラヘルツ領域に展開された周波数コムであるテラヘルツコムを放射させる。
【０００８】
一方、プローブ光を検出用光伝導アンテナに入射させると、超短パルス光による繰り返し超高速スイッチングの結果、検出用光伝導アンテナ内に光励起電流の周波数コム（ＰＣコム）が生成される。このＰＣコムは、モード同期周波数がｆ２で、周波数帯域はテラヘルツコムと同じである。このＰＣコムが励起された検出用光伝導アンテナにテラヘルツコムが入射すると、両者の相互作用により、ビート周波数Δｆ（Δｆ＝ｆ１−ｆ２）を周波数間隔とする２次的な周波数コムがＲＦ帯域に電気信号として発生する。このＲＦ帯域の周波数コム（ＲＦコム）は、テラヘルツコムの周波数スケールをΔｆ／ｆ１でスケーリングしたものである。したがって、このＲＦコムをスペクトラムアナライザで直接観測し、周波数軸をリスケーリングすれば、テラヘルツコムの周波数スペクトルを得ることができる。従来は、機械的ステージによるプローブ光の時間遅延走査をしてテラヘルツ波の時間波形を測定し、ＦＦＴ処理することで周波数スペクトルを測定していたが、多周波ヘテロダイン光伝導検出法を用いると、直接テラヘルツ波の周波数スペクトルを得られるので、高速に測定することができる。
【先行技術文献】
【特許文献】
【０００９】
【特許文献１】特開２００４−２８６１８
【特許文献２】ＷＯ２００６／０９２８７４
【非特許文献】
【００１０】
【非特許文献１】T. Yasui, Y. Kabetani, E. Saneyosh, S. Yokoyama, and T. Araki, "Terahertz frequency comb by multi-frequency-heterodyning photoconductive detection for high-accuracy, high-resolution terahertz spectroscopy", Appl. Phys. Lett., Vol. 88, 241104 (2006)
【発明の概要】
【発明が解決しようとする課題】
【００１１】
しかし、特許文献１に記載のテラヘルツパルス波を用いたＴＯＦ法による膜厚測定では、機械的ステージによる時間遅延走査のためテラヘルツ波の時間波形を測定するのに時間がかかり、その結果膜厚測定にも時間がかかった。
【００１２】
また、特許文献２や非特許文献１に記載のＡＯＳ法、多周波ヘテロダイン光伝導検出法によって測定物までの距離を測定するテラヘルツレーダ装置については、知られていない。
【００１３】
そこで本発明の目的は、ＡＯＳ法や多周波ヘテロダイン光伝導検出法によって測定物までの距離測定を短時間で行うことができるテラヘルツレーダ装置を提供することである。
【課題を解決するための手段】
【００１４】
第１の発明は、モード同期周波数が異なる２台のフェムト秒レーザと、２台のフェムト秒レーザのそれぞれのモード同期周波数が安定化され、かつモード同期周波数の差が一定値に保持されるよう制御するモード同期周波数制御手段と、一方のフェムト秒レーザの出力光をポンプ光として用いてテラヘルツコムを発生させるテラヘルツ波発生手段と、テラヘルツ波発生手段により発生させたテラヘルツコムを任意の方向に走査して測定物に照射するテラヘルツ波照射手段と、他方のフェムト秒レーザの出力光をプローブ光として用い、測定物により反射されたテラヘルツコムと、プローブ光とを入射し、多周波ヘテロダイン光伝導検出により、プローブ光のモード同期周波数とテラヘルツコムのモード同期周波数との差周波数をコム間隔とする周波数コムである電気信号を生成して検出するテラヘルツ波検出手段と、電気信号から単数または複数のコムモードを抽出し、各コムモードで位相を測定する位相測定手段と、を備えていることを特徴とするテラヘルツレーダ装置である。
【００１５】
第２の発明は、第１の発明において、位相測定手段は、位相の変化量により測定物の移動量を測定する手段である、ことを特徴とするテラヘルツレーダ装置である。
【００１６】
第３の発明は、第１の発明から第２の発明において、テラヘルツ波検出手段から出力される電気信号の周波数スペクトルを測定する周波数スペクトル測定手段と、測定した周波数スペクトルを、テラヘルツ吸収または散乱の指紋スペクトルのデータベースと比較して測定物を識別する測定物識別手段と、をさらに備えていることを特徴とするテラヘルツレーダ装置である。
【００１７】
第４の発明は、第１の発明から第３の発明において、位相測定手段はロックインアンプを有し、ロックインアンプを用いて参照信号に同期させることで電気信号からコムモードを抽出することを特徴とするテラヘルツレーダ装置である。
【００１８】
第５の発明は、第４の発明において、モード同期周波数制御手段は、ポンプ光およびプローブ光からモード同期周波数の差周波数のｎ倍周波数の信号を生成して出力する手段を有し、この信号をロックインアンプの参照信号とする、ことを特徴とするテラヘルツレーダ装置である。
【００１９】
第６の発明は、第５の発明において、モード同期周波数制御手段は、参照信号の周波数を掃引可能であることを特徴とするテラヘルツレーダ装置である。
【００２０】
第７の発明は、モード同期周波数が異なる２台のフェムト秒レーザと、２台のフェムト秒レーザのそれぞれのモード同期周波数が安定化され、かつモード同期周波数の差が一定値に保持されるよう制御するモード同期周波数制御手段と、一方のフェムト秒レーザの出力光をポンプ光として用いてテラヘルツコムを発生させるテラヘルツ波発生手段と、テラヘルツ波発生手段により発生させたテラヘルツコムを任意の方向に走査して測定物に照射するテラヘルツ波照射手段と、他方のフェムト秒レーザの出力光をプローブ光として用い、測定物により反射されたテラヘルツパルス波と、プローブ光とを入射し、非同期光サンプリング法により、テラヘルツパルス波とプローブ光が同じタイミングで入射したときのみ、テラヘルツパルス波の電場に比例した電気信号を生成して検出するテラヘルツ波検出手段と、ポンプ光とプローブ光とから時間原点信号であるＳＦＧ（和周波混合）相互相関信号を生成する時間原点信号生成手段と、電気信号を増幅して電気信号の時間波形を、時間原点信号をトリガーとして高速デジタイザーで測定し、時間原点信号を基準とした電気信号の遅れから、測定物までの距離を算出する距離算出手段と、を備えていることを特徴とするテラヘルツレーダ装置である。
【００２１】
第８の発明は、第７の発明において、モード同期周波数制御手段は、モード同期周波数の差を所望の値に設定することができる、ことを特徴とするテラヘルツレーダ装置である。モード同期周波数の差は、２台のフェムト秒レーザの一方または双方のモード同期周波数を変えることで所望の値に設定できる。
【００２２】
第９の発明は、第７の発明または第８の発明において、モード同期周波数制御手段は、出力光をポンプ光として用いる方のフェムト秒レーザのモード同期周波数を掃引可能可能である、ことを特徴とするテラヘルツレーダ装置である。
【００２３】
第１０の発明は、第９の発明において、距離算出手段は、ポンプ光のモード同期周波数を変化させた場合の、電気信号の遅れの変化量から測定物までの距離を算出する手段であることを特徴とするテラヘルツレーダ装置である。
【００２４】
第１１の発明は、第１の発明から第１０の発明において、テラヘルツ波生成手段は、光伝導アンテナまたは非線形光学結晶であることを特徴とするテラヘルツレーダ装置である。
【００２５】
第１２の発明は、第１の発明から第１１の発明において、テラヘルツ波検出手段は、光伝導アンテナまたは電気光学結晶であることを特徴とするテラヘルツレーダ装置である。
【００２６】
第１３の発明は、第１の発明から第１２の発明において、モード同期周波数制御手段は、原子時計を基準として２台のフェムト秒レーザのそれぞれのモード同期周波数を制御する、ことを特徴とするテラヘルツレーダ装置である。
【００２７】
第１４の発明は、第１の発明から第１３の発明において、測定物の距離画像を生成し、その距離画像から測定物の識別を行う手段をさらに備える、ことを特徴とするテラヘルツレーダ装置である。
【発明の効果】
【００２８】
第１の発明によると、多周波ヘテロダイン光伝導検出法を用いることにより、測定物までの距離を位相差によって短時間で測定することができる。また、複数のコムモードを抽出すれば、それぞれのコムモードにおいて位相差を求めることができ、高分解能で高ダイナミックレンジな測距が可能である。また、テラヘルツ波をレーダー波として利用するので、天候の影響を受けることなく測距が可能である。
【００２９】
また、第２の発明のように、測定物が移動体である場合にも、位相の変化量を高速に求めることができるので、測定物の移動量をリアルタイムで計測することができる。
【００３０】
また、第３の発明のように、多周波ヘテロダイン光伝導検出法を用いているため周波数スペクトルを高速に測定することができ、距離の測定と同時に測定物の識別をすることができる。
【００３１】
また、第４の発明のように、コムモードの抽出はロックインアンプを用いて行うことができ、第５の発明のようにロックインアンプの参照信号はポンプ光およびプローブ光から生成することができる。また、第６の発明のように、ロックインアンプの参照信号を掃引することで複数のコムモードを抽出することができる。
【００３２】
また、第７の発明によると、ＡＯＳ法により短時間でテラヘルツパルス波の時間波形を得ることができるため、測定物までの距離を短時間で測定することができる。したがって、測定物が移動体である場合にもリアルタイムで測定物の移動量を計測することができる。また、ＡＯＳ法によりテラヘルツ波の時間波形を時間的に拡大できるので、検出器の応答速度によらず、フェムト秒レーザの出力光のパルス幅のみに依存する高分解測距が可能である。また、テラヘルツ波をレーダー波として利用するので、天候の影響を受けることなく測距が可能である。
【００３３】
また、第８の発明のように、モード同期周波数の差を変えることで、テラヘルツパルス波の時間波形計測におけるサンプリング間隔や、テラヘルツパルス波の時間拡大率などを変更することができる。また、第９の発明のように、ポンプ光のモード同期周波数を変更することで測距レンジを変更することができる。
【００３４】
また、第１０の発明によると、測定物が測距レンジを超える距離に位置し、電気信号の遅れからは相対的な距離（位相差）しか測定できない場合であっても、ポンプ光のモード同期周波数を変化させ、異なるモード同期周波数のテラヘルツパルス波でそれぞれ位相差を測定し、その位相差の変化量を求めれば、測定物までの距離を測定することができる。
【００３５】
また、第１１の発明のように、テラヘルツ波の発生手段として光伝導アンテナまたは非線形光学結晶を用いることができ、第１２の発明のように、テラヘルツ波検出手段として光伝導アンテナまたは電気光学結晶を用いることができる。
【００３６】
また、第１３の発明のように、原子時計を基準としてモード同期周波数を制御することで、フェムト秒レーザのモード同期周波数を精度よく制御することができ、モード同期周波数の差を安定化することができる。
【００３７】
また、第１４の発明のように、本発明のテラヘルツレーダ装置によれば測定物へのテラヘルツ波の照射位置を操作することができるので、測定物の距離画像を得ることができ、その距離画像から測定物を識別することが可能となる。
【図面の簡単な説明】
【００３８】
【図１】実施例１のテラヘルツレーダ装置の構成を示した図。
【図２】モード同期周波数制御手段１２のより詳細な構成を示した図。
【図３】測定した位相の時間変化を示した図。
【図４】実施例２のテラヘルツレーダ装置の構成を示した図。
【図５】測定した変位量と真の変位量との関係を示したグラフ。
【図６】テラヘルツパルス波の時間波形を測定した結果を示す図。
【図７】得られたテラヘルツパルス波の時間波形を重ねて示した図。
【図８】ＡＯＳ法の原理について説明する図。
【図９】多周波ヘテロダイン光伝導検出法の原理について説明する図。
【発明を実施するための形態】
【００３９】
以下、本発明の具体的な実施例について図を参照に説明するが、本発明は実施例に限定されるものではない。
【実施例１】
【００４０】
実施例１のテラヘルツレーダ装置は、テラヘルツ波を用い、位相差から距離を測定する装置であり、同時に測定物の識別を行う装置である。図１は、実施例１のテラヘルツレーダ装置の構成を示す図である。テラヘルツレーダ装置は、２台のフェムト秒レーザ１０、１１と、２台のフェムト秒レーザ１０、１１のモード同期周波数を制御するモード同期周波数制御手段１２と、光伝導アンテナ１３、１４と、放射方向制御手段１５と、電流電圧変換アンプ１６と、位相差計１７と、スペクトラムアナライザ１９と、によって構成されている。
【００４１】
フェムト秒レーザ１０、１１は、モード同期チタンサファイアレーザーであり、中心波長８００ｎｍ、パルス幅１０ｆｓのパルス光を出力する。モード同期周波数（パルスの繰り返し周波数）は、ピエゾ素子による共振器長の制御により可変である。フェムト秒レーザ１０は、モード同期周波数ｆ１であり、その出力光はポンプ光として用いる。フェムト秒レーザ１１は、モード同期周波数ｆ２であり、フェムト秒レーザ１０のモード同期周波数ｆ１とは異なる値である。フェムト秒レーザ１１の出力光はプローブ光として用いる。フェムト秒レーザ１０、１１の出力光は、周波数スペクトル領域ではモード同期周波数の間隔で規則的に多数の周波数モード列（コムモード）が並んだ離散スペクトル構造を有した周波数コムである。
【００４２】
モード同期周波数制御手段１２は、フェムト秒レーザ１０、１１のモード同期周波数ｆ１、ｆ２を安定化させ、モード同期周波数の差Δｆ（Δｆ＝ｆ１−ｆ２）を一定に保持するよう制御する手段であり、また位相差計１７に入力する参照信号を生成する手段である。
【００４３】
図２は、モード同期周波数制御手段１２のより詳細な構成を示す図である。モード同期周波数制御手段１２は、ルビジウム原子時計３０と、周波数シンセサイザ３１と、発振器３２、３３と、超高速フォトディテクタ３４、３５と、ミキサ３６、３７、３８と、ローパスフィルタ３９、４０と、ＰＬＬ回路４１、４２と、ピエゾ素子コントローラ４３、４４と、により構成されている。周波数シンセサイザ３１、発振器３２、３３の出力信号は、ルビジウム原子時計３０の出力信号を基準として制御されている。
【００４４】
モード同期周波数制御手段１２は、以下のような動作により、フェムト秒レーザ１０、１１のモード同期周波数ｆ１、ｆ２を安定化し、Δｆを一定に保持する。
【００４５】
ビームスプリッタ４５により抽出されたフェムト秒レーザ１０の出力光は、超高速フォトディテクタ３４に検出される。検出される信号のうち、１００次の高調波成分である、１００ｆ１を抽出するため、周波数シンセサイザ３１の出力（周波数ｆｓ）をｆ１付近とし、超高速フォトディテクタ３４の出力と周波数シンセサイザ３１の出力をミキサ３６に入力し、ローパスフィルタ３９を通すことで、周波数１００ｆ１−ｆｓの信号を抽出する。そして、ＰＬＬ回路４１において、発振器３２から出力された基準信号を分周したものと、周波数１００ｆ１−ｆｓの信号とを比較して、その誤差信号をピエゾ素子コントローラ４３に入力してフィードバック制御する。１００次の高調波成分１００ｆ１を用い、モード同期周波数ｆ１の揺らぎを拡大して制御しているため、また、周波数シンセサイザ３１、発振器３２の出力信号はルビジウム原子時計３０の出力信号を基準とすることで高精度に安定化されているため、ｆ１の安定化を高精度に行うことができる。
【００４６】
フェムト秒レーザ１１のモード同期周波数ｆ２も、同様にして制御され、安定化される。すなわち、ビームスプリッタ４６によりフェムト秒レーザ１０の出力光を抽出して超高速フォトディテクタ３５により検出し、ミキサ３７とローパスフィルタ４０により周波数１００ｆ２−ｆｓの信号を抽出し、ＰＬＬ回路４２において、発振器３３から出力された基準信号を分周したものと、周波数１００ｆ２−ｆｓの信号とを比較して、その誤差信号をピエゾ素子コントローラ４４に入力してフィードバック制御することにより、ｆ２の安定化を高精度に行う。
【００４７】
また、発振器３２と発振器３３の周波数差がｎΔｆとなるようにし、発振器３２、３３の出力信号をミキサ３８に入力して周波数がｎΔｆの信号を生成する。この周波数ｎΔｆの信号は、位相差計１７の参照信号とする。ｎの値は、発振器３２、３３の出力信号の周波数を調整することで任意の値を選択することができる。発振器３２、３３の出力信号はルビジウム原子時計３０により高精度に安定化されているため、この参照信号もまた、高精度に安定化されている。
【００４８】
光伝導アンテナ１３は、光伝導膜上にボウタイ形状の一対の電極が微小なギャップを隔てて形成された構造である。光伝導アンテナ１３の電極にバイアスを印加して電極間のギャップにフェムト秒レーザ１０からのポンプ光を入射させると、一種の波長変換が生じ、パルス状のテラヘルツ波が発生する。このテラヘルツ波は、周波数スペクトル領域では、基本波成分ｆ１と多数の高調波成分（２ｆ１、３ｆ１、４ｆ１、・・・）が等間隔で並び、テラヘルツ領域にまで展開された周波数コムであり、以下テラヘルツコムと呼ぶ。
【００４９】
光伝導アンテナ１４は、測定物により反射されたテラヘルツコムを、多周波ヘテロダイン光伝導検出法により検出するための装置であり、構造は光伝導アンテナ１３と同様のボウタイ型である。多周波ヘテロダイン光伝導検出法は、テラヘルツコムを電波周波数（ＲＦ）帯まで正確にダウンスケーリングする方法である。光伝導アンテナ１４には、フェムト秒レーザ１１からのプローブ光が入力され、光伝導アンテナ１４内に光励起電流の周波数コム（ＰＣコム）が生成される。そして、このＰＣコムが励起された光伝導アンテナ１４にテラヘルツコムを入射させ、ＰＣコムとテラヘルツコムとを相互作用させることによって、ビート周波数Δｆ（Δｆ＝ｆ１−ｆ２）を周波数間隔とする２次的な周波数コムがＲＦ帯域に電気信号として発生する。このＲＦ帯域の周波数コム（ＲＦコム）は、テラヘルツコムの周波数スケールをΔｆ／ｆ１でスケーリングしたものである。
【００５０】
放射方向制御手段１５は、ミラーやビームスプリッタによって構成され、光伝導アンテナ１３によって発生したテラヘルツコムを任意の方向に走査して測定物Ｘに照射し、測定物Ｘによって反射されたテラヘルツコムを光伝導アンテナ１４に導く手段である。
【００５１】
電流電圧変換アンプ１６は、光伝導アンテナ１４から出力されるＲＦコムを増幅して、位相差計１７とスペクトラムアナライザ１９に出力する。位相差計１７はロックインアンプであり、モード同期周波数制御手段１２からの参照信号に同期する周波数のコムモードをＲＦコムから抽出し、そのコムモードの信号の位相差を測定する。
【００５２】
スペクトラムアナライザ１９は、ＲＦコムの周波数スペクトルを測定し、その周波数スペクトルをもとの周波数スケールに戻し、データベース２１に記憶された各物質が有するテラヘルツ帯域における固有の吸収または散乱スペクトル（指紋スペクトル）と比較することで、測定物Ｘの識別を行う。
【００５３】
なお、テラヘルツレーダ装置内におけるポンプ光、プローブ光、およびテラヘルツコムの伝搬は、空間を伝搬させてもよいし、光ファイバや光導波路を用いて伝搬させてもよい。
【００５４】
次に、実施例１のテラヘルツレーダ装置による測定物Ｘまでの距離測定方法、および測定物Ｘの識別方法について説明する。
【００５５】
まず、モード同期周波数制御手段１２によってフェムト秒レーザ１０、１１のモード同期周波数ｆ１、ｆ２を安定させ、モード同期周波数の差Δｆを一定とする。そして、フェムト秒レーザ１０からのポンプ光を光伝導アンテナ１３に入射させ、テラヘルツコムを発生させる。発生したテラヘルツコムは、放射方向制御手段１５によって導かれ、テラヘルツコムの放射方向を制御して走査し、測定物Ｘにテラヘルツコムを照射する。測定物Ｘにより反射されたテラヘルツコムは、放射方向制御手段１５に導かれて光伝導アンテナ１４に入射する。
【００５６】
一方、フェムト秒レーザ１１からのプローブ光は、光伝導アンテナ１４に入射させ、光伝導アンテナ１４内にＰＣコムを励起させる。このＰＣコムが励起された状態で、光伝導アンテナ１４に測定物Ｘにより反射されたテラヘルツコムを入射させることにより、ＲＦコムを発生させる。ＲＦコムは、電流電圧変換アンプ１６によって増幅された後、位相差計１７とスペクトラムアナライザ１９に入力される。
【００５７】
位相差計１７では、参照信号に同期させることでＲＦコムの特定のコムモードを抽出する。参照信号はモード同期周波数制御手段１２により生成された周波数ｎΔｆの信号を用いる。これにより、周波数がｎΔｆのコムモードを抽出することができる。そして、抽出したコムモードの位相を測定することにより、測定物Ｘまでの相対的な距離、たとえば測定物Ｘの移動量や、測定物Ｘ表面の凹凸形状などを求めることができる。また、参照信号の周波数は可変であり、複数のコムモードを選択、抽出することができるので、たとえば、周波数のオーダーが異なる複数のコムモードを抽出して位相を測定すれば、高ダイナミックレンジな測距が可能となる。また、異なる周波数のコムモード（つまり次数の異なるコムモード、たとえば１６００次と１６０１次のコムモード）を抽出し、それぞれの位相を求めてその差を算出すれば、絶対距離が原理的に測定可能である。
【００５８】
他方、スペクトラムアナライザ１９では、ＲＦコムの周波数スペクトルを測定する。この得られた周波数スペクトルの周波数スケールをｆ１／Δｆ倍することでリスケーリングし、データベース２１に記憶された指紋スペクトルと比較することで、測定物Ｘの材質を識別する。
【００５９】
以上のように、実施例１のテラヘルツレーダ装置によると、さまざまな波長での位相測定により高分解で高ダイナミックレンジな測距が高速で可能であり、また、測距と同時に測定物の識別も行うことができる。また、テラヘルツ波を用いるため、天候の影響を受けることなく高精度な測距を行うことが可能である。
【００６０】
図３は、ｆ１＝８１．８ＭＨｚ、Δｆ＝１０Ｈｚとし、１６００次のコムモード（テラヘルツコムでは１３１ＧＨｚ（波長２．２９ｍｍ）、ＲＦコムでは１６ｋＨｚ）を選択し、測定物Ｘまでの位相差（距離）を１秒ごとに測定した結果である。図３のグラフにおいて横軸は測定を開始してからの経過時間（単位は秒）であり、縦軸は測定開始時の位相を基準とした位相差（単位は度）である。測定物Ｘは５０秒毎に５００μｍずつ階段状に移動させた。
【００６１】
図３を見ると、値に揺らぎはあるが、約５０秒ごとに位相が約７７°変化していることがわかる。波長２．２９ｍｍであるから７７°の位相変化量は約５００μｍであって、測定物Ｘの移動量とおおよそ一致していることから、確かに実施例１のテラヘルツレーダ装置により測距が可能であることがわかる。
【００６２】
次に、図３での条件と同じく１６００次のコムモードを選択して求めた位相変化量と、測定物Ｘの移動量との線形性を評価したところ、次の結果を得た。０．５ｍｍ、１ｍｍ、１．５ｍｍ、２ｍｍの測定物Ｘの移動に伴う位相変化量を各５回測定し、位相変化を光路長と比べて評価した。位相の変動する時間を無視し、その前後での平均値の変化を位相変化量とした。そして、５回測定した位相変化量の平均値と標準偏差を計算した。測定物Ｘの移動量の真値と測定値との誤差の標準偏差を測距精度とすると、１００μｍの測距精度が達成されていることがわかった。
【実施例２】
【００６３】
実施例２のテラヘルツレーダ装置は、テラヘルツパルス波を用い、ＴＯＦ（タイム・オブ・フライト）法により測定物Ｘまでの距離を測定する装置である。図４は、実施例２のテラヘルツレーダ装置の構成を示す図である。実施例２のテラヘルツレーダ装置は、２台のフェムト秒レーザ１１０、１１１と、２台のフェムト秒レーザ１１０、１１１のモード同期周波数を制御するモード同期周波数制御手段１１２と、光伝導アンテナ１１３、１１４と、放射方向制御手段１１５と、電流電圧変換アンプ１１６と、ＳＦＧ強度相互相関計１２２と、高速デジタイザ１２３と、コンピュータ１２４と、レンズ１２５、１２６によって構成されている。
【００６４】
フェムト秒レーザ１１０、１１１は、実施例１のテラヘルツレーダ装置に用いたフェムト秒レーザ１０、１１と同様である。フェムト秒レーザ１１０の出力光はポンプ光として用い、フェムト秒レーザ１１１の出力光はプローブ光として用いる。ポンプ光のモード同期周波数ｆ１と、プローブ光のモード同期周波数ｆ２は、僅かに異なる値である。
【００６５】
モード同期周波数制御手段１１２は、図２に示すモード同期周波数制御手段１２から参照信号を生成するミキサを省略した構成である。モード同期周波数制御手段１１２により、ポンプ光のモード同期周波数ｆ１、およびプローブ光のモード同期周波数ｆ２は高精度に安定化されている。フェムト秒レーザ１０のモード同期周波数ｆ１は掃引可能であり、モード同期周波数の差Δｆ（＝ｆ１−ｆ２）は所望の値に設定することができる。
【００６６】
光伝導アンテナ１１３は、櫛形電極構造を持つ大口径の光伝導アンテナである。フェムト秒レーザ１１０からのポンプ光をレンズ１２５により集光して光伝導アンテナ１１３に入射させることで、高強度で広帯域のテラヘルツパルス波を発生させる。このテラヘルツパルス波は、モード同期周波数ｆ１で、パルス幅はｐｓオーダーである。
【００６７】
光伝導アンテナ１１４は、光伝導アンテナ１３、１４と同様のボウタイ型の電極構造であり、非同期光サンプリング法（ＡＯＳ法）により時間的に拡大されたテラヘルツパルス波の時間波形を検出する手段である。光伝導アンテナ１１４には、レンズ１２６によって集光されてプローブ光が入射し、測定物Ｘにより反射されて戻ってきたテラヘルツパルス波が入射する。光伝導アンテナ１１４は、プローブ光とテラヘルツパルス波が同じタイミングで入射した場合のみ、テラヘルツパルス波の電場に比例した電気信号を出力する。この電気信号の時間波形は、テラヘルツパルス波の電場時間波形を時間的にｆ１／Δｆ倍に拡大したものである。これは、テラヘルツパルス波のモード同期周波数ｆ１と、プローブ光のモード同期周波数ｆ２が僅かに異なるので、光伝導アンテナ１１４においてテラヘルツパルス波とプローブ光の重なるタイミングが、パルスごとにＳ（Ｓ＝Δｆ／（ｆ１＊ｆ２）、Δｆ＝ｆ１−ｆ２）ずつ自動的にずれていくためであり、その重なるタイミングの計測は、テラヘルツパルス波の時間波形を間隔Ｓでサンプリング計測して時間的に拡大したのと等価である。ずれていったテラヘルツパルス波とプローブ光の重なるタイミングが、再びもとに戻るまでに要する時間は、１／Δｆであり、これはテラヘルツパルス波の電場時間波形を得るのに要する時間（測定周期）である。
【００６８】
放射方向制御手段１１５は、光伝導アンテナ１１３から放射されたテラヘルツパルス波を軸外し放物面鏡により平行光にコリメートして測定物Ｘに照射し、測定物Ｘにより反射されたテラヘルツパルス波を導いて光伝導アンテナ１１４に入射させる手段である。テラヘルツパルス波は、ミラーによって放射方向を制御して走査することで測定物Ｘに照射する。測定物Ｘにより反射されたテラヘルツパルス波は、軸外し放物面鏡により集光して光伝導アンテナ１１４に入射させる。
【００６９】
ＳＦＧ強度相互相関計１２２は、和周波光を発生させるための非線形光学結晶と、その和周波光を検出する光電子増倍管と、光電子増倍管からの出力を増幅する電流電圧変換アンプにより構成されていて、ポンプ光およびプローブ光からＳＦＧ強度相互相関信号を生成する手段である。ビームスプリッタ１２７、１２８によって抽出されたポンプ光とプローブ光は、非線形光学結晶に入力され、和周波混合によりＳＦＧ相互相関信号光を発生させる。これを光電子増倍管により検出し、電流電圧変換アンプにより増幅することで、ＳＦＧ強度相互相関信号を得る。ＳＦＧ強度相互相関信号は、高速デジタイザ１２３のトリガー信号、つまり時間原点信号として用いる。
【００７０】
高速デジタイザ１２３は、光伝導アンテナ１１４から出力され、電流電圧変換アンプ１１６により増幅された電気信号の時間波形を測定する。サンプリング速度は、５００ｋサンプル／秒である。
【００７１】
コンピュータ１２４は、高速デジタイザ１２３によって測定されたテラヘルツパルス波の時間波形からパルスの遅れを求め、パルスの遅れから測定物Ｘまでの距離を算出する手段である。
【００７２】
なお、実施例１のテラヘルツレーダ装置の場合と同様に、テラヘルツレーダ装置内におけるポンプ光、プローブ光、およびテラヘルツパルス波の伝搬は、空間を伝搬させてもよいし、光ファイバや光導波路を用いて伝搬させてもよい。
【００７３】
次に、実施例２のテラヘルツレーダ装置による測定物Ｘまでの距離測定方法について説明する。
【００７４】
まず、モード同期周波数制御手段１１２によってフェムト秒レーザ１１０、１１１のモード同期周波数ｆ１、ｆ２を安定させる。そして、フェムト秒レーザ１１０からのポンプ光をレンズ１２５により集光して光伝導アンテナ１１３に入射させ、テラヘルツパルス波を発生させる。発生したテラヘルツパルス波は、放射方向制御手段１１５によって導かれ、テラヘルツパルス波の放射方向を制御して走査し、測定物Ｘにテラヘルツパルス波を照射する。測定物Ｘにより反射されたテラヘルツパルス波は、放射方向制御手段１１５に導かれて光伝導アンテナ１１４に入射する。
【００７５】
一方、フェムト秒レーザ１１１からのプローブ光を光伝導アンテナ１１４に入射させ、プローブ光とテラヘルツパルス波が同じタイミングで入射した場合のみ、テラヘルツパルス波の電場に比例した電気信号を出力する。この電気信号の時間波形は、上述のように、テラヘルツパルス波の電場時間波形を間隔Ｓ（＝Δｆ／（ｆ１＊ｆ２）でサンプリングして時間的にｆ１／Δｆ倍に拡大したものである。光伝導アンテナ１１４から出力された電気信号は、電流電圧変換アンプ１１６によって増幅された後、高速デジタイザ１２３に入力される。
【００７６】
また、ビームスプリッタ１２７、１２８によってポンプ光、プローブ光を抽出してＳＦＧ強度相互相関計１２２に入射させ、ＳＦＧ強度相互相関信号を発生させる。
【００７７】
次に、高速デジタイザ１２３により、ＳＦＧ強度相互相関信号をトリガー信号（時間原点信号）として電流電圧変換アンプ１１６からの電気信号の時間波形を測定する。そして、コンピュータ１２４において、測定された電気信号の時間波形からパルスの遅れ（トリガー信号からのパルスの遅れ）を求め、パルスの遅れから測定物Ｘまでの距離を算出する。たとえば、パルスの時間遅れがｔであれば、時間スケールをリスケーリングしてｔ＊（Δｆ／ｆ１）とし、光速ｃを乗じた値、ｃｔ＊（Δｆ／ｆ１）が、テラヘルツレーダ装置と測定物Ｘ間の往復の光路長であるから、測定物Ｘまでの距離は、ｃｔ＊（Δｆ／ｆ１）＊（１／２）となる。１回の時間波形取得に要する時間は、測定周期１／Δｆであるから、１／Δｆごとに測定物Ｘの距離測定を行うことができる。
【００７８】
また、テラヘルツパルス波のサンプリング間隔Ｓ、時間拡大率、測定周期は、いずれもモード同期周波数の差Δｆに依存しているので、モード同期周波数制御手段１１２によりｆ１、ｆ２を変えることで、サンプリング間隔Ｓ、時間拡大率、測定周期を所望の値に設定することができる。特に、ｆ１を変更することで、測距レンジｃ／ｆ１を変更することができる。
【００７９】
また、異なるｆ１でパルスの遅れを測定し、パルスの遅れの変化量を測定すれば、測定物Ｘの距離が測距レンジを超えている場合であっても、測定物Ｘまでの距離を測定することができる。
【００８０】
以上のように、実施例２のテラヘルツレーダ装置によると、機械的ステージを用いないため、テラヘルツパルス波の時間波形を高速に求めることができ、テラヘルツパルス波を時間的に拡大して計測するので検出器の応答速度によらず高分解な計測が可能であるため、高分解な測距を高速に行うことができる。したがって、測定物が移動体の場合にもリアルタイムで測距を行うことが可能である。また、測定レンジを超える距離に測定物がある場合であっても、測定物までの距離を測定することができる。また、テラヘルツ波を用いるため、天候の影響を受けることなく高精度な測距を行うことが可能である。
【００８１】
図５は、測定物Ｘの真の変位量と、実施例２のテラヘルツレーダ装置を用いて測定した変位量との関係を示したグラフである。ｆ１は８１．８ＭＨｚ、Δｆは１０Ｈｚとした。ここで変位量は、機械的ステージに乗せた測定物Ｘをテラヘルツレーダ装置から１ｍ遠方に配置し、測定物Ｘを１ｃｍずつ動かした時の、初期位置（テラヘルツレーダ装置から１ｍ遠方の位置）からの測定物Ｘの移動量である。なお、測定物Ｘにはアルミ板を用い、初期位置から２０ｃｍの距離まで測定物Ｘを移動させた。図５より、真の変位量と測定した変位量とに良好な線形関係が得られていることがわかる。また、測定物Ｘの真の変位量と測定した変位量との誤差の標準偏差を測距精度とすると、８００μｍの測距精度が得られていることがわかった。
【００８２】
図６は、測定物Ｘを機械的ステージに乗せて４０ｍｍ／ｓの速度で連続的に移動させながら、テラヘルツパルス波の時間波形を測定した結果である。ｆ１は８１．８ＭＨｚ、Δｆは１０Ｈｚとした。図６は１秒ごとの時間波形を示しており、グラフの横軸は時間をサンプリング点数で表わしたものである。パルス１周期は５００００ｐｏｉｎｔであり、グラフでは１００００ｐｏｉｎｔ分、つまりパルス１／５周期分の時間波形を示している。図６をみると、１秒後ごとにパルスのピークが約１１００ｐｏｉｎｔ移動していることがわかる。１１００ｐｏｉｎｔは距離にして約４０ｍｍであり、測定物Ｘの１秒間の移動量におおよそ一致している。このように、実施例２のテラヘルツレーダ装置によると、測定物Ｘが移動体である場合にも、距離の測定が可能である。
【００８３】
また、図７は、測定物Ｘとして階段状の段差（各段差は５０ｍｍ、５０ｍｍ、７０ｍｍ、１７０ｍｍ）が設けられたサンプルを用意し、テラヘルツパルス波を走査してサンプルの各段に照射して時間波形を測定し、得られた時間波形を重ねて示したグラフである。グラフの横軸はテラヘルツレーダ装置からの距離であり、時間波形のピーク位置は、各段までの距離を示している。ｆ１は８１．８ＭＨｚ、Δｆは１０Ｈｚとした。図７のように、ピークの間隔として５０ｍｍ、７０ｍｍ、１７０ｍｍが得られており、これは確かにサンプルの各段差の距離に一致している。
【００８４】
なお、実施例１、２において、テラヘルツ波の発生手段、検出手段として光伝導アンテナを用いているが、非線形光学結晶を用いてテラヘルツ波を発生、検出してもよい。
【００８５】
また、実施例２のテラヘルツレーダ装置において、テラヘルツパルス波の時間波型を元の時間スケールに戻し、フーリエ変換して周波数スペクトルを測定し、テラヘルツ吸収または散乱の指紋スペクトルのデータベースと比較する手段を設ければ、実施例１のテラヘルツレーダ装置のように、測定物の距離測定と平行して測定物の識別を行うことも可能である。
【００８６】
また、実施例１、２において、テラヘルツ波を放射方向制御手段により２次元的に走査すれば、測定物の２次元的な距離画像を得ることもできる。また、得られた測定物の距離画像を画像データベース等で比較することにより測定物を識別する手段を、実施例１、２のテラヘルツレーダ装置に設けてもよい。
【産業上の利用可能性】
【００８７】
本発明のテラヘルツレーダ装置によると、短時間で測定物の距離を測定することができるので、リアルタイムでの測距が可能である。
【符号の説明】
【００８８】
１０、１１、１１０、１１１：フェムト秒レーザ
１２、１１２：モード同期周波数制御手段
１３、１４、１１３、１１４：光伝導アンテナ
１５、１１５：放射方向制御手段
１６、１１６：電流電圧変換アンプ
１７：位相差計
１９：オシロスコープ
１２２：ＳＦＧ強度相互相関計
１２３：高速デジタイザ
１２４：コンピュータ
【図１】
2011080860.tif 000003
【図２】
2011080860.tif 000004
【図３】
2011080860.tif 000005
【図４】
2011080860.tif 000006
【図５】
2011080860.tif 000007
【図８】
2011080860.tif 000008
【図９】
2011080860.tif 000009
【図６】
2011080860.tif 000010
【図７】
2011080860.tif 000011