基于信息擴散理論的邢臺市農(nóng)業(yè)旱災風險分析

張瑞清，李 敏，陳佳飛，商慧蓮，商彥蕊

（河北師范大學 資源與環(huán)境科學學院 河北省環(huán)境演變與生態(tài)建設(shè)重點實驗室，石家莊 050016）

基于信息擴散理論的邢臺市農(nóng)業(yè)旱災風險分析

張瑞清，李 敏，陳佳飛，商慧蓮，商彥蕊

（河北師范大學 資源與環(huán)境科學學院 河北省環(huán)境演變與生態(tài)建設(shè)重點實驗室，石家莊 050016）

以邢臺市2000－2007年的統(tǒng)計數(shù)據(jù)為基礎(chǔ)，應用基于信息擴散理論的模糊數(shù)學方法，建立邢臺市農(nóng)業(yè)旱災風險評價模型，對邢臺旱災風險進行分析。結(jié)果表明：邢臺旱災受災概率在受災指數(shù)為10%時達到峰值，在受災指數(shù)超過30%之后幾乎為零；成災風險值顯示邢臺市在應對旱災方面已取得顯著成效，但該市平均每1.7a遭受到一次受災面積超過10%的旱災，農(nóng)業(yè)旱災仍不容忽視。研究結(jié)果增加了對邢臺市旱情的清楚認識，對區(qū)域防災減災的可持續(xù)發(fā)展具有指導意義。

信息擴散；農(nóng)業(yè)旱災；風險分析；邢臺

近年來，在全球氣候變暖的背景下，干旱災害頻發(fā)，且強度不斷增強，氣象災害引發(fā)的極端事件日益增多，嚴重威脅著受影響地區(qū)的社會經(jīng)濟和農(nóng)業(yè)生產(chǎn)［1］。邢臺市干旱發(fā)生頻率高，波及面積廣，持續(xù)時間長，旱災次數(shù)和旱災面積逐年增加，進行農(nóng)業(yè)旱災研究意義重大。

本文基于信息擴散理論對邢臺市進行農(nóng)業(yè)旱災風險分析，采用針對小區(qū)域小樣本風險評估的模糊數(shù)學方法，得出該區(qū)旱災風險估計值，以期對區(qū)域災害風險研究，提供一個嘗試性案例。

1 基于信息擴散的農(nóng)業(yè)旱災模糊風險分析模型

風險分析的傳統(tǒng)數(shù)學方法是概率統(tǒng)計，即利用研究區(qū)歷史上不同風險水平災害的發(fā)生頻率，得到風險概率分布。但此法進行計算需30個以上樣本才能得到較為滿意的風險評估結(jié)果。實際研究中，面對數(shù)據(jù)缺乏，持有樣本量提供的認識風險信息并不完善，具有模糊不確定性［2］。為了彌補信息不足而考慮優(yōu)化利用樣本模糊信息，采用樣本集值化的模糊數(shù)學處理方法［3］，將一個單值樣本的資料信息擴散到指標論域中所有點，挖掘盡可能多的有用的信息，以提高系統(tǒng)風險識別精度的技術(shù)稱為模糊信息處理技術(shù)，最常用的模型是正態(tài)擴散模型［4－7］。其原理如下：

風險指數(shù)論域U＝｛u1，u2，u3，…，um｝，設(shè)x1，x2，x3，…，xn為風險因素指標在n年的實際測量值，它由可度量的值組成。稱X＝｛x1，x2，x3，…，xn｝為樣本集合。一個單值的樣本xi依式（1）可以將其所帶的信息擴散給指標論域U上的所有點。

其中h為擴散系數(shù)，可根據(jù)樣本中最大值b和最小值a及樣本點個數(shù)n來確定。

計算公式為：

可定義樣本值相對應的模糊子集的隸屬函數(shù)為：

故稱μxi（uj）為樣本點xi的歸一化信息分布。對μxi（uj）進行處理，就可以得到一種效果較好的風險評估結(jié)果。

其物理意義是：由觀測樣本集合X＝｛x1，x2，x3，…，xn｝，經(jīng)過信息擴散推斷出，如果災害指標觀測值只能取U＝｛u1，u2，u3，…，um｝中的一個，在將xi看作樣本代表時，觀測值為uj的樣本點個數(shù)為q（uj）。顯然，q（uj）通常不是一個正整數(shù)，但一定是一個不小于零的數(shù)。

