Analysis of a deformable fracture in permeable material

The response of deformable fractures to changes in fluid pressure controls phenomena ranging from the flow of fluids near wells to the propagation of hydraulic fractures. We developed an analysis designed to simulate fluid flows in the vicinity of asperity‐supported fractures at rest, or fully open fractures that might be propagating. Transitions between at‐rest and propagating fractures can also be simulated. This is accomplished by defining contact aperture as the aperture when asperities on a closing fracture first make contact. Locations on a fracture where the aperture is less than the contact aperture are loaded by both fluid pressure and effective stress, whereas locations where the aperture exceeds the contact aperture are loaded only by fluid pressure. Fluid pressure and effective stress on the fracture are determined as functions of time by solving equations of continuity in the fracture and matrix, and by matching the global displacements of the fracture walls to the local deformation of asperities. The resulting analysis is implemented in a numerical code that can simulate well tests or hydraulic fracturing operations. Aperture changes during hydraulic well tests can be measured in the field, and the results predicted using this analysis are similar to field observations. The hydraulic fracturing process can be simulated from the inflation of a pre‐existing crack, to the propagation of a fracture, and the closure of the fracture to rest on asperities or proppant. Two‐dimensional, multi‐phase fluid flow in the matrix is included to provide details that are obscured by simplifications of the leakoff process (Carter‐type assumptions) used in many hydraulic fracture models. Execution times are relatively short, so it is practical to implement this code with parameter estimation algorithms to facilitate interpretation of field data. Copyright © 2006 John Wiley & Sons, Ltd.     

Modelling fluid flow through a single fracture with different contacts using cellular automata

The representation of the surface of a rock fracture and a numerical method to simulate fluid flow in single fractures are the keys to understanding the hydraulic behaviour of rock fractures. In this paper, a cellular automaton (CA) approach is used to generate the single fracture structure, which is assumed to be composed of contacts and voids. We develop a CA evolution rule to produce a contact area, and randomly model a single rock fracture with different contact ratios to reflect natural fracture properties such as dead voids, islands and tortuous flow path. Then, based on the localisation theory of a CA, a numerical method to simulate fluid flow in single fractures with contacts is developed. In this method, the fracture is discretised into a system composed of cell elements. Different apertures, i.e., zero for contacts and non-zero for voids, are assigned to each cell element. Therefore, the contribution of the cell elements in a contact on a cell’s transmissivity can be ignored completely. The local transmissivity is assumed to conform to the cubic law. The fluid flow in a fracture with different contact situations is then modelled using the method established in this paper. The fluid flow path, flow velocity and fluid head distributions as well as the channel flow in the fracture are well-modelled. The flow behaviour of the fracture strongly depends on the effective fluid flow path.     

Fracture Permeability Alteration due to Chemical and Mechanical Processes: A Coupled High-Resolution Model

Reactive fluid-flow experiments in fractures subjected to normal stress suggest the potential for either increased or decreased permeability resulting from fracture-surface dissolution. We present a computational model that couples mechanical deformation and chemical alteration of fractures subjected to constant normal stress and reactive fluid flow. The model explicitly represents micro-scale roughness of the fracture surfaces and calculates elastic deformation of the rough surfaces using a semi-analytical approach that ensures the surfaces remain in static equilibrium. A depth-averaged reactive transport model calculates chemical alteration of the surfaces, which leads to alteration of the contacting fracture surfaces. The mechanical deformation and chemical alteration calculations are explicitly coupled, which is justified by the disparate timescales required for equilibration of mechanical stresses and reactive transport processes. An idealized analytical representation of dissolution from a single contacting asperity shows that under reaction-limited conditions, contacting asperities can dissolve faster than the open regions of the fracture. Computational simulations in fractures with hundreds of contacting asperities show that the transition from transport-limited conditions (low flow rates) to reaction-rate-limited conditions (high flow rates) causes a shift from monotonically increasing permeability to a more complicated process in which permeability initially decreases and then increases as contacting asperities begin to dissolve. These results are qualitatively consistent with a number of experimental observations reported in the literature and suggest the potential importance of the relative magnitude of mass transport and reaction kinetics on the evolution of fracture permeability in fractures subjected to combined normal stress and reactive fluid flow.     

