JP2013210941A - Traffic flow prediction device, traffic flow prediction method, and traffic flow prediction program - Google Patents

Traffic flow prediction device, traffic flow prediction method, and traffic flow prediction program Download PDF

Info

Publication number
JP2013210941A
JP2013210941A JP2012081996A JP2012081996A JP2013210941A JP 2013210941 A JP2013210941 A JP 2013210941A JP 2012081996 A JP2012081996 A JP 2012081996A JP 2012081996 A JP2012081996 A JP 2012081996A JP 2013210941 A JP2013210941 A JP 2013210941A
Authority
JP
Japan
Prior art keywords
traffic
approximate curve
traffic volume
relationship
volume
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP2012081996A
Other languages
Japanese (ja)
Other versions
JP5898553B2 (en
Inventor
Ryota Horiguchi
良太 堀口
Hisatomo Hanafusa
比佐友 花房
Katsuaki Koide
勝亮 小出
Masato Kobayashi
正人 小林
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
I-TRANSPORT LAB CO Ltd
Original Assignee
I-TRANSPORT LAB CO Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by I-TRANSPORT LAB CO Ltd filed Critical I-TRANSPORT LAB CO Ltd
Priority to JP2012081996A priority Critical patent/JP5898553B2/en
Publication of JP2013210941A publication Critical patent/JP2013210941A/en
Application granted granted Critical
Publication of JP5898553B2 publication Critical patent/JP5898553B2/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Traffic Control Systems (AREA)

Abstract

PROBLEM TO BE SOLVED: To provide a traffic flow prediction device capable of rapidly and highly accurately predicting a traffic flow of vehicles.SOLUTION: A traffic flow prediction device includes: means for obtaining a first approximation curve showing a relationship between traffic density and traffic volume within a zone on the basis of probe traffic information of vehicles travelling in the zone; means for obtaining a second approximation curve on the basis of a result of performing traffic flow simulation; means for computing a magnification factor for matching the first approximation curve to the second approximation curve; means for obtaining traffic volume and link capacity minimizing a distance between a value indicating a relationship between traffic density and traffic volume within a zone that is a product of multiplying a respective value indicating the relationship between the traffic density and traffic volume within the zone when obtaining the first approximation curve, by the magnification factor, and a value indicating a relationship between traffic density and traffic volume within the zone when obtaining the second approximation curve; and means that inputs the obtained traffic volume and link capacity and that repeats processing until a result of executing the traffic flow simulation satisfies a prescribed termination condition.

Description

本発明は、車両の交通流を予測する交通流予測装置、交通流予測方法及び交通流予測プログラムに関する。   The present invention relates to a traffic flow prediction device, a traffic flow prediction method, and a traffic flow prediction program that predict a traffic flow of a vehicle.

現在から比較的長い将来にわたった交通状況を高い精度で予測することができる交通状況推定装置が知られている(例えば、特許文献1参照)。この交通状況推定装置は、交通状況計算部において、道路ネットワークデータ、道路ネットワーク上の信号データ、及びOD交通量を用いた交通流シミュレーションを行なって、交通量の計算値を計算し、将来予測処理部よって、過去の交通量の観測値と現在の交通量の観測値とに基づいて、将来の交通量の予測値を予測し、OD交通量推定部で計算値が現在の交通量の観測値及び予測値に一致するように交通状況計算部で使用するOD交通量を修正するものである。この交通状況推定装置によれば、交通流シミュレーションを実行することにより、交通状況を示す交通量を予測する際に、交通量の計算値が現在の交通量の観測値及び交通量の予測値に一致するように修正されたOD交通量を用いているので交通量を高い精度で予測することができる。   2. Description of the Related Art A traffic situation estimation device that can predict a traffic situation from the present to a relatively long future with high accuracy is known (see, for example, Patent Document 1). In this traffic situation estimation device, the traffic situation calculation unit performs traffic flow simulation using road network data, signal data on the road network, and OD traffic volume, calculates a traffic volume calculation value, and performs future prediction processing. Section predicts the predicted value of future traffic based on the observed value of the past traffic and the observed value of the current traffic, and the calculated value is the observed value of the current traffic by the OD traffic estimation unit. The OD traffic volume used in the traffic situation calculation unit is corrected so as to coincide with the predicted value. According to this traffic situation estimation apparatus, when a traffic volume indicating a traffic situation is predicted by executing a traffic flow simulation, the calculated traffic volume is changed to the current observed traffic volume and the predicted traffic volume. Since the OD traffic volume corrected so as to match is used, the traffic volume can be predicted with high accuracy.

特開2004−118735号公報JP 2004-118735 A

ところで、これまでの交通流シミュレーションは、主要幹線道路などの路線に沿って固定的に配置されるトラフィックカウンターなどによって検出された交通情報データを用いて、トラフィックカウンターが配置された地点の車両の通行量と速度に基づいて交通流シミュレーションを実行し、交通量を予測していた。近年、走行する車両を移動するセンサーに見立て、実際に走行している車両の走行状態を検出して提供情報として用いることが検討されている。このような車両を検出することにより得られた情報は、「プローブ交通情報」といわれ、利用範囲を拡大することが期待されている。   By the way, traffic flow simulations so far have been based on traffic information data detected by traffic counters, etc., which are fixedly arranged along main trunk roads, etc. Based on the volume and speed, a traffic flow simulation was performed to predict the traffic volume. 2. Description of the Related Art In recent years, it has been studied to detect a traveling state of a vehicle that is actually traveling and use it as provided information, considering the traveling vehicle as a moving sensor. Information obtained by detecting such a vehicle is called “probe traffic information” and is expected to expand the range of use.

しかしながら、プローブ交通情報は道路上を走行する全ての車両のうち一部を構成するのみであるため、プローブ交通情報から車両の通行量を直接計測することができず、通行量を調整しながら交通流予測を行うことができないという問題がある。また、従来の交通流シミュレーションで実際の交通状態に則した交通流予測を行うためには、試行錯誤を繰り返しながらシミュレーション実行を何度も実行しなければならないため、予測処理に多大な時間を要するという問題がある。   However, since the probe traffic information only constitutes a part of all the vehicles traveling on the road, the traffic volume of the vehicle cannot be directly measured from the probe traffic information, and the traffic is adjusted while adjusting the traffic volume. There is a problem that flow prediction cannot be performed. In addition, in order to perform traffic flow prediction according to the actual traffic state in the conventional traffic flow simulation, it is necessary to execute the simulation many times while repeating trial and error. There is a problem.

本発明は、このような事情に鑑みてなされたもので、車両の交通流を迅速にかつ高い精度で予測することができる交通流予測装置、交通流予測方法及び交通流予測プログラムを提供することを目的とする。   The present invention has been made in view of such circumstances, and provides a traffic flow prediction device, a traffic flow prediction method, and a traffic flow prediction program capable of predicting a traffic flow of a vehicle quickly and with high accuracy. With the goal.

本発明は、観測値である時刻、緯度経度情報からなるプローブ情報、もしくは時刻、区間番号、区間通過台数、及び区間通過平均所要時間からなるプローブ交通情報から、任意に指定した区域の交通密度と交通量の関係を定量化し、該関係を用いて交通量を推定する交通流予測装置であって、前記区域を走行した車両のプローブ交通情報に基づき、前記区域内の交通密度と交通量の関係を示す第1の近似曲線を求める第1の近似曲線算出手段と、交通流シミュレーションモデルを用いて交通流シミュレーションを行った結果に基づき、前記区域内の交通密度と交通量の関係を示す第2の近似曲線を求める第2の近似曲線算出手段と、前記第1の近似曲線を前記第2の近似曲線に合致させるための拡大率を算定する拡大率算定手段と、前記第1の近似曲線を求めた際の前記区域内の交通密度と交通量の関係を示す値のそれぞれに対して、前記拡大率を乗算して得られた前記区域内の交通密度と交通量の関係を示す値が、前記第2の近似曲線を求めた際の前記区域内の交通密度と交通量の関係を示す値との距離が最小化させるための交通量とリンク容量を求めるパラメータ調整手段と、前記パラメータ調整手段によって求めた前記交通量と前記リンク容量を入力として前記交通流シミュレーションを実行した結果が所定の終了条件を満たすまでパラメータ調整手段による前記交通量と前記リンク容量を求める処理と前記交通流シミュレーションの実行を繰り返す結果判定手段とを備えたことを特徴とする。   The present invention relates to the traffic density of an arbitrarily designated area from probe information consisting of time, latitude and longitude information, or probe traffic information consisting of time, section number, number of sections passing, and section passing average required time. A traffic flow prediction apparatus that quantifies a traffic volume relationship and estimates a traffic volume using the relationship, and based on probe traffic information of a vehicle that has traveled in the zone, a relationship between traffic density and traffic volume in the zone Based on the result of the traffic flow simulation using the first approximate curve calculation means for obtaining the first approximate curve indicating the traffic flow and the traffic flow simulation model, the second shows the relationship between the traffic density in the area and the traffic volume. Second approximate curve calculating means for obtaining the approximate curve, an enlargement ratio calculating means for calculating an enlargement ratio for matching the first approximate curve with the second approximate curve, and the first The relationship between the traffic density in the area and the traffic volume obtained by multiplying the expansion rate for each value indicating the relationship between the traffic density and the traffic volume in the area when the approximate curve is obtained is shown. Parameter adjusting means for determining a traffic volume and a link capacity for minimizing a distance between a value indicating a relationship between the traffic density in the area and the traffic volume when the second approximate curve is obtained; Processing for obtaining the traffic volume and the link capacity by the parameter adjustment means until the result of executing the traffic flow simulation with the traffic volume and the link capacity obtained by the parameter adjustment means as input satisfies a predetermined termination condition, and the traffic flow And a result determination means for repeating the execution of the simulation.

本発明は、観測値である時刻、緯度経度情報からなるプローブ情報、もしくは時刻、区間番号、区間通過台数、及び区間通過平均所要時間からなるプローブ交通情報から、任意に指定した区域の交通密度と交通量の関係を定量化し、該関係を用いて交通量を推定する交通流予測装置が行う交通流予測方法であって、前記区域を走行した車両のプローブ交通情報に基づき、前記区域内の交通密度と交通量の関係を示す第1の近似曲線を求める第1の近似曲線算出ステップと、交通流シミュレーションモデルを用いて交通流シミュレーションを行った結果に基づき、前記区域内の交通密度と交通量の関係を示す第2の近似曲線を求める第2の近似曲線算出ステップと、前記第1の近似曲線を前記第2の近似曲線に合致させるための拡大率を算定する拡大率算定ステップと、前記第1の近似曲線を求めた際の前記区域内の交通密度と交通量の関係を示す値のそれぞれに対して、前記拡大率を乗算して得られた前記区域内の交通密度と交通量の関係を示す値が、前記第2の近似曲線を求めた際の前記区域内の交通密度と交通量の関係を示す値との距離が最小化させるための交通量とリンク容量を求めるパラメータ調整ステップと、前記パラメータ調整ステップによって求めた前記交通量と前記リンク容量を入力として前記交通流シミュレーションを実行した結果が所定の終了条件を満たすまでパラメータ調整ステップによる前記交通量と前記リンク容量を求める処理と前記交通流シミュレーションの実行を繰り返す結果判定ステップとを備えたことを特徴とする。   The present invention relates to the traffic density of an arbitrarily designated area from probe information consisting of time, latitude and longitude information, or probe traffic information consisting of time, section number, number of sections passing, and section passing average required time. A traffic flow prediction method performed by a traffic flow prediction device that quantifies a traffic volume relationship and estimates a traffic volume using the relationship, and based on probe traffic information of a vehicle that has traveled in the zone, traffic in the zone Based on the result of having performed the traffic flow simulation using the traffic flow simulation model and the first approximate curve calculating step for obtaining the first approximate curve indicating the relationship between the density and the traffic volume, the traffic density and the traffic volume in the area A second approximate curve calculating step for obtaining a second approximate curve indicating the relationship between the first approximate curve and an enlargement factor for matching the first approximate curve with the second approximate curve. In the area obtained by multiplying the enlargement ratio to the enlargement ratio calculating step and each of the values indicating the relationship between the traffic density and the traffic volume in the area when the first approximate curve is obtained. The value indicating the relationship between the traffic density and the traffic volume is a traffic volume for minimizing the distance between the traffic density in the area and the value indicating the traffic volume when the second approximate curve is obtained. A parameter adjustment step for obtaining a link capacity; and the traffic volume obtained by the parameter adjustment step until a result of executing the traffic flow simulation with the traffic volume obtained in the parameter adjustment step and the link capacity as an input satisfies a predetermined termination condition; It is characterized by comprising a process for obtaining the link capacity and a result determination step for repeating the execution of the traffic flow simulation.

本発明は、観測値である時刻、緯度経度情報からなるプローブ情報、もしくは時刻、区間番号、区間通過台数、及び区間通過平均所要時間からなるプローブ交通情報から、任意に指定した区域の交通密度と交通量の関係を定量化し、該関係を用いて交通量を推定する交通流予測装置上のコンピュータに、前記区域を走行した車両のプローブ交通情報に基づき、前記区域内の交通密度と交通量の関係を示す第1の近似曲線を求める第1の近似曲線算出ステップと、交通流シミュレーションモデルを用いて交通流シミュレーションを行った結果に基づき、前記区域内の交通密度と交通量の関係を示す第2の近似曲線を求める第2の近似曲線算出ステップと、前記第1の近似曲線を前記第2の近似曲線に合致させるための拡大率を算定する拡大率算定ステップと、前記第1の近似曲線を求めた際の前記区域内の交通密度と交通量の関係を示す値のそれぞれに対して、前記拡大率を乗算して得られた前記区域内の交通密度と交通量の関係を示す値が、前記第2の近似曲線を求めた際の前記区域内の交通密度と交通量の関係を示す値との距離が最小化させるための交通量とリンク容量を求めるパラメータ調整ステップと、前記パラメータ調整ステップによって求めた前記交通量と前記リンク容量を入力として前記交通流シミュレーションを実行した結果が所定の終了条件を満たすまでパラメータ調整ステップによる前記交通量と前記リンク容量を求める処理と前記交通流シミュレーションの実行を繰り返す結果判定ステップとを行わせることを特徴とする。   The present invention relates to the traffic density of an arbitrarily designated area from probe information consisting of time, latitude and longitude information, or probe traffic information consisting of time, section number, number of sections passing, and section passing average required time. Based on the probe traffic information of the vehicle that has traveled in the area, the computer on the traffic flow prediction device that quantifies the traffic volume relationship and estimates the traffic volume using the relationship is used to determine the traffic density and traffic volume in the area. A first approximate curve calculating step for obtaining a first approximate curve indicating the relationship and a result of performing a traffic flow simulation using the traffic flow simulation model is used to indicate a relationship between the traffic density and the traffic volume in the area. A second approximate curve calculation step for obtaining two approximate curves, and an enlargement factor calculation for calculating an enlargement factor for matching the first approximate curve with the second approximate curve The traffic density in the area obtained by multiplying the enlargement rate for each of the step and the value indicating the relationship between the traffic density in the area and the traffic volume when the first approximate curve is obtained The traffic volume and the link capacity for minimizing the distance between the traffic density in the area and the value indicating the traffic volume when the second approximate curve is obtained. A parameter adjustment step to be obtained; and the traffic flow and the link capacity obtained by the parameter adjustment step until a result of executing the traffic flow simulation by inputting the traffic and the link capacity obtained by the parameter adjustment step satisfies a predetermined end condition. And a result determination step that repeats execution of the traffic flow simulation.

本発明によれば、プローブ交通情報から得られる交通状態を用いて交通流予測を行うようにしたため、交通流シミュレーションを必要最小限利用するのみで高精度な交通流予測が可能になるという効果が得られる。   According to the present invention, since the traffic flow prediction is performed using the traffic state obtained from the probe traffic information, it is possible to perform highly accurate traffic flow prediction only by using the traffic flow simulation as much as necessary. can get.

本発明の一実施形態の構成を示すブロック図である。It is a block diagram which shows the structure of one Embodiment of this invention. 図1に示す交通流予測装置1の動作を示すフローチャートである。It is a flowchart which shows operation | movement of the traffic flow prediction apparatus 1 shown in FIG. ネットワークデータの一例を示す説明図である。It is explanatory drawing which shows an example of network data. ネットワークデータを構成するノードとリンクの一例を示す説明図である。It is explanatory drawing which shows an example of the node and link which comprise network data. シミュレーションによって求めた集計QKをプロットした結果から求めた近似曲線と、プローブ交通情報から求めた集計QKをプロットした結果から求めた近似曲線を示す図である。It is a figure which shows the approximate curve calculated | required from the result of plotting the total curve obtained from the result of plotting the total QK calculated | required by simulation, and the total QK calculated | required from probe traffic information. プローブ交通情報から求めた集計QKの値に拡大率α、βを乗算した結果をプロットした図である。It is the figure which plotted the result of multiplying the value of total QK calculated | required from the probe traffic information by expansion rate (alpha) and (beta).

以下、図面を参照して、本発明の一実施形態による交通流予測装置を説明する。図1は同実施形態の構成を示すブロック図である。この図において、符号1は、コンピュータ装置で構成する交通流予測装置である。符号2は、キーボードやポインティングデバイスで構成する入力部である。符号3は、ディスプレイ装置で構成する表示部である。符号4は、交通流を予測するのに必要なデータが記憶されたデータファイルである。   Hereinafter, a traffic flow prediction apparatus according to an embodiment of the present invention will be described with reference to the drawings. FIG. 1 is a block diagram showing the configuration of the embodiment. In this figure, the code | symbol 1 is the traffic flow prediction apparatus comprised with a computer apparatus. Reference numeral 2 denotes an input unit composed of a keyboard and a pointing device. Reference numeral 3 denotes a display unit configured by a display device. Reference numeral 4 denotes a data file in which data necessary for predicting traffic flow is stored.

符号11は、データファイル4及び入力部2から交通流を予測するのに必要なデータを入力するデータ入力部である。符号12は、データ入力部11において入力したデータに基づき交通流シミュレーションを実行するシミュレーション実行部である。ここで、実行する交通流シミュレーションは、公知のものを利用可能であり、例えば、特許文献1に記載の交通流シミュレーションや、交通流シミュレーションモデル(SOUND;Simulation on Urban road Network with Dynamic route choice;東京大学生産技術研究所において開発された広域道路ネットワークシミュレーションモデルである)が適用可能である。   Reference numeral 11 denotes a data input unit that inputs data necessary for predicting traffic flow from the data file 4 and the input unit 2. Reference numeral 12 denotes a simulation execution unit that executes a traffic flow simulation based on the data input in the data input unit 11. Here, as the traffic flow simulation to be executed, known ones can be used. For example, the traffic flow simulation described in Patent Document 1 or a traffic flow simulation model (SOUND; Simulation on Urban road Network with Dynamic route choice; Tokyo) A wide-area road network simulation model developed at the university production technology research institute is applicable.

符号13は、データ入力部11において入力したプローブ交通情報に基づき集計QK(Kはエリア存在台数[台]、Qは集計交通量[台km/時]であり、このQKの関係を示すデータ)を算出する集計QK算出部である。符号14は、シミュレーション実行部12において算出した集計QKと、集計QK算出部13において算出した集計QKとを合致させるための拡大率を算定する拡大率算定部である。符号15は、シミュレーション実行部12における実行結果が所定の条件を満たしたか否かを判定する結果判定部である。符号16は、結果判定部15において、所定の条件を満たさなかった場合、パラメータを調整し、再度シミュレーションを実行させるパラメータ調整部である。   Reference numeral 13 denotes an aggregated QK based on the probe traffic information input in the data input unit 11 (K is the number of existing areas [units], Q is the aggregated traffic volume [units km / hour], and data indicating the relationship of this QK) This is a total QK calculation unit for calculating. Reference numeral 14 denotes an enlargement rate calculation unit that calculates an enlargement rate for matching the total QK calculated by the simulation execution unit 12 with the total QK calculated by the total QK calculation unit 13. Reference numeral 15 denotes a result determination unit that determines whether the execution result in the simulation execution unit 12 satisfies a predetermined condition. Reference numeral 16 denotes a parameter adjustment unit that adjusts parameters and executes simulation again when the result determination unit 15 does not satisfy a predetermined condition.

ここで、図3、図4を参照して、本明細書内で用いる用語について説明する。図3は、ネットワークデータの一例を示す説明図である。図3において、破線で示す対象エリア(例えば、東京都)のゾーンAからゾーンBへ移動する際に主な走行経路として、経路1、経路2、経路3があることを示している。ここで、ゾーンとは、任意に決められた区域のことであり、例えば、市町村等の行政区域に相当する。   Here, the terms used in this specification will be described with reference to FIGS. FIG. 3 is an explanatory diagram showing an example of network data. FIG. 3 shows that there are route 1, route 2, and route 3 as main travel routes when moving from zone A to zone B in a target area (for example, Tokyo) indicated by a broken line. Here, the zone is an arbitrarily determined area and corresponds to an administrative area such as a municipality.

図4は、ネットワークデータを構成するノードとリンクの一例を示す図である。図4に示すように、ネットワークデータは、ノードとリンクによって構成する。ノード(図4の○印)とは、道路の交差点のことであり、リンク(図4の矢印)とは、端部がノードで車両の進行方向が決められている道路である。ノードとリンクには、それぞれを識別する識別情報と、座標情報とからなる。ネットワークデータは、図3の破線で示す対象エリア内の全てのノードとリンクのデータによって構成するものである。   FIG. 4 is a diagram illustrating an example of nodes and links constituting network data. As shown in FIG. 4, network data is configured by nodes and links. Nodes (circles in FIG. 4) are intersections of roads, and links (arrows in FIG. 4) are roads whose ends are nodes and the traveling direction of the vehicle is determined. Each node and link includes identification information for identifying each node and coordinate information. The network data is constituted by data of all nodes and links in the target area indicated by a broken line in FIG.

図1に示す交通流予測装置1は、観測値である時刻、緯度経度情報からなるプローブ交通情報から、任意に指定した集計QK区域(集計QKゾーン)の交通密度と交通量の関係(集計QK)を定量化し、その関係を用いて時間帯別車種別OD交通量を推定するものである。   The traffic flow predicting apparatus 1 shown in FIG. 1 has a relationship between traffic density and traffic volume (total QK) in an arbitrarily specified total QK area (total QK zone) based on probe traffic information including time and latitude / longitude information as observed values. ) Is quantified, and the OD traffic volume by vehicle type is estimated by using the relationship.

