天气诊断分析报告

第二章资料处理和客观分析

§ 2.1资料处理

随着气象观测手段的发展和现代化，得到的资料数量增多，门类拓广，有常规站点观测，有非定点海洋观测，有飞机观测，有气象卫星观测，有定时观测和非定常观测等等，怎样使用这些资料呢？

一般来说，对气象资料的要求有两方面：一是可靠性，二是便于使人所周知。气象测站的分布是不规则的，因此我们只能得到这些不规则点上的气象资料，但是数值预报中的网格点是规则的，因而资料无法直接使用。另外，无论是用穿孔纸带或用电信号的形式将气象电报直接输入电子计算机，都要首先按照专门的程序

进行译码、检查、整理。因为气象电报的内容，是按照气象电码格式编发的，而它的形式又是按照邮电电码格式编发的，从观测、编码、发报，到传递、转换、接收等，在每个工序和环节上，都存在着出错的可能性。因此，我们所接收到的气象电码，不可避免地存在着一些错误或不妥之处。所以，要正确使用这些气象资料，必须经过必要的处理。比较简单的资料处理可分为以下几个方面：

2.1.1记录错情判断

在利用接收到的气象资料之前，首先对资料要作错情判断。一般的做法是根据不同等压面上各种要素值的大小，给出相应略大于其最大值的一个数作为其上限值；也给出相应的略小于其最小值的一个数作为其下限值。例如在我国范围内的各测站，冬季500百帕层上的位势高度最小值不超过500位势什米，我们则取500作为其下限值；最大值不超过600位势什米，我们就取600作为其上限值。然后利用比较大小的子程序，由计算机对每一组数据进行判断，凡是大于上限值或小于下限值的记录，我们就认为它是错误的，予以舍掉，作为缺测记录。另外还可利用气象要素在时间变化上的连续性和空间分布上的连续性，来判断一个气象要素记录是否错误。对于错误记录，可用下面介绍的补缺测记录的方法，另外补一个值。

2.1.2补缺测或漏传记录

一般可把缺测或漏传记录的测点看作是一个网格点，然后由下节介绍的客观分析方法，禾U用周围已有的测站记录，插一个值补上。

2.1.3实测风矢量的分解

气象台站观测到的风场资料，是一个既有大小又有方向的风矢量，为便于该资料的

利用经常将实测风分解为东西和南北两个分量。分别用u、v表示，并规定：u向东为正,

V向北为正。其数量值分别由下式计算：

u Vsin( /180 )

v V cos( /180 )

这里V为实测风速值，为测风报告发布的风向度数。如，实测风为西南风，为240 °，风速10m/s，算得：

u 10 sin(240 /180 ) 8.7(m/s)

v 10 cos(240 /180 ) 5(m/s)

2.1.4风场订正

诊断分析一般都是在有限区域内进行的，多数都采用正方形网格。客观分析后所得到的网格点上的u、v分量值，并不处处与网格区的x,y轴平行，因此还必须进行风向订正。因为只有在基线上的网格点，其东西、南北方向与x,y方向一致，其他网格点上东西、南北方向与x,y方向总有一个偏差角，这显然会给计算带来误差，特别是当计算范围取得较大时，边缘的网格上，这种风向的误差显得更加突出。不进行适当订正是不行的。

如图2.2，先考虑在基线以西的某一网格点A, N是北极，NO为基线，NA和CA分别为经过A 点的经线和纬线，MA和LA分别和这个正方形网格系统的X轴及丫轴平行。假设A点的风速在经纬的分量分别为u和V。而在网格的X,Y方向的分量分别为u '、v ', 由于A点所在的经线不与基线相平行，故u、v分量和u '、v'分量彼此也不平行，而是有一夹角a,a=Z ANO由图可看出，它们之间有如下的换算关系：

u' ucos vcos(90 ) ucos vsin

v' vcos ucos(90 ) vcos usin

(2.1)

(2.2)

图2.1实测风矢量的分解

在基线以西的网格 AO>0所以=sin -1(AO/AN)>0;在基线以东的网格，AO<0所以<0; 在基线上的网格，a= 0,则u ' = u 、v ' = v 。即不必订。

设在某网格点上，u = v= 10米/秒，a= 45°由(2.2)式订正后

u 10 cos45 10 sin45 14(米 /秒) v 10 cos45 10 sin45 0(米 / 秒) 2.1.5平滑和滤波

气象观测资料，总存在着各种各样的误差。比如由气象仪器安装不标准等带来的，

非偶然性误差(器差)和由工作人员在观测、编码、发收报等造成的偶然性误差，以及 将要素值内插到网格点上时，产生的舍入和插值误差等等。无疑，这些误差都将会影响 计算的结果，为了减少误差的影响，通常在计算之前先对原始资料等进行平滑和过滤

(滤

波)，滤掉那些次要的小的天气意义的东西，而保留和突出主要的量，或者，为了研究 的需要滤去资料中某些波长的量，而保留与问题有关的量。现分别介绍如下：

(1) 一维平滑算子

这是最简单的平滑算子。利用同一直线上三点的资料，又称三点平滑算子

S

S

f j

f j 2 f j 1 f j 1 2f j (1 S)f j 2 f j1 f j 1

(2.3)

这里f j 表示第j 点平滑后的值，f j 表示第j 点平滑前的值。S 为平滑系数(可正可负)。 该平滑算子对j 点对称，其权重除j+1，j 和j-1点外均为0。

N

图2.2风场订正示意图

C

对函数f (X )可展成富氏级数，在f j 点，可写成:

ik(x.)

f j (x) C Ae j

比较(2.4)、(2.6)式，平滑后的波相未变，改变的只是波的振幅，平滑后的振幅为:

令R = A/A ，称响应函数，即平滑波幅同原波幅之比，表示平滑后的变化。显然， R = 1,表示平滑后波幅一样，R<1,表示平滑使原波幅衰减。R -0,表示平滑后使原波 动被消失(即波动全被滤掉)。

R>1,表示平滑后使原波幅被放大。由(2.7)式知：

R(k,s) A/ A 1

S(1 cosk x)

或者 R(L,S) 1 2S?sin 2( x/ L)

(2.8)

可见对于固定的网格距，响应函数 R 只与波数k (波长L )以及平滑系数S 有关。

由于0W sin 2( △ x/L) < 1若希望平滑后使原波动衰减，以致滤掉(但不希望出现 反位相情况)，由只须要有0W R < 1,于是由(2.6)式知：0W S W 1/2。倘若希望平滑后使 原波动增幅，即R>1,则必有平滑系数S<0b

作为一特例，取S=1/2，即最大平滑系数，此时响应函数为

(2.4)

这里C 为常数，A 为波动的振幅。k 2 /L 为波数，L 为波长,

为位相。同样在X j !

f j 1(x)

C ik (x ； x

Ae j

)C ik x ik(X j Ae e

)

f j 1(x) C "ik (X ；

x

) —

" ik x ik(x

i )

Ae j

C Ae e

将f j1, f j ,

f j1 代入(2.3) 式得:

f j (x) C (1 S)Ae

ik(Xj

)

2

)ik x

(e

C A1 S(1

cosk 、 ik(X

j

)

x) e

(2.5)

A A1 S(1 cosk x)

(2.7)

图2.3

和X j !点，函数f (X)可写成

e ik x )

(2.6)