The traditional finite element method (FEM) usually exhibits significant numerical dispersion error for solving the Helmholtz equation in relatively high-frequency range, resulting in insufficiently accurate solutions. To address this problem, this paper proposes a novel enriched finite element method (EFEM) based on radial basis functions (RBFs) which are frequently used in meshless numerical techniques. In the proposed method, the partition of unity (PU) framework is retained, and nodal interpolation functions are formed using the RBFs. Furthermore, the linear dependence (LD) problem commonly encountered in many of the PU-based methods using polynomial basis functions (PBFs) is effectively avoided by using the present RBFs. To enrich the approximation space generated by the RBFs, the PBFs are introduced to construct the local enrichment functions. Several typical numerical experiments are conducted in this work. The results indicate that the proposed method can significantly reduce the dispersion error and yield accurate solutions even for relatively high-frequency Helmholtz problems. More importantly, the proposed method can be directly implemented with standard quadrilateral meshes as in FEM. Therefore, the proposed method represents a promising numerical scheme for solving relatively high-frequency Helmholtz problems.