Estimation of free gas saturation from seismic reflection surveys by the genetic algorithm inversion of a P-wave attenuation model

Geophysics ◽  
2012 ◽  
Vol 77 (4) ◽  
pp. R175-R187 ◽  
Author(s):  
Eugene C. Morgan ◽  
Maarten Vanneste ◽  
Isabelle Lecomte ◽  
Laurie G. Baise ◽  
Oddvar Longva ◽  
...  
Keyword(s):  
Genetic Algorithm ◽  
Seismic Reflection ◽  
Wave Attenuation ◽  
In Situ Measurements ◽  
P Wave ◽  
Ratio Method ◽  
Attenuation Model ◽  
Free Gas ◽  
Gas Saturation

Many previously proposed methods of estimating free gas saturation from seismic survey data rely on calibration to invasively collected, in situ measurements. Typically, such in situ measurements are used to parameterize or calibrate rock-physics models, which can then be applied to seismic data to achieve saturation estimates. We tested a technique for achieving estimates of the spatial distribution of gas saturation solely from shipboard seismic surveys. We estimated the quality factor from seismic reflection surveys using the spectral ratio method, and then inverted a mesoscopic-scale P-wave attenuation model to find the parameters that matched the modeled attenuation to our estimates of observed attenuation within the range of seismic frequencies. By using a genetic algorithm for this inversion, we not only searched efficiently for a global solution to the nonlinear set of equations that compose the model, but also constrain the search to a relatively broad set of realistic parameter values. Thus, our estimates do not rely on in situ measurements of these parameters, but on distributions of their possible values, many of which may be referenced from literature. We first tested this method at Blake Ridge, offshore North and South Carolina, where an approximately 400-m-deep gas-saturated zone underlies a field of methane hydrates. The extensive field work and subsequent studies at this site make it ideal for validating our method. We also demonstrated the applicability of our method to shallower deposits by presenting results from Finneidfjord, Norway, where the inversion of the P-wave attenuation model recognizes very small gas saturations.

Mesoscopic P-wave attenuation model to estimation of free gas saturation

2013 ◽  
Author(s):  
Francisco Cabrera ◽  
Flor A. Vivas ◽  
German D. Camacho ◽  
Herling Gonzalez
Keyword(s):  
Wave Attenuation ◽  
P Wave ◽  
Attenuation Model ◽  
Free Gas ◽  
Gas Saturation

scholarly journals Anisotropic P-wave attenuation measured from a multi-azimuth surface seismic reflection survey

2009 ◽  
Vol 57 (5) ◽  
pp. 835-845 ◽  
Author(s):  
Roger A. Clark ◽  
Philip M. Benson ◽  
Andrew J. Carter ◽  
Carlos A. Guerrero Moreno
Keyword(s):  
Seismic Reflection ◽  
Wave Attenuation ◽  
P Wave

scholarly journals Local and regional P-wave spectral attenuation model for Iran

2020 ◽  
Vol 224 (1) ◽  
pp. 241-256
Author(s):  
Ehsan Moradian Bajestani ◽  
Anooshiravan Ansari ◽  
Ehsan Karkooti
Keyword(s):  
Wave Attenuation ◽  
Vertical Component ◽  
P Wave ◽  
P Waves ◽  
Attenuation Model ◽  
Hypocentral Distance ◽  
Anelastic Attenuation ◽  
Curve Versus ◽  
Path Coverage ◽  
Spectral Attenuation

