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

    交通管制在新型冠狀病毒傳播中的作用

    2020-07-31 09:47:22顧天奇莊楚天李晉梁華包淵秋
    上海城市管理 2020年4期
    關鍵詞:新型冠狀病毒

    顧天奇 莊楚天 李晉 梁華 包淵秋

    摘要:在大范圍疫病傳播時,通過交通管制,政府可以有效控制疫區(qū)與其他區(qū)域的聯系,以期減緩、阻斷疫病的傳播。為研究交通管制措施對控制疫病傳播的作用,論文提出一種基于百度遷徙數據、考慮了人口遷徙效果的改進SIR傳染病傳播動力學模型,進行“有無”交通管制的對比分析。研究發(fā)現,在不同的遷徙人口與本地人口的混合程度下,無交通管制情況下新型冠狀病毒染病者峰值為有交通管制情況下峰值的1.61倍~2.69倍。本研究將城市規(guī)劃與交通規(guī)劃中常用的OD矩陣與經典的傳染病動力學模型綜合考慮,建立了融合模型,具有普適性。

    關鍵詞:新型冠狀病毒;傳染病傳播動力學;SIR模型;交通管控;人口遷徙;OD矩陣

    DOI:10.3969/j.issn.1674-7739.2020.04.009

    引言

    傳染性疾病的傳播受生態(tài)、環(huán)境和人口等多方面的因素影響。氣候、溫度、傳播媒介、當地的醫(yī)療衛(wèi)生條件以及人口的遷徙都可以影響疫病傳播效果。研究表明,人口大量聚集的城市和與外界交流頻繁的城市更容易遭到諸如風疹、麻疹,天花、登革熱的侵襲。[1][2]甚至對人口流動性較好的城市群,一些傳染?。ㄈ缧猴L疹)的流行有可能無法徹底根除。[3]

    很顯然在應對疫病傳播時,是沒法寄希望于氣候、溫度、蟲媒等自然因素的。而短期內提高疫區(qū)醫(yī)療衛(wèi)生條件并補充動員大量的醫(yī)護人員,需要政府部門極大的協(xié)調組織、物資儲備調配能力,實際操作中難度也很大。事實上,只有封鎖疫區(qū)實行交通管制是相對簡單快速、具有可操作性的應對措施。但是大規(guī)模的交通管制、對社會經濟發(fā)展有負面的影響。因此研究、預測人口遷徙在疫病傳播中的作用,具有重要的現實意義,也是傳染病動力學研究的一項重要內容?;?0世紀麻疹在英國傳播的數據,O N Bj?rnstad & B T Grenfell提出了一個風險模型,預測流動性不同的區(qū)域暴發(fā)疫病的可能。[3]C Poletto等研究了交通管制在2014年西非埃博拉病毒傳播中的作用,[4]發(fā)現大規(guī)模的跨國交通管制作用有限,僅能有限延緩病毒在鄰國暴發(fā)的時間。此外值得注意的是,相比于傳統(tǒng)人口流動性的研究方法,借助大數據手段追蹤人口遷徙進而研究疫病傳播,具有樣本大、精度高的優(yōu)點,近年來采用的較多。A Wesolowski運用肯尼亞、納米比亞、巴基斯坦三國的手機信令數據,研究了上述國家大范圍季節(jié)性人口遷徙(節(jié)日、寒暑假)在傳染病傳播動力學中的作用。[5]Wesolowski團隊還利用手機信令數據研究巴基斯坦登革熱的傳播與人口遷徙之間的關系,以預測不同區(qū)域疫病的傳播、集中暴發(fā)風險。[6]在最近《柳葉刀》雜志刊登的關于新型冠狀病毒傳播研究中, J T Wu等學者借助騰訊位置大數據對病毒空間傳播進行了預測。[7]上述研究中,跨區(qū)域的交通,通常被認為是影響傳染系數的一個因素,很少作為獨立參數進入傳染動力學模型進行考察。

    2019年底新型冠狀病毒肺炎(下文簡稱COVID-19)在武漢開始傳播,截至2020年3月10日,中國總計報告累計確診病例80 931例,目前已轉入可防可控態(tài)勢。但是國外累計確診也達到了67 760例。截至論文提交時間(3月15日),日本、韓國、伊朗、意大利等國家感染案例每日都在激增,全球的疫病傳播形勢依舊嚴峻。為應對疫情,1月23日武漢“封城”——實行嚴格的進出交通管制措施,截斷了春節(jié)期間以及春節(jié)后的跨市客流。

    如果不采用如此嚴格的交通管制措施,病毒的跨市傳播可能是怎樣的情況?如何量化交通管制在疫情控制中的作用?這些問題將會是下一階段相關疫病評估工作的重點,也是本研究的立足點?;诮浀涞膫魅静恿W預測模型(SIR模型),本研究建立一個考慮了空間人口遷徙的簡單房室模型,模擬COVID-19在湖北省各城市之間傳播的時空分布(傳播時長、感染人數等)。與常規(guī)傳染病動力學研究不同,本研究重點關注跨區(qū)域的人員遷徙在疫病傳播中的作用,從城市、交通規(guī)劃的視角考察空間管制策略的有效性,探索交通規(guī)劃模型與傳染病傳播動力學模型兩者間融合的可能。

    一、模型選擇及數據準備

    (一)傳染病動力學預測模型(SIR model)

    SIR模型是傳染病動力學語境下的房室模式(Compartmental models)中的經典模型,[8]普遍被用來預測傳染病傳播(如總的傳染人數、持續(xù)時間)。SIR模型將(潛在)傳染人群分為三類:易感者(Susceptibles),其數量記為S(t),表示t時刻未染病但有可能被該類疾病傳染的人數;染病者(Infectives),其數量記為I(t),表示t時刻已被感染成為病人而且具有傳染力的人數;移除者(Removals),其數量記為R(t),表示t時刻已從染病者中移出的人數。假設N(t)且不變,則有N(t)=S(t)+I(t)+R(t)。

    此外,t時刻單位時間內,一個病人能傳染的易感者數目與此環(huán)境內易感者總數S(t)成正比,設定比例系數為β,從而在t時刻單位時間內被所有病人傳染的人數為βS(t)I(t)。通過接觸發(fā)生感染的期望時間Tc=β-1。鐘南山、[9]J M READ[10]等不同研究團隊的共識是:1月23日開始,由于醫(yī)療條件、社區(qū)隔離度的加強,β隨時間降低。t時刻,單位時間內從染病者中移出的人數與病人數量成正比,比例系數為γ,單位時間內移出者的數量為γi(t)。感染后治愈(死亡)的期望時間Tr=γ-1。β/γ為基本傳染率(R0),即在不做管控情況下,一個染病者在一個時間步長(通常為一天)引起凈增加的傳染者人數。早期 J M READ等認為武漢COVID-19的R0為3.11,[10]近期J T Wu發(fā)表于柳葉刀的論文中,R0略有下降,為2.68。[7]

    在上述假設條件下,可知當易染者和染病者個體混合時,感染個體增長率為βIS-γI,易染個體下降率為βIS,恢復個體增長率為γI。對特定的時間t,整個過程可用微分方程表達如下:

    基于上述假設,傳染機制的邏輯是傳染者在一定時間內傳播給潛在的感染者,同時傳染者本身不停地被治愈(或病死),因此同一時間內傳染者凈增加還是凈減少,可以簡單理解為取決于感染率和治愈率。

    作為經典的傳染病動力學模型,SIR模型邏輯清晰,計算簡單,但也存在未考慮人口自然增長、未考慮疾病潛伏期等模型固有限制,因此在此基礎上又衍生出SEIR等一系列改進模型。但是SIR及其衍生模型通常不單獨考慮跨區(qū)出行在傳播中的作用,而是將其作為影響傳染率β的次一級因素考慮。伴隨人口的遷徙,從(疫區(qū))遷徙到當前地區(qū)的人員流動造成的疫病傳染以及疫區(qū)人口遷徙到其他區(qū)域造成的疫病傳染均沒法在模型中直接體現,因此SIR及其衍生模型較難直接與城市規(guī)劃、交通規(guī)劃數據、模型建立直接聯系用于確定行政邊界的疫病傳播模擬。基于此,為考慮空間的人口遷徙,本研究參考O N Bj?rnstad & B T Grenfell2007年的研究成果,在SIR中考慮人口遷徙,建立考慮了人口遷徙的SIR預測模型預測COVID-19在不同交通管制情況下的傳播情況。

    (二)疫病暴發(fā)的概率模型和考慮了人口遷徙的SIR預測模型(Mobility-SIR model)

    O N Bj?rnstad & B T Grenfell在2007年提出一種在不同社區(qū)傳播、并考慮時空因素的疫病暴發(fā)風險模型。[3]該模型認為,對某社區(qū)j來說,其產生第一個傳染者的概率(即疫情暴發(fā)概率)可由兩個變量計算,一是疫區(qū)與j社區(qū)的人群接觸的概率,以及發(fā)生接觸后產生傳染的概率。對于前者,O N Bj?rnstad & B T Grenfell認為,j社區(qū)易染者與外來社區(qū)中染病者產生一次接觸的概率ι服從二項分布:[3]

    此概率與易感人群S,其他社區(qū)的傳染人群比例和作為空間聯系函數正相關。此外應注意到,當易感人群S足夠大時候,風險模型可近似認為只與和相關。也就是說只與其他社區(qū)感染情況和j社區(qū)與其他社區(qū)的交通聯系有關。

    為簡化計算,進一步假設函數即為從其他社區(qū)到j社區(qū)的出行量mt,j。從交通規(guī)劃的角度,注意到mt,j事實上指OD矩陣中流入j社區(qū)的列和。至此,受上述風險模型的啟發(fā),本研究在SIR模型中考慮人口空間遷徙,將公式(1)修改如下:

    值得注意的是,SIR模型求解的一個約束條件是S+I+R之和為常數N,而上式中考慮人口空間遷徙的之和并不是常數。為簡化計算,在各社區(qū)(城市)OD的行和與列和(即流出與流入量分別之和)差別不大,且與S相比較小的情況下,上式可以進一步簡化為:

    上式即為考慮了人口空間遷徙的SIR模型。在這里,有幾個新的前提條件,第一,考慮人口遷徙的模型中,N隨時間變化,但是不同時間段的N(t)=S(t)+I(t)+R(t)仍然成立。此外,改進模型不但假設疫病在某個區(qū)域傳播,還假設跨區(qū)出行染病者也與當地易染者存在感染接觸、發(fā)生傳染。以往的SIR模型,交通管制因素往往作為影響傳染系數β的次一級因素,而考慮人口空間遷徙的SIR模型有直接反映不同地區(qū)交通流動的變量mj,k,并在模型中獨立表達,更符合實際情況。

    值得注意的是,上式中的,分別反映了從j向外圍社區(qū),以及從外圍社區(qū)地流向j地的出行量,在數值上等于空間OD矩陣中j的行和與列和。理論上在經典SIR模型之外,只要有一個按時間變化的OD矩陣(區(qū)一級、市一級、跨省、跨國均可),即可以預測相對空間尺度區(qū)域之間、考慮了人口遷徙的疫病傳播。因OD矩陣是城市(交通)規(guī)劃的模型基礎,至此,傳染病動力學模型與傳統(tǒng)城市(交通)規(guī)劃建立起了空間聯系。而下一步,就是確定上式中的Σmj,k,即武漢流入其他城市的出行量。

    (三)人口遷徙數據

    為簡化情況,論文主要研究武漢前往湖北省其他城市的人口遷徙。假設病毒從武漢開始傳播,且與人員流動有關,需要知道武漢出發(fā)、到達其他城市的出行人數——既包括現實實施交通管控下的出行量,也包括未實施管控情況下的預測出行量。

    基于百度遷徙的數據[12]并參考 J T Wu[7]等人的研究成果,標定后的武漢1月1號至3月7日的實際流出、流入人口數據,以及1月24日起(武漢于1月23日開始采取跨市交通管制)無交通管制狀態(tài)下的流出人口數據。

    根據測算,1月1日開始,流入、流出武漢數量持續(xù)走高。1月20日前,流入大于流出,預計是春運回鄉(xiāng)潮所致。而1月20日后,流出人數突增并超過流入,直至1月23日凌晨出城人數達到最高峰。結合當時的嚴峻情況,這也許與疫情嚴峻導致部分外地人口突擊離開武漢有關。1月23日封城之后,人員流動性被凍結,出行量急劇下降。但是應注意一直到1月28日,武漢的流出量才下降到一個較低的水平(5萬人次/日),也就是說武漢頒布“封城”措施的5天后才真正完成了限制跨市出行這一目標,這樣的狀態(tài)一直保持至今(3月份)。事實上,正常春節(jié)后的流出的春運量應該與節(jié)前流入的春運量基本一致,但是由于交通管制,以1月23日為界,之前大量回武漢的人員在春運后沒有流出。特別是2月1日(正月初八)后,預期的春運返程潮并沒有出現,武漢的流入維持在最高峰約1/50、3~4萬人次/日的水平,流出的更少約為2萬人次/日。考慮到其間武漢出城管控極嚴,但還是存在極少量醫(yī)護資源、保障物流進入的情況,流出略少于流入也與實際一致。

    為預測在無跨市交通管控情況下的人口遷徙,基于歷年春運的特征,本研究假設春節(jié)前后流入和流出武漢的人數基本呈現鏡像分布——以春節(jié)前各市流入武漢的流量作為春節(jié)后流出武漢的流量,并參考上一年度春運人口遷徙情況進行標定。例如,2月1日(初八)武漢流出量,為1月22日(臘月廿九)的流入量?;谏鲜黾僭O,對武漢市人口流出量進行預測,1月23日開始,無交通管制下預測的流出量維持在40萬人次每日的水平,并于2月10左右達到高峰(90萬人次每日),之后穩(wěn)定在30~40萬人次的水平,約為有交通管制下實際流出量的15~20倍。

    此外,模型中關于未管控假想狀態(tài)下的人口遷徙采用了相對簡單的假設情況,以春節(jié)前各市流入武漢的流量作為春節(jié)后流出武漢的流量,可能存在一定誤差,在今后研究中可進一步優(yōu)化。

    參考文獻:

    [1]Wesolowski A, Eagle N, Tatem A J, et al. Quantifying the impact of human mobility on malaria[J]. Science, 2012, 338(6104): 267-270.

    [2]Metcalf C J E, Munayco C, Chowell G, et al. Rubella metapopulation dynamics and importance of spatial coupling to the risk of congenital rubella syndrome in Peru[J]. Journal of the Royal Society Interface, 2011, 8(56): 369-376.

    [3]Bj?rnstad O N, Grenfell B T. Hazards, spatial transmission and timing of outbreaks in epidemic metapopulations[J]. Environmental and Ecological Statistics, 2008, 15(3): 265-277.

    [4]Poletto C, Gomes M F, Y Piontti A P, et al. Assessing the impact of travel restrictions on international spread of the 2014 West African Ebola epidemic[J]. Euro surveillance: bulletin Europeen sur les maladies transmissibles= European communicable disease bulletin, 2014, 19(42).

    [5]Wesolowski A, Zu Erbach-Schoenberg E, Tatem A J, et al. Multinational patterns of seasonal asymmetry in human movement influence infectious disease dynamics[J]. Nature communications, 2017, 8(1): 1-9.

    [6]Wesolowski A, Qureshi T, Boni M F, et al. Impact of human mobility on the emergence of dengue epidemics in Pakistan[J]. Proceedings of the National Academy of Sciences, 2015, 112(38): 11887-11892.

    [7]Wu J T, Leung K, Leung G M. Nowcasting and forecasting the potential domestic and international spread of the 2019-nCoV outbreak originating in Wuhan, China: a modelling study[J]. The Lancet, 2020, 395(10225): 689-697.

    [8]Hethcote H W. The mathematics of infectious diseases[J]. SIAM review, 2000, 42(4): 599-653.

    [9]Yang Z, Zeng Z, Wang K. Modified SEIR and AI prediction of the epidemics trend of COVID-19 in China under public health interventions[J]. J Thorac Dis, 2020.

    [10]Read J M, Bridgen J R, Cummings D A, et al. Novel coronavirus 2019-nCoV: early estimation of epidemiological parameters and epidemic predictions[J]. medRxiv, 2020: 2020.01.23.20018549.

    [11]Bj?rnstad O N, Finkenst?dt B F, Grenfell B T. Dynamics of measles epidemics: estimating scaling of transmission rates using a time series SIR model[J]. Ecological monographs, 2002, 72(2): 169-184.

    [12]百度.百度遷徙大數據,2020. http://qianxi.baidu.com/.

    [13]Guan W-J, Ni Z-Y, Hu Y, et al. Clinical characteristics of 2019 novel coronavirus infection in China[J]. medRxiv, 2020: 2020.02.06.20020974.

    Abstract:Infectious diseases transmission is affected by various factors as weather condition, temperature, medical condition as well as human mobility. Among these factors, human mobility is the most easy-to-manipulate one. During the 2020 coronavirus epidemic, the Chinese government implemented temporal travel restriction, tempted to mitigate the spatial coupling between epidemic area and other areas, to interrupt or even block the disease transmission. To estimate the effect of the travel restriction, an advanced SIR epidemic dynamic model based on human mobility is developed based on Baidu Qianxi data. The model is used to estimate the corona-virus transmission provided that no travel restriction was implemented, so that a before and after change outcome could be investigated and the effect of travel restriction policy could be measured. It is found that peak value of infections without travel restriction is 2.69-1.61 times as that of infections with travel restriction, in terms of different contact situations, and the peak time is also delayed to 8 -11 days when travel restriction is implemented. This paper builds an inclusive model for good incorporating OD matrix frequently seen in urban? and transportation planning and classic infectious disease dynamic model.

    Key words:corona-virus; epidemic dynamics; SIR model; travel restriction; human mobility; OD matrix

    責任編輯:許? 丹

    猜你喜歡
    新型冠狀病毒
    2019新型冠狀病毒SEIR模型求解及MATLAB軟件實現
    淺析新型冠狀病毒疫情時期的家庭健身
    大經貿(2020年2期)2020-05-08 08:43:25
    “新冠肺炎”風暴中的心理障礙及防護研究
    防疫期“燃點”
    新型冠狀病毒暴發(fā)流行對護理教師的啟示
    疫情防控中急需關注青少年的全球素養(yǎng)
    從疫情到輿情
    突發(fā)性公共衛(wèi)生事件:政府與公眾的變與不變
    法制與社會(2020年7期)2020-03-31 03:09:19
    突發(fā)公共衛(wèi)生事件中的藥物使用法規(guī)依據
    疫情防控下的高校網上財務管理思考
    會計之友(2020年6期)2020-03-11 13:18:22
    免费观看无遮挡的男女| 少妇裸体淫交视频免费看高清| 99久久精品国产国产毛片| 中文字幕亚洲精品专区| 不卡视频在线观看欧美| 少妇人妻一区二区三区视频| 精品一区在线观看国产| 久久久久精品性色| 国产成人a区在线观看| 国产欧美亚洲国产| 少妇被粗大猛烈的视频| 综合色丁香网| 干丝袜人妻中文字幕| 国产精品人妻久久久久久| 中文天堂在线官网| 中文精品一卡2卡3卡4更新| 国模一区二区三区四区视频| videos熟女内射| 国产在线视频一区二区| 日日啪夜夜撸| 久久久久人妻精品一区果冻| 亚洲国产精品成人久久小说| 久久婷婷青草| 人体艺术视频欧美日本| 在线观看一区二区三区| 午夜福利视频精品| 免费黄网站久久成人精品| 亚洲国产日韩一区二区| 一区二区av电影网| 一本久久精品| 免费少妇av软件| 大陆偷拍与自拍| 精品国产乱码久久久久久小说| 人人妻人人添人人爽欧美一区卜 | 亚洲第一av免费看| 日韩人妻高清精品专区| 国产精品欧美亚洲77777| 在线 av 中文字幕| 久久精品国产亚洲av天美| 久久久久国产网址| 草草在线视频免费看| 国产免费视频播放在线视频| 精品人妻一区二区三区麻豆| h日本视频在线播放| 少妇丰满av| 久久 成人 亚洲| 亚洲精品中文字幕在线视频 | 夫妻性生交免费视频一级片| 最近最新中文字幕免费大全7| 18禁在线无遮挡免费观看视频| 欧美极品一区二区三区四区| 免费大片黄手机在线观看| 中文欧美无线码| 一级二级三级毛片免费看| 亚洲欧美日韩另类电影网站 | 爱豆传媒免费全集在线观看| 久久精品夜色国产| 伦精品一区二区三区| 久久99精品国语久久久| 久久国产亚洲av麻豆专区| 中文字幕亚洲精品专区| 插逼视频在线观看| 国产国拍精品亚洲av在线观看| 大香蕉97超碰在线| 人人妻人人看人人澡| 日韩一本色道免费dvd| 亚洲国产欧美人成| 妹子高潮喷水视频| av福利片在线观看| 亚洲自偷自拍三级| 97精品久久久久久久久久精品| 成年女人在线观看亚洲视频| 在线观看免费日韩欧美大片 | 一级a做视频免费观看| 久久久久精品性色| 内射极品少妇av片p| 欧美亚洲 丝袜 人妻 在线| 丰满迷人的少妇在线观看| 在线免费观看不下载黄p国产| 国产成人午夜福利电影在线观看| 亚洲天堂av无毛| 久久人人爽av亚洲精品天堂 | 十分钟在线观看高清视频www | 麻豆乱淫一区二区| 岛国毛片在线播放| 亚洲av福利一区| 在线 av 中文字幕| 中文资源天堂在线| 亚洲欧美成人综合另类久久久| 国产成人a∨麻豆精品| 欧美高清性xxxxhd video| 亚洲经典国产精华液单| 又粗又硬又长又爽又黄的视频| 精品一区二区三区视频在线| 亚洲性久久影院| 美女主播在线视频| 久久精品国产自在天天线| 欧美另类一区| 妹子高潮喷水视频| 精品国产露脸久久av麻豆| 亚洲自偷自拍三级| 欧美xxxx性猛交bbbb| 女的被弄到高潮叫床怎么办| 国产免费又黄又爽又色| 看免费成人av毛片| 亚洲电影在线观看av| 免费人妻精品一区二区三区视频| 国产精品秋霞免费鲁丝片| a级毛片免费高清观看在线播放| 欧美成人a在线观看| 亚洲人与动物交配视频| 老熟女久久久| 欧美老熟妇乱子伦牲交| 日本黄大片高清| 男女下面进入的视频免费午夜| 少妇的逼好多水| 91久久精品国产一区二区成人| 狂野欧美白嫩少妇大欣赏| 国产色爽女视频免费观看| 三级经典国产精品| 午夜日本视频在线| xxx大片免费视频| 观看美女的网站| 成年人午夜在线观看视频| 99热这里只有精品一区| 在线观看免费视频网站a站| 亚洲中文av在线| 精品人妻偷拍中文字幕| 久久久成人免费电影| 国产精品精品国产色婷婷| 18禁裸乳无遮挡免费网站照片| 婷婷色av中文字幕| 国产av一区二区精品久久 | 高清日韩中文字幕在线| 久久 成人 亚洲| 欧美日韩亚洲高清精品| 两个人的视频大全免费| 亚洲精品国产成人久久av| 香蕉精品网在线| 老司机影院毛片| 777米奇影视久久| 韩国av在线不卡| 日韩国内少妇激情av| 国产伦精品一区二区三区视频9| videossex国产| 成人美女网站在线观看视频| 美女主播在线视频| 国产在视频线精品| 久久久久久久久久成人| 亚州av有码| 久久久久国产精品人妻一区二区| 人人妻人人澡人人爽人人夜夜| 777米奇影视久久| 麻豆精品久久久久久蜜桃| 国产在线视频一区二区| 国产亚洲午夜精品一区二区久久| 国产91av在线免费观看| 午夜福利高清视频| 美女福利国产在线 | 观看av在线不卡| 亚洲国产最新在线播放| 午夜免费男女啪啪视频观看| 国产 一区精品| 国产免费一区二区三区四区乱码| 中国美白少妇内射xxxbb| 国产伦精品一区二区三区视频9| 久久久久网色| 日本黄色日本黄色录像| 2018国产大陆天天弄谢| 狠狠精品人妻久久久久久综合| 你懂的网址亚洲精品在线观看| 国产精品免费大片| 欧美极品一区二区三区四区| 亚洲av国产av综合av卡| 99久久精品国产国产毛片| kizo精华| 精品午夜福利在线看| 久久人人爽人人片av| 大香蕉久久网| 日韩强制内射视频| 国产大屁股一区二区在线视频| 日韩av不卡免费在线播放| 色视频在线一区二区三区| 国产亚洲av片在线观看秒播厂| 久久精品人妻少妇| 性色avwww在线观看| 小蜜桃在线观看免费完整版高清| 十分钟在线观看高清视频www | 老熟女久久久| 亚洲美女黄色视频免费看| av国产久精品久网站免费入址| 免费人成在线观看视频色| 亚洲天堂av无毛| 亚洲av免费高清在线观看| 91狼人影院| 国产成人免费观看mmmm| 国产精品久久久久久久久免| 久久久久久九九精品二区国产| 天美传媒精品一区二区| 秋霞伦理黄片| 久久久久久久久大av| 欧美xxxx性猛交bbbb| 一个人看的www免费观看视频| 精品久久久久久久久av| 直男gayav资源| 男女边摸边吃奶| 精品久久久久久久末码| 欧美激情国产日韩精品一区| 亚洲国产欧美在线一区| 九九爱精品视频在线观看| 国产爱豆传媒在线观看| 一区二区三区四区激情视频| 秋霞伦理黄片| 亚洲国产精品专区欧美| 男人狂女人下面高潮的视频| 亚洲av欧美aⅴ国产| 国产av国产精品国产| 亚洲第一av免费看| 色婷婷av一区二区三区视频| 婷婷色av中文字幕| 亚洲av中文av极速乱| 欧美成人一区二区免费高清观看| 一个人看视频在线观看www免费| 精品午夜福利在线看| 久热久热在线精品观看| 日韩av在线免费看完整版不卡| 男男h啪啪无遮挡| 久久久欧美国产精品| 免费播放大片免费观看视频在线观看| 久久久久性生活片| 成人毛片60女人毛片免费| 日韩成人av中文字幕在线观看| 欧美一级a爱片免费观看看| 欧美日韩精品成人综合77777| 老司机影院成人| 亚洲中文av在线| 男女国产视频网站| 久久精品熟女亚洲av麻豆精品| 99热国产这里只有精品6| a 毛片基地| 亚洲av二区三区四区| 亚洲国产欧美人成| 高清在线视频一区二区三区| 国产黄频视频在线观看| av女优亚洲男人天堂| 国产成人精品福利久久| 免费大片18禁| 欧美变态另类bdsm刘玥| xxx大片免费视频| 99热这里只有精品一区| 哪个播放器可以免费观看大片| 国产淫片久久久久久久久| 三级经典国产精品| 国产高清有码在线观看视频| 久久av网站| 亚洲av福利一区| 熟妇人妻不卡中文字幕| av国产免费在线观看| 亚洲国产色片| 欧美日韩一区二区视频在线观看视频在线| 嫩草影院新地址| 久热这里只有精品99| 少妇人妻久久综合中文| 内地一区二区视频在线| 亚洲成人中文字幕在线播放| 最近的中文字幕免费完整| 欧美日韩一区二区视频在线观看视频在线| 国产av国产精品国产| 午夜免费男女啪啪视频观看| 有码 亚洲区| 亚洲国产精品国产精品| 少妇人妻 视频| 国产中年淑女户外野战色| 精品熟女少妇av免费看| 欧美xxxx黑人xx丫x性爽| 国产成人一区二区在线| 日韩强制内射视频| 小蜜桃在线观看免费完整版高清| 麻豆国产97在线/欧美| 亚洲av在线观看美女高潮| 国产午夜精品久久久久久一区二区三区| 中文乱码字字幕精品一区二区三区| 免费av不卡在线播放| 亚洲精品国产成人久久av| 国产日韩欧美在线精品| 日韩大片免费观看网站| 欧美成人午夜免费资源| 欧美bdsm另类| 晚上一个人看的免费电影| 久久精品人妻少妇| 狂野欧美白嫩少妇大欣赏| 永久网站在线| 精品久久久久久久末码| 久久女婷五月综合色啪小说| av国产精品久久久久影院| 国产又色又爽无遮挡免| 亚洲欧洲日产国产| h视频一区二区三区| 免费不卡的大黄色大毛片视频在线观看| 亚洲精品国产av蜜桃| 亚洲一区二区三区欧美精品| 久久久久国产精品人妻一区二区| 精品人妻偷拍中文字幕| 国产美女午夜福利| 在线观看av片永久免费下载| 国产片特级美女逼逼视频| 精品国产露脸久久av麻豆| 卡戴珊不雅视频在线播放| 中文资源天堂在线| av一本久久久久| 美女中出高潮动态图| 亚洲一级一片aⅴ在线观看| 老熟女久久久| 菩萨蛮人人尽说江南好唐韦庄| 妹子高潮喷水视频| 只有这里有精品99| 人妻系列 视频| 欧美精品一区二区大全| 乱系列少妇在线播放| 肉色欧美久久久久久久蜜桃| 中文字幕精品免费在线观看视频 | 男人爽女人下面视频在线观看| 高清午夜精品一区二区三区| 最后的刺客免费高清国语| 色综合色国产| 亚洲aⅴ乱码一区二区在线播放| 亚洲无线观看免费| 91久久精品国产一区二区三区| 91精品国产九色| 我的老师免费观看完整版| 国精品久久久久久国模美| 乱码一卡2卡4卡精品| 日本与韩国留学比较| 国产黄色视频一区二区在线观看| 午夜福利网站1000一区二区三区| 丝袜脚勾引网站| 亚洲欧美精品自产自拍| 在现免费观看毛片| 国产av码专区亚洲av| 亚洲婷婷狠狠爱综合网| 色网站视频免费| 亚洲婷婷狠狠爱综合网| 人体艺术视频欧美日本| 国产亚洲精品久久久com| 亚洲成人av在线免费| 国产成人a∨麻豆精品| 黄色配什么色好看| 男男h啪啪无遮挡| 女人久久www免费人成看片| 久久99热这里只频精品6学生| 亚洲国产精品999| 伦理电影免费视频| 91狼人影院| 久久人人爽av亚洲精品天堂 | av国产久精品久网站免费入址| 看十八女毛片水多多多| 欧美日韩国产mv在线观看视频 | 国产美女午夜福利| 久久热精品热| a级毛色黄片| av网站免费在线观看视频| 男人爽女人下面视频在线观看| 亚洲欧美一区二区三区黑人 | 亚洲丝袜综合中文字幕| 免费久久久久久久精品成人欧美视频 | 高清视频免费观看一区二区| 青青草视频在线视频观看| 国产伦理片在线播放av一区| 人人妻人人添人人爽欧美一区卜 | 97在线视频观看| 久久久精品免费免费高清| 97在线视频观看| 亚洲人成网站在线观看播放| 搡老乐熟女国产| 欧美zozozo另类| 中国美白少妇内射xxxbb| 一级二级三级毛片免费看| 久久亚洲国产成人精品v| 成人漫画全彩无遮挡| 蜜桃在线观看..| 久久影院123| 91狼人影院| 日韩成人av中文字幕在线观看| 国产探花极品一区二区| 久久 成人 亚洲| 亚洲av二区三区四区| 欧美成人一区二区免费高清观看| 天堂中文最新版在线下载| 久久国内精品自在自线图片| 午夜福利网站1000一区二区三区| 在线免费观看不下载黄p国产| 午夜福利在线在线| 蜜桃亚洲精品一区二区三区| 日韩不卡一区二区三区视频在线| 又爽又黄a免费视频| 美女国产视频在线观看| 精品久久久噜噜| 国产男女内射视频| 人妻 亚洲 视频| 狂野欧美激情性xxxx在线观看| 午夜视频国产福利| 观看美女的网站| 国产欧美亚洲国产| 丰满少妇做爰视频| 在线看a的网站| 简卡轻食公司| 免费高清在线观看视频在线观看| 亚洲精品456在线播放app| 3wmmmm亚洲av在线观看| 新久久久久国产一级毛片| 一级毛片 在线播放| 亚洲第一av免费看| 日韩强制内射视频| 亚洲欧美日韩东京热| 国产成人a∨麻豆精品| 99精国产麻豆久久婷婷| 国产片特级美女逼逼视频| 日韩人妻高清精品专区| 建设人人有责人人尽责人人享有的 | 有码 亚洲区| 国产又色又爽无遮挡免| 免费观看性生交大片5| 国产精品.久久久| 啦啦啦视频在线资源免费观看| 成人免费观看视频高清| 最后的刺客免费高清国语| 大话2 男鬼变身卡| 亚洲欧美日韩另类电影网站 | 国产黄色视频一区二区在线观看| 中文欧美无线码| 中文字幕免费在线视频6| 亚洲av.av天堂| 欧美97在线视频| 国产亚洲最大av| 高清视频免费观看一区二区| 亚洲高清免费不卡视频| 亚洲欧美精品自产自拍| 久久人人爽人人片av| 亚洲天堂av无毛| 亚洲av在线观看美女高潮| 日韩三级伦理在线观看| 欧美高清成人免费视频www| 欧美日韩视频高清一区二区三区二| 欧美激情极品国产一区二区三区 | 日韩av不卡免费在线播放| 色视频www国产| 日韩精品有码人妻一区| 日本欧美视频一区| 在线观看国产h片| www.av在线官网国产| 一级毛片我不卡| 1000部很黄的大片| 人妻系列 视频| 久久精品久久久久久噜噜老黄| 国产 一区 欧美 日韩| 亚洲av福利一区| 亚洲精品一二三| 久久6这里有精品| 熟女电影av网| 亚洲欧美日韩另类电影网站 | 99热这里只有是精品在线观看| 中文欧美无线码| 亚洲图色成人| 有码 亚洲区| 欧美人与善性xxx| 国产免费视频播放在线视频| 三级国产精品欧美在线观看| 美女高潮的动态| 日韩精品有码人妻一区| 啦啦啦在线观看免费高清www| 久久久a久久爽久久v久久| 欧美国产精品一级二级三级 | 日韩精品有码人妻一区| 久久国产精品男人的天堂亚洲 | 亚洲av二区三区四区| 老司机影院成人| 自拍偷自拍亚洲精品老妇| 欧美极品一区二区三区四区| 亚洲国产欧美在线一区| 伦精品一区二区三区| av国产久精品久网站免费入址| 中文天堂在线官网| 特大巨黑吊av在线直播| 中国三级夫妇交换| 18+在线观看网站| 91久久精品国产一区二区成人| 日本-黄色视频高清免费观看| 少妇 在线观看| 日本色播在线视频| 亚洲国产毛片av蜜桃av| 欧美一区二区亚洲| 国产av码专区亚洲av| 午夜激情福利司机影院| 日韩不卡一区二区三区视频在线| 一级毛片aaaaaa免费看小| 天美传媒精品一区二区| 亚洲电影在线观看av| 久久精品国产亚洲网站| 搡女人真爽免费视频火全软件| 国产成人a区在线观看| 成人漫画全彩无遮挡| av福利片在线观看| 精品人妻熟女av久视频| 国产一区二区在线观看日韩| 黄色欧美视频在线观看| 久久久久久久精品精品| 免费观看av网站的网址| 干丝袜人妻中文字幕| 我的女老师完整版在线观看| 另类亚洲欧美激情| 国产免费一区二区三区四区乱码| 男人舔奶头视频| 中文字幕免费在线视频6| 国产亚洲91精品色在线| 日韩,欧美,国产一区二区三区| 国产亚洲精品久久久com| 国产黄频视频在线观看| 最近最新中文字幕大全电影3| 久久久久人妻精品一区果冻| 亚洲成人av在线免费| 少妇 在线观看| 国产精品秋霞免费鲁丝片| a级一级毛片免费在线观看| 高清不卡的av网站| 成年美女黄网站色视频大全免费 | 51国产日韩欧美| 国产一区二区在线观看日韩| 成人一区二区视频在线观看| 啦啦啦视频在线资源免费观看| 午夜福利网站1000一区二区三区| 中文字幕制服av| 国产男女内射视频| 成人免费观看视频高清| 欧美xxxx性猛交bbbb| 亚洲av成人精品一二三区| 80岁老熟妇乱子伦牲交| 国产免费视频播放在线视频| 日韩一本色道免费dvd| 国产91av在线免费观看| 97在线视频观看| 亚洲av二区三区四区| 少妇人妻精品综合一区二区| 少妇人妻久久综合中文| 中国美白少妇内射xxxbb| 一本久久精品| 国产一区二区三区av在线| 国产 一区 欧美 日韩| 18禁在线播放成人免费| 国产av码专区亚洲av| 色视频www国产| 97超视频在线观看视频| 欧美+日韩+精品| 国产精品嫩草影院av在线观看| 美女主播在线视频| 国产伦精品一区二区三区四那| av卡一久久| 久久久午夜欧美精品| 在线观看人妻少妇| 精品少妇久久久久久888优播| 在线看a的网站| 狠狠精品人妻久久久久久综合| 亚洲国产欧美人成| 国产高清不卡午夜福利| 韩国av在线不卡| 国产精品国产三级国产av玫瑰| 麻豆成人av视频| 久久国内精品自在自线图片| 欧美日韩精品成人综合77777| 成人一区二区视频在线观看| 嫩草影院入口| 免费人妻精品一区二区三区视频| 亚洲av成人精品一区久久| 欧美日韩视频高清一区二区三区二| 大话2 男鬼变身卡| 精品99又大又爽又粗少妇毛片| 精品一区在线观看国产| 九九在线视频观看精品| 性色av一级| 大片电影免费在线观看免费| www.av在线官网国产| 狂野欧美白嫩少妇大欣赏| 日日撸夜夜添| 免费看光身美女| 麻豆成人午夜福利视频| 国产伦精品一区二区三区视频9| 中国美白少妇内射xxxbb| 国内精品宾馆在线| 美女主播在线视频| 午夜免费男女啪啪视频观看| 春色校园在线视频观看| 黄片wwwwww| 赤兔流量卡办理| 婷婷色综合www| 国产一区二区在线观看日韩| 亚洲人成网站在线观看播放| 美女cb高潮喷水在线观看| 少妇被粗大猛烈的视频| av视频免费观看在线观看| 国模一区二区三区四区视频| 国产一区二区三区av在线| 在线观看免费视频网站a站| 能在线免费看毛片的网站| 一区在线观看完整版| 欧美变态另类bdsm刘玥| 美女cb高潮喷水在线观看| 国产一区二区在线观看日韩| 国产又色又爽无遮挡免| 一本一本综合久久| 日本wwww免费看| 内射极品少妇av片p| 最近2019中文字幕mv第一页| 99热网站在线观看| 人妻制服诱惑在线中文字幕| 狂野欧美激情性bbbbbb| 国产色爽女视频免费观看|