JP2021142873A - Learning model generation method, probability determination method and determination device - Google Patents

Learning model generation method, probability determination method and determination device Download PDF

Info

Publication number
JP2021142873A
JP2021142873A JP2020042630A JP2020042630A JP2021142873A JP 2021142873 A JP2021142873 A JP 2021142873A JP 2020042630 A JP2020042630 A JP 2020042630A JP 2020042630 A JP2020042630 A JP 2020042630A JP 2021142873 A JP2021142873 A JP 2021142873A
Authority
JP
Japan
Prior art keywords
station
data
error
driving
learning model
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.)
Granted
Application number
JP2020042630A
Other languages
Japanese (ja)
Other versions
JP7328921B2 (en
Inventor
大輔 鈴木
Daisuke Suzuki
大輔 鈴木
綾子 鈴木
Ayako Suzuki
綾子 鈴木
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.)
Railway Technical Research Institute
Original Assignee
Railway Technical Research Institute
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 Railway Technical Research Institute filed Critical Railway Technical Research Institute
Priority to JP2020042630A priority Critical patent/JP7328921B2/en
Publication of JP2021142873A publication Critical patent/JP2021142873A/en
Application granted granted Critical
Publication of JP7328921B2 publication Critical patent/JP7328921B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Train Traffic Observation, Control, And Security (AREA)

Abstract

To provide a technique for predicting error event occurrence possibility at individual stations by use of operation information recorded and accumulated.SOLUTION: When station-specific operation data of a given station is inputted with data including station identification information 311A, which represents either a first class station at which error operation exceeding a fixed position to stop has not occurred during a prescribed period or a second class station at which the error operation has occurred, and station-specific operation data 312, which is statistic data of a plurality of information concerning operation performed when stopping at the station during the prescribed period, as teacher data 310A, a machine learning model 10A, which outputs a prediction value representing the possibility that the given station relevant to the first class station or the second class station, is created when the station-specific data of the given station has been inputted.SELECTED DRAWING: Figure 1

Description

本発明は、停止すべき駅でエラー運転が発生する蓋然性を判定する蓋然性判定方法等に関する。 The present invention relates to a probability determination method for determining the probability that an error operation will occur at a station to be stopped.

鉄道においては、運転情報記録装置(運転状況記録装置ともいわれる)の設置が省令(鉄道に関する技術上の基準を定める省令)で義務付けられている。運転情報記録装置は、鉄道車両に搭載されて、列車の位置や速度、運転士の操作状況や主要機器の動作状況といった列車の運転状況に関する様々な情報を記録する装置である(例えば、特許文献1参照)。 In railways, the installation of a driving information recording device (also called an operating status recording device) is obligatory by a ministerial ordinance (a ministerial ordinance that sets technical standards for railways). The operation information recording device is a device mounted on a railway vehicle and records various information related to the operation status of the train such as the position and speed of the train, the operation status of the driver, and the operation status of the main equipment (for example, patent documents). 1).

特開2008−221902号公報Japanese Unexamined Patent Publication No. 2008-221902

運転情報記録装置により記録された運転情報は、主に、事故等が発生した場合にその発生要因の解析に利用されている。また、記録・蓄積された運転情報をもとに、ヒューマンエラーを含む各種のエラー事象の発生要因を体系的に分析することで、同種のエラー事象の再発防止に役立てる、といった活用もなされている。現状では過去のエラー事象の発生要因の分析が主であるが、今後の課題として、将来のエラー事象の発生可能性の予測といった、記録・蓄積された運転情報の新たな活用の仕方が求められている。例えば、個々の駅における運転士の運転操作に起因するエラー事象の発生可能性の予測、が挙げられる。 The driving information recorded by the driving information recording device is mainly used for analysis of the cause of an accident or the like. In addition, by systematically analyzing the causes of various error events including human error based on the recorded and accumulated driving information, it is also used to help prevent the recurrence of similar error events. .. Currently, the main focus is on analyzing the causes of past error events, but as a future task, new ways of utilizing recorded and accumulated driving information, such as predicting the possibility of future error events, are required. ing. For example, prediction of the possibility of error events caused by the driver's driving operation at each station.

本発明は、上記事情に鑑みてなされたものであり、その目的とするところは、記録・蓄積された運転情報を用いて、個々の駅における運転士の運転操作に起因するエラー事象の発生可能性を予測する技術を提供すること、である。 The present invention has been made in view of the above circumstances, and an object of the present invention is to be able to generate an error event caused by a driver's driving operation at an individual station by using recorded and accumulated driving information. To provide a technique for predicting sex.

上記課題を解決するための第1の発明は、
停止すべき定位置を超過したエラー運転が所定期間の間に発生していない第1分類駅と、発生した第2分類駅とについて、当該所定期間の間に当該駅に停止する際に行われた複数の運転情報の統計データである駅別運転データを取得することと(例えば、図5,図9の教師データ取得部202)、
前記駅別運転データと、当該駅別運転データが前記第1分類駅の情報か前記第2分類駅の情報かを示す駅識別情報とを教師データとして、所与の駅の駅別運転データが入力されたときに、当該所与の駅が前記第1分類駅又は前記第2分類駅に該当する可能性を示す予測値を出力する機械学習モデルを生成することと(例えば、図5,図9の予測モデル生成部204)、
を含む学習モデル生成方法である。
The first invention for solving the above problems is
It is performed when the first-class station and the second-class station where error operation exceeding the fixed position to be stopped did not occur during the predetermined period and stopped at the station during the predetermined period. Acquiring station-specific driving data, which is statistical data of a plurality of driving information (for example, teacher data acquisition unit 202 in FIGS. 5 and 9),
The station-specific operation data of a given station is obtained by using the station-specific operation data and the station identification information indicating whether the station-specific operation data is the information of the first category station or the information of the second category station as teacher data. To generate a machine learning model that outputs a predicted value indicating the possibility that the given station corresponds to the first classification station or the second classification station when input (for example, FIGS. 5 and 5). 9 Prediction model generation unit 204),
It is a learning model generation method including.

第1の発明によれば、記録・蓄積された運転情報を用いて、個々の駅におけるエラー事象の発生可能性を予測することが可能となる。つまり、エラー事象として停止すべき駅の定位置を超過したエラー運転が発生していない第1分類駅と、発生した第2分類駅とのそれぞれの駅別運転データを教師データとして、所与の駅の駅別運転データを入力すると、当該駅が第1分類駅又は第2分類駅に該当する可能性を示す予測値を出力する機械学習モデルを生成することができる。機械学習モデルが出力する予測値は、所与の駅に停止する際に行われた運転の仕方が、エラー運転が発生した駅と発生していない駅とのどちらの駅において行われた運転の仕方にどの程度似ているかを示している。従って、機械学習モデルが出力する予測値から、所与の駅において、エラー運転が発生する可能性を予測することができる。 According to the first invention, it is possible to predict the possibility of error events occurring at individual stations by using the recorded and accumulated driving information. That is, given as teacher data, the operation data for each station of the first-class station where the error operation exceeding the fixed position of the station to be stopped as an error event does not occur and the second-class station where the error operation has occurred is used as the teacher data. By inputting the station-specific operation data of a station, it is possible to generate a machine learning model that outputs a predicted value indicating the possibility that the station corresponds to a first-class station or a second-class station. The predicted value output by the machine learning model is that the driving method performed when stopping at a given station is the driving performed at either the station where the error driving occurred or the station where the error driving did not occur. It shows how similar it is to the way it is done. Therefore, it is possible to predict the possibility of error operation occurring at a given station from the predicted value output by the machine learning model.

第2の発明は、第1の発明において、
前記駅別運転データには、当該駅の定位置に停止するまでの速度パターンの種類数が含まれる、
学習モデル生成方法である。
The second invention is the first invention.
The station-specific operation data includes the number of types of speed patterns until the station stops at a fixed position.
This is a learning model generation method.

第2の発明によれば、駅の定位置に停止するまでの速度パターンの種類数をエラー運転が発生するか否かの要因に含めて予測値を算出する機械学習モデルを生成することができる。 According to the second invention, it is possible to generate a machine learning model that calculates a predicted value by including the number of types of speed patterns until the vehicle stops at a fixed position of a station as a factor of whether or not error driving occurs. ..

第3の発明は、第1又は第2の発明において、
前記駅別運転データには、常用最大ブレーキの使用割合が含まれ、
前記機械学習モデルは、ブレーキ遅れを主因とする前記エラー運転に係る前記予測値を出力するためのモデルである(例えば、第1実施形態における予測モデル10A)、
学習モデル生成方法である。
The third invention is the first or second invention.
The operation data for each station includes the usage rate of the maximum regular brake.
The machine learning model is a model for outputting the predicted value related to the error driving mainly due to the brake delay (for example, the predicted model 10A in the first embodiment).
This is a learning model generation method.

第3の発明によれば、常用最大ブレーキの使用割合をエラー運転が発生するか否かの要因に含めて予測値を算出する機械学習モデルを生成することができる。 According to the third invention, it is possible to generate a machine learning model for calculating a predicted value by including the usage ratio of the maximum normal brake as a factor of whether or not error driving occurs.

