Compositional Geometry and Mass Conservation

A geometrical structure is imposed on compositional data by physical and chemical laws, principally mass conservation. Therefore, statistical or mathematical investigation of possible relations between data values and such laws must be consistent with this structure. This demands that geometrical concepts, such as points that specify both mass and composition in linear space, and lines in projective space that specify composition only, be clearly defined and consistent with mass conservation. Mass thus becomes the norm in composition space in place of the Euclidean norm of ordinary space. Coordinate transformations inconsistent with this geometry are accordingly unnatural and misleading. They are also unnecessary because correlation arising from the constant mass presents no unusual difficulty in the analysis of the underlying quadratic form.     

高斯平面坐标换算到测区平均高程面上的方法

依据高斯投影原理，提出了一种在高斯平面与任决高程平面间进行坐标转换的方法，并设计了相应的计算软件。     

Conditional curvilinear stochastic simulation using pixel-based algorithms

In geology, structures displaying differing local directions of continuity are widespread, a typical example being a flusial depositional system. Conventional pixel-based geostatistical algorithms, may fail to reproduce such curvilinear structures. Conversely, object-based algorithms can reproduce curvilinear shapes but are difficult to condition to dense local data. Local depositional directions as obtained from dipmeter data. 3D seismic data, and geological interpretation represent critical information. An improved pixel-based geostatistical algorithm is proposed to account for such directional information. Case studies demonstrate the potential and limitations of the algorithm.     

论广义逆矩阵阶数增大的递推算法

本文针对加拿大地球物理学家ＥｄｏＮｇｌａｎｄ教授提出的问题，对Ｍ—Ｐ广义逆矩阵的递推算法进行了探分论述探讨的初步结果和对Ｇｒｅｖｉｌｌｅ递推算法的扩充。文中对算法的扩充给出严格的证明。     

三维溶质运移问题的分步广义迎风解法

对对流占优的三维溶质运移问题提出了分步广义迎风解法，首先利用Ｎ．Ｎ．Ｙａ－ｎｅｎｋｏ对水动力弥散方程分步求解的思想，将原来的一个定解问题分解为两个定解问题即对流定解问题和扩散定解问题，对对流定解问题采用广义迎风对偶单元均衡法求解，对扩散定解问题采用一般的Ｇａｌｅｒｋｉｎ有限元法求解，不仅避免了用一般有限元法和有限差分法求解对流占优的地下水水质数学模型时常出现数值弥散和过量问题，而且避免了求节点速度这一步，节省运算步骤，对井点的浓度变化给出了更合适的求解方法。     

Generalized plasticity state parameter‐based model for saturated and unsaturated soils. Part 1: Saturated state

The behavior of granular materials is known to depend on its loose or dense nature, which in turns depends both on density and confining pressure. Many models developed in the past require the use of different sets of constitutive parameters for the same material under different confining pressures. The purpose of this paper is to extend a basic generalized plasticity model for sands proposed by Pastor, Zienkiewicz and Chan by modifying the main ingredients of the model flow—rule, loading–unloading discriminating direction and plastic modulus—to include a dependency on the state parameter. The proposed model is tested against the available experimental data on three different sands, using for each of them a single set of material parameters, finding a reasonably good agreement between experiments and predictions. Copyright © 2010 John Wiley & Sons, Ltd.     

Generalized plasticity state parameter‐based model for saturated and unsaturated soils. Part II: Unsaturated soil modeling

The aim of this paper is to extend the generalized plasticity state parameter‐based model presented in part 1 to reproduce the hydro‐mechanical behavior of unsaturated soils. The proposed model is based on two pairs of stress–strain variables and a suitable hardening law taking into account the bonding—debonding effect of suction and degree of saturation. A generalized state parameter for unsaturated state is proposed to reproduce soil behavior using a single set of material parameters. Generalized plasticity gives a suitable framework to reproduce not only monotonic stress path but also cyclic behavior. The hydraulic hysteresis during a drying—wetting cycle and the void ratio effect on the hydraulic behavior is introduced. Comparison between model simulations and a series of experimental data available, both cohesive and granular, are given to illustrate the accuracy of the enhanced generalized plasticity equation. Copyright © 2010 John Wiley & Sons, Ltd.     

A nonstationary extreme value distribution for analysing the cessation of karst spring discharge

