scholarly journals Numerical modeling of wave processes in fractured porous fluid-saturated media

Author(s):  
М.А. Новиков ◽  
В.В. Лисица ◽  
А.А. Козяев

Одной из актуальных задач современной прикладной геофизики является выделение характерных признаков наличия развитой трещиноватости в пласте по сейсмическим данным. Более того, необходимо выделять флюидозаполненные системы трещин, образующих связанную систему трещин, способную обеспечивать достаточную гидродинамическую проницаемость резервуара. В настоящей статье представлен численный алгоритм расчета волновых полей в трещиноватых пороупругих средах, основанный на конечно-разностной аппроксимации уравнений Био. На основе численных экспериментов показано, что связность трещин, т.е. наличие систем пересекающихся трещин, существенно повышает поглощение сейсмической энергии, что обусловлено возникновением локальных потоков флюида внутри трещин. Приводится детальный частотный анализ затухания сейсмических волн и обусловленной этим дисперсии. One of important problems in modern applied geophysics is to distinguish the specific indications of developed fractures in geological formations using seismic data. In addition, it is necessary to distinguish the fluid-saturated systems of cracks capable of providing a sufficient hydrodynamic permeability of reservoirs. In this paper we propose a numerical algorithm to study wave fields in fractured porous fluid-saturated media on the basis of a finite-difference approximation of Biot's equations. Our numerical results show that the intersecting fractures significantly increase the absorption of seismic energy due to the appearance of fluid flows in cracks. A detailed frequency analysis of seismic wave attenuation is given.

Geophysics ◽  
1981 ◽  
Vol 46 (5) ◽  
pp. 806-808 ◽  
Author(s):  
N. C. Dutta

In the petroleum industry there is considerable current interest in understanding the causes of seismic energy loss. During the past few years, considerable progress has been made in the field of seismic wave attenuation as a result of both controlled laboratory studies of the static and the dynamic properties of rock specimens and theoretical modeling. During the 49th Annual International Meeting of the Society of Exploration Geophysicists in New Orleans, a research workshop on Seismic Loss Mechanisms was organized. The purpose of this workshop was to assess current understanding of the physical processes that cause attenuation of seismic energy. The workshop, organized by Kenneth Larner of Western Geophysical and the author, featured a multispeaker format involving invited papers from both industry and academia. The presentations were followed by discussions lasting approximately 90 minutes which included the speakers, the audience, and a panel of experts in the field of seismology.


Geophysics ◽  
1989 ◽  
Vol 54 (4) ◽  
pp. 524-527 ◽  
Author(s):  
R. Mörig ◽  
H. Burkhardt

Seismic wave attenuation has been a subject of interest during the last 40 years because it may be of use in interpreting seismic data. From this attenuation parameter, more detailed information about the lithology of the subsurface may be deduced if we understand the absorption mechanisms by which dissipation of seismic energy is governed. We are, therefore, studying in the laboratory the effects of different parameters such as porosity, permeability, pore fluid, and saturation state on the absorption of seismic waves in porous rocks over a wide spectrum ranging from seismic to ultrasonic frequencies (Burkhardt et al., 1986).


2021 ◽  
Vol 2 (2) ◽  
pp. 186-195
Author(s):  
Mikhail A. Novikov ◽  
Vadim V. Lisitsa

In our work we investigate the effect of transport and elastic properties anisotropy on seismic attenuation due to fracture-to-fracture wave-induced fluid flow using numerical algorithm of estimation of seismic wave attenuation in anisotropic fractured porous fluid-saturated media. Algorithm is based on numerical solution of anisotropic Biot equations using finite-difference scheme on staggered grid. We perform a set of numerical experiments to model wave propagation in fractured media with anisotropic fractured-filling material providing wave-induced fluid flow within interconnected fractures. Recorded signals are used for numerical estimation of inverse quality factor. Results demonstrate the effect of fracture-filling material anisotropy on seismic wave attenuation.


2018 ◽  
Vol 2018 ◽  
pp. 1-13 ◽  
Author(s):  
Mingwei Zhang ◽  
Shengdong Liu ◽  
Shuzhao Chen ◽  
Yanlong Chen ◽  
Guang Xu ◽  
...  

Based on the energy attenuation characteristics of residual wave in deep rock, a method was developed to determine the microseismic focus energy. Differential energy loss in infinitesimal spreading distance is logically deduced, upon which energy attenuation equation was established. With a logarithmic transformation, a linear relation of the residual seismic energy with distance is formulated. Its intercept was used to determine the microseismic focus energy. The result is compared with that determined by the energy density method. The reliability of the determined focus energy and the impact of the built-in velocity threshold on the residual wave energy computation are discussed. Meanwhile, the energy absorption coefficient used for representing the absorption characteristics of the rock medium in the mining region under study is also clarified. Key findings show that the microseismic focus energy confirmed by the residual wave attenuation is reliable. The result’s accuracy is quite high, especially for the events in deep rock with great homogeneity. The developed focus energy computation method is closely dependent on the integrity of waveform, accuracy of repositioning, and reliability of effective components extraction. The new method has been shown to be effective and practical.


