潘建波，王文川，徐冬梅

（1.華北水利水電學(xué)院 土木與交通學(xué)院，河南 鄭州450011；2.華北水利水電學(xué)院 水利學(xué)院，河南 鄭州450011）

0 引言

中國是農(nóng)業(yè)大國，農(nóng)業(yè)旱災(zāi)是農(nóng)業(yè)生產(chǎn)的重大威脅.根據(jù)國家防汛抗旱總指揮部的統(tǒng)計(jì)，1979—2008年間，中國平均每年農(nóng)作物因旱受災(zāi)面積2.56×107h m2，成災(zāi)10.29×107h m2，因旱造成糧食損失244.6億kg，約占全年糧食總產(chǎn)量的5.7%［1］.特別是近幾年，中國旱情有加重的趨勢，從2008年11月下旬開始，中國冬麥區(qū)降水量明顯偏少，旱情持續(xù)發(fā)展，加之入冬后幾次強(qiáng)降溫天氣過程造成旱凍交加，對冬小麥安全越冬正常發(fā)育造成嚴(yán)重影響.這次干旱持續(xù)時間之長、受旱范圍之廣、受旱程度之重為歷史罕見.截至2009年2月3日，災(zāi)害已經(jīng)讓全國近43%的小麥產(chǎn)區(qū)受旱.河南、安徽、山東、河北、山西、陜西、甘肅等7個主產(chǎn)區(qū)小麥?zhǔn)芎?.53×106h m2，比上年同期增加8.93×106h m2.其中，嚴(yán)重受旱3.79×106h m2，比上年同期增加3.74×106h m2.2010年大旱導(dǎo)致西南五省區(qū)市除四川之外，至少有218萬人返貧，經(jīng)濟(jì)損失超350億元.冬麥區(qū)大部自2010年10月以來累計(jì)降水量不足25 mm，比多年同期偏少5至9成，截至2011年1月28日，全國作物受旱面積5.16×106h m2，257萬人、279萬頭大牲畜因旱出現(xiàn)飲水困難，僅山東、河南兩省受旱面積就達(dá)3.06×106h m2，占全國受旱面積的59%.因此，如何及時制定有效的抗旱減災(zāi)方案，合理地進(jìn)行灌區(qū)抗旱工作，這就需要對某一自然條件和工程條件下的灌區(qū)旱情程度進(jìn)行客觀的評價和分析.因此探討合理有效地旱情等級評價方法，確定其旱情等級，對管理部門抗旱救災(zāi)具有重要的指導(dǎo)意義.

目前，國內(nèi)外關(guān)于灌區(qū)干旱的研究多集中在干旱指標(biāo)的選取上［2］，而對評價的方法則較少涉及［3］.常用的灌區(qū)干旱評價方法多為單指標(biāo)評價法，如降雨頻率法、降雨隨機(jī)模擬法、干濕系數(shù)法、干燥度法、降雨距平法及干旱指數(shù)法等，多指標(biāo)的有模糊綜合評價法和可拓評價法［3－5］.模糊綜合評價是利用最大隸屬度原則確定分級，在模糊概念分級條件下，最大隸屬度原則易導(dǎo)致信息丟失而引起判別結(jié)果出現(xiàn)偏差［6］.作者基于對立統(tǒng)一定理，建立了灌區(qū)旱情等級評價的可變模糊評價新方法，模型能綜合考慮評價系統(tǒng)的線性特征和非線性特征，在具體應(yīng)用時，可根據(jù)問題的實(shí)際情況和問題本身的性質(zhì)進(jìn)行確定，進(jìn)一步提高評價結(jié)果的可靠性.

1 基于可變模糊集理論的對立統(tǒng)一定理

可變模糊集理論是陳守煜教授于2005年建立，并于2010年提出對立統(tǒng)一定理［6－7］，為了論述方便，本文作一簡介.

