篷歉19年9 97月.pdfVIP

1997年9月 系统工程理论与实践 第g期 (18I吗／数据处理之进展一一见于气象四维同化问题 。—、＼ lIV’]I (中国科学院系统科学所 ·北京 100o8d) 0 摘要 对迅速发展的四堆同化方法作了评价。尤其是A．c．Loreoc的论文与方案t对王f保、地庳…散 据处理，有参考价值。 关键词 数据处理 四维同化 AIC准则 数学建檬 DataProcessinginMeteorlolgicalDataAssimilation LuoQiaolin (AcademiaSinica，InstituteotSystemsScience，13eijing100080) Keywords dataprocessing；dataassmilation AICprincplelmathematicalmodeling 1 导言 墨堡旦些旦萼t是塞熬擅面益崖二l：童，即微分方程组求解中它担任给出较合理初的譬 务+简覃说，全国各气象站有气压、温度……的观测值，但微分方程数据值预报，需要在划一的网格点上，给 出初值a 地质探矿中+由少数探井的情况；推断全矿区的情况，卫星定轨中，由雷达观测值，确定卫星的真实轨 道……均属同一应用数学问题。 在 1960年前后+一般采取极简单的方法，如由网格点为中心，以适当半径作一圊，圊内各观测站的取 值的加权平均，作为该网格点的取疽；又如将全国分为若干区，各区内用多项式作回归，再适当考虑边界的 连接。 本文是在与大气所纪立人同志，气象中心朱宗申同志多年合作后的小结。这里，感谢他们及其他有关 同志t对我的指导与帮助。 本文中将介绍一些数据处理的方法，也将介绍在这三十多年中，各种方法兴衰荣辱中，毖藏在这历史 后的铁的法则——数据处理的思想，很多思想出自A—C．Lorenc的工作 介绍一个同化方法，它在应用上没有价值，但在理论上很富启发性-即各网格点上取值由与它最近的 观测点之取值决定 一 四维同化问题．实质上是一个回归问题，即全国各点上气压之取值，是客观存在的一个曲面(回归面)， 我们不知道，但它在各气象站上取值，我们知道，叫观测值，我们的任务是由观测值(集)，去确定出这整十 曲面的情况、构造出这个曲面。 请注意这可是一个非线性非参数回归。 不添加假定，这个问题是没有解答的，请想想这对不对?因为数学上有各种奇形怪状的曲面，一般要假 定它是光滑曲面(或分为有限块，每块是光滑曲面)，这个假定在物理上是合理的，也同实际观测相符合。 因而以后叙述的同化方法，实际上都规定了一个假定t即回归面在怎样一族曲面中选取不同的曲面的 本文于1995年l2月6日收到 第9期 数据处理之进展——见于气象四维同化问题 99 规定。不同的选优准则，导至不同的解。四维同化问题的发展，实删上是这些假定与寻优准则的变化与发 展，有的消失了．有些发展了，新的必起了，有的淘汰了。 2]中把四维同化问题描述为利用一切有教信息，去确定出大气状态在给定时刻的取值，这个定义是 耐^寻味的，请想想与本文最初介绍的定义矛盾不矛盾?有何异同。 2 最优插值与叠代寻优 在6o年代初期，各种统计方法被引进四维同化问题，运渐取代了早期的直观方法，成为主流t本文介 绍其中较突出的两种。 朋维同化问题是由观测值去确定网格点上的取值(一般叫分析值)，即由观测矢量Y’一(t…， )去 计算分析矢量 ’一L ，zz，… )：气象观测一般是多要素的，为菏鹏起见，先讲述单要素的。_v一般比M 要大一个数量级t单靠观测信息去确定网格点上取值，是很难准确的，需靠其它信息或知识来帮助，最方便 最可靠的是由数值预{匪所还来的信息。数值预报是基于流体力学与热力学之上一并综合了^们的经验与上 一 时刻的观测信息-这才是最优插值高于前述菏单直观方法的要点。最优插值实质上就是线性回归t但下 面的公式推