In this paper, we propose and investigate the global dynamics of a SARS-CoV-2 infection model with diffusion and antiviral treatment. The proposed model takes into account the two modes of transmission (virus-to-cell and cell-to-cell), the lytic and nonlytic immune responses. The diffusion into the model is formulated by the regional fractional Laplacian operator. Furthermore, the global asymptotic stability of equilibria is rigorously established by means of a new proposed method constructing Lyapunov functions for a class of partial differential equations (PDEs) with regional fractional Laplacian operator. The proposed method is applied to the classical reaction-diffusion equations with normal diffusion.