首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 250 毫秒
1.
离散裂隙渗流方法与裂隙化渗透介质建模   总被引:4,自引:1,他引:4  
流体渗流模拟的连续介质方法通常适用于多孔地质体,并不一定适用于裂隙岩体,由于裂隙分布及其特征与孔隙差异较大。若流体渗流主要受裂隙的控制,对于一定尺寸的裂隙岩体,多孔介质假设则较难刻划裂隙岩体的渗流特征。离散裂隙渗流方法不但可直接用于模拟裂隙岩体非均质性和各向异性等渗流特征,而且可用其确定所研究的裂隙岩体典型单元体及其水力传导(渗透)张量大小。主要讨论了以下问题:(1)饱和裂隙介质中一般的离散流体渗流模拟;(2)裂隙岩体中的REV(典型单元体)及其水力传导(渗透)张量的确定;(3)利用离散裂隙网络流体渗流模型研究裂隙方向几何参数对水力传导系数和REV的影响;(4)在二维和三维离散裂隙流体渗流模型中对区域大裂隙和局部小裂隙的处理方法。调查结果显示离散裂隙流体渗流数学模型可用来评价不同尺度上的裂隙岩体的水力特征,以及裂隙方向对裂隙化岩体的水力特征有着不可忽视的影响。同时,局部小裂隙、区域大裂隙应当区别对待,以便据其所起的作用及水力特征,建立裂隙化岩体相应的流体渗流模型。  相似文献   

2.
INTRODUCTIONGroundwaterorfluidflowmodelinginfracturedrocksisacomplicatedtheoreticalandappliedtopic.Boththeoreticallyandoperationally ,itisimportantinmanyfieldssuchasgeologicalandhydrogeologicalengineering ,environmentalengineeringandpetroleumengineerin…  相似文献   

3.
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.  相似文献   

4.
裂隙岩体渗透系数确定方法研究   总被引:3,自引:0,他引:3  
裂隙岩体渗透系数以及渗透主方向的确定对研究岩体渗透性大小及各向异性具有重要意义。高放废物地质处置库介质岩体的渗透性能将直接影响其使用安全性。本文运用离散裂隙网络模拟的方法对我国高放废物处置库甘肃北山预选区3#钻孔附近裂隙岩体进行了渗透性质分析。通过对3#钻孔171.5~178.0m段压水试验数据的反演,标定了离散裂隙网络渗流模型中的裂隙渗透参数(导水系数T)。利用标定的离散裂隙网络模型对场区裂隙岩体进行了渗流模拟,确定了该区域裂隙岩体的渗流表征单元体(REV)的尺寸大小以及渗透主值和主渗透方向。运用离散裂隙网络模型计算得出的渗透主值的几何均值与现场压水试验计算结果较接近,证明了计算结果的有效性。  相似文献   

5.
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.  相似文献   

6.
Krzaczek  M.  Nitka  M.  Kozicki  J.  Tejchman  J. 《Acta Geotechnica》2020,15(2):297-324

The paper deals with two-dimensional (2D) numerical modelling of hydro-fracking (hydraulic fracturing) in rocks at the meso-scale. A numerical model was developed to characterize the properties of fluid-driven fractures in rocks by combining the discrete element method (DEM) with computational fluid dynamics (CFD). The mechanical behaviour of the rock matrix was simulated with DEM and the behaviour of the fracturing fluid flow in newly developed and pre-existing fractures with CFD. The changes in the void geometry in the rock matrix were taken into account. The initial 2D hydro-fracking simulation tests were carried out for a rock segment under biaxial compression with one injection slot in order to validate the numerical model. The qualitative effect of several parameters on the propagation of a hydraulic fracture was studied: initial porosity of the rock matrix, dynamic viscosity of the fracking fluid, rock strength and pre-existing fracture. The characteristic features of a fractured rock mass due to a high-pressure injection of fluid were realistically modelled by the proposed coupled approach.

  相似文献   

