New heat flow data from the immediate vicinity of the Kola super-deep borehole: Vertical variation in heat flow confirmed and attributed to advection

We present new heat flow values and other geothermal data in the upper crystalline crust in the immediate vicinity of the 12.4-km deep Kola super-deep borehole, NW Russia. Our results show a systematic vertical increase in geothermal gradient and heat flow density as deep as we could measure (1.6 km). Our results confirm earlier results on vertical heat flow trends of in the uppermost part of the Kola super-deep hole, and imply that the thermal regime is not in steady-state conductive conditions. In an area of 3-km × 5-km measurements were performed in 1–2-km deep boreholes surrounding the Kola super-deep hole and on core samples from these holes. Temperature logs are available from 36 holes. Core data exists from 23 boreholes with a total length of 11.5 km at a vertical resolution of 10 m. We carried out a very detailed study on thermal conductivity with regard to anisotropy, inhomogeneity and temperature dependence. Tensor components of thermal conductivity were determined on 1375 core samples from 21 boreholes in 3400 measurements. Additionally, we measured specific heat capacity, heat generation rate, density, porosity, and permeability on selected subsets of core samples. Heat flow from 19 boreholes varies between 31 and 45 mW m⁻² with an average value of 38 mW m⁻². In most boreholes the vertical heat flow profiles show a considerable variation with depth. This is consistent with observations in the upper part of the Kola super-deep borehole. We conclude that this variation is not caused by technical operations but reflects a natural process. It is considered to be due to a combination of advective, structural and paleoclimatic effects. Preliminary 3-D numerical modeling of heat and flow in the study area provides an indication of relative contributions of each of these factors: advective heat transfer turns out to have a major influence on the vertical variation of heat flow, although transient changes in surface temperature may also cause a significant variation. Heterogeneity of the rocks in the study area is less important.     

Review and classification of structural controls on fluid flow during regional metamorphism

Mechanisms for kilometre-scale, open-system fluid flow during regional metamorphism remain problematic. Debate also continues over the degree of fluid flow channellization during regional metamorphism, and the mechanisms for pervasive fluid flow at depth. The requirements for pervasive long-distance fluid flow are an interconnected porosity and a large regional gradient in fluid pressure and hydraulic head (thermally or structurally controlled) that dominates over local perturbations in hydraulic head due to deformation. In contrast, dynamic or transient porosity interconnection and fluid flow accompanying deformation of heterogeneous rock suites should result in moderately to strongly channellized flow at a range of scales, of which there are many examples in the literature. Classification of fluid flow types based on scale and degree of equilibration between fluid and rock, wallrock permeability, and mode of fluid transport contributes to an understanding of key factors that control fluid flow. Closed-system fluid behaviour, with restricted fluid flow in microcracks or cracks and limited fluid–rock interaction, occurs over a range of strains and crustal depths, but requires low permeabilities and/or small fluid fluxes. Long-distance, open-system fluid flow in channels is favoured in heterogeneous rocks at high strains, moderate (but variable) permeabilities, and moderate to high fluid fluxes. Long-distance, broad, pervasive fluid flow during regional metamorphism requires that the rocks are not accumulating high strains and have high permeabilities, low permeability contrasts, and high fluid fluxes. The ideal situation for such fluid flow is in situations where the rocks are undergoing stress relaxation immediately after a major deformation phase. In the mid-crust, fairly specific conditions are thus required for pervasive fluid flow. During active orogenesis, structurally controlled fluid flow (with focused open systems surrounding regions of closed-system behaviour) predominates in most, but not all, regional metamorphic situations, at a range of scales.     

Determination of the elastic modulus set of foliated rocks from ultrasonic velocity measurements

Ultrasonic measurements of compressional and shear wave velocities under hydrostatic pressure up to 70 MPa were carried out on cylindrical specimens cored across and along the foliation planes. Our measurements revealed that the foliation of the metamorphic rocks induces a clear velocity anisotropy between two orthogonal directions; faster along the foliation plane and slower across the plane in most rock types. All velocity components monotonically increase with the confining pressure, probably due to the closure of microcracks distributed in rock specimens. We determined the complete set of dynamic moduli of foliated metamorphic rocks with two assumptions; transverse isotropy due to the foliation and ellipsoidal seismic energy propagation from a point source. The calculated elastic moduli referring to different directions could be valuable for the design of various engineering structures in planar textured rock mass.     

