共查询到20条相似文献，搜索用时 31 毫秒
1.
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. 相似文献
2.
3.
This paper treats the upscaling of the absolute permeability in a heterogeneous reservoir. By replacing the fine scale permeability tensor with an upscaled, or effective permeability tensor, a modelling error is introduced. An a posteriori error estimate on this modelling error is formulated and tested. An implementation of the theory, based on domain decomposition coupled with a hierarchical representation of the absolute permeability field, is given. As hierarchical basis functions we have chosen the Haar system, which leads to a wavelet representation of the permeability. The wavelet representation offers a natural upscaling technique which resembles the highcut filters commonly used in signal analysis. This procedure represents an adaptive upscaling method. The numerical results show that this method conserves both the dissipation and the mean velocity in the problem fairly well. The a posteriori error estimate on the modelling error coupled with domain decomposition methods constitutes a powerful modelling tool. 相似文献
4.
The impact of organic matter on the flow capacity of shale oil rocks is presumably significant, and the knowledge about the representative size is fundamental for the upscaling studies. The error of the experimentally determined permeability values is comparable with the contribution of kerogen to shale permeability, instead a 2D numerical model is employed to explore the normalised equivalent permeability and the representative elementary area (REA) of shale oil rocks in detail incorporating the effects of kerogen. The discussions on the normalised equivalent permeability and the REA are based on the statistical average and standard deviation from 1000 different runs, respectively. The inorganic permeability heterogeneity is introduced based on the assumption of a lognormal pore size distribution and the Monte Carlo sampling method. The effects of kerogen geometric characteristics are incorporated by putting forward several representative cases for comparison. The effects of the organic permeability contrast (ratio of permeability to the inorganic permeability with no heterogeneity), total organic carbon (TOC, volume fraction), inorganic permeability heterogeneity and kerogen geometric characteristics on the normalised equivalent permeability (ratio of the intrinsic equivalent permeability to inorganic permeability with no heterogeneity) and the REA are discussed comprehensively. This work can provide a better understanding of shale oil rocks at the micrometer scale. 相似文献
5.
A multi-resolution workflow to generate high-resolution models constrained to dynamic data 总被引：2，自引：0，他引：2
Céline Scheidt Jef Caers Yuguang Chen Louis J. Durlofsky 《Computational Geosciences》2011,15(3):545-563
Distance-based stochastic techniques have recently emerged in the context of ensemble modeling, in particular for history
matching, model selection and uncertainty quantification. Starting with an initial ensemble of realizations, a distance between
any two models is defined. This distance is defined such that the objective of the study is incorporated into the geological
modeling process, thereby potentially enhancing the efficacy of the overall workflow. If the intent is to create new models
that are constrained to dynamic data (history matching), the calculation of the distance requires flow simulation for each
model in the initial ensemble. This can be very time consuming, especially for high-resolution models. In this paper, we present
a multi-resolution framework for ensemble modeling. A distance-based procedure is employed, with emphasis on the rapid construction
of multiple models that have improved dynamic data conditioning. Our intent is to construct new high-resolution models constrained
to dynamic data, while performing most of the flow simulations only on upscaled models. An error modeling procedure is introduced
into the distance calculations to account for potential errors in the upscaling. Based on a few fine-scale flow simulations,
the upscaling error is estimated for each model using a clustering technique. We demonstrate the efficiency of the method
on two examples, one where the upscaling error is small, and another where the upscaling error is significant. Results show
that the error modeling procedure can accurately capture the error in upscaling, and can thus reproduce the fine-scale flow
behavior from coarse-scale simulations with sufficient accuracy (in terms of uncertainty predictions). As a consequence, an
ensemble of high-resolution models, which are constrained to dynamic data, can be obtained, but with a minimum of flow simulations
at the fine scale. 相似文献
6.
This paper is aimed towards investigating the filtration law of an incompressible viscous Newtonian fluid through a rigid non-inertial porous medium (e.g. a porous medium placed in a centrifuge basket). The filtration law is obtained by upscaling the flow equations at the pore scale. The upscaling technique is the homogenization method of multiple scale expansions which rigorously gives the macroscopic behaviour and the effective properties without any prerequisite on the form of the macroscopic equations. The derived filtration law is similar to Darcy's law, but the tensor of permeability presents the following remarkable properties: it depends upon the angular velocity of the porous matrix, it verifies Hall–Onsager's relationship and it is a non-symmetric tensor. We thus deduce that, under rotation, an isotropic porous medium leads to a non-isotropic effective permeability. In this paper, we present the results of numerical simulations of the flow through rotating porous media. This allows us to highlight the deviations of the flow due to Coriolis effects at both the microscopic scale (i.e. the pore scale), and the macroscopic scale (i.e. the sample scale). The above results confirm that for an isotropic medium, phenomenological laws already proposed in the literature fails at reproducing three-dimensional Coriolis effects in all types of pores geometry. We show that Coriolis effects may lead to significant variations of the permeability measured during centrifuge tests when the inverse Ekman number Ek−1 is 𝒪(1). These variations are estimated to be less than 5% if Ek−1<0.2, which is the case of classical geotechnical centrifuge tests. We finally conclude by showing that available experimental data from tests carried out in centrifuges are not sufficient to determining the effective tensor of permeability of rotating porous media. Copyright © 2004 John Wiley & Sons, Ltd. 相似文献
7.
We propose a methodology, called multilevel local–global (MLLG) upscaling, for generating accurate upscaled models of permeabilities
or transmissibilities for flow simulation on adapted grids in heterogeneous subsurface formations. The method generates an
initial adapted grid based on the given fine-scale reservoir heterogeneity and potential flow paths. It then applies local–global
(LG) upscaling for permeability or transmissibility [7], along with adaptivity, in an iterative manner. In each iteration of MLLG, the grid can be adapted where needed to reduce
flow solver and upscaling errors. The adaptivity is controlled with a flow-based indicator. The iterative process is continued
until consistency between the global solve on the adapted grid and the local solves is obtained. While each application of
LG upscaling is also an iterative process, this inner iteration generally takes only one or two iterations to converge. Furthermore,
the number of outer iterations is bounded above, and hence, the computational costs of this approach are low. We design a
new flow-based weighting of transmissibility values in LG upscaling that significantly improves the accuracy of LG and MLLG
over traditional local transmissibility calculations. For highly heterogeneous (e.g., channelized) systems, the integration
of grid adaptivity and LG upscaling is shown to consistently provide more accurate coarse-scale models for global flow, relative
to reference fine-scale results, than do existing upscaling techniques applied to uniform grids of similar densities. Another
attractive property of the integration of upscaling and adaptivity is that process dependency is strongly reduced, that is,
the approach computes accurate global flow results also for flows driven by boundary conditions different from the generic
boundary conditions used to compute the upscaled parameters. The method is demonstrated on Cartesian cell-based anisotropic
refinement (CCAR) grids, but it can be applied to other adaptation strategies for structured grids and extended to unstructured
grids. 相似文献
8.
9.
Uncertainty quantification is typically accomplished by simulating multiple geological realizations, which can be very expensive computationally if the flow process is complicated and the models are highly resolved. Upscaling procedures can be applied to reduce computational demands, though it is essential that the resulting coarse-model predictions correspond to reference fine-scale solutions. In this work, we develop an ensemble level upscaling (EnLU) procedure for compositional systems, which enables the efficient generation of multiple coarse models for use in uncertainty quantification. We apply a newly developed global compositional upscaling method to provide coarse-scale parameters and functions for selected realizations. This global upscaling entails transmissibility and relative permeability upscaling, along with the computation of a-factors to capture component fluxes. Additional features include near-well upscaling for all coarse parameters and functions, and iteration on the a-factors, which is shown to improve accuracy. In the EnLU framework, this global upscaling is applied for only a few selected realizations. For 90 % or more of the realizations, upscaled functions are assigned statistically based on quickly computed flow and permeability attributes. A sequential Gaussian co-simulation procedure is incorporated to provide coarse models that honor the spatial correlation structure of the upscaled properties. The resulting EnLU procedure is applied for multiple realizations of two-dimensional models, for both Gaussian and channelized permeability fields. Results demonstrate that EnLU provides P10, P50, and P90 results for phase and component production rates that are in close agreement with reference fine-scale results. Less accuracy is observed in realization-by-realization comparisons, though the models are still much more accurate than those generated using standard coarsening procedures. 相似文献
10.
Identifying the Representative Elementary Volume for Permeability in Heterolithic Deposits Using Numerical Rock Models 总被引：1，自引：0，他引：1
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. 相似文献
11.
The upscaling process of a high-resolution geostatistical reservoir model to a dynamic simulation grid model plays an important role in a reservoir study. Several upscaling methods have been proposed in order to create balance between the result accuracy and computation speed. Usually, a high-resolution grid model is upscaled according to the heterogeneities assuming single phase flow. However, during injection processes, the relative permeability adjustment is required. The so-called pseudo-relative permeability curves are accepted, if their corresponding coarse model is a good representation of the fine-grid model. In this study, an upscaling method based on discrete wavelet transform (WT) is developed for single-phase upscaling based on the multi-resolution analysis (MRA) concepts. Afterwards, an automated optimization method is used in which evolutionary genetic algorithm is applied to estimate the pseudo-relative permeability curves described with B-spline formulation. In this regard, the formulation of B-spline is modified in order to describe the relative permeability curves. The proposed procedure is evaluated in the gas injection case study from the SPE 10th comparative solution project’s data set which provides a benchmark for upscaling problems [1]. The comparisons of the wavelet-based upscaled model to the high-resolution model and uniformly coarsened model show considerable speedup relative to the fine-grid model and better accuracy relative to the uniformly coarsened model. In addition, the run time of the wavelet-based coarsened model is comparable with the run time of the uniformly upscaled model. The optimized coarse models increase the speed of simulation up to 90% while presenting similar results as fine-grid models. Besides, using two different production/injection scenarios, the superiority of WT upscaling plus relative permeability adjustment over uniform upscaling and relative permeability adjustment is presented. This study demonstrates the proposed upscaling workflow as an effective tool for a reservoir simulation study to reduce the required computational time. 相似文献
12.
Although a number of methods for calculating dynamic pseudo-functions have been developed over the years, there is still a lack of understanding as to why a certain method will succeed in some cases but fail in others. In this paper, we describe the results of an assessment of several upscaling methods, namely the Kyte and Berry (KB) method, the Stone method, the Hewett and Archer (HA) method and the Transmissibility-Weighted (TW) method. We have analyzed the equations for deriving the methods and investigated the results of numerical simulations of gas displacing oil, in a variety of models to enable us to gain new insights into these, and related, upscaling methods. In particular, some novel observations on methods based on fluid potential are presented and the issue of using predicted fluid mobilities as a criterion of accuracy of an upscaling method is clarified. 相似文献
13.
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. 相似文献
14.
Qinghua Lei Xiaoguang Wang Jiansheng Xiang John-Paul Latham 《Hydrogeology Journal》2017,25(8):2251-2262
A study about the influence of polyaxial (true-triaxial) stresses on the permeability of a three-dimensional (3D) fractured rock layer is presented. The 3D fracture system is constructed by extruding a two-dimensional (2D) outcrop pattern of a limestone bed that exhibits a ladder structure consisting of a “through-going” joint set abutted by later-stage short fractures. Geomechanical behaviour of the 3D fractured rock in response to in-situ stresses is modelled by the finite-discrete element method, which can capture the deformation of matrix blocks, variation of stress fields, reactivation of pre-existing rough fractures and propagation of new cracks. A series of numerical simulations is designed to load the fractured rock using various polyaxial in-situ stresses and the stress-dependent flow properties are further calculated. The fractured layer tends to exhibit stronger flow localisation and higher equivalent permeability as the far-field stress ratio is increased and the stress field is rotated such that fractures are preferentially oriented for shearing. The shear dilation of pre-existing fractures has dominant effects on flow localisation in the system, while the propagation of new fractures has minor impacts. The role of the overburden stress suggests that the conventional 2D analysis that neglects the effect of the out-of-plane stress (perpendicular to the bedding interface) may provide indicative approximations but not fully capture the polyaxial stress-dependent fracture network behaviour. The results of this study have important implications for understanding the heterogeneous flow of geological fluids (e.g. groundwater, petroleum) in subsurface and upscaling permeability for large-scale assessments. 相似文献
15.
油藏精细地质模型网格粗化算法及其效果 总被引：1，自引：0，他引：1
在前人研究基础上, 根据DP(Dykstra-Parsons)系数能定量评价储层非均质性, 微网格块的渗透率值粗化后, 其等效渗透率的上、下限(Cmin、Cmax)能反映渗透率的各向异性的特点, 提出了一种运算速度快和相对有效的网格粗化算法。该算法能考虑到储层非均质性对不同方向渗透率值的影响, 且求解过程相对简单。应用该方法对鄂尔多斯盆地中部某油藏的陆相储层的精细地质模型进行了网格粗化计算, 然后在粗化后的模型上进行油藏数值模拟研究, 同时针对研究区地质背景和产出流体微可压缩的物性特征, 首次利用流线模拟器对精细地质模型进行了油藏数值模拟研究, 并以此结果为标准, 对该网格粗化算法时效性进行了系统评价。分析表明, 该算法具有较快的计算速度和较高的可靠性, 是解决储层非均质强、物性差的陆相成因油藏精细油藏数值模拟的一种行之有效的手段。 相似文献
16.
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. 相似文献
17.
We review and perform comparison studies for three recent multiscale methods for solving elliptic problems in porous media
flow; the multiscale mixed finite-element method, the numerical subgrid upscaling method, and the multiscale finite-volume
method. These methods are based on a hierarchical strategy, where the global flow equations are solved on a coarsened mesh
only. However, for each method, the discrete formulation of the partial differential equations on the coarse mesh is designed
in a particular fashion to account for the impact of heterogeneous subgrid structures of the porous medium. The three multiscale
methods produce solutions that are mass conservative on the underlying fine mesh. The methods may therefore be viewed as efficient,
approximate fine-scale solvers, i.e., as an inexpensive alternative to solving the elliptic problem on the fine mesh. In addition,
the methods may be utilized as an alternative to upscaling, as they generate mass-conservative solutions on the coarse mesh.
We therefore choose to also compare the multiscale methods with a state-of-the-art upscaling method – the adaptive local–global
upscaling method, which may be viewed as a multiscale method when coupled with a mass-conservative downscaling procedure.
We investigate the properties of all four methods through a series of numerical experiments designed to reveal differences
with regard to accuracy and robustness. The numerical experiments reveal particular problems with some of the methods, and
these will be discussed in detail along with possible solutions. Next, we comment on implementational aspects and perform
a simple analysis and comparison of the computational costs associated with each of the methods. Finally, we apply the three
multiscale methods to a dynamic two-phase flow case and demonstrate that high efficiency and accurate results can be obtained
when the subgrid computations are made part of a preprocessing step and not updated, or updated infrequently, throughout the
simulation.
The research is funded by the Research Council of Norway under grant nos. 152732 and 158908. 相似文献
18.
Subsurface flow models can exhibit strong full-tensor anisotropy due to either permeability or grid nonorthogonality effects.
Upscaling procedures, for example, generate full-tensor effects on the coarse scale even for cases in which the underlying
fine-scale permeability is isotropic. A multipoint flux approximation (MPFA) is often needed to accurately simulate flow for
such systems. In this paper, we present and apply a different approach, nonlinear two-point flux approximation (NTPFA), for
modeling systems with full-tensor effects. In NTPFA, transmissibility (which provides interblock connections) is determined
from reference global flux and pressure fields for a specific flow problem. These fields can be generated using either fully
resolved or approximate global simulations. The use of fully resolved simulations leads to an NTPFA method that corresponds
to global upscaling procedures, while the use of approximate simulations gives a method corresponding to recently developed
local–global techniques. For both approaches, NTPFA algorithms applicable to both single-scale full-tensor permeability systems
and two-scale systems are described. A unified framework is introduced, which enables single-scale and two-scale problems
to be viewed in a consistent manner. Extensive numerical results demonstrate that the global and local–global NTPFA techniques
provide accurate flow predictions over wide parameter ranges for both single-scale and two-scale systems, though the global
procedure is more accurate overall. The applicability of NTPFA to the simulation of two-phase flow in upscaled models is also
demonstrated. 相似文献
19.
Homogenization has proved its effectiveness as a method of upscaling for linear problems, as they occur in single-phase porous media flow for arbitrary heterogeneous rocks. Here we extend the classical homogenization approach to nonlinear problems by considering incompressible, immiscible two-phase porous media flow. The extensions have been based on the principle of preservation of form, stating that the mathematical form of the fine-scale equations should be preserved as much as possible on the coarse scale. This principle leads to the required extensions, while making the physics underlying homogenization transparent. The method is process-independent in a way that coarse-scale results obtained for a particular reservoir can be used in any simulation, irrespective of the scenario that is simulated. Homogenization is based on steady-state flow equations with periodic boundary conditions for the capillary pressure. The resulting equations are solved numerically by two complementary finite element methods. This makes it possible to assess a posteriori error bounds. 相似文献
20.
To physically investigate permeability upscaling, over 13,000 permeability values were measured with four different sample supports (i.e., sample volumes) on a block of Berea Sandstone. At each sample support, spatially exhaustive permeability datasets were measured, subject to consistent flow geometry and boundary conditions, with a specially adapted minipermeameter test system. Here, we present and analyze a subset of the data consisting of 2304 permeability values collected from a single block face oriented normal to stratification. Results reveal a number of distinct and consistent trends (i.e., upscaling) relating changes in key summary statistics to an increasing sample support. Examples include the sample mean and semivariogram range that increase with increasing sample support and the sample variance that decreases. To help interpret the measured mean upscaling, we compared it to theoretical models that are only available for somewhat different flow geometries. The comparison suggests that the nonuniform flow imposed by the minipermeameter coupled with permeability anisotropy at the scale of the local support (i.e., smallest sample support for which data is available) are the primary controls on the measured upscaling. This work demonstrates, experimentally, that it is not always appropriate to treat the local-support permeability as an intrinsic feature of the porous medium, that is, independent of its conditions of measurement. 相似文献