Numerical preservation of long-term dynamics by stochastic two-step methods