JP6806242B2 - 衝突リスク算出プログラム、衝突リスク算出方法および衝突リスク算出装置 - Google Patents

Description

本発明は、衝突リスク算出プログラム、衝突リスク算出方法および衝突リスク算出装置に関する。

従来より、船舶に搭載されるＡＩＳ（Automatic Identification System：船舶自動識別装置）から発信されたＡＩＳデータを用いて、船舶の航行における衝突リスクを算出する種々の手法がある。そして、種々の手法によって算出された衝突リスクをグラフなどで可視化して表示することで、船舶の航行における衝突リスクを直感的に認識させる技術がある。

特開平１１−２７２９９９号公報

しかしながら、上述の種々の手法により算出された衝突リスクのうち、幾何学的な領域として表されるものは、船舶の航行における衝突リスクを直感的に認識させる点では有用であるが、衝突リスクの度合いを定量的に把握させることができないという問題点がある。例えば、幾何学的な領域として表される衝突リスクは、船舶の操船側の視点からすると、衝突回避のためにどのように操船すればよいか直感的に把握できる。その一方、幾何学的な領域として表される衝突リスクは、航行管制側の視点からすると、衝突リスクが存在する海域の存在を把握できても、衝突リスクを定量的に把握できるものではない。

一つの側面では、幾何学的な領域として表される衝突リスクを定量的に把握できる衝突リスク算出プログラム、衝突リスク算出方法および衝突リスク算出装置を提供することを目的とする。

第１の案では、衝突リスク算出プログラムは、コンピュータに、第１の船舶および第２の船舶それぞれの位置および速度に関する進行情報を取得させる。そして、衝突リスク算出プログラムは、コンピュータに、第１の船舶および第２の船舶それぞれの進行情報から、第１の船舶および第２の船舶が将来衝突する可能性がある領域を算出させる。そして、衝突リスク算出プログラムは、コンピュータに、領域を回避するために第１の船舶または第２の船舶が取る操船量に基づく第１のリスク値を算出させる。

本発明の一の実施態様によれば、幾何学的な領域として表される衝突リスクを定量的に把握できるという効果を奏する。

図１は、実施例にかかる支援システムの概略的な構成の一例を示す図である。 図２は、実施例にかかる衝突リスク算出装置の概略的な構成の一例を示す図である。 図３は、実施例にかかるグリッド情報のデータ構成の一例を示す図である。 図４は、実施例にかかる過去航跡情報のデータ構成の一例を示す図である。 図５は、実施例にかかるグリッドごとの進入角および退出角の一例を示す図である。 図６は、実施例にかかるＯＺＴを示す危険領域情報のデータ構成の一例を示すデータの一例を示す図である。 図７は、実施例にかかるＯＺＴの表示の一例を示す図である。 図８は、実施例にかかるリスク算出結果情報の一例を示す図である。 図９は、実施例にかかるＯＺＴにおける領域型リスク値の第１の変換方法を示す図である。 図１０は、実施例にかかるＯＺＴにおける領域型リスク値の第２の変換方法を示す図である。 図１１は、実施例にかかるＯＺＴにおける領域型リスク値の第３の変換方法を示す図である。 図１２は、実施例にかかるＯＺＴにおける領域型リスク値の第３の変換方法を示す図である。 図１３は、実施例にかかる過去航跡情報集計処理の一例を示すフローチャートである。 図１４は、実施例にかかる衝突リスク算出処理の一例を示すフローチャートである。 図１５は、実施例にかかる数値型リスク値算出処理の一例を示すサブルーチンを示す図である。 図１６は、実施例にかかる領域型リスク値算出処理の一例を示すサブルーチンを示す図である。 図１７は、実施例にかかるＣＤＬにおける領域型リスク値の変換方法を示す図である。 図１８は、実施例にかかるＣＤＬの表示の一例を示す図である。 図１９は、衝突リスク算出プログラムを実行するコンピュータを示す図である。

以下に、開示の技術にかかる衝突リスク算出プログラム、衝突リスク算出方法および衝突リスク算出装置の実施例を図面に基づいて詳細に説明する。なお、以下の実施例により本発明が限定されるものではない。そして、各実施例は、処理内容を矛盾させない範囲で適宜組み合わせることが可能である。以下の実施例では、開示の技術を、船舶の航行を支援する支援システムに適用した場合を例に説明する。なお、以下の実施例において、「リスク」、「リスク値」、「衝突リスク」、「衝突リスク値」は、航行管制の対象海域を航行する船舶同士が、当該時刻における進路および速度を取り続けた場合に、所定時間後に衝突する可能性の度合いを示す。「時刻」とは、日付情報および時刻情報を含む。

（支援システムの構成）
最初に、実施例にかかる支援システム１０の一例を説明する。図１は、支援システムの概略的な構成の一例を示す図である。支援システム１０は、船舶の航行を支援する海上交通管制システムである。

図１には、２隻の船舶１１（船舶１（自船）および船舶２（相手船））、陸上施設１３が示されている。なお、船舶１（自船）および船舶２（相手船）は、相対的な関係を示すに過ぎず、船舶２（相手船）から船舶１（自船）を見た場合には、船舶２が自船であり、船舶１が相手船となる。

各船舶１１には、ＡＩＳ（Automatic Identification System）装置１２が搭載されている。例えば、特定の船舶は、法令などにより、ＡＩＳ装置１２の搭載が義務付けられている。特定の船舶は、国際航海に従事する３００総トン以上の全ての船舶、国際航海に従事する全ての旅客船、および、国際航海に従事しない５００総トン以上の全ての船舶が該当する。なお、特定の船舶以外の船舶も、ＡＩＳ装置１２を搭載してもよい。

ＡＩＳ装置１２は、搭載された船舶１１に関する各種の情報を含んだＡＩＳ情報を無線通信で周期的に送信する。ＡＩＳ情報には、例えば、緯度および経度による位置、速度、船名、時刻、船首方向、ＭＭＳＩ（Maritime Mobile Service Identity）、船長、船幅などの情報が含まれている。ＡＩＳ情報は、他の船舶１１や陸上施設１３により受信される。他の船舶１１や陸上施設１３は、受信したＡＩＳ情報の発信元の船舶１１の位置、速度、船名、時刻、船首方向、ＭＭＳＩ、船長、船幅などの各種の情報を把握できる。

陸上施設１３は、例えば、海上の船舶について監視および情報提供する役割を有する海上交通センタや港内交通管制室など、各船舶１１の航行管制を行う施設である。陸上施設１３は、各船舶１１から受信したＡＩＳ情報やレーダで検出された情報などをもとに、各船舶１１の位置を把握し、各船舶１１に対し、海上交通に関する各種の情報を提供する。

（衝突リスク算出装置の構成）
次に、実施例にかかる衝突リスク算出装置２０の構成について説明する。図２は、衝突リスク算出装置の概略的な構成の一例を示す図である。衝突リスク算出装置２０は、陸上施設１３に設けられ、船舶の航行を支援する装置である。例えば、衝突リスク算出装置２０は、陸上施設１３内に配置されたサーバなどのコンピュータ装置に実装される。

衝突リスク算出装置２０は、外部Ｉ／Ｆ（インタフェース）部２１、入力部２２、表示部２３、記憶部３０、制御部４０を有する。

