• <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)與應用
    国产黄频视频在线观看| 亚洲精品一二三| 寂寞人妻少妇视频99o| 精品酒店卫生间| 国产午夜精品久久久久久一区二区三区| 亚洲精品乱久久久久久| 高清欧美精品videossex| 成人亚洲精品av一区二区| 91精品一卡2卡3卡4卡| 一区二区三区高清视频在线| 日韩,欧美,国产一区二区三区| 亚洲精品456在线播放app| 99久久精品一区二区三区| 一级毛片黄色毛片免费观看视频| 草草在线视频免费看| 免费大片18禁| 一级av片app| 国产成人午夜福利电影在线观看| 麻豆久久精品国产亚洲av| 日韩制服骚丝袜av| 久久久久九九精品影院| 欧美97在线视频| 亚洲成人av在线免费| 黄片wwwwww| 观看美女的网站| a级毛色黄片| 亚洲av中文字字幕乱码综合| 日产精品乱码卡一卡2卡三| 久久久精品94久久精品| 国产黄色视频一区二区在线观看| 成人美女网站在线观看视频| 日本-黄色视频高清免费观看| 久久这里只有精品中国| 美女cb高潮喷水在线观看| 久久久久久久国产电影| 黄色一级大片看看| 亚洲av.av天堂| 国产精品久久视频播放| 亚洲av成人精品一区久久| 看十八女毛片水多多多| 久久精品久久久久久噜噜老黄| 最近中文字幕高清免费大全6| 国产精品一二三区在线看| 欧美变态另类bdsm刘玥| 国产一区二区三区av在线| 超碰av人人做人人爽久久| 亚洲一区高清亚洲精品| 久久精品国产亚洲av涩爱| 最近中文字幕2019免费版| 十八禁网站网址无遮挡 | 99久久精品国产国产毛片| 中文在线观看免费www的网站| 国产精品一区二区在线观看99 | 在现免费观看毛片| 天堂网av新在线| 午夜福利在线观看免费完整高清在| 国产精品国产三级国产av玫瑰| 国模一区二区三区四区视频| 国产一区有黄有色的免费视频 | 校园人妻丝袜中文字幕| 成年免费大片在线观看| 亚洲色图av天堂| 欧美bdsm另类| 亚洲va在线va天堂va国产| 在线观看人妻少妇| 国产爱豆传媒在线观看| 亚洲av一区综合| 国产亚洲一区二区精品| 亚洲色图av天堂| 看非洲黑人一级黄片| 欧美zozozo另类| 女的被弄到高潮叫床怎么办| 国产老妇女一区| 色吧在线观看| 大香蕉97超碰在线| 国产伦精品一区二区三区四那| 亚洲精品aⅴ在线观看| 有码 亚洲区| 亚洲精品第二区| 91精品一卡2卡3卡4卡| 亚洲真实伦在线观看| 最近视频中文字幕2019在线8| 久久精品国产亚洲网站| 久久久久久久久久久免费av| 精品久久久久久久久亚洲| 亚洲精品自拍成人| 久久97久久精品| 中文字幕免费在线视频6| 免费无遮挡裸体视频| 国产精品av视频在线免费观看| 亚洲国产精品成人综合色| 国产精品久久久久久久久免| 色综合色国产| 亚洲熟妇中文字幕五十中出| 国产精品一区www在线观看| 日韩欧美一区视频在线观看 | 非洲黑人性xxxx精品又粗又长| 国产成人午夜福利电影在线观看| 欧美日韩精品成人综合77777| 一区二区三区四区激情视频| 又黄又爽又刺激的免费视频.| 成年女人看的毛片在线观看| 久久久久久国产a免费观看| 亚洲欧美一区二区三区国产| 九草在线视频观看| 草草在线视频免费看| 成人无遮挡网站| 观看免费一级毛片| 亚洲精品乱码久久久久久按摩| 综合色av麻豆| 久久久久久久国产电影| 成人毛片60女人毛片免费| 有码 亚洲区| 国产精品国产三级国产专区5o| 内射极品少妇av片p| 国产一区二区在线观看日韩| 老司机影院成人| 啦啦啦韩国在线观看视频| 日韩av在线免费看完整版不卡| 久久综合国产亚洲精品| 精品久久久久久久久久久久久| 国产有黄有色有爽视频| 五月玫瑰六月丁香| 99热全是精品| 欧美不卡视频在线免费观看| 亚洲国产欧美在线一区| 熟女人妻精品中文字幕| 老女人水多毛片| 亚洲不卡免费看| 免费观看精品视频网站| 18禁裸乳无遮挡免费网站照片| 男女边摸边吃奶| 精品午夜福利在线看| 国产精品女同一区二区软件| 精品一区二区三区视频在线| 99热这里只有是精品50| 人妻一区二区av| 亚洲成色77777| 午夜福利在线在线| 亚洲精品色激情综合| 日日干狠狠操夜夜爽| 免费黄色在线免费观看| 丝袜喷水一区| 黄色配什么色好看| 在线观看av片永久免费下载| av免费在线看不卡| 69av精品久久久久久| 国产精品一二三区在线看| 免费人成在线观看视频色| 在线a可以看的网站| 黄片wwwwww| 国产精品女同一区二区软件| 精品99又大又爽又粗少妇毛片| 夫妻性生交免费视频一级片| 狂野欧美激情性xxxx在线观看| 蜜桃亚洲精品一区二区三区| 天堂中文最新版在线下载 | 欧美bdsm另类| 五月天丁香电影| 日韩在线高清观看一区二区三区| 少妇被粗大猛烈的视频| 久久97久久精品| 少妇熟女欧美另类| 免费看av在线观看网站| 国产精品一区www在线观看| 麻豆av噜噜一区二区三区| 久久久精品免费免费高清| 夜夜看夜夜爽夜夜摸| 男人舔女人下体高潮全视频| 天美传媒精品一区二区| 建设人人有责人人尽责人人享有的 | 久久草成人影院| 中文资源天堂在线| 国产极品天堂在线| 一级a做视频免费观看| 成人av在线播放网站| 成人一区二区视频在线观看| 欧美bdsm另类| 久99久视频精品免费| 国产精品99久久久久久久久| 美女黄网站色视频| 特大巨黑吊av在线直播| 亚洲欧美一区二区三区黑人 | 在线观看人妻少妇| 久久久午夜欧美精品| 日韩av在线大香蕉| 国产成人精品一,二区| 国产精品国产三级专区第一集| 久久韩国三级中文字幕| 成人午夜精彩视频在线观看| 国产精品一及| 亚洲一区高清亚洲精品| 边亲边吃奶的免费视频| 美女内射精品一级片tv| 国产午夜精品一二区理论片| 国产乱人视频| 少妇的逼好多水| 亚洲aⅴ乱码一区二区在线播放| 麻豆乱淫一区二区| 精品人妻一区二区三区麻豆| 夜夜看夜夜爽夜夜摸| 蜜臀久久99精品久久宅男| 国产成年人精品一区二区| 成人无遮挡网站| 久久精品国产亚洲av天美| 大香蕉久久网| 日韩欧美精品v在线| 乱人视频在线观看| 少妇熟女欧美另类| 18+在线观看网站| 亚洲激情五月婷婷啪啪| 久久久久久久午夜电影| 亚洲欧美中文字幕日韩二区| 欧美97在线视频| 国产女主播在线喷水免费视频网站 | 啦啦啦中文免费视频观看日本| 亚洲人成网站高清观看| 亚洲熟女精品中文字幕| 一本一本综合久久| 色综合色国产| 国产精品三级大全| 乱码一卡2卡4卡精品| av在线观看视频网站免费| 好男人视频免费观看在线| 国产午夜精品一二区理论片| 亚州av有码| 国产单亲对白刺激| 久久久亚洲精品成人影院| 国产免费视频播放在线视频 | 一夜夜www| 午夜精品一区二区三区免费看| 少妇人妻一区二区三区视频| 三级毛片av免费| 在线a可以看的网站| 青春草亚洲视频在线观看| 在线观看人妻少妇| av在线观看视频网站免费| 成人午夜高清在线视频| 女人被狂操c到高潮| 午夜精品一区二区三区免费看| 日韩,欧美,国产一区二区三区| 亚洲欧美一区二区三区国产| 欧美激情在线99| 欧美 日韩 精品 国产| 国产精品久久久久久av不卡| 欧美 日韩 精品 国产| 伦理电影大哥的女人| 亚洲欧美成人综合另类久久久| 午夜福利视频精品| 亚洲不卡免费看| 永久网站在线| 一级毛片电影观看| 亚洲一区高清亚洲精品| 18+在线观看网站| 久久99精品国语久久久| 成人性生交大片免费视频hd| 免费av不卡在线播放| 乱人视频在线观看| 网址你懂的国产日韩在线| 久久久久久国产a免费观看| 尤物成人国产欧美一区二区三区| 国产午夜精品久久久久久一区二区三区| 亚洲最大成人手机在线| 在线a可以看的网站| 在线观看美女被高潮喷水网站| 精品久久久噜噜| 国内揄拍国产精品人妻在线| 人妻系列 视频| 欧美日韩综合久久久久久| 国产综合精华液| 国产淫语在线视频| 日韩制服骚丝袜av| 久久久久久久国产电影| 2021天堂中文幕一二区在线观| 好男人在线观看高清免费视频| 国产黄片美女视频| 免费观看性生交大片5| 少妇的逼水好多| 亚洲精品色激情综合| 亚洲av不卡在线观看| 久久久午夜欧美精品| 五月天丁香电影| 成人亚洲欧美一区二区av| 一夜夜www| 国产成人午夜福利电影在线观看| 九九久久精品国产亚洲av麻豆| 禁无遮挡网站| 99热这里只有精品一区| 午夜免费男女啪啪视频观看| 国产免费福利视频在线观看| 777米奇影视久久| 三级毛片av免费| 免费看光身美女| 草草在线视频免费看| 男人和女人高潮做爰伦理| 精品久久久久久久人妻蜜臀av| 街头女战士在线观看网站| 搡老妇女老女人老熟妇| 国产精品蜜桃在线观看| 国产成人freesex在线| 国产精品一及| 最近2019中文字幕mv第一页| 特级一级黄色大片| 美女黄网站色视频| 99热这里只有是精品在线观看| 亚洲成色77777| 久久精品综合一区二区三区| 午夜免费观看性视频| 国产成年人精品一区二区| 亚洲精品一区蜜桃| 国产精品久久久久久精品电影| 中文字幕制服av| 国产精品无大码| 99热全是精品| 日产精品乱码卡一卡2卡三| 简卡轻食公司| 精品人妻偷拍中文字幕| 最近2019中文字幕mv第一页| 男女边摸边吃奶| 亚洲国产色片| 天堂中文最新版在线下载 | 嫩草影院新地址| 国产成人午夜福利电影在线观看| 2018国产大陆天天弄谢| 成人综合一区亚洲| 日韩av免费高清视频| 亚洲在线观看片| 国产 亚洲一区二区三区 | 国产v大片淫在线免费观看| 国产三级在线视频| 婷婷色av中文字幕| 亚洲成人中文字幕在线播放| 美女高潮的动态| 国产极品天堂在线| 成年免费大片在线观看| 黄片无遮挡物在线观看| 蜜桃亚洲精品一区二区三区| 汤姆久久久久久久影院中文字幕 | 亚洲精品国产av蜜桃| av福利片在线观看| 色哟哟·www| 麻豆av噜噜一区二区三区| 黄色一级大片看看| 久久久久免费精品人妻一区二区| 午夜免费男女啪啪视频观看| 日韩伦理黄色片| 久久久久网色| 久久鲁丝午夜福利片| 一级片'在线观看视频| 99久国产av精品| 亚洲欧美精品自产自拍| 国产黄色小视频在线观看| av在线老鸭窝| 黄片无遮挡物在线观看| 日本一二三区视频观看| 日韩,欧美,国产一区二区三区| 高清av免费在线| 熟妇人妻不卡中文字幕| av网站免费在线观看视频 | 久久久久国产网址| 搡老妇女老女人老熟妇| 欧美不卡视频在线免费观看| 美女内射精品一级片tv| 高清午夜精品一区二区三区| 久久久色成人| 日韩一本色道免费dvd| 最近2019中文字幕mv第一页| 免费高清在线观看视频在线观看| 美女高潮的动态| 中文欧美无线码| 免费看av在线观看网站| 国产不卡一卡二| 国产91av在线免费观看| 美女国产视频在线观看| 成人美女网站在线观看视频| 国产精品三级大全| 亚洲精华国产精华液的使用体验| 在线观看美女被高潮喷水网站| 国产精品99久久久久久久久| 搡老妇女老女人老熟妇| 白带黄色成豆腐渣| 国产在视频线精品| 日韩欧美 国产精品| 精品人妻熟女av久视频| 免费看a级黄色片| 超碰97精品在线观看| av免费观看日本| 中文字幕亚洲精品专区| 色综合亚洲欧美另类图片| 极品少妇高潮喷水抽搐| 国产男人的电影天堂91| 国产午夜精品一二区理论片| 少妇裸体淫交视频免费看高清| 国产探花极品一区二区| 嫩草影院入口| 国产av码专区亚洲av| 2021少妇久久久久久久久久久| 亚洲av日韩在线播放| 简卡轻食公司| 国产三级在线视频| 国产亚洲一区二区精品| 边亲边吃奶的免费视频| 听说在线观看完整版免费高清| 国产精品嫩草影院av在线观看| 欧美潮喷喷水| 午夜精品在线福利| 观看免费一级毛片| 久久精品久久久久久噜噜老黄| 最近的中文字幕免费完整| 美女脱内裤让男人舔精品视频| 一个人免费在线观看电影| 美女被艹到高潮喷水动态| 国产精品久久久久久精品电影| 国产69精品久久久久777片| 国产免费一级a男人的天堂| 大片免费播放器 马上看| 特级一级黄色大片| 久久久久久国产a免费观看| 免费黄频网站在线观看国产| 又大又黄又爽视频免费| 国产69精品久久久久777片| 国内少妇人妻偷人精品xxx网站| 亚洲欧美中文字幕日韩二区| 午夜日本视频在线| 久热久热在线精品观看| 日日撸夜夜添| 久久这里有精品视频免费| 我要看日韩黄色一级片| 日本熟妇午夜| 一本久久精品| 日韩一本色道免费dvd| 国产亚洲一区二区精品| 亚洲美女视频黄频| 国产色婷婷99| 亚洲av电影在线观看一区二区三区 | 特级一级黄色大片| 精品久久久久久久末码| 国产女主播在线喷水免费视频网站 | 久久久久久久久久人人人人人人| 好男人在线观看高清免费视频| 亚洲激情五月婷婷啪啪| 晚上一个人看的免费电影| 精品一区二区免费观看| 国产淫语在线视频| 久久综合国产亚洲精品| 国产白丝娇喘喷水9色精品| 成人无遮挡网站| 永久免费av网站大全| 国产精品国产三级国产专区5o| 欧美日韩精品成人综合77777| 日日啪夜夜爽| 国产一区二区亚洲精品在线观看| 一个人看视频在线观看www免费| 18禁在线无遮挡免费观看视频| 寂寞人妻少妇视频99o| 97精品久久久久久久久久精品| 亚洲熟女精品中文字幕| 亚洲天堂国产精品一区在线| 成人特级av手机在线观看| 熟妇人妻久久中文字幕3abv| 最后的刺客免费高清国语| 欧美精品一区二区大全| 精品午夜福利在线看| 久久精品久久精品一区二区三区| 三级男女做爰猛烈吃奶摸视频| 好男人在线观看高清免费视频| 久久久午夜欧美精品| 亚洲av电影不卡..在线观看| 亚洲自偷自拍三级| 中国美白少妇内射xxxbb| 丝袜美腿在线中文| av在线蜜桃| 天堂俺去俺来也www色官网 | 国产乱来视频区| 26uuu在线亚洲综合色| 男女那种视频在线观看| 最近视频中文字幕2019在线8| 国产成人精品福利久久| 日韩亚洲欧美综合| 91aial.com中文字幕在线观看| 直男gayav资源| 床上黄色一级片| 欧美激情国产日韩精品一区| 18禁裸乳无遮挡免费网站照片| 国产乱人视频| 国产91av在线免费观看| 一级片'在线观看视频| 国产成人一区二区在线| 嫩草影院入口| 91久久精品电影网| 欧美日韩亚洲高清精品| 亚洲av成人av| 亚洲欧美日韩东京热| 春色校园在线视频观看| 人人妻人人澡欧美一区二区| 高清在线视频一区二区三区| 一级毛片电影观看| 成年版毛片免费区| 能在线免费观看的黄片| 日韩伦理黄色片| 精品久久久噜噜| 久久综合国产亚洲精品| 在线天堂最新版资源| 在线免费观看不下载黄p国产| 久久精品久久久久久噜噜老黄| 午夜精品一区二区三区免费看| 偷拍熟女少妇极品色| 精品国产一区二区三区久久久樱花 | 亚洲国产精品sss在线观看| 国产大屁股一区二区在线视频| 91久久精品国产一区二区三区| 免费无遮挡裸体视频| 国产人妻一区二区三区在| 欧美日韩精品成人综合77777| 99热这里只有是精品50| 观看美女的网站| 色综合色国产| 一区二区三区免费毛片| 成人亚洲欧美一区二区av| 观看美女的网站| 尾随美女入室| 噜噜噜噜噜久久久久久91| 日本wwww免费看| 好男人视频免费观看在线| 久久久久久久久久人人人人人人| 18禁在线无遮挡免费观看视频| 国产男女超爽视频在线观看| 亚洲国产最新在线播放| 一级毛片黄色毛片免费观看视频| 亚洲精品国产av成人精品| 免费观看av网站的网址| 日韩亚洲欧美综合| 男的添女的下面高潮视频| 国产黄色视频一区二区在线观看| 国产黄色免费在线视频| 韩国高清视频一区二区三区| 免费观看a级毛片全部| 久久久久久久国产电影| 国产精品久久久久久精品电影| 国产亚洲av片在线观看秒播厂 | 午夜精品在线福利| 91久久精品电影网| 99热全是精品| 国产 一区精品| 99久久人妻综合| 国产精品99久久久久久久久| 亚洲精品亚洲一区二区| 国产v大片淫在线免费观看| www.色视频.com| 亚洲丝袜综合中文字幕| 午夜亚洲福利在线播放| 伦理电影大哥的女人| 久久国内精品自在自线图片| 国产午夜福利久久久久久| 国产中年淑女户外野战色| 婷婷色综合大香蕉| 欧美成人午夜免费资源| 我的女老师完整版在线观看| 国产免费又黄又爽又色| 丰满乱子伦码专区| 十八禁网站网址无遮挡 | 亚洲美女搞黄在线观看| 日韩欧美三级三区| 人人妻人人看人人澡| 欧美zozozo另类| 免费观看在线日韩| 韩国av在线不卡| 日韩精品有码人妻一区| 日韩视频在线欧美| 偷拍熟女少妇极品色| 国产白丝娇喘喷水9色精品| 亚洲欧美清纯卡通| 亚洲精品自拍成人| 欧美极品一区二区三区四区| 最近最新中文字幕免费大全7| 少妇的逼水好多| 久久6这里有精品| 午夜福利在线观看免费完整高清在| 麻豆乱淫一区二区| 日本wwww免费看| 啦啦啦韩国在线观看视频| 青青草视频在线视频观看| 亚洲av二区三区四区| 成年人午夜在线观看视频 | 国产一级毛片七仙女欲春2| 亚洲欧洲日产国产| 日韩在线高清观看一区二区三区| 综合色丁香网| 亚洲av二区三区四区| 欧美日韩国产mv在线观看视频 | 欧美3d第一页| 日本wwww免费看| 午夜精品一区二区三区免费看| 综合色av麻豆| 岛国毛片在线播放| 午夜福利视频1000在线观看| 日韩一本色道免费dvd| 99久国产av精品| 丰满人妻一区二区三区视频av| 午夜视频国产福利| 久久97久久精品| 欧美三级亚洲精品| 日韩一区二区三区影片| 久久精品夜夜夜夜夜久久蜜豆| 99热这里只有是精品50| 搡老乐熟女国产| 国产白丝娇喘喷水9色精品| 99久国产av精品国产电影| 日日啪夜夜撸| 精品99又大又爽又粗少妇毛片| 五月玫瑰六月丁香| 久久久久免费精品人妻一区二区| 一二三四中文在线观看免费高清| 亚洲国产高清在线一区二区三| 麻豆乱淫一区二区|