共查询到20条相似文献，搜索用时 388 毫秒
1.
The present paper proposes a new family of multiscale finite volume methods. These methods usually deal with a dual mesh resolution, where the pressure field is solved on a coarse mesh, while the saturation fields, which may have discontinuities, are solved on a finer reservoir grid, on which petrophysical heterogeneities are defined. Unfortunately, the efficiency of dual mesh methods is strongly related to the definition of up-gridding and down-gridding steps, allowing defining accurately pressure and saturation fields on both fine and coarse meshes and the ability of the approach to be parallelized. In the new dual mesh formulation we developed, the pressure is solved on a coarse grid using a new hybrid formulation of the parabolic problem. This type of multiscale method for pressure equation called multiscale hybrid-mixed method (MHMM) has been recently proposed for finite elements and mixed-finite element approach (Harder et al. 2013). We extend here the MH-mixed method to a finite volume discretization, in order to deal with large multiphase reservoir models. The pressure solution is obtained by solving a hybrid form of the pressure problem on the coarse mesh, for which unknowns are fluxes defined on the coarse mesh faces. Basis flux functions are defined through the resolution of a local finite volume problem, which accounts for local heterogeneity, whereas pressure continuity between cells is weakly imposed through flux basis functions, regarded as Lagrange multipliers. Such an approach is conservative both on the coarse and local scales and can be easily parallelized, which is an advantage compared to other existing finite volume multiscale approaches. It has also a high flexibility to refine the coarse discretization just by refinement of the lagrange multiplier space defined on the coarse faces without changing nor the coarse nor the fine meshes. This refinement can also be done adaptively w.r.t. a posteriori error estimators. The method is applied to single phase (well-testing) and multiphase flow in heterogeneous porous media. 相似文献
2.
A method for multiscale parameter estimation with application to reservoir history matching is presented. Starting from a
given fine-scale model, coarser models are generated using a global upscaling technique where the coarse models are tuned
to match the solution of the fine model. Conditioning to dynamic data is done by history-matching the coarse model. Using
consistently the same resolution both for the forward and inverse problems, this model is successively refined using a combination
of downscaling and history matching until model-matching dynamic data are obtained at the finest scale. Large-scale corrections
are obtained using fast models, which, combined with a downscaling procedure, provide a better initial model for the final
adjustment on the fine scale. The result is thus a series of models with different resolution, all matching history as good
as possible with this grid. Numerical examples show that this method may significantly reduce the computational effort and/or
improve the quality of the solution when achieving a fine-scale match as compared to history-matching directly on the fine
scale. 相似文献
3.
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. 相似文献
4.
This paper presents a novel mass-conservative mixed multiscale method for solving flow equations in heterogeneous porous media. The media properties (the permeability) contain multiple scales and high contrast. The proposed method solves the flow equation in a mixed formulation on a coarse grid by constructing multiscale basis functions. The resulting velocity field is mass-conservative on the fine grid. Our main goal is to obtain first-order convergence in terms of the mesh size which is independent of local contrast. This is achieved, first, by constructing some auxiliary spaces, which contain global information that cannot be localized, in general. This is built on our previous work on the generalized multiscale finite element method (GMsFEM). In the auxiliary space, multiscale basis functions corresponding to small (contrast-dependent) eigenvalues are selected. These basis functions represent the high-conductivity channels (which connect the boundaries of a coarse block). Next, we solve local problems to construct multiscale basis functions for the velocity field. These local problems are formulated in the oversampled domain, taking into account some constraints with respect to auxiliary spaces. The latter allows fast spatial decay of local solutions and, thus, allows taking smaller oversampled regions. The number of basis functions depends on small eigenvalues of the local spectral problems. Moreover, multiscale pressure basis functions are needed in constructing the velocity space. Our multiscale spaces have a minimal dimension, which is needed to avoid contrast dependence in the convergence. The method’s convergence requires an oversampling of several layers. We present an analysis of our approach. Our numerical results confirm that the convergence rate is first order with respect to the mesh size and independent of the contrast. 相似文献
5.
Luz Angélica Caudillo-Mata Eldad Haber Christoph Schwarzbach 《Computational Geosciences》2017,21(5-6):963-980
In order to reduce the computational cost of the simulation of electromagnetic responses in geophysical settings that involve highly heterogeneous media, we develop a multiscale finite volume method with oversampling for the quasi-static Maxwell’s equations in the frequency domain. We assume a coarse mesh nested within a fine mesh that accurately discretizes the problem. For each coarse cell, we independently solve a local version of the original Maxwell’s system subject to linear boundary conditions on an extended domain, which includes the coarse cell and a neighborhood of fine cells around it. The local Maxwell’s system is solved using the fine mesh contained in the extended domain and the mimetic finite volume method. Next, these local solutions (basis functions) together with a weak-continuity condition are used to construct a coarse-mesh version of the global problem. The basis functions can be used to obtain the fine-mesh details from the solution of the coarse-mesh problem. Our approach leads to a significant reduction in the size of the final system of equations and the computational time, while accurately approximating the behavior of the fine-mesh solutions. We demonstrate the performance of our method using two 3D synthetic models: one with a mineral deposit in a geologically complex medium and one with random isotropic heterogeneous media. Both models are discretized using an adaptive mesh refinement technique. 相似文献
6.
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. 相似文献
7.
Implementation of a Locally Conservative Numerical Subgrid Upscaling Scheme for Two-Phase Darcy Flow 总被引：1，自引：0，他引：1
Todd Arbogast 《Computational Geosciences》2002,6(3-4):453-481
We present a locally mass conservative scheme for the approximation of two-phase flow in a porous medium that allows us to obtain detailed fine scale solutions on relatively coarse meshes. The permeability is assumed to be resolvable on a fine numerical grid, but limits on computational power require that computations be performed on a coarse grid. We define a two-scale mixed finite element space and resulting method, and describe in detail the solution algorithm. It involves a coarse scale operator coupled to a subgrid scale operator localized in space to each coarse grid element. An influence function (numerical Greens function) technique allows us to solve these subgrid scale problems independently of the coarse grid approximation. The coarse grid problem is modified to take into account the subgrid scale solution and solved as a large linear system of equations posed over a coarse grid. Finally, the coarse scale solution is corrected on the subgrid scale, providing a fine grid representation of the solution. Numerical examples are presented, which show that near-well behavior and even extremely heterogeneous permeability barriers and streaks are upscaled well by the technique. 相似文献
8.
We pay a revisit to some classical geomechanics problems using a novel computational multiscale modelling approach. The multiscale approach employs a hierarchical coupling of the finite element method (FEM) and the discrete element method. It solves a boundary value problem at the continuum scale by FEM and derives the material point response from the discrete element method simulation attached to each Gauss point of the FEM mesh. The multiscale modelling framework not only helps successfully bypass phenomenological constitutive assumptions as required in conventional modelling approaches but also facilitates effective cross‐scale interpretation and understanding of soil behaviour. We examine the classical retaining wall and footing problems by this method and demonstrate that the simulating results can be well validated and verified by their analytical solutions. Furthermore, the study sheds novel multiscale insights into these classical problems and offers a new tool for geotechnical engineers to design and analyse geotechnical applications based directly upon particle‐level information of soils. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献
9.
In this paper we study the problem of determining the effective permeability on a coarse scale level of problems with strongly
varying and discontinuous coefficients defined on a fine scale. The upscaled permeability is defined as the solution of an
optimization problem, where the difference between the fine scale and the coarse scale velocity field is minimized. We show
that it is not necessary to solve the fine scale pressure equation in order to minimize the associated cost‐functional. Furthermore,
we derive a simple technique for computing the derivatives of the cost‐functional needed in the fix‐point iteration used to
compute the optimal permeability on the coarse mesh. Finally, the method is illustrated by several analytical examples and
numerical experiments.
This revised version was published online in July 2006 with corrections to the Cover Date. 相似文献
10.
The prediction of fluid flows within hydrocarbon reservoirs requires the characterization of petrophysical properties. Such characterization is performed on the basis of geostatistics and history-matching; in short, a reservoir model is first randomly drawn, and then sequentially adjusted until it reproduces the available dynamic data. Two main concerns typical of the problem under consideration are the heterogeneity of rocks occurring at all scales and the use of data of distinct resolution levels. Therefore, referring to sequential Gaussian simulation, this paper proposes a new stochastic simulation method able to handle several scales for both continuous or discrete random fields. This method adds flexibility to history-matching as it boils down to the multiscale parameterization of reservoir models. In other words, reservoir models can be updated at either coarse or fine scales, or both. Parameterization adapts to the available data; the coarser the scale targeted, the smaller the number of unknown parameters, and the more efficient the history-matching process. This paper focuses on the use of variational optimization techniques driven by the gradual deformation method to vary reservoir models. Other data assimilation methods and perturbation processes could have been envisioned as well. Last, a numerical application case is presented in order to highlight the advantages of the proposed method for conditioning permeability models to dynamic data. For simplicity, we focus on two-scale processes. The coarse scale describes the variations in the trend while the fine scale characterizes local variations around the trend. The relationships between data resolution and parameterization are investigated. 相似文献
11.
Subsurface flows are affected by geological variability over a range of length scales. The modeling of well singularity in
heterogeneous formations is important for simulating flow in aquifers and petroleum reservoirs. In this paper, two approaches
in calculating the upscaled well index to capture the effects of fine scale heterogeneity in near-well regions are presented
and applied. We first develop a flow-based near-well upscaling procedure for geometrically flexible grids. This approach entails
solving local well-driven flows and requires the treatment of geometric effects due to the nonalignment between fine and coarse
scale grids. An approximate coarse scale well model based on a well singularity analysis is also proposed. This model, referred
to as near-well arithmetic averaging, uses only the fine scale permeabilities at well locations to compute the coarse scale
well index; it does not require solving any flow problems. These two methods are systematically tested on three-dimensional
models with a variety of permeability distributions. It is shown that both approaches provide considerable improvement over
a simple (arithmetic) averaging approach to compute the coarse scale well index. The flow-based approach shows close agreement
to the fine scale reference model, and the near-well arithmetic averaging also offers accuracy for an appropriate range of
parameters. The interaction between global flow and near-well upscaling is also investigated through the use of global fine
scale solutions in near-well scale-up calculations. 相似文献
12.
This paper presents an improved discrete element model, which incorporates stochastically distorted contact mechanics, for the simulation of double-twisted hexagonal wire meshes that are commonly used in rockfall protection. First, the characteristics of such meshes are investigated by conducting quasi-static and dynamic experimental tests. Second, the discrete model for the simulation of such meshes is presented. A stochastically distorted contact model is introduced to account for distortions of the wires and hexagons, allowing a more realistic representation of the mechanical response of the mesh from the deformation point of view and the force point of view. Quasi-static tensile tests of a plane net sheet, subjected to a constant strain rate, are used to study the effect of the stochastically distorted contact formulation and to calibrate the numerical model. Finally, the dynamic response of an impacting block on a horizontal mesh sheet is used to compare the numerical predictions against experimental results in order to validate the proposed approach. 相似文献
13.
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. 相似文献
14.
15.
A finite element, variable mesh analysis of unconfined steady-state seepage problems is presented based on a nonlinear programming algorithm. It is shown that the minimization of an objective function which merely represents a measure of the total flux leaving or entering the mesh at the free surface nodes (except those that belong also to pervious boundaries) does not permit a unique definition of the free surface geometry. This problem, which is apparently related to the numerical instabilities often met when using variable mesh approaches, can be eliminated by adding to the objective function a term representing a sort of overall ‘regularity’ condition for the shape of the free surface. The modified solution procedure turns out to be stable and able to provide meaningful results for practical problems even when rather coarse meshes are adopted. 相似文献
16.
Numerical simulation of fluid flow coupled with chemical reactions has been an active field in the hydrogeology community and many formulations have been programmed into different software. In recent years, this subject has attracted increasing interest in the reservoir simulation community, partly for the application of chemical methods for hydrocarbon extraction but also for research on the geological sequestration of CO2. In this paper, an extension to the concept of dual mesh for reactive transport modeling is presented. This approach involves two meshes, a low-resolution mesh to resolve the pressure equation and a high-resolution mesh to transport the species and to calculate the geochemical equilibrium. The main objective is to preserve the fine scale heterogeneities to reach a more accurate field behavior simulation than conventional approach which consist in performing simulations on a coarser mesh. The method is applied to a simulation of CO2 storage in the SPE10 model that keep a high resolution of the heterogeneities. 相似文献
17.
The problem of multiphase phase flow in heterogeneous subsurface porous media is one involving many uncertainties. In particular, the permeability of the medium is an important aspect of the model that is inherently uncertain. Properly quantifying these uncertainties is essential in order to make reliable probabilistic-based predictions and future decisions. In this work, a measure-theoretic framework is employed to quantify uncertainties in a two-phase subsurface flow model in high-contrast media. Given uncertain saturation data from observation wells, the stochastic inverse problem is solved numerically in order to obtain a probability measure on the space of unknown permeability parameters characterizing the two-phase flow. As solving the stochastic inverse problem requires a number of forward model solves, we also incorporate the use of a conservative version of the generalized multiscale finite element method for added efficiency. The parameter-space probability measure is used in order to make predictions of saturation values where measurements are not available, and to validate the effectiveness of the proposed approach in the context of fine and coarse model solves. A number of numerical examples are offered to illustrate the measure-theoretic methodology for solving the stochastic inverse problem using both fine and coarse solution schemes. 相似文献
18.
基于云南省S214思茅-江城二级公路某一边坡,利用有限差分软件FLAC3D中动力分析模块,研究了地震作用下锚杆支护的上覆红黏土岩体边坡的动力响应。基于土体与支护结构相互作用及其协同工作,建立了三维有限差分模型,给出了阻尼和动力作用下边界的选取方法,分析了地震作用下锚杆支护上覆红黏土岩体边坡的动力响应规律,研究了在地震作用下锚杆支护边坡的抗震效果。结果表明：地震作用后各层锚杆轴力和砂浆的剪应力都有所增大,但每层锚杆轴力的增幅都各不相同,锚杆轴力沿全长分布不均匀,且各层锚杆轴力均在红黏土与基岩的交界处最大,剪应力则表现为整体增大的趋势而且最大值向坡面靠近,地震作用前、后剪应力的最小值都是在红黏土与基岩的交界处;地震作用下锚杆的支护很好地限制了边坡的变形,加大了边坡的延性,具有很好的抗震性能;边坡在地震作用下产生了永久位移;坡体内加速度在竖向随高程增加而增大;平台的设置削弱了坡面加速度的增大趋势,起到了一定减震作用。研究结论对锚杆支护边坡的抗震设计与动力分析有一定参考价值。 相似文献
19.
It is becoming easier to combine geographical data and dynamic models to provide information for problem solving and geographical cognition. However, the scale dependencies of the data, model, and process can confuse the results. This study extends traditional scale research in static geographical patterns to dynamic processes and focuses on the combined scale effect of multiscale geographical data and dynamic models. The capacity for topographical expression under the combined scale effect was investigated by taking multiscale topographical data and meteorological processes in Hong Kong as a case study. A meteorological simulation of the combined scale effect was evaluated against data from Hong Kong Observatory stations. The experiments showed that (1) a digital elevation model (DEM) using 3 arc sec data with a 1 km resolution Weather Research and Forecasting (WRF) model gives better topographical expression and meteorological reproduction in Hong Kong; (2) a fine-scale model is sensitive to the resolution of the DEM data, whereas a coarse-scale model is less sensitive to it; (3) better topographical expression alone does not improve weather process simulation; and (4) uncertainty arising from a scale mismatch between the DEM data and the dynamic model may account for 38 % of the variance in certain meteorological variables (e.g., temperature). This case study gives a clear explanation of the significance and implementation of scale matching for multiscale geographical data and dynamic models. 相似文献
20.
脆性颗粒材料的动态多尺度模型研究 总被引：1，自引：0，他引：1
脆性颗粒材料的多尺度模型一般包含微观尺度的基本粒子、细观尺度的颗粒和宏观尺度的颗粒堆积体3个尺度。基于离散元方法（DEM）构建多尺度模型,并将该模型应用于动态加载。首先,对多尺度模型所涉及的两种接触模型和两种黏结模型的参数进行分析,详细讨论微细观模型参数与宏观材料常数之间的联系。然后,选用Hertz-Mindlin接触模型[1]和平行键黏结模型,建造石英砂的动态多尺度模型。通过选择合适的强度和局部阻尼参数发现,模型宏细观尺度上的动态压缩响应与对石英砂的相关试验结果吻合很好。利用多尺度模型和选定的参数,探讨与动态加载密切相关的局部阻尼机制对多尺度模型各个尺度上力学响应的影响。结果表明,阻尼越大则颗粒材料对波的衰减能力越强,但过高的阻尼会使团簇强度和模型的宏观压缩曲线都表现出异常的加载速度效应（后者实际是阻尼引起的微惯性效应）。另外,高阻尼会过度衰减颗粒破碎过程产生的应力波,从而阻碍颗粒破碎。最后,应用改进的动态多尺度模型,对脆性颗粒材料的动态破碎特性进行研究,发现该模型不但能给出与试验相吻合的颗粒级配曲线,还能揭示出颗粒破碎过程中微裂纹分布的空间不均匀性,即颗粒破碎过程中波的产生机制和衰减机制相互作用导致的微裂纹聚团分布的现象。 相似文献