Tensor based geology preserving reservoir parameterization with Higher Order Singular Value Decomposition (HOSVD)

Parameter estimation through robust parameterization techniques has been addressed in many works associated with history matching and inverse problems. Reservoir models are in general complex, nonlinear, and large-scale with respect to the large number of states and unknown parameters. Thus, having a practical approach to replace the original set of highly correlated unknown parameters with non-correlated set of lower dimensionality, that captures the most significant features comparing to the original set, is of high importance. Furthermore, de-correlating system's parameters while keeping the geological description intact is critical to control the ill-posedness nature of such problems. We introduce the advantages of a new low dimensional parameterization approach for reservoir characterization applications utilizing multilinear algebra based techniques like higher order singular value decomposition (HOSVD). In tensor based approaches like HOSVD, 2D permeability images are treated as they are, i.e., the data structure is kept as it is, whereas in conventional dimensionality reduction algorithms like SVD data has to be vectorized. Hence, compared to classical methods, higher redundancy reduction with less information loss can be achieved through decreasing present redundancies in all dimensions. In other words, HOSVD approximation results in a better compact data representation with respect to least square sense and geological consistency in comparison with classical algorithms. We examined the performance of the proposed parameterization technique against SVD approach on the SPE10 benchmark reservoir model as well as synthetic channelized permeability maps to demonstrate the capability of the proposed method. Moreover, to acquire statistical consistency, we repeat all experiments for a set of 1000 unknown geological samples and provide comparison using RMSE analysis. Results prove that, for a fixed compression ratio, the performance of the proposed approach outperforms that of conventional methods perceptually and in terms of least square measure. HighlightsIntroduced new permeability parameterization using High Order Singular Value Decomposition (HOSVD).HOSVD yield reduced computational time as compared to classical SVD based on the same compression ratio.The new methodology improves geological description by capturing all important features (spatial) as compared to SVD methods.The HOSVD method is general for any type of reservoir parameter.HOSVD can be applied in the optimization under the uncertainty paradigm.

[1]  Michael Andrew Christie,et al.  Tenth SPE Comparative Solution Project: a comparison of upscaling techniques , 2001 .

[2]  J. Leeuw,et al.  Principal component analysis of three-mode data by means of alternating least squares algorithms , 1980 .

[3]  A. Journel,et al.  Reservoir Modeling Using Multiple-Point Statistics , 2001 .

[4]  Andrzej Cichocki,et al.  Tensor Decompositions for Signal Processing Applications: From two-way to multiway component analysis , 2014, IEEE Signal Processing Magazine.

[5]  B. Jafarpour Sparsity-Promoting Solution of Subsurface Flow Model Calibration Inverse Problems , 2013 .

[6]  Jef Caers,et al.  Geostatistical Reservoir Modelling using Statistical Pattern Recognition , 2001 .

[7]  Joos Vandewalle,et al.  On the Best Rank-1 and Rank-(R1 , R2, ... , RN) Approximation of Higher-Order Tensors , 2000, SIAM J. Matrix Anal. Appl..

[8]  W. B. Whalley,et al.  The use of fractals and pseudofractals in the analysis of two-dimensional outlines: Review and further exploration , 1989 .

[9]  Yalchin Efendiev,et al.  Local-global multiscale model reduction for flows in high-contrast heterogeneous media , 2012, J. Comput. Phys..

[10]  M. Hinze,et al.  Proper Orthogonal Decomposition Surrogate Models for Nonlinear Dynamical Systems: Error Estimates and Suboptimal Control , 2005 .

[11]  Eduardo Gildin,et al.  Efficient Inference of Reservoir Parameter Distribution Utilizing Higher Order SVD Reparameterization , 2014 .

[12]  C. W. Harper,et al.  A FORTRAN IV program for comparing ranking algorithms in quantitative biostratigraphy , 1984 .

[13]  Jan Dirk Jansen,et al.  Improved Reservoir Management Through Optimal Control and Continuous Model Updating , 2004 .

[14]  G. Mariéthoz,et al.  Multiple-point Geostatistics: Stochastic Modeling with Training Images , 2014 .

[15]  L. De Lathauwer,et al.  Exponential data fitting using multilinear algebra: the decimative case , 2009 .

[16]  Louis J. Durlofsky,et al.  Linearized reduced-order models for subsurface flow simulation , 2010, J. Comput. Phys..

[17]  K. Aziz,et al.  Petroleum Reservoir Simulation , 1979 .

[18]  Tamara G. Kolda,et al.  Tensor Decompositions and Applications , 2009, SIAM Rev..

[19]  Jef Caers,et al.  Uncertainty Quantification in Reservoir Performance Using Distances and Kernel Methods--Application to a West Africa Deepwater Turbidite Reservoir , 2009 .

[20]  Joos Vandewalle,et al.  A Multilinear Singular Value Decomposition , 2000, SIAM J. Matrix Anal. Appl..

[21]  Demetri Terzopoulos,et al.  Multilinear Analysis of Image Ensembles: TensorFaces , 2002, ECCV.

[22]  Behnam Jafarpour,et al.  Hybrid Parameterization for Robust History Matching , 2014 .

[23]  VandewalleJoos,et al.  On the Best Rank-1 and Rank-(R1,R2,. . .,RN) Approximation of Higher-Order Tensors , 2000 .

[24]  D. Oliver,et al.  Recent progress on reservoir history matching: a review , 2011 .

[25]  J. Leeuw,et al.  Some additional results on principal components analysis of three-mode data by means of alternating least squares algorithms , 1987 .

[26]  L. Y. Hu,et al.  Multiple‐point geostatistics for modeling subsurface heterogeneity: A comprehensive review , 2008 .

[27]  Ning Liu,et al.  Inverse Theory for Petroleum Reservoir Characterization and History Matching , 2008 .

[28]  Eduardo Gildin,et al.  Heterogeneous reservoir characterization using efficient parameterization through higher order SVD (HOSVD) , 2014, 2014 American Control Conference.

[29]  Benoy K. Ghose Strech: a subroutine for stretching time series and its use in stratigraphic correlation , 1984 .

[30]  L. Durlofsky,et al.  Efficient real-time reservoir management using adjoint-based optimal control and model updating , 2006 .

[31]  Anand Rangarajan,et al.  Image Denoising Using the Higher Order Singular Value Decomposition , 2013, IEEE Transactions on Pattern Analysis and Machine Intelligence.

[32]  Narendra Ahuja,et al.  Rank-R approximation of tensors using image-as-matrix representation , 2005, 2005 IEEE Computer Society Conference on Computer Vision and Pattern Recognition (CVPR'05).

[33]  Siep Weiland,et al.  Reconstruction and Approximation of Multidimensional Signals Described by Proper Orthogonal Decompositions , 2008, IEEE Transactions on Signal Processing.

[34]  J. Jansen,et al.  Reduced-order optimal control of water flooding using proper orthogonal decomposition , 2006 .

[35]  J. Wade Davis,et al.  Statistical Pattern Recognition , 2003, Technometrics.

[36]  Eduardo Gildin,et al.  Permeability Parametrization Using Higher Order Singular Value Decomposition (HOSVD) , 2013, 2013 12th International Conference on Machine Learning and Applications.