This article deals with the discretization of hybrid dimensional model of Darcy flow in fractured porous media. These models couple the flow in the fractures represented as the surfaces of codimension one with the flow in the surrounding matrix. The convergence analysis is carried out in the framework of Gradient schemes which accounts for a large family of conforming and nonconforming dis-cretizations. The Vertex Approximate Gradient (VAG) scheme and the Hybrid Finite Volume (HFV) scheme are applied to such models and are shown to verify the Gradient scheme framework. Our theoretical results are confirmed by a few numerical experiments performed both on tetrahedral and hexahedral meshes in heterogeneous isotropic and anisotropic media.