共查询到20条相似文献，搜索用时 31 毫秒
1.
The fractal geometry of flow paths in natural fractures in rock and the approach to percolation 总被引：2，自引：0，他引：2
The distributions of contact areas in single, natural fractures in quartz monzonite (Stripa granite) are found to have fractal dimensions which decrease fromD=2.00 to values nearD=1.96 as stress normal to the fractures is increased from 3 MPa up to 85 MPa. The effect of stress on fluid flow is studied in the same samples. Fluid transport through a fracture depends on two properties of the fracture void space geometry. the void aperture; and the tortuosity of the flow paths, determined through the distribution of contact area. Each of these quantities change under stress and contribute to changes observed in the flow rate. A general flow law is presented which separates these different effects. The effects of tortuosity on flow are largely governed by the proximity of the flow path distribution to a percolation threshold. A fractal model of correlated continuum percolation is presented which quantitatively reproduces the flow path geometries. The fractal dimension in this model is fit to the measured fractal dimensions of the flow systems to determine how far the flow systems are above the percolation threshold. 相似文献
2.
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. 相似文献
3.
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. 相似文献
4.
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. 相似文献
5.
FRICTIONAL SLIDING OF PLAGIOCLASE GOUGE UNDER LOWER-CRUST TEMPERATURE AND RELATIVELY LOW EFFECTIVE NORMAL STRESS 
下载免费PDF全文
下载免费PDF全文 The discovery of tremors on the lower crust portion of the San Andreas Fault has attracted more attention on the mechanical properties of the lower crust in recent years, and some experimental studies have been carried out to understand the mechanical behavior. Previous experiments under effective normal stresses of 200MPa have shown that pyroxene and plagioclase mineral separated from the gabbro and their mixtures all show velocity weakening in the lower-crust temperature range, which results in unstable slip when frictional sliding is the dominant deformation mechanism. This work is to examine whether the velocity-weakening behavior of plagioclase gouge also applies to relatively lower effective normal stress. Our experiments were performed under effective normal stress of about 100MPa, with a constant confining pressure control, with pore pressure of 30MPa and temperature of 100℃ to 600℃. We found that the frictional sliding of plagioclase are basically the same with the previous results obtained under effective normal stress of 200MPa, both of which show velocity weakening over the entire temperature range. The only difference is the out-of-trend drop of constitutive parameter a at 600℃ for the lower effective normal stress of 100MPa. It is thus concluded that reducing the effective normal stress has little effect on the sliding stability of plagioclase, and the previous conclusion made for mechanical behavior of the lower crust that unstable slips are possible therein also applies to the lower effective normal stress of 100MPa. 相似文献
6.
The experimental studies done at high temperature and high pressure find that increased temperature can lead to dramatic velocity and strength reductions of most of rocks at high confining pressure[1,2]. What causes this phenomenon? Is it due to dehydrati… 相似文献
7.
Compressional- and shear-wave velocities as a function of confining stress in dry sandstones 总被引：3，自引：0，他引：3
A laboratory study was carried out to investigate the influence of confining stress on compressional- and shear-wave velocities for a set of rock samples from gas-producing sandstone reservoirs in the Cooper Basin, South Australia. The suite of samples consists of 22 consolidated sublitharenites with helium porosity ranging from 2.6% to 16.6%. We used a pulse-echo technique to measure compressional- and shear-wave velocities on dry samples (cylindrical 4.6 × 2 cm) at room temperature and at elevated confining stress (≤ 60 MPa). Compressional- and shear-wave velocities in samples increase non-linearly with confining stress. A regression equation of the form V = A ? Be?DP gives a good fit to the measured velocities with improved prediction of velocities at high confining stresses compared with equations suggested by other studies. The predicted microcrack-closure stresses of the samples show values ranging from 70 MPa to 95 MPa and insignificant correlation with porosity, permeability or clay content. There is a positive correlation between change in velocity with core porosity and permeability, but this association is weak and diminishes with increasing confining stress. Experimental results show that pore geometry, grain-contact type, and distribution and location of clay particles may be more significant than total porosity and clay content in describing the stress sensitivity of sandstones at in situ reservoir effective stress. The stress dependence of Cooper Basin sandstones is very large compared with data from other studies. The implication of our study for hydrocarbon exploration is that where the in situ reservoir effective stress is much less than the microcrack-closure stress of the reservoir rocks, the variation of reservoir effective stress could cause significant changes in velocity of the reservoir rocks. The velocity changes induced by effective stress in highly stress-sensitive rocks can be detected at sonic-log and probably surface-seismic frequencies. 相似文献
8.
9.
The stress state near the Gofukuji Fault, central Japan, was estimated by simultaneously applying deformation rate analysis (DRA), the acoustic emission (AE) method, and AE rate analysis (AERA) to drilling core samples recovered from depths of 327 and 333 m at a site close to the fault. The obtained stress state was for a strike‐slip fault. It is predicted from the estimated stresses that a tensile stress of 6.4 ± 1.1 MPa acted on the borehole wall at 333 m depth during drilling. This is comparable to the typical tensile strength of granodiorite. The tensile stress estimated at 327 m depth is less than the tensile strength. This is consistent with observations of drilling‐induced tensile fracture (DITF) at depths between 329 and 334 m. Ratios of the shear stress to normal stress (normalized shear stress) acting on the fault are calculated to be 0.4–1.0, which are comparable to friction coefficients of rocks measured in laboratories. The large normalized shear stress may represent strength recovery of the Gofukuji Fault. The impendency of earthquake occurrence on this fault is evaluated to be high from the recurrence interval of earthquakes on the fault. 相似文献
10.
Numerical simulations of variable-density flow and solute transport have been conducted to investigate dense plume migration for various configurations of 2D fracture networks. For orthogonal fractures, simulations demonstrate that dispersive mixing in fractures with small aperture does not stabilize vertical plume migration in fractures with large aperture. Simulations in non-orthogonal 2D fracture networks indicate that convection cells form and that they overlap both the porous matrix and fractures. Thus, transport rates in convection cells depend on matrix and fracture flow properties. A series of simulations in statistically equivalent networks of fractures with irregular orientation show that the migration of a dense plume is highly sensitive to the geometry of the network. If fractures in a random network are connected equidistantly to the solute source, few equidistantly distributed fractures favor density-driven transport. On the other hand, numerous fractures have a stabilizing effect, especially if diffusive transport rates are high. A sensitivity analysis for a network with few equidistantly distributed fractures shows that low fracture aperture, low matrix permeability and high matrix porosity impede density-driven transport because these parameters reduce groundwater flow velocities in both the matrix and the fractures. Enhanced molecular diffusion slows down density-driven transport because it favors solute diffusion from the fractures into the low-permeability porous matrix where groundwater velocities are smaller. For the configurations tested, variable-density flow and solute transport are most sensitive to the permeability and porosity of the matrix, which are properties that can be determined more accurately than the geometry and hydraulic properties of the fracture network, which have a smaller impact on density-driven transport. 相似文献
11.
— 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. 相似文献
12.
碳是影响岩石电导率大小的一个重要因素,可能是造成龙门山断层带电导率异常的重要原因之一.为了研究不同的碳含量、矿物颗粒粒径与碳晶体结构对断层带电导率的影响,在干燥、常温、0.2~300 MPa的压力条件下实验研究了人工模拟断层泥样品(石英粉末与含碳粉末混合的样品,简称模拟样品)和采自映秀-北川断层八角庙剖面的天然断层岩样品(简称天然样品)的电导率.实验结果显示,当模拟样品中的含碳粉末连通时,电导率与碳体积率的关系符合逾渗理论模型;而含碳粉末未连通时,电导率随总孔隙度降低而指数性升高.同时模拟样品的电导率也随石英颗粒粒径的变化而发生改变.相比于模拟样品中的含碳粉末主要分布于石英颗粒支撑的孔隙中,天然样品中的碳则主要以碳膜的形式赋存在颗粒边缘,导致碳体积率相同的条件下,模拟样品的电导率小于天然样品.此外,天然样品的电导率(9×10~(-4)S·m~(-1))也要小于野外大地电磁探测的结果(0.03~0.1 S·m~(-1)).在今后的实验中还需要考虑在动态摩擦条件下对含有完整含碳结构的天然样品进行电导率的实验研究. 相似文献
13.
Understanding the Geometry of Connected Fracture Flow with Multiperiod Oscillatory Hydraulic Tests 下载免费PDF全文
下载免费PDF全文 An understanding of the spatial and hydraulic properties of fast preferential flow pathways in the subsurface is necessary in applications ranging from contaminant fate and transport modeling to design of energy extraction systems. One method for the characterization of fracture properties over interwellbore scales is Multiperiod Oscillatory Hydraulic (MOH) testing, in which the aquifer response to oscillatory pressure stimulations is observed. MOH tests were conducted on isolated intervals of wells in siliciclastic and carbonate aquifers in southern Wisconsin. The goal was to characterize the spatial properties of discrete fractures over interwellbore scales. MOH tests were conducted on two discrete fractured intervals intersecting two boreholes at one field site, and a nest of three piezometers at another field site. Fracture diffusivity estimates were obtained using analytical solutions that relate diffusivity to observed phase lag and amplitude decay. In addition, MOH tests were used to investigate the spatial extent of flow using different conceptual models of fracture geometry. Results indicated that fracture geometry at both field sites can be approximated by permeable two‐dimensional fracture planes, oriented near‐horizontally at one site, and near‐vertically at the other. The technique used on MOH field data to characterize fracture geometry shows promise in revealing fracture network characteristics important to groundwater flow and transport. 相似文献
14.
15.
Ultrasonic compressional‐ and shear‐wave velocities have been measured on 34 samples of sandstones from hydrocarbon reservoirs. The sandstones are all of low clay content, high porosity, and cover a wide range of permeabilities. They were measured dry and brine‐saturated under hydrostatic effective stresses of 10, 20 and 40 MPa. For eight of the sandstones, ultrasonic velocity measurements were made at different partial water saturations in the range from dry to fully saturated. The Gassmann–Biot theory is found to account for most of the changes in velocities at high effective stress levels when the dry sandstones are fully saturated with brine, provided the lower velocities resulting when the dry sandstone initially adsorbs small amounts of moisture are used to determine the elastic properties of the ‘dry’ sandstone. At lower effective stress levels, local flow phenomena due to the presence of open microcracks are assumed to be responsible for measured velocities higher than those predicted by the theory. The partial saturation results are modelled fairly closely by the Gassmann–Biot theory, assuming heterogeneous saturation for P‐waves. 相似文献
16.
— Knowledge of rock properties controlling the fluid movement is a basic prerequisite to understand the dynamical processes and the temperature and stress regime of the upper crust. Fracture networks were investigated on different scales to obtain quantitative results of fracture geometry like fracture length, orientation and fracture frequencies. Due to the scale effect, these parameters differ in several orders of magnitude in dependence of the scale of investigation. On the microscopic scale, fluorescent thin sections from cores were analysed and permeability was estimated for 2-D hydraulic networks. On the macroscopic scale, fracture parameters were determined from images of structural borehole measurements. The vicinity of the drill site represents the megascopic scale, where seismic reflectors were assumed as active hydraulic structures for construction of a fracture network. Compiling the fracture densities from all investigated scales and taking into consideration only the networks above the percolation threshold, the fracture length distribution follows a power law with an exponent of ?1.9 ± 0.05. Besides the scale differences of the geometric parameters like fracture density and length and the hydraulic parameters like permeability, the connectivity of the networks seems to be a confining characteristic. This is quantitatively described by the percolation parameter and the mean number of intersections per fracture. When assuming a macroscopic hydraulic system at the percolation threshold for the KTB site, the macroscopic mean fracture length can be estimated to approximately 30 m. This stands in agreement with the hydraulic experiments on site. 相似文献
17.
Cracks play a very important role in many geotechnical issues and in a number of processes in the Earth’s crust. Elastic waves can be used as a remote sensing tool for determining crack density. The effect of varying crack density in crystalline rock on the P- and S-wave velocity and dynamic elastic properties under confining pressure has been quantified. The evolution of P- and S-wave velocity were monitored as a suite of dry Westerly granite samples were taken to 60, 70, 80 and 90 % of the unconfined uniaxial strength of the sample. The damaged samples were then subjected to hydrostatic confining pressure from 2 MPa to 200 MPa to quantify the effect of varying crack density on the P- and S-wave velocity and elastic properties under confining pressure. The opening and propagation of microcracks predominantly parallel to the loading direction during uniaxial loading caused a 0.5 and 6.3 % decrease in the P- and S-wave velocity, respectively. During hydrostatic loading, microcracks are closed at 130 MPa confining pressure. At lower pressures the amount of crack damage in the samples has a small but measureable effect. We observed a systematic 6 and 4 % reduction in P- and S-wave velocity, respectively, due to an increase in the fracture density at 2 MPa confining pressure. The overall reduction in the P- and S-wave velocity decreased to 2 and 1 %, respectively, at 50 MPa. The elastic wave velocities of samples that have a greater amount of microcrack damage are more sensitive to pressure. Effective medium modelling was used to invert elastic wave velocities and infer crack density evolution. Comparing the crack density results with experimental data on Westerly granite samples shows that the effective medium modelling used gave interpretable and reasonable results. Changes in crack density can be interpreted as closure or opening of cracks and crack growth. 相似文献
18.
Poro-acoustoelastic theory has made a great progress in both theoretical and experimental aspects, but with no publications on the joint research from theoretical analyses, experimental measurements, and numerical validations. Several key issues challenge the joint research with comparisons of experimental and numerical results, such as digital imaging of heterogeneous poroelastic properties, estimation of acoustoelastic constants, numerical dispersion at high frequencies and strong heterogeneities, elastic nonlinearity due to compliant pores, and contamination by boundary reflections. Conventional poroacoustoelastic theory, valid for the linear elastic deformation of rock grains and stiff pores, is modified by incorporating a dualporosity model to account for elastic nonlinearity due to compliant pores subject to high-magnitude loading stresses. A modified finite-element method is employed to simulate the subtle effect of microstructures on wave propagation in prestressed digital cores. We measure the heterogeneity of samples by extracting the autocorrelation length of digital cores for a rough estimation of scattering intensity. We conductexperimental measurements with a fluid-saturated sandstone sample under a constant confining pressure of 65 MPa and increasing pore pressures from 5 to 60 MPa. Numerical simulations for ultrasound propagation in the prestressed fluid-saturated digital core of the sample are followed based on the proposed poro-acoustoelastic model with compliant pores. The results demonstrate a general agreement between experimental and numerical waveforms for different stresses, validating the performance of the presented modeling scheme. The excellent agreement between experimental and numerical coda quality factors demonstrates the applicability for the numerical investigation of the stress-associated scattering attenuation in prestressed porous rocks. 相似文献
19.
Tracer experiments conducted in the laboratory on undisturbed core samples (<7.3-cm-diameter) have been a standard method for estimating hydraulic and transport properties of fractured till since the 1980s. This study assesses the relationship between visible fractures on the top and bottom of core samples and the resulting hydraulic and mass transport properties of the core. We hypothesized that more visible fractures would indicate the presence of a well-connected fracture network, leading to greater hydraulic conductivity (K) values and earlier chemical breakthrough times. To test this hypothesis, water flow and bromide (Br-) tracer experiments were performed on 10, 16-cm diameter, 16-cm-tall samples of fractured Dows Formation till from central Iowa. Visually identifiable fractures were present on the top and bottom of every sample. Results indicate that the visual identification of fractures does not predict a connected fracture network, as some samples produced breakthrough curves showing rapid first arrival times and shapes characteristic of solute transport in a fractured medium, while others appeared similar to an unfractured medium. No correlation was found between the number of visible fractures and K (Pearson's r = 0.25), or Br- first arrival time (r = −0.33), but a strong negative correlation between K and first arrival time (r = −0.92). Results indicate that the sample volume was not large enough to reliably contain a connected fracture network. Thus, testing large volumes of till at the field scale coupled with fracture-flow modeling likely represents the best approach for estimating hydraulic and mass transport properties for fractured till. 相似文献
20.
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. 相似文献