Static deformation of two welded monoclinic elastic half-spaces due to a long inclined strike-slip fault

Static deformation of two monoclinic elastic half-spaces in welded contact due to a long inclined strike-slip fault situated in one of the half-spaces is studied analytically and numerically. Closed-form algebraic expressions for the displacement at any point of the medium are obtained. The variation of the displacement at the interface with the horizontal distance from the fault is studied. The effect of anisotropy on the displacement field is examined. It is found that while the anisotropy of the source half-space has a significant effect on the displacement at the interface, the anisotropy of the other half-space has only a marginal effect.     

Construction of three-phase data to model multiphase flow in porous media: Comparing an optimization approach to the finite element approach

Multiphase flow modelling is a major issue in the assessment of groundwater pollution. Three-phase flows are commonly governed by mathematical models that associate a pressure equation with two saturation equations. These equations involve a number of secondary variables that reflect the fluid behaviour in a porous medium. To improve the computational efficiency of multiphase flow simulators, several simplified reformulations of three-phase flow equations have been proposed. However, they require the construction of new secondary variables adapted to the reformulated flow equations. In this article, two different approaches are compared to quantify these variables. A numerical example is given for a typical fine sand.     

Two formulations for dynamic response of a cylindrical cavity in cross‐anisotropic porous media

Two formulations for calculating dynamic response of a cylindrical cavity in cross‐anisotropic porous media based on complex functions theory are presented. The basis of the method is the solution of Biot's consolidation equations in the complex plane. Employing two groups of potential functions for solid skeleton and pore fluid (each group includes three functions), the u–w formulation of Biot's equations are solved. Difference of these two solutions refers to use of two various potential functions. Equations for calculating stress, displacement and pore pressure fields of the medium are mentioned based on each two formulations. Copyright © 2009 John Wiley & Sons, Ltd.     

非均匀介质中交错网格高阶有限差分数值模拟

地震波场的数值模拟一直是地球物理学的一个重要的研究领域,而在数值正演模拟方法的研究中,计算精度和计算效率是评价该方法有效性及优越性的二个关键问题。这里从一阶速度—应力弹性波动方程出发,着重介绍如何构造离散化模型的网格,如何求解空间导数,如何选取边界条件等内容,从而更有效地提高数值计算的精度与计算效率。文中构造了不同类型的介质模型,并在交错网格中,利用高阶有限差分模拟非均匀介质的波场传播。模拟结果表明,该方法实现简单,具有很好地稳定性和较高的精度,能够直观、高效地反映出介质中波场的传播规律。     

A statistical approach to evaluating distance metrics and analog assignments for pollen records

The modern analog technique typically uses a distance metric to determine the dissimilarity between fossil and modern biological assemblages. Despite this quantitative approach, interpretation of distance metrics is usually qualitative and rules for selection of analogs tend to be ad hoc. We present a statistical tool, the receiver operating characteristic (ROC) curve, which provides a framework for identifying analogs from distance metrics. If modern assemblages are placed into groups (e.g., biomes), this method can (1) evaluate the ability of different distance metrics to distinguish among groups, (2) objectively identify thresholds of the distance metric for determining analogs, and (3) compute a likelihood ratio and a Bayesian probability that a modern group is an analog for an unknown (fossil) assemblage. Applied to a set of 1689 modern pollen assemblages from eastern North America classified into eight biomes, ROC analysis confirmed that the squared-chord distance (SCD) outperforms most other distance metrics. The optimal threshold increased when more dissimilar biomes were compared. The probability of an analog vs no-analog result (a likelihood ratio) increased sharply when SCD decreased below the optimal threshold, indicating a nonlinear relationship between SCD and the probability of analog. Probabilities of analog computed for a postglacial pollen record at Tannersville Bog (Pennsylvania, USA) identified transitions between biomes and periods of no analog.     

Generalized Coarse Graining Procedures for Flow in Porous Media

This paper focuses on heterogeneous soil conductivities and on the impact their resolution has on a solution of the piezometric head equation: owing to spatial variations of the conductivity, the flow properties at larger scales differ from those found for experiments performed at smaller scales. The method of coarse graining is proposed in order to upscale the piezometric head equation on arbitrary intermediate scales. At intermediate scales large scale fluctuations of the conductivities are resolved, whereas small scale fluctuations are smoothed by a partialy spatial filtering procedure. The filtering procedure is performed in Fourier space with the aid of a low-frequency cut-off function. We derive the partially upscaled head equations. In these equations, the impact of the small scale variability is modeled by scale dependent effective conductivities which are determined by additional differential equations. Explicit results for the scale dependent conductivity values are presented in lowest order perturbation theory. The perturbation theory contributions are summed up with using a renormalisation group analysis yielding explicit results for the effective conductivity in isotropic media. Therefore, the results are also valid for highly heterogeneous media. The results are compared with numerical simulations performed by Dykaar and Kitanidis (1992). The method of coarse graining combined by a renormalisation group analysis offers a tool to derive exact and explicit expressions for resolution dependent conductivity values. It is, e.g., relevant for the interpretation of measurement data on different scales and for reduction of grid-block resolution in numerical modeling. This revised version was published online in July 2006 with corrections to the Cover Date.     