7.
岩体裂隙系统渗流场与应力场耦合模型   总被引:15,自引:0,他引:15  
岩体系统具有复杂的结构。一般认为,岩体系统是非均质各向异性不连续的多相介质体系。当岩体以裂隙为主,且其分布较密集时,可将岩体系统看作等效连续多相介质体系。本文运用等效连续介质理论,提出了两种岩体裂隙系统渗流场与应力场耦合模型:一是以渗透水压力与隙变形关系、应力与渗透系统数关系为基础,建立渗透系数张量计算公式,进而建立等效效连续介质渗流为数学模型。以裂隙岩体应变张量分析为基础,建立裂隙岩体效应力张量  相似文献   

8.
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  相似文献   

9.
A hybrid discrete‐continuum numerical scheme is developed to study the behavior of a hydraulic fracture crossing natural fractures. The fully coupled hybrid scheme utilizes a discrete element model for an inner domain, within which the hydraulic fracture propagates and interacts with natural fractures. The inner domain is embedded in an outer continuum domain that is implemented to extend the length of the hydraulic fracture and to better approximate the boundary effects. The fracture is identified to propagate initially in the viscosity‐dominated regime, and the numerical scheme is calibrated by using the theoretical plane strain hydraulic fracture solution. The simulation results for orthogonal crossing indicate three fundamental crossing scenarios, which occur for various stress ratios and friction coefficients of the natural fracture: (i) no crossing, that is, the hydraulic fracture is arrested by the natural fracture and makes a T‐shape intersection; (ii) offset crossing, that is, the hydraulic fracture crosses the natural fracture with an offset; and (iii) direct crossing, that is, the hydraulic fracture directly crosses the natural fracture without diversion. Each crossing scenario is associated with a distinct net pressure history. Additionally, the effects of strength contrast and stiffness contrast of rock materials and intersection angle between the hydraulic fracture and the natural fracture are also investigated. The simulations also illustrate that the level of fracturing complexity increases as the number and extent of the natural fractures increase. As a result, we can conclude that complex hydraulic fracture propagation patterns occur because of complicated crossing behavior during the stimulation of naturally fractured reservoirs. Copyright © 2017 John Wiley & Sons, Ltd.  相似文献   

10.
Generally, induced hydraulic fractures are generated by fluid overpressure and are used to increase reservoir permeability through forming interconnected fracture systems. However, in heterogeneous and anisotropic rocks, many hydraulic fractures may become arrested or offset at layer contacts under certain conditions and do not form vertically connected fracture networks. Mechanical layering is an important factor causing anisotropy in sedimentary layers. Hence, in this study, with a shale gas reservoir case study in the Longmaxi Formation in the southeastern Chongqing region, Sichuan Basin, we present results from several numerical models to gain quantitative insights into the effects of mechanical layering on hydraulic fracturing. Results showed that the fractured area caused by hydraulic fracturing indicated a linear relationship with the neighboring layer’s Young’s modulus. An increase of the neighboring layer’s Young’s modulus resulted in better hydraulic fracturing effects. In addition, the contact between two neighboring layers is regarded as a zone with thickness and mechanical properties, which also influences the effects of hydraulic fracturing in reservoirs. The initial hydraulic fracture was unable to propagate into neighboring layers under a relatively low contact’s Young’s modulus. When associated local tensile stresses exceeded the rock strength, hydraulic fractures propagated into neighboring layers. Moreover, with the contact’s Young’s modulus becoming higher, the fractured area increased rapidly first, then slowly and finally became stable.  相似文献   

11.
李传夫  李术才  李树忱 《岩土力学》2006,27(Z1):198-202
常规方法研究非饱和裂隙介质往往采用宏观连续体的概念,然而现场试验证明,非饱和裂隙岩体中的渗流具有相当的非均质性。天然裂隙岩体的非饱和渗流是发生在三维裂隙网络中的多场非等温流,要在不完整的信息基础上刻画和表现介质的不均匀性,通常借助随机模拟的手段。将三维裂隙系统近似为概化的二维非均匀多孔介质的二维平面裂隙。采用模拟退火法将取自各种来源的信息资料通过建立适当的目标函数汇集至模型之中,以此模拟裂隙的空间特征。  相似文献   

