Fangzheng county is the main production areas of rice in Heilongjiang Province,which also serves as one of the major producing counties of the enrichedSelenium(Se)rice.In this study,the spatial heterogeneity and the distribution of the soil total selenium and the relationship between selenium concentration variability and geochemical attributes were studied using a combination method of geostatistics and ArcGIS in Fangzheng county.The results showed that the variation function of soil total selenium was fit to a exponential model and the spatial heterogeneity index(Q)was 0.125,indicating that the variation mainly contributed to the structural factors and there was a strong spatial autocorrelation. The Ordinary Kriging(OK)was used to map the soil total selenium distribution of the whole county.From the map,it showed that the vast majority of the soils(69.3%)were sufficient with selenium,and the total soil selenium content was relatively high in the central and the southwest of the study area while lowvalue areas were mainly located in the northeastern.Among different soil types,peat soil had a relatively high Se content(0.267 mg·kg-1)and bog soil had a relatively low Se content(0.153 mg·kg-1).Considering the land use,paddy soils had a higher Se content(0.230 mg·kg-1)than upland soils(0.217 mg·kg-1).Correlation analysis results showed that soil total organic carbon，clay and silt content were the main factors affecting soil selenium content.