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

    Stochastic Epidemic Model of Covid-19 via the Reservoir-People Transmission Network

    2022-08-24 12:56:42KazemNouriMiladFahimiLeilaTorkzadehandDumitruBaleanu
    Computers Materials&Continua 2022年7期

    Kazem Nouri, Milad Fahimi, Leila Torkzadehand Dumitru Baleanu

    1Department of Mathematics, Faculty of Mathematics, Statistics and Computer Sciences, Semnan University, P.O.Box 35195-363, Semnan, Iran

    2Department of Mathematics, Faculty of Arts and Sciences, Cankaya University, Ankara, Turkey

    3Institute of Space Sciences, Magurele-Bucharest, Romania

    Abstract: The novel Coronavirus COVID-19 emerged in Wuhan, China in December 2019.COVID-19 has rapidly spread among human populations and other mammals.The outbreak of COVID-19 has become a global challenge.Mathematical models of epidemiological systems enable studying and predicting the potential spread of disease.Modeling and predicting the evolution of COVID-19 epidemics in near real-time is a scientific challenge, this requires a deep understanding of the dynamics of pandemics and the possibility that the diffusion process can be completely random.In this paper, we develop and analyze a model to simulate the Coronavirus transmission dynamics based on Reservoir-People transmission network.When faced with a potential outbreak, decision-makers need to be able to trust mathematical models for their decision-making processes.One of the most considerable characteristics of COVID-19 is its different behaviors in various countries and regions, or even in different individuals, which can be a sign of uncertain and accidental behavior in the disease outbreak.This trait reflects the existence of the capacity of transmitting perturbations across its domains.We construct a stochastic environment because of parameters random essence and introduce a stochastic version of theReservoir-Peoplemodel.Then we prove the uniqueness and existence of the solution on the stochastic model.Moreover, the equilibria of the system are considered.Also, we establish the extinction of the disease under some suitable conditions.Finally, some numerical simulation and comparison are carried out to validate the theoretical results and the possibility of comparability of the stochastic model with the deterministic model.

    Keywords: Coronavirus; infectious diseases; stochastic modeling; brownian motion; reservoir-people model; transmission simulation; stochastic differential equation

    1 Introduction

    Mathematical biology is one of the most interesting research areas for applied mathematicians.In these research areas, modeling infectious diseases are considered as an important tool to describe the dynamics of epidemic transmission [1-11].These models play an important role in quantifying the efficient control and preventive measures of infectious diseases.Different scholars used various types of mathematical models in their analysis.

    The ongoing Coronavirus disease 2019 (COVID-19) outbreak, emerged in Wuhan, China at the end of 2019 and increased the attention worldwide [12].COVID-19 has rapidly spread among human populations and other mammals.Additionally, COVID-19 can be capable of spreading from person to person and between cities.Modeling and predicting the evolution of COVID-19 epidemics in near real-time is a scientific challenge [13], this requires a deep understanding of the dynamics of pandemics and the possibility that the diffusion process can be completely random [14].Nowadays, a large number of models are used to simulate the current COVID-19 pandemic.One of the common ways for analyzing the dynamic of COVID-19 is the use of ordinary differential equations (ODEs), but there are some limitations compared to stochastic models.Recently, many works related to the ODEs have been published (see, for example [15-19]).Furthermore, Tahir et al.[20] developed a mathematical model (forMERS) in form of a nonlinear system of differential equations.Chen et al.[21] developed a Bats-Hosts-Reservoir-People (BHRP) transmission network model for simulating the phase-based transmissibility of Coronavirus, which focus on calculatingR0.Additionally, Hashemizadeh et al.[22], presented a numerical solution for themathematicalmodel of the novel Coronavirus by the application of alternative Legendre polynomials to find the transmissibility of COVID-19 based on Reservoir-People network model.Also, Gao et al.[23] presented a numerical solution for the BHRP transmission network model for simulating the transmissibility of Coronavirus.

    The deterministic models make assumptions about the expected value of parameters in the future, but they ignore the variation and fluctuation about the expected value of parameters.One of the main benefits of the stochastic models is that they allow the validity of assumptions to be tasted statistically and produce an estimate with additional degree of realism.However, there are times when stochastic output cannot be thoroughly reviewed and compared to expectations.Stochastic models have advantages but they are computationally quite complex to perform and need careful consideration of outputs.The parameters of deterministic models governing the equations are supposed to be known and therefore the solutions are often unique.This limitation poses a practical problem because nature is intrinsically heterogeneous and the system is only measured at a discrete (and often small) number of locations.Stochastic models can be regarded as a tool to combine deterministic models, statistics and uncertainty within a coherent theoretical framework.So, this idea leads researchers to consider the stochastic models [3, 24-35].On the other hand, one of the most considerable characteristics of COVID-19 is its different behaviors in various countries and regions, or even in different individuals, which can be a sign of uncertain and accidental behavior in the disease outbreak.This trait reflects the existence of the capacity of transmitting perturbations across its domains.In past studies, [10-11, 29, 30, 35-39] have been considered a stochastic form of the COVID-19 outbreak.

    In this paper, we establish a model for investigation of the COVID-19 transmission based on a stochastic version of the Reservoir-People network with additional degree of realism.Also, it is important to note that an implication of our work is that the numerical simulations demonstrate the efficiency of our model and the possibility of comparability of the stochastic model with the deterministic model.We believe that the stochastic model established by this study is useful for researchers and scientists inmaking informed decisions and taking appropriate steps to dominance the COVID-19 disease.We conduct numerical simulations to support the theoretical findings.Numerical simulations demonstrate the efficiency of our model.We hope that the stochastic model established by this study will be useful for researchers and scientists in making informed decisions and taking appropriate steps to dominance the COVID-19 disease.The mathematical analysis of this stochastic model can be investigated in future researches and the simulation results can be extended to other countries involved in the global outbreak.This paper could lead to other studies that have included random and uncertainty of the model and can be more consistent with the reality of Coronavirus transmission.This paper is organized as follows.In Section 2, preliminaries are presented.In Section 3,the stochastic version of the Reservoir-People network model of transmission is presented.Section4 is devoted to present dynamical analysis of solution.Finally, in order to validate our analytical results,numerical simulations are presented in Sections 5 and 6.

    2 Preliminaries

    Let (Ω,{Ft}t≥0,P) be a complete probability space with a filtration {Ft}t≥0satisfying the usual conditions.We also let= {x∈Rd:xi>0, 1≤i≤d},d-dimensional stochastic differential equation can be expressed as follows:

    with initial valueX(t0) =X0, wheref(t,x) is a function in Rddefined on [t0, +∞)×Rd,g(t,x) is ad×mmatrix andB(t) is anm-dimensional standard Brownian motion defined on the probability space (Ω,Ft,P), alsof,gare locally Lipschitz functions inx.

    We define the differential operatorLassociated with Eq.(1) by [39],

    IfLacts on functionVinC2×1(Rd×[t0,∞)), then we obtain

    where

    The generalized It? formula implies that

    dV(t,x) =LV(t,x)dt+Vx(t,x)g(t,x)dB(t).

    3 COVID-19 Model

    In this work, we consider a model of 2019-novel Coronavirus (COVID-19) named as Reservoirpeople model.This model is the simplified and normalized form of Bats-Hosts-Reservoir-People model [21].AssumingSpandEprefer to the number of susceptible and exposed people,Ipdenotes symptomatic infected people,Apdenotes asymptomatic infected people,Rpdenotes removed people including recovered and dead people andWdenotes the COVID in the reservoir (the seafood market).These variables have the unit of number which are considered at timet.Also, the model parameters and their definitions are given in Tab.1.Furthermore, the pictorial diagram indicating of the model is shown in Fig.1.

    Table 1: Model parameters and their definitions

    Figure 1: The model’s flow diagram

    In this paper,we assumenpandmpas the rate of peoplemoving into and moving out from a region per day and the unit of measuring is (individual/day).The incubation period (wp) is the time between infection or contact with the agent and the onset of symptoms or signs of infection and the unit of measuring is (day-1).In this model, the unit of measuring,γpandis (day-1).The Coronavirus typically spreads via droplets from an infected person’s coughs or sneezes.Coronavirus could survive on different types, ranging from 2 h to 14 days.In this model, the unit of measuring ? is (day-1).Furthermore, the unit of measuring transmission rates βpand βwis the number of infections per unit time.Also, we defined asymptomatic COVID-19 positive persons as those who did not present any symptoms.The proportion δpis calculated as the number of asymptomatic COVID-19 infections at initial testing over COVID-19 positive persons.

    The process of normalization is implemented as follows [21]:

    Therefore the six-dimensional integer-order COVID-19 model can be expressed by the following ordinary differential equation:

    3.1 Implementation of Stochastic Description

    We can provide an additional degree of realism by defining the white noise and Brownian motion and introducing a stochastic model.Therefore, we implement this idea by replacing random parameters

    whereBi(t) and σi,i= 1,2,...,6 are the Brownian motions and the intensities of the white noises,respectively.These parameters are selected for the implementation of stochastic environment because of their random essence.So, the following stochastic differential equation for COVID-19 model is obtained

    4 Dynamical Analysis of Solution

    In this section, existence and uniqueness of the solution of the introduced model (2) and the equilibria of the system are presented.

    4.1 Existence and Uniqueness of the Solution

    Theorem 4.1.The coefficients of the stochastic differential Eq.(2) are locally Lipschitz.

    Proof.We define the following kernel γ(sp,ip)=sp.ip.Nowwe should prove the Lipschitz condition for γ(sp,ip).First we set d= (sp,ip) therefore γ(d) =sp.ipand γ(d)∈C1(Ω,R).So we must prove

    Now we choose suitable a subset of Ω and considerN(x,r) as an open ball centered atxand with radius r, such that

    N(x,r) = {y|‖y-x‖≤r}

    Then we continue our proof by applying:d1,d2∈N(x,r)?Ω?R2.Since γ(d) is unbounded therefore, due to (3), we replace K withKxwhich depends onx.

    So we define the path ψ: [0,1]→N(x,r) such that ψ(t) = (1 -t)d1+td2.So.We consider the following increment from d1to d2:

    As a result, there is no global Lipschitz constant K valid on all of Ω= R2in this case, thus γ(d) is in fact locally Lipschitz.

    Lemma 4.1.Suppose thatX?Rn, and let f, g:X?Rnsatisfy Lipschitz conditions, thenf+gsatisfies Lipschitz condition.

    Proof.The proof is presented in [40].

    Therefore, we have local Lipschitz condition for the coefficient ofdB1(t) and similarly, this condition will be proved for the coefficient ofdB2(t),dB3(t),dB4(t),dB5(t) anddB6(t).So our Model contains stochastic differential equations with Lipschitz coefficient and we can obtain a strong solution and for each initial value, the solution will be unique.

    Theorem 4.2.For any given initial value(sp(0),ep(0),ip(0),ap(0),rp(0),w(0))∈R6+, there exists the solution(sp(t),ep(t),ip(t),ap(t),rp(t), w(t))formodel (2) and the solution will remain inR6+with probability one.

    Proof.There exists a local solution on the interval [0,τe) where τeis the explosion time [39].Now,we need to prove τe= +∞.Letm0≥0 be sufficiently large such that each component of initial value(sp(0),ep(0),ip(0),ap(0),rp(0),w(0)) remain within the interval.For each integerm≥m0,we define

    where τmis stopping time.We denote τ+∞=τmand we will show that τ+∞= +∞is valid almost surely.Then the proof goes by contradiction.If this statement is false, then there exist constants T≥0 and η∈(0,1)such that P{τm≤T}≥η for any m≥m0.

    Let f be aC2-function, by definingf:→R such that

    f(sp(t),ep(t),ip(t),ap(t),rp(t),w(t)) =lnsp(t) +lnep(t) +lnip(t) +lnap(t) +lnrp(t) +lnw(t)

    So we obtain

    Therefore we obtain

    where K is a positive constant and we have

    Due to (7) and integrating both sides of (6) from to (τm∧T)and taking expectation deduce thatEf≤f(sp(0),ep(0),ip(0),ap(0),rp(0),w(0)) +KT

    Thus for every ω∈Ωm= {τm≤T} we obtain

    where 1Ωm(ω)is indicator function of Ω.Lettingm→+∞leads to contradiction

    +∞>f(sp(0),ep(0),ip(0),ap(0),rp(0),w(0)) +KT≥+∞.

    So the statement τ+∞= +∞is being satisfied and the proof is complete.

    4.2 Extinction of the Disease

    In this section, we establish sufficient conditions for the extinction of the disease.In order to form these conditions, we consider variation in the number of symptomatic infected people (Ip).Chen et al.[21] have proved that the deterministic model admits the basic reproduction numberR0as follows:

    In deterministic epidemiological models, basic reproduction number is planned to be a criterion of transmissibility of infectious agents.IfR0>1, then the outbreak is expected to continue ifR0<1 then the outbreak is expected to become extinct.In this stochastic COVID-19 model, we consider the effect of high intensity white noises on the prevalence or extinction of the COVID-19.So we present the following theoremto demonstrate that random effect may lead the disease to extinct under suitable conditions, although the disease may be immanent for the deterministic model.Before presenting the main theorem of this section, we refer to the following lemma presented in [39].

    Lemma4.2.Let Z(t):t≥0be a continuous local martingale[41]and<Z(t),Z(t)>beitsquadratic variation.Let 0<c<1 and k be a random integer.Then for almost all ω∈Ω, there exit a random integer k0(ω) such that for all k>k0:

    Theorem 4.3.Let (sp(0),ep(0),ip(0),ap(0),rp(0),w(0))∈R6+be the solution of model (2) with initial value (sp(0),ep(0),ip(0),ap(0),rp(0),w(0))∈R6+, we obtain

    Furthermore, if the assumptions+Υp+mpand δp>1 hold, then I(t) will tend to zero exponentially with probability one.

    Proof.It?’s formula yields that

    So we obtain

    According to Lemma 4.2, we obtain

    where 0<c<1 and k is a random integer.Applying Borel-Cantelli Lemma [41] leads to that for almost all ω∈Ω, there exit a random integer k0(ω) such that for k>k0:

    Therefore

    for 0≤t≤k.Substituting (9) into (8) yields

    So we have

    By dividing both sides of Eq.(10) by t such that k - 1≤t≤k, we achive

    Letting k→∞leads to n→∞and according to the strong law of large numbers [41] to the

    Brownian motion, we have lim= 0

    So we obtain

    Therefore, our claim holds and this completes our proof.

    Remark 4.1.Theorem 4.3 shows that the disease will decay exponentially to zero under suitable conditions.When the noises σ2and σ5are sufficiently large then the disease goes extinct exponentially.The white noises σ2and σ5are related to the fluctuation ofmp(the death rate of people) and δp(the proportion of asymtomatic infection of people), respectively.

    4.3 Equilibria of the System

    Equilibria of the system can be estimated by setting the all derivatives equal to zero.This gives the system

    In order to obtain the fixed points of COVID-19 model (2), we consider the solution of system(11) as follows

    In order to achieve biologically relevant solutions, we find equilibrium points such that this system has positive solutions.There are two equilibrium points of system (11).

    i) Disease free equilibrium pointE1= ((t), 0, 0, 0, 0, 0).In this disease free equilibrium point, the population of five compartments of categorized people tends to zero and the system can be recognized in healthy situation.

    ii) Endemic equilibrium pointE2=where its components are expressed above.

    5 Numerical Simulation

    In this section,we get deep insights into themodel’s dynamical behavior.For numerical simulation,the results of the implementation are given by the known Euler-Maruyama method [42].We consider the following discretization model fort= 0,Δt,2Δt,...,nΔ andk= 1,2,...,n.

    where ΔBi(t)~andZi~N(0,1),i= 1,2,...,6.The source code has been written in MATLAB R2017b.Data given in Tab.2 were calibrated for the cases in China, the first big source of COVID-19.On 21 February 2020 China (and WHO) informed that they have revised their guidance on case classification for COVID-19, removing the classification of clinically diagnosed previously used for Hubei province, and retaining only suspected and confirmed for all areas, the latter requiring laboratory confirmation.Some previously reported clinically diagnosed cases are thus expected to be discarded over the coming days as laboratory testing is conducted and some are found to be COVID-19-negative (see [43]).Thus, in order to use more accurate and reliable data, we use the reported data in China before October 2020.

    Table 2: The parameters values used in Example 5.1

    Example 5.1In order to illustrate the adaptation and compatibility of the presented model with the deterministic model, we used the reported data in Wuhan City, China [21].The parameters values used in the numerical simulation are given in Tab.2.

    The world health organization (WHO) reported an incubation period for COVID-19 between 2-10 days and the United States center for Disease Prevention and Control (CDC) estimates the incubation period for COVID-19 to be between 2-14 days.The mean incubation was 5.2 days (95%confidence interval CI: 4.1 to 7.0).According to the report of WHO on 3 March 2020, globally about 3.4% of reported COVID-19 cases have died.In this simulation, we assume thatwp=and we set the infectious period of the cases as 5.8 days [21].Therefore, we set γp= 0.1724 .Since the Coronavirus can live for hours to days on surfaces, so we set ?= 0.1 for 10 days lifetime of virus.Meanwhile, the transmission rates in our model depend on the epidemiological status and environmental conditions which change with time.The other parameters values used in this simulation are assumed based on the reported data in Wuhan City [21].Also, initial conditions have been chosen as (sp(0),ep(0),ip(0),ap(0),rp(0),w(0)) = (0.1, 0.1, 0.1, 0.1, 0.1, 0.1) and σi= 0.4,i= 1,2,...,6.The strength of the presented numerical method is particularly demonstrated in Fig.2 which shows the adaptation of our model with the deterministic model.

    In this example, if the intensities of the white noises of model (2) increase, the solution of model(2) will fluctuate around the disease free equilibrium.Fig.3 shows the simulations with σi= 0.6,i=1,2,...,6.

    Figure 2: (Continued)

    Figure 2: Comparison between the our presented stochastic model and deterministic model (σi= 0.4)

    Figure 3: (Continued)

    Figure 3: Comparison between the our presented stochastic model and deterministic model (σi= 0.6)

    Table 3: The parameters values used in Example 5.2

    Example 5.2Since the information is changing and due to the lack of complete data on many parameters used in the proposed model, we now present some numerical simulations to illustrate the efficiency of our approach as an experiment for considering sensitivity analysis.Sensitivity analysis is commonly used to determine the robustness of model predictions to parameter values, since there are usually errors in collected data and assumptive parameter values.The parameter values used in this example are given in Tab.3.In this example, we decrease the value of γ and γp, so the infectious period of symptomatic and asymptomatic infection of people will increase.Fig.4 shows that the time interval for occurring extinction of the disease increases when the infectious period of symptomatic and asymptomatic infection of people increases.

    Example5.3In this example we take the intensities of the white noises of the model(2)as σ2= 0.7,σ5= 0.8 and σ1=σ3=σ4=σ6= 0.4 which satisfy the condition of Theorem 4.3.The other parameter values used in this example are given in Tab.4.On the basis of simulation, Fig.5 shows that the number of symptomatic infected peopleipwill decay to zero under suitable conditions presented in Theorem 4.3.Therefore numerical simulation confirms the effect of high intensity white noises on the prevalence or extinction of the COVID-19.

    Figure 4: Effect of increasing the infectious period of symptomatic and asymptomatic infection of people

    Table 4: The parameters values used in Example 5.3

    Figure 5: Extinction of the stochastic model under suitable conditions

    6 Results and Discussion

    In the proposed approach, the results have been obtained using the known Euler-Maruyama method as mentioned above.Moreover, to ensure the accuracy and precision of the obtained results,the comparison of the proposed scheme has been done with the deterministic model.In stochastic models, there are random variations due to uncertainties in the parameter.According to the presented comparisons (see Fig.2 and Fig.3), the results of our stochastic epidemic model are convergent to the deterministic model.The simulation results demonstrate the efficiency of our model, and the possibility of comparability of the stochastic model with the deterministic model.Evaluation results and dynamical analysis confirmed its effectiveness, efficiency, and user satisfaction.In this paper,we found that the infection can be controlled by suitable health decisions and herd immunity as a naturally occurring phenomenon can make the number of infectious cases descending.Furthermore,numerical simulation confirms the effect of high intensity white noises on the prevalence or extinction of the COVID-19.From the simulated findings (please see Fig.4 and Fig.5), it can be noticed that the our stochastic model shows a strong agreement with the deterministic case at different values of white noise.Additionally, graphical investigations confirmed that the infectious population tends to zero under suitable extinction conditions presented in Theorem 4.3.When the noises σ2and σ5are sufficiently large then the disease goes extinct.The white noises σ2and σ5are related to the fluctuationofmp(te death rate of people)andδp(the proportion of asymtomatic infection of people),respectively.

    7 Conclusion and Remarks

    In this paper, we have developed and analyzed an epidemic model for simulating transmissibility of the COVID-19 based on Reservoir-People transmission network.We construct a stochastic version of the Reservoir-People model by using the white noise and Brownian motion.Also, we proved that there is a unique solution for the presented stochastic model.Moreover, the equilibria of the system are considered.Also, we establish the extinction of the disease under some suitable conditions.Finally,the numerical simulation is carried out to demonstrate the efficiency of our model, and the possibility of comparability of the stochastic model with the deterministic model.Furthermore, numerical results demonstrate the validity of suitable conditions for the extinction of the disease.Since mathematical models are powerful tools to understand the dynamics of real world phenomena, particularly the transmission of an infectious disease, therefore this paper could lead to other studies that have includedrandom and uncertainty of the model and can be more consistent with the reality of Coronavirus transmission.The presented resultsmay be fruitful for the existing outbreak in a better way and can be used in taking defensive techniques to decrease the infection.The area of stochasticmodeling has been extended recently, therefore thismodel is an indication for further study in this area.After getting high motivation from this paper and using ideas of modeling approaches in new researches (see [43-52]), hybrid stochastic fractional approach will be considered to investigate the dynamics of the COVID-19,as a scope of future works.

    Funding Statement:The authors received no specific funding for this study.

    Conflicts of Interest:The authors declare that they have no conflicts of interest to report regarding the present study.

    亚洲欧洲日产国产| 精品人妻一区二区三区麻豆| 免费看光身美女| 欧美精品人与动牲交sv欧美| 国产黄频视频在线观看| 精品国产露脸久久av麻豆| h视频一区二区三区| 国产精品.久久久| 国产日韩欧美视频二区| 狠狠精品人妻久久久久久综合| 久久久久久久亚洲中文字幕| 欧美精品国产亚洲| 黑人猛操日本美女一级片| 91在线精品国自产拍蜜月| av卡一久久| 国产 精品1| 如何舔出高潮| 亚洲第一av免费看| 免费在线观看完整版高清| 精品国产露脸久久av麻豆| 婷婷色麻豆天堂久久| 亚洲激情五月婷婷啪啪| freevideosex欧美| 母亲3免费完整高清在线观看 | 国产一区有黄有色的免费视频| av黄色大香蕉| 丝瓜视频免费看黄片| 亚洲av综合色区一区| 性色av一级| 欧美成人精品欧美一级黄| 免费看av在线观看网站| 18禁观看日本| 成人综合一区亚洲| 久久久久久久久久成人| 国产精品一区二区在线观看99| 在线天堂中文资源库| 久久精品人人爽人人爽视色| av黄色大香蕉| 日韩 亚洲 欧美在线| 男女高潮啪啪啪动态图| xxxhd国产人妻xxx| 亚洲av日韩在线播放| 老司机影院成人| av卡一久久| 一区二区三区四区激情视频| √禁漫天堂资源中文www| √禁漫天堂资源中文www| 美女xxoo啪啪120秒动态图| av国产久精品久网站免费入址| 波野结衣二区三区在线| 国产成人午夜福利电影在线观看| 久久狼人影院| 十八禁网站网址无遮挡| 大话2 男鬼变身卡| 最黄视频免费看| 欧美精品一区二区免费开放| av在线观看视频网站免费| 亚洲美女视频黄频| 高清欧美精品videossex| 国产一区有黄有色的免费视频| 狠狠婷婷综合久久久久久88av| 国产成人91sexporn| 丰满少妇做爰视频| 男女午夜视频在线观看 | 久久 成人 亚洲| 亚洲国产成人一精品久久久| 考比视频在线观看| 国产精品99久久99久久久不卡 | 国产精品一区二区在线不卡| 日韩人妻精品一区2区三区| 国产精品久久久av美女十八| 午夜激情av网站| 国产欧美日韩综合在线一区二区| 亚洲欧美成人精品一区二区| 日韩制服丝袜自拍偷拍| 欧美成人午夜免费资源| 亚洲人成网站在线观看播放| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 如何舔出高潮| 国产亚洲av片在线观看秒播厂| av视频免费观看在线观看| 亚洲精品美女久久av网站| 90打野战视频偷拍视频| 高清黄色对白视频在线免费看| 波多野结衣一区麻豆| 黑人高潮一二区| 欧美性感艳星| 久久 成人 亚洲| 欧美xxxx性猛交bbbb| 久久久久久久大尺度免费视频| 欧美激情国产日韩精品一区| 日韩制服骚丝袜av| 曰老女人黄片| 99国产精品免费福利视频| 欧美日韩精品成人综合77777| 高清不卡的av网站| 老熟女久久久| 日本黄色日本黄色录像| 成人免费观看视频高清| 不卡视频在线观看欧美| 亚洲丝袜综合中文字幕| 在线亚洲精品国产二区图片欧美| 9191精品国产免费久久| 亚洲美女视频黄频| 全区人妻精品视频| 日韩熟女老妇一区二区性免费视频| 人人妻人人爽人人添夜夜欢视频| 亚洲精品乱码久久久久久按摩| 亚洲精品久久午夜乱码| 婷婷色av中文字幕| 亚洲av电影在线进入| av免费在线看不卡| 少妇熟女欧美另类| www.av在线官网国产| 亚洲美女视频黄频| 国产老妇伦熟女老妇高清| 亚洲欧美色中文字幕在线| 伊人亚洲综合成人网| 亚洲精品国产色婷婷电影| 热re99久久精品国产66热6| 黄色一级大片看看| 精品午夜福利在线看| 少妇人妻精品综合一区二区| 看十八女毛片水多多多| 成人国语在线视频| 丰满饥渴人妻一区二区三| 国产一区二区激情短视频 | av卡一久久| 国产精品麻豆人妻色哟哟久久| 人成视频在线观看免费观看| 观看美女的网站| 熟妇人妻不卡中文字幕| 国产福利在线免费观看视频| 日韩大片免费观看网站| 久久国产亚洲av麻豆专区| av不卡在线播放| 日本色播在线视频| 国产一区二区三区综合在线观看 | 亚洲一区二区三区欧美精品| 免费黄频网站在线观看国产| 成人二区视频| 飞空精品影院首页| 夜夜爽夜夜爽视频| 亚洲精华国产精华液的使用体验| 在线天堂中文资源库| 一级片免费观看大全| 亚洲美女视频黄频| 日韩中字成人| www.熟女人妻精品国产 | 女性生殖器流出的白浆| 国产极品天堂在线| 宅男免费午夜| 亚洲精品美女久久av网站| 国产精品偷伦视频观看了| 黑人欧美特级aaaaaa片| 建设人人有责人人尽责人人享有的| 午夜免费观看性视频| av播播在线观看一区| 久久人人爽av亚洲精品天堂| 巨乳人妻的诱惑在线观看| 只有这里有精品99| 好男人视频免费观看在线| 国产又爽黄色视频| 国产黄频视频在线观看| 在线 av 中文字幕| 免费大片黄手机在线观看| 日产精品乱码卡一卡2卡三| 色网站视频免费| 亚洲精品中文字幕在线视频| 国产又爽黄色视频| 国产免费福利视频在线观看| 黑人高潮一二区| 久久99热这里只频精品6学生| 久久国内精品自在自线图片| 国产无遮挡羞羞视频在线观看| 久久久国产欧美日韩av| 国产一区二区激情短视频 | 80岁老熟妇乱子伦牲交| 免费观看在线日韩| 日韩精品有码人妻一区| 精品少妇久久久久久888优播| 少妇 在线观看| 久久精品国产a三级三级三级| 亚洲,欧美,日韩| 22中文网久久字幕| 亚洲av男天堂| 国语对白做爰xxxⅹ性视频网站| 韩国精品一区二区三区 | 晚上一个人看的免费电影| 99国产综合亚洲精品| 97在线视频观看| 免费黄网站久久成人精品| 人成视频在线观看免费观看| 亚洲 欧美一区二区三区| 免费少妇av软件| 日本与韩国留学比较| 精品一区在线观看国产| 亚洲精品456在线播放app| 最近中文字幕2019免费版| 国产一区二区在线观看日韩| 国产精品嫩草影院av在线观看| 深夜精品福利| 日韩av免费高清视频| 香蕉精品网在线| 久久人妻熟女aⅴ| 在线观看美女被高潮喷水网站| av黄色大香蕉| 久久99精品国语久久久| 国产精品嫩草影院av在线观看| 久久久久精品人妻al黑| 久久久精品94久久精品| av黄色大香蕉| 男人操女人黄网站| 日韩 亚洲 欧美在线| 国产精品一区二区在线观看99| 一级毛片电影观看| 国产日韩欧美在线精品| 黄片无遮挡物在线观看| 色婷婷久久久亚洲欧美| 久久国产精品大桥未久av| 久久午夜综合久久蜜桃| 精品视频人人做人人爽| 国产精品国产三级国产av玫瑰| 丰满饥渴人妻一区二区三| 免费观看在线日韩| 久久99热这里只频精品6学生| 精品国产一区二区三区四区第35| 久久久久久伊人网av| 丰满饥渴人妻一区二区三| 午夜福利在线观看免费完整高清在| 人人妻人人添人人爽欧美一区卜| 国产精品国产av在线观看| 亚洲国产精品专区欧美| 亚洲国产成人一精品久久久| 久久99蜜桃精品久久| 国产亚洲精品久久久com| 亚洲av电影在线观看一区二区三区| 九九爱精品视频在线观看| 视频在线观看一区二区三区| 丰满少妇做爰视频| 国产欧美亚洲国产| 色吧在线观看| av在线播放精品| av免费在线看不卡| 久久精品久久精品一区二区三区| 性色avwww在线观看| 久久精品国产鲁丝片午夜精品| 国产成人一区二区在线| 国产极品天堂在线| 9色porny在线观看| 久久人人爽人人片av| 精品视频人人做人人爽| 汤姆久久久久久久影院中文字幕| 天美传媒精品一区二区| 黑人巨大精品欧美一区二区蜜桃 | 欧美变态另类bdsm刘玥| 少妇高潮的动态图| 精品国产一区二区久久| 久久午夜综合久久蜜桃| 在线免费观看不下载黄p国产| 在线观看免费视频网站a站| 成人亚洲欧美一区二区av| 国产精品秋霞免费鲁丝片| 美女内射精品一级片tv| 国产成人精品一,二区| 天堂中文最新版在线下载| 色哟哟·www| 三级国产精品片| 伦精品一区二区三区| 国产精品无大码| 亚洲色图综合在线观看| 亚洲av在线观看美女高潮| 香蕉丝袜av| 啦啦啦在线观看免费高清www| 侵犯人妻中文字幕一二三四区| 久久亚洲国产成人精品v| 看非洲黑人一级黄片| 99九九在线精品视频| 天天影视国产精品| 五月伊人婷婷丁香| 国产成人欧美| 国产精品偷伦视频观看了| 丝袜脚勾引网站| 色网站视频免费| av在线播放精品| 宅男免费午夜| 制服人妻中文乱码| 少妇的逼水好多| 在线观看三级黄色| 亚洲第一av免费看| 大码成人一级视频| 日本91视频免费播放| 日韩中字成人| 亚洲国产欧美日韩在线播放| av在线观看视频网站免费| 久久久久精品性色| 高清黄色对白视频在线免费看| 成人无遮挡网站| 国产毛片在线视频| 国产一区二区在线观看av| 国产熟女欧美一区二区| 五月伊人婷婷丁香| 久久 成人 亚洲| 久久精品熟女亚洲av麻豆精品| 九九在线视频观看精品| 久久久久视频综合| 一个人免费看片子| 中国三级夫妇交换| 伊人久久国产一区二区| 国产成人91sexporn| 国产爽快片一区二区三区| 只有这里有精品99| 国产亚洲一区二区精品| 欧美日韩成人在线一区二区| 亚洲欧美日韩另类电影网站| 汤姆久久久久久久影院中文字幕| 国产日韩欧美视频二区| 只有这里有精品99| 热re99久久精品国产66热6| 亚洲激情五月婷婷啪啪| 一级黄片播放器| 母亲3免费完整高清在线观看 | 久久精品国产综合久久久 | 国产精品嫩草影院av在线观看| 国产av精品麻豆| 美女福利国产在线| 亚洲av电影在线进入| 国产成人免费无遮挡视频| 18禁国产床啪视频网站| 黑人欧美特级aaaaaa片| 亚洲在久久综合| 国产高清国产精品国产三级| a级片在线免费高清观看视频| 亚洲av福利一区| videos熟女内射| 国产精品无大码| 性色av一级| av在线观看视频网站免费| 亚洲精品成人av观看孕妇| 久久久亚洲精品成人影院| 亚洲四区av| 肉色欧美久久久久久久蜜桃| 国产亚洲精品久久久com| 一级毛片黄色毛片免费观看视频| 丰满少妇做爰视频| 免费观看在线日韩| 国产片特级美女逼逼视频| 亚洲av日韩在线播放| av不卡在线播放| 免费观看在线日韩| 美女国产高潮福利片在线看| 日韩一本色道免费dvd| 9热在线视频观看99| 久久人妻熟女aⅴ| 欧美人与性动交α欧美软件 | 青春草国产在线视频| a级毛片在线看网站| 欧美精品亚洲一区二区| 亚洲精品久久久久久婷婷小说| 少妇 在线观看| 亚洲av中文av极速乱| 99热6这里只有精品| 欧美精品一区二区大全| 丝袜脚勾引网站| 22中文网久久字幕| 日韩一区二区三区影片| 欧美精品一区二区免费开放| 中文字幕亚洲精品专区| 五月伊人婷婷丁香| 韩国精品一区二区三区 | 精品人妻在线不人妻| 国国产精品蜜臀av免费| 欧美激情 高清一区二区三区| 青春草亚洲视频在线观看| 欧美变态另类bdsm刘玥| 久久人人爽av亚洲精品天堂| www.av在线官网国产| 午夜福利视频精品| 国产有黄有色有爽视频| 亚洲国产色片| 2022亚洲国产成人精品| 丝袜在线中文字幕| 免费观看a级毛片全部| 国产成人欧美| 国产精品秋霞免费鲁丝片| 亚洲国产精品国产精品| 精品久久蜜臀av无| 国产av一区二区精品久久| 精品一区在线观看国产| 91精品三级在线观看| 亚洲,一卡二卡三卡| 性高湖久久久久久久久免费观看| 嫩草影院入口| 欧美xxⅹ黑人| 国产成人aa在线观看| 这个男人来自地球电影免费观看 | 黄色视频在线播放观看不卡| 蜜臀久久99精品久久宅男| 久久精品久久久久久久性| 亚洲熟女精品中文字幕| 乱人伦中国视频| av免费在线看不卡| 成人18禁高潮啪啪吃奶动态图| 18+在线观看网站| 狂野欧美激情性xxxx在线观看| 久久av网站| 视频中文字幕在线观看| 90打野战视频偷拍视频| 国产精品熟女久久久久浪| 男的添女的下面高潮视频| 久久久久国产精品人妻一区二区| 丝袜人妻中文字幕| 蜜桃国产av成人99| 菩萨蛮人人尽说江南好唐韦庄| 亚洲精品国产av蜜桃| www.色视频.com| 国产精品久久久av美女十八| 丰满乱子伦码专区| 深夜精品福利| 18禁在线无遮挡免费观看视频| 国产精品久久久久久久电影| 熟女人妻精品中文字幕| 久久99蜜桃精品久久| 欧美 日韩 精品 国产| 日本wwww免费看| 国产精品麻豆人妻色哟哟久久| 少妇的逼好多水| 热99久久久久精品小说推荐| 一区二区三区四区激情视频| 久久久久久久久久久免费av| 高清欧美精品videossex| 边亲边吃奶的免费视频| 性色av一级| 最新中文字幕久久久久| 在线观看一区二区三区激情| 欧美国产精品一级二级三级| 欧美+日韩+精品| 亚洲色图综合在线观看| 在线天堂最新版资源| 久久久久精品性色| 如日韩欧美国产精品一区二区三区| 成年av动漫网址| 国产精品国产三级国产av玫瑰| 久久精品熟女亚洲av麻豆精品| av天堂久久9| 看免费成人av毛片| 99久久精品国产国产毛片| 一级片免费观看大全| 性色av一级| 精品国产一区二区久久| 成人综合一区亚洲| 久久热在线av| 国产欧美亚洲国产| 国内精品宾馆在线| 在线观看免费日韩欧美大片| 老熟女久久久| 一级a做视频免费观看| 亚洲欧美成人精品一区二区| 在线观看免费高清a一片| 飞空精品影院首页| 精品一区二区免费观看| 亚洲精品久久成人aⅴ小说| 搡老乐熟女国产| 亚洲国产日韩一区二区| 春色校园在线视频观看| 一级,二级,三级黄色视频| 少妇熟女欧美另类| 少妇人妻 视频| 国产xxxxx性猛交| 久久午夜综合久久蜜桃| 国产精品一区二区在线观看99| 老司机影院成人| 免费高清在线观看视频在线观看| 免费黄网站久久成人精品| 精品久久蜜臀av无| 人妻一区二区av| 好男人视频免费观看在线| 亚洲精品aⅴ在线观看| 热re99久久精品国产66热6| 精品一区二区三区四区五区乱码 | 成人国产av品久久久| 飞空精品影院首页| 亚洲成国产人片在线观看| 热99久久久久精品小说推荐| 美女福利国产在线| 男的添女的下面高潮视频| 婷婷成人精品国产| 免费不卡的大黄色大毛片视频在线观看| 赤兔流量卡办理| 一二三四中文在线观看免费高清| 肉色欧美久久久久久久蜜桃| 久久99热6这里只有精品| 日韩不卡一区二区三区视频在线| 男人添女人高潮全过程视频| 看免费av毛片| 免费女性裸体啪啪无遮挡网站| 欧美成人午夜精品| 色婷婷av一区二区三区视频| 午夜免费男女啪啪视频观看| 99久久综合免费| 午夜福利网站1000一区二区三区| 婷婷成人精品国产| av又黄又爽大尺度在线免费看| 国产成人午夜福利电影在线观看| 美女内射精品一级片tv| 1024视频免费在线观看| 日韩制服丝袜自拍偷拍| 欧美国产精品一级二级三级| 女人久久www免费人成看片| 中文乱码字字幕精品一区二区三区| 在线观看免费日韩欧美大片| 日本91视频免费播放| 久久影院123| 国产白丝娇喘喷水9色精品| 国产亚洲精品久久久com| 日本与韩国留学比较| 91成人精品电影| 国产日韩一区二区三区精品不卡| 99国产综合亚洲精品| 一级黄片播放器| 精品一区二区免费观看| 亚洲精华国产精华液的使用体验| 成人无遮挡网站| 最新的欧美精品一区二区| 久久人人97超碰香蕉20202| 国产精品三级大全| av有码第一页| 国产精品99久久99久久久不卡 | 两性夫妻黄色片 | 精品一区二区三区四区五区乱码 | 国产亚洲欧美精品永久| 久久精品国产亚洲av涩爱| av不卡在线播放| 亚洲色图 男人天堂 中文字幕 | 国产成人91sexporn| 免费看不卡的av| 咕卡用的链子| av播播在线观看一区| 国产精品麻豆人妻色哟哟久久| 人人澡人人妻人| 天堂中文最新版在线下载| 妹子高潮喷水视频| 777米奇影视久久| 看免费成人av毛片| 午夜激情av网站| 日本wwww免费看| 久久久欧美国产精品| 亚洲高清免费不卡视频| 国产精品国产三级专区第一集| videos熟女内射| 国产综合精华液| 大话2 男鬼变身卡| 精品99又大又爽又粗少妇毛片| √禁漫天堂资源中文www| 人人妻人人澡人人爽人人夜夜| 亚洲五月色婷婷综合| 最新中文字幕久久久久| 热99久久久久精品小说推荐| av电影中文网址| 国产一区二区三区av在线| 夜夜爽夜夜爽视频| 国产免费视频播放在线视频| 午夜日本视频在线| 国产日韩欧美亚洲二区| 在线观看国产h片| 免费观看在线日韩| 久久久久久人人人人人| 丰满乱子伦码专区| 国产熟女欧美一区二区| 黄色 视频免费看| 自拍欧美九色日韩亚洲蝌蚪91| 国产在线免费精品| 99热6这里只有精品| 国产精品国产三级国产av玫瑰| 十八禁网站网址无遮挡| 97人妻天天添夜夜摸| 国产成人91sexporn| 欧美人与善性xxx| 亚洲欧美中文字幕日韩二区| 亚洲伊人久久精品综合| 国产成人一区二区在线| av女优亚洲男人天堂| 国产麻豆69| 国产成人一区二区在线| 国产精品久久久久久精品电影小说| 日本wwww免费看| 久久av网站| 久久久久精品性色| 亚洲国产精品专区欧美| 精品少妇黑人巨大在线播放| 熟女人妻精品中文字幕| 99久久中文字幕三级久久日本| 捣出白浆h1v1| 水蜜桃什么品种好| 色吧在线观看| 岛国毛片在线播放| 麻豆乱淫一区二区| 啦啦啦在线观看免费高清www| 久久久久久久久久久免费av| 99热这里只有是精品在线观看| 国产成人精品婷婷| 亚洲人成77777在线视频| 国产精品一二三区在线看| xxx大片免费视频| 久久久国产精品麻豆| 在线观看www视频免费| 有码 亚洲区| 全区人妻精品视频| 三上悠亚av全集在线观看| 中国国产av一级| 天天操日日干夜夜撸| 极品少妇高潮喷水抽搐| 亚洲 欧美一区二区三区| 18在线观看网站| 国产精品久久久av美女十八|