# Transmission Model of Hepatitis B Virus with the Migration Effect

Authors: Muhammad Altaf Khan; Saeed Islam; Muhammad Arif; Zahoor ul Haq
Journal: BioMed Research International
 (2013)

---

## Abstract

Hepatitis B is a globally infectious disease. Mathematical modeling of HBV transmission is an interesting research area. In this paper, we present characteristics of HBV virus transmission in the form of a mathematical model. We analyzed the effect of immigrants in the model to study the effect of immigrants for the host population. We added the following flow parameters: “the transmission between migrated and exposed class” and “the transmission between migrated and acute class.” With these new features, we obtained a compartment model of six differential equations. First, we find the basic threshold quantity Ro and then find the local asymptotic stability of disease-free equilibrium and endemic equilibrium. Furthermore, we find the global stability of the disease-free and endemic equilibria. Previous similar publications have not added the kind of information about the numerical results of the model. In our case, from numerical simulation, a detailed discussion of the parameters and their numerical results is presented. We claim that with these assumptions and by adding the migrated class, the model informs policy for governments, to be aware of the immigrants and subject them to tests about the disease status. Immigrants for short visits and students should be subjected to tests to reduce the number of immigrants with disease.


---

## Body

## 1. Introduction

According to World Health Organization, about 350 million people are infected with the hepatitis B virus (HBV)1 and about 170 million people are chronically infected with the hepatitis C virus (HCV)2. The majority of those infected live in developing countries with few incidences in Western countries. Migration is one of the defining issues of our time. For example, more than 5 million Canadians migrate out of the country each year and over 250,000 new immigrants arrive in Canada each year. Countries of the world are increasingly connected through travel and migration, and thus, migration has health implications in one location for both local and global migrations, since infectious diseases do not remain isolated geographically. The UK Hepatitis Foundation estimated in 2007 that the number of hepatitis B cases in the UK doubled in the previous 6 years chiefly due to immigration of infected people, many from the new member states of the European Union where the prevalence of viral hepatitis is higher. The United Arab Emirates mandated that hepatitis C testing is to be done at the time of residence visa renewals and added hepatitis C to the list of diseases such as HIV and hepatitis B as diseases warranting deportation [1]. In China, hepatitis B virus infection is a major public health problem. Hepatitis B is the first one among the diseases with legal management measures. In China, an estimated 93 million people have been infected with the hepatitis B virus [2]. The seroepidemiological survey on HBV infection conducted in 2006 showed that HBsAg carrier rate was 7.18 percent in the overall dynamics of HBV. In this paper, we consider a system of ordinary differential equations which describes the transmission of HBV transmission in China. Several mathematical models have been formulated on the HBV transmission in China. Medley and coauthors used a mathematical model and developed the strategies to eliminate the HBV in New Zealand in 2008 [3, 4]. Anderson and May [5] used a mathematical model which illustrated the effects of carriers on the transmission of HBV. An age structure model was proposed by Zhao et al. [6] to predict the dynamics of HBV transmission and evaluate the long-term effectiveness of the vaccination program in China. Wang et al. [7] proposed and analyzed the hepatitis B virus infection in a diffusion model confined to a finite domain. A hepatitis B virus (HBV) model with spatial diffusion and saturation response of the infection rate is investigated by Xu and Ma [8]. Also, Zou et al. [9] proposed a mathematical model to understand the transmission dynamics and prevalence of HBV in mainland China. A model to describe waning of immunity after sometime has been studied by a number of authors [10–16]. In the context of rapid global migration, there is a potential for any disease to be transferred faster than was previously possible. These implications concerning the movement of HBV and HCV merit far more attention by countries and the international community than they have given the problem to date. This is especially important given that the scope and speed of migration is expected to grow in coming years. In the context of rapid global migration, there is a potential for any disease to be transferred faster than was previously possible. These implications concerning the movement of HBV and HCV merit far more attention by countries and the international community than they have given the problem to date. This is especially important given that the scope and speed of migration is expected to grow in coming years.In this paper, we construct the compartmental model of hepatitis B transmission. We have categorized the model into six compartments: Susceptible-S(t), Exposed-E(t), Acute-A(t), Carrier-C(t), Vaccinated-V(t), and Migrated-M(t) individuals. The migrated class of individuals comes from different parts of the world to the host country, and their interaction occurs in the form of sexual interactions, blood transportation, and transfusion. We modify the model from Pang et al. [10] by adding some new transmission dynamics and introduce the migrated class in the model. Furthermore, some authors [11, 13, 15, 17] show that acute hepatitis B could be found today in newborns of infected mothers. Pang et al. [10] added the vertical transmission term to exposed class from chronic carriers class on the basis of the characteristics of HBV transmission. In this paper, we improved the model of [10], with these new features, by adding the migrated class M(t) and the following parameters: (i)
μ1: the transmission rate from migrated class to exposed class,(ii)
μ2: the transmission rate from migrated to acute class,(iii)
δ: the death rate at the migrated class. The paper is organized as follows. Section 2 is devoted to the mathematical formation of the model. In Section 3 we find the Basic Reproduction Number, the disease-free equilibrium, and endemic equilibrium of the proposed model. Local asymptotic stability of the disease-free equilibrium and endemic equilibrium is discussed in Section 4. In Section 5, we study the global asymptotic stability of disease-free and endemic equilibria using the Lyapunov function. In Section 6, we study the numerical results of the proposed model and present the results in the form of plots for illustrations. The conclusion and references are presented in Section 7.


