一种考虑减速伞气动力影响的风速风向探测方法

摘要

本发明公开了一种风速风向探测方法，将探空仪携带至预定高度及位置，精确投放后，探空仪在减速伞带动下自高空下落，实时测得高空至海平面的风速风向，测量过程中考虑减速伞气动力影响，对该影响产生的风速风向误差进行补偿，提高了精度。应用本发明方法探测范围精确可控，测风高度可根据需求设置，不仅可在陆地等空旷场地测风，还可至海上等从地面释放气球后其不易抵达的地方进行测风。

权利要求

1.一种风速风向探测方法，其特征在于，包括以下步骤：(1)将探空仪从预定高度向下投放，探空仪在减速伞的带动下下落飞行；

(2)获取当前时刻探空仪在x、y、z向的坐标分量xi、yi、zi、速度分量vxi、vyi、vzi以及气动力分量Fqxi、Fqyi、Fqzi，并将x、y、z向的气动力分量合成为气动力Fqi；

(3)计算当前时刻风速与探空仪的相对速度 其中，Cd为减速伞阻力系数，S为减速伞参考面积，ρi为当前时刻空气密度；

(4)结合气动力分量Fqxi、Fqyi、Fqzi计算风速与探空仪当前时刻的相对速度在x、y、z向的分量vxi相对、vyi相对、vzi相对：(5)计算当前时刻探空仪在重力方向上的平衡速度 进而得到该平衡速度在x、y、z向的速度分量

m为质

量，g为重力加速度；

(6)计算当前时刻x、y、z向的实际风速vxi实际、vyi实际、vzi实际：vxi实际＝vxi+vxi相对-vxi平衡vyi实际＝vyi+vyi相对-vyi平衡vzi实际＝vzi+vzi相对-vzi平衡；

(7)根据实际风速计算当前时刻在xy平面的风向θxy和xz平面的风向θxz：

2.根据权利要求1所述的风速风向探测方法，其特征在于，所述步骤(2)具体为：(21)探空仪采集得到前一时刻的坐标xi-1、yi-1、zi-1和当前时刻的坐标xi、yi、zi，同时根据相邻采集时刻的时间间隔Δt，分别计算x、y、z向当前时刻速度vxi、vyi、vzi，计算公式如下：探空仪采集得到当前时刻的坐标xi、yi、zi和下一时刻的坐标xi+1、yi+1、zi+1，根据相邻采集时刻的时间间隔Δt，计算下一时刻速度vxi+1、vyi+1、vzi+1，计算公式如下：(22)计算当前时刻x向加速度分量axi、y向加速度分量ayi及z向加速度分量azi，计算公式如下：(23)确定当前时刻x向合力分量Fxi、y向合力分量Fyi、z向合力分量Fzi，计算公式如下：Fxi＝maxi

Fyi＝mayi

Fzi＝mazi；

(24)根据探空仪重力mg及当前时刻的坐标xi、yi、zi，确定探空仪重力mg当前时刻在x向分量mgxi、y向分量mgyi及z向分量mgzi，计算公式如下：(25)根据当前时刻x方向合力Fxi，当前时刻y方向合力Fyi，当前时刻z方向合力Fzi，以及探空仪重力mg当前时刻在x向分量mgxi、y向分量mgyi、z向分量mgzi，计算气动力当前时刻的x向分量Fqxi、y向分量Fqyi、z向分量Fqzi，计算公式如下：Fqxi＝Fxi+mgxi

Fqyi＝Fyi+mgyi

Fqzi＝Fzi+mgzi；

(26)合成当前时刻气动力

说明书

一种考虑减速伞气动力影响的风速风向探测方法

技术领域

[0001] 本发明属于航天信息识别技术领域，具体为一种考虑减速伞气动力影响的风速风向探测方法。

背景技术

