针对系统噪声及测量噪声统计量不准确的组合导航系统线性高斯状态模型,提出一种组合导航系统的改进变分贝叶斯自适应滤波算法(modified variational Bayesian adaptive filter,MVBAKF)。首先,选择Wishart分布作为已知均值的高斯分布协方差矩阵的共轭先验,并给出测量噪声方差、状态向量及其预测误差协方差矩阵的联合概率分布函数;然后,利用变分贝叶斯方法给出测量噪声方差及状态向量预测误差协方差矩阵的计算公式,进而提出具有迭代性质的MVBAKF算法;最后,进行基于MVBAKF算法的GNSS/SINS组合导航系统仿真实验。结果表明,相对于传统VBAKF算法,MVBAKF算法可较准确地估计测量噪声方差,有效克服系统噪声统计量不准确对滤波精度的影响,进而提高组合导航系统的滤波精度。
考虑到高性能原子钟的高稳定性及对应站钟钟差序列历元间强相关等特点,对接收机钟差附加先验信息约束进行精密单点定位授时测试。基于BDS-3的B1C/B2b观测信息测试结果表明,相比于传统授时模型,附加先验信息的精密授时模型具有更高的授时性能,尤其是短期稳定性提升效果可达60%。此外,BDS-3的B1I/B3I和B1C/B2b组合观测值PPP精度基本相当,均可获得水平方向优于0.5 cm、高程方向约2 cm的定位精度;得益于B1C/B2b更高的观测数据质量,其定位精度较B1I/B3I组合提高13.8%。
为对比北斗三号卫星导航系统新老信号的定位性能,利用连续两周17个MGEX测站观测数据对BDS-3新频点B1C/B2a消电离层组合精密单点定位(PPP)服务进行全面评估。结果表明,单系统BDS-3 B1C/B2a静态PPP在E、N、U三方向的精度分别为1.3 cm、0.8 cm、1.9 cm,仿动态E、N、U三方向的精度分别为2.6 cm、1.8 cm、4.2 cm;与BDS-3 B1I/B3I组合相比,静态E、N两方向精度相当,U方向提升13.6%,仿动态三方向精度相当,差值在1 mm以内。BDS-3 B1C/B2a 静态PPP三维精度40 min可达cm级,仿动态PPP三维精度5 min达到dm级。双系统BDS-3 B1C/B2a+GPS L1/L2静态PPP三维精度20 min达到cm级,仿动态PPP三维精度30 min达到cm级;与双系统BDS-3 B1I/B3I+GPS L1/L2对比,两者定位精度和收敛速度基本一致。综上所述,可以证明B1C/B2a信号的定位性能优于B1I/B3I信号。
收集低成本接收机3种速度区间路段的动态数据,分析多系统互操作信号B1C/L1/E1和B2a/L5/E5a在不同速度区间的数据接收能力、数据质量和导航定位性能。结果表明:1)随着速度加快,互操作信号的数据接收能力逐渐变差,其中BDS-3可以保持较多的可见卫星数,平均在9颗以上;B1C/L1/E1数据完整率较优,但更容易发生数据中断。2)随着速度加快,互操作信号载噪比降低、周跳比增加,而伪距与载波相位噪声及多路径效应受速度影响较弱,其中B1C/L1/E1载波相位噪声与多路径效应值较小。3)载体速度对定位精度的影响没有明显的规律,速度与环境遮挡程度综合影响定位精度,其中多系统互操作信号组合定位较单系统定位性能提升较多;BDS-3在中低速度区间精度较高,PPP精度较其他系统高17%以上。
MODIS反演的可降水量(PWV)为连续的面状数据,精度不高;而GNSS水汽探测技术精度高,但为离散的点数据。本文综合两者优势,提出一种融合反距离加权和傅里叶变换的MODIS水汽校正方法,利用GNSS PWV来校正MODIS PWV产品,并采用香港地区GNSS数据和MODIS水汽产品进行实验验证。该方法首先考虑空间插值问题,使GNSS PWV与MODIS水汽在空间点位上一致。选取不同数量的MODIS像元点,对比6种插值算法的精度,发现选取5个插值参数时插值效果最佳,选择反距离加权法得到的插值水汽产品精度最高,平均偏差为-0.99 mm。其次考虑水汽校正问题,构建3种模型对MODIS PWV进行校正,发现傅里叶变换模型的均方根误差改进率可达70%以上。新模型的MODIS水汽校正方法能够减弱局部数据的影响,更准确地反映该地区的水汽分布情况。
分别运用原始曲线分析法和归一化速率法分析2021-05-22玛多MS7.4地震震中400 km范围内4个地电阻率台站的数据变化。结果表明:1)玛多地震前,4个地电阻率台站资料出现不同程度的异常变化,其中玛曲NW、金银滩EW、甘孜N10°E、拦隆口EW测道原始曲线变幅最明显,表现为小时值突跳变化、日均值加速变化和曲线转折变化、月均值破年变和保留年变形态的趋势变化;2)玛曲NW、金银滩EW、甘孜N10°E、拦隆口EW归一化速率曲线在玛多地震前均出现超阈值现象,但各自归一化速率形态不同;3)综合分析认为,4个台站地电阻率资料异常变化与玛多MS7.4地震的孕育及活动有关。
采用最小二乘与迭代算法,对A10-057绝对重力仪、FG5X-265绝对重力仪及iGrav-053超导重力仪的同期同址观测数据进行处理,得到FG5X-265及A10-057精密测定iGrav-053的格值因子分别为-883.929 4 ns-2/V、-884.134 ns-2/V,标定精度为0.048 26%、0.093 7%,均优于0.1%。同时对比FG5X-265与A10-057在标定工作中各项测量指标与选取测量对数对标定结果的影响,通过列举多场景对格值精度的要求,最终证明A10型绝对重力仪具备标定相对重力仪的能力,为后续相对重力仪标定工作提供参考。
利用滇西北地区(25.0°~26.9°N,99.20°~100.74°E)85个地震台站记录的三分量远震波形资料提取P波接收函数,并选用H-k扫描法获取研究区莫霍面深度和泊松比,再根据接收函数共转换点(CCP)叠加方法绘制NW向和NE向共6条垂直剖面。结果表明,研究区地壳厚度在38.5~52.4 km之间变化,呈纵向差异分布;泊松比范围为0.24~0.34,多呈四象限分布特征;2013年以来云龙、洱源、漾濞等几次5级以上地震均发生在泊松比高梯度带上。从CCP剖面来看,研究区NW向莫霍面起伏相对较大,由北向南莫霍面呈变浅趋势,上隆幅度约为5~10 km;NE向莫霍面起伏相对平稳,平均深度在42 km左右;地壳内分层情况较为复杂,几乎没有较为连续的界面。漾濞地震下方中上地壳存在较宽的负幅值带,推测可能存在低速区并为地震的发生提供了孕震环境。D410间断面平均深度约为420 km,D660间断面平均深度接近660 km,地幔转换带的厚度在210~260 km之间。
为给出湖北襄阳地区典型粘性土、粉土、砂土及碎石土剪切波波速与深度之间的统计公式,基于收集到的102个实测钻孔波速资料,利用多种一元回归模型对襄阳地区常见土层剪切波波速与深度之间的关系进行拟合分析和相关性研究;同时,讨论土体状态对二者关系的影响;最后,利用实测钻孔数据验证统计公式的合理性和可靠性。研究结果表明:1)除人工填土外,襄阳地区其他常规土类剪切波波速与深度之间存在较强的相关性,且离散程度随深度的增加而增大;2)同一土类随密实度或塑性状态的增大,其平均剪切波波速也变大;3)考虑土体状态可以使回归模型的预测结果更加具体、更接近真实值,但对模型拟合优度的变化具有不确定性;4)区域性对土层剪切波波速与深度的统计关系有较大影响。本文研究的襄阳地区各土类剪切波波速回归模型较全国土类模型具有更好的预测精度和误差平稳性,可为南襄盆地及周边区域地震灾害风险普查、场地土动力学性质研究和地震动参数确定等提供参考。
构建三峡地区大气和陆地水质量变化模型,利用GNSS、GRACE-FO和环境质量变化模型计算三峡地区2019~2022年环境负荷引起的地壳垂直形变,探讨基于GRACE-FO数据的球谐函数法和结合质量变化模型的格林函数法在计算环境负荷垂直形变上的差异。研究发现,2种方法计算的陆地水负荷垂直形变在时空分布上具有一定的一致性,但是后者的精度和空间分辨率更高;基于质量变化模型计算的CQWZ站和HBZG站的环境总负荷垂直形变与GNSS实测结果相关系数达到0.88和0.83;2019~2022年三峡库区和洞庭湖流域的地壳垂直形变速率分别为0.5~0.7 mm/a和0.6~1.5 mm/a,原因是2022年该区域水位降低导致地壳回弹上升;三峡库区附近地壳夏季上升、冬季下沉,垂直形变幅度约为-5~18 mm,主要受水库蓄水和大气负荷驱动的影响。
利用SAR偏移量追踪(offset-tracking)技术获取冰川形变作为SBAS-InSAR的补充,采用2种技术计算2018-01~10色东普流域灾前形变,联合分析灾前形变特征及影响因素。结果表明,色东普流域冰川与沟道在2018-10-17冰崩灾害发生前已出现形变;冰川主要形变区形变趋势表现为加速-平缓-加速,7~9月形变量达到-7.69 m;沟道内堆积物长期呈下滑趋势,7月后与冰川均加速形变;气温升高是冰崩碎屑流灾害发生的主导因素。联合SBAS-InSAR与offset-tracking技术能够满足不同形变量级的监测需求,可用于冰崩灾害的早期识别与形变反演,为青藏高原地区冰崩灾害防治提供参考。
利用StaMPS技术,基于2017~2021年覆盖北京首都国际机场的133景升轨和95景降轨Sentinel-1影像,对机场及周边地面沉降、机场内地裂缝稳定性等进行研究。结果显示,升降轨影像监测结果的一致性较好,在机场南部、温榆河以南区域形成沉降漏斗,漏斗内部最大沉降速率超过40 mm/a,沉降漏斗沿温榆河分布。对升降轨影像监测结果进行内符合精度验证,结果表明,90%以上公共像元的形变速率差位于±10 mm/a之间,证明了升降轨影像监测结果的一致性。特征点时序形变结果表明,金盏乡、富豪村和黄金花园小区附近沉降较严重,最大累积沉降量分别超过250 mm、150 mm和110 mm。从机场整体形变及地裂缝两侧时序和剖线结果均可看出,地裂缝两侧地表的差异形变对比明显,且差异分界线与地裂缝位置较吻合。
VP倾斜仪固体潮信号受仪器监测复杂环境限制,多含有大量环境噪声。为获得真实固体潮曲线,提出一种基于灰色关联分析优化改进的自适应噪声完备集合经验模态分解(ICEEMDAN)VP倾斜仪信号降噪模型(GRA-ICEEMDAN)。该方法首先将含噪信号进行ICCEMDAN处理,得到若干个固有模态函数(IMF),并依次排列与标记;然后基于这些IMF分别计算相关系数、互信息、R2、Adj-R2、MSE、SSE、RMSE、MAE、MAPE、样本熵等10个评价指标值,构建IMF可信度评价指标矩阵;最后借助灰色关联分析(GRA)计算各评价指标与不同IMF之间的关联系数和关联度,依据关联度大小对各个IMF进行排序,将排名靠前的IMF进行线性重构,即可完成信号降噪。仿真去噪实验和实测去噪实验均表明,GRA-ICEEMDAN模型优于卡尔曼滤波、70阶低通FIR滤波、Savitzky-Golay等经典降噪模型,能显著区分噪声成分和有效成分,原始信号分解后的重构误差与信号损失极小,可推广至其他仪器的复杂信号降噪中。