首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 46 毫秒
1.
The no-slip boundary condition has usually been assumed to hold for the Reynolds equations (local cubic law) for fluid flow through rough-walled fractures. However, its validity for non-wetting fluid flow, such as prevails in fractured oil reservoirs, has been questioned. A series of experiments with a rough-walled fracture with mean aperture of 760 μm finds a higher flow rate for non-wetting fluid than wetting fluid. A modified Reynolds equation with a slip boundary condition is derived for non-wetting fluid flow through rough-walled fractures. Comparison of the modified Reynolds equation predictions with experimental results confirms that slip was a plausible explanation for a higher flow rate. The amount by which the flow rate for non-wetting fluids exceeds that for wetting fluids is found to depend highly on, and increase with, the degree to which the flowing fluid was non-wetted to thin immobile films on the surfaces. Numerical studies using the modified Reynolds equation indicate that the flow rate of non-wetting fluid became higher than that of wetting fluid as the roughness of the fracture increases. As the aperture becomes smaller, the flow rate ratio of non-wetting fluid to wetting fluid becomes large, leading to the endpoint relative permeability for the non-wetting fluid to exceed 1. The experimental and numerical studies clearly show that as the aperture of the fracture became less than a few hundred microns, the modified Reynolds equation with slip boundary conditions provides a better model for flow of a non-wetting fluid through rough-walled fractures.  相似文献   

2.
To investigate the influence of compression, Poisson effect and turbulence on the fluid flow process and the inversion for fracture surface geometries, we simulate two sets of fractures: one with a defined fracture height standard deviation σ constant and a varying autocorrelation length λ and another with a fixed λ and a changing σ. Under compression, the normal stress closes fractures with a large aperture and thus reduces the effective permeability. However, the Poisson effect, which is induced by the compression, has little influence on the fluid flow properties and does not affect the inversion for fracture height standard deviation or the autocorrelation length. When introducing turbulence, we observe a significant difference between the performance of the Navier–Stokes equation and the local cubic law; compared with the Navier–Stokes equation, the local cubic law overestimates the peak value of the breakthrough time curve and effective permeability, thereby underestimating the mean fracture aperture.  相似文献   

3.
— Experimental observations show that flow through a fracture decreases more rapidly than the cube of the mean aperture (Cook, 1992). In order to provide a possible explanation of these experimental findings, we study creeping flow through a fracture of varying aperture that is symmetric about its midplane, using the power series of the stream function obtained by Van Dyke (1987) for low Reynolds numbers. For the case of sinusoidally-varying walls, a simple expression relating the effective hydraulic aperture of the channel to the mean aperture and to the amplitude and wavelength of the sinusoidal wall profiles is obtained. Comparison is made to previous studies (Kitanidis and Dykaar, 1997) and to finite element calculations, and good agreement is obtained. The effect of fracture closure is then modelled as a decrease of the mean aperture without a change in the roughness. A power law relationship can be obtained between the flowrate and the mean aperture, with an exponent as high as 10, thus providing a potential mechanistic explanation of the experimental findings of Pyrak-Nolte et al. (1987).  相似文献   

4.
《水文科学杂志》2013,58(6):1125-1138
Abstract

Type curves are derived analytically for radial flow in rough horizontal fractures toward a well. The basic assumptions are that there is no turbulent flow near the borehole and the well storage is ignored. The basis of the methodology is to write explicit expressions for the continuity and cubic law flow equations, which are combined using a Boltzmann transformation leading to a simple ordinary differential equation for groundwater movement. Solutions are presented as a set of type curves for different fracture apertures. It is observed that the solutions provide a method of uniquely identifying fracture hydraulic parameters when the fracture is smooth, but pose ambiguity for rough fracture parameter estimations. However, large time portions of these type curves appear as straight lines on semi-logarithmic paper, which provides a unique way for rough fracture parameter determination. Identification of the fracture parameters, namely, the aperture and relative roughness, is possible in a unique manner with the use of these lines and the dimensionless time drawdown concept. The cubic law is the asymptotic behaviour, either for large times or large fracture apertures. Prior to this asymptotic part, there is a non-cubic portion which gives rise to systematic deviations from the cubic law. The technique presented is useful, especially for evaluating pumping tests from a single major fracture isolated by packers.  相似文献   

5.
Modeling fracture porosity development using simple growth laws   总被引：2，自引：0，他引：2  
A model of porosity development has been developed to investigate general relationships between simple fracture aperture growth laws and fracture porosity in evolved fracture arrays in aquifers. The growth of fracture apertures in two-dimensional orthogonal arrays with initially spatially uncorrelated lognormal aperture distributions has been studied, where aperture growth rate is proportional to an exponent of the flow rate through each fracture. The evolved arrays show geometrical phase changes as a function of the aperture growth rate exponent, e, and the standard deviation of the initial aperture distribution, sigma(z). Low values of e and sigma(z) lead to bimodal aperture distributions, where apertures parallel to flow are preferentially enlarged. At moderate values of e and sigma(z), there is a transition to a regime of more complex geometries consisting of networks of channel-like structures of preferentially enlarged apertures. At larger values of e, array-spanning channel-like paths of preferentially enlarged apertures develop, where the tortuosity of the channel-like paths is a linear function of sigma(z). Following an initial growth phase, during which dynamically stable aperture configurations develop, arrays undergo simple amplification. The geometry of the evolved aperture fields is diverse and they can be highly complex; consequently, parameterization and prediction of their evolution in terms of the initial aperture distributions and growth rate laws is not trivial.  相似文献   

