共查询到20条相似文献，搜索用时 531 毫秒
1.
Jean-Laurent Mallet 《Mathematical Geology》2004,36(1):1-32
Interpolating physical properties in the subsurface is a recurrent problem in geology. In sedimentary geology, the geometry of the layers is generally known with a precision much superior to that which one can reasonably expect for the properties. The geometry of the layers is affected by folding and faulting since the time of deposition, whereas the distribution of properties is, to a certain extent, determined at the time of deposition. As a consequence, it may be wise to model first the geometry of the layers and then, simplify the geologic equation by removing the influence of that geometry. Inspired from the work of H. E. Wheeler on Time-Stratigraphy, we define, mathematically, a new space where all horizons are horizontal planes and where faults, if any, have disappeared. We surmise that this new space, however approximative, is better to model physical properties of the subsurface whatever the subsequent interpolation method used. The proposed mathematical framework also provides solutions to complex problems such as determination of strains resulting from tectonic events and up-scaling of permeabilities on structured and unstructured 3D grids. 相似文献
2.
Stochastic continuum modeling of flow and transport in a crystalline rock mass: Fanay-Augères,France, revisited 总被引：1，自引：0，他引：1
A stochastic discrete-fracture model was used by Cacas et al.a,b to interpret flow measurements and transport experiments in a fractured crystalline rock mass at Fanay-Augères. They considered continuum models to be incapable of properly interpreting small-scale measurements or tracer tests in fractured systems, which, in their view, require three-dimensional modeling of numerous discrete channels; in their opinion, continuum modeling applies only to average flow on a relatively large scale. Cacas et al. considered their discrete fracture model to have been validated by its demonstrated ability to reproduce selected experimental results. In this paper, flow and transport at Fanay-Augères are modeled by viewing the fractured rock as a stochastic continuum in a manner originally proposed by Neumanc,d. The stochastic continuum approach obviates the need for detailed information about fracture geometry or assumptions about how individual fractures control flow and transport. All it requires is the delineation of a few dominant features, which can be embedded into the stochastic continuum model as heterogeneous porous slabs. Though a fault zone has been identified at the Fanay-Augères experimental site, it has been modeled neither by Cacas et al. nor in this paper. In fact, in this paper, a larger selection of experimental results than those considered by Cacas et al. are reproduced merely by modeling the rock as a statistically homogeneous continuum in two dimensions. These results demonstrate that a continuum approach may be well suited for the analysis of flow and transport in fractured rock. This does not constitute a validation of the continuum approach, just as the results of Cacas et al. fall short of validating the discrete fracture approach. Instead, the two sets of results illustrate jointly the well-established principle that an open system, especially one as complex as fractured hydrogeologic environments tend to be, cannot be described uniquely on the basis of sparse data and need not be described in great detail to capture its salient behavior by a model.a Cacas MC, Ledoux E, de Marsily G, Barbreau A, Calmels P, Gaillard B, Margritta R (1990a) Modelling fracture flow with a stochastic discrete fracture network: calibration and validation. 1. The flow model. Water Resour Res 26(3):479–489b Cacas MC, Ledoux E, de Marsily G, Barbreau A, Calmels P, Gaillard B, Margritta R (1990b) Modelling fracture flow with a stochastic discrete fracture network: calibration and validation. 2. The transport model. Water Resour Res 26(3):491–500c Neuman SP (1987) Stochastic continuum representation of fractured rock permeability as an alternative to the REV and fracture network concepts, in Rock Mechanics. In: Farmer IW, Daemen JJK, Desai CS, Glass CE, Neuman SP (eds) Proceedings of the 28th U.S. Symposium, Tucson, Arizona. Balkema, Rotterdam, pp 533–561d Neuman SP (1988) A proposed conceptual framework and methodology for investigating flow and transport in Swedish crystalline rocks. SKB Swedish Nuclear Fuel and Waste Management Co., Stockholm, September, Arbetsrapport 88–37, 39 pp 相似文献
3.
László Szarka 《中国地质大学学报(英文版)》2009,20(3)
As I learned it from extensive geo-electromaguetic analogue modeling experiments, some specific nonconventional interpretation parameters, in certain conditions, give more detailed information about the geometry of subsurface resistivity inbomogeneities than the routinely used parameters. In this article, I show several examples, and I present how these early results influenced our later research. An enhanced geometric sensitivity may be due to special array geometry (as we call It "null array"), or it may be due to a narrow and very special frequency range (Le., the so-called "keyhole" range). Nonconventional but physically based interpretation parameters (like the Poyuting vector) or higher order invariants of resistivity or impedance tensors may also give useful additional information about the shape of subsurface bodies. One should be very careful in their application because a large part of these nontroditioual approaches are strongly constrained by measuring errors and geological noise. 相似文献
4.
Layered rock masses can be modelled either as standard, orthotropic continua if the layer bending can be neglected or as Cosserat continua if the influence of layer bending is essential. This paper presents a finite element smeared joint model based on the Cosserat theory. The layers are assumed to be elastic with equal thickness and equal mechanical properties. All the cosserat parameters are expressed through the elastic properties of layers, layer thickness and joint stiffness. Plastic-slip as well as tensile-opening of layer interface (joint) are accounted for in a manner similar to the conventional non-associative plasticity theory. As an application, the behaviour of an excavation in a layered rock mass is examined. The displacement and stress fields given by smeared joint models based on the Cosserat continuum and the conventional anisotropic continuum approaches are compared with those obtained from the discrete joint model. The conventional anisotropic continuum model is found to break-down completely when the effective shear modulus in the direction parallel to layering is low in comparison to the shear modulus of the intact layer, whereas the Cosserat model is found to be capable of accurately reproducing complex load–deflection patterns irrespective of the differences in shear moduli. © 1998 John Wiley & Sons, Ltd. 相似文献
5.
6.
四川盆地侏罗系发育具有分型结构的超致密储层,孔隙度多数在5%之下,渗透率小于0.1m D。长期以来,有关侏罗系油藏是裂缝油藏,岩性油藏,还是致密油,一直存在争议;对于如此致密储层,在没有裂缝的情况下,石油能不能进入储层,是如何进去的,运移路径是什么;是储层先致密后进油,还是先进油而后储层致密化;致密油的运聚机理与富集因素是什么,应采取怎样的勘探思路等等,是勘探和研究中需要面临和解决的基础科学问题。本文运用流体包裹体、场发射和环境扫描电镜、纳米CT和成藏物理模拟等方法,从储层微观介质空间含油性、流体包裹体分布及特征、储层演化与生烃演化等研究入手,开展了侏罗系石油运移、聚集机理和富集规律研究,结果表明四川盆地侏罗系油层以致密油为主体,储层致密化过程中及以后发生石油充注,致密储层基质中各种成岩次生微缝、结晶矿物晶间缝隙、各种微-纳米缝隙等基质缝是石油运移和渗滤有效路径,生烃增压和毛细管压力差是致密油运聚的主要动力,致密储层非达西渗流特征明显,具有启动压力梯度。四川盆地侏罗系致密油具有"源储共控、近源富集、网状赋油、甜点高产"的分布规律,指出在勘探中需要遵循致密油的自身属性和地质规律,按照"整体部署、分层评价、优选甜点"的工作思路,实现创新研究、科学评价和效益勘探。 相似文献
7.
Analysis of strength and moduli of jointed rocks 总被引：1，自引：0，他引：1
This paper deals with two aspects of jointed rock mass behavior, first the finite element modeling of a jointed rock mass as an equivalent continuum, second the comparison of empirical strength criteria of a jointed rock mass. In finite element modeling the jointed rock properties are represented by a set of empirical relationships, which express the properties of the jointed medium as a function of joint factor and the properties of the intact rock. These relationships have been derived from a large set of experimental data of tangent elastic modulus. It is concluded that equivalent continuum analysis gives the best results for both single and multiple jointed rock. The reliability of the analysis depends on the estimation of joint factor, which is a function of the joint orientation, joint frequency and joint strength.Empirical strength criteria for jointed rocks, namely Hoek and Brown, Yudhbir et al., Ramamurthy and Arora, Mohr–Coulomb have been incorporated in a nonlinear finite element analysis of jointed rock using the equivalent continuum approach, to determine the failure stress. The major principal stress at failure, obtained using Ramamurthy's criteria, compares very well with experimental results. 相似文献
8.
The effectiveness of secondary recovery methods in reservoir development studies depends on the knowledge about how fluid-carrying regions (i.e. good-quality rock types) are connected between injection and production wells. To estimate reservoir performance uncertainty, comprehensive simulations on many reservoir model realisations are necessary, which is very CPU consuming and time demanding. Alternatively, we can use much simpler and physically based methods such as percolation approach. Classic percolation assumes connectivity between opposite 2-D faces of a 3-D system; whereas, hydrocarbon production is achieved through active wells that are one-dimensional lines (e.g. vertical, horizontal or deviated wells). The main contribution of this study is to analyse the percolation properties of 3-D continuum percolation models with more realistic well representations during secondary recovery. In particular, the connection of randomly distributed sands (i.e. good-quality rock types) between two lines (representing two wells) located at two corners of the system are modelled by Monte Carlo simulations. Subsequently, the connectivity and conductivity of such a line-to-line well representation is compared with that of face-to-face well representations in the previously published results. The critical percolation properties of those systems as well as the universality concept are also investigated. As there are many rooms for connections in 3-D models, we found that the principal percolation properties will not be altered significantly when the problem with a face-to-face connection is transformed to a line-to-line connection model. 相似文献
9.
Witherspoon, P.A. and Gale, J.E., 1977. Mechanical and hydraulic properties of rocks related to induced seismicity. Eng. Geol., 11(1): 23–55.The mechanical and hydraulic properties of fractured rocks are considered with regard to the role they play in induced seismicity. In many cases, the mechanical properties of fractures determine the stability of a rock mass. The problems of sampling and testing these rock discontinuities and interpreting their non-linear behavior are reviewed. Stick slip has been proposed as the failure mechanism in earthquake events. Because of the complex interactions that are inherent in the mechanical behavior of fractured rocks, there seems to be no simple way to combine the deformation characteristics of several sets of fractures when there are significant perturbations of existing conditions. Thus, the more important fractures must be treated as individual components in the rock mass.In considering the hydraulic properties, it has been customary to treat a fracture as a parallel-plate conduit and a number of mathematical models of fracture systems have adopted this approach. Non-steady flow in fractured systems has usually been based on a two-porosity model, which assumes the primary (intergranular) porosity contributes only to storage and the secondary (fracture) porosity contributes only to the overall conductivity. Using such a model, it has been found that the time required to achieve quasi-steady state flow in a fractured reservoir is one or two orders of magnitude greater than it is in a homogeneous system. In essentially all of this work, the assumption has generally been made that the fractures are rigid.However, it is clear from a review of the mechanical and hydraulic properties that not only are fractures easily deformed but they constitute the main flow paths in many rock masses. This means that one must consider the interaction of mechanical and hydraulic effects. A considerable amount of laboratory and field data is now available that clearly demonstrates this stress-flow behavior. Two approaches have been used in attempting to numerically model such behavior: (1) continuum models, and (2) discrete models. The continuum approach only needs information as to average values of fracture spacing and material properties. But because of the inherent complexity of fractured rock masses and the corresponding decrease in symmetry, it is difficult to develop an equivalent continuum that will simulate the behavior of the entire system. The discrete approach, on the other hand, requires details of the fracture geometry and material properties of both fractures and rock matrix. The difficulty in obtaining such information has been considered a serious limitation of discrete models, but improved borehole techniques can enable one to obtain the necessary data, at least in shallow systems. The possibility of extending these methods to deeper fracture systems needs more investigation. Such data must be considered when deciding whether to use a continuum or discrete model to represent the interaction of rock and fluid forces in a fractured rock system, especially with regard to the problem of induced seismicity. When one is attempting to alter the pressure distribution in a fault zone by injection or withdrawal of fluids, the extent to which this can be achieved will be controlled in large measure by the behavior of the fractures that communicate with the borehole. Since this is essentially a point phenomenon, i.e., the changes will propagate from a relatively small region around the borehole, the use of a discrete model would appear to be preferable. 相似文献
10.
针对致密砂岩气储层复杂的微观孔隙结构进行岩石物理建模,在模型中比较了单一孔隙纵横比、双孔隙模型两种表征孔隙结构的表征方式。岩石物理正演分析表明,两种孔隙结构模型均可解释致密砂岩复杂的速度-孔隙度关系。岩石物理反演结果表明,双孔隙模型对测井横波速度的预测精度更高,说明该模型更适用于表征研究区致密砂岩的孔隙结构,反演的软孔比例参数能够反映地层中孔隙结构的非均匀分布。应用双孔隙模型计算致密砂岩地层岩石骨架的弹性模量,与Krief及Pride等传统经验公式相比,该方法考虑了岩石骨架模量与矿物基质、孔隙度和孔隙结构等微观物性因素的关系,理论上更具有严谨性。对致密砂岩骨架模量计算结果的分析表明,少量微裂隙的存在即能够显著影响致密砂岩骨架的弹性性质,同时孔隙空间中的球形孔隙是致密气的主要赋存空间。并且,通过致密砂岩骨架弹性模量,进一步计算了可用于地层评价的Biot系数等岩石物理参数。致密砂岩骨架模量的预测结果可为Gassmann流体替换理论、BISQ孔隙弹性介质理论等岩石物理方法提供关键参数。 相似文献
11.
Antoine S. Bailly 《GeoJournal》1993,31(3):247-250
After twenty years of work on the geography of representations, how is it still possible to define geography as “the science of space”, ie as direct knowledge of material reality? This conception of the discipline — based on Cartesian precepts of evidence (eg the observer's independent certainty), reductionism (ie a disaggregation into sets of simple elements), causality (ie the presupposition of a linear linkage between cause and effect) and exhaustiveness (ie the certainty that nothing essential has been omitted) — has been thrown into question by the geography of representations' holistic approach. How can our scientific practices be separated from our interior existence with its affective and emotional aspects? Is not scientific action an extension of being? Mustn't the geographer, above and beyond the observation of concrete phenomena, also understand the subtle and complex — at times random and hidden — links which unite human beings and their life-space, be it from the viewpoint of the poet, or of all those who take alternative approaches to geography? What I would like to demonstrate is (1) how in an historically and socially given environment, the individual constructs his own reality in linking together the structural, functional and symbolic; (2) how the representation of the landscape is related to our existential experience; and (3) how the imaginary and the real are connected in each place. 相似文献
12.
Seyed Erfan Saberhosseini Reza Keshavarzi Kaveh Ahangari 《Arabian Journal of Geosciences》2017,10(6):152
One of the crucial consequences of steam assisted gravity drainage (SAGD) process is abnormal reservoir uplifting under thermal steam injection, which can significantly influence the reservoir rock deformation, specifically thin bed reservoirs and causes intensive failures and fractures into the cap rock formations. A thorough understanding of the influences of rock thermo-mechanical properties on reservoir uplifting plays an important role in preventing those aforementioned failures within design and optimization process in SAGD. In addition, coupling of reservoir porous medium and flowing of specific fluid with temperature as an additional degree of freedom with initial pore pressure and in-situ stress condition, are also very challenging parts of geomechanical coupled simulation which would be clearly explained. Thus, a fully coupled thermo-poro-elastic geomechanical model with finite element codes was performed in ABAQUS to investigate the role of rock thermo-mechanical parameters on reservoir vertical uplift during steam injection. It is clearly observed that, any increase in rock thermo-mechanical properties specifically rock’s thermal properties such as specific heat, thermal expansion, and formation’s thermal conductivity, have significant influences on reservoir uplift. So by coupling the temperature as an additional degree of freedom with the coupled pore-fluid stress and diffusion finite element model of SAGD process, the more realistic simulation will be conducted; hence, the errors related to not having heat as an additional degree of freedom will be diminished. In addition, Young’s modulus and specific heat are the rock thermo-mechanical parameters which have the maximum and minimum effects on the reservoir uplift, respectively. 相似文献
13.
Fully coupled simulation of a hydraulic fracture interacting with natural fractures with a hybrid discrete‐continuum method 
下载免费PDF全文
下载免费PDF全文 A hybrid discrete‐continuum numerical scheme is developed to study the behavior of a hydraulic fracture crossing natural fractures. The fully coupled hybrid scheme utilizes a discrete element model for an inner domain, within which the hydraulic fracture propagates and interacts with natural fractures. The inner domain is embedded in an outer continuum domain that is implemented to extend the length of the hydraulic fracture and to better approximate the boundary effects. The fracture is identified to propagate initially in the viscosity‐dominated regime, and the numerical scheme is calibrated by using the theoretical plane strain hydraulic fracture solution. The simulation results for orthogonal crossing indicate three fundamental crossing scenarios, which occur for various stress ratios and friction coefficients of the natural fracture: (i) no crossing, that is, the hydraulic fracture is arrested by the natural fracture and makes a T‐shape intersection; (ii) offset crossing, that is, the hydraulic fracture crosses the natural fracture with an offset; and (iii) direct crossing, that is, the hydraulic fracture directly crosses the natural fracture without diversion. Each crossing scenario is associated with a distinct net pressure history. Additionally, the effects of strength contrast and stiffness contrast of rock materials and intersection angle between the hydraulic fracture and the natural fracture are also investigated. The simulations also illustrate that the level of fracturing complexity increases as the number and extent of the natural fractures increase. As a result, we can conclude that complex hydraulic fracture propagation patterns occur because of complicated crossing behavior during the stimulation of naturally fractured reservoirs. Copyright © 2017 John Wiley & Sons, Ltd. 相似文献
14.
What determines the geometry of spheroidal weathering patterns? Spheroidal weathering rings appear to arise from the leaching of fresh rock, rather than from reprecipitation of dissolved minerals diffusing out of weathered rock. In previous work the authors developed a mathematical model to explain the geometry of spheroidal weathering structures, and showed how some of the features of weathering patterns, otherwise difficult to explain, were natural consequences of the model. In describing the geometrical aspect of these shapes, we hope to understand how natural forces acted upon the unweathered rock and why and under what conditions spheroidal weathering structures were formed. In this paper we review the main features of our model and examine and discuss twinning of spheroidal weathering structures, concluding that the structures were formed during a distinct period in the past. Twinning is an apparent anomaly in that its natural occurrence, as part of the weathering process, cannot be explained by our model without additional ad hoc assumptions. Careful examination of twins reveals an explanation which does not require revision of the model. 相似文献
15.
Current visualization techniques for computational fluid dynamics applications are sophisticated and work well in simple geometries. For complex geometries such as pore spaces, multiple domain boundaries obstruct the view and make the studying of fluid flow fields difficult. To overcome these deficiencies, we use two-sided materials to render the domain boundaries. Using this technique, it is possible to place the camera inside the domain and have a non-obstructed view of the surrounding flow field without losing spatial reference to the domain boundaries. As a result, a larger part of fluid flow visualization is visible. Two-sided material rendering was successfully applied to display still images with Blender Cycles renderer, in a virtual reality environment, and several implementation techniques were explored for using the Visualization Toolkit. 相似文献
16.
The paper describes a method for the radio identification of decameter-wavelength sources based on their continuum spectra and analysis of their coordinates in relatively large error boxes surrounding a specified position on the sky. The distribution of continuum spectra and identifications in other wavelength ranges are analyzed for the resulting radio catalog. Using identifications with the FIRST and NVSS surveys, the statistics of the spectral index-size and spectral index-flux density distributions for steep-spectrum sources have been studied, and a catalog of ultrasteep-spectrum (α相似文献
17.
Dual-continuum (DC) models can be tractable alternatives to explicit approaches for the numerical modelling of multiscale materials with multiphysics behaviours. This work concerns the conceptual and numerical modelling of poroelastically coupled dual-scale materials such as naturally fractured rock. Apart from a few exceptions, previous poroelastic DC models have assumed isotropy of the constituents and the dual-material. Additionally, it is common to assume that only one continuum has intrinsic stiffness properties. Finally, little has been done into validating whether the DC paradigm can capture the global poroelastic behaviours of explicit numerical representations at the DC modelling scale. We address the aforementioned knowledge gaps in two steps. First, we utilise a homogenisation approach based on Levin's theorem to develop a previously derived anisotropic poroelastic constitutive model. Our development incorporates anisotropic intrinsic stiffness properties of both continua. This addition is in analogy to anisotropic fractured rock masses with stiff fractures. Second, we perform numerical modelling to test the DC model against fine-scale explicit equivalents. In doing, we present our hybrid numerical framework, as well as the conditions required for interpretation of the numerical results. The tests themselves progress from materials with isotropic to anisotropic mechanical and flow properties. The fine-scale simulations show that anisotropy can have noticeable effects on deformation and flow behaviour. However, our numerical experiments show that the DC approach can capture the global poroelastic behaviours of both isotropic and anisotropic fine-scale representations. 相似文献
18.
Numerical modelling of rock slides is a versatile approach to understand the failure mechanism and the dynamics of rock slopes. Finite element slope stability analysis of three rock slopes in Garhwal Himalaya, India has been carried out using a two dimensional plane strain approach. Two different modelling techniques have been attempted for this study. Firstly, the slope is represented as a continuum in which the effect of discontinuities is considered by reducing the properties and strength of intact rock to those of rock mass. The equivalent Mohr-Coulomb shear strength parameters of generalised Hoek-Brown (GHB) criterion and modified Mohr-Coulomb (MMC) criterion has been used for this continuum approach. Secondly, a combined continuum-interface numerical method has been attempted in which the discontinuities are represented as interface elements in between the rock walls. Two different joint shear strength models such as Barton-Bandis and Patton’s model are used for the interface elements. Shear strength reduction (SSR) analysis has been carried out using a finite element formulation provided in the PHASE2. For blocky or very blocky rock mass structure combined continuum-interface model is found to be the most suitable one, as this model is capable of simulating the actual field scenario. 相似文献
19.
20.
Analysis of contaminant transport through fractured crystalline rocks has received considerable attention, particularly with regard to subsurface nuclear waste repositories. Most of the studies have employed the dual continuum approach, with the fractures and the rock matrix as the two continuums, assuming that fractures control the overall conductivity of the rock and the porous matrix just provides storage. However, field observations of rock fractures have shown that the real situation can be very complex. Based on some recent investigations, it has been reported that the portion of the rock matrix adjacent to many open fractures is physically and chemically altered. These alterations, referred to as the fracture skin, can have different sorption and diffusion properties compared to those of the undisturbed rock matrix and this may influence the transport of solutes through such formations. In the present study, a numerical model is developed to simulate conservative solute transport in a fractured crystalline rock formation using the triple continuum approach ?? with the fracture, fracture skin and the rock matrix as the three continuums. The model is solved using a fully implicit finite difference scheme. Contaminant migration in the fractured formation with and without skin has been simulated. It is observed that contaminant penetration along the fracture is enhanced at large flow velocities. The effect of flow velocity on conservative solute transport is investigated for different fracture apertures and fracture skin thicknesses. The influence of flow velocity on contaminant transport is demonstrated to be more with change in fracture aperture than with change in skin thickness. 相似文献