![](\"./figs/Dawsonville.png\")
"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"### Load data"
]
},
{
"cell_type": "code",
"execution_count": 7,
"metadata": {},
"outputs": [],
"source": [
"data = np.loadtxt(\"data/dawsonville_slug.txt\")\n",
"to = data[:, 0]\n",
"ho = data[:, 1]"
]
},
{
"cell_type": "markdown",
"metadata": {
"editable": true,
"slideshow": {
"slide_type": ""
},
"tags": []
},
"source": [
"### Parameters and model"
]
},
{
"cell_type": "code",
"execution_count": 9,
"metadata": {},
"outputs": [],
"source": [
"# known parameters\n",
"b = 98 # aquifer thickness in m\n",
"zt = -24 # top of aquifer in m\n",
"zb = zt - b # bottom of aquifer in m\n",
"rw = 0.076 # well radius of Ln-2 Well in m\n",
"rc = 0.076 # casing radius of Ln-2 Well in m\n",
"Q = 10.16 / 1000 # slug volume in m^3 (10.16 l = 0.01016 m^3)"
]
},
{
"cell_type": "code",
"execution_count": 10,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"self.neq 1\n",
"solution complete\n"
]
}
],
"source": [
"ml = tft.ModelMaq(kaq=10, z=[zt, zb], Saq=1e-4, tmin=1e-6, tmax=1e-3, topboundary=\"conf\")\n",
"w = tft.Well(ml, xw=0, yw=0, rw=rw, rc=rc, tsandQ=[(0, -Q)], layers=0, wbstype=\"slug\")\n",
"ml.solve()"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"### Estimate aquifer parameters\n",
"\n",
"We calibrate hydraulic conductivity and specific storage, as in the KGS model (Hyder et al. 1994)."
]
},
{
"cell_type": "code",
"execution_count": 33,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"..............................\n",
"Fit succeeded.\n"
]
}
],
"source": [
"# unknown parameters: kay, Saq\n",
"cal = tft.Calibrate(ml)\n",
"cal.set_parameter(name=\"kaq0\", initial=10, pmin=0, layers=0)\n",
"cal.set_parameter(name=\"Saq0\", initial=1e-4, layers=0)\n",
"cal.seriesinwell(name=\"obs\", element=w, t=to, h=ho)\n",
"cal.fit(report=False)"
]
},
{
"cell_type": "code",
"execution_count": 13,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"| \n", " | layers | \n", "optimal | \n", "std | \n", "perc_std | \n", "pmin | \n", "pmax | \n", "initial | \n", "inhoms | \n", "parray | \n", "
|---|---|---|---|---|---|---|---|---|---|
| kaq0_0_0 | \n", "0 | \n", "0.420912 | \n", "0.018404 | \n", "4.372498 | \n", "0.0 | \n", "inf | \n", "10.0000 | \n", "None | \n", "[[0.4209121893500565]] | \n", "
| Saq0_0_0 | \n", "0 | \n", "0.000017 | \n", "0.000005 | \n", "31.179338 | \n", "-inf | \n", "inf | \n", "0.0001 | \n", "None | \n", "[[1.700498848737957e-05]] | \n", "
| \n", " | k [m/d] | \n", "Ss [1/m] | \n", "RMSE [m] | \n", "
|---|---|---|---|
| timflow | \n", "0.42 | \n", "1.70e-05 | \n", "0.004 | \n", "
| MLU | \n", "0.41 | \n", "1.94e-05 | \n", "0.004 | \n", "