Contact mechanism of a rock fracture subjected to normal loading and its impact on fast closure behavior during initial stage of fluid flow experiment

Fast closure of rock fractures has been commonly observed in the initial stage of fluid flow experiments at environmental temperatures under low or moderate normal stresses. To fully understand the mechanisms that drive this fast closure, the evolution of local stresses acting on contacting asperities on the fracture surfaces prior to fluid flow tests needs to be evaluated. In this study, we modeled numerically the asperity deformation and failure processes during initial normal loading, by adopting both elastic and elastic–plastic deformation models for the asperities on a real rock fracture with measured surface topography data, and estimated their impact on initial conditions for fluid flow test performed under laboratory conditions. Compared with the previous models that simulate the normal contact of a fracture as the approach of two rigid surfaces without deformations, our models of deformable asperities yielded smaller contact areas and higher stresses on contacting asperities at a given normal stress or normal displacement. The results show that the calculated local stresses were concentrated on the contacts of a few major asperities, resulting in crushing of asperity tips. With these higher contact stresses, however, the predicted closure rates by pressure solution are still several orders of magnitude lower than that of the experimental measurements at the initial stage of fluid flow test. This indicates that single pressure solution may not likely to be the principal compaction mechanism for this fast closure, and that the damages on contacting asperities that occur during the initial normal loading stage may play an important role. Copyright © 2015 John Wiley & Sons, Ltd.     

Simulation of fluid percolation in a rough-walled rock fracture

A comparison of experimental and numerical results is presented addressing two-phase immiscible displacement in rough fractures. Quasi-static, capillary displacement in a rough fracture is modeled using the modified invasion percolation approach proposed by Glass et al (1998), and the results are compared against experimental observations obtained from two-phase flow through a rock fracture using x-ray computed-tomography scanning. The model is based on an algorithm seeking the least resistant pathway for the advancement of the invading fluid using the Young-Laplace equation and accounting for local in-plane curvature of the advancing fluid front. The saturation distribution map generated by the model yields good agreement with the experimental phase distribution and presents more realistic phase structures than those obtained from the conventional invasion percolation approach. The improvement in the results obtained with the modified invasion percolation approach is attributed to the contribution of the in-plane curvature term, which captures the effect of regionalized apertures, rather than single-point apertures, on the shape of the invading front. Glass RJ, Nicholl MJ, Yarrington L (1998) A modified invasion percolation model for low-capillary number immiscible displacements in horizontal rough-walled fractures: influence of local in-plane curvature. Water Resour Res 34:3215–3234     

Numerical simulation of two-phase flow in conceptualized fractures

Two-phase flow in fractured rock is an important phenomenon related to a range of practical problems, including non-aqueous phase liquid contamination of groundwater. Although fractured rocks consist of fracture networks, the study of two-phase flow in a single fracture is a pre-requisite. This paper presents a conceptual and numerical model of two-phase flow in a variable fracture. The void space of the fracture is conceptualized as a system of independent channels with position-dependent apertures. Fundamental equations, governing two-phase displacement in each channel, are derived to represent the interface positions and fractional flows in the fracture. For lognormal aperture distributions, simple approximations to fractional flows are obtained in analytical form by assuming void occupancy based on a local capillary allowability criterion. The model is verified by analytical solutions including two-phase flow in a parallel-plate fracture, and used to study the impacts of aperture variation, mobility ratio and fracture orientation on properties of two-phase flow. Illustrative examples indicate that aperture variation may control the distribution of wetting and non-wetting fluids within the fracture plane and hence the ability of the fracture to transmit these fluids. The presence of wetting fluid does little to hinder non-wetting fluid flow in fractures with large aperture variations, whereas a small volume of non-wetting fluid present in the fracture can significantly reduce wetting fluid flow. Large mobility ratios and high fracture slope angles facilitates the migration of non-wetting fluid through fractures.     

