附加系统参数平差法在导线网中的应用

本文探讨了附加系统参数平差法在带有测距系统误差的导线网中的应用，说明了当系统误差比较显著时．采用附加系统参数平差法的可行性、必要性及其优越性。     

Multi-reference evaluation of uncertainty in earth orientation parameter measurements

Uncertainties in polar motion and length-of-day measurements are evaluated empirically using several data series from the space-geodetic techniques of the global positioning system (GPS), satellite laser ranging (SLR), and very long baseline interferometry (VLBI) during 1997–2002. In the evaluation procedure employed here, known as the three-corner hat (TCH) technique, the signal common to each series is eliminated by forming pair-wise differences between the series, thus requiring no assumed values for the “truth” signal. From the variances of the differenced series, the uncertainty of each series can be recovered when reasonable assumptions are made about the correlations between the series. In order to form the pair-wise differences, the series data must be given at the same epoch. All measurement data sets studied here were sampled at noon (UTC); except for the VLBI series, whose data are interpolated to noon and whose UT1 values are also numerically differentiated to obtain LOD. The numerical error introduced to the VLBI values by the interpolation and differentiation is shown to be comparable in magnitude to the values determined by the TCH method for the uncertainties of the VLBI series. The TCH estimates for the VLBI series are corrupted by such numerical errors mostly as a result of the relatively large data intervals. Of the remaining data sets studied here, it is found that the IGS Final combined series has the smallest polar motion and length-of-day uncertainties.     

A new method for indentifying the model error of adjustment system

Some theory problems affecting parameter estimation are discussed in this paper. Influence and transformation between errors of stochastic and functional models is pointed out as well. For choosing the best adjustment model, a formula, which is different from the literatures existing methods, for estimating and identifying the model error, is proposed. On the basis of the proposed formula, an effective approach of selecting the best model of adjustment system is given. Project supported by the Open Research Fund Program of the Key Laboratory of Geospace Environment and Geodesy, Ministry of Education, Wuhan University (No. 905276031-04-10).     

Investigating the propagation mechanism of unmodelled systematic errors on coordinate time series estimated using least squares

The propagation of unmodelled systematic errors into coordinate time series computed using least squares is investigated, to improve the understanding of unexplained signals and apparent noise in geodetic (especially GPS) coordinate time series. Such coordinate time series are invariably based on a functional model linearised using only zero and first-order terms of a (Taylor) series expansion about the approximate coordinates of the unknown point. The effect of such truncation errors is investigated through the derivation of a generalised systematic error model for the simple case of range observations from a single known reference point to a point which is assumed to be at rest by the least squares model but is in fact in motion. The systematic error function for a one pseudo-satellite two-dimensional case, designed to be as simple but as analogous to GPS positioning as possible, is quantified. It is shown that the combination of a moving reference point and unmodelled periodic displacement at the unknown point of interest, due to ocean tide loading, for example, results in an output coordinate time series containing many periodic terms when only zero and first-order expansion terms are used in the linearisation of the functional model. The amplitude, phase and period of these terms is dependent on the input amplitude, the locations of the unknown point and reference point, and the period of the reference point's motion. The dominant output signals that arise due to truncation errors match those found in coordinate time series obtained from both simulated data and real three-dimensional GPS data.     

显式自由表面模型时间分裂格式的一个积分方案

讨论了有关显式自由表面模型的时间分裂格式中的分裂误差。在传统的积分方案中，由于采用蛙跳格式，斜压动量方程具有三层的时间离散结构，但是，正压动量方程在相应的时间层上只有二层的离散结构。因此，正压方程与斜压方程的时间离散结构是不相容的。分析结果表明，这是导致分裂误差的原因。本提出一种能保持时间离散结构相容的正压和斜压方程的时间积分方案。该积分方案的分裂误差很小因此特别适合于需要长期积分的气候研究。该方案在北太平洋的一个自由表面模型中得到了应用。     

Numerical Analysis of Space–Time Finite Element Formulations for Miscible Displacements

A finite element formulation is proposed to approximate a nonlinear system of partial differential equations, composed by an elliptic subsystem for the pressure–velocity and a transport equation (convection–diffusion) for the concentration, which models the incompressible miscible displacement of one fluid by another in a rigid porous media. The pressure is approximated by the classical Galerkin method and the velocity is calculated by a post-processing technique. Then, the concentration is obtained by a Galerkin/least-squares space–time (GLS/ST) finite element method. A numerical analysis is developed for the concentration approximation. Then, stability, convergence and numerical results are presented confirming the a priori error estimates.     

ARTIFICIAL NEURAL NETWORK CORRECTION OF ROUGH ERRORS OF OBSERVATIONS*

In the context of tower measured radiation datasets.following the correction principle meeting a diagnostic equation in data quality control and in terms of a technique for model construction on data and ANN (artificial neural network) retrieval for BP correction of radiation measurements with rough errors available,a BP model is presented.Evidence suggests that the developed model works well and is superior to a convenient multivariate linear regression model,indicating its wide applications.     

饱和岩石的各向异性及非线性黏弹性响应

在MTS伺服压机上对不同饱和状态的砂岩、大理岩标本进行了垂直层理和平行层理两个方向的正弦波加载试验,研究了饱和岩石的各向异性及非线性黏弹性行为. 在最大载荷低于岩石的屈服点时,获得的衰减、杨氏模量、泊松比、波速都表现出显著的各向异性,同时它们还具有较明显的应变振幅效应和频率效应. 随着应变振幅的增大,衰减近似呈线性增长,杨氏模量和泊松比近似呈线性下降. 在0.005-4Hz频段内除衰减不明显依赖频率以外,杨氏模量、泊松比、波速对频率的依赖性都较强; 当频段扩大到15Hz时,这4个参数都表现出较强的频率效应. 饱和岩石的衰减、杨氏模量、泊松比、波速均随饱和液体的黏滞系数增大而增大.     

Precise revision of the garnet-muscovite geothermometer

The garnet-muscovite geothermometer was refined through empirical calibration by using natural rocks metamorphosed under the physical conditions of 238—1306 MPa and 490—700℃. Input temperatures and pressures were determined through simultaneously applying the garnet-biotite geothermometer and the garnet-biotite-plagioclase-quartz barometer, assuming that all FeO in muscovite and garnet be ferrous. Garnet was treated as the asymmetric quaternary solid solution, and muscovite as the symmetric binary solid solution. Input muscovite compositions include Fe atoms between 0.03—0.19 and Mg atoms between 0.04—0.16 on the basis of 11 oxygen atoms, and input garnet compositions include spessartine fractions between 0.01—0.289, grossular fractions between 0.028—0.273, and the Fe/Mg ratio between 3.387-18.986. The resulting garnet-muscovite geothermometer reproduces temperatures within (50℃ compared with the garnet-biotite thermometer. Total random error of ±37℃ of the new thermometer may stem from the pressure uncertainty of ±200 MPa, and uncertainties of ±5% of Fe and Mg components in muscovite, and ±5% of Fe, Mg, Mn and Ca components in garnet, altogether. When there exist 10%, 20%, 30%, 40% and 50% Fe3+ in muscovite, respectively, the computed garnet-muscovite temperatures will be 1—6℃, 2—12℃, 3—16℃, 5—24℃ and 7—29℃, respectively, lower than those obtained when assuming that all FeO be ferrous. The new garnet-muscovite geothermometer can efficiently reflect temperature change of typical prograde sequences and contact aureole rocks, and may be applied to low- to high-grade and low- to high-pressure metamorphic rocks.