首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
逾渗与岩石破裂的计算机模拟研究   总被引:4,自引:0,他引:4  
用计算机模拟方法研究了岩石破裂的逾渗模型,给出了在微裂隙随机分布情况下破裂集团的分布图象、破裂集团的大小分布规律、总破裂集团数、集团平均大小、b值、逾渗分维D及临界情形。所得结果与其他模型和实验基本一致,表明通渗模型是描述岩石破裂扩展演化过程的一种较合适的模型。还发现了岩石破裂过程存在的两个转折点,深入研究这两个转折点附近岩石的宏观性质,可为探索地震孕育演化过程的前兆性质提供依据。  相似文献   

2.
Cavitation,a phenomenon produced by a moving fluid,is ubiquitous in the water environment of the Earth's surface and its related mechanical action in the process of cavitation leads to the widespread erosion of rock in nature.Although the mechanical action of flowing water body that accelerates the rock mass loss and fragmentation of rock(abrasion,erosion,and etching)and other phenomena have been much studied,its acceleration of mineral crystal dissolution is rarely reported.The physical mechanism of effect is not yet clear.The cavitation bubble produced in the cavitation process is at the micron level,and its related mechanical action leading to the accumulation of rock mineral dissolution is manifested in time and space in the process of the chemical element's migration between water and rock minerals.Cavitation erosion may be one of the important driving forces for the migration of geochemical elements within the lithosphere and hydrosphere.In this paper,based on the crystal dissolution stepwave dynamic theory and the theoretical derivation and calculation of Gibbs free energy change of the mineral crystals plastic deformation which caused by the mechanical action of cavitation erosion,we give the possible mechanism of accelerating the transient dissolution of mineral crystals by cavitation erosion—the cavitation bubbles on the surface of the near crystal release the high speed micro-jet and shock wave perpendicular to the surface during the collapsing,in which the water hammer pressure produced by micro-jet at the solid–liquid interface causes instantaneous plastic deformation on the crystal surface under the condition that it is larger than the yield stress of the crystal.Under the influence of the thermal effect of the plastic deformation process and the change of Gibbs free energy(the dislocation elastic strain energy of plastic deformation on the crystal surface may be included),the local instantaneous dissolution rate of the mineral surface is accelerated.The continuous cavitation erosion eventually causes fracture and breaking of the mineral crystal,meanwhile,the Gibbs–Thomson effect may enhance the dissolution of mineral crystals more prominently.At the same time,the correctness of the mechanism is verified qualitatively by the acoustic cavitation experiment with the same erosion mechanism.  相似文献   

3.
4.
Measurements of seismic anisotropy in fractured rock are used at present to deduce information about the fracture orientation and the spatial distribution of fracture intensity. Analysis of the data is based upon equivalent-medium theories that describe the elastic response of a rock containing cracks or fractures in the long-wavelength limit. Conventional models assume frequency independence and cannot distinguish between microcracks and macrofractures. The latter, however, control the fluid flow in many subsurface reservoirs. Therefore, the fracture size is essential information for reservoir engineers. In this study we apply a new equivalent-medium theory that models frequency-dependent anisotropy and is sensitive to the length scale of fractures. The model considers velocity dispersion and attenuation due to a squirt-flow mechanism at two different scales: the grain scale (microcracks and equant matrix porosity) and formation-scale fractures. The theory is first tested and calibrated against published laboratory data. Then we present the analysis and modelling of frequency-dependent shear-wave splitting in multicomponent VSP data from a tight gas reservoir. We invert for fracture density and fracture size from the frequency dependence of the time delay between split shear waves. The derived fracture length matches independent observations from borehole data.  相似文献   

5.
The solutional origin of limestone caves was recognized over a century ago, but the short penetration length of an undersaturated solution made it seem impossible for long conduits to develop. This is contradicted by field observations, where extended conduits, sometimes several kilometers long, are found in karst environments. However, a sharp drop in the dissolution rate of CaCO3 near saturation provides a mechanism for much deeper penetration of reactant. The notion of a “kinetic trigger” – a sudden change in rate constant over a narrow concentration range – has become a widely accepted paradigm in speleogenesis modeling. However, it is based on one-dimensional models for the fluid and solute transport inside the fracture, assuming that the dissolution front is planar in the direction perpendicular to the flow. Here we show that this assumption is incorrect; a planar dissolution front in an entirely uniform fracture is unstable to infinitesimal perturbations and inevitably breaks up into highly localized regions of dissolution. This provides an alternative mechanism for cave formation, even in the absence of a kinetic trigger. Our results suggest that there is an inherent wavelength to the erosion pattern in dissolving fractures, which depends on the reaction rate and flow rate, but is independent of the initial roughness. In contrast to one-dimensional models, two-dimensional simulations indicate that there is only a weak dependence of the breakthrough time on kinetic order; localization of the flow tends to keep the undersaturation in the dissolution front above the threshold for non-linear kinetics.  相似文献   

