共查询到20条相似文献,搜索用时 46 毫秒
1.
Fluid conductivity and elastic properties in fractures depend on the aperture geometry – in particular, the roughness of fracture surfaces. In this study, we have characterized the surface roughness with a log-normal distribution and investigated the transport and flow behaviour of the fractures with varying roughness characteristics. Numerical flow and transport simulations have been performed on a single two-dimensional fracture surface, whose aperture geometry changes with different variances and correlation lengths in each realization. We have found that conventional measurement of hydraulic conductivity alone is insufficient to determine these two parameters. Transient transport measurements, such as the particle breakthrough time, provide additional constraints to the aperture distribution. Nonetheless, a unique solution to the fracture aperture distribution is still under-determined with both hydraulic conductivity and transport measurements. From numerical simulations at different compression states, we have found that the flow and transport measurements exhibit different rates of changes with respect to changes in compression. Therefore, the fracture aperture distribution could be further constrained by considering the flow and transport properties under various compression states. 相似文献
2.
3.
Numerical simulations of variable-density flow and solute transport have been conducted to investigate dense plume migration for various configurations of 2D fracture networks. For orthogonal fractures, simulations demonstrate that dispersive mixing in fractures with small aperture does not stabilize vertical plume migration in fractures with large aperture. Simulations in non-orthogonal 2D fracture networks indicate that convection cells form and that they overlap both the porous matrix and fractures. Thus, transport rates in convection cells depend on matrix and fracture flow properties. A series of simulations in statistically equivalent networks of fractures with irregular orientation show that the migration of a dense plume is highly sensitive to the geometry of the network. If fractures in a random network are connected equidistantly to the solute source, few equidistantly distributed fractures favor density-driven transport. On the other hand, numerous fractures have a stabilizing effect, especially if diffusive transport rates are high. A sensitivity analysis for a network with few equidistantly distributed fractures shows that low fracture aperture, low matrix permeability and high matrix porosity impede density-driven transport because these parameters reduce groundwater flow velocities in both the matrix and the fractures. Enhanced molecular diffusion slows down density-driven transport because it favors solute diffusion from the fractures into the low-permeability porous matrix where groundwater velocities are smaller. For the configurations tested, variable-density flow and solute transport are most sensitive to the permeability and porosity of the matrix, which are properties that can be determined more accurately than the geometry and hydraulic properties of the fracture network, which have a smaller impact on density-driven transport. 相似文献
4.
Fractures have been recently identified in potential host rock for high level nuclear waste disposal, like indurated argilite formations. These fractures appear as potential rapid pathways for radionuclides transport and hydrodynamic properties of the transport inside these systems must thus be characterized. Miscible non-sorbing and sorbing tracers displacements were performed on a 2-D model derived from a real fracture geometry observed in the Tournemire argilite formation with a lattice-gas cellular automaton (LGA). LGA was shown to easily handle the complex geometry of such a fracture. The numerical breakthrough curves obtained were inverted with the 1-D CDE and MIM transport models. Two main conclusions were drawn: (i) at the length scale of the study, the non-sorbing tracer transport in our fracture geometry was more accurately interpreted in terms of the MIM model rather than in terms of the classical CDE model; (ii) in order to correctly model the sorbing tracers migration, the hydrodynamic dispersion coefficient value was found to increase with the increase of the retardation factor. A semi-empirical relation based on the Taylor–Aris theory was then used to describe this dependency. 相似文献
5.
This study formulates and analyzes continuous time random walk (CTRW) models in radial flow geometries for the quantification of non-local solute transport induced by heterogeneous flow distributions and by mobile–immobile mass transfer processes. To this end we derive a general CTRW framework in radial coordinates starting from the random walk equations for radial particle positions and times. The particle density, or solute concentration is governed by a non-local radial advection–dispersion equation (ADE). Unlike in CTRWs for uniform flow scenarios, particle transition times here depend on the radial particle position, which renders the CTRW non-stationary. As a consequence, the memory kernel characterizing the non-local ADE, is radially dependent. Based on this general formulation, we derive radial CTRW implementations that (i) emulate non-local radial transport due to heterogeneous advection, (ii) model multirate mass transfer (MRMT) between mobile and immobile continua, and (iii) quantify both heterogeneous advection in a mobile region and mass transfer between mobile and immobile regions. The expected solute breakthrough behavior is studied using numerical random walk particle tracking simulations. This behavior is analyzed by explicit analytical expressions for the asymptotic solute breakthrough curves. We observe clear power-law tails of the solute breakthrough for broad (power-law) distributions of particle transit times (heterogeneous advection) and particle trapping times (MRMT model). The combined model displays two distinct time regimes. An intermediate regime, in which the solute breakthrough is dominated by the particle transit times in the mobile zones, and a late time regime that is governed by the distribution of particle trapping times in immobile zones. These radial CTRW formulations allow for the identification of heterogeneous advection and mobile-immobile processes as drivers of anomalous transport, under conditions relevant for field tracer tests. 相似文献
6.
The continuous time random walk (CTRW) has both an elegant mathematical theory and a successful record at modeling solute transport in the subsurface. However, there are some interpretation ambiguities relating to the relationship between the discrete CTRW transition distributions and the underlying continuous movement of solute that have not been addressed in existing literature. These include the exact definition of “transition”, and the extent to which transition probability distributions are unique/quantifiable from data. Here, we present some theoretical results which address these uncertainties in systems with an advective bias. Simultaneously, we present an alternative, reduced parameter CTRW formulation for general advective transport in heterogeneous porous media, which models early- and late-time transport by use of random transition times between sparse, imaginary planes normal to flow. We show that even in the context of this reduced-parameter formulation there is nonuniqueness in the definitions of both transition lengths and waiting time distributions, and that neither may be uniquely determined from experimental data. For practical use of this formulation, we suggest Pareto transition time distributions, leading to a two-degree-of-freedom modeling approach. We then demonstrate the power of this approach in fitting two sets of existing experimental data. While the primary focus is the presentation of new results, the discussion is designed to be pedagogical and to provide a good entry point into practical modeling of solute transport with the CTRW. 相似文献
7.
A systematic numerical method has been presented to investigate the constitutive relationships between two-phase flow properties of horizontal fractures and aperture distributions. Based on fractal geometry, single rough-walled fractures are generated numerically by modified successive random addition (SRA) method and then aperture distributions with truncated Gaussian distribution are formed by shear displacement between lower and upper surfaces. (The truncated Gaussian distribution is used to describe aperture evolution under different normal stresses.) According to the assumption of two-dimensional porous media and local parallel plate model, invasion percolation approach is employed to model the two-phase flow displacement (imbibition) in generated horizontal fractures, in which capillary forces are dominant over viscous and gravity forces. For truncated Gaussian distributions, constitutive relationships from numerical simulation are compared to closed-form relationships and a good agreement is obtained. The simulation results indicate strong phase interference with the sum of two phase relative permeability values being less than one in the intermediate saturations. It is found that fracture properties related to residual saturations depend on spatial correlation of aperture distributions. Based on the simulation results, we proposed an empirical relationship between the fracture residual-saturation-rated parameters and the corresponding aperture distributions. 相似文献
8.
9.
Understanding the Geometry of Connected Fracture Flow with Multiperiod Oscillatory Hydraulic Tests
下载免费PDF全文
![点击此处可从《Ground water》网站下载免费的PDF全文](/ch/ext_images/free.gif)
An understanding of the spatial and hydraulic properties of fast preferential flow pathways in the subsurface is necessary in applications ranging from contaminant fate and transport modeling to design of energy extraction systems. One method for the characterization of fracture properties over interwellbore scales is Multiperiod Oscillatory Hydraulic (MOH) testing, in which the aquifer response to oscillatory pressure stimulations is observed. MOH tests were conducted on isolated intervals of wells in siliciclastic and carbonate aquifers in southern Wisconsin. The goal was to characterize the spatial properties of discrete fractures over interwellbore scales. MOH tests were conducted on two discrete fractured intervals intersecting two boreholes at one field site, and a nest of three piezometers at another field site. Fracture diffusivity estimates were obtained using analytical solutions that relate diffusivity to observed phase lag and amplitude decay. In addition, MOH tests were used to investigate the spatial extent of flow using different conceptual models of fracture geometry. Results indicated that fracture geometry at both field sites can be approximated by permeable two‐dimensional fracture planes, oriented near‐horizontally at one site, and near‐vertically at the other. The technique used on MOH field data to characterize fracture geometry shows promise in revealing fracture network characteristics important to groundwater flow and transport. 相似文献
10.
Modeling fracture porosity development using simple growth laws 总被引:2,自引:0,他引:2
A model of porosity development has been developed to investigate general relationships between simple fracture aperture growth laws and fracture porosity in evolved fracture arrays in aquifers. The growth of fracture apertures in two-dimensional orthogonal arrays with initially spatially uncorrelated lognormal aperture distributions has been studied, where aperture growth rate is proportional to an exponent of the flow rate through each fracture. The evolved arrays show geometrical phase changes as a function of the aperture growth rate exponent, e, and the standard deviation of the initial aperture distribution, sigma(z). Low values of e and sigma(z) lead to bimodal aperture distributions, where apertures parallel to flow are preferentially enlarged. At moderate values of e and sigma(z), there is a transition to a regime of more complex geometries consisting of networks of channel-like structures of preferentially enlarged apertures. At larger values of e, array-spanning channel-like paths of preferentially enlarged apertures develop, where the tortuosity of the channel-like paths is a linear function of sigma(z). Following an initial growth phase, during which dynamically stable aperture configurations develop, arrays undergo simple amplification. The geometry of the evolved aperture fields is diverse and they can be highly complex; consequently, parameterization and prediction of their evolution in terms of the initial aperture distributions and growth rate laws is not trivial. 相似文献
11.
Numerical techniques for subsurface flow and transport modeling are often limited by computational limitations including fine mesh and small time steps to control artificial dispersion. Particle-tracking simulation offers a robust alternative for modeling solute transport in subsurface formations. However, the modeling scale usually differs substantially from the rock measurement scale, and the scale-up of measurements have to be made accounting for the pattern of spatial heterogeneity exhibited at different scales. Therefore, it is important to construct accurate coarse-scale simulations that are capable of capturing the uncertainties in reservoir and transport attributes due to scale-up. A statistical scale-up procedure developed in our previous work is extended by considering the effects of unresolved (residual) heterogeneity below the resolution of the finest modeling scale in 3D. First, a scale-up procedure based on the concept of volume variance is employed to construct realizations of permeability and porosity at the (coarse) transport modeling scale, at which flow or transport simulation is performed. Next, to compute various effective transport parameters, a series of realizations exhibiting detailed heterogeneities at the fine scale, whose domain size is the same as the transport modeling scale, are generated. These realizations are subjected to a hybrid particle-tracking simulation. Probabilistic transition time is considered, borrowing the idea from the continuous time random walk (CTRW) technique to account for any sub-scale heterogeneity at the fine scale level. The approach is validated against analytical solutions and general CTRW formulation. Finally, coarse-scale transport variables (i.e., dispersivities and parameterization of transition time distribution) are calibrated by minimizing the mismatch in effluent history with the equivalent averaged models. Construction of conditional probability distributions of effective parameters is facilitated by integrating the results over the entire suite of realizations. The proposed method is flexible, as it does not invoke any explicit assumption regarding the multivariate distribution of the heterogeneity. In contrast to other hierarchical CTRW formulation for modeling multi-scale heterogeneities, the proposed approach does not impose any length scale requirement regarding sub-grid heterogeneities. In fact, it aims to capture the uncertainty in effective reservoir and transport properties due to the presence of heterogeneity at the intermediate scale, which is larger than the finest resolution of heterogeneity but smaller than the representative elementary volume, but it is often comparable to the transport modeling scale. 相似文献
12.
— A new method for characterising the detailed fracture geometry in sandstone cores is presented. This method is based on the impregnation of samples with coloured resin, without significant disturbance of the fractures. The fractures are made clearly visible by the resin, thus allowing the fracture geometry to be examined digitally. In order to model the bulk hydraulic conductivity, the samples are sectioned serially perpendicular to the flow direction. The hydraulic conductivity of individual sections is estimated by summing the contribution of the matrix and each fracture from the digital data. Finally, the hydraulic conductivity of the bulk sample is estimated by a harmonic average in series along the flow path. Results of this geometrical method are compared with actual physical conductivity values measured from fluid experiments carried out prior to sectioning. The predicted conductivity from the fracture geometry parameters (e.g., fracture aperture, fracture width, fracture length and fracture relative roughness all measured using an optical method) is in good agreement with the independent physical measurements, thereby validating the approach. 相似文献
13.
The fractal geometry of flow paths in natural fractures in rock and the approach to percolation 总被引:2,自引:0,他引:2
The distributions of contact areas in single, natural fractures in quartz monzonite (Stripa granite) are found to have fractal dimensions which decrease fromD=2.00 to values nearD=1.96 as stress normal to the fractures is increased from 3 MPa up to 85 MPa. The effect of stress on fluid flow is studied in the same samples. Fluid transport through a fracture depends on two properties of the fracture void space geometry. the void aperture; and the tortuosity of the flow paths, determined through the distribution of contact area. Each of these quantities change under stress and contribute to changes observed in the flow rate. A general flow law is presented which separates these different effects. The effects of tortuosity on flow are largely governed by the proximity of the flow path distribution to a percolation threshold. A fractal model of correlated continuum percolation is presented which quantitatively reproduces the flow path geometries. The fractal dimension in this model is fit to the measured fractal dimensions of the flow systems to determine how far the flow systems are above the percolation threshold. 相似文献
14.
Numerical Study on Behavior of Groundwater Flow in Single Rough Fractures Under Mechanical Effect
下载免费PDF全文
![点击此处可从《Ground water》网站下载免费的PDF全文](/ch/ext_images/free.gif)
Woochang Jeong 《Ground water》2017,55(1):38-50
This paper numerically investigates the characteristics of groundwater flow in spatially correlated variable aperture fractures under the mechanical effect. Spatially correlated aperture distributions are generated using the geostatistical method (i.e., turning bands algorithm in this study). To represent a nonlinear relationship between the effective normal stress and the fracture aperture, a simple mechanical formula is combined with a local flow model. Numerical results indicate that the groundwater flow is significantly affected by the geometry of aperture distribution, varying with the applied effective normal stress as well as the spatial correlation length of aperture distribution. Moreover, using the flow results simulated in this study, two empirical formulae are proposed: (1) the first one (modified Louis formula) is to represent the relationship between the effective normal stress and the effective permeability of fracture and (2) the second one is to represent the relationship between relative roughness and effective permeability. 相似文献
15.
Sung‐Hoon Ji Yong‐Kwon Koh Kristopher L. Kuhlman Moo Yul Lee Jong Won Choi 《Ground water》2013,51(2):298-304
In a series of field experiments, we evaluate the influence of a small water pressure change on fracture aperture during a hydraulic test. An experimental borehole is instrumented at the Korea Atomic Energy Research Institute (KAERI) Underground Research Tunnel (KURT). The target fracture for testing was found from the analyses of borehole logging and hydraulic tests. A double packer system was developed and installed in the test borehole to directly observe the aperture change due to water pressure change. Using this packer system, both aperture and flow rate are directly observed under various water pressures. Results indicate a slight change in fracture hydraulic head leads to an observable change in aperture. This suggests that aperture change should be considered when analyzing hydraulic test data from a sparsely fractured rock aquifer. 相似文献
16.
The objective of this work is to develop a new numerical approach for the three-dimensional modelling of flow and transient solute transport in fractured porous media which would provide an accurate and efficient treatment of 3D complex geometries and inhomogeneities. For this reason, and in order to eliminate as much as possible the number of degrees of freedom, the fracture network, fractures and their intersections, are solved with a coupled 2D–1D model while the porous matrix is solved independently with a 3D model. The interaction between both models is accounted for by a coupling iterative technique. In this way it is possible to improve efficiency and reduce CPU usage by avoiding 3D mesh refinements of the fractures. The approach is based on the discrete-fracture model in which the exact geometry and location of each fracture in the network must be provided as an input. The formulation is based on a multidimensional coupling of the boundary element method-multidomain (BEM-MD) scheme for the flow and boundary element dual reciprocity method-multidomain (BE-DRM-MD) scheme for the transport. Accurate results and high efficiency have been obtained and are reported in this paper. 相似文献
17.
A Simple Model for Deviations from the Cubic Law for a Fracture Undergoing Dilation or Closure 总被引:1,自引:0,他引:1
S. Sisavath A. Al-Yaaruby C. C. Pain R. W. Zimmerman 《Pure and Applied Geophysics》2003,160(5-6):1009-1022
— Experimental observations show that flow through a fracture decreases more rapidly than the cube of the mean aperture (Cook, 1992). In order to provide a possible explanation of these experimental findings, we study creeping flow through a fracture of varying aperture that is symmetric about its midplane, using the power series of the stream function obtained by Van Dyke (1987) for low Reynolds numbers. For the case of sinusoidally-varying walls, a simple expression relating the effective hydraulic aperture of the channel to the mean aperture and to the amplitude and wavelength of the sinusoidal wall profiles is obtained. Comparison is made to previous studies (Kitanidis and Dykaar, 1997) and to finite element calculations, and good agreement is obtained. The effect of fracture closure is then modelled as a decrease of the mean aperture without a change in the roughness. A power law relationship can be obtained between the flowrate and the mean aperture, with an exponent as high as 10, thus providing a potential mechanistic explanation of the experimental findings of Pyrak-Nolte et al. (1987). 相似文献
18.
This article examines the required spatial discretization perpendicular to the fracture-matrix interface (FMI) for numerical simulation of solute transport in discretely fractured porous media. The discrete-fracture, finite-element model HydroGeoSphere ( Therrien et al. 2005 ) and a discrete-fracture implementation of MT3DMS ( Zheng 1990 ) were used to model solute transport in a single fracture, and the results were compared to the analytical solution of Tang et al. (1981) . To match analytical results on the relatively short timescales simulated in this study, very fine grid spacing perpendicular to the FMI of the scale of the fracture aperture is necessary if advection and/or dispersion in the fracture is high compared to diffusion in the matrix. The requirement of such extremely fine spatial discretization has not been previously reported in the literature. In cases of high matrix diffusion, matching the analytical results is achieved with larger grid spacing at the FMI. Cases where matrix diffusion is lower can employ a larger grid multiplier moving away from the FMI. The very fine spatial discretization identified in this study for cases of low matrix diffusion may limit the applicability of numerical discrete-fracture models in such cases. 相似文献
19.
—For small-scale microseismic events, the source sizes provided by shear models are unrealistically large when compared to visual observations of rock fractures near underground openings. A detailed analysis of the energy components in data from a mine-by experiment and from some mines showed that there is a depletion of S-wave energy for events close to the excavations, indicating that tensile cracking is the dominant mechanism in these microseismic events.¶In the present study, a method is proposed to estimate the fracture size from microseismic measurements. The method assumes tensile cracking as the dominant fracture mechanism for brittle rocks under compressive loads and relates the fracture size to the measured microseismic energy. With the proposed method, more meaningful physical fracture sizes can be obtained and this is demonstrated by an example on data from an underground excavation with detailed, high-quality microseismic records. 相似文献
20.
To more accurately predict the migration behavior of pollutants in porous media, we conduct laboratory scale experiments and model simulation. Aniline (AN) is used in one-dimensional soil column experiments designed under various media and hydrodynamic conditions. The advection-dispersion equation (ADE) and the continuous-time random walk (CTRW) were used to simulate the breakthrough curves (BTCs) of the solute transport. The results show that the media and hydrodynamic conditions are two important factors affecting solute transport and are related to the degree of non-Fickian transport. The simulation results show that CTRW can more effectively describe the non-Fickian phenomenon in the solute transport process than ADE. The sensitive parameter in the CTRW simulation process is , which can reflect the degree of non-Fickian diffusion in the solute transport. Understanding the relationship of with velocity and media particle size is conducive to improving the reactive solute transport model. The results of this study provide a theoretical basis for better prediction of pollutant transport in groundwater. 相似文献