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

    連續(xù)API產(chǎn)匯流模型的改進(jìn)及應(yīng)用

    2018-01-04 18:11:20
    水利水電快報(bào) 2017年12期
    關(guān)鍵詞:基流徑流降雨

    (長(zhǎng)江水利委員會(huì) 長(zhǎng)江上游水文水資源勘測(cè)局,重慶 400020)

    試驗(yàn)與研究

    連續(xù)API產(chǎn)匯流模型的改進(jìn)及應(yīng)用

    曹磊陳川建

    (長(zhǎng)江水利委員會(huì) 長(zhǎng)江上游水文水資源勘測(cè)局,重慶 400020)

    結(jié)合長(zhǎng)江流域上游水文特性和實(shí)際預(yù)報(bào)工作,采用連續(xù)的日消退系數(shù)計(jì)算方法、連續(xù)的地下水分割方式、自動(dòng)計(jì)算降雨徑流相關(guān)曲線(P~Pa~R曲線)和瞬時(shí)單位線計(jì)算方法,并結(jié)合降雨徑流相關(guān)計(jì)算流程進(jìn)行產(chǎn)匯流計(jì)算和評(píng)定,得到一套較完整的連續(xù)API降雨徑流率定參數(shù)。通過該方法將API模型歸納成11個(gè)數(shù)字參數(shù),方便利用智能優(yōu)化算法進(jìn)行優(yōu)化率定計(jì)算,提高了API洪水預(yù)報(bào)模型的率定精度。提出的連續(xù)API模型的計(jì)算方法,為降雨徑流相關(guān)模型方案的建立、率定和實(shí)際應(yīng)用提供了一套新思路。

    產(chǎn)匯流計(jì)算;連續(xù)API模型;蓄滿產(chǎn)流;模型參數(shù)率定

    API模型區(qū)間產(chǎn)匯流模型作為一種洪水預(yù)報(bào)模型在大多數(shù)流域得到廣泛使用。該模型的參數(shù)率定涉及到場(chǎng)次洪水的選峰、地下水分割、凈流量計(jì)算、時(shí)段單位線率定等多方面,率定過程較復(fù)雜,主要靠人工結(jié)合計(jì)算機(jī)交互技術(shù),采用單場(chǎng)次洪水進(jìn)行交互分析和數(shù)據(jù)預(yù)處理才能完成,特別是利用人工退水曲線的方式來分割場(chǎng)次凈雨具有一定的隨意性。隨著地下水的分割[1]和還原[2-4]方法的使用,以及利用瞬時(shí)單位線取代時(shí)段單位線,逐漸與計(jì)算機(jī)技術(shù)結(jié)合后,提出了連續(xù)API模型并得到改進(jìn)[5-6],已在區(qū)間流計(jì)算方面得到推廣和應(yīng)用[7],取得了一定的效果。然而,目前還是通過選擇場(chǎng)次洪水人工交互分析得到連續(xù)API模型中退水曲線及場(chǎng)次洪水的分割,以及核心參數(shù)(P~Pa~R曲線)。同時(shí),日消退系數(shù)K值,仍然靠人為來確定,這些都約束了連續(xù)API模型的率定和自動(dòng)計(jì)算的應(yīng)用,沒有真正達(dá)到連續(xù)率定和連續(xù)計(jì)算的目的。本文以蓄滿產(chǎn)流原理為基礎(chǔ),結(jié)合地下水分割和分布水文模型中的植被指數(shù)等,結(jié)合計(jì)算機(jī)技術(shù)提出了一套真正意義上連續(xù)率定、連續(xù)計(jì)算的連續(xù)API模型,包括自動(dòng)進(jìn)行地下水還原計(jì)算、P~Pa~R曲線率定、計(jì)算瞬時(shí)單位線和每月K值的計(jì)算過程。

    1 模型介紹

    API模型是前期降雨指數(shù)(Antecedent Precipitation Index)模型的縮寫,其使用的表達(dá)式是傳統(tǒng)的降雨徑流相關(guān)圖,故又稱為降雨徑流經(jīng)驗(yàn)相關(guān)法。在20世紀(jì)40年代,R.K.林斯雷等提出了以季節(jié)、歷時(shí)、前期降雨指數(shù)API為參數(shù)的五變數(shù)降雨徑流相關(guān)圖,配合時(shí)段單位線技術(shù),成為一種有效的降雨洪水預(yù)報(bào)模型。1969年,西納(Sittner)等提出了模型地下徑流的建議,結(jié)合謝爾曼時(shí)段單位線即構(gòu)成了可模擬流域降雨凈雨過程的API模型。

    API模型屬于多輸入、單輸出靜態(tài)的系統(tǒng)數(shù)學(xué)模型。模型建立之初,主要用于一次洪水徑流量預(yù)報(bào)。在我國(guó),模型經(jīng)簡(jiǎn)化,形成了四變數(shù)關(guān)系線、三變數(shù)關(guān)系圖(P~Pa~R)或兩變數(shù)關(guān)系圖(P+Pa~R)。簡(jiǎn)化后的API模型便于方案的率定和使用。其降雨產(chǎn)流函數(shù)關(guān)系可表示為:R=f(P,Pa,T,季節(jié)) 、R=f(P,Pa,T) 、R=f(P,Pa)或R=f(P,P+Pa)。其中,北方流域使用較多P+Pa~R產(chǎn)流方案,而南方、特別是長(zhǎng)江流域多采用P~Pa~R三軸曲線。API模型的主要參數(shù)有:最大初損(IM)、前期降雨指數(shù)(Pa)、日消退系數(shù)(K)、降雨凈流相關(guān)曲線(P~Pa~R)和流域時(shí)段單位線(UH),計(jì)算時(shí)采用流域平均降雨來推算產(chǎn)匯流,參數(shù)的求解和含義可參考其他相應(yīng)文獻(xiàn)。

    在模型使用過程中,水文預(yù)報(bào)工作者對(duì)該模型進(jìn)行了不斷改進(jìn)[6],包括Pa計(jì)算方式的改進(jìn),并將分塊計(jì)算方式引入到了API模型中[6],考慮到降雨分布不均勻情況下的產(chǎn)流計(jì)算;將API的計(jì)算方式與計(jì)算機(jī)結(jié)合提出API模型的連續(xù)計(jì)算方法[7];將納西(Nash)瞬時(shí)單位線(IUH)引入模型計(jì)算中,利用調(diào)蓄滯時(shí)(Ku)和串聯(lián)水庫(kù)個(gè)數(shù)(N)兩個(gè)參數(shù)來確定流域瞬時(shí)單位線,再根據(jù)流域面積和計(jì)算時(shí)段長(zhǎng)也可將瞬時(shí)單位線轉(zhuǎn)化成時(shí)段單位線,使得單位線的率定更方便。

    2 模型的不足與改進(jìn)

    不管是傳統(tǒng)的API模型還是改進(jìn)后的模型,在進(jìn)行模型率定時(shí),主要是靠人工結(jié)合計(jì)算機(jī)交互功能,經(jīng)過場(chǎng)次樣本選取、基流分割、徑流計(jì)算、單位線推求、樣本評(píng)定5個(gè)主要的復(fù)雜步驟來完成。特別是地下水分割時(shí)具有一定的隨意性和不確定性(如圖1),使得模型的率定和實(shí)際使用過程中出現(xiàn)不一致的情況,這一定程度上制約了API模型的率定與使用,影響了模型計(jì)算精度。

    圖1 地下水不同分割示意

    針對(duì)API模型的不足做了改進(jìn),提出如下連續(xù)API模型的率定和計(jì)算方法。

    2.1 基流分割與還原

    2.1.1 常規(guī)分割方法的局限性

    在日常的API模型參數(shù)率定過程中,通常采用圖1中列出的方法來進(jìn)行基流分割,分割方式方法有一定的隨意性,計(jì)算的凈雨量和實(shí)際凈雨量有時(shí)差異較大。在API模型方案參數(shù)的實(shí)際使用過程中,降雨未結(jié)束洪峰未出現(xiàn)時(shí)(預(yù)報(bào)根據(jù)時(shí)間后的過程),無法確定直接徑流停止點(diǎn)(B點(diǎn))(平割法除外),同時(shí)也很難確定計(jì)算的直接徑流和地下水過程的總流量過程。不少使用者在直接徑流的基礎(chǔ)上加一個(gè)基流定值作為流域出口總出流,相當(dāng)于不管采用哪種分割方式確定的方案都當(dāng)成平割法進(jìn)行處理,這種處理在一定程度上降低了預(yù)報(bào)的精度,計(jì)算流程也不是太合理。

    2.1.2 數(shù)字濾波

    濾波起源于通信,它是將信號(hào)中特定波段頻率濾除得到有用信號(hào)的一種技術(shù)。內(nèi)森(Nathan)和麥克馬洪(McMahon)在1990年首次將數(shù)字濾波技術(shù)應(yīng)用到基流分割中,該方法已成為近年應(yīng)用最廣泛的基流分割方法[8]。分割方程如下:

    (1)

    Qb2=Q2-Qd2

    (2)

    式中,Qd2和Qd1分別為濾波得到的第2時(shí)刻和第1時(shí)刻的地表徑流,m3/s;Q2和Q1分別為第2時(shí)刻和第1時(shí)刻的徑流,m3/s;Qb2為第2時(shí)刻的基流,m3/s;α為濾波系數(shù)。

    這是一個(gè)遞推公式,基于較長(zhǎng)的系列資料進(jìn)行計(jì)算時(shí),特別是枯季開始時(shí),Qd1比較小,可設(shè)置一個(gè)初值,利用公式(1)計(jì)算得到Qd2,再利用公式(2)計(jì)算得到第2時(shí)刻的基流,以此類推計(jì)算得到整個(gè)洪水過程的基流。這種方法在水文中得到不斷的應(yīng)用和推廣,取得了較好的效果[9-11]。

    2.1.3 地下水還原計(jì)算[3]

    1969年美國(guó)天氣局提出API模型中地下水出流采用馬斯京根法線性水庫(kù)(x=0)演算,模型中假定入滲到地下水庫(kù)中的下滲量f與直接徑流Q呈線性關(guān)系,即:

    f=Z(Q-Qg)

    (3)

    地下水匯流用馬法演算,可以得到如下公式:

    Qg2=C0(f1+f2)+C2Qg1

    (4)

    將式(3)帶入式(4)得到如下地下水推算公式:

    (5)

    式中,Qg1為上時(shí)段初地下水流量,m3/s;Q1、Q2為時(shí)段始末的直接徑流量,m3/s;C0、C2為馬法演算參數(shù);Qg2為要計(jì)算的時(shí)段末地下水流量,m3/s;Z為地下水入流比例系數(shù)。

    該方法基于直接徑流進(jìn)行地下水還原計(jì)算,避開了地下水的分割難題,在實(shí)際預(yù)報(bào)過程中采用公式(5)即可推導(dǎo)得出。

    2.2 P~Pa~R參數(shù)計(jì)算改進(jìn)

    根據(jù)新安江三水源模型的蓄水容量曲線和降雨徑流關(guān)系的詮釋,蓄水容量曲線與降雨徑流具有對(duì)應(yīng)關(guān)系,故流域蓄水容量曲線的線形經(jīng)確定后降雨徑流關(guān)系曲線便隨之確定[3]。蓄水容量曲線的線性在濕潤(rùn)地區(qū)廣泛采用拋物線,因此,給定流域最大平均蓄水容量Wm(API模型里為IM)和拋物線的指數(shù)b,根據(jù)不同的前期蓄水量(在蓄滿產(chǎn)流中為W0,在API模型中為Pa),可對(duì)應(yīng)生成一簇P~Pa~R曲線或一條P+Pa~R曲線,把復(fù)雜的關(guān)系曲線族歸結(jié)于IM、b兩個(gè)參數(shù)。

    根據(jù)蓄滿產(chǎn)流的原理,只要有降雨,根據(jù)蓄水容量分布曲線總有部分流域產(chǎn)流。但在實(shí)際工作中發(fā)現(xiàn),在降雨較小的情況下,流域基本無產(chǎn)流或產(chǎn)流較小,連續(xù)API模型中設(shè)置了初始損失降雨P(guān)s,使得在Pa較小時(shí),曲線直接交Y軸與最大初損,Pa越大,初損逐步遞減;有的流域土壤厚度較大,很難蓄滿達(dá)到全流域產(chǎn)流的情況,在連續(xù)API模型中還設(shè)置了將整個(gè)曲線族按照曲線與Y軸交線進(jìn)行旋轉(zhuǎn)的參數(shù)β,得到最后相關(guān)線族。歸納起來,P~Pa~R曲線簇由流域最大缺水(IM)、蓄水容量曲線指數(shù)(b)、流域最大損失(Ps)和旋轉(zhuǎn)角度(β)4個(gè)參數(shù)可推導(dǎo)出來。

    2.3 K值的改進(jìn)

    早期的日消退系數(shù)K采用一個(gè)常值(通常取值0.85左右[4]),為了反映不同季節(jié)不同蒸散發(fā),后來多采用每月一個(gè)值,在一定程度上更趨于合理,但相當(dāng)于把原來單獨(dú)的一個(gè)參數(shù)無形中變成了12個(gè)參數(shù),增加了模型率定的難度。從K值反應(yīng)的蒸散發(fā)的含義出發(fā),力求算法簡(jiǎn)便。K.本旺(Keith Benven)[12-13]歸納了前人的研究成果,整理得到估算潛在蒸散發(fā)的季正弦曲線公式,公式如下:

    (6)

    另外,程根偉[14]等提出物候指數(shù)公式:

    (7)

    式中,Day為時(shí)間日序列的日數(shù),即從年初開始的天數(shù),d;A、B分別為與植物有關(guān)系的系數(shù);91為水文年與自然年的日差值。

    根據(jù)以上理論和API模型的日消退系數(shù)的含義,將該模型的日消退系數(shù)公式定義為:

    (8)

    圖2 日消退系數(shù)變化

    2.4 API模型中Pa數(shù)據(jù)計(jì)算模式

    最初的API模型中Pa值按照?qǐng)龃谓涤赀M(jìn)行計(jì)算。在連續(xù)API模型中,利用計(jì)算機(jī)連續(xù)計(jì)算時(shí),對(duì)于連續(xù)分散的小雨的情況,很難確定降雨開始時(shí)間及Pa值,這也會(huì)給預(yù)報(bào)計(jì)算帶來一定的誤差。連續(xù)API模型沿用原API模型的日模型計(jì)算結(jié)構(gòu),進(jìn)行每日滾動(dòng)計(jì)算日Pa值。每日切換一個(gè)Pa值對(duì)應(yīng)的一條P~Pa~R關(guān)系曲線,每時(shí)段采用累積降雨進(jìn)行查線的方法計(jì)算。

    2.5 改進(jìn)后的結(jié)構(gòu)

    圖3 連續(xù)API模型結(jié)構(gòu)示意

    該流程利用多個(gè)數(shù)據(jù)化參數(shù),概括了連續(xù)API模型結(jié)構(gòu)和率定、計(jì)算的順序,可利用遺傳算法等AI算法進(jìn)行優(yōu)化計(jì)算,得到較為精確的模型參數(shù)。

    表1 連續(xù)API模型優(yōu)化計(jì)算參數(shù)

    注:采用優(yōu)化計(jì)算方法自動(dòng)率定不包括濾波系數(shù)α。

    3 連續(xù)API模型的應(yīng)用

    連續(xù)API模型可以應(yīng)用人工逐步率定,也可以利用優(yōu)化算法進(jìn)行自動(dòng)優(yōu)化計(jì)算。

    3.1 人工交互率定

    人工交互率定可以分為兩個(gè)步驟完成。

    (1) 利用人工進(jìn)行交互率定時(shí),可先通過調(diào)整濾波系數(shù)(α),將整個(gè)過程分成地下徑流和直接徑流,再調(diào)整地下水馬法參數(shù)(Kmas)和地下水入流比例系數(shù)(Z),使得模擬的地下徑流跟分割的地下徑流過程吻合且具有一定的合理性。

    3.2 優(yōu)化計(jì)算

    連續(xù)API模型11個(gè)參數(shù)(自動(dòng)率定為10個(gè)參數(shù))都已數(shù)字化,可以直接利用如遺傳算法等優(yōu)化計(jì)算方法來進(jìn)行優(yōu)化率定。優(yōu)化計(jì)算的適應(yīng)度函數(shù)除了利用預(yù)報(bào)過程的確定性系數(shù)外,還可以結(jié)合計(jì)算過程的洪峰高度和洪峰時(shí)間誤差等系數(shù)統(tǒng)一進(jìn)行約束,達(dá)到洪峰高度和長(zhǎng)系列過程都能滿足要求的最優(yōu)參數(shù)。也可以先通過人工交互的方式調(diào)整好幾個(gè)不敏感的參數(shù)并進(jìn)行固定,其他剩下的參數(shù)再利用優(yōu)化計(jì)算法進(jìn)行優(yōu)化計(jì)算,最終得到合理的優(yōu)化模型參數(shù)。

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

    以長(zhǎng)江上游某流域?yàn)槔?,其集水面積1 441 km2,采用雨洪對(duì)應(yīng)較好的1990年洪水摘錄資料,計(jì)算時(shí)段長(zhǎng)Δt=1 h,利用遺傳算法,設(shè)置50個(gè)基因個(gè)體,對(duì)該連續(xù)API模型進(jìn)行優(yōu)化計(jì)算。采用完全確定性系數(shù)[15]控制和多洪峰(峰時(shí)和峰量)結(jié)合確定性系數(shù)控制(三者權(quán)值分別為20%,20%和60%)兩種方法分別進(jìn)行優(yōu)化計(jì)算,每種方式各計(jì)算10次,最后對(duì)計(jì)算結(jié)果進(jìn)行分析。

    分析得出,按照確定性系數(shù)控制條件每次計(jì)算得到的最低指標(biāo)都在乙級(jí)以上,按照多目標(biāo)權(quán)值控制計(jì)算模式計(jì)算結(jié)果適應(yīng)度函數(shù)平均值有所提高,但整個(gè)長(zhǎng)系列的確定性系數(shù)稍有減小,對(duì)于洪峰的高度控制較好。計(jì)算結(jié)果見圖4,模型參數(shù)如表1所示。

    圖4 連續(xù)API模型優(yōu)化計(jì)算結(jié)果

    利用該表中優(yōu)化后的參數(shù),逐年對(duì)1990~2012年23 a長(zhǎng)系列資料進(jìn)行檢驗(yàn),方案合格率達(dá)到乙級(jí)的年份所占百分比在75%以上,效果良好。

    4 展望與不足

    連續(xù)API模型采用11個(gè)數(shù)字參數(shù)很好地模擬了API模型的率定和實(shí)際使用計(jì)算,也為采用如遺傳算法等優(yōu)化計(jì)算方法進(jìn)行API模型參數(shù)的優(yōu)化計(jì)算提供了可能,一定程度簡(jiǎn)化了API模型率定的過程。通過控制優(yōu)化計(jì)算的控制條件,采用傳統(tǒng)的“確定性系數(shù)”結(jié)合洪峰高度和洪峰時(shí)間的適應(yīng)度函數(shù),提高并優(yōu)化了計(jì)算的精度。

    但由于連續(xù)API模型對(duì)一些參數(shù)進(jìn)行了概化,也存在以下一些不足之處。

    (1) 瞬時(shí)單位線由調(diào)蓄滯時(shí)(Ku)和串聯(lián)水庫(kù)個(gè)數(shù)(N)決定,在整個(gè)汛期過程采用的都是一條單位線,不能對(duì)單位線進(jìn)行分類和分組。

    (2) 公式(3)中,下滲量f和直接徑流成比列關(guān)系只是一個(gè)假設(shè),實(shí)際上f的量值既與土壤下滲特性有關(guān),還與包氣帶調(diào)蓄作用有關(guān)[3]。

    (3) 該方法只能提供三參數(shù)和兩參數(shù)的API模型計(jì)算,未能考慮四變數(shù)和五變數(shù)的模型。

    (4) 連續(xù)計(jì)算后,計(jì)算結(jié)果為某個(gè)流域的整個(gè)汛期產(chǎn)匯流過程,不是單場(chǎng)次洪水樣本,只能使用確定性系數(shù)法來評(píng)定。若使用該連續(xù)API模型還需要進(jìn)行場(chǎng)次評(píng)定,也可在整個(gè)長(zhǎng)系列計(jì)算完成后,選擇每場(chǎng)降雨和對(duì)應(yīng)的洪水過程作為樣本,設(shè)定樣本的開始時(shí)間和結(jié)束時(shí)間對(duì)每場(chǎng)洪水的預(yù)見期、變幅、凈雨、洪峰和峰現(xiàn)時(shí)間等進(jìn)行場(chǎng)次精度評(píng)定和統(tǒng)計(jì)。

    (5) 由于受到優(yōu)化計(jì)算算法的限制,每次得到的參數(shù)不一定為最優(yōu),計(jì)算完成后最好進(jìn)行合理性分析,并進(jìn)行多次計(jì)算,以選擇精度高并且又合理的參數(shù)為最后結(jié)果,在多次計(jì)算結(jié)果差別不大的情況下也可取平均值,需要具體問題具體分析。

    5 結(jié) 語

    隨著計(jì)算機(jī)技術(shù)在洪水預(yù)報(bào)中的不斷應(yīng)用,如何在原API模型中結(jié)合智能優(yōu)化算法實(shí)現(xiàn)自動(dòng)率定和參數(shù)自動(dòng)優(yōu)選,本文提出改進(jìn)后的“連續(xù)API模型”不失為一種有效、實(shí)用的方法。

    [1] 袁作新.流域水文模型[M].北京:中國(guó)水利電力出版社,1990.

    [2] 林三益.水文預(yù)報(bào)[M].北京:中國(guó)水利水電出版社(高等學(xué)校教材第二版),2001.

    [3] 張泉生.水文預(yù)報(bào)[M].南京:河海大學(xué)出版社, 2000.

    [4] 包為民.水文預(yù)報(bào)(第三版)[M].北京:中國(guó)水利水電出版社,2006.

    [5] 劉國(guó)富.連續(xù)API水文模型的研究及應(yīng)用[J].大壩與安全,2008(2):4-8.

    [6] 張恭肅,王成明.對(duì)API模型的改進(jìn)[J].水文.1996(4):20-25.

    [7] 鐘小燕,文磊,余鐘波.連續(xù)API模型子沂河臨沂站徑流預(yù)報(bào)中的應(yīng)用[J].人民長(zhǎng)江,2017,48(13):26-30.

    [8] 熊立華,郭生練.采用非線性水庫(kù)假設(shè)的基流分割方法及應(yīng)用[J].武漢大學(xué)學(xué)報(bào)(工學(xué)版),2005,38(1):27-29.

    [9] 黃國(guó)如.流量過程線的自動(dòng)分割方法探討[J].灌溉排水學(xué)報(bào),2007,26(1):73-76.

    [10]徐磊磊.水文過程的基流分割方法研究進(jìn)展[J].應(yīng)用生態(tài)學(xué)報(bào),2011,22(11):3073-3080.

    [11]陳文艷.數(shù)字濾波法分割基流的論證[J].水力發(fā)電,2014,40(2):37-40.

    [12]Keith Beven.Rainfall-Runoff Modelling Lancaster University[M]. UK:John Wiley & Sons, Ltd, 2012.

    [13]張行南,劉建芬,夏達(dá)忠.降雨-徑流模擬[M].北京:中國(guó)水利水電出版社,2016.

    [14]程根偉,舒棟材.水文預(yù)報(bào)的理論與數(shù)學(xué)模型[M].北京:中國(guó)水利水電出版社,2006.

    [15]GBT 22482-2008水文情報(bào)預(yù)報(bào)規(guī)范.

    2017-10-10

    曹 磊,男,長(zhǎng)江水利委員會(huì)長(zhǎng)江上游水文水資源勘測(cè)局,工程師.

    1006-0081(2017)12-0031-05

    P456.8

    A

    (編輯:朱曉紅)

    猜你喜歡
    基流徑流降雨
    美國(guó)長(zhǎng)期氣象干旱對(duì)河道基流的影響研究
    滄州市2016年“7.19~7.22”與“8.24~8.25”降雨對(duì)比研究
    基于改進(jìn)Tennant法的洋河流域生態(tài)基流估算研究
    紅黏土降雨入滲的定量分析
    Topmodel在布哈河流域徑流模擬中的應(yīng)用
    不同生態(tài)基流對(duì)水利樞紐發(fā)電損益的影響
    探秘“大徑流”
    攻克“大徑流”
    南方降雨不斷主因厄爾尼諾
    增加河流生態(tài)基流的山間河谷型調(diào)蓄模式研究
    看非洲黑人一级黄片| 欧美亚洲日本最大视频资源| 两个人免费观看高清视频| 久久精品人人爽人人爽视色| 欧美激情 高清一区二区三区| 欧美xxⅹ黑人| 一区二区av电影网| 啦啦啦中文免费视频观看日本| 97精品久久久久久久久久精品| av线在线观看网站| 亚洲国产精品999| 制服丝袜香蕉在线| av播播在线观看一区| 免费久久久久久久精品成人欧美视频 | 男人爽女人下面视频在线观看| 久热久热在线精品观看| 九草在线视频观看| 亚洲国产av新网站| av国产久精品久网站免费入址| 香蕉精品网在线| 宅男免费午夜| 国产精品一国产av| 在线观看美女被高潮喷水网站| 少妇 在线观看| 亚洲三级黄色毛片| 天天躁夜夜躁狠狠久久av| 免费不卡的大黄色大毛片视频在线观看| av卡一久久| 在线亚洲精品国产二区图片欧美| 成年人免费黄色播放视频| 两性夫妻黄色片 | 免费看av在线观看网站| 欧美人与善性xxx| 欧美日韩成人在线一区二区| 亚洲精品久久午夜乱码| 狠狠精品人妻久久久久久综合| 捣出白浆h1v1| 日本wwww免费看| 亚洲,一卡二卡三卡| 国产一级毛片在线| 久久久国产精品麻豆| 中文字幕人妻熟女乱码| 少妇猛男粗大的猛烈进出视频| 飞空精品影院首页| 国产精品一区二区在线观看99| 99热国产这里只有精品6| 观看美女的网站| 国产亚洲一区二区精品| 99国产精品免费福利视频| 亚洲精品乱码久久久久久按摩| 2018国产大陆天天弄谢| 人妻系列 视频| 国产精品人妻久久久久久| a级毛片在线看网站| 午夜久久久在线观看| 中文乱码字字幕精品一区二区三区| videos熟女内射| 九九在线视频观看精品| 一本色道久久久久久精品综合| 精品久久久精品久久久| 免费人成在线观看视频色| 男女免费视频国产| 亚洲精品国产色婷婷电影| 欧美xxxx性猛交bbbb| 午夜视频国产福利| 亚洲精品乱码久久久久久按摩| 国产视频首页在线观看| 香蕉精品网在线| 在线天堂中文资源库| 欧美日韩精品成人综合77777| av在线播放精品| 飞空精品影院首页| 欧美最新免费一区二区三区| 欧美精品国产亚洲| 韩国av在线不卡| a级片在线免费高清观看视频| 久久99精品国语久久久| 校园人妻丝袜中文字幕| 极品少妇高潮喷水抽搐| 男女免费视频国产| 欧美国产精品一级二级三级| 婷婷色av中文字幕| 免费高清在线观看日韩| 国产精品免费大片| 韩国av在线不卡| 亚洲性久久影院| 日本av免费视频播放| 最近中文字幕高清免费大全6| 香蕉精品网在线| 女人久久www免费人成看片| 一区二区av电影网| 水蜜桃什么品种好| www.av在线官网国产| 哪个播放器可以免费观看大片| 久久97久久精品| 狠狠精品人妻久久久久久综合| 国产精品一区二区在线观看99| 99热网站在线观看| 亚洲国产精品成人久久小说| 在线精品无人区一区二区三| 人人妻人人澡人人看| 欧美日韩国产mv在线观看视频| 国产亚洲精品第一综合不卡 | 三级国产精品片| 高清视频免费观看一区二区| 亚洲av.av天堂| 春色校园在线视频观看| 亚洲精华国产精华液的使用体验| 国产欧美日韩综合在线一区二区| 大香蕉97超碰在线| 日韩成人伦理影院| 亚洲国产av新网站| 狂野欧美激情性bbbbbb| 人妻一区二区av| 国产熟女午夜一区二区三区| 午夜福利视频精品| 午夜福利乱码中文字幕| 两性夫妻黄色片 | 色5月婷婷丁香| 少妇被粗大猛烈的视频| 热99国产精品久久久久久7| 久热这里只有精品99| freevideosex欧美| 亚洲高清免费不卡视频| 成人综合一区亚洲| 免费人妻精品一区二区三区视频| 嫩草影院入口| 考比视频在线观看| 视频区图区小说| 日韩免费高清中文字幕av| 男人添女人高潮全过程视频| 亚洲色图 男人天堂 中文字幕 | 午夜福利,免费看| 一区在线观看完整版| 亚洲婷婷狠狠爱综合网| h视频一区二区三区| 久久99精品国语久久久| 国产精品久久久久久av不卡| 老司机影院毛片| 亚洲精品乱久久久久久| 51国产日韩欧美| 欧美激情国产日韩精品一区| 国产精品国产三级国产av玫瑰| 精品一区二区三区视频在线| 精品人妻熟女毛片av久久网站| 亚洲欧洲国产日韩| 最近2019中文字幕mv第一页| 国产成人精品在线电影| 国产精品免费大片| 飞空精品影院首页| 国产女主播在线喷水免费视频网站| av天堂久久9| 亚洲国产欧美日韩在线播放| 午夜福利网站1000一区二区三区| av国产精品久久久久影院| 极品少妇高潮喷水抽搐| 国产精品成人在线| 高清欧美精品videossex| 亚洲人成网站在线观看播放| www.色视频.com| 欧美 亚洲 国产 日韩一| 一级爰片在线观看| av免费观看日本| 日韩欧美精品免费久久| 亚洲欧美成人精品一区二区| 亚洲欧洲国产日韩| 丝袜脚勾引网站| 高清视频免费观看一区二区| 777米奇影视久久| 国产乱来视频区| 久久久国产精品麻豆| 国产亚洲欧美精品永久| 18禁观看日本| 日韩av在线免费看完整版不卡| 啦啦啦中文免费视频观看日本| 午夜福利视频在线观看免费| 成人黄色视频免费在线看| 丰满少妇做爰视频| 亚洲在久久综合| 99热网站在线观看| 欧美成人午夜精品| 寂寞人妻少妇视频99o| 制服人妻中文乱码| 国产色婷婷99| 亚洲欧美成人综合另类久久久| www日本在线高清视频| 欧美激情 高清一区二区三区| 亚洲国产av新网站| 激情视频va一区二区三区| 免费av不卡在线播放| 亚洲欧洲日产国产| 亚洲国产欧美日韩在线播放| 男女边摸边吃奶| 高清不卡的av网站| 国产成人午夜福利电影在线观看| 亚洲成国产人片在线观看| 欧美xxxx性猛交bbbb| 亚洲中文av在线| 午夜激情久久久久久久| 日韩制服丝袜自拍偷拍| 国产黄频视频在线观看| 亚洲国产色片| 美女脱内裤让男人舔精品视频| 亚洲欧洲精品一区二区精品久久久 | 国产成人免费观看mmmm| 国产无遮挡羞羞视频在线观看| 人体艺术视频欧美日本| 国产精品一区二区在线不卡| 免费少妇av软件| 欧美xxⅹ黑人| 国产精品国产三级国产av玫瑰| 中文欧美无线码| 90打野战视频偷拍视频| 日本黄大片高清| 国产欧美日韩综合在线一区二区| 激情五月婷婷亚洲| 有码 亚洲区| 国产成人免费无遮挡视频| 在线观看免费高清a一片| 性色av一级| 熟女av电影| 亚洲美女黄色视频免费看| 在线观看美女被高潮喷水网站| 免费看av在线观看网站| 国产在线视频一区二区| 一区在线观看完整版| 国精品久久久久久国模美| 久久久久久久大尺度免费视频| 精品少妇黑人巨大在线播放| 街头女战士在线观看网站| 青春草亚洲视频在线观看| 久久久精品区二区三区| 99国产精品免费福利视频| 哪个播放器可以免费观看大片| 啦啦啦视频在线资源免费观看| 成人综合一区亚洲| 我的女老师完整版在线观看| 成年av动漫网址| 亚洲精品自拍成人| 曰老女人黄片| av.在线天堂| av国产精品久久久久影院| 在线观看免费视频网站a站| 亚洲美女视频黄频| 三级国产精品片| 日本-黄色视频高清免费观看| 青春草亚洲视频在线观看| 少妇精品久久久久久久| 日本wwww免费看| 国产精品一二三区在线看| 你懂的网址亚洲精品在线观看| 综合色丁香网| 另类亚洲欧美激情| 成人国语在线视频| 久久久久久久大尺度免费视频| 母亲3免费完整高清在线观看 | 久久午夜综合久久蜜桃| 久久婷婷青草| 色婷婷av一区二区三区视频| 国产欧美日韩一区二区三区在线| 国产精品一国产av| 久久99热6这里只有精品| 免费女性裸体啪啪无遮挡网站| 亚洲,一卡二卡三卡| 久热这里只有精品99| 亚洲av男天堂| 久热久热在线精品观看| 1024视频免费在线观看| 99久久精品国产国产毛片| 欧美激情极品国产一区二区三区 | 五月开心婷婷网| 国产亚洲最大av| 亚洲第一av免费看| 一本—道久久a久久精品蜜桃钙片| 免费看av在线观看网站| 欧美日韩成人在线一区二区| 国产黄色视频一区二区在线观看| 亚洲一码二码三码区别大吗| 久久精品国产综合久久久 | 丰满少妇做爰视频| 90打野战视频偷拍视频| 欧美日韩视频高清一区二区三区二| 国产xxxxx性猛交| 极品人妻少妇av视频| 国产无遮挡羞羞视频在线观看| 久久久久久久久久久久大奶| 欧美日韩国产mv在线观看视频| 久热这里只有精品99| 国产一区二区在线观看av| 97精品久久久久久久久久精品| 熟女电影av网| 女性被躁到高潮视频| 18+在线观看网站| 免费不卡的大黄色大毛片视频在线观看| 丝袜在线中文字幕| 亚洲精品日韩在线中文字幕| 在现免费观看毛片| 最近中文字幕2019免费版| 一级a做视频免费观看| 午夜精品国产一区二区电影| 国产永久视频网站| 99久国产av精品国产电影| 亚洲国产日韩一区二区| 欧美性感艳星| 久久国产精品男人的天堂亚洲 | 一级黄片播放器| 亚洲成色77777| 成人18禁高潮啪啪吃奶动态图| 亚洲国产毛片av蜜桃av| 日韩中字成人| 欧美成人午夜精品| 欧美日韩综合久久久久久| 黄片播放在线免费| 在线观看免费视频网站a站| 久久久久视频综合| 国产精品蜜桃在线观看| 成人手机av| 国语对白做爰xxxⅹ性视频网站| 亚洲,欧美,日韩| 久久人人爽人人爽人人片va| 高清在线视频一区二区三区| 日韩精品有码人妻一区| 亚洲成人一二三区av| 边亲边吃奶的免费视频| 国产在视频线精品| 久久久久国产网址| 亚洲情色 制服丝袜| 精品国产一区二区三区四区第35| 99热网站在线观看| 久久久久久久国产电影| 一区在线观看完整版| 少妇 在线观看| 精品熟女少妇av免费看| 午夜福利视频精品| 亚洲精品美女久久av网站| 黑人欧美特级aaaaaa片| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产精品三级大全| 久久精品aⅴ一区二区三区四区 | 高清不卡的av网站| 国产白丝娇喘喷水9色精品| 久久久久久久精品精品| 久久精品国产亚洲av涩爱| 国产精品一区二区在线不卡| 国产伦理片在线播放av一区| 丰满迷人的少妇在线观看| 少妇的逼水好多| 国产精品欧美亚洲77777| 国产免费又黄又爽又色| 少妇的逼好多水| 在线观看一区二区三区激情| 最黄视频免费看| 亚洲精品日本国产第一区| 亚洲第一区二区三区不卡| 一区二区三区乱码不卡18| 老司机影院毛片| 黄色毛片三级朝国网站| 女人久久www免费人成看片| 最近中文字幕高清免费大全6| 国产色爽女视频免费观看| 亚洲美女黄色视频免费看| 久久精品夜色国产| 肉色欧美久久久久久久蜜桃| 亚洲精品久久午夜乱码| www.av在线官网国产| 亚洲色图综合在线观看| 不卡视频在线观看欧美| 少妇猛男粗大的猛烈进出视频| 女性生殖器流出的白浆| 午夜老司机福利剧场| 国产av精品麻豆| 日韩欧美精品免费久久| 亚洲欧美一区二区三区黑人 | 亚洲三级黄色毛片| 欧美少妇被猛烈插入视频| av视频免费观看在线观看| www.熟女人妻精品国产 | 卡戴珊不雅视频在线播放| 久久久久国产网址| 精品卡一卡二卡四卡免费| 日本欧美视频一区| 美女主播在线视频| 久久午夜福利片| 一个人免费看片子| 如日韩欧美国产精品一区二区三区| 精品卡一卡二卡四卡免费| 亚洲性久久影院| 亚洲精品乱久久久久久| 亚洲av综合色区一区| 天堂中文最新版在线下载| 久久午夜综合久久蜜桃| 自线自在国产av| 欧美亚洲 丝袜 人妻 在线| 咕卡用的链子| 免费在线观看完整版高清| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲av综合色区一区| av网站免费在线观看视频| 亚洲婷婷狠狠爱综合网| 免费在线观看完整版高清| 久久久国产一区二区| 国产在线一区二区三区精| 91aial.com中文字幕在线观看| 国产老妇伦熟女老妇高清| 看十八女毛片水多多多| 国产亚洲最大av| 免费少妇av软件| 菩萨蛮人人尽说江南好唐韦庄| 建设人人有责人人尽责人人享有的| 免费人妻精品一区二区三区视频| 18禁动态无遮挡网站| 26uuu在线亚洲综合色| 99久久人妻综合| 亚洲国产精品专区欧美| 日韩视频在线欧美| 久久ye,这里只有精品| 亚洲,欧美,日韩| 亚洲,一卡二卡三卡| 久久女婷五月综合色啪小说| 在现免费观看毛片| 这个男人来自地球电影免费观看 | 国产老妇伦熟女老妇高清| 男女啪啪激烈高潮av片| 国产一级毛片在线| 欧美bdsm另类| 国产日韩欧美亚洲二区| 2018国产大陆天天弄谢| 亚洲,欧美精品.| av国产精品久久久久影院| 亚洲成色77777| 欧美人与性动交α欧美软件 | 如何舔出高潮| 桃花免费在线播放| 一区二区日韩欧美中文字幕 | 午夜福利乱码中文字幕| 青青草视频在线视频观看| 黄色一级大片看看| 国语对白做爰xxxⅹ性视频网站| 国产免费福利视频在线观看| 少妇熟女欧美另类| 欧美bdsm另类| 亚洲欧美色中文字幕在线| 男女边吃奶边做爰视频| 不卡视频在线观看欧美| av电影中文网址| 中文天堂在线官网| 亚洲精品中文字幕在线视频| 日本wwww免费看| 在线观看一区二区三区激情| 搡老乐熟女国产| 日本黄色日本黄色录像| 欧美日韩国产mv在线观看视频| 大话2 男鬼变身卡| 天堂俺去俺来也www色官网| av免费观看日本| 国产乱人偷精品视频| 国产亚洲欧美精品永久| 欧美+日韩+精品| 亚洲情色 制服丝袜| 亚洲欧美成人精品一区二区| 欧美国产精品一级二级三级| 最近的中文字幕免费完整| 精品一区二区三卡| 亚洲精品,欧美精品| 日韩一区二区三区影片| 一区二区三区精品91| 韩国高清视频一区二区三区| 丝袜人妻中文字幕| 97精品久久久久久久久久精品| 日韩一本色道免费dvd| 国产黄色视频一区二区在线观看| 国产精品 国内视频| 狂野欧美激情性bbbbbb| 亚洲av中文av极速乱| 免费黄色在线免费观看| 这个男人来自地球电影免费观看 | 久久久亚洲精品成人影院| 久久久久网色| 制服丝袜香蕉在线| 五月天丁香电影| 亚洲第一区二区三区不卡| 欧美人与性动交α欧美软件 | 亚洲人与动物交配视频| 最近最新中文字幕大全免费视频 | 国产高清三级在线| 日本wwww免费看| 熟妇人妻不卡中文字幕| 久久这里只有精品19| 美女中出高潮动态图| 成人影院久久| 女人久久www免费人成看片| 免费日韩欧美在线观看| 黄片无遮挡物在线观看| 日本91视频免费播放| 伦理电影免费视频| 国产成人精品在线电影| 最近2019中文字幕mv第一页| 国产极品粉嫩免费观看在线| 国产xxxxx性猛交| 欧美日韩视频高清一区二区三区二| 黑丝袜美女国产一区| 国产精品久久久久久精品电影小说| 国产一区二区三区av在线| 国产日韩欧美亚洲二区| 制服人妻中文乱码| 大陆偷拍与自拍| 亚洲精品第二区| 最新中文字幕久久久久| 亚洲精品久久成人aⅴ小说| 免费观看性生交大片5| 国产有黄有色有爽视频| 国产日韩欧美视频二区| 亚洲国产精品专区欧美| 蜜桃在线观看..| 少妇猛男粗大的猛烈进出视频| 永久免费av网站大全| 美女视频免费永久观看网站| 国产欧美另类精品又又久久亚洲欧美| 久久av网站| av电影中文网址| 免费高清在线观看视频在线观看| 涩涩av久久男人的天堂| 青春草亚洲视频在线观看| 国产男女超爽视频在线观看| 少妇熟女欧美另类| 日日撸夜夜添| 亚洲国产精品999| av电影中文网址| 在线观看免费日韩欧美大片| 国产深夜福利视频在线观看| www.熟女人妻精品国产 | 曰老女人黄片| 久久精品国产自在天天线| 日本欧美视频一区| 啦啦啦中文免费视频观看日本| 国产成人精品久久久久久| 亚洲国产精品成人久久小说| 日韩av在线免费看完整版不卡| 久久久久国产精品人妻一区二区| 欧美日本中文国产一区发布| 免费播放大片免费观看视频在线观看| 99热全是精品| 日韩熟女老妇一区二区性免费视频| 日韩成人av中文字幕在线观看| 久久久久久久久久成人| 五月玫瑰六月丁香| 久久精品久久精品一区二区三区| 精品亚洲成a人片在线观看| 观看av在线不卡| 中文字幕免费在线视频6| 日韩精品有码人妻一区| 免费黄网站久久成人精品| 久久久久精品久久久久真实原创| 久久久久久久国产电影| 国产一区二区在线观看av| 老女人水多毛片| 久久人人爽人人爽人人片va| 欧美变态另类bdsm刘玥| 亚洲av电影在线观看一区二区三区| 欧美+日韩+精品| 中文字幕av电影在线播放| 国产白丝娇喘喷水9色精品| 国语对白做爰xxxⅹ性视频网站| 日本黄大片高清| 免费黄色在线免费观看| 国产一区二区三区av在线| 九九爱精品视频在线观看| 一本色道久久久久久精品综合| 免费人成在线观看视频色| 国产精品一区二区在线不卡| 久久鲁丝午夜福利片| 国精品久久久久久国模美| 亚洲经典国产精华液单| www日本在线高清视频| 欧美bdsm另类| 午夜免费男女啪啪视频观看| 最后的刺客免费高清国语| 中文欧美无线码| 男人添女人高潮全过程视频| 国产成人a∨麻豆精品| 国产极品天堂在线| 国产色爽女视频免费观看| av国产精品久久久久影院| 天天躁夜夜躁狠狠久久av| av黄色大香蕉| 亚洲精品久久午夜乱码| 国产精品 国内视频| av在线app专区| 天美传媒精品一区二区| 美女大奶头黄色视频| 国产日韩欧美亚洲二区| 久久久久久久精品精品| a级毛片在线看网站| 精品亚洲成国产av| 亚洲 欧美一区二区三区| 亚洲精品自拍成人| 丝袜喷水一区| 日日摸夜夜添夜夜爱| freevideosex欧美| 天天影视国产精品| 天美传媒精品一区二区| 久久免费观看电影| 国产亚洲av片在线观看秒播厂| 午夜免费鲁丝| 久久久亚洲精品成人影院| 在线观看人妻少妇| 国精品久久久久久国模美| 亚洲精品久久久久久婷婷小说| 天天影视国产精品| 99久久综合免费| 中文欧美无线码| 成人午夜精彩视频在线观看| 中文字幕精品免费在线观看视频 | 欧美丝袜亚洲另类| 一级a做视频免费观看|