JP6300344B2 - 海流推定方法及び海流推定システム - Google Patents

Description

本発明は、海上における船舶等を用いて海流の方向、流速を推定する海流推定方法、海流推定システムに関する。

海水の流れには、潮汐作用によって生ずる潮流と、潮流とは異なる原因で地球規模において発生する海流がある。潮流は比較的短い時間内において周期的に変動し、天文学的に予測が可能であるのに対し、海流の時間的変化は一般的には潮流よりも緩やかであるがその予測が困難である。一方、船舶の運航に対しては、この海流は大きな影響を及ぼす。このため、海流ベクトル（海流の向き、流速（強度））、あるいはその分布を認識することは、船舶の運航上極めて重要である。例えば、海流を把握した運航をすれば、運航時間の短縮や燃料消費を低減することもできる。

海流は、例えば観測船や係留ブイ等による計測、人工衛星を用いたリモートセンシング、沿岸に設置された海洋レーダー等によって計測される。こうした離散的な実測データを用いたシミュレーション（計算）によって、海流の全体像が得られ、公開されている。こうした場合においては、海流を実測する測定点を増やし、かつ測定の精度を高めることが正確に海流分布を把握するためには必要である。このように潮流、海流を検出するために、各種の技術が提案されている。特許文献１には、船舶の船位情報と、船舶の形状等を元にしてモデル化した運動モデルを用いた計算によって、船位情報としての船舶の位置の変化と、船首方向から潮流を算出する技術が記載されている。船舶の位置は、例えばＧＰＳ（Ｇｌｏｂａｌ Ｐｏｓｉｔｉｏｎｉｎｇ Ｓｙｓｔｅｍ）等を用いて、また船首方向は、例えばジャイロコンパスを用いて海上においても認識することができる。

特許文献２には、船舶の位置の軌跡と船首方向を同時にディスプレイに表示し、潮流による影響を視覚的に把握させることのできる表示装置が記載されている。

特許文献３には、複数の船舶の各々が潮流を測定し、これらの情報を基地局に集約することによって海流の分布を推定し、潮流が実測された箇所以外の箇所における潮流の情報を推定するシステムが記載されている。潮流を実測した船舶以外の船舶は、この情報によって、自己のいる箇所における潮流を認識した運行をすることができる。

特開２００５−１７２６１８号公報 特開平８−１１４４６２号公報 特開２００１−２６４４３７号公報

特許文献１に記載の技術では、船舶の形状等を元にしたモデル化を行い、流体力学的な計算を行って海流の向き、流速を算出している。このため、計算に用いる前提条件が非常に多く、この前提条件が実際の船舶に全て合致することは稀であり、実際の状況と乖離している場合が多くなる。このため、潮流（海流）ベクトルを高精度で求めることは困難であった。

特許文献２に記載の技術では、潮流を視覚的に認識することは容易であるものの、正確にその向き、流速を算出することが困難であることは明らかである。特許文献３に記載の技術は、個々の船舶において潮流が正確に測定できることを前提としているために、その精度は、結局は個々の船舶における潮流ベクトルの検出の精度に依存する。このため、特許文献１に記載の技術と同様に、潮流（海流）に関する情報を高精度で求めることは困難であった。また、個々の船舶が潮流分布測定だけの目的のために、専用の設備を搭載し協力して行動をする必要があった。

従って、海流に関する情報（向き、流速）を正確に測定することは困難であった。

本発明は、かかる問題点に鑑みてなされたものであり、上記問題点を解決するもので、海流に関する情報（向き、流速）を正確に、容易に推定可能な海流推定方法及び海流推定システムを提供することを目的とする。

本発明は、上記課題を解決すべく、以下に掲げる構成とした。
本発明の請求項１に係る海流推定方法は、海上の船舶の運航に関する状況が特定された運航情報を複数用いることによって海流の向き及び強さを示す海流ベクトルを推定する海流推定方法であって、前記運航情報として海流を推定する海域を通る複数の前記船舶の個船毎の対地速度ベクトル及び前記船舶の船首方位を取得し、個船毎の前記対地速度ベクトル及び前記船首方位に基づいて、複数の前記船舶の航行する前記海域における前記海流ベクトルを算出するとともに、前記運航情報を取得する複数の前記船舶の隻数を増すことにより前記海流ベクトルの推定確度を上げることを特徴とする。
この発明においては、少なくとも複数の船舶の対地速度ベクトル、船首方位に関する情報が含まれる運航情報が複数取得され、これらの情報に基づいて海流ベクトル（海流の向き、強さ）が算出される。ここで、特に複数の船舶における運行情報が用いられる。
本発明の請求項２に係る海流推定方法において、前記運航情報には、前記船舶の対地速度が特定された対地速度情報、及び前記船舶の対地方位が特定された対地方位情報が含まれ、前記対地速度及び前記対地方位から求めた前記対地速度ベクトルを用いて前記海流ベクトルを算出することを特徴とする。
この発明においては、各運航情報には、船舶の対地速度（対地速度ベクトルの絶対値）が特定された対地速度情報と、対地方位（対地速度ベクトルの向き）が特定された対地方位情報が含まれる。これによって各運航情報毎に対地速度ベクトルが特定される。
本発明の請求項３に係る海流推定方法は、前記船首方位が特定された船首方位情報、前記対地速度情報、及び前記対地方位情報を船舶自動識別システム（Ａｕｔｏｍａｔｉｃ Ｉｄｅｎｔｉｆｉｃａｔｉｏｎ Ｓｙｓｔｅｍ：ＡＩＳ）によって取得することを特徴とする。
この発明においては、船首方位情報、対地速度情報、対地方位情報は、特に大型の船舶に設置が義務づけられているＡＩＳによって取得される。

