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

    固原市隆德縣降水特性分析

    2023-02-01 14:33:10許昊麗李春光李超超田軍倉
    寧夏工程技術(shù) 2023年4期
    關(guān)鍵詞:隆德縣降水強(qiáng)度檢驗(yàn)法

    許昊麗,李春光,李超超,田軍倉,王 培

    (1.寧夏大學(xué) 土木與水利工程學(xué)院,寧夏 銀川 750021;2.北方民族大學(xué) 數(shù)學(xué)與信息科學(xué)學(xué)院,寧夏 銀川 750021;3.寧夏回族自治區(qū)水文水資源監(jiān)測預(yù)警中心,寧夏 銀川 750021;4.首創(chuàng)環(huán)保集團(tuán),北京 100000)

    我國幅員遼闊,不同地帶的氣候與降水變化呈現(xiàn)差異性。大氣中溫室氣體含量的升高導(dǎo)致全球氣候變暖,我國又屬于對氣候變化較為敏感的區(qū)域,氣候以及水文要素的變化不僅是導(dǎo)致旱澇等自然災(zāi)害頻繁發(fā)生的關(guān)鍵因素,也是影響農(nóng)業(yè)、生態(tài)系統(tǒng)平衡的重要因素[1]。T.R.Karl 等[2]、田立鑫等[3]指出降水演變分析不僅是水資源潛力分析的基礎(chǔ),也是探討氣候變化的主要手段之一。H.Tabari 等[4]通過M-K檢驗(yàn)法、Sen's 坡度法和線性回歸法對伊朗41 個雨量站的年降水進(jìn)行了時空分布分析。趙恩來等[5]采用線性回歸法、距平分析法、M-K 檢驗(yàn)法、小波分析法等方法,研究了駐馬店市各時間段降水的趨勢性、豐枯性及周期性。盧雅婷等[6]通過M-K 檢驗(yàn)法、線性傾向分析法和Morlet 小波分析法對金沙江流域60 a的降水進(jìn)行了特征分析,從而為該地區(qū)水資源的合理利用及生態(tài)環(huán)境保護(hù)提供了科學(xué)參考。王翠等[7]根據(jù)葉爾羌流域近50 a 的水文氣象資料,分析了極端干旱流域內(nèi)降水的時空分布特征。H.L.Yang等[8]采用創(chuàng)新趨勢分析法對寧夏58 a 來不同強(qiáng)度的降水序列進(jìn)行了時空變化分析。目前,針對降水特征的研究主要集中在大尺度范圍內(nèi),而針對西北地區(qū)小尺度范圍進(jìn)行的降水分析才能夠更好地指導(dǎo)當(dāng)?shù)剞r(nóng)業(yè)發(fā)展及水資源規(guī)劃。

    隆德縣位于寧夏南部山區(qū),人均水資源占有量為400 m3,屬于嚴(yán)重缺水地區(qū),且縣內(nèi)水資源分布與人口、土地資源分布不相適應(yīng)[9]。楊宗孝[10]對隆德縣水資源利用及供需情況進(jìn)行了分析,指出隆德縣水資源時空分布不均,其中農(nóng)業(yè)用水最多,其次為工業(yè)用水,并且未來10 a 內(nèi)該地區(qū)的水資源供需矛盾會進(jìn)一步加劇?;诖?,本文采用滑動平均法、M-K檢驗(yàn)法、Sen's 坡度法、R/S 分析法、小波分析法及多元線性回歸法等方法,從趨勢分析、周期分析、貢獻(xiàn)度分析3 個方面入手,分析了隆德縣近59 a 的降水變化特性和影響時段降水的主要因素,以期能為當(dāng)?shù)氐霓r(nóng)業(yè)生產(chǎn)、水資源開發(fā)及可持續(xù)利用提供參考。

    1 研究區(qū)概況

    寧夏固原市隆德縣總面積為999.45 km2,年平均氣溫為5.1 ℃,多年平均水面蒸發(fā)量為900 mm(圖1)。隆德縣的水資源主要來源于大氣降水,而降水量年內(nèi)分配不均,主要集中在7—9 月,6—9 月的降水量為年內(nèi)最大降水量,約占年降水量的70%;4—6 月為作物生長需水關(guān)鍵期,但這期間的降水量僅占年降水量的25%左右,因此,春旱發(fā)生較為普遍。此外,區(qū)域內(nèi)降水量呈現(xiàn)由東向西遞減的趨勢。降水的時空分布不均性導(dǎo)致隆德縣出現(xiàn)嚴(yán)重的季節(jié)性干旱[8-11]。

    圖1 研究區(qū)概況圖

    2 數(shù)據(jù)來源與研究方法

    2.1 數(shù)據(jù)來源

    文中采用的數(shù)據(jù)來源于1961—2019 年隆德縣逐日降水監(jiān)測資料,該資料由寧夏水文水資源勘測局提供。本文按氣象部門的標(biāo)準(zhǔn)劃分四季,并根據(jù)《寧夏回族自治區(qū)抗旱防汛條例》[12],將研究區(qū)內(nèi)的6—9 月定為汛期階段。

    2.2 研究方法

    2.2.1 趨勢分析方法 歷史趨勢分析主要采用氣候傾向率、M-K 檢驗(yàn)法和Sen's 坡度法進(jìn)行計(jì)算,能夠綜合確定不同時段降水量的變化趨勢。氣候傾向率以10 a 平均變率為依據(jù),通過求解時間序列的線性回歸方程來推斷序列的變化趨勢。周期變化和隨機(jī)波動會對水文時間序列造成影響,對其進(jìn)行線性回歸時,整體變化趨勢易受單個因子的影響,因此,本文輔以M-K 檢驗(yàn)法和Sen's 坡度法,綜合確定了各時段降水量的變化趨勢。其中M-K 檢驗(yàn)法是一種非參數(shù)統(tǒng)計(jì)檢驗(yàn)法,由于其不受異常值和序列分布特征的影響,在水文氣象要素分析中被世界氣象組織推薦并廣泛使用。而Sen's 坡度法采用非參數(shù)估計(jì)量的線性擬合值,能夠有效評估序列趨勢及變化幅度,并且異常值和序列長度對其分析結(jié)果的影響較小,通常可將該方法與M-K 檢驗(yàn)法相結(jié)合,綜合確定序列趨勢。趨勢分析的計(jì)算公式為[13]

    式中,Xj,Xk分別代表時間序列為j,k時的值。將計(jì)算得到的統(tǒng)計(jì)值βi按照從小到大的順序進(jìn)行排列,計(jì)算Sen's 坡度值。若統(tǒng)計(jì)值中值βmed>0,表示序列具有上升趨勢,反之為下降趨勢。

    Hurst 指數(shù)法通過分析序列未來變化趨勢與歷史趨勢的相關(guān)關(guān)系來判定序列的未來變化態(tài)勢,其基本原理見文獻(xiàn)[14]。按照Hurst 指數(shù)值的范圍,可將相關(guān)性分為若干等級(表1)。

    表1 Hurst指數(shù)值分布表

    2.2.2 周期分析方法 Morlet 連續(xù)的復(fù)小波變換分析能夠解析出隨機(jī)分布時間序列的多種周期變化,同時可以定性評估出未來趨勢。小波函數(shù)的變化公式為[15]

    式中:Wf(a,b)為變化函數(shù);a為尺度因子;b為平移因子。

    小波方差的公式為

    此外,通過繪制小波系數(shù)實(shí)部等值線圖可以確定序列的變化周期;通過小波方差圖可以確定序列的主、次周期。

    2.2.3 貢獻(xiàn)度分析方法 降水頻率與降水強(qiáng)度的變化會引起時段內(nèi)降水總量的變化。E.Lu 等[16]采用多元線性回歸法分析了兩種變化對降水總量變化的貢獻(xiàn)程度,計(jì)算公式為

    式中:P為f天內(nèi)的總降水量,mm;i為濕日內(nèi)平均降水強(qiáng)度,mm/d;f為時段內(nèi)日降水量>1 mm 的天數(shù),d;a,b,c為多元線性回歸方程(式(6))的回歸系數(shù),由最小二乘法估計(jì)得出,mm/d,d,mm;Sf,Si分別為f和i對降水總量的貢獻(xiàn)度;αf為f的變化尺度,d;αi為i的變化尺度,mm/d。

    3 結(jié)果分析與討論

    3.1 年降水量趨勢變化分析

    本文選用5 a 滑動平均法和M-K 檢驗(yàn)法(顯著性水平為0.05)對隆德縣1961—2019 年的年降水序列進(jìn)行了分析。隆德縣降水年際變化過程見圖2。由圖2 可知:隆德縣年降水的年際變化波動較大;多年平均降水量為519.67 mm;1961 年出現(xiàn)年最大降水量(825.93 mm),最小降水量(320.90 mm)出現(xiàn)在1982 年,極值比為2.57。由線性擬合結(jié)果可知,隆德縣年內(nèi)降水量總體呈現(xiàn)下降趨勢,下降速率為0.89 mm/a。由5 a 滑動平均曲線可知:年降水量波動幅度較大,1961—1971 年及2003—2007 年降水量呈下降趨勢;1971—1977 年及2006 年之后降水量呈上升趨勢;1977—2003 年降水量呈波動變化,但變化幅度較小,表明降水量較為平穩(wěn)。經(jīng)計(jì)算,M-K 檢驗(yàn)法中的統(tǒng)計(jì)量Z為-0.93。由于|Z|<1.96,表明隆德縣年降水量存在不顯著的減少趨勢。同時βmed=-0.7,表明隆德縣年降水量以7 mm/(10 a)的速率減小,但減小趨勢并不明顯。圖3 為隆德縣年降水量距平的年際變化圖。由圖3 可知,年降水量距平值正負(fù)交替變化,但總體基本持平。

    圖2 隆德縣降水年際變化圖

    圖3 隆德縣年降水量距平的年際變化圖

    本文采用Morlet 小波分析法對降水序列進(jìn)行了分析,得到隆德縣年降水量小波系數(shù)實(shí)部等值線圖(圖4a)。由圖4a 可以得出降水量的時間尺度及其位相分布,其中26~30 a,20~24 a 的時間尺度遍布整個研究時段,其中心尺度分別為28 a 和22 a,在該時間尺度下,降水量豐枯變化規(guī)律顯著;此外,10~14 a 的時間尺度較為明顯,但相較于前2 個時間尺度,其分布具有局限性,主要在1961—1979 年、2002—2019 年有所表現(xiàn);10 a 以下的時間尺度過小,正負(fù)交替頻繁,表明該地區(qū)在小時間尺度下年降水量波動頻繁。

    圖4 隆德縣年降水量小波分析圖

    隆德縣年降水序列小波方差圖如圖4b 所示,其反映了降水時間序列在不同時間尺度下的波動能量變化,圖中極大值點(diǎn)表示降水演化過程中潛在的主周期。由圖4b 可知:該地區(qū)降水量存在多時間尺度的周期變化;有3 個明顯的極大值點(diǎn),分別為28,22,12 a;6 a 和4 a 也對應(yīng)2 個極值點(diǎn),但并不顯著。根據(jù)水文時間序列的互相關(guān)關(guān)系理論,本文研究對象的最大潛在周期為28 a,而時間序列長度為59 a,因此,由小波方差圖確定的周期是準(zhǔn)確可靠的。對于隆德縣1961—2019 年的降水序列,28 a 為該區(qū)域降水變化的第1 主周期,22 a 為第2 主周期,2 個周期的波動決定了降水量變化的主要特性。

    本文分別給出了28,22,12 a 時間尺度下隆德縣年降水小波系數(shù)實(shí)部變化過程圖(圖5),系數(shù)為正表明為豐水期,反之為枯水期。由圖5 可知,隆德縣實(shí)際降水過程是不同時間周期降水過程的疊加。在28 a 時間尺度上,1961—2019 年隆德縣的降水經(jīng)歷了“豐—枯—豐—枯—豐—枯—豐”共7 次的豐枯交替。由此可以預(yù)測,在未來的一段時間內(nèi),該區(qū)域處于偏豐時期。在22 a 的時間尺度上,豐枯交替有所差異,但差異并不顯著,共經(jīng)歷了9 次的豐枯交替。綜上所述,在3 個時間尺度下,本文均可預(yù)見未來一段時間內(nèi)隆德縣處于降水偏豐時期。同時,本文結(jié)合R/S 法對該區(qū)域年降水量未來變化趨勢進(jìn)行了分析,其Hurst 指數(shù)為0.62(大于0.5),表明隆德縣未來降水量的變化對歷史趨勢具有較弱的持續(xù)性,即未來一段時間內(nèi)隆德縣降水量雖呈下降趨勢,但仍處于偏豐期,需注意極端強(qiáng)降水事件的發(fā)生。

    圖5 隆德縣年降水小波系數(shù)實(shí)部變化過程圖

    3.2 季節(jié)降水量趨勢變化分析

    隆德縣屬于典型的中溫帶季風(fēng)區(qū)半濕潤向半干旱過渡性氣候,具有春低溫少雨、夏短暫多雹、秋陰澇霜旱、冬嚴(yán)寒綿長的季節(jié)特征,且各個季節(jié)的降水量差異性較大。圖6 為隆德縣四季降水量箱線圖,其中3—5 月為春季,6—8 月為夏季,9—11 月為秋季,12 月—次年2 月為冬季。1961—2019 年隆德縣四季降水量統(tǒng)計(jì)值見表2。由圖6 和表2 可知,(1)隆德縣年內(nèi)降水量主要集中在夏季,其多年平均值為292.99 mm,約占年平均降水量的56.4 %;通過上、下4 分位之間的間距以及極值差可知,夏季降水量的波動范圍最大,但極值差與多年平均值的比值最?。?.10),且夏季降水量的均值略高于中位數(shù),呈現(xiàn)右偏態(tài)分布。(2)隆德縣冬季多年平均值最小,僅為12.18 mm,占年平均降水量的2.3%;與其他3季相比,冬季降水量的波動范圍最小,集中于0 值附近,極值差與多年平均值的比值較大(1.88)。(3)隆德縣春、秋兩季多年平均值分別為86.92 mm(占年平均降水量的16.7%)和127.58 mm(占年平均降水量的24.6%),其極值差與多年平均值之比分別為1.60 和2.07(最大);春、秋兩季降水量波動范圍較為適中,由于春季降水量均值略小于中位數(shù),可知其呈左偏態(tài)分布,而秋季降水量均值高于中位數(shù),可知其呈右偏態(tài)分布,且右偏態(tài)分布程度較大。(4)由Z值可知,春、冬兩季的降水量具有增加趨勢,且冬季降水量增加趨勢較為明顯,而夏、秋兩季的降水量減少趨勢不明顯。(5)結(jié)合Hurst 指數(shù)可知,未來春季降水更為隨機(jī),夏季降水則對變化趨勢具有較弱的反持續(xù)性,而秋、冬兩季的降水量變化趨勢具有較強(qiáng)的持續(xù)性。

    表2 隆德縣四季降水量統(tǒng)計(jì)表

    圖6 隆德縣四季降水量箱線圖

    3.3 降水趨勢決定因素分析

    一定時段內(nèi)降水序列變化的差異性主要受降水頻率和降水強(qiáng)度的影響。本文依據(jù)式(5),分析了隆德縣降水頻率和降水強(qiáng)度在各時段降水變化趨勢中的貢獻(xiàn)度大小,得到各時段內(nèi)Sf和Si的分布散點(diǎn)圖(圖7)。由圖7 可知:僅有冬季降水貢獻(xiàn)度值在直線Sf=Si以下,其偏離程度較小,表明冬季降水量年際間的差異性受降水頻率的影響居多,但整體降水量的變化主要是受降水頻率和降水強(qiáng)度的雙重影響;其他3 季降水貢獻(xiàn)度值均在直線Sf=Si以上,其中春、秋兩季貢獻(xiàn)度值的偏離程度較小,表明這兩季降水量年際間的變化受降水頻率和降水強(qiáng)度的共同影響;而夏季降水量的貢獻(xiàn)度值偏離程度較大,表明其降水量年際間的變化主要受降水強(qiáng)度的影響。此外,夏季降水量占年降水量的比重較大,因此,隆德縣年降水總量的變化主要取決于降水強(qiáng)度的變化,降水強(qiáng)度的減少導(dǎo)致年降水總量逐漸減少。

    圖7 隆德縣各時段內(nèi)Sf和Si散點(diǎn)圖

    4 結(jié)論

    本文基于隆德縣1961—2019 年逐日降水?dāng)?shù)據(jù),運(yùn)用5 a 滑動平均法、M-K 檢驗(yàn)法、R/S 分析法等方法,通過趨勢分析、周期分析、貢獻(xiàn)度分析等過程,研究了隆德縣降水在年際和季節(jié)尺度上的降水特性及影響不同時段降水量變化的關(guān)鍵因素,得出以下結(jié)論。

    (1)隆德縣近59 年降水量以8.94 mm/(10 a)的速率減少,但減少趨勢并不明顯。各季節(jié)中,春、冬兩季呈現(xiàn)降水量增加趨勢,其中冬季增加趨勢較為明顯,而夏、秋季的降水量呈不明顯減少趨勢,并且除春、夏季之外,年降水以及季節(jié)降水在未來一段時間內(nèi)均延續(xù)歷史趨勢。年內(nèi)降水主要集中在夏季,約占年平均降水量的56.4 %;冬季占比最小,占年平均降水量的比重不足3%。降水幅度中,夏季變化最大,春、秋季次之,冬季最小。

    (2)隆德縣降水存在豐枯交替的周期變化,年際間變化以28 a 為第1 主周期,并伴有20 a 和12 a的周期變化。根據(jù)第1 主周期的變化,隆德縣年降水量在2019—2024 年處于相對偏豐期,應(yīng)注意做好城市防洪和農(nóng)田排澇工作。此外,隆德縣降水量未來呈減小趨勢,需進(jìn)一步考慮雨水的合理利用,充分發(fā)揮隆德縣骨干水庫的調(diào)洪和滯蓄作用,滿足干旱時間段的蓄水要求并實(shí)現(xiàn)水資源可持續(xù)利用。

    (3)隆德縣降水量年際間的差異性主要受降水強(qiáng)度的影響,降水強(qiáng)度的減小是引起年降水量減少的主要原因;而春、秋和冬季降水量變化受降水頻率和降水強(qiáng)度的共同影響,夏季降水量的變化主要受降水強(qiáng)度的影響,同時夏季的極端降水可能出現(xiàn)減小趨勢。

    猜你喜歡
    隆德縣降水強(qiáng)度檢驗(yàn)法
    1961—2021年濉溪縣極端降水變化分析
    記一次寫生
    速讀·上旬(2021年7期)2021-08-05 02:33:55
    不同降水強(qiáng)度下風(fēng)廓線雷達(dá)譜矩特征與測風(fēng)準(zhǔn)確性分析
    2007—2017年大連市降水及降水日數(shù)變化特征分析
    寧夏隆德縣農(nóng)村飲水安全工程建設(shè)經(jīng)驗(yàn)及做法
    隆德縣三通古碑考釋
    西夏研究(2017年4期)2017-08-17 08:53:10
    PCR 檢驗(yàn)法和細(xì)菌培養(yǎng)法用于陰道細(xì)菌檢驗(yàn)的效果
    石羊河流域人工增雨雪的效果評估分析
    城市地理(2015年7期)2015-07-13 04:39:36
    隆德縣農(nóng)村電影放映創(chuàng)新發(fā)展之路
    關(guān)于協(xié)方差的U統(tǒng)計(jì)量檢驗(yàn)法
    99热网站在线观看| 欧美成人精品欧美一级黄| 在线观看一区二区三区激情| 亚洲人成77777在线视频| 国产精品一国产av| 在线观看国产h片| 亚洲综合精品二区| 91成人精品电影| 日韩熟女老妇一区二区性免费视频| 日韩一本色道免费dvd| 中文字幕色久视频| 五月天丁香电影| 国产亚洲精品第一综合不卡| 在线看a的网站| 七月丁香在线播放| 男男h啪啪无遮挡| 晚上一个人看的免费电影| 欧美黄色片欧美黄色片| 日韩大码丰满熟妇| 午夜免费鲁丝| 久久人人爽人人片av| 国产日韩欧美在线精品| 亚洲第一区二区三区不卡| 极品人妻少妇av视频| 亚洲精品av麻豆狂野| 999精品在线视频| 欧美激情 高清一区二区三区| 母亲3免费完整高清在线观看| 久久久久国产一级毛片高清牌| 丰满少妇做爰视频| 狠狠精品人妻久久久久久综合| 高清视频免费观看一区二区| 欧美日韩福利视频一区二区| 99久久99久久久精品蜜桃| 国产欧美日韩一区二区三区在线| 亚洲av日韩精品久久久久久密 | 亚洲人成77777在线视频| 欧美日韩一级在线毛片| 好男人视频免费观看在线| 女的被弄到高潮叫床怎么办| 黄色视频不卡| 菩萨蛮人人尽说江南好唐韦庄| 在线观看免费高清a一片| 18禁动态无遮挡网站| 欧美人与善性xxx| 免费人妻精品一区二区三区视频| 久久综合国产亚洲精品| 中文字幕制服av| 国产一区二区三区综合在线观看| 日韩一卡2卡3卡4卡2021年| 卡戴珊不雅视频在线播放| 老司机靠b影院| 亚洲在久久综合| 国产精品偷伦视频观看了| 青春草亚洲视频在线观看| 下体分泌物呈黄色| 欧美老熟妇乱子伦牲交| 亚洲国产精品一区三区| 久久 成人 亚洲| 无限看片的www在线观看| 婷婷色av中文字幕| videosex国产| 自拍欧美九色日韩亚洲蝌蚪91| 2021少妇久久久久久久久久久| 亚洲国产毛片av蜜桃av| 你懂的网址亚洲精品在线观看| 国产97色在线日韩免费| 国产黄色视频一区二区在线观看| 尾随美女入室| 欧美人与性动交α欧美精品济南到| 色视频在线一区二区三区| 亚洲三区欧美一区| 大香蕉久久网| 亚洲伊人久久精品综合| 久久青草综合色| 纯流量卡能插随身wifi吗| www.av在线官网国产| 99香蕉大伊视频| 卡戴珊不雅视频在线播放| 老熟女久久久| 男女下面插进去视频免费观看| 日韩av免费高清视频| 亚洲国产精品成人久久小说| 深夜精品福利| 久久这里只有精品19| 99九九在线精品视频| 91成人精品电影| 亚洲欧美一区二区三区黑人| 国产亚洲最大av| 日韩一区二区三区影片| 精品酒店卫生间| 亚洲国产av新网站| 中文乱码字字幕精品一区二区三区| 美女午夜性视频免费| 国产成人精品久久二区二区91 | 人人妻,人人澡人人爽秒播 | 欧美日韩成人在线一区二区| 最近中文字幕高清免费大全6| 飞空精品影院首页| 香蕉丝袜av| xxx大片免费视频| 97人妻天天添夜夜摸| 国产成人午夜福利电影在线观看| 在线观看免费午夜福利视频| 亚洲四区av| 日韩精品免费视频一区二区三区| 制服诱惑二区| 日韩免费高清中文字幕av| 日本欧美视频一区| 国产黄频视频在线观看| 午夜日韩欧美国产| 国产精品香港三级国产av潘金莲 | 国产亚洲最大av| 啦啦啦啦在线视频资源| 啦啦啦 在线观看视频| 天天操日日干夜夜撸| 免费黄频网站在线观看国产| 如何舔出高潮| 精品视频人人做人人爽| 欧美激情极品国产一区二区三区| av网站免费在线观看视频| 99热国产这里只有精品6| 777米奇影视久久| 国产成人a∨麻豆精品| 97精品久久久久久久久久精品| 97在线人人人人妻| 国产毛片在线视频| 成人毛片60女人毛片免费| 久久精品久久久久久噜噜老黄| 国产人伦9x9x在线观看| 不卡av一区二区三区| 黑人巨大精品欧美一区二区蜜桃| 欧美人与善性xxx| 最新的欧美精品一区二区| 两个人看的免费小视频| 国产野战对白在线观看| 欧美亚洲 丝袜 人妻 在线| 久久99热这里只频精品6学生| www.自偷自拍.com| 日韩一本色道免费dvd| 综合色丁香网| 国产亚洲av片在线观看秒播厂| 男女下面插进去视频免费观看| 国产高清国产精品国产三级| 久久久久久人人人人人| 欧美在线一区亚洲| av卡一久久| 国产成人系列免费观看| 亚洲精品日韩在线中文字幕| 99久久99久久久精品蜜桃| 精品视频人人做人人爽| bbb黄色大片| 欧美 亚洲 国产 日韩一| h视频一区二区三区| 国产一区二区三区av在线| 精品卡一卡二卡四卡免费| 一级爰片在线观看| 亚洲av中文av极速乱| 国精品久久久久久国模美| 日本91视频免费播放| 老鸭窝网址在线观看| 日本爱情动作片www.在线观看| 亚洲精品一区蜜桃| 亚洲精品,欧美精品| 热re99久久国产66热| 欧美av亚洲av综合av国产av | 中文乱码字字幕精品一区二区三区| 欧美老熟妇乱子伦牲交| 中文字幕色久视频| 国产精品久久久久久精品古装| 免费高清在线观看日韩| 欧美另类一区| 新久久久久国产一级毛片| 天天躁夜夜躁狠狠久久av| 日韩成人av中文字幕在线观看| 极品少妇高潮喷水抽搐| 超色免费av| 大片免费播放器 马上看| 精品一区二区三卡| 99热国产这里只有精品6| 国产精品熟女久久久久浪| 欧美黄色片欧美黄色片| 一级片'在线观看视频| 热re99久久精品国产66热6| 免费高清在线观看日韩| 精品一区二区三区av网在线观看 | 午夜av观看不卡| 久久精品国产综合久久久| 天天影视国产精品| 欧美人与性动交α欧美精品济南到| 国语对白做爰xxxⅹ性视频网站| 成人午夜精彩视频在线观看| 99久久精品国产亚洲精品| 国产精品秋霞免费鲁丝片| 男人爽女人下面视频在线观看| 欧美黑人欧美精品刺激| 在线观看三级黄色| 老汉色∧v一级毛片| 亚洲欧美精品自产自拍| 毛片一级片免费看久久久久| 国产成人一区二区在线| 亚洲欧洲日产国产| 新久久久久国产一级毛片| 丝袜人妻中文字幕| 我要看黄色一级片免费的| 亚洲国产精品999| 亚洲av男天堂| 黄网站色视频无遮挡免费观看| 免费高清在线观看日韩| 亚洲av电影在线观看一区二区三区| www.自偷自拍.com| 欧美精品一区二区免费开放| 2021少妇久久久久久久久久久| 国产精品二区激情视频| 丰满乱子伦码专区| 午夜影院在线不卡| 日本午夜av视频| 美女国产高潮福利片在线看| 亚洲成人av在线免费| 国产片内射在线| 国产在线一区二区三区精| 亚洲一区二区三区欧美精品| 好男人视频免费观看在线| 18禁国产床啪视频网站| 午夜久久久在线观看| 日韩制服丝袜自拍偷拍| 午夜av观看不卡| www日本在线高清视频| 久久精品国产亚洲av涩爱| 大码成人一级视频| 日日啪夜夜爽| 国产欧美日韩综合在线一区二区| 久久久精品区二区三区| 日韩人妻精品一区2区三区| 91精品伊人久久大香线蕉| 赤兔流量卡办理| 日日撸夜夜添| 久久久久久人人人人人| 汤姆久久久久久久影院中文字幕| 不卡av一区二区三区| 岛国毛片在线播放| 丰满迷人的少妇在线观看| 波多野结衣av一区二区av| 国产 一区精品| 下体分泌物呈黄色| 国产av码专区亚洲av| 午夜免费鲁丝| 欧美老熟妇乱子伦牲交| 又粗又硬又长又爽又黄的视频| 三上悠亚av全集在线观看| 精品福利永久在线观看| 国产成人免费无遮挡视频| 人妻 亚洲 视频| 日韩制服骚丝袜av| 中文字幕制服av| 视频在线观看一区二区三区| 下体分泌物呈黄色| 精品免费久久久久久久清纯 | 婷婷色麻豆天堂久久| 中文字幕最新亚洲高清| 蜜桃国产av成人99| 国产成人精品无人区| 精品一区二区三区av网在线观看 | 色视频在线一区二区三区| 青春草国产在线视频| 免费在线观看完整版高清| 91aial.com中文字幕在线观看| 黄色毛片三级朝国网站| 亚洲伊人色综图| 成年美女黄网站色视频大全免费| 一级片'在线观看视频| 99热全是精品| 精品人妻熟女毛片av久久网站| 少妇人妻 视频| 精品久久久精品久久久| 精品国产一区二区三区四区第35| 久久久久久久久久久久大奶| 丰满迷人的少妇在线观看| 男人舔女人的私密视频| 51午夜福利影视在线观看| 久久久精品94久久精品| 女人高潮潮喷娇喘18禁视频| 国产熟女午夜一区二区三区| 一二三四在线观看免费中文在| 国产有黄有色有爽视频| 大片免费播放器 马上看| 中文字幕人妻丝袜一区二区 | 19禁男女啪啪无遮挡网站| 一级毛片电影观看| 日韩欧美一区视频在线观看| 亚洲综合色网址| 国产深夜福利视频在线观看| bbb黄色大片| 超碰97精品在线观看| 久久免费观看电影| 亚洲国产精品成人久久小说| 大话2 男鬼变身卡| 2021少妇久久久久久久久久久| 亚洲成色77777| 黄色怎么调成土黄色| 考比视频在线观看| 免费av中文字幕在线| 免费黄频网站在线观看国产| 精品一区二区三区四区五区乱码 | 丝袜喷水一区| 女人精品久久久久毛片| 韩国精品一区二区三区| 午夜老司机福利片| 亚洲av成人不卡在线观看播放网 | 亚洲激情五月婷婷啪啪| 最近2019中文字幕mv第一页| 伊人亚洲综合成人网| av不卡在线播放| 99精国产麻豆久久婷婷| 久久韩国三级中文字幕| 女性生殖器流出的白浆| 毛片一级片免费看久久久久| 午夜福利免费观看在线| 嫩草影院入口| 免费高清在线观看日韩| 欧美日韩国产mv在线观看视频| 精品酒店卫生间| 日韩成人av中文字幕在线观看| 国产一卡二卡三卡精品 | 无遮挡黄片免费观看| 国产麻豆69| 18禁观看日本| 国产麻豆69| 午夜精品国产一区二区电影| 美女主播在线视频| 精品国产超薄肉色丝袜足j| 男女边摸边吃奶| 亚洲精品日本国产第一区| 天天添夜夜摸| 精品第一国产精品| 天堂俺去俺来也www色官网| netflix在线观看网站| 国产又爽黄色视频| 亚洲精品久久久久久婷婷小说| 91精品伊人久久大香线蕉| 久久人人97超碰香蕉20202| 日韩人妻精品一区2区三区| 女人爽到高潮嗷嗷叫在线视频| 欧美av亚洲av综合av国产av | 人人妻人人添人人爽欧美一区卜| 午夜福利视频在线观看免费| 久久女婷五月综合色啪小说| 高清av免费在线| 色精品久久人妻99蜜桃| 免费av中文字幕在线| 国产亚洲av高清不卡| 成人免费观看视频高清| 成人国语在线视频| 欧美成人午夜精品| 色网站视频免费| 制服人妻中文乱码| 亚洲成人免费av在线播放| 18在线观看网站| 亚洲综合精品二区| 99re6热这里在线精品视频| 国产男人的电影天堂91| 18禁动态无遮挡网站| 老鸭窝网址在线观看| 在线观看免费高清a一片| 久久久精品免费免费高清| 日本wwww免费看| 日韩一区二区三区影片| 免费女性裸体啪啪无遮挡网站| 日本爱情动作片www.在线观看| 男女无遮挡免费网站观看| 日韩不卡一区二区三区视频在线| 精品少妇一区二区三区视频日本电影 | 国产免费现黄频在线看| 9191精品国产免费久久| 亚洲成色77777| 大片免费播放器 马上看| 久久女婷五月综合色啪小说| 日韩免费高清中文字幕av| 欧美日韩亚洲综合一区二区三区_| 精品国产国语对白av| 午夜福利在线免费观看网站| 又黄又粗又硬又大视频| 欧美在线一区亚洲| 久久人人爽人人片av| 十八禁人妻一区二区| 亚洲欧美激情在线| 欧美日韩视频精品一区| 狠狠精品人妻久久久久久综合| 午夜av观看不卡| 水蜜桃什么品种好| 一级爰片在线观看| 免费女性裸体啪啪无遮挡网站| xxx大片免费视频| 九九爱精品视频在线观看| 一区在线观看完整版| 热re99久久精品国产66热6| 又大又爽又粗| 热99国产精品久久久久久7| 一级毛片黄色毛片免费观看视频| 一区二区三区四区激情视频| 日韩视频在线欧美| 日韩一卡2卡3卡4卡2021年| 啦啦啦啦在线视频资源| 国产成人精品福利久久| 久久久精品94久久精品| 99久久人妻综合| 69精品国产乱码久久久| 亚洲欧美一区二区三区国产| 在线天堂最新版资源| 亚洲国产欧美网| 岛国毛片在线播放| 久久青草综合色| 久久毛片免费看一区二区三区| 99久久精品国产亚洲精品| 中文字幕制服av| 男人操女人黄网站| 一级a爱视频在线免费观看| 大话2 男鬼变身卡| 午夜91福利影院| 免费在线观看完整版高清| √禁漫天堂资源中文www| 精品国产一区二区三区四区第35| 亚洲免费av在线视频| av在线app专区| 国产成人免费观看mmmm| 欧美亚洲日本最大视频资源| 久久天躁狠狠躁夜夜2o2o | 久久午夜综合久久蜜桃| av在线app专区| 欧美人与性动交α欧美精品济南到| 少妇人妻精品综合一区二区| 狠狠婷婷综合久久久久久88av| 久久精品久久精品一区二区三区| 国产午夜精品一二区理论片| 精品久久蜜臀av无| 亚洲美女搞黄在线观看| 免费人妻精品一区二区三区视频| 曰老女人黄片| 精品国产露脸久久av麻豆| 国产99久久九九免费精品| 免费看不卡的av| 精品午夜福利在线看| 国产乱人偷精品视频| 成人黄色视频免费在线看| 97精品久久久久久久久久精品| 欧美xxⅹ黑人| 精品国产国语对白av| 亚洲,欧美,日韩| www.精华液| 国产日韩欧美在线精品| 少妇的丰满在线观看| 久久精品久久精品一区二区三区| 男女床上黄色一级片免费看| 两个人看的免费小视频| 久久久久精品人妻al黑| 韩国精品一区二区三区| 免费高清在线观看视频在线观看| 最近的中文字幕免费完整| 一级毛片黄色毛片免费观看视频| 午夜91福利影院| 香蕉丝袜av| 我要看黄色一级片免费的| 久久人妻熟女aⅴ| 国产精品嫩草影院av在线观看| 亚洲av日韩在线播放| 又黄又粗又硬又大视频| 免费女性裸体啪啪无遮挡网站| 91精品伊人久久大香线蕉| 久久久国产精品麻豆| 精品少妇久久久久久888优播| 人人妻人人爽人人添夜夜欢视频| 国产一区二区 视频在线| 欧美黄色片欧美黄色片| 国产成人av激情在线播放| 亚洲久久久国产精品| 成人18禁高潮啪啪吃奶动态图| 亚洲国产毛片av蜜桃av| 各种免费的搞黄视频| 丝袜喷水一区| 亚洲伊人久久精品综合| 国产野战对白在线观看| 中文字幕色久视频| 天天躁狠狠躁夜夜躁狠狠躁| 成人国语在线视频| 国产精品久久久av美女十八| 中国三级夫妇交换| 极品少妇高潮喷水抽搐| 欧美日韩福利视频一区二区| 免费看av在线观看网站| 国产精品偷伦视频观看了| 亚洲熟女毛片儿| 国产精品国产三级专区第一集| 亚洲欧美精品综合一区二区三区| 在线看a的网站| 免费在线观看视频国产中文字幕亚洲 | 国产精品国产三级专区第一集| 韩国av在线不卡| 日韩大码丰满熟妇| 在线免费观看不下载黄p国产| 国产精品 国内视频| 黄色视频在线播放观看不卡| 交换朋友夫妻互换小说| www日本在线高清视频| 亚洲国产精品一区二区三区在线| 99久国产av精品国产电影| 亚洲第一区二区三区不卡| 日韩免费高清中文字幕av| 久久天躁狠狠躁夜夜2o2o | 五月开心婷婷网| 色婷婷久久久亚洲欧美| 免费观看a级毛片全部| 在线观看一区二区三区激情| 综合色丁香网| 国产免费福利视频在线观看| 欧美成人午夜精品| 99热网站在线观看| 18禁国产床啪视频网站| 女人爽到高潮嗷嗷叫在线视频| 1024视频免费在线观看| 久久精品国产a三级三级三级| 日韩av免费高清视频| 欧美人与善性xxx| 夜夜骑夜夜射夜夜干| 久久99精品国语久久久| 久久久国产欧美日韩av| 91老司机精品| 你懂的网址亚洲精品在线观看| 成年人午夜在线观看视频| 日韩人妻精品一区2区三区| 国产精品99久久99久久久不卡 | 亚洲人成网站在线观看播放| 大片电影免费在线观看免费| 亚洲国产毛片av蜜桃av| 亚洲国产av影院在线观看| 嫩草影视91久久| 精品一区二区三区四区五区乱码 | 曰老女人黄片| 麻豆av在线久日| 亚洲欧美日韩另类电影网站| 只有这里有精品99| 蜜桃国产av成人99| 99香蕉大伊视频| 九九爱精品视频在线观看| 99热网站在线观看| 搡老岳熟女国产| 国产成人啪精品午夜网站| 两个人免费观看高清视频| 亚洲精品国产色婷婷电影| 午夜福利,免费看| av在线播放精品| 国产精品一区二区在线观看99| av国产久精品久网站免费入址| 国产av码专区亚洲av| 在线 av 中文字幕| 天天躁日日躁夜夜躁夜夜| 欧美亚洲 丝袜 人妻 在线| 丁香六月天网| 午夜免费鲁丝| 亚洲精品aⅴ在线观看| 18禁观看日本| 国产亚洲最大av| 99精国产麻豆久久婷婷| 亚洲精品在线美女| 日韩中文字幕视频在线看片| 精品人妻熟女毛片av久久网站| 久久久久网色| 午夜老司机福利片| 宅男免费午夜| 日本欧美国产在线视频| 99香蕉大伊视频| 日韩欧美一区视频在线观看| av国产久精品久网站免费入址| 一区二区三区精品91| 午夜影院在线不卡| 亚洲精品av麻豆狂野| 日韩av不卡免费在线播放| 女人高潮潮喷娇喘18禁视频| 欧美日韩视频高清一区二区三区二| 亚洲av日韩在线播放| 成年av动漫网址| 成人影院久久| av.在线天堂| 欧美日韩亚洲高清精品| 久久久精品94久久精品| 香蕉丝袜av| 晚上一个人看的免费电影| 国产精品欧美亚洲77777| 最近中文字幕2019免费版| 婷婷色av中文字幕| 欧美中文综合在线视频| 波野结衣二区三区在线| 啦啦啦啦在线视频资源| 亚洲精品国产一区二区精华液| 国产成人精品在线电影| 久久久久久免费高清国产稀缺| 高清av免费在线| 天天添夜夜摸| 美女视频免费永久观看网站| 亚洲国产中文字幕在线视频| 99国产综合亚洲精品| 老司机靠b影院| 久久久精品94久久精品| 中文字幕另类日韩欧美亚洲嫩草| 久久久久久久大尺度免费视频| 国产精品三级大全| 国产av码专区亚洲av| 女人被躁到高潮嗷嗷叫费观| 亚洲精品国产色婷婷电影| 黄频高清免费视频| 免费黄网站久久成人精品| 丰满迷人的少妇在线观看| 国产女主播在线喷水免费视频网站| 99久久人妻综合| 观看美女的网站| 伊人亚洲综合成人网|