This paper explores a high-dimensional stochastic SIS epidemic model characterized by a mean-reverting, stochastic process. Firstly, we establish the existence and uniqueness of a global solution to the stochastic system. Additionally, by constructing a series of appropriate Lyapunov functions, we confirm the presence of a stationary distribution of the solution under R0s>1. Taking 3D as an example, we analyze the local stability of the endemic equilibrium in the stochastic SIS epidemic model. We introduce a quasi-endemic equilibrium associated with the endemic equilibrium of the deterministic system. The exact probability density function around the quasi-stable equilibrium is determined by solving the corresponding Fokker–Planck equation. Finally, we conduct several numerical simulations and parameter analyses to demonstrate the theoretical findings and elucidate the impact of stochastic perturbations on disease transmission.