地球化学模拟及其在水文地质中的应用

从正向地球化学模拟、反向地球化学模拟和水文地球化学模拟三方面系统概述了地球化学模拟的研究进展；结合笔者近年来的科研实践，对其在水文地质中的应用作了介绍，并着重讨论了地下水化学组分存在形式、矿物饱和指数的计算、水化学组分测量及地下水地球化学演化等方面的应用；最后，指出了这一领域中当前存在的问题和今后的研究方向。     

The use of electron backscatter diffraction and orientation contrast imaging as tools for sulphide textural studies: example from the Greens Creek deposit (Alaska)

The Greens Creek polymetallic massive sulphide deposit is hosted in a typical polyphase deformed lower greenschist facies orogenic setting. The structure of the host rocks is well constrained, exhibiting a series of three superimposed ductile deformations followed by two brittle episodes. The ore is found both in fold hinges where early-formed depositional features are preserved and in fold limbs where primary features are typically strongly modified or obliterated. Samples from both settings have been investigated using electron backscatter diffraction (EBSD) coupled with forescatter orientation contrast (OC) imaging in order to observe the effects of deformation and lower greenschist facies metamorphism on pyrite. Results suggest that colloform pyrite may preserve information relevant to palaeoenvironment, that apparently simple textures are generally more complex, and that pyrite can deform plastically by dislocation glide and creep processes at lower temperatures and/or strain rates than generally accepted. This analysis indicates that EBSD and OC imaging provide powerful tools for observing textural relationships in pyrite that are not shown by more traditional methods. They should become routine tools for pyrite texture analysis.Editorial handling: D. Lentz     

Komatiites and nickel sulfide ores of the Black Swan area,Yilgarn Craton,Western Australia. 4. Platinum group element distribution in the ores,and genetic implications

The Black Swan komatiite sequence, in the Eastern Goldfields province of the Archaean Yilgarn Craton in Western Australia, is a body of dominantly olivine-rich cumulates with lesser volumes of spinifex textured rocks, interpreted as a section through an extensive komatiite lava flow field. The sequence hosts a number of nickel sulfide orebodies, including the Silver Swan massive shoot and the Cygnet and Black Swan disseminated orebodies. The massive sulfide orebodies of the Black Swan Succession are pervasively depleted in all platinum group elements (PGEs), particularly Pt and Pd, despite very high Ni contents. This depletion cannot be explained by R-factor variations, which would also require relatively low Ni tenors. The PGE depletion could be explained in part if the ores are enriched in a monosulfide solid solution (MSS) cumulate component, but requires some additional fractional segregation of sulfide melt upstream from the site of deposition. The Silver Swan orebody shows a remarkably consistent vertical zonation in PGE contents, particularly in Ir, Ru, Rh, Os, which increase systematically from very low levels at the stratigraphic base of the sulfide body to maxima corresponding roughly with the top of a lower layer of the orebody rich in silicate inclusions. Platinum shows the opposite trend, but is somewhat modified by remobilisation during talc carbonate alteration. A similar pattern is also observed in the adjacent White Swan orebody. This zonation is interpreted and modelled as the result of fractional crystallisation of MSS from the molten sulfide pool. The strong IPGE depletion towards the base of the orebody may be a consequence of sulfide liquid crystallisation in an inverted thermal gradient, between a thin rapidly cooling upper rind of komatiite lava and a hot substrate.Electronic Supplementary Material Supplementary material is available in the online version of this article at Editorial handling: Peter Lightfoot     

硫化物矿物氧化反应动力学实验研究

硫化物矿床和含硫化物矿床（如中高硫煤）在开发过程中由于氧化作用向环境释放大量的酸、硫酸盐、重金属及其它有害物质,造成一定范围内的次生地球化学异常和环境污染;而中国大陆硫化物的排放量2010年将达到39．1×10 6 t。因此,人们对硫化物矿物在表生环境下的氧化行为日益重视。硫化物矿物的氧化速率受矿物本身物理化学特性（组成、结构、表面性质）和介质条件（溶解氧DO、Fe 3＋浓度、pH、Eh、温度及微生物）等多种因素的控制。在全面分析近年来国内外文献的基础上,系统介绍了硫化物矿物氧化动力学实验研究方法、实验参数的选择、反应机理、反应产物及反应速率等的研究现状。指出今后应加强对初始阶段的氧化反应速率、反应中间产物、多因素的综合作用、分子水平上的反应机理解释等方面的研究。     