The effects of climate change and population growth in recent decades are leading us to consider their combined and potentially extreme consequences, particularly regarding hydrological processes, which can be modeled using a generalized extreme value (GEV) distribution. Most of the GEV models were based on a stationary assumption for hydrological processes, in contrast to the nonstationary reality due to climate change and human activities. In this paper, we present the nonstationary generalized extreme value (NSGEV) distribution and use it to investigate the risk of Niangziguan Springs discharge decreasing to zero. Rather than assuming the location, scale, and shape parameters to be constant as one might do for a stationary GEV distribution analysis, the NSGEV approach can reflect the dynamic processes by defining the GEV parameters as functions of time. Because most of the GEV model is designed to evaluate maxima (e.g. flooding, represented by positive numbers), and spring discharge cessation is a ?minima’, we deduced an NSGEV model for minima by applying opposite numbers, i.e. negative instead of positive numbers. The results of the model application to Niangziguan Springs showed that the probability of zero discharge at Niangziguan Springs will be 1/80 in 2025, and 1/10 in 2030. After 2025, the rate of decrease in spring discharge will accelerate, and the probability that Niangziguan Springs will cease flowing will dramatically increase. The NSGEV model is a robust method for analysing karst spring discharge. Copyright © 2013 John Wiley & Sons, Ltd.     

Improved algorithms for generalized finite element simulations of three-dimensional hydraulic fracture propagation

This paper reports improvements to algorithms for the simulation of 3-D hydraulic fracturing with the Generalized Finite Element Method (GFEM). Three optimizations are presented and analyzed. First, an improved initial guess based on solving a 3-D elastic problem with the pressure from the previous step is shown to decrease the number of Newton iterations and increase robustness. Second, an improved methodology to find the time step that leads to fracture propagation is proposed and shown to decrease significantly the number of iterations. Third, reduced computational cost is observed by properly recycling the linear part of the coupled stiffness matrix. Two representative examples are used to analyze these improvements. Additionally, a methodology to include the leak-off term is presented and verified against asymptotic analytical solutions. Conservation of mass is shown to be well satisfied in all examples.     

南沙群岛海域长腹剑水蚤(Oithona spp.)的种类组成、数量分布及其与环境因子的关系

中小型桡足类是海洋食物网中重要的中间环节,在碳循环里细菌碳和藻类碳向高营养级转化的过程中扮演重要角色。本文通过分析2013年3—4月、6—7月、9—10月和11—12月南沙群岛海域(111°30′E—112°30′E,4°30′N—11°30′N)的生态环境调查数据,对长腹剑水蚤属(Oithona)的种类组成、数量分布及其与环境因子的关系进行了研究。结果表明:(1)南沙群岛海域共出现长腹剑水蚤15种(包括1个未定种),其中线长腹水蚤(Oithona linearis)、粗长腹剑水蚤(O.robusta)和长刺长腹剑水蚤(O.longispina)为南沙群岛海域首次记录,长腹剑水蚤属种类组成沿水深变化的现象较为明显,该属虽然在南沙群岛海域广泛分布,但单种的出现频率较低,其种类出现的季节性较强;(2)南沙群岛海域的长腹剑水蚤可以划分为常见广布组群Ⅰ和低频组群Ⅱ;(3)长腹剑水蚤年均栖息密度为30ind./m3,季节变化呈单峰型,秋季最高、冬季最低,高数量区主要出现在受南沙西部沿岸流和东部沿岸流影响的近岸海域;(4)在0—750m水深范围内,长腹剑水蚤数量沿水深梯度的垂直变化明显,呈现由浅至深数量明显减少的趋势,长腹剑水蚤主要集中分布在0—75m的水层中,75—500m数量急剧降低至最低后,在500—750m基本维持不变;(5)GAM分析表明,纬度、海水温度、盐度和叶绿素a浓度对长腹剑水蚤数量均有影响,以纬度的影响最为显著,其中长腹剑水蚤属最适温度为28.6—29.2℃,最适盐度为32.6—33.2,最适叶绿素a浓度为5—10μg/L;(6)长腹剑水蚤属更适栖息于低温和低盐,且受沿岸流影响的环境中。因此,其数量分布的季节变化与季风驱动的沿岸流势力强弱密切相关。