變化環(huán)境對華南濕潤區(qū)洪水頻率變化的影響

1、變化環(huán)境對華南濕潤區(qū)變化環(huán)境對華南濕潤區(qū) 洪水頻率變化的影響分析洪水頻率變化的影響分析 陳曉宏陳曉宏1,2 1,2、汪麗娜 、汪麗娜3 3 1 1 中山大學(xué)水資源與環(huán)境研究中心中山大學(xué)水資源與環(huán)境研究中心 2 2 華南地區(qū)水循環(huán)與水安全廣東省普通高校重點實驗室華南地區(qū)水循環(huán)與水安全廣東省普通高校重點實驗室 3 3 華南師范大學(xué)地理科學(xué)學(xué)院華南師范大學(xué)地理科學(xué)學(xué)院 2012.72012.7 研究背景研究背景 全球氣候變化導(dǎo)致災(zāi)害性極端天氣事件的發(fā)生頻率和強度的增加，尤其是全球氣候變化導(dǎo)致災(zāi)害性極端天氣事件的發(fā)生頻率和強度的增加，尤其是 洪水災(zāi)害的頻繁發(fā)生。洪水災(zāi)害的頻繁發(fā)生。 以武江流域為例，進

2、入以武江流域為例，進入2121世紀以來：世紀以來： 2002 2002年出現(xiàn)了近年出現(xiàn)了近3030年一遇的洪水一次，同年出現(xiàn)近年一遇的洪水一次，同年出現(xiàn)近2020年一遇洪水一次；年一遇洪水一次； 20062006年出現(xiàn)了一場近年出現(xiàn)了一場近600600年一遇的洪水，同年又出現(xiàn)年一遇的洪水，同年又出現(xiàn)2020年一遇的洪水；年一遇的洪水；20082008年年 超超2020年一遇的洪水再次訪問該流域。年一遇的洪水再次訪問該流域。 LOGO 近近2020年來武江流域的小概率洪水年來武江流域的小概率洪水 事件時常出現(xiàn)，小概率洪水常態(tài)化了，事件時常出現(xiàn)，小概率洪水常態(tài)化了， 這與統(tǒng)計學(xué)、概率學(xué)上的概念不完

3、全符這與統(tǒng)計學(xué)、概率學(xué)上的概念不完全符 合，所得的結(jié)果難以令人信服。合，所得的結(jié)果難以令人信服。 研究背景研究背景 LOGO 研究思路研究思路 研究意義研究意義 研究思路和研究意義研究思路和研究意義 引發(fā)思考引發(fā)思考 初步初步 分析分析 洪水頻洪水頻 發(fā)現(xiàn)象發(fā)現(xiàn)象 解開矛盾解開矛盾 癥結(jié)癥結(jié) 氣候變氣候變 化化 內(nèi)驅(qū)動力內(nèi)驅(qū)動力外驅(qū)動力外驅(qū)動力 + 進一步進一步 思考思考 洪水特性發(fā)生洪水特性發(fā)生 變化變化 尋找適合當前環(huán)尋找適合當前環(huán) 境下的洪水頻率境下的洪水頻率 計算模型計算模型 v為流域防洪減災(zāi)、洪水資源化利用、水庫調(diào)度規(guī)則的建立、洪澇評估分析，為流域防洪減災(zāi)、洪水資源化利用、水庫調(diào)度規(guī)

4、則的建立、洪澇評估分析， 提供理論支撐。提供理論支撐。 挖掘洪水時間序列的特性。挖掘洪水時間序列的特性。 解讀氣候變化下洪水的異變性。解讀氣候變化下洪水的異變性。 研究氣候變化對洪水頻率理論及計算方法的影響，研究氣候變化對洪水頻率理論及計算方法的影響， 試圖尋找適合當前環(huán)境下的洪水頻率計算模型。試圖尋找適合當前環(huán)境下的洪水頻率計算模型。 研究目標 理論上的小概理論上的小概 率洪水率洪水 實際中的常實際中的常 遇洪水遇洪水 矛盾矛盾 現(xiàn)象現(xiàn)象 揭示隱藏在揭示隱藏在 背后的本質(zhì)背后的本質(zhì) 洪水特性的洪水特性的 挖掘挖掘 洪水異變性洪水異變性 的解析的解析 洪水演進過程時洪水演進過程時 - -頻域的

