Fracture detection and fluid identification based on anisotropic Gassmann equation and linear-slip model

Geophysics ◽  
2019 ◽  
Vol 84 (1) ◽  
pp. R85-R98 ◽  
Author(s):  
Xinpeng Pan ◽  
Guangzhi Zhang
Keyword(s):  
Seismic Data ◽  
Order Approximation ◽  
Low Frequency ◽  
Rock Physics ◽  
Transversely Isotropic ◽  
Slip Model ◽  
Fracture Detection ◽  
Fluid Identification ◽  
Pp Wave ◽  
Linear Slip Model

Detection of fracture and fluid properties from subsurface azimuthal seismic data improves our abilities to characterize the saturated porous reservoirs with aligned fractures. Motivated by the fracture detection and fluid identification in a fractured porous medium, we have developed a feasible approach to perform a rock physics model-based amplitude variation with offset and azimuth (AVOAz) inversion for the fracture and fluid parameters in a horizontal transversely isotropic (HTI) medium using the PP-wave angle gathers along different azimuths. Based on the linear-slip model, we first use anisotropic Gassmann’s equation to derive the expressions of saturated stiffness components and their perturbations of first-order approximation in terms of elastic properties of an isotropic porous background and fracture compliances induced by a single set of rotationally invariant fractures. We then derive a linearized PP-wave reflection coefficient in terms of fluid modulus, dry-rock matrix term, porosity, density, and fracture compliances or quasi-compliances for an interface separating two weakly HTI media based on the Born scattering theory. Finally, we solve the AVOAz inverse problems iteratively constrained by the Cauchy-sparse regularization and the low-frequency regularization in a Bayesian framework. The results demonstrate that the fluid modulus and fracture quasi-compliances are reasonably estimated in the case of synthetic and real seismic data containing moderate noise in a gas-filled fractured porous reservoir.

scholarly journals Elastic-Impedance-Based Fluid/Porosity Term and Fracture Weaknesses Inversion in Transversely Isotropic Media with a Tilted Axis of Symmetry

Geofluids ◽  
2020 ◽  
Vol 2020 ◽  
pp. 1-17
Author(s):  
Xinpeng Pan ◽  
Lin Li ◽  
Guangzhi Zhang ◽  
Yian Cui
Keyword(s):  
Reflection Coefficient ◽  
Seismic Data ◽  
Wave Reflection ◽  
Transversely Isotropic ◽  
Unknown Parameters ◽  
Fluid Identification ◽  
Pp Wave ◽  
Axis Of Symmetry ◽  
Dip Angle ◽  
Elastic Impedance

The rock containing a set of tilted fractures is equivalent to a transversely isotropic (TTI) medium with a tilted axis of symmetry. To implement fluid identification and tilted fracture detection, we propose an inversion approach of utilizing seismic data to simultaneously estimate parameters that are sensitive to fluids and tilted fractures. We first derive a PP-wave reflection coefficient and elastic impedance (EI) in terms of the dip angle, fluid/porosity term, shear modulus, density, and fracture weaknesses, and we present numerical examples to demonstrate how the PP-wave reflection coefficient and EI vary with the dip angle. Based on the information of dip angle of fractures provided by geologic and well data, we propose a two-step inversion approach of utilizing azimuthal seismic data to estimate unknown parameters involving the fluid/porosity term and fracture weaknesses: (1) the constrained sparse spike inversion (CSSI) for azimuthally anisotropic EI data and (2) the estimation of unknown parameters with the low-frequency constrained regularization term. Synthetic and real data demonstrate that fluid and fracture parameters are reasonably estimated, which may help fluid identification and fracture characterization.

Reservoir Prediction Under Control of Sedimentary Facies

2017 ◽  
Vol 25 (03) ◽  
pp. 1750022
Author(s):  
Xiuwei Yang ◽  
Peimin Zhu
Keyword(s):  
Seismic Data ◽  
Low Frequency ◽  
Prediction Method ◽  
Rock Physics ◽  
Real Data ◽  
Sedimentary Facies ◽  
Seismic Inversion ◽  
Rock Properties ◽  
Reservoir Prediction ◽  
Frequency Model

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.

Building tilted transversely isotropic depth models using localized anisotropic tomography with well information

Geophysics ◽  
2010 ◽  
Vol 75 (4) ◽  
pp. D27-D36 ◽  
Author(s):  
Andrey Bakulin ◽  
Marta Woodward ◽  
Dave Nichols ◽  
Konstantin Osypov ◽  
Olga Zdraveva
Keyword(s):  
Seismic Data ◽  
Model Building ◽  
Symmetry Axis ◽  
Transverse Isotropy ◽  
Synthetic Data ◽  
Rock Physics ◽  
Transversely Isotropic ◽  
Data Set ◽  
Anisotropic Models ◽  
Well Data

