• 検索結果がありません。

About the settlement of the parameters in urban expressway simulation model *

N/A
N/A
Protected

Academic year: 2022

シェア "About the settlement of the parameters in urban expressway simulation model * "

Copied!
4
0
0

読み込み中.... (全文を見る)

全文

(1)

1

都市内高速道路シミュレーションモデルにおけるパラメータの設定について*

About the settlement of the parameters in urban expressway simulation model *

森内正寿**・森田綽之***・吉井稔雄****・小根山裕之*****・島崎雅博******

By Masatoshi MORIUCHI**・Hirohisa MORITA***・Toshio YOSHII****・

Hiroyuki ONEYAMA*****・Masahiro SHIMAZAKI******

1.はじめに

首都高速道路公団では、交通渋滞の延伸の表現や 時々刻々と変化する交通状況の評価が可能で、かつ 柔 軟 性 を あ わ せ も つ 「 動 的 シ ミ ュ レ ー シ ョ ン モ デ ル 」TRANDMEX(TRANsportation Dynamic Model on urban EXpressway)を開発してきた。これは、シ ミュレ ーションモデ ル SOUND1)を基本として、首 都高速道路公団の実務へ適用するために、データ作 成労力軽減等の改良を加えたモデルである。

ここで、TRANDMEX において、データの簡略化 を図り、高い再現性を得るために採用したロジック の考 え方、及びパラメー タ(Q-K 曲線、交通容 量 など)の設定方法について述べることとする。

2.TRANDMEXの特徴   TRANDMEXのシ ミュレーション手順 は、図‑1に示すよう にオンランプから車 両が発生し、発生し た車両が移動と経路 選択を繰り返しなが ら目的地であるオフ ランプを目指して 

 

走行するロジックを適用している。 

3.TRANDMEXのロジック及び設定パラメータ

(1)車両の発生

演 算 時 間 を 短 縮 す る た め に 、 同 じ 属 性 の 車 両 数 台を1つの車群として取り扱う。同じ属性は①ラン

プ間 OD、②時間帯、③経路選択特性が同じ場合に

同一として扱う。なお、現在は、シミュレーション 単位時間を3秒、車群の単位を3台とした。これは、

精度的には1 秒1台が優れているが、実務への適用 性を考慮して、演算時間が大幅に短縮できることか ら採用した。

(2)車両の移動 

(a)Q‑K 曲線の適用 

車両速度は、自車(B)と前車(A)の車頭距離 S とそ の走行リンク上の Q-K 曲線から算出する。その速 度にシミュレーション単位時間を乗じて移動距離を 算出し、車両を移動させる。(図-2)

な お 、 車 両 は 一 本 の レ ー ル の 上 を 走 行 す る と 考 える。このため、車線変更については考慮しない。

ただし、織り込み現象による交通容量の低下につい ては、入力条件である Q-K 曲線の交通容量で考慮 することで対応している。

図‑2  車両移動の概念 

*キーワーズ:配分交通、経路選択、ネットワーク交通流

**正員,工修,首都高速道路公団神奈川建設局建設第二部

***フェロー員,工博,日本大学総合科学研究所教授

****正員,工博,高知工科大学社会システム工学科助教授

*****正員,工博,東京大学生産技術研究所

******正員,工修,パシフィックコンサルタンツ(株)