5、響應(yīng)頻域的響應(yīng) 洪水序列隨機性、洪水序列隨機性、 獨立性的辨識獨立性的辨識 P-III型分布使型分布使 用的限制用的限制 極值變量單一輸極值變量單一輸 入的不合理性入的不合理性 隨機性和獨隨機性和獨 立性的前提立性的前提 氣候變化氣候變化 影響影響 Gumbel-Logistic模型模型 洪峰流量洪峰流量 洪水總量洪水總量 洪峰水位洪峰水位 較真實洪水較真實洪水 頻率結(jié)果的頻率結(jié)果的 解答解答 技術(shù)路線技術(shù)路線 降雨量降雨量 研究內(nèi)容研究內(nèi)容 解析了武江流域洪水的動力學(xué)特性 綜合分析了武江流域洪水的異變性 解答了小概率洪水事件和常遇洪水之間的矛 盾癥結(jié)所在，并論證適合氣候變化背景下的 洪水頻率

6、計算模型 武江流域介紹 武江是北江流域第二大的一級支流，水系發(fā)達。武江干流在廣東境內(nèi)比降較 陡，平均比降1.27，流速大，洪水傳播時間快，是典型的彎曲型山區(qū)河流。 進入21世紀以后，全球氣候變化的作用，導(dǎo)致武江流域臺風雨增多。 珠江流域北江水系武江流域圖珠江流域北江水系武江流域圖 LOGO 1 洪水動力學(xué)特性解析 v 解析的方法 G_P算法、Cao氏算法計算嵌入維、Lyapunov指數(shù)的計算、 互信息算法確定延遲、重構(gòu)相空間和R/S分析方法。 6種方法與結(jié)果之間的相互連接: 恰當?shù)那‘數(shù)?+ 嵌入維數(shù)嵌入維數(shù)時間延遲時間延遲 奇異吸引子奇異吸引子重構(gòu)相空間重構(gòu)相空間 Lyapunov指指 數(shù)的

7、計算數(shù)的計算 初值的敏感性初值的敏感性 Cao氏算法氏算法互信息法互信息法 G_P算法算法 R/S分析方法分析方法 分形性分形性 無標度區(qū)的辨識無標度區(qū)的辨識 無標度區(qū)的準確性將直接影響無標度區(qū)的準確性將直接影響 分形維的準確性分形維的準確性 正確辨識正確辨識 洪水時間洪水時間 序列為純序列為純 隨機性？隨機性？ 混沌序列？混沌序列？ 1 洪水動力學(xué)特性解析 LOGO v 解析的時間序列的對象 由于洪水時間序列不同于其他水文時間序列，每年具 有一場或多場洪水過程，因此在本匯報中以兩種洪水時間 序列為研究對象： i)以1955-2007年為時間跨度，每年中選一個洪峰流 量值，由53個洪峰流量數(shù)據(jù)

8、組成一個洪峰流量時間序列; ii)以1955年、1956年、2006年和2007年的洪峰 流量過程為例，分析武江流域的洪水時間序列的動力特性。 v 解析結(jié)果 武江流域53場洪水和年洪峰流量具有混沌性、長程相 關(guān)性和分形性。 1 洪水動力學(xué)特性解析 LOGO v 氣候的變化導(dǎo)致水文時間序列的突變，洪水作為水 文時間序列的要素之一，即具有與之相同的共性， 又有自身的特性。 v 本研究結(jié)合洪水時間序列自身的特征，綜合分析各 項要素的異變性。 從以下幾方面分析洪水的異變性： 1）洪水強度大小的比較及異變性分析 2）洪水形態(tài)的分類及異變性辨識 3）Qmax、Qtotla和Q3d的異變性 2 武江流域洪水

9、異變的辨識 LOGO 2.1.1 強度指標的選擇 2.1.2 分類模型： 2 武江流域洪水異變的辨識 LOGO 2.1.3 分類結(jié)果及異變性的辨識： 洪水 等級洪水分類結(jié)果 IV2006 III1961,1968, 1994,2002 II 1955,1956,1957,1962,1969,1972,1973, 1975,1976,1977,1978,1980,1982, 1983,1985,1992,1993,1995,1997,1998, 1999,2001,2005,2007 I 1958,1959,1960,1963,1964,1965,1966, 1967,1970,1971,197

