• <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)系
    亚洲av电影在线观看一区二区三区| 丰满乱子伦码专区| 亚洲美女视频黄频| 在线观看三级黄色| 亚洲色图综合在线观看| 一级毛片我不卡| 成年人午夜在线观看视频| 午夜激情av网站| 日韩一区二区三区影片| 蜜桃久久精品国产亚洲av| 青春草视频在线免费观看| 人妻 亚洲 视频| 十八禁网站网址无遮挡| 中文字幕人妻丝袜制服| 高清毛片免费看| 国产精品国产三级国产专区5o| 亚洲精品乱码久久久v下载方式| 精品一区二区三卡| 国产爽快片一区二区三区| 亚洲美女搞黄在线观看| 日韩欧美一区视频在线观看| 黄色配什么色好看| 母亲3免费完整高清在线观看 | 少妇丰满av| 色婷婷av一区二区三区视频| 亚洲精品国产色婷婷电影| 人人妻人人添人人爽欧美一区卜| av国产久精品久网站免费入址| 大香蕉97超碰在线| 中文字幕制服av| 高清在线视频一区二区三区| 99re6热这里在线精品视频| 日本免费在线观看一区| 3wmmmm亚洲av在线观看| 男女国产视频网站| 人人妻人人添人人爽欧美一区卜| 五月天丁香电影| 丰满少妇做爰视频| 曰老女人黄片| 欧美日韩亚洲高清精品| 国产精品国产av在线观看| 伦理电影大哥的女人| 高清午夜精品一区二区三区| 日韩中字成人| 国产亚洲欧美精品永久| 久久久久国产网址| 免费黄色在线免费观看| 午夜激情福利司机影院| 欧美精品高潮呻吟av久久| 久久久欧美国产精品| 国产白丝娇喘喷水9色精品| 三级国产精品欧美在线观看| 一区二区av电影网| 免费大片18禁| 熟女人妻精品中文字幕| 在线观看人妻少妇| 久久鲁丝午夜福利片| 亚洲欧美一区二区三区黑人 | 99九九线精品视频在线观看视频| 性色avwww在线观看| 亚洲欧美日韩另类电影网站| 久久久久久久亚洲中文字幕| 免费高清在线观看日韩| 日韩三级伦理在线观看| 丝袜美足系列| 中文字幕制服av| av有码第一页| 国产精品秋霞免费鲁丝片| 91精品国产国语对白视频| 蜜桃国产av成人99| 久久久久久久久久久免费av| 最近最新中文字幕免费大全7| 在线观看三级黄色| 免费av中文字幕在线| 亚洲精品av麻豆狂野| 成年人午夜在线观看视频| 国产淫语在线视频| 国产精品国产三级专区第一集| 九色成人免费人妻av| 热re99久久国产66热| 午夜激情福利司机影院| 亚洲精品久久成人aⅴ小说 | 美女大奶头黄色视频| 久久ye,这里只有精品| 亚洲精品国产av成人精品| 丰满饥渴人妻一区二区三| 少妇精品久久久久久久| 国产一区二区三区综合在线观看 | 国产精品 国内视频| 欧美另类一区| 老司机影院毛片| 午夜视频国产福利| 三上悠亚av全集在线观看| 最近的中文字幕免费完整| 丝袜美足系列| 国产成人a∨麻豆精品| 9色porny在线观看| 亚洲综合色网址| 草草在线视频免费看| 最新的欧美精品一区二区| 久久免费观看电影| 精品熟女少妇av免费看| 黑丝袜美女国产一区| 看免费成人av毛片| 欧美3d第一页| 婷婷色综合www| 一区二区三区免费毛片| 免费观看性生交大片5| 精品少妇内射三级| 丝袜美足系列| 亚洲精品日本国产第一区| 亚洲人与动物交配视频| 国产视频内射| av女优亚洲男人天堂| 国语对白做爰xxxⅹ性视频网站| 成人漫画全彩无遮挡| 精品午夜福利在线看| 国产极品粉嫩免费观看在线 | 欧美97在线视频| 超碰97精品在线观看| 晚上一个人看的免费电影| 欧美变态另类bdsm刘玥| 人成视频在线观看免费观看| 国产精品国产三级国产av玫瑰| 国产黄片视频在线免费观看| 黄色视频在线播放观看不卡| 91成人精品电影| 日本黄色日本黄色录像| 国产精品国产三级专区第一集| 久久久国产一区二区| 亚洲国产最新在线播放| 自线自在国产av| a级片在线免费高清观看视频| 黑丝袜美女国产一区| 2022亚洲国产成人精品| 精品一区在线观看国产| av女优亚洲男人天堂| 老女人水多毛片| 久久精品夜色国产| 久久精品国产鲁丝片午夜精品| 看免费成人av毛片| 少妇丰满av| 欧美日韩视频高清一区二区三区二| 国产精品久久久久久精品电影小说| 纯流量卡能插随身wifi吗| 国产一区有黄有色的免费视频| 看十八女毛片水多多多| 国产国语露脸激情在线看| 国产成人av激情在线播放 | 亚洲欧美色中文字幕在线| 少妇高潮的动态图| 最黄视频免费看| 青青草视频在线视频观看| 97精品久久久久久久久久精品| 国产深夜福利视频在线观看| 18+在线观看网站| 日本免费在线观看一区| 国产成人91sexporn| 久久久久久久精品精品| 丝袜在线中文字幕| 亚洲人成网站在线观看播放| 久久精品国产亚洲av涩爱| 久久99蜜桃精品久久| 欧美亚洲 丝袜 人妻 在线| 中文欧美无线码| 久久久精品94久久精品| 简卡轻食公司| 成人午夜精彩视频在线观看| 久久99一区二区三区| 亚洲av.av天堂| 国产免费视频播放在线视频| 免费黄网站久久成人精品| 一个人看视频在线观看www免费| 亚洲精品美女久久av网站| 母亲3免费完整高清在线观看 | 日本爱情动作片www.在线观看| 亚洲国产精品成人久久小说| 久久久精品免费免费高清| 永久网站在线| 大片电影免费在线观看免费| 色婷婷久久久亚洲欧美| 免费人成在线观看视频色| 久久久久国产网址| 亚洲久久久国产精品| 考比视频在线观看| 国产亚洲一区二区精品| 简卡轻食公司| 久久久精品区二区三区| 美女视频免费永久观看网站| 国产精品一二三区在线看| 人妻少妇偷人精品九色| 少妇人妻精品综合一区二区| 秋霞在线观看毛片| 99热全是精品| 蜜桃国产av成人99| 九九爱精品视频在线观看| 黑人巨大精品欧美一区二区蜜桃 | 国产精品嫩草影院av在线观看| 一个人免费看片子| 蜜桃在线观看..| 夜夜看夜夜爽夜夜摸| 菩萨蛮人人尽说江南好唐韦庄| 嫩草影院入口| 久久精品国产a三级三级三级| 伦理电影大哥的女人| 日韩制服骚丝袜av| 一区二区三区乱码不卡18| 欧美精品亚洲一区二区| 菩萨蛮人人尽说江南好唐韦庄| 亚洲精品日韩av片在线观看| 国产极品粉嫩免费观看在线 | 91精品国产国语对白视频| 日本av手机在线免费观看| 人人妻人人添人人爽欧美一区卜| 啦啦啦中文免费视频观看日本| 久久久久久久精品精品| 亚洲欧美一区二区三区国产| 热99国产精品久久久久久7| 一级毛片我不卡| 久久久精品免费免费高清| 日韩不卡一区二区三区视频在线| 麻豆乱淫一区二区| 在线看a的网站| 日韩成人av中文字幕在线观看| 国产乱人偷精品视频| 中文天堂在线官网| 丝袜美足系列| 久久鲁丝午夜福利片| 国产男人的电影天堂91| 国产欧美日韩一区二区三区在线 | 热99国产精品久久久久久7| 国产毛片在线视频| 色吧在线观看| 在线亚洲精品国产二区图片欧美 | 妹子高潮喷水视频| 国产有黄有色有爽视频| 亚洲精品久久午夜乱码| 赤兔流量卡办理| 九九久久精品国产亚洲av麻豆| 亚洲人与动物交配视频| 日本-黄色视频高清免费观看| 久久影院123| 亚洲av在线观看美女高潮| 精品人妻偷拍中文字幕| 亚洲国产最新在线播放| 五月开心婷婷网| √禁漫天堂资源中文www| av国产精品久久久久影院| 久久久久久人妻| 人成视频在线观看免费观看| 精品久久久精品久久久| 日韩熟女老妇一区二区性免费视频| 午夜激情久久久久久久| 老熟女久久久| 纵有疾风起免费观看全集完整版| 99热网站在线观看| 看免费成人av毛片| 如日韩欧美国产精品一区二区三区 | 欧美一级a爱片免费观看看| 欧美+日韩+精品| 一级毛片黄色毛片免费观看视频| 亚洲精品日韩在线中文字幕| 九九久久精品国产亚洲av麻豆| 最近中文字幕高清免费大全6| 国产日韩欧美视频二区| 免费少妇av软件| 在线天堂最新版资源| 日韩制服骚丝袜av| 黄色毛片三级朝国网站| 亚洲精品456在线播放app| 午夜影院在线不卡| 欧美日韩一区二区视频在线观看视频在线| 亚洲欧洲国产日韩| 久久韩国三级中文字幕| 少妇精品久久久久久久| 中文字幕人妻丝袜制服| 天堂8中文在线网| 日韩视频在线欧美| 日本vs欧美在线观看视频| a级毛色黄片| 夜夜看夜夜爽夜夜摸| 日日摸夜夜添夜夜爱| 日韩一区二区视频免费看| 五月伊人婷婷丁香| 国产日韩欧美视频二区| 视频在线观看一区二区三区| 亚洲色图综合在线观看| 十八禁高潮呻吟视频| 欧美激情国产日韩精品一区| 亚洲精品自拍成人| 一级毛片电影观看| 在线观看www视频免费| 欧美激情国产日韩精品一区| 一本—道久久a久久精品蜜桃钙片| 一级黄片播放器| 岛国毛片在线播放| 极品人妻少妇av视频| 夜夜骑夜夜射夜夜干| 精品视频人人做人人爽| 成人手机av| 王馨瑶露胸无遮挡在线观看| 国产日韩欧美亚洲二区| 久久99一区二区三区| 中文精品一卡2卡3卡4更新| 国产黄色视频一区二区在线观看| 欧美变态另类bdsm刘玥| 汤姆久久久久久久影院中文字幕| a级毛片免费高清观看在线播放| 最近的中文字幕免费完整| 熟女av电影| 免费大片黄手机在线观看| 日韩中字成人| 夜夜骑夜夜射夜夜干| 日韩精品免费视频一区二区三区 | 久久久久视频综合| 女性被躁到高潮视频| 男女无遮挡免费网站观看| 亚洲精品视频女| 免费大片18禁| 18禁观看日本| 永久网站在线| 男女国产视频网站| 国产69精品久久久久777片| 人妻一区二区av| 五月伊人婷婷丁香| 99久久精品国产国产毛片| 欧美变态另类bdsm刘玥| 国产极品天堂在线| 我要看黄色一级片免费的| 久久久欧美国产精品| 日韩精品有码人妻一区| 热99国产精品久久久久久7| 18+在线观看网站| 搡女人真爽免费视频火全软件| 亚洲四区av| 国产精品国产av在线观看| 日韩伦理黄色片| 中文乱码字字幕精品一区二区三区| 日本黄色片子视频| .国产精品久久| 又粗又硬又长又爽又黄的视频| 最近手机中文字幕大全| 久久婷婷青草| 国产午夜精品久久久久久一区二区三区| av在线播放精品| 又黄又爽又刺激的免费视频.| 边亲边吃奶的免费视频| 观看美女的网站| 成年av动漫网址| 国产免费福利视频在线观看| av天堂久久9| 亚洲精品一区蜜桃| 欧美日韩一区二区视频在线观看视频在线| 伦理电影免费视频| 51国产日韩欧美| 成年人免费黄色播放视频| 美女国产视频在线观看| 精品人妻熟女av久视频| 久久精品熟女亚洲av麻豆精品| 街头女战士在线观看网站| 热re99久久国产66热| 91精品一卡2卡3卡4卡| 亚洲怡红院男人天堂| 最后的刺客免费高清国语| 男人添女人高潮全过程视频| 国产黄频视频在线观看| 久久综合国产亚洲精品| 婷婷成人精品国产| 三级国产精品欧美在线观看| 久久ye,这里只有精品| 久热久热在线精品观看| 熟女电影av网| 狂野欧美激情性xxxx在线观看| 国产在线一区二区三区精| 亚洲欧美精品自产自拍| 观看美女的网站| av线在线观看网站| 日本-黄色视频高清免费观看| 亚洲欧洲精品一区二区精品久久久 | 女性生殖器流出的白浆| 男女国产视频网站| 久久精品国产亚洲av天美| 日本免费在线观看一区| 国产精品免费大片| 日日啪夜夜爽| 亚洲精品,欧美精品| av电影中文网址| kizo精华| 我的女老师完整版在线观看| 午夜激情福利司机影院| 亚洲av二区三区四区| 国产精品欧美亚洲77777| 在线观看www视频免费| 我的老师免费观看完整版| 黑人欧美特级aaaaaa片| 中文字幕人妻熟人妻熟丝袜美| 一级,二级,三级黄色视频| 国产av国产精品国产| 青春草视频在线免费观看| 成人黄色视频免费在线看| 一区二区三区乱码不卡18| 91精品一卡2卡3卡4卡| 久久精品久久久久久噜噜老黄| 午夜福利影视在线免费观看| 亚洲性久久影院| 女人精品久久久久毛片| 国产视频内射| 亚洲成色77777| 成人影院久久| 国产亚洲av片在线观看秒播厂| 少妇人妻 视频| 亚洲欧美日韩另类电影网站| 亚洲欧美一区二区三区黑人 | 精品人妻偷拍中文字幕| 毛片一级片免费看久久久久| 久热这里只有精品99| 美女主播在线视频| 国产精品嫩草影院av在线观看| 欧美人与善性xxx| 777米奇影视久久| 99热这里只有精品一区| 22中文网久久字幕| 欧美日韩国产mv在线观看视频| 久久久欧美国产精品| 国产一区二区在线观看av| 女的被弄到高潮叫床怎么办| 母亲3免费完整高清在线观看 | 欧美老熟妇乱子伦牲交| 日韩大片免费观看网站| 婷婷成人精品国产| 精品人妻偷拍中文字幕| av国产久精品久网站免费入址| 亚洲无线观看免费| 亚洲av福利一区| 日日啪夜夜爽| 欧美日韩精品成人综合77777| 亚洲五月色婷婷综合| 欧美最新免费一区二区三区| 久久午夜综合久久蜜桃| 成人综合一区亚洲| 一边亲一边摸免费视频| 国产av一区二区精品久久| 成人漫画全彩无遮挡| 国产成人freesex在线| 国产av国产精品国产| 国产av精品麻豆| 高清欧美精品videossex| 永久免费av网站大全| 免费黄色在线免费观看| 日本av免费视频播放| 三上悠亚av全集在线观看| 青春草国产在线视频| 我要看黄色一级片免费的| 伊人久久精品亚洲午夜| 激情五月婷婷亚洲| 亚洲国产精品国产精品| 在线观看免费日韩欧美大片 | 91精品一卡2卡3卡4卡| 又粗又硬又长又爽又黄的视频| 九九爱精品视频在线观看| 日韩一本色道免费dvd| 好男人视频免费观看在线| 伊人久久精品亚洲午夜| 欧美另类一区| 极品少妇高潮喷水抽搐| freevideosex欧美| 亚洲精品国产av成人精品| 久热这里只有精品99| 成人国产麻豆网| 国产精品嫩草影院av在线观看| 久久鲁丝午夜福利片| 亚洲内射少妇av| 97在线人人人人妻| 日本av免费视频播放| 狂野欧美激情性bbbbbb| av.在线天堂| 在线观看免费高清a一片| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 日韩亚洲欧美综合| 久久鲁丝午夜福利片| 国产片内射在线| 美女中出高潮动态图| 精品一区二区免费观看| 婷婷成人精品国产| 亚洲精品久久久久久婷婷小说| 精品人妻熟女av久视频| 午夜福利在线观看免费完整高清在| 久久 成人 亚洲| 国产成人一区二区在线| 一级毛片电影观看| 在线观看免费视频网站a站| 五月玫瑰六月丁香| 99热这里只有是精品在线观看| 男的添女的下面高潮视频| 人妻少妇偷人精品九色| 免费观看av网站的网址| 亚洲高清免费不卡视频| 日本黄色日本黄色录像| 久久久久人妻精品一区果冻| 少妇被粗大的猛进出69影院 | 69精品国产乱码久久久| 久久精品久久久久久久性| 亚洲av中文av极速乱| 简卡轻食公司| 亚洲激情五月婷婷啪啪| 人妻少妇偷人精品九色| 亚洲成色77777| 免费观看无遮挡的男女| 一个人看视频在线观看www免费| av天堂久久9| 国产成人一区二区在线| 搡老乐熟女国产| av专区在线播放| 国产精品三级大全| 两个人免费观看高清视频| 自拍欧美九色日韩亚洲蝌蚪91| 中文字幕免费在线视频6| 五月伊人婷婷丁香| 美女视频免费永久观看网站| 亚洲精品日韩av片在线观看| 超碰97精品在线观看| 一区二区三区精品91| 久久久久国产精品人妻一区二区| 交换朋友夫妻互换小说| 妹子高潮喷水视频| 日产精品乱码卡一卡2卡三| 99热6这里只有精品| 亚洲激情五月婷婷啪啪| 亚洲欧美一区二区三区国产| 久久免费观看电影| 国产老妇伦熟女老妇高清| 大陆偷拍与自拍| 亚洲国产最新在线播放| 久久99一区二区三区| 欧美3d第一页| 97超视频在线观看视频| 久久婷婷青草| 在线观看免费高清a一片| 国产成人精品无人区| 国产一级毛片在线| 黑人高潮一二区| 国产免费又黄又爽又色| 国产亚洲最大av| 看十八女毛片水多多多| 熟女人妻精品中文字幕| 亚洲精品日韩av片在线观看| 亚洲国产精品国产精品| 成人影院久久| 国产精品久久久久久久电影| 久久99热6这里只有精品| 久久久久久久久大av| 日日摸夜夜添夜夜添av毛片| 久久午夜福利片| 欧美精品人与动牲交sv欧美| 狠狠精品人妻久久久久久综合| 777米奇影视久久| 久久久久久久久久成人| 最后的刺客免费高清国语| 18在线观看网站| 老女人水多毛片| 一级二级三级毛片免费看| 国产av码专区亚洲av| 另类精品久久| 少妇丰满av| 精品人妻在线不人妻| 91精品伊人久久大香线蕉| 视频在线观看一区二区三区| 伊人久久国产一区二区| 9色porny在线观看| 日本-黄色视频高清免费观看| 免费黄色在线免费观看| 十分钟在线观看高清视频www| 欧美激情 高清一区二区三区| 一区二区三区乱码不卡18| 国产精品无大码| 老女人水多毛片| 国产综合精华液| 成人毛片60女人毛片免费| 亚洲国产成人一精品久久久| 丰满乱子伦码专区| 日韩成人伦理影院| 国产av国产精品国产| 91精品国产九色| 全区人妻精品视频| 一本久久精品| 制服丝袜香蕉在线| 午夜激情av网站| 国产精品国产三级专区第一集| 亚洲av综合色区一区| 日日撸夜夜添| 精品视频人人做人人爽| 国产精品免费大片| 青春草视频在线免费观看| 中文字幕最新亚洲高清| 性高湖久久久久久久久免费观看| 精品国产国语对白av| 夜夜看夜夜爽夜夜摸| 欧美精品亚洲一区二区| 18在线观看网站| 国产欧美日韩一区二区三区在线 | 亚洲人与动物交配视频| 只有这里有精品99| 国产一级毛片在线| 少妇被粗大猛烈的视频| 亚洲精品乱码久久久v下载方式| 色吧在线观看| 一区二区日韩欧美中文字幕 | 亚洲国产精品专区欧美| 丝袜喷水一区| 国产亚洲欧美精品永久| av视频免费观看在线观看| 免费大片18禁| 国产白丝娇喘喷水9色精品| 中文字幕精品免费在线观看视频 | 国产极品天堂在线|