Using instantaneous frequency and colored inversion attributes to distinguish and determine the sandstones facies of the Late Ordovician Mamuniyat reservoir, R-field in Murzuq Basin, Libya

2016 ◽  
Vol 4 (4) ◽  
pp. T507-T519 ◽  
Author(s):  
Yousf Abushalah ◽  
Laura Serpa

The Mamuniyat petroleum reservoir in southwestern Libya is comprised of clean sandstones and intercalated shale and sand facies that are characterized by spatial porosity variations. Seismic reflection data from the field exhibit relatively low vertical seismic resolution, side lobes of reflection wavelets, reflection interference, and low acoustic impedance contrast between the reservoir and the units underneath the reservoir, which make mapping those facies a difficult task. In the absence of broadband seismic data, optimizing frequency bands of bandlimited data can be used to suppress pseudoreflectors resulting from side-lobe effects and help to separate the clean sandstone facies of the reservoir. We have optimized the data based on our investigation of seismic frequency bands and used instantaneous frequency analysis to reveal the reflection discontinuity that is mainly associated with the reservoir boundary of the sandstone facies of the clean Mamuniyat reservoir. We also preformed rock-physics diagnostic modeling and inverted the seismic data using spectral-based colored inversion into relative acoustic impedance. The inverted impedance matches the up-scaled impedance from the well data and the inversion of relative acoustic impedance confirms the conclusion that was drawn from the instantaneous frequency results. The interpretation of facies distributions based on the instantaneous frequency was supported by the inversion results and the rock-physics models.

Geophysics ◽  
2015 ◽  
Vol 80 (1) ◽  
pp. R31-R41 ◽  
Author(s):  
Andrea Zunino ◽  
Klaus Mosegaard ◽  
Katrine Lange ◽  
Yulia Melnikova ◽  
Thomas Mejer Hansen

Determination of a petroleum reservoir structure and rock bulk properties relies extensively on inference from reflection seismology. However, classic deterministic methods to invert seismic data for reservoir properties suffer from some limitations, among which are the difficulty of handling complex, possibly nonlinear forward models, and the lack of robust uncertainty estimations. To overcome these limitations, we studied a methodology to invert seismic reflection data in the framework of the probabilistic approach to inverse problems, using a Markov chain Monte Carlo (McMC) algorithm with the goal to directly infer the rock facies and porosity of a target reservoir zone. We thus combined a rock-physics model with seismic data in a single inversion algorithm. For large data sets, the McMC method may become computationally impractical, so we relied on multiple-point-based a priori information to quantify geologically plausible models. We tested this methodology on a synthetic reservoir model. The solution of the inverse problem was then represented by a collection of facies and porosity reservoir models, which were samples of the posterior distribution. The final product included probability maps of the reservoir properties in obtained by performing statistical analysis on the collection of solutions.


Geophysics ◽  
1993 ◽  
Vol 58 (3) ◽  
pp. 419-428 ◽  
Author(s):  
Arthur E. Barnes

Fourier power spectra are often usefully characterized by average measures. In reflection seismology, the important average measures are center frequency, spectral bandwidth, and dominant frequency. These quantities have definitions familiar from probability theory: center frequency is the spectral mean, spectral bandwidth is the standard deviation about that mean, and dominant frequency is the square root of the second moment, which serves as an estimate of the zero‐crossing frequency. These measures suggest counterparts defined with instantaneous power spectra in place of Fourier power spectra, so that they are instantaneous in time though they represent averages in frequency. Intuitively reasonable requirements yield specific forms for these instantaneous quantities that can be computed with familiar complex seismic trace attributes. Instantaneous center frequency is just instantaneous frequency. Instantaneous bandwidth is the absolute value of the derivative of the instantaneous amplitude divided by the instantaneous amplitude. Instantaneous dominant frequency is the square root of the sum of the squares of the instantaneous frequency and instantaneous bandwidth. Instantaneous bandwidth and dominant frequency find employment as additional complex seismic trace attributes in the detailed study of seismic data. Instantaneous bandwidth is observed to be nearly always less than instantaneous frequency; the points where it is larger may mark the onset of distinct wavelets. These attributes, together with instantaneous frequency, are perhaps, of greater use in revealing the time‐varying spectral properties of seismic data. They can help in the search for low frequency shadows or in the analysis of frequency change due to effects of data processing. Instantaneous bandwidth and dominant frequency complement instantaneous frequency and should find wide application in the analysis of seismic reflection data.


