A method, including: obtaining an initial model of a subsurface property; simulating synthetic data from the initial model; obtaining recorded borehole seismic data, wherein the recorded borehole seismic data was obtained with a seismic source or receiver located in a well; and inverting, with a computer, the recorded borehole seismic data by full wavefield inversion, wherein the full wavefield inversion includes comparing the synthetic data to the recorded borehole seismic data, and computing a cost function, obtaining a gradient function from the cost function, wherein the gradient function is related to a change in the objective function with an incremental change in model parameters, using the initial model to compute an illumination function or a resolution function for seismic sources and receivers, and obtaining a conditioned gradient function by conditioning the gradient function with the illumination function or the resolution function.