Response function parameters

This notebook demonstrates the use of response function parameter calculation functions for various aquifer parameters.

from pastas_plugins.responses.parameters import (
    exponential_parameters,
    hantush_parameters,
    kraijenhoff_parameters,
    theis_parameters,
)

0. Aquifer parameters

S = 0.1  # Storativity
K = 10.0  # Hydraulic conductivity [L/T]
D = 50.0  # Aquifer thickness [L]
x = 0.0  # Location in aquifer (center at 0.0)
L = 100.0  # Aquifer length [L]

1. Kraijenhoff van de Leur Parameters

Homogeneous aquifer between two parallel canals

A, a, b = kraijenhoff_parameters(S, K, D, x, L)
print(f"Kraijenhoff Parameters: A={A:.4f}, a={a:.4f}, b={b:.4f}")

Kraijenhoff Parameters: A=2.5000, a=0.2026, b=0.0000

2. Exponential Parameters

Linear reservoir system

c = 200.0  # Drainage resistance [T]

A, a = exponential_parameters(S, c)
print(f"Exponential Parameters: A={A:.4f}, a={a:.4f}")

Exponential Parameters: A=200.0000, a=20.0000

3. Hantush Parameters

Well in a confined aquifer

c = 500.0  # Aquitard resistance [T]
r = 50.0  # Distance from well [L]

A, a, b = hantush_parameters(S, K, D, c, r)
print(f"Hantush Parameters: A={A:.6f}, a={a:.4f}, b={b:.4f}")

Hantush Parameters: A=-0.000773, a=50.0000, b=0.0025

4. Theis Parameters

Well between two parallel canals (y=0)

A, a, b = theis_parameters(S, K, D, x, L)
print(f"Theis Parameters: A={A:.6f}, a={a:.4f}, b={b:.4f}")

Theis Parameters: A=-0.000159, a=0.2026, b=0.0000