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

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

滇西北地区富碱岩体(脉)地质学及岩石化学特征

喜山期富碱岩体广布于滇西北地区 ,通过对其地质学和岩石化学特征的研究 ,认为富碱岩体 (脉 )由正长 (斑 )岩、二长斑岩、花岗斑岩及相关的火山岩组成 ,具有多次活动的特征 ,并形成区内富有特色的岩浆岩系列 ;该系列岩浆岩以富碱高钾为突出特点 ,具有正常的岩浆演化趋势 ,不同岩带岩石地质学和岩石化学方面有区域性变化规律。     

熊耳山地区太古宙超镁铁质岩的地球化学特征及成因探讨

熊耳山地区太古宙超镁铁质岩在空间分布上具有不均一性、成带性和多方向性特征 ,岩石以橄榄岩、辉石岩、角闪石岩为主 ,蚀变强烈。岩石贫钙、镁 ,低铝、K2 O Na2 O,稀土总量较高 ,轻重稀土分馏显著。地球化学特征表明大部分超镁铁质岩体为岩浆侵入成因     

上扬子台地晚震旦世灯影组中葡萄状-雪花状白云岩的成因意义

葡萄状 -雪花状白云岩是白云岩形成之后经淡水淋滤作用而成 ,葡萄内或雪花周围的有机物质不是原生藻 ,而是蓝绿藻经细菌腐解后形成的细菌集合体。该套白云岩的形成与沉积环境关系密切 ,与藻类多少无关。葡萄状白云岩富含有机物质 ,而雪花状白云岩孔洞发育 ,因此葡萄状 -雪花状白云岩既是生油层 ,又是良好的储集层 ,对其重新认识和定位 ,必将为寻找深部油气资源提供极为重要场所     

成矿流体演化与成矿物理化学

成矿流体是富含挥发分 (CO2 、CH4等 )是具有较高含盐度的特殊地质流体。本文讨论了在流体演化过程中挥发分的来源 ,指出主要来自水岩作用、有机质分解及地幔去气和岩浆 ;碱金属及卤素同样具有多来源的性质 ,以海水、含盐系淋滤、建造水为主要来源 ,含盐系重熔可以形成富含碱金属的成矿流体。流体演化过程中氢氧同位素、硫同位素的分馏主要与温度、水岩比值或硫源丰度有关。一个重要的结论是 ,成矿流体的形成主要与地质作用有关 ,而流体来源是次要的。海底热水流体的地球化学特征以高δ3 4 S值、中稀土富集及正铕异常为特征。本文总结了热水流体成矿物理化学条件 ,指出水热流体物相点 :1) 10 80℃ ,7.5× 10 8Pa水溶液与硅酸岩熔浆分熔点 ;2 )水溶液的第二个临界点是气水溶液的超临界点 (374.15℃ ,2 .2 1× 10 7Pa) ;3)水溶液的沸点 (≥ 10 0℃ ,≥ 1× 10 5Pa) ;4)水溶液的冰点 (≤ 0℃ ,1× 10 5Pa) ;5 )H2 O CO2 体系的不混溶温度点 (2 6 6℃ ,2 .15×10 8Pa[1 3 ] 等是重要的成矿相变点。     

青藏高原北羌塘半岛湖新生代粗面玄武岩橄榄石电子探针和激光探针分析

利用电子探针和激光探针剥蚀系统（LA-ICP-MS），对北羌塘新第三纪粗面玄武岩中的橄榄石主元素和微量、稀土元素进行了系统分析。结果表明，本区橄榄石Fo平均为88，属贵橄榄石种属。相对富集Ni、Co和重稀土，而强烈亏损轻稀土及Rb、Sr、Ba、Zr等大离子亲石元素。其稀土元素配分型配分型式与粗面玄武岩全岩稀土配分型式呈互补状态。     

A non-linear stress-stiffness model for geomaterials at small to intermediate strains

A double exponential fitting model (DEFM) capable of expressing the non-linear stress-stiffness relationship of geomaterials has been proposed by Shibuya et al. (1997). The model comprises two material constants; the elastic stiffness at very small strains and the strength, together with other free parameters to determine the complete stress-stiffness relationship. In this paper, the capability of the original function used for DEFM in simulating the tangent stiffness-stress relationship of geomaterials is first discussed. Second, the methods for determining the free model parameters, as well as its conversion to obtain a stress-strain relationship are proposed. The applicability of DEFM to predicting non-linear stress-stiffness relationship is examined in detail in a total of forty-nine fitting cases of compression test data on sedimentary rock, artificial soft rock and soft clay. It is found that the DEFM is effective in expressing the non-linear stress-stiffness relationship of various kinds of geomaterials at small to intermediate strains, say less than 0.5%. The superiority of this model compared to other fitting models currently in use is also demonstrated in some of the fitting cases.     

An experimental study of grain scale melt segregation mechanisms in two common crustal rock types

Creation of pathways for melt to migrate from its source is the necessary first step for transport of magma to the upper crust. To test the role of different dehydration‐melting reactions in the development of permeability during partial melting and deformation in the crust, we experimentally deformed two common crustal rock types. A muscovite‐biotite metapelite and a biotite gneiss were deformed at conditions below, at and above their fluid‐absent solidus. For the metapelite, temperatures ranged between 650 and 800 °C at P_c=700 MPa to investigate the muscovite‐dehydration melting reaction. For the biotite gneiss, temperatures ranged between 850 and 950 °C at P_c=1000 MPa to explore biotite dehydration‐melting under lower crustal conditions. Deformation for both sets of experiments was performed at the same strain rate (ε.) 1.37×10^?5 s^?1. In the presence of deformation, the positive ΔV and associated high dilational strain of the muscovite dehydration‐melting reaction produces an increase in melt pore pressure with partial melting of the metapelite. In contrast, the biotite dehydration‐melting reaction is not associated with a large dilational strain and during deformation and partial melting of the biotite gneiss melt pore pressure builds more gradually. Due to the different rates in pore pressure increase, melt‐enhanced deformation microstructures reflect the different dehydration melting reactions themselves. Permeability development in the two rocks differs because grain boundaries control melt distribution to a greater extent in the gneiss. Muscovite‐dehydration melting may develop melt pathways at low melt fractions due to a larger volume of melt, in comparison with biotite‐dehydration melting, generated at the solidus. This may be a viable physical mechanism in which rapid melt segregation from a metapelitic source rock can occur. Alternatively, the results from the gneiss experiments suggest continual draining of biotite‐derived magma from the lower crust with melt migration paths controlled by structural anisotropies in the protolith.     

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.     

Fluid permeability of sedimentary rocks in a complete stress–strain process

The permeability of sedimentary rocks during triaxial compression tests was investigated to relate it to the complete strain–stress process. It was found that the permeability was not constant, but varied with the stress and strain states in the rocks. Prior to the peak strength, the permeability decreases with increasing load. A dramatic increase in permeability occurs during the strain softening period. In the present study, in situ measurements of fluid flow and pressure in floor strata was carried out in a double longwall mining face in the Yangzhuang colliery. These measurements show that both the strata pressure and the position with respect to the mining face influence the hydrogeologic properties. The permeability increased in the floor strata behind the mining face because those mining induced fractures opened as the strata pressure decreased. To better understand this change in hydraulic behavior around the mining faces, 3-D numerical modeling was carried out. The model provides the general picture of the stress distribution and failure zone both in the floor and roof strata. The field and model results demonstrate the importance of changes in the stress and strain states on the hydrogeology of a site.