对桐庐早白垩世火山—侵入岩成因的几点认识

本文报导了对桐庐火山－侵入岩进一步研究所获得的认识。单颗锆石Ｕ－Ｐｂ年龄测定表明桐庐火山－侵入岩是早白垩世火山活动的产物（１３４．４～１３４．９Ｍａ）。浅部岩浆房中从上往下的流纹质、流纹英安质和石英闪长质岩浆具有相近似的Ｒｂ－Ｓｒ和Ｓｍ－Ｎｄ同位素初始比,（８７Ｓｒ／８６Ｓｒ）ｉ＝０．７０８～０．７０９、εＮｄ（ｔ）＝－５．２７～－６．５６,高场强元素比（如Ｚｒ／Ｎｂ＝１０．１～１３．１）也相近似,表明盆地中各种火山－侵入杂岩是同源的。流纹英安岩中钾长石碎斑周围“组分过冷边”的存在说明岩穹壳内部酸性岩浆的冷凝、结晶要比酸性熔岩流来得慢。盆地中心浅成侵入岩中的中性岩包体的Ｒｂ－Ｓｒ和Ｓｍ－Ｎｄ同位素成分、Ｍｇ／（Ｍｇ＋Ｆｅ）以及Ｎａ／（Ｎａ＋Ｃａ）比表明它们是岩浆房下部石英闪长质岩浆的岩浆熔离作用产物     

平潭和漳州深成杂岩中斜长石捕虏晶与岩浆混合作用

在闽东南平潭和漳州火成杂岩中,广泛发育岩浆混合现象,表现为花岗质岩石中存在辉长－闪长质同深成岩墙和淬冷岩石包体。花岗质和闪长质岩石中普遍存在Ａｎ＞７０％的斜长石捕虏晶。这种成分的斜长石是辉长岩特有的,而不是任一种花岗岩的斑晶。在斜长石捕虏晶四周,都有与花岗岩中的斜长石成分相当的斜长石加生膜。斜长石捕虏晶可分为自形的、熔融状、筛孔状和碎片状等四类,它们各自形成于不同的动力环境     

初议中亚成矿域

给出了中亚成矿域的地质含义与地理范围。此巨型成矿域以发育中、晚元古代及晚古生代成矿作用为其特点。以超大型矿床为重点,并主要围绕金、铜讨论了中亚成矿域6种类型矿床的时空展布。近年取得的进展说明,中亚成矿域在矿产资源优势与潜力方面并不亚于环太平洋成矿域及特提斯成矿域。     

松树沟元古宙蛇绿岩Nd、Sr、Pb同位素地球化学特征

松树沟蛇绿岩中镁铁质岩的Sm-Nd全岩等时年龄为1030±46Ma,∈Nd(t)=+5.7,模式年龄在1422-1271Ma之间,矿物的内部等时年龄为983±140Ma。全岩等时年龄为其形成上限,矿物内部等时年龄为其变质年龄,说明该蛇绿岩形成于中元古代中、晚期。镁铁质岩的∈Nd(t)在+4.2-+6.9之间,是DMM与EMI两个地幔端元的混合产物;206Pb/204Pb在18.06-18.66之间,207Pb/204Pb和208Pb/204Pb值较高,分别在15.55-15.60和37.59-38.38之间,207Pb/204Pb和208Pb/204Pb对206Pb/204Pb的关系表明,岩浆源于具DUPAL异常的源区,(207Pb/204Pb)i相对(208Pb/204Pb)i更偏离NHRL可能是变质流体作用的反映。∈Sr(t)值较高且变化大,可与蛇绿岩类比,87Sr/86Sr比值变化大,可能与海水蚀变作用有关。据此推测,松树沟蛇绿岩代表洋壳残片。     

结晶岩的膨胀性能研究

采集了来自德国大陆深钻（ＫＴＢ）中的一组结晶岩,在实验室用位移记录仪或应变片接通电桥测量了岩石浸入不同流体中应变随时间的变化．结果表明：１．不同的岩石浸入水中,其膨胀速率和膨胀量是不同的,这与岩石的孔隙结构有关．２．同一岩样浸入极性和非极性的流体中,其膨胀性能是不同的,这与水分子的极化特性及其表面张力有关．３．很多情况下,膨胀是各向异性的,尤其是片麻岩．因此体积的增长可根据在３个互相垂直方向上的相对长度变化的总值来计算．可以认为,膨胀的主要因素是极性流体与岩石颗粒作用的结果,膨胀性能是由于颗粒表面的双扩散层造成的．实验证实了地壳岩石变形的机制强烈地受到流体的存在与否以及不同类型流体的影响．流体与岩石间相互作用的研究将对深入探讨地壳内流体的热过程、力学过程、油气运移以及岩石物理性质有重要意义．     