10、4,1979,1981, 1984,1986,1987,1988,1989,1990,1991, 1996,2000,2003,2004 2 武江流域洪水異變的辨識 LOGO 2.2 洪水形態(tài)的分類及異變性辨識洪水形態(tài)的分類及異變性辨識 2.2.1洪水形態(tài)指標選擇：洪水形態(tài)指標選擇： 1）偏度：偏度： 峭度：峭度：峰越尖銳，峭度越大；峰越平坦，峭度越小。峰越尖銳，峭度越大；峰越平坦，峭度越小。 2）時域特征（時域特征（6項指標）：項指標）： 3）頻域特征（頻域特征（5項指標）：項指標）：洪水時間序列的頻域帶個數(shù)、能量最為集中的頻域中洪水時間序列的頻域帶個數(shù)、能量最為集中的頻域中 心值、能量最為

11、集中的時間值、高頻率的頻域跨度值、高頻率的頻域時間持續(xù)值。心值、能量最為集中的時間值、高頻率的頻域跨度值、高頻率的頻域時間持續(xù)值。 2 武江流域洪水異變的辨識 LOGO 2 武江流域洪水異變的辨識 LOGO IIII型雙峰型型雙峰型 I I型雙峰型型雙峰型 2 武江流域洪水異變的辨識 LOGO 年洪峰流量 () 洪水總量 （ ） 最大72小時洪量 () 突變時間顯著水平突變時間顯著水平突變時間顯著水平 1994年*1993年*1993年* 2 武江流域洪水異變的辨識 LOGO 傳統(tǒng)的水文頻率分析（傳統(tǒng)的水文頻率分析（P-III型分布）僅考慮水文時間的型分布）僅考慮水文時間的 某一特征量，對洪水

12、特征的描述不全面，是導(dǎo)致采用某一特征量，對洪水特征的描述不全面，是導(dǎo)致采用P- III型分布計算出的洪水頻率結(jié)果與實際情況存在偏頗的型分布計算出的洪水頻率結(jié)果與實際情況存在偏頗的 原因之一。原因之一。 武江流域洪峰流量具有混沌性、分形性和長程相關(guān)性，武江流域洪峰流量具有混沌性、分形性和長程相關(guān)性， 因此現(xiàn)今洪水時間序列可能不完全服從傳統(tǒng)的因此現(xiàn)今洪水時間序列可能不完全服從傳統(tǒng)的P-III型分型分 布，亦是導(dǎo)致洪水頻率結(jié)果與實際情況存在偏頗的另一布，亦是導(dǎo)致洪水頻率結(jié)果與實際情況存在偏頗的另一 個原因。個原因。 將進一步探討小概率洪水事件與常遇洪水之間的矛盾癥將進一步探討小概率洪水事件與常遇洪水

13、之間的矛盾癥 結(jié)，并尋找適合氣候變化背景下的洪水頻率計算模型。結(jié)，并尋找適合氣候變化背景下的洪水頻率計算模型。 3 解析氣候變化下的洪水頻率值 LOGO 3.1 隨機性和獨立性的解析隨機性和獨立性的解析 v 時間序列的選取時間序列的選取 1 1）以以1955-2007年的年的53場洪水的洪峰流量時間序列為研究對象場洪水的洪峰流量時間序列為研究對象 2）以）以4場調(diào)查到的歷史洪水（場調(diào)查到的歷史洪水（1835年、年、1915年、年、1931年和年和1935年年 ） 與與1955-2007年的年的53場洪水作為輸入，解析場洪水作為輸入，解析57場洪水的洪峰流量的隨場洪水的洪峰流量的隨 機性和獨立性

14、。機性和獨立性。 v 解析方法解析方法 隨機性：通過自相關(guān)系數(shù)值的判斷隨機性：通過自相關(guān)系數(shù)值的判斷 ； 獨立性：輪次分析法。獨立性：輪次分析法。 3 解析氣候變化下的洪水頻率值 LOGO 3.1 隨機性和獨立性的解析隨機性和獨立性的解析 v 解析結(jié)果解析結(jié)果 1）53場洪水的洪峰流量具有隨機性，但由于年洪峰流場洪水的洪峰流量具有隨機性，但由于年洪峰流 量負輪次長的數(shù)學(xué)期望量負輪次長的數(shù)學(xué)期望E(53)=2.3642.0，表明它不具，表明它不具 有獨立性。有獨立性。 2）57場洪水的洪峰流量具有隨機性，但由于年洪峰流場洪水的洪峰流量具有隨機性，但由于年洪峰流 量負輪次長的數(shù)學(xué)期望量負輪次長的數(shù)

