Elastic inverse scattering for fluid variation with time-lapse seismic data

Geophysics ◽  
2015 ◽  
Vol 80 (2) ◽  
pp. WA61-WA67 ◽  
Author(s):  
Zhaoyun Zong ◽  
Xingyao Yin ◽  
Guochen Wu ◽  
Zhiping Wu

Elastic inverse-scattering theory has been extended for fluid discrimination using the time-lapse seismic data. The fluid factor, shear modulus, and density are used to parameterize the reference medium and the monitoring medium, and the fluid factor works as the hydrocarbon indicator. The baseline medium is, in the conception of elastic scattering theory, the reference medium, and the monitoring medium is corresponding to the perturbed medium. The difference in the earth properties between the monitoring medium and the baseline medium is taken as the variation in the properties between the reference medium and perturbed medium. The baseline and monitoring data correspond to the background wavefields and measured full fields, respectively. And the variation between the baseline data and monitoring data is taken as the scattered wavefields. Under the above hypothesis, we derived a linearized and qualitative approximation of the reflectivity variation in terms of the changes of fluid factor, shear modulus, and density with the perturbation theory. Incorporating the effect of the wavelet into the reflectivity approximation as the forward solver, we determined a practical prestack inversion approach in a Bayesian scheme to estimate the fluid factor, shear modulus, and density changes directly with the time-lapse seismic data. We evaluated the examples revealing that the proposed approach rendered the estimation of the fluid factor, shear modulus, and density changes stably, even with moderate noise.

Geophysics ◽  
1994 ◽  
Vol 59 (6) ◽  
pp. 963-972 ◽  
Author(s):  
Bastian Blonk ◽  
Gérard C. Herman

A method is presented for eliminating near‐surface scattered noise from seismic data. Starting from an appropriately chosen background model, a surface‐consistent scattering model is determined using linearized elastodynamic inverse scattering theory. This scattering model does not necessarily equal the actual scatterer distribution, but it enables one to calculate, approximately, the near‐surface scattered part of the data. The method honors at least some of the complexity of the near‐surface scattering process and can be applied in cases where traditional methods, like wavenumber‐frequency filtering techniques and methods for static corrections, are ineffective. From a number of tests on synthetic data, we conclude that the method is rather robust; its main sensitivity is because of errors in the determination of the background Rayleigh‐wave velocity.


2000 ◽  
Vol 3 (06) ◽  
pp. 517-524
Author(s):  
J.K. Hughes

Summary The propagation of elastic waves in rocks is determined by the bulk modulus, shear modulus, and bulk density of the rock. In porous rocks all these properties are affected by the distribution of pore space, the geometry and interconnectivity of the pores, and the nature of the fluid occupying the pore space. In addition, the bulk and shear moduli are also affected by the effective pressure, which is equivalent to the difference between the confining (or lithostatic) pressure and pore pressure. During production of hydrocarbons from a reservoir, the movement of fluids and changes in pore pressure may contribute to a significant change in the elastic moduli and bulk density of the reservoir rocks. This phenomenon is the basis for reservoir monitoring by repeated seismic (or time-lapse) surveys whereby the difference in seismic response during the lifetime of the field can be directly related to changes in the pore fluids and/or pore pressure. Under suitable conditions, these changes in the reservoir during production can be quantitatively estimated by appropriate repeat three-dimensional (3D) seismic surveys which can contribute to understanding of the reservoir model away from the wells. The benefit to reservoir management is a better flow model which incorporates the information derived from the seismic data. What are suitable conditions? There are two primary factors which determine whether the reservoir changes we wish to observe will be detectable in the seismic data:the magnitude of the change in the elastic moduli (and bulk density) of the reservoir rocks as a result of fluid displacement, pressure changes, etc.;the magnitude of the repeatability errors between time-lapse seismic surveys. This includes errors associated with seismic data collection, ambient noise and data processing. The first is the signal component and the second the noise component. Previous reviews of seismic monitoring suggest that for 3D seismic surveys a signal-to-noise (S/N) ratio of 1.0 is sufficient for qualitative estimation of reservoir changes. Higher S/N ratios may allow quantitative estimates. After a brief examination of the rock physics affecting the seismic signal, we examine the second factor, repeatability errors, and use a synthetic seismic model to illustrate some of the factors which contribute to repeatability error. We also use two land 3D surveys over a Middle East carbonate reservoir to illustrate seismic repeatability. The study finds that repeatability errors, while always larger than desired, are generally within limits which will allow production-induced changes in seismic reflectivity to be confidently detected. Introduction Seismic data have been used successfully for many decades in the petroleum industry and have contributed significantly to the discovery of new fields throughout the world. Initially, seismic surveys were primarily an exploration tool, assisting in the identification of potential hydrocarbon structural and stratigraphic traps for drilling targets. With the introduction of 3D seismic surveys in the 1970's, accurate geological structural mapping became possible while the use of new seismic attributes as hydrocarbon indicators improved the success rate of discovery wells. More recently seismic data have also contributed to a better reservoir description away from the wells by making use of the correlation between suitable seismic attributes and petrophysical quantities such as porosity and net to gross, and by incorporating robust geostatistical methods for estimating the static reservoir model. Better seismic acquisition technology, improved seismic processing methods and an overall improvement in signal to noise have led to further 3D seismic surveys over producing fields primarily for better imaging of the reservoir and improved reservoir characterization. The concept of using repeated seismic surveys (time-lapse seismic) for monitoring changes in the reservoir due to production was suggested in the 1980's,1-3 and early tests were done by Arco in the Holt Sand fireflood4 from 1981-83. Over the last few years, the number of publications relating to time-lapse seismic [often referred to as four-dimensional (4D) seismic] has increased dramatically. Prior to time-lapse seismic monitoring, seismic data have been the domain of geologists and geophysicists, but the possibility of monitoring fluid displacements and pressure changes in a producing reservoir, away from the wells, has direct relevance to reservoir engineers and reservoir management. More exciting possibilities have been introduced by the use of time-lapse seismic data in combination with production history matching5 for greater refinement in optimization of the reservoir model. It is important, however, that reliable criteria are used to assess the feasibility of seismic monitoring.6


