We describe a large-scale simulation of the seismic wave propagation from the 2011 Tohoku-Oki earthquake. We combine structure models for topography, sediments, crust and subducting plates to construct a 3D structure model for the northeastern Japan area. Then we simulate the near-field wave propagation by using a FDTD method which is accelerated by the use of multi-GPUs (Graphics Processing Unit). The simulated wave-field exhibits strongly complex pattern reflecting the source complexity and the heterogeneities around the source region.