[0002] 高空定位测风通常是将飞升气球作为随气流移动的质点，用地面设备(经纬仪或雷达)跟踪气球的飞升轨迹，读取其时间间隔的仰角、方位角、斜距，确定其空间位置坐标，或气球携带探空仪，依靠卫星对探空仪定位；根据定位信息，求出气球所经过高度上的平均风速风向，测量高度可高达10公里以上。定位测风结果通常采用滤波和平滑的方法去除定位误差、气球的摆动和风的短周期变化造成的测量误差，获得大气平流的信息。

[0003] 在定位测风方法中通常依靠气球携带探空仪上升，从而根据其位置等信息进行风速、风向解析时，认为气球系统的质量较小，故忽略其惯性，认为气球在水平方向的运动就是水平气流的流动。实际上，飞行过程中气球依然受到气动力的影响，使得方法存在测量误差。另外，在气球带动下探空仪上升过程中气流作用稳定性差，且测量高度不可控。气球本身容易爆炸，导致探测失败。

发明内容

[0004] 本发明的目的是提供一种考虑减速伞气动力影响的风速风向探测方法，探测过程中考虑了减速伞气动力影响，对该影响产生的误差进行了补偿，提高了探测精度。

[0005] 实现本发明目的所采用的具体技术方案如下：

[0006] 一种风速风向探测方法，包括如下步骤：

[0007] (1)将探空仪从预定高度向下投放，探空仪在减速伞的带动下下落飞行；

[0008] (2)获取当前时刻探空仪在x、y、z向的坐标分量xi、yi、zi、速度分量vxi、vyi、vzi以及气动力分量Fqxi、Fqyi、Fqzi，并将x、y、z向的气动力分量合成为气动力Fqi；

[0009] (3)计算当前时刻风速与探空仪的相对速度 其中，Cd为减速伞阻力系数，S为减速伞参考面积，ρi为当前时刻空气密度；

[0010] (4)结合气动力分量Fqxi、Fqyi、Fqzi计算风速与探空仪当前时刻的相对速度在x、y、z向的分量vxi平衡、vyi平衡、vzi平衡：

[0011]

[0012]

[0013]

[0014] (5)计算当前时刻探空仪在重力方向上的平衡速度 进而得到该平衡速度在x、y、z向的速度分量
m为质
量，g为重力加速度；

[0015] (6)计算当前时刻x、y、z向的实际风速vxi实际、vyi实际、vzi实际：

[0016] vxi实际＝vxi+vxi相对-vxi平衡

[0017] vyi实际＝vyi+vyi相对-vyi平衡

[0018] vzi实际＝vzi+vzi相对-vzi平衡

[0019] (7)根据实际风速计算当前时刻在xy平面的风向θxy和xz平面的风向θxz：

[0020]

[0021]

[0022] 本发明相比于传统的风速风向探测方法具有如下优点：

[0023] 本发明探空仪由譬如火箭等探空设备携带至预定高度及位置，精确投放后，在减速伞带动下探空仪自高空下落，实时测得高空至海平面的风速风向，测量过程中考虑了减速伞气动力影响，对该影响产生的误差进行了补偿，提高了精度。探测范围精确可控，通过火箭等高空投放机构，可实现精确的定点定高投放探空仪，保证测试范围精确可控。该方法的测风高度可根据需求控制，同时不仅可在陆地等空旷场地测风，还可至海上等从地面释放气球后其不易抵达的地方进行测风。

附图说明

[0024] 图1是本发明减速伞携带探空仪飞行探空仪受力图；

[0025] 图2是本发明探测方法具体流程图。

具体实施方式

[0026] 为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。此外，下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。

[0027] 本发明高空风速风向探测方法，具体步骤如下：

[0028] (1)将探空仪从预定高度向下投放，探空仪在减速伞的带动下下落飞行。减速伞可采用球型伞、带条伞等等，优选球型伞，阻力小、姿态稳定性好。

[0029] (2)获取当前时刻探空仪在x、y、z向的坐标分量xi、yi、zi、速度分量vxi、vyi、vzi以及气动力分量Fqxi、Fqyi、Fqzi。

