Our goal in this paper is to adapt to the context of the Greater Abidjan Region (Côte d'Ivoire) an existing model reflecting the evolution of the COVID-19 epidemic in Wuhan (China). This model is a deterministic compartmental model which is translated by a system of ordinary differential equations. We are reinvesting this model to obtain the parameters of the epidemic in the great Abidjan. We study some mathematical aspects of this deterministic model.
Using data from the evolution of the COVID-19 epidemic in Abidjan, we obtain constants corresponding to the Greater Abidjan Region epidemic. We are performing numerical simulations using these constants to predict the behavior of the epidemic in the Greater Abidjan Region.
The initial model does not take into account certain disturbances and sudden shocks which could disturb its behavior. In the Greater Abidjan Region, some continuous or sudden events disrupt the behavior of the epidemic. To stick to the realities of the Greater Abidjan Region, we introduce a white noise and jumps that correspond to the different disturbances that can occur. We obtain a stochastic model of COIVD-19 with jumps. We prove that our new model with jumps has a positive global solution. In certain conditions, this solution oscillates in the set of the disease-free equilibrium of the deterministic model. We perform numerical simulations which corroborate our theoretical results.
COVID-19 epidemic model, Disease-free equilibrium, Jump perturbation
This study aims at analyzing the evolution of the COVID-19 epidemic in the great Abidjan over the period of public measures application. For this reasons, we have worked on data related only to that period. For the next step, we plan to make predictions to account the evolution of this epidemic.
In [1], Pierre Magal, et al. develop a mathematic model for COVID-19 epidemic in Wuhan (China) in which the population was composed of four compartments:
• S(t) is the number of susceptible individuals at time t i.e., people who are not infected yet but might become infectious individuals in the future,
• I(t) is the number of asymptomatic infectious individuals at time t i.e., people who have contracted the disease but have not yet developed it at time t,
• R(t) is the number of reported symptomatic infectious individuals i.e., symptomatic infectious with severe symptoms at time t,
• U(t) is the number of unreported symptomatic infectious individuals i.e., symptomatic infectious with mild symptoms at time.
The instant t0=0 corresponds to December 31, 2020. The epidemic starts at instant t0.
At this time, they noted S(t0)=S0>0, I(t0)=I0>0, R(t0)=0 and U(t0)=U0≥0. Times are expressed in days and for t≥t0 they obtained:
⎧⎩⎨⎪⎪⎪⎪⎪⎪⎪⎪dS(t)=[−τS(t)[I(t)+U(t)]]dtdI(t)=[τS(t)[I(t)+U(t)]−νI(t)]dtdR(t)=[ν1I(t)−ηR(t)]dtdU(t)=[ν2I(t)−ηU(t))]dt (1.1)
τ, ν, ν1, ν2 and η are real constants on [0;1]. τ is the transmission rate.
ν is the rate at which the asymptomatic infectious individuals become symptomatic. They assume 1ν is the average time during which asymptomatic infectious are asymptomatic.
f is the fraction of symptomatic infectious that become reported symptomatic infectious. So ν1=fν is the rate at which asymptomatic infectious become reported symptomatic infectious.
1−f is obviously the fraction of symptomatic infectious that become unreported symptomatic infectious. So
ν2=(1−f)ν is the rate at which asymptomatic infectious become unreported symptomatic infectious.
η is the rate at which the symptomatic infectious individuals lose their symptoms by death or healing. They get
1η the average time symptomatic infectious have symptoms.
Pierre Magal, et al. calculated the parameters of their model from the Wuhan epidemic data before public confinement and isolation measures.
A nonlinear adjustment method allowed them to obtain an expression of the CR function giving the cumulative number of reported symptomatic infectious individual at time t. They obtain that CR(t)=CR(t)=χ1exp(χ2t)−χ3.
As CR(t0)=0, they deduce t0=1χ2[ln(χ3)−ln(χ1)].
For experimental reasons, they assume that the average time during which asymptomatic infectious are asymptomatic and the average time symptomatic infectious have symptoms are around 7 days. They get 1ν=1η=17. From their data, they assume that f is in the interval [0.8 ; 1]. They use the value f=0.8 afterwards.
From all these results, they calculated the values of all the parameters corresponding to the epidemic of COVID-19 in Wuhan, notably the transmission rate and the basic reproductive number R0. Theoretically, they computed R0=τS0ν(1+ν2η). The data of Wuhan epidemic permit them to obtain τ=4.44×10−8, and R0=4.13 before the public measures.
With numerical simulations they show the evolution of the epidemic in Wuhan without public confinement and isolation measures. By varying the parameters of the model, they use simulations to predict the peak of the epidemic.
Since January 23, public confinement and isolation measures have been taken in Wuhan. Pierre Magal, et al. considered the effectiveness of its measurements after January 25. They therefore modified their model with the following non constant transmission rate:
τ(t)={4.44×10−8,0,if t∈[t0;25]if t>25
With these public measures, numerical simulations show that the epidemic in Wuhan will be eradicated. We quote them: "As a consequence of our study, we note that public health measures, such as isolation, quarantine, and public closings, greatly reduce the final size of this epidemic, and make the turning point much earlier than without these measures."
In our work, we use the first deterministic model of Pierre Magal, et al. (1.1) because we consider that it could reflect the dynamics of the epidemic in the Greater Abidjan Region as we will see later with the data of this region. We rephrase some starting notations. We study the system from an instant which we will define later. At this instant t=0 we note (S0,I0,R0,U0) the state of the system. N is the total population size before the epidemic. In this model, there is no birth or recruitment. So, we have S0+I0+R0+U0≤N.
In section 1, we prove that the solution of the system (1.1) is bounded in Δ={(x1,x2,x3,x4)∈R4+;x1+x2+x3+x4≤N}.
Obviously M={Ea=(a;0;0;0);0≤a≤N} is the set of the different disease-free equilibrium for the system (1.2). By using the LaSalle's global invariance principle [2,3], we show that all the solutions of (1.3), bounded for t≥0 in Δ, converge in M when t→+∞. These mathematical results comfort us in the possibility of reaching a state without the disease.
In Côte d'Ivoire, the Greater Abidjan Region includes the city of Abidjan and six other districts which are peripheral.
In Côte d'Ivoire, the first patient of COVID-19 was detected in Abidjan on March 11, 2020. At the start of the epidemic, the rest of the country was unaffected by the disease. So the public authorities took the following measures [4].
• The suspension of air flights to countries with more than 100 reported cases (decided on March 20, 2020).
• The closure of all land, sea and air borders from midnight on March 22, 2020 until further notice.
• On March 23, 2020 the state of emergency is established. A curfew is in place from 9 p.m. to 5 a.m. for two weeks.
• From March 29, 2020, no transport is authorized between the Greater Abidjan region and the rest of the country.
• On April 9, 2020, it was decided: the wearing of the mandatory mask in the Greater Abidjan; compulsory confinement of fragile people; reduction of non-essential travels and the number of people in vehicles.
The Greater Abidjan Region is the epicenter of the epidemic with 98 percent of the confirmed cases. Since March 29 the big Abidjan has been isolated from the rest of the country. The rest of the country has no more than 5 reported cases. Since at least April 21, 2020, the rest of the country has no COVID-19 patients. National data in [4] can therefore be assimilated to that of the Big Abidjan (Table 1, Table 2, Table 3, Table 4, Table 5 and Table 6).
The majority of the economic activities of the Ivorian population are activities of the informal sector (small trade, rustic restaurant, street vendors, illegal taxi drivers,...). These activities mainly focus on the city of Abidjan. This did not allow the establishment of total containment in Côte d'Ivoire. In addition to conventional barrier measures (washing hands regularly, coughing in the elbow...), the measures above have been gradually implemented to curb the epidemic. As our study is limited to the Greater Abidjan Region, we will analyze the data from the isolation of the Greater Abidjan Region (March 30, 2020). We will therefore study our model with the eradication measures implemented.
On the other hand, from this isolation date, the epidemic data are more reliable because of the multiplication of screening centers.
These measures, which do not correspond to total containment, favor the transmission of the virus. The transmission rate of the Greater Abidjan Region model is therefore not zero.
In section 2, as was done in [1,5,6], we will determine the constants of the system (1.1) corresponding to the epidemic in the Great Abidjan Region with the effect of the public measures. We consider the period of the effect of the public measures for two reasons: daily statistical data is more significant during this period because of the multiplication of screening centers and we want to know the behavior of the epidemic with these measures. We use the software CurveExpert Pro to fit our data. That permits us to compute the different parameters of the Greater Abidjan Region model and the basic reproductive number . Simulation with MatLab software permits us to find the peak of the epidemic in the Greater Abidjan Region in the case of this deterministic model. Simulations with other parameter values make it possible to observe a more rapid eradication of the epidemic and to propose measures to achieve this.
In Côte d'Ivoire, continuously, certain protective measures such as minimum physical distance of one meter, regular hand washing etc. are not respected in markets and public transport vehicles. Abidjan has many large markets in each of its districts. Our region has many uncontrolled urban transport that does not allow compliance with health measures. Also, the practice of economic activities in the informal sector promotes many contacts between individuals. All these factors continually disturb the transmission rate of the model and could be assimilated to white noise to be introduced into the system (1.1). That's why we add a white noise to the model (1.1) to disrupt the transmission rate.
The isolation of the Greater Abidjan Region is not fully respected despite the efforts of the police. Sudden inflows or outflows of people take place formally with administrative authorizations or illegally by village roads.
Sudden mid-term measures such as the obligation to wear a mask and a curfew extension caused sudden shocks to the model at the first time when they were applied.
The model described by the system (1.1) does not take these factors into account. We need to introduce jumps in the system (1.1) to represent these sudden changes which disturb the transmission rate.
As in [7-9] and from the model (1.1) we develop a model with jumps perturbation that disrupt the disease transmission rate. We get the system. (1.2)
⎧⎩⎨⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪dS(t)=[−τS(t)[I(t)+U(t)]]dt−σS(t)[I(t)+U(t)]dB(t) −∫ZC(z)S(t)[I(t)+U(t)]N˜(dt,dz)dI(t)=[τS(t)[I(t)+U(t)]−νI(t)]dt+σS(t)[I(t)+U(t)]dB(t) +∫ZC(z)S(t)[I(t)+U(t)]N˜(dt,dz)dR(t)=[ν1I(t)−ηR(t)]dtdU(t)=[ν2I(t)−ηU(t))]dt (1.2)
where B is a Brownian motion defined on a stochastic basis (Ω,F,(Ft≥0),P) and σ is constant.
X(t-) means the left limit of X(t), N˜(dt,dz) is a Poisson counting measure with the stationary compensator Π(dz)dt and Π is defined on a measurable subset Z of [0,∞) with Π(Z) < ∞.
In Section 3, we prove that the system (1.2) has a unique global and positive solution. We show how the jumps influence the behavior of the solution around the set of the equilibrium states. Inspired by [10,11] and using MaLab R2019b software, we perform numerical simulations which corroborate our theoretical results.
Here we justify simple mathematical properties that confer biological significance to our model.
N is the size of the population before the epidemic. We have S0+I0+R0+U0≤N. We prove that for any initial state in Δ={(x1,x2,x3,x4)∈R4+;x1+x2+x3+x4≤N}, the solution of the system (1.1) are bounded in Δ={(x1,x2,x3,x4)∈R4+;x1+x2+x3+x4≤N}.
For any given initial value in Δ={(x1,x2,x3,x4)∈R4+;x1+x2+x3+x4≤N}, the set Δ is positively invariant for the system (1.1).
With the system (1.1), we obtain by summing the equation of the system:
d[S(t)+I(t)+R(t)+U(t)]dt =−η[R(t)+U(t)]
If we prove that for t≥t0, R(t)+U(t)≥0 then t↦S(t)+I(t)+R(t)+U(t) will be decreasing and therefore less than or equal to that N.
On the other hand, initiating the solution of system (1.1) in positive orthant R4+, we have :
dS(t)dt∣∣S=0=0, dI(t)dt∣∣I=0=τS(t)U(t)≥0, dR(t)dt∣∣R=0=ν1I(t)≥0 and dU(t)dt∣∣U=0=ν1I(t)≥0.
So the solutions of system (1.1) remain in R4+ at any time t, i.e. S, I, R and U are positive. We thus obtain the desired result.
Starting from any initial condition, the following result guarantees us, under certain conditions, the convergence of any solution in M={Ea=(a;0;0;0);0≤a≤N} the set of the disease-free equilibrium for the system (1.1).
If R0≤1, all the solutions of (1.1) bounded for t≥0 in Δ, converge in M the set of all the disease-free equilibrium for the system (1.1) when t→+∞.
We define V:Δ→R by V(S,I,R,U)=S+I+τS0R+τS0ηU.
V˙(S,I,R,U)=−τS(I+U)+τS(I+U)−νI +τS0ν1I−ητS0R+τS0ην2I−τS0U
=−τSI−τSU+τSI+τSU−νI+τS0ν1I−τS0ηR+τS0ην2I−τS0U
=−τSU+τ(S−S0)U+(τS0ν1+τS0ν2η−ν)I−ητS0R
=−τS0U+(τS0ν1ν+τS0ν2ην−1)νI−ητS0R
=−τS0U+[τS0ν(ν1+ν2η)−1]νI−ητS0R.
Since all the model parameters are non-negative, S≤S0 and ν1≤1, we get :
τS0ν(ν1+ν2η) ≤τS0ν(1+ν2η)=R0 ≤1.
So we have τS0ν(ν1+ν2η)−1≤0. Since all the model parameters are non-negative and S≤S0, we deduce that V˙(S,I,R,U)≤0.
Moreover,
V˙(S,I,R,U)=0
⇔−τS0U+[τS0ν(ν1+ν2η)−1]νI−ητS0R=0
⇔(I,R,U)=(0,0,0)
⇔(S,I,R,U)=(a,0,0,0);0≤a≤N
⇔(S,I,R,U)∈M.
The only compact invariant subset of {x∈Δ;V˙=0} is M. So, by Lasalle's global invariance principle, all the solutions (1.1), bounded for t≥0 in Δ, converge in M when t→+∞.
In practical terms, when R0≤1, any solution of the system (1.1) converges towards a disease-free state.
Determination of the constants in the model of the COVID-19 epidemic in the Greater Abidjan Region with the effect of public measures
In this section, we will try to determine the parameters of the deterministic model corresponding to the dynamics of the epidemic in the Greater Abidjan Region after the implementation of all the public measures that led to its isolation.
As in [1] :
• 1ν, the average time during which asymptomatic infectious are asymptomatic is equal to 7 days and ν=17.
• 1η, the average time during which symptomatic infectious have symptoms is equal to 7 days and η=17.
All this brings us to consider our data 14 days after the isolation of the Greater Abidjan Region to tend towards the dynamics of the epidemic after all the public measures have been taken.
The initial instant t=0 of our study of the epidemic in the big Abidjan with the effect of public measures therefore corresponds to April 12, 2020. In our case we will name (S0,I0,R0,U0) the state of the system at this time. This is the initial state of the model that we are going to study.
We will not use the value of f indicated in [1] by Pierre Magal, et al. We will consider that f=0.05. It is assumed in [12] and [13] that the young age of an individual favors the absence of symptoms or the development of very weak symptoms after COVID-19 contraction. The last population census of Côte d'Ivoire [14] shows that children (0-14 years-old) represent 41.8 percent of the total population and young people (15-34 years-old completed) constitute 35.5 percent of the total population. 77.3 percent of the total population is under 35 years of age. This extreme youth of the population could explain the low number of symptomatic infectious because the symptoms of the disease rarely appear. On the other hand, the low attendance of health facilities (around 47 percent, [15]) may not favor the detection of COVID-19. Screening is only systematic in the event of serious symptoms and symptomatic infectious are hardly detected. There are not enough screening centers, although their number has been increased. For all these reasons, we will assume that around 5 percent of infectious individuals become confirmed cases. f=0.05 and 1−f=0.95.
The parameters ν=17, η=17, f=0.05 and therefore ν1=fν=0.057 and ν2=(1−f)ν=0.957 are thus known.
N=5816060 is the size of the population of the big Abidjan before the epidemic.
We know the number of reported cases at t=0 (April 12, 2020) : R0=R(0)=574−533=41.
As in [1], we define the CR function of the cumulative number of reported symptomatic infectious cases at time t for t≥0 with the function
CR(t)=ν1∫tt0I(s)ds.
where t0 is the solution of CR(t)=0.
The statistical data of the COVID-19 epidemic show that the variations in the number of susceptible people due to the epidemic in the beginning of the epidemic are so small so we can consider S(t0)=N and S0=S(0)=5816060 because as we will see later the distance between 0 and t0 is small. That will permit us to evaluate the other parameters of the model. We can later adjust the value of S0.
As we look at the epidemic from the effect of the measures at date t = 0 then, for t0≤t < 0, CR(t) exists but has no interesting biological significance and is used just for calculations (Figure 1).
To exploit this adjustment function, we consider that For t≥0, CR(t)=a×exp(bt)−c where the values of a, b and c are defined above.
When we use the methods to estimate the parameters presented and justified in [1] in section 4 we obtain that :
t0=1b(ln(c)−ln(a)),
I(0)=cbfν×expb(−t0),
τ=b+νS0×η+bν1+η+b,
U(0)=(1−f)νη+bcbfν×expb(−t0),
and
R0=τS0ν(1+ν2η).
From the previous expressions and the determination of the previous parameters, we calculate the others constants of the model of the COVID-19 epidemic in the great Abidjan with the effect of public measures. We get t0=−15.146, I0=I(0)=5373.7, U0=U(0)=4992.8, τ=1.3019×10−8 and R0=1.0335.
With MatLab software, from the date of (April 12, 2020, t=0), the effect of all the public measures, we propose a simulation of the current situation in the Greater Abidjan Region with the constants previously determined. Let's recap.
The transmission rate is τ=1.3019×10−8. The rate at which the asymptomatic infectious individuals become symptomatic is ν=17. f=0.05 is the fraction of asymptomatic infectious that become reported symptomatic infectious. So ν1=fν=0.057.
The rate at which asymptomatic infectious become unreported symptomatic infectious is
ν2=(1−f)ν=0.957
η=17 is the rate at which the symptomatic infectious individuals lose their symptoms by death or healing. The basic reproductive number for our deterministic model is R0=1.0335.
At t=0, S0=5816060, I0=5373.7, U0=4992.8, τ=1.3019×10−8 and R0=41. We get the following numerical simulations with MatLab sofware where we will not present the trajectory of the susceptible people because of a question of scale (Figure 2).
Now we present the evolution of susceptible individuals (Figure 3).
Under the conditions we have adopted, the epidemic in the Greater Abidjan lives around 1400 days after April 12, 2020, with a peak reaching 230 days after this date.
We can see that at the end of the epidemic, the number of susceptible people in the population of the big Abidjan is around 5,300,000. This drop should not be seen as deaths but rather as cures because the death rate from the disease is very weak (it is less than 1 percent as the statistics indicate. Refer also to the linear graphs). At the end of the epidemic, we have approximately 300,000 cured which does not influence the dynamics of the model.
In support of the public measures taken, we agree with the two new directions taken by the government of Côte d'Ivoire:
• Increase in screening tests,
• Repression and punishment of people who do not respect public measures and protective gestures.
An increase in screening tests results in an increase in the parameter f and therefore in the rate ν1.
We suggest a large increase from f to 0.8. We hope that 80 percent of the infected will be reported cases in order to get them out of the dynamics of the epidemic.
The degree of repression and sanction can result in a coefficient k between 0 and 1 by which the current transmission rate is multiplied. If the repression and the sanctions are high then k is small and the rate of transmission decreases. For our simulation, we opt for k = 0.5, i.e. a reduction by half of the current transmission rate. With the previous values of the other parameters, using the same methods of parameter approximations the following simulation gives R0=0.5131 (Figure 4).
The following simulations shows in the same conditions the evolution of the number of individuals susceptible which knows a small decrease to stabilize at 5,815,600 at the end of the epidemic (Figure 5).
In this section, we study the model (1.2) for where represents the date of April 12, 2020, the date from which we assume that the public measures are effective.
Lyapunov analysis [7,8,16], permits us to prove that the solution of model (1.2) is positive and global.
We assume that the jump coefficient satisfies the following conditions:
\begin{itemize}
• (H1) ∫Z|H(x,z)−H(y,z)∣∣∣∣2Π(dz)≤Lm∣∣∣∣x−y∣∣∣∣2 where
H(x,z)=C(z)S(x−)[I(x−)+U(x−)], with |x|∨|y|≤m
• (H2) C is a bounded function and there exists a strictly positive constant K1 such that C(z) < K1.
We will first prove that the system (1.2) has a unique solution which is global and positive. In a second point we will study the asymptotic behavior of this solution in the case R0≤1. Simulations carried out with Matlab software confirm our results.
As in [16], we show that jump processes can suppress the explosion.
Theorem 6.1: With the assumptions (H1) and (H2), for any given initial value (S(0),I(0),R(0),U(0)∈Δ={(x1,x2,x3,x4)∈R4+;x1+x2+x3+x4≤N}, the equation (1.2) has a unique global solution (S(t),I(t),R(t),U(t))∈Δ for all t≥0 almost surely.
Proof: For any given initial value (S(0),I(0),R(0),U(0))∈Δ there is a unique local solution (S,I,R,U) on a random interval [0,τe[ where τe is the explosion time because of (H1) and the fact the drift and diffusion are locally Lipchitz.
To prove this local solution is global, we need to show that τe=∞ a.s. We will also prove that for any given initial value ((S(0),I(0),R(0),U(0))∈Δ, the equation (1.2) has its solution
((S(t),I(t),R(t),U(t))∈Δ for all t≥0 almost surely.
At first we prove that the solution is positive. In second part, we prove that the solution checks S(t)+I(t)+R(t)+U(t))≤N for all t≥0 almost surely.
The set {(S(t),I(t),R(t),U(t)):t∈[0,τe[ with S(t)≤0 or I(t)≤0 or R(t)≤0 or U(t)≤0} is not empty cause at t=0, (S(0),I(0),R(0),U(0)) is in this set. So we define the following stopping time
τ=inf{t∈[0,τe[:S(t)≤0 or I(t)≤0 or R(t)≤0 or U(t)≤0}.
Let m0 > 0 be so large that max{|S(0)|,|I(0)|,|R(0)|,|U(0)|}≤m0.
For m≥m0, we define the stopping time
τm=inf{t∈[0,τ[:|S(t)|>m or |I(t)|>m or |R(t)|>m or |U(t)|>m}
Clearly, τm is increasing as m↑∞ a.s. Set τ∞=limm→∞τm, whence, τ∞≤τ≤τe a.s.
If we can show that τ∞=∞ a.s is true, then τe=∞=τ a.s then the solution (S(t),I(t),R(t),U(t)) is global and positive a.s.
For x=(x1;x2;x3;x4)∈R4+, we define F a C2-function by F:(x1;x2;x3;x4)↦x21+x22+x23+x24.
Let m≥m0 and T>0 be arbitrary. For t∈[0;τm∧T[,
F(S(t),I(t),R(t),U(t))=S2(t)+I2(t)+R2(t)+U2(t) is well defined.
Applying Itô's formula, we obtain
dF(S(t),I(t),R(t),U(t))=LF(S(t),I(t),R(t),U(t))dt
+σS(t)(I(t)+U(t))(S2(t)−2S(t)+2I(t)−I2(t))dB(t)
+∫ZC(z)S(t−)(I(t−)+U(t−))[2C(z)S(t−)(I(t−)+U(t−))+2I(t−)−2S(t−)]N˜(dt,dz)
where
LF(S(t),I(t),R(t),U(t))=2[τS(t)I(t)(I(t)+U(t)) +ν1I(t)R(t)+ν2I(t)U(t)+4σ2
−τS2(t)(I(t)+U(t))−νI2(t)−ηR2(t)−ηU2(t)]
+∫ZC(z)S(t)(I(t)+U(t))[2C(z)S(t)(I(t)+U(t))+S2(t)−I2(t)+4I(t)−4S(t)]dΠ(z).
For t∈[0;τm∧T[, S(t),I(t),R(t), live in $ [0;m]$ and σ is constant. As τ,ν,ν1,ν2 and η live in [0;1], we get
LF(S(t),I(t),R(t),U(t))≤2[τS(t)I(t)(I(t)+U(t)) +ν1I(t)R(t)+ν2I(t)U(t)+4σ2
+∫ZC(z)S(t−)(I(t−)+U(t−))[2C(z)S(t−)(I(t−)+U(t−))
+S2(t−)−I2(t−)+4I(t−)−4S(t−)]dΠ(z).
We are in conditions where S(t),I(t),R(t), τ,ν,ν1,ν2 and η are dominated by positive constants. So, there is a strictly positive real constant K2 such that
[τS(t)I(t)(I(t)+U(t))+ν1I(t)R(t)+ν2I(t)U(t)++4σ2≤K2.
In other words by using the assumption H2, the fact S(t),I(t),R(t), live in [0;m] for t∈ [0;τm∧T] and the fact
Π is defined on a measurable subset Z with Π(Z) < ∞, we find a strictly positive real constant K3 such that
∫ZC(z)S(t)(I(t)+U(t))[2C(z)S(t)(I(t)+U(t))+S2(t)−I2(t)+4I(t)−4S(t)]dΠ(z) < K3.
Finally
LF(S(t),I(t),R(t),U(t))≤K2+K3=K
Therefore,
∫0τm∧TdF(S(t),I(t),R(t),U(t)) ≤∫0τm∧TKdt
+∫0τm∧TσS(t)(I(t)+U(t))(S2(t)−2S(t)+2I(t)−I2(t))dB(t)
+∫0τm∧T∫ZC(z)S(t−)(I(t−)+U(t−))[2C(z)S(t−)(I(t−)
+U(t−))+2I(t−)−2S(t−)]N˜(dt,dz).
Taking expectation, yields
E(F(S(τm∧T),I(τm∧T),R(τm∧T),U(τm∧T))) ≤F(S(0),I(0),R(0),U(0))+KT=K′ (3.1)
We define for each u≥ 0,
μ(u)=inf{F(x):x=(x1,x2,x3,x4)∈[0;+∞[4with |x1|>u or |x2|>u or |x3|>u or |x4|>u}
Due to the property of the function $ F $ defined by F(x1,x2,x3,x4):=x21+x22+x23+x24, we see that
limu→+∞μ(u)=+∞
We know that if X is a Levy process then, for each t≥0 we have ΔX(t)=0 almost surely where
ΔX(t)=X(t)−X(t−)
S, R, I and $ R $ are Levy processes. So we obtain almost surely from (3.1) :
μ(m)P(τm≤T) ≤E(F(S(τm∧T),I(τm∧T),R(τm∧T),U(τm∧T))I{τm≤T})
≤E(F(S(τm∧T),I(τm∧T),R(τm∧T),U(τm∧T))
≤K′.
We deduce
P(τm≤T)≤K′μ(m).
Letting m→+∞ yields
P(τ∞≤T)=0
Since T is arbitrary, we must have
P(τ∞=+∞)=1 a.s.
So τe=τ=+∞ a.s. and the equation (1.2) has a unique global positive solution X(t) for t≥0.
Now we prove that the solution is dominated.
Let ω be fixed in Ω. t→S(ω)(t)+I(ω)(t)+R(ω)(t)+U(ω)(t) is a deterministic function such that
S(ω)(t)+I(ω)(t)+R(ω)(t)+U(ω)(t)dt=−η[R(ω)(t)+U(ω)(t)]
η is a positive constant. As we have shown above, U(t) and R(t) are positive for t≥0 we can deduce that for t≥0, the deterministic function t→S(ω)(t)+I(ω)(t)+R(ω)(t)+U(ω)(t) is a decreasing function.
Therefore for t≥0, S(ω)(t)+I(ω)(t)+R(ω)(t)+U(ω)(t)≤S(0)+I(0)+R(0)+U(0).
The initial state (S(0),I(0),R(0),U(0) being taken in
Δ={(x1,x2,x3,x4)∈R4+;x1+x2+x3+x4≤N}, therefore we have for
t≥0, S(ω)(t)+I(ω)(t)+R(ω)(t)+U(ω)(t)≤N.
By varying ω in the universe of probability Ω, we thus show that:
for t≥0, S(t)+I(t)+R(t)+U(t)≤N almost surely.
This completes our proof and for any given initial value
(S(0),I(0),R(0),U(0)∈Δ={(x1,x2,x3,x4)∈R4+;x1+x2+x3+x4≤N}, the equation (1.2) has a unique global solution (S(t),I(t),R(t),U(t))∈Δ for all t≥0 almost surely.
Asymptotic behavior around the set of the disease-free equilibrium for the system (1.1)
Considering that R0≤1, we prove that the solution of the stochastic system (1.2) oscillates in the set of the disease-free equilibrium for the deterministic system (1.1).
Theorem 6.2: Let (S(t),I(t),R(t),U(t)), be the solution of the system (1.2) with initial value. If R0≤1 , then there exists ξ < 0 such that
limsupt→∞ 1t E∫t0 [S(x)(I(x)+U(x))+U(x)+I(x)+R(x]dx ≤σ22N+∫Z12NC2(z)Π(dz)ξ
where ξ=min{τN5 ; τS0 ; ν(1−R0) ; ητS0}.
Proof: We know that (S(0),I(0),R(0),U(0))=(S0,I0,R0,U0)
Define a C2-function by
G(S,I,R,U)=12N5S2+S+I+τS0R+τS0ηU
The function G is positive definite on Δ, and
dG(S(t),I(t),R(t),U(t))=LGdt +σN5S2(I+U)dB(t) +∫Z{12N5C(z)S(t−)(I(t−)+U(t−))
×[C(z)S(t−)(I(t−)+U(t−))−2S(t−)−2N5]}N˜(dt,dz) (3.2)
where
LG=−τN5S(t)2(I(t)+U(t))−τS(t)(I(t)+U(t))
+τS(t)(I(t)+U(t))−νI(t)+τS0ν1I(t)−ητS0R(t)+τS0ην2I(t)−τS0U(t)
+σ22N5S2(t)(I( +∫Z[12N5C2(z)S2(t)(I(t)+U(t))2 −C(z)S(t)(I(t)+U(t))]Π(dz).
..
We obtain
LG=−τN5S(t)2(I(t)+U(t))) −τS(t)I(t)−τS(t)U(t)+τS(t)I(t)+τS(t)U(t)
−νI(t)+τS0ν1I(t)−τS0ηR(t) +τS0ην2I(t)−τS0U(t)(t) +σ22N5S2(t)(I(
+∫Z[12N5C2(z)S2(t)(I(t)+U(t))2−C(z)S(t)(I(t)+U(t))]Π(dz).
Thereafter,
LG=−τN5S(t)2(I(t)+U(t))
−τS(t)U(t)+τ(S(t)−S0)U(t)+(τS0ν1+τS0ν2η−ν)I(t)−ητS0R(t)
+σ22N5S2(t)(I(+∫Z[12N5C2(z)S2(t)(I(−C(z)S(t)(I(t)+U(t))]Π(dz),
LG =−τN5S(t)2(I(t)+U(t)))
−τS0U(t)+(τS0ν1ν+τS0ν2ην−1)νI(t)−ητS0R(t)
+σ22N5S2(t)(I(+∫Z[12N5C2(z)S2(t)(I(−C(z)S(t)(I(t)+U(t))]Π(dz),
and
LG =−τN5S(t)2(I(t)+U(t))
−τS0U(t)+[τS0ν(ν1+ν2η)−1]νI(t)−ητS0R(t)
+σ22N5S2(t)(I(+∫Z[12N5C2(z)S2(t)(I(−C(z)S(t)(I(t)+U(t))]Π(dz).
Since all the model parameters are non-negative, S≤S0 and ν1≤1, we get :
τS0ν(ν1+ν2η) ≤τS0ν(1+ν2η)=R0 ≤1.
Since all the model parameters are non-negative and S(t)≤S0, we deduce that
LG≤−τN5S(t)2(I(t)+U(t))
−τS0U(t)+(R0−1)νI(t)−ητS0R(t)
+σ22N5S2(t)(I(+∫Z[12N5C2(z)S2(t)(I(−C(z)S(t)(I(t)+U(t))]Π(dz).
Using that for t≥0, (S(t),I(t),R(t),U(t))∈Δ={(x1,x2,x3,x4)∈R4+;x1+x2+x3+x4≤N} and the fact that the function C is dominated by a positive real K1, there is a real T such that
∫Z−C(z)S(t)(I(t)+U(t))Π(dz) < T.
Hence,
LG≤−τN5S(t)(I(t)+U(t))−τS0U(t)+(R0−1)νI(t)−ητS0R(t)+T+σ22N+∫Z12NC2(z)Π(dz). (3.3)
Integrating both sides of (3.2) from 0 to t, then taking expectation, yields
0≤EG(S(t),I(t),R(t),U(t))=G(S(0),I(0),R(0),U(0))+E∫0tLG(S(x),I(x),R(x),U(x)]dx
By (3.3) we have
E[∫0t{τN5S(x)2(I(x)+U(x))+τS0U(x)+(1−R0)νI(x)+ητS0R(x)}dx]
≤G(S(0),I(0),R(0),U(0))+T+σ22N5S2(t)(I(+∫Z12N5C2(z)S2(t)(I(Π(dz).
Therefore
limsupt→∞1tE∫0t[S(x)2(I(x)+U(x))+U(x)+I(x)+R(x]dx≤σ22N+∫Z12NC2(z)Π(dz)ξ
Where ξ=min{τN5 ; τS0 ; ν(1−R0) ; ητS0}
Remark 6.2: When R0≤1 the solution of (1.2) oscillates more closely in the set M of the disease-free equilibrium for the system (1.1) as the intensity of the noise and the jumps decrease. The case fatality rate is less than 1 percent (see linear graphs). Convergence towards a disease-free state is an encouraging sign. It is a disease-free situation with few deaths.
As in [7,10] using MatLab R2019b software, we add white noise and Poisson jumps to the first example seen in the deterministic case which reflects the reality in the big Abidjan. To be close to the difficulties of lowering the transmission rate in greater Abidjan, we lowered this rate by 10 percent in the first two figures and by 25 percent in the last two. We get basic reproductive numbers less than 1.
We obtain the following figures which corroborate the results of our last theorem (Figure 6).
Now we present the evolution of susceptible individuals (Figure 7 and Figure 8).
Now we present the evolution of susceptible individuals (Figure 9).
In this study, we tried to describe the behavior of the COVID-19 epidemic in the great Abidjan during the period of public measures. We use a stochastic model with white noise and jumps which reflects the disturbances that hamper the implementation of these measures. Total containment which would reduce the transmission rate to zero is impossible in the current context. The application of controls and actions restriction, would greatly reduce the transmission rate τ. Massive screening corresponds to a high value of the parameter f. This means a large value of ν1 and a small value of ν2. This permits to reduce the propagation of the epidemic because this allows a large number of infectious individuals to escape from the dynamics of the epidemic. Therefore, the basic reproductive number R0=τS0ν(1+ν2η) which directly depend on τ and ν2 decreases.
These two factors combined (actions restriction and massive screening) allow a large reduction of R0=τS0ν(1+ν2η), the basic reproductive number of the deterministic model and a rapid drop of the size of the epidemic, in the stochastic model with jumps. The case fatality rate is less than 1 percent (see linear graphs), and the solution of our stochastic model together flow in the set of the disease-free states. Consequently, we tend to get out of the epidemic with few deaths.