• <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
    亚洲精品久久成人aⅴ小说| 亚洲精品久久午夜乱码| 国产男女超爽视频在线观看| 日韩一区二区视频免费看| av在线app专区| 婷婷色综合www| 91精品国产国语对白视频| 国产亚洲一区二区精品| 国产激情久久老熟女| 日韩人妻精品一区2区三区| 香蕉丝袜av| 在线天堂最新版资源| 午夜免费观看性视频| 国产片特级美女逼逼视频| 日本猛色少妇xxxxx猛交久久| 桃花免费在线播放| 一级毛片黄色毛片免费观看视频| 性色avwww在线观看| 久久人人爽人人爽人人片va| 日韩视频在线欧美| 亚洲精品色激情综合| 成年美女黄网站色视频大全免费| 美女中出高潮动态图| 人妻少妇偷人精品九色| 校园人妻丝袜中文字幕| 搡女人真爽免费视频火全软件| 久久久久久久久久久久大奶| 十八禁网站网址无遮挡| 精品人妻在线不人妻| 日韩一区二区视频免费看| 国产无遮挡羞羞视频在线观看| 国产精品久久久久久av不卡| 精品第一国产精品| av天堂久久9| 一级片免费观看大全| 亚洲人成77777在线视频| 大香蕉久久网| 欧美激情国产日韩精品一区| 亚洲丝袜综合中文字幕| 人妻少妇偷人精品九色| 国国产精品蜜臀av免费| 亚洲精品美女久久av网站| 自线自在国产av| 国产精品一二三区在线看| 国产精品久久久久久精品电影小说| 男男h啪啪无遮挡| 欧美少妇被猛烈插入视频| 亚洲欧洲日产国产| 亚洲国产精品专区欧美| 久久免费观看电影| 国产伦理片在线播放av一区| 久久狼人影院| 巨乳人妻的诱惑在线观看| 乱人伦中国视频| 在线亚洲精品国产二区图片欧美| 水蜜桃什么品种好| 精品久久久精品久久久| 一区二区三区精品91| 中文精品一卡2卡3卡4更新| 精品一区二区免费观看| 日韩制服骚丝袜av| 国产免费福利视频在线观看| 亚洲在久久综合| 亚洲色图综合在线观看| 亚洲情色 制服丝袜| xxx大片免费视频| 精品久久久久久电影网| 99热国产这里只有精品6| 国产国语露脸激情在线看| 欧美日韩成人在线一区二区| 伦理电影大哥的女人| 少妇被粗大的猛进出69影院 | 久久青草综合色| 亚洲天堂av无毛| 欧美精品一区二区大全| 国产色婷婷99| 国产淫语在线视频| 亚洲精品456在线播放app| 国产精品熟女久久久久浪| 国产免费视频播放在线视频| 欧美精品av麻豆av| 成人国产av品久久久| 亚洲成人手机| 香蕉国产在线看| 亚洲av电影在线进入| 久久久a久久爽久久v久久| 久久久久久久大尺度免费视频| 国产精品99久久99久久久不卡 | 久久久a久久爽久久v久久| 大香蕉97超碰在线| 国产精品一区二区在线不卡| 成人国产av品久久久| 久久久欧美国产精品| av视频免费观看在线观看| 18+在线观看网站| 一级爰片在线观看| 精品人妻一区二区三区麻豆| 18禁动态无遮挡网站| 国产乱来视频区| 免费人妻精品一区二区三区视频| 国精品久久久久久国模美| 老女人水多毛片| 国产精品久久久久久精品古装| 亚洲国产欧美日韩在线播放| 免费黄色在线免费观看| 99九九在线精品视频| 一本—道久久a久久精品蜜桃钙片| 日韩制服骚丝袜av| 亚洲av国产av综合av卡| 性高湖久久久久久久久免费观看| 欧美激情 高清一区二区三区| 另类亚洲欧美激情| 成年人免费黄色播放视频| 中文字幕另类日韩欧美亚洲嫩草| 纵有疾风起免费观看全集完整版| 亚洲av日韩在线播放| 十八禁网站网址无遮挡| 波野结衣二区三区在线| 午夜福利在线观看免费完整高清在| 一级毛片黄色毛片免费观看视频| 精品少妇黑人巨大在线播放| 国产色爽女视频免费观看| 亚洲av综合色区一区| 大香蕉久久成人网| 国产免费一级a男人的天堂| 大话2 男鬼变身卡| 大陆偷拍与自拍| 亚洲精品一二三| av在线老鸭窝| 久久精品国产鲁丝片午夜精品| 黄色 视频免费看| 黄网站色视频无遮挡免费观看| 日韩 亚洲 欧美在线| 春色校园在线视频观看| 国国产精品蜜臀av免费| 亚洲内射少妇av| 欧美人与性动交α欧美软件 | 最近2019中文字幕mv第一页| 9热在线视频观看99| av不卡在线播放| 免费黄色在线免费观看| 99久国产av精品国产电影| 国产老妇伦熟女老妇高清| av不卡在线播放| 一边摸一边做爽爽视频免费| 秋霞在线观看毛片| 高清毛片免费看| 80岁老熟妇乱子伦牲交| 欧美+日韩+精品| 九草在线视频观看| 国产毛片在线视频| 亚洲精品美女久久av网站| 亚洲精品日本国产第一区| 大片免费播放器 马上看| videossex国产| 亚洲欧美日韩卡通动漫| 美女内射精品一级片tv| 丰满少妇做爰视频| 18禁在线无遮挡免费观看视频| 日日摸夜夜添夜夜爱| 性色avwww在线观看| 人妻少妇偷人精品九色| 亚洲精品成人av观看孕妇| 看免费成人av毛片| 丝袜美足系列| 人妻人人澡人人爽人人| 咕卡用的链子| 国产精品嫩草影院av在线观看| 精品卡一卡二卡四卡免费| 18+在线观看网站| 久久毛片免费看一区二区三区| 久久女婷五月综合色啪小说| 中文天堂在线官网| 汤姆久久久久久久影院中文字幕| 午夜av观看不卡| 欧美日韩精品成人综合77777| 成人国产av品久久久| 久久这里有精品视频免费| av线在线观看网站| 精品一区二区三卡| 九草在线视频观看| 国产精品一二三区在线看| 免费观看性生交大片5| 中文字幕精品免费在线观看视频 | 亚洲精品久久成人aⅴ小说| 欧美成人午夜免费资源| 一级爰片在线观看| 高清av免费在线| 亚洲国产毛片av蜜桃av| av在线app专区| h视频一区二区三区| 久久久久久久亚洲中文字幕| 内地一区二区视频在线| 超碰97精品在线观看| 99热这里只有是精品在线观看| 18禁在线无遮挡免费观看视频| 男女边吃奶边做爰视频| 亚洲欧洲精品一区二区精品久久久 | 成人国产麻豆网| 亚洲av男天堂| 多毛熟女@视频| 欧美少妇被猛烈插入视频| 欧美另类一区| 一边亲一边摸免费视频| 久久精品久久精品一区二区三区| 最近手机中文字幕大全| 久久国内精品自在自线图片| 最近中文字幕高清免费大全6| 在线观看国产h片| 最近最新中文字幕大全免费视频 | 女人久久www免费人成看片| 欧美日韩综合久久久久久| 亚洲一区二区三区欧美精品| 哪个播放器可以免费观看大片| 精品午夜福利在线看| 最近中文字幕高清免费大全6| 少妇猛男粗大的猛烈进出视频| 91精品伊人久久大香线蕉| 五月天丁香电影| 久久久久久久国产电影| 国产又爽黄色视频| 亚洲,欧美精品.| 你懂的网址亚洲精品在线观看| 日韩人妻精品一区2区三区| 亚洲美女视频黄频| 国产精品三级大全| 亚洲国产最新在线播放| a级毛色黄片| 18禁动态无遮挡网站| 国产日韩欧美在线精品| 十八禁网站网址无遮挡| 1024视频免费在线观看| 亚洲综合色网址| 男女无遮挡免费网站观看| 99热网站在线观看| 视频中文字幕在线观看| 国产片特级美女逼逼视频| 久久精品国产亚洲av天美| 在线观看一区二区三区激情| kizo精华| 国产精品久久久久久久电影| 2022亚洲国产成人精品| 少妇被粗大猛烈的视频| 少妇的丰满在线观看| 欧美性感艳星| 男人爽女人下面视频在线观看| 最近中文字幕高清免费大全6| 国产视频首页在线观看| 午夜福利视频精品| 久久久精品区二区三区| 亚洲av.av天堂| av又黄又爽大尺度在线免费看| 制服人妻中文乱码| 久久精品国产综合久久久 | 欧美激情 高清一区二区三区| 免费av不卡在线播放| 丝袜人妻中文字幕| 69精品国产乱码久久久| 亚洲综合色网址| 免费观看性生交大片5| 在现免费观看毛片| 亚洲av综合色区一区| 一级爰片在线观看| 日本黄色日本黄色录像| 久久人妻熟女aⅴ| 免费人妻精品一区二区三区视频| 亚洲av.av天堂| 两个人看的免费小视频| 日韩免费高清中文字幕av| 秋霞在线观看毛片| 男女午夜视频在线观看 | 中文字幕人妻丝袜制服| 一区二区三区精品91| 欧美精品一区二区免费开放| 国产黄频视频在线观看| 欧美精品一区二区大全| 国产黄色免费在线视频| 啦啦啦在线观看免费高清www| 青春草亚洲视频在线观看| 久久精品国产综合久久久 | 视频在线观看一区二区三区| 日韩中文字幕视频在线看片| 亚洲精品av麻豆狂野| 久久久久精品性色| 精品99又大又爽又粗少妇毛片| 高清毛片免费看| 免费看不卡的av| 寂寞人妻少妇视频99o| 亚洲色图 男人天堂 中文字幕 | 日本av手机在线免费观看| av电影中文网址| 国产亚洲精品第一综合不卡 | 国产一区二区在线观看av| 欧美精品av麻豆av| av卡一久久| 欧美成人午夜精品| 男女午夜视频在线观看 | 国产精品欧美亚洲77777| 婷婷色综合www| 高清毛片免费看| www.av在线官网国产| 天天躁夜夜躁狠狠久久av| 九色成人免费人妻av| 日本欧美视频一区| 精品久久久精品久久久| 韩国高清视频一区二区三区| 国产免费视频播放在线视频| 国产精品一二三区在线看| 久久久欧美国产精品| 这个男人来自地球电影免费观看 | 下体分泌物呈黄色| videos熟女内射| 国产男女超爽视频在线观看| 久久这里只有精品19| 咕卡用的链子| 涩涩av久久男人的天堂| 欧美日韩一区二区视频在线观看视频在线| 日本色播在线视频| 久久ye,这里只有精品| 久久99热6这里只有精品| 80岁老熟妇乱子伦牲交| 欧美xxⅹ黑人| 国国产精品蜜臀av免费| 青春草国产在线视频| xxxhd国产人妻xxx| 老司机影院成人| 免费少妇av软件| 男女下面插进去视频免费观看 | 精品99又大又爽又粗少妇毛片| 五月开心婷婷网| 亚洲高清免费不卡视频| 最新中文字幕久久久久| 免费观看性生交大片5| 九色亚洲精品在线播放| 久久久久久久大尺度免费视频| 欧美激情 高清一区二区三区| 日日撸夜夜添| 视频中文字幕在线观看| 亚洲精品第二区| 亚洲精品久久成人aⅴ小说| 成年美女黄网站色视频大全免费| av片东京热男人的天堂| 毛片一级片免费看久久久久| 国产爽快片一区二区三区| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 91在线精品国自产拍蜜月| 性色av一级| 大香蕉久久网| 如日韩欧美国产精品一区二区三区| 色视频在线一区二区三区| 国产69精品久久久久777片| 亚洲成av片中文字幕在线观看 | 国产午夜精品一二区理论片| 中文字幕最新亚洲高清| 精品一区二区三区四区五区乱码 | 看免费成人av毛片| 免费观看a级毛片全部| 制服人妻中文乱码| 少妇熟女欧美另类| 国产av一区二区精品久久| 久久免费观看电影| 亚洲精品乱久久久久久| 美女国产高潮福利片在线看| 午夜免费男女啪啪视频观看| 汤姆久久久久久久影院中文字幕| 亚洲av电影在线观看一区二区三区| 亚洲av综合色区一区| 久久久久精品人妻al黑| 波多野结衣一区麻豆| 久久 成人 亚洲| 丝袜喷水一区| 99国产综合亚洲精品| 爱豆传媒免费全集在线观看| 国产老妇伦熟女老妇高清| 尾随美女入室| 亚洲av中文av极速乱| 黄色怎么调成土黄色| 亚洲精品一区蜜桃| 日本黄大片高清| 国产精品欧美亚洲77777| 欧美xxxx性猛交bbbb| 亚洲,欧美,日韩| www.色视频.com| a 毛片基地| 下体分泌物呈黄色| 久久国产精品男人的天堂亚洲 | 精品酒店卫生间| 国语对白做爰xxxⅹ性视频网站| 午夜久久久在线观看| 街头女战士在线观看网站| 久久人人爽人人爽人人片va| www日本在线高清视频| 精品亚洲乱码少妇综合久久| 久久久国产精品麻豆| 成人国产麻豆网| 18禁在线无遮挡免费观看视频| 久久婷婷青草| 亚洲国产欧美在线一区| 乱码一卡2卡4卡精品| 亚洲美女视频黄频| 老女人水多毛片| 一本—道久久a久久精品蜜桃钙片| 国产爽快片一区二区三区| 色视频在线一区二区三区| 丁香六月天网| www.色视频.com| 亚洲国产欧美日韩在线播放| 久久99精品国语久久久| 亚洲国产色片| 色网站视频免费| 成人漫画全彩无遮挡| 免费看不卡的av| 亚洲一级一片aⅴ在线观看| 侵犯人妻中文字幕一二三四区| 熟女电影av网| 美女主播在线视频| 99re6热这里在线精品视频| 美女脱内裤让男人舔精品视频| 观看av在线不卡| 18禁在线无遮挡免费观看视频| 欧美激情极品国产一区二区三区 | videosex国产| 亚洲婷婷狠狠爱综合网| 亚洲成色77777| 日日摸夜夜添夜夜爱| 亚洲成av片中文字幕在线观看 | 亚洲成人手机| 亚洲成色77777| 精品少妇内射三级| 亚洲成人av在线免费| 成人无遮挡网站| 国产成人精品在线电影| 日韩免费高清中文字幕av| 你懂的网址亚洲精品在线观看| 国产精品熟女久久久久浪| 欧美精品一区二区大全| 99热6这里只有精品| av福利片在线| 在线观看www视频免费| 一级片免费观看大全| 狂野欧美激情性bbbbbb| 日韩一本色道免费dvd| 最近的中文字幕免费完整| 亚洲精品自拍成人| 久久99热6这里只有精品| 亚洲美女黄色视频免费看| 韩国高清视频一区二区三区| 99热国产这里只有精品6| 欧美成人午夜免费资源| 黑人高潮一二区| 精品一区二区三卡| 搡老乐熟女国产| 中国美白少妇内射xxxbb| 久久人人爽av亚洲精品天堂| 欧美激情极品国产一区二区三区 | 国产女主播在线喷水免费视频网站| 曰老女人黄片| freevideosex欧美| 五月伊人婷婷丁香| 大香蕉久久成人网| 亚洲精品成人av观看孕妇| 婷婷色麻豆天堂久久| 最近的中文字幕免费完整| 大香蕉久久成人网| 中文字幕亚洲精品专区| 人人妻人人澡人人看| 欧美成人午夜免费资源| 日韩中文字幕视频在线看片| 日韩欧美一区视频在线观看| 国国产精品蜜臀av免费| 免费在线观看完整版高清| 人人妻人人添人人爽欧美一区卜| 少妇高潮的动态图| 国产精品久久久久成人av| 亚洲第一区二区三区不卡| 国产 精品1| freevideosex欧美| 日韩制服丝袜自拍偷拍| 爱豆传媒免费全集在线观看| 久久婷婷青草| 成人18禁高潮啪啪吃奶动态图| 亚洲精品美女久久久久99蜜臀 | 99视频精品全部免费 在线| 老司机亚洲免费影院| 欧美精品人与动牲交sv欧美| 亚洲五月色婷婷综合| 69精品国产乱码久久久| 黑人巨大精品欧美一区二区蜜桃 | 日韩在线高清观看一区二区三区| av一本久久久久| 伦精品一区二区三区| 一级毛片电影观看| 亚洲国产精品国产精品| 亚洲天堂av无毛| 男人爽女人下面视频在线观看| 少妇被粗大的猛进出69影院 | 亚洲婷婷狠狠爱综合网| 制服丝袜香蕉在线| 人妻系列 视频| av福利片在线| 久久久国产欧美日韩av| 97在线视频观看| 久久亚洲国产成人精品v| 久久久久久久精品精品| 插逼视频在线观看| av女优亚洲男人天堂| 国产精品国产三级专区第一集| 国产精品久久久久久久久免| 我的女老师完整版在线观看| 国产免费又黄又爽又色| 欧美另类一区| 久久鲁丝午夜福利片| 亚洲av免费高清在线观看| 亚洲一区二区三区欧美精品| 18禁观看日本| 国产男人的电影天堂91| 欧美激情国产日韩精品一区| 久久国产亚洲av麻豆专区| 亚洲色图 男人天堂 中文字幕 | 香蕉精品网在线| 精品一区二区三区四区五区乱码 | 人人妻人人澡人人看| 成人综合一区亚洲| 9色porny在线观看| 国产成人精品在线电影| a级毛片黄视频| av播播在线观看一区| 丝瓜视频免费看黄片| 自拍欧美九色日韩亚洲蝌蚪91| kizo精华| 亚洲高清免费不卡视频| 国产av码专区亚洲av| 狂野欧美激情性bbbbbb| 亚洲欧美成人综合另类久久久| 精品一区二区三卡| 久久青草综合色| 亚洲第一区二区三区不卡| 久久99热6这里只有精品| 91午夜精品亚洲一区二区三区| 午夜老司机福利剧场| 国产亚洲午夜精品一区二区久久| 99久久中文字幕三级久久日本| 亚洲精品成人av观看孕妇| 免费大片黄手机在线观看| 亚洲欧洲日产国产| 国产亚洲精品第一综合不卡 | 嫩草影院入口| 亚洲在久久综合| av一本久久久久| 国产日韩欧美在线精品| 亚洲国产av新网站| 亚洲国产欧美日韩在线播放| 街头女战士在线观看网站| 人妻系列 视频| 一级,二级,三级黄色视频| 国产精品久久久久成人av| 在线天堂最新版资源| 九色亚洲精品在线播放| 久久国产精品男人的天堂亚洲 | 成人午夜精彩视频在线观看| 日本av手机在线免费观看| 国产男女超爽视频在线观看| a 毛片基地| 中文字幕人妻熟女乱码| 久久精品熟女亚洲av麻豆精品| 丰满迷人的少妇在线观看| 欧美日韩综合久久久久久| 欧美精品亚洲一区二区| 婷婷色综合大香蕉| 少妇精品久久久久久久| 亚洲av电影在线观看一区二区三区| 久久女婷五月综合色啪小说| a 毛片基地| 免费在线观看黄色视频的| av在线老鸭窝| 观看美女的网站| 夫妻性生交免费视频一级片| 老司机影院成人| 香蕉国产在线看| 极品少妇高潮喷水抽搐| 久久久精品94久久精品| 精品酒店卫生间| 丝袜脚勾引网站| 26uuu在线亚洲综合色| 日本欧美视频一区| 黄色怎么调成土黄色| 狠狠婷婷综合久久久久久88av| 亚洲av综合色区一区| 爱豆传媒免费全集在线观看| 国产在线视频一区二区| 亚洲欧美中文字幕日韩二区| 欧美人与性动交α欧美精品济南到 | 人人澡人人妻人| 成年美女黄网站色视频大全免费| 黄色 视频免费看| 久久人人爽人人爽人人片va| 亚洲中文av在线| 亚洲精品第二区| 亚洲一级一片aⅴ在线观看| 精品少妇黑人巨大在线播放| 亚洲精品日本国产第一区| 日韩制服骚丝袜av| 国产精品国产av在线观看| 婷婷色综合www| 国产精品成人在线| 欧美97在线视频| 免费观看无遮挡的男女| 久久久亚洲精品成人影院| videos熟女内射| 欧美日韩视频高清一区二区三区二| 午夜福利在线观看免费完整高清在| 如何舔出高潮| 男女高潮啪啪啪动态图| 免费女性裸体啪啪无遮挡网站| 蜜桃国产av成人99|