Pumping test analysis with welltestpy Description In this workflow, we analyse two pumping test campaigns on two field sites: "Horkheimer Insel" (Heilbronn, Germany) "Lauswiesen" (Tübingen, Germany) The aim is to estimate parameters of heterogeneity from transient puming test data and to analyse their sensitivities. Target parameters mean of log-transmissivity variance of log-transmissivity length scale of log-transmissivity storage Applied methods The applied methods utilizing effecitive head solutions for the groundwater flow equation under a pumping test condition are described in: Zech, A., Müller, S., Mai, J., Heße, F., and Attinger, S.: Extending Theis' Solution: Using Transient Pumping Tests to Estimate Parameters of Aquifer Heterogeneity, Water Resour. Res., 52, 6156–6170, https://doi.org/10.1002/2015WR018509, 2016. These methods were implemented in welltestpy to automatically interprete pumping test data. The underlying type-curves are implemented in AnaFlow. Data sources The data for the "Horkheimer Insel" field site was manually taken from: Schad H.: Variability of hydraulic parameters in non-uniform porous media: experiments and stochastic modeling at different scales. University Tübingen; 1997. Ph.D. thesis. The pumping data from the "Lauswiesen" field site was kindly provided by Dr. Carsten Leven-Pfister and is made available on a repository of the University of Tübingen: Research Data Portal FDAT Citable as: Leven, C. (2020): Pumping Tests in Wells B1-B5 at the Hydrogeological Research Site Lauswiesen, Eberhard Karls Universität Tübingen, http://hdl.handle.net/10900.1/bfb0b0f7-7065-4a24-8a91-92ad8aa8fc40 Structure The workflow is organized by the following structure: data/ contains the campaign files for both sites in the welltestpy format contains time series for diagnostic plots src/ - contains the scripts to produce the results 00_wtp_plot.py - plotting well-constellation and campaign overviews 01_est_run.sh - bash file running 02_para_estimation.py in parallel 01b_est_run.sh - bash file running 02b_para_estimation.py in parallel 02_para_estimation.py - estimate parameters of heterogeneity from the pumping tests 02b_para_estimation.py - estimate equivalent parameters of homogeneity the pumping tests 03_postpro_results.py - plotting the estimation results for both sites 04_postpro_sensitivity.py - plotting the sensitivity results for both sites 05_est_radial_sens.sh - bash file running 06_rad_sens.py in parallel 06_rad_sens.py - estimate parameter sensitivites depending on the radial distance to the pumping well. when run in serial, results will be plotted. 07_comparison_len_scale.py - generate comparison plot for different length scales 08_check_unconfined_effect.py - generate diagnostic plots results/ - all produced results Python environment Main Python dependencies are stored in requirements.txt: welltestpy==1.0.3 anaflow==1.0.1 spotpy==1.5.9 mpi4py==3.0.2 matplotlib You can install them with pip (potentially in a virtual environment): pip install -r requirements.txt Contact You can contact us via info@geostat-framework.org. License MIT © 2021