次に、交通流予測装置1が時間帯別車種別OD交通量を推定する処理の概要を説明する。基本の交通状態を示すデータとして、ノードとリンクから作られるネットワークデータ、車両が出発地、目的地となるODゾーンデータ、初期値としての時間帯別車種別OD交通量を入力条件とする。Oは出発地ゾーン、Dは目的地ゾーンを示す。また、予め推定対象エリアに任意の集計QKゾーンを展開する。   Next, the outline | summary of the process which the traffic flow prediction apparatus 1 estimates the vehicle classification OD traffic amount according to time zone is demonstrated. As data indicating the basic traffic state, network data created from nodes and links, OD zone data from which a vehicle is a departure point and destination, and OD traffic volume by time zone as an initial value are set as input conditions. O represents a departure zone and D represents a destination zone. In addition, an arbitrary aggregate QK zone is developed in advance in the estimation target area.

まず、集計QKゾーン内を走行したプローブ交通情報によって集計QKを計算し、横軸をK、縦軸をQとした集計QKグラフ上に展開する。その後、各集計QKから、観測値による集計QKの近似曲線を時間帯別に計算するこのとき、Kはエリア存在台数とし、集計QKゾーン内の1時間あたりのエリア存在台数リンク平均旅行時間とリンク長の積の総和から求める。Qは集計交通を示し、集計QKゾーン内の1時間あたりのリンク通過交通量とリンク長の積の総和で求める。   First, the total QK is calculated from the probe traffic information traveling in the total QK zone, and is developed on the total QK graph with the horizontal axis representing K and the vertical axis representing Q. Thereafter, an approximate curve of the total QK based on the observed values is calculated for each time zone from each total QK. At this time, K is the number of existing areas, and the average number of area existing links per hour in the total QK zone and the link length. It is calculated from the sum of products. Q indicates total traffic, and is obtained as the sum of products of link traffic volume per hour and link length in the total QK zone.

