傳染病史料的計量史學研究

周 瑤

傳染病史料的計量史學研究

周 瑤

（黃岡職業(yè)技術學院 醫(yī)藥學院，湖北 黃岡 438002）

應用計量史學的方法，研究中國歷史上瘟疫（即傳染?。┬纬傻挠绊懸蛩丶疤攸c。搜集公元前680年～公元1948年之間的疫情史料，收錄出發(fā)生瘟疫的年份及相關資料，并對其同時存在相關自然因素和社會因素進行計量史學的方法研究和分析。在2628年之間，瘟疫流行年為766年次，占29.15%。經(jīng)二元Logistic回歸分析，瘟疫的流行與氣溫變化（＝0.000，95%的CI：1.516～2.184）、地震（＝0.000，95%的CI：4.274～6.609）、災荒（＝0.016，95%的CI：1.081～2.142）等有明顯相關性，而與戰(zhàn)爭不相關（＝0.056，95%的CI：0.995～1.474）；南方氣候地理條件更易產(chǎn)生瘟疫（X2＝75.787，＝0.000）。瘟疫的防控要結合氣候變化、自然災害、地理環(huán)境，要保護監(jiān)測河流水源、做好邊境檢疫等。

瘟疫流行；傳染病學；計量史學

瘟疫在人類歷史進程中扮演了重要角色。歷史上的瘟疫，即為現(xiàn)代的傳染病。研究瘟疫歷史，并發(fā)現(xiàn)其發(fā)生發(fā)展的規(guī)律，有益于人類了解疾病、對抗疾病，為人類的發(fā)展和健康產(chǎn)生積極的意義。計量史學是一門新興的學科，旨在運用現(xiàn)代數(shù)學、統(tǒng)計學方法和手段，把歷史資料數(shù)據(jù)化，并進行分析，找出規(guī)律，使史料研究精確化和科學化[1-2]。本研究即在此基礎上，對中國歷史上的瘟疫進行分析，得出一些有益結論。

1 資料和方法

1.1 研究資料

疫情史料：《二十四史》《中國古代疫病流行年 表》[3]、全國各省府縣地方志等；地域、戰(zhàn)爭史料：《中國歷史地圖集》[4]《中國戰(zhàn)爭史地圖集》[5]《中國歷代戰(zhàn)爭年表》[6]；地震、氣溫史料：《中國地震歷史資料匯編》[7]等；使用Adobe Flash制作瘟疫發(fā)生發(fā)展時空動畫，輔助判斷疫情傳播的時空關系。

1.2 數(shù)據(jù)錄入及處理

瘟疫：以年為單位，逐年考察當年疫情及相關的情況。

（1）疫情發(fā)生的地域歸屬：結合《中國歷史地圖集》地名，以《基于GIS的中國南北分界的計算和模擬》[8]提出的南北分界線為準，將疫情的流行地域分為三種情況：北方、南方、南北方；（2）疫情流行強度：局限于某一省或州內為“流行”，流行地域大于一省或州為“大流行”；（3）年氣溫升降變化：以竺可楨《中國近五千年來氣候變遷的初步研究》[9]和劉禹等[10]研究結果為依據(jù)，并使用Engauge Digitizer軟件，將氣溫曲線轉換成氣溫數(shù)據(jù)，判斷出當年氣溫的升降幅度；（4）有無戰(zhàn)爭：判斷依據(jù)為《中國戰(zhàn)爭史地圖集》《中國歷代戰(zhàn)爭年表》等；（5）有無地震：《中國地震歷史資料匯編》，以大于5級以上地震為判斷依據(jù)；（6）有無旱災、洪澇等災害：根據(jù)《中國救荒史》[11]及史料記載推斷；（7）瘟疫的發(fā)生發(fā)展：根據(jù)史料記載中時間空間關系判斷，并結合制作的Flash動畫；（8）發(fā)病季節(jié)：根據(jù)史料記載，考慮發(fā)病的潛伏期，以12月～5月為“冬春”季，6月～11月為“夏秋”季。

數(shù)據(jù)處理和統(tǒng)計軟件為Microsoft Excel 2010、SPSS 24.0。統(tǒng)計學方法根據(jù)具體數(shù)據(jù)類型進行統(tǒng)計分析。

2 結果

2.1 古代瘟疫發(fā)生總體情況

結合多種文獻[2,12-20]，收錄自公元前680年至公元1948年的2628年間的數(shù)據(jù)，有記載的瘟疫年份為766年，占比29.15%；其中，流行強度為大流行的年份348年，占比45.43%。

2.2 瘟疫流行的地域規(guī)律

