0.3的指數(shù)作為預(yù)測因子組,采用逐步回歸統(tǒng)計(jì)法,建立荊州汛期降水預(yù)測模型和主汛期降水預(yù)測模型。結(jié)果表明,荊州汛期和主汛期降水預(yù)測模型的"/>
  • <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    荊州汛期降水趨勢模型研究

    2018-08-25 09:26:08鄧艷君趙卓勛張倫瑾
    湖北農(nóng)業(yè)科學(xué) 2018年10期
    關(guān)鍵詞:預(yù)測模型

    鄧艷君 趙卓勛 張倫瑾

    摘要:利用國家氣候中心每月下發(fā)的130項(xiàng)氣候系統(tǒng)監(jiān)測指數(shù)和荊州站1954-2016年的降水資料,逐一分析這些氣候指數(shù)與汛期和主汛期降水距平百分率的相關(guān)性,選取相關(guān)系數(shù)>0.3的指數(shù)作為預(yù)測因子組,采用逐步回歸統(tǒng)計(jì)法,建立荊州汛期降水預(yù)測模型和主汛期降水預(yù)測模型。結(jié)果表明,荊州汛期和主汛期降水預(yù)測模型的相關(guān)系數(shù)分別為0.874和0.914,均明顯大于單個(gè)因子的相關(guān)性。模型預(yù)測2016年汛期和主汛期的降水距平百分率,結(jié)果分別為偏多17.3%和偏多223.2%,與汛期降水距平百分率偏多6.2%和主汛期降水距平百分率偏多30.2%相比,汛期降水預(yù)測模型預(yù)測結(jié)果較好,同屬于略多的等級,主汛期降水預(yù)測模型預(yù)測結(jié)果雖能預(yù)測出偏多的趨勢,但數(shù)值明顯偏大,可能與該模型中預(yù)測因子上年11月印度副高強(qiáng)度指數(shù)有效數(shù)據(jù)較少,系數(shù)偏大有關(guān),可能需要更多數(shù)據(jù)來調(diào)整該項(xiàng)系數(shù),從而提高預(yù)測精度。

    關(guān)鍵詞:汛期降水;氣候監(jiān)測指數(shù);預(yù)測因子;預(yù)測模型

    中圖分類號:P468.0+24 文獻(xiàn)標(biāo)識碼:A 文章編號:0439-8114(2018)10-0054-06

    DOI:10.14088/j.cnki.issn0439-8114.2018.10.013

    Prediction Model of Precipitation Tendency during Flood Season in Jingzhou

    DENG Yan-jun,ZHAO Zhuo-xun,ZHANG Lun-jin

    (Jingzhou Meteorological Bureau of Hubei Province,Jingzhou 434020,Hubei,China)

    Abstract: The 130 climate system monitoring index issued monthly by the National Climate Center and the precipitation data from 1954 to 2016 recorded by Jingzhou climate station were analyzed to find the correlation coefficient between these climate index and the flood season precipitation anomaly percentage as well as the main flood season precipitation anomaly percentage. Using stepwise regression method, the index with the correlation coefficients larger than 0.3 were chosen as predictor group to establish the Jingzhou precipitation forecast models of flood season and main flood season. The statistical results showed that the correlation coefficients of precipitation forecast models of flood season and main flood season were 0.874 and 0.914 respectively, which were significantly higher than the coefficients of single factor. The flood season and main flood season precipitation anomaly percentages of Jingzhou in 2016 were forecasted to be more than the normal 17.3% and 223.2% with those models. Comparing to the actual value of more than 6.2% and 30.2%, the flood season model gave better forecasting result and the main flood season model gave same trend but too large value. The large forecasting error of main flood season model may be related to the shortage of effective data of last November India subtropical high pressure zone index and larger coefficient. More data was needed to tune the relevant coefficient to improve the forecasting accuracy.

    Key words: flood season precipitation; climatic monitoring index; predictive factors; prediction model

    汛期降水對國民經(jīng)濟(jì)特別是農(nóng)業(yè)生產(chǎn)影響極大,其中大范圍的旱澇氣候?yàn)?zāi)害的影響尤為嚴(yán)重[1]。汛期降水是各級政府極為關(guān)注的重要問題,也是氣象學(xué)界長期以來研究的熱點(diǎn)和難點(diǎn)。湖北省荊州市地處亞熱帶季風(fēng)區(qū),汛期降水受中高緯度西風(fēng)環(huán)流、西太平洋副熱帶高壓、臺風(fēng)、東亞季風(fēng)和印度西南季風(fēng)環(huán)流等天氣系統(tǒng)的影響,具有年際變化大、影響因素多等特點(diǎn),容易造成洪澇與嚴(yán)重干旱等災(zāi)害。

    學(xué)者預(yù)測汛期降水方法歸納起來主要有經(jīng)驗(yàn)方法、統(tǒng)計(jì)方法和動(dòng)力學(xué)方法。譚桂容等[2]對華北旱澇年和整個(gè)北半球上各層的環(huán)流異常變化情況進(jìn)行了研究;韓照宇等[3]利用滑動(dòng)相關(guān)普查法分析表明一些環(huán)流因子和夏半年降水有很好的相關(guān)關(guān)系;羅伯良等[4]研究了歐亞中高緯環(huán)流形勢和極渦形態(tài)以及加利福尼亞寒流和黑潮區(qū)的海溫,對夏季降水異常起著重要作用;丑紀(jì)范[5,6]提出基于歷史相似信息對動(dòng)力模式誤差進(jìn)行預(yù)報(bào)的思路,據(jù)此眾多學(xué)者[7-10]進(jìn)行了汛期降水預(yù)測研究;還有學(xué)者[11-13]利用相似年的模式誤差信息,實(shí)現(xiàn)對預(yù)報(bào)年氣候模式預(yù)測結(jié)果訂正的汛期降水動(dòng)力統(tǒng)計(jì)客觀定量化預(yù)測的方法,有效改進(jìn)了模式預(yù)測結(jié)果。本研究旨在建立一種經(jīng)驗(yàn)與統(tǒng)計(jì)相結(jié)合的預(yù)測模型,分析氣候系統(tǒng)監(jiān)測指數(shù)集與汛期降水距平百分率的相關(guān)關(guān)系,建立荊州汛期降水的統(tǒng)計(jì)預(yù)測模型。

    1 研究資料

    氣候系統(tǒng)指數(shù)集產(chǎn)品來源于國家氣候中心每月下發(fā)的特征指數(shù)資料,可在氣象業(yè)務(wù)內(nèi)網(wǎng)[14]和國家氣候中心監(jiān)測診斷主頁(CMDP)[15]自行下載使用。氣候系統(tǒng)監(jiān)測指數(shù)集共包含130項(xiàng)氣候系統(tǒng)指數(shù),其中重新整編計(jì)算的指數(shù)114項(xiàng)(大氣環(huán)流指數(shù)88項(xiàng),海溫指數(shù)26項(xiàng)),收集下載的指數(shù)16項(xiàng),時(shí)間分辨率為月平均,起止時(shí)間為1951年1月至2016年。汛期各時(shí)段降水量選用湖北省荊州區(qū)國家基本氣象站1954-2016年記錄的逐月降水量資料。

    2 方法

    2.1 數(shù)據(jù)整理

    通過Access、Excel等軟件將130項(xiàng)氣候系統(tǒng)監(jiān)測指數(shù)整理成按名稱、年、月導(dǎo)出,選用各月、季度平均及全年平均(或者總和)形式的資料,計(jì)算出荊州站汛期(5-9月)、主汛期(6-8月)降水距平百分率,然后將上年3月到本年2月逐月、季度平均及全年平均氣候系統(tǒng)監(jiān)測指數(shù)與該年汛期(5-9月)、主汛期(6-8月)降水距平百分率一起建立綜合數(shù)據(jù)庫。

    2.2 降水?dāng)?shù)據(jù)分析

    從1954-2016年荊州站汛期(5-9月)降水距平百分率的序列圖(圖1)中可以看出,汛期降水距平百分率呈振蕩式分布,汛期降水最多的年份為1954年,偏多96.59%,汛期降水最少的年份為1966年,偏少54.74%。

    從1954-2016年荊州站主汛期(6-8月)降水距平百分率的序列圖(圖2)中可以看出,主汛期降水距平百分率也呈振蕩式分布,且值域范圍更大,主汛期降水最多年份為1954年,偏多122.7%,主汛期降水最少的年份為1972年,偏少69.9%。

    根據(jù)湖北省月短期氣候預(yù)測質(zhì)量評定辦法,降水趨勢預(yù)測按照6級評分制評定,其預(yù)測用語為特少、偏少、略少、略多、偏多、特多。按該評定辦法統(tǒng)計(jì)分析1954-2016年荊州汛期和主汛期的降水情況,結(jié)果見表1。由表1可知,63年汛期降水中偏多的年份為31年,其中,特多的年份為5年,偏多的年份為11年,略多的年份為15年;偏少的年份為32年,其中,特少的年份為1年,偏少的年份為15年,略少的年份為16年,汛期降水偏多偏少的年份數(shù)基本持平。63年主汛期降水中偏多的年份為25年,其中,特多的年份為8年,偏多的年份為8年,略多的年份為10年;偏少的年份為38年,其中,特少的年份為5年,偏少的年份為20年,略少的年份為13年,主汛期降水偏少的年份數(shù)多于偏多的年份數(shù)。

    2.3 相關(guān)性分析

    運(yùn)用數(shù)理統(tǒng)計(jì)原理和方法,在SPSS統(tǒng)計(jì)軟件中,逐一分析各月、季度平均及全年平均氣候系統(tǒng)監(jiān)測指數(shù)與該年汛期(5-9月)、主汛期(6-8月)降水距平百分率的Pearson相關(guān)性,選出相關(guān)系數(shù)>0.3,雙側(cè)顯著性大于95%(P小于0.01代表在99%的置信度下顯著相關(guān),P小于0.05代表在95%的置信度下相關(guān))的預(yù)報(bào)因子。本著擇優(yōu)選用的原則,剔除相關(guān)性不好、數(shù)據(jù)不完整的因子,選取相關(guān)性好的因子,最終得到與汛期(5-9月)降水百分率相關(guān)性較好的38個(gè)因子(表2),與主汛期(6-8月)降水百分率相關(guān)性較好的28個(gè)因子(表3)。

    由表2可以看出,①所選取的相關(guān)性較好的因子相關(guān)性在0.300~0.404,相關(guān)性最好的因子為上年8月東太平洋副高北界位置,相關(guān)系數(shù)為0.404,其他因子的相關(guān)系數(shù)均小于0.4。②選取的39個(gè)因子(上年11月印度副高強(qiáng)度指數(shù)和上年11月印度副高面積指數(shù)相同,均為0.374,實(shí)為2個(gè)因子)包括33個(gè)大氣環(huán)流指數(shù)、5個(gè)海溫指數(shù)、1個(gè)收集下載的指數(shù),其中副高有22個(gè),海溫有6個(gè),可見副高與汛期(5-9月)降水相關(guān)性最大,海溫相關(guān)性也比較大。③從時(shí)間分布上來看,冬季有6個(gè)因子,上年8月有5個(gè),當(dāng)年2月、上年12月、夏季和全年各有4個(gè),當(dāng)年1月、上年3月和上年11月各有2個(gè),上年4、5、6、9月、春季、秋季各有1個(gè)。表明冬季的因子與汛期降水相關(guān)性最好,上年8月、當(dāng)年2月、上年12月、夏季和全年的因子與汛期降水相關(guān)性較好。

    由表3可以看出,①所選取的相關(guān)性較好的因子相關(guān)性在0.302~0.397,相關(guān)性最好的因子為上年11月印度副高強(qiáng)度指數(shù)和印度副高面積指數(shù),相關(guān)系數(shù)均為0.397,除冬季北非-北大西洋-北美副高脊線位置指數(shù)外,其他因子的相關(guān)系數(shù)均小于0.375。②選取的28個(gè)因子包括22個(gè)大氣環(huán)流指數(shù)、3個(gè)海溫指數(shù)和3個(gè)收集下載的指數(shù),其中副高有12個(gè),海溫有4個(gè),極渦位置、強(qiáng)度有4個(gè),可見副高與主汛期(6-8月)降水相關(guān)性最大,海溫與極渦位置、強(qiáng)度相關(guān)性也比較大。③從時(shí)間分布上來看,冬季有5個(gè)因子,上年8月和上年12月各有3個(gè),當(dāng)年2月、上年5月、上年6月、上年11月和全年各有2個(gè),當(dāng)年1月、上年3月、上年7月、上年9月、上年10月、夏季和秋季各有1個(gè)??梢哉J(rèn)為冬季的因子與主汛期降水相關(guān)性最好,其次是上年8月和上年12月的因子與主汛期降水相關(guān)性較好。

    3 模型建立及結(jié)果檢驗(yàn)

    利用篩選的相關(guān)性較好的39個(gè)因子(相關(guān)系數(shù)>0.3)和荊州站1954-2015年汛期(5-9月)降水距平百分率資料,采用逐步回歸統(tǒng)計(jì)方法,建立荊州站汛期降水距平百分率回歸預(yù)測方程,其中上年11月印度副高強(qiáng)度指數(shù)和上年11月印度副高面積指數(shù)與汛期降水距平百分率的相關(guān)性相同,且數(shù)據(jù)有高度的一致性,僅保留一項(xiàng)作為預(yù)測因子即可,最終建立的荊州汛期降水回歸預(yù)測模型(式1)包括38項(xiàng)氣候系統(tǒng)監(jiān)測指數(shù),模型相關(guān)系數(shù)為0.874,明顯大于單個(gè)因子的相關(guān)系數(shù)。用該汛期降水回歸模型來計(jì)算2016年汛期降水百分率,從而驗(yàn)證汛期降水預(yù)測模型的預(yù)測精度。荊州汛期降水距平百分率的模擬與驗(yàn)證結(jié)果見表4。由表4可以看出,1954-2015年62年模型擬合的結(jié)果中,52個(gè)年份能正確預(yù)測降水偏多偏少趨勢,10個(gè)年份預(yù)測錯(cuò)誤。對降水特多的年份1954、1973、1979、1980、1996年,僅1996年預(yù)測值屬偏多的等級,其他4年準(zhǔn)確預(yù)測降水特多的趨勢。對降水特少的年份1966年,模型預(yù)測結(jié)果為略少,有一定差距。荊州汛期降水預(yù)測模型預(yù)測2016年荊州汛期降水偏多17.3%,相比實(shí)際值6.2%,同屬于略多的等級。

    Y=231.28+4.864X1+0.006X2-0.932X3-4.281X4+3.566X5-13.734X6+8.619X7-12.082X8-1.265X9-0.858X10-0.199X11-8.394X12+9.945X13+2.296X14-2.829X15-0.478X16+4.811X17+7.028X18-4.41X19-42.553X20+44.247X21-0.463X22-1.004X23+3.317X24-5.209X25-3.837X26+3.363X27+14.373X28-4.695X29-4.904X30+0.857X31-0.011X32+17.582X33+10.146X34+25.677X35-5.344X36-19.272X37+0.119X38

    (R=0.874) (1)

    同理,利用1954-2015年主汛期降水資料,建立荊州站主汛期(6-8月)降水距平百分率逐步回歸預(yù)測模型(式2),該預(yù)測模型包含27項(xiàng)氣候系統(tǒng)監(jiān)測指數(shù),模型相關(guān)系數(shù)為0.914,明顯大于單個(gè)因子的相關(guān)系數(shù)。2016年主汛期降水資料用于模型驗(yàn)證,模型模擬與驗(yàn)證的結(jié)果見表5。由表5可以看出,1954-2015年62年間模型擬合的結(jié)果中,55個(gè)年份能正確預(yù)測降水偏多偏少趨勢,僅7個(gè)年份預(yù)測錯(cuò)誤。對降水特多的年份1954、1955、1964、1969、1979、1980、1983、1996年,1964、1983、1996年預(yù)測值屬偏多的等級,其他5年準(zhǔn)確預(yù)測降水特多的趨勢。對降水特少的年份1961、1966、1972、1974、1977年,1972、1974年模型準(zhǔn)確預(yù)測主汛期降水特少的趨勢,1961、1977年模型預(yù)測結(jié)果為偏少,1966年模型預(yù)測結(jié)果為略少。荊州主汛期降水預(yù)測模型預(yù)測2016年荊州主汛期降水偏多223.2%,與實(shí)際值30.2%相比,雖能預(yù)測2016年主汛期降水偏多,但數(shù)值明顯偏大。進(jìn)一步分析預(yù)測模型,發(fā)現(xiàn)模型中預(yù)測因子上年11月印度副高強(qiáng)度指數(shù)的有效數(shù)據(jù)較少,僅1953年和2015年的數(shù)值不為0,其他年份均為0,1954年汛期和主汛期降水均為有記錄以來最多的年份,但該因子在主汛期降水預(yù)測模型中,系數(shù)相比汛期降水預(yù)測模型大兩個(gè)數(shù)量級,從而導(dǎo)致計(jì)算出的2016年主汛期降水明顯偏大,可能需要更多數(shù)據(jù)來調(diào)整該項(xiàng)的系數(shù),從而提高預(yù)測精度。

    Y=753.699-3.97X1-4.541X2+3.237X3-10.339X4-0.02X5-0.007X6-8.368X7+8.799X8-41.907X9-2.679X10-0.061X11+1.149X12-30.323X13-0.019X14+363.333X15+0.014X16-17.318X17-2.51X18+22.87X19+15.147X20+3.609X21-4.323X22-2.866X23-0.027X24+15.59X25+0.763X26+0.063X27

    (R=0.914) (2)

    4 小結(jié)

    1)荊州汛期降水預(yù)測模型選取39個(gè)相關(guān)性>0.3的氣候系統(tǒng)監(jiān)測指數(shù)作為預(yù)測因子,其中副高與汛期(5-9月)降水相關(guān)性最大,海溫相關(guān)性也比較大,相關(guān)性最好的因子為上年8月東太平洋副高北界位置,相關(guān)系數(shù)為0.404,其他因子的相關(guān)系數(shù)均小于0.4。從時(shí)間上來看,冬季的因子與汛期降水相關(guān)性最好。

    2)荊州主汛期降水預(yù)測模型選取28個(gè)相關(guān)性>0.3的氣候系統(tǒng)監(jiān)測指數(shù)作為預(yù)測因子,其中副高與主汛期(6-8月)降水相關(guān)性最大,海溫、極渦位置和強(qiáng)度相關(guān)性也比較大。相關(guān)性最好的因子為上年11月印度副高強(qiáng)度指數(shù)和印度副高面積指數(shù),相關(guān)系數(shù)均為0.397,除冬季北非-北大西洋-北美副高脊線位置指數(shù)外,其他因子的相關(guān)系數(shù)均小于0.375。從時(shí)間上來看,冬季的因子與主汛期降水相關(guān)性最好。

    3)1954-2015年62年汛期降水資料擬合的荊州汛期降水預(yù)測模型,相關(guān)系數(shù)為0.874,52個(gè)年份能正確預(yù)測降水偏多偏少趨勢,10個(gè)年份預(yù)測錯(cuò)誤。該模型預(yù)測的2016年荊州汛期降水趨勢為偏多17.3%,相比實(shí)際值6.2%,同屬于略多的等級。1954-2015年62年主汛期降水資料擬合的荊州主汛期降水預(yù)測模型,相關(guān)系數(shù)高達(dá)0.914,55個(gè)年份能正確預(yù)測降水偏多偏少趨勢,僅7個(gè)年份預(yù)測錯(cuò)誤。該模型預(yù)測2016年荊州主汛期降水偏多223.2%,與實(shí)際值30.2%相比,雖能預(yù)測2016年主汛期降水偏多,但數(shù)值明顯偏大。進(jìn)一步分析預(yù)測模型,發(fā)現(xiàn)模型中預(yù)測因子上年11月印度副高強(qiáng)度指數(shù)的有效數(shù)據(jù)較少,可能需要更多數(shù)據(jù)來調(diào)整該項(xiàng)的系數(shù),提高模型預(yù)測精度。綜合分析可知,本研究所建立的荊州汛期和主汛期降水預(yù)測模型對荊州汛期和主汛期降水趨勢預(yù)測有較好的指導(dǎo)意義。

    參考文獻(xiàn):

    [1] 黃榮輝,蔡榕碩,陳際龍,等.我國旱澇氣候?yàn)?zāi)害的年代際變化及其與東亞氣候系統(tǒng)變化的關(guān)系[J].大氣科學(xué),2006,30(5):730-743.

    [2] 譚桂容,孫照渤,陳海山.華北夏季早澇的環(huán)流特征分析[J].氣象科學(xué),2003,23(2):135-143.

    [3] 韓照宇,劉 榮,王振華.山西省降水與環(huán)流特征量及ELnino的關(guān)系[J].山西氣象,2001(2):22-25.

    [4] 羅伯良,朱明輝.湖南省夏季旱澇預(yù)報(bào)模型研究及試報(bào)[J].氣象,1999,25(10):23-27.

    [5] 丑紀(jì)范.天氣數(shù)值預(yù)報(bào)中使用過去資料的問題[J].中國科學(xué),1974,

    1(6):635-644.

    [6] 丑紀(jì)范.短期氣候預(yù)測的現(xiàn)狀、問題與出路(一)[J].新疆氣象,2003,26(1):1-4.

    [7] 邱崇踐,丑紀(jì)范.天氣預(yù)報(bào)的相似-動(dòng)力方法[J].大氣科學(xué),1989,

    13(1):22-28.

    [8] 黃建平,王紹武.相似-動(dòng)力模式的季節(jié)預(yù)報(bào)試驗(yàn)[J].中國科學(xué)(B輯),1991,35(2):216-224.

    [9] 任宏利,丑紀(jì)范.統(tǒng)計(jì)-動(dòng)力相結(jié)合的相似誤差訂正法[J].氣象學(xué)報(bào),2005,63(6):988-993.

    [10] 鄭志海,任宏利,黃建平.基于季節(jié)氣候可預(yù)報(bào)分量的相似誤差訂正方法和數(shù)值試驗(yàn)[J].物理學(xué)報(bào),2009,58(10):7359-7367.

    [11] 王啟光,封國林,鄭志海,等.長江中下游汛期降水優(yōu)化多因子組合客觀定量化預(yù)測研究[J].大氣科學(xué),2011,35(2):287-297.

    [12] 楊 杰,王啟光,支 蓉,等.動(dòng)態(tài)最優(yōu)多因子組合的華北汛期降水模式誤差估計(jì)及預(yù)報(bào)[J].物理學(xué)報(bào),2011,60(2):837-839.

    [13] 楊 杰,趙俊虎,鄭志海,等.華北汛期降水多因子相似訂正方案與預(yù)報(bào)試驗(yàn)[J].大氣科學(xué),2012,36(1):11-22.

    [14] 百項(xiàng)氣候系統(tǒng)指數(shù)集[DB/OL].http://10.1.64.154/portal/web-link.index_inid=1001.

    [15] 百項(xiàng)氣候系統(tǒng)指數(shù)集[DB/OL].http://cmdp.ncc.cma.gov.cn/Monitoring/cn_index_130.php.

    猜你喜歡
    預(yù)測模型
    基于矩陣?yán)碚撓碌母咝=處熑藛T流動(dòng)趨勢預(yù)測
    東方教育(2016年9期)2017-01-17 21:04:14
    基于支持向量回歸的臺灣旅游短期客流量預(yù)測模型研究
    基于神經(jīng)網(wǎng)絡(luò)的北京市房價(jià)預(yù)測研究
    商情(2016年43期)2016-12-23 14:23:13
    中國石化J分公司油氣開發(fā)投資分析與預(yù)測模型研究
    基于IOWHA法的物流需求組合改善與預(yù)測模型構(gòu)建
    基于小波神經(jīng)網(wǎng)絡(luò)的GDP預(yù)測
    區(qū)域環(huán)境質(zhì)量全局評價(jià)預(yù)測模型及運(yùn)用
    中國市場(2016年30期)2016-07-18 05:15:34
    組合預(yù)測法在汽車預(yù)測中的應(yīng)用
    中國市場(2016年19期)2016-06-06 05:11:02
    H市電動(dòng)汽車充電站項(xiàng)目建設(shè)需求規(guī)模與技術(shù)分析
    商(2016年13期)2016-05-20 10:31:40
    基于預(yù)測模型加擾動(dòng)控制的最大功率點(diǎn)跟蹤研究
    科技視界(2016年1期)2016-03-30 13:37:45
    免费观看精品视频网站| 又大又黄又爽视频免费| av在线天堂中文字幕| 国产在视频线精品| 亚洲av成人精品一区久久| 亚洲经典国产精华液单| 日韩制服骚丝袜av| 性色avwww在线观看| 国产一区二区三区综合在线观看 | 三级经典国产精品| 男人舔女人下体高潮全视频| 少妇猛男粗大的猛烈进出视频 | 国产伦在线观看视频一区| 天天躁夜夜躁狠狠久久av| 一级毛片 在线播放| 久久97久久精品| 干丝袜人妻中文字幕| 国产色爽女视频免费观看| 日韩在线高清观看一区二区三区| 天天一区二区日本电影三级| 大陆偷拍与自拍| 纵有疾风起免费观看全集完整版 | 69人妻影院| 亚洲欧美清纯卡通| 在现免费观看毛片| 51国产日韩欧美| 三级国产精品片| 女的被弄到高潮叫床怎么办| 成人亚洲精品av一区二区| 中文字幕人妻熟人妻熟丝袜美| 欧美97在线视频| 麻豆成人av视频| 人体艺术视频欧美日本| 中文字幕av成人在线电影| 日韩中字成人| 一级片'在线观看视频| 久久久国产一区二区| 噜噜噜噜噜久久久久久91| 日日摸夜夜添夜夜爱| 免费观看的影片在线观看| 久久久久性生活片| 99久久人妻综合| 成人欧美大片| 欧美激情久久久久久爽电影| 亚洲av日韩在线播放| 婷婷色麻豆天堂久久| 国产成人91sexporn| 日韩强制内射视频| 日韩欧美精品免费久久| 亚洲图色成人| 免费观看av网站的网址| 精品久久久久久久久av| 久久久久性生活片| 精品酒店卫生间| 五月伊人婷婷丁香| 六月丁香七月| 午夜精品一区二区三区免费看| 全区人妻精品视频| 国产一区二区亚洲精品在线观看| 麻豆精品久久久久久蜜桃| 18+在线观看网站| 午夜日本视频在线| 久久99热这里只频精品6学生| 久久草成人影院| 久久精品国产亚洲av涩爱| 亚洲国产精品sss在线观看| 视频中文字幕在线观看| 别揉我奶头 嗯啊视频| 人人妻人人澡人人爽人人夜夜 | 国产乱来视频区| 久久亚洲国产成人精品v| 人人妻人人澡欧美一区二区| 国产黄色小视频在线观看| 一级av片app| 日日摸夜夜添夜夜添av毛片| eeuss影院久久| 日日摸夜夜添夜夜爱| 欧美成人午夜免费资源| 我的老师免费观看完整版| 午夜激情久久久久久久| 青春草国产在线视频| 亚洲性久久影院| 少妇人妻一区二区三区视频| 搡老妇女老女人老熟妇| 中文字幕人妻熟人妻熟丝袜美| 成年女人看的毛片在线观看| 一级毛片电影观看| 欧美成人精品欧美一级黄| 特级一级黄色大片| 最近最新中文字幕免费大全7| 高清欧美精品videossex| 国产成人精品一,二区| 身体一侧抽搐| 亚洲av一区综合| 女的被弄到高潮叫床怎么办| 国产成人精品久久久久久| 一级片'在线观看视频| 亚洲人与动物交配视频| 国产精品.久久久| 亚洲国产av新网站| 亚洲国产精品成人久久小说| 国产成人午夜福利电影在线观看| 又大又黄又爽视频免费| 午夜老司机福利剧场| 一级爰片在线观看| 国产真实伦视频高清在线观看| 国产高清有码在线观看视频| 国产 一区 欧美 日韩| 国产成人精品久久久久久| 男女边吃奶边做爰视频| 国产精品蜜桃在线观看| 国产免费福利视频在线观看| 最近视频中文字幕2019在线8| 亚洲精品中文字幕在线视频 | 国产综合精华液| 亚洲国产成人一精品久久久| 亚洲成色77777| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 51国产日韩欧美| 亚洲在线自拍视频| 亚洲国产色片| 插阴视频在线观看视频| 国产视频内射| 久久这里有精品视频免费| 寂寞人妻少妇视频99o| 国产一区二区三区av在线| 亚洲精品乱久久久久久| 亚洲国产欧美人成| 欧美另类一区| 18禁裸乳无遮挡免费网站照片| 国产精品一二三区在线看| 久久久久久久亚洲中文字幕| 午夜免费观看性视频| 亚洲色图av天堂| 中文欧美无线码| 久久久久久国产a免费观看| 一区二区三区四区激情视频| 非洲黑人性xxxx精品又粗又长| 成人综合一区亚洲| 九九在线视频观看精品| 国产三级在线视频| 97在线视频观看| 51国产日韩欧美| 联通29元200g的流量卡| 99视频精品全部免费 在线| 自拍偷自拍亚洲精品老妇| 国产精品人妻久久久影院| 丝袜美腿在线中文| 久久久久性生活片| 黄色欧美视频在线观看| 国产男女超爽视频在线观看| 黄片无遮挡物在线观看| 国产黄a三级三级三级人| 一级毛片我不卡| 99久国产av精品| 久久亚洲国产成人精品v| 一区二区三区高清视频在线| 免费电影在线观看免费观看| 色播亚洲综合网| 婷婷色av中文字幕| 人人妻人人看人人澡| 欧美精品一区二区大全| 亚洲精品久久午夜乱码| 亚洲丝袜综合中文字幕| 欧美bdsm另类| 99久久精品国产国产毛片| 少妇的逼水好多| 精品久久久久久久久亚洲| 性色avwww在线观看| 亚洲av中文字字幕乱码综合| 亚洲综合色惰| 国产一级毛片在线| 国产高潮美女av| 一级毛片久久久久久久久女| 亚洲av成人精品一区久久| 免费看a级黄色片| 国产成人a∨麻豆精品| 丝袜喷水一区| 国产真实伦视频高清在线观看| 性插视频无遮挡在线免费观看| 国产成人a区在线观看| 国产精品国产三级专区第一集| 一个人观看的视频www高清免费观看| 亚洲人成网站高清观看| 国产一区二区在线观看日韩| 免费黄色在线免费观看| 一级毛片黄色毛片免费观看视频| 国产亚洲最大av| 99热6这里只有精品| 国产亚洲最大av| 性插视频无遮挡在线免费观看| 麻豆av噜噜一区二区三区| 国产欧美日韩精品一区二区| 国产老妇女一区| 亚洲国产欧美人成| 国产在视频线在精品| 午夜福利在线在线| 亚洲四区av| 国产午夜精品论理片| 免费电影在线观看免费观看| 中文字幕亚洲精品专区| 国产精品一区二区性色av| 搡老妇女老女人老熟妇| 国产午夜精品一二区理论片| 97超视频在线观看视频| 国产精品蜜桃在线观看| 午夜福利在线在线| 最近手机中文字幕大全| 人妻制服诱惑在线中文字幕| 欧美xxⅹ黑人| 九色成人免费人妻av| 极品教师在线视频| 国产黄色免费在线视频| 精品欧美国产一区二区三| 国产精品99久久久久久久久| 有码 亚洲区| 亚洲久久久久久中文字幕| 中文字幕av成人在线电影| 欧美成人a在线观看| 欧美变态另类bdsm刘玥| 国产精品国产三级国产av玫瑰| 久久久久性生活片| 国产精品综合久久久久久久免费| 97热精品久久久久久| 亚洲精品一二三| 一级毛片 在线播放| 久久精品久久久久久久性| 国产熟女欧美一区二区| 亚洲国产精品成人综合色| 久久这里有精品视频免费| 久久精品夜夜夜夜夜久久蜜豆| 狂野欧美白嫩少妇大欣赏| 波野结衣二区三区在线| 91在线精品国自产拍蜜月| 精品一区二区三卡| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲婷婷狠狠爱综合网| 一个人看视频在线观看www免费| 777米奇影视久久| 晚上一个人看的免费电影| 国产女主播在线喷水免费视频网站 | 久久精品夜色国产| av在线蜜桃| 久久6这里有精品| av专区在线播放| 美女国产视频在线观看| 免费黄网站久久成人精品| 亚洲av中文av极速乱| 能在线免费看毛片的网站| 伦精品一区二区三区| 国产av不卡久久| 亚洲人成网站在线播| 亚洲精品视频女| 一区二区三区乱码不卡18| 内射极品少妇av片p| 国产男人的电影天堂91| 国产亚洲av片在线观看秒播厂 | 国产中年淑女户外野战色| 国产av码专区亚洲av| 天美传媒精品一区二区| av专区在线播放| 亚洲av男天堂| 九九爱精品视频在线观看| 国产精品国产三级专区第一集| 国内精品宾馆在线| 大陆偷拍与自拍| 日韩成人av中文字幕在线观看| 成人综合一区亚洲| 国产成年人精品一区二区| 午夜精品国产一区二区电影 | 亚洲无线观看免费| 亚洲熟妇中文字幕五十中出| 午夜久久久久精精品| 日韩欧美精品v在线| 亚洲成人久久爱视频| 在线a可以看的网站| 一个人免费在线观看电影| 成人一区二区视频在线观看| 亚洲18禁久久av| 亚洲人成网站在线播| 亚洲乱码一区二区免费版| 国语对白做爰xxxⅹ性视频网站| 国产成人freesex在线| 亚洲国产欧美人成| 精品久久久噜噜| 亚洲国产av新网站| 亚洲精品乱码久久久v下载方式| 亚洲成人精品中文字幕电影| 日本黄色片子视频| 91午夜精品亚洲一区二区三区| 国产在视频线精品| 国产av码专区亚洲av| 91久久精品国产一区二区成人| 日韩成人伦理影院| 噜噜噜噜噜久久久久久91| 日韩电影二区| 日韩av免费高清视频| 观看免费一级毛片| 亚洲av男天堂| 国产大屁股一区二区在线视频| 国产不卡一卡二| 国产精品国产三级国产专区5o| 国产高潮美女av| 99久久精品国产国产毛片| 国产精品1区2区在线观看.| 少妇裸体淫交视频免费看高清| a级一级毛片免费在线观看| 亚洲国产欧美在线一区| 男女那种视频在线观看| 亚洲伊人久久精品综合| 99久久人妻综合| eeuss影院久久| 99热6这里只有精品| av福利片在线观看| 久久久久久伊人网av| 草草在线视频免费看| 男人爽女人下面视频在线观看| 秋霞伦理黄片| 精华霜和精华液先用哪个| 亚洲综合精品二区| 精品久久久噜噜| 国产乱人视频| 青青草视频在线视频观看| 久久精品国产亚洲av天美| 国内精品美女久久久久久| 青春草亚洲视频在线观看| 内地一区二区视频在线| 我的老师免费观看完整版| 青青草视频在线视频观看| 国产亚洲av片在线观看秒播厂 | 国产 亚洲一区二区三区 | 亚洲av免费高清在线观看| 国产精品不卡视频一区二区| 在线观看一区二区三区| 在现免费观看毛片| 欧美人与善性xxx| 最近最新中文字幕免费大全7| 麻豆精品久久久久久蜜桃| 国产av码专区亚洲av| 噜噜噜噜噜久久久久久91| 一级爰片在线观看| 亚洲一级一片aⅴ在线观看| 极品少妇高潮喷水抽搐| 日韩欧美三级三区| 一级爰片在线观看| 日韩av免费高清视频| 亚洲三级黄色毛片| 国产精品日韩av在线免费观看| 欧美潮喷喷水| 久久久久久久午夜电影| 国产精品不卡视频一区二区| 国产高清国产精品国产三级 | 欧美97在线视频| 亚洲18禁久久av| 亚洲在线自拍视频| 床上黄色一级片| 大又大粗又爽又黄少妇毛片口| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 青春草视频在线免费观看| 日韩欧美三级三区| 日本一本二区三区精品| 一边亲一边摸免费视频| 赤兔流量卡办理| 亚洲欧洲日产国产| 亚洲欧美成人精品一区二区| 99九九线精品视频在线观看视频| 内射极品少妇av片p| 天美传媒精品一区二区| 免费看光身美女| 少妇裸体淫交视频免费看高清| 国产伦一二天堂av在线观看| 国产av国产精品国产| 欧美激情久久久久久爽电影| 国产精品久久视频播放| 全区人妻精品视频| 麻豆成人午夜福利视频| 日韩一本色道免费dvd| 国产美女午夜福利| 1000部很黄的大片| 亚洲高清免费不卡视频| 极品少妇高潮喷水抽搐| 国产视频内射| 亚洲第一区二区三区不卡| 少妇猛男粗大的猛烈进出视频 | 欧美激情国产日韩精品一区| 国产精品三级大全| 尤物成人国产欧美一区二区三区| 精华霜和精华液先用哪个| 国产精品国产三级国产专区5o| 午夜福利网站1000一区二区三区| 亚洲婷婷狠狠爱综合网| 赤兔流量卡办理| 亚洲18禁久久av| 国产亚洲av片在线观看秒播厂 | 男插女下体视频免费在线播放| 国产精品.久久久| 听说在线观看完整版免费高清| 真实男女啪啪啪动态图| 少妇丰满av| 天堂俺去俺来也www色官网 | 精品久久久久久久久av| 亚洲人成网站高清观看| 国产精品日韩av在线免费观看| 人人妻人人看人人澡| 国产精品女同一区二区软件| 99久久九九国产精品国产免费| 六月丁香七月| 国内精品美女久久久久久| 简卡轻食公司| 亚洲成人精品中文字幕电影| 特大巨黑吊av在线直播| 一级毛片 在线播放| 亚洲在线观看片| 一级a做视频免费观看| 国产精品人妻久久久久久| 综合色丁香网| 免费av不卡在线播放| 最后的刺客免费高清国语| 亚洲欧美一区二区三区黑人 | 久久精品国产自在天天线| 一级毛片 在线播放| 国产91av在线免费观看| 97超视频在线观看视频| 国产在线男女| 少妇裸体淫交视频免费看高清| 亚洲精品久久午夜乱码| 伦理电影大哥的女人| 91精品伊人久久大香线蕉| 国产老妇女一区| 五月天丁香电影| 男人舔奶头视频| 亚洲国产精品成人久久小说| 免费无遮挡裸体视频| 51国产日韩欧美| 国产毛片a区久久久久| 激情五月婷婷亚洲| 亚洲精品一区蜜桃| 蜜臀久久99精品久久宅男| 亚洲国产最新在线播放| 亚洲欧美成人精品一区二区| 国产黄a三级三级三级人| 欧美变态另类bdsm刘玥| 最近中文字幕高清免费大全6| 成人亚洲精品一区在线观看 | 少妇丰满av| 精品久久久久久久久久久久久| 亚洲精品色激情综合| 成人午夜精彩视频在线观看| 18禁在线播放成人免费| 国产精品一区二区在线观看99 | 成年版毛片免费区| 午夜精品在线福利| av在线观看视频网站免费| 男人狂女人下面高潮的视频| 国产伦精品一区二区三区视频9| 日本一本二区三区精品| 22中文网久久字幕| 国产极品天堂在线| 爱豆传媒免费全集在线观看| 一级黄片播放器| 两个人的视频大全免费| 成年人午夜在线观看视频 | 丝瓜视频免费看黄片| 久久精品人妻少妇| av女优亚洲男人天堂| 高清午夜精品一区二区三区| 午夜老司机福利剧场| 99re6热这里在线精品视频| 深爱激情五月婷婷| 人妻一区二区av| 免费大片18禁| 亚洲在久久综合| 久久99热这里只频精品6学生| 国产伦在线观看视频一区| 国产精品精品国产色婷婷| 国产激情偷乱视频一区二区| 国产一区二区亚洲精品在线观看| 国产午夜精品一二区理论片| 久久99热这里只有精品18| av黄色大香蕉| 精品久久久久久电影网| 亚洲美女搞黄在线观看| 精品人妻偷拍中文字幕| 久久午夜福利片| 欧美不卡视频在线免费观看| 国产精品精品国产色婷婷| 人人妻人人看人人澡| 青春草亚洲视频在线观看| 亚洲色图av天堂| 国产老妇伦熟女老妇高清| 2021天堂中文幕一二区在线观| 久久精品国产亚洲av天美| 少妇的逼好多水| 欧美激情国产日韩精品一区| 国产精品久久久久久av不卡| 黄色一级大片看看| 啦啦啦中文免费视频观看日本| 黄色一级大片看看| 精华霜和精华液先用哪个| 超碰av人人做人人爽久久| 天堂av国产一区二区熟女人妻| 久久久精品欧美日韩精品| 国产伦理片在线播放av一区| 只有这里有精品99| 韩国av在线不卡| 国产午夜精品久久久久久一区二区三区| 非洲黑人性xxxx精品又粗又长| 国产精品福利在线免费观看| 国产视频内射| 哪个播放器可以免费观看大片| 国产精品伦人一区二区| 日韩中字成人| 亚洲精品乱码久久久v下载方式| 1000部很黄的大片| 国产精品日韩av在线免费观看| 亚洲欧美成人精品一区二区| 69人妻影院| 国产男女超爽视频在线观看| 最近手机中文字幕大全| 欧美日韩一区二区视频在线观看视频在线 | 国产伦理片在线播放av一区| 热99在线观看视频| 久久久久久久国产电影| 日韩伦理黄色片| 色视频www国产| 人人妻人人澡欧美一区二区| 最近的中文字幕免费完整| 婷婷色综合大香蕉| 啦啦啦韩国在线观看视频| 国精品久久久久久国模美| 成人毛片60女人毛片免费| 欧美 日韩 精品 国产| 亚洲精品色激情综合| 欧美激情国产日韩精品一区| 大香蕉97超碰在线| 亚洲精品日韩av片在线观看| 欧美97在线视频| 国精品久久久久久国模美| 亚洲精品aⅴ在线观看| 成人亚洲精品av一区二区| 国产亚洲91精品色在线| 大香蕉久久网| 国产精品一区www在线观看| eeuss影院久久| 成人欧美大片| 大片免费播放器 马上看| 亚洲av中文av极速乱| 国产在视频线在精品| 18禁在线播放成人免费| 精品熟女少妇av免费看| 国产精品久久久久久久久免| 日韩av免费高清视频| 成年版毛片免费区| 亚洲精品第二区| 亚洲国产精品专区欧美| 日韩成人伦理影院| 午夜精品在线福利| 男人舔女人下体高潮全视频| 赤兔流量卡办理| 成人无遮挡网站| 亚洲乱码一区二区免费版| 国产成人a区在线观看| 成人午夜高清在线视频| 国产精品综合久久久久久久免费| 日本三级黄在线观看| 搡老乐熟女国产| 欧美三级亚洲精品| 久久久久久久大尺度免费视频| 国产成人福利小说| av天堂中文字幕网| 国产69精品久久久久777片| 亚洲av福利一区| 女人久久www免费人成看片| 在线免费观看不下载黄p国产| 七月丁香在线播放| 国产精品99久久久久久久久| av女优亚洲男人天堂| 亚洲av成人精品一二三区| 亚洲精华国产精华液的使用体验| 美女cb高潮喷水在线观看| 五月伊人婷婷丁香| 国产午夜精品一二区理论片| 最近手机中文字幕大全| 五月伊人婷婷丁香| 看黄色毛片网站| 男插女下体视频免费在线播放| 精华霜和精华液先用哪个| 成人午夜精彩视频在线观看| 亚洲va在线va天堂va国产| 亚洲国产精品sss在线观看| 成人欧美大片| 一区二区三区乱码不卡18| 97人妻精品一区二区三区麻豆| 99久久精品一区二区三区| 美女内射精品一级片tv| 91aial.com中文字幕在线观看| 成人美女网站在线观看视频| 毛片一级片免费看久久久久| 91午夜精品亚洲一区二区三区| 老女人水多毛片| 69av精品久久久久久| 男人舔奶头视频| 亚洲在线自拍视频| 毛片一级片免费看久久久久| 久久99蜜桃精品久久| 久久精品熟女亚洲av麻豆精品 | 亚洲欧美日韩东京热| 国产美女午夜福利| 亚洲一级一片aⅴ在线观看| 国产成人精品福利久久| 国产精品国产三级国产专区5o| 又粗又硬又长又爽又黄的视频| 国产男女超爽视频在线观看| 国产一级毛片在线| 成人特级av手机在线观看|