scholarly journals Inversion and interpretation of seismic-derived rock properties in the Duvernay play

2018 ◽  
Vol 6 (2) ◽  
pp. SE1-SE14 ◽  
Author(s):  
Ronald M. Weir ◽  
David W. Eaton ◽  
Larry R. Lines ◽  
Donald C. Lawton ◽  
Eneanwan Ekpo

We have developed an interpretive seismic workflow that incorporates multicomponent seismic inversion, guided by structural mapping, for characterizing low-permeability unconventional reservoirs. The workflow includes the determination of a calibrated time-depth relationship, generation of seismic-derived structural maps, poststack inversion, amplitude-variation-with-offset analysis, and PP-PS joint inversion. The subsequent interpretation procedure combines structural and inversion results with seismic-derived lithologic parameters, such as the Young’s modulus, Poisson’s ratio, and brittleness index. We applied this workflow to a 3D multicomponent seismic data set from the Duvernay play in the Kaybob area in Alberta, Canada. Subtle faults are discernible using isochron maps, horizontal time slices, and seismic stratal slices. Fault-detection software is also used to aid in the delineation of structural discontinuities. We found that seismic-derived attributes, coupled with structural mapping, can be used to map reservoir facies and thus to highlight zones that are most favorable for hydraulic-fracture stimulation. By imaging structural discontinuities and preexisting zones of weakness, seismic mapping also contributes to an improved framework for understanding the induced-seismicity risk.

2019 ◽  
Vol 38 (2) ◽  
pp. 151-160 ◽  
Author(s):  
Ronald Weir ◽  
Don Lawton ◽  
Laurence Lines ◽  
Thomas Eyre ◽  
David Eaton

Simultaneous prestack inversion of multicomponent 3D seismic data integrated with structural interpretation can provide an effective workflow to maximize value for unconventional plays. We outline an integrated workflow for characterizing the Duvernay play in western Canada, an emerging world-class low-permeability unconventional resource fairway. This workflow includes the determination of a time-depth relationship using synthetic seismograms, generation of seismic-derived time- and depth-converted structural maps, and calculation of inversion-based parameters of density and P- and S-wave velocity. The model-based procedure includes poststack (acoustic) inversion, amplitude variation with offset prestack inversion, and joint PP-PS inversion. With these rock properties determined, calculations are made to determine Young's modulus, Poisson's ratio, and brittleness. Faults are mapped based on time slices, isochrons, and correlatable vertical displacements of stratigraphic marker reflections. Significant strike-slip movements are identified by lateral displacement on interpreted geologic features, such as channels and reef edges. Seismic-derived attributes, combined with structural mapping, highlight zones that are conducive to hydraulic fracturing as well as areas unfavorable for development. Mapping of structural discontinuities provides a framework for understanding zones of preexisting weakness and induced-seismicity hazards.


Geophysics ◽  
1993 ◽  
Vol 58 (3) ◽  
pp. 408-418 ◽  
Author(s):  
L. R. Jannaud ◽  
P. M. Adler ◽  
C. G. Jacquin

A method developed for the determination of the characteristic lengths of an heterogeneous medium from the spectral analysis of codas is based on an extension of Aki’s theory to anisotropic elastic media. An equivalent Gaussian model is obtained and seems to be in good agreement with the two experimental data sets that illustrate the method. The first set was obtained in a laboratory experiment with an isotropic marble sample. This sample is characterized by a submillimetric length scale that can be directly observed on a thin section. The spectral analysis of codas and their inversion yields an equivalent correlation length that is in good agreement with the observed one. The second data set is obtained in a crosshole experiment at the usual scale of a seismic survey. The codas are recorded, analysed, and inverted. The analysis yields a vertical characteristic length for the studied subsurface that compares well with the characteristic length measured by seismic and stratigraphic logs.


2017 ◽  
Vol 5 (3) ◽  
pp. SL1-SL8 ◽  
Author(s):  
Ehsan Zabihi Naeini ◽  
Russell Exley

