scholarly journals Attribute-guided well-log interpolation applied to low-frequency impedance estimation

Geophysics ◽  
2008 ◽  
Vol 73 (6) ◽  
pp. R83-R95 ◽  
Author(s):  
Thomas Mejer Hansen ◽  
Klaus Mosegaard ◽  
Radmila Pedersen-Tatalovic ◽  
Anette Uldall ◽  
Nils Lange Jacobsen

Several approaches exist to use trends in 3D seismic data, in the form of seismic attributes, to interpolate sparsely sampled well-log measurements between well locations. Kriging and neural networks are two such approaches. We have applied a method that finds a relation between seismic attributes (such as two-way times, interval velocities, reflector roughness) and rock properties (in this case, acoustic impedance) from information at well locations. The relation is designed for optimum prediction of acoustic impedances away from well sites, and this is accomplished through a combination of cross validation and the Tikhonov-regularized least-squares method. The method is fast, works well even for highly underdetermined problems, and has general applicability. We apply it to two case studies in which we estimate 3D cubes of low-frequency impedance, which is essential for producing good porosity models. We show that the method is superior to traditional least squares: Numerous blind tests show that estimated low-frequency impedance away from well locations can be determined with an accuracy very close to estimations obtained at well locations.

Author(s):  
V. B. Olaseni ◽  
Y. S. Onifade ◽  
J. O. Airen ◽  
L. Adeoti

A statistically driven spectral method was carried out on 3D seismic data and well logs in ‘’VIC’’ Field within the Niger Delta with the aim of deriving reservoir properties and delineating stratigraphic features using edge detection attributes like coherence so as to have a better and clearer view of subsurface structure of a reservoir interval that possesses hydrocarbon using Spectral method. A suite of data consisting of seismic sections and composite logs comprising Gamma-ray, Resistivity, Spontaneous Potential, Sonic Time and Porosity logs (density and Neutron) were utilized to identify reservoir interval on log signature across wells 4 and 5 and the reservoir interval obtained was between 11,164 feet and 11,196 feet. Edge detection attribute like coherence was computed from the amplitude data in time domain and transformed to frequency domain using Fourier Transform tool in MATLAB. In order to display well log in time, well to seismic tie was carried out using check shot data which was used as time to depth relationship. The analysis of the spectral domain shows distinct bright spots that vary with measured depth due to variation in fluid and formation properties. The results led to an enhancement of seismic data interpretation in the field of study due to a spectral technique method that was applied to calculate the frequency slices. The results indicate that the spectral domain in coherence attributes revealed better geological features and the reservoir character such as faults and fractures. Frequency domain gives better geological maps as it is used to filter data, which means it is an enhancement of hidden features in time domain and gives a smoother variation of the features that has low frequency values. A reservoir with low frequency values is a sandy environment showing stratigraphy features. Hence, the reservoir is suspected to be a channel fill reservoir. This implies that Spectral domain (frequency) defines major geological areas of the ‘’VIC’’ field and gives much clearer image of the reservoir features within the field than in time domain.  


2021 ◽  
pp. 1-69
Author(s):  
Marwa Hussein ◽  
Robert R. Stewart ◽  
Deborah Sacrey ◽  
Jonny Wu ◽  
Rajas Athale

Net reservoir discrimination and rock type identification play vital roles in determining reservoir quality, distribution, and identification of stratigraphic baffles for optimizing drilling plans and economic petroleum recovery. Although it is challenging to discriminate small changes in reservoir properties or identify thin stratigraphic barriers below seismic resolution from conventional seismic amplitude data, we have found that seismic attributes aid in defining the reservoir architecture, properties, and stratigraphic baffles. However, analyzing numerous individual attributes is a time-consuming process and may have limitations for revealing small petrophysical changes within a reservoir. Using the Maui 3D seismic data acquired in offshore Taranaki Basin, New Zealand, we generate typical instantaneous and spectral decomposition seismic attributes that are sensitive to lithologic variations and changes in reservoir properties. Using the most common petrophysical and rock typing classification methods, the rock quality and heterogeneity of the C1 Sand reservoir are studied for four wells located within the 3D seismic volume. We find that integrating the geologic content of a combination of eight spectral instantaneous attribute volumes using an unsupervised machine-learning algorithm (self-organizing maps [SOMs]) results in a classification volume that can highlight reservoir distribution and identify stratigraphic baffles by correlating the SOM clusters with discrete net reservoir and flow-unit logs. We find that SOM classification of natural clusters of multiattribute samples in the attribute space is sensitive to subtle changes within the reservoir’s petrophysical properties. We find that SOM clusters appear to be more sensitive to porosity variations compared with lithologic changes within the reservoir. Thus, this method helps us to understand reservoir quality and heterogeneity in addition to illuminating thin reservoirs and stratigraphic baffles.


Geophysics ◽  
2018 ◽  
Vol 83 (4) ◽  
pp. V243-V252
Author(s):  
Wail A. Mousa