Q事實上就是各uj點上樣本點數(shù)的總和，從理論上來講，必有Q＝n，但由于數(shù)值計算四舍五入的誤差，Q與n之間略有差別。

p（uj）就是樣本落在uj處頻率值（概率值），可以作為概率的估計值。對于某一指標X＝｛x1，x2，x3，…，xn｝，通?？蓪⒎治龅闹笜酥祒i取為論域U中某一個uj，那么其超越概率值應為：

P＝｛P1，P2，P3，…，Pm｝稱為旱災因子風險值。如果論域指標是災害損失值，則成為旱災損失風險值。

2 指標選取及數(shù)據(jù)來源

對災害風險進行定量估算，必須選取可行性指標。災害損失（貨幣度量）隨社會發(fā)展逐年增加，不符合數(shù)據(jù)的一致性要求，選取受災指數(shù)、成災指數(shù)等指標，相對而言有較好的一致性，可當作隨機變量進行頻率分析［8］。研究數(shù)據(jù)來自邢臺統(tǒng)計年鑒和河北農(nóng)村統(tǒng)計年鑒（1999－2008年），旱災災情及氣候數(shù)據(jù)見表1。本文采用旱災受災指數(shù)（受災率）、成災指數(shù)（成災率）作為風險指數(shù)進行旱災風險分析，定義：受災（成災）率＝旱災受災（成災）面積/總播種面積

表1 邢臺2000－2007年農(nóng)業(yè)旱災災情統(tǒng)計及年降水量

3 實例分析

3.1 邢臺市概況

邢臺市位于河北省南部，北緯36°50′－37°47′和東經(jīng)113°52′－115°49′之間，面積12 486km2，截止到2007年農(nóng)業(yè)人口5 180 178人，占市總?cè)丝跀?shù)的74.6%。年末耕地總面積占全市面積的52%，是河北省主要農(nóng)業(yè)地市之一。本區(qū)屬暖溫帶半濕潤大陸性季風氣候，春旱風大，夏熱多雨，秋涼時短，冬寒少雪。自西而東為山地、丘陵、平原階梯排列，地貌、土壤類型齊全，因水資源短缺、降水時間和空間上分布不均勻?qū)е滤临Y源組合條件較差，是典型的旱災區(qū)，農(nóng)業(yè)生產(chǎn)受旱災影響最為嚴重［9］。該市兼有華北平原和太行山區(qū)兩大地形區(qū)的典型特征，以邢臺市為例進行旱災風險研究，對華北平原和太行山區(qū)的農(nóng)業(yè)生產(chǎn)系統(tǒng)可持續(xù)發(fā)展有積極意義。

3.2 邢臺市農(nóng)業(yè)旱災特點

3.2.1 降水影響顯著 區(qū)域降水對農(nóng)作物的產(chǎn)量影響巨大，降水的不穩(wěn)定性和降水量的不足，使農(nóng)業(yè)干旱可以發(fā)生在作物生長的任何階段［10］。2000－2007年邢臺市全年平均降水量516mm，由于降水季節(jié)變率大，降水量、降水時間與農(nóng)作物的生長期需水量不匹配，造成作物生育期缺水嚴重。由表1看出，各年份降水量的高低與成災率的大小大體相一致，處于年平均降水量以上的年份成災率較低，反之較高。個別年份的差異跟人工灌溉程度等因素有關(guān)。

3.2.2 春旱多發(fā)突出 受所處氣候帶影響，春季多風少雨，最易形成旱災。據(jù)資料統(tǒng)計，受旱的年內(nèi)變化顯示，春旱最多，發(fā)生頻率達84%～94%，初夏旱頻率為25%～45%，伏旱頻率10%～30%，秋旱頻率為20%～33%［11］。春季是該區(qū)主要作物生長周期的關(guān)鍵期，春旱的高發(fā)嚴重威脅本區(qū)糧食產(chǎn)量，因此，春旱是造成經(jīng)濟損失的關(guān)鍵因素之一，也是邢臺旱情最突出的特點之一。