Tilted transverse isotropy (TTI) is increasingly recognized as a more geologically plausible description of anisotropy in sedimentary formations than vertical transverse isotropy (VTI). Although model-building approaches for VTI media are well understood, similar approaches for TTI media are in their infancy, even when the symmetry-axis direction is assumed known. We describe a tomographic approach that builds localized anisotropic models by jointly inverting surface-seismic and well data. We present a synthetic data example of anisotropic tomography applied to a layered TTI model with a symmetry-axis tilt of 45 degrees. We demonstrate three scenarios for constraining the solution. In the first scenario, velocity along the symmetry axis is known and tomography inverts for Thomsen’s [Formula: see text] and [Formula: see text] parame-ters. In the second scenario, tomography inverts for [Formula: see text], [Formula: see text], and velocity, using surface-seismic data and vertical check-shot traveltimes. In contrast to the VTI case, both these inversions are nonunique. To combat nonuniqueness, in the third scenario, we supplement check-shot and seismic data with the [Formula: see text] profile from an offset well. This allows recovery of the correct profiles for velocity along the symmetry axis and [Formula: see text]. We conclude that TTI is more ambiguous than VTI for model building. Additional well data or rock-physics assumptions may be required to constrain the tomography and arrive at geologically plausible TTI models. Furthermore, we demonstrate that VTI models with atypical Thomsen parameters can also fit the same joint seismic and check-shot data set. In this case, although imaging with VTI models can focus the TTI data and match vertical event depths, it leads to substantial lateral mispositioning of the reflections.

Seismic inversion of a carbonate buildup: A case study

2017 ◽  
Vol 5 (4) ◽  
pp. T641-T652 ◽  
Author(s):  
Mark Sams ◽  
Paul Begg ◽  
Timur Manapov
Keyword(s):  
Seismic Data ◽  
Probability Distributions ◽  
Low Frequency ◽  
Rock Physics ◽  
Seismic Inversion ◽  
Accurate Estimation ◽  
Amplitude Analysis ◽  
Simultaneous Inversion ◽  
Band Limited ◽  
Using Data

The information within seismic data is band limited and angle limited. Together with the particular physics and geology of carbonate rocks, this imposes limitations on how accurately we can predict the presence of hydrocarbons in carbonates, map the top carbonate, and characterize the porosity distribution through seismic amplitude analysis. Using data for a carbonate reef from the Nam Con Son Basin, Vietnam, the expectations based on rock-physics analysis are that the presence of gas can be predicted only when the porosity at the top of the carbonate is extremely high ([Formula: see text]), but that a fluid contact is unlikely to be observed in the background of significant porosity variations. Mapping the top of the carbonate (except when the top carbonate porosities are low) or a fluid contact requires accurate estimates of changes in [Formula: see text]. The seismic data do not independently support such an accurate estimation of sharp changes in [Formula: see text]. The standard approach of introducing low-frequency models and applying rock-physics constraints during a simultaneous inversion does not resolve the problems: The results are heavily biased by the well control and the initial interpretation of the top carbonate and fluid contact. A facies-based inversion in which the elastic properties are restricted to values consistent with the facies predicted to be present removes the well bias, but it does not completely obviate the need for a reasonably accurate initial interpretation in terms of prior facies probability distributions. Prestack inversion improves the quality of the facies predictions compared with a poststack inversion.

Azimuthally anisotropic elastic impedance inversion for fluid indicator driven by rock physics

Geophysics ◽  
2017 ◽  
Vol 82 (6) ◽  
pp. C211-C227 ◽  
Author(s):  
Xinpeng Pan ◽  
Guangzhi Zhang ◽  
Xingyao Yin
Keyword(s):  
Fractured Reservoirs ◽  
Low Frequency ◽  
Rock Physics ◽  
Real Data ◽  
Inversion Method ◽  
Anisotropic Rock ◽  
Weak Anisotropy ◽  
Fluid Identification ◽  
Anisotropic Elastic ◽  
Elastic Impedance

The normal-to-tangential fracture compliance ratio is usually used as a fracture fluid indicator (FFI) for fluid identification in fractured reservoirs. With a new parameterization for fracture weaknesses, we have defined a new FFI based on azimuthally anisotropic elastic impedance (EI) inversion and fractured anisotropic rock-physics models. First, we derived a new azimuthally anisotropic EI equation with a similar expression for the isotropic and anisotropic EI parts to remove the exponential correction of EI that is attributable to weak anisotropy. Then, we built a fractured anisotropic rock-physics model used for the estimation of well-log parameters for the normal and tangential fracture weaknesses, which built the initial background low-frequency trend of fracture weaknesses. Finally, based on the azimuthally anisotropic EI inversion method with the Cauchy-sparse and low-frequency information regularization, we estimated an FFI applied to fluid identification in fractured reservoirs. Tests on the synthetic and real data demonstrate that the anisotropic parameters related to fracture weaknesses can be estimated reasonably and stably and that our method appears to provide an alternative available for fluid identification in fractured reservoirs.

