SUMMARYIn acoustic wavefield modelling, the 2.5-D approximation, in which the model parameter is considered to be invariant in the strike direction and the source is assumed to be a point excitation so that the wavefield has 3-D features, is an economical and realistic scheme for seismic full-waveform inversion. On the other hand, according to Parseval's theorem, the least-squares inversion of time-domain seismic data is equivalent to that carried out with frequency-domain data, but the latter has some advantages such as its efficiency and flexibility in non-linear waveform (or diffraction) inversion with singleor multiple-frequency crosshole data. These inverse techniques require the computation of the response of a specified source and the Fréchet derivatives in the frequency domain, both of which can be calculated via the Green's function.In this paper, a damping method was applied to the computation of the 2.5-D Green's function for arbitrary acoustic media in the wavenumber-frequency domain. The 2.5-D acoustic and time-damped wave equation was solved with the finite-element method (FEM). The computations for homogeneous and heterogeneous media were performed in the frequency range 100-300 Hz. It was found that the time damper is important to the solution. A linear damper, as in time-domain modelling, can be chosen, but it depends on the frequency, the wavenumber and the width of the absorbing zone in the wavenumber-frequency-domain modelling. An appropriate time damper and a reasonable number of k y samples can be determined in terms of the range of frequency and the critical values of the wavenumber of the media. By comparing the numerical solutions with the exact analytic solution for the homogeneous case and the semi-analytic solution for two semi-infinite media in contact, it is shown that the method can yield satisfactory solutions in the frequency domain and the time domain. Finally, a synthetic experiment of crosshole seismic full-waveform inversion is used to illustrate the method.