Model‐based inversion of amplitude‐variations‐with‐offset data using a genetic algorithm

Geophysics ◽  
1995 ◽  
Vol 60 (4) ◽  
pp. 939-954 ◽  
Author(s):  
Subhashis Mallick

I cast the inversion of amplitude‐variation‐with‐offset (AVO) data into the framework of Bayesian statistics. Under such a framework, the model parameters and the physics of the forward problem are used to generate synthetic data. These synthetic data are then matched with the observed data to obtain an a‐posteriori probability density (PPD) function in the model space. The genetic algorithm (GA) uses a directed random search technique to estimate the shape of the PPD. Unlike the classical inversion methods, GA does not depend upon the choice of an initial model and is well suited for the AVO inversion. For the single‐layer AVO inversion where the amplitudes from a single reflection event are inverted, GA estimates the normal incidence reflection coefficient [Formula: see text] and the contrast of the Poisson’s ratio (Δσ) to reasonable accuracy, even when the signal‐to‐noise ratio is poor. Comparisons of single‐layer amplitude inversion using synthetic data demonstrate that GA inversion obtains more accurate results than does the least‐squares fit to the approximate reflection coefficients as is usually practiced in the industry. In the multilayer AVO waveform inversion, all or a part of the prestack data are inverted. Inversion of this type is nonunique for the estimation of the absolute values of velocities, Poisson’s ratios, and densities. However, by applying simplified approximations to the P‐wave reflection coefficient, I verify that [Formula: see text], the contrast in the acoustic impedance (ΔA), and the gradient in the reflection coefficient (G), can be estimated from such an inversion. From the GA estimated values of [Formula: see text], ΔA, and G, and from reliable estimates of velocity and Poisson’s ratio at the start time of the input data, an inverted model can be generated. I apply this procedure to marine data and demonstrate that the the synthetics computed from such an inverted model match the input data to reasonable accuracy. Comparison of the log data from a nearby well shows that the GA inversion obtains both the low‐ and the high‐frequency trends (within the bandwidth of seismic resolution) of the P‐wave acoustic impedance. In addition to P‐wave acoustic impedance, GA also obtains an estimate of the Poisson’s ratio, an extremely important parameter for the direct detection of hydrocarbons from seismic data.

Geophysics ◽  
1984 ◽  
Vol 49 (10) ◽  
pp. 1637-1648 ◽  
Author(s):  
W. J. Ostrander

The P-wave reflection coefficient at an interface separating two media is known to vary with angle of incidence. The manner in which it varies is strongly affected by the relative values of Poisson’s ratio in the two media. For moderate angles of incidence, the relative change in reflection coefficient is particularly significant when Poisson’s ratio differs greatly between the two media. Theory and laboratory measurements indicate that high‐porosity gas sands tend to exhibit abnormally low Poisson’s ratios. Embedding these low‐velocity gas sands into sediments having “normal” Poisson’s ratios should result in an increase in reflected P-wave energy with angle of incidence. This phenomenon has been observed on conventional seismic data recorded over known gas sands.


1974 ◽  
Vol 64 (2) ◽  
pp. 473-491
Author(s):  
Harold M. Mooney

abstract We consider a version of Lamb's Problem in which a vertical time-dependent point force acts on the surface of a uniform half-space. The resulting surface disturbance is computed as vertical and horizontal components of displacement, particle velocity, acceleration, and strain. The goal is to provide numerical solutions appropriate to a comparison with observed wave forms produced by impacts onto granite and onto soil. Solutions for step- and delta-function sources are not physically realistic but represent limiting cases. They show a clear P arrival (larger on horizontal than vertical components) and an obscure S arrival. The Rayleigh pulse includes a singularity at the theoretical arrival time. All of the energy buildup appears on the vertical components and all of the energy decay, on the horizontal components. The effects of Poisson's ratio upon vertical displacements for a step-function source are shown. For fixed shear velocity, an increase of Poisson's ratio produces a P pulse which is larger, faster, and more gradually emergent, an S pulse with more clear-cut beginning, and a much narrower Rayleigh pulse. For a source-time function given by cos2(πt/T), −T/2 ≦ T/2, a × 10 reduction in pulse width at fixed pulse height yields an increase in P and Rayleigh-wave amplitudes by factors of 1, 10, and 100 for displacement, velocity and strain, and acceleration, respectively. The observed wave forms appear somewhat oscillatory, with widths proportional to the source pulse width. The Rayleigh pulse appears as emergent positive on vertical components and as sharp negative on horizontal components. We show a theoretical seismic profile for granite, with source pulse width of 10 µsec and detectors at 10, 20, 30, 40, and 50 cm. Pulse amplitude decays as r−1 for P wave and r−12 for Rayleigh wave. Pulse width broadens slightly with distance but the wave form character remains essentially unchanged.