6.
相对于常规砂岩,致密砂岩在岩石物理性质、力学性质等方面具有明显差异,并呈现出很强的非均质性.岩石物理模型能将储层参数与地震特性信息联系起来,因此可以作为致密砂岩储层参数与地震特性信息转换的桥梁.常规的岩石物理模型通常只考虑单一因素(例如非均匀性,单一孔隙,单一尺度等),建立的岩石物理模板并不适用于致密砂岩.本文针对高饱和气、微裂隙发育、非均质性强的致密砂岩储层,利用Voigt-Reuss-Hill模型计算混合矿物的弹性模量,采用微分等效介质(DEM)模型描述含裂隙、孔隙岩石的骨架弹性模量,基于Biot-Rayleigh波动方程构建了岩石物理弹性模板,给出了致密砂岩储层弹性参数与物性的关系.基于测井数据和实验数据对岩石物理弹性模板进行校正,并将校正后的岩石物理弹性模板结合叠前地震资料应用于川西地区储层孔隙度与裂隙含量预测.结果显示,反演裂隙含量、孔隙度与储层试气报告、测井孔隙度基本吻合,表明该模板能够较合理地应用于致密砂岩储层孔隙度及裂隙含量解释中.  相似文献   

7.
Underground fractures play an important role in the storage and movement of hydrocarbon fluid. Fracture rock physics has been the useful bridge between fracture parameters and seismic response. In this paper, we aim to use seismic data to predict subsurface fractures based on rock physics. We begin with the construction of fracture rock physics model. Using the model, we may estimate P-wave velocity, S-wave velocity and fracture rock physics parameters. Then we derive a new approximate formula for the analysis of the relationship between fracture rock physics parameters and seismic response, and we also propose the method which uses seismic data to invert the elastic and rock physics parameters of fractured rock. We end with the method verification, which includes using well-logging data to confirm the reliability of fracture rock physics effective model and utilizing real seismic data to validate the applicability of the inversion method. Tests show that the fracture rock physics effective model may be used to estimate velocities and fracture rock physics parameters reliably, and the inversion method is resultful even when the seismic data is added with random noise. Real data test also indicates the inversion method can be applied into the estimation of the elastic and fracture weaknesses parameters in the target area.  相似文献   

8.
A systematic investigation of the effect of configurations of stochastically distributed fracture networks on hydraulic behaviour for fractured rock masses could provide either quantitative or qualitative correlation between the structural configuration of the fracture network and its corresponding hydraulic behaviour, and enhance our understanding of appropriate application of groundwater flow and contaminant transport modelling in fractured rock masses. In this study, the effect of block sizes, intersection angles of fracture sets, standard deviations of fracture orientation, and fracture densities on directional block hydraulic conductivity and representative elementary volume is systematically investigated in two dimensions by implementing a numerical discrete fracture fluid flow model and incorporating stochastically distributed fracture configurations. It is shown from this investigation that the configuration of a stochastically distributed fracture network has a significant quantitative or qualitative effect on the hydraulic behaviour of fractured rock masses. Compared with the deterministic fracture configurations that have been extensively dealt with in a previous study, this investigation is expected to be more practical and adequate, since fracture geometry parameters are inherently stochastically distributed in the field. Moreover, the methodology and approach presented in this study may be generally applied to any fracture system in investigating the hydraulic behaviours from configurations of the fracture system while establishing a ‘bridge’ from the discrete fracture network flow modelling to equivalent continuum modelling in fractured rock masses. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

9.
岩石分形断裂的统计理论   总被引:9,自引:0,他引:9       下载免费PDF全文
本文研究了岩石分形断裂的统计理论.认为岩石断裂是岩石内部微裂缝成核、形成和传播的非平衡过程,导出了描述这个过程的微分方程.运用位错塞积成核机理,得到了岩石发生分形断裂的微裂缝大小的统计分布函数和平均长大速率.从理论上得出了岩石非弹性体积应变与微裂缝数目成正比,从而解释了实验室研究的结果.另外发现岩石的分维是决定岩石脆性的一个重要因素,分维值愈小,岩石愈脆.但岩石的分维值与岩石强度无明显关系.最后,从微裂缝大小的统计分布函数出发,导出了b值与岩石分维值的关系.  相似文献   

