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

    松花江流域極端氣象水文要素分析

    2021-05-10 00:50:24章曉夢(mèng)蕭峻瓊次旦多杰
    關(guān)鍵詞:松花江流域洪量嫩江

    章曉夢(mèng),郝 潔,鞠 琴*,蕭峻瓊,次旦多杰,陳 璽

    (1. 河海大學(xué)水文水資源學(xué)院,江蘇 南京 210098;2. 水利部交通運(yùn)輸部國(guó)家能源局南京水利科學(xué)研究院,江蘇 南京 210000;3. 西藏水文水資源勘測(cè)局林芝水文水資源分局, 西藏 林芝 860000;4. 長(zhǎng)江水利委員會(huì)水文局,湖北 武漢 430010)

    松花江流域內(nèi)各地區(qū)有顯著的氣候差異,且降水空間分布不均勻[1]。任國(guó)玉等[2]的研究指出,東北偏高緯度地區(qū)的降水量呈增加或者不明顯變化趨勢(shì)。曹小磊等[3]研究結(jié)果表明氣溫呈顯著升高趨勢(shì),降水量和徑流量呈不顯著下降趨勢(shì),影響天然徑流的主要?dú)庀笠厥墙邓?。袁杰等[4]研究結(jié)果表明松花江各流域年徑流量均有下降的趨勢(shì)。已有學(xué)者針對(duì)極端氣候事件定義方法、極端降水事件規(guī)律特點(diǎn)進(jìn)行研究[5-8]。目前針對(duì)松花江流域氣候變化的研究大部分集中在單個(gè)氣象、水文要素,如降水、氣溫或徑流的研究,缺少對(duì)極端氣象水文要素變化的分析以及其成因的研究。

    本文擬以1961—2018年間松花江流域39個(gè)氣象站點(diǎn)的逐日降水觀測(cè)資料,以及1961—2018年間松花江流域佳木斯、大賚、扶余3個(gè)水文站點(diǎn)的日徑流量資料為研究對(duì)象,采用趨勢(shì)分析法、Mann-kendall非參數(shù)檢驗(yàn)方法、Mann-kendall突變檢驗(yàn)法、克里金插值法、皮爾遜相關(guān)系數(shù)等方法進(jìn)行松花江流域極端降水和極端徑流的時(shí)空變化規(guī)律和變異點(diǎn)診斷,分析松花江流域及各子流域極端氣象水文特征以及極端降水與極端徑流事件的關(guān)聯(lián)性,為該流域極端氣象水文事件的預(yù)測(cè)及防災(zāi)減災(zāi)規(guī)劃提供科學(xué)依據(jù)。

    1 資料與方法

    1.1 研究區(qū)域概況及數(shù)據(jù)來(lái)源

    松花江流域位于中國(guó)東北的北部地區(qū),其經(jīng)緯度范圍為北緯41°42′至51°38′、東經(jīng)119°52′至132°31′,流域內(nèi)有內(nèi)蒙古、黑龍江、吉林三個(gè)省份,總流域面積為55.68×104km2。將松花江流域分為松花江干流流域、二松流域、嫩江流域三個(gè)子流域進(jìn)行分析,選擇流域內(nèi)39個(gè)氣象站點(diǎn)及3個(gè)具有代表性的水文站點(diǎn),3個(gè)水文站分別為嫩江的大賚站、第二松花江的扶余站以及松花江干流的佳木斯站。松花江流域站點(diǎn)分布圖見(jiàn)圖1。

    注:Ⅰ:松花江干流流域,Ⅱ:第二松花江流域,Ⅲ:嫩江流域。圖1 松花江流域站點(diǎn)分布圖Fig.1 Site distribution map in Songhua River Basin

    所用降水?dāng)?shù)據(jù)為中國(guó)氣象數(shù)據(jù)網(wǎng)提供的1961—2018年的逐日降水實(shí)測(cè)資料。對(duì)數(shù)據(jù)資料進(jìn)行篩選和質(zhì)控處理,并對(duì)存在降水?dāng)?shù)據(jù)缺失的日期進(jìn)行了剔除,最終提取出所選39個(gè)站點(diǎn)的相關(guān)信息,得到了完整連續(xù)的逐日降水實(shí)測(cè)數(shù)據(jù)。徑流實(shí)測(cè)數(shù)據(jù)為佳木斯、扶余、大賚水文站1961年1月1日—2018年12月31日的逐日徑流數(shù)據(jù)。

    1.2 研究方法

    世界氣象組織氣候委員會(huì)等組織成立了氣候變化監(jiān)測(cè)和指標(biāo)專家組,并定義了27個(gè)典型氣候指數(shù),包括16個(gè)氣溫指數(shù)和11個(gè)降水指數(shù)[9],目前這些指數(shù)已被廣泛應(yīng)用于極端氣候研究。本文選取其中11個(gè)降水指數(shù)(表1)進(jìn)行研究。極端徑流事件采用洪量作為標(biāo)準(zhǔn),用最大1日洪量、最大3日洪量、最大7日洪量、最大15日洪量、最大30日洪量作為標(biāo)準(zhǔn)[10](表2)。

    運(yùn)用趨勢(shì)分析法、Mann-kendall非參數(shù)檢驗(yàn)方法[11]對(duì)極端降水指數(shù)和極端徑流指數(shù)進(jìn)行趨勢(shì)分析,采用克里金插值繪制極端降水指數(shù)的空間分布圖,采用Mann-kendall突變檢驗(yàn)法[11]得到各極端指數(shù)的突變點(diǎn),利用皮爾遜相關(guān)系數(shù)分析極端降水與極端徑流的關(guān)聯(lián)性。

    2 結(jié)果分析

    2.1 極端降水要素

    2.1.1 極端降水要素的空間分布特征

    松花江流域的極端降水指數(shù)空間分布特征見(jiàn)圖2。強(qiáng)度指數(shù)中,降水強(qiáng)度(SDII)分布較均勻,最大降水強(qiáng)度9.14 mm/d出現(xiàn)在長(zhǎng)春,最小降水強(qiáng)度7.49 mm/d出現(xiàn)在富錦,其余強(qiáng)度指數(shù)在松花江干流流域由西向東先增加再減小,第二松花江流域由西南向東北部逐漸減小,嫩江子流域由西南向東北呈先增加后減小趨勢(shì),Rx1day、Rx5day這兩個(gè)強(qiáng)度指數(shù)的最大值均出現(xiàn)在流域最南部,分別為67.99 mm和107.69 mm。

    表1 極端降水指數(shù)定義

    表2 極端徑流指數(shù)定義

    絕對(duì)指數(shù)中,在嫩江子流域,R1和R10指數(shù)由東南向西北呈增加趨勢(shì),R20指數(shù)由東南向西北先增大后減?。辉诘诙苫ń退苫ń闪鞲浇?,三個(gè)絕對(duì)指數(shù)的分布特征相同,呈中南部大,西部、東部小的分布特點(diǎn)。

    持續(xù)指數(shù)中,持續(xù)干旱指數(shù)和持續(xù)濕潤(rùn)指數(shù)是長(zhǎng)時(shí)間內(nèi)持續(xù)降雨情況的反映,其中持續(xù)干旱指數(shù)反映長(zhǎng)時(shí)間內(nèi)無(wú)降雨的情況,是研究干旱的方法之一[12]。持續(xù)干期(CDD)與持續(xù)濕期(CWD)、年總降水量(PRCPTOT)呈完全相反的空間分布,CDD指數(shù)在嫩江子流域呈南部大、北部小的趨勢(shì),在二松子流域及松花江干流子流域呈南部小、北部大的分布趨勢(shì),最大值出現(xiàn)在嫩江流域南部的富錦,約達(dá)75.07 d,最小值出現(xiàn)在第二松花江流域南部的東崗,為27.84 d。

    相對(duì)指數(shù)在嫩江流域從東南到西北呈增大趨勢(shì),在松花江干流和第二松花江附近從西到東先增大后減小,從南到北先減小后增大。

    極端降水指數(shù)存在區(qū)域性差異,除CDD和SDII外,其余指數(shù)的分布特征均呈現(xiàn)出第二松花江流域?yàn)楦咧祬^(qū),松花江干流流域次之,嫩江流域?yàn)榈椭祬^(qū)。說(shuō)明在松花江流域,嫩江流域降水量、降水日數(shù)、降水強(qiáng)度較小,第二松花江流域極端降水事件強(qiáng)度高且持續(xù)性強(qiáng),應(yīng)成為重點(diǎn)防護(hù)對(duì)象。

    2.1.2 極端降水要素的空間趨勢(shì)分布

    極端降水要素在不同空間的趨勢(shì)性反映了其變化趨勢(shì)的總體規(guī)律,可借助M-K趨勢(shì)檢驗(yàn)對(duì)松花江流域各極端降水指標(biāo)的變化趨勢(shì)進(jìn)行檢驗(yàn),其空間分布情況如圖3所示。松花江流域極端降水事件表現(xiàn)出明顯的區(qū)域差異性。

    強(qiáng)度指數(shù)是一個(gè)衡量是否發(fā)生洪澇災(zāi)害的重要指標(biāo),23個(gè)站點(diǎn)的Rx1day指數(shù)表現(xiàn)出上升趨勢(shì),其中5個(gè)站點(diǎn)通過(guò)了α=0.05的顯著性檢驗(yàn);有21個(gè)站點(diǎn)的Rx5day指數(shù)表現(xiàn)出上升趨勢(shì),18個(gè)站點(diǎn)表現(xiàn)出下降趨勢(shì);25個(gè)站點(diǎn)的SDII指數(shù)表現(xiàn)出上升趨勢(shì),14個(gè)站點(diǎn)表現(xiàn)出下降趨勢(shì)。Rx1day、Rx5day、SDII三個(gè)強(qiáng)度指數(shù)顯著增加的站點(diǎn)大多分布在嫩江子流域南部和松花江干流子流域西部,且最大增幅均出現(xiàn)在明水站。

    圖2 極端降水指數(shù)空間分布Fig.2 Spatial distribution of extreme precipitation index

    絕對(duì)指數(shù)中,R10指數(shù)趨于增加站點(diǎn)數(shù)為27個(gè),僅有安達(dá)站變化趨勢(shì)通過(guò)了α=0.05的顯著性檢驗(yàn);R20指數(shù)趨于增加站點(diǎn)數(shù)為25個(gè),僅有安達(dá)站和白城站變化趨勢(shì)通過(guò)了α=0.05的顯著性檢驗(yàn);有29個(gè)站點(diǎn)的R1指數(shù)表現(xiàn)出上升趨勢(shì),其中5個(gè)站點(diǎn)呈顯著上升趨勢(shì)。與下降趨勢(shì)相比,絕對(duì)指數(shù)在松花江流域的上升趨勢(shì)更加明顯,但變化趨勢(shì)并不顯著,大多呈微弱上升趨勢(shì)。

    持續(xù)指數(shù)中,僅有5個(gè)站點(diǎn)的CDD指數(shù)呈微弱上升趨勢(shì),34個(gè)站點(diǎn)呈下降趨勢(shì),呈顯著下降趨勢(shì)站點(diǎn)主要分布在嫩江流域中部和松花江干流流域北部;有19個(gè)站點(diǎn)的CWD指數(shù)表現(xiàn)出上升趨勢(shì),20個(gè)站點(diǎn)表現(xiàn)出下降趨勢(shì);年總降水量(PRCPTOT)趨于增加的站點(diǎn)數(shù)為32個(gè),僅有7個(gè)站點(diǎn)PRCPTOT指數(shù)趨于減小,大多集中在嫩江流域南部,代表松花江流域整體降水量有所增加。

    相對(duì)指數(shù)中,19個(gè)站點(diǎn)的R95pTOT指數(shù)表現(xiàn)出上升趨勢(shì),其中4個(gè)站點(diǎn)通過(guò)了α=0.05的顯著性檢驗(yàn),22個(gè)站點(diǎn)的R99pTOT指數(shù)表現(xiàn)出上升趨勢(shì),其中4個(gè)站點(diǎn)呈顯著上升趨勢(shì)。由此可見(jiàn),同一個(gè)極端降水指標(biāo)在不同時(shí)間段和不同區(qū)域的變化趨勢(shì)均不同,體現(xiàn)了極端降水事件在時(shí)間和空間上的復(fù)雜性和差異性。

    2.1.3 極端降水要素的突變分析

    掌握徑流演變規(guī)律,除分析其總體趨勢(shì)外,還必須判斷并檢驗(yàn)突變發(fā)生的時(shí)間、次數(shù)以及變化幅度[13]。水文學(xué)或氣候?qū)W上的突變是指系統(tǒng)從一種穩(wěn)定狀態(tài)跳躍到另一種穩(wěn)定狀態(tài)的現(xiàn)象[14]。借助M-K突變檢驗(yàn)方法對(duì)松花江流域三個(gè)子流域1961—2018年的極端降水指數(shù)進(jìn)行突變分析,突變分析結(jié)果見(jiàn)表3。

    在松花江干流流域,強(qiáng)度指數(shù)和絕對(duì)指數(shù)無(wú)突變點(diǎn),持續(xù)指數(shù)CDD、CWD、PCRPTOT的突變年份分別為1998年、1972年、1963年,相對(duì)指數(shù)R95pTOT、R99pTOT的突變年份為1964和1966年左右;二松流域僅有持續(xù)干期CDD指數(shù)在1967年有明顯突變;嫩江流域的強(qiáng)度指數(shù)中的Rx1day、Rx5day指數(shù)分別在1970年、1963年發(fā)生突變;整個(gè)松花江流域各極端降水指數(shù)中僅有Rx5day和CDD指數(shù)分別在1964年和1993年左右開(kāi)始突變。松花江流域及各子流域在突變點(diǎn)后的指數(shù)序列均為突變減小的過(guò)程,各指數(shù)的突變時(shí)間表現(xiàn)出明顯的空間差異性。

    圖3 極端降水指數(shù)空間趨勢(shì)分布Fig.3 Spatial trend distribution of extreme precipitation index

    表3 極端降水指數(shù)突變檢驗(yàn)

    2.2 極端徑流要素

    不同時(shí)間尺度下年最大洪量的變化趨勢(shì)和突變年份相近,見(jiàn)表4,佳木斯站、扶余站、大賚站1961—2018年各極端徑流指數(shù)M-K檢驗(yàn)的臨界值|Z|均小于1.96,未通過(guò)α=0.05的顯著性檢驗(yàn),在佳木斯站和大賚站呈微弱下降趨勢(shì),在扶余站Rx1day、Rx3day指數(shù)呈微弱下降趨勢(shì),Rx7day、Rx15day、Rx30day指數(shù)呈微弱上升趨勢(shì),各指數(shù)總體變化趨勢(shì)較為穩(wěn)定。圖4以年最大7日洪量(Rx7day)為例分析各水文站變化趨勢(shì)和突變情況。

    三個(gè)水文控制站中佳木斯站和大賚站極端徑流量較大,扶余站的徑流量較小。從圖4可以看出,佳木斯站RX7day指數(shù)M-K檢驗(yàn)圖的UF和UB曲線在0.05顯著性水平閾值線內(nèi)有多個(gè)交點(diǎn),其中1967年為突變點(diǎn),在1976年左右UF曲線超過(guò)了置信曲線,指數(shù)值有一個(gè)突變減小的過(guò)程,其余交點(diǎn)可以看出RX7day有突變減小的趨勢(shì),但減小趨勢(shì)尚未達(dá)到0.05顯著性水平;指數(shù)值年際變化較大,在20世紀(jì)60年代至80年代初呈下降趨勢(shì),20世紀(jì)80年代至2000年左右呈上升趨勢(shì),21世紀(jì)初呈現(xiàn)減少再增加的趨勢(shì)。

    扶余站RX7day指數(shù)M-K檢驗(yàn)圖的UF和UB曲線在0.05顯著性水平閾值線內(nèi)有多個(gè)交點(diǎn),但UF曲線從未超過(guò)置信曲線,代表RX7day指數(shù)在扶余站無(wú)突變點(diǎn),指數(shù)值呈波動(dòng)性變化,總體變化趨勢(shì)不明顯。

    表4 極端徑流指數(shù)趨勢(shì)及突變

    圖4 最大7日降水量(Rx7day)突變分析Fig.4 Mutation test of Rx7day

    大賚站RX7day指數(shù)M-K檢驗(yàn)圖的UF和UB曲線在0.05顯著性水平閾值線內(nèi)有多個(gè)交點(diǎn),其中突變點(diǎn)為1964年,在1975年左右UF曲線超過(guò)了置信曲線,徑流量有一個(gè)突變減小的過(guò)程。在置信區(qū)間外UF曲線和UB曲線交于1983年,說(shuō)明在該點(diǎn)RX7day指數(shù)有增加的趨勢(shì),但變化趨勢(shì)沒(méi)有達(dá)到突變的水平,因此不是突變點(diǎn);指數(shù)值年代際變化幅度較大,20世紀(jì)50年代至80年代初、21世紀(jì)以后相對(duì)減少,20世紀(jì)80年代至2000年相對(duì)增加。

    嫩江流域和松花江全流域的強(qiáng)度指數(shù)均在20世紀(jì)60年代發(fā)生突變,其控制站大賚站和佳木斯站的年最大洪量也于20世紀(jì)60年代發(fā)生突變。在二松流域,1953年豐滿水庫(kù)開(kāi)始進(jìn)行人工調(diào)節(jié)之后,對(duì)大洪水洪峰的削減效果十分顯著,對(duì)防洪起到了很大作用,因此扶余站極端徑流指數(shù)的年代際變化不明顯,且無(wú)突變點(diǎn)。

    2.3 極端降水和極端徑流事件的關(guān)聯(lián)性

    根據(jù)極端徑流指數(shù)所定義的年最大1、3、7、15、30日洪量,取其中相對(duì)較長(zhǎng)時(shí)間尺度,分析松花江流域在時(shí)間尺度為7日、15日、30日時(shí)極端降水和徑流事件兩者間的關(guān)聯(lián)性,通過(guò)皮爾遜相關(guān)系數(shù)對(duì)松花江流域不同時(shí)間尺度下的年最大降水量與年最大洪量進(jìn)行相關(guān)性分析。由表5可知,在松花江流域,各極端降水指數(shù)和極端徑流指數(shù)的相關(guān)系數(shù)均通過(guò)99%置信度檢驗(yàn),表現(xiàn)出顯著的正相關(guān)關(guān)系,表明這兩者之間存在一定的關(guān)聯(lián)性,極端降水量的變化必然帶來(lái)極端徑流量的變化。在各控制站所對(duì)應(yīng)流域范圍內(nèi),年最大30天降水量與年最大30天洪量表現(xiàn)出最為顯著的關(guān)聯(lián)性,其相關(guān)系數(shù)分別為0.678、0.713、0.807。時(shí)間尺度為30天的極端降水和極端徑流的關(guān)聯(lián)性最為明顯,說(shuō)明松花江流域面積較大,支流較多,匯流時(shí)間較長(zhǎng),所以隨著時(shí)間尺度的增大,其極端降水與極端徑流的相關(guān)性更加顯著。

    表5 極端降水與極端徑流相關(guān)系數(shù)

    3 結(jié)論

    本文基于1961—2018年松花江流域39個(gè)氣象站點(diǎn)逐日降水資料和3個(gè)水文站點(diǎn)逐日流量資料,對(duì)極端降水要素和極端徑流要素的變化趨勢(shì)、時(shí)空分布特征及關(guān)聯(lián)性進(jìn)行了分析。

    1)松花江流域極端降水指數(shù)值變化規(guī)律呈現(xiàn)出區(qū)域差異性,除持續(xù)干期CDD和降水強(qiáng)度SDII外,其余指數(shù)的分布特征均呈現(xiàn)出第二松花江流域?yàn)楦咧祬^(qū),松花江干流流域次之,嫩江流域?yàn)榈椭祬^(qū)。各控制水文站極端徑流指數(shù)變化趨勢(shì)均未通過(guò)95%顯著性檢驗(yàn),且僅扶余站極端徑流指數(shù)未出現(xiàn)突變點(diǎn)。

    2)反映強(qiáng)降水量的極端降水指數(shù)在大部分地區(qū)表現(xiàn)為增加趨勢(shì),但極端徑流指數(shù)卻大多表現(xiàn)為減少趨勢(shì)或不顯著上升趨勢(shì),未顯示出同步變化趨勢(shì)。降水量與徑流量之間關(guān)系密切,但復(fù)雜的氣象活動(dòng)和大壩、水庫(kù)等水利工程的修建都是出現(xiàn)該現(xiàn)象的重要影響因素,例如二松流域上的豐滿水庫(kù)會(huì)使徑流過(guò)程逐漸改變,從而影響極端徑流事件的演變規(guī)律,導(dǎo)致在極端降水增加的情況下,極端徑流指數(shù)卻呈減少趨勢(shì)。

    3)不同時(shí)間尺度下各極端降水指數(shù)與極端徑流指數(shù)的相關(guān)系數(shù)均通過(guò)99%顯著性檢驗(yàn),由于松花江流域面積較大,支流較多,匯流時(shí)間較長(zhǎng),當(dāng)時(shí)間尺度為30天時(shí),極端降水與極端徑流表現(xiàn)出最為顯著的關(guān)聯(lián)性。

    猜你喜歡
    松花江流域洪量嫩江
    基于SPA 的北江流域峰量關(guān)系研究
    陜西水利(2023年12期)2023-12-19 03:28:32
    自然災(zāi)害、農(nóng)戶傳統(tǒng)風(fēng)險(xiǎn)管理手段與天氣指數(shù)保險(xiǎn)需求
    推進(jìn)嫩江市霍龍門鎮(zhèn)社會(huì)全面振興 做到“六個(gè)堅(jiān)持”實(shí)施“六大戰(zhàn)略”
    活力(2021年4期)2021-07-28 05:34:10
    百年嫩江
    吉林省松花江流域產(chǎn)業(yè)生態(tài)系統(tǒng)適應(yīng)性分析
    遼河干流主要控制站近75年最大洪峰及洪量變化特征分析研究
    基于一維水動(dòng)力模型的洪水頂托影響分析
    農(nóng)業(yè)面源氨氮污染減排體系的構(gòu)建及減排控制對(duì)策
    嫩江情
    梁喜才:嫩江水畔的領(lǐng)頭雁
    免费看日本二区| 日本撒尿小便嘘嘘汇集6| 777久久人妻少妇嫩草av网站| 99精品久久久久人妻精品| 日韩欧美在线乱码| 制服人妻中文乱码| 中文字幕久久专区| 欧美乱色亚洲激情| 国产精品一区二区免费欧美| 久久精品国产99精品国产亚洲性色| 中文字幕熟女人妻在线| 每晚都被弄得嗷嗷叫到高潮| 亚洲av中文字字幕乱码综合| 香蕉国产在线看| 91麻豆精品激情在线观看国产| 激情在线观看视频在线高清| 成年女人毛片免费观看观看9| 午夜亚洲福利在线播放| 久久久久亚洲av毛片大全| 国产精品 欧美亚洲| 91成年电影在线观看| 日韩免费av在线播放| 中文字幕av在线有码专区| 黄色片一级片一级黄色片| 91成年电影在线观看| 国产精品98久久久久久宅男小说| 欧洲精品卡2卡3卡4卡5卡区| 亚洲性夜色夜夜综合| 中文字幕人成人乱码亚洲影| 国产精品乱码一区二三区的特点| 久久精品影院6| 美女大奶头视频| 国产片内射在线| 午夜福利视频1000在线观看| АⅤ资源中文在线天堂| 国产精品影院久久| 国产成+人综合+亚洲专区| 桃色一区二区三区在线观看| 亚洲aⅴ乱码一区二区在线播放 | 啦啦啦免费观看视频1| 757午夜福利合集在线观看| 精品一区二区三区av网在线观看| 1024手机看黄色片| 国产高清视频在线观看网站| 日本免费一区二区三区高清不卡| 国产aⅴ精品一区二区三区波| 精品欧美国产一区二区三| 亚洲精品粉嫩美女一区| 老司机靠b影院| www日本在线高清视频| 性欧美人与动物交配| 国产精品久久久久久人妻精品电影| 黄色女人牲交| 午夜视频精品福利| 中文字幕av在线有码专区| 国产区一区二久久| 香蕉丝袜av| 97人妻精品一区二区三区麻豆| 亚洲国产日韩欧美精品在线观看 | 国内久久婷婷六月综合欲色啪| 真人一进一出gif抽搐免费| 午夜亚洲福利在线播放| 最好的美女福利视频网| 一本大道久久a久久精品| www.熟女人妻精品国产| 1024手机看黄色片| 18禁美女被吸乳视频| 男女做爰动态图高潮gif福利片| 999久久久国产精品视频| 亚洲av中文字字幕乱码综合| 91麻豆av在线| 91国产中文字幕| 亚洲av成人av| 波多野结衣高清无吗| 欧美日韩瑟瑟在线播放| 91大片在线观看| 高清毛片免费观看视频网站| 日韩欧美三级三区| av天堂在线播放| 一区二区三区国产精品乱码| 久久久久久亚洲精品国产蜜桃av| 最近最新免费中文字幕在线| 久久精品综合一区二区三区| 最好的美女福利视频网| 精品一区二区三区视频在线观看免费| 色综合婷婷激情| 欧美3d第一页| 国产精华一区二区三区| 日本免费a在线| 国产高清视频在线观看网站| 亚洲 国产 在线| 色老头精品视频在线观看| 国产99久久九九免费精品| 老司机午夜福利在线观看视频| 亚洲黑人精品在线| 亚洲狠狠婷婷综合久久图片| 美女免费视频网站| 变态另类成人亚洲欧美熟女| 成人精品一区二区免费| 欧美极品一区二区三区四区| 国产欧美日韩一区二区精品| av有码第一页| 国产亚洲精品第一综合不卡| 欧美激情久久久久久爽电影| 丰满人妻熟妇乱又伦精品不卡| 久久99热这里只有精品18| 亚洲国产高清在线一区二区三| 一区二区三区国产精品乱码| 日本一本二区三区精品| 真人一进一出gif抽搐免费| 一夜夜www| 午夜免费观看网址| av在线播放免费不卡| 国产精品一及| 亚洲国产欧洲综合997久久,| 五月玫瑰六月丁香| 亚洲av日韩精品久久久久久密| 美女黄网站色视频| 国产一区二区在线观看日韩 | 国产一级毛片七仙女欲春2| 日本成人三级电影网站| 老司机靠b影院| 婷婷精品国产亚洲av在线| 99久久久亚洲精品蜜臀av| 亚洲国产精品sss在线观看| 久久精品国产亚洲av香蕉五月| ponron亚洲| 午夜视频精品福利| 亚洲国产高清在线一区二区三| 国产亚洲精品一区二区www| 无限看片的www在线观看| 夜夜看夜夜爽夜夜摸| 好男人电影高清在线观看| 久久婷婷成人综合色麻豆| 国产免费男女视频| 一进一出抽搐gif免费好疼| 啦啦啦免费观看视频1| 少妇人妻一区二区三区视频| 国产野战对白在线观看| 人妻夜夜爽99麻豆av| 在线观看免费视频日本深夜| 757午夜福利合集在线观看| 女人爽到高潮嗷嗷叫在线视频| 天堂av国产一区二区熟女人妻 | 两个人看的免费小视频| 色av中文字幕| 夜夜夜夜夜久久久久| 女警被强在线播放| 久久性视频一级片| av视频在线观看入口| 国产伦人伦偷精品视频| 国产一区二区三区在线臀色熟女| 长腿黑丝高跟| 免费高清视频大片| 亚洲精品色激情综合| 亚洲成人中文字幕在线播放| 亚洲av成人不卡在线观看播放网| 精品一区二区三区视频在线观看免费| 亚洲成人中文字幕在线播放| 激情在线观看视频在线高清| 国产精品电影一区二区三区| 久久久久久久久免费视频了| 久久草成人影院| 最新在线观看一区二区三区| av福利片在线| 亚洲五月婷婷丁香| 黄片大片在线免费观看| 精品免费久久久久久久清纯| 欧美日韩亚洲综合一区二区三区_| 欧美黑人欧美精品刺激| 久久天堂一区二区三区四区| 丁香欧美五月| 日韩精品青青久久久久久| 日本黄色视频三级网站网址| 在线国产一区二区在线| 亚洲第一电影网av| 午夜成年电影在线免费观看| 99热6这里只有精品| 一本久久中文字幕| 久久天躁狠狠躁夜夜2o2o| 淫秽高清视频在线观看| 亚洲激情在线av| 在线永久观看黄色视频| 美女午夜性视频免费| 国产三级在线视频| 97碰自拍视频| 国产野战对白在线观看| 国内揄拍国产精品人妻在线| 男女做爰动态图高潮gif福利片| 国产av在哪里看| 亚洲全国av大片| 亚洲熟女毛片儿| 一二三四在线观看免费中文在| 老司机深夜福利视频在线观看| 好看av亚洲va欧美ⅴa在| 搡老熟女国产l中国老女人| 9191精品国产免费久久| 欧美高清成人免费视频www| 亚洲中文日韩欧美视频| 黄色毛片三级朝国网站| 91老司机精品| 高潮久久久久久久久久久不卡| 夜夜躁狠狠躁天天躁| 91大片在线观看| 国产aⅴ精品一区二区三区波| 人妻丰满熟妇av一区二区三区| 黄色a级毛片大全视频| 人人妻人人看人人澡| 日日夜夜操网爽| 国产99久久九九免费精品| 免费观看精品视频网站| 人妻丰满熟妇av一区二区三区| 中文字幕久久专区| 久久久久久久久中文| 亚洲一区二区三区不卡视频| 国产免费男女视频| 国产精品1区2区在线观看.| 在线观看免费视频日本深夜| av福利片在线观看| 亚洲熟女毛片儿| 国产高清视频在线播放一区| 午夜免费成人在线视频| 午夜福利18| 两个人免费观看高清视频| 亚洲成人国产一区在线观看| 亚洲国产欧洲综合997久久,| 久久久久久国产a免费观看| 日本成人三级电影网站| 丁香六月欧美| 9191精品国产免费久久| 欧美性长视频在线观看| 国产精品99久久99久久久不卡| 成人高潮视频无遮挡免费网站| 国产精品98久久久久久宅男小说| 日韩中文字幕欧美一区二区| 少妇粗大呻吟视频| 天天躁狠狠躁夜夜躁狠狠躁| 一a级毛片在线观看| 黄片小视频在线播放| www.999成人在线观看| 亚洲人成77777在线视频| 在线观看美女被高潮喷水网站 | 搡老熟女国产l中国老女人| 日韩精品中文字幕看吧| 麻豆成人av在线观看| 午夜福利高清视频| av国产免费在线观看| 88av欧美| 一个人免费在线观看的高清视频| 制服丝袜大香蕉在线| 亚洲五月婷婷丁香| 热99re8久久精品国产| 亚洲第一欧美日韩一区二区三区| 欧美中文综合在线视频| 国产精品电影一区二区三区| 久久精品国产综合久久久| www日本在线高清视频| 亚洲精品在线美女| 神马国产精品三级电影在线观看 | 精品久久久久久久末码| 亚洲av第一区精品v没综合| 十八禁人妻一区二区| 色av中文字幕| 真人做人爱边吃奶动态| 三级男女做爰猛烈吃奶摸视频| 久久久国产欧美日韩av| 床上黄色一级片| 成人手机av| 精品日产1卡2卡| 国产免费男女视频| 啦啦啦免费观看视频1| 两个人视频免费观看高清| 午夜福利成人在线免费观看| 欧美不卡视频在线免费观看 | 亚洲欧美激情综合另类| 亚洲av中文字字幕乱码综合| 一二三四在线观看免费中文在| 女人高潮潮喷娇喘18禁视频| 中文字幕熟女人妻在线| 狂野欧美白嫩少妇大欣赏| 日本三级黄在线观看| 国产精品国产高清国产av| 伊人久久大香线蕉亚洲五| 亚洲精品久久成人aⅴ小说| 最近最新中文字幕大全免费视频| 国产亚洲欧美98| 黑人欧美特级aaaaaa片| 亚洲午夜精品一区,二区,三区| 日本五十路高清| 国产精品免费一区二区三区在线| 一级作爱视频免费观看| www日本在线高清视频| 亚洲精品美女久久久久99蜜臀| 免费在线观看成人毛片| 国产精品一区二区三区四区免费观看 | 欧美成人性av电影在线观看| av有码第一页| 少妇的丰满在线观看| 亚洲 欧美一区二区三区| 国内少妇人妻偷人精品xxx网站 | 麻豆国产97在线/欧美 | 一区二区三区国产精品乱码| 久久久久亚洲av毛片大全| 香蕉丝袜av| 免费看日本二区| 日韩中文字幕欧美一区二区| 国产精品爽爽va在线观看网站| 亚洲成a人片在线一区二区| 亚洲专区国产一区二区| 色综合站精品国产| 欧洲精品卡2卡3卡4卡5卡区| 少妇粗大呻吟视频| 两个人视频免费观看高清| 人妻丰满熟妇av一区二区三区| 国产成人系列免费观看| 后天国语完整版免费观看| 精品国产乱子伦一区二区三区| 午夜福利18| 丝袜人妻中文字幕| 成人18禁高潮啪啪吃奶动态图| 国产一区二区三区视频了| 国产高清视频在线观看网站| 久久这里只有精品19| 亚洲午夜理论影院| 国产不卡一卡二| 色播亚洲综合网| www.熟女人妻精品国产| 亚洲一区中文字幕在线| 黄片小视频在线播放| 高清毛片免费观看视频网站| 亚洲精品国产一区二区精华液| 欧美日韩精品网址| 一区福利在线观看| 99riav亚洲国产免费| 成人18禁高潮啪啪吃奶动态图| 我的老师免费观看完整版| 91成年电影在线观看| 无人区码免费观看不卡| 丁香六月欧美| 国产黄片美女视频| 亚洲精品在线观看二区| 亚洲国产欧美人成| 女警被强在线播放| 黑人操中国人逼视频| 久久人人精品亚洲av| 精品国产乱码久久久久久男人| 首页视频小说图片口味搜索| 亚洲欧美日韩东京热| 两个人免费观看高清视频| 日韩大码丰满熟妇| 高潮久久久久久久久久久不卡| 老熟妇仑乱视频hdxx| 成在线人永久免费视频| 可以免费在线观看a视频的电影网站| 色老头精品视频在线观看| 亚洲av美国av| 一本一本综合久久| 国产探花在线观看一区二区| 亚洲专区中文字幕在线| 国产黄a三级三级三级人| 操出白浆在线播放| 中文资源天堂在线| 巨乳人妻的诱惑在线观看| 啪啪无遮挡十八禁网站| 日韩欧美国产在线观看| 妹子高潮喷水视频| 国产亚洲欧美在线一区二区| x7x7x7水蜜桃| 一本综合久久免费| avwww免费| 国模一区二区三区四区视频 | 午夜成年电影在线免费观看| 三级男女做爰猛烈吃奶摸视频| 国产精品一区二区三区四区久久| 日韩有码中文字幕| 女人爽到高潮嗷嗷叫在线视频| 欧美在线黄色| 一边摸一边抽搐一进一小说| 欧美一级a爱片免费观看看 | 国产精华一区二区三区| 成人欧美大片| 女人高潮潮喷娇喘18禁视频| 国产高清有码在线观看视频 | 少妇熟女aⅴ在线视频| 午夜影院日韩av| 国产亚洲精品久久久久5区| 成熟少妇高潮喷水视频| 欧美不卡视频在线免费观看 | 日本撒尿小便嘘嘘汇集6| 免费一级毛片在线播放高清视频| 每晚都被弄得嗷嗷叫到高潮| 亚洲精品美女久久av网站| 午夜激情福利司机影院| 国产1区2区3区精品| 日本五十路高清| 一个人免费在线观看电影 | bbb黄色大片| 免费在线观看日本一区| 妹子高潮喷水视频| 日韩欧美三级三区| 别揉我奶头~嗯~啊~动态视频| 精品国内亚洲2022精品成人| 日韩欧美 国产精品| 亚洲狠狠婷婷综合久久图片| 首页视频小说图片口味搜索| 在线十欧美十亚洲十日本专区| 五月伊人婷婷丁香| 久久精品91无色码中文字幕| 91老司机精品| 91字幕亚洲| 搡老妇女老女人老熟妇| 欧洲精品卡2卡3卡4卡5卡区| 99国产精品一区二区三区| 在线视频色国产色| а√天堂www在线а√下载| 亚洲第一电影网av| 久久久久久人人人人人| 露出奶头的视频| 久久这里只有精品19| 亚洲专区字幕在线| 日韩三级视频一区二区三区| 一本精品99久久精品77| 精品人妻1区二区| 色播亚洲综合网| 中文字幕久久专区| 少妇被粗大的猛进出69影院| 99热这里只有精品一区 | 99热只有精品国产| e午夜精品久久久久久久| 婷婷精品国产亚洲av| 亚洲中文av在线| 欧美极品一区二区三区四区| 欧美一区二区精品小视频在线| cao死你这个sao货| 一区二区三区国产精品乱码| 婷婷六月久久综合丁香| 1024手机看黄色片| 白带黄色成豆腐渣| 99在线人妻在线中文字幕| 性欧美人与动物交配| 一区二区三区激情视频| 一个人免费在线观看电影 | 亚洲精品国产一区二区精华液| 亚洲电影在线观看av| 国产一区二区三区在线臀色熟女| 国产又黄又爽又无遮挡在线| 国产亚洲精品第一综合不卡| 一本久久中文字幕| 国产午夜精品论理片| 亚洲成人国产一区在线观看| 成人国语在线视频| 999精品在线视频| 男女之事视频高清在线观看| 亚洲av成人精品一区久久| 无遮挡黄片免费观看| 国产欧美日韩一区二区精品| 亚洲欧美日韩东京热| 色综合站精品国产| 日韩免费av在线播放| 欧美成人性av电影在线观看| 亚洲av中文字字幕乱码综合| 国产精品香港三级国产av潘金莲| 日本一区二区免费在线视频| 精品少妇一区二区三区视频日本电影| 久久人妻福利社区极品人妻图片| 成人一区二区视频在线观看| 午夜影院日韩av| 午夜精品在线福利| 成人国产综合亚洲| 日韩三级视频一区二区三区| 熟女少妇亚洲综合色aaa.| 日本黄大片高清| 欧美日韩乱码在线| 黑人欧美特级aaaaaa片| 久久国产乱子伦精品免费另类| 青草久久国产| 色老头精品视频在线观看| 午夜福利免费观看在线| 黑人欧美特级aaaaaa片| 一区福利在线观看| 国产成人av教育| 亚洲国产精品合色在线| 久久精品影院6| 婷婷丁香在线五月| 一区二区三区高清视频在线| 免费人成视频x8x8入口观看| 天堂av国产一区二区熟女人妻 | 男女之事视频高清在线观看| 窝窝影院91人妻| 少妇被粗大的猛进出69影院| aaaaa片日本免费| 在线a可以看的网站| 黄色视频不卡| 国内久久婷婷六月综合欲色啪| 男女午夜视频在线观看| 好男人在线观看高清免费视频| av免费在线观看网站| 露出奶头的视频| 免费在线观看日本一区| 色综合亚洲欧美另类图片| 久久草成人影院| 波多野结衣高清无吗| 狂野欧美白嫩少妇大欣赏| 每晚都被弄得嗷嗷叫到高潮| 中亚洲国语对白在线视频| 看黄色毛片网站| 国产91精品成人一区二区三区| 久久久久国产一级毛片高清牌| 午夜福利高清视频| 在线观看免费日韩欧美大片| 曰老女人黄片| 制服丝袜大香蕉在线| 国产成人精品久久二区二区免费| 国产高清视频在线观看网站| 99热这里只有精品一区 | 免费在线观看视频国产中文字幕亚洲| 亚洲av成人av| 免费观看精品视频网站| 国产精品免费一区二区三区在线| 国产激情久久老熟女| 人人妻人人澡欧美一区二区| 麻豆国产av国片精品| 熟女电影av网| 亚洲五月天丁香| 岛国在线观看网站| 久久中文字幕人妻熟女| 国产1区2区3区精品| 一二三四社区在线视频社区8| 国产高清videossex| 啦啦啦免费观看视频1| 亚洲免费av在线视频| 国产三级中文精品| e午夜精品久久久久久久| 一级a爱片免费观看的视频| 特大巨黑吊av在线直播| 黄色女人牲交| 最新美女视频免费是黄的| 亚洲熟妇中文字幕五十中出| 国产激情偷乱视频一区二区| 亚洲av电影在线进入| 99久久99久久久精品蜜桃| 丝袜美腿诱惑在线| 欧美性长视频在线观看| 丝袜美腿诱惑在线| 国产亚洲精品av在线| 国内毛片毛片毛片毛片毛片| 欧美乱色亚洲激情| 大型av网站在线播放| 777久久人妻少妇嫩草av网站| 亚洲aⅴ乱码一区二区在线播放 | 免费高清视频大片| 成人欧美大片| 舔av片在线| 欧美高清成人免费视频www| 亚洲av成人av| 一区二区三区国产精品乱码| 91成年电影在线观看| 欧美 亚洲 国产 日韩一| 日日干狠狠操夜夜爽| 老司机深夜福利视频在线观看| 我要搜黄色片| 日韩国内少妇激情av| 亚洲一码二码三码区别大吗| 免费在线观看视频国产中文字幕亚洲| 性色av乱码一区二区三区2| 亚洲国产精品合色在线| 免费人成视频x8x8入口观看| 国产乱人伦免费视频| 国产精品久久视频播放| 国产亚洲精品久久久久久毛片| 丝袜人妻中文字幕| 999久久久精品免费观看国产| 操出白浆在线播放| 国产亚洲精品第一综合不卡| 蜜桃久久精品国产亚洲av| 美女黄网站色视频| 黄色毛片三级朝国网站| 中文字幕熟女人妻在线| 国产午夜精品论理片| 极品教师在线免费播放| 欧美 亚洲 国产 日韩一| 亚洲精品中文字幕一二三四区| 老司机深夜福利视频在线观看| 搡老妇女老女人老熟妇| 国产免费男女视频| 岛国在线免费视频观看| 在线十欧美十亚洲十日本专区| 亚洲七黄色美女视频| www日本在线高清视频| 搡老妇女老女人老熟妇| 日本熟妇午夜| 日日干狠狠操夜夜爽| 成人国产综合亚洲| 国产精品 国内视频| 久热爱精品视频在线9| 久久精品成人免费网站| 全区人妻精品视频| 青草久久国产| 麻豆久久精品国产亚洲av| 欧美人与性动交α欧美精品济南到| 久久婷婷成人综合色麻豆| 18禁黄网站禁片午夜丰满| 神马国产精品三级电影在线观看 | 欧美三级亚洲精品| 午夜精品久久久久久毛片777| 老熟妇仑乱视频hdxx| 精品一区二区三区四区五区乱码| 2021天堂中文幕一二区在线观| 久久草成人影院| 午夜a级毛片| 中文在线观看免费www的网站 | 亚洲国产精品久久男人天堂| 午夜免费观看网址| 亚洲av成人精品一区久久| 麻豆国产97在线/欧美 |