外部Ｉ／Ｆ部２１は、例えば、他の装置と各種の情報を送受信するインタフェースである。外部Ｉ／Ｆ部２１は、陸上施設１３に設けられたアンテナなどの無線通信装置１３Ａを介して、各船舶１１と無線通信が可能とされており、各船舶１１と各種の情報を送受信する。例えば、外部Ｉ／Ｆ部２１は、無線通信装置１３Ａを介して、各船舶１１からＡＩＳ情報を受信する。

入力部２２は、各種の情報を入力する入力デバイスである。入力部２２としては、マウスやキーボードなどの操作の入力を受け付ける入力デバイスが挙げられる。入力部２２は、各種の情報の入力を受け付ける。例えば、入力部２２は、各種の処理の開始を指示する操作入力を受け付ける。入力部２２は、受け付けた操作内容を示す操作情報を制御部４０に入力する。

表示部２３は、各種情報を表示する表示デバイスである。表示部２３としては、ＬＣＤ（Liquid Crystal Display）やＣＲＴ（Cathode Ray Tube）などの表示デバイスが挙げられる。表示部２３は、各種情報を表示する。例えば、表示部２３は、操作画面など各種の画面を表示する。

記憶部３０は、ＨＤＤ（Hard Disk Drive）、ＳＳＤ（Solid State Drive）、光もしくは光磁気ディスクなどの外部記憶装置である。なお、記憶部３０は、ＲＡＭ（Random Access Memory）、フラッシュメモリ、ＮＶＳＲＡＭ（Non Volatile Static Random Access Memory）などのデータを書き換え可能な半導体メモリであってもよい。

記憶部３０は、制御部４０で実行されるＯＳ（Operating System）や各種プログラムを記憶する。例えば、記憶部３０は、後述する過去航跡情報集計処理や衝突リスク算出処理を実行するプログラムを記憶する。さらに、記憶部３０は、制御部４０で実行されるプログラムで用いられる各種データを記憶する。例えば、記憶部３０は、ＡＩＳ蓄積データ３１、グリッド情報３２、過去航跡情報３３、危険領域情報３４、リスク算出結果情報３５を記憶する。ＡＩＳ蓄積データ３１、グリッド情報３２、過去航跡情報３３、危険領域情報３４、リスク算出結果情報３５のそれぞれは、一例としてテーブルのデータ形式である。しかし、これに限られず、ＡＩＳ蓄積データ３１、グリッド情報３２、過去航跡情報３３、危険領域情報３４、リスク算出結果情報３５のそれぞれは、ＣＳＶ（Comma Separated Values）形式など、その他のデータ形式であってもよい。

ＡＩＳ蓄積データ３１は、各船舶１１から受信されたＡＩＳ情報を蓄積したデータである。グリッド情報３２は、陸上施設１３が航行管制の対象海域を所定サイズごとのグリッドに分割したグリッドに関する各種の情報を記憶したデータである。例えば、グリッド情報３２には、グリッドを識別する識別情報と、グリッドの領域の境界の位置の情報が記憶されている。

図３は、実施例にかかるグリッド情報のデータ構成の一例を示す図である。図３に示すように、グリッド情報３２は、「グリッドＩＤ」、「グリッド範囲」などの項目を有する。なお、図３に示したグリッド情報３２の各項目は、一例であり、その他の項目を有してもよい。

「グリッドＩＤ」は、グリッドを識別する識別情報を記憶する領域である。グリッドには、それぞれを識別する識別情報としてグリッドＩＤ（IDentifier）が付与される。「グリッドＩＤ」には、グリッドに付与されたグリッドＩＤが格納される。「グリッド範囲」は、グリッドの領域の各頂点の位置の緯度および経度を記憶する領域である。

ここで、グリッドのサイズについて説明する。船舶の場合、衝突を避けるための回避行動は、少なくとも約３０秒以上の回頭として顕れるとする。例えば、後述する衝突リスクについて、リスクの変化が読み取れる可能性が高い１０秒で評価するものとする。船舶の一般的な航行速度を約１０〜１２［ｋｎ（ノット）］とすると、１０秒間での航行距離が約５０〜６０［ｍ］となる。本実施例では、衝突リスクを安定して評価するため、衝突リスクを評価する周期で船舶の位置を求めた場合、船舶が隣接していないグリッドに位置することが避けられるようグリッドのサイズを定める。例えば、グリッドは、最小で一辺の幅が１００［ｍ］の矩形領域とする。なお、グリッドは、一辺の幅を１００〜２００［ｍ］としてもよい。また、グリッドの形状は、矩形に限られない。例えば、３角形や６角形などの多角形であってもよい。また、対象範囲を、複数の多角形の組合せによりグリッドに分割してもよい。

（過去航跡情報）
図４は、実施例にかかる過去航跡情報のデータ構成の一例を示す図である。過去航跡情報３３は、各グリッドを通過した船舶に関する各種の情報を記憶したデータである。図４に示すように、過去航跡情報３３は、「グリッドＩＤ」、「進入角」、「退出角」、「速度差分」、「船舶ＩＤ」などの項目を有する。なお、図４に示した過去航跡情報３３の各項目は、一例であり、その他の項目を有してもよい。

「グリッドＩＤ」は、グリッドのグリッドＩＤを記憶する領域である。「進入角」は、グリッドを通過した船舶のグリッドへの進入角を記憶する領域である。「退出角」は、グリッドを通過した船舶のグリッドからの退出角を記憶する領域である。「速度差分」は、船舶が各グリッドに進入した際の船舶の速度と、船舶が各グリッドから退出した際の船舶の速度との差分である。

「進入角」および「退出角」は、所定の方向を基準（０度）とした角度とする。例えば、「進入角」および「退出角」は、北の方向を基準として右回りの角度とする。なお、「退出角」は、グリッドから退出した角度であってもよく、「進入角」に対するグリッドから退出した角度の角度差を用いてもよい。「船舶ＩＤ」は、船舶を識別する識別情報を記憶する領域である。「船舶ＩＤ」には、グリッドを「進入角」、「退出角」、「速度差分」で通過した船舶のＭＭＳＩが記憶される。

（危険領域情報）
図６は、実施例にかかるＯＺＴを示す危険領域情報のデータ構成の一例を示すデータの一例を示す図である。ＯＺＴとは、Obstacle Zone by Targetの略であり、船舶衝突の危険領域を円で幾何的に示す１つの指標である。図６に示すように、危険領域情報３４は、「時刻」、「ＭＭＳＩ１」、「ＭＭＳＩ２」、「円中心点（緯度）」、「円中心点（経度）」、「円半径（ｍ）」などの項目を有する。なお、図６に示した危険領域情報３４の各項目は、一例であり、その他の項目を有してもよい。

「時刻」は、「ＭＭＳＩ１」で特定される船舶１（自船）および「ＭＭＳＩ２」で特定される船舶２（相手船）に対応する危険領域が「円中心点（緯度）」、「円中心点（経度）」、「円半径（ｍ）」となる時刻である。