## 2. Model Formulation

In this section, we present the mathematical formulation of the compartmental model of hepatitis B, which consists of a system of differential equations. The model is based on the characteristics of HBV transmission. We divide the total population into six compartments, that is, Susceptible individualsS(t), exposed E(t), Acute A(t), Carrier C(t), Immunity class V(t), and Migrated M(t). The flow diagram (Figure 1) and the system is given in the following:
(1)dS(t)dt=δπ(1-ηC(t))-δS(t)-β(A(t)+κC(t))S(t)+δoV(t)-pS(t),dE(t)dt=β(A(t)+κC(t))S(t)-δE(t)+δπηC(t)-γ1E(t)+μ1M(t),dA(t)dt=γ1E(t)-δA(t)-qγ2A(t)-(1-q)γ1A(t)+μ2M(t),dC(t)dt=qγ2A(t)-δC(t)-γ3C(t),dV(t)dt=γ3C(t)+(1-q)γ1A(t)-δoV(t)-δV(t)+δ(1-π)+pS(t),dM(t)dt=-(μ1+μ2)M(t)-δM(t).
Subject to the initial conditions,
(2)S(0)≥0,E(0)≥0,A(0)≥0,C(0)≥0,V(0)≥0,M(0)≥0.Figure 1
The complete flow diagram of hepatitis B virus transmission model.The proportion of failure immunization is shown byπ. δ represent both the death and birth rate. At the γ1 rate the exposed individuals become infectious and move to the Acute class. γ2 is the rate at which the individuals move to the carrier class, γ3 is the flow of carrier to vaccinated class, β shows the transmission coefficient, κ represents the carrier infectiousness to acute infection, q is the proportion of acute individuals that become carrier, δo represent the loss of immunity rate and the individual become the susceptible again, p represents the vaccination of susceptible individuals, μ1 represents the rate of flow from migrated class to exposed class, and γ2 is the rate of transmission from migrated class to acute class. η is the unimmunized children born to carrier mothers, δ(1-π) measures the successful immunization of newborn babies, and the term δπ(1-ηC(t)) shows that the newborns are unimmunized and become susceptible again.We assume that the total population size is equal to 1, and just for simplifications,S(t) is the susceptible, E(t) the exposed, A(t) the acute, C(t) the carrier, V(t) the immunity, and M(t) is the migrated class representing the state variables in our proposed population model. The sum of the total population is S(t)+E(t)+A(t)+C(t)+V(t)+M(t)=1, holds. We just add (1), and we can easily get. We ignore the fifth equation in system (1); so, the new models become
(3)dS(t)dt=δπ(1-ηC(t))-δS(t)-β(A(t)+κC(t))S(t)dS(t)dt=+δo(1-(S(t)+E(t)+A(t)+C(t)+M(t)))-pS(t),dE(t)dt=β(A(t)+κC(t))S(t)-δE(t)+δπηC(t)-γ1E(t)+μ1M(t),dA(t)dt=γ1E(t)-δA(t)-qγ2A(t)-(1-q)γ1A(t)+μ2M(t),dC(t)dt=qγ2A(t)-δC(t)-γ3C(t),dM(t)dt=-(μ1+μ2)M(t)-δM(t).
Let
(4)Γ={(S,E,A,C,M)∈ℜ+5,∣S(t)≤δπ+δoδ+δo+p,(S+E+A+C+M≤δπ+δoδ+δo)}.
Here, Γ is a positively invariant set. All the solutions lie inside Γ which is our main focus.


