共查询到20条相似文献，搜索用时 15 毫秒
1.
Sedimentary rocks have structures on all length scales from the millimeter to the kilometer. These structures are generally associated with variations in rock permeability. These need to be modeled if we are to make predictions about fluid flow through the rock. However, existing computers are not powerful enough for us to be able to represent all scales of heterogeneity explicitly in our fluid flow models—hence, we need to upscale. Small cell renormalization is a fast method for upscaling permeability, derived from an analogue circuit of resistors. However, it assumes that the small scale permeability distribution is known. In practice, this is unlikely. The only information available about small scale properties is either qualitative, derived from the depositional setting of the reservoir, or local to the wells as a result of coring or logging. The influence of small scale uncertainty on large scale properties is usually modelled by the Monte Carlo method. This is time-consuming and inaccurate if not enough realisations are used. This paper describes a new implementation of renormalization, which enables the direct upscaling of uncertain small-scale permeabilities to produce the statistical properties of the equivalent coarse grid. This is achieved by using a perturbation expansion of the resistor-derived equation. The method is verified by comparison with numerical simulations using the Monte Carlo method. The prediction of expected large-scale permeability and its standard deviation are shown to be accurate for small cell standard deviations of up to 40% of the mean cell value, using just the first nonzero term of the perturbation expansion. Inclusion of higher order terms allows larger standard deviations to be modeled accurately. Evaluation of cross-terms allows correlations of actual cell values, over and above the background structure of mean cell values. The perturbation method is significantly faster than conventional Monte Carlo simulation. It needs just two calculations whereas the Monte Carlo method needs many thousands of realisations to be generated and renormalized to converge. This results in significant savings in computer time. 相似文献
2.
Certain rock properties that depend on intergranular fracture and frictional sliding appear to be independent of rock type. This relationship is true for the rock-on-rock frictional sliding coefficient. The generalization has been widely applied to geomechanical modelling of upper crustal strength. Porous sandstones can be relatively weak and poorly cohesive, hence susceptible to deformation involving grain fragmentation and pore collapse. The critical state theory is commonly applied to describe such behaviour. Previous work showed that the yield surface is substantially independent of rock type when mean stress and differential stress are normalized by the grain crushing pressure, implying that the critical state line is rock type-independent and equivalent to the frictional sliding criterion. We test these hypotheses using previously published data for a range of porous sandstones augmented by new experimental results on Hollington and Berea sandstones deformed to large strains to define the critical state line over a wide range of pressures for each rock type. Results confirm the rock type-independence of the critical state line and show that it is nearly equivalent to frictional sliding. These relationships point to a simple procedure for estimating approximately the mechanical properties of sandstones based only on petrographic characteristics. 相似文献
3.
Relations of Rock Structure and Composition to Petrophysical and Geomechanical Rock Properties: Examples from Permocarboniferous Red-Beds 总被引：2，自引：0，他引：2
Summary. This study relates textural properties to physical and mechanical properties of coarse grained sedimentary rocks of Permocarboniferous age. As an equivalent to rock texture the principle of geomechanical order is applied. The geomechanical order describes a rock as a function of its structural and compositional order which are derived from petrological analyses. Our results indicate that rock properties like density and porosity are stronger dependent on the structural order, while strength properties additionally depend on the compositional order. The ultrasonic wave velocity responds to both structural and compositional properties. These evidences imply that the geomechanical order is not an independent parameter but a variable function of structural or compositional features, which needs specification for correlation purposes to distinct physical and mechanical rock properties. 相似文献
4.
砂泥岩间互地层等效岩石力学参数计算模型及其应用 总被引：4，自引：3，他引：1
以地层岩石的变形分析为基础,从能量的角度出发,利用应变能理论及能量守恒定律推导出了等效岩石力学参数的一般计算公式,并建立了测井资料计算模型;利用该模型计算了库车坳陷克拉A气田目的层系巴什基奇克组的等效岩石力学参数,并以计算结果为基础进行了应力场及储层裂缝的数值模拟;最后讨论了等效岩石力学参数测井资料计算模型的优点及适用条件。模拟结果显示,储层裂缝集中发育在背斜高点及单斜构造部位,与气田目前的开发实践基本一致,表明了所建模型的合理性;等效岩石力学参数适用于区域性的横向应力场及储层裂缝预测,不适于预测局部或垂向上的储层裂缝分布规律。 相似文献
5.
Uncertainty quantification for geomechanical and reservoir predictions is in general a computationally intensive problem, especially if a direct Monte Carlo approach with large numbers of full-physics simulations is used. A common solution to this problem, well-known for the fluid flow simulations, is the adoption of surrogate modeling approximating the physical behavior with respect to variations in uncertain parameters. The objective of this work is the quantification of such uncertainty both within geomechanical predictions and fluid-flow predictions using a specific surrogate modeling technique, which is based on a functional approach. The methodology realizes an approximation of full-physics simulated outputs that are varying in time and space when uncertainty parameters are changed, particularly important for the prediction of uncertainty in vertical displacement resulting from geomechanical modeling. The developed methodology has been applied both to a subsidence uncertainty quantification example and to a real reservoir forecast risk assessment. The surrogate quality obtained with these applications confirms that the proposed method makes it possible to perform reliable time–space varying dependent risk assessment with a low computational cost, provided the uncertainty space is low-dimensional. 相似文献
6.
Mortar Upscaling for Multiphase Flow in Porous Media 总被引：1，自引：0，他引：1
In mortar space upscaling methods, a reservoir is decomposed into a series of subdomains (blocks) in which independently constructed numerical grids and possibly different physical models and discretization techniques can be employed in each block. Physically meaningful matching conditions are imposed on block interfaces in a numerically stable and accurate way using mortar finite element spaces. Coarse mortar grids and fine subdomain grids provide two-scale approximations. In the resulting effective solution flow is computed in subdomains on the fine scale while fluxes are matched on the coarse scale. In addition the flexibility to vary adaptively the number of interface degrees of freedom leads to more accurate multiscale approximations. This methodology has been implemented in the Center for Subsurface Modeling's multiphysics multiblock simulator IPARS (Integrated Parallel Accurate reservoir Simulator). Computational experiments demonstrate that this approach is scalable in parallel and it can be applied to non-matching grids across the interface, multinumerics and multiphysics models, and mortar adaptivity. Moreover unlike most upscaling approaches the underlying systems can be treated fully implicitly. 相似文献
7.
Extension of plasticity theory to debonding,grain dissolution,and chemical damage of calcarenites 
下载免费PDF全文
下载免费PDF全文 The mechanical properties of calcarenites are known to be significantly affected by water saturation: both stiffness and strength decrease for wetting in the short term and for chemical dissolution in the long term. Both processes mainly affect bonds among grains: immediately after inundation depositional bonds fall in suspension, whereas diagenetic bonds dissolve more slowly. In this paper, the authors started from the micro‐structural analysis of the weathering processes to conceive a strain hardening hydro‐chemo‐mechanical coupled elastoplastic constitutive model. The concept of extended hardening rules is here enriched: weathering functions have been determined by employing a micro to macro simplified upscaling procedure. Chemical damage is incorporated into the formulation by means of a scalar damage function. Its evolution is also described by using a multiscale approach. A new term is added to the strain rate tensor in order to incorporate the dissolution induced chemical deformations developing once the soft rock is turned into a granular material. A calibration procedure for the constitutive parameters is suggested, and the model is validated by using both coupled and uncoupled chemo‐mechanical experimental test results. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献
8.
Modern geostatistical techniques allow the generation of high-resolution heterogeneous models of hydraulic conductivity containing millions to billions of cells. Selective upscaling is a numerical approach for the change of scale of fine-scale hydraulic conductivity models into coarser scale models that are suitable for numerical simulations of groundwater flow and mass transport. Selective upscaling uses an elastic gridding technique to selectively determine the geometry of the coarse grid by an iterative procedure. The geometry of the coarse grid is built so that the variances of flow velocities within the coarse blocks are minimum. Selective upscaling is able to handle complex geological formations and flow patterns, and provides full hydraulic conductivity tensor for each block. Selective upscaling is applied to a cross-bedded formation in which the fine-scale hydraulic conductivities are full tensors with principal directions not parallel to the statistical anisotropy of their spatial distribution. Mass transport results from three coarse-scale models constructed by different upscaling techniques are compared to the fine-scale results for different flow conditions. Selective upscaling provides coarse grids in which mass transport simulation is in good agreement with the fine-scale simulations, and consistently superior to simulations on traditional regular (equal-sized) grids or elastic grids built without accounting for flow velocities. 相似文献
9.
The aim of upscaling is to determine equivalent homogeneous parameters at a coarse-scale from a spatially oscillating fine-scale parameter distribution. To be able to use a limited number of relatively large grid-blocks in numerical oil reservoir simulators or groundwater models, upscaling of the permeability is frequently applied. The spatial fine-scale permeability distribution is generally obtained from geological and geostatistical models. After upscaling, the coarse-scale permeabilities are incorporated in the relatively large grid-blocks of the numerical model. If the porous rock may be approximated as a periodic medium, upscaling can be performed by the method of homogenization. In this paper the homogenization is performed numerically, which gives rise to an approximation error. The complementarity between two different numerical methods – the conformal-nodal finite element method and the mixed-hybrid finite element method – has been used to quantify this error. These two methods yield respectively upper and lower bounds for the eigenvalues of the coarse-scale permeability tensor. Results of 3D numerical experiments are shown, both for the far field and around wells. 相似文献
10.
Large‐scale poroelastic fractured reservoirs modeling using the fast multipole displacement discontinuity method 下载免费PDF全文
下载免费PDF全文 An effective approach to modeling the geomechanical behavior of the network and its permeability variation is to use a poroelastic displacement discontinuity method (DDM). However, the approach becomes rather computationally intensive for an extensive system of cracks, particularly when considering coupled diffusion/deformation processes. This is because of additional unknowns and the need for time‐marching schemes for the numerical integration. The Fast Multipole Method (FMM) is a technique that can accelerate the solution of large fracture problems with linear complexity with the number of unknowns both in memory and CPU time. Previous works combining DDM and FMM for large‐scale problems have accounted only for elastic rocks, neglecting the fluid leak‐off from the fractures into the matrix and its influence on pore pressure and stress field. In this work we develop an efficient geomechanical model for large‐scale natural fracture networks in poroelastic reservoirs with fracture flow in response to injection and production operations. Accuracy and computational performance of the proposed method with those of conventional poroelastic DDM are compared through several case studies involving up to several tens of thousands of boundary elements. The results show the effectiveness of the FMM approach to successfully evaluate field‐scale problems for the design of exploitation strategies in unconventional geothermal and petroleum reservoirs. An example considering faults reveals the impact of reservoir compartmentalization because of sealing faults for both geomechanical and flow variables under elastic and poroelastic rocks. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献
11.
Numerical computations of rock dissolution and geomechanical effects for CO2 geological storage 下载免费PDF全文
下载免费PDF全文 The paper is motivated by the long‐term safety analysis of the CO2 geological storage. We present a methodology for the assessment of the geomechanical impact of progressive rock dissolution. The method is based on the use of X‐ray tomography and the numerical dissolution technique. The influence of evolution of the microstructure on the macroscopic properties of the rock is analysed by using periodic homogenization method. The numerical computations show progressive degradation of all components of the stiffness (orthotropic) tensor. Moreover, the evolution of associated mass transfer properties (as tortuosity and conductivity tensors), by using the periodic homogenization method, is also calculated. The correlation between the mechanical parameters and the transfer properties during the dissolution process is presented. The results show that the highest increase of the hydraulic conductivity (in direction Y) is not associated with the highest decrease of Young modulus in this direction. Moreover, the highest decrease of Young modulus (in the direction X) is not associated with percolation in this direction. Finally, an incremental law to calculate settlement, in case of a rock with evolving microstructure, is proposed. The solution of the macroscopic settlement problem under constant stress and drained conditions showed that the geomechanical effects of the rock dissolution are rather limited. Copyright © 2014 John Wiley & Sons, Ltd. 相似文献
12.
Massimiliano Ferronato Giuseppe Gambolati Carlo Janna Pietro Teatini 《国际地质力学数值与分析法杂志》2008,32(6):633-657
The stress variation induced by gas/oil production may activate pre‐existing regional faults. This may enhance the expected land subsidence due to the generation of mechanically weak points close to the producing field. A class of elasto‐plastic interface elements (IE), specifically designed to address the mechanical behaviour of faults over a regional scale, is integrated into a finite element (FE) geomechanical model and used to investigate the role exerted by active faults in anthropogenic land subsidence. The importance of regional faults depends on a variety of factors including depth of the depleted reservoir, fault number, orientation and size, geomechanical properties of porous medium, pore pressure drawdown induced by fluid production, etc. With the aid of some representative examples, a useful indication is provided as to where and how fault activation may influence both magnitude and extent of the land subsidence bowl above producing gas/oil reservoirs, pointing to a generally limited impact on the ground surface. The simulation of a real faulted gas reservoir in a complex 3‐D setting shows that the proposed IE can be simply and efficiently incorporated into a FE geomechanical model, thus improving the quality of the stress and displacement prediction. Copyright © 2007 John Wiley & Sons, Ltd. 相似文献
13.
混凝土喷层支护节理岩体等效力学模型及其应用 总被引：1，自引：0，他引：1
考虑喷层在节理面局部的限剪和限裂作用,提出喷层支护节理岩体流变模型,在此基础上推导出喷层支护节理岩体的本构关系,建立了一种喷层支护节理岩体等效力学模型。采用该模型进行结构有限元计算,喷层、节理均不需要离散网格,进行等效模拟;当实际节理裂隙的数量与方向或喷层的厚度变化时,也不需要变化计算网格,从而可降低前处理工作量,为大型复杂结构的计算分析提供较大便利。最后,通过算例和实际工程验证了该模型的有效性--不仅具有离散模型便于模拟喷层真实作用机制的优点,而且具有等效模拟简单易行的优点,具有较好的工程应用前景。 相似文献
14.
15.
空间尺度转换是近年来区域生态水文研究领域的一个基本研究问题。其需要主要是源于模型的输入数据与所能提供的数据空间尺度不一致以及模型所代表的地表过程空间尺度与所观测的地表过程空间尺度不吻合。综述了目前区域生态水文模拟研究中常用的空间尺度转换研究方法,包括向上尺度转换和向下尺度转换。详细论述了2种向下尺度转换方法: 统计学经验模型和动态模型。前者是通过将GCM大尺度数据与长期的历史观测数据比较从而建立统计学相关模型, 然后利用这个统计学经验模型进行向下的空间尺度转换. 然而动态模型并不直接对GCM数据进行向下尺度的转换,而是对与GCM进行动态耦合的区域气候模型(RCM) 的输出数据进行空间尺度转换. 通常后者所获得的数据精度要比前者高,但是一个主要缺点就是并不是全球所有的研究区域都有对应的RCM。还详细论述了2种向上尺度转换方法: 统计学经验模型和斑块模型。前者是建立一个能代表小尺度信息在大尺度上分布的密度分布概率函数, 然后利用这个函数在所需的大尺度上进行积分而求得大尺度所需的信息。而后者是根据相似性最大化原则将大尺度划分为若干个可操作的小尺度斑块,然后将计算的每个小尺度斑块的信息平均化得到大尺度所需的信息。通常在计算这种斑块化的小尺度信息的时候,对每个小尺度也会采用统计学经验模型来计算代表整个斑块小尺度的信息。建议用斑块模型与统计学经验模型相集合的方法来实现向上的空间尺度转换 相似文献
16.
Anna Nissen Eirik Keilegavlen Tor Harald Sandve Inga Berre Jan Martin Nordbotten 《Computational Geosciences》2018,22(2):451-467
In simulation of fluid injection in fractured geothermal reservoirs, the characteristics of the physical processes are severely affected by the local occurence of connected fractures. To resolve these structurally dominated processes, there is a need to develop discretization strategies that also limit computational effort. In this paper, we present an upscaling methodology for geothermal heat transport with fractures represented explicitly in the computational grid. The heat transport is modeled by an advection-conduction equation for the temperature, and solved on a highly irregular coarse grid that preserves the fracture heterogeneity. The upscaling is based on different strategies for the advective term and the conductive term. The coarse scale advective term is constructed from sums of fine scale fluxes, whereas the coarse scale conductive term is constructed based on numerically computed basis functions. The method naturally incorporates the coupling between solution variables in the matrix and in the fractures, respectively, via the discretization. In this way, explicit transfer terms that couple fracture and matrix solution variables are avoided. Numerical results show that the upscaling methodology performs well, in particular for large upscaling ratios, and that it is applicable also to highly complex fracture networks. 相似文献
17.
A new method for upscaling fine scale permeability fields to general quadrilateral-shaped coarse cells is presented. The procedure, referred to as the conforming scale up method, applies a triangle-based finite element technique, capable of accurately resolving both the coarse cell geometry and the subgrid heterogeneity, to the solution of the local fine scale problem. An appropriate averaging of this solution provides the equivalent permeability tensor for the coarse scale quadrilateral cell. The general level of accuracy of the technique is demonstrated through application to a number of flow problems. The real strength of the conforming scale up method is demonstrated when the method is applied in conjunction with a flow-based gridding technique. In this case, the approach is shown to provide results that are significantly more accurate than those obtained using standard techniques. 相似文献
18.
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. 相似文献
19.
A Study of Stress Change and Fault Slip in Producing Gas Reservoirs Overlain by Elastic and Viscoelastic Caprocks 总被引：3，自引：2，他引：1
Geomechanical simulations were conducted to study the effects of reservoir depletion on the stability of internal and boundary faults in gas reservoirs overlain by elastic and viscoelastic salt caprocks. The numerical models were of a disk-shaped gas reservoir with idealized geometry; they mimic the structure of a gas field in the northern Netherlands which has experienced induced seismicity during gas production. The geomechanical simulations identified the area of the internal fault most sensitive to fault reactivation as coinciding with the epicenters of the largest seismic events associated with gas production. Depletion-induced shear slip is initiated at the depth of the reservoir, in the fault areas where the vertical fault throw ranges between 0.5 and 1.5 times the reservoir thickness. The extent of reactivated area differs depending on whether the caprock is viscoelastic or elastic: when it is viscoelastic, there is more down-dip shear displacement. High initial horizontal stresses in the rock salt and lower stresses in the elastic side-seal and the reservoir promote unloading of the internal and reservoir-bounding faults even before the reservoir is depleted. Particularly prone to fault reactivation are the fault zones along the interface between the reservoir rock and the salt caprock, which may already be critically stressed before depletion and are likely to be reactivated early during gas production. Stress relaxation and associated geomechanical changes affecting fault stability and ground surface deformation may continue long after production has ceased, due to the viscous behavior of the salt. 相似文献
20.
Hema J. Siriwardane Raj K. Gondle Duane H. Smith 《International Journal of Coal Geology》2009,77(1-2):188-202
Geologic sequestration in deep unmineable coal seams and enhanced coalbed methane production is a promising choice, economically and environmentally, to reduce anthropogenic gases such as carbon dioxide in the atmosphere. Unmineable coal seams are typically known to adsorb large amounts of carbon dioxide in comparison to the sizeable amounts of sorbed methane, which raises the potential for large scale sequestration projects. During the process of sequestration, carbon dioxide is injected into the coalbed and desorbed methane is produced. The coal matrix is believed to shrink when a gas is desorbed and swell when a gas is sorbed, sometimes causing profound changes in the cleat porosity and permeability of the coal seam. These changes may have significant impact on the reservoir performance. Therefore, it is necessary to understand the combined influence of swelling and shrinkage, and geomechanical properties including elastic modulus, cleat porosity, and permeability of the reservoir.The present paper deals with the influence of swelling and shrinkage on the reservoir performance, and the geomechanical response of the reservoir system during the process of geologic sequestration of carbon dioxide and enhanced coalbed methane production in an actual field project located in northern New Mexico. A three-dimensional swelling and shrinkage model was developed and implemented into an existing reservoir model to understand the influence of geomechanical parameters, as well as swelling and shrinkage properties, on the reservoir performance. Numerical results obtained from the modified simulator were compared to available measured values from that site and previous studies. Results show that swelling and shrinkage, and the combination of geomechanical and operational parameters, have a significant influence on the performance of the reservoir system. 相似文献