4.4 Article

Randomized Tensor Decomposition for Large-Scale Data Assimilation Problems for Carbon Dioxide Sequestration

Journal

MATHEMATICAL GEOSCIENCES
Volume 54, Issue 7, Pages 1139-1163

Publisher

SPRINGER HEIDELBERG
DOI: 10.1007/s11004-022-10005-1

Keywords

Inverse problem; CO2 storage; Randomized linear algebra; High-order singular value decomposition

Funding

  1. School of Energy Resources
  2. Department of Geology and Geophysics of University of Wyoming

Ask authors/readers for more resources

A computationally efficient method using randomized tensor decomposition is developed in this study to reduce model parameters and observations for efficient data assimilation in low-dimensional spaces.
Data assimilation methods are commonly used to predict petrophysical properties of deep saline aquifers for carbon dioxide sequestration studies. However, data assimilation is usually computationally challenging for large-scale geological models with a large number of geophysical observations. A computationally efficient method is developed for large-scale data assimilation problems. In the proposed method, model parameters and observations are reduced by randomized tensor decomposition, and the inversion is performed in lower-dimensional spaces. Tensors are multiway arrays that generalize matrices to multiple dimensions and therefore provide a natural way to represent three-dimensional geological models and geophysical measurements. To exploit the inherent multidimensional structure of tensors, high-order singular value decomposition performs orthogonal decomposition of the tensor in the high-order space instead of using a flattened matrix representation. However, for complex geological models, the dimension of tensors becomes very large, making the decomposition computationally unfeasible. To overcome this limitation, the dimension of the tensors is first reduced by a randomized linear algebra algorithm which guarantees that most of the information of the original tensors is preserved with high probability in the reduced tensor, and then high-order singular value decomposition is applied on the smaller tensors with relatively small computational cost. By using randomized tensor decomposition for model and data re-parameterization, data assimilation is efficiently performed in the low-dimensional model and data space. The proposed method is demonstrated and validated on the three-dimensional data set synthesized based on the Johansen formation, a potential carbon dioxide storage site offshore in Norway.

Authors

I am an author on this paper
Click your name to claim this paper and add it to your profile.

Reviews

Primary Rating

4.4
Not enough ratings

Secondary Ratings

Novelty
-
Significance
-
Scientific rigor
-
Rate this paper

Recommended

No Data Available
No Data Available