第4の発明は、第3の発明において、
前記駅別運転データには、前記定位置から所定距離手前の位置での走行速度のバラツキに関する情報が含まれる、
学習モデル生成方法である。
The fourth invention is the third invention.
The station-specific driving data includes information on variations in traveling speed at a position a predetermined distance before the fixed position.
This is a learning model generation method.

第4の発明によれば、定位置から所定距離手前の位置での走行速度のバラツキに関する情報をエラー運転が発生するか否かの要因に含めて予測値を算出する機械学習モデルを生成することができる。“速度のバラツキ”とは、例えば、速度の標準偏差である。 According to the fourth invention, a machine learning model for calculating a predicted value is generated by including information on a variation in traveling speed at a position a predetermined distance before a fixed position as a factor of whether or not error driving occurs. Can be done. The "velocity variation" is, for example, the standard deviation of velocity.

第5の発明は、第1又は第2の発明において、
前記取得することは、当該駅の定位置が編成車両数に応じて異なるか否かに関する定位置可変情報を更に取得し、
前記生成することは、前記教師データに前記定位置可変情報を更に含めて前記機械学習モデルを生成する(例えば、第2実施形態における予測モデル10B)、
学習モデル生成方法である。
The fifth invention is the first or second invention.
The acquisition further acquires the fixed position variable information regarding whether or not the fixed position of the station differs depending on the number of trains.
The generation generates the machine learning model by further including the fixed position variable information in the teacher data (for example, the prediction model 10B in the second embodiment).
This is a learning model generation method.

第5の発明によれば、駅の定位置が編成車両数に応じて異なるか否かに関する定位置可変情報をエラー運転が発生するか否かの要因に含めて予測値を算出する機械学習モデルを生成することができる。“定位置可変情報”には、例えば、編成車両数が「n」の列車(n両編成の列車)の停止位置目標を示すn両標の有無や、停止位置目標の数等がある。 According to the fifth invention, a machine learning model that calculates a predicted value by including fixed position variable information regarding whether or not the fixed position of a station differs according to the number of trains in the factor of whether or not error driving occurs. Can be generated. The "fixed position variable information" includes, for example, the presence / absence of an n-car mark indicating a stop position target of a train having an "n" number of trains (a train of an n-car train), the number of stop position targets, and the like.

第6の発明は、
第1〜第5の何れかの発明の学習モデル生成方法によって生成された機械学習モデルに対して、所与の駅に係る駅別運転データを入力することで得られる前記予測値に基づいて、当該駅で前記エラー運転が発生する蓋然性を判定する(例えば、図5,図9の蓋然性判定部210)、
蓋然性判定方法である。
The sixth invention is
Based on the predicted value obtained by inputting station-specific operation data for a given station to the machine learning model generated by the learning model generation method of any one of the first to fifth inventions. Determining the probability that the error operation will occur at the station (for example, the probability determination unit 210 in FIGS. 5 and 9),
This is a probability determination method.

他の発明として、
停止すべき定位置を超過したエラー運転が所定期間の間に発生していない第1分類駅と、発生した第2分類駅とについて、当該所定期間の間に当該駅に停止する際に行われた複数の運転情報の統計データである駅別運転データと、当該駅別運転データが前記第1分類駅の情報か前記第2分類駅の情報かを示す駅識別情報とを教師データとして、所与の駅に係る駅別運転データが入力されたときに、当該所与の駅が前記第1分類駅又は前記第2分類駅に該当する可能性を示す予測値を出力するように学習させた機械学習モデルに、所与の駅に係る駅別運転データを入力して、当該所与の駅の前記予測値を得る予測処理を実行する予測処理実行手段(例えば、図5,図9の予測部208)と、
前記予測処理で得られた前記予測値に基づいて、前記入力した所与の駅で前記エラー運転が発生する蓋然性を判定する蓋然性判定手段(例えば、図5,図9の蓋然性判定部210)と、
を備える判定装置を構成してもよい。
As another invention
It is performed when the first category station where the error operation exceeding the fixed position to be stopped does not occur during the predetermined period and the second category station where it occurs stops at the station during the predetermined period. Station-specific driving data, which is statistical data of a plurality of driving information, and station identification information indicating whether the station-specific driving data is the information of the first classification station or the information of the second classification station are used as teacher data. When the station-specific operation data for a given station is input, the trained trainer is to output a predicted value indicating the possibility that the given station corresponds to the first classification station or the second classification station. Prediction processing execution means (for example, prediction of FIGS. 5 and 9) for inputting station-specific operation data for a given station into a machine learning model and executing a prediction process for obtaining the predicted value of the given station. Part 208) and
With the probability determination means (for example, the probability determination unit 210 in FIGS. 5 and 9) for determining the probability that the error driving will occur at the input given station based on the prediction value obtained in the prediction process. ,
A determination device including the above may be configured.

第6の発明等によれば、機械学習モデルに対して所与の駅の駅別運転データを入力することで得られる予測値に基づいて、当該所与の駅においてエラー運転が発生する蓋然性を判定することができる。つまり、機械学習モデルが出力する予測値は、所与の駅に停止する際に行われる運転の仕方が、エラー運転が発生した駅と発生していない駅とのどちらの駅において行われた運転の仕方にどの程度似ているかを示している。従って、例えば、所与の駅においては未だエラー運転が発生していないが、機械学習モデルが出力した予測値がエラー運転が発生した第2分類駅に該当する可能性が高いことを示す場合には、近い将来にエラー運転が発生する可能性が高い、と判定することができる。 According to the sixth invention and the like, based on the predicted value obtained by inputting the station-specific operation data of a given station into the machine learning model, the probability that an error operation will occur at the given station is determined. Can be determined. In other words, the predicted value output by the machine learning model is that the driving method performed when stopping at a given station is the driving performed at either the station where the error driving occurred or the station where the error driving did not occur. It shows how similar it is to. Therefore, for example, when an error operation has not yet occurred at a given station, but the predicted value output by the machine learning model indicates that it is highly likely that the station corresponds to the second category station where the error operation has occurred. Can determine that there is a high possibility that error operation will occur in the near future.

第1実施形態における予測モデルの生成の説明図。The explanatory view of the generation of the prediction model in 1st Embodiment. 駅別運転データの生成の説明図。Explanatory drawing of generation of operation data for each station. 予測モデルを用いたエラー運転が発生する蓋然性の判定の説明図。Explanatory drawing for determination of probability that error driving occurs using a prediction model. 判定結果の一例。An example of the judgment result. 判定装置の機能構成図。Functional configuration diagram of the judgment device. 第2実施形態における予測モデルの生成の説明図。The explanatory view of the generation of the prediction model in 2nd Embodiment. 予測モデルを用いたエラー運転が発生する蓋然性の判定の説明図。Explanatory drawing for determination of probability that error driving occurs using a prediction model. 判定結果の一例。An example of the judgment result. 判定装置の機能構成図。Functional configuration diagram of the judgment device.

以下、図面を参照して本発明の好適な実施形態について説明する。なお、以下に説明する実施形態によって本発明が限定されるものではなく、本発明を適用可能な形態が以下の実施形態に限定されるものでもない。また、図面の記載において、同一要素には同一符号を付す。 Hereinafter, preferred embodiments of the present invention will be described with reference to the drawings. It should be noted that the embodiments described below do not limit the present invention, and the embodiments to which the present invention can be applied are not limited to the following embodiments. Further, in the description of the drawings, the same elements are designated by the same reference numerals.

[概要]
本実施形態は、鉄道車両に搭載される運転情報記録装置に記録された運転情報を用いて、駅の停止すべき定位置を超過したエラー運転が発生する蓋然性を判定するものである。定位置を超過したエラー運転には、定位置の過走(オーバーラン)と、停止せずに駅を通過したこととを含む。運転情報は、時刻に対応付けて、列車の位置や速度、ノッチ数を含むブレーキの操作状況や動作状況等を時系列に記録したデータである。運転情報は、運転情報記録装置に記録されたデータであることとしたが、同項目のデータであれば、運転情報記録装置に記録されるデータ以外のデータも運転情報として使用することができる。
[Overview]
In the present embodiment, the driving information recorded in the driving information recording device mounted on the railway vehicle is used to determine the probability that an error driving exceeding a fixed position to be stopped at the station will occur. Error driving that exceeds a fixed position includes overrunning at a fixed position and passing through a station without stopping. The operation information is data in which the position and speed of the train, the operation status of the brake including the number of notches, the operation status, and the like are recorded in chronological order in association with the time. The operation information is the data recorded in the operation information recording device, but if it is the data of the same item, data other than the data recorded in the operation information recording device can be used as the operation information.

具体的には、過去の所定期間における複数の駅それぞれに停止した際の運転情報に基づき、駅に停止する際の運転士の運転操作とエラー運転との関係を示す機械学習モデルである予測モデルを生成する。この予測モデルを用いることで、所与の駅に停止した際の運転情報から、当該駅において近い将来にエラー運転が行われる可能性を予測する。以下、エラー運転の要因が異なる2つの実施形態を説明する。 Specifically, it is a prediction model that is a machine learning model that shows the relationship between the driver's driving operation and error driving when stopping at a station based on the driving information when stopping at each of a plurality of stations in the past predetermined period. To generate. By using this prediction model, it is possible to predict the possibility of error driving at a given station in the near future from the driving information when the train stops at a given station. Hereinafter, two embodiments in which the factors of error operation are different will be described.

