The method proposed and implemented in (Pozdniakov et al. (1998)) for imaging scatterers/diffractors in 2D situations modified for 3D heterogeneous media and is applied for mapping of 3D cracked areas, decompression areas and so on within consolidated granitoid subsurface block. This block is chosen as one of the possible sites for subsurface nuclear waste deposit. The approach bases on the use of asymmetric samples from multi offset data in order to avoid regular reflections and to concentrate on scattered/diffracted waves. This is possible because of the almost circular dispersion index for scattering via small intrusions. In the paper results of real data processing are presented. Next, this approach is applied in order introduce new seismic attribute - intensity of scattered waves. It was computed for some real oil field (Nepsko-Botoubinskoe) in East Siberia region. Regression analysis is performed in order to justify reliability of this seismic attribute as predication factor for estimation of collector capacity.