The article presents a numerical model designed for determining groundwater dynamics and water balance of the catchments of the Oziąbel (Czarna Woda) river and the Wołczyński Strumień river in Wołczyn region. Hydrogeological mapping and modelling research covered the area of 238.9 km2. As a result of measurements performed in 2008-2009, flows were determined in major rivers and water table positions were measured at 26 points. In the major part of the area described, the water table, lying at the depth of 1.5-18.7 m, has unconfined character, and the aquifer is built of Neogene (Quaternary) sands and gravels. In the area under study, groundwaters are drawn from 6 wells with total withdrawal of 6133 m3/d. The numerical modelling was performed with the use of Visual Modflow 3.1.0 software. The area was partitioned by a discretization grid with a step size