Nanoscale ◽  
2021 ◽  
Author(s):  
Lingling Bai ◽  
Yifan Gao ◽  
Junhao Peng ◽  
Xing Xie ◽  
Minru Wen ◽  
...  

Zero (or near-zero) Poisson's ratio (ZPR) materials have important applications in the field of precision instrument because one of its faces is stable and will not be affected by strain....


Geophysics ◽  
1994 ◽  
Vol 59 (9) ◽  
pp. 1352-1361 ◽  
Author(s):  
James W. Spencer ◽  
Michael E. Cates ◽  
Don D. Thompson

In this study, we investigate the elastic moduli of the empty grain framework (the “frame” moduli) in unconsolidated sands and consolidated sandstones. The work was done to improve the interpretation of seismic amplitude anomalies and amplitude variations with offset (AVO) associated with hydrocarbon reservoirs. We developed a laboratory apparatus to measure the frame Poisson’s ratio and Young’s modulus of unconsolidated sands at seismic frequencies (0.2 to 155 Hz) in samples approximately 11 cm long. We used ultrasonic pulse velocity measurements to measure the frame moduli of consolidated sandstones. We found that the correlation coefficient between the frame Poisson’s ratio [Formula: see text] and the mineral Poisson’s ratio [Formula: see text] is 0.84 in consolidated sandstones and only 0.28 in unconsolidated sands. The range of [Formula: see text] values in unconsolidated sands is 0.115 to 0.237 (mean = 0.187, standard deviation = 0.030), and [Formula: see text] cannot be estimated without core or log analyses. Frame moduli analyses of core samples can be used to calibrate the interpretation of seismic amplitude anomalies and AVO effects. For use in areas without core or log analyses, we developed an empirical relation that can be used to estimate [Formula: see text] in unconsolidated sands and sandstones from [Formula: see text] and the frame P‐wave modulus.


Author(s):  
Haohao Zhang ◽  
Jun Lu ◽  
Benchi Chen ◽  
Xuejun Ma ◽  
Zhidong Cai

Abstract The considerable depth and complicated structure of the Tahe Oilfield in the Tuofutai area of China make it very difficult to delineate its Ordovician carbonate fracture-cavity reservoir. The resolution of conventional ground seismic data is inadequate to satisfy current exploitation requirements. To further improve the understanding of the carbonate fracture-cavity reservoir of the Tahe Oilfield and to provide predictions of reservoir properties that are more accurate, a walkaround 3D-3C vertical seismic profiling (VSP) survey was conducted. First, we preprocessed raw VSP data and developed a method of joint PP- and PSV-wave prestack time migration. In contrast to ground seismic imaging profiles, VSP imaging profiles have a higher resolution and wider spectrum range that provide more detailed strata information. Then, using the joint PP- and PSV-wave prestack inversion method, we obtained the PP- and PSV-wave impedance and Poisson's ratio parameters of the Ordovician carbonate reservoir. Compared with the P-wave impedance of the ground seismic inversion, we found the VSP inversion results had higher accuracy, which enabled clearer identification of the internal characteristics of the carbonate reservoir. Finally, coupled with the Poisson's ratio attribute, we predicted the distribution of favorable reservoirs and interwell connectivity. The prediction results were verified using both logging and production data. The findings of this study demonstrate the applicability of the proposed technical method for the exploration of deep carbonate fracture-cavity reservoirs.


2014 ◽  
Vol 104 (8) ◽  
pp. 081902 ◽  
Author(s):  
Baolin Wang ◽  
Jiangtao Wu ◽  
Xiaokun Gu ◽  
Hanqing Yin ◽  
Yujie Wei ◽  
...  

2016 ◽  
Vol 4 (4) ◽  
pp. T613-T625 ◽  
Author(s):  
Qizhen Du ◽  
Bo Zhang ◽  
Xianjun Meng ◽  
Chengfeng Guo ◽  
Gang Chen ◽  
...  

