1. A method for building a velocity model for pre-stack depth migration for geological exploration comprising: receiving seismic, gravity and electromagnetic data relating to a geologic subsurface; and simultaneously joint inverting the seismic, gravity and electromagnetic data to produce the velocity model.