一个确定海-气界面气体交换率理论途径.pdfVIP

114 中国物理海洋学现状与展望 一个确定海一气界面气体交换率的理论途径 管玉平 (中国科学院南海海洋研究所热带海洋环境动力学重点实验室，广州510301) 大气CO。含量及其变化对气候的影响是当今全球变化研究中最重要的问 题之一。海洋对大气COz变化所起的作用(即源或汇的问题)是研究海洋在全 球气候变化中的作用的主要内容之一。要直接定量海域是大气COz的源或汇 的问题，需要高精度测定CO：在该界面的交换速率及其影响因素。而目前国际 上在测定海一气界面CO：的交换速率时采用的经验参数法不确定性因素很高， 需要在理论和测量方法上取得突破(戴民汉等，2001)。本文探讨理论上的一个 可能途径。 界面现象广泛存在于自然界中，然而，人们对于地球上最大的自然(气一液) 界面——海洋与大气的交界面——海一气界面上所发生的过程认识得还很不 够。实际海一气界面是很复杂的，海水流动、表面波动、风生波破碎和浪沫等都 影响海一气界面附近的结构，这给观测不同海表状况下确定海一气界面层带来困 难(高登义，1994)，同时给理论研究提出了许多挑战性课题。由于液一气(汽)界 面上即液体的表面性质和扩散系数等是化工、冶金、环保、石油等诸多工业领域 在设计、开发和生产过程中必不可少的基础物性参数，在实际工程中有着广泛 而重要的应用，尤其是化工领域对液一气(汽)界面处物质、热量的传输理论研究 早就备受关注。 研究表明，水一气界面附近湍流是气体传输的主要控制因素，但目前对于 水一气界面处气体传输机制还不能给出明确的结论性意见(贾复，2000)。所以， 海一气界面处发生了什么是我们感兴趣的，因为海一气间物质和能量的交换是在 海一气界面上进行的，换言之，海一气相互作用是首先在海一气界面上发生然后逐 渐展开的。海一气界面是液体的三维空间结构在与界面或表面垂直的方向上中 断或突然变化的体系，必然会形成一系列特殊的物理、化学性质而呈现出独特 的界面现象，如毛细现象、吸附、扩散与迁移、蒸发与凝结等等。虽然这些物质 和能量的传输过程很复杂，但最终都是通过分子进行的。 ’ 大尺度海洋环流与海气相互作用 115 由于气一液界面非常薄，使得现有的实验技术在界面现象研究的某些方面 存在一定的困难，例如难以对其分子微观结构和动力学行为进行研究，难以准 确测定紧邻界面处的热力学参数等。计算机分子模拟技术是经过几十年发展 的数值方法，它能够从分子水平上深入地研究界面现象，还能通过大量微观状 态的统计获得界面的宏观热力学性质，把微观与宏观联系起来，是研究界面现 象的有力工具，已成为物理、化学、生物、工程等许多学科领域的主要研究手段。 分子模拟研究方法的基本思想是为所研究的模型体系提供大量的微观状 态，在所有可能存在的微观状态中，体系中的每个分子都有一定的位置、取向和 动量。根据平衡态(或非平衡态)统计力学原理，体系的各种宏观性质如热力学 性质、传递性质、结构性质、光谱性质等都是相应微观量的系综平均值或时间平 均值。因此，只要在分子模拟过程中取得的微观状态足够多，就能够得到准确 可靠的模拟数据。 分子模拟有两种方法：一是MonteCarlo(MC)模拟方法，按照一定的概率 分布构作不同的微观状态，体系的宏观性质为大量微观状态的系综平均值；二 是分子动力学(MolecularDynamics，MD)模拟方法，所有的粒子先被赋予起始 位置和动量，通过求解体系中每个分子的经典牛顿运动方程，得到不同时刻体 系的微观状态，体系的宏观性质为大量微观状态的时间平均值。两者共同之处 在于都需要描述体系中分子之间相互作用的位能函数，区别是MC模拟方法通 常只用于平衡性质的模拟计算而MD模拟方法不仅可以模拟计算流体的热力 学平衡性质而且可模拟流体的传递性质。 从宏观物质系统是由大量微观粒子组成这一事实出发，统计物理学认为物 质的宏观性质是大量微观粒子运动的平均效果，宏观物理量是相应的微观物理 量的统计平均值。不难想像，海气相互作用所引起的现象最初是由海一气界面 上大量分子忙忙碌碌的结果。这似乎必须模拟十分庞大的粒子系统才能准确 再现真实系统的行为，但既使是用当今超级并行计算机来研究粒子数在1024量 级的体系仍不现实。那么，如何由少量粒子的模拟来得到宏观体系的性质呢? 对