6.
A quasi-three-dimensional particle tracking model is developed to characterize the spatial and temporal effects of advection, molecular diffusion, Taylor dispersion, fracture wall deposition, matrix diffusion, and co-transport processes on two discrete plumes (suspended monodisperse or polydisperse colloids and dissolved contaminants) flowing through a variable aperture fracture situated in a porous medium. Contaminants travel by advection and diffusion and may sorb onto fracture walls and colloid particles, as well as diffuse into and sorb onto the surrounding porous rock matrix. A kinetic isotherm describes contaminant sorption onto colloids and sorbed contaminants assume the unique transport properties of colloids. Sorption of the contaminants that have diffused into the matrix is governed by a first-order kinetic reaction. Colloids travel by advection and diffusion and may attach onto fracture walls; however, they do not penetrate the rock matrix. A probabilistic form of the Boltzmann law describes filtration of both colloids and contaminants on fracture walls. Ensemble-averaged breakthrough curves of many fracture realizations are used to compare arrival times of colloid and contaminant plumes at the fracture outlet. Results show that the presence of colloids enhances contaminant transport (decreased residence times) while matrix diffusion and sorption onto fracture walls retard the transport of contaminants. Model simulations with the polydisperse colloids show increased effects of co-transport processes.  相似文献   

7.
— A numerical fracture flow simulation based on the lubrication approximation is used to investigate the influence of roughness on the flow inside a rough fracture, at low Reynolds number. Facing surfaces are described as self-affine topographies with identical roughness magnitude. Resolution of the Reynolds equation is achieved using two distinct numerical schemes, with consistency. Fracture closure is studied assuming perfect plastic contact between facing surfaces. Long-range correlations are shown to exist in the local aperture field due to the fracture geometry and subsequently in the local fluxes inside the fracture. Flow channeling is the result of these correlations in terms of spatial distribution of the flow, and is responsible for either flow-enhancing or flow-inhibiting behavior of the fracture. Matching between the two surfaces at scales larger than a mismatch scale is studied. The mismatch scale is the upper limit scale for the local apertures scale invariance. It appears to control flow channeling and the related dispersion of the possible behaviors over a large statistics of fractures with identical statistical features. Hydraulic anisotropy of a given fracture is investigated: the dependence of the fracture transmittivity on the pressure drop orientation is proved to be sinusoidal, with an amplitude that is controlled by the mismatch scale.  相似文献   

8.
利用国内外287次5级以上的实际地震记录,研究了地震动有效幅值的不对称性特征及其规律。分别用有效幅值不对称比和有效幅值作用次数不对称比来描述地震动的不对称性,统计结果具有规律性,并提出了其估计公式。结果表明:地震动有效幅值不对称性和有效幅值作用次数不对称性特征明显,幅值不对称比的平均值不受震级影响,但与截断比呈非线性递增关系,截断比小于1/3时不对称比在1.07～1.20之间,截断比大于1/3时不对称比在1.13～1.96之间;有效幅值作用次数不对称比的平均值不受震级影响,但与截断比呈非线性递增关系,截断比小于1/3时作用次数不对称比在1.04～1.31之间变化,截断比大于1/3时作用次数不对称比在1.22～1.80之间。  相似文献   

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

10.
Local cubic law (LCL) is one of the most commonly applied physical laws for flow in single fractures (SF) and fractured media. The foundation of LCL is Darcian flow. This experimental study examines if LCL is valid for flow in a single rough fracture and how the fracture roughness and Reynolds number (Re) affect flow. Similar to the Moody diagram for flow in pipes, a diagram for flow in a single rough fracture has been generated to relate the friction coefficient with Re and the roughness. Under the experimental condition of this study, flow appears to be substantially different from Darcian flow. The flow law of q ∝ enJm appears to be valid for describing the flow scheme where q, e, and J are the unit width flux, the average aperture, and the hydraulic gradient. The value of the power index m is found to be around 0·83 ～ 0·98, less than what has been used in Darcian flow (m = 1). The power index n is around 11·2 and 13·0, much greater than the n value used in the LCL (n = 3), and it increases with the average velocity. The Moody type of diagram shows that the friction factor for flow in SFs is influenced by Re and the roughness. It decreases with Re when Re is small, and becomes less sensitive to Re when Re is large enough. It also increases with the roughness. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

11.
  相似文献   

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

