Control del clima del invernadero Modelo en vivo

Diriges el clima de un invernadero de vidrio de una hectárea que cultiva tomates. El sol, el aire exterior y el propio cultivo alteran cada diez minutos la temperatura, la humedad y el CO₂; tus palancas son las ventanas cenitales, la consigna de calefacción, una pantalla de sombreo/ahorro de energía, la nebulización de alta presión y la dosificación de CO₂. Cada cifra sale de un balance de energía y de masa que puedes inspeccionar más abajo.

Qué aprenderás

Simulador

Tiempo 0 min
Temperatura del aire del invernadero 20.0 °C · Humedad relativa 75% · Déficit de presión de vapor (DPV) 0.58 kPa · Concentración de CO₂ 420 ppm · Temperatura exterior 18.4 °C · Cultivo en confort☀ 18.4 °C · 58% · 371 W/m²✓✓✓✓✓✓✓✓✓✓🌡 20.0 °C · 💧 75% · ♨ 0 W/m² · 🌿 0%CO₂ 420 ppmVPD 0.58 kPa♨ Calor 0.00 kWh/m²💧 Nebulización 0.00 L/m²CO₂ dosificado 0 g/m²Horas de estrés: calor · frío · humedad🔥 0.0 h · ❄ 0.0 h · 🍄 0.0 h
  • Cultivo en confort
  • Cultivo con estrés leve
  • Cultivo en riesgo (calor, frío o enfermedad)
  • Gotas de nebulización
  • Tuberías de calefacción en marcha

Controles

El intercambio de aire crece con la apertura y el viento. Elimina calor, humedad y CO₂ por igual.

La caldera calienta hasta 180 W/m² siempre que el aire esté más frío que este valor.

Cerrada: bloquea el 60 % de la luz y reduce la pérdida de calor por la cubierta de 6,2 a 4,0 W/m²K, pero también dificulta el intercambio de aire y atrapa la humedad debajo.

Hasta 0,4 L/m² por hora. Las gotas al evaporarse enfrían el aire; con aire húmedo mojan las hojas.

Hasta 8 g/m² por hora. Solo merece la pena mientras las ventanas cenitales están casi cerradas.

Indicadores

Temperatura del aire del invernadero
20.0°C
normal
Humedad relativa
75%
normal
Déficit de presión de vapor (DPV)
0.58kPa
normal
Fotosíntesis del cultivo (relativa)
0%
crítico
Temperatura exterior18.4 °C
Humedad exterior58 %
Luz solar exterior371 W/m²
Concentración de CO₂420 ppm
Margen sobre el punto de rocío4.6 K
Potencia de calefacción0 W/m²
Transpiración del cultivo0.00 L/m²·h
Renovaciones de aire1 1/h
Energía de calefacción usada0.00 kWh/m²
Agua de nebulización usada0.00 L/m²
CO₂ dosificado0 g/m²
CO₂ fijado por el cultivo0.0 g/m²
Horas de estrés por calor (> 35 °C)0.0 h
Horas de estrés por frío (< 10 °C)0.0 h
Horas de riesgo de enfermedad (HR ≥ 90 % u hojas mojadas)0.0 h
Grados-día de crecimiento (base 10 °C)0.00 °C·d

Tendencia

Temperatura del aire del invernadero: — °C45.0-5.0

Escenarios de crisis

Nivel 1 · Tarde de ola de calor

Son las diez de la mañana, cielo despejado, y el pronóstico anuncia 33 °C en el exterior a media tarde. Las ventanas están abiertas solo un tercio y el invernadero ya se está calentando. Evita que el cultivo sufra estrés por calor durante la tarde sin malgastar agua.

  • Tiempo de estrés por calor por encima de 35 °C ≤ 15 min
  • Agua de nebulización ≤ 1,8 L/m²
  • DPV medio ≤ 2,2 kPa
  • Tiempo con hojas mojadas o muy húmedo ≤ 30 min

Nivel 2 · Noche fría y despejada

Las seis de la tarde, cielo sin una nube y se esperan −4 °C antes del amanecer. Siguen los ajustes de día: calefacción a 18 °C, pantalla abierta y ventanas entreabiertas. Mantén el cultivo por encima de 10 °C toda la noche con el mínimo combustible y devuélvele la luz al amanecer.

  • Tiempo por debajo de 10 °C ≤ 15 min
  • Energía de calefacción ≤ 1,4 kWh/m²
  • Fotosíntesis media tras el amanecer ≥ 12 %

Nivel 3 · Día nublado y húmedo

Un día gris y húmedo: 9 °C fuera, el aire exterior casi saturado y el invernadero cerrado desde anoche a 14 °C. El cultivo sigue transpirando y la humedad sube hacia condiciones de moho gris. Mantén bajas las horas húmedas y el aire al menos a 15,5 °C sin agotar tu presupuesto de calefacción, y no tapes la poca luz que hay.

  • Horas de riesgo de enfermedad ≤ 2 h
  • Energía de calefacción ≤ 0,8 kWh/m²
  • El aire nunca baja de 15,5 °C
  • Fotosíntesis media ≥ 32 %

