In conventional analysis of the regional groundwater flow, the model parameters T, S are identified by a direct comparison of the computed head values with the observed head data. However, we should keep in mind that the computed values are macroscopic solutions, while the observed data include local random-components caused by the heterogeneity of aquifer. This study aims to propose a new method for the parameter identification, which takes full advantage of the statistical property on uncertainty of actual head behavior.