1 UVOD Nestabilnost prometnega toka na avtocestah povzroča zgostitve v prometu in nastajanje zastojev, ki povzročajo veliko ekonomsko izgubo [1]. Zato je veliko truda usmerjenega v razvoj novih metod za stabiliziranje prometnega toka [1-3]. S tem ciljem v nadaljevanju empirično izpeljemo novo različico osnovnega zakona prometa, ki opisuje vpliv gostote prometa na hitrost in pretok vozil na avtocesti [2-4]. Na podlagi značilnosti ustreznih diagramov nato opišemo nestabilni režim prometnega toka in razložimo možnost za njegovo stabiliziranje. Za ta namen predlagamo optimalno prilagajanje omejitve hitrosti obstoječi gostoti prometa, kar nas vodi do opredelitve pripadajočega kontrolnega zakona. Z njim je opredeljeno, kako je treba prilagajati omejitev hitrosti, če želimo zagotoviti optimalno pretočnost in hkrati s tem stabilen prometni tok. 2 NOVI OSNOVNI ZAKON PROMETA Z naraščanjem gostote vozil 𝜌 prometni tok 𝑄 najprej narašča do neke maksimalne vrednosti 𝑄𝑚𝑎𝑥, nato pa začne znova upadati. V prvem režimu prostega toka je stanje prometa stabilno, medtem ko je v drugem režimu zgoščenega toka nestabilno [1-3]. Ker je položaj maksimuma odvisen od vrednosti omejitve hitrosti 𝑣𝑜, je mogoče doseči stabilni režim s primerno prilagoditvijo te vrednosti gostoti vozil na cesti. Vrednost 𝑣𝑜 je smiselno postaviti tako, da se pri dani gostoti vzpostavi stabilno stanje, ki ustreza maksimalnemu pretoku vozil. Do ustreznega zakona za vodenje pridemo z analizo osnovnih diagramov prometa v nadaljnjem besedilu. Obstoječa literatura o dinamiki prometa [1-3] vsebuje veliko različic fizikalnega opisa lastnosti prometa in ustreznih diagramov. Najosnovnejša je različica, ki podaja odvisnost povprečne hitrosti vozil 𝑣 od gostote vozil ρ. Iz izkušenj vemo, da povprečna hitrost z naraščajočo gostoto vozil upada; hkratne meritve obeh spremenljivk pa to potrjujejo tudi kvantitativno. Na sliki 1 so z oznako –o– prikazani eksperimentalni podatki iz reference [3] za dvopasovno avtocesto z dovoljeno hitrostjo 120 km/h, ki zaradi prisotnosti tovornjakov ustreza nazivni omejitvi hitrosti 𝑣𝑜 = 110 km/h. Na podlagi značilnosti tega grafa so že prej opredelili [3-str.112] prvi osnovni zakon prometa, ki opisuje odvisnost hitrosti od gostote z nastavkom: 𝑣(𝜌) = 𝑣𝑜[1 − (𝜌/𝜌𝑘) α]𝛽. V njem pomeni 𝜌𝑘 največjo mogočo gostoto in α ter β empirično opredeljiva parametra. Čeprav ta nastavek kar dobro opiše najbolj Prejet 30. maj, 2017 Odobren 30. junij, 2017 168 ŠVEGL, GRABEC Slika 1: Odvisnost hitrosti 𝑣 od gostote 𝜌 za omejitev hitrosti 𝑣𝑜 = 110 km/h. (–o–) eksperimentalni podatki iz reference [3- str. 259], (…) st-model, (---) p-model. izstopajoče lastnosti prometa, pa takšna opredelitev nastavka ne ponuja ugodne analitične možnosti za vključitev v celotni opis nekaterih znanih lastnosti prometa, ki vodijo do nestabilnosti prometnega toka pri velikih gostotah. Zato je bilo že veliko raziskav usmerjenih v opredelitev osnovnega zakona prometa tudi s pomočjo statistične analize dinamike prometa. Kljub uspehu na tem področju [2,3], pa je takšen pristop precej kompleksen in težko prilagodljiv za naš namen. Zato smo se odločili za srednjo pot ter poskusili priti do opisa prvega osnovnega zakona tako, da upoštevamo, kako vozniki prilagajajo hitrost svojih vozil obstoječemu stanju prometa na avtocesti [2]. Kot osnovo za izpeljavo smo tudi mi uporabili eksperimentalne podatke, o prvem osnovnem zakonu prometa, ki so predstavljeni na sliki 1 z grafom (–o–). Značilnosti prometa na avtocesti so v glavnem opredeljene z omejitvijo hitrosti 𝑣𝑜, in primerno hitrostjo: 𝑤 = 𝑣𝑜(𝑟 – 𝑟𝑘) / 𝑡. V tem izrazu opisuje 𝑡 reakcijski čas voznika, 𝑟 pomeni razdaljo med vozili, 𝑟𝑘 pa je razdalja, ki je malo večja kot povprečna dolžina vozil in ustreza popolnoma zgoščenemu stanju, kjer se promet ustavi. Ta razdalja je bistveno manjša od razdalje 𝑟𝑜 = 𝑡 𝑣𝑜 , ki jo vozilo opravi v reakcijskem času 𝑡 pri največji dovoljeni hitrosti 𝑣𝑜. Hitrost 𝑤 se imenuje »primerna« zato, ker omogoča vozniku v kritičnih razmerah ustaviti vozilo na območju proste ali varnostne razdalje med vozili 𝑟 – 𝑟𝑘 v intervalu reakcijskega časa 𝑡 . Razdalja med vozili določa gostoto 𝜌 = 1/𝑟 , medtem ko 𝑟𝑘 določa največjo mogočo gostoto 𝜌𝑘 = 1/𝑟𝑘 pri kateri promet obstane. Gostota in primerna hitrost sta torej povezani z enačbo 𝑤 = (1/𝜌 – 1/𝜌𝑘) / 𝑡 . V prostem toku pri majhni gostoti prometa vozniki v glavnem prilagajajo hitrost mejni vrednosti opredeljeni z omejitvijo hitrosti 𝑣𝑜, medtem ko pri zgoščenem prometu z zastoji, vozijo približno s primerno hitrostjo 𝑤. V skladu s temi lastnostmi lahko grobo opišemo odvisnost hitrosti od gostote z enačbo: 𝑣 = min(𝑣𝑜 , (1/𝜌 – 1/𝜌𝑘) / 𝑡) (1) Ta izraz pomeni stopničasti model (st), ker je prehod iz ene vrednosti v drugo opisan s stopničasto funkcijo: min(𝑣𝑜, 𝑤) . Pripadajoči graf je prikazan na sliki 1 s točkasto črto (…) ob eksperimentalnih podatkih, ki so predstavljeni s krivuljo (-o-) [3,4]. Stopničasti model kar dobro opisuje nekaj osnovnih lastnosti prometa pri veliki ali pa pri majhni gostoti [2]. Očitno odstopanje pa je opazno pri prehodu iz omejitve hitrosti 𝑣𝑜 v primerno hitrost 𝑤, ko je 𝑣𝑜 = 𝑤 = (1/𝜌𝑝 – 1/𝜌𝑘) / 𝑡. Iz te povezave dobimo izraz za gostoto prehoda 𝜌𝑝 = 1/(𝑟𝑜 + 𝑟𝑘). Ker označuje vrednost 𝜌𝑝 prehod iz prostega toka pri 𝜌 < 𝜌𝑝 v zgoščeni tok pri 𝜌 > 𝜌𝑝, je vrednost 𝜌𝑝 pomemben parameter za obravnavanje problemov, povezanih s stabiliziranjem prometnega toka. Zato, da bi prilagodili stopničasti model eksperimentalnim podatkom predpostavimo, da se prehod iz ene v drugo omejevalno vrednost hitrosti ne zgodi pri točno določeni vrednosti 𝜌𝑝, temveč v intervalu okoli nje. To predpostavko lahko upoštevamo, če vsako omejevalno komponento opišemo s primerno izbrano potenco inverzne vrednosti pripadajoče hitrosti, nato pa oba omejevalna pogoja preprosto seštejemo pri opredelitvi izraza za ustrezno hitrost: 1/𝑣𝑝 = 1/𝑣𝑜 𝑝 + 1/𝑤𝑝, ki nas nato vodi do izraza za opis odvisnosti hitrosti od gostote. Tukaj potenca 𝑝 opisuje kako nenadno opravijo vozniki prehod iz ene omejevalne hitrosti v drugo; čim večji je 𝑝, tem ožji je interval prehoda. Vrednost parametra 𝑝 najpreposteje določimo na podlagi eksperimentalnih podatkov. V ta namen preoblikujemo izraz, ki opisuje omejevalni pogoj v obliko: 𝑣 = 𝑣𝑜 (1 + (𝑣𝑜/𝑤) 𝑝 )−1/𝑝. Iz lastnosti uporabljenih funkcij izhaja, da hitrost 𝑣 ne presega nobene od omejevalnih komponent 𝑣𝑜, 𝑤. Ker je primerna hitrost 𝑤 odvisna od gostote vozil 𝜌, je od 𝜌 odvisna tudi hitrost: 𝑣 = 𝑣(𝜌). Ustrezna funkcija, ki opisuje tako imenovani potenčni, oziroma p-model, je nato opisana z izrazom: 𝑣(𝜌) = 𝑣𝑜 {1 + [ 𝑟𝑜 1 𝜌⁄ + 1 𝜌𝑘⁄ ] 𝑝 } −1/𝑝 (2) Ta funkcija, ki opisuje prvi osnovni zakon prometa, je podobna prej omenjenemu empirično pridobljenemu nastavku in je navidezno zapletena; vendar je dejansko dokaj preprosta, ker vključuje samo potence in osnovne aritmetične operacije, ki dajejo možnost opredelitve inverzne funkcije in tudi preproste nadaljnje analize značilnosti prometa. Prednost tega modela je, da vključuje zgolj en prilagodljivi ali prosti parameter p, medtem ko so 𝑣𝑜, 𝑟𝑜 = 𝑡𝑣𝑜 in 𝜌𝑘 = 1/𝑟𝑘 objektivno določeni z lastnostmi vozil, voznikov in same ceste. Veliko predhodno objavljenih oblik osnovnih zakonov ne kaže teh ugodnih značilnosti [1-3]. MOŽNOST STABILIZIRANJA PROMETNEGA TOKA S PRILAGAJANJEM OMEJITVE HITROSTI 169 Slika 2: Odvisnost hitrosti 𝑣 od gostote vozil 𝜌 za omejitev hitrosti 50 km/h ≤ 𝑣𝑜 ≤ 150 km/h. (-o-) eksperimentalni podatki iz reference [3], (---) p - model. Parametre v p-modelu izberemo tako, da ustrezna krivulja sovpada z eksperimentalno pridobljeno, ki je objavljena v vodilni literaturi s tega področja [3] in prikazana s krivuljo (-o-) na sliki 1. Pripadajoče vrednosti parametrov modela so: 𝑟𝑘 = 7m, 𝑡 = 1,2s in 𝑝 = 2,5. Diagram, prikazan na Sliki 1 s polno črto kaže, da naš p-model presenetljivo dobro opisuje osnovne značilnosti prometa na avtocesti, oziroma njegov prvi osnovni zakon pridobljen z meritvami. Za svojo nadaljnjo analizo uporabimo p-model pri različnih vrednostih omejitve hitrosti 𝑣𝑜 na intervalu od 50 do 150 km/h z razmikom 𝑑𝑣𝑜 = 10 km/h in prikažemo ustrezne diagrame na sliki 2. Z enačbo: 𝜌 = 1/𝑟 pretvorimo sliko 2 v sliko 3, ki prikazuje povezavo med razdaljo 𝑟 in hitrostjo 𝑣. V prostem toku, ko je razdalja med vozili 𝑟 velika, je hitrost približno enaka omejitvi hitrosti 𝑣𝑜. Z manjšanjem razdalje med vozili pa se hitrost skoraj linearno zmanjšuje z 𝑟, kakor je določeno s primerno hitrostjo 𝑤 = (𝑟 – 𝑟𝑘)/𝑡. Slika 3: Odvisnost hitrosti 𝑣 od razdalje 𝑟 za omejitev hitrosti na intervalu 50 km/h ≤ 𝑣𝑜 ≤ 150 km/h. Slika 4: Odvisnost prometnega toka 𝑄 od gostote 𝜌 pri omejitvi hitrosti 𝑣𝑜 = 110 km/h. (–o–) eksperimentalni podatki iz reference [3], (…) st-model, (---) p-model. Če pomnožimo enačbo (2) z gostoto 𝜌, dobimo drugi osnovni zakon prometa, ki opisuje odvisnost prometnega toka 𝑄 od gostote 𝜌: 𝑄(𝜌) = 𝑣𝑜𝜌 {1 + [ 𝑟𝑜 1 𝜌⁄ + 1 𝜌𝑘⁄ ] 𝑝 } −1/𝑝 (3) Sliki 1 in 2 se nato spremenita v sliki 4 in 5. Slika 4 kaže tri grafe, ki ustrezajo eksperimentalnim podatkom, st-modelu, in p-modelu. Kakor v prvi verziji zakona tudi v tej, drugi zelo dobro sovpadata diagrama, ki opisujeta eksperimentalne podatke in p-model, medtem ko st- model s svojo trikotno obliko nakazuje predvsem osnovne lastnosti naraščanja in upadanja prometnega toka z naraščajočo gostoto vozil [2]. Iz izraza za gostoto prehoda 𝜌𝑝 = 1/(𝑟𝑜 + 𝑟𝑘) sledi, da nastopi maksimum prometnega toka 𝑄𝑚𝑎𝑥 v st-modelu pri gostoti: 𝜌𝑝 = 𝜌𝑘 1+ 𝑡𝑣𝑜𝜌𝑘 (4) Ta enačba kaže da se z manjšanjem vrednosti dovoljene hitrosti 𝑣𝑜 veča vrednost gostote 𝜌𝑝, oziroma veča širina stabilnega območja prometa. To je pomembna lastnost, ki že nakazuje, kako je treba spreminjati omejitev hitrosti, če želimo zagotoviti stabilne razmere v prometu! Do podobne lastnosti pridemo tudi z analizo p-modela. Iz enačbe 3 dobimo z upoštevanjem pogoja za ekstrem 𝑑𝑄/𝑑𝜌 = 0 naslednji izraz za gostoto, ki ustreza maksimumu toka: 𝜌𝑚𝑎𝑥 = 𝜌𝑘 1+[𝑡𝑣𝑜𝜌𝑘] 𝑝 (𝑝+1)⁄ (5) Ta enačba se do potence 𝑝 (𝑝 + 1)⁄ ≈ 0,7 v imenovalcu ujema z enačbo 4, kjer ima vrednost 1 , zato so tudi zaključki, ki sledijo iz nje podobni, malo se razlikujejo samo vrednosti ocenjene iz obeh enačb. 170 ŠVEGL, GRABEC Slika 5: Odvisnost prometnega toka 𝑄 od gostote 𝜌 za omejitev hitrosti 50 km/h ≤ 𝑣𝑜 ≤ 150 km/h. Iz izraza za gostoto, ki ustreza ekstremu, lahko dobimo s pomočjo izraza 𝑄(𝜌) = 𝜌 𝑣(𝜌) nato še izraz za samo vrednost maksimalnega toka. Iz st-modela dobimo za maksimalni tok izraz: 𝑄𝑚𝑎𝑥 = (1 – 𝜌𝑝/ 𝜌𝑘) / 𝑡, ki kaže, da optimalna vrednost skoraj linearno pada s širino stabilnega režima. Podobne lastnosti prometa, kakor jih dobimo iz analize st-modela, dobimo tudi z analizo p-modela, ki ga opisuje enačba 3, vendar so pri tem dobljene malo drugačne vrednosti spremenljivk, ki opisujejo ekstrem. Ker je enačba 3 za hitro oceno lastnosti malo preobširna, je vpliv omejitve hitrosti na potek 𝑄(𝜌) in položaj maksimumov prikazan z naborom diagramov na sliki 5. Pri veliki vrednosti 𝑣𝑜 nastopi maksimum 𝑄𝑚𝑎𝑥 pri nižji vrednosti gostote 𝜌𝑚𝑎𝑥 (zgornji diagram) kot pa pri nizki vrednosti 𝑣𝑜 (spodnji diagram). To pomeni, da je prometni tok manj stabilen pri visokih vrednostih dovoljene hitrosti 𝑣𝑜 kakor pa pri nizkih vrednostih. V skladu s tem moramo zmanjšati vrednost dovoljene hitrosti, če hočemo razširiti območje stabilnosti. Ta navidezno presenetljiv zaključek se nam lahko zdi paradoksalen, vendar izvira iz analize eksperimentalnih podatkov [3]! 4 STABILIZIRANJE PROMETNEGA TOKA Z namenom, da pridemo do napotka za optimalno prilagajanje omejitve hitrosti obstoječi gostoti prometa 𝜌, določimo najprej za st-model inverzno funkcijo funkciji, podani v enačbi 4: 𝑣𝑜 = (𝜌𝑘/𝜌 – 1)/(𝑡𝜌𝑘) (6) Ta enačba nam kaže kako moramo vrednost dovoljene hitrosti 𝑣𝑜 zmanjševati z naraščajočo gostoto, če želimo doseči pri dani gostoti prometa optimalni pretok. Hitrost pri maksimumu pada z naraščajočo gostoto podobno kakor primerna hitrost 𝑤. Isto velja tudi za p-model, le da je funkcija, ki sledi iz enačbe 5, za malenkost drugačna: Slika 6: Odvisnost omejitve hitrosti 𝑣𝑜 od gostote 𝜌𝑚𝑎𝑥. 𝑣𝑜 = (𝜌𝑘/𝜌 – 1) (𝑝+1) 𝑝⁄ /(𝑡𝜌𝑘) (7) Enačbi 6 in 7 kažeta iskani optimalni zakon za prilagajanje omejitve hitrosti k obstoječi gostoti prometa 𝜌 , ki sledi iz st- oziroma p-modela. Ustrezni funkciji sta prikazani na sliki 6 in kažeta prej omenjeno presenetljivo lastnost izrazitega padanja vrednosti omejitve hitrosti 𝑣𝑜 s povečevanjem gostote 𝜌 . Funkcija, podana v enačbi 7 je uporabna kot osnova za stabiliziranje prometnega toka z optimalnim prilagajanjem omejitve hitrosti gostoti prometa. Za ta namen moramo najprej izmeriti gostoto vozil na cesti in nato njej prilagoditi tisto vrednost omejitve hitrosti, pri kateri prometni tok doseže maksimum. Za stabiliziranje prometnega toka mora biti nato vrednost optimalne omejitve hitrosti posredovana voznikom z opozorilnimi tablami prometne infrastrukture. Smiselno pa je posredovati voznikom po mobilnih telefonih tudi informacijo o razmerah v samem prometu in kaj bodo s prilagoditvijo spremenjeni vrednosti omejitve hitrosti dosegli. V skladu s podatki, objavljenimi v literaturi, lahko pričakujemo, da bi taka kontrola pripomogla k stabilnosti prometa [1,2]. Učinek lahko pričakujemo v bolj homogenem toku z manj zgostitvami in zastoji, v skladu s tem pa naj bi se skrajšali tudi potovalni časi. V centrih za nadzor in vodenje prometa se po navadi zanimajo tudi za maksimalno vrednost homogenega toka 𝑄𝑚𝑎𝑥 pri dani omejitvi hitrosti, ki jo imenujemo kapaciteta ali zmogljivost ceste. Za ta namen lahko uporabijo diagram na sliki 7, ki prikazuje zmogljivost ceste v odvisnosti od omejitve hitrosti 𝑄𝑚𝑎𝑥(𝑣𝑜). Zanimiva je tudi odvisnost maksimalnega toka od dane gostote 𝑄𝑚𝑎𝑥(𝜌), ki je prikazana na sliki 8. Pripadajoča funkcija je za st-model opisana kar s premico: 𝑄𝑚𝑎𝑥(𝜌) = (1 − 𝜌/𝜌𝑘)/𝑡 med začetno in končno točko. Oba podatka sta koristna predvsem pri preusmerjanju prometa ob kritičnih dogodkih na avtocesti. Za nadzornike prometa je zanimiva tudi odvisnost povprečnega toka < 𝑄 > od omejitve hitrosti, ker je < MOŽNOST STABILIZIRANJA PROMETNEGA TOKA S PRILAGAJANJEM OMEJITVE HITROSTI 171 𝑄 > najlažje merljiv in najpogosteje uporabljani podatek o prometu. Pripadajoči graf kaže slika 9. Slika 7: Odvisnost zmogljivosti ceste 𝑄𝑚𝑎𝑥 od omejitve hitrosti 𝑣𝑜. Slika 8: Odvisnost zmogljivosti ceste 𝑄𝑚𝑎𝑥 od ustrezne gostote 𝜌𝑚𝑎𝑥. Graf na sliki 9 je podoben grafu, ki prikazuje zmogljivost ceste na sliki 7, vendar je vrednost povprečnega toka samo približno polovica zmogljivosti. V strokovni literaturi, ki se nanaša na to temo, je pogosto omenjena povezava med hitrostjo in tokom, ki je prikazana na sliki 10. Ta povezava je pomembna, ker nam posreduje dodatne podatke za oceno prehoda iz stabilnega v nestabilni režim na podlagi meritev toka 𝑄 [1-3]. Osnovna spremenljivka za stabiliziranje toka je gostota 𝜌. Po navadi ni izmerjena neposredno, temveč je izračunana iz povprečnega toka < 𝑄 > in povprečne hitrosti < 𝑣 > z enačbo 𝜌 = < 𝑄 > / < 𝑣 >. Meritve povprečnih vrednosti < 𝑄 > in < 𝑣 > zahtevajo nekaj časa, zato takšni podatki omogočajo samó kontrolo z določenim časovnim zamikom, kar pa lahko zmanjša učinek sáme kontrole. Temu problemu se lahko izognemo z napovedovanjem gostote in omejitve hitrosti za kratek časovni interval vnaprej [4,5]. Metoda napovedovanja daje možnost dolgoročne napovedi za optimalno omejitev hitrosti, ki je nato lahko iz Slika 9: Odvisnost < 𝑄 > od omejitve hitrosti vo. Slika 10. Odvisnost hitrosti v od toka 𝑄 za omejitve hitrosti 50km/h ≤ 𝑣𝑜 ≤ 150 km/h. kontrolnih centrov prometa po različnih komunikacijskih kanalih posredovana voznikom [6]. 5 SKLEP Naš osnovni zakon prometa, ki je podan z enačbo (2), se zelo dobro ujema z merskimi podatki. Z njegovo analizo pridobimo osnovo za opredelitev stabilnega in optimalnega prometnega režima in tudi podatke za stabiliziranje prometa, ki temelji na prilagajanju omejitve hitrosti 𝑣𝑜 obstoječi gostoti prometa 𝜌. Predhodno razvita metoda za napovedovanje prometa [5,6] in tukaj izpeljana ocena za optimalno omejitev hitrosti 𝑣𝑜(𝜌) [7], dajeta podlago za razvoj novih metod za optimalno kontrolo prometa. Za ta namen mora biti predlagana metoda najprej preverjena in potrjena v razmerah na avtocesti, nato pa mora biti še cestna infrastruktura posodobljena z ustreznimi prometnimi znaki in komunikacijskimi napravami. 172 ŠVEGL, GRABEC ZAHVALE To delo je dobilo podporo EraSme MOTRAC in COST Akcije TU1102 ARTS.