We have developed justified, verified and implemented for supercomputers with parallel architecture the reliable instrument for the studying complicated processes of waves' propagation in realistic 3D heterogeneous multiscale models of geological media - the special finite-difference method with local mesh refinement in time and space. On this base, it becomes possible to deal with realistic models and acquisitions and take into account both small-scale heterogeneities of a reservoir and mesoscale variations in overburden. This software was applied for simulation of seismic waves' propagation through realistic synthetic models developed for some East Siberian oil field. We have found that orientation of fracture corridors and fluid saturation of reservoir microstructure has very specific impact in synthetic images of scattered waves, which can be used to develop predictive criteria in real life data processing and interpretation. These criteria are verified by comparison of predictions with well log data (fracture orientation) and permeability (fluid saturation) of a collector by test results.