The objectives of this study were to reproduce the dynamic change process of groundwater in the period from 1970s to 1990s, during which local groundwater was in relatively low level, and to estimate the groundwater recharge process for the specified region by method of numerical simulation. The region which is located at the area in front of an estuary of a watershed (area: 10×12 km2 ) was adopted as the study area. Based on the 2 D numerical method in addition to the physical conditions of the soil vertical profile, a quasi-3 D numerical model was developed by considering water movement in the vertical direction using ground - water flow equation. The results of numerical simulation indicate that model efficiency is more than 80 % ; the groundwater recharge process for each aquifer was estimated based on numerical results of 1966-2005. The result of this study is looking forward to provide the reference method for the studies on numerical simulation and assessment of groundwater resources for regional groundwater.