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

    無定河流域徑流侵蝕功率時空變化特征

    2020-04-27 06:07:22李占斌任宗萍
    水土保持研究 2020年1期

    王 偉, 李占斌,2, 楊 瑞, 賈 路, 任宗萍, 李 鵬

    (1.西安理工大學 省部共建西北旱區(qū)生態(tài)水利國家重點實驗室, 西安 710048; 2.中國科學院 水利部 水土保持研究所黃土高原土壤侵蝕與旱地農(nóng)業(yè)國家重點實驗室, 陜西 楊凌 712100; 3 延安大學 建筑工程學院, 陜西 延安 716000)

    流域尺度上,徑流不僅是泥沙的輸移媒介,而且是能量傳遞、轉(zhuǎn)化的紐帶和載體。因此,徑流侵蝕能量的傳遞過程、空間尺度效應及其與輸沙之間的響應關系是流域侵蝕研究的薄弱環(huán)節(jié),亟待對其進行定量研究[1]。研究人員[2]嘗試采用徑流能耗理論對侵蝕產(chǎn)沙過程進行了模擬和分析,得到了較好的結果。Yang C T[3]審查了常規(guī)泥沙輸送方程推導中所用的基本假設,提出徑流功率理論。同時為了進一步確定總含沙量,楊志達[4]提出了包括有關變量的函數(shù)關系。崔文濱等[5]應用這一理論對不同坡度和流量條件下的侵蝕產(chǎn)沙特征進行了分析和比較,發(fā)現(xiàn)徑流功率理論能夠較好地表征侵蝕產(chǎn)沙。以上研究說明,侵蝕能量理論能夠較好的描述土壤侵蝕。尺度問題在流域水沙模擬中具有非常重要意義,不同時空尺度的水沙響應關系十分復雜,還需要科研人員不斷地進行研究和探索。魯克新等[6-7]基于次暴雨洪水過程中徑流深、洪峰流量模數(shù)兩個水文參數(shù),提出了徑流侵蝕功率的概念,表明徑流侵蝕功率在次暴雨尺度可以更好地表征坡面和流域尺度的水力侵蝕動力。程圣東[8]首次提出了計算年功率、月功率的理論和計算方法。研究表明,在不同尺度坡面侵蝕產(chǎn)沙中,表明徑流侵蝕功率與坡面次暴雨侵蝕模數(shù)之間存在極顯著的冪函數(shù)相關關系,徑流侵蝕功率可以較好地表征坡面次暴雨水力侵蝕動力[9]。

    本文基于水蝕動力過程的徑流侵蝕能量理論,將其運用至年尺度,首先采用Mann-Kendall法和線性回歸法分析無定河流域1956—2010年年徑流量和年徑流侵蝕功率的時空變化規(guī)律;利用Mann-Kendall突變點檢驗,識別出徑流的突變年份,對比分析突變年份前后的年徑流侵蝕功率和年徑流量的時空變化規(guī)律,并建立無定河流域不同空間尺度的年徑流侵蝕功率—輸沙相關模型。

    1 研究區(qū)概況

    無定河是黃河河口至龍門區(qū)間最大的支流,發(fā)源于陜西白于山(37°00′00″—39°0′0″N,107°47′48″—110°24′21″E)。無定河流域面積30 261 km2,流域干流全長491 km,主河道平均比降為1.97%,植被覆蓋率低,墾殖指數(shù)高達60%,溝壑縱橫,地形破碎、林木稀少,水土流失嚴重,為泥沙的主要來源區(qū)[10]。

    2 數(shù)據(jù)和方法

    2.1 數(shù)據(jù)來源與處理

    本文所用的水文資料來源于黃河水文年鑒中1956—2010年的水文資料,選取無定河流域5個代表水文測站李家河站、綏德站、趙石窯站、丁家溝站和白家川站(圖1),各代表站分別控制小理河流域、大理河流域、無定河趙石窯站以上流域、無定河丁家溝站以上流域、無定河流域,控制流域面積分別為807 km2,3 893 km2,15 325 km2,23 422 km2,29 662 km2。

    圖1 無定河流域位置及水文站點分布

    2.2 研究方法

    目前,在研究徑流侵蝕方面主要有剪切力模型和侵蝕能量模型兩大類[11-14],本文采用侵蝕能量模型中的徑流侵蝕功率理論。

    通過以下公式計算出徑流侵蝕產(chǎn)沙的侵蝕動力指標:

    (1)

    E=Con·F·V

    (2)

    式中:W為徑流總量(m3);A為流域面積(m2);Qm為洪峰流量(m3/s);A′為與Qm對應的流域出口斷面的過水面積(m2);V為流域出口斷面與Qm對應的平均流速(m/s);ρ為水的密度(kg/m3);g為重力加速度(m/s2);F為作用力(N)。

    可以看出,指標E具有功率的量綱。本文將徑流年內(nèi)分配試作一次徑流過程[15],通過使用徑流侵蝕功率的計算公式,計算出年徑流侵蝕功率,從年尺度上研究流域徑流侵蝕功率的時空變化特征。

    則以下式中:各項含義分別為:

    Ea=Q′m·H

    (3)

    本文采用Mann-Kendall非參數(shù)秩次相關檢驗法[16]來進行水沙變化趨勢分析??蓸嬙旖y(tǒng)計量:

    (4)

    式中:

    (5)

    (6)

    統(tǒng)計量U稱為Kendall秩次相關系數(shù),當n增加時,U將很快收斂于標準正態(tài)分布。給定顯著性水平α,其雙尾檢驗臨界值為Uα/2,當|U|Uα/2時,系列趨勢顯著,如U>0,系列呈上升趨勢,如U<0,系列呈下降趨勢。

    線性回歸法通過建立年徑流序列yt與相應的時序t之間的線性回歸方程來檢驗時間序列的線性變化趨勢[17-18]。線性回歸方程為:

    yt=at+b

    (7)

    式中:yt為實測序列;t為時序(t=1,2,…,n;n為序列長度);a為斜率,表征時間序列的平均趨勢變化率;b為截距。

    Mann-Kendall突變點診斷法被用于進行突變年份的檢驗。當Mann-Kendall法用于檢驗序列突變性時[19],需構造一個秩序列d:

    (8)

    式中:

    (9)

    在時間序列隨機獨立的假定下,dk的均值和方差可由下邊兩個式子計算:

    (10)

    (11)

    定義統(tǒng)計變量

    (12)

    按時間序列逆序,再重復上述過程,同時使UBL=-UFK(L=n+1-k),由UFK繪制出曲線C1,由UBL繪制出曲線C2。若UFK或UBL的值超過臨界直線,表明序列上升或者下降趨勢顯著。如果C1或C2出現(xiàn)交點,且交點在臨界線之內(nèi),那么交點對應的時刻便是突變開始的時間。

    3 結果與分析

    3.1 流域水沙趨勢分析

    通過運用M-K非參數(shù)秩次相關檢驗,計算5個代表水文站徑流變化趨勢特征,見表1。可以看出,無定河流域空間內(nèi)各小流域徑流量均發(fā)生了顯著變化。徑流量秩檢驗系數(shù)都為負數(shù),徑流量呈下降趨勢,且絕對值隨流域面積逐漸增大:白家川>丁家溝>趙石窯>綏德>李家河,無定河干流水文站顯著水平0.001,支流水文站顯著水平0.05,干流下降趨勢比支流更顯著。

    表1 年徑流量變化趨勢M-K非參數(shù)秩次相關檢驗

    使用M-K檢驗法對無定河流域白家川站、丁家溝站、趙石窯站、綏德站、李家河站等5個水文站的年徑流系列進行突變點檢驗,檢驗結果見圖2。5個水文站的徑流變化呈減小趨勢,趨勢明顯超過α=0.05。由圖中選取第一個交點作為突變點進行分析,可以看出,白家川站年徑流系列突變點在1985年左右,丁家溝站年徑流系列突變點在1979年左右,趙石窯站年徑流系列突變點在1978年左右、綏德站年徑流系列突變點在1974年左右,李家河站年徑流系列突變點在1971年左右;說明大理河流域最早受人類活動影響,大理河一級支流小理河流域受人類活動影響更早。

    圖2 年徑流量M-K統(tǒng)計量時程變化

    通過統(tǒng)計無定河流域5個代表水文站的水沙特征,同時統(tǒng)計出5個水文站的年徑流量在各自突變年份前后的統(tǒng)計特征值,見表2和表3:通過表2中可以看出,無定河流域內(nèi)控制面積逐漸增大的5個水文站李家河站、綏德站、趙石窯站、丁家溝站、白家川站的年徑流量、年輸沙量的平均值、變異系數(shù)等統(tǒng)計特征都隨著水文站控制的流域面積增大而有減少,說明流域面積越小徑流與輸沙年系列變化越不穩(wěn)定,流域面積越大徑流與輸沙的年系列變化幅度越小,變化越穩(wěn)定。研究表明,在小尺度上人類活動的影響可能更加明顯、強烈。

    表2 無定河流域各水文站水沙特征

    表3 突變前后徑流量和徑流侵蝕功率

    由表3可以分析出,突變年份后無定河流域白家川站、丁家溝站、趙石窯站、綏德站、李家河站等5個水文站徑流量的平均值、變異系數(shù)等統(tǒng)計特征值基本小于突變年份之前,說明突變年份后的徑流系列更為穩(wěn)定,同時與表2對比,發(fā)現(xiàn)突變年份之前的水文要素統(tǒng)計特征指標值比總體統(tǒng)計特征指標值大,突變年份之后的水文要素統(tǒng)計特征指標值比總體統(tǒng)計特征指標值小。

    圖3 1959-2010年長時間序列年徑流侵蝕功率

    圖4 徑流侵蝕功率—輸沙關系

    3.2 流域徑流侵蝕功率的時空變化特征分析

    通過使用徑流侵蝕功率理論,計算出了無定河流域5個水文站1956—2010年的年徑流侵蝕功率,見圖3。由圖3可以看出,5個站點的徑流侵蝕功率時間序列都呈下降減小趨勢;通過統(tǒng)計發(fā)現(xiàn)白家川站、丁家溝站、趙石窯站、綏德站、李家河站的5個水文站控制流域的徑流侵蝕功率的標準差等統(tǒng)計值,在各自徑流量突變年份之前都大于突變年份后,突變年份之前流域年徑流侵蝕功率平均標準差高于突變年份后,這說明突變后的侵蝕功率更為穩(wěn)定。

    通過對白家川站、丁家溝站、趙石窯站、綏德站和李家河站5個水文站1959—2010年長時間序列徑流侵蝕功率取平均值,以及在各自突變年份之前取平均值和突變年份之后取平均值,可以分析出,無定河流域白家川站、丁家溝站、趙石窯站、綏德站和李家河站的平均徑流侵蝕功率分別為0.73,1.03,0.68,2.50,18.73。隨著控制流域尺度的增加呈遞減小趨勢,無定河支流大理河流域綏德站和小理河流域李家河站的平均徑流侵蝕功率大于干流;流域平均徑流侵蝕功率隨著控制流域尺度的增加呈遞減小趨勢。突變年份之后年平均徑流侵蝕功率比突變年份前平均減少1.05×10-4m4/(s·km2),且支流徑流侵蝕功率大于干流,干流侵蝕功率變化基本趨于穩(wěn)定,突變年份前無定河流域白家川站、丁家溝站、趙石窯站、綏德站和李家河站的年平均徑流侵蝕功率隨著流域尺度增大而逐漸減小的變化規(guī)律比較明顯,突變年份后無定河流域白家川站、丁家溝站、趙石窯站、綏德站和李家河站的年平均徑流侵蝕功率均有所減少,隨著流域尺度增大變化規(guī)律比較不明顯。

    3.3 徑流侵蝕功率—輸沙相關模型分析

    通過對無定河流域5個代表水文站的年輸沙量與年徑流侵蝕功率線性擬合,結果見圖4。根據(jù)曲線分析,無定河流域白家川站、丁家溝站、趙石窯站、綏德站、李家河站徑流侵蝕功率與輸沙量相關關系如圖4,可以看出,徑流侵蝕功率與輸沙量呈線性關系,輸沙量隨著徑流侵蝕功率的增大而增大,5個代表水文站大部分徑流侵蝕功率與輸沙量的關系點均密集分布在相關線左右,可以看出徑流侵蝕功率與輸沙量關系沒有出現(xiàn)系統(tǒng)性的偏離,線性相關方程R2分別為0.85,0.53,0.92,0.90,0.89,說明了徑流侵蝕功率理論在流域年尺度上能更好的表征徑流侵蝕輸沙關系。

    4 討 論

    徑流侵蝕功率理論最早是由魯克新等人針對岔巴溝流域次暴雨洪水尺度上提出的,該理論外推到其他時空尺度需要進一步研究。研究和表征徑流侵蝕輸沙的理論有多種,如剪切力模型和侵蝕能量模型。本文在流域年尺度上使用了徑流侵蝕功率理論,建立了無定河流域不同空間尺度的年徑流侵蝕功率—輸沙相關模型,可以較好的表征了徑流侵蝕輸沙關系,但是由于國內(nèi)外在降雨侵蝕力理論上有非常豐富的研究成果[20-29],不同時空尺度皆取得了一系列重要成果,但徑流侵蝕功率和降雨侵蝕力兩者的相關研究較為薄弱,需要進一步對其進行深入研究。

    通過突變點檢驗,可以發(fā)現(xiàn)白家川站年徑流系列有一個突變年份在1985年左右,所以將白家川站的年徑流系列分為1957—1985年,1986—2010年兩個時段進行對比分析,分析氣候變化與人類活動對無定河流域徑流變化的影響大小。

    通過表4可以看出,無定河流域?qū)崪y徑流量在1986—2010年的年均值比基準年減少4.85億m3,其中降水減少影響占-66.96%,人類活動影響減少占166.96%;分析減少的原因可以看出,在1986—2010年降雨量實際上是在增加的,水土保持等人類活動對徑流的減少具有重要的影響,從而可以看出人類活動的影響占主導地位。同時,無定河流域?qū)崪y輸沙量在1981—2010年的年均值比基準年減少1.03億t,其中降水減少影響占2.68%,人類活動影響減少占97.32%;分析減少的原因可以看出,在1981—2010年降雨量和水土保持等人類活動對輸沙的減少均具有重要的影響,但人類活動的影響占主導地位。研究表明,水土保持等人類活動可能是造成無定河流域徑流侵蝕功率變化的主導因素。

    表4 白家川站控制流域徑流輸沙變化

    5 結 論

    (1) 無定河流域1956—2010年年徑流量有顯著減小趨勢;徑流與輸沙均隨流域面積的增大流域逐漸增大;徑流突變年份在1971—1985年,突變年份后徑流離散程度變小。

    (2) 無定河流域1956—2010年年5個代表水文站年徑流侵蝕功率呈下降減小趨勢,線性回歸系數(shù)均為負數(shù),突變年份之前徑流侵蝕功率平均標準差高于突變年份后;年平均徑流侵蝕功率隨流域面積的增大而減小;突變年份之后年徑流侵蝕功率平均減少1.05×10-4m4/(s·km2),且突變年份之前徑流侵蝕功率遞減趨勢與總體保持一致。

    (3) 無定河流域5個代表水文站徑流侵蝕功率—輸沙相關模型相關關系均顯著(p<0.01),線性相關方程R2分別為0.85,0.53,0.92,0.90,0.89,徑流侵蝕功率理論在流域年尺度上能較好的表征徑流侵蝕輸沙關系。水土保持等人類活動可能是造成無定河流域徑流侵蝕功率變化的主導因素。

    一个人免费看片子| 波野结衣二区三区在线| 久久国内精品自在自线图片| 香蕉国产在线看| 赤兔流量卡办理| 午夜福利影视在线免费观看| 国产xxxxx性猛交| 天堂中文最新版在线下载| 人妻系列 视频| 只有这里有精品99| 久久这里有精品视频免费| 久久精品国产综合久久久| 国产精品一国产av| 99久国产av精品国产电影| 少妇人妻 视频| 高清视频免费观看一区二区| 99香蕉大伊视频| 波多野结衣av一区二区av| 国产精品久久久久久av不卡| 哪个播放器可以免费观看大片| 国产成人一区二区在线| 天堂俺去俺来也www色官网| 精品少妇内射三级| 超色免费av| 亚洲中文av在线| 国产日韩欧美在线精品| 久久综合国产亚洲精品| 18在线观看网站| 1024香蕉在线观看| av国产精品久久久久影院| 亚洲一码二码三码区别大吗| kizo精华| 999精品在线视频| av网站在线播放免费| 亚洲国产精品999| 妹子高潮喷水视频| 久久久国产一区二区| 日韩欧美一区视频在线观看| 亚洲婷婷狠狠爱综合网| 一本色道久久久久久精品综合| 在线看a的网站| 一个人免费看片子| 国产精品久久久久久精品电影小说| 亚洲av福利一区| 交换朋友夫妻互换小说| 赤兔流量卡办理| 伊人久久大香线蕉亚洲五| 亚洲 欧美一区二区三区| 日本色播在线视频| 国产av码专区亚洲av| 国产精品麻豆人妻色哟哟久久| 国产老妇伦熟女老妇高清| 日韩大片免费观看网站| 一级片免费观看大全| 18+在线观看网站| 妹子高潮喷水视频| 青青草视频在线视频观看| 国产 一区精品| 国产有黄有色有爽视频| 亚洲av综合色区一区| 国产精品欧美亚洲77777| 精品福利永久在线观看| 欧美97在线视频| 亚洲成人手机| 91国产中文字幕| 国产老妇伦熟女老妇高清| 久久久久网色| 美女视频免费永久观看网站| 女人高潮潮喷娇喘18禁视频| 黄网站色视频无遮挡免费观看| 免费在线观看视频国产中文字幕亚洲 | 日韩制服骚丝袜av| 老鸭窝网址在线观看| 精品国产露脸久久av麻豆| 久久av网站| 国产精品久久久av美女十八| 欧美 日韩 精品 国产| 黄色视频在线播放观看不卡| 春色校园在线视频观看| 我的亚洲天堂| 91久久精品国产一区二区三区| 七月丁香在线播放| 日韩三级伦理在线观看| 国产成人91sexporn| 只有这里有精品99| 18在线观看网站| 精品国产露脸久久av麻豆| 欧美激情极品国产一区二区三区| 亚洲,欧美,日韩| kizo精华| 国产精品香港三级国产av潘金莲 | 丰满乱子伦码专区| 2018国产大陆天天弄谢| 777米奇影视久久| 欧美另类一区| 日本av手机在线免费观看| 久久影院123| 黑人猛操日本美女一级片| 国产激情久久老熟女| 桃花免费在线播放| 成人国语在线视频| 欧美中文综合在线视频| 久久久精品94久久精品| √禁漫天堂资源中文www| 99久久中文字幕三级久久日本| 久久久久国产一级毛片高清牌| 日本wwww免费看| 搡女人真爽免费视频火全软件| 精品一区二区免费观看| 99热国产这里只有精品6| 99香蕉大伊视频| 国产精品秋霞免费鲁丝片| 国产精品成人在线| 亚洲熟女精品中文字幕| 99国产综合亚洲精品| 天堂中文最新版在线下载| 免费高清在线观看日韩| 丝袜美腿诱惑在线| 最新的欧美精品一区二区| 老鸭窝网址在线观看| 男女下面插进去视频免费观看| 国产av码专区亚洲av| 99精国产麻豆久久婷婷| 欧美最新免费一区二区三区| 女人久久www免费人成看片| 国产探花极品一区二区| 人妻少妇偷人精品九色| 久久综合国产亚洲精品| 国产精品欧美亚洲77777| 各种免费的搞黄视频| 成人国语在线视频| 女人精品久久久久毛片| 美女国产视频在线观看| 国产精品秋霞免费鲁丝片| 又大又黄又爽视频免费| 久久久久精品久久久久真实原创| 美女主播在线视频| av又黄又爽大尺度在线免费看| www.精华液| 黄色怎么调成土黄色| 9191精品国产免费久久| 欧美国产精品va在线观看不卡| 亚洲国产精品成人久久小说| 亚洲一区中文字幕在线| 最近的中文字幕免费完整| 青春草国产在线视频| 天天躁夜夜躁狠狠躁躁| 精品亚洲成a人片在线观看| 久久久a久久爽久久v久久| 秋霞伦理黄片| a级片在线免费高清观看视频| kizo精华| 大陆偷拍与自拍| 制服诱惑二区| 国产亚洲最大av| 91久久精品国产一区二区三区| 国产1区2区3区精品| 日韩免费高清中文字幕av| 亚洲五月色婷婷综合| 999精品在线视频| 日日啪夜夜爽| 国产成人精品一,二区| 免费看不卡的av| 日本黄色日本黄色录像| 五月开心婷婷网| 免费人妻精品一区二区三区视频| 久久久亚洲精品成人影院| 最黄视频免费看| 国产精品亚洲av一区麻豆 | a级毛片在线看网站| av在线观看视频网站免费| 亚洲美女视频黄频| 十分钟在线观看高清视频www| 王馨瑶露胸无遮挡在线观看| 精品99又大又爽又粗少妇毛片| 亚洲视频免费观看视频| 日韩av不卡免费在线播放| 亚洲国产欧美在线一区| 久久精品国产综合久久久| 伊人久久国产一区二区| 日韩电影二区| 免费看不卡的av| 在线天堂中文资源库| 亚洲欧美精品自产自拍| 极品人妻少妇av视频| 最近中文字幕高清免费大全6| 又黄又粗又硬又大视频| 韩国高清视频一区二区三区| 在线观看一区二区三区激情| 亚洲,欧美,日韩| 日韩一本色道免费dvd| 在线亚洲精品国产二区图片欧美| 亚洲精品成人av观看孕妇| 妹子高潮喷水视频| 亚洲欧美色中文字幕在线| 两个人免费观看高清视频| 免费大片黄手机在线观看| 老司机影院毛片| 电影成人av| 美女国产视频在线观看| 久久影院123| 久久久久久久精品精品| 啦啦啦中文免费视频观看日本| 国产精品久久久久久av不卡| 又黄又粗又硬又大视频| 侵犯人妻中文字幕一二三四区| 国产精品免费视频内射| 麻豆乱淫一区二区| 一区二区三区激情视频| 狠狠精品人妻久久久久久综合| 老女人水多毛片| 午夜免费男女啪啪视频观看| 日韩,欧美,国产一区二区三区| 美女脱内裤让男人舔精品视频| 亚洲欧美一区二区三区黑人 | 日韩成人av中文字幕在线观看| 国产成人a∨麻豆精品| 国产精品熟女久久久久浪| 欧美成人午夜精品| 黑人巨大精品欧美一区二区蜜桃| 巨乳人妻的诱惑在线观看| 国产精品国产av在线观看| 老汉色av国产亚洲站长工具| 国产av码专区亚洲av| 免费人妻精品一区二区三区视频| 国产麻豆69| 2018国产大陆天天弄谢| 精品少妇久久久久久888优播| 夜夜骑夜夜射夜夜干| 日韩 亚洲 欧美在线| 国产一区二区三区综合在线观看| 久久鲁丝午夜福利片| 国产男女超爽视频在线观看| 成人毛片a级毛片在线播放| 亚洲色图 男人天堂 中文字幕| 久久久久视频综合| 日本vs欧美在线观看视频| 在线观看美女被高潮喷水网站| 曰老女人黄片| av在线app专区| 欧美日韩亚洲高清精品| 超碰97精品在线观看| 免费少妇av软件| 日本黄色日本黄色录像| 一级,二级,三级黄色视频| av免费在线看不卡| 久久久久久久精品精品| 午夜激情av网站| 亚洲伊人色综图| 日本黄色日本黄色录像| 寂寞人妻少妇视频99o| 91aial.com中文字幕在线观看| 人人澡人人妻人| 国产精品久久久久久精品电影小说| videos熟女内射| 亚洲成av片中文字幕在线观看 | 91aial.com中文字幕在线观看| 人人澡人人妻人| 亚洲图色成人| 看十八女毛片水多多多| 成人国产av品久久久| 国产在视频线精品| 黄色 视频免费看| 日韩在线高清观看一区二区三区| 欧美在线黄色| 在线天堂最新版资源| 成年人免费黄色播放视频| 亚洲欧美精品自产自拍| 精品一区二区三区四区五区乱码 | 99久久综合免费| 久久久久国产精品人妻一区二区| 丝袜喷水一区| 久久精品国产亚洲av涩爱| 中文天堂在线官网| 精品国产一区二区久久| 午夜久久久在线观看| 王馨瑶露胸无遮挡在线观看| 精品国产超薄肉色丝袜足j| 国产成人精品一,二区| 国产成人精品在线电影| 日韩一本色道免费dvd| 大片电影免费在线观看免费| 亚洲国产最新在线播放| 麻豆乱淫一区二区| 少妇精品久久久久久久| 亚洲人成77777在线视频| 亚洲av国产av综合av卡| av片东京热男人的天堂| 夜夜骑夜夜射夜夜干| 欧美成人午夜免费资源| 中文字幕人妻丝袜一区二区 | 久久久久国产网址| freevideosex欧美| 日韩视频在线欧美| 国产福利在线免费观看视频| 久久这里只有精品19| 高清不卡的av网站| 超碰97精品在线观看| 亚洲,欧美精品.| 午夜日本视频在线| 国产精品 欧美亚洲| 精品少妇一区二区三区视频日本电影 | 视频在线观看一区二区三区| 十分钟在线观看高清视频www| 免费观看无遮挡的男女| 国产精品久久久久久久久免| 黄片小视频在线播放| 永久网站在线| 视频区图区小说| 免费观看性生交大片5| 免费黄色在线免费观看| av在线老鸭窝| 九草在线视频观看| 最新的欧美精品一区二区| 国产成人精品婷婷| 欧美人与性动交α欧美软件| 男人舔女人的私密视频| 欧美+日韩+精品| 日韩av在线免费看完整版不卡| 色吧在线观看| av国产久精品久网站免费入址| 亚洲av综合色区一区| 少妇的逼水好多| 国产成人精品无人区| av免费在线看不卡| 国产免费视频播放在线视频| 婷婷成人精品国产| 国产又色又爽无遮挡免| 天天操日日干夜夜撸| 国产97色在线日韩免费| 黄色毛片三级朝国网站| 国产成人欧美| 宅男免费午夜| 国产精品国产三级专区第一集| 一级片'在线观看视频| 中文字幕精品免费在线观看视频| 国产野战对白在线观看| 国产探花极品一区二区| 18在线观看网站| 久久ye,这里只有精品| 最近2019中文字幕mv第一页| 免费少妇av软件| 狠狠婷婷综合久久久久久88av| 免费播放大片免费观看视频在线观看| 久久综合国产亚洲精品| 啦啦啦啦在线视频资源| 另类亚洲欧美激情| 午夜日本视频在线| 丰满乱子伦码专区| 国产精品免费视频内射| 亚洲精品一二三| 人成视频在线观看免费观看| 亚洲一级一片aⅴ在线观看| 9191精品国产免费久久| av线在线观看网站| 90打野战视频偷拍视频| 免费在线观看黄色视频的| 国产av国产精品国产| 欧美 日韩 精品 国产| 91精品三级在线观看| 99香蕉大伊视频| 国产1区2区3区精品| 女人高潮潮喷娇喘18禁视频| 国产有黄有色有爽视频| 菩萨蛮人人尽说江南好唐韦庄| 国产精品99久久99久久久不卡 | 亚洲国产精品国产精品| 亚洲成av片中文字幕在线观看 | 我的亚洲天堂| 韩国高清视频一区二区三区| 亚洲成色77777| 国产 一区精品| 一区二区三区四区激情视频| 亚洲少妇的诱惑av| 黄频高清免费视频| 久久久国产一区二区| 水蜜桃什么品种好| 80岁老熟妇乱子伦牲交| 亚洲av中文av极速乱| 中国国产av一级| 久久人人97超碰香蕉20202| 成人亚洲欧美一区二区av| 欧美精品国产亚洲| 免费黄色在线免费观看| 欧美+日韩+精品| 日韩精品有码人妻一区| 亚洲综合精品二区| 天天躁狠狠躁夜夜躁狠狠躁| 午夜福利在线观看免费完整高清在| 国产精品秋霞免费鲁丝片| 91精品伊人久久大香线蕉| 亚洲av在线观看美女高潮| 日本免费在线观看一区| 中文精品一卡2卡3卡4更新| 99热全是精品| 国产精品免费大片| 国产成人av激情在线播放| 香蕉精品网在线| 日本-黄色视频高清免费观看| 欧美另类一区| 欧美 日韩 精品 国产| 国产精品一区二区在线观看99| 精品国产一区二区久久| 好男人视频免费观看在线| 久久精品熟女亚洲av麻豆精品| 亚洲熟女精品中文字幕| 婷婷色综合www| tube8黄色片| 亚洲激情五月婷婷啪啪| 91在线精品国自产拍蜜月| 熟妇人妻不卡中文字幕| 亚洲成人av在线免费| 亚洲欧洲日产国产| 天堂俺去俺来也www色官网| 最新的欧美精品一区二区| 亚洲第一区二区三区不卡| 大片电影免费在线观看免费| 国产片内射在线| 国产黄色视频一区二区在线观看| 熟女av电影| 免费在线观看黄色视频的| 最近手机中文字幕大全| 80岁老熟妇乱子伦牲交| 美女主播在线视频| 99久久精品国产国产毛片| 欧美+日韩+精品| 中文精品一卡2卡3卡4更新| 天天躁夜夜躁狠狠躁躁| 男女啪啪激烈高潮av片| 人人妻人人添人人爽欧美一区卜| 久久久久久久大尺度免费视频| 久久精品人人爽人人爽视色| 免费不卡的大黄色大毛片视频在线观看| 高清欧美精品videossex| 国产淫语在线视频| 亚洲天堂av无毛| 欧美日韩国产mv在线观看视频| 日韩一区二区视频免费看| av免费观看日本| 国语对白做爰xxxⅹ性视频网站| 一区二区日韩欧美中文字幕| 亚洲天堂av无毛| av不卡在线播放| 国产成人一区二区在线| 日韩大片免费观看网站| 一边摸一边做爽爽视频免费| 免费观看在线日韩| 最近中文字幕高清免费大全6| 精品人妻一区二区三区麻豆| 这个男人来自地球电影免费观看 | 亚洲欧洲国产日韩| 下体分泌物呈黄色| 国产一区二区三区综合在线观看| 一本色道久久久久久精品综合| 亚洲成人一二三区av| 1024视频免费在线观看| 好男人视频免费观看在线| 成人免费观看视频高清| 只有这里有精品99| 久久韩国三级中文字幕| 久久久a久久爽久久v久久| 国产黄色免费在线视频| 精品少妇内射三级| 久久久久久久亚洲中文字幕| 男女边摸边吃奶| 精品少妇久久久久久888优播| 一区二区日韩欧美中文字幕| 啦啦啦在线观看免费高清www| 国产精品久久久久久久久免| 亚洲欧美精品自产自拍| 国产人伦9x9x在线观看 | 视频区图区小说| 激情视频va一区二区三区| 久久免费观看电影| 777久久人妻少妇嫩草av网站| 制服丝袜香蕉在线| 又黄又粗又硬又大视频| 亚洲男人天堂网一区| 最新中文字幕久久久久| www日本在线高清视频| 人人妻人人爽人人添夜夜欢视频| 一本—道久久a久久精品蜜桃钙片| 在线 av 中文字幕| 高清视频免费观看一区二区| 大片电影免费在线观看免费| 久久精品久久久久久久性| 国产男女超爽视频在线观看| 欧美中文综合在线视频| 三上悠亚av全集在线观看| 免费大片黄手机在线观看| av天堂久久9| 国产一区有黄有色的免费视频| 如日韩欧美国产精品一区二区三区| 不卡av一区二区三区| 亚洲精华国产精华液的使用体验| 中文精品一卡2卡3卡4更新| 国产熟女午夜一区二区三区| 亚洲欧洲精品一区二区精品久久久 | 国产伦理片在线播放av一区| 亚洲精品乱久久久久久| 熟妇人妻不卡中文字幕| 国产黄色免费在线视频| 高清欧美精品videossex| 精品酒店卫生间| 国产成人av激情在线播放| 日韩电影二区| 亚洲av男天堂| 亚洲av在线观看美女高潮| 日韩精品有码人妻一区| 亚洲成人一二三区av| 亚洲欧美一区二区三区黑人 | 精品一品国产午夜福利视频| 日本vs欧美在线观看视频| 久久久久久久大尺度免费视频| 天天操日日干夜夜撸| 日日爽夜夜爽网站| 精品人妻熟女毛片av久久网站| 亚洲av电影在线观看一区二区三区| av视频免费观看在线观看| 成人影院久久| 久久亚洲国产成人精品v| av一本久久久久| 老女人水多毛片| 一区二区av电影网| 天堂俺去俺来也www色官网| 天天操日日干夜夜撸| 超碰成人久久| 在线天堂最新版资源| 国产成人a∨麻豆精品| 精品人妻熟女毛片av久久网站| 欧美中文综合在线视频| 一级爰片在线观看| 亚洲av男天堂| 天天躁夜夜躁狠狠躁躁| 18在线观看网站| 亚洲国产色片| 在线看a的网站| 亚洲欧美清纯卡通| 亚洲人成77777在线视频| 91在线精品国自产拍蜜月| 免费在线观看完整版高清| 色视频在线一区二区三区| 亚洲精华国产精华液的使用体验| av又黄又爽大尺度在线免费看| www.自偷自拍.com| 国产精品一国产av| 在线观看www视频免费| 午夜日本视频在线| 在现免费观看毛片| 久久久久久久久久久免费av| 99久久综合免费| 黄色一级大片看看| av国产精品久久久久影院| 欧美老熟妇乱子伦牲交| 多毛熟女@视频| 久久精品熟女亚洲av麻豆精品| 可以免费在线观看a视频的电影网站 | 成人国产av品久久久| 一区二区av电影网| 成年av动漫网址| 精品一区二区三区四区五区乱码 | 夫妻性生交免费视频一级片| 亚洲 欧美一区二区三区| 黄片小视频在线播放| 免费久久久久久久精品成人欧美视频| 久久国产精品男人的天堂亚洲| 日韩一卡2卡3卡4卡2021年| 亚洲美女搞黄在线观看| 69精品国产乱码久久久| 久久久久久免费高清国产稀缺| 天天躁夜夜躁狠狠久久av| 久久久久网色| 老鸭窝网址在线观看| 国产成人aa在线观看| 大话2 男鬼变身卡| 欧美国产精品va在线观看不卡| 中文字幕人妻丝袜制服| 丝袜美足系列| 亚洲人成77777在线视频| av在线播放精品| 国产成人精品福利久久| 久久影院123| 国产精品不卡视频一区二区| 亚洲av电影在线观看一区二区三区| 国产淫语在线视频| 国产亚洲欧美精品永久| 久久狼人影院| 亚洲内射少妇av| 日韩一区二区视频免费看| 一边亲一边摸免费视频| 欧美精品高潮呻吟av久久| 亚洲国产欧美在线一区| 日韩电影二区| 男女国产视频网站| 在线观看免费视频网站a站| 一级毛片电影观看| 亚洲,欧美精品.| 最近中文字幕高清免费大全6| 老汉色av国产亚洲站长工具| 人人妻人人添人人爽欧美一区卜| 中文字幕人妻丝袜一区二区 | 国产精品人妻久久久影院| 精品亚洲成国产av| 日韩一卡2卡3卡4卡2021年| 精品亚洲成a人片在线观看| 日韩大片免费观看网站| 高清在线视频一区二区三区| 日本vs欧美在线观看视频| 高清在线视频一区二区三区| 久久久久视频综合| 欧美97在线视频| 国产激情久久老熟女| 国产又色又爽无遮挡免| 久久97久久精品| 黑丝袜美女国产一区|