14-17 May 2018
New Orleans
US/Central timezone

Multiscale Data Assimilation of Spatially Distributed Information

16 May 2018, 15:58
New Orleans

New Orleans

Oral 20 Minutes MS 2.08: Recent Advances in Multiscale Methods and Uncertainty Quantification Parallel 8-H


Rafael Moraes (TU Delft)


In data assimilation problems, various types of data are naturally linked to different spatial scales (e.g. seismic and electromagnetic data), and these scales are usually not coincident to the subsurface simulation model scale. Alternatives like down/upscaling [1] 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 gradient via a MS strategy [2, 3]. The fine-scale uncertain parameters are directly updated and the MS grid [4, 5] is constructed in a resolution that meets the observed data resolution. The advantages of the technique are demonstrated in the assimilation of data represented in 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 numerical experiments show that the MS strategy provides advantages on the computational side – expensive operations are only performed at the coarse scale; on the accuracy – the matched uncertain parameter distribution is closer to the “truth” model; and in the optimization performance – faster convergence behavior and smaller OF value are obtained. As conclusion, the newly developed method is capable of providing superior results when compared to strategies that rely on the downscaling of the observed data, addressing the scale dissimilarity via a robust, consistent MS strategy.


[1] Aanonsen, S. I. (2005, January). Efficient history matching using a multiscale technique. In SPE Reservoir Simulation Symposium. Society of Petroleum Engineers.
[2] Moraes, R. J., Rodrigues, J. R., Hajibeygi, H., & Jansen, J. D. (2017). Multiscale gradient computation for flow in heterogeneous porous media. Journal of Computational Physics, 336, 644-663.
[3] Moraes, R., Rodrigues, J. R. P., Hajibeygi, H., & Jansen, J. D. (2017, February). Multiscale Gradient Computation for Multiphase Flow in Porous Media. In SPE Reservoir Simulation Conference. Society of Petroleum Engineers.
[4] Hou, T. Y., & Wu, X. H. (1997). A multiscale finite element method for elliptic problems in composite materials and porous media. Journal of computational physics, 134(1), 169-189.
[5] Jenny, P., Lee, S. H., & Tchelepi, H. A. (2003). Multi-scale finite-volume method for elliptic problems in subsurface flow simulation. Journal of Computational Physics, 187(1), 47-67.

Acceptance of Terms and Conditions Click here to agree

Primary author

Rafael Moraes (TU Delft)


Hadi Hajibeygi (TU Delft) Prof. Jan Dirk Jansen (TU Delft)

Presentation Materials

There are no materials yet.