首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
Based on a simulation of three-dimensional fracture networks and a superposition principle of liquid dissipation energy for fractured rock masses, a model of the fracture permeability tensor is proposed. An elastic constitutive model of rock fractures, considering fracture closure and dilation during shearing, is also proposed, based on the dilation angle of the fracture. Algorithms of flow-path searching and calculation of the effective flow coefficients for fracture networks are presented, together with a discussion on the influence of geometric parameters of the fractures (trace length, spacing, aperture, orientation and the number of fracture sets) on magnitude, anisotropy of hydraulic permeability and the size of a representative elementary volume (REV). The anisotropy of hydraulic permeability of fractured rock masses is mainly affected by orientation and the number of fracture sets, and the REV size is mainly influenced by trace length, spacing and the number of fracture sets. The results of studies on REV size and the influence of in-situ stress on hydraulic conductivity of the rock mass on the slope of Jinping-I hydropower station, China, are presented using the developed models and methods. The simulation results agreed well with the results obtained from field water-pressure measurements, with an error of less than 10 %.  相似文献   

2.
Using a range of realistic 3D numerical lithofacies (dm-scale) models of ripple laminated sandstone intercalated with mudstone we evaluate how single-phase permeability varies as a function of sample support. The models represent a range of mudstone content which is typical for tidal deposits. Furthermore, the spatial distribution of flow barriers (i.e. mudstone) is not random, but governed by sedimentological rules giving a variable anisotropy ratio as a function of mudstone content. Both vertical and horizontal permeability are found to vary at small sample volumes, but these fluctuations reduce as the sample volume increases. The vertical permeability increases while the horizontal permeability is nearly constant as a function of sample support for small mudstone contents. For higher mudstone content, the horizontal permeability decreases while the vertical permeability is nearly constant as a function of sample support. We propose a criterion, based on a normalised standard deviation, to determine the Representative Elementary Volume (REV). The size of the REV is dependent on both the property measured (vertical and horizontal permeability) and the correlation lengths of the lithological elements (i.e. lithofacies). Based on this we identify three flow upscaling regimes that each require a different method for upscaling: (1) layered systems where the arithmetic and harmonic averages are appropriate, (2) systems close to the percolation threshold where a percolation model should be used, and (3) discontinuous systems where an effective medium method provides the best estimate of permeability. The work gives, by using numerical experiments on a range of heterogeneous systems, a new insight in determination of the REV for permeability at the lithofacies scale and its relation to sedimentological parameters.  相似文献   

3.
Discontinuities such as fractures or cracks are common in geo-materials. Connectivity is a comprehensive parameter that includes the influence of length, orientation and density of the discontinuities. The main objectives of this paper are to define connectivity based on the statistical geometrical parameters of a network of random discontinuities and to quantify the hydraulic properties of the network using connectivity. An analytical method is first proposed to evaluate connectivity through the density of degree of freedom (DDOF). The analytically estimated DDOF is verified using numerical analysis. DDOF increases with the increases in the mean length, relative orientation and density of the random discontinuity network. Connectivity becomes better and the permeability becomes larger with increasing DDOF. The representative element volume (REV) of a random discontinuity network can be estimated based on connectivity. When a discontinuity network is poorly connected, a REV does not exist when the size of the discontinuity network is smaller than 10 times of the mean length of the discontinuities. When a discontinuity network is normally connected, a REV exists when the REV size is smaller than 10 times of the mean length of the discontinuities. When the discontinuity network is highly connected, the REV size is only 2–4 times of the mean length of the discontinuities. The results indicate that DDOF is a good indicator to quantify the connectivity of a random discontinuity network and the relationship between connectivity and permeability.  相似文献   

4.
Katagiri  Jun  Kimura  Sho  Noda  Shohei 《Acta Geotechnica》2020,15(8):2195-2203
Acta Geotechnica - This paper is a study of determination of representative elementary volume (REV) size suitable for pore-scale flow simulation (PFS) and evaluation of permeability anisotropy for...  相似文献   

