登革熱的傳播研究

劉敏 郭翔宇 徐樂希

摘 ?要：本文主要針對登革熱的傳播進行了相關研究。首先我們估算了2017年廣州爆發月份，通過建立的模型以及對題目中的相關數據進行插值擬合處理，從而計算出了c值，并且用布雷圖指數檢驗進行了驗證。我們又收集整理了2017年廣州登革熱的發病情況，并且根據福州市1694例登革熱的臨床分析，估算出登革熱感染者的平均患病天數是9.4d，并且通過SIR模型以及MATLAB的畫圖驗證了取值的合理性。然后根據廣東的案例數和相關資料，我們估算出2017年廣州爆發月份為8月，之后將題目中的數據進行擬合插值處理，在不考慮蚊子的雌雄比例的情況下，通過MATLAB估算出了c值為0.047。若考慮到雌雄之間的比例關系，參考相關文獻后我們選取雌雄比例3：1進行計算，此時估算的c值為0.062。

關鍵詞：登革熱;SIR模型;布雷圖指數檢驗

1.引言

登革熱（俗稱斷骨熱）是一種由登革熱病毒引起的急性發熱傳染病，其通過帶有登革熱病毒的雌性伊蚊叮咬而傳染給人類，主要癥狀是發熱、頭痛、胃痛、肌肉痛或關節痛，臨床表現為高熱、頭痛、肌肉、骨關節劇烈酸痛等[1]。由于蚊媒在登革熱的傳播中占據主導地位，研究表明蚊子在不同的溫度和相對濕度下對登革熱病毒的傳播能力有所不同[2]。因此，溫度和相對濕度在決定登革熱是否爆發的過程中起到了重要的作用[3]。

2.2 數據處理

首先我們根據文獻查閱到廣州登革熱爆發的月份在8、9、10月，由于廣州作為廣東省會，病例數約占廣東省70%，因此我們用廣東省的數據來確定廣州的爆發月份。由爆發的定義知：某地區某病在短時間內（一般以小時、天、周或月計算）發病數突然增多時稱爆發，所以我們根據廣東的每月案例的增長率，來判斷爆發的月份。從數據來看，6月登革熱案例有19例，7月登革熱案例有23例，8月登革熱案例有139例，9月登革熱案例有332例，10月登革熱案例有701例。7月登革熱案例的增長率為0.21，8月登革熱案例的增長率為5.04，9月登革熱案例的增長率為1.38，10月登革熱案例的增長率為1.11，8月是增長率最高的月份。從時間橫向比較來看，8月為爆發月份也是合理的。

根據收集整理的2017年的平均氣溫和平均相對濕度，我們利用MATLAB分析每個月數據的插值，之后用MATLAB進行數據處理得到 。若 不考慮雌雄蚊子的比例，我們得出c=0.047，若考慮雌雄蚊子比例3：1，我們得出的c=0.062。

2.3 模型檢驗

可知要想求得C值，必須先要求出M的值。根據M定義—載體數與人類數量的比例得：必須求出廣州市每月蚊子數目和廣州市人口數。而布雷圖指數正是求解蚊子數量的重要指標。

通過查閱廣州市疫情防控中心給定資料得到5-12月份每周給定的平均布雷圖指數，再通過求平均的方法計算得出每月平均布雷圖指數。

根據布雷圖指標定義，必須求得2017年廣州家庭戶數。但由于家庭戶數數據缺失，（可以不要）經過查閱廣州市2018年統計局統計年鑒得出，廣州市2017年人口數n人為1449.84萬人，平均每戶2.75人，易得出全市家庭戶數n戶為5272145戶。題目給予簡化不予考慮短時間人口變化，因此每月人口及家庭數可當做一個定值處理。

利用上表所給數據，進行MATLAB插值，求得5-12月份 與 ;由公式： 。在廣州市2017年人口數n人一定情況下，此通過M定義，我們得出5-12月份每個月的M。明確M，我們只需要對公式： 簡單變形，即可得出c值為0.04962。

通過簡單計算，發現c值與所建模型所求c值有一定波動，但大體相似，可間接說明模型所求c值的正確性。

綜合考慮廣州市各個區之間，不同的地區的氣象因素對登革熱的爆發情況不盡相同，可能由于不同地區登革熱的爆發規律與本地區的特殊地理位置、氣候條件、政策原因等有關，所以考慮到登革熱的爆發時間可能存在滯后效應，因此預測的爆發月份與實際廣州爆發疫情的情況是趨于一致的。

參考文獻

[1] ?王藝瑋，陳鳴浩，朱瑾郁，林支桂. 大型滅蚊行動對2014年廣州市登革熱疫情的影響[J]. 揚州大學學報（自然科學版），2017，20（02）：10-12+41.

[2] ?江毅民，嚴子鏘，胡志剛，徐建敏，李成玲，梁雪瑩. 誘蚊誘卵指數作為蚊蟲密度監測指標適用性的探討[J]. 熱帶醫學雜志，2015，15（11）：1555-1557.

[3] ?嚴子鏘，胡志剛，江毅民，伍浩榮，潘志民. 標準間指數在白紋伊蚊自然種群動態監測中的應用研究[J]. 中華衛生殺蟲藥械，2006（05）：352-354.