DEM(Discrete Element Method) is a newly developed and prospective numerical simulation method. However, it takes tremendous time to determine the neighborhood relation among elements. In this paper an improved algorithm (Region Method) is introduced to reduce the time. The analysis of the complication of the algoritm and the comparison with the previous method (Window Method) are also presented. It illustrates the obvious advantage of the new algorithm.