設(shè)論域U中的任意元素u的對立模糊概念（事物、現(xiàn)象）或u對立的基本模糊屬性，以A～與表示.在參考連續(xù)統(tǒng)區(qū)間［1，0］（對）與［0，的任一點(diǎn)上，對立模糊屬性的相對隸屬度分別為且

～

～在pl和pr之間必定存在中介點(diǎn)pm，如圖1所示，該點(diǎn)的對立模糊概念（事物、現(xiàn)象）或?qū)α⒒灸：龑傩缘南鄬﹄`屬度相等，即則pm為連續(xù)統(tǒng)區(qū)間［1，0］（對的漸變式質(zhì)變點(diǎn).因?yàn)閺膒l點(diǎn)變化到pr點(diǎn)，其間必定經(jīng)過0.5的中介點(diǎn)pm.由于所以Pm同時具有.在線段在線段pmpr間則中介 點(diǎn) pm必 為 連 續(xù) 統(tǒng) 區(qū) 間 ［1，0］（對的漸變式質(zhì)變點(diǎn)，即對立統(tǒng)一矛盾性質(zhì)的漸變式轉(zhuǎn)化點(diǎn).

圖1 對立統(tǒng)一定理的連續(xù)統(tǒng)示意圖Fig.1 The continuum sketch map of theorem of unity of opposites

2 基于對立統(tǒng)一定理的灌區(qū)農(nóng)業(yè)旱情等級評價新方法

2.1 計(jì)算待評灌區(qū)農(nóng)業(yè)對象u對旱情等級h的指標(biāo)相對隸屬度

設(shè)待評灌區(qū)對象u，根據(jù)已知的多個級別h（h＝1，2，…，c.為級別總數(shù)），多個指標(biāo)（i＝1，2，…，m，為評價指數(shù)總數(shù)）的指標(biāo)標(biāo)準(zhǔn)區(qū)間矩陣可表示為

式中：i＝1，2，…，m；h＝1，2，…，c.并設(shè)定1級干旱程度強(qiáng)于2級，2級干旱程度強(qiáng)于3級，c級最弱.

在實(shí)際干旱程度定量分析中指標(biāo)標(biāo)準(zhǔn)值區(qū)間［aih，bih］有兩種情況：①aih＜bih，稱遞增系列，其指標(biāo)特征值越大，干旱程度越強(qiáng)；②aih＞bih，稱遞減系列，其指標(biāo)特征值越小，干旱程度越強(qiáng).

根據(jù)實(shí)際情況的物理概念分析可知，在已知級別h指標(biāo)i特征值的相對隸屬度為1的矩陣：

的條件下，位于Mih與Mih＋（）1的指標(biāo)特征值ui，對級別h與（h＋1）構(gòu)成相對的對立模糊概念，根據(jù)對立統(tǒng)一定理有

式中：μh（ui）、μh＋1（ui）分別為待評對象u指標(biāo)i對級別h與（h＋1）的相對隸屬度.

設(shè)待評灌區(qū)對象u指標(biāo)i的特征值ui落入h與（h＋1）級相對隸屬度為1的M矩陣的區(qū)間［Mih，Mi（h＋1）］內(nèi)，則ui對h 級的相對隸屬度為

當(dāng)ui＝Mih時，μh（ui）＝1，由式（3）μh＋1（ui）＝0；當(dāng)ui＝Mi（h＋1）時，μh（ui）＝0，由式（4）μh＋1（ui）＝1；當(dāng)ui＝ （Mih＋Mi（h＋1））／2時，μh（ui）＝0.5，由式（4）μh＋1（ui）＝0.5.所以，應(yīng)用式（4）和式（5）可以計(jì)算待評對象u指標(biāo)i的特征值ui對級別h與h＋1的相對隸屬度.顯然，對于小于h，大于（h＋1）的指標(biāo)i的相對隸屬度均等于零.

2.2 計(jì)算待評灌區(qū)農(nóng)業(yè)對象u對旱情等級h的綜合相對隸屬度

根據(jù)式（5）和式（4）計(jì)算得到的相對隸屬度μh（ui）、μh＋1（ui）是指標(biāo)i對級別h 的單指標(biāo)相對隸屬度.灌區(qū)旱情等級評價是多指標(biāo)綜合評價問題，指標(biāo)具有不同的權(quán)重.設(shè)已知指標(biāo)權(quán)重向量為

式中：p為距離參數(shù)；p＝1為海明距離；p＝2為歐式距離.

若待評灌區(qū)對象u對級別h的多指標(biāo)綜合相對隸屬度以v（u）表示，則其算式為［7］

式中：α為優(yōu)化準(zhǔn)則參數(shù)；α＝1為最小一乘方優(yōu)化準(zhǔn)則；α＝2為最小二乘方優(yōu)化準(zhǔn)則.式（10）滿足：

（1）當(dāng)pi與pl重合時，μh（ui）＝1，dh（pl，pi）＝0，則vh（u）＝1，vh＋1（u）＝0；

（2）當(dāng)pi與pr重合時，μh＋1（ui）＝1，dh（pi，pr）＝0，則vh（u）＝0，vh＋1（u）＝1；

（3）當(dāng)pi與pm重合時，μh（ui）＝μh＋1（ui）＝0.5，dh（pl，pi）＝dh（pi，pr），則vh（u）＝vh＋1（u）＝0.5；

（4）當(dāng)pi位于區(qū)段plpm時，有dh（pl，pi）＜dh（pi，pr），則vh（u）＞vh＋1（u）；

（5）當(dāng)pi位于區(qū)段pmpr時，有dh（pl，pi）＞dh（pi，pr），則vh（u）＜vh＋1（u）；

因此，式（8）滿足對立統(tǒng)一定理，用于計(jì)算待評灌區(qū)對象u對級別h的相對隸屬度.

若采用優(yōu)化準(zhǔn)則參數(shù)α＝1，距離參數(shù)p＝1，則式（9）變?yōu)?/p>

式（9）是一個線性公式.如果考慮灌區(qū)綜合評價為非線性系統(tǒng)，可采用p＝2即歐氏距離，式（9）變?yōu)?/p>

這是一個非線性公式.因此綜合評價模型（9），只要改變距離參數(shù)p，既可用于線性系統(tǒng)評價，也可用于非線性系統(tǒng)評價.

當(dāng)α＝2時，無論p＝1還是p＝2，式（9）都是非線性公式.考慮到α＝2對距離比值dh（pl，pi）／dh（pi，pr）具有放大或縮小效應(yīng)，所以在綜合評價中宜采用α＝1.所以，可變模糊綜合評價模型式（9）在數(shù)學(xué)上滿足對立統(tǒng)一定理，說明它在哲學(xué)上符合唯物辯證法的對立統(tǒng)一規(guī)律.既可以評價線性的系統(tǒng)又可以評價非線性的系統(tǒng).這是該模型的特點(diǎn)和優(yōu)點(diǎn).

根據(jù)文獻(xiàn)［8］的級別特征值公式，可計(jì)算待評灌區(qū)的旱情級別特征值

3 應(yīng)用實(shí)例

關(guān)中西部四大灌區(qū)位于陜西省關(guān)中渭河平原西部，由馮家山灌區(qū)、寶雞峽灌區(qū)、羊毛灣灌區(qū)和石頭河灌區(qū)組成，涉及西安、咸陽、寶雞三個市，是陜西省主要的糧食、果林、油料、和蔬菜產(chǎn)地之一.關(guān)中地區(qū)屬于干旱和半干旱氣候帶，水資源匱乏，加上本區(qū)內(nèi)蓄水工程少，水資源調(diào)蓄能力差，供水保證率低，旱災(zāi)出現(xiàn)頻繁.探討合理有效地旱情等級評價方法，確定其旱情等級，對管理部門抗旱救災(zāi)具有重要的指導(dǎo)意義.為驗(yàn)證本文方法的可行性和有效性，引用文獻(xiàn)［3］中的數(shù)據(jù)，選取降雨指標(biāo)、徑流指標(biāo)、水庫蓄水指標(biāo)、地下水指標(biāo)為評價指標(biāo).劃為重度干旱、中度干旱、輕度干旱、不發(fā)生干旱4個等級表述灌區(qū)干旱程度［9］.灌區(qū)旱情等級劃分標(biāo)準(zhǔn)如表1所示，關(guān)中地區(qū)四大灌區(qū)評價因素指標(biāo)值如表2所示.

表1 灌區(qū)旱情等級劃分標(biāo)準(zhǔn)Tab.1 The division standard of drought grade of irrigation district

表2 關(guān)中地區(qū)四大灌區(qū)評價因素指標(biāo)值Tab.2 The index value of evaluation f actor of for irrigation district in Guanzhong area

根據(jù)表1、2可得關(guān)中地區(qū)四大灌區(qū)的旱情等級劃分標(biāo)準(zhǔn)區(qū)間矩陣和評價因素指標(biāo)特征值矩陣：

根據(jù)文獻(xiàn)［7］Mih與aih、bih的關(guān)系式可以為

根據(jù)式（12）可得級別h指標(biāo)i特征值的相對隸屬度等于1的矩陣：

根據(jù)評價因素指標(biāo)特征值矩陣X可知馮家山地下水指標(biāo)u4＝71，指標(biāo)（4）落入指標(biāo)相對隸屬度等于1的矩陣M的1、2級之間，即落入?yún)^(qū)間［M41，M42］＝［0.70，0.75］.應(yīng)用公式（5），i＝1，h＝1，得μ1（u4）＝（0.75－0.71）／（0.75－0.70）＝0.8，由對立統(tǒng)一定理公式（4）μ2（u4）＝0.2.根據(jù)對立統(tǒng)一定理，小于h大于（h＋1）級的指標(biāo)相對隸屬度為0，即μ3（u4）＝0、μ4（u4）＝0.則得馮家山地下水指標(biāo)u4的相對隸屬度向量為μ?（u4）＝（0.8，0.2，0，0）.類似可得到馮家山干旱指標(biāo)ui（i＝1，2，3，4，）對級別h（h＝1，2，3，4）的相對隸屬度矩陣為

為了與文獻(xiàn)［4］進(jìn)行對比分析，采用文獻(xiàn)［4］給出的4個灌區(qū)的評價指標(biāo)權(quán)重向量.

應(yīng)用可變模糊綜合模型公式（9），可求得馮家山對級別h（h＝1，2，3，4）的線性和非線性的綜合相對隸屬度，見表3，同理也可求出寶雞峽、石頭河和羊毛灣灌區(qū)的線性和非線性的綜合相對隸屬度，其結(jié)果如表3所示.根據(jù)級別特征值公式（12），可計(jì)算灌區(qū)旱情等級的線性和非線性模型的級別特征值見表3.

表3 線性和非線性模型的綜合相對隸屬度及級別特征值Tab.3 The comprehensive relative membership degree and grade characteristic value of linear and nonlinear model

4 結(jié)果分析與討論

根據(jù)表3中線性評價模型和非線性評價模型的級別特征值，可得出馮家山、寶雞峽、石頭河和羊毛灣灌區(qū)的綜合平均級別特征值分別為1.045、1、1.2、1.47.評價4個灌區(qū)均屬于1級，即重度干旱.但重度干旱的強(qiáng)弱有所不同，寶雞峽最強(qiáng)，馮家山次之，羊毛灣最弱.文獻(xiàn)［3］采用非對稱貼近度原則評價旱情等級，馮家山、寶雞峽、石頭河和羊毛灣灌區(qū)的等級變量特征值分別為3.6282，3.5985，3.5294，3.6237，評價4個灌區(qū)均屬于重度干旱與本文相同.但從給出的級別特征值來看，其重度干旱的強(qiáng)弱排序應(yīng)該為馮家山、羊毛灣、寶雞峽、石頭河，與本文不同.結(jié)合表1、2中數(shù)據(jù)進(jìn)行分析可知，馮家山有一項(xiàng)指標(biāo)地下水u4＝0.71屬于中度干旱，但非常接近重度干旱的指標(biāo)0.70，其它三項(xiàng)指標(biāo)屬于重度干旱；寶雞峽的4項(xiàng)干旱指標(biāo)特征值全部屬于重度干旱；石頭河有一項(xiàng)指標(biāo)水庫蓄水u3＝0.31屬于中度干旱，但接近中度干旱區(qū)間［0.20，0.50］中間值，其它三項(xiàng)指標(biāo)屬于重度干旱；羊毛灣有一項(xiàng)指標(biāo)地下水u4＝0.88屬于輕度干旱，其它三項(xiàng)指標(biāo)屬于重度干旱.因此，寶雞峽應(yīng)該是重度干旱的強(qiáng)度最大，羊毛灣重度干旱的強(qiáng)度最小，與本文量化分析的結(jié)果一致，說明本文的計(jì)算結(jié)果更符合實(shí)際情況.

5 結(jié)束語

作者基于可變模糊集理論的對立統(tǒng)一定理，建立了灌區(qū)旱情等級評價的可變模糊評價新方法，提出的方法滿足唯物辯證法的對立統(tǒng)一規(guī)律，理論嚴(yán)謹(jǐn)，概念清晰，計(jì)算簡便.模型綜合考慮灌區(qū)旱情評價系統(tǒng)的線性特征和非線性特征，均衡了兩者之間的差異，令結(jié)果更具說服力和可靠性.實(shí)例應(yīng)用結(jié)果表明，該方法不但能準(zhǔn)確給出旱情的等級，而且能清楚地給出其強(qiáng)弱排序，這對進(jìn)行抗旱的先后順序具有重要的指導(dǎo)意義.與文獻(xiàn)［4］和實(shí)際情況的對比分析，說明了本文的計(jì)算結(jié)果更符合實(shí)際情況.

［1］ 水利部辦公廳，水利部發(fā)展研究中心.水利改革發(fā)展30年回顧與展望［M］.北京：中國水利水電出版社.2010.

［2］ 王曉紅，胡鐵松，吳鳳燕，等.灌區(qū)農(nóng)業(yè)干旱評估指標(biāo)分析及應(yīng)用［J］.中國農(nóng)村水利水電，2003（7）：4－6.

［3］ 孫廷容，黃強(qiáng)，張洪波，等.基于粗集權(quán)重的改進(jìn)可拓評價方法在灌區(qū)干旱評價中的應(yīng)用［J］.農(nóng)業(yè)工程學(xué)報，2006，22（4）：70－74.

［4］ 劉曉寧，賈忠華.模糊綜合評價在關(guān)中灌區(qū)干旱評價中的應(yīng)用［J］.水資源與水工程學(xué)報，2005，16（2）：62－66.

［5］ 王斌，付強(qiáng)，王忠波，等.基于降水隨機(jī)模擬的旱作區(qū)農(nóng)業(yè)旱情等級劃分方法［J］.中國農(nóng)村水利水電，2011（3）：167－170.

［6］ 陳守煜.可變模糊集理論與模型及其應(yīng)用［M］.大連：大連理工大學(xué)出版社，2009.

［7］ 陳守煜.基于可變模糊集的辯證法三大規(guī)律數(shù)學(xué)定理及其應(yīng)用［J］.大連理工大學(xué)學(xué)報，2010，50（5）：838－844.

［8］ 陳守煜.工程模糊集理論與應(yīng)用［M］.北京：國防工業(yè)出版社，1998.

［9］ 成福云.關(guān)中灌區(qū)干旱指標(biāo)體系的建立與抗旱措施研究［D］.西安：西安理工大學(xué)水利水電學(xué)院，2003.