[第1実施形態]
第1実施形態は、運転エラーとして、運転士のブレーキ遅れを主因として停止位置を誤るブレーキ時機エラーが発生する可能性を予測する実施形態である。
[First Embodiment]
The first embodiment is an embodiment that predicts the possibility of a braking time error in which the stop position is erroneous, mainly due to the driver's braking delay, as a driving error.

図1は、予測モデルの生成を説明する図である。図1に示すように、第1実施形態の予測モデル10Aは、複数の教師データ310Aを用いた機械学習を行うことで生成される。1つの教師データ310Aは、1つの駅に関するデータであり、駅識別情報311Aと、駅別運転データ312と、駅設備データ313とを含む。駅識別情報311Aは、当該駅が、所定期間の間にエラー運転が発生していない第1分類駅であるか、エラー運転が発生した第2分類駅であるかを示す情報である。第1実施形態におけるエラー運転は、運転士のブレーキ遅れを主因とするブレーキ時機エラーである。駅別運転データ312は、列車が当該駅に停止した際の運転士による運転操作に関するデータであり、所定期間の間に当該駅に停止した際の複数の運転情報の統計データである。 FIG. 1 is a diagram illustrating the generation of a prediction model. As shown in FIG. 1, the prediction model 10A of the first embodiment is generated by performing machine learning using a plurality of teacher data 310A. One teacher data 310A is data related to one station, and includes station identification information 311A, station-specific operation data 312, and station equipment data 313. The station identification information 311A is information indicating whether the station is a first-class station in which error driving has not occurred or a second-class station in which error driving has occurred during a predetermined period. The error driving in the first embodiment is a braking time error mainly caused by the driver's braking delay. The station-specific operation data 312 is data related to the operation by the driver when the train stops at the station, and is statistical data of a plurality of operation information when the train stops at the station during a predetermined period.

図2は、駅別運転データ312の生成を説明する図である。図2に示すように、駅別運転データ312は、所定期間の間に当該駅に停止した際の複数の運転情報を統計演算することで生成される。1つの運転情報は、運転士の1回の運転に相当し、運転日や運転士名とともに、当該駅に停車するまでの各時刻に、列車の速度と、位置と、ブレーキノッチとを対応付けたデータである。つまり、運転情報は、当該運転士による対象駅に停車するまでの、列車の位置や速度、ブレーキノッチ等の時系列データである。複数の運転情報に対する統計演算として、例えば、ブレーキノッチの使用状況、所定時点や所定地点での速度の平均値や標準偏差を求めることで、駅別運転データ312を生成する。 FIG. 2 is a diagram illustrating the generation of operation data 312 for each station. As shown in FIG. 2, the station-specific operation data 312 is generated by statistically calculating a plurality of operation information when the station stops at the station during a predetermined period. One driving information corresponds to one driving of the driver, and associates the speed, position, and brake notch of the train with the driving date and the driver's name at each time until the train stops at the station. This is the data. That is, the driving information is time-series data such as the position and speed of the train and the brake notch until the train stops at the target station by the driver. As a statistical calculation for a plurality of driving information, for example, the usage status of the brake notch, the average value of the speeds at a predetermined time point or a predetermined point, and the standard deviation are obtained to generate the operation data 312 for each station.

駅別運転データ312は、停車前5秒間の常用最大ブレーキの使用割合[%]、停車前5秒間のブレーキ3ノッチ以上追加の割合[%]、停車前5秒間のブレーキ全緩めの使用割合[%]、停車前5秒間のブレーキノッチ移動回数の平均値[回]、停車前5秒間のブレーキノッチ移動回数の標準偏差[回]、停車前5秒間の追加ブレーキ量の平均値[ノッチ]、停車前5秒間の追加ブレーキ量の標準偏差[ノッチ]、停止位置の20m手前地点での速度の平均値[km/h]、停止位置の20m手前地点での速度の標準偏差[km/h]、停止位置の200m手前地点での速度の平均値[km/h]、停止位置の200m手前地点での速度の標準偏差[km/h]、の11個の項目を含む。 The operation data 312 by station shows the usage rate of the maximum normal brake for 5 seconds before the stop [%], the additional rate of 3 notches or more of the brake for 5 seconds before the stop [%], and the usage rate of the brake fully loosened for 5 seconds before the stop [ %], Average value of brake notch movement number for 5 seconds before stopping [times], standard deviation of brake notch movement number for 5 seconds before stopping [times], average value of additional brake amount for 5 seconds before stopping [notch], Standard deviation of the amount of additional braking for 5 seconds before stopping [notch], average value of speed 20m before the stop position [km / h], standard deviation of speed 20m before the stop position [km / h] , The average value of the speed at the point 200 m before the stop position [km / h], and the standard deviation of the speed at the point 200 m before the stop position [km / h].

図1に戻り、教師データ310Aに含まれる駅設備データ313は、列車が停止する際の運転士の運転操作に影響を与える可能性がある当該駅の設備や構造に関するデータであり、4両標の有無、6両標の有無、8両標の有無、10両標の有無、ホームにおける停止位置目標の数[個]、場内信号機の設置箇所の勾配[‰]、ホーム後端(進入側)の勾配[‰]、ホーム中央の勾配[‰]、ホーム先端(進出側)の勾配[‰]、速度パターン数[個]、の10個の項目を含む。 Returning to FIG. 1, the station equipment data 313 included in the teacher data 310A is data on the equipment and structure of the station that may affect the driving operation of the driver when the train stops, and is a four-car marker. Presence / absence of 6-car mark, presence / absence of 8-car mark, presence / absence of 10-car mark, number of stop position targets at the platform [pieces], gradient of the installation location of the on-site signal [‰], rear end of the platform (entrance side) Includes 10 items: gradient [‰] in the center of the platform, gradient [‰] in the center of the platform, gradient [‰] at the tip of the platform (advanced side), and number of speed patterns [pieces].

“n両標”とは、n両編成の列車の停止位置目標を示す標識のことである。このn両標は、列車の編成車両数に応じて停止位置目標が異なる駅に設けられる。つまり、4両標の有無、6両標の有無、8両標の有無、10両標の有無、及び、ホームにおける停止位置目標の数[個]の項目は、定位置である停止位置目標が編成車両数に応じて異なるか否かに関する定位置可変情報に相当する。 The "n-car mark" is a sign indicating a stop position target of an n-car train. These n-car markers are provided at stations with different stop position targets depending on the number of trains. That is, the presence / absence of the 4-car mark, the presence / absence of the 6-car mark, the presence / absence of the 8-car mark, the presence / absence of the 10-car mark, and the number [pieces] of the stop position targets at the platform are the fixed positions of the stop position targets. It corresponds to the fixed position variable information regarding whether or not it differs depending on the number of trains.

“速度パターン数”とは、当該駅に停止する際に取り得る大まかな速度パターン(位置に対する列車の速度のグラフ)の種類数である。例えば、当該駅に停止する際の複数の運転情報それぞれから、駅構内に進入してから停止するまでの速度パターンを作成する。そして、作成した複数の速度パターンをパターンの形状が類似するグループに分類し、分類したグループ数を速度パターン数とする。速度パターンの形状は、主に、信号現示に応じて決まるため、速度パターンの数は、場内信号機を含む信号機の数及びその信号現示に応じて決まることになる。例えば、信号機が場内信号機のみであり、且つ、場内信号機が取り得る現示が、青のみ、青+黄、黄のみ、黄+黄の4種類ならば、速度パターン数は「4」となる。 The "number of speed patterns" is the number of types of rough speed patterns (graph of train speed with respect to position) that can be taken when stopping at the station. For example, a speed pattern from entering the station yard to stopping is created from each of a plurality of driving information when stopping at the station. Then, the created plurality of speed patterns are classified into groups having similar pattern shapes, and the number of classified groups is defined as the number of speed patterns. Since the shape of the speed pattern is mainly determined according to the signal display, the number of speed patterns is determined according to the number of traffic lights including the on-site signal and the signal display. For example, if the traffic light is only the on-site traffic light and the indications that the on-site traffic light can take are four types of blue only, blue + yellow, yellow only, and yellow + yellow, the number of speed patterns is "4".

予測モデル10Aは、教師データ310Aに含まれるデータのうち、駅別運転データ312及び駅設備データ313を入力とし、駅識別情報311Aを出力とした機械学習を行うことで生成する。 The prediction model 10A is generated by performing machine learning with station-specific operation data 312 and station equipment data 313 as inputs and station identification information 311A as output among the data included in the teacher data 310A.

本実施形態では、機械学習の種類はロジスティック回帰分析とするが、ニューラルネットワーク等の他の機械学習であってもよい。つまり、教師データ310Aの駅別運転データ312及び駅設備データ313を説明変数とし、駅識別情報311Aを目的変数としたロジスティック回帰分析を行うことで、予測モデル10Aを生成する。説明変数は、駅別運転データ312である11個の項目と、駅設備データ313である10個の項目との合計21個の項目とする。目的変数は、第1分類駅であることを「0」、第2分類駅であることを「1」とする二値変数である。ロジスティック回帰分析の結果となる予測モデル式は、式(1)に示すモデル式とする。

