In this talk we introduce an energy-preserinvg numerical scheme for a class of stochastic Lotka-Volterra systems, as well as a class of discrete-gradient-projection methods that simultaneously preserve all invariant Hamiltonians for the underlying stochastic Poisson systems. We prove they can also preserve the Casimir functions of the systems under certain conditions. Mean-square convergence orders of the methods are investigated.