Method for rapidly computing updates to frequency-domain seismic wave fields by utilizing a matrix perturbation approach. The method speeds up model (e.g., velocity) parameter estimation by iterative inversion of measured seismic data. The method applies to the line search where the optimal size of the model update is estimated by testing different size updates to see which one generates the minimum objective function. By treating the model update as a perturbation, perturbation theory is used to relate the model perturbation to a corresponding wavefield perturbation. Thus, the Helmholtz equation is solved only once per iteration cycle.