彭高輝 宋艷紅 王銀鴻 馬建琴

摘要：根據(jù)河南省17個(gè)氣象站1961-2013年逐月平均氣溫和降水資料，基于關(guān)聯(lián)規(guī)則，采用加權(quán)溫度干旱均一化指標(biāo)，分析河南省旱澇空間分布特點(diǎn)。結(jié)果表明：豫中、豫東、豫西、豫南、豫北地區(qū)各城市間旱澇關(guān)聯(lián)性較強(qiáng);空間分布上，干旱具有明顯的地域性且分布不均，從西北向東南方向逐漸減弱，主要發(fā)生在黃河流域及其以北地區(qū)，而洪澇呈團(tuán)狀分布，主要發(fā)生在黃河流域以南及偏北地區(qū)。

關(guān)鍵詞：旱澇特征;關(guān)聯(lián)規(guī)則;降水溫度均一化指標(biāo);河南省

中圖分類號(hào)：P338⁺.6 文獻(xiàn)標(biāo)志碼：A

河南省位于我國中東部、黃河中下游，屬于亞熱帶濕潤性季風(fēng)氣候與溫帶半濕潤性季風(fēng)氣候的過渡地帶，氣候復(fù)雜多樣、降水時(shí)空分布不均[1]，自然災(zāi)害頻繁發(fā)生，干旱和洪澇是主要?dú)庀鬄?zāi)害，往往造成嚴(yán)重的經(jīng)濟(jì)損失，其中對(duì)農(nóng)業(yè)影響尤為嚴(yán)重[2-3]。據(jù)1450-1949年的史料統(tǒng)計(jì)[4]，500a間發(fā)生全省性的干旱將近200a，其中大旱37a、特大旱6a。1949-2003年，旱災(zāi)面積平均每年為107.29萬hm²，澇災(zāi)面積平均每年為86.29萬hm²，幾乎每年都有旱澇災(zāi)害發(fā)生[5]。2003-2010年，干旱造成的年均糧食受災(zāi)面積為12.33萬hm²，洪澇造成的受災(zāi)面積為1.15萬hm²[6]。近幾年旱澇災(zāi)害仍持續(xù)發(fā)生，僅2016年7月，河南省洪澇災(zāi)害造成的農(nóng)作物受災(zāi)面積就達(dá)14.9萬hm²，折合成人民幣約為13.9億元[7]。因此，研究河南省旱澇時(shí)空分布，掌握旱澇整體分布規(guī)律，提前做好管理和規(guī)劃，具有重要的現(xiàn)實(shí)意義。

目前，對(duì)旱澇時(shí)空分布特征的研究大多采用氣象干旱指標(biāo)方法，如標(biāo)準(zhǔn)降水指數(shù)法[8]、降水距平百分率法[9]、Z指標(biāo)法[10]等，這些方法均基于降水要素，未考慮地理空間的關(guān)聯(lián)性。本文先對(duì)研究區(qū)域進(jìn)行劃分，然后利用數(shù)據(jù)挖掘中的關(guān)聯(lián)規(guī)則理論，構(gòu)造降水量和平均氣溫綜合指標(biāo)，通過對(duì)區(qū)域范圍內(nèi)的歷史災(zāi)害信息進(jìn)行挖掘，尋找災(zāi)害事件間的關(guān)聯(lián)性，進(jìn)而判斷災(zāi)害發(fā)生的類型和空間分布特征[11]。

1 研究方法

依據(jù)中國氣象數(shù)據(jù)網(wǎng)（http;//data.cma.cn/site/in-dex.html）提供的1961-2013年河南省17個(gè)氣象站點(diǎn)（見圖1）的月平均氣溫和月降水量數(shù)據(jù)，運(yùn)用關(guān)聯(lián)規(guī)則和加權(quán)溫度干旱均一化指標(biāo)[12]方法對(duì)河南省旱澇時(shí)空分布進(jìn)行研究。

1.1 降水溫度均一化指標(biāo)

