An automated strategy for calculation of phase diagram sections and retrieval of rock properties as a function of physical conditions

We formulate an algorithm for the calculation of stable phase relations of a system with constrained bulk composition as a function of its environmental variables. The basis of this algorithm is the approximate representation of the free energy composition surfaces of solution phases by inscribed polyhedra. This representation leads to discretization of high variance phase fields into a continuous mesh of smaller polygonal fields within which the composition and physical properties of the phases are uniquely determined. The resulting phase diagram sections are useful for understanding the phase relations of complex metamorphic systems and for applications in which it is necessary to establish the variations in rock properties such as density, seismic velocities and volatile‐content through a metamorphic cycle. The algorithm has been implemented within a computer program that is general with respect to both the choice of variables and the number of components and phases possible in a system, and is independent of the structure of the equations of state used to describe the phases of the system.     

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.     

Verification of the dualism ordering method by the canonical correspondence analysis: fish community samples

Two ordination techniques (Romaniszyn diagram, RD, and Canonical correspondence analysis, CCA) were applied to the same data set (fish assemblage composition along two Paraná river tributaries and associated environmental variables), to identify spatial community patterns and further to evaluate the performance of RD. The major patterns in these Brazilian streams were successfully revealed by both methods (i.e. clear separation of both tributaries on the basis of fish assemblage composition and environmental factors), but more complex situations may limit the application of RD in comparison with CCA.     

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.     

Melting of β-quartz up to 2.0 GPa and thermodynamic optimization of the silica liquidus up to 6.0 GPa

Melting relations of β-quartz were experimentally determined at 1.0 GPa (1900±20 °C), 1.5 GPa (2033±20 °C), and 2.0 GPa (2145±20 °C) using a new high-pressure assembly in a piston–cylinder apparatus and substantial differences were found with data previously reported. The new melting data of β-quartz were combined and optimized with all available thermodynamic, volumetric, and phase equilibria data for β-cristobalite, β-quartz and coesite to produce a P–T liquidus diagram for silica valid up to 6.0 GPa. Using the new optimized thermodynamic parameters, the invariant point β-cristobalite+β-quartz+liquid and β-quartz+coesite+liquid were determined to lie at 1687±17 °C and 0.457 GPa, and 2425±25 °C and 5.00 GPa, respectively.     

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.     

由三元系相图计算活度的方法Ⅱ.有二元同分熔点化合物的体系

利用熔化自由能 ,由具有同分熔点二元化合物的三元系相图计算活度。计算的范围遍及三元相图的液相面、液相线、共熔线和三元低共熔点     

GPS在地面沉降监测中的应用探讨

利用GPS进行地面沉降监测能够解决传统测量方法工作量大、周期长的问题。该文结合实例，论述了GPS在地面沉降监测中的应用，分析了GPS测量方法的优点和误差，提出了GPS测量地面沉降的适用条件。     

机械卡盘设计中的计算误差分析与应用

通过对普通机构二维力因素影响的分析，利用修正公式计算求值，以修正单因素影响最大值计算所产生的误差。利用这一计算方法，对CD—3型钻机卡盘进行了计算，得到了正确的设计结果，说明用修正公式对多因素影响采用单因素影响最大值计算进行修正是可行的。