共查询到20条相似文献，搜索用时 15 毫秒
1.
Knut-Andreas Lie Jostein R. Natvig Stein Krogstad Yahan Yang Xiao-Hui Wu 《Computational Geosciences》2014,18(3-4):357-372
A Dirichlet–Neumann representation method was recently proposed for upscaling and simulating flow in reservoirs. The DNR method expresses coarse fluxes as linear functions of multiple pressure values along the boundary and at the center of each coarse block. The number of flux and pressure values at the boundary can be adjusted to improve the accuracy of simulation results and, in particular, to resolve important fine-scale details. Improvement over existing approaches is substantial especially for reservoirs that contain high-permeability streaks or channels. As an alternative, the multiscale mixed finite-element (MsMFE) method was designed to obtain fine-scale fluxes at the cost of solving a coarsened problem, but can also be used as upscaling methods that are flexible with respect to geometry and topology of the coarsened grid. Both methods can be expressed in mixed-hybrid form, with local stiffness matrices obtained as “inner products” of numerically computed basis functions with fine-scale sub-resolution. These basis functions are determined by solving local flow problems with piecewise linear Dirichlet boundary conditions for the DNR method and piecewise constant Neumann conditions for MsMFE. Adding discrete pressure points in the DNR method corresponds to subdividing faces in the coarse grid and hence increasing the number of basis functions in the MsMFE method. The methods show similar accuracy for 2D Cartesian cases, but the MsMFE method is more straightforward to formulate in 3D and implement for general grids. 相似文献
2.
3.
Fast 3D Reservoir Simulation and Scale Up Using Streamtubes 总被引：1，自引：0，他引：1
This paper presents an implementation of a semianalytical method for oil recovery calculation in heterogeneous reservoirs that is both fast and accurate. The method defines streamline paths based on a conventional single-phase incompressible flow calculation. By calculating the time-of-flight for a particle along a streamline and assigning a volumetric flux to each streamline, the cumulative pore volume of a streamtube containing the streamline can be calculated. Subsequently, the streamtube geometries are kept constant and the effects of the time varying mobility distribution in two-phase flow are accounted for by varying the flow rate in each streamtube, based on fluid resistance changes along the streamtube. Oil recovery calculations are then done based on the 1D analytical Buckley–Leverett solution. This concept makes the method extremely fast and easy to implement, making it ideal to simulate large reservoirs generated by geostatiscal methods. The simulation results of a 3D heterogeneous reservoir are presented and compared with those of other simulators. The results shows that the new simulator is much faster than a traditional finite difference simulator, while having the same accuracy. The method also naturally handles the upscaling of absolute and relative permeability. We make use of these upscaling abilities to generate a coarse curvilinear grid that can be used in conventional simulators with a great advantage over conventional upscaled Cartesian grids. This paper also shows an upscaling example using this technique. 相似文献
4.
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. 相似文献
5.
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. 相似文献
6.
Erick R. Burns Larry R. Bentley Rene Therrien Clayton V. Deutsch 《Hydrogeology Journal》2010,18(6):1357-1373
An upscaling algorithm has been developed that generates an irregular coarse grid that preserves flow connectivity by applying a rule-based upscaling algorithm to a fine-scale facies distribution. The algorithm is demonstrated using stochastically generated paleo-fluvial facies distributions. First, an irregular grid honoring the channel facies is created, followed by computation of effective anisotropic parameters for all coarse-grid cells. For the apparent layer-cake geometry of overbank deposits seen in outcrop, two local upscaling methods are compared: (1) the layered system approximation and (2) the mode. To assess upscaling performance, flow simulations for the original and upscaled grids are compared. The horizontal layered approximation (arithmetic mean) performs poorly, over-predicting lateral connectivity where even infrequent disconnection becomes important. Performance of the mode as an upscaling algorithm depends on the probability that a coarse-grid cell will be dominated by a single facies, and it performs surprisingly well because the upscaled grid-generation algorithm honors the channels, informing the upscaling process. Lastly, the irregular coarse grid was compared to a uniform coarse grid, showing superior performance with the irregular grid. The reduction in grid size achieved by irregular-grid generation will be a function of the geometrical complexity of the geologic objects to be honored. 相似文献
7.
Benjamin Faigle Rainer Helmig Ivar Aavatsmark Bernd Flemisch 《Computational Geosciences》2014,18(5):625-636
A sequential solution procedure is used to simulate compositional two-phase flow in porous media. We employ a multiphysics concept that adapts the numerical complexity locally according to the underlying processes to increase efficiency. The framework is supplemented by a local refinement of the simulation grid. To calculate the fluxes on such grids, we employ a combination of the standard two-point flux approximation and a multipoint flux approximation where the grid is refined. This is then used to simulate a large-scale example related to underground CO2 storage. 相似文献
8.
Multiscale mixed/mimetic methods on corner-point grids 总被引：1，自引：0，他引：1
Multiscale simulation is a promising approach to facilitate direct simulation of large and complex grid models for highly
heterogeneous petroleum reservoirs. Unlike traditional simulation, approaches based on upscaling/downscaling, multiscale methods
seek to solve the full flow problem by incorporating subscale heterogeneities into local discrete approximation spaces. We
consider a multiscale formulation based on a hierarchical grid approach, where basis functions with subgrid resolution are
computed numerically to correctly and accurately account for subscale variations from an underlying (fine-scale) geomodel
when solving the global flow equations on a coarse grid. By using multiscale basis functions to discretise the global flow
equations on a (moderately sized) coarse grid, one can retain the efficiency of an upscaling method and, at the same time,
produce detailed and conservative velocity fields on the underlying fine grid. For pressure equations, the multiscale mixed
finite-element method (MsMFEM) has been shown to be a particularly versatile approach. In this paper, we extend the method
to corner-point grids, which is the industry standard for modelling complex reservoir geology. To implement MsMFEM, one needs
a discretisation method for solving local flow problems on the underlying fine grids. In principle, any stable and conservative
method can be used. Here, we use a mimetic discretisation, which is a generalisation of mixed finite elements that gives a
discrete inner product, allows for polyhedral elements, and can (easily) be extended to curved grid faces. The coarse grid
can, in principle, be any partition of the subgrid, where each coarse block is a connected collection of subgrid cells. However,
we argue that, when generating coarse grids, one should follow certain simple guidelines to achieve improved accuracy. We
discuss partitioning in both index space and physical space and suggest simple processing techniques. The versatility and
accuracy of the new multiscale mixed methodology is demonstrated on two corner-point models: a small Y-shaped sector model
and a complex model of a layered sedimentary bed. A variety of coarse grids, both violating and obeying the above mentioned
guidelines, are employed. The MsMFEM solutions are compared with a reference solution obtained by direct simulation on the
subgrid. 相似文献
9.
The use of limited global information in multiscale simulations is needed when there is no scale separation. Previous approaches entail fine-scale simulations in the computation of the global information. The computation of the global information is expensive. In this paper, we propose the use of approximate global information based on partial upscaling. A requirement for partial homogenization is to capture long-range (non-local) effects present in the fine-scale solution, while homogenizing some of the smallest scales. The local information at these smallest scales is captured in the computation of basis functions. Thus, the proposed approach allows us to avoid the computations at the scales that can be homogenized. This results in coarser problems for the computation of global fields. We analyze the convergence of the proposed method. Mathematical formalism is introduced, which allows estimating the errors due to small scales that are homogenized. The proposed method is applied to simulate two-phase flows in heterogeneous porous media. Numerical results are presented for various permeability fields, including those generated using two-point correlation functions and channelized permeability fields from the SPE Comparative Project (Christie and Blunt, SPE Reserv Evalu Eng 4:308–317, 2001). We consider simple cases where one can identify the scales that can be homogenized. For more general cases, we suggest the use of upscaling on the coarse grid with the size smaller than the target coarse grid where multiscale basis functions are constructed. This intermediate coarse grid renders a partially upscaled solution that contains essential non-local information. Numerical examples demonstrate that the use of approximate global information provides better accuracy than purely local multiscale methods. 相似文献
10.
This paper extends the multipoint flux-approximation (MPFA) control-volume method to quadrilateral grids for which the adjacent cells do not necessarily share corners. Examples are grids with faults and locally refined grids. This paper gives a derivation of the method for such grids. The difference between two-point flux-approximation (TPFA) results and MPFA results for faults and local grid refinements is demonstrated for synthetic problems. Further, the results are compared with results from uniform fine-grid simulations. The effect of repeated fault patterns as well as anisotropy is investigated. Large errors may be found for the TPFA method for flow through a series of faults in an anisotropic medium. Finally, a comparison is done for a reservoir field application. 相似文献
11.
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. 相似文献
12.
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.
K.-A. Lie O. Møyner J. R. Natvig A. Kozlova K. Bratvedt S. Watanabe Z. Li 《Computational Geosciences》2017,21(5-6):981-998
For the past 10 years or so, a number of so-called multiscale methods have been developed as an alternative approach to upscaling and to accelerate reservoir simulation. The key idea of all these methods is to construct a set of prolongation operators that map between unknowns associated with cells in a fine grid holding the petrophysical properties of the geological reservoir model and unknowns on a coarser grid used for dynamic simulation. The prolongation operators are computed numerically by solving localized flow problems, much in the same way as for flow-based upscaling methods, and can be used to construct a reduced coarse-scale system of flow equations that describe the macro-scale displacement driven by global forces. Unlike effective parameters, the multiscale basis functions have subscale resolution, which ensures that fine-scale heterogeneity is correctly accounted for in a systematic manner. Among all multiscale formulations discussed in the literature, the multiscale restriction-smoothed basis (MsRSB) method has proved to be particularly promising. This method has been implemented in a commercially available simulator and has three main advantages. First, the input grid and its coarse partition can have general polyhedral geometry and unstructured topology. Secondly, MsRSB is accurate and robust when used as an approximate solver and converges relatively fast when used as an iterative fine-scale solver. Finally, the method is formulated on top of a cell-centered, conservative, finite-volume method and is applicable to any flow model for which one can isolate a pressure equation. We discuss numerical challenges posed by contemporary geomodels and report a number of validation cases showing that the MsRSB method is an efficient, robust, and versatile method for simulating complex models of real reservoirs. 相似文献
15.
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. 相似文献
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.
The paper is devoted to the upscaling method appropriate for single-phase flow in media with discontinuous permeability distribution.
The suggested algorithm is a modification of the iterative adaptive local–global upscaling developed by Chen and coauthors.
The key feature of this method is a consistency between local and coarse global calculated characteristics. In this work,
we apply a modified procedure to determine the boundary conditions used in the local fine-scale computation. To increase the
accuracy of these boundary conditions on each iteration, we involve an additional preliminary step based on the results of
coarse scale calculations from the previous iteration. Numerical tests show an essential improvement of the accuracy of upscaled
flow rates for most of the realizations of statistical permeability distribution. Although the developed method is universal,
its efficiency increases with increasing of permeability contrast. 相似文献
18.
Thomas H. Robey 《Mathematical Geology》1995,27(6):709-729
Geostatistical techniques allow simulation of properties such as porosity or conductivity on a fine scale. Typically, porous media flow modeling is performed at a coarser scale. Upscaling properties from the fine scale to the coarser scale introduces potential errors which are constrained by the degree of homogeneity of the cell or element. Adaptive grid techniques can be used to minimize the heterogeneity in the individual cells or elements, thus minimizing potential upscaling errors. A geostatistical adaptive grid (GAG) algorithm based on local minimization of heterogeneity is introduced. Local minimization allows greater control over the type of distortion permitted. Comparisons are made with a general elastic grid adjustment (GEGA) algorithm based on global minimization of heterogeneity. Several sample problems are used to test and demonstrate the two approaches. 相似文献
19.
Use of Border Regions for Improved Permeability Upscaling 总被引：1，自引：0，他引：1
A procedure for the improved calculation of upscaled grid block permeability tensors on Cartesian grids is described and applied. The method entails the use of a border region of fine-scale cells surrounding the coarse block for which the upscaled permeability is to be computed. The implementation allows for the use of full-tensor permeability fields on the fine and coarse scales. Either periodic or pressure–no flow boundary conditions are imposed over the extended local domain (target block plus border regions) though averaged quantities, used to compute the upscaled permeability tensor, are computed only over the target block region. Flow and transport results using this procedure are compared to those from standard methods for different types of geological and simulation models. Improvement using the new approach is consistently observed for the cases considered, though the degree of improvement varies for different models and flow quantities. 相似文献
20.
We present a new nonlinear monotone finite volume method for diffusion equation and its application to two-phase flow model. We consider full anisotropic discontinuous diffusion or permeability tensors on conformal polyhedral meshes. The approximation of the diffusive flux uses the nonlinear two-point stencil which provides the conventional seven-point stencil for the discrete diffusion operator on cubic meshes. We show that the quality of the discrete flux in a reservoir simulator has great effect on the front behavior and the water breakthrough time. We compare two two-point flux approximations (TPFA), the proposed nonlinear TPFA and the conventional linear TPFA, and multipoint flux approximation (MPFA). The new nonlinear scheme has a number of important advantages over the traditional linear discretizations. Compared to the linear TPFA, the nonlinear TPFA demonstrates low sensitivity to grid distortions and provides appropriate approximation in case of full anisotropic permeability tensor. For nonorthogonal grids or full anisotropic permeability tensors, the conventional linear TPFA provides no approximation, while the nonlinear flux is still first-order accurate. The computational work for the new method is higher than the one for the conventional TPFA, yet it is rather competitive. Compared to MPFA, the new scheme provides sparser algebraic systems and thus is less computational expensive. Moreover, it is monotone which means that the discrete solution preserves the nonnegativity of the differential solution. 相似文献