← Back to All 38 Engineering Tools
⚡ Flagship Tool 38 • Decarbonized Geothermal CEA

Aquifer Thermal Energy Storage (ATES) & Heat Pump Sizing Simulator

Model commercial greenhouse seasonal doublet well hydraulics, Darcy drawdown, VDI 4640 seasonal heat pump COP, and German BEW (Bundesförderung für effiziente Wärmenetze) / US IRA 48C decarbonization economics.

🏛️ 1. Greenhouse & Peak Thermal Demand

Cultivation Area (Hectares) 2.5 ha
Peak Heating Demand (W/m²) 140 W/m²
Peak Cooling Demand (W/m²) 110 W/m²
Annual Full-Load Heating Hours 1,650 h/yr

🌊 2. Aquifer Hydrogeology (VDI 4640)

Aquifer Saturated Thickness (H_aq) 25 m
Hydraulic Conductivity (k_f × 10⁻⁴ m/s) 2.5 × 10⁻⁴ m/s (Sand/Gravel)
Warm Well Storage Temp (T_warm) 19.0 °C
Cold Well Storage Temp (T_cold) 7.5 °C
Well Separation Distance (L_well) 130 m

💶 3. Heat Pump & Energy Economics

Heat Pump Supply Water Temp 42 °C (Low-Temp Rail)
Electricity Tariff (€/kWh) 0.14 €/kWh
Natural Gas Baseline (€/kWh) 0.085 €/kWh
BEW / IRA Clean Energy Subsidy 45% Non-Dilutive Grant
Pumping Flow Rate
145
m³/h (Doublet Peak)
Heating COP (VDI 4640)
5.12
W_th / W_el (Seasonal)
Well Drawdown
3.4
m (Dupuit-Theis)
Thermal Recovery (η_th)
84.2
% (No Breakthrough)
Annual CO₂ Reduction
842
Tons CO₂e / Year
Annual OPEX Savings
€164,200
Net Annual Profit
Subsidized Payback
3.4
Years (BEW 45% Grant)

📐 Hydrogeological & Thermodynamic Formulation

Groundwater flow and seasonal heat recovery obey the coupled Darcy groundwater potential and energy conservation laws:

Pumping Flow: V_dot = Q_thermal / (rho_w * C_p * Delta_T_well) [m³/h]
Darcy Well Drawdown: s_w = V_dot / (2 * pi * k_f * H_aq) * ln(R_influence / r_well) [m]
Carnot COP: COP_th = eta_carnot * (T_supply_K / (T_supply_K - T_evap_K))
Thermal Radius: R_thermal = sqrt((C_w * V_seasonal) / (pi * H_aq * C_aquifer)) [m]
VDI 4640 Compliance Status:
✅ Excellent: Well distance (130m) exceeds minimum safe thermal distance (94m). Zero thermal breakthrough risk over 25-year operational lifecycle.

📚 Theoretical Hydrogeology, VDI 4640 Standards & Numerical Verification

Aquifer Thermal Energy Storage (ATES) utilizes saturated geological formations (predominantly unconsolidated sand and gravel aquifers confined by clay aquicludes) as seasonal thermal energy buffers. During the summer cooling cycle, chilled groundwater ($6\text{--}9^\circ\text{C}$) is extracted from the Cold Well, routed through plate heat exchangers to absorb excess greenhouse sensible heat and fan-coil loads, and re-injected into the Warm Well at $16\text{--}22^\circ\text{C}$. In winter, the flow direction reverses: warm water is pumped from the Warm Well into the evaporator of industrial ammonia ($\text{NH}_3$/R717) or hydrocarbon heat pumps, lifting heat to $45\text{--}55^\circ\text{C}$ for pipe-rail heating networks while returning chilled water to recharge the cold reservoir.

1. Dupuit-Theis Hydraulic Drawdown & Well Interference

For a fully penetrating well in a confined aquifer of thickness $H_{\text{aq}}$ and hydraulic conductivity $k_f$, steady-state radial drawdown $s_w$ at well radius $r_w$ under volumetric extraction $\dot{V}$ is governed by the Dupuit-Theis formulation:

s_w = \frac{\dot{V}}{2 \pi \cdot k_f \cdot H_{\text{aq}}} \cdot \ln\left(\frac{R_{\text{influence}}}{r_w}\right) \quad [\text{m}]

To prevent catastrophic hydraulic jamming and borehole screen siltation, Darcy entry velocity must not exceed $v_{\text{crit}} = \sqrt{k_f} / 30 \text{ m/s}$ (Sichardt criteria).

2. Thermal Plume Geometry & Breakthrough Prevention (VDI 4640 Part 4)

Under radial advective transport, the seasonal thermal radius $R_{\text{th}}$ of injected warm or cold water within a porous matrix having volumetric heat capacity $c_{\text{aq}} = (1-\phi)c_{\text{solid}} + \phi c_{\text{water}}$ is calculated as:

R_{\text{th}} = \sqrt{\frac{c_w \cdot V_{\text{seasonal}}}{\pi \cdot H_{\text{aq}} \cdot c_{\text{aq}}}} \quad [\text{m}]

According to VDI 4640 guidelines, well separation distance $L_{\text{well}}$ must strictly satisfy $L_{\text{well}} \ge 2.5 \cdot R_{\text{th}} + v_{\text{regional}} \cdot \tau_{\text{seasonal}}$ to prevent early thermal breakthrough, which degrades seasonal coefficient of performance (COP) over multi-decade lifecycles.

3. Frequently Asked Questions (FAQ)

Q1: What are the primary hydrogeological requirements for an agricultural ATES system?

An unconsolidated sand or gravel aquifer with transmissivity $T = k_f \cdot H_{\text{aq}} > 50\text{ m}^2/\text{day}$, natural hydraulic gradient $< 0.005$ to avoid plume drift, low dissolved iron ($< 0.2\text{ mg/L}$) and manganese to prevent well screen clogging, and confining upper/lower clay aquicludes to isolate potable drinking water aquifers.

Q2: How does ATES compare to conventional closed-loop borehole heat exchangers (BHE)?

For commercial facilities exceeding 1 hectare (heating loads $> 1.5\text{ MW}_{\text{th}}$), ATES achieves $40\text{--}60\%$ lower capital expenditure per kilowatt installed compared to drilling hundreds of 100m U-tube boreholes, while offering significantly higher flow rates ($100\text{--}300\text{ m}^3/\text{h}$) and direct free cooling capacities.

Q3: How do German BEW (Bundesförderung effiziente Wärmenetze) grants apply?

Under BEW Modul 2, innovative deep and shallow geothermal installations, industrial heat pumps, and aquifer doublet wells qualify for up to $40\text{--}45\%$ non-dilutive capital grants, shortening commercial payback periods from 6.8 years to under 3.5 years.

Q4: What maintenance is required to prevent well clogging and geochemical scaling?

Nitrogen gas blanketing in buffer tanks, pressurized closed-loop hydraulic loops keeping water above $\text{CO}_2$ outgassing pressure, seasonal regenerative backwashing, and continuous redox potential monitoring prevent calcium carbonate precipitation and ferric hydroxide scaling.