考虑非达西效应的酸蚀裂缝流场数值模拟

酸压改造后的酸蚀裂缝是储层流体的主要流动通道,研究流体在酸蚀裂缝中的流动规律是十分必要的。本文介绍了目前裂隙流动研究中的4种常用控制方程及各自的适用条件,确定了Navier-Stokes方程为酸蚀裂缝流场数值模拟的控制方程。为了分析酸蚀裂缝流场分布规律以及非达西效应对裂缝导流能力的影响,利用逆向工程技术对两种不同刻蚀形态的酸蚀岩样壁面进行了实体重构,并采用有限元数值模拟方法进行了不同流量下的酸蚀裂缝流动试验。结果表明：酸刻蚀形态和裂缝接触关系对流场分布的影响较大。隙宽分布平滑的均匀刻蚀裂缝,其流态稳定、曲折度低,但由于裂缝开度较窄且比表面较高,边界层效应明显,流体流动时产生的黏滞阻力较高;沟槽裂缝曲面粗糙,隙宽分布复杂,流体流动时流态不稳定,曲折度高,在高流量下会产生明显的涡流流动,增大惯性阻力。裂缝通道存在缩颈现象会引起流体绕流和多次加减速,产生额外的压力损耗。随着模拟流量的增大,压降与流量会逐渐偏离线性关系而呈现出非达西流动现象。酸蚀裂缝壁面越粗糙,产生非达西效应时的临界流量和临界雷诺数就越小,在相同流量下非达西效应就越强烈,导流能力下降速率更快。     

Fault zone characteristics,fracture systems and permeability implications of Middle Triassic Muschelkalk in Southwest Germany

Fault zone structure and lithology affect permeability of Triassic Muschelkalk limestone-marl-alternations in Southwest Germany, a region characterized by a complex tectonic history. Field studies of eight fault zones provide insights into fracture system parameters (orientation, density, aperture, connectivity, vertical extension) within fault zone units (fault core, damage zone). Results show decreasing fracture lengths with distances to the fault cores in well-developed damage zones. Fracture connectivity at fracture tips is enhanced in proximity to the slip surfaces, particularly caused by shorter fractures. Different mechanical properties of limestone and marl layers obviously affect fracture propagation and thus fracture system connectivity and permeability. Fracture apertures are largest parallel and subparallel to fault zones and prominent regional structures (e.g., Upper Rhine Graben) leading to enhanced fracture-induced permeabilities. Mineralized fractures and mineralizations in fault cores indicate past fluid flow. Permeability is increased by the development of hydraulically active pathways across several beds (non-stratabound fractures) to a higher degree than by the formation of fractures interconnected at fracture tips. We conclude that there is an increase of interconnected fractures and fracture densities in proximity to the fault cores. This is particularly clear in more homogenous rocks. The results help to better understand permeability in Muschelkalk rocks.     

Hydraulic Characteristics of Rough Fractures in Linear Flow under Normal and Shear Load