5.
Grand Mesa is an erosional remnant on the southern margin of the Piceance basin (Colorado, USA) that appears to host topographically driven groundwater flow in low permeability strata via a pervasive network of vertical extensional fractures. The vertical fractures cut more than 1 km of clay-rich lithology ranging in age from Upper Cretaceous through Eocene, and likely formed from horizontal dilation, cooling, and erosional unloading associated with 2.8 km of regional uplift and 1.5 km of incision by the Gunnison and Colorado rivers. The vertical fractures create anisotropy in which vertical permeability exceeds horizontal permeability. This enhances vertical flow and depth of penetration of groundwater, favors local flow regimes over regional flow, and results in groundwater temperatures that are elevated by up to 10°C over mean surface temperatures at the location of springs. The uplift and cooling that formed the fractures may also have produced domains of abnormally underpressured pore fluids and natural gas within blocks of low permeability bedrock bounded by the fractures. Pore pressures inside these blocks may be in disequilibrium with the groundwater flow system due to ongoing stress release, and the long time scales required for pressure equilibration in the low permeability strata.  相似文献   

6.
Three different types of permeability tests were conducted on 23 intact and singly jointed rock specimens, which were cored from rock blocks collected from a rock cavern under construction in Singapore. The studied rock types belong to inter-bedded meta-sandstone and meta-siltstone with very low porosity and high uniaxial compressive strength. The transient pulse water flow method was employed to measure the permeability of intact meta-sandstone under a confining pressure up to 30 MPa. It showed that the magnitude order of meta-sandstone’s intrinsic permeability is about 10?18 m2. The steady-state gas flow method was used to measure the permeability of both intact meta-siltstone and meta-sandstone in a triaxial cell under different confining pressures spanning from 2.5 to 10 MPa. The measured permeability of both rock types ranged from 10?21 to 10?20 m2. The influence of a single natural joint on the permeability of both rock types was studied by using the steady-state water flow method under different confining pressures spanning from 1.25 to 5.0 MPa, including loading and unloading phases. The measured permeability of both jointed rocks ranged from 10?13 to 10?11 m2, where the permeability of jointed meta-siltstone was usually slightly lower than that of jointed meta-sandstone. The permeability of jointed rocks decreases with increasing confining pressure, which can be well fitted by an empirical power law relationship between the permeability and confining pressure or effective pressure. The permeability of partly open cracked specimens is lower than that of open cracked specimens, but it is higher than that of the specimen with a dominant vein for the meta-sandstone under the same confining pressure. The permeability of open cracked rock specimens will partially recover during the unloading confining pressure process. The equivalent crack (joint) aperture is as narrow as a magnitude order of 10?6 m (1 μm) in the rock specimens under confining pressures spanning from 1.25 to 5.0 MPa, which represent the typical ground stress conditions in the cavern. The in situ hydraulic conductivity measurements conducted in six boreholes by the injection test showed that the in situ permeability of rock mass varies between 10?18 and 10?11 m2. The lower bound of the in situ permeability is larger than that of the present laboratory-tested intact rock specimens, while the upper bound of the in situ permeability is less than that of the present laboratory-tested jointed rock specimens. The in situ permeability test results were thus compatible with our present laboratory permeability results of both intact and jointed rock specimens.  相似文献   

7.
申林方  王志良  李邵军 《岩土力学》2015,36(Z2):689-694
为了研究土体的细观渗流特性,假设土体是完全饱和且在渗流过程中水分的流动始终处于层流状态。考虑宏观统计参数（孔隙率、渗透率及有效黏滞系数等）的影响,基于表征体元（REV）尺度的格子博尔兹曼（Boltzmann）方法,建立了压力作用下土体细观渗流的数值模型。采用D2Q9模型考虑水分流动的离散速度分布,宏观边界条件为左右侧面为不透水边界 ,上下边界设置不同的密度来控制压力边界,在微观边界条件上采用非平衡态外推格式。编制相应的计算程序,将计算区域内的多孔介质材料设置成流体（孔隙率 1.0,渗透率 ）,验证了经典的Poiseuille流。此外,结合算例分别讨论了土体在压力作用下孔隙率、渗透率及渗透压力等影响因素与渗流速度的相互关系,研究表明该数值方法与Darcy定律得到的计算结果较为吻合。因此,基于REV尺度的格子Boltzmann方法可以有效地模拟土体的渗流机制,为进一步研究土体渗流特性提供了一种新的研究手段。  相似文献   

8.
Pumping tests carried out in the fissured layer of a granitic hard-rock aquifer, interpreted at the observation wells by means of the analytical solution of Neuman and at the pumping wells with that of Gringarten show the existence of a strong vertical anisotropy of this layer of the aquifer; the horizontal permeability is clearly and systematically higher than the vertical one. These results agree perfectly with the geological observations, the fissured layer of the weathered granite profile showing the existence of many sub-horizontal fissures. It confirms that, within the fissured layer, the permeability of sub-horizontal fissures due to the weathering process dominates over that of sub-vertical fissures of tectonic origin. To cite this article: J.-C. Maréchal et al., C. R. Geoscience 335 (2003).To cite this article: J.-C. Maréchal et al., C. R. Geoscience 335 (2003).  相似文献   