## 3. Basic Reproduction Number/Threshold Quantity

The basic reproduction number or the threshold quantityℜ0 for the proposed model gives an average number of secondary infection when an infection is introduced in a purely susceptible population. We use the idea developed by [18], and also for detail see [19]. We have(5)ℱ=[0βSoβκSo000000],V=[δ+γ10-δπη-γ1δ+qγ2+(1-q)γ100-qγ2δ+γ3],V-1=[(δ+qγ2+(1-q)γ1)(δ+γ3)qγ2δπηδπη(δ+qγ2+(1-q)γ1)γ1(δ+γ3)(δ+γ1)(δ+γ3)δπηγ1qγ1γ2qγ2(δ+γ1)(δ+γ1)(δ+qγ2+(1-q)γ1)](δ+γ1)(δ+qγ2+(1-q)γ1)(δ+γ3)-δπηqγ1γ2,where Q1=(δ+qγ2+(1-q)γ1) and(6)FV-1=[0βSoδπηqγ2(δ+γ1)(δ+qγ2+(1-q)γ1)(δ+γ3)-δπηqr1r2κβSo(δ+γ1(1-q)+qγ2)(δ+γ1)(δ+qγ2+(1-q)γ1)(δ+γ3)-δπηqr1r2000000].So, the reproduction number given by ρ(FV-1) is
(7)ℜo=κβSo(δ+γ1(1-q)+qγ2)+δπηqr1r2(δ+γ1)(δ+qγ2+(1-q)γ1)(δ+γ3).
Here, So shows the disease-free equilibrium (DFE), and Do=(So,0,0,0,0), giving So=(δπ+δo)/(δ+δo+p). The endemic equilibrium point T*=(S*,E*,A*,C*,M*) for system (3), whose endemic equilibrium is given in the following subsection.Endemic Equilibria. To find the endemic equilibria of the system (3), by setting S=S*, E=E*, A=A*, C=C*, and M=M*, equating left side of the system (3) equal to zero, we obtained
(8)S*=Q1(δ+γ1)(δ+γ3)(1-R0)+κβSoQ1γ1β((δ+γ3)+κqγ2)-(δ+γ3)[(δ+γ1)+μ1γ1]M*γ1β((δ+γ3)+κqγ2)A*,E*=Q1A*-μ2M*γ1,C*=qγ2(δ+γ3)A*.


## 4. Local Stability Analysis

