Digital rock physics applied to squirt flow

Geophysics ◽  
2021 ◽  
pp. 1-40
Author(s):  
Simón Lissa ◽  
Matthias Ruf ◽  
Holger Steeb ◽  
Beatriz Quintal

We present a workflow for computing the seismic wave moduli dispersion and attenuation due to squirt flow in a numerical model derived from a micro X-Ray Computed Tomography image of cracked (through thermal treatment) Carrara marble sample. To generate the numerical model, the image is processed, segmented and meshed. The finite-element method is adopted to solve the linearized, quasi-static Navier-Stokes equations describing laminar flow of a compressible viscous fluid inside the cracks coupled with the quasi-static Lamé-Navier equations for the solid phase. We compute the effective P- and S-wave moduli in the three Cartesian directions for a model in dry conditions (saturated with air) and for a smaller model fully saturated with glycerin and having either drained or undrained boundary conditions. For the model saturated with glycerin, the results show significant and frequency-dependent P- and S-wave attenuation and the corresponding dispersion caused by squirt flow. Squirt flow occurs in response to fluid pressure gradients induced in the cracks by the imposed deformations. Our digital rock physics (DRP) workflow can be used to interpret laboratory measurements of attenuation using images of the rock sample.

2021 ◽  
Author(s):  
Simón Lissa ◽  
Matthias Ruf ◽  
Holger Steeb ◽  
Beatriz Quintal

<p>Seismic waves are affected by rock properties such as porosity, permeability, grain material and by their heterogeneities as well as by the fluid properties saturating the rocks. Consequently, seismic methods are a valuable tool for the indirect characterization of rocks. For example, at the microscale, the presence of compliant pores (cracks or grain contacts) in fluid-saturated rocks can cause strong seismic attenuation and velocity dispersion. In this case, the deformation caused by a passing wave induces a fluid pressure gradient between compressed compliant pores and much less compressed pores (stiff isometric pores or cracks having a different orientation than the most compressed ones) if they are hydraulically connected. The consequent fluid pressure diffusion (FPD) dissipates seismic energy due to viscous friction in the fluid.</p><p>Digital rock physics (DRP) aims to reproduce experimental measurements using numerical simulation in models derived from high resolution rock images. We developed a DRP workflow to calculate the frequency dependent seismic moduli dispersion and attenuation in fluid-saturated models derived from micro X-Ray Computed Tomography (µXRCT) images. Filtering, segmentation and meshing procedures are applied on sub-volumes of different rock images to create 3D numerical models. We apply our workflow to calculate seismic moduli attenuation due to FPD at the microscale (squirt flow). We consider a µXRCT image of a cracked (through thermal treatment) Carrara marble sample. A detailed visualization of the fluid pressure as well as of the energy dissipation rate in the 3D model helps to understand the squirt flow attenuation process at different frequencies.</p>


2021 ◽  
Vol 11 (5) ◽  
pp. 2113-2125
Author(s):  
Chenzhi Huang ◽  
Xingde Zhang ◽  
Shuang Liu ◽  
Nianyin Li ◽  
Jia Kang ◽  
...  

AbstractThe development and stimulation of oil and gas fields are inseparable from the experimental analysis of reservoir rocks. Large number of experiments, poor reservoir properties and thin reservoir thickness will lead to insufficient number of cores, which restricts the experimental evaluation effect of cores. Digital rock physics (DRP) can solve these problems well. This paper presents a rapid, simple, and practical method to establish the pore structure and lithology of DRP based on laboratory experiments. First, a core is scanned by computed tomography (CT) scanning technology, and filtering back-projection reconstruction method is used to test the core visualization. Subsequently, three-dimensional median filtering technology is used to eliminate noise signals after scanning, and the maximum interclass variance method is used to segment the rock skeleton and pore. Based on X-ray diffraction technology, the distribution of minerals in the rock core is studied by combining the processed CT scan data. The core pore size distribution is analyzed by the mercury intrusion method, and the core pore size distribution with spatial correlation is constructed by the kriging interpolation method. Based on the analysis of the core particle-size distribution by the screening method, the shape of the rock particle is assumed to be a more practical irregular polyhedron; considering this shape and the mineral distribution, the DRP pore structure and lithology are finally established. The DRP porosity calculated by MATLAB software is 32.4%, and the core porosity measured in a nuclear magnetic resonance experiment is 29.9%; thus, the accuracy of the model is validated. Further, the method of simulating the process of physical and chemical changes by using the digital core is proposed for further study.


Author(s):  
Mohammad Ebadi ◽  
Denis Orlov ◽  
Ivan Makhotin ◽  
Vladislav Krutko ◽  
Boris Belozerov ◽  
...  

Minerals ◽  
2019 ◽  
Vol 9 (11) ◽  
pp. 669
Author(s):  
Rongrong Lin ◽  
Leon Thomsen

With a detailed microscopic image of a rock sample, one can determine the corresponding 3-D grain geometry, forming a basis to calculate the elastic properties numerically. The issues which arise in such a calculation include those associated with image resolution, the registration of the digital numerical grid with the digital image, and grain anisotropy. Further, there is a need to validate the numerical calculation via experiment or theory. Because of the geometrical complexity of the rock, the best theoretical test employs the Hashin–Shtrikman result that, for an aggregate of two isotropic components with equal shear moduli, the bulk modulus is uniquely determined, independent of the micro-geometry. Similarly, for an aggregate of two isotropic components with a certain combination of elastic moduli defined herein, the Hashin–Shtrikman formulae give a unique result for the shear modulus, independent of the micro-geometry. For a porous, saturated rock, the solid incompressibility may be calculated via an “unjacketed” test, independent of the micro-geometry. Any numerical algorithm proposed for digital rock physics computation should be validated by successfully confirming these theoretical predictions. Using these tests, we validate a previously published staggered-grid finite difference damped time-stepping algorithm to calculate the static properties of digital rock models.


2022 ◽  
Author(s):  
Omar Alfarisi ◽  
Djamel Ouzzane ◽  
Mohamed Sassi ◽  
TieJun Zhang

<p><a></a>Each grid block in a 3D geological model requires a rock type that represents all physical and chemical properties of that block. The properties that classify rock types are lithology, permeability, and capillary pressure. Scientists and engineers determined these properties using conventional laboratory measurements, which embedded destructive methods to the sample or altered some of its properties (i.e., wettability, permeability, and porosity) because the measurements process includes sample crushing, fluid flow, or fluid saturation. Lately, Digital Rock Physics (DRT) has emerged to quantify these properties from micro-Computerized Tomography (uCT) and Magnetic Resonance Imaging (MRI) images. However, the literature did not attempt rock typing in a wholly digital context. We propose performing Digital Rock Typing (DRT) by: (1) integrating the latest DRP advances in a novel process that honors digital rock properties determination, while; (2) digitalizing the latest rock typing approaches in carbonate, and (3) introducing a novel carbonate rock typing process that utilizes computer vision capabilities to provide more insight about the heterogeneous carbonate rock texture.<br></p>


2021 ◽  
pp. 105008
Author(s):  
Eric J. Goldfarb ◽  
Ken Ikeda ◽  
Richard A. Ketcham ◽  
Maša Prodanović ◽  
Nicola Tisato

Sign in / Sign up

Export Citation Format

Share Document