Nonlinear seismic imaging is designed to specify the velocity-depth model (VDM) and anisotropy parameters on the residual moveout corrections after depth migration pre-stack and also to build a primary VDM in time and RMS velocity of the reflecting horizons. The simultaneous refinement of all or selected parameters given a priori (borehole) information for isotropic or poor anisotropic subsurfaces (VTI\TTI).