Lab experiment and mathematical simulation Modular three dimensional finite difference groundwater (MODFLOW) were performed in a soil tank to simulate the hydrogeochemical interaction between lake and typical unconfined aquifer. Results show that the velocity decreases exponentially with the transect distance on seepage face. The maximal velocity occurs at the top point of seepage face. The obliquity of seepage face has a great influence on the maximum and distribution of seepage velocity. With the increase of the obliquity of seepage face, the maximal velocity decreases quickly and the velocity distribution becomes much more even. Most of groundwater flow and pollutant flux discharges through a narrow portion near the top of seepage face. The flow and mass concentrated in the narrow portion increase with the decrease of the obliquity of seepage face. These will benefit to design a reasonable and economical scenario to manage lakeshore and to control the pollution of lake water near lakeshore.