• <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.

    人体艺术视频欧美日本| 亚洲精品乱码久久久久久按摩| 一卡2卡三卡四卡精品乱码亚洲| 国产欧美另类精品又又久久亚洲欧美| 99热这里只有精品一区| 欧美人与善性xxx| 日韩亚洲欧美综合| 久久久久久久久久黄片| 伊人久久精品亚洲午夜| 成人二区视频| 有码 亚洲区| 日本黄色视频三级网站网址| 欧美成人免费av一区二区三区| 亚洲伊人久久精品综合 | 亚洲国产精品成人久久小说| 神马国产精品三级电影在线观看| 联通29元200g的流量卡| 久久精品国产鲁丝片午夜精品| 国产黄色视频一区二区在线观看 | 老司机影院成人| 日韩制服骚丝袜av| 免费av不卡在线播放| 精品久久久久久久久久久久久| 日韩高清综合在线| 国产精品无大码| 午夜视频国产福利| 在线免费观看不下载黄p国产| 国国产精品蜜臀av免费| 午夜福利在线在线| 欧美激情久久久久久爽电影| 欧美成人精品欧美一级黄| 亚洲最大成人av| 国产三级在线视频| 欧美+日韩+精品| 国产片特级美女逼逼视频| 精品久久久久久电影网 | 99在线人妻在线中文字幕| 91久久精品电影网| 久久久a久久爽久久v久久| 午夜精品在线福利| 国产成人精品婷婷| 精品国产一区二区三区久久久樱花 | 日本免费在线观看一区| 午夜福利成人在线免费观看| 男女视频在线观看网站免费| 色综合亚洲欧美另类图片| 久久精品91蜜桃| 亚洲无线观看免费| 国产av一区在线观看免费| 好男人在线观看高清免费视频| 美女xxoo啪啪120秒动态图| 国产伦精品一区二区三区视频9| 国产免费一级a男人的天堂| 91狼人影院| 国产精品麻豆人妻色哟哟久久 | 国产黄色小视频在线观看| 亚洲av.av天堂| 日产精品乱码卡一卡2卡三| av福利片在线观看| 久久久精品大字幕| videossex国产| 哪个播放器可以免费观看大片| av专区在线播放| 亚洲aⅴ乱码一区二区在线播放| 国产av在哪里看| 久久精品熟女亚洲av麻豆精品 | 成年女人看的毛片在线观看| 中文字幕亚洲精品专区| 日本爱情动作片www.在线观看| 亚洲精品日韩在线中文字幕| 日韩成人伦理影院| 蜜桃久久精品国产亚洲av| videossex国产| 国产av不卡久久| 亚洲精品,欧美精品| 看片在线看免费视频| 最近手机中文字幕大全| 国产午夜精品久久久久久一区二区三区| 亚洲丝袜综合中文字幕| videossex国产| 男女下面进入的视频免费午夜| 亚洲av日韩在线播放| 在线播放无遮挡| av黄色大香蕉| 久久久久性生活片| av卡一久久| 变态另类丝袜制服| 亚洲精品久久久久久婷婷小说 | 欧美激情久久久久久爽电影| 午夜精品一区二区三区免费看| 免费在线观看成人毛片| 日韩 亚洲 欧美在线| 在线免费观看不下载黄p国产| www.色视频.com| 欧美精品一区二区大全| 欧美人与善性xxx| 色尼玛亚洲综合影院| 最近的中文字幕免费完整| 听说在线观看完整版免费高清| 国产一区二区三区av在线| 精华霜和精华液先用哪个| 亚洲成人av在线免费| 91精品一卡2卡3卡4卡| 久热久热在线精品观看| 91久久精品国产一区二区成人| 日韩精品有码人妻一区| 国产极品精品免费视频能看的| 97超碰精品成人国产| 亚洲精品乱码久久久v下载方式| 午夜爱爱视频在线播放| 天美传媒精品一区二区| 国国产精品蜜臀av免费| 男女下面进入的视频免费午夜| 日本黄大片高清| 久久久久国产网址| 亚洲av中文字字幕乱码综合| 老女人水多毛片| 国产精品蜜桃在线观看| 少妇裸体淫交视频免费看高清| 久久这里只有精品中国| 日韩一本色道免费dvd| 在线a可以看的网站| 日韩在线高清观看一区二区三区| 欧美成人精品欧美一级黄| 欧美色视频一区免费| 嫩草影院入口| 亚洲av电影在线观看一区二区三区 | 纵有疾风起免费观看全集完整版 | 一边摸一边抽搐一进一小说| 一边亲一边摸免费视频| 久久精品国产自在天天线| 成年女人看的毛片在线观看| 亚洲,欧美,日韩| 又黄又爽又刺激的免费视频.| 99热这里只有是精品在线观看| 国内精品美女久久久久久| 国产午夜精品一二区理论片| 大又大粗又爽又黄少妇毛片口| 国产欧美日韩精品一区二区| 欧美色视频一区免费| 国产乱人视频| 欧美极品一区二区三区四区| 欧美性猛交黑人性爽| 精品人妻一区二区三区麻豆| 搞女人的毛片| 午夜a级毛片| av黄色大香蕉| 久久综合国产亚洲精品| 国产女主播在线喷水免费视频网站 | 亚洲av成人精品一区久久| 亚洲最大成人手机在线| 亚洲va在线va天堂va国产| 女人十人毛片免费观看3o分钟| 最近中文字幕高清免费大全6| 最近视频中文字幕2019在线8| 高清av免费在线| 青春草亚洲视频在线观看| 国产精品三级大全| 亚洲内射少妇av| 国产精品三级大全| 精品久久国产蜜桃| 国国产精品蜜臀av免费| 国产一区二区在线av高清观看| av又黄又爽大尺度在线免费看 | 色吧在线观看| 综合色av麻豆| 国产视频首页在线观看| 久久久成人免费电影| 国产 一区精品| 国产 一区精品| 成人三级黄色视频| 日韩欧美精品v在线| 亚洲欧美精品专区久久| 日韩一本色道免费dvd| 亚洲欧美精品自产自拍| av女优亚洲男人天堂| 午夜免费激情av| 久久久久网色| 一本久久精品| 日本欧美国产在线视频| 精品欧美国产一区二区三| 亚洲欧美日韩高清专用| 99久久精品国产国产毛片| 日韩亚洲欧美综合| 麻豆乱淫一区二区| 国产视频内射| 国产精品永久免费网站| 蜜桃久久精品国产亚洲av| 搞女人的毛片| 日韩欧美在线乱码| 免费一级毛片在线播放高清视频| 午夜精品一区二区三区免费看| 精品国内亚洲2022精品成人| 成人午夜精彩视频在线观看| 可以在线观看毛片的网站| 小蜜桃在线观看免费完整版高清| 午夜免费男女啪啪视频观看| 国产人妻一区二区三区在| 99久久精品国产国产毛片| 嫩草影院入口| 免费不卡的大黄色大毛片视频在线观看 | 国产精品爽爽va在线观看网站| 亚洲av成人av| 国产国拍精品亚洲av在线观看| 又粗又硬又长又爽又黄的视频| 色吧在线观看| 亚洲国产精品专区欧美| 我的女老师完整版在线观看| 亚洲av免费在线观看| 成人三级黄色视频| 免费观看a级毛片全部| 免费观看a级毛片全部| 成人三级黄色视频| 婷婷色综合大香蕉| 七月丁香在线播放| 亚洲av男天堂| 丰满人妻一区二区三区视频av| 天堂√8在线中文| 久久久精品94久久精品| 一级二级三级毛片免费看| 一夜夜www| 性色avwww在线观看| 亚洲成人久久爱视频| 秋霞在线观看毛片| 免费看日本二区| 99久久人妻综合| 真实男女啪啪啪动态图| 日韩一区二区视频免费看| 91精品一卡2卡3卡4卡| 白带黄色成豆腐渣| 欧美三级亚洲精品| 蜜臀久久99精品久久宅男| 久久99热这里只有精品18| 精品久久国产蜜桃| 99九九线精品视频在线观看视频| 一二三四中文在线观看免费高清| 春色校园在线视频观看| a级毛色黄片| 嘟嘟电影网在线观看| 欧美高清性xxxxhd video| 精品久久久噜噜| 国产精品人妻久久久久久| 亚洲av二区三区四区| 精品一区二区三区人妻视频| 在线观看av片永久免费下载| 免费观看a级毛片全部| 菩萨蛮人人尽说江南好唐韦庄 | 变态另类丝袜制服| 青青草视频在线视频观看| 日韩 亚洲 欧美在线| 国产精品嫩草影院av在线观看| av在线亚洲专区| 国产黄片美女视频| 久久久午夜欧美精品| 男女下面进入的视频免费午夜| 亚洲av成人精品一二三区| 国产视频首页在线观看| 一夜夜www| 欧美成人一区二区免费高清观看| 久久精品熟女亚洲av麻豆精品 | 久久精品国产自在天天线| 在线免费十八禁| 黄色一级大片看看| 国产一区二区亚洲精品在线观看| 级片在线观看| 亚洲精品成人久久久久久| 国模一区二区三区四区视频| 久久久久国产网址| 97在线视频观看| 亚洲国产成人一精品久久久| 日韩一本色道免费dvd| 男女视频在线观看网站免费| 欧美激情国产日韩精品一区| 久久精品综合一区二区三区| 天堂网av新在线| 国产 一区精品| 好男人在线观看高清免费视频| 九草在线视频观看| av在线天堂中文字幕| 国产色爽女视频免费观看| 精品人妻视频免费看| 亚洲欧美日韩高清专用| 免费观看在线日韩| 久久午夜福利片| 中国美白少妇内射xxxbb| 国产一区亚洲一区在线观看| 国产精品永久免费网站| www日本黄色视频网| 久久精品夜色国产| 精品一区二区三区人妻视频| 亚洲婷婷狠狠爱综合网| 亚洲国产精品合色在线| 国产精品女同一区二区软件| 欧美变态另类bdsm刘玥| 久久精品影院6| 亚洲不卡免费看| 国产三级在线视频| 最近2019中文字幕mv第一页| 三级国产精品欧美在线观看| 午夜老司机福利剧场| 搞女人的毛片| 狂野欧美激情性xxxx在线观看| 久久久a久久爽久久v久久| 日本色播在线视频| 久久99热6这里只有精品| 高清午夜精品一区二区三区| 欧美成人a在线观看| 欧美激情国产日韩精品一区| av在线天堂中文字幕| 国产精品电影一区二区三区| 亚洲精品色激情综合| 免费看日本二区| 日日啪夜夜撸| 国产精品美女特级片免费视频播放器| 国产一级毛片七仙女欲春2| 丝袜美腿在线中文| 欧美人与善性xxx| 亚洲欧洲日产国产| 日韩视频在线欧美| 亚洲自偷自拍三级| 国产综合懂色| 亚洲欧美成人综合另类久久久 | 免费播放大片免费观看视频在线观看 | 卡戴珊不雅视频在线播放| 亚洲欧美成人精品一区二区| 91久久精品国产一区二区三区| 日本一本二区三区精品| 亚洲综合色惰| 国产精品永久免费网站| 国产伦精品一区二区三区视频9| 亚洲一区高清亚洲精品| 人人妻人人澡欧美一区二区| 欧美xxxx性猛交bbbb| 夜夜看夜夜爽夜夜摸| 久久6这里有精品| 亚洲精华国产精华液的使用体验| 免费看a级黄色片| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 少妇的逼水好多| 熟妇人妻久久中文字幕3abv| 精品酒店卫生间| 亚洲真实伦在线观看| 青春草亚洲视频在线观看| 午夜激情欧美在线| 欧美日本视频| 久久久午夜欧美精品| a级毛片免费高清观看在线播放| 午夜爱爱视频在线播放| 国产激情偷乱视频一区二区| 久久热精品热| 久久久久久国产a免费观看| 人妻少妇偷人精品九色| 国产激情偷乱视频一区二区| 成年女人永久免费观看视频| a级一级毛片免费在线观看| 亚洲av免费在线观看| 国产在线一区二区三区精 | 在线免费十八禁| 永久网站在线| 老司机影院成人| 在线天堂最新版资源| 能在线免费看毛片的网站| 看十八女毛片水多多多| 变态另类丝袜制服| 精品人妻一区二区三区麻豆| 国产精品久久久久久精品电影| 久久久精品大字幕| 亚洲最大成人中文| 久久精品久久久久久噜噜老黄 | 亚洲人成网站在线观看播放| eeuss影院久久| 国产高清国产精品国产三级 | 久久久久久久久大av| 久久久精品94久久精品| 国产精品久久久久久精品电影小说 | 国产亚洲av嫩草精品影院| 中文在线观看免费www的网站| 亚洲无线观看免费| 国产午夜精品久久久久久一区二区三区| 91久久精品国产一区二区成人| 一级黄片播放器| 99久久精品一区二区三区| 亚洲av中文av极速乱| 丝袜美腿在线中文| 乱码一卡2卡4卡精品| 午夜福利高清视频| 乱人视频在线观看| 蜜桃亚洲精品一区二区三区| 尾随美女入室| 亚洲国产欧美在线一区| 直男gayav资源| 蜜臀久久99精品久久宅男| 免费搜索国产男女视频| 天堂影院成人在线观看| 国产午夜福利久久久久久| 国产一区亚洲一区在线观看| 少妇人妻精品综合一区二区| 久久鲁丝午夜福利片| 天堂√8在线中文| 99久国产av精品| 色网站视频免费| 日韩强制内射视频| 久久精品久久久久久久性| 三级男女做爰猛烈吃奶摸视频| 色综合站精品国产| 中文字幕免费在线视频6| 日本免费一区二区三区高清不卡| 毛片一级片免费看久久久久| 午夜精品一区二区三区免费看| 国产一区二区三区av在线| 国产人妻一区二区三区在| 白带黄色成豆腐渣| 亚洲av电影在线观看一区二区三区 | 精品不卡国产一区二区三区| 搡老妇女老女人老熟妇| 日韩视频在线欧美| 五月伊人婷婷丁香| 在线免费观看不下载黄p国产| 成人亚洲欧美一区二区av| 嫩草影院精品99| 亚洲在线观看片| 精品国产露脸久久av麻豆 | 国内揄拍国产精品人妻在线| 99热6这里只有精品| 在线a可以看的网站| 99九九线精品视频在线观看视频| 国产亚洲av片在线观看秒播厂 | 日本黄色视频三级网站网址| 国产黄色小视频在线观看| 精品国内亚洲2022精品成人| 黄片无遮挡物在线观看| 一区二区三区免费毛片| 午夜视频国产福利| 亚洲精品乱久久久久久| 午夜激情福利司机影院| 国产伦一二天堂av在线观看| 中文欧美无线码| 视频中文字幕在线观看| 国产v大片淫在线免费观看| 熟妇人妻久久中文字幕3abv| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 成人三级黄色视频| 最后的刺客免费高清国语| 免费人成在线观看视频色| 国产 一区精品| av在线播放精品| 又爽又黄无遮挡网站| h日本视频在线播放| АⅤ资源中文在线天堂| 日本黄色片子视频| 白带黄色成豆腐渣| 最近视频中文字幕2019在线8| 一卡2卡三卡四卡精品乱码亚洲| 国产av码专区亚洲av| or卡值多少钱| 色尼玛亚洲综合影院| 大话2 男鬼变身卡| 观看美女的网站| 亚洲人成网站高清观看| 深夜a级毛片| 免费观看的影片在线观看| 国产亚洲精品av在线| 国产成人精品婷婷| 国产精品久久久久久久久免| 人妻夜夜爽99麻豆av| 岛国在线免费视频观看| 真实男女啪啪啪动态图| 精品国产三级普通话版| 亚洲欧美中文字幕日韩二区| 国产v大片淫在线免费观看| 亚洲欧美精品综合久久99| 欧美不卡视频在线免费观看| 亚洲一级一片aⅴ在线观看| 久久久色成人| av在线老鸭窝| 国产一级毛片七仙女欲春2| 长腿黑丝高跟| 人妻系列 视频| 又粗又硬又长又爽又黄的视频| 亚洲婷婷狠狠爱综合网| 久久精品国产鲁丝片午夜精品| 亚洲国产高清在线一区二区三| 亚洲av福利一区| 狂野欧美白嫩少妇大欣赏| 国产精品国产高清国产av| 美女国产视频在线观看| 麻豆成人午夜福利视频| 亚洲成av人片在线播放无| 日韩欧美三级三区| 一级av片app| 久久久久久九九精品二区国产| 亚洲欧美日韩卡通动漫| 啦啦啦韩国在线观看视频| 美女大奶头视频| 日本欧美国产在线视频| 国内精品一区二区在线观看| 少妇的逼好多水| 日本三级黄在线观看| 精华霜和精华液先用哪个| 能在线免费观看的黄片| 黄色一级大片看看| 高清av免费在线| 啦啦啦啦在线视频资源| 色网站视频免费| 国产精品,欧美在线| av国产久精品久网站免费入址| 亚洲精品456在线播放app| 国内精品一区二区在线观看| 国产不卡一卡二| 日韩亚洲欧美综合| 97超碰精品成人国产| 免费电影在线观看免费观看| 亚洲欧美精品综合久久99| 久久亚洲精品不卡| 蜜桃久久精品国产亚洲av| 亚洲精品亚洲一区二区| 一二三四中文在线观看免费高清| 观看免费一级毛片| av在线观看视频网站免费| 久久久精品94久久精品| 一区二区三区乱码不卡18| 美女被艹到高潮喷水动态| 在线观看av片永久免费下载| 精品免费久久久久久久清纯| 九九爱精品视频在线观看| 欧美精品一区二区大全| 性插视频无遮挡在线免费观看| 全区人妻精品视频| 岛国在线免费视频观看| 一级毛片aaaaaa免费看小| www日本黄色视频网| 一本一本综合久久| 高清在线视频一区二区三区 | 乱码一卡2卡4卡精品| 国产爱豆传媒在线观看| 毛片一级片免费看久久久久| 亚洲人成网站高清观看| 哪个播放器可以免费观看大片| 18禁在线无遮挡免费观看视频| 免费人成在线观看视频色| 色噜噜av男人的天堂激情| 亚洲成人久久爱视频| 欧美+日韩+精品| 久久久精品94久久精品| 真实男女啪啪啪动态图| 搡女人真爽免费视频火全软件| 免费看日本二区| 亚洲真实伦在线观看| 久久久国产成人免费| 国产熟女欧美一区二区| 国产成人a∨麻豆精品| 成人美女网站在线观看视频| 久久久久久九九精品二区国产| av免费在线看不卡| 中文乱码字字幕精品一区二区三区 | 久久精品影院6| 国产69精品久久久久777片| 国产免费一级a男人的天堂| 国产黄片视频在线免费观看| 精品久久久久久久人妻蜜臀av| 男女国产视频网站| 波多野结衣巨乳人妻| 久久久精品94久久精品| 美女被艹到高潮喷水动态| 精品久久久久久久久av| 国产一区有黄有色的免费视频 | 欧美日韩综合久久久久久| 校园人妻丝袜中文字幕| 亚洲精品一区蜜桃| 成人亚洲欧美一区二区av| 国产亚洲最大av| 国产在线一区二区三区精 | 汤姆久久久久久久影院中文字幕 | 人妻制服诱惑在线中文字幕| 亚洲av一区综合| 熟妇人妻久久中文字幕3abv| 观看美女的网站| 国产精品99久久久久久久久| 一夜夜www| 伦理电影大哥的女人| 99国产精品一区二区蜜桃av| 国产精品.久久久| 国产黄片视频在线免费观看| 少妇熟女aⅴ在线视频| 亚洲av一区综合| 欧美日韩国产亚洲二区| 欧美性猛交黑人性爽| 欧美高清成人免费视频www| 九九热线精品视视频播放| 欧美97在线视频| 国产一区二区亚洲精品在线观看| 亚洲av中文字字幕乱码综合| 最近视频中文字幕2019在线8| 国产成人一区二区在线| 免费观看性生交大片5| 国产精品.久久久| 亚洲综合精品二区| 寂寞人妻少妇视频99o| 免费黄色在线免费观看| 人人妻人人澡欧美一区二区| 老司机福利观看| 国产精品一及| 免费人成在线观看视频色| 嫩草影院精品99| 天天一区二区日本电影三级| 校园人妻丝袜中文字幕| 日韩,欧美,国产一区二区三区 | av线在线观看网站| 免费电影在线观看免费观看| 成人美女网站在线观看视频| 国产亚洲91精品色在线| 国产一区有黄有色的免费视频 | av.在线天堂| 美女cb高潮喷水在线观看| 国产乱人视频|