Prediction of the non-linear flow in porous media is still a major scientific and engineering challenge, despite major technological advances in both theoretical and computational thermodynamics in the past two decades. Specifically, essential controls on non-linear flow in porous media are not yet definitive. The principal aim of this paper is to develop a meaningful and reasonable quantitative model that manifests the most important fundamental controls on low velocity non-linear flow. By coupling a new derivative with fractional order, referred to conformable derivative, Swartzendruber equation and modified Hertzian contact theory as well as fractal geometry theory, a flow velocity model for porous media is proposed to improve the modeling of Non-linear flow in porous media. Predictions using the proposed model agree well with available experimental data. Salient results presented here include (1) the flow velocity decreases as effective stress increases; (2) rock types of “softer” mechanical properties may exhibit lower flow velocity; (3) flow velocity increases with the rougher pore surfaces and rock elastic modulus. In general, the proposed model illustrates mechanisms that affect non-linear flow behavior in porous media.