Development of optimal models of porous media by combining static and dynamic data: the porosity distribution.

This paper is part of a project, the goal of which is the development of the optimal spatial distributions of the porosity and permeability of a large-scale porous medium by using complementary static and dynamic data for the medium. The data include limited measurements of the porosity, which the method honors (preserves) in the optimal model and utilizes its correlation function, together with the first-arrival (FA) times, at a certain number of receivers, of seismic waves that have propagated in the medium and the time dependence of the pressure of a fluid flowing in the medium. The method uses the simulated-annealing (SA) technique in order to develop the optimal model. In the present paper we utilize the porosity and FA times data in order to develop the optimal spatial distribution of the porosity. This is accomplished by combining the SA method with a simulator that solves for the numerical solution of the acoustic-wave equation from which the FA times are estimated, limited porosity, and FA times data. We show that the optimal model not only honors the data, but also provides accurate estimates of the porosities in the rest of the porous medium. The efficiency of the computations is discussed in detail.

[1]  S. Torquato,et al.  Reconstructing random media , 1998 .

[2]  Berk Nf,et al.  Scattering properties of the leveled-wave model of random morphologies. , 1991 .

[4]  D. Janaki Ram,et al.  Parallel Simulated Annealing Algorithms , 1996, J. Parallel Distributed Comput..

[5]  Hilfer,et al.  Stochastic reconstruction of sandstones , 2000, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[6]  M. A. Dablain,et al.  The application of high-order differencing to the scalar wave equation , 1986 .

[7]  S. Patz,et al.  Probing porous media with gas diffusion NMR. , 1999, Physical review letters.

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

[9]  Jack K. Cohen,et al.  Mathematics of Multidimensional Seismic Imaging, Migration, and Inversion , 2001 .

[10]  Rintoul,et al.  Reconstruction of the Structure of Dispersions , 1997, Journal of colloid and interface science.

[11]  S. Torquato,et al.  Reconstructing random media. II. Three-dimensional media from two-dimensional cuts , 1998 .

[12]  J. Thovert,et al.  Grain reconstruction of porous media: application to a low-porosity Fontainebleau sandstone. , 2001, Physical review. E, Statistical, nonlinear, and soft matter physics.

[13]  Miguel Bosch,et al.  The optimization approach to lithological tomography: Combining seismic data and petrophysics for porosity prediction , 2004 .

[14]  Thomas Bohlen,et al.  Paralel 3-D viscoelastic finite difference seismic modelling , 2002 .

[15]  Muhammad Sahimi,et al.  Computer simulation of gas generation and transport in landfills II: Dynamic conditions , 2006 .

[16]  Christian Meise On the convergence of parallel simulated annealing , 1998 .

[17]  Pierre M. Adler,et al.  The formation factor of reconstructed porous media , 1992 .

[18]  Pierre M. Adler,et al.  Flow in simulated porous media , 1990 .

[19]  M. Sahimi,et al.  Generation of long-range correlations in large systems as an optimization problem. , 2006, Physical review. E, Statistical, nonlinear, and soft matter physics.

[20]  Berk Nf,et al.  Scattering properties of a model bicontinuous structure with a well defined length scale. , 1987 .

[21]  Daniel R. Greening,et al.  Parallel simulated annealing techniques , 1990 .

[22]  J. Quiblier A new three-dimensional modeling technique for studying porous media , 1984 .

[23]  Muhammad Sahimi,et al.  Multiresolution Wavelet Scale Up of Unstable Miscible Displacements in Flow Through Heterogeneous Porous Media , 2004 .

[24]  S. Torquato Random Heterogeneous Materials , 2002 .

[25]  M. Sahimi Flow phenomena in rocks : from continuum models to fractals, percolation, cellular automata, and simulated annealing , 1993 .

[26]  Alessandro Bevilacqua,et al.  A Methodological Approach to Parallel Simulated Annealing on an SMP System , 2002, J. Parallel Distributed Comput..

[27]  John W. Cahn,et al.  Phase Separation by Spinodal Decomposition in Isotropic Systems , 1965 .

[28]  Pierre M. Adler,et al.  The effective mechanical properties of reconstructed porous media , 1996 .

[29]  Barry Hilary Valentine Topping,et al.  Parallel simulated annealing for structural optimization , 1999 .

[30]  Teubner,et al.  Transport properties of heterogeneous materials derived from Gaussian random fields: Bounds and simulation. , 1995, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[31]  Xin-Quan Ma,et al.  Simultaneous inversion of prestack seismic data for rock properties using simulated annealing , 2002 .

[32]  W. V. Pinczewski,et al.  Velocity‐porosity relationships, 1: Accurate velocity model for clean consolidated sandstones , 2003 .

[33]  Mrinal K. Sen,et al.  Nonlinear one-dimensional seismic waveform inversion using simulated annealing , 1991 .

[34]  Guido Kneib,et al.  Accurate and efficient seismic modeling in random media , 1993 .

[35]  M. Sahimi,et al.  Shape of a wave front in a heterogenous medium. , 2006, Physical review letters.

[36]  Kevin P. Dorrington,et al.  Genetic‐algorithm/neural‐network approach to seismic attribute selection for well‐log prediction , 2004 .

[37]  Roberts,et al.  Structure-property correlations in model composite materials. , 1996, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[38]  M. Teubner LEVEL SURFACES OF GAUSSIAN RANDOM-FIELDS AND MICROEMULSIONS , 1991 .

[39]  Hong Ye,et al.  Speed-up simulated annealing by parallel coordinates , 2006, Eur. J. Oper. Res..

[40]  M. Sahimi,et al.  Self-affine fractal distributions of the bulk density, elastic moduli, and seismic wave velocities of rock. , 2005, Physical review. E, Statistical, nonlinear, and soft matter physics.

[41]  Gerard T. Schuster,et al.  Wave-equation traveltime inversion , 1991 .

[42]  Jonathan Carter,et al.  Using genetic algorithms for reservoir characterisation , 2001 .