The fractional order derivative was introduced to the seepage flow research to establish the relaxation models of non-Newtonian viscoelastic fluids in dual porous media. The flow characteristics of non-Newtonian viscoelastic fluids through a dual porous medium were studied by using the Hankel transform, the discrete Laplace transform of sequential fractional derivatives and the generalized Mittag-Leffler function. Exact solutions were obtained for arbitrary fractional order derivative. The long-time and short-time asymptotic solutions for an infinite formation were also resulted. The pressure transient behavior of non-Newtonian viscoelastic fluids flow through an infinite dual porous media was studied by using Stehfest's inversion method of the numerical Laplace transform. It shows that the characteristics of the fluid flow are appreciably affected by the order of the fractional derivative.
