The current level of development of numerical methods and high-performance computer systems opens the way to obtain detailed information about the structure of geological objects using 3D seismic study. A universally recognized necessary component that ensures the successful development of modern high-tech technologies for acquiring, processing, and interpreting geophysical data is the complete digital models of geological objects - their digital counterparts. It is on this basis that a detailed assessment of the resolution and information content of the proposed methods and their comparison with the already known processing and interpretation algorithms using the example of a specific geological object becomes possible. Besides, the presence of such digital models allows you to determine the optimal acquisition system, focused on the study of specific features of the object being studied and the selection of the most appropriate graph for processing the data obtained. In this paper, the primary attention is paid to the construction of a realistic three-dimensional geological model with a family of faults, as well as fracture corridors and clusters of cavities. After constructing such an inhomogeneous multi-scale model, we perform a finite-difference numerical simulation of 3D seismic waves’ propagation. The data obtained are processed using the original procedures for extracting scattered/diffracted waves with the subsequent construction of images of the corresponding small-scale objects, which generate these waves. The results obtained are using for verification of the algorithms of scattering and diffraction imaging as well as full waveform inversion.