Figure 2021142873
In the present embodiment, the type of machine learning is logistic regression analysis, but other machine learning such as a neural network may be used. That is, the prediction model 10A is generated by performing logistic regression analysis using the station-specific operation data 312 and the station equipment data 313 of the teacher data 310A as explanatory variables and the station identification information 311A as the objective variable. The explanatory variables are a total of 21 items, including 11 items of station-specific operation data 312 and 10 items of station equipment data 313. The objective variable is a binary variable in which the first classification station is "0" and the second classification station is "1". The prediction model formula that is the result of the logistic regression analysis is the model formula shown in formula (1).
Figure 2021142873

式(1)に示すモデル式は、21個の説明変数のうち、停車前5秒間の常用最大ブレーキの使用割合[%](説明変数Xとする)、停止位置の200m手前地点での速度の平均値[km/h](説明変数Xとする)、停止位置の200m手前地点での速度の標準偏差[km/h](説明変数Xとする)、速度パターン数[個](説明変数Xとする)、の4個の説明変数X〜Xが、目的変数に影響を与えていることを示している。92箇所の駅において行われた1年間の運転情報をもとに本実施形態の機械学習の実験を行った結果、一例として、説明変数X〜Xの係数a〜a、及び、定数項aは、例えば、a=3.69,a=0.37,a=−0.72,a=0.90,a=−22.79、であった。式(1)のモデル式により算出される予測値Pは、「0.0〜1.0」の範囲の数値であり、駅が第2分類駅である可能性、つまり、エラー運転が発生する確率を表す。 In the model formula shown in the formula (1), out of the 21 explanatory variables, the usage rate [%] (explanatory variable X 1 ) of the maximum regular brake for 5 seconds before the stop, and the speed at a point 200 m before the stop position. [Km / h] (explanatory variable X 2 ), standard deviation of velocity [km / h] (explanatory variable X 3 ) at a point 200 m before the stop position, number of velocity patterns [pieces] ( the explanatory variable X 4), four explanatory variables X 1 to X 4 of the shows that are affecting the objective variable. As a result of conducting the machine learning experiment of this embodiment based on the one-year operation information conducted at 92 stations, as an example, the coefficients a 1 to a 4 of the explanatory variables X 1 to X 4 and the coefficients a 1 to a 4 of the explanatory variables X 1 to X 4 and The constant term a 0 was, for example, a 1 = 3.69, a 2 = 0.37, a 3 = −0.72, a 4 = 0.90, a 0 = -22.79. The predicted value P calculated by the model formula of the formula (1) is a numerical value in the range of "0.0 to 1.0", and there is a possibility that the station is a second classification station, that is, an error operation occurs. Represents a probability.

図3は、予測モデル10Aを用いた、所与の駅においてエラー運転が発生する蓋然性の判定を説明する図である。図3に示すように、判定対象の駅(図3の例では、“Z駅”)に関するデータである駅データ330を、予測モデル10Aに入力する。駅データ330は、駅別運転データと、駅設備データとを含む。駅別運転データに含まれる、説明変数Xである停車前5秒間の常用最大ブレーキの使用割合[%]、説明変数Xである停止位置の200m手前地点での速度の平均値[km/h]、説明変数Xである停止位置の200m手前地点での速度の標準偏差[km/h]と、駅設備データに含まれる、説明変数Xである速度パターン数と、の4個の説明変数X〜Xの値を、式(1)で表されるモデル式に入力して、予測値Pを算出する。 FIG. 3 is a diagram for explaining the determination of the probability that error driving will occur at a given station using the prediction model 10A. As shown in FIG. 3, station data 330, which is data related to the station to be determined (“Z station” in the example of FIG. 3), is input to the prediction model 10A. The station data 330 includes station-specific operation data and station equipment data. Explanatory variable X 1 included in the operation data for each station, the usage rate [%] of the maximum regular brake for 5 seconds before the stop, and explanatory variable X 2 the average value of the speed 200 m before the stop position [km / h], the explanatory variables X 3 is a stop position velocity and standard deviation [km / h] of at 200m before the point, contained in the station facility data, the speed pattern number and an explanatory variable X 4, 4 pieces of the The values of the explanatory variables X 1 to X 4 are input to the model formula represented by the formula (1) to calculate the predicted value P.

予測値Pは、第2分類駅である可能性を示すので、予測値Pが「1」に近いほど、第2分類駅に該当する可能性が高く、逆に、「0」に近いほど、第1分類駅に該当する可能性が高いといえる。従って、例えば予測値Pの閾値を「0.5」として、予測値Pが「0.5」未満ならば、判定対象の駅は第1分類駅に該当し、予測値Pが「0.5」以上ならば、判定対象の駅は第2分類駅に該当すると判定することができる。勿論、2段階で判定するのではなく、予測値Pが「0」に近いか「1」に近いかで、第1分類駅に該当する可能性が高い/低い、第2分類駅に該当する可能性が高い/低い、どちらともいえない、の計5段階で判定することとしてもよい。 Since the predicted value P indicates the possibility of being a second-class station, the closer the predicted value P is to "1", the more likely it is to correspond to the second-class station, and conversely, the closer it is to "0", the more likely it is. It can be said that there is a high possibility that it corresponds to the first category station. Therefore, for example, if the threshold value of the predicted value P is "0.5" and the predicted value P is less than "0.5", the station to be determined corresponds to the first classification station, and the predicted value P is "0.5". If it is the above, it can be determined that the station to be determined corresponds to the second category station. Of course, instead of judging in two stages, it corresponds to the second classification station, which is likely / low to correspond to the first classification station depending on whether the predicted value P is close to "0" or "1". It may be judged in a total of 5 stages of high / low possibility and neither.

そして、第2分類駅に該当すると判定したけれども、所定の対象期間の間にエラー運転が発生していないならば、判定対象の駅は、今後エラー運転が発生する可能性が高い、と判定する。これは、過去の運転情報からは第2分類駅に該当する、つまり、エラー運転が発生する可能性が高い駅と予測されるが、所定期間の間はエラー運転が発生していない駅であり、この場合には、今後エラー運転が発生する可能性が高いと考えられるからである。なお、予測値Pが「1」に近いほど今後エラー運転が発生する可能性が高いと判定して、その可能性の度合を予測することとしてもよい。 Then, if it is determined that the station corresponds to the second classification station, but no error driving has occurred during the predetermined target period, it is determined that the station to be determined is likely to have an error driving in the future. .. This is a station that corresponds to the second category station from the past driving information, that is, a station that is likely to cause error driving, but no error driving occurs during a predetermined period. In this case, it is highly likely that error operation will occur in the future. It should be noted that the closer the predicted value P is to "1", the higher the possibility that error driving will occur in the future, and the degree of the possibility may be predicted.

図4は、第1実施形態の判定方法による判定結果(実験結果)の一例である。92箇所の駅それぞれに停止した際の運転情報を教師データ310Aとして、ロジスティック回帰分析により、式(1)に示したモデル式で表される予測モデル10Aを生成した。92箇所の駅の内訳は、第1分類駅が46箇所、第2分類駅が46箇所である。次いで、教師データ310Aとして用いた92箇所の駅それぞれを判定対象の駅として、その駅データ330を生成した予測モデル10Aに入力して予測値Pを算出した。 FIG. 4 is an example of a determination result (experimental result) by the determination method of the first embodiment. Using the driving information when the trains stopped at each of the 92 stations as teacher data 310A, a prediction model 10A represented by the model formula shown in the formula (1) was generated by logistic regression analysis. Of the 92 stations, there are 46 first-class stations and 46 second-class stations. Next, each of the 92 stations used as the teacher data 310A was set as the station to be determined, and the station data 330 was input to the generated prediction model 10A to calculate the predicted value P.

図4は、このように算出した92箇所の駅それぞれの予測値Pをプロットした図である。左側に、第1分類駅の予測値Pのプロットを示し、右側に、第2分類駅の予測値Pのプロットを示している。縦軸は、予測値Pである。また、第1分類駅のうち、所定期間以降に最初にエラー運転が発生した1つの駅についてのみ、“白抜きの三角”のプロットで示している。 FIG. 4 is a diagram in which the predicted values P of each of the 92 stations calculated in this way are plotted. The plot of the predicted value P of the first classification station is shown on the left side, and the plot of the predicted value P of the second classification station is shown on the right side. The vertical axis is the predicted value P. Further, among the first-class stations, only one station in which the error operation first occurs after a predetermined period is shown by the plot of "white triangle".

図4によれば、第1分類駅の8割弱の予測値Pが「0.5」未満となっており、第1分類駅に該当すると予測されている。また、第1分類駅のうち、所定期間以降に最初にエラー運転が発生した駅の予測値Pは「0.59」であり、第2分類駅に該当すると予測されている。また、第2分類駅の8割弱の予測値Pが「0.5」以上となっており、第2分類駅に該当すると予測されている。 According to FIG. 4, the predicted value P of less than 80% of the first-class stations is less than "0.5", and it is predicted that the stations correspond to the first-class stations. Further, among the first-class stations, the predicted value P of the station where the error driving first occurs after the predetermined period is "0.59", and it is predicted that the station corresponds to the second-class station. In addition, the predicted value P of less than 80% of the second classification station is "0.5" or more, and it is predicted that it corresponds to the second classification station.

