• <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法對圖們江干流含沙量年際變化的分析
    av.在线天堂| 国产精品三级大全| 免费在线观看黄色视频的| 日韩一本色道免费dvd| 欧美激情高清一区二区三区 | kizo精华| 久久久精品免费免费高清| 91成人精品电影| 日本欧美视频一区| 一区二区日韩欧美中文字幕| 成年av动漫网址| 亚洲国产成人一精品久久久| 久久久久久久精品精品| 国产一区二区在线观看av| 老鸭窝网址在线观看| 日韩人妻精品一区2区三区| 精品福利永久在线观看| 午夜日本视频在线| 亚洲精品日韩在线中文字幕| 美女脱内裤让男人舔精品视频| 深夜精品福利| 看十八女毛片水多多多| 国产色婷婷99| 久久久久久久国产电影| 99久国产av精品国产电影| 久久影院123| 爱豆传媒免费全集在线观看| 男的添女的下面高潮视频| 黄色配什么色好看| 日韩精品免费视频一区二区三区| 日日啪夜夜爽| 一级黄片播放器| 久久韩国三级中文字幕| 久久青草综合色| 中文字幕av电影在线播放| 久久午夜综合久久蜜桃| av一本久久久久| 校园人妻丝袜中文字幕| 99久久综合免费| 99久久精品国产国产毛片| 男人爽女人下面视频在线观看| 免费黄色在线免费观看| 黄色 视频免费看| 国产在线一区二区三区精| 免费黄频网站在线观看国产| 一区福利在线观看| 亚洲,欧美精品.| 精品少妇黑人巨大在线播放| 国产欧美亚洲国产| 国产av精品麻豆| 大话2 男鬼变身卡| 纵有疾风起免费观看全集完整版| 国产av码专区亚洲av| 欧美xxⅹ黑人| 久久久亚洲精品成人影院| 久久精品亚洲av国产电影网| 成人黄色视频免费在线看| 性少妇av在线| 国产精品 欧美亚洲| 女人精品久久久久毛片| 久久久久视频综合| 久久久精品94久久精品| 青青草视频在线视频观看| 久久午夜福利片| 下体分泌物呈黄色| 婷婷成人精品国产| 亚洲美女搞黄在线观看| 国产成人精品一,二区| 久久精品aⅴ一区二区三区四区 | 亚洲国产精品一区二区三区在线| 国产精品偷伦视频观看了| 亚洲天堂av无毛| 国产国语露脸激情在线看| 日韩三级伦理在线观看| 一级毛片电影观看| 久久99热这里只频精品6学生| 中文乱码字字幕精品一区二区三区| 一边摸一边做爽爽视频免费| 久久久久久人人人人人| 女的被弄到高潮叫床怎么办| 国产日韩欧美在线精品| 两个人免费观看高清视频| 99国产精品免费福利视频| 看免费成人av毛片| 天堂中文最新版在线下载| 日产精品乱码卡一卡2卡三| 久久精品久久久久久噜噜老黄| 爱豆传媒免费全集在线观看| 中国国产av一级| 久久精品国产综合久久久| 久久精品夜色国产| 免费黄网站久久成人精品| 精品国产一区二区三区四区第35| 尾随美女入室| 九草在线视频观看| 最近最新中文字幕大全免费视频 | 999精品在线视频| 高清在线视频一区二区三区| 日韩电影二区| 久久久久国产网址| 制服丝袜香蕉在线| 免费观看a级毛片全部| 国产欧美日韩综合在线一区二区| 欧美日韩亚洲国产一区二区在线观看 | 亚洲国产成人一精品久久久| 亚洲国产毛片av蜜桃av| 18+在线观看网站| 看免费成人av毛片| 大话2 男鬼变身卡| 亚洲色图 男人天堂 中文字幕| 2018国产大陆天天弄谢| 国产成人免费无遮挡视频| 亚洲一区中文字幕在线| 国产成人欧美| 国产日韩欧美在线精品| 国产成人精品福利久久| 精品人妻偷拍中文字幕| 成人国产麻豆网| 免费av中文字幕在线| videossex国产| 日本91视频免费播放| 成年动漫av网址| 免费大片黄手机在线观看| 中文字幕另类日韩欧美亚洲嫩草| av在线观看视频网站免费| 女人精品久久久久毛片| 日本av手机在线免费观看| 国产精品熟女久久久久浪| 三上悠亚av全集在线观看| 久久婷婷青草| 九草在线视频观看| 精品人妻偷拍中文字幕| 菩萨蛮人人尽说江南好唐韦庄| 亚洲精品国产色婷婷电影| 国产精品久久久久成人av| 亚洲成人一二三区av| 免费黄频网站在线观看国产| 久久久久精品人妻al黑| 精品福利永久在线观看| 亚洲av成人精品一二三区| 成人18禁高潮啪啪吃奶动态图| 国产精品久久久久成人av| www.自偷自拍.com| 尾随美女入室| 你懂的网址亚洲精品在线观看| 这个男人来自地球电影免费观看 | 精品久久久久久电影网| 欧美精品av麻豆av| 久久精品国产综合久久久| 久久久精品区二区三区| 精品一区二区三卡| 欧美日本中文国产一区发布| 美女大奶头黄色视频| 色婷婷av一区二区三区视频| 国产精品 欧美亚洲| 亚洲色图综合在线观看| a 毛片基地| 极品少妇高潮喷水抽搐| 熟女av电影| 国产免费福利视频在线观看| 男人操女人黄网站| 国产人伦9x9x在线观看 | 看十八女毛片水多多多| 国产精品偷伦视频观看了| 亚洲国产欧美在线一区| 国精品久久久久久国模美| 久久久久国产精品人妻一区二区| 免费少妇av软件| 久久女婷五月综合色啪小说| 久久99蜜桃精品久久| 中文字幕色久视频| 色视频在线一区二区三区| 日韩 亚洲 欧美在线| 精品人妻一区二区三区麻豆| 国产精品二区激情视频| 男女下面插进去视频免费观看| 一区二区av电影网| 成人影院久久| 欧美人与性动交α欧美软件| 国产无遮挡羞羞视频在线观看| 国产又爽黄色视频| 伦理电影大哥的女人| 街头女战士在线观看网站| av一本久久久久| 欧美精品亚洲一区二区| 国产色婷婷99| 久久久国产精品麻豆| 国产极品天堂在线| 精品午夜福利在线看| 久久97久久精品| 国产精品国产三级国产专区5o| 免费少妇av软件| 欧美最新免费一区二区三区| 看非洲黑人一级黄片| www.熟女人妻精品国产| 我要看黄色一级片免费的| 十八禁高潮呻吟视频| 国产一区二区三区综合在线观看| 91国产中文字幕| 搡老乐熟女国产| 男人爽女人下面视频在线观看| 精品亚洲成a人片在线观看| 国产高清不卡午夜福利| 亚洲色图 男人天堂 中文字幕| 精品久久蜜臀av无| 一区二区三区四区激情视频| 欧美日韩精品网址| 日本91视频免费播放| 大香蕉久久网| 婷婷色麻豆天堂久久| 成人国产av品久久久| 日韩在线高清观看一区二区三区| 熟女电影av网| 亚洲激情五月婷婷啪啪| 久久久国产一区二区| 午夜福利影视在线免费观看| 国产一级毛片在线| 午夜福利一区二区在线看| 久久久久久久久久人人人人人人| 少妇被粗大猛烈的视频| 久久午夜综合久久蜜桃| 欧美日韩精品网址| 亚洲欧美色中文字幕在线| 久久国产精品大桥未久av| 在线观看免费日韩欧美大片| 欧美精品av麻豆av| 成年人午夜在线观看视频| 久久久久久人妻| 亚洲成av片中文字幕在线观看 | 啦啦啦在线观看免费高清www| 久久国产亚洲av麻豆专区| 99热全是精品| 亚洲三级黄色毛片| 丝瓜视频免费看黄片| 在线观看人妻少妇| 国产一区二区激情短视频 | 最黄视频免费看| 国产一区二区 视频在线| 色婷婷久久久亚洲欧美| 纯流量卡能插随身wifi吗| 亚洲久久久国产精品| 久久99精品国语久久久| 熟女少妇亚洲综合色aaa.| 亚洲婷婷狠狠爱综合网| 国产无遮挡羞羞视频在线观看| 少妇的逼水好多| 香蕉丝袜av| 国产麻豆69| 啦啦啦在线免费观看视频4| 青春草国产在线视频| 国产午夜精品一二区理论片| 天天操日日干夜夜撸| 亚洲精品国产一区二区精华液| av国产久精品久网站免费入址| 在线亚洲精品国产二区图片欧美| kizo精华| 免费高清在线观看视频在线观看| 一边亲一边摸免费视频| 婷婷色综合www| 欧美av亚洲av综合av国产av | 最近中文字幕2019免费版| 国产成人91sexporn| 亚洲,一卡二卡三卡| 狠狠精品人妻久久久久久综合| 亚洲精品美女久久久久99蜜臀 | 女的被弄到高潮叫床怎么办| 欧美变态另类bdsm刘玥| 电影成人av| 国产精品三级大全| 黄色视频在线播放观看不卡| 丝袜喷水一区| 欧美激情高清一区二区三区 | 久久这里只有精品19| 啦啦啦视频在线资源免费观看| 久热久热在线精品观看| 亚洲成av片中文字幕在线观看 | 成人18禁高潮啪啪吃奶动态图| 99久久精品国产国产毛片| 成年av动漫网址| 色网站视频免费| 中国三级夫妇交换| 不卡av一区二区三区| 日韩精品免费视频一区二区三区| 国产福利在线免费观看视频| av网站在线播放免费| 黄色视频在线播放观看不卡| 两个人免费观看高清视频| 欧美精品一区二区大全| 午夜久久久在线观看| 在现免费观看毛片| h视频一区二区三区| 亚洲欧美精品自产自拍| 国产精品三级大全| 亚洲欧美精品综合一区二区三区 | 一区二区三区精品91| √禁漫天堂资源中文www| 侵犯人妻中文字幕一二三四区| 赤兔流量卡办理| 亚洲美女黄色视频免费看| 观看美女的网站| 日韩制服丝袜自拍偷拍| 亚洲中文av在线| 国产97色在线日韩免费| 国产精品久久久久成人av| 天天影视国产精品| 女人久久www免费人成看片| 欧美日韩一级在线毛片| 一二三四中文在线观看免费高清| 蜜桃国产av成人99| 成年人免费黄色播放视频| 亚洲熟女精品中文字幕| 91在线精品国自产拍蜜月| 亚洲av综合色区一区| 国产一区二区三区综合在线观看| 亚洲天堂av无毛| 亚洲成人手机| 巨乳人妻的诱惑在线观看| 永久网站在线| av一本久久久久| 妹子高潮喷水视频| 午夜福利一区二区在线看| 青草久久国产| 好男人视频免费观看在线| 日本欧美视频一区| 1024视频免费在线观看| 欧美另类一区| 免费黄色在线免费观看| 久久久欧美国产精品| 观看av在线不卡| 亚洲婷婷狠狠爱综合网| 在线观看免费视频网站a站| 亚洲一区二区三区欧美精品| 青青草视频在线视频观看| 性色av一级| 久久精品人人爽人人爽视色| 精品久久久久久电影网| av国产精品久久久久影院| 国产精品久久久av美女十八| 一级片免费观看大全| 国产视频首页在线观看| 少妇熟女欧美另类| 大香蕉久久网| 久久久久久伊人网av| 99久国产av精品国产电影| tube8黄色片| 久久久久久久久久久久大奶| 男女国产视频网站| 天堂俺去俺来也www色官网| 久久精品人人爽人人爽视色| 欧美精品国产亚洲| 建设人人有责人人尽责人人享有的| 国产在线视频一区二区| 亚洲av.av天堂| 一二三四在线观看免费中文在| kizo精华| 最新的欧美精品一区二区| 青春草视频在线免费观看| 99re6热这里在线精品视频| 亚洲熟女精品中文字幕| 久久久久久久久免费视频了| 国产精品国产三级国产专区5o| 在线天堂中文资源库| 哪个播放器可以免费观看大片| 少妇猛男粗大的猛烈进出视频| 999久久久国产精品视频| 高清黄色对白视频在线免费看| 夫妻午夜视频| 久久久久久人人人人人| 午夜福利在线观看免费完整高清在| 色哟哟·www| 中文字幕亚洲精品专区| 免费观看av网站的网址| 国产亚洲午夜精品一区二区久久| 国产亚洲一区二区精品| 久久久久国产网址| 一本—道久久a久久精品蜜桃钙片| 不卡av一区二区三区| 人妻 亚洲 视频| 国产精品一国产av| 亚洲国产欧美网| 国产精品久久久av美女十八| 最新的欧美精品一区二区| 美女xxoo啪啪120秒动态图| 亚洲国产欧美日韩在线播放| 90打野战视频偷拍视频| 亚洲一级一片aⅴ在线观看| 国语对白做爰xxxⅹ性视频网站| 欧美变态另类bdsm刘玥| 青草久久国产| 一级毛片我不卡| 国产淫语在线视频| 一区二区三区精品91| 午夜福利视频精品| 丝瓜视频免费看黄片| 精品人妻偷拍中文字幕| 看免费av毛片| 久久人人97超碰香蕉20202| 欧美精品一区二区大全| 黑人欧美特级aaaaaa片| 欧美 日韩 精品 国产| 亚洲精品美女久久久久99蜜臀 | 999久久久国产精品视频| 国产一区二区激情短视频 | av网站在线播放免费| 一级黄片播放器| 久久鲁丝午夜福利片| 爱豆传媒免费全集在线观看| 免费黄色在线免费观看| 免费在线观看完整版高清| 美女主播在线视频| 国产精品熟女久久久久浪| 亚洲欧美成人精品一区二区| 又大又黄又爽视频免费| 亚洲成人手机| 欧美激情 高清一区二区三区| 亚洲av成人精品一二三区| 中文乱码字字幕精品一区二区三区| 国产精品国产av在线观看| 26uuu在线亚洲综合色| 欧美人与性动交α欧美软件| 久久毛片免费看一区二区三区| 两性夫妻黄色片| 一二三四中文在线观看免费高清| 天天影视国产精品| 国产 一区精品| xxxhd国产人妻xxx| 欧美在线黄色| 免费高清在线观看视频在线观看| 亚洲色图综合在线观看| 波多野结衣av一区二区av| 欧美日韩成人在线一区二区| 中文字幕av电影在线播放| 日韩制服丝袜自拍偷拍| 午夜福利网站1000一区二区三区| 久久青草综合色| 少妇精品久久久久久久| 边亲边吃奶的免费视频| 我的亚洲天堂| 日韩视频在线欧美| 免费在线观看黄色视频的| 一级a爱视频在线免费观看| 国产欧美亚洲国产| 国产亚洲午夜精品一区二区久久| 18禁动态无遮挡网站| 边亲边吃奶的免费视频| 午夜福利影视在线免费观看| 成年女人在线观看亚洲视频| 免费播放大片免费观看视频在线观看| 啦啦啦啦在线视频资源| 欧美xxⅹ黑人| 黄色毛片三级朝国网站| 青春草视频在线免费观看| 国产麻豆69| 久热久热在线精品观看| 日韩av不卡免费在线播放| 久久午夜福利片| 岛国毛片在线播放| 国产高清国产精品国产三级| tube8黄色片| 亚洲av成人精品一二三区| 国产精品久久久久久精品古装| 欧美黄色片欧美黄色片| 精品久久久久久电影网| 日韩中文字幕欧美一区二区 | 国产精品免费视频内射| 欧美日韩综合久久久久久| 国产成人免费无遮挡视频| 国产欧美亚洲国产| 国产综合精华液| 亚洲精华国产精华液的使用体验| 中文字幕色久视频| 日韩制服丝袜自拍偷拍| 中文欧美无线码| 一二三四在线观看免费中文在| 国产免费福利视频在线观看| 国产免费一区二区三区四区乱码| 欧美+日韩+精品| 一级片'在线观看视频| 黑人巨大精品欧美一区二区蜜桃| a级毛片在线看网站| 亚洲欧美清纯卡通| 久久久国产精品麻豆| 久久热在线av| 久久精品aⅴ一区二区三区四区 | 国产精品久久久久久av不卡| 精品少妇黑人巨大在线播放| 午夜福利,免费看| 亚洲精品乱久久久久久| 久久精品国产亚洲av天美| 波多野结衣一区麻豆| 国产成人aa在线观看| 日本免费在线观看一区| 日韩成人av中文字幕在线观看| 国产精品国产av在线观看| 国产男女内射视频| 免费黄频网站在线观看国产| 色婷婷av一区二区三区视频| 晚上一个人看的免费电影| 亚洲av国产av综合av卡| 婷婷色综合www| www.av在线官网国产| 天堂8中文在线网| 80岁老熟妇乱子伦牲交| 久久精品久久久久久噜噜老黄| 青春草国产在线视频| 免费高清在线观看视频在线观看| 黄网站色视频无遮挡免费观看| 国产精品国产三级专区第一集| 日韩 亚洲 欧美在线| 一级,二级,三级黄色视频| 成年人午夜在线观看视频| 肉色欧美久久久久久久蜜桃| 午夜福利影视在线免费观看| 国产精品av久久久久免费| 午夜福利视频精品| 老汉色av国产亚洲站长工具| 久久综合国产亚洲精品| 久久毛片免费看一区二区三区| 国产精品三级大全| 人妻 亚洲 视频| 国产片内射在线| 赤兔流量卡办理| 亚洲欧洲精品一区二区精品久久久 | 国产精品一国产av| 天堂俺去俺来也www色官网| 丝袜人妻中文字幕| 一边摸一边做爽爽视频免费| 亚洲av国产av综合av卡| 久久久久国产精品人妻一区二区| 亚洲激情五月婷婷啪啪| 日韩 亚洲 欧美在线| 亚洲av日韩在线播放| 不卡av一区二区三区| 亚洲国产精品一区二区三区在线| 欧美 亚洲 国产 日韩一| 久久99一区二区三区| www日本在线高清视频| 男女下面插进去视频免费观看| 欧美日韩精品网址| 人成视频在线观看免费观看| 久久久久精品性色| 欧美xxⅹ黑人| 99久国产av精品国产电影| 午夜福利在线观看免费完整高清在| 中文字幕精品免费在线观看视频| 午夜福利在线观看免费完整高清在| 亚洲精品国产av成人精品| 久久国产精品大桥未久av| 99九九在线精品视频| 国产av码专区亚洲av| 人妻人人澡人人爽人人| 大陆偷拍与自拍| 中文字幕精品免费在线观看视频| 天堂中文最新版在线下载| 老司机影院成人| 中文天堂在线官网| 亚洲精品第二区| 青草久久国产| 久久久久久人人人人人| 日韩伦理黄色片| av在线app专区| 亚洲av.av天堂| 天天躁狠狠躁夜夜躁狠狠躁| 久久久久国产一级毛片高清牌| 久久狼人影院| 国产亚洲最大av| 久久亚洲国产成人精品v| 制服人妻中文乱码| 91精品三级在线观看| 在线观看免费高清a一片| 免费高清在线观看日韩| 亚洲欧美清纯卡通| 永久免费av网站大全| 97精品久久久久久久久久精品| 国产黄频视频在线观看| 女的被弄到高潮叫床怎么办| 国产熟女午夜一区二区三区| 国产爽快片一区二区三区| 国产成人91sexporn| 免费少妇av软件| 国产精品 欧美亚洲| 亚洲欧美成人精品一区二区| 人人妻人人添人人爽欧美一区卜| 亚洲精品美女久久久久99蜜臀 | 精品国产一区二区三区久久久樱花| 又黄又粗又硬又大视频| 高清黄色对白视频在线免费看| 日韩 亚洲 欧美在线| 久久免费观看电影| 麻豆乱淫一区二区| 在线观看免费日韩欧美大片| 美女视频免费永久观看网站| 精品国产一区二区久久| 久久久久久久亚洲中文字幕| 亚洲精品国产一区二区精华液| 免费看不卡的av| 91在线精品国自产拍蜜月| 波野结衣二区三区在线| 久久国内精品自在自线图片| 一级片免费观看大全| 人妻系列 视频| 少妇被粗大猛烈的视频| av线在线观看网站| 日本猛色少妇xxxxx猛交久久| 国产午夜精品一二区理论片| 午夜免费鲁丝| 欧美最新免费一区二区三区| av不卡在线播放| 国产黄色免费在线视频| 国产精品久久久久久精品古装| 国产国语露脸激情在线看| 亚洲国产看品久久| 免费av中文字幕在线| 不卡视频在线观看欧美| 午夜福利影视在线免费观看|