A divergence-free moving particle semi-implicit method is introduced for free-surface flow through porous media. Numerical incompressibility is conserved by solving additional pressure Poisson equation (PPE). Depending on current particle coordinates, a porosity-based factor is introduced to incorporate the effect of solid volume inside the porous domain. A hybrid formulation containing specified boundary condition and PPE is utilized on free-surface particles. The current framework is tested for four different problems. The first problem shows the effect of the proposed factor in vertical flow through a rectangular porous block and its representative volume change for different phases. Second and third problems validate the numerical model for dam break through a rectangular block of homogeneous porous media. In the fourth problem, flow through a trapezoidal porous block consisting of different porous media with variable effective porosity and permeability is simulated.