Note
Click here to download the full example code
Estimate homogeneous parametersΒΆ
Here we estimate transmissivity and storage from a pumping test campaign with the classical theis solution.
import welltestpy as wtp
campaign = wtp.load_campaign("Cmp_UFZ-campaign.cmp")
estimation = wtp.estimate.Theis("Estimate_theis", campaign, generate=True)
estimation.run()
Out:
Initializing the Shuffled Complex Evolution (SCE-UA) algorithm with 5000 repetitions
The objective function will be minimized
Starting burn-in sampling...
Initialize database...
['csv', 'hdf5', 'ram', 'sql', 'custom', 'noData']
* Database file '/home/docs/checkouts/readthedocs.org/user_builds/welltestpy/checkouts/v1.0.2/examples/Estimate_theis/2020-09-02_14-53-19_db.csv' created.
Burn-in sampling completed...
Starting Complex Evolution...
ComplexEvo loop #1 in progress...
ComplexEvo loop #2 in progress...
ComplexEvo loop #3 in progress...
ComplexEvo loop #4 in progress...
ComplexEvo loop #5 in progress...
ComplexEvo loop #6 in progress...
ComplexEvo loop #7 in progress...
ComplexEvo loop #8 in progress...
ComplexEvo loop #9 in progress...
ComplexEvo loop #10 in progress...
ComplexEvo loop #11 in progress...
ComplexEvo loop #12 in progress...
ComplexEvo loop #13 in progress...
ComplexEvo loop #14 in progress...
ComplexEvo loop #15 in progress...
ComplexEvo loop #16 in progress...
ComplexEvo loop #17 in progress...
ComplexEvo loop #18 in progress...
ComplexEvo loop #19 in progress...
THE POPULATION HAS CONVERGED TO A PRESPECIFIED SMALL PARAMETER SPACE
SEARCH WAS STOPPED AT TRIAL NUMBER: 2572
NUMBER OF DISCARDED TRIALS: 0
NORMALIZED GEOMETRIC RANGE = 0.000368
THE BEST POINT HAS IMPROVED IN LAST 100 LOOPS BY 100000.000000 PERCENT
*** Final SPOTPY summary ***
Total Duration: 1.97 seconds
Total Repetitions: 2572
Minimal objective value: 77.2517
Corresponding parameter setting:
mu: -9.21583
lnS: -9.10167
******************************
Best parameter set:
mu=-9.215829527790921, lnS=-9.101669397418332
In addition, we run a sensitivity analysis, to get an impression of the impact of each parameter
estimation.sensitivity()
Out:
Initializing the Fourier Amplitude Sensitivity Test (FAST) with 260 repetitions
Starting the FAST algotrithm with 260 repetitions...
Creating FAST Matrix
Initialize database...
['csv', 'hdf5', 'ram', 'sql', 'custom', 'noData']
* Database file '/home/docs/checkouts/readthedocs.org/user_builds/welltestpy/checkouts/v1.0.2/examples/Estimate_theis/2020-09-02_14-53-30_sensitivity_db.csv' created.
*** Final SPOTPY summary ***
Total Duration: 0.28 seconds
Total Repetitions: 260
Minimal objective value: 547.005
Corresponding parameter setting:
mu: -9.04744
lnS: -9.80937
Maximal objective value: 2.78196e+06
Corresponding parameter setting:
mu: -15.855
lnS: -11.2826
******************************
260
Parameter First Total
mu 0.655260 0.939710
lnS 0.058136 0.397956
260
Total running time of the script: ( 0 minutes 12.407 seconds)