The reflection of plane waves in a poroelastic half-space saturated with inviscid fluid

This paper discusses surface displacements, surface strain, rocking, and energy partitioning during reflection-of-plane waves in a fluid-saturated poroelastic half-space. The medium is modeled by Biot's theory, and is assumed to be saturated with inviscid fluid. A linear porosity-modulus relation based on experimental data on sandstones is used to determine the material parameters for Biot's model. Numerical results in terms of angle of incident waves and Poisson's ratio are illustrated for various porosities and degrees of solid frame stiffness. The results show that the amount of solid frame stiffness controls the response of a fluid-saturated porous system. A poroelastic medium with essentially dry-frame stiffness behaves like an elastic medium, and the influence of pore fluid increases as dry-frame stiffness is reduced. The effects of a second P-wave become noticeable in poroelastic media with low dry-frame stiffness.     

The role of moisture cycling in the weathering of a quartz chlorite schist in a tropical environment: findings of a laboratory simulation

Long‐term weathering of a quartz chlorite schist via wetting and drying was studied under a simulated tropical climate. Cubic rock samples (15 mm × 15 mm × 15 mm) were cut from larger rocks and subjected to time‐compressed climatic conditions simulating the tropical wet season climate at the Ranger Uranium Mine in the Northern Territory, Australia. Fragmentation, moisture content and moisture uptake rate were monitored over 5000 cycles of wetting and drying. To determine the impact of climatic variables, five climatic regimes were simulated, varying water application, temperature and drying. One of the climatic regimes reproduced observed temperature and moisture variability at the Ranger Uranium Mine, but over a compressed time scale. It is shown that wetting and drying is capable of weathering quartz chlorite schist with changes expected over a real time period of decades. While wetting and drying alone does produce changes to rock morphology, the incorporation of temperature variation further enhances weathering rates. Although little fragmentation occurred in experiments, significant changes to internal pore structure were observed, which could potentially enhance other weathering mechanisms. Moisture variability is shown to lead to higher weathering rates than are observed when samples are subjected only to leaching. Finally, experiments were conducted on two rock samples from the same source having only subtle differences in mineralogy. The samples exhibited quite different weathering rates leading to the conclusion that our knowledge of the role of rock type and composition in weathering is insufficient for the accurate determination of weathering rates. Copyright © 2005 John Wiley & Sons, Ltd.     

储层岩石电性特征的逾渗网络模型研究(英文)

建立能够反映储层孔隙结构、流体特征的逾 渗网络模型,通过数值模拟研究了储层孔隙尺寸、 孔隙形状、连通性、微孔隙的发育状况等对I-Sw曲 线影响的定性规律,研究了地层水矿化度对岩石电 阻率的影响。分析讨论了不同因素对电阻率影响的 相对强度。最后,通过曲线拟合定量地研究了不同 因素对I-Sw曲线的影响规律。研究表明,在上述影 响因素中连通性和微孔隙对I-Sw曲线的影响很大, 其它因素的影响强度较小。地层水矿化度对岩石电 阻率绝对值的影响很大。“非阿尔奇”现象普遍存在, 在低渗透(低连通性)储层岩石中更为明显。     

Cause of winter gravity wave spectrum saturation

Gravity waves play a significant role in establishing the large-scale circulation and structure of the middle atmosphere. Through gravity wave saturation proc-esses, such motions are believed to cause turbulence, resulting in divergence of momentum flux and the diffusion of heat and constituents in the meso-sphere[1,2]. The mechanisms that contribute signifi-cantly to the gravity wave saturation are thought to be the dynamical and convective instabilities[3]. However, it is difficult to distin…     

致密砂岩的形成机制及其地质意义——以塔里木盆地英南2井为例

砂质岩的孔隙和喉道被网格状粘土矿物和次生加大矿物充填成微细孔喉状结构时,形成具有较高毛细管压力的致密砂岩。它与地层水发生水锁效应,可大大降低渗透率,成为致密砂岩盖层。当含水饱和度在50%以下时,束缚水饱和度比较低,致密砂岩储层可以产气;当含水饱和度在50%～90%区间时,具有较高的束缚水饱和度,相对渗透率非常低,它既不能产气也不能产水,反映为渗透率瓶颈区(具有盖层性质);当含水饱和度大于90%以上时,致密砂岩储层仅微量产水。塔里木盆地英南2井侏罗系气藏盖层由致密砂岩构成,不含水时的气体渗透率在(0.027～0.081)×10-3μm2,不能构成封堵;当含水饱和度达到60%以上时,相对渗透率几乎为零,构成有效盖层。     

