This paper presents a methodology that combines the hybrid snapshot simulation (Bai and Wang in Int J Comput Methods 18:2050029, 2020) and the discrete empirical interpolation method (DEIM) to reduce the computational cost of constructing a DEIM-based nonlinear reduced order model (ROM) while preserving its accuracy. In distinct contrast to its traditional counterpart, the time span of the hybrid snapshot simulation is divided into multiple intervals, and a majority of the intervals are simulated by local ROMs, while a fraction by the full order model (FOM). To tackle the challenge associated with limited FOM data for ROM-DEIM construction, a new approach is proposed to reconstruct and enrich the snapshot data of nonlinear model terms using solution data in ROM intervals. A formulation and procedure based on the cell-centered finite volume method (FVM) scheme are also developed to take into account two types of nonlinearities in ROM-DEIM: the componentwise and the transport-related, non-componentwise. A global ROM-DEIM is produced immediately at the end of the snapshot simulation and can be used to accelerate online simulation for different scenarios. The proposed methodology demonstrates excellent computational performance. Specifically, the computational time of the snapshot simulation is reduced by \(\sim\) 50% without compromising ROM-DEIM accuracy (relative error less than 0.8% compared with FOM). The results prove feasibility of combining hybrid snapshot simulation and DEIM for robust and efficient ROM-DEIM development.