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

    病媒生物空間抽樣方法應(yīng)用研究進(jìn)展及對蚊蟲監(jiān)測的啟示

    2022-02-27 00:44:08余可憶吳健平周毅彬姚申君
    關(guān)鍵詞:釘螺蚊蟲異質(zhì)性

    余可憶 吳健平 周毅彬 姚申君*

    (1.地理信息科學(xué)教育部重點(diǎn)實(shí)驗(yàn)室,上海 200241;2.華東師范大學(xué)地理科學(xué)學(xué)院,上海 200241;3.自然資源部超大城市自然資源時空大數(shù)據(jù)分析應(yīng)用重點(diǎn)實(shí)驗(yàn)室,上海 200241;4.上海市疾病預(yù)防控制中心傳染病防治所,上海 200336)

    媒介生物一般指的是作為疾病傳播途徑且能將傳染病傳播給人類的生物。作為媒介生物的一種,蚊蟲不僅刺叮騷擾人群,給人們帶來較差的生活體驗(yàn),而且可以傳播登革熱、黃熱病、瘧疾、淋巴絲蟲病、流行性乙型腦炎等蚊媒傳染病,是需要重點(diǎn)監(jiān)測的媒介生物(閆冬明等,2020;田波等,2022)。過去100年來,全球氣溫升高加快了蚊媒疾病的快速傳播;蓬勃發(fā)展的經(jīng)濟(jì)建設(shè)與快速的城市化進(jìn)程使蚊媒疾病的傳播模式更為復(fù)雜(Anoopkumaretal.,2022)。蚊蟲傳染病的頻繁暴發(fā)(Yueetal.,2014;牟笛等,2017),使其成為當(dāng)下亟待解決的公共衛(wèi)生問題之一,相關(guān)防治工作也受到越來越多的關(guān)注與重視(冷培恩等,2019)。

    通過對媒介蚊蟲密度的監(jiān)測以指導(dǎo)開展控制,可以有效阻斷蚊媒傳染病的傳播,從而遏制其進(jìn)一步蔓延。我國現(xiàn)有關(guān)于蚊蟲監(jiān)測的研究多著重于對監(jiān)測工具、監(jiān)測方法、監(jiān)測結(jié)果的闡述,以及對其結(jié)果進(jìn)行季節(jié)性特征分析(冷培恩等,2019;劉美德等,2021;王純玉等,2022;王巧燕等,2022)。近年來,學(xué)者開始關(guān)注監(jiān)測點(diǎn)的抽樣問題,并逐漸意識到空間自相關(guān)和空間異質(zhì)性的重要影響,如有學(xué)者以社區(qū)劃分調(diào)查網(wǎng)格進(jìn)行監(jiān)測(張佳一等,2021),依據(jù)居民區(qū)、公園等不同的生境設(shè)置監(jiān)測點(diǎn)(王巧燕等,2022),根據(jù)地形、植被、水文和土地覆蓋/利用等進(jìn)行環(huán)境分層實(shí)現(xiàn)抽樣(Rouxetal.,2013;夏斯偉等,2021)。然而蚊蟲監(jiān)測的空間抽樣研究尚在起步階段,現(xiàn)有蚊蟲監(jiān)測點(diǎn)的布樣主要以選擇蚊蟲孳生的最佳生境或疫情高發(fā)區(qū)最為常見(周毅彬等,2021),需要更為詳細(xì)的規(guī)范,深入理解蚊蟲空間分布特征,實(shí)現(xiàn)精準(zhǔn)防控。亟需學(xué)者在理論、方法與應(yīng)用等多方面進(jìn)行探索,以豐富蚊蟲監(jiān)測的理論與實(shí)踐。

    除蚊蟲以外,鼠類和釘螺也是兩種較為重要的病媒生物,危害著世界人民的公共衛(wèi)生安全。鼠類不僅危害人類健康,也在一定程度上影響了農(nóng)林畜牧業(yè)的發(fā)展(黃英凱等,2015)。湖北釘螺是日本血吸蟲的唯一中間宿主。由血吸蟲引起的血吸蟲病(Donohueetal.,2017)威脅著我國長江流域及其以南地區(qū)人民健康(郭蘇影等,2021)??臻g抽樣方法已被廣泛應(yīng)用于鼠類和釘螺監(jiān)測,對及時且精確地做好滅鼠工作、阻斷血吸蟲病的傳播具有重要的意義(李忠武等,2013;雷正龍等,2015)。

    本文旨在闡述空間抽樣理論在病媒生物監(jiān)測中的應(yīng)用,并歸納適用于蚊蟲監(jiān)測的空間抽樣框架以供相關(guān)研究借鑒。下文將先簡單介紹空間抽樣理論,再論述空間抽樣方法在蚊蟲、鼠類和釘螺監(jiān)測上的應(yīng)用,最后提出適用于蚊蟲監(jiān)測的空間抽樣方案設(shè)想。

    1 空間抽樣與統(tǒng)計推斷

    經(jīng)典抽樣方法的前提是樣本單元之間是互相獨(dú)立同分布的,但絕大多數(shù)分布在空間的研究對象(如蚊蟲等媒介生物),通常具有空間相關(guān)性以及空間異質(zhì)性。因此,從具有特殊分布性質(zhì)的研究對象中采用經(jīng)典抽樣方法得到的樣本,往往是非獨(dú)立性、非均質(zhì)性且具有屬性值隨著空間的不同劃分而變化的特點(diǎn),會對抽樣統(tǒng)計結(jié)果產(chǎn)生一系列的有偏估計(王勁峰等,2019)。如空間自相關(guān)會對用樣本簡單平均值估計區(qū)域均值以及超總體均值產(chǎn)生影響,導(dǎo)致樣本均值方差的偏小或偏大(Haining,1988;Cressie,1991;Griffithetal.,1994;Haining,2003;Wangetal.,2012);當(dāng)研究對象具有空間異質(zhì)性,采用分層抽樣但某些層沒有樣本時,樣本簡單平均值的數(shù)學(xué)期望將有偏于總體均值(Wangetal.,2011;Huetal.,2013;Wangetal.,2013a)。因此,在抽樣調(diào)查具有特殊空間分布的對象時,更適合采用將空間數(shù)據(jù)特性納入數(shù)理統(tǒng)計范圍內(nèi)的空間抽樣方法(Atkinson,1991;Stehmanetal.,2003;Wangetal.,2012)。

    空間抽樣在經(jīng)典抽樣方法的基礎(chǔ)上更加強(qiáng)調(diào)空間的重要性。與經(jīng)典抽樣不同的是,空間抽樣的對象具有地理空間坐標(biāo),如果改變被抽取空間對象的位置,將產(chǎn)生完全不同的抽樣調(diào)查總體,從而影響抽樣調(diào)查結(jié)果(姜成晟等,2009)??臻g抽樣調(diào)查對象由地理空間坐標(biāo)和屬性值共同組成,因此在空間抽樣過程中需考慮研究對象的空間分布特征:是否獨(dú)立同分布、是否具有空間自相關(guān)特征、是否具有分層異質(zhì)性特征、抑或既具有空間自相關(guān)又具有分層異質(zhì)性特征(王勁峰等,2019)。20世紀(jì)末至今,學(xué)界普遍認(rèn)為由空間異質(zhì)性特征對空間進(jìn)行分層采樣可以實(shí)現(xiàn)更有效的空間采樣和統(tǒng)計推斷(Lietal.,2008;Wangetal.,2010)。在考慮到空間異質(zhì)性、相關(guān)性以及小樣本多報告單元的問題,王勁峰等于2013年正式提出了空間三明治抽樣(Wangetal.,2013b)。隨著空間抽樣的理論與模型逐步完善,國內(nèi)外學(xué)者將其廣泛運(yùn)用于林業(yè)(董士偉,2018)、農(nóng)業(yè)(王迪等,2021)、土壤(陳超等,2021)、生物種群(Liuetal.,2018)、土地利用(馬煒,2019)調(diào)查等多個領(lǐng)域。

    在大多數(shù)研究中,空間抽樣的目的是對總體或超總體進(jìn)行估算(如蚊蟲的數(shù)量、鼠密度、釘螺密度等),方法是從總體中以某種方式抽取一定數(shù)量的樣本單元(如使用簡單隨機(jī)抽樣或分層抽樣等),然后將獲取的樣本匯總得到樣本的統(tǒng)計量(如算數(shù)平均等),利用上述統(tǒng)計量并結(jié)合統(tǒng)計推斷方法,計算得到可以表征總體的統(tǒng)計量,完成對總體的估算(王勁峰等,2019)。上述過程主要涉及確定地學(xué)對象的屬性值、抽樣方案(樣本量與布樣情況)和統(tǒng)計量(王勁峰等,2019),有著眾多的選擇與組合。如何選擇與研究對象性質(zhì)相符合的抽樣方法以及統(tǒng)計量,以期使用較少的樣本得到精度較高的統(tǒng)計結(jié)果,是所有學(xué)者在進(jìn)行空間抽樣前需要解決的問題。為此,王勁峰等(2019)采用屬性值、抽樣方案和統(tǒng)計量所組成的“三位一體”的方式對空間抽樣與統(tǒng)計推斷加以描述,提供了一個空間抽樣理論框架,根據(jù)研究對象的空間分布特性,空間抽樣方法以及統(tǒng)計推斷方式可有不同選擇。如當(dāng)研究對象既具有空間相關(guān)性又具有空間異質(zhì)性,且采用分層抽樣時,在每層樣本足夠多的情況下,可以采用異質(zhì)表面均值估計(Mean of Surface with Nonhomogeneity,MSN)及其抽樣(Wangetal.,2009;Huetal.,2011);當(dāng)有些層的樣本缺失時,即樣本有偏的情況下,應(yīng)采用Bshade統(tǒng)計量(Wangetal.,2011;Huetal.,2013;Xuetal.,2013);若僅有一個樣本點(diǎn)時,則需要使用單點(diǎn)面積估計(Single Point Areal Estimation,SPA),并結(jié)合輔助信息進(jìn)行統(tǒng)計推斷(Wangetal.,2013c;王勁峰等,2019)。

    后文將按照不同的主流空間抽樣方法,對蚊蟲、鼠類和釘螺三種病媒生物的監(jiān)測應(yīng)用展開綜述。

    2 空間分層抽樣

    空間分層抽樣是較為常用且簡便易行的一種空間抽樣方法,大多數(shù)用于存在空間異質(zhì)性的研究對象,可在一定程度上減少抽樣樣本數(shù),進(jìn)而降低抽樣成本。由于釘螺和蚊蟲的空間分布特征與地理環(huán)境密切相關(guān),而地理環(huán)境均存在著一定的空間異質(zhì)性,因此多數(shù)學(xué)者均采用了空間分層抽樣對上述兩種病媒生物進(jìn)行數(shù)量監(jiān)測。本節(jié)將先介紹空間分層抽樣原理,再對釘螺與蚊蟲監(jiān)測進(jìn)行總結(jié)梳理。

    2.1 原理

    空間分層抽樣應(yīng)先按照一定的規(guī)則將總體劃分為多層,使得屬性值相近的點(diǎn)被劃分為同一層,即層內(nèi)方差最小、層間方差盡可能大(Wangetal.,2010;Wangetal.,2016)。這里的規(guī)則視不同的研究對象而定,須遵循該對象的空間分布特征。

    再根據(jù)一定的分配原則,將總樣本量分配到各個層內(nèi);最后根據(jù)各層內(nèi)的樣本量,在各層內(nèi)部進(jìn)行簡單隨機(jī)抽樣。分配原則包括:各層平均分配;按各層單元數(shù)比例分配;按某層標(biāo)準(zhǔn)離散方差Sh與該層單元數(shù)Nh乘積比例分配,當(dāng)各層的樣本單元成本Ch不同時,層內(nèi)最優(yōu)抽樣如公式1所示(Cochran,1977)。上述三種分配原則的抽樣效率依次提高。

    (1)

    2.2 釘螺

    由于釘螺孳生與地理環(huán)境密切相關(guān),土壤、植被等的空間異質(zhì)性導(dǎo)致了距離相鄰的采樣點(diǎn)的釘螺數(shù)可能相差很大(邱娟,2012;張文馨,2016),因此空間分層抽樣成為了釘螺調(diào)查的首選。通過該方法,可將數(shù)理統(tǒng)計意義上同一生境下的研究區(qū)進(jìn)行歸類,使得層內(nèi)同質(zhì)性較高、層間異質(zhì)性較高,再對各層進(jìn)行合理抽樣。在確保無偏采樣的同時,獲得更具有代表性的樣本,從而對此研究區(qū)的釘螺數(shù)量做出較為準(zhǔn)確的估計。

    張文馨(2015)選擇江西省新建縣恒湖農(nóng)場茶葉港的一塊50 m×50 m草洲作為試驗(yàn)樣地,通過比較3種不同抽樣方法(簡單隨機(jī)抽樣、系統(tǒng)抽樣和空間分層抽樣)所需樣本量以及抽樣的相對和絕對誤差,發(fā)現(xiàn)依據(jù)高程作為輔助因子進(jìn)行空間分層抽樣所需的樣本量最少,絕對誤差最小,優(yōu)于傳統(tǒng)抽樣方法。商伊迪(2017)也采用了空間分層抽樣策略對鄱陽湖的釘螺進(jìn)行調(diào)查。由于灰化苔草和南荻的高度、蓋度均與釘螺孳生具有明顯空間關(guān)系,將上述植被因子作為輔助變量,通過K均值進(jìn)行聚類,將其分為7層,再依據(jù)奈曼分配計算總樣本量以及各層樣本量,最后對各層進(jìn)行隨機(jī)抽樣。汪訓(xùn)平(2018)提出了一種基于植物豐度的空間分層抽樣策略,將與釘螺孳生密切相關(guān)的植物豐度作為分層依據(jù),采用Hammond McCullagh方程、依據(jù)植物豐度的變異系數(shù),確定每層最優(yōu)采樣點(diǎn)的數(shù)量,最后使用多向插值離散度的空間采樣布局方案,以覆蓋最大化為原則在各層中進(jìn)行合理布樣。

    2.3 蚊蟲

    由于蚊蟲密度不僅與植被覆蓋(陰涼處)、積水分布、居民密度及人類活動(使用殺蟲劑等)有關(guān),也會受溫度、降雨量、濕度等氣象因素影響(仲潔等,2015)。因此即使是相同的土地利用類型,也會存在一定的空間異質(zhì)性,即對蚊蟲進(jìn)行監(jiān)測時也可使用空間分層抽樣(Ageepetal.,2009;Rouxetal.,2013;Seddaetal.,2019)。

    Ageep等(2009)對蘇丹北部阿拉伯按蚊的時空分布展開研究,主要采用了空間分層隨機(jī)抽樣。先利用多光譜QuickBird數(shù)據(jù)以及高分辨率數(shù)字高程模型,依據(jù)物體的形狀、紋理等對研究區(qū)土地利用類型進(jìn)行分層合并,得到5個層復(fù)合定義的土地利用類型。之后每個月在上述層內(nèi)隨機(jī)抽取40個像元作為樣本,每個樣本面積為1公頃。Roux等(2013)設(shè)計了一種基于環(huán)境分層的成年按蚊調(diào)查和棲息地特征的采樣策略。首先需要在研究區(qū)中隨機(jī)抽樣確定500個地點(diǎn),作為后續(xù)聚類的候選采樣點(diǎn)。接著擬定在地形、植被、水文和土地覆蓋/利用等大類下的共43個不同的環(huán)境變量。然后對上述環(huán)境變量進(jìn)行主成分分析以及混合群的因子分析,確定影響分層的決定性變量。最后利用500個采樣候選點(diǎn)以及決定性變量,對其進(jìn)行K均值聚類,最終獲得6個預(yù)期集群(層)。Sedda等(2019)提出了一個生態(tài)采樣設(shè)計框架,該框架主要討論了樣本大小、如何分層以及采樣點(diǎn)的地理分配3個要素。首先,通過計算已知協(xié)方差函數(shù)的Log Gaussian Cox Process模型,計算上述模型參數(shù)的總方差以及預(yù)測蚊蟲數(shù)的標(biāo)準(zhǔn)誤差,得到最優(yōu)樣本大小。然后,利用二次判別分析算法,對研究區(qū)進(jìn)行分層,并得到最佳層數(shù)以及其地理界限。最后,依據(jù)每一層的樣本點(diǎn)的數(shù)量在該層內(nèi)進(jìn)行隨機(jī)抽樣或自適應(yīng)抽樣。

    國內(nèi)外對釘螺與蚊蟲進(jìn)行數(shù)量監(jiān)測時大多使用了空間分層抽樣,或提供了一個基于空間分層抽樣的抽樣模型。大多數(shù)學(xué)者均沿用了空間分層抽樣的理念,但對如何進(jìn)行分層、如何確定最佳樣本數(shù)以及如何對各層進(jìn)行抽樣等問題進(jìn)行了更深入的研究與創(chuàng)新,使其更加適用于釘螺調(diào)查與蚊蟲密度監(jiān)測。

    3 空間三明治抽樣

    空間分層抽樣方法主要針對單報告單元,當(dāng)報告單元較多時,對每個報告單元進(jìn)行分別抽樣,將會產(chǎn)生樣本量大、費(fèi)用高的問題。空間三明治抽樣模型(Wangetal.,2002;Wangetal.,2013b)解決了在總體分異條件下小樣本多報告單元抽樣問題,可以用較少的樣本量實(shí)現(xiàn)多報告單元的匯總。另外,由于空間三明治抽樣對于樣本圖層中采樣點(diǎn)的空間分布沒有要求,因此可以最大限度地利用現(xiàn)有先驗(yàn)驗(yàn)數(shù)據(jù)。

    此抽樣主要分為以下3個步驟(Wangetal.,2002;Wangetal.,2013b):首先將研究對象按照層間方差最大,層內(nèi)方差最小進(jìn)行分層,形成知識層,如圖1所示,可按照研究區(qū)域的屬性對層進(jìn)行劃分。接著,將樣本按照知識層進(jìn)行分配,計算出各層的樣本均值以及樣本均值方差。最后將知識層的均值以及均值方差推算至各報告單元,得到各報告單元的總體情況。

    圖1 空間三明治抽樣(Wang et al.,2013b)Fig.1 Spatial sandwich sampling

    從上述步驟可以看出,空間三明治抽樣實(shí)際上是空間分層抽樣的變體,采用了兩次分層抽樣的技術(shù),第一次是從樣本圖層到知識層,第二次是從知識層到報告圖層。

    該方法主要用于對研究區(qū)內(nèi)鼠密度進(jìn)行抽樣調(diào)查與制圖。如任周鵬(2011)以三峽庫區(qū)現(xiàn)有鼠密度監(jiān)測數(shù)據(jù),并選擇土壤、地貌以及NDVI(normalized difference vegetation index,歸一化植被指數(shù))作為估計研究區(qū)各縣平均鼠密度的輔助因子,將NDVI的分類結(jié)果作為三明治模型的知識層,完成三峽庫區(qū)蓄水前后各縣平均鼠密度的估計,并采用以行政區(qū)為單元和以格網(wǎng)為單元的兩種報告方式進(jìn)行制圖。Liu等(2018)在中國吉林省進(jìn)行鼠密度調(diào)查時發(fā)現(xiàn),土壤類型、土地利用方式以及植被類型對鼠類密度分布有顯著影響,且研究區(qū)域異質(zhì)性較高。因此通過地理探測器模型,在評估各種自然因素和人為因素對鼠類分布的影響及其交互效應(yīng)的情況下,得到空間分層分異的情況。通過三明治模型分層計算知識層的鼠密度并匯總至報告層。

    與傳統(tǒng)抽樣方法相比,空間三明治抽樣方法的主要優(yōu)勢是考慮到媒介生物密度存在空間自相關(guān)性與空間異質(zhì)性,并且在使用現(xiàn)有監(jiān)測點(diǎn)數(shù)據(jù)的情況下,解決總體分異條件下小樣本多報告單元抽樣問題,統(tǒng)計推斷得到任意單元內(nèi)的鼠類密度,完成對鼠密度的多尺度匯總。

    4 對蚊蟲監(jiān)測的啟示

    通過上述綜述可以發(fā)現(xiàn),空間抽樣理論體系已構(gòu)建完整,很多成熟空間抽樣方法已被廣泛使用于媒介生物監(jiān)測,但國內(nèi)有關(guān)蚊蟲的監(jiān)測大多采用棲息地分層的方式,并且大多是基于主觀經(jīng)驗(yàn)或定性進(jìn)行分層(Ageepetal.,2009),但相同類型的棲息地也有可能存在一定的空間異質(zhì)性,而不同類型的棲息地之間也有可能有著一定的空間自相關(guān)。因此未確定可表征研究區(qū)環(huán)境特征的變量、未檢測樣本在時空上的獨(dú)立性以及忽略了相同類型棲息地的異質(zhì)性,都有可能使監(jiān)測結(jié)果產(chǎn)生偏差。

    在蚊蟲的監(jiān)測過程中引入空間分層抽樣與空間三明治抽樣可提高抽樣的準(zhǔn)確性,在一定程度上解決上述問題。在使用空間分層抽樣過程中,首先應(yīng)確定分層的規(guī)則,規(guī)則的優(yōu)劣會對分層抽樣結(jié)果的準(zhǔn)確度產(chǎn)生一定影響。此規(guī)則往往和影響蚊蟲數(shù)量的環(huán)境因子、氣象因素、人口密度等相關(guān)聯(lián),可以通過先驗(yàn)知識或混合群的因子分析得到。在分層前,可使用主成分分析,確保各因素間相互獨(dú)立。之后,可使用地理探測器(Wangetal.,2012),探究上述因素是否對蚊蟲密度有重要影響??墒褂肒均值聚類或模糊C均值聚類法進(jìn)行分層,使得屬性值相近的點(diǎn)被劃分為同一層;再根據(jù)一定的分配原則,將總樣本量分配到各個層內(nèi);最后在各層內(nèi)部進(jìn)行簡單隨機(jī)抽樣。在得到抽樣結(jié)果后,可使用地理加權(quán)回歸(geographically weighted regression,GWR)(Brunsdonetal.,1998;Fotheringhametal.,1999)對各采樣點(diǎn)的影響因子(自變量)與蚊蟲數(shù)量(因變量)進(jìn)行回歸分析,擬合得到回歸方程。從上述回歸方程可以解釋蚊蟲數(shù)量的關(guān)鍵變量或?qū)ρ芯繀^(qū)內(nèi)各地蚊蟲數(shù)量進(jìn)行預(yù)測。而空間三明治抽樣方法的應(yīng)用,則可以最大限度地利用先驗(yàn)數(shù)據(jù),完成蚊蟲多報告單元的匯總。

    5 總結(jié)與展望

    蚊蟲作為一種可以傳播多種傳染病的媒介生物,其監(jiān)測工作對于控制蚊蟲數(shù)量、預(yù)防相關(guān)傳染病具有重要作用。關(guān)于蚊蟲密度監(jiān)測方法、每種方法的適用范圍、操作步驟以及密度的計算方法,國家標(biāo)準(zhǔn)中均有詳細(xì)的規(guī)定,但缺少與蚊蟲監(jiān)測點(diǎn)布樣相關(guān)的規(guī)范。如今已有部分學(xué)者開始意識到空間分布特征對于蚊蟲監(jiān)測點(diǎn)布樣的影響,已開始積極引入空間抽樣理論,在考慮樣本的空間相關(guān)性和異質(zhì)性的同時,增加對異質(zhì)性區(qū)域的監(jiān)測,減少對同質(zhì)區(qū)域的重復(fù)采樣,提高了抽樣效率和質(zhì)量(張杰等,2016)。然而國內(nèi)有關(guān)蚊蟲的監(jiān)測尚處于初步探索階段,可根據(jù)特征選擇合適的空間抽樣方法,通過借鑒、融合、發(fā)展已有方法,在蚊蟲監(jiān)測的理論與實(shí)踐上進(jìn)行更深一步的創(chuàng)新探究。

    猜你喜歡
    釘螺蚊蟲異質(zhì)性
    辛巴狗幽默日常
    2020年洪澇災(zāi)害后湖北省釘螺擴(kuò)散情況調(diào)查
    帶您了解蚊蟲驅(qū)避劑
    2016—2021年池州市貴池區(qū)新發(fā)和復(fù)現(xiàn)釘螺情況分析
    夏天來了 實(shí)用防蚊蟲妙招請收好
    基于可持續(xù)發(fā)展的異質(zhì)性債務(wù)治理與制度完善
    一種快速鑒別湖北釘螺(Oncomelania hupensis)雌雄性的方法
    快樂大講堂③蚊蟲叮咬止癢法
    現(xiàn)代社區(qū)異質(zhì)性的變遷與啟示
    1949年前譯本的民族性和異質(zhì)性追考
    最近手机中文字幕大全| 国产69精品久久久久777片| 国产视频内射| 观看免费一级毛片| 性高湖久久久久久久久免费观看| 男女边摸边吃奶| 日韩精品有码人妻一区| 卡戴珊不雅视频在线播放| 美女内射精品一级片tv| 国产精品久久久久久av不卡| 丰满迷人的少妇在线观看| 美女主播在线视频| 精品亚洲成国产av| 99热这里只有精品一区| 中文字幕免费在线视频6| 欧美成人午夜免费资源| 男人添女人高潮全过程视频| 亚洲不卡免费看| 久久精品国产亚洲av天美| 亚洲人与动物交配视频| 国产 精品1| 少妇熟女欧美另类| 人妻系列 视频| 久久99热6这里只有精品| 亚洲成人一二三区av| 国产精品伦人一区二区| 美女视频免费永久观看网站| 国产在线免费精品| 久久99热这里只有精品18| 精品少妇久久久久久888优播| 国产精品爽爽va在线观看网站| 尤物成人国产欧美一区二区三区| 国产精品国产三级国产专区5o| 99久国产av精品国产电影| 国产精品久久久久久av不卡| 一区二区三区四区激情视频| av在线播放精品| 日韩成人av中文字幕在线观看| 大片免费播放器 马上看| 免费黄色在线免费观看| 国产精品人妻久久久久久| 日韩人妻高清精品专区| 免费久久久久久久精品成人欧美视频 | 久久久午夜欧美精品| 九色成人免费人妻av| 国产在视频线精品| 男人添女人高潮全过程视频| 这个男人来自地球电影免费观看 | 欧美日韩视频精品一区| 高清av免费在线| 婷婷色综合www| 日日啪夜夜爽| 人人妻人人澡人人爽人人夜夜| 国产一区二区在线观看日韩| 免费观看性生交大片5| 欧美3d第一页| 日韩欧美一区视频在线观看 | 久久久久久久久大av| 久久久成人免费电影| 久久精品国产亚洲网站| 大片电影免费在线观看免费| 国产精品女同一区二区软件| 99久国产av精品国产电影| 一级黄片播放器| 97在线人人人人妻| 国产成人一区二区在线| 大陆偷拍与自拍| 黄色日韩在线| 性高湖久久久久久久久免费观看| 久久亚洲国产成人精品v| 联通29元200g的流量卡| 91狼人影院| 精品国产一区二区三区久久久樱花 | 亚洲精品,欧美精品| 中文天堂在线官网| 国产精品福利在线免费观看| 日本色播在线视频| 亚洲精品国产成人久久av| 日韩亚洲欧美综合| 欧美丝袜亚洲另类| 热re99久久精品国产66热6| 亚洲av在线观看美女高潮| 日韩在线高清观看一区二区三区| 精品酒店卫生间| 成人毛片a级毛片在线播放| 男女啪啪激烈高潮av片| 日日啪夜夜撸| 看免费成人av毛片| 久久人妻熟女aⅴ| 纵有疾风起免费观看全集完整版| 国产精品av视频在线免费观看| 麻豆国产97在线/欧美| 久久人妻熟女aⅴ| 伊人久久国产一区二区| 欧美一级a爱片免费观看看| 亚洲成人av在线免费| 日本与韩国留学比较| 又大又黄又爽视频免费| 亚州av有码| 日韩强制内射视频| 啦啦啦中文免费视频观看日本| 六月丁香七月| 少妇熟女欧美另类| 高清欧美精品videossex| av国产精品久久久久影院| 精品久久久精品久久久| 国产免费一区二区三区四区乱码| 夜夜看夜夜爽夜夜摸| 中文天堂在线官网| 日本午夜av视频| 91精品一卡2卡3卡4卡| 成人亚洲欧美一区二区av| 国产乱人视频| 亚洲精品乱码久久久久久按摩| 日韩一区二区三区影片| 韩国av在线不卡| videos熟女内射| www.av在线官网国产| 亚洲精品aⅴ在线观看| 少妇人妻一区二区三区视频| 日韩三级伦理在线观看| 午夜福利视频精品| 狂野欧美白嫩少妇大欣赏| 国产乱人偷精品视频| 建设人人有责人人尽责人人享有的 | 97超碰精品成人国产| 精品国产一区二区三区久久久樱花 | 汤姆久久久久久久影院中文字幕| 在线观看三级黄色| 美女内射精品一级片tv| 性色av一级| 成人美女网站在线观看视频| 老师上课跳d突然被开到最大视频| 国产精品秋霞免费鲁丝片| 大片电影免费在线观看免费| 日韩亚洲欧美综合| 久久国产精品男人的天堂亚洲 | 九草在线视频观看| 免费观看av网站的网址| 五月伊人婷婷丁香| 少妇被粗大猛烈的视频| 国产成人精品一,二区| 久久毛片免费看一区二区三区| 日日摸夜夜添夜夜添av毛片| 成人毛片60女人毛片免费| 国产精品人妻久久久影院| 王馨瑶露胸无遮挡在线观看| 久久精品国产亚洲av天美| 伦理电影大哥的女人| 精品国产三级普通话版| 亚洲精品久久午夜乱码| 国产无遮挡羞羞视频在线观看| 日韩欧美一区视频在线观看 | 内地一区二区视频在线| 国产久久久一区二区三区| 亚洲人成网站高清观看| 在线观看一区二区三区激情| 日韩电影二区| 国产精品一区二区三区四区免费观看| 91午夜精品亚洲一区二区三区| 亚洲精品视频女| 人妻系列 视频| 亚洲,一卡二卡三卡| 王馨瑶露胸无遮挡在线观看| 美女cb高潮喷水在线观看| 一级爰片在线观看| 国产高潮美女av| 97在线视频观看| 在线 av 中文字幕| 精品亚洲成a人片在线观看 | 国产成人91sexporn| 王馨瑶露胸无遮挡在线观看| 最近最新中文字幕大全电影3| 久久久久久久久久人人人人人人| 亚洲精品视频女| 亚洲一级一片aⅴ在线观看| 欧美日韩视频高清一区二区三区二| a 毛片基地| 超碰av人人做人人爽久久| 成年人午夜在线观看视频| 新久久久久国产一级毛片| 亚洲成人av在线免费| 精品一区二区免费观看| 天天躁日日操中文字幕| 欧美丝袜亚洲另类| 国产伦精品一区二区三区四那| 色网站视频免费| 看免费成人av毛片| kizo精华| 久久久亚洲精品成人影院| 免费不卡的大黄色大毛片视频在线观看| 久久人人爽人人爽人人片va| 欧美三级亚洲精品| 各种免费的搞黄视频| 在现免费观看毛片| 九九爱精品视频在线观看| 国产白丝娇喘喷水9色精品| 欧美精品人与动牲交sv欧美| 中国三级夫妇交换| 免费观看无遮挡的男女| 一二三四中文在线观看免费高清| 国产高清不卡午夜福利| 国产午夜精品久久久久久一区二区三区| 高清黄色对白视频在线免费看 | 夫妻性生交免费视频一级片| 中文字幕av成人在线电影| 国产亚洲午夜精品一区二区久久| 毛片一级片免费看久久久久| 国产v大片淫在线免费观看| 国产精品国产av在线观看| 永久网站在线| 少妇人妻 视频| 国国产精品蜜臀av免费| 日日摸夜夜添夜夜爱| 欧美精品亚洲一区二区| 午夜激情福利司机影院| 日韩强制内射视频| 亚洲电影在线观看av| 这个男人来自地球电影免费观看 | 中文欧美无线码| 久久久成人免费电影| 最近中文字幕高清免费大全6| 熟女人妻精品中文字幕| 国产日韩欧美在线精品| 人妻夜夜爽99麻豆av| av免费在线看不卡| 一区二区三区四区激情视频| 国产91av在线免费观看| www.av在线官网国产| 亚洲精品日韩av片在线观看| 我要看日韩黄色一级片| 免费大片18禁| 国产成人91sexporn| 91久久精品国产一区二区成人| 国产久久久一区二区三区| 国产精品久久久久久久电影| 中文欧美无线码| 久久人人爽人人爽人人片va| 久久精品久久精品一区二区三区| 国产精品熟女久久久久浪| 亚洲av综合色区一区| 国产精品99久久久久久久久| 国产大屁股一区二区在线视频| 亚洲国产高清在线一区二区三| 成人国产av品久久久| 九草在线视频观看| 亚洲怡红院男人天堂| 一本一本综合久久| 国产在线一区二区三区精| 亚洲经典国产精华液单| 五月伊人婷婷丁香| 卡戴珊不雅视频在线播放| 国产精品嫩草影院av在线观看| 久久久久精品久久久久真实原创| 91久久精品国产一区二区三区| 日本黄色片子视频| 精品酒店卫生间| 大码成人一级视频| 欧美3d第一页| 精品国产露脸久久av麻豆| 日韩中字成人| 一区二区三区四区激情视频| 国产色爽女视频免费观看| 两个人的视频大全免费| 视频中文字幕在线观看| 国产v大片淫在线免费观看| 国产成人精品福利久久| 久久精品国产自在天天线| 色视频在线一区二区三区| 嘟嘟电影网在线观看| 久久久久久久久久成人| 亚洲国产精品专区欧美| 久久 成人 亚洲| 成年av动漫网址| 精品酒店卫生间| 一级毛片 在线播放| 精品久久久精品久久久| 亚洲国产毛片av蜜桃av| h视频一区二区三区| 偷拍熟女少妇极品色| av不卡在线播放| 日日啪夜夜爽| 国产av精品麻豆| 99久久中文字幕三级久久日本| 女性被躁到高潮视频| 午夜视频国产福利| 3wmmmm亚洲av在线观看| 免费观看性生交大片5| 成年免费大片在线观看| 亚洲国产日韩一区二区| 亚洲欧美日韩无卡精品| av国产精品久久久久影院| 国产精品伦人一区二区| 一个人看视频在线观看www免费| 亚洲人成网站在线播| 妹子高潮喷水视频| 日韩,欧美,国产一区二区三区| av国产免费在线观看| 亚洲色图综合在线观看| 中文字幕人妻熟人妻熟丝袜美| 亚洲一区二区三区欧美精品| 亚洲国产欧美人成| 午夜福利影视在线免费观看| 午夜免费观看性视频| 国产精品福利在线免费观看| 久久人人爽av亚洲精品天堂 | 伦精品一区二区三区| 久久精品人妻少妇| 欧美日韩视频高清一区二区三区二| 乱码一卡2卡4卡精品| 夜夜看夜夜爽夜夜摸| 欧美少妇被猛烈插入视频| 高清不卡的av网站| 国产 一区精品| 成人18禁高潮啪啪吃奶动态图 | 99热这里只有精品一区| 成人高潮视频无遮挡免费网站| 26uuu在线亚洲综合色| 嫩草影院入口| 中文字幕久久专区| 22中文网久久字幕| 中国国产av一级| 国产精品伦人一区二区| 久久鲁丝午夜福利片| 日本wwww免费看| 国产精品熟女久久久久浪| 日韩一区二区视频免费看| 亚洲欧美日韩无卡精品| 色哟哟·www| 欧美日韩综合久久久久久| 亚洲av成人精品一二三区| av卡一久久| 国产精品麻豆人妻色哟哟久久| 在线看a的网站| 亚洲精品国产成人久久av| 国产精品国产三级国产专区5o| 黄色配什么色好看| 亚洲综合色惰| 国产精品不卡视频一区二区| 成人亚洲精品一区在线观看 | 国内揄拍国产精品人妻在线| 婷婷色综合www| 国产av一区二区精品久久 | xxx大片免费视频| 国产成人91sexporn| 亚洲精品亚洲一区二区| 色视频在线一区二区三区| 亚洲精品乱码久久久v下载方式| 欧美日韩一区二区视频在线观看视频在线| 国产亚洲91精品色在线| 最近的中文字幕免费完整| 女性生殖器流出的白浆| 91久久精品国产一区二区三区| 高清毛片免费看| 亚洲欧美精品专区久久| 在线免费十八禁| 久久国产精品大桥未久av | 亚州av有码| 亚洲精品456在线播放app| 三级国产精品欧美在线观看| 日日摸夜夜添夜夜爱| 国产在线男女| 日本wwww免费看| 国产成人aa在线观看| 岛国毛片在线播放| 日本猛色少妇xxxxx猛交久久| 干丝袜人妻中文字幕| 亚洲欧美成人精品一区二区| 亚洲精品456在线播放app| 国产欧美日韩一区二区三区在线 | 亚州av有码| 麻豆精品久久久久久蜜桃| 性高湖久久久久久久久免费观看| 99热网站在线观看| 欧美丝袜亚洲另类| 男人添女人高潮全过程视频| 一区二区三区精品91| 欧美zozozo另类| 精品一品国产午夜福利视频| 又大又黄又爽视频免费| a 毛片基地| 久久久欧美国产精品| 在线观看一区二区三区激情| 波野结衣二区三区在线| 十八禁网站网址无遮挡 | 亚洲色图av天堂| 在线观看国产h片| 99热这里只有是精品50| 国产午夜精品久久久久久一区二区三区| 王馨瑶露胸无遮挡在线观看| 人人妻人人爽人人添夜夜欢视频 | 99视频精品全部免费 在线| 欧美另类一区| 中文字幕制服av| 日本黄大片高清| 噜噜噜噜噜久久久久久91| 欧美丝袜亚洲另类| 国产无遮挡羞羞视频在线观看| 亚洲国产色片| 国产成人午夜福利电影在线观看| 九九久久精品国产亚洲av麻豆| 久久久久国产精品人妻一区二区| 三级经典国产精品| 少妇人妻久久综合中文| 日韩亚洲欧美综合| 久久久久国产网址| 精品一区在线观看国产| 蜜桃久久精品国产亚洲av| 日韩人妻高清精品专区| 一区二区三区免费毛片| 尾随美女入室| 成人无遮挡网站| 久久久精品94久久精品| 亚洲av不卡在线观看| 91午夜精品亚洲一区二区三区| 国产免费视频播放在线视频| 久久这里有精品视频免费| 日本wwww免费看| 性色av一级| 纵有疾风起免费观看全集完整版| 91久久精品国产一区二区成人| 亚洲不卡免费看| 建设人人有责人人尽责人人享有的 | 女的被弄到高潮叫床怎么办| 看免费成人av毛片| a级毛片免费高清观看在线播放| 亚洲av不卡在线观看| 亚洲av国产av综合av卡| 高清欧美精品videossex| 免费av不卡在线播放| 偷拍熟女少妇极品色| 少妇精品久久久久久久| 国产精品无大码| 99热6这里只有精品| 国产免费福利视频在线观看| 亚洲高清免费不卡视频| 中文资源天堂在线| 黑人高潮一二区| 男人狂女人下面高潮的视频| 天堂中文最新版在线下载| 男女国产视频网站| 欧美日韩亚洲高清精品| 日本黄色日本黄色录像| 日日啪夜夜撸| 18禁在线无遮挡免费观看视频| 免费人成在线观看视频色| 成人二区视频| 午夜福利影视在线免费观看| 国产伦精品一区二区三区四那| 久久久成人免费电影| 蜜桃亚洲精品一区二区三区| 制服丝袜香蕉在线| 欧美zozozo另类| 中国美白少妇内射xxxbb| 成人18禁高潮啪啪吃奶动态图 | 国产精品成人在线| 最近2019中文字幕mv第一页| 国产精品一区二区性色av| 国语对白做爰xxxⅹ性视频网站| 国产又色又爽无遮挡免| 2018国产大陆天天弄谢| 丝袜喷水一区| 亚洲av男天堂| 久久久久久久久久久丰满| 最新中文字幕久久久久| 国产精品久久久久久久久免| 91aial.com中文字幕在线观看| 久久精品国产亚洲网站| 22中文网久久字幕| 春色校园在线视频观看| 人妻制服诱惑在线中文字幕| 日韩中文字幕视频在线看片 | 亚洲综合色惰| 欧美成人一区二区免费高清观看| 久热这里只有精品99| 爱豆传媒免费全集在线观看| 欧美3d第一页| 婷婷色麻豆天堂久久| 国语对白做爰xxxⅹ性视频网站| 黄色欧美视频在线观看| 日韩不卡一区二区三区视频在线| 美女cb高潮喷水在线观看| 91精品国产国语对白视频| 精品一品国产午夜福利视频| 色婷婷av一区二区三区视频| 国产精品伦人一区二区| 亚洲美女黄色视频免费看| 色视频www国产| 亚洲真实伦在线观看| 乱系列少妇在线播放| 三级国产精品欧美在线观看| 久久久色成人| 在线精品无人区一区二区三 | 亚洲精品456在线播放app| 少妇人妻一区二区三区视频| 国产无遮挡羞羞视频在线观看| 国产黄频视频在线观看| 久久99热6这里只有精品| 男女边摸边吃奶| 大香蕉97超碰在线| 我要看黄色一级片免费的| 国产欧美日韩一区二区三区在线 | 国产在视频线精品| 色综合色国产| 高清不卡的av网站| 人人妻人人爽人人添夜夜欢视频 | 国产黄色免费在线视频| 大香蕉97超碰在线| 国产在线视频一区二区| 成人高潮视频无遮挡免费网站| 国产视频首页在线观看| 狂野欧美激情性xxxx在线观看| 嘟嘟电影网在线观看| 免费高清在线观看视频在线观看| 男人爽女人下面视频在线观看| av福利片在线观看| 国产成人免费无遮挡视频| 人妻制服诱惑在线中文字幕| 天美传媒精品一区二区| 一个人看视频在线观看www免费| 国产精品精品国产色婷婷| 欧美97在线视频| 啦啦啦在线观看免费高清www| 涩涩av久久男人的天堂| 亚洲人成网站在线播| 人妻制服诱惑在线中文字幕| 国产精品女同一区二区软件| 小蜜桃在线观看免费完整版高清| 国产成人精品婷婷| 亚洲精品中文字幕在线视频 | 久久青草综合色| 午夜免费观看性视频| 99久久中文字幕三级久久日本| 欧美成人午夜免费资源| 观看av在线不卡| 亚洲成人一二三区av| 欧美日韩视频精品一区| 久久精品久久久久久久性| 尤物成人国产欧美一区二区三区| 内射极品少妇av片p| 在线观看国产h片| 麻豆成人av视频| 国产亚洲91精品色在线| 亚洲伊人久久精品综合| 欧美成人精品欧美一级黄| 日韩伦理黄色片| 国产91av在线免费观看| 国产伦精品一区二区三区视频9| 亚洲精品国产av蜜桃| 九九久久精品国产亚洲av麻豆| 亚洲av在线观看美女高潮| 搡女人真爽免费视频火全软件| 免费观看a级毛片全部| 国产亚洲精品久久久com| 国产av码专区亚洲av| 国产在视频线精品| 欧美三级亚洲精品| 久久国内精品自在自线图片| 网址你懂的国产日韩在线| 18禁在线无遮挡免费观看视频| 波野结衣二区三区在线| 欧美成人一区二区免费高清观看| h日本视频在线播放| 亚洲精品国产av成人精品| 欧美三级亚洲精品| 国产亚洲精品久久久com| 插逼视频在线观看| 精品久久久久久久末码| 高清午夜精品一区二区三区| 99久久综合免费| 波野结衣二区三区在线| 成人18禁高潮啪啪吃奶动态图 | 十八禁网站网址无遮挡 | 午夜福利影视在线免费观看| 国产在线一区二区三区精| 国产成人精品福利久久| 18禁裸乳无遮挡动漫免费视频| 性色av一级| 国产大屁股一区二区在线视频| 大码成人一级视频| 久久这里有精品视频免费| 日韩强制内射视频| 欧美人与善性xxx| 91在线精品国自产拍蜜月| 波野结衣二区三区在线| 久久久久网色| av播播在线观看一区| 精品亚洲乱码少妇综合久久| xxx大片免费视频| 在线播放无遮挡| 国产成人免费无遮挡视频| 欧美xxⅹ黑人| 乱系列少妇在线播放| 亚洲欧美精品专区久久| 免费看日本二区| 美女国产视频在线观看| 一级毛片aaaaaa免费看小| 偷拍熟女少妇极品色| 久久久久久久久久人人人人人人| 亚洲国产最新在线播放| 日韩欧美 国产精品| 欧美精品国产亚洲| 嘟嘟电影网在线观看| 国产成人精品福利久久| 卡戴珊不雅视频在线播放| 成人亚洲欧美一区二区av| 亚洲精品乱码久久久久久按摩| 久久久久精品性色| 亚洲一区二区三区欧美精品| 香蕉精品网在线| 亚洲国产av新网站| 免费少妇av软件| 一级二级三级毛片免费看| 久久鲁丝午夜福利片| 一本久久精品|