In this paper the problems of constructing sonar images of the seabed according to measurements of the multibeam side scan sonar are considered. The inverse problem for the non-stationary equation of radiation transfer with the diffuse reflection conditions at the boundary which consists in finding the discontinuity lines of the bottom scattering coefficient is investigated. A numerical algorithm for solving the inverse problem is developed, and an analysis of the quality of reconstructing the boundaries of inhomogeneities of the seabed is carried out, depending on the number of views and the width of a radiation pattern and the sounding range.