Two-dimensional depth-averaged advection-dispersion equation was simulated using FEM. In the straight rectangular channel, the advection-dispersion processes are simulated so that these results can be compared with analyti-cal solutions for the transverse line injection and the point injection. In the straight domain the standard Galerkin method with the linear basis function is found to be inadequate to the advection-dispersion analysis compared to the upwind finite element scheme. The experimental data in the S-curved channel were compared with the result by the numerical model using SUPG(Streamline upwind Petrov-Galerkin) method.