The boson peak, which represents an excess of vibrational states compared to Debye’s prediction at low frequencies, has been studied extensively, and yet, its nature remains controversial. In this study, we focus on understanding the nature of the boson peak based on the spatial heterogeneity of modulus fluctuations using a simple model system of a highly jammed two-dimensional granular material. Despite the simplicity of our system, we find that the boson peak in our two-dimensional system shows a shape very similar to that of three-dimensional molecular glasses when approaching their boson peak frequencies. Our finding indicates a strong connection between the boson peak and the spatial heterogeneity of shear modulus fluctuations.