Summary A hydro-mechanical testing system, which is capable of measuring both the flow rates and the normal and shear displacement of a rock fracture, was built to investigate the hydraulic behaviour of rough tension fractures. Laboratory hydraulic tests in linear flow were conducted on rough rock fractures, artificially created using a splitter under various normal and shear loading. Prior to the tests, aperture distributions were determined by measuring the topography of upper and lower fracture surfaces using a laser profilometer. Experimental variograms of the initial aperture distributions were classified into four groups of geostatistical model, though the overall experimental variograms could be well fitted to the exponential model. The permeability of the rough rock fractures decayed exponentially with respect to the normal stress increase up to 5 MPa. Hydraulic behaviours during monotonic shear loading were significantly affected by the dilation occurring until the shear stress reached the peak strength. With the further dilation, the permeability of the rough fracture specimens increased more. However, beyond shear displacement of about 7 to 8 mm, permeability gradually reached a maximum threshold value. The combined effects of both asperity degradation and gouge production, which prohibited the subsequent enlargement of mean fracture aperture, mainly caused this phenomenon. Permeability changes during cyclic shear loading showed somewhat irregular variations, especially after the first shear loading cycle, due to the complex interaction from asperity degradations and production of gouge materials. The relation between hydraulic and mechanical apertures was analyzed to investigate the valid range of mechanical apertures to be applied to the cubic law. Received June 12, 2001; accepted February 26, 2002 Published online September 2, 2002     

剪胀对复杂裂隙岩体渗透性影响的离散元分析

裂隙岩体流固耦合问题是目前国内外研究热点之一,采用离散元软件UDEC对裂隙岩体发生节理剪胀的渗透性变化规律进行了模拟分析。基于现场调查的裂隙信息统计生成裂隙网络岩体模型。 通过固定垂直应力、不断增加应力比RS(RS=水平应力/垂直应力)使岩体出现剪胀,采用库伦滑移节理模式对岩体在剪胀过程中的渗透性变化情况进行模拟。结果发现：当应力比较小(RS3.1)时,节理水力隙宽、流速、渗透系数等参数都随着应力比的增加表现出明显的降低; 而当岩体出现剪胀现象之后(应力比大于3.1),发生剪切滑移和剪胀现象的节理控制着裂隙岩体的总体渗流行为,与不考虑节理剪胀的计算结果相比,岩体渗透能力出现了显著增长。这一结果表明,剪胀对裂隙岩体渗透性的影响是显著而不可忽视的。     

Numerical modelling of fluid flow tests in a rock fracture with a special algorithm for contact areas

The fluid flow in rock fractures during shear processes has been an important issue in rock mechanics and is investigated in this paper using finite element method (FEM), considering evolutions of aperture and transmissivity with shear displacement histories under different normal stress and normal stiffness conditions as measured during laboratory coupled shear-flow tests. The distributions of fracture aperture and its evolution during shearing were calculated from the initial aperture, based on the laser-scanned sample surface roughness results, and shear dilations measured in the laboratory tests. Three normal loading conditions were adopted in the tests: simple normal stress and mixed normal stress and normal stiffness to reflect more realistic in situ conditions. A special algorithm for treatment of the contact areas as zero-aperture elements was used to produce more accurate flow field simulations, which is important for continued simulations of particle transport but often not properly treated in literature. The simulation results agree well with the measured hydraulic apertures and flow rate data obtained from the laboratory tests, showing that complex histories of fracture aperture and tortuous flow fields with changing normal loading conditions and increasing shear displacements. With the new algorithm for contact areas, the tortuous flow fields and channeling effects under normal stress/stiffness conditions during shearing were more realistically captured, which is not possible if traditional techniques by assuming very small aperture values for the contact areas were used. These findings have an important impact on the interpretation of the results of coupled hydro-mechanical experiments of rock fractures, and on more realistic simulations of particle transport processes in fractured rocks.     

Investigating the role of hydromechanical coupling on flow and transport in shallow fractured-rock aquifers