2012 ◽  
Vol 12 (1) ◽  
pp. 193-225 ◽  
Author(s):  
N. Anders Petersson ◽  
Björn Sjögreen

AbstractWe develop a stable finite difference approximation of the three-dimensional viscoelastic wave equation. The material model is a super-imposition of N standard linear solid mechanisms, which commonly is used in seismology to model a material with constant quality factor Q. The proposed scheme discretizes the governing equations in second order displacement formulation using 3N memory variables, making it significantly more memory efficient than the commonly used first order velocity-stress formulation. The new scheme is a generalization of our energy conserving finite difference scheme for the elastic wave equation in second order formulation [SIAM J. Numer. Anal., 45 (2007), pp. 1902-1936]. Our main result is a proof that the proposed discretization is energy stable, even in the case of variable material properties. The proof relies on the summation-by-parts property of the discretization. The new scheme is implemented with grid refinement with hanging nodes on the interface. Numerical experiments verify the accuracy and stability of the new scheme. Semi-analytical solutions for a half-space problem and the LOH.3 layer over half-space problem are used to demonstrate how the number of viscoelastic mechanisms and the grid resolution influence the accuracy. We find that three standard linear solid mechanisms usually are sufficient to make the modeling error smaller than the discretization error.


1966 ◽  
Vol 6 (03) ◽  
pp. 217-227 ◽  
Author(s):  
Hubert J. Morel-Seytoux

Abstract The influence of pattern geometry on assisted oil recovery for a particular displacement mechanism is the object of investigation in this paper. The displacement is assumed to be of unit mobility ratio and piston-like. Fluids are assumed incompressible and gravity and capillary effects are neglected. With these assumptions it is possible to calculate by analytical methods the quantities of interest to the reservoir engineer for a great variety of patterns. Specifically, this paper presentsvery briefly, the methods and mathematical derivations required to obtain the results of engineering concern, andtypical results in the form of graphs or formulae that can be used readily without prior study of the methods. Results of this work provide checks for solutions obtained from programmed numerical techniques. They also reveal the effect of pattern geometry and, even though the assumptions of piston-like displacement and of unit mobility ratio are restrictive, they can nevertheless be used for rather crude but quick, cheap estimates. These estimates can be refined to account for non-unit mobility ratio and two-phase flow by correlating analytical results in the case M=1 and the numerical results for non-Piston, non-unit mobility ratio displacements. In an earlier paper1 it was also shown that from the knowledge of closed form solutions for unit mobility ratio, quantities called "scale factors" could be readily calculated, increasing considerably the flexibility of the numerical techniques. Many new closed form solutions are given in this paper. INTRODUCTION BACKGROUND Pattern geometry is a major factor in making water-flood recovery predictions. For this reason many numerical schemes have been devised to predict oil recovery in either regular patterns or arbitrary configurations. The numerical solutions, based on the method of finite difference approximation, are subject to errors often difficult to evaluate. An estimate of the error is possible by comparison with exact solutions. To provide a variety of checks on numerical solutions, a thorough study of the unit mobility ratio displacement process was undertaken. To calculate quantities of interest to the reservoir engineer (oil recovery, sweep efficiency, etc.), it is necessary to first know the pressure distribution in the pattern. Then analytical procedures are used to calculate, in order of increasing difficulty: injectivity, breakthrough areal sweep efficiency, normalized oil recovery and water-oil ratio as a function of normalized PV injected. BACKGROUND Pattern geometry is a major factor in making water-flood recovery predictions. For this reason many numerical schemes have been devised to predict oil recovery in either regular patterns or arbitrary configurations. The numerical solutions, based on the method of finite difference approximation, are subject to errors often difficult to evaluate. An estimate of the error is possible by comparison with exact solutions. To provide a variety of checks on numerical solutions, a thorough study of the unit mobility ratio displacement process was undertaken. To calculate quantities of interest to the reservoir engineer (oil recovery, sweep efficiency, etc.), it is necessary to first know the pressure distribution in the pattern. Then analytical procedures are used to calculate, in order of increasing difficulty: injectivity, breakthrough areal sweep efficiency, normalized oil recovery and water-oil ratio as a function of normalized PV injected.


Sign in / Sign up

Export Citation Format

Share Document