Ensemble-Based Assimilation of Nonlinearly Related Dynamic Data in Reservoir Models Exhibiting Non-Gaussian Characteristics

Devesh Kumar, Sanjay Srinivasan

Research output: Contribution to journalArticle

Abstract

Inverse modeling techniques for estimating reservoir parameters (e.g., transmissivity, permeability) utilize some dynamic (secondary) information (e.g., hydraulic head or production data at well locations). Ensemble based data assimilation methods are one such class of inverse modeling techniques. Ensemble Kalman filter (EnKF) in particular is built around the basic framework of linearly updating a state vector ψ f to ψ a , based on observed dynamic data. Components of the state vector are reservoir model parameters such as transmissivity, permeability, storativity, porosity, hydraulic head, and phase-saturation. Although EnKF is able to update a large number of parameters successively as data become available, it has some shortcomings. It is optimal only in the cases where multivariate joint distribution describing the relationship between components in the state vector is Gaussian. The relationship between the model parameters and the observed data space is quantified in terms of a covariance function, which is adequate to describe a multi-Gaussian distribution. These assumptions and simplifications may result in updated models that are inconsistent with the prior geology (exhibiting non-Gaussian characteristics) and that, in turn, may yield inaccurate predictions of reservoir performance. The aim of this research work is to propose a novel method for data assimilation which is free from the Gaussian assumptions. This new method can be used to assimilate dynamic data into reservoir models using an ensemble-based approach. Model updates are performed after transforming the model parameters into indicator variables. The indicator transform is insensitive to non-linear operations, and it thus allows us to overcome some of the limitations of the traditional EnKF.

Original languageEnglish (US)
Pages (from-to)75-107
Number of pages33
JournalMathematical Geosciences
Volume51
Issue number1
DOIs
StatePublished - Jan 10 2019

Fingerprint

Ensemble
Ensemble Kalman Filter
Kalman filter
Inverse Modeling
Data Assimilation
hydraulic head
transmissivity
data assimilation
Hydraulics
Permeability
Update
Model
permeability
F-vector
Geology
Covariance Function
Multivariate Distribution
research work
Inaccurate
Porosity

All Science Journal Classification (ASJC) codes

  • Mathematics (miscellaneous)
  • Earth and Planetary Sciences(all)

Cite this

@article{1d2647b09485477d8a4f5622e6f1d006,
title = "Ensemble-Based Assimilation of Nonlinearly Related Dynamic Data in Reservoir Models Exhibiting Non-Gaussian Characteristics",
abstract = "Inverse modeling techniques for estimating reservoir parameters (e.g., transmissivity, permeability) utilize some dynamic (secondary) information (e.g., hydraulic head or production data at well locations). Ensemble based data assimilation methods are one such class of inverse modeling techniques. Ensemble Kalman filter (EnKF) in particular is built around the basic framework of linearly updating a state vector ψ f to ψ a , based on observed dynamic data. Components of the state vector are reservoir model parameters such as transmissivity, permeability, storativity, porosity, hydraulic head, and phase-saturation. Although EnKF is able to update a large number of parameters successively as data become available, it has some shortcomings. It is optimal only in the cases where multivariate joint distribution describing the relationship between components in the state vector is Gaussian. The relationship between the model parameters and the observed data space is quantified in terms of a covariance function, which is adequate to describe a multi-Gaussian distribution. These assumptions and simplifications may result in updated models that are inconsistent with the prior geology (exhibiting non-Gaussian characteristics) and that, in turn, may yield inaccurate predictions of reservoir performance. The aim of this research work is to propose a novel method for data assimilation which is free from the Gaussian assumptions. This new method can be used to assimilate dynamic data into reservoir models using an ensemble-based approach. Model updates are performed after transforming the model parameters into indicator variables. The indicator transform is insensitive to non-linear operations, and it thus allows us to overcome some of the limitations of the traditional EnKF.",
author = "Devesh Kumar and Sanjay Srinivasan",
year = "2019",
month = "1",
day = "10",
doi = "10.1007/s11004-018-9762-x",
language = "English (US)",
volume = "51",
pages = "75--107",
journal = "Mathematical Geosciences",
issn = "1874-8961",
publisher = "Springer Netherlands",
number = "1",

}