図６に示す危険領域情報３４は、陸上施設１３が航行管制の対象海域内を航行する各船舶１１から受信したＡＩＳ情報に含まれるＭＭＳＩ１の船舶１およびＭＭＳＩ２の船舶２の各時刻における位置を含む。そして、危険領域情報３４は、船舶１および船舶２が各時刻における針路および速度を取り続けた場合に、所定時間後に２船が衝突する可能性が所定値より高い領域を「円中心点（緯度）」および「円中心点（経度）」を中心とする「円半径（ｍ）」の円で表す。ここで、「各時刻」とは、例えば、１秒間隔で補完された情報である。あるいは、「各時刻」とは、例えば、１秒ごとに更新された情報であってもよい。また、「各時刻から所定時間後」とは、例えば５分後である。

（ＯＺＴの表示）
図７は、実施例にかかるＯＺＴの表示の一例を示す図である。図７は、図６に示す危険領域情報３４を視覚的に表示したものである。例えば、図７に示すように、船舶１および船舶２が、現在における針路および速度を取り続けた場合に、例えば５分後に２船が衝突する可能性が閾値以上となる領域が円の危険領域ＤＺ１で表される。図７では、危険領域ＤＺ１は、１つのみを示すが、これに限られず、航行管制の対象海域を航行する船舶の数に応じて、複数となる。

（リスク算出結果情報）
図８は、実施例にかかるリスク算出結果情報の一例を示す図である。図８に示すように、リスク算出結果情報３５は、「時刻」、「船舶１についての情報」「船舶２についての情報」「リスク値１」「リスク値２」「総合リスク値」などの項目を有する。なお、図８に示したリスク算出結果情報３５の各項目は、一例であり、その他の項目を有してもよい。

「船舶１についての情報」および「船舶２についての情報」は、船舶１および船舶２それぞれから受信したＡＩＳ情報に基づく。例えば「船舶１についての情報」は、船舶１に関する「ＭＭＳＩ」、「緯度」、「経度」、「ＳＯＧ」、「ＣＯＧ」、「船長」、「船幅」などの情報を含む。

船舶１に関する「ＭＭＳＩ」は、リスク算出結果情報３５において対応付けられている「時刻」における船舶１のＭＭＳＩである。船舶１に関する「緯度」「経度」は、リスク算出結果情報３５において対応付けられている「時刻」における船舶１の位置を示す緯度情報および経度情報である。船舶１に関する「ＳＯＧ」は、リスク算出結果情報３５において対応付けられている「時刻」における船舶１の対地速力（Speed Over Ground）である。船舶１に関する「ＣＯＧ」は、リスク算出結果情報３５において対応付けられている「時刻」における船舶１の対地針路（Course Over Ground）である。船舶１に関する「船長」は、船舶１の最大長である。船舶１に関する「船幅」は、船舶１の最大幅である。なお、図８において、「船舶２についての情報」も「船舶１についての情報」と同様である。

また、リスク算出結果情報３５における「リスク値１」は、船舶１および船舶２の衝突リスクを数値形式で表す数値型リスク値である。連続的な数値形式で表現できる数値型リスク値には、例えば、ＤＣＰＡ（Distance to Closest Point of Approach）、ＴＣＰＡ（Time to Closest Point of Approach）などがある。また、連続的もしくは離散的な数値形式で表現できる数値型リスク値には、例えば、RiskLevelなどがある。リスク算出結果情報３５における「リスク値１」は、船舶１および船舶２の衝突リスクを数値形式で表す数値型リスク値のうちの１つまたは複数のリスク値を含む。例えば、図８において、１行目の「時刻」“2013/1/10”の行に対応する「リスク値１」は、“r011”である。「リスク値１」は、第２のリスク値の一例である。

また、リスク算出結果情報３５における「リスク値２」は、船舶１および船舶２の衝突リスクを、幾何学的な形式で表す領域型リスク値である。領域型リスク値には、ＯＺＴ（Obstacle Zone by Target）、ＣＤＬ（Collision Danger Line）などがある。ＯＺＴは、２つの船舶の同時存在確率に基いて円の危険領域を算出する手法である。本実施例では、「リスク値２」は、ＯＺＴであるとする。例えば、図８において、１行目の「時刻」“2013/1/10”の行に対応する「リスク値２」は、“r012”である。「リスク値２」は、第１のリスク値の一例である。

ここで、「リスク値２」は、元来は幾何学的な形式で表されたリスク値であるが、数値形式へ変換されたリスク値となっている。幾何学的な形式で表されたリスク値から数値形式で表されたリスク値へ変換する変換方法は、例えば、次の変換方法１〜３がある。本実施例では、変換方法１〜３のいずれか１、あるいは２つ以上を重み付けした結果を用いるものとする。変換方法１〜３については後述する。

なお、「リスク値１」および「リスク値２」は、取りうる値の幅が様々であるため、正規化により、値が取り得る範囲（例えば０以上１以下などの下限値および上限値）を予め揃えておく。また、リスク算出結果情報３５における「総合リスク値」は、１または複数の「リスク値１」と、１または複数の「リスク値２」とを重み付けしたリスク値である。例えば、図８において、１行目の「時刻」“2013/1/10”の行に対応する「総合リスク値」は、“R011”である。「総合リスク値」は、第３のリスク値の一例である。

なお、「リスク値１」、「リスク値２」、「総合リスク値」を算出する際の重み付けは、次の各手法を用いることができる。第１の手法は、重みを均一とする単純平均である。第２の手法は、人による主観的なリスク値の正解データがある場合における多変量解析（重回帰分析など）あるいは機械学習等の手法を用いて正解に合致させるような重みを付与するアンサンブル学習である。第３の手法は、正解データがない場合における教師なし機械学習手法（クラスタリングなど）である。

このように、幾何学的な領域（例えば衝突の危険性がある危険領域）の形式で表される領域型リスク値を数値形式で表されるリスク値へ変換した上で、数値形式で表される数値型リスク値とともに重み付けされた「総合リスク値」を、船舶の衝突リスク値とする。これにより、数値型リスク値および領域型リスク値の各特長を有する、より精度が高いリスク値を得ることができる。

（第１の変換方法）
図９は、実施例にかかるＯＺＴにおける領域型リスク値の第１の変換方法を示す図である。リスク算出部４４は、図９の速度ベクトルＡ１１で示されるように、各時刻における船舶１（自船）の将来予想航路を生成する。同様に、リスク算出部４４は、図９の速度ベクトルＡ１２で示されるように、船舶２（相手船）の将来予想航路を作成する。そして、リスク算出部４４は、船舶１（自船）および船舶２（相手船）の組合せに対して危険領域を算出する。図９に示す例では、危険領域は、危険領域ＤＺ１１、ＤＺ１２、ＤＺ１３の３つ算出される。そして、リスク算出部４４は、これらの危険領域ＤＺ１１、ＤＺ１２、ＤＺ１３を回避するために船舶１（自船）が取る、現在の針路に対する最小旋回角θ１を算出する。第１の変換方法では、最小旋回角θ１が、変換後のリスク値である。最小旋回角θ１は、最大１８０°として正規化されている。最小旋回角θ１は、遠方で船舶１（自船）を横切る船舶２（相手船）よりも、間近で船舶１（自船）を横切る船舶２（相手船）の方が、角度が大きくなり高リスクと判定されるという特長を有する。最小旋回角θ１は、操船量の一例である。なお、図９では、最小旋回角θ１は、右旋回角度であるが、左旋回角度、または、右旋回角度および左旋回角度のうち最小角となるいずれかであってもよい。