降水溫度均一化指標(biāo)是降水標(biāo)準(zhǔn)化變量與溫度標(biāo)準(zhǔn)化變量之差：式中：I_s為降水溫度均一化指標(biāo);R為年平均降水量;R為多年平均降水量;σ_R為降水量均方差;σ_T為氣溫均方差;T為年平均氣溫;T為多年平均氣溫。

張強(qiáng)等[13]對(duì)3種干旱指標(biāo)進(jìn)行比較，考慮溫度對(duì)旱澇影響的比重，將標(biāo)準(zhǔn)化降水指標(biāo)與降水溫度均一化指標(biāo)相結(jié)合，提出了新的干旱指標(biāo)：式中：Z為標(biāo)準(zhǔn)化降水指數(shù)（見文獻(xiàn)[12]）;α_i為系數(shù)（依據(jù)降水與氣溫成顯著負(fù)相關(guān)的特點(diǎn)[14]求得具體值）;i為河南地區(qū)的劃分序號(hào)。

河南地區(qū)劃分及經(jīng)計(jì)算驗(yàn)證得到的系數(shù)α_i見表1，降水溫度均一化指標(biāo)的旱澇等級(jí)見文獻(xiàn)[12]，整理后的旱澇等級(jí)見表2。

1.2 基于地區(qū)劃分的關(guān)聯(lián)規(guī)則

關(guān)聯(lián)規(guī)則是數(shù)據(jù)挖掘中的典型研究方法，最早由Agxawal等[15]于1993年提出，其核心概念是支持度和置信度，定義如下。

設(shè)I={i₁，i₂，…，i_n}是所有項(xiàng)目的集合，其中的元素稱為項(xiàng);D為事務(wù)T的集合，事務(wù)T為項(xiàng)的集合，且T互1。每個(gè)事務(wù)有唯一的標(biāo)識(shí)TID。設(shè)I₁是I中1個(gè)項(xiàng)的集合，若I₁■T，則稱事務(wù)T包含I₁。

關(guān)聯(lián)規(guī)則是形如I₁→I₂的蘊(yùn)含式，其中I₁，I₂■I，且I₁∩I₂=φ。規(guī)則I₁→I₂在事務(wù)數(shù)據(jù)庫D中的支持度（support）是事務(wù)數(shù)據(jù)庫中包含I₁和I₂事務(wù)數(shù)與所有事務(wù)數(shù)之比，記為support（I₁→I₂）：

support（I₁→I₂）=|{T：I₁∪I₂■T，T∈D}|/|D|（3）

I₁→I₂在事務(wù)集中的置信度（confidence）是包含I₁和I₂的事務(wù)數(shù)與包含I₁的事務(wù)數(shù)之比，記為confidence（I₁→I₂）：

在集合D中找出具有用戶給定的滿足min support和min confidence的關(guān)聯(lián)規(guī)則，稱為強(qiáng)關(guān)聯(lián)規(guī)則。關(guān)聯(lián)規(guī)則挖掘主要有2個(gè)階段：①從集合D中找出大于等于用戶指定最小支持度的項(xiàng)目集;②利用頻繁項(xiàng)目集生成所需要的關(guān)聯(lián)規(guī)則。

1.3 空間分析方法

對(duì)于空間分布特征采用克里金插值法進(jìn)行分析。克里金插值法又稱空間局部插值法，以空間自相關(guān)為基礎(chǔ)，利用樣本數(shù)量和變異函數(shù)的結(jié)構(gòu)性，對(duì)區(qū)域化變量進(jìn)行無偏估計(jì)、最優(yōu)估計(jì)[16]。其表達(dá)式為

Z（s）=μ（s）+ε（s）（5）式中：Z（s）為待空間化的氣象值：μ（s）為已確定站點(diǎn)的氣象數(shù)據(jù)值;ε（s）為自相關(guān)隨機(jī)差。

2 結(jié)果與分析

2.1 旱澇關(guān)聯(lián)性判斷

2.1.1 干旱關(guān)聯(lián)性判斷

結(jié)合關(guān)聯(lián)性存在的理論基礎(chǔ)及區(qū)域干旱發(fā)生的規(guī)律，同時(shí)依據(jù)數(shù)據(jù)挖掘理論中的關(guān)聯(lián)規(guī)則算法[17]，并按式（3）和式（4）分別計(jì)算其支持度與置信度，結(jié)果見表3～表7，其中L_k為頻繁項(xiàng)目集，X_m-1為頻繁項(xiàng)目子集。

