Three dimensional simulation of fluid flow in X‐ray CT images of porous media

A numerical scheme is developed in order to simulate fluid flow in three dimensional (3‐D) microstructures. The governing equations for steady incompressible flow are solved using the semi‐implicit method for pressure‐linked equations (SIMPLE) finite difference scheme within a non‐staggered grid system that represents the 3‐D microstructure. This system allows solving the governing equations using only one computational cell. The numerical scheme is verified through simulating fluid flow in idealized 3‐D microstructures with known closed form solutions for permeability. The numerical factors affecting the solution in terms of convergence and accuracy are also discussed. These factors include the resolution of the analysed microstructure and the truncation criterion. Fluid flow in 2‐D X‐ray computed tomography (CT) images of real porous media microstructure is also simulated using this numerical model. These real microstructures include field cores of asphalt mixes, laboratory linear kneading compactor (LKC) specimens, and laboratory Superpave gyratory compactor (SGC) specimens. The numerical results for the permeability of the real microstructures are compared with the results from closed form solutions. Copyright © 2004 John Wiley & Sons, Ltd.     

Isotope hydrological study of mean transit time in the granitic Strengbach catchment (Vosges massif,France): application of the FlowPC model with modified input function

Measurements of ¹⁸O concentrations in precipitation, soil solution, spring and runoff are used to determine water transit time in the small granitic Strengbach catchment (0·8 km²; 883–1146 m above sea level) located in the Vosges Mountains of northeastern France. Water transit times were calculated by applying the exponential, exponential piston and dispersion models of the FlowPC program to isotopic input (rainfall) and output (spring and stream water) data sets during the period 1989–95. The input function of the model was modified compared with the former version of the model and estimated by a deterministic approach based on a simplified hydrological balance. The fit between observed and calculated output data showed marked improvements compared with results obtained using the initial version of the model. An exponential piston version of the model applied to spring water indicates a 38·5 month mean transit time, which suggests that the volume in the aquifer, expressed in water depth, is 2·4 m. A considerable thickness (>45 m) of fractured bedrock may be involved for such a volume of water to be stored in the aquifer. Copyright © 2005 John Wiley & Sons, Ltd.     

广西老厂铅锌矿床地质特征及其成因探讨

本文是在专题研究基础上,探讨广西老厂铅锌矿床地质特征和矿床成因。该矿床位于海洋山复式背斜南段的老厂穹窿构造核部,铅锌矿化受构造、地层和岩性控制。矿体呈脉状,产于NE向组张扭性断裂构造中;围岩为寒武系白云岩,板岩、砂岩;矿石类型以锌为主,矿物成分简单,以闪锌矿、方铅矿与石英、白云石为主,具充填、交代矿石结构、构造,矿石平均品位含Zn9.94%、Pb1.50%;围岩蚀变主要为硅化、绢云母化和碳酸盐化。据地层岩石元素丰度值、稀土元素配分和稳定同位素数据、资料,探讨了成矿物质来源于围岩,热液主要来自大气降水。因此,认为本矿床实属产于寒武系轻变质岩中的层控矿型大气降水热液铅锌矿床。     

冰碛扇型砂金矿床特征及其成矿初探

成生于冰碛扇内，经一定的生物、物理、化学成矿作用，以砂金形态分布而富集成矿的金矿称为冰碛扇型砂金矿床。它具有沿造山带一定标高范围成群、成带分布，沉积物为近源且半胶结，砂金呈面状分布，并以粒金、块金为主，具胶状、浑圆状、包块状构造形态，分布极不均匀等特征。冰碛扇型砂金矿床的成矿作用与河成砂金矿床有本质的区别。它的成矿作用模式是：造山带下地壳富含的活化金通过新构造运动活化的断裂运移地表，经地下流体和地表流体的迁移聚集到冰碛扇这一封闭稳定的生物、物理、化学障环境中，主要经高效聚金微生物有机胶体成矿作用沉淀、再生加大形成砂金，逐步富集成矿床。而红色磨拉石建造中风化离解的Ｆｅ，Ｍｎ物质对沉淀环境ｐＨ值起着一定的调节作用。冰碛扇型砂金矿床在我国西部造山带内广泛分布，具有一定的工业价值，是砂金矿床中一重要类型，应予以重视     

Tephrochronology and tephrostratigraphy of two Pleistocene continental fossiliferous successions from central Italy

The stratigraphical context of two Middle Pleistocene fossiliferous palaeosols from Central Italy (Abruzzo and Tuscany) have been studied. Small mammals and molluscs occur in both palaeosols, which are covered by tephra layers that were analysed using an interdisciplinary approach. Application of fission‐track dating to apatites separated from the Case Picconetto tephra (Pescara, Abruzzo), yielded an age of 0.48 ± 0.04 Ma, indistinguishable from those previously determined for the Campani Quarry (Lower Valdarno, Tuscany) (0.46 ± 0.05 Ma and 0.48 ± 0.05 Ma). Geochemical and petrographic investigations indicate that these tephra originated from different volcanoes, the Alban Hills Volcanic Complex and the Vico Volcano (Latium) respectively. Small mammal and mollusc assemblages indicate different palaeoclimatic and palaeoenvironmental conditions for the Case Picconetto and Campani Quarry palaeosols. Warm and humid conditions can be inferred for the Campani Quarry site, whereas open and cold conditions can be inferred for Case Picconetto. On the basis of faunal data, fission‐track dates and attribution of tephra to specific volcanic eruptions, we suggest a correlation of these faunas with marine oxygen isotope stage 14 (Case Picconetto) and with marine oxygen isotope stage 11 (Campani Quarry), respectively. Copyright © 2003 John Wiley & Sons, Ltd.     

