In this paper, a stochastic SIRD model of Ebola with double saturated incidence rates and vaccination is considered. Firstly, the existence and uniqueness of a global positive solution are obtained. Secondly, by constructing suitable Lyapunov functions and using Khasminskii's theory, we show that the stochastic model has a unique stationary distribution. Moreover, the extinction of the disease is also analyzed. Finally, numerical simulations are carried out to portray the analytical results. MSC: Primary 37H10; 60H10; secondary 92C60; 92D30