The spatial distribution pattern of pine wilt disease was analyzed with geostatistics. The results showed that pine wilt disease was an aggregated distribution within China, at the scale of the entire country, the semivariogram function of spatial distribution can be described by a spherical model with a spatial dependence range of 10.1888°. At the provincial scale, the distribution patterns of Jiangsu, Anhui, Guangdong and Zhejiang provinces could also be described by the spherical model with spatial dependence range of 3.4464°, 2.1581°, 3.0800°, and 2.4376° respectively. In addition,2 aggregated distribution areas were found. One was centered in Nanjing with a range of (30.5-32.5°N, 117.7-120.5°E), which include almost all of the area of Jiangsu and Anhui provinces as well as the north-west position of Zhejiang province. The second aggregated distribution area was in Guangdong Province (22.5 -24.0°N, 113.0-114.5°E). The natural distribution model of the early stage of pine wilt disease in China can be formulated as R=5.4743e0.4139t (R2=09204).