多尺度岩石渗透特性的试验研究

岩体赋存于一定地质环境中,经过漫长的地质构造作用,形成由连续性状的岩石材料和非连续性状的各种地质结构面组成的复杂地质体。渗透率是影响这种复杂岩体渗透特性的重要参数之一,渗透率选择的正确与否直接关系到裂隙岩体渗流分析结果的可靠性与工程稳定性。采用室内试验的方法,结合达西定律确定裂隙岩体的渗透率。通过自行研发的5种尺寸压力室测定不同尺寸岩石试样的渗透率,得到由不同裂隙、孔洞组合的岩石试样在稳定流、非稳定流下的渗透率,并由此分析得到多尺度岩石的渗透特性。从多组试验结果中分析可得,随着岩石试样的尺寸增加,岩石的渗流通道增多,从而渗透率增大,但增大到一定程度后会趋于稳定,这一渗透规律对实际工程有指导意义。     

The controls and selected peculiarities of the origin of contact metamorphic zonation

At least sis or perhaps seven types of contact metamorphism may be distinguished in nature (see Table 1). They differ from each other by a set of metamorphic facies in the exocontact aureoles, also by thermodynamic conditions of metamorphism. The manifestation of some types of metamorphism depends chiefly on magma temperature and composition, on the initial temperature of the country rocks (prior to contact metamorphism) and on the depth of the intrusive formation. The movement of magma through the intrusive channel (chamber) and the kinetic peculiarities of metamorphism exert additional influence on metamorphic conditions.Since the temperature elevation with depth corresponds in a general way to the increase in the lithostatic pressure, the maximum temperature levels attained at the direct intrusive contact must differ for various pressure levels (or depths). The two series of contact metamorphic rooks may be distinguished by their pressure: 1. common hornfelsic rock low-pressure series (types of metamorphism 1–4, Table 1), and 2. thermal transformed gneissous rock series metamorphosed under moderately high pressure (types 5–6, Table 1). The initial temperature of the country rock is low in the first series. It may apparently be discounted in certain cases. The usual non-abyssal contact low-pressure metamorphism (the proper contact metamorphism) depends essentially on magma temperature and composition, while in the second series, the initial temperature of the country rocks may be rather high. This and the effect of the intrusive heat and high lithostatic pressure results in rock transformation under conditions of the regional metamorphism of the facies (under moderately high pressure).Since superposition of the local temperature field on a regional high-temperature field in the neighbourhood of intrusives (in deep conditions) results in an increase in the thickness of the contact aureole, Ingersoll's criterion (the ratio of the thickness of the aureole to the thickness of the intrusive body) may be helpful for distinguishing abyssal from non-abyssal contact metamorphism. The values of this criterion cannot exceed 0.2–0.3 for the non-abyssal contact metamorphism.     

不同温压下岩石弹性波速度、衰减及各向异性与组构的关系

结合岩石组构分析 ,阐述了岩石弹性波传播速度和衰减以及它们的各向异性与岩石组构之间的关系。在不同温压条件下对具有很强晶格优选方位的岩石样品的研究表明 ,随着围压的增加 ,波速和Q值均增大 ,但是在相互正交的 3个方向上 (垂直或平行于层理面及线理方向 )增大的速度并不相同 ,这与微裂隙的逐渐闭合密切相关。观测到的波速和Q值的各向异性具有不同的形成机理 ,波速各向异性主要与定向分布的微裂隙和主要矿物的晶格优选方位等构造因素有关 ;高围压下Q值各向异性与速度各向异性正好相反 ,可能是由于定向排列的矿物晶体沿不同方向其边界之间接触程度不同造成的。对岩石组构的研究不仅可以揭示岩体的变形机制、变形的动力学过程及其有关的热力学信息 ,还可以对宏观岩石的各种物理性质 ,尤其是力学特性 ,从微观机理上加以解释。文中特别强调了岩石组构分析对研究岩石物理性质的各向异性具有十分重要的意义。     

Combining compositional zoning and foliation intersection axes (FIAs) in garnet to quantitatively determine early P-T-t paths in multiply deformed and metamorphosed schists: north central Massachusetts,USA