[0030] (21)探空仪采集得到地心空间(WGS-84)坐标系内前一时刻的坐标xi-1、yi-1、zi-1和当前时刻的坐标xi、yi、zi，同时根据获得两组坐标的时间间隔Δt，分别计算x、y、z方向当前时刻速度vxi、vyi、vzi，计算公式如下：

[0031]

[0032]

[0033]

[0034] 式中，vxi为当前时刻x方向速度，xi-1为前一时刻x方向位置坐标，xi为当前时刻x方向位置坐标，vyi为当前时刻y方向速度，yi-1为前一时刻y方向位置坐标，yi为当前时刻y方向位置坐标，vzi为当前时刻z方向速度，zi-1为前一时刻z方向位置坐标，zi为当前时刻z方向位置坐标，Δt为采集两组坐标的时间间隔。

[0035] 根据探空仪采集得到WGS-84坐标系内当前时刻的坐标xi、yi、zi和下一时刻的坐标xi+1、yi+1、zi+1，同时根据获得两组坐标的时间间隔Δt，计算WGS-84坐标系内下一时刻速度vxi+1、vyi+1、vzi+1，计算公式如下：

[0036]

[0037]

[0038]

[0039] 式中，vxi+1为下一时刻x方向速度，xi+1为下一时刻x方向位置坐标，vyi+1为下一时刻y方向速度，yi+1为下一时刻y方向位置坐标，vzi+1为下一时刻z方向速度，zi+1为下一时刻z方向位置坐标，Δt为采集两组坐标的时间间隔。

[0040] (22)根据获得的当前时刻速度vxi、vyi、vzi，下一时刻速度vxi+1、vyi+1、vzi+1，当前时刻的坐标xi、yi、zi及下一时刻的坐标xi+1、yi+1、zi+1，计算得到当前时刻x方向加速度axi、当前时刻y方向加速度ayi及当前时刻z方向加速度azi，计算公式如下：

[0041]

[0042]

[0043]

[0044] 式中，axi表示当前时刻x方向加速度，ayi表示当前时刻y方向加速度，azi表示当前时刻z方向加速度。

[0045] (23)根据计算得到的当前时刻加速度axi、ayi、azi，确定当前时刻x方向合力Fxi，当前时刻y方向合力Fyi，当前时刻z方向合力Fzi，计算公式如下：

[0046] Fxi＝maxi

[0047] Fyi＝mayi

[0048] Fzi＝mazi

[0049] 式中，Fxi为当前时刻x方向合力，Fyi为当前时刻y方向合力，Fzi为当前时刻z方向合力，m为探空仪质量。

[0050] (24)根据探空仪重力mg及当前时刻的坐标xi、yi、zi，确定探空仪重力mg当前时刻在x方向分量mgxi、当前时刻在y方向分量mgyi及当前时刻在z方向分量mgzi，计算公式如下：

[0051]

[0052]

[0053]

[0054] 式中，mgxi为探空仪重力当前时刻在x方向分量，mgyi为探空仪重力当前时刻在y方向分量及mgzi为探空仪重力当前时刻在z方向分量。

[0055] (25)根据当前时刻x方向合力Fxi，当前时刻y方向合力Fyi，当前时刻z方向合力Fzi，以及探空仪重力mg当前时刻在x方向分量mgxi、当前时刻在y方向分量mgyi、当前时刻在z方向分量mgzi，计算气动力当前时刻在x方向的分量Fqxi，气动力当前时刻在y方向的分量Fqyi，气动力当前时刻在z方向的分量Fqzi，计算公式如下：

[0056] Fqxi＝Fxi+mgxi

[0057] Fqyi＝Fyi+mgyi

[0058] Fqzi＝Fzi+mgzi

[0059] 式中，Fqxi表示气动力当前时刻在x方向的分量，Fqyi表示气动力当前时刻在y方向的分量，Fqzi表示气动力当前时刻在z方向的分量。

[0060] (26)根据气动力当前时刻在x方向的分量Fqxi，气动力当前时刻在y方向的分量Fqyi，气动力当前时刻在z方向的分量Fqzi，计算当前时刻气动力Fqi，计算公式如下：

