Numerical simulation in the frequency-space domain has inherent advantages, such as： it is possible to simulate wave propagation from multiple sources simultaneously; there are no cumulative errors; only the interesting frequencies can be selected; and it is more suitable for wave propagation in viscoelastic media. The only obstacle to using the method is the requirement of huge computer storage. We extend the compressed format for storing the coefficient matrix. It can reduce the required computer storage dramatically. We get the optimal coefficients by least-squares method to suppress the numerical dispersion and adopt the perfectly matched layer （PML） boundary conditions to eliminate the artificial boundary reflections. Using larger grid intervals decreases computer storage requirements and provides high computational efficiency. Numerical experiments demonstrate that these means are economic and effective, providing a good basis for elastic wave imaging and inversion.