15、學(xué)期望E(57)=2.15382.0，表明該時間，表明該時間 序列亦不具有獨立性。序列亦不具有獨立性。 3 解析氣候變化下的洪水頻率值解析氣候變化下的洪水頻率值 LOGO 當采用當采用P-III型分布計算洪水頻率時，由于模型輸入的時間序列的非型分布計算洪水頻率時，由于模型輸入的時間序列的非 獨立性，導(dǎo)致獨立性，導(dǎo)致P-III型分布計算出的洪水頻率結(jié)果與實際存在偏差。型分布計算出的洪水頻率結(jié)果與實際存在偏差。 Gumbel認為現(xiàn)實狀態(tài)下的數(shù)據(jù)難以具有獨立性，并且實際研究表明認為現(xiàn)實狀態(tài)下的數(shù)據(jù)難以具有獨立性，并且實際研究表明 許多水文極值現(xiàn)象都服從許多水文極值現(xiàn)象都服從Gumbe1分布，由于分布

16、，由于Gumbel-logistic模型的模型的 限制條件較少，使用范圍更大，本研究中采用該模型計算洪水特征限制條件較少，使用范圍更大，本研究中采用該模型計算洪水特征 量的頻率量的頻率/重現(xiàn)期值。重現(xiàn)期值。 3 解析氣候變化下的洪水頻率值解析氣候變化下的洪水頻率值 3.2 氣候變化下的洪水頻率計算模型的選取氣候變化下的洪水頻率計算模型的選取 LOGO 3 解析氣候變化下的洪水頻率值解析氣候變化下的洪水頻率值 3.2 氣候變化下的洪水頻率計算模型的選取氣候變化下的洪水頻率計算模型的選取 LOGO 3 解析氣候變化下的洪水頻率值解析氣候變化下的洪水頻率值 3.2 氣候變化下的洪水頻率計算模型的選取

17、氣候變化下的洪水頻率計算模型的選取 項目項目結(jié)果結(jié)果柯爾摩哥洛夫檢驗的臨界值（柯爾摩哥洛夫檢驗的臨界值（Dn） 53場洪水場洪水 (n=53) 洪峰流量洪峰流量0.1224 D(53)=0.17981 洪水總量洪水總量0.1119 57場洪水場洪水 (n=57) 洪峰流量洪峰流量0.1236 D(57)=0.17231 洪峰水位洪峰水位0.0536 LOGO 3 解析氣候變化下的洪水頻率值解析氣候變化下的洪水頻率值 3.2 氣候變化下的洪水頻率計算模型的選取氣候變化下的洪水頻率計算模型的選取 項目項目年份年份 重現(xiàn)期（年）重現(xiàn)期（年） Gumbel-Logistic模型模型 53場洪水場洪水

18、1961年年 3.09 1968年年 10.56 1994年年 25.76 2002年年 10.94 2006年年 153 LOGO 3 解析氣候變化下的洪水頻率值解析氣候變化下的洪水頻率值 3.2 氣候變化下的洪水頻率計算模型的選取氣候變化下的洪水頻率計算模型的選取 項目項目 年份年份 重現(xiàn)期（年）重現(xiàn)期（年） P-III型分布型分布Gumbel-Logistic模型模型 57場洪水場洪水 1853年年 10037.23 1915年年 3017.52 1931年年 2012.86 1935年年 138.5 1961年年 107.27 1968年年 137.64 1994年年 3016.15

19、2002年年 2213.21 2006年年 599131.05 LOGO 小結(jié)：小結(jié)： 氣候變化的作用下，對于不具有獨立性的時間序列，如果采用氣候變化的作用下，對于不具有獨立性的時間序列，如果采用P- III型曲線分布的統(tǒng)計方法，不僅不能客觀描述洪水頻率的結(jié)果，型曲線分布的統(tǒng)計方法，不僅不能客觀描述洪水頻率的結(jié)果， 且計算出的重現(xiàn)期難以從數(shù)理統(tǒng)計學(xué)的角度理解。且計算出的重現(xiàn)期難以從數(shù)理統(tǒng)計學(xué)的角度理解。 因此，本研究從洪水具有的特性，以及氣候變化作用下展示出的因此，本研究從洪水具有的特性，以及氣候變化作用下展示出的 異變性出發(fā)，根據(jù)異變性出發(fā)，根據(jù)Gumbel-Logistic模型，結(jié)合兩變量

