• <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)放縮方法
    女人被躁到高潮嗷嗷叫费观| 男女做爰动态图高潮gif福利片 | 男女床上黄色一级片免费看| 国产黄a三级三级三级人| 久久久精品欧美日韩精品| 精品免费久久久久久久清纯| 99国产精品一区二区蜜桃av| 成年人黄色毛片网站| 女警被强在线播放| 女人被狂操c到高潮| 国产精品久久电影中文字幕| 精品国产乱子伦一区二区三区| 在线国产一区二区在线| 亚洲色图av天堂| 久久婷婷成人综合色麻豆| 午夜福利18| tocl精华| 国产国语露脸激情在线看| 91老司机精品| 日本免费一区二区三区高清不卡 | 在线观看舔阴道视频| 纯流量卡能插随身wifi吗| 搡老岳熟女国产| www.999成人在线观看| 男人的好看免费观看在线视频 | 久久精品国产99精品国产亚洲性色 | 日本黄色视频三级网站网址| 天堂动漫精品| 可以在线观看的亚洲视频| 高潮久久久久久久久久久不卡| 午夜两性在线视频| 国产欧美日韩精品亚洲av| 午夜免费激情av| 伦理电影免费视频| 亚洲电影在线观看av| 国产精品久久久人人做人人爽| 69av精品久久久久久| 狂野欧美激情性xxxx| 欧美日本视频| 欧美在线黄色| 18美女黄网站色大片免费观看| 国产成人啪精品午夜网站| 亚洲国产精品合色在线| 色综合亚洲欧美另类图片| 免费搜索国产男女视频| 日本在线视频免费播放| 一区二区三区激情视频| 天堂动漫精品| 午夜精品久久久久久毛片777| 久久青草综合色| 免费无遮挡裸体视频| 免费看a级黄色片| 午夜福利免费观看在线| 亚洲自拍偷在线| 国产成人系列免费观看| 夜夜爽天天搞| 免费久久久久久久精品成人欧美视频| av免费在线观看网站| 久久亚洲精品不卡| 给我免费播放毛片高清在线观看| 一级毛片女人18水好多| 国产一区二区三区视频了| 99热只有精品国产| 日韩有码中文字幕| 搡老熟女国产l中国老女人| 人人妻,人人澡人人爽秒播| 欧美日韩中文字幕国产精品一区二区三区 | 精品无人区乱码1区二区| 麻豆一二三区av精品| 日韩免费av在线播放| 国产一区二区在线av高清观看| 国产熟女午夜一区二区三区| 精品国产国语对白av| 亚洲少妇的诱惑av| 一个人观看的视频www高清免费观看 | 美女国产高潮福利片在线看| 国产精品,欧美在线| 大香蕉久久成人网| 国产99久久九九免费精品| 亚洲国产欧美日韩在线播放| 极品教师在线免费播放| 老司机福利观看| 97人妻天天添夜夜摸| 日韩欧美一区视频在线观看| 两性夫妻黄色片| 欧美一级a爱片免费观看看 | 国产成人欧美| 国产精品日韩av在线免费观看 | 男人舔女人下体高潮全视频| 99riav亚洲国产免费| 久久影院123| 成人特级黄色片久久久久久久| 国产精品,欧美在线| 香蕉国产在线看| 人妻丰满熟妇av一区二区三区| 999精品在线视频| 日韩三级视频一区二区三区| 国产97色在线日韩免费| 午夜激情av网站| 日韩欧美一区二区三区在线观看| 99在线视频只有这里精品首页| 欧美最黄视频在线播放免费| 天天躁狠狠躁夜夜躁狠狠躁| 国产精品久久视频播放| 看免费av毛片| 国产精品久久久久久亚洲av鲁大| 亚洲七黄色美女视频| www国产在线视频色| 熟妇人妻久久中文字幕3abv| 男人舔女人下体高潮全视频| 亚洲片人在线观看| 午夜精品国产一区二区电影| 亚洲第一电影网av| 波多野结衣一区麻豆| 1024视频免费在线观看| 亚洲精品久久国产高清桃花| 国产亚洲av嫩草精品影院| 中国美女看黄片| 久久精品91无色码中文字幕| 欧美一级a爱片免费观看看 | 国产一卡二卡三卡精品| 91国产中文字幕| 欧美大码av| 激情视频va一区二区三区| 91精品三级在线观看| 精品欧美国产一区二区三| 亚洲在线自拍视频| 久久久久亚洲av毛片大全| 黄色视频,在线免费观看| 亚洲免费av在线视频| av电影中文网址| 人人妻人人澡欧美一区二区 | 9热在线视频观看99| 国产主播在线观看一区二区| 波多野结衣一区麻豆| 久久 成人 亚洲| 午夜免费成人在线视频| 男女下面插进去视频免费观看| aaaaa片日本免费| 精品卡一卡二卡四卡免费| 国产男靠女视频免费网站| 久久中文字幕人妻熟女| 一卡2卡三卡四卡精品乱码亚洲| 欧美日韩精品网址| 亚洲熟妇熟女久久| 久久久久国产一级毛片高清牌| 9热在线视频观看99| 亚洲av美国av| 女人精品久久久久毛片| 亚洲熟女毛片儿| 亚洲精品美女久久久久99蜜臀| 久99久视频精品免费| netflix在线观看网站| 国产xxxxx性猛交| 亚洲精华国产精华精| 亚洲天堂国产精品一区在线| 999久久久国产精品视频| 欧美激情极品国产一区二区三区| 免费在线观看亚洲国产| 亚洲午夜精品一区,二区,三区| 在线观看一区二区三区| 日韩大码丰满熟妇| 日日夜夜操网爽| www.自偷自拍.com| а√天堂www在线а√下载| 成年女人毛片免费观看观看9| 咕卡用的链子| 国产主播在线观看一区二区| 精品欧美国产一区二区三| 日日爽夜夜爽网站| 亚洲国产精品sss在线观看| 男女之事视频高清在线观看| 多毛熟女@视频| 一边摸一边抽搐一进一出视频| 亚洲欧美日韩无卡精品| 国产精品久久视频播放| 啦啦啦观看免费观看视频高清 | 久久这里只有精品19| 亚洲va日本ⅴa欧美va伊人久久| 国产av一区二区精品久久| 成年人黄色毛片网站| 热re99久久国产66热| 一级a爱片免费观看的视频| 久久久久久久精品吃奶| 九色国产91popny在线| 国产一区二区三区视频了| 亚洲av成人av| 亚洲欧美精品综合久久99| 久久精品亚洲熟妇少妇任你| 国产免费av片在线观看野外av| 91国产中文字幕| 一级毛片高清免费大全| 日韩精品免费视频一区二区三区| tocl精华| 又黄又粗又硬又大视频| 宅男免费午夜| 免费在线观看亚洲国产| 精品福利观看| 亚洲精品美女久久av网站| 99香蕉大伊视频| 在线国产一区二区在线| 中文字幕人成人乱码亚洲影| 日韩大码丰满熟妇| 亚洲中文日韩欧美视频| 日韩欧美一区视频在线观看| 久久性视频一级片| 黑人操中国人逼视频| 涩涩av久久男人的天堂| 免费看十八禁软件| 亚洲欧美激情综合另类| 在线播放国产精品三级| 国产精品综合久久久久久久免费 | 成人精品一区二区免费| 亚洲va日本ⅴa欧美va伊人久久| 亚洲国产精品合色在线| 18禁观看日本| 中文字幕人成人乱码亚洲影| 日韩欧美一区二区三区在线观看| 黄片大片在线免费观看| 色综合婷婷激情| 亚洲少妇的诱惑av| 国产伦人伦偷精品视频| 九色国产91popny在线| 久久久久久人人人人人| 午夜亚洲福利在线播放| 黑人欧美特级aaaaaa片| 亚洲国产精品sss在线观看| 久久精品91无色码中文字幕| 丝袜美足系列| 久久狼人影院| 亚洲色图 男人天堂 中文字幕| 国内精品久久久久精免费| 午夜福利高清视频| 最近最新中文字幕大全免费视频| 国产成人av激情在线播放| 丝袜在线中文字幕| 日本免费a在线| 法律面前人人平等表现在哪些方面| 亚洲精品一卡2卡三卡4卡5卡| 亚洲av成人av| e午夜精品久久久久久久| 免费在线观看视频国产中文字幕亚洲| 男女床上黄色一级片免费看| 精品一区二区三区四区五区乱码| 精品一区二区三区视频在线观看免费| 久久久久国产精品人妻aⅴ院| 亚洲av成人av| 成人18禁高潮啪啪吃奶动态图| 男女午夜视频在线观看| 一级作爱视频免费观看| 激情视频va一区二区三区| 亚洲七黄色美女视频| 人人澡人人妻人| 一级a爱片免费观看的视频| 别揉我奶头~嗯~啊~动态视频| 亚洲一区中文字幕在线| 久久中文字幕人妻熟女| 色综合婷婷激情| 午夜福利,免费看| 精品午夜福利视频在线观看一区| 神马国产精品三级电影在线观看 | 亚洲九九香蕉| 亚洲精品av麻豆狂野| 美女免费视频网站| 99久久精品国产亚洲精品| 亚洲中文字幕日韩| 国产精品 国内视频| 午夜日韩欧美国产| 99国产精品免费福利视频| 国产精品久久电影中文字幕| 久久国产乱子伦精品免费另类| 成年人黄色毛片网站| 欧美性长视频在线观看| 亚洲精品美女久久av网站| www.精华液| 人人妻,人人澡人人爽秒播| 中国美女看黄片| 午夜福利影视在线免费观看| 免费高清视频大片| 久久影院123| 美女免费视频网站| 久久精品国产99精品国产亚洲性色 | 亚洲欧美日韩另类电影网站| 亚洲一码二码三码区别大吗| 午夜福利影视在线免费观看| 90打野战视频偷拍视频| 日韩大尺度精品在线看网址 | 国产日韩一区二区三区精品不卡| 精品人妻在线不人妻| 精品一区二区三区四区五区乱码| 午夜视频精品福利| 波多野结衣一区麻豆| 女同久久另类99精品国产91| 久久精品91无色码中文字幕| 黄色成人免费大全| 免费人成视频x8x8入口观看| 91在线观看av| 国产精品精品国产色婷婷| 啦啦啦观看免费观看视频高清 | 色播在线永久视频| 岛国视频午夜一区免费看| 亚洲人成伊人成综合网2020| 日韩国内少妇激情av| 中文字幕人成人乱码亚洲影| a在线观看视频网站| 国产精品野战在线观看| 亚洲少妇的诱惑av| 一个人观看的视频www高清免费观看 | 亚洲九九香蕉| 亚洲免费av在线视频| 老熟妇乱子伦视频在线观看| 99re在线观看精品视频| 99riav亚洲国产免费| 夜夜躁狠狠躁天天躁| 黄色 视频免费看| www.999成人在线观看| 69精品国产乱码久久久| 午夜日韩欧美国产| 美女午夜性视频免费| www.www免费av| 美国免费a级毛片| 亚洲激情在线av| 成年人黄色毛片网站| 亚洲av片天天在线观看| 一a级毛片在线观看| 99在线视频只有这里精品首页| 精品一区二区三区av网在线观看| 久久人人97超碰香蕉20202| 香蕉丝袜av| 国产亚洲欧美精品永久| 9热在线视频观看99| 亚洲成人国产一区在线观看| 真人一进一出gif抽搐免费| 成人精品一区二区免费| 亚洲少妇的诱惑av| 日日干狠狠操夜夜爽| 高清在线国产一区| 国产亚洲av高清不卡| 亚洲精品美女久久久久99蜜臀| 亚洲人成伊人成综合网2020| 黑人巨大精品欧美一区二区mp4| 国产欧美日韩精品亚洲av| 亚洲男人的天堂狠狠| 午夜精品国产一区二区电影| 日本a在线网址| 久久青草综合色| 国产又爽黄色视频| 国产亚洲欧美精品永久| 婷婷六月久久综合丁香| 色尼玛亚洲综合影院| 亚洲国产看品久久| 熟女少妇亚洲综合色aaa.| 99re在线观看精品视频| 午夜成年电影在线免费观看| 精品人妻在线不人妻| 国产成人啪精品午夜网站| 免费少妇av软件| 在线av久久热| 99re在线观看精品视频| 日本在线视频免费播放| 丰满的人妻完整版| 制服人妻中文乱码| 精品乱码久久久久久99久播| 午夜久久久久精精品| 久久精品国产亚洲av香蕉五月| 欧美日韩亚洲国产一区二区在线观看| 中文字幕人成人乱码亚洲影| svipshipincom国产片| 男人舔女人下体高潮全视频| 在线播放国产精品三级| 88av欧美| 丁香六月欧美| 国产av精品麻豆| 午夜日韩欧美国产| 黄频高清免费视频| 久久精品亚洲精品国产色婷小说| 高清黄色对白视频在线免费看| 欧美性长视频在线观看| a级毛片在线看网站| 9热在线视频观看99| 亚洲黑人精品在线| 少妇 在线观看| av超薄肉色丝袜交足视频| 免费无遮挡裸体视频| 人成视频在线观看免费观看| 一级黄色大片毛片| 国产成人影院久久av| 亚洲欧美日韩另类电影网站| 久久午夜亚洲精品久久| 天天躁夜夜躁狠狠躁躁| 如日韩欧美国产精品一区二区三区| 天天一区二区日本电影三级 | 变态另类成人亚洲欧美熟女 | 自线自在国产av| 精品一区二区三区四区五区乱码| 亚洲欧美日韩无卡精品| 老司机福利观看| 一a级毛片在线观看| 欧美丝袜亚洲另类 | 亚洲精品在线观看二区| 国产精品日韩av在线免费观看 | 神马国产精品三级电影在线观看 | 国产精品香港三级国产av潘金莲| 99riav亚洲国产免费| 99精品在免费线老司机午夜| 精品久久久久久,| 欧美成人一区二区免费高清观看 | 亚洲成人久久性| a级毛片在线看网站| 国产亚洲精品av在线| 在线永久观看黄色视频| 老司机午夜福利在线观看视频| 国产xxxxx性猛交| 午夜久久久在线观看| 免费av毛片视频| 淫妇啪啪啪对白视频| 国产一区二区三区综合在线观看| 精品国产国语对白av| 搡老妇女老女人老熟妇| 侵犯人妻中文字幕一二三四区| 国产成人影院久久av| 亚洲一区二区三区色噜噜| 国产欧美日韩综合在线一区二区| 无遮挡黄片免费观看| 亚洲精品国产一区二区精华液| 日韩欧美一区视频在线观看| 精品一区二区三区av网在线观看| 亚洲无线在线观看| 女性被躁到高潮视频| 母亲3免费完整高清在线观看| 久久国产精品影院| 一区在线观看完整版| 久久精品国产清高在天天线| 午夜久久久久精精品| 国产欧美日韩精品亚洲av| 看黄色毛片网站| 叶爱在线成人免费视频播放| 国产精品国产高清国产av| 色精品久久人妻99蜜桃| 国产午夜福利久久久久久| 一进一出抽搐动态| 久久国产乱子伦精品免费另类| 午夜福利成人在线免费观看| 亚洲,欧美精品.| 中文字幕久久专区| 日韩视频一区二区在线观看| 曰老女人黄片| 在线观看免费日韩欧美大片| 神马国产精品三级电影在线观看 | 国产精品秋霞免费鲁丝片| 男女午夜视频在线观看| 久久伊人香网站| 中文字幕人妻丝袜一区二区| 国产成人精品久久二区二区免费| 久久人妻熟女aⅴ| 手机成人av网站| 在线观看午夜福利视频| 男人舔女人的私密视频| 国产单亲对白刺激| 国产一卡二卡三卡精品| 一级毛片高清免费大全| 亚洲av熟女| 久热爱精品视频在线9| 两性夫妻黄色片| 久久欧美精品欧美久久欧美| 美女扒开内裤让男人捅视频| 岛国视频午夜一区免费看| 久久久久久久久中文| 最近最新中文字幕大全免费视频| 夜夜看夜夜爽夜夜摸| 欧美日韩一级在线毛片| 日本欧美视频一区| 国产aⅴ精品一区二区三区波| АⅤ资源中文在线天堂| 午夜福利影视在线免费观看| 久久久久国内视频| 99久久久亚洲精品蜜臀av| 久久久久久亚洲精品国产蜜桃av| 日韩欧美一区视频在线观看| 日日夜夜操网爽| 搡老妇女老女人老熟妇| 亚洲电影在线观看av| 亚洲人成伊人成综合网2020| 精品人妻在线不人妻| 怎么达到女性高潮| 黄频高清免费视频| 在线观看www视频免费| 免费在线观看亚洲国产| 久久精品aⅴ一区二区三区四区| 精品第一国产精品| 国产精品爽爽va在线观看网站 | 国产精品免费视频内射| 99香蕉大伊视频| 搡老岳熟女国产| 欧美久久黑人一区二区| 嫩草影院精品99| 国产精品98久久久久久宅男小说| 麻豆国产av国片精品| 波多野结衣av一区二区av| 成人18禁在线播放| 国产一区二区在线av高清观看| 国产成+人综合+亚洲专区| 久久精品91无色码中文字幕| 自线自在国产av| 亚洲男人的天堂狠狠| 亚洲熟妇熟女久久| 亚洲一区中文字幕在线| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲成人精品中文字幕电影| 正在播放国产对白刺激| 亚洲av成人av| 超碰成人久久| 青草久久国产| 日本五十路高清| svipshipincom国产片| 久久久久久国产a免费观看| 美女国产高潮福利片在线看| 一级作爱视频免费观看| 国产精品香港三级国产av潘金莲| 国产激情欧美一区二区| 亚洲一码二码三码区别大吗| 国产高清视频在线播放一区| 亚洲avbb在线观看| 琪琪午夜伦伦电影理论片6080| 色婷婷久久久亚洲欧美| 99香蕉大伊视频| 国产又色又爽无遮挡免费看| 看免费av毛片| 人成视频在线观看免费观看| 精品电影一区二区在线| 校园春色视频在线观看| 亚洲中文av在线| 国产乱人伦免费视频| 国产国语露脸激情在线看| 午夜福利欧美成人| 亚洲 欧美一区二区三区| 亚洲欧美一区二区三区黑人| 999久久久精品免费观看国产| 夜夜夜夜夜久久久久| 在线免费观看的www视频| 欧美乱码精品一区二区三区| 国产精华一区二区三区| 国产野战对白在线观看| 国产精品一区二区免费欧美| 性色av乱码一区二区三区2| 色哟哟哟哟哟哟| www国产在线视频色| 一二三四社区在线视频社区8| 熟妇人妻久久中文字幕3abv| 国产色视频综合| 精品久久久久久,| 午夜免费鲁丝| 搡老妇女老女人老熟妇| 女人被躁到高潮嗷嗷叫费观| 午夜福利成人在线免费观看| 欧美日韩中文字幕国产精品一区二区三区 | 最好的美女福利视频网| 久久午夜综合久久蜜桃| 黄频高清免费视频| 一进一出抽搐gif免费好疼| 法律面前人人平等表现在哪些方面| videosex国产| 国产视频一区二区在线看| 欧美黄色淫秽网站| 久久精品影院6| 午夜激情av网站| 热99re8久久精品国产| 午夜福利一区二区在线看| 波多野结衣巨乳人妻| www.www免费av| 天天躁狠狠躁夜夜躁狠狠躁| 精品国内亚洲2022精品成人| 一级片免费观看大全| 亚洲少妇的诱惑av| 淫妇啪啪啪对白视频| cao死你这个sao货| 国产精品久久久久久亚洲av鲁大| 久久人妻av系列| 狠狠狠狠99中文字幕| 91字幕亚洲| 欧美国产日韩亚洲一区| 久久久精品欧美日韩精品| 亚洲熟女毛片儿| 精品午夜福利视频在线观看一区| 琪琪午夜伦伦电影理论片6080| 国产精品亚洲一级av第二区| 在线av久久热| 国产一区二区在线av高清观看| 久久国产精品影院| 欧美中文日本在线观看视频| 国产午夜精品久久久久久| 国产精品一区二区精品视频观看| 极品人妻少妇av视频| 午夜福利,免费看| 成人永久免费在线观看视频| 制服诱惑二区| 久久精品国产亚洲av香蕉五月| av电影中文网址| 99精品久久久久人妻精品| 国产精品久久久人人做人人爽| 咕卡用的链子| 久久香蕉激情| 亚洲视频免费观看视频| 无限看片的www在线观看| 久9热在线精品视频| 免费人成视频x8x8入口观看| 国产男靠女视频免费网站| 中文字幕久久专区| 1024视频免费在线观看| 国产在线观看jvid| 国产精品野战在线观看| 男人的好看免费观看在线视频 | 久久久国产欧美日韩av| 黑人欧美特级aaaaaa片| 母亲3免费完整高清在线观看|