scholarly journals Modeling compaction effects on the elastic properties of clay-water composites

Geophysics ◽  
2012 ◽  
Vol 77 (5) ◽  
pp. D171-D183 ◽  
Author(s):  
Bernardo Moyano ◽  
Kyle T. Spikes ◽  
Tor Arne Johansen ◽  
Nazmul Haque Mondol

Modeling the elastic properties of clay-bearing rocks (shales) requires thorough knowledge of the mineral constituents, their elastic properties, pore space microstructure, and orientations of clay platelets. Information about these variables and their complex interrelationships is rarely available for real rocks. We theoretically modeled the elastic properties of synthetic clay-water composites compacted in the laboratory, including estimates of pore space topology and percolation behavior. The mineralogy of the samples was known exactly, and the focus was on two monomineralic samples comprised of kaolinite and smectite. We used differential effective medium theory (DEM) and analysis of scanning electron microscope (SEM) images of the compacted kaolinite and smectite samples. Percolation behavior was included through calculations of critical porosities from measurements of the liquid limits of the individual clay powders. Quantitative analysis of the SEM images showed that the large scale ([Formula: see text]) pore space of the smectite composite had more rounded pores (mean aspect ratio [Formula: see text]) than the kaolinite composite (mean pore’s aspect ratio [Formula: see text]). However, models that used only these large-scale pore shapes could not explain the compressional and shear velocity measurements. DEM simulations with a single pore aspect ratio showed that bulk and shear moduli are controlled by different pore shapes. Conversely, modeling results that combined critical porosity and dual porosity models into DEM theory compared well with the measured bulk and shear moduli of compacting kaolinite and smectite composites. The methods and results we used could be used to model unconsolidated clay-bearing rocks of more complex mineralogy.

Geophysics ◽  
2012 ◽  
Vol 77 (3) ◽  
pp. WA157-WA168 ◽  
Author(s):  
Osni Bastos de Paula ◽  
Marina Pervukhina ◽  
Dina Makarynska ◽  
Boris Gurevich

Modeling dispersion and attenuation of elastic waves in fluid-saturated rocks due to squirt flow requires the knowledge of a number of geometrical parameters of the pore space, in particular, the characteristic aspect ratio of the pores. These parameters are usually inferred by fitting measurements on saturated rocks to model predictions. To eliminate such fitting and thus make the model more predictive, we propose to recover the geometrical parameters of the pore space from the pressure dependency of elastic moduli on dry samples. Our analysis showed that the pressure dependency of elastic properties of rocks (and their deviation from Gassmann’s prediction) at ultrasonic frequencies is controlled by the squirt flow between equant, stiff, and so-called intermediate pores (with aspect ratios between [Formula: see text]). Such intermediate porosity is expected to close at confining pressures of between 200 and 2000 MPa, and thus cannot be directly obtained from ultrasonic experiments performed at pressures below 50 MPa. However, the presence of this intermediate porosity is inferred from the significant linear trend in the pressure dependency of elastic properties of the dry rock and the difference between the bulk modulus of the dry rock computed for spherical pores and the measured modulus at 50 MPa. Moreover, we can infer the magnitude of the intermediate porosity and its characteristic aspect ratio. Substituting these parameters into the squirt model, we have computed elastic moduli and velocities of the water-saturated rock and compared these predictions against laboratory measurements of these velocities. The agreement is good for a number of clean sandstones, but not unexpectedly worse for a broad range of shaley sandstones. Our predictions showed that dispersion and attenuation caused by the squirt flow between compliant and stiff pores may occur in the seismic frequency band. Confirmation of this prediction requires laboratory measurements of elastic properties at these frequencies.


2021 ◽  
Author(s):  
Michele Paulatto ◽  
Joanna Morgan ◽  
Kajetan Chrapkiewicz ◽  
Emilie Hooft ◽  
Doug Toomey ◽  
...  

