GNSS约束的2022年泸定M6.8地震滑动分布及同震应力变化

2022年9月5日四川泸定县发生M6.8级强震,震中位于鲜水河断裂东南末端,野外地质调查初步结果显示本次地震并未发现明显的地表破裂迹象.本文基于震后科考GNSS流动观测和震中周边连续站观测资料,解算并提取震中90 km范围内31个测站的静态水平向同震位移.结果显示:GNSS同震形变场空间分布呈现明显的左旋走滑特征,观测到的最大水平向同震位移达23 cm,震中距50 km范围内同震位移量普遍大于1 cm.基于GNSS观测资料反演得到的同震滑动分布显示泸定地震地表破裂主要集中在磨西至田湾之间,主破裂深度2~8 km,最大滑动量~1.96 m,地震矩9.25×10^(18)N·m,对应矩震级M_W6.6.静态库仑应力结果显示本次地震增强了震源破裂区周边活动断层的库仑应力,并触发了大量余震,余震主要发生在库仑应力增强区域.结合震间闭锁分布、历史地震及库仑应力变化,我们认为未来需要密切关注与磨西断裂交接的安宁河、大凉山断裂以及康定—磨西段的地震危险性.

GNSS系统时间偏差参数分析

分析各GNSS在系统时差参数设计上的特点和变化,并以国际权度局(BIPM)和时间频率公报发布的数据为参考对GNSS播发的系统时差进行评估分析。结果表明,BDS播发的3种系统时差均存在十几ns的偏差,通过BDS间接获得的GLONASS与GPS、Galileo与GPS的系统时差准确度优于GLONASS和Galileo系统电文发播的同种时差参数。

GNSS辅助风云三号卫星MERSI近红外通道的大气可降水量反演方法

大气水汽是对流层中的重要参数之一,已被广泛应用于短临天气预警和长期气候监测等领域。我国风云三号(FY-3)系列卫星搭载的中分辨率光谱成像仪可用于大气水汽监测,但反演大气可降水量(PWV)时存在大气透过率参数低估、水汽与大气透过率回归系数经验选取的缺陷,无法满足在短临降雨监测、数值同化等高精度PWV应用方面的需求。针对该问题,本文提出一种GNSS辅助FY-3卫星的高精度PWV反演方法。本文方法引入高精度实测GNSS PWV作为大气透过率计算模型的回归拟合参数,辅助FY-3 L1级数据精确估计PWV和大气透过率的模型回归系数;同时,该方法顾及季节和高程因素对FY-3-L1 PWV反演的影响,分季节反演PWV并引入数字高程模型修正由于部分大气透过率参数低估导致的FY-3 L1 PWV相对不准确的现象。选取中国区域2013—2014年FY-3A卫星的L1数据和中国地壳运动观测网络的260个GNSS测站数据进行试验。结果表明,本文提出的GNSS辅助FY-3系列卫星PWV反演方法优于传统方法(FY-3A-L2 PWV),其整体精度改善率为74.5%,可得到更加可靠、稳健性更强的PWV格网产品,对于其在短临降雨监测和数值同化等方面具有重要意义。

因子图框架下里程计辅助GNSS/INS组合导航算法

在复杂观测环境下,GNSS/INS组合导航系统的GNSS信号易受干扰从而导致INS独立导航精度迅速下降。针对上述问题,本文基于因子图的里程计辅助GNSS/INS组合导航算法,利用里程计观测信息结合非完整性约束构建航向速度约束方程,同时采用能多次线性化计算和多次迭代的因子图优化方法进行参数估计。实际车载试验解算结果表明,在GNSS信号良好时,基于因子图方法比滤波方法具有更快的收敛时间,收敛速度提高了近10倍;在GNSS信号发生中断时,添加里程计辅助后组合导航系统在东向和北向分别提升了83%和89%。与传统的滤波融合手段相比,本文采用因子图优化后在东向和北向的定位精度分别有63%、70%的改善。

