Three bark beetles, Tomicus minor, T. yunnanensis and T. brevipilosus, are distributed in southwest China and mainly infest Pinus yunnanensis, causing serious damage to forestry production. In order to study the effect of the relationship of coexistence among three sympatric Tomicus species on their spatial distribution, the spatial distribution patterns of different degrees of infestation in the trunk were analyzed with geostatistical and traditional aggregation index methods. The results showed that the population density of T. yunnanensis was significant higher than that of T. minor in the heavily infested stand, contrary to the lightly infested stand. The results from the traditional index method showed that the spatial distribution patterns of three sympatric Tomicus species were all aggregated in the different degrees of infestation stands. The aggregation distribution of T. minor and T. yunnanensis resulted from their aggregation behavior and environment. The results from the geostatistical method showed that all presented the aggregation distribution, in addition to T. brevipilosus in the heavily infested stand. The level of spatial dependence of all populations ranged from 4.01 to 7.45 m, except the T. minor population in the heavily infested stand. The semivariogram models of T. minor and T. yunnanensis switched between a spherical model and a gaussian model. In short, the coexistence of three sympatric Tomicus species in different population densities did not affect their spatial distribution patterns, but impacted the semivariogram models and model parameters.