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

    北洛河流域TOPMODEL模型基流分割定量評估

    2020-10-12 14:27:14王萌萌魏占民張瑞強
    人民黃河 2020年4期
    關(guān)鍵詞:數(shù)字濾波基流徑流

    王萌萌 魏占民 張瑞強

    摘 要:為探究水文模擬結(jié)果對總徑流、地表徑流和地下徑流(基流)的影響,首先采用Mann-Kendall檢驗方法識別出頭水文站水文序列突變點,然后構(gòu)建北洛河流域TOPMODEL水文模型,最后采用數(shù)字濾波方法分別

    對模擬水文過程和實測水文過程進行基流分割,并定量分析水文模型模擬結(jié)果對地表徑流和地下徑流的影響。結(jié)果表明:①北洛河流域徑流突變年份為1994年;②率定期和驗證期納什效率系數(shù)、可決系數(shù)和相對誤差絕對值分別為0.69與0.63、0.87與0.82以及18.5%與13.4%;③率定期基流指數(shù)大于驗證期且大于0.75,這意味著北洛河流域基流量對河川徑流的補給作用顯著;④模型在非汛期徑流的模擬精度優(yōu)于汛期的,能夠為非汛期的基流分割提供依據(jù);⑤不同時期基流量水文過程線起伏差異較大,但其形狀基本與徑流過程線保持一致。

    關(guān)鍵詞:突變檢驗;TOPMODEL水文模型;基流分割;北洛河流域

    中圖分類號:P333 ? 文獻標志碼:A

    doi:10.3969/j.issn.1000-1379.2020.04.004

    Abstract: In order to explore the influence of hydrological simulation results to the runoff, surface runoff and underground runoff (basic flow), the TOPMODEL model of Beiluo River basin was constructed based on the identification of the locations of change points by Mann-Kendall mutation test method in Zhuangtou Hydrological Station. The basic flow segmentation of simulated and measured values by digital filtering method could quantitatively analyze the influence of hydrological model to the surface runoff and underground runoff simulation. The results show that a) the runoff sudden change of Beiluo River basin appeared in 1994. b) Nash efficiency coefficient (NSE), correlation coefficient (R2) and relative error (Re) are 0.69 and 0.63, 0.87 and 0.82, and 18.5% and 13.4% respectively in the calibration and validation periods. c) Basic flow index (BFI) in the calibration periods is greater than the validation period and is greater than 0.75, which means that the underground runoff in the Beiluo River basin has a significant effect on the recharge of river runoff. d) The simulation accuracy of the model in the non-flood period is better than that in the flood season, which can provide a basis for the non-flood period base flow segmentation. e) The fluctuation degree of basic flow process line varies greatly in different periods, but its shape is basically consistent with the runoff process line.

    Key words: mutation test; TOPMODEL hydrological model; basic flow segmentation; Beiluo River basin

    1 引 言

    河川基流是指由地下水補給河川的水量,是枯水季節(jié)河川徑流的主要來源,在維持河川生態(tài)系統(tǒng)健康中具有重要作用[1]。同時,科學(xué)合理地分割流域徑流過程,對流域水資源規(guī)劃、水資源分析計算、水資源配置等具有重要的現(xiàn)實依據(jù)和理論基礎(chǔ)。此外,在干旱半干旱流域,基流對于維持水源的穩(wěn)定性與持續(xù)性具有重要作用。近年來,國內(nèi)外學(xué)者針對基流分割問題開展了大量研究,并獲得了良好成果[2-8]。李瑞等[2]以黃河上游藥水河流域為研究對象,采用兩種基流分割方法(滑動最小值方法和數(shù)字濾波方法)對典型年徑流序列進行了基流分割,結(jié)果表明滑動最小值方法較數(shù)字濾波方法的基流分割值偏小,但均大于全國平均值(0.26);豆林等[3]采用3種基流分割算法(PART法、數(shù)字濾波法和滑動最小值法)對黃土高原區(qū)6個子流域的多年徑流進行了基流分割,結(jié)果表明3種分割方法獲取的基流分割值差異顯著,其中數(shù)字濾波法計算結(jié)果與流域?qū)嶋H徑流過程相一致;焦瑋等[4]以錫林河流域控制站——錫林浩特水文站1963—2012年實測徑流量為基礎(chǔ)數(shù)據(jù),采用數(shù)字濾波法、HYSEP法、滑動最小值法以及加里寧法進行了基流分割,并評估其適應(yīng)性,結(jié)果表明不同算法估算的年基流指數(shù)差異顯著,其中數(shù)字濾波算法和HYSEP法計算結(jié)果較好,最接近流域地下徑流和地表徑流過程,且濾波參數(shù)取0.85時分割的基流結(jié)果最好;黃文等[5]對數(shù)字濾波法、滑動最小值法和時間步長法三類共9種基流分割方法在黑河上游的適用性進行了分析,結(jié)果表明時間步長法中的局部最小值法最適合在黑河上游流域進行基流分割;高東東等[6]采用不同的基流分割方法,對貢嘎山黃崩溜溝流域的徑流過程進行了分割,并對比分析其結(jié)果的適應(yīng)性,得出采用數(shù)字濾波法能夠得到穩(wěn)定的基流分割結(jié)果,同時豐水期地表徑流量相對降水的滯后時間較枯水期短,而基流量相對降水的滯后時間在豐水期、平水期和枯水期差異不明顯;Tularam G. A.等[7]采用不同基流分割方法分割基流和地表徑流,并以最適合的基流分割方法(指數(shù)平滑方法)為例,采用連續(xù)損失、比例損失以及體積徑流系數(shù)3個指標作為基流分割結(jié)果評價指標,定量評估基流分割損失,結(jié)果表明體積徑流系數(shù)適用于不同子流域基流分割損失的評估,且體積徑流系數(shù)的誤差對連續(xù)損失和比例損失無影響;Zhang J.等[8]運用4種不同非示蹤劑基流分離方法,對澳大利亞東部的5個子流域進行了基流分割并對比分析其結(jié)果,得出不同模型參數(shù)的選取對數(shù)字濾波方法基流分割的結(jié)果具有顯著影響,其中運用自動基流分割方法估算的衰退系數(shù)優(yōu)于特定參數(shù)取值,年徑流序列絕對偏差平均減小20%,與不同的參數(shù)化方案相比,不同非示蹤劑基流分離方法之間的差異較小。

    流域水文模型是模擬和分析流域水文過程的重要工具,在特定流域構(gòu)建分布式或概念式水文模型能夠較好地反映其水文變化特征[9-10]。楊軍軍等[9]構(gòu)建了湟水流域SWAT模型并進行了參數(shù)不確定性分析,結(jié)果表明構(gòu)建的SWAT模型對水文過程的模擬精度較高,參數(shù)的不確定性是導(dǎo)致模型率定結(jié)果精度不一的主要原因;宮興龍等[10]通過在半干旱區(qū)的洮兒河流域建立TOPMODEL模型,模擬了各水文要素的時空變化特征,結(jié)果表明由次洪模擬徑流相對誤差絕對值的平均值為 4.45%、洪峰流量相對誤差絕對值的平均值為 5.00%、確定性系數(shù)為 0.76,意味著該模型能夠在干旱區(qū)得到很好的應(yīng)用。綜上可知,數(shù)字濾波分割法能夠較好地分割基流和地表徑流,且不同的水文模型在特定流域均能夠較好地反映流域的水文過程,但未從側(cè)面評估水文模型模擬結(jié)果對基流分割的影響。本文以干旱半干旱區(qū)的北洛河流域為研究對象,為真實反映流域的水文過程特征,以流域控制站——頭水文站1960—2010年徑流資料為基礎(chǔ)數(shù)據(jù),采用Mann-Kendell突變檢驗方法和有序聚類方法識別出徑流序列的突變點,其次構(gòu)建北洛河流域TOPMODEL水文模型,采用遺傳算法(GA)對模型參數(shù)進行率定并模擬流域近似的水文過程;為了評估模擬值與實測值的差異,采用數(shù)字濾波算法對模擬水文過程和實測水文過程進行基流分割,最后定量評估水文模型模擬結(jié)果對總徑流、地表徑流和基流的影響。

    2 流域概況

    北洛河發(fā)源于陜西省定邊縣白于山南麓的草梁山脈,自西北至東南流入渭河。流域面積約為2.69萬km2,流域?qū)儆跍貛О敫珊岛痛箨懶约撅L(fēng)氣候區(qū),四季分明,春季干燥、夏季旱澇相間、秋季濕潤、冬季寒冷干燥。降水主要集中在汛期(6—9月),占全年降水量的76%左右,年平均降雨量為440.3 mm,降水空間分布差異顯著,自東南向西北遞減,年均氣溫為7.5 ℃。干流全長680 km。流域地理位置及空間分布見圖1。

    4 結(jié)果與分析

    4.1 徑流突變檢驗結(jié)果

    Mann-Kendall方法因?qū)λ鶛z驗樣本不必服從特定的頻率分布且不受少量變異值的影響,故被廣泛用于降水、徑流、蒸發(fā)等水文氣象序列的突變檢驗[20]。本文為了劃定流域徑流的天然期(無人類活動影響)和變異期(人類活動影響),首先采用Mann-Kendall方法對頭水文站1960—2010年年徑流序列進行突變點檢驗,結(jié)果見圖2(a),可知徑流突變年份出現(xiàn)在1994年、1996年以及2004年,為了進一步確定突變點位置,使用有序聚類方法對徑流序列進行進一步檢驗,結(jié)果見圖2(b),可知徑流突變年份為1994年。由此可以認為,1994年后北洛河徑流受人類活動影響較大,年尺度徑流模擬過程受人類活動影響顯著,不利于水文模型構(gòu)建。

    4.2 TOPMODEL模型徑流模擬結(jié)果

    4.2.1 子流域地形指數(shù)分布特征與模型參數(shù)的確定

    圖3為北洛河流域不同子流域的地形指數(shù)—面積分布。由圖3可知,各子流域地形指數(shù)—面積分布特征差異較小,均呈現(xiàn)出“幾”字形,當?shù)匦沃笖?shù)值為7.88時,各子流域的面積占比分別為45.30%、42.40%、42.20%、44.90%、46.30%,意味著不同子流域土壤蓄水能力差異不明顯,受人類活動影響較小。

    依據(jù)頭水文站徑流突變點檢測結(jié)果(1994 年為突變年),選取 1986—1988年為模型率定期,1989—1990年為模型驗證期。采用遺傳算法(GA)對模型參數(shù)進行率定,各子流域模型參數(shù)率定結(jié)果見表1(SRmax、SR0分別為根系帶最大蓄水量、根系帶初始含水量;Td為重力排水的時間滯時系數(shù);Ry為河道匯流的有效速度;CHV為河網(wǎng)寬度函數(shù)的有效地表匯流速率)。

    4.2.2 模型率定期和驗證期模擬結(jié)果

    表2為北洛河流域模型率定期和驗證期適應(yīng)性評價結(jié)果,由表2可知,率定期和驗證期納什效率系數(shù)(NSE)、可決系數(shù)(R2)和相對誤差絕對值(|Re|)分別為0.69與0.63、0.87與0.82、18.5%與13.4%,表明率定期和驗證期徑流模擬結(jié)果較好,能夠真實反映流域水文過程,可用于北洛河流域模擬結(jié)果對總徑流、地表徑流和基流影響的定量評估。

    由圖4可知,率定期和驗證期模型模擬值與實測日徑流值擬合效果較好,從不同擬合效果可以看出,小流量時模型能夠很好地捕捉到流量過程,而大流量時捕捉結(jié)果較差,但總體來說徑流模擬結(jié)果較優(yōu)。

    4.2.3 汛期和非汛期水文過程模擬結(jié)果

    為了進一步分析年內(nèi)汛期(6—9月)和非汛期(10月—次年5月)徑流模擬的精度,采用數(shù)理統(tǒng)計方法對率定期和驗證期模擬和實測水文過程進行對比分析,結(jié)果見圖5。由圖5可知,汛期模擬流量的置信區(qū)間寬度大于實測流量的,但模擬流量平均值低于實測流量值平均值;非汛期模擬流量置信區(qū)間近似于實測流量值且異常值較少,同時平均流量值集中在25.68 m3/s;此外,非汛期無論是模擬值還是實測值其置信區(qū)間均小于汛期,意味著非汛期的徑流模擬精度高于汛期,這進一步驗證了小流量時模型能夠很好地捕捉到流量過程,而大流量時捕捉結(jié)果較差。

    4.3 數(shù)字濾波基流分割結(jié)果分析

    4.3.1 率定期和驗證期基流分割結(jié)果

    依據(jù)所構(gòu)建的北洛河流域TOPMODEL模型徑流模擬結(jié)果,為了對比分析徑流模擬結(jié)果,以頭水文站1986—1990年實測日徑流序列為依據(jù),采用數(shù)字濾波方法分別對徑流模擬序列和實測徑流序列進行基流分割,從而分割出地表徑流和地下徑流(基流)。此外,選取單一的基流指數(shù)(BFI)作為評價指標,對率定期和驗證期基流分割結(jié)果進行評估,結(jié)果見表3。由表3可知,率定期基流指數(shù)(BFI)均大于驗證期且均大于0.75,這表明北洛河流域地下徑流對河道徑流的補給作用顯著,這也進一步反映了流域降雨過程主要集中在汛期。此外,采用TOPMODEL模型模擬的流量在率定期和驗證期均大于實測流量,但總的來說相差微小。

    為了深入了解數(shù)字濾波方法所得基流分割過程的特點,采用數(shù)字濾波算法分別對模擬序列和實測序列的率定期與驗證期進行基流分割,結(jié)果見圖6。由圖6可知:①不同時期基流過程線起伏程度差異較大,但其形狀基本與徑流過程線一致;②模擬序列和實測序列基流分割過程線變化趨勢基本一致,但實測序列在汛期的分割過程線更符合流域?qū)嶋H徑流過程;③汛期受集中降水作用,率定期和驗證期地表徑流的增加會導(dǎo)致土壤下滲量增大,基流量也隨之增多,其趨勢與徑流過程相同。

    4.3.2 基流分割過程對徑流量的影響

    圖7為對不同時期采用數(shù)字濾波算法估算的地表徑流和地下徑流(基流)的影響量。由圖7可知:①率定期模擬徑流量大于實測徑流量,兩者之間的差值流量為1 964.60 m3/s,模擬地表徑流量小于實測地表徑流量,流量差為632.49 m3/s,模擬基流量大于實測基流量,流量差值為1 332.11 m3/s;②驗證期模擬徑流量小于實測徑流量,模擬地表徑流和基流量均小于實測地表徑流量和基流量,其流量差值分別為826.91、1 022.73、195.82 m3/s;這主要歸因于驗證期模型模擬精度較率定期差,進而導(dǎo)致驗證期模擬水文過程被低估,進而導(dǎo)致驗證期模擬的徑流成分(基流、地表徑流)較實測水文過程小??偟膩碚f,率定期基流分割效果優(yōu)于驗證期。

    5 結(jié) 論

    (1)通過Mamn-Kendall突變檢驗和有序聚類方法確定頭水文站徑流量的突變年份為1994年,在此基礎(chǔ)上選取TOPMODEL模型的率定期為1986—1988年,驗證期為1989—1990年。

    (2)構(gòu)建的北洛河流域TOPMODEL模型模擬精度較高,率定期和驗證期的NSE分別為0.69、0.63,R2分別為0.87、0.82, Re 均小于20%。

    (3)非汛期徑流模擬精度優(yōu)于汛期,這為進行基流分割奠定了良好基礎(chǔ)。

    (4)不同時期基流過程線起伏程度差異較大,但其形狀基本與徑流過程線一致;率定期模擬徑流量和基流量均大于實測徑流量,但模擬地表徑流量小于實測地表徑流量;驗證期模擬流量大于實測流量,模擬地表徑流和基流量均大于實測地表徑流量和基流量,其流量差值分別為826.91、1 022.73、195.82 m3/s。

    參考文獻:

    [1] 雷泳南,張曉萍,張建軍,等.自動基流分割法在黃土高原水蝕風(fēng)蝕交錯區(qū)典型流域適用性分析[J].中國水土保持科學(xué),2011,9(6):57-64.

    [2] 李瑞,張士鋒.兩種自動基流分割方法在干旱半干旱地區(qū)的對比研究[J].地理科學(xué)進展,2017,36(7):864-872.

    [3] 豆林,黃明斌.自動基流分割方法在黃土區(qū)流域的應(yīng)用研究[J].水土保持通報,2010,30(3):107-111.

    [4] 焦瑋,朱仲元,宋小園,等.基流分割方法在錫林河流域適用性分析[J].干旱區(qū)研究,2017,34(1):26-35.

    [5] 黃文,李占玲,徐宗學(xué),等.基流分割方法在寒旱區(qū)流域的適用性分析[J].水文,2018,38(3):21-28.

    [6] 高東東,吳勇,陳盟,等.貢嘎山森林系統(tǒng)小流域基流分割與降雨入滲補給計算[J].長江流域資源與環(huán)境,2015,24(6):949-955.

    [7] TULARAM G A, ILAHEE M. Exponential Smoothing Method of Basic flow Separation and Its Impact on Continuous Loss Estimates[J]. American Journal of Environmental Sciences, 2008, 4(2): 136-144.

    [8] ZHANG J, ZHANG Y, SONG J, et al. Evaluating Relative Merits of Four Basic Flow Separation Methods in Eastern Australia[J]. Journal of Hydrology, 2017, 549(1): 252-263.

    [9] 楊軍軍,高小紅,李其江,等.湟水流域SWAT模型構(gòu)建及參數(shù)不確定性分析[J].水土保持研究,2013,20(1):82-88.

    [10] 宮興龍,付強,王斌,等.聯(lián)合改進TOPMODEL和PDSI的半干旱區(qū)干旱評估模型構(gòu)建[J].農(nóng)業(yè)工程學(xué)報,2018,34(6):137-144.

    [11] 郭巧玲,楊云松,魯學(xué)綱.黑河流域1957—2008年徑流變化特性分析[J].水資源與水工程學(xué)報,2011,22(3):77-81.

    [12] 羅開盛,陶福祿.基于SWAT的西北干旱區(qū)縣域水文模擬:以臨澤縣為例[J].生態(tài)學(xué)報,2018,38(23):8593-8603.

    [13] 熊立華,郭生練.分布式流域水文模型[M].北京:中國水利水電出版社,2004:142-143.

    [14] 李抗彬,沈冰,宋孝玉,等.TOPMODEL模型在半濕潤地區(qū)徑流模擬分析中的應(yīng)用及改進[J].水利學(xué)報,2015,46(12):1453-1459.

    [15] ECKHARD T K. A Comparison of Basic flow Indices, Which were Calculated with Seven Different Basic flow Separation Methods [J]. Journal of Hydrology, 2008, 352(1): 168-173.

    [16] 李芳,靳少波,鄒松兵,等.黃河源區(qū)白河基流分割方法適用性分析[J].冰川凍土,2018,40(3):679-689.

    [17] 符傳君,黃國如,陳永勤. 用TOPMODEL模型模擬流域枯水徑流[J].應(yīng)用基礎(chǔ)與工程科學(xué)學(xué)報,2007(4):509-516.

    [18] 劉吉開,萬甜,程文,等.未來氣候情境下渭河流域陜西段非點源污染負荷響應(yīng)[J].水土保持通報,2018,38(4):82-86.

    [19] 黎云云,暢建霞,王於琪,等.渭河流域河川徑流對氣候變化的時空響應(yīng)機理[J].應(yīng)用基礎(chǔ)與工程科學(xué)學(xué)報,2018,26(3):502-514.

    [20] 周帥.SWAT模型參數(shù)不確定性分析及動態(tài)性評估[D].西安:西安理工大學(xué),2018:17-22.

    【責(zé)任編輯 翟戌亮】

    猜你喜歡
    數(shù)字濾波基流徑流
    高精度傳感器溫度漂移補償分析與研究
    美國長期氣象干旱對河道基流的影響研究
    數(shù)字濾波在語音信號降噪中的應(yīng)用
    電子制作(2019年24期)2019-02-23 13:22:40
    基于改進Tennant法的洋河流域生態(tài)基流估算研究
    Topmodel在布哈河流域徑流模擬中的應(yīng)用
    不同生態(tài)基流對水利樞紐發(fā)電損益的影響
    濾波器長度對濾波結(jié)果的影響研究
    探秘“大徑流”
    攻克“大徑流”
    增加河流生態(tài)基流的山間河谷型調(diào)蓄模式研究
    亚洲精品成人久久久久久| 婷婷色综合大香蕉| 日本与韩国留学比较| 欧美又色又爽又黄视频| 天天一区二区日本电影三级| 黄色一级大片看看| 中国美白少妇内射xxxbb| 搡老熟女国产l中国老女人| 亚洲av五月六月丁香网| 又黄又爽又刺激的免费视频.| 99热精品在线国产| 成人高潮视频无遮挡免费网站| 亚洲人成网站在线播放欧美日韩| 久久精品人妻少妇| 国产一区二区三区av在线 | 日韩,欧美,国产一区二区三区 | 长腿黑丝高跟| 卡戴珊不雅视频在线播放| 国内精品宾馆在线| 欧美成人免费av一区二区三区| 亚洲人成网站高清观看| 色5月婷婷丁香| 国产精品亚洲一级av第二区| 蜜臀久久99精品久久宅男| 五月伊人婷婷丁香| 听说在线观看完整版免费高清| 卡戴珊不雅视频在线播放| 久99久视频精品免费| 亚洲av不卡在线观看| 精品久久久久久久久久免费视频| 国产精品乱码一区二三区的特点| 亚洲欧美精品自产自拍| 男女下面进入的视频免费午夜| 国产高清视频在线观看网站| 亚洲婷婷狠狠爱综合网| 夜夜夜夜夜久久久久| 色尼玛亚洲综合影院| 欧美一级a爱片免费观看看| 久久久久久久久久成人| 亚洲中文字幕日韩| 99久久精品热视频| 免费大片18禁| 午夜激情欧美在线| 可以在线观看毛片的网站| 亚洲欧美日韩东京热| 成年女人看的毛片在线观看| 国产 一区精品| 波多野结衣高清作品| 在线观看美女被高潮喷水网站| 最近手机中文字幕大全| 晚上一个人看的免费电影| 又爽又黄a免费视频| 久久久久久久久中文| 人人妻人人看人人澡| 亚洲精品456在线播放app| 91av网一区二区| 别揉我奶头 嗯啊视频| 长腿黑丝高跟| 精品福利观看| 中文字幕av在线有码专区| 黄色欧美视频在线观看| 中出人妻视频一区二区| 国产中年淑女户外野战色| 日韩欧美精品免费久久| 狠狠狠狠99中文字幕| 国产精品久久久久久精品电影| 国产私拍福利视频在线观看| av天堂在线播放| 色哟哟哟哟哟哟| 黄色配什么色好看| 精品少妇黑人巨大在线播放 | 欧美3d第一页| 国产蜜桃级精品一区二区三区| 午夜影院日韩av| 两性午夜刺激爽爽歪歪视频在线观看| 国产美女午夜福利| 在线观看一区二区三区| 又黄又爽又免费观看的视频| 51国产日韩欧美| 亚洲精品粉嫩美女一区| 亚洲国产精品国产精品| 天堂动漫精品| 久久久久久久亚洲中文字幕| 毛片女人毛片| 久久这里只有精品中国| 欧美激情久久久久久爽电影| 中出人妻视频一区二区| 成人午夜高清在线视频| 亚洲美女黄片视频| 免费人成在线观看视频色| 亚洲婷婷狠狠爱综合网| eeuss影院久久| 久久亚洲精品不卡| 日韩欧美 国产精品| 国产高清不卡午夜福利| 欧美不卡视频在线免费观看| 久久人妻av系列| 精品日产1卡2卡| 99热这里只有是精品在线观看| 18禁在线播放成人免费| 亚洲av美国av| 成人特级黄色片久久久久久久| 亚洲中文日韩欧美视频| 久久久久国产精品人妻aⅴ院| 一区二区三区高清视频在线| 丰满乱子伦码专区| 1024手机看黄色片| 成人美女网站在线观看视频| 人人妻人人澡人人爽人人夜夜 | 免费一级毛片在线播放高清视频| av.在线天堂| 91av网一区二区| 久久久久久久午夜电影| 国产精品一区www在线观看| 日本三级黄在线观看| 国产精品电影一区二区三区| 国产色爽女视频免费观看| 亚洲五月天丁香| 国产国拍精品亚洲av在线观看| 精品午夜福利视频在线观看一区| 两个人视频免费观看高清| 国产蜜桃级精品一区二区三区| 久久精品久久久久久噜噜老黄 | 熟妇人妻久久中文字幕3abv| 国产大屁股一区二区在线视频| 99久国产av精品| 老师上课跳d突然被开到最大视频| 又黄又爽又刺激的免费视频.| 晚上一个人看的免费电影| 久久久久性生活片| 国产久久久一区二区三区| 夜夜夜夜夜久久久久| 国产精品永久免费网站| 在线看三级毛片| 午夜激情欧美在线| 简卡轻食公司| 国产片特级美女逼逼视频| 日本撒尿小便嘘嘘汇集6| 如何舔出高潮| 别揉我奶头 嗯啊视频| 又黄又爽又刺激的免费视频.| 免费一级毛片在线播放高清视频| 人妻制服诱惑在线中文字幕| 女同久久另类99精品国产91| 国产高清不卡午夜福利| 99视频精品全部免费 在线| 蜜臀久久99精品久久宅男| 一a级毛片在线观看| 久久久久久久久久黄片| 欧美成人一区二区免费高清观看| aaaaa片日本免费| 极品教师在线视频| 内射极品少妇av片p| 国产精品无大码| 亚洲内射少妇av| 亚洲精品国产av成人精品 | 日韩大尺度精品在线看网址| 啦啦啦韩国在线观看视频| 久久久久久久久中文| 成人特级黄色片久久久久久久| 午夜福利18| 亚州av有码| 午夜福利在线观看吧| 如何舔出高潮| 伦精品一区二区三区| 日日啪夜夜撸| 久久精品综合一区二区三区| 在线免费十八禁| 午夜精品一区二区三区免费看| 国产老妇女一区| 99久久精品热视频| 丝袜美腿在线中文| 免费看av在线观看网站| 日日摸夜夜添夜夜添av毛片| 日韩精品青青久久久久久| 91久久精品国产一区二区三区| 亚洲国产色片| 成人一区二区视频在线观看| 亚洲精品国产av成人精品 | 精品久久久久久成人av| 亚洲专区国产一区二区| 日产精品乱码卡一卡2卡三| 丰满乱子伦码专区| 国产亚洲av嫩草精品影院| 久久久久久久亚洲中文字幕| 大又大粗又爽又黄少妇毛片口| 亚洲欧美精品自产自拍| 一级a爱片免费观看的视频| 亚洲最大成人av| 麻豆精品久久久久久蜜桃| 午夜免费男女啪啪视频观看 | 天天躁日日操中文字幕| 午夜激情欧美在线| 国产伦精品一区二区三区四那| 日本黄大片高清| 国产精品,欧美在线| av卡一久久| 国内久久婷婷六月综合欲色啪| 国产成人freesex在线 | 中文字幕精品亚洲无线码一区| 国产欧美日韩精品一区二区| 国产成人精品久久久久久| 久久久成人免费电影| 亚洲七黄色美女视频| 极品教师在线视频| 天堂av国产一区二区熟女人妻| 亚洲乱码一区二区免费版| 俺也久久电影网| 久久久久免费精品人妻一区二区| 午夜a级毛片| 午夜爱爱视频在线播放| 熟女电影av网| 卡戴珊不雅视频在线播放| 亚洲在线自拍视频| 国产蜜桃级精品一区二区三区| 乱系列少妇在线播放| 亚洲最大成人中文| 97在线视频观看| 直男gayav资源| 亚洲精品日韩在线中文字幕 | 欧美国产日韩亚洲一区| 国产午夜精品久久久久久一区二区三区 | 日产精品乱码卡一卡2卡三| 日日摸夜夜添夜夜添av毛片| 亚洲一区二区三区色噜噜| 最近视频中文字幕2019在线8| 久久久久九九精品影院| 真实男女啪啪啪动态图| 九九在线视频观看精品| 国产亚洲91精品色在线| 99热6这里只有精品| 一个人看视频在线观看www免费| 国产三级在线视频| 成人午夜高清在线视频| 亚洲天堂国产精品一区在线| 日韩人妻高清精品专区| 校园人妻丝袜中文字幕| 中文在线观看免费www的网站| 日韩精品青青久久久久久| 欧美在线一区亚洲| 免费观看精品视频网站| 全区人妻精品视频| 色哟哟·www| 午夜激情福利司机影院| 日韩在线高清观看一区二区三区| 久久精品人妻少妇| 一区二区三区四区激情视频 | 亚洲中文日韩欧美视频| 精品一区二区免费观看| 国产中年淑女户外野战色| 国产亚洲精品久久久久久毛片| 亚洲av一区综合| 黄色欧美视频在线观看| 亚洲成人精品中文字幕电影| 久久人人爽人人片av| 在线观看66精品国产| 嫩草影院精品99| 国产综合懂色| 五月伊人婷婷丁香| 最近2019中文字幕mv第一页| 欧美区成人在线视频| 国产亚洲91精品色在线| 午夜福利成人在线免费观看| 久久热精品热| 成人av在线播放网站| 小说图片视频综合网站| 亚洲精品影视一区二区三区av| 少妇的逼水好多| 日韩欧美三级三区| 亚洲天堂国产精品一区在线| 九九久久精品国产亚洲av麻豆| 国产69精品久久久久777片| 精品一区二区三区人妻视频| 日本一二三区视频观看| av福利片在线观看| 日日干狠狠操夜夜爽| 看黄色毛片网站| 精品乱码久久久久久99久播| 亚洲国产精品久久男人天堂| avwww免费| 国产69精品久久久久777片| ponron亚洲| 午夜老司机福利剧场| 久久久久国产网址| 中国国产av一级| 99热精品在线国产| 国产av不卡久久| 天堂影院成人在线观看| 国产高清激情床上av| 欧美色欧美亚洲另类二区| 日韩欧美 国产精品| 午夜老司机福利剧场| 久久久久国产精品人妻aⅴ院| 中国国产av一级| 亚洲丝袜综合中文字幕| 国产精品人妻久久久影院| 亚洲不卡免费看| 国产三级在线视频| 日本黄色片子视频| 人人妻,人人澡人人爽秒播| 好男人在线观看高清免费视频| 波野结衣二区三区在线| 日韩av在线大香蕉| 亚洲欧美日韩东京热| 中国国产av一级| 国产精品野战在线观看| 男女边吃奶边做爰视频| 国产精品野战在线观看| 97在线视频观看| 成人综合一区亚洲| 国产精品人妻久久久影院| 国产av一区在线观看免费| 深夜a级毛片| 在线看三级毛片| 国产乱人视频| 中文字幕精品亚洲无线码一区| 国产女主播在线喷水免费视频网站 | 久久久久久久久中文| 天堂√8在线中文| 成人三级黄色视频| 午夜老司机福利剧场| 99久久九九国产精品国产免费| 日韩欧美免费精品| 一个人看视频在线观看www免费| 国产一区二区三区在线臀色熟女| 久久久欧美国产精品| 亚洲熟妇熟女久久| 国产欧美日韩精品一区二区| 男女边吃奶边做爰视频| 午夜精品一区二区三区免费看| 99热网站在线观看| 国产乱人偷精品视频| 精品国产三级普通话版| 精品99又大又爽又粗少妇毛片| 日本五十路高清| 美女 人体艺术 gogo| 色哟哟哟哟哟哟| 国产亚洲av嫩草精品影院| 欧美最新免费一区二区三区| 性插视频无遮挡在线免费观看| 日韩国内少妇激情av| 成年女人永久免费观看视频| 中文字幕av在线有码专区| 国产一区二区亚洲精品在线观看| 波多野结衣巨乳人妻| 97超视频在线观看视频| 女的被弄到高潮叫床怎么办| 91av网一区二区| 久久精品夜色国产| 97碰自拍视频| 韩国av在线不卡| 三级国产精品欧美在线观看| 亚洲国产精品成人久久小说 | 久久久成人免费电影| 久久精品国产亚洲av天美| 国产一区二区在线观看日韩| 午夜福利成人在线免费观看| 久久精品国产亚洲av香蕉五月| 精品99又大又爽又粗少妇毛片| 婷婷色综合大香蕉| 亚洲精品乱码久久久v下载方式| 五月伊人婷婷丁香| 蜜桃久久精品国产亚洲av| 成人国产麻豆网| 欧美高清性xxxxhd video| 精品久久久久久久久亚洲| 性欧美人与动物交配| 日本与韩国留学比较| 日本黄色片子视频| 国产一区二区在线观看日韩| 日本三级黄在线观看| 九九久久精品国产亚洲av麻豆| 久久精品国产亚洲av天美| 国产伦精品一区二区三区视频9| 色尼玛亚洲综合影院| 丰满的人妻完整版| 久久久精品大字幕| 床上黄色一级片| av专区在线播放| 丰满乱子伦码专区| 亚洲精华国产精华液的使用体验 | aaaaa片日本免费| 人人妻人人澡欧美一区二区| 欧美成人精品欧美一级黄| 亚洲欧美日韩无卡精品| 色av中文字幕| 草草在线视频免费看| 亚洲第一电影网av| 亚洲av免费在线观看| 亚洲在线观看片| 国产一区亚洲一区在线观看| 免费一级毛片在线播放高清视频| 黄色视频,在线免费观看| 中国美白少妇内射xxxbb| 国产精品女同一区二区软件| 国内精品一区二区在线观看| 国模一区二区三区四区视频| 免费观看在线日韩| 麻豆成人午夜福利视频| 日韩精品青青久久久久久| 国产视频内射| 少妇丰满av| 欧美高清性xxxxhd video| 丰满的人妻完整版| 免费观看精品视频网站| 激情 狠狠 欧美| 国产成人一区二区在线| 成人午夜高清在线视频| 九九在线视频观看精品| 国产真实乱freesex| 久久欧美精品欧美久久欧美| 久久草成人影院| 日韩精品有码人妻一区| 老师上课跳d突然被开到最大视频| 国产成年人精品一区二区| 国产精品电影一区二区三区| 欧美性感艳星| 又黄又爽又免费观看的视频| av在线播放精品| 91在线观看av| 亚洲国产欧美人成| 日本一本二区三区精品| 亚洲人成网站高清观看| 91午夜精品亚洲一区二区三区| 最近手机中文字幕大全| 久久久久国产精品人妻aⅴ院| 久久久久久大精品| 国产又黄又爽又无遮挡在线| av天堂在线播放| 国产成人aa在线观看| 永久网站在线| 亚洲国产精品成人综合色| 久久鲁丝午夜福利片| 久久综合国产亚洲精品| 熟妇人妻久久中文字幕3abv| 亚洲最大成人av| 久久久久久久久久久丰满| 精华霜和精华液先用哪个| 一级av片app| 亚洲自拍偷在线| 国产视频内射| av中文乱码字幕在线| 国产人妻一区二区三区在| 寂寞人妻少妇视频99o| 性插视频无遮挡在线免费观看| 久久鲁丝午夜福利片| 一级毛片电影观看 | 久久鲁丝午夜福利片| 国产精品爽爽va在线观看网站| 人妻丰满熟妇av一区二区三区| 国内揄拍国产精品人妻在线| 免费看av在线观看网站| 日韩精品青青久久久久久| 欧美色视频一区免费| 你懂的网址亚洲精品在线观看 | 蜜臀久久99精品久久宅男| 成人av在线播放网站| 不卡一级毛片| 春色校园在线视频观看| 免费av不卡在线播放| 精品久久久久久久久av| 人妻丰满熟妇av一区二区三区| 99热只有精品国产| 国产在视频线在精品| 国产v大片淫在线免费观看| 色噜噜av男人的天堂激情| 丝袜喷水一区| 97在线视频观看| 两性午夜刺激爽爽歪歪视频在线观看| 99riav亚洲国产免费| 精品人妻偷拍中文字幕| 欧美成人精品欧美一级黄| 免费电影在线观看免费观看| 免费一级毛片在线播放高清视频| 狠狠狠狠99中文字幕| 在线a可以看的网站| 亚洲欧美中文字幕日韩二区| 精品人妻视频免费看| 久久精品国产鲁丝片午夜精品| АⅤ资源中文在线天堂| 国产亚洲精品av在线| 一个人观看的视频www高清免费观看| 国产熟女欧美一区二区| 久久精品国产亚洲网站| 黄片wwwwww| 国产精品伦人一区二区| 国产精品一及| 午夜视频国产福利| 欧美一区二区亚洲| 插逼视频在线观看| 在线看三级毛片| 男女做爰动态图高潮gif福利片| 亚洲人成网站高清观看| 偷拍熟女少妇极品色| 日本成人三级电影网站| 欧美xxxx性猛交bbbb| 国产不卡一卡二| 亚洲在线观看片| 最近视频中文字幕2019在线8| 一个人看的www免费观看视频| 亚洲国产色片| 3wmmmm亚洲av在线观看| 啦啦啦观看免费观看视频高清| 18+在线观看网站| 精品免费久久久久久久清纯| 国产极品精品免费视频能看的| 免费在线观看成人毛片| 久久久久久伊人网av| 久久久久久久久中文| 天天躁日日操中文字幕| 中文字幕av成人在线电影| 一卡2卡三卡四卡精品乱码亚洲| 国产精品一及| 色5月婷婷丁香| 最后的刺客免费高清国语| 久久精品人妻少妇| 别揉我奶头 嗯啊视频| 国国产精品蜜臀av免费| 国产毛片a区久久久久| 99久国产av精品国产电影| 日日啪夜夜撸| 色尼玛亚洲综合影院| 精品久久久久久久久av| 伦理电影大哥的女人| 国产人妻一区二区三区在| 黄片wwwwww| 欧美最黄视频在线播放免费| av视频在线观看入口| 亚洲第一区二区三区不卡| 欧美3d第一页| 国产精品人妻久久久影院| 欧美激情国产日韩精品一区| 日本一二三区视频观看| 精品人妻视频免费看| 精品久久久久久久久久免费视频| 一个人免费在线观看电影| 永久网站在线| 欧美成人精品欧美一级黄| 久久精品影院6| 欧美激情久久久久久爽电影| 国产人妻一区二区三区在| 在线观看一区二区三区| 午夜视频国产福利| 一区二区三区免费毛片| 亚洲av中文av极速乱| 国产乱人偷精品视频| 中国国产av一级| 高清午夜精品一区二区三区 | 久久久久免费精品人妻一区二区| 男插女下体视频免费在线播放| 午夜福利18| 黄色欧美视频在线观看| 日日摸夜夜添夜夜爱| 12—13女人毛片做爰片一| 精品99又大又爽又粗少妇毛片| 亚洲国产欧美人成| 三级男女做爰猛烈吃奶摸视频| 麻豆精品久久久久久蜜桃| 免费不卡的大黄色大毛片视频在线观看 | 亚洲精品色激情综合| 日韩国内少妇激情av| av女优亚洲男人天堂| 日韩中字成人| 国产在线男女| 久久天躁狠狠躁夜夜2o2o| 免费看a级黄色片| 久久久久久久亚洲中文字幕| 成人午夜高清在线视频| 免费观看在线日韩| 热99re8久久精品国产| 少妇人妻一区二区三区视频| 少妇裸体淫交视频免费看高清| 国产在视频线在精品| 国产精品免费一区二区三区在线| 人妻丰满熟妇av一区二区三区| 97在线视频观看| 亚洲性久久影院| 精品免费久久久久久久清纯| 天天一区二区日本电影三级| 两个人视频免费观看高清| 在线国产一区二区在线| 国产精品三级大全| 日韩一区二区视频免费看| 天堂√8在线中文| 成人精品一区二区免费| 99九九线精品视频在线观看视频| 国产精品乱码一区二三区的特点| 搡老妇女老女人老熟妇| 亚洲人成网站高清观看| 国产精品爽爽va在线观看网站| 搞女人的毛片| videossex国产| 国产精品亚洲美女久久久| 欧美性猛交╳xxx乱大交人| 亚洲国产日韩欧美精品在线观看| 18禁在线播放成人免费| 久久久久国产精品人妻aⅴ院| 69av精品久久久久久| 你懂的网址亚洲精品在线观看 | 国产精品野战在线观看| 女生性感内裤真人,穿戴方法视频| 小蜜桃在线观看免费完整版高清| 超碰av人人做人人爽久久| 久久久久国产网址| 国产精品久久视频播放| 国产熟女欧美一区二区| 亚洲国产精品成人久久小说 | 久久亚洲精品不卡| 国产色婷婷99| 嫩草影院精品99| 国产午夜精品论理片| 日本a在线网址| 成人永久免费在线观看视频| 一夜夜www| 中文亚洲av片在线观看爽|