中國歷史隨著經(jīng)濟文化中心由北向南轉移，瘟疫流行也隨著變化。公元前680年至1127年，有記載南、北地域的249年次，其中發(fā)生在北方145年次，南方104年次；公元1127年至1948年，有記載分南北地域的有320年次，其中發(fā)生在北方72年次，南方248年次。數(shù)據(jù)統(tǒng)計學處理表明，隨經(jīng)濟文化中心轉移，瘟疫也相應的發(fā)生轉移（X2＝75.787，＝0.000）。

部分瘟疫有沿長江、黃河、淮河向下傳播的特點：結合制作的Flash動畫觀察時間——地域的疫情發(fā)生發(fā)展過程，部分流行情況總結如下表1：

表1 部分有時間－地域特點的瘟疫蔓延經(jīng)過

雖然不能全面總結瘟疫發(fā)展的特點，但可以看出：在古代交通不發(fā)達的時候，戰(zhàn)爭能加速瘟疫的傳播；很多瘟疫在1～2年內傳遍大半個中國；有些瘟疫有逐水傳播的特點，即沿長江、黃河、淮河向下游傳播。提示水域微生物監(jiān)測有重要意義。如圖1：

圖1 Flash動畫顯示部分瘟疫沿長江向下游傳播（以1706～1709年為例）

各省份頻次：將各省份瘟疫年次進行統(tǒng)計計數(shù)，結果如下圖2：

圖2各省自公元前680年～公元1948年瘟疫發(fā)生的頻次（n＝766）

總體印象是長江、黃河中下游區(qū)域的瘟疫發(fā)生率高，當然也可能與經(jīng)濟發(fā)達、人口眾多有關。其中，環(huán)太湖區(qū)域有記載的瘟疫年份有104次之多，是一個特別高發(fā)的地區(qū)。

2.3 瘟疫的發(fā)生與地震發(fā)生的規(guī)律

統(tǒng)計公元前680年～公元1948年共2628年之間，發(fā)生5.0級以上地震共460年次，有瘟疫288年次，無瘟疫172年次。統(tǒng)計學結果表明，地震的發(fā)生更易引發(fā)瘟疫的發(fā)生（X2＝289.679，＝0.000）。

2.4 瘟疫的發(fā)生與氣溫變化的規(guī)律

將竺可楨的氣溫曲線[9]和劉禹等根據(jù)祁連圓柏樹輪研究的氣溫曲線[10]進行Engauge Digitizer軟件處理，得出氣溫歷年近似數(shù)據(jù)，再與瘟疫發(fā)生進行相關研究，繪圖如下圖3：

圖3 竺可楨氣溫曲線與瘟疫頻次（世紀）關系圖

從每世紀發(fā)生瘟疫頻次與曲線對比關系圖可以看出，當氣溫變化較大時，更容易發(fā)生瘟疫，分別體現(xiàn)在6、12、13、14、16、17、18、19世紀。因此，氣溫進行數(shù)據(jù)處理時，取其變化斜率作為判斷依據(jù)。

2.5 瘟疫與各種災害之間的關系

統(tǒng)計自公元前７世紀至19世紀之間發(fā)生的水災、旱災等[11]，并與瘟疫發(fā)生的頻次進行分析，見表2：

將上表數(shù)據(jù)單獨與瘟疫頻次進行線性回歸分析，值均為0.000，進行多重回歸分析時，對應的災情類別VIF分別為9.850、9.909、13.039、13.081、6.071、4.988、4.982，存在多重共線性問題。

2.6 瘟疫發(fā)生的季節(jié)規(guī)律，以及瘟疫發(fā)生的周期規(guī)律

有記載起始時間為冬春季的155年，占56.98%，夏秋季的117年，占43.01%。

按瘟疫年份發(fā)生的時間順序，進行runs test（Z＝－24.822，＝0.000），說明瘟疫發(fā)生在時間順序上存在一定規(guī)律。再按60年一周期（甲子）進行頻次統(tǒng)計，如圖4：

表2 以世紀為單位災情與瘟疫的關系

圖4 按60年一周期瘟疫發(fā)生頻次統(tǒng)計圖（n＝766）

總體印象是有一定周期規(guī)律的，即約10年左右一個周期，或者4～5年左右一個周期，并非為“庚子年”更容易發(fā)生瘟疫。

2.7 瘟疫的影響因素二元Logistic分析

將瘟疫發(fā)生與否、地震有無、戰(zhàn)爭有無、氣溫聚變（相鄰年氣溫變化＞0. 5%）是否、是否是救荒年[11]等，進行1、0數(shù)據(jù)轉換，進行二元Logistic分析，如下表3。

