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

    基于柵格新安江模型的秦淮河流域洪水模擬及實(shí)時(shí)校正研究

    2021-10-28 11:12:02孫文宇李致家龔珺夫黃鵬年
    人民珠江 2021年10期
    關(guān)鍵詞:新安江秦淮河場(chǎng)次

    孫文宇,王 容,姚 成,何 健,李致家,龔珺夫,黃鵬年

    (1.河海大學(xué)水文水資源學(xué)院,江蘇 南京 210098;2.水利部水文水資源監(jiān)測(cè)預(yù)報(bào)中心,北京 100053;3.江蘇省水文水資源勘測(cè)局,江蘇 南京 210029;4.南京信息工程大學(xué)水文與水資源工程學(xué)院,江蘇 南京 210044)

    新安江模型是20世紀(jì)60年代中國著名水文學(xué)家趙人俊教授及其團(tuán)隊(duì)提出的流域水文模型,是中國工程水文學(xué)發(fā)展過程中重要成果之一[1]。該模型采用蓄滿產(chǎn)流機(jī)制,適用于濕潤(rùn)地區(qū)或半濕潤(rùn)半干旱地區(qū)。李致家等[2]將新安江模型與TOPKAPI模型應(yīng)用于安徽省呈村流域,取得了良好的模擬效果;姚成[3]構(gòu)建了基于柵格的分布式新安江模型并將其應(yīng)用于錢塘江支流密賽流域;Yang等[4]采用分布式新安江模型,對(duì)中國西南喀斯特流域洪水預(yù)報(bào)進(jìn)行研究;大量的實(shí)踐應(yīng)用表明了新安江模型良好的模擬效果及廣泛的適用范圍。

    實(shí)時(shí)校正是中小河流水文模擬及實(shí)時(shí)洪水預(yù)報(bào)中不可或缺的一環(huán)。Liu等[5]將LSTM-KNN法應(yīng)用于中國多個(gè)氣候區(qū)的流域?qū)崟r(shí)洪水預(yù)報(bào)中,取得了良好的模擬效果;徐杰等[6]利用K最近鄰(KNN)算法、誤差自回歸(AR)方法和反饋模擬3種實(shí)時(shí)校正方法來提高新安江模型洪水預(yù)報(bào)精度,成功應(yīng)用于陜西省陳河流域;黃一昕等[7]梳理了洪水預(yù)報(bào)實(shí)時(shí)校正技術(shù)的研究進(jìn)展,并對(duì)其發(fā)展方向進(jìn)行了展望。水文模型與實(shí)時(shí)校正方法相結(jié)合,既遵循流域本身的水文規(guī)律,又能及時(shí)修正人為因素影響帶來的系統(tǒng)性偏差。

    本研究采用新安江模型和柵格新安江模型對(duì)秦淮河上游流域進(jìn)行場(chǎng)次洪水模擬,再利用KNN法和反饋模擬進(jìn)行實(shí)時(shí)校正,取得了良好的次洪模擬結(jié)果,既能給秦淮河水系中下游水力學(xué)模型模擬研究提供上邊界條件,也為秦淮河及相似流域洪水預(yù)報(bào)工作提供參考。

    1 研究區(qū)概況和數(shù)據(jù)來源

    秦淮河位于長(zhǎng)江下游南京段南岸,總體流向自東南往西北,發(fā)源于溧水河、句容河。兩河于方山西南處匯集,北上至河定橋處分流,東向?yàn)榍鼗葱潞?,到大勝關(guān)入江口處匯入長(zhǎng)江;北向稱古秦淮河,古秦淮河又于東水關(guān)分流為內(nèi)秦淮河及外秦淮河,并于西水關(guān)重匯,最終由三汊河入江口匯入長(zhǎng)江。流域地屬北亞熱帶季風(fēng)氣候區(qū),氣候濕潤(rùn),四季分明,常于6、7月份因冷暖氣團(tuán)相遇產(chǎn)生鋒面低壓和靜止鋒,導(dǎo)致連日陰雨,即梅雨期。流域整體呈蒲扇形,長(zhǎng)寬各50 km左右,四周為丘陵山區(qū),腹地為低洼圩區(qū)。

    對(duì)秦淮河水系水文特性及洪水模擬研究于南京市的居民生命財(cái)產(chǎn)安全有重大意義。早在2004年,王建群等[8]便嘗試將流域數(shù)字水文模型應(yīng)用于秦淮河流域,論證了不同土地利用變化對(duì)產(chǎn)匯流的影響;此后,易建軍[9]進(jìn)一步對(duì)秦淮河流域洪水預(yù)報(bào)與調(diào)度模型進(jìn)行研究,將秦淮河劃分為山區(qū)和圩區(qū),分別建立了產(chǎn)匯流模型;因人類活動(dòng)對(duì)秦淮河流域的影響較大,李倩[10]著重分析了城市化空間格局變化的水文效應(yīng),Yuan等[11]分析了城市面積及格局變化對(duì)秦淮河流域洪水過程的影響,袁玉、王懷志等[12-13]將城市群的“防洪圈”概化為圩垸,進(jìn)行了秦淮河流域圩垸式防洪模式洪水模擬。

    由于中下游河段比降極小,下游建有武定門閘和秦淮新河閘,終點(diǎn)處為長(zhǎng)江入江口,故而水位受人工調(diào)度及長(zhǎng)江潮汐頂托影響較大,易出現(xiàn)徑流反向的情況;上游至句容河與溧水河交匯處,受到人工調(diào)度及潮汐頂托的影響相對(duì)較小,故能以交匯處設(shè)立的前垾村(秦淮)水文站為出口點(diǎn)(圖1),采用新安江模型和柵格新安江模型對(duì)出口點(diǎn)以上流域的汛期洪水流量過程進(jìn)行模擬研究[14]。

    圖1 秦淮河流域概況

    研究所用徑流數(shù)據(jù)為秦淮河流域前垾村(秦)站2014—2018年的汛期逐時(shí)徑流量,由《中華人民共和國水文年鑒——長(zhǎng)江流域水文資料》中的實(shí)測(cè)流量成果表計(jì)算得到,摘取了2014年7月到2018年8月汛期內(nèi)的10場(chǎng)洪水進(jìn)行研究。模型中使用的高分辨率數(shù)字高程資料源自地理信息數(shù)據(jù)云網(wǎng)站(http://www.gscould.cn)。

    2 研究方法

    2.1 新安江模型

    新安江模型(XAJ Model)采用三層蒸散發(fā)模式計(jì)算流域蒸散發(fā),產(chǎn)流計(jì)算采用蓄滿產(chǎn)流機(jī)制,匯流計(jì)算分地表徑流、壤中水和地下徑流匯流。主要參數(shù)有蒸散發(fā)折算系數(shù)K、自由水蓄水容量SM、河網(wǎng)水流消退系數(shù)CS等。新安江模型對(duì)流域的劃分較為簡(jiǎn)單,主要以下墊面條件、降水分布情況、站點(diǎn)布設(shè)位置等主要因素,和泰森多邊形法、等雨量線法、反距離加權(quán)法等方法對(duì)流域進(jìn)行分塊處理。新安江模型對(duì)數(shù)據(jù)資料的要求相對(duì)較少,前期準(zhǔn)備和計(jì)算過程都更為便利,在濕潤(rùn)半濕潤(rùn)流域應(yīng)用方便且精度也能得到保證。關(guān)于模型的更多細(xì)節(jié),參見文獻(xiàn)[15]。

    2.2 柵格新安江模型

    柵格新安江模型(Grid-XAJ Model)以流域內(nèi)的每個(gè)DEM柵格作為計(jì)算單元,根據(jù)柵格間的演算次序,將各個(gè)柵格上劃分的地表徑流、壤中流和地下徑流3種水源演算至流域出口。柵格模型考慮了各個(gè)柵格降雨、地貌特征、土壤類型、植被覆蓋等條件的差異性,考慮了柵格間進(jìn)行匯流演算時(shí)上游來水對(duì)蓄滿產(chǎn)流的影響,也考慮了河道降水及河道排水網(wǎng)絡(luò)對(duì)柵格單元產(chǎn)匯流計(jì)算的影響。關(guān)于柵格新安江模型的更多細(xì)節(jié),參見文獻(xiàn)[16]。

    2.3 實(shí)時(shí)校正方法

    根據(jù)實(shí)測(cè)的水位或者流量資料,通過反饋機(jī)制修正預(yù)報(bào)值的過程稱為實(shí)時(shí)校正。中小河流的洪水模擬預(yù)報(bào)易受人為因素?cái)_動(dòng),預(yù)見期相對(duì)較短,具有采用實(shí)時(shí)校正方法的必要性和可行性。本研究選用了KNN法和反饋法2種水文模型中常用的實(shí)時(shí)校正方法。

    2.3.1KNN法

    KNN法是一種歷史匹配方法,通過比較樣本特征向量間的歐氏距離來匹配數(shù)個(gè)近鄰樣本序列,用反距離權(quán)重法對(duì)預(yù)報(bào)序列誤差值進(jìn)行校正(圖2)。KNN法步驟包括歷史樣本庫構(gòu)建、臨近樣本篩選、預(yù)報(bào)值校正和滾動(dòng)校正。主要包含2個(gè)參數(shù):參與誤差估計(jì)的近鄰樣本數(shù)目k、當(dāng)前預(yù)報(bào)誤差有關(guān)聯(lián)的歷史預(yù)報(bào)誤差個(gè)數(shù)s。k取值范圍5~100,取值過小則誤差校正依賴少數(shù)樣本,曲線震蕩,取值過大會(huì)削弱有價(jià)值樣本的作用,曲線過分平滑;s取值范圍2~10,取值過小則僅適用于洪峰等洪水突變部位,取值過大則僅適用于洪水狀態(tài)變化緩慢的部位。

    圖2 KNN法原理示意

    2.3.2反饋法

    反饋法首先計(jì)算實(shí)測(cè)及預(yù)報(bào)序列相鄰差值,再由相鄰2個(gè)時(shí)段的實(shí)測(cè)流量差值之和與預(yù)報(bào)流量差值之和的比值計(jì)算FACT因子,最后根據(jù)預(yù)報(bào)序列相鄰流量差值判斷漲退水段,分別進(jìn)行反饋校正計(jì)算。反饋法更多應(yīng)用細(xì)節(jié),參見文獻(xiàn)[17]。

    實(shí)時(shí)校正的效果與預(yù)見期的選取密切相關(guān)。若預(yù)見期足夠短,實(shí)時(shí)校正進(jìn)行過程中,有著足夠接近時(shí)段的實(shí)測(cè)數(shù)據(jù)和模擬數(shù)據(jù),校正后的模擬值幾乎等于實(shí)測(cè)值,但預(yù)見期必須有足夠的長(zhǎng)度,才能起到預(yù)報(bào)、預(yù)警、輔助決策的作用;選取較長(zhǎng)的預(yù)見期,雖然有利于預(yù)警調(diào)度等工作,但容易導(dǎo)致校正結(jié)果誤差增大,無法及時(shí)響應(yīng)洪峰段過程線的變化,校正后的洪峰會(huì)過分偏離實(shí)測(cè)值。

    3 在秦淮河流域上游段的應(yīng)用

    采用柵格新安江模型及新安江模型程序?qū)η皥ゴ?秦淮)水文站的10場(chǎng)洪水實(shí)測(cè)數(shù)據(jù)進(jìn)行模擬,計(jì)算時(shí)段長(zhǎng)取1 h,其中前7場(chǎng)用于率定,后3場(chǎng)用于檢驗(yàn)。其中柵格新安江模型部分參數(shù)由輸入的土壤類型及植被覆蓋數(shù)據(jù)直接估計(jì),或從物理意義上建立聯(lián)系計(jì)算得出[18],如SM、C、WM、WUM、WLM等,根據(jù)相應(yīng)柵格的下墊面及植被季節(jié)性分布情況有不同取值,柵格之間各有差異。部分參數(shù)假定它們的取值流域空間分布均勻,流域內(nèi)統(tǒng)一賦常用值或估計(jì)值,包括K、CI、CG、CS等。均以最高確定性系數(shù)為目標(biāo),采用SCE-UA算法,自動(dòng)率定得到一組最優(yōu)參數(shù)值。流域內(nèi)統(tǒng)一賦值參數(shù)的率定值分別見表1、2,模型計(jì)算結(jié)果特征值見表3,部分典型場(chǎng)次的模擬結(jié)果見圖3。

    a)20150614場(chǎng)次洪水過程

    表1 柵格新安江模型參數(shù)值

    表2 新安江模型參數(shù)值

    表3 次洪模擬特征值

    從柵格新安江模型的模擬效果來看,10場(chǎng)洪水中有7場(chǎng)洪水確定性系數(shù)大于0.80,其中20140702、20150602、20160701、20170923四場(chǎng)洪峰和洪量相對(duì)誤差小于20%,模擬效果良好[19];20160701場(chǎng)次洪水過程雨量洪量都較大,模擬結(jié)果中的洪峰計(jì)算值過大,結(jié)合實(shí)際情況分析,存在人為調(diào)度等因素,削減了實(shí)際洪峰;20160713場(chǎng)次洪水過程的雨量較小,但該場(chǎng)次暴雨中心主要位于城鎮(zhèn)、工業(yè)園區(qū)較為集中的土橋江寧,城市不透水面占比大,導(dǎo)致洪峰陡漲陡落;20180815場(chǎng)次洪水洪峰起漲前存在短暫的負(fù)流量,受下游感潮河段水力聯(lián)系影響,模擬結(jié)果的洪水過程線整體偏大,增大了洪量和洪峰相對(duì)誤差。

    新安江模型在前垾村流域10場(chǎng)洪水的模擬結(jié)果中,有5場(chǎng)洪水模擬確定性系數(shù)超過0.80,其中20150623場(chǎng)次洪峰、洪量模擬較好,而該場(chǎng)次降雨在實(shí)測(cè)洪峰出現(xiàn)前后各集中發(fā)生1次,模型模擬計(jì)算時(shí)2次降雨洪水過程之間出現(xiàn)短暫退水,導(dǎo)致模擬過程線出現(xiàn)雙峰,峰現(xiàn)時(shí)間誤差增大;模擬結(jié)果偏差的2014702、20150614、20160713三場(chǎng)次洪過程,洪峰相對(duì)誤差均超過40%,需要配合相應(yīng)的實(shí)時(shí)校正系統(tǒng),以實(shí)現(xiàn)一定預(yù)見期內(nèi)更為準(zhǔn)確的洪水流量過程模擬。

    4 實(shí)時(shí)校正

    中下游閘壩的人為調(diào)度因水力聯(lián)系對(duì)上游有顯著影響,也一定程度上為模擬預(yù)報(bào)工作帶來了更多的變化因素。汛期調(diào)度工作的主觀性和多變性增加了預(yù)報(bào)工作的難度。采用實(shí)時(shí)校正程序,可以有效補(bǔ)足原始預(yù)報(bào)模型對(duì)人為因素影響考慮的欠缺。

    采用實(shí)時(shí)校正程序,依次選取1~6 h作為預(yù)見期[20],取KNN法參數(shù)k=8,s=3[21],對(duì)柵格新安江模型和新安江模型模擬的10場(chǎng)次洪進(jìn)行校正。其中6 h預(yù)見期(e=6 h)對(duì)應(yīng)的確定性系數(shù)和峰現(xiàn)時(shí)間誤差見表4,不同預(yù)見期下2種校正方法模擬結(jié)果的洪峰相對(duì)誤差見圖4。

    表4 實(shí)時(shí)校正后確定性系數(shù)及峰現(xiàn)時(shí)間誤差(e=6 h)

    a)20140702場(chǎng)次

    經(jīng)實(shí)時(shí)校正后,2種模型的模擬效果都明顯提高,主要反映在較大的確定性系數(shù)(表4)和總體減小的洪峰相對(duì)誤差(圖4)上,峰現(xiàn)時(shí)間誤差總體也相對(duì)減小,6 h預(yù)見期實(shí)時(shí)校正后,以3 h作為峰現(xiàn)時(shí)間許可誤差,柵格新安江模型峰現(xiàn)時(shí)間誤差合格率達(dá)0.8,新安江模型峰現(xiàn)時(shí)間誤差合格率達(dá)0.6。但類似20160701場(chǎng)次,洪水過程多變,短時(shí)間內(nèi)出現(xiàn)多次洪峰,進(jìn)行預(yù)見期較大的實(shí)時(shí)校正后可能會(huì)進(jìn)一步放大峰現(xiàn)時(shí)間誤差。實(shí)時(shí)校正的效果與預(yù)見期時(shí)長(zhǎng)密切相關(guān),10場(chǎng)次洪的實(shí)時(shí)校正后洪峰相對(duì)誤差及峰現(xiàn)時(shí)間誤差隨著預(yù)見期減小而呈現(xiàn)減小的趨勢(shì),即預(yù)見期越短,實(shí)時(shí)校正的模擬效果就越好。分析實(shí)時(shí)校正后洪峰相對(duì)誤差較大的20150602場(chǎng)次,該場(chǎng)次降雨歷時(shí)短,強(qiáng)度大,洪水漲洪快,在較長(zhǎng)預(yù)見期的實(shí)時(shí)校正過程中,實(shí)際洪峰已達(dá)洪峰處時(shí),校正程序依賴的數(shù)據(jù)系列在相應(yīng)預(yù)見期之前的漲洪段,從而錯(cuò)誤地抬高洪峰,增大洪峰誤差。預(yù)見期縮短時(shí),可利用的數(shù)據(jù)序列更加接近洪峰段,實(shí)時(shí)校正洪峰誤差明顯減??;而新安江模型在該場(chǎng)次的模擬結(jié)果中,本身計(jì)算洪峰偏小,抵消了誤差,所以模擬結(jié)果較好。

    實(shí)時(shí)校正前,柵格新安江模型平均確定性系數(shù)略高,洪峰相對(duì)誤差及峰現(xiàn)時(shí)間誤差較小,模擬效果優(yōu)于新安江模型;實(shí)時(shí)校正后,柵格新安江模型利用2種實(shí)時(shí)校正方法獲得的模擬結(jié)果確定性系數(shù)也較新安江模型更高,但新安江模型實(shí)時(shí)校正后的洪峰誤差較小。

    KNN法和反饋法都能實(shí)現(xiàn)良好的實(shí)時(shí)校正效果,均能有效提高洪水模擬預(yù)報(bào)精度。對(duì)比2種方法,從校正后的確定性系數(shù)、洪峰相對(duì)誤差及流量過程線形態(tài)來看,反饋模擬得到的結(jié)果確定性系數(shù)整體更高,洪峰相對(duì)誤差稍小,且流量過程線更平滑自然。但KNN法能根據(jù)實(shí)際需要調(diào)節(jié)參數(shù)以適應(yīng)工作重點(diǎn),或者通過參數(shù)優(yōu)選的方法確定具體流域的最優(yōu)參數(shù)組,應(yīng)用更為靈活。

    5 結(jié)論及展望

    采用新安江模型和柵格新安江模型,結(jié)合KNN法及反饋模擬2種實(shí)時(shí)校正方法,對(duì)秦淮河上游流域2014—2018年10場(chǎng)次洪過程進(jìn)行模擬研究,主要結(jié)論如下。

    a)在秦淮河上游前垾村流域的洪水模擬工作中,柵格新安江模型及新安江模型經(jīng)由KNN法和反饋法進(jìn)行實(shí)時(shí)校正后,模擬預(yù)報(bào)在1~6 h預(yù)見期內(nèi)均能達(dá)到0.90以上的確定性系數(shù)及25%以內(nèi)的洪峰相對(duì)誤差,峰現(xiàn)時(shí)間誤差合格情況良好,可應(yīng)用于前垾村流域的洪水模擬及預(yù)報(bào)工作。

    b)預(yù)見期選取對(duì)實(shí)時(shí)校正效果影響顯著,縮短預(yù)見期可以使實(shí)時(shí)校正程序獲得更接近的數(shù)據(jù)序列,大大提高實(shí)時(shí)校正的模擬精度;預(yù)見期較大時(shí),實(shí)時(shí)校正會(huì)放大部分漲洪較快場(chǎng)次洪水的洪峰模擬值。

    c)城市不透水面積增大、暴雨中心集中在城鎮(zhèn)區(qū)域?qū)涤旰樗^程存在影響,導(dǎo)致產(chǎn)匯流速度加快,洪峰陡漲陡落,峰形尖瘦,致災(zāi)風(fēng)險(xiǎn)高。

    基于對(duì)研究結(jié)果的分析和思考,城市地區(qū)不透水區(qū)域、地下管網(wǎng)排水系統(tǒng)等因素對(duì)徑流過程存在的影響是后續(xù)研究工作的目標(biāo)。在新安江模型中添加相應(yīng)的計(jì)算模塊進(jìn)行表達(dá),如城市管網(wǎng)排水系統(tǒng)產(chǎn)匯流調(diào)節(jié)模塊、暴雨中心位置權(quán)重分析模塊等,對(duì)秦淮河及相似中小河流流域水文模擬預(yù)報(bào)研究具有一定實(shí)用價(jià)值。

    猜你喜歡
    新安江秦淮河場(chǎng)次
    長(zhǎng)江上游高洪水期泥沙輸移特性
    基于運(yùn)行場(chǎng)次用時(shí)誤差的載人設(shè)備故障預(yù)警可視化研究
    排考場(chǎng)次分配方法及其SQL實(shí)現(xiàn)
    夜航
    新安江霧
    照相機(jī)(2017年10期)2017-11-22 05:51:11
    秦淮河水冷,戰(zhàn)事幾回傷
    中國三峽(2016年11期)2017-01-15 14:00:08
    南京城與秦淮河
    中國三峽(2016年11期)2017-01-15 14:00:04
    新安江上游橫江地表水中砷、硒、汞含量變化研究
    新安江上游流域農(nóng)村污水處理技術(shù)應(yīng)用
    地鐵觀影指南
    電影故事(2015年33期)2015-09-06 01:05:30
    av福利片在线观看| 色噜噜av男人的天堂激情| 成人高潮视频无遮挡免费网站| 黄色一级大片看看| 麻豆久久精品国产亚洲av| 动漫黄色视频在线观看| 婷婷精品国产亚洲av| 久久精品91蜜桃| 久久欧美精品欧美久久欧美| 欧美日韩精品成人综合77777| 久久人人精品亚洲av| 久久久久久伊人网av| 18禁黄网站禁片免费观看直播| 在线观看一区二区三区| 亚洲精品成人久久久久久| 精品久久久久久,| 精品一区二区免费观看| 99在线人妻在线中文字幕| 免费高清视频大片| 一区二区三区免费毛片| x7x7x7水蜜桃| 午夜精品一区二区三区免费看| 久久久久久久亚洲中文字幕| 午夜福利成人在线免费观看| 欧美国产日韩亚洲一区| 联通29元200g的流量卡| 午夜免费男女啪啪视频观看 | 成人精品一区二区免费| 中国美白少妇内射xxxbb| 国产精品1区2区在线观看.| 99久久久亚洲精品蜜臀av| 国产精品福利在线免费观看| 婷婷六月久久综合丁香| 亚洲内射少妇av| 最近视频中文字幕2019在线8| 免费搜索国产男女视频| 国产人妻一区二区三区在| 99久国产av精品| 欧美一级a爱片免费观看看| 国产精品野战在线观看| 国产伦一二天堂av在线观看| 色综合色国产| 三级国产精品欧美在线观看| 亚洲成人中文字幕在线播放| 丰满的人妻完整版| 床上黄色一级片| 美女xxoo啪啪120秒动态图| 欧美性猛交╳xxx乱大交人| 亚洲精品影视一区二区三区av| 一本一本综合久久| 女人被狂操c到高潮| 在线播放国产精品三级| 国产精品野战在线观看| 婷婷亚洲欧美| 成人av一区二区三区在线看| 18禁黄网站禁片免费观看直播| 99热只有精品国产| 一级毛片久久久久久久久女| 久久久色成人| 午夜福利成人在线免费观看| 午夜视频国产福利| 女同久久另类99精品国产91| 免费看a级黄色片| 国产综合懂色| 日本免费a在线| 欧美xxxx性猛交bbbb| 国内精品久久久久精免费| .国产精品久久| 美女大奶头视频| 久久久精品欧美日韩精品| 伦精品一区二区三区| 老师上课跳d突然被开到最大视频| 国产精品99久久久久久久久| 日韩中文字幕欧美一区二区| 久久热精品热| 黄片wwwwww| 免费电影在线观看免费观看| 欧美日韩中文字幕国产精品一区二区三区| 91在线精品国自产拍蜜月| 桃红色精品国产亚洲av| 一级av片app| 99久久成人亚洲精品观看| 日日摸夜夜添夜夜添小说| 老熟妇乱子伦视频在线观看| 午夜视频国产福利| 天天一区二区日本电影三级| 国产精品不卡视频一区二区| 精品福利观看| 在线观看av片永久免费下载| 日韩精品有码人妻一区| 97热精品久久久久久| 天堂影院成人在线观看| 91久久精品电影网| 一夜夜www| 精品午夜福利视频在线观看一区| 欧美三级亚洲精品| 少妇的逼好多水| 欧美xxxx性猛交bbbb| 在线观看舔阴道视频| 午夜影院日韩av| 日本爱情动作片www.在线观看 | 精品人妻一区二区三区麻豆 | 婷婷精品国产亚洲av在线| 极品教师在线免费播放| 亚洲久久久久久中文字幕| 亚洲av二区三区四区| 99热网站在线观看| 最近在线观看免费完整版| 日韩一区二区视频免费看| 欧美又色又爽又黄视频| 偷拍熟女少妇极品色| 极品教师在线视频| 欧美最新免费一区二区三区| 精品久久国产蜜桃| 内地一区二区视频在线| 中文亚洲av片在线观看爽| 久久久久久久久久黄片| 免费在线观看成人毛片| 老熟妇乱子伦视频在线观看| 99久久中文字幕三级久久日本| 久久九九热精品免费| 国产精品电影一区二区三区| 欧美日韩精品成人综合77777| 国产精品久久视频播放| 国产精品电影一区二区三区| 久久精品国产自在天天线| 精品久久久久久久末码| 看片在线看免费视频| 亚洲乱码一区二区免费版| 久久6这里有精品| 国产av一区在线观看免费| 夜夜看夜夜爽夜夜摸| 亚洲熟妇熟女久久| 如何舔出高潮| 久久香蕉精品热| 少妇的逼好多水| 国产亚洲欧美98| 黄色一级大片看看| 88av欧美| 蜜桃亚洲精品一区二区三区| 身体一侧抽搐| 国产精品爽爽va在线观看网站| 老司机午夜福利在线观看视频| 国产精品98久久久久久宅男小说| 国国产精品蜜臀av免费| 久久久国产成人精品二区| 国产乱人伦免费视频| 亚洲一级一片aⅴ在线观看| 一进一出抽搐动态| 一个人看的www免费观看视频| 国产私拍福利视频在线观看| 亚洲欧美激情综合另类| 中文亚洲av片在线观看爽| 欧美一区二区国产精品久久精品| 免费黄网站久久成人精品| 夜夜看夜夜爽夜夜摸| a级毛片免费高清观看在线播放| 中出人妻视频一区二区| 精品不卡国产一区二区三区| 九九热线精品视视频播放| 成年女人毛片免费观看观看9| 小说图片视频综合网站| 欧美一区二区国产精品久久精品| 亚洲美女搞黄在线观看 | 亚洲av成人av| 夜夜夜夜夜久久久久| 白带黄色成豆腐渣| 欧美一区二区国产精品久久精品| 久久久久精品国产欧美久久久| 国产免费一级a男人的天堂| 国内毛片毛片毛片毛片毛片| 午夜精品一区二区三区免费看| 免费看a级黄色片| 中国美女看黄片| 国产欧美日韩精品一区二区| 久久亚洲精品不卡| 麻豆国产av国片精品| 中国美白少妇内射xxxbb| 久久婷婷人人爽人人干人人爱| 真人一进一出gif抽搐免费| 国产精品不卡视频一区二区| 天堂动漫精品| 成人亚洲精品av一区二区| 日本黄大片高清| 制服丝袜大香蕉在线| 又爽又黄a免费视频| 亚洲专区中文字幕在线| 欧美在线一区亚洲| 成人综合一区亚洲| 在线观看一区二区三区| 久久精品国产亚洲网站| av在线蜜桃| 日本爱情动作片www.在线观看 | 免费在线观看影片大全网站| 男女视频在线观看网站免费| 日日摸夜夜添夜夜添av毛片 | 国产激情偷乱视频一区二区| www日本黄色视频网| 男人的好看免费观看在线视频| 18+在线观看网站| 国产高清激情床上av| 一进一出抽搐动态| 在线观看66精品国产| 欧美一区二区国产精品久久精品| 欧洲精品卡2卡3卡4卡5卡区| 最好的美女福利视频网| 午夜精品久久久久久毛片777| 日韩,欧美,国产一区二区三区 | 婷婷色综合大香蕉| bbb黄色大片| 国产乱人伦免费视频| 欧美一区二区国产精品久久精品| 国产精品永久免费网站| 午夜免费成人在线视频| 真人一进一出gif抽搐免费| 免费看光身美女| 天堂动漫精品| av黄色大香蕉| 欧美bdsm另类| 很黄的视频免费| .国产精品久久| 嫩草影院精品99| 如何舔出高潮| 日本成人三级电影网站| 亚洲乱码一区二区免费版| 九色国产91popny在线| 欧美激情国产日韩精品一区| 午夜福利在线在线| 国内毛片毛片毛片毛片毛片| 一夜夜www| 韩国av在线不卡| 国产中年淑女户外野战色| 久久久久久久亚洲中文字幕| 日本成人三级电影网站| 草草在线视频免费看| 一个人观看的视频www高清免费观看| 一级av片app| 伦理电影大哥的女人| 少妇丰满av| 亚洲精华国产精华精| 久久精品国产鲁丝片午夜精品 | 国产一区二区三区在线臀色熟女| 国内精品久久久久久久电影| 日韩欧美在线乱码| 人人妻人人澡欧美一区二区| 欧美一区二区亚洲| 成年女人毛片免费观看观看9| 国内精品一区二区在线观看| 色在线成人网| 美女cb高潮喷水在线观看| 天天一区二区日本电影三级| 国产一区二区三区在线臀色熟女| 99久久精品一区二区三区| 久久久成人免费电影| 很黄的视频免费| 99视频精品全部免费 在线| 超碰av人人做人人爽久久| 亚洲在线观看片| 国产主播在线观看一区二区| 99久国产av精品| 日韩一本色道免费dvd| 午夜福利在线观看吧| 尾随美女入室| 国产一区二区亚洲精品在线观看| av黄色大香蕉| 亚洲美女视频黄频| 一级av片app| 在线观看av片永久免费下载| 最近在线观看免费完整版| 黄色日韩在线| 女人十人毛片免费观看3o分钟| 欧美最黄视频在线播放免费| 日韩欧美三级三区| 午夜久久久久精精品| 美女免费视频网站| 亚洲成a人片在线一区二区| 精品一区二区三区人妻视频| 18禁裸乳无遮挡免费网站照片| 国产伦一二天堂av在线观看| 听说在线观看完整版免费高清| 99热这里只有是精品50| 国产精品,欧美在线| 极品教师在线视频| 日本色播在线视频| 美女高潮的动态| 悠悠久久av| 18禁裸乳无遮挡免费网站照片| 成人国产综合亚洲| 搡老妇女老女人老熟妇| 亚洲国产日韩欧美精品在线观看| 男女边吃奶边做爰视频| 中文字幕久久专区| 97超视频在线观看视频| 久久久久精品国产欧美久久久| 五月伊人婷婷丁香| 婷婷丁香在线五月| 亚洲精品成人久久久久久| 欧美色欧美亚洲另类二区| 国产爱豆传媒在线观看| 久久久国产成人免费| 麻豆成人av在线观看| 狠狠狠狠99中文字幕| 免费av观看视频| 国产精品三级大全| 天堂av国产一区二区熟女人妻| 97超视频在线观看视频| 国产精品一区二区免费欧美| 国产私拍福利视频在线观看| 中文在线观看免费www的网站| 久久精品国产自在天天线| 美女高潮的动态| 男人舔奶头视频| 真人一进一出gif抽搐免费| 国产熟女欧美一区二区| 欧美成人a在线观看| 欧美精品啪啪一区二区三区| 亚洲欧美日韩卡通动漫| 黄色配什么色好看| 免费大片18禁| 老司机深夜福利视频在线观看| 噜噜噜噜噜久久久久久91| 日日夜夜操网爽| 免费不卡的大黄色大毛片视频在线观看 | 免费av毛片视频| 精品久久久久久成人av| 亚洲美女视频黄频| 高清在线国产一区| 我的女老师完整版在线观看| 午夜久久久久精精品| 黄色女人牲交| 黄色丝袜av网址大全| eeuss影院久久| 色视频www国产| 久久精品人妻少妇| 51国产日韩欧美| 国产精品综合久久久久久久免费| 国产成人aa在线观看| 永久网站在线| 国产真实乱freesex| 三级国产精品欧美在线观看| 99精品久久久久人妻精品| 女同久久另类99精品国产91| 午夜福利视频1000在线观看| 国产一级毛片七仙女欲春2| 中亚洲国语对白在线视频| 欧美日韩瑟瑟在线播放| 欧美激情在线99| 欧美激情久久久久久爽电影| 久久久久久久精品吃奶| 亚洲经典国产精华液单| 色精品久久人妻99蜜桃| 国产熟女欧美一区二区| 亚洲av免费高清在线观看| 午夜福利在线观看免费完整高清在 | 婷婷色综合大香蕉| 免费黄网站久久成人精品| 色噜噜av男人的天堂激情| 亚洲av不卡在线观看| 日本色播在线视频| 精品99又大又爽又粗少妇毛片 | 亚洲三级黄色毛片| 亚洲国产欧美人成| 波多野结衣巨乳人妻| 色综合站精品国产| 国产午夜精品论理片| 亚洲最大成人av| 亚洲成av人片在线播放无| 中文字幕熟女人妻在线| 欧美+日韩+精品| 男女之事视频高清在线观看| 国产伦在线观看视频一区| 乱人视频在线观看| 免费看日本二区| 国产午夜精品久久久久久一区二区三区 | 亚洲欧美清纯卡通| 亚洲中文字幕一区二区三区有码在线看| 亚洲va在线va天堂va国产| 亚洲专区国产一区二区| 99久久精品国产国产毛片| 欧美日韩黄片免| 国产亚洲av嫩草精品影院| 免费观看在线日韩| 99在线视频只有这里精品首页| 不卡一级毛片| 热99re8久久精品国产| 午夜福利在线在线| 中出人妻视频一区二区| 国产伦精品一区二区三区视频9| 亚洲第一电影网av| 国产精品免费一区二区三区在线| 中文字幕久久专区| 中文字幕熟女人妻在线| 色5月婷婷丁香| 嫩草影院新地址| 老司机深夜福利视频在线观看| 亚洲七黄色美女视频| 成人高潮视频无遮挡免费网站| 69av精品久久久久久| 国产乱人伦免费视频| 午夜福利在线观看吧| 熟女电影av网| 久久国内精品自在自线图片| 亚洲专区中文字幕在线| 啦啦啦啦在线视频资源| 精品久久久久久久久亚洲 | 一进一出好大好爽视频| av黄色大香蕉| 少妇的逼好多水| 午夜免费男女啪啪视频观看 | 嫁个100分男人电影在线观看| 日韩精品中文字幕看吧| 长腿黑丝高跟| 国内久久婷婷六月综合欲色啪| 日本与韩国留学比较| 国产精品嫩草影院av在线观看 | 最好的美女福利视频网| 日韩高清综合在线| 黄色丝袜av网址大全| 国产精品福利在线免费观看| 美女cb高潮喷水在线观看| 免费看光身美女| 亚洲av日韩精品久久久久久密| 久久国产精品人妻蜜桃| www日本黄色视频网| 成人性生交大片免费视频hd| 内射极品少妇av片p| 久久精品影院6| 又爽又黄a免费视频| 亚洲成人中文字幕在线播放| 51国产日韩欧美| 欧美日韩黄片免| 精品久久久久久久末码| 特级一级黄色大片| 国产一区二区亚洲精品在线观看| 日本在线视频免费播放| 久久国产乱子免费精品| 国产中年淑女户外野战色| 亚洲成人免费电影在线观看| 两个人视频免费观看高清| 91狼人影院| 99国产精品一区二区蜜桃av| 亚洲第一区二区三区不卡| av女优亚洲男人天堂| 欧美+亚洲+日韩+国产| 久久久国产成人精品二区| 热99在线观看视频| 又黄又爽又刺激的免费视频.| 国产成人福利小说| 欧美高清成人免费视频www| 日韩欧美免费精品| 精品福利观看| 99在线视频只有这里精品首页| 午夜福利在线在线| 亚洲久久久久久中文字幕| 97超视频在线观看视频| 午夜日韩欧美国产| 一个人观看的视频www高清免费观看| 亚洲成人免费电影在线观看| 无人区码免费观看不卡| 久久精品国产自在天天线| 老女人水多毛片| 国产一区二区激情短视频| 成人av一区二区三区在线看| 免费观看的影片在线观看| 日韩一本色道免费dvd| 免费看日本二区| 婷婷六月久久综合丁香| 99精品在免费线老司机午夜| 极品教师在线免费播放| 在线观看66精品国产| 麻豆成人午夜福利视频| 美女大奶头视频| 日韩欧美国产在线观看| 欧美日韩中文字幕国产精品一区二区三区| 亚洲最大成人中文| 久久国内精品自在自线图片| av在线观看视频网站免费| 九九爱精品视频在线观看| 免费不卡的大黄色大毛片视频在线观看 | 中国美白少妇内射xxxbb| 色播亚洲综合网| 亚洲一区二区三区色噜噜| 97超级碰碰碰精品色视频在线观看| avwww免费| 国产色婷婷99| 成人特级av手机在线观看| 欧美日韩瑟瑟在线播放| 国产精品嫩草影院av在线观看 | 午夜爱爱视频在线播放| 精华霜和精华液先用哪个| 亚洲不卡免费看| 亚洲国产欧洲综合997久久,| 亚洲,欧美,日韩| 99riav亚洲国产免费| 欧美+亚洲+日韩+国产| 一卡2卡三卡四卡精品乱码亚洲| 99热只有精品国产| 亚洲性久久影院| 色精品久久人妻99蜜桃| 亚洲 国产 在线| 一卡2卡三卡四卡精品乱码亚洲| 亚洲18禁久久av| 国产成人一区二区在线| 欧美成人一区二区免费高清观看| 男女啪啪激烈高潮av片| 日韩欧美在线二视频| 九色成人免费人妻av| 天天一区二区日本电影三级| 久久午夜亚洲精品久久| 午夜福利18| 老女人水多毛片| 亚洲性夜色夜夜综合| 在线观看66精品国产| 51国产日韩欧美| 国产高清激情床上av| 美女被艹到高潮喷水动态| 一区二区三区四区激情视频 | 中国美女看黄片| 久久久久性生活片| 悠悠久久av| 简卡轻食公司| 成人毛片a级毛片在线播放| 一本久久中文字幕| 欧美三级亚洲精品| 国产精品嫩草影院av在线观看 | 亚洲综合色惰| 观看免费一级毛片| 99在线人妻在线中文字幕| a级一级毛片免费在线观看| 国国产精品蜜臀av免费| 91在线观看av| 91午夜精品亚洲一区二区三区 | 欧美绝顶高潮抽搐喷水| 日本爱情动作片www.在线观看 | 国产精品国产高清国产av| 狠狠狠狠99中文字幕| 亚洲最大成人中文| 校园人妻丝袜中文字幕| 国产免费男女视频| 岛国在线免费视频观看| h日本视频在线播放| 高清毛片免费观看视频网站| 尤物成人国产欧美一区二区三区| 黄片wwwwww| 直男gayav资源| 99热只有精品国产| 99精品在免费线老司机午夜| 久久久久久久午夜电影| 欧美日本视频| 欧美日韩乱码在线| 看免费成人av毛片| 中国美女看黄片| 999久久久精品免费观看国产| 成人性生交大片免费视频hd| 日日啪夜夜撸| 两人在一起打扑克的视频| 日韩欧美精品v在线| 蜜桃亚洲精品一区二区三区| 男人舔女人下体高潮全视频| 又爽又黄a免费视频| 精品午夜福利在线看| 色综合婷婷激情| www.www免费av| 亚洲真实伦在线观看| 久久久久久久亚洲中文字幕| 岛国在线免费视频观看| 男女啪啪激烈高潮av片| 午夜久久久久精精品| 少妇裸体淫交视频免费看高清| 亚洲第一区二区三区不卡| 特大巨黑吊av在线直播| 女生性感内裤真人,穿戴方法视频| 亚洲av二区三区四区| 国产日本99.免费观看| 国产精品免费一区二区三区在线| 日韩国内少妇激情av| 人妻少妇偷人精品九色| 午夜福利欧美成人| 国产又黄又爽又无遮挡在线| 国产精品av视频在线免费观看| 女生性感内裤真人,穿戴方法视频| 中文字幕av在线有码专区| 一夜夜www| 波野结衣二区三区在线| 少妇裸体淫交视频免费看高清| 熟女电影av网| 99热6这里只有精品| 草草在线视频免费看| 日本一本二区三区精品| 国产精品爽爽va在线观看网站| 亚洲av五月六月丁香网| 久9热在线精品视频| 亚洲乱码一区二区免费版| 国产av一区在线观看免费| 成年女人永久免费观看视频| 国产毛片a区久久久久| 69人妻影院| av.在线天堂| x7x7x7水蜜桃| 日韩精品有码人妻一区| 欧美潮喷喷水| 日本爱情动作片www.在线观看 | 搞女人的毛片| 男人的好看免费观看在线视频| 亚洲 国产 在线| 亚洲avbb在线观看| 少妇的逼好多水| 少妇熟女aⅴ在线视频| 国内精品宾馆在线| 亚洲人成网站在线播| 成人国产麻豆网| 日日夜夜操网爽| av视频在线观看入口| 国产精品一区二区性色av| 国产精品伦人一区二区|