Base: el modelo detrás de las cifras

Todas las relaciones que usa el simulador, con su fuente. Las constantes marcadas como supuestos son calibraciones ilustrativas.

La temperatura del aire sigue un balance de energía: entra el sol, sale calor por la cubierta y con el aire de ventilación, se gasta calor en evaporar agua y se suma la calefacción.
C·dT/dt = τ·I·(1 − 0.6·0.75·screen) − λ·(E_crop + E_fog) + Q_heat + (U·Ac/Af + ρc_p·Q_v)·(T_out − T) + G_soil·(T_soil − T)[1]Supuesto: el tamaño del invernadero, la capacidad térmica, las ventanas, la pantalla, la nebulización y los coeficientes del cultivo son valores ilustrativos de un invernadero de vidrio moderno; la física es un único volumen de aire bien mezclado. Las instalaciones reales calibran cada constante.
La pérdida de calor por la cubierta es U × superficie × diferencia de temperatura; una pantalla cerrada baja U, un cielo nocturno despejado la sube.
U = 6.2 − (6.2 − 4.0)·screen W/m²K (single glass → with thermal blanket); ×(1 + 0.25·clear sky) at night[1]Supuesto: el tamaño del invernadero, la capacidad térmica, las ventanas, la pantalla, la nebulización y los coeficientes del cultivo son valores ilustrativos de un invernadero de vidrio moderno; la física es un único volumen de aire bien mezclado. Las instalaciones reales calibran cada constante.
La ventilación es la suma de las fugas y del flujo, impulsado por el viento, a través de las ventanas abiertas.
Q_v = (N·H/3600 + vent·C_w·max(1, u_wind))·(1 − 0.4·screen), N = 0.75 h⁻¹[1][4]Supuesto: el tamaño del invernadero, la capacidad térmica, las ventanas, la pantalla, la nebulización y los coeficientes del cultivo son valores ilustrativos de un invernadero de vidrio moderno; la física es un único volumen de aire bien mezclado. Las instalaciones reales calibran cada constante.
La presión de vapor de saturación aumenta con fuerza con la temperatura; el DPV indica cuánto le falta al aire para saturarse.
e_s(T) = 0.6108·exp(17.27·T / (T + 237.3)) kPa; VPD = e_s(T) − e_a; RH = e_a / e_s(T)[2]
Balance de vapor de agua: la transpiración y la nebulización aportan humedad; la ventilación y la condensación en el vidrio frío la retiran.
H·dχ/dt = E_crop + E_fog − Q_v·(χ − χ_out) − k_c·(Ac/Af)·(1 − 0.7·screen)·max(0, e_a − e_s(T_cover)), χ = e·M_w/(R·T)[6][2]Supuesto: el tamaño del invernadero, la capacidad térmica, las ventanas, la pantalla, la nebulización y los coeficientes del cultivo son valores ilustrativos de un invernadero de vidrio moderno; la física es un único volumen de aire bien mezclado. Las instalaciones reales calibran cada constante.
La transpiración tiene una parte debida a la radiación y otra impulsada por el aire seco (DPV).
E_crop = A·(1 − e^(−k·LAI))·I_in/λ + B·LAI·VPD (A 0.45, k 0.64, LAI 3, B 0.02 kg/m²·h·kPa by day)[3]Supuesto: el tamaño del invernadero, la capacidad térmica, las ventanas, la pantalla, la nebulización y los coeficientes del cultivo son valores ilustrativos de un invernadero de vidrio moderno; la física es un único volumen de aire bien mezclado. Las instalaciones reales calibran cada constante.
La nebulización solo enfría en la medida en que se evapora: unos 2,45 MJ por litro, más con aire seco que con aire húmedo.
E_fog = fog·0.4 kg/m²·h × min(1, VPD/0.8); each evaporated kg removes λ = 2.45 MJ of sensible heat; the rest wets leaves[1][2]Supuesto: el tamaño del invernadero, la capacidad térmica, las ventanas, la pantalla, la nebulización y los coeficientes del cultivo son valores ilustrativos de un invernadero de vidrio moderno; la física es un único volumen de aire bien mezclado. Las instalaciones reales calibran cada constante.
Balance de CO₂: entra por la dosificación; salen la ventilación y la absorción del cultivo.
H·dC/dt = dose − Q_v·(C − C_out) − P_net; dosing ≈ 5.5 g/m²·h holds ≈ 1000 ppm at 1 air change per hour[1]
La fotosíntesis responde a la luz, el CO₂, la temperatura y el DPV; la respuesta al CO₂ sigue una regla empírica publicada para cultivos de fruto en invernadero.
P = P_ref·(1 − e^(−I_in/300))·exp(150/400 − 150/C)·f(T: 1 at 17–27 °C)·f(VPD) (from the CO₂ rule dP/P = 1.5 %·(1000/C)² per 100 ppm)[5][7]Supuesto: el tamaño del invernadero, la capacidad térmica, las ventanas, la pantalla, la nebulización y los coeficientes del cultivo son valores ilustrativos de un invernadero de vidrio moderno; la física es un único volumen de aire bien mezclado. Las instalaciones reales calibran cada constante.
Contadores de estrés del cultivo: horas con demasiado calor, demasiado frío o con humedad suficiente para el moho gris.
heat hours: T > 35 °C · cold hours: T < 10 °C · disease-risk hours: RH ≥ 90 % or wet leaves · GDD = Σ max(0, T − 10)·Δt[1][7][8][10]
Clima exterior: curva diaria de temperatura, punto de rocío constante, radiación solar reducida por las nubes y fluctuaciones aleatorias de nubosidad y viento.
T_out = T̄ + A·cos(2π(h − 15)/24); e_out = e_s(T_dew); I = I_clear·(1 − 0.75·N^3.4); cloud N and wind fluctuate (normal, seeded)[9][2]Supuesto: el tamaño del invernadero, la capacidad térmica, las ventanas, la pantalla, la nebulización y los coeficientes del cultivo son valores ilustrativos de un invernadero de vidrio moderno; la física es un único volumen de aire bien mezclado. Las instalaciones reales calibran cada constante.
Otras constantes que usa el modelo.
height 4.5 m · cover/floor 1.2 · τ 0.8 · heat capacity 30 kJ/m²K · heating ≤ 180 W/m² · screen blocks 60 % of light (75 % of that as heat) and 40 % of air exchange, and hides 70 % of the cold glass (condensation) · vents 0.016 m³/s·m² per m/s wind · soil 2 W/m²K · fog ≤ 0.4 L/m²·h · CO₂ ≤ 8 g/m²·h · P_ref 4 g CO₂/m²·h · condensation 1.9·10⁻⁸ kg/m²·s·PaSupuesto: el tamaño del invernadero, la capacidad térmica, las ventanas, la pantalla, la nebulización y los coeficientes del cultivo son valores ilustrativos de un invernadero de vidrio moderno; la física es un único volumen de aire bien mezclado. Las instalaciones reales calibran cada constante.

