(57)【要約】  センサの未設置区間における渋滞起点、及び、渋滞末尾の特定精度を向上させる。渋滞推定装置は、車両ごとに検出された位置情報を基に、渋滞領域を定めるための渋滞末尾および渋滞起点を特定する渋滞起点・末尾特定部と、前記渋滞末尾から上流側のセンサで測定された交通量と、前記渋滞起点から下流側のセンサで測定された交通量と、前記位置情報を基に、時刻ごとの前記渋滞末尾への流入量および前記時刻ごとの前記渋滞起点からの流出量を推定する流入・流出量推定部と、前記流入量を基に前記渋滞末尾を更新し、前記流出量を基に前記渋滞末尾起点を更新する渋滞領域更新部を備える。
2018180688.tif 000002
【特許請求の範囲】
【請求項１】
車両ごとに検出された位置情報を基に、渋滞領域を定めるための渋滞末尾、及び、渋滞起点を特定する、渋滞起点・末尾特定手段と、
前記渋滞末尾から上流側のセンサで測定された交通量と、前記渋滞起点から下流側のセンサで測定された交通量と、前記位置情報とを基に、時刻ごとの前記渋滞末尾への流入量および前記時刻ごとの前記渋滞起点からの流出量を推定する、流入・流出量推定手段と、
前記流入量を基に前記渋滞末尾を更新し、前記流出量を基に前記渋滞起点を更新する、渋滞領域更新手段と、
を備える、渋滞推定装置。
【請求項２】
前記流入・流出量推定手段は、
前記上流側のセンサで測定された交通量と、前記上流側のセンサを通過した車両が前記渋滞末尾に到着する割合を示す到着率とを基に、到達流入量を算出し、
前記到達流入量と、ある時刻より前に前記上流側のセンサで測定された交通量のうち、ある時刻に前記渋滞末尾に到着する残留交通量とを基に、前記流入量を算出する、
請求項１に記載の渋滞推定装置。
【請求項３】
前記流入・流出量推定手段は、
ある時刻より前に前記上流側のセンサと前記渋滞末尾との間に残留した車両台数と、ある時刻に前記渋滞末尾に到達可能な車両の割合とを基に、前記残留交通量を算出する、
請求項２に記載の渋滞推定装置。
【請求項４】
前記流入・流出量推定手段は、
前記下流側のセンサで測定された交通量と前記渋滞起点から離脱した車両の割合を示す離脱率とに基づいて離脱交通量を算出し、
前記離脱交通量と、ある時刻に前記渋滞起点を離脱し、前記ある時刻より後に前記下流側のセンサに到着する離脱残留交通量とを基に、前記流出量を算出する、
請求項１乃至３のいずれか１つに記載の渋滞推定装置。
【請求項５】
前記流入・流出量推定手段は、
ある時刻より後に前記渋滞起点と前記下流側のセンサとの間に残留した残留台数と、ある時刻より後に前記下流側のセンサに到着する割合とを基に、前記離脱残留交通量を算出する、
請求項４に記載の渋滞推定装置。
【請求項６】
前記渋滞起点・末尾特定手段は、
前記位置情報を基に、前記車両の平均速度又は交通密度を算出して前記渋滞末尾および前記渋滞起点を特定する、
請求項１乃至５のいずれか１つに記載の渋滞推定装置。
【請求項７】
前記渋滞領域更新手段は、
前記流入量と前記流出量との差を基に、前記渋滞領域の延伸を更新する、
請求項１乃至６のいずれか１つに記載の渋滞推定装置。
【請求項８】
車両ごとに検出された位置情報を基に、渋滞領域を定めるための渋滞末尾、及び、渋滞起点を特定し、
前記渋滞末尾から上流側のセンサで測定された交通量と、前記渋滞起点から下流側のセンサで測定された交通量と、前記位置情報とを基に、時刻ごとの前記渋滞末尾への流入量および前記時刻ごとの前記渋滞起点からの流出量を推定し、
前記流入量を基に前記渋滞末尾を更新し、前記流出量を基に前記渋滞起点を更新する、
渋滞推定方法。
【請求項９】
車両ごとに検出された位置情報を基に、渋滞領域を定めるための渋滞末尾、及び、渋滞起点を特定し、
前記渋滞末尾から上流側のセンサで測定された交通量と、前記渋滞起点から下流側のセンサで測定された交通量と、前記位置情報とを基に、時刻ごとの前記渋滞末尾への流入量および前記時刻ごとの前記渋滞起点からの流出量を推定し、
前記流入量を基に前記渋滞末尾を更新し、前記流出量を基に前記渋滞起点を更新する、
ことをコンピュータに実行させる、プログラムを格納した記録媒体。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、渋滞推定装置等に関する。
【背景技術】
【０００２】
道路管理者にとって、道路の各地点における交通状況を計測して監視することは道路の安全管理、道路の品質維持において重要である。特に、高速道路において、渋滞が発生する起点及び渋滞の末尾を特定することは交通量の制御又は道路の増設、拡幅計画の策定にとっても重要である。そのため、さまざまな計測機器（センサ）が道路に設置され、それらの情報を収集し活用する交通管制システムが構築されている。【０００３】
計測機器には、ループコイルによる車両検知器、ＣＣＴＶ（ｃｌｏｓｅｄ−ｃｉｒｃｕｉｔ  ｔｅｌｅｖｉｓｉｏｎ）カメラによって撮像映像を画像処理し、車体を認識する車両検知器、超音波を用いた車両検知器、あるいは、光学式車両検知器等が用いられる。【０００４】
しかしながら、これらの設置型センサによる交通状況の監視には、設置型センサだけでなく、その設置型センサの情報を処理する処理装置、さらに処理装置と交通管制システムとの通信回線が必要であり設置コストがかかる。そのため、全ての道路区間に対して設置型センサを配備することは困難であり、渋滞が発生しやすい道路区間など重要な道路区間にのみ配備される場合が多い。一方で、重要な道路区間以外でも、事故又は非定常な渋滞の発生時には安全管理上、状況把握が必要となるが、現状では監視が難しい。【０００５】
また、非設置型センサを用いない交通状況測定として、ＧＰＳ（Ｇｌｏｂａｌ  Ｐｏｓｉｔｉｏｎｉｎｇ  Ｓｙｓｔｅｍ）車載器を用いた交通状況測定が検討されている。ＧＰＳによる交通状況測定の利点は、道路管理者の費用負担が小さいこと、設置型センサの無い地点でも交通状況が測定できることである。【０００６】
一方で、ＧＰＳによる交通状況測定にも以下のようなデメリットがある。例えば、ＧＰＳ車載器を搭載した車両が道路を走行しないと交通状況が測定できない。このため、交通状況の測定点が時間と地域によって大きな偏りが発生する。また、ＧＰＳの誤差の影響を受けるため、特に衛星が把握しづらい都心部では測定誤差が増大しやすい。さらにＧＰＳ搭載車両の走行速度および移動時間が測定できたとしても、その道路区間全体としての情報、例えば平均走行速度、通過台数又は密度などの情報を知ることができない。【０００７】
このような状況を鑑み、上記の設置型センサによる交通状況測定と、ＧＰＳによる交通状況測定を組み合わせることで時間経過に応じて変化する渋滞領域の推定技術が検討されている。【０００８】
特許文献１は、主に一般道の交差点間で流入・流出車両が測定できない分岐路が存在する場合に、交差点間に滞留している車両台数を推定する方法を開示している。特許文献１の推定方法では、センサが設置された二つの交差点間において、流入側交差点の通過台数と流出側交差点の通過台数の差と、未知の分岐路への流出及び分岐路からの流入の割合を示す補正係数を用いて、測定できていない道路区間に存在する車両台数を推定している。【０００９】
さらに、特許文献１は、補正係数を更新する際にＧＰＳ情報を用いて補正係数の精度を向上させる方法も開示する。前記の通り、補正係数は未測定の分岐路との流出入割合を表しており、特許文献１に記載の方法ではこの補正係数は特定の周期で過去の測定データを用いて更新される。この更新は流入側交差点での通過台数と流出側交差点での交差点の比に基づいて更新されるが、流出側の交差点において測定される交通量は、流入側においては同時刻の計測データではなく、より過去の時刻で計測された交通量と相関がある。車両の移動時間を加味した上で流入側・流出側での測定時間帯を合わせ、その交通流において、未測定の分岐路との流出入を推定することが望ましい。そのため、ＧＰＳ情報を用いて車両の移動時間を測定し、その時間に基づいて流出側通過台数と流入側通過台数の時刻を適正に設定している。【００１０】
また、特許文献２も同じく設置型センサで測定された情報と、ＧＰＳによって測定された情報を組み合わせて未測定の地点の交通状況を推定する方法を開示している。特許文献２では速度を測定・推定対象としている。具体的には、ＧＰＳによる測定速度を教師データととらえ、設置型センサの測定速度を説明変数、ＧＰＳによる測定速度を目的変数とした回帰分析を行うことで設置型センサの情報から精度の高い速度情報を任意の時間で推定する方法を開示している。また、この方法は説明変数に他の地点の設置型センサ情報、目的変数にセンサ未設置区間のＧＰＳ測定情報を用いることもできる。【００１１】
特許文献３も、設置型センサで測定された渋滞末尾情報に加えて、ＧＰＳによって測定された情報を組み合わせて、渋滞の末尾の測定精度を向上させる方法を開示する。特許文献３の方法は、ＧＰＳによって測定された車両の位置情報を収集し、その位置間隔が一定以下となる車両が渋滞領域に存在していると推定する。そして特許文献３の方法は、推定した渋滞領域のもっとも上流の地点を渋滞末尾ととらえることによって、設置型センサの無い区間においても渋滞末尾の変化を測定する。
【先行技術文献】
【特許文献】
【００１２】
【特許文献１】国際公開第２０１５／０４５６９５号
【特許文献２】米国特許出願公開第２０１０／０２８６８９９号明細書
【特許文献３】特開２００１−１８８９８７号公報
【発明の概要】
【発明が解決しようとする課題】
【００１３】
特許文献１が開示する技術によって、ＧＰＳで測定された道路移動時間を考慮した上で、未測定道路区間における車両の残留割合が推定できる。本方式を高速道路に適用する際には、オンランプ・オフランプを分岐路として捉えることで同様に適用することが出来る。しかしながら、特許文献１の方法は、センサ設置点で測定された単位時間当たりの交通量が、分岐路における流出入を除いてセンサ未設置区間においても均一に割り当てられるという前提がある。そのため、山間部の高速道路など、センサ未設置区間が長距離にわたる場合、車両ごとの平均速度の違いによって地点ごとに交通量が均一ではなくなり、推定精度が低くなるため、渋滞の推定精度も低下する。【００１４】
特許文献２の技術では、ＧＰＳによる測定と合うように未測定地点の推定を行う回帰モデルを構築する。しかし、このモデル構築はＧＰＳ情報を教師データとするため、十分なＧＰＳ情報がない場合には推定できない地点数が増加し、渋滞の特定が困難となる。【００１５】
特許文献３の技術も同じくＧＰＳによる測定のみに基づいてセンサの未設置区間における渋滞変化を推定しているため、ＧＰＳデータが時間的・空間的に細粒度に収集させる環境でなければ渋滞の推定精度が低下する。【００１６】
本発明の目的は、センサの未設置区間であっても渋滞領域の起点、及び、末尾を特定する精度を向上させることが可能な渋滞推定装置等を提供することにある。
【課題を解決するための手段】
【００１７】
本発明の渋滞推定装置の一態様は、車両ごとに検出された位置情報を基に、渋滞領域を定めるための渋滞末尾および渋滞起点を特定する渋滞起点・末尾特定部と、前記渋滞末尾から上流側のセンサで測定された交通量と、前記渋滞起点から下流側のセンサで測定された交通量と、前記位置情報を基に、時刻ごとの前記渋滞末尾への流入量および前記時刻ごとの前記渋滞起点からの流出量を推定する流入・流出量推定部と、前記流入量を基に前記渋滞末尾を更新し、前記流出量を基に前記渋滞末尾起点を更新する渋滞領域更新部とを備える。【００１８】
本発明の渋滞推定方法の一態様は、車両ごとに検出された位置情報を基に、渋滞領域を定めるための渋滞末尾および渋滞起点を特定し、前記渋滞末尾から上流側のセンサで測定された交通量と、前記渋滞起点から下流側のセンサで測定された交通量と、前記位置情報を基に、時刻ごとの前記渋滞末尾への流入量および前記時刻ごとの前記渋滞起点からの流出量を推定し、前記流入量を基に前記渋滞末尾を更新し、前記流出量を基に前記渋滞末尾起点を更新する。【００１９】
本発明の記録媒体に格納されたプログラムの一態様は、車両ごとに検出された位置情報を基に、渋滞領域を定めるための渋滞末尾および渋滞起点を特定し、前記渋滞末尾から上流側のセンサで測定された交通量と、前記渋滞起点から下流側のセンサで測定された交通量と、前記位置情報を基に、時刻ごとの前記渋滞末尾への流入量および前記時刻ごとの前記渋滞起点からの流出量を推定し、前記流入量を基に前記渋滞末尾を更新し、前記流出量を基に前記渋滞末尾起点を更新することをコンピュータに実行させる。
【発明の効果】
【００２０】
本発明は、センサの未設置区間であっても渋滞領域の起点、及び、末尾を特定する精度を向上させることができる。
【図面の簡単な説明】
【００２１】
【図１】第１の実施形態の渋滞推定装置の構成を示すブロック図である。
【図２】第１の実施形態の渋滞推定装置の動作を示すフローチャート図である。
【図３】第１の実施形態のセンサ情報の一例を示すデータシートである。
【図４】第１の実施形態のＧＰＳ情報の一例を示すデータシートである。
【図５】第１の実施形態におけるセンサ未設置区間の初期渋滞領域を示す説明図である。
【図６】第１の実施形態の渋滞情報の一例を示すデータシートである。
【図７】第１の実施形態における渋滞推定装置等をコンピュータで実現したハードウエア構成を示すブロック図である。
【発明を実施するための形態】
【００２２】
第１の実施形態の渋滞推定装置について、図面を用いて説明する。図１は、第１の実施形態の渋滞推定装置の構成を示すブロック図である。図１に示す渋滞推定装置１００は、情報収集部１０１、ＧＰＳ情報記憶部１０２、センサ情報記憶部１０３、道路情報記憶部１０４、渋滞起点・末尾特定部１０５、流入・流出量推定部１０６、渋滞領域更新部１０７、渋滞情報記憶部１０８、及び、情報提供部１０９を備える。
【００２３】
情報収集部１０１は、車両に搭載されたＧＰＳ車載器や携帯電話端末などから、ＧＰＳを用いて測定された交通情報（以降、「ＧＰＳ情報」とも称する）を収集し、ＧＰＳ情報記憶部１０２に記憶させる。ＧＰＳ情報は、例えば、走行車両のタイムスタンプ付の位置情報である。
【００２４】
情報収集部１１０はまた、道路に設置されているセンサで測定された交通情報（以降、「センサ情報」とも称する）を収集し、センサ情報記憶部１２３に記憶させる。センサ情報は、例えば、センサが設置された道路区間（以後、設置道路区間とも称す）で測定された車両速度、交通量、断面交通量、又は、交通密度等である。設置道路区間は、例えば、高速道路、交差点、又は、重要道路などである。
【００２５】
交通量は、道路のある地点を単位時間あたりに通過する車両台数（通過台数）を表す。断面交通量は、２地点の間に１つの断面を想定し、その断面を通過する交通量を表す。交通密度は、道路の単位区間に在る自動車等の車両台数、又は、空間・時間占有率を表す。
【００２６】
なお、情報収集部１０１が、直接ＧＰＳ搭載車両から車両位置情報を受信し、ＧＰＳ情報を生成してもよい。あるいは、情報収集部１０１が、外部装置（図示せず）によって生成されたＧＰＳ情報を収集して、ＧＰＳ情報記憶部１０２に記憶させてもよい。
【００２７】
さらに情報収集部１０１は、センサで測定された車両速度等をセンサから直接受信してもよく、センサに直結する一次処理装置などから直接収集してもよい。また、交通管制装置などの他の装置が収集・加工した情報を収集して記憶させてもよい。
【００２８】
道路情報記憶部１０４は、道路ネットワークにおける道路の接続情報、道路の区間距離など道路構造に関する情報を記憶する。
【００２９】
渋滞起点・末尾特定部１０５は、車両ごとに検出された位置情報を基に、渋滞領域を定めるための渋滞末尾、及び、渋滞起点を特定する。具体的には、ＧＰＳ情報およびセンサ情報に基づいて、初期渋滞領域の渋滞起点および渋滞末尾を特定する。渋滞起点・末尾特定部１０５は、特定した渋滞起点および渋滞末尾を記憶装置（図示せず）に記憶させる。なお、ＧＰＳ情報記憶部１０２又はセンサ情報記憶部１０３に渋滞起点又は渋滞末尾を記憶させてもよい。
【００３０】
流入・流出量推定部１０６は、渋滞末尾から上流側のセンサで測定された交通量と、渋滞起点から下流側のセンサで測定された交通量と、位置情報を基に、時刻ごとの渋滞末尾への流入量を推定し、時刻ごとの渋滞起点からの流出量を推定する。流入・流出量推定部１０６は、推定した流入量および流出量を、記憶装置（図示せず）に記憶させる。流入・流出量推定部１０６の詳細については、渋滞推定装置の動作にて具体的に説明する。
【００３１】
渋滞領域更新部１０７は、流入・流出量推定部１０６によって推定された流入量を基に渋滞末尾を更新し、推定された流出量を基に渋滞末尾起点を更新する。また、渋滞領域更新部１０７は、流入・流出量推定部１０６が推定した流入量と流出量との差を基に、渋滞領域の延伸を更新する。渋滞領域更新部１０７は、更新した渋滞領域の情報を渋滞情報記憶部１０８に記憶させる。
【００３２】
渋滞情報記憶部１０８は、渋滞領域更新部１０７が推定した渋滞の時間変化情報を記憶する。情報提供部１０９は、渋滞情報記憶部１０８に記憶された渋滞情報を外部装置（図示せず）に提供する。
【００３３】
次に、第１の実施形態の渋滞推定装置の動作について図面を用いて説明する。はじめに、渋滞推定装置における全体処理の動作について説明する。
【００３４】
図２は、第１の実施形態の渋滞推定装置の動作を示すフローチャートである。
【００３５】
渋滞起点・末尾特定部１０５は、車両ごとに検出された位置情報を基に、渋滞領域を定めるための渋滞末尾、及び、渋滞起点を特定する（ステップＳ２０１）。
【００３６】
流入・流出量推定部１０６は、渋滞末尾から上流側のセンサで測定された交通量と、渋滞起点から下流側のセンサで測定された交通量と、位置情報を基に、時刻ごとの渋滞末尾への流入量を推定し（ステップＳ２０２）、時刻ごとの渋滞起点からの流出量を推定する（ステップＳ２０３）。
【００３７】
渋滞領域更新部１０７は、推定した流入量、流出量と、センサ設置点の間における道路の車線数、及び、平均車両長に基づいて、渋滞末尾の延伸を計算し、渋滞情報を更新する（ステップＳ２０４）。この時、渋滞情報を更新した結果、渋滞が解消している場合（ステップＳ２０５のＹｅｓ）、処理を終了する。渋滞が解消していない場合（ステップＳ２０５のＮｏ）、ステップＳ２０２の処理を継続する。
【００３８】
以降に、各ステップの処理の詳細について説明する。
【００３９】
渋滞起点・末尾特定部１０５による、初期渋滞領域の渋滞末尾、渋滞起点を特定する処理（ステップＳ２０１）について説明する。
【００４０】
初期渋滞領域を特定する特定処理には、センサ情報とＧＰＳ情報を用いる特定処理、又は、ＧＰＳ情報を用いる特定処理がある。まず、センサ情報とＧＰＳ情報を用いる特定処理について説明する。
【００４１】
渋滞起点・末尾特定部１０５は、センサ情報記憶部１０３のセンサ情報を用いて、センサ情報によって検知可能な渋滞地点（以後、センサ渋滞地点と示す）を特定する。センサ渋滞地点は、例えば、センサで測定された車両の平均速度が一定以下となる地点、あるいは、交通密度が一定以上となる地点である。
【００４２】
図３は、センサ情報記憶部１０３に記憶されたセンサ情報の一例である。図３に示すセンサ情報には、センサ_ＩＤ、時刻、交通量、車両密度、地点速度が含まれる。センサＩＤは、道路に設置された固定センサを識別するための識別子である。時刻は、センサの測定時刻である。
【００４３】
交通量は、センサによってある時刻に測定された単位時間あたりの車両台数（センサを通過した台数）である。車両密度は、単位長さあたりの車両台数である。地点速度は、センサ設置点における車両の平均車両速度である。
【００４４】
図４は、ＧＰＳ情報記憶部１０２に記憶されたＧＰＳ情報の一例である。図４に示すＧＰＳ情報には、ＧＰＳ_ＣＡＲ_ＩＤ、タイムスタンプ、位置情報、速度が含まれる。ＧＰＳ_ＣＡＲ_ＩＤは、ＧＰＳ搭載車両を識別するための識別子である。タイムスタンプには、ＧＰＳ搭載車両が車両位置をＧＰＳ測位した日時が含まれる。位置情報は、ＧＰＳ測位時の車両位置である。速度は、ＧＰ搭載車両の車両速度である。なお、車両軌跡は、例えば、あるセンサ設置点間から次のセンサ設置点までの間におけるＧＰＳ搭載車両の位置情報の集合として表すことができる。
【００４５】
渋滞起点・末尾特定部１０５は、センサ渋滞地点が連続して出現する場合、その連続したセンサ渋滞地点の区間を組み合わせ、渋滞領域（センサ）として導出する。さらに渋滞起点・末尾特定部１０５は、この渋滞領域（センサ）を通過したＧＰＳ搭載車両の車両軌跡（軌跡情報）をＧＰＳ情報記憶部１０２より取得する。渋滞起点・末尾特定部１０５は、取得した車両軌跡を基に、ＧＰＳ搭載車両の平均速度が一定値以下となっている領域を渋滞領域（センサ＋ＧＰＳ）として導出する。
【００４６】
さらに、渋滞起点・末尾特定部１０５は、渋滞領域（センサ＋ＧＰＳ）より上流側あるいは下流側でＧＰＳ搭載車両の平均速度が一定値以下となっている領域があれば、当該領域を渋滞領域（センサ＋ＧＰＳ）に追加して更新する。なお、車両の走行方向は上流から下流に向かうものとする。
【００４７】
渋滞起点・末尾特定部１０５は、この更新後の渋滞領域（センサ＋ＧＰＳ）を初期渋滞領域として特定する。初期渋滞領域の端（ＧＰＳ搭載車両の平均速度が一定値以下となる位置）は、それぞれ初期渋滞領域の渋滞起点、及び、初期渋滞領域の渋滞末尾となる。
【００４８】
なお、センサ情報を用いずＧＰＳ情報を用いて初期渋滞領域を特定することも可能である。渋滞起点・末尾特定部１０５は、ＧＰＳ搭載車両の車両軌跡をＧＰＳ情報記憶部１０２より取得し、渋滞領域（センサ）を用いずに、取得した車両軌跡を基にＧＰＳ搭載車両の平均速度が一定値以下となっている領域を渋滞領域（ＧＰＳ）として導出する。渋滞起点・末尾特定部１０５は、渋滞領域（ＧＰＳ）を初期渋滞領域として特定する。初期渋滞領域の端は、それぞれ初期渋滞領域の渋滞起点、及び、初期渋滞領域の渋滞末尾となる。
【００４９】
次に、流入・流出量推定部１０６による、渋滞末尾への流入量の推定（ステップＳ２０２）、渋滞起点からの流出量の推定（ステップＳ２０３）について説明する。
【００５０】
流入・流出量推定部１０６は、各時刻に渋滞末尾に到着する単位時間あたりの到着台数を、渋滞末尾への流入量として推定する。また、各時刻に渋滞起点から離脱する単位時間あたりの離脱台数を、渋滞起点からの流出量として推定する。
【００５１】
図５は、第１の実施形態におけるセンサ未設置区間の初期渋滞領域を示す説明図である。図５の例では、上流センサ設置点と下流センサ設置点の間に初期渋滞領域が存在している。図中、上流センサ設置点から渋滞末尾までと下流センサ設置点から渋滞起点までの間は、未渋滞領域となっている。
【００５２】
次に、流入量の推定について説明する。渋滞末尾の近傍にある上流センサで測定された交通量（センサ設置点における単位時間あたりの車両台数）は、そのままの交通量で渋滞末尾に到着するとは限らない。上流センサ設置点と渋滞末尾との距離が長いほど、上流センサ（上流側のセンサとも称する）を通過した車両の速度変化の影響を受けるようになり、上流センサ設置点から渋滞末尾までの距離に応じて車両が渋滞末尾に到着する割合（到着率）が異なることになる。
【００５３】
また、前回、上流センサで測定されて渋滞末尾に到着しなかった残留車両が、渋滞末尾に到着することも考慮する必要がある。
【００５４】
このため、流入・流出量推定部１０６は、下記の式（１）によって各時刻における渋滞末尾への流入量を推定する。
2018180688.tif 000003
ここで、
2018180688.tif 000004
は、時刻ｔにおける渋滞末尾への流入量（単位時間あたりの到着台数）、
2018180688.tif 000005
は、時刻tに上流センサで測定された交通量（単位時間あたりのセンサ通過台数）、
2018180688.tif 000006
は、時刻ｔ−１における上流センサ設置点から渋滞末尾までの距離、
ｖは、上流センサ設置点から渋滞末尾まで走行したＧＰＳ搭載車両の平均速度、
2018180688.tif 000007
は、時刻ｔ−１に上流センサを通過した車両が、渋滞末尾に到達するまでの所要時間、
Δｔは、上流センサの計測時間の幅を示す。
【００５５】
また、式（１）の第二項の
2018180688.tif 000008
は、時刻tより前に上流センサで測定された交通量（単位時間あたりのセンサ通過台数）のうち、時刻tに渋滞末尾に到着する流入量（単位時間あたりの到着台数）を示す。
【００５６】
式（１）において、ある計測時間において、計測開始から当該所要時間より後に上流センサを通過した車両は、時刻ｔの計測時間の幅内に渋滞末尾に辿りつくことはできない。そのため、上流センサによる計測時間の幅からこれを引き、上流センサの計測時間の幅で除算することで到着率を求めている。時刻ｔに上流センサで測定された交通量に当該到着率を乗算することで渋滞末尾に到達する到達流入量（単位時間あたりの到着台数）を推定することができる。
【００５７】
また、式（１）の
2018180688.tif 000009
は、下記の式（２）で計算される。
2018180688.tif 000010
上記式（２）の
2018180688.tif 000011
は、時間ステップｔ−ｉにおいて未渋滞領域に残留した残留台数を示している。
【００５８】
その次の項は、時刻tからi遡った時間において未渋滞領域に残留した車両のうち、時刻tに渋滞末尾へ到達可能な車両の割合を示している。分母は未渋滞領域に残留した車両が使用可能な時間幅、分子はその使用可能な時間幅から、渋滞末尾まで走行した車両の平均時間を引いたものを示す。
【００５９】
流入・流出量推定部１０６は、式（１）、式（２）から各時刻に渋滞末尾に到着する単位時間あたりの到着台数を、渋滞末尾への流入量として推定する。流入・流出量推定部１０６は、推定した渋滞末尾への流入量を記憶装置（図示せず）に記憶させる。
【００６０】
次に、流出量の推定について説明する。センサ未設置区間に存在する渋滞起点からの流出量についても同様に近傍の下流センサ（下流側のセンサとも称する）で測定された交通量をそのまま使用することはできない。渋滞起点から下流センサ設置点までの距離に応じて車両が渋滞起点から離脱する割合（離脱率）が異なるからである。このため、流入・流出量推定部１０６は、式（３）によって時刻ｔにおける渋滞起点からの流出量（単位時間あたりの離脱台数）を推定する。
2018180688.tif 000012
ここで、
2018180688.tif 000013
は、時刻ｔにおける渋滞起点からの流出量（単位時間あたりの離脱台数）、
2018180688.tif 000014
は、時刻ｔに下流センサで測定された交通量（単位時間あたりのセンサ通過台数）、
2018180688.tif 000015
は、時刻ｔ−１における渋滞起点から下流センサまでの距離、
2018180688.tif 000016
は、時刻ｔ−１に渋滞起点を離脱した車両が下流センサに到達するまでの所要時間、
ｖは、渋滞起点から下流センサまで走行したＧＰＳ搭載車両の平均速度、
Δｔは、下流センサの計測時間の幅を示す。
【００６１】
式（３）の第一項は、時刻ｔに下流センサで測定された交通量（単位時間あたりのセンサ通過台数）に、渋滞起点から離脱してきた車両の離脱率を乗算することで渋滞起点から離脱する離脱交通量（単位時間あたりの離脱台数）を推定している。
【００６２】
上記式（３）の第二項の
2018180688.tif 000017
は、時刻ｔに渋滞起点を離脱し、時刻ｔより後に下流センサで測定された交通量（単位時間あたりのセンサ通過台数）を表している。
【００６３】
また、式（３）の
2018180688.tif 000018
は下記の式で計算される。
2018180688.tif 000019
上記式の
2018180688.tif 000020
は、時間ステップｔ＋ｉにおいて未渋滞領域に残留した残留台数を示している。
【００６４】
その次の項は、時刻tからi経過した時間において未渋滞領域へ離脱した車両のうち、時刻t＋１に下流センサに到着可能な車両の割合を示している。この項に上述の残量台数を乗じた値が離脱残留交通量となる。分母は未渋滞領域に残留した車両が使用可能な時間幅、分子はその使用可能な時間幅から下流センサ走破するのにかかる平均時間を引いたものを示す。
【００６５】
流入・流出量推定部１０６は、式（３）、式（４）から各時刻に渋滞起点から離脱する単位時間あたりの離脱台数を、渋滞起点からの流出量として推定する。流入・流出量推定部１０６は、推定した渋滞末尾への流入量を記憶装置（図示せず）に記憶させる。
【００６６】
渋滞領域更新部１０７による、初期渋滞領域の渋滞末尾の更新、渋滞起点を更新する処理（ステップＳ２０４）について説明する。
【００６７】
渋滞領域更新部１０７は、流入・流出量推定部１０６が推定した流入量を基に初期渋滞領域の渋滞末尾を更新し、推定した流出量を基に、初期渋滞領域の渋滞起点を更新する。具体的には、渋滞領域更新部１０７は、時間当たりの到着台数と離脱台数のそれぞれに対し、車線数で除算すると車線あたりの到着台数、離脱台数の増減分が算出され、この値に平均車両長＋平均車間距離を乗算することで、渋滞末尾、渋滞起点のそれぞれにおける単位時間あたりの渋滞長の延伸が計算できる。
【００６８】
また、渋滞領域更新部１０７は、流入・流出量推定部１０６が推定した流入量と流出量との差を基に初期渋滞領域の時間変化を推定する。具体的には、渋滞領域更新部１０７は、時間当たりの到着台数と離脱台数の差分を求める。この差分を車線数で除算すると車線あたりの滞留台数の増減分が算出され、この値に平均車両長＋平均車間距離を乗算することで、初期渋滞領域の渋滞長の延伸が計算できる。
【００６９】
渋滞推定装置１００は、このようにして計算した渋滞領域を渋滞情報記憶部１０８に保存する。図６は、渋滞情報の一例を示すデータシートである。図６に示す渋滞情報は、渋滞ＩＤ、渋滞起点、渋滞末尾、渋滞長、時刻を含む。渋滞ＩＤは、渋滞を識別するための識別子である。渋滞起点は、渋滞ＩＤの渋滞領域における渋滞の起点位置である。渋滞末尾は、渋滞ＩＤの渋滞領域における渋滞の末尾位置である。渋滞長は、渋滞起点と渋滞末尾で規定される渋滞領域の長さである。時刻は、渋滞長を推定した時刻である。
【００７０】
渋滞推定装置１００は、例えば、道路が高速道路である場合、渋滞起点、渋滞末尾を当該高速道路に設置された「キロポスト」情報として渋滞情報記憶部１０８に記憶させてもよい。
【００７１】
情報提供部１０９は、渋滞情報記憶部１０８の更新された渋滞領域を取得して、ユーザへ提供する。
【００７２】
本実施形態によれば、センサの未設置区間であっても渋滞区間の起点、及び、末尾を特定する精度を向上させることができる。渋滞推定装置の渋滞起点・末尾特定部は、車両ごとに検出された位置情報を基に、渋滞領域を定めるための渋滞末尾、及び、渋滞起点を特定する。さらに、流入・流出量推定部は、渋滞末尾から上流側のセンサで測定された交通量と、渋滞起点から下流側のセンサで測定された交通量と、位置情報を基に、時刻ごとの渋滞末尾への流入量を推定し、時刻ごとの渋滞起点からの流出量を推定する。流入・流出量推定部は、推定された流入量を基に渋滞末尾を更新し、推定された流出量を基に渋滞末尾起点を更新するからである。
【００７３】
上記実施形態における道路は、高速道路又は山間部の道路に適用可能である。
【００７４】
なお、ＧＰＳ情報記憶部１０２に保存されるＧＰＳ情報は、どの時間帯にどの区間を車両が走行したかを示す走行履歴であればよく、ＧＰＳの観測結果に基づいて取得されることに限定されない。例えば、車載器のＭＡＣ（Ｍｅｄｉａ  Ａｃｃｅｓｓ  Ｃｏｎｔｒｏｌ）アドレスを用いて車両を追跡した結果を取得してもよいし、あるいは、通信端末にユーザから入力された情報から取得してもよいし、ユーザから地点と時間を示す情報を何らかの手法で取得してもよい。
（ハードウエア構成）
図７は、第１の実施形態の渋滞推定装置をコンピュータで実現したハードウエア構成を示す図である。第１の実施形態の渋滞推定装置の各構成要素は、機能単位のブロックを示している。各構成要素の一部又は全部は、例えば、図７に示すようなコンピュータ６００とプログラムとの任意の組み合わせにより実現される。コンピュータ６００は、一例として、以下のような構成を含む。
ＣＰＵ（Ｃｅｎｔｒａｌ  Ｐｒｏｃｅｓｓｉｎｇ  Ｕｎｉｔ）６０１、
ＲＯＭ（Ｒｅａｄ  Ｏｎｌｙ  Ｍｅｍｏｒｙ）６０２、
ＲＡＭ（Ｒａｎｄｏｍ  Ａｃｃｅｓｓ  Ｍｅｍｏｒｙ）６０３、
ＲＡＭ６０３にロードされるプログラム６０４、
プログラム６０４を格納する記憶装置６０５、
記録媒体６０６の読み書きを行うドライブ装置６０７、
通信ネットワーク６０９と接続する通信インターフェース６０８、
データの入出力を行う入出力インターフェース６１０、
各構成要素を接続するバス６１１
渋滞推定装置の各構成要素は、これらの機能を実現するプログラム６０４をＣＰＵ６０１が取得して実行することで実現される。各構成要素の機能を実現するプログラム６０４は、例えば、予め記憶装置６０５又はＲＡＭ６０３に格納されており、必要に応じてＣＰＵ６０１が読み出す。なお、プログラム６０４は、通信ネットワーク６０９を介してＣＰＵ６０１に供給されてもよいし、予め記録媒体６０６に格納されており、ドライブ装置６０７が当該プログラムを読み出してＣＰＵ６０１に供給してもよい。
【００７５】
渋滞推定装置の実現方法には、様々な変形例がある。例えば、渋滞推定装置は、各構成要素にそれぞれ別個のコンピュータ６００とプログラムとの任意の組み合わせにより実現されてもよい。また、渋滞推定装置が備える複数の構成要素が、一つのコンピュータ６００とプログラムとの任意の組み合わせにより実現されてもよい。
【００７６】
また、渋滞推定装置の各構成要素の一部又は全部は、その他の汎用または専用の回路、プロセッサ等やこれらの組み合わせによって実現される。これらは、単一のチップによって構成されてもよいし、バスを介して接続される複数のチップによって構成されてもよい。また、コンピュータ６００の代わりにＦＰＧＡ（Ｆｉｅｌｄ−Ｐｒｏｇｒａｍｍａｂｌｅ  Ｇａｔｅ  Ａｒｒａｙ）のようなプログラマブルロジックデバイスを用いてもよい。
【００７７】
さらに、渋滞推定装置の各構成要素の一部又は全部は、上述した回路等とプログラムとの組み合わせによって実現されてもよい。
【００７８】
また、渋滞推定装置の各構成要素の一部又は全部が複数のコンピュータや回路等により実現される場合には、複数のコンピュータや回路等は、集中配置されてもよいし、分散配置されてもよい。例えば、コンピュータや回路等は、クライアントアンドサーバシステム、クラウドコンピューティングシステム等、各々が通信ネットワークを介して接続される形態として実現されてもよい。
【００７９】
以上、実施形態を参照して本願発明を説明したが、本願発明は上記実施形態に限定されるものではない。本願発明の構成や詳細には、本願発明のスコープ内で当業者が理解し得る様々な変更をすることができる。
【００８０】
また、図面中の矢印の方向は、一例を示すものであり、ブロック間の信号の向きを限定するものではない。
【００８１】
この出願は、２０１７年３月３１日に出願された日本出願特願２０１７−０７０４５１を基礎とする優先権を主張し、その開示の全てをここに取り込む。
【符号の説明】
【００８２】
１００  渋滞推定装置
１０１  情報収集部
１０２  ＧＰＳ情報記憶部
１０３  センサ情報記憶部
１０４  道路情報記憶部
１０５  渋滞起点・末尾特定部
１０６  流入・流出量推定部
１０７  渋滞領域更新部
１０８  渋滞情報記憶部
１０９  情報提供部
６００  コンピュータ
６０１  ＣＰＵ
６０２  ＲＯＭ（Ｒｅａｄ  Ｏｎｌｙ  Ｍｅｍｏｒｙ）
６０３  ＲＡＭ
６０４  プログラム
６０５  記憶装置
６０６  記録媒体
６０７  ドライブ装置
６０８  通信インターフェース
６０９  通信ネットワーク
６１０  入出力インターフェース
６１１  バス
【図１】
2018180688.tif 000021
【図２】
2018180688.tif 000022
【図３】
2018180688.tif 000023
【図４】
2018180688.tif 000024
【図５】
2018180688.tif 000025
【図６】
2018180688.tif 000026
【図７】
2018180688.tif 000027
【国際調査報告】
2018180688.tif 000028
2018180688.tif 000029