超稠油油藏转小井距蒸汽驱条件及时机

蒸汽驱是稠油油藏经过蒸汽吞吐开采后进一步提高原油采收率的主要热采阶段。但不是说任何时候转驱都可取得较好的效果。笔者系统研究了影响蒸汽驱开采效果的因素：原油粘度、油层有效厚度、纯总厚度比、起始含油饱和度和油层非均质性等,确定出适宜蒸汽驱开采的油藏条件。在分析稠油油藏加密吞吐后油藏条件变化的基础上．利用油藏工程数值模拟技术研究确定出超稠油油藏转小井距蒸汽驱的合理时机。     

Geochemistry and geothermometry of volcanic rocks from Serra Branca, Iberian Pyrite Belt, Portugal

Volcanic rocks from Serra Branca, Iberian Pyrite Belt, Portugal, consist of calc-alkaline felsic and intermediate rocks. The latter are massive andesites, whereas the former include four dacitic to rhyolitic lithologies, distinguishable on spiderdiagrams and binary plots of immobile elements. Zircon thermometry indicates that two felsic suites may have formed from different magmas produced at distinct temperatures, with only limited fractionation within each suite. Alternatively, all the felsic rocks can be related through fractionation of a single magma if the lower zircon saturation temperature obtained for one suite merely results from Zr dilution, mostly reflecting silicification.

The relatively high magma temperatures at Serra Branca ease the classification of felsic rocks based on their HFSE contents and also indicate volcanogenic massive sulfide deposit favorability. This contrasts with other areas of the Belt that register lower magma temperatures and are subsequently barren. However, magma temperatures may have not been high enough to cause complete melting of refractory phases in which HFSE reside during crustal fusion of an amphibolite protolith, implying difficult discrimination of tectonic environments for the felsic rocks. The intermediate rocks were possibly formed by mixing between basaltic magmas and crustal material, compatible with volcanism in an attenuated continental lithosphere setting.     

Effect of loading rate on the volume change behavior of unsaturated residual soil

Residual soils occur in most countries of the world but those that occur in greater area and depth are usually found in the tropical humid areas. This research examines the effects of loading rate of applied or net mean stress on the volume change behavior of an unsaturated granitic residual soil. Three loading rates were adopted: fast, slow and extra slow loading rate. It was observed that the loading rate of applied or net mean stress has a pronounced effect on the void ratio and degree of saturation but has an insignificant effect on the water content of residual soils subjected to constant matric suction.     

The Aguablanca Ni–(Cu) sulfide deposit, SW Spain: geologic and geochemical controls and the relationship with a midcrustal layered mafic complex

The Aguablanca Ni–(Cu) sulfide deposit is hosted by a breccia pipe within a gabbro–diorite pluton. The deposit probably formed due to the disruption of a partially crystallized layered mafic complex at about 12–19 km depth and the subsequent emplacement of melts and breccias at shallow levels (<2 km). The ore-hosting breccias are interpreted as fragments of an ultramafic cumulate, which were transported to the near surface along with a molten sulfide melt. Phlogopite Ar–Ar ages are 341–332 Ma in the breccia pipe, and 338–334 Ma in the layered mafic complex, and are similar to recently reported U–Pb ages of the host Aguablanca Stock and other nearby calc-alkaline metaluminous intrusions (ca. 350–330 Ma). Ore deposition resulted from the combination of two critical factors, the emplacement of a layered mafic complex deep in the continental crust and the development of small dilational structures along transcrustal strike-slip faults that triggered the forceful intrusion of magmas to shallow levels. The emplacement of basaltic magmas in the lower middle crust was accompanied by major interaction with the host rocks, immiscibility of a sulfide melt, and the formation of a magma chamber with ultramafic cumulates and sulfide melt at the bottom and a vertically zoned mafic to intermediate magmas above. Dismembered bodies of mafic/ultramafic rocks thought to be parts of the complex crop out about 50 km southwest of the deposit in a tectonically uplifted block (Cortegana Igneous Complex, Aracena Massif). Reactivation of Variscan structures that merged at the depth of the mafic complex led to sequential extraction of melts, cumulates, and sulfide magma. Lithogeochemistry and Sr and Nd isotope data of the Aguablanca Stock reflect the mixing from two distinct reservoirs, i.e., an evolved siliciclastic middle-upper continental crust and a primitive tholeiitic melt. Crustal contamination in the deep magma chamber was so intense that orthopyroxene replaced olivine as the main mineral phase controlling the early fractional crystallization of the melt. Geochemical evidence includes enrichment in SiO₂ and incompatible elements, and Sr and Nd isotope compositions (⁸⁷Sr/⁸⁶Sr_i 0.708–0.710; ¹⁴³Nd/¹⁴⁴Nd_i 0.512–0.513). However, rocks of the Cortegana Igneous Complex have low initial ⁸⁷Sr/⁸⁶Sr and high initial ¹⁴³Nd/¹⁴⁴Nd values suggesting contamination by lower crustal rocks. Comparison of the geochemical and geological features of igneous rocks in the Aguablanca deposit and the Cortegana Igneous Complex indicates that, although probably part of the same magmatic system, they are rather different and the rocks of the Cortegana Igneous Complex were not the direct source of the Aguablanca deposit. Crust–magma interaction was a complex process, and the generation of orebodies was controlled by local but highly variable factors. The model for the formation of the Aguablanca deposit presented in this study implies that dense sulfide melts can effectively travel long distances through the continental crust and that dilational zones within compressional belts can effectively focus such melt transport into shallow environments.Electronic supplementary material Supplementary material is available in the online version of this article at and is accessible for authorized users.     

