综述：植物叶片气孔系统的仿生热力学模型：能量转化与CO?代谢

《Agriculture》：Bionic Thermodynamic Model of the Stomatal System of a Plant Leaf: Energy Transformation and CO2 Metabolism

　　

植物叶片气体交换机制及其与能量转化的关系在定量层面仍未得到充分研究，这限制了对CO?气体交换过程的评估。本研究将植物叶片的气孔系统作为微尺度、宏尺度或纳尺度热力学系统进行分析，并应用基于仿生学原理的建模方法。构建了一个理想化热力学循环，用以评估热能向机械功的转化。热效率系数的理论上限估计为η_t ≈ 0.003，在叶片与环境温度差约为1 °C时，机械能通量可达0.6 W/m2，即约占叶片吸收太阳辐射通量的0.3%。研究结果表明，即使效率较低，这种能量转化也足以影响气体交换的强度。基于该模型，还估算了植物CO?吸收潜力与冠层面积的关系。研究得出结论：植物叶片的气孔系统可被解释为一种理论上的仿生能量转化模型，适用于CO?交换过程分析以及仿生微尺度、宏尺度或纳尺度系统的开发。

本文对植物叶片气孔系统的仿生热力学模型进行了系统分析。以下按论文主体结构对各部分内容进行总结。

## 1. 引言

全球气温上升是气候变化的主要指标，与温室气体(GHG)排放密切相关，其中CO?占比最大。通过国际协议推动清洁技术和能源效率提升虽有助于减排，但仍不足以应对气候变化，因此CO?捕获与封存技术也在持续发展。自然界中从大气中去除CO?最有效且能耗最低的机制仍是植物的自然光合作用。热能通过循环发生的热力学过程转化为机械能，常规热机在高温高压下运行，循环效率直接取决于这些参数。根据热力学第二定律，无补偿地将低势热能转化为其他形式的能量几乎不可能，因此植物叶片中仅几度的温度差能否实现这种转化需要专门研究。在植物中，同化过程直接依赖于高效的气体交换，CO?通过叶片气孔系统被吸收，水蒸气和代谢产物在叶片与环境之间运输，这些过程受气孔导度、叶肉扩散特性、叶片解剖结构以及环境因素引起的梯度调控。然而，克服气体通过这些通道运动阻力的机械力的来源仍未被解释，这表明叶片内部存在将低势热能转化为机械能的机制。叶片吸收的大部分未被生化过程利用的太阳辐射以热量形式耗散，直射太阳辐射可使叶片温度显著超过环境气温，但蒸腾等生物调节机制限制了长期的临界过热。叶片与环境温度之间不断更新的差异为热能向机械能转化创造了先决条件。植物叶片特定的解剖结构通过气孔和叶肉组织中的细胞间隙、微米级、毫米级和纳米级通道系统实现高效气体运动。已有研究指出气孔系统中发生的热力学过程允许将低势热能转化为机械能。本研究构建了一个理想化的等压-绝热循环，得到了热效率系数的解析表达式，并将其与冠层尺度的CO?吸收评估联系起来。该模型不取代气孔调控理论，而是通过定量评估低势热能向机械能的转化对现有理论进行补充。

## 2. 仿生叶片的结构与物理模型

在技术装置中，热能通过基于循环热力学过程的热机转化为机械功，类似的循环热质交换原理可应用于植物叶片的气孔系统，其中温度、压力和体积的变化为能量转化创造了条件。该热力学过程在性质上类似于喷气发动机的工作循环，涉及一系列体积变化以及温度和压力的变动，从而引起植物体内物质运动和能量转化。植物叶片在白天日照和温度脉动的影响下发生能量转化过程，这些过程与CO?同化和气体交换直接相关。由于生物和物理过程同时发生，准确描述叶片中热能转化为机械能的热力学过程存在困难，因此可将叶片的气孔和叶肉通道系统作为微尺度、宏尺度和纳尺度的热力学结构进行分析。叶片的光合机构由上表皮和下表皮之间的叶肉组成，叶肉包括栅栏组织和海绵组织，栅栏组织适应光合作用并拥有大部分叶绿体，海绵组织形成细胞间隙系统。叶片厚度取决于光照条件，例如欧洲山毛榉(Fagus sylvatica)叶片厚度范围为90至210 μm。叶片两侧由带有角质层的表皮保护，该结构对水蒸气和气体几乎不可渗透，叶片内部空气连同CO?仅通过气孔吸收。气孔系统的几何特征，包括气孔开度面积、密度和内部通道尺寸，直接决定了气体交换腔体中体积和压力变化的限度。每平方毫米叶片面积上气孔数量可达450个，气孔开度面积在不同植物中变化范围为0.17至239 μm2，气孔总面积占叶片表面的0.52%–5.28%。约90%–95%的CO?通过气孔进入叶片，仅5%–10%通过表皮和角质层扩散。气孔的行为直接影响气体交换，以优化碳获取与水分散失之间的平衡。尽管气孔大小和密度对气体交换的影响已被建模，但其与叶片中发生的热-机械能转化的关系在很大程度上仍未得到探索，正是这些形态特征决定了气孔系统作为微尺度热机的运行限度。

