This study was conducted in Larix principis-rupprechtii forests, located in GuanDi Mountain, Shanxi Province, China. Two plots（32 m×32 m） were selected which represented disturbed stand(plot A) and undisturbed(plot B) respectively. Based on theory and methodology of spatial pattern analysis in geostatistics, samples were taken at three soil depths (0~10 cm, 10~20 cm and 20~30 cm.) using soil core (7.0 cm in diameter) . The spatial heterogeneity of soil characteristics (soil moisture , pH ,total N )was analyzed. The result showed that soil moisture and its spatial heterogeneity decreased in disturbed stand, remained mainly spatial variation at small scales. Harvesting increased the means of soil pH (pH =6.86,plot A).Compared with soil pH at same soil depth in different stand, the proportion of randomly variation increased in disturbed stand and spatial autocorrelation variation became weaker. After cutting disturbance, the soil total nitrogen decreased significantly. Semivariogram analysis demonstrated that harvesting modified spatial pattern of soil total nitrogen , the spatial heterogeneity decreased at large scales.The spatial heterogeneity intensity of soil moisture, pH and total nitrogen became lower along with increase of soil depths.