20、的優(yōu)勢，多模型，結(jié)合兩變量的優(yōu)勢，多 角度、較全面刻畫出洪水頻率在氣候變化作用下的真實結(jié)果。角度、較全面刻畫出洪水頻率在氣候變化作用下的真實結(jié)果。 3 解析氣候變化下的洪水頻率值解析氣候變化下的洪水頻率值 LOGO 4.1 分析對象的選擇分析對象的選擇 v 樣本的選取樣本的選取 1)以年最大抽樣法得到樣本，即以年最大抽樣法得到樣本，即AM系列。系列。 2)超定量法取樣：超定量法取樣： i）平均每年以）平均每年以1個超定量洪峰納入分析序列中獲得的樣本，即個超定量洪峰納入分析序列中獲得的樣本，即 POT1。 i i ）平均每年以）平均每年以3個超定量洪峰納入分析序列中獲得的樣本，即個超定量洪峰納入

21、分析序列中獲得的樣本，即 POT3。 v 時間序列的分段時間序列的分段 1)1951-1980年為一個時間段，年為一個時間段， 2)1981-2010年為另一個時間段。年為另一個時間段。 4 解析氣候變化下的洪水頻率值解析氣候變化下的洪水頻率值 LOGO 4.1 分析對象的選擇分析對象的選擇 v 獲得的序列：獲得的序列： 1）AM_1(1951-1980年年)、AM_2 (1981-2010年年); 2）POT1_1(1951-1980年年)、POT1_2 (1981-2010年年); 3）POT3_1(1951-1980年年) 、POT3_2 (1981-2010年年); v 模型的輸入：模

22、型的輸入： P-III型分布：型分布：以最大日徑流量或最大以最大日徑流量或最大3天流量天流量 Gumbel-Logistic模型：模型：以最大日徑流量和最大以最大日徑流量和最大3天流量天流量 4 解析氣候變化下的洪水頻率值解析氣候變化下的洪水頻率值 LOGO 4 珠江流域其他站點洪水頻率值解析珠江流域其他站點洪水頻率值解析 4.2 AM_1序列的頻率值 站名站名 K-S檢驗結(jié)果檢驗結(jié)果獨立性檢驗獨立性檢驗 隨機性檢驗隨機性檢驗 QdQ3dQdQ3dQdQ3d 遷江遷江0.0910.1161.52 三岔三岔0.1310.08321.75 南寧南寧0.0950.0912.3332 Qd LOGO

23、4 珠江流域其他站點洪水頻率值解析珠江流域其他站點洪水頻率值解析 具有獨立性的最大日徑流量三岔站具有獨立性的最大日徑流量三岔站: 結(jié)論：結(jié)論：P-III型分布計算結(jié)果與型分布計算結(jié)果與Gumbel-Logistic模型計算結(jié)果較為接近。模型計算結(jié)果較為接近。 4.2 AM_1序列的頻率值序列的頻率值 LOGO 4 珠江流域其他站點洪水頻率值解析珠江流域其他站點洪水頻率值解析 不具有獨立性的最大日徑流量遷江站和南寧站不具有獨立性的最大日徑流量遷江站和南寧站 ： 結(jié)論：結(jié)論：P-III型分布計算結(jié)果與型分布計算結(jié)果與Gumbel-Logistic模型計算結(jié)果差異較大。模型計算結(jié)果差異較大。 4.2

24、 AM_1序列的頻率值序列的頻率值 4 珠江流域其他站點洪水頻率值解析珠江流域其他站點洪水頻率值解析 具有獨立性的最大具有獨立性的最大3天徑流量遷江站和南寧站天徑流量遷江站和南寧站 ： 結(jié)論：結(jié)論：1)P-III型分布計算結(jié)果與型分布計算結(jié)果與Gumbel-Logistic模型計算結(jié)果相接近。模型計算結(jié)果相接近。 2)存在差異是由于存在差異是由于Gumbel-Logistic模型以二元輸入為背景進行計算的。模型以二元輸入為背景進行計算的。 4.2 AM_1序列的頻率值序列的頻率值 4 珠江流域其他站點洪水頻率值解析珠江流域其他站點洪水頻率值解析 4.3 AM_2序列的頻率值序列的頻率值 站名站

