Seawater intrusion is one of the most serious environmental problems in coastal areas of many countries, and is a hot issue world widely. The coast of Laizhou Bay is one of the most typical areas suffered from seawater intrusion. Seawater intrusion resulted in groundwater salinization, decreased the residents’ living quality and damaged the eco-environment seriously. Therefore, forecasting the trend of seawater intrusion, assessing the risk of seawater intrusion and proposed improvement measurements are the urgency tasks for scientific water resources management and sustainable development. 
　　　On the basis of summary of predecessors’ research achievements, this paper utilized the numerical simulation technique and spatial analysis approach of Geographic Information System to simulate seawater intrusion and assess the risk of it via SEAWAT4.0 and ArcGIS9.3. The main achievements are as follows:
　　　 (1) On the basis of analyzing hydrogeological conditions, the zonal hydrogeological concept model was constructed. 
　　　 (2) A variable density numerical groundwater flow and solute transport model was constructed to simulate seawater intrusion and identified and verified by SEAWAT4.0.
　　　 (3) Based on the analysis of the precipitation, groundwater withdrawal and sea level rise, four scenarios were set up for simulating the evolution of seawater intrusion using the previously verified model. The results showed the range of seawater intrusion was gradually decreasing with the pass of time, whereas seawater intrusion became more severely in the vicinity of Sanshandao and Diaolongju. Moreover, precipitation is the most important factor for effecting seawater intrusion; groundwater withdrawal is the secondly most important factor, sea level rise has little effects on seawater intrusion.
　　　 (4) Cl- was selected as the assessment index and used to classify the risk level. And then the risk assessment was made under the different scenarios designed previously. The results showed the high, medium and low risk areas were successively zonal distribution along the coastal line, and presented the character of “The middle is small, and the two points are big”; moreover, the area of risk zones was decreasing as time passing.