Smoothing and Change Point Detection for Gamma Ray Count Data

Gamma ray detectors are used to measure the natural radioactivity of rocks. For a number of boreholes drilled at a site the gamma ray detector is lowered into each borehole and counts of gamma ray emissions at different depths are recorded as the instrument is gradually raised to ground level. The profile of gamma counts can be informative about the geology at each location. The raw count data are highly variable, and in this paper we describe the use of adaptive smoothing techniques and change point models in order to identify changes in the geology based on the gamma logs. We formulate all our models for the data in the framework of the class of generalized linear models, and describe computational methods for Bayesian inference and model selection for generalized linear models that improve on existing techniques. Application is made to gamma ray data from the Castelreagh Waste Management Centre which served as a hazardous waste disposal facility for the Sydney region between March 1974 and August 1998. Understanding the geological structure of this site is important for further modelling the transport of pollutants beneath the waste disposal area.     

基于TOUGHREACT的成矿过程化学反应数值模拟——以虎头崖铅锌多金属矿床为例

虎头崖铅锌多金属矿床是青海省祁漫塔格成矿带中典型的矽卡岩型矿床。本文基于前人对该矿区地质特征、矿物学和地球化学特征等方面的研究成果，利用TOUGHREACT数值模拟软件，针对石英-晚期硫化物成矿阶段进行成矿过程化学反应数值模拟研究。通过建立地质概念模型和网格化物化参数的设置，在设置的时间步长内对多相流体流动和热对流传导、溶质运移以及地球化学反应过程进行模拟计算，研究了矿区内不同温压条件下成矿元素化学平衡浓度变化的总体趋势以及对方铅矿和闪锌矿矿物溶解沉淀情况的影响。结果表明：当温度从250℃降至155℃时，Pb²⁺的化学平衡浓度出现骤减；当温度从300℃降至100℃时，Zn²⁺的化学平衡浓度出现骤减；由于矿物的沉淀需要消耗大量相关的成矿元素，方铅矿和闪锌矿最佳成矿温度范围分别为155℃~250℃和100℃~300℃；不同温度和压力对成矿元素化学平衡浓度影响不同，温度是控制方铅矿和闪锌矿沉淀化学反应的核心因素，压力的影响较小。模拟结果与传统地质解译结果相对符合，证明了模拟的有效性。     

高温高压下水岩界面反应与金的迁移富集

高温高压下水岩界面反应与金的迁移富集王玉荣阚小凤王志强（中国科学院广州地球化学研究所，广州５１０６４０）关键词水岩界面反应自发电池实验地球化学本实验以哈达门沟［１］及东坪金矿［２］为地质背景，研究钾长石岩中金矿化富集的物理化学条件。该区金矿既与早期变...     

表生环境下闪锌矿氧化反应模拟实验

采用混合流反应器装置,进行了表生环境条件下闪锌矿氧化动力学模拟实验研究,选择Fe³⁺和O₂为氧化剂,实验条件为:温度20～55℃,pH=1.0～7.8、氧化剂Fe³⁺浓度1.0×10^-5～1.0×10^-2 mol/L,氧气流量0.5 L/min。实验结果表明,在Fe³⁺为氧化剂时,闪锌矿氧化速率随着Fe³⁺浓度的增加、温度升高、pH值降低而增加,且氧化过程中Zn和Cd的释放速率大致相同,其反应的活化能分别为Ea(Zn) 41.75 kJ/mol或Ea(Cd) 42.51 kJ/mol,表明闪锌矿氧化速率主要受矿物表面反应控制;而以氧气作为氧化剂时,随着pH值的变化闪锌矿的氧化机理发生变化,在pH值小于5.95范围内,闪锌矿氧化速率随着pH值的增加而降低,进入碱性范围后,随着pH值增加,闪锌矿氧化速率反而增加。对实验数据进行双对数处理并建立了闪锌矿的氧化速率公式。     

多孔介质中含均相/非均相反应传质过程数值模拟

结合化学反应方程式,并应用多孔多相介质溶混污染物输运过程的数值模型,对多孔多相介质中含均相/非均相化学反应传质过程进行了数值模拟。化学反应主要包含均相快速/慢速和非均相快速/慢速等5种化学反应过程,溶质输运行为的控制机制主要考虑对流、扩散及降解、吸附等。基于原有的隐式特征线Galerkin离散化的有限元方法,求解模型控制方程的边值初值问题,求解过程中把均相化学反应物质中按照反应物和生成物分开,非均相反应物质按照固相和液相分开,对均相反应物及非均相液相物质浓度耦合求解,而均相生成物和非均相固相物质独立求解。使方程组按照其不同类型进行分类,同时可减少未知数的个数。对于含有非线性内状态变量的右端项进行迭代求解。数值例题结果验证了所提出的数值方法的有效性、计算精度和稳定性。     