<p>The lack of direct seismological evidence for large molten magma chambers is considered to be one of the most important arguments in support of the mush paradigm. However, most published melt fraction estimates based on interpretation of seismological data are associated with large uncertainties because of two limitations: i) inherent limits to resolution of seismic tomography and, ii) trade-offs in the constitutive relationships that tie seismic properties to melt fraction. Low-velocity volumes associated with magma storage are particularly difficult to image with conventional travel-time tomography due to limited resolution and wavefront healing, resulting in blurred images and a high velocity bias. We tackle these limitation by applying full waveform inversion to active source seismic data collected over the Kolumbo submarine volcano (Greece). We recover a previously undetected Vp anomaly of –50% beneath the volcano and interpret this as a shallow magmatic intrusion. Extension of this approach to the wider Santorini volcanic system is ongoing. Concurrently, we are tackling the second limitation, which is the result of the dependence of elastic properties on the microgeometry of the melt. Seismological melt estimates rely on the assumption that the melt pore space can be represented by simple geometrical shapes, usually ellipsoids, with a given aspect ratio. Since the aspect ratio is poorly constrained, this results in a trade-off between melt fraction and melt geometry. We have adapted a method for the homogenisation of the elastic properties of multi-phase composites and applied it to calculating the elastic properties of partially molten rocks starting from the melt microstructure determined by X-ray CT scanning. The microgeometry of the mush can be inferred from the study of glomerocrysts: crystal mush inclusions with quenched interstitial melt that are carried to the surface by erupted lava. After the sample is digitized and segmented into its constitutive phases (crystals, melt, vesicles), the average elastic properties are determined by numerical homogenisation which consists of numerically simulating the deformation of the sample under load and predicting its elastic response. The results are compared to a semi analytical solution for ellipsoidal inclusions. We apply this approach to a plutonic nodule from St Kitts and show that the melt microstructure leads to an elastic response equivalent to that of ellipsoidal melt inclusions with an aspect ratio of 0.1 (oblate spheroids). This equivalent aspect ratio is used to refine melt estimates for Montserrat, Santorini and Kolumbo volcano.</p>


Author(s):  
C. A. Callender ◽  
Wm. C. Dawson ◽  
J. J. Funk

The geometric structure of pore space in some carbonate rocks can be correlated with petrophysical measurements by quantitatively analyzing binaries generated from SEM images. Reservoirs with similar porosities can have markedly different permeabilities. Image analysis identifies which characteristics of a rock are responsible for the permeability differences. Imaging data can explain unusual fluid flow patterns which, in turn, can improve production simulation models.Analytical SchemeOur sample suite consists of 30 Middle East carbonates having porosities ranging from 21 to 28% and permeabilities from 92 to 2153 md. Engineering tests reveal the lack of a consistent (predictable) relationship between porosity and permeability (Fig. 1). Finely polished thin sections were studied petrographically to determine rock texture. The studied thin sections represent four petrographically distinct carbonate rock types ranging from compacted, poorly-sorted, dolomitized, intraclastic grainstones to well-sorted, foraminiferal,ooid, peloidal grainstones. The samples were analyzed for pore structure by a Tracor Northern 5500 IPP 5B/80 image analyzer and a 80386 microprocessor-based imaging system. Between 30 and 50 SEM-generated backscattered electron images (frames) were collected per thin section. Binaries were created from the gray level that represents the pore space. Calculated values were averaged and the data analyzed to determine which geological pore structure characteristics actually affect permeability.


2020 ◽  
Vol 6 (1) ◽  
pp. 50-56
Author(s):  
Francesco Baino ◽  
Elisa Fiume

AbstractPorosity is known to play a pivotal role in dictating the functional properties of biomedical scaffolds, with special reference to mechanical performance. While compressive strength is relatively easy to be experimentally assessed even for brittle ceramic and glass foams, elastic properties are much more difficult to be reliably estimated. Therefore, describing and, hence, predicting the relationship between porosity and elastic properties based only on the constitutive parameters of the solid material is still a challenge. In this work, we quantitatively compare the predictive capability of a set of different models in describing, over a wide range of porosity, the elastic modulus (7 models), shear modulus (3 models) and Poisson’s ratio (7 models) of bioactive silicate glass-derived scaffolds produced by foam replication. For these types of biomedical materials, the porosity dependence of elastic and shear moduli follows a second-order power-law approximation, whereas the relationship between porosity and Poisson’s ratio is well fitted by a linear equation.