本発明の請求項４に係る海流推定方法において、前記運航情報には、前記運航情報が取得された際の前記船舶の位置が特定された位置情報が含まれ、当該位置情報に基づいて前記海域における前記海流ベクトルを算出することを特徴とする。
この発明においては、船舶の位置に基づいて複数の運航情報が選択されて使用される。例えば、海流を推定すべき海域を航行する２隻以上の船舶、近接して航行する２隻以上の船舶等、が選択され、その海域の海流ベクトルの推定（算出）に使用される。ここで、「複数の運航情報」は、必ずしも複数の船舶におけるものでないことは前記と同様である。
本発明の請求項５に係る海流推定方法において、前記運航情報には、前記運航情報が取得された際の前記船舶の位置が特定された位置情報、及び前記位置情報が取得された時刻が特定された時刻情報が含まれ、前記位置情報及び前記時刻情報に基づいて前記対地速度ベクトルを算出することを特徴とする。
この発明においては、対地速度情報と対地方位情報を用いる代わりに、運航情報における位置情報と時刻情報を用い、船舶の位置の時間変化から、対地速度ベクトルが算出される。
本発明の請求項６に係る海流推定方法は、前記船首方位が特定された船首方位情報、及び前記位置情報を船舶自動識別システム（Ａｕｔｏｍａｔｉｃ Ｉｄｅｎｔｉｆｉｃａｔｉｏｎ Ｓｙｓｔｅｍ：ＡＩＳ）によって取得することを特徴とする。
この発明においては、位置情報も船首方位情報と同様に、ＡＩＳによって取得される。

本発明の請求項７に係る海流推定方法は、前記海流ベクトルの算出に用いる情報を、前記運航情報及び/又は前記運航情報の取得状況に基づいて補正することを特徴とする。
この発明においては、海流ベクトルの算出に用いる情報（運航情報によって得られた対地速度ベクトル、船首方位等や、これら以外に海流ベクトルの算出に用いるパラメータ）が、より高精度に海流ベクトルが算出できるように補正されてから使用される。
本発明の請求項８に係る海流推定方法は、前記海流ベクトルに関する情報をパラメータとし、各前記運航情報毎に運航情報別目的関数を作成し、複数の前記運航情報毎の前記運航情報別目的関数の総和である目的関数を最適化する前記海流ベクトルを算出することを特徴とする。
この発明においては、最適化法によって海流ベクトルが推定される。ここでは、複数の運航情報における対地速度ベクトル、船首方位と海流ベクトルを用いて定義される目的関数の最適値（例えば最小値）を与える海流ベクトルが、推定される海流ベクトルであるとされる。目的関数は、運航情報毎に定義された運航情報別目的関数の総和として表される。
本発明の請求項９に係る海流推定方法は、前記海流ベクトルに関する情報をパラメータとして用い各前記運航情報を用いて運航情報別目的関数を作成し、複数の前記運航情報毎の前記運航情報別目的関数に重み係数を乗じた値の総和である目的関数を最適化する前記海流ベクトルを算出し、前記運航情報毎の前記重み係数の設定によって前記補正を行うことを特徴とする。
この発明においては、前記の補正は、運航情報毎の重み係数の設定によって行われる。
本発明の請求項１０に係る海流推定方法は、前記重み係数を、前記運航情報が取得された位置と前記海流を推定する位置との差、前記運航情報が取得された時刻と前記海流を推定する時刻との差、前記運航情報を取得した船舶の種類、あるいは前記運航情報を用いて過去に行われた海流ベクトルの算出結果、のうちの少なくとも１つに基づき設定することを特徴とする。
この発明においては、重み係数の設定は、位置、時刻、船舶の種類、過去の算出履歴に基づいて設定される。これによって、使用される運航情報の優先度が設定される。

本発明の請求項１１に係る海流推定システムは、
請求項１から請求項１０までのいずれか１項に記載の海流推定方法に基づいて海流ベクトルを推定する海流推定システムであって、複数の前記船舶の前記運航情報を取得する情報取得部と、前記取得した前記運航情報を記憶する記憶手段と、前記海域を特定の時刻に航行する複数の前記船舶をリストアップするリストアップ手段と、リストアップされた前記船舶の前記運航情報を前記記憶手段から読み取る読取手段と、読み取られた前記運航情報に基づいて、前記海流ベクトルを算出する海流ベクトル算出部とを具備することを特徴とする。
この発明においては、情報取得部が、複数の船舶の運航情報を入手し、この運行情報は記憶手段に記憶される。これらの船舶のうち、リストアップ手段によって対象となる海域を特定の時刻に航行する船舶がリストアップされ、これに対応する運行情報が記憶手段から読み取られる。海流ベクトル算出部は、このように読み取られた運行情報に基づいて海流ベクトル（海流の向き、強さ）を算出する。
本発明の請求項１２に係る海流推定システムにおいて、前記情報取得部は、予め予測された海象情報を取得し、前記海流ベクトル算出部は、推定された前記海流ベクトルを前記海象情報に基づいて補正することを特徴とする。
この発明においては、例えば外部から入手した海象情報によって潮流を認識し、推定された海流ベクトルから、この潮流成分を除いた本来の海流が推定される。
本発明の請求項１３に係る海流推定システムにおいて、前記海流ベクトル算出部は、複数の前記海域における前記海流ベクトルを算出することにより、海流分布を推定することを特徴とする。
この発明においては、複数の海域における海流ベクトルを推定することにより、海流分布が推定される。
本発明の請求項１４に係る海流推定システムにおいて、前記海流ベクトル算出部は、予め予測された複数の前記海域における海流分布に対して、算出された前記海流ベクトルを同化させ、新たな海流分布を推定することを特徴とする。
この発明においては、計算により予め予測された海流分布が存在する場合に、この海流分布に対して、推定された海流ベクトルのデータを同化させることにより、より高精度の海流分布が算出される。
本発明の請求項１５に係る海流推定システムは、推定された前記海流分布を元にして、航路の指示を海上の船舶に対して行う航路指示部を具備することを特徴とする。
この発明においては、航路指示部が、推定された海流分布を元にして算出された最適な航路の指示を行う。

本発明の海流推定方法、海流推定システムは以上のように構成されているので、海流ベクトルを高精度で推定することができる。
ここで、対地速度ベクトルと船首方位は、船舶の運航に関する状況を特定する運航情報によって求めることができるため、船舶の通常の設備と運航情報を利用して海流ベクトルの推定を特に容易に行うことができる。特に、正確な対地速度情報、対地方位情報、船首方位情報は、大型の船舶に搭載が義務付けられているＡＩＳによって容易に入手することができる。これにより、特別の運航情報取得手段を用いることなく、上記の算出を高精度かつ容易に行うことができる。
また、位置情報、時刻情報も、対地速度情報及び船首方位情報と共に使用することによって、より高精度で海流ベクトルの算出を行うことができる。この際、位置情報も、ＡＩＳによって正確な情報を容易に入手することができる。位置情報と時刻情報に基づいて対地速度ベクトルを正確に算出し、これによってより高精度に海流ベクトルを算出することもできる。また、位置情報によって海流ベクトルを得た船舶の航行する海域を特定することができる。
また、運航情報を用いて目的関数を作成し、最適化法によって海流ベクトルを推定することにより、より広範囲のデータを用いることによって、より信頼性の高い海流ベクトルを推定することができる。この際、海流ベクトルを算出すべき位置、時刻等を考慮して運航情報別目的関数の重み付けを設定すれば、位置、時刻等に応じて最も確度の高い海流ベクトルを推定することができる。
また、複数の海域における海流ベクトルを推定し、広域の海流分布も推定することができる。この際、海象情報を用いて、得られた海流分布から潮流成分を除いた本来の海流分布を得ることもできる。
また、予め推定される低精度の海流分布がある場合には、推定された海流ベクトルのデータをこれに同化させ、より高精度の海流分布を得ることができる。こうして推定された海流分布を用いて、最適化された航路の指示を行うこともできる。

