An approach to reconstructing a velocity depth model of the near surface using full wave inversion is proposed, implemented, and tested. We have shown previously that, although the involvement of multiples formed on the free surface increases the method stability, it decreases its resolving capacity. Therefore, at the first stage we engaged the entire wave field in processing, but afterwards we suppressed the free-surface related multiples. The results obtained demonstrate the potential of using this method to reconstruct quite complex near-surface structures, even if they have trap intrusions.
