Finite-difference modeling of wave propagation and diffusion in poroelastic media

Geophysics ◽  
2009 ◽  
Vol 74 (4) ◽  
pp. T55-T66 ◽  
Author(s):  
Fabian Wenzlau ◽  
Tobias M. Müller
Keyword(s):  
Numerical Modeling ◽  
Finite Difference ◽  
Diffusion Processes ◽  
Spatial Scales ◽  
Rock Physics ◽  
Seismic Attenuation ◽  
Rock Properties ◽  
Modeling Tools ◽  
Reservoir Rocks ◽  
Poroelastic Media

Numerical modeling of seismic waves in heterogeneous, porous reservoir rocks is an important tool for interpreting seismic surveys in reservoir engineering. Various theoretical studies derive effective elastic moduli and seismic attributes from complex rock properties, involving patchy saturation and fractured media. To confirm and further develop rock-physics theories for reservoir rocks, accurate numerical modeling tools are required. Our 2D velocity-stress, finite-difference scheme simulates waves within poroelastic media as described by Biot’s theory. The scheme is second order in time, contains high-order spatial derivative operators, and is parallelized using the domain-decomposition technique. A series of numerical experiments that are compared to exact analytical solutions allow us to assess the stability conditions and dispersion relations of the explicit poroelastic finite-differ-ence method. The focus of the experiments is to model wave-induced flow accurately in the vicinity of mesoscopic heterogeneities such as cracks and gas inclusions in partially saturated rocks. For that purpose, a suitable numerical setup is applied to extract seismic attenuation and dispersion from quasi-static experiments. Our results confirm that finite-difference modeling is a valuable tool to simulate wave propa-gation in heterogeneous poroelastic media, provided the temporal and spatial scales of the propagating waves and of the induced fluid-diffusion processes are resolved properly.

Numerical Modeling of Wave Propagation in Poroelastic Media with a Staggered-Grid Finite-Difference Method

2013 ◽  
Vol 275-277 ◽  
pp. 612-617
Author(s):  
Wen Sheng Zhang ◽  
Li Tong
Keyword(s):  
Numerical Modeling ◽  
Wave Propagation ◽  
Finite Difference ◽  
Finite Difference Method ◽  
Difference Method ◽  
Staggered Grid ◽  
Poroelastic Medium ◽  
Numerical Computations ◽  
Poroelastic Media ◽  
Biot’S Equations

In this paper, wave propagation in poroelastic medium is simulated with a staggered-grid finite-difference method. The formulation is discretized based on the second-order Biot’s equations rather than the corresponding velocity-stress form. In order to eliminate boundary reflections, the PML method is applied. Numerical computations are implemented and the results show the correctness and effectiveness of the schemes presented in this paper.

Seismic moduli dispersion and attenuation obtained using Digital Rock Physics

2021 ◽  
Author(s):  
Simón Lissa ◽  
Matthias Ruf ◽  
Holger Steeb ◽  
Beatriz Quintal
Keyword(s):  
Seismic Waves ◽  
Numerical Models ◽  
Fluid Pressure ◽  
Rock Physics ◽  
Viscous Friction ◽  
Seismic Attenuation ◽  
Rock Properties ◽  
Digital Rock Physics ◽  
Digital Rock ◽  
Dispersion And Attenuation

<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>

scholarly journals Elastic-Electrical Rock-Physics Template for the Characterization of Tight-Oil Reservoir Rocks

Lithosphere ◽  
2021 ◽  
Vol 2021 (Special 3) ◽  
Author(s):  
Mengqiang Pang ◽  
Jing Ba ◽  
José M. Carcione ◽  
Erik H. Saenger
Keyword(s):  
Complex Structure ◽  
Rock Physics ◽  
Clay Content ◽  
Rock Properties ◽  
Tight Oil ◽  
Log Data ◽  
Reservoir Rocks ◽  
Porosity And Permeability ◽  
Parallel Network ◽  
Log Interpretation

Abstract Tight-oil reservoirs have low porosity and permeability, with microcracks, high clay content, and a complex structure resulting in strong heterogeneities and poor connectivity. Thus, it is a challenge to characterize this type of reservoir with a single geophysical methodology. We propose a dual-porosity-clay parallel network to establish an electrical model and the Hashin-Shtrikman and differential effective medium equations to model the elastic properties. Using these two models, we compute the rock properties as a function of saturation, clay content, and total and microcrack porosities. Moreover, a 3D elastic-electrical template, based on resistivity, acoustic impedance, and Poisson’s ratio, is built. Well-log data is used to calibrate the template. We collect rock samples and log data (from two wells) from the Songliao Basin (China) and analyze their microstructures by scanning electron microscopy. Then, we study the effects of porosity and clay content on the elastic and electrical properties and obtain a good agreement between the predictions, log interpretation, and actual production reports.

scholarly journals On the generation of internal waves by river plumes in subcritical initial conditions

2021 ◽  
Vol 11 (1) ◽  
Author(s):  
R. Mendes ◽  
J. C. B. da Silva ◽  
J. M. Magalhaes ◽  
B. St-Denis ◽  
D. Bourgault ◽  
...  
Keyword(s):  
Numerical Modeling ◽  
Internal Waves ◽  
Coastal Waters ◽  
Initial Conditions ◽  
Spatial Scales ◽  
River Plumes ◽  
Wide Range ◽  
Near Shore ◽  
Generation Mechanisms ◽  
Douro River

AbstractInternal waves (IWs) in the ocean span across a wide range of time and spatial scales and are now acknowledged as important sources of turbulence and mixing, with the largest observations having 200 m in amplitude and vertical velocities close to 0.5 m s−1. Their origin is mostly tidal, but an increasing number of non-tidal generation mechanisms have also been observed. For instance, river plumes provide horizontally propagating density fronts, which were observed to generate IWs when transitioning from supercritical to subcritical flow. In this study, satellite imagery and autonomous underwater measurements are combined with numerical modeling to investigate IW generation from an initial subcritical density front originating at the Douro River plume (western Iberian coast). These unprecedented results may have important implications in near-shore dynamics since that suggest that rivers of moderate flow may play an important role in IW generation between fresh riverine and coastal waters.

Digital Rock Typing DRT Algorithm Formulation with Optimal Supervised Semantic Segmentation

2022 ◽  
Author(s):  
Omar Alfarisi ◽  
Djamel Ouzzane ◽  
Mohamed Sassi ◽  
TieJun Zhang
Keyword(s):  
Carbonate Rock ◽  
Chemical Properties ◽  
Rock Physics ◽  
Semantic Segmentation ◽  
Rock Properties ◽  
Digital Rock Physics ◽  
Digital Rock ◽  
Rock Types ◽  
Fluid Saturation ◽  
Rock Typing

<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>

Sign in / Sign up

Export Citation Format

Share Document