Adjoint Numerical Method for a Multiphysical Inverse Problem of Two-Phase Well Testing in Petroleum Reservoirs

Abstract

Inverse problem solution is an integral part of data interpretation for well testing in petroleum reservoirs. In case of two-phase well tests with water injection, forward problem is based on the multiphase flow model in porous media and solved numerically. The inverse problem is based on a misfit or likelihood objective function. Adjoint methods have proved robust and efficient for gradient calculation of the objective function in this type of problems. However, if time-lapse electrical resistivity measurements during the well test are included in the objective function, both the forward and inverse problems become multiphysical, and straightforward application of the adjoint method is problematic. In this paper we present a novel adjoint algorithm for the inverse problems considered. It takes into account the structure of cross dependencies between flow and electrical equations and variables, as well as specifics of the equations (mixed parabolic-hyperbolic for flow and elliptic for electricity), numerical discretizations and grids, and measurements in the inverse problem. Decomposition is proposed for the adjoint problem which makes possible step-wise solution of the electric adjoint equations, like in the forward problem, after which a cross-term is computed and added to the right-hand side of the flow adjoint equations at this timestep. The overall procedure provides accurate gradient calculation for the multiphysical objective function while preserving robustness and efficiency of the adjoint methods. Example cases of the adjoint gradient calculation are presented and compared to straightforward difference-based gradient calculation in terms of accuracy and efficiency.