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

    基于“三年化疫”理論探討百日咳發(fā)病與前期氣象因素的相關性并建立預測模型*

    2015-06-01 12:33:28軒,賀
    西部中醫(yī)藥 2015年11期
    關鍵詞:百日咳北京地區(qū)時段

    張 軒,賀 娟

    北京中醫(yī)藥大學基礎醫(yī)學院,北京 100029

    基于“三年化疫”理論探討百日咳發(fā)病與前期氣象因素的相關性并建立預測模型*

    張 軒,賀 娟△

    北京中醫(yī)藥大學基礎醫(yī)學院,北京 100029

    目的:探討北京地區(qū)百日咳發(fā)病與前期(1~3年前)氣象因素的關聯(lián)性,并建立BP人工神經(jīng)網(wǎng)絡醫(yī)療氣象預測模型。方法:收集北京地區(qū)1970—2004年35年的氣象資料和百日咳發(fā)病資料,采用BP人工神經(jīng)網(wǎng)絡方法,分別從1年前、2年前、3年前3個不同時間維度,建立百日咳發(fā)病的氣象預警模型。結果:1)百日咳發(fā)病與前期(1~3年前)氣象因素具有相關性,其中關系最密切的氣象因素是1~3年前的平均相對濕度;2)利用前期氣象因素皆可成功建立百日咳發(fā)病的預測模型,以1年前氣象因素建模的預測效果最佳。結論:1~3年前的氣候可能會影響某些傳染病的發(fā)病,今后對傳染病發(fā)病的研究應注意考慮前期的氣象變化。

    百日咳;BP人工神經(jīng)網(wǎng)絡;前期氣象變化;三年化疫;五運六氣

    百日咳是由百日咳桿菌引起的急性呼吸道傳染病,臨床表現(xiàn)以陣發(fā)性痙攣性咳嗽,以及咳嗽終止時伴有雞鳴樣吸氣吼聲為特征,多發(fā)生于兒童[1]。近年研究[2-4]開始重視從氣象學角度探討引起百日咳發(fā)病的氣候條件,旨在為其防控提供氣象決策依據(jù)。但目前大多數(shù)研究仍側重于分析同期氣象變化,而鮮有關于前期氣候對百日咳發(fā)病影響的報道。

    五運六氣學說是中醫(yī)理論的重要組成部分,文字數(shù)量占《黃帝內(nèi)經(jīng)·素問》三分之一還要多,有舉足輕重的地位[5]。該學說中的“三年化疫”理論認為疫病的出現(xiàn)不但與當時的氣候有關,而且與近3年的氣候變化有關系[6]。這一觀點延伸了對前期氣象研究的時間范圍,提示可以選擇1~3年前的氣象資料分析其與傳染病發(fā)病的關聯(lián)性。本課題既往的研究結果發(fā)現(xiàn)某些傳染病的發(fā)生與前期的氣象因素有關:如痢疾[7]和傷寒副傷寒[8]均與3年前的氣候變化存在確切相關性;通過前1年的氣象變化可以預測當年感染性腹瀉高發(fā)月份的病例數(shù)[9];利用2年前的氣象因素可建立流行性腮腺炎發(fā)病的回歸預測模型,且擬合效果最佳[10]。

    因此,本研究試圖利用北京地區(qū)1970—2004年35年的傳染病數(shù)據(jù)庫和氣象數(shù)據(jù)庫,采用統(tǒng)計學方法,探討百日咳發(fā)病與1~3年前氣象變化的相關性,并建立氣象預警模型。

    1 資料與方法

    1.1 研究地點 北京地區(qū)屬典型的暖溫帶半濕潤大陸性季風氣候,四季分明。春季干旱,夏季炎熱多雨,秋季天高氣爽,冬季寒冷干燥。風向有明顯的季節(jié)變化。有學者[11]研究發(fā)現(xiàn)北京全年六氣階段的實際氣候特點與中醫(yī)六氣模式“厥陰風木、少陰君火、少陽相火、太陰濕土、陽明燥金、太陽寒水特點”基本符合,認為北京地區(qū)的氣候特點比較適合運氣理論的研究。因此,將研究地點選定在北京地區(qū),并采用六氣分段的方式處理數(shù)據(jù),以期較為客觀地分析氣象因素與百日咳發(fā)病之間的關系。

    1.2 資料來源 北京市疾病預防控制中心提供1970—2004年記錄存檔的百日咳逐月發(fā)病數(shù)據(jù),以及北京市觀象臺提供的同期7項基本氣象因子的數(shù)據(jù)資料,包括日平均氣溫、日平均風速、日降水量、日平均相對濕度、日平均水汽壓、日低云量、日舒適度。

    1.3 研究方法

    1.3.1 數(shù)據(jù)處理 本研究選擇以立春為起始點[12],將每年按24節(jié)氣均分為六個時段,即初之氣(立春~春分)、二之氣(清明~小滿)、三之氣(芒種~大暑)、四之氣(立秋~秋分)、五之氣(寒露~小雪)、終之氣(大雪~大寒)。以六氣時段為單位,統(tǒng)計35年間各時段內(nèi)氣溫、風速、降水量、相對濕度、水汽壓、低云量、舒適度7項氣象因子的平均值;以及1970—2003年百日咳逐年的六段發(fā)病數(shù)(因缺少2005年1月的發(fā)病數(shù)據(jù),故無法統(tǒng)計2004年終之氣的發(fā)病數(shù))。

    1.3.2 統(tǒng)計學方法 數(shù)據(jù)采用SPSS 17.0統(tǒng)計軟件處理。由于人工神經(jīng)網(wǎng)絡具有獨特的信息存儲方式、良好的容錯性、大規(guī)模的非線性并行處理方式以及強大的自組織、自學習和自適應能力[13],故近年來該方法被廣泛應用于傳染病預測及與氣象因素關系的分析中[14-17]。針對數(shù)據(jù)的非線性特點,本研究將采用BP人工神經(jīng)網(wǎng)絡建立百日咳發(fā)病的氣象預測模型。模型1(1年前):選擇1973—2003年百日咳六氣時段的發(fā)病數(shù)及1972—2002年7個氣象因子的六段均值數(shù)進行統(tǒng)計,各自包含186組數(shù)據(jù)。將各氣象因子作為輸入變量,百日咳發(fā)病數(shù)作為輸出變量,采用統(tǒng)計學中常用且默認的7∶3比例隨機抽取數(shù)據(jù)進行訓練及預測,建立預測模型,評價其預測精度,并分析1年前的各氣象因子對百日咳發(fā)病的影響程度。模型2(2年前):選擇1973—2003年百日咳六氣時段的發(fā)病數(shù)及1971—2001年7個氣象因子的六段均值數(shù)進行統(tǒng)計,各自包含186組數(shù)據(jù)。如上法建立預測模型,評價其預測精度,并分析2年前的各氣象因子對百日咳發(fā)病的影響程度。模型3(3年前):選擇1973—2003年百日咳六氣時段的發(fā)病數(shù)及7個氣象因子的六段均值數(shù)進行統(tǒng)計,各自包含186組數(shù)據(jù)。如上法建立預測模型,評價其預測精度,并分析3年前的各氣象因子對百日咳發(fā)病的影響程度。

    2 結果

    2.1 各氣象因素的六氣時段變化趨勢 在六氣時段中,有6個氣象因子(即:氣溫、降水量、相對濕度、水汽壓、低云量和舒適度)的波峰值均出現(xiàn)在三之氣(芒種~大暑),其中相對濕度在三之氣和四之氣(立秋~秋分)均呈現(xiàn)波峰情況;6個氣象因子的波谷值均出現(xiàn)在終之氣(大雪~大寒),其次為初之氣(立春~春分)。只有風速在三之氣的變化并不明顯,其最高值出現(xiàn)在二之氣(清明~小滿),其次為初之氣,最低值出現(xiàn)在四之氣,見圖1—7。

    圖1 北京地區(qū)1970—2004年氣溫的六氣時段變化趨勢

    圖2 北京地區(qū)1970—2004年風速的六氣時段變化趨勢

    圖3 北京地區(qū)1970—2004年降水量的六氣時段變化趨勢

    圖4 北京地區(qū)1970—2004年相對濕度的六氣時段變化趨勢

    圖5 北京地區(qū)1970—2004年水汽壓的六氣時段變化趨勢

    圖6 北京地區(qū)1970—2004年低云量的六氣時段變化趨勢

    圖7 北京地區(qū)1970—2004年舒適度的六氣時段變化趨勢

    2.2 百日咳六氣時段的發(fā)病特點 在六氣時段中,百日咳在三之氣呈現(xiàn)高發(fā)趨勢,在終之氣呈現(xiàn)低發(fā)趨勢。三之氣的發(fā)病與其他各時段的發(fā)病情況比較差異有統(tǒng)計學意義(χ2=11.920,P<0.05),即確定三之氣為百日咳的高發(fā)病時段,見圖8。

    圖8 北京地區(qū)1970—2004年百日咳六氣時段發(fā)病趨勢

    2.3 百日咳發(fā)病與1年前氣象因素的相關性及氣象預測模型 結果顯示,該模型訓練樣本占數(shù)據(jù)的73.7%,測試樣本占數(shù)據(jù)的26.3%。訓練相對錯誤值為0.330,測試的相對錯誤值為0.365,均小于1,說明模型建立成功,預測精度約為66%。其中,與百日咳發(fā)病密切相關的氣象因素(即標準化重要性超過50%)有:1年前的相對濕度、1年前的風速、1年前的低云量、1年前的降水量。其中,相對濕度的貢獻度最大,達100%,見表1。

    表1 1年前各氣象因素的重要性

    2.4 百日咳發(fā)病與2年前氣象因素的相關性及氣象預測模型 結果顯示,該模型訓練樣本占數(shù)據(jù)的72.0%,測試樣本占數(shù)據(jù)的28.0%。訓練相對錯誤值為0.524,測試的相對錯誤值為0.444,均小于1,說明模型建立成功,預測精度約為52%。其中,與百日咳發(fā)病密切相關的氣象因素(即標準化重要性超過50%)有:2年前的相對濕度和2年前的風速。其中,相對濕度的貢獻度最大,達100%,見表2。

    2.5 百日咳發(fā)病與3年前氣象因素的相關性及氣象預測模型 結果顯示,該模型訓練樣本占數(shù)據(jù)的72.0%,測試樣本占數(shù)據(jù)的28.0%。訓練相對錯誤值為0.787,測試的相對錯誤值為0.576,均小于1,說明模型建立成功,預測精度約為32%。其中,與百日咳發(fā)病密切相關的氣象因素(即標準化重要性超過50%)有:3年前的相對濕度、3年前的風速、3年前的低云量、3年前的降水量。其中,相對濕度的貢獻度最大,達100%,見表3。

    表2 2年前各氣象因素的重要性

    表3 3年前各氣象因素的重要性

    3 討論

    本研究顯示,北京地區(qū)百日咳高發(fā)病的三之氣時段,正是多種氣象因素(包括氣溫、降水量、相對濕度、水汽壓、低云量、舒適度等)高峰出現(xiàn)的時期;而BP人工神經(jīng)網(wǎng)絡統(tǒng)計結果顯示,北京地區(qū)百日咳發(fā)病與1~3年前的相對濕度、風速、降水量、低云量等氣象因素具有密切相關性,其中相對濕度是最重要的氣象因子。本研究結果與課題之前對百日咳發(fā)病與同期氣象因素相關性研究[18]的結果對比,發(fā)現(xiàn)無論在當年,還是1~3年前,相對濕度都是誘發(fā)百日咳發(fā)病的最關鍵因素。在利用BP人工神經(jīng)網(wǎng)絡建立醫(yī)療氣象預測模型方面,本研究顯示,利用3年前的氣象因素建模,其訓練相對錯誤值為0.787,測試的相對錯誤值為0.576,推斷模型預測精度約為32%,效果并不理想。利用2年前的氣象因素建模,其訓練相對錯誤值為0.524,測試的相對錯誤值為0.444,推斷預測精度約為52%,效果尚可。利用1年前的氣象因素建模,其訓練相對錯誤值為0.330,測試的相對錯誤值為0.365,推斷預測精度約為66%,效果優(yōu)于前2個模型。既往的研究結果顯示,利用當年的氣象因素建模,其訓練相對錯誤值為0.855,測試的相對錯誤值為0.805,推斷預測精度約為18%,效果較差。兩者對比發(fā)現(xiàn),利用前期氣象因素建立百日咳發(fā)病的預測模型,效果優(yōu)于利用同期氣象因素建立的模型。其中,利用1年前氣象因素建模的預測效果最佳。該結果提示,北京地區(qū)百日咳的發(fā)病可能與上一年的氣候變化密切相關。

    本研究的創(chuàng)新之處在于:1)突破同期氣象與傳染病發(fā)病的研究范疇[19-21],基于中醫(yī)運氣“三年化疫”的特色理論,從前期氣象變化進行分析,并將前期時間范圍定義在1~3年前,利用前期氣象因素建立百日咳發(fā)病的預測模型,從時間上促進預警系統(tǒng)關口的前移,凸顯了中醫(yī)運氣理論預測的超前性特征。2)針對傳染病數(shù)據(jù)非正態(tài)的分布特點,以及氣象因子之間的交叉作用,選擇較為合適的統(tǒng)計方法,即BP人工神經(jīng)網(wǎng)絡建立預測模型,并可同時分析關系密切的氣象因素,與其他研究中常用的單因素相關分析[22]、多元回歸[23]、指數(shù)平滑[24]、時間序列[25]等方法不同。

    但是,本研究尚具有一定的局限性:由于資料有限,研究僅選擇了7個基本氣象因子作為研究對象,而氣候的變化是多個氣象因子綜合作用的結果,各氣象因子之間存在著不同的影響權重。在以后研究中,筆者認為應擴大入選的氣象指標,并引入氣象綜合參數(shù)等概念,參考能反映氣象綜合變動的數(shù)學模型,用于分析不同氣候特點對百日咳發(fā)病的影響,以期建立精度更高的預測模型。

    [1] 楊紹基,任紅.傳染病學[M].北京:人民衛(wèi)生出版社,2011: 194.

    [2] 周國琪,楊柳,鄒純樸.運氣學司天之氣與溫州地區(qū)氣象及百日咳發(fā)病規(guī)律相關性的研究[J].中國中醫(yī)基礎醫(yī)學雜志,2008,14(11):805-806.

    [3] 鄒純樸,韓淑杰,王建,等.東北三地猩紅熱、百日咳發(fā)病率與氣象因素相關性分析[J].2008,10(11):40-41.

    [4] 王清,果麗平,鄧健,等.邯鄲市百日咳發(fā)病率與氣象因素之間的關系研究[J].醫(yī)學動物防制,2014,30(3):246-248.

    [5] 蘇穎.五運六氣探微[M].北京:人民衛(wèi)生出版社,2014:3.

    [6] 蘇穎.中醫(yī)運氣學[M].北京:中國中醫(yī)藥出版社,2009:93.

    [7] 劉忠第,張軒,賀娟.基于“3年化疫”理論探討氣象因素與北京地區(qū)痢疾發(fā)病的關系[J].北京中醫(yī)藥大學學報,2013,36(4):258-262.

    [8] 張軒,祖麗胡瑪爾·艾尼瓦爾,劉忠第,等.基于運氣理論探討傷寒副傷寒與氣象因素的相關性并建立預測模型[J].北京中醫(yī)藥大學學報,2013,36(12):797-801.

    [9] 湯巧玲,付幫澤,黃羚,等.北京市感染性腹瀉與氣象及六氣的關聯(lián)性研究[J].中醫(yī)雜志,2013,54(7):570-572.

    [10]張軒,賀娟.北京地區(qū)1990—2004年流行性腮腺炎發(fā)病與氣象因素相關性研究[J].中華中醫(yī)藥雜志,2014,29(5):1740-1745.

    [11]張德山,高思華,李鴻濤,等.北京地區(qū)中醫(yī)六氣的氣候分析[J].中醫(yī)研究,2010,23(9):1-3.

    [12]劉玉慶,高思華,等.從北京市60年氣象資料看王冰注釋運氣計算模式的謬誤[J].北京中醫(yī)藥大學學報,2010,12(33): 813-816.

    [13]孫振球.醫(yī)學統(tǒng)計學[M].北京:人民衛(wèi)生出版社,2007:478.

    [14]馬師雷,賀娟,高思華.從六氣角度探討病毒性肝炎與氣象因子的相關性及疾病預測模型[J].中華中醫(yī)藥雜志,2012,27(4):823-827.

    [15]胡婧媛,蔣夢姣,景元書,等.基于神經(jīng)網(wǎng)絡的氣象條件對瀘州市肺結核發(fā)病率預測[J].科技通報,2013,29(5):19-23.

    [16]鄧斌,周志剛,馬澤粦,等.肺結核病與氣象因素關系的BP神經(jīng)網(wǎng)絡模型研究[J].國際醫(yī)藥衛(wèi)生導報,2008,14(1):17-20.

    [17]曲波,黃德生,郭海強,等.流行性腦脊髓膜炎與氣象因素關系的BP神經(jīng)網(wǎng)絡模型研究[J].中國醫(yī)科大學學報,2006,35(2):158-159.

    [18]湯巧玲,付幫澤,劉宏偉,等.以百日咳為例探討六氣、氣象與傳染性疾病的關聯(lián)性研究方法[J].中華中醫(yī)藥雜志,2013,28(5):1301-1305.

    [19]黎新宇,王全意,賈蕾,等.北京市傷寒發(fā)病率與氣象因素關系研究[J].首都公共衛(wèi)生,2007,1(2):91-93.

    [20]M a SL,Tang Q L,Li u H W,et al,等.Correl at i on anal ysi s for t he at t ack of baci l l ary dysent ery and m et eorol ogi cal fact ors based on t he Chi nese m edi ci ne t heory ofYunQi and t he m edi cal-m et eorol ogi calforecast m odel[J].Chi nese JournalofInt egrat i ve M edi ci ne,2013,19(3):182-186.

    [21]賈蕾,黎新宇,劉桂榮,等.北京氣象因素與細菌性痢疾發(fā)病關系的研究[J].現(xiàn)代預防醫(yī)學,2007,34(13):2470-2471.

    [22]鄭能雄,林云欽,鄭高,等.福州市主要傳染病與氣象因素關系的研究[J].醫(yī)學動物防制,2004,20(5):275-278.

    [23]安慶玉,吳雋,王曉立,等.氣象因素變化與大連市腸道傳染病發(fā)病時間分布關系的研究[J].中國預防醫(yī)學雜志,2012,13(4):288-291.

    [24]史蕓萍,馬家奇.指數(shù)平滑法在流行性腮腺炎預測預警中的應用[J].中國疫苗和免疫,2010,16(3):233-237.

    [25]王春平,王志鋒,單杰,等.隨機時間序列分析法在傳染病預測中的應用[J].中國醫(yī)藥統(tǒng)計,2006,13(3):229-232.

    Exploration on the Correlation between the Incidence of Pertussis and Previous Meteorological Variables Based on the Theory of″Plague Transformation in Three Years″and Building Prediction Models

    ZHANG Xuan,HE Juan△
    Basic Medicine School of Beijing University of Chinese Medicine,Beijing 100029,China

    Objective:To explore the correlation between pertussis incidence in Beijing and previous(one year to three years ago)meteorological variables,and establish medical meteorological prediction models of BP artificial neural network.Methods:The data of climate change from 1970 to 2004 and pertussis incidence in Beijing were collected and meteorological prediction models were established by adopting the method of BP artificial neural network from three different time dimensions including one year,two years and three years ago respectively.Results:1)Pertussis incidence is related to previous(one year to three years ago)meteorological variables,,among them,closely related meteorological factor was average relative humidity(one year to three years ago);2)Prediction model of pertussis incidence could be successfully established by using previous meteorological variables,especially prediction models established on the foundation of meteorological one year ago showing the best prediction effects.Conclusion:Climates one to three years ago might be affect the incidence of some infectious diseases,and the study on the incidences of infectious diseases should pay attention to meteorological changes early.

    pertussis;BP artificial neural network;previous meteorological variables;plague transformation in three years;five movements and six climates

    R516.7

    A

    1004-6852(2015)11-0038-05

    2014-112-27

    國家自然科學基金項目(編號81072896)。

    張軒(1987—),女,在讀博士研究生。研究方向:五運六氣學說。

    △通訊作者:賀娟(1964—),女,博士學位,博士研究生導師,教授。研究方向:五運六氣學說。

    猜你喜歡
    百日咳北京地區(qū)時段
    百日咳要咳一百天嗎
    阿奇霉素在小兒百日咳的應用
    小心百日咳
    健康博覽(2019年3期)2019-12-02 05:42:16
    四個養(yǎng)生黃金時段,你抓住了嗎
    咳不停,是不是得了百日咳
    北京地區(qū)居民鎘攝入量評估
    傍晚是交通事故高發(fā)時段
    分時段預約在PICC門診維護中的應用與探討
    1949—1966年北京地區(qū)貞操觀的變革——解放與進步
    分時段預約掛號的實現(xiàn)與應用
    亚洲欧美精品自产自拍| 9热在线视频观看99| 欧美精品一区二区免费开放| 少妇精品久久久久久久| 一级黄片播放器| 中文字幕高清在线视频| 嫁个100分男人电影在线观看 | 久久av网站| 手机成人av网站| 国产精品秋霞免费鲁丝片| 国产又色又爽无遮挡免| 一二三四在线观看免费中文在| 久久午夜综合久久蜜桃| 纵有疾风起免费观看全集完整版| 日韩电影二区| 国产精品一二三区在线看| 国产日韩一区二区三区精品不卡| 美女主播在线视频| 成在线人永久免费视频| 亚洲av成人不卡在线观看播放网 | 少妇的丰满在线观看| 国产不卡av网站在线观看| 日韩大码丰满熟妇| 赤兔流量卡办理| 水蜜桃什么品种好| 国产精品一区二区免费欧美 | 国产一区二区激情短视频 | 欧美精品啪啪一区二区三区 | h视频一区二区三区| 欧美成人午夜精品| 亚洲情色 制服丝袜| 色婷婷av一区二区三区视频| www.999成人在线观看| 久久久久久免费高清国产稀缺| 国产在视频线精品| 国产激情久久老熟女| 日韩av不卡免费在线播放| 99久久人妻综合| a 毛片基地| 国产免费视频播放在线视频| 一本—道久久a久久精品蜜桃钙片| 黄频高清免费视频| 日韩 欧美 亚洲 中文字幕| 国产精品免费视频内射| 黄色毛片三级朝国网站| 各种免费的搞黄视频| 18在线观看网站| 人体艺术视频欧美日本| 久久狼人影院| 纯流量卡能插随身wifi吗| 高清欧美精品videossex| 亚洲欧洲日产国产| 亚洲精品乱久久久久久| 免费看十八禁软件| 精品久久久久久久毛片微露脸 | 制服人妻中文乱码| 亚洲精品久久久久久婷婷小说| 久久久国产精品麻豆| 久久精品aⅴ一区二区三区四区| 亚洲中文字幕日韩| 久久国产亚洲av麻豆专区| 又粗又硬又长又爽又黄的视频| 日韩av在线免费看完整版不卡| 人人妻人人爽人人添夜夜欢视频| 丝袜在线中文字幕| 大片免费播放器 马上看| 欧美日韩亚洲国产一区二区在线观看 | 亚洲激情五月婷婷啪啪| 国产精品香港三级国产av潘金莲 | 十八禁网站网址无遮挡| 免费看av在线观看网站| av在线播放精品| 各种免费的搞黄视频| 亚洲欧美一区二区三区黑人| 国产淫语在线视频| 色婷婷久久久亚洲欧美| av一本久久久久| av天堂久久9| 午夜91福利影院| 777久久人妻少妇嫩草av网站| 婷婷丁香在线五月| 两人在一起打扑克的视频| 成人亚洲欧美一区二区av| 欧美亚洲日本最大视频资源| xxx大片免费视频| 久久99精品国语久久久| 久久精品国产亚洲av高清一级| 性高湖久久久久久久久免费观看| 免费看不卡的av| av天堂在线播放| 秋霞在线观看毛片| 99九九在线精品视频| 久久综合国产亚洲精品| 国产成人啪精品午夜网站| 少妇人妻久久综合中文| 亚洲色图 男人天堂 中文字幕| 国产女主播在线喷水免费视频网站| 国产精品九九99| 精品人妻在线不人妻| 日韩av免费高清视频| 美女主播在线视频| 免费看十八禁软件| 成人亚洲精品一区在线观看| 一区二区日韩欧美中文字幕| 黄色 视频免费看| 亚洲欧美精品综合一区二区三区| 热99国产精品久久久久久7| 91精品国产国语对白视频| av天堂久久9| 自拍欧美九色日韩亚洲蝌蚪91| 日本一区二区免费在线视频| 精品国产乱码久久久久久小说| 欧美日韩视频精品一区| 别揉我奶头~嗯~啊~动态视频 | 亚洲国产精品999| 久久精品aⅴ一区二区三区四区| 在线观看免费视频网站a站| 一个人免费看片子| 999久久久国产精品视频| 一级黄片播放器| 一边摸一边做爽爽视频免费| 日韩一区二区三区影片| 亚洲熟女毛片儿| 国产日韩欧美视频二区| 一边摸一边抽搐一进一出视频| 久久国产亚洲av麻豆专区| 国产精品一区二区在线不卡| www.精华液| 中文精品一卡2卡3卡4更新| 自拍欧美九色日韩亚洲蝌蚪91| 国产成人91sexporn| 美女主播在线视频| 99香蕉大伊视频| 国产极品粉嫩免费观看在线| 国产精品秋霞免费鲁丝片| 久久久精品国产亚洲av高清涩受| 婷婷色综合大香蕉| 亚洲熟女精品中文字幕| 视频在线观看一区二区三区| 性色av乱码一区二区三区2| 性高湖久久久久久久久免费观看| 大香蕉久久网| 欧美日本中文国产一区发布| 捣出白浆h1v1| 91精品三级在线观看| 欧美人与善性xxx| 亚洲视频免费观看视频| 欧美变态另类bdsm刘玥| 国产精品久久久av美女十八| av线在线观看网站| 精品人妻1区二区| 丁香六月天网| av天堂在线播放| 91成人精品电影| 好男人电影高清在线观看| 黄色 视频免费看| 亚洲精品在线美女| 欧美成人午夜精品| 亚洲欧美激情在线| 精品国产超薄肉色丝袜足j| 一边摸一边做爽爽视频免费| 亚洲欧美一区二区三区久久| 午夜福利视频在线观看免费| 精品国产超薄肉色丝袜足j| 巨乳人妻的诱惑在线观看| 久久国产亚洲av麻豆专区| 国产成人av激情在线播放| 亚洲精品av麻豆狂野| 午夜两性在线视频| 亚洲精品成人av观看孕妇| 丝袜在线中文字幕| 黑人巨大精品欧美一区二区蜜桃| 亚洲伊人色综图| 亚洲成人免费av在线播放| 国产成人系列免费观看| 亚洲av在线观看美女高潮| 亚洲精品国产av蜜桃| 在线观看人妻少妇| 亚洲少妇的诱惑av| 亚洲国产精品国产精品| 精品亚洲成a人片在线观看| 久久久久久久国产电影| 三上悠亚av全集在线观看| 自线自在国产av| 日本91视频免费播放| 十八禁高潮呻吟视频| 乱人伦中国视频| 成年av动漫网址| 老熟女久久久| av在线app专区| 国产日韩欧美在线精品| 亚洲国产看品久久| av视频免费观看在线观看| 精品人妻1区二区| 999久久久国产精品视频| 日本91视频免费播放| 久久毛片免费看一区二区三区| 欧美老熟妇乱子伦牲交| 亚洲成色77777| 纵有疾风起免费观看全集完整版| 在线观看免费日韩欧美大片| 精品熟女少妇八av免费久了| 欧美黄色淫秽网站| 男女国产视频网站| 91精品国产国语对白视频| 一本久久精品| 亚洲午夜精品一区,二区,三区| 亚洲专区国产一区二区| 侵犯人妻中文字幕一二三四区| 美女视频免费永久观看网站| 三上悠亚av全集在线观看| av在线老鸭窝| 久久毛片免费看一区二区三区| 男女床上黄色一级片免费看| 看十八女毛片水多多多| 色网站视频免费| 99国产精品一区二区蜜桃av | 国产精品二区激情视频| 亚洲精品美女久久av网站| 国产成人精品在线电影| 大型av网站在线播放| 观看av在线不卡| 爱豆传媒免费全集在线观看| 欧美精品av麻豆av| 真人做人爱边吃奶动态| 国产老妇伦熟女老妇高清| 国产无遮挡羞羞视频在线观看| 一个人免费看片子| 欧美黄色片欧美黄色片| 叶爱在线成人免费视频播放| 亚洲一区二区三区欧美精品| 欧美日韩视频精品一区| av片东京热男人的天堂| 亚洲欧洲精品一区二区精品久久久| 黄色a级毛片大全视频| 美女福利国产在线| 欧美人与性动交α欧美软件| 男女国产视频网站| 黄色a级毛片大全视频| 一本大道久久a久久精品| 国产成人免费无遮挡视频| 国产成人免费无遮挡视频| 成人18禁高潮啪啪吃奶动态图| 两个人看的免费小视频| 99九九在线精品视频| 亚洲专区国产一区二区| 天天躁夜夜躁狠狠躁躁| 亚洲第一av免费看| 国产一卡二卡三卡精品| 国产熟女欧美一区二区| 亚洲美女黄色视频免费看| 性少妇av在线| 国产一区二区 视频在线| 丰满少妇做爰视频| 精品国产超薄肉色丝袜足j| 国产主播在线观看一区二区 | 777米奇影视久久| 亚洲,一卡二卡三卡| 国产老妇伦熟女老妇高清| 中文字幕亚洲精品专区| 青春草视频在线免费观看| 黄色一级大片看看| 欧美精品人与动牲交sv欧美| 日本wwww免费看| 国产精品 欧美亚洲| 日韩免费高清中文字幕av| 亚洲中文字幕日韩| 国产真人三级小视频在线观看| 中文字幕人妻丝袜一区二区| 高清黄色对白视频在线免费看| 亚洲欧洲国产日韩| 国产精品免费大片| 男男h啪啪无遮挡| 久久精品国产a三级三级三级| 午夜福利免费观看在线| 亚洲国产欧美一区二区综合| 少妇人妻 视频| 夫妻午夜视频| 亚洲成人手机| 国产日韩欧美视频二区| 国产又爽黄色视频| 国产成人精品在线电影| 青草久久国产| 男女边吃奶边做爰视频| 夜夜骑夜夜射夜夜干| 国产高清视频在线播放一区 | 国产在线一区二区三区精| 男女无遮挡免费网站观看| 免费高清在线观看日韩| 国产精品偷伦视频观看了| 欧美黄色淫秽网站| 国产激情久久老熟女| 美女扒开内裤让男人捅视频| 国产亚洲午夜精品一区二区久久| 操出白浆在线播放| 91国产中文字幕| 大话2 男鬼变身卡| 国产伦人伦偷精品视频| 99精国产麻豆久久婷婷| 欧美在线一区亚洲| 午夜福利,免费看| 欧美变态另类bdsm刘玥| 波多野结衣av一区二区av| 99久久人妻综合| 自线自在国产av| 亚洲国产欧美日韩在线播放| 九色亚洲精品在线播放| 一本久久精品| 最近手机中文字幕大全| 午夜日韩欧美国产| 中文字幕最新亚洲高清| 高清av免费在线| 高清av免费在线| 咕卡用的链子| 国产黄色视频一区二区在线观看| 国产爽快片一区二区三区| av欧美777| 久久人妻福利社区极品人妻图片 | 天天躁日日躁夜夜躁夜夜| 国产黄色视频一区二区在线观看| 国产在线一区二区三区精| 1024香蕉在线观看| 日本a在线网址| 美女福利国产在线| 日本色播在线视频| 最新在线观看一区二区三区 | 少妇猛男粗大的猛烈进出视频| 国产亚洲一区二区精品| 搡老岳熟女国产| 夫妻性生交免费视频一级片| 国产黄频视频在线观看| 久久久国产一区二区| 丰满少妇做爰视频| 三上悠亚av全集在线观看| 欧美激情 高清一区二区三区| 在线看a的网站| 一本—道久久a久久精品蜜桃钙片| 91精品伊人久久大香线蕉| 欧美日韩黄片免| 亚洲精品一二三| 免费人妻精品一区二区三区视频| √禁漫天堂资源中文www| 又黄又粗又硬又大视频| av在线app专区| 女人久久www免费人成看片| 久久免费观看电影| 成人国语在线视频| 日韩免费高清中文字幕av| 丝瓜视频免费看黄片| 亚洲国产日韩一区二区| cao死你这个sao货| 这个男人来自地球电影免费观看| 国产黄色免费在线视频| 国产成人影院久久av| 丝袜人妻中文字幕| 免费看十八禁软件| 欧美精品一区二区大全| 国产精品av久久久久免费| 美女脱内裤让男人舔精品视频| 国产色视频综合| 亚洲av成人精品一二三区| 日韩一卡2卡3卡4卡2021年| 99热国产这里只有精品6| 精品熟女少妇八av免费久了| 99香蕉大伊视频| 精品国产一区二区久久| 欧美日韩福利视频一区二区| 日本一区二区免费在线视频| 少妇人妻久久综合中文| 成人国产av品久久久| 精品一区二区三区av网在线观看 | 亚洲图色成人| 黄色视频不卡| 又大又黄又爽视频免费| 午夜福利视频精品| 久久人人97超碰香蕉20202| 成年人免费黄色播放视频| 婷婷色综合大香蕉| 啦啦啦在线免费观看视频4| 少妇粗大呻吟视频| 国产三级黄色录像| 精品国产一区二区久久| 9191精品国产免费久久| 国产伦人伦偷精品视频| 嫁个100分男人电影在线观看 | 伊人久久大香线蕉亚洲五| 久久精品人人爽人人爽视色| 黄色片一级片一级黄色片| 大片免费播放器 马上看| 国产精品麻豆人妻色哟哟久久| 成年美女黄网站色视频大全免费| 国产97色在线日韩免费| 国产精品熟女久久久久浪| 精品一区二区三卡| 99国产精品一区二区蜜桃av | 又大又黄又爽视频免费| 久久精品久久精品一区二区三区| 另类精品久久| 女性生殖器流出的白浆| 最新的欧美精品一区二区| 亚洲欧洲精品一区二区精品久久久| 久久久久国产精品人妻一区二区| 蜜桃在线观看..| 久9热在线精品视频| 男男h啪啪无遮挡| 少妇粗大呻吟视频| 精品一品国产午夜福利视频| 777久久人妻少妇嫩草av网站| 一本大道久久a久久精品| 黄色a级毛片大全视频| 久久99一区二区三区| 日韩欧美一区视频在线观看| 日本黄色日本黄色录像| 日本色播在线视频| 国产1区2区3区精品| 在线观看人妻少妇| 亚洲av欧美aⅴ国产| 欧美精品av麻豆av| 老司机靠b影院| 国产1区2区3区精品| 嫩草影视91久久| 肉色欧美久久久久久久蜜桃| 成年女人毛片免费观看观看9 | 国产一区二区三区av在线| 亚洲欧美日韩高清在线视频 | 十八禁人妻一区二区| 人人妻人人添人人爽欧美一区卜| 悠悠久久av| 免费女性裸体啪啪无遮挡网站| 国产亚洲午夜精品一区二区久久| 男人爽女人下面视频在线观看| 美女中出高潮动态图| 中文欧美无线码| 久久精品久久精品一区二区三区| 午夜福利视频在线观看免费| 国产av精品麻豆| 国产男女内射视频| 久久精品久久久久久久性| www.av在线官网国产| 七月丁香在线播放| 老司机午夜十八禁免费视频| 99国产精品99久久久久| 国产一区有黄有色的免费视频| 亚洲av国产av综合av卡| 丁香六月欧美| 欧美精品高潮呻吟av久久| 国产精品成人在线| 99九九在线精品视频| 两个人免费观看高清视频| 久9热在线精品视频| 美国免费a级毛片| 黄片小视频在线播放| av片东京热男人的天堂| 纵有疾风起免费观看全集完整版| 欧美激情极品国产一区二区三区| 日韩大片免费观看网站| 大片免费播放器 马上看| 亚洲精品国产区一区二| 久久中文字幕一级| 亚洲国产av新网站| 一级毛片我不卡| 久久久欧美国产精品| 一区二区av电影网| 水蜜桃什么品种好| 国产亚洲欧美在线一区二区| 大码成人一级视频| 老司机靠b影院| 中国美女看黄片| 咕卡用的链子| 免费看不卡的av| 国产一区二区在线观看av| a级毛片黄视频| 一二三四社区在线视频社区8| 极品人妻少妇av视频| 下体分泌物呈黄色| 国产精品二区激情视频| a 毛片基地| 另类亚洲欧美激情| 精品一区二区三区四区五区乱码 | 中文乱码字字幕精品一区二区三区| 91精品三级在线观看| 18禁裸乳无遮挡动漫免费视频| 久热爱精品视频在线9| 久久久国产精品麻豆| 国产男女内射视频| 晚上一个人看的免费电影| 人成视频在线观看免费观看| 久久精品久久精品一区二区三区| 亚洲国产精品国产精品| 免费看十八禁软件| 好男人视频免费观看在线| 日韩 欧美 亚洲 中文字幕| 免费观看a级毛片全部| 不卡av一区二区三区| 久久久欧美国产精品| 两人在一起打扑克的视频| 97精品久久久久久久久久精品| 色94色欧美一区二区| 久久久久精品国产欧美久久久 | 一级片'在线观看视频| 久久精品久久久久久久性| 狠狠婷婷综合久久久久久88av| 青青草视频在线视频观看| 国产又爽黄色视频| 操美女的视频在线观看| 亚洲av综合色区一区| 满18在线观看网站| 丁香六月欧美| 另类亚洲欧美激情| 少妇人妻久久综合中文| 国产精品一二三区在线看| 香蕉丝袜av| 亚洲精品久久午夜乱码| 成人午夜精彩视频在线观看| 一级,二级,三级黄色视频| 女人爽到高潮嗷嗷叫在线视频| 深夜精品福利| 国产激情久久老熟女| 欧美激情高清一区二区三区| 亚洲精品国产色婷婷电影| 国产日韩欧美亚洲二区| 久久精品成人免费网站| 午夜日韩欧美国产| 91麻豆精品激情在线观看国产 | 国产精品一区二区精品视频观看| 性高湖久久久久久久久免费观看| 国产亚洲午夜精品一区二区久久| av天堂在线播放| 久久性视频一级片| 亚洲欧美中文字幕日韩二区| 日韩熟女老妇一区二区性免费视频| 国产精品 欧美亚洲| 国产精品麻豆人妻色哟哟久久| 久久精品久久久久久噜噜老黄| 一级黄色大片毛片| 欧美xxⅹ黑人| 好男人视频免费观看在线| 99国产精品免费福利视频| 一级毛片电影观看| 在线观看一区二区三区激情| 精品国产一区二区久久| 国产成人欧美| 后天国语完整版免费观看| 精品福利观看| 亚洲国产精品一区二区三区在线| 国产免费一区二区三区四区乱码| 中文字幕人妻丝袜一区二区| 在线观看国产h片| 一区二区日韩欧美中文字幕| 又粗又硬又长又爽又黄的视频| 国产野战对白在线观看| 精品国产一区二区久久| 国产高清视频在线播放一区 | 亚洲黑人精品在线| 婷婷色综合www| 国产精品av久久久久免费| 99九九在线精品视频| 午夜福利免费观看在线| 高清不卡的av网站| 老汉色av国产亚洲站长工具| 国产日韩一区二区三区精品不卡| 两个人免费观看高清视频| 高清av免费在线| 日本欧美视频一区| 亚洲成国产人片在线观看| 亚洲少妇的诱惑av| 99久久精品国产亚洲精品| 男人舔女人的私密视频| 操出白浆在线播放| 每晚都被弄得嗷嗷叫到高潮| 日本av免费视频播放| 妹子高潮喷水视频| 亚洲美女黄色视频免费看| 久久久精品区二区三区| 成人手机av| 热re99久久精品国产66热6| 久久99精品国语久久久| 黄网站色视频无遮挡免费观看| 十八禁人妻一区二区| 99热全是精品| 91国产中文字幕| 一区二区三区激情视频| 久久精品亚洲av国产电影网| 九色亚洲精品在线播放| av网站免费在线观看视频| 亚洲欧美精品综合一区二区三区| 亚洲激情五月婷婷啪啪| 国产亚洲av高清不卡| 只有这里有精品99| 婷婷成人精品国产| 亚洲色图 男人天堂 中文字幕| 午夜两性在线视频| 亚洲人成电影观看| 久久久久久人人人人人| 亚洲伊人久久精品综合| 人人妻,人人澡人人爽秒播 | 成人三级做爰电影| 在线观看人妻少妇| 少妇人妻久久综合中文| 美女视频免费永久观看网站| 久久天躁狠狠躁夜夜2o2o | 亚洲视频免费观看视频| 午夜久久久在线观看| 一个人免费看片子| 两人在一起打扑克的视频| 亚洲国产欧美日韩在线播放| 国产成人影院久久av| 老司机在亚洲福利影院| 你懂的网址亚洲精品在线观看| 黑人欧美特级aaaaaa片| 国产不卡av网站在线观看| 啦啦啦在线观看免费高清www| 午夜福利视频在线观看免费| 操美女的视频在线观看| 老熟女久久久|