（第２の変換方法）
図１０は、実施例にかかるＯＺＴにおける領域型リスク値の第２の変換方法を示す図である。第１の変換方法同様に、リスク算出部４４は、図１０の速度ベクトルＡ２１で示されるように、各時刻における船舶１（自船）の将来予想航路を作成する。同様に、リスク算出部４４は、図１０の速度ベクトルＡ２２で示されるように、船舶２（相手船）の将来予想航路を作成する。そして、リスク算出部４４は、船舶１（自船）および船舶２（相手船）の組合せに対して危険領域を算出する。図１０に示す例では、危険領域は、危険領域ＤＺ２１、ＤＺ２２、ＤＺ２３の３つ算出される。そして、リスク算出部４４は、これらの危険領域ＤＺ２１、ＤＺ２２、ＤＺ２３を回避するために船舶１（自船）が取る、現在の針路に対する最小旋回角θ２を算出する。そして、危険領域ＤＺ２１、ＤＺ２２、ＤＺ２３までの距離を、船舶１（自船）と危険領域ＤＺ２１を結ぶ接線の接点と、船舶１（自船）までの距離Ｄ２とする。そして、船舶１（自船）の各時刻における速度で距離Ｄ２を航行するためにかかる時間Ｔ２を算出する。第２の変換方法では、最小旋回角θ２を時間Ｔ２で除算したリスク値θ２／Ｔ２が、変換後のリスク値である。リスク値θ２／Ｔ２は、幾何学的な領域をリスク値に変換する際に、危険領域ＤＺ２１、ＤＺ２２、ＤＺ２３に到達するまでの時間Ｔ２を考慮することで、船舶１（自船）の旋回操作量だけでなく、衝突回避までの時間をリスク値に反映することができる。最小旋回角θ２は、最大１８０°として正規化されている。これにより、例えば、船舶１（自船）と危険領域との位置関係が同一であっても、速度が高い方が高リスクと判定されるという特長を有する。最小旋回角θ２および船舶１（自船）の各時刻における速度は、操船量の一例である。なお、図１０では、最小旋回角θ２は、右旋回角度であるが、左旋回角度、または、右旋回角度および左旋回角度のうち最小角となるいずれかであってもよい。

（第３の変換方法）
図１１および図１２は、実施例にかかるＯＺＴにおける領域型リスク値の第３の変換方法を示す図である。リスク算出部４４は、図１１の速度ベクトルＡ３１で示されるように、各時刻における船舶１（自船）の将来予想航路を作成する。同様に、リスク算出部４４は、図１１の速度ベクトルＡ３２で示されるように、船舶２（相手船）の将来予想航路を作成する。そして、リスク算出部４４は、図１１の速度ベクトルＡ３１出で示される船舶１（自船）の針路の直線と、危険領域の両方に接する円Ｃ３の半径を船舶１（自船）が危険領域ＤＺ３１〜ＤＺ３３を回避するために操舵する際に旋回する最小旋回半径Ｒｒｅｑとして算出する。また、リスク算出部４４は、最大舵角旋回試験における縦距Ｌ３１、旋回圏（２×Ｌ３２）を算出し、旋回圏（２×Ｌ３２）の半径であるＬ３２を船舶１（自船）の旋回半径Ｒｍｉｎとする。ここで、縦距Ｌ３１は、例えば船舶１（自船）の船長の４．５倍以下である。また、旋回圏（２×Ｌ３２）は、例えば船舶１（自船）の船長の４．５倍以下である。旋回半径Ｒｍｉｎは、例えばＩＭＯ（International Maritime Organization）基準（２００２年）に準拠する。第３の変換方法では、船舶１（自船）の最小旋回半径Ｒｒｅｑを旋回半径Ｒｍｉｎで除算したリスク値Ｒｒｅｑ／Ｒｍｉｎが、変換後のリスク値である。なお、リスク値Ｒｒｅｑ／Ｒｍｉｎは、最小旋回半径Ｒｒｅｑが旋回半径Ｒｍｉｎ以上となる場合に、リスク値が最大値であると正規化されている。最小旋回角θ２および船舶１（自船）の各時刻における速度は、操船量の一例である。なお、図１１では、最小旋回半径Ｒｒｅｑは、右旋回角度であるが、左旋回角度、または、右旋回角度および左旋回角度のうち最小旋回半径をとるいずれの旋回角度であってもよい。

制御部４０は、衝突リスク算出装置２０を制御するデバイスである。制御部４０としては、ＣＰＵ（Central Processing Unit）、ＭＰＵ（Micro Processing Unit）などの処理装置や、ＡＳＩＣ（Application Specific Integrated Circuit）、ＦＰＧＡ（Field Programmable Gate Array）などの集積回路を採用できる。制御部４０は、各種の処理手順を規定したプログラムや制御データを格納するための内部メモリを有し、これらによって種々の処理を実行する。制御部４０は、各種のプログラムが動作することにより各種の処理部として機能する。例えば、制御部４０は、取得部４１、過去航跡情報算出部４２、危険領域情報算出部４３、リスク算出部４４、出力部４５を有する。

取得部４１は、各種の情報を取得する。例えば、取得部４１は、各船舶の位置および速度に関する進行情報をそれぞれ取得する。例えば、取得部４１は、各船舶の進行情報として、無線通信装置１３Ａを介して、各船舶１１からＡＩＳ情報を取得する。取得部４１は、取得したＡＩＳ情報をＡＩＳ蓄積データ３１に格納する。なお、各船舶の速度は、ＡＩＳ情報に記憶された速度を用いてもよく、各船舶の時刻ごとの位置の変化から算出してもよい。本実施例では、衝突リスク算出装置２０でＡＩＳ情報を受信する場合を説明するが、ＡＩＳ情報は、ストレージ装置など外部の記憶装置に記憶されていてもよい。この場合、取得部４１は、外部の記憶装置から、各船舶１１のＡＩＳ情報を取得する。

過去航跡情報算出部４２は、グリッドごとに、当該グリッドを通過した船舶の航行の特徴を示す情報を算出する。例えば、過去航跡情報算出部４２は、ＡＩＳ蓄積データ３１から、グリッドごとに、当該グリッドを通過した各船舶の進行方向を求める。図５は、実施例にかかるグリッドごとの進入角および退出角の一例を示す図である。例えば、過去航跡情報算出部４２は、ＡＩＳ蓄積データ３１を参照して、図５に示すように、グリッドごとに、当該グリッドを通過した各船舶の時刻ごとの位置を求め、進行方向として、各船舶のグリッドへの進入角を算出する。また、過去航跡情報算出部４２は、ＡＩＳ蓄積データ３１を参照して、グリッドごとに、当該グリッドを通過した各船舶の速度を求める。速度は、グリッドを通過した際の平均速度であってもよく、グリッドへ進入した際の速度であってもよい。

なお、各船舶からＡＩＳ情報が送信される周期が異なる場合、過去航跡情報算出部４２は、船舶ごとにＡＩＳ情報の位置や速度から補間により各時刻の位置や速度を求めてもよい。例えば、過去航跡情報算出部４２は、船舶ごとに、１秒間隔で補完された（あるいは１秒ごとに更新された）各時刻の位置や速度を算出する。そして、過去航跡情報算出部４２は、算出した各船舶１１の各時刻の位置や速度を、各船舶のＩＤ（例えばＭＭＳＩなど）に対応付けて、図４に示す過去航跡情報３３に格納する。