The response of an elastic half‐space under a momentary shear line impulse

The response of an ideal elastic half‐space to a line‐concentrated impulsive vector shear force applied momentarily is obtained by an analytical–numerical computational method based on the theory of characteristics in conjunction with kinematical relations derived across surfaces of strong discontinuities. The shear force is concentrated along an infinite line, drawn on the surface of the half‐space, while being normal to that line as well as to the axis of symmetry of the half‐space. An exact loading model is introduced and built into the computational method for this shear force. With this model, a compatibility exists among the prescribed applied force, the geometric decay of the shear stress component at the precursor shear wave, and the boundary conditions of the half‐space; in this sense, the source configuration is exact. For the transient boundary‐value problem described above, a wave characteristics formulation is presented, where its differential equations are extended to allow for strong discontinuities which occur in the material motion of the half‐space. A numerical integration of these extended differential equations is then carried out in a three‐dimensional spatiotemporal wavegrid formed by the Cartesian bicharacteristic curves of the wave characteristics formulation. This work is devoted to the construction of the computational method and to the concepts involved therein, whereas the interpretation of the resultant transient deformation of the half‐space is presented in a subsequent paper. Copyright © 2003 John Wiley & Sons, Ltd.     

Modelling of deformation response and chemo‐mechanical coupling in chalk

A constitutive relation is derived for describing the mechanical response of chalk. The approach is based on a phenomenological framework which employs chemo‐plasticity. The properties of the material are assumed to be affected by the physico‐chemical processes that occur through the interaction between the skeleton and the pore fluid. The underlying mechanism is discussed by invoking a micromechanical analysis. The performance of the framework is illustrated by examining the evolution of mechanical characteristics in the presence of different pore fluids. Copyright © 2006 John Wiley & Sons, Ltd.     

Use of multi‐platform,multi‐temporal remote‐sensing data for calibration of a distributed hydrological model: an application in the Arno basin,Italy

Images from satellite platforms are a valid aid in order to obtain distributed information about hydrological surface states and parameters needed in calibration and validation of the water balance and flood forecasting. Remotely sensed data are easily available on large areas and with a frequency compatible with land cover changes. In this paper, remotely sensed images from different types of sensor have been utilized as a support to the calibration of the distributed hydrological model MOBIDIC, currently used in the experimental system of flood forecasting of the Arno River Basin Authority. Six radar images from ERS‐2 synthetic aperture radar (SAR) sensors (three for summer 2002 and three for spring–summer 2003) have been utilized and a relationship between soil saturation indexes and backscatter coefficient from SAR images has been investigated. Analysis has been performed only on pixels with meagre or no vegetation cover, in order to legitimize the assumption that water content of the soil is the main variable that influences the backscatter coefficient. Such pixels have been obtained by considering vegetation indexes (NDVI) and land cover maps produced by optical sensors (Landsat‐ETM). In order to calibrate the soil moisture model based on information provided by SAR images, an optimization algorithm has been utilized to minimize the regression error between saturation indexes from model and SAR data and error between measured and modelled discharge flows. Utilizing this procedure, model parameters that rule soil moisture fluxes have been calibrated, obtaining not only a good match with remotely sensed data, but also an enhancement of model performance in flow prediction with respect to a previous calibration with river discharge data only. Copyright © 2006 John Wiley & Sons, Ltd.     

胶东“玲珑-焦家式”金矿资源潜力与找矿

“玲珑焦家式”金矿的成岩成矿,从新太古代至中生代存在着继承性的成因关系,成矿具有长期性、多源性特点,但中生代燕山早期岩浆活动是主导成矿因素。区域化探成果显示胶东金矿直接矿源岩———郭家岭花岗岩金的丰度比原始矿源岩———早前寒武纪变质沉积岩系及侵入岩系金的丰度低1/2。郭家岭花岗岩面积在2 km2以上岩体定量估计金亏损总量与胶东金矿总量预测的资源量相当,从一个侧面展现了胶东金矿资源的潜力。望儿山金矿的找矿实践展示,焦家式金矿在1 000 m以下尚有巨大潜力,整个胶东金矿产区在3 000 m以上的资源潜力应在5000 t以上。     

上扬子地台震旦系铅锌矿床类型及找矿方向

上扬子地台含铅锌矿地层主要为震旦系。震旦系在上扬子地台主要发育上统灯影组,次有下统陡山沱组。将上扬子地台震旦系铅锌矿床成因划分为一级为叠生矿床,二级为层控矿床,三级为热水沉积改造矿床和热水沉积再造矿床四级类型,该类型按铅锌比例进一步划分为富铅型、富锌型、铅锌共生型。研究区内最重要的铅锌矿床成因类型为南郑马元、会东大梁子铅锌矿为代表的富锌型。陕西南郑、紫阳—镇坪地区,湖北武当、神农架、黄陵地区,以及地台西缘的安宁河、甘洛—小江、峨边—寻甸等地区的震旦系是寻找铅锌矿的最有利地段。