Quantitative interpretation (QI) is an important part of successful exploration, appraisal, and development activities. Seismic amplitude variation with offset (AVO) provides the primary signal for the vast majority of QI studies allowing the determination of elastic properties from which facies can be determined. Unfortunately, many established AVO-based seismic inversion algorithms are hindered by not fully accounting for inherent subsurface facies variations and also by requiring the addition of a preconceived low-frequency model to supplement the limited bandwidth of the input seismic. We apply a novel joint impedance and facies inversion applied to a North Sea prospect using broadband seismic data. The focus was to demonstrate the significant advantages of inverting for each facies individually and iteratively determine an optimized low-frequency model from facies-derived depth trends. The results generated several scenarios for potential facies distributions thereby providing guidance to future appraisal and development decisions.


Geophysics ◽  
2020 ◽  
Vol 86 (1) ◽  
pp. C1-C18
Author(s):  
Xinpeng Pan ◽  
Lin Li ◽  
Shunxin Zhou ◽  
Guangzhi Zhang ◽  
Jianxin Liu

The characterization of fracture-induced tilted transverse isotropy (TTI) seems to be more suitable to actual scenarios of geophysical exploration for fractured reservoirs. Fracture weaknesses enable us to describe fracture-induced anisotropy. With the incident and reflected PP-wave in TTI media, we have adopted a robust method of azimuthal amplitude variation with offset (AVO) parameterization and inversion for fracture weaknesses in a fracture-induced reservoir with TTI symmetry. Combining the linear-slip model with the Bond transformation, we have derived the stiffness matrix of a dipping-fracture-induced TTI medium characterized by normal and tangential fracture weaknesses and a tilt angle. Integrating the first-order perturbations in the stiffness matrix of a TTI medium and scattering theory, we adopt a method of azimuthal AVO parameterization for PP-wave reflection coefficient for the case of a weak-contrast interface separating two homogeneous weakly anisotropic TTI layers. We then adopt an iterative inversion method by using the partially incidence-angle-stacked seismic data with different azimuths to estimate the fracture weaknesses of a TTI medium when the tilt angle is estimated based on the image well logs prior to the seismic inversion. Synthetic examples confirm that the fracture weaknesses of a TTI medium are stably estimated from the azimuthal seismic reflected amplitudes for the case of moderate noise. A field data example demonstrates that geologically reasonable results of fracture weaknesses can be determined when the tilt angle is treated as the prior information. We determine that the azimuthal AVO inversion approach provides an available tool for fracture prediction in a dipping-fracture-induced TTI reservoir.


2014 ◽  
Vol 2 (2) ◽  
pp. SC37-SC46
Author(s):  
Jessa-lyn Lee

Using a combined amplitude-variation-with-offset (AVO) and inversion workflow, a 2D seismic data set was used to predict the lateral lithology changes in a sand bar deposit cut by shale-filled channels. Drilling into a shale channel drastically affects the success of a well, so understanding the spatial extent of such a feature is considered important to the economics of the development program. Based on results from extensive pre- and poststack modeling, a combination of AVO and inversion attributes provided the best chance at lithology prediction. Using this method, a channel map was created for a small study area. It is important to identify the risks and uncertainties intrinsic in the processes that are being applied as well as the effect of the overlying geology. This example in particular showed how vital it is to understand the geology in a specific area, know the technical limits of the data being used, and adapt workflows accordingly.


2020 ◽  
Vol 17 (6) ◽  
pp. 929-939
Author(s):  
Daiane R Rosa ◽  
Juliana M C Santos ◽  
Rafael M Souza ◽  
Dario Grana ◽  
Denis J Schiozer ◽  
...  

