• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    基于二次重現(xiàn)期的多變量洪水風(fēng)險(xiǎn)評(píng)估*

    2015-06-15 16:30:41陳子燊
    湖泊科學(xué) 2015年2期
    關(guān)鍵詞:洪量洪峰歷時(shí)

    黃 強(qiáng),陳子燊

    (中山大學(xué)水資源與環(huán)境系,廣州 510275)

    基于二次重現(xiàn)期的多變量洪水風(fēng)險(xiǎn)評(píng)估*

    黃 強(qiáng),陳子燊**

    (中山大學(xué)水資源與環(huán)境系,廣州 510275)

    由于洪水是一種具有多個(gè)特征屬性的隨機(jī)事件,頻率分析成為洪水風(fēng)險(xiǎn)評(píng)估的一種有效手段,多變量重現(xiàn)期與設(shè)計(jì)值的定義與計(jì)算則是洪水頻率分析中的重點(diǎn)和難點(diǎn).本文通過構(gòu)造洪水歷時(shí)、洪峰與洪量的聯(lián)合分布,介紹了一種新的多變量重現(xiàn)期定義——二次重現(xiàn)期,并探討了“或”重現(xiàn)期、“且”重現(xiàn)期和二次重現(xiàn)期對(duì)安全與危險(xiǎn)域識(shí)別的差異性,以及在洪水風(fēng)險(xiǎn)管理與工程設(shè)計(jì)中的合理性與可靠性.傳統(tǒng)的“或”和“且”多變量重現(xiàn)期對(duì)安全與危險(xiǎn)域的識(shí)別存在局限性,利用Kendall函數(shù)定義的二次重現(xiàn)期則提供了更加合理的安全與風(fēng)險(xiǎn)域識(shí)別,避免了對(duì)安全事件與危險(xiǎn)事件的錯(cuò)誤判定,更有利于指導(dǎo)洪水風(fēng)險(xiǎn)的管理.在給定的二次重現(xiàn)期條件下,依據(jù)出現(xiàn)概率最大原則推算的歷時(shí)、洪峰與洪量設(shè)計(jì)值組合可以滿足工程設(shè)計(jì)以較低成本承受較大風(fēng)險(xiǎn)的追求,相比于單變量設(shè)計(jì)值,考慮了洪水多個(gè)屬性聯(lián)合特征的多變量設(shè)計(jì)值提供了更加全面和可靠的參考信息.

    多變量洪水特征;極值分布;安全與危險(xiǎn)域;Kendall函數(shù);二次重現(xiàn)期;多變量設(shè)計(jì)值

    洪水作為一種自然災(zāi)害事件,每年造成的經(jīng)濟(jì)損失數(shù)以億計(jì),如何控制洪水并減少損失一直是災(zāi)害風(fēng)險(xiǎn)管理的重點(diǎn).由于洪水事件發(fā)生的隨機(jī)性,頻率分析成為洪水風(fēng)險(xiǎn)評(píng)估的一種有效手段,重現(xiàn)期與設(shè)計(jì)值的定義與計(jì)算則是頻率分析中的重點(diǎn).同一個(gè)洪水事件前后兩次發(fā)生的平均間隔時(shí)間即為該事件的重現(xiàn)期,若把該洪水事件造成的破壞和影響(以下稱效應(yīng))作為洪水風(fēng)險(xiǎn)管理與工程設(shè)計(jì)的標(biāo)準(zhǔn),重現(xiàn)期更大的洪水事件則可視為該標(biāo)準(zhǔn)下的危險(xiǎn)事件.傳統(tǒng)的單變量頻率分析并不能滿足洪水事件多個(gè)特征屬性的特點(diǎn),隨著Copula函數(shù)在水文分析中的應(yīng)用,多變量的洪水頻率分析已成為近年來的研究熱點(diǎn)[1-4].與單變量重現(xiàn)期相比,多變量重現(xiàn)期的定義則要復(fù)雜得多,“或”和“且”定義法[5]是目前最常用的兩種多變量重現(xiàn)期定義方法,但在安全事件與危險(xiǎn)事件的判定上,兩者都存在局限性[6].Salvadori等[6]在“或”和“且”多變量重現(xiàn)期定義法的基礎(chǔ)上,基于安全與危險(xiǎn)域的概念提出了新的多變量重現(xiàn)期定義——二次重現(xiàn)期,相關(guān)概念與定義已在海岸工程設(shè)計(jì)的研究中有了初步的應(yīng)用[7-8].

    本文以東江下游博羅站點(diǎn)為例,介紹多變量重現(xiàn)期新的定義與計(jì)算方法,通過分析傳統(tǒng)重現(xiàn)期與二次重現(xiàn)期在安全與危險(xiǎn)域識(shí)別的差異性,探討二次重現(xiàn)期在洪水風(fēng)險(xiǎn)應(yīng)用的合理性與可靠性.本文的目的在于理論方法的介紹與探討,為洪水風(fēng)險(xiǎn)管理與工程設(shè)計(jì)的研究提供一種新的思路和方法.

    1 洪水特征屬性

    洪水事件具有多個(gè)特征屬性,其中最重要的是歷時(shí)、洪峰和洪量.一場洪水的選取關(guān)鍵在于辨別洪水開始和結(jié)束的時(shí)間,當(dāng)流量開始增加并超過了基流量時(shí)即可認(rèn)為是洪水開始發(fā)生,當(dāng)流量又退落到基流量時(shí)即意味著洪水的結(jié)束.然而,由于每一場洪水的基流量都不一樣,實(shí)際上要準(zhǔn)確地辨認(rèn)出基流量是十分困難的,一般可以將流量突然并持續(xù)增加的時(shí)刻看作是洪水開始的時(shí)間,將流量突然下降并退落至起漲流量或穩(wěn)定流量的時(shí)刻看作是洪水結(jié)束的時(shí)間[9].

    洪水的特征屬性可以按照流量過程線來提取,如圖1所示,流量從突然增加到退落至穩(wěn)定所經(jīng)歷的時(shí)間為洪水歷時(shí)D(d),其間最大的流量即為洪峰流量Q(m3/s),所增加的水量為洪水總量(洪量)V(斜線部分面積,日m3/s).歷時(shí)和洪量可分別根據(jù)下式計(jì)算[9]:

    D=te-ts

    (1)

    (2)

    式中,ts和te分別為洪水開始和結(jié)束的時(shí)間,q為日流量觀測序列.

    圖1 洪水過程及相應(yīng)的特征屬性Fig.1 The flood hydrograph and characteristics

    2 邊緣分布與聯(lián)合分布

    洪水頻率分析中的樣本序列主要通過年最大抽樣的方法從每年發(fā)生的洪水事件中選取,選取的年最大洪水事件可認(rèn)為是極值事件,根據(jù)極值理論,極值事件應(yīng)服從極值分布[10].極值理論中的單元極值分布為廣義極值分布(GEV),其分布函數(shù)為:

    (3)

    式中,μ、σ和k分別為位置、尺度和形狀參數(shù).

    與傳統(tǒng)的多變量分布函數(shù)相比,Copula函數(shù)具有簡單靈活的結(jié)構(gòu)特點(diǎn),并且不受邊緣分布的限制,是一種可以適用于多個(gè)變量間不同相關(guān)結(jié)構(gòu)的多變量分布函數(shù)[11].Archimedean copula是水文分析中最常用的Copula函數(shù)族,其基本形式為:

    (4)

    式中,uj∈[0, 1](j>1)為邊緣分布,φθ為Archimedean copula生成元,θ為參數(shù).當(dāng)C滿足下式時(shí),則C是一個(gè)多變量極值分布[13].

    (5)

    可以證明Gumbel-Hougaard copula滿足上式,并且是Archimedean copula中唯一的多變量極值copula[12-13],其結(jié)構(gòu)形式為:

    (6)

    僅當(dāng)邊緣分布和Copula函數(shù)均具有極值性質(zhì)時(shí),構(gòu)造的聯(lián)合分布才是極值分布[13],可用于極值事件的頻率分析中.

    3 多變量重現(xiàn)期與分位值

    3.1 傳統(tǒng)多變量重現(xiàn)期

    在傳統(tǒng)的多變量重現(xiàn)期定義中,E∨對(duì)應(yīng)的重現(xiàn)期為“或”重現(xiàn)期T∨,E∧對(duì)應(yīng)的重現(xiàn)期則為“且”重現(xiàn)期T∧.年最大洪水事件的3變量“或”重現(xiàn)期和“且”重現(xiàn)期可分別由以下兩式計(jì)算:

    (7)

    (8)

    假設(shè)E的安全率為Psafe(危險(xiǎn)率則為1-Psafe),在Psafe=p的條件下,以兩變量的情形進(jìn)行闡述,“或”重現(xiàn)期和“且”重現(xiàn)期對(duì)安全與危險(xiǎn)域的識(shí)別如圖2所示.在相同安全率的情況下,“或”重現(xiàn)期和“且”重現(xiàn)期標(biāo)準(zhǔn)下不同的分位值組合對(duì)應(yīng)的安全與危險(xiǎn)域都不一樣.

    圖2 “或”重現(xiàn)期(a)和“且”重現(xiàn)期(b)的安全與危險(xiǎn)域識(shí)別(灰色部分為危險(xiǎn)域)Fig.2 Identification of safe and dangerous regions for “OR” return period(a) and “AND” return period(b) (the grey part represents the dangerous one)

    一般地,設(shè)計(jì)的重現(xiàn)期越大,相應(yīng)的危險(xiǎn)率就越小,設(shè)置的災(zāi)害防御標(biāo)準(zhǔn)遭遇的風(fēng)險(xiǎn)就越低.就此而言,“或”重現(xiàn)期和“且”重現(xiàn)期識(shí)別的安全與危險(xiǎn)域顯然具有局限性[6].如圖3a所示,對(duì)于“或”重現(xiàn)期,A點(diǎn)處的重現(xiàn)期大于B點(diǎn)(p1>p2),而A點(diǎn)識(shí)別的危險(xiǎn)域卻有一部分在B點(diǎn)識(shí)別的危險(xiǎn)域之外(深灰色部分),即當(dāng)事件E*發(fā)生時(shí),若以B點(diǎn)的組合值作為災(zāi)害防御的標(biāo)準(zhǔn)則可以認(rèn)為E*是安全的,E*造成的效應(yīng)在防御

    的標(biāo)準(zhǔn)之下;若以重現(xiàn)期更大的A點(diǎn)組合值作為災(zāi)害防御的標(biāo)準(zhǔn)則E*顯然是危險(xiǎn)事件.高風(fēng)險(xiǎn)的防御標(biāo)準(zhǔn)識(shí)別的安全事件卻被更低風(fēng)險(xiǎn)的防御標(biāo)準(zhǔn)識(shí)別為危險(xiǎn)事件,這顯然是矛盾的.對(duì)于“且”重現(xiàn)期(圖3b),A點(diǎn)處的重現(xiàn)期小于B點(diǎn)(p1

    圖3 “或”重現(xiàn)期(a)和“且”重現(xiàn)期(b)對(duì)安全危險(xiǎn)域識(shí)別的局限性Fig.3 Limitations on the identification of safe and dangerous regions for “OR” return period(a) and “AND” return period(b)

    3.2 分位值曲面

    3.3 二次重現(xiàn)期

    (9)

    TK稱為二次重現(xiàn)期(Secondary return period)或Kendall重現(xiàn)期,KC為Kendall分布函數(shù)[16]:

    (10)

    (11)

    對(duì)于三維Gumbel-Hougaard copula,Kendall分布函數(shù)可以表達(dá)為:

    (12)

    由公式(7)、(8)和(9)以及C的非遞減性可知,單變量重現(xiàn)期、“或”重現(xiàn)期、“且”重現(xiàn)期和二次重現(xiàn)期四者的關(guān)系[14,17]為:T∨≤T≤TK≤T∧,當(dāng)且僅當(dāng)uj=0時(shí),等號(hào)成立.

    圖4 二次重現(xiàn)期的安全與危險(xiǎn)域識(shí)別 (灰色部分為危險(xiǎn)域)Fig.4 Identification of safe and dangerous regions for secondary return period (the grey part represents the dangerous one)

    3.4 多變量設(shè)計(jì)值

    (13)

    f(d,q,v)=c(ud,uq,uv)fd(d)fq(q)fv(v)

    (14)

    式中,c為三維Gumbel-Hougaard copula的概率密度函數(shù),fd(d)、fq(q)和fv(v)則分別為歷時(shí)、洪峰和洪量的概率密度函數(shù).

    綜上所述,二次重現(xiàn)期條件下的洪水多變量分位值估算步驟可以總結(jié)為:

    (2) 利用式(4),計(jì)算所有使C(ud,uq,uv)=p成立的概率組合(ud,uq,uv);

    (3) 由式(15)計(jì)算出使f(d,q,v)達(dá)到最大值的一組(ud,uq,uv);

    (4) 最后分別再根據(jù)邊緣分布的反函數(shù)推求出歷時(shí)、洪峰和洪量的分位值(dml,qml,vml),dml=F-1(ud)、qml=F-1(uq)、vml=F-1(uv).

    4 案例研究

    博羅站位于東江下游,是東江流域重要的水文控制站點(diǎn)之一.本文選用博羅站1953-2010年的日流量觀測數(shù)據(jù)作為研究資料,依據(jù)年最大極值理論,首先提取每年的最大日流量作為該年極端洪水事件的洪峰流量,以洪峰流量為基準(zhǔn),再按流量過程線提取并計(jì)算出該場洪水的歷時(shí)與洪量,總共提取了58場洪水作為博羅站點(diǎn)多變量洪水頻率分析的樣本序列.歷時(shí)與洪峰、洪量樣本兩兩間的Kendall秩相關(guān)系數(shù)表明洪水特征屬性間存在較強(qiáng)的相關(guān)性,其中,洪峰與洪量的相關(guān)性最大(0.665),歷時(shí)與洪量次之(0.623),歷時(shí)與洪峰則相對(duì)較小(0.355).

    為了檢驗(yàn)GEV和Gumbel-Hougaard copula構(gòu)造洪水事件聯(lián)合分布的合理性,選用水文分析中常用的對(duì)數(shù)正態(tài)分布(GNO)、廣義邏輯斯特分布(GLO)和皮爾遜三型分布(PE3),Clayton和Frank copula[12]分別對(duì)歷時(shí)、洪峰和洪量樣本進(jìn)行擬合優(yōu)度檢驗(yàn),并與GEV和Gumbel-Hougaard copula進(jìn)行對(duì)比.選用的分布參數(shù)均采用極大似然法進(jìn)行估計(jì),擬合優(yōu)度檢驗(yàn)則依據(jù)Akaike信息準(zhǔn)則(AIC)[7],AIC值越小,代表擬合優(yōu)度越高.選用的備選分布都通過了K-S檢驗(yàn),表明選用的分布都可用于樣本的擬合.各種分布的擬合優(yōu)度AIC值如表1所示,雖然擬合優(yōu)度檢驗(yàn)的結(jié)果并沒有顯示GEV和Gumbel-Hougaard copula對(duì)所有歷時(shí)、洪峰和洪量樣本以及它們間組合的擬合都是最優(yōu)的(GEV對(duì)歷時(shí)和Gumbel-Hougaard copula對(duì)歷時(shí)與洪量組合的擬合為次優(yōu)),但從總體角度和出于極值分布的考慮,我們選用GEV和Gumbel-Hougaard copula分別作為邊緣分布和多變量分布來構(gòu)造洪水特征屬性間的聯(lián)合分布.圖5和圖6所示的GEV和Gumbel-Hougaard copula對(duì)歷時(shí)、洪峰和洪量樣本的擬合效果亦表明兩者用于構(gòu)造洪水聯(lián)合分布是合理的.

    表1 邊緣分布與copula函數(shù)的擬合優(yōu)度AIC值Tab.1 The AIC values for different marginal distributions and copulas

    圖5 GEV對(duì)洪水歷時(shí)、洪峰和洪量樣本的擬合優(yōu)度Fig.5 The goodness of fit of GEV for flood duration, peak and volume

    圖6 Gumbel-Hougaard copula對(duì)洪水歷時(shí)、洪峰和洪量樣本的擬合優(yōu)度Fig.6 The goodness of fit of Gumbel-Hougaard copula for flood duration, peak and volume

    歷時(shí)、洪峰和洪量保證率分別為96%、98%、99%和99.5%時(shí)的單變量重現(xiàn)期、“或”重現(xiàn)期、“且”重現(xiàn)期和二次重現(xiàn)期計(jì)算值如表2所示,在同一保證率水平下(保證率不為0),二次重現(xiàn)期與“且”重現(xiàn)期相差較小,與“或”重現(xiàn)期則相差較大;保證率越小,二次重現(xiàn)期與“或”重現(xiàn)期和“且”重現(xiàn)期之間的差值則越小,隨著保證率的增大,三者的差值逐漸增大.根據(jù)3.4所述步驟計(jì)算的25、50、100和200年重現(xiàn)期歷時(shí)、洪峰和洪量的多變量設(shè)計(jì)值以及單獨(dú)推算的單變量設(shè)計(jì)值如表3所示,單獨(dú)推算的歷時(shí)與洪量設(shè)計(jì)值均比多變量推算的大,而多變量洪峰設(shè)計(jì)值則要大于單變量設(shè)計(jì)值,并且隨著重現(xiàn)期的增大,兩者的差值逐漸減小.一場洪水造成的效應(yīng)并不是由單個(gè)特征屬性決定的,而是多個(gè)特征屬性共同作用的結(jié)果,如高洪峰伴隨大洪量的洪水造成的影響和破壞要比高洪峰伴隨小洪量的洪水大,長歷時(shí)高洪峰的洪水造成的影響和破壞要比短歷時(shí)高洪峰的洪水大.再者,歷時(shí)、洪峰與洪量之間的組合以及出現(xiàn)的概率又取決于它們間的相關(guān)結(jié)構(gòu)和聯(lián)合分布.因此,相比于單變量設(shè)計(jì)值,考慮了洪水多個(gè)屬性聯(lián)合特征的多變量設(shè)計(jì)值為防洪工程設(shè)計(jì)提供了更加全面和可靠的參考信息.

    表2 同一保證率水平下的單變量、“或”、“且”和二次重現(xiàn)期計(jì)算值Tab.2 The calculation values of univariate, “OR, “AND” and secondary return periods in the same probability level

    表3 25、50、100和200年重現(xiàn)期歷時(shí)、洪峰和洪量的單變量與多變量設(shè)計(jì)值Tab.3 The design values for flood duration, peak and volume in a multivariate framework under the secondary return periods of 25, 50, 100 and 200 years, respectively

    5 結(jié)論

    1) 傳統(tǒng)的“或”和“且”多變量重現(xiàn)期對(duì)安全與危險(xiǎn)域的識(shí)別具有局限性,重現(xiàn)期與危險(xiǎn)域范圍大小的矛盾會(huì)造成對(duì)安全事件與危險(xiǎn)事件的錯(cuò)誤識(shí)別.

    2) 根據(jù)C(ud,uq,uv)=p臨界條件并利用Kendall函數(shù)定義的二次重現(xiàn)期,具有同一臨界水平下任意歷時(shí)、洪峰與洪量組合對(duì)安全與危險(xiǎn)域識(shí)別的一致性,重現(xiàn)期越大,對(duì)應(yīng)的危險(xiǎn)域則越小,避免了對(duì)安全事件與危險(xiǎn)事件的錯(cuò)誤識(shí)別,更有利于指導(dǎo)洪水風(fēng)險(xiǎn)的管理.

    3) 給定二次重現(xiàn)期條件下的洪水事件可以具有不同的歷時(shí)、洪峰與洪量組合,在工程設(shè)計(jì)中,與每個(gè)特征屬性同時(shí)取較大值相比,依據(jù)出現(xiàn)概率最大的原則推算的設(shè)計(jì)值組合可以滿足以較低成本承受較大風(fēng)險(xiǎn)的追求.相比于單變量設(shè)計(jì)值,考慮了洪水多個(gè)屬性聯(lián)合特征的多變量設(shè)計(jì)值為防洪工程設(shè)計(jì)提供了更加全面和可靠的參考信息.

    [1] 郭生練,閆寶偉,肖 義等.Copula函數(shù)在多變量水文分析計(jì)算中的應(yīng)用及研究進(jìn)展.水文,2008,28(3):1-6.

    [2] Zhang L, Singh VP. Bivariate flood frequency analysis using the copula method.JournalofHydrologicEngineering, 2006, 11(2):150-164.

    [3] Zhang L, Singh VP. Trivariate flood frequency analysis using the Gumbel-Hougaard copula.JournalofHydrologicEngineering, 2007, 12(4):431-439.

    [4] 侯蕓蕓,宋松柏,趙麗娜等.基于Copula函數(shù)的3變量洪水頻率研究.西北農(nóng)林科技大學(xué)學(xué)報(bào):自然科學(xué)版,2010,38(2):219-228.

    [5] Shiau JT. Fitting drought duration and severity with two-dimensional copulas.WaterResourcesManagement, 2006, 20:795-815.

    [6] Salvadori G, De Michele C, Durante F. On the return period and design in a multivariate framework.HydrologyandEarthSystemSciences, 2011, 15:3293-3305.

    [7] Corbella S, Stretch DD. Multivariate return periods of sea storms for coastal erosion risk assessment.NaturalHazardsandEarthSystemSciences, 2012, 12:2699-2708.

    [8] Salvadori G, Tomasicchio GR, Alessandro FD. Multivariate approach to design coastal and off-shore structures.JournalofCoastalResearch, 2013, 65:386-391.

    [9] Yue S, Ouarda TBMJ, Bobée Betal. The Gumbel mixed model for flood frequency analysis.JournalofHydrology, 1999, 226(1):88-100.

    [10] Salvadori G, de Michele C, Kottegoda Netal. Extremes in Nature: An approach using copulas. Dordrecht: Springer, 2007:15-30.

    [11] Favre AC, Adlouni SE, Perreault Letal. Multivariate hydrological frequency analysis using copulas.WaterResourcesResearch, 2004, 40: W01101. doi:10.1029/2003WR002456.

    [12] Nelson RB. An introduction to copulas: Second edition. New York: Springer, 2006:115-143.

    [13] Salvadori G, de Michele C. Multivariate multiparameter extreme value models and return periods: A copula approach.WaterResourceResearch, 2010, 46: W10501. doi:10.1029/2009WR009040.

    [14] Salvadori G, de Michele C. Frequency analysis via copulas: Theoretical aspects and applications to hydrological events.WaterResourcesResearch, 2004, 40: W12511. doi:10.1029/2004WR003133.

    [15] Chebana F, Ouarda TBMJ. Multivariate quantiles in hydrological frequency analysis.Environmetrics, 2011, 22:63-78.

    [16] Genest C, Quessy JF, Rémillard B. Goodness-of-fit procedures for copula model based on the probability integral transformation.ScandinavianJournalofStatistics, 2006, 33(2):337-366.

    [17] Graler B, van den Berg MJ, Vandenberghe Setal. Multivariate return periods in hydrology: a critical and practical review focusing on synthetic design hydrograph estimation.HydrologyandEarthSystemSciences, 2013, 17: 1281-1296.

    Multivariate flood risk assessment based on the secondary return period

    HUANG Qiang & CHEN Zishen

    (DepartmentofWaterResourcesandEnvironment,SunYat-senUniversity,Guangzhou510275,P.R.China)

    Frequency analysis is a useful tool for flood risk assessment, but the definition and calculation of return periods and design values in a multivariate framework are difficult tasks. In this paper, by constructing the joint distribution of flood duration, peak discharge and volume, we introduced the definitions of multivariate “OR”, “AND” and secondary return periods and discuss their differences in identifying the safe and dangerous regions in a critical level. The rationality and reliability of information provided from them in flood risk management and engineering design were then analyzed, respectively. The traditional “OR-AND” approach in multivariate return periods definition is limited in the identification of safe and dangerous regions, while the secondary return period based on the Kendall measure would be more rational. Consequently, mistakes in the recognition of safe and dangerous event can be avoided by the use of second return period, which, obviously, is better for the flood risk management. For a given second return period, the combination of flood duration, peak discharge and volume design values calculated from the most-likely approach, may have an advantage to meet the requirement that with low cost to bear greater risks in engineering design. Furthermore, the information provided from multivariate design values would be more considerate and reliable than the univariate ones.

    Multivariate flood characteristics; extreme value distribution; safe and dangerous regions; Kendall function; secondary return period; multivariate design values

    *國家自然科學(xué)基金項(xiàng)目(41371498)和廣東水利創(chuàng)新基金項(xiàng)目(2009-41)聯(lián)合資助.2014-02-27收稿;2014-08-11收修改稿.黃強(qiáng)(1989~),男,博士研究生;E-mail:huangq52@mail2.sysu.edu.cn.

    **通信作者;E-mail:eesczs@mail.sysu.edu.cn.

    猜你喜歡
    洪量洪峰歷時(shí)
    基于SPA 的北江流域峰量關(guān)系研究
    陜西水利(2023年12期)2023-12-19 03:28:32
    量詞“只”的形成及其歷時(shí)演變
    常用詞“怠”“惰”“懶”的歷時(shí)演變
    對(duì)《紅樓夢(mèng)》中“不好死了”與“……好的”的歷時(shí)考察
    古今字“兌”“說”“悅”“敚”歷時(shí)考察
    遼河干流主要控制站近75年最大洪峰及洪量變化特征分析研究
    淡定!
    解禁洪峰
    地火(2014年4期)2014-03-01 01:55:30
    適用于電算的設(shè)計(jì)洪水過程線放縮方法
    我要搜黄色片| 国内精品美女久久久久久| 日韩大尺度精品在线看网址| 日日摸夜夜添夜夜添av毛片| 小说图片视频综合网站| 久久99蜜桃精品久久| 九九久久精品国产亚洲av麻豆| 国产国拍精品亚洲av在线观看| 我要搜黄色片| 蜜桃亚洲精品一区二区三区| 国产精品麻豆人妻色哟哟久久 | 中文字幕久久专区| 国产高潮美女av| 美女黄网站色视频| 春色校园在线视频观看| 午夜精品一区二区三区免费看| 欧美高清成人免费视频www| 国产高清不卡午夜福利| 色噜噜av男人的天堂激情| 免费观看人在逋| 久久精品国产亚洲网站| 久久久久久久久中文| 两个人的视频大全免费| 欧美日韩在线观看h| 我的女老师完整版在线观看| 91麻豆精品激情在线观看国产| 国内精品一区二区在线观看| 亚洲激情五月婷婷啪啪| 色噜噜av男人的天堂激情| 国产精品爽爽va在线观看网站| 最近手机中文字幕大全| 99热6这里只有精品| 中出人妻视频一区二区| 大又大粗又爽又黄少妇毛片口| 中国美白少妇内射xxxbb| 久久久久久九九精品二区国产| 国产av在哪里看| 国产成年人精品一区二区| 日日干狠狠操夜夜爽| 99热这里只有是精品在线观看| 九九爱精品视频在线观看| 小蜜桃在线观看免费完整版高清| 精品少妇黑人巨大在线播放 | 久久久久国产网址| 变态另类成人亚洲欧美熟女| 99久久久亚洲精品蜜臀av| 美女国产视频在线观看| 成人特级av手机在线观看| 欧美日韩乱码在线| 免费看美女性在线毛片视频| 国产黄色小视频在线观看| 精品久久久久久久末码| 能在线免费观看的黄片| 校园春色视频在线观看| 国产色婷婷99| 国产精品乱码一区二三区的特点| 久久精品91蜜桃| 亚洲乱码一区二区免费版| 国产精品免费一区二区三区在线| 狂野欧美激情性xxxx在线观看| 亚洲精品影视一区二区三区av| 国产麻豆成人av免费视频| av黄色大香蕉| 亚洲经典国产精华液单| 国产三级在线视频| 日韩欧美三级三区| 听说在线观看完整版免费高清| 欧美xxxx黑人xx丫x性爽| 99热这里只有精品一区| 在线观看免费视频日本深夜| 久久精品影院6| 麻豆国产av国片精品| 青春草视频在线免费观看| 12—13女人毛片做爰片一| 国产一区亚洲一区在线观看| 秋霞在线观看毛片| 日韩,欧美,国产一区二区三区 | 一级二级三级毛片免费看| 久久久久性生活片| 日韩,欧美,国产一区二区三区 | 99久久无色码亚洲精品果冻| 欧美一级a爱片免费观看看| 亚洲电影在线观看av| 午夜福利成人在线免费观看| 亚洲av二区三区四区| 久久精品久久久久久噜噜老黄 | 国产在视频线在精品| 国产老妇女一区| 色吧在线观看| 国产中年淑女户外野战色| 91狼人影院| 亚洲国产精品sss在线观看| 亚洲欧美日韩卡通动漫| 又黄又爽又刺激的免费视频.| 日本黄色视频三级网站网址| 变态另类丝袜制服| 一本一本综合久久| 蜜桃久久精品国产亚洲av| 国产精品一区二区三区四区免费观看| 美女大奶头视频| 精品久久久久久久久av| 久久精品国产亚洲网站| 亚洲经典国产精华液单| 哪里可以看免费的av片| 最近的中文字幕免费完整| 中国美白少妇内射xxxbb| 美女大奶头视频| 国产成人a∨麻豆精品| 国产亚洲av片在线观看秒播厂 | 国产探花在线观看一区二区| 亚洲成人精品中文字幕电影| 一区二区三区免费毛片| 99热精品在线国产| 国产真实伦视频高清在线观看| 女人十人毛片免费观看3o分钟| 别揉我奶头 嗯啊视频| 女同久久另类99精品国产91| 国产精品日韩av在线免费观看| 99热6这里只有精品| 特级一级黄色大片| 我的女老师完整版在线观看| 韩国av在线不卡| 亚洲最大成人中文| 亚洲av成人av| 国产69精品久久久久777片| 少妇熟女aⅴ在线视频| 亚洲欧美成人综合另类久久久 | 久久久精品大字幕| av在线观看视频网站免费| 91在线精品国自产拍蜜月| 色尼玛亚洲综合影院| 精品熟女少妇av免费看| 深爱激情五月婷婷| 内地一区二区视频在线| 22中文网久久字幕| 欧美激情在线99| 国产久久久一区二区三区| 国产精品国产高清国产av| 精品一区二区三区人妻视频| 国产综合懂色| 成人高潮视频无遮挡免费网站| 在线a可以看的网站| 亚洲精品日韩av片在线观看| 麻豆国产97在线/欧美| 亚洲国产精品成人久久小说 | 国产精品久久久久久精品电影| 亚洲五月天丁香| 在线播放国产精品三级| 少妇高潮的动态图| 搡女人真爽免费视频火全软件| 日日摸夜夜添夜夜爱| 99久国产av精品| 内射极品少妇av片p| 成人三级黄色视频| 日本免费a在线| 中出人妻视频一区二区| 欧美最黄视频在线播放免费| 永久网站在线| 日日撸夜夜添| 亚洲婷婷狠狠爱综合网| 亚洲不卡免费看| 九草在线视频观看| 免费观看a级毛片全部| 一区二区三区四区激情视频 | 男的添女的下面高潮视频| 国产精品,欧美在线| videossex国产| 欧美色欧美亚洲另类二区| 免费av观看视频| 人人妻人人澡欧美一区二区| 五月玫瑰六月丁香| 日韩视频在线欧美| 免费av观看视频| 国产白丝娇喘喷水9色精品| 天天躁夜夜躁狠狠久久av| 97超视频在线观看视频| 日韩欧美精品免费久久| 午夜久久久久精精品| 久久鲁丝午夜福利片| 色哟哟哟哟哟哟| 久久婷婷人人爽人人干人人爱| 一级二级三级毛片免费看| 可以在线观看的亚洲视频| 国产探花在线观看一区二区| 一区二区三区高清视频在线| 丰满人妻一区二区三区视频av| 能在线免费看毛片的网站| 国产av一区在线观看免费| 91狼人影院| 亚洲精品乱码久久久久久按摩| а√天堂www在线а√下载| 老师上课跳d突然被开到最大视频| 大型黄色视频在线免费观看| 一边摸一边抽搐一进一小说| 国产精品99久久久久久久久| 有码 亚洲区| 欧美成人免费av一区二区三区| 色哟哟哟哟哟哟| 69人妻影院| 成人综合一区亚洲| 精华霜和精华液先用哪个| 给我免费播放毛片高清在线观看| 日韩,欧美,国产一区二区三区 | 亚洲精品久久国产高清桃花| 久久这里有精品视频免费| 在线观看一区二区三区| 欧美色欧美亚洲另类二区| 变态另类成人亚洲欧美熟女| 69人妻影院| 联通29元200g的流量卡| 国产精品一及| 久久久久性生活片| 成人毛片60女人毛片免费| 两个人视频免费观看高清| 中出人妻视频一区二区| 成人毛片60女人毛片免费| 一级毛片久久久久久久久女| 亚洲在线观看片| 91av网一区二区| 亚洲四区av| 一个人看视频在线观看www免费| 美女黄网站色视频| 欧美高清性xxxxhd video| 日本免费一区二区三区高清不卡| 99久久精品国产国产毛片| 一个人看视频在线观看www免费| 亚洲国产欧美在线一区| 国产乱人偷精品视频| 国产三级中文精品| 亚洲精品影视一区二区三区av| 欧美高清成人免费视频www| 国产精品一及| 久久久国产成人精品二区| 日本黄色视频三级网站网址| 亚洲婷婷狠狠爱综合网| 伊人久久精品亚洲午夜| 淫秽高清视频在线观看| 国产又黄又爽又无遮挡在线| eeuss影院久久| 成人av在线播放网站| 欧美不卡视频在线免费观看| 青春草亚洲视频在线观看| 老师上课跳d突然被开到最大视频| 国产三级在线视频| 久久久色成人| 国产精品.久久久| 狂野欧美激情性xxxx在线观看| 老女人水多毛片| 18禁黄网站禁片免费观看直播| 麻豆成人午夜福利视频| 天堂中文最新版在线下载 | 深爱激情五月婷婷| 久久人妻av系列| 91午夜精品亚洲一区二区三区| 国产成人精品婷婷| 亚洲欧美清纯卡通| 99久久人妻综合| 变态另类成人亚洲欧美熟女| 变态另类丝袜制服| 日本撒尿小便嘘嘘汇集6| 伊人久久精品亚洲午夜| 亚洲av免费在线观看| 非洲黑人性xxxx精品又粗又长| 99久久九九国产精品国产免费| 国产成年人精品一区二区| 少妇熟女欧美另类| 国产成人freesex在线| 日韩三级伦理在线观看| 久久99热6这里只有精品| 人人妻人人澡人人爽人人夜夜 | 精品久久久久久久人妻蜜臀av| 久久精品久久久久久噜噜老黄 | 成熟少妇高潮喷水视频| 三级国产精品欧美在线观看| 亚洲欧美日韩卡通动漫| 亚洲av中文字字幕乱码综合| 久久精品夜夜夜夜夜久久蜜豆| 国产人妻一区二区三区在| 欧美激情久久久久久爽电影| 色哟哟·www| 欧美高清性xxxxhd video| 日本黄色视频三级网站网址| 爱豆传媒免费全集在线观看| 可以在线观看的亚洲视频| 91狼人影院| 国产成人91sexporn| 99热6这里只有精品| 亚洲人成网站在线播放欧美日韩| 欧美性猛交黑人性爽| 亚洲第一电影网av| 亚洲真实伦在线观看| 18禁裸乳无遮挡免费网站照片| 在线播放无遮挡| 国产成年人精品一区二区| 国产精品久久久久久久久免| 一进一出抽搐gif免费好疼| 91久久精品国产一区二区成人| 联通29元200g的流量卡| 亚洲美女视频黄频| 日本撒尿小便嘘嘘汇集6| 熟女电影av网| 欧美日韩乱码在线| 大又大粗又爽又黄少妇毛片口| 成人国产麻豆网| 人妻系列 视频| 亚洲一区二区三区色噜噜| 最近中文字幕高清免费大全6| 日韩一区二区视频免费看| 久久久成人免费电影| 麻豆精品久久久久久蜜桃| 网址你懂的国产日韩在线| 亚洲精品粉嫩美女一区| 又爽又黄无遮挡网站| av.在线天堂| 偷拍熟女少妇极品色| 99久久成人亚洲精品观看| 又粗又爽又猛毛片免费看| 亚洲国产精品成人综合色| 亚洲三级黄色毛片| 亚洲欧美精品专区久久| 日韩 亚洲 欧美在线| 三级国产精品欧美在线观看| 久久久色成人| 亚洲中文字幕一区二区三区有码在线看| 欧美性猛交╳xxx乱大交人| 亚洲aⅴ乱码一区二区在线播放| 日韩强制内射视频| 国产v大片淫在线免费观看| 国产午夜福利久久久久久| 亚洲av第一区精品v没综合| 非洲黑人性xxxx精品又粗又长| 亚洲性久久影院| 亚洲自偷自拍三级| 久久午夜亚洲精品久久| 国产精品不卡视频一区二区| 亚洲精品色激情综合| 国产中年淑女户外野战色| 黄色日韩在线| 亚洲欧洲日产国产| 内地一区二区视频在线| 两个人的视频大全免费| 最近中文字幕高清免费大全6| 国产大屁股一区二区在线视频| 国产一级毛片七仙女欲春2| 听说在线观看完整版免费高清| 国产91av在线免费观看| 身体一侧抽搐| 美女xxoo啪啪120秒动态图| 精品久久久久久久久久久久久| 亚洲精品亚洲一区二区| 亚洲av第一区精品v没综合| 国产精品永久免费网站| 1000部很黄的大片| 国产精品女同一区二区软件| 99热这里只有精品一区| 亚洲一区二区三区色噜噜| 青青草视频在线视频观看| 乱码一卡2卡4卡精品| 日韩精品有码人妻一区| 色尼玛亚洲综合影院| 婷婷六月久久综合丁香| 国产乱人偷精品视频| 99视频精品全部免费 在线| 色播亚洲综合网| 免费无遮挡裸体视频| 少妇猛男粗大的猛烈进出视频 | 成人漫画全彩无遮挡| 中文字幕av在线有码专区| 国产一区二区激情短视频| 国产激情偷乱视频一区二区| 嫩草影院精品99| 亚洲成av人片在线播放无| 噜噜噜噜噜久久久久久91| 国产精品人妻久久久久久| 日韩强制内射视频| 伦精品一区二区三区| 一个人免费在线观看电影| 亚洲在久久综合| 亚洲av第一区精品v没综合| 天天一区二区日本电影三级| 天美传媒精品一区二区| 国产白丝娇喘喷水9色精品| 毛片一级片免费看久久久久| 高清午夜精品一区二区三区 | 只有这里有精品99| 极品教师在线视频| 国产高清激情床上av| 久久精品国产亚洲av天美| 婷婷六月久久综合丁香| 日本一本二区三区精品| 国产av不卡久久| 啦啦啦韩国在线观看视频| 精品久久国产蜜桃| 日韩成人伦理影院| 一区二区三区四区激情视频 | 久久久国产成人免费| 国产爱豆传媒在线观看| 桃色一区二区三区在线观看| 婷婷色综合大香蕉| 国产乱人偷精品视频| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲欧美日韩无卡精品| 亚洲久久久久久中文字幕| 久久久成人免费电影| 国产色爽女视频免费观看| 亚洲熟妇中文字幕五十中出| 人人妻人人看人人澡| 床上黄色一级片| 亚洲国产精品国产精品| 亚洲精品日韩av片在线观看| 久久热精品热| 欧美zozozo另类| 天堂影院成人在线观看| 国产三级在线视频| 国产精品嫩草影院av在线观看| www.av在线官网国产| 国产一区二区在线av高清观看| 国产久久久一区二区三区| 只有这里有精品99| 国产精品嫩草影院av在线观看| 成年版毛片免费区| 日韩强制内射视频| 黄色欧美视频在线观看| 热99在线观看视频| 久久精品久久久久久久性| 人妻夜夜爽99麻豆av| 一区二区三区免费毛片| 伦精品一区二区三区| 国产一区二区在线av高清观看| 中文字幕人妻熟人妻熟丝袜美| 村上凉子中文字幕在线| 可以在线观看的亚洲视频| 日本黄色视频三级网站网址| 少妇猛男粗大的猛烈进出视频 | 国产黄色小视频在线观看| av黄色大香蕉| 日韩av不卡免费在线播放| 欧美xxxx黑人xx丫x性爽| 熟妇人妻久久中文字幕3abv| 变态另类丝袜制服| 成人二区视频| 日本免费一区二区三区高清不卡| 久久国内精品自在自线图片| 婷婷精品国产亚洲av| 亚洲av第一区精品v没综合| 亚洲成av人片在线播放无| 综合色av麻豆| 九草在线视频观看| 国产亚洲av嫩草精品影院| 大型黄色视频在线免费观看| 91精品国产九色| 国产精品久久久久久精品电影| 国产色婷婷99| 青青草视频在线视频观看| 看十八女毛片水多多多| 亚洲欧美日韩高清专用| 国产爱豆传媒在线观看| 极品教师在线视频| 免费不卡的大黄色大毛片视频在线观看 | 亚洲av熟女| 插逼视频在线观看| 狂野欧美白嫩少妇大欣赏| 日韩欧美精品v在线| 99热这里只有精品一区| 国产亚洲91精品色在线| 中国美白少妇内射xxxbb| 毛片一级片免费看久久久久| 国产女主播在线喷水免费视频网站 | 成人特级av手机在线观看| 日韩国内少妇激情av| 日韩欧美精品免费久久| 亚洲精品自拍成人| 亚洲久久久久久中文字幕| 黑人高潮一二区| 人人妻人人看人人澡| 免费看美女性在线毛片视频| 欧美区成人在线视频| 麻豆成人午夜福利视频| 日本色播在线视频| 国产成人a区在线观看| 丰满乱子伦码专区| 少妇高潮的动态图| 久久久国产成人免费| 国产精品无大码| 级片在线观看| 狂野欧美激情性xxxx在线观看| 少妇高潮的动态图| 国产淫片久久久久久久久| 国产午夜精品久久久久久一区二区三区| 精品久久久久久成人av| 床上黄色一级片| 一个人看视频在线观看www免费| 最后的刺客免费高清国语| 久久久久久九九精品二区国产| 嘟嘟电影网在线观看| 亚洲一级一片aⅴ在线观看| 18+在线观看网站| 69av精品久久久久久| 永久网站在线| 久久精品国产99精品国产亚洲性色| 亚洲av男天堂| 国产成人午夜福利电影在线观看| 卡戴珊不雅视频在线播放| 亚洲欧美精品自产自拍| 干丝袜人妻中文字幕| 色吧在线观看| 色5月婷婷丁香| 最后的刺客免费高清国语| 亚洲精品乱码久久久久久按摩| 99热这里只有精品一区| 亚洲欧洲日产国产| 亚洲色图av天堂| 在现免费观看毛片| 欧美日韩在线观看h| 观看免费一级毛片| 我要看日韩黄色一级片| 美女脱内裤让男人舔精品视频 | 成人毛片60女人毛片免费| 久久久国产成人精品二区| 亚洲一级一片aⅴ在线观看| 久久亚洲国产成人精品v| 97热精品久久久久久| 亚洲精品亚洲一区二区| 直男gayav资源| 色哟哟·www| 欧美不卡视频在线免费观看| 可以在线观看的亚洲视频| avwww免费| 国产精品野战在线观看| 欧美成人a在线观看| 欧美三级亚洲精品| 久久精品国产99精品国产亚洲性色| 爱豆传媒免费全集在线观看| 熟妇人妻久久中文字幕3abv| 午夜精品在线福利| 精华霜和精华液先用哪个| 中文字幕久久专区| 99国产精品一区二区蜜桃av| 国内少妇人妻偷人精品xxx网站| 精品久久久久久成人av| 欧美性猛交黑人性爽| 在线观看66精品国产| 国产中年淑女户外野战色| 美女国产视频在线观看| av在线亚洲专区| 高清毛片免费观看视频网站| 欧美日本亚洲视频在线播放| 国产成人a区在线观看| 搡老妇女老女人老熟妇| 亚洲精品亚洲一区二区| 别揉我奶头 嗯啊视频| 综合色av麻豆| 在线天堂最新版资源| 99热全是精品| 久久久久久久久大av| 色5月婷婷丁香| 国内精品宾馆在线| 午夜福利在线在线| 婷婷色av中文字幕| 午夜a级毛片| 中国美女看黄片| 国产人妻一区二区三区在| 精品人妻偷拍中文字幕| av女优亚洲男人天堂| 久久人人精品亚洲av| 直男gayav资源| 变态另类丝袜制服| 91精品一卡2卡3卡4卡| 国产高清有码在线观看视频| h日本视频在线播放| 国产成人影院久久av| 亚洲精品亚洲一区二区| 亚洲av电影不卡..在线观看| 国产一区亚洲一区在线观看| 日韩成人伦理影院| 好男人视频免费观看在线| 国产高潮美女av| 蜜桃亚洲精品一区二区三区| 丰满的人妻完整版| 国产在线男女| 一区二区三区四区激情视频 | 床上黄色一级片| 久久亚洲国产成人精品v| 毛片一级片免费看久久久久| 99久久久亚洲精品蜜臀av| 免费看日本二区| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 欧美xxxx黑人xx丫x性爽| 亚洲成人av在线免费| 麻豆国产av国片精品| 久久99精品国语久久久| 校园春色视频在线观看| 亚洲自偷自拍三级| 久久久欧美国产精品| 久久6这里有精品| 亚洲电影在线观看av| 国产中年淑女户外野战色| 欧美日韩精品成人综合77777| 一本精品99久久精品77| 国产精品日韩av在线免费观看| 两性午夜刺激爽爽歪歪视频在线观看| 免费搜索国产男女视频| 亚洲中文字幕一区二区三区有码在线看| 久久人人爽人人片av| 欧美日韩精品成人综合77777| 男女那种视频在线观看| 欧美一区二区亚洲| av福利片在线观看| 在线免费十八禁| 欧美日韩在线观看h| 老女人水多毛片| 深爱激情五月婷婷| 三级经典国产精品| 麻豆国产av国片精品| 国产精品无大码|