Precise orbit determination of BeiDou constellation: method comparison

Chinese BeiDou navigation satellite system is in official service as a regional constellation with five geostationary earth orbit (GEO) satellites, five inclined geosynchronous satellite orbit (IGSO) satellites and four medium earth orbit (MEO) satellites. There are mainly two methods for precise orbit determination of the BeiDou constellation found in the current literatures. One is the independent single-system method, where only BeiDou observations are used without help from other GNSS systems. The other is the two-step GPS-assisted method where in the first step, GPS data are used to resolve some common parameters, such as station coordinates, receiver clocks and zenith tropospheric delay parameters, which are then introduced as known quantities in BeiDou processing in the second step. We conduct a thorough performance comparison between the two methods. Observations from the BeiDou experimental tracking stations and the IGS Multi-GNSS Experiment network from January 1 to March 31, 2013, are processed with the Positioning and Navigation Data Analyst (PANDA) software. The results show that for BeiDou IGSO and MEO satellites, the two-step GPS-assisted method outperforms the independent single-system method in both internal orbit overlap precision and external satellite laser ranging validation. For BeiDou GEO satellites, the two methods show close performances. Zenith tropospheric delays estimated from the first method are very close to those estimated from GPS precise point positioning in the second method, with differences of several millimeters. Satellite clock estimates from the two methods show similar performances when assessing the stability of the BeiDou on board clocks.     

实时卫星钟差基准精化方法

在基于精密单点定位(PPP)的授时方法中,卫星钟差产品的高精度时间基准至关重要.针对实时卫星钟差产品时间基准不够稳定的问题,本文采用一组具有原子钟外部输入的国际全球卫星导航系统(GNSS)服务(IGS)跟踪站建立了顾及原子钟变化特性的基准精化方法.该方法首先采用阿伦方差对不同的IGS跟踪站外接原子钟进行稳定度分析,挑选出一组稳定度高的原子钟用以精化时间基准.在此基础上,利用阿伦方差分析各台原子钟的噪声参数特征,并确定不同原子钟之间的权比关系.最终,建立时间基准改正量的随机模型,并计算出精化后的时间基准.通过实例验证表明:与IGS事后精密钟差产品定义的时间基准比较,改正后的实时钟差基准单天内的标准差(STD)优于0.1 ns,相比于改正前最高提升了93%.同时,基准改正后的天内万秒稳达到10^-15量级,实现了一个量级的提高.此外,通过相对钟差精度的分析,表明钟差基准修正不影响PPP的定位精度.     

HHT-EEMD用于IGS站高程时间序列分析

利用经验模态分解（EMD）和整体经验模态分解（EEMD）方法,将BJFS站2000~2015年高程时间序列进行分解,发现其不仅存在1 a、0.5 a、0.25 a、2个月、1个月以及长周期等周期项,还存在以往方法很难探测出来的近似2 a周期信号。与EMD分解结果对比,整体经验模态分解可以明显减弱模式混叠现象。对各分量进行Hilbert 变换（HHT）,得到时间-频率-能量的Hilbert 频谱图。结果表明,年周期和2 a周期变化是高程运动的主要贡献项。利用小波变换方法对比验证EEMD的分解结果表明,与小波分析相比,EEMD重构信号与高程序列差异的RMS更小,证明了HHT-EEMD方法在数据资料分析过程中的有效性。对环境负载及GRACE负载造成的测站位移进行功率谱分析得出,环境负载确实会造成IGS站高程时间序列的1 a、0.5 a以及季节性运动,GRACE负载还验证了2 a信号的存在。     

岩石颗粒破碎的尺寸效应

固体颗粒破碎强度的尺寸效应是一种普遍存在的现象,冰块、岩石颗粒、陶瓷和混凝土块等的破碎强度都表现出随颗粒直径增加而减小的现象,分形模型为解释固体颗粒破碎强度的尺寸效应提供了可行的方法。本文采用Steacy和Sammis分形模型模拟了岩石颗粒压碎特征,分析岩石颗粒破碎后的颗粒分布规律,给出颗粒破碎分维的确定方法,建立颗粒压碎强度与粒径的理论关系,颗粒破碎强度与颗粒粒径的关系用分维D表示为fdD-3。已有的颗粒破碎分布的数据表明,岩石颗粒破碎的分维大约为2.50～2.60,颗粒破碎强度符合用分维表示的尺寸效应。     

