首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
Analysis of contaminant transport through fractured crystalline rocks has received considerable attention, particularly with regard to subsurface nuclear waste repositories. Most of the studies have employed the dual continuum approach, with the fractures and the rock matrix as the two continuums, assuming that fractures control the overall conductivity of the rock and the porous matrix just provides storage. However, field observations of rock fractures have shown that the real situation can be very complex. Based on some recent investigations, it has been reported that the portion of the rock matrix adjacent to many open fractures is physically and chemically altered. These alterations, referred to as the fracture skin, can have different sorption and diffusion properties compared to those of the undisturbed rock matrix and this may influence the transport of solutes through such formations. In the present study, a numerical model is developed to simulate conservative solute transport in a fractured crystalline rock formation using the triple continuum approach ?? with the fracture, fracture skin and the rock matrix as the three continuums. The model is solved using a fully implicit finite difference scheme. Contaminant migration in the fractured formation with and without skin has been simulated. It is observed that contaminant penetration along the fracture is enhanced at large flow velocities. The effect of flow velocity on conservative solute transport is investigated for different fracture apertures and fracture skin thicknesses. The influence of flow velocity on contaminant transport is demonstrated to be more with change in fracture aperture than with change in skin thickness.  相似文献   

2.
Quantification of mass and heat transport in fractured porous rocks is important to areas such as contaminant transport, storage and release in fractured rock aquifers, the migration and sorption of radioactive nuclides from waste depositories, and the characterization of engineered heat exchangers in the context of enhanced geothermal systems. The large difference between flow and transport characteristics in fractures and in the surrounding matrix rock means models of such systems are forced to make a number of simplifications. Analytical approaches assume a homogeneous system, numerical approaches address the scale at which a process is operating, but may lose individual important processes due to averaging considerations. Numerical stability criteria limit the contrasts possible in defining material properties. Here, a hybrid analytical–numerical method for transport modeling in fractured media is presented. This method combines a numerical model for flow and transport in a heterogeneous fracture and an analytical solution for matrix diffusion. By linking the two types of model, the advantages of both methods can be combined. The methodology as well as the mathematical background are developed, verified for simple geometries, and applied to fractures representing experimental field conditions in the Grimsel rock laboratory.  相似文献   

3.
A new semi-analytical solution for the transport of a conservative contaminant species in a fractured medium having a regular two- or three-dimensional fracture network is presented. The application of the technique and some of the practical implications arising from an examination of contaminant migration in fractured systems is discussed. Particular consideration is given to the effects of Darcy velocity, fracture spacing, matrix porosity, dispersivity and the mass of contaminant available for transport. The implications of uncertainty with respect to fracture opening size and ground-water velocity is also discussed and it is shown that provided one can obtain a reasonable estimate of the hydraulic gradient and hydraulic conductivity for the rock mass, uncertainty regarding the magnitude of the opening size and the groundwater velocity does not have a significant effect on predicted contaminant migration for the class of problems being considered.  相似文献   

4.
Numerical modeling of stress effects on solute transport in fractured rocks   总被引:1,自引:0,他引:1  
The effects of stress/deformation on fluid flow and contaminant transport in fractured rocks is one of the major concerns for performance and safety assessments of many subsurface engineering problems, especially radioactive waste disposal and oil/gas reservoir fields. However, very little progress has been made to study this issue due to difficulties in both experiments and numerical modeling. The objective of this study is to systematically investigate the influence of stress on solute transport in fractured rocks for the first time, considering different stress and hydraulic pressure conditions. A hybrid approach combining discrete element method (DEM) for stress-flow simulations and a particle tracking algorithm is developed. The impact of matrix diffusion (diffusion of molecular size solutes in and out of the rock matrix, and sorption onto the surface of micropores in rock matrix) is also included. The numerical results show that stress not only significantly changes the solute residence time through the fracture networks, but also changes the solute travel paths. Matrix diffusion plays a dominant role in solute transport when the hydraulic gradient is small, which is often encountered in practice.  相似文献   

