• <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)蓄模式研究
    久久亚洲国产成人精品v| 亚洲第一av免费看| 久久久久久久久大av| 日本欧美视频一区| 国产伦精品一区二区三区视频9| 最近最新中文字幕大全电影3| 夜夜骑夜夜射夜夜干| 精品久久国产蜜桃| 久久精品国产自在天天线| 直男gayav资源| 中国三级夫妇交换| 免费大片黄手机在线观看| 免费看日本二区| 国产精品一及| 人人妻人人爽人人添夜夜欢视频 | 国产亚洲av片在线观看秒播厂| 精品一区二区三区视频在线| 舔av片在线| 亚洲第一av免费看| 丰满迷人的少妇在线观看| 国产av码专区亚洲av| 成年免费大片在线观看| 国产精品av视频在线免费观看| 成人亚洲精品一区在线观看 | videossex国产| 男人舔奶头视频| 免费av不卡在线播放| 久久av网站| 精品久久久精品久久久| 国产男女超爽视频在线观看| av在线播放精品| 99视频精品全部免费 在线| 99久久精品国产国产毛片| 亚洲av欧美aⅴ国产| 国产精品蜜桃在线观看| 亚洲精品乱久久久久久| 永久免费av网站大全| 日韩,欧美,国产一区二区三区| 女的被弄到高潮叫床怎么办| 久久久色成人| 亚洲国产色片| 亚洲精品久久午夜乱码| 国产乱人视频| 国产毛片在线视频| 久久久久精品性色| 久久精品久久久久久噜噜老黄| 美女中出高潮动态图| 狠狠精品人妻久久久久久综合| 下体分泌物呈黄色| 我要看黄色一级片免费的| 高清av免费在线| 国产在线男女| 国产午夜精品久久久久久一区二区三区| 特大巨黑吊av在线直播| 久久99热这里只频精品6学生| 精品少妇黑人巨大在线播放| 国产中年淑女户外野战色| 91在线精品国自产拍蜜月| xxx大片免费视频| 亚洲综合色惰| 国产av一区二区精品久久 | 日日撸夜夜添| 日本免费在线观看一区| 麻豆精品久久久久久蜜桃| a级毛片免费高清观看在线播放| 国产爱豆传媒在线观看| 久久久久久久久久人人人人人人| 久久毛片免费看一区二区三区| 18禁在线无遮挡免费观看视频| 五月伊人婷婷丁香| 大陆偷拍与自拍| 99热这里只有精品一区| 欧美日韩视频高清一区二区三区二| 国产熟女欧美一区二区| av在线老鸭窝| 久久精品国产a三级三级三级| 国产精品久久久久久精品古装| 丰满迷人的少妇在线观看| 国产精品一区二区性色av| 亚洲图色成人| 卡戴珊不雅视频在线播放| 网址你懂的国产日韩在线| h视频一区二区三区| av线在线观看网站| 日韩视频在线欧美| h视频一区二区三区| 中文字幕久久专区| 人妻一区二区av| 另类亚洲欧美激情| 久久久久久久久久人人人人人人| 插逼视频在线观看| 插逼视频在线观看| 国产男女超爽视频在线观看| 永久免费av网站大全| av在线观看视频网站免费| 久久久久久九九精品二区国产| 天天躁夜夜躁狠狠久久av| 国产av码专区亚洲av| 夫妻性生交免费视频一级片| av网站免费在线观看视频| 亚洲av中文av极速乱| 久久精品国产自在天天线| 免费久久久久久久精品成人欧美视频 | 亚洲av免费高清在线观看| 欧美xxxx黑人xx丫x性爽| 六月丁香七月| 久久久久国产精品人妻一区二区| 在线免费十八禁| 插阴视频在线观看视频| 蜜桃久久精品国产亚洲av| 夜夜爽夜夜爽视频| 亚洲av在线观看美女高潮| 亚洲丝袜综合中文字幕| 国产黄片美女视频| 欧美亚洲 丝袜 人妻 在线| 亚洲精品色激情综合| 国产av国产精品国产| 水蜜桃什么品种好| 久久精品熟女亚洲av麻豆精品| 久久国产精品大桥未久av | 亚洲欧美一区二区三区黑人 | 日日摸夜夜添夜夜添av毛片| 美女脱内裤让男人舔精品视频| 一区二区三区精品91| 爱豆传媒免费全集在线观看| 91精品国产国语对白视频| 亚洲国产精品成人久久小说| 最近2019中文字幕mv第一页| 少妇人妻精品综合一区二区| av一本久久久久| 夜夜骑夜夜射夜夜干| 精品久久久噜噜| 亚洲欧美一区二区三区黑人 | 日本黄色日本黄色录像| 日韩人妻高清精品专区| 精品一区二区免费观看| 精品久久国产蜜桃| 在现免费观看毛片| 亚洲国产精品国产精品| 国产亚洲91精品色在线| 一本色道久久久久久精品综合| 丰满乱子伦码专区| av国产久精品久网站免费入址| 如何舔出高潮| 黄色一级大片看看| 少妇丰满av| 午夜福利视频精品| 欧美极品一区二区三区四区| 99re6热这里在线精品视频| 国产精品国产三级专区第一集| 一本—道久久a久久精品蜜桃钙片| 能在线免费看毛片的网站| 久久精品国产a三级三级三级| 亚洲经典国产精华液单| 日本午夜av视频| 五月玫瑰六月丁香| 色婷婷久久久亚洲欧美| 黄色欧美视频在线观看| 国产成人精品久久久久久| 在线观看三级黄色| 国产亚洲精品久久久com| av在线老鸭窝| 99视频精品全部免费 在线| 2022亚洲国产成人精品| 精品一区在线观看国产| 亚洲精品成人av观看孕妇| 国产精品人妻久久久久久| 自拍偷自拍亚洲精品老妇| 少妇裸体淫交视频免费看高清| 国产视频内射| 日韩国内少妇激情av| 国产 一区精品| a级毛片免费高清观看在线播放| 人妻少妇偷人精品九色| 午夜免费观看性视频| 日韩欧美 国产精品| 国产亚洲最大av| 亚洲欧美日韩另类电影网站 | 青春草视频在线免费观看| 国产一区有黄有色的免费视频| 美女脱内裤让男人舔精品视频| 爱豆传媒免费全集在线观看| 一本一本综合久久| 日韩电影二区| 永久免费av网站大全| 精品久久国产蜜桃| 亚洲欧美一区二区三区国产| 久久毛片免费看一区二区三区| 久久久午夜欧美精品| 在线观看美女被高潮喷水网站| 蜜桃亚洲精品一区二区三区| 日韩国内少妇激情av| 一区二区三区精品91| 九色成人免费人妻av| 久久99热这里只频精品6学生| 成年人午夜在线观看视频| 观看免费一级毛片| 欧美日韩综合久久久久久| 亚洲美女视频黄频| 亚洲高清免费不卡视频| 亚洲精品一区蜜桃| 国产欧美亚洲国产| 亚洲成人av在线免费| 久久精品久久久久久噜噜老黄| 久久精品国产自在天天线| 亚洲色图综合在线观看| 乱系列少妇在线播放| 人妻 亚洲 视频| 国产免费一级a男人的天堂| 精品午夜福利在线看| av不卡在线播放| 久久精品熟女亚洲av麻豆精品| 欧美精品人与动牲交sv欧美| 日韩强制内射视频| 人妻系列 视频| 国产在线男女| 成人毛片a级毛片在线播放| 日本午夜av视频| 校园人妻丝袜中文字幕| 最近2019中文字幕mv第一页| 麻豆精品久久久久久蜜桃| av一本久久久久| 国产精品一区www在线观看| 18+在线观看网站| 91午夜精品亚洲一区二区三区| 精品久久久久久电影网| 超碰av人人做人人爽久久| 水蜜桃什么品种好| 亚洲国产精品999| 成人黄色视频免费在线看| 91精品国产国语对白视频| 欧美成人精品欧美一级黄| 国产探花极品一区二区| 久久国产乱子免费精品| 国产高清有码在线观看视频| 三级国产精品片| 精品一品国产午夜福利视频| 交换朋友夫妻互换小说| 国产成人91sexporn| 欧美成人一区二区免费高清观看| 如何舔出高潮| 国产av精品麻豆| 色吧在线观看| 99久久精品国产国产毛片| 成人综合一区亚洲| 成人毛片60女人毛片免费| 亚洲美女黄色视频免费看| 久久久久久久精品精品| 91精品国产九色| 岛国毛片在线播放| 这个男人来自地球电影免费观看 | 人人妻人人添人人爽欧美一区卜 | 久久久久久人妻| 在线精品无人区一区二区三 | 麻豆乱淫一区二区| 久热这里只有精品99| 国产乱人视频| 欧美性感艳星| 蜜桃亚洲精品一区二区三区| 伦精品一区二区三区| 午夜福利网站1000一区二区三区| 日产精品乱码卡一卡2卡三| 婷婷色麻豆天堂久久| 免费大片18禁| 国产精品国产三级国产专区5o| 80岁老熟妇乱子伦牲交| 只有这里有精品99| 好男人视频免费观看在线| 精品一区二区三卡| 91久久精品电影网| 亚洲自偷自拍三级| 亚洲婷婷狠狠爱综合网| tube8黄色片| h日本视频在线播放| 午夜福利高清视频| 成人漫画全彩无遮挡| 亚洲人成网站在线播| 国产男女超爽视频在线观看| 一级片'在线观看视频| 狂野欧美白嫩少妇大欣赏| 国产爽快片一区二区三区| 成人高潮视频无遮挡免费网站| 国产有黄有色有爽视频| 七月丁香在线播放| 综合色丁香网| 国产美女午夜福利| av在线app专区| 成人午夜精彩视频在线观看| 网址你懂的国产日韩在线| 中国三级夫妇交换| 国产人妻一区二区三区在| 久久久a久久爽久久v久久| 日本-黄色视频高清免费观看| 黄片wwwwww| 亚洲aⅴ乱码一区二区在线播放| 精品亚洲乱码少妇综合久久| 亚洲av二区三区四区| 七月丁香在线播放| 一级毛片我不卡| 又黄又爽又刺激的免费视频.| 国产精品一区二区性色av| 久久久国产一区二区| 久久久久国产精品人妻一区二区| 国产精品人妻久久久久久| 欧美亚洲 丝袜 人妻 在线| 欧美成人a在线观看| 午夜免费男女啪啪视频观看| 久久国产乱子免费精品| 久久午夜福利片| 国产淫语在线视频| 久久99蜜桃精品久久| 王馨瑶露胸无遮挡在线观看| 国产色婷婷99| 高清黄色对白视频在线免费看 | 少妇高潮的动态图| 亚洲国产精品专区欧美| 日韩 亚洲 欧美在线| 免费看光身美女| 一本—道久久a久久精品蜜桃钙片| 99热这里只有是精品在线观看| 五月开心婷婷网| 99久久中文字幕三级久久日本| 国产片特级美女逼逼视频| 久久久久久久久久成人| 亚洲欧美成人综合另类久久久| 国产黄片视频在线免费观看| 欧美日韩国产mv在线观看视频 | 成人国产麻豆网| 六月丁香七月| 赤兔流量卡办理| 亚洲美女黄色视频免费看| 久久6这里有精品| tube8黄色片| 日韩一本色道免费dvd| 少妇的逼好多水| 日韩 亚洲 欧美在线| 中文资源天堂在线| 少妇人妻久久综合中文| 97精品久久久久久久久久精品| 九九在线视频观看精品| 国产一区亚洲一区在线观看| 在线 av 中文字幕| 青春草视频在线免费观看| 亚洲精品aⅴ在线观看| 爱豆传媒免费全集在线观看| 国内精品宾馆在线| 国产精品国产av在线观看| 免费播放大片免费观看视频在线观看| 熟女人妻精品中文字幕| 日韩国内少妇激情av| xxx大片免费视频| 美女高潮的动态| av在线蜜桃| 观看美女的网站| 又黄又爽又刺激的免费视频.| 2022亚洲国产成人精品| 一区二区三区乱码不卡18| av视频免费观看在线观看| 少妇丰满av| 精品一区二区三区视频在线| 91狼人影院| 欧美精品国产亚洲| 纯流量卡能插随身wifi吗| 大话2 男鬼变身卡| 91精品伊人久久大香线蕉| 久久这里有精品视频免费| 日日撸夜夜添| 亚洲av综合色区一区| 日本黄大片高清| 综合色丁香网| 一级爰片在线观看| 国产精品福利在线免费观看| 久久 成人 亚洲| 国产在线男女| 日本wwww免费看| 欧美精品亚洲一区二区| 三级经典国产精品| 国产精品免费大片| 免费看日本二区| 精品国产一区二区三区久久久樱花 | 欧美 日韩 精品 国产| 国产免费福利视频在线观看| 男男h啪啪无遮挡| 乱码一卡2卡4卡精品| 国产精品偷伦视频观看了| 国产av国产精品国产| 日韩欧美一区视频在线观看 | 街头女战士在线观看网站| 啦啦啦视频在线资源免费观看| 一级毛片aaaaaa免费看小| 一区在线观看完整版| 国产精品久久久久成人av| 人妻一区二区av| 在线观看三级黄色| 久热久热在线精品观看| 看十八女毛片水多多多| 啦啦啦在线观看免费高清www| 免费看av在线观看网站| 国产亚洲精品久久久com| 18+在线观看网站| 一级毛片 在线播放| 超碰97精品在线观看| 黄色日韩在线| 久久人人爽人人爽人人片va| 亚洲国产成人一精品久久久| 欧美精品人与动牲交sv欧美| 美女内射精品一级片tv| 欧美一级a爱片免费观看看| 最黄视频免费看| 天堂俺去俺来也www色官网| 你懂的网址亚洲精品在线观看| 精品少妇久久久久久888优播| 亚洲国产毛片av蜜桃av| 男人舔奶头视频| 亚洲国产最新在线播放| 免费看不卡的av| 中文欧美无线码| 亚洲国产最新在线播放| 麻豆成人av视频| 又粗又硬又长又爽又黄的视频| 国产v大片淫在线免费观看| 22中文网久久字幕| 亚洲精品国产色婷婷电影| 亚洲色图综合在线观看| 久久婷婷青草| 日韩三级伦理在线观看| 亚洲精品国产av蜜桃| 欧美精品一区二区免费开放| 国产 精品1| 国产男女超爽视频在线观看| 丰满迷人的少妇在线观看| 国产色婷婷99| 国产成人91sexporn| 中文字幕免费在线视频6| 久久99蜜桃精品久久| av不卡在线播放| 亚洲成色77777| 激情 狠狠 欧美| 夜夜爽夜夜爽视频| 十分钟在线观看高清视频www | 少妇人妻久久综合中文| 美女视频免费永久观看网站| 久久99精品国语久久久| 狂野欧美白嫩少妇大欣赏| 人妻夜夜爽99麻豆av| 嫩草影院入口| 人人妻人人看人人澡| 亚洲自偷自拍三级| 久久久久久久大尺度免费视频| 久久99精品国语久久久| 国产淫片久久久久久久久| 黄色日韩在线| 人人妻人人添人人爽欧美一区卜 | 男的添女的下面高潮视频| 天天躁日日操中文字幕| 99re6热这里在线精品视频| 多毛熟女@视频| 又黄又爽又刺激的免费视频.| av国产免费在线观看| 热99国产精品久久久久久7| 黄色日韩在线| 天美传媒精品一区二区| 欧美变态另类bdsm刘玥| 精品午夜福利在线看| videos熟女内射| 在线观看av片永久免费下载| 插阴视频在线观看视频| 日本黄色日本黄色录像| 国产黄片视频在线免费观看| 亚洲性久久影院| 青春草亚洲视频在线观看| av在线app专区| 欧美精品一区二区大全| 欧美zozozo另类| 少妇精品久久久久久久| 最新中文字幕久久久久| 超碰97精品在线观看| 午夜福利在线在线| 中文乱码字字幕精品一区二区三区| 少妇高潮的动态图| 欧美少妇被猛烈插入视频| 内地一区二区视频在线| av免费观看日本| 青春草国产在线视频| 久久久亚洲精品成人影院| 成年人午夜在线观看视频| 精品熟女少妇av免费看| 午夜老司机福利剧场| 成人高潮视频无遮挡免费网站| 久久久国产一区二区| 亚洲国产精品国产精品| 在线观看免费日韩欧美大片 | 国产大屁股一区二区在线视频| 亚洲欧美精品自产自拍| 成人18禁高潮啪啪吃奶动态图 | 99热这里只有精品一区| 亚洲色图av天堂| 免费久久久久久久精品成人欧美视频 | 精品少妇黑人巨大在线播放| 26uuu在线亚洲综合色| 日本av手机在线免费观看| 精品99又大又爽又粗少妇毛片| av.在线天堂| 嫩草影院新地址| tube8黄色片| 亚洲精品视频女| av国产精品久久久久影院| 亚洲av成人精品一区久久| 日本av免费视频播放| 午夜福利影视在线免费观看| 国产精品不卡视频一区二区| 日韩精品有码人妻一区| 丝袜脚勾引网站| 波野结衣二区三区在线| 日韩成人av中文字幕在线观看| 精品一区二区三卡| 亚洲av.av天堂| 亚洲美女搞黄在线观看| 免费高清在线观看视频在线观看| 午夜老司机福利剧场| 高清视频免费观看一区二区| 一边亲一边摸免费视频| 国产日韩欧美亚洲二区| 五月伊人婷婷丁香| 亚洲美女黄色视频免费看| 欧美xxxx黑人xx丫x性爽| 黄色一级大片看看| 日韩成人av中文字幕在线观看| a级一级毛片免费在线观看| 亚洲一级一片aⅴ在线观看| 亚洲成人一二三区av| 我的老师免费观看完整版| 午夜老司机福利剧场| 欧美区成人在线视频| 日日摸夜夜添夜夜爱| 亚洲欧美成人精品一区二区| 精品久久久久久久久av| h视频一区二区三区| 夫妻午夜视频| 大又大粗又爽又黄少妇毛片口| av又黄又爽大尺度在线免费看| 成人综合一区亚洲| 少妇的逼好多水| 成人美女网站在线观看视频| 日韩中文字幕视频在线看片 | 在线观看av片永久免费下载| 一本久久精品| 亚洲欧美精品自产自拍| 毛片一级片免费看久久久久| 国产精品一区二区在线观看99| 日韩免费高清中文字幕av| av福利片在线观看| 亚洲va在线va天堂va国产| 五月开心婷婷网| 美女cb高潮喷水在线观看| 天天躁夜夜躁狠狠久久av| 一级av片app| 亚洲成人中文字幕在线播放| 国产v大片淫在线免费观看| 欧美精品一区二区免费开放| 国产无遮挡羞羞视频在线观看| 久久久亚洲精品成人影院| 99久久精品国产国产毛片| 国产av精品麻豆| 99九九线精品视频在线观看视频| 免费人成在线观看视频色| 精品午夜福利在线看| 国产大屁股一区二区在线视频| 日本午夜av视频| 国产精品一二三区在线看| 国产黄色免费在线视频| 99热全是精品| 亚洲精品自拍成人| 51国产日韩欧美| 午夜福利高清视频| 国产精品嫩草影院av在线观看| 一级毛片黄色毛片免费观看视频| 99久久综合免费| 成人免费观看视频高清| 美女福利国产在线 | 亚洲av国产av综合av卡| 老司机影院成人| 中文天堂在线官网| 国产乱来视频区| 亚洲成人中文字幕在线播放| 欧美高清性xxxxhd video| 亚州av有码| 小蜜桃在线观看免费完整版高清| 久久影院123| 欧美老熟妇乱子伦牲交| 欧美性感艳星| 亚洲,欧美,日韩| 综合色丁香网| 伦精品一区二区三区| 免费观看性生交大片5| 黄色配什么色好看| 高清不卡的av网站| 久久女婷五月综合色啪小说| 国内精品宾馆在线| 中文天堂在线官网| 成年免费大片在线观看| av在线观看视频网站免费| 国产v大片淫在线免费观看| 精品久久久久久久久亚洲| 久久久久久人妻| 七月丁香在线播放| 亚洲第一av免费看| 噜噜噜噜噜久久久久久91| 黄色配什么色好看| av国产久精品久网站免费入址| 国产成人精品久久久久久| 日韩 亚洲 欧美在线| 久久99蜜桃精品久久| 老熟女久久久| 亚洲内射少妇av|