Search

In data assimilation problems, various types of data are naturally linked to different spatial resolutions (e.g., seismic and electromagnetic data), and these scales are usually not coincident to the subsurface simulation model scale. Alternatives like upscaling/downscaling of the data and/or the simulation model can be used, but with potential loss of important information. Such alternatives introduce additional uncertainties which are not in the nature of the problem description, but the result of the post processing of the data or the geo-model. To address this issue, a novel multiscale (MS) data assimilation method is introduced. The overall idea of the method is to keep uncertain parameters and observed data at their original representation scale, avoiding upscaling/downscaling of any quantity. The method relies on a recently developed mathematical framework to compute adjoint gradients via a MS strategy in an algebraic framework. The fine-scale uncertain parameters are directly updated and the MS grid is constructed in a resolution that meets the observed data resolution. This formulation therefore enables a consistent assimilation of data represented at a coarser scale than the simulation model. The misfit objective function is constructed to keep the MS nature of the problem. The regularization term is represented at the simulation model (fine) scale, whereas the data misfit term is represented at the observed data (coarse) scale. The computational aspects of the method are investigated in a simple synthetic model, including an elaborate uncertainty quantification step, and compared to upscaling/downscaling strategies. The experiment shows that the MS strategy provides several potential advantages compared to more traditional scale conciliation strategies: (1) expensive operations are only performed at the coarse scale; (2) the matched uncertain parameter distribution is closer to the "truth"; (3) faster convergence behavior occurs due to faster gradient computation; and (4) better uncertainty quantification results are obtained. The proof-of-concept example considered in this paper sheds new lights on how one can reduce uncertainty within fine-scale geo-model parameters with coarse-scale data, without the necessity of upscaling/downscaling the data nor the geo-model. The developments demonstrate how to consistently formulate such a gradient-based MS data assimilation strategy in an algebraic framework which allows for implementation in available computational platforms. [ABSTRACT FROM AUTHOR]

Field-scale simulation of flow in porous media in presence of incomplete mixing demands for high-resolution computational grids, much beyond the scope of state-of-the-art simulators. Hence, the upscaling-based Todd and Longstaff (TL) approach is typically used, where coarse grid cells are employed with effective mixing fluid properties and parameters found by matching results obtained with fully resolved reference simulations. Dynamic local grid refinement (DLGR) techniques, on the other hand, only employ fine-scale grid resolution where the fully mixed assumption is not valid. The rest of the domain is then solved at coarser resolutions, where the fully mixed assumption is valid. Here, we assess the accuracy and the robustness of DLGR- and TL-based simulations of miscible displacements in homogeneous and heterogeneous porous media. Due to the intrinsic uncertainty within the unstable displacement nature of the studied incomplete mixing processes, the performance of the methods is also investigated based on a range of acceptable solutions rather than relying only on a single reference one. Systematic numerical results illustrate that the DLGR method is much more robust and accurate than the upscaling-based TL approach, and employs only a small fraction of fine-scale reference grids. Especially, the TL upscaling results (though history matched with computationally expensive fine-scale results) are very sensitive to the change of the simulation parameters. Based on this study, we propose a dynamic multilevel simulation strategy for efficient and reliable large-scale simulation of the complex incomplete mixing processes. [ABSTRACT FROM AUTHOR]

We propose a two-stage preconditioner for accelerating the iterative solution by a Krylov subspace method of Biot's poroelasticity equations based on a displacement-pressure formulation. The spatial discretization combines a finite element method for mechanics and a finite volume approach for flow. The fully implicit backward Euler scheme is used for time integration. The result is a 2 × 2 block linear system for each timestep. The preconditioning operator is obtained by applying a two-stage scheme. The first stage is a global preconditioner that employs multiscale basis functions to construct coarse-scale coupled systems using a Galerkin projection. This global stage is effective at damping low-frequency error modes associated with long-range coupling of the unknowns. The second stage is a local block-triangular smoothing preconditioner, which is aimed at high-frequency error modes associated with short-range coupling of the variables. Various numerical experiments are used to demonstrate the robustness of the proposed solver. [ABSTRACT FROM AUTHOR]

The MultiScale Finite Volume (MSFV) method is known to produce non-monotone solutions. The causes of the non-monotone solutions are identified and connected to the local flux across the boundaries of primal coarse cells induced by the basis functions. We propose a monotone MSFV (m-MSFV) method based on a local stencil-fix that guarantees monotonicity of the coarse-scale operator, and thus, the resulting approximate fine-scale solution. Detection of non-physical transmissibility coefficients that lead to non-monotone solutions is achieved using local information only and is performed algebraically. For these 'critical' primal coarse-grid interfaces, a monotone local flux approximation, specifically, a Two-Point Flux Approximation (TPFA), is employed. Alternatively, a local linear boundary condition can be used for the dual basis functions to reduce the degree of non-monotonicity. The local nature of the two strategies allows for ensuring monotonicity in local sub-regions, where the non-physical transmissibility occurs. For practical applications, an adaptive approach based on normalized positive off-diagonal coarse-scale transmissibility coefficients is developed. Based on the histogram of these normalized coefficients, one can remove the large peaks by applying the proposed modifications only for a small fraction of the primal coarse grids. Though the m-MSFV approach can guarantee monotonicity of the solutions to any desired level, numerical results illustrate that employing the m-MSFV modifications only for a small fraction of the domain can significantly reduce the non-monotonicity of the conservative MSFV solutions. [ABSTRACT FROM AUTHOR]