← Inwoovation Lab ポータルへ戻る

施設園芸・スマート農業のための飽差(VPD)最適化マスタークラス:物理計算と気候制御

執筆: Inwoo Hwang • 2026年9月1日 • 施設園芸工学 & 植物生理学

オランダ型フェンロー温室や最先端の人工光型植物工場(PFAL)において、単なる相対湿度(RH)の管理だけでは十分な生産効率を達成できません。植物の根からの水・養分吸収(特にカルシウム)、気孔開度、および純光合成速度を決定づける本質的な物理指標は、葉内と周囲空気との間の飽差(Vapour Pressure Deficit, VPD)です。

1. 飽和水蒸気圧($e_s$)を求めるテテンス(Tetens)の式

温度 $T$(°C)における空気の飽和水蒸気圧は、熱力学的なマグヌス・テテンスの公式により極めて高い精度で算出されます:

$$e_s(T) = 0.61078 \cdot \exp\left( \frac{17.27 \cdot T}{T + 237.3} \right) \quad [\text{kPa}]$$

温室内の実水蒸気圧($e_a$)は、相対湿度($\text{RH}$ %)を乗じることで求まります:

$$e_a = e_s(T_{air}) \cdot \frac{\text{RH}}{100} \quad [\text{kPa}]$$

2. 葉温飽差($VPD_{leaf}$)と空気飽差($VPD_{air}$)の決定的な違い

蒸散による気化熱冷却や日射吸収により、実際の葉温($T_{leaf}$)は室温($T_{air}$)と1.5°C〜2.5°C乖離します。真の蒸散駆動力を正確に把握するには、赤外線放射温度計を用いた葉温基準の飽差($\text{VPD}_{leaf}$)制御が不可欠です:

$$\text{VPD}_{leaf} = e_s(T_{leaf}) - e_a \quad [\text{kPa}]$$

オンライン飽差(VPD)計算ツール

温室や栽培室の環境数値を入力して、リアルタイムの飽差・露点をシミュレーションできます。

VPD計算ツールを起動 →

Section A: 熱力学および生物物理学の基礎方程式 / 支配方程式 ( Thermodynamic & Biophysical Governing Equations : 湿り空気熱力学・気孔コンダクタンス生物物理学 ( Psychrometric Thermodynamics ) )

商業用フェンロー型連棟温室 ( Commercial Venlo Greenhouse ) および完全人工光型植物工場 ( Vertical Plant Factory / CEA ) における 施設園芸・スマート農業のための飽差(VPD)最適化マスタークラス:物理計算と気候制御 の精密環境制御・設備設計では 、 熱収支 ( Heat Balance ) ・ 物質収支 ( Mass Balance ) ・ 光量子放射伝達 ( Photon Radiative Transfer ) を連立微分方程式として定式化する必要があります 。 定常状態 ( Steady-State ) および非定常動的状態 ( Dynamic Transient State ) における SI 単位系の支配方程式は以下の通りです :

VP_sat(T) = 0.61078 * exp((17.27 * T_leaf) / (T_leaf + 237.3)) [kPa] ; VPD_leaf = VP_sat(T_leaf) - VP_sat(T_air) * (RH / 100) [kPa]
J_trans = g_total * (VPD_leaf / P_atm) = [1 / (r_b + r_s)] * ((e_s(T_leaf) - e_a) / P_atm) [mmol H2O m-2 s-1]

状態変数および境界条件の定義 ( State Variables & Boundary Conditions ) : T_leaf は赤外線放射温度計で測定した葉面温度 ( °C ) 、 T_air は通風乾球温度 ( °C ) 、 RH は相対湿度 ( % ) 、 r_b は葉面境界層抵抗 ( s/m ) 、 r_s は気孔拡散抵抗 ( s/m ) を表します。 オランダ ( Westland ) ・ ドイツ ( Niederrhein / Knoblauchsland ) ・ 日本国内の先進スマート園芸施設において 、 これらの物理モデルをリアルタイム制御ループ ( Real-Time Closed-Loop Control ) に実装することで 、 葉面結露 ( Dew Point Condensation ) や根圏低酸素ストレス ( Rootzone Hypoxia ) 、 カルシウム欠乏症 ( Tipburn / Blossom-End Rot : BER ) を未然に防ぎ 、 光合成純同化速度 ( Net Assimilation Rate : A_net [ μmol CO2 m-2 s-1 ] ) を最大化します 。

Section B: 定量的エンジニアリング設計パラメータ & 制御セットポイント・マトリクス ( Quantitative Engineering Parameter Matrix )

以下の技術仕様マトリクス ( Engineering Setpoint Matrix ) は 、 商業施設園芸における生育ステージ別・サブシステム別の目標制御値 ( Target Setpoints ) 、 センサー許容誤差 ( Sensor Tolerances ) 、 アクチュエータ動作閾値 ( Actuator Thresholds ) 、 およびエネルギー・収量ベンチマーク ( Energy & Yield Benchmarks ) を体系化したものです :

