In the petroleum engineering and chemical industries, fluids engaging in displacement often have non-Newtonian properties, even though many former studies assume constant viscosities in the defending fluid. In this study, the computational fluid dynamics approach was performed in a two-dimensional model with uniformly distributed disks. This arrangement helps reveal the phenomenon and mechanics of how non-Newtonian characteristics of defending fluid affect two-phase displacement in porous media. Both global (in the whole medium) and regional (in the pore throat) studies revealed that shear-thinning makes capillary force and the pressure in the invading fluid decisive and leads to a uniform pattern. Meanwhile, the shear-thickening causes fingering due to the pressure drop in the defending fluid that becomes decisive. Cases of increasing injection rates were investigated to verify their ability to improve efficiency. The results verified that increased injection rates are effective in shear-thinning cases but energy-intensive when it comes to costs in shear-thickening cases. Finally, the viscosity ratio and capillary number (M-Ca) diagram were extended by plotting non-Newtonian cases as lines to consider viscosity variation. An estimation method was presented, which calculates the characteristic viscosity and locates non-Newtonian cases on an M-Ca diagram. This work can serve as a reference for enhanced oil recovery method development and microfluidic manipulation.