Generalized Multiscale Inversion for Heterogeneous Problems

In this work, we propose a generalized multiscale inversion algorithm for heterogeneous problems that aims at solving an inverse problem on a computational coarse grid. Previous inversion techniques for multiscale problems seek a coarse-grid media properties, e.g., permeability and conductivity, and by doing so, they assume that there exists a homogenized representation of the underlying fine-scale permeability field on a coarse grid. Generally such assumptions do not hold for highly heterogeneous fields, e.g., fracture media or channelized fields, where the width of channels are very small compared to the coarse-grid sizes. In these cases, grid refinement can lead to many degrees of freedom, and thus unattractive to apply. The proposed algorithm is based on the Generalized Multiscale Finite Element Method (GMsFEM), which uses local spectral problems to identify non-localized features, i.e., channels (high-conductivity inclusions that connect the boundaries of the coarse-grid block). The inclusion of these features in the coarse space enables one to achieve a good accuracy. The approach is valid under the assumption that the solution can be well represented in a reduced-dimensional space by multiscale basis functions. In practice, these basis functions are non-observable as we do not identify the fine-scale features of the permeability field. Our inversion algorithm finds the discretization parameters of the resulting system. By doing so, we identify the appropriate coarse-grid parameters representing the permeability field instead of fine-grid permeability field. We illustrate the approach by numerical results for fractured media.

[1]  Victor M. Calo,et al.  Randomized Oversampling for Generalized Multiscale Finite Element Methods , 2014, Multiscale Model. Simul..

[2]  Thomas Y. Hou,et al.  A Multiscale Finite Element Method for Elliptic Problems in Composite Materials and Porous Media , 1997 .

[3]  Mårten Gulliksson,et al.  A separating oscillation method of recovering the G-limit in standard and non-standard homogenization problems , 2016 .

[4]  H. Kazemi,et al.  NUMERICAL SIMULATION OF WATER-OIL FLOW IN NATURALLY FRACTURED RESERVOIRS , 1976 .

[5]  Yalchin Efendiev,et al.  Domain Decomposition Preconditioners for Multiscale Flows in High Contrast Media: Reduced Dimension Coarse Spaces , 2010, Multiscale Model. Simul..

[6]  Yu-Shu Wu,et al.  A multiple-porosity method for simulation of naturally fractured petroleum reservoirs , 1988 .

[7]  Bangti Jin,et al.  Multilevel Markov Chain Monte Carlo Method for High-Contrast Single-Phase Flow Problems , 2014, 1402.5068.

[8]  Yalchin Efendiev,et al.  Multiscale finite element methods for high-contrast problems using local spectral basis functions , 2011, J. Comput. Phys..

[9]  Bangti Jin,et al.  Inverse Problems , 2014, Series on Applied Mathematics.

[10]  Yalchin Efendiev,et al.  Coupling of multiscale and multi-continuum approaches , 2017, GEM - International Journal on Geomathematics.

[11]  Yalchin Efendiev,et al.  Generalized multiscale finite element methods (GMsFEM) , 2013, J. Comput. Phys..

[12]  J. E. Warren,et al.  The Behavior of Naturally Fractured Reservoirs , 1963 .

[13]  G. I. Barenblatt,et al.  Basic concepts in the theory of seepage of homogeneous liquids in fissured rocks [strata] , 1960 .

[14]  Todd Arbogast,et al.  Derivation of the double porosity model of single phase flow via homogenization theory , 1990 .

[15]  H. Engl,et al.  Regularization of Inverse Problems , 1996 .

[16]  G. A. Pavliotis,et al.  Multiscale modelling and inverse problems , 2010, 1009.2943.

[17]  Ivan P. Gavrilyuk,et al.  Lagrange multiplier approach to variational problems and applications , 2010, Math. Comput..

[18]  Christina Frederick,et al.  Numerical methods for multiscale inverse problems , 2014 .

[19]  Yalchin Efendiev,et al.  Multiscale model reduction for shale gas transport in fractured media , 2015, Computational Geosciences.

[20]  P. Maass,et al.  Sparsity regularization for parameter identification problems , 2012 .

[21]  Guanglian Li,et al.  Low-Rank Approximation to Heterogeneous Elliptic Problems , 2018, Multiscale Model. Simul..

[22]  Peng Chen,et al.  Uncertainty propagation using infinite mixture of Gaussian processes and variational Bayesian inference , 2015, J. Comput. Phys..

[23]  Yalchin Efendiev,et al.  Adaptive multiscale model reduction with Generalized Multiscale Finite Element Methods , 2016, J. Comput. Phys..

[24]  T. Narasimhan,et al.  On fluid reserves and the production of superheated steam from fractured, vapor‐dominated geothermal reservoirs , 1982 .

[25]  E. Somersalo,et al.  Statistical and computational inverse problems , 2004 .

[26]  Barbara Kaltenbacher,et al.  Regularization Methods in Banach Spaces , 2012, Radon Series on Computational and Applied Mathematics.

[27]  Y. Efendiev,et al.  A Domain Decomposition Preconditioner for Multiscale High-Contrast Problems , 2011 .

[28]  Albert Tarantola,et al.  Inverse problem theory - and methods for model parameter estimation , 2004 .