A new scheme for numerical modelling of flow and transport processes in 3D fractured porous media

The objective of this work is to develop a new numerical approach for the three-dimensional modelling of flow and transient solute transport in fractured porous media which would provide an accurate and efficient treatment of 3D complex geometries and inhomogeneities. For this reason, and in order to eliminate as much as possible the number of degrees of freedom, the fracture network, fractures and their intersections, are solved with a coupled 2D–1D model while the porous matrix is solved independently with a 3D model. The interaction between both models is accounted for by a coupling iterative technique. In this way it is possible to improve efficiency and reduce CPU usage by avoiding 3D mesh refinements of the fractures. The approach is based on the discrete-fracture model in which the exact geometry and location of each fracture in the network must be provided as an input. The formulation is based on a multidimensional coupling of the boundary element method-multidomain (BEM-MD) scheme for the flow and boundary element dual reciprocity method-multidomain (BE-DRM-MD) scheme for the transport. Accurate results and high efficiency have been obtained and are reported in this paper.     

Tracer measurements in groyne fields for the quantification of mean hydraulic residence times and of the exchange with the stream

Tracer experiments were carried out in groyne fields (GF) of the river Elbe near Havelberg (Germany) in order to estimate the hydraulic connectivity with the river channel. The characteristic times of the five groyne fields, which were estimated from the exponentially declining tracer curves in 43 runs, ranged between 15 min and 69 min and did not correlate with the water level. Methodological investigations show that single point injection and two measurements (in the outflowing water and in the dominant region) are sufficient to provide robust in‐situ tracer curves. Using simplified mathematical simulations with connected stirred tanks, the conditions are investigated for the development of breaks in tracer curves and for the occurrence of significant errors in the estimation of intrinsic residence times. It was shown that an initial uniform dye distribution is not mandatory for the estimation during steady states. In special cases, point injections are more advisable. Moreover, the mean hydraulic residence time was found to be not equivalent to the estimated characteristic time. In fact, it is mostly overestimated by tracer experiments. The degree of overestimation depends on mixing and the volumetric proportions between the different parts of the GF and can be calculated from measured dye concentration differences. For example, an overestimation of 32% was calculated for a groyne field with a commonly found circulation flow pattern.     

A conceptual model for kimberlite emplacement by solitary interfacial mega-waves on the core mantle boundary

If convection in the Earth's liquid outer core is disrupted, degrades to turbulence and begins to behave in a chaotic manner, it will destabilize the Earth's magnetic field and provide the seeds for kimberlite melts via turbulent jets of silicate rich core material which invade the lower mantle. These (proto-) melts may then be captured by extreme amplitude solitary nonlinear waves generated through interaction of the outer core surface with the base of the mantle. A pressure differential behind the wave front then provides a mechanism for the captured melt to ascend to the upper mantle and crust so quickly that emplacement may indirectly promote a type of impact fracture cone within the relatively brittle crust. These waves are very rare but of finite probability. The assumption of turbulence transmission between layers is justified using a simple three-layer liquid model. The core derived melts eventually become frozen in place as localised topographic highs in the Mohorovicic discontinuity (Moho), or as deep rooted intrusive events. The intrusion's final composition is a function of melt contamination by two separate sources: the core contaminated mantle base and subducted Archean crust. The mega-wave hypothesis offers a plausible vehicle for early stage emplacement of kimberlite pipes and explains the age association of diamondiferous kimberlites with magnetic reversals and tectonic plate rearrangements.     

雅鲁藏布江蛇绿岩中超高压矿物硅尖晶石的研究

从西藏雅鲁藏布江蛇绿岩带的罗布莎蛇绿岩的铬铁矿中，发现一个由70-80种奇异矿物组成的地幔矿物群，其中包括一种成分特殊的尖晶石类矿物。该种尖晶石呈包裹体分布在毒砂中，28粒该矿物的平均化学成分：Na2O 1.58％，MgO7．52％，Al2O3 36．59％，SiO2 44．45％，FeO 8．72％，并含少量CaO和TiO2。经激光拉曼谱仪测试，一部分颗粒具有Franclinite（锌铁尖晶石ZnFe2O4）拉曼谱。根据尖晶石结构和化学成分，可以得出两种分子式：（Mg0．52Na0．14Fe0.32Al0.74）1.72（Si2.00Al1.20）3.20O8和（Mg0.52Na0.14Fe0.32Si0.50）1.48相（Si1.50Al1.94）3．44O8。两种分子式都表明阳离子Si呈六配位占据尖晶石八面体晶格位置。Si离子呈六配位的硅酸盐。实验证明具有超高压性质，来自相当于过渡带400-670km的深部。表明西藏雅鲁藏布江蛇绿岩（古大洋岩石圈）的岩浆活动达到过渡带。可能是地幔柱活动将硅尖晶石类等超高压矿物搬运到上地幔浅部的。     

