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

    沭河流域季節(jié)實(shí)測徑流變化特性分析

    2020-06-12 08:11:44陸建宇張育德
    海河水利 2020年3期
    關(guān)鍵詞:沭河年際徑流量

    陸建宇,張育德,杜 娟

    (中水北方勘測設(shè)計(jì)研究有限責(zé)任公司,天津 300222)

    河川徑流是一種重要的水資源,徑流量的變化將直接影響水資源的開發(fā)利用[1]。在氣候與人類影響下,近些年來,我國各大流域河川徑流特性均已發(fā)生了不同程度的改變[2],這不僅會(huì)影響到人類社會(huì)的健康發(fā)展,還會(huì)給區(qū)域工農(nóng)業(yè)生產(chǎn)、水資源管理及生態(tài)環(huán)境等帶來諸多不利影響[3]。因此,分析實(shí)測徑流變化特性具有重要的現(xiàn)實(shí)意義。

    沭河屬淮河流域沂沭泗水系,汛期雨強(qiáng)較大,易發(fā)生區(qū)域性洪澇災(zāi)害[4]。目前,關(guān)于沭河徑流方面的研究成果已有一些報(bào)道。其中,郝敏[5]等研究了沭河徑流年際變化規(guī)律及年內(nèi)分配特征,發(fā)現(xiàn)年徑流量存在減小趨勢,徑流年內(nèi)分配漸漸趨于平均;張魯[6]探討了沂沭河流域洪水管理措施和洪水資源開發(fā)利用方案,擬定了經(jīng)大官莊水利樞紐向外流域的引水規(guī)模與調(diào)水方案;陸建宇[7]等采用多種生態(tài)徑流計(jì)算方法,估算了沭河流域生態(tài)徑流與生態(tài)需水,發(fā)現(xiàn)多數(shù)時(shí)段內(nèi)的河流生態(tài)需水都可被滿足。此外,盡管還有一些其他方面的研究成果[8,9],但涉及徑流變化特性分析的成果仍然很少。鑒于此,筆者嘗試以流域控制站——大官莊站為研究對象,借助多種時(shí)間統(tǒng)計(jì)及分析方法,著重研究該流域春、夏、秋、冬季實(shí)測徑流變化特性,以期為合理規(guī)劃與利用流域水資源提供一定的科學(xué)指導(dǎo)。

    1 沭河流域概況

    沭河流域位于山東省南部及江蘇省北部,是山東省臨沂市的母親河[7]。大官莊水文站位于沭河流域上游,是流域主要控制性站點(diǎn),控制流域面積4 529 km2,干流河長 205.5 km[5]。沭河在大官莊站以下分為2 支,即新沭河、老沭河,其中新沭河上建有節(jié)制閘、老沭河上筑有人民勝利堰,所以大官莊站處流量為以上2處閘堰流量之和[7];大官莊站以上基本是山丘區(qū),流域汛期洪水經(jīng)常表現(xiàn)出漲落迅速、量大峰高的特征;大官莊站以下為狹長平坦的帶狀河谷平原。

    沭河流域?yàn)榈湫捅睖貛Т箨懶约撅L(fēng)氣候,四季特征明顯,春季風(fēng)大干燥,夏季濕熱多雨,秋季舒適涼爽,冬季寒冷少雨。該流域常年降水比較豐富,年均降水量806.0 mm,主要分布在6—9 月,占全年的70%以上,7 月是年內(nèi)最大暴雨易發(fā)月份,其月降水量占到全年的30%左右[6];降水是沭河河川徑流的主要來源,降水時(shí)空分布和地形起伏變化直接決定了流域水資源的年際、年內(nèi)的豐枯變化,同時(shí)山洪極易暴發(fā)在夏、秋兩季,這十分不利于流域水資源綜合規(guī)劃、有序開發(fā)和高效利用[5]。

    2 資料及方法

    本次搜集并整理了沭河流域內(nèi)新沭河節(jié)制閘、老沭河人民勝利堰的實(shí)測徑流數(shù)據(jù),匯總得到沭河流域主要控制站大官莊站逐月實(shí)測徑流系列,研究時(shí)段為1958—2016年,總計(jì)59 a。

    2.1 累積濾波器法

    累積濾波器法的關(guān)鍵在于累積均值,該方法能夠定性地反映實(shí)測徑流序列的變化趨勢[10],其原理如下:

    式中:Sk為第k年累積均值;xi為第i年值;xˉ為多年均值;k為序列長度,k=1,2,3,…,n。

    2.2 Man-Kendall突變檢驗(yàn)法

    Man-Kendall 法是一種非參數(shù)統(tǒng)計(jì)檢驗(yàn)與分析方法[11],該方法不需要樣本數(shù)據(jù)服從一定的分布特征,可用于徑流序列的定量化趨勢分析和突變點(diǎn)檢驗(yàn)。在檢測徑流序列突變點(diǎn)時(shí),需構(gòu)造ri(i=1,2,···,n),統(tǒng)計(jì)量Sn、UFn

    的計(jì)算公式為:

    式中:mi為rj>ri(2 ≤i≤j)的累積數(shù);n為檢驗(yàn)樣本系列長度;N為樣本系列總長度。

    同時(shí),逆序排列ri(i=1,2,···n),按式(3)重復(fù)上述步驟,即得UBn':

    在檢測突變點(diǎn)時(shí),若UFn和UBn'

    在a=0.01置信區(qū)間內(nèi)有交點(diǎn),則交點(diǎn)對應(yīng)的時(shí)間即為實(shí)測徑流序列的突變點(diǎn)[12]。

    2.3 Morlet小波分析法

    Morlet 小波是高斯包絡(luò)下的單頻率復(fù)正弦函數(shù)[13],其波形與徑流序列變化過程相近,因而常借助Morlet 小波識別徑流序列的多時(shí)間尺度特征。Mor?let小波函數(shù)形式如下:

    式中:?(t)為子小波;i為虛數(shù);w0為常數(shù),常取6.2。

    同時(shí),小波方差圖中Var(a)的峰值對應(yīng)時(shí)間尺度b的數(shù)值,可確定徑流序列的主周期[13],相應(yīng)函數(shù)形式如下:

    式中:Wf(a,b)為小波變換系數(shù);a,b分別為尺度參數(shù),時(shí)間參數(shù)。

    3 計(jì)算結(jié)果與分析

    3.1 季節(jié)徑流年際過程

    沭河流域季節(jié)實(shí)測徑流量年際過程(如圖1 所示)及其年代均值過程(見表1)顯示:①年際變化上,大官莊站春、夏、秋、冬季實(shí)測徑流量年際波動(dòng)較大,多年平均徑流量依次為6 718 萬、61 922 萬、23 291 萬和7 883 萬m3,對應(yīng)徑流量極值比分別為1 021.4、41.1、95.3 和 51.6。同時(shí),各季徑流量分別以445.86萬、-8 960.9萬、793.58萬、1 301.9萬m3/10 a的線性速度呈增加、減少、增加、增加變化,其中夏季徑流減少、冬季徑流變幅較顯著。②年代變化上,大官莊站春季各年代實(shí)測徑流量在波動(dòng)中總體呈增加態(tài)勢,由1950 s的4 511萬m3增加到2 010 s的9 962萬m3;與1960、1970 s相比,夏季1980—2010 s實(shí)測徑流量年代均值大幅減少,其中1980 s最小,僅為29 472萬m3;秋季實(shí)測徑流量年代變化上為“增-減-增-減”過程,其中1980 s最小僅為10 322萬m3,2000 s最大達(dá)37 474 萬m3,各年代時(shí)段徑流量極值比為3.63;1960—1980 s,冬季徑流年代均值由10 415 萬m3一直持續(xù)減少到1 910萬m3,而在1980 s以后又持續(xù)增加,增至2010 s的16 110萬m3。整體上,沭河流域各季節(jié)年際、年代實(shí)測徑流量相差較大,現(xiàn)階段春、冬季實(shí)測徑流較多年均值偏豐,夏季偏枯,秋季則與多年均值持平,各季徑流總體變化不平穩(wěn)。

    表1 沭河流域季節(jié)實(shí)測徑流量年代均值統(tǒng)計(jì) 萬m3

    圖1 沭河流域季節(jié)實(shí)測徑流量年際變化過程

    3.2 季節(jié)徑流變化趨勢分析

    基于沭河大官莊站1958—2016 年逐月實(shí)測徑流資料,繪制各季實(shí)測徑流量的累積均值過程線,如圖2 所示。由圖2 可看出:①春季實(shí)測徑流量總體表現(xiàn)為增加過程,1958—1964 年急劇增加,1965—1971 年大幅減少,之后在連續(xù)波動(dòng)中維持小幅增加,但增加趨勢相對較弱;②夏季實(shí)測徑流量在1960—1966 年間不斷增加,雖然經(jīng)歷了1967—1971 年的短期減少和1972—1975 年的暫時(shí)增加,但1975 年以后又保持顯著減少變化;③秋季實(shí)測徑流量總體表現(xiàn)為增加變化,其中1958—1965 年為增加過程、1966—1970 和 1971—1992 年分別為增加和減少變化、1993 年以來總體為增加態(tài)勢;4 冬季實(shí)測徑流量總體呈增加趨勢,1958—1962 年急劇增加,1963—1968 年大幅減少,之后在小幅突增后持續(xù)減少,但2003 年以后轉(zhuǎn)為明顯增加態(tài)勢。總體上,沭河流域季節(jié)徑流量具有明顯的豐枯階段性變化特征,同時(shí)其累積均值過程能較好地體現(xiàn)各季徑流量年際變化趨勢和階段性特征。

    圖2 沭河流域季節(jié)實(shí)測徑流量累積均值過程

    為定量分析流域季節(jié)實(shí)測徑流量變化趨勢的顯著水平,借助Man-Kendall 趨勢檢驗(yàn)法[14]對各季徑流量做趨勢檢驗(yàn),并將其檢驗(yàn)結(jié)果與累積濾波器法分析對比,結(jié)果見表2。

    表2 沭河流域季節(jié)實(shí)測徑流量變化趨勢分析

    從表2 可以看出,大官莊站季節(jié)徑流量年際波動(dòng)十分顯著,Cv值均維持在0.72 以上,同時(shí)Man-Kendall 趨勢檢驗(yàn)法的檢驗(yàn)結(jié)果與累積濾波器法的分析結(jié)果一致。其中,春季徑流量的統(tǒng)計(jì)量Z值為1.007,表現(xiàn)為不顯著增加變化;夏季徑流量的Z值為-1.936,存在顯著減少變化(90%顯著性水平);秋季徑流量的Z值為0.209,呈不顯著增加趨勢;冬季徑流量的Z值為2.001,表現(xiàn)為十分顯著的增加態(tài)勢(95%顯著性水平)。

    3.3 季節(jié)徑流Man-Kendall法突變分析

    沭河流域季節(jié)徑流量Man-Kendall 法突變檢驗(yàn)結(jié)果,如圖3所示。由圖3可以看出:①在a=0.05臨界區(qū)間內(nèi),大官莊站季節(jié)實(shí)測徑流量的UF、UB存在多個(gè)相交點(diǎn),即徑流突變點(diǎn),部分位于1980年以前,這主要是由在沭河流域內(nèi)實(shí)施“導(dǎo)沭整沂”工程和“蓄泄兼籌”工程期間在流域上游修筑大量水庫及其他引水、分洪、調(diào)水等水利工程和水土保持工程引起的[15]。②春季徑流量突變點(diǎn)位于1959、1961、1963、2004 年,夏季位于1966、1969、1971、2006、2009 年,秋季位于1962、1964、1999、2002 年,冬季位于2003年;沭河上游已建成大型水庫4座,其中沂水縣沙溝水庫總庫容1.04 億m3、1959 年蓄水,莒縣青峰嶺水庫總庫容4.101億m3、1960年蓄水,莒縣小仕陽水庫總庫容1.254億m3、1958—1959年蓄水,莒南縣陡山水庫總庫容 2.92 億 m3、1959 年 7 月竣工;結(jié)合文獻(xiàn)[16]可知,各雨量站的長系列實(shí)測年降水中,年最大、年最小降水量多出現(xiàn)在2001—2003年。③1951—1966年,春季徑流量的UF>0,說明該時(shí)段徑流量有增加趨勢,1967—1977年短暫為負(fù)值,1972—1980年轉(zhuǎn)為正值,之后持續(xù)保持負(fù)值,2008 年以后以正值持續(xù)至2016 年,說明春季徑流量增減變化十分頻繁;1960—1966年,夏季的UF>0,說明該時(shí)段徑流量呈短暫增加變化,之后一直為負(fù)值,持續(xù)至2016年,說明夏季徑流量又持續(xù)減少;1960—1965年,秋、冬季的UF>0,說明該時(shí)段徑流量暫時(shí)呈增加趨勢,之后一直為負(fù)值,持續(xù)至2007年,說明秋、冬季徑流量又持續(xù)減少,2007 年以后以正值持續(xù)至2016 年,說明秋、冬季徑流量近期持續(xù)增加。整體上,沭河流域?qū)崪y徑流突變特征受水利工程、降水變化共同影響,各季徑流量的UF的正、負(fù)變化相對頻繁,但僅夏冬季實(shí)測徑流量的UF長期超出a=0.05的臨界區(qū)間,說明夏冬季實(shí)測徑流量波動(dòng)較大、趨勢特征較為顯著,這也與前述趨勢分析結(jié)果是相符的。

    圖3 沭河流域季節(jié)實(shí)測徑流量Man-Kendall法突變檢驗(yàn)

    3.4 季節(jié)徑流Morlet小波周期分析

    圖4—5 為沭河流域季節(jié)實(shí)測徑流量Morlet 小波實(shí)部等值線圖和小波方差圖。圖4中實(shí)線表示的小波實(shí)部系數(shù)為正,說明對應(yīng)徑流量屬偏多狀態(tài);虛線表示的小波實(shí)部系數(shù)為負(fù),說明對應(yīng)徑流量屬偏少狀態(tài);黑色實(shí)線加粗部分對應(yīng)著徑流量正發(fā)生豐枯轉(zhuǎn)變[11]。由圖4—5可以看出:①春季徑流量存在4、12、17 a 3個(gè)震蕩周期,其中最劇烈的時(shí)間尺度為12 a,主要分布在1960—2000 年;夏季徑流量存在4、12 a 2 個(gè)周期特征,其中12 a 為最主要的震蕩周期;秋季徑流量存在2、5、11 a 3個(gè)層次的周期變化,其中11 a 尺度最穩(wěn)定,且經(jīng)歷了多個(gè)“枯—豐—枯”的循環(huán)交替過程;冬季徑流量存在7、20 a 2 個(gè)時(shí)間尺度,其中7 a為主周期。②各季徑流量具有十分明顯的階段性變化。在各季實(shí)測徑流主周期上,春季表現(xiàn)為 3 個(gè)豐水段即 1960—1964、1972—1978 和2006—2012 年,3 個(gè)枯水段即 1965—1971、1979—1985 和 2000—2005 年;夏季豐水段為 1960—1966、1972—1975 年,枯水期體現(xiàn)在1967—1971、1976—1982、2001—2005 年;秋季豐水段為 1962—1965、1971—1975、1993—1997 年,枯水段為1966—1970、1998—2003 年;冬季豐水段為1964—1967、1968—1971、1975—1977、2004—2006、2010—2013 年 ,枯 水 段 為 1958—1963、1972—1974、1977—1980、2000—2003、2007—2009 年。2016 年以后,春季實(shí)測徑流量較少,夏季徑流量相對偏枯,秋季徑流量較少,冬季徑流量相對偏枯,近期年徑流量也將處于偏少狀態(tài)。

    圖4 沭河流域季節(jié)實(shí)測徑流量Morlet小波實(shí)部等值線

    圖5 沭河流域季節(jié)實(shí)測徑流量Morlet小波方差

    4 小結(jié)

    基于多種時(shí)間序列分析方法,采用大官莊站1958—2016年逐月實(shí)測徑流資料,分析了沭河流域各季實(shí)測徑流量變化特性,得出主要結(jié)論如下。

    (1)沭河流域春、夏、秋和冬季實(shí)測徑流主周期分別為12、12、11 和7 a,同時(shí)各季徑流量存在顯著階段性特征。其中,春季徑流量在1960—1964、1972—1978、2006—2012 年偏豐,在 1965—1971、1979—1985、2000—2005 年偏枯;夏季徑流量在1960—1966、1972—1975 年偏多,在 1967—1971、1976—1982、2001—2005 年偏少;秋季豐水段為1962—1965、1971—1975、1993—1997 年,枯水段為1966—1970、1998—2003 年;冬季徑流量在1964—1967、1968—1971、1975—1977、2004—2006、2010—2013 年豐沛,在 1958—1963、1972—1974、1977—1980、2000—2003、2007—2009年相對偏少。

    (2)沭河流域春、夏、秋、冬季實(shí)測徑流量依次以445.86 萬、-8 960.9 萬、793.58 萬、1 301.9 萬 m3/10 a的速率呈增加、減少、增加、增加變化;春季實(shí)測徑流量突變點(diǎn)依次位于1959、1961、1963、2004 年,夏季位于 1966、1969、1971、2006、2009 年,秋季位于1962、1964、1999、2002年,冬季位于2003年;上述研究表明,流域各季實(shí)測徑流量變化十分復(fù)雜,現(xiàn)階段春、冬季實(shí)測徑流量較多年均值偏豐,夏季偏枯,秋季與多年均值持平;同時(shí),夏季實(shí)測徑流量的減少雖然可能會(huì)降低流域山洪災(zāi)害沖擊強(qiáng)度及影響范圍,但也將在一定程度上增加流域工農(nóng)業(yè)生產(chǎn)和人們生活的用水壓力,致使區(qū)域水資源供需矛盾更加尖銳和突出。

    猜你喜歡
    沭河年際徑流量
    美麗的沭河公園
    我家門前沭水流
    北緯30°中層頂區(qū)域鈉與鐵原子層的結(jié)構(gòu)和年際變化
    蘆 葦
    水文比擬法在計(jì)算河川徑流量時(shí)的修正
    亞洲夏季風(fēng)的年際和年代際變化及其未來預(yù)測
    與北大西洋接壤的北極海冰和年際氣候變化
    SCS模型在紅壤土坡地降雨徑流量估算中的應(yīng)用
    老沭河口頭站二0一二年水情分析
    基于M-K法對圖們江干流含沙量年際變化的分析
    www.av在线官网国产| 99久久人妻综合| 免费播放大片免费观看视频在线观看| 欧美日韩综合久久久久久| 欧美一级a爱片免费观看看| 国产精品国产三级国产专区5o| 久久久色成人| 日韩伦理黄色片| 美女脱内裤让男人舔精品视频| 晚上一个人看的免费电影| 亚洲综合精品二区| 日本爱情动作片www.在线观看| 一二三四中文在线观看免费高清| 免费观看无遮挡的男女| 亚洲aⅴ乱码一区二区在线播放| 亚洲高清免费不卡视频| 2021少妇久久久久久久久久久| 2022亚洲国产成人精品| 一区在线观看完整版| 色吧在线观看| 久久久精品94久久精品| 一级毛片久久久久久久久女| 日韩成人av中文字幕在线观看| 久久 成人 亚洲| 免费大片18禁| 亚洲欧美成人综合另类久久久| 亚洲欧美精品自产自拍| 自拍偷自拍亚洲精品老妇| 日韩不卡一区二区三区视频在线| 日日撸夜夜添| 亚洲内射少妇av| 久久久成人免费电影| 亚洲国产精品专区欧美| 日本午夜av视频| 亚洲欧美日韩无卡精品| 国产高清不卡午夜福利| 欧美最新免费一区二区三区| 精品一品国产午夜福利视频| 亚洲美女搞黄在线观看| 我要看黄色一级片免费的| 青春草亚洲视频在线观看| 日本欧美视频一区| 中文在线观看免费www的网站| 亚洲不卡免费看| 久久99热6这里只有精品| 美女内射精品一级片tv| 综合色丁香网| 国产精品国产三级专区第一集| 国产成人精品久久久久久| 国产伦理片在线播放av一区| 久久这里有精品视频免费| 免费少妇av软件| 99热这里只有精品一区| 在线观看一区二区三区激情| 大又大粗又爽又黄少妇毛片口| 国产精品不卡视频一区二区| 国产探花极品一区二区| 老熟女久久久| 国产片特级美女逼逼视频| 婷婷色综合www| 最近2019中文字幕mv第一页| 日韩亚洲欧美综合| 亚洲四区av| 51国产日韩欧美| 最近最新中文字幕大全电影3| 美女中出高潮动态图| 不卡视频在线观看欧美| 日本黄大片高清| 国产成人91sexporn| 深爱激情五月婷婷| 中文字幕精品免费在线观看视频 | 亚洲欧美一区二区三区国产| videos熟女内射| 国产白丝娇喘喷水9色精品| 免费人妻精品一区二区三区视频| 国产一区亚洲一区在线观看| 王馨瑶露胸无遮挡在线观看| 免费看av在线观看网站| 国产淫语在线视频| 午夜精品国产一区二区电影| 黄色视频在线播放观看不卡| 国产 一区 欧美 日韩| 国产成人午夜福利电影在线观看| 欧美三级亚洲精品| 精品国产乱码久久久久久小说| 日韩强制内射视频| 99久久精品国产国产毛片| 欧美成人精品欧美一级黄| 中文字幕久久专区| 亚洲欧美日韩东京热| 日韩av在线免费看完整版不卡| www.色视频.com| 综合色丁香网| 国产精品一区www在线观看| 一个人免费看片子| 尾随美女入室| 人人妻人人澡人人爽人人夜夜| 欧美+日韩+精品| 内射极品少妇av片p| 只有这里有精品99| 国产亚洲91精品色在线| 啦啦啦视频在线资源免费观看| 婷婷色av中文字幕| 亚洲av中文字字幕乱码综合| 免费大片黄手机在线观看| 1000部很黄的大片| 色婷婷av一区二区三区视频| 深爱激情五月婷婷| 一区二区三区免费毛片| 免费av中文字幕在线| 精品久久久久久久末码| 国产v大片淫在线免费观看| 精品一区二区三卡| 九草在线视频观看| 男女啪啪激烈高潮av片| 99热国产这里只有精品6| 美女xxoo啪啪120秒动态图| 日本av免费视频播放| 爱豆传媒免费全集在线观看| 亚洲婷婷狠狠爱综合网| 国产爽快片一区二区三区| 国产毛片在线视频| 日韩av不卡免费在线播放| 精品久久久久久久久av| av不卡在线播放| 午夜日本视频在线| 久久久精品94久久精品| 乱系列少妇在线播放| 日韩欧美一区视频在线观看 | 中文字幕人妻熟人妻熟丝袜美| 亚洲性久久影院| 高清午夜精品一区二区三区| 亚洲激情五月婷婷啪啪| 亚洲精品色激情综合| 99热网站在线观看| 成人亚洲精品一区在线观看 | 美女脱内裤让男人舔精品视频| 午夜福利在线在线| 18禁裸乳无遮挡动漫免费视频| 久久97久久精品| 久热久热在线精品观看| 熟女av电影| 日韩人妻高清精品专区| 国产亚洲91精品色在线| 亚洲欧美清纯卡通| 在线观看人妻少妇| 一级毛片黄色毛片免费观看视频| 自拍偷自拍亚洲精品老妇| 一级毛片电影观看| 久久久久精品性色| 国产在线男女| 亚洲激情五月婷婷啪啪| 国产男女内射视频| 麻豆成人午夜福利视频| 国产高清三级在线| 伦理电影大哥的女人| 国产精品偷伦视频观看了| 久久久久视频综合| av国产精品久久久久影院| 色婷婷久久久亚洲欧美| 新久久久久国产一级毛片| 日韩国内少妇激情av| 肉色欧美久久久久久久蜜桃| 亚洲伊人久久精品综合| 国产亚洲最大av| 午夜福利影视在线免费观看| 亚洲图色成人| 亚洲欧美成人精品一区二区| 寂寞人妻少妇视频99o| 日韩大片免费观看网站| 国产精品女同一区二区软件| 久久久久网色| 大又大粗又爽又黄少妇毛片口| 国产成人精品久久久久久| 各种免费的搞黄视频| 亚洲综合精品二区| 免费黄网站久久成人精品| 色5月婷婷丁香| 高清在线视频一区二区三区| 国产综合精华液| 在线免费十八禁| 丝袜喷水一区| 麻豆成人av视频| 国产极品天堂在线| 国产亚洲91精品色在线| 亚洲欧美中文字幕日韩二区| 欧美变态另类bdsm刘玥| 日本黄大片高清| 久久久国产一区二区| 天美传媒精品一区二区| 国产中年淑女户外野战色| 亚洲国产精品国产精品| 成人高潮视频无遮挡免费网站| 蜜臀久久99精品久久宅男| 久久人人爽人人爽人人片va| 99久久综合免费| 久久久久国产精品人妻一区二区| 天天躁日日操中文字幕| 欧美高清性xxxxhd video| 国产在线免费精品| 又大又黄又爽视频免费| 久久影院123| 丝袜脚勾引网站| 黄色欧美视频在线观看| 高清黄色对白视频在线免费看 | 中国三级夫妇交换| 天天躁日日操中文字幕| 婷婷色综合www| 久久久精品94久久精品| 日本猛色少妇xxxxx猛交久久| 毛片女人毛片| 久久精品熟女亚洲av麻豆精品| 永久网站在线| 七月丁香在线播放| 国内精品宾馆在线| 国产精品国产三级国产专区5o| h视频一区二区三区| 肉色欧美久久久久久久蜜桃| 久久精品国产自在天天线| 777米奇影视久久| 啦啦啦中文免费视频观看日本| 中文乱码字字幕精品一区二区三区| 国产中年淑女户外野战色| 亚洲第一av免费看| 高清视频免费观看一区二区| 免费看日本二区| 亚洲国产精品999| 汤姆久久久久久久影院中文字幕| 在线观看国产h片| 日韩强制内射视频| 中文字幕亚洲精品专区| 老司机影院毛片| 毛片女人毛片| 人人妻人人看人人澡| 美女主播在线视频| 一本久久精品| 久久久精品免费免费高清| 亚洲中文av在线| 99热网站在线观看| 国产 精品1| 秋霞伦理黄片| 久久国产乱子免费精品| 亚洲aⅴ乱码一区二区在线播放| 老师上课跳d突然被开到最大视频| 纵有疾风起免费观看全集完整版| 亚洲内射少妇av| 在线免费观看不下载黄p国产| tube8黄色片| 全区人妻精品视频| 99久久精品热视频| 王馨瑶露胸无遮挡在线观看| 欧美精品亚洲一区二区| 一级毛片久久久久久久久女| 欧美3d第一页| 日韩一区二区视频免费看| 午夜福利高清视频| 婷婷色综合大香蕉| 久久韩国三级中文字幕| 中文字幕免费在线视频6| 国产男女超爽视频在线观看| 国产男人的电影天堂91| 国产美女午夜福利| 九草在线视频观看| 一级片'在线观看视频| 麻豆精品久久久久久蜜桃| 免费av不卡在线播放| 国产亚洲欧美精品永久| 51国产日韩欧美| 日韩,欧美,国产一区二区三区| 汤姆久久久久久久影院中文字幕| 国产成人a∨麻豆精品| 国产色爽女视频免费观看| 日韩av在线免费看完整版不卡| 免费观看无遮挡的男女| 日日啪夜夜爽| 插逼视频在线观看| 有码 亚洲区| 久久精品国产亚洲网站| 人人妻人人看人人澡| 亚洲精品,欧美精品| 色视频在线一区二区三区| 欧美日韩亚洲高清精品| www.av在线官网国产| 哪个播放器可以免费观看大片| 91在线精品国自产拍蜜月| 亚洲av不卡在线观看| 成人高潮视频无遮挡免费网站| 九九久久精品国产亚洲av麻豆| 丰满少妇做爰视频| 免费观看性生交大片5| 国产高清国产精品国产三级 | a级一级毛片免费在线观看| 国产免费福利视频在线观看| 国产亚洲91精品色在线| 国产精品99久久99久久久不卡 | 在线观看免费视频网站a站| 亚洲精品国产av蜜桃| 国产成人一区二区在线| 1000部很黄的大片| 两个人的视频大全免费| 国产精品一区二区在线不卡| 欧美 日韩 精品 国产| 免费av不卡在线播放| 97超视频在线观看视频| 欧美日韩国产mv在线观看视频 | 国产爱豆传媒在线观看| av国产精品久久久久影院| 中文字幕人妻熟人妻熟丝袜美| 久久99精品国语久久久| 午夜老司机福利剧场| 国产成人午夜福利电影在线观看| 在线观看一区二区三区| videossex国产| 国产色婷婷99| 汤姆久久久久久久影院中文字幕| 18禁在线播放成人免费| 久热久热在线精品观看| 五月开心婷婷网| 97超视频在线观看视频| 久热久热在线精品观看| 亚洲av日韩在线播放| 日本爱情动作片www.在线观看| 在线观看人妻少妇| 国产精品一区二区性色av| 最近中文字幕高清免费大全6| 欧美一级a爱片免费观看看| 免费观看a级毛片全部| 老司机影院成人| 日韩成人伦理影院| 国产伦精品一区二区三区视频9| 我要看黄色一级片免费的| 欧美一级a爱片免费观看看| 亚洲精品久久久久久婷婷小说| 少妇裸体淫交视频免费看高清| 国产成人a∨麻豆精品| 国产男女内射视频| 免费黄色在线免费观看| av黄色大香蕉| 久久久久久九九精品二区国产| 一级毛片电影观看| 香蕉精品网在线| 色网站视频免费| 精品久久国产蜜桃| 日韩一区二区三区影片| 好男人视频免费观看在线| 一级毛片久久久久久久久女| 国产精品久久久久久久久免| 欧美日韩一区二区视频在线观看视频在线| 久久综合国产亚洲精品| 国产淫语在线视频| av国产免费在线观看| 日韩免费高清中文字幕av| 国产精品国产av在线观看| 99热网站在线观看| 99热这里只有是精品在线观看| 成人亚洲欧美一区二区av| 日韩三级伦理在线观看| 亚洲成人一二三区av| 亚洲经典国产精华液单| 久久久精品94久久精品| 最近最新中文字幕大全电影3| 街头女战士在线观看网站| 亚洲精品国产av成人精品| 不卡视频在线观看欧美| 免费观看性生交大片5| 丝瓜视频免费看黄片| 久久久久久久久久人人人人人人| 国产精品一二三区在线看| 欧美人与善性xxx| 日本vs欧美在线观看视频 | 久久韩国三级中文字幕| 欧美97在线视频| 久久99精品国语久久久| 99久久精品国产国产毛片| 亚洲,欧美,日韩| 高清av免费在线| 国产91av在线免费观看| 1000部很黄的大片| 天美传媒精品一区二区| 日韩欧美一区视频在线观看 | 国产真实伦视频高清在线观看| 中文字幕精品免费在线观看视频 | 91狼人影院| 在线观看三级黄色| 国产久久久一区二区三区| 欧美亚洲 丝袜 人妻 在线| 午夜福利视频精品| 一级毛片aaaaaa免费看小| 亚洲av中文av极速乱| 午夜视频国产福利| 国产 一区精品| 亚洲精品国产av成人精品| 久久久久久久久久成人| 大片电影免费在线观看免费| 国产精品久久久久久精品古装| 99精国产麻豆久久婷婷| 欧美国产精品一级二级三级 | 国产精品一区二区在线不卡| 国产欧美日韩一区二区三区在线 | 晚上一个人看的免费电影| 最近中文字幕2019免费版| 有码 亚洲区| 99热6这里只有精品| 伦精品一区二区三区| 大片电影免费在线观看免费| 亚洲不卡免费看| 国产亚洲5aaaaa淫片| 亚洲欧美一区二区三区国产| 黄色视频在线播放观看不卡| 久久av网站| 99热网站在线观看| 亚洲四区av| 国产精品一区二区三区四区免费观看| 丰满人妻一区二区三区视频av| 久久久久性生活片| 你懂的网址亚洲精品在线观看| 我要看日韩黄色一级片| 欧美xxxx性猛交bbbb| 欧美最新免费一区二区三区| 国产高清国产精品国产三级 | 久久热精品热| 毛片一级片免费看久久久久| 80岁老熟妇乱子伦牲交| 国产伦精品一区二区三区视频9| 久久久久性生活片| 精品人妻偷拍中文字幕| 国内精品宾馆在线| 一级毛片黄色毛片免费观看视频| 91精品伊人久久大香线蕉| 亚洲精华国产精华液的使用体验| 日韩三级伦理在线观看| 欧美成人精品欧美一级黄| av网站免费在线观看视频| 一级毛片久久久久久久久女| 久久久久久人妻| 国产精品.久久久| 亚洲精品色激情综合| 一级黄片播放器| 亚洲中文av在线| 免费久久久久久久精品成人欧美视频 | 日本黄色日本黄色录像| 欧美日韩精品成人综合77777| 成人国产av品久久久| 精品国产三级普通话版| 舔av片在线| 丝袜脚勾引网站| 亚洲精品aⅴ在线观看| 中国国产av一级| 国产精品久久久久久精品古装| 国产成人精品福利久久| 一本—道久久a久久精品蜜桃钙片| 成年人午夜在线观看视频| 国产午夜精品一二区理论片| 欧美精品一区二区免费开放| 天美传媒精品一区二区| 香蕉精品网在线| 久久久久久久久久人人人人人人| 国产精品国产三级国产av玫瑰| 国内精品宾馆在线| 久久久欧美国产精品| 一区二区三区精品91| 欧美zozozo另类| 国产成人精品婷婷| 免费看av在线观看网站| 久久99蜜桃精品久久| 日本黄色片子视频| 亚洲精品久久久久久婷婷小说| 亚洲欧美日韩另类电影网站 | 大话2 男鬼变身卡| 久久国内精品自在自线图片| 最近中文字幕2019免费版| 毛片一级片免费看久久久久| 看免费成人av毛片| 免费大片18禁| 欧美激情国产日韩精品一区| 我要看日韩黄色一级片| kizo精华| 国产成人91sexporn| 免费黄网站久久成人精品| 精品久久久久久久久亚洲| 国产一级毛片在线| 中文欧美无线码| 国产亚洲一区二区精品| 亚洲av中文av极速乱| videossex国产| 欧美最新免费一区二区三区| 久久国产亚洲av麻豆专区| 亚洲中文av在线| 91精品国产国语对白视频| 大香蕉久久网| 久久青草综合色| 国产成人freesex在线| 我要看日韩黄色一级片| 久热这里只有精品99| 超碰av人人做人人爽久久| 亚洲精华国产精华液的使用体验| 狂野欧美激情性xxxx在线观看| 在线观看免费视频网站a站| 一边亲一边摸免费视频| 久久婷婷青草| 亚洲内射少妇av| 成年女人在线观看亚洲视频| 大码成人一级视频| 久久精品国产亚洲av天美| 久久国产乱子免费精品| 妹子高潮喷水视频| 国产 一区精品| 免费看av在线观看网站| 午夜福利高清视频| 99久久精品国产国产毛片| 国产免费一级a男人的天堂| 色婷婷久久久亚洲欧美| 国产成人精品久久久久久| 丰满少妇做爰视频| 亚洲av二区三区四区| 啦啦啦视频在线资源免费观看| 99re6热这里在线精品视频| 成人国产av品久久久| 一区二区三区免费毛片| 小蜜桃在线观看免费完整版高清| 中文精品一卡2卡3卡4更新| 少妇精品久久久久久久| 日本vs欧美在线观看视频 | 亚洲人成网站在线播| 2022亚洲国产成人精品| 亚洲国产精品成人久久小说| 制服丝袜香蕉在线| 日韩亚洲欧美综合| 欧美日韩一区二区视频在线观看视频在线| 嫩草影院入口| 欧美日韩综合久久久久久| 男男h啪啪无遮挡| 黄色配什么色好看| 少妇裸体淫交视频免费看高清| 人人妻人人看人人澡| 精品久久久久久久久av| av专区在线播放| 美女国产视频在线观看| 蜜桃久久精品国产亚洲av| 少妇被粗大猛烈的视频| 妹子高潮喷水视频| 亚洲熟女精品中文字幕| 久久精品久久久久久久性| 亚洲精品视频女| 国产精品久久久久久精品电影小说 | 欧美 日韩 精品 国产| 精品一区二区三卡| 777米奇影视久久| 免费人成在线观看视频色| 国产成人精品婷婷| 夜夜爽夜夜爽视频| 色哟哟·www| 国产av一区二区精品久久 | 久久99热这里只频精品6学生| 在线精品无人区一区二区三 | 男女无遮挡免费网站观看| 国产亚洲欧美精品永久| 亚洲三级黄色毛片| 夫妻午夜视频| 国产精品不卡视频一区二区| 婷婷色麻豆天堂久久| 视频区图区小说| 精品99又大又爽又粗少妇毛片| 26uuu在线亚洲综合色| 蜜桃久久精品国产亚洲av| 一级毛片 在线播放| 男人添女人高潮全过程视频| 国产精品欧美亚洲77777| 尾随美女入室| 亚洲色图综合在线观看| 亚洲成色77777| 欧美日韩国产mv在线观看视频 | 男女下面进入的视频免费午夜| 亚洲国产成人一精品久久久| 最黄视频免费看| 国产久久久一区二区三区| 精品少妇黑人巨大在线播放| 亚洲欧美一区二区三区黑人 | 国产高清不卡午夜福利| 欧美xxxx黑人xx丫x性爽| 九九在线视频观看精品| 2021少妇久久久久久久久久久| 亚洲电影在线观看av| 国产亚洲精品久久久com| 亚洲精品aⅴ在线观看| 一级黄片播放器| 亚洲美女黄色视频免费看| 婷婷色av中文字幕| 欧美日韩精品成人综合77777| 欧美极品一区二区三区四区| 搡女人真爽免费视频火全软件| 亚洲精品456在线播放app| 国产精品成人在线| 秋霞在线观看毛片| 亚洲精品视频女| 秋霞在线观看毛片| 免费不卡的大黄色大毛片视频在线观看| 国产亚洲av片在线观看秒播厂| 亚州av有码| 欧美变态另类bdsm刘玥| 美女xxoo啪啪120秒动态图| 亚洲不卡免费看| 国产成人精品久久久久久| 国产人妻一区二区三区在| 午夜福利在线在线| 中文字幕免费在线视频6| 丰满人妻一区二区三区视频av| 久久这里有精品视频免费| 国产精品免费大片| 伦精品一区二区三区| 毛片一级片免费看久久久久| 欧美日韩一区二区视频在线观看视频在线| 国产人妻一区二区三区在| 九草在线视频观看| 久久久久久久久大av|