1887

Abstract

Summary

We introduce a novel approach to compute the permeability on very large CT images of core samples. The LIR-tree is used for spatial partitioning of the geometry. The images are coarsened in areas where the velocity and pressure do not vary much while keeping the original resolution near the solid surfaces. In addition, solid regions do not occupy computational memory. Variables are arranged in a way such that each cell is able to satisfy the Stokes-equations independently from its neighbor cells. Pressure and velocity are discretized on staggered grids but instead of using one velocity variable on the cell faces we introduce two velocity variables. They discretize the two one-sided limits at the center of the cell surface. The discretization of the momentum and mass conservation equations yields one small linear system per cell. This structure allows to use the block Gauß-Seidel-algorithm as iterative solver.

We compare our method to three other solvers on a large complex Berea sandstone dataset. The method is 3–6 times faster, scales well for up to 32 processors and has very low memory requirements. An additional benefit of this work is that the range of permeabilities from the benchmark can be narrowed down significantly.

Loading

Article metrics loading...

/content/papers/10.3997/2214-4609.20140834
2014-06-16
2024-04-25
Loading full text...

Full text loading...

References

  1. Andrä, H. et al.
    [2013a] Digital rock physics benchmarks—Part I: Imaging and segmentation. Computers & Geosciences, 50, 25–32, ISSN 0098-3004, benchmark problems, datasets and methodologies for the computational geosciences.
    [Google Scholar]
  2. [2013b] Digital rock physics benchmarks—part II: Computing effective properties. Computers & Geosciences, 50, 33–43, ISSN 0098-3004, benchmark problems, datasets and methodologies for the computational geosciences.
    [Google Scholar]
  3. Cheng, L., Wiegmann, A. and Rief, S.
    [2013] SIMPLE-FFT for flow computations in low porosity μCT images. 5th International Conference on Porous Media & Annual Meeting, Czech Republic, 21–24 May, 2013.
    [Google Scholar]
  4. Kanit, T., Forest, S., Galliet, I., Mounoury, V. and Jeulin, D.
    [2003] Determination of the size of the representative volume element for random composites: statistical and numerical approach. International Journal of Solids and Structures, 40(13–14), 3647–3679, ISSN 0020-7683.
    [Google Scholar]
  5. Keehm, Y.
    [2003] Computational rock physics: Transport properties in porous media and applications. Ph.D. thesis, Stanford University.
    [Google Scholar]
  6. Linden, S.
    [2013] The LIR Space Partitioning System applied to the Stokes Equations. Ph.D. thesis, University of Kaiserslautern.
    [Google Scholar]
  7. Linden, S., Wiegmann, A. and Hagen, H.
    [2013] The LIR Space Partitioning System applied to Cartesian Grids, to appear at Eighth International Conference on Mathematical Methods for Curves and Surfaces, Oslo, June 28–July 3, 2012.
    [Google Scholar]
  8. Patankar, S.V.
    [1980] Numerical heat transfer and fluid flow. Taylor & Francis.
    [Google Scholar]
  9. Wiegmann, A. and Kenneth, P.B.
    [1997] The Explicit-Jump Immersed Interface Method: Finite Difference Methods For PDE With Piecewise Smooth Solutions. SIAM J. Numer. Anal, 37, 827–862.
    [Google Scholar]
http://instance.metastore.ingenta.com/content/papers/10.3997/2214-4609.20140834
Loading
/content/papers/10.3997/2214-4609.20140834
Loading

Data & Media loading...

This is a required field
Please enter a valid email address
Approval was a Success
Invalid data
An Error Occurred
Approval was partially successful, following selected items could not be processed due to error