Multiphase flow reconstruction in oil pipelines by capacitance tomography using simulated annealing

A highly optimized simulated annealing (SA) algorithm is applied to reconstruct permittivity images of real two-phase gas-oil flows through a cylindrical vessel using electrical capacitance tomography (ECT). ECT yields low-accuracy images but is robust, inexpensive and much faster than many other tomography processes. This non-intrusive method essentially measures non-conductive system distributions and is applied in oil industry processes such as mixing or stirring vessels, fluidized bed reactors, separator tanks and pipelines carrying multiphase flows. A forward problem is solved at each step of an iterative algorithm to solve the inverse problem using simulated annealing (SA). Comparisons with linear methods like The Projected Landweber technique are discussed. In this paper we introduce a finite volume discretization with local mesh refinements in a cylindrical configuration close to the electrodes in order to improve resolution in the calculation of capacitances, and to avoid problems with resolution at the centre of cylindrical container when finite differences are used. This discretization has the advantage of a conservative formulation as used in finite element methods and features the flexibility of mesh refinement close to the electrodes. Thus, improvement of local accuracy is achieved without increasing prohibitively the number of mesh points. Performance of the forward problem resolution is compared with finite element based methods and experimental data. We show that the non linear version of SA provides better reconstructions of threephase flows than the Landweber method.

[1]  Kris Vasudevan,et al.  Simulated annealing statics computation using an order-based energy function , 1991 .

[2]  Maurice Beck,et al.  Experimental evaluation of capacitance tomographic flow imaging systems using physical models , 1994 .

[3]  F. J. Dickin,et al.  Tomographic imaging of industrial process equipment : techniques and applications , 1992 .

[4]  Wuqiang Yang,et al.  Development of capacitance tomographic imaging systems for oil pipeline measurements , 1995 .

[5]  Wuqiang Yang,et al.  Electrical Capacitance Tomography — from Design to Applications , 1995 .

[7]  C. D. Gelatt,et al.  Optimization by Simulated Annealing , 1983, Science.

[8]  Lihui Peng,et al.  Image reconstruction algorithms for electrical capacitance tomography , 2003 .

[9]  Mrinal K. Sen,et al.  Global Optimization Methods in Geophysical Inversion , 1995 .

[10]  Ø. Isaksen,et al.  A review of reconstruction techniques for capacitance tomography , 1996 .

[11]  Dieter Mewes,et al.  Recent developments and industrial/research applications of capacitance tomography , 1996 .

[12]  Jan-Erik Nordtvedt,et al.  A new reconstruction algorithm for process tomography , 1993 .

[13]  M. Sambridge,et al.  Genetic algorithms in seismic waveform inversion , 1992 .

[14]  Brian S. Hoyle,et al.  Electrical capacitance tomography for flow imaging: system model for development of image reconstruction algorithms and design of primary sensors , 1992 .

[15]  N. Metropolis,et al.  Equation of State Calculations by Fast Computing Machines , 1953, Resonance.

[16]  Wuqiang Yang,et al.  An image-reconstruction algorithm based on Landweber's iteration method for electrical-capacitance tomography , 1999 .

[17]  Wuqiang Yang,et al.  Optimization of an iterative image reconstruction algorithm for electrical capacitance tomography , 1999 .

[18]  Supachai Bangliang,et al.  The use of simultaneous iterative reconstruction technique for electrical capacitance tomography , 2000 .

[19]  M. S. Beck,et al.  8-electrode capacitance system for two-component flow identification. I. Tomographic flow imaging , 1989 .

[20]  Wuqiang Yang,et al.  Principles and Industrial Applications of Electrical Capacitance Tomography , 1997 .

[21]  Per Christian Hansen,et al.  Rank-Deficient and Discrete Ill-Posed Problems , 1996 .

[22]  Xiaoping Wu.  A 3-D finite-element algorithm for DC resistivity modelling using the shifted incomplete Cholesky conjugate gradient method , 2003 .