2021 ◽  
pp. 565-577
Author(s):  
Nowfal A. Nassir ◽  
Ahmed S. AL- Banna ◽  
Ghazi H. Al-Sharaa

The estimation of rock petrophysical parameters is an essential matter to characterize any reservoir. This research deals with the evaluation of  effective porosity (Pe), shale volume (Vsh) and water saturation (Sw) of reservoirs at Kumait and Dujalia fields, which were analyzed from well log and seismic data. The absolute acoustic impedance (AI) and relative acoustic impedance (RAI) were derived from a model which is based on the inversion of seismic 3-D post-stack data. NahrUmr formation’s sand reservoirs are identified by the RAI section of the study area. Nahr Umr sand-2 unit in Kumait field is the main reservoir; its delineation depends on the available well logs and AI sections information. The results of well logging interpretation showed a decrease of Sw and Vsh and an increase of effective porosity in the oil reservoir area, which coincides with the decrease of AI values. The existence of the water reservoir in Du-2 well revealed a convergence of the results of AI and effective porosity with those of Kumait wells , along with and some differential results of Sw and Vsh values that may be related to changes in lithology and fluid density.


Geophysics ◽  
2008 ◽  
Vol 73 (1) ◽  
pp. V1-V9 ◽  
Author(s):  
Chun-Feng Li ◽  
Christopher Liner

Although the passage of singularity information from acoustic impedance to seismic traces is now well understood, it remains unanswered how routine seismic processing, mode conversions, and multiple reflections can affect the singularity analysis of surface seismic data. We make theoretical investigations on the transition of singularity behaviors from acoustic impedances to surface seismic data. We also perform numerical, wavelet-based singularity analysis on an elastic synthetic data set that is processed through routine seismic processing steps (such as stacking and migration) and that contains mode conversions, multiple reflections, and other wave-equation effects. Theoretically, seismic traces can be approximated as proportional to a smoothed version of the [Formula: see text] derivative of acoustic impedance,where [Formula: see text] is the vanishing moment of the seismic wavelet. This theoretical approach forms the basis of linking singularity exponents (Hölder exponents) in acoustic impedance with those computable from seismic data. By using wavelet-based multiscale analysis with complex Morlet wavelets, we can estimate singularity strengths and localities in subsurface impedance directly from surface seismic data. Our results indicate that rich singularity information in acoustic impedance variations can be preserved by surface seismic data despite data-acquisition and processing activities. We also show that high-resolution detection of singularities from real surface seismic data can be achieved with a proper choice of the scale of the mother wavelet in the wavelet transform. Singularity detection from surface seismic data thus can play a key role in stratigraphic analysis and acoustic impedance inversion.


Geophysics ◽  
2020 ◽  
Vol 85 (5) ◽  
pp. M73-M83 ◽  
Author(s):  
Leonardo Azevedo ◽  
Dario Grana ◽  
Leandro de Figueiredo

Accurate subsurface modeling and characterization require the prediction of facies and rock properties within the reservoir model. This is commonly achieved by inverting geophysical data, such as seismic reflection data, using a two-step approach either in the discrete or the continuous domain. We have adopted an iterative simultaneous method, namely, stochastic perturbation optimization, to invert seismic reflection data jointly for facies and rock properties. Facies first are simulated according to a Markov chain model, and then rock properties are generated with stochastic sequential simulation and cosimulation conditioned to each facies. Elastic and seismic data are computed by applying a rock-physics model to the realizations of petrophysical properties and a seismic convolutional model. The similarity between observed and synthetic seismic data is used to update the solution by perturbing facies and rock properties until convergence. Coupling the discrete and continuous domains ensures a consistent perturbation of the reservoir models throughout the iterations. We have evaluated the method in a 1D synthetic example for the estimation of facies and porosity from zero-offset seismic data assuming a linear rock-physics model to demonstrate the validity of the method. Then, we apply the method to a real 3D data set from the North Sea for the joint estimation of facies and petrophysical properties from prestack seismic data. The results show spatially consistent rock and fluid inverted models in which the predicted facies reproduce the vertical ordering as observed in the well data.


