Surrogate models for seismological inverse problems: Constructing fast approximations to misfit functions