唐 莉，張永波，祝雪萍，王權(quán)威

(太原理工大學(xué)水利科學(xué)與工程學(xué)院，山西太原030024)

0 引 言

合理確定汛期分期不僅是分期設(shè)計(jì)洪水和分期汛限水位的基礎(chǔ)，對(duì)洪泛區(qū)管理、水庫防洪及實(shí)現(xiàn)洪水安全利用等也都具有重要意義[1]。汛期具有隨機(jī)性、確定性、模糊性等變化規(guī)律[2]，目前常用的汛期分期方法有[3- 5]：數(shù)理統(tǒng)計(jì)法、Fisher最優(yōu)分割法、模糊集分析法、變點(diǎn)分析法、分形分析法等。這些方法存在考慮因素單一，或計(jì)算繁瑣、理論依據(jù)薄弱、閾值選取主觀性強(qiáng)等缺陷，不能簡便合理地進(jìn)行汛期分期[6]。汛期分期是一個(gè)可變集，并且非汛期過渡到汛期再逐漸過渡到非汛期以及前汛期過渡到主汛期再到后汛期，符合唯物辯證法哲學(xué)規(guī)律。即，具有量變到質(zhì)變的過渡性[7]。本文采用以可變集理論為基礎(chǔ)，且數(shù)學(xué)概念清晰、理論性強(qiáng)的量變質(zhì)變判別模式，結(jié)合主成分分析法對(duì)李仙江流域進(jìn)行多類型指標(biāo)下的汛期分期，并將多種方法分期結(jié)果進(jìn)行了對(duì)比分析。

1 主成分分析法

汛期分期是基于多種指標(biāo)的有序聚類分析，各指標(biāo)對(duì)樣本的重要性有所不同，本文采用主成分分析法計(jì)算各指標(biāo)權(quán)重。具體步驟如下：

(1)構(gòu)造樣本陣X，對(duì)樣本陣X進(jìn)行標(biāo)準(zhǔn)化得到標(biāo)準(zhǔn)化陣

(1)

(2)計(jì)算標(biāo)準(zhǔn)化陣Y的相關(guān)系數(shù)矩陣

(2)

(3)計(jì)算特征值λi與特征向量ei。對(duì)于λi(i=1，2，…，m)由|λIm-R|=0求出，并對(duì)其進(jìn)行排序，即λ1≥λ2≥…≥λm≥0；同時(shí)根據(jù)每一個(gè)特征值λi求出其特征向量ei。

(4)計(jì)算主成分Fi貢獻(xiàn)率Ci及累計(jì)貢獻(xiàn)率Ai。即

(3)

(4)

當(dāng)λk對(duì)應(yīng)的Ci≥85%時(shí)，則其有k(k≤m)個(gè)與λ1，λ2，…，λk對(duì)應(yīng)的主成分。

(5)計(jì)算主成分載荷

通過式(5)得到lij，再由式(6)得到各主成分的線性模型

Fi=e1iX1+e2iX2+…emiXm(i=1，2，…，k)

(6)

(6)綜合得分模型為

其中，ai反映了各指標(biāo)變量在主成分中的綜合重要性。

(7)各指標(biāo)變量的權(quán)重為

(8)

在構(gòu)造好樣本陣X后，可用PASS軟件分析得到主成分貢獻(xiàn)率、累計(jì)貢獻(xiàn)率與主成分載荷(即成分荷載)。

2 量變質(zhì)變判別模式的汛期分期原理與步驟

2.1 對(duì)立差異函數(shù)

令

(9)

D:U→[-1，1]u|→D(u)∈[-1，1]

(10)

圖1 對(duì)立差異函數(shù)示意

2.2 模式原理與步驟

(1)構(gòu)建n個(gè)指標(biāo)特征值矩陣

其中，xim為時(shí)序m指標(biāo)i的特征值。

(2)確定時(shí)序研究對(duì)象ut指標(biāo)i的特征值xim落入汛前期與主汛期或主汛期與后汛期相對(duì)隸屬度為1的標(biāo)準(zhǔn)值區(qū)間矩陣[Ni1，Ni2]，已知矩陣X，可得標(biāo)準(zhǔn)區(qū)間矩陣Y=(Nih)，其中i=1，2，3…n;h=1，2則