交通技術本部  道路部

    (東京都新宿区西新宿2-7-1,TEL:03-3344-0074,

E-mail:[email protected]

B

Bの移動前の車頭距離 L Bの移動距離 D Bの移動後の車頭距離 S

B B

A A A

車両Aの移動

時刻

時刻 t+⊿t

1スキャンタイム

  V=Q/K=f(K)/K=f(1/S)/(1/S)

  L=D+S、D=V×(Δt/3600)

車両の発生(オンランプ)

首都高上の車両の移動

JCT部での経路選択

目的オフランプ へ到着したか

目的オフランプへの到着 NO

YES

図‑1  シミュレーション手順 

(2)

2

(b)Q‑K 曲線の設定

  各リンクの Q-K 曲線については、首都高速上の 車両感知器の5分間データをシステマティックにア ップデートして得られたものを用いており、様々な 交通状況を包括した平均的な値と位置づけられる。

この Q-K 曲線は、図-3 に示すように自由流側・

渋滞流側いずれも2次曲線で近似されているが、演 算時間の短縮を考慮して交通挙動の影響の少ない自 由流側を直線式に置き換えた。

図‑3  Q‑K 曲線(リンクごとに設定) 

(c)大型車混入率の考慮

  図-3 の Q-K 曲線は実台数により構成されている が、交通現況を再現する上で大型車の影響を無視す ることはできない。この影響をシミュレーションで 考慮するため、区間を通過する日平均大型車混入率 で Q-K 曲線を相似拡大(図-4 の(1)の補正)し、時 間帯別にOD交通量を補正した。ここでは、大型車 換算係数は2.0を用いた。

(d)渋滞時の交通容量の低減

交通容量が発現する付近の Q-K が不連続かつ不 安定な状態であり、自由流域における最大通過交通 量が渋滞時におけるものより大きいこと、一度渋滞 現象が生じた場合に回復しがたいことを表現するた めに、ボトルネック箇所において、渋滞状態が 20 分以上継続した場合、交通容量を5%減少させ、さ らに渋滞状況が 20 分継続した場合はさらに5%減 少させるロジックを導入した。(図-4の(2)の補正)

図‑4  Q‑K 曲線の補正 

(2)経路選択

(a)2経路選択モデル

  経路選択モデルは、一般的に用いられる方法とし てロジットモデルをベースとしたDial配分2)があ る。しかし、この方法ではIIA(重複した経路を独 立した経路として経路選択を行う問題)とエフィシ ェントパス(経路が複数存在しても、その内の一部 が経路選択対象から外れる現象)の問題がある。こ れらの問題に対応するために経路選択用のネットワ ークを別途準備する方法などもある(SOUND では、

走行用ネットワークとは別に、経路選択用のネット ワークでロジットモデルを用いて経路比率を算出し ている)が、作業効率上労力が大きい。そこで、静 的配分のネットワークをそのまま用い、作業労力軽 減を図ることができる2経路選択モデルを採用した。

  このモデルは、JCT 等の分流部手前で、目的オフ ランプまでの右左折した場合のそれぞれの最短経路 時間から分流比率を算出している。

  たとえば、浜崎橋 JCT から目的オフランプへ向 かう経路が複数存在する場合、そのうち左折した場 合 の 最 短 経 路 ( ① ) 、 右 折 し た 場 合 の 最 短 経 路

( ② ) に 着 目 し て 、 そ の 2 経 路 の 旅 行 時 間 を 図-5 の式へ適用することにより分流比率を算出している。

図‑5  経路選択方法   

 

江戸橋JCT 5号池袋線

4号新宿線

3号渋谷線

2号目黒線

1号羽田線 浜崎橋JCT

至目的ランプ

至出発ランプ

 

交通量

交通密度 y=  xar 2+bx+c

(1)

(2)

(1):大型車混入率による拡大 (2):渋滞による低減

Pi:経路iの選択確率

n:経路数(2経路)

θ:トレードオフパラメータ

ti:経路iの旅行時間

tj:経路jの旅行時間

e

-θti

Pi= 

Σe

-θtj

j=1

区間ごとに設定 交通量

交通密 y=ax2+bx+c

y=a'x

(3)

3

(b)トレードオフパラメータ

経 路 選 択 確 率 の ト レ ー ド オ フ パ ラ メ ー タ θ に つ い て は 、 第 22 回 首 都 高 速 道 路 交 通 起終 点 調 査 3)

(以降 22OD と記す)における代表的な JCT 間の 時間帯別経路選択比率と車両感知器データから得ら れる時間帯別経路別旅行時間に基づいて経路選択確 率の残差平方和が最小となる 0.00835(1/秒)を用 いた。

こ こ で 、 θ の 感 度 分 析 で は 、 大 き く 設 定 す る 場 合はあまり影響はなかった。しかし、θを小さくす ると、都心環状線等で周回交通が発生し、現実の交 通状況が再現されない場合があった。