この判定例のように、教師データ310として予測モデル10Aの生成に用いた駅を、当該予測モデル10Aを用いた判定対象の駅とすることで、これらの駅の中から、今後、運転士のブレーキ遅れを主因とするエラー運転が発生する可能性が高い駅を抽出することができる。つまり、第1分類駅であるが、予測モデル10Aから予測される予測値が「1」に近い(例えば、予測値Pが「0.5」以上)駅を、今後エラー運転が発生する可能性が高い駅として抽出することができる。 As in this determination example, the station used to generate the prediction model 10A as the teacher data 310 is set as the station to be determined using the prediction model 10A. It is possible to extract stations where error driving is likely to occur mainly due to brake delay. That is, although it is the first classification station, there is a possibility that an error operation will occur in the future at a station where the predicted value predicted from the prediction model 10A is close to "1" (for example, the predicted value P is "0.5" or more). Can be extracted as a high station.

図5は、判定装置1Aの機能構成図である。この判定装置1Aは、上述したように、対象駅に停止した際の運転情報を用いて、駅に停止する際の運転士の運転操作とエラー運転との関係を示す予測モデル10Aを生成し、この予測モデル10Aを用いて、所与の判定対象の駅においてエラー運転が発生する蓋然性を判定する装置である。図5に示すように、判定装置1Aは、入力部102と、表示部104と、音出力部106と、通信部108と、処理部200と、記憶部300とを備えて構成され、一種のコンピュータシステムとして実現される。なお、判定装置1Aは、1台のコンピュータで実現してもよいし、複数台のコンピュータを接続して構成することとしてもよい。 FIG. 5 is a functional configuration diagram of the determination device 1A. As described above, the determination device 1A uses the driving information when the train stops at the target station to generate a prediction model 10A showing the relationship between the driver's driving operation and the error driving when the train stops at the target station. Using this prediction model 10A, it is a device that determines the probability that an error operation will occur at a given station to be determined. As shown in FIG. 5, the determination device 1A includes an input unit 102, a display unit 104, a sound output unit 106, a communication unit 108, a processing unit 200, and a storage unit 300, and is a kind of determination device 1A. Realized as a computer system. The determination device 1A may be realized by one computer, or may be configured by connecting a plurality of computers.

入力部102は、例えばキーボードやマウス、タッチパネル、各種スイッチ等の入力装置で実現され、なされた操作に応じた操作信号を処理部200に出力する。表示部104は、例えば液晶ディスプレイやタッチパネル等の表示装置で実現され、処理部200からの表示信号に基づく各種表示を行う。音出力部106は、例えばスピーカ等の音出力装置で実現され、処理部200からの音信号に基づく各種音出力を行う。通信部108は、例えば無線通信モジュールやルータ、モデム、有線用の通信ケーブルのジャックや制御回路等で実現される通信装置であり、所与の通信ネットワークに接続して外部装置とのデータ通信を行う。 The input unit 102 is realized by an input device such as a keyboard, a mouse, a touch panel, and various switches, and outputs an operation signal corresponding to the performed operation to the processing unit 200. The display unit 104 is realized by a display device such as a liquid crystal display or a touch panel, and performs various displays based on the display signal from the processing unit 200. The sound output unit 106 is realized by a sound output device such as a speaker, and outputs various sounds based on the sound signal from the processing unit 200. The communication unit 108 is a communication device realized by, for example, a wireless communication module, a router, a modem, a jack or a control circuit of a wired communication cable, and is connected to a given communication network to perform data communication with an external device. conduct.

処理部200は、CPU(Central Processing Unit)やFPGA(Field Programmable Gate Array)等の演算装置や演算回路で実現されるプロセッサーであり、記憶部300に記憶されたプログラムやデータ、入力部102や通信部108からの入力データ等に基づいて、判定装置1Aの全体制御を行う。また、処理部200は、機能的な処理ブロックとして、教師データ取得部202、予測モデル生成部204、駅データ取得部206、予測部208、蓋然性判定部210、を有する。処理部200が有するこれらの各機能部は、処理部200がプログラムを実行することでソフトウェア的に実現することも、専用の演算回路で実現することも可能である。本実施形態では、前者のソフトウェア的に実現することとして説明する。 The processing unit 200 is a processor realized by an arithmetic unit such as a CPU (Central Processing Unit) or an FPGA (Field Programmable Gate Array) or an arithmetic circuit, and is a program or data stored in the storage unit 300, an input unit 102, or a communication unit. Overall control of the determination device 1A is performed based on the input data from the unit 108 and the like. Further, the processing unit 200 has a teacher data acquisition unit 202, a prediction model generation unit 204, a station data acquisition unit 206, a prediction unit 208, and a probability determination unit 210 as functional processing blocks. Each of these functional units included in the processing unit 200 can be realized by software by the processing unit 200 executing a program, or can be realized by a dedicated arithmetic circuit. In the present embodiment, the former will be described as being realized by software.

教師データ取得部202は、停止すべき駅の定位置を超過したエラー運転が所定期間の間に発生していない第1分類駅と、発生した第2分類駅とについて、当該所定期間の間に当該駅に停止する際に行われた複数の運転情報の統計データである駅別運転データ312を含む教師データ310Aを取得する。教師データ310Aは、駅別運転データ312のほか、駅識別情報311Aと、駅設備データ313とを含む。駅識別情報311Aが示すエラー運転は、運転士のブレーキ遅れを主因として停止位置を誤るブレーキ時機エラーである。駅別運転データ312には、停車前5秒間の常用最大ブレーキの使用割合(説明変数X)、停止位置の200m手前地点での速度の平均値(説明変数X)、定位置から所定距離手前の速度のバラツキに関する情報である停止位置の200m手前地点での速度の標準偏差(説明変数X)が含まれ、駅設備データ313には、当該駅の定位置である停止位置目標に停止するまでの速度パターンの種類数(説明変数X)が含まれる(図1参照)。 The teacher data acquisition unit 202 describes the first-class station in which the error operation exceeding the fixed position of the station to be stopped has not occurred during the predetermined period and the second-class station in which the error operation has occurred during the predetermined period. The teacher data 310A including the station-specific driving data 312, which is statistical data of a plurality of driving information performed when the train stops at the station, is acquired. The teacher data 310A includes station identification information 311A and station equipment data 313, in addition to station-specific operation data 312. The error driving indicated by the station identification information 311A is a braking time error in which the stop position is erroneous mainly due to the driver's braking delay. The operation data 312 by station includes the usage ratio of the maximum regular brake for 5 seconds before the stop (explanatory variable X 1 ), the average value of the speed 200 m before the stop position (explanatory variable X 2 ), and the predetermined distance from the fixed position. The standard deviation of the speed (explanatory variable X 3 ) at the point 200 m before the stop position, which is information on the variation of the speed in front, is included, and the station equipment data 313 stops at the stop position target which is the fixed position of the station. The number of types of speed patterns (explanatory variable X 4 ) is included (see FIG. 1).

予測モデル生成部204は、駅別運転データと、駅設備データと、当該駅別運転データ及び駅設備データが第1分類駅の情報か前記第2分類駅の情報かを示す駅識別情報とを教師データとして、所与の駅の駅別運転データが入力されたときに、当該所与の駅が第1分類駅又は第2分類駅に該当する可能性を示す予測値を出力する機械学習モデルを生成する。 The prediction model generation unit 204 generates station-specific operation data, station equipment data, and station identification information indicating whether the station-specific operation data and station equipment data are information on a first-class station or information on a second-class station. A machine learning model that outputs predicted values indicating the possibility that a given station corresponds to a first-class station or a second-class station when the station-specific operation data of a given station is input as teacher data. To generate.

つまり、予測モデル生成部204は、教師データ取得部202により取得された教師データ310Aを用いた機械学習を行って、所与の判定対象の駅の駅別運転データ及び駅設備データを含む駅データ330が入力されたときに、当該判定対象の駅が第1分類駅又は第2分類駅に該当する可能性を示す予測値Pを出力する機械学習モデルである予測モデル10Aを生成する(図1参照)。機械学習としてロジスティック回帰分析を行った場合には、予測モデル10Aとして、式(1)に一例を示すモデル式が生成される。予測モデル生成部204が生成した予測モデル10Aに関するデータは、予測モデルデータ320として記憶される。予測モデルデータ320は、予測モデル10Aを定義するデータであり、具体的には、式(1)に一例を示すモデル式の係数a〜aのデータである(図1参照)。 That is, the prediction model generation unit 204 performs machine learning using the teacher data 310A acquired by the teacher data acquisition unit 202, and station data including station-specific operation data and station equipment data of a given determination target station. When 330 is input, a prediction model 10A, which is a machine learning model that outputs a prediction value P indicating the possibility that the station to be determined corresponds to the first classification station or the second classification station, is generated (FIG. 1). reference). When logistic regression analysis is performed as machine learning, a model formula showing an example in formula (1) is generated as the prediction model 10A. The data related to the prediction model 10A generated by the prediction model generation unit 204 is stored as the prediction model data 320. The prediction model data 320 is data that defines the prediction model 10A, and specifically, is data of coefficients a 0 to a 4 of the model formula showing an example in the formula (1) (see FIG. 1).

駅データ取得部206は、所与の駅に係る駅別運転データ及び駅設備データを含む駅データ330を取得する(図3参照)。 The station data acquisition unit 206 acquires station data 330 including station-specific operation data and station equipment data related to a given station (see FIG. 3).

