改进粒子滤波器的室内地磁匹配精度研究

针对粒子滤波算法在室内地磁分布的连续性变化中会发生滤波发散,加之算法中重采样步骤使该现象恶化,直接导致室内定位过程中目标丢失的问题。该文受Hausdorff距离度量法中点集匹配思想的启发,提出了一种利用初定位误差作为最大粒子匹配距离的约束思想,通过缩小匹配空间的区域范围的方法来解决传统算法的发散问题。最后实验表明,此方法解决了粒子滤波器在地磁匹配中出现的滤波发散现象,同时使得定位结果更加快速收敛。     

基于有效区域约束的GPU-CPU协同影像快拼方法

提出了一种以有效区域约束的Voronoi图为拼接线网络,利用GPU-CPU协同处理航空影像快速拼接方法。首先,基于成像有效区域生成有效区域约束的Voronoi图拼接线网络,解决了传统Voronoi图拼接线网络在低重叠度条件下拼接后影像局部区域不被覆盖的问题,然后利用GPU-CPU协同处理将正射纠正嵌入到影像拼接的过程中,并且只对每张影像的有效区域进行纠正,再通过选择配置优化和存储层次性优化进一步提高拼接效率。实验表明,对237张高分辨率航空影像进行快速正射纠正和拼接,本文算法较传统先纠正再拼接的方法效率提高近20倍,同时保证很高的拼接精度,可以满足应急测绘要求。     

全球CGCS2000坐标框架的构建理论研究

全球坐标参考框架是发展全球卫星导航系统最具重要性的基础设施之一,也是我国空间信息基础设施建设的迫切需求。研究了构建全球CGCS2000（China Geodetic Coordinate System 2000）坐标参考框架的数学模型,研究了采用内约束理论构建全球CGCS2000框架理论和算法,最后证明了构建全球CGCS2000框架的二步法和一步法数学模型在一定条件下等价。内部约束算法根据框架转换参数及其速率变化量最小标准重新定义框架约束条件,经过坐标相似变化得到最佳网形,并保证最佳网形和初始网形的拟合度达到最高,最后得到最优网形中测站坐标估值,并能保持不同观测技术参考框架的内在纯粹性,发挥不同观测技术优势。本文研究结果可完善中国坐标框架构建基础理论体系,对CGCS2000框架的全球化实现提供科学方法。     

附不等式约束平差模型的一种快速搜索算法

大地测量中常存在一些先验不等式约束信息,充分利用它们可以保证参数解的唯一性和稳定性。然而,现有的不等式约束平差算法主要是基于优化理论,算法通常比较复杂,需要选取有效约束或建立罚函数。在最小二乘平差准则基础上,把不等式约束看成是一个可行域,借助Fisher函数在可行域中快速搜索使误差平方和达到最小的最优解,推导出了可行解为最优解的充分必要条件。建立了基于Wolfe-Powell算法的非精确快速搜索算法,从而减小了搜索算法的计算量,得到了一种新的不等式约束平差计算方法。该算法的平差准则与最小二乘平差准则一致,不需要矩阵求逆运算,可适用于维数较大的平差问题解算。     

九龙站在西南低涡/晴朗天气背景下观测资料的对比分析

本文基于中国气象局成都高原气象研究所2017年西南涡加密观测资料,对比分析了西南低涡主要源地九龙站在西南低涡/晴朗这两种不同天气背景下,各个基本物理量的日变化特征及存在的差异,认为:西南涡天气背景下,边界层高度较低,而晴朗天气状况下,边界层大气湍流运动强烈,对流混合伸展高度非常惊人,能超过3000m。西南涡天气背景下,夜间大气较为暖湿,气压显著偏低,上层大气风速较大,便于动量下传。晴朗天气背景下,大气较为干燥,温度日较差偏大,大气压较高,白天湍流混合强。越靠近地面差异越显著,差值随高度的增加而逐渐减小。     

GRAPES全球奇异向量方法改进及试验分析

基于总能量模的奇异向量扰动常用于构造集合预报的初始条件。以建立GRAPES（Global and Regional Assimilation PrEdiction System）全球集合预报系统为目的,基于前期研发的GRAPES全球模式奇异向量方法,在GRAPES全球切线性模式和伴随模式2.0版的框架下,开展了引入线性化边界层方案来改善奇异向量结构,并提高奇异向量计算效率的研究。通过连续试验,从奇异向量的扰动能量结构、扰动能量谱及扰动空间分布等方面,综合分析改进GRAPES全球奇异向量的结构及演变特征。试验结果表明,改进后的GRAPES奇异向量方法有效抑制了之前扰动能量在近地面层不合理的快速增长,同时,奇异向量最优扰动的结构更客观地体现了中高纬度区域大气初始条件中的斜压不稳定扰动及其演变,如在初始时刻奇异向量扰动能量主要位于对流层中层,并呈现出随高度向西倾斜的大气斜压特征;经过线性化演变,扰动能量向较大水平尺度转移,并在垂直结构上表现出向对流层高层上传及向对流层低层下传的特征等。针对GRAPES奇异向量迭代求解中伴随模式计算耗时为主的情况,改进伴随模式中广义共轭余差方案的调用方式,并采用大内存存储法来提高其计算效率,进而将奇异向量总计算时间缩短了25%。总之,改进后的GRAPES奇异向量方法,可应用于构建面向业务应用的GRAPES全球集合预报系统。      