危険領域情報算出部４３は、航行管制の対象海域内を航行する例えば船舶１および船舶２の各時刻における位置を含み、船舶１および船舶２が各時刻における針路および速度を取り続けた場合に、所定時間後に２船が衝突する可能性が所定値より高い領域を算出する。危険領域情報算出部４３は、所定時間後に２船が衝突する可能性が所定値より高い領域を「円中心点（緯度）」「円中心点（経度）」を中心とする「円半径（ｍ）」の円で表す。危険領域情報算出部４３は、ＡＩＳ蓄積データ３１をもとに、図６に示すＯＺＴに基づく危険領域情報３４を生成する。また、危険領域情報算出部４３は、危険領域情報３４に基づいて、表示部２３を介して、図７に示すようにＯＺＴを視覚的に表示する。

リスク算出部４４は、ＡＩＳ蓄積データ３１から、「船舶１についての情報」と、「船舶２についての情報」とを取得する。そして、リスク算出部４４は、数値的に表される１または複数の「リスク値１」と、幾何的に表される１または複数のリスク値を数値へ変換した「リスク値２」とを算出する。そして、リスク算出部４４は、「リスク値１」および「リスク値２」を重み付けすることにより、「総合リスク値」を算出する。

出力部４５は、各種の出力を行う。例えば、出力部４５は、リスク算出部４４により算出された衝突リスクが閾値以上高い場合、警告を出力する。例えば、出力部４５は、衝突の危険性が高い旨を画面、衝突の危険性が高い船舶１１のＡＩＳ装置１２、外部装置に出力する。これにより、出力部４５は、衝突の危険性が高いことを通知できる。

（過去航跡情報集計処理）
図１３は、実施例にかかる過去航跡情報集計処理の一例を示すフローチャートである。過去航跡情報集計処理は、所定のタイミング、例えば、後述する衝突リスク算出処理の前のタイミングや、処理開始を指示する所定操作を受け付けたタイミングで実行される。

図１３に示すように、ステップＳ１１では、過去航跡情報算出部４２は、ＡＩＳ蓄積データ３１から、船舶ごとに、補間などにより１秒ごとに各時刻の位置や速度を算出する。ステップＳ１２では、過去航跡情報算出部４２は、陸上施設１３が航行管制の対象海域を緯度および経度に応じた領域（例えばグリッド）に分割する。

ステップＳ１３では、過去航跡情報算出部４２は、ある１つの領域について、航行した全ての船舶を探索する。ステップＳ１４では、過去航跡情報算出部４２は、当該領域を通過した各船舶の進入角、退出各、速度を算出する。ステップＳ１５では、過去航跡情報算出部４２は、領域ごとの各船舶の進入角、退出各、速度を、過去航跡情報３３に格納する。ステップＳ１６では、過去航跡情報算出部４２は、全ての領域について、ステップＳ１３〜ステップＳ１５の処理が完了したか否かを判定する。過去航跡情報算出部４２は、全ての領域について、ステップＳ１３〜ステップＳ１５の処理が完了した場合（ステップＳ１６Ｙｅｓ）、過去航跡情報集計処理を終了する。一方、過去航跡情報算出部４２は、全ての領域について、ステップＳ１３〜ステップＳ１５の処理が完了していない場合（ステップＳ１６Ｎｏ）、ステップＳ１３へ処理を移す。

（衝突リスク算出処理）
図１４は、実施例にかかる衝突リスク算出処理の一例を示すフローチャートである。衝突リスク算出処理は、所定のタイミング、例えば、所定周期もしくは処理開始を指示する所定操作を受け付けたタイミングで実行される。

図１４に示すように、ステップＳ２１では、リスク算出部４４は、ＡＩＳ蓄積データ３１から、船舶ごとに、補間などにより１秒ごとに各時刻の位置や速度を算出する。ステップＳ２２では、リスク算出部４４は、ステップＳ２２の処理による補完済みのデータを対象に、リスク算出を実行する。具体的には、ステップＳ２２に続くステップＳ２３Ａにおいて、リスク算出部４４は、連続的もしくは離散的な数値形式で表現できる数値型リスク値を算出する数値型リスク値算出モジュールを呼び出し、数値型リスク値算出処理を実行する。また、ステップＳ２２に続くステップＳ２３Ｂにおいて、リスク算出部４４は、幾何的形式で表現できる領域型リスク値を算出する領域型リスク値算出モジュールを呼び出し、領域型リスク値算出処理を実行する。

なお、図１４に示す衝突リスク算出処理のフローチャートでは、説明の簡略化のため、数値型リスク値を算出する数値型リスク値算出モジュールの呼び出しステップは、ステップＳ２３Ａの１つのみであり、１種類の数値型リスク値を算出する場合を示す。複数種類の数値型リスク値が算出される場合には、数値型リスク値ごとの数値型リスク値算出モジュールが並列に呼び出されて実行される。同様に、図１４に示す衝突リスク算出処理のフローチャートでは、領域型リスク値を算出する領域型リスク算出モジュールの呼び出しステップは、ステップＳ２３Ｂの１つのみであり、１種類の領域型リスク値を算出する場合を示す。複数種類の領域型リスク値が算出される場合には、領域型リスク値ごとの領域型リスク算出モジュールが並列に呼び出されて実行される。

ステップＳ２３Ａでは、リスク算出部４４は、算出対象の数値型リスク値を算出するための数値型リスク値算出モジュールをサブルーチンとして呼び出して、算出対象の数値型リスク値を算出する。ステップＳ２３Ａの詳細は、図１５を参照して後述する。

他方、ステップＳ２３Ｂでは、リスク算出部４４は、算出対象の領域型リスク値を算出するための領域型リスク値算出モジュールをサブルーチンとして呼び出して、算出対象の領域型リスク値を算出する。ステップＳ２３Ｂの詳細は、図１６を参照して後述する。

ステップＳ２３ＡおよびステップＳ２３Ｂに続き、ステップＳ２４では、リスク算出部４４は、ステップＳ２３Ａによる数値型リスク値算出結果を「リスク値１」とし、ステップＳ２３Ｂによる領域形式リスク値算出結果を「リスク値２」とする。そして、リスク算出部４４は、算出の基礎となった「船舶１についての情報」および「船舶２についての情報」とともに、リスク算出結果情報３５に格納する。なお、複数種類の数値型リスク値が算出される場合には、数値型リスク値ごとの「リスク値１」が算出され、複数種類の領域型リスク値が算出される場合には、領域型リスク値ごとの「リスク値２」が算出され、リスク算出結果情報３５に格納される。

ステップＳ２５では、リスク算出部４４は、ステップＳ２４でリスク算出結果情報３５に格納した「リスク値１」および「リスク値２」に対し正規化を行い、重み算出方法に従って「リスク値１」および「リスク値２」から「総合リスク値」を算出する。例えば、リスク算出部４４は、「リスク値１」および「リスク値２」それぞれを０以上１以下の値の範囲となるように正規化し、「リスク値１」および「リスク値２」の平均を「総合リスク値」として算出する。

