中國區域運氣學說適宜性研究(Ⅱ)
——寒熱氣候

王毅榮，伏芬琪，張尚祖，金榮花，陳祺，方峰，羅漢，龐朝云，丁瑞津

(1.中國氣象局蘭州干旱氣象研究所,甘肅省干旱氣候變化減災重點實驗室,中國氣象局干旱氣候變化與減災重點開放實驗室，甘肅 蘭州 730020；2.甘肅省人工影響天氣辦公室，甘肅 蘭州 730020；3.蘭州資源環境職業技術學院，甘肅 蘭州 730022；4.甘肅中醫藥大學，甘肅 蘭州 730019；5.國家氣象中心，北京 100081)

運氣學說最早記載于《黃帝內經》，創造性提出人體發病與自然氣候之間存在著密切關系。同時，運氣學說也是一套完整的醫學氣候預測模型，利用該模型進行簡單實用的推演，使得人們能夠預測年、季等不同時間尺度寒暑熱火等邪氣演變，進而對自然界的氣象變化和疾病的流行作出預測。如在《黃帝內經·氣交變大論篇》中，岐伯指出：“炎暑流行，金肺受邪；寒氣流行，邪害心火。”“民病瘧，少氣咳喘，血溢血泄注下，嗌燥耳聾，中熱肩背熱。”或“民病身熱煩心，躁悸，陰厥上下中寒、譫妄心痛，甚則腹大脛腫，喘咳，寢汗出憎風。”在近2 000多年的醫療實踐中，運氣學說展示出了極高的實用價值[1-3],近年來對運氣學說(理論)的研究大多集中于局地或應用[4-7]，加之我國各地氣候復雜多樣[8-10]，對我國全區時空代表性不強，有必要從更大時空尺度上加以研究及改進。在六淫邪氣中，寒熱邪氣是最為頻發的一種[11-12]，基于我國寒熱邪氣時空演變特征，對運氣學說寒熱氣候的適宜性予以研究，得到一些很有價值的結果，證實了運氣學說繼承和發揚的科學必要，體現中醫宏觀整體思維方法應對復雜系統的獨特優勢[13-14]。

1 材料與方法

1.1 材料

本文選取了分布在我國的，氣象記錄年代延續最長的83個氣象觀測站，收集并整理了55年(1951—2005年)間近地面大氣溫度逐日實測數據資料。

1.2 方法

為提取氣象要素場時空變化特征，采用自然正交分解函數(empirical orthogonal function,EOF)和旋轉自然正交分解函數(rotate empirical orthogonal function,REOF)方法對全國長期近地面溫度時空場特征進行分解刻畫，該方法有著十分明顯的優勢[15]；為分析寒熱過程的復雜性，計算大氣溫度序列分形維數D (Fractal Dimension)[16]，其大小反映了次區域寒熱變化過程線形態復雜程度。

以大寒日為年度起始分界日，對83站(氣象觀測站)55年的氣溫(剔去線性趨勢)年距平序列資料進行EOF和REOF分解。以氣溫的正負距平劃分熱寒，將正距平定義為熱年(偏熱)，負距平定義為寒年(偏冷)。

大運是主管每年全年氣候變化的歲運，又名中運、歲運，大運有太過與不及。每年的年號凡逢天干甲和己年，大運都屬于土運，從土運開始依次各年按天干五行順序排列。

將公元記年轉為甲子記年，如1951年記為辛卯年，依次按天干地支組合排列。歲運的計算依據《素問·天元紀大論》“甲己之歲土運統之，乙庚之歲金運統之，丙辛之歲水運統之，丁壬之歲木運統之，戊癸之歲火運統之。”“歲火太過炎暑流行，歲水太過寒氣流行；歲火不及寒乃大行，歲金不及炎火乃行。”(《素問·天氣交變大論》)和歲運太過與不及和五行生克關系，定義值年大運與寒熱氣候(見表1)，為便于計算，將歲運寒熱序列中熱年記為1，寒記為-1。1951—2005年間的寒熱序列構成了一個容量為55(55年)的取值為1或-1數字序列，簡稱為值年熱寒序列。

表1 值年大運與對應熱寒

2 結果

2.1 寒熱一致性

55年氣溫(已剔去線性趨勢)距平的第一空間型EOF1方差貢獻51%，明顯高于其它空間型，EOF1及其主分量(PC1)反映出全國氣溫場演變的最主要特征。年氣溫距平EOF1揭示寒熱時空演變最主要特征是，我國寒熱邪氣最容易在黃淮流域地區體現出來，該區域是我國心肺最易受邪的地區；PC1曲線反應了我國寒熱的年際趨勢變化，在1951—2005年間，1988年之前以寒邪為主，之后以熱邪為主；寒熱邪氣演變存在寒邪減弱、熱邪加重的趨勢，則心火受邪不斷減輕，肺金受邪趨向加重。

