CN106290969A - 一种考虑减速伞气动力影响的风速风向探测方法 - Google Patents
一种考虑减速伞气动力影响的风速风向探测方法 Download PDFInfo
- Publication number
- CN106290969A CN106290969A CN201510242729.4A CN201510242729A CN106290969A CN 106290969 A CN106290969 A CN 106290969A CN 201510242729 A CN201510242729 A CN 201510242729A CN 106290969 A CN106290969 A CN 106290969A
- Authority
- CN
- China
- Prior art keywords
- current time
- component
- sonde
- wind speed
- speed
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Landscapes
- Indicating Or Recording The Presence, Absence, Or Direction Of Movement (AREA)
Abstract
本发明公开了一种风速风向探测方法,将探空仪携带至预定高度及位置,精确投放后,探空仪在减速伞带动下自高空下落,实时测得高空至海平面的风速风向,测量过程中考虑减速伞气动力影响,对该影响产生的风速风向误差进行补偿,提高了精度。应用本发明方法探测范围精确可控,测风高度可根据需求设置,不仅可在陆地等空旷场地测风,还可至海上等从地面释放气球后其不易抵达的地方进行测风。
Description
技术领域
本发明属于航天信息识别技术领域,具体为一种考虑减速伞气动力影响的风速风向探测方法。
背景技术
高空定位测风通常是将飞升气球作为随气流移动的质点,用地面设备(经纬仪或雷达)跟踪气球的飞升轨迹,读取其时间间隔的仰角、方位角、斜距,确定其空间位置坐标,或气球携带探空仪,依靠卫星对探空仪定位;根据定位信息,求出气球所经过高度上的平均风速风向,测量高度可高达10公里以上。定位测风结果通常采用滤波和平滑的方法去除定位误差、气球的摆动和风的短周期变化造成的测量误差,获得大气平流的信息。
在定位测风方法中通常依靠气球携带探空仪上升,从而根据其位置等信息进行风速、风向解析时,认为气球系统的质量较小,故忽略其惯性,认为气球在水平方向的运动就是水平气流的流动。实际上,飞行过程中气球依然受到气动力的影响,使得方法存在测量误差。另外,在气球带动下探空仪上升过程中气流作用稳定性差,且测量高度不可控。气球本身容易爆炸,导致探测失败。
发明内容
本发明的目的是提供一种考虑减速伞气动力影响的风速风向探测方法,探测过程中考虑了减速伞气动力影响,对该影响产生的误差进行了补偿,提高了探测精度。
实现本发明目的所采用的具体技术方案如下:
一种风速风向探测方法,包括如下步骤:
(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:
本发明相比于传统的风速风向探测方法具有如下优点:
本发明探空仪由譬如火箭等探空设备携带至预定高度及位置,精确投放后,在减速伞带动下探空仪自高空下落,实时测得高空至海平面的风速风向,测量过程中考虑了减速伞气动力影响,对该影响产生的误差进行了补偿,提高了精度。探测范围精确可控,通过火箭等高空投放机构,可实现精确的定点定高投放探空仪,保证测试范围精确可控。该方法的测风高度可根据需求控制,同时不仅可在陆地等空旷场地测风,还可至海上等从地面释放气球后其不易抵达的地方进行测风。
附图说明
图1是本发明减速伞携带探空仪飞行探空仪受力图;
图2是本发明探测方法具体流程图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。此外,下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
本发明高空风速风向探测方法,具体步骤如下:
(1)将探空仪从预定高度向下投放,探空仪在减速伞的带动下下落飞行。减速伞可采用球型伞、带条伞等等,优选球型伞,阻力小、姿态稳定性好。
(2)获取当前时刻探空仪在x、y、z向的坐标分量xi、yi、zi、速度分量vxi、vyi、vzi以及气动力分量Fqxi、Fqyi、Fqzi。
(21)探空仪采集得到地心空间(WGS-84)坐标系内前一时刻的坐标xi-1、yi-1、zi-1和当前时刻的坐标xi、yi、zi,同时根据获得两组坐标的时间间隔Δt,分别计算x、y、z方向当前时刻速度vxi、vyi、vzi,计算公式如下:
式中,vxi为当前时刻x方向速度,xi-1为前一时刻x方向位置坐标,xi为当前时刻x方向位置坐标,vyi为当前时刻y方向速度,yi-1为前一时刻y方向位置坐标,yi为当前时刻y方向位置坐标,vzi为当前时刻z方向速度,zi-1为前一时刻z方向位置坐标,zi为当前时刻z方向位置坐标,Δt为采集两组坐标的时间间隔。
根据探空仪采集得到WGS-84坐标系内当前时刻的坐标xi、yi、zi和下一时刻的坐标xi+1、yi+1、zi+1,同时根据获得两组坐标的时间间隔Δt,计算WGS-84坐标系内下一时刻速度vxi+1、vyi+1、vzi+1,计算公式如下:
式中,vxi+1为下一时刻x方向速度,xi+1为下一时刻x方向位置坐标,vyi+1为下一时刻y方向速度,yi+1为下一时刻y方向位置坐标,vzi+1为下一时刻z方向速度,zi+1为下一时刻z方向位置坐标,Δt为采集两组坐标的时间间隔。
(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,计算公式如下:
式中,axi表示当前时刻x方向加速度,ayi表示当前时刻y方向加速度,azi表示当前时刻z方向加速度。
(23)根据计算得到的当前时刻加速度axi、ayi、azi,确定当前时刻x方向合力Fxi,当前时刻y方向合力Fyi,当前时刻z方向合力Fzi,计算公式如下:
Fxi=maxi
Fyi=mayi
Fzi=mazi
式中,Fxi为当前时刻x方向合力,Fyi为当前时刻y方向合力,Fzi为当前时刻z方向合力,m为探空仪质量。
(24)根据探空仪重力mg及当前时刻的坐标xi、yi、zi,确定探空仪重力mg当前时刻在x方向分量mgxi、当前时刻在y方向分量mgyi及当前时刻在z方向分量mgzi,计算公式如下:
式中,mgxi为探空仪重力当前时刻在x方向分量,mgyi为探空仪重力当前时刻在y方向分量及mgzi为探空仪重力当前时刻在z方向分量。
(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
式中,Fqxi表示气动力当前时刻在x方向的分量,Fqyi表示气动力当前时刻在y方向的分量,Fqzi表示气动力当前时刻在z方向的分量。
(26)根据气动力当前时刻在x方向的分量Fqxi,气动力当前时刻在y方向的分量Fqyi,气动力当前时刻在z方向的分量Fqzi,计算当前时刻气动力Fqi,计算公式如下:
式中,Fqi表示当前时刻气动力。
(3)根据当前时刻气动力Fqi及已知的减速伞阻力系数Cd,减速伞参考面积S,当前时刻大气密度ρi,计算风速与探空仪当前时刻的相对速度vi相对,计算公式如下:
式中,vi相对表示风速与探空仪的相对速度,Cd为探空仪减速伞阻力系数,S为探空仪减速伞参考面积,ρi为当前时刻空气密度。
(4)根据风速与探空仪当前时刻的相对速度vi相对,气动力当前时刻在x方向的分量Fqxi,气动力当前时刻在y方向的分量Fqyi,气动力当前时刻在z方向的分量Fqzi,当前时刻气动力Fqi,计算风速与探空仪的相对速度当前时刻在x方向的分量vxi相对,风速与探空仪的相对速度当前时刻在y方向的分量vyi相对,风速与探空仪的相对速度当前时刻在z方向的分量vzi相对,计算公式如下:
式中,vxi相对表示风速与探空仪的相对速度当前时刻在x方向的分量,vyi相对表示风速与探空仪的相对速度当前时刻在y方向的分量,vzi相对表示风速与探空仪的相对速度当前时刻在z方向的分量。
(5)根据探空仪重力mg、减速伞阻力系数Cd,减速伞参考面积S,当前时刻大气密度ρi,计算如果无风探空仪在重力方向上当前时刻的平衡速度vi平衡,计算公式如下:
式中,vi平衡表示如果无风探空仪在重力方向上当前时刻的平衡速度。
根据当前时刻平衡速度vi平衡以及当前时刻位置坐标xi、yi、zi,计算当前时刻平衡速度在x方向的分量vxi平衡,当前时刻平衡速度在y方向的分量vyi平衡,当前时刻平衡速度在z方向的分量vzi平衡,计算公式如下:
式中,vxi平衡表示当前时刻平衡速度在x方向的分量,vyi平衡表示当前时刻平衡速度在y方向的分量,vzi平衡表示当前时刻平衡速度在z方向的分量。
(6)根据当前时刻速度vxi、vyi、vzi,当前时刻平衡速度在x方向的分量vxi平衡,当前时刻平衡速度在y方向的分量vyi平衡,当前时刻平衡速度在z方向的分量vzi平衡,风速与探空仪的相对速度当前时刻在x方向的分量vxi相对,风速与探空仪的相对速度当前时刻在y方向的分量vyi相对,风速与探空仪的相对速度当前时刻在z方向的分量vzi相对,计算当前时刻x方向实际风速vxi实际,当前时刻y方向实际风速vyi实际,当前时刻z方向实际风速vzi实际的计算公式如下:
vxi实际=vxi+vxi相对-vxi平衡
vyi实际=vyi+vyi相对-vyi平衡
vzi实际=vzi+vzi相对-vzi平衡
式中,vxi实际表示当前时刻x方向实际风速,vyi实际表示当前时刻y方向实际风速,vzi实际表示当前时刻z方向实际风速。
(7)根据当前时刻x方向实际风速vxi实际,当前时刻y方向实际风速vyi实际,当前时刻z方向实际风速vzi实际,计算xy平面内的当前时刻风向θxy和xz平面内的当前时刻风向θxz的计算公式如下:
式中,θxy为xy平面内的当前时刻风向和θxz为xz平面内的当前时刻风向。
循环进行步骤(1)~(7)即可实现在一段时间内考虑减速伞气动力影响的风速和风向的探测。
本领域的技术人员容易理解,以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (2)
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)合成当前时刻气动力
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510242729.4A CN106290969B (zh) | 2015-05-12 | 2015-05-12 | 一种考虑减速伞气动力影响的风速风向探测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510242729.4A CN106290969B (zh) | 2015-05-12 | 2015-05-12 | 一种考虑减速伞气动力影响的风速风向探测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106290969A true CN106290969A (zh) | 2017-01-04 |
CN106290969B CN106290969B (zh) | 2019-01-18 |
Family
ID=57630813
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510242729.4A Active CN106290969B (zh) | 2015-05-12 | 2015-05-12 | 一种考虑减速伞气动力影响的风速风向探测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106290969B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107132590A (zh) * | 2017-04-06 | 2017-09-05 | 中国科学院大气物理研究所 | 利用下投探空测量大气垂直风的方法和装置 |
CN109917493A (zh) * | 2019-04-03 | 2019-06-21 | 西安邮电大学 | 一种高空风能资源测量装置及方法 |
CN114675055A (zh) * | 2022-03-29 | 2022-06-28 | 北京航空航天大学 | 一种基于惯性系统的探空仪风速测量误差补偿方法 |
CN116306333A (zh) * | 2022-12-09 | 2023-06-23 | 中国能源建设集团广东省电力设计研究院有限公司 | 一种高空风能捕获装置空气动力学评估方法及系统 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE2939958A1 (de) * | 1978-10-02 | 1980-04-10 | Tracor | Einrichtung und verfahren zur wetteruntersuchung |
EP1327580A3 (en) * | 2002-01-15 | 2004-05-26 | Kawasaki Jukogyo Kabushiki Kaisha | Method and system for setting hull parameter of airship and method of adjusting ascension rate of the same |
CN101592741A (zh) * | 2009-07-03 | 2009-12-02 | 中国人民解放军理工大学气象学院 | 高空气象探测的北斗探测单元测风方法及其装置 |
CN102323627A (zh) * | 2010-10-29 | 2012-01-18 | 中国科学院大气物理研究所 | 一种多通道气象探空系统 |
-
2015
- 2015-05-12 CN CN201510242729.4A patent/CN106290969B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE2939958A1 (de) * | 1978-10-02 | 1980-04-10 | Tracor | Einrichtung und verfahren zur wetteruntersuchung |
EP1327580A3 (en) * | 2002-01-15 | 2004-05-26 | Kawasaki Jukogyo Kabushiki Kaisha | Method and system for setting hull parameter of airship and method of adjusting ascension rate of the same |
CN101592741A (zh) * | 2009-07-03 | 2009-12-02 | 中国人民解放军理工大学气象学院 | 高空气象探测的北斗探测单元测风方法及其装置 |
CN102323627A (zh) * | 2010-10-29 | 2012-01-18 | 中国科学院大气物理研究所 | 一种多通道气象探空系统 |
Non-Patent Citations (2)
Title |
---|
ESTELCARDELLACH, ET AL.: "Mediterranean Balloon Experiment: ocean wind speed sensing from the stratosphere, using GPS reflections", 《REMOTE SENSING OF ENVIRONMENT》 * |
张国舫 等: "GPS探空仪两种测风计算方法的比较", 《气象水文海洋仪器》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107132590A (zh) * | 2017-04-06 | 2017-09-05 | 中国科学院大气物理研究所 | 利用下投探空测量大气垂直风的方法和装置 |
CN107132590B (zh) * | 2017-04-06 | 2020-02-21 | 中国科学院大气物理研究所 | 利用下投探空测量大气垂直风的方法和装置 |
CN109917493A (zh) * | 2019-04-03 | 2019-06-21 | 西安邮电大学 | 一种高空风能资源测量装置及方法 |
CN114675055A (zh) * | 2022-03-29 | 2022-06-28 | 北京航空航天大学 | 一种基于惯性系统的探空仪风速测量误差补偿方法 |
CN116306333A (zh) * | 2022-12-09 | 2023-06-23 | 中国能源建设集团广东省电力设计研究院有限公司 | 一种高空风能捕获装置空气动力学评估方法及系统 |
CN116306333B (zh) * | 2022-12-09 | 2023-10-20 | 中国能源建设集团广东省电力设计研究院有限公司 | 一种高空风能捕获装置空气动力学评估方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN106290969B (zh) | 2019-01-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107966156B (zh) | 一种适用于运载火箭垂直回收段的制导律设计方法 | |
CN107844128B (zh) | 一种基于复合比例导引的高超声速飞行器巡航段制导方法 | |
CN103472503B (zh) | 探空仪及基于ins的高空风探测方法 | |
CN106290969A (zh) | 一种考虑减速伞气动力影响的风速风向探测方法 | |
CN102607639A (zh) | 基于bp神经网络的大攻角飞行状态下大气数据测量方法 | |
CN105486307B (zh) | 针对机动目标的视线角速率估计方法 | |
CN108061477B (zh) | 一种导引头与惯组系统间相对安装误差校正方法 | |
CN107101636B (zh) | 一种使用卡尔曼滤波器辨识多旋翼动力学模型参数的方法 | |
CN104568652B (zh) | 一种高精度测量临近空间大气密度的方法及测量装置 | |
CN112526872B (zh) | 带大落角约束中制导和末制导交接与制导信息的处理方法 | |
CN108709956B (zh) | 基于落球定位信息的大气参数测量方法和设备 | |
CN106123917B (zh) | 考虑外杆臂效应的捷联惯导系统罗经对准方法 | |
CN106017460B (zh) | 一种地形辅助惯导紧组合的水下潜器导航定位方法 | |
CN110286390A (zh) | 一种指定路径风速测量方法、装置及测风雷达标定方法 | |
CN111611648B (zh) | 助推-滑翔式飞行器最大飞行距离能力评估方法 | |
CN109855649A (zh) | 一种水下滑翔机平台运动轨迹被动确定方法 | |
CN109190248A (zh) | 一种用于滑翔飞行器的滑翔射程解析方法及解析系统 | |
CN110007318B (zh) | 风场干扰下基于卡尔曼滤波的单无人机判断gps欺骗的方法 | |
CN108827345A (zh) | 一种基于杆臂挠曲变形补偿的机载武器传递对准方法 | |
CN115248038A (zh) | 一种发射系下的sins/bds组合导航工程算法 | |
CN103486904A (zh) | 一种简易制导弹药的拟速度追踪制导方法 | |
CN116992700B (zh) | 一种物流无人机导航精度确定的方法及设备 | |
CN105973237B (zh) | 基于实际飞行数据插值的仿真动态轨迹解析生成方法 | |
CN109445283B (zh) | 一种用于欠驱动浮空器在平面上定点跟踪的控制方法 | |
De Bruijn et al. | Measuring low-altitude winds with a hot-air balloon and their validation with Cabauw tower observations |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |