We simulate advection-diffusion flow of tracers in fractured rock geometries. Voronoi tessellation generates polygonal patterns, which we then use to introduce fractures and obtain fractured geometry. The rock geometry is discretized using a scalable, in-house developed discretization software. Lattice-Boltzmann and random-walk particle-tracking methods are employed to obtain flow field and recover tracer behavior, respectively. Tracers are allowed to cross semi-permeable interface between fractures and matrix. In addition, tracers can have variable partitioning coefficients. The implemented numerical framework allows simulating field-scale tracer experiments designed to estimate residual oil saturation. Use of HPC platform is necessary to perform such simulations in three dimensions.


Article metrics loading...

Loading full text...

Full text loading...


  1. An, C., Yan, B., Alfi, M., Mi, L., Killough, J. and Heidari, Z.
    [2017] Estimating spatial distribution of natural fractures by changing NMR relaxation with magnetic nanoparticles. J. Pet. Sci. Eng.157, 273–287.
    [Google Scholar]
  2. Bechtold, M., Vanderborght, J., Ippisch, O. and Vereecken, H.
    [2011] Efficient random walk particle tracking algorithm for advective-dispersive transport in media with discontinuous dispersion coeffcients and water contents. Water Resour. Res., 47, W10526.
    [Google Scholar]
  3. Cunningham, K., Sukop, M., Huang, H., Alvarez, P., Curran, A., Renken, R. and Joann, D.
    [2009] Prominence of ichnologically-influenced macroporosity in the karst Biscayne aquifer: Stratiform “super-K” zones. Geol. Soc. Am. Bull.121, 164–180.
    [Google Scholar]
  4. Daneyko, A., Hlushkou, D., Baranau, V., Khirevich, S., Seidel-Morgenstern, A. and Tallarek, U.
    [2015] Computational investigation of longitudinal diffusion, eddy dispersion, and trans-particle mass transfer in bulk, random packings of coreshell particles with varied shell thickness and shell diffusion coefficient. J. Chrom. A, 1407, 139–156.
    [Google Scholar]
  5. Ginzburg, I., Verhaeghe, F. and d’Humières, D.
    [2008] Two-Relaxation-Time Lattice Boltzmann Scheme: About Parametrization, Velocity, Pressure and Mixed Boundary Conditions. Commun. Comput. Phys., 3, 427–478.
    [Google Scholar]
  6. Khirevich, S., Ginzburg, I. and Tallarek, U.
    [2015] Coarse- and fine-grid numerical behavior of MRT/TRT lattice-Boltzmann schemes in regular and random sphere packings. J. Comput. Phys., 281, 708–742.
    [Google Scholar]
  7. Raghavachary, S.
    [2002] Fracture generation on polygonal meshes using Voronoi polygons. In ACM SIGGRAPH, 187–187.
    [Google Scholar]
  8. Sanni, M., Al-Abbad, M., Kokal, S., Dugstad, Ø., Hartvig, S. and Huseby, O.
    [2016] Pushing the Envelope of Residual Oil Measurement: A Field Case Study of a New Class of Inter-Well Chemical Tracers. SPE Annual Technical Conference and Exhibition, 26–28 September, Dubai, UAE.
    [Google Scholar]
  9. Szymczak, P. and Ladd, A.
    [2003] Boundary conditions for stochastic solutions of the convection diffusion equation. Phys. Rev. E68, 036704.
    [Google Scholar]

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