Seismic reservoir characterization of Utica-Point Pleasant shale with efforts at fracability evaluation — Part 2: A case study

2018 ◽  
Vol 6 (2) ◽  
pp. T325-T336 ◽  
Author(s):  
Ritesh Kumar Sharma ◽  
Satinder Chopra ◽  
James Keay ◽  
Hossein Nemati ◽  
Larry Lines
Keyword(s):  
Seismic Data ◽  
Reservoir Characterization ◽  
Low Frequency ◽  
Rock Physics ◽  
Well Log ◽  
X Ray Diffraction ◽  
Impedance Inversion ◽  
Seismic Reservoir ◽  
Petrophysical Modeling

The Utica Formation in eastern Ohio possesses all the prerequisites for being a successful unconventional play. Attempts at seismic reservoir characterization of the Utica Formation have been discussed in part 1, in which, after providing the geologic background of the area of study, the preconditioning of prestack seismic data, well-log correlation, and building of robust low-frequency models for prestack simultaneous impedance inversion were explained. All these efforts were aimed at identification of sweet spots in the Utica Formation in terms of organic richness as well as brittleness. We elaborate on some aspects of that exercise, such as the challenges we faced in the determination of the total organic carbon (TOC) volume and computation of brittleness indices based on mineralogical and geomechanical considerations. The prediction of TOC in the Utica play using a methodology, in which limited seismic as well as well-log data are available, is demonstrated first. Thereafter, knowing the nonexistence of the universally accepted indicator of brittleness, mechanical along with mineralogical attempts to extract the brittleness information for the Utica play are discussed. Although an attempt is made to determine brittleness from mechanical rock-physics parameters (Young’s modulus and Poisson’s ratio) derived from seismic data, the available X-ray diffraction data and regional petrophysical modeling make it possible to determine the brittleness index based on mineralogical data and thereafter be derived from seismic data.

Modeling processed seismic data to improve seismic facies prediction

2015 ◽  
Vol 3 (4) ◽  
pp. SAC91-SAC98 ◽  
Author(s):  
Adrian Pelham
Keyword(s):  
Seismic Data ◽  
Confusion Matrix ◽  
Low Frequency ◽  
Rock Physics ◽  
Seismic Facies ◽  
Reservoir Modeling ◽  
Elastic Parameters ◽  
Frequency Model ◽  
Lithology Prediction ◽  
Well Data

Interpreters need to screen and select the most geologically robust inversion products from increasingly larger data volumes, particularly in the absence of significant well control. Seismic processing and inversion routines are devised to provide reliable elastic parameters ([Formula: see text] and [Formula: see text]) from which the interpreter can predict the fluid and lithology properties. Seismic data modeling, for example, the Shuey approximations and the convolution inversion models, greatly assist in the parameterization of the processing flows within acceptable uncertainty limits and in establishing a measure of the reliability of the processing. Joint impedance facies inversion (Ji-Fi®) is a new inversion methodology that jointly inverts for acoustic impedance and seismic facies. Seismic facies are separately defined in elastic space ([Formula: see text] and [Formula: see text]), and a dedicated low-frequency model per facies is used. Because Ji-Fi does not need well data from within the area to define the facies or depth trends, wells from outside the area or theoretical constraints may be used. More accurate analyses of the reliability of the inversion products are a key advance because the results of the Ji-Fi lithology prediction may then be quantitatively and independently assessed at well locations. We used a novel visual representation of a confusion matrix to quantitatively assess the sensitivity and uncertainty in the results when compared with facies predicted from the depth trends and well-elastic parameters and the well-log lithologies observed. Thus, using simple models and the Ji-Fi inversion technique, we had an improved, quantified understanding of our data, the processes that had been applied, the parameterization, and the inversion results. Rock physics could further transform the elastic properties to more reservoir-focused parameters: volume of shale and porosity, volumes of facies, reservoir property uncertainties — all information required for interpretation and reservoir modeling.

Incorporation of geology with rock physics enables subsalt poststack inversion: A case study in the Gulf of Mexico

Geophysics ◽  
2011 ◽  
Vol 76 (5) ◽  
pp. WB53-WB65 ◽  
Author(s):  
Huyen Bui ◽  
Jennifer Graham ◽  
Shantanu Kumar Singh ◽  
Fred Snyder ◽  
Martiris Smith
Keyword(s):  
Seismic Data ◽  
Acoustic Impedance ◽  
Low Frequency ◽  
Rock Physics ◽  
Seismic Inversion ◽  
Inversion Method ◽  
Reservoir Properties ◽  
Blind Test ◽  
Frequency Model ◽  
Layer Sequence