SUMMARY A robust frequency-dependent local and regional P-wave attenuation model is estimated for continental paths in the Iranian Plateau. In order to calculate the average attenuation parameters, 46 337 vertical-component waveforms related to 9267 earthquakes, which are recorded at the Iranian Seismological Center (IRSC) stations, have been selected in the distance range 10–1000 km. The majority of the event's magnitudes are less than 4.5. This collection of records provides high spatial ray path coverage. Results indicate that the shape of attenuation P-wave curve versus distance is not uniform and has three distinct sections with hinges at 90 and 175 km. A trilinear model for attenuation of P-wave amplitude in the frequency range 1–10 Hz is proposed in this study. Fourier spectral amplitudes are found to decay as R−1.2 (where R is hypocentral distance), corresponding to geometric spreading within 90 km from the source. There is a section from 90 to 175 km, where the attenuation is described as R0.8, and the attenuation is described well beyond 175 km by R−1.3. Moreover, the average quality factor for Pg and Pn waves (QPg and QPn), related to anelastic attenuation is obtained as Qpg= (54.2 ± 2.6)f(1.0096±0.07) and Qpn= (306.8 ± 7.4)f (0.51±0.05). There is a good agreement between the results of the model and observations. Also, the attenuation model shows compatibility with the recent regional studies. From the results it turns out that the amplitude of P waves attenuates more rapidly in comparison with the global models in local distances.

P-wave attenuation measurements from laboratory resonance and sonic waveform data

Geophysics ◽  
1989 ◽  
Vol 54 (1) ◽  
pp. 76-81 ◽  
Author(s):  
D. Goldberg ◽  
B. Zinszner
Keyword(s):  
Wave Attenuation ◽  
Compressional Wave ◽  
P Wave ◽  
Laboratory Measurements ◽  
Frequency Range ◽  
Waveform Data ◽  
Sonic Log ◽  
In Situ Conditions ◽  
The Mean

We computed compressional‐wave velocity [Formula: see text] and attenuation [Formula: see text] from sonic log waveforms recorded in a cored, 30 m thick, dolostone reservoir; using cores from the same reservoir, laboratory measurements of [Formula: see text] and [Formula: see text] were also obtained. We used a resonant bar technique to measure extensional and shear‐wave velocities and attenuations in the laboratory, so that the same frequency range as used in sonic logging (5–25 kHz) was studied. Having the same frequency range avoids frequency‐dependent differences between the laboratory and in‐situ measurements. Compressional‐wave attenuations at 0 MPa confining pressure, calculated on 30 samples, gave average [Formula: see text] values of 17. Experimental and geometrical errors were estimated to be about 5 percent. Measurements at elevated effective pressures up to 30 MPa on selected dolostone samples in a homogeneous interval showed mean [Formula: see text] and [Formula: see text] to be approximately equal and consistently greater than 125. At effective stress of 20 MPa and at room temperature, the mean [Formula: see text] over the dolostone interval was 87, a minimum estimate for the approximate in‐situ conditions. We computed compressional‐wave attenuation from sonic log waveforms in the 12.5–25 kHz frequency band using the slope of the spectral ratio of waveforms recorded 0.914 m and 1.524 m from the source. Average [Formula: see text] over the interval was 13.5, and the mean error between this value and the 95 percent confidence interval of the slope was 15.9 percent. The laboratory measurements of [Formula: see text] under elevated pressure conditions were more than five times greater than the mean in‐situ values. This comparison shows that additional extrinsic losses in the log‐derived measurement of [Formula: see text], such as scattering from fine layers and vugs or mode conversion to shear energy dissipating radially from the borehole, dominate the apparent attenuation.

Physical modeling and analysis of P-wave attenuation anisotropy in transversely isotropic media

Geophysics ◽  
2007 ◽  
Vol 72 (1) ◽  
pp. D1-D7 ◽  
Author(s):  
Yaping Zhu ◽  
Ilya Tsvankin ◽  
Pawan Dewangan ◽  
Kasper van Wijk
Keyword(s):  
Attenuation Coefficient ◽  
Symmetry Axis ◽  
Wave Attenuation ◽  
Transversely Isotropic ◽  
P Wave ◽  
Velocity Variation ◽  
Ratio Method ◽  
Fracture Detection ◽  
Wide Range ◽  
Attenuation Anisotropy