3.2.3 地域空間差異 大氣降水、農(nóng)田灌溉和土壤保水能力三大要素組合的地域分布不同，導致干旱發(fā)生頻率、受災面積存在明顯的地域差異。該市中低山區(qū)旱災面積較小，成災率為30%～50%；東部低平原為較旱區(qū)，成災率為20%～30%。山麓平原為少旱區(qū)，成災率在20%以下，且大旱較少（表2）。

3.3 邢臺市農(nóng)業(yè)旱災風險分析

設(shè)Xj（j＝1，2，…，8）分別為2000年，2001年，…，2007年邢臺市農(nóng)作物旱災受災面積與總播種面積的比值，即旱災受災率Xj。同理xj（j＝1，2，…，8）分別為2000年，2001年，…，2007年邢臺市農(nóng)作物旱災成災面積與總播種面積的比值，即旱災成災率xj（見表1）。旱災受災（成災）指數(shù)是衡量相對損失的指標，較客觀地反映出旱災的受災風險情況。由于樣本數(shù)量較少，故選用基于信息擴散理論的模糊風險分析模型，考慮計算精度需要，取受災（成災）指數(shù)論域U＝｛u1，u2，…，um｝＝｛0，0.01，…，1｝，將Xj和xj按上述公式（1）－（9）進行計算，可以計算出旱災受災（成災）風險估計值，并做出旱災受災指數(shù)的概率密度曲線及超越概率密度曲線如圖1和圖2所示。

圖1 邢臺市旱災受災率概率密度曲線

表2 邢臺市農(nóng)業(yè)旱災風險分析結(jié)果

由圖1得到任何受損區(qū)間上的概率值，從而推斷出不同災情水平發(fā)生的可能性大小。邢臺旱災受災率在10%的可能性最大，超過30%受災可能性的發(fā)生概率很小。圖2曲線有一個從凸到凹的趨勢，說明對應的概率密度曲線有一個極大值點，與圖1所反映的情況相符，該曲線坡度較陡，表明旱災發(fā)生的情況相對集中，確定在一定范圍內(nèi)，風險性較低。表2受災指數(shù)為10%的風險估計值為0.601 43，即≥10%的概率為0.601 43。就是說，在邢臺平均每1.7a就要遭受到一次受災面積超過10%的旱災。同理，受災面積超過20%的概率為0.094 13，即平均每10.6a就有一次受災面積為20%的旱災。這說明，邢臺出現(xiàn)受災面積到10%的災情，幾乎是每年一次，而達到20%以上受災面積的旱災十年一遇，邢臺旱災的發(fā)生頻率和強度對農(nóng)業(yè)的影響是比較嚴重的，是危害農(nóng)業(yè)生產(chǎn)的主要自然災害之一。將計算的成災風險值與受災風險值對照，風險性的概率要低許多，說明邢臺在應對旱災采取的抗旱措施已取得了顯著成效，進一步分析其對旱災適應性有更深遠的意義。將表3成災指數(shù)概率分布與圖3旱澇保收率進行對照，不同年份成災指數(shù)概率的大小與旱澇保守率的整體走勢比較一致。無論風險性的高低變化，圖3顯示旱澇保收率的整體趨勢是逐年提高，其影響因素綜合復雜，包括自然和人為因素的諸多方面，對其因素的分析還有待進一步探討。

表3 2000－2007年邢臺市成災指數(shù)概率分布

圖3 2000－2007年邢臺市旱澇保收率變化情況

4 結(jié)論

運用評估模型對邢臺市2000－2007年8a的統(tǒng)計數(shù)據(jù)進行風險評估，得出旱災影響的時間頻率和損失程度的概率分布，簡明直觀，便于對邢臺旱情有更清楚的認識。

邢臺旱災幾乎連年發(fā)生，受損程度主要集中在10%～20%，加之凍害、風災等其它災害的影響，對于以農(nóng)業(yè)為主要產(chǎn)業(yè)的邢臺市十分不利，農(nóng)業(yè)旱災不容忽視。要求有關(guān)部門，加強旱情監(jiān)測，抓好農(nóng)田水利設(shè)施建設(shè)，合理利用水資源，建立完善的農(nóng)業(yè)防災減災體系，將災害損失降到最低限度。