12.
Numerical modeling of stress-permeability coupling in rough fractures   总被引:2,自引:0,他引:2  
A numerical model is described for coupled flow and mechanical deformation in fractured rock. The mechanical response of rock joints to changes in hydraulic pressure is strongly influenced by the geometric characteristics of the joint surfaces. The concept of this work is to combine straightforward finite element solutions with complex and realistic fracture surface geometry in order to reproduce the non-linear stress-deformation-permeability coupling that is commonly observed in fractures. Building on the numerous studies that have expanded the understanding of the key parameters needed to describe natural rough-walled fractures, new methods have been developed to generate a finite element mesh representing discrete fractures with realistic rough surface geometries embedded in a rock matrix. The finite element code GeoSys/Rockflow was then used to simulate the coupled effects of hydraulic stress, mechanical stress, and surface geometry on the evolving permeability of a single discrete fracture. The modeling concept was experimentally verified against examples from the literature. Modeling results were also compared to a simple interpenetration model.  相似文献   

13.
裂隙岩体渗透系数确定方法综述   总被引:4,自引:2,他引:4  
总结了近年来裂隙岩体渗透系数的确定方法,主要有现场水力试验法、裂隙测量法和离散裂隙网络渗流数值试验法。每种方法都有自己的适用性和测量尺度。裂隙岩体渗透系数存在尺度效应,针对不同尺度的研究对象,应尽量选取与渗流模型网格剖分尺度匹配的测量方法。裂隙的延伸具有方向性,测试点的布设应合理科学,不同测量方法需要结合起来才能得到裂隙岩体真实的渗透系数。  相似文献   

14.
15.
Microcracks in brittle rocks affect not only the local mechanical properties, but also the poroelastic behavior and permeability. A continuum coupled hydro-mechanical modeling approach is presented using a two-scale conceptual model representing realistic rock material containing micro-fractures. This approach combines a microcrack-based continuous damage model within generalized Biot poroelasticity, in which the tensors of macroscopic elastic stiffness, Biot effective stress coefficient and of overall permeability are directly related to microcrack growth. Heterogeneity in both mechanical and hydraulic properties evolves from an initially random distribution of damage to produce localized failure and fluid transmission. A significant advantage of the approach is the ability to accurately predict the evolution of realistic fracturing and associated fluid flow in permeable rocks where pre-existing fractures exert significant control. The model is validated for biaxial failure of rock in compression and replicates typical pre- and post-peak strength metrics of stress drop, AE event counts, permeability evolution and failure modes. The model is applied to the simulation of hydraulic fracturing in permeable rocks to examine the effects of heterogeneities, permeability and borehole pressurization rate on the initiation of fracturing. The results indicate that more homogenous rocks require higher hydraulic pressure to initiate fracturing and breakdown. Moreover, both the fracturing initiation pressure and breakdown pressure decrease with permeability but increase with borehole pressurization rate, and the upper and lower limit of the initiation pressure are seen to be given by the impermeable (Hubbert–Willis) and permeable (Haimson–Fairhurst) borehole wall solutions, respectively. The numerical results are shown to be in good agreement with the experimental observations and theoretical results. This coupled damage and flow modeling approach provides an alternative way to solve a variety of complicated hydro-mechanical problems in practical rock engineering with the process coupling strictly enforced.  相似文献   

16.
In this paper a new analytical model is proposed to determine the permeability tensor for fractured rock masses based on the superposition principle of liquid dissipation energy. This model relies on the geometrical characteristics of rock fractures and the corresponding fracture network, and demonstrates the coupling effect between fluid flow and stress/deformation. This model empirically considers the effect of pre‐peak shear dilation and shear contraction on the hydraulic behavior of rock fractures and can be used to determine the applicability of the continuum approach to hydro‐mechanical coupling analysis. Results of numerical analysis presented in this paper show that the new model can effectively describe the permeability of fractured rock masses, and can be applied to the coupling analysis of seepage and stress fields. Copyright © 2008 John Wiley & Sons, Ltd.  相似文献   

