滇西北保山核桃坪铅锌矿床成矿流体来源的同位素地球化学证据

滇西北核桃坪铅锌矿是保山地块内沿澜沧江断裂带发育的一个大型铅锌多金属矿床,矿体严格受近SN向的断裂带及层间破碎带控制,赋存于上寒武统核桃坪组和沙河厂组大理岩化灰岩及矽卡岩中。文章在分析矿床地质特征的基础上,通过对V1矿体的氢、氧、碳、硫、硅和铅同位素地球化学的示踪研究,探讨了成矿溶液中水、碳、硅、硫和铅的来源以及成矿流体的演化问题。研究表明,V1矿体中矿石的石英δ18O水值变化范围为6·1‰~7·6‰,平均为6·7‰,其δDV-SMOW值变化范围为-100‰~-108‰,平均为-104‰,说明成矿流体主要来自深部岩浆分异水,并在后期成矿作用过程中有地层建造水和大气降水的加入;矿石中方解石的δ13CV-PDB和δ18O水值分别为-6·6‰~-5·9‰和5·0‰~5·2‰,反映成矿流体的碳、氧具有深部流体与围岩发生交代反应的特征。硅、硫和铅同位素具深部岩浆或地幔来源的特点,但遭受后期壳源物质的混染。据此提出该矿床是与深部岩浆热液或幔源流体贯入有关的构造热液脉型矿床,其成矿作用受控于深部过程的驱动和控制,可为指导地质找矿提供依据。     

胶东三甲金矿床流体包裹体特征

三甲金矿是胶东牟平-乳山金成矿带内重要的石英脉型金矿,金主要产于黄铁矿和多金属硫化物石英脉中。流体包裹体研究表明,三甲金矿蚀变岩石和各成矿阶段金矿石中的流体包裹体主要有三种类型：H2O-CO2包裹体、富CO2包裹体和H2O溶液包裹体。早期乳白色石英中主要赋存原生的H2O-CO2包裹体;成矿期黄铁矿石英脉和多金属硫化物石英脉中的富CO2包裹体主要为原生,随机分布,气液比变化较大,常与早期H2O溶液包裹体共生且均一温度接近,显示不混溶流体包裹体组合特征;在成矿晚期的石英和方解石中主要发育原生H2O溶液包裹体。显微测温结果显示,成矿前（第1阶段）H2O-CO2包裹体的完全均一温度（Tb．TOT,至液相）为280℃至416℃,成矿期（第Ⅱ和Ⅲ阶段）富CO2包裹体的完全均一温度为210—330℃,同期的H2O溶液包裹体均一温度为253～377℃,成矿后（第Ⅳ阶段）H2O溶液包裹体的均一温度为176—207℃。成矿流体为低盐度的CO2-H2O-NaCl型热液,成矿应力场转变导致的流体减压沸腾作用可能是三甲金矿金沉淀成矿的主要原因。     

Subsurface Water Masses in the Central North Pacific Transition Region: The Repeat Section along the 18° Meridian

A repeat hydrographic section has been maintained over two decades along the 180° meridian across the subarctic-subtropical transition region. The section is naturally divided into at least three distinct zones. In the Subarctic Zone north of 46°N, the permanent halocline dominates the density stratification, supporting a subsurface temperature minimum (STM). The Subarctic Frontal Zone (SFZ) between 42°–46°N is the region where the subarctic halocline outcrops. To the south is the Subtropical Zone, where the permanent thermocline dominates the density stratification, containing a pycnostad of North Pacific Central Mode Water (CMW). The STM water colder than 4°C in the Subarctic Zone is originated in the winter mixed layer of the Bering Sea. The temporal variation of its core temperature lags 12–16 months behind the variations of both the winter sea surface temperature (SST) and the summer STM temperature in the Bering Sea, suggesting that the thermal anomalies imposed on the STM water by wintertime air-sea interaction in the Bering Sea spread over the western subarctic gyre, reaching the 180° meridian within a year or so. The CMW in this section originates in the winter mixed layer near the northern edge of the Subtropical Zone between 160°E and 180°. The CMW properties changed abruptly from 1988 to 1989; its temperature and salinity increased and its potential density decreased. It is argued that these changes were caused by the climate regime shift in 1988/1989 characterized by weakening of the Aleutian Low and the westerlies and increase in the SST in the subarctic-subtropical transition region. This revised version was published online in August 2006 with corrections to the Cover Date.     

Integration of local–global upscaling and grid adaptivity for simulation of subsurface flow in heterogeneous formations

