JP6864483B2 - Disturbance guesser - Google Patents

Disturbance guesser Download PDF

Info

Publication number
JP6864483B2
JP6864483B2 JP2017007807A JP2017007807A JP6864483B2 JP 6864483 B2 JP6864483 B2 JP 6864483B2 JP 2017007807 A JP2017007807 A JP 2017007807A JP 2017007807 A JP2017007807 A JP 2017007807A JP 6864483 B2 JP6864483 B2 JP 6864483B2
Authority
JP
Japan
Prior art keywords
disturbance
ship
information
vector
route
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
JP2017007807A
Other languages
Japanese (ja)
Other versions
JP2018114897A (en
Inventor
謙志 柳原
謙志 柳原
啓延 植竹
啓延 植竹
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tokyo Keiki Inc
Original Assignee
Tokyo Keiki Inc
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Tokyo Keiki Inc filed Critical Tokyo Keiki Inc
Priority to JP2017007807A priority Critical patent/JP6864483B2/en
Publication of JP2018114897A publication Critical patent/JP2018114897A/en
Application granted granted Critical
Publication of JP6864483B2 publication Critical patent/JP6864483B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Description

本発明は、航路を航行する船体に対する外乱を推測する装置に関する。 The present invention relates to a device for estimating a disturbance to a hull navigating a route.

従来、船舶の航行方法として、出発地点と目的地点とを複数の直線区間により接続した航路を指定し、この航路から離れないように船体を制御する航路制御という技術が知られている。この航路制御によれば、船体は、ある区間においては設定された直線航路に沿って航行するように制御される。このため、この直線航路が風や潮流などの外乱に逆らう航路となっている場合には、航路は最短距離であるものの、外乱を考慮して航行する場合と比較して、この区間を航行するのに掛かる時間が大きくなってしまうことがある。したがって、この航路制御は、エネルギー効率の観点からは必ずしも最適ではなかった。 Conventionally, as a navigation method of a ship, there is known a technique called a route control in which a route connecting a departure point and a destination point by a plurality of straight sections is specified and the hull is controlled so as not to depart from this route. According to this route control, the hull is controlled to navigate along a set straight route in a certain section. For this reason, if this straight route is a route that opposes disturbances such as wind and tidal currents, the route is the shortest distance, but compared to the case of navigating in consideration of disturbances, this section is navigated. It may take a long time. Therefore, this route control was not always optimal from the viewpoint of energy efficiency.

このような航路制御に対して、予め測定された風や潮流の測定結果から、エネルギー効率の観点から最適な航路を計画する航路計画という技術が知られている。また、この航路計画に関連する技術として、対象船舶の機器と船体運動に関してデータ収集を行う運航モニタリングシステムと、対象船舶の実海域における性能の解析を行う就航解析システムと、船舶の航行海域の気象、海象、海流の情報を基に最適な航路を計算する最適航路計算システムとから構成される船舶の運航支援システムにおける船舶の運航支援方法、が知られている(例えば、特許文献1参照)。 For such route control, a technique called a route plan for planning an optimum route from the viewpoint of energy efficiency is known from the measurement results of wind and tidal current measured in advance. In addition, as technologies related to this route planning, an operation monitoring system that collects data on the equipment and hull movement of the target ship, an in-service analysis system that analyzes the performance of the target ship in the actual sea area, and the weather in the navigation area of the ship. , A ship operation support method in a ship operation support system including an optimum route calculation system for calculating an optimum route based on information on sea conditions and sea currents (see, for example, Patent Document 1).

特許第4970346号公報Japanese Patent No. 4970346

しかしながら、航路計画において利用される、予め計測された風や潮流は地域的に大局的なものであり、船体に対する外乱として影響する局所的な風や潮流は計測されておらず、このため、航路計画において考慮される船体に対する外乱は精度が十分ではなく、実際の外乱とは大きくことなることがある、という問題がある。 However, the pre-measured winds and tides used in the route planning are regional global, and the local winds and tides that affect the hull as a disturbance are not measured. There is a problem that the disturbance to the hull considered in the plan is not accurate enough and may differ greatly from the actual disturbance.

本発明の実施形態は、上述した問題点を解決するためになされたものであり、船体に対する外乱をより精度良く推測することができる外乱推測装置を提供することを目的とする。 An embodiment of the present invention has been made to solve the above-mentioned problems, and an object of the present invention is to provide a disturbance estimation device capable of estimating a disturbance to a hull with higher accuracy.

上述した課題を解決するため、本実施形態の外乱推測装置は、予め設定された設定航路を航行する対象船に対する外乱ベクトルを推測する外乱推測装置であって、前記対象船とは異なる少なくとも1つ以上の他船から、該他船において算出された該他船に対する外乱ベクトルと該他船の位置とを含む船舶情報を取得する取得部と、前記船舶情報に基づいて、前記設定航路上において前記対象船の前方に位置する所定地点の外乱ベクトルを推測する外乱推測部とを備える。 In order to solve the above-mentioned problems, the disturbance estimation device of the present embodiment is a disturbance estimation device that estimates a disturbance vector for a target ship navigating a preset set route, and is at least one different from the target ship. The acquisition unit that acquires ship information including the disturbance vector with respect to the other ship and the position of the other ship calculated by the other ship from the above other ships, and the said on the set route based on the ship information. It is equipped with a disturbance estimation unit that estimates the disturbance vector at a predetermined point located in front of the target ship.

本発明の実施形態によれば、船体に対する外乱をより精度良く推測することができる。 According to the embodiment of the present invention, the disturbance to the hull can be estimated more accurately.

実施形態に係る対象船及び他船を示す概略図である。It is a schematic diagram which shows the target ship and other ship which concerns on embodiment. 実施形態に係る外乱推測装置のハードウェア構成を示すブロック図である。It is a block diagram which shows the hardware configuration of the disturbance guessing apparatus which concerns on embodiment. 外乱推測装置の機能構成を示すブロック図である。It is a block diagram which shows the functional structure of a disturbance guessing apparatus. 外乱推測装置の全体処理の動作を示すフローチャートである。It is a flowchart which shows the operation of the whole processing of a disturbance guessing apparatus. 狭域外乱推測処理の動作を示すフローチャートである。It is a flowchart which shows the operation of the narrow region disturbance estimation processing. 狭域範囲内の他船を示す概略図である。It is the schematic which shows the other ship in a narrow area range. 広域外乱推測処理の動作を示すフローチャートである。It is a flowchart which shows the operation of wide area disturbance estimation processing. 広域範囲内の他船を示す概略図である。It is a schematic diagram which shows the other ship in a wide area. 交点における外乱ベクトルの算出方法を示す概略図である。It is a schematic diagram which shows the calculation method of the disturbance vector at an intersection. 航路推測処理の動作を示すフローチャートである。It is a flowchart which shows the operation of the route estimation processing. 外乱ベクトルの補間方法を示す概略図である。It is a schematic diagram which shows the interpolation method of a disturbance vector. 推測ベクトルの算出方法を示す概略図である。It is a schematic diagram which shows the calculation method of a guess vector. 航路の算出方法を示す概略図である。It is a schematic diagram which shows the calculation method of a route. 対象船の補償角を示す概略図である。It is a schematic diagram which shows the compensation angle of a target ship. 補償角により変針された対象船を示す概略図である。It is the schematic which shows the target ship which changed the needle by the compensation angle.

以下、図面を参照しながら、本発明の実施形態について説明する。 Hereinafter, embodiments of the present invention will be described with reference to the drawings.

(対象船及び外乱推測装置の構成)
まず、本実施形態に係る対象船及び他船について説明する。図1は、実施形態に係る対象船及び他船を示す概略図である。
(Configuration of target ship and disturbance estimation device)
First, the target ship and other ships according to the present embodiment will be described. FIG. 1 is a schematic view showing a target ship and another ship according to the embodiment.

図1に示すように、本実施形態において、航路上に生じる外乱を算出する対象とする対象船Tは、予め設定された航路に含まれる1区間である航路区間Sを航行しているものとし、また、航路Sの周辺には対象船Tとは異なる複数の他船O1〜O4が航行しているものとする。 As shown in FIG. 1, in the present embodiment, it is assumed that the target ship T for which the disturbance generated on the route is calculated is navigating the route section S, which is one section included in the preset route. Further, it is assumed that a plurality of other vessels O1 to O4 different from the target vessel T are navigating around the route S.

次に、本実施形態に係る外乱推測装置の構成について説明する。図2は、実施形態に係る外乱推測装置のハードウェア構成を示すブロック図である。図3は、外乱推測装置の機能構成を示すブロック図である。 Next, the configuration of the disturbance estimation device according to the present embodiment will be described. FIG. 2 is a block diagram showing a hardware configuration of the disturbance estimation device according to the embodiment. FIG. 3 is a block diagram showing a functional configuration of the disturbance estimation device.

外乱推測装置10は、対象船T及び他船O1〜O4にそれぞれ備えられるものとするが、以下の説明においては、対象船Tに備えられる外乱推測装置10について説明する。なお、本実施形態において、対象船Tは自船であるが、外乱推測装置10は、船舶に備えられていなくとも良く、対象船T及び他船O1〜O4と通信可能となっていれば、例えば、陸地に設置されるサーバであっても良い。また、対象船T及び他船O1〜O4のそれぞれは、船体のsurge速度を検出するスピードログ、船体の船首方位を検出するジャイロコンパス、GPS等の衛星測位システム(GNSS)からの船体位置を検出するGNSSセンサを含むセンサ類を備えるものとする。 The disturbance estimation device 10 is provided on the target ship T and the other ships O1 to O4, respectively. In the following description, the disturbance estimation device 10 provided on the target ship T will be described. In the present embodiment, the target ship T is the own ship, but the disturbance estimation device 10 does not have to be provided in the ship, and if it is possible to communicate with the target ship T and other ships O1 to O4, For example, it may be a server installed on land. In addition, each of the target ship T and the other ships O1 to O4 detects the hull position from a satellite positioning system (GNSS) such as a speed log for detecting the hull speed, a gyrocompass for detecting the bow orientation of the hull, and GPS. It shall be provided with sensors including a GNSS sensor.

図2に示すように、外乱推測装置10は、ハードウェアとして、CPU(Central Processing Unit)11、メモリ12、記憶装置13、入出力I/F14、アンテナ15を備える。CPU11は、メモリ12と協働して後述する各種機能を実行する。メモリ12は、RAM(Random Access Memory)などの揮発性の主記憶装置である。記憶装置13は、CPU11及びメモリ12による各種機能の実行に伴う各種データを記憶する外部記憶装置である。入出力I/F14は、CPU11がアンテナ15と通信するためのインターフェイスである。アンテナ15は、後述する船舶情報を電波として送信または受信する。なお、外乱推測装置10は、入出力I/F14を介して、対応する船体に備えられたセンサ類による測定値を取得するものとする。 As shown in FIG. 2, the disturbance estimation device 10 includes a CPU (Central Processing Unit) 11, a memory 12, a storage device 13, an input / output I / F 14, and an antenna 15 as hardware. The CPU 11 cooperates with the memory 12 to execute various functions described later. The memory 12 is a volatile main storage device such as a RAM (Random Access Memory). The storage device 13 is an external storage device that stores various data associated with the execution of various functions by the CPU 11 and the memory 12. The input / output I / F 14 is an interface for the CPU 11 to communicate with the antenna 15. The antenna 15 transmits or receives ship information described later as radio waves. The disturbance estimation device 10 shall acquire the measured values by the sensors provided on the corresponding hull via the input / output I / F14.

また、外乱推測装置10は図3に示すように、機能として、外乱推定部101、送出部102、取得部103、時間判定部104、外乱推測部105、外乱補間部106、推測ベクトル算出107、推測航路算出部108、補償角算出部109を備え、外乱推測部105は、第1外乱推測部105aと第2外乱推測部105bとを含む。 Further, as shown in FIG. 3, the disturbance estimation device 10 has functions such as a disturbance estimation unit 101, a transmission unit 102, an acquisition unit 103, a time determination unit 104, a disturbance estimation unit 105, a disturbance interpolation unit 106, and an estimation vector calculation 107. The estimation route calculation unit 108 and the compensation angle calculation unit 109 are provided, and the disturbance estimation unit 105 includes a first disturbance estimation unit 105a and a second disturbance estimation unit 105b.

外乱推定部101は、対象船T、即ち自船に影響を及ぼす外乱を推定する。ここで、外乱は風及び潮流であり、外乱推定部101は、対象船Tに備えられたセンサ類から得られた測定値に基づいて、既存の手法により外乱ベクトルを推定する。送出部102は、外乱推定部101により推定された外乱ベクトルと、外乱ベクトルの算出時刻である時刻情報と、対象船Tの位置を示す位置情報とが対応付けられた船舶情報をアンテナ15に送信させる。取得部103は、他船O1〜O4より送信され、アンテナ15により受信された船舶情報と、航行ベクトルとしての対象船Tの航行方位及び対水速度とを取得して記憶装置13に記憶する。 The disturbance estimation unit 101 estimates the target ship T, that is, the disturbance that affects the own ship. Here, the disturbance is a wind and a tidal current, and the disturbance estimation unit 101 estimates the disturbance vector by an existing method based on the measured values obtained from the sensors provided on the target ship T. The sending unit 102 transmits to the antenna 15 ship information in which the disturbance vector estimated by the disturbance estimation unit 101, the time information which is the calculation time of the disturbance vector, and the position information indicating the position of the target ship T are associated with each other. Let me. The acquisition unit 103 acquires the ship information transmitted from the other ships O1 to O4 and received by the antenna 15 and the navigation direction and the water speed of the target ship T as a navigation vector and stores them in the storage device 13.

時間判定部104は、取得部103により取得された船舶情報に含まれる時刻情報に基づいて、これに対応付けられた外乱情報が有効であるか否かを判定する。第1外乱推測部105aは航路区間Sを含む狭域範囲内の他船から受信された船舶情報に基づいて、航路区間S上の一地点における外乱ベクトルを推測する。また、第2外乱推測部105bは、航路区間Sを含む広域範囲内の他船から受信された船舶情報に基づいて、航路区間S上の一地点における外乱ベクトルを推測する。 Based on the time information included in the ship information acquired by the acquisition unit 103, the time determination unit 104 determines whether or not the disturbance information associated with the time information is valid. The first disturbance estimation unit 105a estimates the disturbance vector at one point on the route section S based on the ship information received from another ship in the narrow area including the route section S. Further, the second disturbance estimation unit 105b estimates the disturbance vector at one point on the route section S based on the ship information received from another ship in the wide area including the route section S.

外乱補間部106は、航路区間S上において外乱ベクトルが推測された地点のうち、対象船Tの航行方向に連続する地点間にある補間地点における外乱を補間する。推測ベクトル算出部107は、外乱ベクトルと対象船Tの航行ベクトルとに基づいて、外乱が生じた場合の対象船Tのベクトルである推測ベクトルを算出する。 The disturbance interpolation unit 106 interpolates the disturbance at the interpolation point between the points where the disturbance vector is estimated on the route section S and which are continuous in the navigation direction of the target ship T. The guess vector calculation unit 107 calculates a guess vector, which is a vector of the target ship T when a disturbance occurs, based on the disturbance vector and the navigation vector of the target ship T.

推測航路算出部108は、推測ベクトルに基づいて、外乱が生じた場合の対象船Tの航路である推測航路を算出する。補償角算出部109は、対象船Tが推測航路を航行した場合に生じる航路区間S上にある到達地点からのずれを補償するための補償角を算出する。 The estimated route calculation unit 108 calculates the estimated route, which is the route of the target ship T when a disturbance occurs, based on the estimated vector. The compensation angle calculation unit 109 calculates the compensation angle for compensating for the deviation from the arrival point on the route section S that occurs when the target ship T navigates the estimated route.

(外乱推測装置の全体動作)
次に、外乱推測装置の動作について説明する。図4は、外乱推測装置の全体処理の動作を示すフローチャートである。なお、図4において、全ての他船から船舶情報がすでに送信されているものとし、また、全体処理は所定の周期毎に実行されているものとする。
(Overall operation of disturbance estimation device)
Next, the operation of the disturbance estimation device will be described. FIG. 4 is a flowchart showing the operation of the entire processing of the disturbance estimation device. In FIG. 4, it is assumed that ship information has already been transmitted from all other ships, and that the entire process is executed at predetermined intervals.

図4に示すように、まず、取得部103は、全ての他船O1〜O4からそれぞれ送信された船舶情報の全てと対象船Tの航行方位及び対水速度とを取得し(S101)、これらを記憶装置13に記憶する(S102)。次に、時間判定部104は、記憶装置13に記憶されたすべての船舶情報のうち、未選択の船舶情報があるか否かを判断する(S103)。 As shown in FIG. 4, first, the acquisition unit 103 acquires all the ship information transmitted from all the other ships O1 to O4, the navigation direction of the target ship T, and the water speed (S101). Is stored in the storage device 13 (S102). Next, the time determination unit 104 determines whether or not there is unselected ship information among all the ship information stored in the storage device 13 (S103).

未選択の船舶情報がある場合(S103,YES)、時間判定部104は、未選択の船舶情報を1つ選択し(S104)、選択した船舶情報に含まれている算出時刻が現時刻から所定時間以内であるか否かを判断する(S105)。 When there is unselected ship information (S103, YES), the time determination unit 104 selects one unselected ship information (S104), and the calculated time included in the selected ship information is predetermined from the current time. It is determined whether or not it is within the time (S105).

算出時刻が現時刻から所定時間以内である場合(S105,YES)、時間判定部104は、選択中の船舶情報を有効と判定する(S106)。記憶装置13には、リアルタイム、即ち全体処理の周期内に得られた船舶情報だけではなく、過去に得られた船舶情報も記憶されているため、時間判定部104は、このように現時刻から予め設定された所定時間内に算出された船舶情報のみを有効とする。次に、時間判定部104は、記憶装置13に記憶された全ての船舶情報を選択したか否かを判断する(S107)。 When the calculated time is within a predetermined time from the current time (S105, YES), the time determination unit 104 determines that the selected ship information is valid (S106). Since the storage device 13 stores not only the ship information obtained in real time, that is, within the cycle of the entire processing, but also the ship information obtained in the past, the time determination unit 104 thus starts from the current time. Only the ship information calculated within the preset predetermined time is valid. Next, the time determination unit 104 determines whether or not all the ship information stored in the storage device 13 has been selected (S107).

全ての船舶情報が選択された場合(S107,YES)、狭域外乱推測処理が実行され(S108)、広域外乱推測処理が実行され(S109)、推測航路算出処理が実行され(S110)、取得部103は、再度、全ての他船O1〜O4からそれぞれ送信された船舶情報を全て取得する(S101)。なお、狭域外乱推測処理、広域外乱推測処理、推測航路算出処理については後述する。 When all the ship information is selected (S107, YES), the narrow area disturbance estimation process is executed (S108), the wide area disturbance estimation process is executed (S109), the estimated route calculation process is executed (S110), and the acquisition is performed. Unit 103 again acquires all the ship information transmitted from all the other ships O1 to O4 (S101). The narrow area disturbance estimation process, wide area disturbance estimation process, and estimated route calculation process will be described later.

一方、全ての船舶情報が選択されていない場合(S107,NO)、時間判定部104は、再度、未選択の船舶情報を1つ選択する(S104)。 On the other hand, when all the ship information is not selected (S107, NO), the time determination unit 104 selects one unselected ship information again (S104).

また、ステップS105において、算出時刻が現時刻から所定時間以内ではない場合(S105,NO)、時間判定部104は、記憶装置13に記憶された全ての船舶情報を選択したか否かを判断する(S107)。 Further, in step S105, when the calculated time is not within a predetermined time from the current time (S105, NO), the time determination unit 104 determines whether or not all the ship information stored in the storage device 13 has been selected. (S107).

また、ステップS103において、未選択の船舶情報がない場合(S103,NO)、取得部103は、再度、全ての他船O1〜O4からそれぞれ送信された船舶情報を全て取得する(S101)。 Further, in step S103, when there is no unselected ship information (S103, NO), the acquisition unit 103 again acquires all the ship information transmitted from all the other ships O1 to O4 (S101).

(狭域外乱推測処理)
次に、狭域外乱推測処理について説明する。図5は、狭域外乱推測処理の動作を示すフローチャートである。図6は、狭域範囲内の他船を示す概略図である。
(Narrow region disturbance estimation processing)
Next, the narrow region disturbance estimation process will be described. FIG. 5 is a flowchart showing the operation of the narrow region disturbance estimation process. FIG. 6 is a schematic view showing another ship in a narrow area.

図5に示すように、まず、第1外乱推測部105aは、時間判定部104に有効と判定された船舶情報のうち、第1外乱推測部105aにより未選択の船舶情報があるか否かを判定する(S201)。 As shown in FIG. 5, first, the first disturbance estimation unit 105a determines whether or not there is ship information unselected by the first disturbance estimation unit 105a among the ship information determined to be valid by the time determination unit 104. Judgment (S201).

未選択の船舶情報がある場合(S201,YES)、第1外乱推測部105aは、有効と判定され、未選択である船舶情報を選択し(S202)、この船舶情報に含まれる位置情報により示される他船位置とこの他船位置から直近にある航路区間Sにおける一地点との距離を離間距離として算出する(S203)。更に第1外乱推測部105aは、この離間距離が予め設定された所定距離D1以内であり、且つ他船位置が航路区間Sにおいて対象船Tの前方にあるか否かを判断する(S204)。この判断によれば、図6に示すように、狭域範囲N内に船舶情報を送信した他船があるか否かが判定され、本実施形態においては、他船O3及び他船O4が狭域範囲N内にある他船に該当する。 When there is unselected ship information (S201, YES), the first disturbance estimation unit 105a determines that it is valid, selects unselected ship information (S202), and indicates by the position information included in this ship information. The distance between the position of the other ship and one point in the route section S closest to the position of the other ship is calculated as the separation distance (S203). Further, the first disturbance estimation unit 105a determines whether or not the separation distance is within the preset predetermined distance D1 and the position of the other ship is in front of the target ship T in the route section S (S204). According to this determination, as shown in FIG. 6, it is determined whether or not there is another ship that has transmitted the ship information within the narrow area range N, and in the present embodiment, the other ship O3 and the other ship O4 are narrow. It corresponds to another ship within the range N.

離間距離が所定距離D1以内であり、且つ他船位置が対象船Tの前方にある場合(S204,YES)、第1外乱推測部105aは、選択中の船舶情報に含まれる外乱ベクトルに基づいて、航路区間S上の一地点の外乱ベクトルを算出する(S205)。ここで、第1外乱推測部105aは、船舶情報に含まれる外乱ベクトルを、これに対応する位置情報に示される他船位置に対して航路区間S上において直近にある地点の外乱ベクトルとすることによって、航路区間S上の一地点の外乱ベクトルを算出する。これによって、図6に示すように、他船O3,O4それぞれの外乱ベクトルが航路区間S上の外乱ベクトルa2,a3となる。なお、図6における外乱ベクトルa1は、外乱推定部101により推定された対象船Tの外乱ベクトルである。次に、第1外乱推測部105aは、有効と判定された全ての船舶情報を選択したか否かを判断する(S206)。 When the separation distance is within the predetermined distance D1 and the position of the other ship is in front of the target ship T (S204, YES), the first disturbance estimation unit 105a is based on the disturbance vector included in the selected ship information. , Calculate the disturbance vector at one point on the route section S (S205). Here, the first disturbance estimation unit 105a sets the disturbance vector included in the ship information as the disturbance vector of the point closest to the other ship position shown in the corresponding position information on the route section S. Calculates the disturbance vector at one point on the route section S. As a result, as shown in FIG. 6, the disturbance vectors of the other vessels O3 and O4 become the disturbance vectors a2 and a3 on the route section S. The disturbance vector a1 in FIG. 6 is a disturbance vector of the target ship T estimated by the disturbance estimation unit 101. Next, the first disturbance estimation unit 105a determines whether or not all the ship information determined to be valid has been selected (S206).

全ての有効な船舶情報が選択された場合(S206,YES)、第1外乱推測部105aは狭域外乱推測処理を終了する。 When all the valid ship information is selected (S206, YES), the first disturbance estimation unit 105a ends the narrow region disturbance estimation process.

一方、全ての有効な船舶情報が選択されていない場合(S206,NO)、第1外乱推測部105aは、再度、有効と判定され、未選択である船舶情報を選択する(S202)。 On the other hand, when all the valid ship information is not selected (S206, NO), the first disturbance estimation unit 105a again determines that it is valid and selects the unselected ship information (S202).

また、ステップS204において、離間距離が所定距離D1以内でないか、または他船位置が対象船Tの前方にない場合(S204,NO)、第1外乱推測部105aは、有効と判定された全ての船舶情報を選択したか否かを判断する(S206)。 Further, in step S204, when the separation distance is not within the predetermined distance D1 or the position of the other ship is not in front of the target ship T (S204, NO), the first disturbance estimation unit 105a is all determined to be effective. It is determined whether or not the ship information is selected (S206).

また、ステップS201において、未選択の船舶情報がない場合(S201,NO)、第1外乱推測部105aは、狭域外乱推測処理を終了する。 Further, in step S201, when there is no unselected ship information (S201, NO), the first disturbance estimation unit 105a ends the narrow region disturbance estimation process.

(広域外乱推測処理)
次に、広域外乱推測処理の動作について説明する。図7は、広域外乱推測処理の動作を示すフローチャートである。図8は、広域範囲内の他船を示す概略図である。
(Wide area disturbance estimation processing)
Next, the operation of the wide area disturbance estimation process will be described. FIG. 7 is a flowchart showing the operation of the wide area disturbance estimation process. FIG. 8 is a schematic view showing another ship within a wide area.

図7に示すように、まず、第2外乱推測部105bは、対象船Tの前後の航路を含む広域範囲を設定し(S301)、その位置情報により示される位置が広域範囲内であり、且つ有効と判定された船舶情報を全て候補に追加する(S302)。図8に示すように、広域範囲Bは、航路区間Sに直交する対象船Tの左右方向の距離D2と、対象船Tの前方、後方の距離D3、D4により規定される範囲であり、距離D2は狭域範囲を規定する距離D1より大きく、また、距離D3は距離D4より大きくすると良い。本実施形態においては、他船O1〜O4が広域範囲内の他船に該当する。 As shown in FIG. 7, first, the second disturbance estimation unit 105b sets a wide range including the routes before and after the target ship T (S301), and the position indicated by the position information is within the wide range. All the ship information determined to be valid is added to the candidates (S302). As shown in FIG. 8, the wide area range B is a range defined by the distance D2 in the left-right direction of the target ship T orthogonal to the route section S and the distances D3 and D4 in front of and behind the target ship T. It is preferable that D2 is larger than the distance D1 that defines the narrow range, and the distance D3 is larger than the distance D4. In this embodiment, the other vessels O1 to O4 correspond to other vessels within a wide range.

次に、第2外乱推測部105bは、候補に2つ以上の船舶情報があるか否かを判断する(S303)。 Next, the second disturbance estimation unit 105b determines whether or not the candidate has two or more ship information (S303).

候補に2つ以上の船舶情報がある場合(S303,YES)、第2外乱推測部105bは、候補に始点として選択されていない船舶情報があるか否かを判定する(S304)。 When the candidate has two or more ship information (S303, YES), the second disturbance estimation unit 105b determines whether or not the candidate has ship information that has not been selected as the starting point (S304).

候補に始点として未選択の船舶情報がある場合(S304,YES)、第2外乱推測部105bは、候補における始点として未選択の船舶情報の1つを始点として選択し(S305)、候補に始点として、または選択中の始点に対する終点として選択されていない船舶情報があるか否かを判定する(S306)。 When the candidate has unselected ship information as the starting point (S304, YES), the second disturbance estimation unit 105b selects one of the unselected ship information as the starting point in the candidate as the starting point (S305), and the candidate has the starting point. It is determined whether or not there is ship information that has not been selected as, or as an end point with respect to the selected start point (S306).

候補に始点または選択中の始点に対する終点として未選択の船舶情報がある場合(S306,YES)、第2外乱推測部105bは、候補における始点または選択中の始点に対する終点として未選択の船舶情報の1つを終点として選択し(S307)、始点として選択した船舶情報の位置情報による位置と終点として選択した船舶情報の位置情報による位置とを接続する線分と、航路区間Sとが対象船Tの前方において交差するか否かを判断する(S308)。本実施形態においては、図8に示すように、他船O1と他船O4とを接続する線分、他船O2と他船O3とを接続する線分、及び他船O3と他船O4とを接続する線分が、対象船Tの前方において航路区間Sと交差する。 When the candidate has unselected vessel information as a start point or an end point for the selected start point (S306, YES), the second disturbance estimation unit 105b is the unselected vessel information as the end point for the start point in the candidate or the selected start point. One is selected as the end point (S307), and the line segment connecting the position based on the position information of the ship information selected as the start point and the position based on the position information of the ship information selected as the end point, and the route section S are the target ships T. It is determined whether or not to intersect in front of (S308). In the present embodiment, as shown in FIG. 8, the line segment connecting the other ship O1 and the other ship O4, the line segment connecting the other ship O2 and the other ship O3, and the other ship O3 and the other ship O4 The line segment connecting the two intersects the route section S in front of the target ship T.

線分と航路区間Sとが対象船Tの前方において交差する場合(S308,YES)、第2外乱推測部105bは、他船間を接続する線分と航路区間Sとの交点における外乱ベクトルを後述する算出方法により算出し(S309)、再度、候補に始点として選択されていない船舶情報があるか否かを判定する(S304)。 When the line segment and the route section S intersect in front of the target ship T (S308, YES), the second disturbance estimation unit 105b sets the disturbance vector at the intersection of the line segment connecting between other ships and the route section S. It is calculated by the calculation method described later (S309), and again, it is determined whether or not there is ship information that has not been selected as the starting point in the candidates (S304).

一方、線分と航路区間Sとが対象船Tの前方において交差しない場合(S308,NO)、第2外乱推測部105bは、再度、候補に始点として選択されていない船舶情報があるか否かを判定する(S304)。 On the other hand, when the line segment and the route section S do not intersect in front of the target ship T (S308, NO), the second disturbance estimation unit 105b again determines whether or not the candidate has ship information that has not been selected as the starting point. Is determined (S304).

また、ステップS306において、候補に始点としてまたは選択中の始点に対する終点として未選択の船舶情報がない場合(S306,NO)、第2外乱推測部105bは、再度、候補に始点として選択されていない船舶情報があるか否かを判定する(S304)。 Further, in step S306, when there is no unselected ship information as the start point in the candidate or as the end point with respect to the selected start point (S306, NO), the second disturbance estimation unit 105b is not selected as the start point in the candidate again. It is determined whether or not there is ship information (S304).

また、ステップS304において、候補に始点として未選択の船舶情報がない場合(S304,NO)、第2外乱推測部105bは、広域外乱推測処理を終了する。 Further, in step S304, when there is no unselected ship information as a starting point in the candidate (S304, NO), the second disturbance estimation unit 105b ends the wide area disturbance estimation process.

また、ステップS302において、候補に2つ以上の船舶情報がない場合(S303,NO)、第2外乱推測部105bは、広域外乱推測処理を終了する。 Further, in step S302, when the candidate does not have two or more ship information (S303, NO), the second disturbance estimation unit 105b ends the wide area disturbance estimation process.

(交点における外乱ベクトルの算出方法)
次に、上述した広域外乱推測処理における外乱ベクトルの算出方法について説明する。図9は、交点における外乱ベクトルの算出方法を示す概略図である。
(Calculation method of disturbance vector at intersection)
Next, a method of calculating the disturbance vector in the above-mentioned wide area disturbance estimation process will be described. FIG. 9 is a schematic diagram showing a method of calculating the disturbance vector at the intersection.

2つの船舶情報における位置情報に示される位置(他船位置)間を接続する線分と航路区間Sとの交点における外乱ベクトルは、2つの船舶情報のそれぞれにおける外乱ベクトルと、一方の他船位置から交点までの距離と他方の他船位置から交点までの距離との比とに基づいて算出される。具体的には、図9に示すように、線分の一方の端点としての他船位置をAとし、他方の端点としての他船位置をBとし、交点をPとし、更にAとPとの距離をn、BとPとの距離をmとした場合、以下の式により交点における外乱ベクトルが算出される。 The disturbance vector at the intersection of the line segment connecting between the positions (other ship positions) shown in the position information in the two ship information and the route section S is the disturbance vector in each of the two ship information and the position of one other ship. It is calculated based on the ratio of the distance from to the intersection to the distance from the other ship's position to the intersection. Specifically, as shown in FIG. 9, the position of the other ship as one end point of the line segment is A, the position of the other ship as the other end point is B, the intersection is P, and A and P Assuming that the distance is n and the distance between B and P is m, the disturbance vector at the intersection is calculated by the following equation.

Figure 0006864483
Figure 0006864483

(推測航路算出処理)
次に、推測航路算出処理について説明する。図10は、航路推測処理の動作を示すフローチャートである。図11は、外乱ベクトルの補間方法を示す概略図である。図12は、推測ベクトルの算出方法を示す概略図である。図13は、航路の算出方法を示す概略図である。図14は、対象船の補償角を示す概略図である。図15は、補償角により変針された対象船を示す概略図である。
(Estimated route calculation process)
Next, the estimated route calculation process will be described. FIG. 10 is a flowchart showing the operation of the route estimation process. FIG. 11 is a schematic view showing a method of interpolating the disturbance vector. FIG. 12 is a schematic diagram showing a method of calculating the guess vector. FIG. 13 is a schematic view showing a method of calculating a route. FIG. 14 is a schematic view showing the compensation angle of the target ship. FIG. 15 is a schematic view showing a target ship whose needle has been changed by the compensation angle.

図10に示すように、まず、外乱補間部 106は、すでに算出された複数の外乱ベクトル間を、図11に示すように航路区間S上で線形補完することにより、航路区間S上の外乱ベクトルを更に詳細に推測する(S401)。図11において、a1は外乱推定部101により推定された対象船Tの外乱ベクトルを示し、a2及びa3は狭域外乱推測処理により推測された外乱ベクトルを示し、b1〜b3は広域外乱推測処理により推測された外乱ベクトルを示し、c1〜c5は線形補完された外乱ベクトルを示す。 As shown in FIG. 10, first, the disturbance interpolation unit 106 linearly complements the already calculated plurality of disturbance vectors on the route section S as shown in FIG. 11, thereby causing the disturbance vector on the route section S. Is inferred in more detail (S401). In FIG. 11, a1 shows the disturbance vector of the target ship T estimated by the disturbance estimation unit 101, a2 and a3 show the disturbance vector estimated by the narrow region disturbance estimation process, and b1 to b3 show the disturbance vector estimated by the wide area disturbance estimation process. The inferred disturbance vector is shown, and c1 to c5 indicate the linearly complemented disturbance vector.

次に、推測ベクトル算出部107は、航路区間S上の外乱ベクトルそれぞれについて、この外乱ベクトルと取得部103により取得された対象船Tの航行ベクトルとに基づいて、図12に示すように推測ベクトルE1〜Enを順次算出する(S402)。この推測ベクトルは、個々の外乱ベクトルと対象船Tの航行ベクトルとを合成することによって算出される。 Next, the guess vector calculation unit 107 makes a guess vector for each of the disturbance vectors on the route section S, based on the disturbance vector and the navigation vector of the target ship T acquired by the acquisition unit 103, as shown in FIG. E1 ~ En are sequentially calculated (S402). This guess vector is calculated by synthesizing the individual disturbance vectors and the navigation vector of the target ship T.

推測ベクトルの算出後、推測航路算出部108は、複数の推測ベクトルを、図13に示すように順次累積することによって、推測航路を算出する(S403)。具体的には、対象船Tの現在地点である位置P0を起点として推測ベクトルを順次累積することで、通過地点をP1〜Pn−1、到達地点をPnとし、外乱の影響を加味した対象船Tの航路である推測航路を算出する。 After calculating the guess vector, the guess route calculation unit 108 calculates the guess route by sequentially accumulating a plurality of guess vectors as shown in FIG. 13 (S403). Specifically, by sequentially accumulating the estimation vectors starting from the position P0, which is the current position of the target ship T, the passing points are set to P1 to Pn-1, the reaching points are set to Pn, and the target ship is affected by the disturbance. The estimated route, which is the route of T, is calculated.

次に、補償角算出部109は、推測航路に基づいて、図14に示すように、航路区間Sと、位置P0と位置Pnとを接続する線分とが成す角度を補償角度Cとして算出する(S404)。この補償角度Cによれば、図15に示すように、対象船Tを、現行の航行方向に対して−Cだけ変針させることによって、外乱の影響による到達地点のずれを解消させつつ、対象船Tに外乱を考慮した航路を航行させることができる。 Next, the compensation angle calculation unit 109 calculates the angle formed by the route section S and the line segment connecting the position P0 and the position Pn as the compensation angle C, as shown in FIG. 14, based on the estimated route. (S404). According to this compensation angle C, as shown in FIG. 15, by changing the needle of the target ship T by -C with respect to the current navigation direction, the target ship while eliminating the deviation of the arrival point due to the influence of the disturbance. It is possible to make T navigate a route in consideration of disturbance.

以上説明したように、他船において算出された外乱ベクトルに基づいて航路区間上の外乱ベクトルを推測することによって、従来と比較してより精度高く外乱を推測することができ、延いては、対象船の航路を航行時間及び消費エネルギーの観点からより最適なものとすることができる。 As explained above, by estimating the disturbance vector on the route section based on the disturbance vector calculated by the other ship, it is possible to estimate the disturbance with higher accuracy than before, and by extension, the target. The route of the ship can be optimized from the viewpoint of navigation time and energy consumption.

本発明の実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら新規な実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれるとともに、特許請求の範囲に記載された発明とその均等の範囲に含まれる。 The embodiments of the present invention are presented as examples and are not intended to limit the scope of the invention. These novel embodiments can be implemented in various other embodiments, and various omissions, replacements, and changes can be made without departing from the gist of the invention. These embodiments and modifications thereof are included in the scope and gist of the invention, and are also included in the scope of the invention described in the claims and the equivalent scope thereof.

10 外乱推測装置
103 取得部
104 時間判定部
105 外乱推定部
106 外乱補間部
108 推測航路算出部
109 補償角算出部
10 Disturbance estimation device 103 Acquisition unit 104 Time determination unit 105 Disturbance estimation unit 106 Disturbance interpolation unit 108 Guessing route calculation unit 109 Compensation angle calculation unit

Claims (7)

予め設定された設定航路を航行する対象船に対する外乱ベクトルを推測する外乱推測装置であって、
前記対象船とは異なる少なくとも1つ以上の他船から、該他船において算出された該他船に対する外乱ベクトルと該他船の位置とを含み、該他船から送信された船舶情報を取得する取得部と、
前記船舶情報に基づいて、前記設定航路上において前記対象船の前方に位置する所定地点の外乱ベクトルを推測する外乱推測部と
を備える外乱推測装置。
It is a disturbance estimation device that infers the disturbance vector for the target ship navigating the preset set route.
Acquires different from at least one or more other vessels, viewed including the position of the disturbance vector and said other vessels for the another ship calculated in said other ship, ship information transmitted from the other ship and the target ship Acquisition department and
A disturbance estimation device including a disturbance estimation unit that estimates a disturbance vector at a predetermined point located in front of the target ship on the set route based on the ship information.
前記外乱推測部は、前記船舶情報のうち、該船舶情報に含まれる位置が前記設定航路に対して前記対象船の左右方向それぞれに第1距離を有する第1範囲内である船舶情報を第1船舶情報として選択し、該第1船舶情報に含まれる外乱ベクトルを該第1船舶情報に含まれる位置から最も近接する前記設定航路上の一地点の外乱ベクトルとすることを特徴とする請求項1に記載の外乱推測装置。 The disturbance estimation unit first obtains ship information in the ship information in which the position included in the ship information is within the first range having a first distance in each of the left and right directions of the target ship with respect to the set route. Claim 1 is selected as ship information, and the disturbance vector included in the first ship information is set as the disturbance vector at one point on the set route closest to the position included in the first ship information. Disturbance estimation device described in. 前記取得部は、少なくとも2つ以上の他船から前記船舶情報を取得し、
前記外乱推測部は、前記船舶情報のうち、該船舶情報に含まれる位置が前記設定航路に対して前記対象船の左右方向それぞれに前記第1距離より大きい第2距離を有する第2範囲内である船舶情報を第2船舶情報として複数選択し、該複数の第2船舶情報それぞれに含まれる位置同士を接続した線分と前記設定航路との交点における外乱ベクトルを、前記線分の端点である位置が含まれる複数の第2船舶情報それぞれに含まれる外乱ベクトルに基づいて算出することを特徴とする、請求項2に記載の外乱推測装置。
The acquisition unit acquires the ship information from at least two or more other ships, and obtains the ship information.
The disturbance estimation unit is within a second range of the ship information in which the position included in the ship information has a second distance larger than the first distance in each of the left and right directions of the target ship with respect to the set route. A plurality of certain ship information is selected as the second ship information, and the disturbance vector at the intersection of the line segment connecting the positions included in each of the plurality of second ship information and the set route is the end point of the line segment. position and calculates, based on the disturbance vectors contained in each of the plurality of second vessels that contain disturbance guess device according toMotomeko 2.
前記船舶情報は、前他船に対する外乱ベクトルの算出時刻を示す時刻情報を更に含み、
前記時刻情報に基づいて、該時刻情報を含む船舶情報が有効であるか否かを判定する時間判定部を更に備え、
前記外乱推測部は、有効と判定された船舶情報に基づいて前記設定航路上における所定地点の外乱ベクトルを推測することを特徴とする、請求項1〜請求項3のいずれか一項に記載の外乱推測装置。
The ship information further includes time information indicating the calculation time of the disturbance vector with respect to the other ship in front.
Further, a time determination unit for determining whether or not the ship information including the time information is valid based on the time information is provided.
The disturbance estimation unit according to any one of claims 1 to 3, wherein the disturbance estimation unit estimates a disturbance vector at a predetermined point on the set route based on the ship information determined to be valid. Disturbance guesser.
前記外乱推測部は、前記設定航路上における複数の推測地点における外乱ベクトルを推測し、
前記設定航路上における前記複数の推測地点間にある地点の外乱ベクトルを、該複数の推測地点において外乱ベクトルに基づく線形補完により算出する外乱補間部を更に備えることを特徴とする、請求項1〜請求項4のいずれか一項に記載の外乱推測装置。
The disturbance estimation unit estimates disturbance vectors at a plurality of estimation points on the set route, and estimates the disturbance vectors.
Claims 1 to 1, further comprising a disturbance interpolation unit that calculates a disturbance vector at a point between the plurality of estimated points on the set route by linear interpolation based on the disturbance vector at the plurality of estimated points. The disturbance estimation device according to any one of claims 4.
前記外乱推測部により推測された外乱ベクトルと、前記外乱補間部により算出された外乱ベクトルとに基づいて、前記設定航路上の外乱を考慮した航路である推測航路を算出する推測航路算出部を更に備えることを特徴とする請求項5に記載の外乱推測装置。 Based on the disturbance vector estimated by the disturbance estimation unit and the disturbance vector calculated by the disturbance interpolation unit, the estimated route calculation unit that calculates the estimated route, which is the route considering the disturbance on the set route, is further added. The disturbance estimation device according to claim 5, further comprising. 前記推測航路における前記対象船の起点と到達地点とを接続する線分と前記設定航路とが成す補償角を算出する補償角算出部を更に備えることを特徴とする請求項6に記載の外乱推測装置。 The disturbance estimation according to claim 6, further comprising a compensation angle calculation unit for calculating the compensation angle formed by the line segment connecting the starting point and the arrival point of the target ship in the estimated route and the set route. apparatus.
JP2017007807A 2017-01-19 2017-01-19 Disturbance guesser Active JP6864483B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2017007807A JP6864483B2 (en) 2017-01-19 2017-01-19 Disturbance guesser

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2017007807A JP6864483B2 (en) 2017-01-19 2017-01-19 Disturbance guesser

Publications (2)

Publication Number Publication Date
JP2018114897A JP2018114897A (en) 2018-07-26
JP6864483B2 true JP6864483B2 (en) 2021-04-28

Family

ID=62983758

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2017007807A Active JP6864483B2 (en) 2017-01-19 2017-01-19 Disturbance guesser

Country Status (1)

Country Link
JP (1) JP6864483B2 (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
TWI714040B (en) * 2019-03-27 2020-12-21 財團法人船舶暨海洋產業研發中心 A vessel navigation system and navigation method thereof

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0820391A (en) * 1994-07-05 1996-01-23 Mitsubishi Heavy Ind Ltd Steering system for ship
JP3760159B2 (en) * 2003-04-15 2006-03-29 欣三 井上 Navigation safety evaluation method, navigation safety evaluation system, and navigation safety evaluation program
JP5956233B2 (en) * 2012-04-25 2016-07-27 Jmuディフェンスシステムズ株式会社 Route holding control device and ship
JP5442071B2 (en) * 2012-07-04 2014-03-12 ジャパンマリンユナイテッド株式会社 Ship maneuvering control apparatus, automatic ship maneuvering control system, ship maneuvering control method, and program
JP6038087B2 (en) * 2014-09-22 2016-12-07 三菱重工業株式会社 Ship operation plan optimization system and ship operation plan optimization method

Also Published As

Publication number Publication date
JP2018114897A (en) 2018-07-26

Similar Documents

Publication Publication Date Title
US9541403B2 (en) Method and device for displaying route
US20150330803A1 (en) Information display device and method
CN108362281A (en) A kind of Long baselines underwater submarine matching navigation method and system
US11772758B2 (en) Apparatus, method, and recording medium for autonomous ship navigation
JP2008145312A (en) Optimum route search method
JPWO2016017358A1 (en) Ship characteristic estimation device and automatic steering device
GB2474715A (en) Aiding navigation of a marine vessel in a tidal region
JP6567665B2 (en) A method for estimating each drift (floating) vector at all points in a ship's route
JP6864483B2 (en) Disturbance guesser
CN110543179A (en) water surface target avoiding method based on three-dimensional sparse array sound source orientation recognition
CN113176588B (en) Distance measuring method, device and system and readable storage medium
KR20180076936A (en) Optimum Ocean Route guide system and method using low cost estimation device
JP2013079814A (en) Fish shoal prediction position display system, underwater detection device, storage device, fish shoal prediction device, fish shoal prediction position display program, and fish shoal prediction position display method
Jin et al. A compensation algorithm with motion constraint in DVL/SINS tightly coupled positioning
KR20170061534A (en) Method for ship sailing considering sea weather data
JP6014884B1 (en) Sensor laying apparatus, sound source object estimation system and method
CN114096887A (en) Method for generating a three-dimensional environment model using GNSS measurements
US20210018319A1 (en) Search support apparatus, search support method, and computer-readable recording medium
EP3851367A1 (en) Submarine machine system and work method
WO2023037500A1 (en) Information processing device, determining method, program, and storage medium
JP6874992B2 (en) Depth measuring device and nautical chart image display device
JP7362343B2 (en) Underwater work system
JP6300344B2 (en) Ocean current estimation method and ocean current estimation system
WO2022091677A1 (en) Seacraft monitoring system, seacraft monitoring method, information processing device, and program
JP2023040680A (en) Information processing device, control method, program, and storage medium

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20191001

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20200728

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20200811

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20201001

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20210315

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20210402

R150 Certificate of patent or registration of utility model

Ref document number: 6864483

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250