JP7719645B2 - Congestion rate prediction device and program - Google Patents
Congestion rate prediction device and programInfo
- Publication number
- JP7719645B2 JP7719645B2 JP2021113717A JP2021113717A JP7719645B2 JP 7719645 B2 JP7719645 B2 JP 7719645B2 JP 2021113717 A JP2021113717 A JP 2021113717A JP 2021113717 A JP2021113717 A JP 2021113717A JP 7719645 B2 JP7719645 B2 JP 7719645B2
- Authority
- JP
- Japan
- Prior art keywords
- congestion
- congestion rate
- section
- vehicle
- processor
- 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
Links
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Description
本発明は、混雑率予測装置及びプログラムに関する。 The present invention relates to a congestion rate prediction device and program.
従来、鉄道などの乗り物を運用する運用者などは、乗り物の快適さなどサービスレベルの維持及び向上のため、乗り物の混雑の状態を知ることで、輸送設備及び運用などの改善策に取り組んでいる。当該運用者は、乗り物の混雑の状態を知るために混雑率(乗車率)などの指標を用いている。また、当該運用者などは、乗り物の混雑の状態を、混雑率の計測及びその統計分析などにより混雑区間及び混雑時間帯などの大まかな傾向として求めている。計測技術の発展により一部の区間や一部の乗り物について混雑率をリアルタイムに計測することが可能となってきているが、大規模な路線全体の混雑率をリアルタイムに推定する技術及び列車の車両毎など詳細な混雑率を推定する技術のため、更には未来の時刻の混雑率の変化を効率よく逐次予測するためには従来とは異なる新規な手法を用いることが有用であると考えられる。 Traditionally, operators of transportation such as railways have worked to improve transportation facilities and operations by knowing the congestion state of their vehicles in order to maintain and improve service levels, such as comfort on board. To understand the congestion state of their vehicles, these operators use indicators such as congestion rate (occupancy rate). Furthermore, these operators determine the congestion state of vehicles as general trends, such as congested sections and congested time periods, by measuring congestion rates and conducting statistical analysis of the measured values. While advances in measurement technology have made it possible to measure congestion rates in real time for some sections or some vehicles, it is believed that new, unconventional methods would be useful for estimating the congestion rate of entire large-scale routes in real time, estimating detailed congestion rates for individual train cars, and efficiently predicting changes in congestion rates at future times.
本発明の実施形態が解決しようとする課題は、乗り物の混雑の状態を新規な方法で予測することができる混雑率予測装置及びプログラムを提供することである。 The problem that embodiments of the present invention aim to solve is to provide a congestion rate prediction device and program that can predict vehicle congestion conditions using a novel method.
実施形態の混雑率予測装置は、処理部を備える。処理部は、1又は複数の区間により構成される第1の区間での乗り物の混雑の状態の統計データと、前記第1の区間での前記乗り物の混雑の状態の第1の計測値を用いて、前記第1の区間での予測対象の乗り物の混雑の状態を予測する。 The congestion rate prediction device of the embodiment includes a processing unit. The processing unit predicts the congestion state of a vehicle to be predicted in a first section, which is composed of one or more sections, using statistical data on the congestion state of vehicles in the first section and a first measurement value of the congestion state of vehicles in the first section.
本発明は、乗り物の混雑の状態を新規な方法で予測することができる。 The present invention provides a novel way to predict vehicle congestion.
以下、いくつかの実施形態に係る混雑予測システムについて図面を用いて説明する。なお、以下の実施形態の説明に用いる各図面は、各部の縮尺を適宜変更している場合がある。また、以下の実施形態の説明に用いる各図面は、説明のため、構成を省略して示している場合がある。また、各図面及び本明細書中において、同一の符号は同様の要素を示す。 The following describes congestion prediction systems according to several embodiments using drawings. Note that the scale of each part in the drawings used to explain the following embodiments may be changed as appropriate. Also, for the sake of explanation, each drawing used to explain the following embodiments may omit components. Also, in each drawing and throughout this specification, the same reference numerals indicate similar elements.
〔第1実施形態〕
図1は、第1実施形態に係る混雑予測システム1及び混雑予測システム1に含まれる構成要素の要部構成の一例を示すブロック図である。混雑予測システム1は、鉄道の各列車の混雑率を求めるシステムである。なお、混雑率は、例えば、列車又は車両の定員数に対する乗客数の割合である。ただし、混雑率は、他の定義であっても良い。混雑予測システム1は、一例として、混雑予測サーバー100、運行管理システム200、計測装置300及び混雑率管理システム400を含む。
First Embodiment
FIG. 1 is a block diagram showing an example of the main configuration of a congestion prediction system 1 according to the first embodiment and the components included in the congestion prediction system 1. The congestion prediction system 1 is a system that calculates the congestion rate of each train on a railway. The congestion rate is, for example, the ratio of the number of passengers to the capacity of the train or vehicle. However, the congestion rate may be defined in other ways. The congestion prediction system 1 includes, as an example, a congestion prediction server 100, an operation management system 200, a measurement device 300, and a congestion rate management system 400.
混雑予測サーバー100、運行管理システム200、計測装置300及び混雑率管理システム400は、ネットワークNWに接続する。ネットワークNWは、典型的にはインターネットを含む通信網である。ネットワークNWは、典型的にはWAN(wide area network)を含む通信網である。ネットワークNWは、イントラネットなどのプライベートネットワークを含む通信網であっても良い。ネットワークNWは、LAN(local area network)を含む通信網であっても良い。また、ネットワークNWは、無線回線でも良いし有線回線でも良く、無線回線と有線回線とが混在していても良い。また、ネットワークNWは、専用線又は公衆携帯電話網などを含む通信網であっても良い。 The congestion prediction server 100, the traffic management system 200, the measurement device 300, and the congestion rate management system 400 are connected to a network NW. The network NW is typically a communication network that includes the Internet. The network NW is typically a communication network that includes a WAN (wide area network). The network NW may be a communication network that includes a private network such as an intranet. The network NW may be a communication network that includes a LAN (local area network). The network NW may also be a wireless line or a wired line, or may be a combination of wireless and wired lines. The network NW may also be a communication network that includes a dedicated line or a public mobile phone network.
混雑予測サーバー100は、鉄道の各列車の混雑率を求める装置である。混雑予測サーバー100は、1つの装置であっても良いし、複数の装置からなる装置であっても良い。混雑予測サーバー100は、一例として、プロセッサー110、ROM(read-only memory)120、RAM(random-access memory)130、補助記憶装置140及び通信I/F(interface)150を含む。そして、バス160などが、これら各部を接続する。なお、混雑予測サーバー100は、混雑率予測装置の一例である。 The congestion prediction server 100 is a device that calculates the congestion rate of each railway train. The congestion prediction server 100 may be a single device or a device consisting of multiple devices. As an example, the congestion prediction server 100 includes a processor 110, a ROM (read-only memory) 120, a RAM (random-access memory) 130, an auxiliary storage device 140, and a communication I/F (interface) 150. A bus 160 and the like connect these components. The congestion prediction server 100 is an example of a congestion rate prediction device.
プロセッサー110は、混雑予測サーバー100の動作に必要な演算及び制御などの処理を行うコンピューターの中枢部分であり、各種演算及び処理などを行う。プロセッサー110は、例えば、CPU(central processing unit)、MPU(micro processing unit)、SoC(system on a chip)、DSP(digital signal processor)、GPU(graphics processing unit)、ASIC(application specific integrated circuit)、PLD(programmable logic device)又はFPGA(field-programmable gate array)などである。あるいは、プロセッサー110は、これらのうちの複数を組み合わせたものである。また、プロセッサー110は、これらにハードウェアアクセラレーターなどを組み合わせたものあっても良い。プロセッサー110は、ROM120又は補助記憶装置140などに記憶されたファームウェア、システムソフトウェア及びアプリケーションソフトウェアなどのプログラムに基づいて、混雑予測サーバー100の各種の機能を実現するべく各部を制御する。また、プロセッサー110は、当該プログラムに基づいて後述する処理を実行する。なお、当該プログラムの一部又は全部は、プロセッサー110の回路内に組み込まれていても良い。また、プロセッサー110は、処理部の一例である。 The processor 110 is the central part of the computer, performing calculations, control, and other processes necessary for the operation of the congestion prediction server 100, and performs various calculations and processes. The processor 110 may be, for example, a central processing unit (CPU), a microprocessing unit (MPU), a system on a chip (SoC), a digital signal processor (DSP), a graphics processing unit (GPU), an application-specific integrated circuit (ASIC), a programmable logic device (PLD), or a field-programmable gate array (FPGA). Alternatively, the processor 110 may be a combination of several of these. The processor 110 may also be a combination of these with a hardware accelerator. The processor 110 controls each component to realize the various functions of the congestion prediction server 100 based on programs such as firmware, system software, and application software stored in the ROM 120 or the auxiliary storage device 140. The processor 110 also executes the processes described below based on these programs. Note that some or all of these programs may be incorporated into the circuitry of the processor 110. Additionally, the processor 110 is an example of a processing unit.
ROM120及びRAM130は、プロセッサー110を中枢としたコンピューターの主記憶装置である。
ROM120は、専らデータの読み出しに用いられる不揮発性メモリである。ROM120は、上記のプログラムのうち、例えばファームウェアなどを記憶する。また、ROM120は、プロセッサー110が各種の処理を行う上で使用するデータなども記憶する。
RAM130は、データの読み書きに用いられるメモリである。RAM130は、プロセッサー110が各種の処理を行う上で一時的に使用するデータを記憶するワークエリアなどとして利用される。RAM130は、典型的には揮発性メモリである。
The ROM 120 and RAM 130 are main memory devices of the computer with the processor 110 at its core.
The ROM 120 is a non-volatile memory used exclusively for reading data. The ROM 120 stores, for example, firmware among the above programs. The ROM 120 also stores data used by the processor 110 when it performs various processes.
The RAM 130 is a memory used for reading and writing data. The RAM 130 is used as a work area for storing data that is temporarily used when the processor 110 performs various processes. The RAM 130 is typically a volatile memory.
補助記憶装置140は、プロセッサー110を中枢としたコンピューターの補助記憶装置である。補助記憶装置140は、例えばEEPROM(electric erasable programmable read-only memory)、HDD(hard disk drive)又はフラッシュメモリなどである。補助記憶装置140は、上記のプログラムのうち、例えば、システムソフトウェア及びアプリケーションソフトウェアなどを記憶する。また、補助記憶装置140は、プロセッサー110が各種の処理を行う上で使用するデータ、プロセッサー110での処理によって生成されたデータ及び各種の設定値などを記憶する。 The auxiliary storage device 140 is an auxiliary storage device for a computer with the processor 110 at its core. The auxiliary storage device 140 is, for example, an EEPROM (electric erasable programmable read-only memory), a HDD (hard disk drive), or flash memory. The auxiliary storage device 140 stores the above programs, such as system software and application software. The auxiliary storage device 140 also stores data used by the processor 110 when performing various processes, data generated by the processes in the processor 110, various setting values, etc.
また、補助記憶装置140は、混雑データベース141も記憶する。混雑データベース141は、後述の分析処理によって求められた混雑率などのデータを管理及び記憶するデータベースである。混雑予測サーバー100は、混雑率の予測を行う際に混雑データベース141を参照する。 The auxiliary storage device 140 also stores a congestion database 141. The congestion database 141 is a database that manages and stores data such as congestion rates determined by the analysis process described below. The congestion prediction server 100 references the congestion database 141 when predicting congestion rates.
通信I/F150は、混雑予測サーバー100がネットワークNWなどを介して通信するためのインターフェースである。 The communication I/F 150 is an interface that allows the congestion prediction server 100 to communicate via a network NW, etc.
バス160は、コントロールバス、アドレスバス及びデータバスなどを含み、混雑予測サーバー100の各部で授受される信号を伝送する。 Bus 160 includes a control bus, address bus, and data bus, and transmits signals exchanged between each part of the congestion prediction server 100.
運行管理システム200は、鉄道の運行管理などを行うシステムである。また、運行管理システム200は、鉄道の運行に関する各種情報を記憶する。運行管理システム200は、例えば、サーバーなどを含む。 The traffic management system 200 is a system that manages train operations. The traffic management system 200 also stores various information related to train operations. The traffic management system 200 includes, for example, a server.
計測装置300は、鉄道列車の各車両の混雑率を、センサーなどを用いて計測する装置である。当該センサーは、例えば、応荷重装置又はカメラなどである。応荷重装置は、例えば各車両に搭載される。応荷重装置は、車両の搭載荷重を計測するセンサーである。計測装置300は、当該搭載荷重から、所定の計算式などを用いて混雑率を求める。また、カメラは、例えば、軌条の外側に設置され、車両を側面外側などから撮影する。あるいは、カメラは、車両の内部などに設置され、車両の内部を撮影する。計測装置300は、カメラによって撮影された画像から乗客がどの程度乗っているかを画像認識することで、混雑率を計測する。計測装置300によって計測される混雑率は、例えば、リアルタイムの計測値である。計測装置300は、計測した混雑率を混雑率管理システム400などに送信する。混雑率管理システム400は、受信した混雑率を列車別、車両別、駅区間別に記憶する。また、混雑率管理システム400は、混雑率の計測時点の駅又は駅区間の状況を示す状況データも関連付けて記憶する。当該状況は、例えば、日時、曜日、天気、天気予報、ダイヤ乱れなどの運行情報並びに近隣でのイベントの開催有無及び当該イベントの規模などである。また、混雑率管理システム400は、混雑率の対象の列車の始発駅及び行き先、普通又は急行などの種別、車両の種別及び編成数、並びにその他の列車情報も関連付けて記憶する。なお、混雑率管理システム400は、運行情報及び列車情報を例えば運行管理システム200から取得する。 The measuring device 300 is a device that measures the congestion rate of each car of a railway train using a sensor or the like. The sensor is, for example, a load-adjusting device or a camera. The load-adjusting device is mounted, for example, on each car. The load-adjusting device is a sensor that measures the load on the car. The measuring device 300 calculates the congestion rate from the load using a predetermined calculation formula or the like. The camera is, for example, installed on the outside of the rail and takes pictures of the car from the outside of the side. Alternatively, the camera is installed inside the car and takes pictures of the interior of the car. The measuring device 300 measures the congestion rate by image recognition of the number of passengers on board from the images captured by the camera. The congestion rate measured by the measuring device 300 is, for example, a real-time measurement value. The measuring device 300 transmits the measured congestion rate to the congestion rate management system 400 or the like. The congestion rate management system 400 stores the received congestion rates by train, car, and station section. The congestion rate management system 400 also associates and stores situation data indicating the status of the station or station section at the time the congestion rate is measured. This status includes, for example, operation information such as the date, time, day of the week, weather, weather forecast, and schedule disruptions, as well as whether or not an event is being held nearby and the scale of that event. The congestion rate management system 400 also associates and stores the starting station and destination of the train for which the congestion rate is being measured, its type (such as local or express), the type and number of carriages, and other train information. The congestion rate management system 400 obtains operation information and train information from, for example, the traffic management system 200.
なお、混雑率は人によって目視で計測されたものなどであっても良い。この場合、担当者などが、人によって計測された混雑率を混雑率管理システム400に入力する。混雑率管理システム400は、当該混雑率を記憶する。 The congestion rate may also be measured visually by a person. In this case, a person in charge or the like inputs the manually measured congestion rate into the congestion rate management system 400. The congestion rate management system 400 stores the congestion rate.
混雑率管理システム400は、前述のように計測装置300などによって計測された混雑率の計測値を記憶及び管理する。混雑率管理システム400は、例えば、サーバーなどを含む。 The congestion rate management system 400 stores and manages the congestion rate measurements taken by the measuring device 300, etc., as described above. The congestion rate management system 400 includes, for example, a server.
以下、第1実施形態に係る混雑予測システム1の動作を図2及び図3などに基づいて説明する。なお、動作説明における処理の内容は一例であって、同様な結果を得ることが可能な様々な処理を適宜に利用できる。図2及び図3は、混雑予測サーバー100のプロセッサー110による処理の一例を示すフローチャートである。プロセッサー110は、例えば、ROM120又は補助記憶装置140などに記憶されたプログラムに基づいて図2及び図3の処理を実行する。
プロセッサー110は、例えば、図2及び図3の処理を、混雑予測サーバー100の起動又は混雑予測用のアプリケーションの起動などに基づいて開始する。
The operation of the congestion prediction system 1 according to the first embodiment will be described below with reference to Figures 2 and 3. Note that the processing content in the operational description is an example, and various processing that can achieve similar results can be used as appropriate. Figures 2 and 3 are flowcharts showing an example of processing by the processor 110 of the congestion prediction server 100. The processor 110 executes the processing of Figures 2 and 3 based on a program stored in, for example, the ROM 120 or the auxiliary storage device 140.
The processor 110 starts the processing of Figures 2 and 3, for example, based on the activation of the congestion prediction server 100 or the activation of an application for congestion prediction.
図2に示す処理は、混雑率の予測に用いるための各種データを混雑率などの過去の混雑率データを用いて求める処理(以下「分析処理」という。)を含む処理である。
図2のステップST11において混雑予測サーバー100のプロセッサー110は、分析処理を開始するか否かを判定する。なお、分析処理は、ステップST12~ステップST17の処理を示す。プロセッサー110は、例えば、所定の条件を満たした場合に分析処理を行うと判定する。所定の条件は、前回分析処理を行ってから所定の期間経過したこと、及び混雑率の予測値と混雑率の計測値との差が大きい状態が高い頻度で起きることなどである。混雑率の予測値と混雑率の計測値との差が大きい状態が高い頻度で起きている状態とは、例えば、当該差の二乗若しくは絶対値の平均値が所定より大きい状態、当該差の二乗若しくは絶対値が所定より大きい頻度が所定以上である状態、又はこれらと同値の状態などである。所定の条件は、例えば、分析処理を行うように指示する入力があったことである。当該入力は、例えば、混雑予測サーバー100のコンソールなどを用いた操作入力に基づく。あるいは、当該入力は、他の装置などから送信されて通信I/F150へ入力されるコマンドなどに基づく。当該コマンドは、分析処理を行うように指示することを示す。
The process shown in FIG. 2 includes a process (hereinafter referred to as "analysis process") for obtaining various data to be used for predicting the congestion rate using past congestion rate data such as the congestion rate.
In step ST11 of FIG. 2 , the processor 110 of the congestion prediction server 100 determines whether to start analysis processing. The analysis processing refers to the processing of steps ST12 to ST17. The processor 110 determines to perform analysis processing when, for example, a predetermined condition is met. The predetermined condition may be that a predetermined period of time has elapsed since the previous analysis processing, or that a state in which the difference between the predicted congestion rate and the measured congestion rate is large occurs frequently. A state in which a state in which the difference between the predicted congestion rate and the measured congestion rate is large frequently may be, for example, a state in which the average value of the square or absolute value of the difference is greater than a predetermined value, a state in which the frequency at which the square or absolute value of the difference is greater than a predetermined value is greater than or equal to a predetermined value, or a state equivalent to these. The predetermined condition may be, for example, that an input has been made to instruct analysis processing to be performed. The input may be based, for example, on an operation input using a console of the congestion prediction server 100. Alternatively, the input may be based on a command transmitted from another device and input to the communication I/F 150. The command indicates an instruction to perform analysis processing.
ステップST12においてプロセッサー110は、分析処理に用いるデータ(以下「分析用データ」という。)を運行管理システム200及び混雑率管理システム400などから取得する。分析用データは、例えば、路線データ、時刻表、各列車及び各車両の各駅における混雑率の過去の所定の日時から現在までのデータなどを含む。プロセッサー110は、例えば、路線データ及び時刻表などのデータを運行管理システム200から取得し、混雑率のデータを混雑率管理システム400から取得する。
プロセッサー110は、分析用データを送信するように指示する要求を運行管理システム200及び混雑率管理システム400に送信するように通信I/F150に対して指示する。この送信の指示を受けて通信I/F150は、当該要求を運行管理システム200及び混雑率管理システム400に送信する。送信された当該要求は、運行管理システム200及び混雑率管理システム400によって受信される。そして、運行管理システム200及び混雑率管理システム400は、当該要求を受信したことに応じて分析用データを混雑予測サーバー100に送信する。当該データは、混雑予測サーバー100の通信I/F150によって受信される。
In step ST12, the processor 110 acquires data to be used for the analysis process (hereinafter referred to as "analysis data") from the traffic control system 200, the congestion rate management system 400, etc. The analysis data includes, for example, route data, timetables, and congestion rate data for each train and each vehicle at each station from a predetermined past date and time to the present. The processor 110 acquires, for example, data such as route data and timetables from the traffic control system 200, and congestion rate data from the congestion rate management system 400.
The processor 110 instructs the communication I/F 150 to send a request to the traffic management system 200 and the congestion rate management system 400 to send analysis data. Upon receiving this transmission instruction, the communication I/F 150 sends the request to the traffic management system 200 and the congestion rate management system 400. The sent request is received by the traffic management system 200 and the congestion rate management system 400. Then, in response to receiving the request, the traffic management system 200 and the congestion rate management system 400 send analysis data to the congestion prediction server 100. The data is received by the communication I/F 150 of the congestion prediction server 100.
ステップST13においてプロセッサー110は、分析用データを用いて、各駅区間の列車平均混雑率を求める(算出する)。なお、駅区間の列車平均混雑率とは、駅を発車してから次の駅に着くまでの列車の混雑率の平均である。例えば、A駅-B駅区間の列車平均混雑率は、A駅を発車してからB駅に着くまでの列車の混雑率の平均である。A駅-B駅区間の混雑率は、例えば、列車がA駅を発車又は通過する時点の混雑率でも良いしB駅に到着又は通過する時点での混雑率でも良いしA駅とB駅の間の位置に到着又は通過する時点での混雑率でも良い。なお、B駅は、A駅の次の駅であるとする。また、列車平均混雑率は、所定の期間における列車の混雑率の平均である。したがって、複数車両からなる列車の列車平均混雑率は、当該列車の各車両の混雑率を平均した値と同じである。プロセッサー110は、列車平均混雑率を、例えば、出発時刻ごとに求める。出発時刻は、ダイヤ上の時刻であっても良いし、実際に計測した時刻であっても良い。一例として、プロセッサー110は、8時10分発の列車の列車平均混雑率と、8時14分発の列車の列車平均混雑率を別々に求める。あるいは、プロセッサー110は、出発時刻の時間帯ごとに列車平均混雑率を求めても良い。一例として、プロセッサー110は、8時00分から8時15分までの時間帯に発車する列車の列車平均混雑率と、8時15分から8時30分までの時間帯に発車する列車の列車平均混雑率を求める。また、プロセッサー110は、出発時刻以外にも様々な条件ごとに列車平均混雑率を求める。なお、列車平均混雑率などの条件分けに用いる条件を以下「混雑条件」という。プロセッサー110は、混雑条件として以下に挙げた条件の一部又は全部を用いても良く、また、以下に挙げていない条件を用いても良い。当該条件は、例えば以下のようなものを挙げることができる。なお、プロセッサー110は、混雑条件による条件分けをせずに列車平均混雑率を求めても良い。
・曜日
・平日、土曜日及び日祝のいずれであるかなど。
・列車情報
・運行情報
・天気
・天気予報
・連休であるか否か。
・季節又は時期
・月日などの暦
・駅の近隣でイベントが開催されているか否か。またそのイベントの規模。
In step ST13, the processor 110 uses the analysis data to determine (calculate) the average train congestion rate for each station section. The average train congestion rate for a station section is the average train congestion rate from the time a train departs from a station until it arrives at the next station. For example, the average train congestion rate for the Station A-Station B section is the average train congestion rate from the time a train departs from Station A until it arrives at Station B. The congestion rate for the Station A-Station B section may be, for example, the congestion rate at the time a train departs from or passes through Station A, the congestion rate at the time it arrives at or passes through Station B, or the congestion rate at the time it arrives at or passes a location between Station A and Station B. It is assumed that Station B is the station next to Station A. The average train congestion rate is also the average train congestion rate for a specified period. Therefore, the average train congestion rate for a train consisting of multiple cars is the same as the average congestion rate of each car of the train. The processor 110 determines the average train congestion rate, for example, for each departure time. The departure time may be a time listed on the timetable or an actually measured time. As an example, the processor 110 separately calculates the average train congestion rate for the train departing at 8:10 and the average train congestion rate for the train departing at 8:14. Alternatively, the processor 110 may calculate the average train congestion rate for each time slot of the departure time. As an example, the processor 110 calculates the average train congestion rate for trains departing between 8:00 and 8:15 and the average train congestion rate for trains departing between 8:15 and 8:30. The processor 110 also calculates the average train congestion rate for various conditions other than the departure time. The conditions used to classify the average train congestion rate and other conditions are hereinafter referred to as "congestion conditions." The processor 110 may use some or all of the conditions listed below as congestion conditions, or may use conditions not listed below. Examples of such conditions include the following. The processor 110 may also calculate the average train congestion rate without classifying the conditions based on congestion conditions.
- Day of the week - whether it is a weekday, Saturday, Sunday or public holiday, etc.
-Train information, operation information, weather, weather forecast, whether it is a holiday or not.
- Season or time of year, date, etc. - Whether or not an event is being held near the station, and the scale of the event.
一例として、プロセッサー110は、平日の8時14分発の列車の列車平均混雑率と、日曜日の8時14分発の列車の列車平均混雑率を別々に求める。別の一例として、プロセッサー110は、晴れの日の春の平日の8時14分発の列車の列車平均混雑率と雨の日の春の平日の8時14分発の列車の列車平均混雑率と晴れの日の夏の平日の8時14分発の列車の列車平均混雑率とを別々に求める。 As one example, processor 110 separately calculates the average train congestion rate for trains departing at 8:14 on weekdays and the average train congestion rate for trains departing at 8:14 on Sundays. As another example, processor 110 separately calculates the average train congestion rate for trains departing at 8:14 on sunny spring weekdays, the average train congestion rate for trains departing at 8:14 on rainy spring weekdays, and the average train congestion rate for trains departing at 8:14 on sunny summer weekdays.
また、プロセッサー110は、求めた列車平均混雑率を、当該列車平均混雑率の条件分けに用いた混雑条件を示す情報と関連付けて、混雑データベース141に記憶する。 In addition, the processor 110 associates the calculated average train congestion rate with information indicating the congestion conditions used to classify the average train congestion rate, and stores it in the congestion database 141.
ステップST14においてプロセッサー110は、分析用データを用いて号車別混雑分布を求める。号車別混雑分布は、列車についての車両(号車)ごとの混雑率を示す。プロセッサー110は、ステップST13で求める列車平均混雑率と同様に、条件分けをして号車別混雑分布を求めても良いし、条件分けせずに号車別分布を求めても良い。また、プロセッサー110は、ステップST13とは異なる混雑条件で号車別分布を求めても良い。号車別混雑分布は、例えば、各車両について、車両の平均の混雑率が列車平均混雑率の何倍であるかを示す。したがって、号車別混雑分布の値が1である車両の混雑率は、混雑条件が同一である列車平均混雑率と同じ値となる。 In step ST14, processor 110 uses the analysis data to determine the congestion distribution by car. The congestion distribution by car indicates the congestion rate for each car (car number) of the train. Processor 110 may determine the congestion distribution by car by dividing the conditions, as with the average train congestion rate determined in step ST13, or may determine the distribution by car without dividing the conditions. Processor 110 may also determine the distribution by car under congestion conditions different from those in step ST13. The congestion distribution by car indicates, for each car, for example, how many times the average congestion rate of the car is compared to the average train congestion rate. Therefore, the congestion rate of a car for which the value of the congestion distribution by car is 1 will be the same as the average train congestion rate for the same congestion conditions.
また、プロセッサー110は、求めた号車別混雑分布を、当該号車別混雑分布の条件分けに用いた混雑条件を示す情報と関連付けて、混雑データベース141に記憶する。 The processor 110 also associates the calculated congestion distribution by car with information indicating the congestion conditions used to classify the congestion distribution by car, and stores it in the congestion database 141.
なお、プロセッサー110は、列車平均混雑率又は号車別混雑分布のうち、いずれかについて、今回の処理では求めずに、前回求めたものを用いても良い。この場合、プロセッサー110は、前回求めた列車平均混雑率又は号車別混雑分布を混雑データベース141から取得する。 In addition, the processor 110 may use the previously calculated average train congestion rate or the previously calculated congestion distribution by car, rather than calculating it in the current process. In this case, the processor 110 obtains the previously calculated average train congestion rate or previously calculated congestion distribution by car from the congestion database 141.
ステップST15においてプロセッサー110は、号車別混雑率を求める。号車別混雑率は、各車両の混雑率である。プロセッサー110は、例えば、列車平均混雑率に各車両の号車別混雑分布をかける。これにより、プロセッサー110は、混雑条件ごとの各車両の混雑率を求めることができる。号車別混雑率は、列車平均混雑率に各車両の号車別混雑分布をかけたものであるので、列車平均混雑率と同様に条件分けされたものである。 In step ST15, processor 110 calculates the congestion rate by car. The congestion rate by car is the congestion rate of each car. Processor 110, for example, multiplies the average train congestion rate by the congestion distribution by car of each car. This allows processor 110 to calculate the congestion rate of each car for each congestion condition. The congestion rate by car is calculated by multiplying the average train congestion rate by the congestion distribution by car of each car, and is therefore categorized by conditions in the same way as the average train congestion rate.
図4は、号車別混雑分布をグラフ化した図である。図4は、特定の混雑条件の列車についての号車別混雑率の一例を示している。当該列車は、1号車から8号車までの8両の車両を含むものとする。図4は、当該列車の1号車から8号車までのそれぞれの号車別混雑率を棒グラフで示している。また、折れ線Laは、当該棒グラフを折れ線グラフとして示したものである。図4は、当該列車の列車平均混雑率を示す混雑率Caも示している。図4に示す号車別混雑率では、1号車2号車、4号車及び5号車が列車平均混雑率よりも混雑率が低い。すなわち、これらの車両の号車別混雑分布の値は1未満である。そして、3号車、6号車、7号車及び8号車は、列車平均混雑率よりも混雑率が高い。すなわち、これらの車両の号車別混雑分布の値は1超である。 Figure 4 is a graph of the congestion distribution by car. Figure 4 shows an example of the congestion rate by car for a train under specific congestion conditions. The train includes eight cars, from car 1 to car 8. Figure 4 shows the congestion rate by car for each of cars 1 to 8 of the train as a bar graph. Line La shows the bar graph as a line graph. Figure 4 also shows the congestion rate Ca, which indicates the average congestion rate of the train. In the congestion rate by car shown in Figure 4, cars 1, 2, 4, and 5 have lower congestion rates than the train average congestion rate. In other words, the congestion distribution values by car for these cars are less than 1. Furthermore, cars 3, 6, 7, and 8 have higher congestion rates than the train average congestion rate. In other words, the congestion distribution values by car for these cars are greater than 1.
ステップST16においてプロセッサー110は、変動表を作成する。変動表は、例えば、列車の各車両についての、駅区間ごとの号車別混雑率の変動を示す表である。プロセッサー110は、例えば、号車別混雑率を用いて変動表を作成する。変動表は、号車別混雑率と同様に条件分けされたものであっても条件分けされていないものでも良い。プロセッサー110は、ある条件の変動表を作成する場合、当該ある条件に該当する号車別混雑率を用いて変動表を作成する。プロセッサー110は、一例としてデータクラスタリングにより変動表を作成する。ある条件の変動表中の特定の駅区間の特定の号車の混雑率は、条件、駅区間及び号車が同じ号車別混雑率の平均である。例えば、条件が晴れの日である変動表の9時00分のB駅-C駅区間の3号車の混雑率は、晴れの日の9時00分のB駅-C駅区間の3号車の号車別混雑率の平均である。当該号車別混雑率は、晴れの日の9時00分であればよく、季節などの変動表の条件に含まれない条件は何でも良い。あるいは、プロセッサー110は、号車別混雑率に代えて、混雑率管理システム400に記憶された混雑率を用いても良い。変動表が条件分けされている場合、一例として、休日におけるA駅7時50分始発の列車の各駅区間の平均混雑率などが変動表から分かる。なお、変動表は、駅区間ごとの列車の混雑の状態の変化を示す変動情報の一例である。また、変動表は、第1の区間と第2の区間の相関関係を示すデータの一例である。 In step ST16, processor 110 creates a fluctuation table. The fluctuation table is, for example, a table showing fluctuations in congestion rates by car for each station section for each car of a train. Processor 110 creates the fluctuation table, for example, using congestion rates by car. The fluctuation table may be categorized by conditions, as with congestion rates by car, or may not be categorized by conditions. When creating a fluctuation table for a certain condition, processor 110 creates the fluctuation table using congestion rates by car that meet the certain condition. As an example, processor 110 creates the fluctuation table using data clustering. The congestion rate of a specific car in a specific station section in a fluctuation table for a certain condition is the average congestion rate by car for the same condition, station section, and car. For example, in a fluctuation table where the condition is a sunny day, the congestion rate of car 3 in the section between stations B and C at 9:00 is the average congestion rate by car for car 3 in the section between stations B and C at 9:00 on a sunny day. The congestion rate for each car may be 9:00 on a sunny day, and may be any condition not included in the conditions of the fluctuation table, such as season. Alternatively, the processor 110 may use the congestion rate stored in the congestion rate management system 400 instead of the congestion rate for each car. When the fluctuation table is categorized by conditions, for example, the average congestion rate for each station section for the first train departing from Station A at 7:50 on a holiday can be determined from the fluctuation table. The fluctuation table is an example of fluctuation information that shows changes in the train congestion state for each station section. The fluctuation table is also an example of data that shows the correlation between the first section and the second section.
図5は、変動表をグラフ化した図である。図5は、特定の混雑条件の列車の特定の車両についての変動表の一例を示している。当該列車は、A駅、B駅、C駅、D駅、E駅、F駅、G駅の順に進んでいく列車であるとする。図5は、当該特定の車両の各駅区間における号車別混雑率を棒グラフで示している。また、折れ線Lbは、当該棒グラフを折れ線グラフとして示したものである。 Figure 5 is a graph of the fluctuation table. Figure 5 shows an example of a fluctuation table for a specific car of a train under specific congestion conditions. Assume that the train in question travels through stations A, B, C, D, E, F, and G in that order. Figure 5 uses a bar graph to show the congestion rate by car for that specific car in each station section. Furthermore, line Lb is a line graph that shows the bar graph.
ステップST17においてプロセッサー110は、各駅区間について混雑条件ごとに時系列モデルを作成する。時系列モデルは、例えば駅区間の列車の混雑率の時間変化を号車ごとに示す時間の関数である。プロセッサー110は、特定の混雑条件における特定の駅区間の時系列モデルを作成する場合、当該駅区間の当該混雑条件の号車別混雑率を用いる。プロセッサー110は、当該号車別混雑率から、最小二乗法、AR(autoregressive model)、ARMA(autoregressive moving average model)、ARIMA(autoregressive integrated moving average model)又は状態空間モデルなどを使用して時系列モデルを求める。あるいは、プロセッサー110は、各号車別混雑率を時系列順に直線で結んだものを時系列モデルとしても良い。なお、時系列モデルは、駅区間での列車の混雑の状態の統計データの一例である。 In step ST17, processor 110 creates a time series model for each congestion condition for each station section. The time series model is, for example, a function of time that indicates the change in train congestion rate for each car in the station section over time. When processor 110 creates a time series model for a specific station section under specific congestion conditions, it uses the congestion rate for each car under the congestion conditions in that station section. From the congestion rate for each car, processor 110 obtains the time series model using least squares, AR (autoregressive model), ARMA (autoregressive moving average model), ARIMA (autoregressive integrated moving average model), or a state space model. Alternatively, processor 110 may use a straight line connecting the congestion rates for each car in chronological order as the time series model. The time series model is an example of statistical data on the train congestion state in the station section.
ステップST18においてプロセッサー110は、分析処理で求めた号車別混雑率、変動表及び時系列モデルなどの各データを混雑データベース141などに記憶する。プロセッサー110は、ステップST18の処理の後、ステップST11へと戻る。 In step ST18, the processor 110 stores the data obtained by the analysis process, such as the congestion rate by car, fluctuation table, and time series model, in the congestion database 141, etc. After processing step ST18, the processor 110 returns to step ST11.
図3に示す処理は、号車別混雑率及び変動表などを用いて混雑率を予測する処理を含む処理である。図6を用いて図3の処理について説明する。
図6は、A駅からG駅までの各駅区間における複数の列車の混雑率を示した図である。当該混雑率は、例えば、計測装置300によって計測される。混雑率A1~混雑率F6までの各点の高さが混雑率を示す。ただし、図6に混雑率A1~混雑率F6の符号を全て図示すると図が見づらくなるため、一部の符号は省略して図示していない。なお、混雑率を示す各点の符号は、左側のアルファベットと右側の数字の2つの部分を含んで成り立っている。一例として「A1」であれば左側のアルファベットの「A」と、右側の数字の「1」の2つの部分を含む。左側のアルファベットは、駅区間のうち、列車が先に到着する側の駅を示す。例えば、AであればA駅-B駅区間、BであればB駅-C駅区間、CであればC駅-D駅区間である。また、右側の数字は、駅区間を通る列車の順番を示す。なお、駅区間を列車が通る時点は、例えば、当該駅区間内の特定の位置を列車が通る時点を示す。当該位置は、例えば、駅区間の端の駅であっても良いし間の位置であっても良い。一例として、列車がA駅-B駅区間を通る時点とは、列車がA駅を発車若しくは通過する時点、B駅に到着若しくは通過する時点、又はA駅とB駅の間の特定の位置に到着若しくは通過する時点である。なお、駅区間を列車が通る時点は、実際に計測された時刻に限らず、ダイヤ上の時刻であっても良い。例えば、混雑率A1は、A駅-B駅区間を通る1番目の列車の混雑率を示す。この場合、混雑率A2は、A駅-B駅区間を通る2番目の列車の混雑率を示す。すなわち、混雑率A1に対応する列車がA駅-B駅区間を通過した次にA駅-B駅区間を通る列車は混雑率A2に対応する列車である。また例えば、混雑率C3は、C駅-D駅区間を通る3番目の列車の混雑率を示す。なお、右側の数字が1である符号に対応する列車が特定の駅区間を通るより前に当該駅区間を通る列車があっても良い。また、右側の数字が同じ場合、同じ列車を示す。したがって、混雑率A1~混雑率F1は、全て同じ列車の混雑率を示す。ある特定の列車のA駅-B駅区間の混雑率が混雑率A1、当該特定の列車のB駅-C駅区間の混雑率が混雑率B1、当該特定の列車のC駅-D駅区間の混雑率が混雑率C1、…といった具合である。なお、図6に示す混雑率は、特定の号車の混雑率を代表的に示している。特定の号車は、何号車であっても良いが、一例として混雑率A1が3号車の混雑率を示したものであれば、混雑率A1~混雑率F6は全て3号車の混雑率を示したものである。
The process shown in Fig. 3 includes a process for predicting the congestion rate by car number using a congestion rate table, etc. The process shown in Fig. 3 will be described with reference to Fig. 6 .
Figure 6 is a diagram showing the congestion rates of multiple trains in each station section from Station A to Station G. The congestion rates are measured, for example, by a measurement device 300. The height of each point, from congestion rate A1 to congestion rate F6, indicates the congestion rate. However, because showing all the symbols for congestion rates A1 to F6 in Figure 6 would make the diagram difficult to read, some symbols have been omitted. Note that the symbols for each point indicating the congestion rate consist of two parts: an alphabet on the left and a number on the right. For example, "A1" includes two parts: the alphabet "A" on the left and the number "1" on the right. The alphabet on the left indicates the station in the station section where the train arrives first. For example, A indicates the Station A-Station B section, B indicates the Station B-Station C section, and C indicates the Station C-Station D section. The numbers on the right indicate the order in which trains pass through the station section. Note that the time at which a train passes through a station section indicates, for example, the time at which the train passes a specific location within the station section. The location may be, for example, a station at the end of the station section or a location in between. As an example, the time when a train passes through the Station A-Station B section refers to the time when the train departs from or passes through Station A, the time when the train arrives at or passes through Station B, or the time when the train arrives at or passes through a specific location between Station A and Station B. The time when a train passes through a station section is not limited to an actually measured time, but may also be a time on a timetable. For example, congestion rate A1 indicates the congestion rate of the first train passing through the Station A-Station B section. In this case, congestion rate A2 indicates the congestion rate of the second train passing through the Station A-Station B section. In other words, after the train corresponding to congestion rate A1 passes through the Station A-Station B section, the train passing through the Station A-Station B section next is the train corresponding to congestion rate A2. For example, congestion rate C3 indicates the congestion rate of the third train passing through the Station C-Station D section. Note that there may be a train passing through the station section before the train corresponding to a code with a number 1 on the right side passes through the specific station section. Furthermore, the same numbers on the right side indicate the same train. Therefore, congestion rates A1 to F1 all indicate the congestion rates of the same train. The congestion rate of a specific train in the section between station A and station B is congestion rate A1, the congestion rate of that specific train in the section between station B and station C is congestion rate B1, the congestion rate of that specific train in the section between station C and station D is congestion rate C1, and so on. The congestion rates shown in Figure 6 representatively indicate the congestion rate of a specific car. The specific car can be any car number, but as an example, if congestion rate A1 indicates the congestion rate of car 3, congestion rates A1 to F6 all indicate the congestion rate of car 3.
ここでは、混雑率C1~混雑率C4を用いて他の混雑率を予測する場合を例に図3の処理を説明する。この場合、C駅-D駅区間は、第1の区間の一例である。なお、第1の区間は、1つの駅区間で構成されていても良いし、複数の駅区間によって構成されていても良い。また、1つの駅区間の駅と駅の間に別の駅があっても良い。また、駅区間は、駅と駅の間に限らない。例えば、駅に代えて操車場、信号場又は停車場などであっても良い。 Here, the processing in Figure 3 will be explained using congestion rates C1 to C4 as an example to predict other congestion rates. In this case, the Station C-D section is an example of a first section. Note that the first section may consist of a single station section, or multiple station sections. There may also be another station between stations in a single station section. Also, station sections are not limited to sections between stations. For example, instead of stations, they may be marshalling yards, signal stations, or train stations.
図3のステップST21においてプロセッサー110は、混雑率の予測を開始するか否かを判定する。例えば、プロセッサー110は、混雑率管理システム400が計測装置300から新たに列車の最新の混雑率を取得したことに応じて、当該列車の混雑率実測値を用いた混雑率の予測を開始すると判定する。プロセッサー110は、混雑率の予測を開始すると判定しないならば、ステップST21においてNoと判定してステップST21の処理を繰り返す。対して、プロセッサー110は、混雑率の予測を開始すると判定するならば、ステップST21においてYesと判定してステップST22へと進む。なお、混雑率管理システム400が新たに取得した最新の混雑率は、ここでは混雑率C4であるとする。 In step ST21 of FIG. 3, processor 110 determines whether to start predicting the congestion rate. For example, when the congestion rate management system 400 acquires the latest congestion rate of a train from the measurement device 300, processor 110 determines to start predicting the congestion rate using the actual measured congestion rate value of the train. If processor 110 does not determine to start predicting the congestion rate, it determines No in step ST21 and repeats the processing of step ST21. On the other hand, if processor 110 determines to start predicting the congestion rate, it determines Yes in step ST21 and proceeds to step ST22. Here, the latest congestion rate newly acquired by the congestion rate management system 400 is assumed to be congestion rate C4.
ステップST22においてプロセッサー110は、混雑率の予測に使用する混雑率データを混雑率管理システム400などから取得する。プロセッサー110は、ここで、リアルタイムの混雑率などを取得することができる。一例として、プロセッサー110は、当該混雑率データとして、混雑率C1~混雑率C4を取得する。また、プロセッサー110は、各混雑率データに関連付けて記憶されている状況データも取得する。混雑率C1~混雑率C4などの、ステップST22において取得される同一駅区間の複数の混雑率データは、駅区間での列車の混雑の状態の第1の計測値の一例である。 In step ST22, processor 110 acquires congestion rate data to be used for predicting the congestion rate from the congestion rate management system 400 or the like. Here, processor 110 can acquire real-time congestion rates, etc. As an example, processor 110 acquires congestion rates C1 to C4 as the congestion rate data. Processor 110 also acquires situation data stored in association with each congestion rate data. Multiple congestion rate data for the same station section acquired in step ST22, such as congestion rates C1 to C4, are an example of a first measurement value of the train congestion state in the station section.
ステップST23においてプロセッサー110は、混雑率の予測に使用する時系列モデル及び変動表などを混雑データベース141から取得する。混雑率の予測に使用される時系列モデルは、ステップST22で取得された状況データが示す状況と一致する混雑条件の時系列モデルである。また、混雑率の予測に使用される時系列モデルは、ステップST22で取得された混雑率データの駅区間と同一の駅区間の時系列モデルである。そして、混雑率の予測に使用される変動表は、ステップST22で取得された状況データが示す状況と一致する混雑条件の変動表である。また、混雑率の予測に使用される変動表は、混雑率を予測する対象の列車についての変動表である。 In step ST23, the processor 110 acquires from the congestion database 141 a time series model and a fluctuation table to be used to predict the congestion rate. The time series model used to predict the congestion rate is a time series model of congestion conditions that match the situation indicated by the situation data acquired in step ST22. The time series model used to predict the congestion rate is a time series model of the same station section as the station section of the congestion rate data acquired in step ST22. The fluctuation table used to predict the congestion rate is a fluctuation table of congestion conditions that match the situation indicated by the situation data acquired in step ST22. The fluctuation table used to predict the congestion rate is a fluctuation table for the train for which the congestion rate is to be predicted.
ステップST24においてプロセッサー110は、ステップST22で取得した混雑率データとステップST23で取得した時系列モデルの差分を求める。図7を用いて、計測された混雑率データ、時系列モデル及び差分の関係を説明する。
図7は、混雑率の時間変化の一例を示すグラフである。図7のグラフは、横軸が時間(時刻)、縦軸が混雑率である。図7には、ステップST22で取得された混雑率データの一例として点Paをプロットしている。ただし、図7に示す点Paは、混雑率C1~混雑率C4とは別の例である。曲線Lcは、ステップST23で取得される時系列モデルの一例を示す。また、差分Daは、点Paと曲線Lcとの混雑率の差を示す。なお、差分Daは、統計データと計測値の差分の一例を示す。
In step ST24, the processor 110 calculates the difference between the congestion rate data acquired in step ST22 and the time series model acquired in step ST23. The relationship between the measured congestion rate data, the time series model, and the difference will be described with reference to FIG. 7.
FIG. 7 is a graph showing an example of a change in congestion rate over time. In the graph of FIG. 7, the horizontal axis represents time (hour) and the vertical axis represents congestion rate. In FIG. 7, point Pa is plotted as an example of congestion rate data acquired in step ST22. However, point Pa shown in FIG. 7 is a different example from congestion rates C1 to C4. Curve Lc represents an example of a time series model acquired in step ST23. Furthermore, difference Da represents the difference in congestion rate between point Pa and curve Lc. Note that difference Da represents an example of the difference between statistical data and a measured value.
ステップST25においてプロセッサー110は、ステップST24で求めた差分から、差分の予測値を求める。例えば、プロセッサー110は、最小二乗法、AR、ARMA、ARIMA又は状態空間モデルなどを使用して差分を時間の関数(以下「差分関数」という。)で示す。図8を用いて差分と差分関数の関係を説明する。
図8は、差分の時間変化の一例を示すグラフである。図8のグラフは、縦軸が混雑率の差分、横軸が時間である。点Pbは、差分Daをグラフにプロットしたものである。曲線Ldは、差分Daから求められた差分関数の一例を示す。差分関数を用いることで、最新の点Pbよりも先の時間の差分を予測することが可能である。
In step ST25, the processor 110 calculates a predicted value of the difference from the difference calculated in step ST24. For example, the processor 110 expresses the difference as a function of time (hereinafter referred to as a "difference function") using the least squares method, AR, ARMA, ARIMA, a state space model, or the like. The relationship between the difference and the difference function will be explained using FIG. 8.
Fig. 8 is a graph showing an example of how the difference changes over time. In the graph of Fig. 8, the vertical axis represents the difference in congestion rate, and the horizontal axis represents time. Point Pb is the difference Da plotted on the graph. Curve Ld shows an example of a difference function calculated from the difference Da. By using the difference function, it is possible to predict the difference for a time beyond the latest point Pb.
ステップST26においてプロセッサー110は、差分関数を用いて、図7の最新の点Paよりも先の時間の混雑率を予測する。曲線Leは、混雑率の時間変化の予測値を示す時間の関数(以下「時系列予測関数」という。)であり、時系列モデルと差分関数を足した関数である。例えば、時系列モデルを関数Fa(t)、差分関数を関数Fb(t)、時系列予測関数を関数Fc(t)とすると、Fc(t)=Fa(t)+Fb(t)である。プロセッサー110は、時系列予測関数を用いることで、任意の時間における混雑率Fc(t)を求めることができる。プロセッサー110は、例えば、駅区間を列車が通る時刻での混雑率を、時系列予測関数から求める。プロセッサー110は、一例として、時系列予測関数を用いて混雑率C5及び混雑率C6を求める。なお、プロセッサー110は、駅区間の混雑率を求める場合、当該駅区間のうち列車が先に到達する側の駅を当該列車が発車する時刻の混雑率を、時系列予測関数を用いて求める。一例として、プロセッサー110は、C駅-D駅区間の混雑率を求める場合、C駅を発車する時刻の混雑率を、時系列予測関数を用いて求める。あるいは、プロセッサー110は、駅区間の特定の位置を列車が通過又は到達する時刻の混雑率を、時系列予測関数を用いて求めても良い。なお、時系列予測関数は、駅区間での列車の混雑の状態の時間変化を予測する関数の一例である。 In step ST26, processor 110 uses a differential function to predict the congestion rate for a time beyond the latest point Pa in Figure 7. Curve Le is a function of time (hereinafter referred to as the "time series prediction function") that indicates the predicted value of the change in congestion rate over time, and is a function obtained by adding a time series model and a differential function. For example, if the time series model is function Fa(t), the differential function is function Fb(t), and the time series prediction function is function Fc(t), then Fc(t) = Fa(t) + Fb(t). Processor 110 can calculate the congestion rate Fc(t) at any time by using the time series prediction function. Processor 110, for example, calculates the congestion rate at the time a train passes through the station section from the time series prediction function. As an example, processor 110 calculates congestion rates C5 and C6 using the time series prediction function. When calculating the congestion rate for a station section, processor 110 uses a time series prediction function to calculate the congestion rate at the time the train departs from the station that the train will arrive at first in the station section. As an example, when calculating the congestion rate for the Station C-Station D section, processor 110 uses a time series prediction function to calculate the congestion rate at the time the train departs from Station C. Alternatively, processor 110 may use a time series prediction function to calculate the congestion rate at the time the train passes through or arrives at a specific position in the station section. The time series prediction function is an example of a function that predicts changes over time in the train congestion state in a station section.
ステップST27においてプロセッサー110は、対象混雑率とステップST23で取得した変動表を用いて、対象混雑率とは異なる駅区間の混雑率を推定する。対象混雑率は、ステップST27の混雑率の推定に用いる既知の混雑率である。対象混雑率は、例えば、ステップST22で取得された混雑率及びステップST26で求められた混雑率である。プロセッサー110は、1つの対象混雑率(以下「混雑率Xa」という。)と、混雑率Xaと列車及び号車が同じ変動表(以下「対象変動表」という。)を用いて、当該対象混雑率と同じ列車及び号車の、当該対象混雑率と異なる駅区間の混雑率を推定する。対象変動表のうち、混雑率Xaと同じ駅区間の混雑率を混雑率Yaとする。プロセッサー110は、混雑率Xaと混雑率Yaがどれだけ異なった値であるかを求める。一例として、プロセッサー110は、混雑率Xaと混雑率Yaがどれだけ異なった値であるかを示す値として、混雑率Xaと混雑率Yaの差αを求める。例えば、差αは、下式(1)で求めることができる。
α=Xa-Ya (1)
ここで、推定対象の混雑率を混雑率Xbとする。また、対象変動表のうち、混雑率Xbと同じ駅区間の混雑率を混雑率Ybとする。この場合、プロセッサー110は、例えば、下式(2)を用いて混雑率Xbを求める。
Xb=Yb+α (2)
また、プロセッサー110は、混雑率Xaと混雑率Yaがどれだけ異なった値であるかを示す値として、混雑率Xaと混雑率Yaの比βを求める。例えば、比βは、下式(3)で求めることができる。
β=Xa÷Ya (3)
この場合、プロセッサー110は、例えば、下式(4)を用いて混雑率Xbを求める。
Xb=Yb×β (4)
In step ST27, the processor 110 uses the target congestion rate and the fluctuation table acquired in step ST23 to estimate the congestion rate of a station section different from the target congestion rate. The target congestion rate is a known congestion rate used to estimate the congestion rate in step ST27. The target congestion rate is, for example, the congestion rate acquired in step ST22 and the congestion rate calculated in step ST26. The processor 110 uses one target congestion rate (hereinafter referred to as "congestion rate Xa") and a fluctuation table (hereinafter referred to as "target fluctuation table") with the same train and car as the congestion rate Xa to estimate the congestion rate of a station section different from the target congestion rate for the same train and car. In the target fluctuation table, the congestion rate of the station section that is the same as the congestion rate Xa is set to congestion rate Ya. The processor 110 calculates the difference between the congestion rate Xa and the congestion rate Ya. As an example, the processor 110 calculates the difference α between the congestion rates Xa and Ya as a value indicating how much the congestion rates Xa and Ya differ from each other. For example, the difference α can be calculated using the following formula (1):
α=Xa-Ya (1)
Here, the congestion rate to be estimated is defined as congestion rate Xb. Furthermore, the congestion rate in the target variation table for the same station section as congestion rate Xb is defined as congestion rate Yb. In this case, the processor 110 calculates congestion rate Xb using, for example, the following formula (2):
Xb=Yb+α (2)
The processor 110 also calculates the ratio β of the congestion rate Xa to the congestion rate Ya as a value indicating how much the congestion rate Xa and the congestion rate Ya differ from each other. For example, the ratio β can be calculated using the following equation (3).
β = Xa ÷ Ya (3)
In this case, the processor 110 calculates the congestion rate Xb using, for example, the following equation (4).
Xb = Yb × β (4)
一例として、プロセッサー110は、ステップST26で求めた混雑率C5と、混雑率C5と列車及び号車が同じ変動表とを用いて、混雑率A5、混雑率B5、混雑率D5~混雑率F6を求める。この場合、混雑率C5が混雑率Xaで、混雑率A5、混雑率B5、混雑率D5~混雑率F6が混雑率Xbである。また、この場合、C駅-D駅区間が第1の区間の一例であり、A駅-B駅区間、B駅-C駅区間、D駅-E駅区間、E駅-F駅区間及びF駅-G駅区間のそれぞれが第2の区間の一例である。 As an example, processor 110 uses the congestion rate C5 calculated in step ST26 and a fluctuation table that has the same train and car number as congestion rate C5 to calculate congestion rates A5, B5, D5 to F6. In this case, congestion rate C5 is congestion rate Xa, and congestion rates A5, B5, D5 to F6 are congestion rate Xb. In this case, the Station C-D station section is an example of a first section, and the Station A-B station section, the Station B-C station section, the Station D-E station section, the Station E-F station section, and the Station F-G station section are each an example of a second section.
ステップST28においてプロセッサー110は、ステップST27及びステップST28で予測した混雑率を混雑データベース141に記憶する。 In step ST28, the processor 110 stores the congestion rate predicted in steps ST27 and ST28 in the congestion database 141.
ステップST29においてプロセッサー110は、ステップST27及びステップST28で予測した混雑率を各種の装置に送信するように通信I/F150に対して指示する。当該各種の装置は、混雑率の推定値を用いて何らかの処理を行う装置である。当該各種の装置は、例えば、駅又は列車などに設置されたコンピューターなどである。当該コンピューターは、例えば、受信した混雑率を駅員又は列車の乗員などに報知する。当該コンピューターは、例えば、液晶ディスプレイなどのディスプレイに表示するなどして、混雑率を報知する。あるいは、当該各種の装置は、混雑率の推定値をユーザーに提供するサービスを提供するサーバーなどの装置である。当該サーバーは、ユーザーが使用するPC(personal computer)又はスマートホンなどの端末などに、受信した混雑率を送信する。
プロセッサー110は、ステップST29の処理の後、ステップST21へと戻る。
In step ST29, the processor 110 instructs the communication I/F 150 to transmit the congestion rate predicted in steps ST27 and ST28 to various devices. The various devices are devices that perform some kind of processing using the congestion rate estimate. The various devices are, for example, computers installed in stations or on trains. The computers, for example, notify station staff or train crew members of the received congestion rate. The computers notify the congestion rate by displaying it on a display such as an LCD display. Alternatively, the various devices are devices such as servers that provide a service that provides users with estimated congestion rates. The server transmits the received congestion rate to a terminal, such as a personal computer (PC) or smartphone, used by the user.
After processing step ST29, processor 110 returns to step ST21.
第1実施形態の混雑予測システム1によれば、混雑予測サーバー100は、時系列モデル及び混雑率データを用いて混雑率を予測する。このように混雑予測サーバー100は、新規な方法で混雑率を予測することができる。また、このような方法を用いることで、混雑予測サーバー100は、従来よりも高精度に高速に混雑率の予測することができると考えられる。したがって、混雑予測サーバー100は、大規模な路線全体の混雑率をリアルタイムに推定することができる。また、混雑予測サーバー100は、列車の車両毎など詳細な混雑率を推定することができる。未来の時刻の混雑率の変化を効率よく逐次予測することができる。また、混雑予測サーバー100は、大規模な路線ネットワークにおけるリアルタイムな混雑予測を効率よく、列車の車両単位で求めることができる。 According to the congestion prediction system 1 of the first embodiment, the congestion prediction server 100 predicts the congestion rate using a time series model and congestion rate data. In this way, the congestion prediction server 100 can predict the congestion rate using a novel method. Furthermore, by using such a method, it is believed that the congestion prediction server 100 can predict the congestion rate more accurately and quickly than conventional methods. Therefore, the congestion prediction server 100 can estimate the congestion rate of an entire large-scale route in real time. Furthermore, the congestion prediction server 100 can estimate detailed congestion rates, such as for each train car. It can efficiently and sequentially predict changes in congestion rate at future times. Furthermore, the congestion prediction server 100 can efficiently obtain real-time congestion predictions for each train car in a large-scale route network.
また、第1実施形態の混雑予測システム1によれば、混雑予測サーバー100は、時系列予測関数を求め、時系列予測関数を用いて混雑率を予測する。このように混雑予測サーバー100は、新規な方法で混雑率を予測することができる。また、このような方法を用いることで混雑予測サーバー100は、従来よりも高精度に高速に混雑率の予測することができると考えられる。 Furthermore, according to the congestion prediction system 1 of the first embodiment, the congestion prediction server 100 calculates a time series prediction function and predicts the congestion rate using the time series prediction function. In this way, the congestion prediction server 100 can predict the congestion rate using a novel method. Furthermore, by using such a method, it is believed that the congestion prediction server 100 can predict the congestion rate more accurately and quickly than conventional methods.
また、第1実施形態の混雑予測システム1によれば、混雑予測サーバー100は、変動表を用いて混雑率を予測する。このように混雑予測サーバー100は、新規な方法で混雑率を予測することができる。また、このような方法を用いることで、混雑予測サーバー100は、従来よりも高精度に高速に混雑率の予測することができると考えられる。 Furthermore, according to the congestion prediction system 1 of the first embodiment, the congestion prediction server 100 predicts the congestion rate using a fluctuation table. In this way, the congestion prediction server 100 can predict the congestion rate using a novel method. Furthermore, by using such a method, it is believed that the congestion prediction server 100 can predict the congestion rate more accurately and quickly than conventional methods.
〔第2実施形態〕
第2実施形態の混雑予測システムは、旅客フローを推定して列車平均混雑率を求める。
図9は第2実施形態に係る混雑予測システム1b及び混雑予測システム1bに含まれる構成要素の要部構成の一例を示すブロック図である。
混雑予測システム1bは、第1実施形態の混雑予測システム1に加えて改札システム500を含む。改札システム500は、ネットワークNWに接続する。
Second Embodiment
The congestion prediction system of the second embodiment estimates passenger flow to determine the average train congestion rate.
FIG. 9 is a block diagram showing an example of the main configuration of a congestion prediction system 1b according to the second embodiment and components included in the congestion prediction system 1b.
The congestion prediction system 1b includes a ticket gate system 500 in addition to the congestion prediction system 1 of the first embodiment. The ticket gate system 500 is connected to a network NW.
改札システム500は、自動改札機及びサーバーなどを含む。改札システム500は、旅客の改札内への入出場などの管理などを行うシステムである。改札システム500は、各旅客の入出場情報を記録し、混雑予測サーバー100に送信する。混雑予測サーバー100の通信I/F150は、当該入出場情報を受信する。プロセッサー110は、当該入出場情報を混雑データベース141に記憶する。なお、入出場情報は、旅客の入場駅、入場時間、出場駅及び出場時間などを示す情報含む。入出場情報は、経由駅を示す情報を含んでも良い。なお、入出場情報は、旅客がどこから前記乗り物に乗り、どこで前記乗り物から降りたかを示す乗降情報の一例である。また、旅客は、列車のユーザーの一例である。 The ticket gate system 500 includes an automatic ticket gate and a server. The ticket gate system 500 is a system that manages passenger entry and exit through ticket gates. The ticket gate system 500 records the entry and exit information of each passenger and transmits it to the congestion prediction server 100. The communication I/F 150 of the congestion prediction server 100 receives the entry and exit information. The processor 110 stores the entry and exit information in the congestion database 141. The entry and exit information includes information indicating the passenger's entry station, entry time, exit station and exit time. The entry and exit information may also include information indicating intermediate stations. The entry and exit information is an example of boarding and alighting information that indicates where the passenger boarded the vehicle and where he or she disembarked. The passenger is an example of a train user.
以下、第2実施形態に係る混雑予測システム1bの動作を図10及び図3などに基づいて説明する。図10は、混雑予測サーバー100のプロセッサー110による処理の一例を示すフローチャートである。プロセッサー110は、例えば、ROM120又は補助記憶装置140などに記憶されたプログラムに基づいて図10の処理を実行する。
プロセッサー110は、例えば、図10の処理を、混雑予測サーバー100の起動又は混雑予測用のアプリケーションの起動などに基づいて開始する。
The operation of the congestion prediction system 1b according to the second embodiment will be described below with reference to Fig. 10, Fig. 3, etc. Fig. 10 is a flowchart showing an example of processing by the processor 110 of the congestion prediction server 100. The processor 110 executes the processing of Fig. 10 based on a program stored in, for example, the ROM 120 or the auxiliary storage device 140.
The processor 110 starts the process of FIG. 10, for example, based on the activation of the congestion prediction server 100 or the activation of an application for congestion prediction.
プロセッサー110は、第1実施形態の図2の処理に代えて図10の処理を行う。また、プロセッサー110は、第1実施形態と同様に図3の処理を行う。第2実施形態の混雑予測システム1bの動作については、第1実施形態と異なる部分について説明する。 The processor 110 performs the processing shown in FIG. 10 instead of the processing shown in FIG. 2 in the first embodiment. The processor 110 also performs the processing shown in FIG. 3, as in the first embodiment. Regarding the operation of the congestion prediction system 1b in the second embodiment, differences from the first embodiment will be explained below.
プロセッサー110は、図10のステップST12の処理の後、ステップST31へと進む。
ステップST31においてプロセッサー110は、混雑データベース141から所定の期間分の入出場情報を取得する。ここで、プロセッサー110は、所定の期間の入出場情報のうち、運行が平常でない分を除いた、運行が平常である分の入出場情報を取得することが好ましい。なお、プロセッサー110は、例えば、どの期間の運行が平常であり、どの期間の運行が平常でないかを示す情報を運行管理システム200から取得する。そして、プロセッサー110は、当該情報を用いることで、入出場情報から運行が平常で無い分を除く。
After processing step ST12 in FIG. 10, processor 110 proceeds to step ST31.
In step ST31, the processor 110 acquires entry/exit information for a predetermined period from the congestion database 141. Here, it is preferable that the processor 110 acquires entry/exit information for the predetermined period for which operation is normal, excluding the portion of the entry/exit information for which operation is abnormal. Note that the processor 110 acquires, for example, information indicating which periods of operation are normal and which periods of operation are abnormal from the traffic management system 200. Then, the processor 110 uses this information to exclude the portion of the entry/exit information for which operation is abnormal.
ステップST32においてプロセッサー110は、ステップST31で取得した入出場情報を用いて、各駅区間を通る人数を時間帯別に推定する。当該入出場情報は、ステップST31で運行が平常でない分が除かれている場合には、運行が平常である分の入出場情報である。例えば、プロセッサー110は、ある駅STA1から入場して別のある駅STA2で出場した旅客がある駅区間SS1を通ったと推定できる。ただし、駅STA1から駅STA2までの経路が複数ある場合、駅STA1から入場して駅STA2で出場した旅客の全てが駅区間SS1を通ったわけではない場合がある。旅客がどの経路を通って駅STA1から駅STA2まで行くかは、予め割合が定められている。一例として、駅STA1から駅STA2までの経路が経路RO1~経路RO3の3種類ある場合、経路RO1を通る旅客が50%、経路RO2を通る旅客が30%、経路RO3を通る旅客が20%といったように割合が定められている。プロセッサー110は、この割合に従い、駅STA1から乗って駅STA2で降りた旅客のうち、何人が駅区間SS1を通ったか人数を推定する。一例として、経路RO1のみが駅区間SS1を含む場合、プロセッサー110は、駅STA1から入場して駅STA2で出場した旅客のうち、50%が駅区間SS1を通ったと推定する。また、プロセッサー110は、旅客が駅STA1に入場した時間から駅区間SS1を通った時間も推定する。一例として、プロセッサー110は、時刻t1から時刻t2までに入場した旅客は、時刻t3に駅区間SS1を通ると推定する。時刻t1及び時刻t2と時刻t3との関係は、例えば、時刻表などに基づき予め定められている。あるいは、プロセッサー110は、時刻表などから時刻t3を算出しても良い。プロセッサー110は、入場駅と出場駅との組み合わせそれぞれについて、同様に駅区間SS1を通る人数を推定する。そして、プロセッサー110は、当該組み合わせごとの駅区間SS1を通る人数を通る時間帯ごとに足し合わせることで、時間帯別の駅区間SS1を通る人数を推定する。また、プロセッサー110は、駅区間SS1以外の駅区間についても同様に時間帯別の通る人数を推定する。 In step ST32, processor 110 uses the entry/exit information acquired in step ST31 to estimate the number of people passing through each station section by time period. If the entry/exit information is for when operation is not normal in step ST31, it is the entry/exit information for when operation is normal. For example, processor 110 can estimate that a passenger who entered at station STA1 and exited at another station STA2 passed through station section SS1. However, if there are multiple routes from station STA1 to station STA2, not all passengers who entered at station STA1 and exited at station STA2 may have passed through station section SS1. The proportion of passengers who take which route from station STA1 to station STA2 is predetermined. As an example, if there are three routes from station STA1 to station STA2, route RO1 to route RO3, the ratios are set such that 50% of passengers take route RO1, 30% take route RO2, and 20% take route RO3. Based on these ratios, processor 110 estimates the number of passengers who boarded at station STA1 and disembarked at station STA2 and passed through station section SS1. As an example, if only route RO1 includes station section SS1, processor 110 estimates that 50% of passengers who entered station STA1 and exited at station STA2 passed through station section SS1. Processor 110 also estimates the time that passengers passed through station section SS1 from the time they entered station STA1. As an example, processor 110 estimates that passengers who entered between time t1 and time t2 will pass through station section SS1 at time t3. The relationship between time t1, time t2, and time t3 is determined in advance based on, for example, a timetable. Alternatively, processor 110 may calculate time t3 from a timetable. Processor 110 similarly estimates the number of people passing through station section SS1 for each combination of entry station and exit station. Processor 110 then adds up the number of people passing through station section SS1 for each combination for each time period, thereby estimating the number of people passing through station section SS1 for each time period. Processor 110 also similarly estimates the number of people passing through station sections other than station section SS1 for each time period.
ステップST33においてプロセッサー110は、ステップST32で推定した時間帯別の通る人数を用いて、時間帯別の列車平均混雑率を求める。例えば、駅区間SS1のある時間帯の通る人数がNP1であるとする。また、駅区間SS1を当該時間帯に通る列車の定員を足し合わせた数をCA1とする。この場合、当該時間帯の列車平均混雑率は、例えば、下式(5)で求めることができる。
(列車平均混雑率)=NP1÷CA1 (5)
プロセッサー110は、各駅区間の各時間帯について同様に(5)式などを用いて列車平均混雑率を求める。プロセッサー110は、ステップST33の処理の後、ステップST14へと進む。なお、第2実施形態では、プロセッサー110は、第1実施形態のステップST13で求める列車平均混雑率に代えて、ステップST33で求める列車平均混雑率を用いて各処理を行う。
In step ST33, the processor 110 calculates the average train congestion rate for each time period using the number of passengers for each time period estimated in step ST32. For example, assume that the number of passengers passing through station section SS1 during a certain time period is NP1. Also, assume that the total number of passengers on trains passing through station section SS1 during that time period is CA1. In this case, the average train congestion rate for that time period can be calculated, for example, using the following formula (5):
(Average train congestion rate) = NP1 ÷ CA1 (5)
The processor 110 similarly calculates the average train congestion rate for each time period in each station section using equation (5) or the like. After processing step ST33, the processor 110 proceeds to step ST14. Note that in the second embodiment, the processor 110 performs each process using the average train congestion rate calculated in step ST33 instead of the average train congestion rate calculated in step ST13 of the first embodiment.
第2実施形態の混雑予測システム1bは、第1実施形態と同様の効果が得られる。 The congestion prediction system 1b of the second embodiment achieves the same effects as the first embodiment.
また、第2実施形態の混雑予測システム1bによれば、混雑予測サーバー100は、入出場情報を用いて混雑率を求める。これにより、混雑予測サーバー100は、混雑率をより精度良く求めることができる。 Furthermore, according to the congestion prediction system 1b of the second embodiment, the congestion prediction server 100 calculates the congestion rate using entry/exit information. This allows the congestion prediction server 100 to calculate the congestion rate with greater accuracy.
〔第3実施形態〕
第3実施形態の混雑予測システムは、走行中の列車の混雑率を用いて変動表を補正する。
第3実施形態の混雑予測システム1bの構成は、第2実施形態と同様であるので説明を省略する。なお、第3実施形態の混雑予測システムの構成は、第1実施形態と同様であっても良い。
Third Embodiment
The congestion prediction system of the third embodiment corrects the fluctuation table using the congestion rate of a running train.
The configuration of the congestion prediction system 1b of the third embodiment is the same as that of the second embodiment, and therefore description thereof will be omitted. Note that the configuration of the congestion prediction system of the third embodiment may be the same as that of the first embodiment.
以下、第3実施形態に係る混雑予測システム1bの動作を図10及び図11などに基づいて説明する。図11は、混雑予測サーバー100のプロセッサー110による処理の一例を示すフローチャートである。プロセッサー110は、例えば、ROM120又は補助記憶装置140などに記憶されたプログラムに基づいて図11の処理を実行する。
プロセッサー110は、例えば、図11の処理を、混雑予測サーバー100の起動又は混雑予測用のアプリケーションの起動などに基づいて開始する。
The operation of the congestion prediction system 1b according to the third embodiment will be described below with reference to Fig. 10 and Fig. 11. Fig. 11 is a flowchart showing an example of processing by the processor 110 of the congestion prediction server 100. The processor 110 executes the processing of Fig. 11 based on a program stored in the ROM 120 or the auxiliary storage device 140, for example.
The processor 110 starts the process of FIG. 11, for example, based on the activation of the congestion prediction server 100 or the activation of an application for congestion prediction.
プロセッサー110は、第2実施形態の図3の処理に代えて図11の処理を行う。また、第3実施形態のプロセッサー110は、第2実施形態と同様に図10の処理を行う。第3実施形態の混雑予測システム1bの動作については、第2実施形態と異なる部分について説明する。 The processor 110 performs the processing shown in FIG. 11 instead of the processing shown in FIG. 3 in the second embodiment. Furthermore, the processor 110 in the third embodiment performs the processing shown in FIG. 10, similar to the second embodiment. Regarding the operation of the congestion prediction system 1b in the third embodiment, differences from the second embodiment will be explained below.
プロセッサー110は、図11のステップST26の処理の後、ステップST41へと進む。
ステップST41においてプロセッサー110は、混雑率の予測対象の列車について、始発駅から現在走行中の位置までのうち、計測装置300によって計測済みの混雑率を混雑率管理システム400などから取得する。一例として、混雑率D5を予測する場合には、混雑率D5に対応する列車が予測対象の列車である。なお、ステップST41で取得される混雑率を以下「実測混雑率」という。実測混雑率は、予測対象の乗り物の混雑の状態の第2の計測値の一例である。
After processing step ST26 in FIG. 11, processor 110 proceeds to step ST41.
In step ST41, the processor 110 acquires, from the congestion rate management system 400 or the like, the congestion rate measured by the measurement device 300 for the train whose congestion rate is to be predicted, from the starting station to the current running position. As an example, when predicting congestion rate D5, the train corresponding to congestion rate D5 is the train to be predicted. The congestion rate acquired in step ST41 will be referred to as the "measured congestion rate" hereinafter. The measured congestion rate is an example of a second measurement value of the congestion state of the vehicle to be predicted.
ステップST42においてプロセッサー110は、実測混雑率を用いて対象変動表を補正する。例えば、プロセッサー110は、対象変動表の各駅区間の混雑率について、実測混雑率との差又は比などを求める。そして、プロセッサー110は、当該差又は比を用いて最小二乗法、AR、ARMA、ARIMA又は状態空間モデルなどを使用して、予測対象の列車がまだ到達していない駅区間それぞれについての差又は比を予測する。そして、プロセッサー110は、予測した差又は比を用いて変動表を補正する。プロセッサー110は、ステップST42の処理の後、ステップST27へと進む。プロセッサー110は、補正済みの変動表を用いてステップST27の処理を行う。 In step ST42, processor 110 corrects the target fluctuation table using the measured congestion rate. For example, processor 110 calculates the difference or ratio between the congestion rate of each station section in the target fluctuation table and the measured congestion rate. Then, processor 110 uses the difference or ratio to predict the difference or ratio for each station section that the train to be predicted has not yet arrived, using least squares, AR, ARMA, ARIMA, or a state space model. Then, processor 110 corrects the fluctuation table using the predicted difference or ratio. After processing step ST42, processor 110 proceeds to step ST27. Processor 110 performs the processing of step ST27 using the corrected fluctuation table.
第3実施形態の混雑予測システム1bは、第2実施形態と同様の効果が得られる。 The congestion prediction system 1b of the third embodiment achieves the same effects as the second embodiment.
また、第3実施形態の混雑予測システム1bによれば、混雑予測サーバー100は、計測装置300によって計測された、混雑率の予測対象の列車の混雑率を用いて変動表を補正する。これにより、混雑予測サーバー100は、混雑率をより精度良く求めることができる。 Furthermore, according to the congestion prediction system 1b of the third embodiment, the congestion prediction server 100 corrects the fluctuation table using the congestion rate of the train for which the congestion rate is to be predicted, measured by the measurement device 300. This allows the congestion prediction server 100 to calculate the congestion rate with greater accuracy.
〔第4実施形態〕
第4実施形態の混雑予測システムは、鉄道の運行に異常が発生した場合に、平常時とは異なる方法で混雑率を推定する。また、第4実施形態の混雑予測システムは、異常時の混雑率から平常時の混雑率へ復旧する見込み時間を推定する。
異常が発生した路線は、輸送手段の運行が大きく制約される場合が多い。また、異常の原因及び障害対策方法が多様である。さらに、旅客への対応も情報提供の内容によって多様となる。このため、異常が発生した路線に対する混雑状態の予測は、他の路線の混雑状態の予測に比べて難しい。したがって、例えば、異常が発生した路線に対しては、運転再開見込みの情報提供を優先することが好ましい。そして、異常が発生した路線に繋がる路線又は近隣の路線などへの混雑の影響、混雑状態の異常の予測及び平常状態への復旧時間の予測を第4実施形態の混雑予測システムを使用して推測することが考えられる。ただし、異常が発生した路線に対して第4実施形態の混雑予測システムを用いることも可能である。
Fourth Embodiment
The congestion prediction system of the fourth embodiment estimates the congestion rate when an abnormality occurs in railway operation using a method different from that used during normal times. The congestion prediction system of the fourth embodiment also estimates the expected time for the congestion rate during the abnormal time to return to the normal congestion rate.
In many cases, transportation operations are significantly restricted on routes where an abnormality has occurred. Furthermore, there are a variety of causes of the abnormality and methods for dealing with the problem. Furthermore, passenger responses vary depending on the content of the information provided. For this reason, predicting congestion on routes where an abnormality has occurred is more difficult than predicting congestion on other routes. Therefore, for example, it is preferable to prioritize providing information on the expected resumption of service for routes where an abnormality has occurred. It is also possible to use the congestion prediction system of the fourth embodiment to estimate the impact of congestion on routes connected to the route where an abnormality has occurred or neighboring routes, predict the abnormal congestion status, and predict the time it will take to return to normal. However, the congestion prediction system of the fourth embodiment can also be used for routes where an abnormality has occurred.
第4実施形態の混雑予測システム1bの構成は、第3実施形態と同様であるので説明を省略する。 The configuration of the congestion prediction system 1b in the fourth embodiment is the same as that in the third embodiment, so a detailed description will be omitted.
以下、第4実施形態に係る混雑予測システム1bの動作を図10~図12などに基づいて説明する。図12は、混雑予測サーバー100のプロセッサー110による処理の一例を示すフローチャートである。プロセッサー110は、例えば、ROM120又は補助記憶装置140などに記憶されたプログラムに基づいて図12の処理を実行する。
プロセッサー110は、例えば、図12の処理を、混雑予測サーバー100の起動又は混雑予測用のアプリケーションの起動などに基づいて開始する。
The operation of the congestion prediction system 1b according to the fourth embodiment will be described below with reference to Figures 10 to 12. Figure 12 is a flowchart showing an example of processing by the processor 110 of the congestion prediction server 100. The processor 110 executes the processing of Figure 12 based on a program stored in, for example, the ROM 120 or the auxiliary storage device 140.
The processor 110 starts the process of FIG. 12, for example, based on the activation of the congestion prediction server 100 or the activation of an application for congestion prediction.
第4実施形態のプロセッサー110は、第3実施形態と同様に図10及び図11の処理を行う。また、プロセッサー110は、第3実施形態の図10の処理に加えて図12の処理を行う。また、第4実施形態の混雑予測システム1bの動作については、第3実施形態と異なる部分について説明する。 The processor 110 of the fourth embodiment performs the processing shown in Figures 10 and 11, similar to that of the third embodiment. Furthermore, the processor 110 performs the processing shown in Figure 12 in addition to the processing shown in Figure 10 of the third embodiment. Regarding the operation of the congestion prediction system 1b of the fourth embodiment, differences from the third embodiment will be described below.
ステップST51においてプロセッサー110は、予測対象の路線とは異なる路線又は予測対象の路線の運行に異常が発生したか否かを判定する。例えば、運行管理システム200は、各路線の運行に異常が発生した場合、運行に異常が発生したことを通知する情報を混雑予測サーバー100に送信する。当該情報は、混雑予測サーバー100の通信I/F150によって受信される。プロセッサー110は、当該情報が受信されたことに応じて路線の運行に異常が発生したと判定する。あるいは、プロセッサー110は、混雑率の計測値の変化から、異常が発生したことを判定しても良い。 In step ST51, the processor 110 determines whether an abnormality has occurred in the operation of a route different from the route being predicted or the route being predicted. For example, when an abnormality occurs in the operation of each route, the operation management system 200 transmits information notifying the congestion prediction server 100 that an abnormality has occurred in the operation. The information is received by the communication I/F 150 of the congestion prediction server 100. In response to receiving the information, the processor 110 determines that an abnormality has occurred in the operation of the route. Alternatively, the processor 110 may determine that an abnormality has occurred from a change in the measured value of the congestion rate.
図13は、混雑率の時間変化の一例を示すグラフである。図13のグラフは、縦軸が混雑率、横軸が時間である。図13には、混雑率の計測値を点Pcとしてプロットしている。また、折れ線Lfは、点Pcを結んだ折れ線である。そして、曲線Lgは、正常時の混雑率の時間変化を示す。 Figure 13 is a graph showing an example of how congestion rates change over time. In the graph in Figure 13, the vertical axis represents congestion rates and the horizontal axis represents time. In Figure 13, measured values of congestion rates are plotted as points Pc. The broken line Lf connects the points Pc. The curve Lg shows how congestion rates change over time under normal conditions.
点Pcは、例えば、特定の駅区間の列車又は車両の混雑率である。この場合、折れ線Lfは、例えば、特定の駅区間における列車又は車両の混雑率の時間変化を示す。またこの場合、曲線Lgは、例えば時系列モデルである。
あるいは、点Pcは、特定の列車又は車両の時間ごとの混雑率である。この場合、折れ線Lfは、例えば、特定の列車又は車両の混雑率の時間変化を示す。またこの場合、曲線Lgは、例えば変動表をグラフ化したものである。
なお、曲線Lgは、通常時用の変動表を用いて予測された混雑率のグラフであっても良い。
Point Pc represents, for example, the congestion rate of trains or vehicles in a specific station section. In this case, polygonal line Lf represents, for example, the change in the congestion rate of trains or vehicles in a specific station section over time. In this case, curve Lg represents, for example, a time series model.
Alternatively, point Pc represents the hourly congestion rate of a specific train or vehicle. In this case, polygonal line Lf represents, for example, the change in the congestion rate of a specific train or vehicle over time. In addition, curve Lg represents, for example, a graph of a fluctuation table.
The curve Lg may be a graph of the congestion rate predicted using a fluctuation table for normal times.
また、図13には差分混雑率CDも示している。差分混雑率CDは、同じ時間の点Pcと曲線Lgとの差又は曲線Lfと曲線Lgとの差である。図13では、代表的に、1つの点Pcについての差分混雑率CDを示している。 Figure 13 also shows the differential congestion rate CD. The differential congestion rate CD is the difference between point Pc and curve Lg at the same time, or the difference between curve Lf and curve Lg. Figure 13 representatively shows the differential congestion rate CD for one point Pc.
図14は、差分混雑率CDの時間変化の一例を示すグラフである。点Pdは、点Pcと曲線Lgとの差である差分混雑率CDをプロットしたものである。また、線Lhは、点Pdを結んだ折れ線である。したがって、線Lhは、差分混雑率CDの時間変化を示す。また、図14には、閾値CDaも示している。閾値CDaは、異常が発生しているか否かの境界を示す閾値である。
プロセッサー110は、例えば、最新の混雑率の計測値に基づく差分混雑率CDが閾値CDa以上であるならば、路線の運行に異常が発生したと判定する。図13及び図14には、差分混雑率CDが閾値CDa以上となった時刻として、時刻Taを示している。
FIG. 14 is a graph showing an example of the change over time in the differential congestion rate CD. Point Pd is a plot of the differential congestion rate CD, which is the difference between point Pc and curve Lg. Line Lh is a polygonal line connecting points Pd. Therefore, line Lh shows the change over time in the differential congestion rate CD. FIG. 14 also shows a threshold value CDa. Threshold value CDa is a threshold value that indicates the boundary between whether or not an abnormality has occurred.
For example, if the differential congestion rate CD based on the latest congestion rate measurement value is equal to or greater than a threshold value CDa, the processor 110 determines that an abnormality has occurred in the operation of the route. 13 and 14 show time Ta as the time when the differential congestion rate CD becomes equal to or greater than the threshold value CDa.
プロセッサー110は、路線の運行に異常が発生したと判定しないならば、ステップST51においてNoと判定してステップST51の処理を繰り返す。対して、プロセッサー110は、路線の運行に異常が発生したと判定するならば、ステップST51においてYesと判定してステップST52へと進む。 If processor 110 does not determine that an abnormality has occurred in the operation of the route, it determines No in step ST51 and repeats the processing of step ST51. On the other hand, if processor 110 determines that an abnormality has occurred in the operation of the route, it determines Yes in step ST51 and proceeds to step ST52.
ステップST52においてプロセッサー110は、図11のステップST26における混雑率の予測方法を通常時用のものから異常時用のものに変更する。異常時用の混雑率の予測方法を用いる場合、プロセッサー110は、例えば、ステップST26の処理において求めた混雑率の予測値をn1倍する。あるいは、プロセッサー110は、ステップST26の処理において求めた混雑率にm1ポイント足す。なお、n1は、1以上の数である。m1は、0以上の数である。また、混雑率には上限があっても良い。混雑率に上限を設ける場合、プロセッサー110は、混雑率をn1倍にすることで混雑率が上限以上となる場合、又は混雑率にm1ポイント足すことで混雑率が上限以上となる場合、計算結果に拘らず混雑率を上限の値にする。混雑率の上限は、例えば、満員状態を示す。すなわち、混雑率が上限である列車及び車両は満員状態となることが予想される。なお、駅区間ごとにn1又はm1の値が異なっていても良い。また、異常が発生した路線及び駅区間によってn1又はm1の値が異なっていても良い。また、混雑率の予測対象の路線及び駅区間によってn1又はm1の値が異なっていても良い。また、プロセッサー110は、ステップST26の処理において求めた混雑率に基づく差分混雑率をn1倍するかm1足しても良い。そして、プロセッサー110は、当該差分混雑率から混雑率を求める。
図13に示す混雑率が平坦となっている部分は、混雑率が上限となっている部分である。混雑率が平坦となっている部分は、時刻Tcから時刻Tdの範囲である。
In step ST52, the processor 110 changes the congestion rate prediction method used in step ST26 of FIG. 11 from the normal congestion rate prediction method to the abnormal congestion rate prediction method. When using the abnormal congestion rate prediction method, the processor 110, for example, multiplies the congestion rate prediction value calculated in step ST26 by n1. Alternatively, the processor 110 adds m1 points to the congestion rate calculated in step ST26. Here, n1 is a number equal to or greater than 1. Here, m1 is a number equal to or greater than 0. The congestion rate may also have an upper limit. When setting an upper limit for the congestion rate, if multiplying the congestion rate by n1 or adding m1 points to the congestion rate results in the congestion rate being equal to or greater than the upper limit, the processor 110 sets the congestion rate to the upper limit regardless of the calculation result. The upper limit for the congestion rate indicates, for example, a full-occupancy state. In other words, trains and cars with a maximum congestion rate are expected to be full. The values of n1 and m1 may differ for each station section. The value of n1 or m1 may also differ depending on the line and station section where the abnormality occurred. The value of n1 or m1 may also differ depending on the line and station section for which the congestion rate is to be predicted. The processor 110 may also multiply the differential congestion rate based on the congestion rate calculated in the processing of step ST26 by n1 or add m1 to it. The processor 110 then calculates the congestion rate from the differential congestion rate.
13, the portion where the congestion rate is flat is the portion where the congestion rate is at its upper limit, and is the range from time Tc to time Td.
図12のステップST53においてプロセッサー110は、混雑率が正常化し始めたか否かを判定する。運行の異常を原因として混雑率は平常時より増加するが、運行が正常化するにつれて混雑率の増加量は減少する。プロセッサー110は、例えば、計測装置300で計測された混雑率に基づき、混雑率の増加量が減少し始めたことを検知して、混雑率が正常化し始めたと判定する。あるいは、プロセッサー110は、運行管理システム200から混雑率が正常化し始めたことを示す情報を受信したことに応じて混雑率が正常化し始めたと判定する。混雑率が正常化し始めたことを示す情報は、一例として、列車の運転再開を示す情報である。 In step ST53 of FIG. 12, processor 110 determines whether the congestion rate has started to return to normal. Although the congestion rate increases from normal due to an operational abnormality, the rate of increase in the congestion rate decreases as operation returns to normal. Processor 110 determines that the congestion rate has started to return to normal by detecting, for example, that the rate of increase in the congestion rate has started to decrease based on the congestion rate measured by measuring device 300. Alternatively, processor 110 determines that the congestion rate has started to return to normal in response to receiving information from traffic management system 200 indicating that the congestion rate has started to return to normal. One example of information indicating that the congestion rate has started to return to normal is information indicating that train operations have resumed.
図14には、閾値CDbを示している。閾値CDbは、混雑率が正常化し始めたか否かの境界を示す閾値である。プロセッサー110は、例えば、差分混雑率CDが閾値CDb以下となったならば混雑率が正常化し始めたと判定する。 Figure 14 shows the threshold CDb. The threshold CDb is a threshold that indicates the boundary between whether the congestion rate has started to return to normal. For example, if the differential congestion rate CD becomes equal to or less than the threshold CDb, the processor 110 determines that the congestion rate has started to return to normal.
プロセッサー110は、混雑率が正常化し始めたと判定しないならば、ステップST53においてNoと判定してステップST53の処理を繰り返す。対して、プロセッサー110は、混雑率が正常化し始めたと判定するならば、ステップST53においてYesと判定してステップST54へと進む。 If processor 110 does not determine that the congestion rate has begun to return to normal, it determines No in step ST53 and repeats the processing of step ST53. On the other hand, if processor 110 determines that the congestion rate has begun to return to normal, it determines Yes in step ST53 and proceeds to step ST54.
ステップST54においてプロセッサー110は、図11のステップST26における混雑率の予測方法を異常時用のものから正常化時用のものに変更する。正常化時用の予測方法は、例えば、異常時から正常時への過渡状態に用いるための方法である。正常化時用の混雑率の予測方法を用いる場合、プロセッサー110は、ステップST26の処理において求めた混雑率の予測値をn2倍する。あるいは、プロセッサー110は、ステップST26処理において求めた混雑率にm2ポイント足す。なお、n2は、1以上の数である。m2は、0以上の数である。n2及びm2は、正常化し始めてからの経過時間t11が長いほど減少する値である。t11が0である場合、例えば、n1=n2、m1=m2である。一例として、n2及びm2は、
n2=1+(n1-1)×e(-t11/τ1) (6)
m2=m1×e(-t11/τ2) (7)
である。ここで、eは、ネイピア数である。また、τ1及びτ2は、時定数である。なお、n2を示す関数は、(6)式に限らない。また、m2を示す関数は、(7)式に限らない。また、プロセッサー110は、ステップST26の処理において求めた混雑率に基づく差分混雑率をn2倍するかm2足しても良い。そして、プロセッサー110は、当該差分混雑率から混雑率を求める。なお、図14の減衰曲線Liは、ステップST26の処理において求めた混雑率に基づく差分混雑率をn2倍するかm2足した差分混雑率の時間変化を示す。
また、図14に示す曲線Liは、n2又はm2を用いて求めた混雑率に基づく差分混雑率CDの時間変化の一例を示す曲線である。
In step ST54, the processor 110 changes the congestion rate prediction method in step ST26 of FIG. 11 from that for abnormal times to that for normal times. The normal time prediction method is, for example, a method to be used in a transient state from an abnormal time to a normal time. When using the congestion rate prediction method for normal times, the processor 110 multiplies the predicted value of the congestion rate obtained in the processing of step ST26 by n2. Alternatively, the processor 110 adds m2 points to the congestion rate obtained in the processing of step ST26. Note that n2 is a number equal to or greater than 1. m2 is a number equal to or greater than 0. n2 and m2 are values that decrease as the elapsed time t11 from the start of normalization becomes longer. When t11 is 0, for example, n1 = n2 and m1 = m2. As an example, n2 and m2 are
n2=1+(n1-1)×e (-t11/τ1) (6)
m2=m1×e (-t11/τ2) (7)
where e is Napier's constant. τ1 and τ2 are time constants. Note that the function representing n2 is not limited to equation (6). Note that the function representing m2 is not limited to equation (7). The processor 110 may also multiply the differential congestion rate based on the congestion rate calculated in the processing of step ST26 by n2 or add m2 to it. The processor 110 then calculates the congestion rate from the differential congestion rate. Note that the attenuation curve Li in FIG. 14 shows the time change of the differential congestion rate calculated by multiplying the differential congestion rate based on the congestion rate calculated in the processing of step ST26 by n2 or adding m2 to it.
A curve Li shown in FIG. 14 is a curve showing an example of a change over time in the differential congestion rate CD based on the congestion rate calculated using n2 or m2.
あるいは、プロセッサー110は、i番目の点Pdからの経過時間をt12とし、n2に代えてn3を、m2に代えてm3を用いても良い。一例として、n3及びm3は、
n3=1+(n1-1)×e(-t12/τ3) (8)
m3=m1×e(-t12/τ4) (9)
である。
Alternatively, the processor 110 may set the elapsed time from the i-th point Pd to t12, and use n3 instead of n2 and m3 instead of m2.
n3=1+(n1-1)×e (-t12/τ3) (8)
m3=m1×e (-t12/τ4) (9)
is.
プロセッサー110は、例えば図14に示すような点Pdを2つ以上用いて時定数τを求める。例えば、混雑率が正常化し始めてから最初の点Pdの差分混雑率CDをR1、2番目の点Pdの差分混雑率CDをR2、3番目の点Pdの差分混雑率CDをR3、…とする。また、(i-1)番目の点Pdの差分混雑率CDをR(i-1)、i番目の点Pdの差分混雑率CDをRiとする。なお、iは、1以上N以下の任意の整数である。Nは、時定数τを求めるために用いることができる点Pdの個数である。また、i番目の点Pdの時刻と(i-1)番目の点Pdの時刻との差はΔTiである。なお、ΔTiを総称してΔTというものとする。ΔTiの値は、一定ではなくiごとに異なる。あるいは、プロセッサー110は、ΔTiがiにかかわらず一定となるように補正しても良い。時定数τは、例えば下式(10)で表すことができる。なお、時定数τ1~時定数τ4を総称して時定数τというものとする。また、εは、推定誤差である。
(R(i-1)/Ri) ∝ τ+ε (10)
プロセッサー110は、推定誤差εを、i=1~Nの範囲で平均値処理することで圧縮することにより求める。あるいは、プロセッサー110は、平均値処理をすることで、状態変化の追随を考慮して移動平均をとっても良い。
The processor 110 calculates the time constant τ using two or more points Pd, such as those shown in FIG. 14 . For example, let R1 be the differential congestion rate CD at the first point Pd after the congestion rate begins to normalize, R2 be the differential congestion rate CD at the second point Pd, R3 be the differential congestion rate CD at the third point Pd, and so on. Let R(i-1) be the differential congestion rate CD at the (i-1)th point Pd, and Ri be the differential congestion rate CD at the i-th point Pd. Note that i is an integer between 1 and N. N is the number of points Pd that can be used to calculate the time constant τ. The difference between the time of the i-th point Pd and the time of the (i-1)th point Pd is ΔTi. Note that ΔTi is collectively referred to as ΔT. The value of ΔTi is not constant but varies for each i. Alternatively, the processor 110 may correct ΔTi so that it is constant regardless of i. The time constant τ can be expressed, for example, by the following equation (10): The time constants τ1 to τ4 are collectively referred to as the time constant τ, and ε is the estimation error.
(R(i-1)/Ri) ∝ τ+ε (10)
The processor 110 obtains the estimation error ε by compressing it through averaging in the range of i = 1 to N. Alternatively, the processor 110 may perform averaging to obtain a moving average in consideration of tracking state changes.
Riは、例えば、下式(11)で表すことができる。
Ri=(R(i-1))×e(-ΔTi/τ) (11)
したがって、プロセッサー110は、時定数τを例えば下式(12)で求める。なお、ここで求められるτは、例えば時定数τ3又は時定数τ4である。また、lnは、自然対数を示す。
τi=ΔTi×ln(Ri/(R(i-1))) (12)
なお、時定数τiは、Ri及びR(i-1)を用いて求めた時定数τを示す。ここで、i=1の場合について考える。i=1の場合、R(i-1)は、R0である。R0は、R1に対応する点Pdの1つ前の点Pdの差分混雑率CDを示す。なお、R0は、図14には示していない。プロセッサー110は、R0に対応する点Pdに対応する点Pcの混雑率が上限である場合、誤差が大きくなるため、R2の値が分かるまで時定数τを求めなくても良い。そして、プロセッサー110は、R2の値を求めるなどしてR2の値が分かったならばi=2の場合の(12)式により時定数τを求める。なお、点Pdに対応する点Pcとは、曲線Lgの値を引いて差分混雑率CDを求める前の点Pcである。したがって、点Pdと点Pdに対応する点Pcとは時間は同じである。プロセッサー110は、時定数τを求めるまでは、混雑率の予測方法を異常時用のものから正常化時用のものに変更しなくても良い。
Ri can be expressed, for example, by the following formula (11).
Ri=(R(i-1))×e (-ΔTi/τ) (11)
Therefore, the processor 110 calculates the time constant τ using, for example, the following equation (12): Here, τ is, for example, the time constant τ3 or the time constant τ4. Furthermore, ln indicates the natural logarithm.
τi = ΔTi × ln(Ri / (R(i-1))) (12)
Note that the time constant τi denotes the time constant τ calculated using Ri and R(i-1). Here, consider the case where i = 1. When i = 1, R(i-1) is R0. R0 indicates the differential congestion rate CD of the point Pd immediately preceding the point Pd corresponding to R1. Note that R0 is not shown in FIG. 14. If the congestion rate of the point Pc corresponding to the point Pd corresponding to R0 is at the upper limit, the error will be large, so the processor 110 does not need to calculate the time constant τ until the value of R2 is known. Then, once the value of R2 is known, for example, the processor 110 calculates the time constant τ using equation (12) for i = 2. Note that the point Pc corresponding to point Pd is the point Pc before the differential congestion rate CD is calculated by subtracting the value of the curve Lg. Therefore, the time between point Pd and the point Pc corresponding to point Pd is the same. The processor 110 does not need to change the congestion rate prediction method from that for abnormal times to that for normal times until the time constant τ is calculated.
混雑率の計測値を用いて求めた差分混雑率は、様々な誤差が含まれている。時定数τの推定値τeは、
τe=E[τi+εi]=E[τi]+E[εi] (13)
と表すことができる。ここで、E[x]はxの期待値を示す。なお、τiとεiは独立とする。また、εiは、時定数τiについての推定誤差である。また、εiは、各測定値などに含まれる誤差をまとめたものである。プロセッサー110は、次の式によりτeを求めても良い。
τe≒(1/N)・Σ(ΔTi・ln(Ri/R(i-1)))+(1/N)・Σ(εi) (14)
ここで、右辺の第2項(1/N)・Σ(εi)は、0に近いため無視しても良い。すなわち、当該第2項を無視する場合、プロセッサー110は、N個のτiの平均値をτeとする。プロセッサー110は、時定数の推定値τeを時定数τとして用いても良い。
The differential congestion rate calculated using the measured congestion rate value contains various errors. The estimated value τe of the time constant τ is
τe=E[τi+εi]=E[τi]+E[εi] (13)
Here, E[x] represents the expected value of x. Note that τi and εi are independent. εi is the estimation error for the time constant τi. εi is a sum of the errors included in each measurement value, etc. The processor 110 may calculate τe using the following formula:
τe≒(1/N)・Σ(ΔTi・ln(Ri/R(i-1)))+(1/N)・Σ(εi) (14)
Here, the second term on the right-hand side, (1/N)·Σ(εi), may be ignored because it is close to 0. In other words, when the second term is ignored, the processor 110 sets the average value of N τi as τe. The processor 110 may use the estimated value τe of the time constant as the time constant τ.
ステップST55においてプロセッサー110は、運行が正常になるまでの時間ΔTb及び運行が正常になる時刻Tbを推定する。時刻Tbは、例えば、ステップST54でn2をかけるかm2を足すことによって求めた混雑率に基づく差分混雑率が閾値CDa以下となる時刻である。運行が正常になるまでの時間ΔTbは、現在時刻から時刻Tbまでの時間である。 In step ST55, processor 110 estimates the time ΔTb until operation returns to normal and the time Tb at which operation returns to normal. Time Tb is, for example, the time at which the differential congestion rate based on the congestion rate calculated in step ST54 by multiplying by n2 or adding m2 falls below threshold CDa. The time ΔTb until operation returns to normal is the time from the current time to time Tb.
あるいは、プロセッサー110は、Ri及びR(i-1)から時刻Tbを推定しても良い。例えば、プロセッサー110は、i番目の点Pdと(i-1)番目の点Pdを結んだ直線SLiの傾きmを求める。傾きmは、例えば、
m=((Ri-R(i-1))/ΔTi) (15)
で求めることができる。また、プロセッサー110は、i番目の点Pdの時刻Tiと点Qの時刻TQとの差の時間Txiを求める。なお、点Qは、直線SLiが閾値CDaを示す直線と交わる点である。換言すると、点Qは、直線SLi上の、差分混雑率CDがCDaである点である。時間Txiは、例えば、
Txi=ΔTi-Ri/m (16)
で求めることができる。プロセッサー110は、例えば、時刻Tiから時刻Tbまでの時間Tbiを求める。時間Tbiは、例えば
Tbi=Txi×γ1 (17)
で求めることができる。ここで、係数γ1は、予め定められた所定の係数である。また、プロセッサー110は、時刻Tbを推定する。時刻Tbは、例えば、
Tb=Ti+Tbi (18)
で求めることができる。
Alternatively, the processor 110 may estimate the time Tb from Ri and R(i-1). For example, the processor 110 calculates the slope m of the straight line SLi connecting the i-th point Pd and the (i-1)-th point Pd. The slope m can be calculated, for example, as follows:
m=((Ri-R(i-1))/ΔTi) (15)
The processor 110 also calculates the time Txi, which is the difference between the time Ti of the i-th point Pd and the time TQ of the point Q. The point Q is the point where the line SLi intersects with the line indicating the threshold CDa. In other words, the point Q is the point on the line SLi where the differential congestion rate CD is CDa. The time Txi can be calculated, for example, as follows:
Txi=ΔTi−Ri/m (16)
The processor 110 calculates, for example, the time Tbi from the time Ti to the time Tb. The time Tbi can be calculated, for example, by the following equation: Tbi=Txi×γ1 (17)
Here, the coefficient γ1 is a predetermined coefficient. The processor 110 also estimates the time Tb. The time Tb can be calculated, for example, as follows:
Tb=Ti+Tbi (18)
It can be found by:
あるいは、プロセッサー110は、時間Txiに代えて時間Tx(i-1)を、時間Tbiに代えて時間Tb(i-1)を求めても良い。時間Tx(i-1)は、(i-1)番目の点Pdの時刻Tiと点Qの時刻TQとの差である。時間Tb(i-1)は、(i-1)番目の点Pdの時刻T(i-1)から時刻Tbまでの時間である。時間Tx(i-1)は、
Tx(i-1)=ΔT(i-1)-(R(i-1))/m (19)
で求めることができる。また、時間Tb(i-1)は、例えば、
Tb(i-1)=Tx(i-1)×γ2 (20)
で求めることができる。ここで、係数γ2は、予め定められた所定の係数である。また、プロセッサー110は、時刻Tbを推定する。時刻Tbは、例えば、
Tb=T(i-1)+Tb(i-1) (21)
で求めることができる。
Alternatively, the processor 110 may obtain time Tx(i-1) instead of time Txi, and time Tb(i-1) instead of time Tbi. Time Tx(i-1) is the difference between time Ti at the (i-1)th point Pd and time TQ at point Q. Time Tb(i-1) is the time from time T(i-1) to time Tb at the (i-1)th point Pd. Time Tx(i-1) is
Tx(i-1)=ΔT(i-1)-(R(i-1))/m (19)
The time Tb(i-1) can be calculated by, for example,
Tb(i-1)=Tx(i-1)×γ2 (20)
Here, the coefficient γ2 is a predetermined coefficient. The processor 110 also estimates the time Tb. The time Tb can be calculated, for example, as follows:
Tb=T(i-1)+Tb(i-1) (21)
It can be found by:
また、時間ΔTbは、時刻Tbから現在時刻を引くことで求めることができる。
なお、時刻Tb及び時間ΔTbを求める各式は上に示したものに限定しない。
Furthermore, the time ΔTb can be obtained by subtracting the current time from the time Tb.
The equations for determining the time Tb and the time ΔTb are not limited to those shown above.
また、プロセッサー110は、運行が正常になるまでの時間ΔTb及び時刻Tbを、例えば、ステップST29と同様に各種の装置に送信する。各種の装置は、例えば、ディスプレイに表示するなどして、運行が正常になるまでの時間ΔTb及び時刻Tbを報知する。各種の装置は、例えば、ユーザーが使用する端末などに運行が正常になるまでの時間ΔTb及び時刻Tbを送信する。なお、プロセッサー110は、ステップST51とステップST55において共通の閾値CDaを用いている。しかしながら、プロセッサー110は、ステップST51とステップST55で異なる閾値を用いても良い。 In addition, the processor 110 transmits the time ΔTb and time Tb until operation returns to normal to various devices, for example, as in step ST29. The various devices notify the time ΔTb and time Tb until operation returns to normal, for example, by displaying them on a display. The various devices transmit the time ΔTb and time Tb until operation returns to normal to, for example, a terminal used by a user. Note that the processor 110 uses a common threshold CDa in steps ST51 and ST55. However, the processor 110 may use different thresholds in steps ST51 and ST55.
ステップST56においてプロセッサー110は、運行が正常になったか否かを判定する。プロセッサー110は、例えば、混雑率が正常化し始めてから所定の時間経過したならば、運行が正常になったと判定する。あるいは、プロセッサー110は、運行管理システム200から運行が正常になったことを示す情報を受信したことに応じて運行が正常になったと判定する。あるいは、プロセッサー110は、計測装置300で計測された混雑率が所定の正常範囲内となって事に応じて運行が正常になったと判定する。プロセッサー110は、運行が正常になったと判定しないならば、ステップST56においてNoと判定してステップST57へと進む。対して、プロセッサー110は、運行が正常になったと判定するならば、ステップST56においてYesと判定してステップST57へと進む。 In step ST56, processor 110 determines whether operation has returned to normal. For example, processor 110 determines that operation has returned to normal if a predetermined time has passed since the congestion rate began to return to normal. Alternatively, processor 110 determines that operation has returned to normal in response to receiving information from traffic management system 200 indicating that operation has returned to normal. Alternatively, processor 110 determines that operation has returned to normal in response to the congestion rate measured by measuring device 300 falling within a predetermined normal range. If processor 110 does not determine that operation has returned to normal, it determines No in step ST56 and proceeds to step ST57. On the other hand, if processor 110 determines that operation has returned to normal, it determines Yes in step ST56 and proceeds to step ST57.
ステップST57においてプロセッサー110は、時系列モデル及び変動表を正常化時用のものから通常時用のものに変更する。すなわち、プロセッサー110は、図11のステップST26において求めた混雑率をそのまま用いる。プロセッサー110は、図12のステップST57の処理の後、ステップST51へと戻る。 In step ST57, processor 110 changes the time series model and fluctuation table from those for normalization to those for normal times. That is, processor 110 uses the congestion rate calculated in step ST26 of FIG. 11 as is. After processing step ST57 of FIG. 12, processor 110 returns to step ST51.
第4実施形態の混雑予測システム1bは、第3実施形態と同様の効果が得られる。 The congestion prediction system 1b of the fourth embodiment achieves the same effects as the third embodiment.
また、第4実施形態の混雑予測システム1bによれば、混雑予測サーバー100は、運行に異常が発生した場合、異常時用の混雑率の予測方法を用いて混雑率を予測する。これにより、混雑予測サーバー100は、運行に異常が発生している場合の混雑率をより精度良く求めることができる。また、混雑予測サーバー100は、ダイヤ乱れ等により混雑状態が過渡的に異常な状態になった場合に、異常な混雑状態から平常状態に収束する時間の予測ができる。 Furthermore, according to the congestion prediction system 1b of the fourth embodiment, when an abnormality occurs in operation, the congestion prediction server 100 predicts the congestion rate using a congestion rate prediction method for abnormal situations. This allows the congestion prediction server 100 to more accurately determine the congestion rate when an abnormality occurs in operation. Furthermore, when the congestion state becomes transiently abnormal due to a train schedule disruption or the like, the congestion prediction server 100 can predict the time it will take for the abnormal congestion state to converge to a normal state.
上記の実施形態は、以下のような変形も可能である。
上記の実施形態では、プロセッサー110は、ステップST24~ステップST26において差分を用いて混雑率を予測した。しかしながら、プロセッサー110は、差分に代えて比を用いて混雑率を予測しても良い。
The above embodiment can be modified as follows.
In the above embodiment, the processor 110 predicts the congestion rate using the difference in steps ST24 to ST26. However, the processor 110 may predict the congestion rate using the ratio instead of the difference.
上記の実施形態においては混雑率管理システム400が実行する処理の一部又は全部を混雑予測サーバー100又は運行管理システム200が実行しても良い。上記の実施形態において運行管理システム200が実行する処理の一部又は全部を混雑予測サーバー100又は混雑率管理システム400が実行しても良い。上記の実施形態において混雑予測サーバー100が実行する処理の一部又は全部を運行管理システム200又は混雑率管理システム400が実行しても良い。 In the above embodiment, some or all of the processing performed by the congestion rate management system 400 may be performed by the congestion prediction server 100 or the traffic management system 200. In the above embodiment, some or all of the processing performed by the traffic management system 200 may be performed by the congestion prediction server 100 or the congestion rate management system 400. In the above embodiment, some or all of the processing performed by the congestion prediction server 100 may be performed by the traffic management system 200 or the congestion rate management system 400.
上記の実施形態において混雑率管理システム400が記憶する情報の一部又は全部を混雑予測サーバー100又は運行管理システム200が記憶しても良い。上記の実施形態において運行管理システム200が記憶する情報の一部又は全部を混雑予測サーバー100又は混雑率管理システム400が記憶しても良い。上記の実施形態において混雑予測サーバー100が記憶する情報の一部又は全部を運行管理システム200又は混雑率管理システム400が記憶しても良い。 In the above embodiments, some or all of the information stored by the congestion rate management system 400 may be stored by the congestion prediction server 100 or the traffic management system 200. In the above embodiments, some or all of the information stored by the traffic management system 200 may be stored by the congestion prediction server 100 or the congestion rate management system 400. In the above embodiments, some or all of the information stored by the congestion prediction server 100 may be stored by the traffic management system 200 or the congestion rate management system 400.
混雑予測システム1は、混雑率に代えて、混雑率以外の乗り物の混雑の状態を示す値を用いても良い。 Instead of the congestion rate, the congestion prediction system 1 may use a value other than the congestion rate that indicates the state of vehicle congestion.
混雑予測システム1は、平均値に代えて中央値などを用いても良い。 The congestion prediction system 1 may use a median instead of the average value.
上記の実施形態の混雑予測システム1は、鉄道の混雑率を予測するシステムである。しかしながら、実施形態の混雑予測システムは、鉄道以外の乗り物の混雑率を予測するシステムであっても良い。例えば、実施形態の混雑予測システムは、路線バス、トロリーバス、ガイドウェイバス若しくはその他のバス、船舶、航空機、エレベーター又はその他の乗り物の混雑率を予測する。 The congestion prediction system 1 of the above embodiment is a system that predicts the congestion rate of a railway. However, the congestion prediction system of the embodiment may also be a system that predicts the congestion rate of a vehicle other than a railway. For example, the congestion prediction system of the embodiment predicts the congestion rate of a route bus, trolley bus, guideway bus, or other bus, a ship, an airplane, an elevator, or other vehicle.
混雑予測サーバー100は、複数の装置に分かれていても良い。そして、当該複数の装置は、上記の実施形態における処理を分担して処理しても良い。 The congestion prediction server 100 may be divided into multiple devices. These multiple devices may then share the processing responsibilities of the above embodiments.
プロセッサー110は、上記実施形態においてプログラムによって実現する処理の一部又は全部を、回路のハードウェア構成によって実現するものであっても良い。 The processor 110 may implement some or all of the processing implemented by a program in the above embodiments using a hardware circuit configuration.
実施形態の処理を実現するプログラムは、例えば装置に記憶された状態で譲渡される。しかしながら、当該装置は、当該プログラムが記憶されない状態で譲渡されても良い。そして、当該プログラムが別途に譲渡され、当該装置へと書き込まれても良い。このときのプログラムの譲渡は、例えば、リムーバブルな記憶媒体に記録して、あるいはインターネット又はLANなどのネットワークを介したダウンロードによって実現できる。 The program that realizes the processing of the embodiment is transferred, for example, in a state stored in the device. However, the device may also be transferred without the program stored therein. The program may then be transferred separately and written to the device. In this case, the program may be transferred, for example, by recording it on a removable storage medium or by downloading it via a network such as the Internet or a LAN.
以上、本発明の実施形態を説明したが、例として示したものであり、本発明の範囲を限定するものではない。本発明の実施形態は、本発明の要旨を逸脱しない範囲において種々の態様で実施可能である。 The above describes an embodiment of the present invention, but it is presented as an example and does not limit the scope of the present invention. The embodiment of the present invention can be implemented in various forms without departing from the spirit of the present invention.
1,1b 混雑予測システム
100 混雑予測サーバー
110 プロセッサー
120 ROM
130 RAM
140 補助記憶装置
141 混雑データベース
150 通信I/F
160 バス
200 運行管理システム
300 計測装置
400 混雑率管理システム
500 改札システム
1, 1b Congestion prediction system 100 Congestion prediction server 110 Processor 120 ROM
130 RAM
140 Auxiliary storage device 141 Congestion database 150 Communication I/F
160 Bus 200 Operation control system 300 Measuring device 400 Congestion rate control system 500 Ticket gate system
Claims (9)
前記処理部は、
前記統計データと前記第1の計測値の差分を用いて、前記第1の区間での前記予測対象の乗り物の混雑の状態の時間変化を予測する関数を求め、前記関数を用いて前記第1の区間での前記予測対象の乗り物の混雑の状態を予測し、
前記第1の区間とは異なる第2の区間での前記予測対象の乗り物の混雑の状態を、前記予測対象の乗り物の前記第1の区間での混雑の状態を用い、前記第1の区間と前記第2の区間の相関関係として区間ごとの前記乗り物の混雑の状態の変化を示す変動情報を用いて予測し、
前記変動情報は、前記予測対象の乗り物としての列車の各車両について、駅区間ごとの号車別混雑率の変動を示す変動表である、
混雑率予測装置。 a processing unit that predicts a state of vehicle congestion in a first section that is a prediction target, using statistical data on a state of vehicle congestion in a first section that is composed of one or more sections and a first measurement value of the state of vehicle congestion in the first section ;
The processing unit
using a difference between the statistical data and the first measurement value to obtain a function for predicting a time change in the congestion state of the vehicle to be predicted in the first section, and predicting the congestion state of the vehicle to be predicted in the first section using the function;
predicting a congestion state of the vehicle to be predicted in a second section different from the first section using the congestion state of the vehicle to be predicted in the first section and using fluctuation information indicating a change in the congestion state of the vehicle for each section as a correlation between the first section and the second section;
The fluctuation information is a fluctuation table showing fluctuations in congestion rate by car for each station section for each car of the train as the vehicle to be predicted.
Congestion rate prediction device.
前記変動表の各駅区間の混雑率と実測混雑率との差又は比を求め、当該差又は比を用いて予測対象の列車がまだ到達していない駅区間それぞれについての差又は比を予測し、Calculating the difference or ratio between the congestion rate of each station section in the fluctuation table and the actually measured congestion rate, and predicting the difference or ratio for each station section that the train to be predicted has not yet arrived using the difference or ratio;
予測した差又は比を用いて前記変動表を補正する、請求項1又は2に記載の混雑率予測装置。The congestion rate prediction device according to claim 1 or 2, wherein the variation table is corrected using the predicted difference or ratio.
1又は複数の区間により構成される第1の区間での乗り物の混雑の状態の統計データと、前記第1の区間での前記乗り物の混雑の状態の第1の計測値を用いて、前記第1の区間での予測対象の乗り物の混雑の状態を予測する処理部として機能させるコンピュータプログラムであって、
前記処理部は、
前記統計データと前記第1の計測値の差分を用いて、前記第1の区間での前記予測対象の乗り物の混雑の状態の時間変化を予測する関数を求め、前記関数を用いて前記第1の区間での前記予測対象の乗り物の混雑の状態を予測し、
前記第1の区間とは異なる第2の区間での前記予測対象の乗り物の混雑の状態を、前記予測対象の乗り物の前記第1の区間での混雑の状態を用い、前記第1の区間と前記第2の区間の相関関係として区間ごとの前記乗り物の混雑の状態の変化を示す変動情報を用いて予測し、
前記変動情報は、前記予測対象の乗り物としての列車の各車両について、駅区間ごとの号車別混雑率の変動を示す変動表である、
コンピュータプログラム。 The processor in the congestion rate prediction device is
A computer program that functions as a processing unit that predicts a state of vehicle congestion in a first section that is a target of prediction, using statistical data on a state of vehicle congestion in a first section that is composed of one or more sections and a first measurement value of the state of vehicle congestion in the first section, the computer program comprising:
The processing unit
using a difference between the statistical data and the first measurement value to obtain a function for predicting a time change in the congestion state of the vehicle to be predicted in the first section, and predicting the congestion state of the vehicle to be predicted in the first section using the function;
predicting a congestion state of the vehicle to be predicted in a second section different from the first section using the congestion state of the vehicle to be predicted in the first section and using fluctuation information indicating a change in the congestion state of the vehicle for each section as a correlation between the first section and the second section;
The fluctuation information is a fluctuation table showing fluctuations in congestion rate by car for each station section for each car of the train as the vehicle to be predicted.
Computer program .
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2021113717A JP7719645B2 (en) | 2021-07-08 | 2021-07-08 | Congestion rate prediction device and program |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2021113717A JP7719645B2 (en) | 2021-07-08 | 2021-07-08 | Congestion rate prediction device and program |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2023009987A JP2023009987A (en) | 2023-01-20 |
| JP7719645B2 true JP7719645B2 (en) | 2025-08-06 |
Family
ID=85119156
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2021113717A Active JP7719645B2 (en) | 2021-07-08 | 2021-07-08 | Congestion rate prediction device and program |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP7719645B2 (en) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2015009604A (en) | 2013-06-27 | 2015-01-19 | 株式会社日立製作所 | Train congestion prediction system and train congestion prediction method |
| JP2016168876A (en) | 2015-03-11 | 2016-09-23 | 株式会社東芝 | Congestion prediction apparatus and congestion prediction method |
| JP2017194863A (en) | 2016-04-21 | 2017-10-26 | 株式会社Nttドコモ | Demand forecasting device |
| JP2018140683A (en) | 2017-02-27 | 2018-09-13 | 三菱電機株式会社 | Train operation management device, train operation management method, and train operation management program |
| JP2020177354A (en) | 2019-04-16 | 2020-10-29 | 大日本印刷株式会社 | Predictor, prediction system, computer program, prediction method and generation method |
| JP2021049813A (en) | 2019-09-24 | 2021-04-01 | 株式会社日立製作所 | Train boarding rate management system and train boarding rate management method |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP3539844B1 (en) * | 2016-11-08 | 2025-10-08 | Hitachi, Ltd. | Transportation system, schedule proposal system, and train operations system |
-
2021
- 2021-07-08 JP JP2021113717A patent/JP7719645B2/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2015009604A (en) | 2013-06-27 | 2015-01-19 | 株式会社日立製作所 | Train congestion prediction system and train congestion prediction method |
| JP2016168876A (en) | 2015-03-11 | 2016-09-23 | 株式会社東芝 | Congestion prediction apparatus and congestion prediction method |
| JP2017194863A (en) | 2016-04-21 | 2017-10-26 | 株式会社Nttドコモ | Demand forecasting device |
| JP2018140683A (en) | 2017-02-27 | 2018-09-13 | 三菱電機株式会社 | Train operation management device, train operation management method, and train operation management program |
| JP2020177354A (en) | 2019-04-16 | 2020-10-29 | 大日本印刷株式会社 | Predictor, prediction system, computer program, prediction method and generation method |
| JP2021049813A (en) | 2019-09-24 | 2021-04-01 | 株式会社日立製作所 | Train boarding rate management system and train boarding rate management method |
Also Published As
| Publication number | Publication date |
|---|---|
| JP2023009987A (en) | 2023-01-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109789885B (en) | Traffic system, schedule suggestion system and vehicle operation system | |
| JP7399627B2 (en) | Diagram creation device, diagram creation method, and automobile control system | |
| EP3437957B1 (en) | Train traffic control system and train traffic control method | |
| Soza-Parra et al. | Public transport reliability across preferences, modes, and space | |
| Zhang et al. | Bus dwell time estimation and prediction: A study case in Shanghai-China | |
| El-Geneidy et al. | Limited-stop bus service: an evaluation of an implementation strategy | |
| JP6068634B2 (en) | Route operation result prediction time generation system, route operation result prediction time generation method, and route operation result prediction time generation program | |
| JP2014091480A (en) | Operation arrangement device and method | |
| JP6454222B2 (en) | Data processing system and data processing method | |
| JP2017091008A (en) | Data processing method and data processing system | |
| JP7158288B2 (en) | Traffic state estimation system and traffic state estimation method | |
| JP2018140683A (en) | Train operation management device, train operation management method, and train operation management program | |
| JP2021098425A (en) | Operation forecasting system and method for the same | |
| JP5969957B2 (en) | Train congestion degree prediction system and train congestion degree prediction method | |
| JP6393531B2 (en) | Train selection support system, train selection support method, and train selection support program | |
| US12272241B2 (en) | Artificial intelligence for responsive operation for vehicle fleet management and method thereof | |
| JP7719645B2 (en) | Congestion rate prediction device and program | |
| Baali et al. | Long enough but not too long: a posteriori determination of the dwell time margins from high-resolution passenger flow data | |
| JP7425680B2 (en) | Navigation device and navigation method | |
| JP7358767B2 (en) | Information processing device, notification method and program | |
| JP2019160119A (en) | Transportation transit determination device and transportation transit determination method | |
| JP2022502777A (en) | How to tune a multimodal transportation network | |
| JP7787725B2 (en) | Seating information processing device, seating information providing system, route search system, program, and seating information processing method | |
| Kunimatsu et al. | The prediction of passenger flow under transport disturbance using accumulated passenger data | |
| Ponkratov et al. | Evaluation of overcrowded mass transit demand-supply interaction considering the fail to board probability and waiting time increase |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20240424 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20250129 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20250204 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20250403 |
|
| 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: 20250708 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20250725 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 7719645 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |