Abstract In this paper, we propose a complex stochastic hepatitis B virus (HBV) epidemic model with vaccination strategy, where random fluctuations of the transmission dynamics of HBV is driven by Black-Karasinski process, for the first time. It is shown that Black-Karasinski process is a both biologically and mathematically reasonable assumption compared with existing stochastic modeling approaches. For the deterministic model, the basic reproduction number R₀, possible equilibria, and related asymptotic stability are studied. Then for the stochastic model, the existence and global positivity of the solution are proved. We further derive two stochastic critical values R₀^S and R₀^E related to R₀ to characterize the long-term behavior of HBV, and it turns out that (i) the stochastic model has a stationary distribution if R₀^S1; (ii) the infected individuals will go extinct exponentially fast when R₀^E1; (iii) R₀^S=R₀^E=R₀ if there is no environmental noise. Our results reveal that random fluctuations introduced will facilitate HBV prevalence. Moreover, by analyzing the stable structure of the model, we provide a complete classification and explicit approximation for the local density function of the stationary distribution. Finally, some numerical examples are performed to support our theoretical findings. The techniques and methods of analysis in this paper can be applied to many complex high-dimensional epidemic models motivated by Black-Karasinski process.
Han et al. (Wed,) studied this question.