安徽沿江地区中生代岩浆岩的基本特点

安徽沿江地区位于长江中下游岩浆带中段,中生代岩浆活动强烈,类型复杂,特点突出。它具有夹心饼干式分带现象。中间为富碱的高钾钙碱性岩系和橄榄安粗岩系;两侧为钙碱性岩系。中间的岩浆岩具有高钾富碱、高Ｓｒ,低Ｍｇ、Ｃｒ、Ｎｉ,εＮｄ变化大但Ｉｓｒ变化小且较低以及岩浆演化的多样性为特点。     

哀牢山—金沙江富碱侵入岩及其与裂谷构造关系初步研究

沿哀牢山—金沙江分布的富碱侵入岩,年龄为50-30Ma,化学成分是高钾(K_2O=5.07-8.30％)富碱(K_2O+Na_2O=8.17-10.08％),LREE富集(∑Ce/∑Y=7.5-9.4),铕一般无异常(δEu=0.98-1.06),δS~(34)值小(1.11-3.27)和锶初始比低(i=0.706)等,以及成因演化受第三纪裂谷控制。     

中国东北地区新生代火山岩的年代学研究

用常规K-Ar、~(40)Ar-~(39)Ar、~(14)C等方法系统测定了东北地区新生代火山岩的年龄,获得近百个有效年龄数据。其精度和准确度都很高,与地质背景相当吻合。对于玄武岩中超镁铁质岩色体也进行了年龄测定尝试、发现包体中含继承氩,并可能存在~(36)Ar过剩。在阶段加热过程中,样品中氩的释放有两个高峰,一个在低温阶段(450—620℃),另一个在高温阶段(>1000℃),约50％的放射成因~(40)Ar在低温阶段释放。     

A semianalytical approach for solving first-order perturbation equations of dissolution-timescale reactive infiltration instability problems in fluid-saturated rocks

This paper presents a semianalytical approach for solving first-order perturbation (FOP) equations, which are used to describe dissolution-timescale reactive infiltration instability (RII) problems in fluid-saturated rocks. The proposed approach contains two parts because the chemical dissolution reaction divides the whole problem domain into two subdomains. In the first part, the interface-condition substitution strategy is used to derive the analytical expressions of purely mathematical solutions for the FOP equations in the upstream subdomain, where the dissolution chemical reaction is ceased and the FOP equations are weakly coupled. In the second part, the finite element method (FEM) is used to derive the analytical expressions of numerical solutions for the FOP equations in the downstream subdomain, where the dissolution chemical reaction needs to be considered and the FOP equations are strongly coupled so that it is impossible to derive purely mathematical solutions for them. Particular attention is paid to the development of the element-by-element forward marching strategy, which is associated with the use of the FEM for solving this new kind of scientific problem. The related analytical results demonstrated that (1) both the dynamic characteristic of a reactive infiltration system and the dimensionless wavenumber can have pronounced influences on the distribution of the FOP dimensionless acid concentration within the entire domain of the dissolution-timescale RII problems in fluid-saturated rocks and (2) the FOP dimensionless acid concentration distribution exhibits two significantly different patterns in the upstream and downstream subdomains of the dissolution-timescale RII system.     

Development of a GPGPU-parallelized hybrid finite-discrete element method for modeling rock fracture

The hybrid finite-discrete element method (FDEM) is widely used for engineering applications, which, however, is computationally expensive and needs further development, especially when rock fracture process is modeled. This study aims to further develop a sequential hybrid FDEM code formerly proposed by the authors and parallelize it using compute unified device architecture (CUDA) C/C++ on the basis of a general-purpose graphics processing unit (GPGPU) for rock engineering applications. Because the contact detection algorithm in the sequential code is not suitable for GPGPU parallelization, a different contact detection algorithm is implemented in the GPGPU-parallelized hybrid FDEM. Moreover, a number of new features are implemented in the hybrid FDEM code, including the local damping technique for efficient geostatic stress analysis, contact damping, contact friction, and the absorbing boundary. Then, a number of simulations with both quasi-static and dynamic loading conditions are conducted using the GPGPU-parallelized hybrid FDEM, and the obtained results are compared both quantitatively and qualitatively with those from either theoretical analysis or the literature to calibrate the implementations. Finally, the speed-up performance of the hybrid FDEM is discussed in terms of its performance on various GPGPU accelerators and a comparison with the sequential code, which reveals that the GPGPU-parallelized hybrid FDEM can run more than 128 times faster than the sequential code if it is run on appropriate GPGPU accelerators, such as the Quadro GP100. It is concluded that the GPGPU-parallelized hybrid FDEM developed in this study is a valuable and powerful numerical tool for rock engineering applications.