natural science 学会　講演概要

領域１【物理から生物】
アリの行列とバスのダンゴ運転

友枝 明保†,1,2、西成 活裕2,3
1 明治大学先端数理科学インスティテュート・〒213-8571, 神奈川県川崎市多摩区東三田1-1-1
2 東京大学工学系研究科 航空宇宙工学専攻・〒113-8656, 東京都文京区本郷7-3-1
3 独立行政法人 科学技術振興機構 さきがけ・〒113-8656, 東京都文京区本郷7-3-1
(Tel : 81-3-5841-6631; Fax : 81-3-5841-8561)
(†Email address:tomoeda◎mail.ecc.u-tokyo.ac.jp)

概要: 確率セルオートマトンモデルを粒子と場の相互作用が存在する系に拡張し，アリの並進運動ダイナミクスとバスの運行システムに応用した．これらの二つの数理モデルの数値シミュレーションからいずれの場合も密度が上がるに連れて速度が単調減少するのではなく，速度が最大となる最適な密度が存在することがわかった．さらに，二つの数理モデルを渋滞現象という視点で捉えると数学的に完全に等価であることもわかった．

Keywords: 渋滞学，数理モデル，セルオートマトン，確率過程

I. 研究の背景

今日の社会では「渋滞」と呼ばれる現象がいたるところで見られる．例えば，高速道路での車の長い列，花火大会が終わった後の人の混雑，バスや電車のダンゴ運転などが挙げられる．一方，自然界においても，魚・鳥の群れやアリの行列など生物の集団運動としての渋滞現象も観測されている．これらの多岐に渡る一見全く関係の無さそうな渋滞現象を統一的に扱う最先端の学問が「渋滞学」である[1]．このような様々な渋滞現象に共通することは，すべて排除体積効果を持った自分自身で能動的に動くことのできる自己駆動粒子の多体相互作用で渋滞を捉えることができる，ということである．自己駆動粒子はその粒子が持っている知性や心理のため，個々の振る舞いが物理学の大原則である作用=反作用の法則や慣性の法則に従わず，微分方程式による運動の記述が困難となってしまう．そこで，このような現象を記述する際に有効となるのが，個々の粒子の振る舞いをルールによって記述するセルオートマトンモデルである．セルオートマトンモデルの利点は，そのルールによる記述のみならず，時間・空間に関する全ての物理量が離散化されており， コンピューターシミュレーションが容易であることや，特にAsymmetric Simpler Exclusion Process (ASEP) [2] や Zero Range Process (ZRP) [3] と呼ばれる定常状態における粒子の分布を厳密に求めることができる可解確率セルオートマトンモデルが存在することである．

本稿では，確率セルオートマトンモデルを粒子と場の相互作用系に拡張し，アリの並進運動とバスの運行システムに応用した数理モデルについて紹介する．セルオートマトンモデルを用いる準備として，空間をLセルにわけ，各セルにラベルi(i=1,2,･･･L)で番号付けをする．さらに粒子は一次元の道を一方向にだけ進むものとし，時刻tの粒子の変数をUti，場の変数をσtiとする．

II. 確率セルオートマトンモデル

1. アリの並進運動モデル[4]

アリ同士はフェロモンを用いてお互いのコミュニケーションを実現しているため，アリ粒子とフェロモン場の相互作用系を考え，以下のようなダイナミクスでモデル化する．

(ステージ１)
アリがある時刻tでセルiにいたとする(Uti=1)．もし前にアリがいる(Uti+1=1)ならば，セルiのアリは動かない(体積排除効果)．そして前にアリがいないときは，前に進もうとするがフェロモン(σti)のある/なしに応じて進むホップ確率pが以下のように変わる．


このように2種類の確率q,Qを導入し，フェロモンがある方が動きやすいので，一般的にq<Qとおける．

(ステージ２)
ステージ１においてアリがいるセルにはすべてフェロモンを生成する．つまりもし，U^{t+1}_i=1ならば，σ^{t+1}_i=1とする． また，アリがいないセルのフェロモンは確率fで蒸発するとする．すなわち，U^{t+1}_i=0かつσ^{t}_i=1ならば，確率fでσ^{t+1}_i=0とする．

2. バスの運行システムモデル[5]

バスの運行システムでは，バスとバス停で待つ乗客の相互作用系を考え，以下のようなダイナミクスでモデル化する．

(ステージ１)
アリと同様，バスがある時刻tでセルiにいた場合(U^{t}_i=1)，もし前にバスがいる(U^{t}_{i+1}=1)ならば，セルiのバスは動けず，その一方， 前にバスがいないときは，前に進もうとするが待つ乗客のいる/いない(σ^{t}_i)に応じて進むホップ確率pが(1)式のように変わると考えることができる． 待つ乗客がいれば乗り降りに時間がかかるためq<Qが従う．

(ステージ２)
ステージ１においてバスがいるセルでは乗客をバスに乗せ，待っている乗客はいなくなる(U^{t+1}_i=1ならば，σ^{t+1}_i=0)． 一方，バスのいないセルには確率fで乗客が到着するとする(Ut+1i=0かつσti=0ならば，確率fでσt+1i=1)．

III. 各数理モデルのシミュレーション結果

上述した二つのモデルそれぞれに対して数値シミュレーションを行った．境界は周期境界条件とした．アリの並進運動モデルの結果が図１であり，バスの運行システムモデルの結果が図２である．各図は粒子の密度に対して平均速度をシミュレーションしたものである．

