专利名称：一种中高层大气风场反演方法
一种中高层大气风场反演方法技术领域
本发明属于激光雷达探测技术领域，尤其涉及一种中高层大气风场反演方法。
技术背景
中高层大气风场参数是空间物理学研究的重要参数之一，结合空间物理的理论探索，能够揭示空间的物理现象以及内含的相互作用和因果关系，同时中高层大气风场参数在气象研究、天气预报、大气环境监测和国防高技术等方面都具有广泛的应用。
激光多普勒遥感测量已经有几十年的发展过程，所采用的多普勒测量技术包括相干(外差)技术和非相干技术。其中，非相干技术对光学系统的要求不高、容易实现和加工， 其应用范围还扩展到分子散射的探测。非相干系统利用大气粒子产生的多普勒频移直接得到径向速度分布，是目前测量大气风场的最有利的工具，并且非相干系统具有测量精度高、 分辨率高和三维风场信息等特点。目前国际普遍采用直接接收(非相干)激光雷达作为大气风场测量的手段。
多普勒测风激光雷达径向单波束探测得到的是水平风场在该探测方向的径向分量，在适当的扫描方式下，多波束组合即可反演得到三维风场。目前，多普勒测风激光雷达采用具有固定仰角和相同间隔方位角的四个波束的扫描方式，即采用四波束进行风场反演，原理如图1所示。
以正东方向为X轴正方向，正北方向为Y轴正方向，天顶的方向为Z轴正方向建立坐标系，激光的发射仰角为θ，分别探测四个径向的风场，从而反演东西南北四个方向的径向风速，径向风速是指实际风速在激光束指向的这个方向上的分量值。一般规定正北方向为0°，东、南、西按顺时针依次为90°、180°和360°，获得径向风速VH。
Vn = Vdl -i = N, S,E,ff (1)
式(1)中Vdi为多普勒频移，λ为激光波长，根据几何关系可得
VrN =VyCosO+ V2SmO VrE =Vx cos θ+ V2sm θ Vrs=-Vy cos θ+ V2sm θ Vrw =-Fx cos ^+ Fz sin ^⑵
经过推导可得大气风场在χ轴、y轴、ζ轴的各自分量为
ν -VV -V V +V +V +Vy _ v rE v rW γ _ ν rN ν rS γ __ rE ^ r rW ^ r rN ^ r rS(3)2 cos 6 r 2cos6>4 sin 6
则水平风速大小Vh和方向γ分别为2cos6> 丨
y = arctan (Vx/Vy) + π {1-sign [ (Vy+1 Vy |) ‘ VjjVy ^ 0
在传统的大气风场反演方法中，需要通过多普勒测风激光雷达反演大气层中四个
Vh =J(Vx)2HVy)2-Kw)2+ (K,-VJ2T ⑷4方向的径向风速，在反演径向风速的过程中，首先向大气层某一预设方向发射激光信号，激光信号遇到大气分子产生瑞利后向散射，瑞利后向散射大气回波信号被接收，之后根据该方向的瑞利后向散射大气回波信号确定该方向上的多普勒频移，根据公式(1)确定大气层中该方向的径向风速，在获取大气层中四个预设方向的径向风速后，根据公式C3)和(4)计算水平风速大小和方向。
在实施过程中，要借助于一台多普勒测风激光雷达依次探测大气层中四个预设方向的径向风速，这导致大气风场测量时间较长，降低了时间分辨率，时间分辨率是指在同一区域进行的相邻两次风场观测或预报结果的最小时间间隔。发明内容
有鉴于此，本发明的目的在于提供一种中高层大气风场反演方法，可以缩短大气风场测试时间，提高时间分辨率。
为实现上述目的，本发明提供如下技术方案
一种中高层大气风场反演方法，包括
确定三个探测方向，三个探测方向包括第一探测方向、第二探测方向和第三探测方向，所述第一探测方向的天顶角为、其方位角为Q1，所述第二探测方向的天顶角为 Φ2、其方位角为θ 2，所述第三探测方向的天顶角为0°，其中，所述Q1* θ 2相差90°，所述牝和φ2相同、大于0°且小于90° ；
分别反演所述第一探测方向上的径向风速Vri、第二探测方向上的径向风速Vrt和第三探测方向上的径向风速V,3 ；
在预先构建的三维坐标系中，根据所述Vri、Vrt和Vrt计算大气风场的水平风速Vh、 以及所述大气风场的方向与所述三维坐标系中X轴正方向之间的角度γ。
优选的，在上述方法中，所述预先构建的三维坐标系，以正东方向为X轴正方向、 以正北方向为Y轴正方向、以指向天顶方向为Z轴正方向，则
Vh =机2 + V；,
y = arctan(NjNy) + π {1-sign[ (Vy+1 Vy |) .VJ}，其中，^rVr2 cos φχ - F1 cos φ2sin φ2 sin θ2 cos φχ - sin φχ sin θχ cos φ2χ sin φ2 cos θ2 cos φχ - sin φχ cos θχ cos φ2 sin φ2 cos θ2 cos φχ - sin φχ cos θχ cos φ2
(Vr3 cosφι - Vrl)(sinφ2 cosθ2 cosφι - sinφι cosθλ cosφ2) + (Vr2 cosφι - Vrl cosφ2)(sinφι cosθλ) (sinφχ sin θχ)(sin φ2 cosθ2 cos病-sinφχ cosθχ cosφ2)~ (sinφ2 sin θ2 cos病-sin φχ sinθχ cosφ2)(sinφχ cosθχ)
(Vr3 cosφι - Vrl)(sinφ2 cosθ2 cosφι - sinφι cosO1 cosφ2) + (Vr2 cosφι - Vrl cosφ2)(sinφι cosO1)V =-少 (sinφχ sinO1 )(sinφ2 cosO1 cosφχ -sinφχ cosO1 cosφ2)- (sinφ2 sinO1 cosφχ -sinφχ sinO1 cosφ2)(sinφχ cosO1)ο
优选的，在上述方法中，反演任意一个探测方向的径向风速的过程，包括
向大气层中一个探测方向发射激光信号；
接收瑞利后向散射大气回波信号；
计算所述探测方向的多普勒频移vd，vd = [R-1 (v0+vd, Ta) -Tl-1 (v0)]，其中，R^1 (v0+vd, Ta)为频率响应函数R(V(1+Vd，Ta)的反函数，IV1 (Vtl)为IY (Vtl)的反函数，
R(v0+vd,Ta) = Tlil (v。+vd，Ta)/TK2 (v。+vd，Ta)，其中，Tiu(VfVdJa)为瑞利后向散射大气回波信号经过F-P标准具第一信号通道的透过率，TE2(v0+vd, Ta)为瑞利后向散射大气回波信号经过F-P标准具第二信号通道的透过率，
Tl(V0) = a5IJ(a4I J，其中，a5为多普勒测风激光雷达接收机中第五探测器的校准常数，I。为所述第五探测器探测到的光强，&为所述多普勒测风激光雷达接收机中第四探测器的校准常数，Ils为所述第四探测器探测到的光强；
根据所述探测方向的多普勒频移，计算所述探测方向的径向风速。
优选的，在上述方法中，计算所述测量方向的多普勒频移过程中，所述温度Ta为模式温度。
优选的，在上述方法中，计算所述测量方向的多普勒频移过程中，所述温度Ta为实测温度；
确定所述实测温度的过程包括
测量接收到的来自于所述第三探测方向的瑞利后向散射大气回波信号的光子数 N (ζ)；
根据所述光子数计算所述大气层的大气密度P (ζ) , Ρ(ζ) = Ν[Ζ^(Ζ2ο)，其中，%N(Z0)Z0为参考高度，P (ζ0)为所述大气层中参考高度处的密度，N(Ztl)为所述大气层中参考高度处的光子数；
根据所述大气密度计算所述大气层中被探测高度处的实际温度T(Z)，
Γ⑷=叩。+其中，T(Z0)为所述参考高度处P(Z) M(Z0) R ^ p(z)大气层的温度，M(Ztl)为所述大气层中参考高度处的摩尔质量，M(Z)为所述大气层中被探测高度处的摩尔质量，R为气体常数，g(z)为重力加速度。
优选的，在上述方法中，所述Ct1和Φ2为30°。
由此可见，本发明的有益效果为本发明公开的中高层大气风场反演方法中，在对大气风场进行反演过程中，仅通过多普勒测风激光雷达对三个测量方向进行多普勒频移测试，之后分别计算各个方向的径向风速，最终根据三个测试方向的径向风速确定大气风场的水平风速和方向，由于整个大气风场反演过程中只需要测量三个测量方向的多普勒频移，与现有技术中测量四个测量方向的多普勒频移相比，缩短了测量时间，从而提高了时间分辨率。
另外，本发明在反演径向风速的过程中，采用被探测大气层的实测温度来计算多普勒频移，可以进一步提高大气风场数据的准确性。
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为现有的四波束扫描示意图2为本发明公开的一种中高层大气风场反演方法的流程图3为本发明公开的三波束扫描示意图4为本发明公开的反演一个探测方向的径向风速的方法的流程图5为现有的多普勒测风激光雷达接收机的结构示意图6为本发明公开的一种三维坐标系中激光束的分解图7为本发明中确定被探测大气层实测温度的方法的流程图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
本发明公开了一种中高层大气风场反演方法，可以缩短大气风场测试时间，提高时间分辨率。其基本原理为在对大气风场进行反演过程中，仅通过多普勒测风激光雷达对三个测量方向进行多普勒频移测试，之后分别计算各个方向的径向风速，最终根据三个测试方向的径向风速确定大气风场的水平风速和方向，由于整个大气风场反演过程中只需要测量三个测量方向的多普勒频移，与现有技术中测量四个测量方向的多普勒频移相比，缩短了测量时间，从而提高了时间分辨率。
参见图2，图2为本发明公开的一种中高层大气风场反演方法的流程图。包括
步骤Sl 确定三个探测方向。
三个探测方向包括第一探测方向、第二探测方向和第三探测方向。其中，第一探测方向的天顶角为、其方位角为Q1，第二探测方向的天顶角为Φ2、其方位角为θ2，第三探测方向的天顶角小3为0°。并且，Q1* θ 2相差90°，(^和Φ2相同、大于0°且小于 90°。即第三探测方向为垂直向上，第一探测方向和第二探测方向为斜向上，并且第一探测方向在水平面上的投影与第二探测方向在水平面上的投影正交。该三波束的示意图请参见图3。
步骤S2 分别反演三个探测方向的径向风速。
利用多普勒测风激光雷达对三个探测方向进行径向风速反演的过程基本一致。其中，反演一个探测方向的径向风速的过程如图4所示，包括
步骤S21 向大气层中一个探测方向发射激光信号。
步骤S22 接收瑞利后向散射大气回波信号。
当激光信号遇到大气分子后会产生瑞利后向散射，多普勒测风激光雷达可以接收该瑞利后向散射大气回波信号。
步骤S23 计算该探测方向的多普勒频移vd。
vd = [R-1 (v0+vd, Ta) -Tl-1 (v0)]，其中，vQ为多普勒频移，vQ为激光发射频率， R-1 (v0+vd, Ta)为频率响应函数R(V(1+Vd，Ta)的反函数，IV1 (Vtl)为Tl(V0)的反函数。
频率响应函数R(v。+vd，Ta) = Tei (v0+vd, Ta)/TE2(v0+vd, Ta)，其中，Tei (v0+vd, Ta)为瑞利后向散射大气回波信号经过F-P标准具第一信号通道的透过率，TE2(v0+vd, Ta)为瑞利后向散射大气回波信号经过F-P标准具第二信号通道的透过率。F-P标准具两个信号通道的透过率可以表示为
Tri (v0 +vdJJ = HXv0 +vd) fL(v0 +vd) fRay(v0 +vd,Ta)
其中，fL(v0+vd)为激光发射谱线，h, (v0+vd)为F-P标准具透过率函数，fEay(v0+vd, Τ)为瑞利后向散射增宽谱线，Ta为温度，i = 1，2为两个信号通道，《 ，，表示卷积。
TL(v0) = Mw^a4IJ，其中，a5为多普勒测风激光接收机中第五探测器的校准常数，ILe为该第五探测器探测到的光强，a4为多普勒测风激光雷达的接收机中第四探测器的校准常数，Ils为该第四探测器探测到的光强。多普勒测风激光雷达接收机的结构如图5所示，主要包括第一感测器101、第二探测器102、第三探测器103、第四探测器104和第五探测器105，其中，第五探测器105用于检测参考光的光强Ile，第四探测器104用于检测参考光透过F-P标准具后的光强U。
步骤S24 根据该探测方向的多普勒频移计算径向风速。
某个探测方向的径向风速由该探测方向的多普勒频移确定，径向风速与多普勒频移的关系为厂 =Vift^i = 1，2，3，其中，i表示三个探测方向，λ为激光波长。即某个探测方向的径向风速为该探测方向的多普勒频移与激光波长的乘积的1/2。
对每个探测方向分别执行步骤S21至S24的步骤后，就可以确定三个探测方向上的径向风速。
步骤S3 在预先构建的三维坐标系中，根据三个探测方向的径向风速计算大气风场的水平风速、以及大气风场的方向与三维坐标系中X轴正方向之间的角度。
确定三个探测方向的径向风速后，将其以向量的形式分别在预先建立的三维坐标系中进行表示，其中，向量的大小表示径向风速的大小、向量的方向表示探测方向。之后，基于三角函数公式计算大气风场的水平风速、以及大气风场的方向与该三维坐标系中X轴正半轴之间的角度。实施中，还可以进一步根据X轴正半轴所代表的实际方向和角度Y确定大气风场的方向。
图6为本发明公开的一种三维坐标系中激光束的示意图，该三维坐标系中，以正东方向为X轴正方向、以正北方向为Y轴正方向、以指向天顶方向为Z轴正方向。
第一探测方向巧的单位矢量为g = sin病cos《J+sin病sin《-) + cos病·系。式中，θ 为第一探测方向Α的方位角，Ct1为第一探测方向巧的天顶角J为χ方向的单位矢量， 为 y方向的单位矢量，石为ζ方向的单位矢量。
以此类推，
第二探测方向r2的单位矢量为g = sin武cos^J+sin武sin代_) + cos武4。其中，θ 2为第二探测方向r2的方位角，Φ2为第二探测方向r2的天顶角。
第三探测方向r3的单位矢量为= Sin疼cos .7+sin疼sin +cost系。其中θ 3为第三探测方向r3的方位角，Φ 3为第三探测方向r3的天顶角。由于第三探测方向r3为垂直向上，因此，92和φ3均为0°。
假定水平风矢量为歹=Vx -1+ Vy-J+ V2-k ,则三个径向风速为6/8页Vrl = Vx sin 於 Cos^1 + Vy sin 病 SinOl +Vz cos ^
\ Vr2 = Vx sin φ2 cos02 + Vy sin ^2 sin ^2 + Vz cos 於2Vr3 = Vx smφ3 cos03 +Vy sin sin ^3 + Vz cos
可求得r/Vr2 cos φλ - F1 cos φ2sin φ2 sin θ2 cos φλ - sin φλ sin θλ cos φ2
Vx =——x sin φ2 cos O1 cos φλ - sin φλ cos Ox cos φ2 sin φ2 cos O1 cos φλ - sin φλ cos Ox cos φ2(F3 cosφλ - F1 )(sinφ2 cosO1 cosφλ - sinφλ cosOx cosφ2) + (F2 cosφλ - F1 cosφ2)(sinφλ cosOx) (sinφλ sin θλ)(sin φ2 cosO1 cosφλ - sinφλ cosOx cosφ2) - (sinφ2 sin θ2 cosφλ - sin φλ sinθλ cosφ2)(sinφλ cosOx)
(Vr3 cosφλ - Vrl)(sinφ2 cosO1 cosφλ - sinφλ cosθχ cosφ2) + (Vr2 cosφλ - Vrl cosφ2)(sinφλ cosO1)V —-少 (sinφχ sinOx)(sinφ2 cosO1 cosφχ -sinφχ cosθχ cosφ2)~ (sinφ2 sinO1 cosφχ -sinφχ sinOx cosφ2)(sinφχ cosθχ)
大气层风场的水平风速大小、以及风场方向与三维坐标系中X轴正半轴之间的夹角分别为
Vh=^V1x +V1y
y = arctan (NjN1) + Ji {1-sign [ (Vy+1 Vy |) · Vj} Vy ^ 0
当然，三维坐标系可以采用多种形式，而并不限于上述公开的形式。优选的，以指向天顶的方向为Z轴的正半轴，以正南、正北、正东和正西四个方向中相邻的两个方向作为 X轴的正半轴和Y轴的正半轴，此时，可以降低确定\和Vy的计算量。
本发明上述公开的中高层大气风场反演方法中，在对大气风场进行反演过程中， 仅通过多普勒测风激光雷达对三个测量方向进行多普勒频移测试，之后分别计算各个方向的径向风速，最终根据三个测试方向的径向风速确定大气风场的水平风速和方向，由于整个大气风场反演过程中只需要测量三个测量方向的多普勒频移，与现有技术中测量四个测量方向的多普勒频移相比，缩短了测量时间，从而提高了时间分辨率。
实施中，在计算某个探测方向的多普勒频移的过程中，温度Ta可以采用模式温度， 该模式温度由本领域技术人员惯用的标准温度，但是大气层的温度受到很多因素的影响， 因此模式温度与被探测大气层的实际温度可能存在一定差异。为了进一步提高大气风场数据的准确性，本发明中优选被探测大气层的实测温度作为计算多普勒频移的参量。
参见图7，图7为本发明中确定被探测大气层实测温度的方法的流程图。包括
步骤SlOl 测量接收到的来自于第三探测方向的瑞利后向散射大气回波信号的光子数N(Z)。
在多普勒测风激光雷达接收到来自于第三探测方法的瑞利后向散射大气回波信号之后，通过光子计数探测器测量该回波信号的光子数，记为N(Z)。
步骤S102 根据该光子数计算所述大气层的大气密度P (ζ)。
激光雷达方程为
N(z) = ^4^exp[-2[a(zyz](5)ζO
在式(5)中，C。为系统常数，P (ζ)为大气密度，σ (ζ)为大气消光系数，N(ζ)为9
在中高空，大气的消光系数很小，则式(5)可以简化为TV⑷= ^4^,取某高度~Z为参考高度，可得
CrlP(Zn)N(Z0) =(6)Z
联立式(5)和式(6)可求得大气密度廓线
彻=(V)N(Z0)Z0
在式(7)中，ζ为要探测的高度，Z0为设定的参考高度，P (Z0)为大气层中参考高度处的密度，N(Ztl)为从参考高度处接收的光子数，N(Z)为从要探测高度处接收的光子数。
步骤S103 根据大气密度计算大气层中被探测高度处的实际温度T(Z)。
理想气体方程为
pV = mRT/M = nRT (8)
在式(8)中，ρ为状态参量压强，V为体积，m为质量，M为摩尔质量，R为气体常量， T为绝对温度，η为物质的量。
ccP(z) ^ \ pgdz(9)ζ
式(9)表示高度ζ处的气压等于从该高度到大气上界的单位面积气柱所受的重力。
将式(7)代入式⑶和式(9)，可以求得大气温度廓线T(Z)
T(Z) = T(Z0)pi^ M(z)(10)P(Z) M(Z0) R ^ p(z)v )
其中，T(Ztl)为参考高度处大气层的温度，M(Z0)为大气层中参考高度处的摩尔质量，M(Z)为大气层中被探测高度处的摩尔质量，R为气体常数，g(z)为重力加速度。
执行上述步骤SlOl至S103后，可以确定大气层被探测高度处的实测温度，将该实测温度用于计算多普勒频移的过程，可以提高多普勒频移的精度，最终提高大气风场反演数据的准确性。
实施中，当第一探测方向的天顶角Φ1和第二探测方向的天顶角Φ2很小时，会增大多普勒测风激光雷达中望远镜的安装难度，当第一探测方向的天顶角小工和第二探测方向的天顶角Φ2很大时，对多普勒测风激光雷达中激光器的能量要求较高。因此，本发明中可以设置牝和Φ2位于30° 60°之间，一方面降低了望远镜的安装难度，另一方面也降低了对激光器能量的要求。
优选的，将Ct1和Φ2设置为30°。
另外，可以调整第一探测方向和第二探测方向，使得第一探测方向和第二探测方10向在水平面的投影为正南、正北、正东和正西中相邻的两个，此时，Vx和Vy的计算公式被简化，从而为了降低计算大气风场的水平风速和方向的运算量。
例如，在图3中，在三维坐标系以正南方向为X轴正方向、以正东方向为Y轴正方向、以指向天顶的方向为Z轴正方向，同时，第一探测方向的方位角Q1Sf、天顶角牝为 30°，第二探测方向的方位角92为90°、天顶角，第三探测方向的方位角93为 0°、天顶角(^为0°。Vx=IVrlSvr3
求得，K-Ks
本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。
对所公开的实施例的上述说明，使本领域专业技术人员能够实现或使用本发明。 对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的，本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下，在其它实施例中实现。因此，本发明将不会被限制于本文所示的这些实施例，而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。
权利要求
1.一种中高层大气风场反演方法，其特征在于，包括确定三个探测方向，所述三个探测方向包括第一探测方向、第二探测方向和第三探测方向，所述第一探测方向的天顶角为、其方位角为Q1，所述第二探测方向的天顶角为 Φ2、其方位角为θ 2，所述第三探测方向的天顶角为0°，其中，所述Q1* θ 2相差90°，所述牝和φ2相同、大于0°且小于90° ；分别反演所述第一探测方向上的径向风速Vri、第二探测方向上的径向风速Vrt和第三探测方向上的径向风速Vrt;在预先构建的三维坐标系中，根据所述Vri、Vrt和Vrt计算大气风场的水平风速Vh、以及所述大气风场的方向与所述三维坐标系中X轴正方向之间的角度Y。
2.根据权利要求1所述的方法，其特征在于，所述预先构建的三维坐标系，以正东方向为X轴正方向、以正北方向为Y轴正方向、以指向天顶方向为Z轴正方向，则K=^jK2+ K,Y = arctan (Vx/Vy)+ JI {1-sign [ (Vy+1 Vy |) .VJ}，其中，F2 cos φχ - F1 cos φ2sin φ2 sin θ2 cos φχ - sin φχ sin θχ cos φ2Vx =---1 sin φ2 cos θ2 cos φλ - sin φλ cos O1 cos φ2 sin φ2 cos θ2 cos φλ - sin φλ cos O1 cos φ2(Vr3 cosφι - Vrl)(sinφ2 cosθ2 cosφι - sinφι cosO1 cosφ2) + (Vr2 cosφι - Vrl cosφ2)(sinφι cosO1) (sin φι sin O1 )(sin φ2 cosθ2 cos病—sinφι cosO1 cosφ2)- (sinφ2 sin θ2 cos病—sin φι sinO1 cos φ2)(sinφι cosO1)(Vr3 cosφι - Vrl)(sinφ2 cosθ2 cosφι - sinφι cosO1 cosφ2) + (Vr2 cosφι - Vrl cosφ2)(sinφι cosO1)V =-少 (sinφχ sinO1 )(sinφ2 cosO1 cosφχ -sinφχ cosO1 cosφ2)- (sinφ2 sinO1 cosφχ -sinφχ sinO1 cosφ2)(sinφχ cosO1)ο
3.根据权利要求1或2所述的方法，其特征在于，反演任意一个探测方向的径向风速的过程，包括向大气层中一个探测方向发射激光信号； 接收瑞利后向散射大气回波信号；计算所述探测方向的多普勒频移vd，vd = [R-1 (v0+vd, Ta) -Tl-1 (v0)]，其中，IT1 (v0+vd, Ta) 为频率响应函数R(V(1+Vd，Ta)的反函数，IV1 (Vtl) % Tl(V0)的反函数，R(v0+vd,Ta) = Tlil (v。+vd，Ta)/TK2 (v。+vd，Ta),其中，Tiu(Vc^VdJa)为瑞利后向散射大气回波信号经过F-P标准具第一信号通道的透过率，TE2(v0+vd, Ta)为瑞利后向散射大气回波信号经过F-P标准具第二信号通道的透过率，TL(v0) =Al^ZOi4IJ,其中， 为多普勒测风激光雷达接收机中第五探测器的校准常数，ILe为所述第五探测器探测到的光强，a4为所述多普勒测风激光雷达接收机中第四探测器的校准常数，Ils为所述第四探测器探测到的光强；根据所述探测方向的多普勒频移，计算所述探测方向的径向风速。
4.根据权利要求3所述的方法，其特征在于，计算所述测量方向的多普勒频移过程中， 所述温度Ta为模式温度。
5.根据权利要求3所述的方法，其特征在于，计算所述测量方向的多普勒频移过程中， 所述温度Ta为实测温度；确定所述实测温度的过程包括测量接收到的来自于所述第三探测方向的瑞利后向散射大气回波信号的光子数根据所述光子数计算所述大气层的大气密度P (ζ)，P(z) = N(')Z fZ2o)，其中，Z0为参N(Z0)Z0考高度，P (z0)为所述大气层中参考高度处的密度，N(Ztl)为所述大气层中参考高度处的光子数；根据所述大气密度计算所述大气层中被探测高度处的实际温度T(Z)， Γ㈡=作。)孕# +华「。^ii/z,其中，T(Z0)为所述参考高度处大气M(Z0) R ^ ρ(ζ)层的温度，M(Ztl)为所述大气层中参考高度处的摩尔质量，M(Z)为所述大气层中被探测高度处的摩尔质量，R为气体常数，g(z)为重力加速度。
6.根据权利要求1所述的方法，其特征在于，所述Ct1和Φ2为30°。
全文摘要
本发明公开了一种中高层大气风场反演方法，包括确定三个探测方向，所述三个探测方向包括第一探测方向、第二探测方向和第三探测方向，所述第一探测方向的天顶角为φ1、其方位角为θ1，所述第二探测方向的天顶角为φ2、其方位角为θ2，所述第三探测方向的天顶角为0°，所述θ1和θ2相差90°，所述φ1和φ2相同、大于0°且小于90°；分别反演所述第一探测方向上的径向风速Vr1、第二探测方向上的径向风速Vr2和第三探测方向上的径向风速Vr3；在预先构建的三维坐标系中，根据所述Vr1、Vr2和Vr3计算大气风场的水平风速Vh、以及所述大气风场的方向与所述三维坐标系中X轴正方向之间的角度γ。本发明公开的中高层大气风场反演方法，缩短了测量时间，从而提高了时间分辨率。
文档编号G01S17/95GK102508222SQ201110369049
公开日2012年6月20日 申请日期2011年11月18日 优先权日2011年11月18日
发明者夏海云, 孙东松, 王国成, 窦贤康, 胡冬冬, 薛向辉 申请人:中国科学技术大学