排序方式: 共有6条查询结果,搜索用时 0 毫秒
1
1.
We present a generic, semi-automated algorithm for generating non-uniform coarse grids for modeling subsurface flow. The method is applicable to arbitrary grids and does not impose smoothness constraints on the coarse grid. One therefore avoids conventional smoothing procedures that are commonly used to ensure that the grids obtained with standard coarsening procedures are not too rough. The coarsening algorithm is very simple and essentially involves only two parameters that specify the level of coarsening. Consequently the algorithm allows the user to specify the simulation grid dynamically to fit available computer resources, and, e.g., use the original geomodel as input for flow simulations. This is of great importance since coarse grid-generation is normally the most time-consuming part of an upscaling phase, and therefore the main obstacle that has prevented simulation workflows with user-defined resolution. We apply the coarsening algorithm to a series of two-phase flow problems on both structured (Cartesian) and unstructured grids. The numerical results demonstrate that one consistently obtains significantly more accurate results using the proposed non-uniform coarsening strategy than with corresponding uniform coarse grids with roughly the same number of cells. 相似文献
2.
How contact metamorphism can trigger global climate changes: Modeling gas generation around igneous sills in sedimentary basins 总被引:2,自引:0,他引:2
Ingrid Aarnes Henrik Svensen Yuri Y. Podladchikov 《Geochimica et cosmochimica acta》2010,74(24):7179-7195
Large volumes of greenhouse gases such as CH4 and CO2 form by contact metamorphism of organic-rich sediments in aureoles around sill intrusions in sedimentary basins. Thermogenic gas generation and dehydration reactions in shale are treated numerically in order to quantify basin-scale devolatilization. We show that aureole thicknesses, defined as the zone of elevated metamorphism relative to the background level, vary within 30-250% of the sill thickness, depending on the temperature of the host-rock and intrusion, besides the sill thickness. In shales with total organic carbon content of >5 wt.%, CH4 is the dominant volatile (85-135 kg/m3) generated through organic cracking, relative to H2O-generation from dehydration reactions (30-110 kg/m3). Even using conservative estimates of melt volumes, extrapolation of our results to the scale of sill complexes in a sedimentary basin indicates that devolatilization can have generated ∼2700-16200 Gt CH4 in the Karoo Basin (South Africa), and ∼600-3500 Gt CH4 in the Vøring and Møre basins (offshore Norway). The generation of volatiles is occurring on a time-scale of 10-1000 years within an aureole of a single sill, which makes the rate of sill emplacement the time-constraining factor on a basin-scale. This study demonstrates that thousands of gigatons of potent greenhouse gases like methane can be generated during emplacement of Large Igneous Provinces in sedimentary basins. 相似文献
3.
J?rg E. Aarnes Vegard Kippe Knut-Andreas Lie Alf Birger Rustad 《Computational Geosciences》2002,6(1):1-4
Introduction
Upscaling of Geological Models for Reservoir Simulation: Issues and Approaches 相似文献4.
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. 相似文献
5.
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. 相似文献
6.
Christakos Konstantinos Furevik Birgitte R. Aarnes Ole Johan Breivik Øyvind Tuomi Laura Byrkjedal Øyvind 《Ocean Dynamics》2020,70(1):57-75
Ocean Dynamics - Accurate predictions of surface ocean waves in coastal areas are important for a number of marine activities. In complex coastlines with islands and fjords, the quality of wind... 相似文献
1