A multiscale method for subsurface inverse modeling: Single-phase transient flow

2019

Comput Geosci

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.

Section: Introductionmentioning

confidence: 99%

A multiscale method for data assimilation

2019

Comput Geosci

“…These are efficient simulation strategies capable of solving the flow problem in a coarser grid, while being capable of accurately representing fine scale heterogeneities. Fu et al (2010Fu et al ( , 2011) developed a MS adjoint method for single-phase flow in porous media. The most efficient data assimilation methods are gradient-based, with gradients obtained with the adjoint method Oliver et al (2008).…”

Section: Introductionmentioning

confidence: 99%

A Multiscale Method For Data Assimilation

2018

Proceedings

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 down/upscaling of the data and/or the simulation model can be used, but with potential loss of important information. 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 down/upscaling of any quantity. The method relies on a recently developed mathematical framework to compute adjoint gradients via a MS strategy. The fine-scale uncertain parameters are directly updated and the MS grid is constructed in a resolution that meets the observed data resolution. The advantages of the technique are demonstrated in the 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 performance of the method is demonstrated in synthetic models and compared to down/upscaling strategies. The experiments show that the MS strategy provides advantages 1) on the computational side -expensive operations are only performed at the coarse scale; 2) with respect to accuracy -the matched uncertain parameter distribution is closer to the "truth"; and 3) in the optimization performance -faster convergence behaviour due to faster gradient computation. In conclusion, the newly developed method is capable of providing superior results when compared to strategies that rely on the up/downscaling of the response/observed data, addressing the scale dissimilarity via a robust, consistent MS strategy.

“…In order to address the computational efficiency and dissimilarity of scales in data assimilation, the present work focuses on the computation of MS adjoint gradients applied to CAHM studies. In this context, Fu et al (2011) show that a MSFV adjoint formulation can be efficiently applied to inverse problems of single-phase flow in heterogeneous porous media. Furthermore, Moraes et al (2016) proposes a generic framework for the computation of MS gradient information for single-phase flow models.…”

Section: Introductionmentioning

confidence: 99%

Multiscale Gradient Computation for Multiphase Flow in Porous Media

Rodrigues

Day 3 Wed, February 22, 2017

2017

A multiscale gradient computation method for multiphase flow in heterogeneous porous media is developed. The method constructs multiscale primal and dual coarse grids, imposed on the given fine-scale computational grid. Local multiscale basis functions are computed on (dual-) coarse blocks, constructing an accurate map (prolongation operator) between coarse-and fine-scale systems. While the expensive operations involved in computing the gradients are performed at the coarse scale, sensitivities with respect to uncertain parameters (e.g., grid block permeabilities) are expressed in the fine scale via the partial derivatives of the prolongation operator. Hence, the method allows for updating of the geological model, rather than the dynamic model only, avoiding upscaling and the inevitable loss of information. The formulation and implementation are based on automatic differentiation (AD), allowing for convenient extensions to complex physics. An IMPES coupling strategy for flow and transport is followed, in the forward simulation. The flow equation is computed using a multiscale finite volume (MSFV) formulation and the transport equation is computed at the fine scale, after reconstruction of mass conservative velocity field. To assess the performance of the method, a synthetic multiphase flow test case is considered. The multiscale gradients are compared against those obtained from a fine-scale reference strategy. Apart from its computational efficiency, the benefits of the method include flexibility to accommodate variables expressed at different scales, specially in multiscale data assimilation and reservoir management studies.