PC1與歲運寒熱序列間相關系數為0.272(在a=0.05水平上相關顯著)，表明了五運概括了我國寒熱演變的主要規律。大運年推算的寒熱年與實測一致則記為正確，則全國83站55年正確率為60.2%，89%的站點正確率在56%以上；若按氣象評分標準(預測與實況趨勢一致或相反趨勢距平不超過0.5 ℃)計，準確率達到86.6%。寒熱正確率存在地域差異，大致是華南地區和西北東部地區相對較高，西南等地相對較低。

五運寒熱與東部區(1區)相關性最高(表2)，我國中部區(2區)相對較低，西部區(3區)幾乎不相關；與六氣相關程度很低。

表2 RPC主分量與五運六氣熱寒序列相關系數

2.2 寒熱復雜性

在五運寒熱年基礎上，綜合六氣“子午之歲，上見少陰；丑未之歲，上見太陰；寅申之歲，上見少陽；卯酉之歲，上見陽明；辰戌之歲，上見太陽；巳亥之歲，上見厥陰。厥陰之上，風氣主之；少陰之上，熱氣主之；太陰之上，濕氣主之；少陽之上，相火主之；陽明之上，燥氣主之；太陽之上，寒氣主之。”之后，構成五運六氣寒熱年序列與氣溫PC1間相關系數反而下降，55年正確率為52.3%，正確率在56%以上站點降到20%。可見，寒熱演變就全國整體而言六氣不夠適宜，體現出寒熱演變的區域復雜性。

分維數D值大小反映出變化過程線形態復雜程度，主要反映的是相鄰峰谷差的平均變化。平均峰谷差越大，則過程線的起伏越大，過程線越復雜，表現在過程線局部峰谷變化頻繁、峰谷落差大。采用rk=2k(k=0,1,2,3)作為時間序列度量的特征線度計算，得到全國1951—2005年間候氣溫分維數D，總體上D值由東南向西北遞減，越過黃河河套向西地域分布變得凌亂，表明我國寒熱年際變化由東南向西北變化幅度加劇、寒熱行為愈趨復雜。

2.3 寒熱地域性

從年氣溫距平的EOF1概括了51%方差(1951—2005年)，也反映出我國氣溫變化還存在著地域性差異。因此，利用REOF分解將空間相關程度高的區域劃分在一起，前3個空間型統計學顯著，解釋了總方差的80%。依據空間模態分布中數值大小反映出區域內相關程度，大值區域內相關程度最高。第1空間模REOF1解釋了38.24%方差，數值大值區在我國東北地區，構成了一個高相關典型區，簡稱為東部區；第2空間模REOF2解釋了31.22%方差，大值區在我國中西部地區，構成了一個高相關典型區，簡稱中部區；第3空間模REOF3解釋了11.02%方差，大值區在我國最西部，構成了一個高相關區域，簡稱西部區。

2.4 寒熱邪氣演變

以沈陽、西安和烏魯木齊等分別代表1～3區，研究各次區域寒熱時間分布特征。從55年總體氣候態來看，沈陽、西安和烏魯木齊代表站的氣溫曲線形態十分相似，各地暑熱冬寒(冷)特征突出，4月初—10月底以熱邪為主(平均氣溫高于歷年平均)，其余時段寒邪為主。寒熱變幅烏魯木齊最大，西安變幅最小，可見東北區、中部區、西部區都是夏季肺金易受邪，冬季心火易受邪；其中烏魯木齊區域熱寒邪氣都較重，西安區域寒熱邪氣相對較輕，沈陽介于二者之間。

圖1給出1～3區REOF的時間系數(RPC)序列(含六階主旨函數)及其小波分析，六階主旨函數與原序列在a=0.01水平上相關顯著。

圖1 沈陽(a)、西安(b)和烏魯木齊(c)區域時間系數RPC序列(左)及其小波分析(右)

圖中可見，東部區(圖1a)寒熱階段交替出現、波動較大；振蕩在1960年代以15年周期為主，2000年代周期不明顯，其余時段以準5年周期為主；1988年之前以寒邪為主，1970年代最盛，心火易受邪；1988年之后熱邪為主，肺金易受邪。可見，寒熱邪氣演變的短周期、長周期和階段性與“五運”周期一致性程度高。

中部區(圖1b)大致寒熱2階段、波動較小；振蕩以2～3年為主，在1980年代以前15～16年長周期也比較明顯；1993年之前寒邪流行，心火易受邪；1993年之后熱邪為主，肺金易受邪。可見，寒熱邪氣演變的短周期、長周期和階段性與“五運”周期比較接近。

西部區(圖1c)寒熱邪氣交替頻繁、波動較大，振蕩1970年代以10年期為主，其余時段以準8年周期為主；1950年代寒邪最重，1970年代、1980年代中初期寒邪較重，期間心火易受邪；其余時段以熱邪為主，肺金易受邪。可見，寒熱邪氣在1970年代振蕩和階段性與“五運”周期接近，其余時段不一致。