## 3. 植物叶片气体交换系统中的能量转换循环

CO?、O?和水蒸气通过叶片气体交换通道的运动遇到内部阻力，这种运动需要能量转化为机械功，以产生驱动力确保气体运输。根据叶片形态结构和气体交换建模，已确定叶片与大气之间的气体交换可由机械功促进，该机械功由按热机原理运行的生物过程产生，即叶片与环境之间的气体交换机械力由叶片中运行的发动机（热机的生物原型）产生的机械功所创造。叶片吸收太阳辐射能产生的热能在细胞微结构中转化为机械能，这使评估植物中可能存在微尺度热机成为可能。该过程不可避免地依赖于热力学第二定律。植物叶片中所有生物和能量代谢过程都与局部叶片温度密切相关。在典型晴天，叶片经历脉动的短波辐射通量Q_sp和脉动的对流热通量Q_conv，这些波动的热通量导致叶片中积累热量Q_ac的脉动，能量平衡表达式为：±Q_ac ± Q_sp ± Q_conv = 0。叶片中热积累的脉动过程引起温度波动，进而决定气体交换腔体中的压力变化。实验模型表明，温度仅增加1 °C（例如19 °C到20 °C），压力可增加到347 Pa。在热力学循环分析中，假设恒定均匀的气体质量参与循环，即位于叶片气体腔体中的绝对干燥空气；腔体体积视为恒定(V = const)；循环期间环境压力保持恒定(p₀ = const)。随着叶片温度变化，气体腔体中空气质量体积也发生变化。分析中引入"虚拟膜"概念，将叶片气体交换腔体条件性分为两部分：腔体V₁中的气体体积直接取决于叶片组织温度，另一侧腔体V₂中的气体通过气孔开口向外界环境运动。该热力学循环在p–V坐标系中表示为通过六个特征状态的闭合路径。循环阶段及其对应生物过程包括：1–2等压冷却至T_min，V₁腔体体积减小（叶片因对流和蒸腾冷却）；2–3环境空气通过气孔通道进入，压力升至大气压（吸入含CO?空气）；3–4多变压缩伴随温度升高（叶片组织吸收太阳辐射而升温）；4–5等压加热至T_max，压力超过大气压（叶片组织因吸收太阳辐射而变暖）；5–6气体通过气孔排入环境（排出O?和水蒸气）；6–1冷却并返回初始状态。循环过程中叶片气孔腔体和气孔中的气体完成机械功L，其大小与循环中吸收的热量相对应，由循环1–2–3–4–5–6–1的面积量化。

### 3.1 生物系统中热效率的理论基础

为分析热转化为机械能的可能性，基于工程热力学中广泛使用的理想化方法对过程建模。基于理想卡诺循环，热机（生物原型）的热效率η_t取决于叶片组织温度及其变化：η_t = 1 ? (T_min/T_max)，其中T_min/T_max为一个循环中叶片组织的最低和最高温度。在极小的温度差（如1–3 °C）下，热效率系数η_t保持极低水平（不超过约1%），这与低势热转化的原理一致。η_t对温度绝对值依赖性弱，但对循环中的温度变化(ΔT)极其敏感，表明温度梯度的控制是叶片腔体产生热功的决定性因素。

### 3.2 气孔系统中理想化热力学循环的分析