海底硬岩钻探的现状与前景分析

海底硬岩是世界大洋钻探实现“莫霍钻”目标必然钻遇的地层。世界大洋钻探实施的航次中已多次钻遇海底硬岩,不仅采获了岩心样品,也在钻进施工中发现了诸多问题。本文搜集了近年来世界大洋钻探实施的309、312、335、360和384等有关海底硬岩钻探航次的资料,在对航次情况简要介绍的基础上,重点阐述了目前海底硬岩钻探在钻头、井壁稳定以及取心工具等方面遇到的挑战以及采取的应对措施,并对今后海底硬岩钻探的发展进行了展望。     

碳酸锰氧化实验和热力学分析

通过对碳酸锰氧化实验的结果分析和讨论，建立了一种适合于锰氧化过程的非平衡热力学模型，并利用此模型对硫酸锰氧化过程作了热力学可行性分析，研究结果表明，碳酸锰氧化过程是矿物组合的变化过程，也是锰自身价态的转为过程，整个过程既存在平衡状态，又存在非平衡状态，是一个不可逆的，多相的复杂反应。     

Petrochemical constraints on the models of Cretaceous-Eocene tectonic evolution of the Eastern Pontic chain (Turkey)

Sixteen selected samples from the Upper Cretaceous volcanic belt of the Eastern Pontids have been analysed for major elements, Rb, Sr and Zr. On the basis of the K₂O versus SiO₂ distribution, two groups of rocks have been distinguished, one with calc-alkaline affinity and a second group with shoshonitic character. The calc-alkaline rocks have porphyritic texture with clinopyroxene, plagioclase and orthopyroxene as phenocryst and in the groundmass. The orthopyroxene is lacking in the shoshonites where plagioclase, clinopyroxene and, in the more evolved terms, amphibole and biotite are the main phenocryst minerals. The shoshonitic rocks have higher $\text{K}{}_2\text{O}\text{Na}{}_2\text{O}$ ratio, K₂O, P₂O₅ and Rb, contents with respect to the calc-alkaline samples. The TiO₂ content is invariably low, never exceeding approximately 1%. The occurrence of volcanic rocks ranging in composition from calc-alkaline to shoshonitic in the Upper Cretaceous volcanic belt of the Eastern Pontids suggests that the Upper Cretaceous volcanic cycle reached its mature stage before the onset of the Eocene calc-alkaline volcanism which is believed to be neither genetically nor tectonically related with the Upper Cretaceous volcanism.     

A garnet population in Yellowknife schist, Canada

Abstract Data are presented on a garnet population in a specimen of garnet-biotite-plagioclase-quartz schist from the cordierite zone of an Archaean thermal dome in the Southern Slave Province of the Canadian Shield. Garnet crystals are bounded by planar dodecahedral faces and by trapezohedral faces which on the 10-μm scale are corrugated. Crystal distribution, as revealed by dissection of a small cubic volume of rock, is random. The size distribution is normal, with a mean diameter of 0.81 mm and a standard deviation of 0.32 mm. In the largest crystal of the population (mean radius 0.83 mm), [Mn] = 100 Mn/(Fe + Mg + Mn + Ca) decreases from 14.5 at the centre to 7.5 and then increases in the outer margin to 8.5; [Fe] increases continuously from 67 at the centre to 77 at the surface; [Mg] increases from 12.5 to 13.5 and then falls sharply to 11; [Ca] remains unchanged at 4.0 and then drops to 3.3. Progressively smaller crystals have progressively lower [Mn] and higher [Fe] concentrations at their centres, while all crystals have the same margin composition. Growth vectors extending from given concentration contours to crystal surfaces are of equal length regardless of the size of the crystal in which the vector is located. A garnet-forming model is presented in which reaction was initiated by a rise in temperature. Nucleation sites were randomly selected. The nucleation rate increased with time and then declined. Crystal faces advanced at a constant linear rate, which implies an increase in volume proportional to surface area. Initially, the composition of garnet deposited on crystal surfaces was determined by van Laar equations of equilibrium, which demanded the withdrawal of Mn and Fe from within chlorite crystals. This transfer reaction was then accompanied by an ion exchange reaction which moved Mn and Fe to garnet surfaces from biotite, in exchange for Mg. The exchange reaction provides an explanation for the high overall concentration of Mn and Fe in garnet and for the observed Mn and Mg reversals in the margins of crystals. The increase of garnet volume in the garnet population is found to be parabolic, i.e. Vαα⁵.