Dynamics of a stochastic HBV infection model with drug therapy and immune response

Hepatitis B is a disease that damages the liver, and its control has become a public health problem that needs to be solved urgently. In this paper, we investigate analytically and numerically the dynamics of a new stochastic HBV infection model with antiviral therapies and immune response represented by CTL cells. Through using the theory of stochastic differential equations, constructing appropriate Lyapunov functions and applying Itô's formula, we prove that the disease-free equilibrium of the stochastic HBV model is stochastically asymptotically stable in the large, which reveals that the HBV infection will be eradicated with probability one. Moreover, the asymptotic behavior of globally positive solution of the stochastic model near the endemic equilibrium of the corresponding deterministic HBV model is studied. By using the Milstein's method, we provide the numerical simulations to support the analysis results, which shows that sufficiently small noise will not change the dynamic behavior, while large noise can induce the disappearance of the infection. In addition, the effect of inhibiting virus production is more significant than that of blocking new infection to some extent, and the combination of two treatment methods may be the better way to reduce HBV infection and the concentration of free virus.