[0061]

[0062] 式中，Fqi表示当前时刻气动力。

[0063] (3)根据当前时刻气动力Fqi及已知的减速伞阻力系数Cd，减速伞参考面积S，当前时刻大气密度ρi，计算风速与探空仪当前时刻的相对速度vi相对，计算公式如下：

[0064]

[0065] 式中，vi相对表示风速与探空仪的相对速度，Cd为探空仪减速伞阻力系数，S为探空仪减速伞参考面积，ρi为当前时刻空气密度。

[0066] (4)根据风速与探空仪当前时刻的相对速度vi相对，气动力当前时刻在x方向的分量Fqxi，气动力当前时刻在y方向的分量Fqyi，气动力当前时刻在z方向的分量Fqzi，当前时刻气动力Fqi，计算风速与探空仪的相对速度当前时刻在x方向的分量vxi相对，风速与探空仪的相对速度当前时刻在y方向的分量vyi相对，风速与探空仪的相对速度当前时刻在z方向的分量vzi相对，计算公式如下：

[0067]

[0068]

[0069]

[0070] 式中，vxi相对表示风速与探空仪的相对速度当前时刻在x方向的分量，vyi相对表示风速与探空仪的相对速度当前时刻在y方向的分量，vzi相对表示风速与探空仪的相对速度当前时刻在z方向的分量。

[0071] (5)根据探空仪重力mg、减速伞阻力系数Cd，减速伞参考面积S，当前时刻大气密度ρi，计算如果无风探空仪在重力方向上当前时刻的平衡速度vi平衡，计算公式如下：

[0072]

[0073] 式中，vi平衡表示如果无风探空仪在重力方向上当前时刻的平衡速度。

[0074] 根据当前时刻平衡速度vi平衡以及当前时刻位置坐标xi、yi、zi，计算当前时刻平衡速度在x方向的分量vxi平衡，当前时刻平衡速度在y方向的分量vyi平衡，当前时刻平衡速度在z方向的分量vzi平衡，计算公式如下：

[0075]

[0076]

[0077]

[0078] 式中，vxi平衡表示当前时刻平衡速度在x方向的分量，vyi平衡表示当前时刻平衡速度在y方向的分量，vzi平衡表示当前时刻平衡速度在z方向的分量。

[0079] (6)根据当前时刻速度vxi、vyi、vzi，当前时刻平衡速度在x方向的分量vxi平衡，当前时刻平衡速度在y方向的分量vyi平衡，当前时刻平衡速度在z方向的分量vzi平衡，风速与探空仪的相对速度当前时刻在x方向的分量vxi相对，风速与探空仪的相对速度当前时刻在y方向的分量vyi相对，风速与探空仪的相对速度当前时刻在z方向的分量vzi相对，计算当前时刻x方向实际风速vxi实际，当前时刻y方向实际风速vyi实际，当前时刻z方向实际风速vzi实际的计算公式如下：

[0080] vxi实际＝vxi+vxi相对-vxi平衡

[0081] vyi实际＝vyi+vyi相对-vyi平衡

[0082] vzi实际＝vzi+vzi相对-vzi平衡

[0083] 式中，vxi实际表示当前时刻x方向实际风速，vyi实际表示当前时刻y方向实际风速，vzi实际表示当前时刻z方向实际风速。

[0084] (7)根据当前时刻x方向实际风速vxi实际，当前时刻y方向实际风速vyi实际，当前时刻z方向实际风速vzi实际，计算xy平面内的当前时刻风向θxy和xz平面内的当前时刻风向θxz的计算公式如下：

[0085]

[0086]

[0087] 式中，θxy为xy平面内的当前时刻风向和θxz为xz平面内的当前时刻风向。

[0088] 循环进行步骤(1)～(7)即可实现在一段时间内考虑减速伞气动力影响的风速和风向的探测。

[0089] 本领域的技术人员容易理解，以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。