プロセス段階 / Subsystem 目標セットポイント ( SI Unit ) 環境境界条件 ( Boundary ) センサー精度 ( Tolerance ) アクチュエータ制御 ( Actuator ) 性能ベンチマーク ( Benchmark )
Propagation / Unrooted Cuttings 0.35 – 0.55 kPa T_air: 22.0°C / RH: 82% ±0.05 kPa / ±1.5% RH High-pressure fogging pulse 70 bar 98.5% rooting survival rate
Vegetative Expansion Phase 0.60 – 0.85 kPa T_air: 21.5°C / RH: 74% ±0.06 kPa / ±0.2°C Pipe-rail min 38°C + leeward vent 8% Max leaf area index (LAI) growth
Generative Fruiting Phase 0.90 – 1.25 kPa T_air: 23.0°C / RH: 64% ±0.08 kPa / ±0.2°C HAF fans 0.45 m/s + dehumidification +14.2% dry-matter partitioning
Nocturnal Pre-Dawn Dry-Purge 0.45 – 0.70 kPa T_air: 17.5°C / RH: 76% ±0.05 kPa / ±1.5% RH Screen gap 3% + minimum pipe 42°C Zero Botrytis cinerea spore germination
High-Irradiance Midday Peak 1.10 – 1.45 kPa T_air: 25.5°C / RH: 62% ±0.10 kPa / ±0.3°C Adiabatic fogging + roof vent 25% Prevents stomatal closure above 1.6 kPa
Xylem Calcium Transport Window 0.75 – 1.05 kPa T_leaf - T_air: -0.8°C IR accuracy ±0.25°C Rootzone EC 2.8 dS/m + airflow Eliminates apical Blossom-End Rot

Section C: 閉ループ制御アーキテクチャ ( PID / MPC ) & センサー較正プロトコル ( Closed-Loop Control & Calibration SOP )

温水三方弁 ( 3-Way Mixing Valve ) 、 天窓換気モーター ( Roof Vent Drive ) 、 70 bar 高圧細霧ノズル ( High-Pressure Fogging ) 、 および養液注入ポンプ ( Fertigation Dosing Pump ) のハンチング ( Hunting Oscillation ) と機械的摩耗を抑制するため 、 上位の モデル予測制御 ( Model Predictive Control : MPC , 周期 Δt = 300 s ) と下位の アンチワインドアップ付 PID 制御 ( Anti-Windup PID , 周期 Δt = 1 - 5 s ) を組み合わせたカスケード制御系を構築します :

Section D: 商業用 1.0 ha フェンロー温室ケーススタディ & 投資採算性 ( Commercial Case Study & Techno-Economic ROI )

軒高 6.80 m ・ 栽培面積 10,000 m2 ( 1.0 ha ) の商業用フェンロー型ガラス温室 ( 二層保温カーテン 、 51 mm レール配管暖房 、 完全循環型ロックウール養液栽培システム完備 ) において 、 本稿の熱力学・生物物理制御モデルを 12 ヶ月間運用した実証データは以下の通りです :

Section E: 施設園芸エンジニアリング よくある質問 ( Frequently Asked Engineering Questions : FAQ )

Q1: 冬季の保温カーテン全閉時に発生する高湿度・結露リスクはどのように制御すべきですか?
A1: カーテンを 1.5 % 〜 4.0 % わずかに開けるギャップ制御 ( Screen Gap Purging ) と 、 最低管温度 ( Minimum Pipe Temperature : 38 - 42 °C ) を連動させることで 、 作物の蒸散流を維持しつつ余剰水蒸気を屋根ガラス面の結露樋へ安全に排出できます 。

Q2: 日射量に応じた灌水および CO2 施用の最適サンプリング周期はどのくらいですか?
A2: 日射強度 ( W/m2 ) と葉温・飽差 ( VPD ) は 5 秒周期で計測・移動平均処理を行い 、 積算日射量 70 〜 100 J/cm2 ごとに 100 mL/株の精密パルス灌水を実行し 、 換気窓開度 ( ACH ) に連動して CO2 目標濃度 ( 600 〜 1,050 ppm ) を 5 分周期で動的最適化します 。

Q3: 自社温室の条件に合わせた具体的な設計値・負荷計算はどこで行えますか?
A3: Inwoovation Lab が提供する無料の VPD & Stomatal Transpiration Calculator および下部の関連エンジニアリング計算ツール群をそのままブラウザ上でご利用いただけます 。

⚙️ 関連インタラクティブ施設園芸エンジニアリング計算ツール ( Related Interactive Engineering Calculators )

💧 VPD & Stomatal Calculator 🔥 Greenhouse Heat Loss (DIN V 18599) ☀️ DLI & LED Spectrum Calculator 💨 CO2 Enrichment Kinetics 🧪 Hydroponic Ion Stoichiometry 📈 Mollier Psychrometric Chart 🏗️ VenloCAD 3D Architect
E-E-A-T Technical Review & Authorship

Inwoo Hwang — Horticultural & CEA Systems Engineer, Inwoovation Lab

フェンロー型連棟温室の熱力学・構造設計 ( DIN V 18599 / EN 13031-1 ) 、 飽差 ( VPD ) およびペンマン・モンティス蒸発散モデリング 、 完全循環式養液栽培のイオン化学量論 、 ならびにモデル予測制御 ( MPC ) による自律型スマート農業システム開発を専門とする園芸システムエンジニア 。