We present a novel fast synthetic aperture radar (SAR) image simulation method based on the echo matrix cell algorithm including multiple scattering. To improve the efficiency of SAR image simulation while ensuring the fidelity of the simulated results, we first discretized the target facets set in the SAR beams footprint into lattice targets using the range-Doppler (RD) imaging geometry model and provided the basis for simulating electromagnetic wave transmission. Based on the simulation of electromagnetic waves transmission, we used the ray tracing algorithm to calculate the multiple backscattering field including multi-polarimetric information and various material properties. Then, based on the echo matrix cell algorithm, we set the echo matrix cell as the subfield of the target backscattering field and designed the CUDA kernel function to implement a computation parallelization for SAR echo generation. All the echo matrix cells are traversed in parallel to obtain the total backscattering field of the target, reproducing the time-varying characteristic of the backscatter coefficient for each lattice target within the synthetic aperture time. The echo signal is processed using the RD imaging algorithm to obtain the simulated SAR image. Finally, we select some targets including aircraft carrier and airplane models for simulation tests. The computation efficiency is improved over 170-fold by comparing the computations of the proposed method and CPU single-thread. We also performed some qualitative and quantitative evaluations on the fidelity of the simulated SAR images. The experimental results verify the effectiveness of the proposed method.