Under consideration are the numerical methods for simulation of a fluid flow in fractured porous media. The fractures are taken into account explicitly by using a discrete fracture model. The formulated single-phase filtering problem is approximated by an implicit finite element method on unstructured grids that resolve fractures at the grid level. The systems of linear algebraic equations (SLAE) are solved by the iterative methods of domain decomposition in the Krylov subspaces using the KRYLOVlibrary of parallel algorithms. The results of solving some model problem are presented. A study is conducted of the efficiency of the computational implementation for various values of contrast coefficients which significantly affect the condition number and the number of iterations required for convergence of the method.