Frawieście Falowe Epidemii: Model SIR w Przestrzeni |… |…

Siedzimy w miejscu przyzwyczajonym

Klasyczny model SIR dzieli populację na kompartmenty Zakażonych, Chorych i Wyzdrowiałych i śledzi, jak przepływają one do siebie w czasie, bez uwzględnienia geografii – zakłada się, że każdy mieszka z każdym równie. Rzeczywiste epidemie nie zachowują się w ten sposób: zaczynają się w jednym punkcie i rozprzestrzeniają na zewnątrz, miasto po mieście, ulica po ulicy. Dodanie wymiaru przestrzennego przekształca zwykłe równania różniczkowe w równania reakcji-difuzji cząstkowe, a rezultat to prawdziwa fala podróżująca infekcji rozprzestrzeniająca się przez populację z przewidywalną prędkością.

demo na żywo · powiązana symulacja● LIVE

Terminy reakcji Kermacka-McKendrick

W każdym punkcie przestrzeni lokalne dynamiki nadal przestrzegają oryginalnych równań z 1927 roku Kermacka-McKendrick – nowe zachorowania proporcjonalne do iloczynu gęstości populacji podatnej i zainfekowanej, a wyzdrowienia proporcjonalne do gęstości populacji zainfekowanej:

dS/dt = -β * S * I + D * Laplacian(S) dI/dt = β * S * I - γ * I + D * Laplacian(I) R/dt = γ * I

β = współczynnik transmisji γ = współczynnik leczenia D = współczynnik dyfuzji przestrzennej R0 = β / γ (liczba podstawowego namnażania)

dS/dt = -beta * S * I  +  D * Laplacian(S)
dI/dt =  beta * S * I  -  gamma * I  +  D * Laplacian(I)
dR/dt =  gamma * I

beta  = transmission rate    gamma = recovery rate    D = spatial diffusion coefficient
R0 = beta / gamma   (basic reproduction number)

Dlaczego przód porusza się z prędkością stałą

Kiedy uformowana jest fale rozchodzącej się, ustala się ona w kształcie, który przemieszcza się bez zmiany – rozwiązanie typu fali podróżującej. Minimalna prędkość, jaką taka fala może utrzymać, ma czyste, zamknięte wyrażenie dla liniowej równania Fishera-KPP: c ≈ 2·√(D·r), gdzie r jest efektywnym lokalnym tempem wzrostu infekcji ((β·S₀ − γ) w pobliżu przedniej krawędzi) a D jest współczynnikiem dyfuzji. Podwojenie współczynnika dyfuzji zwiększa prędkość frontu o około 41% (√2), podczas gdy podwojenie netto tempa wzrostu ma ten sam wpływ – zależność pierwiastkowa oznacza, że zmniejszanie surowej transmisji jest generalnie bardziej efektywnym mechanizmem niż próba spowolnienia fizycznego mieszania, choć realne interwencje zazwyczaj atakują oba.

R0, odporność zbiorowa i kurczący się front

Podstawowy współczynnik rozrodu R₀ = β/γ określa, czy wybuch epidemii jest możliwy – poniżej 1, każdy lokalny wybuch uśmierza się przed powstaniem frontu; powyżej 1, powstaje fala i się ona rozwija. Wraz z pochłanianiem podatnych osób za rosnącym frontem, efektywny współczynnik rozrodu R_eff = R₀·S/N spada, ponieważ liczba podatnych osób dostępnych dla każdego zainfekowanego człowieka maleje. Odporność zbiorowa to próg, w którym R_eff spada do 1, przy czym odsetek populacji odpornej wynosi około 1 − 1/R₀: w tym momencie fala przestaje się utrzymać i gaśnie, nawet jeśli wielu osób pozostaje niezakażonych, ponieważ odsetek odpornych osobników rozdziera ciągi transmisji, na których polega front.

Dlaczego rzeczywiste wybuchy nie są gładkimi okręgami

Model jednorodnej dyfuzji przewiduje uporządkowany rozszerzający się okrąg, ale faktyczne fronty epidemii są widocznie nieregularne z powodu transmisji zachodzącej natychmiastowo na duże odległości – pojedynczy zarażony podróżnik może wywołać nowy front setki kilometrów dalej, znacznie wcześniej niż lokalna fala dyfuzyjna. Ta dalekosiężna interakcja, nakładająca się na lokalną dyfuzję, jest powodem, dla którego modele rzeczywistych epidemi (takich jak historyczne zarazy lub współczesne rozprzestrzenianie się grypy) lepiej pasują do mieszanki krótkotrwałej dyfuzji i sieci długodystansowych skoków wzdłuż tras transportowych, zamiast dyfuzji samej w sobie; wyjaśnia to również, dlaczego ograniczenia podróży celowo ukierunkowują się na termin dalekiego zasięgu, a lokalna dyfuzja jest w porównaniu trudniejsza do zatrzymania.

Frequently asked questions

Dlaczego infekcja rozprzestrzenia się falą zamiast natychmiast wszędzie?

Wynika to z tego, że transmisja jest lokalna – zarażona osoba głównie zaraża osoby w pobliżu, które są podatne na chorobę, a nie całe populacje jednocześnie. Ten lokalny związek, matematycznie wyrażony jako człon dyfuzji dodany do równań SIR, generuje falę, która wywija się na zewnątrz od pierwotnych przypadków z prędkością w przybliżeniu stałą, zamiast infekcji pojawiać się równomiernie na mapie.

Co determinuje prędkość ruchu fali frontu?

Dla postaci SIR o reakcji-dyfuzji, minimalna zrównoważona prędkość fali jest proporcjonalna do pierwiastka kwadratowego z współczynnika dyfuzji pomnożonego przez lokalną szybkość wzrostu (beta*S - gamma). To skalowanie pierwiastkiem kwadratowym oznacza, że prędkość reaguje bardziej powoli na zmiany w którymkolwiek z tych parametrów niż mogłoby się intuicyjnie spodziewać.

Czy odporność zbiorowa oznacza, że fala epidemii zatrzymuje się natychmiast?

Oznacza to, że fala nie może utrzymać się, gdy liczba przenoszeń (efektywny współczynnik reprodukcji) spada do 1, ale fala nie zanika natychmiast – zwalnia i pozostała populacja chorych maleje w czasie, podobnie jak dowolny oscylator z dodanym tłumieniem, który wycisza się, a nie zatrzymuje się nagle.

Wypróbuj na żywo

Wszystko powyżej działa bezpośrednio w Twojej przeglądarce — otwórz Epidemic Wavefront i zmieniaj parametry podczas działania. Nic nie jest instalowane ani przesyłane na serwer, cały model działa w jednej karcie.

▶ Otwórz symulację Epidemic Wavefront

Frawieście Falowe Epidemii: Kiedy SIR Zyskuje Mapę