(57)【要約】【課題】  海面により反射する電波を分析し、海流の流速の変化を高感度で検出し、津波を高感度で検出できる海洋レーダ装置を得る。【解決手段】  観測する海面を区切るセル毎に海面からの反射電波を分析し、前記反射電波のドップラ分布を前記セル毎に算出する信号処理装置と、地震の発生位置および規模から予測された津波の流速の海面上の分布をもとに、前記津波の流速が所定の流速から一定の範囲内である、前記津波の前縁部に位置する複数の前記セルを選択し、この選択した複数の前記セルからなる積分区間を設定する積分区間設定手段と、前記信号処理装置が算出したドップラ分布を前記積分区間について積分し、積分結果を算出する積分手段と、前記積分結果から、前記積分区間で振幅が最大となるドップラ周波数を求め、このドップラ周波数の変化を検出する検出手段とを備えるようにした。【選択図】          図２
2014098686.tif 000002
【特許請求の範囲】
【請求項１】
観測する海面を区切るセル毎に海面からの反射電波を分析し、前記反射電波のドップラ分布を前記セル毎に算出する信号処理装置と、
地震の発生位置および規模から予測された津波の流速の海面上の分布をもとに、前記津波の流速が所定の流速から一定の範囲内である、前記津波の前縁部に位置する複数の前記セルを選択し、この選択した複数の前記セルからなる積分区間を設定する積分区間設定手段と、
前記信号処理装置が算出したドップラ分布を前記積分区間について積分し、積分結果を算出する積分手段と、
前記積分結果から、前記積分区間で振幅が最大となるドップラ周波数を求め、このドップラ周波数の変化を検出する検出手段と
を備えた海洋レーダ装置。
【請求項２】
前記積分区間設定手段は、前記積分区間を、前記津波の流速が第１の流速から一定の範囲内の複数の前記セルからなる第１の積分区間と、前記津波の流速が第１の流速とは異なる第２の流速から一定の範囲内の複数の前記セルからなる第２の積分区間として設定し、
前記積分手段は、入力されたドップラ分布を前記第１の積分区間と前記第２の積分区間とについて積分し積分結果をそれぞれ算出し、
前記検出手段は、前記積分手段が算出した、前記第１の積分区間についての積分結果と前記第２の積分区間についての積分結果のそれぞれから、前記第１の積分区間と前記第２の積分区間のそれぞれで振幅が最大となるドップラ周波数を求め、求めたドップラ周波数の前記第１の積分区間と前記第２の積分区間との間の変化を検出する
請求項１に記載の海洋レーダ装置。
【請求項３】
前記検出手段は、前記積分結果から、前記積分区間で振幅が最大となるドップラ周波数を求め、このドップラ周波数の時間変化を検出することを特徴とする請求項１に記載の海洋レーダ装置。
【請求項４】
通常流予測手段が通常海流の予測速度から算出した、海面からの反射電波の通常海流によるドップラ周波数が入力され、前記信号処理装置が前記セル毎に算出したドップラ分布を、前記通常流予測手段が算出した通常海流によるドップラ周波数によりセル毎に周波数方向にシフトさせ、セル毎の補正されたドップラ分布を前記積分手段に出力する周波数シフト手段を備えた
請求項１ないし３のいずれか１項に記載の海洋レーダ装置。
【請求項５】
前記信号処理装置がセル毎に算出した前記ドップラ分布から、所定の周期の変動成分を抽出するドップラ変動周期分離手段と、
前記信号処理装置がセル毎に算出した前記ドップラ分布を、前記ドップラ変動周期分離手段が抽出した前記所定の周期の変動成分の周波数によりをセル毎に周波数方向にシフトさせ、セル毎の補正されたドップラ分布を前記積分手段に出力する周波数シフト手段を備えた
請求項１ないし３のいずれか１項に記載の海洋レーダ装置。
【請求項６】
観測する海面を区切るセル毎に海面からの反射電波を分析し、前記反射電波のドップラ分布を前記セル毎に算出する信号処理装置と、
地震の発生位置および規模から予測された津波の流速の海面上の分布をもとに、前記津波の波面に平行な方向に並んだ複数のセルを選択し、この選択した複数の前記セルからなる積分区間を設定する積分区間設定手段と、
前記信号処理装置が算出したドップラ分布を前記積分区間について積分し、積分結果を算出する積分手段と、
前記積分結果から、前記積分区間で振幅が最大となるドップラ周波数を求め、このドップラ周波数の変化を検出する検出手段と
を備えた海洋レーダ装置。
【請求項７】
観測する海面を区切るセル毎に海面からの反射電波を分析し、前記反射電波のドップラ分布を前記セル毎に算出する信号処理装置と、
地震の発生位置および規模から予測された津波の流速の海面上の分布をもとに、前記津波の流速が所定の流速から一定の範囲内である、前記津波の前縁部に位置する複数の前記セルを選択し、この選択した複数の前記セルからなる積分区間を設定する積分区間設定手段と、
前記信号処理装置が算出したドップラ分布を前記積分区間について積分し、積分結果を算出する積分手段と、
前記積分結果から、前記積分区間で振幅が最大となる第１のドップラ周波数を求め、この第１のドップラ周波数の変化を検出する第１の検出手段と、
前記反射電波のドップラ分布から、前記セル毎に振幅が最大となる第２のドップラ周波数を求め、この第２のドップラ周波数の変化を検出する第２の検出手段と、
を備えた海洋レーダ装置。
【請求項８】
前記信号処理装置は、前記セル毎の前記反射電波のドップラ分布を、
第１の受信位置で受信した第１の反射電波を分析し算出した第１のドップラ分布と
第２の受信位置で受信した第２の反射電波を分析し算出した第２のドップラ分布として算出し、
前記積分手段は、前記積分区間設定手段の設定した積分区間について、前記信号処理装置が算出した前記第１のドップラ分布と前記第２のドップラ分布とをそれぞれ積分し積分結果を算出し、
前記検出手段は、前記積分手段が算出した前記第１のドップラ分布の積分結果と前記第２のドップラ分布の積分結果のそれぞれから前記積分区間で分布が最大となるドップラ周波数を求め、この第１のドップラ分布の積分結果から求めたドップラ周波数の変化とこの第２のドップラ偏移の周波数分布の積分結果から求めたドップラ周波数の変化とをそれぞれ検出する
請求項１に記載の海洋レーダ装置。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、観測海域に電波を送信し、その反射信号を分析することで津波を検出する海洋レーダ装置に関する。
【背景技術】
【０００２】
津波は、その構成する海水の量が膨大なものであることなどから、大きな被害をもたらすものとなりうる。このため、防災の対象となる地域では、津波の到来前に十分な対処する時間を得るため、できるだけ遠方から津波を検出する必要がある。津波を探知する方法としては、固定的に設けられたレーダ（例えば、特許文献１参照）、波浪計、監視カメラなどが有り、また、特に遠方の津波を探知する方法としては、地震発生情報を基に津波の発生が予測される地域を航空機により監視する方法なども知られている（例えば、特許文献２参照）。これらの手段を比較すると、天候に左右されにくく、連続した時間の観測が可能で、監視対象を広い範囲にできるなどの理由で、防災の対象付近の固定した位置に設けたレーダ装置により監視を行なうことが他の方法に比べて有利な面が多い。
【０００３】
津波の周期は短くとも数分、長い場合には数時間になることもあり、波長は数ｋｍから数百ｋｍに及ぶ。また、津波は、孤立波として伝わり、その水平方向の流速成分は、津波の前縁で激しく変化するという特徴を持つ。このため、レーダにより津波を監視する方法としては、海洋に電波を送信し、その反射波を分析することで、海面の距離分解能等に対応したセル毎に流速を検出し、その流速のセル同士の速度偏差などにより流速の急激な変化を検出することで津波の検出を行なっている（例えば、特許文献１参照）。この為、津波検出の感度は、電波の海面でのレーダ電波の反射信号のドップラをいかに高感度で精度良く探知できるかに左右される。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開平８−２９２２７３号公報（段落００１９〜００５４）
【特許文献２】特開２００９−２２９４２４公報（第１〜第３図）
【発明の概要】
【発明が解決しようとする課題】
【０００５】
従来の海洋レーダ装置では、津波を遠距離から探知するためには、海面から反射する電波を受信し、海流の流速の変化を検出することを高感度で行なう必要がある。なお、ここでは、海流とは、海面の海水の流れを意味するものとする。反射信号の検出を高感度で行なうためには、長時間の信号を集積することが一般的であるが、長時間の信号を集積することは、津波の到来前に十分な対処する時間を得ることの必要性から制限される。
【０００６】
この発明は、上記のような問題を解決するためのもので、海面により反射する電波を分析し、海流の流速の変化を高感度で検出し、津波を高感度で検出できる海洋レーダ装置を得ることを目的とする。
【課題を解決するための手段】
【０００７】
この発明に係る海洋レーダ装置は、上記のような課題を解決するため、観測する海面を区切るセル毎に海面からの反射電波を分析し、前記反射電波のドップラ分布を前記セル毎に算出する信号処理装置と、地震の発生位置および規模から予測された津波の流速の海面上の分布をもとに、前記津波の流速が所定の流速から一定の範囲内である、前記津波の前縁部に位置する複数の前記セルを選択し、この選択した複数の前記セルからなる積分区間を設定する積分区間設定手段と、前記信号処理装置が算出したドップラ分布を前記積分区間について積分し、積分結果を算出する積分手段と、前記積分結果から、前記積分区間で振幅が最大となるドップラ周波数を求め、このドップラ周波数の変化を検出する検出手段とを備えるようにしたものである。
【発明の効果】
【０００８】
この発明に係る海洋レーダ装置は、津波の前縁部の、津波の流速が所定の流速から一定の範囲内である複数のセルを積分区間として海水の流速を測定し分析したドップラ分布データを積分するため、各セルのドップラ成分の振幅が最大となるドップラ周波数は積分する範囲で積分することにより積み上がり強調される。このため、積分により積み上がったドップラ周波数の変化を検出することで、海水の流速の変化を高感度で検出することができる。
【図面の簡単な説明】
【０００９】
【図１】本発明の海洋レーダ装置の検出する津波の発生、伝搬の概念図である。
【図２】本発明の実施の形態１に係る海洋レーダ装置の一例を表す系統図である。
【図３】本発明の実施の形態１に係る海洋レーダ装置の信号処理装置の出力を表す図である。
【図４】本発明の実施の形態１に係る海洋レーダ装置の海流によるドップラ分析の概念を示す図である。
【図５】本発明の実施の形態１に係る海洋レーダ装置の津波検出の概念を表す図である。
【図６】本発明の実施の形態１に係る海洋レーダ装置の津波検出における積分区間の概念を表す図である。
【図７】本発明の実施の形態１に係る海洋レーダ装置の津波検出における積分効果の概念を表す図である。
【図８】本発明の実施の形態１に係る海洋レーダ装置の処理の流れを表す図である。
【図９】本発明の実施の形態２に係る海洋レーダ装置の一例を表す系統図である。
【図１０】本発明の実施の形態３における海洋レーダ装置の分析するドップラ分布の概念図である。
【図１１】本発明の実施の形態３に係る海洋レーダ装置の一例を示す系統図である。
【図１２】本発明の実施の形態３に係る海洋レーダ装置のドップラ周波数の補正の概念を示す図である。
【図１３】本発明の実施の形態３に係る海洋レーダ装置の一例を示す系統図である。
【図１４】本発明の実施の形態３に係る海洋レーダ装置の一例を示す系統図である。
【図１５】本発明の実施の形態４に係る海洋レーダ装置の一例を示す系統図である。
【図１６】本発明の実施の形態４に係る海洋レーダ装置のドップラ変動周期分離手段の構成の一例を示す系統図である。
【図１７】本発明の実施の形態４に係る海洋レーダ装置のドップラ変動周期分離手段の構成の一例を示す系統図である。
【図１８】本発明の実施の形態４に係る海洋レーダ装置の一例を示す系統図である。
【図１９】本発明の実施の形態４に係る海洋レーダ装置の一例を示す系統図である。
【図２０】本発明の実施の形態４に係る海洋レーダ装置の一例を示す系統図である。
【発明を実施するための形態】
【００１０】
実施の形態１．
津波は、図１に示すように、地震による短時間での広範に及ぶ地殻変動により震源地の上に存在する海水が海面隆起となることで発生し、海上を伝播するものである。海洋レーダ装置１００は、図１に示す防災対象に対し、海洋から伝わる津波に対する対策を行なうため、防災対象の付近の場所に固定して設置し、海洋に対して電波を送信し、反射波を受信、分析することにより防災対象に向けて到来する津波を検出する。
【００１１】
図２は、本発明の実施の形態１に係る海洋レーダ装置の一例を示す系統図である。図２に示す海洋レーダ装置１００において、送信機１は津波を検出しようとする海洋（図２には記載していない）に向けてＦＭ変調を行ったＣＷ（Ｃｏｎｔｉｎｕｏｕｓ  Ｗａｖｅ）信号を送信し、反射した電波を受信機２で受信する。受信機２は、複数の受信アンテナ（図２には記載していない）で反射した電波を受信し、それぞれ送信機１の送信するＦＭ変調を行ったＣＷ（Ｃｏｎｔｉｎｕｏｕｓ  Ｗａｖｅ）信号と混合し、複数の受信アンテナに対応した複数のチャネルのビート信号を生成する。このようにして生成した複数のチャネルのビート信号は、それぞれが信号の反射点とレーダとの間の距離に対応した周波数の信号となる。受信機２では、それぞれのビート信号にＡ／Ｄ変換を行ない、複数のチャネルのデジタル信号を信号処理装置３に出力する。
【００１２】
信号処理装置３は、複数のチャネルのデジタル信号を入力して、方位分離手段６では、ＤＢＦ（デジタルビームホーミング）により方位ごとに合成して分離し、さらに、距離分離手段７では、ＦＦＴを行い、信号を周波数すなわち距離により分離する。このようにして、受信信号を、方位分離手段６のＤＢＦおよび、距離分離手段７のＦＦＴの分解能に応じた、方位および距離のセル毎のデータに分離する。ドップラ分離手段８では、方位および距離のセル毎のデータに分離したデータを複数回収集し、ＦＦＴを行うことで、海流の流速により生じたドップラ成分に分離する。なお、距離と方位で区切られたセルでは、セルの中でもその場所により海流の速度に微妙な違いがある。このセルの中の海流の速度の違いのため、それぞれのセルで観測される反射波のドップラ周波数は、１つの周波数に揃わない。それぞれのセルの中では、海流がそれぞれの速度で流れる面積が、速度に対して分布することに対応して、受信信号がそれぞれのドップラ周波数を持つ信号成分の振幅は、ドップラ周波数に対し分布を持つ。この、ドップラ周波数に対する振幅の分布を、以後ドップラ分布とよぶ。また、それぞれのセルにおけるドップラ周波数の観測値は、ドップラ分布で振幅が最大（ピーク）となる周波数で定義される。また、それぞれのセルにおける流速の観測値は、レーダ電波の反射波に対して上記のドップラ周波数の観測値をもたらす速度で定義される。信号処理装置３は、距離と方位で区切られたセル毎の海流のドップラ分布となったデータを検出装置４に出力し、検出装置４で津波の検出を行なう。
【００１３】
津波予測装置５は、検出装置４で津波を検出する際に使用する津波の流速予測データを作成する。津波予測装置５にある津波予測処理９は、地震発生位置および地震規模を入力すると、津波の伝搬計算等による津波シミュレータ１０によるシミュレーションを実施するかまたは、予め予想される地震発生地点についてシミュレーションを行なった結果を纏めた津波データベース１１を検索するなどにより津波の現れる位置の海流の流速の予測を求める。津波の到来方向、速度および形状は事前に津波シミュレーション等を使用して算出し、その結果を震源の位置、震度によって複数のケースに分けて記憶しておく事が出来る。そのため、津波データベース１１では、事前に複数のケースの津波の到来方向、速度および形状を算出して準備しておき、地震発生時に震源、震度等の情報が得られた際、この各ケースの中で震源、震度等の条件が最も近しいものを引用することが出来、各時刻における波面の位置を予測することができる。津波予測処理９では、このようにして津波の現れる位置の海流の流速の予測を求め、求めた海流の流速を、信号処理装置３がドップラ分布を出力する距離と方位で区切られたセル毎のデータに変換し、検出装置４に出力する。検出装置４では、信号処理装置３から入力するドップラ分布から津波を検出する際に、津波予測装置５から入力する津波の現れる位置の海流の流速の予測を使用する。
【００１４】
ここで、本発明の実施の形態１に係る海洋レーダ装置の検出する対象となる津波の特徴とその検出原理を示す。図３は、信号処理装置３の出力するドップラ分布の概念を示したものである。図３（ａ）の扇形の部分全体は、測定地点から測定する海面を表す。信号処理装置３は、図３（ａ）に示すように、海面を海洋レーダ装置の探知可能な、測定地点からの距離、方位で分割したセルに分割し、セル毎にドップラ分布を算出して出力する。このため、１つ１つのセルは、図３（ｂ）に示すようなドップラ分布を持つ。ところで、海面に送信した電波は、ブラッグ散乱により、（１）で示すブラッグ散乱周波数ｆＢに対して、±ｆＢで示すドップラ周波数を持つ成分が反射され、観測される。
【００１５】
【数１】
2014098686.tif 000003
【００１６】
但し、ｇは重力加速度、λは、海洋レーダ装置の使用している電波の波長を示す。このため、レーダ装置で観測されるドップラ分布は、１つ１つのセルについて、図３（ｂ）のように、±ｆＢ付近にピーク（振幅の最大値）を持つ分布となる。更に、海面による反射には、海流による影響が加わり、海洋レーダ装置で観測されるドップラ分布は、海流の流速による（２）に示すドップラ周波数ｆｓによりシフトされるため、観測される電波は、±ｆＢ＋ｆｓにそれぞれピークを持つ。
【００１７】
【数２】
2014098686.tif 000004
【００１８】
ここで、ｖは、海面の流速を示す。
【００１９】
次に図４および図５により、海面の流速の変化を分析することの概念を示す。図４は、海洋レーダ装置の測定する範囲に津波の波面が無い場合の流速の変化の分析の例を表す。図４（ａ）は、図３（ａ）と同様に、海洋レーダ装置の測定する海面を表す。図４（ａ）に示すように、レーダの信号処理装置３は、海面を海洋レーダ装置の探知可能な、測定地点からの距離、方位で分割したセルに分割し、セル毎にドップラ分布を算出して出力する。図４（ｂ）は、測定地点方向に進行する津波を検出するため、信号処理装置３の出力するドップラ分布を距離方向のセルについて並べ、図４（ｃ）は更に、図４（ｂ）を距離とドップラの関係をコンタープロット（等高線）で表すことでその分布量の多い部分を示し、距離に対するドップラ分布の変化を表したものである。図４（ａ）の様に津波が存在しない場合は、海流の流速によるドップラ分布は急激な変化を示さない。このため、信号処理装置３の出力するドップラ分布を距離方向のセルについて並べた結果は、図４（ｂ）および図４（ｃ）で示すとおり、ドップラ分布のピークの、振幅が最大となるドップラ周波数が直線状に並ぶ。
【００２０】
津波が存在する場合、津波の前縁部では、海流の流速が急激に変化する。図５（ａ）は、図４（ａ）と同様に、海洋レーダ装置の測定する海面を表す。図５（ａ）に示すように、レーダの信号処理装置３は、海面を海洋レーダ装置の探知可能な、測定地点からの距離、方位で分割したセルに分割し、セル毎にドップラ分布を算出して出力する。図５（ｂ）は、測定地点方向に進行する津波を検出するため、信号処理装置３の出力するドップラ分布を距離方向のセルについて並べ、図５（ｃ）は更に、図５（ｂ）を距離とドップラの関係をコンタープロット（等高線）で表すことでその分布量の多い部分を示し、距離に対するドップラ分布の変化を表したものである。図５（ａ）の様に、距離方向のセルの並びを津波の前縁部が通過する場合は、図５（ｂ）および図５（ｃ）では、Ａで示したとおり、津波の前縁部の通過する部分について、流速の急激な変化に従いドップラ分布のピークの、ドップラ分布の振幅が最大となるドップラ周波数が距離と共に激しく変化する。図４（ｂ）と図５（ｂ）、図４（ｃ）と図５（ｃ）を比較することにより明らかな通り、津波が到来すると、図５（ｂ）および図５（ｃ）の様に津波の前縁部でドップラ分布のピークが激しく変化するため、この変化を検出することにより、津波の検出を行なう。
【００２１】
ところで、津波の前縁部でドップラ分布が激しく変化することを検出するためには、ドップラ分布のピークの、ドップラ分布の振幅が最大となるドップラ周波数を２箇所以上で検出し、そのピーク同士が異なることを検出する必要がある。ドップラ分布のピークの検出を高感度で行なうためには、検出に使用する信号をより多く集積するため、広い面積からの反射信号を集積するなどが必要となる。広い面積からの反射信号を集めることは、多数のセルを集めて積分区間とし、積分区間内のセルの信号を積分することに相当する。ところで、流速が激しく変化する津波の前縁部においては、隣接する範囲のセルであっても、流速が大きく異なるセルも流速がほぼ等しいセルも存在する。積分区間内に流速が大きく異なるセルを多数含んでしまうと、積分することにより流速が相殺し合い、振幅のピークを強調する効果が得られない。このため、ドップラ分布のピークの検出を高感度で行なうためには、流速がほぼ等しいと推定されるセルを集めて積分区間を構成し、そのようにして構成した積分区間について積分した結果からピークの検出を行なう必要がある。
【００２２】
図６は、津波による海面のドップラ分布の変化を検出する原理を示したものである。図６（ａ）に示すように、津波の進行方向にほぼ垂直方向に波面が形成されるが、前述の通り、津波の前縁部では海面の流速が激しく変化する。このため、津波の進行方向に沿う方向に並んだセルについて比較すると観測されるドップラ周波数の変化が激しいが、波面に沿う方向に並んだセルについて比較すると海流の変化は少なく、観測されるドップラ周波数は一定している。
【００２３】
ここで、図６（ａ）の様に、波面に平行な方向に並んだセルを集めて第１の積分区間３１および、第２の積分区間３２を設けると、波面に平行な方向については、海流の変化は少ないため、第１の積分区間３１を構成するセルで観測されるドップラ分布は互いに同じ周波数にピークを持つ。また、同様に、第２の積分区間３２を構成するセルで観測されるドップラ分布についても、互いに同じ周波数にピークを持つ。また、津波の進行方向に沿う方向では海面の流速が激しく変化するため、第１の積分区間３１を構成するセルで観測されるドップラ分布と第２の積分区間３２を構成するセルで観測されるドップラ分布を比較すると、図６（ｂ）の様に、異なる位置にピークを持つ。
【００２４】
このことから、図７に示すように、それぞれが同じ位置にドップラ分布のピークを持つセルで構成された第１の積分区間３１でドップラ分布を積分することにより、ドップラのピークは積み上がり強調される。同様に、それぞれが同じ位置にドップラ分布のピークを持つセルで構成された第２の積分区間３２でドップラ分布を積分することにより、ドップラのピークは積み上がり強調される。それぞれの積分区間を構成するセルのピークが異なるため、第１の積分区間３１と第２の積分区間３２では、積み上がり強調されるピークの位置は、異なるものとなる。このため、観測されるドップラ分布について、第１の積分区間３１で積分を行なった結果と、第２の積分区間３２で積分を行なった結果を比較することで、海流の変化を高感度で検出することができる。
【００２５】
上記原理を踏まえ、図２で示される検出装置４の動作を説明する。検出装置４は信号処理装置３からドップラ分布を入力し、メモリ１２に記憶する。他方、検出装置４では、積分区間設定手段１３により、津波予測装置５の出力を基に波面に平行な方向に並んだセルを選択し、第１の積分区間３１および、第２の積分区間３２として設定する。積分区間設定手段１３は、信号処理装置３がドップラ分布を出力する距離と方位で区切られたセル毎のデータとして津波予測装置５から海流の流速の予測値を入力する。入力した流速予測値をセル毎に相関を行い、同じ積分区間内で流速が一定の範囲内であって、流速が同じと予測されるセル同士を集めることで、それぞれ津波の波面に沿う方向に広がりを持つ積分区間とする。その積分区間を、流速が同じと予想されるセル同士を集めた第１の積分区間３１と、流速が、第１の積分区間とは異なる流速と同じと予想されるセル同士を集めた第２の積分区間３２とを設定する。
【００２６】
例として、海流の流速の極大値に対して、流速が極大値から一定の範囲内であるセルを選択して集めて第１の積分区間３１を設定し、海流の流速の極小値に対して、流速が極小値から一定の範囲内であるセルを選択して集めて第２の積分区間３２を設定するなど、ドップラ分布で違いが判別できる２つの流速について、セル内の流速がそれぞれの流速になると予想される２つの積分区間を設定する。ここで、セルを集めて積分区間を設定する際に規準とする流速の差についての「一定の範囲」としては、例えば、信号処理装置３のドップラ分離手段８の速度分解能や、受信機２と信号処理装置３の周波数精度を含めた海流の速度の分解能、検出装置４で津波を検出した結果から類推できる値などを設定する。
【００２７】
積分手段（空間積分手段１４）は、距離、方位でセルに分割した各セルのドップラ分布をメモリ１２より入力し、入力したドップラ分布を第１の積分区間３１および第２の積分区間３２それぞれの積分区間について積分し、その積分結果を検出手段１５に出力する。検出手段１５は、ドップラ分布を第１の積分区間３１および第２の積分区間３２それぞれの積分区間について積分した積分結果からそれぞれドップラ分布のピークを検出し、それぞれのピークを比較する。このとき、積分区間設定手段１３の設定した第１の積分区間３１および第２の積分区間３２を津波の前縁が通過すると、それぞれの積分区間におけるドップラ分布のピークが積み上がり強調され、その変化を高感度で検出することができる。
【００２８】
尚、本発明の実施の形態１に係る発明の構成では、積分手段（空間積分手段１４）と検出手段１５を第１の検出系とすると、それ以外にも、第２の検出系として、検出手段１６でも津波の検出を行い、更に、第３の検出系として、時間積分手段１７と検出手段１８でも津波の検出を行なう。第２の検出系に係る検出手段１６は、距離、方位でセルに分割した各セルのドップラ分布をメモリ１２より入力し、入力したドップラ分布のピークの、最大となるドップラ周波数をそれぞれのセル間で比較し、海流の流速の変化を検出する。
【００２９】
また、第３の検出系では、時間積分手段１７は、距離、方位でセルに分割した各セルのドップラ分布をメモリ１２より入力し、入力したドップラ分布をそれぞれのセルについてそれぞれ一定時間蓄積し、積分することでドップラ分布のピークを強調し、検出手段１８でそれぞれのセル間で比較し、海流の流速の変化を検出する。検出手段１５、検出手段１６、検出手段１８のそれぞれで検出された津波は、検出結果統合手段１９で同一のものについては１つにまとめ、メモリ２０に記憶する。メモリ２０に記憶した津波検出内容は、図２には表示されていないが、表示や警報の発信など各種の処理に使用される。このように、第１ないし第３の検出系を用いることで、それぞれの検出系の長所を組合せ、広範囲の津波を高感度で、高い信頼度で検出することができる。
【００３０】
図８は、本発明の実施の形態１に係る海洋レーダ装置の動作の一例を表す。最初に、津波を検出する対象となる海域に対して、送信機１により電波を送信し、その反射波を受信機２で受信することで測定を行なう（ＳＴ８０１）。次に、受信機２の受信信号を信号処理装置３で分析し、方位および距離のセル毎のデータに分離し、さらに、各セル毎のデータについて、ＦＦＴにより、流速により生じたドップラ成分の分布を算出する（ＳＴ８０２）。ＳＴ８０２により算出した各セル毎のドップラ成分の分布をＳＴ８２１以下による検出、ＳＴ８３１以下による検出、ＳＴ８４１以下による検出の３通りの方法で検出を行なう。
【００３１】
ＳＴ８２１以下では、セル毎のドップラ成分の分布の空間積分を使用するが、積分区間として、津波予測装置５による予測結果を使用する。地震が発生すると、地震発生位置、規模などのデータが津波予測装置に入力される（ＳＴ８１１）。津波予測装置５は、入力されたデータを基に津波の流速分布を求めるが、このとき、津波シミュレータ１０によりリアルタイムのシミュレーションを行なうかまたは、津波データベース１１に事前に算出して蓄積した地震発生地点と津波到達時間のデータを津波データベース１１から最も条件が近しいケースを求め、シミュレーション結果とデータベースの検索結果のいずれかまたは両方の結果を、信号処理装置３がドップラ分布を出力する距離と方位で区切られたセル毎のデータに変換し、検出装置４に出力する（ＳＴ８１２）。
【００３２】
津波予測装置５から津波の流速分布を入力した検出装置４では、積分区間設定手段１３が、例えば流速の極大値と極小値など、異なる２つの流速について、流速が同じと予測される方位、距離範囲を纏めてそれぞれ第１の積分区間３１、第２の積分区間３２とし、空間積分手段１４に出力する（ＳＴ８２１）。空間積分手段１４では、ＳＴ８０２で算出したセル毎のドップラ分布の内、第１の積分区間３１を構成するセルのドップラ分布について、ドップラ成分毎に足し合わせることで積分を行なう。同様に、空間積分手段１４では、ＳＴ８０２で算出したセル毎のドップラ分布の内、第２の積分区間３２を構成するセルについてもドップラ成分毎に足し合わせることで積分を行なう。このようにして、第１の積分区間３１、第２の積分区間３２それぞれについてのドップラ分布の積分を行ない、検出手段１５に出力する（ＳＴ８２２）。
【００３３】
検出手段１５では、ＳＴ８２２で算出した第１の積分区間３１、第２の積分区間３２それぞれについて、積分区間毎に、ドップラ分布の積分結果を、積分区間について設定したスレショルドと比較し、検出された信号のピークを求める。検出に使用するスレショルドについては、積分結果について、前述の（１）式で示されるブラッグ散乱周波数ｆＢについて、±ｆＢの前後のドップラ成分をレファレンス信号とし、ＣＦＡＲ（Ｃｏｎｓｔａｎｔ  Ｆａｌｓｅ  Ａｌａｒｍ  Ｒａｔｅ）の原理により算出する。さらに、求めたピークの位置すなわち海流の速度を第１の積分区間３１と第２の積分区間３２とで比較し、その比較結果により急激な流速変化が有れば、津波を検出とする（ＳＴ８２３）。
【００３４】
なお、地震の震源の位置、震度、伝搬条件などにより、複数の津波が予想される場合がある。津波予測装置５は、複数の津波が予想される場合、複数の予測をそれぞれ検出装置４に出力する。検出装置４の積分区間設定手段１３では、津波予測装置５の出力する複数の予測それぞれに対して第１の積分区間３１および第２の積分区間３２を算出する。また、空間積分手段１４は、複数の予測に対して、それぞれ第１の積分区間３１および第２の積分区間３２についての積分を行い、検出手段１５は、複数の予測に対して、それぞれ第１の積分区間３１の積分結果と第２の積分区間３２の積分結果の比較を行い、津波を検出する。
【００３５】
ＳＴ８２１以下では、セル毎のドップラ成分の分布の空間積分を使用するが、それと並行して、ＳＴ８３１では、セル毎のドップラ分布を積分せずに津波の検出を行なう。検出手段１６では、信号処理装置３の算出したセル毎のドップラ成分の分布について、セル毎にドップラ分布のピークの、振幅が最大になるドップラ周波数を求め、さらに、そのピークの周波数すなわち海流の速度を隣接するセル同士で比較し、その比較結果により急激な流速変化が有れば、津波を検出とする（ＳＴ８３１）。
【００３６】
さらに、津波の前縁部は数ｋｍの広がりを持ち、津波の周期は数分から長い場合には数時間であることから、津波に係る流速データは比較的長い時間に渡って安定したものであり、このことから、ＳＴ８４１以下では、セル毎に得られるドップラ分布を時間方向に積分し、観測したデータを積み上げることにより、感度の向上を図る。時間積分手段１７では、信号処理装置３の算出したセル毎のドップラ成分の分布について、セル毎に所定の時間、ドップラ成分毎に足し合わせることで積分を行なう（ＳＴ８４１）。なお、所定の時間は、津波の想定される周期より十分短い時間で設定される。時間積分手段１７の出力するセル毎の積分結果について、検出手段１８では、セル毎にドップラ分布のピークを求め、さらに、そのピークの位置すなわち海流の速度を隣接するセル同士で比較し、その比較結果により急激な流速変化が有れば、津波を検出とする（ＳＴ８４２）。
【００３７】
検出結果統合手段１９では、検出手段１５、検出手段１６、検出手段１８それぞれに検出した津波について、同一の検出結果を１つに纏める（ＳＴ８５１）。さらに、連続して津波を探知する、一定時間ごとに津波を探知する、または、緊急時に津波を探知するなどの運用に従い、ＳＴ８０１からＳＴ８５１までを繰り返す（ＳＴ８５２）。
【００３８】
なお、上記実施の形態１では、空間積分手段１４で積分を行なう積分区間については、第１の積分区間３１と第２の積分区間３２の２つの積分区間を使用する例を示したが、海流の速度が異なる３つまたはそれ以上の積分区間を設け、その３つ以上の積分区間について積分した結果を比較しても同様の効果を得ることができる。さらに、積分区間を１つのみ設定し、その積分区間のドップラ分布を時間と共に観測し、ドップラ分布のピークの、振幅が最大となるドップラ周波数が変化することから津波を検出することも可能である。また、検出手段１６および検出手段１８についても同様に、隣接するセル同士で比較するかわりに、同一セルについてドップラ分布のピークの、振幅が最大となるドップラ周波数を時間と共に観測し、ドップラ分布のピークの、振幅が最大となるドップラ周波数が変化することから津波を検出することも可能である。
【００３９】
また、上記実施の形態１では、津波予測装置５の出力する海流の流速の予測については津波予測処理９で流速の予測を、信号処理装置３がドップラ分布を出力する距離と方位で区切られたセル毎のデータに変換した上で検出装置４に出力する例を示したが、津波予測処理９は津波シミュレータ１０および津波データベース１１固有の座標で算出した海流の流速の予測を検出装置４に出力し、検出装置４では積分区間設定手段１３などで、流速の予測を、信号処理装置３がドップラ分布を出力する距離と方位で区切られたセル毎のデータに変換する構成としても良い。また、津波予測処理９は、津波の流速の予測を、任意の座標系で津波の流速の海面上の分布として検出装置４に出力し、検出装置４では積分区間設定手段１３などで、流速の予測を、信号処理装置３がドップラ分布を出力する距離と方位で区切られたセル毎のデータに変換する構成としても良い。さらに、上記実施の形態１における信号処理装置がドップラ分布を算出するセルについても、海面を海洋レーダ装置の探知可能な、測定地点からの距離、方位で分割したセルで説明したが、セルは、海面を細かく分割するものであれば、その分割の方法が津波の検出に影響を与えるものではなく、任意の方法で海面を分割したセルにより処理を行ってもかまわない。
【００４０】
以上により、本発明の第１の実施の形態に係る海洋レーダ装置では、津波の前縁部に、海水の流速が同じであるセルを集めた積分区間を１つまたは複数設定し、設定した積分区間についてドップラ分布を積分するため、積分する範囲のドップラ分布のピークを積分することで強調される。このため、強調されたドップラ成分のピークの周波数の変化を検出することで、海流の速度の変化を高感度で検出することができる。また、セル毎のドップラ分布のピークの周波数を直接比較することによる検出と組み合わせること、および、さらに、セル毎のドップラ分布を時間積分し、積分結果を比較することによる検出と組み合わせることにより、広範囲の津波を高感度で、高い信頼度で検出することができる。
【００４１】
なお、上記実施の形態１では、送信機１はＦＭ変調を行ったＣＷ信号を送信し、受信機２は、受信した電波それぞれに送信機１の送信するＦＭ変調を行ったＣＷ信号と混合し、信号の反射した距離と対応する周波数の信号に変換し、信号処理装置３は、送信信号と混合した受信信号の周波数を分離することで受信した信号の距離による分離を行なうＦＭ−ＣＷ（または、ＦＭ−ＩＣＷ）方式を使用しているが、レーダとしては、パルスレーダを使用しても良い。パルスレーダを使用する場合、送信機１はパルス信号を送信し、受信機２は反射信号を受信し、信号処理装置３では、距離分離手段７において、周波数の替りに反射信号の反射時間により距離分離を行なう以外は、同様である。
【００４２】
実施の形態２．
図９は、本発明の実施の形態２に係る海洋レーダ装置の一例を示す系統図である。実施の形態１は、１台の受信機により海面の流速を測定し、津波を検出する構成だが、実施の形態２では、異なる場所に設けた２台の受信機により受信したデータを分析し、合成する。この構成の違いにより、実施の形態１では、レーダの距離方向の海流の流速の成分によるドップラを測定し、実施の形態２では、２つのレーダの距離方向を合成した２次元の海流の流速の成分によるドップラを測定する。このため、実施の形態１の構成では、レーダの距離方向の流速成分の大きい津波は良く検出するが、例えば、レーダの距離方向と垂直方向に伝搬する津波などレーダの距離方向の流速成分が小さい津波については、検出能力が制限される。これに対して、実施の形態２の構成では、津波を検出するに際して、方向による制限が少なく、全方向に対する検出能力が向上する。図９において、送信機１は、図２と同様であるため、説明を省略する。また、津波予測装置５ａは、津波予測処理９ａが、津波の現れる位置の海流の流速の予測を求め、検出装置４ａに出力する際に、求めた海流の流速を、信号処理装置３ａが受信機２ａと受信機２ｂに対応してドップラ分布を出力する距離と方位で区切られた２種類のセル毎のデータに変換する以外は図２の津波予測装置５と同様である。
【００４３】
図９において、受信機２ａおよび、受信機２ｂは、それぞれ図２における受信機２と同等の構成であり、それぞれを異なる場所に設置されている。このため、受信機２ａ、受信機２ｂで受信した信号は、海面の同一の地点については、別の角度から測定した信号となる。信号処理装置３ａは、図２における信号処理装置３を受信機２ａの信号を処理する部分、受信機２ｂの信号を処理する部分の２系統にしたものである。
【００４４】
受信機２ａは、複数の受信アンテナ（図９には記載していない）で受信した電波それぞれに送信機１の送信するＦＭ変調を行ったＣＷ（Ｃｏｎｔｉｎｕｏｕｓ  Ｗａｖｅ）信号と混合し、ビート信号を生成する。このようにして生成したそれぞれのビート信号は、信号の反射点とレーダとの間の距離に対応した周波数の信号となる。受信機２ａでは、それぞれのビート信号にＡ／Ｄ変換を行ない、複数のチャネルのデジタル信号を信号処理装置３ａに出力する。
【００４５】
受信機２ｂは、複数の受信アンテナ（図９には記載していない）で受信した電波それぞれに送信機１の送信するＦＭ変調を行ったＣＷ（Ｃｏｎｔｉｎｕｏｕｓ  Ｗａｖｅ）信号と混合し、ビート信号を生成する。このようにして生成したそれぞれのビート信号は、信号の反射点とレーダとの間の距離に対応した周波数の信号となる。受信機２ｂでは、それぞれのビート信号にＡ／Ｄ変換を行ない、複数のチャネルのデジタル信号を信号処理装置３ａに出力する。
【００４６】
受信機２ａの受信した信号は、信号処理装置３ａの方位分離手段６ａで、ＤＢＦ（デジタルビームホーミング）により方位ごとに分離し、さらに、距離分離手段７ａでＦＦＴを行い、信号を周波数すなわち距離により分離する。このようにして、受信信号を、方位分離手段６ａのＤＢＦおよび、距離分離手段７ａのＦＦＴの分解能に応じた大きさの、方位および距離毎のセル毎のデータに分離した後、ドップラ分離手段８ａでは、方位および距離毎のセル毎のデータに分離したデータを複数回収集し、ＦＦＴを行うことで、海流の流速により生じたドップラ成分に分離し、その結果を検出装置４ａに出力する。
【００４７】
受信機２ｂの受信した信号は、信号処理装置３ａの方位分離手段６ｂで、ＤＢＦ（デジタルビームホーミング）により方位ごとに分離し、さらに、距離分離手段７ｂでＦＦＴを行い、信号を周波数すなわち距離により分離する。このようにして、受信信号を、方位分離手段６ｂのＤＢＦおよび、距離分離手段７ｂのＦＦＴの分解能に応じた大きさの、方位および距離毎のセル毎のデータに分離した後、ドップラ分離手段８ｂでは、方位および距離毎のセル毎のデータに分離したデータを複数回収集し、ＦＦＴを行うことで、海流の流速により生じたドップラ成分に分離し、その結果を検出装置４ａに出力する。
【００４８】
検出装置４ａは信号処理装置３ａから受信機２ａおよび受信機２ｂのそれぞれで受信した信号のドップラ分布を入力し、メモリ１２ａに記憶する。他方、検出装置４ａでは、津波予測装置５ａから受信機２ａおよび受信機２ｂのそれぞれの受信信号に対応するセル毎の海流の流速予測値を基に、積分区間設定手段１３ａが、検出装置４と同様に、波面に沿った方向に並んだセルを算出し、その範囲のセルを集めて、第１の積分区間３１および、第２の積分区間３２を受信機２ａおよび受信機２ｂのそれぞれの受信信号に対して設定する。
【００４９】
積分手段（空間積分手段１４ａ）は、距離、方位でセルに分割した各セルの受信機２ａおよび受信機２ｂのそれぞれの受信信号についてのドップラ分布をメモリ１２ａより入力し、入力したドップラ分布を第１の積分区間３１および第２の積分区間３２それぞれの積分区間について積分し、その積分結果を受信機２ａおよび受信機２ｂのそれぞれで受信した信号毎に検出手段１５ａに出力する。検出手段１５ａは、受信機２ａおよび受信機２ｂのそれぞれで受信した信号について、ドップラ分布を第１の積分区間３１および第２の積分区間３２それぞれの積分区間について積分した積分結果からそれぞれドップラ分布のピークの、振幅が最大となるドップラ周波数を検出し、それぞれのピークのドップラ周波数を比較し、津波を検出する。検出した津波は、同じセルについて受信機２ａと受信機２ｂで受信したそれぞれの速度を合成し、２次元の速度成分を持った津波のデータとする。
【００５０】
検出手段１６ａは、距離、方位でセルに分割した、受信機２ａおよび受信機２ｂのそれぞれで受信した信号の各セルのドップラ分布をメモリ１２ａより入力し、入力したドップラ分布をそれぞれのセル間で比較し、受信機２ａおよび受信機２ｂのそれぞれで受信した信号毎に海流の流速の変化を検出する。検出した津波は、同じセルについて受信機２ａと受信機２ｂで受信したそれぞれの速度を合成し、２次元の速度成分を持った津波のデータとする。
【００５１】
時間積分手段１７ａは、距離、方位でセルに分割した受信機２ａおよび受信機２ｂのそれぞれで受信した信号の各セルのドップラ分布をメモリ１２ａより入力し、入力したドップラ分布をそれぞれのセルについて受信機２ａおよび受信機２ｂのそれぞれの受信信号毎に、それぞれ一定時間蓄積し、積分することでドップラ分布のピークを強調し、検出手段１８ａで受信機２ａおよび受信機２ｂのそれぞれの受信信号毎に、それぞれのセル間で比較し、海流の流速の変化を受信機２ａおよび受信機２ｂのそれぞれの受信信号毎に検出する。検出した津波は、同じセルについて受信機２ａと受信機２ｂで受信したそれぞれの速度を合成し、２次元の速度成分を持った津波のデータとする。
【００５２】
検出手段１５ａ、検出手段１６ａ、検出手段１８ａのそれぞれで検出された津波は、検出結果統合手段１９ａで同一のものについては１つにまとめ、メモリ２０ａに記憶する。メモリ２０ａに記憶した津波検出内容は、図９には表示されていないが、表示や警報の発信など各種の処理に使用される。
【００５３】
なお、図９では、全方向に向かう津波を検出するため、受信機２ａおよび受信機２ｂの２台の受信機を用いる構成だが、３台以上の受信機により受信し、その検出結果を合成しても良い。
【００５４】
また、図９では、１台の送信機１から送信した電波の海面による反射を受信機２ａおよび受信機２ｂで受信する構成としたが、それぞれの受信機に対して例えば送信機１ａと送信機１ｂなどを設けることにより、送信機１ａで送信した電波の海面による反射信号を受信機２ａで受信し、送信機１ｂで送信した電波の海面による反射信号を受信機２ｂで受信する構成としてもよい。
【００５５】
以上により、本発明の第２の実施の形態に係る海洋レーダ装置では、受信機２ａおよび受信機２ｂなどの異なる場所に設置された２台以上のレーダで受信した受信データから津波を検出し、検出結果を合成することにより、実施の形態１の効果に加え、津波を検出するに際して、方向による制限が改善され、全方向に対する検出能力が向上するという効果がある。
【００５６】
実施の形態３．
海面における海水の流れには、地震により発生する津波と、地震の有る無しにかかわらず発生している流れ（以下、通常海流と言う）がある。通常海流には、暖流、寒流、潮汐により発生する潮流、および風により発生する吹走流等が存在する。このため、レーダ装置により観測され、信号処理装置３で分析される海水の流れの速度は、津波による流れと上記の通常海流による流れが重ね合わせされたものとなる。津波による流速をｖｔ、通常海流の流速をｖｄとすると、海面の流速ｖは、（３）式のようになる。
【００５７】
【数３】
2014098686.tif 000005
【００５８】
海面における海水の流速と、レーダ装置により観測されるドップラ周波数については、前述の（２）式の関係が有るため、津波によるドップラ周波数の成分をｆｔ、通常海流によるドップラ周波数の成分をｆｄとすると、信号処理装置３で算出されるドップラ周波数ｆSとの関係は、（４）式のようになる。
【００５９】
【数４】
2014098686.tif 000006
【００６０】
図１０は、通常海流の流速の合計ｖｄによるドップラ周波数ｆｄの影響を示すものである。ｆｄ＝０であれば、海洋レーダ装置で観測されるドップラ分布のピーク±ｆＢ＋ｆｓは、±ｆＢ＋ｆｔとなり、ブラッグ散乱周波数±ｆＢの周囲で、津波による周波数ｆｔの変動が加わった周波数となる。一方、ｆｄ≠０の場合は、海洋レーダ装置で観測されるドップラ分布のピーク±ｆＢ＋ｆｓは、±ｆＢ＋ｆｄの周囲で、津波による周波数ｆｔの変動が加わった周波数となる。
【００６１】
実施の形態１および実施の形態２に係る海洋レーダ装置では、積分区間設定手段１３、１３ａは、津波による流速が一定と予想されるセルで積分区間を構成する。このため、津波以外の通常海流の流速の合計速度ｖｄが積分区間内の各セルで異なる場合は、（４）式のｆｄが積分区間内のセル毎に異なることとなり、レーダ装置で観測され、信号処理装置３で算出されるドップラ周波数ｆｓの分布もそれぞれのセルで異なるものとなる。このような場合は、レーダ装置で観測されるドップラ周波数ｆｓを積分区間内で積分しても、そのピークが積み上がるという効果を十分に得ることが困難になる場合がある。
【００６２】
実施の形態３に係る海洋レーダ装置では、積分区間設定手段１３の設定する積分区間内の各セルについて、通常海流の速度ｖｄによるドップラの影響ｆｄを予測し、レーダ装置で観測されるドップラ周波数ｆｓを津波以外の通常海流の流速によるドップラの影響の予測値だけ差し引き補正する。通常海流の流速ｖｄによるドップラ分布への影響ｆｄを差し引いた上で、空間積分手段１４がドップラ分布を積分することで、通常海流の合計速度ｖｄが積分区間内の各セルで異なる場合でも、そのピークが積み上がるという効果を得るものである。
【００６３】
図１１は、実施の形態３に係る海洋レーダ装置の構成を示す図である。また、図１２は、実施の形態３に係る海洋レーダ装置の、ドップラ周波数の補正の概念を示す図である。図１１において、海洋レーダ装置は、送信機１、受信機２、信号処理装置３、検出装置４ｂ、津波予測装置５、および通常流予測装置４１を備えている。送信機１、受信機２、信号処理装置３、津波予測装置５は、図２と同様である。また、検出装置４ｂの各構成要素についても、図２と同一符号のものは、同様の構成を表す。
【００６４】
通常流予測装置４１は、潮流予測処理４２、吹走流予測処理４３、暖流・寒流予測装置４４、通常流合成処理４５を有し、通常海流の速度ｖｄをセル毎に予測する。潮流予測処理４２は、潮汐表、気圧情報をもとに各セルについての潮流の速度を予測する。これは、潮流は、潮位変化などによる水圧により発生するためである。また、吹走流予測処理４３は、風速、吹続時間、吹走距離をもとに各セルについての吹走流の速度を予測する。これは、吹走流は、海上の風の応力により発生するためである。暖流、寒流予測処理４４は、海上保安庁の有する海流情報などの海流情報により暖流、寒流の速度を予測する。通常流合成処理４５は、成因の異なる潮流、吹走流、暖流、寒流の速度を合成し、津波以外の通常海流の速度ｖｄをセル毎に予測する。例えば、図１２（ａ）に示す第１の積分区間３１を構成する各セルをそれぞれセル１、セル２、・・・セルｎとすると、通常流合成処理４５では、それぞれのセル１、セル２、・・・セルｎでの津波以外の通常海流の流速ｖｄ１、ｖｄ２、・・・ｖｄｎを算出する。
【００６５】
検出装置４ｂでは、それぞれのセルで観測されるドップラ分布は、ドップラ分布シフト手段５１を通して空間積分手段１４に入力される。空間積分手段１４は、ドップラ分布シフト手段５１を通して入力されるドップラ分布を積分し、検出手段１５は、空間積分手段１４の積分結果から津波の検出を行なう。ドップラ分布シフト手段５１は、通常流予測装置４２の予測した速度ｖｄを基に津波以外の通常海流の流速によるドップラｆｄを算出し、メモリ１２上のドップラ分布に対してそれぞれ−ｆｄの周波数方向のシフトを行い、セル毎に速度ｖｄの影響を差し引く。例えば、図１２（ａ）の第１の積分区間３１を構成する各セルであるセル１、セル２、・・・セルｎについては、図１２（ｂ）で示すように、それぞれのセルで観測されるドップラ分布は、津波以外の津波以外の通常海流の流速の合計ｖｄ１、ｖｄ２、・・・ｖｄｎの影響により、それぞれ、ｆｄ１、ｆｄ２、・・・ｆｄｎだけ周波数方向にシフトされたものとなっている。
【００６６】
図１２（ｃ）は、ドップラ分布シフト手段５１が、それぞれのセルで観測されるドップラ分布を、−ｆｄ１、−ｆｄ２、・・・−ｆｄｎだけ周波数方向にシフトすることにより、津波以外の通常海流の影響を差し引いた結果を示したものである。セル１、セル２、・・・セルｎは、それぞれ、津波によるドップラ周波数ｆｔ１、ｆｔ２、・・・ｆｔｎの予測値が同じセルである。セル１、セル２、・・・セルｎでは、ドップラシフト手段５１により周波数方向にシフトし補正した結果、ドップラ周波数がほぼ同じ値になる。空間積分手段１４は、ドップラ分布シフト手段５１で周波数方向にシフトを行なったドップラ分布を積分するが、上記のように、補正し津波の流速によるほぼ同じ値となったドップラ周波数を積分するため、積分手段（空間積分手段１４）の算出する積分結果において、ドップラ分布のピークは積み上がり強調される。検出手段１５は、空間積分手段１４の出力である積分結果の変化により津波を検出する。空間積分手段１４の出力である積分結果では、ドップラ分布のピークが強調されているため、検出手段１５は、津波を高感度で検出することができる。
【００６７】
なお、図１１では、検出装置４ｂは、ドップラ分布シフト手段５１によりメモリ１２上のドップラ分布のシフトを行い、空間積分手段１４に出力し、空間積分手段１４により、ドップラ分布シフト手段５１からの入力を積分するが、図１３に示す構成の様に、検出装置４ｃは、ドップラ分布シフト手段５１によりメモリ１２上のドップラ分布のシフトを行ない、その結果をメモリ５２に出力し、メモリ５２上のドップラ分布を空間積分手段１４により積分する構成としても良い。
【００６８】
図１３の通常流予測装置４１は、図１１における通常流予測装置４１と同一の構成である。図１３のドップラ分布シフト手段５１は、図１１におけるドップラ分布シフト手段５１と同一の構成であり、ドップラ分布から津波以外の通常海流の影響を除いた結果を、直接空間積分手段１４ではなく、メモリ５２に出力する。直接空間積分手段１４は、ドップラ分布シフト手段５１からではなく、メモリ５２より入力される、津波以外の通常海流の影響を除いたドップラ分布を積分する。
【００６９】
さらに、検出手段１６による津波の検出、時間積分手段１７と検出手段１８による津波の検出も、メモリ１２からではなく、メモリ５２上のドップラ分布のシフトを行ったドップラ分布により津波の検出に係る処理を行う構成とする。このような構成とした場合、検出手段１６による津波の検出、時間積分手段１７と検出手段１８による津波の検出、および空間積分手段１４と検出手段１５による津波の検出ではそれぞれメモリ５２上の、ドップラ分布から津波以外の通常海流の影響を差し引いたドップラ分布により津波の検出に係る処理を行う。このため、メモリ１２上のドップラ分布から津波以外の通常海流の影響を除いたドップラ分布を使用するため、検出手段１６、検出手段１８における津波の検出においても、津波以外の潮流、吹走流などの流速のセル毎の違いや時間変動に影響されずに検出を行なうことができる。
【００７０】
さらに、通常流予測装置４１については、本発明に係る海洋レーダ装置あるいは別途設置された海洋レーダ装置等で計測された海面の流速のデータをデータベース化し、データ取得時の季節、気温、風速、気圧などの条件が一致するデータをデータベースより抽出することで、津波以外の通常海流の速度ｖｄをセル毎に予測しても良い。
【００７１】
加えて、実施の形態３に係る海洋レーダ装置の、ドップラ分布シフト手段５１により、ドップラ分布から、通常流予測装置４１で求めた通常海流の影響を差し引く構成は、図９のように、受信機２ａおよび受信機２ｂなどの異なる場所に設置された２台以上のレーダで受信した受信データから津波を検出し、検出結果を合成する構成に対しても適用することが可能である。図９では、受信機２ａで受信された電波は、信号処理装置３ａ中の方位分離手段６ａ、距離分離手段７ａ、速度分離手段８ａで処理され、空間積分手段１４ａで積分される。これらの処理は、図１４において、受信機２で受信された電波が、信号処理装置３中の方位分離手段６、距離分離手段７、速度分離手段８で処理され、空間積分手段１４で積分されるそれぞれの処理と全く同等である。また、図９における受信機２ｂで受信された電波についても同様である。このため、図９における受信機２ａ、受信機２ｂで受信された電波のそれぞれについても、図１４における、受信機２で受信された電波と全く同様に、ドップラ分布シフト手段５１により、ドップラ分布から、通常流予測装置４１で求めた通常海流の影響を除くことができる。
【００７２】
実施の形態４．
津波の周期は、短いもので数分から長いものは数時間である。それに対して、潮流の周期は、１２時間および２４時間程度であり、津波の周期とは異なる。実施の形態４では、ノイズとなる津波と津波以外の潮流の影響をドップラ周波数の変動の周期の違いにより除去する。このようにして、ノイズとなる潮流の影響が各セル毎に異なることによる、空間積分手段１４により積分への影響を除き、検出手段１８における津波の検出においても、潮流の流速のセル毎の違いや時間変動に影響されずに検出を行なうことができるようにするものである。
【００７３】
図１５は、実施の形態４に係る海洋レーダ装置および検出装置の構成を示す図である。図１５において、海洋レーダ装置は、送信機１、受信機２、信号処理装置３、検出装置４ｅ、および津波予測装置５を備えている。送信機１、受信機２、信号処理装置３、津波予測装置５は、図２と同様である。また、検出装置４ｅの各構成要素についても、図２と同一符号のものは、同様の構成を表す。
【００７４】
検出装置４ｅでは、図１１における通常流予測装置４１のかわりに、ドップラ変動周期分離手段５３が、メモリ１２上のドップラ分布を元に、潮流などによるドップラｆｄをセル毎に推定する。ドップラ分布シフト手段５１ｂは、ドップラ変動周期分離手段５３が推定した潮流などによるドップラｆｄを使用し、メモリ１２上のドップラ分布に対してそれぞれセル毎に−ｆｄの周波数シフトを行うことで、潮流などの長い周期で速度が変動する流れによる影響を除いたドップラ分布を算出する。空間積分手段１４は、ドップラシフト手段５１ｂが潮流などの影響を除いたドップラ分布を積分する。
【００７５】
図１６は、ドップラ変動周期分離手段５３の構成を示す。ドップラ変動周期分離手段５３は、図１６の様に、海流のドップラ算出手段６１、ドップラ変動解析手段６２、流速変動成分抽出手段６３、潮流のドップラ算出手段６４で構成され、メモリ１２上のドップラ分布を分析し、潮流の流速をセル毎に予測する。海流のドップラ算出手段６１は、セル毎に図１０のドップラ分布からその２つのピークｆＢ＋ｆｓ、−ｆＢ＋ｆｓを求め、その平均値を求めることで海流の流速によるドップラｆｓを算出する。
【００７６】
ドップラ変動解析手段６２は、海流のドップラ算出手段６１が算出した海流の流速によるドップラｆｓをセル毎に、例えば数日間などの十分長い時間蓄積し、ＦＦＴにより海流の流速の周期毎のスペクトルを算出する。流速変動成分抽出手段６３は、ドップラ変動解析手段６２が算出した海流の流速の周期毎のスペクトルから、その周期により潮流の成分を抽出する。抽出する基準としては、想定する津波の周期より長い周期の成分を抽出する、潮汐の最短の周期である半日や、その倍数となる周期の成分を抽出するなどがある。抽出成分のドップラ算出手段６４は、流速変動成分抽出手段６３が抽出した潮流のスペクトル成分から逆ＦＦＴにより潮流によるドップラ成分ｆｄをセル毎に算出する。
【００７７】
抽出成分のドップラ算出手段６４が算出したドップラｆｄすなわちドップラ変動周期分離手段５３が算出したドップラｆｄにより、実施の形態３と同様に、ドップラ分布シフト手段５１は、ドップラ分布を周波数方向にシフトさせることにより、各セルのドップラ分布から潮流などの影響を除去する。空間積分手段１４では、ドップラ分布シフト手段５１で周波数シフトを行なった結果を積分することにより、ドップラ分布のピークが積み上がり強調され、検出手段１５により、津波を高感度で検出することができる。
【００７８】
なお、ドップラ変動周期分離手段５３の構成として、図１６では、潮流などの影響を算出する為に、フーリエ変換による周波数解析を行なう構成としたが、図１７に示すように、潮流の推定の際に一般によく用いられる調和分解を行なう構成としても良い。図１７の構成では、ドップラ変動周期分離手段５３は、海流のドップラ算出手段６１、調和定数算出手段６２ａ、潮流成分抽出手段６３ａ、および潮流算出手段６４ａにより構成される。
【００７９】
調和分解は、１５日から１年間にわたりセル毎の潮流のドップラ周波数ｆｄを蓄積し、蓄積したドップラ周波数ｆｄの変動を（５）式にあてはめて、各成分に分解し、潮流などの周期的な変動を行なう成分を抽出するものである。
【００８０】
【数５】
2014098686.tif 000007
【００８１】
調和定数算出手段６２ａは、ドップラ周波数ｆｄの変動を最小２乗法などにより（５）式にあてはめて、それぞれの調和定数を算出する。潮流成分抽出手段６３ａは、算出された調和定数から、潮流などの津波と異なる周期を持つ成分を抽出する。潮流算出手段６４ａは、潮流成分抽出手段６３ａが抽出した成分に関する調和定数のみから、（５）式により抽出した潮流などの成分に関する流速を算出する。
【００８２】
ドップラ変動周期分離手段５３としては、図１６および図１７に示す他、ドップラ周波数の周期的変動の違いにより分離する方法であれば、デジタルフィルタなどを使用して分離を行なっても良い。
【００８３】
また、図１５の構成では、検出装置４ｅは、ドップラ変動周期分離手段５３により、津波と周期の異なるドップラ成分を抽出するが、図１８のように、検出装置４ｆは、ドップラ変動分離手段５３により津波と周期の異なるドップラ成分を抽出し、周期の違いでは分離が困難な場合のある暖流、寒流、吹走流などについては、実施の形態３と同様の吹走流予測装置４３と暖流・寒流予測処理４４とにより予測を行う構成としても良い。図１８では、検出装置４ｆは、通常流合成装置４５ａによりドップラ変動分離手段５３の抽出結果と吹走流予測装置４３と暖流・寒流予測処理４４の予測結果とを合成し、ドップラ分布シフト手段５１ｂによりドップラ分布から津波以外の通常海流の影響を差し引く。空間積分手段１４は、津波以外の通常海流の流速によるドップラ周波数の影響の予測値が差し引かれたドップラ分布を積分する。このように構成すると、空間積分手段１４の出力では、ドップラ分布のピークが積み上がり強調され、検出手段１５により、津波を高感度で検出することができる。
【００８４】
さらに、図１５の構成では、検出装置４ｅは、ドップラ分布シフト手段５１によりメモリ１２上のドップラ分布のシフトを行い、空間積分手段１４でドップラ分布の積分を行うが、図１９の構成の様に、検出装置４ｇは、ドップラ分布シフト手段５１によりメモリ１２上のドップラ分布のシフトを行なった結果をメモリ５２に入力し、メモリ５２上のドップラ分布を空間積分手段１４により積分する構成としても良い。
【００８５】
さらに、検出手段１６による津波の検出、時間積分手段１７と検出手段１８による津波の検出も、メモリ１２からではなく、メモリ５２上のドップラ分布のシフトを行ったドップラ分布により津波の検出に係る処理を行う構成とする。このような構成とした場合、検出手段１６による津波の検出、時間積分手段１７と検出手段１８による津波の検出、および空間積分手段１４と検出手段１５による津波の検出ではそれぞれメモリ５２上の、ドップラ分布から津波以外の通常海流の影響を差し引いたドップラ分布により津波の検出に係る処理を行う。このため、メモリ５２上の、ドップラ分布から津波以外の通常海流の影響を除いたドップラ分布を使用するため、検出手段１６、検出手段１８における津波の検出においても、津波以外の通常海流の流速のセル毎の違いや時間変動に影響されずに検出を行なうことができる。
【００８６】
加えて、実施の形態４に係る海洋レーダ装置の、ドップラ分布シフト手段５１ｂにより、ドップラ分布から、ドップラ変動周期分離手段５３で求めた潮流の影響を差し引く構成は、図２０のように、受信機２ａおよび受信機２ｂなどの異なる場所に設置された２台以上のレーダで受信した受信データから津波を検出し、検出結果を合成する構成に対しても適用することが可能である。図２０では、受信機２ａで受信された電波は、信号処理装置３ａ中の方位分離手段６ａ、距離分離手段７ａ、速度分離手段８ａで処理され、空間積分手段１４ａで積分される。これらの処理は、図１５において、受信機２で受信された電波が、信号処理装置３中の方位分離手段６、距離分離手段７、速度分離手段８で処理され、空間積分手段１４で積分されるそれぞれの処理と全く同等である。また、図２０における受信機２ｂで受信された電波についても同様である。このため、図２０における受信機２ａ、受信機２ｂで受信された電波のそれぞれについても、図１５における、受信機２で受信された電波と全く同様に、ドップラ分布シフト手段５１ｂにより、ドップラ分布から、ドップラ変動周期分離手段５３で求めた潮流の影響を除くことができる。
【００８７】
以上のように、実施の形態４に係る海洋レーダ装置は、ドップラ周波数の周期の違いにより津波と津波以外の通常海流の影響によるドップラ周波数ｆｄを算出し、ドップラ分布シフト手段５１で周波数方向にシフトを行なった結果を積分することにより、ドップラ分布のピークが積み上がり強調され、検出手段１５により、津波を高感度で検出することができる。
【符号の説明】
【００８８】
１、１ａ、１ｂ  送信機
２、２ａ、２ｂ  受信機
３、３ａ  信号処理装置
４、４ａ、４ｂ、４ｃ、４ｄ、４ｅ、４ｆ  検出装置
５  津波予測装置
６、６ａ、６ｂ  方位分離手段
７、７ａ、７ｂ  距離分離手段
８、８ａ、８ｂ  ドップラ分離手段
９  津波予測処理
１０  津波シミュレータ
１１  津波データベース
１２、１２ａ  メモリ
１３  積分区間設定手段
１４、１４ａ  空間積分手段
１５、１５ａ、１６，１６ａ、１８、１８ａ  検出手段
１７、１７ａ  時間積分手段
１９、１９ａ  検出結果統合手段
２０、２０ａ  メモリ
３１、３２  積分区間
４１  通常流予測装置
４２  潮流予測装置
４３  吹走流予測装置
４４  暖流・寒流予測装置
４５、４５ａ  通常流合成処理
５１、５１ａ  ドップラ分布シフト手段
５２  メモリ
５３、５３ａ  ドップラ変動周期分離手段
１００  海洋レーダ装置
【図１】
2014098686.tif 000008
【図２】
2014098686.tif 000009
【図３】
2014098686.tif 000010
【図４】
2014098686.tif 000011
【図５】
2014098686.tif 000012
【図６】
2014098686.tif 000013
【図７】
2014098686.tif 000014
【図８】
2014098686.tif 000015
【図９】
2014098686.tif 000016
【図１０】
2014098686.tif 000017
【図１１】
2014098686.tif 000018
【図１２】
2014098686.tif 000019
【図１３】
2014098686.tif 000020
【図１４】
2014098686.tif 000021
【図１５】
2014098686.tif 000022
【図１６】
2014098686.tif 000023
【図１７】
2014098686.tif 000024
【図１８】
2014098686.tif 000025
【図１９】
2014098686.tif 000026
【図２０】
2014098686.tif 000027