We propose an expanded mixed finite element method (EMFEM) with sharp recognition for such intrinsic natures of a Markov process as the probability conservation, the Chapman-Kolmogorov equation, skewness and excess kurtosis for the fractional-order advection diffusion equation (FADE) oriented in financial engineering. We improve the result in [11] for purely fractional diffusion model by modifying the domain of one bilinear form so that it is well defined for the appearance of the advection term, and the result in [13] for second order elliptic models by using lower-order finite element spaces to reduce computation costs. The solvability and optimal convergence rates of the EMFEM are proved by showing the coerciveness and the inf-sup condition for the modified bilinear forms over the lower-order finite element spaces, and thus the systematically mathematical theory of the EMFEM for general FADEs are developed. Numerical experiments are conducted to confirm these theoretical findings.