H2O-CH4-H2-CO2-CO-O2在高温、高压和高密度条件下的状态方程

为了计算高温，高压和高密度流体的热力学性质，提出了一个具有19个参数的维里型状态方程，其中参数与温度间的函数关系采用由Sutherland位能函数导出的维里系数近似式。除临界点附近以外，在已报道的pVT数据所覆盖的大部分超临界区域内，该方程均可适用。用该方程对H2O，CH4，H2，CO2，CO和O2等流体pVT关系的计算结果令人满意，其中pVT上限分别为：91-610GPa,1.6-11.0cm^3/mol,4000-5000K。计算体积的平均偏 小于0．8％，最大偏差小于5．4％。     

离散裂隙渗流方法与裂隙化渗透介质建模

流体渗流模拟的连续介质方法通常适用于多孔地质体，并不一定适用于裂隙岩体，由于裂隙分布及其特征与孔隙差异较大。若流体渗流主要受裂隙的控制，对于一定尺寸的裂隙岩体，多孔介质假设则较难刻划裂隙岩体的渗流特征。离散裂隙渗流方法不但可直接用于模拟裂隙岩体非均质性和各向异性等渗流特征，而且可用其确定所研究的裂隙岩体典型单元体及其水力传导（渗透）张量大小。主要讨论了以下问题：（1）饱和裂隙介质中一般的离散流体渗流模拟；（2）裂隙岩体中的REV（典型单元体）及其水力传导（渗透）张量的确定；（3）利用离散裂隙网络流体渗流模型研究裂隙方向几何参数对水力传导系数和REV的影响；（4）在二维和三维离散裂隙流体渗流模型中对区域大裂隙和局部小裂隙的处理方法。调查结果显示离散裂隙流体渗流数学模型可用来评价不同尺度上的裂隙岩体的水力特征，以及裂隙方向对裂隙化岩体的水力特征有着不可忽视的影响。同时，局部小裂隙、区域大裂隙应当区别对待，以便据其所起的作用及水力特征，建立裂隙化岩体相应的流体渗流模型。     

得尔布干成矿带中段八大关——新峰山成矿地质条件分析

得尔布干成矿带中段八大关—新峰山地区断裂构造发育 ,北东向得尔布干深大断裂是重要的控岩控矿构造。与成矿有关岩浆岩主要为燕山期的中性—中酸性侵入岩及次火山岩。赋矿层位主要有上库力组、塔木兰沟组及谢尔塔拉组。在分析典型矿床的成矿地质特征基础上 ,进一步指出找矿方向     

华北克拉通北部内、外缘中生代流体成矿作用研究--以大井、蔡家营典型矿床为例

华北克拉通北部内、外缘中生代形成的有色金属矿种繁多 ,以外缘富 Sn、内缘富 Mo,Cu、Pb、 Zn遍布内、外缘为主要特色。通过内、外缘典型矿床的研究和综合对比 ,认为外缘的富 Sn是与富 Sn的基底有关 ,内缘的富 Mo是与结晶基底地层和岩浆热液关系密切。而内、外缘 Cu、 Pb、 Zn的成矿与中生代构造转折 ,造成大规模陆相火山活动和成矿流体的卸载密切相关     

川滇密西西比河谷型铅锌矿床成矿流体来源研究:流体Na-Cl-Br体系的证据

应用流体包裹体滤液分析方法，测定了川滇7个MVT铅锌矿床成矿流体的Na,Cl,Br含量，结果表明成矿流体的x(Na/Br)和x(Cl/Br)的平均值分别为185和73，并与高度蒸发浓缩的残留海水的x(Na/Br)和x(Cl/Br)相近。成矿流体的Na,Cl含量呈正相关，在lgNa-lgCl图解中呈线性分布。根据这些事实以及矿床地质特征，认为原始成矿卤水起源于蒸发浓缩的残留海水，原始含矿卤水与富含有机质的大气降水混合导致矿质沉淀而成矿。     

Groundwater-flow modeling in the Yucatan karstic aquifer, Mexico

The current conceptual model of the unconfined karstic aquifer in the Yucatan Peninsula, Mexico, is that a fresh-water lens floats above denser saline water that penetrates more than 40 km inland. The transmissivity of the aquifer is very high so the hydraulic gradient is very low, ranging from 7–10 mm/km through most of the northern part of the peninsula. The computer modeling program AQUIFER was used to investigate the regional groundwater flow in the aquifer. The karstified zone was modeled using the assumption that it acts hydraulically similar to a granular, porous medium. As part of the calibration, the following hypotheses were tested: (1) karstic features play an important role in the groundwater-flow system; (2) a ring or belt of sinkholes in the area is a manifestation of a zone of high transmissivity that facilitates the channeling of groundwater toward the Gulf of Mexico; and (3) the geologic features in the southern part of Yucatan influence the groundwater-flow system. The model shows that the Sierrita de Ticul fault, in the southwestern part of the study area, acts as a flow barrier and head values decline toward the northeast. The modeling also shows that the regional flow-system dynamics have not been altered despite the large number of pumping wells because the volume of water pumped is small compared with the volume of recharge, and the well-developed karst system of the region has a very high hydraulic conductivity. Electronic Publication     

Origin of Gold—Bearing Fluid and Its Initiative Localization Mechanism in Xiadian Gold Deposit,Shandong Province