由此可見，三個次區域寒熱演變時段內存在的準5年、10年和15年振蕩周期，與運氣學說中“五運”基本周期(5年)是相一致，同時寒熱也存在著周期不明顯或周期遠離5年的階段。

準2～3年和8年較之于“五運”基本周期(5年)存在著“過或不及”，需要進一步訂正；運氣學說中用“六氣從化、主客加臨、司天在泉”等思想予以修正。

給出1～3區REOF的時間系數(RPC)序列和修訂五運六氣逐年修正寒熱(熱為1，寒為-1)對應關系。對于東部區，將值年大運中陰火寒年修訂為熱年(木生火，歲火不及修訂為歲火太過)；對于中部區，土生金、土克水，故將值年大運中陰金、陰水年熱修訂為寒；對于西部區，本地金使弱金加強，金生水使弱水加強，故將值年大運中陰金和陰水年熱修訂為寒年；由此，相關程度顯著水平再次提高，中部地區提高最為明顯，成為比較適宜(表3)。五運六氣與各區相關程度明顯差異表明了，五運理論東部區最適宜，運氣學說成書與流行于中原地區有一定的道理。

表3 RPC主分量與修正后五運、六氣熱寒序列相關系數

即按六氣從化、主客加臨、司天在泉等予以修正，五行觀點認為五方主氣不同，如西方燥金、東方風木、南方火熱、北方寒水、中央濕土等，而五行之間又存在生克制化關系，結合逢陽的干年為太過，逢陰的干年為不及。太過即主歲的歲運旺盛而有余，為本氣流行；不及即主歲的歲運衰少而不足，克己之氣流行；《素問·天氣交變大論》“歲火太過炎暑流行，歲水太過寒氣流行；歲火不及寒乃大行，歲金不及炎火乃行。”和“子午之歲，上見少陰；丑未之歲，上見太陰；寅申之歲，上見少陽；卯酉之歲，上見陽明；辰戌之歲，上見太陽；巳亥之歲，上見厥陰。厥陰之上，風氣主之；少陰之上，熱氣主之；太陰之上，濕氣主之；少陽之上，相火主之；陽明之上，燥氣主之；太陽之上，寒氣主之。”及“少陰司天，熱淫所勝，怫熱至，火行其政……少陽司天，火淫所勝，則溫氣流行。”(《素問·至真要大論》)訂正之后，顯著水平再次提高，相關程度提高到顯著水平，通過a=0.05水平信度檢驗，達到適宜水平。

3 結論和討論

3.1 地域差異

從55年間我國寒熱演變存在東部區、中部區和西部區等3類高相關典型區，我國寒熱邪氣演變由東南向西變化幅度加劇、或寒或熱行為愈趨復雜；變幅西部區最大，中部區變幅最小，東北區介于二者之間。

3.2 周期特征

東部區寒熱階段交替出現、波動較大；振蕩在1960年代以15年周期為主，2000年代周期不明顯，其余時段以準5年周期為主；1988年之前以寒邪為主，1970年代最盛，心火易受邪；1988年之后熱邪為主，肺金易受邪；寒熱邪氣演變的短周期、長周期和階段性與“五運”周期一致性程度高。

中部區振蕩以2～3年為主，在1980年代以前15～16年長周期也比較明顯；1993年之前寒邪流行，心火易受邪；1993年之后熱邪為主，肺金易受邪；寒熱邪氣演變的短周期、長周期和階段性與“五運”周期比較接近。

西部區振蕩1970年代以10年期為主，其余時段以準8年周期為主；1950年代寒邪最重，1970年代、1980年代中初寒邪較重，期間心火易受邪；其余時段以熱邪為主，肺金易受邪；寒熱邪氣在1970年代振蕩和階段性與“五運”周期接近，其余時段不一致。

三個次區域寒熱演變時段內存在的準5年、10年和15年振蕩周期，與運氣學說中“五運”基本周期(5年)是相一致，同是寒熱也存在著周期不明顯或周期遠離5年的階段。

3.3 “運氣”修正

五運推演寒熱準確率較高，全國83站55年正確率為60.2%，89%的站點正確率在56%以上；單純六氣推演在我國大部地方不夠適合。準2～3年和8年較之于“五運”基本周期(5年)存在著“過或不及”，利用“六氣從化、主客加臨、司天在泉”等思想予以修正后效果大幅改善，都成為適宜區，其中全國東部區最適宜，中部區次之，西部最差。

值年大運學說中的五年周期與我國熱寒氣候演變準5年振蕩周期十分一致，是運氣學說概括了我國熱寒氣候演變最主要、最基本規律的內在根源；寒熱氣候演變中的準2～3年和8年較之于“五運”基本周期(5年)存在著“過或不及”，是“六氣從化、主客加臨、司天在泉”等思想存在的必然要求。予以修正后效果大為改善，在我國的地域適宜性大幅提高，但存在著難以量化和缺少氣候趨勢分析的不足，需要大力發掘加以提高。