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

    1956 年~2020 年黃河源區(qū)徑流變化規(guī)律分析

    2022-07-08 15:45:56藍云龍關(guān)銅壘徐東坡劉佳嘉
    陜西水利 2022年6期
    關(guān)鍵詞:趨勢性瑪曲源區(qū)

    藍云龍,黎 曙,李 霞,關(guān)銅壘,徐東坡,劉佳嘉

    (1.黃河水利委員會西寧水文水資源勘測局,青海 西寧 810008;2.中國水利水電科學(xué)研究院,北京100038;3. 黃河勘測規(guī)劃設(shè)計研究院有限公司,河南 鄭州 450003)

    1 引言

    隨著全球氣候變暖,氣候變異性越來越突出,影響水文循環(huán)過程,使降水、蒸散發(fā)、徑流、冰川和凍土等發(fā)生變化,對水資源管理、規(guī)劃及配置提出了嚴(yán)峻的考驗。由于氣候變化以及人類活動頻繁,黃河流域徑流量減少趨勢甚是明顯[1]。黃河源區(qū)(唐乃亥水文站以上區(qū)域)是我國重點保護的水源涵養(yǎng)區(qū),一方面,氣候變化導(dǎo)致源區(qū)內(nèi)的凍土退化和冰川消融;另一方面,黃河源區(qū)人類活動更加頻繁,使得源區(qū)植被、草原、濕地等發(fā)生變化,進而使得源區(qū)徑流大小變化明顯,徑流變化作為水源涵養(yǎng)量的一個影響因素,對水源涵養(yǎng)量的影響不容忽視[2]。因此,研究黃河源區(qū)徑流的變化規(guī)律對黃河流域生態(tài)保護和高質(zhì)量發(fā)展、增加源區(qū)水源涵養(yǎng)量具有重要的參考價值。

    目前關(guān)于徑流變化特征的研究主要利用實測數(shù)據(jù)分析其趨勢性、突變性及周期性。秦年秀等人利用長江流域代表性水文站宜昌站、漢口站和大通站,分別選取1882 年~2000 年、1870 年~2000 年和1950 年~2000 年的月平均流量數(shù)據(jù)系列,通過非參數(shù)Mann-Kendall 法來檢驗長江徑流的趨勢變化[3]。劉嘉琦等人選取1950 年~2009 年長江流域大通水文站逐日流量系列,針對年、枯季、洪季、汛后及逐月平均流量序列,通過Pettitt 突變點檢驗法、Mann-Kendall 趨勢性檢驗法進行突變點分析和分段及整體序列的演變規(guī)律分析[4]。張嵐婷等人利用大渡河流域內(nèi)年徑流量數(shù)據(jù),為研究1951 年~2012 年序列的變化趨勢、周期性及突變點,選取Mann-Kendall 檢驗法、小波分析法、有序聚類法、累積距平法及Pettitt 突變檢驗法等方法對序列進行處理分析[5]。

    黃河源區(qū)徑流變化規(guī)律對研究黃河源區(qū)水文循環(huán)特征及水文要素的預(yù)測具有重要意義。鄭紅星等人為研究黃河源區(qū)瑪曲站、唐乃亥站的月天然徑流序列演變規(guī)律,選取1952 年~1997 年數(shù)據(jù)分析得出瑪曲站年內(nèi)分配的不均勻性、集中度以及相對變化幅度均略高于唐乃亥站,而其絕對變化幅度較小[6]。劉希勝等人為黃河源徑流變化規(guī)律,選取過程線法、集中度和集中期方法分析了黃河沿站、吉邁站、瑪曲站及唐乃亥站4個水文站1960 年~2012 年徑流的集中度和趨勢性,研究發(fā)現(xiàn)吉邁以上徑流量年際變化呈不顯著增加,而吉邁以下呈減少趨勢[7]。蔡宜晴等人為研究三江源區(qū)1956 年~2012 年的年徑流量演變規(guī)律,采用線性回歸、滑動平均、Mann-Kendall 檢驗、累積距平、滑動t 檢驗和Morlet 小波分析等方法研究發(fā)現(xiàn)黃河源區(qū)吉邁站呈上升趨勢,唐乃亥站呈下降趨勢,吉邁站在2004 年發(fā)生顯著性均值突變[8]。周帥等人利用綜合Mann-Kendall 方法和累積距平方法診斷出黃河源區(qū)1960 年~2010 年徑流量呈減少趨勢,徑流突變點位置在1989 年[9]。趙延存等人利用黃河源區(qū)1956 年~2019 年徑流序列,采用啟發(fā)式分割法、Mann-Kendall 檢驗法和集合經(jīng)驗?zāi)B(tài)分解(EEMD)方法對源區(qū)年徑流序列分析得到徑流序列在1990 年和2017 年發(fā)生變異,呈現(xiàn)不顯著下降趨勢,主周期為2.56 年[10]。綜上所述,各學(xué)者分析黃河源區(qū)選取的徑流數(shù)據(jù)系列及水文站點各不相同,多選取黃河沿、吉邁、瑪曲及唐乃亥干流水文站,時間序列多在2012 年以前。因此,本文仍選取黃河沿站、瑪曲站、唐乃亥站作為對象,水文系列選取1956 年~2020 年,采用Mann-Kendall 趨勢性分析法、Pettitt 突變分析法和小波分析法對黃河源區(qū)各站點年徑流量的趨勢性、突變性及周期性進行分析。

    2 研究區(qū)概況和資料

    黃河源區(qū)為唐乃亥水文站以上的匯水區(qū)域,面積12.2 萬km2,占黃河流域總面積的17%,是黃河的主要產(chǎn)流區(qū)[2]。黃河源區(qū)位于青藏高原東北部,見圖1,涉及青海,甘肅,四川三省,源區(qū)內(nèi)地勢總體呈西高東低[11]。黃河源區(qū)多年(1960 年~2019 年)平均年降水量為564.9mm,多集中在6 月~9 月[12]。黃河沿站、瑪曲站、唐乃亥站3 個水文站1956 年~2020 年徑流數(shù)據(jù)來源于水文統(tǒng)計年鑒《黃河流域水文資料第一冊(黃河上游區(qū)上段)》,其中黃河沿站1968 年~1975 年數(shù)據(jù)缺失,瑪曲站1956 年~1958 年及1989 年數(shù)據(jù)缺失,根據(jù)相鄰站點相關(guān)性進行差值。

    圖1 黃河源區(qū)分區(qū)圖

    3 研究方法

    3.1 Mann-Kendall 趨勢分析法

    Mann-Kendall 趨勢性分析方法(簡稱MK 檢驗)是一種非參數(shù)檢驗,在時間序列趨勢分析中該方法是世界氣象組織推薦并已廣泛使用的非參數(shù)檢驗方法,在水文分析中各學(xué)者利用該方法分析降水、徑流、氣溫及水質(zhì)等要素時間序列的趨勢變化[13]。在MK 檢驗中,原假設(shè)H0表示數(shù)據(jù)集X的樣本獨立同分布,沒有趨勢性存在,備擇假設(shè)H1表示數(shù)據(jù)集X中存在一個單調(diào)的趨勢性變化[14]。

    MK 檢驗所構(gòu)造的統(tǒng)計量為:

    其中,

    式中:xk和xi為樣本數(shù)據(jù)值,n為樣本容量,根據(jù)xk-xi的正負,sgn(xk-xi)分別為1、0 和-1。如果-Z1-a/2≤Zc≤Z1-a/2,原假設(shè)H0即被接受。另外MK 檢驗的另一個指標(biāo)Kendall 傾斜度β,用此指標(biāo)來量化單調(diào)趨勢,β即為:

    式中:1<j<i<n。當(dāng)β>0 時,反映上升的趨勢,反之則為下降的趨勢。

    3.2 Pettitt 突變分析法

    Pettitt 檢驗法是一種基于Mann-Whitney 非參數(shù)檢驗的方法,該方法可以對水文氣象要素序列進行突變分析獲得突變點,量化突變點在統(tǒng)計意義上的顯著水平[15]。其中Mann-Whitney 非參數(shù)統(tǒng)計量為[16]:

    根據(jù)統(tǒng)計量可計算:

    一般情況下,認為當(dāng)p≤0.05 時認為數(shù)據(jù)中存在突變點。

    3.3 小波分析法

    小波分析法將不同頻率成分的時間序列按尺度分解為低頻和高頻,再利用小波系數(shù)計算得到水文時間序列的周期性,選擇合理有效的Morlet 小波,對年徑流系列數(shù)據(jù)進行分析[17]。Morlet 小波表達式為:

    Morlet 小波的頻域函數(shù)為:

    采用上述Mann-Kendall 趨勢性分析方法分析1956 年~2020 年黃河沿、瑪曲和唐乃亥年徑流量的趨勢性,采用Pettitt 檢驗法分析突變性和突變點,采用小波分析方法中的Morlet 小波進行周期性分析。

    4 結(jié)果與分析

    4.1 趨勢性分析

    1956 年~2020 年黃河源區(qū)年徑流量變化過程見圖2,根據(jù)MK 趨勢檢驗結(jié)果可知,黃河沿站MK 檢驗統(tǒng)計量Zc值為-0.59,大于-1.96(0.05 顯著性水平),且傾斜度β為-0.03<0,因此,黃河沿站年徑流量呈現(xiàn)不顯著的減少趨勢;瑪曲站MK 趨勢檢驗可知,檢驗統(tǒng)計量Zc值為0.02,且小于1.96,傾斜度β為0.01>0,則呈現(xiàn)不顯著的增加趨勢;唐乃亥站檢驗統(tǒng)計量Zc值為0.03,傾斜度β為0.01>0,則呈現(xiàn)不顯著的增加趨勢。

    圖2 1956 年~2020 年黃河源區(qū)年徑流量變化趨勢圖

    4.2 突變性分析

    對1956 年~2020 年黃河源區(qū)年徑流量進行Pettitt 突變檢驗,檢驗結(jié)果見圖3,黃河沿站Pettitt 突變檢驗p值為0.17>0.05,則黃河沿站年徑流量在時間序列中沒有發(fā)生顯著的突變;瑪曲站Pettitt 突變檢驗p值為0.57>0.05,沒有發(fā)生顯著的突變;唐乃亥站Pettitt 突變檢驗p值為0.38>0.05,沒有發(fā)生顯著的突變。

    圖3 1956 年~2020 年黃河源區(qū)年徑流量Pettitt 突變檢驗結(jié)果圖

    4.3 周期性分析

    對1956 年~2020 年黃河源區(qū)年徑流量進行小波分析研究其周期性見圖4,黃河沿站第一主周期的周期大小為18 年,第二主周期的周期大小為6 年;瑪曲站第一主周期的周期大小為24 年,第二主周期的周期大小為7 年;唐乃亥站第一主周期的周期大小為24 年,第二主周期的周期大小為9 年。分別將周期尺度對應(yīng)的小波實部繪制出小波實部過程線(圖5、圖6)。黃河沿站在18 年的第一主周期上,大約經(jīng)歷了3個豐枯轉(zhuǎn)換期(圖5(a));在6 年的第二主周期上,大約經(jīng)歷了9 個豐枯轉(zhuǎn)換期(圖6(a))?,斍驹?4 年的第一主周期上,大約經(jīng)歷了3 個豐枯轉(zhuǎn)換期(圖5(b));在7 年的第二主周期上,大約經(jīng)歷了9 個豐枯轉(zhuǎn)換期(圖6(b))。唐乃亥站在24 年的第一主周期上,大約經(jīng)歷了2 個豐枯轉(zhuǎn)換期(圖5(c));在9 年的第二主周期上,大約經(jīng)歷了7 個豐枯轉(zhuǎn)換期(圖6(c))。

    圖4 1956 年~2020 年黃河源區(qū)年徑流量小波方差圖

    圖5 1956 年~2020 年黃河源區(qū)年徑流量第一主周期小波實部過程圖

    圖6 1956 年~2020 年黃河源區(qū)年徑流量第二主周期小波實部過程圖

    4.4 討論

    黃河源區(qū)黃河沿站1956 年~2020 年平均年徑流量為10.90 億m3,瑪曲站為144.52 億m3,唐乃亥站為203.72 億m3,可知黃河沿-瑪曲是黃河源區(qū)主要的產(chǎn)流區(qū),其徑流占黃河源區(qū)唐乃亥出口處徑流的66%。主要原因為西太平洋暖濕氣流從久治、瑪曲、若爾蓋一帶進入黃河源區(qū),黃河沿-瑪曲區(qū)域的多年平均降水量高于其他分區(qū),進而影響徑流量大小[18]。

    1956 年~2020 年黃河源區(qū)黃河沿站年徑流量呈現(xiàn)不顯著的減少趨勢,瑪曲站、唐乃亥站呈現(xiàn)不顯著的增加趨勢?;谒臍庀笳军c1956 年~2018 年實測降水和氣溫資料,采用考慮相關(guān)系數(shù)的距離平方反比法對降水和氣溫實測數(shù)據(jù)進行空間展布,得到3 個分區(qū)(黃河沿以上、黃河沿-瑪曲、瑪曲-唐乃亥)逐年的降水量[19]和氣溫[20]。采用Mann-Kendall 趨勢性分析法對三個分區(qū)進行分析,發(fā)現(xiàn)黃河沿以上年降水呈顯著增加趨勢,黃河沿-瑪曲、瑪曲-唐乃亥呈不顯著增加趨勢;而3 個區(qū)域年均氣溫均呈現(xiàn)顯著性增長趨勢,但黃河沿以上增加更快。判斷黃河沿以上蒸發(fā)增加的速度超過降水增加的速度,而黃河沿-瑪曲和瑪曲-唐乃亥區(qū)間蒸發(fā)增加的速度低于降水增加的速度,故導(dǎo)致不同站點徑流變化呈現(xiàn)不同的特征。

    5 結(jié)論

    為深入分析黃河源區(qū)1956 年~2020 年徑流量變化規(guī)律,選取Mann-Kendall 趨勢性分析法、Pettitt 突變分析法和小波分析法進行分析。由結(jié)果可知,1956 年~2020 年黃河源區(qū)中黃河沿年徑流量呈不顯著減少趨勢,瑪曲、唐乃亥呈現(xiàn)不顯著增加的趨勢。分析原因可能是因為降水與氣溫共同作用。在整個系列年中3 個站點無顯著突變點。3 個站小波分析均得到兩個主周期,黃河沿第一主周期的周期大小為18 年,第二主周期的周期大小為6 年;瑪曲第一主周期的周期大小為24 年,第二主周期的周期大小為7 年;唐乃亥第一主周期的周期大小為24 年,第二主周期的周期大小為9 年。

    猜你喜歡
    趨勢性瑪曲源區(qū)
    冬小麥蒸散源區(qū)代表性分析
    瑪曲:美麗的黃河首曲大草原
    農(nóng)產(chǎn)品流通競爭力的趨勢性和區(qū)域分異性
    渭河源區(qū)徑流量變化特征及趨勢分析
    基于M—K方法及成分提取的海洋表層熱狀分析
    瑪曲行吟
    中國三峽(2016年10期)2017-01-15 13:59:44
    在瑪曲
    紅豆(2016年9期)2016-05-14 11:11:49
    中國經(jīng)濟發(fā)展新常態(tài)十大趨勢性變化及因應(yīng)策略
    全球化(2015年9期)2015-02-28 12:40:07
    近57年來和豐縣氣溫和降水量的趨勢性及突變特征近
    問候瑪曲
    江南詩(2013年4期)2013-10-26 03:43:25
    av又黄又爽大尺度在线免费看| 亚洲五月色婷婷综合| 亚洲专区中文字幕在线 | 国产精品一区二区在线不卡| 高清欧美精品videossex| 如何舔出高潮| 亚洲av成人不卡在线观看播放网 | 观看av在线不卡| 亚洲av成人精品一二三区| 操美女的视频在线观看| 一级毛片 在线播放| www.精华液| 桃花免费在线播放| 日本vs欧美在线观看视频| 男人操女人黄网站| 国产一区二区在线观看av| 日本色播在线视频| 国产麻豆69| 最新在线观看一区二区三区 | 久久免费观看电影| 只有这里有精品99| 美女主播在线视频| 国产精品亚洲av一区麻豆 | 国产一区二区激情短视频 | 亚洲av国产av综合av卡| 青春草亚洲视频在线观看| 精品少妇久久久久久888优播| 一本大道久久a久久精品| 侵犯人妻中文字幕一二三四区| 久久久久久久国产电影| 国产一区二区三区综合在线观看| 香蕉国产在线看| 中国国产av一级| 久热爱精品视频在线9| 亚洲第一av免费看| 欧美日韩亚洲综合一区二区三区_| 人人妻人人澡人人看| 久久精品久久精品一区二区三区| 国产精品蜜桃在线观看| videosex国产| 最近最新中文字幕免费大全7| 亚洲婷婷狠狠爱综合网| 亚洲国产毛片av蜜桃av| 在线看a的网站| 性少妇av在线| av不卡在线播放| 亚洲伊人色综图| 国产不卡av网站在线观看| 在线天堂最新版资源| 男女边摸边吃奶| 日本色播在线视频| 国产成人啪精品午夜网站| 国产日韩一区二区三区精品不卡| 久久影院123| 精品少妇久久久久久888优播| 黄色毛片三级朝国网站| 99久国产av精品国产电影| 丝袜脚勾引网站| 国产精品国产三级国产专区5o| 中文字幕亚洲精品专区| 黑人巨大精品欧美一区二区蜜桃| 一级爰片在线观看| 熟女av电影| 日韩熟女老妇一区二区性免费视频| 国产高清不卡午夜福利| 一级片免费观看大全| 亚洲av福利一区| 亚洲精品国产av蜜桃| 日韩电影二区| 丝袜人妻中文字幕| 人人澡人人妻人| 免费看av在线观看网站| 看十八女毛片水多多多| 少妇人妻 视频| 99久久综合免费| 97在线人人人人妻| 男的添女的下面高潮视频| 亚洲国产中文字幕在线视频| 一级a爱视频在线免费观看| 赤兔流量卡办理| 最新的欧美精品一区二区| 久久影院123| 另类精品久久| 国产高清不卡午夜福利| www.精华液| 精品少妇一区二区三区视频日本电影 | 1024香蕉在线观看| 又大又黄又爽视频免费| 不卡av一区二区三区| 日韩 欧美 亚洲 中文字幕| 亚洲第一区二区三区不卡| 久久久久国产一级毛片高清牌| 国产片特级美女逼逼视频| 成年动漫av网址| 日韩av在线免费看完整版不卡| 青春草亚洲视频在线观看| 蜜桃国产av成人99| 国产av码专区亚洲av| 爱豆传媒免费全集在线观看| 男人爽女人下面视频在线观看| 国产av一区二区精品久久| 亚洲精品视频女| 丰满迷人的少妇在线观看| 欧美亚洲日本最大视频资源| 老司机影院毛片| 亚洲精品乱久久久久久| 国产亚洲av高清不卡| 高清视频免费观看一区二区| av片东京热男人的天堂| 最近的中文字幕免费完整| 宅男免费午夜| 操美女的视频在线观看| 久久精品亚洲熟妇少妇任你| 一区福利在线观看| 99精品久久久久人妻精品| 好男人视频免费观看在线| 丁香六月天网| 欧美日韩成人在线一区二区| 亚洲国产欧美网| 日本91视频免费播放| 国产精品av久久久久免费| 街头女战士在线观看网站| 亚洲国产欧美日韩在线播放| 在线观看一区二区三区激情| 国产精品二区激情视频| 美女大奶头黄色视频| 99热国产这里只有精品6| 一级片'在线观看视频| 免费观看性生交大片5| 街头女战士在线观看网站| 久久精品国产亚洲av涩爱| 午夜福利在线免费观看网站| 国产伦人伦偷精品视频| 日韩伦理黄色片| 欧美国产精品一级二级三级| 久久 成人 亚洲| 一级黄片播放器| 91成人精品电影| 日韩欧美一区视频在线观看| 五月开心婷婷网| 天堂中文最新版在线下载| 只有这里有精品99| 大码成人一级视频| 91精品伊人久久大香线蕉| 亚洲五月色婷婷综合| 国产亚洲av片在线观看秒播厂| 老汉色av国产亚洲站长工具| 亚洲一区二区三区欧美精品| 成人国产av品久久久| 另类精品久久| 亚洲欧美日韩另类电影网站| 97精品久久久久久久久久精品| 另类精品久久| 婷婷成人精品国产| 麻豆乱淫一区二区| 亚洲国产欧美日韩在线播放| 国产亚洲av高清不卡| 新久久久久国产一级毛片| a级片在线免费高清观看视频| 欧美激情极品国产一区二区三区| 天天操日日干夜夜撸| 视频区图区小说| 成人影院久久| 国产精品久久久久成人av| 婷婷色av中文字幕| 亚洲av电影在线进入| 亚洲综合色网址| 黄网站色视频无遮挡免费观看| 性高湖久久久久久久久免费观看| 在线观看三级黄色| 纯流量卡能插随身wifi吗| 中文字幕人妻丝袜制服| 成人18禁高潮啪啪吃奶动态图| 男人舔女人的私密视频| 国产成人精品无人区| 黄片播放在线免费| 精品亚洲成a人片在线观看| av网站免费在线观看视频| 综合色丁香网| 在线观看免费视频网站a站| av免费观看日本| 高清欧美精品videossex| 亚洲精品一二三| 视频区图区小说| 国产老妇伦熟女老妇高清| 国产不卡av网站在线观看| 少妇人妻 视频| 欧美日韩精品网址| 亚洲av电影在线进入| 多毛熟女@视频| 久久久久久久大尺度免费视频| 欧美日韩亚洲综合一区二区三区_| 男女国产视频网站| 日日啪夜夜爽| 大片电影免费在线观看免费| 悠悠久久av| 亚洲第一区二区三区不卡| 青春草亚洲视频在线观看| 精品亚洲乱码少妇综合久久| 香蕉国产在线看| 夜夜骑夜夜射夜夜干| 亚洲,欧美精品.| 99久久精品国产亚洲精品| 免费黄色在线免费观看| 五月开心婷婷网| 亚洲国产中文字幕在线视频| 男女午夜视频在线观看| 亚洲天堂av无毛| 99久国产av精品国产电影| 精品久久蜜臀av无| 成年av动漫网址| 欧美精品人与动牲交sv欧美| 精品久久久精品久久久| 久久亚洲国产成人精品v| 亚洲伊人久久精品综合| 热re99久久精品国产66热6| 人体艺术视频欧美日本| 久久99热这里只频精品6学生| 日本91视频免费播放| 青草久久国产| 亚洲一区中文字幕在线| 精品人妻熟女毛片av久久网站| 一本—道久久a久久精品蜜桃钙片| 亚洲熟女毛片儿| 在线天堂最新版资源| 91老司机精品| 国产精品国产av在线观看| 国产一区二区在线观看av| 韩国av在线不卡| 最近中文字幕2019免费版| 中文字幕最新亚洲高清| 精品视频人人做人人爽| 一级毛片 在线播放| 妹子高潮喷水视频| 一级,二级,三级黄色视频| 王馨瑶露胸无遮挡在线观看| 亚洲,欧美,日韩| av又黄又爽大尺度在线免费看| 老司机靠b影院| 中文天堂在线官网| 欧美 亚洲 国产 日韩一| 亚洲国产日韩一区二区| 国产高清国产精品国产三级| 亚洲国产看品久久| 99热国产这里只有精品6| 久久久久国产精品人妻一区二区| 亚洲欧美一区二区三区黑人| 中文字幕高清在线视频| 免费在线观看完整版高清| 国产成人欧美在线观看 | 亚洲成人国产一区在线观看 | 在现免费观看毛片| 国产成人精品在线电影| 国产99久久九九免费精品| 亚洲精品久久成人aⅴ小说| 亚洲美女视频黄频| 国产精品蜜桃在线观看| 亚洲美女黄色视频免费看| 亚洲国产精品成人久久小说| 久久亚洲国产成人精品v| 欧美人与性动交α欧美精品济南到| 啦啦啦 在线观看视频| 精品国产乱码久久久久久男人| av福利片在线| 午夜免费鲁丝| 精品国产超薄肉色丝袜足j| 久久毛片免费看一区二区三区| 婷婷色麻豆天堂久久| 亚洲情色 制服丝袜| 精品少妇一区二区三区视频日本电影 | 777久久人妻少妇嫩草av网站| 亚洲国产看品久久| 母亲3免费完整高清在线观看| 日本vs欧美在线观看视频| 婷婷色综合www| 一级爰片在线观看| 久久久亚洲精品成人影院| 亚洲av电影在线观看一区二区三区| 在线观看人妻少妇| 99久久精品国产亚洲精品| 亚洲精品国产区一区二| 色精品久久人妻99蜜桃| 亚洲国产精品成人久久小说| 99国产综合亚洲精品| 亚洲色图综合在线观看| 在线精品无人区一区二区三| 精品第一国产精品| 亚洲欧美成人精品一区二区| 精品午夜福利在线看| 丝瓜视频免费看黄片| 黄片无遮挡物在线观看| 国产日韩欧美视频二区| 久久天躁狠狠躁夜夜2o2o | 巨乳人妻的诱惑在线观看| 日本wwww免费看| 另类精品久久| 亚洲美女黄色视频免费看| 国产精品一区二区精品视频观看| 午夜福利网站1000一区二区三区| 欧美人与性动交α欧美软件| 狠狠婷婷综合久久久久久88av| 一区二区三区四区激情视频| 亚洲精品美女久久久久99蜜臀 | 亚洲精品第二区| 亚洲七黄色美女视频| 国产免费一区二区三区四区乱码| 久久女婷五月综合色啪小说| 男男h啪啪无遮挡| 久久毛片免费看一区二区三区| 亚洲综合色网址| 精品亚洲成a人片在线观看| 国产欧美日韩综合在线一区二区| 欧美 日韩 精品 国产| 在线观看免费日韩欧美大片| 日韩中文字幕欧美一区二区 | 久久午夜综合久久蜜桃| 成人午夜精彩视频在线观看| 男女高潮啪啪啪动态图| 人妻 亚洲 视频| av网站在线播放免费| 无遮挡黄片免费观看| 又大又黄又爽视频免费| 欧美 日韩 精品 国产| 成人手机av| 国产精品免费视频内射| 国产成人精品久久二区二区91 | 女性生殖器流出的白浆| 一二三四中文在线观看免费高清| 亚洲中文av在线| 成人18禁高潮啪啪吃奶动态图| 80岁老熟妇乱子伦牲交| 一级毛片电影观看| kizo精华| 中国国产av一级| 国产精品国产av在线观看| 午夜久久久在线观看| 精品一区二区三卡| 国产男女内射视频| 久久精品aⅴ一区二区三区四区| 国产精品人妻久久久影院| 丁香六月天网| 日韩中文字幕视频在线看片| 男人添女人高潮全过程视频| 欧美日韩亚洲国产一区二区在线观看 | 丝袜在线中文字幕| 99香蕉大伊视频| 国产一级毛片在线| 美女主播在线视频| 熟妇人妻不卡中文字幕| 菩萨蛮人人尽说江南好唐韦庄| 18禁观看日本| 国产国语露脸激情在线看| 欧美亚洲日本最大视频资源| 亚洲av中文av极速乱| 一二三四中文在线观看免费高清| 欧美在线黄色| 99精品久久久久人妻精品| 亚洲精品国产区一区二| 欧美精品一区二区免费开放| 激情视频va一区二区三区| 十八禁人妻一区二区| 少妇 在线观看| 1024视频免费在线观看| 国产精品国产三级国产专区5o| 午夜日本视频在线| 青春草国产在线视频| 国产精品国产三级专区第一集| 97精品久久久久久久久久精品| 日韩一本色道免费dvd| 日韩 亚洲 欧美在线| 一级片免费观看大全| 黄色视频在线播放观看不卡| 1024视频免费在线观看| 在线观看www视频免费| 青草久久国产| 在线亚洲精品国产二区图片欧美| 精品一区二区三区四区五区乱码 | 午夜福利乱码中文字幕| 欧美中文综合在线视频| 80岁老熟妇乱子伦牲交| 欧美久久黑人一区二区| 精品人妻熟女毛片av久久网站| 久久久久久久精品精品| 亚洲中文av在线| 中文字幕高清在线视频| 中文字幕另类日韩欧美亚洲嫩草| 日本欧美国产在线视频| 男的添女的下面高潮视频| 九色亚洲精品在线播放| 超色免费av| 激情视频va一区二区三区| 午夜av观看不卡| 晚上一个人看的免费电影| av福利片在线| 精品人妻一区二区三区麻豆| 国产一区亚洲一区在线观看| 欧美黑人精品巨大| 9色porny在线观看| 啦啦啦 在线观看视频| 久久韩国三级中文字幕| 日韩熟女老妇一区二区性免费视频| 久久99热这里只频精品6学生| 国产成人午夜福利电影在线观看| 国产精品人妻久久久影院| av在线老鸭窝| 欧美中文综合在线视频| 国产毛片在线视频| av国产精品久久久久影院| 黄色一级大片看看| 美国免费a级毛片| 51午夜福利影视在线观看| a级毛片黄视频| 1024视频免费在线观看| 老汉色∧v一级毛片| 99热国产这里只有精品6| 亚洲图色成人| 久久久久视频综合| 久热这里只有精品99| 精品一区在线观看国产| 久久精品国产亚洲av涩爱| 啦啦啦视频在线资源免费观看| 国产免费现黄频在线看| 亚洲四区av| 亚洲国产欧美一区二区综合| 另类精品久久| 极品人妻少妇av视频| 国产野战对白在线观看| 亚洲欧洲日产国产| 尾随美女入室| 热re99久久国产66热| 欧美激情高清一区二区三区 | 哪个播放器可以免费观看大片| 久久久久久人人人人人| 国产成人91sexporn| 18禁观看日本| 制服诱惑二区| 一二三四在线观看免费中文在| 国产欧美亚洲国产| 男女高潮啪啪啪动态图| 久久青草综合色| kizo精华| 啦啦啦在线观看免费高清www| 999久久久国产精品视频| 久久免费观看电影| 一本大道久久a久久精品| 欧美激情极品国产一区二区三区| 黑人猛操日本美女一级片| 人人妻,人人澡人人爽秒播 | 国产在线免费精品| 亚洲国产精品一区三区| 精品国产一区二区三区四区第35| a级片在线免费高清观看视频| 丰满乱子伦码专区| 国产欧美亚洲国产| 2021少妇久久久久久久久久久| 久久免费观看电影| 18禁动态无遮挡网站| 操美女的视频在线观看| 秋霞伦理黄片| 久久久久久久精品精品| 叶爱在线成人免费视频播放| 国产精品久久久久久精品电影小说| 超色免费av| 成年av动漫网址| 少妇的丰满在线观看| 精品国产一区二区三区四区第35| 欧美久久黑人一区二区| 国产男女超爽视频在线观看| 成年人午夜在线观看视频| 人体艺术视频欧美日本| 丰满少妇做爰视频| 成人亚洲欧美一区二区av| 亚洲色图综合在线观看| 国产欧美亚洲国产| 高清不卡的av网站| 亚洲中文av在线| av在线播放精品| 欧美黑人欧美精品刺激| 黄片播放在线免费| 亚洲久久久国产精品| 久久午夜综合久久蜜桃| 欧美老熟妇乱子伦牲交| 成人漫画全彩无遮挡| 午夜免费男女啪啪视频观看| 久久ye,这里只有精品| 欧美变态另类bdsm刘玥| 久久久久精品国产欧美久久久 | 亚洲av福利一区| 天天躁狠狠躁夜夜躁狠狠躁| 成年女人毛片免费观看观看9 | svipshipincom国产片| 久久国产精品男人的天堂亚洲| 国产成人免费观看mmmm| 男女午夜视频在线观看| 成人国产av品久久久| 日韩 欧美 亚洲 中文字幕| 爱豆传媒免费全集在线观看| 亚洲欧美色中文字幕在线| 成人国产av品久久久| 人人妻人人爽人人添夜夜欢视频| 亚洲国产精品一区三区| 国产成人一区二区在线| 99九九在线精品视频| 电影成人av| 女的被弄到高潮叫床怎么办| 777米奇影视久久| 一本一本久久a久久精品综合妖精| 中文字幕高清在线视频| 母亲3免费完整高清在线观看| 悠悠久久av| 可以免费在线观看a视频的电影网站 | 菩萨蛮人人尽说江南好唐韦庄| 国产精品香港三级国产av潘金莲 | 国产精品久久久久久精品古装| 精品一区二区免费观看| 久久毛片免费看一区二区三区| 久久 成人 亚洲| 国产成人精品久久二区二区91 | 1024视频免费在线观看| 女性被躁到高潮视频| 国产精品偷伦视频观看了| 日日啪夜夜爽| 国产高清国产精品国产三级| 国产男女超爽视频在线观看| 久久久久国产精品人妻一区二区| 午夜老司机福利片| 午夜福利一区二区在线看| 日日摸夜夜添夜夜爱| 满18在线观看网站| 国产成人欧美| 亚洲国产欧美在线一区| 久久久久精品人妻al黑| 女性被躁到高潮视频| av免费观看日本| 久久国产精品男人的天堂亚洲| 色视频在线一区二区三区| 免费女性裸体啪啪无遮挡网站| 韩国高清视频一区二区三区| 国产成人91sexporn| 热99久久久久精品小说推荐| 亚洲婷婷狠狠爱综合网| 丝袜人妻中文字幕| av在线app专区| 搡老乐熟女国产| 人体艺术视频欧美日本| 国产精品一二三区在线看| 精品一区二区免费观看| 亚洲精华国产精华液的使用体验| 国产成人a∨麻豆精品| 久久久亚洲精品成人影院| 老熟女久久久| 制服人妻中文乱码| 99精国产麻豆久久婷婷| 韩国高清视频一区二区三区| 亚洲精品乱久久久久久| 亚洲精品美女久久久久99蜜臀 | 久久婷婷青草| 国产免费视频播放在线视频| 日韩av不卡免费在线播放| 99久久人妻综合| 国产精品一区二区在线不卡| 欧美激情高清一区二区三区 | 日韩大片免费观看网站| 欧美人与善性xxx| 免费日韩欧美在线观看| 精品国产乱码久久久久久小说| 秋霞伦理黄片| 99久久99久久久精品蜜桃| 国产极品天堂在线| av在线老鸭窝| 一区二区三区四区激情视频| 母亲3免费完整高清在线观看| 丝瓜视频免费看黄片| 久久ye,这里只有精品| 一级毛片黄色毛片免费观看视频| 蜜桃国产av成人99| 国产男女内射视频| 欧美最新免费一区二区三区| 一区二区三区激情视频| 各种免费的搞黄视频| 久久鲁丝午夜福利片| 曰老女人黄片| 极品人妻少妇av视频| 午夜福利,免费看| 五月天丁香电影| 人人妻人人澡人人看| 宅男免费午夜| 老司机亚洲免费影院| 国产国语露脸激情在线看| 1024香蕉在线观看| 日本欧美国产在线视频| 99热国产这里只有精品6| 久久久欧美国产精品| 欧美人与性动交α欧美软件| 99精品久久久久人妻精品| 国产av一区二区精品久久| 精品久久久精品久久久| 哪个播放器可以免费观看大片| 一区二区日韩欧美中文字幕| 咕卡用的链子| 国产 一区精品| 国产淫语在线视频| 99九九在线精品视频| 如日韩欧美国产精品一区二区三区| 美女午夜性视频免费| 91精品伊人久久大香线蕉| 只有这里有精品99| 男人添女人高潮全过程视频| 亚洲婷婷狠狠爱综合网| 香蕉国产在线看| 日韩一区二区视频免费看| 人妻人人澡人人爽人人| 欧美人与性动交α欧美精品济南到| 考比视频在线观看| 国产精品久久久久成人av|