土壤热传导方程解析解和那曲地区.PDFVIP

第 60 卷第 3 期 气 象 学 报 Vol. 60, No. 3 2002 年 6 月 ACTA METEOROLOGICA SINICA June 2002 土壤热传导方程解析解和那曲地区 土壤热扩散率研究X 高志球 卞林根 张雅斌 王金星 ( 中国气象科学研究院 北京, 100081) 姜冬梅 ( 山东省淄博气象局, 淄博 , 255048) 摘 要 文中用 Laplace 变换推导了土壤热传导方程的解析解和包含 热对流项的 土壤热传 导方程的解 析解。用青藏 高 原 8 个土壤湿度、温度廓线观测站 1998 年 9 月 4 日到 10 日实测资料基础上, 根据谐波 方法和 Laplace 变换方法得到 了土壤热传导 方程的解析解, 计算了这些站的总体土 壤热扩 散率; 用包含 热对流 项的土 壤热传导 方程的 解析解 计 算了土壤热扩散率。结果表明: 对于一 个深 度从 0. 04~ 0. 20 m 的浅 薄土壤 层, 总体 土壤热 扩散 率的 值为 0. 30 @ 10- 6 ~ 0. 98@ 10- 6 m2/ s, 土壤热扩散率的值为 0. 15@ 10- 6 ~ 0. 72@ 10- 6m2 / s。由谐波方法得到的 总体土壤热扩散率 比由 Laplace 变换的值稍大; 总体土壤热 扩散率总是比土壤热扩散率大。 关键词: 热传导方程, 解析解, 土壤温度廓线, 土壤热扩散率。 热通量和土壤表面的温度, 这些都迫切需要深入认 1 引 言 识土壤的热性质, 如热传导率、总体热扩散率以及体 青藏高原在亚洲季风 中扮演了重要角色[ 1, 2] , 积热容量等, 土 壤的体 积热 容量由 土壤 的成 分决 青藏高原的热动力作用是由高原表面地-气相互作 定[ 9] , 热传导率和总体热扩散率通过体积热容量联 用决定的。详细分析有关物理过程( 如辐射收支、感 系在一起, 所以仅需要确定其中一个。通常情况下, 热和潜热输送等) 、深入研究有关的物理参数( 空气 只研究总体热扩散率, 因为它描述热传导随温度边 动力粗糙度、土壤热传导率、土壤热扩散率等) 有利 界条件的瞬时过程。事实上, 土壤的热传输包含了 于更好地分析能量分配[ 3, 4] , 从而为利用一维、中尺 传导过程和多孔的对流过程[ 10] 。研究人员都将土 度或大尺度模式模拟地-气相互作用、中尺度, 大尺 壤热传导视为一个包含对流的总体过程, 常常见到 度天气过程提供这个特殊地区的参数化方案[ 5] 。为 确定总体土壤热扩散率或总体土壤热传导率的有关 了更深入地了解青藏高原的能量和水循环, 多