与仅考虑极端温度的卡诺模型不同，所构建的理想化循环还考虑中间过程阶段，从而能更准确地确定实际可达的热效率。理想化循环由两条绝热线和两条等压线组成，假设工质为空气并具有理想气体性质。等压过程中植物叶片气孔发动机从环境获得的热量q₁ = c_p·(T₄ ? T₃)，向环境释放的热量q₂ = c_p·(T₁ ? T₂)。通过绝热过程的参数关系推导，最终热效率表达式为η_t = 1 ? (p₁/p₄)^(k?1)/k。温度变化与叶片和环境之间更高效的能量及气体交换相关。根据计算结果，ΔT = 0.5 °C时效率系数约为0.0025，ΔT = 10 °C时达到约0.035。即使温度差很小（ΔT = 1 °C），叶片发动机理论上也能产生激活CO?交换的机械力。效率主要取决于温度梯度而非绝对温度值。假设理论效率系数η_t = 0.003，考虑自然条件下落在叶片上的太阳辐射通量（20–200 W/m2），可有0.06至0.6 W/m2转化为机械能，约占叶片吸收太阳辐射通量的0.3%。尽管理论效率低，但表明即使极微小的温度差也能引发足以激活气体交换的机械效应。直接实验确认该循环面临三个障碍：研究对象是发生复杂CO?同化过程的活体器官，侵入式测量会干扰该过程；同化过程中CO?从物理状态转化为生化状态，涉及众多相互关联的生物和物理因素；测量对象本身是由微米级、毫米级和纳米级通道构成的气体填充结构，现有技术手段无法进行直接的压力和温度测量，因此需采用间接测量方法验证模型。

## 4. 植物叶片CO?气体吸收评估方法

气孔发动机解释了温度脉动如何产生决定叶片水平气体交换强度的机械力。为进一步将这一局部过程与整体结果（植物在整个冠层面积和活跃生长期内吸收的CO?量）联系起来，采用以下评估方法。植物在白天与环境进行能量和气体交换时，CO?气体的消耗、从环境中吸收以及在后续生化过程中的同化是代谢过程的重要方面，与植物发育和生长动态及碳封存密切相关。年度CO?气体吸收量计算公式为：M_CO?^y = A_tree · q_CO?，其中M_CO?^y为CO?吸收量(CO? eq. kg/年)，A_tree为树木冠层面积(m2)，q_CO?为CO?吸收量(CO? eq. kg/m2/年)。每日CO?气体吸收量计算公式为：M_CO?^d = M_CO?^y / τ_ABy，其中M_CO?^d为CO?吸收量(CO? eq. kg/天)，τ_ABy为树木每年进行CO?吸收的平均天数。

## 5. 植物叶片热机对CO?吸收和减少环境污染的意义

CO?封存量直接取决于活跃叶面积，生长快速且叶片茂盛的植物最为高效。光合作用在叶片中进行，吸收环境中的CO?并向环境释放O?作为代谢产物。植物和树木不会迅速降低环境中的CO?浓度，但树木具有突出的长寿特性，在较长的生命周期中可消耗大量温室气体CO?。平均叶面积的CO?封存约为每天0.002至0.004克，该值可能随具体环境条件变化。选择适应当地条件和气候的植物种类有助于提高CO?封存效率。在冠层水平上，CO?吸收由多个相互关联的因素决定：光通量决定叶片吸收的能量，进而决定叶片与环境温度之差；温度具有双重效应，较大的温差增加热力学分量，但过高的叶片温度抑制光合作用的生化能力；水分可利用性是限制因素，干旱时气孔关闭以保存水分，无论热循环如何，气体交换停止；植物种类决定气孔密度和开度，以及冠层密度和生长速率；光合活性决定叶片腔体与环境空气之间的CO?浓度梯度。城市树木对减少空气污染、降低环境温度有显著影响，可减缓烟雾形成过程并去除气体和颗粒污染物。具有气孔并进行能量和气体交换的植物可直接从环境中去除污染气体（CO、NO_X、O?和SO?）。研究表明，树木平均每年从环境空气中去除约0.002%（0.34 g/m2）的CO、0.8%（1.24 g/m2）的NO?、0.3%（1.09 g/m2·年）的SO?、0.3%（3.07 g/m2·年）的O?。冠层面积约50 m2的树木每年CO?吸收量为11至18 kg。若移除一棵健康的冠层面积为50 m2的成熟树木，环境将失去重要的CO?减排来源。按比吸收率0.26 kg CO? eq./m2/年计算，冠层面积50 m2的树木每年CO?消耗量为13 kg，而冠层面积为3 m2的幼树每年仅0.78 kg。树木相互遮荫会减少到达叶片的辐射通量，降低叶片与环境之间的温度差，从而降低气孔发动机的热效率，减少环境中CO?气体的消耗。

## 6. 讨论