新疆及邻区现今GNSS变形特征与地震关系研究

基于GNSS观测数据,采用球面最小二乘配置方法计算了新疆及邻区的应变率张量特征,在统计分析研究区域应变率分布以及1900—2022年M≥6.0地震分布之间的关系的基础上,研究了GNSS应变率特征对强震地点的指示意义。通过构建三维弹性块体模型反演得到的区内主要断裂的运动变形特征,结合震源机制解结果,对比分析了断裂运动变形特征与不同区域强震类型之间的关系。结果表明:天山西段—帕米尔地区、阿尔金断裂带邻近地区呈现第二应变率高值特征,其中帕米尔构造结附近高值特性最明显;最大剪应变率方向在南天山西段—帕米尔地区主要为NE-NEE向,反映了该区以倾滑变形为主的特征。研究区域的M≥6.0地震主要分布在应变率高值区及其边缘,特别是帕米尔构造结东部地区的强震集中非常明显。区内断裂运动性质具有明显的分类特征,除整体以挤压运动为主外,NE走向断裂带以左旋走滑运动为主,NW走向的断裂以右旋走滑运动为主。NE走向的柯坪、迈丹和那拉提断裂带与NW走向的克孜勒陶、塔拉斯—费尔干纳断裂带汇集的南天山西段和帕米尔西构造结东部地区强震密集分布,强震的震源机制类型与断层运动方式较为一致。

基于自适应高斯渐进滤波的工程车GNSS/INS紧组合定位

研究了量测野值影响下的工程车GNSS/INS紧组合定位问题,提出了一种基于自适应高斯渐进滤波的车辆定位方法。首先,为降低量测野值对滤波器的破坏风险,利用假设检验方法对量测野值进行检测和剔除;其次,对于野值漏检测引起的定位性能下降的问题,设计了自适应的高斯渐进滤波方法来补偿量测的不确定性;特别地,利用线性化误差与系统估计误差的变化关系,对渐进量测更新方式进行了改进,从而实现对线性化误差的间接补偿。最后,通过工程车GNSS/INS紧组合定位实验进行结果分析,验证了所提方法的可靠性和优越性。

欧美月球GNSS规划现状分析综述

月球是地球最重要的天然卫星,当前国际上正在迎来新一轮月球探索高潮,数十个机构和商业团队正在规划月球探索任务,并设想在未来实现航天员长期驻月,围绕月球的“太空竞赛”刚刚开始。月球GNSS(基于现有的地球GNSS以及新的环月卫星通信导航基础设施的月球卫星通信导航定位技术)是空间基准科研的基础,能够提供航天器着陆定位以及月面(及其覆盖空间)定位、导航与授时等服务,同时可以将月球作为试验场,将导航工具包扩展到更远的目的地(如火星)。对欧美近期发布的月球GNSS规划进行了整理归纳,其中包括美国月球GNSS接收机实验(LuGRE)计划和欧洲月光(MoonLight)计划,以及美国中远期月球通信中继和导航系统(LCRNS)计划,这些计划可以为我国开展月球GNSS规划提供参考。

基于高度角随机模型的GNSS外辐射源雷达定位算法

针对GNSS外辐射源雷达定位时不同卫星对定位的误差贡献不同的问题,提出基于高度角随机模型的定位算法,理论分析目标位置估计量的克拉美罗界和统计特性,并计算卫星位置误差和地面站位置误差对定位误差的影响。仿真结果表明:所提出算法对不同GNSS卫星的直射、反射路径的伪距差测量值误差进行了合理分配,定位性能达到了克拉美罗界,且不会因为选星方案的改变而大幅恶化。对地面站位置误差和卫星位置误差的分析表明:地面站位置的标准差小于10 cm、卫星位置的标准差小于1 km时对定位总误差的贡献可以忽略不计。

基于GNSS的西秦岭北缘断裂带现今活动特征研究