Fractured-rock aquifers display spatially and temporally variable hydraulic conductivity generally attributed to variable fracture intensity and connectivity. Empirical evidence suggests fracture aperture and hydraulic conductivity are sensitive to in situ stress. This study investigates the sensitivity of fractured-rock hydraulic conductivity, groundwater flow paths, and advection-dominated transport to variable shear and normal fracture stiffness magnitudes for a range of deviatoric stress states. Fracture aperture and hydraulic conductivity are solved for analytically using empirical hydromechanical coupling equations; groundwater flow paths and ages are then solved for numerically using groundwater flow and advection-dispersion equations in a traditional Toth basin. Results suggest hydraulic conductivity alteration is dominated by fracture normal closure, resulting in decreasing hydraulic conductivity and increasing groundwater age with depth, and decreased depth of long flow paths with decreasing normal stiffness. Shear dilation has minimal effect on hydraulic conductivity alteration for stress states investigated here. Results are interpreted to suggest that fracture normal stiffness influences hydraulic conductivity of hydraulically active fractures and, thus, affects flow and transport in shallow (<1 km) fractured-rock aquifers. It is suggested that observed depth-dependent hydraulic conductivity trends in fractured-rock aquifers throughout the world may be partly a manifestation of hydromechanical phenomena.     

Correlations Developed for Estimation of Hydraulic Parameters of Rough Fractures Through the Simulation of JRC Flow Channels

The hydro-mechanical response of fractured rock masses is complex, due partly to the presence of fractures at different scales. Surface morphology has a significant influence on fluid flow behaviour of a fracture. Different empirical correlations and statistical models have been proposed to estimate the equivalent hydraulic aperture and determine the pressure drop along a fracture. However, the existing models suffer from not being adequately generalised to be applicable to a wide range of real fracture surfaces. To incorporate the effect of profile roughness in the hydro-mechanical behaviour of fractured rock masses, the joint roughness coefficient (JRC) is the most widely used empirical approach. However, the average JRC of two fracture walls in fluid flow analysis, as is a common practice, appears to be inappropriate. It will be shown how different combinations of pairs of JRCs could lead to a similar JRC value. Also, changing the position of the top and bottom walls of a fracture can significantly change the hydraulic response of the fracture while the average JRC is identical in both cases. In this paper, correlations are developed which are based on the simulation of JRCs using estimated fluid flow parameters of 2D fractures can be estimated. In order to widen the application range of the correlations, JRC flow channels were generated: these are 2D channels with their top and bottom walls being made from two of the JRC profiles. To estimate the JRC of linear profiles, a correlation developed between JRC and a newly developed Riemannian roughness parameter, D _R1, is proposed. Considering ten JRC profiles, a total of 100 JRC flow channels were generated. In order to only investigate the effect of surface roughness on fluid flow, the minimum closure between the top and bottom walls of JRC flow channels were considered to be constant. Three cases with minimum closures of 0.01, 0.05 and 0.10 cm were considered in this study. All JRC flow channels were subjected to fluid analysis using FLUENT software. Based on these results, correlations were developed between the geometrical and hydraulic properties of flow channels. Analysis of several real fractures demonstrated the applicability of these correlations.     

An experimental study of single and two-phase fluid flow through fractured granite specimens

Triaxial tests on the two-phase flow of air and water through fractured granite specimens were performed to discover whether the two-phase fluid flow within rock fractures was laminar or turbulent. The two-phase flow characterization was carried out based on the macroscopic two-phase steady state flow model and the homogeneous steady state flow model. Rock specimens with a single natural fracture (joint roughness coefficient, JRC < 10) were tested using two-phase, high pressure triaxial rig. Experimental results show that the estimated Reynolds numbers for various inlet fluid pressures are well below 1000. The findings of this study reveal that both single and two-phase flow through rock fractures (JRC < 10) can be characterized as laminar flows at moderate inlet fluid pressures. However, for single-phase air flow, an increase in inlet air pressures may result in the formation of turbulent flow.     

Fracture closure modes during flowback from hydraulic fractures