表3 瘟疫發(fā)生原因二元Logistic分析結果

從上表統(tǒng)計學處理結果來看，地震、氣溫驟變、需要救荒均是導致瘟疫發(fā)生的原因。需要救荒的年份，顯然當年存在各種自然災害。雖然戰(zhàn)爭沒有成為瘟疫發(fā)生的原因，但=0.056非常接近=0.05，可能原因是戰(zhàn)爭的當年與瘟疫的發(fā)生存在一定時間差。通過動畫演示和流行時期內戰(zhàn)爭情況分析，戰(zhàn)爭能加快傳播瘟疫。例如，研究瘟疫流行Flash動畫發(fā)現(xiàn)，207年襄陽之戰(zhàn)、208年赤壁之戰(zhàn)之后，隨著戰(zhàn)爭，首發(fā)于襄陽等地的瘟疫在魏國和吳國蔓延了數(shù)年。

3 討論

瘟疫，或稱疫癘，是中醫(yī)對傳染病的稱呼，既包含病毒感染，也可能是細菌、立克次體等感染。比如公元318年：“比歲有疫病，天行豌豆，癍瘡狀如火燒瘡，皆戴白漿，隨決隨生，不治，數(shù)日必死?！睆呐R床表現(xiàn)來看，應該是水痘或天花之類的疾病。有些病，如公元405年：“十月，大疫，發(fā)赤斑乃愈。”可見，有些傳染病，出過疹子后，即能痊愈，有似麻疹或其他出疹性疾病。有些瘟疫，比如“吊腳痧”，則為細菌性疾病霍亂。

3.1 管理和控制傳染源，適當預判

疫情是否可以預測預判？從研究結果中可以看出，瘟疫伴隨著地震、旱災、水災等災害后出現(xiàn)。俗語“大災之后必有大疫”，所以每次洪水、地震過后，政府及防疫部門均進行大規(guī)模消殺，對大災后的大疫進行了有效的防控。

結果表明，氣溫異常變化是瘟疫發(fā)生的因素，如果出現(xiàn)氣溫異常，也應該進行傳染病的預防，這個結論與吳曉旭等[21]研究結論一致。王穎，梁丁元，李垚等[22]研究也表明，很多病毒性傳染病是氣候敏感性疾病。而歷史文獻來看，多次出現(xiàn)“夏四月”而后瘟疫流行，可見氣候異常是瘟疫發(fā)生的“前奏”。

積極治療和管理病患，能有效控制傳染源?！叭缮辏橥ㄌ煜退幰辕熋窦?。”還有一些特殊的手段處理傳染源：1207年，“余親見泰和六年丙寅，征南師旅大舉，至明年軍回，是歲瘴癘殺人，莫知其數(shù)。昏瞀懊依，十死八九，皆火之化也?！?，可見古人已經(jīng)使用火化處理尸體，1910年東北鼠疫流行期間，伍連德處理方法亦是如此。

隔離傳染源，限制傳染源的活動范圍，做好邊境檢疫工作，古人也做得非常好?！坝篮湍?，多疾疫。舊制，朝臣家有時疾，染有三人以上者，身雖無病，百日不得入宮。”看來古人很早就認識這個道理，還形成了制度沿用下來?！皶r疫疙瘩腫毒病者，古方書所論不見其說……自天眷、皇統(tǒng)間，生于嶺北，次于太原，后于燕薊。山野頗罹此患，至今不絕，互相傳染，多至死亡。”根據(jù)描述，此病可能為鼠疫（腺），因蒙古軍與西域戰(zhàn)爭而傳入。提示邊境檢疫工作的重要。

3.2 阻斷傳播徐徑，監(jiān)測水源、管好飲用水

限制人口聚集和流動、控制人畜共患病、水域水源監(jiān)測、戴口罩等均是控制傳播途徑的有效方法。

1228年，“春，大疫，比屋相枕藉，安吉尤甚，戶減十五六。烹魚者，率從腹中得人指發(fā)?！?315年，“會昌民素不知井飲，汲于河流，故多疾癘?！?520年，“瘟疫大行，得疾者，親友不相訪問，染之即不起?！比诵蠊不疾〉拿枋鲇校?681年，“晉寧疫，人牛多斃；曲陽大疫?！?698年，“秋，房縣人牛俱疫。”人牛共患病，可能是天花牛痘之類。一些傳染病有沿河流流向發(fā)展的趨勢，并且不是一個偶然現(xiàn)象，且鄰長江、黃河、淮河流域的省份疫情發(fā)生明顯比其他省份要高?，F(xiàn)在全民使用自來水，經(jīng)河流水域傳播傳染病已經(jīng)成為歷史。