Quantitative compositional and microstructural analysis of garnet porphyroblasts in kyanite–staurolite–garnet grade rocks from the northeastern flank of the Pelham dome, north central Massachusetts, distinguishes the effects of Acadian deformation and metamorphism from extensive overprinting Alleghanian shearing. The P–T conditions and the metamorphic path during the Acadian were determined using samples preserving well defined stages in a lengthy tectonic history revealed by a succession of five foliation intersection axis trends preserved within porphyroblasts (FIAs). This Acadian succession extends at least 120 km to the north into rocks where no evidence has been found of an Alleghanian overprint. For each sample where garnet first nucleated during one of these stages in the tectonic history, the PT of core growth was determined by plotting the intersection of the Mn, Fe and Ca isopleths calculated for the core composition on a P–T pseudosection for that sample using THERMOCALC. Combining the PT data from all these samples indicates that the temperature and pressure increased throughout Acadian orogenesis, causing episodic garnet growth. During the Alleghanian, locally intense shearing, especially against the margin of the Pelham dome, formed the dominant schistosity, which truncated all foliations defined by inclusion trails in porphyroblasts and obliterated all remains of Acadian deformation and metamorphism in the rock matrix. Shearing was accompanied by near complete homogenization of the compositional zoning in garnet porphyroblasts and an associated apparent increase in the temperature of the matrix to 700°C in those rocks lying directly adjacent to the Pelham dome, and resulted from the rocks of the Northfield syncline being thrust a large distance southwards over the gneisses in the dome.     

Directional dependence of P- and S-wave propagation and polarization in foliated rocks from the Kola superdeep well: Evidence from laboratory measurements and calculations based on TOF neutron diffraction

We have measured P- and S-wave velocities on two amphibolite and two gneiss samples from the Kola superdeep borehole as a function of pressure (up to 600 MPa) and temperature (up to 600 °C). The velocity measurements include compressional (Vp) and shear wave velocities (Vs1, Vs2) propagating in three orthogonal directions which were in general not parallel to inherent rock symmetry axes or planes. The measurements are accompanied by 3D-velocities calculations based on lattice preferred orientation (LPO) obtained by TOF (Time Of Flight) neutron diffraction analysis which allows the investigation of bulk volumes up to several cubic centimetres due to the high penetration depth of neutrons. The LPO-based numerical velocity calculations give important information on the different contribution of the various rock-forming minerals to bulk elastic anisotropy and on the relations of seismic anisotropy, shear wave splitting, and shear wave polarization to the structural reference frame (foliation and lineation). Comparison with measured velocities obtained for the three propagation directions that were not in accordance with the structural frame of the rocks (foliation and lineation) demonstrate that for shear waves propagating through anisotropic rocks the vibration directions are as important as the propagation directions. The study demonstrates that proper measurement of shear wave splitting by means of two orthogonal polarized sending and receiving shear wave transducers is only possible when their propagation and polarization directions are parallel and normal to foliation and lineation, respectively.     

Three-dimensional discrete element simulation of indirect tensile behaviour of a transversely isotropic rock

This paper presents the development of a three-dimensional discrete element model using flat-joint and smooth-joint contact models to investigate the effect of anisotropy on the tensile behaviour of slate, a transversely isotropic rock, under Brazilian testing from both macro and microscales. The effect of anisotropy is further realised by exploring the influence of foliation orientations (β and ψ) on the tensile strength, fracture pattern, microcracking and stress distribution of the transversely isotropic rock. The variation of tensile strength with foliation orientation is presented. The cross-weak-plane fracture growth observed in laboratory is reproduced, and the criterion for which to form is also given from the aspect of foliation orientation. Furthermore, the proportional variations of microcracks well account for the effects of foliation orientation on the tensile strength and failure pattern. Finally, it is found that the existence of weak planes increases both the heterogeneity and the anisotropy of stress distributions within the transversely isotropic rock, with the degree of influence varying with the foliation orientation.     

The initiation and development of metamorphic foliation in the Otago Schist, Part 1: competitive oriented growth of white mica