Flowback analysis recently has been considered as a potential tool for assessing induced fractures through corresponding pressure analyses and rate transient analysis. In this paper, we study fracture closure mechanisms during the flowback of fracturing fluid after hydraulic fracturing treatments. Although it is known that flowback can be significantly affected by the geometry of the fractures and closure stress, there has not been any effort to understand the problem from the geomechanical point of view; rather, available studies assume that a fracture closes uniformly with constant fracture compressibility. The coupled geomechanics and fluid flow model presented in this paper help to elucidate how geostresses may affect fracturing fluid recovery under different conditions. We perform a scaling analysis to formulate the occurrence of different fracture closure modes and then use numerical analyses to verify our scaling parameters. The factors governing the flowback process include the mechanical and petrophysical properties of the rock as well as preexisting discontinuities such as natural fractures. Three different closure modes for fracture networks are described and numerically verified. The occurrence of each mode can dramatically affect fracturing fluid recovery. The role of fluid leakoff into the formation, fractures conductivity, and drawdown strategy are examined for each fracture closure scenario.     

岩体裂隙网络等效渗透系数方向性的数值计算

采用渗流力学理论并结合Monte Carlo方法描述岩体裂隙的随机分布,研究渗流模型的尺寸效应并确定表征单元体积（REV）,得到了3种开口度分布形式的等效渗透系数椭圆曲线,建立了等效渗透系数方向性的判别标准。离散裂隙网络（DFN）模型假定流体只在岩体裂隙内部流动,而不通过岩体本身渗流。基于二维离散元程序UDEC并进行二次开发,建立DFN模型,通过改变流体的流动方向,得到不同流动方向下岩体裂隙网络的等效渗透系数,并分析不同的开口度分布形式对岩体裂隙网络等效渗透系数方向性的影响。计算结果表明,表征单元体积存在的条件是等效渗透系数保持稳定且渗透椭圆比较光滑。等效渗透系数的方向性受开口度分布形式的影响很大：当开口度-长度关联分布时,等效渗透系数各向异性;当开口度对数正态分布时,等效渗透系数各向同性;当开口度恒定分布时,等效渗透系数的特性介于二者之间。变化系数（CV）是否大于5%是判定岩体裂隙网络渗透系数是否具有方向性的判别标准。     

Fracture flow modelling based on satellite images of the Wajid Sandstone, Saudi Arabia

Large-scale geological features have been identified by satellite imagery and global positioning system data in the Wajid Sandstone in Saudi Arabia. The main objective is to evaluate the importance of fractures for the overall flow behaviour in this fractured rock aquifer and to estimate in-situ hydraulic apertures. Data on fractures and lineaments were available for three outcrops. By applying a “cut-out” routine on the fracture endpoint data of these fracture trace windows, three deterministic discrete fracture networks (DFN), with an area of 100 m?×?100 m, could be generated. These were used to simulate the fracture flow and to determine the hydraulic conductivity tensors. Using additional data on hydraulic pumping tests and matrix conductivities, in-situ hydraulic apertures could be determined. Average in-situ hydraulic apertures range from 1,300 to 1,700 µm. Observations from the field support these results. In addition, a hydraulic conductivity ratio between the matrix and fracture system was used to identify the contribution of the DFN to the overall fluid transport. A ratio of 10.4 was determined, which indicates that the effective flow behaviour in the Wajid Sandstone aquifer is not entirely dominated by the fracture system, though evidently strongly controlled by it.     

Natural fracture profiles,fractal dimension and joint roughness coefficients

Summary Thirteen natural rock profiles (Barton and Choubey, 1977) are analyzed for their fractal properties. Most of the profiles were found to approximate fractal curves but some also showed features of specific wavelengths and amplitudes superimposed on fractal characteristics. The profiles showed fractal dimensions from 1.1 to 1.5 covering a range of selfsimilar and self-affine curves. The analysis results suggest a negative correlation between fractal dimension,D, and amplitude,A. Joint roughness coefficients (JRC) show a positive correlation with amplitude,A, and a negative correlation with fractal dimension,D. A numerical model of fracture closure is used to investigate the effects of different profile characteristics (D, A and sample size) on the nature of dilation and contact area, using the natural profiles and synthetic fractional Brownian motion profiles. Smooth profiles (low JRC, highD, lowA) display many small contact regions whereas rough fractures (high JRC, lowD, highA) display few large contact areas. The agreement with published experimental data supports the suggested correlations between JRC and the fractal parameters,A andD. It is suggested that observed scale effects in JRC and joint dilation can be explained by small differential strain discontinuities across fractures, which originate at the time of fracture formation.     

