共查询到20条相似文献，搜索用时 15 毫秒
1.
The representation of the surface of a rock fracture and a numerical method to simulate fluid flow in single fractures are the keys to understanding the hydraulic behaviour of rock fractures. In this paper, a cellular automaton (CA) approach is used to generate the single fracture structure, which is assumed to be composed of contacts and voids. We develop a CA evolution rule to produce a contact area, and randomly model a single rock fracture with different contact ratios to reflect natural fracture properties such as dead voids, islands and tortuous flow path. Then, based on the localisation theory of a CA, a numerical method to simulate fluid flow in single fractures with contacts is developed. In this method, the fracture is discretised into a system composed of cell elements. Different apertures, i.e., zero for contacts and non-zero for voids, are assigned to each cell element. Therefore, the contribution of the cell elements in a contact on a cell’s transmissivity can be ignored completely. The local transmissivity is assumed to conform to the cubic law. The fluid flow in a fracture with different contact situations is then modelled using the method established in this paper. The fluid flow path, flow velocity and fluid head distributions as well as the channel flow in the fracture are well-modelled. The flow behaviour of the fracture strongly depends on the effective fluid flow path. 相似文献
2.
岩石破裂行为的实验研究 总被引：6，自引：0，他引：6
岩样的破裂行为、破坏过程和参数测试是裂隙断裂构造研究的基础和依据。实验岩石样品采自四川东部和新疆北部地区,为测试准确起见,对岩样进行了应力等值线的有限元法计算。通过单轴和三轴实验的岩样破坏观察和应力应变曲线对比,将岩石的破裂行为、应力应变划分为四个阶段,即裂隙压密阶段、弹性变形阶段、微观劈裂阶段和宏观破裂阶段。基于单轴抗压实验岩石劈裂—破裂—碎裂发展过程的微观分析,可以看出宏观破裂主要是沿岩样原有的隐裂隙、临界裂隙发育的,许多新裂隙则主要是在宏观破裂阶段产生的。 相似文献
3.
Friction Factor of Water Flow Through Rough Rock Fractures 总被引：2，自引：1，他引：1
Fluid flow through rock joints occurs in many rock engineering applications. As the fluid flows through rough-walled rock fractures, pressure head loss occurs due to friction drag of the wall and local aperture changes. In this study, the friction factor was experimentally investigated by performing flow tests through sandstone fractures with joint roughness coefficient ranging from 5.5 to 15.4 under changing normal stress from 0.5 to 3.5 MPa. According to the experimental results, the friction factor was formulated as a function of two-independent variables—Reynolds number and relative roughness. Relative roughness is defined as the ratio of maximum asperity height to equivalent hydraulic aperture. The experimental results show that the proposed predictor of the friction factor fits the data with a coefficient of determination R 2 > 0.93. Sensitivity analyses indicate that in general, the proposed friction factor increases with the relative roughness of confined fractures. The large difference of friction factor induced by relative roughness occurs when the Reynolds number is lower than unity, especially for Re < 0.2. For Reynolds numbers greater than unity, the difference of friction factor induced by relative roughness is smaller. Inclusion of joint roughness in calculating the friction resistance to fluid flow in rough rock joints and the influence of normal stress to the joints is a major step towards more accurate predictions for fluid flow in underground joint networks. This study provides a significant improvement in fundamental understanding of fluid flow in the jointed strata. 相似文献
4.
An experimental study of single and two-phase fluid flow through fractured granite specimens 总被引：3，自引：2，他引：1
P. G. Ranjith 《Environmental Earth Sciences》2010,59(7):1389-1395
Triaxial tests on the two-phase flow of air and water through fractured granite specimens were performed to discover whether
the two-phase fluid flow within rock fractures was laminar or turbulent. The two-phase flow characterization was carried out
based on the macroscopic two-phase steady state flow model and the homogeneous steady state flow model. Rock specimens with
a single natural fracture (joint roughness coefficient, JRC < 10) were tested using two-phase, high pressure triaxial rig.
Experimental results show that the estimated Reynolds numbers for various inlet fluid pressures are well below 1000. The findings
of this study reveal that both single and two-phase flow through rock fractures (JRC < 10) can be characterized as laminar
flows at moderate inlet fluid pressures. However, for single-phase air flow, an increase in inlet air pressures may result
in the formation of turbulent flow. 相似文献
5.
We derive a macroscopic model for single-phase, incompressible, viscous fluid flow in a porous medium with small cavities called vugs. We model the vuggy medium on the microscopic scale using Stokes equations within the vugular inclusions, Darcy's law within the porous rock, and a Beavers–Joseph–Saffman boundary condition on the interface between the two regions. We assume periodicity of the medium and obtain uniform energy estimates independent of the period. Through a two-scale homogenization limit as the period tends to zero, we obtain a macroscopic Darcy's law governing the medium on larger scales. We also develop some needed generalizations of the two-scale convergence theory needed for our bimodal medium, including a two-scale convergence result on the Darcy–Stokes interface. The macroscopic Darcy permeability is computable from the solution of a cell problem. An analytic solution to this problem in a simple geometry suggests that: (1) flow along vug channels is primarily Poiseuille with a small perturbation related to the Beavers–Joseph slip, and (2) flow that alternates from vug to matrix behaves as if the vugs have infinite permeability. 相似文献
6.
This paper aims at developing a pore-scale numerical model for non-Darcy fluid flow through rough-walled fractures. A simple general relationship between the local hydraulic conductivity and the flow velocity is proposed. A new governing equation for non-Darcy fluid flow through rough-walled fractures is then derived by introducing this relationship into the Reynolds equation. Based on the non-linear finite element method, a self-developed code is used to simulate the non-Darcy fluid flow through fractures. It is found that the macroscopic results obtained by the numerical simulation agree well with the experimental results. Furthermore, some interesting experimental observations can be reproduced. 相似文献
7.
A calibration procedure, in which coupled effects of microscopic parameters are considered, is proposed to determine the values of the microscopic parameters in the Discrete Element Method (DEM) for Fujian quartz sand. Laboratory tri-axial tests are conducted to be compared with the DEM simulations and the effects of end restraint in the laboratory tests are eliminated through a digital image measurement system. Sensitivities of the macroscopic behaviour of the specimen to the microscopic parameters are analyzed through DEM simulations. Four coupled effects of the microscopic parameters on the macroscopic behaviour are investigated through a graphic method and then considered in the calibration procedure. 相似文献
8.
Cracks play a major role in most rocks submitted to crustal conditions. Mechanically, cracks make the rock much more compliant. They also make it much easier for fluid to flow through any rock body. Relying on Fracture Mechanics and Statistical Physics, we introduce a few key concepts, which allow to understand and quantify how cracks do modify both the elastic and transport properties of rocks. The main different schemes, which can be used to derive the elastic effective moduli of a rock, are presented. It is shown from experimental results that an excellent approximation is the so-called non-interactive scheme. The main consequences of the existence of cracks on the elastic waves is the development of elastic anisotropy (due to the anisotropic distribution of crack orientations) and the dispersion effect (due to microscopic local fluid flow). At a larger scale, macroscopic fluid flow takes place through the crack network above the percolation threshold. Two macroscopic fluid flow regimes can be distinguished: the percolative regime close to the percolation threshold and the connected regime well above it. Experimental data on very different rock types show both of these behaviors. 相似文献
9.
基于微观力学特性的脆性岩石变形过程模拟 总被引：1，自引：0，他引：1
为了建立能够描述岩石宏观变形及其组成部分的变形模拟方法,考虑微缺陷的不均匀性对岩石宏观变形产生的影响,将岩石视为由空隙部分材料和骨架部分材料两部分组成。首先,基于岩石各组成部分材料变形分析,利用真应变描述方法及非线性变形的特点建立空隙部分材料的变形力学分析方法,同时,引入统计损伤以及固体力学理论建立骨架部分材料的变形力学分析方法。然后,通过岩石及其组成部分的变形力学分析建立岩石宏观变形力学分析方法,从而获得能够模拟岩石变形破坏全过程的统计损伤本构模型,并给出了模型参数的确定方法。最后,将文中模型与既有模型的理论曲线与试验曲线进行了对比分析,且基于文中模型对岩石及其组成部分变形进行了讨论。研究结果表明：文中模型不仅能反映既有模型所描述的主要变形特征,还能克服其难以较好地反映初期宏观变形非线性的不足;同时,该模型还能描述岩石组成部分的变形过程,不仅揭示了岩石宏观变形与其组成部分变形之间的关系,也揭示了空隙部分材料变形是引起岩石宏观变形非线性的根本原因。表明本文模型与方法具有一定的合理性与优越性。 相似文献
10.
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.
相似文献11.
对交叉裂隙渗流传质特性的定量描述是研究整个裂隙网络渗透传质特性的基础。为真实模拟水流及溶质在三维交叉裂隙中的运移过程,首先通过三维轮廓仪获取天然岩石裂隙表面的形貌数据,再应用三维重构技术生成相应的三维交叉裂隙模型,随后求解Navier-Stokes方程,假定溶质运移满足Fick定律,模拟水流和溶质在三维交叉裂隙中的运移过程。通过对比粗糙裂隙模型与平行平板模型的模拟结果发现：粗糙度对流体的分布及流动状态存在显著的影响;不同进、出口工况下的流体流动及溶质运移状态亦表明：裂隙交叉的几何形貌会显著地影响溶质混合行为。这些结果表明,目前被广泛采用的平行平板模型在评估岩体内特别是交叉口的物质运移特性时将导致较大的偏差,在将来的研究中有必要针对裂隙交叉口的几何特征建立修正的模型以提高评估的准确性。 相似文献
12.
An embedded fracture modeling framework for simulation of hydraulic fracturing and shear stimulation 总被引：1，自引：0，他引：1
Jack H. Norbeck Mark W. McClure Jonathan W. Lo Roland N. Horne 《Computational Geosciences》2016,20(1):1-18
A numerical modeling framework is described that is able to calculate the coupled processes of fluid flow, geomechanics, and rock failure for application to general engineering problems related to reservoir stimulation, including hydraulic fracturing and shear stimulation. The numerical formulation employs the use of an embedded fracture modeling approach, which provides several advantages over more traditional methods in terms of computational complexity and efficiency. Specifically, the embedded fracture modeling strategy avoids the usual requirement that the discretization of the fracture domain conforms to the discretization of the rock volume surrounding the fractures. As fluid is exchanged between the two domains, conservation of mass is guaranteed through a coupling term that appears as a simple source term in the governing mass balance equations. In this manner, as new tensile fractures nucleate and propagate subject to mechanical effects, numerical complexities associated with the introduction of new fracture control volumes are largely negated. In addition, the ability to discretize the fractures and surrounding rock volume independently provides the freedom to choose an acceptable level of discretization for each domain separately. Three numerical examples were performed to demonstrate the utility of the embedded fracture model for application to problems involving fluid flow, mechanical deformation, and rock failure. The results of the numerical examples confirm that the embedded fracture model was able to capture accurately the complex and nonlinear evolution of reservoir permeability as new fractures propagate through the reservoir and as fractures fail in shear. 相似文献
13.
This paper is aimed towards investigating the filtration law of an incompressible viscous Newtonian fluid through a rigid non-inertial porous medium (e.g. a porous medium placed in a centrifuge basket). The filtration law is obtained by upscaling the flow equations at the pore scale. The upscaling technique is the homogenization method of multiple scale expansions which rigorously gives the macroscopic behaviour and the effective properties without any prerequisite on the form of the macroscopic equations. The derived filtration law is similar to Darcy's law, but the tensor of permeability presents the following remarkable properties: it depends upon the angular velocity of the porous matrix, it verifies Hall–Onsager's relationship and it is a non-symmetric tensor. We thus deduce that, under rotation, an isotropic porous medium leads to a non-isotropic effective permeability. In this paper, we present the results of numerical simulations of the flow through rotating porous media. This allows us to highlight the deviations of the flow due to Coriolis effects at both the microscopic scale (i.e. the pore scale), and the macroscopic scale (i.e. the sample scale). The above results confirm that for an isotropic medium, phenomenological laws already proposed in the literature fails at reproducing three-dimensional Coriolis effects in all types of pores geometry. We show that Coriolis effects may lead to significant variations of the permeability measured during centrifuge tests when the inverse Ekman number Ek−1 is 𝒪(1). These variations are estimated to be less than 5% if Ek−1<0.2, which is the case of classical geotechnical centrifuge tests. We finally conclude by showing that available experimental data from tests carried out in centrifuges are not sufficient to determining the effective tensor of permeability of rotating porous media. Copyright © 2004 John Wiley & Sons, Ltd. 相似文献
14.
Mining induced subsidence can significantly affect mining costs where major surface facilities and natural environment need to be protected. Overburden grout injection is a technology used to control coal mine subsidence by injecting the mine waste material extracted from the coal back into the inter-burden rock during longwall mining. The flowing slurry is here categorised as a nonlinear viscous cohesive (Bingham plastic) fluid. During longwall mining the grout slurry is pumped into the separated beds of the rock mass through a central vertical borehole, which is drilled deep into the inter-burden rock strata above the coal seam. However, a blockage can occur in the injection system when the slurry velocity falls below a certain critical threshold velocity, indicating a material phase change from cohesive-viscous to cohesive-frictional. In situ field injection tests through boreholes have been simulated at a smaller scale at the CSIRO laboratory in Brisbane by pumping the slurry through a radial disk (gap = 4 mm) from its centre. Laboratory experiments indicate a general, nonlinear, cohesive, viscous, frictional model for shear behaviour of the slurry, in which the material shear parameters are functions of the disk radial distance. Complete dimensional and dimensionless analytical solutions have been developed based on an approach related to Bingham–Herschel–Bulkley fluid mechanics. The derived formulae include relations for minimum pump pressure, local pressure and pressure gradient, wall shear stress, volume rate, velocity and velocity gradient. The theoretical results match the experimental measurements. The experiments covered slurries with maximum particle sizes of 0.5 to 2 mm with about 50% being larger than 100 µm. The viscosities at the various solids concentrations were measured with a standard torsion viscometer. This study differs from the previous research in several distinct aspects, namely, consideration of the variable shear parameters rather than fixed values, inclusion of total nonlinear behaviour, and implementation of a friction function to mimic behaviour of the deposited and consolidating stiff slurry, which can cause a significant pressure rise as a result of the increased shear resistance. 相似文献
15.
Fracture flow modelling based on satellite images of the Wajid Sandstone, Saudi Arabia 总被引：1，自引：0，他引：1
Conny Zeeb Daniel Göckus Paul Bons Hussain Al Ajmi Randolf Rausch Philipp Blum 《Hydrogeology Journal》2010,18(7):1699-1712
Large-scale geological features have been identified by satellite imagery and global positioning system data in the Wajid Sandstone in Saudi Arabia. The main objective is to evaluate the importance of fractures for the overall flow behaviour in this fractured rock aquifer and to estimate in-situ hydraulic apertures. Data on fractures and lineaments were available for three outcrops. By applying a “cut-out” routine on the fracture endpoint data of these fracture trace windows, three deterministic discrete fracture networks (DFN), with an area of 100 m?×?100 m, could be generated. These were used to simulate the fracture flow and to determine the hydraulic conductivity tensors. Using additional data on hydraulic pumping tests and matrix conductivities, in-situ hydraulic apertures could be determined. Average in-situ hydraulic apertures range from 1,300 to 1,700 µm. Observations from the field support these results. In addition, a hydraulic conductivity ratio between the matrix and fracture system was used to identify the contribution of the DFN to the overall fluid transport. A ratio of 10.4 was determined, which indicates that the effective flow behaviour in the Wajid Sandstone aquifer is not entirely dominated by the fracture system, though evidently strongly controlled by it. 相似文献
16.
Size and stress gradient effects on fracture around cavities 总被引：1，自引：2，他引：1
B. J. Carter 《Rock Mechanics and Rock Engineering》1992,25(3):167-186
Summary The effect of hole size on fracture initiation around cavities in rock is examined through physical model tests on Tyndall limestone. Twenty five models, employing nine different hole sizes, ranging from 3.2 to 62 mm in diameter were tested in uniaxial compression. Strain gauges placed both on the face of the model block and inside the cylindrical cavity were used to record the state of fracture.Fracture initiation data for primary, remote and sidewall fractures were determined using the stress-strain curves from the respective strain gauge locations. For all three fracture types, the initiation stress decreases with increasing hole size and appears to approach a horizontal asymptote for large cavity sizes. The fracture mechanics formulations and the stress averaging method give the best fit to the primary fracture initiation except for very small cavity sizes. The stress averaging method also provides a very good fit to the sidewall fracture initiation data. 相似文献
17.
裂隙岩体流固耦合问题是目前国内外研究热点之一,采用离散元软件UDEC对裂隙岩体发生节理剪胀的渗透性变化规律进行了模拟分析。基于现场调查的裂隙信息统计生成裂隙网络岩体模型。 通过固定垂直应力、不断增加应力比RS(RS=水平应力/垂直应力)使岩体出现剪胀,采用库伦滑移节理模式对岩体在剪胀过程中的渗透性变化情况进行模拟。结果发现:当应力比较小(RS3.1)时,节理水力隙宽、流速、渗透系数等参数都随着应力比的增加表现出明显的降低; 而当岩体出现剪胀现象之后(应力比大于3.1),发生剪切滑移和剪胀现象的节理控制着裂隙岩体的总体渗流行为,与不考虑节理剪胀的计算结果相比,岩体渗透能力出现了显著增长。这一结果表明,剪胀对裂隙岩体渗透性的影响是显著而不可忽视的。 相似文献
18.
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. 相似文献
19.
The fluid flow in rock fractures during shear processes has been an important issue in rock mechanics and is investigated in this paper using finite element method (FEM), considering evolutions of aperture and transmissivity with shear displacement histories under different normal stress and normal stiffness conditions as measured during laboratory coupled shear-flow tests. The distributions of fracture aperture and its evolution during shearing were calculated from the initial aperture, based on the laser-scanned sample surface roughness results, and shear dilations measured in the laboratory tests. Three normal loading conditions were adopted in the tests: simple normal stress and mixed normal stress and normal stiffness to reflect more realistic in situ conditions. A special algorithm for treatment of the contact areas as zero-aperture elements was used to produce more accurate flow field simulations, which is important for continued simulations of particle transport but often not properly treated in literature. The simulation results agree well with the measured hydraulic apertures and flow rate data obtained from the laboratory tests, showing that complex histories of fracture aperture and tortuous flow fields with changing normal loading conditions and increasing shear displacements. With the new algorithm for contact areas, the tortuous flow fields and channeling effects under normal stress/stiffness conditions during shearing were more realistically captured, which is not possible if traditional techniques by assuming very small aperture values for the contact areas were used. These findings have an important impact on the interpretation of the results of coupled hydro-mechanical experiments of rock fractures, and on more realistic simulations of particle transport processes in fractured rocks. 相似文献
20.
岩体裂隙网络渗流广泛存在于地下工程中,对地下工程的建设和运行安全具有重要的影响。因此,研究裂隙网络渗流有着重要的理论和实际意义。本文根据立方定律和Forchheimer方程推导所得的交叉裂隙渗流模型,运用数值模拟和人工神经网络方法,对平面交叉裂隙渗流模型非线性参数与模型几何条件的关系进行探究。通过数值模拟,获得了平面交叉裂隙非线性渗流模型的参数;运用人工神经网络遗传算法,探究了交叉裂隙几何条件与交叉裂隙渗流模型中非线性系数之间的关系,证明了平面交叉裂隙非线性渗流模型适用于描述交叉裂隙渗流规律,验证了神经网络方法预测非线性系数的可行性和准确性。同时,还对比分析了运用拟合数值表达式和人工神经网络两种方法的特点。 相似文献