Aleatoriedad: un generador mulberry32 con semilla; distribuciones usadas: uniforme, exponencial (CDF inversa), normal (Box–Muller), Poisson (Knuth). La semilla se muestra y se puede compartir.

Fuentes

  1. ANSI/ASAE EP406.4 JAN2003 (R2008) — Heating, Ventilating and Cooling Greenhouses (eq. 1–4, Table 1 U-values, Table 2 infiltration: new glass or fiberglass 0.50–1.0 h⁻¹, §5.8.5 CO₂ injection, §7.1.3 temperature limits); since withdrawn, available from ASABE as a historical document — American Society of Agricultural and Biological Engineers, 2003
  2. R. G. Allen, L. S. Pereira, D. Raes, M. Smith — Crop evapotranspiration, FAO Irrigation and Drainage Paper 56 (eq. 11 saturation vapour pressure, eq. 14 dew point, λ = 2.45 MJ/kg) — FAO, Rome, 1998
  3. M. Baille, A. Baille, J. C. Laury — A simplified model for predicting evapotranspiration rate of nine ornamental species vs. climate factors and leaf area — Scientia Horticulturae 59:217–232, 1994
  4. T. Boulard, A. Baille — Modelling of air exchange rate in a greenhouse equipped with continuous roof vents — Journal of Agricultural Engineering Research 61(1):37–48, 1995
  5. E. M. Nederhoff — Effects of CO₂ concentration on photosynthesis, transpiration and production of greenhouse fruit vegetable crops (PhD thesis; CO₂ rule) — Wageningen Agricultural University, 1994
  6. B. H. E. Vanthoor, C. Stanghellini, E. J. van Henten, P. H. B. de Visser — A methodology for model-based greenhouse design: Part 1, a greenhouse climate model for a broad range of designs and climates — Biosystems Engineering 110:363–377, 2011
  7. R. R. Shamshiri et al. — Review of optimum temperature, humidity, and vapour pressure deficit for microclimate evaluation and control in greenhouse cultivation of tomato: a review — International Agrophysics 32:287–302, 2018
  8. Identifying and Managing Three Major Disease Threats in Greenhouse Tomatoes: Leaf Mold, Gray Mold, and Bacterial Wilt (SPES-820) — Virginia Cooperative Extension
  9. F. Kasten, G. Czeplak — Solar and terrestrial radiation dependent on the amount and type of cloud — Solar Energy 24:177–189, 1980
  10. G. S. McMaster, W. W. Wilhelm — Growing degree-days: one equation, two interpretations — Agricultural and Forest Meteorology 87:291–300, 1997

Quién se dedica a esto

Modelo educativo: no sirve para decisiones operativas. Las instalaciones reales calibran cada constante con sus propios equipos y datos.