在最小支持度為7.5%和最小置信度為57.1%的閾值條件下，可知表3～表7均為強(qiáng)關(guān)聯(lián)規(guī)則，表明豫中、豫東、豫西、豫南、豫北地區(qū)各城市間干旱聯(lián)動(dòng)性較強(qiáng)門

2.1.2 洪澇關(guān)聯(lián)性判斷

同樣，對(duì)各區(qū)域城市進(jìn)行洪澇關(guān)聯(lián)規(guī)則挖掘，統(tǒng)計(jì)結(jié)果見表8～表12。

在最小支持度為5.7%和最小置信度為41.7%的閾值條件下，可知表8～表12均為強(qiáng)關(guān)聯(lián)規(guī)則，表明豫中、豫東、豫西、豫南、豫北地區(qū)各城市間洪澇聯(lián)動(dòng)性較強(qiáng)。

2.2 旱澇災(zāi)害的空間關(guān)聯(lián)性分析

2.2.1 干旱災(zāi)害的空間關(guān)聯(lián)性分析

圖2為1961-2013年河南省干旱頻率空間分布情況。可以看出，河南省空間上的干旱化趨勢(shì)為自豫西至豫東逐漸減弱，其中豫西地區(qū)的三門峽市和盧氏縣干旱化最為嚴(yán)重。究其原因，豫西黃土丘陵區(qū)屬大陸季風(fēng)性氣候區(qū)，具有秋季晴朗日照長、冬季寒冷雨雪少、春季干旱風(fēng)沙多、夏季氣溫高、地溫回升快、蒸發(fā)量大、空氣干燥的特點(diǎn)[18]。針對(duì)該區(qū)氣候特點(diǎn)，文獻(xiàn)提出推廣節(jié)水抗旱技術(shù)可大幅提高水資源利用效率，有效減少旱災(zāi)損失。

2.2.2 洪澇災(zāi)害的空間關(guān)聯(lián)性分析

從1961-2013年河南省洪澇頻率空間分布情況（見圖3）可以看出，洪澇頻率呈明顯的團(tuán)狀分布，且極其不均勻，較為嚴(yán)重的洪澇主要發(fā)生在信陽和孟津。究其原因，信陽市地勢(shì)南高北低，由西南向東北傾斜，偏北地區(qū)的異常季風(fēng)使得雨帶停滯在該區(qū)，降雨量偏多，從而形成洪澇災(zāi)害[20-21]。孟津縣地處豫西丘陵地區(qū)，地形復(fù)雜，屬亞熱帶和溫帶的過渡地帶，季風(fēng)環(huán)流影響明顯，夏季降雨量變化大，易發(fā)生洪澇。針對(duì)該區(qū)特點(diǎn)，文獻(xiàn)[22]提出采取合理的工程措施與非工程措施減輕洪澇災(zāi)害。

3 結(jié)論

（1）河南省干旱化地域差異明顯，趨勢(shì)為自豫西至豫東逐漸減弱，其中豫西地區(qū)的三門峽市和盧氏縣干旱化最為嚴(yán)重。

（2）河南地區(qū)的洪澇呈團(tuán)狀分布，較為嚴(yán)重的洪澇主要發(fā)生在信陽和孟津。

（3）豫中、豫東、豫西、豫南、豫北地區(qū)各城市間旱澇關(guān)聯(lián)性較強(qiáng)。

參考文獻(xiàn)：

[1]趙路偉，徐剛.河南省1961-2014年氣溫和降水量的時(shí)空變化特征[J].南水北調(diào)與水利科技，2016，14（3）：17-23.

[2]羅浩，包為民，王友恒.河南省旱澇態(tài)勢(shì)時(shí)空特征分析[J].人民黃河，2016，38（11）：39-43.

[3]張延偉，史本林，朱孔來.1960—2010年河南省極端氣候事件變化趨勢(shì)[J].人民黃河，2016，38（8）：10-13.