25、名 K-S檢驗結(jié)果檢驗結(jié)果獨立性檢驗獨立性檢驗 隨機性檢驗隨機性檢驗 最大日徑流最大日徑流 量量 最大最大3天流天流 量量 最大日徑流最大日徑流 量量 最大最大3天流天流 量量 最大日徑流最大日徑流 量量 最大最大3天流天流 量量 百色百色0.11680.08081.82 長壩長壩0.13060.070521.7 高道高道0.10060.139422 4 珠江流域其他站點洪水頻率值解析珠江流域其他站點洪水頻率值解析 具有獨立性的最大日徑流量高道站和長壩站具有獨立性的最大日徑流量高道站和長壩站 ： 結(jié)論：結(jié)論：P-III型分布計算結(jié)果與型分布計算結(jié)果與Gumbel-Logistic模型計算結(jié)果相

26、接近。模型計算結(jié)果相接近。 4.3 AM_2序列的頻率值序列的頻率值 4 珠江流域其他站點洪水頻率值解析珠江流域其他站點洪水頻率值解析 百色站百色站 不具有獨立性的最大日徑流量不具有獨立性的最大日徑流量(a)： 具有獨立性的最大具有獨立性的最大3天流量天流量(b)： 結(jié)論：結(jié)論：圖圖(a)中中P-III型分布計算結(jié)果與型分布計算結(jié)果與Gumbel-Logistic模型計算結(jié)果相差較模型計算結(jié)果相差較 大；圖大；圖(b)中的結(jié)果相接近。中的結(jié)果相接近。 4.3 AM_2序列的頻率值序列的頻率值 4 珠江流域其他站點洪水頻率值解析 4.4 POT1_1序列的頻率值 站名站名 K-S檢驗結(jié)果檢驗結(jié)果

27、獨立性檢驗獨立性檢驗 隨機性檢驗隨機性檢驗 最大日徑流最大日徑流 量量 最大最大3天流天流 量量 最大日徑流最大日徑流 量量 最大最大3天流天流 量量 最大日徑流最大日徑流 量量 最大最大3天流天流 量量 新和新和0.1670.10522 長壩長壩0.1850.1261.22 石角石角0.0880.0881.5562 4 珠江流域其他站點洪水頻率值解析珠江流域其他站點洪水頻率值解析 具有獨立性的最大日徑流量的新和站具有獨立性的最大日徑流量的新和站 ： 結(jié)論：結(jié)論：P-III型分布計算結(jié)果與型分布計算結(jié)果與Gumbel-Logistic模型計算結(jié)果相接近。模型計算結(jié)果相接近。 4.4 POT1_

28、1序列的頻率值序列的頻率值 4 珠江流域其他站點洪水頻率值解析珠江流域其他站點洪水頻率值解析 不具有獨立性的最大日徑流量的長壩站和石角站：不具有獨立性的最大日徑流量的長壩站和石角站： 結(jié)論：結(jié)論：圖中圖中P-III型分布計算結(jié)果不可靠。型分布計算結(jié)果不可靠。 4.4 POT1_1序列的頻率值序列的頻率值 4 珠江流域其他站點洪水頻率值解析珠江流域其他站點洪水頻率值解析 具有獨立性的最大具有獨立性的最大3天流量的長壩站和石角站：天流量的長壩站和石角站： 結(jié)論：結(jié)論：P-III型分布計算結(jié)果與型分布計算結(jié)果與Gumbel-Logistic模型計算結(jié)果相接近。模型計算結(jié)果相接近。 4.4 POT1_

29、1序列的頻率值序列的頻率值 4 珠江流域其他站點洪水頻率值解析珠江流域其他站點洪水頻率值解析 4.5 POT1_2序列的頻率值序列的頻率值 站名站名 K-S檢驗結(jié)果檢驗結(jié)果獨立性檢驗獨立性檢驗 隨機性檢驗隨機性檢驗 最大日徑流最大日徑流 量量 最大最大3天流天流 量量 最大日徑流最大日徑流 量量 最大最大3天流天流 量量 最大日徑流最大日徑流 量量 最大最大3天流天流 量量 百色百色0.25450.10842.172 南寧南寧0.13490.144122 4 珠江流域其他站點洪水頻率值解析珠江流域其他站點洪水頻率值解析 具有獨立性的最大日徑流量的南寧站：具有獨立性的最大日徑流量的南寧站： 結(jié)論