基于1999~2007年、2009~2013年、2013~2017年和2017~2021年共4期GNSS地壳运动速度场构建GNSS速度剖面,分析西秦岭北缘断裂带各次级断裂现今地壳差异运动特征,并进一步建立负位错模型反演获取各次级断裂现今闭锁程度与滑动亏损速率。结果表明:1)西秦岭北缘断裂带现今整体呈现出以左旋走滑为主、兼逆冲挤压的活动特征。2)1999~2007年,锅麻滩段和天水-宝鸡段整体处于强闭锁状态,漳县段和鸳凤段仅部分区段闭锁且闭锁程度较低,闭锁深度也较浅;2009~2013年,西秦岭北缘断裂带中段闭锁程度有所减弱,鸳凤段东部与天水-宝鸡段东部由较强闭锁转变为蠕滑状态,锅麻滩段和天水-宝鸡段中西部依然呈强闭锁状态;2013~2017年,西秦岭北缘断裂带整体处于较强闭锁状态,仅天水-宝鸡段东部呈蠕滑状态;2017~2021年,锅麻滩段中西部、漳县段、鸳凤段东部和天水-宝鸡段中部处于较强闭锁状态,鸳凤段中西部与天水-宝鸡段西部则呈蠕滑状态。3)各次级断裂滑动亏损速率整体分布特征与断裂闭锁程度基本一致。西秦岭北缘断裂带锅麻滩段和天水-宝鸡段中部现今处于较强闭锁状态,具有较强的应变能累积背景。

基于GNSS的2023年2月23日塔吉克斯坦M_(S)7.2地震前异常资料变化特征研究

为了进一步探究中强地震发生前GNSS坐标时序中是否存在中短期异常现象,通过分析2023年2月23日塔吉克斯坦M_(S)7.2地震前后中国大陆构造环境监测网络(简称:“陆态网络”)塔什库尔干、布伦口GNSS基准站观测资料的变化特征,发现该地震震前塔什库尔干、布伦口均出现较为明显的GNSS基准站时间序列异常现象,震前EW向存在趋势性W向运动,震后E向加速;震前NS向趋势性S向运动,临震前状态快速改变,布伦口在震后还显示出7.3 mm的EW向永久性形变,研究结果表明GNSS在地震发生短临阶段具有重要的预测意义。与2015年12月7日塔吉克斯坦M S7.4地震进行对比发现,两次M S≥7.0地震对塔什库尔干、布伦口基准站的影响差异较大。2015年12月7日塔吉克斯坦M S7.4地震为左旋走滑,影响方向主要为NE向,而2023年2月23日塔吉克斯坦M_(S)7.2地震为右旋走滑,影响方向主要为SE向。震前基线均有挤压缩短的前兆异常显示。梳理站点异常出现情况,发现异常持续时间越长、异常出现次数越多,发震震级可能越大。

卡斯卡迪亚慢滑移信息的GNSS时空探测

针对断层慢滑移时空分布信息探测困难的问题,基于全球卫星导航系统(GNSS)坐标序列,提出一种利用多通道奇异谱分析(MSSA)探测卡斯卡迪亚消减带慢滑移事件时空分布的方法:根据慢滑移分量的振荡特点确定窗口长度,采用时间迟滞矩阵对协方差阵进行增广,以提高对异常信息的识别能力;对振幅归一化获取空间响应,通过快速傅里叶变换分析慢滑移频谱特征。结果表明,卡斯卡迪亚消减带2007-01-05—2007-02-23和2008-04-22—2008-06-14发生了2次慢滑移事件,2次事件呈东南、西北区域的反向运动特征,且东南部响应程度明显高于西北地区;频谱特征显示,慢滑移信息表现为低频特征,且主要为随机游走噪声。研究结果可为板内地震前兆性信息探测提供参考。

融合BiLSTM和VMD的GNSS坐标时间序列重构

