Reservoir description by using a Piecewise Constant Level Set Method

We consider the permeability estimation problem in two-phase porous media flow. We try to identify the permeability field by utilizing both the production data from wells as well as inverted seismic data. The permeability field is assumed to be piecewise constant, or can be approximated well by a piecewise constant function. A variant of the level set method, called Piecewise Constant Level Set Method is used to represent the interfaces between the regions with different permeability levels. The inverse problem is solved by minimizing a functional, and TV norm regularization is used to deal with the ill-posedness. We also use the operator-splitting technique to decompose the constraint term from the fidelity term. This gives us more flexibility to deal with the constraint and helps to stabilize the algorithm.