31 Mayıs 2020 18:04
Özgür Ertunç, M. Pınar Mengüç ve Reyhan Diz-Küçükkaya*
(Bu makale, sarkac.org'dan alınmıştır)
Bu yazı araştırmanın ilk sonuçlarını aktarmaktadır. Araştırma henüz bilimsel dergilerde hakemlik sürecinden geçmemiştir. Türkiye’de normalleşme sürecinde Covid-19 salgınına karşı alınmış olan tedbirler kaldırıldığında farklı senaryolara göre ikinci bir salgın dalgasının ne zaman ve ne büyüklükte olabileceği tahmin edilmektedir. Semptomsuz olduğu veya hastalığını bildirmediği için resmi kayıtlara geçmemiş hasta sayıları hesaba katılabilse tahminlerin doğruluğu artacaktır. Kestirimler Türkiye için toplam sayılar üzerinden yapıldığından şehirlerarası seyahat yasağının kalkması ile artacak bulaşıklık da beklenen ikinci dalgayı büyütebilir. Ancak bilinmeyen yaygın bir bağışıklık oluşmuşsa ya da toplumun edindiği korunma alışkanlıkları tedbirlerin kalkmasından sonra da etkili olarak devam ederse yeni hasta ve ölüm sayılarında 1 Haziran öncesinde sürmekte olan azalma 1 Haziran’da başlayan normalleşmeye rağmen kalıcı olabilir. |
Günümüzün en önemli sorunu olan Covid-19 salgınının yayılımının önceden tahmin edilmesi ve önlem senaryolarının hızlı tespiti hem ülkelerin hem de organizasyonların hızlı ve etkin bir şekilde önlem alabilmesi için anahtar rol oynuyor. Doğru yapılan simülasyonların öngörüsü sayesinde, devlet kurumları, şirketler ve toplum salgına etkin ve hızlı bir şekilde tepki verip, salgının en az kayıpla geçmesini sağlayabilir. İçinde bulunduğumuz Covid-19 salgınında da, toplumsal yaşamın salgını yok edecek şekilde düzenlenmesi, sağlık hizmetlerinin etkin planlanması, salgın süresince lojistik ihtiyaçların giderilmesi, kamu ve özel kurumların operasyonlarını planlaması için simülasyonlar her ülkede yapılıyor ve önlemler tartışılıyor.
Dünya Sağlık Örgütü’nün (DSÖ) pandemi tanımına göre (Şekil 1) Mayıs sonu itibariyle ülke olarak tepe sonrası (post peak) fazına geçiyoruz. Brezilya, Amerika, Hindistan ve Rusya gibi ülkeler ise 6. fazda bulunuyor. Yani salgın dünyada eş zamanlı yayılmıyor ve benzer şiddette yaşanmıyor.
6. fazdan sonra erken normalleşme, yurt dışından gelen bulaştırıcı hastalar ve virüsün mutasyonla değişmesi ikincil bir dalgaya yol açabilir. Dalgaların arasında aylar olabileceği biliniyor. İçinde bulunduğumuz normalleşme çalışmalarının, salgının gelişimine etkisini görebilmek şu an için acil bir konu. Bu yazıda detaylarını anlatacağımız çalışma, normalleşme ve akabinde yaşanabilecek salgın dalgaların etkisini göstermek ve önlemleri tartışmak için yapıldı. Aynı zamanda, bu çalışma yöneticilere ilerde alacakları kararlara bilimsel olarak yaklaşabilmeleri ve etkin bir yol ve yön stratejisi hazırlamaları için kurgulandı.
Aşağıda detaylarını görebileceğiniz analizde iki senaryo düşünülüyor. Birinci senaryoda normalleşme sürecine Haziran başı, Temmuz başı veya Ağustos başında girildiği takdirde olacaklar tahmin ediliyor. Simülasyonlar beklendiği gibi her üç durumda da ikinci dalganın oluşacağını fakat Haziran’ın normalleşme için hala erken olduğunu gösteriyor. Temmuz veya Ağustos’ta normalleşme sürecine girdiğimiz takdirde ikinci dalga çok daha kontrol edilebilir olacak.
İkinci senaryoda ise Haziran’da normalleşme sürecine girildiği düşünülüyor. Normalleşme sonucunda günlük hasta sayısı Mart’ın ikinci yarısındaki seviyeye ulaştığında sert önlemlerin 1 hafta, 2 hafta ve 3 hafta gecikmeyle alınmasının sonuçları tahmin ediliyor. Simülasyonlar önlem almada 1 haftadan fazla gecikilirse, Mart ayında başlayan salgından daha büyük boyutlu bir salgının yaşanacağını gösteriyor. Bugüne kadar farklı ülkelerde geciken önlemlerin bilançosunu az çok gördük, buradaki hesaplar da benzer bir tabloyu ortaya koyuyor.
Türkiye’deki salgının şimdiki halini (şimdigörü) ve içinde bulunduğumuz normalleşme sürecinin yakın zamandaki etkilerini görebilmek için SEIR (Susceptible Exposed Infective Removed) [4][5][6] türevi bir salgın modeli kullanıyoruz. Bu model bulaşabilen (Susceptible), bulaşan ama bulaştırıcı olmayan (Exposed), bulaştıran (Infective), iyileşen ve vefat eden kişi sayılarının (Removed) zamana bağlı olarak hesabına izin veriyor.
Bu sisteme eklenen ek denklemlerle, hastaneye gelenler, toplam tespit edilenler, hastaneye gelip de iyileşenler, korunanlar ve vefat edenler zamana bağlı olarak hesaplanıyor. Şimdilik 9 türevsel denklemden oluşan modeldeki zamana bağlı integralleri sayısal olarak hesaplayarak salgının gelişimini öngörebiliyoruz. Sağlık sisteminin tanı, karantina ve filyasyon protokollerine göre bu denklem sistemi genişletilebilir.
Bu modelin doğru çalışabilmesi için, hastalığın seyrine dair bazı bilgiler gerekli ve bu verilerin detayı arttıkça modelin doğruluğu da artıyor. Modelin ihtiyaç duyduğu, en önemli parametrelerin bir kısmı hastalığın seyrine dair olanlar. Hastalar bulaşmadan a-gün sonrasında bulaştırıcı oluyorlar. Ancak ilk semptomların ortaya çıkması için, bundan sonra b-gün bir zaman geçmesi gerekiyor. Bu iki zamanın toplamı ise inkübasyon (kuluçka) süresi c olarak tanımlanıyor [6]. Hasta b süresi ne kadar uzun olursa, fark etmeden o süre boyunca çevresine bu hastalığı yayabiliyor. Zaten bu salgının bu kadar yayılmasındaki en önemli etkenlerden birisi de bu. Semptomlar belirdikten sonra, hasta d-gün kadar bulaştırıcı olmaya devam ediyor. Ancak zamanla bulaştırıcılığı azalıyor. Karantinaya alınmadığında toplamda e=b+d gün kadar bulaştırıcı oluyor.
Şimdiye kadar yapılan çalışmalardan belirlenen bu parametreler Şekil 2’de gösteriliyor. Hastaların vefat-vaka oranları da sağlık sisteminin tanı ve tedavi protokollerine bağlı. Sağlık Bakanlığı’nın verilerinden elde edilen vefat-vaka oranlarında bu oranın, Nisan’ın ortasından itibaren yaklaşık %2,5±0,8 (%95 güven aralığı) olduğu tespit edilmiş ve kullanılmıştır.
Şekilde bir hasta için anlık bulaştırma katsayısı BK(t), bulaşıcılığın devam ettiği e gün sonunda 0’a iniyor. BK(t) burada görüldüğü gibi doğrusal olarak azalmayabilir [8], bizim için önemli olan kişinin bir günde ortalama kaç kişiye bulaştırdığını gösteren ortalama bulaştırma katsayısı (OBK). OBK’yı hesaplamak için salgın sürecinde sıkça duyduğumuz ve bir kişinin hastalığı boyunca ortalama kaç kişiye bulaştırdığını gösteren temel üreme katsayısını (basic reproduction number-R0) ve e’yi bilmemiz yeterli.
Bulaşmayı en fazla etkileyen faktörler ve simülasyon tabanlı yapılan optimizasyonun tahmin ettiği parametreler şöyle:
Bilinmeyen parametrelerin bulunması için bir global optimizasyon algoritması olan genetik algoritma kullandık Sağlık Bakanlığı’nın yayınladığı tespit edilen vaka ve vefat sayılarına, simülasyonu en doğru yakınsatan parametre kümelerini bu algoritma yardımı ile bir Pareto kümesi olarak tespit ettik. Pareto kümesi, aslında birçok olası çözümden oluşuyor.
Şekil 3’te yukarıda anlattığımız modeli Sağlık Bakanlığı verileriyle karşılaştırıyoruz. Sol panelde sokağa çıkma yasaklarına göre düzenlenen ve hesaplarda kullanılan bulaştırabilen birey sayısı grafiği var. Orta ve sağ panelde ise günlük tespit edilen vaka ve günlük vefat sayıları için Pareto kümesindeki parametreler kullanılarak bulunan çözümler görülüyor.
Oluşturulan modele göre bu çözümlerin herhangi biri olası bir çözüm. Orta ve sağ panellerde modeli Sağlık Bakanlığı verileriyle karşılaştırdığımızda bazı çözümlerin vefat sayılarını daha iyi hesapladığını, bazı çözümlerin ise tespit edilen vakaları daha iyi hesaplayabildiğini görüyoruz. Çözümlerde görülen dalgalanma, uygulanan sokağa çıkma yasaklarının etkisi.
Bu çözümlerde sağlık sisteminin tespit ettiği her hastanın izole edildiği (karantinaya alındığı) varsayımını yaptığımızı ekleyelim. Ancak Türkiye’de ağır olmayan hastaların evde kendi imkanları ile izole edilerek tedavisinin yapıldığı ve durumunun takip edildiği biliniyor.
Şekil’4'te bu hesaplardan ortaya çıkan günlük bulaşmış ama bulaştırıcı olmayan birey sayısı (exposed), günlük bulaştırıcı birey sayısı ve toplam iyileşmiş birey sayıları görülüyor. Hesaplanan iyileşmiş sayıları ve Sağlık Bakanlığı verileri arasındaki uyumsuzluğun (Şekil 4-sağ panel) tedavi süreçlerinin modelde kabul edilen 7-8 günden daha uzun olabileceğinden kaynaklandığı düşünüyoruz.
Sağ panelde kırmızı ile gösterilen çözümlerdeki 550 bin ve 650 bin arasında değişen toplam iyileşen sayısını aynı zamanda bağışıklık geliştirmiş birey sayısı olarak kabul ediyoruz. Simülasyon, tespit edilmediği halde hastalığı geçirip iyileşen (bağışıklık kazanmış olabilecek) hasta sayısını tespit edilen sayının 3,5-4,5 katı arasında olabileceğini öngörüyor. Şekil 3 ve 4 dikkatlice incelendiğinde, günlük tespit edilmeden dolaşan bulaştırıcı hasta sayısının (Şekil 4-orta panel), tespit edilen vaka sayısın (Şekil 3-orta panel) hastalığın en tepe noktasında 30 katına ulaşabileceği görülebilir. Bu asemptomatik hasta oranının çok yüksek olduğunu gösterir ve bu durum genellikle buz-dağına benzetilir.
Yukarıda gösterdiğimiz simülasyonlar (Şekil 4), Haziran başında, salgının başladığı Mart ayının ortasındaki bulaştırıcı hasta kadar bireyin serbest bulaştırıcı olduğunu gösteriyor. Bu kişiler, Haziran başında gerçekleşecek normalleşme ile ikinci bir dalganın oluşmasının ana nedeni olacaklar. Daha önce serbest bulaştırıcı sayısının salgının gidişatı için çok önemli olduğunu söylemiştik, bu sayı ne kadar büyük olursa dalga da o kadar büyük oluyor. Önlemler daha uzatılırsa serbest bulaştırıcı sayısı düşüyor, ikinci dalganın şiddeti de azalıyor.
İkinci dalganın oluşumunu göstermek ve alınacak önlemlerin bu dalga üstüne etkisini göstermek için iki senaryo düşündük. Her iki senaryoda da normalleşme sürecine girildiğinde bulaşabilen bireylerin sayısını (bağışıklık kazanmış olduğu tahmin edilenler dışındaki) tüm nüfus olarak alıyoruz, ancak öğrenilmiş hijyen/mesafelenme alışkanlıkları ve okulların kapalı olmasını dikkate alarak ortalama bulaştırma katsayısının (OBK), Mart ayı başındakinin yarısı olduğunu varsayıyoruz.
Normalleşme stratejilerinin Haziran, Temmuz veya Ağustos ayları başında gerçekleştirilmiş olacağı üç değişik durum düşünüldü. Her normalleşmeden sonra tespit edilen vaka sayısı salgın başladıktan hemen sonraki (yaklaşık Mart 2020 ortası) seviyesine eriştiğinde, bir hafta gecikmeyle Nisan ayındaki gibi sokağa çıkma yasağı uygulanacağı öngörüldü.
Şekil 5’te görüldüğü gibi Haziran başında gerçekleşecek bir normalleşme, Haziran’ın ikinci haftasından sonra iyice belirgin hale gelen ikinci bir dalgayı oluşturuyor. Devlet, bir hafta gecikme ile Nisan ayındaki gibi yasakları devreye sokarsa Mayıs sonundaki seviyeye ulaşmak için Eylül 2020’yi beklemek gerekeceği görülüyor (Şekil 7’de daha detaylı görülebilir). Buna ek olarak Mayıs sonunda 5 bin seviyesindeki vefat sayısı 8 bin -11 bin aralığına sıçrıyor. Oluşan ikinci dalga, başka ülkeler için yapılan benzeri çalışmalarla paralellik gösteriyor.
Şekil 5 Haziran, Temmuz ve Ağustos başı başlatılacak normalleşmenin yıl boyunca etkisi. Haziran ve Temmuz için iki hafta sonra, Ağustos için 3 hafta sonra oluşan hasta sayısında artışa istinaden, bir hafta gecikmeyle sokağa çıkma yasakları getirilmiş olduğu düşünülerek hesaplar yapılmıştır. Şimdiye kadar açıklanan veri kırmızı dairelerle gösteriliyor, simülasyonlar Mayıs sonuna kadar yapılan (siyah) hesaplamaların verdiği güvenle diğer oluşturuldu.
Eğer normalleşme sürecine Haziran başında değil, Temmuz veya Ağustos ayında girilirse ve hasta sayısının artmasıyla önlemler yeniden devreye sokulursa, sırasıyla Eylül ve Ekim aylarında normalleşme yaşanabilir.
Şekil 6’da, serbest bulaştırıcı ve tespit edilip karantinaya alınan hasta sayılarını belirgin hale getirmek için veri daha dar bir aralıkta gösteriliyor. Görüldüğü gibi asıl sorun, sağlık sisteminin fark edemediği serbest dolaşan hastalar. Temmuz başındaki sayılara bakıldığında, tespit edilen hasta sayısı 20 civarında olduğunda dahi serbest dolaşan ve bulaştıran sayısı 200-1000 arasında olabiliyor. Bu simülasyonlar, Haziran’da yapılacak normalleşmenin sonuçlarına hazırlıklı olmamız gerektiğini ve tespit edilen her vakanın aslında sistemin görüntüleyemediği birçok vakaya eş değer olduğunu bilerek yönetsel kararlar vermemiz gerektiğini gösteriyor.
Haziran başında yaşanacak bir normalleşme için, tespit edilen günlük vaka sayısı yaklaşık 15 Haziran’da Martın ikinci yarısındaki seviyeye geri geliyor. 15 Haziran’dan sonra bir hafta, iki hafta ve üç hafta gecikmeyle Nisan ayında uygulanan sokağa çıkma yasakları uygulandığı takdirde Şekil 7’deki sonuçları elde ediyoruz. Simülasyonlar önlem almada 1 haftadan fazla gecikilirse, Mart ayında başlayan salgından daha büyük boyutlu bir salgının yaşanacağını gösteriyor. Üç hafta gecikme ile getirilecek sokağa çıkma yasakları, vefat sayısını 20 bin bandına taşıma potansiyeline sahip. Böyle bir durum toplumda ve ekonomide öngörülemez etkiler yaratabilir.
Her ne kadar tersini istesek de, bu çalışma, takipçi dalgaların kaçınılmaz olduğunu gösteriyor. Bu dalgaların oluşmasının ana nedeni normalleşme sürecine girerken serbest bulaştırıcı hasta sayısının yeni bir dalgayı tetiklemeye yetecek boyutta olması. Bu dalgalar, sadece Türkiye’de değil tüm dünyada bekleniyor.
Sağlık konusunda yaşanacak bu dalgalanmalar, ekonomik dalgalanmaları hemen tetikleyebilir, etkileri orta vadede güvenlik ve adalette, uzun vadede ise eğitimde görülmeye başlayabilir. İçinde bulunduğumuz durum, fırtınalı bir denizde limana varmaya çalışan bir balıkçı teknesini andırıyor. Bu fırtınalı deniz salgının yarattığı küresel belirsizlik ortamı, balıkçı teknesi ise bizim gibi devletler. Güvenli limana, en az hasarla varmak için, sadece sağlıkta değil, tüm ülkenin ve dünyanın işleyişinde olacak değişikliklerin simülasyonlar ve diğer bilimsel araştırmaların ışığı ile öngörülüp, planlamaların ona göre yapılması sayesinde ulaşılabilir.
Bu çalışmadaki hesaplar, Türkiye’ye özel klinik veriler, filyasyon çalışmalarından ve toplum genelinde yapılacak testlerden elde edilecek asemptomatik hasta oranları ile desteklendiğinde daha doğru sonuçlar verecektir.
Burada yapılan simülasyonlar sonunda hem yöneticiler hem de toplum için aşağıda listelenen sonuçlar çıkıyor:
Özgür Ertunç (Özyeğin Üniversitesi, Makina Mühendisliği Bölümü)
M. Pınar Mengüç (Bilim Akademisi üyesi, Özyeğin Üniversitesi, Enerji Çevre ve Ekonomi Merkezi Direktörü)
Reyhan Diz-Küçükkaya (Bilim Akademisi üyesi, İstanbul Üniversitesi Fen Fakültesi, Moleküler Biyoloji ve Genetik Bölümü)
Bu yazı 30 Mayıs'ta sarkac.org'da yayınlandı
© Tüm hakları saklıdır.