We provide the detailed asymptotic behavior for first-order aggregation models of heterogeneous oscillators. Due to the dissimilarity of natural frequencies, one could expect that all relative distances converge to definite positive value and furthermore that each oscillator converges to a possibly different stationary point. In order to establish the desired results, we introduce a novel method, called dimension reduction method that can be applied to a specific situation when the degree of freedom of the natural frequency is one. In this way, we would say that although a small perturbation is allowed, convergence toward an equilibrium of the gradient flow is still guaranteed. Several first-order aggregation models are provided as concrete examples by using the dimension reduction method to study the structure of the equilibrium, and numerical simulations are conducted to support theoretical results.