本発明の実施の形態に係る海流推定方法（船舶が２隻以下の場合）の原理を説明する図である。 本発明の実施の形態に係る海流推定方法（船舶が３隻の場合）の原理を説明する図である。 本発明の実施の形態に係る海流推定システムの構成の一例を模式的に示す図である。

以下、本発明を実施するための形態となる海流推定システム、海流推定方法について説明する。図１は、この海流推定システムにおいて海流を推定（算出）する方法の原理を説明する図である。まず、ここでは、２隻以下の船舶を用いる場合について説明する。図１（ａ）は、１隻の船舶１の海上における運航状況を上側から見た図であり、単位時間あたりにおける船舶１の動きを模式的に示している。ただし、図１においては、潮流も海流に含め、全体的な海水の流れを海流と呼称することにする。

ここで、海流がない場合には、船舶１は常にその船首方位（船首の向く方向）に沿って対水速度（水面に対する速度）ベクトルＶ_Ｗをもって進行し、実際の船舶１の速度（対地速度：絶対座標（地上）に対する速度）ベクトルＶ_Ｓは、対水速度ベクトルＶ_Ｗと等しくなることは明らかである。しかしながら、海流ベクトルＶ_Ｃが零でない場合には、対地速度ベクトルＶ_Ｓは、対水速度ベクトルＶ_Ｗと海流ベクトルＶ_Ｃの和となる。この際、特に船首方位と垂直な成分をもつ海流ベクトルＶ_Ｃが存在した場合には、図１（ａ）に示されるように、船舶１の船首方位は実際の進行方向とずれる。ここで、船舶１自身に搭載されたＧＰＳ信号を用いた機器により、船舶１自身の位置情報は高精度で認識することができ、この位置情報の時間変化によって、対地速度ベクトルＶ_Ｓ（実際の速度の向き、絶対値）も高精度で検出することが可能である。ここで、対地速度ベクトルＶ_Ｓを検出することは、対地速度ベクトルＶ_Ｓの絶対値と、その向きである対地方位が検出されることを意味する。対地方位は、図１（ａ）に示されるように、例えば北方向を基準とした時計回りの角度θ_Ｓとして定義される。また、船舶１の船首方位も、同様に角度θ_Ｗとして定義すれば、同様に高精度で認識することができる。船首方位θ_Ｗは、対水速度ベクトルＶ_Ｗの向となる。なお、以下の説明においては、対地速度ベクトルの絶対値を対地速度、対地速度ベクトルの向き（θ_Ｓ）を対地方位、対水速度ベクトルの絶対値を対水速度と呼称する。海流ベクトルＶ_Ｃが零の場合（海流が存在しない場合）には、対地速度と対水速度は等しく、対地方位と船首方位は等しくなる。

ここで、図１（ａ）より、対水速度ベクトルＶ_Ｗも既知となれば、対地速度ベクトルＶ_Ｓとの差分によって海流ベクトルＶ_Ｃを算出することができる。これに対して、船舶に装着した流速計等を用いて対水速度ベクトルＶ_Ｗを検出することは不可能ではない。しかしながら、ＧＰＳ信号等を用いて自己の位置を認識し、その時間変化より正確に認識される対地速度ベクトルＶ_Ｓと異なり、対水速度ベクトルＶ_Ｗを測定するための基準となるのは海水であるため、これを対地速度ベクトルＶ_Ｓと同程度に高精度で検出することは困難である。このため、前記の通り、対水速度ベクトルＶ_Ｗについては、その方向のみがＶ_Ｓと同様に船首方位θ_Ｗとして高精度で認識され、その絶対値Ｖ_Ｗ０を高精度で検出することは困難である。

この場合、図１（ｂ）において、船首方位θ_Ｗと対地速度ベクトルＶ_Ｓが正確に定まっても、海流ベクトルＶ_Ｃの始点（対水速度ベクトルＶ_Ｗの終点）が定まらないため、海流ベクトルＶ_Ｃを一義的に求めることができない。例えば、図１（ｂ）においては、ベクトルＶ_Ｃ１、Ｖ_ＣＰ、Ｖ_Ｃ２のいずれもが海流ベクトルＶ_Ｃの候補となる。ここで、ベクトルＶ_ＣＰは船首方向（対水速度ベクトルＶ_Ｗ）と垂直である。この結果より、船首方向と垂直な方向の海流ベクトルＶ_Ｃの成分は、ベクトルＶ_Ｃ１、Ｖ_Ｃ２においてもベクトルＶ_ＣＰとなることは明らかである。すなわち、船舶１の対地速度ベクトルＶ_Ｓと船首方位θ_Ｗが認識できた場合、海流ベクトルＶ_Ｃの向き、絶対値を正確に算出することは困難であるが、海流ベクトルＶ_Ｃの船首方位と垂直な方向成分は一義的に算出することができる。船舶の操舵においては、海流ベクトルＶ_Ｃの船首方位と垂直な成分が特に大きな影響を与えることは明らかである。このため、対水速度ベクトルＶ_Ｗを用いずに、船舶１の対地速度ベクトルＶ_Ｓと船首方位θ_Ｗのみを用い、図１（ｂ）の作図を行うことによってＶ_ＣＰを算出する上記の方法は有効である。

