共查询到20条相似文献，搜索用时 31 毫秒
1.
2.
裂隙岩体中非饱和渗流与运移的概念模型及数值模拟 总被引：12，自引：2，他引：12
探讨了裂隙岩体中非饱和地下水渗流与溶质运移的几种概念模型的构造及数值模拟问题 ,如裂隙网络模型、连续体模型、等效连续体模型、双孔隙度 (单渗透率 )模型、双渗透率模型、多组份连续体模型等。在裂隙岩体中 ,非饱和地下水的渗流可能只局限于岩体中的岩石组份、或裂隙网络 ,也可能在裂隙和岩石中同时发生 ;对前一种情形只需考虑单一连续体中的流动 ,而后一种情况则需要包括地下水在岩石和裂隙之间的交换。岩体中的裂隙网络往往是溶质运移的主要通道 ;但当溶质在裂隙与岩石之间的渗透和扩散是重要的运移机制时 ,就需要考虑岩石与裂隙界面处的溶质交换。为了模拟岩石与裂隙之间地下水和溶质的交换 ,就需要了解岩石与裂隙之间相互作用的模式和范围 ,使得这类问题的概念模型较单一连续体模型多了一层不确定性、其数值模拟也变得更为困难。因为在实际问题中不易、甚至根本不能判别非饱和渗流的实际形态 ,具体采用哪种模型主要取决于分析的目的和对现场数据的掌握程度。不论哪种模型都会受到模型及参数不确定性的影响 ,因此必须考虑与其他辅助模型的比较. 相似文献
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.
Kai Witthüser Dalinaidu Arnepalli Devendra N. Singh 《Geotechnical and Geological Engineering》2006,24(2):325-334
Contaminant transport through fractured rock mass is predominated by diffusion. This is due to the continuous interaction
of the mobile water present in the fracture network and relatively immobile pore water, which is adsorbed on the surface and
in the rock matrix itself. Even though the advective flow through the fracture network is high, besides sorption of rock mass,
the diffusive exchange into the rock mass leads to significant retardation of contaminant transport. Hence, for describing
contaminant transport in fractured rock mass, more precisely, the effect of retardation attributed to the matrix diffusion
must be taken in account. With this in view, a methodology, which can be employed for determination of the diffusion characteristics
of the rock mass, has been developed and its details are presented in this paper. Validation of the methodology has been demonstrated
with the help of Archie’s law. 相似文献
5.
Current designs for nuclear-waste repositories rely primarily upon subsurface geologic barriers for long-term containment. Because water and air are generally considered to be the mechanisms most likely to transport radioactivity to the surface environment, flow and transport models are important tools in repository assessment. Most models assume that the geologic medium can be treated as a continuum. A substantial body of recent work has focused on applying these models to difficult-to-solve problems, such as the simulation of variably dense or variably saturated flow and transport, large and complex flow systems, sharp solute concentration fronts, and fractured rock systems. The complex chemical interactions between the transport fluid and solid particles within the system have been analyzed using geochemical flow models, most of which assume that the system is at chemical equilibrium. The role of colloids in contaminant transport is a relatively new area of research. The large-scale effects of small-scale variability within the geologic system have been the subject of intense investigation. Inherent limitations of the continuum approach have prompted the design of models in which the flow occurs in discrete fractures. The difficulty and complexity of simulating transport has led to the development of network transport models, which represent the flow field as a series of 1-D path segments. The widespread use of models for prediction and analysis has prompted investigations of their reliability and relative merits. 相似文献
6.
Hybrid analytical and finite element numerical modeling of mass and heat transport in fractured rocks with matrix diffusion 总被引：1，自引：0，他引：1
Christopher I. McDermott Robert Walsh Ralph Mettier Georg Kosakowski Olaf Kolditz 《Computational Geosciences》2009,13(3):349-361
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. 相似文献
7.
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. 相似文献
8.
对交叉裂隙渗流传质特性的定量描述是研究整个裂隙网络渗透传质特性的基础。为真实模拟水流及溶质在三维交叉裂隙中的运移过程,首先通过三维轮廓仪获取天然岩石裂隙表面的形貌数据,再应用三维重构技术生成相应的三维交叉裂隙模型,随后求解Navier-Stokes方程,假定溶质运移满足Fick定律,模拟水流和溶质在三维交叉裂隙中的运移过程。通过对比粗糙裂隙模型与平行平板模型的模拟结果发现：粗糙度对流体的分布及流动状态存在显著的影响;不同进、出口工况下的流体流动及溶质运移状态亦表明：裂隙交叉的几何形貌会显著地影响溶质混合行为。这些结果表明,目前被广泛采用的平行平板模型在评估岩体内特别是交叉口的物质运移特性时将导致较大的偏差,在将来的研究中有必要针对裂隙交叉口的几何特征建立修正的模型以提高评估的准确性。 相似文献
9.
Sung-Hoon Ji Kyung Woo Park Doo-Hyun Lim Chunsoo Kim Kyung Su Kim William Dershowitz 《Hydrogeology Journal》2012,20(7):1341-1353
The development and implementation of a hybrid discrete fracture network/equivalent porous medium (DFN/EPM) approach to groundwater flow at the Gyeong-Ju low- and intermediate-level radioactive waste (LILW) disposal site in the Republic of Korea is reported. The geometrical and hydrogeological properties of fractured zones, background fractures and rock matrix were derived from site characterization data and implemented as a DFN. Several DFN realizations, including the deterministic fractured zones and the stochastic background fractures, whose statistical properties were verified by comparison with in-situ fracture and hydraulic test data, were suggested, and they were then upscaled to continuums using a fracture tensor approach for site-scale flow simulations. The upscaled models were evaluated by comparison to in-situ pressure monitoring data, and then used to simulate post-closure hydrogeology for the LILW facility. Simulation results demonstrate the importance of careful characterization and implementation of fractured zones. The study highlighted the importance of reducing uncertainty regarding the properties and variability of natural background fractures, particularly in the immediate vicinity of repository emplacement. 相似文献
10.
Witherspoon, P.A. and Gale, J.E., 1977. Mechanical and hydraulic properties of rocks related to induced seismicity. Eng. Geol., 11(1): 23–55.The mechanical and hydraulic properties of fractured rocks are considered with regard to the role they play in induced seismicity. In many cases, the mechanical properties of fractures determine the stability of a rock mass. The problems of sampling and testing these rock discontinuities and interpreting their non-linear behavior are reviewed. Stick slip has been proposed as the failure mechanism in earthquake events. Because of the complex interactions that are inherent in the mechanical behavior of fractured rocks, there seems to be no simple way to combine the deformation characteristics of several sets of fractures when there are significant perturbations of existing conditions. Thus, the more important fractures must be treated as individual components in the rock mass.In considering the hydraulic properties, it has been customary to treat a fracture as a parallel-plate conduit and a number of mathematical models of fracture systems have adopted this approach. Non-steady flow in fractured systems has usually been based on a two-porosity model, which assumes the primary (intergranular) porosity contributes only to storage and the secondary (fracture) porosity contributes only to the overall conductivity. Using such a model, it has been found that the time required to achieve quasi-steady state flow in a fractured reservoir is one or two orders of magnitude greater than it is in a homogeneous system. In essentially all of this work, the assumption has generally been made that the fractures are rigid.However, it is clear from a review of the mechanical and hydraulic properties that not only are fractures easily deformed but they constitute the main flow paths in many rock masses. This means that one must consider the interaction of mechanical and hydraulic effects. A considerable amount of laboratory and field data is now available that clearly demonstrates this stress-flow behavior. Two approaches have been used in attempting to numerically model such behavior: (1) continuum models, and (2) discrete models. The continuum approach only needs information as to average values of fracture spacing and material properties. But because of the inherent complexity of fractured rock masses and the corresponding decrease in symmetry, it is difficult to develop an equivalent continuum that will simulate the behavior of the entire system. The discrete approach, on the other hand, requires details of the fracture geometry and material properties of both fractures and rock matrix. The difficulty in obtaining such information has been considered a serious limitation of discrete models, but improved borehole techniques can enable one to obtain the necessary data, at least in shallow systems. The possibility of extending these methods to deeper fracture systems needs more investigation. Such data must be considered when deciding whether to use a continuum or discrete model to represent the interaction of rock and fluid forces in a fractured rock system, especially with regard to the problem of induced seismicity. When one is attempting to alter the pressure distribution in a fault zone by injection or withdrawal of fluids, the extent to which this can be achieved will be controlled in large measure by the behavior of the fractures that communicate with the borehole. Since this is essentially a point phenomenon, i.e., the changes will propagate from a relatively small region around the borehole, the use of a discrete model would appear to be preferable. 相似文献
11.
Stochastic continuum modeling of flow and transport in a crystalline rock mass: Fanay-Augères,France, revisited 总被引：1，自引：0，他引：1
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.
Groundwater flow and solute transport through fractured rock is highly responsive to the hydraulic anisotropy and heterogeneity that are specific to every major fracture. A major fracture is modeled as the combination of some primal master fractures and several splay fractures that branch out from primal master fractures: step structures (or jog parts). Step structures are commonly observed along a major fracture on various scales. Master fractures were formed and developed by shear movement while some splay fractures were formed by extension normal to their wall. This difference in fracturing process may lead to a permeability difference between master fractures and splay fractures which seems to be one of the major factors controlling flow and solute transport through the fracture networks due to its hydraulic anisotropic and heterogeneous features. This study is composed of two major components: (1) identification and characterization of a step structure from borehole data; (2) evaluation of effect of some idealized step structures on breakthrough curve by numerical simulations. The fracture data of four 1000-m boreholes were used to make clear fracture patterns in the Tono area of Japan. Some major fractures were identified using stereographic projection technique. On the basis of these results, several idealized models of a major fracture having a step was constructed for the numerical study. The obtained results from numerical simulations clearly imply that geometry of step structure plays an important role in flow and transport through the fracture networks. 相似文献
13.
Groundwater flow in fractured rocks is modeled using a coupled model based on the domain decomposition method. In the model, the fractured porous medium is divided into two non-overlapping sub-domains. One is the rock matrix, in which the medium is described using a continuum model. The other consists of deep fractures and fissure zones, where the medium is described using a discrete fracture network (DFN) model. The two models are coupled through the continuity of the hydraulic heads and fluxes on the common boundaries. The coupled model is used to simulate groundwater flow in a hydropower station. The results show that the model simulates groundwater levels that are in agreement with the measured groundwater levels. Furthermore, the model’s parameters relating to deep fractures and fissure zones are verified by comparing three different models (the continuum model, coupled model, and DFN model). The results show that the coupled model can capture and duplicate the hydrogeological conditions in the study domain, whereas the continuum model overestimates and the DFN model underestimates the measured hydraulic heads. A sensitivity analysis shows that fracture aperture has a considerable effect on the groundwater level. So, when the fracture aperture is large, the coupled model or DFN model is more appropriate than the continuum model in the fracture domain. 相似文献
14.
This paper proposes a three-dimensional coupled hydrothermal model for fractured rock based on the finite-discrete element method to simulate fluid flow and heat transport. The 3D coupled hydrothermal model is composed of three main parts: a heat conduction model for the rock matrix, a heat transfer model for the fluid in the fractures (including heat conduction and heat convection), and a heat exchange model between the rock matrix and the fluid in the fractures. Four examples with analytical solutions are provided to verify the model. A heat exchange experiment of circulating water in a cylindrical granite sample with one fracture is simulated. The simulation results agree well with the experimental results. The effects of the fracture aperture, fluid viscosity, and pressure difference on the heat exchange between the fluid and rock are studied. Finally, an application concerned with heat transport and fluid flow in fractured rock is presented. The simulation results indicate that the 3D fully coupled hydrothermal model can capture the fluid flow and temperature evolution of rocks and fluids. 相似文献
15.
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. 相似文献
16.
Fracture networks play a more significant role in conducting fluid flow and solute transport in fractured rock masses, comparing with that of the rock matrix. Accurate estimation of the permeability of fracture networks would help researchers and engineers better assess the performance of projects associated with fluid flow in fractured rock masses. This study provides a review of previous works that have focused on the estimation of equivalent permeability of two-dimensional (2-D) discrete fracture networks (DFNs) considering the influences of geometric properties of fractured rock masses. Mathematical expressions for the effects of nine important parameters that significantly impact on the equivalent permeability of DFNs are summarized, including (1) fracture-length distribution, (2) aperture distribution, (3) fracture surface roughness, (4) fracture dead-end, (5) number of intersections, (6) hydraulic gradient, (7) boundary stress, (8) anisotropy, and (9) scale. Recent developments of 3-D fracture networks are briefly reviewed to underline the importance of utilizing 3-D models in future research. 相似文献
17.
考虑裂隙的连通率、间距、孔隙基质和裂隙材料在表征单元（REV）中的体积分数,并假定双重孔隙-裂隙介质的等效内摩擦角保持常数,而等效黏聚力是固有黏聚力、等效塑性应变、基质吸力、溶质浓度及温度的函数,提出了一种在热-水-应力-迁移耦合条件下确定表征单元内任一平面上等效的黏聚力及内摩擦角的方法。针对一个假定的位于非饱和双重孔隙-裂隙岩体中的高放废物地质处置模型进行了数值模拟及分析。结果表明,基质吸力对于等效黏聚力的增强作用大于等效塑性应变和溶质浓度的减弱作用,使得等效黏聚力得到了提高,故减少了围岩中的塑性区;由此岩体应力、孔（裂）隙水压力及流速、孔（裂）隙溶质浓度的分布及量值也发生相应的改变。 相似文献
18.
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. 相似文献
19.
Two-dimensional reactive transport modeling of the Maqarin Eastern Springs site, a natural analogue for the alteration of a fractured limestone by high-pH Portland cement waters, has been performed using the CrunchFlow code. These 2D calculations included transport by advection–dispersion–diffusion along a single fracture and diffusion in the wall rock. Solute transport was coupled to mineral dissolution and precipitation. A limited sensitivity analysis evaluated the effect of different values of primary mineral surface areas, flow velocity and sulfate concentration of the inflowing high-pH solution.Major secondary minerals include ettringite–thaumasite, C–S–H/C–A–S–H and calcite. C–S–H/C–A–S–H precipitation is controlled by the dissolution of primary silicates. Ettringite precipitation is controlled by diffusion of sulfate and aluminum from the wall rock to the fracture, with aluminum provided by the dissolution of albite. Calcite precipitation is controlled by diffusion of carbonate from the wall rock. Extents of porosity sealing along the fracture and in the fracture-wall rock interface depend on assumptions regarding flow velocity and composition of the high-pH solution. The multiple episodes of fracture sealing and reactivation evidenced in the fracture infills were not included in the simulations. Results can qualitatively reproduce the reported decrease in porosity in the fractures and in the wall rock next to the fractures. Instances of porosity increase next to fractures caused by carbonate dissolution were not reproduced by the calculations. 相似文献
20.
米尺度裂隙岩体模型水流-传热试验的数值模拟分析 总被引：1，自引：0，他引：1
为了研究高放射性核废物地下处置库近场的水流-传热耦合问题,采用国内高放废物地下处置库预选场址--甘肃北山地区的花岗岩石块体,加工组合成米尺度的规则裂隙岩体模型,设置边界热源和裂隙水流,试验模拟裂隙水水流与传热之间的相互作用。作为该室内模型试验的前期理论研究,采用等效孔隙介质数值模型,着重分析了裂隙开度、裂隙流量和热源功率对流场和温度场的影响。在设定条件下,计算分析表明：热传导和裂隙水水流由热源作用初期的不耦合很快转化为耦合;不流动的裂隙水主要表现为热存储和热传导,而流动的裂隙水还引起流动传热和水与岩石之间的对流换热,使岩体温度场明显不同于单纯热传导的情况;如果保持裂隙水流量不变,则裂隙开度的变化对水流-传热影响不大;如果保持裂隙水流速不变,则裂隙开度的变化对水流-传热影响显著;热源功率越大,通过裂隙水的热流量越大,裂隙水压强越大,而当温度超过100 ℃时,裂隙水会因汽化而压强显著增大;加热7 d时,热量的输入和输出几乎相等,裂隙水流带走的热量接近热源供给的热量,模型系统基本达到了热平衡。 相似文献