本研究由俄罗斯托木斯克理工大学能源工程研究所的Nikolay V. Baranovskiy和Alexander E. Nee共同完成，其论文《Mathematical Simulation of Heat Transfer in Surface Soil Layer at Forest Fire Influence in Three-Dimension Statement》发表于2015年的《EPJ Web of Conferences》期刊第82卷。

研究背景 本研究属于野火生态效应与传热学交叉领域，具体聚焦于森林火灾对地表土壤层的热影响。森林火灾不仅直接烧毁植被，其产生的剧烈热量也会向下传递至土壤，深刻影响土壤的物理化学性质、水分状况以及其中的微生物群落（微生物生物地理群落， microbiogeocenosis）。以往的研究在评估这种热影响时，常采用简化的一维模型来模拟土壤中的热传导过程。然而，实际的火灾火源（ignition centre）往往是局部的、非均匀的，其热传递过程具有显著的空间（三维）特征。虽然已有一些实验研究试图量化这种影响，但野外或受控燃烧实验通常耗费巨大资源与时间，且存在较高的安全风险。相比之下，数值模拟方法提供了一种高效且安全的替代研究途径。先前的研究表明，即使是相对简单的数学模型，其数值结果也能与实验数据较好地吻合。因此，本研究旨在建立并求解一个三维数学模型，以更真实地模拟局部火源影响下土壤层中的非稳态热传导过程，从而更精确地评估森林火灾对表层土壤的热影响范围与程度。

详细研究流程 本研究是一项纯粹的数值模拟研究，不涉及实体实验对象，其核心流程围绕数学模型的建立、求解与结果分析展开，具体可分为以下几个步骤：

1. 物理与数学模型构建：这是研究的核心基础。研究者将地表系统抽象为一个四层结构，自上而下分别为：空气层、森林可燃物层、富含有机质的表层土壤（厚度5厘米）、黏土层（厚度5厘米）。将土壤视为两层系统是基于不同层次物质成分（有机质与黏土）的热物理性质差异。模型假设系统内的热量传递仅通过热传导机制进行，忽略了空气层中的气相过程。森林可燃物层中考虑了其热分解反应，作为一个内热源。整个计算区域如图1所示。

2. 控制方程与定解条件确立：针对上述四个区域，分别建立了描述其内部温度场变化的导热微分方程。

   • 对于黏土层、有机质土层和空气层，其控制方程为经典的非稳态导热方程。
   • 对于森林可燃物层，其控制方程在导热方程的基础上增加了一项源项，用以描述可燃物热分解反应释放的热量，其形式为 ( q_p k_p \rho_3 \xi \exp(-E_1/(RT_3)) )，其中包含了反应热效应 ( q_p )、指前因子 ( k_p )、活化能 ( E_1 ) 和干燥有机质体积分数 ( \xi ) 等动力学参数。
   • 此外，还为可燃物层建立了一个描述干燥有机质体积分数消耗的动力学方程。
   • 初始条件设定了各层在时间为零时的均匀温度（土壤温度、有机层温度、可燃物层温度、空气温度）以及可燃物中干燥有机质的初始体积分数。
   • 边界条件方面：在空气层外边界设定了第三类边界条件（对流换热）；在垂直于z轴的计算区域底部边界也设定了第三类边界条件，以考虑与更深处土壤的热交换；在各层之间的交界面上，设定了第四类边界条件（温度连续、热流密度连续），以确保热量传递的连续性。

3. 数值求解方法实施：由于所建立的偏微分方程组在三维空间上具有非线性特征，无法获得解析解，因此采用了数值方法进行求解。具体步骤如下：

   • 离散化：采用有限差分法（finite difference method）在均匀网格上对控制方程和定解条件进行离散。
   • 算法选择：对于导热方程的离散形式，采用了萨马尔斯基（Samarskiy）的局部一维格式进行处理，将多维问题转化为一系列一维问题。这些一维差分方程随后使用追赶法（marching method，通常指Thomas算法）进行求解。
   • 非线性处理：森林可燃物层方程中的非线性源项（与温度指数相关）通过简单迭代法（method of simple iteration）进行处理。
   • 参数设定：研究使用了具体的物性参数和反应动力学参数进行计算，例如各层的密度（ρ）、比热容（c）、热导率（λ），以及可燃物热分解的反应热、指前因子、活化能等。这些参数为模拟提供了具体的物理基础。