9.
采用渗流力学理论并结合Monte Carlo方法描述岩体裂隙的随机分布,研究渗流模型的尺寸效应并确定表征单元体积（REV）,得到了3种开口度分布形式的等效渗透系数椭圆曲线,建立了等效渗透系数方向性的判别标准。离散裂隙网络（DFN）模型假定流体只在岩体裂隙内部流动,而不通过岩体本身渗流。基于二维离散元程序UDEC并进行二次开发,建立DFN模型,通过改变流体的流动方向,得到不同流动方向下岩体裂隙网络的等效渗透系数,并分析不同的开口度分布形式对岩体裂隙网络等效渗透系数方向性的影响。计算结果表明,表征单元体积存在的条件是等效渗透系数保持稳定且渗透椭圆比较光滑。等效渗透系数的方向性受开口度分布形式的影响很大：当开口度-长度关联分布时,等效渗透系数各向异性;当开口度对数正态分布时,等效渗透系数各向同性;当开口度恒定分布时,等效渗透系数的特性介于二者之间。变化系数（CV）是否大于5%是判定岩体裂隙网络渗透系数是否具有方向性的判别标准。  相似文献   

10.
Grain boundary processes contribute significantly to electronic and ionic transports in materials within Earth’s interior. We report a novel experimental study of grain boundary conductivity in highly strained olivine aggregates that demonstrates the importance of misorientation angle between adjacent grains on aggregate transport properties. We performed electrical conductivity measurements of melt-free polycrystalline olivine (Fo90) samples that had been previously deformed at 1200 °C and 0.3 GPa to shear strains up to γ?=?7.3. The electrical conductivity and anisotropy were measured at 2.8 GPa over the temperature range 700–1400 °C. We observed that (1) the electrical conductivity of samples with a small grain size (3–6 µm) and strong crystallographic preferred orientation produced by dynamic recrystallization during large-strain shear deformation is a factor of 10 or more larger than that measured on coarse-grained samples, (2) the sample deformed to the highest strain is the most conductive even though it does not have the smallest grain size, and (3) conductivity is up to a factor of ~?4 larger in the direction of shear than normal to the shear plane. Based on these results combined with electrical conductivity data for coarse-grained, polycrystalline olivine and for single crystals, we propose that the electrical conductivity of our fine-grained samples is dominated by grain boundary paths. In addition, the electrical anisotropy results from preferential alignment of higher-conductivity grain boundaries associated with the development of a strong crystallographic preferred orientation of the grains.  相似文献   

11.
Shear-wave splitting measurements from local and teleseismic earthquakes are used to investigate the seismic anisotropy in the upper mantle beneath the Rwenzori region of the East African Rift system. At most stations, shear-wave splitting parameters obtained from individual earthquakes exhibit only minor variations with backazimuth. We therefore employ a joint inversion of SKS waveforms to derive hypothetical one-layer parameters. The corresponding fast polarizations are generally rift parallel and the average delay time is about 1 s. Shear phases from local events within the crust are characterized by an average delay time of 0.04 s. Delay times from local mantle earthquakes are in the range of 0.2 s. This observation suggests that the dominant source region for seismic anisotropy beneath the rift is located within the mantle. We use finite-frequency waveform modeling to test different models of anisotropy within the lithosphere/asthenosphere system of the rift. The results show that the rift-parallel fast polarizations are consistent with horizontal transverse isotropy (HTI anisotropy) caused by rift-parallel magmatic intrusions or lenses located within the lithospheric mantle—as it would be expected during the early stages of continental rifting. Furthermore, the short-scale spatial variations in the fast polarizations observed in the southern part of the study area can be explained by effects due to sedimentary basins of low isotropic velocity in combination with a shift in the orientation of anisotropic fabrics in the upper mantle. A uniform anisotropic layer in relation to large-scale asthenospheric mantle flow is less consistent with the observed splitting parameters.  相似文献   

12.
  相似文献   

