共查询到20条相似文献，搜索用时 913 毫秒
1.
Microcrack-based coupled damage and flow modeling of fracturing evolution in permeable brittle rocks
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. 相似文献
2.
Yu-Shu Wu Guoping Lu Keni Zhang Lehua Pan Gudmundur S. Bodvarsson 《Hydrogeology Journal》2007,15(3):553-572
Characterizing percolation patterns in unsaturated fractured rock has posed a greater challenge to modeling investigations than comparable saturated zone studies due to the heterogeneous nature of unsaturated media and the great number of variables impacting unsaturated flow. An integrated modeling methodology has been developed for quantitatively characterizing percolation patterns in the unsaturated zone of Yucca Mountain, Nevada (USA), a proposed underground repository site for storing high-level radioactive waste. The approach integrates moisture, pneumatic, thermal, and isotopic geochemical field data into a comprehensive three-dimensional numerical model for analyses. It takes into account the coupled processes of fluid and heat flow and chemical isotopic transport in Yucca Mountain’s highly heterogeneous, unsaturated fractured tuffs. Modeling results are examined against different types of field-measured data and then used to evaluate different hydrogeological conceptualizations through analyzing flow patterns in the unsaturated zone. In particular, this model provides clearer understanding of percolation patterns and flow behavior through the unsaturated zone, both crucial issues in assessing repository performance. The integrated approach for quantifying Yucca Mountain’s flow system is demonstrated to provide a practical modeling tool for characterizing flow and transport processes in complex subsurface systems. 相似文献
3.
Stylolites are rough surfaces, formed by localized rock dissolution, and prevalent in carbonates and other sedimentary rocks. Their impact on porosity and permeability, and capacity to accommodate compactive strain, are well documented. This paper presents a meso-scale field study on sedimentary stylolites in carbonates, characterizing large-scale distributions of stylolites, including measurements conducted on longer than kilometer-long stylolites. Our field study suggests that on large scales connections between stylolites become important. Since connectivity, and also lack of connectivity, are expected to play a significant role in strain accommodation and hydraulic rock properties, we suggest that large-scale analysis may require a new characterization scheme for “stylolite populations”, based on their connectivity. We therefore divide sedimentary stylolite populations into three end-member types, which are correlated with the three possibilities for percolation of such systems: isolated stylolites (with zero percolation/connectivity), long-parallel stylolites (with 2-dimensional percolation/connectivity), and interconnected stylolite networks (with 3-dimensional percolation/connectivity). New statistical parameters and measures are devised and used to quantitatively characterize the different population types. Schematic mechanistic models are then offered to explain the evolution of the three end-member connectivity-classes. In addition we discuss the effect on fluid flow of the different population types. 相似文献
4.
〖HT2H〗低渗透砂岩油层相对渗透率曲线特征、影响因素及其对开发的影响 总被引：4，自引：0，他引：4
低渗透储层油、水两相渗透率特征及其影响因素的研究,是认识、分析低渗透储层油、水两相渗流机理的重要途径,是合理开发油田的基础。本文在总结低渗透油藏油水相对渗透率曲线特征的基础上,对以火山岩、凝灰岩岩屑为主的岩屑砂岩相对渗透率曲线特征进行研究,并分析了其影响因素,研究表明低渗透储层的岩石学特征,碎屑颗粒的成份、结构,其表面的粗糙程度,及其粘土矿物的产状对相对渗透率曲线有明显的影响。这一研究成果对合理开发此类低渗透油田有着十分重要的意义。 相似文献
5.
Po-Yu Chuang Yeeping Chia Ya-Hsuan Liou Mao-Hua Teng Ching-Yi Liu Tsai-Ping Lee 《Hydrogeology Journal》2016,24(7):1651-1662
Recent advances in borehole geophysical techniques have improved characterization of cross-hole fracture flow. The direct detection of preferential flow paths in fractured rock, however, remains to be resolved. In this study, a novel approach using nanoscale zero-valent iron (nZVI or ‘nano-iron’) as a tracer was developed for detecting fracture flow paths directly. Generally, only a few rock fractures are permeable while most are much less permeable. A heat-pulse flowmeter can be used to detect changes in flow velocity for delineating permeable fracture zones in the borehole and providing the design basis for the tracer test. When nano-iron particles are released in an injection well, they can migrate through the connecting permeable fracture and be attracted to a magnet array when arriving in an observation well. Such an attraction of incoming iron nanoparticles by the magnet can provide quantitative information for locating the position of the tracer inlet. A series of field experiments were conducted in two wells in fractured rock at a hydrogeological research station in Taiwan, to test the cross-hole migration of the nano-iron tracer through permeable connected fractures. The fluid conductivity recorded in the observation well confirmed the arrival of the injected nano-iron slurry. All of the iron nanoparticles attracted to the magnet array in the observation well were found at the depth of a permeable fracture zone delineated by the flowmeter. This study has demonstrated that integrating the nano-iron tracer test with flowmeter measurement has the potential to characterize preferential flow paths in fractured rock. 相似文献
6.
The effectiveness of secondary recovery methods in reservoir development studies depends on the knowledge about how fluid-carrying regions (i.e. good-quality rock types) are connected between injection and production wells. To estimate reservoir performance uncertainty, comprehensive simulations on many reservoir model realisations are necessary, which is very CPU consuming and time demanding. Alternatively, we can use much simpler and physically based methods such as percolation approach. Classic percolation assumes connectivity between opposite 2-D faces of a 3-D system; whereas, hydrocarbon production is achieved through active wells that are one-dimensional lines (e.g. vertical, horizontal or deviated wells). The main contribution of this study is to analyse the percolation properties of 3-D continuum percolation models with more realistic well representations during secondary recovery. In particular, the connection of randomly distributed sands (i.e. good-quality rock types) between two lines (representing two wells) located at two corners of the system are modelled by Monte Carlo simulations. Subsequently, the connectivity and conductivity of such a line-to-line well representation is compared with that of face-to-face well representations in the previously published results. The critical percolation properties of those systems as well as the universality concept are also investigated. As there are many rooms for connections in 3-D models, we found that the principal percolation properties will not be altered significantly when the problem with a face-to-face connection is transformed to a line-to-line connection model. 相似文献
7.
Bedoya-Gonzalez Diego Kessler Timo Rinder Thomas Schafmeister Maria-Theresia 《Hydrogeology Journal》2022,30(6):1717-1736
Underground hard coal mining usually disrupts the mechanical equilibrium of rock sequences, creating fractures within minor permeable rocks. The present study employs a dual-continuum model to assess how both fractured and porous sandstone media influence the percolation process in postmining setups. To test the approach, the software TOUGH2 was employed to simulate laminar fluid flow in the unsaturated zone of the Ibbenbüren Westfield mining area. Compared to other coal mining districts in Germany, this area is delineated by the topography and local geology, leading to a well-defined hydrogeological framework. Results reveal good agreement between the calculated and measured mine water discharge for the years 2008 and 2017. The constructed model was capable of reproducing the bimodal flow behavior of the adit by coupling a permeable fractured continuum with a low-conductivity rock matrix. While flow from the fractured continuum results in intense discharge events during winter months, the rock matrix determines a smooth discharge limb in summer. The study also evaluates the influence of individual and combined model parameters affecting the simulated curve. A detailed sensitivity analysis displayed the absolute and relative permeability function parameters of both continua among the most susceptible variables. However, a strong a priori knowledge of the value ranges for the matrix continuum helps to reduce the model ambiguity. This allowed for calibration of some of the fractured medium parameters for which sparse or variable data were available. However, the inclusion of the transport component and acquisition of more site-specific data is recommended to reduce their uncertainty.
相似文献8.
The Southern Alps of New Zealand is an actively rising mountain belt which displays a thermal anomaly adjacent to the Alpine Fault, the Australian-Pacific plate boundary. Extensive fluid movement occurs in this uplift zone, resulting in metallic vein mineralization. Gold mineralization is confined to greenschist facies rocks, while younger veins in amphibolite facies rocks near the Alpine Fault are enriched in copper. Transport and deposition of metals in this complex hydrothermal system is governed by interaction between rising metamorphic fluids and downward-percolating meteoric fluid. Metamorphic fluids have equilibrated with graphitic schist country rock and are relatively reduced. Infiltration and mixing of meteoric water increases oxygen activity and decreases sulphur activity in the fluid. Oxidised meteoric water heats up and dissolves Cu during downward percolation. This Cu is deposited as the fluid becomes more reduced. Hence, there is a progressive increase in copper content in the middle portions of the hydrothermal system, especially in the more permeable highly fractured rocks near the Alpine Fault. 相似文献
9.
Fractured rock has often been conceptualized as a dual-continuum system for many practical applications. This study proposes a systematic approach to deal with multiphase flow in a dual-continuum system. Considering that fluid flow occurs in pore volumes (including fracture apertures), we first develop a so-called pore-space conservation equation for deformed fractured rock and then combine this equation with fluid mass balance to derive governing equations for multiphase flow associated with rock deformation. Constitutive relationships are also presented for describing stress dependence of hydraulic properties and effective mechanical parameters for bulk rock body (as a function of the corresponding parameters for fracture and matrix continua). Finally, we applied the developed approach to a CO2 geological sequestration problem to demonstrate the usefulness of the approach. 相似文献
10.
裂隙岩体渗流概念模型研究 总被引：31，自引：4，他引：27
裂隙岩体中的渗流和传统的多孔介质渗流在机理上存在本质的差别,这种差别主要表现为裂隙岩体在各种尺度上存在的非均质性。模拟裂隙岩体渗流的主要困难在于描述这种非均质性。目前的概念模型,包括等效连续体模型、离散裂隙网络模型和混合模型使用了不同的技术来预测裂隙岩体中的渗流。这些模型基于不同的假设和概念框架,有着各自不同的优缺点。在实际应用时,应当根据研究域的具体特点和所要解决的问题的要求对其选择,此外,还讨论了单裂隙的概念模型。 相似文献
11.
12.
One-dimensional advection-dispersion models predict that characteristic δ18 O vs. distance and δ18 O vs. δ13 C profiles should be produced during isothermal metamorphic fluid flow under equilibrium conditions. However, the patterns of isotopic resetting in rocks that have experienced fluid flow are often different from the predictions. Two-dimensional advection-dispersion simulations in systems with simple geometries suggest that such differences may be as a result of fluid channelling and need not indicate disequilibrium, high dispersivities, or polythermal flow. The patterns of isotopic resetting are a function of: (1) the permeability contrast between more permeable layers ('channels') and less permeable layers ('matrix'); (2) the width and spacing of the channels; (3) the width and spacing of discrete fractures; and (4) the orientation of the pressure gradient with respect to layering. In fractured systems, the efficiency of isotopic transport depends on the fracture aperture and the permeability of the surrounding rock. Resetting initially occurs along and immediately adjacent to the fractures, but with time isotopic resetting because of flow through the rock as a whole increases in importance. Application of the one-dimensional advection-dispersion equations to metamorphic fluid flow systems may yield incorrect estimates of fluid fluxes, intrinsic permeabilities, dispersivities, and permeability contrasts unless fluid flow occurred through zones of high permeability that were separated by relatively impermeable layers. 相似文献
13.
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. 相似文献
14.
Connectivity is an important measure for assessing flow transport in rock, especially through fractures. In this paper, rock
fracture systems are modelled by a discrete fracture model simulated by a marked point process. A connectivity index is then
introduced to quantify the connectivity between any two points in space. Monte Carlo simulation is used to evaluate the connectivity
index for stationary cases and relationships between the connectivity index and the parameters of the discrete fracture model
are analysed. The average number of intersections per fracture, Xf, and the fracture intensity, P12 (P32), are calculated and the relationships between these parameters and the connectivity index are investigated, concluding that
Xf is the more suitable parameter for the classification of rock mass flow properties. The relationships between the connectivity
index and the percolation state of the fractured medium are also discussed. An edge correction is briefly discussed and a
practical example is used to demonstrate the method of computing the connectivity index. 相似文献
15.
16.
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. 相似文献
17.
A comparison of experimental and numerical results is presented addressing two-phase immiscible displacement in rough fractures. Quasi-static, capillary displacement in a rough fracture is modeled using the modified invasion percolation approach proposed by Glass et al (1998), and the results are compared against experimental observations obtained from two-phase flow through a rock fracture using x-ray computed-tomography scanning. The model is based on an algorithm seeking the least resistant pathway for the advancement of the invading fluid using the Young-Laplace equation and accounting for local in-plane curvature of the advancing fluid front. The saturation distribution map generated by the model yields good agreement with the experimental phase distribution and presents more realistic phase structures than those obtained from the conventional invasion percolation approach. The improvement in the results obtained with the modified invasion percolation approach is attributed to the contribution of the in-plane curvature term, which captures the effect of regionalized apertures, rather than single-point apertures, on the shape of the invading front. Glass RJ, Nicholl MJ, Yarrington L (1998) A modified invasion percolation model for low-capillary number immiscible displacements in horizontal rough-walled fractures: influence of local in-plane curvature. Water Resour Res 34:3215–3234 相似文献
18.
To predict the macroscopic properties (e.g., transport, electromagnetic, and mechanical properties) of porous media, it is necessary to have a three‐dimensional (3D) representation of porous media. We reconstruct the geologically realistic 3D structure of Fontainebleau sandstone based on the two‐dimensional (2D) thin sections by using the multiple‐point statistics method. For this method, the size of template is an important parameter that reflects the perceived scale of spatial structure of porous media. In this paper, we take advantage of entropy method to obtain the appropriate size of the template, which is proven to be correct and feasible. The reconstruction method proposed by us combines successive 2D MPS simulations as well as 3D MPS simulation, which takes account into the pore structure information (e.g., heterogeneity and connectivity) both intralayer and interlayer. This reconstruction method is tested on Fontainebleau sandstone for which 3D images from micro‐CT scanning are available. Applying local percolation theory analysis, this new approach can depict the expected patterns of geological heterogeneities. In addition, it also can well reproduce a high degree of connectivity of the pore space better than other reconstruction methods based on lower‐order statistics. Copyright © 2011 John Wiley & Sons, Ltd. 相似文献
19.
Numerical simulation of fluid flow coupled with chemical reactions has been an active field in the hydrogeology community and many formulations have been programmed into different software. In recent years, this subject has attracted increasing interest in the reservoir simulation community, partly for the application of chemical methods for hydrocarbon extraction but also for research on the geological sequestration of CO2. In this paper, an extension to the concept of dual mesh for reactive transport modeling is presented. This approach involves two meshes, a low-resolution mesh to resolve the pressure equation and a high-resolution mesh to transport the species and to calculate the geochemical equilibrium. The main objective is to preserve the fine scale heterogeneities to reach a more accurate field behavior simulation than conventional approach which consist in performing simulations on a coarser mesh. The method is applied to a simulation of CO2 storage in the SPE10 model that keep a high resolution of the heterogeneities. 相似文献
20.
An Approach for Modeling Rock Discontinuous Mechanical Behavior Under Multiphase Fluid Flow Conditions 总被引：1，自引：0，他引：1
Peng-Zhi Pan Jonny Rutqvist Xia-Ting Feng Fei Yan 《Rock Mechanics and Rock Engineering》2014,47(2):589-603
In this paper, the two computer codes TOUGH2 and RDCA (for “rock discontinuous cellular automaton”) are integrated for coupled hydromechanical analysis of multiphase fluid flow and discontinuous mechanical behavior in heterogeneous rock. TOUGH2 is a well-established code for geohydrological analysis involving multiphase, multicomponent fluid flow and heat transport; RDCA is a numerical model developed for simulating the nonlinear and discontinuous geomechanical behavior of rock. The RDCA incorporates the discontinuity of a fracture independently of the mesh, such that the fracture can be arbitrarily located within an element, while the fluid pressure calculated by TOUGH2 can be conveniently applied to fracture surfaces. We verify and demonstrate the coupled TOUGH–RDCA simulator by modeling a number of simulation examples related to coupled multiphase flow and geomechanical processes associated with the deep geological storage of carbon dioxide—including modeling of ground surface uplift, stress-dependent permeability, and the coupled multiphase flow and geomechanical behavior of fractures intersecting the caprock. 相似文献