Abstract: Based on data collected in the field survey,spatial variation of benzoapyrene in soil was studied by means of geostatistical analyst with ArcGIS,and relationships of its contamination level with the natural environment and industry distribution of the region were analyzed.Results indicate that spatial heterogeneity of benzoapyrene concentration exists in the soil of the study area,and is especially distinct at the spatial scale of 240 km,with nugget/sill ratio C₀/（C₀+C） being 44%,suggesting that the spatial variation is made up of structural factors（e.g.soil organic matter） and random factors（e.g.distribution of pollution sources）.The steel industry,including coking plant,is the main emission sources of（benzoa）pyrene,and benzoapyrene contaminated soils are distributed around the steel industry developed cities and their vicinities.