4. 模拟计算与结果分析：在设定的初始和边界条件下，运行数值计算程序，获取不同时刻（如t=60秒，300秒，600秒）整个三维计算区域内的温度场分布。研究者通过二维截面图（如x=11米和x=8米处的截面）可视化温度分布（图2，图3），并提取了特定方向（如z轴方向）上的温度数据进行定量分析。特别地，研究还将三维模型的模拟结果与对应的二维模型结果进行了对比分析，以评估模型维度简化带来的差异。

主要研究结果 数值模拟成功获得了局部火源影响下，“黏土-有机层-可燃物-空气”系统中三维温度场的时空演化数据。

1. 短期火源影响下的热穿透深度：模拟结果显示，在火灾影响的初期（例如60秒），火源产生的热量对土壤的影响范围较浅。从“土壤-可燃物”界面算起，显著的温升（温度变化明显）仅能到达约1至1.5厘米的深度。在此深度范围内，富含有机质的土壤层温度可升高超过100°C。如此高的温度可能导致土壤水分发生高温蒸发，并且在最顶部的厘米级土壤层中，可能发生炭化（charring）和植被残留物的热解（pyrolysis）。

2. 长期火源影响下的热穿透深度：当火源影响时间延长（例如300秒），热影响的深度显著增加。此时，显著的温升可以深入到“土壤-可燃物”界面以下4至5厘米处。在深度2至2.5厘米的有机质土层中，温度仍可超过100°C，导致更深处土壤水分的蒸发。同时，顶部约两厘米厚的土壤层将经历更剧烈的炭化和植被残留物热解过程。

3. 二维与三维模型对比结果：研究的一个重要发现是，在t=600秒时刻，沿z轴方向（深度方向）在表层土壤中选取一系列坐标点对比温度值（如表1所示），二维模型与三维模型的计算结果表现出极高的一致性，两者之间的差异不超过0.1 K。这表明，对于此类研究森林火灾对土壤热影响的问题，使用计算成本更低的二维模型进行模拟是可行的，其结果能够较好地近似三维情况。这一发现为开发用于地理监测系统（geomonitoring systems）的高效算法提供了依据。

研究结论与意义 本研究首次提出了一个用于评估森林火灾热模式对地表土壤层影响的简单三维数学模型。通过数值模拟，明确了在局部火源作用下，土壤温度的显著变化主要局限于富含有机质的表层。这种急剧的温度升高会对土壤微生物生物地理群落的功能和土壤本身的性质产生重大影响。 研究所获得的温度场分布数据，可以直接用于界定森林火灾火源对土壤生态功能（特别是微生物活动）产生影响的空间边界。这些成果可被整合进包含森林火灾预测模型的地理监测系统中，为火灾后的生态影响评估和管理决策提供定量化的科学参考。

研究亮点 1. 模型维度创新：首次针对森林火灾土壤热影响问题，建立了明确的三维非稳态热传导数学模型，相较于以往常用的一维模型，能更真实地反映局部火源的空间热传递特征。 2. 系统分层细化：将土壤-可燃物-空气系统精细化为四个具有不同热物性的层次，并考虑了可燃物层的热分解反应作为内热源，使模型更贴近实际物理场景。 3. 关键的对比验证：通过系统对比二维与三维模型的模拟结果，得出了“在此类问题中二维模型近似性良好”的重要结论。这一结论不仅验证了先前许多基于二维模型研究的合理性，更重要的是为后续需要快速计算的应用场景（如实时监测、风险评估）提供了使用简化模型的科学依据，具有重要的方法论意义和应用价值。 4. 明确的生态学指向：研究始终将物理模拟结果与潜在的生态学后果（如水分蒸发、土壤炭化、有机物热解、微生物群落影响）相联系，体现了跨学科研究的特色，使纯粹的物理数学模型输出具有明确的生态学解释和应用目标。

其他有价值的内容 研究在讨论部分提到了未来工作的可能性，例如可以使用更详细的数学模型来设定土壤中未受扰动时的垂直温度剖面作为初始条件，这暗示了模型进一步复杂化和精细化的方向。此外，文中引用的参考文献涵盖了从数学建模、土壤科学到火灾生态学的多个方面，为有兴趣的读者提供了深入探索该主题的完整学术脉络。