Low-rank Kalman filtering for efficient state estimation of subsurface advective contaminant transport models
KAUST DepartmentApplied Mathematics and Computational Science Program
Computational Transport Phenomena Lab
Earth Fluid Modeling and Prediction Group
Earth Science and Engineering Program
Environmental Science and Engineering Program
Physical Science and Engineering (PSE) Division
Permanent link to this recordhttp://hdl.handle.net/10754/562146
MetadataShow full item record
AbstractAccurate knowledge of the movement of contaminants in porous media is essential to track their trajectory and later extract them from the aquifer. A two-dimensional flow model is implemented and then applied on a linear contaminant transport model in the same porous medium. Because of different sources of uncertainties, this coupled model might not be able to accurately track the contaminant state. Incorporating observations through the process of data assimilation can guide the model toward the true trajectory of the system. The Kalman filter (KF), or its nonlinear invariants, can be used to tackle this problem. To overcome the prohibitive computational cost of the KF, the singular evolutive Kalman filter (SEKF) and the singular fixed Kalman filter (SFKF) are used, which are variants of the KF operating with low-rank covariance matrices. Experimental results suggest that under perfect and imperfect model setups, the low-rank filters can provide estimates as accurate as the full KF but at much lower computational effort. Low-rank filters are demonstrated to significantly reduce the computational effort of the KF to almost 3%. © 2012 American Society of Civil Engineers.
CitationEl Gharamti, M., Hoteit, I., & Sun, S. (2012). Low-Rank Kalman Filtering for Efficient State Estimation of Subsurface Advective Contaminant Transport Models. Journal of Environmental Engineering, 138(4), 446–457. doi:10.1061/(asce)ee.1943-7870.0000484
SponsorsThis publication utilized work supported in part by funds from the KAUST GCR Collaborative Fellow program.
Showing items related by title, author, creator and subject.
Particle Kalman Filtering: A Nonlinear Framework for Ensemble Kalman FiltersHoteit, Ibrahim; Luo, Xiaodong; Pham, Dinh-Tuan; Moroz, Irene M. (AIP Publishing, 2010-09-20) [Conference Paper]Optimal nonlinear filtering consists of sequentially determining the conditional probability distribution functions (pdf) of the system state, given the information of the dynamical and measurement processes and the previous measurements. Once the pdfs are obtained, one can determine different estimates, for instance, the minimum variance estimate, or the maximum a posteriori estimate, of the system state. It can be shown that, many filters, including the Kalman filter (KF) and the particle filter (PF), can be derived based on this sequential Bayesian estimation framework. In this contribution, we present a Gaussian mixture-based framework, called the particle Kalman filter (PKF), and discuss how the different EnKF methods can be derived as simplified variants of the PKF. We also discuss approaches to reducing the computational burden of the PKF in order to make it suitable for complex geosciences applications. We use the strongly nonlinear Lorenz-96 model to illustrate the performance of the PKF.
Particle Kalman Filtering: A Nonlinear Bayesian Framework for Ensemble Kalman Filters*Hoteit, Ibrahim; Luo, Xiaodong; Pham, Dinh-Tuan (Monthly Weather Review, American Meteorological Society, 2012-02) [Article]This paper investigates an approximation scheme of the optimal nonlinear Bayesian filter based on the Gaussian mixture representation of the state probability distribution function. The resulting filter is similar to the particle filter, but is different from it in that the standard weight-type correction in the particle filter is complemented by the Kalman-type correction with the associated covariance matrices in the Gaussian mixture. The authors show that this filter is an algorithm in between the Kalman filter and the particle filter, and therefore is referred to as the particle Kalman filter (PKF). In the PKF, the solution of a nonlinear filtering problem is expressed as the weighted average of an “ensemble of Kalman filters” operating in parallel. Running an ensemble of Kalman filters is, however, computationally prohibitive for realistic atmospheric and oceanic data assimilation problems. For this reason, the authors consider the construction of the PKF through an “ensemble” of ensemble Kalman filters (EnKFs) instead, and call the implementation the particle EnKF (PEnKF). It is shown that different types of the EnKFs can be considered as special cases of the PEnKF. Similar to the situation in the particle filter, the authors also introduce a resampling step to the PEnKF in order to reduce the risk of weights collapse and improve the performance of the filter. Numerical experiments with the strongly nonlinear Lorenz-96 model are presented and discussed.
GaukharSH/mlenkf: Julia Implementation of Multilevel Ensemble Kalman Filtering (MLEnKF) with Local Kalman GainsHoel, Håkon; Shaimerdenova, Gaukhar; Tempone, Raul (Github, 2019-12-22) [Software]Julia Implementation of Multilevel Ensemble Kalman Filtering (MLEnKF) with Local Kalman Gains