一方、交通流シミュレーションモデルを用いて対象時間帯の交通流シミュレーションを行う。その結果を用いて、あるリンクの各車両の車種、出発地、目的地から、リンクあたりの時間帯別車種別OD選択確率を算出する。同時に、時間帯別車種別OD交通量を増減させた場合の交通流シミュレーションを複数ケース行い、各集計QKゾーンにおいて計算される集計QKから、シミュレーションによる集計QKの近似曲線を時間帯別に求める。さらに、各集計QKゾーンについて、実測値による集計QKの近似曲線を、シミュレーションによる集計QKの近似曲線に合致させるための時間帯別拡大率を求める。このとき、各集計QKゾーンについて、時間帯別拡大率によって拡大された実測値による集計QKと、シミュレーションによる集計QKの位置関係が最小化するときの時間帯別車種別OD交通量、及びリンクパラメータ値(リンク容量)を得る。   On the other hand, the traffic flow simulation of the target time zone is performed using the traffic flow simulation model. Using the result, the vehicle type OD selection probability per time zone per link is calculated from the vehicle type, departure point, and destination of each vehicle on a certain link. At the same time, a plurality of traffic flow simulations are performed when the OD traffic volume by vehicle type is increased or decreased by time zone, and an approximate curve of the aggregate QK by simulation is obtained for each time zone from the aggregate QK calculated in each aggregate QK zone. Further, for each total QK zone, an enlargement factor for each time zone for matching the approximate curve of the total QK based on the actual measurement value with the approximate curve of the total QK based on the simulation is obtained. At this time, for each aggregated QK zone, the aggregated QK based on the actual measurement value expanded by the expansion rate for each time zone and the OD traffic volume by vehicle type for each time zone when the positional relationship between the aggregated QK by simulation is minimized, and the link parameter Get the value (link capacity).