対水速度ベクトルＶ_Ｗの方向（船首方位）のみが既知でありその絶対値（対水速度）が不明である場合、海流ベクトルＶ_Ｃの他成分、例えば船首方位における成分を算出することは不可能であり、海流ベクトルＶ_Ｃの向き、絶対値を一義的に求めることは不可能である。しかしながら、対地速度ベクトルＶ_Ｓ１、船首方位θ_Ｗ１である船舶１の他に、船舶２（対地速度ベクトルＶ_Ｓ２、船首方位θ_Ｗ２）が航行しており、その対地速度ベクトルＶ_Ｓ２、船首方位θ_Ｗ２が正確に検出され、かつ海流ベクトルＶ_Ｃが船舶１、２に対して共通である場合には、船舶１、２の対水速度ベクトルＶ_Ｗ１、Ｖ_Ｗ２の絶対値が不明であっても、図１（ｃ）に示されるように、海流ベクトルＶ_Ｃは一義的に定まる。例えば、船舶１と船舶２の距離が近接していれば、これらの間で海流ベクトルＶ_Ｃは共通と考えることができる。この場合には海流ベクトルＶ_Ｃ（海流の向き、強度）を正確に算出することができる。あるいは、船舶１、２が離間していても、これらが航行する海域において海流ベクトルＶ_Ｃが等しいと考えられる場合にも、同様である。

具体的には、この算出方法は、以下の通りとなる。ここで、図１（ａ）に示されるようにｘ、ｙ軸をとり、対水速度ベクトルＶ_Ｗを（Ｖ_Ｗｘ、Ｖ_Ｗｙ）とすると、これは、対水速度をＶ_Ｗ０として、Ｖ_Ｗ０×（ｃｏｓθ_Ｗ、ｓｉｎθ_Ｗ）と書ける。ここで、Ｖ_Ｗ０が不定であり、θ_Ｗは既知となる。対地速度ベクトルＶ_Ｓを（Ｖ_Ｓｘ、Ｖ_Ｓｙ）とすると、これも既知となる。算出すべき海流ベクトルＶ_Ｃを（Ｖ_Ｃｘ、Ｖ_Ｃｙ）とすると、これらの間には、（１）（２）式が成立する。

ここで、（１）（２）式より、不定であるＶ_Ｗ０を消去すると、（３）式が得られる。

船舶１、２の各々について（３）式が成立し、Ｖ_Ｃｘ、Ｖ_Ｃｙは両者に共通であると仮定すると、船舶１、２についてのＶ_Ｓｘ、Ｖ_Ｓｙ、θ_Ｗについては添字１、２を付して、（４）式が成立する。

これにより、海流ベクトルＶ_Ｃ＝（Ｖ_Ｃｘ、Ｖ_Ｃｙ）は（５）式によって算出される。

このため、θ_Ｗ１≠θ_Ｗ２の場合、すなわち、船舶１と船舶２の船首方位が異なる場合には、（Ｖ_Ｃｘ、Ｖ_Ｃｙ）を算出することができる。北方向を基準とした海流の方位（θ_Ｗ、θ_Ｓと同様に定義）はｔａｎ^−１（Ｖ_Ｃｘ／Ｖ_Ｃｙ）であり、海流の絶対値（強さ）は（Ｖ_Ｃｘ ^２＋Ｖ_Ｃｙ ^２）^１／２となる。

上記の通り、海流ベクトルＶ_Ｃの算出に必要なのは、各船舶の対地速度ベクトルＶ_Ｓ（Ｖ_Ｓｘ、Ｖ_Ｓｙ）と、船首方位θ_Ｗである。対地速度ベクトルＶ_Ｓの算出に必要なのは、対地速度が特定された対地速度情報と、対地方位θ_Ｓが特定された対地方位情報である。更に、船首方位θ_Ｗが特定された船首方位情報が必要となる。前記の通り、自船におけるこれらの情報は、ＧＰＳ等を用いて得ることができる。

一方、「１９７４年の海上における人命の安全に関する条約（ＳＯＬＡＳ７４）」によって、大型の船舶、例えば国際航海に従事する３００総トン以上の船舶には、船舶自動識別システム（Ａｕｔｏ Ｉｄｅｎｔｉｆｉｃａｔｉｏｎ Ｓｙｓｔｅｍ：ＡＩＳ）搭載が義務付けられている。ＡＩＳによれば、船舶は、自船に関するＡＩＳ情報をＶＨＦ帯の電波で発信する義務が発生し、他船のＡＩＳ情報も受信して入手することができる。これにより、船舶の航行の安全を図ることができる。ＡＩＳの詳細についての詳細は、例えばｈｔｔｐ：／／ｗｗｗ６．ｋａｉｈｏ．ｍｌｉｔ．ｇｏ．ｊｐ／ｏｓａｋａｗａｎ／ｏｔｈｅｒｓ／ａｉｓ／ａｉｓ．ｈｔｍ等に記載されている。ＡＩＳ情報のうち、船舶の状態を示す動的情報としては、船舶の位置情報、時刻（世界標準時）、対地方位（対地針路）、対地速度、船首方位、回頭率（船首方位の時間変化率）等がある。

このため、前記の対地速度情報、対地方位情報、船首方位情報は、ＡＩＳ情報に含まれる。対地速度ベクトルＶ_Ｓは、ＡＩＳ情報における対地速度、対地方位から求まる。あるいは、ＡＩＳ情報における位置情報の時間変化から、対地速度ベクトルを算出することもできる。すなわち、対地速度ベクトルＶ_ＳをＡＩＳ情報から一義的に算出することができる。このため、ＡＩＳを搭載した船舶であれば、自船の対地速度ベクトルＶ_Ｓ、船首方位θ_Ｗを認識することができる。これにより、海流ベクトルＶ_Ｃの船首方位と垂直な成分（Ｖ_ＣＰ）を算出することが可能となる。

また、ＡＩＳを搭載した船舶は、他船のＡＩＳ情報（対地速度情報、対地方位情報、船首方位情報）を受信することもできる。このため、この船舶は、図１（ｃ）における対地速度ベクトルＶ_{Ｓ１（２）}、船首方向θ_{Ｗ１（２）}を同様に認識することができる。更に、ＡＩＳの受信機能をもつ地上局においても、ＡＩＳを搭載する複数の船舶におけるこれらの情報を認識することができる。このため、海流ベクトルＶ_Ｃの算出を、船舶上、地上局のどちらにおいても行うことができる。

なお、上記の計算では、図１（ｃ）に示されるように、２隻の船舶１、２における対地速度ベクトル、船首方位情報が使用された。ここで、船舶１、２は異なる船舶である必要はなく、例えばある時間内において１隻の船舶が船首方位を変えて航行し、この時間内にこの船舶が航行した海域における海流ベクトルＶ_Ｃが変わらない場合には、この船舶における異なる時刻あるいは位置における対地速度ベクトル、船首方位情報を上記の２隻の対地速度ベクトル、船首方位情報の代わりに用いることもできる。

