孟雪芹, 高春芳, 曹亞文

(安徽師范大學(xué) 數(shù)學(xué)與統(tǒng)計(jì)學(xué)院, 安徽 蕪湖 241000)

旱澇災(zāi)害屬于氣象災(zāi)害之一, 在所有的自然災(zāi)害中, 旱澇災(zāi)害發(fā)生頻率最高、影響范圍最廣、持續(xù)時(shí)間最長, 這不僅嚴(yán)重制約著我國社會經(jīng)濟(jì)的發(fā)展和穩(wěn)定, 也對人民的生命財(cái)產(chǎn)安全造成威脅.海陸因素使得我國有顯著的季風(fēng)氣候, 而旱澇災(zāi)害形成的直接原因就是季風(fēng)氣候?qū)е麓髿猸h(huán)境異常, 其中我國長江中下游地區(qū)和華北地區(qū)受旱澇災(zāi)害的影響最為嚴(yán)重[1].望江縣位于長江中下游, 因此也是旱澇災(zāi)害頻發(fā)的地區(qū).

為了預(yù)防和減少旱澇災(zāi)害帶來的危害, 國內(nèi)許多學(xué)者對旱澇災(zāi)害的致災(zāi)因素、發(fā)生規(guī)律及發(fā)展趨勢進(jìn)行了研究.1937年徐近之[2]利用統(tǒng)計(jì)分析法研究拉薩旱澇的現(xiàn)狀與成因, 這也是國內(nèi)對旱澇研究的起點(diǎn).1953年安振[3]提出興修水利保留水土解決旱澇問題, 得到了國內(nèi)專家對旱澇災(zāi)害研究的重視, 在該研究上取得了很大進(jìn)步.近80年來國內(nèi)對旱澇災(zāi)害問題的研究經(jīng)歷了三個(gè)階段: 以旱澇災(zāi)害現(xiàn)象描述分析為主的萌芽階段; 以研究旱澇災(zāi)害的分布規(guī)律及變化趨勢為主的發(fā)展階段; 以旱澇災(zāi)害的預(yù)測和防控舉措為主的成熟階段[4].

1 模型、變量與數(shù)據(jù)處理

1.1 旱澇主要致災(zāi)因素的選取

以各氣象因素?cái)?shù)據(jù)為依據(jù), 通過置信區(qū)間[5]與誤差分析法[6]得到每個(gè)觀測因素比重, 從而可選出導(dǎo)致旱澇災(zāi)害最主要的致災(zāi)因素.望江縣每年都有不同程度的旱澇災(zāi)害, 因此該地區(qū)每年的氣象數(shù)據(jù)對旱澇災(zāi)害的研究都具有一定價(jià)值.以氣壓、氣溫、水汽壓、相對濕度、降雨量為氣象觀測因素, 2006年至2017年望江縣境內(nèi)的氣象統(tǒng)計(jì)數(shù)據(jù)如圖1～5 所示.

圖1 氣壓數(shù)據(jù)

圖2 氣溫?cái)?shù)據(jù)

圖3 水汽壓數(shù)據(jù)

圖4 相對濕度數(shù)據(jù)

圖5 降雨量數(shù)據(jù)

對以上氣象數(shù)據(jù)資料進(jìn)行整理, 可得到多屬性決策矩陣, 記為X=(xij)144×5.為了保留更多的數(shù)據(jù)特征, 本文選用75%置信區(qū)間, 得到望江縣2006年至2017年每月各氣象因素?cái)?shù)據(jù)的置信區(qū)間, 如圖6～10 所示.

圖6 氣壓數(shù)據(jù)的置信區(qū)間

圖7 氣溫?cái)?shù)據(jù)的置信區(qū)間

圖8 水汽壓數(shù)據(jù)的置信區(qū)間

圖9 相對濕度數(shù)據(jù)的置信區(qū)間

對中間值矩陣X*與誤差估計(jì)值矩陣 ΔX*進(jìn)行規(guī)范化處理, 可得中間值規(guī)范化矩陣與誤差估計(jì)值規(guī)范化矩陣