ステップＳ２６では、リスク算出部４４は、ステップＳ２５により算出した「総合リスク値」をリスク算出結果情報３５に格納する。ステップＳ２６が終了すると、リスク算出部４４は、衝突リスク算出処理を終了する。

（数値型リスク値算出処理）
図１５は、実施例にかかる数値型リスク値算出処理の一例を示すサブルーチンを示す図である。図１５は、図１４に示す衝突リスク算出処理のステップＳ２３Ａの処理において、リスク算出部４４が、算出対象の数値型リスク値に対応する数値型リスク値算出モジュールを呼び出して実行し、「リスク値１」を算出する処理である。

図１５に示すように、ステップＳ２３Ａ−１では、リスク算出部４４は、例えば図１に示すように船舶１（自船）に対して、一定距離以内（例えば所定の近傍のグリッド）の船舶２（相手船）を抽出する。ステップＳ２３Ａ−２では、リスク算出部４４は、ステップＳ２３Ａ−１で抽出した船舶１（自船）および船舶２（相手船）の将来予想航路を、過去航跡情報３３に蓄積された情報をもとにした予測方法により算出する。例えば、船舶の針路および速度の予測を行う際には、複数の将来予測航路が考えられるが、船舶が取り得る針路について、各針路および各速度を、それぞれが生起する確率で重み付けした和の針路および速度を、船舶が取り得る将来予測航路とする。

ステップＳ２３Ａ−３では、リスク算出部４４は、定義された数値型リスク値の計算式に基づいて、ステップＳ２３Ａ−２で算出した将来予測航路から、算出対象の数値型リスク値を算出する。具体的には、リスク算出部４４は、算出対象の数値型リスク値が、例えば船舶１（自船）と船舶２（相手船）の間のリスク値を示すＤＣＰＡである場合には、船舶１（自船）に対する船舶２（相手船）のＤＣＰＡを定義する計算式に基づいて、ＤＣＰＡを算出する。ステップＳ２３Ａ−４では、リスク算出部４４は、ステップＳ２３Ａ−３で算出した数値型リスク値を、船舶１（自船）に対する船舶２（相手船）の「リスク値１」として得る。

ステップＳ２３Ａ−５では、リスク算出部４４は、ステップＳ２３Ａ−１〜２３Ａ−４の処理が、「リスク値１」を算出する対象期間の全ての時刻について完了したか否かを判定する。リスク算出部４４は、ステップＳ２３Ａ−１〜２３Ａ−４の処理が、「リスク値１」を算出する対象期間の全ての時刻について完了した場合（ステップＳ２３Ａ−５Ｙｅｓ）、ステップＳ２３Ａ−６へ処理を移す。一方、リスク算出部４４は、ステップＳ２３Ａ−１〜２３Ａ−４の処理が、「リスク値１」を算出する対象期間の全ての時刻について完了していない場合（ステップＳ２３Ａ−５Ｎｏ）、ステップＳ２３Ａ−１へ処理を移す。

ステップＳ２３Ａ−６では、リスク算出部４４は、ステップＳ２３Ａ−１〜２３Ａ−４の一連の処理を繰り返すことにより算出された「リスク値１」を、次のステップＳ２４の処理へ出力する。ステップＳ２３Ａ−６が終了すると、リスク算出部４４は、数値型リスク値算出処理を終了する。

なお、上述の数値型リスク値算出処理では、リスク算出部４４は、複数の将来予測航路を確率で重み付けして算出した将来予測航路を用いて「リスク値１」を算出するとするが、これに限られるものではない。すなわち、リスク算出部４４は、複数の将来予測航路ごとに「リスク値１」を算出し、複数の将来予測航路の各確率で重み付けして「リスク値１」を算出してもよい。これにより、より柔軟性があるリスク値を算出することが可能となる。

また、リスク算出部４４は、ステップＳ２３Ａ−２の処理において、ステップＳ２３Ａ−１で抽出した船舶１（自船）および船舶２（相手船）の将来予想航路を、過去航跡情報３３に蓄積された情報に基づかない予測方法により算出してもよい。この場合、図２に示す過去航跡情報算出部４２、図４に示す過去航跡情報３３、図１３に示す過去航跡情報集計処理は、省略できる。

（領域型リスク値算出処理）
図１６は、実施例にかかる領域型リスク値算出処理の一例を示すサブルーチンを示す図である。図１６は、図１４に示す衝突リスク算出処理のステップＳ２３Ｂの処理において、リスク算出部４４が、算出対象の領域型リスク値に対応する領域型リスク値算出モジュールを呼び出して実行し、「リスク値２」を算出する処理である。

図１６に示すように、ステップＳ２３Ｂ−１では、リスク算出部４４は、例えば図１に示すように船舶１（自船）に対して、一定距離以内（例えば所定の近傍のグリッド）の船舶２（相手船）を抽出する。ステップＳ２３Ｂ−２では、リスク算出部４４は、ステップＳ２３Ｂ−１で抽出した船舶１（自船）および船舶２（相手船）の将来予想航路を、過去航跡情報３３に蓄積された情報をもとにした予測方法により算出する。将来予想航路の算出方法は、図１５に示す数値型リスク値算出処理と同様である。

ステップＳ２３Ｂ−３では、リスク算出部４４は、危険領域情報算出部４３により、定義された危険領域（ＯＺＴ）の計算式に基づいて、ステップＳ２３Ｂ−２で算出した将来予測航路から、危険領域情報３４（図６参照）を生成する。具体的には、リスク算出部４４は、危険領域情報算出部４３によって、算出対象の領域型リスク値が領域型リスク値（ＯＺＴ）である場合に、ＯＺＴの定義に基づいて、船舶１（自船）と船舶２（相手船）のＯＺＴを算出する。ステップＳ２３Ｂ−４では、リスク算出部４４は、ステップＳ２３Ｂ−３で算出したＯＺＴを、船舶１（自船）と船舶２（相手船）の危険領域（「リスク値２」）を得る。

ステップＳ２３Ｂ−５では、領域型リスク値である危険領域を、数値型リスク値に変換するために、例えば、上述の「第１の変換方法」による最小旋回角θ１を算出する。

ステップＳ２３Ｂ−６では、リスク算出部４４は、ステップＳ２３Ｂ−１〜２３Ｂ−５の処理が、「リスク値２」を算出する対象期間の全ての時刻について完了したか否かを判定する。リスク算出部４４は、ステップＳ２３Ｂ−１〜２３Ｂ−５の処理が、「リスク値２」を算出する対象期間の全ての時刻について完了した場合（ステップＳ２３Ｂ−６Ｙｅｓ）、ステップＳ２３Ｂ−７へ処理を移す。一方、リスク算出部４４は、ステップＳ２３Ｂ−１〜２３Ｂ−５の処理が、「リスク値２」を算出する対象期間の全ての時刻について完了していない場合（ステップＳ２３Ｂ−６Ｎｏ）、ステップＳ２３Ｂ−１へ処理を移す。

ステップＳ２３Ｂ−７では、リスク算出部４４は、ステップＳ２３Ｂ−１〜２３Ｂ−５の一連の処理を繰り返すことにより算出された「リスク値２」を、次のステップＳ２４の処理へ出力する。ステップＳ２３Ｂ−７が終了すると、リスク算出部４４は、領域型リスク値算出処理を終了する。

