• <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)與應用
    亚洲成人手机| 久久人妻福利社区极品人妻图片| 操出白浆在线播放| 欧美国产精品一级二级三级| 免费日韩欧美在线观看| 啦啦啦免费观看视频1| 大香蕉久久网| 欧美人与性动交α欧美精品济南到| 亚洲精品国产色婷婷电影| 好看av亚洲va欧美ⅴa在| 悠悠久久av| 精品国产一区二区三区四区第35| 天天操日日干夜夜撸| 人人妻人人澡人人看| 国产片内射在线| 欧美精品啪啪一区二区三区| 午夜影院日韩av| 久久香蕉国产精品| 99久久国产精品久久久| 怎么达到女性高潮| 国产成人欧美| 日韩三级视频一区二区三区| 精品熟女少妇八av免费久了| 91大片在线观看| 丰满迷人的少妇在线观看| 韩国av一区二区三区四区| 一进一出抽搐gif免费好疼 | 伊人久久大香线蕉亚洲五| 一区福利在线观看| 日韩三级视频一区二区三区| 国产精品美女特级片免费视频播放器 | 日韩 欧美 亚洲 中文字幕| 国产一区二区三区在线臀色熟女 | xxxhd国产人妻xxx| 美女视频免费永久观看网站| 深夜精品福利| 一区在线观看完整版| 国产aⅴ精品一区二区三区波| 又黄又爽又免费观看的视频| 午夜福利一区二区在线看| 亚洲精品中文字幕一二三四区| 日本撒尿小便嘘嘘汇集6| 91精品国产国语对白视频| 国产一区二区三区视频了| 高清毛片免费观看视频网站 | 国产高清videossex| 亚洲国产精品合色在线| 国产高清激情床上av| 亚洲欧美一区二区三区黑人| 少妇的丰满在线观看| 精品国产国语对白av| 中文字幕制服av| 91精品三级在线观看| 欧美日韩成人在线一区二区| 午夜免费鲁丝| 大型黄色视频在线免费观看| 黄色视频不卡| 日本vs欧美在线观看视频| 少妇猛男粗大的猛烈进出视频| 国产一区二区三区视频了| 日本精品一区二区三区蜜桃| 在线观看免费午夜福利视频| 99久久99久久久精品蜜桃| 嫩草影视91久久| 女人久久www免费人成看片| 男人的好看免费观看在线视频 | 国产精品1区2区在线观看. | 夫妻午夜视频| 一区二区三区激情视频| 久久香蕉国产精品| 国产精品美女特级片免费视频播放器 | 在线观看免费日韩欧美大片| 精品国内亚洲2022精品成人 | 一级片免费观看大全| 欧美日韩亚洲国产一区二区在线观看 | 十八禁高潮呻吟视频| 女性生殖器流出的白浆| 十分钟在线观看高清视频www| 热re99久久国产66热| 国产在线一区二区三区精| 视频区欧美日本亚洲| 一级作爱视频免费观看| 丁香六月欧美| 亚洲美女黄片视频| 18禁美女被吸乳视频| 色婷婷av一区二区三区视频| 欧美精品啪啪一区二区三区| 黄片小视频在线播放| 久久九九热精品免费| 黑人欧美特级aaaaaa片| 久久人人97超碰香蕉20202| 五月开心婷婷网| 日韩有码中文字幕| 国产成人免费无遮挡视频| 91字幕亚洲| 成人亚洲精品一区在线观看| 99久久综合精品五月天人人| 黑人巨大精品欧美一区二区mp4| 人妻 亚洲 视频| 亚洲欧美一区二区三区黑人| 国产精品一区二区免费欧美| 成人国产一区最新在线观看| 亚洲精品中文字幕一二三四区| 欧美成人免费av一区二区三区 | 国产成人精品在线电影| 亚洲精品国产一区二区精华液| 激情在线观看视频在线高清 | 久久青草综合色| 一进一出抽搐gif免费好疼 | 18禁裸乳无遮挡动漫免费视频| 欧美人与性动交α欧美软件| 亚洲视频免费观看视频| 欧美精品一区二区免费开放| 亚洲av成人一区二区三| 午夜福利欧美成人| 国产精品偷伦视频观看了| 久久人人爽av亚洲精品天堂| 免费看a级黄色片| 免费观看人在逋| 成人三级做爰电影| 身体一侧抽搐| 欧美日韩乱码在线| 成熟少妇高潮喷水视频| 亚洲欧美精品综合一区二区三区| 十八禁网站免费在线| 人人妻人人澡人人看| 美国免费a级毛片| 亚洲国产毛片av蜜桃av| 一区二区三区精品91| 看免费av毛片| 久久国产精品男人的天堂亚洲| 69av精品久久久久久| 高清av免费在线| 国产精品99久久99久久久不卡| www.自偷自拍.com| 成在线人永久免费视频| 夫妻午夜视频| 伊人久久大香线蕉亚洲五| 欧美日韩乱码在线| 久久久久精品国产欧美久久久| 精品国产超薄肉色丝袜足j| 久久香蕉国产精品| 可以免费在线观看a视频的电影网站| 最近最新中文字幕大全电影3 | 午夜影院日韩av| 满18在线观看网站| 91老司机精品| 韩国精品一区二区三区| 国产精品久久久久久人妻精品电影| svipshipincom国产片| 在线观看www视频免费| 黄色视频不卡| 免费久久久久久久精品成人欧美视频| 少妇粗大呻吟视频| 亚洲黑人精品在线| 丰满迷人的少妇在线观看| 成人特级黄色片久久久久久久| 亚洲精品久久成人aⅴ小说| 满18在线观看网站| 精品久久久久久电影网| 精品国产国语对白av| 日韩 欧美 亚洲 中文字幕| 亚洲人成电影观看| av不卡在线播放| 成年动漫av网址| 亚洲五月婷婷丁香| 好看av亚洲va欧美ⅴa在| 午夜福利乱码中文字幕| 青草久久国产| 亚洲国产欧美网| 国产精品影院久久| 99re6热这里在线精品视频| 捣出白浆h1v1| 免费人成视频x8x8入口观看| 国产欧美日韩综合在线一区二区| 国产精品久久电影中文字幕 | 中文亚洲av片在线观看爽 | 18禁裸乳无遮挡动漫免费视频| 夜夜躁狠狠躁天天躁| 一级毛片精品| a级毛片黄视频| 下体分泌物呈黄色| 丝瓜视频免费看黄片| 午夜福利乱码中文字幕| 精品人妻熟女毛片av久久网站| cao死你这个sao货| 可以免费在线观看a视频的电影网站| 国产精品99久久99久久久不卡| 欧美日本中文国产一区发布| 欧美成人免费av一区二区三区 | 电影成人av| 人人澡人人妻人| 婷婷精品国产亚洲av在线 | 国产高清视频在线播放一区| cao死你这个sao货| 欧美日韩国产mv在线观看视频| 精品一品国产午夜福利视频| 国产成人av教育| 女性被躁到高潮视频| 欧美激情 高清一区二区三区| 操出白浆在线播放| 麻豆乱淫一区二区| 午夜亚洲福利在线播放| 免费在线观看日本一区| 另类亚洲欧美激情| 狂野欧美激情性xxxx| 天堂俺去俺来也www色官网| 天堂√8在线中文| 亚洲精品乱久久久久久| 国产高清视频在线播放一区| 99riav亚洲国产免费| 黑人操中国人逼视频| av欧美777| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美精品亚洲一区二区| 国产有黄有色有爽视频| 日日摸夜夜添夜夜添小说| 99久久人妻综合| 精品福利永久在线观看| 日韩免费高清中文字幕av| 午夜福利,免费看| aaaaa片日本免费| 一边摸一边做爽爽视频免费| 久久中文字幕人妻熟女| 欧美日韩亚洲国产一区二区在线观看 | 一区二区日韩欧美中文字幕| 久久这里只有精品19| 在线观看66精品国产| 久久久久久久午夜电影 | 少妇 在线观看| 国产欧美亚洲国产| 久久人妻av系列| 极品少妇高潮喷水抽搐| 十八禁人妻一区二区| 热re99久久国产66热| 热99久久久久精品小说推荐| 色综合欧美亚洲国产小说| 一进一出抽搐gif免费好疼 | 黄片播放在线免费| 国产不卡一卡二| 久久香蕉激情| 亚洲五月色婷婷综合| 亚洲人成77777在线视频| 男女高潮啪啪啪动态图| 嫩草影视91久久| 18禁裸乳无遮挡免费网站照片 | 99在线人妻在线中文字幕 | 天堂动漫精品| 黄色丝袜av网址大全| 女人高潮潮喷娇喘18禁视频| 波多野结衣一区麻豆| 一边摸一边抽搐一进一出视频| 久久久水蜜桃国产精品网| 国产单亲对白刺激| 午夜免费成人在线视频| 国产欧美日韩综合在线一区二区| 国产视频一区二区在线看| 极品教师在线免费播放| 精品国产乱子伦一区二区三区| 亚洲精品在线美女| 久久久水蜜桃国产精品网| 一级片'在线观看视频| 免费日韩欧美在线观看| 成年人免费黄色播放视频| 亚洲国产中文字幕在线视频| 欧美精品高潮呻吟av久久| 久99久视频精品免费| 黑人巨大精品欧美一区二区mp4| 丝袜美腿诱惑在线| 在线观看免费高清a一片| 亚洲精品一二三| 欧美成人免费av一区二区三区 | 午夜精品在线福利| 国产亚洲av高清不卡| 午夜亚洲福利在线播放| 久久精品国产综合久久久| 欧美老熟妇乱子伦牲交| 国产精品久久电影中文字幕 | 中文字幕人妻熟女乱码| 中文亚洲av片在线观看爽 | 国产一区二区三区视频了| 交换朋友夫妻互换小说| 老司机亚洲免费影院| 亚洲精品粉嫩美女一区| 搡老乐熟女国产| 亚洲欧美色中文字幕在线| 黄色片一级片一级黄色片| 日韩中文字幕欧美一区二区| 亚洲精华国产精华精| 好看av亚洲va欧美ⅴa在| 啦啦啦 在线观看视频| 免费av中文字幕在线| 天天躁狠狠躁夜夜躁狠狠躁| 日韩一卡2卡3卡4卡2021年| 国产无遮挡羞羞视频在线观看| 亚洲国产欧美一区二区综合| 欧美精品啪啪一区二区三区| 精品免费久久久久久久清纯 | 99国产精品一区二区蜜桃av | 中国美女看黄片| 757午夜福利合集在线观看| 一个人免费在线观看的高清视频| 一区二区三区国产精品乱码| 亚洲熟妇熟女久久| 国产成人av激情在线播放| 精品人妻在线不人妻| 久久久久视频综合| 色在线成人网| 欧美日韩亚洲国产一区二区在线观看 | 国产免费现黄频在线看| 丁香六月欧美| 国产精品九九99| 制服人妻中文乱码| 另类亚洲欧美激情| 亚洲精品国产色婷婷电影| 国产不卡一卡二| 18禁裸乳无遮挡动漫免费视频| 欧美+亚洲+日韩+国产| 久久九九热精品免费| 又紧又爽又黄一区二区| 这个男人来自地球电影免费观看| 又紧又爽又黄一区二区| 国产黄色免费在线视频| 亚洲精品成人av观看孕妇| 满18在线观看网站| 一边摸一边抽搐一进一出视频| 日韩欧美免费精品| 国产亚洲精品一区二区www | 精品福利观看| a级毛片在线看网站| 超碰成人久久| 免费女性裸体啪啪无遮挡网站| 少妇猛男粗大的猛烈进出视频| 精品人妻熟女毛片av久久网站| 久久国产精品人妻蜜桃| 欧美老熟妇乱子伦牲交| 亚洲情色 制服丝袜| 亚洲视频免费观看视频| 91大片在线观看| av视频免费观看在线观看| 少妇的丰满在线观看| 精品亚洲成国产av| 可以免费在线观看a视频的电影网站| 女人高潮潮喷娇喘18禁视频| 在线av久久热| 免费黄频网站在线观看国产| 成年女人毛片免费观看观看9 | 国产极品粉嫩免费观看在线| 免费在线观看影片大全网站| 在线观看日韩欧美| 国产精品综合久久久久久久免费 | 欧美日本中文国产一区发布| 亚洲精品av麻豆狂野| 99热只有精品国产| 怎么达到女性高潮| 国产三级黄色录像| 黄色毛片三级朝国网站| 国产精品一区二区在线观看99| 真人做人爱边吃奶动态| 一个人免费在线观看的高清视频| 咕卡用的链子| 极品教师在线免费播放| 老司机深夜福利视频在线观看| 国产有黄有色有爽视频| 国产成+人综合+亚洲专区| 俄罗斯特黄特色一大片| 又黄又粗又硬又大视频| 亚洲aⅴ乱码一区二区在线播放 | 80岁老熟妇乱子伦牲交| 欧美大码av| 国产欧美日韩一区二区三| 午夜久久久在线观看| 在线观看免费视频网站a站| 777米奇影视久久| 国产xxxxx性猛交| 精品午夜福利视频在线观看一区| 午夜福利乱码中文字幕| av天堂久久9| 日韩欧美一区视频在线观看| 99久久综合精品五月天人人| 日韩 欧美 亚洲 中文字幕| 少妇的丰满在线观看| 国产精品久久久av美女十八| 免费久久久久久久精品成人欧美视频| 他把我摸到了高潮在线观看| 久久久久久久久免费视频了| 黄色视频不卡| 久久精品国产综合久久久| 成在线人永久免费视频| av不卡在线播放| 亚洲一码二码三码区别大吗| 成年版毛片免费区| 麻豆av在线久日| 国产亚洲av高清不卡| 99riav亚洲国产免费| 麻豆乱淫一区二区| 免费观看a级毛片全部| 老汉色∧v一级毛片| 中文字幕人妻熟女乱码| 日本五十路高清| 国产精品亚洲av一区麻豆| 午夜影院日韩av| 中出人妻视频一区二区| 80岁老熟妇乱子伦牲交| 精品久久久久久久久久免费视频 | 久久国产精品大桥未久av| 韩国av一区二区三区四区| 中文字幕另类日韩欧美亚洲嫩草| 欧美国产精品va在线观看不卡| 亚洲欧美一区二区三区黑人| 黄色成人免费大全| 最新的欧美精品一区二区| 亚洲欧美日韩另类电影网站| 久久久久国内视频| 12—13女人毛片做爰片一| 18禁黄网站禁片午夜丰满| 亚洲精品在线观看二区| 色综合欧美亚洲国产小说| 国产高清激情床上av| 搡老乐熟女国产| 日本精品一区二区三区蜜桃| 亚洲成人手机| 男女免费视频国产| 国产亚洲av高清不卡| 香蕉国产在线看| 天堂俺去俺来也www色官网| 免费观看a级毛片全部| 亚洲 国产 在线| 免费在线观看亚洲国产| 精品一区二区三卡| 欧美国产精品一级二级三级| 亚洲熟妇中文字幕五十中出 | 国产精品影院久久| 欧美成人免费av一区二区三区 | 中文字幕色久视频| 亚洲成国产人片在线观看| 国产野战对白在线观看| 法律面前人人平等表现在哪些方面| 一进一出好大好爽视频| 黄色怎么调成土黄色| 在线天堂中文资源库| 一区二区三区国产精品乱码| 亚洲欧美一区二区三区黑人| 国产精品亚洲一级av第二区| 久99久视频精品免费| 新久久久久国产一级毛片| 欧美黄色淫秽网站| 人妻 亚洲 视频| 精品久久久久久久久久免费视频 | 男女免费视频国产| 亚洲熟女毛片儿| 中文字幕人妻熟女乱码| 天天躁狠狠躁夜夜躁狠狠躁| 国产片内射在线| 国产欧美日韩精品亚洲av| 午夜日韩欧美国产| 人成视频在线观看免费观看| 久久久精品区二区三区| 精品国产一区二区久久| 国产xxxxx性猛交| 伦理电影免费视频| av有码第一页| www日本在线高清视频| 丝袜人妻中文字幕| 亚洲专区中文字幕在线| 99久久综合精品五月天人人| 日韩欧美三级三区| 免费一级毛片在线播放高清视频 | 欧美日本中文国产一区发布| 国精品久久久久久国模美| 人人妻人人澡人人爽人人夜夜| 国产在视频线精品| 久久久久国内视频| 日韩制服丝袜自拍偷拍| 欧美一级毛片孕妇| 久久精品国产清高在天天线| 中文字幕高清在线视频| 最近最新中文字幕大全免费视频| 国产视频一区二区在线看| 99精品久久久久人妻精品| 免费黄频网站在线观看国产| 久久ye,这里只有精品| 国产一区二区激情短视频| 老司机深夜福利视频在线观看| 亚洲精品自拍成人| 久久久久久人人人人人| 91麻豆av在线| 亚洲久久久国产精品| 久久精品国产99精品国产亚洲性色 | 欧美精品av麻豆av| 久久精品人人爽人人爽视色| 9热在线视频观看99| 亚洲色图 男人天堂 中文字幕| 男女午夜视频在线观看| 国产精品久久电影中文字幕 | 久久久久精品国产欧美久久久| 日韩三级视频一区二区三区| 高清欧美精品videossex| 最近最新免费中文字幕在线| 欧美 日韩 精品 国产| 亚洲成国产人片在线观看| 亚洲精品一二三| 丝袜美足系列| 国产精品久久久久久人妻精品电影| 国产精品久久久久久精品古装| 日本vs欧美在线观看视频| 人妻 亚洲 视频| 久久 成人 亚洲| 亚洲一码二码三码区别大吗| 欧美亚洲日本最大视频资源| 久久草成人影院| 久热这里只有精品99| 日日爽夜夜爽网站| 天天操日日干夜夜撸| 老司机靠b影院| 可以免费在线观看a视频的电影网站| 亚洲成人国产一区在线观看| 咕卡用的链子| 黄片大片在线免费观看| 国产野战对白在线观看| 成人18禁高潮啪啪吃奶动态图| 日韩熟女老妇一区二区性免费视频| 欧美性长视频在线观看| 91国产中文字幕| 亚洲美女黄片视频| 两个人看的免费小视频| 国产精品自产拍在线观看55亚洲 | 99久久国产精品久久久| 人妻 亚洲 视频| 99国产精品99久久久久| 久久亚洲真实| 欧美+亚洲+日韩+国产| 午夜成年电影在线免费观看| 久久久国产精品麻豆| 九色亚洲精品在线播放| 国产精品秋霞免费鲁丝片| 久久久精品免费免费高清| 成人精品一区二区免费| 夜夜爽天天搞| 老熟妇仑乱视频hdxx| 熟女少妇亚洲综合色aaa.| 成年人免费黄色播放视频| 成年人午夜在线观看视频| 午夜精品在线福利| 啦啦啦免费观看视频1| 久久青草综合色| 69av精品久久久久久| 三上悠亚av全集在线观看| 亚洲综合色网址| 亚洲精品粉嫩美女一区| 国产精品99久久99久久久不卡| 午夜老司机福利片| 国产欧美日韩一区二区三| 搡老熟女国产l中国老女人| 一区二区三区精品91| 波多野结衣av一区二区av| 丰满迷人的少妇在线观看| 欧美激情高清一区二区三区| 亚洲成人免费电影在线观看| 两人在一起打扑克的视频| 在线免费观看的www视频| 亚洲自偷自拍图片 自拍| 欧美精品一区二区免费开放| 国产一区在线观看成人免费| 免费一级毛片在线播放高清视频 | 久久人妻熟女aⅴ| 丝袜在线中文字幕| 校园春色视频在线观看| 国产午夜精品久久久久久| 国产精品一区二区精品视频观看| 极品人妻少妇av视频| 黄片小视频在线播放| 免费在线观看视频国产中文字幕亚洲| 男人操女人黄网站| 美女午夜性视频免费| 真人做人爱边吃奶动态| 这个男人来自地球电影免费观看| 亚洲成人免费av在线播放| 最新美女视频免费是黄的| 男女之事视频高清在线观看| 中文字幕制服av| 欧美国产精品一级二级三级| 99精国产麻豆久久婷婷| 国产99白浆流出| 久久久国产一区二区| 久久香蕉国产精品| 国产亚洲av高清不卡| 亚洲一区高清亚洲精品| 日本精品一区二区三区蜜桃| 亚洲午夜精品一区,二区,三区| 在线观看免费午夜福利视频| 日韩欧美一区二区三区在线观看 | 女人久久www免费人成看片| 又紧又爽又黄一区二区| 国产成人欧美在线观看 | 国产成+人综合+亚洲专区| 一区在线观看完整版| 人妻 亚洲 视频| 国产在视频线精品| 99精品久久久久人妻精品| 在线永久观看黄色视频| 搡老熟女国产l中国老女人| 一区在线观看完整版| 午夜福利一区二区在线看| 成人永久免费在线观看视频| 制服诱惑二区| 80岁老熟妇乱子伦牲交| 中文字幕精品免费在线观看视频| av福利片在线| 国产91精品成人一区二区三区| 成人永久免费在线观看视频| 欧美日韩一级在线毛片| 一进一出好大好爽视频| 欧美性长视频在线观看| 国产成人一区二区三区免费视频网站|