I have a dataset containing high frequency measurements of the hydraulic head at different locations in a 5x6km area. Also known are the precipitation, evaporation, and pumped discharges from a number of wells (unfortunately only monthly totals for the latter). Using the Menyanthes software, timeseries analysis is performed returning the stationary drawdown at each location as a result of the groundwater withdrawal.
A very schematic crosssection of the pumping wells and measurement wells is given below. The red lines indicate pumping wells, blue lines are some of measurement locations, in reality there are 10 pumping wells and some 40 measurement wells. The thick dotted lines are faults, hatched lines indicate an aquitard, while crosshatching indicates an aquiclude. I'm mainly interested in the top layer as that is where the drinking water is extracted from, but I included the deeper layers for completeness.
From the calculated drawdown at each location I want to calculate (or at least estimate) the transmissivity of each aquifer so I can hopefully tell a bit more about the resistance caused by the faults (by smeared out clay layers mostly). I know of the Thiem and Theis models (see for example this article) for the drawdown in equilibrium and fluctuating situations, but are there any other models around that might be helpful to determine the geohydrological properties of the different compartments?
I realise this question probably doesn't have a single conclusive answer, but I think it is relevant for this site anyway as it tries to compare the available models/equations for evaluation of a definitely on-topic problem.