共查询到20条相似文献，搜索用时 937 毫秒
1.
Contact mechanism of a rock fracture subjected to normal loading and its impact on fast closure behavior during initial stage of fluid flow experiment 
下载免费PDF全文
下载免费PDF全文 Fast closure of rock fractures has been commonly observed in the initial stage of fluid flow experiments at environmental temperatures under low or moderate normal stresses. To fully understand the mechanisms that drive this fast closure, the evolution of local stresses acting on contacting asperities on the fracture surfaces prior to fluid flow tests needs to be evaluated. In this study, we modeled numerically the asperity deformation and failure processes during initial normal loading, by adopting both elastic and elastic–plastic deformation models for the asperities on a real rock fracture with measured surface topography data, and estimated their impact on initial conditions for fluid flow test performed under laboratory conditions. Compared with the previous models that simulate the normal contact of a fracture as the approach of two rigid surfaces without deformations, our models of deformable asperities yielded smaller contact areas and higher stresses on contacting asperities at a given normal stress or normal displacement. The results show that the calculated local stresses were concentrated on the contacts of a few major asperities, resulting in crushing of asperity tips. With these higher contact stresses, however, the predicted closure rates by pressure solution are still several orders of magnitude lower than that of the experimental measurements at the initial stage of fluid flow test. This indicates that single pressure solution may not likely to be the principal compaction mechanism for this fast closure, and that the damages on contacting asperities that occur during the initial normal loading stage may play an important role. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献
2.
Lifang Zou Boris G. Tarasov Arcady V. Dyskin Deepak P. Adhikary Elena Pasternak Weiya Xu 《Rock Mechanics and Rock Engineering》2013,46(1):67-81
This paper presents the results of laboratory experiments conducted to study the impact of stress on fracture deformation and permeability of fractured rocks. The physical models (laboratory specimens) consisted of steel cubes simulating a rock mass containing three sets of orthogonal fractures. The laboratory specimens were subjected to two or three cycles of hydrostatic loading/unloading followed by the measurement of displacement and permeability. The results show a considerable difference in both deformation and permeability trends between the first loading and the subsequent loading/unloading cycles. However, the micrographs of the contact surfaces taken before and after the tests show that the standard deviation of asperity heights of measured surfaces are affected very little by the loadings. This implies that both deformation and permeability are rather controlled by the highest surface asperities which cannot be picked up by the conventional roughness characterization technique. We found that the dependence of flow rate on mechanical aperture follows a power law with the exponent n smaller or larger than three depending upon the loading stage. Initially, when the maximum height of the asperities is high, the exponent is slightly smaller than 3. The first loading, however, flattens these asperities. After that, the third loading and unloading yielded the exponent of around 4. Due to the roughness of contact surfaces, the flow route is no longer straight but tortuous resulting in flow length increase. 相似文献
3.
Hide Yasuhara Naoki Kinoshita Hiroaki Ohfuji Dae Sung Lee Shinichiro Nakashima Kiyoshi Kishida 《Applied Geochemistry》2011,26(12):2074-2088
Examining the evolution of fracture permeability under stressed and temperature-elevated conditions, a series of flow-through experiments on a single rock fracture in granite has been conducted under confining pressures of 5 and 10 MPa, under differential water pressures ranging from 0.04 to 0.5 MPa, and at temperatures of 20–90 °C, for several hundred hours in each experiment. Measurements of fluid and dissolved mass fluxes, and post-experimental microscopy, were conducted to constrain the progress of mineral dissolution and/or precipitation and to examine its effect on transport properties. Generally, the fracture aperture monotonically decreased with time at room temperature, and reached a steady state in relatively short periods (i.e., <400 h). However, once the temperature was elevated to 90 °C, the aperture resumed decreasing and kept decreasing throughout the rest of the experimental periods. This reduction may result from the removal of the mineral mass from the bridging asperities within the fracture. Post-experimental observations by scanning electron microscopy, coupled with energy dispersive X-ray spectroscopy (SEM-EDX), revealed the formation of several kinds of secondary minerals such as silica and calcite. However, the precipitated minerals seemed to have had little influence on the flow characteristics within the fracture, because the precipitation was limited to quite local and small areas. The evolving rates and ultimate magnitudes of the fracture aperture are likely to be controlled by the stress exerted over the contacting asperities and temperatures, and by the prescribed flow conditions. Thus, this complex behavior should be attributed to the coupled chemically- and mechanically-induced effect. A coupled chemo–mechano conceptual model, accounting for pressure and free-face dissolutions, is presented in this paper to follow the evolution of the fracture permeability observed in the flow-through experiments. This model addresses the two dissolution processes at the contacting asperities and the free walls within the fractures, and is also capable of describing multi-mineral dissolution behavior. The model shows that the evolution of a fracture aperture (or related permeability) and of element concentrations may be followed with time under arbitrary temperature and pressure conditions. The model predictions for the evolving fracture aperture and elements concentrations show a relatively good agreement with the experimental measurements, although it is not possible to replicate the abrupt reduction observed in the early periods of the experiments, which is likely to be due to an unaccounted mechanism of more stress-mediated fracture compaction driven by the fracturing of the propping asperities. 相似文献
4.
Natural fractures are characterized by rough surfaces and complex fluid flows. A large distribution of apertures (residual voids) within their walls and the presence of contact points (in situ normal loads) produce heterogeneous flows (channeling). The resulting permeabilities, porosities or fluid–rock exchange surfaces cannot be realistically modeled by parallel and smooth plate models. Four natural fractures are sampled at different depths and degrees of alteration in the Soultz sandstone and granite (EPS1 drillhole, Soultz-sous-Forêts, Bas-Rhin, France). The fracture surfaces are measured with mechanical profilometry and maps of asperity heights (XYZ). Resulting local apertures (XYe) are then calculated. A statistical study of the surface profiles (XZ) show that the fractures are more or less rough and tortuous according to the types of alteration. Altered samples are characterized by smoother surfaces of fractures. Such differences imply that (i) the average fracture aperture is not representative for the whole fracture and that (ii) the different local apertures should be integrated in hydraulic and mechanical models. A hydraulic model (finite difference calculations) of fluid flow, taking into account the elastic closure (Hertz contact theory) of fractures with depth, is used. Maps of contact points and relative local loads within the fracture planes are compared to flow maps. They show different channeling of fluid flows. Strongly altered fractures are characterized by homogeneous fluxes despite the presence of numerous contact zones during the closure of fracture. By contrast, fresh fractures develop, increasing fluid flow channels with depth.Fracture closure (increasing normal stress) does not systematically increase the channeling of fluid flow. There is evidence for a general smoothing out of the irregularities of the fracture walls due to precipitation of secondary minerals, indicating that the cubic law can be commonly valid, also at great crustal depth but this validity depends on the degree of fracture alteration. Mineralogical and geochemical observations, thus, should be taken into account to perform more accurate permeability calculations and models of fluid circulation in fracture networks. 相似文献
5.
This paper presents the development of a discrete fracture model of fully coupled compressible fluid flow, adsorption and geomechanics to investigate the dynamic behaviour of fractures in coal. The model is applied in the study of geological carbon dioxide sequestration and differs from the dual porosity model developed in our previous work, with fractures now represented explicitly using lower-dimensional interface elements. The model consists of the fracture-matrix fluid transport model, the matrix deformation model and the stress-strain model for fracture deformation. A sequential implicit numerical method based on Galerkin finite element is employed to numerically solve the coupled governing equations, and verification is completed using published solutions as benchmarks. To explore the dynamic behaviour of fractures for understanding the process of carbon sequestration in coal, the model is used to investigate the effects of gas injection pressure and composition, adsorption and matrix permeability on the dynamic behaviour of fractures. The numerical results indicate that injecting nonadsorbing gas causes a monotonic increase in fracture aperture; however, the evolution of fracture aperture due to gas adsorption is complex due to the swelling-induced transition from local swelling to macro swelling. The change of fracture aperture is mainly controlled by the normal stress acting on the fracture surface. The fracture aperture initially increases for smaller matrix permeability and then declines after reaching a maximum value. When the local swelling becomes global, fracture aperture starts to rebound. However, when the matrix permeability is larger, the fracture aperture decreases before recovering to a higher value and remaining constant. Gas mixtures containing more carbon dioxide lead to larger closure of fracture aperture compared with those containing more nitrogen. 相似文献
6.
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. 相似文献
7.
A numerical model is presented to describe the evolution of fracture aperture (and related permeability) mediated by the competing chemical processes of pressure solution and free‐face dissolution/precipitation; pressure (dis)solution and precipitation effect net‐reduction in aperture and free‐face dissolution effects net‐increase. These processes are incorporated to examine coupled thermo‐hydro‐mechano‐chemo responses during a flow‐through experiment, and applied to reckon the effect of forced fluid injection within rock fractures at geothermal and petroleum sites. The model accommodates advection‐dominant transport systems by employing the Lagrangian–Eulerian method. This enables changes in aperture and solute concentration within a fracture to be followed with time for arbitrary driving effective stresses, fluid and rock temperatures, and fluid flow rates. This allows a systematic evaluation of evolving linked mechanical and chemical processes. Changes in fracture aperture and solute concentration tracked within a well‐constrained flow‐through test completed on a natural fracture in novaculite (Earth Planet. Sci. Lett. 2006, in press) are compared with the distributed parameter model. These results show relatively good agreement, excepting an enigmatic abrupt reduction in fracture aperture in the early experimental period, suggesting that other mechanisms such as mechanical creep and clogging induced by unanticipated local precipitation need to be quantified and incorporated. The model is applied to examine the evolution in fracture permeability for different inlet conditions, including localized (rather than distributed) injection. Predictions show the evolution of preferential flow paths driven by dissolution, and also define the sense of permeability evolution at field scale. Copyright © 2006 John Wiley & Sons, Ltd. 相似文献
8.
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. 相似文献
9.
Hydraulic Characteristics of Rough Fractures in Linear Flow under Normal and Shear Load 总被引：19，自引：3，他引：19
Summary
A hydro-mechanical testing system, which is capable of measuring both the flow rates and the normal and shear displacement
of a rock fracture, was built to investigate the hydraulic behaviour of rough tension fractures. Laboratory hydraulic tests
in linear flow were conducted on rough rock fractures, artificially created using a splitter under various normal and shear
loading. Prior to the tests, aperture distributions were determined by measuring the topography of upper and lower fracture
surfaces using a laser profilometer. Experimental variograms of the initial aperture distributions were classified into four
groups of geostatistical model, though the overall experimental variograms could be well fitted to the exponential model.
The permeability of the rough rock fractures decayed exponentially with respect to the normal stress increase up to 5 MPa.
Hydraulic behaviours during monotonic shear loading were significantly affected by the dilation occurring until the shear
stress reached the peak strength. With the further dilation, the permeability of the rough fracture specimens increased more.
However, beyond shear displacement of about 7 to 8 mm, permeability gradually reached a maximum threshold value. The combined
effects of both asperity degradation and gouge production, which prohibited the subsequent enlargement of mean fracture aperture,
mainly caused this phenomenon. Permeability changes during cyclic shear loading showed somewhat irregular variations, especially
after the first shear loading cycle, due to the complex interaction from asperity degradations and production of gouge materials.
The relation between hydraulic and mechanical apertures was analyzed to investigate the valid range of mechanical apertures
to be applied to the cubic law.
Received June 12, 2001; accepted February 26, 2002 Published online September 2, 2002 相似文献
10.
The possibility that gradients in concentration may develop within single pores and fractures, potentially giving rise to scale-dependent mineral dissolution rates, was investigated with experimentally validated reactive transport modeling. Three important subsurface mineral phases that dissolve at widely different rates, calcite, plagioclase, and iron hydroxide, were considered. Two models for analyzing mineral dissolution kinetics within a single pore were developed: (1) a Poiseuille Flow model that applies laboratory-measured dissolution kinetics at the pore or fracture wall and couples this to a rigorous treatment of both advective and diffusive transport within the pore, and (2) a Well-Mixed Reactor model that assumes complete mixing within the pore, while maintaining the same reactive surface area, average flow rate, geometry, and multicomponent chemistry as the Poiseuille Flow model. For the case of a single fracture, a 1D Plug Flow Reactor model was also considered to quantify the effects of longitudinal versus transverse mixing. Excellent agreement was obtained between results from the Poiseuille Flow model and microfluidic laboratory experiments in which pH 4 and 5 solutions were flowed through a single 500 μm diameter by 4000 μm long cylindrical pore in calcite. The numerical modeling and time scale analysis indicated that rate discrepancies arise primarily where concentration gradients develop under two necessary conditions: (1) comparable rates of reaction and advective transport, and (2) incomplete mixing via molecular diffusion. For plagioclase and iron hydroxide, the scaling effects are negligible at the single pore and fracture scale because of their slow rates. In the case of calcite, where dissolution rates are rapid, scaling effects can develop at high flow rates from 0.1 to 1000 cm/s and for fracture lengths less than 1 cm. Under more normal flow conditions where flow is usually slower than 0.001 cm/s, however, mixing via molecular diffusion is effective in homogenizing the concentration field, thus eliminating any discrepancies between the Poiseuille Flow and the Well-Mixed Reactor model. The analysis suggests that concentration gradients are unlikely to develop within single pores and fractures under typical geological/hydrologic conditions, implying that the discrepancy between laboratory and field rates must be attributed to other factors. 相似文献
11.
Summary The mechanical and hydromechanical behaviour of isolated rock joints is of prime importance for a correct understanding of
the behaviour of jointed rock masses. This paper focuses on the mechanical behaviour of a fracture under normal stress (fracture
closure), using approaches based on both experimentation and modelled analysis. Experimental closure tests were carried out
by positioning four displacement transducers around a fracture, leading to results which tended to vary as a function of transducer
location. Such variations can be explained by the non-constant void space distribution between both walls of the fracture.
The present study focuses on the importance of transducer location in such a test, and on the significant role played, in
terms of mechanical response, by the morphology of the fracture surfaces.
An analytical mechanical model is then developed, which takes into account the deformation of surface asperities and of the
bulk material surrounding the fracture; it also includes the effects of mechanical interaction between contact points. The
model is validated by simulating the behaviour which is very similar to experimental observations. Various parametric studies
(scale effect, spatial distribution of contact points) are then carried out. The study of scale effects reveals a decrease
in the normal stiffness with increasing fracture size. Finally, analysis of the role of various mechanical parameters has
shown that the most influential of these is Young’s modulus corresponding to the bulk material surrounding the joint. Many
applications, such as geothermal fluid recovery from fractures, could benefit from these results.
Correspondence: Antoine Marache, Université Bordeaux 1, GHYMAC, Av. des Facultés, 33405 Talence Cedex, France 相似文献
12.
Fractured-rock aquifers display spatially and temporally variable hydraulic conductivity generally attributed to variable fracture intensity and connectivity. Empirical evidence suggests fracture aperture and hydraulic conductivity are sensitive to in situ stress. This study investigates the sensitivity of fractured-rock hydraulic conductivity, groundwater flow paths, and advection-dominated transport to variable shear and normal fracture stiffness magnitudes for a range of deviatoric stress states. Fracture aperture and hydraulic conductivity are solved for analytically using empirical hydromechanical coupling equations; groundwater flow paths and ages are then solved for numerically using groundwater flow and advection-dispersion equations in a traditional Toth basin. Results suggest hydraulic conductivity alteration is dominated by fracture normal closure, resulting in decreasing hydraulic conductivity and increasing groundwater age with depth, and decreased depth of long flow paths with decreasing normal stiffness. Shear dilation has minimal effect on hydraulic conductivity alteration for stress states investigated here. Results are interpreted to suggest that fracture normal stiffness influences hydraulic conductivity of hydraulically active fractures and, thus, affects flow and transport in shallow (<1 km) fractured-rock aquifers. It is suggested that observed depth-dependent hydraulic conductivity trends in fractured-rock aquifers throughout the world may be partly a manifestation of hydromechanical phenomena. 相似文献
13.
The response of deformable fractures to changes in fluid pressure controls phenomena ranging from the flow of fluids near wells to the propagation of hydraulic fractures. We developed an analysis designed to simulate fluid flows in the vicinity of asperity‐supported fractures at rest, or fully open fractures that might be propagating. Transitions between at‐rest and propagating fractures can also be simulated. This is accomplished by defining contact aperture as the aperture when asperities on a closing fracture first make contact. Locations on a fracture where the aperture is less than the contact aperture are loaded by both fluid pressure and effective stress, whereas locations where the aperture exceeds the contact aperture are loaded only by fluid pressure. Fluid pressure and effective stress on the fracture are determined as functions of time by solving equations of continuity in the fracture and matrix, and by matching the global displacements of the fracture walls to the local deformation of asperities. The resulting analysis is implemented in a numerical code that can simulate well tests or hydraulic fracturing operations. Aperture changes during hydraulic well tests can be measured in the field, and the results predicted using this analysis are similar to field observations. The hydraulic fracturing process can be simulated from the inflation of a pre‐existing crack, to the propagation of a fracture, and the closure of the fracture to rest on asperities or proppant. Two‐dimensional, multi‐phase fluid flow in the matrix is included to provide details that are obscured by simplifications of the leakoff process (Carter‐type assumptions) used in many hydraulic fracture models. Execution times are relatively short, so it is practical to implement this code with parameter estimation algorithms to facilitate interpretation of field data. Copyright © 2006 John Wiley & Sons, Ltd. 相似文献
14.
Eren Ucar Eirik Keilegavlen Inga Berre Jan Martin Nordbotten 《Computational Geosciences》2018,22(4):993-1007
Simulating the deformation of fractured media requires the coupling of different models for the deformation of fractures and the formation surrounding them. We consider a cell-centered finite-volume approach, termed the multi-point stress approximation (MPSA) method, which is developed in order to discretize coupled flow and mechanical deformation in the subsurface. Within the MPSA framework, we consider fractures as co-dimension one inclusions in the domain, with the fracture surfaces represented as line pairs in 2D (face pairs in 3D) that displace relative to each other. Fracture deformation is coupled to that of the surrounding domain through internal boundary conditions. This approach is natural within the finite-volume framework, where tractions are defined on surfaces of the grid. The MPSA method is capable of modeling deformation, considering open and closed fractures with complex and nonlinear relationships governing the displacements and tractions at the fracture surfaces. We validate our proposed approach using both problems, for which analytical solutions are available, and more complex benchmark problems, including comparison with a finite-element discretization. 相似文献
15.
Veena U. Joshi 《Journal of the Geological Society of India》2009,74(6):711-722
Quartz sand grains obtained from a deeply gullied topography along the banks of two tributaries of River Pravara in Godavari
Basin, Maharashtra have been examined with a scanning electron microscope (SEM) in order to make the environmental interpretations
of these deposits. The sediments reveal features resulting from mechanical grinding as well as from chemical alteration. Conchoidal
fractures, cleavage planes, grooves, v-shaped indentations etc are the mechanical features documented on the grains whereas
solution pits of varying sizes and intensity, precipitation surfaces, oriented v-pits, solution crevasses and etching are
the features of diagenetic origin. Few sand grains show the evidence of wind transported sediments. Several evidence indicate
that the samples have undergone digenetic changes. Few grains exhibit the features of intense chemical breakdown. The overall
assemblages of the grain surface features suggest that the samples have been subjected to subaqueous transport for a considerable
period of time. The chemical features such as etching, solution pits or semi circular arcuate steps that are found in abundance
in these grains are due to the dissolution of the sediments in a low energy fluviatile environment, such as in floodplain
region. 相似文献
16.
Numerical Simulation of 3D Hydraulic Fracturing Based on an Improved Flow-Stress-Damage Model and a Parallel FEM Technique 总被引：4，自引：4，他引：0
L. C. Li C. A. Tang G. Li S. Y. Wang Z. Z. Liang Y. B. Zhang 《Rock Mechanics and Rock Engineering》2012,45(5):801-818
The failure mechanism of hydraulic fractures in heterogeneous geological materials is an important topic in mining and petroleum engineering. A three-dimensional (3D) finite element model that considers the coupled effects of seepage, damage, and the stress field is introduced. This model is based on a previously developed two-dimensional (2D) version of the model (RFPA2D-Rock Failure Process Analysis). The RFPA3D-Parallel model is developed using a parallel finite element method with a message-passing interface library. The constitutive law of this model considers strength and stiffness degradation, stress-dependent permeability for the pre-peak stage, and deformation-dependent permeability for the post-peak stage. Using this model, 3D modelling of progressive failure and associated fluid flow in rock are conducted and used to investigate the hydro-mechanical response of rock samples at laboratory scale. The responses investigated are the axial stress–axial strain together with permeability evolution and fracture patterns at various stages of loading. Then, the hydraulic fracturing process inside a rock specimen is numerically simulated. Three coupled processes are considered: (1) mechanical deformation of the solid medium induced by the fluid pressure acting on the fracture surfaces and the rock skeleton, (2) fluid flow within the fracture, and (3) propagation of the fracture. The numerically simulated results show that the fractures from a vertical wellbore propagate in the maximum principal stress direction without branching, turning, and twisting in the case of a large difference in the magnitude of the far-field stresses. Otherwise, the fracture initiates in a non-preferred direction and plane then turns and twists during propagation to become aligned with the preferred direction and plane. This pattern of fracturing is common when the rock formation contains multiple layers with different material properties. In addition, local heterogeneity of the rock matrix and macro-scale stress fluctuations due to the variability of material properties can cause the branching, turning, and twisting of fractures. 相似文献
17.
Groundwater flow is an important control on subsurface evaporite (salt) dissolution. Salt dissolution can drive faulting and associated subsidence on the land surface and increase salinity in groundwater. This study aims to understand the groundwater flow system of Gypsum Canyon watershed in the Paradox Basin, Utah, USA, and whether or not groundwater-driven dissolution affects surface deformation. The work characterizes the groundwater flow and solute transport systems of the watershed using a three-dimensional (3D) finite element flow and transport model, SUTRA. Spring samples were analyzed for stable isotopes of water and total dissolved solids. Spring water and hydraulic conductivity data provide constraints for model parameters. Model results indicate that regional groundwater flow is to the northwest towards the Colorado River, and shallow flow systems are influenced by topography. The low permeability obtained from laboratory tests is inconsistent with field observed discharges, supporting the notion that fracture permeability plays a significant role in controlling groundwater flow. Model output implies that groundwater-driven dissolution is small on average, and cannot account for volume changes in the evaporite deposits that could cause surface deformation, but it is speculated that dissolution may be highly localized and/or weaken evaporite deposits, and could lead to surface deformation over time. 相似文献
18.
This paper presents a composite element algorithm of coupled normal stress and fluid flow process for fractured rock mass, developed from the composite element method (CEM). The coupled relation between the fracture flow and normal stress makes use of the “filled model”, which examines the asperities in the fracture as a layer of granular medium having high porosity and being clipped by the two parallel plates. The existence of fractures is not considered in the mesh generation, but it will be considered explicitly in the mapped composite element. The coupled normal stress and fluid flow process has been simulated by applying a cross iterative algorithm between the two fields. The proposed algorithm considers not only the flow through the fractures, but also the flow exchange between fractures and the surrounding rock blocks. In addition, it can be used for both the filled and non-filled fractures. The verification of the proposed algorithm has been conducted through the illustration of three examples by comparison with the conventional finite element method (FEM), from which the advantages and reliability of the proposed algorithm have been shown clearly. 相似文献
19.
构造-流体-成矿体系的反应-输运-力学耦合模型和动力学模拟 总被引：12，自引：0，他引：12
建立了一个综合的构造流体成矿体系的反应输运力学耦合动力学模型。利用有限元方法求解岩石变形、断裂作用和断裂网络统计动力学、流体流动、有机和无机地球化学反应及成岩成矿作用、压力溶液和其它压实力学、热迁移的方程组 ,可以对构造流体成矿体系的动力学演化过程进行 1～ 3维数值模拟。模拟的主要内容是在各种过程耦合作用下描述构造流体成矿体系的主要变量的时空演化 :( 1)与成矿流体的形成和性质有关的变量 ,如地层中矿物 (包括成矿物质 )的溶解速率、流体中各组分的浓度与饱和度、流体温度、压力、离子强度等 ;( 2 )与构造变形和流体运移有关的各变量 ,如应力与变形速率、岩石孔隙度、构造 (断裂 )渗透率等 ;( 3 )与沉淀成矿有关的变量 ,如矿物 (金属矿物和脉石矿物 )的成核速率、各矿物的沉淀量等 ;( 4 )上述各有关变量间的时空耦合关系 ,如断裂渗透率时空演化与流体流动、汇聚和成矿的耦合关系等。以湖南沃溪金锑钨矿床为例 ,应用该模型和方法对成矿动力学过程和动力学机制进行了初步的模拟与分析。 相似文献
20.
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. 相似文献