上記の算出において使用する情報は、ＡＩＳ情報に限定されない。例えば、海運企業毎の船舶の運行状況をモニターした情報であり、その中に対地速度情報、対地方位情報、船首方位情報が含まれるものであれば、同様に用いることができることは明らかである。また、例えば、１隻の船舶における異なる複数の時刻の情報を用いることもできる。すなわち、上記の算出を行うためには、船舶の対地速度情報、対地方位情報、船首方位情報が含まれる任意の運航情報（海上の船舶の運航に関する状況が特定された情報）を複数用いることができる。１隻の船舶における異なる複数の時刻の情報の場合には、この運航情報は、ＡＩＳ情報のように、船舶から外部に発信された情報である必要もない。

上記においては、船舶が２隻（あるいは運航情報が２つ）あり、それぞれの海域における海流ベクトルＶ_Ｃが厳密に等しいものとしたために、海流ベクトルＶ_Ｃを一義的に算出することができた。しかしながら、実際には、両者の海域間で海流ベクトルＶ_Ｃは厳密には等しくない場合も多い。こうした場合には、３隻以上の船舶のデータ（３つ以上の運航情報）を用い、最も確度の高い海流ベクトルを推定することができる。また、２隻の船舶（２つの運航情報）を用いた場合には上記の通り海流ベクトルＶ_Ｃは一義的に定まるものの、算出に利用されるＶ_{Ｓ１（２）}、θ_{Ｗ１（２）}の測定誤差に起因した誤差が発生する。このような誤差を低減し、推定される海流ベクトルＶ_Ｃの確度を高くするためには、これより多くの船舶（運航情報）のデータを用いる方が好ましいことは明らかである。この場合における海流ベクトルＶ_Ｃの算出方法（推定方法）について以下に説明する。

まず、海流ベクトルＶ_Ｃが船舶１、２、３の３隻について厳密に同一となっている場合には、図１（ｃ）と同様の関係が船舶１、２、３のうちの任意の２隻間で成立するため、上記の各ベクトルの関係は図２（ａ）に示される通りとなる。この図においては、各々の対水速度ベクトルＶ_Ｗの終点、あるいは各船舶の船首方位を向く直線（対水速度ベクトルＶ_Ｗの方向）の交点が点Ｃ_１で一致すると共に、各々の対地速度ベクトルＶ_Ｓの終点は点Ｃ_２で一致する。図２（ａ）においては、海流ベクトルＶ_Ｃは、点Ｃ_１を始点、点Ｃ_２を終点とするベクトルとなり一義的に定まる。

しかしながら、厳密には図２（ａ）の関係が成立していない場合には、図２（ｂ）に示されるように、海流ベクトルＶ_Ｃの終点となる点Ｃ_２（ｘ_１、ｘ_２）を定めた場合に、各船舶の船首方位を向く直線は共通の交点をもたない。このため、海流ベクトルＶ_Ｃの始点となるべきＣ_１の位置を定めることができず、海流ベクトルＶ_Ｃを厳密に算出することができない。この場合には、最適化法として知られるように、海流ベクトルＶ_Ｃの始点とすべき点Ｐ（ｘ_０、ｙ_０）におけるｘ_０、ｙ_０を、ある目的関数が最小値をとるような値として決定することができる。

この最適化において、図２（ｃ）に示されるように、点Ｐは、船舶１の対水速度ベクトルＶ_Ｗ１に対応する直線Ｌ_１、船舶２の対水速度ベクトルＶ_Ｗ２に対応する直線Ｌ_２、船舶３の対水速度ベクトルＶ_Ｗ３に対応する直線Ｌ_３との間の距離（それぞれｈ_１、ｈ_２、ｈ_３とする）がいずれも短くなる点として定義することができる。図２（ｃ）に示されるように、この点Ｐは、直線Ｌ_１、Ｌ_２、Ｌ_３で囲まれる３角形の中に存在する。

推定される海流ベクトルＶ_Ｃの始点となる点Ｐの座標（ｘ_０、ｙ_０）としては、例えば、点Ｐと各直線Ｌ_１、Ｌ_２、Ｌ_３との間の距離の２乗和を目的関数とし、この値が最小となるように設定することができる。直線Ｌがａｘ＋ｂｙ＝１で示される場合に、点Ｐ（ｘ_０、ｙ_０）と直線Ｌとの間の距離ｈは、（６）式で与えられる。ここで、点Ｐと直線Ｌとの間の距離とは、直線Ｌと直交し点Ｐを通る直線と直線Ｌとの交点と点Ｐとの間の距離である。

このため、点Ｐから各直線Ｌ_ｉまでの距離ｈをｈ_ｉとし、ｈ_ｉの２乗を運航情報別目的関数とし、最小とすべき目的関数ｆを、運航情報別目的関数ｈ_ｉ ^２の総和（ｈ_１ ^２＋ｈ_２ ^２＋ｈ_３ ^２）とすることができる。この場合、直線Ｌ_ｉ（ｉ＝１、２、３）がａ_ｉｘ＋ｂ_ｉｙ＝１で表されるものとして、目的関数ｆは（７）式で与えられる。ここで、各船舶の対地速度ベクトルＶ_Ｓ１、Ｖ_Ｓ２、Ｖ_Ｓ３は既知であるために、図２（ｂ）における各船舶の位置は、対地速度ベクトルＶ_Ｓ１、Ｖ_Ｓ２、Ｖ_Ｓ３の終点が１点Ｃ_２となるようにして定まる。また、各船舶の船首方位θ_Ｗ１、θ_Ｗ２、θ_Ｗ３も既知であるため、直線Ｌ_ｉが船舶ｉの位置を通りこの船首方位θ_ｉを向く直線となるように、ａ_ｉ、ｂ_ｉ（ｉ＝１、２、３）は一義的に定まる。

目的関数ｆを最小化するｘ_０、ｙ_０は、∂ｆ／∂ｘ_０＝０、∂ｆ／∂ｙ_０＝０として定まる。その結果、Ｐの座標ｘ_０、ｙ_０は、（８）、（９）式で与えられる。ここで、Ａ、Ｂ、Ｃ、Ｄ、Ｅは（１０）〜（１４）式で与えられる。すなわち、これによってｘ_０、ｙ_０を算出することができ、推定される海流ベクトルＶ_Ｃは、点Ｐから点Ｃ_２に向かうベクトルと推定されるため、（Ｖ_Ｃｘ、Ｖ_Ｃｙ）＝（ｘ_１−ｘ_０、ｙ_１−ｙ_０）となる。

