反演地表质量变化的附有方差约束的径向点质量方法

针对目前径向点质量方法在反演GRACE (Gravity Recovery and Climate Experiment)时变地球重力场信号时条带噪声和泄漏误差较大的问题,本文提出了一种附有方差约束的径向点质量方法(Variance Constrained Radial Point-Mass Method,VCPM).该方法利用质量块的先验方差信息设计正则化矩阵对质量变化反演进行空间约束,同时为了解决单一的正则化参数不能适应全球质量变化大小不均匀的问题,构建了迭代的正则化处理策略.为了检验VCPM方法的可靠性,本文反演了2003-01-2014-11全球地表质量变化,并与GRACE球谐位系数法、GLDAS2.1水文模型、CSR (University of Texas Center for Space Research,德克萨斯大学空间研究中心)与JPL (Jet Propulsion Laboratory,喷气推进实验室) RL06 Mascon数据在全球的均方根误差、长期趋势、周年振幅以及局部的冰冻圈、陆地流域质量变化进行了对比分析.结果显示VCPM与CSR Mascon在全球、陆地和海洋区域的RMSE按面积加权平均值分别为2.12 cm、4.16 cm和1.25 cm,在全球和陆地区域均低于CSR Mascon与JPL Mascon的相应RMSE (2.22 cm和4.62 cm),在海洋区域仅略大于CSR Mascon与JPL Mascon的RMSE (1.18 cm);在8个陆地流域和内陆海域的时间序列对比中,VCPM与两种Mascon数据在扣除季节项和趋势项后的相关系数在绝大多数区域仍能达到0.86以上;另外相对于传统点质量方法和球谐位系数法,VCPM能有效抑制南极和格陵兰岛等区域的泄漏误差.这表明本文所提出的VCPM方法可以有效去除条带误差、抑制海陆信号泄漏且大幅提高反演结果的精度和空间分辨率,从而为计算高精度高分辨率地表质量变化提供了一种可供选择的有效途径.

     

Scaling of response spectrum and duration for aftershocks

Aftershocks have the potential to cause collapse of a structure that has been already damaged by the preceding main shock. Seismic safety of a structure should therefore be ascertained through a damage analysis using the anticipated main shock and few larger-aftershock motions. Simulation of aftershock motions needs characterization of the seismic hazard due to aftershocks, and therefore it will be useful to develop a conditional scaling model that can predict the response spectrum of an anticipated aftershock motion consistent with the design spectrum of the main shock motion anticipated at the same station. In this study an attempt is made to develop a conditional scaling model for the pseudo spectral velocity spectrum via linear regression analysis on the aftershock and main shock recordings for the 1999 Chi–Chi earthquake. It is shown that it may be possible to obtain a simpler and approximate version of the conditional model from an unconditional model. Damage-causing potential of a ground motion also depends on its strong motion duration (SMD) and therefore a conditional scaling model is developed for SMD of the aftershock motion in several narrow frequency-bands. The model is developed for the larger-aftershock motions and it is shown that a reasonable replacement of such a model may be obtainable directly from an unconditional model. Finally, a simple weighted averaging scheme is proposed to obtain the composite SMD from the SMDs for different frequency bands by using the pseudo spectral acceleration spectrum of the motion.     

The Data Depth-weight-kernel Estimation of the Model Error of Satellite Orbit Determination

It is an objective fact that there exists error in the satellite dynamic model and it will be transferred to satellite orbit determination algorithm, forming a part of the connotative model error. Mixed with the systematic error and random error of the measurements, they form the unitive model error and badly restrict the precision of the orbit determination. We deduce in detail the equations of orbit improvement for a system with dynamic model error, construct the parametric model for the explicit part of the model and nonparametric model for the error that can not be explicitly described. We also construct the partially linear orbit determination model, estimate and fit the model error using a two-stage estimation and a kernel function estimation, and finally make the corresponding compensation in the orbit determination. Beginning from the data depth theory, a data depth weight kernel estimator for model error is proposed for the sake of promoting the steadiness of model error estimation. Simulation experiments of SBSS are performed. The results show clearly that the model error is one of the most important effects that will influence the precision of the orbit determination. The kernel function method can effectively estimate the model error, with the window width as a major restrict parameter. A data depth-weight-kernel estimation, however, can improve largely the robustness of the kernel function and therefore improve the precision of orbit determination.     

基于离散元方法的落石碰撞破碎研究