予測部208は、機械学習モデルに、所与の駅に係る駅別運転データ及び駅設備データを入力して、当該所与の駅の予測値を得る予測処理を実行する。つまり、予測部208は、駅データ取得部206により取得された駅データ330を、予測モデル生成部204により予測モデル10Aとして生成されたモデル式(1)に入力することで、判定対象の駅が第1分類駅又は第2分類駅に該当する可能性を示す予測値Pを算出する。 The prediction unit 208 inputs station-specific operation data and station equipment data related to a given station into the machine learning model, and executes a prediction process for obtaining a predicted value of the given station. That is, the prediction unit 208 inputs the station data 330 acquired by the station data acquisition unit 206 into the model formula (1) generated as the prediction model 10A by the prediction model generation unit 204, so that the station to be determined can be determined. The predicted value P indicating the possibility of corresponding to the first classification station or the second classification station is calculated.

蓋然性判定部210は、予測処理で得られた予測値に基づいて、所与の駅においてエラー運転が発生する蓋然性を判定する。つまり、蓋然性判定部210は、予測部208により算出された予測値Pが閾値(例えば、「0.5」)以上であり、且つ、判定対象の駅において所定期間の間にエラー運転が発生していないならば、判定対象の駅において近い将来にエラー運転が発生する可能性が高い、と判定する。 The probability determination unit 210 determines the probability that an error driving will occur at a given station based on the predicted value obtained in the prediction process. That is, in the probability determination unit 210, the predicted value P calculated by the prediction unit 208 is equal to or higher than the threshold value (for example, “0.5”), and an error operation occurs at the station to be determined during a predetermined period. If not, it is determined that there is a high possibility that an error operation will occur at the station to be determined in the near future.

記憶部300は、ROM(Read Only Memory)やRAM(Random Access Memory)等のIC(Integrated Circuit)メモリやハードディスク等の記憶装置で実現され、処理部200が判定装置1Aを統合的に制御するためのプログラムやデータ等を記憶しているとともに、処理部200の作業領域として用いられ、処理部200が実行した演算結果や、入力部102や通信部108からの入力データ等が一時的に格納される。本実施形態では、記憶部300には、判定プログラム302と、教師データ取得部202により取得された教師データ310Aと、予測モデル生成部204により生成された予測モデルデータ320Aと、駅データ取得部206により取得された駅データ330と、蓋然性判定部210より判定された判定結果データ340Aと、が記憶される。 The storage unit 300 is realized by a storage device such as an IC (Integrated Circuit) memory such as a ROM (Read Only Memory) or a RAM (Random Access Memory) or a hard disk, and the processing unit 200 controls the determination device 1A in an integrated manner. Is stored as a work area of the processing unit 200, and the calculation result executed by the processing unit 200 and the input data from the input unit 102 and the communication unit 108 are temporarily stored. NS. In the present embodiment, the storage unit 300 contains the determination program 302, the teacher data 310A acquired by the teacher data acquisition unit 202, the prediction model data 320A generated by the prediction model generation unit 204, and the station data acquisition unit 206. The station data 330 acquired by the above and the determination result data 340A determined by the probability determination unit 210 are stored.

[第2実施形態]
第2実施形態は、運転エラーとして、運転士の停止位置の勘違いを主因として停止位置を誤る停目両数エラーが行われる可能性を予測する。なお、以下の第2実施形態においては、第1実施形態と異なる点を中心に説明する。また、上述の第1実施形態と同一の構成要素については同一符号を付して詳細な説明を省略する。
[Second Embodiment]
In the second embodiment, as a driving error, it is predicted that there is a possibility that a stop position error may occur in which the stop position is erroneous mainly due to the driver's misunderstanding of the stop position. In the second embodiment below, the points different from those of the first embodiment will be mainly described. Further, the same components as those in the above-described first embodiment are designated by the same reference numerals, and detailed description thereof will be omitted.

図6は、予測モデルの生成を説明する図である。図6に示すように、第2実施形態の予測モデル10Bは、第1実施形態と同様に、複数の教師データ310Bを用いた機械学習を行うことで生成される。教師データ310Bに含まれる駅識別情報311Bが示すエラー運転は、運転士の停止位置の勘違いを主因として停止位置を誤る停目両数エラーであり、第1実施形態におけるブレーキ時機エラーとは異なる。 FIG. 6 is a diagram illustrating the generation of the prediction model. As shown in FIG. 6, the prediction model 10B of the second embodiment is generated by performing machine learning using a plurality of teacher data 310B as in the first embodiment. The error driving indicated by the station identification information 311B included in the teacher data 310B is a stop position error in which the stop position is erroneous mainly due to the driver's misunderstanding of the stop position, and is different from the brake time error in the first embodiment.

予測モデル10Bは、教師データ310Bに含まれるデータのうち、駅別運転データ312及び駅設備データ313を入力とし、駅識別情報311Bを出力とした機械学習を行うことで生成する。つまり、教師データ310の駅別運転データ312及び駅設備データ313を説明変数とし、駅識別情報311Bを目的変数としたロジスティック回帰分析を行うことで、予測モデル10Bを生成する。ロジスティック回帰分析の結果となる予測モデル式は、式(2)に示すモデル式とする。

Figure 2021142873
The prediction model 10B is generated by performing machine learning with station-specific operation data 312 and station equipment data 313 as inputs and station identification information 311B as output among the data included in the teacher data 310B. That is, the prediction model 10B is generated by performing logistic regression analysis using the station-specific operation data 312 and the station equipment data 313 of the teacher data 310 as explanatory variables and the station identification information 311B as the objective variable. The prediction model formula that is the result of the logistic regression analysis is the model formula shown in formula (2).
Figure 2021142873

式(2)に示すモデル式は、21個の説明変数のうち、停止位置の20m手前地点での速度の平均値[km/h](説明変数Xとする)、4両標の有無(説明変数Xとする)、6両標の有無(説明変数Xとする)、速度パターン数[個](説明変数Xとする)、の4個の説明変数X〜Xが、目的変数に影響を与えていることを示している。72箇所の駅において行われた1年間の運転情報をもとに本実施形態の機械学習の実験を行った結果、一例として、説明変数X〜Xの係数a〜a、及び、定数項aは、例えば、a=0.80,a=1.42,a=1.15,a=1.07,a=−18.30、であった。式(2)のモデル式により算出される予測値Pは、「0.0〜1.0」の範囲の数値であり、駅が第2分類駅である可能性、つまり、エラー運転が発生する確率を表す。第2実施形態は、主因が異なるエラー運転の発生可能性を予測するものであり、教師データ310Bに含まれる駅識別情報311Bが第1実施形態と異なるから、第1実施形態の予測モデル10Aとは異なる予測モデル10Bが生成される。 Model formula shown in Equation (2), among the 21 pieces of explanatory variables, the average value of the velocity at 20m before the point of the stop position [km / h] (as explanatory variables X 5), whether the four-car mark ( the explanatory variable X 6), the presence or absence of 6-Car mark (the explanatory variable X 7), the speed pattern number [pieces] (explanatory variable X 8), 4 pieces of explanatory variables X 5 to X 8 of, It shows that it affects the objective variable. As a result of conducting the machine learning experiment of this embodiment based on the one-year operation information conducted at 72 stations, as an example, the coefficients a 5 to a 8 of the explanatory variables X 5 to X 8 and the coefficients a 5 to a 8 of the explanatory variables X 5 to X 8 and The constant term a 9 was, for example, a 5 = 0.80, a 6 = 1.42, a 7 = 1.15, a 8 = 1.07, a 9 = -18.30. The predicted value P calculated by the model formula of the formula (2) is a numerical value in the range of "0.0 to 1.0", and there is a possibility that the station is a second classification station, that is, an error operation occurs. Represents a probability. The second embodiment predicts the possibility of error driving having different main causes, and since the station identification information 311B included in the teacher data 310B is different from the first embodiment, it is different from the prediction model 10A of the first embodiment. Generates a different prediction model 10B.

図7は、予測モデル10Bを用いた、所与の駅においてエラー運転が発生する蓋然性の判定を説明する図である。図7に示すように、判定対象の駅(図7の例では、“Z駅”)に関するデータである駅データ330を、予測モデル10Bに入力する。つまり、駅データ330の駅別運転データに含まれる、説明変数Xである停止位置の20m手前地点での速度の平均値[km/h]と、駅設備データに含まれる、説明変数Xである4両標の有無と、説明変数Xである6両標の有無と、説明変数Xである速度パターン数と、の4個の説明変数X〜Xの値を、式(2)で表されるモデル式に入力して、予測値Pを算出する。 FIG. 7 is a diagram for explaining the determination of the probability that error driving will occur at a given station using the prediction model 10B. As shown in FIG. 7, station data 330, which is data related to the station to be determined (“Z station” in the example of FIG. 7), is input to the prediction model 10B. That is, the average value [km / h] of the speed at the point 20 m before the stop position, which is the explanatory variable X 5 included in the station data 330 for each station , and the explanatory variable X 6 included in the station equipment data. and the presence of four-car mark is, the presence of six-car mark an explanatory variable X 7, speed pattern number and an explanatory variable X 8, the value of four explanatory variables X 5 to X 8 in the formula ( The predicted value P is calculated by inputting into the model formula represented by 2).

