Reservoir description using a binary level set model

  • Authors:

  • Lars Kristian Nielsen;Hongwei Li;Xue-Cheng Tai;Sigurd Ivar Aanonsen;Magne Espedal

  • Affiliations:

  • Univ. of Bergen, Dept. of Math., 5020, Bergen and Univ. of Bergen, CIPR-Centre for Integrated Petroleum Res., 5020, Bergen and Laksevåg Upper Secondary/High Sch., Bergen, Norway;University of Bergen, CIPR-Centre for Integrated Petroleum Research, 5020, Bergen, Norway and Capital Normal University, Department of Mathematics, 10037, Beijing, People’s Republic of Chin ...;University of Bergen, Department of Mathematics, 5020, Bergen, Norway and Nanyang Technological University, Division of Mathematical Sciences, School of Physical and Mathematical Sciences, 637616, ...;University of Bergen, Department of Mathematics, 5020, Bergen, Norway and University of Bergen, CIPR-Centre for Integrated Petroleum Research, 5020, Bergen, Norway;University of Bergen, Department of Mathematics, 5020, Bergen, Norway and University of Bergen, CIPR-Centre for Integrated Petroleum Research, 5020, Bergen, Norway

  • Venue:
  • Computing and Visualization in Science
  • Year:
  • 2009

Quantified Score

Hi-index 0.00

Visualization

Abstract

We consider the inverse problem of permeability estimation for two-phase flow in porous media. In the parameter estimation process we utilize both data from the wells (production data) and spatially distributed data (from time-lapse seismic data). The problem is solved by approximating the permeability field by a piecewise constant function, where we allow the discontinuity curves to have arbitrary shape with some forced regularity. To achieve this, we have utilized level set functions to represent the permeability field and applied an additional total variation regularization. The optimization problem is solved by a variational augmented Lagrangian approach. A binary level set formulation is used to determine both the curves of discontinuities and the constant values for each region. We do not need any initial guess for the geometries of the discontinuities, only a reasonable guess of the constant levels is required.