(57)【特許請求の範囲】
【請求項１】
複数の送信機および受信機を具備し、送信機が受信機から離れて位置するバイスタティックレーダシステムを用いて海洋の表面流ベクトルをマッピングする方法であって、
前記受信機でエコーを生成するために、前記送信機から散乱領域内の波に向けて送信されたレーダ信号を散乱させるステップと、
前記散乱領域への方位角（φ）を決定するステップと、
前記送信機から前記受信機までのレーダエコー時間遅延を測定するために、送信後の時間に対してサンプリングするステップと、
前記散乱領域の位置を決定するステップと、
前記散乱領域における流速を決定するステップとを具備し、
前記散乱領域の位置を決定するステップは、
一定時間遅延（Ｄ）が前記測定されたレーダエコー時間遅延と等しくなるように、前記散乱領域が通過し、前記送信機および受信機が共焦点となる楕円の長軸（Ａ）を、次式を用いて決定するステップと、
【数８】
0004043878.tif 000018
ここで、Ｆは、分離している前記送信機および受信機の距離、すなわち焦点間距離である。
前記楕円の短軸（Ｂ）を次式を用いて決定するステップと、
【数９】
0004043878.tif 000019
前記散乱領域から局所座標系の原点までの角度（Ψ）の正弦および余弦を、次式を用いて決定するステップと、
【数１０】
0004043878.tif 000020
および、
【数１１】
0004043878.tif 000021
前記角度の正弦および余弦ないし前記楕円の長軸および短軸から局所座標系における散乱領域の位置を決定するステップとを具備し、
前記散乱領域における流速を決定するステップは、後方散乱レーダシステムを発展させたコンピュータプログラムを用いて次式で表されるドップラーシフト（ｆD）を算出するステップを具備し、
【数６】
0004043878.tif 000022
更に、算出されたドップラーシフト（ｆD）を用いて、次式で表される散乱領域を通過する双曲線に沿った流速（Ｖh）を決定するステップを具備することを特徴とするバイスタティックレーダシステムを用いて表面流ベクトルをマッピングする方法。
【数７】
0004043878.tif 000023
ここで、λは、前記レーダ信号の波長、ｇは、重力加速度（９．８０６ｍ／ｓ2）、Ｖr  は、後方散乱レーダシステムを発展させた前記コンピュータプログラムを用いた擬似後方散乱過程と、前記複数の送信機および受信機によって探査され、計測された海洋状況とから抽出された擬似放射流速であり、更に、θは、前記送信機および前記散乱領域を結ぶ線分と、前記受信機および前記散乱領域を結ぶ線分との間の角であるバイスタティック角である。
【請求項２】
前記散乱領域への前記方位角（φ）を決定するステップは、多重信号分類（ＭＵＳＩＣ）方位発見アルゴリズムを用いることを特徴とする請求項１記載の表面流ベクトルをマッピングする方法。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は一般的なレーダシステム、特に、水量の多い領域の表面流をマッピングするためのシステムおよびそのマッピング方法に関する。
【０００２】
【従来の技術】
ＭＦ、ＨＦ、ＶＨＦおよび低ＵＨＦ帯で動作する低周波の後方散乱レーダシステムは、回流や川に沿った流れ等の水面の流速をマッピングするために、広く用いられている。
【０００３】
流速を測定するために、もっとも一般的に用いられているのは、ドップラー効果を利用する後方散乱システムである。後方散乱ドップラーシフトは、以下の２つの効果をもたらす。（ｉ）ブラッグ散乱、すなわち、波長が正確にレーダ波長の半分である波の動きを捉えることができる。（ii）レーダの伝搬方向の表面波を輸送する内在的な流れの構成を捉えることができる。後方散乱レーダにおいて、波のドップラーは、単一の周波数で生じ、測定されたドップラーシフトから所望の放射状の流速を得るまで減衰することがよく知られている。このように、波は、内在的な流れをトレースするのに役立つ。つまり、これらの波についての情報は、探知されるものではなく、レーダターゲットに供するためだけに必要となるものである。
【０００４】
散乱領域の範囲および距離は、送信レーダ信号と受信レーダ信号との時間遅延により得られる。例えば、米国特許５３６１０７２号により具体的に記載されている。範囲処理の後、それぞれの範囲の領域についてのエコー時間系は、ドップラースペクトルといくつかの受信アンテナまたは要素とを得るために、フーリエ変換される。流れを抽出するために用いられるブラッグピークは、単一かつ分離している。それぞれのドップラースペクトルビン（流れの放射速度として定義できる）において、エコーの方位角は、アルゴリズムを決定する方位を用いて受信アンテナの指向性より得られる。好適な方位決定のアルゴリズムとしては、例えば、米国特許５９９０８３４号に記載された多重信号分類（ＭＵＳＩＣ）方位発見アルゴリズムがある。このように、後方散乱レーダは、放射流速が距離と方位角との関数としてマップされた極座標系の測定結果を得る。
【０００５】
単一の後方散乱レーダシステムは、２次元の水平流速ベクトルである放射成分だけ測定することができる。それゆえ、通常２つの後方散乱レーダシステムが１対として用いられ、海岸に沿って数１０キロメートルごとに独立して動作させる。公知であるジオメトリと相互に観測される散乱領域の位置とに基いて、対象区域にオーバーラップした２つの結果の放射速度成分が結合されて、合計速度ベクトルマップが生成される。このように、従来システムの欠点の１つは、水平速度マッピングのために、後方散乱システムによる多重化が必要であり、コストがかかることである。
【０００６】
従来システムの他の欠点としては、後方散乱システムは、合計ベクトルあるいは２つの後方散乱レーダシステムが１つのラインに結合する付近を分離することができない。これは、双方が同じ速度成分として計測されるためである。このことは、海岸沿い、川を跨ぐ箇所、湾口において測定される領域が欠けてしまうこととなり、特に問題である。
【０００７】
水面速度をマッピングするための代替方法として、１つまたはそれ以上の受信機が離れた所にあり、２つまたはそれ以上の送信機で動作されるバイスタティックレーダシステムが用いられる。合成エコーは、それぞれ独立した送信機−受信機対から同時に放射された信号に基いて設定されるため、散乱領域における非平行な速度成分を構築することができる。この方法の主な欠点は、個々の送信機および受信機を地理的に同期させることがコスト高となることである。従前のバイスタティックシステムは、高安定のセシウムまたはルビジウムの時間標準を採用するか、あるいは、個々の位置におけるコヒーレント信号を維持するための過制御されたフェーズロックループを採用していた。このような方法は、高価であることに加えて、海洋探査システムが日常晒される過酷な環境には、不向きである。
【０００８】
後方散乱システムに比べて、バイスタティックレーダシステムを使用することの、さらなる根本的問題は、バイスタティックレーダシステムにおいては、送信機−受信機対について楕円の共焦点を見つけるために一定の時間遅延が発生することである。この結果、波の動作によるドップラーシフトは、もはや後方散乱レーダシステムのように一定ではなくなり、一定の時間遅延の生じた外形の位置が変化してしまう。
【０００９】
【発明が解決しようとする課題】
したがって、本発明は、かかる従来技術の問題点を解決するべくなされたもので、より広い領域に渡って、従来の後方散乱レーダシステムではシステムのジオメトリが通常モニタリング不可能な領域において正確にマッピングすることが可能なバイスタティックレーダシステムおよび表面流ベクトルをマッピングする方法を提供することを目的とする。また、本発明は、頑強かつ安価であるバイスタティックレーダシステムおよび表面流ベクトルをマッピングする方法を提供することを他の目的とする。さらに、現存する海洋探査レーダシステムに容易かつ安価にアップグレードすることができるバイスタティックレーダシステムおよび表面流ベクトルをマッピングする方法を提供することを他の目的とする。本発明のバイスタティックレーダシステムおよび表面流ベクトルをマッピングする方法によって、従来技術の課題を解決し、上記効果を得ることができる。
【００１０】
【課題を解決するための手段】
本発明は、バイスタティックレーダシステムおよび多量の水の表面流をマッピングするための方法を供給する。本発明は、広い領域に渡って、従来の後方散乱レーダシステムでは、有効にマッピングできなかった領域の流速を安定して積算することを提供する。
【００１１】
本発明の１つの目的は、従来の後方散乱レーダシステムを用いると、不安定で、通常はモニタリングが不可能な領域において、より広範囲で、正確なバイスタティックレーダシステムを用いて、海洋流の情報を供給することである。
【００１２】
本発明のさらなる目的は、複雑性、コストおよびＧＰＳ時間信号を用いたバイスタティックレーダシステムの送信機および受信機の変調信号を同期するための放射スペクトル資源の使用を減らすことである。
【００１３】
本発明のさらなる目的は、従来の後方散乱レーダシステムにおいて用いられる時間遅延または範囲、方位および速度の関係をバイスタティックレーダシステムの楕円／双曲線のジオメトリを用いた所望の流速マップに変換するための方法およびコンピュータプログラムを供給することである。
【００１４】
本発明の一側面によると、バイスタティックレーダシステムは、レーダ信号を送受信するための多数の送信機および受信機を具備する。好ましくは、１つの送信機は、少なくとも１つの受信機から離れて位置しており、送信機および受信機は、それぞれローカルオシレータを具備する。このローカルオシレータは、前記複数の送信機および受信機が、前記複数の送信機および受信機の間で干渉を生じさせるためにＧＰＳ信号に組み合わされる。さらに好ましくは、前記送信機および受信機は、海洋状況を探査し、計測するように構成され、バイスタティックレーダシステムは、前記海洋状況における情報を引き出すための信号処理手段をさらに具備する。送信機および受信機によって探査され、計測された海洋状況は、例えば、表面流速ベクトルを具備する。
【００１５】
前記信号処理手段は、前記バイスタティックレーダシステムを用いて測定されたドップラーシフト（ｆD）を用いて、次式で表される散乱領域内の流速を決定するように構成される。
【数１４】
0004043878.tif 000002
ここで、λは、前記レーダ信号の波長、ｇは、重力加速度（９．８０６ｍ／ｓ2）、およびθは、前記送信機および前記散乱領域を結ぶ線分と、前記受信機および前記散乱領域を結ぶ線分との間の角であるバイスタティック角である。前記流速（Ｖh）は、前記散乱領域を通過し、前記送信機および受信機を共焦点とする楕円に垂直な双曲線に沿って決定され、前記楕円は、測定されたレーダエコー時間遅延に等しい一定時間遅延（Ｄ）を有する。具体的には、ドップラーシフト（ｆD）は、バイスタティックレーダシステムを用いて直接的に計測される。他の具体例としては、前記信号処理手段は、前記バイスタティックレーダシステムを用いて測定されたドップラーシフト（ｆD）を用いて、次式で表される散乱領域内の流速を決定する。
【数１５】
0004043878.tif 000003
ここで、λは、前記レーダ信号の波長、ｇは、重力加速度（９．８０６ｍ／ｓ2）、θは、前記送信機および前記散乱領域を結ぶ線分と、前記受信機および前記散乱領域を結ぶ線分との間の角であるバイスタティック角、およびＶhは、前記散乱領域を通過する双曲線に沿った前記流速である。
【００１６】
好ましくは、前記表面流速ベクトルは、前記表面流速ベクトルとほぼ同じドップラースペクトル領域上の速度を有する波動から独立しているように構成される。さらに好ましくは、前記バイスタティックレーダシステムは、前記送信機から離れた位置にある前記受信機が結合する線に沿った領域に合計流ベクトルを供給するように構成される。
【００１７】
本発明の他の側面によると、複数の送信機および受信機を具備し、少なくとも１つの送信機が少なくとも１つの受信機から離れて位置するレーダシステムを用いて表面流ベクトルをマッピングする方法である。本方法は、（ｉ）前記受信機でエコーを生成するために、前記送信機から散乱領域内の波に向けて送信されたレーダ信号を散乱させるステップと、（ii）方位決定アルゴリズムを用いて前記散乱領域への方位角（φ）を決定するステップと、（iii）前記送信機から前記受信機までのレーダエコー時間遅延を測定するために、送信後の時間に対してサンプリングするステップと、（iv）前記散乱領域の位置を決定するステップと、（ｖ）前記散乱領域における流速を決定するステップとを具備する。
【００１８】
前記散乱領域における流速を決定するステップは、送信機および散乱領域および受信機を結ぶ線の間の角であるバイスタティック角、レーダエコーのドップラーシフト（ｆD）、および散乱領域の位置に基いている。前述されるように、ドップラーシフトは、バイスタティック送信機−受信機対のドップラーシフトから直接的に測定されるか、受信機の１つから測定された仮想の放射後方散乱の速度成分より算出される。
【００１９】
散乱領域の位置を決定するためのステップは、（ｉ）一定時間遅延（Ｄ）が前記測定されたレーダエコー時間遅延と等しくなるように、前記散乱領域が通過し、前記送信機および受信機が共焦点となる楕円の長軸（Ａ）を決定するステップと、（ii）前記楕円の短軸（Ｂ）を決定するステップと、（iii）前記散乱領域から局所座標系の原点までの角度（Ψ）の正弦および余弦を決定するステップと、（iv）前記角度の正弦および余弦ないし楕円の長軸および短軸から局所座標系における散乱領域の位置を決定するステップとを具備する。
【００２０】
さらに、他の側面においては、本発明は、上記の方法を実現するために、コンピュータシステムとともに使用するためのコンピュータプログラムを提供する。本コンピュータプログラムは、読み込み可能な記憶媒体およびコンピュータプログラム構造が組み込まれ、前記コンピュータ記憶媒体は、複数の送信機および受信機を具備し、少なくとも１つの送信機が少なくとも１つの受信機から離れて位置するレーダシステムを用いて表面流ベクトルをマッピングするために、所定の仕様に基いて前記コンピュータシステムを指示するプログラムモジュールを具備する。また、本コンピュータプログラムは、（ｉ）レーダ制御サブルーチン＝プログラムモジュール、（ii）方位角決定サブルーチン＝プログラムモジュール、（iii）エコー時間遅延サブルーチン＝プログラムモジュール、（iv）位置決定サブルーチン＝プログラムモジュール、および（ｖ）流速サブルーチン＝プログラムモジュールを具備する。
【００２１】
本発明のコンピュータプログラムおよび方法は、以下のどれか１つまたは全ての利点を有する。
【００２２】
（ｉ）   複雑性、コストおよびＧＰＳ時間信号を用いたバイスタティックレーダシステムの送信機および受信機の変調信号を同期するための放射スペクトル資源の使用を減らすことができる。
【００２３】
（ii）   送信機および／または受信機を増加させて発展させることによって流れマッピングのための従来の後方散乱レーダシステムの範囲領域と正確さを迅速かつ比較的安価に改良することができる。
【００２４】
（iii）  コンピュータもエアコンも必要なく、分散配置には好適なソーラーエネルギーを利用することが可能な小型の送信システムをデザインすることができる。
【００２５】
（iv）   結果として、方位角を正確に測定する能力に影響する複雑なアンテナパターンを破壊することのない、ブイ、沖合の構造物、または建物の屋根上に取り付けることが可能な垂直なホイップ送信アンテナから信号を放射することができる。
【００２６】
【発明の実施の形態】
本発明は、大洋、海、湾、港、川および湖のような水量の多い領域の表面の流れをマップするためのバイスタティックレーダシステムおよびレーダシステムを用いて表面流ベクトルをマッピングする方法である。
【００２７】
本発明に係るバイスタティックレーダシステムは、図１に示されている。図１はレーダ信号を送受信するための複数の送信機１０２および複数の受信機１０６を含むバイスタティックレーダシステム１００の模式的ブロック図である。少なくとも１つの送信機１０２は、少なくとも１つの受信機１０６から離れて位置し、本実施形態のようにＧＰＳ信号を使用することによって送信機１０２と受信機１０６とが干渉することが必要である。明確のために、レーダシステムの多くについての詳細は、広く知られており、本発明に関係しないため、省略されている。
【００２８】
図１は、バイスタッティックレーダシステム１００のハードウェアの構成および信号処理過程を示している。角が角張ったブロックをハードウェアの構成として表示し、角の丸いブロックを受信したレーダデータを処理するキー信号処理過程として表示しており、デジタルコンピュータ上のソフトウェアプログラムによってリアルタイムで実行させることができる。
【００２９】
図１に示されるように、バイスタッティックレーダシステム１００は、海洋あるいは水の表面状態を調査することができるいくつかの独立した受信機１０６を有し、受信部１０４と協同する送信機１０２を具備している。受信機１０６の受信アンテナ１０８は、送信機１０２の送信アンテナ１１０から送信され、散乱地点あるいは散乱領域（図示せず）において波によって散乱したレーダ信号のエコーを受信する。信号処理機１１２は、海洋の状態に関する情報を決定し、表示するためにレーダエコーから導き出された情報を処理する。このような情報は、表面流速ベクトルを具備している。バイスタティックモードにおいて、図２に示されるように、送信アンテナ１１０によって、地理的に離れている少なくとも１つの受信アンテナ１０８へ指示される。つまり、送信アンテナ１１０および少なくとも１つの受信アンテナ１０８は、送信アンテナ１１０を繋ぐ線と散乱領域１１４との間の角度および受信アンテナ１０８を繋ぐ線と散乱領域１１４との間の角度が０にならないように位置する。好ましくは、バイスタッティックレーダシステム１００によってカバーされる領域を最大限にするために、すべての受信機１０６または受信アンテナ１０８は、送信アンテナ１１０から離れているとともに、お互いの受信アンテナ１０８から離れているように構成される。送信機１０２と受信機１０６との間は、物理的に接続されておらず、電信による接合や通信路によって連結される。バイスタッティックレーダシステム１００が必要とする周波数およびマップ対象に依存する両者の距離は、約２〜２００ｋｍにまで及ぶ。後方散乱システム１１６においては、比較と例示のために、送信機１０２と送信アンテナ１１０とは、受信機１０６と同じ位置に配置されている。
【００３０】
送信機１０２および受信機１０６の電子的構成は、それぞれのアンテナ１０８，１１０の近くにある必要はない。また、受信機１０６は、信号処理機１１２の近くにある必要もない。例えば、受信機１０６から離れた場所にある施設に、ユーザあるいはオペレータが海洋の状態に関する情報にアクセスする場所に信号処理機１１２を配置するのが好ましい。信号処理機１１２は、受信機１０６に陸線または無線通信路を経て接続可能である。同様に、送信機１０２および受信機１０６の電子的構成は、アンテナ１０８，１１０が設置され、ショートケーブルまたはワイヤを経てアンテナ１０８，１１０に連結された囲みの中、すなわち建物内に設置可能である。
【００３１】
好適には、送信機１０２および送信アンテナ１１０は、例えば、光電池、燃料電池またはバッテリーのような電源を自ら具備している。好ましくは、送信機１０２は、ボート、はしけ、プラットフォームまたはブイのような、耐水性の基盤上に設置され、受信アンテナ１０８の一部または全部がリングで囲まれている。例えば、送信機１０２は、港または湾１１８内のブイ上にあり、受信アンテナ１０８は、図２に示すように、湾の周りの地面１２０上に建てられた建物のように、他のブイ上または岸辺に設置される。
【００３２】
送信機１０２は、ローカルオシレータ（図示せず）を有するデジタルシンセサイザ１２２を具備し、ローカルオシレータは、送信機１０２により送信されたレーダ信号を生成するためにキャリア周波数、すなわちキャリア波を発生させ、キャリア波を変調するための信号を発生させる。バイスタティック構成のためのレーダ信号は、例えば、米国特許５３６１０７２号に記載されるように連続的な２０ＭＨｚのキャリア波が変調されて形成される。変調されたレーダ信号は、デジタルシンセサイザ１２２から送信アンプまたはパワーアンプ１２４を通過し、送信アンテナ１１０から放射される。送信アンテナ１１０は、通常、ブロードビームを供給するために組み込まれ、１２０°から３６０°の方位角を有し、海洋の状態に関する情報が必要とされる海洋または水の広範囲に照射される。
【００３３】
受信機１０６は、レーダエコーを受信し、デコードするための受信モジュール１２６を具備し、受信したエコーから情報を引き出すためのいくつかの信号処理過程または信号処理機を具備する。受信モジュール１２６は、受信したエコーとローカルオシレータ（図示せず）で発生した信号とを混合し、エコーをデコードする。一般に、信号処理機は、ハードウェアおよびソフトウェアの要素を有している。好ましくは、受信モジュール１２６は、それぞれアナログ−デジタル（Ａ／Ｄ）変換器（図示せず）を具備し、信号処理機によってデジタル処理できるようにデジタル時系列に出力する。デジタル信号処理機１２８は、時間遅延領域の位置を決定し、レーダエコーおよび雑音の時系列をそれぞれの領域枠ごとに出力する。コヒーレント積分やコヒーレント処理として知られる処理を行うドップラー信号処理機１３０は、潜在的なターゲットエコーを２次元空間の領域およびドップラー周波数に順序立てる。最終処理過程は、方位角処理機１３２で行われ、エコー方位角を決定する前に、有用な信号を明確化し、背景雑音や散乱から抽出する。これら３つの信号処理機１２８，１３０，１３２は、それぞれ米国特許５３６１０７２号および５９９０８３４号に詳述されており、これらを参照することにより、具体化される。コヒーレントレーダ処理のために、放射信号と受信信号との間の正確なタイミングが必要となるため、送信機１０２および受信機１０６のそれぞれは、受信モジュール１２６で生成された信号を送信機１０２において変調信号を同期させるためのＧＰＳ衛星受信機１３４をさらに具備する。これらＧＰＳ衛星受信機１３４は、小型のデバイスであり、それ自身に衛星アンテナ１３６を装備し、複数のベンダー、例えば、カリフォルニア・サニーベイルにあるTrimble Navigation Limited社のThunderbolt GPS Disciplined Clockから商業的に利用できる。送信機１０２および受信機１０６の変調信号を同期するためのＧＰＳタイミング信号の利用方法の１つとしては、例えば、２００１年８月２８日出願の米国出願６０／３１５５６７号に記載されており、これを参照して具体化される。ローカルフェーズ固定オシレータ（ＰＬＯ）は、送信機１０２および受信機１０６内の共通の基準を供給するためにＧＰＳ衛星受信機１３４からのＧＰＳ信号を固定する。ＰＬＯは、一機関に渡って１００万につき１０12〜１０13パーツの正確さで構築される電子的フライホイールの役割を果たしている。
【００３４】
信号処理機１１２は、本発明が海洋の状態の情報を決定し、表示するために、レーダエコーから抽出された情報を処理するようにプログラムされたパソコンのような汎用のデジタルコンピュータを具備する。また、信号処理機１１２は、情報が処理されるいくつかの要素または過程を包含し、方位角決定過程１３８、バイスタティック散乱地点決定過程１４０、バイスタティック流速過程１４２および合計流速ベクトル過程１４４を具備する。これらの過程およびその作用については、以下で図１および図３を参照しつつ説明する。
【００３５】
第１の過程、すなわち、方位角決定過程１３８において、方位範囲および散乱領域１１４ごとのドップラーは、受信機１０６ごとから受信されるデジタルターゲットエコー信号から決定される。好ましくは、この決定は、米国特許５９９０８３４号に記載の多重信号分類（ＭＵＳＩＣ）方位発見アルゴリズムのような、好適な方位決定アルゴリズムを用いてなされる。本過程、すなわち、方位角決定過程１３８は、擬似後方散乱過程と呼ばれている。というのは、送信機と受信機とが同じ所に位置する後方散乱レーダシステムにおいて、本過程１３８からの出力は、放射表面流速と範囲と方位とで構成される極座標系のマッピングとして表される。それぞれのドップラーシフト（あるいはスペクトルビン）は、後方散乱レーダシステムに向かって、あるいは、離れるように動く、散乱領域１１４の速度成分に関係する。これは、放射速度Ｖrとして定義される。１つの具体例において、擬似後方散乱過程１３８の出力は、後方散乱レーダシステムの放射表面流速マップに一致するバイステティック合計流速ベクトルマップに変換される。この放射表面流速マップから合計流速ベクトルマップへの変換の数学的および物理的理解は、以下に図３を参照しつつ説明する。図３に示されるように、本発明のバイスタティックレーダシステム１００の送信機１０２は、座標ｘ，ｙ＝１４．１，１４．１にある円１４６に位置している。受信機１０６の位置は、座標ｘ，ｙ＝０，０の原点にある円１４８によって示される。これは、図示して比較するためのモノスタティック、すなわち、この受信機１０６が一緒に作動する後方散乱レーダシステムに形成される後方散乱送信機の位置ともなる。後方散乱レーダシステムにとって、一定の時間遅延の外形は、ターゲットのエコー信号の送受信間における円１４８で示される後方散乱レーダシステムの位置と同軸の円１５０，１５２として形成される。対して、バイスタティックレーダシステム１００にとって、一定時間遅延の外形は、円１４６の送信機１０２と円１４８の受信機１０６とを共焦点とする一群の楕円１５４，１５６で定義される。その結果、楕円群の焦点間距離はＦ、すなわち、送信機１０２および受信機１０６間の距離で示され、図３においては、２０単位の距離に等しい。
【００３６】
後方散乱レーダシステムにおいては、測定された３つの出力があり、そのアナログバイスタティックな対象物は、所望のバイスタティック合計流速ベクトルマップを定義するために探知される。その３つとは、すなわち、（ｉ）散乱領域１１４の範囲、（ii）散乱領域１１４に対する方位角φ、（iii）散乱領域１１４における水の（流れの）放射速度Ｖr、である。バイスタティックレーダシステムにおいて、これに一致する出力、すなわち、対象物は、（ｉ）散乱領域１１４の位置（水平位置を定義するために２つの座標が必要である）および（ii）散乱領域１１４でのバイスタティックレーダから適用可能な流速成分である。
【００３７】
バイスタティック構成においては、局所座標系（図示せず）を採用し、その原点は、送信機１０２と受信機１０６との間の線分Ｆの中点である。この局所座標系のｘ’軸は、円１４６の送信機１０２の方向に取られ、ｙ’軸は、これに半時計回りに垂直な方向に取られる。
【００３８】
局所座標系におけるＲ，ｘ’，ｙ’で表される散乱領域１１４の位置は、以下の数学的ステップを用いることにより、バイスタティック散乱地点決定過程１４０において決定される。
【００３９】
測定されたレーダエコー時間遅延Ｄは、散乱領域１１４を通過する楕円を示し、受信機１０６から散乱領域１１４までの距離Ｐ、散乱領域１１４から送信機１０２までの距離Ｑ、および送信機１０２と受信機１０６との距離Ｆに関係し、次式で表される。
【数１６】
0004043878.tif 000004
【００４０】
時間遅延Ｄが与えられるときの散乱領域１１４が通過する楕円の長軸Ａは、次式で表される。
【数１７】
0004043878.tif 000005
【００４１】
時間遅延Ｄが与えられるときの散乱領域１１４が通過する楕円の短軸Ｂは、ｙ’軸に沿って存在し、次式で表される。
【数１８】
0004043878.tif 000006
【００４２】
散乱領域１１４の方位角φが測定されるとともに、方向決定アルゴリズムを用いた方位角決定過程１３８から知られることにより、後述する式の分母Ｄｅｎは、次式で表される。
【数１９】
0004043878.tif 000007
【００４３】
局所座標軸系の原点ｘ’＝０，ｙ’＝０からＲにおける散乱領域１１４までの角度の余弦は、次式で与えられる。
【数２０】
0004043878.tif 000008
【００４４】
局所座標軸系の原点ｘ’＝０，ｙ’＝０からＲにおける散乱領域１１４までの角度の正弦は、次式で与えられる。
【数２１】
0004043878.tif 000009
【００４５】
散乱領域１１４の所望の位置は、局所座標軸系において次式で与えられる。
【数２２】
0004043878.tif 000010
【００４６】
流速Ｖhは、散乱領域１１４を通過する楕円に垂直な双曲線に沿った方向を有しており、以下に示す数学的ステップを用いたバイスタティック流速過程１４２において決定される。
【００４７】
図３の円１４８にある送信機−受信機対から見えるエコー信号の後方散乱ドップラーシフトｆDは、レーダおよび流速の放射成分の波長に関係している。注目すべきは、ドップラーシフトｆDは、前述のドップラー信号処理機１３０において既に計測されていることである。後方散乱構成からのエコー信号のドップラーシフトｆDは、次式で表現される。
【数２３】
0004043878.tif 000011
ここで、λは、既知のレーダ波長、ｇは、重力加速度（９．８０６ｍ／ｓ2）、そして、Ｖrは、レーダに向かってくる、または、レーダから離れていく流速の放射成分である。±の記号は、ＨＦレーダによって見られる２つのスペクトルブラッグピークを示す。レーダ波長が向かってくる波（＋）と、離れていく波（−）とに分けられているためである。双方の波は、重複する要素を供給し、精度を向上させて、波の情報を抽出するために用いられる。
【００４８】
バイスタティック構成における円１４８にある受信機１０６および円１４６にある送信機１０２から見えるエコー信号のバイスタティックドップラーシフトｆDは、次式で表される。
【数２４】
0004043878.tif 000012
ここで、λは、レーダ波長、ｇは、重力加速度（９．８０６ｍ／ｓ2）、θは、送信機１０２および散乱領域１１４を結ぶ線分と受信機１０６および散乱領域１１４を結ぶ線分とのなす角であるバイスタティック角、および、Ｖhは、バイスタティックレーダシステム１００によって測定可能な流速である。
【００４９】
後方散乱ドップラーシフトｆDを示す式８や式９の第１項によると、ドップラーシフトｆDは、波動に依存している。つまり、ブラッグ波による散乱は、送信機１０２と受信機１０６とを分ける基準線に向かって進むことにより、一定時間遅延Ｄが増大し、基準線から離れて進むことにより、一定時間遅延Ｄが減少する。この項は、バイスタティック構成の複雑性を示唆している。波動に依存するドップラーシフトｆDは、後方散乱構成にとって、もはや一定ではなく、範囲領域、すなわち、一定時間遅延Ｄの楕円に沿って位置を変える。このように、流れに依存するドップラーシフトから波動に依存するドップラーシフトを分離することが必要である。従来の明確でない海洋の状態を観測するバイスタティックレーダシステムを用いると、この分離は、困難である。バイスタティックシステムが測定できる速度成分は、後方散乱ジオメトリと対照的にＶhと呼ばれる。なぜなら、その速度成分は、散乱地点における楕円に垂直な双曲線に沿った方向を有しているからである。このベクトルＶhは、図３の点Ｒにおいて例示されている。式９の有効性を確認するために、送信機１０２および受信機１０６がともに移動する様子に注目する。構成がモノスタティック、すなわち、後方散乱レーダと同様になる。つまり、バイスタティック角が０に向かうにつれて、式９は、式８に近づいていく。
【００５０】
双曲線に沿った流速Ｖhは、次式を用いて測定されたドップラーシフトから取り出される。
【数２５】
0004043878.tif 000013
ここで、λは、レーダ波長、ｇは、重力加速度（９．８０６ｍ／ｓ2）、θは、バイスタティック角、そしてＶhは、散乱領域１１４を通過する双曲線に沿った流速である。
【００５１】
次に、局所的なｘ’，ｙ’座標系に関してのＶhの方向は、以下のように決定される。想起されるバイスタティックドップラーシフトは、ドップラー信号処理機１３０において既に測定されている。以下のステップでは、必要なバイスタティック角を決定する方法について示す。
【００５２】
次式は、ベクトルＰ，Ｑおよびそれらの単位ベクトルｐ’，ｑ’を送信機１０２および受信機１０６の位置からそれぞれ散乱地点に向かった方向に取って、定義する。
【数２６】
0004043878.tif 000014
【００５３】
散乱地点における楕円から外向きの標準ベクトルＮおよび単位標準ベクトルは、次式によって定義される。
【数２７】
0004043878.tif 000015
【００５４】
単位標準ベクトルのこれらの集合から、式９において速度成分を決定するために必要なθ／２の余弦は、次式を用いて明らかとなる。
【数２８】
0004043878.tif 000016
【００５５】
θ／２の余弦を式９に代入して、バイスタティックベクトルＶhの大きさおよび局所座標系のｘ’，ｙ’に沿った成分（これをｕ’，ｖ’と定義する）を用いて表すと次式が与えられる。
【数２９】
0004043878.tif 000017
【００５６】
これらは、後方散乱、すなわち、送信機−受信機ジオメトリによって定義された局所座標系の点ｘ’，ｙ’における、バイスタティックレーダシステム１００によって測定可能な、最後の流速成分である。多くのドップラースペクトルビンに基いて、双曲線に沿った速度成分Ｖhのマップは、測定された海洋のエコーデータから構築することができる。
【００５７】
この送信機−受信機対に指向性を有するマップで定義される点ｘ’，ｙ’におけるこれらの構成ｕ’，ｖ’は、よく知られている数学的テクニックによって、変換、回転させることにより、容易に他の座標系にすることができる。例えば、海岸線が、図３の太い破線１５８に沿って指向されるとき、海岸線に指向性を有する系を取ることがより好ましい。また、他の例として、地球を基準とした緯度−経度（北緯−東経）による系を選択してもよい。
【００５８】
上記に概説したようなバイスタティックレーダシステム１００から抽出された流れのマッピングに基いた、双曲線に沿った速度Ｖhは、水平の流れを完全に表現してはいない。単一の後方散乱レーダを有する後方散乱レーダシステムから抽出可能でもない。両システムは、マップ上のそれぞれの位置における２次元の流ベクトル１つの構成としてのみ計測する。それゆえ、合計流速ベクトルマップは、最終過程、すなわち、合計流速ベクトル過程１４４において決定される。
【００５９】
従来の後方散乱レーダシステムは、放射流速マップを作成するために、２つの異なる方向から同じ地点の表面を見るための少なくとも２つの完全な後方散乱レーダを必要とする。さらに、上述のように、２つの後方散乱レーダシステムが含まれる「ベースライン」およびその近くの流ベクトルを決定することができない。例えば、第２の後方散乱レーダが図３の海岸線を示す太い破線１５８に沿って位置するとき、例えば、位置（２０，０）において、後方散乱レーダシステムは、２つのベクトル間の角度が０となる海岸線付近の重要な領域をカバーした流ベクトルマップを作成することができない。
【００６０】
これに対して、本発明のバイスタティックレーダシステム１００は、円１４８にある送信機−受信機対に、円１４６にある送信機１０２が増設されることにより、正確な合計流ベクトルを生成し、海岸線沿いの領域を含むバイスタティックレーダシステム１００によってカバーされる全領域を積算することができる。円は、一定時間遅延の外形（すなわち、範囲領域）を定義し、そこで、後方散乱形状に動作する円１４８での送信機−受信機対からの放射マップが放射ベクトルＶrを生成する。これと同時に動作する、円１４６にある送信機１０２および円１４８にある受信機１０６のバイスタティック構成からの測定は、同じ海岸地点にある双曲線に沿ったベクトルＶhを生成する。示されるように、これらのベクトルは、もはや平行ではない。実際、それらは、互いに少なからぬ角度を有している。このように、本発明に係るバイスタティックレーダシステム１００は、２つの送信機１０２、単一の受信機１０６および信号処理機１１３を具備し、海岸線沿いの領域を含むシステムによってカバーされる全領域において合計速度ベクトルマップを供給するのに有効である。
【００６１】
本発明に係る表面流速ベクトルをマッピングするための方法を図４を参照しつつ説明する。図４は本発明の一実施形態に係るバイスタティックレーダシステムを動作させるための方法についてのステップを示すフローチャートである。本方法では、レーダ信号は、送信機１０２から放射され、散乱領域で反射し、受信機１０６においてエコーが生成される（ステップ１６０）。散乱領域の方位角φは、ＭＵＳＩＣ方位発見アルゴリズムや他の好適な方位角決定アルゴリズムを用いて、決定される（ステップ１６２）。送信機１０２から受信機１０６までのレーダエコー時間遅延を測定するために、送信後の時間に対してレーダエコーがサンプリングされる。散乱領域の位置が決定され（ステップ１６６）、散乱領域にある表面流速は、上記式１０を用いて決定、算出される（ステップ１６８）。散乱領域にある流速を決定するステップ１６８において、レーダエコーのためのドップラーシフトｆDを決定することも含めてもよい。具体例の１つとして、このドップラーシフトｆDは、バイスタティックレーダシステムの分離されている送信機１０２および受信機１０６を用いて直接的に計測され、ドップラーシフトｆDは、上記式８を用いて方位角φおよび受信機１０６に関する放射流速から計算される。
【００６２】
本発明の実施形態に係る散乱領域の位置を決定するための方位を図５を参照しつつ説明する。図５に示すように、本方法は、次の手順で達成される。（ｉ）一定時間遅延Ｄが測定されたレーダエコー時間遅延と等しくなるように、散乱領域が通過し、送信機および受信機が共焦点となる楕円の長軸Ａを決定する（ステップ１７０）。（ii）楕円の短軸Ｂを決定する（ステップ１７２）。（iii）散乱領域から局所座標系の原点までの角度Ψの正弦および余弦を決定する（ステップ１７４）。そして、（iv）前記角度Ψの正弦および余弦ないし楕円の長軸および短軸から、局所座標系における散乱領域の位置を決定する（ステップ１７６）。長軸を決定するステップ１７０は、上記式２を用いることにより達成される。短軸を決定するステップ１７２は、上記式３を用いることにより達成され、角度Ψの正弦および余弦を決定するステップ１７４は、式４および式５を用いることで達成される。最後に、局所座標系における散乱領域の位置を決定するステップ１７６は、上記式７を用いることで達成される。
【００６３】
上述のそれぞれのステップを処理していくためのコンピュータプログラムについて図６を参照しつつ説明する。図６は本発明の実施形態に係るコンピュータプログラムの模式的ブロックダイヤグラムである。コンピュータプログラム１７８は、（ｉ）レーダ制御サブルーチン＝プログラムモジュール１８０、（ii）方位角決定サブルーチン＝プログラムモジュール１８２、（iii）エコー時間遅延サブルーチン＝プログラムモジュール１８４、（iv）位置決定サブルーチン＝プログラムモジュール１８６、および（ｖ）流速サブルーチン＝プログラムモジュール１８８を具備する。レーダ制御プログラムモジュール１８０は、送信機１０２から散乱領域に向けた散乱レーダ信号により、受信機１０６においてエコーを生成するためにバイスタティックレーダシステム１００を動作し、制御させるようなプログラムコードを有する。方位角決定プログラムモジュール１８２は、方位角決定アルゴリズムを用いて、散乱領域１１４についての方位角φを決定するためのプログラムコードを有する。エコー時間遅延プログラムモジュール１８４は、エコーの時間遅延を測定するためのプログラムコードを有する。位置決定プログラムモジュール１８６は、局所座標系内における散乱領域１１４の位置を決定するためのプログラムコードを有する。流速プログラムモジュール１８８は、散乱領域１１４における流速を決定するためのプログラムコードを有する。
【００６４】
１つの具体例として、流速プログラムモジュール１８８は、バイスタティックレーダシステム１００の地理的に離れた位置にある送信機１０２および受信機１０６を用いて直接測定されたドップラーシフトを用いたプログラムコードを有する。あるいは、流速プログラムモジュール１８８は、後方散乱における受信機１０６に関して測定された方位角φおよび放射流速Ｖrより、ドップラーシフトを計算するプログラムをさらに有してもよい。
【００６５】
以上の記述は、発明のより好ましい具体例を図および記述によって表したものである。これは、本発明を実施形態に開示された正確な形に限定するものではなく、また排他するものでもない。そして、上記より多くの変更およびバリエーションが想起され得ることは、明らかである。本実施形態は、本発明の原理を最もよく説明するために、選択し、記載したものであって、実用的な応用例は、本実施形態によって、他の当業者でも特定の使用を考慮することにより、好適となる本発明の様々な形態や様々な変化形態を最もよく利用することができる。例えば、２つ以上の送信機は、単一の受信機に組み合わせて動作させてもよい。同様に、複数の受信機を単一の送信機と組み合わせて動作させてもよく、この場合、無線封止である方が好ましい軍事活動において利点を有する。これらの可能性の全ては、本発明の他の実施形態を構成する。このように、請求の範囲によって定義される本発明の範囲は、前述の実施形態に限定されるものではない。
【００６６】
【発明の効果】
本発明に係るバイスタティックレーダシステム１００は、円１４８にある送信機−受信機対に、円１４６にある送信機１０２が増設されることにより、正確な合計流ベクトルを生成し、海岸線沿いの領域を含むバイスタティックレーダシステムによってカバーされる全領域を積算することができる。したがって、海岸線沿いの領域を含むシステムによってカバーされる全領域において合計速度ベクトルマップを供給するのに有効である。
【図面の簡単な説明】
【図１】レーダ信号を送受信するための複数の送信機および複数の受信機を含むバイスタティックレーダシステムの模式的ブロック図である。
【図２】本発明の一実施形態に係るバイスタティックレーダシステムの大量の水の回りにある送信機および受信機の位置を示す図である。
【図３】本発明の一実施形態に係る測定可能な流ベクトルと後方散乱およびバイスタティックレーダシステム双方のジオメトリとの関係を示す図である。
【図４】本発明の一実施形態に係るバイスタティックレーダシステムを動作させるための方法についてのステップを示すフローチャートである。
【図５】本発明の一実施形態に係る散乱領域の位置を決定するための方法についてのステップを示すフローチャートである。
【図６】本発明の実施形態に係るコンピュータプログラムの模式的ブロックダイヤグラムである。
【符号の説明】
１００…バイスタティックレーダシステム
１０２…送信機
１０６…受信機
１１２…信号処理機
１１４…散乱領域
【図１】
0004043878.tif 000024
【図２】
0004043878.tif 000025
【図３】
0004043878.tif 000026
【図４】
0004043878.tif 000027
【図５】
0004043878.tif 000028
【図６】
0004043878.tif 000029