We extend FD-CSI method (Abubakar et al., 2009) to reconstruct density and velocity.
FD-CSI needs to solve full forward problem only once, making it highly efficient.
Multiplicative constraint is incorporated to improve the inversion results.
FD-CSI method is more efficient than the conventional nonlinear CG methods.
Two-phase inversion strategy can alleviate cross-talk effects between parameters.