<p>Granular soil is a nonhomogeneous and discontinuous porous medium with a large amount of unique amorphous particles. The gradation and particle morphological characteristics are essential factors that affect the hydraulic and mechanical properties of granular soils in the numerical simulations of suffusion. The current study on the discrete element modeling of fractal graded soil with realistic particle shapes is not sufficiently detailed. In this paper, a new modeling method of granular soil is proposed, which enables the generation of granular soil specimens with predesigned fractal gradations and numerous irregular particle shapes. The proposed methodology is validated by generating five graded specimens with different particle shapes. A series of numerical simulations of fine migration are performed to evaluate the influence of the gradations and particle morphologies on the ability to prevent the loss of fine particles. The simulation results indicate that the capacity of the soil matrix to prevent the fines loss decreases with the increasing particle size ratio and increases with the decreasing fractal dimension. Furthermore, the angularity of skeletal particles plays a critical role in enhancing the suffusion resistance. This study contributes to advancing the understanding of the internal stability of graded soils, facilitating the filter layer design and accurate prediction of the erodible particles during suffusion.</p>