Generation of a Probabilistic Facial Model on the Basis of Lithology Logs, Well Logs, Seismic Data and Existing Methods of Machine Learning and Classification

2021 ◽  
Author(s):  
Henadz Zaitsau ◽  
Valeri Shumilyak ◽  
Alexander Konyushenko

Abstract The main topic of an article is machine learning and classification (neural net) use for prognostic lithological model creation. Moreover, research preceding stages such as attribute analysis, seismic inversion, seismogeological modeling and briefly the results of lithological and petrophysical investigations are described in this art

2021 ◽  
Author(s):  
Pavlo Kuzmenko ◽  
Rustem Valiakhmetov ◽  
Francesco Gerecitano ◽  
Viktor Maliar ◽  
Grigori Kashuba ◽  
...  

Abstract The seismic data have historically been utilized to perform structural interpretation of the geological subsurface. Modern approaches of Quantitative Interpretation are intended to extract geologically valuable information from the seismic data. This work demonstrates how rock physics enables optimal prediction of reservoir properties from seismic derived attributes. Using a seismic-driven approach with incorporated prior geological knowledge into a probabilistic subsurface model allowed capturing uncertainty and quantifying the risk for targeting new wells in the unexplored areas. Elastic properties estimated from the acquired seismic data are influenced by the depositional environment, fluid content, and local geological trends. By applying the rock physics model, we were able to predict the elastic properties of a potential lithology away from the well control points in the subsurface whether or not it has been penetrated. Seismic amplitude variation with incident angle (AVO) and azimuth (AVAZ) jointly with rock-derived petrophysical interpretations were used for stochastical modeling to capture the reservoir distribution over the deep Visean formation. The seismic inversion was calibrated by available well log data and by traditional structural interpretation. Seismic elastic inversion results in a deep Lower Carboniferous target in the central part of the DDB are described. The fluid has minimal effect on the density and Vp. Well logs with cross-dipole acoustics are used together with wide-azimuth seismic data, processed with amplitude control. It is determined that seismic anisotropy increases in carbonate deposits. The result covers a set of lithoclasses and related probabilities: clay minerals, tight sandstones, porous sandstones, and carbonates. We analyzed the influence of maximum angles determination for elastic inversion that varied from 32.5 to 38.5 degrees. The greatest influence of the far angles selection is on the density. AI does not change significantly. Probably the 38,5 degrees provides a superior response above the carbonates. It does not seem to damage the overall AVA behavior, which result in a good density outcome, as higher angles of incidence are included. It gives a better tie to the wells for the high density layers over the interval of interest. Sand probability cube must always considered in the interpretation of the lithological classification that in many cases may be misleading (i.e. when sand and shale probabilities are very close to each other, because of small changes in elastic parameters). The authors provide an integrated holistic approach for quantitative interpretation, subsurface modeling, uncertainty evaluation, and characterization of reservoir distribution using pre-existing well logs and recently acquired seismic data. This paper underpins the previous efforts and encourages the work yet to be fulfilled on this subject. We will describe how quantitative interpretation was used for describing the reservoir, highlight values and uncertainties, and point a way forward for further improvement of the process for effective subsurface modeling.


2021 ◽  
pp. 1-57
Author(s):  
Chen Liang ◽  
John Castagna ◽  
Marcelo Benabentos

Sparse reflectivity inversion of processed reflection seismic data is intended to produce reflection coefficients that represent boundaries between geological layers. However, the objective function for sparse inversion is usually dominated by large reflection coefficients which may result in unstable inversion for weak events, especially those interfering with strong reflections. We propose that any seismogram can be decomposed according to the characteristics of the inverted reflection coefficients which can be sorted and subset by magnitude, sign, and sequence, and new seismic traces can be created from only reflection coefficients that pass sorting criteria. We call this process reflectivity decomposition. For example, original inverted reflection coefficients can be decomposed by magnitude, large ones removed, the remaining reflection coefficients reconvolved with the wavelet, and this residual reinverted, thereby stabilizing inversions for the remaining weak events. As compared with inverting an original seismic trace, subtle impedance variations occurring in the vicinity of nearby strong reflections can be better revealed and characterized when only the events caused by small reflection coefficients are passed and reinverted. When we apply reflectivity decomposition to a 3D seismic dataset in the Midland Basin, seismic inversion for weak events is stabilized such that previously obscured porous intervals in the original inversion, can be detected and mapped, with good correlation to actual well logs.


Geophysics ◽  
2021 ◽  
pp. 1-67
Author(s):  
Luanxiao Zhao ◽  
Caifeng Zou ◽  
Yuanyuan Chen ◽  
Wenlong Shen ◽  
Yirong Wang ◽  
...  