Geophysics ◽  
2007 ◽  
Vol 72 (2) ◽  
pp. G21-G29 ◽  
Author(s):  
Timothy C. Johnson ◽  
Partha S. Routh ◽  
Warren Barrash ◽  
Michael D. Knoll

Ground-penetrating radar (GPR) attenuation-difference tomography is a useful tool for imaging the migration of electrically anomalous tracer or contaminant plumes. Attenuation-difference tomography uses the difference in the trace amplitudes of tomographic data sets collected at different times to image the distribution of bulk-conductivity changes within the medium. The most common approach for computing the tomographic sensitivities uses ray theory, which is well understood and leads to efficient computations. However, ray theory requires the assumption that waves propagate at infinite frequency, and thus sensitivities are distributed along a line between the source and receiver. The infinite-frequency assumption in ray theory leads to a significant loss of resolution (both spatially and in terms of amplitude) of the recovered image. We use scattering theory to approximate the sensitivity of electromagnetic (EM) wave amplitude to changes in bulk conductivity within the medium.These sensitivities occupy the first Fresnel zone, account for the finite frequency nature of propagating EM waves, and are valid when velocity variations within the medium do not cause significant ray bending. We evaluate the scattering theory sensitivities by imaging a bromide tracer plume as it migrates through a coarse alluvial aquifer over two successive days. The scattering theory tomograms display a significant improvement in resolution over the ray-based counterparts, as shown by a direct comparison of the tomograms and also by a comparison of the vertical fluid conductivity distribution measured in a monitoring well, located within the tomographic plane. By improving resolution, the scattering theory sensitivities increase the utility of GPR attenuation-difference tomography for monitoring the movement of electrically anomalous plumes. In addition, the improved accuracy of information gathered through attenuation-difference tomography using scattering theory is a positive step toward future developments in using GPR data to help characterize the distribution of hydrogeologic properties.


2021 ◽  
Vol 11 (1) ◽  
Author(s):  
O. N. Senkov ◽  
D. B. Miracle

AbstractTwo classical criteria, by Pugh and Pettifor, have been widely used by metallurgists to predict whether a material will be brittle or ductile. A phenomenological correlation by Pugh between metal brittleness and its shear modulus to bulk modulus ratio was established more than 60 years ago. Nearly four decades later Pettifor conducted a quantum mechanical analysis of bond hybridization in a series of intermetallics and derived a separate ductility criterion based on the difference between two single-crystal elastic constants, C12–C44. In this paper, we discover the link between these two criteria and show that they are identical for materials with cubic crystal structures.


Author(s):  
Michele Di Cristo ◽  
Giacomo Milan

Abstract We investigate the reciprocity gap functional method, which has been developed in the inverse scattering theory, in the context of electrical impedance tomography. In particular, we aim to reconstruct an inclusion contained in a body, whose conductivity is different from the conductivity of the surrounding material. Numerical examples are given, showing the performance of our algorithm.


1983 ◽  
Vol 24 (6) ◽  
pp. 1502-1508 ◽  
Author(s):  
Helena M. A. de Castro ◽  
Walter F. Wreszinski

Sign in / Sign up

Export Citation Format

Share Document