In order to simulate drifting and run-up ice floes driven by tsunamis, we developed a quasi-3D Discrete Element Method (DEM) in which phenomena unique to granular solids, such as the arch action, jam and pile up of such solids, can be simulated with a small computation load. The validity of the simulation method was verified by hydraulic experiments using synthetic ice. Some numerical simulations of drifting and run-up ice floes driven by a tsunami were performed to investigate the fundamental characteristics of ice run-up and the applicability to actual seas. It was confirmed that the simulation results could estimate risk areas where ice pile up and jam form as well as hazardous areas showing the degree of damage to buildings. Such calculation was expected to be useful for the compilation of hazard maps and the development of disaster mitigation plans in the future.