(57)【要約】  性能データ取得部１１１１は、複数の気象海象の環境の各々に関し、当該環境下を船舶が航行する場合の、航行速度と単位時間当たり燃料消費量との関係を示す性能データを取得する。気象海象データ取得部１１１３は、過去の複数の時期の各々に関し、当該時期における、複数の海域の各々の気象海象を示す気象海象データを取得する。算出部１１４は、性能データと気象海象データとを用いて、過去の複数の時刻の各々に関し、当該時刻を出発時刻として船舶が指定された航路に従い、指定された航行速度または単位時間当たり燃料消費量をもたらすように航行した場合に当該航行に要する時間または燃料消費量を算出する。統計処理部１１５は、算出部１１４により算出された時間または燃料消費量を標本とする母集団の確率分布を推定する。
2018069989.tif 000002
【特許請求の範囲】
【請求項１】
複数の気象海象の環境の各々に関し、当該環境下を船舶が航行する場合の、航行速度と相関する第１の変数と、単位時間当たり燃料消費量と相関する第２の変数との関係を示す性能データを取得する性能データ取得部と、
過去の複数の時期の各々に関し、当該時期における、１以上の海域の各々の気象海象を示す気象海象データを取得する気象海象データ取得部と、
前記船舶が前記１以上の海域を通過する航路を航行する場合の、前記第１の変数に関する指定された値を示す変数データを取得する変数データ取得部と、
前記性能データと前記気象海象データとを用いて、過去の複数の時刻の各々に関し、当該時刻を出発時刻として前記船舶が前記航路を前記変数データ取得部が取得した変数データが示す値をもたらすように航行した場合の前記第２の変数に関する値または当該航行に要する燃料消費量を示す値を算出する算出部と、
過去の複数の時刻の各々に関し前記算出部により算出された値を統計処理する統計処理部と
を備えるデータ処理装置。
【請求項２】
前記統計処理部は、過去の複数の時刻の各々に関し前記算出部により算出された値を標本とする母集団の確率分布を推定する
請求項１に記載のデータ処理装置。
【請求項３】
前記統計処理部により推定された確率分布を用いて、前記変数データ取得部により取得された変数データが示す値に応じた航行速度で前記船舶が前記航路を航行するために所定の確率で十分な燃料消費量を特定する燃料消費量特定部を備える
請求項２に記載のデータ処理装置。
【請求項４】
複数の気象海象の環境の各々に関し、当該環境下を船舶が航行する場合の、航行速度と相関する第１の変数と、単位時間当たり燃料消費量と相関する第２の変数との関係を示す性能データを取得する性能データ取得部と、
過去の複数の時期の各々に関し、当該時期における、１以上の海域の各々の気象海象を示す気象海象データを取得する気象海象データ取得部と、
前記船舶が前記１以上の海域を通過する航路を航行する場合の、前記第２の変数に関する指定された値を示す変数データを取得する変数データ取得部と、
前記性能データと前記気象海象データとを用いて、過去の複数の時刻の各々に関し、当該時刻を出発時刻として前記船舶が前記航路を前記変数データ取得部が取得した変数データが示す値をもたらすように航行した場合の前記第１の変数に関する値または当該航行に要する時間を示す値を算出する算出部と、
過去の複数の時刻の各々に関し前記算出部により算出された値を統計処理する統計処理部と
を備えるデータ処理装置。
【請求項５】
前記統計処理部は、過去の複数の時刻の各々に関し前記算出部により算出された値を標本とする母集団の確率分布を推定する
請求項４に記載のデータ処理装置。
【請求項６】
前記統計処理部により推定された確率分布を用いて、前記船舶が前記航路の航行を所定の確率で完了するために要する時間を特定する時間特定部を備える
請求項５に記載のデータ処理装置。
【請求項７】
前記変数データ取得部は、複数の変数データを取得し、
前記統計処理部は、前記複数の変数データの各々に関し、前記算出部により算出された値を標本とする母集団の確率分布を推定し、
前記時間特定部は、前記統計処理部により前記複数の変数データの各々に関し推定された確率分布を用いて、前記船舶が前記航路の航行を所定の確率で所定の時間内に完了するための第２の変数を特定する
請求項６に記載のデータ処理装置。
【請求項８】
前記統計処理部により推定された確率分布を用いて、前記船舶が前記航路の航行を所定の時間内に完了する確率を特定する確率特定部を備える
請求項５に記載のデータ処理装置。
【請求項９】
前記変数データ取得部は、複数の変数データを取得し、
前記統計処理部は、前記複数の変数データの各々に関し、前記算出部により算出された値を標本とする母集団の確率分布を推定し、
前記確率特定部は、前記統計処理部により前記複数の変数データの各々に関し推定された確率分布を用いて、前記船舶が前記航路の航行を所定の確率で所定の時間内に完了するための第２の変数を特定する
請求項８に記載のデータ処理装置。
【請求項１０】
前記船舶の前記航路に従う航行における積載量を示す積載量データを取得する積載量データ取得部を備え、
前記性能データ取得部は、複数の積載量の各々に関し、当該積載量の積載物を積載した前記船舶が航行する場合の前記第１の変数と前記第２の変数との関係を示す前記性能データを取得し、
前記算出部は、前記積載量データが示す積載量に関する前記性能データを用いて、前記算出を行う
請求項１乃至９のいずれか１項に記載のデータ処理装置。
【請求項１１】
コンピュータに、
複数の気象海象の環境の各々に関し、当該環境下を船舶が航行する場合の、航行速度と相関する第１の変数と、単位時間当たり燃料消費量と相関する第２の変数との関係を示す性能データを取得する処理と、
過去の複数の時期の各々に関し、当該時期における、１以上の海域の各々の気象海象を示す気象海象データを取得する処理と、
前記船舶が前記１以上の海域を通過する航路を航行する場合の、前記第１の変数に関する指定された値を示す変数データを取得する処理と、
前記性能データと前記気象海象データとを用いて、過去の複数の時刻の各々に関し、当該時刻を出発時刻として前記船舶が前記航路を前記変数データが示す値をもたらすように航行した場合の前記第２の変数に関する値または当該航行に要する燃料消費量を示す値を算出する処理と、
過去の複数の時刻の各々に関し算出した前記第２の変数に関する値または燃料消費量を示す値の統計値を算出する処理と
を実行させるためのプログラム。
【請求項１２】
コンピュータに、
複数の気象海象の環境の各々に関し、当該環境下を船舶が航行する場合の、航行速度と相関する第１の変数と、単位時間当たり燃料消費量と相関する第２の変数との関係を示す性能データを取得する処理と、
過去の複数の時期の各々に関し、当該時期における、１以上の海域の各々の気象海象を示す気象海象データを取得する処理と、
前記船舶が前記１以上の海域を通過する航路を航行する場合の、前記第１の変数に関する指定された値を示す変数データを取得する処理と、
前記性能データと前記気象海象データとを用いて、過去の複数の時刻の各々に関し、当該時刻を出発時刻として前記船舶が前記航路を前記変数データが示す値をもたらすように航行した場合の前記第２の変数に関する値または当該航行に要する燃料消費量を示す値を算出する処理と、
過去の複数の時刻の各々に関し算出した前記第２の変数に関する値または燃料消費量を示す値の統計値を算出する処理と
を実行させるためのプログラムを持続的に記録するコンピュータに読み取り可能な記録媒体。
【請求項１３】
コンピュータに、
複数の気象海象の環境の各々に関し、当該環境下を船舶が航行する場合の、航行速度と相関する第１の変数と、単位時間当たり燃料消費量と相関する第２の変数との関係を示す性能データを取得する処理と、
過去の複数の時期の各々に関し、当該時期における、１以上の海域の各々の気象海象を示す気象海象データを取得する処理と、
前記船舶が前記１以上の海域を通過する航路を航行する場合の、前記第２の変数に関する指定された値を示す変数データを取得する処理と、
前記性能データと前記気象海象データとを用いて、過去の複数の時刻の各々に関し、当該時刻を出発時刻として前記船舶が前記航路を前記変数データが示す値をもたらすように航行した場合の前記第１の変数に関する値または当該航行に要する時間を示す値を算出する処理と、
過去の複数の時刻の各々に関し算出した前記第１の変数に関する値または時間を示す値の統計値を算出する処理と
を実行させるためのプログラム。
【請求項１４】
コンピュータに、
複数の気象海象の環境の各々に関し、当該環境下を船舶が航行する場合の、航行速度と相関する第１の変数と、単位時間当たり燃料消費量と相関する第２の変数との関係を示す性能データを取得する処理と、
過去の複数の時期の各々に関し、当該時期における、１以上の海域の各々の気象海象を示す気象海象データを取得する処理と、
前記船舶が前記１以上の海域を通過する航路を航行する場合の、前記第２の変数に関する指定された値を示す変数データを取得する処理と、
前記性能データと前記気象海象データとを用いて、過去の複数の時刻の各々に関し、当該時刻を出発時刻として前記船舶が前記航路を前記変数データが示す値をもたらすように航行した場合の前記第１の変数に関する値または当該航行に要する時間を示す値を算出する処理と、
過去の複数の時刻の各々に関し算出した前記第１の変数に関する値または時間を示す値の統計値を算出する処理と
を実行させるためのプログラムを持続的に記録するコンピュータに読み取り可能な記録媒体。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、船舶の航行に要する時間または燃料消費量を推定するための技術に関する。
【背景技術】
【０００２】
船舶の航行において、航行に要する時間は、設定された到着予定時刻までに目的港に到着できるか否かを決する重要な要素である。また、航行に要する燃料消費量は、航行に要するコストを大きく左右する重要な要素である。航行に要する時間を削減するために航行速度を上げると、一般的に航行に要する燃料消費量は増加する。【０００３】
船舶の航行に要する時間および燃料消費量は、個々の船舶の性能、航行時の積載物の重量等に加え、船舶が航行中に遭遇する気象海象の環境によって左右される。船舶の航行に要する時間および燃料消費量に影響を与える気象海象の環境とは、風向、風速、波向、波高、波周期、潮向、潮速等の組み合わせである。【０００４】
船舶が航行中に遭遇する気象海象の環境を考慮して、航行に要する時間および燃料消費量の観点から最適航路を探索する技術が提案されている。例えば、特許文献１には、対象海域の気象海象の予報データを用いて対象海域内における最適航路を算出するシステムが提案されている。特許文献１に記載のシステムは、予報データに加えて、対象海域を航行している複数の船舶の各々が測定した気象海象を示す測定データを用いて最適航路を算出する。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２０１４−０１３１４５号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
既述のように、船舶の航行に要する時間および燃料消費量は船舶が航行時に遭遇する気象海象の環境によって左右される。ただし、未来の気象海象は正確に予測することはできない。また、予測する気象海象の時期が現在時刻から離れる程、予測の精度は低下する。従って、船舶の航行に要する時間および燃料消費量を正確に推定することはできない。【０００７】
特許文献１に記載のシステムは、気象海象の予報データに加え、実際に現在、対象海域を航行している船舶により実測された気象海象の測定データを用いることにより、予測データのみを用いる場合と比較して、より高い精度で未来の気象海象を予測する。しかしながら、測定データは現在の気象海象を示し、未来の気象海象を示すわけではない。従って、特許文献１に記載のシステムは、上述した、予測する気象海象の時期が現在時刻から離れる程、予測の精度が低下する、という問題を解決することはできない。【０００８】
上記の事情に鑑み、本発明は、予測困難な気象海象によって左右される船舶の航行に要する時間または燃料消費量を合理的に推定することを可能とする手段を提供する。
【課題を解決するための手段】
【０００９】
上記課題を解決するため、本発明は、複数の気象海象の環境の各々に関し、当該環境下を船舶が航行する場合の、航行速度と相関する第１の変数と、単位時間当たり燃料消費量と相関する第２の変数との関係を示す性能データを取得する性能データ取得部と、過去の複数の時期の各々に関し、当該時期における、１以上の海域の各々の気象海象を示す気象海象データを取得する気象海象データ取得部と、前記船舶が前記１以上の海域を通過する航路を航行する場合の、前記第１の変数に関する指定された値を示す変数データを取得する変数データ取得部と、前記性能データと前記気象海象データとを用いて、過去の複数の時刻の各々に関し、当該時刻を出発時刻として前記船舶が前記航路を前記変数データ取得部が取得した変数データが示す値をもたらすように航行した場合の前記第２の変数に関する値または当該航行に要する燃料消費量を示す値を算出する算出部と、過去の複数の時刻の各々に関し前記算出部により算出された値を統計処理する統計処理部とを備えるデータ処理装置を第１の態様として提供する。【００１０】
上記の第１の態様にかかるデータ処理装置において、前記統計処理部は、過去の複数の時刻の各々に関し前記算出部により算出された値を標本とする母集団の確率分布を推定する、という構成が第２の態様として採用されてもよい。【００１１】
上記の第２の態様にかかるデータ処理装置において、前記統計処理部により推定された確率分布を用いて、前記変数データ取得部により取得された変数データが示す値に応じた航行速度で前記船舶が前記航路を航行するために所定の確率で十分な燃料消費量を特定する燃料消費量特定部を備える、という構成が第３の態様として採用されてもよい。【００１２】
また、本発明は、複数の気象海象の環境の各々に関し、当該環境下を船舶が航行する場合の、航行速度と相関する第１の変数と、単位時間当たり燃料消費量と相関する第２の変数との関係を示す性能データを取得する性能データ取得部と、過去の複数の時期の各々に関し、当該時期における、１以上の海域の各々の気象海象を示す気象海象データを取得する気象海象データ取得部と、前記船舶が前記１以上の海域を通過する航路を航行する場合の、前記第２の変数に関する指定された値を示す変数データを取得する変数データ取得部と、前記性能データと前記気象海象データとを用いて、過去の複数の時刻の各々に関し、当該時刻を出発時刻として前記船舶が前記航路を前記変数データ取得部が取得した変数データが示す値をもたらすように航行した場合の前記第１の変数に関する値または当該航行に要する時間を示す値を算出する算出部と、過去の複数の時刻の各々に関し前記算出部により算出された値を統計処理する統計処理部とを備えるデータ処理装置を第４の態様として提供する。【００１３】
上記の第４の態様にかかるデータ処理装置において、前記統計処理部は、過去の複数の時刻の各々に関し前記算出部により算出された値を標本とする母集団の確率分布を推定する、という構成が第５の態様として採用されてもよい。【００１４】
上記の第５の態様にかかるデータ処理装置において、前記統計処理部により推定された確率分布を用いて、前記船舶が前記航路の航行を所定の確率で完了するために要する時間を特定する時間特定部を備える、という構成が第６の態様として採用されてもよい。【００１５】
上記の第６の態様にかかるデータ処理装置において、前記変数データ取得部は、複数の変数データを取得し、前記統計処理部は、前記複数の変数データの各々に関し、前記算出部により算出された値を標本とする母集団の確率分布を推定し、前記時間特定部は、前記統計処理部により前記複数の変数データの各々に関し推定された確率分布を用いて、前記船舶が前記航路の航行を所定の確率で所定の時間内に完了するための第２の変数を特定する、という構成が第７の態様として採用されてもよい。【００１６】
上記の第５の態様にかかるデータ処理装置において、前記統計処理部により推定された確率分布を用いて、前記船舶が前記航路の航行を所定の時間内に完了する確率を特定する確率特定部を備える、という構成が第８の態様として採用されてもよい。【００１７】
上記の第８の態様にかかるデータ処理装置において、前記変数データ取得部は、複数の変数データを取得し、前記統計処理部は、前記複数の変数データの各々に関し、前記算出部により算出された値を標本とする母集団の確率分布を推定し、前記確率特定部は、前記統計処理部により前記複数の変数データの各々に関し推定された確率分布を用いて、前記船舶が前記航路の航行を所定の確率で所定の時間内に完了するための第２の変数を特定する、という構成が第９の態様として採用されてもよい。【００１８】
上記の第１乃至第９のいずれかの態様にかかるデータ処理装置において、前記船舶の前記航路に従う航行における積載量を示す積載量データを取得する積載量データ取得部を備え、前記性能データ取得部は、複数の積載量の各々に関し、当該積載量の積載物を積載した前記船舶が航行する場合の前記第１の変数と前記第２の変数との関係を示す前記性能データを取得し、前記算出部は、前記積載量データが示す積載量に関する前記性能データを用いて、前記算出を行う、という構成が第１０の態様として採用されてもよい。【００１９】
また、本発明は、コンピュータに、複数の気象海象の環境の各々に関し、当該環境下を船舶が航行する場合の、航行速度と相関する第１の変数と、単位時間当たり燃料消費量と相関する第２の変数との関係を示す性能データを取得する処理と、過去の複数の時期の各々に関し、当該時期における、１以上の海域の各々の気象海象を示す気象海象データを取得する処理と、前記船舶が前記１以上の海域を通過する航路を航行する場合の、前記第１の変数に関する指定された値を示す変数データを取得する処理と、前記性能データと前記気象海象データとを用いて、過去の複数の時刻の各々に関し、当該時刻を出発時刻として前記船舶が前記航路を前記変数データが示す値をもたらすように航行した場合の前記第２の変数に関する値または当該航行に要する燃料消費量を示す値を算出する処理と、過去の複数の時刻の各々に関し算出した前記第２の変数に関する値または燃料消費量を示す値の統計値を算出する処理とを実行させるためのプログラムを第１１の態様として提供する。【００２０】
また、本発明は、コンピュータに、複数の気象海象の環境の各々に関し、当該環境下を船舶が航行する場合の、航行速度と相関する第１の変数と、単位時間当たり燃料消費量と相関する第２の変数との関係を示す性能データを取得する処理と、過去の複数の時期の各々に関し、当該時期における、１以上の海域の各々の気象海象を示す気象海象データを取得する処理と、前記船舶が前記１以上の海域を通過する航路を航行する場合の、前記第１の変数に関する指定された値を示す変数データを取得する処理と、前記性能データと前記気象海象データとを用いて、過去の複数の時刻の各々に関し、当該時刻を出発時刻として前記船舶が前記航路を前記変数データが示す値をもたらすように航行した場合の前記第２の変数に関する値または当該航行に要する燃料消費量を示す値を算出する処理と、過去の複数の時刻の各々に関し算出した前記第２の変数に関する値または燃料消費量を示す値の統計値を算出する処理とを実行させるためのプログラムを持続的に記録するコンピュータ読み取り可能な記録媒体を第１２の態様として提供する。【００２１】
また、本発明は、コンピュータに、複数の気象海象の環境の各々に関し、当該環境下を船舶が航行する場合の、航行速度と相関する第１の変数と、単位時間当たり燃料消費量と相関する第２の変数との関係を示す性能データを取得する処理と、過去の複数の時期の各々に関し、当該時期における、１以上の海域の各々の気象海象を示す気象海象データを取得する処理と、前記船舶が前記１以上の海域を通過する航路を航行する場合の、前記第２の変数に関する指定された値を示す変数データを取得する処理と、前記性能データと前記気象海象データとを用いて、過去の複数の時刻の各々に関し、当該時刻を出発時刻として前記船舶が前記航路を前記変数データが示す値をもたらすように航行した場合の前記第１の変数に関する値または当該航行に要する時間を示す値を算出する処理と、過去の複数の時刻の各々に関し算出した前記第１の変数に関する値または時間を示す値の統計値を算出する処理とを実行させるためのプログラムを第１３の態様として提供する。【００２２】
また、本発明は、コンピュータに、複数の気象海象の環境の各々に関し、当該環境下を船舶が航行する場合の、航行速度と相関する第１の変数と、単位時間当たり燃料消費量と相関する第２の変数との関係を示す性能データを取得する処理と、過去の複数の時期の各々に関し、当該時期における、１以上の海域の各々の気象海象を示す気象海象データを取得する処理と、前記船舶が前記１以上の海域を通過する航路を航行する場合の、前記第２の変数に関する指定された値を示す変数データを取得する処理と、前記性能データと前記気象海象データとを用いて、過去の複数の時刻の各々に関し、当該時刻を出発時刻として前記船舶が前記航路を前記変数データが示す値をもたらすように航行した場合の前記第１の変数に関する値または当該航行に要する時間を示す値を算出する処理と、過去の複数の時刻の各々に関し算出した前記第１の変数に関する値または時間を示す値の統計値を算出する処理とを実行させるためのプログラムを持続的に記録するコンピュータ読み取り可能な記録媒体を第１４の態様として提供する。
【発明の効果】
【００２３】
本発明によれば、船舶が航行する海域における過去の気象海象に基づき、その海域を船舶が航行した場合に気象海象が航行に与える影響を考慮した航行に要する航行速度または時間、もしくは、単位時間当たり燃料消費量または燃料消費量の統計値が算出される。従って、本発明により算出される統計値を用いて、予測困難な気象海象によって左右される船舶の航行に要する時間または燃料消費量の合理的な推定が可能となる。
【図面の簡単な説明】
【００２４】
【図１】一実施形態にかかるシステムの構成を示した図。
【図２】一実施形態にかかる端末装置のハードウェアとして用いられるコンピュータの構成を示した図。
【図３】一実施形態にかかるサーバ装置のハードウェアとして用いられるコンピュータ２０の構成を示した図。
【図４】一実施形態にかかる端末装置の機能構成を示した図。
【図５】一実施形態にかかる端末装置が用いる性能データ群の構成を示した図。
【図６】一実施形態にかかる性能データが示す航行速度と単位時間当たり燃料消費量との関係を示したグラフ群。
【図７】一実施形態にかかる端末装置が用いる航路テーブルの構成を示した図。
【図８】一実施形態にかかる端末装置が用いる気象海象テーブル群の構成を示した図。
【図９】一実施形態にかかる端末装置が表示する画面を示した図。
【図１０】一実施形態にかかる端末装置が用いる作業テーブルの構成を示した図。
【図１１】一実施形態にかかる端末装置が用いる標本テーブルの構成を示した図。
【図１２Ａ】一実施形態にかかる端末装置が行う処理のフローを示した図。
【図１２Ｂ】一実施形態にかかる端末装置が行う処理のフローを示した図。
【図１３】一実施形態にかかる端末装置が表示する画面を示した図。
【図１４】一実施形態にかかる端末装置が表示する画面を示した図。
【図１５】一変形例にかかる端末装置が表示する画面を示した図。
【図１６】一変形例にかかる端末装置が用いる航行速度テーブルの構成を示した図。
【図１７】一変形例にかかる端末装置が用いる燃料消費量テーブルの構成を示した図。
【図１８】一変形例にかかる端末装置が表示する画面を示した図。
【図１９】一変形例にかかる端末装置が表示する画面を示した図。
【図２０】一変形例にかかる端末装置が行う処理のフローを示した図。
【図２１】一変形例にかかる端末装置が表示する画面を示した図。
【図２２】一変形例にかかる端末装置が行う処理のフローを示した図。
【図２３】一変形例にかかる端末装置が表示する画面を示した図。
【発明を実施するための形態】
【００２５】
［実施形態］
以下に本発明の一実施形態にかかるシステム１を説明する。システム１は、ユーザに、船舶の航行に要する時間および燃料消費量に関する情報を提供するシステムである。システム１は、過去の気象海象に基づき、船舶が指定された航路に従う航行に要する時間と燃料消費量を推定し、推定した時間および燃料消費量をユーザに提示するシステムである。以下、システム１において行われる航行に要する時間と燃料消費量の推定のための処理を「推定処理」という。
【００２６】
図１はシステム１の構成を示した図である。システム１は、ネットワーク９を介してデータ通信が可能な端末装置１１とサーバ装置１２を備える。なお、図１には端末装置１１が１つのみ示されているが、端末装置１１の数はシステム１の利用者の数に応じて任意に変化する。
【００２７】
端末装置１１のハードウェアは、例えば一般的な端末装置用のコンピュータである。図２は、端末装置１１のハードウェアとして用いられるコンピュータ１０の構成を示した図である。コンピュータ１０は、各種データを記憶するメモリ１０１、メモリ１０１に記憶されているプログラムに従う各種データ処理を行うプロセッサ１０２、外部の装置との間でデータ通信を行うＩＦ（Interface）である通信ＩＦ１０３、ユーザに対し画像を表示する液晶ディスプレイ等の表示装置１０４、ユーザの操作を受け付けるキーボード等の操作装置１０５を備える。なお、コンピュータ１０に内蔵される表示装置１０４に代えて、もしくは加えて、コンピュータ１０に接続される外付けの表示装置が用いられてもよい。また、コンピュータ１０に内蔵される操作装置１０５に代えて、もしくは加えて、コンピュータ１０に接続される外付けの操作装置が用いられてもよい。
【００２８】
サーバ装置１２のハードウェアは、例えば一般的なサーバ装置用のコンピュータである。図３は、サーバ装置１２のハードウェアとして用いられるコンピュータ２０の構成を示した図である。コンピュータ２０は、各種データを記憶するメモリ２０１、メモリ２０１に記憶されているプログラムに従う各種データ処理を行うプロセッサ２０２、外部の装置との間でデータ通信を行う通信ＩＦ２０３を備える。
【００２９】
図４は、端末装置１１の機能構成を示した図である。すなわち、端末装置１１のハードウェアを構成するコンピュータ１０は、端末装置１１用のプログラムに従うデータ処理を実行することにより、図４に示す構成部を備える装置として動作する。以下、図４に示される端末装置１１の各構成部を説明する。
【００３０】
取得部１１１は、サーバ装置１２等の外部の装置から送信されてくるデータおよびユーザにより入力されたデータを取得する。外部の装置から送信されてくるデータを取得する場合、取得部１１１は通信ＩＦ２０３により実現される。また、ユーザにより入力されたデータを取得する場合、取得部１１１は操作装置１０５により実現される。
【００３１】
記憶部１１２は、取得部１１１により取得されたデータを記憶する。記憶部１１２はメモリ１０１により実現される。
【００３２】
記憶部１１２には、複数の船舶の各々の性能を示す性能データを含む性能データ群と、複数の航路の各々を示す航路データを含む航路テーブルと、複数の海域の各々における過去の気象海象を示す気象海象データを含む気象海象テーブル群が記憶されている。
【００３３】
図５は、性能データ群の構成を示した図である。性能データ群は複数の船舶の各々に関する性能データの集まりである。性能データは、船舶が航行する場合の航行速度と単位時間当たり燃料消費量の関係を示すデータである。船舶が航行する場合の航行速度と単位時間当たり燃料消費量の関係は、船舶の積載量（本実施形態において、積載量は最大積載量に対する比率（％）で示されるものとする）と、船舶が航行する際に遭遇する気象海象とによって変化する。従って、個々の船舶に関する性能データは、様々な積載量の各々に関し準備され、さらに、様々な気象海象の各々に関し準備されている。
【００３４】
図６は、或る船舶の或る積載量に関する性能データが示す航行速度と単位時間当たり燃料消費量との関係を例示したグラフ群である。図６のグラフ群の横軸は航行速度（ｋｎｏｔ）、縦軸は単位時間当たり燃料消費量（ｔｏｎ／ｄａｙ）である。
【００３５】
図６には、ビューフォート（ＢＦ）０〜８の各々に応じた９本のグラフが含まれている。ビューフォートは本来、風力の階級であるが、本実施形態においては船舶の航行に影響を与える気象海象の環境の穏やかさを示す指標として用いられ、数が小さい程、穏やかな気象海象を示す。例えば、ビューフォート１は風速１．０ｍ／ｓ、波高０．１ｍ、波周期１．２秒、等の気象海象のパラメータの組み合わせで示される。なお、風向、波向等はいずれのビューフォートに関しても、船舶の進行方向に対し正面から向かう方向が想定されている。すなわち、図６に含まれる９本のグラフの各々は、代表的な気象海象の環境下において船舶が航行を行った場合における、航行速度と単位時間当たり燃料消費量との関係を示している。
【００３６】
なお、図６はビューフォート０〜８に応じた９つの気象環境に応じた航行速度と単位時間当たり燃料消費量の関係を示しているが、実際の性能データは、より多数の様々な気象環境（例えば、風速、風向、波高、波向、波周期等の様々な組み合わせ）に応じた航行速度と単位時間当たり燃料消費量の関係を示す。
【００３７】
図７は、航路テーブルの構成を示した図である。航路テーブルは複数の航路の各々に関するデータレコード（航路データ）の集まりである。航路テーブルは、データフィールド［航路名］、［港名］、［停泊時間］、［海域名］、［距離］を有する。データフィールド［航路名］には、航路を識別する航路名が格納される。データフィールド［港名］には、航路上の複数の港（出発港および目的港を含む）の港名を示すデータが航路上の順序で格納される。データフィールド［停泊時間］には、対応する港における停泊時間を示すデータが格納される。データフィールド［海域名］には、対応する港から航路上の次の港までの航行区間において、航路が経由する海域の海域名を示すデータが、航路上の順序で格納される。データフィールド［距離］には、航路のうち対応する海域を通過する部分の距離を示すデータが格納される。
【００３８】
性能データ群および航路テーブルは、例えば、ユーザによって端末装置１１に直接入力されたデータである。なお、端末装置１１が性能データ群および航路テーブルを取得する方法はユーザによる入力に限られない。例えば、図１に図示されない装置において性能データ群および航路テーブルが生成され、端末装置１１が当該装置から性能データ群および航路テーブルを受信してもよい。
【００３９】
図８は、気象海象テーブル群の構成を示した図である。気象海象テーブル群は複数の海域の各々に応じた気象海象テーブルの集まりである。気象海象テーブルは、所定時間長の期間の各々に応じたデータレコードの集まりである。気象海象テーブルは、データフィールド［期間］、［気象海象］を有する。データフィールド［期間］には、期間を示すデータが格納される。データフィールド［気象海象］には、サブフィールド［風向］、［風速］、［波向］、［波高］、［波周期］、［潮向］、［潮速］が含まれる。これらのサブフィールドには、気象海象テーブルに応じた海域の、データフィールド［期間］に格納されているデータが示す過去の期間における、風向、風速等（サブフィールドの名称に応じた気象海象のパラメータ）を示すデータが格納される。
【００４０】
気象海象テーブル群に格納されるデータは、端末装置１１がサーバ装置１２から受信したデータである。サーバ装置１２には、時間の経過に伴い新たな気象海象データが蓄積される。端末装置１１は、例えば所定時間の経過毎にサーバ装置１２から未受信の気象海象データを受信し、気象海象テーブル群に追加する。
【００４１】
図４を参照しつつ、端末装置１１の構成の説明を続ける。取得部１１１は、性能データ群を取得する性能データ取得部１１１１、航路テーブルを取得する航路データ取得部１１１２、気象海象テーブル群に格納される気象海象データを取得する気象海象データ取得部１１１３を有する。
【００４２】
さらに、取得部１１１は、後述する設定画面においてユーザが端末装置１１に対し入力する航行速度または単位時間当たり燃料消費量を示す変数データを取得する変数データ取得部１１１４、設定画面においてユーザが端末装置１１に対し入力する積載量を示す積載量データを取得する積載量データ取得部１１１５を有する。
【００４３】
なお、取得部１１１が取得するデータは性能データ取得部１１１１〜積載量データ取得部１１１５により取得するデータに限られない。例えば、取得部１１１は、推定処理において使用される出発時刻を示すデータ等の各種データを取得する。
【００４４】
表示部１１３は、ユーザに対し各種情報を表示する。表示部１１３はプロセッサ１０２の制御下で動作する表示装置１０４により実現される。図９は、表示部１１３により表示される設定画面を示した図である。設定画面は、ユーザが、推定処理において用いられる各種パラメータを設定するための画面である。
【００４５】
設定画面の領域Ａ１には、船舶名の入力欄が表示される。ユーザが領域Ａ１の入力欄の右側のボタンをクリックすると、設定画面には船舶名の選択肢のリストが表示される。当該リストに表示される船舶名は、性能データ群（図５参照）に含まれる性能データの各々に応じた船舶名である。ユーザは、当該リストから希望する船舶名を選択することで、入力欄に船舶名を容易に入力することができる。
【００４６】
設定画面の領域Ａ２には、「航行速度一定モード」と「燃料消費量一定モード」のいずれかを選択するためのラジオボタンと、各モードに応じた変数の入力欄が表示される。航行速度一定モードは、船舶が航行速度を一定に保つように航行する航行方法である。燃料消費量一定モードは、船舶が単位時間当たり燃料消費量を一定に保つように航行する航行方法である。
【００４７】
ユーザは、領域Ａ２の航行速度一定モードに応じたラジオボタンをチェックした場合、航行速度の入力欄に航行速度を入力することができる。また、ユーザは、領域Ａ２の燃料消費量一定モードに応じたラジオボタンをチェックした場合、単位時間当たり燃料消費量の入力欄に単位時間当たり燃料消費量を入力することができる。なお、領域Ａ２の入力欄に入力されるデータが、変数データ取得部１１１４により取得される変数データである。
【００４８】
設定画面の領域Ａ３には、航路名の入力欄と世界地図が表示される。ユーザが領域Ａ３の入力欄の右側のボタンをクリックすると、設定画面には航路名の選択肢のリストが表示される。当該リストに表示される航路名は、航路テーブル（図７参照）に含まれるデータレコードの各々に応じた航路名である。ユーザは、当該リストから希望する航路名を選択することで、入力欄に航路名を容易に入力することができる。ユーザにより航路名の入力が行われると、領域Ａ３の世界地図には、入力された航路名に応じた航路が図示される。
【００４９】
設定画面の領域Ａ４には、出発港の港名の入力欄、出発時刻（出発港からの出港時刻）の入力欄、目的港の港名の入力欄、到着時刻（目的港への着港時刻）の入力欄が表示される。ユーザが出発港または目的港の入力欄の右側のボタンをクリックすると、設定画面には港名の選択肢のリストが表示される。当該リストに表示される港名は、領域Ａ３において入力された航路名に応じた航路データ（図７参照）が示す航路上の港の港名である。ユーザは、当該リストから希望する港名を選択することで、入力欄に港名を容易に入力することができる。
【００５０】
領域Ａ２において航行速度一定モードが選択されている場合、領域Ａ２の航行速度の入力欄に航行速度が入力されると、領域Ａ４の到着時刻の入力欄には自動的に到着時刻が入力される。端末装置１１は、領域Ａ２の航行速度の入力欄に航行速度が入力されると、まず、領域Ａ３において入力されている航路名に応じた航路データ（図７参照）が示す航路上の港間の距離のうち、領域Ａ４の出発港および目的港の入力欄に入力されている港名に応じた港間の区間（以下、「指定区間」という）の距離を、領域Ａ２の航行速度の入力欄に入力された航行速度で除算して、指定区間の航行に要する時間を算出する。続いて、端末装置１１は、領域Ａ４の出発時刻の入力欄に入力されている時刻から、指定区間の航行時間と、航路データが示す指定区間上の港の各々における停泊時間の合計が経過した後の時刻を、到着時刻として特定する。端末装置１１は、このように特定した到着時刻を、領域Ａ４の到着時刻の入力欄に表示する。
【００５１】
また、領域Ａ２において航行速度一定モードが選択されている場合、領域Ａ４の到着時刻の入力欄に到着時刻が入力されると、領域Ａ２の航行速度には自動的に航行速度が入力される。端末装置１１は、領域Ａ４の到着時刻の入力欄に到着時刻が入力されると、まず、領域Ａ４の到着時刻の入力欄に入力されている時刻から領域Ａ４の出発時刻の入力欄に入力されている出発時刻までの時間を全所要時間として算出する。続いて、端末装置１１は、領域Ａ３において入力されている航路名に応じた航路データ（図７参照）が示す指定区間上の港の各々における停泊時間の合計を、全所領時間から減算し、指定区間の航行時間を算出する。続いて、端末装置１１は、航路データが示す指定区間の距離を、指定区間の航行時間で除算して、航行速度を算出する。端末装置１１は、このように算出した航行速度を、領域Ａ２の航行速度の入力欄に表示する。
【００５２】
設定画面の領域Ａ５には、過去の期間を指定するための複数のチェックボックスが表示される。領域Ａ５において指定される期間は、推定処理において、航行に要する時間または燃料消費量の標本（後述）を算出するために用いられる仮想的な出発時刻の取り得る範囲を示す。以下、領域Ａ４の出発時刻の入力欄に入力された出発時刻を「指定出発時刻」と呼び、領域Ａ５において指定された期間内で標本の算出のために設定される仮想的な出発時刻を「仮想出発時刻」と呼ぶ。
【００５３】
設定画面の領域Ａ６には、積載量の入力欄が表示される。ユーザは、領域Ａ６の入力欄に積載量を入力することができる。なお、領域Ａ６の入力欄に入力されるデータが、積載量データ取得部１１１５により取得される積載量データである。
【００５４】
ユーザは、設定画面において必要なパラメータの入力を完了すると、領域Ａ７の「ＯＫ」ボタンをクリックする。この操作に応じて、端末装置１１は推定処理を行う。推定処理は、端末装置１１（図４参照）が備える算出部１１４、統計処理部１１５、燃料消費量特定部１１６、時間特定部１１７、確率特定部１１８により行われる。再び図４を参照しつつ、端末装置１１が備える構成部の説明を続ける。
【００５５】
算出部１１４は、性能データ群（図５参照）に含まれる性能データと、気象海象テーブル群（図８参照）に含まれる気象海象データとを用いて、過去の複数の時刻の各々に関し、当該時刻を出発時刻として指定された船舶が指定された航路を、指定された航行速度または単位時間当たり燃料消費量をもたらすように航行した場合に、当該航行に要する時間および燃料消費量の少なくとも一方を算出する。
【００５６】
統計処理部１１５は、過去の複数の時刻の各々に関し算出部１１４により算出された時間または燃料消費量を標本とする母集団の確率分布を推定する。
【００５７】
燃料消費量特定部１１６は、統計処理部１１５により推定された燃料消費量に関する確率分布を用いて、所定の確率で全航路の航行に十分と推定される燃料消費量、すなわち、船舶が保持すべき燃料量を特定する。
【００５８】
時間特定部１１７は、統計処理部１１５により推定された時間に関する確率分布を用いて、船舶が指定された航路に従う航行を所定の確率で完了するために要する時間を特定する。
【００５９】
確率特定部１１８は、統計処理部１１５により推定された時間に関する確率分布を用いて、船舶が指定された航路に従う航行を所定の時間内に完了する確率を特定する。
【００６０】
算出部１１４、統計処理部１１５、燃料消費量特定部１１６、時間特定部１１７、確率特定部１１８の各々が行う処理の詳細は、後述の推定処理の説明において合わせて説明する。
【００６１】
サーバ装置１２（図１参照）は、複数の海域の各々に関し、当該海域の過去の気象海象を示す気象海象データを蓄積し、外部の装置からの要求に応じて、蓄積している気象海象データの配信を行うサーバ装置である。サーバ装置１２はデータ配信を行う一般的なサーバ装置であるため、その機能構成の説明を省略する。
【００６２】
続いて、端末装置１１が行う推定処理を説明する。推定処理が行われる間、記憶部１１２には、算出部１１４により生成されるデータを一時的に格納するための作業テーブルと標本テーブルが一時的に記憶される。
【００６３】
算出部１１４は、設定画面（図９参照）の領域Ａ５においてユーザが指定した期間内の様々な仮想出発時刻の各々に関し、船舶が当該仮想出発時刻に航行を開始した場合の、航行に要する時間および燃料消費量を標本として算出する。作業テーブルは、仮想出発時刻の各々に応じて一時的に生成され、算出部１１４が１組の時間および燃料消費量の標本の算出を終えると破棄される。
【００６４】
図１０は、作業テーブルの構成を示した図である。作業テーブルは、データフィールド［港間航路］、［海域名］、［区間番号］、［距離］、［開始時刻］、［航行速度］、［単位時間当たり燃料消費量］、［終了時刻］、［燃料消費量］を有する。
【００６５】
データフィールド［港間航路］には、航路の互いに隣接する港間の部分（以下、「港間航路」という）の出発港と目的港の港名を示すデータが格納される。データフィールド［海域名］には、港間航路が通過する海域の各々を識別する海域名を示すデータが、港間航路上の順序で格納される。データフィールド［区間番号］には、港間航路のうち海域を通過する部分を出発港に近い側から所定長の距離（例えば、１０ｍｉｌｅ）毎に区分した区間の各々を識別する区間番号が格納される。データフィールド［距離］には、対応する区間の距離を示すデータが格納される。なお、各海域の最後の区間を除き、各区間の距離は所定長（例えば、１０ｍｉｌｅ）である。
【００６６】
データフィールド［開始時刻］には、対応する区間の航行の開始時刻を示すデータが格納される。データフィールド［航行速度］には、対応する区間を船舶が航行する際の航行速度が格納される。データフィールド［単位時間当たり燃料消費量］には、対応する区間を船舶が航行する際の単位時間当たり燃料消費量が格納される。データフィールド［終了時刻］には、対応する区間の航行の終了時刻を示すデータが格納される。データフィールド［燃料消費量］には、対応する区間を船舶が航行する際の燃料消費量が格納される。
【００６７】
図１１は、標本テーブルの構成を示した図である。標本テーブルは、算出部１１４が、複数の仮想出発時刻の各々に関し算出した航行に要する時間および燃料消費量の標本を格納するためのテーブルである。標本テーブルは、仮想出発時刻の各々に応じたデータレコードの集まりであり、データフィールド［仮想出発時刻］、［仮想到着時刻］、［燃料消費量］を有する。
【００６８】
データフィールド［仮想出発時刻］には、仮想出発時刻を示すデータが格納される。データフィールド［仮想到着時刻］には、船舶が仮想出発時刻に航行を開始した場合に目的港に到着する時刻に、仮想出発時刻から指定出発時刻までの時間を加えた時刻（以下、「仮想到着時刻」という）を示すデータが格納される。すなわち、仮想到着時刻は、指定出発時刻に航行を開始した船舶が、仮想出発時刻に航行を開始した場合と同様の気象海象に遭遇したと仮定した場合に、目的港に到着する時刻である。データフィールド［燃料消費量］には、船舶が仮想出発時刻に航行を開始した場合に目的港に到着するまでに要する燃料消費量を示すデータが格納される。
【００６９】
推定処理は、大きく、算出部１１４が航行に要する時間および燃料消費量の標本を多数、生成する処理と、算出部１１４により生成された標本を用いて統計処理部１１５が確率分布を推定する処理と、統計処理部１１５により推定された確率分布を用いて燃料消費量特定部１１６、時間特定部１１７、および確率特定部１１８が燃料消費量、時刻、および確率を各々特定する処理とで構成される。
【００７０】
図１２Ａおよび図１２Ｂは、算出部１１４が、航行に要する時間および燃料消費量の標本を多数、生成する処理のフローを示した図である。設定画面（図９参照）において、ユーザが領域Ａ７の「ＯＫ」ボタンをクリックすると、算出部１１４は標本テーブル（図１１参照）を生成する（ステップＳ００１）。なお、ステップＳ００１において生成される標本テーブルはまだデータレコードを含まない空のテーブルである。
【００７１】
続いて、算出部１１４は、設定画面の領域Ａ３の入力欄に入力されている航路名に応じた航路データを航路テーブル（図７参照）から読み出す（ステップＳ００２）。
【００７２】
続いて、算出部１１４は、設定画面の領域Ａ５においてチェックされているチェックボックスに応じた期間内の日から古い順に１日を選択し、領域Ａ４の入力欄に入力されている出発時刻（指定出発時刻）の日付を、選択した日付で置き換えた時刻を、仮想出発時刻として設定する（ステップＳ００３）。
【００７３】
続いて、算出部１１４は、ステップＳ００２において読み出した航路データが示す航路のうち、設定画面の領域Ａ４の入力欄に入力されている出発港から目的港に至る部分に関する作業テーブル（図１０参照）を生成する（ステップＳ００４）。ステップＳ００４において、算出部１１４は、作業テーブルのデータフィールド［海域名］、［区間番号］、［距離］に対応するデータを格納する。また、ステップＳ００４において、算出部１１４は、作業テーブルの最初の区間に応じたデータフィールド［開始時刻］に、ステップＳ００３において設定した仮想出発時刻を示すデータを格納する。すなわち、仮想出発時刻が、最初の区間の航行の開始時刻として設定される。
【００７４】
続いて、算出部１１４は、作業テーブルのデータフィールド［区間番号］に格納されている区間番号により識別される区間の中から、作業テーブルにおける上から順に１つの区間を選択する（ステップＳ００５）。以下、ステップＳ００５において選択された区間を「選択区間」という。
【００７５】
続いて、算出部１１４は、気象海象テーブル群（図８参照）から、選択区間を含む海域に応じた気象海象テーブルを読み出し、読み出した気象海象テーブルから、選択区間の航行の開始時刻を含む期間に応じた気象海象データを読み出す（ステップＳ００６）。
【００７６】
続いて、算出部１１４は、性能データ群（図５参照）から、設定画面の領域Ａ１の入力欄に入力されている船舶名、領域Ａ６の入力欄に入力されている積載量、ステップＳ００６において読み出した気象海象データが示す気象海象、の組み合わせに応じた性能データを読み出す（ステップＳ００７）。
【００７７】
続いて、算出部１１４は、設定画面の領域Ａ２において、航行速度一定モードと燃料消費量一定モードのいずれが選択されているかを判定する（ステップＳ００８）。
【００７８】
航行速度一定モードが選択されている場合（ステップＳ００８；「１」）、算出部１１４は、作業テーブルの、選択区間に応じたデータフィールド［航行速度］に、設定画面の領域Ａ２の入力欄に入力されている航行速度を示すデータを格納する。すなわち、算出部１１４は、選択区間の航行速度として、ユーザにより指定された航行速度を設定する（ステップＳ００９）。
【００７９】
続いて、算出部１１４は、ステップＳ００７において読み出した性能データが示す航行速度と単位時間当たり燃料消費量の関係に従い、ステップＳ００９において設定した航行速度に応じた単位時間当たり燃料消費量を特定する。算出部１１４は、作業テーブルの、選択区間に応じたデータフィールド［単位時間当たり燃料消費量］に、特定した単位時間当たり燃料消費量を示すデータを格納する。すなわち、算出部１１４は、選択区間の航行に要する単位時間当たり燃料消費量として、ユーザにより設定された航行速度に応じた単位時間当たり燃料消費量を設定する（ステップＳ０１０）。
【００８０】
燃料消費量一定モードが選択されている場合（ステップＳ００８；「２」）、算出部１１４は、作業テーブルの、選択区間に応じたデータフィールド［単位時間当たり燃料消費量］に、設定画面の領域Ａ２の入力欄に入力されている単位時間当たり燃料消費量を示すデータを格納する。すなわち、算出部１１４は、選択区間の航行に要する単位時間当たり燃料消費量として、ユーザにより指定された単位時間当たり燃料消費量を設定する（ステップＳ０１１）。
【００８１】
続いて、算出部１１４は、ステップＳ００７において読み出した性能データが示す航行速度と単位時間当たり燃料消費量の関係に従い、ステップＳ０１１において設定した単位時間当たり燃料消費量に応じた航行速度を特定する。算出部１１４は、作業テーブルの、選択区間に応じたデータフィールド［航行速度］に、特定した航行速度を示すデータを格納する。すなわち、算出部１１４は、選択区間の航行速度として、ユーザにより設定された単位時間当たり燃料消費量に応じた航行速度を設定する（ステップＳ０１２）。
【００８２】
ステップＳ０１０またはＳ０１２の後、算出部１１４は、選択区間の距離を選択区間の航行速度で除算し、選択区間の航行に要する時間を算出する（ステップＳ０１３）。続いて、算出部１１４は、選択区間の開始時刻に、ステップＳ０１３において算出した時間を加算して、選択区間の航行の終了時刻を特定する（ステップＳ０１４）。算出部１１４は、作業テーブルの、選択区間に応じたデータフィールド［終了時刻］に、特定した終了時刻を示すデータを格納する。
【００８３】
続いて、算出部１１４は、選択区間の単位時間当たり燃料消費量に対し、ステップＳ０１３において算出した時間を乗算して、選択区間の航行に要する燃料消費量を算出する（ステップＳ０１５）。算出部１１４は、作業テーブルの、選択区間に応じたデータフィールド［燃料消費量］に、算出した燃料消費量を示すデータを格納する。
【００８４】
続いて、算出部１１４は、選択区間が、作業テーブルの最後の区間、すなわち全航路の最後の区間であるか否かを判定する（ステップＳ０１６）。選択区間が、全航路の最後の区間でない場合（ステップＳ０１６；「Ｎｏ」）、算出部１１４は、選択区間が、港間航路の最後の区間であるか否かを判定する（ステップＳ０１７）。
【００８５】
選択区間が、港間航路の最後の区間である場合（ステップＳ０１７；「Ｙｅｓ」）、算出部１１４は、ステップＳ００２において読み出した航路データが示す、選択区間を含む港間航路の目的港における停泊時間を、選択区間の航行の終了時刻に加算して、選択区間の次の区間の航行の開始時刻を特定する。算出部１１４は、作業テーブルの、選択区間の次の区間に応じたデータフィールド［開始時刻］に、特定した開始時刻を示すデータを格納する。すなわち、算出部１１４は、港間航路の最後の区間の航行の終了時刻から、港における停泊時間が経過した時刻を次の区間の航行の開始時刻として設定する（ステップＳ０１８）。
【００８６】
選択区間が、港間航路の最後の区間でない場合（ステップＳ０１７；「Ｎｏ」）、算出部１１４は、選択区間の航行の終了時刻を示すデータを、作業テーブルの、選択区間の次の区間に応じたデータフィールド［開始時刻］に格納する。すなわち、算出部１１４は、港間航路の途中の区間の航行の終了時刻を、次の区間の航行の開始時刻として設定する（ステップＳ０１９）。
【００８７】
ステップＳ０１８またはＳ０１９の後、算出部１１４は処理をステップＳ００５に戻し、作業テーブルの次の区間を選択した後、ステップＳ００６以降の処理を繰り返す。
【００８８】
ステップＳ０１６の判定において、選択区間が、全航路の最後の区間であると判定した場合（ステップＳ０１６；「Ｙｅｓ」）、算出部１１４は、選択区間、すなわち全航路の最後の区間の航行の終了時刻に対し、ステップＳ００３において設定した仮想出発時刻から、設定画面の領域Ａ４の入力欄に入力されている出発時刻（指定出発時刻）までの時間を加算して、仮想到着時刻を特定する（ステップＳ０２０）。
【００８９】
続いて、算出部１１４は、作業テーブルのデータフィールド［燃料消費量］に格納されているデータが示す燃料消費量を全て加算して、全航路の航行に要する燃料消費量を算出する（ステップＳ０２１）。
【００９０】
続いて、算出部１１４は、標本テーブルに新しいデータレコードを１つ追加して、追加したデータレコードのデータフィールド［仮想出発時刻］にステップＳ００３において設定した仮想出発時刻を示すデータを、データフィールド［仮想到着時刻］にステップＳ０２０において特定した仮想到着時刻を示すデータを、また、データフィールド［燃料消費量］にステップＳ０２１において算出した燃料消費量を示すデータを、各々格納する（ステップＳ０２２）。
【００９１】
続いて、算出部１１４は、作業テーブルを破棄する（ステップＳ０２３）。続いて、算出部１１４は、ステップＳ００３において設定した仮想出発時刻が、設定画面の領域Ａ５においてチェックされているチェックボックスに応じた期間内の最後の日を選択して設定した仮想出発時刻、すなわち最後の仮想出発時刻であるか否かを判定する（ステップＳ０２４）。
【００９２】
ステップＳ００３において設定した仮想出発時刻が最後の仮想出発時刻でない場合（ステップＳ０２４；「Ｎｏ」）、算出部１１４は処理をステップＳ００３に戻し、設定画面の領域Ａ５においてチェックされているチェックボックスに応じた期間内の次の日を選択して新たな仮想出発時刻を設定した後、ステップＳ００４以降の処理を繰り返す。
【００９３】
ステップＳ００３において設定した仮想出発時刻が最後の仮想出発時刻である場合（ステップＳ０２４；「Ｙｅｓ」）、算出部１１４は処理を終了する。
【００９４】
上述した算出部１１４の処理が終了すると、統計処理部１１５は、記憶部１１２に記憶されている標本テーブル（図１１参照）に格納されているデータを用いて、確率分布の推定を行う。
【００９５】
統計処理部１１５が推定する確率分布は、設定画面の領域Ａ２において選択されているモードに応じて異なる。燃料消費量一定モードが選択されている場合、統計処理部１１５は、標本テーブルのデータフィールド［燃料消費量］に格納されているデータが示す燃料消費量の各々を標本とする母集団の確率分布（以下、「燃料消費量に関する確率分布」という）と、標本テーブルのデータフィールド［仮想到着時刻］に格納されているデータが示す仮想到着時刻の各々を標本とする母集団の確率分布（以下、「時間に関する確率分布」という）を推定する。一方、航行速度一定モードが選択されている場合、仮想到着時刻は一定となるため、統計処理部１１５は燃料消費量に関する確率分布のみを推定し、時間に関する確率分布の推定は行わない。
【００９６】
本実施形態において、統計処理部１１５は標本テーブルに格納されているデータが示す標本の母集団が正規分布に従うものとして燃料消費量に関する確率分布または時間に関する確率分布を推定する。ただし、統計処理部１１５が正規分布以外の確率分布を推定してもよい。なお、与えられた標本の母集団の確率分布を推定する方法は既知であるため、その説明は省略する。
【００９７】
図１３および図１４は、統計処理部１１５により推定された確率分布のグラフを含む結果表示画面を示した図である。図１３は設定画面において航行速度一定モードが選択された場合の結果表示画面を示し、図１４は設定画面において燃料消費量一定モードが選択された場合の結果表示画面を示している。
【００９８】
図１３に示される航行速度一定モードの結果表示画面の領域Ｂ１には、統計処理部１１５により推定された燃料消費量に関する確率分布のグラフが表示される。図１３に例示のグラフは、船舶が全航路の航行に要する燃料消費量は、９９．７３％の確率で１５６８．２ｔｏｎ（平均から＋３σの値）以下であることを示している。
【００９９】
燃料消費量特定部１１６は、統計処理部１１５により推定された燃料消費量に関する確率分布を用いて、９９．７３％の確率で全航路の航行に十分と推定される燃料消費量、すなわち、全航路の航行のために船舶が保持すべき燃料量を特定する。航行速度一定モードの結果表示画面の領域Ｂ１に表示されるグラフには、燃料消費量特定部１１６により特定された燃料量が表示される。
【０１００】
航行速度一定モードの結果表示画面の領域Ｂ２には、航行速度、到着時刻、燃料消費量マージンが表示される。航行速度一定モードの結果表示画面に表示される航行速度および到着時刻は、設定画面の領域Ａ２においてユーザが指定した航行速度と当該航行速度に応じて自動設定された到着時刻、もしくは、領域Ａ４においてユーザが設定した到着時刻と当該到着時刻に応じて自動設定された航行速度である。
【０１０１】
燃料消費量マージンは、仮に船舶が航行中に遭遇する気象海象が基準となるビューフォート（例えば、ビューフォート０）であった場合に、全航行に要した燃料消費量をＸ（ｔｏｎ）とし、確率分布に従い９９．７３％の確率で要すると推定される燃料消費量（図１３の例では１５６８．２ｔｏｎ）をＹ（ｔｏｎ）とした場合、（Ｙ−Ｘ）／Ｘ×１００（％）で算出される値である。
【０１０２】
統計処理部１１５は、推定した燃料消費量に関する確率分布を用いて、上記の算出式に従い、燃料消費量マージンを算出する。航行速度一定モードの結果表示画面の領域Ｂ２には、このように統計処理部１１５により算出された燃料消費量マージンが表示される。
【０１０３】
図１４に示される燃料消費量一定モードの結果表示画面の領域Ｃ１には、統計処理部１１５により推定された時間に関する確率分布のグラフが表示される。また、領域Ｃ２には、統計処理部１１５により推定された燃料消費量に関する確率分布のグラフが表示される。
【０１０４】
図１４の領域Ｃ１に例示のグラフは、船舶が９９．７３％の確率で２０１６年１０月１２日１３時１５分（平均から＋３σの値）までに目的港に到着することを示している。時間特定部１１７は、統計処理部１１５により推定された時間に関する確率分布を用いて、ユーザが設定画面の領域Ａ４において指定した出発時刻に船舶が航行を開始した場合に、９９．７３％の確率で目的港に到着している時刻（図１４の例では２０１６年１０月１２日１３時１５分）を特定する。燃料消費量一定モードの結果表示画面の領域Ｃ１に表示されるグラフには、このように時間特定部１１７により特定された時刻が表示される。
【０１０５】
図１４の領域Ｃ２に例示のグラフは、船舶が全航路の航行に要する燃料消費量は、９９．７３％の確率で１５２６．８ｔｏｎ（平均から＋３σの値）以下であることを示している。
【０１０６】
燃料消費量特定部１１６は、統計処理部１１５により推定された燃料消費量に関する確率分布を用いて、９９．７３％の確率で全航路の航行に十分であると推定される燃料消費量、すなわち、全航路の航行のために船舶が保持すべき燃料量を特定する。航行速度一定モードの結果表示画面の領域Ｃ２に表示されるグラフには、このように燃料消費量特定部１１６により特定された燃料量が表示される。
【０１０７】
燃料消費量一定モードの結果表示画面の領域Ｃ３には、単位時間当たり燃料消費量、燃料消費量マージン、到着確率が表示される。燃料消費量一定モードの結果表示画面に表示される単位時間当たり燃料消費量は、設定画面の領域Ａ２においてユーザが指定した単位時間当たり燃料消費量である。燃料消費量一定モードの結果表示画面に表示される燃料消費量マージンは、航行速度一定モードの結果表示画面に表示される燃料消費量マージンと同様に、統計処理部１１５により算出される値である。燃料消費量一定モードの結果表示画面に表示される到着確率は、設定画面の領域Ａ４においてユーザが指定した到着時刻までに目的港に到着する確率である。
【０１０８】
確率特定部１１８は、統計処理部１１５により推定された時間に関する確率分布を用いて、ユーザが設定画面の領域Ａ４において指定した出発時刻に船舶が航行を開始した場合に、ユーザが設定画面の領域Ａ４において指定した到着時刻までに船舶が目的港に到着する確率（図１４の例では６５．２％）を特定する。燃料消費量一定モードの結果表示画面の領域Ｃ３には、このように確率特定部１１８により特定された確率が到着確率として表示される。
【０１０９】
上述したように、システム１によれば、ユーザは船舶が航行を行う際に所定の確率で必要となる時間および燃料消費量を知ることができる。また、システム１によれば、ユーザは指定した出発時刻に船舶が航行を開始した場合、指定した到着時刻までに目的港に到着する確率（到着確率）を知ることができる。システム１がユーザに提示する時間、燃料消費量、到着確率は、過去の気象海象に基づき推定された確率分布に従い特定されたものである。従って、システム１がユーザに提示する時間、燃料消費量、到着確率は、気象海象が船舶の航行に与える影響が考慮された合理的な推定値である。従って、ユーザは、システム１により提供される情報を、船舶の航行における航行速度または単位時間当たり燃料消費量の決定、到着時刻が指定された目的地への船舶の航行の是非の判断等に利用することができる。
【０１１０】
［変形例］
上述した実施形態は、本発明の技術的思想の範囲内で様々に変形することができる。以下にそれらの変形の例を示す。
【０１１１】
（１）上述した実施形態において、端末装置１１は燃料消費量一定モードが指定された場合、ユーザにより指定された１つの目的地に船舶が所定の確率（例えば、９９．７３％）で到着できる時刻を到着時刻として特定し表示する。これに代えて、端末装置１１が、指定された出発港を指定された出発時刻に出発した船舶が指定された単位時間当たり燃料消費量で航行を行った場合に、複数の港の各々に関し、当該港に所定の確率（例えば、９９．７３％）で到着できる時刻を到着時刻として特定し表示する構成が採用されてもよい。
【０１１２】
図１５はこの変形例にかかる端末装置１１が表示する画面を示した図である。図１５の画面の領域Ｄ１において、ユーザは船舶名、出発港の港名、出発時刻、単位時間当たり燃料消費量、積載量を入力する。その後、ユーザが領域Ｄ２の「ＯＫ」ボタンをクリックすると、領域Ｄ３には、複数の目的港の各々に関し、９９．７３％の確率で到着できる時刻が到着時刻として表示される。この変形例において、統計処理部１１５は、出発港から目的港の各々に至る複数の航路の各々に関し、上述した実施形態における燃料消費量一定モード時の処理と同様の処理を行い、確率分布を推定する。また、時間特定部１１７は、複数の航路の各々に関し、上述した実施形態における処理と同様の処理を行い、統計処理部１１５により推定された確率分布を用いて、領域Ｄ３に表示される到着時刻を特定する。
【０１１３】
（２）上述した実施形態において、積載量データは最大積載量に対する積載量の比率（％）を示す。これに代えて、積載量データが積載物の重量（ｔｏｎ）を示してもよい。また、積載量データが船舶の喫水を示してもよい。
【０１１４】
（３）上述した実施形態において、算出部１１４が性能データに従い航行速度または単位時間当たり燃料消費量を特定する単位である区間として、港間航路のうち海域を通過する部分を出発港に近い側から所定長の距離（例えば、１０ｍｉｌｅ）毎に区分した区間が用いられる。区間の決定方法はこれに限られない。例えば、燃料消費量一定モードにおいて、船舶が所定長の時間（例えば、１５分）に航行する航路の部分が、算出部１１４が性能データに従い航行速度を特定する単位である区間として用いられてもよい。この場合、港間航路における各区間の位置および距離は、船舶が遭遇する気象海象によって変化する。従って、複数の仮想出発時刻の各々に応じて、港間航路における各区間の位置および距離は異なることになる。
【０１１５】
（４）上述した実施形態において、燃料消費量特定部１１６は、全航路の航行に要する燃料消費量を特定する。これに代えて、もしくは加えて、燃料消費量特定部１１６が、補油量を特定する構成が採用されてもよい。この変形例において、取得部１１１は、補油港に到着した時の船舶の残油量を示す残油量データを取得する。燃料消費量特定部１１６は、補油港から目的港までの航行に要する燃料消費量を特定し、特定した燃料消費量から、残油量データが示す残油量を現在して、補油量を算出する。表示部１１３は、燃料消費量特定部１１６により算出された補油量を結果表示画面に表示する。この変形例によれば、ユーザは補油港において補油すべき燃料油の量を容易に知ることができる。
【０１１６】
（５）上述した実施形態において、性能データは、航行速度と単位時間当たり燃料消費量の関係を表形式で示す。これに代えて、性能データが、航行速度と単位時間当たり燃料消費量の関係を計算式で示す構成が採用されてもよい。
【０１１７】
（６）上述した実施形態においては、船舶の停泊中の燃料消費量は考慮されない。これに代えて、端末装置１１が船舶の停泊中の燃料消費量を考慮する構成が採用されてもよい。この変形例において、取得部１１１は、停泊中の単位時間当たり燃料消費量を示すデータを取得する。算出部１１４は、経由する港の各々における停泊時間を、取得部１１１が取得したデータが示す停泊中の単位時間当たり燃料消費量に乗算して、停泊中の燃料消費量を算出する。算出部１１４は、このように算出した停泊中の燃料消費量を、航行中の燃料消費量に加算することで、停泊中の燃料消費量を考慮した、全航路の航行に要する燃料消費量を算出する。この変形例によれば、より実際に近い燃料消費量の推定が行われる。
【０１１８】
（７）上述した実施形態において、算出部１１４は１日単位で異なる仮想出発時刻を設定し、到着時刻または燃料消費量の標本を算出する。算出部１１４が設定する仮想出発時刻の時間間隔は１日に限られず、例えば、１２時間単位で異なる仮想出発時刻が設定されてもよい。
【０１１９】
（８）上述した実施形態において、ユーザが指定する変数、すなわち、航行速度（航行速度一定モード）または単位時間当たり燃料消費量（燃料消費量一定モード）は、全航路を通して一定である。これに代えて、ユーザが、例えば港間航路または海域の各々に関し異なる変数を指定可能としてもよい。
【０１２０】
（９）上述した実施形態において、船舶の航行における航行速度と単位時間当たり燃料消費量の関係に影響を与えるパラメータとして、気象海象および積載量が考慮される。船舶の航行における航行速度と単位時間当たり燃料消費量の関係に影響を与えるパラメータとして、トリム等の他のパラメータが考慮されてもよい。この場合、性能データ取得部１１１１は、トリム等のパラメータに応じた性能データを取得し、算出部１１４は指定されたトリム等のパラメータに応じた性能データを用いて、各区間の航行における航行速度または燃料消費量を算出する。
【０１２１】
（１０）上述した実施形態において、性能データは、航行速度と単位時間当たり燃料消費量の関係を示す。性能データは、航行速度と相関する第１の変数と、単位時間当たり燃料消費量と相関する第２の変数との関係を示す限り、これらの変数は限定されない。例えば、性能データが、航行速度と馬力の関係、プロペラ回転数と負荷の関係等を示してもよい。
【０１２２】
（１１）上述した実施形態において、端末装置１１が燃料消費量および到着時刻に関する推定処理を行う対象の航路は、ユーザにより指定された航路の出発港から目的港までの全航路である。これに代えて、端末装置１１が、船舶の現在位置から目的港までの航路に関する推定処理を行ってよい。この場合、端末装置１１の取得部１１１は、例えば船舶に搭載されているＧＮＳＳ（Global Navigation Satellite System）等により特定される船舶の現在位置を示す位置データを取得する。そして、端末装置１１は、取得した位置データが示す船舶の現在位置から目的港までの航路（残航）に関し、推定処理を行う。
【０１２３】
（１２）上述した実施形態において、算出部１１４は、複数の仮想出発時刻の各々に応じた、航行に要する燃料消費量（航行速度一定モードおよび燃料消費量一定モード）または到着時刻（燃料消費量一定モード）の標本を算出する。これに代えて、もしくは加えて、算出部１１４が、複数の仮想出発時刻の各々に応じた、単位時間当たり燃料消費量（航行速度一定モード）または航行速度（燃料消費量一定モード）の標本を算出してもよい。この場合、統計処理部１１５は、単位時間当たり燃料消費量の確率分布（航行速度一定モード）または航行速度（燃料消費量一定モード）の確率分布を推定することになる。
【０１２４】
（１３）上述した実施形態においては、推定処理が行われる度に作業テーブル（図１０参照）が生成され、推定処理が終了すると作業テーブルは破棄される。すなわち、算出部１１４は、各区間に関する航行速度（燃料消費量一定モード）または単位時間当たり燃料消費量（航行速度一定モード）を推定処理が行われる度に繰り返し算出する。これに代えて、算出部１１４が様々な仮想出発時刻の各々に関し、様々な単位時間当たり燃料消費量に応じた各区間の航行速度（燃料消費量一定モード）または様々な航行速度に応じた各区間の単位時間当たり燃料消費量（航行速度一定モード）を算出し、算出した結果を記憶部１１２に記憶させておいてもよい。
【０１２５】
図１６は、この変形例において算出部１１４により算出された各区間の航行速度を格納するテーブル（以下、「航行速度テーブル」という）の構成を例示した図である。
【０１２６】
図１６に例示の航行速度テーブルには、或る積載量の或る船舶が、仮想出発時刻「２０１０年１月１日１２：００」に、単位時間当たり燃料消費量「５０（ｔｏｎ／ｄａｙ）」にて或る航路に従った航行を開始した場合の、３時間で航行される区間毎に算出部１１４が算出した航行速度が格納されている。航行速度テーブルは、データフィールド［期間］、［区間の始点］、［区間の終点］、［航行速度］を備える。データフィールド［期間］には航行の期間を示すデータが格納される。この場合、データフィールド［期間］に格納されるデータが示す期間の時間長は常に３時間である。
【０１２７】
データフィールド［区間の始点］には、データフィールド［期間］に格納されるデータが示す期間に船舶が航行する区間の始点の位置を示すデータが格納される。また、データフィールド［区間の終点］には、データフィールド［期間］に格納されるデータが示す期間に船舶が航行する区間の終点の位置を示すデータが格納される。
【０１２８】
データフィールド［航行速度］には、算出部１１４が、データフィールド［期間］に格納されるデータが示す期間と、データフィールド［区間の始点］および［区間の終点］に格納されるデータが示す海域とに応じた気象海象データに基づき、性能データに従い算出した航行速度を示すデータが格納される。
【０１２９】
記憶部１１２には、図１６に示すデータ構成の航行速度テーブルが、様々な船舶、積載量、航路、仮想出発時刻、単位時間当たり燃料消費量の組み合わせに関し格納される。
【０１３０】
図１７は、この変形例において算出部１１４により算出された各区間の単位時間当たり燃料消費量を格納するテーブル（以下、「燃料消費量テーブル」という）の構成を例示した図である。
【０１３１】
図１７に例示の燃料消費量テーブルには、或る積載量の或る船舶が、仮想出発時刻「２０１０年１月１日１２：００」に、航行速度「１５（ｋｎｏｔ）」にて或る航路に従った航行を開始した場合の、３時間で航行される区間毎に算出部１１４が算出した単位時間当たり燃料消費量が格納されている。燃料消費量テーブルは、データフィールド［期間］、［区間の始点］、［区間の終点］、［単位時間当たり燃料消費量］を備える。燃料消費量テーブルのデータフィールド［期間］、［区間の始点］、［区間の終点］に格納されるデータは、航行速度テーブルの同じ名称のデータフィールドに格納されるデータと同じである。
【０１３２】
データフィールド［単位時間当たり燃料消費量］には、算出部１１４が、データフィールド［期間］に格納されるデータが示す期間と、データフィールド［区間の始点］および［区間の終点］に格納されるデータが示す海域とに応じた気象海象データに基づき、性能データに従い算出した単位時間当たり燃料消費量を示すデータが格納される。
【０１３３】
記憶部１１２には、図１７に示すデータ構成の燃料消費量テーブルが、様々な船舶、積載量、航路、仮想出発時刻、航行速度の組み合わせに関し格納される。
【０１３４】
算出部１１４は、記憶部１１２に記憶されている複数の航行速度テーブルから、設定画面（図９参照）の領域Ａ５においてユーザにより指定された期間内の仮想出発時刻に応じた航行速度テーブルを抽出する。続いて、算出部１１４は、抽出した航行速度テーブルに格納されている航行速度を用いて、航行時間および目的港への到着時刻を算出する。統計処理部１１５は、このように算出部１１４により算出された到着時刻を標本とする母集団の確率分布（時間に関する確率分布）を推定する。
【０１３５】
また、算出部１１４は、記憶部１１２に記憶されている複数の航行速度テーブルから、設定画面（図９参照）の領域Ａ５においてユーザにより指定された期間内の仮想出発時刻に応じた燃料消費量テーブルを抽出する。続いて、算出部１１４は、抽出した燃料消費量テーブルに格納されている単位時間当たり燃料消費量を用いて、全航路の航行に要する燃料消費量を算出する。統計処理部１１５は、このように算出部１１４により算出された燃料消費量を標本とする母集団の確率分布（燃料消費量に関する確率分布）を推定する。
【０１３６】
（１４）端末装置１１の統計処理部１１５が、例えば上述した変形例において記憶部１１２に記憶される航行速度テーブル（図１６参照）または燃料消費量テーブル（図１７参照）に格納されるデータを用いて、指定された単位当たり燃料消費量で船舶が航行した場合の航行速度に関する統計値、または、指定された航行速度で船舶が航行した場合の単位時間当たり燃料消費量に関する統計値を算出してもよい。
【０１３７】
図１８は、この変形例において、燃料消費量一定モードが指定されている場合に、表示部１１３が表示する画面を例示した図である。図１８に例示の画面には、ユーザにより指定された船舶が、ユーザにより指定された積載量および単位時間当たり燃料消費量で、ユーザにより指定された航路を航行する場合の、航行速度に関する統計値が表示される。
【０１３８】
図１８の左側に示されるグラフは、ユーザにより指定された条件下で船舶が航行する場合の、１月から１２月までの各月における船舶の航行速度の平均値を、ユーザにより指定された基準となるビューフォートの環境下で船舶が航行した場合の航行速度と比較した値（ｋｎｏｔ）を示している。例えば、算出部１１４は、ユーザに指定された船舶、積載量および基準となるビューフォートに応じた性能データに従い、ユーザに指定された単位時間当たり燃料消費量に対応する航行速度を基準となる航行速度Ｗとして特定する。
【０１３９】
続いて、統計処理部１１５は、ユーザに指定された船舶、航路、積載量、および単位時間当たり燃料消費量の組み合わせに応じた複数の航行速度テーブルに含まれるデータレコードの中から、１月に応じたデータレコードを抽出し、抽出したデータレコードに格納されるデータが示す航行速度の平均値Ｖを算出する。統計処理部１１５は、１月に関し算出した航行速度の平均値Ｖから、算出部１１４により算出された基準となる航行速度Ｗを減算した値（Ｖ−Ｗ）を、１月に関する航行速度マージン（ｋｎｏｔ）として算出する。統計処理部１１５は、２月〜１２月に関しても同様に航行速度マージンを算出する。図１８の左側に示されるグラフは、このように算出された航行速度マージンを示している。
【０１４０】
図１８の右側には、ユーザにより指定された１以上の月に関する航行速度の中央値、平均値が表示される。統計処理部１１５は、ユーザに指定された船舶、航路、積載量、および単位時間当たり燃料消費量の組み合わせに応じた複数の航行速度テーブルに含まれるデータレコードの中から、ユーザにより指定された月に応じたデータレコードを抽出し、抽出したデータレコードに格納されるデータが示す航行速度の中央値および平均値を算出する。図１８の右側に示される中央値、平均値は、このように算出された統計値を示している。
【０１４１】
図１９は、この変形例において、航行速度一定モードが指定されている場合に、表示部１１３が表示する画面の他の例を示した図である。図１９に例示の画面には、ユーザにより指定された船舶が、ユーザにより指定された積載量および航行速度で、ユーザにより指定された航路を航行する場合の、単位時間当たり燃料消費量に関する統計値が表示される。
【０１４２】
図１９の左側に示されるグラフは、ユーザにより指定された条件下で船舶が航行する場合の、１月から１２月までの各月における船舶の単位時間当たり燃料消費量の平均値を、ユーザにより指定された基準となるビューフォートの環境下で船舶が航行した場合の単位時間当たり燃料消費量と比較した値（％）を示している。例えば、算出部１１４は、ユーザに指定された船舶、積載量および基準となるビューフォートに応じた性能データに従い、ユーザに指定された航行速度に対応する単位時間当たり燃料消費量を基準となる単位時間当たり燃料消費量Ｍとして特定する。
【０１４３】
続いて、統計処理部１１５は、ユーザに指定された船舶、航路、積載量、および航行速度の組み合わせに応じた複数の燃料消費量テーブルに含まれるデータレコードの中から、１月に応じたデータレコードを抽出し、抽出したデータレコードに格納されるデータが示す単位時間当たり燃料消費量の平均値Ｎを算出する。統計処理部１１５は、１月に関し算出した単位時間当たり燃料消費量の平均値Ｎから算出部１１４により算出された基準となる単位時間当たり燃料消費量Ｍを減算した値（Ｎ−Ｍ）を、基準となる単位時間当たり燃料消費量Ｍで除算し１００を乗じた値（（Ｎ−Ｍ）／Ｍ×１００）を、１月に関する単位時間当たり燃料消費量マージン（％）として算出する。算出部１１４は、２月〜１２月に関しても同様に単位時間当たり燃料消費量マージンを算出する。図１９の左側に示されるグラフは、このように算出された単位時間当たり燃料消費量マージンを示している。
【０１４４】
図１９の右側には、ユーザにより指定された１以上の月に関する単位時間当たり燃料消費量の中央値、平均値が表示される。統計処理部１１５は、ユーザに指定された船舶、航路、積載量、および航行速度の組み合わせに応じた複数の燃料消費量テーブルに含まれるデータレコードの中から、ユーザにより指定された月に応じたデータレコードを抽出する。続いて、統計処理部１１５は、抽出したデータレコードに格納されるデータが示す単位時間当たり燃料消費量の各々に関し、上述した単位時間当たり燃料消費量マージンを算出する。統計処理部１１５は、このように算出した単位時間当たり燃料消費量マージンの中央値および平均値を算出する。図１９の右側に示される中央値、平均値は、このように算出された統計値を示している。
【０１４５】
（１５）端末装置１１が、ユーザにより指定された船舶、航路、積載量、航行時間、到着確率に応じた単位時間当たり燃料消費量を特定してもよい。到着確率とは、到着時刻までに船舶が航路の目的港に到着している確率である。この変形例において、ユーザは、船舶、航路、積載量、出発時刻、到着時刻、到着確率を指定する。
【０１４６】
図２０は、この変形例における端末装置１１の動作のフローを示した図である。変数データ取得部１１１４は、まず、単位時間当たり燃料消費量の初期値「０．１ｔｏｎ／ｄａｙ」を、燃料消費量一定モードにおける単位時間当たり燃料消費量として設定する（ステップＳ１０１）。
【０１４７】
続いて、算出部１１４は、ユーザにより指定された船舶、航路、積載量と、ステップＳ１０１において設定された単位時間当たり燃料消費量とに関し、上述した実施形態における燃料消費量一定モードにおける処理と同様の処理を行い、複数の仮想出発時刻の各々に応じた航行に要する時間を標本として算出する（ステップＳ１０２）。
【０１４８】
続いて、統計処理部１１５は、ステップＳ１０２において算出部１１４により算出された標本の母集団の確率分布（時間に関する確率分布）を推定する（ステップＳ１０３）。続いて、時間特定部１１７は、ステップＳ１０３において推定された確率分布に従い、ユーザにより指定された到着確率で目的港に到着している時刻を特定する（ステップＳ１０４）。
【０１４９】
続いて、記憶部１１２は、ステップＳ１０４において時間特定部１１７が特定した時刻を、当該時刻の特定に用いられた単位時間当たり燃料消費量および確率分布に対応付けて記憶する（ステップＳ１０５）。
【０１５０】
続いて、時間特定部１１７は、その時点で設定されている単位時間当たり燃料消費量が船舶の単位時間当たり燃料消費量の最大値であるか否かを判定する（ステップＳ１０６）。設定されている単位時間当たり燃料消費量が最大値でない場合（ステップＳ１０６；Ｎｏ）、変数データ取得部１１１４はその時点で設定されている単位時間当たり燃料消費量を所定量（例えば、０．１ｔｏｎ／ｄａｙ）だけ増加する（ステップＳ１０７）。その後、端末装置１１は、ステップＳ１０２以降の処理を繰り返す。
【０１５１】
設定されている単位時間当たり燃料消費量が最大値である場合（ステップＳ１０６；Ｙｅｓ）、時間特定部１１７は記憶部１１２に記憶されている多数の時刻（ステップＳ１０４において特定された時刻）のうちユーザに指定された到着時刻に最も近い時刻に対応付けて記憶されている単位時間当たり燃料消費量を特定する（ステップＳ１０８）。このように特定される単位時間当たり燃料消費量が、ユーザに指定された到着確率でユーザに指定された到着時刻までに船舶が目的港に到着するための単位時間当たり燃料消費量である。
【０１５２】
表示部１１３は、ステップＳ１０５において記憶部１１２が記憶したデータを用いて、例えば図２１に示す画面（以下、「結果表示画面」という）を表示する。結果表示画面の領域Ｅ１には、ユーザにより設定された到着確率等が表示される。領域Ｅ２には、ユーザにより指定された条件で船舶が航行するための単位時間当たり燃料消費量と燃料消費量が表示される。
【０１５３】
なお、領域Ｅ２に表示される単位時間当たり燃料消費量はステップＳ１０８において特定された単位時間当たり燃料消費量である。また、領域Ｅ２に表示される燃料消費量は、領域Ｅ２に表示される単位時間当たり燃料消費量に航行時間、すなわち、ユーザにより指定された出発時刻から到着時刻までの時間から港における停泊時間を差し引いた時間を乗じて算出される。
【０１５４】
領域Ｅ３には、単位時間当たり燃料消費量と到着時刻の関係を示すグラフが表示される。領域Ｅ３のグラフは、ステップＳ１０５において記憶部１１２が記憶したデータが示す単位時間当たり燃料消費量と時刻の関係を示している。領域Ｅ４には、領域Ｅ２に表示される単位時間当たり燃料消費量で船舶が航行した場合の到着時刻の確率分布を示すグラフが表示される。領域Ｅ４のグラフは、ステップＳ１０８において特定された単位時間当たり燃料消費量に対応付けて記憶されている確率分布を示している。
【０１５５】
上述した端末装置１１の処理においては、時間特定部１１７により特定される時刻がユーザにより指定された到着時刻と比較される。これに代えて、確率特定部１１８により特定される確率がユーザにより指定された到着確率と比較されてもよい。この場合に端末装置１１が行う処理のフローを図２２に示す。図２２に示される処理のうち、図２０に示される処理と同じものには同じステップ番号が付されている。以下、図２２に示される処理のうち、図２０に示される処理と異なる処理を説明する。
【０１５６】
図２０のフローのステップＳ１０４に代えて、確率特定部１１８は、ステップＳ１０３において推定された確率分布に従い、ユーザにより指定された到着時刻までに目的港に到着している確率を特定する（ステップＳ２０４）。
【０１５７】
続いて、記憶部１１２は、ステップＳ２０４において確率特定部１１８が特定した確率を、当該確率の特定に用いられた単位時間当たり燃料消費量および確率分布に対応付けて記憶する（ステップＳ２０５）。
【０１５８】
続いて、確率特定部１１８は、その時点で設定されている単位時間当たり燃料消費量が船舶の単位時間当たり燃料消費量の最大値であるか否かを判定する（ステップＳ１０６）。設定されている単位時間当たり燃料消費量が最大値でない場合（ステップＳ１０６；Ｎｏ）、変数データ取得部１１１４はその時点で設定されている単位時間当たり燃料消費量を所定量（例えば、０．１ｔｏｎ／ｄａｙ）だけ増加する（ステップＳ１０７）。その後、端末装置１１は、ステップＳ１０２以降の処理を繰り返す。
【０１５９】
設定されている単位時間当たり燃料消費量が最大値である場合（ステップＳ１０６；Ｙｅｓ）、確率特定部１１８は記憶部１１２に記憶されている多数の確率（ステップＳ２０４において特定された確率）のうちユーザに指定された到着確率に最も近い確率に対応付けて記憶されている単位時間当たり燃料消費量を特定する（ステップＳ２０８）。このように特定される単位時間当たり燃料消費量が、ユーザに指定された到着確率でユーザに指定された到着時刻までに船舶が目的港に到着するための単位時間当たり燃料消費量である。
【０１６０】
表示部１１３は、ステップＳ２０５において記憶部１１２が記憶したデータを用いて、例えば図２３に示す結果表示画面を表示する。図２３の結果表示画面は、図２１の結果表示画面と比較し、領域Ｅ３のグラフが異なっている。この場合、領域Ｅ３には、単位時間当たり燃料消費量と到着確率の関係を示すグラフが表示される。このグラフは、ステップＳ２０５において記憶部１１２が記憶したデータが示す単位時間当たり燃料消費量と確率の関係を示している。
【０１６１】
この変形例によれば、ユーザは到着時刻までに目的地に指定した到着確率で到着するために船舶の航行において必要となる最小限の単位時間当たり燃料消費量を知ることができる。
【０１６２】
（１６）上述した実施形態における単位時間当たり燃料消費量に代えて、負荷が用いられてもよい。この場合、端末装置１１は、例えば記憶部１１２に、単位時間当たり燃料消費量と負荷との間の関係を示す対応表または換算式を記憶し、必要に応じて単位時間当たり燃料消費量と負荷の間の変換を行う。
【０１６３】
（１７）上述した実施形態において、端末装置１１が扱う時刻は、出発時刻からの時間で表現されてもよい。例えば、算出部１１４は時間に関する標本の算出において、到着時刻に代えて、航行時間を算出してもよい。到着時刻は出発時刻から航行時間および停泊時間が経過した時刻であるため、航行時間と到着時刻のいずれが用いられてもよい。従って、時間特定部１１７が到着時刻を特定する代わりに、航行時間を特定してもよい。また、確率特定部１１８が、ユーザにより指定された到着時刻までに船舶が目的港に到着する確率を特定する代わりに、ユーザにより指定された航行時間（出発時刻から到着時刻までの時間から停泊時間を減じた時間）内に船舶が目的港に到着する確率を特定してもよい。
【０１６４】
（１８）上述した実施形態においては、端末装置１１およびサーバ装置１２は一般的なコンピュータがプログラムに従った処理を実行することにより実現される。これに代えて、端末装置１１およびサーバ装置１２の少なくとも一部が、いわゆる専用装置として構成されてもよい。
【０１６５】
（１９）本発明にかかるプログラム、すなわち、上述した実施形態において、コンピュータ１０により実行されて端末装置１１を実現するためのプログラムとして例示されたプログラムは、例えば、光記録媒体、半導体メモリなどのコンピュータ読取り可能な記録媒体に記憶された状態で提供されてもよいし、インターネット等の通信ネットワークを介して提供されてもよい。本発明にかかるプログラムが記録媒体に記憶された状態で提供される場合、コンピュータ１０が当該プログラムを記録媒体から読み取り用いる。また、本発明にかかるプログラムが通信ネットワークを介して提供される場合、コンピュータ１０が当該プログラムを配信元の装置から受信して用いる。
【符号の説明】
【０１６６】
１…システム、１０…コンピュータ、１１…端末装置、１２…サーバ装置、２０…コンピュータ、１０１…メモリ、１０２…プロセッサ、１０３…通信ＩＦ、１０４…表示装置、１０５…操作装置、１１１…取得部、１１２…記憶部、１１３…表示部、１１４…算出部、１１５…統計処理部、１１６…燃料消費量特定部、１１７…時間特定部、１１８…確率特定部、２０１…メモリ、２０２…プロセッサ、２０３…通信ＩＦ、１１１１…性能データ取得部、１１１２…航路データ取得部、１１１３…気象海象データ取得部、１１１４…変数データ取得部、１１１５…積載量データ取得部
【図１】
2018069989.tif 000003
【図２】
2018069989.tif 000004
【図３】
2018069989.tif 000005
【図４】
2018069989.tif 000006
【図５】
2018069989.tif 000007
【図６】
2018069989.tif 000008
【図７】
2018069989.tif 000009
【図８】
2018069989.tif 000010
【図９】
2018069989.tif 000011
【図１０】
2018069989.tif 000012
【図１１】
2018069989.tif 000013
【図１２Ａ】
2018069989.tif 000014
【図１２Ｂ】
2018069989.tif 000015
【図１３】
2018069989.tif 000016
【図１４】
2018069989.tif 000017
【図１５】
2018069989.tif 000018
【図１６】
2018069989.tif 000019
【図１７】
2018069989.tif 000020
【図１８】
2018069989.tif 000021
【図１９】
2018069989.tif 000022
【図２０】
2018069989.tif 000023
【図２１】
2018069989.tif 000024
【図２２】
2018069989.tif 000025
【図２３】
2018069989.tif 000026
【国際調査報告】
2018069989.tif 000027
2018069989.tif 000028
2018069989.tif 000029
2018069989.tif 000030