5.
A double-porosity model is used to describe the oscillatory gas motion and associated contaminant transport induced by cyclical variations in the barometric pressure at the surface of a fractured porous medium. Flow along the fractures and within the permeable matrix blocks is locally one-dimensional. The interaction between fractures and blocks includes seepage of fluid as well as diffusion of contaminant. To guard against artificial numerical diffusion, the FRAM filtering remedy and methodology of Chapman is used in calculating the advective fluxes along fractures and within blocks. The entire system of equations, including the fracture-matrix interaction terms, is solved by a largely implicit non-iterative algorithm which remains stable and conservative even when the computational time step is large compared to the cross-block transit time of pressure waves. The numerical accuracy is tested by comparison with exact solutions for oscillatory and unidirectional flows, some of which include diffusion interaction between the fracture and the matrix. The method is used to estimate the rate of vertical transport of radioactive gases through the rubblized chimney produced by an underground nuclear explosion.  相似文献   

6.
7.
A two-dimensional boundary integral method to analyse the flow of contaminant in fractured media having a two- or three-dimensional orthogonal fracture network is presented. The method assumes that the fractures provide the paths of least resistance for transport of contaminants while the matrix, because of its low permeability, acts as ‘storage blocks’ into which the contaminant diffuses. Laplace transform is used to eliminate the time variable in the governing equation in order to facilitate the formulation of a boundary integral equation in the Laplace transform space. Conventional boundary element techniques are applied to solve for the contaminant concentrations at specified locations in the spatial domain. The concentration in the time domain is then obtained by using an efficient inversion technique developed by Talbot. The method is able to analyse the behaviour of waste repositories which have diminishing concentration due to the mass transport of the contaminant into the surrounding fractured media.  相似文献   

8.
裂隙渗流会引起裂隙周围岩体中的温度场变化,在低温岩体中其影响更为明显;此外,裂隙水与周围低温岩石介质发生热交换会引起裂隙中的水冰相变过程发生,而裂隙水冻结将阻碍裂隙渗流,引起裂隙渗流场的变化。因此,低温下的裂隙岩体水-热相互作用是一个强耦合过程。考虑裂隙中的水冰相变过程和渗流作用,建立了低温冻结条件下裂隙岩体水-热耦合模型;以冻结法施工为例,考察了低温冻结过程中裂隙水渗流对裂隙冻结交圈的影响。研究结果表明:由于裂隙渗流的存在,距裂隙较远处岩石先冻结,裂隙冻结所需时间远大于周围岩石;裂隙宽度和裂隙水压力差都会影响冻结交圈时间,裂隙越宽、水压力差越大,裂隙冻结需要时间越长;随着冻结时间的推进,裂隙水渗流速度逐渐降低,当裂隙冻结后裂隙渗流停止。最后通过构建随机裂隙网络模型,利用所建立的水-热耦合模型考察了裂隙网络渗流对冻结交圈的影响,说明了在冻结法施工中考虑裂隙的重要性。  相似文献   

9.
10.
The unsaturated fractured volcanic deposits at Yucca Mountain in Nevada, USA, have been intensively investigated as a possible repository site for storing high-level radioactive waste. Field studies at the site have revealed that there exist large variabilities in hydrological parameters over the spatial domain of the mountain. Systematic analyses of hydrological parameters using a site-scale three-dimensional unsaturated zone (UZ) flow model have been undertaken. The main objective of the sensitivity analyses was to evaluate the effects of uncertainties in hydrologic parameters on modeled UZ flow and contaminant transport results. Sensitivity analyses were carried out relative to fracture and matrix permeability and capillary strength (van Genuchten α) through variation of these parameter values by one standard deviation from the base-case values. The parameter variation resulted in eight parameter sets. Modeling results for the eight UZ flow sensitivity cases have been compared with field observed data and simulation results from the base-case model. The effects of parameter uncertainties on the flow fields were evaluated through comparison of results for flow and transport. In general, this study shows that uncertainties in matrix parameters cause larger uncertainty in simulated moisture flux than corresponding uncertainties in fracture properties for unsaturated flow through heterogeneous fractured rock.  相似文献   