In this section, we find the local stability of disease-free and endemic equilibria. First, we show the local asymptotical stability of DFE equilibrium, and then we find the local asymptotical stability of endemic equilibrium. Now, we show the local stability of DFE about the pointDo=(So,0,0,0,0) in the following theorem.Theorem 1.
ForR0≤1, the disease-free equilibrium of the system (3) about an equilibrium point Do=(So,0,0,0,0) is locally asymptotically stable if ((Q1(δ+δo+p)(δ+γ1)>βγ1(δπ+δo)); otherwise, the disease-free equilibrium of the system (3) is unstable for R0>1.Proof.
To show the local stability of the system (3), about the point Do, we set the left-hand side of the system (3) equating to zero, and we obtain the following Jacobian matrix Jo(ζ):(9)Jo(ζ)=(-δ-δo-pδo-βSo-δo-δπη-βκSo-δo00-(δ+γ1)βSoβκSo+δπημ10γ1-Q10μ200qγ2-(δ+γ3)00000-(μ1+μ2+δ)).By the elementary row operation, we get the following matrix:(10)Jo(ζ)=(-(δ+δo+p)-δo-βSo-δo-δπη-βκSo00-(δ+γ1)βSoβκSo+δπημ100-Q1(δ+γ1)+γ1βSoγ1(βκSo+δπη)(δ+γ1)μ2+μ1γ1000T100000-(μ1+μ2+δ)),where Q1=(δ+qγ2+(1-q)γ1) and T1=-(-Q1(δ+γ1)+γ1βSo)-qγ1γ2(βκSo+δπη). The characteristic equation to the previous Jacobian matrix is given by
(11)(λ+δ+δo+p)(λ+μ1+μ2+δ)(λ3+a1λ2+a2λ+a3)=0.
The first two eigenvalues -(δ+δo+p) and -(μ1+μ2+δ) have negative real parts. For the rest of the eigenvalues, we get
(12)a1=(δ+δo+p)(Q1(δ+γ3)+1+Q1)(1-R0)+β(δπ+δo)((γ1(δ+γ3)+qγ1γ2)-(Q1κ+γ1))×(δ+δo+p)-1,a2=(Q1(δ+γ1)(δ+δo+p)-βγ1(δπ+δo))×{(δ+δo+p)[(δ+γ1)+(δ+γ1)2]+(δ+γ1)[Q1(δ+γ1)(δ+δo+p)-βγ1(δπ+δo)]+qγ1γ2(δ+γ1)[βκ(δπ+δo)+δπη(δ+δ+p)][(δ+γ1)+(δ+γ1)2]}×(δ+δo+p)-2+{qγ1γ2(δ+γ1)(βκ(δπ+δo)+δπη(δ+δo+p))}×(δ+δo+p)-1,a3=(δ+γ1)(δ+δo+p)(Q1(δ+δo+p)(δ+γ1)-βγ1(δπ+δo))×[qγ1γ2(βκ(δπ+δo)+δπη)+(Q1(δ+γ1)(δ+δo+p)-βγ1(δπ+δo))].
By the Routh-Hurwitz criteria, a1>0, a3>0, and a1a2>a3. Here, a1>0 when R0≤1 and ((Q1(δ+δo+p)(δ+γ1)>βγ1(δπ+δo)). Also, a2>0 and a3>0, and then a1a2>a3. So, according to the Routh-Hurwitz criteria, the Jacobian matrix has negative real parts if and only if R0≤1. Thus by Routh-Hurwitz criteria, the DFE of the system (3) is locally asymptotically stable about the point Do=(So,0,0,0,0). The proof is completed.The stability of the disease-free equilibrium forR0≤1 means that the disease dies out from the population. Next we show that the endemic equilibrium of the system (3) is locally asymptotically stable for R0>1. When the disease-free equilibrium is locally asymptotically stable for R0≤1, then the endemic equilibrium does not exist, but we are interested to know about the properties of an endemic equilibrium for R0>1.


### 4.1. Stability of Endemic Equilibrium (EE)

In this subsection, we find the local asymptotic stability of EE aboutD*=(S*,E*,A*,C*,M*), and we prove the local stability of endemic equilibrium in the following.Theorem 2.
ForR0>1, the endemic equilibrium D* of system (3) is locally asymptotical stable, if the following conditions hold:
(13)γ1Z2δoS*2>γ1β2κ,G1>G2;
otherwise, the system is unstable.Proof.
Here, we prove the that the system (3) about the equilibrium point D* is locally stable, and for this, we obtain the Jacobian matrix J*(ζ) of the system (3), in the following:(14)J*(ζ)=(-(δ+δo+p+β(A*+κC*))δo-βS*-δo-δπη-βκS*-δo-δoβ(A*+κC*)-(δ+γ1)βS*βκS*+δπημ10γ1-Q10μ200qγ2-(δ+γ3)00000-(δ+μ1+μ2)).By elementary row operation and after simplification, we get the following Jacobian matrix:(15)J*(ζ)=(-Z1-δo-βS*-δoZ3-δo0-Z1(δ+γ1)-Z2δoZ1βS*-Z2(βS*+δo)Z1Z4+Z2Z3Z1μ1-Z2δo00Z5γ1(Z1Z4+Z2Z3)Z6000Z700000-(μ1+μ2+δ)),where
(16)Z1=δ+δo+p+β(A*+κC*),Z2=β(A*+κC*),Z3=-(δπη+βκS*+δo),Z4=(βκS*+δπη),Z5=-Q1(-Z1(δ+γ1)-Z2δo)+γ1(Z4βS*-Z2(βS*+δo)),Z6=-μ2(-Z1(δ+γ1)-Z2δo)+γ1(Z1μ1-Z2δo),Z7=-Z5(δ+γ3)-Z5γ1(Z1Z4+Z2Z3).
The eigenvalue λ1=-(μ1+μ2+δ)<0, λ2=-Z1, and using the value of Z1 and further C*, we get λ2=-(δ+δo+p+β((δ+γ3)+κqγ2)A*)<0. λ3=-(Z1(δ+γ1)+Z2δo)<0, as Z1>0 and Z2>0. λ4=Z5, λ4<0, if and only if Z5<0. After the simplifications, we get
(17)β*A*(γ1Z2δo-γ1β2κS*2)+{γ1βδπηβ***M*-[(γ1βδπηβ**+β*Q1)Z1(δ+γ1)+Q1Z2δoβ*]}>0,
where β*=γ1β((δ+γ3)+κqγ2), β**=Q1(δ+γ1)(δ+γ3)(1-R0)+κβSoQ1, and {γ1βδπηβ***M*>[(γ1βδπηβ**+β*Q1)Z1(δ+γ1)+Q1Z2δoβ*]}, say G1=γ1βδπηβ***M* and G2=[(γ1βδπηβ**+β*Q1)Z1(δ+γ1)+Q1Z2δoβ*]β***=(δ+γ3)[(δ+γ1)+γ1μ1]. So, λ4 has negative real part if (γ1Z2δo>γ1β2κS*2). For λ5=Z7, we obtained negative real parts, and by using the Z5 which is positive under the conditions described in λ4, and Z1>0, Z2>0, we get the negative real parts. Thus, all the eigenvalues have negative real parts, so by the Routh-Hurwitz criteria the endemic equilibrium point D* is locally asymptotically stable when R0>1. The proof is completed.


## 4.1. Stability of Endemic Equilibrium (EE)

In this subsection, we find the local asymptotic stability of EE aboutD*=(S*,E*,A*,C*,M*), and we prove the local stability of endemic equilibrium in the following.Theorem 2.
ForR0>1, the endemic equilibrium D* of system (3) is locally asymptotical stable, if the following conditions hold:
(13)γ1Z2δoS*2>γ1β2κ,G1>G2;
otherwise, the system is unstable.Proof.
Here, we prove the that the system (3) about the equilibrium point D* is locally stable, and for this, we obtain the Jacobian matrix J*(ζ) of the system (3), in the following:(14)J*(ζ)=(-(δ+δo+p+β(A*+κC*))δo-βS*-δo-δπη-βκS*-δo-δoβ(A*+κC*)-(δ+γ1)βS*βκS*+δπημ10γ1-Q10μ200qγ2-(δ+γ3)00000-(δ+μ1+μ2)).By elementary row operation and after simplification, we get the following Jacobian matrix:(15)J*(ζ)=(-Z1-δo-βS*-δoZ3-δo0-Z1(δ+γ1)-Z2δoZ1βS*-Z2(βS*+δo)Z1Z4+Z2Z3Z1μ1-Z2δo00Z5γ1(Z1Z4+Z2Z3)Z6000Z700000-(μ1+μ2+δ)),where
(16)Z1=δ+δo+p+β(A*+κC*),Z2=β(A*+κC*),Z3=-(δπη+βκS*+δo),Z4=(βκS*+δπη),Z5=-Q1(-Z1(δ+γ1)-Z2δo)+γ1(Z4βS*-Z2(βS*+δo)),Z6=-μ2(-Z1(δ+γ1)-Z2δo)+γ1(Z1μ1-Z2δo),Z7=-Z5(δ+γ3)-Z5γ1(Z1Z4+Z2Z3).
The eigenvalue λ1=-(μ1+μ2+δ)<0, λ2=-Z1, and using the value of Z1 and further C*, we get λ2=-(δ+δo+p+β((δ+γ3)+κqγ2)A*)<0. λ3=-(Z1(δ+γ1)+Z2δo)<0, as Z1>0 and Z2>0. λ4=Z5, λ4<0, if and only if Z5<0. After the simplifications, we get
(17)β*A*(γ1Z2δo-γ1β2κS*2)+{γ1βδπηβ***M*-[(γ1βδπηβ**+β*Q1)Z1(δ+γ1)+Q1Z2δoβ*]}>0,
where β*=γ1β((δ+γ3)+κqγ2), β**=Q1(δ+γ1)(δ+γ3)(1-R0)+κβSoQ1, and {γ1βδπηβ***M*>[(γ1βδπηβ**+β*Q1)Z1(δ+γ1)+Q1Z2δoβ*]}, say G1=γ1βδπηβ***M* and G2=[(γ1βδπηβ**+β*Q1)Z1(δ+γ1)+Q1Z2δoβ*]β***=(δ+γ3)[(δ+γ1)+γ1μ1]. So, λ4 has negative real part if (γ1Z2δo>γ1β2κS*2). For λ5=Z7, we obtained negative real parts, and by using the Z5 which is positive under the conditions described in λ4, and Z1>0, Z2>0, we get the negative real parts. Thus, all the eigenvalues have negative real parts, so by the Routh-Hurwitz criteria the endemic equilibrium point D* is locally asymptotically stable when R0>1. The proof is completed.


## 5. Global Stability of DFE

In this section, we present the global stability of disease-free equilibrium DFE of the system (3). For different biological model, the Lyapunov function was used by [20, 21] for the global stability. For our model, we define and construct Lyapunov function in the following for the global stability of DFE. Further the global stability of endemic equilibrium we use the Lyapunov function and find its global asymptotical stability.Theorem 3.
ForR0≤1, the disease-free equilibrium of the system (3) is stable globally asymptotically, if S=So and unstable for R0>1.Proof.
Here, we define the Lyapunov function for the global stability of disease-free equilibrium, given by(18)V(t)=[d1(S-So)+d2E+d3A+d4C+d5M].
Differentiating the previous function with respect to t and using the system (3),
(19)V′(t)=d1S′+d2E′+d3A′+d4C′+d5M′,V′(t)=d1[δπ(1-ηC(t))-δS(t)-β(A(t)+κC(t))S(t)+δo(1-(S(t)+E(t)+A(t)+C(t)+M(t)))-pS(t)]+d2[β(A(t)+κC(t))S(t)-δE(t)+δπηC(t)-γ1E(t)+μ1M(t)]+d3[γ1E(t)-δA(t)-qγ2A(t)-(1-q)γ1A(t)+μ2M(t)]+d4[qγ2A(t)-δC(t)-γ3C(t)]+d5[-(μ1+μ2)M(t)-δM(t)],
where di, i=1,2,…5, are some positive constants to be chosen later.
After the arrangement, we obtain(20)V′(t)=[d2-d1](β(A+κC))S+[d2-d1]δπηC+[d3γ1-d2(δ+γ1)]E+[d2μ1-d5δ]M+[d4qγ2-d3Q1]A+d1δπ-d1δS+d1δo-d1δo(S+E+A+C+M)-d1pS-d4(δ+γ3)C.
Choosing the constants, d1=d2=γ1,  d3=(δ+γ1),  d4=Q1(δ+γ1)/qγ2, and d5=γ1μ1/δ.
After the simplification, we get(21)V′(t)=-(S-So)-γ1δoE-γ1δoA-(γ1δo+(δ+γ3)(δ+γ1)Q1qγ2)C-(γ1δo+(δo+γ1)μ2)M,
where So=(δπ+δo)/(δ+δo+p). V′(t)=0 if and only if S=So and E=A=C=M=0. Also, V′(t) is negative for (S>So). So, by [22], the DFE is globally asymptotically stable in Γ. The proof is completed.


### 5.1. Global Stability of Endemic Equilibrium

In this subsection, we show the global asymptotical stability of the system (3). To do this, we state and prove the following theorem for the global stability of endemic equilibrium.Theorem 4.
ForR0>1, system (3) is globally asymptotically stable, if S=S* and δo>μ1, and unstable for R0≤1.Proof.
To prove that system (3) is globally asymptotically stable, we define the Lyapunov in the following:
(22)L(t)=(μ1+μ2+μ3)γ1γ1+δ(S-S*)+(μ1+μ2+μ3)γ1γ1+δE+(μ1+μ2+μ3)A+(μ1+μ2+μ3)Q1qγ2C+μ2M.
Taking the derivative with respect to time t, using the system (3),
(23)L′(t)=(μ1+μ2+μ3)γ1γ1+δ×[δπ(1-ηC(t))-δS(t)-β(A(t)+κC(t))S(t)+δo(1-(S(t)+E(t)+A(t)+C(t)+M(t)))-pS(t)]+(μ1+μ2+μ3)γ1γ1+δ[β(A(t)+κC(t))S(t)-δE(t)+δπηC(t)-γ1E(t)+μ1M(t)]+(μ1+μ2+μ3)[γ1E(t)-δA(t)-qγ2A(t)-(1-q)γ1A(t)+μ2M(t)]+(μ1+μ2+μ3)Q1qγ2[qγ2A(t)-δC(t)-γ3C(t)]+μ2[-(μ1+μ2+δ)M].
Simplifying, we obtained
(24)L′(t)=-(μ1+μ2+μ3)γ1(γ1+δ)(δ+δo+p)(S-S*)-(μ1+μ2+μ3)γ1δo(γ1+δ)[E+A+C]-(μ1+μ2+μ3)γ1(γ1+δ)[δo-μ1]M.
The endemic equilibrium of the system (3) is globally asymptotically stable for R0>1, if S=S* and δo>μ1. So, the endemic equilibrium of the system (3) is globally asymptotically stable. The proof is completed.


## 5.1. Global Stability of Endemic Equilibrium

In this subsection, we show the global asymptotical stability of the system (3). To do this, we state and prove the following theorem for the global stability of endemic equilibrium.Theorem 4.
ForR0>1, system (3) is globally asymptotically stable, if S=S* and δo>μ1, and unstable for R0≤1.Proof.
To prove that system (3) is globally asymptotically stable, we define the Lyapunov in the following:
(22)L(t)=(μ1+μ2+μ3)γ1γ1+δ(S-S*)+(μ1+μ2+μ3)γ1γ1+δE+(μ1+μ2+μ3)A+(μ1+μ2+μ3)Q1qγ2C+μ2M.
Taking the derivative with respect to time t, using the system (3),
(23)L′(t)=(μ1+μ2+μ3)γ1γ1+δ×[δπ(1-ηC(t))-δS(t)-β(A(t)+κC(t))S(t)+δo(1-(S(t)+E(t)+A(t)+C(t)+M(t)))-pS(t)]+(μ1+μ2+μ3)γ1γ1+δ[β(A(t)+κC(t))S(t)-δE(t)+δπηC(t)-γ1E(t)+μ1M(t)]+(μ1+μ2+μ3)[γ1E(t)-δA(t)-qγ2A(t)-(1-q)γ1A(t)+μ2M(t)]+(μ1+μ2+μ3)Q1qγ2[qγ2A(t)-δC(t)-γ3C(t)]+μ2[-(μ1+μ2+δ)M].
Simplifying, we obtained
(24)L′(t)=-(μ1+μ2+μ3)γ1(γ1+δ)(δ+δo+p)(S-S*)-(μ1+μ2+μ3)γ1δo(γ1+δ)[E+A+C]-(μ1+μ2+μ3)γ1(γ1+δ)[δo-μ1]M.
The endemic equilibrium of the system (3) is globally asymptotically stable for R0>1, if S=S* and δo>μ1. So, the endemic equilibrium of the system (3) is globally asymptotically stable. The proof is completed.


## 6. Numerical Simulations

In this section, we present the numerical simulation of the proposed model (3), by using the Runge-Kutta order four scheme. For different values of the parameters, the numerical results are presented in Figures 2, 3, 4, 5, 6, 7, and 8. The variation of migration parameters μ1 and μ2, with different values, is presented. In the numerical solution of the model, the parameters and their values are presented in Table 1. In our simulation, the susceptible individuals are shown by dashed line, the exposed individuals by dotted dashed, the acute individuals by bold line, the carrier by dashed line and the migrated individuals by red bold line. Figures 2 to 8 represent the compartmental population of hepatitis B individuals with migration effect. The values presented in Table 1 are fixed except for μ1 and μ2. In Figure 2, by the values for μ1=0.9 and μ2=0.9, we see that the population of exposed and acute individuals is increasing. In Figure 3, we set μ1=0.8 and μ2=0.9, and the population of exposed and acute individuals decreased. Decreasing the values of μ1 and μ2, we obtain different results; see Figures 2 to 8. When we decrease the population of immigrants who have the HBV virus, we see the decrease in the population of exposed and acute individuals. The parameters presented in Table 1 were used by different authors, for example, the natural death rate equally birth rate by [4], the rate at which the latent individuals become infectious by [3], the rate at which the individuals move to the carrier class by [3], the rate at which the individuals move to the carrier class [3], the transmission coefficient β by [24], δo the loss of immunity rate by [13], the value of η by [13], the value of κ by [10], and π and q from [23]. We assume the values for the parameters γ3, p, μ1, and μ2 in our simulation.Table 1
Parameter values used in numerical simulations.
Notation
Parameter description
Range
Source
δ
Natural death rate equally birth rate
0.0143
[4]
π
The failure immunization
0-1
[23]
γ
1
The rate at which the latent individuals become infectious
6
[3]
γ
2
The rate at which the individuals move to the carrier class
4
[3]
γ
3
The rate of flow from carrier to the vaccinated class
0.34
Assumed
β
The transmission coefficient
0.8
[24]
q
The proportion of individuals become carrier
0.005
[23]
δ
o
Represent the loss of immunity
0.06–0.03
[13]
p
Represent the vaccination of susceptible
0.3
Assumed
μ
1
The rate of flow from Migrated class to exposed class
0.23
Assumed
μ
2
The rate of flow from Migrated class to acute class
0.56
Assumed
η
Unimmunized children born to carrier mothers
0.7
[13]
κ
The infectiousness of carriers related to acute infection
0-1
[10]Figure 2
The plot shows the HBV transmission model of hepatitis B, withμ1=0.90 and μ2=0.90.Figure 3
The plot shows the HBV transmission model of hepatitis B, withμ1=0.80 and μ2=0.90.Figure 4
The plot shows the HBV transmission model of hepatitis B, withμ1=0.70 and μ2=0.80.Figure 5
The plot shows the HBV transmission model of hepatitis B, withμ1=0.50 and μ2=0.60.Figure 6
The plot shows the HBV transmission model of hepatitis B, withμ1=0.20 and μ2=0.30.Figure 7
The plot shows the HBV transmission model of hepatitis B, withμ1=0.20 and μ2=0.40.Figure 8
The plot shows the HBV transmission model of hepatitis B, withμ1=0.10 and μ2=0.20.


## 7. Conclusion

A compartmental model of HBV transmission virus has been presented. A mathematical model has been obtained by adding(1) the migrated class, (2) HBV transmission rate between migrated class and exposed class, (3) transmission rate between migrated class and acute class, and (4) death rate of individuals in the migrated class. By adding these new features, we have obtained a compartmental model of HBV with migration effect. First, we obtained the basic reproduction number for the proposed model. The disease-free equilibrium is locally as well as globally asymptotically stable for R0≤1. We obtained the local and global asymptotical stability for the endemic equilibrium. For the reproduction number R0>1, the endemic equilibrium is, locally as well as globally, asymptotically stable. Furthermore, we have solved the compartment model numerically, and the results are presented in Figures 2 to 8. By changing the values of μ1 and μ2, different results have been obtained. It is concluded that when the value of μ1 and μ2 decreases, the population of (exposed, acute, and carrier) individuals also decrease. The proportion of infected individuals decreases when the proportion of migrated individuals (who have the HBV virus) decreases. So, the number of infected individuals is directly proportional to the number of migrated individuals.Last yet not the least, the authors of this work have agreed to devise in a course of time a more advanced model on restraining HBV transmission through migration.



---

*Source: 150681-2013-06-24.xml*