Hostname: page-component-7bb8b95d7b-dvmhs Total loading time: 0 Render date: 2024-10-03T16:39:38.914Z Has data issue: false hasContentIssue false

Binary Level Set Methods for Dynamic Reservoir Characterization by Operator Splitting Scheme

Published online by Cambridge University Press:  03 June 2015

Changhui Yao*
Affiliation:
Department of Mathematics, Zhengzhou University, No. 100 of Science Road, Zhengzhou 450001, Henan, China
*
*Corresponding author. Email: chyao@lsec.cc.ac.cn (C. H. Yao)
Get access

Abstract

In this paper, operator splitting scheme for dynamic reservoir characterization by binary level set method is employed. For this problem, the absolute permeability of the two-phase porous medium flow can be simulated by the constrained augmented Lagrangian optimization method with well data and seismic time-lapse data. By transforming the constrained optimization problem in an unconstrained one, the saddle point problem can be solved by Uzawas algorithms with operator splitting scheme, which is based on the essence of binary level set method. Both the simple and complicated numerical examples demonstrate that the given algorithms are stable and efficient and the absolute permeability can be satisfactorily recovered.

Type
Research Article
Copyright
Copyright © Global-Science Press 2012

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)

References

[1]Aanonsen, S. I. and Eydinov, D., A multiscale method for distributed parameter estimation with application to reservoir history matching, Comput. Geosci., 10(1) (2006), pp. 97117.Google Scholar
[2]Bai, Z.-Z., Huang, Y.-M., Ng, M. K. and Yang, X., Preconditioned iterative methods for algebraic systems from multiplicative half-quadratic regularization image restorations, Numer. Math. Theor. Meth. Appl., 3 (2010), pp. 461474.Google Scholar
[3]Barkve, T.Cominelli, A.Gosselin, O., Kolasinski, M., Reme, H., Aanonsen, S. I., AND Aavatsmark, I., Effect of scale dependent data correlation in an integrated history matching loop combining production data and 4d seismeic data, Technical report, presented at the SPE Reservoir Simulation Symposium held in Houston, Texas, Feb.2003.Google Scholar
[4]Chan, T. and Vese, L. A., Active contours without edges, IEEE Tran. image Proc., 10 (2001), pp. 266277.CrossRefGoogle ScholarPubMed
[5]Chan, T. F. and Tai, X.-C., Augmented lagrangian and total variation methods for recovering discontinuous coefficients from elliptic equations, in Bristeau, M., Etgen, G., Fitzgibbon, W., Lions, J. L., Periaux, J. and Wheeler, M. F., editors, Computational Science for the 21st Century, pp. 597607, John Wiley & Sons, 1997.Google Scholar
[6]Gosselin, O., Aavatsmark, I., Iaanonsen, S., Cominelli, A. and Barkve, T., Inte-grationn of 4d data in the history match loop by investigating scale dependent correlations in the acustic impedence cube, Technical report, Proc. 8th European Conference on Mathematics of Oil Recovery, Freiberg, German, 2002.Google Scholar
[7]Guan, W. and Hu, H. S., The parameter averaging technique in finite-difference modeling of elastic waves in combined structures with solid, fluid and porous subregions, Commun. Com-put. Phys., 10 (2011), pp. 695715.Google Scholar
[8]Ito, K. and Kunisch, K., The augmented Lagrangian method for parameter estimation in elliptic systems, SIAM J. Control Optim., 28 (1990), pp. 113136.Google Scholar
[9]Ito, K., Li, Z.-L. and Qiao, Z., The sensitivity analysis for the flow past obstacles problem with respect to the Reynolds number, Adv. Appl. Math. Mech., 4 (2012), pp. 2135.Google Scholar
[10]Ito, K. and Qiao, Z., A high order compact MAC finite difference scheme for the stokes equations: augmented variable approach, J. Comput. Phys., 227(17) (2008), pp. 81778190.Google Scholar
[11]Li, F. and Ng, M. K., Kernel density estimation based multiphase fuzzy region competition method for texture image segmentation, Commun. Comput. Phys., 8 (2010), pp. 623641.Google Scholar
[12]Li, H.-W. and Tai, X.-C., Piecewise constant level set method for multiphase motion, Int. J. Numer. Anal. Model., 4(2) (2007), pp. 291305.Google Scholar
[13]Lie, J., Lysaker, M. and Tai, X.-C., A binary level set model and some applications for Mumford-Shah image segmentation, IEEE Trans. Image Proc., 15(4) (2006), pp. 11711181.CrossRefGoogle ScholarPubMed
[14]Lie, J., Lysaker, M. and Tai, X.-C., A variant of the level set method and applications to image segmentation, Math. Comput., 75(255) (2006), pp. 11551174.Google Scholar
[15]Lu, T., Neittaanmaäki, P. and Tai, X.-C., A parallel splitting up method and its application to Navier-Stokes equations, Appl. Math. Lett., 4(2) (1991), pp. 2529.Google Scholar
[16]Lu, T., Neittaanmaäki, P. and Tai, X.-C., A parallel splitting-up method for partial differential equations and its applications to Navier-Stokes equations, RAIRO Modél. Math. Anal. Numér., 26(6) (1992), pp. 673708.Google Scholar
[17]Nielsen, L. K., Li, H.-W., Tai, X.-C., Aanonsen, S. I. and Espedal, M., Reservoir description using a binary level set model, Comput. Vis. Sci., 13(1) (2009), pp. 4158.Google Scholar
[18]Nielsen, L. K., Tai, X.-C., Aanonsen, S. I. and Espedal, M., Identifying reservoir facies using a binary level method, SPE, 2007, accepted.Google Scholar
[19]Nielsen, L. K., Tai, X.-C., Aanonsen, S. I. and Espedal, M., Reservoir decription using a binary level set model, CAM-report-05-50, UCLA, Applied Mathematics, 2005.Google Scholar
[20]Nielsen, L. K., Tai, X.-C., Aanonsen, Sigurd Ivar and Espedal, Magne, A binary level set model for elliptic inverse problems with discontinuous coefficients, Int. J. Numer. Anal. Model., 4(1) (2007), pp. 7499.Google Scholar
[21]Rada, L. and Chen, K., A new variational model with dual level set functions for selective segmentation, Commun. Comput. Phys., 12 (2012), pp. 261283.Google Scholar
[22]Reynolds, A. C., Li, R. and Oliver, D. S., History matching of three-phase flow production data, SPE Journal, 8(4), December 2003.Google Scholar
[23]Ren, K., Recent developments in numerical techniques for transport-based medical imaging methods, Commun. Comput. Phys., 8 (2010), pp. 150.Google Scholar
[24]Song, P. and Sun, S. Y., Contaminant flow and transport simulation in cracked porous media using locally conservative schemes, Adv. Appl. Math. Mech., 4 (2012), pp. 389421.Google Scholar
[25]Tai, X.-C., Christiansen, O., Lin, P. and Skjaelaaen, I., A remark on the MBO scheme and some piecewise constant level set methods, CAM-report-05-24, UCLA, Applied Mathematics, 2005.Google Scholar
[26]Wang, L.-L. and Gu, Y., Efficient dual algorithms for image segmentation using TV-Allen-Cahn type models, Commun. Comput. Phys., 9 (2011), pp. 859877.CrossRefGoogle Scholar
[27]Wang, S.-G., Zhao, G.-Z., Xu, L.-B., Guo, D. and Sun, S.-Y., Optimization for automatic history matching, Int. J. Numer. Anal. Model., 2 (2005), pp. 131137.Google Scholar
[28]Zhu, J.-P. and Chen, Y. M., Multilevel grid method for history matching multidimensional multiphase reservoir models, Appl. Numer. Math., 10(2) (1992), pp. 159174.Google Scholar