Boshan is an industrial city in central Shandong Province. Groundwater is the unique source for the urban water supply. The main water source is fracture-karst water in Middle Ordovician carbonate rocks. This paper discusses the main pollution sources and the path of pollution, the main contaminants in groundwater and the characteristics of their spatial distribution and evaluation of groundwater quality with the fuzzy mathematic method. Geostatistical methods, such as the trend analysis and kriging methods, were used to simulate the distribution of contaminants in groundwater. The grey system method and characteristic finite element solution of the advective-dispersive equation were adopted to forecast pollution in the future. Finally some proposals for protection of groundwater environment in the Boshan area are put forward.