30、：結(jié)論：P-III型分布計算結(jié)果與型分布計算結(jié)果與Gumbel-Logistic模型計算結(jié)果相接近。模型計算結(jié)果相接近。 4.5 POT1_2序列的頻率值序列的頻率值 4 珠江流域其他站點洪水頻率值解析珠江流域其他站點洪水頻率值解析 百色站百色站 不具有獨立性的最大日徑流量不具有獨立性的最大日徑流量(a)： 具有獨立性的最大具有獨立性的最大3天流量天流量(b)： 結(jié)論：結(jié)論：圖圖(a)中中P-III型分布計算結(jié)果與型分布計算結(jié)果與Gumbel-Logistic模型計算結(jié)果相差較模型計算結(jié)果相差較 大；圖大；圖(b)中的結(jié)果相接近。中的結(jié)果相接近。 4.5 POT1_2序列的頻率值序列的頻率值

31、4 珠江流域其他站點洪水頻率值解析珠江流域其他站點洪水頻率值解析 4.6 POT3_1序列的頻率值序列的頻率值 站名站名 K-S檢驗結(jié)果檢驗結(jié)果獨立性檢驗獨立性檢驗 隨機性檢驗隨機性檢驗 最大日徑流最大日徑流 量量 最大最大3天流天流 量量 最大日徑流最大日徑流 量量 最大最大3天流天流 量量 最大日徑流最大日徑流 量量 最大最大3天流天流 量量 蔗香蔗香0.1010.09522.15 柳州柳州0.1090.0732.2112 高道高道0.1010.0842.4712 4 珠江流域其他站點洪水頻率值解析珠江流域其他站點洪水頻率值解析 具有獨立性的最大日徑流量的蔗香站：具有獨立性的最大日徑流量的

32、蔗香站： 結(jié)論：結(jié)論：P-III型分布計算結(jié)果與型分布計算結(jié)果與Gumbel-Logistic模型計算結(jié)果相接近。模型計算結(jié)果相接近。 4.6 POT3_1序列的頻率值序列的頻率值 4 珠江流域其他站點洪水頻率值解析珠江流域其他站點洪水頻率值解析 柳州站柳州站 不具有獨立性的最大日徑流量不具有獨立性的最大日徑流量(a)： 具有獨立性的最大具有獨立性的最大3天流量天流量(b)： 結(jié)論：結(jié)論：圖圖(a)中中P-III型分布計算結(jié)果與型分布計算結(jié)果與Gumbel-Logistic模型計算結(jié)果相差較模型計算結(jié)果相差較 大；圖大；圖(b)中的結(jié)果相接近。中的結(jié)果相接近。 4.6 POT3_1序列的頻率值

33、序列的頻率值 LOGO 4 珠江流域其他站點洪水頻率值解析珠江流域其他站點洪水頻率值解析 高道站 不具有獨立性的最大日徑流量(a)： 具有獨立性的最大3天流量(b)： 結(jié)論：圖(a)中P-III型分布計算結(jié)果與Gumbel-Logistic模型計算結(jié)果相差較大； 圖(b)中的結(jié)果相接近。 4.6 POT3_1序列的頻率值 LOGO 4 珠江流域其他站點洪水頻率值解析 4.7 POT3_2序列的頻率值 站名站名 K-S檢驗結(jié)果檢驗結(jié)果獨立性檢驗獨立性檢驗 隨機性檢驗隨機性檢驗 最大日徑流最大日徑流 量量 最大最大3天流天流 量量 最大日徑流最大日徑流 量量 最大最大3天流天流 量量 最大日徑流最

34、大日徑流 量量 最大最大3天流天流 量量 蔗香蔗香0.09980.069722.1 貴港貴港0.06290.06722.222.11 LOGO 4 珠江流域其他站點洪水頻率值解析 結(jié)論：P-III型分布計算結(jié)果與Gumbel-Logistic模型計算結(jié)果相接近。 4.7 POT3_2序列的頻率值 具有獨立性的最大日徑流量的蔗香站： LOGO 4.7 POT3_2序列的頻率值 4 珠江流域其他站點洪水頻率值解析 不具有獨立性的最大日徑流量的貴港站： P-III型分布計算的結(jié)果與Gumbel-Logistic模型計算的存在很大的差 異。其差異值主要體現(xiàn)在2005年流量為15700m3/s（在右圖中表現(xiàn)為離群點 據(jù)），P-III型分布認為該流量對應(yīng)的重現(xiàn)期為180.53年，而Gumbel- Logist