ま た 、 ト レ ード オ フ パ ラ メ ータ θ は 運転 者 が 認 知する情報の正確さを表すものであり、交通状況情 報の提供が高度化する将来において、この値を小さ く設定することは考えられないため、適用にあたっ ては問題ないと判断した。

(c)経路選択時間間隔

経 路 選 択 時 間間 隔 は 、 首 都 高速 の 道 路情 報 提 供 更新間隔である1分を採用した。この間隔を長く設 定すると、交通が集中し渋滞が過大に発生し現況再 現性が下がる傾向にあった。

(d)交通状況依存層と非依存層の設定

経 路 選 択 行 動 に つ い て 、 所 要時 間 の 情報 提 供 に 対 す る 反 応 の 違 い 等 に よ っ て 、 ① 交 通 状 況 依 存 層

(提供情報を受けて経路選択を行う)、②交通状況 非依存層(自由流時で選択確率に従って経路選択を 行う)に分類した。

首都高速道路公団で実施されたアンケート調査 4)

結果(約 40%)、および吉井ら 1)の推計結果(40%

〜60%)を参考に50:50と設定した。

4.現況再現におけるパラメータの設定方法

(1)交通容量の管理

現 況 再 現 を 行 う に あ た り 、 パ ラ メ ー タ の 設 定 と し て 最 も 影 響 を 受 け る の が 、 ボ ト ル ネ ッ ク 箇 所 の Q-K曲線における交通容量の設定である。

こ こ で 、 車 両 感 知 器 か ら 近 似 式 で 定 義 さ れ る 交 通容量は、実際に車両感知器では観測され難い大き な値を示すことが多く、過大な交通容量を設定して しまうこととなる。そこで、ボトルネック部では、

実態データから速度が低下している時間帯における 平均通過交通量を算出し、ボトルネック部の交通容 量として適用した。この交通容量管理を行うために、

Q-K 曲線の上部をカットした図-6 に示す台形型の Q-K曲線を採用した。

ここ で、 図-7 は、ボ トルネ ッ ク部に お ける通 過 交通量と速度の関係を示しており、渋滞発生時より 渋滞発生中の通過交通量が 10%程度減少している ことを示している。この状況を Q-K プロット図に 示すと図-8 に示すような Q-K 曲線となる。(点線 は、渋滞中の Q-K 曲線で、交通容量が低下してい ることを現している)

図‑6  ボトルネック箇所の Q‑K 曲線 

図‑7  ボトルネック箇所の交通量と速度 

図‑8  ボトルネック箇所の Q‑K プロット図 

(2)交通容量の設定

  上記のように管理されたボトルネック箇所の交通 容量を以下の手順で設定した。

①  車両感知器で首都高速道路上のボトルネックと なっている箇所を抽出。(27箇所抽出)

②  抽 出 し た ボ ト ル ネ ッ ク 箇 所 の 最 大 通 過 交 通 量

①ボトルネック部(葛西→新木場)

0 100 200 300 400 500 600 700

0:00 1:00 2:00 3:00 4:00 5:00 6:00 7:00 8:00 9:00 10:00 11:00 12:00 13:00 14:00 15:00 16:00 17:00 18:00 19:00 20:00 21:00 22:00 23:00

時間帯 交通量(台/5分)

0 20 40 60 80 100 120 速度(km/h)

交通量 速度 6,431(7,589)台/時 5,799(6,843)台/時

(  )内の数値は、乗用車換算値

 

交通量Q

交通密度K 実態Q-Kにお

ける交通容量 実測における 通過交通量

得られる実測データ

Q-Kプロット(葛西→新木場)

0 100 200 300 400 500 600 700

0 50 100 150 200

密度(台/km)

交通量(/5分)

04:00~08:00 09:00~22:00

(4)

4

(1時間)を車両感知器データから算出。

③  各ボトルネック箇所の交通容量を最大通過交通 量に設定。

ここで、設定した交通容量でシミュレーションを 実施した場合でも、再現性が不十分であったため、

これらの交通容量の微調整を行うこととした。この とき、首都高速道路のネットワークが複雑であるた め、各ボトルネックが互いに干渉しあって、1 箇所 の交通容量をわずかに変化させただけで全線の交通 状況が大きく変化することがわかった。 

