Abstract
The use of 3D Computed Tomography (CT) throws light on the interior structure of the core samples. That is why the 3D CT has become increasingly widespread in digital rock physics. However, the CT does not provide information about the elastic properties of the rocks. So far, the most common approach to recover these properties is static laboratory measurements, when a sample is subject to a variety of static loads, followed by the measurement of displacements. This way is rather expensive, time-consuming and sometimes result in the destruction of a sample. Instead, we propose to implement numerical experiments by solving a 3D elastic static problem to perform upscaling of the CT data and to compute effective elastic parameters. In fact, to find these parameters we need to resolve a huge system of linear algebraic equations, which is a troublesome task. We have decided not to use a direct solver, but to apply a new iterative relaxation technique by considering a dynamic elastic problem with a special choice of relaxation parameters. The approach proposed needs parallelization. We have come to the conclusion to use the new feature of the latest Fortran extension known as CoArray Fortran (CAF) and to compare the three ways of parallel implementation as applied to this problem: MPI, MPI+OpenMP and CAF. Our experiments have proved that CAF and MPI approximately demonstrate the same performance, but the CAF possesses a clearer structure and is easy for coding.
Access this chapter
Tax calculation will be finalised at checkout
Purchases are for personal use only
References
Aboudi, J.: Mechanics of Composite Materials: A Unified Micromechanical Approach. Elsevier Science, Amsterdam (1991)
Arns, C.H., Knackstedt, M.A., Pinczewski, W.V., Garboczi, E.J.: Computation of linear elastic properties from microtomographic images: methodology and agreement between theory and experiment. Geophysics 67(5), P1395–P1405 (2002)
Chivers, I., Sleightholme, J.: Introduction to Programming with Fortran, 3rd edn. Springer, London (2015). https://doi.org/10.1007/b137984
Christensen, R.: Introduction to Mechanics of Composite Materials, 1st edn. Wiley, New York (1979)
Fredrich, J.T., DiGiovanni, A.A., Noble, D.R.: Predicting macroscopic transport properties using microscopic image data. J. Geophys. Res. 111, B03201 (2006). https://doi.org/10.1029/2005JB003774
ISO/IEC 1539–1:2010, Fortran - Part 1: Base language, International Standard (2010)
Keehm, Y.: Computationa lrock physics: transport properties in porous media and applications. Ph.D. dissertation. Stanford University (2003)
Kostin, V., Lisitsa, V., Reshetova, G., Tcheverda, V.: Parallel algorithm with modulus structure for simulation of seismic wave propagation in 3D multiscale multiphysics media. In: Malyshkin, V. (ed.) PaCT 2017. LNCS, vol. 10421, pp. 42–57. Springer, Cham (2017). https://doi.org/10.1007/978-3-319-62932-2_4
Kostin, V., Lisitsa, V., Reshetova, G., Tcheverda, V.: Simulation of seismic waves propagation in multiscale media: impact of cavernous/fractured reservoirs. In: Jónasson, K. (ed.) PARA 2010. LNCS, vol. 7133, pp. 54–64. Springer, Heidelberg (2012). https://doi.org/10.1007/978-3-642-28151-8_6
Landau, L.D., Lifshitz, E.M.: Mechanics. Nauka, Moscow (1988)
Reshetova, G., Khachkova, T.: Parallel numerical method to estimate the effective elastic moduli of rock core samples from 3D tomographic images. In: Dimov, I., Faragó, I., Vulkov, L. (eds.) FDM 2018. LNCS, vol. 11386, pp. 452–460. Springer, Cham (2019). https://doi.org/10.1007/978-3-030-11539-5_52
Saenger, E.H.: Numerical methods to determine effective elastic properties. Int. J. Eng. Sci. 46, 598–605 (2008)
Schoenberg, M., Muir, F.: A calculus for finely layered anisotropic media. Geophysics 54(5), 581–589 (1989)
Sendetski, J. (ed.): Composition Materials, vol. 2. Mechanics of Composition Meterials [Russian translation]. Mir, Moscow (1978)
Shermergor, T.: The theory of elasticity of microinhomogeneous media [Russian translation]. Nauka, Moscow (1977)
Zhang, W., Dai, G., Wang, F., Sun, S., Bassir, H.: Using strain energy-based prediction of effective elastic properties in topology optimization of material microstructures. Acta. Mech. Sin. 23(1), 77–89 (2007)
Vishnevsky, D., Lisitsa, V., Tcheverda, V., Reshetova, G.: Numerical study of the interface errors of finite-difference simulations of seismic waves. Geophysics 79, T219–T232 (2014)
Virieux, J.: P-SV wave propagation in heterogeneous media: velocity-stress finite-difference method. Geophysics 51, 889–901 (1986)
Zhan, X., Lawrence, M.S., Wave, S., Toksoz, M.N., Morgan, F.D.: Porescale modeling of electrical and fluid transport in Berea sandstone. Geophysics 75(5), F135–F142 (2010)
Joint Supercomputer Center of RAS. http://old.jscc.ru/eng/index.shtml
Novosibirsk Supercomputer Center of SB RAS. http://www.sscc.icmmg.nsc.ru
Irkutsk Supercomputer Center of SB RAS. http://hpc.icc.ru
Acknowledgements
The reported study was funded by RFBR according to the research project â„–19-01-00347.
The research is carried out using the equipment of the shared research facilities of HPC computing resources at the Joint Supercomputer Center of RAS [20], the Siberian Supercomputer Center [21] and Irkutsk Supercomputer Center of SB RAS [22].
Author information
Authors and Affiliations
Corresponding author
Editor information
Editors and Affiliations
Rights and permissions
Copyright information
© 2019 Springer Nature Switzerland AG
About this paper
Cite this paper
Reshetova, G., Cheverda, V., Khachkova, T. (2019). Numerical Experiments with Digital Twins of Core Samples for Estimating Effective Elastic Parameters. In: Voevodin, V., Sobolev, S. (eds) Supercomputing. RuSCDays 2019. Communications in Computer and Information Science, vol 1129. Springer, Cham. https://doi.org/10.1007/978-3-030-36592-9_24
Download citation
DOI: https://doi.org/10.1007/978-3-030-36592-9_24
Published:
Publisher Name: Springer, Cham
Print ISBN: 978-3-030-36591-2
Online ISBN: 978-3-030-36592-9
eBook Packages: Computer ScienceComputer Science (R0)