大气中传输扩散计算中的数值方法（Ⅰ）：综述.pdfVIP

第五届全国风工程及工业空气动力学学术会论文集 大气中传输扩散计算中的数值方法((I):综述 桑建国 北京大学地球物理系，北京.100871 引 言 在空气污染气象学中大气污染物输送扩散计算中运用最普遍的公式是高斯(正态)烟云 模式。这类模式有明显的优点: 1.输入参数小而且规范化。模式的输入参数只有源强Q、有效源高H、平均风速U 和水平、垂直扩散参数ay和a}。其中源强为源的物理状态·与大气条件无关;有效源高 决定于大气条件和源的物理状态;完全决定于大气条件的参数只有三个，其中风速 U是 常数不随空间时间变化.而ay和。:只是距离、的函数，函数关系和大气条件、ijA.IR 有关，这类函数己经规范化，可在有关的手册和导则中查到。 2.输出结果和输入参数有明显的对应关系。例如浓度 C和源强的线性关系，和风速 u的反比关系，和有效源高的负指数关系等·Q7和6:的大小很直接地在浓度分布图中反 映出来. 3.便于平行比较。对于规一化的地面浓度CU/Q或 C/Q，按不同大气条件可做出标 准的分布图表，任何模式计算结果都可与这类标准图表对比，以判断参数选择和计算结 果的准确性。 a.计算便捷。高斯模式是一代数函数.其计算的快捷是其他模式无法比的。在需要 快速评价的问题中，例如核事故、有毒气体泄漏等事件中，可在极短时间内给出污染物 的分布.为应急响应争得时间。 但高斯模式的缺欠也正是因为它的简单。例如，大气中的风速 u从来不可能是一常 数。诸如复杂地形引起的中小尺度环流:高架源计算中遇到的风速风向的垂直切变;中 长距离传输中风速分布的水平不均匀和日变化等。 为适应复杂地形上的计算，根据湍流统计理论和计算流体力学方法己发展出各类复 杂模式。可统称为数值模式，以区别简单的高斯模式。但高斯模式仍然不可能完全被取 代，这是因为:1.有些模式过于复杂.其运行过程完全不透明，对于这类黑匣子的输出 结果.判断其可靠性的方法之一。是和高斯模式结果做平行比较，以避免出现太离谱的 结果;2.无论地形如何复杂，在适当短距离内，高斯模式总会在一定程序上近似成立:3.涉 及到计算长期平均浓度。还是要采用高斯型的模式。这里也包含着一种观念:复杂地形 引起的流场距其长期平均值的偏差，会在平均过程中得到消除。 数值模式有很多种类，其选取要根据计算区域的地形和大气条件而定。这些条件包 I一 21 { 括源的物理特性:譬如持续性:连续源还是瞬时源;几何形状:点源、面源或复合源: 数量:单源还是多源;高度:高架源还是地而源;计算域的水平范围以及地形的复杂程 度等等。 流场模式 流场模式提供出平均风场Vj(x,Y,z,1)，其中i=1,2,3分别代表u,v,w。有些模式还可 提供温度场和湍流场 (湍流动能)，根据湍流参数化方法.可得出湍流扩散系数K, 流场模式可分为两大类。 1.1 诊断模式 诊断模式或称客观分析，是根据风场观测资料.经过内插和调整，得出三维风场结 构。便如美国Livermore实验室的MATHEW模式。 1.1.1 资料来源。风场观测数据可取自于:i,常规气象观测，地面站和高空站。如果计算 范围在几百公里以上，在测站稠密地区，例如中国东部，可得到足够的资料运行此模式。 ii.大型试验中的布点观测。模式计算结果的好坏，取决于布点的密度及分布的合理性。 值得提出的是，高空风布点的分布密度永远是不足的，因此需要把地面风向上外推 (例 如幂指数外推)。在某些中尺度环流计算中，外推会带来严重问题，这将在后面讨论。 1.1.2调整过程。目前最常用的调整方法是采用变分法.目的是得出合理的垂直速度值. 例如陡峭的小尺度地形.w的取值范围为10-10Pms，对于较大较平坦的地形取值则要 相应减少。这个调整包含着半经验因素。 如果计算的范围超过几百公里，就没有必要再做调整。从而使模式计算简化。 1.1.3 模式优点。i.客观性。给定一组资料.给定调整参数，模式就会提供一个流场.这 个流场与运行者的经验、水平无关，因此称它为客观分析。相比之下，运行 同一动力模 式，因为运行者对大气运动的理解不同，可能得出不同结果.ii.真实性。计算域的大气流 场是多种因素构成的。是多种尺度大气运动系统共同作用的结果;这些系统的来源可能 是移动性的 (即外面移过来的).也可