A non-transitory computer-readable medium embodying a asset of instructions, which, when executed by one or more processors, cause the one or more processors to perform a computer-implemented method for an automated identification of surfaces in connection with building a geological-hydrodynamic model of an oil-and-gas deposit based on seismic data, the computer-implemented method comprising:
- 1) determining a model grid with a specified pitch and maximal allowable size on coordinate grid;
- 2) specifying control parameters of surfaces relative to the grid pitch;
- 3) determining a set of grid surfaces with respect to its pitch;
- 4) determining a domain in which the surface is defined and adding grid nodes where the surface is not defined;
- 5) specifying the parameters of the volumes of seismic data p, q;
- specifying the coordinates X, Y of the newly added node;
- specifying the coordinates of the adjacent node x, y, in which the surface is defined;
- and specifying the surface depth parameter in this node z(x,y);
- 6) arranging the depths in the given grid node in increasing order and setting the depth of the i-th surface in this node equal to the i-th depth;
- 7) determining the surfaces for processing;
- 8) determining the functionals and gradients of surfaces for processing;
- 9) determining the size of the current step along the gradient;
- 10) for each surface to be processed, determining what number of steps is required for the minimal value of the functional to be attained;
- in that case,
- 11) when the minimal value of the functional for each processed surface is less than the step along the gradient, making the displacement along the gradient by the found number of steps;
- 12) saving the size of this displacement in the memory, and determining the maximum of the displacements made;
- 13) temporarily excluding the surfaces for which a decrease in the functional was obtained from processing;
- 14) determining the surfaces for which the functional was not decreased, and once this is done, reiterating steps 7) to 13);
- 15) rearranging the depths of the surfaces determined in the given grid node in increasing order and setting the depth of the i-th surface in this node equal to the i-th depth;
- 16) determining the values of the current grid pitches from the new data on surfaces;
- and
- 17) reiterating all steps of the method until a model with parameters specified at steps 1 and 2 is obtained.