すなわち、この計算においては、海流ベクトルＶ_Ｃに関する情報であるｘ_０、ｙ_０をパラメータとし、使用した運航情報（対地速度情報、対地方位情報、船首方位情報）毎に運航情報別目的関数を算出し、その総和として目的関数を定義し、海流ベクトルＶ_Ｃは、この目的関数を最適化（上記の場合には最小化）するものとして算出（推定）される。

図２の例では、船舶が３隻の場合について記載したが、船舶が４隻以上、あるいは３つ以上の運航情報を用いる場合でも、同様に海流ベクトルＶ_Ｃを算出することができることは明らかである。

上記の式（１）〜（１４）を用いて海流ベクトルＶ_Ｃを算出する場合においては、運航情報（対地速度情報、対地方位情報、船首方位情報）で特定された対地速度ベクトル、船首方位を用いて海流ベクトルＶ_Ｃを算出した。しかしながら、海流ベクトルＶ_Ｃをより高精度に算出するために、運航情報で特定された対地速度ベクトル、船首方位を補正し、これらを用いて海流ベクトルＶ_Ｃを算出することもできる。また、この補正は、対地速度ベクトル、船首方位に対してだけでなく、海流ベクトルＶ_Ｃの算出に用いる他の情報（例えば運航情報別目的関数に関する情報）についても行うことができる。

例えば、対地速度情報、船首方位情報は、一定の時間間隔で取得される場合が多く、この時間間隔は、船舶毎に異なる場合もある。このため、複数の船舶のデータを用いる場合には、全ての船舶において、海流ベクトルを推定すべきある定まった時刻の対地速度情報、船首方位情報が存在するとは限らない。こうした場合には、各船舶のＡＩＳ情報における時刻情報を元にして、同一時刻のデータとなるように補正された対地速度ベクトル、船首方位等を上記の計算に用いることができる。この補正においては、例えば、海流ベクトルを推定すべき時刻と異なる時刻において取得された対地速度ベクトル、船首方位等のデータを内挿、又は外挿すること等によって、新たに対地速度ベクトル、船首方位等のデータを作成する。この新たなデータを、海流ベクトルを推定すべき時刻におけるデータとして扱うことができる。船首方位については、ＡＩＳ情報における回頭率のデータを用いることもできる。

また、例えば、船舶１が船舶２と船舶３の中間付近の海域を航行しており、船舶２と船舶３の船首方位θ_Ｗ２、θ_Ｗ３がほぼ同一であり、船舶１の船首方位θ_Ｗ１はこれとは異なる場合について考える。この場合には、船舶２の対地速度ベクトルＶ_Ｓ２と船舶３の対地速度ベクトルＶ_Ｓ３を内挿することにより、仮に船舶１と同一の海域を船首方位θ_Ｗ３（θ_Ｗ２）で航行した場合における対地速度ベクトルＶ_Ｓ２３を算出することができる。この場合には、船舶１が航行する海域における海流ベクトルＶ_Ｃを算出するためには、この内挿された対地速度ベクトルＶ_Ｓ２３、船首方位θ_Ｗ３（θ_Ｗ２）と、船舶１において実測された対地速度ベクトルＶ_Ｓ１、船首方位θ_Ｗ１を用いて、海流ベクトルＶ_Ｃの算出を行うことができる。すなわち、時刻情報に基づいて対地速度ベクトルＶ_Ｓを補正するだけでなく、位置情報に基づいて補正した対地速度ベクトルＶ_Ｓを使用することにより、より正確に海流ベクトルＶ_Ｃを推定することも可能となる。あるいは、時刻情報、位置情報の両方に基づいて補正した対地速度ベクトルＶ_Ｓ、船首方位θ_Ｗを用いてもよい。

以上の２つの例においては、より高精度に海流ベクトルＶ_Ｃを算出するために対地速度ベクトルＶ_Ｓ、船首方位θ_Ｗが補正された。これに対して、運航情報別目的関数に関する情報を補正する場合について以下に説明する。

目的関数ｆ（（７）式）を用いた上記の例では、各船舶毎に運航情報別目的関数ｈ_ｉ ^２を算出し、その総和を目的関数ｆとした。この際、例えば、船舶１と船舶２は海流ベクトルＶ_Ｃを算出すべき海域の近傍で航行しており、船舶３はこの海域からやや離れた海域にいる場合もある。こうした場合には、各船舶によるデータに重み付けをし、船舶３によるデータ（ｈ_３）の重み付けを船舶１、２によるデータ（ｈ_１、ｈ_２）よりも低くすることが好ましい。具体的には、船舶ｉの重み係数をｗ_ｉとして、（７）式の目的関数ｆの代わりに、（１５）式で定められる目的関数ｇを最小とするようにｘ_０、ｙ_０を定めればよい。この例の場合には、ｗ_１、ｗ_２＞＞ｗ_３とすればよい。前記の目的関数ｆ（（７）式）は、ｗ_ｉ≡１（ｉ＝１、２、３）とし、運航情報別の重み付けを行わなかった場合に相当し、この重み係数の設定が、運航情報別目的関数に関する情報の補正に相当する。

この場合には、上記のＡ〜Ｅは、（１０）〜（１４）式の代わりに下記の（１６）〜（２０）式で与えられるＡ〜Ｅを用い、ｘ_０、ｙ_０は、上記と同様に、（８）（９）式で与えられる。

重み係数ｗ_ｉの設定として、例えば本来の海流ベクトルＶ_Ｃを算出したい海域に近い位置にいる船舶について大きな値、この海域から遠い位置にいる船舶について小さな値となるように、海流ベクトルＶ_Ｃを算出したい海域からの船舶の距離の関数としてｗ_ｉを設定することができる。すなわち、この重み係数ｗ_ｉも、船舶ｉのＡＩＳ情報（位置情報）から定めることができる。