11.
A stochastic discrete-fracture model was used by Cacas et al.a,b to interpret flow measurements and transport experiments in a fractured crystalline rock mass at Fanay-Augères. They considered continuum models to be incapable of properly interpreting small-scale measurements or tracer tests in fractured systems, which, in their view, require three-dimensional modeling of numerous discrete channels; in their opinion, continuum modeling applies only to average flow on a relatively large scale. Cacas et al. considered their discrete fracture model to have been validated by its demonstrated ability to reproduce selected experimental results. In this paper, flow and transport at Fanay-Augères are modeled by viewing the fractured rock as a stochastic continuum in a manner originally proposed by Neumanc,d. The stochastic continuum approach obviates the need for detailed information about fracture geometry or assumptions about how individual fractures control flow and transport. All it requires is the delineation of a few dominant features, which can be embedded into the stochastic continuum model as heterogeneous porous slabs. Though a fault zone has been identified at the Fanay-Augères experimental site, it has been modeled neither by Cacas et al. nor in this paper. In fact, in this paper, a larger selection of experimental results than those considered by Cacas et al. are reproduced merely by modeling the rock as a statistically homogeneous continuum in two dimensions. These results demonstrate that a continuum approach may be well suited for the analysis of flow and transport in fractured rock. This does not constitute a validation of the continuum approach, just as the results of Cacas et al. fall short of validating the discrete fracture approach. Instead, the two sets of results illustrate jointly the well-established principle that an open system, especially one as complex as fractured hydrogeologic environments tend to be, cannot be described uniquely on the basis of sparse data and need not be described in great detail to capture its salient behavior by a model.a Cacas MC, Ledoux E, de Marsily G, Barbreau A, Calmels P, Gaillard B, Margritta R (1990a) Modelling fracture flow with a stochastic discrete fracture network: calibration and validation. 1. The flow model. Water Resour Res 26(3):479–489b Cacas MC, Ledoux E, de Marsily G, Barbreau A, Calmels P, Gaillard B, Margritta R (1990b) Modelling fracture flow with a stochastic discrete fracture network: calibration and validation. 2. The transport model. Water Resour Res 26(3):491–500c Neuman SP (1987) Stochastic continuum representation of fractured rock permeability as an alternative to the REV and fracture network concepts, in Rock Mechanics. In: Farmer IW, Daemen JJK, Desai CS, Glass CE, Neuman SP (eds) Proceedings of the 28th U.S. Symposium, Tucson, Arizona. Balkema, Rotterdam, pp 533–561d Neuman SP (1988) A proposed conceptual framework and methodology for investigating flow and transport in Swedish crystalline rocks. SKB Swedish Nuclear Fuel and Waste Management Co., Stockholm, September, Arbetsrapport 88–37, 39 pp  相似文献   

12.
裂隙岩体中非饱和渗流与运移的概念模型及数值模拟   总被引:12,自引:2,他引:12  
探讨了裂隙岩体中非饱和地下水渗流与溶质运移的几种概念模型的构造及数值模拟问题 ,如裂隙网络模型、连续体模型、等效连续体模型、双孔隙度 (单渗透率 )模型、双渗透率模型、多组份连续体模型等。在裂隙岩体中 ,非饱和地下水的渗流可能只局限于岩体中的岩石组份、或裂隙网络 ,也可能在裂隙和岩石中同时发生 ;对前一种情形只需考虑单一连续体中的流动 ,而后一种情况则需要包括地下水在岩石和裂隙之间的交换。岩体中的裂隙网络往往是溶质运移的主要通道 ;但当溶质在裂隙与岩石之间的渗透和扩散是重要的运移机制时 ,就需要考虑岩石与裂隙界面处的溶质交换。为了模拟岩石与裂隙之间地下水和溶质的交换 ,就需要了解岩石与裂隙之间相互作用的模式和范围 ,使得这类问题的概念模型较单一连续体模型多了一层不确定性、其数值模拟也变得更为困难。因为在实际问题中不易、甚至根本不能判别非饱和渗流的实际形态 ,具体采用哪种模型主要取决于分析的目的和对现场数据的掌握程度。不论哪种模型都会受到模型及参数不确定性的影响 ,因此必须考虑与其他辅助模型的比较.  相似文献   

