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

    時(shí)間步長(zhǎng)和產(chǎn)匯流方法對(duì)小流域次洪模擬的影響

    2019-09-10 08:42:45趙娜趙雪花祝雪萍
    人民長(zhǎng)江 2019年1期
    關(guān)鍵詞:洪量產(chǎn)流步長(zhǎng)

    趙娜 趙雪花 祝雪萍

    摘要:為研究半干旱半濕潤(rùn)流域復(fù)雜產(chǎn)流機(jī)制下不同時(shí)間步長(zhǎng)和產(chǎn)匯流方法對(duì)次洪模擬的影響,選取北張店流域?yàn)檠芯繉?duì)象,對(duì)1993~2005年具有代表性的12場(chǎng)次洪降雨徑流資料,利用HEC-HMS模型分別以10,15,20 min為時(shí)間步長(zhǎng),采用2種方案進(jìn)行模擬計(jì)算并對(duì)比分析結(jié)果。方案1為初損穩(wěn)滲法+Clark單位線(xiàn)法;方案2為格林-安普特法+SCS單位線(xiàn)法;河道匯流演算均為馬斯京根法,基流計(jì)算均為指數(shù)消退法,結(jié)果表明:時(shí)間步長(zhǎng)的選取對(duì)次洪模擬結(jié)果影響較大,兩種方案均在15 min時(shí)間步長(zhǎng)下模擬精度最高,表明時(shí)間步長(zhǎng)的選取并非是越短越精確。時(shí)間步長(zhǎng)為15 min時(shí)方案1的確定性系數(shù)(DC)均值為0.813,方案2的DC均值為0.764,且方案1在洪峰、洪量和峰現(xiàn)時(shí)差模擬方面均優(yōu)于方案2。

    關(guān)?鍵?詞:時(shí)間步長(zhǎng); 產(chǎn)匯流; 次洪模擬; HEC-HMS模型; 半干旱半濕潤(rùn)流域; 北張店流域; 山西省

    中圖法分類(lèi)號(hào): P33?文獻(xiàn)標(biāo)志碼: ADOI:10.16232/j.cnki.1001-4179.2019.01.012

    我國(guó)北方多為半干旱半濕潤(rùn)地區(qū),蓄滿(mǎn)產(chǎn)流與超滲產(chǎn)流并存[1]。由于受季風(fēng)氣候的影響,雨強(qiáng)大、歷時(shí)短,一旦遭遇洪水,具有來(lái)勢(shì)猛、水量集中、陡漲陡落等特點(diǎn),因此半干旱半濕潤(rùn)地區(qū)洪水預(yù)報(bào)一直都是研究的重點(diǎn)和難點(diǎn)。根據(jù)研究區(qū)特性選擇合適的水文計(jì)算方法和時(shí)間尺度已受到了廣泛的關(guān)注。DA Zema等[2]應(yīng)用HEC-HMS模型分析了不同產(chǎn)流方法對(duì)模擬結(jié)果的影響,結(jié)果表明產(chǎn)流方法對(duì)模擬洪峰大小和峰現(xiàn)時(shí)差影響較大;Hessel R[3]應(yīng)用LISEM模型探討了時(shí)間步長(zhǎng)對(duì)模擬結(jié)果的影響,結(jié)果表明隨著時(shí)間步長(zhǎng)的增大,流域降雨強(qiáng)度平均化,模擬洪峰流量減小;Ren Liliang等[4]將不同水文模型應(yīng)用在半干旱半濕潤(rùn)地區(qū),結(jié)果表明時(shí)間尺度對(duì)模型模擬的影響較大;林峰等[5]的研究成果表明模擬時(shí)間步長(zhǎng)增加,模擬峰現(xiàn)時(shí)間延遲、峰值和洪量均減小;Jin Hua等[6]將HEC-HMS模型應(yīng)用于絳河流域,結(jié)果表明該模型在該流域適用性較好,但峰現(xiàn)時(shí)差較大,還需進(jìn)一步對(duì)產(chǎn)匯流方法和時(shí)間步長(zhǎng)進(jìn)行分析以提高模擬精度。

    北張店流域地處黃河流域,屬于典型的半干旱半濕潤(rùn)地區(qū),產(chǎn)匯流機(jī)制十分復(fù)雜。本文以北張店流域?yàn)檠芯繀^(qū)域,基于1993~2005年具有代表性的12場(chǎng)次洪降雨徑流資料,應(yīng)用HEC-HMS分布式模型模擬研究區(qū)次降雨洪水過(guò)程,探討了不同水文計(jì)算方法在多時(shí)間步長(zhǎng)下對(duì)次洪過(guò)程模擬精度的影響。研究成果對(duì)提高半干旱半濕潤(rùn)地區(qū)的洪水預(yù)報(bào)精度有重要參考價(jià)值。

    1?研究區(qū)域與數(shù)據(jù)處理

    1.1?研究區(qū)域

    北張店流域位于山西省長(zhǎng)治市中西部,屯留縣和沁縣的交界處,地理位置為112°26′E~112°37′E,36°15′N(xiāo)~36°32′N(xiāo),流域總面積270 km2,土壤為褐土。區(qū)域內(nèi)地勢(shì)北西南三面高、東部及中部低,水流由北西南向中部匯集后向東出境,最高海拔1 544 m,最低海拔976 m,平均海拔1 276.5 m。研究區(qū)屬于暖溫帶季風(fēng)區(qū),具有典型的半干旱半濕潤(rùn)地區(qū)氣候特征。氣候特點(diǎn)是四季分明、冬長(zhǎng)夏短、春季略長(zhǎng)于秋季、日照充足。年平均氣溫9.0℃,最高溫度37.4℃,最低溫度為-22.0℃;最大凍深70 cm,年無(wú)霜期155 d,年均日照時(shí)數(shù)2 418 h,多年平均蒸發(fā)量1 775.8 mm。北張店流域降水量年際變化較大,全流域多年平均降水量為564 mm,最大年降水量872.2 mm,最小年降水量299.1 mm,流域內(nèi)降雨屬暴雨型,強(qiáng)度大、歷時(shí)短,最大暴雨多集中在7~8月。

    1.2?研究區(qū)數(shù)據(jù)預(yù)處理

    根據(jù)研究區(qū)的地理位置,在地理空間數(shù)據(jù)云(http://www.gscloud.cn)下載GDEM DEM 30 m分辨率數(shù)字高程數(shù)據(jù),利用ArcGIS平臺(tái)提取北張店流域的DEM數(shù)據(jù)?;谘芯繀^(qū)的DEM數(shù)據(jù),運(yùn)用內(nèi)嵌于ArcGIS軟件中的HEC-GeoHMS拓展模塊,通過(guò)填洼、計(jì)算流向、匯流等,提取出水系并將流域劃分為15個(gè)子流域(如圖1所示),最終生成HEC-HMS模型可以使用的數(shù)據(jù)。

    1.3?氣象數(shù)據(jù)庫(kù)的建立

    利用可視化數(shù)據(jù)存儲(chǔ)系統(tǒng)HEC-DSSVue建立研究區(qū)雨量數(shù)據(jù)庫(kù)和流量數(shù)據(jù)庫(kù)。將流域內(nèi)12個(gè)雨量站1993~2005年洪水期降雨徑流數(shù)據(jù)輸入數(shù)據(jù)庫(kù),并利用HEC-DSSVue的數(shù)學(xué)函數(shù)功能,將時(shí)間序列數(shù)據(jù)按10,15,20 min步長(zhǎng)分別插值,最終得到每場(chǎng)降雨雨量和每場(chǎng)洪水實(shí)測(cè)流量數(shù)據(jù)序列。

    2?模型簡(jiǎn)介和方案選取

    2.1HEC-HMS模型

    HEC-HMS水文模型系統(tǒng)是美國(guó)陸軍工程師兵團(tuán)(US Army Corps of Engineer,USACE)水文工程中心(Hydrologic Engineering Center,HEC)開(kāi)發(fā)的新一代軟件產(chǎn)品,被廣泛應(yīng)用于模擬樹(shù)狀流域中降雨-徑流的過(guò)程[7-8]。該水文模型基本原理是在考慮降雨空間分布不均勻性和下墊面條件非均勻性的基礎(chǔ)上,將研究流域劃分成若干子流域(或網(wǎng)格單元),在每一個(gè)子流域內(nèi),根據(jù)其下墊面條件選定參數(shù)進(jìn)行降雨徑流過(guò)程模擬計(jì)算,然后將各子流域計(jì)算結(jié)果演算到流域出口處的流域降雨徑流過(guò)程即為HEC-HMS模型的降雨徑流概化過(guò)程。

    HEC-HMS水文模型主要包括氣象模塊、流域模塊、模擬控制模塊和時(shí)間序列模塊。流域模塊中主要包括凈雨計(jì)算、直接徑流、河道洪水演算和基流4部分。每一部分都可以根據(jù)研究區(qū)情況來(lái)選取不同的計(jì)算方法,通過(guò)設(shè)置各個(gè)子流域的參數(shù)來(lái)計(jì)算。

    2.2?方案的選取

    本文產(chǎn)流計(jì)算模塊選取了國(guó)內(nèi)外較常用的初損常數(shù)法[9-10]和理論基礎(chǔ)充分、模擬結(jié)果較好的格林-安普特法[11]。初損常數(shù)方法只有一個(gè)參數(shù)即土壤最大穩(wěn)滲能力?fc,需要率定。而初始損失量Ia?要取決于研究區(qū)的地形地貌、土壤類(lèi)型、植被條件以及土地使用類(lèi)型等。本文按照初始損失量為總降雨量的10%~20%給出每個(gè)子流域初始參數(shù)值[15]。格林-安普特法是一種概念性下滲方法,它通過(guò)求解降雨徑流過(guò)程中的降雨損失量,來(lái)反求流域產(chǎn)流量。

    直接徑流計(jì)算模塊選取了Clark單位線(xiàn)法和SCS單位線(xiàn)法。河道匯流演算大多數(shù)采用馬斯京根法和運(yùn)動(dòng)波法[12-13]。其中馬斯京根法以所用到的參數(shù)少、方便率定,在國(guó)內(nèi)外水文模擬中廣泛應(yīng)用,效果較好[14],2種方案均選用其進(jìn)行河道匯流演算。基流計(jì)算方法有常數(shù)月變化法、指數(shù)消退法及線(xiàn)性水庫(kù)法。由于指數(shù)消退法很好地解釋了流域中的水量在無(wú)雨期間向流域河流進(jìn)行水量補(bǔ)給的過(guò)程,所以2種方案均選用指數(shù)消退法。模型組合方案見(jiàn)表1。

    3?結(jié)果與分析

    模擬選用8場(chǎng)降雨徑流資料進(jìn)行參數(shù)率定、4場(chǎng)實(shí)測(cè)洪水資料用于模型檢驗(yàn)。所選取的12場(chǎng)次降雨洪水的洪峰流量最小值是11.7 m3/s,最大值是158 m3/s。采用權(quán)重因子法對(duì)降水?dāng)?shù)據(jù)進(jìn)行插值處理,分別以10,15,20 min為模擬時(shí)間步長(zhǎng)。根據(jù)《水文情報(bào)預(yù)報(bào)規(guī)范》(GB/T22482-2008)[16]對(duì)洪峰流量相對(duì)誤差?REp、洪量相對(duì)誤差REv、峰現(xiàn)時(shí)間誤差△T和確定性系數(shù)DC共4個(gè)指標(biāo)對(duì)模型模擬精度進(jìn)行評(píng)價(jià)?。

    REp=qs-q0q0×100%?(1)

    REv=Qs-Q0Q0×100%??(2)

    △T=Ts-T0?(3)

    DC=1-ni=1[qs(i)-q0(i)]2

    ni=1[q0(i)-?q0?]2??(4)

    式中,?qs,q0,?q0?分別為計(jì)算洪峰流量、實(shí)測(cè)洪峰流量、實(shí)測(cè)平均洪峰流量;Qs,Q0分別為計(jì)算洪量和實(shí)測(cè)洪量;Ts,T0分別為計(jì)算

    峰現(xiàn)時(shí)間和實(shí)測(cè)峰現(xiàn)時(shí)間;i?代表時(shí)間步長(zhǎng)序列。

    3.1?不同時(shí)間步長(zhǎng)2種方案的模擬分析

    3.1.1?方案1的模擬結(jié)果分析

    次洪模擬計(jì)算步長(zhǎng)分別為10,15 min和20 min,方案1模擬誤差指標(biāo)絕對(duì)值的平均值見(jiàn)表2,圖2為12場(chǎng)次洪模擬的?DC?值。由表2、圖2可知:

    (1) 從整體模擬結(jié)果來(lái)看,時(shí)間步長(zhǎng)為15 min下的模擬結(jié)果4個(gè)指標(biāo)精度均高于10 min和20 min的,10 min步長(zhǎng)模擬的精度整體高于20 min的。

    (2) 從?DC?值來(lái)看,12場(chǎng)次降雨洪水中有9場(chǎng)15 min步長(zhǎng)模擬的?DC?值高于10 min和20 min的,而且根據(jù)《水文情報(bào)預(yù)報(bào)規(guī)范》(GB/T22482-2008)的規(guī)定,15 min步長(zhǎng)模擬的DC值達(dá)到甲級(jí)(?DC?>0.90)預(yù)報(bào)精度的有3場(chǎng),占25%,達(dá)到乙級(jí)(0.90>?DC?≧?0.70?)的有7場(chǎng),占58.33%;達(dá)到丙級(jí)(0.70>?DC?≧0.50)的有2場(chǎng),占16.7%。10 min步長(zhǎng)模擬的?DC?值達(dá)到甲級(jí)預(yù)報(bào)精度的有1場(chǎng),占8.3%,達(dá)到乙級(jí)的有5場(chǎng),占41.7%;達(dá)到丙級(jí)的有2場(chǎng),占16.7%。20 min步長(zhǎng)模擬的?DC?值達(dá)到甲級(jí)預(yù)報(bào)精度的有0場(chǎng),達(dá)到乙級(jí)的有5場(chǎng),占41.7%;達(dá)到丙級(jí)的有4場(chǎng),占?33.33%?。其中20 min步長(zhǎng)模擬的DC值浮動(dòng)較大,最大值達(dá)到了0.8以上,最小值小于0.15。

    (3) 從峰現(xiàn)時(shí)差來(lái)看,15 min步長(zhǎng)模擬的峰現(xiàn)時(shí)差絕對(duì)值的平均值在0.5 h以?xún)?nèi),10 min和20 min模擬的峰現(xiàn)時(shí)差絕對(duì)值的平均值在2 h以?xún)?nèi)。

    (4) 20 min步長(zhǎng)模擬的洪峰相對(duì)誤差絕對(duì)值的均值最大,比10 min的大4.03%,比15 min的大?7.38%?,而三者洪量相對(duì)誤差絕對(duì)值均值差別不大。

    因此,綜合考慮以上4個(gè)指標(biāo),方案1計(jì)算時(shí)段為15 min的模擬效果優(yōu)于10 min和20 min的。

    3.1.2?方案2的模擬結(jié)果分析

    不同時(shí)間步長(zhǎng)下方案2模擬誤差指標(biāo)絕對(duì)值的平均值見(jiàn)表3,圖3為方案2的12場(chǎng)次洪模擬的?DC?值。由表3、圖3可知:

    (1) 從整體模擬結(jié)果來(lái)看,時(shí)間步長(zhǎng)為15 min下的模擬結(jié)果4個(gè)指標(biāo)精度均高于10 min和20 min的。

    (2) 從?DC?值來(lái)看,12場(chǎng)次降雨洪水中有10場(chǎng)15 min步長(zhǎng)模擬的?DC?值高于10 min和20 min的,15 min步長(zhǎng)模擬的?DC?值達(dá)到甲級(jí)預(yù)報(bào)精度的有1場(chǎng),占8.3%,達(dá)到乙級(jí)的有8場(chǎng),占66.7%,達(dá)到丙級(jí)的有3場(chǎng),占25%。10 min步長(zhǎng)模擬的?DC?值達(dá)到甲級(jí)預(yù)報(bào)精度的有0場(chǎng),達(dá)到乙級(jí)的有4場(chǎng),占33.33%;達(dá)到丙級(jí)的有5場(chǎng),占41.7%。20 min步長(zhǎng)模擬的?DC?值達(dá)到甲級(jí)預(yù)報(bào)精度的有0場(chǎng),達(dá)到乙級(jí)的有4場(chǎng),占?33.33%?;達(dá)到丙級(jí)的有3場(chǎng),占33.33%。其中10 min步長(zhǎng)模擬的?DC?值浮動(dòng)較大,最高達(dá)到0.75以上,最低值小于0.2。

    (3) 從峰現(xiàn)時(shí)差來(lái)看,15 min步長(zhǎng)模擬的峰現(xiàn)時(shí)差絕對(duì)值的平均值在0.75 h以?xún)?nèi),10 min的在2 h以?xún)?nèi),20 min的在3 h以?xún)?nèi)。

    (4) 20 min步長(zhǎng)模擬的洪峰流量相對(duì)誤差絕對(duì)值均值比10 min的大2.96%,比15 min的大4.93%,洪量相對(duì)誤差絕對(duì)值均值也最大。因此,無(wú)論從洪峰相對(duì)誤差、洪量相對(duì)誤差、確定性系數(shù)以及峰現(xiàn)時(shí)差中的單個(gè)因素分析還是綜合分析,方案2計(jì)算時(shí)段為15 min的模擬效果是最優(yōu)的,其次是10 min。

    根據(jù)不同時(shí)間步長(zhǎng)下方案1和方案2的研究結(jié)果可知,2種方案都是在計(jì)算時(shí)段是15 min的模擬效果最好。其原因主要是研究區(qū)屬于半干旱半濕潤(rùn)地區(qū),蓄滿(mǎn)產(chǎn)流與超滲產(chǎn)流兩種產(chǎn)流機(jī)制并存。當(dāng)降雨強(qiáng)度大、歷時(shí)較短時(shí)以超滲產(chǎn)流為主,時(shí)間步長(zhǎng)選取較長(zhǎng)就會(huì)導(dǎo)致降雨強(qiáng)度被平均化, 超滲機(jī)制難以體現(xiàn)。而且,所形成的洪水過(guò)程線(xiàn)趨于坦化,在一定程度上影響了模型模擬效果。時(shí)間步長(zhǎng)選取較短時(shí),降雨強(qiáng)度的均化作用就將減弱,超滲作用相對(duì)增強(qiáng),增大了地面徑流量,導(dǎo)致所形成的洪水過(guò)程線(xiàn)過(guò)于集中,降低了模擬精度。

    3.2?最優(yōu)時(shí)間步長(zhǎng)下2種方案的結(jié)果對(duì)比

    計(jì)算時(shí)段是15 min時(shí)2種方案次洪模擬結(jié)果是最優(yōu)的,模擬結(jié)果見(jiàn)表4,圖4為2場(chǎng)典型的洪水過(guò)程線(xiàn)。由表4、圖4可知:2種方案下模擬的12場(chǎng)洪水洪峰流量、洪量相對(duì)誤差都在20%范圍以?xún)?nèi),峰現(xiàn)時(shí)差都在3 h以?xún)?nèi),確定性系數(shù)都達(dá)到0.6以上。綜合考慮洪峰、洪量相對(duì)誤差、峰現(xiàn)時(shí)差以及確定性系數(shù)4個(gè)指標(biāo),2種方案15 min步長(zhǎng)模擬的結(jié)果都合格,合格率100%。在綜合考慮以上4個(gè)指標(biāo)后,以合格率大小作為標(biāo)準(zhǔn),2種方案下的模擬均達(dá)到甲級(jí)(合格率≧85%)預(yù)報(bào)標(biāo)準(zhǔn)。具體表現(xiàn)為:

    (1) 從洪峰流量來(lái)看,方案1相對(duì)誤差絕對(duì)值均值比方案2低3.05%,有9場(chǎng)洪水的相對(duì)誤差絕對(duì)值低于方案2,占75%。方案1最大相對(duì)誤差值為?-16.15%?,最小值為0.50%,方案2相對(duì)誤差最大達(dá)到了-19.70%,最小為0.48%,表明方案1洪峰流量的模擬更為精確。

    (2) 從洪量來(lái)看,方案1相對(duì)誤差絕對(duì)值均值比方案2低1.7%,有7場(chǎng)洪水的相對(duì)誤差絕對(duì)值值低于方案2,占58.33%。方案1最大相對(duì)誤差值最大為19.30%,最小為0.74%,方案2相對(duì)誤差值最大為19.16%,最小值為0%,整體來(lái)看方案1洪量的模擬更為精確。

    (3) 從峰現(xiàn)時(shí)差來(lái)看,方案2只有3場(chǎng)洪水的相對(duì)誤差絕對(duì)值低于方案1,占25%,整體相對(duì)誤差絕對(duì)值的均值比方案1高16 min,而且方案1中洪號(hào)為19980714,19980821,20031011的洪水峰現(xiàn)時(shí)間與實(shí)測(cè)洪峰同步,表明方案1峰現(xiàn)時(shí)差的模擬比方案2精確。

    (4) 從?DC值看,方案1中DC值均值比方案2高0.049,其中方案1有7場(chǎng)洪水的DC?值高于方案2,占?58.33%?,方案1達(dá)到乙級(jí)預(yù)報(bào)精度以上的洪水比方案2多1場(chǎng),占8.33%。

    因此,綜合考慮以上4個(gè)指標(biāo),計(jì)算時(shí)段為15 min時(shí)方案1的模擬精度高于方案2。

    半干旱半濕潤(rùn)地區(qū),降雨多為強(qiáng)度大、歷時(shí)短的類(lèi)型,加上時(shí)空分布不均勻,易形成陡漲陡落、峰高量小的洪水過(guò)程,HEC-HMS模型對(duì)該類(lèi)型的洪水模擬結(jié)果較好,如19900910場(chǎng)洪水(圖 4(a)),兩種方案模擬的精度均較高。對(duì)于20030826場(chǎng)洪水(圖4(b)),退水期降雨量雖然較少,但是實(shí)測(cè)退水段流量仍然較大、退水時(shí)間較長(zhǎng),退水過(guò)程可能主要是以蓄滿(mǎn)產(chǎn)流為主,對(duì)于退水段流量較大且速度較慢的這類(lèi)型洪水,該模型2種方案下的模擬結(jié)果也較好,說(shuō)明HEC-HMS模型在北張店流域適用性較好。根據(jù)15 min時(shí)間步長(zhǎng)下2種方案的次洪模擬結(jié)果可知,方案1的模擬精度高于方案2。這是由于:

    (1) 北張店流域?qū)儆诎敫珊蛋霛駶?rùn)地區(qū),再加上受季風(fēng)氣候的影響,降雨多強(qiáng)度大、歷時(shí)短,多以超滲產(chǎn)流為主。方案1產(chǎn)流計(jì)算為初損穩(wěn)滲法,是計(jì)算超滲產(chǎn)流的典型方法,所需率定的參數(shù)少,計(jì)算方法簡(jiǎn)單合理,模擬結(jié)果較為精確。方案2所選用的格林-安普特法計(jì)算方法較復(fù)雜,而且由于缺少相關(guān)資料導(dǎo)致參數(shù)的選取中會(huì)出現(xiàn)一定程度的偏差,所以模擬精度相對(duì)較低。

    (2) 在直接徑流計(jì)算方法中,方案1選取了Clark單位線(xiàn)法,方案2選取了SCS單位線(xiàn)法。其中SCS單位線(xiàn)法計(jì)算簡(jiǎn)單,是通過(guò)劃分的子流域,計(jì)算各子流域平均洪峰滯時(shí)。而Clark單位線(xiàn)考慮了天然河道對(duì)洪水的平移與調(diào)蓄作用,所以模擬的峰現(xiàn)時(shí)差精度更高。

    4?結(jié) 論

    (1) 半干旱半濕潤(rùn)地區(qū)降雨強(qiáng)度大、歷時(shí)短,且時(shí)空分布極不均勻的這些特性在很大程度上決定了洪水特性,時(shí)間步長(zhǎng)的長(zhǎng)短對(duì)洪峰流量、洪量和峰現(xiàn)時(shí)差都有很大的影響。

    (2) 分別以10,15,20 min為計(jì)算時(shí)段步長(zhǎng), 應(yīng)用HEC-HMS模型分別在2種方案下對(duì)北張店流域的12場(chǎng)洪水進(jìn)行次洪模擬,分析研究了不同時(shí)間步長(zhǎng)對(duì)次洪過(guò)程模擬的影響。結(jié)果表明,時(shí)間步長(zhǎng)為15 min時(shí)模擬精度最高,并非是時(shí)間步長(zhǎng)越短,次洪模擬效果越好,洪水模擬中存在最優(yōu)時(shí)間步長(zhǎng)。因此,時(shí)間步長(zhǎng)選取較長(zhǎng)或較短都會(huì)影響模型模擬精度。

    (3) 半干旱半濕潤(rùn)地區(qū)產(chǎn)流方式較為復(fù)雜,蓄滿(mǎn)產(chǎn)流與超滲產(chǎn)流并存,同一場(chǎng)洪水可能前期為超滲產(chǎn)流 ,后期變?yōu)樾顫M(mǎn)產(chǎn)流。

    (4) 產(chǎn)匯流計(jì)算方法對(duì)洪峰、洪量、峰現(xiàn)時(shí)差以及確定性系數(shù)DC值均有較大的影響。表明在之后運(yùn)用HEC-HMS水文模型對(duì)研究區(qū)進(jìn)行次洪模擬時(shí),應(yīng)根據(jù)研究區(qū)下墊面特征所確定的相關(guān)水文參數(shù)來(lái)選擇合理的水文計(jì)算方法進(jìn)行優(yōu)選搭配,從而使構(gòu)建的模型更準(zhǔn)確地反映出研究區(qū)復(fù)雜下墊面特征下的降雨徑流機(jī)制,提高次洪模擬的精度。

    參考文獻(xiàn):

    [1]李巧玲,李致家,陳利者,等.半濕潤(rùn)半干旱流域降雨徑流關(guān)系及下墊面相似性[J].河海大學(xué)學(xué)報(bào):自然科學(xué)版,2015,43(2):95-99.

    [2]Zema DA, Labate A, Martino D, et al. Comparing DifferentInfiltration Methods Of The Hec-Hms Model: The Case Study of The Mésima Torrent (Southern Italy)[J].Land Degradation & Development, 2017,28,(1):294-308.

    [3]Hessel R. Effects of grid cell size and time step length on simulation results of the Limburg soil erosion model (LISEM)[J].Hydrological Processes,2010, 19(15):3037-3049.

    [4]Ren L L, Zhang W, Li C H, et al.Comparison of RunoffParameterization Schemes with Spatial Heterogeneity across Different Temporal Scales in Semi humid and Semiarid Regions[J].Journal of Hydrologic Engineering,2008, 13(5):400-409.

    [5]林峰,陳瑩,陳興偉,等.基于HEC-HMS模型的時(shí)間步長(zhǎng)對(duì)次洪模擬的影響分析[J].山地學(xué)報(bào),2011,29(1):50-54.

    [6]Jin H, Liang R, Wang Y, et al. Flood-Runoff in Semi-Arid and Sub-Humid Regions, a Case Study: A Simulation of Jianghe Watershed in Northern China[J]. Water,2015,7(9):5155-5172.

    [7]鄧霞,董曉華,薄會(huì)娟.目標(biāo)函數(shù)對(duì)HEC-HMS模型參數(shù)率定的影響研究[J].水電能源科學(xué),2010,28(8):23-25.

    [8]吳博,鄭秀清,劉姍姍,等.HEC-HMS水文模型在流域洪水模擬中的應(yīng)用——以山西石樓縣義牒河流域?yàn)槔齕J].人民長(zhǎng)江,2017,48(23):1-4.

    [9]李軍.半干旱半濕潤(rùn)地區(qū)產(chǎn)匯流過(guò)程機(jī)理研究[D].北京:中國(guó)科學(xué)院大學(xué),2015.

    [10]Rostami N,Rostamizad G.Comparison methods Green & Ampt and Initial and constant rate loss for estimate precipitation loss to simulation runoff[J]. International Journal of Agriculture Research & Review, 2011,1(1):33-37.

    [11]Reza Kabiri, Andrew Chan, Ramani Bai. Comparison of SCS andGreen-Ampt Methods in Surface Runoff-Flooding Simulation for Klang Watershed in Malaysia[J]. Open Journal of Modern Hydrology, 2013,3(3):102-114.

    [12]王中根,劉昌明,左其亭,等.基于DEM的分布式水文模型構(gòu)建方法[J].地理科學(xué)進(jìn)展,2002,21(5):430-439.

    [13]穆艾塔爾·賽地,阿不都·沙拉木,丁建麗,等.HEC-HMS水文模型在數(shù)據(jù)稀缺山區(qū)流域中的應(yīng)用——以烏魯木齊河流域?yàn)槔齕J].水土保持通報(bào),2015,35(6):140-143.

    [14]趙人俊.流域水文模擬——新安江模型與陜北模型[M].北京:水利電力出版社,1984.

    [15]Haan C T, Johnson H P, Brakensiek D L. Hydrologic modeling of small watersheds[M].American Society of Agricultural Engineers, 1982:342-344

    [16]GB/T22482-2008 水文情報(bào)預(yù)報(bào)規(guī)范[S].北京:中國(guó)標(biāo)準(zhǔn)出版社,2008.

    引用本文:趙?娜,趙雪花,祝雪萍.時(shí)間步長(zhǎng)和產(chǎn)匯流方法對(duì)小流域次洪模擬的影響[J].人民長(zhǎng)江,2019,50(1):64-69.

    Influence of time-steps and runoff generation and confluencecalculation methods on flood simulation in small watersheds

    ZHAO Na,ZHAO Xuehua,ZHU Xueping

    (College of Water Resources Science and Engineering, Taiyuan University of Technology, Taiyuan 030024,China)

    Abstract:In order to study the influence of different time steps and runoff generation and confluence calculation methods on flood simulation in semi-humid and semi-arid regions, we simulated 12 representational flood events in Beizhangdian Basin from 1993 to 2005. The HEC-HMS model was applied to simulate the hydrological process under different time steps (10,15,20 min) and two schemes. The first scheme was the initial-loss and stable seepage method&Clark unit line method, the second one was the Green and Ampt method&SCS unit line method, and both schemes used the Muskingum method for river channel routing and the exponential regression method for base flow calculation. The results show that the time step has a great influence on the simulation results and both schemes had the highest simulation accuracy under 15min time step, indicating that the shorting time step does not always lead to accuracy improvement. When the time step was 15min, the average deterministic coefficient of the first scheme was 0.813, and the second scheme was 0.764. The first scheme is superior to the second one in terms of flood peak, flood volume and flood peak time.

    Key words:?time steps; runoff generation and concentration; flood simulation; HEC-HMS model;semi-humid and semi-arid region;Beizhangdian Watershed;Shanxi Province

    猜你喜歡
    洪量產(chǎn)流步長(zhǎng)
    基于SPA 的北江流域峰量關(guān)系研究
    陜西水利(2023年12期)2023-12-19 03:28:32
    產(chǎn)流及其研究進(jìn)展
    基于Armijo搜索步長(zhǎng)的BFGS與DFP擬牛頓法的比較研究
    不同坡面單元人工降雨產(chǎn)流試驗(yàn)與分析
    北京山區(qū)側(cè)柏林地坡面初始產(chǎn)流時(shí)間影響因素
    遼河干流主要控制站近75年最大洪峰及洪量變化特征分析研究
    地表粗糙度對(duì)黃土坡面產(chǎn)流機(jī)制的影響
    基于逐維改進(jìn)的自適應(yīng)步長(zhǎng)布谷鳥(niǎo)搜索算法
    一種新型光伏系統(tǒng)MPPT變步長(zhǎng)滯環(huán)比較P&O法
    適用于電算的設(shè)計(jì)洪水過(guò)程線(xiàn)放縮方法
    国产片内射在线| 色综合欧美亚洲国产小说| 精品国产一区二区三区久久久樱花| 欧美老熟妇乱子伦牲交| 在线观看人妻少妇| 国产精品1区2区在线观看. | 国产伦理片在线播放av一区| 丝袜喷水一区| 欧美+亚洲+日韩+国产| 欧美成人午夜精品| 中文字幕精品免费在线观看视频| 黑人欧美特级aaaaaa片| 久久久精品94久久精品| 最新在线观看一区二区三区| e午夜精品久久久久久久| 老司机午夜十八禁免费视频| 2018国产大陆天天弄谢| 1024视频免费在线观看| 捣出白浆h1v1| 在线天堂中文资源库| 国产精品影院久久| 国产精品久久久久久精品电影小说| 亚洲国产欧美在线一区| 热99re8久久精品国产| 一区二区av电影网| 99精品久久久久人妻精品| 国产免费福利视频在线观看| 久久国产精品大桥未久av| 午夜福利在线免费观看网站| 精品国产一区二区三区久久久樱花| 日韩人妻精品一区2区三区| 国产国语露脸激情在线看| 国产极品粉嫩免费观看在线| aaaaa片日本免费| 亚洲免费av在线视频| 亚洲综合色网址| 91麻豆av在线| 午夜久久久在线观看| 欧美国产精品一级二级三级| 国产精品一区二区在线不卡| 50天的宝宝边吃奶边哭怎么回事| 国产成人精品久久二区二区免费| 一区二区三区国产精品乱码| 色视频在线一区二区三区| 91大片在线观看| av一本久久久久| 亚洲欧美日韩高清在线视频 | 一个人免费看片子| 国产亚洲精品第一综合不卡| 国产真人三级小视频在线观看| 狠狠精品人妻久久久久久综合| 精品亚洲成a人片在线观看| 国产野战对白在线观看| 国产真人三级小视频在线观看| 久久毛片免费看一区二区三区| 高清视频免费观看一区二区| 日韩视频一区二区在线观看| 99国产综合亚洲精品| 极品教师在线免费播放| 国产免费福利视频在线观看| 国产成人精品久久二区二区91| 亚洲av国产av综合av卡| 日本av免费视频播放| 黄片大片在线免费观看| 欧美成狂野欧美在线观看| 女警被强在线播放| 色婷婷av一区二区三区视频| 国产单亲对白刺激| 免费久久久久久久精品成人欧美视频| 日本av手机在线免费观看| 精品少妇内射三级| 国产精品国产高清国产av | 在线十欧美十亚洲十日本专区| 高清欧美精品videossex| 97人妻天天添夜夜摸| 午夜两性在线视频| 成人av一区二区三区在线看| 91字幕亚洲| 亚洲av成人不卡在线观看播放网| 日本五十路高清| 国产成人系列免费观看| 日韩成人在线观看一区二区三区| 男女午夜视频在线观看| 日韩视频在线欧美| 欧美黑人精品巨大| 亚洲av美国av| 精品国产乱码久久久久久男人| 极品教师在线免费播放| 夜夜爽天天搞| 大香蕉久久成人网| tocl精华| 亚洲五月婷婷丁香| 别揉我奶头~嗯~啊~动态视频| 老司机午夜十八禁免费视频| 夜夜夜夜夜久久久久| 欧美另类亚洲清纯唯美| av一本久久久久| 高潮久久久久久久久久久不卡| 真人做人爱边吃奶动态| 在线观看舔阴道视频| 免费日韩欧美在线观看| 欧美日韩精品网址| 老汉色∧v一级毛片| 精品国产超薄肉色丝袜足j| 十八禁高潮呻吟视频| 久久精品熟女亚洲av麻豆精品| 久久精品人人爽人人爽视色| 亚洲国产精品一区二区三区在线| 精品国产一区二区三区久久久樱花| 免费观看a级毛片全部| 搡老乐熟女国产| 一二三四在线观看免费中文在| 无遮挡黄片免费观看| 女人久久www免费人成看片| 国产精品二区激情视频| 国产男女超爽视频在线观看| 国产成人精品无人区| 视频在线观看一区二区三区| 大香蕉久久成人网| 首页视频小说图片口味搜索| av天堂久久9| 香蕉久久夜色| 欧美亚洲日本最大视频资源| 在线亚洲精品国产二区图片欧美| 纯流量卡能插随身wifi吗| 亚洲全国av大片| 天天影视国产精品| 一进一出好大好爽视频| 亚洲精品粉嫩美女一区| 日本a在线网址| 成年人黄色毛片网站| av网站在线播放免费| 黄色a级毛片大全视频| 啦啦啦中文免费视频观看日本| 亚洲专区字幕在线| 亚洲精品在线美女| 91成人精品电影| 亚洲国产看品久久| 成人亚洲精品一区在线观看| 精品国产乱码久久久久久男人| 这个男人来自地球电影免费观看| 午夜老司机福利片| 黄色a级毛片大全视频| 美国免费a级毛片| 一级黄色大片毛片| 人人妻,人人澡人人爽秒播| 久久精品熟女亚洲av麻豆精品| 免费在线观看完整版高清| 黑丝袜美女国产一区| 国产黄色免费在线视频| 午夜老司机福利片| 男女高潮啪啪啪动态图| 午夜福利视频精品| 777久久人妻少妇嫩草av网站| 91字幕亚洲| 欧美黑人精品巨大| 大陆偷拍与自拍| 黄片大片在线免费观看| 国产精品秋霞免费鲁丝片| 亚洲精品国产色婷婷电影| 99国产精品99久久久久| 麻豆乱淫一区二区| 亚洲第一av免费看| 国产欧美日韩精品亚洲av| 欧美精品一区二区大全| 久9热在线精品视频| 国产成人啪精品午夜网站| 女人爽到高潮嗷嗷叫在线视频| 热re99久久精品国产66热6| 在线观看人妻少妇| 亚洲精品美女久久av网站| 国产91精品成人一区二区三区 | 国产一区有黄有色的免费视频| 免费久久久久久久精品成人欧美视频| 中文字幕制服av| 久久天躁狠狠躁夜夜2o2o| 国产在线免费精品| 美女高潮喷水抽搐中文字幕| 国产亚洲欧美精品永久| 夫妻午夜视频| 怎么达到女性高潮| a级毛片黄视频| 国产男靠女视频免费网站| 99精品欧美一区二区三区四区| 国产亚洲精品久久久久5区| 精品熟女少妇八av免费久了| 久久久久精品国产欧美久久久| 无限看片的www在线观看| 狠狠精品人妻久久久久久综合| 国产精品久久久人人做人人爽| 精品国产一区二区三区久久久樱花| 欧美日韩亚洲综合一区二区三区_| 国产精品偷伦视频观看了| 欧美一级毛片孕妇| 啦啦啦在线免费观看视频4| 亚洲熟女毛片儿| 日韩视频在线欧美| 久久久国产精品麻豆| 1024视频免费在线观看| 99热网站在线观看| 窝窝影院91人妻| 午夜两性在线视频| 国产精品香港三级国产av潘金莲| 日韩熟女老妇一区二区性免费视频| 电影成人av| bbb黄色大片| av天堂久久9| 亚洲五月色婷婷综合| 国产xxxxx性猛交| 久久国产精品大桥未久av| 久久午夜亚洲精品久久| 一本—道久久a久久精品蜜桃钙片| 脱女人内裤的视频| 亚洲欧美激情在线| 80岁老熟妇乱子伦牲交| 免费观看人在逋| 18禁黄网站禁片午夜丰满| 成人影院久久| 女性被躁到高潮视频| 十八禁高潮呻吟视频| 在线天堂中文资源库| 香蕉国产在线看| 伦理电影免费视频| 我要看黄色一级片免费的| 男人舔女人的私密视频| 国产高清激情床上av| 99在线人妻在线中文字幕 | 亚洲成av片中文字幕在线观看| av超薄肉色丝袜交足视频| 国产亚洲午夜精品一区二区久久| 中文字幕精品免费在线观看视频| 老汉色∧v一级毛片| 亚洲视频免费观看视频| 在线看a的网站| 国产单亲对白刺激| 午夜福利影视在线免费观看| 国产精品影院久久| 国产精品九九99| 在线观看免费日韩欧美大片| 亚洲精品国产色婷婷电影| 在线观看66精品国产| 精品视频人人做人人爽| 中文字幕另类日韩欧美亚洲嫩草| 亚洲伊人久久精品综合| 一区二区三区国产精品乱码| 国产无遮挡羞羞视频在线观看| 欧美激情极品国产一区二区三区| 美女主播在线视频| 亚洲欧美日韩高清在线视频 | 男女之事视频高清在线观看| 夜夜夜夜夜久久久久| 久久 成人 亚洲| 每晚都被弄得嗷嗷叫到高潮| 99re6热这里在线精品视频| 桃红色精品国产亚洲av| 夫妻午夜视频| 搡老熟女国产l中国老女人| 日韩制服丝袜自拍偷拍| 亚洲av电影在线进入| aaaaa片日本免费| av天堂久久9| 国产日韩欧美视频二区| 蜜桃在线观看..| 青草久久国产| 欧美久久黑人一区二区| 国产亚洲av高清不卡| 另类亚洲欧美激情| 亚洲av国产av综合av卡| 男女高潮啪啪啪动态图| 精品少妇一区二区三区视频日本电影| 久久精品熟女亚洲av麻豆精品| 纵有疾风起免费观看全集完整版| 国产成人欧美在线观看 | 一本综合久久免费| 国产男靠女视频免费网站| 天天添夜夜摸| 人人妻人人澡人人看| 一区二区三区国产精品乱码| 嫩草影视91久久| 国产精品自产拍在线观看55亚洲 | 菩萨蛮人人尽说江南好唐韦庄| 美女福利国产在线| 亚洲 国产 在线| 欧美日韩中文字幕国产精品一区二区三区 | 色婷婷av一区二区三区视频| 热re99久久国产66热| 欧美亚洲日本最大视频资源| 男人舔女人的私密视频| 亚洲av电影在线进入| 男女下面插进去视频免费观看| 日本撒尿小便嘘嘘汇集6| netflix在线观看网站| 亚洲少妇的诱惑av| 久久人妻福利社区极品人妻图片| 久久久久久亚洲精品国产蜜桃av| 日本欧美视频一区| 国产高清激情床上av| 一二三四在线观看免费中文在| 女人久久www免费人成看片| 自拍欧美九色日韩亚洲蝌蚪91| 又紧又爽又黄一区二区| 色94色欧美一区二区| 国产激情久久老熟女| 动漫黄色视频在线观看| 女性生殖器流出的白浆| 久久久久视频综合| 中文亚洲av片在线观看爽 | videos熟女内射| 最近最新免费中文字幕在线| av国产精品久久久久影院| 狠狠精品人妻久久久久久综合| 国产精品一区二区在线观看99| 国产亚洲精品第一综合不卡| 欧美激情极品国产一区二区三区| 99久久精品国产亚洲精品| 欧美 亚洲 国产 日韩一| 欧美黑人精品巨大| cao死你这个sao货| 2018国产大陆天天弄谢| 激情在线观看视频在线高清 | 91字幕亚洲| 狠狠狠狠99中文字幕| 色尼玛亚洲综合影院| 热re99久久精品国产66热6| 99在线人妻在线中文字幕 | e午夜精品久久久久久久| 精品国产国语对白av| 韩国精品一区二区三区| 国内毛片毛片毛片毛片毛片| 美女扒开内裤让男人捅视频| 亚洲中文字幕日韩| 欧美激情极品国产一区二区三区| 制服诱惑二区| 大片免费播放器 马上看| 波多野结衣av一区二区av| 男女下面插进去视频免费观看| 国产精品一区二区在线不卡| 少妇粗大呻吟视频| 亚洲九九香蕉| 午夜久久久在线观看| 中文亚洲av片在线观看爽 | av视频免费观看在线观看| 久久亚洲真实| 一本—道久久a久久精品蜜桃钙片| 黑人巨大精品欧美一区二区mp4| 国产又色又爽无遮挡免费看| 日韩人妻精品一区2区三区| 成人影院久久| 老司机靠b影院| 久久中文字幕一级| 伦理电影免费视频| 日韩免费av在线播放| 最近最新中文字幕大全免费视频| 高清在线国产一区| 丝瓜视频免费看黄片| 欧美激情高清一区二区三区| 久久久精品94久久精品| 久久久久久久大尺度免费视频| 日本黄色日本黄色录像| 亚洲欧美一区二区三区黑人| 18禁裸乳无遮挡动漫免费视频| 亚洲va日本ⅴa欧美va伊人久久| h视频一区二区三区| 久久中文字幕一级| 成人国语在线视频| 91九色精品人成在线观看| 嫩草影视91久久| 电影成人av| 亚洲精品av麻豆狂野| 国产一区有黄有色的免费视频| 最近最新中文字幕大全免费视频| 国产高清国产精品国产三级| 十八禁高潮呻吟视频| 日韩大码丰满熟妇| 一区二区三区国产精品乱码| 自线自在国产av| 成在线人永久免费视频| 欧美黄色淫秽网站| 激情在线观看视频在线高清 | 超碰97精品在线观看| 久久中文字幕一级| 婷婷丁香在线五月| 九色亚洲精品在线播放| 女人久久www免费人成看片| 大片电影免费在线观看免费| 精品少妇一区二区三区视频日本电影| 宅男免费午夜| 国产麻豆69| 黑人巨大精品欧美一区二区mp4| 大香蕉久久网| 亚洲综合色网址| 超碰97精品在线观看| 一夜夜www| av天堂久久9| 日本五十路高清| 国产黄色免费在线视频| 日日摸夜夜添夜夜添小说| 精品一区二区三卡| 精品熟女少妇八av免费久了| 中文字幕精品免费在线观看视频| 亚洲成人国产一区在线观看| 亚洲免费av在线视频| 国产精品一区二区精品视频观看| 19禁男女啪啪无遮挡网站| 亚洲国产中文字幕在线视频| 久久婷婷成人综合色麻豆| 香蕉丝袜av| 日本av手机在线免费观看| 99re在线观看精品视频| 久久久国产成人免费| 欧美日韩av久久| 成人免费观看视频高清| 法律面前人人平等表现在哪些方面| 久久久久久久国产电影| 最近最新免费中文字幕在线| 亚洲国产欧美网| 宅男免费午夜| 亚洲免费av在线视频| 日韩一区二区三区影片| 久久久国产精品麻豆| 国产精品久久电影中文字幕 | 一区二区三区乱码不卡18| 国产欧美日韩一区二区三| 精品一区二区三区四区五区乱码| 亚洲精品成人av观看孕妇| 欧美日韩中文字幕国产精品一区二区三区 | 国产精品免费一区二区三区在线 | 亚洲精品久久午夜乱码| 自线自在国产av| 桃红色精品国产亚洲av| 啦啦啦视频在线资源免费观看| 日韩制服丝袜自拍偷拍| 99re6热这里在线精品视频| 欧美黄色淫秽网站| 亚洲国产欧美日韩在线播放| 老鸭窝网址在线观看| 国产精品98久久久久久宅男小说| 久久久久久久久免费视频了| 色在线成人网| 久久久国产成人免费| 麻豆成人av在线观看| 黑丝袜美女国产一区| 午夜成年电影在线免费观看| 99国产极品粉嫩在线观看| 免费人妻精品一区二区三区视频| 91九色精品人成在线观看| 99re6热这里在线精品视频| av又黄又爽大尺度在线免费看| 日韩免费高清中文字幕av| 午夜福利一区二区在线看| 欧美亚洲 丝袜 人妻 在线| 黄色毛片三级朝国网站| av欧美777| 菩萨蛮人人尽说江南好唐韦庄| 妹子高潮喷水视频| 人成视频在线观看免费观看| 香蕉久久夜色| 欧美黑人欧美精品刺激| 自拍欧美九色日韩亚洲蝌蚪91| 女人高潮潮喷娇喘18禁视频| 一边摸一边抽搐一进一小说 | 美国免费a级毛片| 国产精品久久久人人做人人爽| 久久这里只有精品19| 妹子高潮喷水视频| 亚洲精华国产精华精| 黄色怎么调成土黄色| av在线播放免费不卡| 国产精品秋霞免费鲁丝片| 国产一区二区三区视频了| 日韩人妻精品一区2区三区| 午夜91福利影院| 十八禁网站免费在线| 亚洲精品国产色婷婷电影| 日本欧美视频一区| av福利片在线| 午夜两性在线视频| 夜夜爽天天搞| 日韩免费高清中文字幕av| 啦啦啦 在线观看视频| 美女扒开内裤让男人捅视频| 无遮挡黄片免费观看| 午夜老司机福利片| 12—13女人毛片做爰片一| 国产精品 国内视频| 亚洲自偷自拍图片 自拍| 久久精品人人爽人人爽视色| 精品国产乱子伦一区二区三区| 久久精品aⅴ一区二区三区四区| 国产区一区二久久| 男人舔女人的私密视频| 亚洲色图 男人天堂 中文字幕| 亚洲五月色婷婷综合| 啪啪无遮挡十八禁网站| 老司机深夜福利视频在线观看| 性高湖久久久久久久久免费观看| 久久影院123| 久久久欧美国产精品| 欧美另类亚洲清纯唯美| 男女午夜视频在线观看| 中文欧美无线码| 国产精品亚洲一级av第二区| 亚洲欧美一区二区三区黑人| 啦啦啦在线免费观看视频4| 国产不卡av网站在线观看| 亚洲avbb在线观看| 精品一区二区三卡| av福利片在线| 变态另类成人亚洲欧美熟女 | 天天操日日干夜夜撸| 天天躁日日躁夜夜躁夜夜| 女性生殖器流出的白浆| 日本av手机在线免费观看| 亚洲男人天堂网一区| 制服人妻中文乱码| 亚洲国产看品久久| 国产色视频综合| 亚洲五月色婷婷综合| 久久天躁狠狠躁夜夜2o2o| 国产精品1区2区在线观看. | 国产一卡二卡三卡精品| 亚洲国产欧美在线一区| 巨乳人妻的诱惑在线观看| 亚洲免费av在线视频| 天天影视国产精品| 亚洲国产欧美在线一区| 日本欧美视频一区| 国产精品免费大片| 久久狼人影院| 中文字幕色久视频| 久久久久久久国产电影| 后天国语完整版免费观看| 亚洲三区欧美一区| 成人手机av| 国产一区二区三区视频了| 国产欧美日韩一区二区精品| 美女视频免费永久观看网站| 最新的欧美精品一区二区| 一二三四在线观看免费中文在| 亚洲精品自拍成人| 亚洲国产欧美网| 国产精品九九99| 王馨瑶露胸无遮挡在线观看| 亚洲性夜色夜夜综合| 一本久久精品| 狠狠婷婷综合久久久久久88av| 中文字幕精品免费在线观看视频| 亚洲av日韩在线播放| 一边摸一边抽搐一进一出视频| 母亲3免费完整高清在线观看| 曰老女人黄片| 人成视频在线观看免费观看| 欧美日韩av久久| 亚洲一区二区三区欧美精品| 国产精品1区2区在线观看. | 久久av网站| 日日夜夜操网爽| 成人免费观看视频高清| 亚洲av片天天在线观看| 成人亚洲精品一区在线观看| 大型av网站在线播放| 咕卡用的链子| 久久精品国产亚洲av高清一级| 9热在线视频观看99| 视频在线观看一区二区三区| 国产精品亚洲一级av第二区| 中文字幕色久视频| 午夜两性在线视频| 国产精品成人在线| 12—13女人毛片做爰片一| 日本wwww免费看| 少妇精品久久久久久久| 老司机福利观看| 亚洲第一av免费看| 国产高清videossex| 日本精品一区二区三区蜜桃| 国产免费视频播放在线视频| 18禁美女被吸乳视频| 亚洲欧美精品综合一区二区三区| 桃红色精品国产亚洲av| 亚洲国产欧美日韩在线播放| 99久久国产精品久久久| 天天影视国产精品| 亚洲欧洲精品一区二区精品久久久| 一本一本久久a久久精品综合妖精| 亚洲色图综合在线观看| 国产真人三级小视频在线观看| 成人18禁高潮啪啪吃奶动态图| 日韩三级视频一区二区三区| 国产成人精品久久二区二区免费| 欧美日韩亚洲综合一区二区三区_| 欧美性长视频在线观看| 久久天躁狠狠躁夜夜2o2o| 国产成人欧美| 亚洲免费av在线视频| 久热这里只有精品99| 色视频在线一区二区三区| 亚洲精品粉嫩美女一区| 高清欧美精品videossex| 色精品久久人妻99蜜桃| 免费在线观看影片大全网站| 日本撒尿小便嘘嘘汇集6| 亚洲国产看品久久| 制服人妻中文乱码| 黑丝袜美女国产一区| 老熟妇乱子伦视频在线观看| 欧美精品一区二区免费开放| 91字幕亚洲| av网站在线播放免费| 1024香蕉在线观看| 丝袜美足系列| 大型av网站在线播放| 久久婷婷成人综合色麻豆| 色婷婷av一区二区三区视频| 午夜福利欧美成人|