そして、第1実施形態と同様に、例えば予測値Pの閾値を「0.5」として、予測値Pが「0.5」未満ならば、判定対象の駅は第1分類駅に該当し、予測値Pが「0.5」以上ならば、判定対象の駅は第2分類駅に該当すると判定する。 Then, as in the first embodiment, for example, if the threshold value of the predicted value P is "0.5" and the predicted value P is less than "0.5", the station to be determined corresponds to the first classification station. If the predicted value P is "0.5" or more, it is determined that the station to be determined corresponds to the second category station.

図8は、第2実施形態の判定方法による判定結果(実験結果)の一例である。72箇所の駅それぞれに停止した際の運転情報を教師データ310Bとして、ロジスティック回帰分析により、式(2)に示したモデル式で表される予測モデル10Bを生成した。72箇所の駅の内訳は、第1分類駅が36箇所、第2分類駅が36箇所である。次いで、教師データ310Aとして用いた72箇所の駅それぞれを判定対象の駅として、その駅データ330を生成した予測モデル10Bに入力して予測値Pを算出した。 FIG. 8 is an example of a determination result (experimental result) by the determination method of the second embodiment. Using the driving information when the trains stopped at each of the 72 stations as the teacher data 310B, the prediction model 10B represented by the model formula shown in the formula (2) was generated by logistic regression analysis. Of the 72 stations, there are 36 first-class stations and 36 second-class stations. Next, each of the 72 stations used as the teacher data 310A was set as the station to be determined, and the station data 330 was input to the generated prediction model 10B to calculate the predicted value P.

図8は、このように算出した72箇所の駅それぞれの予測値Pをプロットした図である。左側に、第1分類駅の予測値Pのプロットを示し、右側に、第2分類駅の予測値Pのプロットを示している。縦軸は、予測値Pである。また、第1分類駅のうち、所定期間以降に最初にエラー運転が発生した1つの駅についてのみ、“白抜きの三角”のプロットで示している。 FIG. 8 is a diagram in which the predicted values P of each of the 72 stations calculated in this way are plotted. The plot of the predicted value P of the first classification station is shown on the left side, and the plot of the predicted value P of the second classification station is shown on the right side. The vertical axis is the predicted value P. Further, among the first-class stations, only one station in which the error operation first occurs after a predetermined period is shown by the plot of "white triangle".

図8によれば、第1分類駅の8割弱の予測値Pが「0.5」未満となっており、第1分類駅に該当すると予測されている。また、第1分類駅のうち、所定期間以降に最初にエラー運転が発生した駅の予測値Pは「0.71」であり、第2分類駅に該当すると予測されている。また、第2分類駅の8割弱の予測値Pが「0.5」以上であり、第2分類駅に該当すると予測されている。 According to FIG. 8, the predicted value P of less than 80% of the first-class stations is less than “0.5”, and it is predicted that the stations correspond to the first-class stations. Further, among the first-class stations, the predicted value P of the station where the error driving first occurs after the predetermined period is "0.71", and it is predicted that the station corresponds to the second-class station. In addition, the predicted value P of less than 80% of the second classification station is "0.5" or more, and it is predicted that it corresponds to the second classification station.

この判定例のように、教師データ310として予測モデル10Bの生成に用いた駅を、当該予測モデル10Bを用いた判定対象の駅とすることで、これらの駅の中から、今後、運転士の停止位置目標の勘違いを主因とするエラー運転が発生する可能性が高い駅を抽出することができる。つまり、第1分類駅であるが、予測モデル10Bから予測される予測値が「1」に近い(例えば、予測値Pが「0.5」以上)駅を、今後エラー運転が発生する可能性が高い駅として抽出することができる。 As in this determination example, the station used to generate the prediction model 10B as the teacher data 310 is set as the station to be determined using the prediction model 10B. It is possible to extract stations that are likely to cause error driving mainly due to misunderstanding of the stop position target. That is, although it is the first classification station, there is a possibility that an error operation will occur in the future at a station where the predicted value predicted from the prediction model 10B is close to "1" (for example, the predicted value P is "0.5" or more). Can be extracted as a high station.

図9は、判定装置1Bの機能構成図である。この判定装置1Bは、上述したように、対象駅に停止した際の運転情報を用いて、駅に停止する際の運転士の運転操作とエラー運転との関係を示す予測モデル10Bを生成し、この予測モデル10Bを用いて、所与の判定対象の駅においてエラー運転が発生する蓋然性を判定する装置である。 FIG. 9 is a functional configuration diagram of the determination device 1B. As described above, the determination device 1B uses the driving information when the train stops at the target station to generate a prediction model 10B showing the relationship between the driver's driving operation and the error driving when the train stops at the target station. Using this prediction model 10B, it is a device that determines the probability that an error operation will occur at a given station to be determined.

判定装置1Bにおいて、教師データ取得部202は、教師データ310Bを取得する。教師データ310Bに含まれる駅識別情報311Bが示すエラー運転は、運転士の停止位置の勘違いを主因として停止位置を誤る停目両数エラーである。 In the determination device 1B, the teacher data acquisition unit 202 acquires the teacher data 310B. The error driving indicated by the station identification information 311B included in the teacher data 310B is a stop position error in which the stop position is erroneous mainly due to the driver's misunderstanding of the stop position.

予測モデル生成部204は、教師データ取得部202により取得された教師データ310Bを用いた機械学習を行って、予測モデル10Bを生成する(図6参照)。機械学習としてロジスティック回帰分析を行った場合には、予測モデル10Bとして、式(2)に一例を示すモデル式が生成される。予測モデル生成部204が生成した予測モデル10Bに関するデータは、予測モデルデータ320Bとして記憶される。予測モデルデータ320Bは、予測モデル10Bを定義するデータであり、具体的には、式(2)に一例を示すモデル式の係数a〜aのデータである。 The prediction model generation unit 204 performs machine learning using the teacher data 310B acquired by the teacher data acquisition unit 202 to generate the prediction model 10B (see FIG. 6). When logistic regression analysis is performed as machine learning, a model formula showing an example in formula (2) is generated as the prediction model 10B. The data related to the prediction model 10B generated by the prediction model generation unit 204 is stored as the prediction model data 320B. Predictive model data 320B is data defining the predictive model 10B, specifically, the data of the coefficient a 5 ~a 9 model type, one example of which is shown in Equation (2).

予測部208は、駅データ取得部206により取得された駅データ330を、予測モデル生成部204により予測モデル10Bとして生成されたモデル式(2)に入力することで、判定対象の駅が第1分類駅又は第2分類駅に該当する可能性を示す予測値Pを算出する。 The prediction unit 208 inputs the station data 330 acquired by the station data acquisition unit 206 into the model formula (2) generated as the prediction model 10B by the prediction model generation unit 204, so that the station to be determined is the first station. The predicted value P indicating the possibility of corresponding to the classification station or the second classification station is calculated.

[作用効果]
このように、第1実施形態及び第2実施形態によれば、記録・蓄積された運転情報を用いて、個々の駅における将来のエラー事象の発生可能性を予測することが可能となる。つまり、エラー事象として停止すべき駅の定位置を超過したエラー運転が発生していない第1分類駅と、発生した第2分類駅とのそれぞれの駅別運転データを教師データとして、所与の駅の駅別運転データを入力すると、当該駅が第1分類駅又は第2分類駅に該当する可能性を示す予測値を出力する機械学習モデルを生成することができる。機械学習モデルが出力する予測値は、所与の駅に停止する際に行われた運転の仕方が、エラー運転が発生した駅と発生していない駅とのどちらの駅において行われた運転の仕方にどの程度似ているかを示している。従って、機械学習モデルが出力する予測値から、所与の駅において、エラー運転が発生する可能性を予測することができる。
[Action effect]
As described above, according to the first embodiment and the second embodiment, it is possible to predict the possibility of future error events occurring at individual stations by using the recorded and accumulated operation information. That is, given as teacher data, the operation data for each station of the first category station where the error operation exceeding the fixed position of the station to be stopped as an error event has not occurred and the second category station where the error operation has occurred is given. By inputting the station-specific operation data of a station, it is possible to generate a machine learning model that outputs a predicted value indicating the possibility that the station corresponds to a first-class station or a second-class station. The predicted value output by the machine learning model is that the driving method performed when stopping at a given station is the driving performed at either the station where the error driving occurred or the station where the error driving did not occur. It shows how similar it is to the way it is done. Therefore, it is possible to predict the possibility of error operation occurring at a given station from the predicted value output by the machine learning model.

なお、本発明の適用可能な実施形態は上述の実施形態に限定されることなく、本発明の趣旨を逸脱しない範囲で適宜変更可能なのは勿論である。 It should be noted that the applicable embodiments of the present invention are not limited to the above-described embodiments, and of course, they can be appropriately changed without departing from the spirit of the present invention.

1A,1B…判定装置
200…処理部
202…教師データ取得部
204…予測モデル生成部
206…駅データ取得部
208…予測部
210…蓋然性判定部
300A,300B…記憶部
302…判定プログラム
310A,310B…教師データ
320A,320B…予測モデルデータ
330…駅データ
340A,340B…判定結果データ
10A,10B…予測モデル
1A, 1B ... Judgment device 200 ... Processing unit 202 ... Teacher data acquisition unit 204 ... Prediction model generation unit 206 ... Station data acquisition unit 208 ... Prediction unit 210 ... Probability judgment unit 300A, 300B ... Storage unit 302 ... Judgment program 310A, 310B ... Teacher data 320A, 320B ... Prediction model data 330 ... Station data 340A, 340B ... Judgment result data 10A, 10B ... Prediction model