なお、上述の領域型リスク値算出処理では、リスク算出部４４は、複数の将来予測航路を確率で重み付けして算出した将来予測航路を用いて「リスク値２」を算出するとするが、これに限られるものではない。すなわち、リスク算出部４４は、複数の将来予測航路ごとに「リスク値２」を算出し、複数の将来予測航路の各確率で重み付けして「リスク値２」を算出してもよい。これにより、より柔軟性があるリスク値を算出することが可能となる。

また、リスク算出部４４は、ステップＳ２３Ｂ−２の処理において、ステップＳ２３Ｂ−１で抽出した船舶１（自船）および船舶２（相手船）の将来予想航路を、過去航跡情報３３に蓄積された情報に基づかない予測方法により算出してもよい。この場合、数値型リスク値算出処理と同様に、図２に示す過去航跡情報算出部４２、図４に示す過去航跡情報３３、図１３に示す過去航跡情報集計処理は、省略できる。

種々の衝突リスク値は、衝突リスクごとに定義が異なることから、同一のＡＩＳデータに基づいたとしても、ある衝突リスク値では高い衝突リスクを示す一方、他の衝突リスク値では低い衝突リスクを示すといった場合が生じる。そこで、種々の複数の衝突リスクを適切に重み付けし、単一の手法に基づく単一の衝突リスクよりもより精度が高い新たな衝突リスクを算出する方法が考えられる。

しかし、複数の衝突リスクを重み付けすることは、数値形式で表される衝突リスク同士であれば可能である一方、数値形式で表される衝突リスクと幾何学的な領域（例えば衝突の危険性がある危険領域）の形式で表される衝突リスクとの間では難しい。このように、形式が異なる衝突リスクを含む衝突リスクの算出が困難であるという問題がある。

そこで、上述の実施例における衝突リスク算出装置２０は、領域型リスク値を数値形式へ変換する。また、上述の実施例における衝突リスク算出装置２０は、数値型リスク値と、領域型リスク値を数値形式へ変換したリスク値とを重み付けした総合リスク値を算出し、総合リスク値により船舶同士の衝突リスクを表す。よって、上述の実施例によれば、幾何学的な領域として表される衝突リスクを定量的に把握できる。また、上述の実施例によれば、種々の数値型リスク値および領域型リスク値それぞれの特長を併せ持った総合リスク値により、「衝突直前にリスクが最大になる」という衝突リスク値が備えるべき性質を満たしつつ、船舶衝突危険性をより高い精度で表すことができる。

また、上述の実施例では、船舶１および船舶２それぞれの将来進路を複数算出し、それぞれの生起確率で複数の将来進路を重み付けして、船舶１および船舶２それぞれの将来進路を算出する。よって、船舶１および船舶２それぞれの将来進路の予測精度を向上させ、延いては、数値型リスク値、領域型リスク値、総合リスク値の精度を向上させることができる。

また、上述の実施例では、領域型リスク値を、危険領域を回避するために船舶１または船舶２が取る旋回角に基づく値に変換する。よって、自船に他船が接近するほど旋回角に基づく値が大きくなる、すなわち、領域型リスク値の変換後の値が大きくなることから、領域型リスク値のリスクの度合いを、より適切に表すことができる。

また、上述の実施例における衝突リスク算出装置２０は、海上交通管制システムなどで船舶の衝突リスクの可視化結果を用いる際に、船舶のリアルタイムな位置を地図上に表示するとともに、各船舶が取りうる将来予測航路、各船舶の衝突リスクのグラフを表示する。よって、上述の実施例によれば、衝突リスクが高まっている船舶が発生した場合には、警報を発したり、画面上で該当の船舶をハイライト表示し、オペレータに表示部２３などで確認させ、必要に応じて船舶１１へ衝突危険性の連絡を実施させることができる。

［他の実施例］
さて、これまで開示の技術に関する実施例について説明したが、開示の技術は上述した実施例以外にも、種々の異なる形態にて実施されてよい。そこで、以下では、開示の技術に含まれる他の実施例を説明する。

（１）危険領域について
上述の実施例では、危険領域は、ＯＺＴであるが、これに限られず、ＣＤＬ（Collision Danger Line）であってもよい。図１７は、実施例にかかるＣＤＬにおける領域型リスク値の変換方法を示す図である。図１８は、実施例にかかるＣＤＬの表示の一例を示す図である。

Ｎ分以内に船舶２が速度ベクトルＡ４２の針路により領域ＤＺ４２に至る場合、船舶１が速度ベクトルＡ４２と平行かつ同じ大きさの速度ベクトルＡ４１（船舶１の針路から左旋回角度θ４１）から速度ベクトルＡ４１’（船舶１の針路から右旋回角度θ４１’）迄の針路を取る。このことにより、船舶１（自船）が危険領域ＤＺ４に進入せず、船舶１（自船）と船舶２（相手船）がＮ分以内は距離が閾値より大となり、航行の安全が保たれると見なされる。この危険領域ＤＺ４が、ＣＤＬである。そして、ＣＤＬは、船舶１（自船）と船舶２（相手船）がＮ分以内に距離が閾値以下となり、航行の安全が保たれない領域ＤＺ４１および領域ＤＺ４２を、船舶１（自船）の速度ベクトルＡ４１の先端が入らないようにするべき危険領域ＤＺ４とする手法である。なお、図１７に示すように、危険領域ＤＺ４は、円および円に接する２本の接線により構成される。

例えば、図１８に示すように、船舶１（自船）および船舶２（相手船）が、現在における針路および速度を取り続けた場合に、例えばＮ分以内に船舶１（自船）および船舶２（相手船）の距離が閾値以下になる領域が危険領域ＤＺ４である。

なお、ＯＺＴ、ＣＤＬに限られず、幾何学的な領域（例えば衝突の危険性がある危険領域）の形式で表される衝突リスク値であれば、いずれの衝突リスク値であっても開示の技術を適用することができる。

（２）時刻について
上述の実施例では、危険領域を示す時刻は、現在時刻から所定時間後、例えば５分後である。しかし、これに限らず、危険領域を示す時刻は、現在時刻であってもよい。

（３）衝突リスク算出装置２０の配置について
上述の実施例では、衝突リスク算出装置２０は、陸上施設１３内に配置されたサーバなどのコンピュータ装置に実装される。しかし、これに限られず、衝突リスク算出装置２０は、船舶１１の衝突防止を目的とした装置として、船舶１１自体へ搭載されてもよい。

（４）過去の事例解析
上述の実施例では、衝突リスク算出装置２０は、現在時刻から所定時間後の船舶同士の衝突リスク値を算出して表示などする。しかし、これに限られず、衝突リスク算出装置２０は、過去データから、事故には至らなかったが危険な事例を抽出してもよい。これにより、衝突リスク算出装置２０は、過去の事例解析や危険箇所マップ作成の支援を行い、船舶衝突の事故防止に資することができる。