Ni1=(minxim，maxxim)，Ni2=(maxxim，minxim)

(11)

式中，maxxim，minxim分別為時(shí)序集xim的最大、最小特征值。

(3)計(jì)算xim落入[Ni1，Ni2]區(qū)間的相對(duì)隸屬度

μA(μt)=(xim-Ni1)/(Ni2-Ni1)

(12)

(4)根據(jù)主成分分析法確定各指標(biāo)的權(quán)重wi。

(5)計(jì)算ut綜合相對(duì)隸屬度

(13)

式中，α為優(yōu)化準(zhǔn)則參數(shù)；P為距離參數(shù)。α=1為最小一乘方優(yōu)化準(zhǔn)則，α=2為最小二乘方優(yōu)化；考慮到α=2對(duì)距離具有放大或縮小效應(yīng)，在汛期分期應(yīng)用中選取α=1。若采用優(yōu)化準(zhǔn)則參數(shù)α=1，距離參數(shù)為海明距離，即P=1；則式(13)變?yōu)?/p>

(14)

若汛期分期為非線性系統(tǒng)，可采用歐氏距離，即P=2，式(13)變?yōu)?/p>

(15)

(16)

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

李仙江流域國境內(nèi)的流域面積為19 309 km2，屬云南省西南縱谷區(qū)南部，位于我國云貴高原西南邊緣、橫斷山脈南段，屬南亞熱帶高原季風(fēng)氣候，受地形、季風(fēng)、低緯度的影響，具有復(fù)雜多變的氣候特征。降水年內(nèi)分配極不均勻，流域的多年平均汛期為6～11月。

3.1 基本資料及指標(biāo)選取

選取汛期6月～11月為分期的研究論域，以旬為基本單位，將整個(gè)論域劃分為18個(gè)旬。以流域內(nèi)多年汛期逐日降雨資料為基本資料，取能反映水庫控制流域汛期內(nèi)暴雨變化特征的5個(gè)因子作為影響因子，即多年旬平均雨量、旬最大1 d雨量(mm)、旬最大3 d雨量(mm)、旬最大7 d雨量(mm)。洪水指標(biāo)選?。貉嗄暄骄髁?m3/s)[9]。為了對(duì)不同類型指標(biāo)對(duì)分期結(jié)果的影響進(jìn)行對(duì)比，本文分研究考慮綜合指標(biāo)的汛期分期、研究僅考慮降雨指標(biāo)的汛期分期和研究僅考慮洪水指標(biāo)的汛期分期3種情況分析。

3.2 綜合因素分析

根據(jù)汛期劃分指標(biāo)特征值，得到從汛前期過渡到主汛期和從主汛期過渡到后汛期的暴雨指標(biāo)指標(biāo)特征值矩陣X1(具體略)；根據(jù)矩陣X1和式(11)可得從汛前期過渡到主汛期和從主汛期過渡到后汛期暴雨指標(biāo)的兩個(gè)標(biāo)準(zhǔn)區(qū)間矩陣

用SPSS軟件對(duì)標(biāo)準(zhǔn)化的樣本數(shù)據(jù)進(jìn)行主成分分析，可得表1及表2。

表1 解釋的總方差

表2 成分1多年旬平均指標(biāo)矩陣

由表1可知，主成分1特征累積方差貢獻(xiàn)率≥85%。因此，第1個(gè)主成分能夠反映各指標(biāo)的大部分信息。

綜合得分模型F=0.470X1+0.472X2+0.474X3+0.473X4+0.328X5。

各指標(biāo)的權(quán)重ω=(ω1，ω2，ω3，ω4，ω5)=(0.212，0.213，0.214，0.213，0.148)。

應(yīng)用式(12)與矩陣Y11和Y12，可得從汛前期過渡到主汛期和從主汛期過渡到后汛期各降雨指標(biāo)的相對(duì)隸屬度矩陣U1t(略)；采用主成分分析法所得權(quán)重，根據(jù)式(13)可得到在線性模型和非線性模型下的相對(duì)優(yōu)屬度。

當(dāng)α=1，P=1時(shí)v1(u1t)=(0.21 0.03 0.73

0.84 0.96 0.93 0.84 0.81 0.70 0.64 0.52

0.44 0.46 0.29 0.39 0.15 0.05 0)；