公婆泉盆地下白垩统新民堡群层序地层和铀成矿特点

本文论述了公婆泉盆地下白垩统新民堡群的层序地层特点———层序分布和时代、界面特征和内部构成、沉积相,以及层间氧化带和铀矿化特点。认为该盆地在公婆泉地段具有一定的找铀矿远景,但北骆驼泉地段铀成矿条件不佳,那林疏怀地段不具备铀成矿远景。     

北京房山岩体形成过程中的岩浆混合作用证据

北京房山侵入岩体主要由石英二长闪长岩单元和三种不同结构的花岗闪长岩单元组成,岩体中广泛分布的镁铁质微粒包体、富集暗色矿物的线状条带以及各岩相单元的非渐变接触关系等宏观特征,指示岩体的形成与岩浆混合作用有密切联系。122个斜长石成分剖面电子探针分析表明,它们众数成分为 An=29的更长石。但是,其成分变化范围极宽(An=14～56),晶体内不同部位的成分有很大差别,出现正常环带的 An 值较低(14～42),内环带的 An 值较高(24～56)的特殊情况,指示斜长石结晶环境的改变,暗示具有富钙基性岩浆的加入。从岩体边部到中心,斜长石 An 牌号逐渐降低,但相邻岩相单元之间有较大的重叠,表明不同岩石单元中的斜长石成分主要受控于岩浆成分的变化,与分离结晶作用无关。斜长石的韵律性环带、特别是内环带的发育则是岩浆混合作用机制的有力证据。基性岩浆的反复注入不仅改变寄主岩浆的成分,而且也控制了它的温度、定位机制和岩石的结构特征。因此,房山岩体的形成过程还包含了镁铁质岩浆多次注入的历史,并且从外到里基性岩浆的作用逐渐减弱。     

同向圆射流混合特性实验研究

利用声学多普勒流速仪(Acoustic Doppler Velocimetry,ADV)自动测速系统和平面激光诱导荧光(Plane Laser Induced Florescence,PLIF)系统对同向圆射流的速度场和浓度场进行了系统的研究,并与前人的工作进行了比较和讨论。时均特性测量结果表明,射流主体段下游各断面的速度和浓度横向分布是自相似的,且近似可用高斯分布来表示;随着射流向下游传播,射流中心流速逐渐衰减,在强射流区和弱射流区遵循不同的规律;紊动特性测量结果表明,下游各断面的紊动强度横向分布是自相似的。     

多孔介质中水敏性应用的研究综述

水敏现象在石油、岩土和环境工程等领域广泛存在。有效地控制水敏现象的发生，不仅具有科学意义，还有一定的经济价值。本文在广泛查阅国内外研究成果的基础上，总结了水敏性在石油、岩土和环境工程方面应用的研究进展，并在此基础上提出了水敏性研究的发展方向。     

舟坝水电站水库蓄水对乐山市地震活动的影响及对策

在对铜街子水电站蓄水后库坝区及附近出现的地震活动特点进行分析的基础上,对另外一个舟坝水电站库坝区的地质结构和历史地震活动情况进行探讨,认为舟坝水电站蓄水产生诱发地震的可能性较大,五渡—利店断裂附近的马边、沐川、沙湾等地的诱发地震可能为5级以上,犍为县城附近可能出现的地震为4级左右。针对这些问题提出了减轻诱发地震造成损失的相应对策。     

Lithium isotope fractionation in the southern Cascadia subduction zone

We present lithium (Li) abundances and isotope compositions for a suite of anhydrous olivine tholeiites (HAOTs) and hydrous basalt-andesitic (BA) lavas from the Mt. Shasta and Medicine Lake regions, California. The values of δ⁷Li vary from + 0.9‰ to + 6.4‰ and correlate inversely with distance from the trench. These data are consistent with continuous isotope fractionation of Li during dehydration of the subducted oceanic lithosphere, an interpretation corroborated by uniformly high pre-eruptive H₂O contents in basaltic andesites accompanied by high Li, Rb, Sr, Ba and Pb abundances. The subduction-derived component that was added to these hydrous magmas is shown to be very similar beneath both Mt. Shasta and Medicine Lake volcanoes despite characteristically distinct Li isotope compositions in the magmas themselves. More evolved andesites and dacites from Mt. Shasta have δ⁷Li from + 2.8 to + 6.9‰ which is identical with the range obtained for HAOTs and BA lavas from Mt. Shasta. Therefore, Li isotopes do not provide evidence for any other crustal component admixed to Mt. Shasta andesites or dacites during magmatic differentiation and magma mixing in the crust.