To increase the calculation speed of the computational fluid dynamics (CFD) based simulation for the gas-solid flow in fluidized beds, a Markov chain model (MCM) is developed to simulate the particle movement in a two-dimensional (2D) circulating fluidized bed (CFB) riser. As a statistic model, the MCM takes the CFD-discrete element method (DEM) results as samples to calculate the transition probability matrix of particles, which can describe the macroscopic regularities of the particle movement and guide the numerical simulation of the particle motion based on the Monte Carlo method. The particle distribution snapshots, residence time distribution (RTD) and mixing obtained from both the CFD-DEM and the MCM are compared. Results show that the computational speed of the MCM is faster than that of the CFD-DEM by about 2 orders of magnitude, and the difference of the mean particle residence time calculated by the two models is less than 2%. Besides, the MCM can well describe the time-averaged particle mixing compared to the CFD-DEM.