We propose a methodology, called multilevel local–global (MLLG) upscaling, for generating accurate upscaled models of permeabilities or transmissibilities for flow simulation on adapted grids in heterogeneous subsurface formations. The method generates an initial adapted grid based on the given fine-scale reservoir heterogeneity and potential flow paths. It then applies local–global (LG) upscaling for permeability or transmissibility [7], along with adaptivity, in an iterative manner. In each iteration of MLLG, the grid can be adapted where needed to reduce flow solver and upscaling errors. The adaptivity is controlled with a flow-based indicator. The iterative process is continued until consistency between the global solve on the adapted grid and the local solves is obtained. While each application of LG upscaling is also an iterative process, this inner iteration generally takes only one or two iterations to converge. Furthermore, the number of outer iterations is bounded above, and hence, the computational costs of this approach are low. We design a new flow-based weighting of transmissibility values in LG upscaling that significantly improves the accuracy of LG and MLLG over traditional local transmissibility calculations. For highly heterogeneous (e.g., channelized) systems, the integration of grid adaptivity and LG upscaling is shown to consistently provide more accurate coarse-scale models for global flow, relative to reference fine-scale results, than do existing upscaling techniques applied to uniform grids of similar densities. Another attractive property of the integration of upscaling and adaptivity is that process dependency is strongly reduced, that is, the approach computes accurate global flow results also for flows driven by boundary conditions different from the generic boundary conditions used to compute the upscaled parameters. The method is demonstrated on Cartesian cell-based anisotropic refinement (CCAR) grids, but it can be applied to other adaptation strategies for structured grids and extended to unstructured grids.     

Accurate local upscaling with variable compact multipoint transmissibility calculations

We propose a new single-phase local upscaling method that uses spatially varying multipoint transmissibility calculations. The method is demonstrated on two-dimensional Cartesian and adaptive Cartesian grids. For each cell face in the coarse upscaled grid, we create a local fine grid region surrounding the face on which we solve two generic local flow problems. The multipoint stencils used to calculate the fluxes across coarse grid cell faces involve the six neighboring pressure values. They are required to honor the two generic flow problems. The remaining degrees of freedom are used to maximize compactness and to ensure that the flux approximation is as close as possible to being two-point. The resulting multipoint flux approximations are spatially varying (a subset of the six neighbors is adaptively chosen) and reduce to two-point expressions in cases without full-tensor anisotropy. Numerical tests show that the method significantly improves upscaling accuracy as compared to commonly used local methods and also compares favorably with a local–global upscaling method.     

Evaluation of change in subsurface thermal environment due to groundwater flow in the Tokyo Lowland,Japan

Information on the distribution of subsurface temperature and hydraulic heads at 24 observation wells in and around the Tokyo Lowland, the eastern part of the Tokyo Metropolitan area, were examined to make clear the relationship between groundwater and the subsurface thermal environment in the urban area. Minimums in temperature–depth logs due to subsurface temperature increasing at shallow parts were recognized in 21 wells. This fact shows subsurface temperature is affected by ground surface warming in almost all of this area. Deeper than minimums, where the effects of surface warming became relatively small, regional variation is observed as follows: high temperatures are shown in the central part to the southern part, and low temperatures shown in the inland to eastern part. The high temperature area corresponds to an area where the lower boundary of groundwater flow is relatively shallow. This area corresponds also to an area with severe land subsidence resulting from excessive groundwater pumping. It is considered that this high temperature area is formed by the effects of upward groundwater flow affected by hydrogeological conditions and pumping. On the other hand, a comparison between past data (1956–1967) and present data (2001–2003) revealed widespread decreasing temperature in the inland area. This is explained by downward groundwater flow based on an analysis of temperature–depth logs. This fact suggests that subsurface temperature is not only increasing from the effects of surface warming but also decreasing from the effects of groundwater environment change due to pumping.     

幔源CO_2释出机理、脱气模式及成藏机制研究进展

针对幔源CO2如何从地幔岩浆中脱出并进入沉积地层中形成CO2气藏聚集这一关键问题,总结了国内外研究进展和前缘方向。研究表明,地幔深部的碱性玄武岩浆和碱性岩浆才是深部流体和CO2等挥发份大量赋存、渗滤和释出的场所。浅成侵入岩、次火山岩和火山通道等是CO2释放和聚集的有利位置,岩浆期后和岩浆衰弱期的热液活动阶段是CO2大量释放和聚集的有利时期。幔源CO2进入沉积盆地中具有3种脱气模式,即沿岩石圈断裂直接脱气模式、热流底辟体脱气模式和壳内岩浆房-基底断裂组合脱气模式。CO2的固有物化性质决定其运移相态多样,具有运移和聚集过程同步的特征。只有在满足大量的化学消耗及地层水或原油的溶解和耗散之后才能形成CO2有效聚集。幔源CO2成藏和分布主要受岩浆气源体和气源断裂体系的控制。今后,在超临界CO2及其对油气运移聚集的作用、CO2与深大断裂及火山岩的关系、CO2脱气运移机制、CO与常规烃类油气的耦合差异成藏机制等方面仍需要进一步的研究和探索。     

Fluid sources for the La Guitarra epithermal deposit (Temascaltepec district, Mexico): Volatile and helium isotope analyses in fluid inclusions