Geophysics ◽  
2021 ◽  
pp. 1-60
Author(s):  
Francesco Turco ◽  
Leonardo Azevedo ◽  
Dario Grana ◽  
Gareth J. Crutchley ◽  
Andrew R. Gorman

Quantitative characterization of gas hydrate systems on continental margins from seismic data is challenging, especially in regions where no well logs are available. However, probabilistical seismic inversion provides an effective means for constraining the physical properties of subsurface strata in such settings and analyzing the variability related to the results. We apply a workflow for the characterization of two deep-water gas hydrate reservoirs east of New Zealand, where high concentrations of gas hydrate have been inferred previously. We estimate porosity and gas hydrate saturation in the reservoirs from multi-channel seismic data through a two-step procedure based on geostatistical seismic and Bayesian petrophysical inversion built on a rock physics model for gas hydrate-bearing marine sediments. We found that the two reservoirs together host between 2.45 × 105 m3 and 1.72 × 106 m3 of gas hydrate, with the best estimate at 9.68 × 105 m3. This estimate provides a first-order assessment for further gas hydrate evaluations in the region. The two-step statistically based seismic inversion method is an effective approach for characterizing gas hydrate systems from long-offset seismic reflection data.


2021 ◽  
pp. 1-29
Author(s):  
Papia Nandi ◽  
Patrick Fulton ◽  
James Dale

As rising ocean temperatures can destabilize gas hydrate, identifying and characterizing large shallow hydrate bodies is increasingly important in order to understand their hazard potential. In the southwestern Gulf of Mexico, reanalysis of 3D seismic reflection data reveals evidence for the presence of six potentially large gas hydrate bodies located at shallow depths below the seafloor. We originally interpreted these bodies as salt, as they share common visual characteristics on seismic data with shallow allochthonous salt bodies, including high-impedance boundaries and homogenous interiors with very little acoustic reflectivity. However, when seismic images are constructed using acoustic velocities associated with salt, the resulting images were of poor quality containing excessive moveout in common reflection point (CRP) offset image gathers. Further investigation reveals that using lower-valued acoustic velocities results in higher quality images with little or no moveout. We believe that these lower acoustic values are representative of gas hydrate and not of salt. Directly underneath these bodies lies a zone of poor reflectivity, which is both typical and expected under hydrate. Observations of gas in a nearby well, other indicators of hydrate in the vicinity, and regional geologic context, all support the interpretation that these large bodies are composed of hydrate. The total equivalent volume of gas within these bodies is estimated to potentially be as large as 1.5 gigatons or 10.5 TCF, considering uncertainty for estimates of porosity and saturation, comparable to the entire proven natural gas reserves of Trinidad and Tobago in 2019.


2018 ◽  
Vol 123 (12) ◽  
pp. 10,810-10,830
Author(s):  
Michael Dentith ◽  
Huaiyu Yuan ◽  
Ruth Elaine Murdie ◽  
Perla Pina-Varas ◽  
Simon P. Johnson ◽  
...  

Author(s):  
A. Ogbamikhumi ◽  
T. Tralagba ◽  
E. E. Osagiede