Metals ◽  
2021 ◽  
Vol 11 (6) ◽  
pp. 968
Author(s):  
Fumitada Iguchi ◽  
Keisuke Hinata

The elastic properties of 0, 10, 15, and 20 mol% yttrium-doped barium zirconate (BZY0, BZY10, BZY15, and BZY20) at the operating temperatures of protonic ceramic fuel cells were evaluated. The proposed measurement method for low sinterability materials could accurately determine the sonic velocities of small-pellet-type samples, and the elastic properties were determined based on these velocities. The Young’s modulus of BZY10, BZY15, and BZY20 was 224, 218, and 209 GPa at 20 °C, respectively, and the values decreased as the yttrium concentration increased. At high temperatures (>20 °C), as the temperature increased, the Young’s and shear moduli decreased, whereas the bulk modulus and Poisson’s ratio increased. The Young’s and shear moduli varied nonlinearly with the temperature: The values decreased rapidly from 100 to 300 °C and gradually at temperatures beyond 400 °C. The Young’s modulus of BZY10, BZY15, and BZY20 was 137, 159, and 122 GPa at 500 °C, respectively, 30–40% smaller than the values at 20 °C. The influence of the temperature was larger than that of the change in the yttrium concentration.


Geophysics ◽  
2013 ◽  
Vol 78 (3) ◽  
pp. D169-D179 ◽  
Author(s):  
Zijian Zhang ◽  
De-hua Han ◽  
Daniel R. McConnell

Hydrate-bearing sands and shallow nodular hydrate are potential energy resources and geohazards, and they both need to be better understood and identified. Therefore, it is useful to develop methodologies for modeling and simulating elastic constants of these hydrate-bearing sediments. A gas-hydrate rock-physics model based on the effective medium theory was successfully applied to dry rock, water-saturated rock, and hydrate-bearing rock. The model was used to investigate the seismic interpretation capability of hydrate-bearing sediments in the Gulf of Mexico by computing elastic constants, also known as seismic attributes, in terms of seismic interpretation, including the normal incident reflectivity (NI), Poisson’s ratio (PR), P-wave velocity ([Formula: see text]), S-wave velocity ([Formula: see text]), and density. The study of the model was concerned with the formation of gas hydrate, and, therefore, hydrate-bearing sediments were divided into hydrate-bearing sands, hydrate-bearing sands with free gas in the pore space, and shallow nodular hydrate. Although relations of hydrate saturation versus [Formula: see text] and [Formula: see text] are different between structures I and II gas hydrates, highly concentrated hydrate-bearing sands may be interpreted on poststack seismic amplitude sections because of the high NI present. The computations of elastic constant implied that hydrate-bearing sands with free gas could be detected with the crossplot of NI and PR from prestack amplitude analysis, and density may be a good hydrate indicator for shallow nodular hydrate, if it can be accurately estimated by seismic methods.


2021 ◽  
Author(s):  
Vladislav Vasilevich Alekseev ◽  
Denis Mihaylovich Orlov ◽  
Dmitry Anatolevich Koroteev

Abstract The approaches of building and methods of using the digital core are currently developing rapidly. The use of these methods makes it possible to obtain petrophysical information by non-destructive methods quickly. Digital rock physics includes two main stages: constructing models and modeling various physical processes on the obtained models. Our work proposes using deep learning methods for mineral and pore space segmentation instead of classical methods such as threshold image processing. Deep neural networks have long been able to show their advantages in many areas of computer vision. This paper proposes and tests methods that help identify different minerals in images from a scanning electron microscope. We used images of rocks of the Achimov formation, which are arkoses, as samples. We tested various deep neural networks such as LinkNet, U-Net, ResUNet, and pix2pix and identified those that performed best in segmentation.


Sign in / Sign up

Export Citation Format

Share Document