3D reconstruction of porous media from a 2D section and comparisons of transport and elastic properties

Morteza Elahi Naraghi, Kyle Spikes, Sanjay Srinivasan

Research output: Contribution to journalArticle

3 Citations (Scopus)

Abstract

High-resolution, 3D microstructure images of rocks can be used to compute the transport and elastic properties of those samples by use of a digital rock-physics approach. Those properties are complex functions of the pore-size distribution, geometry, and morphology, and necessitate the use of accurate 3D volumes. Because of the limited availability of 3D microstructure images of rocks, several attempts to construct 3D images from 2D images have been made. In this study, we propose a new stochastic method to reconstruct a 3D image of the rock by use of only a 2D section of the imagedrock sample. Our method is derived from a simple observation that the pores are gradually deformed from one section to the next. Therefore, the first step is to generate multiple independent realizations by performing multiple-point-statistics (MPS)-based stochastic simulations. These simulations represent independent 2D scans through the rock volume. Next, a succession of images is generated spanning two adjacent independent sections. These images consist of gradually morphed features from one section to those in the next independent section. Juxtaposing these 2D images results in the reconstructed 3D image. We implement the algorithm on a Berea sandstone rock for which the 3D high image was available for comparison. We calculate the spatial connectivity in the third direction, and confirmed that the proposed method can retrieve the connectivity in the third direction accurately. We also compute transport and elastic properties from the reconstructed image and from the original image to verify that this method reproduces the appropriate spatial statistics and pore-size distribution, geometry, and morphology. We also compare the numerical results with laboratory measurements performed on the sample. The results obtained by use of the reconstructed image reveal that the numerically calculated properties are similar to the measured values. We compare the mismatch of transport and elastic properties with the original measurements with that of the previous reconstruction algorithm. These comparisons show that the proposed simulation method has the same accuracy as previous ones. However, the proposed method is much more computationally efficient than the other algorithms, which are dependent on simulation of all 2D layers, mainly because of the faster MPS algorithm and the fact that the simulation is being performed only for the independent layers instead of all the layers. The proposed methodology is an accurate method to reconstruct a 3D representative sample of a rock given only one 2D thin section. The algorithm is computationally efficient and faster than the previously introduced algorithm, and can easily be used to characterize samples for which 3D images are difficult to obtain in terms of both time and expense.

Original languageEnglish (US)
Pages (from-to)342-352
Number of pages11
JournalSPE Reservoir Evaluation and Engineering
Volume20
Issue number2
DOIs
StatePublished - May 1 2017

Fingerprint

elastic property
Porous materials
porous medium
Rocks
Statistics
Pore size
rock
Microstructure
Geometry
simulation
comparison
Sandstone
connectivity
microstructure
Physics
Availability
geometry
method
thin section
physics

All Science Journal Classification (ASJC) codes

  • Fuel Technology
  • Energy Engineering and Power Technology
  • Geology

Cite this

@article{b2800ff377b64881bc88633e5c26ced7,
title = "3D reconstruction of porous media from a 2D section and comparisons of transport and elastic properties",
abstract = "High-resolution, 3D microstructure images of rocks can be used to compute the transport and elastic properties of those samples by use of a digital rock-physics approach. Those properties are complex functions of the pore-size distribution, geometry, and morphology, and necessitate the use of accurate 3D volumes. Because of the limited availability of 3D microstructure images of rocks, several attempts to construct 3D images from 2D images have been made. In this study, we propose a new stochastic method to reconstruct a 3D image of the rock by use of only a 2D section of the imagedrock sample. Our method is derived from a simple observation that the pores are gradually deformed from one section to the next. Therefore, the first step is to generate multiple independent realizations by performing multiple-point-statistics (MPS)-based stochastic simulations. These simulations represent independent 2D scans through the rock volume. Next, a succession of images is generated spanning two adjacent independent sections. These images consist of gradually morphed features from one section to those in the next independent section. Juxtaposing these 2D images results in the reconstructed 3D image. We implement the algorithm on a Berea sandstone rock for which the 3D high image was available for comparison. We calculate the spatial connectivity in the third direction, and confirmed that the proposed method can retrieve the connectivity in the third direction accurately. We also compute transport and elastic properties from the reconstructed image and from the original image to verify that this method reproduces the appropriate spatial statistics and pore-size distribution, geometry, and morphology. We also compare the numerical results with laboratory measurements performed on the sample. The results obtained by use of the reconstructed image reveal that the numerically calculated properties are similar to the measured values. We compare the mismatch of transport and elastic properties with the original measurements with that of the previous reconstruction algorithm. These comparisons show that the proposed simulation method has the same accuracy as previous ones. However, the proposed method is much more computationally efficient than the other algorithms, which are dependent on simulation of all 2D layers, mainly because of the faster MPS algorithm and the fact that the simulation is being performed only for the independent layers instead of all the layers. The proposed methodology is an accurate method to reconstruct a 3D representative sample of a rock given only one 2D thin section. The algorithm is computationally efficient and faster than the previously introduced algorithm, and can easily be used to characterize samples for which 3D images are difficult to obtain in terms of both time and expense.",
author = "Naraghi, {Morteza Elahi} and Kyle Spikes and Sanjay Srinivasan",
year = "2017",
month = "5",
day = "1",
doi = "10.2118/180489-PA",
language = "English (US)",
volume = "20",
pages = "342--352",
journal = "SPE Reservoir Evaluation and Engineering",
issn = "1094-6470",
publisher = "Society of Petroleum Engineers (SPE)",
number = "2",

}