Anisotropic attenuation can provide sensitive attributes for fracture detection and lithology discrimination. This paper analyzes measurements of the P-wave attenuation coefficient in a transversely isotropic sample made of phenolic material. Using the spectral-ratio method, we estimate the group (effective) attenuation coefficient of P-waves transmitted through the sample for a wide range of propagation angles (from [Formula: see text] to [Formula: see text]) with the symmetry axis. Correction for the difference between the group and phase angles and for the angular velocity variation help us to obtain the normalized phase attenuation coefficient [Formula: see text] governed by the Thomsen-style attenuation-anisotropy parameters [Formula: see text] and [Formula: see text]. Whereas the symmetry axis of the angle-dependent coefficient [Formula: see text] practically coincides with that of the velocity function, the magnitude of the attenuation anisotropy far exceeds that of the velocity anisotropy. The quality factor [Formula: see text] increases more than tenfold from the symmetry axis (slow direction) to the isotropy plane (fast direction). Inversion of the coefficient [Formula: see text] using the Christoffel equation yields large negative values of the parameters [Formula: see text] and [Formula: see text]. The robustness of our results critically depends on several factors, such as the availability of an accurate anisotropic velocity model and adequacy of the homogeneous concept of wave propagation, as well as the choice of the frequency band. The methodology discussed here can be extended to field measurements of anisotropic attenuation needed for AVO (amplitude-variation-with-offset) analysis, amplitude-preserving migration, and seismic fracture detection.

Identification of gas hydrate dissociation from wireline-log data in the Shenhu area, South China Sea

Geophysics ◽  
2012 ◽  
Vol 77 (3) ◽  
pp. B125-B134 ◽  
Author(s):  
Xiujuan Wang ◽  
Myung Lee ◽  
Shiguo Wu ◽  
Shengxiong Yang
Keyword(s):  
Wave Velocity ◽  
Gas Hydrate ◽  
Seismic Data ◽  
Pore Space ◽  
P Wave ◽  
Well Log ◽  
Hydrate Dissociation ◽  
Free Gas ◽  
P Wave Velocity

Wireline logs were acquired in eight wells during China’s first gas hydrate drilling expedition (GMGS-1) in April–June of 2007. Well logs obtained from site SH3 indicated gas hydrate was present in the depth range of 195–206 m below seafloor with a maximum pore-space gas hydrate saturation, calculated from pore water freshening, of about 26%. Assuming gas hydrate is uniformly distributed in the sediments, resistivity calculations using Archie’s equation yielded hydrate-saturation trends similar to those from chloride concentrations. However, the measured compressional (P-wave) velocities decreased sharply at the depth between 194 and 199 mbsf, dropping as low as [Formula: see text], indicating the presence of free gas in the pore space, possibly caused by the dissociation of gas hydrate during drilling. Because surface seismic data acquired prior to drilling were not influenced by the in situ gas hydrate dissociation, surface seismic data could be used to identify the cause of the low P-wave velocity observed in the well log. To determine whether the low well-log P-wave velocity was caused by in situ free gas or by gas hydrate dissociation, synthetic seismograms were generated using the measured well-log P-wave velocity along with velocities calculated assuming both gas hydrate and free gas in the pore space. Comparing the surface seismic data with various synthetic seismograms suggested that low P-wave velocities were likely caused by the dissociation of in situ gas hydrate during drilling.

scholarly journals In-situ measurements of atmospheric hydrofluorocarbons (HFCs) and perfluorocarbons (PFCs) at the Shangdianzi regional background station, China

2012 ◽  
Vol 12 (5) ◽  
pp. 11151-11173
Author(s):  
B. Yao ◽  
M. K. Vollmer ◽  
L. X. Zhou ◽  
S. Henne ◽  
S. Reimann ◽  
...  
Keyword(s):  
In Situ Measurements ◽  
Ratio Method ◽  
High Concentration ◽  
Urban Sector ◽  
Mixing Ratios ◽  
Regional Background ◽  
Hfc 134A ◽  
The Mean ◽  
Measurement Period

