The paper studies a model of fluid flow in a fractured porous medium in which fractures are distributed uniformly over the volume. This model includes a nonlinear equation containing several terms with fractional derivatives in the sense of Caputo of order belonging to the interval 1,2. The relevance of studying this problem is determined by its practical significance in the oil industry, since most of the world’s oil reserves are in these types of reservoirs. The uniqueness of the solution to the problem in a differential form and its dependence on the initial data and the right-hand side of the equation is proved. A numerical method is proposed based on the use of the finite difference approximation for integer and fractional time derivatives and the finite element method in the spatial direction. A change of variables is introduced to reduce the order of the fractional derivatives. Furthermore, the fractional derivative is approximated by using the L1-method. The stability and convergence of the proposed numerical method are rigorously proved. The theoretical order of convergence is confirmed by the results of numerical tests for a problem of fluid flow in fractured porous media with a known exact solution.