Evolution of fluid pressure and fracture propagation during contact metamorphism

Abstract Rock fracture enhances permeability and provides pathways through which fluids migrate. During contact metamorphism, fluids contained in isolated pores and fractures expand in response to temperature increases caused by the dissipation of heat from magmas. Heat transport calculations and thermomechanical properties of water-rich fluids demonstrate (1) that thermal energy is a viable mechanism to produce and maintain pore fluid pressure (P_f) in a contact metamorphic aureole; (2) that the magnitude of P_f generated is sufficient to propagate fractures during the prograde thermal history (cause hydrofracture) and enhance permeability; and (3) that P_f-driven fracture propagation is episodic with time-scales ranging from years to thousands of years. Because P_f dissipation is orders of magnitude faster than P, _f buildup, P_f oscillations and cyclical behaviour are generated as thermal heating continues. The P_f cycle amplitude depends on the initial fracture length, geometry and the rock's resistance to failure whereas the frequency of fracture depends on the rate of heating. Consequently, oscillation frequency also varies spatially with distance from the heat source. Time series of fluid pressures caused by this process suggest that cyclical fracture events are restricted to an early time period of the prograde thermal event near the intrusive contact. In the far field, however, individual fracture events have a lower frequency but continue to occur over a longer time interval. Numerous fracture cycles are possible within a single thermal event. This provides a provisional explanation for multiple generations of veins observed in outcrop. P _f cycling and oscillations may explain several petrological features. If pore fluids are trapped at various positions along a pressure cycle, the large amplitude of P_f variations for small fractures may account for different pressures recorded by fluid inclusions analysed from a single sample. P_f oscillations, during a single thermal episode, also drive chemical reactions which can produce complex mineral textures and assemblages for discontinuous reactions and/or zoning patterns for continuous reactions. These can mimic polymetamorphic or disequilibrium features. Temporal aspects of fracture propagation and permeability enhancement also constrain the likely timing of fluid flow and fluid-mineral interactions. These data suggest that fluid flow and fluid-mineral reactions are likely to be restricted to an early period in the prograde thermal history, characterized by high P_f coincident with relatively high temperatures, fracture propagation and consequent increases in permeability. This early prograde hydration event is followed by diffusional peak metamorphic reactions. This relationship is evident in the complex mineralogical textures common in some metamorphosed rocks.     

Cross-borehole flow analysis to characterize fracture connections in the Melechov Granite, Bohemian-Moravian Highland, Czech Republic

Application of the cross-borehole flow method, in which short pumping cycles in one borehole are used to induce time-transient flow in another borehole, demonstrated that a simple hydraulic model can characterize the fracture connections in the bedrock mass between the two boreholes. The analysis determines the properties of fracture connections rather than those of individual fractures intersecting a single borehole; the model contains a limited number of adjustable parameters so that any correlation between measured and simulated flow test data is significant. The test was conducted in two 200-m deep boreholes spaced 21?m apart in the Melechov Granite in the Bohemian-Moravian Highland, Czech Republic. Transient flow was measured at depth stations between the identified transmissive fractures in one of the boreholes during short-term pumping and recovery periods in the other borehole. Simulated flows, based on simple model geometries, closely matched the measured flows. The relative transmissivity and storage of the inferred fracture connections were corroborated by tracer testing. The results demonstrate that it is possible to assess the properties of a fracture flow network despite being restricted to making measurements in boreholes in which a local population of discrete fractures regulates the hydraulic communication with the larger-scale aquifer system.