Seismic prediction of fluid and lithofacies distributions is of great interest to reservoir characterization, geological model building, and flow unit delineation. Inferring fluids and lithofacies from seismic data under the framework of machine learning is commonly subject to issues of limited features, imbalanced data sets, and spatial constraints. As a consequence, an XGBoost based workflow, which takes feature engineering, data balancing, and spatial constraints into account, is proposed to predict the fluid and lithofacies distribution by integrating well-log and seismic data. The constructed feature set based on simple mathematical operations and domain knowledge outperforms the benchmark group consisting of conventional elastic attributes of P-impedance and Vp/Vs ratio. A radial basis function characterizing the weights of training samples according to the distances from the available wells to the target region is developed to impose spatial constraints on the model training process, significantly improving the prediction accuracy and reliability of gas sandstone. The strategy combining the synthetic minority oversampling technique (SMOTE) and spatial constraints further increases the F1 score of gas sandstone and also benefits the overall prediction performance of all the facies. The application of the combined strategy on prestack seismic inversion results generates a more geologically reasonable spatial distribution of fluids, thus verifying the robustness and effectiveness of the proposed workflow.


2014 ◽  
Vol 54 (2) ◽  
pp. 532
Author(s):  
Yazeed Altowairqi ◽  
Reza Rezaee ◽  
Milovan Urosevic

Unconventional resources such as shale gas have been an extremely important exploration and production target. To understand the seismic responses of the shale gas plays, the use of rock physical relationship is important, which is constrained with geology and formation-evaluation analysis. Since organic-rich shale seismic properties remains poorly understood, seismic inversion can be used to identify the organic-rich shale from barren shale. This approach helps identify and map spatial distributions and of the organic rich shales. This study shows the acoustic impedance (AI), which is the product of compressional velocity and density, decreases nonlinearly with increasing total organic carbon (TOC) content. TOC is obtained using Roc-Eval pyrolysis for more than 120 core shale samples for the Perth Basin. By converting the AI data to TOC precent on the seismic data, we therefore can map lateral distribution, thickness, and variation in TOC profile. This extended abstract presents a case study of the northern Perth Basin 3D seismic with application of different approaches of seismic inversion and multi-attribute analysis with the rock physical relationships.


2020 ◽  
Vol 8 (4) ◽  
pp. T917-T925
Author(s):  
Bo Zhang ◽  
Yahua Yang ◽  
Yong Pan ◽  
Hao Wu ◽  
Danping Cao

The accuracy of seismic inversion is affected by the seismic wavelet and time-depth relationship generated by the process of the seismic well tie. The seismic well tie is implemented by comparing the synthetic seismogram computed from well logs and the poststack seismogram at or nearby the borehole location. However, precise waveform matching between the synthetic seismogram and the seismic trace does not guarantee an accurate tie between the elastic properties contained represented by the seismic data and well logs. We have performed the seismic well tie using the impedance log and the impedance inverted from poststack seismic data. We use an improved dynamic time warping to align the impedance log and impedance inverted from seismic data. Our workflow is similar to the current procedure of the seismic well tie except that the matching is implemented between the impedance log and the inverted impedance. The current seismic well-tie converges if there is no visible changes for the wavelets and time-depth relationship in the previous and current tying loops. Similarly, our seismic well tie converges if there are no visible changes for the wavelets, inverted impedance, and time-depth relationship in the previous and current tying loops. The real data example illustrates that more accurate inverted impedance is obtained by using the new wavelet and time-depth relationship.


Geophysics ◽  
2020 ◽  
Vol 85 (4) ◽  
pp. WA185-WA200
Author(s):  
Yuqing Chen ◽  
Gerard T. Schuster

We present a wave-equation inversion method that inverts skeletonized seismic data for the subsurface velocity model. The skeletonized representation of the seismic traces consists of the low-rank latent-space variables predicted by a well-trained autoencoder neural network. The input to the autoencoder consists of seismic traces, and the implicit function theorem is used to determine the Fréchet derivative, i.e., the perturbation of the skeletonized data with respect to the velocity perturbation. The gradient is computed by migrating the shifted observed traces weighted by the skeletonized data residual, and the final velocity model is the one that best predicts the observed latent-space parameters. We denote this as inversion by Newtonian machine learning (NML) because it inverts for the model parameters by combining the forward and backward modeling of Newtonian wave propagation with the dimensional reduction capability of machine learning. Empirical results suggest that inversion by NML can sometimes mitigate the cycle-skipping problem of conventional full-waveform inversion (FWI). Numerical tests with synthetic and field data demonstrate the success of NML inversion in recovering a low-wavenumber approximation to the subsurface velocity model. The advantage of this method over other skeletonized data methods is that no manual picking of important features is required because the skeletal data are automatically selected by the autoencoder. The disadvantage is that the inverted velocity model has less resolution compared with the FWI result, but it can serve as a good initial model for FWI. Our most significant contribution is that we provide a general framework for using wave-equation inversion to invert skeletal data generated by any type of neural network. In other words, we have combined the deterministic modeling of Newtonian physics and the pattern matching capabilities of machine learning to invert seismic data by NML.