13.
For the safe disposal of high-level radioactive waste, different host rocks are currently being considered. The favorable properties of clay are low permeability, some retention capacity concerning radionuclides, and the ability to self-seal cracks and fissures, e.g. by swelling or time-dependent compaction creep. In Switzerland, the Jurassic Opalinus Clay is envisaged as a potential host rock which—at Mont Terri—is subdivided into the sandy, shaly, and carbonate-rich facies, the latter being less abundant. For long-term safety assessments, the understanding of the relations of properties (e.g. mineralogical composition and microstructure) and performance (e.g. mechanical behavior) of clays and claystones is essential. In the case of the sandy Opalinus Clay, the mechanical strength increases with increasing carbonate content, because carbonates form the matrix. The mineralogical investigation of a set of sandy facies samples proved a significantly larger carbonate content (20–40 mass %) when compared to the shaly facies (10–20 mass %). The carbonates of the shaly Opalinus Clay, on the other hand, are mostly localized fossils aligned parallel to the bedding, acting as predetermined breaking points. Image analysis of SEM images of polished sections proved the determined microstructural differences. In addition, carbonate particles of the sandy facies are mostly isometric, whereas carbonate particles of the shaly facies cover a greater range of shapes. The mechanical tests were accompanied by investigations of the p- and s-wave velocities, which revealed that the anisotropy of the sandy facies is less pronounced than sedimentological analyses would suggest. The mechanical strength, which, for the first time, presents results of real triaxial tests of the sandy facies. The samples of the sandy facies exhibit a failure strength of σ eff,B, approximately twice as high as was found for the shaly facies considering the deformation axis parallel to the bedding. Similar values were obtained when measuring perpendicularly to the bedding.  相似文献   

14.
基于三维裂隙网络的裂隙岩体表征单元体研究   总被引：1，自引：0，他引：1  
安玉华  王清 《岩土力学》2012,33(12):3775-3780
岩体表征单元体（REV）的研究一直是岩土工程中的重要问题。裂隙岩体REV是其内部复杂结构在尺寸效应上的体现,归因于岩体内裂隙系统的随机性。在三维网络模型的基础上通过计算体积节理数建立了确定岩体REV的方法。这种方法考虑了岩体内裂隙的三维空间密度的特征,采用体积节理数这一参数的收敛状况确定岩体REV的大小。建立了实际岩体的模拟裂隙系统,用一定大小的立方体来分割裂隙系统,统计不同立方体内的体积节理数,并将均值作为最终裂隙岩体的参数。将此参数作为裂隙岩体的属性,通过判别其收敛状况获得了岩体REV的尺寸。考虑了岩体不连续面的特征,建立了REV的分析模型,得出工程岩体的REV的尺寸是裂隙平均迹长的4倍的结论。  相似文献   

15.
夏露  刘晓非  于青春 《岩土力学》2010,31(12):3991-3996
表征单元体（REV,或典型单元体）的存在是应用连续介质方法研究岩体问题的前提,因此,裂隙岩体表征单元体研究是岩体力学及岩体水力学中的一个基础问题,表征单元体的存在与否可以从不同物理参数的角度进行分析。从岩体块体化程度的角度对其存在性进行了较系统的分析,把岩体中由裂隙圈闭形成的孤立块体的总体积在岩体中所占体积百分比定义为裂隙岩体的块体化程度。介绍了GeneralBlock软件的基本分析过程,利用软件建立了三维裂隙岩体模型。根据国际岩石力学学会的裂隙分级,针对5个延展性级别、7个间距级别建立了35种岩体结构模型。利用GeneralBlock对这35种岩体进行了块体识别,分别讨论了这35种岩体结构的块体化程度随着模型范围变化的波动情况。由于所建立的岩体结构模型是随机的,对每种结构的每个研究范围进行了5次随机实现。结果表明,对所有的35种岩体,表征单元体的尺寸在4～8倍间距之间,不超过8倍间距。  相似文献   

16.
A hydrodynamic threshold between Darcian and non-Darcian flow conditions was found to occur in cubes of Key Largo Limestone from Florida, USA (one cube measuring 0.2 m on each side, the other 0.3 m) at an effective porosity of 33% and a hydraulic conductivity of 10 m/day. Below these values, flow was laminar and could be described as Darcian. Above these values, hydraulic conductivity increased greatly and flow was non-laminar. Reynolds numbers (Re) for these experiments ranged from <0.1 to 7. Non-laminar flow conditions observed in the hydraulic conductivity tests were observed at Re close to 1. Hydraulic conductivity was measured on all three axes in a permeameter designed specifically for samples of these sizes. Positive identification of vertical and horizontal axes as well as 100% recovery for each sample was achieved. Total porosity was determined by a drying and weighing method, while effective porosity was determined by a submersion method. Bulk density, total porosity and effective porosity of the Key Largo Limestone cubes averaged 1.5 g/cm3, 40 and 30%, respectively. Two regions of anisotropy were observed, one close to the ground surface, where vertical flow dominated, and the other associated with a dense-laminar layer, below which horizontal flow dominated.  相似文献   

