Estimate steady homogeneous parameters

Here we estimate transmissivity from the quasi steady state of a pumping test campaign with the classical thiem solution.

import welltestpy as wtp

campaign = wtp.load_campaign("Cmp_UFZ-campaign.cmp")
estimation = wtp.estimate.Thiem("Estimate_thiem", 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.1/examples/Estimate_thiem/2020-09-02_14-47-54_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...
THE POPULATION HAS CONVERGED TO A PRESPECIFIED SMALL PARAMETER SPACE
SEARCH WAS STOPPED AT TRIAL NUMBER: 1366
NUMBER OF DISCARDED TRIALS: 0
NORMALIZED GEOMETRIC RANGE = 0.000991
THE BEST POINT HAS IMPROVED IN LAST 100 LOOPS BY 100000.000000 PERCENT

*** Final SPOTPY summary ***
Total Duration: 0.62 seconds
Total Repetitions: 1366
Minimal objective value: 0.0672645
Corresponding parameter setting:
mu: -9.21029
******************************

Best parameter set:
mu=-9.210293556543638

since we only have one parameter, we need a dummy parameter to estimate sensitivity

estimation.gen_setup(dummy=True)
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.1/examples/Estimate_thiem/2020-09-02_14-47-56_sensitivity_db.csv' created.

*** Final SPOTPY summary ***
Total Duration: 0.16 seconds
Total Repetitions: 260
Minimal objective value: 74.8884
Corresponding parameter setting:
mu: -9.1824
dummy: 0.179636
Maximal objective value: 2.35929e+06
Corresponding parameter setting:
mu: -15.976
dummy: 0.690593
******************************

260
Parameter First Total
mu 0.788866 0.978326
dummy 0.002011 0.054925
260