The 3D shape, size and orientation data for white mica grains sampled along two transects of increasing metamorphic grade in the Otago Schist, New Zealand, reveal that metamorphic foliation, as defined by mica shape‐preferred orientation (SPO), developed rapidly at sub‐greenschist facies conditions early in the deformation history. The onset of penetrative strain metamorphism is marked by the rapid elimination of poorly oriented large clastic mica in favour of numerous new smaller grains of contrasting composition, higher aspect ratios and a strong preferred orientation. The metamorphic mica is blade shaped with long axes defining the linear aspect of the foliation and intermediate axes a partial girdle about the lineation. Once initiated, foliation progressively intensified by an increase in the aspect ratio, size and alignment of grains, although highest grade samples within the chlorite zone record a decrease in aspect ratio and reduction in SPO strength despite continued increase in grain size. These trends are interpreted in terms of progressive competitive anisotropic growth of blade‐shaped grains so that the fastest growth directions and blade lengths tend to parallel the extension direction during deformation. The competitive nature of mica growth is indicated by the progressive increase in size and resultant decrease in number of metamorphic mica with increasing grade, from c. 1000 relatively small mica grains per square millimetre of thin section at lower grades, to c. 100 relatively large grains per square millimetre in higher grade samples. Reversal of SPO intensity and grain aspect ratio trends in higher grade samples may reflect a reduction in the strain rate or reduction in the deviatoric component of the stress field.     

Anisotropic Variation Law of Rock Permeability with the Burial Depth of Limestone

Permeability tensors of both macrofracture and microfracture systems were measured progressively along the depth of limestone formations at severed sites. It was found that the principal permeability values KX, Ky and Kz in these permeability tensors all decrease simultaneously and logarithmically with depth. However, the limestone aquifers are composed of an upper region where the larger permeability ellipsoid is upright or prolate and characterized by Kz>Kx and Kz>Ky, a transitional zone, and a lower zone whose smaller permeability ellipsoid is horizontal or oblate and characterized by KZ相似文献   

A field study (Massachusetts, USA) of the factors controlling the depth of groundwater flow systems in crystalline fractured-rock terrain

Groundwater movement and availability in crystalline and metamorphosed rocks is dominated by the secondary porosity generated through fracturing. The distributions of fractures and fracture zones determine permeable pathways and the productivity of these rocks. Controls on how these distributions vary with depth in the shallow subsurface (<300 m) and their resulting influence on groundwater flow is not well understood. The results of a subsurface study in the Nashoba and Avalon terranes of eastern Massachusetts (USA), which is a region experiencing expanded use of the fractured bedrock as a potable-supply aquifer, are presented. The study logged the distribution of fractures in 17 boreholes, identified flowing fractures, and hydraulically characterized the rock mass intersecting the boreholes. Of all fractures encountered, 2.5% are hydraulically active. Boreholes show decreasing fracture frequency up to 300 m depth, with hydraulically active fractures showing a similar trend; this restricts topographically driven flow. Borehole temperature profiles corroborate this, with minimal hydrologically altered flow observed in the profiles below 100 m. Results from this study suggest that active flow systems in these geologic settings are shallow and that fracture permeability outside of the influence of large-scale structures will follow a decreasing trend with depth.     

A shear‐dilation‐based model for evaluation of hydraulically stimulated naturally fractured reservoirs

The role of shear dilation as a mechanism of enhancing fluid flow permeability in naturally fractured reservoirs was mainly recognized in the context of hot dry rock (HDR) geothermal reservoir stimulation. Simplified models based on shear slippage only were developed and their applications to evaluate HDR geothermal reservoir stimulation were reported. Research attention is recently focused to adjust this stimulation mechanism for naturally fractured oil and gas reservoirs which reserve vast resources worldwide. This paper develops the overall framework and basic formulations of this stimulation model for oil and gas reservoirs. Major computational modules include: natural fracture simulation, response analysis of stimulated fractures, average permeability estimation for the stimulated reservoir and prediction of an average flow direction. Natural fractures are simulated stochastically by implementing ‘fractal dimension’ concept. Natural fracture propagation and shear displacements are formulated by following computationally efficient approximate approaches interrelating in situ stresses, natural fracture parameters and stimulation pressure developed by fluid injection inside fractures. The average permeability of the stimulated reservoir is formulated as a function of discretized gridblock permeabilities by applying cubic law of fluid flow. The average reservoir elongation, or the flow direction, is expressed as a function of reservoir aspect ratio induced by directional permeability contributions. The natural fracture simulation module is verified by comparing its results with observed microseismic clouds in actual naturally fractured reservoirs. Permeability enhancement and reservoir growth are characterized with respect to stimulation pressure, in situ stresses and natural fracture density applying the model to two example reservoirs. Copyright © 2002 John Wiley & Sons, Ltd.     

