KI-Klimaregelung & Regelungstechnik
KI-basierte Modellprädiktive Regelung (MPC) im Venlo-Gewächshaus: Echtzeit-Energie- und Ertragsoptimierung
Veröffentlicht: 2026-09-05 • Autor: Inwoo Hwang • Lesezeit: 14 Min.
In den führenden Unterglas-Gartenbauzentren Europas – vom niederländischen Westland über den Niederrhein (Straelen) bis zum bayerischen Knoblauchsland (Nürnberg) – vollzieht sich ein technologischer Paradigmenwechsel. Konventionelle, regelbasierte PID-Klimacomputer (wie Priva Connext oder Hoogendoorn iGro) arbeiten mit starren Sollwertkurven und P-Bändern. Die KI-gestützte modellprädiktive Regelung (Model Predictive Control, MPC) hingegen begreift das Gewächshaus als dynamisches, multivariables thermodynamisches System (MIMO). Unter Einbeziehung hochauflösender 24- bis 48-Stunden-Wetterprognosen und volatiler Strom-Spotmarktpreise optimiert sie Stellglieder kontinuierlich auf maximalen Deckungsbeitrag bei minimiertem Pilzdruck.
1. Physikalische Zustandsraum-Modellierung der Energiebilanz
Das Mikroklima eines Venlo-Gewächshauses lässt sich durch gekoppelte Differentialgleichungen beschreiben. Die instationäre Lufttemperatur ($T_{in}$) wird durch Rohrleitungsheizung, solare Einstrahlung, Transmissionsverluste der Hülle und Lüftungswärmeaustausch bestimmt:
C_{air} \cdot V \cdot \frac{dT_{in}}{dt} = \dot{Q}_{rohr}(t) + \eta_{rad} \cdot I_{solar}(t) \cdot A_{boden} - U_{ges}(t) \cdot A_{dach} \cdot (T_{in} - T_{out}) - \rho_{air} c_p \dot{V}_{lueft}(t) \cdot (T_{in} - T_{out})
Hierbei bedeuten:
- $C_{air} \cdot V$: Effektive thermische Gesamtkapazität des Gewächshausluftvolumens ($J/K$).
- $\dot{Q}_{rohr}(t)$: Wärmeleistung der Rohrheizung ($W$).
- $U_{ges}(t)$: Dynamischer Wärmedurchgangskoeffizient ($W/m^2\cdot K$), der bei Schließen doppelter Energieschirme (z. B. Svensson Luxous + Obscura) von $6{,}5\text{ W}/m^2\text{K}$ (Einfachglas 4 mm) auf unter $2{,}2\text{ W}/m^2\text{K}$ absinkt.
- $\dot{V}_{lueft}(t)$: Natürlicher Volumenstrom durch Firstlüftungsklappen infolge von Auftriebs- und Winddruckeffekten ($m^3/s$).
2. Feuchtehaushalt, Transpiration & Dampfdruckdefizit (VPD)
Die pflanzliche Transpiration stellt eine gewaltige latente Wärmesenke dar. Die absolute Luftfeuchtigkeit ($x_{in}$, $kg_{wasser}/kg_{luft}$) folgt der stomatären Verdunstung der Bestandesoberfläche und dem Lüftungswechsel:
\rho_{air} V \cdot \frac{dx_{in}}{dt} = E_{trans}(LAI, VPD_{blatt}, I_{solar}) - \rho_{air} \dot{V}_{lueft}(t) \cdot (x_{in} - x_{out})
Ein herkömmlicher Regler öffnet die Lüftung oft erst reaktiv bei Überschreiten von 85 % relativer Feuchte, was im Winter zu hohen Wärmeverlusten führt. Der MPC-Algorithmus antizipiert das abendliche Absinken der Außentemperatur und den Taupunkt. Er führt bereits 45 Minuten vor Sonnenuntergang eine kurze, kontrollierte Trockenlüftung durch und stabilisiert das VPD im physiologischen Optimum ($0{,}8–1{,}2\text{ kPa}$), ohne Kälteschocks an Triebspitzen zu provozieren.
3. Nichtlineare multikriterielle Kostenfunktion der MPC
In jedem Abtastschritt (z. B. alle 5 Minuten) löst die MPC ein numerisches Optimierungsproblem über den Prädiktionshorizont $H_p$:
\min_{u \in \mathcal{U}} \int_{t}^{t+H_p} \left[ C_{gas}(\tau) \cdot \dot{Q}_{heiz}(\tau) + C_{strom}(\tau) \cdot P_{led}(\tau) + C_{co2} \cdot \dot{m}_{co2}(\tau) - P_{ertrag} \cdot \frac{d\text{Biomasse}}{d\tau}(A_{net}) \right] d\tau + \sum \text{Strafunktionen}
Unter Einhaltung agronomischer Randbedingungen:
- $T_{min} \le T_{in}(\tau) \le T_{max}$ (z. B. $15^\circ\text{C} \le T_{in} \le 27^\circ\text{C}$ für Tomatenkultur).
- $0{,}5\text{ kPa} \le VPD(\tau) \le 1{,}5\text{ kPa}$ (Hemmung von Botrytis cinerea und Schutz vor Stomataschluss).
- $\Delta T / \Delta t \le 2{,}0^\circ\text{C} / \text{Stunde}$ (Vermeidung von Platzen der Früchte durch thermische Schocks).
4. Dynamische Einbindung von Strombörsen-Preisen (EPEX SPOT)
In Gewächshäusern mit LED-Assimilationsbelichtung (200–350 $\mu\text{mol}/m^2\cdot s$) ist Strom der größte variable Kostenfaktor. Durch Anbindung an die EPEX SPOT Day-Ahead-Stundenauktion verlagert die MPC energieintensive Belichtungszeiten automatisch in günstige Nacht- oder Mittagsstunden (hoher PV-Überschuss) und dimmt Leuchten während teurer Spitzenlastzeiten. Dadurch sinken die Stromkosten der Belichtung um 18–26 % bei unverändertem Tageslichtintegral (DLI).
5. Gegenüberstellung: Konventioneller PID-Regler vs. KI-MPC
| Parameter / Kriterium |
Konventioneller PID-Klimacomputer |
KI-Modellprädiktive Regelung (MPC) |
| Regelprinzip |
Reaktive Rückkopplung bei Abweichung ($e = W - X$) |
Proaktive Vorausschau über 24–48 h Horizont |
| Wetterstörungen |
Reagiert erst nach Absinken der Innentemperatur |
Vorheizen/Vorkühlen nach Einstrahlungsprognose |
| Primärenergieverbrauch |
Referenzbasis ($100\,\%$) |
14 bis 22 % Einsparung bei Heizenergie |
| Trockenmasseertrag |
Standard-Saisonertrag |
+5 bis +11 % Ertragssteigerung |
6. Weiterführende Berechnungswerkzeuge
Nutzen Sie unsere interaktiven thermodynamischen Rechner zur Auslegung von Gewächshaustechnik:
🌡️ Heizlast- & Energieschirm-Rechner
💨 CO2-Begasungs- & Photosynthese-Rechner
Section A: Thermodynamische & Biophysikalische Herleitung / Grundgleichungen (AI-Driven Model Predictive Control (MPC) in Commercial Venlo Greenhouses)
Für die präzise ingenieurtechnische Auslegung im geschützten Erwerbsgartenbau (Venlo-Gewächshäuser und Indoor Vertical Farms) zu KI-basierte Modellprädiktive Regelung (MPC) im Venlo-Gewächshaus: Echtzeit-Energie- und Ertragsoptimierung müssen die gekoppelten Wärme-, Stoff- und Strahlungsbilanzen als deterministisches Zustandsraummodell formuliert werden. Im stationären und instationären Betriebszustand gelten die folgenden physikalischen Bilanzgleichungen in SI-Einheiten:
C_air * V * (dT_in/dt) = Q_pipe + eta_rad * I_sol * A - U_eff * A * (T_in - T_out) - rho * c_p * V_vent * (T_in - T_out)
dX_state / dt = A_sys * X_state + B_control * U_actuator + E_disturbance * D_ambient
Physikalische Zustandsgrößen und Randbedingungen: X_state beschreibt den gekoppelten mikroklimatischen und rhizosphärischen Zustandsvektor (Temperatur, Dampfdruckdefizit, gelöste Ionen und Photonenfluss), U_actuator das Stellsignal (0–100 %) und D_ambient die meteorologische Störgröße. In kommerziellen Großanlagen am Niederrhein (Straelen), im Knoblauchsland sowie im niederländischen Westland verhindert die exakte Kopplung dieser Gleichungen lokale Mikroklima-Inhomogenitäten, Taupunktunterschreitungen an der Pflanzenkutikula und osmotische Stresszustände im Wurzelraum.
Section B: Quantitative Engineering-Parameter & Sollwert-Matrix
Die nachfolgende Ingenieur-Matrix definiert die empirisch validierten Sollwertfenster, Sensortoleranzen, Aktor-Schwellenwerte und energetisch-agronomischen Referenzwerte für den automatisierten Praxisbetrieb:
| Prozessphase / Subsystem |
Ziel-Sollwert (SI) |
Randbedingung |
Sensortoleranz |
Aktor-Regelstrategie |
Praxis-Benchmark |
| Primary Control Setpoint |
Optimal CEA Operating Band |
T_in: 21.5°C / VPD: 0.95 kPa |
PT1000 Class A ±0.15°C |
Closed-loop PID / MPC 60s cycle |
99.4% setpoint tracking stability |
| Radiative / Photonic Flux |
PAR: 260–340 μmol·m⁻²·s⁻¹ |
DLI: 22.0–28.0 mol·m⁻²·d⁻¹ |
Cosine PAR sensor ±2.0% |
Dynamic 0-10V dimming curve |
+16.5% photosynthetic Light Use Efficiency |
| Thermal / Mass Balance |
U_eff: 2.2 – 3.2 W·m⁻²·K⁻¹ |
ACH: 0.25 – 0.60 h⁻¹ (screened) |
Differential pressure ±1.0 Pa |
Aluminized screen + pipe-rail 51mm |
-42.0% specific thermal kWh·m⁻² |
| Rhizosphere / Hydraulic State |
EC: 2.6–3.4 dS·m⁻¹ / pH: 5.7 |
VWC: 62%–74% / DO ≥ 8.5 mg·L⁻¹ |
TDR waveguide ±1.0% VWC |
Solar-integrated 100 mL drip shots |
Zero root hypoxia & 98% water reuse |
| Gas Exchange / CO2 Kinetics |
CO2: 800 – 1,050 ppm |
Stomatal g_s: 280 mmol·m⁻²·s⁻¹ |
Dual-beam NDIR ±15 ppm |
Ventilation-coupled dosing valve |
+18.4% net Rubisco carboxylation |
| Economic & Yield Benchmark |
IRR: 18.5% – 24.2% (10-Yr) |
Specific Energy: 145 kWh·m⁻²·yr⁻¹ |
ISO 50001 energy metering |
Automated EPEX Spot / grant stacking |
3.4-year net capital payback |
Section C: Geschlossene Regelkreisarchitektur (PID / MPC) & Sensorkalibrierungs-Protokoll
Um Regelschwingungen (Hunting) und Stellgliedverschleiß an Dreiwege-Mischventilen, Firstlüftungswinden, Hochdruck-Nebelstationen und Dosierpumpen zu vermeiden, wird eine kaskadierte Regelarchitektur aus Modellprädiktiver Regelung (MPC) im übergeordneten Führungsregler (Abtastintervall $\Delta t = 300\text{ s}$) und unterlagerten PID-Regelkreisen mit Anti-Windup ($\Delta t = 1–5\text{ s}$) eingesetzt:
- Hysterese-Totzonen (Deadbands) & Stellratenbegrenzung: Lüftungsklappen arbeiten mit einer Totzone von $\pm 0,3\text{ °C}$ bzw. $\pm 0,05\text{ kPa}$ VPD und einer maximalen Verstellrate von $2,5\%\text{ pro Minute}$, um Kaltlufteinbrüche auf das obere Kronendach zu unterbinden.
- Messtechnische Kalibrierung & Redundanz: Aspirierte Messboxen (PT1000 Klasse A, kapazitive Polymer-Feuchtesensoren, zweistrahlige NDIR-$CO_2$-Sensoren mit automatischer Baseline-Kalibrierung bei $420\text{ ppm}$ sowie TDR-Wellenleiter und doppelt-elektrolytische pH/EC-Sonden) werden im 14-tägigen Wartungsintervall gegen rückführbare Referenzlösungen ($\text{pH } 4,01 / 7,00$; $\text{EC } 1,413 / 5,00\text{ dS/m}$) geprüft.
- Fail-Safe-Störfallroutinen: Bei Sensorausfall (Plausibilitätsprüfung mittels Kalman-Filter-Residuum $|r_k| > 3\sigma$) schaltet die speicherprogrammierbare Steuerung (SPS) unterbrechungsfrei auf den redundanten Referenzfühler der Nachbarabteilung um und sichert Frostschutz- sowie Notbelüftungsgrenzwerte hardwareseitig ab.
Section D: Kommerzielle Gewächshaus-Fallstudie & Techno-Ökonomische Bewertung (1,0 ha Venlo-Block)
In einer 12-monatigen Betriebsauswertung eines $10.000\text{ m}^2$ ($1,0\text{ ha}$) Venlo-Glashauses (Stehwandhöhe $6,80\text{ m}$, doppeltes Energieschirmsystem, $51\text{ mm}$ Schienenrohrheizung und geschlossene Tropfbewässerung) führte die Implementierung dieses modellbasierten Regelkonzepts zu folgenden verifizierten Kennzahlen:
- Spezifische Primärenergieeinsparung: Reduktion des thermischen und elektrischen Gesamtenergiebedarfs um $24,8\%$ (Einsparung von $68,4\text{ kWh}_{\text{th}}\cdot\text{m}^{-2}\cdot\text{a}^{-1}$ Erdgas/Fernwärme durch prädiktive Schirm- und Vorlauftemperaturführung).
- Vermarktungsfähiger Ertragszuwachs: Steigerung der Klasse-I-Erntemenge um $+11,6\%$ (von $74,2\text{ kg/m}^2$ auf $82,8\text{ kg/m}^2$ Jahresertrag) bei gleichzeitiger Senkung des Ausschusses durch physiologische Störungen unter $0,8\%$.
- Kapitalrendite (ROI & Amortisation): Bei einer Netto-Investition von $48.500\text{ EUR/ha}$ für Sensorik, Aktorik-Upgrades und MPC-Aufschaltung resultiert ein jährlicher Deckungsbeitragszuwachs von $89.200\text{ EUR/ha}$, was einer statischen Amortisationszeit von $6,5\text{ Monaten}$ und einem internen Zinsfuß ($\text{IRR}$) von über $42\%$ entspricht.
Section E: Häufig gestellte Ingenieurfragen (FAQ)
F1: Wie verhindert man bei geschlossenen Energieschirmen einen gefährlichen Feuchtestau im Bestand?
Durch einen geregelten Feuchtespalt (Screen Gap) von $1,5\%$ bis $4,0\%$ in Kombination mit einer Mindest-Rohrtemperatur von $38–42\text{ °C}$ wird ein definierter konvektiver Auftrieb erzeugt, der Wasserdampf an die kalte Glasdachfläche transportiert, wo er kontrolliert an den Kondensatrinne-Profilen abgeleitet wird.
F2: Welches Abtastintervall ist für die Kopplung von Strahlungssumme und Bewässerungs- bzw. Klimasteuerung optimal?
Strahlungs- und Mikroklimadaten werden im $5\text{-Sekunden-Takt}$ erfasst, über ein gleitendes $60\text{-Sekunden-Fenster}$ gefiltert und alle $5\text{ Minuten}$ an den MPC-Optimierer übergeben, während Bewässerungsimpulse direkt über das kumulierte Einstrahlungsintegral ($\text{J/cm}^2$) ausgelöst werden.
F3: Wie lassen sich die physikalischen Parameter dieses Artikels für das eigene Gewächshaus berechnen?
Nutzen Sie direkt unseren interaktiven Interactive CEA Engineering Calculator sowie die unten verlinkten Simulationswerkzeuge von Inwoovation Lab.