温度控制下的煤层气储层物性动态演化研究

为研究热对煤储层的改造作用,采用不同煤级具有代表性的12件煤样。考虑煤级差异、温度作用差异以及核磁共振(NMR)在研究孔隙流体中的技术优势等,对煤样进行了显微组分、镜质体反射率、孔渗测试以及孔隙结构的压汞和NMR实验,研究了煤储层核磁共振特征、T2谱响应特征、孔径分布、核磁共振孔隙度及核磁渗透率、可动流体等,同时对比了不同变质程度煤储层孔隙结构特征及其对煤层气储层物性的控制作用。温度变化下的煤储层物性研究结果表明:热的作用对各煤级煤储层孔隙度都有改善作用,但是对低煤级影响最大,其次为中煤级,最后为高煤级煤储层。温度对不同煤级煤岩的渗透性影响迥异。高挥发份烟煤受热后渗透率增加了一个数量级。低挥发份烟煤受热后渗透率则下降了一个数量级,这应与煤岩热膨胀性有关。无烟煤渗透率总体变化不明显,说明热作用对高煤级煤储层改造意义不大。     

Impact of GLONASS pseudorange inter-channel biases on satellite clock corrections

GLONASS carrier phase and pseudorange observations suffer from inter-channel biases (ICBs) because of frequency division multiple access (FDMA). Therefore, we analyze the effect of GLONASS pseudorange inter-channel biases on the GLONASS clock corrections. Different Analysis Centers (AC) eliminate the impact of GLONASS pseudorange ICBs in different ways. This leads to significant differences in the satellite and AC-specific offsets in the GLONASS clock corrections. Satellite and AC-specific offset differences are strongly correlated with frequency. Furthermore, the GLONASS pseudorange ICBs also leads to day-boundary jumps in the GLONASS clock corrections for the same analysis center between adjacent days. This in turn will influence the accuracy of the combined GPS/GLONASS precise point positioning (PPP) at the day-boundary. To solve these problems, a GNSS clock correction combination method based on the Kalman filter is proposed. During the combination, the AC-specific offsets and the satellite and AC-specific offsets can be estimated. The test results show the feasibility and effectiveness of the proposed clock combination method. The combined clock corrections can effectively weaken the influence of clock day-boundary jumps on combined GPS/GLONASS kinematic PPP. Furthermore, these combined clock corrections can improve the accuracy of the combined GPS/GLONASS static PPP single-day solutions when compared to the accuracy of each analysis center alone.     

Assessment of code bias variations of BDS triple-frequency signals and their impacts on ambiguity resolution for long baselines

GPS Solutions - Carrier phase ambiguity resolution over long baselines is challenging in BDS data processing. This is partially due to the variations of the hardware biases in BDS code signals and...     

An efficient solution of real-time data processing for multi-GNSS network

Global navigation satellite systems (GNSS) are acting as an indispensable tool for geodetic research and global monitoring of the Earth, and they have been rapidly developed over the past few years with abundant GNSS networks, modern constellations, and significant improvement in mathematic models of data processing. However, due to the increasing number of satellites and stations, the computational efficiency becomes a key issue and it could hamper the further development of GNSS applications. In this contribution, this problem is overcome from the aspects of both dense linear algebra algorithms and GNSS processing strategy. First, in order to fully explore the power of modern microprocessors, the square root information filter solution based on the blocked QR factorization employing as many matrix–matrix operations as possible is introduced. In addition, the algorithm complexity of GNSS data processing is further decreased by centralizing the carrier-phase observations and ambiguity parameters, as well as performing the real-time ambiguity resolution and elimination. Based on the QR factorization of the simulated matrix, we can conclude that compared to unblocked QR factorization, the blocked QR factorization can greatly improve processing efficiency with a magnitude of nearly two orders on a personal computer with four 3.30 GHz cores. Then, with 82 globally distributed stations, the processing efficiency is further validated in multi-GNSS (GPS/BDS/Galileo) satellite clock estimation. The results suggest that it will take about 31.38 s per epoch for the unblocked method. While, without any loss of accuracy, it only takes 0.50 and 0.31 s for our new algorithm per epoch for float and fixed clock solutions, respectively.     

Hybrid continuum–discrete simulation of granular impact dynamics

Acta Geotechnica - Granular impact—the dynamic intrusion of solid objects into granular media—is widespread across scientific and engineering applications including geotechnics....