兴隆金矿区韧性剪切带特征及其含矿性研究

兴隆金矿区韧性剪切带受兴隆沟深断裂影响,其韧性剪切变质变形控制着本区金矿(化)体的空间分布。韧性剪切带由片麻岩、糜棱岩化岩、糜棱岩及超糜棱岩组成,矿物组分和化学成分发生了明显变化,促使金的活化带出,在适当温压条件下充填并沉淀于叶理构造中而形成韧性剪切带型金矿。     

Composition and structure of minerals from Paleoproterozoic sedimentary rocks in the section of the Kola Superdeep Borehole and their homologs in outcrops

Geochemical and structural features of microcline, quartz, and carbonates from the sedimentary Luchlompolo Formation (III SF) of the Pechenga volcano sedimentary complex were compared in outcrops at a depth of 20–400 m in the section of satellite Borehole IX and of 4.8 km in the Kola Superdeep Borehole (SD-3). The rocks occupying an identical position in the section of this formation are homologous and characterized by similar composition and spectroscopic properties of rock-forming minerals. A certain variance in defect structure is revealed in minerals from shallow and deep levels. The crystal structure of microcline and dolomite from deep levels of the Kola Superdeep is distinguished by a larger amount of defects. The isomorphic replacement of Ca with Mn ions in dolomite increases with depth, probably owing to more intense metamorphism and further long-term standing at elevated pressure and temperature. The concentration of Al admixture in quartz decreases with depth, and isomorphic substitution of Al for Si shifts toward compensation with the alkali element in response to the same factors. The decrease in the number of natural paramagnetic Al centers in quartz and decrease in oxygen vacancies with depth is regarded as an effect of the current PT conditions.     

Geotechnical Evaluation of Dam Foundation with Special Reference to In Situ Permeability: A Case Study

Since every dam site is unique and different from the rest, special geological and geotechnical investigation for each dam is a key factor in the success of the project. The Peygham-chay is a 55 m in height is an embankment dam with a clay core that is situated to the northwest of Iran and is under construction at present. The dam is founded on a valley filled by 12 m alluvium deposits that thicknesses increase toward left abutment up to 60 m. The dam site investigated based on surface discontinuity surveying, drilled boreholes logging data, in situ and laboratory tests. Detailed and extensive surface and subsurface explorations have been carried out at the dam site to characterize the geotechnical properties of foundation. The geology of the site consists of a series of cretaceous basic lavas intrusion composed with andesite, basalt and diabase that cropped out on the right abutments and is strongly fractured due to numerous joints and faults. Permeability of fissured rock mass is strongly depends on joint properties (degree of jointing, opening, continuity, filling material and weathering). In this paper in situ permeability of foundation was evaluated based on water pressure tests (lugeon) results and secondary permeability index (SPI). The impact of water pressure on joints characteristics and consequently flow rate are investigated and considered in the results. The overall groutability of rock mass was evaluated by comparison of RQD as an index of rock jointing degree and SPI values. The results indicated that permeability of rock mass in more than 66% of the tests is very low and not required for treatment. Meanwhile in the rest of sections the construction of preferred sealing system in necessary.     

低应力与覆水环境下单裂隙粉砂质泥岩渗流特性

为探究不同外部环境因素影响下浅层粉砂质泥岩边坡裂隙渗流特性,采用自主研发的岩体裂隙渗流试验装置,对含6种不同裂隙面粗糙度（JRC）的粉砂质泥岩裂隙试样进行渗流试验,研究了不同低围压和覆水深度下粉砂质泥岩裂隙渗流特性。结果表明：不同覆水深度及JRC下围压与粉砂质泥岩裂隙渗透系数均呈反相关,两者之间关系可用幂函数表征,且渗透系数的降低过程可分为快速降低（围压为0～30 kPa）和缓慢降低（围压为30～50 kPa）两个阶段,CT扫描结果验证了围压增大使得粉砂质泥岩裂隙开度减小是渗透系数随围压增大而减小的主要原因。随围压的增大或覆水深度的减小,不同JRC粉砂质泥岩裂隙渗透系数的离散程度逐渐减小。当围压增至最大,同时覆水深度最小时,JRC对裂隙渗透系数的影响将会被消除。不同围压下,粉砂质泥岩裂隙渗透系数与覆水深度呈正相关,且两者关系可用指数函数表征。推导出了粉砂质泥岩裂隙渗流非线性Izabsh模型,该模型能较好地反映低应力及低流速下粉砂质泥岩裂隙渗流量与压力梯度之间的非线性变化关系,但随围压的增大,该模型的相关性有一定程度的减小。     

