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

    流域水文氣象變化趨勢及相關性分析

    2016-09-01 07:54:54曹小磊周祖昊王子茹中國水利水電科學研究院流域水循環(huán)模擬與調控國家重點實驗室北京0008大連理工大學水資源與防洪研究所遼寧大連6086沈陽環(huán)境科學研究院遼寧沈陽006
    東北水利水電 2016年8期
    關鍵詞:松花江流域氣象要素徑流

    曹小磊,周祖昊,許 翼,王子茹(.中國水利水電科學研究院 流域水循環(huán)模擬與調控國家重點實驗室,北京 0008;.大連理工大學水資源與防洪研究所,遼寧 大連 6086;.沈陽環(huán)境科學研究院,遼寧 沈陽 006)

    流域水文氣象變化趨勢及相關性分析

    曹小磊1,2,周祖昊1,許翼3,王子茹2
    (1.中國水利水電科學研究院 流域水循環(huán)模擬與調控國家重點實驗室,北京 100038;
    2.大連理工大學水資源與防洪研究所,遼寧 大連 116086;3.沈陽環(huán)境科學研究院,遼寧 沈陽 110016)

    采用Mann-Kendall趨勢檢驗和線性回歸分析法,對松花江流域及水資源二級區(qū)1956-2010年的水文氣象要素進行趨勢分析,并采用Pearson相關系數(shù)方法,分析了研究區(qū)天然徑流與氣象要素之間的相關關系。趨勢分析結果表明:松花江流域的氣溫呈顯著上升趨勢,降水呈不顯著的減少趨勢,日照時數(shù)、相對濕度和風速呈顯著的減小趨勢,潛在蒸散發(fā)呈顯著增加趨勢,天然徑流深呈不顯著的減少趨勢。相關性分析結果表明,影響天然徑流變化的主要因素是降水量,蒸散發(fā)變化也對徑流有較為顯著的影響作用,相對濕度減少和溫度升高也在一定程度上通過對蒸散發(fā)的影響間接作用于徑流變化。

    水文氣象要素;Mann-Kendall;趨勢分析;相關性分析

    0 引言

    氣候變化改變了地表的水熱條件,使全球水資源供給條件發(fā)生變化,進而影響地表的生態(tài)水文過程,以及社會經(jīng)濟發(fā)展和生態(tài)環(huán)境建設[1]。由于區(qū)域條件的差異,氣候變化情況具有不確定性,根據(jù)IPCC報告[2],高緯度地區(qū)屬于氣候變化的敏感區(qū)域。寒區(qū)生態(tài)系統(tǒng)對全球變化的高度敏感性導致生態(tài)系統(tǒng)加速退化[3],在氣候變化及生態(tài)系統(tǒng)變化驅動作用下,寒區(qū)的水文循環(huán)過程發(fā)生了顯著改變。松花江流域是我國氣候變化敏感區(qū),近年來,受氣候變化影響,松花江流域干旱洪澇頻發(fā),水土流失嚴重[4]。目前針對松花江流域氣候變化及其影響的研究大部分集中在單個氣象、水文要素或局部區(qū)域的研究[4-9],缺乏對全流域氣象、水文要素變化以及水文要素與各氣象因子之間相關關系的系統(tǒng)全面的研究分析?;?956—2010年長系列氣象水文資料,對松花江流域的氣象、水文要素的變化趨勢進行了分析,并探討了影響流域徑流變化的主要氣象因素,以期為松花江流域的水資源可持續(xù)利用、生態(tài)建設和社會經(jīng)濟發(fā)展提供科學依據(jù)。

    1 研究區(qū)概況

    松花江流域地處我國東北地區(qū)的北部,流域面積55.68萬km2。松花江有南、北兩源,南源為第二松花江,流域面積7.82萬 km2;北源為嫩江,流域面積28.3萬km2[6],兩源于三岔河附近匯合向東而流注入黑龍江。松花江流域地處北溫帶季風氣候區(qū),多年平均降水量在500 mm左右,東南部山區(qū)降水可達700~900 mm,而干旱的西部地區(qū)只有400 mm,流域多年平均水資源總量為960.9億m3,多年平均氣溫在-5~3℃之間。

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

    氣象資料來自中國氣象科學數(shù)據(jù)共享服務網(wǎng),包括松花江流域內(nèi)的46個氣象站和周邊18個氣象站的1956—2010年的逐日降水、氣溫、日照時數(shù)、相對濕度和風速數(shù)據(jù)。采用距離平方反比法(RDS)對松花江流域內(nèi)的氣象數(shù)據(jù)進行空間插值,統(tǒng)計得到各水資源分區(qū)的平均值。徑流資料來源于松遼流域水資源綜合評價的成果,具體收集了松花江流域1956—2000年45年的天然徑流深數(shù)據(jù),并依據(jù)1999—2010年的《松遼流域水資源公報》將天然徑流數(shù)據(jù)延長至2010年。

    采用線性回歸分析法和Mann-Kendall趨勢分析法[10],分析松花江流域的氣象要素和水文要素的年際變化趨勢,在顯著性水平α=0.05時,Mann-Kendall檢驗統(tǒng)計量的臨界值為±1.96;采用pearson相關系數(shù)方法[11],分析天然徑流與各氣象要素的相關關系。

    3 結果分析

    3.1松花江流域氣象要素變化趨勢分析

    對松花江流域及水資源二級區(qū)1956—2010年的降水、平均氣溫、日照時數(shù)、相對濕度和風速的變化趨勢,進行Mann-Kendall趨勢線檢驗和線性回歸分析,如表1所示。松花江流域降水呈不顯著的減少趨勢,平均減少幅度為0.80 mm/a;嫩江、第二松花江和松花江(三岔河口以下)的降水均呈不顯著減少趨勢,平均減少幅度分別為0.61 mm/年、0.52 mm/a和1.21 mm/a。松花江流域氣溫呈顯著上升趨勢,平均升幅為0.39℃/10 a;嫩江、第二松花江和松花江(三岔河口以下)的氣溫均呈顯著上升趨勢,平均上升幅度分別為0.42℃/1 0a、0.38℃/10 a和0.27℃/10 a。松花江流域日照時數(shù)、相對濕度和風速均呈顯著的減少趨勢,各水資源分區(qū)的日照時數(shù)、相對濕度和風速也均為顯著減少趨勢,只是受地理位置、下墊面條件等因素影響,減少程度有所不同。3.2松花江流域潛在蒸散發(fā)趨勢分析

    表1 松花江流域氣象要素趨勢檢驗結果

    蒸散發(fā)是能量平衡和水量平衡的重要組成部分,是水文循環(huán)的重要環(huán)節(jié),潛在蒸散發(fā)是實際蒸散量的理論上限,通常也是計算實際蒸散量的基礎[12]。根據(jù)插值得到的日氣象資料和彭曼公式[13],計算松花江流域的多年平均潛在蒸散發(fā)為659.53 mm,其中嫩江、第二松花江和松花江(三岔河口以下)的潛在蒸散發(fā)多年平均值分別為667.98 mm、 670.32 mm和645.55 mm。松花江流域潛在蒸散發(fā)的年際變化如圖1所示,Mann-Kendall趨勢檢驗和線性回歸分析結果顯示其呈顯著增加趨勢,平均增幅為0.75 mm/a,見表2。從水資源分區(qū)來看,嫩江的潛在蒸散發(fā)呈顯著增加趨勢,平均增幅為0.95 mm/a,第二松花江和松花江(三岔河口以下)的潛在蒸散發(fā)呈不顯著增加趨勢,平均增幅為0.47 mm/a和0.53 mm/a。 3.3松花江流域徑流趨勢分析

    圖1 松花江流域潛在蒸散發(fā)年際變化

    表2 松花江流域潛在蒸散發(fā)趨勢檢驗結果

    松花江流域的多年平均天然徑流深為142.62 mm,徑流系數(shù)為0.26,利用Mann-Kendall趨勢檢驗法和線性回歸分析方法,對天然徑流的變化趨勢進行分析,結果如圖2所示,松花江流域1956—2010年的徑流呈波動減少趨勢,平均減少幅度為0.53 mm/a,Mann-Kendall趨勢檢驗得到的統(tǒng)計值為-1.29,即不顯著減少趨勢。

    圖2 松花江流域徑流年際變化過程

    表3給出了松花江流域1956—2010年徑流的平均值、徑流系數(shù)、變差系數(shù)和趨勢檢驗結果。嫩江、第二松花江和松花江(三岔河口以下)多年平均天然徑流深分別為94.30 mm、223.16 mm和185.50 mm。從徑流變差系數(shù)來看,嫩江最大,表明嫩江徑流年際變化較大,而第二松花江和松花江(三岔河口以下)徑流變化則相對較小。Man-Kendall趨勢檢驗結果表明:嫩江的天然徑流呈顯著減少趨勢,平均減少幅度為0.48 mm/a;第二松花江和松花江(三岔河口以下)的天然徑流則呈不顯著減少趨勢,平均減少幅度分別為0.11 mm/a 和0.76 mm/a。3.4松花江流域徑流與氣象要素的相關性分析

    表3 松花江流域徑流統(tǒng)計量和趨勢檢驗結果

    分析松花江流域1956—2010年的天然徑流和降水的相關關系,結果見表4,徑流與降水存在顯著的正相關關系,相關系數(shù)為0.93??梢姡苫ń饔蚪鼛资陙韽搅鳒p少,主要是由于降水減少引起的。從水資源分區(qū)來看,嫩江、第二松花江、松花江(三岔河口以下)徑流和降水均呈顯著的正相關關系,相關系數(shù)分別為0.89,0.92和0.93。

    對松花江流域的徑流和潛在蒸散發(fā)進行相關性分析,結果如表4所示,徑流與潛在蒸散發(fā)呈顯著的負相關關系,相關系數(shù)為-0.72。

    表4 松花江流域徑流與氣象要素的相關系數(shù)

    可見,松花江流域蒸散發(fā)變化對徑流的影響雖然不及降水,但也較為顯著,是徑流變化的因素之一。潛在蒸散發(fā)增加的趨勢表明地區(qū)的蒸散發(fā)能力呈逐漸增大趨勢,更多的降水將以蒸散發(fā)的形式返回到大氣中,導致地區(qū)產(chǎn)流系數(shù)的下降。從水資源分區(qū)來看,嫩江、第二松花江、松花江(三岔河口以下)徑流和潛在蒸散發(fā)呈顯著的負相關關系,相關系數(shù)分別為-0.64,-0.59和-0.74。

    區(qū)域蒸散發(fā)又受氣溫、日照、相對濕度和風速4個要素的影響,各氣象要素通過對蒸散發(fā)的影響間接作用于徑流。通過分析氣溫、日照、相對濕度、風速和潛在蒸散發(fā)的相關關系,可以進一步明確各氣象要素對徑流的影響程度。

    表5 松花江流域潛在蒸散發(fā)與氣象要素的相關系數(shù)

    由表5的分析結果可以看出,松花江流域的氣溫、日照、相對濕度、風速和潛在蒸散發(fā)均存在較好的相關性,其中,蒸散發(fā)與相對濕度的相關性最高,相關系數(shù)為-0.78;其次是溫度,相關系數(shù)為0.52??梢?,溫度和相對濕度是影響松花江流域蒸散發(fā)變化的主要氣象因素,同時也在一定程度上引起徑流變化。日照時數(shù)和風速對潛在蒸散發(fā)的影響作用相對較弱,相關系數(shù)分別為0.36和-0.30,受其他氣象因素的影響作用,松花江流域風速與潛在蒸散發(fā)呈負相關關系。水資源分區(qū)各氣象要素與潛在蒸散發(fā)的相關關系基本與松花江流域相同。

    4 結論

    分析松花江流域近195—2010年的水文氣象要素變化趨勢,以及影響松花江流域徑流變化的主要因素,得出:

    松花江流域的氣溫呈上升趨勢,降水、日照時數(shù)、相對濕度和風速呈下降趨勢,潛在蒸散發(fā)呈增加趨勢,天然徑流深呈下降趨勢。其中,溫度的平均上升幅度為0.39℃/10 a,降水的平均減少幅度為0.80 mm/a,潛在蒸散發(fā)的平均增加幅度為0.75 mm/a,徑流的平均減少幅度為0.53 mm/a。

    通過對松花江流域徑流、潛在蒸發(fā)和氣象要素之間的相關性分析,表明影響天然徑流下降的主要原因是受降水減少的影響,此外,蒸散發(fā)量增加也是松花江流域天然徑流減少的主要原因之一;其他四個氣象要素主要通過對蒸散發(fā)的影響間接作用于地區(qū)的徑流量,而氣象要素和潛在蒸散發(fā)之間的相關性分析表明,引起松花江流域蒸

    發(fā)變化的主要原因是相對濕度和氣溫。

    [1]王寶鑒,宋連春,張強,等.石羊河流域水資源對氣候變暖的響應及對生態(tài)環(huán)境的影響[J].地球科學進展,2007,2(7):730-737.

    [2]IPCC.Climatechange2007:Synthesisreport[M].Cambridge:Cambridge University Press,2007.

    [3]王根緒,李元壽,王一博,等.長江源區(qū)高寒生態(tài)與氣候變化對河流徑流過程的影響分析[J].冰川凍土,2007,29(2):159-169.

    [4]俞方圓,鄭粉莉,李志,等.近51年松花江流域溫度變化趨勢分析[J].干旱地區(qū)農(nóng)業(yè)研究,2011,29(5):242-252.

    [5]李想,李維京,趙振國.我國松花江流域和遼河流域降水的長期變化規(guī)律和未來趨勢分析[J].應用氣象學報,2005,16(5):593-599.

    [6]王彥君,王隨繼,蘇騰.1955-2010年松花江流域不同區(qū)段徑流量變化影響因素定量評估[J].地理科學進展,2014,33(1):65-75.

    [7]潘健,唐莉華.松花江流域上游徑流變化及其影響研究[J].水力發(fā)電學報,2013,32(5):58-63,69.

    [8]謝今范,韋小麗,張晨琛,等.第二松花江流域實際蒸散發(fā)的時空變化特征和影響因素[J].生態(tài)學雜志,2013,32(12):3336-3343.

    [9]溫姍姍,姜彤,李修倉,等.1961-2010年松花江流域實際蒸散發(fā)時空變化及影響要素分析[J].氣候變化研究進展,2014,10(2):79-86.

    [10]李占玲,王小娟.黑河中上游氣象要素單調變化趨勢分析[J].地理研究,2011,30(11):2059-2066.

    [11]Pearson K.On some novel properties of partial and multiple correlationcoefficientsinauniverseofmanifoldcharacteristics [J].Biometrika,1916,11(3):231-238.

    [12]高歌,陳德亮,任國玉,等.1956~2000年中國潛在蒸散量變化趨勢[J].地理研究,2006,25(3):378-387.

    [13]Penman H.L.Natural evaportranspiration from open water,bare soil and grass[J].Pro R Soc Lond A,1948,193:120-145.流域上下游各水文站設計洪水成果點繪在地區(qū)綜合圖上進行分析,洪峰、洪量均值隨集水面積的增大而增大,Cv值隨集水面積的增大而減小,各設計成果在地區(qū)分布上是合理的。另外,通過將桓仁、回龍山站的設計洪水成果點繪在實測峰量關系圖上進行分析,設計峰量值與實測點據(jù)的分布趨勢較為協(xié)調。因此,本次計算的設計洪水參數(shù)是合理的。5.3成果對比分析

    將這次計算的桓仁、回龍山、太平哨水電站壩址設計洪水成果與1972年太平哨初設審定成果比較,均值減小,Cv值增大,各頻率設計洪水成果均有所減小。其中洪峰流量設計成果減小幅度在0%~8%之間,3 d洪量成果減小幅度在4%~10%之間,分析其原因,主要為1971年以后渾江流域雖然發(fā)生了1986,1995,2010年等大水,但其量級相對于1888,1960年等歷史洪水還有一定差距,對前4位大洪水的排位無影響,而洪水系列的延長導致大洪水重現(xiàn)期增加,致使頻率曲線中的大洪水點距左移,且本次延長系列中大多數(shù)年份洪水量級不大,導致均值減小,從而使本次洪峰、洪量設計

    值較太平哨初設成果整體減小。

    6 結語

    通過對桓仁、回龍山、太平哨水電站壩址設計洪水進行復核分析,各壩址設計洪峰、洪量成果較審定成果有一定程度的減小,因此,從工程安全角度分析,太平哨初設審定的洪水成果仍然是安全的??紤]各水電站工程均已建成,為保持工程設計洪水成果相對穩(wěn)定、且從工程安全考慮,桓仁、回龍山、太平哨壩址洪水仍可沿用既往審定成果。

    [參考文獻]

    [收稿日期]2016-04-20

    [1]常景坤,鄧秋良等.受上游水庫群影響的設計洪水分析研究[J].中國水運,2014(11):240-241.

    [2]L278-2002,水利水電工程水文計算規(guī)范[S].北京:中國水利水電出版社,2002.

    [3]SL44-2006,水利水電工程設計洪水計算規(guī)范[S].北京:中國水利水電出版社,2006.

    TV12

    B

    1002-0624(2016)08-0029-03

    2016-03-07

    猜你喜歡
    松花江流域氣象要素徑流
    成都電網(wǎng)夏季最大電力負荷變化特征及其與氣象要素的關系
    自然災害、農(nóng)戶傳統(tǒng)風險管理手段與天氣指數(shù)保險需求
    吉林省松花江流域產(chǎn)業(yè)生態(tài)系統(tǒng)適應性分析
    沈陽市1951—2013年氣候變化特征及其區(qū)域蒸發(fā)的響應分析
    基于一維水動力模型的洪水頂托影響分析
    北京市朝陽區(qū)大氣污染物時空分布特征及與氣象要素的關系研究
    探測環(huán)境變化對臨沭站氣象要素的影響
    農(nóng)業(yè)面源氨氮污染減排體系的構建及減排控制對策
    Topmodel在布哈河流域徑流模擬中的應用
    探秘“大徑流”
    国产av不卡久久| ponron亚洲| 搡老岳熟女国产| av欧美777| 成人三级黄色视频| 伊人久久大香线蕉亚洲五| 久久久久免费精品人妻一区二区 | 三级毛片av免费| 叶爱在线成人免费视频播放| 国产精品亚洲一级av第二区| 欧美最黄视频在线播放免费| 欧美成狂野欧美在线观看| 亚洲av第一区精品v没综合| 欧美日韩亚洲综合一区二区三区_| 国产欧美日韩一区二区三| 悠悠久久av| 午夜激情av网站| 91九色精品人成在线观看| 欧美另类亚洲清纯唯美| 大香蕉久久成人网| 欧美激情 高清一区二区三区| 一区二区日韩欧美中文字幕| 99re在线观看精品视频| 日韩国内少妇激情av| 中出人妻视频一区二区| 男人的好看免费观看在线视频 | АⅤ资源中文在线天堂| 亚洲,欧美精品.| 亚洲欧美一区二区三区黑人| 好男人在线观看高清免费视频 | 国产高清videossex| 国产亚洲欧美98| 国产精品美女特级片免费视频播放器 | 国产又色又爽无遮挡免费看| 又紧又爽又黄一区二区| 亚洲av成人一区二区三| 精品一区二区三区av网在线观看| 免费观看精品视频网站| 欧美激情久久久久久爽电影| 一二三四在线观看免费中文在| 国产精品亚洲一级av第二区| 色在线成人网| 99re在线观看精品视频| 亚洲精品美女久久久久99蜜臀| or卡值多少钱| а√天堂www在线а√下载| 丝袜人妻中文字幕| 久久香蕉激情| xxxwww97欧美| 欧美性猛交╳xxx乱大交人| 国产v大片淫在线免费观看| 精品熟女少妇八av免费久了| 亚洲中文字幕日韩| 高清在线国产一区| 国产精品影院久久| 一级毛片精品| 黄色丝袜av网址大全| 久久精品91无色码中文字幕| 日韩欧美一区视频在线观看| 两性午夜刺激爽爽歪歪视频在线观看 | 日韩欧美 国产精品| 男男h啪啪无遮挡| 久9热在线精品视频| 大型av网站在线播放| 99re在线观看精品视频| 国产精品1区2区在线观看.| 免费观看精品视频网站| 成人18禁高潮啪啪吃奶动态图| 久久香蕉激情| 国产黄色小视频在线观看| 天堂√8在线中文| 国产又色又爽无遮挡免费看| 久久久国产欧美日韩av| 欧美在线黄色| 亚洲片人在线观看| 亚洲av电影在线进入| 国产亚洲av嫩草精品影院| 熟女少妇亚洲综合色aaa.| 久久精品亚洲精品国产色婷小说| 国产野战对白在线观看| 午夜福利免费观看在线| 99热6这里只有精品| 亚洲在线自拍视频| 黑人巨大精品欧美一区二区mp4| 欧美黑人精品巨大| 亚洲精华国产精华精| 欧美日韩亚洲国产一区二区在线观看| 国语自产精品视频在线第100页| 91成人精品电影| АⅤ资源中文在线天堂| 久久久久九九精品影院| www.精华液| 美女国产高潮福利片在线看| 色综合亚洲欧美另类图片| 免费搜索国产男女视频| 精品欧美一区二区三区在线| 黑人操中国人逼视频| 亚洲全国av大片| 变态另类丝袜制服| 午夜精品在线福利| 男女那种视频在线观看| 国产视频一区二区在线看| 久久婷婷人人爽人人干人人爱| 国产成人精品无人区| 欧美中文日本在线观看视频| 亚洲午夜精品一区,二区,三区| 丁香欧美五月| 国产精品久久久久久人妻精品电影| 91麻豆av在线| 国产精品1区2区在线观看.| 婷婷六月久久综合丁香| 中亚洲国语对白在线视频| 99热这里只有精品一区 | 美女 人体艺术 gogo| 国产激情欧美一区二区| www.自偷自拍.com| 色综合站精品国产| 后天国语完整版免费观看| 国产精品一区二区免费欧美| 中文亚洲av片在线观看爽| 久久天堂一区二区三区四区| 男人操女人黄网站| 久久久久久久久免费视频了| 国产高清视频在线播放一区| 天堂动漫精品| 中文字幕精品亚洲无线码一区 | 1024香蕉在线观看| 麻豆成人午夜福利视频| 精品国产国语对白av| 成人18禁高潮啪啪吃奶动态图| 亚洲中文字幕日韩| 日本精品一区二区三区蜜桃| bbb黄色大片| 怎么达到女性高潮| av中文乱码字幕在线| 亚洲狠狠婷婷综合久久图片| 色综合欧美亚洲国产小说| 成人国产综合亚洲| 亚洲第一青青草原| bbb黄色大片| 91九色精品人成在线观看| 亚洲真实伦在线观看| 听说在线观看完整版免费高清| 女人被狂操c到高潮| 熟妇人妻久久中文字幕3abv| 国产91精品成人一区二区三区| 日韩大码丰满熟妇| 亚洲国产看品久久| 少妇裸体淫交视频免费看高清 | 欧洲精品卡2卡3卡4卡5卡区| 一级a爱视频在线免费观看| 他把我摸到了高潮在线观看| 国产欧美日韩一区二区精品| 搞女人的毛片| 午夜免费观看网址| 99热6这里只有精品| 亚洲自偷自拍图片 自拍| 一级a爱片免费观看的视频| 两个人视频免费观看高清| 校园春色视频在线观看| 女同久久另类99精品国产91| 两个人免费观看高清视频| 俺也久久电影网| 国产精品 欧美亚洲| 国产亚洲精品久久久久久毛片| 国产精品,欧美在线| 亚洲色图 男人天堂 中文字幕| 国产99久久九九免费精品| 一边摸一边做爽爽视频免费| 亚洲精品在线美女| 国产精品久久电影中文字幕| 午夜免费激情av| 欧美日韩乱码在线| 亚洲一码二码三码区别大吗| 中文字幕另类日韩欧美亚洲嫩草| 国产成人精品无人区| 亚洲精品中文字幕一二三四区| 亚洲精品久久国产高清桃花| 精品卡一卡二卡四卡免费| 成人亚洲精品一区在线观看| 香蕉丝袜av| 久久久国产欧美日韩av| 精品少妇一区二区三区视频日本电影| 99re在线观看精品视频| 亚洲,欧美精品.| 啦啦啦 在线观看视频| 99热只有精品国产| 久久精品国产清高在天天线| 国产精品久久久久久精品电影 | 精品久久久久久久毛片微露脸| 久热这里只有精品99| 制服人妻中文乱码| 欧美一级a爱片免费观看看 | 久久天躁狠狠躁夜夜2o2o| tocl精华| 无人区码免费观看不卡| e午夜精品久久久久久久| 人人妻人人看人人澡| 国产精品亚洲美女久久久| 少妇 在线观看| 日韩免费av在线播放| 亚洲专区中文字幕在线| 美女大奶头视频| 黄色视频,在线免费观看| 法律面前人人平等表现在哪些方面| 亚洲电影在线观看av| 成人免费观看视频高清| 亚洲男人的天堂狠狠| 侵犯人妻中文字幕一二三四区| 亚洲熟妇中文字幕五十中出| 久久草成人影院| 一边摸一边做爽爽视频免费| 日本免费a在线| 久久国产精品人妻蜜桃| 精品人妻1区二区| 丝袜人妻中文字幕| 亚洲午夜理论影院| 国产精品二区激情视频| 成熟少妇高潮喷水视频| 国产高清视频在线播放一区| 亚洲一码二码三码区别大吗| 亚洲av中文字字幕乱码综合 | 国产人伦9x9x在线观看| 国产野战对白在线观看| 免费在线观看黄色视频的| 不卡av一区二区三区| 日韩欧美三级三区| 老司机福利观看| 日韩中文字幕欧美一区二区| 久久精品国产亚洲av高清一级| 高潮久久久久久久久久久不卡| 亚洲自拍偷在线| 18禁裸乳无遮挡免费网站照片 | 国产精品香港三级国产av潘金莲| 女警被强在线播放| 国产亚洲精品久久久久久毛片| 黑丝袜美女国产一区| 91大片在线观看| 久久国产乱子伦精品免费另类| www.www免费av| 91在线观看av| 亚洲精品美女久久av网站| 亚洲av电影不卡..在线观看| www.自偷自拍.com| 后天国语完整版免费观看| 欧美性长视频在线观看| 哪里可以看免费的av片| 精品久久久久久久久久久久久 | 国产精华一区二区三区| 99国产精品一区二区三区| 亚洲精品色激情综合| 久久精品国产99精品国产亚洲性色| 夜夜看夜夜爽夜夜摸| 亚洲性夜色夜夜综合| 久久午夜综合久久蜜桃| 亚洲国产看品久久| 久久久国产成人精品二区| 老司机靠b影院| 欧美黑人巨大hd| 一级毛片高清免费大全| 怎么达到女性高潮| 午夜免费观看网址| 国产激情久久老熟女| 精品电影一区二区在线| 嫁个100分男人电影在线观看| 国产精品电影一区二区三区| 国产精品二区激情视频| 国产aⅴ精品一区二区三区波| 亚洲真实伦在线观看| 两性夫妻黄色片| 国产精品国产高清国产av| 国产一区二区在线av高清观看| 国产精品 欧美亚洲| 2021天堂中文幕一二区在线观 | 亚洲三区欧美一区| 少妇裸体淫交视频免费看高清 | 亚洲av成人一区二区三| 亚洲精品美女久久av网站| 亚洲国产中文字幕在线视频| av中文乱码字幕在线| 男人操女人黄网站| 在线观看www视频免费| 日日干狠狠操夜夜爽| 可以在线观看毛片的网站| 亚洲片人在线观看| 91麻豆av在线| 搡老岳熟女国产| 热99re8久久精品国产| 精品福利观看| 亚洲七黄色美女视频| 一区二区三区国产精品乱码| 国产黄a三级三级三级人| 丰满人妻熟妇乱又伦精品不卡| 最好的美女福利视频网| 国产日本99.免费观看| 成人特级黄色片久久久久久久| 最近最新中文字幕大全免费视频| 国产成+人综合+亚洲专区| 久久久水蜜桃国产精品网| 亚洲天堂国产精品一区在线| 欧美黑人精品巨大| 国产激情久久老熟女| 精品久久久久久,| 国产真实乱freesex| 不卡av一区二区三区| 久久国产乱子伦精品免费另类| 成年版毛片免费区| 成人国产一区最新在线观看| 国产av一区二区精品久久| 亚洲精品中文字幕一二三四区| 久久精品国产99精品国产亚洲性色| 亚洲一区二区三区不卡视频| 波多野结衣av一区二区av| 国产黄a三级三级三级人| 美女 人体艺术 gogo| 窝窝影院91人妻| 侵犯人妻中文字幕一二三四区| 国产欧美日韩一区二区三| 最近最新中文字幕大全免费视频| 久久精品国产99精品国产亚洲性色| 黄色成人免费大全| 欧美黑人精品巨大| 青草久久国产| 一级毛片女人18水好多| 丝袜在线中文字幕| 一夜夜www| 99国产精品99久久久久| 国产黄a三级三级三级人| 亚洲人成网站高清观看| 美女高潮喷水抽搐中文字幕| 久久国产精品人妻蜜桃| 深夜精品福利| 人人澡人人妻人| 久久精品国产亚洲av香蕉五月| 国产成人欧美在线观看| 亚洲精品在线观看二区| 非洲黑人性xxxx精品又粗又长| 国产一级毛片七仙女欲春2 | 757午夜福利合集在线观看| 精品久久久久久久毛片微露脸| 国产精品野战在线观看| 久久香蕉国产精品| 男男h啪啪无遮挡| 美女高潮到喷水免费观看| 怎么达到女性高潮| 美女高潮喷水抽搐中文字幕| 高潮久久久久久久久久久不卡| 欧美黄色片欧美黄色片| 亚洲人成网站高清观看| 亚洲精品美女久久久久99蜜臀| 日韩欧美国产在线观看| 久久久久久久久久黄片| 精品第一国产精品| 亚洲国产欧洲综合997久久, | 亚洲无线在线观看| 精品高清国产在线一区| 欧美性猛交黑人性爽| 国产精品电影一区二区三区| 久久久久久久午夜电影| 中文字幕高清在线视频| 91成年电影在线观看| 脱女人内裤的视频| 久久久国产精品麻豆| 国产精品野战在线观看| 亚洲人成77777在线视频| 窝窝影院91人妻| 19禁男女啪啪无遮挡网站| 啦啦啦观看免费观看视频高清| 国产私拍福利视频在线观看| 淫妇啪啪啪对白视频| 天堂动漫精品| 久久 成人 亚洲| 亚洲精品美女久久久久99蜜臀| 国产精品一区二区免费欧美| 亚洲熟妇熟女久久| 日日夜夜操网爽| 99精品欧美一区二区三区四区| 国产aⅴ精品一区二区三区波| 亚洲成人国产一区在线观看| 精品一区二区三区av网在线观看| 久久久水蜜桃国产精品网| 可以在线观看毛片的网站| 亚洲精品国产精品久久久不卡| 日韩大尺度精品在线看网址| 亚洲无线在线观看| 又黄又粗又硬又大视频| av视频在线观看入口| 日本熟妇午夜| 色播亚洲综合网| 亚洲精品一卡2卡三卡4卡5卡| 午夜福利在线观看吧| 在线观看免费视频日本深夜| 在线观看免费日韩欧美大片| 每晚都被弄得嗷嗷叫到高潮| 黄色片一级片一级黄色片| 18禁观看日本| 日韩一卡2卡3卡4卡2021年| 国产精品影院久久| 精品欧美一区二区三区在线| 国产精品久久久人人做人人爽| 久久午夜综合久久蜜桃| 久久久水蜜桃国产精品网| 麻豆一二三区av精品| 国产精品爽爽va在线观看网站 | 亚洲av美国av| 亚洲 欧美一区二区三区| 波多野结衣高清作品| 国内精品久久久久精免费| 一二三四社区在线视频社区8| 黑人欧美特级aaaaaa片| 99在线人妻在线中文字幕| 久久国产亚洲av麻豆专区| 90打野战视频偷拍视频| 在线观看舔阴道视频| 女性生殖器流出的白浆| 亚洲国产精品久久男人天堂| 日韩中文字幕欧美一区二区| 国产精品爽爽va在线观看网站 | 国产亚洲精品久久久久久毛片| 免费女性裸体啪啪无遮挡网站| 中文字幕久久专区| 欧美av亚洲av综合av国产av| 成人午夜高清在线视频 | 久久狼人影院| 免费在线观看黄色视频的| 国产精品乱码一区二三区的特点| 波多野结衣高清作品| 怎么达到女性高潮| 亚洲中文av在线| 国产激情偷乱视频一区二区| 亚洲狠狠婷婷综合久久图片| www日本在线高清视频| 欧美性猛交╳xxx乱大交人| 日日摸夜夜添夜夜添小说| 欧美性猛交黑人性爽| 女人爽到高潮嗷嗷叫在线视频| 亚洲专区中文字幕在线| 亚洲国产精品成人综合色| 精品福利观看| 后天国语完整版免费观看| 亚洲 国产 在线| 欧美色欧美亚洲另类二区| 亚洲成人精品中文字幕电影| 久久伊人香网站| 久久中文看片网| 亚洲欧美日韩高清在线视频| a级毛片a级免费在线| 国产亚洲精品一区二区www| 久久久久久久精品吃奶| 亚洲精品在线美女| 免费在线观看黄色视频的| 午夜影院日韩av| 精品久久久久久成人av| 日本精品一区二区三区蜜桃| 精品国产国语对白av| 最近最新中文字幕大全免费视频| 欧美成人午夜精品| 国产真实乱freesex| 2021天堂中文幕一二区在线观 | 宅男免费午夜| 一区二区三区国产精品乱码| 亚洲av熟女| 久热爱精品视频在线9| 后天国语完整版免费观看| 久久久久亚洲av毛片大全| 欧美成人免费av一区二区三区| 国产精品久久久久久精品电影 | 国产黄色小视频在线观看| 国产黄片美女视频| 欧美人与性动交α欧美精品济南到| 午夜免费成人在线视频| 亚洲av电影不卡..在线观看| 男女床上黄色一级片免费看| 在线国产一区二区在线| 最近最新中文字幕大全免费视频| 欧美三级亚洲精品| av片东京热男人的天堂| 少妇被粗大的猛进出69影院| 国产激情偷乱视频一区二区| 亚洲av熟女| 亚洲电影在线观看av| 中文字幕精品亚洲无线码一区 | 非洲黑人性xxxx精品又粗又长| 亚洲第一av免费看| 黄片播放在线免费| 叶爱在线成人免费视频播放| 又紧又爽又黄一区二区| 国产真实乱freesex| 欧美av亚洲av综合av国产av| 观看免费一级毛片| 国产精品自产拍在线观看55亚洲| 亚洲av五月六月丁香网| 国产av又大| 村上凉子中文字幕在线| 丰满人妻熟妇乱又伦精品不卡| 啦啦啦 在线观看视频| 搞女人的毛片| 国产精华一区二区三区| 欧美日本亚洲视频在线播放| 日韩欧美国产在线观看| 亚洲人成77777在线视频| 欧美日韩精品网址| 最近在线观看免费完整版| 国产男靠女视频免费网站| 99热这里只有精品一区 | 亚洲熟女毛片儿| 午夜老司机福利片| 男人舔女人下体高潮全视频| 50天的宝宝边吃奶边哭怎么回事| 成在线人永久免费视频| 又紧又爽又黄一区二区| 亚洲人成伊人成综合网2020| 久久99热这里只有精品18| 精品国内亚洲2022精品成人| 久久精品国产亚洲av高清一级| 动漫黄色视频在线观看| 成人国语在线视频| 亚洲成人国产一区在线观看| 91九色精品人成在线观看| 麻豆av在线久日| 午夜视频精品福利| 91麻豆av在线| 亚洲熟妇熟女久久| 美女国产高潮福利片在线看| 99精品久久久久人妻精品| 亚洲国产精品sss在线观看| 成人18禁高潮啪啪吃奶动态图| 免费高清在线观看日韩| 国产精品亚洲一级av第二区| 非洲黑人性xxxx精品又粗又长| 草草在线视频免费看| 真人做人爱边吃奶动态| 精品久久久久久久末码| 国产精品影院久久| 中文字幕最新亚洲高清| or卡值多少钱| 一级a爱视频在线免费观看| 久99久视频精品免费| 美女免费视频网站| 欧美最黄视频在线播放免费| cao死你这个sao货| 男人舔女人的私密视频| 成人18禁在线播放| 亚洲人成77777在线视频| 日本成人三级电影网站| 亚洲成a人片在线一区二区| 国产免费av片在线观看野外av| 一区二区三区国产精品乱码| 久久精品91无色码中文字幕| 午夜福利在线在线| 国产片内射在线| 国产精品自产拍在线观看55亚洲| 熟女电影av网| 亚洲成a人片在线一区二区| 国产精品免费一区二区三区在线| 亚洲专区字幕在线| 后天国语完整版免费观看| www日本黄色视频网| 脱女人内裤的视频| 母亲3免费完整高清在线观看| 大型黄色视频在线免费观看| 久久久久精品国产欧美久久久| 亚洲精品国产精品久久久不卡| 老熟妇乱子伦视频在线观看| 久99久视频精品免费| 婷婷六月久久综合丁香| 亚洲午夜理论影院| 久久性视频一级片| 国产一区在线观看成人免费| 免费一级毛片在线播放高清视频| 久久天躁狠狠躁夜夜2o2o| 国产男靠女视频免费网站| 两个人看的免费小视频| 久久中文字幕人妻熟女| 亚洲 欧美 日韩 在线 免费| 欧美中文日本在线观看视频| 哪里可以看免费的av片| 国产精品98久久久久久宅男小说| 免费在线观看亚洲国产| 日韩欧美三级三区| 久久久久国产精品人妻aⅴ院| 国产极品粉嫩免费观看在线| 国产高清激情床上av| 黑人操中国人逼视频| 青草久久国产| 男人舔女人下体高潮全视频| 看黄色毛片网站| 久久精品成人免费网站| 国产野战对白在线观看| 看免费av毛片| 最近最新免费中文字幕在线| 午夜福利在线在线| 我的亚洲天堂| 国产成人精品久久二区二区91| 久久精品人妻少妇| 日本在线视频免费播放| 精品国产乱码久久久久久男人| 老汉色∧v一级毛片| 成人亚洲精品av一区二区| a级毛片a级免费在线| 精品国产超薄肉色丝袜足j| 欧美成人免费av一区二区三区| 欧美日韩中文字幕国产精品一区二区三区| 亚洲欧美精品综合一区二区三区| 免费在线观看亚洲国产| 亚洲久久久国产精品| 好看av亚洲va欧美ⅴa在| 久久久久九九精品影院| 亚洲第一av免费看| 窝窝影院91人妻| 久久精品成人免费网站| 一级a爱片免费观看的视频| 中文字幕高清在线视频| 成人手机av| 午夜福利欧美成人|