3.3 做好易感人群的保護

我國開展的計劃免疫接種，避免了很多傳染的流行，很多疾病的癥狀及預后明顯得到控制?！抖徽疃ㄕ摗芬粫斗N痘論》篇中：“宋仁宗時，有丞相王旦，初生諸子俱苦于痘……”詳細講了宋仁宗時期種痘的方法?！坝致劮N痘法，起于明朝隆慶年間寧國府太平縣，姓氏失考，得之異人丹傳之家，由此蔓延天下?！?1657年，“春正月，痘疫時，民家延醫(yī)種痘。”比如“大頭瘟”，歷史上萬歷十年記載：“腫項，人見病及哭者即死。”“大頭瘟”即流行性腮腺炎，現(xiàn)代普遍進行了計劃免疫接種，病情輕，傳播范圍小。

可見，歷史上的瘟疫防控，有許多措施和方法值得我們借鑒，中醫(yī)學在預防接種、處理傳染源、阻斷傳染途徑上，均有較大貢獻，需要挖掘。隨著現(xiàn)代科技的進步，很多影響瘟疫形成和傳播的因素，已經(jīng)被有效的阻斷。計量史學在史料研究中的應用，有利于中醫(yī)學的發(fā)展，找到更多的文化自信，也有助于課程思政建設。

[1]霍俊江.計量史學研究入門[M].北京:北京大學出版社,2013.

[2]楊奕望.計量史學方法在中醫(yī)史學研究中的 探索[J].中華中醫(yī)藥雜志,2017,32(8):3381－3383.

[3]張志斌.中國古代疫病流行年表[M].福建科學技術出版社,2006.

[4]譚其驤.中國歷史地圖集[M].北京:中國地圖 出版社,1996.

[5]中國人民軍事博物館.中國戰(zhàn)爭史地圖集[M].北京:星球地圖出版社,2007.

[6]中國軍事史編寫組.中國歷代戰(zhàn)爭年表[M]. 北京:解放軍出版社,2002.

[7]王會安,聞黎明,等.中國地震歷史資料匯編[M].北京:科學出版社,1983.

[8]譚忠厚.基于GIS的中國南北分界的計算和 模擬[D].蘭州大學,2011.

[9]竺可楨.中國近五千年來氣候變遷的初步研究[J].中國科學,1973(2):168-189．

[10]劉禹,安芷生,Hans W,等.青藏高原中東部 過去2485年以來溫度變化的樹輪記錄[J].中國科學 (D輯:地球科學),2009,39(2):166-176．

[11]鄧云特.中國救荒史[M].北京:商務印書館, 2011.

[12]龔勝生,劉楊,張濤.先秦兩漢時期疫災地理 研究[J].中國歷史地理論叢,2010,35(3):96-112．

[13]譚曉東.歷史上傳染病疫情特征與處理[J]. 人民論壇,2020(10):53-57.

[14]吳雨.民國時期云貴川地區(qū)疫災流行與公共衛(wèi)生意識的變遷研究[D].華中師范大學,2019.

[15]齊曉鈺.民國時期京津冀地區(qū)疫災流行與 公共衛(wèi)生意識的變遷研究[D].華中師范大學,2019.

[16]石晴.民國時期云南常見傳染病流行與防治研究（1937-1949）[D].云南農(nóng)業(yè)大學,2017.

[17]廖揚.民國時期成都地區(qū)鼠疫防控研究（1912年 -1949年）[D].四川師范大學,2017.

[18]張琳.民國時期東北地區(qū)疫災流行與公共 衛(wèi)生意識的變遷研究[D].華中師范大學,2016.

[19]石國寧.民國時期兩廣地區(qū)疫災流行與公共衛(wèi)生意識的變遷研究[D].華中師范大學,2016.

[20]蘇東敏.民國時期閩浙地區(qū)疫災流行與公共衛(wèi)生意識的變遷研究[D].華中師范大學,2016.

[21]吳曉旭,田懷玉,周森,等.全球變化對人類傳染病發(fā)生與傳播的影響[J].中國科學:地球科學,2013,43 (11) :1743-1759.

[22]王穎,梁丁元,李垚,等.歷史傳染病疫情的環(huán)境與氣候特征初探及對新冠肺炎疫情的思考[J].環(huán)境 科學研究,2020,33(7) :1555-1561.

R183

A

1672-1047（2022）06-0111-06

10.3969/j.issn.1672-1047.2022.06.29

2022-11-02

周 瑤，男，內蒙古托縣人，副主任醫(yī)師。研究方向：醫(yī)學與教育。

[責任編輯：曾華]