Obtaining time sections required for structural constructions is one of the main procedures at the stage of preliminary digital processing and interpretation of seismic exploration materials. The optimality of obtaining time sections will ultimately consist in finding the values of normal moveout adjustments that best approximate the true values of corrections for target reflected waves. A description is given of the methodology developed by the authors for predicting the lithological composition of clay-carbonate rocks based on the use of well logging (WL) data and velocity analysis of CDP seismic records. When choosing a particular velocity model, the question arises about the nature of the interfaces between the media: are seismic boundaries associated with velocity jumps or with a gradual transition of some rocks and their velocities to others? This issue remains poorly understood to date. An algorithm for approximating discontinuous velocity functions with continuous functions is presented, brought to numerical implementation on a computer. The section describes algorithms for solving forward and inverse problems of transforming depth and time kinematic sections for a number of velocity models.