Effect of suction on the mechanical behaviour of iron ore rock

The effect of suction on the behaviour of iron ore has been studied from both physical and mechanical points of view. The porosity and the suction phenomena have been analysed using different experimental techniques. Uniaxial compressive tests on partially saturated samples have shown that the suction is responsible for strength and cohesion improvement. Considering the theory of partially saturated porous soils of Coussy and Dangla (Mécanique des sols non saturés (2002 edn). Hermès Science: 2002; 390), we have proposed a constitutive law for partially saturated iron ore. The real increase in the apparent cohesion due to the capillary attraction forces is overestimated if the yield function is written in terms of effective stresses. The effect of the capillary cohesion has been modelled with a function in the expression of the apparent cohesion of the yield function. The effect of suction on the mechanical behaviour has been represented in the effective stresses space and in the total stresses space like the Alonso model (Géotechnique 1990; 40 :405–430). Copyright © 2005 John Wiley & Sons, Ltd.     

Experimental determination of argon solubility in silicate melts: An assessment of the effects of liquid composition and temperature

The argon solubility of 38 liquids in the system Na₂O-CaO-MgO-Al₂O₃-SiO₂ (NCMAS) has been determined at 1873 K and 1 bar, the argon concentration of presaturated glasses being measured using a static mass spectrometer. For compositions in the subsystem diopside (CaMgSi₂O₆), nepheline (NaAlSiO₄), albite (NaAlSi₃O₈), anorthite (CaAl₂Si₂O₈), argon solubility is generally a linear function of the relative proportion of each end member, solubility being lowest in diopside melt (1.53 10⁻⁵ cm³ STP · g⁻¹ · bar⁻¹) and highest in albite melt (2.88 10⁻⁴ cm³ STP · g⁻¹ · bar⁻¹). For the tectosilicate joins studied (SiO₂-Na₂Al₂O₄, SiO₂-CaAl₂O₄, SiO₂-MgAl₂O₄) solubility decreases with decreasing silica content in all cases, being highest for Na-bearing liquids and lowest for Mg-bearing liquids at constant molar silica content. Where comparison is possible our results are in good agreement with data from the literature. When our data are considered in isolation we find that argon solubility shows an excellent correlation with calculated ionic porosity. The covariation of argon solubility and liquid density is also reasonable, that with molar volume less convincing and that with polymerization state (as defined by the ratio of the number of nonbridging oxygens and tetrahedral network forming cations; NBO/T) nonexistent. However, when our data are combined with those from the literature no well constrained correlation between argon solubility and ionic porosity is apparent. Based upon this observation and consideration of the temperature dependence of noble gas solubility it is concluded that ionic porosity is not a universally applicable parameter which may be used to predict noble gas solubility as a function of composition, temperature and pressure. Two new models for calculating argon solubility are proposed, both employing the notion of partial molar argon solubilities. The first uses oxide components, for which partial molar argon solubility is directly proportional to partial molar ionic porosity calculated at 1873 K, irrespective of the temperature of experimental equilibration. The second model, which offers the best fit to the available data, employs tetrahedral units rather than oxides as the proposed melt components. This latter model successfully accounts for reported argon solubilities in simple Al-free systems, in simple Al-bearing systems and in natural liquids. This is interpreted to infer that argon is incorporated in large sites in the liquid structure (such as the space within rings of n-tetrahedra) although further work is required to understand the quantitative links between melt structure and noble gas solubility.     

Development of a soil moisture‐based distributed hydrologic model for determining hydrologically based critical source areas