辽宁绥中地区海陆风的多普勒测风激光雷达观测与特征提取

2019年3月,利用相干多普勒测风激光雷达首次在辽东湾西部绥中地区进行了风廓线测量试验。根据研究区域海岸线走向采用风向的十六分位法定义局地海风和陆风,分析和提取海陆风特征验证了多普勒测风激光雷达在春季季风间断期间观测海陆风的可行性,并计算和分析了大气边界层湍流能量的变化以及回流水平变化等特性。结果表明：1）绥中地区春季存在明显的海陆风环流特征,测风激光雷达观测海陆风出现的时间与地面自动气象站观测的数据较为一致,符合海陆风日的定义。2）海陆风日发生时,水平局地回流指数（RF）较小,1.2 km以下的RF值小于0.5,使得污染物循环累积,较易形成雾霾天气;但是海风时大气边界层的高度可达1 km以上,有利于低层大气污染物向高层扩散,减轻低层大气污染。研究结果为该地边界层参数化方案的设计和污染的防治提供了参考依据。     

基于WRF模式的风电场短期风速集成预报方法研究

风能始源于大气的运动,具有很大的随机性和间歇性。风速预测是风电场风功率预测的基础,其准确性具有重要的意义。对于复杂地形条件下,风速的预报一直是各国研究的难点和重点。为了提高风电场短期风速预报的准确性,本研究采用多种边界层参数化方案来集成预报风速,将各单一边界层参数化方案预报的风速及相应的实测风速数据,应用随机森林算法建立集成预报模型,对风电场的短期风速进行集成预报研究。试验结果表明,采用集成预报风速方法,预报的风速误差相比于单一边界层参数化方案预报的风速误差明显减小,对研究区域的风速、风向等气象要素有着较好的模拟效果,能够有效提高风速预报的准确率。     

中国西北地区大气边界层高度变化特征——基于探空资料与ERA-Interim再分析资料

利用中国西北地区2015年9月至2016年8月38个站点L波段探空观测、2016年7月加密探空观测和ERA-Interim边界层高度资料,对比分析了西北地区大气边界层高度变化特征。观测资料表明,在中国西北地区,08:00(北京时,下同)冬季边界层高度最高; 20:00春季边界层高度最高,边界层高度从西部到东部有显著降低的趋势。ERA-Interim资料基本能表现出边界层高度的区域分布,但相对于探空观测得到的边界层高度,除夏季20:00外,ERA-Interim再分析资料边界层高度均偏低。全年平均而言,08:00(20:00)偏低160 m(170 m),其中在08:00(20:00),冬季(春季)偏低最显著。08:00边界层高度与低层稳定度、近地层温度和风速相关更加显著; 20:00边界层高度与低层稳定度和相对湿度相关更加显著。2016年7月加密观测资料对比表明,ERA-Interim资料的对流(中性)边界层高度显著偏高;低层稳定度、相对湿度偏小,风速偏大可能是造成边界层高度偏高的原因; ERA-Interim资料的稳定边界层高度偏低,与低层稳定度和近地层温度偏低相关,但其影响因素相对更加复杂。     

边界层参数化对海南岛海风环流结构模拟的影响

利用WRF V3.7详细分析了应用8种边界层参数化方案（YSU、MYNN2.5、MYNN3、ACM2、BouLac、UW、SH、GBM）所模拟的2014年5月25日海南岛海风环流结构的差异,其中YSU、ACM2和SH为非局地闭合方案,MYNN2.5、MYNN3、BouLac、UW和GBM为局地闭合方案。结果表明:对于海风环流水平结构的模拟,15时,YSU、ACM2、BouLac、UW和SH模拟的北部海风较强,SH和GBM的内陆风速偏大。温度与海风发展强度相对应,MYNN2.5与MYNN3模拟的岛屿温度偏低,海陆温差小,海风相对较弱。对于海风环流垂直结构的模拟,09时海风开始,但强度较小,且存在残余陆风,向内陆传播距离较短,YSU、MYNN2.5和SH方案的海风相对较强。12时,海风已呈现出较为清晰的环流结构,YSU和ACM2的海风厚度及向内陆传播距离相对强于其它方案,MYNN3的环流结构则不太明显,且向内陆推进距离短,海风相对较弱。15时,海风发展强盛,MYNN2.5和MYNN3方案模拟的海风垂直强度较小,ACM2方案的海风垂直环流特征最为明显。18时,海风的强度和扰动均有所减弱,ACM2、BouLac和UW的整体海风相对强于其它方案。21时海风已基本转为陆风,BouLac与UW的陆风环流结构最为清晰。位温、水汽及海风垂直环流强度的发展变化与海风的演变过程基本一致。造成ACM2模拟海风偏强的原因是其边界层垂直混合偏强,形成了足够的湍流混合强度所致。对于边界层高度的模拟,ACM2的边界层顶最高,这与此方案所模拟的海风强度偏大相吻合,其它方案的边界层高度与海风强度并不完全一致。