GNSS观测时间序列包含复杂的非线性构造运动,如地面质量荷载、模型残差、周围环境因素等。由于环境因素的复杂性,季节性信号可能具备准周期时变的特征,传统的时间序列分析模型很难模型化。因此,可以采用一种双向长短期记忆(Bidirectional Long Short-Term Memory,BiLSTM)循环神经网络与变分模态分解(Variational Mode Decomposition,VMD)联合的信号重构方法。首先利用VMD强大的分解能力将GNSS信号进行频域剖分并将其分为多项子信号和噪声项,再基于BiLSTM强大的学习能力对GNSS信号进行训练建模。结果表明,BiLSTM+VMD模型能充分挖掘信号的时频域特征,提高信号重构的精度和稳定性,GNSS N、E、U三分量重构结果均方根误差(Root Mean Squared Error,RMSE)都表现出不同程度的降低,尤其水平方向效果更为显著,相比EMD与VMD方法,E方向离散度分别降低了61%和19%,N方向离散度分别降低了20%和14%。这为GNSS观测时间序列中信号提取与模型参数估计提供了一个有价值的模型。

GNSS/SINS/视觉导航鲁棒算法

全球导航卫星系统(Global Navigation Satellite System, GNSS)、捷联惯性导航系统(Strapdown Inertial Navigation System, SINS)和视觉传感器优势互补,3者信息融合可获得高精度、无漂移的导航定位信息.针对GNSS/SINS/视觉融合导航易受运动速度、光照变化、遮挡等影响导致定位精度和鲁棒性降低问题,本文在图优化框架的代价函数中加入SoftLOne鲁棒核函数,设置量测值粗差检验程序,降低离群点带来的负面影响.进一步,对量测值计算残差进行卡方检验,对超限残差降权处理,提高系统精度和鲁棒性.实验结果表明,本文算法较不施加鲁棒核函数、不采用异常值剔除策略和卡方检验的传统算法,以及加入其他鲁棒核函数的算法精度更高、鲁棒性更好,能够较大程度提升GNSS/SINS/视觉导航定位精度和鲁棒性,在大尺度环境下,未出现较大漂移误差,绝对位姿均方根误差0.735 m,绝对位姿误差标准差0.336 m.

vbICA方法用于GNSS坐标序列共模误差提取研究

全球导航卫星系统(global navigation satellite system,GNSS)坐标序列精度主要受共模误差(common mode error,CME)影响。为提高GNSS坐标序列精度,采用变分贝叶斯独立分量分析方法(variational Bayesian independent component analysis,vbICA)提取实验区20个GNSS测站坐标序列的CME,并对比分析vbICA与主成分分析(principal component analysis,PCA)和独立分量分析(independent component analysis,ICA)的滤波性能。结果表明,vbICA滤波效果优于PCA和ICA;经vbICA滤波后,E、N、U方向坐标残差序列均方根(RMS)平均降低36.57%、31.63%、10.97%,距离相关系数平均降低60.53%、56.84%、25.80%;扣除CME后,GNSS速度场估计更加可靠和精确,可有效提高GNSS坐标序列精度,为地球动力学研究提供可靠的数据支撑。

基于图优化的GNSS/双目视觉/惯性SLAM系统开发及应用

为提高机器人室外长航时定位精度,提出一种基于图优化的全球导航卫星系统(GNSS)/双目视觉/惯性同时定位与建图(SLAM)系统开发及应用。将空间中的线特征作为几何约束的补充,集成至前端的特征提取及后端的位姿优化线程,提升位姿解算精度。同时,以因子图构建联合优化的图结构,并推导出全局观测误差模型。近200 m的BullDog-CX机器人巡检结果表明,所提算法相比于VINSFusion和PL-VINS分别取得约12.6%及3.4%的定位精度提升,为室外机器人长航时导航提供了一种可行方案。

GNSS海洋反射信号接收机设计与分析