17.
ABSTRACT

Strong seismic anisotropy is observed in many subduction zones. This effect is attributed partly to subducting oceanic crust that is transformed into blueschist facies rocks. Because blueschist facies constituents such as glaucophane, epidote, and phengite show strong anisotropic elasticity, seismic anisotropy in subducting oceanic crust can be attributed to the lattice preferred orientation (LPO) of these minerals. We studied the deformation fabrics and seismic properties of phengite-rich epidote–glaucophane schists from the Franciscan Complex of Ring Mountain, California. The samples are composed mainly of glaucophane, epidote, and phengite. Some samples contain abundant phengite, the maximum being 40%. The LPOs of glaucophane showed that the [001] axes are aligned subparallel to lineation, and both (110) poles and [100] axes are aligned subnormal to foliation. The epidote [001] axes are aligned subnormal to foliation, with both (110) and (010) poles aligned subparallel to lineation. The LPOs of phengite are characterized by the maxima of [001] axes subnormal to foliation, and both (110) and (010) poles and [100] axes are aligned in a girdle subparallel to foliation. The phengite showed substantially strong seismic anisotropy (AVP = 42%, max.AVS = 37%). The glaucophane schist with abundant phengite showed significantly stronger seismic anisotropy (AVP = 30%, max.AVS = 23%) than the epidote–glaucophane schist (AVP = 13%, max.AVS = 9%). When the subduction angle of phengite-rich glaucophane schist is considered, the polarization direction of the fast S-waves for vertically propagating S-waves changed to a nearly trench-parallel direction for the subduction angle of 45?60°, and the S-wave anisotropy became stronger for vertically propagating S-waves with increasing subduction angles. Our data showed that phengite-rich blueschist facies rock can therefore contribute to the strong trench-parallel seismic anisotropy occurring at the subducting oceanic crust and at the slab–mantle interface in many subduction zones.  相似文献   

18.
低延展性裂隙岩体REV存在性研究   总被引：1，自引：0，他引：1  
张莉丽  王云  张辛  王允  于青春 《现代地质》2013,27(3):743-749
采用自行开发的离散元软件FractureToKarst,根据国际岩石力学协会对岩石的分类,讨论了低延展性裂隙岩体典型单元体(REV)的存在性及大小,并给出确定其存在性的一般方法。根据研究结果可知,低延展性裂隙岩体REV存在的条件为裂隙平均间距在0.6 m以下,即极密间距、很密间距、密间距和中等间距裂隙岩体存在REV;裂隙平均间距大于0.6 m,即宽间距、很宽间距和极宽间距裂隙岩体不存在REV。 关键词：裂隙岩体;渗流;典型单元体(REV);离散元  相似文献   

19.
We consider flow and upscaling of flow properties from pore scale to Darcy scale, when the pore-scale geometry is changing. The idea is to avoid having to solve for the pore evolution at the pore scale, because this results in unmanageable complexity. We propose to use stochastic modeling to parametrize plausible modifications of the pore geometry and to construct distributions of permeability parametrized by Darcy-scale variables. To localize the effects of, e.g., clogging, we introduce an intermediate scale of pore-network models. We use local Stokes solvers to calibrate the throat permeability.  相似文献   

20.
A numerical procedure to determine the equivalent permeability tensor of a fractured rock is presented, using a stochastic REV (Representative Elementary Volume) concept that uses multiple realizations of stochastic DFN (Discrete Fracture Network) models. Ten square DFN models are generated using the Monte Carlo simulations of the fracture system based on the data obtained from a site characterization program at Sellafield, Cumbria, UK. Smaller models with varying sizes of from 0.25 m×0.25 m to 10 m×10 m are extracted from the generated DFN models and are used as two-dimensional geometrical models for calculation of equivalent permeability tensor. The DFN models are also rotated in 30º intervals to evaluate the tensor characteristics of calculated directional permeability. Results show that the variance of the calculated permeability values decreases significantly as the side lengths of the DFN models increase, which justifies the existence of a REV. The REV side length found in this analysis is about 5 m and 8 m with 20% and 10% acceptable variations, respectively. The calculated directional permeability values at the REV size have tensor characteristic that is confirmed by a close approximation of an ellipse in a polar plot of the reciprocal of square roots of the directional permeability.
Ki-Bok MinEmail: Phone: +46-8-7907919Fax: +46-8-7906810
  相似文献   

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

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