自然災害風險分析對數(shù)據(jù)的要求往往很高，尤其小區(qū)域的數(shù)據(jù)缺乏，無法使用傳統(tǒng)的概率統(tǒng)計方法，本文采用的數(shù)據(jù)量雖小，利用模型得出風險分析結(jié)果與實際情況相符，有其實用性，對邢臺農(nóng)業(yè)旱災的災害研究理論體系的完善有所幫助，對區(qū)域防災減災的可持續(xù)發(fā)展也有一定意義。

［1］何斌，武建軍，呂愛鋒.農(nóng)業(yè)干旱風險研究進展［J］.地理科學進展，2010，29（5）：557－559.

［2］黃崇福.自然災害風險評價理論與實踐［M］.北京：科學出版社，2005：180－183.

［3］梁保松，曹殿立.模糊數(shù)學及其應用［M］.北京：科學出版社，2007：15－38.

［4］王積全，李維德.基于信息擴散理論的干旱區(qū)農(nóng)業(yè)旱災風險分析：以甘肅省民勤縣為例［J］.中國沙漠，2007，27（5）：827－829.

［5］孫才志，張翔.基于信息擴散技術(shù)的遼寧省農(nóng)業(yè)旱災風險評價［J］.農(nóng)業(yè)系統(tǒng)科學與綜合研究，2008，24（4）：507－510.

［6］張俊香，李平日，黃光慶，等.基于信息擴散理論的中國沿海特大臺風暴潮災害風險分析［J］.熱帶地理，2007，27（1）：11－14.

［7］田敏，聶素芬，劉菲，等.昆明市東川區(qū)農(nóng)業(yè)旱災風險區(qū)識別研究［J］.災害學，2011，26（1）：41－42.

［8］盧峰本，黃瀅.基于信息擴散理論的水產(chǎn)養(yǎng)殖災害風險評價［J］.廣西氣象，2006，27（2）：22－23.

［9］商慧蓮.農(nóng)業(yè)旱災系統(tǒng)恢復及其驅(qū)動影響因素初探：以河北省邢臺市為例［D］.石家莊：河北師范大學，2010.

［10］商彥蕊，史培軍.自然災害系統(tǒng)脆弱性研究［M］.西安：西安地圖出版社，2004：124－126.

［11］李玉平，張義文.邢臺市土地生態(tài)環(huán)境保護與整治［J］.邢臺學院學報，2005，20（4）：98－99.

Risk Assessment of Agriculture Drought Disaster in Xingtai City Based on Information Diffusion Technique

ZHANG Rui－qing，LI Min，CHEN Jia－fei，SHANG Hui－lian，SHANG Yan－rui
（School of Geograph，Hebei Normal University，Hebei Key Laboratory of Environment Change and Ecological Construction，Shijiazhuang，Hebei050016，China）

By using the statistical data of Xingtai City in 2000－2007and the fuzzy mathematical method based on information diffusion principle，risk models of agricultural drought disaster of Xingtai City was established，the drought disaster risk of this city was analyzed.The results show that the damage probability of drought disaster in Xingtai City reaches maximum when the index is 10%and it is almost zero after the index is over 30%；the value of disaster risk shows that Xingtai City has made remarkable achievement in response to drought disaster，but this city suffered a drought that the drought－affected area is more than 10%on average 1.7a，agricultural drought still cannot be ignored.The results can provide the clear knowledge about drought situation of Xingtai City，and significant guidance for the sustainable development of regional disaster prevention and mitigation.

information diffusion；agriculture drought disaster；risk analysis；Xingtai City

P954

A

1005－3409（2011）06－0212－04

2011－06－07

2011－09－22

國家自然科學基金“區(qū)域農(nóng)業(yè)旱災適應性評價模型與風險防范模式研究”（41171402）；河北省高校重點學科建設(shè)項目

張瑞清（1985－ ），女，河北邯鄲人，碩士研究生，研究方向：資源評價與可持續(xù)利用。E－mail：zhangruiqing850729＠126.com

商彥蕊（1963－），女，河北蠡縣人，博士，教授，碩士生導師，主要研究方向：資源開發(fā)與減災。E－mail：shangyanrui＠126.com