2013 ◽  
Vol 347-350 ◽  
pp. 2447-2451
Author(s):  
Guo He Li ◽  
Xiang Yue ◽  
Wei Jiang Wu ◽  
Jiang Hui Zhao

In order to set up universal and non-linear map of variables, a full binary tree is constructed as mathematical model. Leaf nodes of the full binary tree are linear combination of input variables, and used as inputs of next nodes. On the basis of weighting two inputs by selector for inner node, the inputs are again linearly combined and used as output for next node. The inputs and outputs of all the inner nodes are constructed in turn as the same, and the output of root node is the output of mathematical model, implementing segment-linear approximation. With the means of machine learning of particle swarm optimization for data from some areas, all the coefficients of mathematical model are achieved for the special. The mathematical model is applied to seismic inversion to interpret stratum by seismic data, approving it very practical.


Geophysics ◽  
2021 ◽  
pp. 1-50
Author(s):  
Ahmad Mustafa ◽  
Motaz Alfarraj ◽  
Ghassan AlRegib

Seismic inversion plays a very useful role in detailed stratigraphic interpretation of migrated seismic volumes by enabling the estimation of reservoir properties over the complete volume. Traditional and machine learning-based seismic inversion workflows are limited to inverting each seismic trace independently of other traces to estimate impedance profiles, leading to lateral discontinuities in the presence of noise and large geological variations in the seismic data. In addition, machine learning-based approaches suffer the problem of overfitting if there is only a small number of wells on which the model is trained. We propose a two-pronged strategy to overcome these problems. We present a temporal convolutional network that models seismic traces temporally. We further inject spatial context for each trace into its estimations of the impedance profile. To counter the problem of limited labeled data, we also present a joint learning scheme whereby multiple datasets are simultaneously used for training, sharing beneficial information among each other. This results in the improvement in generalization performance on all datasets. We present a case study of acoustic impedance inversion using the open-source SEAM and Marmousi 2 datasets. Our evaluations show that our proposed approach is able to estimate impedance in the presence of noisy seismic data and a limited number of well logs with greater robustness and spatial consistency. We compare and contrast our approach to other learning-based seismic inversion methodologies in the literature. On SEAM, we are able to obtain an average MSE of 0.0476, the lowest among all other methodologies.


2021 ◽  
Author(s):  
Khalid Obaid ◽  
Muhammad Aamir ◽  
Tarek Yehia Nafie ◽  
Omar Aly ◽  
Widad Krissat ◽  
...  

Abstract Rock physics/seismic inversion is a powerful tool that deliver information about intra-wells rocks elastic attributes and reservoir properties such as porosity, saturation and rock lithology classification. In principle, inversion is like an engine that should be fueled by proper input quality of both seismic and well data. As for the well data, sonic and density logs measure the rock properties a few inches from the borehole. Reliability of sonic transit-time and bulk density logs can be affected by large and rapid variation in the diameter and shape of the borehole cross-section, as well as the process of drilling fluid invasion. The basic assumption for acoustic well logs editing and conditioning is to use other recorded logs (not affected by bad-hole conditions) in a Multivariate-Regression Algorithm. In addition, Fluid Substitution was implemented to correct for the mud invasion that affects the acoustic and elastic properties based on the PVT data for fluid properties computation. The logs were then quality checked by multiple cross-plotting comparisons to the standard Rock-physics trends templates. As for seismic data, there are several factors affecting the quality of surface seismic data including the presence of residual noise and multiples contamination that caused improper amplitude balancing. Optimizing the seismic data processing for the inversion studies require reviewing and conditioning the seismic gathers and pre-stack volumes, guided by a deterministic seismic-to-well tie analysis after every major stage of the processing sequence. The applied processes are mainly consisting of Curvelet domain noise attenuation to attenuate residual noise. This was followed by high resolution Radon anti-multiple to attenuate residual surface multiples and Extended interbed multiple prediction to attenuate interbed multiples. In addition, Offset dependent amplitude and spectral balancing were applied to maintain the seismic amplitudes fidelity. This paper will illustrate a case from Abu Dhabi where data conditioning results improved the Hydrocarbon saturated carbonates vs brine saturated carbonate and the lithology predictions, leading to optimizing field development plans and drilling operations.


Sign in / Sign up

Export Citation Format

Share Document