Quantifying Damage, Saturation and Anisotropy in Cracked Rocks by Inverting Elastic Wave Velocities

Crack damage results in a decrease of elastic wave velocities and in the development of anisotropy. Using non-interactive crack effective medium theory as a fundamental tool, we calculate dry and wet elastic properties of cracked rocks in terms of a crack density tensor, average crack aspect ratio and mean crack fabric orientation from the solid grains and fluid elastic properties. Using this same tool, we show that both the anisotropy and shear-wave splitting of elastic waves can be derived. Two simple crack distributions are considered for which the predicted anisotropy depends strongly on the saturation, reaching up to 60% in the dry case. Comparison with experimental data on two granites, a basalt and a marble, shows that the range of validity of the non-interactive effective medium theory model extends to a total crack density of approximately 0.5, considering symmetries up to orthorhombic. In the isotropic case, Kachanov's (1994) non-interactive effective medium model was used in order to invert elastic wave velocities and infer both crack density and aspect ratio evolutions. Inversions are stable and give coherent results in terms of crack density and aperture evolution. Crack density variations can be interpreted in terms of crack growth and/or changes of the crack surface contact areas as cracks are being closed or opened respectively. More importantly, the recovered evolution of aspect ratio shows an exponentially decreasing aspect ratio (and therefore aperture) with pressure, which has broader geophysical implications, in particular on fluid flow. The recovered evolution of aspect ratio is also consistent with current mechanical theories of crack closure. In the anisotropic cases—both transverse isotropic and orthorhombic symmetries were considered—anisotropy and saturation patterns were well reproduced by the modelling, and mean crack fabric orientations we recovered are consistent with in situ geophysical imaging. Our results point out that: (1) It is possible to predict damage, anisotropy and saturation in terms of a crack density tensor and mean crack aspect ratio and orientation; (2) using well constrained wave velocity data, it is possible to extrapolate the contemporaneous evolution of crack density, anisotropy and saturation using wave velocity inversion as a tool; 3) using such an inversion tool opens the door in linking elastic properties, variations to permeability.     

近源地震动峰值加速度衰减关系影响因素分析

本文收集了丰富的强震资料, 以峰值加速度为例, 采用简单且体现近场峰值加速度PGA震级饱和和距离饱和特性的衰减模型, 研究了表征震级与距离饱和效应的R0(M)的性质。R0与震级相关, 同时与震源性质、地震波频谱有关。在单个地震的R0(M)的求取中, 由于R0与系数d几乎呈线性关系, 所以要求单个地震R0的值, 必须先根据理论约束确定d的大小。在检验衰减方程的预测效果时, 不仅要判断衰减曲线是否反映了实测资料的平均变化趋势, 而且要判断实测资料是否绝大多数落在84%及16%概率水平的预测曲线之内(之间).