10.
根据裂缝双侧向测井响应的特征,裂缝分3种状态:低角度裂缝、倾斜裂缝和高角度裂缝,给出了判断裂缝状态的方法;建立了裂缝的双侧向测井响应与裂缝孔隙度、裂缝内流体电导率、基岩电导率之间的函数关系;得出人工估算和计算机精确计算裂缝孔隙度的两套反演公式及相应计算方法.  相似文献   

11.
基于贝叶斯线性AVAZ的TTI介质裂缝参数反演   总被引:2,自引:0,他引:2       下载免费PDF全文
裂缝储层岩石物理参数的准确获得对地下裂缝预测具有重要意义,而叠前方位AVA地震反演是获得裂缝岩石物理参数的有效手段.假设地下岩石为倾斜横向各向同性(TTI)介质,本文从裂缝岩石物理等效模型的构建出发,从测井数据中估计出纵横波相对反射系数和裂缝柔度参数.通过推导含裂缝柔度的方位各向异性反射系数公式,基于贝叶斯反演框架建立了P波线性AVAZ反演方法.合成地震数据应用表明基于贝叶斯理论的TTI介质裂缝柔度反演方法具有一定抗噪性,可以降低裂缝柔度估测的不确定性,为地下裂缝预测提供有力的依据.  相似文献   

12.

Uncertainty in input fracture geometric parameters during analysis of the stability of jointed rock slopes is inevitable and therefore the stochastic discrete fracture network (DFN) — distinct element method (DEM) is an efficient modeling tool. In this research, potentially unstable conditions are detected in the right abutment of the Karun 4 dam and downstream of the dam body as a case study. Two extreme states with small and relatively large block sizes are selected and a series of numerical DEM models are generated using a number of validated DFN models. Stability of the rock slope is assessed in both static and dynamic loading states. Based on the design basis earthquake (DBE) and maximum credible earthquake (MCE) expected in the dam site, histories of seismic waves are applied to analyze the stability of the slope in dynamic earthquake conditions. The results indicate that a MCE is likely to trigger sliding of rock blocks on the rock slope major joint. Furthermore, the dynamic analysis also shows a local block failure by the DBE, which can consequently lead to slope instability over the long term. According to the seismic behavior of the two models, larger blocks are prone to greater instability and are less safe against earthquakes.

  相似文献   

13.
The stochastic continuum (SC) representation is one common approach for simulating the effects of fracture heterogeneity in groundwater flow and transport models. These SC reservoir models are generally developed using geostatistical methods (e.g., kriging or sequential simulation) that rely on the model semivariogram to describe the spatial variability of each continuum. Although a number of strategies for sampling spatial distributions have been published in the literature, little attention has been paid to the optimization of sampling in resource- or access-limited environments. Here we present a strategy for estimating the minimum sample spacing needed to define the spatial distribution of fractures on a vertical outcrop of basalt, located in the Box Canyon, east Snake River Plain, Idaho. We used fracture maps of similar basalts from the published literature to test experimentally the effects of different sample spacings on the resulting semivariogram model. Our final field sampling strategy was based on the lowest sample density that reproduced the semivariogram of the exhaustively sampled fracture map. Application of the derived sampling strategy to an outcrop in our field area gave excellent results, and illustrates the utility of this type of sample optimization. The method will work for developing a sampling plan for any intensive property, provided prior information for a similar domain is available; for example, fracture maps or ortho-rectified photographs from analogous rock types could be used to plan for sampling of a fractured rock outcrop.  相似文献   

14.
A field demonstration was performed at Edwards Air Force Base to assess bioaugmentation for treatment of a well‐characterized tetrachloroethene (PCE) dense nonaqueous phase liquid (DNAPL) source area in fractured rock. Groundwater recirculation was employed to deliver remedial amendments, including bacteria, to facilitate reductive dechlorination and enhance DNAPL dissolution. An active treatment period of 9 months was followed by a 10‐month posttreatment rebound evaluation. Dechlorination daughter products were observed in both the shallow and deep fracture zones following treatment. In the shallow fracture zone, the calculated DNAPL mass removed was approximately equal to the DNAPL mass estimated using partitioning tracer testing, and no rebound in chlorinated ethenes or ethene was observed during the posttreatment period. A maximum DNAPL dissolution enhancement factor of 5 was observed in the shallow fracture zone. In the deep fracture zone, only approximately 45% of the DNAPL mass—as estimated via partitioning tracer testing—was removed and rebound in the total molar chlorinated ethenes + ethene was observed. The difference in behavior between the shallow and deep fracture zones was attributed to DNAPL architecture and the fracture flow field.  相似文献   