次に、交通流シミュレーションモデルについて説明する。ここでは、前述のSOUNDを例にして説明する。SOUNDは、以下の特徴を有している。
(1)渋滞現象のダイナミズムを考慮しており、過飽和の交通状況を再現できる。
(2)静・的/動的な経路選択モデルを内包しており、ITS(高度道路交通システム)における情報提供や動的経路誘導などの運用策の評価が可能である。
(3)車種などの各種の属性を付与した個別の車両を扱うので、対象車両を限定した交通運用施策の評価が可能である。
(4)リンク毎に与えた交通量一密度(Q−K)特性を用いて車両移動の計算をするマクロなモデルであり、計算負荷が小さいため、大規模なネットワークに適用可能である。
Next, a traffic flow simulation model will be described. Here, the above-described SOUND will be described as an example. SOUND has the following characteristics.
(1) Considering the dynamism of the traffic jam phenomenon, it is possible to reproduce oversaturated traffic conditions.
(2) It incorporates a static / dynamic / dynamic route selection model, and it is possible to evaluate operational measures such as information provision and dynamic route guidance in ITS (Intelligent Transport System).
(3) Since individual vehicles with various attributes such as vehicle types are handled, it is possible to evaluate traffic operation measures with limited target vehicles.
(4) It is a macro model that calculates vehicle movement using the traffic density density (QK) characteristic given to each link. Since the calculation load is small, it can be applied to a large-scale network.

シミュレーションでは、各種の施策や交通運用策を適用した場合の、リンク交通量や速度、路線別平均旅行時間、総走行距離などの走行効率に関する指標を予測できる。対象とする交通運用策は次の通りである。
(a)道路網整備。
(b)ロードプライシングや需要分散、流入規制などのTDM政策。
(c)右折禁止や一方通行などの進行方向に関する時間帯別交通規制。
(d)信号制御パラメータの変更。
(e)イベントや災害時における一時的な通行止め、あるいは流入制限。
(f)情報提供による経路変更の促進や迂回誘導。
In the simulation, it is possible to predict indexes related to travel efficiency, such as link traffic volume and speed, average travel time by route, and total travel distance when various measures and traffic operation measures are applied. The target traffic management measures are as follows.
(A) Road network development.
(B) TDM policies such as road pricing, demand distribution and inflow regulation.
(C) Traffic regulations by time zone regarding the direction of travel such as prohibition of right turn and one-way traffic
(D) Change of signal control parameters.
(E) Temporary closure or inflow restrictions during events or disasters.
(F) Promotion of route changes and guidance for detours by providing information.

次に、シミュレーションモデルの入出力について説明する。始めに、シミュレーションエンジンについて説明する。シミュレーションエンジンは、道路ネットワークの接続情報と幾何形状、任意の分間隔で指定した時間帯別OD交通量、信号制御パラメータ等の道路条件と交通需要のもとでの交通状況を再現する。   Next, input / output of the simulation model will be described. First, the simulation engine will be described. The simulation engine reproduces traffic conditions based on road conditions and traffic demands such as road network connection information and geometry, OD traffic volume by time zone specified at arbitrary minute intervals, signal control parameters, and the like.

シミュレーションは一定時間毎に逐次計算を行うピリオデイツクスキャン方式で、任意の集計時間帯毎に、リンクの右左折直進別通過交通量、リンクの右左折直進別平均旅行時間、リンク上の滞留台数、あるいは渋滞長を出力する。指定経路の平均旅行時間や総走行距離、総走行時間などの指標は、これらのリンク毎の指標を用いて算出する。   The simulation is a periodical scan method that performs sequential calculation at regular intervals, and for each total time zone, the amount of traffic that passes through the right and left turns of the link, the average travel time of the right and left turns of the link, and the number of vehicles staying on the link Or the traffic jam length is output. Indices such as the average travel time, total travel distance, and total travel time of the designated route are calculated using the indexes for each link.

SOUNDの入力項目は、以下の(1)〜(8)に示す通りである。
(1)シミュレーション設定
シミュレーションの実行に必要な情報として、シミュレーション対象時間帯、単位スキャン時間、車両発生間隔、結果集計間隔、ネットワークの規模、パケットサイズ、経路選択行動のための情報更新間隔、乱数シード値を入力する。
The SOUND input items are as shown in (1) to (8) below.
(1) Simulation settings Information necessary for execution of simulation includes simulation target time zone, unit scan time, vehicle generation interval, result aggregation interval, network size, packet size, information update interval for route selection behavior, random number seed Enter a value.

(2)ネットワークデータ
主に広域シミュレーションを目的としたSOUNDでは、デジタル道路地図(DRM)ベースのネットワーク構造に対応できる。必要となる情報は、ノード位置、リンク形状・区間長、ノード・リンク接続情報、車線構成(本線車線/右左折付加車線)、通行禁止リンクペアである。
(2) Network data SOUND mainly for the purpose of wide area simulation can cope with a digital road map (DRM) based network structure. Necessary information includes node position, link shape / section length, node / link connection information, lane configuration (main lane / left / right turn additional lane), and a traffic-restricted link pair.