Influence of permeability in modeling of reservoir triggered seismicity in Koyna region,western India

The Koyna region located in the west coast of India is a classic example of reservoir triggered seismicity (RTS) that started soon after the impoundment of the Koyna reservoir in 1962. Previous studies have shown that RTS can be explained in terms of stress and pore pressure changes due to poroelastic response of the rock matrix. The permeability of rock matrix is a key parameter for pore pressure diffusion which is mainly responsible for generation of stress perturbation related to seismicity. Based on the poroelastic theory, we employ 2-D finite element models to simulate the evolution of pore pressure up to 5 years after the reservoir impoundment in 1962, using a range in permeability, 10^?16–10^?14 m². Constraints on material properties of Deccan basalt and granitic rocks were taken from available studies. The results show the formation of pore pressure front and its propagation with depth and time since the reservoir impoundment as a function of permeability. While a permeability of 10^?16 m² does not produce any significant change in pore pressure, a ten-fold increase in permeability produces significant changes up to a depth of 2 km only beneath the reservoir after 5 years of impoundment. Permeability values between 10^?15 m² and 10^?14 m² are required to induce critical pore pressure changes in the range 0.1–1 MPa up to depth of 10 km, capable of triggering earthquakes in a critically stressed region. Studies on core samples of granitic basement rock down to a depth of 1522 m in the Koyna region provide evidences of fracture zones that may contribute to water channelization. Direct measurements of material properties through the ongoing deep drilling programme would help to develop more realistic models of RTS.     

Physical and mechanical properties of selected amphibolite core samples from the Kola Superdeep Borehole KSDB-3

Laboratory tests of physical and mechanical properties were performed on amphibolites representing two units from the Kola Superdeep Borehole KSDB‐3 (Russia) ? the Proterozoic Karelian Complex (depths of 3043 m, 3530 m and 4389 m) and the Archean Kola–Belomorian Complex (depths of 7951 m, 8942 m and 9904 m). Obtained grain density and bulk density values are nearly identical for all tested samples, but porosities increased slowly with depth. Marked differences in strength properties were found between amphibolite samples from the Karelian Complex and amphibolite samples from the Kola–Belomorian Complex. The uniaxial compressive strength of the samples from the lower part of the borehole is significantly smaller than in samples from the upper part; a similar trend was found for triaxial tests. The values of the deformation modulus, Young's modulus and Poisson's ratio obtained from samples of the Karelian Complex are considerably higher than those from the Kola–Belomorian Complex. Observed mechanical properties are explained by the variable grain‐size distribution and by the spatial arrangement of main rock‐forming minerals. These properties are also influenced by recovery of rocks from great depths.     

Permeability and P-wave velocity change in granitic rocks under freeze–thaw cycles

An extensive experimental investigation of microstructural changes in granites under freeze–thaw cycles using permeability and P-wave velocity measurements is described. Two types of natural granite rocks are considered and tested under dry and saturated conditions. The specimens were subjected to 200 heating–cooling cycles (??20°C/?+?20°C); each cycle had a duration of 24 h. The results indicate that the ageing process decreases the permeability and the P-wave velocity. The reduction in P-wave velocity is likely to be due to microcrack development (material damage). In dry samples, the microcracks result from the repeated differential contraction–dilatation of the mineral components. In water-saturated samples, there is an additional effect of freezing and thawing of water in the porous network. The decrease in permeability in the dry samples is due to partial closure of existing microcracks. In water-saturated samples, there was no increase in the permeability. A physically acceptable explanation is that new microcracks are not necessarily connected with those that already exist. Therefore the physicochemical process resulting from water–rock interactions also affects the permeability. This phenomenon reduces fluid flow in the material.