15.
—?The potential for large excavation-induced seismic events may be recognised, even if the timing of an event may be inherently unpredictable. In this case, modelling the wave propagation from a potential event could allow the dynamic motions around an excavation to be projected, and for areas of danger to be anticipated. However, the above and other potential applications require accurate models of wave interaction with the openings, as well as with the fractured rock which surrounds such excavations. This paper considers real recorded waveforms and how well these waveforms are modelled by explicit mechanical models of the source, the medium and the excavation. Models of experiments at three different scales of the problem are presented: small and large amplitude waveforms recorded around a deep-level mining tunnel in a synthetic rockburst experiment; waveforms from laboratory experiments of waves through plates of steel representing fractures; waveforms from active pulses in an acoustic emission experiment in a small volume of fractured rock at the surface of an underground excavation. The results show that elastic wave propagation around an excavation was a first approximation for small amplitude waves, but was less successful for modelling large amplitude waves and more fractured rock. Fractures in the models were represented explicitly with displacement discontinuities. Waveforms through known fracture geometries were particularly well-reproduced, and indicate the importance of fracture stiffness, the in situ stress state, and stress-dependence of the fractures in such models. Overall, the models are sufficiently successful at representing recorded behaviour, to be encouraging for the goal of representing accurate wave motions around excavations.  相似文献   

16.
—For small-scale microseismic events, the source sizes provided by shear models are unrealistically large when compared to visual observations of rock fractures near underground openings. A detailed analysis of the energy components in data from a mine-by experiment and from some mines showed that there is a depletion of S-wave energy for events close to the excavations, indicating that tensile cracking is the dominant mechanism in these microseismic events.¶In the present study, a method is proposed to estimate the fracture size from microseismic measurements. The method assumes tensile cracking as the dominant fracture mechanism for brittle rocks under compressive loads and relates the fracture size to the measured microseismic energy. With the proposed method, more meaningful physical fracture sizes can be obtained and this is demonstrated by an example on data from an underground excavation with detailed, high-quality microseismic records.  相似文献   

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

18.
According to the requirement of the project “Establishment of the Physical Model of Earthquake Precursor Fields”, this paper elucidates the train of thinking for research on the project and some scientific problems which must be studied; the elucidation emphasizes that the core of this project is to study the conditions and processes of the generation of strong earthquakes. The paper first outlines the origin and development of the “strong-body earthquake-generating model” proposed by the author in the 1980’s; and then proves the reasonableness of the model from three aspects, namely: deep structures, mechanical analysis and rock fracture experiments. By studying the tomographic image for the northern part of North China, it can be seen that the sources of strong earthquakes are all distributed in high-velocity bodies, or in the contact zone between high-velocity and low-velocity bodies but nearer to the high-velocity body. It has been affirmed through studies of the mechanical models of hard and soft inclusions that the existence of a hard inclusion is an important condition for the high concentration of large amounts of strain energy. A lot of theoretical and experimental studies have been made to investigate the conditions for rock instability; the results have consistently indicated that rock instability, sudden fracture and stress drop would be possible only if the stiffness of the source body is greater than the environmental stiffness.  相似文献   

19.
We present a sequence of purely advective transport models that demonstrate the influence of small-scale geometric inhomogeneities on contaminant transport in fractured crystalline rock. Special weight is placed on the role of statistically generated variable fracture apertures. The fracture network geometry and the aperture distribution are based on information from an in situ radionuclide retardation experiment performed at Grimsel test site (Swiss Alps). The obtained breakthrough curves are fitted with the advection dispersion equation and continuous-time random walks (CTRW). CTRW is found to provide superior fits to the late-arrival tailing and is also found to show a good correlation with the velocity distributions obtained from the hydraulic models. The impact of small-scale heterogeneities, both in fracture geometry and aperture, on transport is shown to be considerable.  相似文献   

20.
A conceptual model for describing effective saturation in fractured hard rock is presented. The fracture network and the rock matrix are considered as an equivalent continuum medium where each fracture is conceptualized as a porous medium of granular structure and the rock matrix is assumed to be impermeable. The proposed model is based on the representation of a rough‐walled fracture by an equivalent porous medium, which is described using classical constitutive models. A simple closed‐form equation for the effective saturation is obtained when the van Genuchten model is used to describe saturation inside fractures and fractal laws are assumed for both aperture and number of fractures. The relative hydraulic conductivity for the fractured rock is predicted from a simple relation derived by Liu and Bodvarsson. The proposed constitutive model contains three independent parameters, which may be obtained by fitting the proposed effective saturation curve to experimental data. Two of the model parameters have physical meaning and can be identified with the reciprocal of the air entry pressure values in the fractures of minimum and maximum apertures. Effective saturation and relative hydraulic conductivity curves match fairly well the simulated constitutive relations obtained by Liu and Bodvarsson. Copyright © 2008 John Wiley & Sons, Ltd.  相似文献   

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

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