同様に、本来の海流ベクトルＶ_Ｃを算出すべき時刻に近い時刻に取得された対地速度情報、船首方位情報をもつ運航情報について大きな値、海流ベクトルＶ_Ｃを算出すべき時刻から離れた時刻のデータをもつ運航情報について小さな値となるように、ｗ_ｉをこの時間差の関数とすることもできる。前記の通り、各船舶の位置情報、対地速度情報、船首方位情報が取得された時刻に関する情報（時刻情報）は、ＡＩＳ情報から読みとることができるため、この設定も、ＡＩＳ情報を元にして行うことができる。あるいは、重み係数ｗ_ｉを、位置、時刻の両方の関数として設定することもできる。あるいは、運航情報における回頭率によってｗ_ｉを設定する、例えば回頭率が大きな場合には、船舶の運航に海流の影響が出にくいと判断してｗ_ｉを小さくすることができる。運航情報に対応した船舶の種類によってｗ_ｉを設定する、例えば風の影響を受けやすい船舶のデータは信憑性が低いと判断してｗ_ｉを小さくすることができる。また、同じ運航情報を用いた過去の計算結果の履歴を元にｗ_ｉを設定する、例えば過去の計算において図２（ｃ）のｈｉが大きくなった運航情報に対しても、ｗ_ｉを小さくすることができる。このように、海流に関するデータとしての信憑性が低いと判断されるデータに対しては重み係数ｗ_ｉを小さく設定することにより、海流ベクトルＶ_Ｃの算出の精度を高めることができる。

また、例えば２隻の船舶のデータが用いられ、この２隻によるデータの数が大きく異なる場合、データの数が多い方の船舶のデータの重み係数ｗ_ｉを小さくすることにより、２隻の船舶のデータを実質的に対等に扱う、等の処理を行うこともできる。

すなわち、重み係数の設定は、運航情報が取得された位置と海流を推定する海域の位置との差、運航情報が取得された時刻と海流を推定する時刻との差、単一の船舶において取得され用いられた運航情報の数、等に基づいて行うことができる。他の基準によってこの設定をすることもできる。また、前記の例のように補正された対地速度ベクトル等を用い、更にこの重み係数の設定を行うこともできる。

上記の通りに海流ベクトルＶ_Ｃを算出（推定）する海流推定システム全体の構成について図３に示す。図３においては、船舶１〜５が海流１００に乗って航行している。各船舶はＡＩＳを搭載しており、船舶１〜５の各々は、自船のＡＩＳ情報を発信する。このＡＩＳ情報は、他船、あるいは地上にある情報取得部１０のアンテナ１１で受信することができる。情報取得部１０は、受信したＡＩＳ情報を、サーバ（海流ベクトル算出部）２０にネットワーク等を介して送信する。サーバ２０は、受信したＡＩＳ情報を記憶装置２１にデータベースとして記憶し、これを用いて上記の計算を行う。航路指示部２２は、最適な航路の指示を船舶１〜５に対して行う。なお、情報取得部（ＡＩＳ情報取得部）１０に接続されたアンテナ１１と、サーバ（海流ベクトル算出部）２０に接続されたアンテナ１１とでは、通信容量、通信に使用される波長や信号方式等に応じ、それぞれにおいて最適なものを用いることができる。この点については、船舶１〜５において用いられるアンテナについても同様である。

例えば、サーバ２０は、ある特定の海域近辺を特定の時刻近くに航行する船舶をリストアップし、そのＡＩＳ情報を記憶装置２１から読み取り、この海域のこの時刻における海流ベクトルＶ_Ｃを推定することができる。この時刻として最新（現在）の時刻を設定することによって、最新の海流情報を推定することができ、これを各船舶に発信することができる。

また、この操作を複数の海域、複数の時刻において行うことにより、現在までの海流の時間変化も推定することができる。これにより、今後の海流を予測することも可能となる。

また、サーバ２０は、海流ベクトルＶ_Ｃを推定する際に用いる対地速度ベクトル、船首方位に関する前記の位置情報、時刻情報等に基づいた補正も行うことができる。

図３の構成においては、サーバ（海流ベクトル算出部）２０、記憶装置２１を各船舶とは別に地上に設置し、このために情報取得部１０も地上に設置するものとした。しかしながら、これらは通常使用されるコンピュータを用いて構成することができ、単純な構成であるため、これらを各船舶に搭載することも可能である。この場合には、自船のＡＩＳは、自船のＡＩＳ情報を取得するための情報取得部として機能する。他船のＡＩＳ情報は、情報取得部１０が入手し、上記の処理を行うことができる。ここで算出された海流ベクトルに関する情報を、他船、あるいは地上に向かって発信することも可能である。

この際、サーバ２０は、単なる海流ベクトルの推定（算出）だけでなく、この結果を用いた各種の処理も行うことができる。まず、前記の通り、海水の流れには、潮汐作用によって生ずる潮流と、潮流とは異なる原因で地球規模において発生する海流（本来の海流）があり、ここで測定された海流ベクトルは、本来の海流に潮流が加わった流れである。ここで、潮流は周期的であり予測が容易であるために、潮流に関するデータは、海象情報として外部から入手することが可能であり、このデータは予め記憶装置２１に記憶させておくことができる。この場合、この潮流のデータを用いて前記の海流ベクトルを補正し、本来の海流に対応するベクトルを算出することができる。

また、海上の複数の海域における海流を推定することにより、広域における海流の分布を推定することができる。時刻毎に推定された海流の分布をデータベース化して記憶装置２１に記憶させ、かつ各船舶にもこれを発信することができる。各船舶は常時ＡＩＳ情報を発信しているために、更新されたＡＩＳ情報をサーバ２０が受信することにより、こうして推定された海流分布を更新して記憶装置２１に記憶させることができる。この際、海象情報を元にして潮流がない場合に対応して算出された本来の海流分布を記憶装置２１に記憶すれば、本来の海流分布の時間変化も容易に追跡することができる。

この際、予め算出された低精度の海流分布のデータがある場合、これを記憶装置２１が記憶し、上記の海流推定方法によって推定された海流ベクトルをこの海流分布に同化させ、より高精度の海流分布を推定すること（データ同化）もできる。こうした同化方法は、例えば、特開２０１０−２２３６３９号公報に記載されている。こうした同化処理をサーバ２０が記憶装置２１を用いて行うことも可能であり、この高精度とされた海流分布に関する情報をサーバ２０が外部に発信することもできる。