そ こ で 、 微 調 整 に お い て は 、 各 ボ ト ル ネ ッ ク 別 に感度を整理し(表‑1)、調整を図った。 

表‑1  各渋滞箇所の影響 

5.現況再現性の検証

  モデルの検証は、22OD のランプ間交通量データ を入力交通量として 24 時間連続したシミュレーシ ョンを行い、当日の車両感知器データと比較を行っ た。交通容量の微調整を図った結果、各時間帯の混 雑量(図-9)、渋滞箇所(図-10)は概ね一致し、1日の 交通状況を再現することができたと考える。よって、

TRANDMEX における基本的なロジック及びパラメ

ータの設定については確立できたと考える。

図−9  時間帯別混雑量の変化 6.おわりに

  実務を行うにあたっての今後の課題として、①時 間帯別将来ランプ間 OD表の作成方法、②将来路線 の Q-K 曲線の設定方法、③車線運用変化に対応し た Q-K 曲線の設定方法が挙げられ、これらについ

て今後検討していく必要があると考えている。 

図−10  現況再現結果(速度ランク図) 

参考文献 

1) 吉井稔雄、桑原雅夫、森田綽之;都市内高速道 路における過飽和ネットワークシミュレーショ ンモデルの開発,交通工学Vol30,No1,1995 2)  土木学会;交通ネットワークの均衡分析,第4

章,1997 

3)  首都高速道路公団;第 22 回首都高速道路交通 起終点調査報告書,平成7年9月実施 

4) 首都高速道路公団;首都高速道路の利用行動の 特性に関する調査,pp.139〜158,平成5年3月

①7時台

(シミュレーション結果)

(22OD当日の車両感知器データ)

②17時台

(シミュレーション結果)

(22OD当日の車両感知器データ)

混雑量

(首都高速道路全線)

0 20 40 60 80 100 120 140

3 4 5 6 7 8 9 1011 121314 1516 171819 20 21 22 23 0 1 2 時間帯

混雑量(km・時)

実測(22OD)

シミュレーション結果

芝公園ランプ合流後

3号渋谷線合流後

北の丸合流後

1号羽田線合流後

飯倉トンネル

4号新宿線合流後

5号池袋線合流後 -

用賀料金所

代々木ランプ合流後

志村料金所

9号深川線合流後

6号三郷線合流後

船堀橋ランプ合流後

中央環状線合流後

大井料金所

北上野料金所

羽田ランプ合流後

三ツ沢料金所

花之木ランプ合流後

11号台場線合流後

新木場ランプ分流手

羽田ランプ合流後

外苑ランプ合流後

飯田橋ランプ合流後

池尻ランプ合流後

※渋滞箇所の交通容量を低下させると渋滞が減る方向:「-」

※渋滞箇所の交通容量を低下させると渋滞が増える方向:「+」

渋滞箇所

6号向島線合流後

3号 狩場線

9号深川線合流後

中央 環状線湾岸線 1号

横羽線 2号 三ツ沢線 11号

台場線 1号 上野線

6号 三郷線川口線 5号

池袋線 6号 向島線

7号 小松川線

9号 深川線 1号

羽田線 2号 目黒線

3号 渋谷線

4号 新宿線

参照

関連したドキュメント

© Metropolitan Expressway Company Limited

(2015) presents an integrated model to explore the spatial distribution of outdoor recreation trips of urban residents based on their case study in Salford, UK. Using

Software In the Loop Simulation 機能検証 計測制御技術 計測制御技術 〈 エンジニアリングもご提供いたします 〉 Hardware In the Loop Simulation システム検証

Aid policies by the World Bank and others on urban lending in Asia dramatically changed in mid 1980s from individual poverty-oriented projects to the institutional development

Aerodynamic drag accounts for a large portion of the total resistive force when the solar car is running in the race. It is a significant power loss at any speed above

Using ETC Panel Data to Analyze Behavioral Changes in Response to a Toll Revision on an Urban Expressway. 後藤

Cruz Joven—the practices of piety (exercitia pietatis) introduced to Timor-Leste during the Indonesian era—was brought to the settlement to commemorate 500 years of the pres- ence

This sub-model, which employed the step-down formula and the overflow formula, was found to overestimate storm water interaction discharge when compared