In this paper, the numerical properties of a nonlinear age-structured hepatitis B virus model with saturated incidence and spatial diffusion are studied. Applying linearly implicit Euler method in time integration, a numerical scheme which can preserve the biological meanings is constructed. The convergence of the numerical solution in finite time is explored. In stability analysis, a threshold is proposed, which is called numerical basic reproduction number and denoted by R0h. It is proved that the numerical solution is locally asymptotically stable at the disease-free equilibrium when R0h<1. Moreover, it is proved the numerical basic reproduction number converges to the exact basic reproduction number of the model with first order accuracy. Furthermore, it is shown that a numerical space independent equilibrium exists and is asymptotically stable if R0h>1, which implies the threshold stability of the model can be preserved by numerical solution proposed. Eventually, our conclusions are tested through numerical experiments.