RESERVOIR DESCRIPTION BY USING A PIECEWISE CONSTANT LEVEL SET METHOD  

作　　者：Hongwei Li Center for Integrated Petroleum Research,University of Bergen,Norway Department of Mathematics,Capital Normal University,Beijing 100037,China Xuecheng Tai Department of Mathematics,University of Bergen,Norway Division of Mathematical Sciences,School of Physical and Mathematical Sciences,Nanyang Technological University,Singapore Sigurd Ivar Aanonsen Center for Integrated Petroleum Research,University of Bergen,Norway Department of Mathematics,University of Bergen,Norway 

出　　处：《Journal of Computational Mathematics》2008年第3期365-377,共13页计算数学（英文）

基　　金：the Norwegian Research Council,Petromaks Programme

摘　　要：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.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.

关 键 词：Inverse problem Level set method Piecewise constant Operator splitting Reservoir description 

分 类 号：P618.13[天文地球—矿床学]