The recent vast deployment of linear seismic arrays provides a good opportunity to conduct full waveform inversion (FWI) for high-resolution subsurface images. Full waveform inversion is considered to be the last course of data analysis in seismology and has been widely explored during the last decade. However, the application of FWI into real data is computationally expensive. To make the powerful and promising FWI computationally feasible, in this talk I will present the SEM-FK hybrid method for the numerical modelling of teleseismic incidence with local heterogeneities. The SEM-FK hybrid method will be the necessary tool for imaging the detailed subsurface structures beneath dense seismic arrays. Besides that, I will also introduce a new software package TomoQuake, which can be used to determine earthquake location, velocity, anisotropy, attenuation, and boundary topography by conducting traveltime inversions. We expect that TomoQuake can generate good starting models for seismic array imaging with the advanced full waveform inversion technology.