基于分形理論的湘江流域洪水分期研究

2010-05-08 12:59:36石月珍鄭仰奇

水土保持通報 2010年5期

石月珍，李淼，鄭仰奇

(1.長沙理工大學水利工程學院，湖南長沙410114;2.湖南省水沙科學與水災害防治重點實驗室，湖南長沙410114)

湘江流域地處濕潤多雨地區，降雨充沛。然而，由于降水時空分布極不均勻，70%的降雨量集中在汛期，大量的水資源在春夏以洪水的形式體現，而且年際豐枯變化很大，容易形成非澇即旱，旱澇交替發生的局面。毛德華等[1]的研究指出，1949—1998年間，湘江流域中下游區平均1.7 a發生一次洪澇災害。而1998年、1999年和2003年秋冬季曾遇特大干旱，湘江長沙段創歷史最低水位，自來水廠取水頻頻告急，2009年10月6日湘潭水文站最低水位25.06 m，更是創歷史新低。湘江流域水資源供需矛盾十分尖銳，嚴重阻礙了長株潭城市群兩型社會配套改革試驗區的建設。很多學者和專家提出了洪水資源化的管理理念并進行研究，取得了不少成果[2-4]。對于湘江流域，在保證防洪安全的前提下，利用流域蓄、滯洪工程對汛期洪水實行分期調度將是解決水資源供需矛盾，提高洪水資源利用效率的重要舉措。因此，對湘江流域汛期洪水進行合理分期就顯得尤為重要。

洪水具有季節性變化規律，在洪水資源利用過程中，必須考慮這種特性。目前洪水分期的劃分一般靠經驗方法確定，將洪峰量級大體一致的時段劃分為一個分期，這種方法帶有主觀性。洪水流量隨季節、時間變化的過程是自然界中的一種復雜現象，在這種復雜現象的背后隱藏著特定的規律性，這種規律性在通常尺度不變規則下不易直觀觀測到，在分形理論特定尺度變換規則下，其規律性則顯而易見。且洪水點據系列具有隨機性、非線性、確定性和相似性，洪水點據系列的這些特性與分形理論研究的對象一致，因此可應用分形理論進行洪水分期的研究。

目前國內外在分形研究中普遍的做法是在對研究對象自相似性的定性認識基礎上，根據lnNN(ε)—ln(ε)關系曲線中是否存在直線段即無標度區來判斷是否為分形的[5-9]。本文的研究即是在上述宗旨下進行的。

1 分形維數的計算

樣本集具有標度不變性是以分形維數來定量的，即分形維數是一個常數。在一個洪水分期中如果分形維數Dc是相近的，那么就稱為同一洪水分期[10]。一般選用容量維進行分維計算[11]，即對于一個時間序列且對任意ε，在尺度小于ε的不規則集下，觀察ε→0時其度量 N(ε)的變化。分形維數Dc可用下式表示:

顯然當ε≈0時，下式成立:

點繪lnN(ε)—ln(ε)圖形，如果曲線中間直線段存在，則可認為此序列為分形，直線段的斜率為Dc，即計算出分形維數。

Smally等[12]在計算新赫布里底群島地震譜系列的時間分維時，引入了一個新的量NN，即相對度量:

式中:N(ε)——絕對度量;NT——劃分的總時段數，NT≈T/ε，T——研究時段總長;ε——時段步長。作出 lnNN(ε)—ln(ε)關系圖，若直線段的斜率為b，則容量維Dm為:

式中:d——拓撲維，洪水流量散點分布在Q—t二維平面上，故d取2，則容量維Dm=2-b。

2 湘江流域洪水分期

2.1 傳統方法確定洪水分期

湘潭水文站是國家重要水文站和湘江干流控制站，集水面積 81 638 km2。本文以湘潭水文站1997—2007年汛期實測日流量資料(4月1日至9月10日)作為分析對象，取歷年汛期日最大流量作為研究系列樣本。

如前所述，傳統洪水分期的劃分一般靠經驗確定，即將洪水量級大體一致的時段劃分為一個分期。具體做法為:以湘潭水文站1997—2007年歷年汛期實測日最大流量數據為基礎，繪制以時間為橫軸，流量值為縱軸的流量過程散點圖，如圖1所示。

圖1 湘潭水文站歷年汛期實測日最大流量散點

從圖2—4可以看出，lnNN(ε)—ln(ε)曲線中有直線段存在，可初步斷定洪水點據序列為分形。此外，還可看出，汛期各分期 lnNN(ε)—ln(ε)關系曲線簇中，A，B曲線的直線段斜率大體一致，而C曲線中直線段的斜率與A，B曲線的直線段斜率相差較大。因此，可認為B曲線對應的時期作為一個分期。據此，本文將湘潭水文站汛期洪水分期劃分為3段。表1列出了湘潭水文站汛期洪水分期結果。按照傳統經驗統計分析，由圖1可以看出，洪水分期分界點大概在5月10日，7月31日前后。如果洪水期劃分為3個洪水分期，那么第一個分期從4月1至5月9日，稱為前汛期;第二分期從5月10至7月31日，稱為主汛期;第三分期從8月1日至9月10日，稱為后汛期。