また、図示した各装置の各構成要素は機能概念的なものであり、必ずしも物理的に図示の如く構成されていることを要しない。すなわち、各装置の分散・統合の具体的状態は図示のものに限られず、その全部または一部を、各種の負荷や使用状況などに応じて、任意の単位で機能的または物理的に分散・統合して構成することができる。例えば、危険領域情報算出部４３と、リスク算出部４４は、統合されてもよい。あるいは、例えば、危険領域情報算出部４３は、船舶１（自船）および船舶２（相手船）の将来予想進路を算出する将来予想進路算出部と、険領域情報を算出する険領域情報算出部とに分散されてもよい。また、リスク算出部４４は、将来進路算出部、リスク算出部とに分散されてもよい。

また、各処理部にて行なわれる各処理機能は、その全部または任意の一部が、ＣＰＵおよび該ＣＰＵにて解析実行されるプログラムにて実現され、あるいは、ワイヤードロジックによるハードウェアとして実現され得る。

（衝突リスク算出プログラム）
また、上述の実施例で説明した各種の処理は、予め用意されたプログラムをパーソナルコンピュータやワークステーションなどのコンピュータシステムで実行することによって実現することもできる。そこで、以下では、上述の実施例と同様の機能を有するプログラムを実行するコンピュータシステムの一例を説明する。図１９は、衝突リスク算出プログラムを実行するコンピュータを示す図である。

図１９に示すように、コンピュータ３００は、ＣＰＵ３１０、ＨＤＤ（Hard Disk Drive）３２０、ＲＡＭ（Random Access Memory）３４０を有する。これら３１０〜３４０の各部は、バス４００を介して接続される。

ＨＤＤ３２０には上述の実施例の各処理部と同様の機能を発揮する衝突リスク算出プログラム３２０ａが予め記憶される。例えば、上述の実施例の取得部４１、過去航跡情報算出部４２、危険領域情報算出部４３、リスク算出部４４、出力部４５と同様の機能を発揮する衝突リスク算出プログラム３２０ａを記憶させる。なお、衝突リスク算出プログラム３２０ａについては、各機能を適宜モジュール分割してもよい。

また、ＨＤＤ３２０は、各種データを記憶する。例えば、ＨＤＤ３２０は、ＯＳや各種データを記憶する。

そして、ＣＰＵ３１０が、衝突リスク算出プログラム３２０ａをＨＤＤ３２０から読み出して実行することで、実施例の各処理部と同様の動作を実行する。すなわち、衝突リスク算出プログラム３２０ａは、実施例の取得部４１、過去航跡情報算出部４２、危険領域情報算出部４３、リスク算出部４４、出力部４５と同様の動作を実行する。

なお、上述した衝突リスク算出プログラム３２０ａについては、必ずしも最初からＨＤＤ３２０に記憶させることを要しない。例えば、コンピュータ３００に挿入されるフレキシブルディスク（ＦＤ）、Compact Disk Read Only Memory（ＣＤ−ＲＯＭ）、Digital Versatile Disk（ＤＶＤ）、光磁気ディスク、ＩＣカードなどの「可搬用の物理媒体」にプログラムを記憶させておく。そして、コンピュータ３００がこれらからプログラムを読み出して実行するようにしてもよい。

さらには、公衆回線、インターネット、ＬＡＮ、ＷＡＮなどを介してコンピュータ３００に接続される「他のコンピュータ（もしくはサーバ）」などにプログラムを記憶させておく。そして、コンピュータ３００がこれらからプログラムを読み出して実行するようにしてもよい。

１ 船舶（自船）
２ 船舶（相手船）
１０ 支援システム
１１ 船舶
１２ ＡＩＳ装置
１３ 陸上施設
１３Ａ 無線通信装置
２０ 衝突リスク算出装置
２１ 外部Ｉ／Ｆ部
２２ 入力部
２３ 表示部
３０ 記憶部
３１ ＡＩＳ蓄積データ
３２ グリッド情報
３３ 過去航跡情報
３４ 危険領域情報
３５ リスク算出結果情報
４０ 制御部
４１ 取得部
４２ 過去航跡情報算出部
４３ 危険領域情報算出部
４４ リスク算出部
４５ 出力部

Claims (8)

1. コンピュータに、
   第１の船舶および第２の船舶それぞれの位置および速度に関する進行情報を取得し、
   前記第１の船舶および前記第２の船舶それぞれの前記進行情報から、前記第１の船舶および前記第２の船舶が将来衝突する可能性がある領域を算出し、
   前記領域を回避するために前記第１の船舶または前記第２の船舶が取る最小旋回角に基づく第１のリスク値を算出する
   処理を実行させることを特徴とする衝突リスク算出プログラム。
2. さらに、前記コンピュータに、
   前記第１の船舶および前記第２の船舶が将来衝突する可能性を数値的に示す第２のリスク値を算出し、
   前記第１のリスク値および前記第２のリスク値を重み付けして前記第１の船舶および前記第２の船舶が将来衝突する可能性を示す第３のリスク値を算出する
   処理を実行させることを特徴とする請求項１に記載の衝突リスク算出プログラム。
3. さらに、前記コンピュータに、
   過去に航行した船舶の進行情報から、前記第１の船舶および前記第２の船舶それぞれの複数の将来進路を算出し、それぞれの生起確率で前記複数の将来進路を重み付けして、前記第１の船舶および前記第２の船舶それぞれの将来進路を算出し、
   前記将来進路をもとに、前記第１のリスク値または前記第２のリスク値を算出する
   処理を実行させることを特徴とする請求項２に記載の衝突リスク算出プログラム。
4. 前記第１のリスク値は、前記最小旋回角である
   ことを特徴とする請求項１〜３のいずれか１項に記載の衝突リスク算出プログラム。
5. 前記第１のリスク値は、前記最小旋回角を、前記最小旋回角を取る前記第１の船舶または前記第２の船舶が前記領域との接点に至るまでにかかる時間で除算した値である
   ことを特徴とする請求項１〜３のいずれか１項に記載の衝突リスク算出プログラム。
6. 前記第１のリスク値は、前記領域を回避するために前記第１の船舶または前記第２の船舶の一方の船舶が行う操舵の際の最小旋回半径を、前記一方の船舶の最大舵角旋回試験における旋回圏の半径で除算した値である
   ことを特徴とする請求項１〜３のいずれか１項に記載の衝突リスク算出プログラム。
7. コンピュータが、
   第１の船舶および第２の船舶それぞれの位置および速度に関する進行情報を取得し、
   前記第１の船舶および前記第２の船舶それぞれの前記進行情報から、前記第１の船舶および前記第２の船舶それぞれの将来進路を算出し、
   前記将来進路から、前記第１の船舶および前記第２の船舶が将来衝突する可能性がある領域を算出し、
   前記領域を回避するために前記第１の船舶または前記第２の船舶が取る最小旋回角に基づく第１のリスク値を算出する
   処理を実行することを特徴とする衝突リスク算出方法。
8. 第１の船舶および第２の船舶それぞれの位置および速度に関する進行情報を取得する取得部と、
   前記第１の船舶および前記第２の船舶それぞれの前記進行情報から、前記第１の船舶および前記第２の船舶それぞれの将来進路を算出する将来進路算出部と、
   前記将来進路から、前記第１の船舶および前記第２の船舶が将来衝突する可能性がある領域を算出する領域情報算出部と、
   前記領域を回避するために前記第１の船舶または前記第２の船舶が取る最小旋回角に基づく第１のリスク値を算出するリスク算出部と
   を有することを特徴とする衝突リスク算出装置。

Images