3D reconstruction of porous media from a 2D section and comparisons of transport and elastic properties. / Naraghi, Morteza Elahi; Spikes, Kyle; Srinivasan, Sanjay.

In: SPE Reservoir Evaluation and Engineering, Vol. 20, No. 2, 01.05.2017, p. 342-352.

Research output: Contribution to journalArticle

TY - JOUR

T1 - 3D reconstruction of porous media from a 2D section and comparisons of transport and elastic properties

AU - Naraghi, Morteza Elahi

AU - Spikes, Kyle

AU - Srinivasan, Sanjay

PY - 2017/5/1

Y1 - 2017/5/1

N2 - High-resolution, 3D microstructure images of rocks can be used to compute the transport and elastic properties of those samples by use of a digital rock-physics approach. Those properties are complex functions of the pore-size distribution, geometry, and morphology, and necessitate the use of accurate 3D volumes. Because of the limited availability of 3D microstructure images of rocks, several attempts to construct 3D images from 2D images have been made. In this study, we propose a new stochastic method to reconstruct a 3D image of the rock by use of only a 2D section of the imagedrock sample. Our method is derived from a simple observation that the pores are gradually deformed from one section to the next. Therefore, the first step is to generate multiple independent realizations by performing multiple-point-statistics (MPS)-based stochastic simulations. These simulations represent independent 2D scans through the rock volume. Next, a succession of images is generated spanning two adjacent independent sections. These images consist of gradually morphed features from one section to those in the next independent section. Juxtaposing these 2D images results in the reconstructed 3D image. We implement the algorithm on a Berea sandstone rock for which the 3D high image was available for comparison. We calculate the spatial connectivity in the third direction, and confirmed that the proposed method can retrieve the connectivity in the third direction accurately. We also compute transport and elastic properties from the reconstructed image and from the original image to verify that this method reproduces the appropriate spatial statistics and pore-size distribution, geometry, and morphology. We also compare the numerical results with laboratory measurements performed on the sample. The results obtained by use of the reconstructed image reveal that the numerically calculated properties are similar to the measured values. We compare the mismatch of transport and elastic properties with the original measurements with that of the previous reconstruction algorithm. These comparisons show that the proposed simulation method has the same accuracy as previous ones. However, the proposed method is much more computationally efficient than the other algorithms, which are dependent on simulation of all 2D layers, mainly because of the faster MPS algorithm and the fact that the simulation is being performed only for the independent layers instead of all the layers. The proposed methodology is an accurate method to reconstruct a 3D representative sample of a rock given only one 2D thin section. The algorithm is computationally efficient and faster than the previously introduced algorithm, and can easily be used to characterize samples for which 3D images are difficult to obtain in terms of both time and expense.

AB - High-resolution, 3D microstructure images of rocks can be used to compute the transport and elastic properties of those samples by use of a digital rock-physics approach. Those properties are complex functions of the pore-size distribution, geometry, and morphology, and necessitate the use of accurate 3D volumes. Because of the limited availability of 3D microstructure images of rocks, several attempts to construct 3D images from 2D images have been made. In this study, we propose a new stochastic method to reconstruct a 3D image of the rock by use of only a 2D section of the imagedrock sample. Our method is derived from a simple observation that the pores are gradually deformed from one section to the next. Therefore, the first step is to generate multiple independent realizations by performing multiple-point-statistics (MPS)-based stochastic simulations. These simulations represent independent 2D scans through the rock volume. Next, a succession of images is generated spanning two adjacent independent sections. These images consist of gradually morphed features from one section to those in the next independent section. Juxtaposing these 2D images results in the reconstructed 3D image. We implement the algorithm on a Berea sandstone rock for which the 3D high image was available for comparison. We calculate the spatial connectivity in the third direction, and confirmed that the proposed method can retrieve the connectivity in the third direction accurately. We also compute transport and elastic properties from the reconstructed image and from the original image to verify that this method reproduces the appropriate spatial statistics and pore-size distribution, geometry, and morphology. We also compare the numerical results with laboratory measurements performed on the sample. The results obtained by use of the reconstructed image reveal that the numerically calculated properties are similar to the measured values. We compare the mismatch of transport and elastic properties with the original measurements with that of the previous reconstruction algorithm. These comparisons show that the proposed simulation method has the same accuracy as previous ones. However, the proposed method is much more computationally efficient than the other algorithms, which are dependent on simulation of all 2D layers, mainly because of the faster MPS algorithm and the fact that the simulation is being performed only for the independent layers instead of all the layers. The proposed methodology is an accurate method to reconstruct a 3D representative sample of a rock given only one 2D thin section. The algorithm is computationally efficient and faster than the previously introduced algorithm, and can easily be used to characterize samples for which 3D images are difficult to obtain in terms of both time and expense.

UR - http://www.scopus.com/inward/record.url?scp=85019408314&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85019408314&partnerID=8YFLogxK

U2 - 10.2118/180489-PA

DO - 10.2118/180489-PA

M3 - Article

AN - SCOPUS:85019408314

VL - 20

SP - 342

EP - 352

JO - SPE Reservoir Evaluation and Engineering

JF - SPE Reservoir Evaluation and Engineering

SN - 1094-6470

IS - 2

ER -