A workflow for recovering fracture network characteristics from seismic data is considered. First, the presented discrete fracture modeling technique properly describes fracture models on the seismic scale. The key procedure of the workflow is 3D diffraction imaging based on the spectral decomposition of different combinations of selective images. Selective images are obtained by the prestack asymmetric migration procedure, whereas spectral decomposition occurs in the Fourier domain with respect to the spatial dip and the azimuth angles. At the final stage, we performed a topological analysis based on the construction of a merge tree from the obtained diffraction images. The results of the topological algorithm are modeling parameters for the discrete fractures. To analyze the effectiveness of our workflow, a statistical comparison of the recovered parameters and true model parameters was conducted. We used the Kolmogorov-Smirnov test for the statistical analysis of the fracture lengths, whereas the behavior of the Morisita index indicates the statistical distribution of the modeled fracture corridors. Numerical examples with synthetic realistic models provide a detailed, reliable reconstruction of the statistical characteristics of the fracture corridors.