Field ‘K’ is a mature field in the coastal swamp onshore Niger delta, which has been producing since 1960. As a huge producing field with some potential for further sustainable production, field monitoring is therefore important in the identification of areas of unproduced hydrocarbon. This can be achieved by comparing production data with the corresponding changes in acoustic impedance observed in the maps generated from base survey (initial 3D seismic) and monitor seismic survey (4D seismic) across the field. This will enable the 4D seismic data set to be used for mapping reservoir details such as advancing water front and un-swept zones. The availability of good quality onshore time-lapse seismic data for Field ‘K’ acquired in 1987 and 2002 provided the opportunity to evaluate the effect of changes in reservoir fluid saturations on time-lapse amplitudes. Rock physics modelling and fluid substitution studies on well logs were carried out, and acoustic impedance change in the reservoir was estimated to be in the range of 0.25% to about 8%. Changes in reservoir fluid saturations were confirmed with time-lapse amplitudes within the crest area of the reservoir structure where reservoir porosity is 0.25%. In this paper, we demonstrated the use of repeat Seismic to delineate swept zones and areas hit with water override in a producing onshore reservoir.


2021 ◽  
Author(s):  
Piotr Krzywiec ◽  
Łukasz Słonka ◽  
Quang Nguyen ◽  
Michał Malinowski ◽  
Mateusz Kufrasa ◽  
...  

<p>In 2016, approximately 850 km of high-resolution multichannel seismic reflection data of the BALTEC survey have been acquired offshore Poland within the transition zone between the East European Craton and the Paleozoic Platform. Data processing, focused on removal of multiples, strongly overprinting geological information at shallower intervals, included SRME, TAU-P domain deconvolution, high resolution parabolic Radon demultiple and SWDM (Shallow Water De-Multiple). Entire dataset was Kirchhoff pre-stack time migrated. Additionally, legacy shallow high-resolution multichannel seismic reflection data acquired in this zone in 1997 was also used. All this data provided new information on various aspects of the Phanerozoic evolution of this area, including Late Cretaceous to Cenozoic tectonics and sedimentation. This phase of geological evolution could be until now hardly resolved by analysis of industry seismic data as, due to limited shallow seismic imaging and very strong overprint of multiples, essentially no information could have been retrieved from this data for first 200-300 m. Western part of the BALTEC dataset is located above the offshore segment of the Mid-Polish Swell (MPS) – large anticlinorium formed due to inversion of the axial part of the Polish Basin. BALTEC seismic data proved that Late Cretaceous inversion of the Koszalin – Chojnice fault zone located along the NE border of the MPS was thick-skinned in nature and was associated with substantial syn-inversion sedimentation. Subtle thickness variations and progressive unconformities imaged by BALTEC seismic data within the Upper Cretaceous succession in vicinity of the Kamień-Adler and the Trzebiatów fault zones located within the MPS documented complex interplay of Late Cretaceous basin inversion, erosion and re-deposition. Precambrian basement of the Eastern, cratonic part of the study area is overlain by Cambro-Silurian sedimentary cover. It is dissected by a system of steep, mostly reverse faults rooted in most cases in the deep basement. This fault system has been regarded so far as having been formed mostly in Paleozoic times, due to the Caledonian orogeny. As a consequence, Upper Cretaceous succession, locally present in this area, has been vaguely defined as a post-tectonic cover, locally onlapping uplifted Paleozoic blocks. New seismic data, because of its reliable imaging of the shallowest substratum, confirmed that at least some of these deeply-rooted faults were active as a reverse faults in latest Cretaceous – earliest Paleogene. Consequently, it can be unequivocally proved that large offshore blocks of Silurian and older rocks presently located directly beneath the Cenozoic veneer must have been at least partly covered by the Upper Cretaceous succession; then, they were uplifted during the widespread inversion that affected most of Europe. Ensuing regional erosion might have at least partly provided sediments that formed Upper Cretaceous progradational wedges recently imaged within the onshore Baltic Basin by high-end PolandSPAN regional seismic data. New seismic data imaged also Paleogene and younger post-inversion cover. All these results prove that Late Cretaceous tectonics substantially affected large areas located much farther towards the East than previously assumed.</p><p>This study was funded by the Polish National Science Centre (NCN) grant no UMO-2017/27/B/ST10/02316.</p>


Sign in / Sign up

Export Citation Format

Share Document