図１：アリの並進運動モデル(Q=0.75, q=0.25)
図１：アリの並進運動モデル(Q=0.75, q=0.25)

図２：バスの運行システムモデル(Q=0.9, q=0.5)
図２：バスの運行システムモデル(Q=0.9, q=0.5)

図１，２それぞれの速度密度関係図から，どちらの場合も，密度が増加するに従って速度が単調に減少するわけではなく，ある密度で速度の最大値をとる確率fが存在することがわかる．つまり，場の物理量の生成/消滅の確率によっては，粒子の密度がある程度増えても渋滞せず，逆に速度を上げて粒子が動くことができる状況が起こりえるのである．

IV. 結論

本稿では確率セルオートマトンモデルを用いた二つの応用数理モデルを紹介した．これらの二つの数理モデルのシミュレーション結果からいずれの場合も，密度がある程度増えても速度が減少せず，逆に速度が上がり，速度を最大にする最適な密度が存在するという結果が得られた．

これらの二つのダイナミクスは異なる自然/社会現象を記述したものであるが，渋滞現象という見地に立ち，フェロモンの存在のある/なしを待っている乗客のいない/いると読み替えることで，数学的に完全に等価なモデルとなっていることもわかった．数理モデルを用いて様々な分野の現象を統一的に取り扱うことで，バスのダンゴ運転解消のアイデアがアリの行列行進ダイナミクスのアイデアから生まれる可能性も大いにあり，数理的にも社会的にも極めて興味深い結果が得られた．

参考文献

[1] 渋滞学，西成活裕，新潮選書
[2] B. Derrida, et al., J. Phys. A. Math. Gen., 26 (1993), 1493.
[3] F. Spitzer, Adv. In. Math. 5 (1970), 246.
[4] K. Nishinari, et al., Phys. Rev. E, 67 (2003), 036120.
[5] A. Tomoeda, et al., Physica A, 384 (2007), 600.

natural science 学会　目次

領域 領域代表 領域アドバイザー
領域１
【物理から生物】		 沢田康次
（東北工業大学学長）		 武田曉（東京大名誉教授、東北大名誉教授）
領域２
【物理から環境】		 今村文彦
（東北大学大学院工学研究科 災害制御研究センター）		 湯本良次（東北学院榴ヶ岡高校 副校長）
領域３
【生物から環境】		 本川達雄
（東京工業大学 生命理工学部）		 菅原研（東北学院大学 准教授）、千葉聡（東北大学大学院生命科学研究科 准教授）

発表題目

領域１ アリの行列とバスのダンゴ運転 ○友枝明保1,2、西成活裕2,3 1．明治大学先端数理科学インスティテュート、2．東京大学工学系研究科、3．独立行政法人科学技術振興機構さきがけ
ボトルネックにおける人の渋滞軽減方法 ○栁澤大地1,2、西成活裕1,3 1．東京大学大学院工学系研究科、2．(独)日本学術振興会特別研究員DC1、3．(独)科学技術振興機構さきがけ
視覚情報の遮断によるヒトの視覚運動機構の変化 ○田村友里恵1,3、林叔克2,3、沢田康次2 1．東北学院大学教養学部、2．東北工業大学、3．NPO法人naturalscience
エチゼンクラゲの工業材料への有効利用研究 勝然 勇也・三浦 拓也・諏訪 勇麻・阿部 令・佐々木 卓博・鈴木 真・遠藤 雄飛・小野 浩明・東 兼規・山田 有理・伊藤 大将・小川原 佳希・今野 優也・佐々木 祐香・門馬 裕里・山内 透也 宮城県石巻工業高等学校
ヒトの運動制御システムの解析 ～ヒトによる倒立振子の制御～ ○佐瀬一弥1,3、林叔克2,3 1．東北大学工学部、2．東北工業大学、3．NPO法人naturalscience
領域２ 津波と地球科学 今村文彦，越村俊一，後藤和久，菅原大助，萱場真太郎，高橋潤，橋本貴之，川俣秀樹 東北大学大学院工学研究科附属災害制御研究センター
津波と地球科学―恐竜絶滅の隕石衝突と地球史上最大の津波― ○後藤和久1、菅原大助2 １ ． 東北大学大学院工学研究科、 ２ ． 東北大学大学院理学研究科
地層に残る津波の証拠―仙台平野を襲った平安時代の津波― ○菅原大助1 、今村文彦2 、後藤和久2 1．東北大学大学院工学研究科、2．東北大学大学院理学研究科
領域３ 逃げるカタツムリと追いかけるヘビ ○細将貴 東北大学大学院生命科学研究科
小数の働きアリが示す活動のリズム ○結城麻衣1,4、菅原研1、林叔克2,4、菊地友則3、辻和希3 1．東北学院大学教養学部、2．東北工業大学、3．琉球大学農学部亜熱帯動物学講座、4．NPO法人naturalscience
女王蟻によるパトロール行動のモデル化とパトロールの最適化問題 ○八重樫和之1,5、林叔克2,5、菊地友則3、菅原研4、辻和希3 1．東北大学工学部、2．東北工業大学、3．琉球大学農学部亜熱帯動物学講座、4．東北学院大学教養学部、5．NPO法人naturalscience
↑ このページのトップへ

■ 学都「仙台・宮城」サイエンス・デイ２０１０

２００９年の開催報告

■ 学都「仙台・宮城」サイエンス・デイ２００９

■ 2009年７月５日（日）
サイエンス・フェスティバル


２００８年の開催報告

■ NPO法人 natural science について