Systems and methods are disclosed. The method includes obtaining a sparse
seismic velocity model for a subterranean
region of interest. The sparse
seismic velocity model includes populated positions and missing positions, where each populated position includes a velocity. The method further includes determining an interpolated
seismic velocity model by determining the velocity at each missing position using the velocity at one or more of the populated positions and interpolation. The method still further includes determining an intermediate dense seismic velocity model by updating the velocity at each missing position using, at least in part, the interpolated seismic velocity model and a
projections onto convex sets method and determining a dense seismic velocity model using, at least in part, the intermediate dense seismic velocity model.