Method, device and software package for locally determining the shape of geological horizons

Process to determine locally the shape of geological horizons, in which there is constructed a continuous function Sij,k(t) by interpolation or approximation of the discrete seismic traces of a three-dimensional seismic matrix, said Sij,k(t)'s function being designated as a “continous local seismic trace”, comprising the following steps:

• a). using as optimum offset of two adjacent continuous local seismic traces, the value of offset rendering maximal their correlation function, with optimal offset not being necessarily a whole multiple of the vertical sampling interval;
• b). taking as conditional neighborhood of a “central” continuous local seismic trace Sij,k(t) the sub-neighborhood consisting in adjacent traces Spq,k(t) corresponding to optimum offsets hij,pq,k associated with correlations Rij,pq,k(hij,pq,k) greater than a predetermined threshold comprised between 0 and 1;
• c). defining for each continuous local seismic trace Spq,k(t) of the conditional neighborhood, a value of residual relative to said “central” continuous local seismic trace Sij,k(t) comprising parametric coefficients;
• d). determining the parametric coefficients aij,k and bij,k by minimization of a set of residuals on the conditional neighborhood.