Abstract: This paper studies the main pollution sources to the groundwater of fracture-karst in Middle Odovician carbonaterockes, which is the unique sources for the urban water supply,and the main contaminants in groundwater, its characteristics of space-distribution, estimates groundwater quality by fuzzy mathematic method. Trend surface analysis and Kriging method in Geostatistical methods are taken to simulate contaminants distribution in space, grey system method and characteristic finite element solution of advantages and disadvantages of these methods are compared. Finally some proposal for protection of ground water environment in Boshan are suggested.