13.
Precipitation-dissolution reactions are important for a number of applications such as isotopic tracer transport in the subsurface. Analytical solutions have been developed for tracer transport in both single-fracture and multiple-fracture systems associated with these reactions under transient and steady-state transport conditions. These solutions also take into account advective transport in fractures and molecular diffusion in the rock matrix. For studying distributions of disturbed tracer concentration (the difference between actual concentration and its equilibrium value), effects of precipitation-dissolution reactions are mathematically equivalent to a “decay” process with a decay constant proportional to the corresponding bulk reaction rate. This important feature significantly simplifies the derivation procedure by taking advantage of the existence of analytical solutions for tracer transport associated with radioactive decay in fractured rock. It is also useful for interpreting tracer breakthrough curves, because the impact of a decay process is relatively easy to analyze. Several illustrative examples are presented, which show that the results are sensitive to fracture spacing, matrix diffusion coefficient (fracture surface area), and bulk reaction rate (or “decay” constant), indicating that the relevant flow and transport parameters may be estimated by analyzing tracer signals.  相似文献   

14.
15.
Modelling of radionuclide transport in fractured media is a primary task for safety evaluation of a deep nuclear waste repository. A performance assessment (PA) model has been derived from site characterization data with the aim of improving confidence for quantifying transport of sorbing radionuclides at a safety time scale of 106 y. The study was conducted on a 200?×?200?×?200 m semi-synthetic fractured block, providing a realistic system derived from the Äspö Hard Rock Laboratory (Sweden) dataset. The block includes 5,632 fractures ranging from 0.5 to 100 m in length and a heterogeneous matrix structure (fracture coating, gouge, mylonite, altered and non-altered diorite). The PA model integrates steady-state flow conditions and transport of released radionuclides during the safety time scale. An original simulation method was developed involving Eulerian flow and transport within fracture planes with a mixed hybrid finite element scheme and a semi-analytical source term to account for heterogeneous matrix diffusion. Total mass flux of radionuclides (conservative to strongly sorbing) was computed. A method to simplify the system was demonstrated, leading to a major path of 12 fractures. This is mainly due to the low connectivity of the fracture network. Matrix diffusion and sorption proved to have major impact on block retention properties for PA conditions.  相似文献   

16.
分析了运移参数的变化对GCL中污染物运移的影响。从衬垫底部浓度比、对流和扩散通量、污染物累积运移量等方面,对污染物在GCL中的运移进行了计算。结果表明,减小渗透系数和扩散系数可以有效减缓衬垫底部浓度比的增加速度。渗透系数对对流通量的大小、扩散通量的变化速度的影响较大;扩散系数主要影响最大扩散通量。当渗透系数较大时,可以忽略扩散累积运移量,用对流累积运移量代替总累积运移量引起的误差很小。当渗透系数较小时,忽略扩散运移量可能产生较大的误差,尤其是在扩散系数也较大的情况下误差更加明显。   相似文献   

