- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
2.背景场(background) 在对模式进行同化以前,对模式状态的一种初估状态.其中用到了一定的物理约束和模式状态在时间上连续的假设.可以是气候态或是其他的不重要的状态,也可以由先前的分析场产生. 二.资料同化的基本概念 3.分析场(analysis) 在一个给定的时间上,能够准确反映真实状态的结果.一般是对模式进行同化以后得到的状态. 二.资料同化的基本概念 模式 初始条件 (初值场) 预报结果 数 值 (天气,气候)预报 资料同化问题 数值天气预报实际上是求解控制大气运动的偏微分方程组的初边值问题。因此,利用数值模式做天气预报时,预报误差来源于初边值的误差和模式误差两个方面。资料同化的目的之一就是提高初值的质量,进而改进预报效果 。 资 料 同 化 数值模式+观测资料 = 更准确的初值 资料同化是指充分利用观测资料(常规、非常规),结合预报模式,尽可能准确地给出大气或海洋运动真实状态的估计。 1) 简单分析 2) 统计插值或最优插值 变分分析 资料同化发展的三个主要阶段 简单分析主要在20世纪50年代使用,那时计算机技术还相当落后,计算能力仍然非常低下。简单分析又主要分为主观分析和客观分析两种,它们可以说是资料同化最早的基础。当时,采用主观分析技术的最小化过程通常只起着“平滑”的作用,其结果非常不可靠,因为受人为因素影响较大。由于这个原因,当时客观分析更受到重视。 简 单 分 析 在20世纪60年代至70年代,统计学方法被引入到气象资料同化中来。由此,一些统计插值(或最优插值)方法被用来把观测资料同化到预报模式中去。后来,这些最优插值分析方法得到了广泛使用。 统计插值或最优插值 到了20世纪80年代至90年代,人们的注意力已逐渐转移到变分方法上来,尤其是四维伴随处理。这一技术不仅在气象和海洋观测资料的同化中得到广泛应用,而且也在敏感性分析和参数估计中发挥了重要作用。 变 分 分 析 资料同化的历史和现状 最优插值(OI)(60-80年代) 三维变分(3D-VAR) ECMWF(1996.1)、NCEP(原NMC) (1991.6) 四维变分(4D-VAR) ECMWF(1997.11) 区域模式 北欧八国(HIRLAM,1999, 3D-VAR) 资料同化目的 获得给定时刻的大气或海洋等的“真实”状态的分析值。 为天气与气候数值模式提供初值 为科学研究提供高质量的资料 环境监测等 大气、海洋、水文等的数值模拟和预报 河口发育的变分同化 1.3资料同化方法 资料同化的有效信息 观测值,背景场,两者各自的误差; 大气发展、变化遵循的物理、动力的规律。 如何进行资料同化? 现行的同化方法可分为两类:序列同化和变分同化。序列同化主要有卡尔曼滤波,集合卡尔曼滤波等 归结为最优化问题。 资料同化的基本原理 资料同化的基本原理 变分同化方法原理简介 代价函数(或罚函数)J : J(x) = ( x - xb )T B-1 ( x - xb ) + (R(x)-Yobs)T O-1 (R(x)-Yobs) 1 2 1 2 1 2 the background or the first guess N-component vector of analysis variables N forecast error covariance matrix N M-component vector of observations observation operator M observational error covariance matrix M 怎样变分同化? x0k+1 = x0k - Sk J(x0k)T { k where Sk a symmetric matrix, an approximation to the inverse Hessian. It only depends on J and . is the step size. J k Limited-memory quasi-Newton method: x01 = xb The minimization: To find an optimal estimate of the i
文档评论(0)