1、JOURNALOFCHENGDUUNIVERSITYOFINFORMATIONTECHNOLOGYJun.20232023年6 月息工程报Vol.38No.3大学学成都3信第38 卷第3期文章编号：2 0 9 6-16 18(2 0 2 3)0 3-0 2 7 1-0 6基于气象无人机观测的测风优化研究凌荣巧，卢会国，蒋娟萍，冯岩莹（成都信息工程大学电子工程学院，四川成都6 10 2 2 5）摘要：气象无人机测风主要依托于皮托-静压管法，通过皮托管和静压管的气压差测得空速，通过CGPS等设备测得地速，根据地速、空速和风速的三角矢量关系间接测得空速。现有的测空速模型采用标准大气压下的空气密度，不
2、考虑无人机在飞行时的空气密度，这与标准大气不相符。本文通过2 0 2 1年12 月11日在四川省都江堰市的无人机飞行数据来分析空速和地速的变化。将原始空速和风速与考虑空气密度得到后的修正空速和风速对比发现：在不考虑空气密度的情况下，空速的绝对误差在-1.35-1.2 5m/s，相对误差在-6%左右；然而，传递给风速的平均相对误差却达到2 0%左右，最大达到40%。而引入空气密度可以消除这部分误差，从而提高测风精度。关键词：无人机；皮托-静压管；测风；空气密度中图分类号：TP79文献标志码：Adoi:10.16836/ki.jcuit.2023.03.0040引言实时准确的气象探测资料是开展气象
3、预报预测、防灾减灾、空中任务保障的重要一环1-2 ,气象探测资料包括温度、湿度、风向、风速、气压等气象数据，气象数据的采集是气象观测的重要工作3。目前，地基观测设备资料的空间性和连续性有很大的不足，在大范围内获取气象要素只能借助探空设备4-5。随着航空和电子技术的发展应用，利用无人机测量气象要素的趋势已经越来越明显。无人机能够持续、不间断地进行气象要素采集，具有机动性强、飞行高度高、续航时间长等优点，在气象观测应用领域将得到更多的关注6 。不过，由于无人机在飞行过程中经常会受到风的影响7 ，使无人机中关于风场参数测定的困难度很大8 ,所以如何通过无人机精确地测定风速风向已是目前无人机气象监测领
4、域的重要研发方向9 。对于空中测风，建立了高空气象探测站，每天早晚两次通过放飞探空气球来获得其他高度层上的风场信息，但是其探空资料十分有限，所以逐渐通过无人机来测风。马舒庆等10 先后提出了水平空速归零法和解析测风法来测风,屈耀红等11 在2 0 0 9 年提出了航位推算法。目前常用的测风方法是皮托-静压管法，是澳大利亚气象局DrGregHolland12在19 9 2 年提出的。针对于皮托-静压管测风法，为进一步提高其精度，减小其空速误差、地速误差、角度误差和探测“野值”等收稿日期：2 0 2 2-0 6-0 6基金项目：国家自然科学基金资助项目（42 0 7 512 9）误差来源，国内外也
5、进行了大量的研究13。任金彬等14 对无人机皮托-静压管测风做了误差分析,认为空速误差是引起风速误差的关键，而引起空速误差的关键是静压误差，可以改变无人机和皮托-静压管的外形和形状以及修正攻角和马赫数的影响来减小。周伟静等15 借助无人机皮托-静压管测风原理，建立了基于磁偏角对测风影响的数学模型,并通过仿真、实验进行相关的验证。王彦杰等16 提出了一种针对气象无人机飞行角度的优化测风方案，对测量总压、空速与地速的夹角进行了修正。金永奇等17 提出了引人加速度的无人机皮托-静压管法的测风模型，分析了无人机在飞行时加速度对测风的影响。Borup等18 利用配备了全球导航卫星系统并辅助以惯性导航系统
6、（G NSS-INS）的固定翼无人机，提出了一种指数稳定的非线性风速观测器，该观测器可以估计皮托-静压管系数并对无人机持续的机动性没有要求，研究结果表明，该观测器计算的风速在无模型误差的情况下收敛于真实风速。本文通过对现有皮托-静压管测风模型进行改进，考虑实时的空气密度，空气密度又和气温与气压有关系，对比引人实时空气密度前后风速误差和空速误差，从而说明引人空气密度的必要性。1现有测风模型当前，固定翼无人机测风主要是使用皮托-静压管法，借助于皮托-静压管测得其空速，再借助GPS等设备测得地速，被测的空速、地速和风速呈矢量三角形,再通过解三角形得到风速19 。将地速、空速和风大272息都第3 8
7、卷报程成信学学速三者映射到大地坐标系上如图1所示，得到风速的计算公式17 2v.U.COST(1)2+V2Wa其中，是风速，。是地速，。是空速，是地速与空速的夹角。图1中T是地速相对于地理北的夹角，通过GPS可获得；T2是空速相对于地理北的夹角，可通过陀螺仪获得；而地速与空速的夹角为T1和T2之差。N空速风速地速T2t图1地速、空速和风速矢量图皮托-静压管又叫空速管，由两部分组成，即总压管（皮托管）和静压管，负责测量无人机飞行时的总压和静压，一般安装于无人机的正前方。由于气象无人机飞行时的巡航速度相对于载客民用飞机飞行时的速度要低得多，马赫数小于0.3。所以，气体可认为是不可压缩的，此时可应用
8、伯努利方程：P,=Po+pu(2)式中P,为总压;P。为静压;p为空气密度常数;u是空速理想气体状态方程为：P=pRT(3)结合伯努利方程和理想气体状态方程可得到空速模型为：V2RTAP(4)V=P其中,R为气体常数，对于理想气体,R一般取2 8 7/J(KgK);T为静温（K);P=P,-Po,即总压与静压之差;P。和P均为静压。无人机的地速是相对于大地运动的速度，无人机的空速是相对于空气运动的速度，二者之间存在一定的夹角，如果在静风的情况下，夹角T为0,即地速等于空速。因此,风速的作用才使无人机的地速和空速不相等，通过测量地速和空速就可以间接测得风速。由地速、空速和风速的三者矢量关系可知,
9、风速误差主要来源于地速测量和空速测量，由图1可知，误差包括大小和方向，主要考虑其速度大小上的误差。皮托静压管测风可用于风洞中做风速测量标准传感器使用,风向的测量则需要在一个圆周范围内每隔一定的间距安装一个皮托-静压管,反演出风向，所以风向的误差较小。地速采用GPS测得，GPS的信号主要有3种码，即C/A码（粗捕获码）、P码（经码）和Y码。民用一般采用粗捕获码，其定位误差在2 9.3 2.9 3m，测速精度0.3m/s,对风速的传递误差为0.32 0 。如果风速为6m/s时，地速误差导致的空速误差很小,所以对地速误差所带来的误差不予考虑。而在空速模型计算中，空气密度是一个常数，在中国国家标准规定
10、的大气压下，空气密度为1.2 2 5kg/m，即无人机在飞行过程中空气密度都是取的这一数值，从而忽略了空气密度对测量空速的影响。而实际上，无人机飞行在不同的高度上,周围的空气密度是不同的。所以，需要考虑无人机所处的实时空气密度2 12考虑空气密度的测风模型无人机处于不同的高度时，其空气密度不同。空气密度是指在一定的温度和压力下，单位体积空气所具有的质量。参考理想气体状态方程，空气密度p的表达式为：PRT(5)由式(5），得知空气密度取决于气温和气压。在大气层中的对流层中，气温随海拔高度的增加而降低，高度每上升1km，气温就下降6.5；空气压力是指单位面积向上延伸到大气上边界的垂直空气柱的重量，
11、所以，随着海拔的升高气压随之降低。海拔高度上升，气压下降，气温也会降低，气压和气温的变化决定了空速的变化,其中气压是静压，气压根据静压管测得，而气温弃用标准大气压下的温度，采用实时的温度，根据式（2）、（3）、（5）可得到修正后的风速模型：2RT,AP(6)V=PT,是实时气温，单位为热力学温度K，而气温传感器的所测温度为摄氏温度，其转换公式为：T,=（t+273.15）,P是静压,单位为hPa。由此可知,空气密度p对于空速测量存在着一定的影响，气压和气温能反映空气密度的变化，现有的空速模型只是把气温看作是一个常数，由此空速测量造成了一定的系统误差，此误差会传递给风速，最终影响风速测量的结果。
12、考虑实时空气密度的模型，消除了系统误差。273凌荣巧人机观测的测风优化研究第3期3结果分析考虑不同海拔高度的空气密度对测量空速的影响，根据无人机测风的实际情况进行数据分析，对比原有测得的空速和风速，可看到改进后的空速模型对减小误差的效果。本文采用的无人机是纵横CW-15大鹏固定翼无人机，其主要参数如表1所示，并且搭载了可以测量气温、相对湿度、气压等气象要素的集成传感器，其安装效果如图2 所示。试验的地点位于四川省成都市都江堰市，当地的海拔高度为7 9 3m，处于东经103.7北纬31.1，试验的时间段2 0 2 1年12 月11日上午11：19-12：11，每隔1s传输一次飞行数据。表1无人机
13、的相关参数模块参数机身长度/mm2060高度/mm620翼展/mm3540最大起飞重量/kg17智能电池/mAh27000智能电池充电器/W21200无线电链路/MHz840845携带箱尺寸/mm1370（长）6 0 5（宽）52 5（高）图2纵横CW-15无人机搭载传感器3.1空速和风速的评判依据无人机以圆和椭圆的轨迹在空中飞行，从地面飞行至海拔12 0 0 m左右的空中进行实时测验。获得原始的地速、空速和风速的数据，在测空速的模型中考虑实时的空气密度从而得到修正的空速和修正的风速，分别对二者的误差进行分析。用探空设备来测风的方式比较有限，通常分为三类：一是通过国家标准气象站在每天两个时间点
14、放飞的探空气球来获得对应时间段在每个垂直高度上的探空数据，但是只能获得国家标准气象站附近区域上空的数据，而国家标准气象站的数量有限,其地理位置和探空气球的放飞时间制约了其应用的广泛性；二是通过气象雷达卫星等探空设备来反演某区域的风场情况，气象雷达设备的站点较为稀疏，数量有限，空间分辨率较低，不能做到精细化的分析；三是通过附近气象站或便携式自组网气象站来获得测风数据，但是只能测量到地面的风的情况。所以，以上3种方式对空中风的测量都提供不了参考依据。本文所采用的皮托一静压管测风可用于风洞中做风速测量标准传感器使用,其风向的测量是在无人机上安装多个皮托-静压管进而来反演出风向，只是在测量空中风时，本
15、文引人了实时的空气密度，进而在数值上对原始的风速测量有影响。因而，本文只能通过引人实时空气密度前后的风速变化来进行对比验证。选定相对误差和绝对误差来评估改进后的测风模型。绝对误差是被测量值与真值之差，即：绝对误差=测量值-真值,所以绝对误差既有量纲又有符号2 2 。O但是真值是一个理想概念，一般无法获得，实际应用中通常用非常接近被测量的实际值来代替真值。本文采用修正后的空速和风速作为真值，原始的空速和风速作为测量值。因此，绝对误差表明了被测量物体的测量值与实际值之间的偏离程度和方向。测量质量不仅与实际偏离的值有关,还和其测量的精度有关,因为一个量的准确度不仅与绝对误差有关，还和这个量本身的大小
16、有关,所以还需引人相对误差。相对误差是绝对误差和测量真值的百分比，相对误差只有大小和方向，没有单位。在本文中其测量真值仍然是修正后的空速和风速。绝对误差的绝对值和相对误差的百分比越大，表示改进后的测风模型越好。3.2空速和风速的定量计算分析3.2.1无人机飞行状态分析2 0 2 1年12 月11日的飞行数据，当天西南风，风速大小为2 级，温度在511,无人机从海拔7 9 3m的山地起飞，稳定飞行在海拔11501200m，最终降落在原地。上午11:19-12：11无人机飞行的三维轨迹如图3所示。图3中时间序列是指：每接收到一条无人机飞行数据的时间节点，通过时间序列、海拔高度以及经纬度可以实时掌握无人机的飞行状态。无人机在起飞和降落阶段为简单的直线飞行，而在平稳飞行阶段主要以绕椭圆飞行为主。图4为无人机在飞行路径上各个点的风矢量图，由于飞行轨迹有重合，所以风矢量产生重叠。虽然在同一轨迹上飞行，但是测量风速的时间有差异，其空气流动的方位和大小也有不同，所以在同一个点可能会出现风向不一致的情况。大274息第38 卷程都报成信学学400035003000250020001500100050031