本研究提出的基于热转化为机械功的叶片气孔结构热力学模型，允许通过工程学棱镜对叶片生理过程进行新的解释。虽然计算的热机效率(η_t ≈ 0.003)较低，但其作为支持气体运输的组分之一效果可能显著。文献强调植物叶片中即使微小的温度变化对蒸腾强度、气孔导度和光合速率也有显著影响。所提出的热力学模型与已建立的气孔调控解释不矛盾，可从三个方面关联：第一，决定气孔开度的保卫细胞膨压作为循环的边界条件，开度面积决定流动阻力；气孔完全关闭时，循环退化为无有用功的等容储热。第二，两种机制在不同时间尺度上运行：保卫细胞的渗透和膨压调节为分钟量级，而模型中的温度脉动为十分之几秒至几秒，因此对单一循环而言气孔导度可视为准恒定。第三，两者关系是相互的：叶片温度升高促进气孔开放作为蒸发冷却机制，较大的开度面积反过来降低流动阻力并增加循环的质量通量，形成正反馈回路。物种差异通过不同的叶片升温间接反映在模型中，唯一输入是叶片温度与环境温度之差。对温带树冠的热成像测量显示，郁闭冠层物种（如椴树或山毛榉）的叶片在晴天正午比气温高约4.5–4.8 °C，而冠层较疏或蒸腾更旺盛的物种（如橡树）约高2.4 °C，这对应热效率约两倍的差异。

### 6.1 对城市绿地和环境质量的实践意义

所提出的气孔发动机模型和CO?吸收评估对城市绿地规划具有直接的实际意义。城市树木通过蒸腾和遮荫降低局部空气温度，调节微气候并改善空气质量。该模型表明，相互遮荫的树木叶片与环境温度差较小，气孔发动机以较低热效率运行，因此在规划种植时，不仅总冠层面积重要，树木间距和冠层分层同样关键。城市中树木的位置还影响气流和污染物的扩散，特别是在街道峡谷中。冠层覆盖、污染与热岛减少之间的关系表明结构化绿化政策的重要性。

### 6.2 模型的局限性和进一步研究方向

该模型是理论性和理想化的，评估时需考虑若干局限性。第一，热力学循环基于理想化假设：工质为干燥理想气体，腔体体积和环境压力恒定，生物因素（蒸腾、光合动力学、细胞代谢）未直接包含在循环中；循环也视为封闭循环，尽管真实振荡系统中存在与环境的质量交换，该假设仅在单次温度脉动持续时间内有效。模型未直接参数化叶片相对于太阳的角度和冠层分层，这些因素通过到达叶片的辐射通量范围间接纳入。第二，模型结果——热效率和机械能通量——是理论输出而非直接测量值，活体叶片微米级、毫米级和纳米级通道中的直接测量无法用现有技术实现，因此必须采用间接测量方法验证模型。模型假设的各个要素——温度差、压力变化、腔体体积变化和通过气孔的质量流——已在文献中得到实验确认，但尚无研究在同一叶片上同时测量所有参数。同步验证需要四项测量：使用高分辨率红外热成像和微热电偶以至少10 Hz频率记录叶片组织温度脉动；使用压力探针技术间接评估叶片内部腔体的压力变化；使用微计算机断层扫描在测量温度的同时监测细胞间隙体积变化；高频气体交换测量以区分扩散和对流流动分量，其中分离对流分量是模型验证的关键标准，因为扩散传输机制已被公认。未来工作应建立气孔系统的物理仿生模拟器，并在受控条件下评估热-机械能转化，同时将理论循环模型与实验测量的气体交换参数和温度脉动数据相关联。

## 7. 结论

植物叶片气体交换系统的特性允许从理论上将其解释为按热力学循环原理运行的微尺度热机的生物原型。模型预测叶片与环境之间ΔT≈1–10 °C的温度差可产生有助于气体通过气孔系统通道运动的机械力。计算的热效率η_t ≈ 0.003是在理想化循环假设下获得的理论上限，表明这种能量转化发生在低势热区域。先前进行的叶片旋转实验提供了机制存在的间接证据，活体叶片的旋转角在相同条件下超过死叶片和非生物模拟近四倍。直接确认需要同步测量活体叶片中的温度脉动、腔体压力和体积变化以及通过气孔的气体流量，特别是分离对流流动模式与扩散模式。所提出的模型为仿生学和微能量学研究方向提供了可验证的假设，即生命系统中的被动能量转化。CO?吸收估算基于文献数据通过理论计算获得，在应用于城市绿地规划或微气候管理之前，必须在冠层尺度上通过测量进行确认。