The composition of quartz inclusions and trace elements in ore indicate that gold-bearing fluid in the Xiadian gold deposit,Shandong Province,stemmed from both mantle and magma,belonging to a composite origin.Based on theoretical analysis and high temperature and high pressure experimental studies,gold-bearing fluid initiative localization mechanism and the forming environment of ore-host rocks are discussed in the present paper.The composite fluid extracted gold from rocks because of its expanding and injecting forces and injecting forces and flew through ore-conducive structures,leading to the breakup of rocks.The generation of ore-host faults and the precipitation of gold-bearing fluid occurred almost simultaneously.This study provides fur-ther information about the relationships between gold ore veins and basic-ultrabasic vein rocks and intermediate vein rocks,the spatial distribution of gold ore veins and the rules governing the migration of ore fluids.     

Constraints on the application of carbon isotope thermometry in high- to ultrahigh-temperature metamorphic terranes

Nine marble horizons from the granulite facies terrane of southern India were examined in detail for stable carbon and oxygen isotopes in calcite and carbon isotopes in graphite. The marbles in Trivandrum Block show coupled lowering of δ¹³C and δ¹⁸O values in calcite and heterogeneous single crystal δ¹³C values (? 1 to ? 10‰) for graphite indicating varying carbon isotope fractionation between calcite and graphite, despite the granulite facies regional metamorphic conditions. The stable isotope patterns suggest alteration of δ¹³C and δ¹⁸O values in marbles by infiltration of low δ¹³C–δ¹⁸O‐bearing fluids, the extent of alteration being a direct function of the fluid‐rock ratio. The carbon isotope zonation preserved in graphite suggests that the graphite crystals precipitated/recrystallized in the presence of an externally derived CO₂‐rich fluid, and that the infiltration had occurred under high temperature and low f_O2 conditions during metamorphism. The onset of graphite precipitation resulted in a depletion of the carbon isotope values of the remaining fluid+calcite carbon reservoir, following a Rayleigh‐type distillation process within fluid‐rich pockets/pathways in marbles resulting in the observed zonation. The results suggest that calcite–graphite thermometry cannot be applied in marbles that are affected by external carbonic fluid infiltration. However, marble horizons in the Madurai Block, where the effect of fluid infiltration is not detected, record clear imprints of ultrahigh temperature metamorphism (800–1000 °C), with fractionations reaching <2‰. Zonation studies on graphite show a nominal rimward lowering δ¹³C on the order of 1 to 2‰. The zonation carries the imprint of fluid deficient/absent UHT metamorphism. Commonly, calculated core temperatures are > 1000 °C and would be consistent with UHT metamorphism.     

An automated strategy for calculation of phase diagram sections and retrieval of rock properties as a function of physical conditions

We formulate an algorithm for the calculation of stable phase relations of a system with constrained bulk composition as a function of its environmental variables. The basis of this algorithm is the approximate representation of the free energy composition surfaces of solution phases by inscribed polyhedra. This representation leads to discretization of high variance phase fields into a continuous mesh of smaller polygonal fields within which the composition and physical properties of the phases are uniquely determined. The resulting phase diagram sections are useful for understanding the phase relations of complex metamorphic systems and for applications in which it is necessary to establish the variations in rock properties such as density, seismic velocities and volatile‐content through a metamorphic cycle. The algorithm has been implemented within a computer program that is general with respect to both the choice of variables and the number of components and phases possible in a system, and is independent of the structure of the equations of state used to describe the phases of the system.     

Physical properties of ultrahigh-pressure metamorphic rocks from the Sulu terrain, eastern central China: implications for the seismic structure at the Donghai (CCSD) drilling site

About 30 samples representing major lithologies of Sulu ultrahigh-pressure (UHP) metamorphic rocks were collected from surface exposures and exploration wells, and compressional (Vp) and shear wave (Vs) velocities and their directional dependence (anisotropy) were determined over a range of constant confining pressures up to 600 MPa and temperatures ranging from 20 to 600 °C. Samples range in composition from acidic to ultramafic. P- and S-wave velocities measured at 600 MPa vary from 5.08 to 8.64 km/s and 2.34 to 4.93 km/s, respectively. Densities are in the range from 2.60 to 3.68 g/cm³. To make a direct tie between seismic measurements (refraction and reflection) and subsurface lithologies, the experimental velocity data (corresponding to shallow depths) were used to calculate velocity profiles for the different lithologies and profiles of reflection coefficients at possible lithologic interfaces across the projected 5000-m Chinese Continental Scientific Drilling Program (CCSD) crustal segment. Comparison of calculated in situ velocities with respective intrinsic velocities suggests that the in situ velocities at shallow depths are lowered by an increased abundance of open microcracks. The strongly reflective zone beneath the Donghai drill site can be explained by the impedance contrasts between the different lithologies. Contacts between eclogite/peridotite and felsic rocks (gt-gneiss, granitic gneiss), in particular, may give rise to strong seismic reflections. In addition, shear-induced (lattice preferred orientation (LPO)-related) seismic anisotropy can increase reflectivity. For the explanation of the high velocity bodies (>6.4 km/s) around 1000 m and below 3200-m depth, large proportions of eclogite/peridotite (about 40 and 30 vol.%, respectively) are needed.