A stable explicit depth wavefield extrapolation is obtained using [Formula: see text] iterative reweighted least-squares (IRLS) frequency-space ([Formula: see text]-[Formula: see text]) finite-impulse response digital filters. The problem of designing such filters to obtain stable images of challenging seismic data is formulated as an [Formula: see text] IRLS minimization. Prestack depth imaging of the challenging Marmousi model data set was then performed using the explicit depth wavefield extrapolation with the proposed [Formula: see text] IRLS-based algorithm. Considering the extrapolation filter design accuracy, the [Formula: see text] IRLS minimization method resulted in an image with higher quality when compared with the weighted least-squares method. The method can, therefore, be used to design high-accuracy extrapolation filters.


2017 ◽  
Vol 25 (03) ◽  
pp. 1750022
Author(s):  
Xiuwei Yang ◽  
Peimin Zhu

Acoustic impedance (AI) from seismic inversion can indicate rock properties and can be used, when combined with rock physics, to predict reservoir parameters, such as porosity. Solutions to seismic inversion problem are almost nonunique due to the limited bandwidth of seismic data. Additional constraints from well log data and geology are needed to arrive at a reasonable solution. In this paper, sedimentary facies is used to reduce the uncertainty in inversion and rock physics modeling; the results not only agree with seismic data, but also conform to geology. A reservoir prediction method, which incorporates seismic data, well logs, rock physics and sedimentary facies, is proposed. AI was first derived by constrained sparse spike inversion (CSSI) using a sedimentary facies dependent low-frequency model, and then was transformed to reservoir parameters by sequential simulation, statistical rock physics and [Formula: see text]-model. Two numerical experiments using synthetic model and real data indicated that the sedimentary facies information may help to obtain a more reasonable prediction.


Geophysics ◽  
2002 ◽  
Vol 67 (6) ◽  
pp. 1877-1885 ◽  
Author(s):  
Xin‐Quan Ma

A new prestack inversion algorithm has been developed to simultaneously estimate acoustic and shear impedances from P‐wave reflection seismic data. The algorithm uses a global optimization procedure in the form of simulated annealing. The goal of optimization is to find a global minimum of the objective function, which includes the misfit between synthetic and observed prestack seismic data. During the iterative inversion process, the acoustic and shear impedance models are randomly perturbed, and the synthetic seismic data are calculated and compared with the observed seismic data. To increase stability, constraints have been built into the inversion algorithm, using the low‐frequency impedance and background Vs/Vp models. The inversion method has been successfully applied to synthetic and field data examples to produce acoustic and shear impedances comparable to log data of similar bandwidth. The estimated acoustic and shear impedances can be combined to derive other elastic parameters, which may be used for identifying of lithology and fluid content of reservoirs.


2014 ◽  
Vol 2 (2) ◽  
pp. 18-26
Author(s):  
Sonny Inichinbia ◽  
◽  
Peter O. Sule ◽  
Halidu Hamza ◽  
Aminu L. Ahmed

2020 ◽  
Vol 8 (2) ◽  
pp. 168
Author(s):  
Nyeneime O. Etuk ◽  
Mfoniso U. Aka ◽  
Okechukwu A. Agbasi ◽  
Johnson C. Ibuot

Seismic attributes were evaluated over Edi field, offshore Western Niger Delta, Nigeria, via 3D seismic data. Manual mappings of the horizons and faults on the in-lines and cross-lines of the seismic sections were done. Various attributes were calculated and out put on four horizons corresponding to the well markers at different formations within the well were identified. The four horizons identified, which includes: H1, H2, H3 and H4 were mapped and interpreted across the field. The operational agenda was thru picking given faults segments on the in–line of seismic volume. A total of five faults coded as F1, F2, F3, F4 and F5, F1 and F5 were the major fault and were observed as extending through the field. Structural and horizon mappings were used to generate time structure maps. The maps showed the various positions and orientations of the faults. Different attributes which include: root mean square amplitude, instantaneous phase, gradient magnitude and chaos were run on the 3D seismic data. The amplitude and incline magnitude maps indicate direct hydrocarbon on the horizon maps; this is very important in the drilling of wells because it shows areas where hydrocarbons are present in the subsurface. The seismic attributes revealed information, which was not readily apparent in the raw seismic data.   


Geophysics ◽  
2009 ◽  
Vol 74 (3) ◽  
pp. R15-R24 ◽  
Author(s):  
Taeyoung Ha ◽  
Wookeen Chung ◽  
Changsoo Shin

Waveform inversion faces difficulties when applied to real seismic data, including the existence of many kinds of noise. The [Formula: see text]-norm is more robust to noise with outliers than the least-squares method. Nevertheless, the least-squares method is preferred as an objective function in many algorithms because the gradient of the [Formula: see text]-norm has a singularity when the residual becomes zero. We propose a complex-valued Huber function for frequency-domain waveform inversion that combines the [Formula: see text]-norm (for small residuals) with the [Formula: see text]-norm (for large residuals). We also derive a discretized formula for the gradient of the Huber function. Through numerical tests on simple synthetic models and Marmousi data, we find the Huber function is more robust to outliers and coherent noise. We apply our waveform-inversion algorithm to field data taken from the continental shelf under the East Sea in Korea. In this setting, we obtain a velocity model whose synthetic shot profiles are similar to the real seismic data.


Sign in / Sign up

Export Citation Format

Share Document