The La Guitarra deposit (Temascaltepec district, South-Central Mexico), belongs to the low/intermediate sulfidation epithermal type, has a polymetallic character although it is currently being mined for Ag and Au. The mineralization shows a polyphasic character and formed through several stages and sub-stages (named I, IIA, IIB, IIC, IID, and III). The previous structural, mineralogical, fluid inclusion and stable isotope studies were used to constrain the selection of samples for volatile and helium isotope analyses portrayed in this study. The N₂/Ar overall range obtained from analytical runs on fluid inclusion volatiles, by means of Quadrupole Mass Spectrometry (QMS), is 0 to 2526, and it ranges 0 to 2526 for stage I, 0 to 1264 for stage IIA, 0 to 1369 for stage IIB, 11 to 2401 for stage IIC, 19 to 324 for stage IID, and 0 to 2526 for stage III. These values, combined with the CO₂/CH₄ ratios, and N₂-He-Ar and N₂-CH₄-Ar relationships, suggest the occurrence of fluids from magmatic, crustal, and shallow meteoric sources in the forming epithermal vein deposit. The helium isotope analyses, obtained by means of Noble Gas Mass Spectrometry, display R/Ra average values between 0.5 and 2, pointing to the occurrence of mantle-derived helium that was relatively diluted or “contaminated” by crustal helium. These volatile analyses, when correlated with the stable isotope data from previous works and He isotope data, show the same distribution of data concerning sources for mineralizing fluids, especially those corresponding to magmatic and crustal sources. Thus, the overall geochemical data from mineralizing fluids are revealed as intrinsically consistent when compared to each other.The three main sources for mineralizing fluids (magmatic, and both deep and shallow meteoric fluids) are accountable at any scale, from stages of mineralization down to specific mineral associations. The volatile and helium isotope data obtained in this paper suggest that the precious metal-bearing mineral associations formed after hydrothermal pulses of predominantly oxidized magmatic fluids, and thus it is likely that precious metals were carried by fluids with such origin. Minerals from base-metal sulfide associations record both crustal and magmatic sources for mineralizing fluids, thus suggesting that base metals could be derived from deep leaching of crustal rocks. At the La Guitarra epithermal deposit there is no evidence for an evolution of mineralizing fluids towards any dominant source. Rather than that, volatile analyses in fluid inclusions suggest that this deposit formed as a pulsing hydrothermal system where each pulse or set of pulses accounts for different compositions of mineralizing fluids.The positive correlation between the relative content of magmatic fluids (high N₂/Ar ratios) and H₂S suggests that the necessary sulfur to carry mostly gold as bisulfide complexes came essentially from magmatic sources. Chlorine necessary to carry silver and base metals was found to be abundant in inclusion fluids and although there is no evidence about its source, it is plausible that it may come from magmatic sources as well.     

Albite vein formation during exhumation of high-pressure terranes: a case study from alpine Corsica

The formation of late‐stage veins can yield valuable information about the movement and composition of fluids during uplift and exhumation of high‐pressure terranes. Albite veins are especially suited to this purpose because they are ubiquitously associated with the greenschist facies overprint in high‐pressure rocks. Albite veins in retrogressed metabasic rocks from high‐pressure ophiolitic units of Alpine Corsica (France) are nearly monomineralic, and have distinct alteration haloes composed of actinolite + epidote + chlorite + albite. Estimated P–T conditions of albite vein formation are 478 ± 31 °C and 0.37 ± 0.14 GPa. The P–T estimates and petrographic constraints indicate that the albite veins formed after the regional greenschist facies retrogression, in response to continued decompression and exhumation of the terrane. Stable isotope geochemistry of the albite veins, their associated alteration haloes and unaltered hostrocks indicates that the vein‐forming fluid was derived from the ophiolite units and probably from the metabasalts within each ophiolite slice. That the vein‐forming fluid was locally derived means that a viable source of fluid to form the veins was retained in the rocks during high‐pressure metamorphism, indicating that the rocks did not completely dehydrate. This conclusion is supported by the observation of abundant lawsonite at the highest metamorphic grades. Fluids were liberated during retrogression via decompression dehydration reactions such as those that break down hydrous high‐pressure minerals like lawsonite. Albite precipitation into veins is sensitive to the solubility and speciation of Al, which is more pressure sensitive than other factors which might influence albite vein formation such as silica saturation or Na:K fluid ratios. Hydraulic fracturing in response to fluid generation during decompression was probably the main mechanism of vein formation. The associated pressure decrease with fracturing and fluid decompression may also have been sufficient to change the solubility of Al and drive albite precipitation in fracture systems.     

Computer simulation studies of Crystallization under confinement conditions

Crystallization under confinement conditions is a very important process in geochemistry and geophysics. Computer simulations of fluids in nanometer scale pore spaces can provide a unique microscopic insight into the structure, dynamics and forces arising from the crystallization process. We discuss in this paper molecular dynamics computer simulations of crystallization in pores of nanometer dimensions. The crystallization pressure due to the freezing of a model of Argon in a nanopore is computed using molecular dynamics simulations. We also investigate the influence of pore geometry in determining the dynamics of confined fluids, as well as mass separation in binary mixtures. It turns out that the pore geometry reveals itself as an important variable, leading to 1) new mechanisms for fast diffusion in confined spaces, and 2) accumulation of solute in specific regions inside the pore.