A simple grid cell‐based distributed hydrologic model was developed to provide spatial information on hydrologic components for determining hydrologically based critical source areas. The model represents the critical process (soil moisture variation) to run‐off generation accounting for both local and global water balance. In this way, it simulates both infiltration excess run‐off and saturation excess run‐off. The model was tested by multisite and multivariable evaluation on the 50‐km² Little River Experimental Watershed I in Georgia, U.S. and 2 smaller nested subwatersheds. Water balance, hydrograph, and soil moisture were simulated and compared to observed data. For streamflow calibration, the daily Nash‐Sutcliffe coefficient was 0.78 at the watershed outlet and 0.56 and 0.75 at the 2 nested subwatersheds. For the validation period, the Nash‐Sutcliffe coefficients were 0.79 at the watershed outlet and 0.85 and 0.83 at the 2 subwatersheds. The per cent bias was less than 15% for all sites. For soil moisture, the model also predicted the rising and declining trends at 4 of the 5 measurement sites. The spatial distribution of surface run‐off simulated by the model was mainly controlled by local characteristics (precipitation, soil properties, and land cover) on dry days and by global watershed characteristics (relative position within the watershed and hydrologic connectivity) on wet days when saturation excess run‐off was simulated. The spatial details of run‐off generation and travel time along flow paths provided by the model are helpful for watershed managers to further identify critical source areas of non‐point source pollution and develop best management practices.     

A methodology for regionalizing 3‐D effective porosity at watershed scale in crystalline aquifers

An innovative approach for regionalizing the 3‐D effective porosity field is presented and applied to two large, overexploited, and deeply weathered crystalline aquifers located in southern India. The method derives from earlier work on regionalizing a 2‐D effective porosity field in that part of an aquifer where the water table fluctuates, which is now extended over the entire aquifer using a 3‐D approach. A method based on geological and geophysical surveys has also been developed for mapping the weathering profile layers (saprolite and fractured layers). The method for regionalizing 3‐D effective porosity combines water table fluctuation and groundwater budget techniques at various cell sizes with the use of satellite‐based data (for groundwater abstraction), the structure of the weathering profile, and geostatistical techniques. The approach is presented in detail for the Kudaliar watershed (983 km²) and tested on the 730 km² Anantapur watershed. At watershed scale, the effective porosity of the aquifer ranges from 0.5% to 2% in Kudaliar and between 0.3% and 1% in Anantapur, which agrees with earlier works. Results show that (a) depending on the geology and on the structure of the weathering profile, the vertical distribution of effective porosity can be very different and that the fractured layers in crystalline aquifers are not necessarily characterized by a rapid decrease in effective porosity and (b) that the lateral variations in effective porosity can be larger than the vertical ones. These variations suggest that within a same weathering profile, the density of open fractures and/or degree of weathering in the fractured zone may significantly vary from a place to another. The proposed method provides information on the spatial distribution of effective porosity that is of prime interest in terms of flux and contaminant transport in crystalline aquifers. Implications for mapping groundwater storage and scarcity are also discussed, which should help in improving groundwater resource management strategies.     

Zircon in amphibolites from Naxos,Aegean Sea,Greece: origin,significance and tectonic setting

We report U–Pb zircon ages of c. 700–550 Ma, 262–220 Ma, 47–38 Ma and 15–14 Ma from amphibolites on Naxos Island in the Aegean extensional province of Greece. The zircon has complex internal structures. Based on cathodoluminescence response, zoning and crosscutting relationships a minimum of four zircon growth stages are identified: inherited core, magmatic core, inner metamorphic (?) rim and an outer metamorphic rim. Trace element compositions of the amphibolites suggest igneous differentiation and crustal assimilation. Zircon solubility as a function of saturation temperatures, Zr content and melt composition indicates that the zircon did not originally crystallize in the mafic bodies but was inherited from felsic precursor rocks, and subsequently assimilated into the mafic intrusives during emplacement. Zircon inheritance is corroborated by the complex, xenocrystic nature of the zircon in one sample. Ages of c. 700–550 Ma and 262–220 Ma are assigned to inherited zircon. Available geochemical data suggest that the 15–14 Ma metamorphic rims grew in situ in the amphibolites, corresponding to a high‐grade metamorphic event at this time. However, the geochemical data cannot conclusively establish if the c. 40 Ma zircon rims also grew in situ, or whether they were inherited along with the xenocrystic cores. Two scenarios for emplacement of the mafic intrusives are discussed: (i) Intrusion during late‐Triassic to Jurassic ocean basin development of the Aegean realm, in which case the 40 Ma zircon rims would have grown in situ, and (ii) emplacement in the Miocene as a result mafic underplating during large‐scale extension. In this case, only the 15–14 Ma metamorphic outer rims would have formed in situ in the amphibolitic host rocks.