Abstract Time-lapse (4D) seismic inversion aims to predict changes in elastic rock properties, such as acoustic impedance, from measured seismic amplitude variations due to hydrocarbon production. Possible approaches for 4D seismic inversion include two classes of method: sequential independent 3D inversions and joint inversion of 4D seismic differences. We compare the standard deterministic methods, such as coloured and model-based inversions, and the probabilistic inversion techniques based on a Bayesian approach. The goal is to compare the sequential independent 3D seismic inversions and the joint 4D inversion using the same type of algorithm (Bayesian method) and to benchmark the results to commonly applied algorithms in time-lapse studies. The model property of interest is the ratio of the acoustic impedances, estimated for the monitor, and base surveys at each location in the model. We apply the methods to a synthetic dataset generated based on the Namorado field (offshore southeast Brazil). Using this controlled dataset, we can evaluate properly the results as the true solution is known. The results show that the Bayesian 4D joint inversion, based on the amplitude difference between seismic surveys, provides more accurate results than sequential independent 3D inversion approaches, and these results are consistent with deterministic methods. The Bayesian 4D joint inversion is relatively easy to apply and provides a confidence interval of the predictions.


2020 ◽  
Vol 4 ◽  
pp. 27-29
Author(s):  
A.V. Novoyavchev ◽  
◽  
A.A. Kleimenov ◽  
M.Yu. Tokarev ◽  
K.M. Myatchin ◽  
...  

Animals ◽  
2020 ◽  
Vol 11 (1) ◽  
pp. 50
Author(s):  
Jennifer Salau ◽  
Jan Henning Haas ◽  
Wolfgang Junge ◽  
Georg Thaller

Machine learning methods have become increasingly important in animal science, and the success of an automated application using machine learning often depends on the right choice of method for the respective problem and data set. The recognition of objects in 3D data is still a widely studied topic and especially challenging when it comes to the partition of objects into predefined segments. In this study, two machine learning approaches were utilized for the recognition of body parts of dairy cows from 3D point clouds, i.e., sets of data points in space. The low cost off-the-shelf depth sensor Microsoft Kinect V1 has been used in various studies related to dairy cows. The 3D data were gathered from a multi-Kinect recording unit which was designed to record Holstein Friesian cows from both sides in free walking from three different camera positions. For the determination of the body parts head, rump, back, legs and udder, five properties of the pixels in the depth maps (row index, column index, depth value, variance, mean curvature) were used as features in the training data set. For each camera positions, a k nearest neighbour classifier and a neural network were trained and compared afterwards. Both methods showed small Hamming losses (between 0.007 and 0.027 for k nearest neighbour (kNN) classification and between 0.045 and 0.079 for neural networks) and could be considered successful regarding the classification of pixel to body parts. However, the kNN classifier was superior, reaching overall accuracies 0.888 to 0.976 varying with the camera position. Precision and recall values associated with individual body parts ranged from 0.84 to 1 and from 0.83 to 1, respectively. Once trained, kNN classification is at runtime prone to higher costs in terms of computational time and memory compared to the neural networks. The cost vs. accuracy ratio for each methodology needs to be taken into account in the decision of which method should be implemented in the application.


2006 ◽  
Vol 06 (04) ◽  
pp. 373-384
Author(s):  
ERIC BERTHONNAUD ◽  
JOANNÈS DIMNET

Joint centers are obtained from data treatment of a set of markers placed on the skin of moving limb segments. Finite helical axis (FHA) parameters are calculated between time step increments. Artifacts associated with nonrigid body movements of markers entail ill-determination of FHA parameters. Mean centers of rotation may be calculated over the whole movement, when human articulations are likened to spherical joints. They are obtained using numerical technique, defining point with minimal amplitude, during joint movement. A new technique is presented. Hip, knee, and ankle mean centers of rotation are calculated. Their locations depend on the application of two constraints. The joint center must be located next to the estimated geometric joint center. The geometric joint center may migrate inside a cube of possible location. This cube of error is located with respect to the marker coordinate systems of the two limb segments adjacent to the joint. Its position depends on the joint and the patient height, and is obtained from a stereoradiographic study with specimen. The mean position of joint center and corresponding dispersion are obtained through a minimization procedure. The location of mean joint center is compared with the position of FHA calculated between different sequential steps: time sequential step, and rotation sequential step where a minimal rotation amplitude is imposed between two joint positions. Sticks are drawn connecting adjacent mean centers. The animation of stick diagrams allows clinical users to estimate the displacements of long bones (femur and tibia) from the whole data set.


Sign in / Sign up

Export Citation Format

Share Document