17.
In fractured-rock aquifers, the geometric and hydraulic properties of the fractures commonly have a dominant influence on transport. Tracer tests are often used to estimate directly the gross transport properties of a fractured rock mass. The prospects for understanding characteristics of the heterogeneities in a fractured porous medium were explored from evidence provided by tracer experiments. The approach was to simulate flow and transport on a large set of prescribed fracture networks in a two-dimensional homogeneous permeable medium, thus generating synthetic tracer test data. The fracture orientation, aperture, spacing and network geometry were systematically altered from one case to the next. A classification scheme was devised for the tracer breakthrough curves using principal component analysis and this classification was linked to the fracture pattern properties. Even under highly simplified and controlled conditions, quite different fracture patterns can produce very similar breakthrough curves. The classification scheme thus demonstrates that a single breakthrough curve cannot reveal the fracture geometry with any precision. However, the scheme provided a methodology for rejecting geometric properties that do not belong to the fracture pattern under investigation, thus reducing the uncertainty in fracture geometry.  相似文献   

18.
岩体是由裂隙和其间的岩块组成。当岩块处于近乎隔水状态时,岩体中的渗流可看作为裂隙网络流。本文通过分析岩体裂隙网络渗流场、应力场以及它们之间的相互力学关系,提出了岩体渗流场与应力场耦合的裂隙网络模型及数值计算方法。  相似文献   

19.
岩溶地区隧道裂隙水突出力学机制研究   总被引:16,自引:4,他引:12  
岩溶地下水是诱发隧道发生突(涌)水地质灾害的主导因素之一,岩溶裂隙水对隧道围岩的危害越来越成为岩溶地区隧道建设中的热点研究问题之一。基于目前岩溶裂隙水突出机制研究现状,运用岩溶地质学、工程水力学和断裂力学相关理论分析岩溶地区隧道水岩相互作用机制,探讨了水岩相互作用对岩溶地区隧道施工发生突(涌)水的影响,揭示岩溶地区隧道裂隙水突出前后过程的力学机制。研究表明:岩溶隧道裂隙水突出是裂隙岩体在岩溶水及水压的持续作用下受施工外力干扰发生劈裂的结果,岩溶水和水压在裂隙岩体突水破坏过程中的力学作用主要体现在4个方面,即突水蓄势期岩溶水对裂隙岩体的软化溶蚀作用、水压对裂隙岩体的劈裂作用,突水失稳期水流的冲刷扩径作用、水压对突水量的动力控制作用。基于上述分析,以断裂力学、弹塑性力学理论为基础,提出岩溶隧道岩溶裂隙水突出的最小岩石防突厚度概念,推导了其半解析解表达式,并为工程实例所验证,其结果可为高风险岩溶地区隧道突水理论与防治措施提供参考。  相似文献   

20.
In situ measurements of deformations, stresses, and closure of fractures, affecting water inflow following coal mining, are challenging due to the inaccessibility of fractured rock. In this paper, the authors studied the closure process of the fractured rock mass with the cover stress re-establishment based on a theoretical analysis and a scale model testing. A quantitative analysis is used to study the fracture distribution in the fractured zone. A function to describe a fracture aperture distribution in the fractured zone is proposed, which takes into account the curvature and thickness of the fractured rock. The theoretical analysis and a scale model testing both indicate that the cover stress re-establishment with mining distance increasing and the relationship between the fracture closure and cover stress re-establishment both satisfy a logarithmic function. The scale model test also shows the following features: (1) the fracture ratio (which is the fracture area divided by the total area of fracture and intact rock with a unit width in the vertical or horizontal direction) in the lower part of the fractured rock mass is greater than that in the upper part; (2) the initially fast decreased of fracture ratios is then followed by a slower decrease during the cover stress re-establishment process; (3) in the upper part of the rock mass, the vertical directional fractures with small apertures are being closed with cover stress re-establishment, which indicates an increase in the water resistance reducing the seepage from these parts of the fractured zone. This study improves the general understanding of the fracture closure process and cover stress re-establishment in the fractured rock mass after coal mining ceased, and provides a theoretical basis for water resource protection in case of underground coal mining.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号