また、サーバ２０が船舶の航路上の海流分布を推定できれば、これに基づいた船舶の最適航路を算出することもできる。この場合、例えば、対水速度の絶対値Ｖ_Ｗ０を定めた場合に、図１（ａ）における対地速度ベクトルＶ_Ｓが目的地まで最短距離を結ぶように船首方位θ_Ｗを設定することができる。これにより、同じＶ_Ｗ０の場合に航行距離を最も短くし、船舶の燃料消費を低減することが可能となる。すなわち、サーバ２０は、船舶の航行に関する指示も行うことができる。図３における航路指示部２２は、このように推定された前記海流分布を元にして、最適な航路の指示を海上の船舶１〜５に対して行う。なお、最適な航路指示は、海流ベクトルＶ_Ｃの算出に当っての運航情報の提供に係わらない任意の船舶を含めて行うことができる。

上記の通り、上記の海流推定システム、海流推定方法は、海上の船舶を用いて海水の流れを推定するために用いられるが、例えば、広い湖沼や河川等においても、同様にその流れを推定することができることは明らかである。また、船舶としても、ＡＩＳを搭載する大型船だけでなく、小型の模型船舶等を用いても同様の測定（推定）ができることは明らかである。

１〜５ 船舶
１０ 情報取得部
１１ アンテナ
２０ サーバ（海流ベクトル算出部）
２１ 記憶装置
２２ 航路指示部
１００ 海流
θ_Ｗ、θ_Ｗ１、θ_Ｗ２ 船首方位
Ｖ_Ｓ、Ｖ_Ｓ１、Ｖ_Ｓ２ 対地速度ベクトル
Ｖ_Ｃ 海流ベクトル

Claims (15)

1. 海上の船舶の運航に関する状況が特定された運航情報を複数用いることによって海流の向き及び強さを示す海流ベクトルを推定する海流推定方法であって、
   前記運航情報として海流を推定する海域を通る複数の前記船舶の個船毎の対地速度ベクトル及び前記船舶の船首方位を取得し、個船毎の前記対地速度ベクトル及び前記船首方位に基づいて、複数の前記船舶の航行する前記海域における前記海流ベクトルを算出するとともに、前記運航情報を取得する複数の前記船舶の隻数を増すことにより前記海流ベクトルの推定確度を上げることを特徴とする海流推定方法。
2. 前記運航情報には、前記船舶の対地速度が特定された対地速度情報、及び前記船舶の対地方位が特定された対地方位情報が含まれ、前記対地速度及び前記対地方位から求めた前記対地速度ベクトルを用いて前記海流ベクトルを算出することを特徴とする請求項１に記載の海流推定方法。
3. 前記船首方位が特定された船首方位情報、前記対地速度情報、及び前記対地方位情報を船舶自動識別システム（Ａｕｔｏｍａｔｉｃ Ｉｄｅｎｔｉｆｉｃａｔｉｏｎ Ｓｙｓｔｅｍ：ＡＩＳ）によって取得することを特徴とする請求項２に記載の海流推定方法。
4. 前記運航情報には、前記運航情報が取得された際の前記船舶の位置が特定された位置情報が含まれ、当該位置情報に基づいて前記海域における前記海流ベクトルを算出することを特徴とする請求項１から請求項３までのいずれか１項に記載の海流推定方法。
5. 前記運航情報には、前記運航情報が取得された際の前記船舶の位置が特定された位置情報、及び前記位置情報が取得された時刻が特定された時刻情報が含まれ、前記位置情報及び前記時刻情報に基づいて前記対地速度ベクトルを算出することを特徴とする請求項１に記載の海流推定方法。
6. 前記船首方位が特定された船首方位情報、及び前記位置情報を船舶自動識別システム（Ａｕｔｏｍａｔｉｃ Ｉｄｅｎｔｉｆｉｃａｔｉｏｎ Ｓｙｓｔｅｍ：ＡＩＳ）によって取得することを特徴とする請求項５に記載の海流推定方法。
7. 前記海流ベクトルの算出に用いる情報を、前記運航情報及び/又は前記運航情報の取得状況に基づいて補正することを特徴とする請求項１から請求項６のうちの１項に記載の海流推定方法。
8. 前記海流ベクトルに関する情報をパラメータとし、各前記運航情報毎に運航情報別目的関数を作成し、
   複数の前記運航情報毎の前記運航情報別目的関数の総和である目的関数を最適化する前記海流ベクトルを算出することを特徴とする請求項１から請求項７までのいずれか１項に記載の海流推定方法。
9. 前記海流ベクトルに関する情報をパラメータとして用い各前記運航情報を用いて運航情報別目的関数を作成し、
   複数の前記運航情報毎の前記運航情報別目的関数に重み係数を乗じた値の総和である目的関数を最適化する前記海流ベクトルを算出し、前記運航情報毎の前記重み係数の設定によって前記補正を行うことを特徴とする請求項７に記載の海流推定方法。
10. 前記重み係数を、前記運航情報が取得された位置と前記海流を推定する位置との差、前記運航情報が取得された時刻と前記海流を推定する時刻との差、前記運航情報を取得した船舶の種類、あるいは前記運航情報を用いて過去に行われた海流ベクトルの算出結果、のうちの少なくとも１つに基づき設定することを特徴とする請求項９に記載の海流推定方法。
11. 請求項１から請求項１０までのいずれか１項に記載の海流推定方法に基づいて海流ベクトルを推定する海流推定システムであって、複数の前記船舶の前記運航情報を取得する情報取得部と、前記取得した前記運航情報を記憶する記憶手段と、前記海域を特定の時刻に航行する複数の前記船舶をリストアップするリストアップ手段と、リストアップされた前記船舶の前記運航情報を前記記憶手段から読み取る読取手段と、読み取られた前記運航情報に基づいて、前記海流ベクトルを算出する海流ベクトル算出部とを具備することを特徴とする海流推定システム。
12. 前記情報取得部は、予め予測された海象情報を取得し、前記海流ベクトル算出部は、推定された前記海流ベクトルを前記海象情報に基づいて補正することを特徴とする請求項１１に記載の海流推定システム。
13. 前記海流ベクトル算出部は、複数の前記海域における前記海流ベクトルを算出することにより、海流分布を推定することを特徴とする請求項１１又は請求項１２に記載の海流推定システム。
14. 前記海流ベクトル算出部は、予め予測された複数の前記海域における海流分布に対して、算出された前記海流ベクトルを同化させ、新たな海流分布を推定することを特徴とする請求項１３に記載の海流推定システム。
15. 推定された前記海流分布を元にして、航路の指示を海上の船舶に対して行う航路指示部を具備することを特徴とする請求項１３又は請求項１４に記載の海流推定システム。

Images