One of the main goals of seismic inversion is to obtain high-resolution relative and absolute impedance for reservoir properties prediction. We aim to study whether the results from seismic inversion of subsalt data are sufficiently robust for reliable reservoir characterization. Approximately [Formula: see text] of poststack, wide-azimuth, anisotropic (vertical transverse isotropic) wave-equation migration seismic data from 50 Outer Continental Shelf blocks in the Green Canyon area of the Gulf of Mexico were inverted in this study. A total of four subsalt wells and four subsalt seismic interpreted horizons were used in the inversion process, and one of the wells was used for a blind test. Our poststack inversion method used an iterative discrete spike inversion method, based on the combination of space-adaptive wavelet processing to invert for relative acoustic impedance. Next, the dips were estimated from seismic data and converted to a horizon-like layer sequence field that was used as one of the inputs into the low-frequency model. The background model was generated by incorporating the well velocities, seismic velocity, seismic interpreted horizons, and the previously derived layer sequence field in the low-frequency model. Then, the relative acoustic impedance volume was scaled by adding the low-frequency model to match the calculated acoustic impedance logs from the wells for absolute acoustic impedance. Finally, the geological information and rock physics data were incorporated into the reservoir properties assessment for sand/shale prediction in two main target reservoirs in the Miocene and Wilcox formations. Overall, the poststack inversion results and the sand/shale prediction showed good ties at the well locations. This was clearly demonstrated in the blind test well. Hence, incorporating rock physics and geology enables poststack inversion in subsalt areas.

APPLICATION OF THE EMPIRICAL MODE DECOMPOSITION METHOD TO GROUND-ROLL NOISE ATTENUATION IN SEISMIC DATA

2013 ◽  
Vol 31 (4) ◽  
pp. 619 ◽  
Author(s):  
Luiz Eduardo Soares Ferreira ◽  
Milton José Porsani ◽  
Michelângelo G. Da Silva ◽  
Giovani Lopes Vasconcelos
Keyword(s):  
Empirical Mode Decomposition ◽  
Seismic Data ◽  
Low Frequency ◽  
High Amplitude ◽  
Seismic Line ◽  
Seismic Processing ◽  
Mode Decomposition ◽  
Ground Roll ◽  
Fortran 90 ◽  
Emd Method

ABSTRACT. Seismic processing aims to provide an adequate image of the subsurface geology. During seismic processing, the filtering of signals considered noise is of utmost importance. Among these signals is the surface rolling noise, better known as ground-roll. Ground-roll occurs mainly in land seismic data, masking reflections, and this roll has the following main features: high amplitude, low frequency and low speed. The attenuation of this noise is generally performed through so-called conventional methods using 1-D or 2-D frequency filters in the fk domain. This study uses the empirical mode decomposition (EMD) method for ground-roll attenuation. The EMD method was implemented in the programming language FORTRAN 90 and applied in the time and frequency domains. The application of this method to the processing of land seismic line 204-RL-247 in Tacutu Basin resulted in stacked seismic sections that were of similar or sometimes better quality compared with those obtained using the fk and high-pass filtering methods.Keywords: seismic processing, empirical mode decomposition, seismic data filtering, ground-roll. RESUMO. O processamento sísmico tem como principal objetivo fornecer uma imagem adequada da geologia da subsuperfície. Nas etapas do processamento sísmico a filtragem de sinais considerados como ruídos é de fundamental importância. Dentre esses ruídos encontramos o ruído de rolamento superficial, mais conhecido como ground-roll . O ground-roll ocorre principalmente em dados sísmicos terrestres, mascarando as reflexões e possui como principais características: alta amplitude, baixa frequência e baixa velocidade. A atenuação desse ruído é geralmente realizada através de métodos de filtragem ditos convencionais, que utilizam filtros de frequência 1D ou filtro 2D no domínio fk. Este trabalho utiliza o método de Decomposição em Modos Empíricos (DME) para a atenuação do ground-roll. O método DME foi implementado em linguagem de programação FORTRAN 90, e foi aplicado no domínio do tempo e da frequência. Sua aplicação no processamento da linha sísmica terrestre 204-RL-247 da Bacia do Tacutu gerou como resultados, seções sísmicas empilhadas de qualidade semelhante e por vezes melhor, quando comparadas as obtidas com os métodos de filtragem fk e passa-alta.Palavras-chave: processamento sísmico, decomposição em modos empíricos, filtragem dados sísmicos, atenuação do ground-roll.

Sign in / Sign up

Export Citation Format

Share Document