In this thesis, we investigate a model for quantum gravity on finite noncommutative spaces using the topological recursion method originated from random matrix theory. More precisely, we consider a particular type of finite noncommutative geometries, in the sense of Connes, called spectral triples of type (1,0), introduced by Barrett. A random spectral triple of type (1,0) has a fixed fermion space, and the moduli space of its Dirac operator D={H,⋅},H∈HN, encoding all the possible geometries over the fermion space, is the space of Hermitian matrices HN. A distribution of the form e−S(D)dD is considered over the moduli space of Dirac operators. We specify the form of the action functional S(D) such that the topological recursion for a repulsive particles system, introduced by Borot, Eynard and Orantin, holds for the large N topological expansion of the n-point correlators Wn(x1,⋯,xn) of our model. In addition, we get the large N topological expansion of the free energy F=logZN and the n-point correlators of the model in terms of the enumerative combinatorics of the stuffed maps, introduced by Borot, whose elementary 2-cells may have the topology of a disk or of a cylinder. One can compute all the stable coefficients Wng(x1,⋯,xn),2g−2+n3˘e0,n≥1,g≥0 of the large N topological expansion of the n-point correlators Wn(x1,⋯,xn) of the model using the topological recursion formula, provided the leading order terms W10(x) and W20(x1,x2) are known. We show that, for our model, the leading order term W10(x) satisfies a quadratic algebraic equation y2+Q(x)y−P(x)=0. The spectral curve Σ of the model is a genus zero complex algebraic curve, given by the pre-mentioned quadratic equation. We find explicit linear (resp. quadratic) expressions for the coefficients of the polynomial Q(x) (resp. P(x)) in terms of the moments of the jump discontinuity of W10(x). We plane to investigate the spectral curve (Σ,ω10,ω20) of the model in more detail