當(dāng)α=1，P=2時(shí)v1(u1t)=(0.30 0.08 0.70

0.79 0.90 0.91 0.84 0.81 0.70 0.64 0.53

0.44 0.45 0.29 0.39 0.15 0.05 0)；

0.82 0.93 0.92 0.84 0.81 0.70 0.64 0.53

0.44 0.45 0.29 0.39 0.15 0.05 0)。

0.43 0.63 0.86 0.84 0.68 0.62 0.40 0.28

0.05 -0.11 -0.09 -0.42 -0.22 -0.69

-0.90 -1]。

應(yīng)用量變與質(zhì)變判別模式對(duì)基于綜合指標(biāo)的汛期演變過程分析如下：

3.3 降雨與洪水因素分析

0.85 0.64 0.59 0.37 0.30 0.08 -0.10

-0.32 -0.18 0.15 0.46 0.71 1 0.77 0.52

0.25 0.06 -0.04 -0.09 -0.12 -0.20

-0.41 -0.49 -59)。

3.4 結(jié)果比較和分析

基于不同指標(biāo)的對(duì)立差異度時(shí)程變化狀況比較見圖2。

圖2 基于不同指標(biāo)的對(duì)立差異度時(shí)程變化狀況

結(jié)合各項(xiàng)指標(biāo)對(duì)立差異度與圖2可知：

(1)依據(jù)綜合指標(biāo)。前汛期、主汛期與后汛期可分別定為6月1日到6月22日、6月23日到9月17日、9月18日到11月30日。

(2)依據(jù)降雨指標(biāo)。前汛期、主汛期與后汛期可分別定為6月1日到6月21日、6月21日到9月18日、9月19日到11月30日

(3)依據(jù)洪水指標(biāo)。前汛期、主汛期與后汛期可分別定為6月1日到6月30日、7月1日到9月20日、9月21日到11月30日。

表3給出基于暴雨因子、洪水因子、綜合因子以及基于數(shù)理統(tǒng)計(jì)法、模糊集合分析法、Fisher最優(yōu)分割法對(duì)研究流域汛期分期結(jié)果的比較。

圖2與表3顯示：

(1)基于洪水指標(biāo)的分期滯后于基于暴雨指標(biāo)的分期，而基于綜合指標(biāo)的分期基本介于兩者之間，9月中旬～10月上旬出現(xiàn)分歧。其原因可能是本文資料有限，所選洪水因素指標(biāo)單一，存在小幅誤差?？傮w上體現(xiàn)了洪水是暴雨與地表下墊面等多種因素綜合作用的反映，洪水要滯后于暴雨，這與實(shí)際情況是相符的，一定程度上反映了汛期分期的物理背景。

表3 基于不同指標(biāo)的汛期分析結(jié)果

(2)本文基于研究結(jié)果與數(shù)理統(tǒng)計(jì)法、模糊集合分析法、Fisher最優(yōu)分割法的比較說明了本文的分析成果可靠。

4 結(jié)論與建議

(1)汛期分期結(jié)果與分期指標(biāo)有密切關(guān)系，實(shí)際應(yīng)用的應(yīng)綜合考慮多種指標(biāo)，讓分期結(jié)果更客觀；另外，也可以利用不同指標(biāo)的分期結(jié)果，提供防洪安全或挖掘洪水資源潛力。

(2)本文采用的方法分期結(jié)果能反映指標(biāo)的物理意義，有較強(qiáng)的理論背景，數(shù)學(xué)概念清楚；計(jì)算方便，結(jié)果直觀，分期合理。

(3)我國南北方汛期時(shí)序具有明顯的差異性，北方采用Fisher最優(yōu)分割的結(jié)果多為三段式，然而南方采用Fisher最優(yōu)分割的結(jié)果可為四段式，本文采用的理論方法視汛期為三段式，具有一定局限性，本研究根據(jù)對(duì)立差異度的意義，提出基于量變質(zhì)變模型的過渡期劃分標(biāo)準(zhǔn)，雖與Fisher最優(yōu)分割法[9]結(jié)果表現(xiàn)出較高的一致性，但也需要在今后的研究中不斷地驗(yàn)證與完善。