Ensemble-Based Assimilation of Nonlinearly Related Dynamic Data in Reservoir Models Exhibiting Non-Gaussian Characteristics. / Kumar, Devesh; Srinivasan, Sanjay.

In: Mathematical Geosciences, Vol. 51, No. 1, 10.01.2019, p. 75-107.

Research output: Contribution to journalArticle

TY - JOUR

T1 - Ensemble-Based Assimilation of Nonlinearly Related Dynamic Data in Reservoir Models Exhibiting Non-Gaussian Characteristics

AU - Kumar, Devesh

AU - Srinivasan, Sanjay

PY - 2019/1/10

Y1 - 2019/1/10

N2 - Inverse modeling techniques for estimating reservoir parameters (e.g., transmissivity, permeability) utilize some dynamic (secondary) information (e.g., hydraulic head or production data at well locations). Ensemble based data assimilation methods are one such class of inverse modeling techniques. Ensemble Kalman filter (EnKF) in particular is built around the basic framework of linearly updating a state vector ψ f to ψ a , based on observed dynamic data. Components of the state vector are reservoir model parameters such as transmissivity, permeability, storativity, porosity, hydraulic head, and phase-saturation. Although EnKF is able to update a large number of parameters successively as data become available, it has some shortcomings. It is optimal only in the cases where multivariate joint distribution describing the relationship between components in the state vector is Gaussian. The relationship between the model parameters and the observed data space is quantified in terms of a covariance function, which is adequate to describe a multi-Gaussian distribution. These assumptions and simplifications may result in updated models that are inconsistent with the prior geology (exhibiting non-Gaussian characteristics) and that, in turn, may yield inaccurate predictions of reservoir performance. The aim of this research work is to propose a novel method for data assimilation which is free from the Gaussian assumptions. This new method can be used to assimilate dynamic data into reservoir models using an ensemble-based approach. Model updates are performed after transforming the model parameters into indicator variables. The indicator transform is insensitive to non-linear operations, and it thus allows us to overcome some of the limitations of the traditional EnKF.

AB - Inverse modeling techniques for estimating reservoir parameters (e.g., transmissivity, permeability) utilize some dynamic (secondary) information (e.g., hydraulic head or production data at well locations). Ensemble based data assimilation methods are one such class of inverse modeling techniques. Ensemble Kalman filter (EnKF) in particular is built around the basic framework of linearly updating a state vector ψ f to ψ a , based on observed dynamic data. Components of the state vector are reservoir model parameters such as transmissivity, permeability, storativity, porosity, hydraulic head, and phase-saturation. Although EnKF is able to update a large number of parameters successively as data become available, it has some shortcomings. It is optimal only in the cases where multivariate joint distribution describing the relationship between components in the state vector is Gaussian. The relationship between the model parameters and the observed data space is quantified in terms of a covariance function, which is adequate to describe a multi-Gaussian distribution. These assumptions and simplifications may result in updated models that are inconsistent with the prior geology (exhibiting non-Gaussian characteristics) and that, in turn, may yield inaccurate predictions of reservoir performance. The aim of this research work is to propose a novel method for data assimilation which is free from the Gaussian assumptions. This new method can be used to assimilate dynamic data into reservoir models using an ensemble-based approach. Model updates are performed after transforming the model parameters into indicator variables. The indicator transform is insensitive to non-linear operations, and it thus allows us to overcome some of the limitations of the traditional EnKF.

UR - http://www.scopus.com/inward/record.url?scp=85053469549&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85053469549&partnerID=8YFLogxK

U2 - 10.1007/s11004-018-9762-x

DO - 10.1007/s11004-018-9762-x

M3 - Article

AN - SCOPUS:85053469549

VL - 51

SP - 75

EP - 107

JO - Mathematical Geosciences

JF - Mathematical Geosciences

SN - 1874-8961

IS - 1

ER -