13.
With the notion that interface and boundary layer phenomena play an important part in those geophysical processes which, by observation appear to be related to the earth's internal boundaries between the solid and liquid phases of its core and mantle, constitutional supercooling suggests itself as a mechanism capable of generating and maintaining inhomogeneities in concentration and density at the boundaries of the liquid core. The mechanism of constitutional supercooling requires a slow overgrowth of mantle and core, and, it implies that this growth process is associated with a selective partitioning of certain impurities shared in different concentrations by the liquid core and the solid phases of mantle and inner core. It can lead to the formation of regular (quasi-periodic) corrugations of the core-mantle and the inner-outer core boundaries with amplitudes of the order of 1 km. Mass redistributions, off-setting continually regenerated concentration and density inhomogeneities, provide a mechanism for core motion in the form of concentration currents. A regular distribution of corrugations or humps may give rise to (zonal) patterns of closed loops of concentration currents either in layers adjacent to the solid-liquid interfaces, or in loops extending through the entire outer core. The development of regular flow patterns should be enhanced if, referable to one particular constituent of the liquid phase, some parts of the solid-liquid interfaces acted as sources, others as sinks.  相似文献   

14.
The concepts of homogenization theory are employed to derive a macro-scale brine transport equation for micro-heterogeneous porous medium of layered structure under assumptions of validity of classical Darcy’s law and Fick’s law at the local scale. Derived macro-scale model is analogous to the so-called phase field equations. The obtained results are verified with direct numerical experiment.  相似文献   

15.
This paper numerically investigates the characteristics of groundwater flow in spatially correlated variable aperture fractures under the mechanical effect. Spatially correlated aperture distributions are generated using the geostatistical method (i.e., turning bands algorithm in this study). To represent a nonlinear relationship between the effective normal stress and the fracture aperture, a simple mechanical formula is combined with a local flow model. Numerical results indicate that the groundwater flow is significantly affected by the geometry of aperture distribution, varying with the applied effective normal stress as well as the spatial correlation length of aperture distribution. Moreover, using the flow results simulated in this study, two empirical formulae are proposed: (1) the first one (modified Louis formula) is to represent the relationship between the effective normal stress and the effective permeability of fracture and (2) the second one is to represent the relationship between relative roughness and effective permeability.  相似文献   

16.
Acceptance limits of the structural response of walls for low‐rise concrete housing were developed. Proposed values are applicable within a performance‐based seismic design framework. Acceptance limits are based on performance indicators of structural response–allowable story drift ratios, width of residual cracks and residual damage index, and expected damage of walls. Cracking limits were defined from parameters obtained at the unloading stage of walls (i.e., residual cracking stage). The residual cracking stage may be used for structural damage evaluation and cost estimation of structural rehabilitation after an earthquake has occurred. The performance indicators proposed herein were derived from test observations and measured response of 39 RC walls' specimens during shaking table and quasistatic testing, as well as from limiting values and results of previous studies. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

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

18.
The east-west asymmetry has been analyzed using the 1946–2000 corona observation data for the green line at the Kislovodsk station. A positive east-west asymmetry has been revealed throughout the entire observation period except for three years, 1994–1996. Time variations of the east-west asymmetry for certain periods coincide with those for flares. Seasonal variations of the east-west asymmetry revealed earlier by other authors are not confirmed. If seasonal variations of the east-west asymmetry in the corona exist, their amplitudes are smaller than or comparable to the instrumental errors, errors caused by atmospheric variations during the observation period and to differences between the corona-intensity measuring systems used at different observatories.  相似文献   

19.
Fluid conductivity and elastic properties in fractures depend on the aperture geometry – in particular, the roughness of fracture surfaces. In this study, we have characterized the surface roughness with a log-normal distribution and investigated the transport and flow behaviour of the fractures with varying roughness characteristics. Numerical flow and transport simulations have been performed on a single two-dimensional fracture surface, whose aperture geometry changes with different variances and correlation lengths in each realization. We have found that conventional measurement of hydraulic conductivity alone is insufficient to determine these two parameters. Transient transport measurements, such as the particle breakthrough time, provide additional constraints to the aperture distribution. Nonetheless, a unique solution to the fracture aperture distribution is still under-determined with both hydraulic conductivity and transport measurements. From numerical simulations at different compression states, we have found that the flow and transport measurements exhibit different rates of changes with respect to changes in compression. Therefore, the fracture aperture distribution could be further constrained by considering the flow and transport properties under various compression states.  相似文献   

20.
Using a magnetic dynamo model, suggested by Kazantsev (J. Exp. Theor. Phys. 1968, vol. 26, p. 1031), we study the small-scale helicity generation in a turbulent electrically conducting fluid. We obtain the asymptotic dependencies of dynamo growth rate and magnetic correlation functions on magnetic Reynolds numbers. Special attention is devoted to the comparison of a longitudinal correlation function and a function of magnetic helicity for various conditions of asymmetric turbulent flows. We compare the analytical solutions on small scales with numerical results, calculated by an iterative algorithm on non-uniform grids. We show that the exponential growth of current helicity is simultaneous with the magnetic energy for Reynolds numbers larger than some critical value and estimate this value for various types of asymmetry.  相似文献   

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

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