Land seismic surveys are often carried out in regions with a complex structure of the upper part of the section. The main complicating factors are significant differences in relief heights, inhomogeneity of the low-velocity zone and the presence of high-velocity layers of permafrost. An important stage of processing, in this case, is the construction of a velocity model of the nearsurface area based on the times of the first arrivals of seismic waves and its further accounting when constructing seismic images. The presence of boundaries in the resulting model may be necessary when choosing a floating or final datum or embedding in a depth model. The article discusses approaches to automating the procedure for constructing a near-surface model within the framework of the reflected waves method. First, a method for automatic picking of the first break times was implemented based on convolutional neural networks. Testing on real data has shown that the use of neural networks provides a more robust acquisition of first break times compared to standard approaches implemented in processing packages. Secondly, a method is proposed for constructing a layered near-surface velocity model based on the first break times. We use ray seismic tomography to build a smooth velocity model, then we convert it to a layered model. Testing on synthetic data simulating the geological conditions of Western Siberia has shown the possibility of building the layered near-surface model with good accuracy. Finally, an example of the implementation of methods in the form of modules in a processing package for their further use for processing real data is shown.
