基于气象无人机观测的测风优化研究.pdf

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 0 2 1年12 月11日在四川省都江堰市的无人机飞行数据来分析空速和地速的变化。将原始空速和风速与考虑空气密度得到后的修正空速和风速对比发现：在不考虑空气密度的情况下，空速的绝对误差在-1.35-1.2 5m/s，相对误差在-6%左右；然而，传递给风速的平均相对误差却达到2 0%左右，最大达到40%。而引入空气密度可以消除这部分误差，从而提高测风精度。关键词：无人机；皮托-静压管；测风；空气密度中图分类号：TP79文献标志码：Adoi:10.16836/ki.jcuit.2023.03.0040引言实时准确的气象探测资料是开展气象预报预测、防灾减灾、空中任务保障的重要一环1-2 ,气象探测资料包括温度、湿度、风向、风速、气压等气象数据，气象数据的采集是气象观测的重要工作3。目前，地基观测设备资料的空间性和连续性有很大的不足，在大范围内获取气象要素只能借助探空设备4-5。随着航空和电子技术的发展应用，利用无人机测量气象要素的趋势已经越来越明显。无人机能够持续、不间断地进行气象要素采集，具有机动性强、飞行高度高、续航时间长等优点，在气象观测应用领域将得到更多的关注6 。不过，由于无人机在飞行过程中经常会受到风的影响7 ，使无人机中关于风场参数测定的困难度很大8 ,所以如何通过无人机精确地测定风速风向已是目前无人机气象监测领域的重要研发方向9 。对于空中测风，建立了高空气象探测站，每天早晚两次通过放飞探空气球来获得其他高度层上的风场信息，但是其探空资料十分有限，所以逐渐通过无人机来测风。马舒庆等10 先后提出了水平空速归零法和解析测风法来测风,屈耀红等11 在2 0 0 9 年提出了航位推算法。目前常用的测风方法是皮托-静压管法，是澳大利亚气象局DrGregHolland12在19 9 2 年提出的。针对于皮托-静压管测风法，为进一步提高其精度，减小其空速误差、地速误差、角度误差和探测“野值”等收稿日期：2 0 2 2-0 6-0 6基金项目：国家自然科学基金资助项目（42 0 7 512 9）误差来源，国内外也进行了大量的研究13。任金彬等14 对无人机皮托-静压管测风做了误差分析,认为空速误差是引起风速误差的关键，而引起空速误差的关键是静压误差，可以改变无人机和皮托-静压管的外形和形状以及修正攻角和马赫数的影响来减小。周伟静等15 借助无人机皮托-静压管测风原理，建立了基于磁偏角对测风影响的数学模型,并通过仿真、实验进行相关的验证。王彦杰等16 提出了一种针对气象无人机飞行角度的优化测风方案，对测量总压、空速与地速的夹角进行了修正。金永奇等17 提出了引人加速度的无人机皮托-静压管法的测风模型，分析了无人机在飞行时加速度对测风的影响。Borup等18 利用配备了全球导航卫星系统并辅助以惯性导航系统（G NSS-INS）的固定翼无人机，提出了一种指数稳定的非线性风速观测器，该观测器可以估计皮托-静压管系数并对无人机持续的机动性没有要求，研究结果表明，该观测器计算的风速在无模型误差的情况下收敛于真实风速。本文通过对现有皮托-静压管测风模型进行改进，考虑实时的空气密度，空气密度又和气温与气压有关系，对比引人实时空气密度前后风速误差和空速误差，从而说明引人空气密度的必要性。1现有测风模型当前，固定翼无人机测风主要是使用皮托-静压管法，借助于皮托-静压管测得其空速，再借助GPS等设备测得地速，被测的空速、地速和风速呈矢量三角形,再通过解三角形得到风速19 。将地速、空速和风大272息都第3 8 卷报程成信学学速三者映射到大地坐标系上如图1所示，得到风速的计算公式17 2v.U.COST(1)2+V2Wa其中，是风速，。是地速，。是空速，是地速与空速的夹角。图1中T是地速相对于地理北的夹角，通过GPS可获得；T2是空速相对于地理北的夹角，可通过陀螺仪获得；而地速与空速的夹角为T1和T2之差。N空速风速地速T2t图1地速、空速和风速矢量图皮托-静压管又叫空速管，由两部分组成，即总压管（皮托管）和静压管，负责测量无人机飞行时的总压和静压，一般安装于无人机的正前方。由于气象无人机飞行时的巡航速度相对于载客民用飞机飞行时的速度要低得多，马赫数小于0.3。所以，气体可认为是不可压缩的，此时可应用伯努利方程：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,即地速等于空速。因此,风速的作用才使无人机的地速和空速不相等，通过测量地速和空速就可以间接测得风速。由地速、空速和风速的三者矢量关系可知,风速误差主要来源于地速测量和空速测量，由图1可知，误差包括大小和方向，主要考虑其速度大小上的误差。皮托静压管测风可用于风洞中做风速测量标准传感器使用,风向的测量则需要在一个圆周范围内每隔一定的间距安装一个皮托-静压管,反演出风向，所以风向的误差较小。地速采用GPS测得，GPS的信号主要有3种码，即C/A码（粗捕获码）、P码（经码）和Y码。民用一般采用粗捕获码，其定位误差在2 9.3 2.9 3m，测速精度0.3m/s,对风速的传递误差为0.32 0 。如果风速为6m/s时，地速误差导致的空速误差很小,所以对地速误差所带来的误差不予考虑。而在空速模型计算中，空气密度是一个常数，在中国国家标准规定的大气压下，空气密度为1.2 2 5kg/m，即无人机在飞行过程中空气密度都是取的这一数值，从而忽略了空气密度对测量空速的影响。而实际上，无人机飞行在不同的高度上,周围的空气密度是不同的。所以，需要考虑无人机所处的实时空气密度2 12考虑空气密度的测风模型无人机处于不同的高度时，其空气密度不同。空气密度是指在一定的温度和压力下，单位体积空气所具有的质量。参考理想气体状态方程，空气密度p的表达式为：PRT(5)由式(5），得知空气密度取决于气温和气压。在大气层中的对流层中，气温随海拔高度的增加而降低，高度每上升1km，气温就下降6.5；空气压力是指单位面积向上延伸到大气上边界的垂直空气柱的重量，所以，随着海拔的升高气压随之降低。海拔高度上升，气压下降，气温也会降低，气压和气温的变化决定了空速的变化,其中气压是静压，气压根据静压管测得，而气温弃用标准大气压下的温度，采用实时的温度，根据式（2）、（3）、（5）可得到修正后的风速模型：2RT,AP(6)V=PT,是实时气温，单位为热力学温度K，而气温传感器的所测温度为摄氏温度，其转换公式为：T,=（t+273.15）,P是静压,单位为hPa。由此可知,空气密度p对于空速测量存在着一定的影响，气压和气温能反映空气密度的变化，现有的空速模型只是把气温看作是一个常数，由此空速测量造成了一定的系统误差，此误差会传递给风速，最终影响风速测量的结果。考虑实时空气密度的模型，消除了系统误差。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无人机的相关参数模块参数机身长度/mm2060高度/mm620翼展/mm3540最大起飞重量/kg17智能电池/mAh27000智能电池充电器/W21200无线电链路/MHz840845携带箱尺寸/mm1370（长）6 0 5（宽）52 5（高）图2纵横CW-15无人机搭载传感器3.1空速和风速的评判依据无人机以圆和椭圆的轨迹在空中飞行，从地面飞行至海拔12 0 0 m左右的空中进行实时测验。获得原始的地速、空速和风速的数据，在测空速的模型中考虑实时的空气密度从而得到修正的空速和修正的风速，分别对二者的误差进行分析。用探空设备来测风的方式比较有限，通常分为三类：一是通过国家标准气象站在每天两个时间点放飞的探空气球来获得对应时间段在每个垂直高度上的探空数据，但是只能获得国家标准气象站附近区域上空的数据，而国家标准气象站的数量有限,其地理位置和探空气球的放飞时间制约了其应用的广泛性；二是通过气象雷达卫星等探空设备来反演某区域的风场情况，气象雷达设备的站点较为稀疏，数量有限，空间分辨率较低，不能做到精细化的分析；三是通过附近气象站或便携式自组网气象站来获得测风数据，但是只能测量到地面的风的情况。所以，以上3种方式对空中风的测量都提供不了参考依据。本文所采用的皮托一静压管测风可用于风洞中做风速测量标准传感器使用,其风向的测量是在无人机上安装多个皮托-静压管进而来反演出风向，只是在测量空中风时，本文引人了实时的空气密度，进而在数值上对原始的风速测量有影响。因而，本文只能通过引人实时空气密度前后的风速变化来进行对比验证。选定相对误差和绝对误差来评估改进后的测风模型。绝对误差是被测量值与真值之差，即：绝对误差=测量值-真值,所以绝对误差既有量纲又有符号2 2 。O但是真值是一个理想概念，一般无法获得，实际应用中通常用非常接近被测量的实际值来代替真值。本文采用修正后的空速和风速作为真值，原始的空速和风速作为测量值。因此，绝对误差表明了被测量物体的测量值与实际值之间的偏离程度和方向。测量质量不仅与实际偏离的值有关,还和其测量的精度有关,因为一个量的准确度不仅与绝对误差有关，还和这个量本身的大小有关,所以还需引人相对误差。相对误差是绝对误差和测量真值的百分比，相对误差只有大小和方向，没有单位。在本文中其测量真值仍然是修正后的空速和风速。绝对误差的绝对值和相对误差的百分比越大，表示改进后的测风模型越好。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.1231.11103.7731.131.095103.76/03.765纬度31.08103.755103.7531.07103.745经度（a）飞行轨迹-时间序列12501200115011001050100095090085031.1231.11103.7731.1103.76103.76531.09纬度31.08103.755103.75经度31.07103.745（b）飞行轨迹-海拔高度图3无人机三维轨迹飞行图31.11531.1131.10531.131.09531.0931.08531.0831.075163.745103.75103.755103.76103.7655103.77103.775图42021年12 月11日上午无人机飞行路径上的风矢量3.2.2空速的误差分析考虑无人机飞行时的空气密度，改进空速模型,并且通过地速、空速和风速三者的矢量关系得到其对应的风速，改进前后的空速对比和风速对比如图5和图6所示。通过图5可知，在考虑了实时的空气密度后，修正空速相比于原始空速都有所增加,并且二者速度变化的一致性都相当吻合，这是由于随着海拔高度上升，气温和气压降低，空速随之增加。稳定飞行阶段时原始空速基本稳定在18 m/s左右，而改进后这一数值达到了19 m/s左右；原始空速的平均速度是18.3m/s，修正后的空速为19.5m/s，飞机在起飞和降落阶段速度变化较为剧烈，一方面是由于飞机加速度的作用,另一方面是由于飞机一直做的是绕圈飞行，在风速的作用下，不同的转弯过程（顺风转弯或者逆风转弯）会产生不同的加速度从而影响飞行时的速度。由图7 可知，无人机在飞行时的空速绝对误差大多分布在-1.35-1.25m/s,相对误差保持在-6%左右。2120.520(s/u)/车19.519区18.51817.51716.55001000 150020002500300035004000时间序列图52021年12 月11日上午无人机空速改进前后对比76(s/u)/4321050010001500200025003000 35004000时间序列图6 2 0 2 1年12 月11日上午风速改进前后对比120010008006004002001.4 1.351.3 1.25-1.21.151.11.0510.950.9（a）空速绝对误差150010005000.068-0.066-0.064-0.062-0.060.058-0.056-0.054-0.052-0.05(b)空速相对误差图7空速绝对误差和相对误差直方图3.2.3风速的误差分析综上可知，考虑了空气密度后修正所测得的空速，其绝对误差和相对误差都比较小，但是其产生的空速误差将传递给风速的测量。由于无人机飞行是近地绕圈飞行，并且飞行的速度不是很快，逆风飞行时其地速减小,顺风飞行时其地速增加,所以修正后的风速较原始风速有所波动。由图6 可知，在时间序列12 0 0-1500和2 10 0-2 40 0 时段，无人机正从西往东飞行，处于顺风的状态,地速增加，但是空速没有什么变化,所以风速的幅值下降，但其变化趋势与同一时刻的原始风速一致。稳定飞行阶段所测得的原始平均风速为4.86m/s，修正后的风速的平均值为5.30 m/s。由图8可知,风速绝对误差分布在-1.41.4m/s,相对误差分布在-2 5%2 5%。由此可知，忽略空气密度所引起的风速测量相对误差比较大，平均的相对误差都在2 0%左右，有时甚至可以超过40%275凌荣巧（机观测的测风优化研究第3期100080060040020001.5-1-0.500.511.5（a）风速绝对误差12001000800F60040020004-0.80.60.40.200.20.40.60.81(b）风速相对误差图：风速绝对误差和相对误差直方图4总结与讨论本文针对无人机皮托-静压管现有测风模型进行改进，考虑将无人机在飞行时的空气密度加人到测风模型中。空气密度与当时所处的气压与气温有关，把实时的气压与气温考虑进测风模型中，修正其对应的空速。由于地速、空速和风速构成矢量三角形，空速的改变会引起风速的变化，从而通过矢量三角形关系来修正风速。修正后的空速较原始空速有所增大，但是其变化趋势完全一致，平均速度从18.3m/s上升到19.5m/s,其绝对误差在-1.35-1.2 5m/s,所带来的相对误差在-6%左右。然而，传递给风速的测量变化较大，但是其风速变化的趋势与原始风速基本一致，其绝对误差在分布在-1.4 1.4m/s,相对误差在20%左右，误差较大的时刻甚至达到了40%左右。所以，考虑空气密度对无人机测风是不可忽视的。总体而言，考虑空气密度对无人机皮托-静压管测风可以提高测风精度，但是还应该考虑基于此方法所带来的其他误差，如地速误差、角度误差，改进测风模型，提高测风精度。参考文献：1朱嘉慧，王海江，李静，等.Meteo-particle模型在ADS-B风场反演中的性能研究J.成都信息工程大学学报,2 0 2 1,36(5)：47 9-48 4.2Alaoui-Sosse S,Durand P,Medina P,et al.BORE-AL-A Fixed-Wing Unmanned Aerial System forthe Measurement of Wind and Turbulence in theAtmospheric Boundary Layer J.Journal of At-mospheric and Oceanic Technology,2022,39(3):387-402.3宁志远，刘厚凤.大气边界层的国内外研究现状J.中国环境管理干部学院学报,2 0 17,2 7(2)：22-25.4Roseman C A,Argrow B M.Weather Hazard RiskQuantification for sUAS Safety Risk ManagementJ.Journal of Atmospheric and Oceanic Technol-0gy,2020,37(7):1251-1268.5何建新，张福贵，周红根，等.龙卷风探测雷达研制及业务化应用研究J.气象科技进展,2 0 2 1，11(4):54.6沈怀荣，邵琼玲，王盛军等编著.无人机气象探测技术M.北京：清华大学出版社,2 0 10.7祁月皎，何建新，王旭.多普勒天气雷达二维理想均匀风场的数值模拟J.成都信息工程学院学报,2 0 14,2 9(2):12 7-132.8Pinto J O,O Sullivan D,Taylor S,et al.The Statusand Future of Small Uncrewed Aircraft Systems(UAS)in Operational MeteorologyJ.Bulletin ofthe American Meteorological Society,2021,102(11):2121-2136.9董天天.基于无人机的大气边界层气象要素探测平台设计D.南京：南京信息工程大学,2 0 17.10马舒庆，汪改，潘毅.微型无人驾驶飞机探空初步试验研究J.南京气象学院学报19 9 7(2)：30-36.11屈耀红，凌琼，闫建国，等.无人机DR/GPS/RP导航中风场估计仿真J.系统仿真学报，2009,21(7):1822-1825.12Holland GJ.The Aerosonde robotic aircraft:Anew paradigm for environmental observations J.Bulletin of the American Meteorological Society,2001,82(5):889-901.13胡昊辉.基于六旋翼无人机的测风系统研究D.长沙：湖南大学,2 0 2 0.14任金彬，沈怀荣.无人机皮托-静压管测风误差分析J.装备指挥技术学院学报,2 0 0 4(4):45-48.15周伟静，沈怀荣.磁偏角对无人机皮托-静压管测风的影响分析J.装备指挥技术学院学报，2006(4):97-101.16王彦杰，司长彬.基于气象无人机飞行角度的优化测风模型研究J.传感器世界，2 0 12，(8):16-20.17金永奇，周树道，卫克晶，等.引入加速度的无大276息报都成程第38 卷信学学人机皮托-静压管法测风模型J.探测与控制学报,2 0 12,34(6):7 2-7 5.18Borup,Kasper T,Fossen,et al.A Nonlinear Mod-el-Based Wind Velocity Observer for UnmannedAerial Vehicles J.IFAC-PapersOnLine,2011,49(18):276-283.19周伟静，沈怀荣.一种基于小型无人机的风场测量方法J.测试技术学报,2 0 0 9,2 3(4)：2 9 7-30 2.20李天文.GPS原理及应用M.北京：科学出版社,2 0 0 3.21程春龙，潘晓春.空气密度时序变化特征及其对风压计算的影响J.黑龙江科学,2 0 19,10(22):1-3.22王川.电子测量技术与仪器M.北京：北京理工大学出版社,2 0 19.Optimization Study of Wind Measurement based onMeteorological Drone ObservationLING Rongqiao,LU Huiguo,JJIANG Juanping,FENG Yanying(College of Electronic Engineering,Chengdu University of Information Technology,Chengdu 610225,China)Abstract:Meteorological UAV wind measurement is mainly based on the Pitot-Static tube method,which measures theair speed through the Pitot tube and static tube air pressure difference.According to the triangular vector relationship ofground speed measured by GPS and other equipments,air speed and wind speed,air speedcan be measuredindirectly.The existing airspeed measurement model uses the air density at standard atmospheric pressure,and does not considerthe air density of the UAV during flight.In this paper,we analyze the variation of airspeed and ground speed by theflight data in Dujiangyan,Sichuan Province on December 11,2021.Comparing the original airspeed and wind speedwith the corrected airspeed and wind speed considering the air density,it is found that:the absolute error of airspeedwhich ignores the influence of air density is between-1.35-1.25 m/s and the relative error is around-6%;howev-er,the average relative error passed to the wind speed reaches around 20%and the maximum reaches 40%.And theintroduction of air density can eliminate this part of error,thus improving the wind measurement accuracy.Keywords:UAV;Pitot-Static tube;wind measurement;air density