17.
Coupled hydro-mechanical (HM) processes are significant in geological engineering such as oil and gas extraction, geothermal energy, nuclear waste disposal and for the safety assessment of dam foundations and rock slopes, where the geological media usually consist of fractured rock masses. In this study, we developed a model for the analysis of coupled hydro-mechanical processes in porous rock containing dominant fractures, by using the numerical manifold method (NMM). In the current model, the fractures are regarded as different material domains from surrounding rock, i.e., finite-thickness fracture zones as porous media. Compared with the rock matrix, these fractured porous media are characterized with nonlinear behavior of hydraulic and mechanical properties, involving not only direct (poroelastic) coupling but also indirect (property change) coupling. By combining the potential energy associated with mechanical responses, fluid flow and solid–fluid interactions, a new formulation for direct HM coupling in porous media is established. For indirect coupling associated with fracture opening/closure, we developed a new approach implicitly considering the nonlinear properties by directly assembling the corresponding strain energy. Compared with traditional methods with approximation of the nonlinear constitutive equations, this new formulation achieves a more accurate representation of the nonlinear behavior. We implemented the new model for coupled HM analysis in NMM, which has fixed mathematical grid and accurate integration, and developed a new computer code. We tested the code for direct coupling on two classical poroelastic problems with coarse mesh and compared the results with the analytical solutions, achieving excellent agreement, respectively. Finally, we tested for indirect coupling on models with a single dominant fracture and obtained reasonable results. The current poroelastic NNM model with a continuous finite-thickness fracture zone will be further developed considering thin fractures in a discontinuous approach for a comprehensive model for HM analysis in fractured porous rock masses.  相似文献   

18.
岩体内裂隙等非连续结构面对岩体的强度及变形等力学特性有着显著的影响,研究岩体裂隙起裂、扩展、相互作用和贯通机制,对工程岩体力学行为的表征和工程性能的评价十分重要。本文基于连续介质力学模型的离散元方法,通过考虑裂隙分布、模型加载条件及其与裂隙产状的关系,建立了一系列裂隙力学计算模型,研究了不同模型裂隙扩展演化特征和岩体破裂机制,分析了岩体裂隙扩展规律及其对岩体破坏路径和强度的影响,研究结果表明:(1)裂隙岩体模型加载条件下的破坏起裂点、最终贯通破坏特征及损伤分布受控于裂隙的产状及其与最大主压应力取向角度大小及围压大小。(2)裂隙弱面走向与最大主压应力取向斜交时,裂隙弱面在加载条件下其端部裂隙扩展、贯通破坏表现比较明显,反之,当裂隙弱面走向与最大主压应力取向一致时,裂隙弱面被动影响裂隙模型内新生裂隙的萌生、扩展和贯通模式,自身未出现新的扩展破坏。(3)裂隙数目的增多和围压的增大会显著增加模型内部剪切裂缝的数量和模型破坏后的破碎程度,模型内部的损伤区域主要围绕破裂面呈滑移线型交叉分布,非破裂面区域损伤呈条带状X型分布。(4)裂隙弱面走向与最大主压应力取向斜交时,裂隙对岩体模型强度的弱化程度高于裂隙弱面走向与最大主压应力取向一致的情况,而裂隙模型破坏后的残余强度则正好相反。  相似文献   

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

20.
增强型地热系统(EGS)是在干热岩技术基础上提出来的一个清洁能源概念,水力压裂建立人工热储是开采地下干热岩热能的有效方法之一。利用TOUGH2系列软件对增强型地热系统进行模拟,具体介绍了对水压致裂过程中裂隙网络模拟的处理方法。裂隙中的水流可以采用不同的概念模型,最为常见的模型包括双空隙率、双渗透率、多重相互作用连续统一体(MINC)以及有效连续统一体(ECM),这些模型明确了对离散的裂隙和基质的模拟方法。应根据基质的渗透性和裂隙的性质灵活地选择裂隙处理方法,也可将不同方法结合起来使用。提出了几种有效的混合模拟方案,对将来高温岩体地热开发具有重要意义。  相似文献   

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

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