针对全球导航定位系统(global navigation satellite system,GNSS)海洋反射信号遥感技术的需求,文章设计了一种GNSS海洋反射信号接收机,用以对GNSS直射信号和反射信号进行接收。多普勒延迟映射结构是GNSS海洋反射信号接收机中的关键部件,采用直射信号通道闭环、反射信号通道开环的方式,反射信号通道利用直射信号通道获得的扩频码和多普勒频率,在半个C/A码片时间内进行相关积分,得到反射信号的相关峰波形。由于海面反射的GNSS信号是多个不同反射信号的集合,合成码的相关波形和直射信号不同,其位置和形状包含了海洋学参数信息,通过此相关函数波形即可解算出海面风场信息。GNSS海洋反射信号接收机技术可以为海洋环境监测、基础测绘、灾害监测等领域提供技术支持。

新工科背景下GNSS课程项目式教学模式实践与改革

结合项目式教学的特点,以成都理工大学GNSS原理及应用课程建设为例,对项目式教学模式建设中的项目设计、教学过程设计、教学考核评价等方面的改革进行探索,总结实践经验,以期提高教学质量。

月球导航星座轨道的地球GNSS可见性评估

随着月球探索进入一个新的时代,为确保未来月球任务的实时高精度定位,大幅提高登月航天器的自主性,且减少对地球基站的依赖,选取合适的轨道构建月球导航星座,并且实现月球导航星座与地球GNSS的通信链路同步尤为重要。为评估月球导航星座轨道中卫星接收GNSS信号的性能,首先对椭圆形月球冻结轨道(ELFO)、近直线晕轨道(NRHO)、顺行圆形轨道(PCO)和低月球轨道(LLO)4种轨道进行仿真。然后基于天线方向图等指标仿真了GNSS卫星信号,并依据主旁瓣波束宽度和载噪比等仿真结果评估了4种不同轨道卫星对GNSS的可见性。结果表明,NRHO和ELFO对GNSS星座有较好的可视效果,最高可见时间占比达99.57%,保障了绝大部分时间内能够稳定接收GNSS信号,有助于实现月球导航星座轨道卫星的轨道和时钟校正,并保证轨道的定位性能。

GNSS-I/MR河流边界及水位测量理论与仿真

从几何构型和理论模型2个角度探究全球导航卫星系统干涉/多径反射测量(GNSS-I/MR)技术提取河流边界和测量河流水位的可行性。基于河流/河岸介电常数的差异性和河岸相对于河流存在斜度的事实,定义对河流边界敏感的反-直比和实验变化率2个观测量。在希尔伯特变换估计的载噪比包络和相位的基础上,提出基于一元二次方程求解反-直比和线性拟合实验变化率的估计方法,以及3-阈值最大类间方差识别河流/河岸和循环拟合提取河流边界的方法。基于河流表面的实验变化率测量河流水位,通过搭建仿真平台和设计2种河流场景,仿真验证所提方法的有效性。仿真结果表明:在设定参数下,利用介电常数检测场景1和场景2河流边界的均方根误差分别为1.47 m和1.13 m;实验变化率检测场景1和场景2河流边界的均方根误差分别为0.49 m和0.35 m;场景1和场景2的测量高度和真实值的误差分别为0.09 m和0.13 m。

一种基于智能算法的GNSS高程拟合方法

广义回归神经网络(GRNN)是一种新型的前馈神经网络模型,具有训练次数少、耗时短、非线性参数的预报能力较强等优点。但GRNN唯一的调节参数SPREAD不能自动获取限制其进一步的应用。针对该缺陷,本文采用果蝇优化算法(FOA)与GRNN相结合构建FOAGRNN模型对GRNN进行优化,自动获取调节参数的值。为了检验FOAGRNN模型的GNSS高程拟合精度,进行了实验分析。实验结果证明了FOAGRNN模型的GNSS高程拟合精度可达6mm。为进一步检验FOAGRNN模型的优越性,采用与平面拟合模型、二次曲面拟合模型进行对比。实验结果表示FOAGRNN模型的拟合精度要优于平面拟合模型和二次曲面拟合模型,证明了FOAGRNN模型在数据样本较少的情况下,其GNSS高程拟合精度仍然可以达到较高精度。