共查询到20条相似文献，搜索用时 31 毫秒
1.
An efficient and accurate numerical model for multicomponent compressible single-phase flow in fractured media is presented. The discrete-fracture approach is used to model the fractures where the fracture entities are described explicitly in the computational domain. We use the concept of cross flow equilibrium in the fractures. This will allow large matrix elements in the neighborhood of the fractures and considerable speed up of the algorithm. We use an implicit finite volume (FV) scheme to solve the species mass balance equation in the fractures. This step avoids the use of Courant–Freidricks–Levy (CFL) condition and contributes to significant speed up of the code. The hybrid mixed finite element method (MFE) is used to solve for the velocity in both the matrix and the fractures coupled with the discontinuous Galerkin (DG) method to solve the species transport equations in the matrix. Four numerical examples are presented to demonstrate the robustness and efficiency of the proposed model. We show that the combination of the fracture cross-flow equilibrium and the implicit composition calculation in the fractures increase the computational speed 20–130 times in 2D. In 3D, one may expect even a higher computational efficiency. 相似文献
2.
《Advances in water resources》2005,28(3):235-255
The main objective of this work is to compare three different models for modelling of flow and solute transport in fractured porous media, in terms of their predictions of the flow and solute transport field variables. The three models are: the equivalent continuum model, the dual porosity model and the discrete fracture/non-homogeneous model. Though it is clear that the three models are based on different assumptions for their validity, it is not clear in which cases two or all of them would give similar results, since there are no such reported comparisons in the open literature.The three methods are compared for two different geometries: a rectangular porous domain with two parallel fractures and a square porous domain with regular mesh of three parallel fractures and another three fractures perpendicular to the first ones. The results helped to draw some conclusions in respect to the similarity of potentials as well as fluxes for the different methods for each of the two geometries.In this research the boundary element dual reciprocity method–multi domain scheme (BE DRM–MD) has been used and its implementation has been described. This numerical scheme has been used for the first time to solve a dual-porosity model. The scheme showed satisfactory accuracy and high flexibility in preparation of the discrete fracture/non-homogeneous meshes. 相似文献
3.
The presence of fractures in fluid‐saturated porous rocks is usually associated with strong seismic P‐wave attenuation and velocity dispersion. This energy dissipation can be caused by oscillatory wave‐induced fluid pressure diffusion between the fractures and the host rock, an intrinsic attenuation mechanism generally referred to as wave‐induced fluid flow. Geological observations suggest that fracture surfaces are highly irregular at the millimetre and sub‐millimetre scale, which finds its expression in geometrical and mechanical complexities of the contact area between the fracture faces. It is well known that contact areas strongly affect the overall mechanical fracture properties. However, existing models for seismic attenuation and velocity dispersion in fractured rocks neglect this complexity. In this work, we explore the effects of fracture contact areas on seismic P‐wave attenuation and velocity dispersion using oscillatory relaxation simulations based on quasi‐static poroelastic equations. We verify that the geometrical and mechanical details of fracture contact areas have a strong impact on seismic signatures. In addition, our numerical approach allows us to quantify the vertical solid displacement jump across fractures, the key quantity in the linear slip theory. We find that the displacement jump is strongly affected by the geometrical details of the fracture contact area and, due to the oscillatory fluid pressure diffusion process, is complex‐valued and frequency‐dependent. By using laboratory measurements of stress‐induced changes in the fracture contact area, we relate seismic attenuation and dispersion to the effective stress. The corresponding results do indeed indicate that seismic attenuation and phase velocity may constitute useful attributes to constrain the effective stress. Alternatively, knowledge of the effective stress may help to identify the regions in which wave induced fluid flow is expected to be the dominant attenuation mechanism. 相似文献
4.
The Kirkham piezometer tube method of measuring hydraulic conductivity of saturated soil uses a numerical constant (shape factor) to account for the relative geometrical parameters of the cavity and the soil domain. Published shape factors have mostly been determined using an electrical analogue. In this paper the application of the finite element method to shape factor determination is considered. Convergence is studied and methods of extrapolation evaluated. A determination is made of the size of finite region required to represent adequately the radially unbounded domain assumed for the piezometer method. Computed shape factors were found to be in good agreement with analogue data. 相似文献
5.
In general, the accuracy of numerical simulations is determined by spatial and temporal discretization levels. In fractured porous media, the time step size is a key factor in controlling the solution accuracy for a given spatial discretization. If the time step size is restricted by the relatively rapid responses in the fracture domain to maintain an acceptable level of accuracy in the entire simulation domain, the matrix tends to be temporally over-discretized. Implicit sub-time stepping applies smaller sub-time steps only to the sub-domain where the accuracy requirements are less tolerant and is most suitable for problems where the response is high in only a small portion of the domain, such as within and near the fractures in fractured porous media. It is demonstrated with illustrative examples that implicit sub-time stepping can significantly improve the simulation efficiency with minimal loss in accuracy when simulating flow and transport in fractured porous media. The methodology is successfully applied to density-dependent flow and transport simulations in a Canadian Shield environment, where the flow and transport is dominated by discrete, highly conductive fracture zones. 相似文献
6.
交错网格高阶差分方法是一种在保持效率的前提下提高弹性波模拟精度的有效方法.本文将可变空间网格与变化的时间步长技术引入到交错网格高阶差分弹性波模拟中,提出一种空间网格可任意奇数倍变化与时间步长任意变化的交错网格高阶差分弹性波模拟方法.一系列数值试验表明,该方法能够在保证模拟精度的同时,通过有效降低空间与时间维度上的过采样来显著提高弹性波模拟的效率.同时,该方法还能够精细刻画含孔缝洞介质以及横向变化剧烈介质的局部细微结构,减小弹性波模拟误差,提高介质细微结构处的弹性波传播模拟精度. 相似文献
7.
K.F. TYAPKIN 《Geophysical Prospecting》1980,28(6):935-944
According to data presented by YU.A. Kosigin, 84% of all endogenous mineral deposits are in some direct or indirect relation with the fractures of the earth's crust. Therefore the discovery and the study of the spatial disposition of the fractures is the most important object for geophysicists. Abyssal fractures are of particular interest. By geophysical methods one can find the geometrical parameters of abyssal fractures such as their extent, the depth of formation, the breadth of the zones, and the amplitude of the relative displacement of separated blocks. The methods determining these parameters are widely known. A calculation of the difference in the levels of the erosion cuts of the blocks are of particular interest under shield conditions. A method to calculate this difference by gravity interpretation of “step” anomalies and by using the gradient model of the earth's crust is proposed. A comparision of the results of gravimetry and magnetometry with those of deep seismic soundings shows that the fractures of the first and the second order in the Ukrainian shield cut the earth's crust and part of the upper mantle so they can be qualified as abyssal fractures. In the Ukrainian shield the spatial regularities of the fractures determined by geophysical methods can be used for the prognosis of the ore deposits. 相似文献
8.
The effect of fracture network geometry on free convection in fractured rock is studied using numerical simulations. We examine the structural properties of fracture networks that control the onset and strength of free convection and the patterns of density-dependent flow. Applicability of the equivalent porous medium approach (EPM) is also tested, and recommendations are given, for which situations the EPM approach is valid. To date, the structural properties of fracture networks that determine free convective flow are examined only in few, predominantly simplified regular fracture networks. We consider fracture networks containing continuous, discontinuous, orthogonal and/or inclined discrete fractures embedded in a low-permeability rock matrix. The results indicate that bulk permeability is not adequate to infer the occurrence and magnitude of free convection in fractured rock. Fracture networks can inhibit or promote convection depending on the fracture network geometry. Continuous fracture circuits are the crucial geometrical feature of fracture networks, because large continuous fracture circuits with a large vertical extent promote convection. The likelihood of continuous fracture circuits and thus of free convection increases with increasing fracture density and fracture length, but individual fracture locations may result in great deviances in strength of convection between statistically equivalent fracture networks such that prediction remains subject to large uncertainty. 相似文献
9.
《Earth and Planetary Science Letters》2006,241(3-4):454-465
An understanding of fluid flow through natural fractures in rocks is important in many areas, such as in the hydrocarbon and water industries, and in the safe design of disposal sites for domestic, industrial and nuclear waste. It is often impractical to obtain this information by field or laboratory scale measurements, so numerical modelling of fluid flow must be carried out using synthetic fractures with rough fracture surfaces that are representative of the natural rock fractures. Clearly there are two practical requirements; (i) the development of a method for analysing natural rock fractures to obtain their characteristic parameters, and (ii) the development of techniques for creating high quality synthetic fractures using these parameters. We have implemented these practical requirements in two new software packages. The first, ParaFrac allows the analysis and parameterisation of fracture surfaces and apertures. The second, SynFrac, enables the numerical synthesis of fracture surfaces and apertures with basic prescribed parameters. Synthetic fractures are created using, (i) a new model, which takes full account of the complex matching properties of fracture surfaces using two new parameters, a minimum matching fraction and a transition length and (ii) an improved method of partially correlated random number generation. This model more closely captures the often complex matching properties of real rock fractures than previous more simplified models. 相似文献
10.
11.
Generalized non‐hyperbolic approximation for qP‐wave relative geometrical spreading in a layered transversely isotropic medium with a vertical symmetry axis 
下载免费PDF全文
下载免费PDF全文 Compensation for geometrical spreading along the ray‐path is important in amplitude variation with offset analysis especially for not strongly attenuative media since it contributes to the seismic amplitude preservation. The P‐wave geometrical spreading factor is described by a non‐hyperbolic moveout approximation using the traveltime parameters that can be estimated from the velocity analysis. We extend the P‐wave relative geometrical spreading approximation from the rational form to the generalized non‐hyperbolic form in a transversely isotropic medium with a vertical symmetry axis. The acoustic approximation is used to reduce the number of parameters. The proposed generalized non‐hyperbolic approximation is developed with parameters defined by two rays: vertical and a reference rays. For numerical examples, we consider two choices for parameter selection by using two specific orientations for reference ray. We observe from the numerical tests that the proposed generalized non‐hyperbolic approximation gives more accurate results in both homogeneous and multi‐layered models than the rational counterpart. 相似文献
12.
Computed tomography scan imaging of natural open fractures in a porous rock; geometry and fluid flow 总被引：2，自引：0，他引：2
Computed tomography scan imaging techniques have been used on core samples to investigate the effect of natural open fractures on reservoir flow in the Snøhvit Gas Condensate Field. Firstly, computed tomography (CT) scanning was used to describe the 3D geometrical properties of the fracture network including orientation and fracture density. Two types of fractures were observed: F1 fractures are short and stylolite related and F2 fractures are longer, cross-cutting the core and without any obvious relationship to stylolites. Secondly, monitoring of single and two phase flow experiments on samples containing these two types of natural open fractures was performed under 10 and 80 bar net confining pressure while using CT scanning. 1-phase miscible flooding experiment shows approximately 3 times higher flooding velocity in an open F2 fracture than in the matrix. 2-phase flooding by gas injection into a 100% water saturated core gave early gas breakthrough due to flow in the fracture and thereafter very little water production. The flow experiments showed that the presence of open fractures has a significant local effect on fluid flow even in a case with relatively high matrix porosity (200–300 mD). The sample containing F1 fractures showed a complex flow pattern influenced both by open fractures and stylolites. The CT scan data enables an exact representation of the fracture network in core scale simulation models and therefore improves the understanding of fracture influence on flow in a fractured porous medium. CT scanning of core samples provides an effective tool for integrating geology and fluid flow properties of a porous fractured medium. 相似文献
13.
Effects of the geometry of two‐dimensional fractures on their hydraulic aperture and on the validity of the local cubic law 下载免费PDF全文
下载免费PDF全文 Flow through rough fractures is investigated numerically in order to assess the validity of the local cubic law for different fracture geometries. Two‐dimensional channels with sinusoidal walls having different geometrical properties defined by the aperture, the amplitude, and the wavelength of the walls' corrugations, the corrugations asymmetry, and the phase shift between the two walls are considered to represent different fracture geometries. First, it is analytically shown that the hydraulic aperture clearly deviates from the mean aperture when the walls' roughness, the phase shift, and/or the asymmetry between the fracture walls are relatively high. The continuity and the Navier–Stokes equations are then solved by means of the finite element method and the numerical solutions compared to the theoretical predictions of the local cubic law. Reynolds numbers ranging from 0.066 to 66.66 are investigated so as to focus more particularly on the effect of flow inertial effects on the validity of the local cubic law. For low Reynolds number, typically less than 15, the local cubic law properly describes the fracture flow, especially when the fracture walls have small corrugation amplitudes. For Reynolds numbers higher than 15, the local cubic law is valid under the conditions that the fracture presents a low aspect ratio, small corrugation amplitudes, and a moderate phase lag between its walls. 相似文献
14.
In this study, a numerical manifold method (NMM) model is developed to analyze flow in porous media with discrete fractures in a non-conforming mesh. This new model is based on a two-cover-mesh system with a uniform triangular mathematical mesh and boundary/fracture-divided physical covers, where local independent cover functions are defined. The overlapping parts of the physical covers are elements where the global approximation is defined by the weighted average of the physical cover functions. The mesh is generated by a tree-cutting algorithm. A new model that does not introduce additional degrees of freedom (DOF) for fractures was developed for fluid flow in fractures. The fracture surfaces that belong to different physical covers are used to represent fracture flow in the direction of the fractures. In the direction normal to the fractures, the fracture surfaces are regarded as Dirichlet boundaries to exchange fluxes with the rock matrix. Furthermore, fractures that intersect with Dirichlet or Neumann boundaries are considered. Simulation examples are designed to verify the efficiency of the tree-cutting algorithm, the calculation's independency from the mesh orientation, and accuracy when modeling porous media that contain fractures with multiple intersections and different orientations. The simulation results show good agreement with available analytical solutions. Finally, the model is applied to cases that involve nine intersecting fractures and a complex network of 100 fractures, both of which achieve reasonable results. The new model is very practical for modeling flow in fractured porous media, even for a geometrically complex fracture network with large hydraulic conductivity contrasts between fractures and the matrix. 相似文献
15.
Volcanic vents rooted on extension fractures and their geodynamic implications in the Ethiopian Rift
Tesfaye Korme Jean Chorowicz Bernard Collet Franco F. Bonavia 《Journal of Volcanology and Geothermal Research》1997,79(3-4)
Most of the extension fractures located in the Ethiopian rift are related to rift dynamics, and they have consistently similar orientations for hundreds of kilometers. The orientation of extension fracture gives the local extension direction which is perpendicular to its horns (end segments of an extension fracture). We have established clear geometrical relationships between tectonics and volcanism in the Main Ethiopian Rift by interpreting high-resolution images obtained from Landsat-TM and SPOT satellites and by quantitatively analyzing the geometry of extension fractures, elongated vents and linear volcanic clusters. Applying the relationships obtained, we show that extension fractures in the rift served as channels for magma rising to the surface, and that extension fractures underlie most of the elongated volcanic vents and linear volcanic clusters. The geometry of extension fractures beneath volcanic edifices can be deduced from the shape of elongated vents and pattern of linear volcanic clusters rooted on them. Utilizing the orientations of extension fractures directly observed and those inferred from elongated volcanic vents and linear clusters, we found that the extension direction of the Main Ethiopian Rift is northwest–southeast and that the direction has been rotated clockwise for about 20° in the time interval 2.83 to 0.023 Ma. The Recent axis of rift opening is oriented N40° and located closer to the southeastern escarpment giving an asymmetric geometry to the Main Ethiopian Rift (MER). 相似文献
16.
Numerical investigation of alternative fracture stiffness measures and their respective scaling behaviours 下载免费PDF全文
下载免费PDF全文 We study the mechanical deformation of fractures under normal stress, via tangent and specific fracture stiffnesses, for different length scales using numerical simulations and analytical insights. First, we revisit an equivalent elastic layer model that leads to two expressions: the tangent stiffness is the sum of an “intrinsic” stiffness and the normal stress, and the specific stiffness is the tangent stiffness divided by the fracture aperture at current stress. Second, we simulate the deformation of rough fractures using a boundary element method where fracture surfaces represented by elastic asperities on an elastic half‐space follow a self‐affine distribution. A large number of statistically identical “parent” fractures are generated, from which sub‐fractures of smaller dimensions are extracted. The self‐affine distribution implies that the stress‐free fracture aperture increases with fracture length with a power law in agreement with the chosen Hurst exponent. All simulated fractures exhibit an increase in the specific stiffness with stress and an average decrease with increase in length consistent with field observations. The simulated specific and tangent stiffnesses are well described by the equivalent layer model provided the “intrinsic” stiffness slightly decreases with fracture length following a power law. By combining numerical simulations and the analytical model, the effect of scale and stress on fracture stiffness measures can be easily separated using the concept of “intrinsic” stiffness. We learn that the primary reason for the variability in specific stiffness with length comes from the fact that the typical aperture of the self‐affine fractures itself scales with the length of the fractures. 相似文献
17.
Weak Zone Related Seismic Cycles in Progressive Failure Leading to Collapse in Brittle Crust 总被引：2，自引：0，他引：2
Until quite recently our understanding of the basic mechanical process responsible for earthquakes and faulting was not well known. It can be argued that this was partly a consequence of the complex nature of fracture in crust and in part because evidence of brittle phenomena in the natural laboratory of the earth is often obliterated or obscured by other geological processes. While it is well understood that the spatial and temporal complexity of earthquakes and the fault structures emerge from geometrical and material built-in heterogeneities, one important open question is how the shearing becomes localized into a band of intense fractures. Here we address these questions through a numerical approach of a tectonic plate by considering rockmass heterogeneity both in microscopic scale and in mesoscopic scale. Numerical simulations of the progressive failure leading to collapse under long-range slow driving forces in the far-field show earthquake-like rupture behavior. En Echelon crack-arrays are reproduced in the numerical simulation. It is demonstrated that the underlying fracturing induced acoustic emissions (or seismic events) display self-organized criticality—from disorder to order. The seismic cycles and the geometric structures of the fracture faces, which are found greatly depending on the material heterogeneity (especially on the macroscopic scale), agree with that observed experimentally in real brittle materials. It is concluded that in order to predict a main shock, one must have extremely detailed knowledge on very minor features of the earth's crust far from the place where the earthquake originated. If correct, the model proposed here seemingly provides an explanation as to why earthquakes to date are not predicted so successfully. The reason is not that we do not understand earthquake mechanisms very well but that we still know little about our earth's crust. 相似文献
18.
A systematic numerical method has been presented to investigate the constitutive relationships between two-phase flow properties of horizontal fractures and aperture distributions. Based on fractal geometry, single rough-walled fractures are generated numerically by modified successive random addition (SRA) method and then aperture distributions with truncated Gaussian distribution are formed by shear displacement between lower and upper surfaces. (The truncated Gaussian distribution is used to describe aperture evolution under different normal stresses.) According to the assumption of two-dimensional porous media and local parallel plate model, invasion percolation approach is employed to model the two-phase flow displacement (imbibition) in generated horizontal fractures, in which capillary forces are dominant over viscous and gravity forces. For truncated Gaussian distributions, constitutive relationships from numerical simulation are compared to closed-form relationships and a good agreement is obtained. The simulation results indicate strong phase interference with the sum of two phase relative permeability values being less than one in the intermediate saturations. It is found that fracture properties related to residual saturations depend on spatial correlation of aperture distributions. Based on the simulation results, we proposed an empirical relationship between the fracture residual-saturation-rated parameters and the corresponding aperture distributions. 相似文献
19.
— A new method for characterising the detailed fracture geometry in sandstone cores is presented. This method is based on the impregnation of samples with coloured resin, without significant disturbance of the fractures. The fractures are made clearly visible by the resin, thus allowing the fracture geometry to be examined digitally. In order to model the bulk hydraulic conductivity, the samples are sectioned serially perpendicular to the flow direction. The hydraulic conductivity of individual sections is estimated by summing the contribution of the matrix and each fracture from the digital data. Finally, the hydraulic conductivity of the bulk sample is estimated by a harmonic average in series along the flow path. Results of this geometrical method are compared with actual physical conductivity values measured from fluid experiments carried out prior to sectioning. The predicted conductivity from the fracture geometry parameters (e.g., fracture aperture, fracture width, fracture length and fracture relative roughness all measured using an optical method) is in good agreement with the independent physical measurements, thereby validating the approach. 相似文献
20.
This paper discusses the resolution of geometrical characteristics of pressurized fractures from tiltmeter data. The quasi-static
deformation and tilt field induced by such fractures can be modeled by superposition of displacement discontinuity (DD) singularities.
Despite the relatively common use of such measurements to infer fracture characteristics, there is a widespread misunderstanding
of what can be accurately determined, depending on the relative distance between the tiltmeter array and the fracture. We
investigate in detail the resolution of the dimensions and orientation of hydraulic fractures or faults from tilt measurements.
In particular, we formally prove that at a distance larger than about twice the characteristic length of the fracture, elastostatic
measurements such as those measured by tiltmeters are not able to resolve independently all the dimensions of the fracture,
although the fracture volume can be robustly inverted from the data. The resolution of fracture orientation is also discussed
using an analysis based on a spatial Fourier Transform of the tilt field. The relative angle between the plane where the measurements
are located and the fracture plane plays a major role in the accuracy of this estimation. In an illustrative field example,
where the measurements are located in the far-field of the fracture deformation field, we show how a single DD singularity
can be used to model tiltmeter data and efficiently obtain the fracture orientation and volume. 相似文献