Biological processes, such as gene networks can be described either as a discrete or continuous Markov processes or a combination of both based on a chemical kinetics approach. Solutions of such mathematical descriptions emerge as probability distributions. In this work we focus on the continuous Markov process regime and propose a reduction framework that deals with the existence of multiple time scales. Trajectories that sample the probability distribution of continuous Markov processes, governed by Fokker-Plank equations, can be obtained as solutions of the Chemical Langevin equation (CLE) or systems of CLEs. CLEs are stochastic differential equations (SDEs) and like ordinary or partially differential equations, their solution can be a stiff numerical problem whenever the underlying physical system exhibits multiple time scales. The problem of multiple time scales can be addressed through model reduction of the system of CLEs. Although extensive work has been done in this field concerning deterministic deferential equations, the extension to a system of SDEs is not straightforward. In this work we develop a three step systematic framework that appropriately reduces the initial system of CLEs to subsystems with similar time scales. Specifically, by applying an appropriate linear transformation to the original system, the initial state vector of the CLE system is decomposed to fast and slow varying subvariables. This allows for a non-stiff description by treating each set differently. A set of sufficient and necessary conditions arises, which has to be met to ensure the existence of the appropriate transformation. The second step is to treat each of the two subsets independently. We extend and apply the method of adiabatic elimination to the systems under consideration. Fast variables are assumed to relax to a pseudo-stationary density under the assumption that the slow variables remain constant. Slow variables are approximated through a Fokker-Plank equation which governs the probability density of only the slow variables. Ultimately the distribution of the slow variables can be sampled through the solutions of a system of CLEs that correspond only to the slow subspace. The final step is to compute the approximated solution of the initial CLEs system by simply multiplying the two independent probability densities. The rigorousness of the proposed framework will be examined through illustrative chemical kinetics examples that lie in the continuous Markov process regime.