Three-term amplitude-variation with offset (AVO) inversion generally suffers from instability when there is limited prior geologic or petrophysical constraints. Two-term AVO inversion shows higher instability compared with three-term AVO inversion. However, density, which is important in the fluid-type estimation, cannot be recovered from two-term AVO inversion. To reliably predict the P- and S-waves and density, we have developed a robust two-step joint PP- and PS-wave three-term AVO-inversion method. Our inversion workflow consists of two steps. The first step is to estimate the P- and S-wave reflectivities using Stewart’s joint two-term PP- and PS-AVO inversion. The second step is to treat the P-wave reflectivity obtained from the first step as the prior constraint to remove the P-wave velocity related-term from the three-term Aki-Richards PP-wave approximated reflection coefficient equation, and then the reduced PP-wave reflection coefficient equation is combined with the PS-wave reflection coefficient equation to estimate the S-wave and density reflectivities. We determined the effectiveness of our method by first applying it to synthetic models and then to field data. We also analyzed the condition number of the coefficient matrix to illustrate the stability of the proposed method. The estimated results using proposed method are superior to those obtained from three-term AVO inversion.


Geophysics ◽  
2004 ◽  
Vol 69 (1) ◽  
pp. 164-179 ◽  
Author(s):  
Shaoming Lu ◽  
George A. McMechan

The elastic properties of hydrated sediments are not well‐known, which leads to inaccuracy in the evaluation of the amount of gas hydrate worldwide. Elastic impedance inversion is useful in estimating the elastic properties of sediments containing gas hydrate, or free gas trapped beneath the gas hydrate, from angle‐dependent P‐wave reflections. We reprocess the multichannel U.S. Geological Survey seismic line BT‐1 from the Blake Ridge off the east coast of North America to obtain migrated common‐angle aperture data sets, which are then inverted for elastic impedance. Two new algorithms to estimate P‐impedance and S‐impedance from the elastic impedance are developed and evaluated using well‐log data from Ocean Drilling Program (ODP) Leg 164; these new algorithms are stable, even in the presence of modest noise in the data. The Vs/Vp ratio, Poisson's ratio, and Lamé parameter terms λρ and λ/μ are estimated from the P‐impedance and S‐impedance. The hydrated sediments have high elastic impedance, high P‐impedance, high S‐impedance, high λρ, slightly higher Vs/Vp ratio, slightly lower Poisson's ratio, and slightly lower λ/μ values compared to those of the surrounding unhydrated sediments. The sediments containing free gas have low elastic impedance, low P‐impedance, nonanomalous background S‐impedance, high Vs/Vp ratio, low Poisson's ratio, low λρ, and low λ/μ values. We conclude that some parameters such as Vs/Vp ratio, Poisson's ratio, and λ/μ, although they help identify the free‐gas charged layers, cannot differentiate between the hydrated sediments and nonhydrated sediments when gas hydrate concentration is low, and cannot differentiate between the hydrated sediments and free‐gas charged sediments when the gas hydrate concentration is high. Three distinct layers of gas hydrate are interpreted as being caused by gas hydrates with gas of different molecular weights, with correspondingly different stability zones in depth. Free gas appears to be present below the two deeper gas‐hydrate layers, but not below the shallowest one because the lack of a trapping structure. The gas hydrate has an average concentration of ∼3–5.5% by volume, and is highest (9%) at the base of the lower gas hydrate stability zone. The free‐gas concentration ranges from 1 to 8% by volume, and is most developed beneath the local topographic high of the ocean bottom.


Geophysics ◽  
1982 ◽  
Vol 47 (5) ◽  
pp. 819-824 ◽  
Author(s):  
Harsh K. Gupta ◽  
Ronald W. Ward ◽  
Tzeu‐Lie Lin

Analysis of P‐ and S‐waves from shallow microearthquakes in the vicinity of The Geysers geothermal area, California, recorded by a dense, telemetered seismic array operated by the U.S. Geological Survey (USGS) shows that these phases are easily recognized and traced on record sections to distances of 80 km. Regional average velocities for the upper crust are estimated to be [Formula: see text] and [Formula: see text] for P‐ and S‐waves, respectively. Poisson’s ratio is estimated at 23 locations using Wadati diagrams and is found to vary from 0.13 to 0.32. In general, the Poisson’s ratio is found to be lower at the locations close to the steam production zones at The Geysers and Clear Lake volcanic field to the northeast. The low Poisson ratio corresponds to a decrease in P‐wave velocity in areas of high heat flow. The decrease may be caused by fracturing of the rock and saturation with gas or steam.


Sign in / Sign up

Export Citation Format

Share Document