中間值指標(biāo)權(quán)重與誤差估計(jì)值指標(biāo)權(quán)重計(jì)算公式分別為

然后利用客觀權(quán)重區(qū)間公式[7]

即可得到五個(gè)觀測因素的客觀權(quán)重區(qū)間分別為

下面采用0.1～0.9 標(biāo)度區(qū)間法對專家進(jìn)行問卷調(diào)查, 再對專家調(diào)查結(jié)果進(jìn)行整理可得區(qū)間數(shù)模糊互補(bǔ)判斷矩陣[8].利用誤差分析法計(jì)算得到五個(gè)觀測因素的主觀權(quán)重區(qū)間分別為

組合權(quán)重區(qū)間公式為

其中λ1,λ2分別為相應(yīng)的客、主觀的權(quán)重系數(shù).由于主客觀權(quán)重是公平競爭的平等關(guān)系[9], 故取λ1=λ2=0.5, 即可得五個(gè)觀測因素的組合權(quán)重區(qū)間為

對組合權(quán)重區(qū)間結(jié)果進(jìn)行分析可知, 氣溫與降雨量的權(quán)重總占比高達(dá)71.54%～87.71%, 所以氣溫和降雨量是望江縣旱澇災(zāi)害的主要致災(zāi)因素.

1.2 望江縣十二年間氣象數(shù)據(jù)每月的組合綜合評價(jià)區(qū)間

即可得到望江縣2006年至2017年氣象數(shù)據(jù)的客觀屬性綜合評價(jià)區(qū)間如圖11所示.

圖11 客觀綜合評價(jià)區(qū)間

再根據(jù)客觀屬性綜合評價(jià)區(qū)間的計(jì)算方法可求得望江縣2006年至2017年氣象數(shù)據(jù)的主觀屬性綜合評價(jià)區(qū)間如圖12所示.

圖12 主觀綜合評價(jià)區(qū)間

由主、客觀綜合評價(jià)區(qū)間的組合公式

其中λ1=λ2=0.5, 可得到望江縣2006年至2017年氣象數(shù)據(jù)的綜合評價(jià)區(qū)間如圖13所示.

圖13 綜合評價(jià)區(qū)間

1.3 望江縣十二年間氣象數(shù)據(jù)每月綜合評價(jià)區(qū)間的排序

根據(jù)上述得到的綜合評價(jià)區(qū)間, 用可能度的區(qū)間排序[11]對其進(jìn)行計(jì)算.利用可能度公式

得到排序向量, 結(jié)果如圖14所示.

圖14 氣象數(shù)據(jù)的序

根據(jù)序的計(jì)算公式可知: 序值越大, 干旱越嚴(yán)重; 序值越小, 洪澇越嚴(yán)重.由此, 從圖14中可得出以下結(jié)論: 望江縣2006年至2017年十二年間干旱一般發(fā)生在12月至2月之間, 其中普遍存在于1月, 而洪澇一般發(fā)生在6月至11月, 其中普遍存在于8月、9月和10月, 并且從2015年開始, 旱澇情況逐漸趨于一個(gè)范圍, 且干旱時(shí)間變長.

2 結(jié)論

根據(jù)國家防汛抗旱總指揮部辦公室文件——《干旱評估標(biāo)準(zhǔn)(試行)》[12]中干旱標(biāo)準(zhǔn)條件可知, 本文對望江縣旱澇災(zāi)害的評估與實(shí)際情況相符合.所以對望江縣2006年至2017年旱澇災(zāi)害的評估方案是合理準(zhǔn)確的, 將來可以提前做出一些防控措施.

本文是將望江縣氣象局的氣象數(shù)據(jù)信息作為望江縣綜合站點(diǎn)來研究的, 除此之外, 我們還可以對望江縣每個(gè)站點(diǎn)的氣象數(shù)據(jù)進(jìn)行細(xì)分, 進(jìn)一步分析每個(gè)站點(diǎn)區(qū)域的旱澇災(zāi)害情形, 得出的結(jié)果將會更加精確.