Computer-implemented method for constructing a geological model of a subterranean formation formed from at least two sedimentary layers, on the basis of N' 0 measurements of the geometry of the limits of said layers, said N' 0 measurements of a said limit being carried out in a number N' 0 of cells of a two-dimensional grid G 0 , said grid G 0 being made up of N 0 cells, characterized in that : - for each limit of said layers, said geometry of said layer limit is determined in each cell of said grid G 0 by means of at least the following steps: a) nk grids G k are constructed from said grid G 0 , each of the grids G k , with k varying from 1 to nk, having the same spatial limits as the grid G 0 and having a number of cells lower than that of a grid G k-1 ; b) each of the grids G k , such that k is comprised between 1 and nk, is allocated via interpolation of said N' 0 measurements of said geometry in the corresponding cells of said grid G k ; c) said geometry of said layer limit is determined in each cell of said grid G 0 by determining said geometry of said layer limit in each cell of each of the nk grids G k , with k varying from 1 to nk, by applying, successively, from the grid G nk to said grid G 0 , an iterative method to minimize an objective function dependent on a property relative to said geometry of said layer limit, and by initializing each of the grids G k , with k varying from 0 to nk-1, with said geometry of said layer limit determined on the grid G k+1 ; and - on the basis of said geometry, of each of said layer limits, determined in each cell of said grid G 0 , said geological model is constructed. - for each limit of said layers, said geometry of said layer limit is determined in each cell of said grid G 0 by means of at least the following steps: a) nk grids G k are constructed from said grid G 0 , each of the grids G k , with k varying from 1 to nk, having the same spatial limits as the grid G 0 and having a number of cells lower than that of a grid G k-1 ; b) each of the grids G k , such that k is comprised between 1 and nk, is allocated via interpolation of said N' 0 measurements of said geometry in the corresponding cells of said grid G k ; c) said geometry of said layer limit is determined in each cell of said grid G 0 by determining said geometry of said layer limit in each cell of each of the nk grids G k , with k varying from 1 to nk, by applying, successively, from the grid G nk to said grid G 0 , an iterative method to minimize an objective function dependent on a property relative to said geometry of said layer limit, and by initializing each of the grids G k , with k varying from 0 to nk-1, with said geometry of said layer limit determined on the grid G k+1 ; and a) nk grids G k are constructed from said grid G 0 , each of the grids G k , with k varying from 1 to nk, having the same spatial limits as the grid G 0 and having a number of cells lower than that of a grid G k-1 ; b) each of the grids G k , such that k is comprised between 1 and nk, is allocated via interpolation of said N' 0 measurements of said geometry in the corresponding cells of said grid G k ; c) said geometry of said layer limit is determined in each cell of said grid G 0 by determining said geometry of said layer limit in each cell of each of the nk grids G k , with k varying from 1 to nk, by applying, successively, from the grid G nk to said grid G 0 , an iterative method to minimize an objective function dependent on a property relative to said geometry of said layer limit, and by initializing each of the grids G k , with k varying from 0 to nk-1, with said geometry of said layer limit determined on the grid G k+1 ; and - on the basis of said geometry, of each of said layer limits, determined in each cell of said grid G 0 , said geological model is constructed. <EOS>