(3)リンク交通特性パラメータ
SOUNDでは、リンクに交通特性を与えて、車両密度の管理を逐次的に行う。交通特性を決定するパラメータは、本線容量[pcu/時/車線]、自由流走行速度[km/時]、ジャム密度[pcu/km]、下流端での右左折直進別飽和交通流率[pcu/有効青1時間/車網である
(3) Link traffic characteristic parameter In SOUND, the traffic density is given to the link, and the vehicle density is sequentially managed. Parameters for determining traffic characteristics are main line capacity [pcu / hour / lane], free-flow traveling speed [km / hour], jam density [pcu / km], right-left turn straight saturated traffic flow rate [pcu] at the downstream end. / Effective blue 1 hour / Vehicle network

(4)セントロイド(ゾーン)情報
シミュレーションでは、交通発生集中点として、ネットワーク端点および街区レベルのゾーン毎にセントロイドを設定するため、道路ネットワークとセントロイドの接続情報として、端点ノートと端点セントロイドの対応、ゾーン内、あるいはゾーン境界のリンクとゾーンセントロイドとの対応を入力する。
(4) Centroid (zone) information In the simulation, since the centroid is set for each zone at the network end point and the block level as the traffic generation concentration point, the end point note and the end point centroid are used as the connection information between the road network and the centroid. The correspondence between the zone centroid and the zone centroid is entered.

(5)OD交通量
SOUNDでは任意の分単位でOD交通量を発生させる。すなわち、 一定時間毎の、車種別、経路選択行動別、その他の属性別OD交通量を入力とする。
(5) OD traffic volume SOUND generates OD traffic volume in arbitrary minutes. In other words, the OD traffic volume by vehicle type, route selection action, and other attributes at regular intervals is input.

(6)経路選択層
SOUNDは複数の経路選択行動を指定できる。前述のOD交通量は、それぞれどの経路選択行動をとるかが設定されている。経路選択行動のパラメータとして、確率的経路選択/最小コスト経路選択、一般化コスト式(経路距離、経路旅行時間、右左折ペナルティの線形和)、ロジットの感度パラメータ(確率的経路選択の場合)を入力する。
(6) Route selection layer SOUND can designate a plurality of route selection actions. Each of the above-mentioned OD traffic volumes sets which route selection action is taken. As parameters for route selection, stochastic route selection / minimum cost route selection, generalized cost formula (linear sum of route distance, route travel time, right / left turn penalty), logit sensitivity parameter (in case of stochastic route selection) input.

(7)信号制御パラメータ
SOUNDでは、交差点に信号を設置することができる。入力情報は、信号制御機番号と設置交差点番号の対応、一連の現示ステップ、各ステップのスプリット時間である。
(7) Signal control parameters In SOUND, signals can be installed at intersections. The input information is the correspondence between the signal controller number and the installation intersection number, a series of display steps, and the split time of each step.

(8)交通規制データ
シミュレーションで評価する対象となる施策のうち、一時的な通行止め規制および車線規制による流入制御を、シミュレーションヘのイベントとして入力することができる。通行止め規制は、対象リンク、対象車種、および時間帯を入力し、車線規制は、対象リンク、閉鎖車線数(容量値によるコントロール)、および時間帯を入力する。
(8) Traffic regulation data Among the measures to be evaluated in the simulation, inflow control based on temporary traffic restriction and lane regulation can be input as an event to the simulation. In the traffic restriction, the target link, the target vehicle type, and the time zone are input, and in the lane control, the target link, the number of closed lanes (control by the capacity value), and the time zone are input.

SOUNDの出力項目は以下の(1)、(2)に示す通りである。
(1)シミュレーション計算実行中に出力するもの
シミュレーション計算中の交通状況を確認するため、1スキャン毎のリンク上の車両位置、各集計時間毎のリンク上の滞留長、各集計時間毎のリンク平均旅行速度を表示部3上のイメージとして表示する。
The output items of SOUND are as shown in (1) and (2) below.
(1) What is output during simulation calculation To confirm the traffic status during simulation calculation, the vehicle position on the link for each scan, the staying length on the link for each counting time, and the link average for each counting time The travel speed is displayed as an image on the display unit 3.

(2)シミュレーション計算終了後に出力するもの
シミュレーション計算終了後は、各集計時間毎に、車種別・進行方向別リンク通過交通量、進行方向別リンク平均旅行時間、リンク上の滞留台数、セントロイド上の滞留台数等の指標を出力する。
(2) What is output after completion of simulation calculation After the completion of simulation calculation, the amount of traffic that passes through the link by vehicle type and direction of travel, the average travel time of link by direction of travel, the number of vehicles staying on the link, Outputs indicators such as the number of staying

平均旅行時間はリンク毎に集計されるため、ある経路や区間の平均旅行時間は、出発時刻からリンク旅行時間をたどることで算出できる。同様に、ある経路や区間、エリアにおける総走行距離・総走行時間などの指標も、リンクの通過交通量とリンク長あるいは平均旅行時間から求めることができる。2次的に求められる指標として、経路/区間平均旅行時間、経路/区間/エリア総走行距離、経路/区間/エリア総走行時間、あるいは総遅れ時間がある。さらに、走行車両の一部をプローブ交通情報として指定することで、プローブ車両が走行したリンクおよびリンク流入時刻の履歴等のサンプル車両走行軌跡を結果として出力することができる。   Since the average travel time is tabulated for each link, the average travel time for a certain route or section can be calculated by following the link travel time from the departure time. Similarly, indicators such as the total travel distance and total travel time in a certain route, section, and area can also be obtained from the link traffic volume and the link length or the average travel time. The second-order index includes route / section average travel time, route / section / area total travel distance, route / section / area total travel time, or total delay time. Furthermore, by designating a part of the traveling vehicle as the probe traffic information, it is possible to output a sample vehicle traveling locus such as a link on which the probe vehicle has traveled and a history of the link inflow time.

図1に示すシミュレーション実行部12において実行するシミュレーションモデルは、SOUNDに限るものではなく、前述の出力指標が得られるシミュレーションモデルであればよい。   The simulation model executed in the simulation execution unit 12 shown in FIG. 1 is not limited to SOUND, and any simulation model can be used as long as the above-described output index can be obtained.

次に、交通流シミュレーションモデル(SOUND)における経路選択について説明する。各車両は常に目的地までの「経路誘導ネットワーク」を参照している。これは、現在走行中のリンクに対して、次に走行するリンクを選ぶ確率が与えられているものである。この確率は一定時間毎に、適切な経路選択モデルによって更新される。経路選択モデルには、必ず最短コストとなる経路を選択するものと、ロジットモデルにより確率選択するものが用意されている。ロジットモデルの選択確率はDialのアルゴリズムを利用して求める。これはあるODペアについてn本の経路があり(例えば、図3に示す例では3本)、各経路のコストが与えられるとき、その経路の選択確率を求めることができる。各車両パケットはリンクに流入した時点で、経路誘導ネットワークを参照し、次に選択するリンクを得ることで、そのリンク下流端での進行方向を決定することができる。   Next, route selection in the traffic flow simulation model (SOUND) will be described. Each vehicle always refers to a “route guidance network” to the destination. This is given the probability of selecting the link to be traveled next for the currently traveling link. This probability is updated by a suitable route selection model at regular intervals. There are two types of route selection models: one that always selects the route with the shortest cost and one that selects the probability using the logit model. The logit model selection probability is determined using the Dial algorithm. In this case, there are n paths for a certain OD pair (for example, three paths in the example shown in FIG. 3), and when the cost of each path is given, the selection probability of the path can be obtained. When each vehicle packet flows into the link, the direction of travel at the downstream end of the link can be determined by referring to the route guidance network and obtaining the next link to be selected.

次に、図2を参照して、図1に示す交通流予測装置1の動作を説明する。図2は、図1に示す交通流予測装置1の動作を示すフローチャートである。まず、データ入力部11は、データファイル4から交通流を予測するのに必要なデータを入力する(ステップS1)。ここで入力するデータは以下の8種類のデータ((1)〜(8)のデータ)である。
(1)プローブ情報:位置情報(緯度、経度)、時刻からなる車両情報
(2)ネットワークデータ:ノードデータ、道路データ(道路線形・道路属性データ)、リンクデータ
(3)リンクパラメータ:リンク容量、ジャム密度、自由流速度、飽和交通流率
(4)交通規制情報:信号情報、車線規制情報、リンク規制情報
(5)ゾーンデータ:ゾーンを定義する情報
(6)経路選択パラメータの初期値(θ):経路選択行動モデルの感度パラメータ
(7)評価値E’のための収束判定値(Th1)
(8)評価値Eのための収束判定値(Th2)、施行回数、前回評価値E(t−1)からの改善率=前回ABS(E(t−1)−E/E
なお、入力データであるプローブ情報に代えて、プローブ情報を集計して得られる時刻、区間番号、区間通過台数、及び区間通過平均所要時間からなるプローブ交通情報を入力するようにしてもよい。
Next, the operation of the traffic flow prediction apparatus 1 shown in FIG. 1 will be described with reference to FIG. FIG. 2 is a flowchart showing the operation of the traffic flow prediction apparatus 1 shown in FIG. First, the data input unit 11 inputs data necessary for predicting traffic flow from the data file 4 (step S1). The data input here are the following eight types of data (data (1) to (8)).
(1) Probe information: vehicle information including position information (latitude, longitude), time (2) network data: node data, road data (road alignment / road attribute data), link data (3) link parameters: link capacity, Jam density, free flow velocity, saturated traffic flow rate (4) Traffic regulation information: signal information, lane regulation information, link regulation information (5) Zone data: Information defining a zone (6) Initial value of path selection parameter (θ ): Sensitivity parameter of route selection behavior model (7) Convergence determination value (Th1) for evaluation value E ′
(8) Convergence judgment value (Th2) for evaluation value E, number of executions, improvement rate from previous evaluation value E (t-1) = previous ABS (E (t-1) -E / E
Instead of the probe information that is input data, probe traffic information that includes the time obtained by aggregating the probe information, the section number, the number of sections that have passed through the section, and the section passage average required time may be input.

次に、シミュレーション実行部12は、データ入力部11において入力した入力データに基づき、乱数系列及びOD交通量を増減させながら交通流シミュレーションの実行を数回繰り返す(ステップS2)ことにより、ゾーン間rsの時間帯hにおける発生交通量Qrsh、ゾーン間時間帯交通量Qrshのうち、時刻tにあるリンクkを通った車両の量Qrsh ktを得る。また、シミュレーション実行部12は、シミュレーション値の集計QKプロットを行い、この集計QKの近似曲線((1)式)を求める。

Figure 2013210941
Next, the simulation execution unit 12 repeats the traffic flow simulation several times while increasing / decreasing the random number sequence and the OD traffic volume based on the input data input in the data input unit 11 (step S2). Among the generated traffic volume Q rsh and the inter-zone time period traffic volume Q rsh in the time zone h, the amount of vehicles Q rsh kt passing through the link k at time t is obtained. In addition, the simulation execution unit 12 performs an aggregate QK plot of simulation values and obtains an approximate curve (equation (1)) of the aggregate QK.
Figure 2013210941

一方、集計QK算出部13は、データ入力部11において入力したプローブ交通情報(観測値)から集計QKを算出する(ステップS3)。そして、集計QK算出部13は、プローブ交通情報の集計QKプロットを行い、この集計QKの近似曲線((2)式)を求める。

Figure 2013210941
On the other hand, the total QK calculation unit 13 calculates the total QK from the probe traffic information (observed values) input in the data input unit 11 (step S3). Then, the total QK calculation unit 13 performs a total QK plot of the probe traffic information and obtains an approximate curve (equation (2)) of the total QK.
Figure 2013210941

次に、拡大率算定部14は、(1)式と(2)式から(3)式、(4式)を満たす拡大率α、βを(5)式、(6)式)によって算定する(ステップS4)。ここで、x’は、拡大後のxの値、y’は、拡大後のyの値である。

Figure 2013210941
Figure 2013210941
Figure 2013210941
Figure 2013210941
Next, the enlargement ratio calculation unit 14 calculates the enlargement ratios α and β satisfying the expressions (3) and (4) from the expressions (1) and (2) by the expressions (5) and (6). (Step S4). Here, x ′ is the value of x after enlargement, and y ′ is the value of y after enlargement.
Figure 2013210941
Figure 2013210941
Figure 2013210941
Figure 2013210941

ここで、図5、図6を参照して拡大率について説明する。図5は、シミュレーションによって求めた集計QKをプロットした結果から求めた(1)式の近似曲線と、プローブ交通情報から求めた集計QKをプロットした結果から求めた(2)式の近似曲線を示す図である。図6は、この2つの近似曲線から前述の拡大率α、βを求め、プローブ交通情報から求めた集計QKの値に拡大率α、βを乗算した結果をプロットした図である。プローブ交通情報は、走行する車両のうちの一部の車両のみであるため、単に集計QKを求めると、図5に示すように、シミュレーションによって求めた近似曲線に比べて、山が小さい曲線となるが、プローブ交通情報に基づき求めた集計QKの値を相似拡大する拡大率を乗算することにより、シミュレーションの集計QKの値より高精度化することが可能となる。   Here, the enlargement ratio will be described with reference to FIGS. FIG. 5 shows the approximate curve of the formula (1) obtained from the result of plotting the total QK obtained by the simulation and the approximate curve of the formula (2) obtained from the result of plotting the total QK obtained from the probe traffic information. FIG. FIG. 6 is a diagram in which the above-described enlargement rates α and β are obtained from these two approximate curves, and the result of multiplying the total QK value obtained from the probe traffic information by the enlargement rates α and β is plotted. Since the probe traffic information is only a part of the traveling vehicles, if the total QK is simply obtained, the curve is smaller than the approximate curve obtained by simulation as shown in FIG. However, by multiplying the value of the total QK obtained based on the probe traffic information by the enlargement factor for similar enlargement, it is possible to make the accuracy higher than the value of the total QK of the simulation.

次に、結果判定部15は、終了条件を満たしたか否かを判定し(ステップS5)、終了条件を満たしていれば処理を終了する。この判定は、(7)式によって得られる評価値E’がデータ入力部11において入力した収束判定値(Th1)以下であるか、または所定の施行回数に達したか否かに基づいて終了条件を満たしたか否かを判定する。

Figure 2013210941
Next, the result determination unit 15 determines whether or not the end condition is satisfied (step S5), and ends the process if the end condition is satisfied. This determination is based on whether the evaluation value E ′ obtained by the expression (7) is equal to or less than the convergence determination value (Th1) input in the data input unit 11 or whether a predetermined number of executions has been reached. It is determined whether or not
Figure 2013210941

次に、結果判定部15により終了条件を満たしていないと判定された場合、パラメータ調整部16は、所定の条件を満たすまでOD交通量とリンク容量を調整する(ステップS6)。ここでいう所定の条件とは、(8)式が評価値Eのための収束判定値(Th2)を下回る、指定された施行回数に達する、改善率がしきい値を下回る、のいずれかの条件である。

Figure 2013210941
Next, when it is determined by the result determination unit 15 that the termination condition is not satisfied, the parameter adjustment unit 16 adjusts the OD traffic volume and the link capacity until a predetermined condition is satisfied (step S6). The predetermined condition here is any one of (8) below the convergence judgment value (Th2) for the evaluation value E, reaching the specified number of enforcement, and the improvement rate below the threshold value. It is a condition.
Figure 2013210941

(8)式において、(7)式と異なる点は、→uzτ sim(→はuの頭に付く、以下、同様)を→uzτと置き換えた点である。ここで、→uzτは、交通流シミュレーションによって求めた値ではなく、パラメータ調整部16が、以下に示す演算によって求めた値である。まず、→uzτは、(9)式によってQτ,Kτを求めることによって得られる。

Figure 2013210941
In the equation (8), the point different from the equation (7) is that → u sim (→ is attached to the head of u, hereinafter the same) is replaced with → u . Here, → u is not a value obtained by the traffic flow simulation but a value obtained by the parameter adjustment unit 16 by the following calculation. First, → u is obtained by obtaining Q τ and K τ by the equation (9).
Figure 2013210941

次に、(9)式のnjτは、(10)式によって求められ、(9)式のTjτは、(11)式によって求めることができる。(11)式においてqjτは、(12)式によって求めることができる。

Figure 2013210941
Figure 2013210941
Figure 2013210941
Next, n jτ in equation (9) can be obtained from equation (10), and T jτ in equation (9) can be obtained from equation (11). In equation (11), q can be obtained by equation (12).
Figure 2013210941
Figure 2013210941
Figure 2013210941

パラメータ調整部16は、(8)式の評価値Eが前述の条件を満たすまで、(10)式、(12)式で用いるαの値(リンク容量)と(12)式のQrsh(OD交通量)を増減させて調整を行う。このとき、αとQrshを増減させて最適化を図る手法は、公知の最適化手法を用いることが可能であるので、ここでは詳細な説明を省略する。パラメータ調整部16は、この処理によって求められた最適化されたαとQrshをシミュレーション実行部12へ出力する。 The parameter adjustment unit 16 determines the value of α (link capacity) used in Equations (10) and (12) and Q rsh (OD) in Equation (12) until the evaluation value E in Equation (8) satisfies the above-described condition. Adjust by increasing or decreasing (traffic volume). At this time, as a method of optimizing by increasing or decreasing α and Q rsh , a known optimization method can be used, and thus detailed description thereof is omitted here. The parameter adjustment unit 16 outputs the optimized α and Q rsh obtained by this processing to the simulation execution unit 12.

これを受けて、シミュレーション実行部12は、パラメータ調整部16から出力されたα(リンク容量)の値とQrsh(OD交通量)の値を入力して交通流シミュレーションを実行する(ステップS7)。そして、再度結果判定部15は、(7)式によって得られる評価値E’が収束判定値(Th1)以下であるか、または所定の施行回数に達したか否かに基づいて終了条件を満たしたか否かを判定する(ステップS5)。そして、結果判定部15は、終了条件を満たしていれば処理を終了し、満たしていなければステップS6、S7の処理を繰り返し行い、ステップS5における終了条件を満たすまで処理を繰り返す。 In response to this, the simulation execution unit 12 inputs the value of α (link capacity) and the value of Q rsh (OD traffic volume) output from the parameter adjustment unit 16 and executes a traffic flow simulation (step S7). . Then, the result determination unit 15 again satisfies the termination condition based on whether the evaluation value E ′ obtained by the expression (7) is equal to or less than the convergence determination value (Th1) or whether the predetermined number of executions has been reached. It is determined whether or not (step S5). Then, the result determination unit 15 ends the process if the end condition is satisfied, repeats the processes in steps S6 and S7 if not satisfied, and repeats the process until the end condition in step S5 is satisfied.

従来のOD交通量推定では、日通過交通量から、1日単位のOD交通量を推定するものが一般的である。すなわち、OD交通量の時間変動を考慮できないという問題がある。また、時間帯別OD交通量を推定するものも存在する(例えば、特許文献1)が、評価値を得るために交通流シミュレーションを毎回実行しているため、非常に計算時間が長くなるという課題を抱えている。   In the conventional OD traffic volume estimation, the daily OD traffic volume is generally estimated from the daily passage traffic volume. That is, there is a problem that time fluctuation of the OD traffic volume cannot be considered. Moreover, although there exists what estimates the OD traffic volume according to time zone (for example, patent document 1), since traffic flow simulation is performed every time in order to obtain an evaluation value, the subject that calculation time becomes very long Have

以上説明したように、本実施形態では、プローブ交通情報から得られる交通状態を用いて時間帯別のOD交通量及びリンク容量パラメータを求めるようにしたため、得られる値の精度を向上させることができる。また、既存のモデルでは目標値が交通量であるのに対し、本実施形態は、プローブ交通情報から集計されたエリア内存在台数(Q)と総通過交通量(K)から得られる集計QK(交通密度―交通量)を目標値としているため、実際の交通状態に則した交通流予測を行うことができる。また、交通流シミュレーションを利用するタイミングを初期計算時とパラメータ収束計算がある目的を達成した際の再評価値計算時のみに留めたため、必要最小限の利用回数で交通流予測が可能になる。   As described above, in this embodiment, since the OD traffic volume and link capacity parameter for each time zone are obtained using the traffic state obtained from the probe traffic information, the accuracy of the obtained value can be improved. . In addition, while the target value is the traffic volume in the existing model, in the present embodiment, the total QK () obtained from the total number of traffic in the area (Q) and the total passing traffic volume (K) calculated from the probe traffic information. Since traffic density (traffic volume) is a target value, traffic flow prediction according to actual traffic conditions can be performed. In addition, since the timing of using the traffic flow simulation is limited to the initial calculation and the re-evaluation value calculation when the purpose of the parameter convergence calculation is achieved, the traffic flow prediction can be performed with the minimum necessary number of uses.

なお、図1における処理部の機能を実現するためのプログラムをコンピュータ読み取り可能な記録媒体に記録して、この記録媒体に記録されたプログラムをコンピュータシステムに読み込ませ、実行することにより交通流予測処理を行ってもよい。なお、ここでいう「コンピュータシステム」とは、OSや周辺機器等のハードウェアを含むものとする。また、「コンピュータシステム」は、ホームページ提供環境(あるいは表示環境)を備えたWWWシステムも含むものとする。また、「コンピュータ読み取り可能な記録媒体」とは、フレキシブルディスク、光磁気ディスク、ROM、CD−ROM等の可搬媒体、コンピュータシステムに内蔵されるハードディスク等の記憶装置のことをいう。さらに「コンピュータ読み取り可能な記録媒体」とは、インターネット等のネットワークや電話回線等の通信回線を介してプログラムが送信された場合のサーバやクライアントとなるコンピュータシステム内部の揮発性メモリ(RAM)のように、一定時間プログラムを保持しているものも含むものとする。   Note that a traffic flow prediction process is performed by recording a program for realizing the function of the processing unit in FIG. 1 on a computer-readable recording medium, causing the computer system to read and execute the program recorded on the recording medium. May be performed. Here, the “computer system” includes an OS and hardware such as peripheral devices. The “computer system” includes a WWW system having a homepage providing environment (or display environment). The “computer-readable recording medium” refers to a storage device such as a flexible medium, a magneto-optical disk, a portable medium such as a ROM and a CD-ROM, and a hard disk incorporated in a computer system. Further, the “computer-readable recording medium” refers to a volatile memory (RAM) in a computer system that becomes a server or a client when a program is transmitted via a network such as the Internet or a communication line such as a telephone line. In addition, those holding programs for a certain period of time are also included.

また、上記プログラムは、このプログラムを記憶装置等に格納したコンピュータシステムから、伝送媒体を介して、あるいは、伝送媒体中の伝送波により他のコンピュータシステムに伝送されてもよい。ここで、プログラムを伝送する「伝送媒体」は、インターネット等のネットワーク(通信網)や電話回線等の通信回線(通信線)のように情報を伝送する機能を有する媒体のことをいう。また、上記プログラムは、前述した機能の一部を実現するためのものであってもよい。さらに、前述した機能をコンピュータシステムにすでに記録されているプログラムとの組み合わせで実現できるもの、いわゆる差分ファイル(差分プログラム)であってもよい。   The program may be transmitted from a computer system storing the program in a storage device or the like to another computer system via a transmission medium or by a transmission wave in the transmission medium. Here, the “transmission medium” for transmitting the program refers to a medium having a function of transmitting information, such as a network (communication network) such as the Internet or a communication line (communication line) such as a telephone line. The program may be for realizing a part of the functions described above. Furthermore, what can implement | achieve the function mentioned above in combination with the program already recorded on the computer system, what is called a difference file (difference program) may be sufficient.

車両の交通流を高精度で予測することが不可欠な用途に適用できる。   It can be applied to applications where it is essential to predict the traffic flow of vehicles with high accuracy.

1・・・交通流予測装置、2・・・入力部、3・・・表示部、4・・・データファイル、11・・・データ入力部、12・・・シミュレーション実行部、13・・・集計QK算出部、14・・・拡大率算定部、15・・・結果判定部、16・・・パラメータ調整部   DESCRIPTION OF SYMBOLS 1 ... Traffic flow prediction apparatus, 2 ... Input part, 3 ... Display part, 4 ... Data file, 11 ... Data input part, 12 ... Simulation execution part, 13 ... Total QK calculation unit, 14 ... enlargement rate calculation unit, 15 ... result determination unit, 16 ... parameter adjustment unit

Claims (3)

観測値である時刻、緯度経度情報からなるプローブ情報、もしくは時刻、区間番号、区間通過台数、及び区間通過平均所要時間からなるプローブ交通情報から、任意に指定した区域の交通密度と交通量の関係を定量化し、該関係を用いて交通量を推定する交通流予測装置であって、
前記区域を走行した車両のプローブ交通情報に基づき、前記区域内の交通密度と交通量の関係を示す第1の近似曲線を求める第1の近似曲線算出手段と、
交通流シミュレーションモデルを用いて交通流シミュレーションを行った結果に基づき、前記区域内の交通密度と交通量の関係を示す第2の近似曲線を求める第2の近似曲線算出手段と、
前記第1の近似曲線を前記第2の近似曲線に合致させるための拡大率を算定する拡大率算定手段と、
前記第1の近似曲線を求めた際の前記区域内の交通密度と交通量の関係を示す値のそれぞれに対して、前記拡大率を乗算して得られた前記区域内の交通密度と交通量の関係を示す値が、前記第2の近似曲線を求めた際の前記区域内の交通密度と交通量の関係を示す値との距離が最小化させるための交通量とリンク容量を求めるパラメータ調整手段と、
前記パラメータ調整手段によって求めた前記交通量と前記リンク容量を入力として前記交通流シミュレーションを実行した結果が所定の終了条件を満たすまでパラメータ調整手段による前記交通量と前記リンク容量を求める処理と前記交通流シミュレーションの実行を繰り返す結果判定手段と
を備えたことを特徴とする交通流予測装置。
Relationship between traffic density and traffic volume in an arbitrarily specified area from probe information consisting of observed time, latitude and longitude information, or probe traffic information consisting of time, section number, number of sections passing, and section average travel time Is a traffic flow prediction device that estimates traffic volume using the relationship,
First approximate curve calculation means for obtaining a first approximate curve indicating a relationship between traffic density and traffic volume in the area based on probe traffic information of a vehicle traveling in the area;
A second approximate curve calculating means for obtaining a second approximate curve indicating a relationship between the traffic density and the traffic volume in the area based on the result of the traffic flow simulation using the traffic flow simulation model;
An enlargement ratio calculating means for calculating an enlargement ratio for matching the first approximate curve with the second approximate curve;
The traffic density and the traffic volume in the area obtained by multiplying each of the values indicating the relationship between the traffic density and the traffic volume in the area when the first approximate curve is obtained by multiplying the expansion rate. Is a parameter adjustment for determining the traffic volume and link capacity for minimizing the distance between the traffic density in the area and the value indicating the traffic volume when the second approximate curve is obtained. Means,
Processing for obtaining the traffic volume and the link capacity by the parameter adjustment means until the result of executing the traffic flow simulation with the traffic volume and the link capacity obtained by the parameter adjustment means as input satisfies a predetermined termination condition, and the traffic A traffic flow prediction apparatus comprising: a result determination unit that repeatedly executes a flow simulation.
観測値である時刻、緯度経度情報からなるプローブ情報、もしくは時刻、区間番号、区間通過台数、及び区間通過平均所要時間からなるプローブ交通情報から、任意に指定した区域の交通密度と交通量の関係を定量化し、該関係を用いて交通量を推定する交通流予測装置が行う交通流予測方法であって、
前記区域を走行した車両のプローブ交通情報に基づき、前記区域内の交通密度と交通量の関係を示す第1の近似曲線を求める第1の近似曲線算出ステップと、
交通流シミュレーションモデルを用いて交通流シミュレーションを行った結果に基づき、前記区域内の交通密度と交通量の関係を示す第2の近似曲線を求める第2の近似曲線算出ステップと、
前記第1の近似曲線を前記第2の近似曲線に合致させるための拡大率を算定する拡大率算定ステップと、
前記第1の近似曲線を求めた際の前記区域内の交通密度と交通量の関係を示す値のそれぞれに対して、前記拡大率を乗算して得られた前記区域内の交通密度と交通量の関係を示す値が、前記第2の近似曲線を求めた際の前記区域内の交通密度と交通量の関係を示す値との距離が最小化させるための交通量とリンク容量を求めるパラメータ調整ステップと、
前記パラメータ調整ステップによって求めた前記交通量と前記リンク容量を入力として前記交通流シミュレーションを実行した結果が所定の終了条件を満たすまでパラメータ調整ステップによる前記交通量と前記リンク容量を求める処理と前記交通流シミュレーションの実行を繰り返す結果判定ステップと
を備えたことを特徴とする交通流予測方法。
Relationship between traffic density and traffic volume in an arbitrarily specified area from probe information consisting of observed time, latitude and longitude information, or probe traffic information consisting of time, section number, number of sections passing, and section average travel time Is a traffic flow prediction method performed by a traffic flow prediction device that estimates traffic using the relationship,
A first approximate curve calculating step for obtaining a first approximate curve indicating a relationship between traffic density and traffic volume in the area based on probe traffic information of a vehicle traveling in the area;
A second approximate curve calculating step for obtaining a second approximate curve indicating a relationship between the traffic density and the traffic volume in the area based on the result of the traffic flow simulation using the traffic flow simulation model;
An enlargement factor calculating step for calculating an enlargement factor for matching the first approximate curve with the second approximate curve;
The traffic density and the traffic volume in the area obtained by multiplying each of the values indicating the relationship between the traffic density and the traffic volume in the area when the first approximate curve is obtained by multiplying the expansion rate. Is a parameter adjustment for determining the traffic volume and link capacity for minimizing the distance between the traffic density in the area and the value indicating the traffic volume when the second approximate curve is obtained. Steps,
Processing for obtaining the traffic volume and the link capacity in the parameter adjustment step until the result of executing the traffic flow simulation with the traffic volume and the link capacity obtained in the parameter adjustment step as inputs is satisfied a predetermined end condition, and the traffic A traffic flow prediction method comprising: a result determination step that repeats execution of a flow simulation.
観測値である時刻、緯度経度情報からなるプローブ情報、もしくは時刻、区間番号、区間通過台数、及び区間通過平均所要時間からなるプローブ交通情報から、任意に指定した区域の交通密度と交通量の関係を定量化し、該関係を用いて交通量を推定する交通流予測装置上のコンピュータに、
前記区域を走行した車両のプローブ交通情報に基づき、前記区域内の交通密度と交通量の関係を示す第1の近似曲線を求める第1の近似曲線算出ステップと、
交通流シミュレーションモデルを用いて交通流シミュレーションを行った結果に基づき、前記区域内の交通密度と交通量の関係を示す第2の近似曲線を求める第2の近似曲線算出ステップと、
前記第1の近似曲線を前記第2の近似曲線に合致させるための拡大率を算定する拡大率算定ステップと、
前記第1の近似曲線を求めた際の前記区域内の交通密度と交通量の関係を示す値のそれぞれに対して、前記拡大率を乗算して得られた前記区域内の交通密度と交通量の関係を示す値が、前記第2の近似曲線を求めた際の前記区域内の交通密度と交通量の関係を示す値との距離が最小化させるための交通量とリンク容量を求めるパラメータ調整ステップと、
前記パラメータ調整ステップによって求めた前記交通量と前記リンク容量を入力として前記交通流シミュレーションを実行した結果が所定の終了条件を満たすまでパラメータ調整ステップによる前記交通量と前記リンク容量を求める処理と前記交通流シミュレーションの実行を繰り返す結果判定ステップと
を行わせることを特徴とする交通流予測プログラム。
Relationship between traffic density and traffic volume in an arbitrarily specified area from probe information consisting of observed time, latitude and longitude information, or probe traffic information consisting of time, section number, number of sections passing, and section average travel time To the computer on the traffic flow prediction device that estimates the traffic volume using the relationship,
A first approximate curve calculating step for obtaining a first approximate curve indicating a relationship between traffic density and traffic volume in the area based on probe traffic information of a vehicle traveling in the area;
A second approximate curve calculating step for obtaining a second approximate curve indicating a relationship between the traffic density and the traffic volume in the area based on the result of the traffic flow simulation using the traffic flow simulation model;
An enlargement factor calculating step for calculating an enlargement factor for matching the first approximate curve with the second approximate curve;
The traffic density and the traffic volume in the area obtained by multiplying each of the values indicating the relationship between the traffic density and the traffic volume in the area when the first approximate curve is obtained by multiplying the expansion rate. Is a parameter adjustment for determining the traffic volume and link capacity for minimizing the distance between the traffic density in the area and the value indicating the traffic volume when the second approximate curve is obtained. Steps,
Processing for obtaining the traffic volume and the link capacity in the parameter adjustment step until the result of executing the traffic flow simulation with the traffic volume and the link capacity obtained in the parameter adjustment step as inputs is satisfied a predetermined end condition, and the traffic A traffic flow prediction program characterized by causing a result determination step to repeatedly execute a flow simulation.
JP2012081996A 2012-03-30 2012-03-30 Traffic flow prediction device, traffic flow prediction method, and traffic flow prediction program Expired - Fee Related JP5898553B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2012081996A JP5898553B2 (en) 2012-03-30 2012-03-30 Traffic flow prediction device, traffic flow prediction method, and traffic flow prediction program

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2012081996A JP5898553B2 (en) 2012-03-30 2012-03-30 Traffic flow prediction device, traffic flow prediction method, and traffic flow prediction program

Publications (2)

Publication Number Publication Date
JP2013210941A true JP2013210941A (en) 2013-10-10
JP5898553B2 JP5898553B2 (en) 2016-04-06

Family

ID=49528684

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2012081996A Expired - Fee Related JP5898553B2 (en) 2012-03-30 2012-03-30 Traffic flow prediction device, traffic flow prediction method, and traffic flow prediction program

Country Status (1)

Country Link
JP (1) JP5898553B2 (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105513357A (en) * 2015-12-30 2016-04-20 银江股份有限公司 Method for calibrating traffic flow basic map parameter on the basis of microwave data
KR101742043B1 (en) 2016-11-15 2017-05-31 한국과학기술정보연구원 Apparatus and method for travel mode choice prediction
KR20200007577A (en) * 2018-07-13 2020-01-22 서울시립대학교 산학협력단 Traffic estimation method using average travel time information and traffic estimation device
CN114516329A (en) * 2020-11-17 2022-05-20 丰田自动车株式会社 Vehicle adaptive cruise control system, method, and computer readable medium
CN115457764A (en) * 2022-08-24 2022-12-09 华南理工大学 Road section traffic density estimation method, device and medium based on vehicle track data

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112017436B (en) * 2020-09-09 2021-09-28 中国科学院自动化研究所 Method and system for predicting urban traffic travel time

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000222669A (en) * 1999-01-28 2000-08-11 Mitsubishi Electric Corp Traffic flow estimating device and traffic flow estimating method
US20020198694A1 (en) * 2001-06-22 2002-12-26 Qi Yang Traffic data management and simulation system
JP2005084722A (en) * 2003-09-04 2005-03-31 Toshiba Corp Road traffic condition analysis apparatus and prediction apparatus
JP2011048754A (en) * 2009-08-28 2011-03-10 I-Transport Lab Co Ltd Device, program and method for analyzing traffic situation

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000222669A (en) * 1999-01-28 2000-08-11 Mitsubishi Electric Corp Traffic flow estimating device and traffic flow estimating method
US20020198694A1 (en) * 2001-06-22 2002-12-26 Qi Yang Traffic data management and simulation system
JP2005084722A (en) * 2003-09-04 2005-03-31 Toshiba Corp Road traffic condition analysis apparatus and prediction apparatus
JP2011048754A (en) * 2009-08-28 2011-03-10 I-Transport Lab Co Ltd Device, program and method for analyzing traffic situation

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105513357A (en) * 2015-12-30 2016-04-20 银江股份有限公司 Method for calibrating traffic flow basic map parameter on the basis of microwave data
KR101742043B1 (en) 2016-11-15 2017-05-31 한국과학기술정보연구원 Apparatus and method for travel mode choice prediction
KR20200007577A (en) * 2018-07-13 2020-01-22 서울시립대학교 산학협력단 Traffic estimation method using average travel time information and traffic estimation device
KR102116908B1 (en) 2018-07-13 2020-06-01 서울시립대학교 산학협력단 Traffic estimation method using average travel time information and traffic estimation device
CN114516329A (en) * 2020-11-17 2022-05-20 丰田自动车株式会社 Vehicle adaptive cruise control system, method, and computer readable medium
JP2022080304A (en) * 2020-11-17 2022-05-27 トヨタ自動車株式会社 Vehicle adaptive cruise control system, method, computer program for executing the method, and computer readable medium
CN115457764A (en) * 2022-08-24 2022-12-09 华南理工大学 Road section traffic density estimation method, device and medium based on vehicle track data
CN115457764B (en) * 2022-08-24 2023-07-18 华南理工大学 Road section traffic density estimation method, device and medium based on vehicle track data

Also Published As

Publication number Publication date
JP5898553B2 (en) 2016-04-06

Similar Documents

Publication Publication Date Title
Liu et al. A cellular automata traffic flow model combined with a BP neural network based microscopic lane changing decision model
EP3035314B1 (en) A traffic data fusion system and the related method for providing a traffic state for a network of roads
Ton et al. Evaluating a data-driven approach for choice set identification using GPS bicycle route choice data from Amsterdam
JP5898553B2 (en) Traffic flow prediction device, traffic flow prediction method, and traffic flow prediction program
JP2013235326A (en) Traffic flow prediction device, traffic flow prediction method, and traffic flow prediction program
Minelli et al. Evaluation of connected vehicle impact on mobility and mode choice
Tang et al. Markov Chains based route travel time estimation considering link spatio-temporal correlation
Moshahedi et al. A network-wide anticipatory control of an urban network using macroscopic fundamental diagram
KR100772511B1 (en) Method for searching the shortest route based on traffic prediction and apparatus thereof
Kawasaki et al. Real-time monitoring of dynamic traffic states by state-space model
Lee et al. Development of reinforcement learning-based traffic predictive route guidance algorithm under uncertain traffic environment
Rodriguez-Vega et al. Dynamic density and flow reconstruction in large-scale urban networks using heterogeneous data sources
Pan et al. Modified volume-delay function based on traffic fundamental diagram: A practical calibration framework for estimating congested and uncongested conditions
Yang et al. Connected vehicle enhanced vehicle routing with intersection turning cost estimation
Ng et al. Development of Macroscopic Cell‐Based Logistic Lane Change Prediction Model
Yu et al. Genetic-fuzzy logic ramp metering control for ramp metering of motorway
Guastella et al. Traffic Modeling with SUMO: a Tutorial
Darda Joint calibration of a microscope traffic simulator and estimation of origin-destination flows
Abedini A Machine-Learning Framework for Clustering and Calibration of Roadway Performance Models with Application in the Large-Scale Traffic Assignment
Allström Highway Traffic State Estimation and Short-term Prediction
Zhu A connected vehicle based coordinated adaptive navigation system
Roocroft et al. Data-driven traffic assignment through density-based road-specific congestion function estimation
Zahedian Introducing a Graph-Based Neural Network for Networkwide Traffic Volume Estimation
Pandey Optimal dynamic pricing for managed lanes with multiple entrances and exits
Kabir Data Driven Method to Assess Safety and Energy on Freeway and Intersection Performance Using Probe Vehicle

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20141211

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20150819

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20150825

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20151019

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20160304

R150 Certificate of patent or registration of utility model

Ref document number: 5898553

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

LAPS Cancellation because of no payment of annual fees