Skip navigation
Skip navigation

Finite element modelling of the effective elastic properties of partially saturated rocks

Makarynska, Dina; Gurevich, Boris; Ciz, Radim; Arns, Christoph; Knackstedt, Mark

Description

Simulation of effective physical properties from microtomographic 3D images of porous structures allows one to relate properties of rocks directly to their microstructure. A static FEM code has been previously used to estimate effective elastic properties of fully saturated monomineralic (quartz) rock under wet and dry conditions. We use the code to calculate elastic properties under partially saturated conditions. The numerical predictions are compared to the Gassmann theory combined with...[Show more]

dc.contributor.authorMakarynska, Dina
dc.contributor.authorGurevich, Boris
dc.contributor.authorCiz, Radim
dc.contributor.authorArns, Christoph
dc.contributor.authorKnackstedt, Mark
dc.date.accessioned2015-12-10T22:12:18Z
dc.identifier.issn0098-3004
dc.identifier.urihttp://hdl.handle.net/1885/49586
dc.description.abstractSimulation of effective physical properties from microtomographic 3D images of porous structures allows one to relate properties of rocks directly to their microstructure. A static FEM code has been previously used to estimate effective elastic properties of fully saturated monomineralic (quartz) rock under wet and dry conditions. We use the code to calculate elastic properties under partially saturated conditions. The numerical predictions are compared to the Gassmann theory combined with Wood's formula (GW) for a mixture of pore fluids, which is exact for a monomineralic macroscopically homogeneous porous medium. Results of the numerical simulations performed for two Boolean sphere pack distributions show significant deviation from the GW limit and depend on the spatial distribution of fluids. This is shown to be a numerical artefact caused by incomplete equilibration of fluid pressure, which is primarily due to insufficient spatial resolution. To investigate the effect of pore-size and pore geometry, we perform FEM simulations for a model with regular pore geometry, where all pore channels have the same size and shape. Accuracy of these simulations increases with the total cross-section area of the channels and the size of individual channels. For the case where the total cross-section of the channels is large enough (on the same order as total porosity), there is a minimum of 4 voxels per channel diameter required for adequate fluid pressure equilibration throughout the pore space. Increasing the spatial resolution of the digital models reduces the discrepancy between the simulations and theory, but unfortunately increases the memory and CPU requirements of the simulations.
dc.publisherPergamon Press
dc.sourceComputers and Geosciences
dc.subjectKeywords: Computer simulation; Elasticity; Finite element method; Porosity; Porous materials; Pressure; Gassmann theory; Partial saturation; Partially saturated rocks; Pore fluids; Poroelasticity; Rocks; Computer simulation; Elasticity; Finite element method; Poros Effective elastic properties of rocks; Finite element method; Gassmann theory; Partial saturation; Poroelasticity
dc.titleFinite element modelling of the effective elastic properties of partially saturated rocks
dc.typeJournal article
local.description.notesImported from ARIES
local.identifier.citationvolume34
dc.date.issued2008
local.identifier.absfor040201 - Exploration Geochemistry
local.identifier.ariespublicationu9210271xPUB188
local.type.statusPublished Version
local.contributor.affiliationMakarynska, Dina, Curtin University of Technology
local.contributor.affiliationGurevich, Boris, Curtin University of Technology
local.contributor.affiliationCiz, Radim, CSIRO Petroleum
local.contributor.affiliationArns, Christoph, College of Physical and Mathematical Sciences, ANU
local.contributor.affiliationKnackstedt, Mark, College of Physical and Mathematical Sciences, ANU
local.description.embargo2037-12-31
local.bibliographicCitation.startpage647
local.bibliographicCitation.lastpage657
local.identifier.doi10.1016/j.cageo.2007.06.009
dc.date.updated2015-12-09T07:51:31Z
local.identifier.scopusID2-s2.0-40949096395
local.identifier.thomsonID000255472800007
CollectionsANU Research Publications

Download

File Description SizeFormat Image
01_Makarynska_Finite_element_modelling_of_2008.pdf1.25 MBAdobe PDF    Request a copy


Items in Open Research are protected by copyright, with all rights reserved, unless otherwise indicated.

Updated:  19 May 2020/ Responsible Officer:  University Librarian/ Page Contact:  Library Systems & Web Coordinator