Claims (7)

停止すべき定位置を超過したエラー運転が所定期間の間に発生していない第1分類駅と、発生した第2分類駅とについて、当該所定期間の間に当該駅に停止する際に行われた複数の運転情報の統計データである駅別運転データを取得することと、
前記駅別運転データと、当該駅別運転データが前記第1分類駅の情報か前記第2分類駅の情報かを示す駅識別情報とを教師データとして、所与の駅の駅別運転データが入力されたときに、当該所与の駅が前記第1分類駅又は前記第2分類駅に該当する可能性を示す予測値を出力する機械学習モデルを生成することと、
を含む学習モデル生成方法。
It is performed when the first-class station and the second-class station where error operation exceeding the fixed position to be stopped did not occur during the predetermined period and stopped at the station during the predetermined period. Acquiring station-specific driving data, which is statistical data of multiple driving information,
The station-specific operation data of a given station is obtained by using the station-specific operation data and the station identification information indicating whether the station-specific operation data is the information of the first classification station or the information of the second classification station as teacher data. To generate a machine learning model that outputs a predicted value indicating the possibility that the given station corresponds to the first classification station or the second classification station when input.
Learning model generation method including.
前記駅別運転データには、当該駅の定位置に停止するまでの速度パターンの種類数が含まれる、
請求項1に記載の学習モデル生成方法。
The station-specific operation data includes the number of types of speed patterns until the station stops at a fixed position.
The learning model generation method according to claim 1.
前記駅別運転データには、常用最大ブレーキの使用割合が含まれ、
前記機械学習モデルは、ブレーキ遅れを主因とする前記エラー運転に係る前記予測値を出力するためのモデルである、
請求項1又は2に記載の学習モデル生成方法。
The operation data for each station includes the usage rate of the maximum regular brake.
The machine learning model is a model for outputting the predicted value related to the error driving mainly due to the brake delay.
The learning model generation method according to claim 1 or 2.
前記駅別運転データには、前記定位置から所定距離手前の位置での走行速度のバラツキに関する情報が含まれる、
請求項3に記載の学習モデル生成方法。
The station-specific driving data includes information on variations in traveling speed at a position a predetermined distance before the fixed position.
The learning model generation method according to claim 3.
前記取得することは、当該駅の定位置が編成車両数に応じて異なるか否かに関する定位置可変情報を更に取得し、
前記生成することは、前記教師データに前記定位置可変情報を更に含めて前記機械学習モデルを生成する、
請求項1又は2に記載の学習モデル生成方法。
The acquisition further acquires the fixed position variable information regarding whether or not the fixed position of the station differs depending on the number of trains.
The generation generates the machine learning model by further including the fixed position variable information in the teacher data.
The learning model generation method according to claim 1 or 2.
請求項1〜5の何れか一項に記載の学習モデル生成方法によって生成された機械学習モデルに対して、所与の駅に係る駅別運転データを入力することで得られる前記予測値に基づいて、当該駅で前記エラー運転が発生する蓋然性を判定する、
蓋然性判定方法。
Based on the predicted value obtained by inputting station-specific operation data for a given station to the machine learning model generated by the learning model generation method according to any one of claims 1 to 5. To determine the probability that the error driving will occur at the station.
Probability determination method.
停止すべき定位置を超過したエラー運転が所定期間の間に発生していない第1分類駅と、発生した第2分類駅とについて、当該所定期間の間に当該駅に停止する際に行われた複数の運転情報の統計データである駅別運転データと、当該駅別運転データが前記第1分類駅の情報か前記第2分類駅の情報かを示す駅識別情報とを教師データとして、所与の駅に係る駅別運転データが入力されたときに、当該所与の駅が前記第1分類駅又は前記第2分類駅に該当する可能性を示す予測値を出力するように学習させた機械学習モデルに、所与の駅に係る駅別運転データを入力して、当該所与の駅の前記予測値を得る予測処理を実行する予測処理実行手段と、
前記予測処理で得られた前記予測値に基づいて、前記入力した所与の駅で前記エラー運転が発生する蓋然性を判定する蓋然性判定手段と、
を備える判定装置。
It is performed when the first category station where the error operation exceeding the fixed position to be stopped does not occur during the predetermined period and the second category station where it occurs stops at the station during the predetermined period. Station-specific driving data, which is statistical data of a plurality of driving information, and station identification information indicating whether the station-specific driving data is the information of the first classification station or the information of the second classification station are used as teacher data. When the station-specific operation data for a given station is input, the trained trainer is to output a predicted value indicating the possibility that the given station corresponds to the first classification station or the second classification station. A prediction processing execution means for inputting station-specific operation data for a given station into a machine learning model and executing a prediction process for obtaining the predicted value of the given station.
Based on the predicted value obtained in the prediction process, a probability determining means for determining the probability that the error driving will occur at the input given station, and a probability determining means.
Judgment device including.
JP2020042630A 2020-03-12 2020-03-12 Learning model generation method, probability determination method, and determination device Active JP7328921B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2020042630A JP7328921B2 (en) 2020-03-12 2020-03-12 Learning model generation method, probability determination method, and determination device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2020042630A JP7328921B2 (en) 2020-03-12 2020-03-12 Learning model generation method, probability determination method, and determination device

Publications (2)

Publication Number Publication Date
JP2021142873A true JP2021142873A (en) 2021-09-24
JP7328921B2 JP7328921B2 (en) 2023-08-17

Family

ID=77765749

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2020042630A Active JP7328921B2 (en) 2020-03-12 2020-03-12 Learning model generation method, probability determination method, and determination device

Country Status (1)

Country Link
JP (1) JP7328921B2 (en)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2013055777A (en) * 2011-09-02 2013-03-21 Mitsubishi Electric Corp Control data collection estimation device and control data collection estimation method
JP2018137857A (en) * 2017-02-21 2018-08-30 三菱重工エンジニアリング株式会社 Vehicle control device, vehicle control method, program
JP2019204170A (en) * 2018-05-21 2019-11-28 Kddi株式会社 Information processing device and information processing program

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2013055777A (en) * 2011-09-02 2013-03-21 Mitsubishi Electric Corp Control data collection estimation device and control data collection estimation method
JP2018137857A (en) * 2017-02-21 2018-08-30 三菱重工エンジニアリング株式会社 Vehicle control device, vehicle control method, program
JP2019204170A (en) * 2018-05-21 2019-11-28 Kddi株式会社 Information processing device and information processing program

Also Published As

Publication number Publication date
JP7328921B2 (en) 2023-08-17

Similar Documents

Publication Publication Date Title
Arbabzadeh et al. A data-driven approach for driving safety risk prediction using driver behavior and roadway information data
Kuang et al. A tree-structured crash surrogate measure for freeways
Ramyar et al. Identification of anomalies in lane change behavior using one-class SVM
Jafarian et al. Application of fuzzy fault tree analysis for evaluation of railway safety risks: an evaluation of root causes for passenger train derailment
Touran et al. A collision model for safety evaluation of autonomous intelligent cruise control
Wen et al. Characterizing car-following behaviors of human drivers when following automated vehicles using the real-world dataset
Bulla-Cruz et al. Event-based road safety assessment: A novel approach towards risk microsimulation in roundabouts
CN113436432A (en) Method for predicting short-term traffic risk of road section by using road side observation data
Wang et al. The impact of vehicle moving violations and freeway traffic flow on crash risk: An application of plugin development for microsimulation
Amini et al. Development of a conflict risk evaluation model to assess pedestrian safety in interaction with vehicles
Mohammadi et al. Reliable, effective, and sustainable urban railways: A model for optimal planning and asset management
CN109784586B (en) Prediction method and system for danger emergence condition of vehicle danger
Zou et al. Multivariate analysis of car-following behavior data using a coupled hidden Markov model
Pasha et al. Safety and delays at level crossings in the United States: Addressing the need for multi-objective resource allocation
Bahmankhah et al. Micro driving behaviour in different roundabout layouts: Pollutant emissions, vehicular jerk, and traffic conflicts analysis
Wang et al. Estimating Rear‐End Accident Probabilities with Different Driving Tendencies at Signalized Intersections in China
US20230054889A1 (en) Automatically tracking a trajectory of a vehicle
AlRajie Investigation of using microscopic traffic simulation tools to predict traffic conflicts between right-turning vehicles and through cyclists at signalized intersections
Katrakazas et al. A simulation study of predicting conflict-prone traffic conditions in real-time
JP2021142873A (en) Learning model generation method, probability determination method and determination device
Cao et al. Trajectory data-based severe conflict prediction for expressways under different traffic states
Kitali et al. Using crash modification factors to appraise the safety effects of pedestrian countdown signals for drivers
JP7328920B2 (en) Learning model generation method, probability determination method, and determination device
Mahmud et al. Traditional approaches to Traffic Safety Evaluation (TSE): Application challenges and future directions
CN107958274A (en) A kind of structural fire protection safety index computational methods based on big data sorting algorithm

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20220607

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20230413

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20230425

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20230516

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: 20230801

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20230804

R150 Certificate of patent or registration of utility model

Ref document number: 7328921

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150