Mean-reverting SDEs of the form dXt=κ(θ−Xt)dt+σ(Xt)dWt pull Xt toward θ at rate κ; since E[∫σdW]=0, moments are found by taking expectations and solving the resulting ODE, not by simulating the noise.
Key facts:
Xt=θ+(X0−θ)e−κt+σ∫0te−κ(t−s)dWs(OU/Vasicek, constant σ),
E[Xt]=θ+(X0−θ)e−κt,Var(Xt)=2κσ2(1−e−2κt).
As t→∞: Xt→N(θ,σ2/2κ).
For CIR, dXt=κ(θ−Xt)dt+σXtdWt: mean is identical to OU/Vasicek,
Var(Xt)=X0κσ2(e−κt−e−2κt)+θ2κσ2(1−e−κt)2.
Feller condition 2κθ≥σ2 keeps Xt≥0.