[ENHANCEMENT] Hantush Well Model rfunc about pastas HOT 1 CLOSED

Interesting idea @ursusmutus, I'm curious if you're able to get good results using that approach!

In the past I've split wells in different aquifers over different HantushWellModel response functions. So all wells in each aquifer have the same response scaled with distance, but different aquifers can respond differently. From a hydrological point of view, and the current capabilities in Pastas I would think that separate response functions for each aquifer would be the more appropriate choice, but I'm curious to see how your approach performs. Would you be willing to share a comparison between a model using your approach and an approach where each aquifer is modeled separately?

For grouping wells in multi-aquifer systems using a single response function, it would make more sense to use a response function that accounts for that explicitly, but that is not currently implemented in Pastas. We have recently been talking about implementing a Pastas-plugins structure that would allow more custom StressModels to be applied. An example would be an multi-layer analytic element model using TTim. The analytic model is then used to compute a response to pumping. This is a kind of hybrid form of time series modeling and more traditional groundwater modeling. This still requires quite some testing and development however, so this won't be available for some time.

I'm converting this issue into a discussion for now, as I think this requires a bit more thinking and testing before we have a concrete idea of what we want to implement in Pastas.

from pastas.