Inversion of Geothermal Reservoir Models Using the Adjoint Method and Randomized Low-Rank Matrix Approximation Algorithms

ResearchSpace Repository

Show simple item record

dc.contributor.advisor O'Sullivan, MJ en
dc.contributor.advisor O'Sullivan, JP en Bjarkason, Elvar en 2019-04-04T01:29:12Z en 2018 en
dc.identifier.uri en
dc.description.abstract Solving the inverse problem of finding suitable parameters for a numerical model describing a high-enthalpy geothermal system is time-consuming, as it involves running numerous nonlinear simulations until model outputs match observations. In geothermal modelling the most popular inversion algorithms are basic implementations of the Levenberg-Marquardt method, which use derivatives of model outputs to improve model parameters iteratively. However, their computational cost scales poorly with the number of parameters Nm, since they use at least Nm Å1 nonlinear simulations and finite differencing to estimate model derivatives. Choosing Nm large is, therefore, impractical with these methods, which limits their flexibility. For reducing the computational cost of inverting highly-parameterized geothermal models, this study considered adjoint and direct methods, as well as randomized matrix algorithms. As demonstrated, the adjoint and direct methods can, by themselves, reduce the cost of standard derivative-based methods for inverting geothermal natural-state and production models. This is especially the case for natural-state models since they require solving expensive transient to steady-state simulations but evaluating model derivatives with the adjoint or direct methods only requires solving linear problems for the final natural-state time-step. For improved efficiency of highly-parameterized inversions, we proposed combining randomized low-rank matrix algorithms with adjoint and direct methods to generate approximate Levenberg-Marquardt updates. As discussed, the proposed randomized Levenberg-Marquardt variants are suited to high performance parallel computing as they do not require serial iterations to generate model updates and can, therefore, outperform previously proposed approximate Levenberg-Marquardt methods, which use iterative methods. To improve the accuracy and performance of the proposed inversion methods we developed new randomized low-rank matrix algorithms which improve on previous algorithms. These improved algorithms include generalized subspace iteration and block Krylov algorithms, and single-pass streaming algorithms. These algorithms can be applied to a wide range of problems which involve forming low-rank matrices. Finally, we advised and demonstrated how these new randomized low-rank algorithms can be used to accelerate inversion and uncertainty quantification of highly-parameterized geothermal models. The results suggest that the proposed randomized methods are useful for inversion of geothermal models, and this may also be the case for other inverse problems. en
dc.publisher ResearchSpace@Auckland en
dc.relation.ispartof PhD Thesis - University of Auckland en
dc.relation.isreferencedby UoA99265148713802091 en
dc.rights Items in ResearchSpace are protected by copyright, with all rights reserved, unless otherwise indicated. Previously published items are made available in accordance with the copyright policy of the publisher. en
dc.rights.uri en
dc.rights.uri en
dc.title Inversion of Geothermal Reservoir Models Using the Adjoint Method and Randomized Low-Rank Matrix Approximation Algorithms en
dc.type Thesis en Engineering Science en The University of Auckland en Doctoral en PhD en
dc.rights.holder Copyright: The author en
dc.rights.accessrights en
pubs.elements-id 767812 en
pubs.record-created-at-source-date 2019-04-04 en

Files in this item

Find Full text

This item appears in the following Collection(s)

Show simple item record Except where otherwise noted, this item's license is described as


Search ResearchSpace

Advanced Search