[4]郝秀平，張振偉，馬建琴，等.基于標(biāo)準(zhǔn)降水指數(shù)的河南省干旱時(shí)空演變規(guī)律分析[J].水電能源科學(xué)，2013（9）：4-7.

[5]朱業(yè)玉，程炳巖，王記芳.河南旱澇災(zāi)害的演變特征分析[J].災(zāi)害學(xué)，2006，21（3）：93-97.

[6]葛笑瑩，張靜雯，朱麗.1981—2011年河南省旱澇災(zāi)害的變化特征分析[J].現(xiàn)代農(nóng)業(yè)科技，2015（20）：194-195.

[7]大河網(wǎng).河南災(zāi)情造成經(jīng)濟(jì)損失48億，家庭損失13.7億元[EB/OL].（2016-07-25）[2017-09-10].http：//hn.ehinaso.com/finance/detail/20160725/1000200032826601469407512297827428_1.html.

[8]KUMAR R，MUSUUZA J L，VAN LOON A F，et al. Multi-scale Evaluation of the Standardized Precipitation Index as aGroundwater Drought Indicator[J]. Hydrology & EarthSystem Sciences Discussions，2016，12（8） ：7405-7436.

[9]CHENG Z. Establishment of Statistical Model for PrecipitationPrediction in the Flood Season in China[J]. Meteorological &Environmental Research，2011 （3） ： 1-3.

[10]WANG Y，ZHAO F，YAO Y，et al. Analysis on DroughtCharacteristics of Shiyang River Basin Based on Z Index[J].Journal of Catastrophology，2013，28（2）：100-106.

[11]李漢巨，梁浩波.基于關(guān)聯(lián)規(guī)則的自然災(zāi)害預(yù)測(cè)系統(tǒng)[J].計(jì)算機(jī)系統(tǒng)應(yīng)用，2017，26（7）：50-55.

[12]王殿武.區(qū)域旱澇特征分析及災(zāi)害預(yù)測(cè)技術(shù)研究[M].北京：中國水利水電出版社，2011：82-88.

[13]張強(qiáng)，鞠笑生，李淑華.三種干旱指標(biāo)的比較和新指標(biāo)的確定[J].氣象科技，1998（2）：48-52.

[14]吳洪寶.我國東南部夏季干旱指數(shù)研究[J].應(yīng)用氣象學(xué)報(bào)，2000，11（2）：137-144.

[15]AGRAWAL R，IMIELINSKI T，SWAMI A N. Mining Asso-ciation Rules Between Sets of Items in Large Databases[C]//Proceedings of the 1993 ACM-SIGMOD InternationalConference on Management of Data （SIGMOD' 93）. NewYork：ACM Press，1993：207-216.

[16]BEARMAN N，F(xiàn)ISHER P. Data for Paper： Using Sound toRepresent Spatial Data in ArcGIS[J]. Computers & Geo-sciences，2017，46（3）：157-163.

[17]彭高輝，馬建琴，賀淵博.基于可公度理論的關(guān)聯(lián)規(guī)則算法在干旱中的預(yù)測(cè)研究[J].南水北調(diào)與水利科技，2013，11（2）：151-153.

[18]張繼珍.類型學(xué)在豫西鄉(xiāng)村聚落更新與發(fā)展中的應(yīng)用研究[D].長沙：湖南大學(xué)，2010：10-11.

[19]楊栓群，高九思.豫西丘陵地區(qū)農(nóng)業(yè)節(jié)水抗旱技術(shù)研究[J].農(nóng)業(yè)災(zāi)害研究，2012，2（1）：74-75.

[20]陳軼.豫南木本園林植物季相特征及植物造景初探：以信陽市羅山縣植物配置為例[D].鄭州：河南農(nóng)業(yè)大學(xué)，2009：10-12.

[21]孔海江，吳勝安.淮河流域和長江中下游流域旱澇年的劃分[J].人民黃河，2012，34（3）：30-33.

[22]李文田，余國忠.信陽市洪澇災(zāi)害成災(zāi)機(jī)制與減災(zāi)對(duì)策研究[J].地域研究與開發(fā)，2013，32（5）：109-113.