AnaFlow provides several analytical and semi-analytical solutions for the groundwater-flow equation.

You can install the latest version with the following command:

pip install anaflow

You can find the documentation under https://anaflow.readthedocs.io.

In the following the well known Theis function is called an plotted for three different time-steps.

import numpy as np
from matplotlib import pyplot as plt
from anaflow import theis


time = [10, 100, 1000]
rad = np.geomspace(0.1, 10)

head = theis(time=time, rad=rad, transmissivity=1e-4, storage=1e-4, rate=-1e-4)

for i, step in enumerate(time):
    plt.plot(rad, head[i], label="Theis(t={})".format(step))

plt.legend()
plt.show()

The following functions are provided directly

• thiem Thiem solution for steady state pumping
• theis Theis solution for transient pumping
• ext_thiem_2d extended Thiem solution in 2D from Zech 2013
• ext_theis_2d extended Theis solution in 2D from Mueller 2015
• ext_thiem_3d extended Thiem solution in 3D from Zech 2013
• ext_theis_3d extended Theis solution in 3D from Mueller 2015
• neuman2004 transient solution from Neuman 2004
• neuman2004_steady steady solution from Neuman 2004
• grf "General Radial Flow" Model from Barker 1988
• ext_grf the transient extended GRF model
• ext_grf_steady the steady extended GRF model
• ext_thiem_tpl extended Thiem solution for truncated power laws
• ext_theis_tpl extended Theis solution for truncated power laws
• ext_thiem_tpl_3d extended Thiem solution in 3D for truncated power laws
• ext_theis_tpl_3d extended Theis solution in 3D for truncated power laws
• ext_thiem_int extended Thiem solution for the Integral variogram model
• ext_theis_int extended Theis solution for the Integral variogram model
• ext_thiem_int_3d extended Thiem solution in 3D for the Integral variogram model
• ext_theis_int_3d extended Theis solution in 3D for the Integral variogram model

We provide routines to calculate the laplace-transformation as well as the inverse laplace-transformation of a given function

• get_lap Get the laplace transformation of a function
• get_lap_inv Get the inverse laplace transformation of a function

• NumPy >= 1.14.5
• SciPy >= 1.1.0
• pentapy >= 1.1.0

You can contact us via info@geostat-framework.org.

MIT © 2019 - 2025