Abstract. In-situ measurements of atmospheric hydrofluorocarbons (HFCs) and perfluorocarbons (PFCs) have been conducted at the Shangdianzi (SDZ) Global Atmosphere Watch (GAW) regional background station, China, from May 2010 to May 2011. The time series for 5 HFCs and 4 PFCs periodically showed high concentration events while background conditions occurred for 36% (HFC-32) to 83% (PFC-218) of all measurements. The mean mixing ratios during background conditions for HFC-23, HFC-32, HFC-125, HFC-134a, HFC-152, CF4, PFC-116, PFC-218 and PFC-318 were 24.5, 5.86, 9.97, 66.0, 9.77, 79.1, 4.22, 0.56, 1.28 ppt (parts per trillion, 10−12, molar), respectively. The background mixing ratios for the compounds at SDZ are consistent with those obtained at mid to high latitude sites in the Northern Hemisphere, except for HFC-32 and PFC-318 for which background mixing ratios were not reported in recent years. All HFCs and PFCs show positive trends at rates of 0.7, 1.4, 1.6, 4.1, 1.1, 0.43, 0.05, 0.01, 0.04 ppt yr−1 for HFC-23, HFC-32, HFC-125, HFC-134a, HFC-152, CF4, PFC-116, PFC-218 and PFC-318, respectively. North-easterly winds were connected with small contributions to atmospheric HFCs and PFCs loadings, whereas south-westerly advection (urban sector) showed increased loadings. Chinese emissions were estimated by a tracer ratio method using CO as tracer with rather well known emissions. The emissions, as derived from our measurement period, were 4.4 ± 0.7, 6.9 ± 0.9, 2.5 ± 0.3, 9.0 ± 1.3, 2.2 ± 0.4, 2.1 ± 0.3, 0.24 ± 0.06, 0.07 ± 0.04, 0.45 ± 0.09 kt yr−1 for HFC-23, HFC-32, HFC-125, HFC-134a, HFC-152, CF4, PFC-116, PFC-218, and PFC-318, respectively. The lower HFC-23 emissions compared to earlier studies may be a result of the HFC-23 abatement measures taken as part of the Clean Development Mechanism (CDM) project that started in 2005.

scholarly journals Seismic wave attenuation in the uppermost glacier ice of Storglaciären, Sweden

2010 ◽  
Vol 56 (196) ◽  
pp. 249-256 ◽  
Author(s):  
Alessio Gusmeroli ◽  
Roger A. Clark ◽  
Tavi Murray ◽  
Adam D. Booth ◽  
Bernd Kulessa ◽  
...  
Keyword(s):  
Seismic Wave ◽  
Wave Attenuation ◽  
Seismic Refraction ◽  
Spectral Ratio ◽  
P Wave ◽  
Ratio Method ◽  
Polar Ice ◽  
Seismic Line ◽  
Seismic Wave Attenuation ◽  
Glacier Ice

AbstractWe conducted seismic refraction surveys in the upper ablation area of Storglaciären, a small valley glacier located in Swedish Lapland. We estimated seismic-wave attenuation using the spectral-ratio method on the energy travelling in the uppermost ice with an average temperature of approximately −1 °C. Attenuation values were derived between 100 and 300 Hz using the P-wave quality factor, QP, the inverse of the internal friction. By assuming constant attenuation along the seismic line we obtained mean QP = 6 ± 1. We also observed that QP varies from 8 ± 1 to 5 ± 1 from the near-offset to the far-offset region of the line, respectively. Since the wave propagates deeper at far offsets, this variation is interpreted by considering the temperature profile of the study area; far-offset arrivals sampled warmer and thus more-attenuative ice. Our estimates are considerably lower than those reported for field studies in polar ice (∼500–1700 at −28°C and 50–160 at −10°C) and, hence, are supportive of laboratory experiments that show attenuation increases with rising ice temperature. Our results provide new in situ estimates of QP for glacier ice and demonstrate a valuable method for future investigations in both alpine and polar ice.

Sign in / Sign up

Export Citation Format

Share Document