纵向弥散作用与渗透介质非均质性定量关系的模拟研究

在地下水污染模拟预报中,弥散参数是很难确定的一个模型参数。因实验室小尺度弥散规律一般不能用于大尺度弥散过程,而野外示踪试验却耗资大、周期长,限制了其实用性。文中利用随机数值模拟手段、基于随机理论的蒙特卡罗方法及序贯高斯模拟技术来生成渗透系数随机场,并研究渗透系数对数场的方差、相关长度以及变异函数类型在不同尺度上对纵向弥散度的影响,进而建立纵向弥散度与随机分布渗透系数场的方差和相关长度的统计定量关系,并与Gelhar理论计算结果进行比较。数值模拟结果表明,经过一定迁移距离后纵向弥散度与随机分布渗透系数对数场的方差和相关长度具有良好的线性统计关系,与Gelhar理论公式表达的关系类型类似。但对于较大的方差,纵向弥散度模拟结果明显大于Gelhar理论计算值,而对于较大相关长度在迁移距离不很大时,纵向弥散度模拟结果明显小于Gelhar理论计算值。本研究可为野外大尺度地下水污染预报模型中水动力弥散参数的确定提供方法借鉴。     

中国土地利用变化现代过程时空特征的研究──基于卫星遥感数据

为实现土地利用变化现代过程的卫星遥感研究,须建立时间序列完整的空间型遥感数据平台。在此基础上,通过建模分析,可以有效地研究土地利用变化的现代过程,其中包括土地利用变化的时间动态特征和空间动态特征。为探讨基于卫星遥感数据的中国土地利用变化现代过程研究的方法论,采用陆地卫星TM分层抽样方法,对全国土地利用类型区域组合、土地利用程度变化和土地利用动态度进行了分析。     

不连续面在双重介质热-水-力三维耦合分析中的有限元数值实现

在建立双重介质热-水-力耦合微分控制方程的基础上,提出了裂隙岩体热-水-力耦合的三维力学模型,对不同介质分别建立以节点位移、水压力和温度为求解量的三维有限元格式,开发了双重介质热-水-力耦合分析的的三维有限元计算程序,在有限元数值分析中不连续面应力计算采用等厚度空间8节点节理单元进行离散,而不连续面渗流和热能计算时采用平面4节点等参单元进行离散,这样保证了不同介质之间的水量、热量交换和两类模型接触处节点水头、温度和位移相等。通过高温岩体地热开发算例,揭示了在热-水-力耦合作用下不连续面处于低应力区,其张开度随运行时间的延长呈非线性增加,非稳定渗流阶段不连续面显著地控制着渗流场的整体分布,它的水头远高于拟连续岩体介质的水头,而进入稳定渗流阶段不连续面的控渗作用不明显,由于高温岩体地热开发系统中存在大规模的热量补给,不连续面对岩体温度场分布的影响并不显著。     

湘黔地区下寒武统黑色岩系中镍-钼矿床黄铁矿的成因

重点研究了贵州遵义中南村和湖南张家界三岔镍-钼多金属矿床的黄铁矿,矿床成因为热水喷流沉积型.黄铁矿样品S/Fe比值的平均值为2.0093-2.048,成分特征均属于铁亏损型,说明其形成温度低.Co/Ni变化范围为0.077~4.5,均值为2.197.据Co/Ni的比值可以判断本区的黄铁矿主要为热液成因,矿化物质来源为热水喷流体系内的热液.成分图解标型显示,黑色岩系中高ω(Co+Ni)/ω(Fe)的黄铁矿和高ω(As+Se+Te)/ω(S)的黄铁矿均大量出现,可能与热水喷流沉积的温差变化范围较大有关.而并非正常沉积的产物.S/Se比值为950.8-1059.9,说明生成环境均一程度较高,而且海水温度比较高.由中南村镍-钼矿床的Se/Te值也能判明其为热液成因.关于黄铁矿的研究对本区镍-钼矿床的寻找有重要意义.