Abstract: Numerical solver using a uniform grid is popular due to its simplicity and low computational cost, but would be unfeasible in the presence of tiny structures in large-scale media. It is necessary to use a nonuniform grid, where upsampling the wavefield from the coarse grid to the fine grid is essential for reducing artifacts. In this paper, we suggest a local refinement scheme using the Fourier interpolation, which is superior to traditional interpolation methods since it is theoretically exact if the input wavefield is band limited. Traditional interpolation methods would fail at high upsampling ratios (say 50); in contrast, our scheme still works well in the same situations, and the upsampling ratio can be any positive integer. A high upsampling ratio allows us to greatly reduce the computational burden and memory demand in the presence of tiny structures and large-scale models, especially for 3D cases.