GNSS系统时间偏差参数分析

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

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

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

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

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

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

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的西秦岭北缘断裂带现今活动特征研究

基于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向。震前基线均有挤压缩短的前兆异常显示。梳理站点异常出现情况,发现异常持续时间越长、异常出现次数越多,发震震级可能越大。

融合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观测时间序列中信号提取与模型参数估计提供了一个有价值的模型。

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海洋反射信号接收机设计与分析

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

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

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

一种基于改进免疫算法的GNSS快速选星方法

针对多GNSS系统组合导航选星的实时性和可靠性需求,提出一种基于改进免疫算法的快速选星方法。该方法通过选取最优的迭代次数和初始化种群数目,可以在获取理想卫星几何结构的同时提高接收机运算效率。利用CUT0测站1 d实测数据对所提算法进行验证,结果表明,当初始化种群数目选在20附近、最大迭代次数选在120附近时,选星效果最好;在单BDS系统定位时,平均GDOP差值与遍历法的比值仅为0.57%,耗时比遍历法缩短53.27%;在GPS/BDS/Galileo/GLONASS/QZSS组合系统下,其平均GDOP与遍历法相差仅0.122,平均耗时(1.06 s)比遍历法缩短96.47%。该算法适用于多GNSS系统组合定位,具有良好的普适性。

基于GNSS载波相位的守时系统远程性能评估方法研究

为测试在建守时系统性能,提出基于GNSS载波相位的远程评估方法,并开展实验验证。以国家标准时间UTC(NTSC)为参考,基于西安和北京地区4个测站连续15 d的GNSS观测数据,分别采用PPP和载波相位共视(PCV)法解算得到待测守时系统相对UTC(NTSC)的时间偏差,然后对待测守时系统的相对频率偏差、频率稳定度等关键指标进行分析。结果表明,PCV和PPP时间传递链路的系统误差小,时间传递准确度高,采用PCV和PPP时间传递可准确测量相距较远的守时系统的频率稳定度、相对频率偏差等特性,可为下一步国产守时系统性能远程测试打下坚实基础。

基于GNSS的红河断裂中段现今变形特征

采用红河断裂带中段加密布设的GNSS连续站观测资料以及前人发布的速度场结果,获取红河断裂带中段及邻区现今地壳运动速度场。根据川滇地区活动断裂分布建立研究区三维有限元模型,以红河断裂带中段及邻区GNSS速度场为约束,获得红河断裂带中段不同段落的现今滑动速率和区域应变率场。结果表明,红河断裂带弥渡-元江段右旋走滑速率为1.2±0.6 mm/a,挤压速率为0.6±0.5 mm/a;红河断裂带元江-元阳段右旋走滑速率为1.8±0.7 mm/a,挤压速率为1.5±0.6 mm/a。应变率结果显示,红河断裂带中段及邻区以剪切变形为主,最大剪应变率高值区位于小江断裂带附近,最大幅度约为62×10^(-9)/a,红河断裂剪切变形相对较弱;面应变率显示,红河断裂元江-元阳段挤压变形较为显著,挤压应变率值约为10×10^(-9)/a,该段落处于强闭锁状态,未来地震危险性值得关注。

基于GNSS观测的陆地水储量研究进展

介绍了基于质量负荷理论利用GNSS观测反演陆地水储量的方法,阐述了利用GNSS研究陆地水储量的进展、应用和发展趋势,以及GNSS与GRACE、水文模型等的联合研究。随着GNSS观测精度的不断提升和观测网络的迅猛发展,利用GNSS观测研究陆地水储量取得了重要进展。GNSS已成为陆地水储量监测的重要手段,在中国大陆地区利用GNSS研究陆地水储量具有巨大的潜力和前景,但仍需要在时间序列信号精细识别、反演算法优化和多源数据融合等方面深入研究。

贵州喀斯特地貌地区GNSS站坐标时间序列特性研究

使用GAMIT/GLOBK10.71软件解算贵州省北斗卫星导航定位基准站网(GZCORS)25个基准站数据,获得各站单天解坐标时间序列。利用主成分分析法提取站坐标共模误差,利用谱分析法分析共模误差基本特性,采用极大似然估计法确定站点最优噪声模型及其运动速度。结果表明,GZCORS坐标序列共模误差中包含有周期项,N方向最大振幅周期分别出现在0.2周/a、1.2周/a、3.2周/a和4.2周/a;GZCORS站点最优噪声模型以WN+FN和WN+GM为主,剔除共模误差后,36%的测站分量噪声特性发生变化;剔除共模误差后,坐标时间序列噪声水平明显降低,各坐标分量速度参数的估计精度均有明显提升,其中N、E、U分量分别提高52%、56%和50%。