落石是山区公路边坡常见的地质灾害。落石撞击破碎是落石灾害中常见的现象,落石破碎致使运动轨迹改变往往会增大落石的危害性,但其通常未被考虑到落石的防护设计中。因此,落石破碎机理的研究对落石防护措施设计有重要的指导意义。本文运用离散元方法（离散元开源软件Esys-Particle）模拟落石铅直撞击地面的过程。在模型中,落石由大量的颗粒组成,并且相邻的颗粒由可断裂的黏结材料黏接;地面由一层固定的颗粒组成,从而模拟地面的摩擦、弹性变形性质。模型模拟出落石反弹、破碎和粉碎等过程,分析黏结材料的杨氏模量、颗粒间凝聚力和内摩擦角对落石撞击地面过程的影响,得出落石撞击过程中破坏率和动能的变化过程。研究发现：当杨氏模量和凝聚力的比值较小时,落石将会发生反弹;当杨氏模量和凝聚力的比值增大到一定程度时,落石将会破碎成不同大小的块体;当杨氏模量和凝聚力的比值超过一定限度时,落石将会破碎成非常细小的块体;内摩擦角的变化对撞击过程影响非常小。     

前方交会法测定建筑物主体倾斜

通过建筑物倾斜观测实例,介绍了应用前方交会法进行倾斜观测的方法、精度要求及注意事项。     

地球重力场与KBR系统频谱关系的建立与分析

以近极圆形轨道的SST ll为研究目标,利用解析的方法分析卫星间精密测距与地球重力场的频谱关系,得出了一般情况下300km高度的SST,星间距离每增加100km,恢复能力提高约10阶次;对于500km卫星高度的SST,星间距离每增加100km,恢复能力提高约5阶次的结论。研究成果可以为我国发展SST探测地球重力场的卫星计划提供参考。     

Robust biased estimation and its applications in geodetic adjustments

In this paper, a number of robust biased estimators (e.g. ordinary robust ridge estimator, robust principal components estimator, robust combined principal components estimator, robust single-parametric principal components estimator, robust root-root estimator) are established by means of a unified expression of biased estimators and based on the principle of equivalent weight. The most attractive advantage of these new estimators is that they can not only overcome the ill-conditioning of the normal equation but also have the ability to resist outliers. A numerical example is used to illustrate that these new estimators are much better than the least-squares estimator and various biased estimators even when both ill-conditioning and outliers exist. Received: 14 November 1995/Accepted: 11 February 1998     

Jason-2高度计数据在中国四大海域的波形重跟踪与分析

在Brown-Hayne模型基础上,利用迭代最小二乘算法实现了回波波形参数的重跟踪,同时利用重心偏移算法(OCOG)、阈值算法对Jason-2 卫星的波形数据在中国渤海、黄海、东海、南海4个海域开展了波形重跟踪试验。试验表明,在4个海域,Brown-Hayne模型法获得的波形重跟踪结果相比较于阈值法、OCOG算法与法国空间局(CNES)公布结果最为一致,相互较差的标准差约2 cm左右。50%阈值算法获得的波形重跟踪结果相比较于其它阈值精度更优,与CNES比较其标准差在3~6 cm,但存在4~10 cm的系统偏差。OCOG算法获得的结果存在较大误差,但其结果可作为Brown-Hayne模型解算的初始值。当卫星地面轨迹靠近大陆时,波形重跟踪获得的海面高改正量会较宽阔海域增大,其中南海区域的海面高改正量最大,可达1 m,其原因可能是卫星地面轨迹经过了西沙群岛,雷达回波波形受到了一定的影响。     

POS与DEM辅助机载SAR多普勒参数估计

由于机载SAR多普勒参数受载机不稳定飞行和复杂地形影响明显,目前仍缺少有效的矢量估计法。本文对POS和DEM数据及其误差对机载多普勒参数的影响进行了分析,以距离-共面方程构建的SAR几何关系为基础,在POS和DEM数据支持下,建立了一种多普勒参数矢量估计法,并通过仿真对提出的方法进行了验证和分析。理论和仿真表明,本文方法在高山地形、不稳飞行及大斜视角等复杂条件下,亦能够为机载SAR瞬时多普勒中心提供高精度估计。     

基于分形方法的蚀变矿物多尺度特征分析

针对矿化蚀变信息常常受其他地物信息的干扰与影像空间分辨率的制约而表现得非常微弱的问题,该文以Landsat7号ETM+30m空间分辨率的影像为数据基础,重采样得到一系列不同空间分辨率的影像;运用ETM+波段3/波段1比值法获得铁染蚀变相对含量灰度图,采用分形模型计算出不同分辨率下蚀变矿物的分维值,进而获得分辨率—分维值的一元线性回归方程;最后建立了基于分形理论的蚀变矿物多尺度分析模型,得到蚀变矿物空间分布特征变化的参考尺度。结果表明,在一定尺度范围内(30~280m),蚀变矿物显示出空间尺度特征。