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

    基于LSTM神經(jīng)網(wǎng)絡(luò)的烏魯木齊市流感樣病例的預(yù)測(cè)研究

    2019-12-02 02:34:16龔風(fēng)云王凱
    科技視界 2019年31期
    關(guān)鍵詞:氣象因素預(yù)測(cè)

    龔風(fēng)云 王凱

    【摘 要】目的:分析烏魯木齊市流感樣病例月發(fā)病數(shù)的變化趨勢(shì),建立長(zhǎng)短期記憶(LSTM)模型,對(duì)流感樣病例例數(shù)數(shù)進(jìn)行預(yù)測(cè),為烏魯木齊市流感的預(yù)防與控制提供科學(xué)依據(jù)。方法:利用2015年1月-2018年3月烏魯木齊市的每月氣象數(shù)據(jù)、流感樣病例監(jiān)測(cè)數(shù)據(jù),采用單變量LSTM模型和多變量LSTM模型對(duì)烏魯木齊市流感樣病例例數(shù)的時(shí)間序列進(jìn)行預(yù)測(cè),使用RMSE和MAE值評(píng)價(jià)不同方法的預(yù)測(cè)精度。結(jié)果:?jiǎn)巫兞縇STM模型和多變量LSTM模型的RMSE值分別是66.17和56.91;MAE值分別是60.42和39.07。與單變量LSTM模型相比,多變量的LSTM模型預(yù)測(cè)效果較好。結(jié)論:本研究所建立的多變量LSTM模型能較好地預(yù)測(cè)ILI病例數(shù)的發(fā)病趨勢(shì),為流感監(jiān)測(cè)和預(yù)防控制提供依據(jù)。

    【關(guān)鍵詞】LSTM;流感樣病例;氣象因素;預(yù)測(cè)

    中圖分類號(hào): TP393.0;TP183 文獻(xiàn)標(biāo)識(shí)碼: A 文章編號(hào): 2095-2457(2019)31-0020-003

    DOI:10.19694/j.cnki.issn2095-2457.2019.31.009

    Prediction of influenza-like cases in urumqi based on LSTM neural network

    GONG Feng-yun1 WANG Kai2*

    (1.College of Applied Mathematics,Xinjiang University of Finance and Economis, Urumqi Xinjiang 830012, China;

    2.College of Medical Engineering and Technology,Xinjiang Medical University, Urumqi Xinjiang 830011, China)

    【Abstract】[Objective]To analyze the change trend of the monthly incidence of influenza-like cases in Urumqi, establish long-short term memory(LSTM) model,and predict the number of influenza-like cases, so as to provide scientific basis for the prevention and control of influenza in Urumqi.[Methods]Monthly meteorological data and influenza-like case monitoring data from January 2015 to March 2018 in Urumqi were used to predict the time series of influenza-like cases in Urumqi using univariate LSTM model and multivariate LSTMmodel, RMSE and MAE values were used to evaluate the prediction accuracy of different methods.[Results]RMSE values of univariate LSTM model and multivariate LSTM model are 66.17 and 56.91 respectively.The MAE values are 60.42 and 39.07 respectively.Compared with the single-variable LSTM model, the multivariable LSTM model has a better predictive effect. [Conclusion] The multivariate LSTM model established in this study can better predict the incidence trend of ILI cases, providing a basis for influenza surveillance, preventionand control.

    【Key words】LSTM; Influenza-like cases; Meteorological factors; Prediction

    流行性感冒簡(jiǎn)稱流感,是由流感病毒引起的一種傳染性呼吸道疾病[1]。流感病毒主要通過空氣傳播,如咳嗽、噴嚏等[2]。據(jù)統(tǒng)計(jì),全球每年約有10%~20%的人群因流感的季節(jié)性流行而感染流感,其中300萬(wàn)~500萬(wàn)人為重癥病例,25萬(wàn)~50萬(wàn)人為死亡病例[3]。流感的流行不僅對(duì)人群健康造成嚴(yán)重威脅,也對(duì)社會(huì)產(chǎn)生了巨大的經(jīng)濟(jì)損失[4]。在美國(guó),流感每年的總成本超過100億美元,而未來(lái)的流感大流行估計(jì)直接和間接成本高達(dá)數(shù)千億美元[5]。Yang J等[6]一項(xiàng)覆蓋全國(guó)范圍的研究顯示,個(gè)人流感門診病例的平均費(fèi)用為$155,流感住院病例的平均費(fèi)用為$1511。在疾病高峰期,診所和醫(yī)院不堪重負(fù)。接種流感疫苗和阻礙傳播途徑是預(yù)防流感有效的方法[7]。新疆烏魯木齊市是流感多發(fā)地之一,該市作為國(guó)家級(jí)流感監(jiān)測(cè)點(diǎn),對(duì)流感樣病例(influenza-like illness,ILI)進(jìn)行了長(zhǎng)期持續(xù)的監(jiān)測(cè)[8]。為了幫助政府、醫(yī)院、診所、制藥公司和其他公司有效地控制流感爆發(fā),并及時(shí)限制傳播途徑,本研究使用LSTM模型對(duì)烏魯木齊市ILI例數(shù)數(shù)據(jù)進(jìn)行訓(xùn)練擬合、分析及預(yù)測(cè),并比較預(yù)測(cè)精度。

    1 資料與方法

    1.1? ILI定義

    發(fā)熱(體溫≥38℃)、伴咳嗽或咽痛之一者[9]。

    1.2 資料來(lái)源

    2015年1月—2018年3月烏魯木齊市的每月ILI例數(shù)來(lái)源于“中國(guó)疾病預(yù)防控制信息系統(tǒng)流感監(jiān)測(cè)信息子系統(tǒng)”。氣象數(shù)據(jù)資料來(lái)自烏魯木齊市氣象局,主要包括月平均氣溫(℃)、月降水量(mm)、月平均氣壓(hpa)、月平均最大氣壓(hpa)、月平均最小氣壓(hpa)、月平均相對(duì)濕度(%)、月平均風(fēng)速(m/s)、月日照時(shí)數(shù)(h)。

    1.3 長(zhǎng)短期記憶人工神經(jīng)網(wǎng)絡(luò)(long-short term memory,LSTM)

    LSTM是一種改進(jìn)的時(shí)間循環(huán)神經(jīng)網(wǎng)絡(luò)(recurrent neural network,RNN)[10],它的提出解決了傳統(tǒng)循環(huán)神經(jīng)網(wǎng)絡(luò)在學(xué)習(xí)過程中由于輸入序列過長(zhǎng)帶來(lái)的梯度消失問題[11]。一個(gè)LSTM單元包括細(xì)胞狀態(tài)(cell state)、輸入門(input gate)、遺忘門(forget gate)、輸出門(forget gate),且LSTM神經(jīng)元的三個(gè)門的激活函數(shù)均為 Sigmoid[12]。

    1.4 模型預(yù)測(cè)效果的比較

    本文通過計(jì)算測(cè)試集的均方根誤差(Root Mean Square Error,RMSE)和平均絕對(duì)誤差(Mean Absolute Error,MAE)來(lái)定量地評(píng)價(jià)模型擬合及預(yù)測(cè)效果。該計(jì)算公式如下:

    RMSE=■

    MAE=■■|Y■-Y■|

    式中,Y1t和Y2t均分別為t時(shí)刻的觀測(cè)值和模型輸出值,n為數(shù)據(jù)點(diǎn)個(gè)數(shù)。

    1.5 統(tǒng)計(jì)分析

    使用Python3.0和R3.6.1軟件建立LSTM模型,其中各變量之間的相關(guān)性分析均采用了統(tǒng)計(jì)學(xué)中的Spearman相關(guān)。檢驗(yàn)水準(zhǔn)α=0.05。

    2 結(jié)果

    2.1 描述性分析

    2015年1月-2018年3月烏魯木齊市的ILI病例總數(shù)為5572例,月平均發(fā)病數(shù)約為143例。其中,月ILI例數(shù)、月平均氣溫(℃)、月降水量(mm)、月平均氣壓(hpa)、月平均最大氣壓(hpa)、月平均最小氣壓(hpa)、月平均相對(duì)濕度(%)、月平均風(fēng)速(m/s)、月日照時(shí)數(shù)(h)分別用fre、temp、rain、press、press_max、press_min、humidity、wind_speed、sunshine_hours表示。該時(shí)間序列圖顯示,流感樣病例數(shù)呈現(xiàn)明顯的季節(jié)性,冬春季發(fā)病數(shù)明顯多于夏秋季。月平均氣溫、月降水量、月平均風(fēng)速和月日照時(shí)數(shù)在冬春季偏低,夏秋季偏高,而月平均氣壓、月平均最大氣壓、月平均最小氣壓和月平均相對(duì)濕度在冬春季偏高,夏秋季偏低。見圖1。2015年1月—2018年3月烏魯木齊市的月平均氣溫、月降水量、月平均氣壓、月平均最大氣壓、月平均最小氣壓、月平均相對(duì)濕度、月平均風(fēng)速、月日照時(shí)數(shù)的均值分別為6.9℃、29.3mm、912.5hPa、917.6hpa、906.7hpa、59.7%、2.0m/s、222.4h。見表1。

    2.2 LSTM模型分析

    2.2.1 單變量LSTM模型

    由于LSTM不需要對(duì)時(shí)間序列進(jìn)行平穩(wěn)性和白噪聲檢驗(yàn),可以直接建立模型。所以,本文首先對(duì)數(shù)據(jù)進(jìn)行預(yù)處理,將時(shí)間序列數(shù)據(jù)轉(zhuǎn)化為學(xué)習(xí)數(shù)據(jù),并且對(duì)數(shù)據(jù)進(jìn)行歸一化處理。然后將數(shù)據(jù)集分為訓(xùn)練集和測(cè)試集,將2015年1月-2016年1月的數(shù)據(jù)作為訓(xùn)練集,2016年2月-2018年3月的數(shù)據(jù)作為測(cè)試集。最后建立LSTM模型,本文設(shè)置模型的隱藏層數(shù)為4,輸入層和輸出層分別為1和1;時(shí)間步長(zhǎng)(time_step)設(shè)置為1。

    圖1 2015年1月-2018年3月烏魯木齊市的每月ILI例數(shù)和氣象因素分布

    2.2.2 多變量LSTM模型

    由于單變量LSTM的預(yù)測(cè)結(jié)果并不太好,所以本文對(duì)單變量LSTM的輸入層與隱藏層進(jìn)行了改變,進(jìn)而構(gòu)建了一個(gè)可以輸入多個(gè)序列的多變量LSTM模型。本文設(shè)置的多變量LSTM模型的隱藏層數(shù)為50;輸入層和輸出層分別為8和1;每批次訓(xùn)練樣本數(shù)(batch_size)設(shè)置為72;時(shí)間步長(zhǎng)(time_step)設(shè)置為1;在 輸入特征時(shí)需要將tensor轉(zhuǎn)成三維作為L(zhǎng)STM cell的輸入。結(jié)果見圖2。

    圖2 多變量LSTM預(yù)測(cè)結(jié)果

    2.2.3 模型評(píng)價(jià)

    從表2可以看出,與單變量LSTM模型相比,納入多個(gè)變量的LSTM模型預(yù)測(cè)的RMSE和MAE略小。其中,多變量的LSTM模型預(yù)測(cè)的RMSE降低了9.26,預(yù)測(cè)的MAE降低了21.35。表明多變量LSTM模型可以提高模型的預(yù)測(cè)精度。

    表2 單變量LSTM模型和多變量LSTM模型的預(yù)測(cè)精度比較

    3 討論

    本研究采用時(shí)間序列分析方法,分析了2015年 1月~2018年3月新疆烏魯木齊市ILI例數(shù)的發(fā)病趨勢(shì)。時(shí)間序列模型通過使用不同的特性可以分為3種類型[5]。第一類模型是自回歸模型,它使用過去的患者數(shù)量作為特征(“xs”),并預(yù)測(cè)未來(lái)的患者數(shù)量作為響應(yīng)(y)。典型的例子包括自回歸綜合移動(dòng)平均(ARIMA)模型和向量自回歸模型(VAR)。第二類模型使用其他參數(shù)(如溫度、濕度等)代替過去的流感數(shù)據(jù)作為回歸模型(如線性回歸、隨機(jī)森林等)的特征。著名的例子是“谷歌流感趨勢(shì)”,它使用搜索引擎查詢數(shù)據(jù)作為特征和線性回歸模型。第三類模型是第一類和第二類的組合。它使用過去流感患者的數(shù)量作為特征(如第一種類型)和回歸模型(如第二種類型)。本研究采用了第三種模型類型,長(zhǎng)短期記憶(LSTM)模型來(lái)預(yù)測(cè)流感的爆發(fā)。

    從描述性分析結(jié)果來(lái)看,2015年1月—2018年3月烏魯木齊市的ILI發(fā)病數(shù)呈現(xiàn)冬春季高發(fā)的特點(diǎn),該特點(diǎn)不同于南方地區(qū)流感流行趨勢(shì),南方春季或初夏流感活動(dòng)較為活躍,這與南北方氣候差異相關(guān)。已有研究表明,地理位置、氣候條件的差異是造成北方地區(qū)與南方地區(qū)流感流行時(shí)間不同的主要因素[13]。由于烏魯木齊市位于中國(guó)西北地區(qū),是溫帶大陸性氣候[14]。溫帶大陸性氣候特點(diǎn)主要包括冬夏2季時(shí)間較長(zhǎng),春秋2季時(shí)間較短,且冬季寒冷,夏季炎熱,年降水量少,比較干旱[15]。所以其發(fā)病呈冬春季高發(fā)的特點(diǎn)。

    從LSTM模型結(jié)果分析來(lái)看,本文將氣象因素中的多變量LSTM模型概念引入到ILI發(fā)病數(shù)的預(yù)測(cè)中,為一般基于單變量 LSTM流感樣病例預(yù)測(cè)模型加入了更多氣象因素特征,不再僅僅將ILI例數(shù)看作為一個(gè)簡(jiǎn)單的數(shù)學(xué)概念上的時(shí)間序列。適當(dāng)使用多變量LSTM模型可能有助于預(yù)測(cè)目前和近期的流感傳播。由于本文現(xiàn)有的流感數(shù)據(jù)有限,所以,流感預(yù)測(cè)模型的準(zhǔn)確性還有待提高。但是,與單變量LSTM模型相比,納入多個(gè)變量的LSTM模型預(yù)測(cè)的RMSE和MAE略小,說明多變量LSTM模型可以用于ILI發(fā)病數(shù)的預(yù)測(cè),同時(shí)為公共衛(wèi)生人員提供一個(gè)相對(duì)準(zhǔn)確的參考點(diǎn),進(jìn)而對(duì)預(yù)防和控制流感的流行提出相關(guān)性建議。

    然而,本文的研究也有一些不足。預(yù)測(cè)模型的準(zhǔn)確性還有待提高,造成這個(gè)結(jié)果的可能原因是數(shù)據(jù)量較少,研究下一步將擴(kuò)大樣本量,再次對(duì)比單變量LSTM模型與多變量LSTM在預(yù)測(cè)疾病發(fā)病數(shù)上的精度。盡管多變量LSTM模型在預(yù)測(cè)精度上沒有明顯的優(yōu)勢(shì),但是多變量LSTM模型也為預(yù)測(cè)ILI例數(shù)提供了新的方法。

    【參考文獻(xiàn)】

    [1]譚婭文,萬(wàn)海同,何昱,杜海霞,楊潔紅,彭學(xué)謙,周惠芬.中藥抗流感病毒的作用及機(jī)制研究進(jìn)展[J].中國(guó)現(xiàn)代應(yīng)用藥學(xué),2019(16):2095-2099.

    [2]WHO,“Influenza(Seasonal),”http://www.who.int/mediacentre/factsheets/fs211/en/index.html.

    [3]崔永梅,杜中強(qiáng),侯瑩.孕婦感染甲型H3N2流感死亡病例的病因?qū)W分析[J].中國(guó)藥物與臨床,2019,19(5):834-835.

    [4]李文娟,王大燕.我國(guó)流感疾病負(fù)擔(dān)相關(guān)研究進(jìn)展[J/OL].中國(guó)人獸共患病學(xué)報(bào):1-6[2019-08-31].http://kns.cnki.net/kcms/detail/35.1284.R.20190723.1447.004.html.

    [5]Jie Z , Kazumitsu N . A comparative study on predicting influenza outbreaks[J]. BioScience Trends, 2017, 11(5):533-541.

    [6]Yang J,Jit M,Leung KS, et al.The economic burden of influenza-associated outpatient visits and hospitalizations in China: a retrospective survey [J]. Infect Dis Poverty,2015,4:44.

    [7]張慕麗,彭質(zhì)斌,鄭建東, 等.中國(guó)兒童流感疾病負(fù)擔(dān)和疫苗應(yīng)用現(xiàn)狀[J].中華實(shí)用兒科臨床雜志,2019,34(2):91-97.

    [8]高楓,阿不都熱依木,樊旭成,等.2016-2017年烏魯木齊市流感監(jiān)測(cè)結(jié)果分析[J].醫(yī)學(xué)信息,2018,31(19):134-136.

    [9]中華人民共和國(guó)衛(wèi)生部.全國(guó)流感監(jiān)測(cè)方案(2010年版)[J].國(guó)際呼吸雜志,2011,31(2):85-88.

    [10]張春露.基于Tensorflow的LSTM在太原空氣質(zhì)量AQI指數(shù)中的分析與預(yù)測(cè)[D].中北大學(xué),2019.

    [11]裴大衛(wèi),朱明.基于多因子與多變量長(zhǎng)短期記憶網(wǎng)絡(luò)的股票價(jià)格預(yù)測(cè)[J].計(jì)算機(jī)系統(tǒng)應(yīng)用,2019,28(8):30-38.

    [12]楊意豪,王梅,左銘.基于深度LSTM的甲亢疾病發(fā)展預(yù)測(cè)及應(yīng)用系統(tǒng)[J].智能計(jì)算機(jī)與應(yīng)用,2019,9(4):128-131.

    [13]蒲玉嬌. 2012一2017年烏魯木齊市流感流行特征與病原學(xué)監(jiān)測(cè)分析[D].新疆:新疆醫(yī)科大學(xué),2018:36-37.

    [14]LI J , LI XM . Response of stomatal conductance of two tree species to vapor pressure deficit in three climate zones[J]. J Arid Land, 2014, 6(6):771-781.

    [15]陶燕,顧天毅,王硯,等.蘭州市城關(guān)區(qū)流行性感冒與氣象因素的時(shí)間序列[J].蘭州大學(xué)學(xué)報(bào)(自然科學(xué)版),2018,54(1):137-142.

    猜你喜歡
    氣象因素預(yù)測(cè)
    無(wú)可預(yù)測(cè)
    黃河之聲(2022年10期)2022-09-27 13:59:46
    選修2-2期中考試預(yù)測(cè)卷(A卷)
    選修2-2期中考試預(yù)測(cè)卷(B卷)
    選修2—2期中考試預(yù)測(cè)卷(A卷)
    不可預(yù)測(cè)
    不必預(yù)測(cè)未來(lái),只需把握現(xiàn)在
    氣象因素對(duì)油菜品種蘇油211生長(zhǎng)發(fā)育和產(chǎn)量的影響
    氣象因素對(duì)油菜品種蘇油211生長(zhǎng)發(fā)育和產(chǎn)量的影響
    中國(guó)科技博覽(2016年1期)2016-04-25 11:41:02
    淺談氣象與農(nóng)業(yè)生產(chǎn)及自然災(zāi)害的關(guān)系
    久久影院123| 国产成人一区二区在线| 2022亚洲国产成人精品| 午夜久久久在线观看| 午夜91福利影院| 飞空精品影院首页| 99re6热这里在线精品视频| 男女啪啪激烈高潮av片| 国产精品免费大片| 欧美国产精品va在线观看不卡| 亚洲成人一二三区av| 久久毛片免费看一区二区三区| 婷婷成人精品国产| 91精品三级在线观看| 日韩av不卡免费在线播放| 日韩中文字幕视频在线看片| 一本久久精品| 男女午夜视频在线观看 | 国产亚洲午夜精品一区二区久久| 熟女人妻精品中文字幕| 国产不卡av网站在线观看| 老司机影院成人| av在线app专区| av有码第一页| 免费观看a级毛片全部| 伦理电影免费视频| 黄色视频在线播放观看不卡| 男人添女人高潮全过程视频| 99精国产麻豆久久婷婷| 国产欧美日韩综合在线一区二区| 久久国产亚洲av麻豆专区| 在线观看www视频免费| 亚洲欧美成人精品一区二区| 男女啪啪激烈高潮av片| 精品第一国产精品| 麻豆精品久久久久久蜜桃| 考比视频在线观看| 亚洲精品,欧美精品| 欧美另类一区| 韩国av在线不卡| 精品人妻在线不人妻| 亚洲国产精品一区二区三区在线| 国产男女内射视频| 亚洲欧美一区二区三区国产| 18禁国产床啪视频网站| 成人手机av| 成人影院久久| 亚洲精品中文字幕在线视频| 国产一区二区激情短视频 | 一区二区三区四区激情视频| 国产69精品久久久久777片| 亚洲精品一二三| 丝袜人妻中文字幕| 啦啦啦中文免费视频观看日本| 丰满饥渴人妻一区二区三| 少妇的逼水好多| 母亲3免费完整高清在线观看 | 天堂8中文在线网| 建设人人有责人人尽责人人享有的| 日韩欧美一区视频在线观看| 成人毛片60女人毛片免费| 日韩中文字幕视频在线看片| 国产 精品1| 亚洲精品久久成人aⅴ小说| 亚洲av欧美aⅴ国产| 大片电影免费在线观看免费| 黄色视频在线播放观看不卡| kizo精华| 汤姆久久久久久久影院中文字幕| 欧美精品av麻豆av| 欧美丝袜亚洲另类| 人人妻人人澡人人爽人人夜夜| 欧美bdsm另类| 精品久久久精品久久久| 在线观看免费视频网站a站| 在线看a的网站| 国产男人的电影天堂91| 精品一区二区三区视频在线| 日韩精品免费视频一区二区三区 | 咕卡用的链子| 午夜福利视频在线观看免费| 人人妻人人澡人人爽人人夜夜| 99热网站在线观看| 人妻少妇偷人精品九色| 精品人妻熟女毛片av久久网站| 黄片无遮挡物在线观看| 国产成人a∨麻豆精品| 熟妇人妻不卡中文字幕| 在线天堂最新版资源| 国产成人免费观看mmmm| 精品一区二区三区视频在线| 十八禁网站网址无遮挡| 日本vs欧美在线观看视频| 丰满少妇做爰视频| 丝袜喷水一区| 亚洲一区二区三区欧美精品| 免费观看无遮挡的男女| 国产精品麻豆人妻色哟哟久久| 欧美xxⅹ黑人| 9色porny在线观看| 男人舔女人的私密视频| 最近最新中文字幕大全免费视频 | 精品少妇黑人巨大在线播放| 少妇人妻久久综合中文| a级毛片黄视频| 午夜免费观看性视频| av网站免费在线观看视频| 黑人猛操日本美女一级片| a级毛色黄片| 亚洲久久久国产精品| 国产熟女午夜一区二区三区| 亚洲综合精品二区| 日韩欧美一区视频在线观看| 国产欧美日韩综合在线一区二区| 高清欧美精品videossex| 国产免费又黄又爽又色| 国产精品国产三级专区第一集| 中国国产av一级| 国产精品蜜桃在线观看| 亚洲av电影在线进入| 亚洲欧美成人综合另类久久久| 91精品国产国语对白视频| av线在线观看网站| 大话2 男鬼变身卡| 久久精品久久精品一区二区三区| 少妇精品久久久久久久| 91午夜精品亚洲一区二区三区| 国产亚洲欧美精品永久| 亚洲成人一二三区av| 久久久国产精品麻豆| 日韩大片免费观看网站| 日韩一区二区三区影片| 插逼视频在线观看| 国产精品欧美亚洲77777| 国产精品一二三区在线看| 不卡视频在线观看欧美| 精品一区二区三区四区五区乱码 | 性高湖久久久久久久久免费观看| 成人国语在线视频| 性色av一级| 成人亚洲精品一区在线观看| 久久久久久久大尺度免费视频| 看免费成人av毛片| 中文精品一卡2卡3卡4更新| 少妇熟女欧美另类| 国产成人精品在线电影| 日本猛色少妇xxxxx猛交久久| 亚洲精品第二区| 亚洲精品美女久久久久99蜜臀 | 亚洲三级黄色毛片| 日本色播在线视频| 亚洲国产最新在线播放| 久久久欧美国产精品| 国产一区二区在线观看日韩| 亚洲国产色片| 亚洲成av片中文字幕在线观看 | 大片电影免费在线观看免费| 纯流量卡能插随身wifi吗| 久久久国产精品麻豆| 在线观看三级黄色| 女人被躁到高潮嗷嗷叫费观| 成人毛片a级毛片在线播放| 搡老乐熟女国产| 国产精品蜜桃在线观看| 男女下面插进去视频免费观看 | 久久久国产精品麻豆| 在线观看三级黄色| 999精品在线视频| 色网站视频免费| 国产亚洲欧美精品永久| 欧美亚洲日本最大视频资源| 精品国产露脸久久av麻豆| 欧美国产精品一级二级三级| 中文字幕免费在线视频6| 中文字幕另类日韩欧美亚洲嫩草| 午夜免费鲁丝| 美女中出高潮动态图| 亚洲美女黄色视频免费看| 黄色毛片三级朝国网站| 国产色婷婷99| 26uuu在线亚洲综合色| 亚洲四区av| 亚洲性久久影院| 色婷婷av一区二区三区视频| 午夜激情av网站| 国产欧美亚洲国产| 国产激情久久老熟女| 欧美亚洲日本最大视频资源| 美女大奶头黄色视频| 午夜免费鲁丝| 久热这里只有精品99| 男人操女人黄网站| 黄色视频在线播放观看不卡| 亚洲成人av在线免费| 91久久精品国产一区二区三区| 中文字幕亚洲精品专区| 亚洲欧美精品自产自拍| 日日摸夜夜添夜夜爱| 国产一区有黄有色的免费视频| 最新的欧美精品一区二区| 男男h啪啪无遮挡| 国产高清三级在线| 国产 一区精品| av电影中文网址| 蜜桃在线观看..| 人人妻人人澡人人看| 国产精品99久久99久久久不卡 | 亚洲欧美成人精品一区二区| 在线观看三级黄色| 久久精品国产a三级三级三级| 中文天堂在线官网| 国产在视频线精品| 亚洲精品美女久久久久99蜜臀 | 精品少妇黑人巨大在线播放| 一级毛片电影观看| 在线观看一区二区三区激情| 草草在线视频免费看| 欧美xxxx性猛交bbbb| 精品酒店卫生间| 国产亚洲午夜精品一区二区久久| 久久人人97超碰香蕉20202| 午夜福利网站1000一区二区三区| 99九九在线精品视频| 久久久欧美国产精品| 国产成人91sexporn| 大香蕉97超碰在线| 精品视频人人做人人爽| 波多野结衣一区麻豆| 国产精品熟女久久久久浪| 捣出白浆h1v1| 成人国语在线视频| 丝袜脚勾引网站| 国产探花极品一区二区| 中文字幕人妻熟女乱码| xxxhd国产人妻xxx| 久久毛片免费看一区二区三区| 极品人妻少妇av视频| 亚洲人与动物交配视频| 国产av精品麻豆| 国产精品一区二区在线不卡| 日本欧美视频一区| 哪个播放器可以免费观看大片| 男人爽女人下面视频在线观看| 久久精品国产亚洲av涩爱| 高清av免费在线| 午夜精品国产一区二区电影| 黄片无遮挡物在线观看| 中文乱码字字幕精品一区二区三区| 看免费成人av毛片| 美女xxoo啪啪120秒动态图| 天天影视国产精品| 草草在线视频免费看| 韩国av在线不卡| 国产精品久久久久久久电影| 99久久综合免费| 午夜福利,免费看| xxx大片免费视频| 欧美3d第一页| 午夜日本视频在线| 丝袜脚勾引网站| 欧美最新免费一区二区三区| 日日爽夜夜爽网站| 91精品三级在线观看| 亚洲内射少妇av| 99久久人妻综合| av播播在线观看一区| 欧美 亚洲 国产 日韩一| 久久久久国产网址| 久久综合国产亚洲精品| kizo精华| www.熟女人妻精品国产 | 91久久精品国产一区二区三区| 精品亚洲成国产av| 99香蕉大伊视频| 狂野欧美激情性bbbbbb| 久久国产精品大桥未久av| 成人漫画全彩无遮挡| 国产色爽女视频免费观看| 日本av免费视频播放| 国产 精品1| 人体艺术视频欧美日本| 国产1区2区3区精品| 欧美老熟妇乱子伦牲交| 丰满乱子伦码专区| 成年女人在线观看亚洲视频| 国产成人91sexporn| 欧美亚洲日本最大视频资源| 日韩av不卡免费在线播放| 婷婷色综合www| 在线天堂中文资源库| www.熟女人妻精品国产 | 久久热在线av| 亚洲欧美成人精品一区二区| 亚洲av在线观看美女高潮| 制服诱惑二区| 亚洲av综合色区一区| 母亲3免费完整高清在线观看 | 亚洲精品第二区| 永久网站在线| 春色校园在线视频观看| 十八禁网站网址无遮挡| 99热国产这里只有精品6| 男人爽女人下面视频在线观看| 国产精品不卡视频一区二区| 国产一区二区在线观看日韩| 亚洲av欧美aⅴ国产| 黄片无遮挡物在线观看| 啦啦啦视频在线资源免费观看| 女的被弄到高潮叫床怎么办| 人人妻人人澡人人爽人人夜夜| 亚洲内射少妇av| 欧美 日韩 精品 国产| 亚洲久久久国产精品| 自拍欧美九色日韩亚洲蝌蚪91| 欧美成人精品欧美一级黄| 中文字幕人妻丝袜制服| 久久久精品94久久精品| 久久久久精品久久久久真实原创| 欧美最新免费一区二区三区| av不卡在线播放| 巨乳人妻的诱惑在线观看| av国产精品久久久久影院| 欧美亚洲 丝袜 人妻 在线| 免费看光身美女| 婷婷色av中文字幕| 咕卡用的链子| 国内精品宾馆在线| 搡老乐熟女国产| 少妇被粗大猛烈的视频| 天天操日日干夜夜撸| 久久亚洲国产成人精品v| 男女午夜视频在线观看 | 性色avwww在线观看| 视频中文字幕在线观看| 一本—道久久a久久精品蜜桃钙片| 伊人亚洲综合成人网| 国产毛片在线视频| 亚洲少妇的诱惑av| 少妇 在线观看| 亚洲国产最新在线播放| 免费观看性生交大片5| 女性被躁到高潮视频| av在线老鸭窝| 人人妻人人澡人人爽人人夜夜| 免费大片黄手机在线观看| 亚洲av电影在线观看一区二区三区| 777米奇影视久久| 久久免费观看电影| 男男h啪啪无遮挡| 成人国产麻豆网| 母亲3免费完整高清在线观看 | 岛国毛片在线播放| 日韩电影二区| 国产一区有黄有色的免费视频| 成人二区视频| 热re99久久精品国产66热6| 午夜av观看不卡| 制服丝袜香蕉在线| 岛国毛片在线播放| 久热这里只有精品99| 亚洲国产色片| 亚洲精品色激情综合| 日本午夜av视频| 免费人成在线观看视频色| 天天躁夜夜躁狠狠躁躁| 欧美日韩av久久| a 毛片基地| 国产精品久久久久久精品电影小说| 老司机影院成人| 久久人人爽人人片av| a 毛片基地| 日韩一区二区三区影片| 麻豆乱淫一区二区| 亚洲成人一二三区av| 欧美成人午夜精品| 青春草国产在线视频| 日韩一区二区三区影片| 另类精品久久| 亚洲少妇的诱惑av| 久久女婷五月综合色啪小说| 香蕉国产在线看| 欧美国产精品一级二级三级| 中文字幕精品免费在线观看视频 | 精品国产一区二区三区久久久樱花| 69精品国产乱码久久久| 久久午夜综合久久蜜桃| 老司机影院成人| 亚洲国产毛片av蜜桃av| 久久精品人人爽人人爽视色| 亚洲色图 男人天堂 中文字幕 | 啦啦啦中文免费视频观看日本| 亚洲精品成人av观看孕妇| 精品国产一区二区三区四区第35| 全区人妻精品视频| 日韩不卡一区二区三区视频在线| 国产亚洲最大av| 水蜜桃什么品种好| 亚洲国产看品久久| 久久女婷五月综合色啪小说| 欧美bdsm另类| 在线观看国产h片| 91精品国产国语对白视频| 最近的中文字幕免费完整| 亚洲欧美色中文字幕在线| 久久av网站| 国产精品 国内视频| 久久久久精品久久久久真实原创| 精品久久蜜臀av无| av在线播放精品| 黄色视频在线播放观看不卡| 老熟女久久久| 国产不卡av网站在线观看| 久久鲁丝午夜福利片| 国产女主播在线喷水免费视频网站| 成人黄色视频免费在线看| 在线观看三级黄色| 亚洲,一卡二卡三卡| 久久久久国产精品人妻一区二区| 亚洲美女搞黄在线观看| 免费黄频网站在线观看国产| 99re6热这里在线精品视频| 国产又色又爽无遮挡免| 在线观看人妻少妇| 人人妻人人澡人人爽人人夜夜| 免费在线观看黄色视频的| 久久青草综合色| √禁漫天堂资源中文www| 亚洲熟女精品中文字幕| 国产免费福利视频在线观看| 人人妻人人爽人人添夜夜欢视频| 午夜影院在线不卡| 国产精品一二三区在线看| 一本大道久久a久久精品| 伦理电影免费视频| av免费在线看不卡| 精品国产一区二区久久| 中文精品一卡2卡3卡4更新| 国产一区二区在线观看av| 一级毛片 在线播放| 午夜91福利影院| tube8黄色片| 国产av码专区亚洲av| 国产av码专区亚洲av| 国产毛片在线视频| 最黄视频免费看| 国产黄色免费在线视频| 下体分泌物呈黄色| 欧美激情国产日韩精品一区| 一边亲一边摸免费视频| 妹子高潮喷水视频| 国产精品 国内视频| 欧美 亚洲 国产 日韩一| 免费大片黄手机在线观看| 国产视频首页在线观看| 亚洲丝袜综合中文字幕| 亚洲精品乱久久久久久| 国产亚洲精品第一综合不卡 | 色哟哟·www| 亚洲精品久久成人aⅴ小说| 自拍欧美九色日韩亚洲蝌蚪91| 一边摸一边做爽爽视频免费| 亚洲中文av在线| 中文字幕av电影在线播放| 久久综合国产亚洲精品| av线在线观看网站| 日产精品乱码卡一卡2卡三| 9191精品国产免费久久| 色94色欧美一区二区| 最近最新中文字幕大全免费视频 | 777米奇影视久久| 在线亚洲精品国产二区图片欧美| 欧美精品一区二区大全| 美国免费a级毛片| 久久女婷五月综合色啪小说| 激情五月婷婷亚洲| 国产日韩欧美在线精品| 婷婷色综合大香蕉| 香蕉精品网在线| 男人操女人黄网站| 久久精品夜色国产| 国产成人精品久久久久久| 亚洲综合精品二区| 欧美日韩一区二区视频在线观看视频在线| 久久青草综合色| 熟女av电影| 一级毛片电影观看| 日日爽夜夜爽网站| 狂野欧美激情性bbbbbb| 亚洲成av片中文字幕在线观看 | 在线亚洲精品国产二区图片欧美| 满18在线观看网站| 午夜免费鲁丝| 亚洲第一av免费看| 香蕉国产在线看| 夜夜骑夜夜射夜夜干| 午夜免费男女啪啪视频观看| 亚洲国产av新网站| 侵犯人妻中文字幕一二三四区| 亚洲中文av在线| 欧美另类一区| 精品亚洲成国产av| 国产极品天堂在线| 男男h啪啪无遮挡| 国产精品国产av在线观看| 99香蕉大伊视频| 美女国产高潮福利片在线看| 大香蕉久久网| 男女边吃奶边做爰视频| 丝瓜视频免费看黄片| 1024视频免费在线观看| 夫妻午夜视频| 建设人人有责人人尽责人人享有的| 在线观看免费视频网站a站| 在线看a的网站| 99热国产这里只有精品6| 嫩草影院入口| 99re6热这里在线精品视频| 国产男女内射视频| 亚洲图色成人| 亚洲一区二区三区欧美精品| av国产久精品久网站免费入址| 伦精品一区二区三区| 大码成人一级视频| 九九爱精品视频在线观看| 亚洲色图 男人天堂 中文字幕 | 色婷婷久久久亚洲欧美| 免费女性裸体啪啪无遮挡网站| 亚洲国产最新在线播放| 捣出白浆h1v1| 两个人看的免费小视频| 亚洲丝袜综合中文字幕| av免费在线看不卡| 妹子高潮喷水视频| 亚洲激情五月婷婷啪啪| 欧美 亚洲 国产 日韩一| 国产毛片在线视频| 男女高潮啪啪啪动态图| 久久狼人影院| 18禁裸乳无遮挡动漫免费视频| 中文字幕亚洲精品专区| 久久久久久久久久久久大奶| 一区二区三区乱码不卡18| 岛国毛片在线播放| 青春草亚洲视频在线观看| 精品一品国产午夜福利视频| 亚洲国产精品专区欧美| 欧美 日韩 精品 国产| 国产精品偷伦视频观看了| 国产成人aa在线观看| 国产av精品麻豆| 国产综合精华液| 国产精品人妻久久久久久| 午夜久久久在线观看| 久久久a久久爽久久v久久| 精品国产乱码久久久久久小说| 妹子高潮喷水视频| 午夜日本视频在线| 水蜜桃什么品种好| 亚洲精品aⅴ在线观看| 人妻人人澡人人爽人人| 国产免费现黄频在线看| 伦理电影免费视频| 成人午夜精彩视频在线观看| 一级片'在线观看视频| 男女午夜视频在线观看 | xxxhd国产人妻xxx| 精品99又大又爽又粗少妇毛片| 国产高清不卡午夜福利| 黑人高潮一二区| 又黄又粗又硬又大视频| 亚洲精品一二三| 亚洲美女视频黄频| 亚洲欧美日韩另类电影网站| 看十八女毛片水多多多| 国产精品三级大全| 亚洲少妇的诱惑av| 菩萨蛮人人尽说江南好唐韦庄| 亚洲精品中文字幕在线视频| 亚洲av国产av综合av卡| 国产乱来视频区| 最新的欧美精品一区二区| 水蜜桃什么品种好| 男女啪啪激烈高潮av片| 国产男人的电影天堂91| 一边摸一边做爽爽视频免费| 亚洲精品视频女| 久久久久国产精品人妻一区二区| 1024视频免费在线观看| 亚洲精品第二区| 亚洲精品久久久久久婷婷小说| 少妇 在线观看| 亚洲av中文av极速乱| 人妻少妇偷人精品九色| 交换朋友夫妻互换小说| 丝袜脚勾引网站| 亚洲精品视频女| 一边摸一边做爽爽视频免费| 这个男人来自地球电影免费观看 | 制服丝袜香蕉在线| 国产亚洲精品久久久com| 男女免费视频国产| 丝袜脚勾引网站| 日韩av在线免费看完整版不卡| 亚洲av免费高清在线观看| 欧美日韩一区二区视频在线观看视频在线| 人妻人人澡人人爽人人| 亚洲第一区二区三区不卡| 老司机影院成人| 美女大奶头黄色视频| 97人妻天天添夜夜摸| 男女边摸边吃奶| 在线免费观看不下载黄p国产| freevideosex欧美| 久久人人爽人人片av| 久久99蜜桃精品久久|