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

    2011—2015年新疆維吾爾自治區(qū)涂陽(yáng)肺結(jié)核時(shí)空格局分析

    2020-06-09 07:22:40尹哲賀湘焱李奇鳳劉亞潔張燕李德洋加依娜拉茲別克古麗娜扎爾艾克拜爾曹明芹
    中國(guó)防癆雜志 2020年6期
    關(guān)鍵詞:新疆模型

    尹哲 賀湘焱 李奇鳳 劉亞潔 張燕 李德洋 加依娜·拉茲別克 古麗娜扎爾·艾克拜爾 曹明芹

    結(jié)核病傳播具有空間聚集性分布的特點(diǎn),在經(jīng)濟(jì)“邊緣化和脆弱地區(qū)”高發(fā),根據(jù)全國(guó)第五次結(jié)核病流行病學(xué)抽樣調(diào)查,西部地區(qū)涂陽(yáng)肺結(jié)核患病率是中部地區(qū)的1.7倍,是東部地區(qū)的2.4倍[1]。新疆維吾爾自治區(qū)(簡(jiǎn)稱“新疆”)地處我國(guó)西北部,經(jīng)濟(jì)發(fā)展相對(duì)落后,加之衛(wèi)生資源稀缺等因素造成基礎(chǔ)疫情嚴(yán)重,防治困難,年平均涂陽(yáng)肺結(jié)核登記率為39.31/10萬(wàn)[2],是我國(guó)結(jié)核病發(fā)病比較嚴(yán)重的地區(qū)之一[3]。

    伴隨著空間流行病學(xué)理論、技術(shù)的快速發(fā)展,地理信息系統(tǒng)(geographical information system,GIS)和空間分析軟件應(yīng)用日益廣泛,利用地理統(tǒng)計(jì)學(xué)的知識(shí)從空間角度揭示結(jié)核病的分布與傳播規(guī)律,已經(jīng)成為廣大科研工作者熱衷研究的方向??死锝鸩逯捣ㄍㄟ^(guò)變異函數(shù)來(lái)研究具有空間相關(guān)性和依賴性的地理分布現(xiàn)象,主要應(yīng)用于采礦業(yè)、農(nóng)業(yè)、地質(zhì)等領(lǐng)域,隨著地質(zhì)統(tǒng)計(jì)學(xué)理論的不斷成熟,應(yīng)用范圍持續(xù)擴(kuò)展[4-5]。

    本研究在2011—2015年新疆涂陽(yáng)肺結(jié)核發(fā)病風(fēng)險(xiǎn)時(shí)空分布探索的基礎(chǔ)上,根據(jù)克里金插值方法最優(yōu)無(wú)偏的特點(diǎn),采用普通克里金法和經(jīng)驗(yàn)貝葉斯克里金法建立模型,探索涂陽(yáng)肺結(jié)核發(fā)病風(fēng)險(xiǎn)時(shí)空分布的相關(guān)性和變異性,評(píng)估新疆不同地區(qū)肺結(jié)核疫情。

    資料和方法

    一、資料來(lái)源

    2011—2015年新疆肺結(jié)核發(fā)病報(bào)告數(shù)據(jù)來(lái)自《中國(guó)疾病預(yù)防控制信息系統(tǒng)》的子系統(tǒng)《傳染病報(bào)告信息管理系統(tǒng)》。2011—2015年新疆涂陽(yáng)肺結(jié)核報(bào)告發(fā)病例數(shù)分別為11 813、14 059、11 149、10 483、10 196例,5年間共計(jì)57 700例。各年份區(qū)(縣)人口數(shù)據(jù)來(lái)自新疆統(tǒng)計(jì)年鑒。省級(jí)行政區(qū)域電子地圖(1∶10萬(wàn))從國(guó)家基礎(chǔ)地理信息系統(tǒng)獲取,各區(qū)(縣)的經(jīng)緯度數(shù)據(jù)由Google earth獲取,以2011年為標(biāo)準(zhǔn)。

    二、分析方法

    1. 涂陽(yáng)肺結(jié)核標(biāo)準(zhǔn)化發(fā)病比(standardized morbidity ratio,SMR)估算:2011—2015年新疆結(jié)核病發(fā)病數(shù)和發(fā)病率數(shù)據(jù)資料經(jīng)整理、計(jì)算,構(gòu)建各區(qū)(縣)涂陽(yáng)肺結(jié)核SMR數(shù)據(jù)庫(kù)。以各區(qū)(縣)年度報(bào)告的全疆人口構(gòu)成和涂陽(yáng)肺結(jié)核發(fā)病率作為標(biāo)準(zhǔn),計(jì)算各個(gè)區(qū)(縣)涂陽(yáng)肺結(jié)核SMR,公式如下:

    公式(1)[6]

    其中,yij表示i區(qū)(縣)(1≤i≤98)在第j年(2011≤j≤2015)的涂陽(yáng)肺結(jié)核報(bào)告例數(shù),Eij表示i區(qū)(縣)在第j年的期望報(bào)告例數(shù),由當(dāng)?shù)厝丝跀?shù)乘以全疆第j年的涂陽(yáng)肺結(jié)核報(bào)告率得到。

    2. 空間自相關(guān)分析:采用全局莫蘭指數(shù)(MoranI)值反映空間自相關(guān)關(guān)系,公式如下:

    公式(2)[7]

    其中,zi是第i個(gè)區(qū)(縣)肺結(jié)核SMR與其平均值的偏差,wij是i和j之間的空間連接矩陣,n是要素總數(shù),即觀測(cè)區(qū)域數(shù)等于98,S0是98個(gè)區(qū)(縣)空間權(quán)重的聚合。MoranI取值范圍為-1~1,MoranI值越趨近1,代表肺結(jié)核SMR呈聚集分布;MoranI值越趨近-1,代表肺結(jié)核發(fā)病率分布具有負(fù)相關(guān)性;MoranI值為零,表明整體呈隨機(jī)分布。假設(shè)檢驗(yàn)使用Z-score得分,檢驗(yàn)水準(zhǔn)α=0.05;Z>1.96或<-1.96時(shí),α<0.05。

    3. 克里金法:普通克里金法可對(duì)周?chē)臏y(cè)量值進(jìn)行加權(quán)以得出未測(cè)量位置的一個(gè)線性估計(jì)預(yù)測(cè)系統(tǒng),適用于區(qū)域化變量存在空間相關(guān)性的情況,能夠給出估計(jì)誤差,充分考慮空間變量相關(guān)性,有效彌補(bǔ)數(shù)據(jù)集存在的聚類影響,插值精度高。公式如下:

    公式(3)[8]

    其中,Z(Si) 是第i個(gè)位置處的實(shí)際涂陽(yáng)肺結(jié)核SMR,λi是第i個(gè)位置處的未知權(quán)重,S0是預(yù)估點(diǎn),N是測(cè)量點(diǎn)數(shù)。在普通克里金法中,權(quán)重λi取決于涂陽(yáng)肺結(jié)核SMR實(shí)際點(diǎn)與待測(cè)點(diǎn)位置的距離,以及待測(cè)點(diǎn)周?chē)膶?shí)際涂陽(yáng)肺結(jié)核SMR之間空間關(guān)系的擬合模型。

    經(jīng)驗(yàn)貝葉斯克里金法是地質(zhì)統(tǒng)計(jì)學(xué)建模技術(shù)之一。它與其他克里金模型不同之處在于其使用內(nèi)在隨機(jī)函數(shù)進(jìn)行空間插值。除了考慮半變異函數(shù)模型基本參數(shù)的不確定性之外,可以快速擬合許多理論半變異函數(shù)模型[9],從中選擇合適的理論模型進(jìn)行最優(yōu)擬合。根據(jù){z(xi);i=1,2,…,N}已知觀測(cè)點(diǎn)實(shí)際涂陽(yáng)肺結(jié)核SMR,進(jìn)一步求得未知測(cè)量點(diǎn)x0(x0∈D)處的貝葉斯克里金估計(jì)值,即:

    公式(4)[10]

    采用均方根誤差(root mean squared error,RMSE)來(lái)評(píng)價(jià)兩種模型的擬合效果,表示實(shí)際值與模型估計(jì)值之間的平均距離,RMSE值越小,誤差越小,模型效果越好。

    三、統(tǒng)計(jì)學(xué)處理

    采用Excel 2010軟件對(duì)2011—2015年新疆肺結(jié)核發(fā)病報(bào)告數(shù)據(jù)進(jìn)行整理。通過(guò)SPSS 25.0軟件進(jìn)行統(tǒng)計(jì)學(xué)分析。使用ArcGIS 10.2軟件實(shí)現(xiàn)空間分布圖、全局相關(guān)性分析、克里金插值分析及其結(jié)果的可視化。采用均方根誤差評(píng)價(jià)兩種模型的擬合效果。檢驗(yàn)水準(zhǔn)均為α=0.05。

    結(jié) 果

    一、2011—2015年新疆涂陽(yáng)肺結(jié)核報(bào)告SMR時(shí)空分布格局

    2011—2015年新疆涂陽(yáng)肺結(jié)核SMR較高的地區(qū)主要分布于南疆的喀什、和田、阿克蘇3個(gè)地區(qū),具體表現(xiàn)在英吉沙縣、澤普縣、于田縣、民豐縣、洛浦縣等13個(gè)區(qū)(縣),見(jiàn)表1。從空間上看,涂陽(yáng)肺結(jié)核SMR水平存在從南向北逐漸降低的趨勢(shì),南疆SMR水平一直居高不下,北疆、東疆的SMR處于較低水平。5年間,各地區(qū)涂陽(yáng)肺結(jié)核SMR均呈現(xiàn)小幅波動(dòng)、整體下降的趨勢(shì)。南疆涂陽(yáng)肺結(jié)核SMR水平較為平穩(wěn),但遠(yuǎn)高于北疆和東疆,在高位上下波動(dòng)。北疆和東疆始終保持在較低水平,以SMR水平降低為主要趨勢(shì),2015年呈現(xiàn)較為明顯的升高趨勢(shì),見(jiàn)圖1~5。

    二、涂陽(yáng)肺結(jié)核SMR全局自相關(guān)分析

    利用全局MoranI指數(shù)對(duì)2011—2015年新疆區(qū)(縣)涂陽(yáng)肺結(jié)核疫情進(jìn)行全局空間自相關(guān)分析。結(jié)果顯示,各年度全疆涂陽(yáng)肺結(jié)核SMR分布存在空間正相關(guān)關(guān)系,且差異有統(tǒng)計(jì)學(xué)意義(MoranI值均>0,P值均<0.001)。其中,整體表現(xiàn)出的聚集性在2013年最高(MoranI=0.376),在2014年最低(MoranI=0.248),詳見(jiàn)表2。

    三、克里金插值模型構(gòu)建涂陽(yáng)肺結(jié)核SMR分布格局

    對(duì)2011—2015年新疆涂陽(yáng)肺結(jié)核SMR數(shù)據(jù)庫(kù)進(jìn)行空間探索分析,變異云圖提示數(shù)據(jù)存在空間相關(guān)性,可以進(jìn)行插值分析。分別建立普通克里金和經(jīng)驗(yàn)貝葉斯克里金模型,生成新疆涂陽(yáng)肺結(jié)核發(fā)病SMR的克里金估計(jì)圖。交叉驗(yàn)證比較,2011—2015年普通克里金模型RMSE范圍在0.379~0.522,經(jīng)驗(yàn)貝葉斯克里金模型RMSE范圍在0.382~0.484,兩種模型擬合度較高。與新疆涂陽(yáng)肺結(jié)核SMR分布圖比較,符合實(shí)際分布規(guī)律,且更為直觀地表現(xiàn)出SMR在空間地理維度上的變化趨勢(shì),提示涂陽(yáng)肺結(jié)核SMR在新疆南部較高,北部較低,中部介于二者之間,南北疆之間存在一個(gè)較為明顯的分隔帶,走向基本與天山山脈走向一致,見(jiàn)圖6~15。

    表1 2011—2015年新疆肺結(jié)核SMR的分布

    表2 2011—2015年新疆涂陽(yáng)肺結(jié)核SMR全局

    討 論

    一、新疆涂陽(yáng)肺結(jié)核SMR時(shí)空格局的探討

    本研究對(duì)2011—2015年新疆98個(gè)區(qū)(縣)涂陽(yáng)肺結(jié)核SMR進(jìn)行分析,以區(qū)(縣)級(jí)為尺度,為新疆尋找更有針對(duì)性的肺結(jié)核精準(zhǔn)防控地區(qū)。全局自相關(guān)分析結(jié)果提示,涂陽(yáng)肺結(jié)核SMR分布非隨機(jī),具有較強(qiáng)的地區(qū)變化規(guī)律。由北向南、由東向西SMR逐漸增高,呈階梯式分布。涂陽(yáng)肺結(jié)核發(fā)病高SMR地區(qū)主要分布在南疆的喀什、和田、阿克蘇3個(gè)地區(qū),具體表現(xiàn)在英吉沙縣、澤普縣、于田縣等10余個(gè)區(qū)(縣)。5年間,新疆涂陽(yáng)肺結(jié)核SMR平均值在0.9~2.0之間上下波動(dòng),2011—2014年在波動(dòng)中持續(xù)下降,2015年有小幅上漲。從宏觀上看,新疆SMR整體呈現(xiàn)一個(gè)好轉(zhuǎn)的趨勢(shì),高SMR地區(qū)范圍逐漸縮小,符合新疆結(jié)核病發(fā)展分布趨勢(shì)。但新疆肺結(jié)核防控工作任務(wù)仍然十分艱巨,需要進(jìn)一步強(qiáng)化和完善各項(xiàng)防治政策,落實(shí)各項(xiàng)防治措施。

    圖1~5 2011—2015年新疆涂陽(yáng)肺結(jié)核SMR時(shí)空分布

    圖6~10 2011—2015年新疆涂陽(yáng)肺結(jié)核SMR普通克里金估計(jì)值

    圖11~15 2011—2015年新疆涂陽(yáng)肺結(jié)核SMR經(jīng)驗(yàn)貝葉斯克里金估計(jì)值

    二、克里金插值在估計(jì)疾病風(fēng)險(xiǎn)中的應(yīng)用價(jià)值

    克里金插值法可對(duì)所研究的對(duì)象提供一種最佳線性無(wú)偏估計(jì),利用半方差函數(shù)圖將統(tǒng)計(jì)相關(guān)性的強(qiáng)度作為距離函數(shù)來(lái)測(cè)量。本研究用涂陽(yáng)肺結(jié)核SMR數(shù)據(jù)庫(kù)分別建立普通克里金和經(jīng)驗(yàn)貝葉斯克里金模型。建模發(fā)現(xiàn)SMR分布與實(shí)際分布基本一致,誤差服從正態(tài)分布,估計(jì)圖、誤差圖均提示兩個(gè)模型擬合度較為理想。進(jìn)行交叉驗(yàn)證比較,采用均方根誤差指標(biāo)評(píng)價(jià)兩模型擬合度,從估計(jì)值和真實(shí)值的離散程度上來(lái)說(shuō),經(jīng)驗(yàn)貝葉斯克里金法(RMSE范圍在0.382~0.484)略優(yōu)于普通克里金法(RMSE范圍在0.379~0.522),擬合度更加符合實(shí)際分布。

    克里金插值估計(jì)圖可直觀展現(xiàn)新疆涂陽(yáng)肺結(jié)核SMR階梯式變化,兩模型均提示在北疆與南疆之間存在一條與天山山脈走向較為相似的分隔帶。天山山脈將新疆分為南北兩大部分,南、北疆在海拔、氣候、經(jīng)濟(jì)、衛(wèi)生等方面差異明顯。南疆地區(qū)全年降水量一般不足100 mm,遠(yuǎn)低于北疆,且沙漠化現(xiàn)象嚴(yán)重[11]。西部大開(kāi)發(fā)以來(lái),隨著新亞歐大陸橋的開(kāi)通,北疆經(jīng)濟(jì)迅猛發(fā)展,南北疆經(jīng)濟(jì)差距不斷擴(kuò)大,“十二五”期間,天山北坡經(jīng)濟(jì)帶建立,作為我國(guó)14個(gè)重點(diǎn)開(kāi)發(fā)區(qū)之一,一直是新疆的重點(diǎn)開(kāi)發(fā)區(qū),是新疆經(jīng)濟(jì)實(shí)力最強(qiáng)、城鎮(zhèn)最集中分布的區(qū)域,這也進(jìn)一步擴(kuò)大了地區(qū)經(jīng)濟(jì)差異[12]。研究表明,環(huán)境因素、經(jīng)濟(jì)因素對(duì)結(jié)核病的發(fā)生發(fā)展有著實(shí)質(zhì)性影響,國(guó)內(nèi)外肺結(jié)核疫情報(bào)告中多有報(bào)道,經(jīng)濟(jì)水平低下地區(qū)、自然條件惡劣地區(qū)的結(jié)核病的發(fā)病率較高[13]。

    克里金插值作為一種局部估計(jì)方法,對(duì)估計(jì)值的整體空間相關(guān)性考慮不夠,其能保證數(shù)據(jù)的估計(jì)局部最優(yōu),卻不能保證數(shù)據(jù)的總體最優(yōu),因?yàn)榭死锝鸸乐档姆讲畋仍紨?shù)據(jù)的方差要小。因此,當(dāng)數(shù)據(jù)點(diǎn)較少且分布不均時(shí)可能會(huì)出現(xiàn)較大的估計(jì)誤差。此外克里金插值法作為內(nèi)插方法,為減小估計(jì)方差而對(duì)真實(shí)觀測(cè)數(shù)據(jù)的離散性進(jìn)行了平滑處理,雖然可以得到由于平滑而更美觀的等值線圖,但一些有意義的異常帶也可能被平滑作用而平滑掉[14]。在本研究的數(shù)據(jù)上也有體現(xiàn),具體表現(xiàn)為經(jīng)過(guò)克里金處理后,新疆涂陽(yáng)肺結(jié)核SMR的估計(jì)值跟實(shí)際值比較起來(lái)有較為明顯的下降,波動(dòng)范圍縮小。但克里金法通過(guò)有限的離散點(diǎn)數(shù)據(jù),研究了整個(gè)區(qū)域的分布特征,直觀地反映了變量的空間結(jié)構(gòu)性;并初步提供了未來(lái)新疆地區(qū)結(jié)核病防范的重點(diǎn)區(qū)域,為及時(shí)調(diào)整或完善防控措施提供線索和依據(jù),對(duì)未來(lái)新疆結(jié)核病的防治有著重要意義。

    猜你喜歡
    新疆模型
    一半模型
    走進(jìn)新疆
    在新疆(四首)
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計(jì)的漸近分布
    3D打印中的模型分割與打包
    FLUKA幾何模型到CAD幾何模型轉(zhuǎn)換方法初步研究
    新疆多怪
    絲綢之路(2014年9期)2015-01-22 04:24:46
    新疆對(duì)外開(kāi)放山峰
    新疆對(duì)外開(kāi)放一類口岸
    美女扒开内裤让男人捅视频| 欧美绝顶高潮抽搐喷水| 久久久久国内视频| 国产成年人精品一区二区| 天天躁狠狠躁夜夜躁狠狠躁| 国产精品电影一区二区三区| 成人亚洲精品av一区二区| 国产熟女午夜一区二区三区| 亚洲七黄色美女视频| 久久天堂一区二区三区四区| 一区二区三区高清视频在线| 国产免费男女视频| 天堂影院成人在线观看| 午夜精品久久久久久毛片777| 听说在线观看完整版免费高清| 国产真人三级小视频在线观看| 国产乱人伦免费视频| 精品国内亚洲2022精品成人| 亚洲人成伊人成综合网2020| 啪啪无遮挡十八禁网站| 波多野结衣高清作品| 久久久久亚洲av毛片大全| 日日干狠狠操夜夜爽| 国产成人系列免费观看| 亚洲专区国产一区二区| 日日爽夜夜爽网站| 欧美黄色淫秽网站| 成人av在线播放网站| а√天堂www在线а√下载| 国产高清视频在线观看网站| 美女午夜性视频免费| 国产精品久久久久久久电影 | 悠悠久久av| 欧美绝顶高潮抽搐喷水| 亚洲欧美日韩高清在线视频| 国产69精品久久久久777片 | 我要搜黄色片| 黄片小视频在线播放| 99国产极品粉嫩在线观看| 亚洲成人免费电影在线观看| 欧美 亚洲 国产 日韩一| 伦理电影免费视频| 国产熟女午夜一区二区三区| 1024视频免费在线观看| 舔av片在线| 一进一出好大好爽视频| 神马国产精品三级电影在线观看 | 日韩三级视频一区二区三区| 久久亚洲精品不卡| 国产不卡一卡二| 黄色 视频免费看| 小说图片视频综合网站| 国产黄片美女视频| 曰老女人黄片| netflix在线观看网站| 精品第一国产精品| 亚洲avbb在线观看| 成人国语在线视频| xxxwww97欧美| 99riav亚洲国产免费| 日韩欧美免费精品| av有码第一页| 国产成年人精品一区二区| 亚洲精品在线观看二区| 亚洲精品在线美女| 老司机深夜福利视频在线观看| 国产探花在线观看一区二区| 日本三级黄在线观看| 久久久久国产一级毛片高清牌| 日韩有码中文字幕| 俺也久久电影网| 亚洲精华国产精华精| 亚洲一码二码三码区别大吗| 一级片免费观看大全| 99国产精品一区二区蜜桃av| 日本a在线网址| 香蕉丝袜av| 亚洲狠狠婷婷综合久久图片| 黄色视频,在线免费观看| 成在线人永久免费视频| 亚洲男人天堂网一区| 亚洲av美国av| 1024香蕉在线观看| 两个人免费观看高清视频| 国产在线观看jvid| 中文亚洲av片在线观看爽| 香蕉久久夜色| 日韩国内少妇激情av| 午夜日韩欧美国产| 男女那种视频在线观看| 国产探花在线观看一区二区| 免费一级毛片在线播放高清视频| 久久人妻av系列| 国产精品亚洲av一区麻豆| 在线看三级毛片| 久久久久久亚洲精品国产蜜桃av| 最近在线观看免费完整版| 中文字幕熟女人妻在线| 视频区欧美日本亚洲| 国产精品一区二区三区四区久久| 黄色女人牲交| 欧美黑人精品巨大| 一a级毛片在线观看| 久久久国产欧美日韩av| 亚洲 国产 在线| 久久久久久人人人人人| 在线a可以看的网站| 亚洲av五月六月丁香网| 欧美绝顶高潮抽搐喷水| 男男h啪啪无遮挡| 成人欧美大片| 老司机深夜福利视频在线观看| 中文在线观看免费www的网站 | 亚洲精品一卡2卡三卡4卡5卡| 国产伦一二天堂av在线观看| 免费一级毛片在线播放高清视频| 久久人妻av系列| 国产精品 欧美亚洲| 最新在线观看一区二区三区| 中文字幕精品亚洲无线码一区| 国产成人一区二区三区免费视频网站| 一级作爱视频免费观看| 好男人在线观看高清免费视频| 久久久久九九精品影院| 久久久久国内视频| 午夜福利在线观看吧| 老鸭窝网址在线观看| 国产亚洲精品久久久久5区| 久久中文字幕一级| 叶爱在线成人免费视频播放| 国产精品爽爽va在线观看网站| 变态另类丝袜制服| 精品久久久久久久末码| 毛片女人毛片| 狂野欧美白嫩少妇大欣赏| 最近在线观看免费完整版| 桃红色精品国产亚洲av| 欧美成人一区二区免费高清观看 | 久久久久久久久久黄片| 五月玫瑰六月丁香| 亚洲熟妇熟女久久| 国产欧美日韩精品亚洲av| 1024香蕉在线观看| 99热只有精品国产| 国产成年人精品一区二区| 99精品欧美一区二区三区四区| 99热只有精品国产| 深夜精品福利| 亚洲熟妇熟女久久| 亚洲中文av在线| 久久久国产精品麻豆| 一个人免费在线观看的高清视频| 亚洲av电影在线进入| 成在线人永久免费视频| 嫩草影院精品99| 久久久精品欧美日韩精品| 成在线人永久免费视频| 俄罗斯特黄特色一大片| 91av网站免费观看| 国产黄色小视频在线观看| 亚洲av第一区精品v没综合| 叶爱在线成人免费视频播放| 国产精品av久久久久免费| 免费看美女性在线毛片视频| av超薄肉色丝袜交足视频| 久99久视频精品免费| 久久精品国产99精品国产亚洲性色| 国产成人啪精品午夜网站| 日本一本二区三区精品| 波多野结衣巨乳人妻| 99国产精品99久久久久| 国产精品98久久久久久宅男小说| 欧美三级亚洲精品| 好男人电影高清在线观看| 美女午夜性视频免费| 久久中文看片网| 制服丝袜大香蕉在线| 欧美在线一区亚洲| 中国美女看黄片| 国产真实乱freesex| 国产精品久久电影中文字幕| 久久人妻av系列| 日日摸夜夜添夜夜添小说| 欧美 亚洲 国产 日韩一| 精品福利观看| 久久香蕉精品热| 欧美+亚洲+日韩+国产| av国产免费在线观看| 亚洲片人在线观看| 天堂动漫精品| 一个人免费在线观看电影 | 国产欧美日韩精品亚洲av| or卡值多少钱| 18禁美女被吸乳视频| 亚洲精华国产精华精| 在线看三级毛片| 丰满人妻熟妇乱又伦精品不卡| 国产精品一区二区免费欧美| 亚洲午夜理论影院| 精品福利观看| 好男人电影高清在线观看| 亚洲欧美日韩无卡精品| 黄片小视频在线播放| 听说在线观看完整版免费高清| 亚洲精华国产精华精| 欧美人与性动交α欧美精品济南到| 免费在线观看成人毛片| www.熟女人妻精品国产| 国产av麻豆久久久久久久| 一本大道久久a久久精品| 国产精品久久久久久久电影 | 热99re8久久精品国产| 成年免费大片在线观看| 午夜成年电影在线免费观看| 国产成人av激情在线播放| 亚洲av中文字字幕乱码综合| 亚洲国产精品合色在线| 亚洲电影在线观看av| 免费看a级黄色片| 一二三四社区在线视频社区8| 免费观看人在逋| 91国产中文字幕| 老司机午夜福利在线观看视频| 啦啦啦韩国在线观看视频| 少妇熟女aⅴ在线视频| 男人舔女人的私密视频| 午夜免费观看网址| 最近视频中文字幕2019在线8| 日韩高清综合在线| 亚洲国产中文字幕在线视频| 日韩欧美在线乱码| 国产高清视频在线播放一区| 九色成人免费人妻av| 麻豆成人午夜福利视频| 一本精品99久久精品77| 搞女人的毛片| 三级国产精品欧美在线观看 | 村上凉子中文字幕在线| 日韩欧美精品v在线| 十八禁网站免费在线| 51午夜福利影视在线观看| 嫁个100分男人电影在线观看| 欧美色欧美亚洲另类二区| 国内久久婷婷六月综合欲色啪| 日日摸夜夜添夜夜添小说| 国产精品一区二区三区四区免费观看 | 母亲3免费完整高清在线观看| 麻豆成人午夜福利视频| 午夜福利成人在线免费观看| 日韩大尺度精品在线看网址| 欧美乱色亚洲激情| 国产99白浆流出| 啦啦啦免费观看视频1| 长腿黑丝高跟| 亚洲av熟女| 在线观看免费午夜福利视频| 熟女少妇亚洲综合色aaa.| 精品国内亚洲2022精品成人| 少妇熟女aⅴ在线视频| 两个人的视频大全免费| av欧美777| a级毛片a级免费在线| 麻豆av在线久日| 国产69精品久久久久777片 | 中出人妻视频一区二区| 最新美女视频免费是黄的| 女人爽到高潮嗷嗷叫在线视频| 在线观看免费视频日本深夜| 女人高潮潮喷娇喘18禁视频| 日韩欧美三级三区| 黄频高清免费视频| 好男人在线观看高清免费视频| 精品欧美一区二区三区在线| 国产探花在线观看一区二区| 无限看片的www在线观看| 老汉色∧v一级毛片| 少妇人妻一区二区三区视频| 国产欧美日韩精品亚洲av| 欧美成人性av电影在线观看| 日韩有码中文字幕| 国产片内射在线| 成人永久免费在线观看视频| 午夜视频精品福利| 免费av毛片视频| 国产伦人伦偷精品视频| 午夜精品久久久久久毛片777| svipshipincom国产片| 黑人操中国人逼视频| 亚洲精品中文字幕在线视频| 婷婷六月久久综合丁香| 亚洲性夜色夜夜综合| 真人一进一出gif抽搐免费| 美女扒开内裤让男人捅视频| 国产久久久一区二区三区| 天堂影院成人在线观看| 淫妇啪啪啪对白视频| 日日爽夜夜爽网站| 国产亚洲av高清不卡| 亚洲国产精品999在线| 亚洲中文av在线| 香蕉国产在线看| 日本在线视频免费播放| 欧美乱色亚洲激情| 18禁观看日本| 亚洲av美国av| 亚洲国产欧美一区二区综合| 性欧美人与动物交配| 国产av麻豆久久久久久久| 熟女电影av网| avwww免费| 国内毛片毛片毛片毛片毛片| 亚洲成人久久爱视频| 欧美大码av| 又黄又爽又免费观看的视频| 神马国产精品三级电影在线观看 | 精品久久蜜臀av无| 亚洲av片天天在线观看| 悠悠久久av| 日韩欧美三级三区| 国产蜜桃级精品一区二区三区| av有码第一页| 久久午夜亚洲精品久久| 精品久久久久久久久久久久久| 天天一区二区日本电影三级| 男人舔奶头视频| 久久精品综合一区二区三区| 香蕉国产在线看| 黄色丝袜av网址大全| 亚洲午夜精品一区,二区,三区| 美女免费视频网站| 免费看美女性在线毛片视频| 女警被强在线播放| 国内精品久久久久精免费| 久久精品91蜜桃| 欧美午夜高清在线| 国产精品 国内视频| 亚洲欧美精品综合久久99| 国产伦在线观看视频一区| 欧美又色又爽又黄视频| 国产一区二区激情短视频| www.www免费av| svipshipincom国产片| 欧美激情久久久久久爽电影| 18禁黄网站禁片免费观看直播| 我要搜黄色片| 两人在一起打扑克的视频| 亚洲国产欧美网| 久久精品人妻少妇| videosex国产| 国产精品美女特级片免费视频播放器 | 国产精品日韩av在线免费观看| 亚洲精品久久国产高清桃花| svipshipincom国产片| 亚洲国产精品合色在线| 美女免费视频网站| 日日干狠狠操夜夜爽| 国产精品亚洲av一区麻豆| 在线观看舔阴道视频| 99国产极品粉嫩在线观看| 一个人免费在线观看的高清视频| 欧美日韩一级在线毛片| 丁香六月欧美| 久久久久久人人人人人| 国产欧美日韩精品亚洲av| 1024手机看黄色片| 高潮久久久久久久久久久不卡| 999久久久精品免费观看国产| 88av欧美| 欧美人与性动交α欧美精品济南到| 欧美一级a爱片免费观看看 | 最近最新中文字幕大全免费视频| 久久精品国产99精品国产亚洲性色| 日本熟妇午夜| 午夜a级毛片| 亚洲精华国产精华精| 国产片内射在线| 一区二区三区激情视频| 日本精品一区二区三区蜜桃| 国产av又大| 一级毛片女人18水好多| 欧美一区二区精品小视频在线| 视频区欧美日本亚洲| av片东京热男人的天堂| 不卡av一区二区三区| avwww免费| 18美女黄网站色大片免费观看| 制服人妻中文乱码| 少妇粗大呻吟视频| 露出奶头的视频| 久久精品国产99精品国产亚洲性色| 超碰成人久久| 一区福利在线观看| 国产成人影院久久av| 淫妇啪啪啪对白视频| 91字幕亚洲| 91av网站免费观看| 一卡2卡三卡四卡精品乱码亚洲| 色老头精品视频在线观看| 他把我摸到了高潮在线观看| 亚洲五月婷婷丁香| 日日夜夜操网爽| 日本五十路高清| 制服丝袜大香蕉在线| 成人三级做爰电影| 小说图片视频综合网站| 国产成年人精品一区二区| 欧美日韩福利视频一区二区| 亚洲av成人一区二区三| 亚洲精品中文字幕在线视频| avwww免费| 母亲3免费完整高清在线观看| 亚洲人成电影免费在线| 在线国产一区二区在线| 搡老岳熟女国产| 久久精品91无色码中文字幕| 黄色a级毛片大全视频| 久久久久亚洲av毛片大全| 日本a在线网址| 国产又黄又爽又无遮挡在线| 亚洲 国产 在线| 亚洲欧美一区二区三区黑人| 伦理电影免费视频| 国产高清视频在线播放一区| 丝袜美腿诱惑在线| 午夜成年电影在线免费观看| 动漫黄色视频在线观看| 欧美日韩中文字幕国产精品一区二区三区| 亚洲免费av在线视频| 一二三四在线观看免费中文在| 波多野结衣高清作品| aaaaa片日本免费| 欧美不卡视频在线免费观看 | 特大巨黑吊av在线直播| 亚洲国产高清在线一区二区三| 我要搜黄色片| 久久午夜综合久久蜜桃| 在线a可以看的网站| 麻豆国产97在线/欧美 | 看黄色毛片网站| 国产一区二区在线观看日韩 | 成人18禁在线播放| 亚洲成人久久性| 性色av乱码一区二区三区2| 久久人妻av系列| 欧美又色又爽又黄视频| 亚洲黑人精品在线| 怎么达到女性高潮| 亚洲美女黄片视频| 久久久久九九精品影院| 99久久99久久久精品蜜桃| 18禁观看日本| 香蕉丝袜av| 午夜久久久久精精品| 欧美+亚洲+日韩+国产| 99久久无色码亚洲精品果冻| 国产探花在线观看一区二区| 美女大奶头视频| 成人手机av| 免费在线观看视频国产中文字幕亚洲| 免费电影在线观看免费观看| 国产视频内射| 国产亚洲精品久久久久久毛片| 亚洲性夜色夜夜综合| 正在播放国产对白刺激| 国产三级中文精品| 久久亚洲真实| 成人精品一区二区免费| 久久精品国产亚洲av香蕉五月| 亚洲人成电影免费在线| 激情在线观看视频在线高清| 午夜激情福利司机影院| 久久久久国产精品人妻aⅴ院| 亚洲18禁久久av| 精品欧美一区二区三区在线| 久久久久久免费高清国产稀缺| 国产精品av久久久久免费| 中文在线观看免费www的网站 | 国产激情久久老熟女| 18禁黄网站禁片免费观看直播| 精品无人区乱码1区二区| 此物有八面人人有两片| 亚洲成人久久爱视频| 国产精品久久久久久久电影 | 亚洲av中文字字幕乱码综合| avwww免费| 免费在线观看成人毛片| 亚洲成人免费电影在线观看| 日韩欧美免费精品| 啦啦啦观看免费观看视频高清| 亚洲成人中文字幕在线播放| tocl精华| 欧美精品亚洲一区二区| 精品久久久久久久毛片微露脸| 91老司机精品| 色在线成人网| 一级作爱视频免费观看| 在线a可以看的网站| 叶爱在线成人免费视频播放| 亚洲在线自拍视频| 国产一区在线观看成人免费| 一进一出抽搐动态| 国产亚洲av高清不卡| 久久人妻福利社区极品人妻图片| xxx96com| av欧美777| 国产69精品久久久久777片 | 成年版毛片免费区| 中文字幕熟女人妻在线| 久久久久久久久免费视频了| 啦啦啦免费观看视频1| 无人区码免费观看不卡| 夜夜看夜夜爽夜夜摸| 最近最新中文字幕大全免费视频| 亚洲av中文字字幕乱码综合| 欧美午夜高清在线| 欧美一级毛片孕妇| 大型黄色视频在线免费观看| 午夜免费激情av| 欧美不卡视频在线免费观看 | 在线观看舔阴道视频| a级毛片a级免费在线| 成人精品一区二区免费| 88av欧美| www.熟女人妻精品国产| 国产在线观看jvid| 免费在线观看亚洲国产| 麻豆久久精品国产亚洲av| 欧美成人一区二区免费高清观看 | 丁香六月欧美| 亚洲成人国产一区在线观看| 精品欧美一区二区三区在线| www.自偷自拍.com| 欧美色视频一区免费| 黑人巨大精品欧美一区二区mp4| 亚洲精品中文字幕一二三四区| 人人妻人人看人人澡| 国产欧美日韩精品亚洲av| 最近最新中文字幕大全电影3| 国产不卡一卡二| 国产成人欧美在线观看| 日韩精品青青久久久久久| 国产欧美日韩一区二区精品| 日本成人三级电影网站| 亚洲aⅴ乱码一区二区在线播放 | 亚洲国产欧洲综合997久久,| 在线播放国产精品三级| 中文字幕人成人乱码亚洲影| 亚洲一区高清亚洲精品| 国产av在哪里看| 亚洲人与动物交配视频| 麻豆久久精品国产亚洲av| 18美女黄网站色大片免费观看| 久久久久国产精品人妻aⅴ院| 丝袜美腿诱惑在线| netflix在线观看网站| 露出奶头的视频| 国产99久久九九免费精品| 色综合亚洲欧美另类图片| 日本三级黄在线观看| 国产又色又爽无遮挡免费看| 国产成人欧美在线观看| 人人妻人人澡欧美一区二区| 国产成人av教育| 日本免费一区二区三区高清不卡| 亚洲成av人片免费观看| 国产亚洲精品久久久久久毛片| 最近最新免费中文字幕在线| 激情在线观看视频在线高清| 男女之事视频高清在线观看| 久久精品亚洲精品国产色婷小说| 久久久精品欧美日韩精品| 波多野结衣巨乳人妻| 日韩三级视频一区二区三区| 淫妇啪啪啪对白视频| 美女黄网站色视频| 制服人妻中文乱码| 国产私拍福利视频在线观看| 国产人伦9x9x在线观看| 国产高清激情床上av| 美女扒开内裤让男人捅视频| 一a级毛片在线观看| 欧美成人免费av一区二区三区| 亚洲国产精品999在线| 性欧美人与动物交配| 亚洲自拍偷在线| 国产成人欧美在线观看| 日韩欧美精品v在线| 久久久久久亚洲精品国产蜜桃av| av天堂在线播放| 日本撒尿小便嘘嘘汇集6| 日韩精品青青久久久久久| 热99re8久久精品国产| www.自偷自拍.com| 精品无人区乱码1区二区| 日韩欧美国产一区二区入口| 精品无人区乱码1区二区| 男人舔奶头视频| aaaaa片日本免费| 久久人人精品亚洲av| 狂野欧美激情性xxxx| 国产在线观看jvid| 午夜福利成人在线免费观看| 国内揄拍国产精品人妻在线| 亚洲国产精品久久男人天堂| 国产精品久久久人人做人人爽| 正在播放国产对白刺激| 美女黄网站色视频| 欧美久久黑人一区二区| 亚洲成人国产一区在线观看| 黑人欧美特级aaaaaa片| 亚洲国产欧洲综合997久久,| 国产探花在线观看一区二区| 国产一区二区在线av高清观看| 真人做人爱边吃奶动态| 精品久久久久久久毛片微露脸|