(57)【特許請求の範囲】
【請求項１】
コンピュータを、
シミュレーション時刻を時々刻々と計時する計時手段、
目的駅と移動経路とを設定した仮想旅客を所定駅に時系列に発生させるとともに、前記計時されているシミュレーション時刻に従って所定のダイヤデータに沿った列車が運行されているとして当該列車に前記仮想旅客を乗車させて、前記所定駅での列車待ち旅客１人１人を管理する列車待ち旅客管理手段、
として機能させるためのプログラムであって、
前記列車待ち旅客管理手段が、前記ダイヤデータのうち前記仮想駅を含む仮想列車の運行を中断させる中断区間を設定する運行中断設定手段を有し、当該運行中断区間における列車が運行されないこととして仮想旅客を管理し、
前記列車待ち旅客管理手段により管理されている仮想旅客１人１人について、前記運行中断設定手段による中断後の運転再開想定時刻を設定する再開想定時刻設定手段、
前記計時手段により計時されているシミュレーション時刻が予め定められた所定時刻に到達した際に、前記列車待ち旅客管理手段により管理されている仮想旅客１人１人について、１）前記運行中断設定手段により中断が設定されなかった場合に、設定されている移動経路で設定されている目的駅までの移動に要する通常時の到着所要推定時間と、２）設定されている運転再開想定時刻までの待ち時間と、３）前記運行中断区間を迂回する迂回経路で設定されている目的駅までの移動に要する迂回時の到着所要推定時間とを用いた所定の確率算出処理を実行して前記迂回経路の選択確率を算出し、算出された確率に従って当該仮想旅客の経路を選択して当該仮想旅客の移動経路を更新する移動経路選択更新手段、
前記仮想旅客１人１人に設定されている運転再開想定時刻を可変する運転再開想定時刻可変手段、
として更に前記コンピュータを機能させるためのプログラム。
【請求項２】
前記移動経路選択更新手段が、１）通常時の到着所要推定時間と待ち時間の合算時間とを全所要時間とした場合の通常時の効用値と、２）迂回時の到着所要推定時間を全所要時間とした場合の迂回時の効用値とのそれぞれについて、全所要時間を変数とする所定の効用値算出演算を行って算出するとともに、それぞれの場合の効用値を変数とする所定の確率算出演算を行って前記迂回経路の選択確率を算出するように前記コンピュータを機能させるための請求項１に記載のプログラム。
【請求項３】
前記移動経路選択更新手段が、前記算出された迂回時の効用値が通常時の効用値より所定の閾値以上又は超過した場合に前記迂回経路の選択確率を１／２以上とするように前記迂回経路の選択確率を可変して算出するように前記コンピュータを機能させるための請求項２に記載のプログラム。
【請求項４】
前記計時手段により計時されているシミュレーション時刻が、運転再開予定時刻を仮想的に報知する時刻として定められた所定のシミュレーション時刻に到達したことを検出する報知時刻到達検出手段、
前記報知時刻到達検出手段による検出がなされた場合に、前記仮想旅客それぞれに設定されている運転再開想定時刻を前記運転再開予定時刻に更新する運転再開時刻更新手段、
として前記コンピュータを機能させるとともに、
前記移動経路選択更新手段が、前記報知時刻到達検出手段による検出がなされた場合に、前記仮想旅客１人１人についての前記迂回経路の選択確率の算出及び算出された確率に従った経路選択を実行するように前記コンピュータを機能させるための請求項１〜３の何れか一項に記載のプログラム。
【請求項５】
前記列車待ち旅客管理手段により管理されている仮想旅客１人１人について、前記運行中断区間を迂回する迂回経路の認知有無を初期設定する迂回経路認知有無初期設定手段として前記コンピュータを機能させるとともに、
前記移動経路選択更新手段が、前記列車待ち旅客管理手段により管理されている仮想旅客のうち、迂回経路の認知有りと設定されている仮想旅客のみを対象として前記迂回経路の選択確率の算出及び算出された確率に従った経路選択を実行するように前記コンピュータを機能させ、
さらに、
前記計時手段により計時されているシミュレーション時刻が、迂回経路を仮想的に報知する時刻として定められた所定のシミュレーション時刻に到達したことを検出する迂回経路報知時刻到達検出手段、
前記迂回経路報知時刻到達検出手段による検出がなされた場合に、迂回経路の認知無しと設定されている仮想旅客を認知有りに更新する迂回経路認知更新手段、
として前記コンピュータを機能させ、
前記移動経路選択更新手段が、前記迂回経路報知到達時刻検出手段による検出がなされた場合に、前記仮想旅客１人１人についての前記迂回経路の選択確率の算出及び算出された確率に従った経路選択を実行するように前記コンピュータを機能させる、
ための請求項１〜４の何れか一項に記載のプログラム。
【請求項６】
シミュレーション時刻を時々刻々と計時する計時手段と、
目的駅と移動経路とを設定した仮想旅客を所定駅に時系列に発生させるとともに、前記計時されているシミュレーション時刻に従って所定のダイヤデータに沿った列車が運行されているとして当該列車に前記仮想旅客を乗車させて、前記所定駅での列車待ち旅客１人１人を管理する列車待ち旅客管理手段と、
を備えたシミュレーション装置であって、
前記列車待ち旅客管理手段が、前記ダイヤデータのうち前記仮想駅を含む仮想列車の運行を中断させる中断区間を設定する運行中断設定手段を有し、当該運行中断区間における列車が運行されないこととして仮想旅客を管理し、
前記列車待ち旅客管理手段により管理されている仮想旅客１人１人について、前記運行中断設定手段による中断後の運転再開想定時刻を設定する再開想定時刻設定手段と、
前記計時手段により計時されているシミュレーション時刻が予め定められた所定時刻に到達した際に、前記列車待ち旅客管理手段により管理されている仮想旅客１人１人について、１）前記運行中断設定手段により中断が設定されなかった場合に、設定されている移動経路で設定されている目的駅までの移動に要する通常時の到着所要推定時間と、２）設定されている運転再開想定時刻までの待ち時間と、３）前記運行中断区間を迂回する迂回経路で設定されている目的駅までの移動に要する迂回時の到着所要推定時間とを用いた所定の確率算出処理を実行して前記迂回経路の選択確率を算出し、算出された確率に従って当該仮想旅客の経路を選択して当該仮想旅客の移動経路を更新する移動経路選択更新手段と、
前記仮想旅客１人１人に設定されている運転再開想定時刻を可変する運転再開想定時刻可変手段と、
を備えたシミュレーション装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、プログラム及びシミュレーション装置に関する。
【背景技術】
【０００２】
鉄道において、事故や故障、災害等によって列車の運行に乱れた生じた場合、安全な運行を遂行するために列車ダイヤを変更する運転整理が行われるが、運転整理案の作成の際には、旅客需要をより正確に予測する必要がある。そこで、旅客の行動を予測するための旅客行動モデルとして、列車運行と旅客数との動的な相互作用を考慮したモデルが提案されている。この旅客行動モデルでは、旅客が列車に乗車するか否かの判断や、列車に乗車している旅客が降車するか否かの判断を、ロジットモデルにもとづいて行う。そして、作成した運転整理案をこの旅客行動モデルを用いて評価し、より良い運転整理案の作成を支援することを可能としている（例えば、特許文献１参照）。
【特許文献１】特開２００１−５５１４５号公報
【発明の開示】
【発明が解決しようとする課題】
【０００３】
しかしながら、上述の特許文献１の旅客行動モデルは、乗車するか否かの乗車と降車するか否かの乗車行動とをロジットモデルを用いて判定する手法のみが開示されている。ダイヤ乱れ時には、旅客は単純に乗車／降車するか否かといった判断以外の他の行動を取る可能性も存在し、その判断には複雑な要因がある。具体的には、旅客は、通常時には、自らの行動を、計画ダイヤに従った運行がなされている前提で決定するが、ダイヤ乱れ時には、列車の運休や遅れ、臨時列車の運行といった列車の運行状況を把握しきれないため、駅員による見込みの情報や自身の予測といった曖昧な情報をもとに決定する。迂回経路を選択するといった行動も、その迂回経路の存在の有無やその迂回経路で目的駅まで移動した場合の所要時間も要因の一つとなってくる。更には、旅客は１人１人目的駅や出現時刻が異なるため、旅客１人１人の行動を推定することが、精度向上に不可欠な技術といえる。本発明は、上記事情に鑑みてなされたものであり、ダイヤ乱れ時の旅客の行動を精度良く予測することを目的としている。
【課題を解決するための手段】
【０００４】
上記課題を解決するための第１の発明は、
コンピュータ（例えば、図３のシミュレーション装置１）を、
シミュレーション時刻を時々刻々と計時する計時手段（例えば、図９のステップＡ３３）、
目的駅と移動経路とを設定した仮想旅客を所定駅に時系列に発生させるとともに、前記計時されているシミュレーション時刻に従って所定のダイヤデータに沿った列車が運行されているとして当該列車に前記仮想旅客を乗車させて、前記所定駅での列車待ち旅客１人１人を管理する列車待ち旅客管理手段（例えば、図９のステップＡ５〜Ａ３１）、
として機能させるためのプログラム（例えば、図３のシミュレーションプログラム５１０）であって、
前記列車待ち旅客管理手段が、前記ダイヤデータのうち前記仮想駅を含む仮想列車の運行を中断させる中断区間を設定する運行中断設定手段（例えば、図９のステップＡ１）を有し、当該運行中断区間における列車が運行されないこととして仮想旅客を管理し、
前記列車待ち旅客管理手段により管理されている仮想旅客１人１人について、前記運行中断設定手段による中断後の運転再開想定時刻を設定する再開想定時刻設定手段（例えば、図９のステップＡ５）、
前記計時手段により計時されているシミュレーション時刻が予め定められた所定時刻に到達した際に、前記列車待ち旅客管理手段により管理されている仮想旅客１人１人について、１）前記運行中断設定手段により中断が設定されなかった場合に、設定されている移動経路で設定されている目的駅までの移動に要する通常時の到着所要推定時間と、２）設定されている運転再開想定時刻までの待ち時間と、３）前記運行中断区間を迂回する迂回経路で設定されている目的駅までの移動に要する迂回時の到着所要推定時間とを用いた所定の確率算出処理を実行して前記迂回経路の選択確率を算出し、算出された確率に従って当該仮想旅客の経路を選択して当該仮想旅客の移動経路を更新する移動経路選択更新手段（例えば、図９のステップＡ２１）、
前記仮想旅客１人１人に設定されている運転再開想定時刻を可変する運転再開想定時刻可変手段（例えば、図９のステップＡ１７）、
として更に前記コンピュータを機能させるためのプログラムである。
【０００５】
また、第６の発明は、
シミュレーション時刻を時々刻々と計時する計時手段（例えば、図３のＣＰＵ１０）と、
目的駅と移動経路とを設定した仮想旅客を所定駅に時系列に発生させるとともに、前記計時されているシミュレーション時刻に従って所定のダイヤデータに沿った列車が運行されているとして当該列車に前記仮想旅客を乗車させて、前記所定駅での列車待ち旅客１人１人を管理する列車待ち旅客管理手段（例えば、図３のＣＰＵ１０）と、
を備えたシミュレーション装置（例えば、図３のシミュレーション装置１）であって、
前記列車待ち旅客管理手段が、前記ダイヤデータのうち前記仮想駅を含む仮想列車の運行を中断させる中断区間を設定する運行中断設定手段を有し、当該運行中断区間における列車が運行されないこととして仮想旅客を管理し、
前記列車待ち旅客管理手段により管理されている仮想旅客１人１人について、前記運行中断設定手段による中断後の運転再開想定時刻を設定する再開想定時刻設定手段（例えば、図３のＣＰＵ１０）と、
前記計時手段により計時されているシミュレーション時刻が予め定められた所定時刻に到達した際に、前記列車待ち旅客管理手段により管理されている仮想旅客１人１人について、１）前記運行中断設定手段により中断が設定されなかった場合に、設定されている移動経路で設定されている目的駅までの移動に要する通常時の到着所要推定時間と、２）設定されている運転再開想定時刻までの待ち時間と、３）前記運行中断区間を迂回する迂回経路で設定されている目的駅までの移動に要する迂回時の到着所要推定時間とを用いた所定の確率算出処理を実行して前記迂回経路の選択確率を算出し、算出された確率に従って当該仮想旅客の経路を選択して当該仮想旅客の移動経路を更新する移動経路選択更新手段（例えば、図３のＣＰＵ１０）と、
前記仮想旅客１人１人に設定されている運転再開想定時刻を可変する運転再開想定時刻可変手段（例えば、図３のＣＰＵ１０）と、
を備えたシミュレーション装置である。
【０００６】
この第１又は第６の発明によれば、仮想旅客を所定駅に時系列に発生させるとともに、シミュレーション時刻に従って所定のダイヤデータに沿って運行される列車に仮想旅客を乗車させて、所定駅での列車待ち旅客１人１人を管理するシミュレーションにおいて、対象駅を含む仮想列車の運行を中断させる中断区間が設定され、設定された運行中断区間における列車が運行されないこととして仮想旅客が管理される。そして、シミュレーション時刻が予め定められた所定時刻に達すると、所定駅での列車待ち旅客１人１人について、１）通常時の到着所要推定時間と、２）運転再開想定時刻までの待ち時間と、３）迂回時の到着所要推定時間とを用いた所定の確率算出処理を実行して迂回経路の選択確率が算出され、算出された選択確率に従ってこの仮想旅客の経路が選択されて移動経路が更新される。つまり、所定駅で列車待ちの旅客それぞれについて、移動経路として、そのまま所定駅で待つか、運行中断区間を迂回するかの何れかが選択される。また、仮想旅客毎に設定される運転再開想定時刻が異なるため、算出される迂回経路の選択確率は仮想旅客毎に異なる。このように、仮想旅客１人１人について迂回経路の選択確率を算出して移動経路を更新することで、仮想旅客１人１人の行動を忠実に再現した精度の良いシミュレーションが可能となり、その結果、運行中断といったダイヤ乱れ時の旅客需要をより正確に予測することが可能となる。
【０００７】
第２の発明は、第１の発明のプログラムであって、
前記移動経路選択更新手段が、１）通常時の到着所要推定時間と待ち時間の合算時間とを全所要時間とした場合の通常時の効用値と、２）迂回時の到着所要推定時間を全所要時間とした場合の迂回時の効用値とのそれぞれについて、全所要時間を変数とする所定の効用値算出演算を行って算出するとともに、それぞれの場合の効用値を変数とする所定の確率算出演算を行って前記迂回経路の選択確率を算出するように前記コンピュータを機能させるためのプログラムである。
【０００８】
この第２の発明によれば、１）通常時の効用値と、２）迂回時の効用値とのそれぞれが、全所要時間を変数とする所定の効用値算出演算を行って算出されるとともに、それぞれの場合の効用値を変数とする所定の確率算出処理を行って、仮想旅客それぞれの迂回経路の選択確率が算出される。
【０００９】
また、第３の発明は、第２の発明のプログラムであって、
前記移動経路選択更新手段が、前記算出された迂回時の効用値が通常時の効用値より所定の閾値以上又は超過した場合に前記迂回経路の選択確率を１／２以上とするように前記迂回経路の選択確率を可変して算出するように前記コンピュータを機能させるためのプログラムである。
【００１０】
この第３の発明によれば、迂回時の効用値が通常時の効用値より所定の閾値以上又は超過した場合に迂回経路の選択確率を１／２以上とするように、迂回経路の選択確率を可変して算出される。運行中断時の旅客の行動としては、通常時の効用値と迂回時の効用値とがほぼ等しい場合、迂回するよりもそのまま所定駅で運転再開を待つ確率の方が高い。このため、迂回時の効用値と通常時の効用値との間に所定の閾値を設けることで、より現実に即したシミュレーションが実現される。
【００１１】
第４の発明は、第１〜第３の何れかの発明のプログラムであって、
前記計時手段により計時されているシミュレーション時刻が、運転再開予定時刻を仮想的に報知する時刻として定められた所定のシミュレーション時刻に到達したことを検出する報知時刻到達検出手段（例えば、図３のＣＰＵ１０；図９のステップＡ１５）、
前記報知時刻到達検出手段による検出がなされた場合に、前記仮想旅客それぞれに設定されている運転再開想定時刻を前記運転再開予定時刻に更新する運転再開時刻更新手段（例えば、図３のＣＰＵ１０；図９のステップＡ１７）、
として前記コンピュータを機能させるとともに、
前記移動経路選択更新手段が、前記報知時刻到達検出手段による検出がなされた場合に、前記仮想旅客１人１人についての前記迂回経路の選択確率の算出及び算出された確率に従った経路選択を実行するように前記コンピュータを機能させるためのプログラムである。
【００１２】
この第４の発明によれば、シミュレーション時刻が、運転再開予定時刻を仮想的に報知する時刻として定められた所定のシミュレーション時刻に到達したことが検出されると、仮想旅客それぞれに設定されている運転再開想定時刻がこの運転再開予定時刻に更新されるとともに、仮想旅客１人１人について、迂回経路の選択確率の算出及び算出された確率に従った経路選択が実行される。旅客の立場から見れば、運転再開時刻を知らない場合、自身が予測している運転再開想定時刻をもとに、迂回か待ちかのどちらかを選択する。そして、駅構内のアナウンス等で運転再開予定時刻が報知されると、この運転再開予定時刻をもとに、迂回か待ちかのどちらかを、再度選択する。つまり、運転再開予定時刻を知っているか否かによって迂回経路の選択確率が異なり、旅客自身が予測している運転再開想定時刻と運転再開予定時刻とのずれが大きいほど、この選択確率のずれも大きくなる。このため、運転再開予定時刻の報知がなされると、各仮想旅客の運転再開想定時刻をこの運転再開予定時刻に更新し、再度、迂回経路の選択確率を算出して経路選択を行うことで、より現実に即した精度の良いシミュレーションが可能となる。
【００１３】
第５の発明は、第１〜第４の何れかの発明のプログラムであって、
前記列車待ち旅客管理手段により管理されている仮想旅客１人１人について、前記運行中断区間を迂回する迂回経路の認知有無を初期設定する迂回経路認知有無初期設定手段（例えば、図３のＣＰＵ１０；図９のステップＡ５）として前記コンピュータを機能させるとともに、
前記移動経路選択更新手段が、前記列車待ち旅客管理手段により管理されている仮想旅客のうち、迂回経路の認知有りと設定されている仮想旅客のみを対象として前記迂回経路の選択確率の算出及び算出された確率に従った経路選択を実行するように前記コンピュータを機能させ、
さらに、
前記計時手段により計時されているシミュレーション時刻が、迂回経路を仮想的に報知する時刻として定められた所定のシミュレーション時刻に到達したことを検出する迂回経路報知時刻到達検出手段（例えば、図３のＣＰＵ１０；図９のステップＡ１１）、
前記迂回経路報知時刻到達検出手段による検出がなされた場合に、迂回経路の認知無しと設定されている仮想旅客を認知有りに更新する迂回経路認知更新手段（例えば、図３のＣＰＵ１０；図９のステップＡ１３）、
として前記コンピュータを機能させ、
前記移動経路選択更新手段が、前記迂回経路報知到達時刻検出手段による検出がなされた場合に、前記仮想旅客１人１人についての前記迂回経路の選択確率の算出及び算出された確率に従った経路選択を実行するように前記コンピュータを機能させる、
ためのプログラムである。
【００１４】
この第５の発明によれば、仮想旅客のうち、迂回経路の認知有りと設定されている仮想旅客のみを対象として、迂回経路の選択確率の算出及び算出された確率に従った経路選択が実行される。そして、計時されているシミュレーション時刻が、迂回経路を仮想的に報知する時刻として定められた所定のシミュレーション時刻に到達したことが検出されると、迂回経路の認知無しと設定されている仮想旅客が認知有りに更新されるとともに、仮想旅客１人１人についての迂回経路の選択確率の算出及び算出された確率に従った経路選択が実行される。旅客にとっては、迂回経路を知らない場合、必然的に運転再開を待つことになり、迂回経路が報知されると、その時点で初めて迂回か待ちかを選択することになる。このため、迂回経路の報知がなされると、迂回経路の認知無しと設定されている旅客について認知有りに更新し、これらの旅客それぞれについて迂回経路の選択確率を算出して経路選択を行うことで、より現実に即した精度の良いシミュレーションが可能となる。
【発明の効果】
【００１５】
本発明によれば、仮想旅客を所定駅に時系列に発生させるとともに、シミュレーション時刻に従って所定のダイヤデータに沿って運行される列車に仮想旅客を乗車させて、所定駅での列車待ち旅客１人１人を管理するシミュレーションにおいて、対象駅を含む仮想列車の運行を中断させる中断区間が設定され、設定された運行中断区間における列車が運行されないこととして仮想旅客が管理される。そして、シミュレーション時刻が予め定められた所定時刻に達すると、所定駅での列車待ち旅客１人１人について、１）通常時の到着所要推定時間と、２）運転再開想定時刻までの待ち時間と、３）迂回時の到着所要推定時間とを用いた所定の確率算出処理を実行して迂回経路の選択確率が算出され、算出された選択確率に従ってこの仮想旅客の経路が選択されて移動経路が更新される。つまり、所定駅で列車待ちの旅客それぞれについて、移動経路として、そのまま所定駅で待つか、運行中断区間を迂回するかの何れかが選択される。また、仮想旅客毎に設定される運転再開想定時刻が異なるため、算出される迂回経路の選択確率は仮想旅客毎に異なる。このように、仮想旅客１人１人について迂回経路の選択確率を算出して移動経路を更新することで、仮想旅客１人１人の行動を忠実に再現した精度の良いシミュレーションが可能となり、その結果、運行中断といったダイヤ乱れ時の旅客需要をより正確に予測することが可能となる。
【発明を実施するための最良の形態】
【００１６】
以下、図面を参照して、本発明の好適な実施形態を説明する。
【００１７】
［原理］
先ず、本実施形態の原理を説明する。鉄道において、ダイヤ乱れ時の旅客の選択行動としては、大きくは、（１）運転再開を待つ、（２）迂回経路を選択する、の何れかとなる。これは、旅客が予定していた「移動経路」、或いはこの移動経路を迂回する「迂回経路」の何れかを選択する二者択一の経路選択モデルとして扱うことができる。
【００１８】
図１に、本実施形態の経路選択モデルの概要を示す。同図では、対象駅である「Ａ駅」に出現した旅客Ｒが、目的駅である「Ｂ駅」に向うケースを示している。「Ａ駅」から「Ｂ駅」に向う経路として、旅客Ｒが予定している「移動経路」と、この移動経路を迂回して対象駅から目的駅に向う経路である「迂回経路」とがある。そして、旅客Ｒの移動経路に相当する「Ａ駅」と「Ｂ駅」との間の区間の列車運行が中断された場合、この旅客Ｒは、（１）運転が再開された後、予定していた「移動経路」でＢ駅に向う、（２）「迂回経路」でＢ駅に向う、の何れかを選択することになる。
【００１９】
本実施形態において、この経路選択モデルは、公知の非集計ロジットモデルで構築される。非集計ロジットモデルの効用関数及び選択確率は、それぞれ、式（１），（２）で与えられる。
【数１】
0004747144.tif 000002
【数２】
0004747144.tif 000003
【００２０】
式（１）において、「Ｕｍ」は経路ｍの効用、「Ｘｋｍ」は説明変数、「βｋ」はパラメータ、「α」は定数項である。また、式（２）において、「Ｕ１，Ｕ２」は、それぞれ経路１，２の効用、「Ｐｎ」は、経路ｎ（ｎ＝１，２）の選択確率である。
【００２１】
ところで、正常ダイヤ時の経路選択モデルでは、二つの経路の効用が全く等しい場合、選択確率Ｐはともに５０％となるように定められる。しかし、ダイヤ乱れ時には、二つの経路の効用が等しいならば、旅客は、当初予定していた「移動経路」を選択する、すなわち「待ち」を選択する確率が高いと考えられる。そこで、本実施形態では、この概念を反映させるために、「迂回」の効用Ｕａと「待ち」の効用Ｕｕとの差に閾値を設定し、双方の効用差がこの閾値を超えたときに「迂回」の選択確率Ｐａが５０％を超えるようなモデル構造としている。
【００２２】
従って、本実施形態の経路選択モデルにおいて、「迂回」の効用Ｕａは式（３）で与えられ、「待ち」の効用Ｕｕは式（４）で与えられる。そして、「迂回」の選択確率Ｐａは式（５）で与えられる。
【数３】
0004747144.tif 000004
【数４】
0004747144.tif 000005
【数５】
0004747144.tif 000006
【００２３】
式（３）によれば、「迂回」の効用関数Ｕａには、説明変数Ｘとして、全所要時間Ｘ１ａと、他社線ダミーＸ２とが含まれている。全所要時間Ｘ１ａは、迂回経路を経由した場合の対象駅から目的駅までの所要時間である。他社線ダミーＸ２は、迂回経路に他の鉄道会社（他社）の路線が含まれているか否かであり、含まれているならば「１」、含まれていないならば「０」となる。また、式（４）によれば、「待ち」の効用関数Ｕｕには、説明変数Ｘとして、全所要時間Ｘ１ｕと、待ち時間の逆数Ｘ３とが含まれている。全所要時間Ｘ１ｕは、通常経路を選択した場合の所要時間であり、運転再開までの待ち時間と、通常経路を経由した場合の対象駅から目的駅までの所要時間との和である。待ち時間の逆数Ｘ３は、上述の効用差の閾値に相当する。
【００２４】
図２に、待ち時間と効用差の閾値との関係の一例を示す。同図では、横軸を運転再開までの待ち時間、縦軸を効用差の閾値としたグラフを示している。同図によれば、待ち時間が短いほど閾値が大きくなっているとともに、待ち時間がある値（同図では、１０分程度）以下となると、閾値が急激に大きくなっている。これは、「待ち時間が短いほど、迂回をしない」傾向を表現しているとともに、待ち時間が短い（例えば、１０分程度）ならば、「迂回」を選択する確率が急激に低下することを表現している。
【００２５】
［構成］
図３は、本実施形態のシミュレーション装置１の内部構成を示すブロック図である。このシミュレーション装置１は、上述の経路選択モデルを用いて旅客それぞれの選択行動をシミュレーションする装置であり、ＣＰＵ１０と、入力部２０と、表示部３０と、通信部４０と、記憶部５０とが、バスＢにより互いにデータ通信可能に接続されて構成された公知のコンピュータシステムを用いて実現される。
【００２６】
ＣＰＵ１０は、記憶部５０に記憶されるプログラムやデータ、入力部２０から入力されたデータ等に基づいて、シミュレーション装置１を構成する各部への指示やデータの転送を行い、シミュレーション装置１の全体制御等を行う。また、本実施形態では、ＣＰＵ１０は、シミュレーションプログラム５１０に従ったシミュレーション処理を行う。このシミュレーション処理では、対象駅から目的駅に向う旅客それぞれについて、この対象駅から目的駅までの経路における列車運行を中断とした場合の選択行動を、上述の経路選択モデルをもとに決定する。
【００２７】
具体的には、先ず、シミュレーション条件を、例えばユーザによる操作入力に従って設定する。このシミュレーション条件には、シミュレーションの対象とする駅（対象駅）及び目的駅や、列車運行を中断する区間、運転中断時刻及び運転再開時刻（運転再開予定時刻）、運行中断の理由のほか、運転再開時刻や迂回経路の情報を対象駅に滞留している旅客（列車待ちの旅客）に提供（報知）する時刻等が含まれている。運行中断区間は、対象駅を含む区間であれば何れでも良いが、本実施形態では、対象駅と目的駅との間の区間とする。
【００２８】
設定されたシミュレーション条件は、条件テーブル５３１に格納される。図４に、条件テーブル５３１のデータ構成の一例を示す。同図によれば、条件テーブル５３１は、対象駅５３１ａと、目的駅５３１ｂと、運行中断区間５３１ｃと、運転中断時刻５３１ｄと、運転再開時刻５３１ｅと、運転中断理由５３１ｆと、迂回経路５３１ｇと、運転再開時刻の提供時刻５３１ｈと、迂回経路の提供時刻５３１ｉとを格納している。
【００２９】
運行中断区間５３１ｃは、対象駅と目的駅との間の区間である。迂回経路５３１ｇは、運行中断区間を迂回して対象駅から目的駅に向う経路である。運転再開時刻の提供時刻５３１ｈは、対象駅に滞留している旅客に対して運転再開時刻の情報を提供する時刻であり、運転中断時刻から運転再開時刻の間に設定される。迂回経路の提供時刻５３１ｉは、対象駅に滞留している旅客に対して迂回経路の情報を提供する時刻であり、運転中断時刻から運転再開時刻までの間に設定される。なお、この運転再開時刻及び迂回経路の提供は、対象駅における駅構内のアナウンス等を想定している。
【００３０】
次いで、ＣＰＵ１０は、設定したシミュレーション条件のもと、対象駅から目的駅に向う旅客それぞれの選択行動を決定する。すなわち、運転中断時刻から運転再開時刻までの間、旅客出現確率テーブル５２３に従って、対象駅に旅客を出現させる。
【００３１】
旅客出現確率テーブル５２３は、旅客の出現確率を定義したデータテーブルである。図５に、旅客出現確率テーブル５２３のデータ構成の一例を示す。同図によれば、旅客出現確率テーブル５２３は、出現駅と目的駅との組合せ毎に生成され、該当する出現駅５２３ａ及び目的駅５２３ｂとともに、時間帯５２３ｃ毎に旅客の出現確率５２３ｄを対応付けて格納している。出現確率５２３ｄは、目的駅に向う旅客の単位時間（例えば、１分）あたりの出現人数である。この出現確率は、例えば各駅に設置されている自動改札機等で収集された旅客の駅構内への入退場時刻のデータ等から決定されたり、各旅客の利用列車や時間帯毎の各駅間の乗車人数といった調査データをもとに決定される値である。
【００３２】
ＣＰＵ１０は、シミュレーション条件として定められた出現駅及び目的駅の組合せに該当する旅客出現確率テーブル５２３において定義された現在時刻に該当する時間帯の出現確率に従って、対象駅に旅客を出現させる。そして、出現させた旅客について、移動経路、迂回経路の認知有無及び運転再開想定時刻を設定する。移動経路は、出現駅から目的駅に向う経路であって、本実施形態では運行中断区間そのものとなる。迂回経路の認知有無は、例えば所定の確率に従って決定する。運転再開想定時刻は、当該旅客が予測する再開時刻であり、例えば所定の確率分布に従って決定する。
【００３３】
出現させた旅客についてのデータは、旅客データ５３２に格納される。図６に、旅客データ５３２のデータ構成の一例を示す。同図によれば、旅客データ５３２は、旅客毎に生成され、識別番号である旅客ＩＤ５３２ａと、出現駅５３２ｂと、出現時刻５３２ｃと、目的駅５３２ｄと、移動経路５３２ｅと、迂回経路認知有無５３２ｆと、運転再開想定時刻５３２ｇと、選択経路５３２ｈと、経路選択フラグ５３２ｉと、経路確定フラグ５３２ｊとを格納している。
【００３４】
迂回経路認知有無５３２ｆは、該当する旅客が迂回経路を認知しているか否かを示す。この迂回経路認知有無は、出現時には、例えば所定確率に従って「有り」或いは「無し」に設定され、迂回経路の提供がなされると「有り」に更新される。運転再開想定時刻５３２ｇは、該当する旅客が想定している運転再開時刻である。この運転再開想定時刻は、出現時には、該当する旅客自身が予測した時刻が設定されており、運転再開時刻の提供がなされると、この提供された時刻に変更される。なお、旅客自身が予測する再開時刻は、運転中断の理由に応じて異なることにしても良い。選択経路５３２ｈは、該当する旅客が選択した経路である。この選択経路は、出現時には「待ち」に設定され、経路選択を行うと、選択された経路（すなわち、「待ち」或いは「迂回」）に更新される。経路選択フラグ５３２ｉは、経路選択を行ったか否かを示すフラグである。この経路選択フラグは、出現時には選択を行っていないことを示す「０」に設定され、１回目の経路選択を行うと「１」に更新される。経路確定フラグ５３２ｊは、選択した経路を確定したか否かを示すフラグである。この経路確定フラグは、出現時には確定していないことを示す「０」に設定され、経路選択を行って選択経路が確定されると「１」に更新される。
【００３５】
続いて、ＣＰＵ１０は、出現させた旅客それぞれについて、上述の経路選択モデルにもとづく経路選択を行う。図７は、旅客についての経路選択の概要を示す図である。経路選択のタイミングは、該当する旅客が迂回経路を認知しているか否かによって異なる。
【００３６】
同図（ａ）は、迂回経路の認知有りの旅客について示している。同図（ａ）に示すように、迂回経路の認知有りの旅客については、先ず、出現させたタイミング（時刻ｔｘ）で１回目の経路選択を行う。その結果、「迂回」を選択した場合には、経路選択を「確定」し、該旅客を対象駅から迂回経路の列車に乗車させる。一方、「待ち」を選択した場合には、該旅客を対象駅に滞留させる。その後、運転再開時刻の提供がなされたタイミング（時刻ｔｙ）で２回目の経路選択を行う。その結果、「待ち」を選択したならば、該旅客を対象駅に滞留させたままとし、「迂回」を選択したならば、該旅客を対象駅から迂回経路の列車に乗車させて対象駅の滞留旅客から除く。そして、該旅客の経路選択を「確定」する。
【００３７】
同図（ｂ）は、迂回経路の認知無しの旅客について示している。同図（ｂ）に示すように、迂回経路の認知無しの旅客については、出現したタイミング（時刻ｔｘ）では経路選択を行わず、迂回経路の提供がなされたタイミング（時刻ｔｚ）で１回目の経路選択を行う。これは、旅客が迂回経路を認知していない場合には、必然的に選択肢が「待ち」のみとなるためである。その結果、「迂回」を選択した場合には、経路選択を「確定」し、対象駅から迂回経路の列車に乗車させて滞留旅客から除く。一方、「待ち」を選択した場合には、該旅客を対象駅に滞留させる。その後、再開時刻の提供がなされたタイミング（時刻ｔｙ）で２回目の経路選択を行う、その結果、「待ち」を選択したならば、該旅客を対象駅に滞留させたままとし、「迂回」を選択したならば、該旅客を対象駅から迂回経路の列車に乗車させて対象駅の滞留旅客から除く。そして、該旅客の選択経路を「確定」する。
【００３８】
旅客の経路選択は、上述の経路選択モデルに従って行う。すなわち、式（３）に従って、「迂回」に対する効用Ｕａを算出する。このとき、全所要時間Ｘ１ａは、経路選択を行う時点の時刻（対象時刻）ｔを基準とした、迂回経路を経由した場合の対象駅から目的駅までの所要時間とする。他社線ダミーＸ２は、迂回経路に他社線の路線が含まれているならば「１」とし、含まれていないならば「０」とする。
【００３９】
また、式（４）に従って、「待ち」に対する効用Ｕｕを算出する。このとき、全所要時間Ｘ１ｕは、上記対象時刻における運転再開までの待ち時間Δｔ１と、通常経路を経由した場合の対象駅から目的駅までの所要時間Δｔ２との和とする。待ち時間Δｔ１は、対象時刻ｔから該旅客の再開想定時刻までの時間である。この再開想定時刻は、運転再開時刻の提供がなされていないならば、該旅客が予測するとして定められた時刻であり、運転再開時刻の提供がなされているならば、この運転再開時刻である。そして、所要時間Δｔ２は、上記再開想定時刻を基準とした、通常経路を経由した場合の通常ダイヤに基づく対象駅から目的駅までの所要時間である。ここで、通常ダイヤは通常ダイヤデータ５２４に格納されている。この通常ダイヤデータ５２４には、自社の各路線のダイヤデータのほか、迂回経路に含まれている他社の各路線のダイヤデータも含まれている。
【００４０】
次いで、算出した効用Ｕｕ，Ｕａをもとに、式（５）に従って、「迂回」の選択確率Ｐａを算出する。そして、選択確率Ｐａをもとに、経路を選択する。
【００４１】
ここで、式（３）〜（５）の各関数式は、経路選択モデル式データ５２１に格納されている。また、各効用関数におけるパラメータβの値は、パラメータテーブル５２２に格納されている。図８に、パラメータテーブル５２２のデータ構成の一例を示す。これらのパラメータの値は、所定のアンケート調査結果をもとに推計されたものである。
【００４２】
図３に戻り、入力部２０は、例えばキーボードやマウス、タッチパネル及び各種スイッチ等によって実現される入力装置であり、操作入力に応じた入力信号をＣＰＵ１０に出力する。表示部３０は、例えばＬＣＤ（Liquid Crystal Display）やＥＬＤ（Electronic Luminescent Display）等によって実現される表示装置であり、ＣＰＵ１０から入力される表示信号に基づく各種画面を表示する。通信部４０は、例えば、無線通信モジュールやルータ、モデム、ＴＡ、有線用の通信ケーブルのジャックや制御回路等によって実現される通信装置であり、外部機器との間でデータ通信を行う。
【００４３】
記憶部５０は、ＣＰＵ１０がシミュレーション装置１を統合的に制御するための諸機能を実現するためのシステムプログラムや、シミュレーションを実現するためのプログラムやデータ等を記憶しているとともに、ＣＰＵ１０の作業領域として用いられ、ＣＰＵ１０が各種プログラムに従って実行した演算結果や、入力部２０からの入力信号等が一時的に格納される。本実施形態では、記憶部５０は、プログラムとしてシミュレーションプログラム５１０を記憶しているとともに、データとして、経路選択モデル式データ５２１と、パラメータテーブル５２２と、旅客出現確率テーブル５２３と、通常ダイヤデータ５２４と、条件テーブル５３１と、旅客データ５３２とを記憶している。
【００４４】
［処理の流れ］
図９は、ＣＰＵ１０が実行するシミュレーション処理の流れを説明するためのフローチャートである。同図によれば、ＣＰＵ１０は、先ず、ユーザからの操作入力等に従ってシミュレーションの条件を決定する（ステップＡ１）。その後、決定したシミュレーション条件のもと、運転中断時刻ｔａから運転再開時刻ｔｂまでの間のシミュレーションを開始する。
【００４５】
すなわち、運転中断時刻ｔａを、シミュレーション時刻である仮想時刻ｔに設定する（ステップＡ３）。続いて、旅客出現処理を行って、対象駅に旅客を出現させる（ステップＡ５）。その後、対象駅に滞留している旅客それぞれを対象としたループＡの処理を行う。ループＡでは、先ず、対象旅客の経路確定フラグＦ２を判断し、「０」ならば（ステップＡ７：ＮＯ）、更に、対象旅客の迂回経路の認知有無を判断する。仮想旅客の迂回経路の認知有無が“無し”ならば（ステップＡ９：ＮＯ）、仮想時刻ｔが、迂回経路の提供時刻に達しているか否かを判断し、達しているならば（ステップＡ１１：ＹＥＳ）、対象旅客の迂回経路の認知有無を“有り”に設定する（ステップＡ１３）。
【００４６】
続いて、仮想時刻ｔが運転再開時刻の提供時刻に達したかを判断し、達しているならば（ステップＡ１５：ＹＥＳ）、対象旅客の再開想定時刻を運転再開時刻に設定する（ステップＡ１７）。そして、経路選択処理を行って、対象旅客の経路選択を行う（ステップＡ２１）。
【００４７】
一方、仮想時刻ｔが運転再開時刻の提供時刻に達していないならば（ステップＡ１５：ＮＯ）、対象旅客の経路選択フラグＦ１を判断し、「０」ならば（ステップＡ１９：ＹＥＳ）、経路選択処理を行って、対象旅客の経路選択を行う（ステップＡ２１）。
【００４８】
経路選択処理を行うと、ＣＰＵ１０は、該処理で選択した経路が「迂回」ならば（ステップＡ２３：ＹＥＳ）、対象旅客の経路確定フラグＦ２を「１」に設定し、選択経路を「迂回」として確定する（ステップＡ２５）。一方、選択した経路が「待ち」ならば（ステップＡ２３：ＮＯ）、続いて、対象旅客の経路選択フラグＦ１を判定する。経路選択フラグＦ１が「１」ならば（ステップＡ２７：ＹＥＳ）、今回の経路選択が２回目の経路選択であると判断し、経路確定フラグＦ２を「１」に設定して、選択経路を「待ち」として確定する（ステップＡ２９）。一方、経路選択フラグＦ１が「０」ならば（ステップＡ２７：ＮＯ）、今回の経路選択が１回目の経路選択であると判断し、経路選択フラグＦ１を「１」に設定する（ステップＡ３１）。ループＡはこのように行われる。
【００４９】
全ての旅客を対象としたループＡの処理を終了すると、ＣＰＵ１０は、仮想時刻ｔを所定の時分（例えば、１分）だけ進める（ステップＡ３３）。その結果、仮想時刻ｔが運転再開時刻に達していないならば（ステップＡ３５：ＮＯ）、ステップＡ５に戻る。一方、仮想時刻ｔが運転再開時刻に達したならば（ステップＡ３５：ＹＥＳ）、シミュレーションを終了し、シミュレーション結果を出力する（ステップＡ３７）。例えば、旅客データ５３２を参照し、各旅客の選択経路（「待ち」或いは「迂回」）を一覧表示させたり、選択経路が「待ち」である旅客の人数、すなわち運転再開時刻における対象駅の滞留旅客人数を算出して表示させる。以上の処理を行うと、ＣＰＵ１０は、本シミュレーション処理を終了する。
【００５０】
［作用・効果］
このように、本実施形態によれば、対象駅に出現した旅客それぞれについて、ロジットモデルである経路選択モデルに従って、そのまま運転再開を待つ「待ち」、及び運行中断区間を迂回する迂回経路を選択する「迂回」の何れかが選択される。経路選択モデルでは、「迂回」に対する効用Ｕａ及び「待ち」に対する効用Ｕｕをもとに迂回の選択確率Ｐａを算出し、算出した選択確率Ｐａに従って「迂回」及び「待ち」のいずれかを選択する。「待ち」の効用Ｕｕにおける説明変数Ｘである全所要時間Ｘ１ｕには、運転再開までの待ち時間が含まれるが、旅客毎に想定する運転再開時刻（運転再開想定時刻）が異なるため、「迂回」の選択確率Ｐａは旅客毎に異なる。このように、仮想旅客それぞれについて「迂回」の選択確率Ｐａを算出して経路を選択することで、仮想旅客１人１人の行動を忠実に再現した精度の良いシミュレーションが可能となり、その結果、運行中断といったダイヤ乱れ時の旅客需要をより正確に予測することが可能となる。
【００５１】
なお、本発明の適用可能な実施形態は上述の実施形態に限定されることなく、本発明の趣旨を逸脱しない範囲で適宜変更可能なのは勿論である。
【００５２】
（Ａ）サーバ・クライアントシステム
例えば、上述の実施形態では、本発明を、１台のシミュレーション装置に適用することにしたが、これを、サーバ装置及びクライアント装置から構成されるサーバ／クライアント型のシミュレーションシステムで適用することにしても良い。この場合、ユーザは、クライアント装置においてシミュレーション条件を入力し、入力されたシミュレーション条件は、クライアント装置からサーバ装置へ送信される。すると、サーバ装置は、クライアント装置から送信されてきたシミュレーション条件をもとづいてシミュレーションを行い、シミュレーション結果として、決定した各旅客の選択経路（「待ち」或いは「迂回」）の一覧や、対象駅の滞留旅客人数等を含むデータを、クライアント装置に送信する。そして、クライアント装置では、サーバ装置から送信されてきたシミュレーション結果を、例えば表示等して出力する。
【図面の簡単な説明】
【００５３】
【図１】経路選択モデルの概要図。
【図２】待ち時間と効用差の閾値との関係例。
【図３】シミュレーション装置の内部構成図。
【図４】条件テーブルのデータ構成例。
【図５】旅客出現確率テーブルのデータ構成例。
【図６】旅客データのデータ構成例。
【図７】経路選択のタイミングの概要図。
【図８】パラメータテーブルのデータ構成例。
【図９】シミュレーション処理のフローチャート。
【符号の説明】
【００５４】
１  シミュレーション装置
１０  ＣＰＵ、２０  入力部、３０  表示部、４０  通信部
５０  記憶部
５１０  シミュレーションプログラム
５２１  経路選択モデル式データ、５２２  パラメータテーブル
５２３  旅客出現確率テーブル、５２４  通常ダイヤデータ
５３１  条件テーブル、５３２  旅客データ
【図１】
0004747144.tif 000007
【図２】
0004747144.tif 000008
【図３】
0004747144.tif 000009
【図４】
0004747144.tif 000010
【図５】
0004747144.tif 000011
【図６】
0004747144.tif 000012
【図７】
0004747144.tif 000013
【図８】
0004747144.tif 000014
【図９】
0004747144.tif 000015