Elastic full-waveform inversion (EFWI) is a high-resolution technique for inverting velocity structures, such as P-wave and S-wave velocities, enabling the indirect estimation of subsurface material properties like porosity and saturation through petrophysical relationships. However, conventional EFWI faces challenges such as low signal-to-noise ratio (SNR) in near-surface land seismic data and heterogeneous Poisson’s ratios in complex environments, which hinder convergence and introduce crosstalk noise during the inversion process. In this work, we propose an integrated approach that combines a neural network-based seismic inversion with a petrophysical inversion framework to enhance the estimation of soil saturation and porosity properties. We introduce an advanced strategy called elastic neural network full-waveform inversion (ENFWI), which employs convolutional neural networks (CNNs) and automatic differentiation (AD) to invert the <inline-formula> <tex-math notation="LaTeX">$V_{p}$ </tex-math></inline-formula> and <inline-formula> <tex-math notation="LaTeX">$V_{s}$ </tex-math></inline-formula> models. This method utilizes two independent CNN architectures to generate <inline-formula> <tex-math notation="LaTeX">$V_{p}$ </tex-math></inline-formula> and <inline-formula> <tex-math notation="LaTeX">$V_{s}$ </tex-math></inline-formula> gradients, effectively suppressing crosstalk noise and providing robustness against low SNR conditions. Subsequently, we develop a petrophysical iterative inversion framework using Hertz-Mindlin (HM) contact theory to update porosity and saturation properties based on the inverted velocity models and resistivity constraint term. The AD algorithm integrates multiple equations to account for the <inline-formula> <tex-math notation="LaTeX">$V_{p}/V_{s}$ </tex-math></inline-formula> ratio and resistivity as constraints. Two typical synthetic models show that the proposed framework inverts reliable saturation, porosity, and velocity models. Furthermore, field data illustrate the framework’s capability to accurately identify a water-bearing interlayer and the Wadi basement interface. This work offers a reliable method for imaging structural distributions and estimating properties in near-surface critical zone (CZ) applications.

Accurate reservoir characterization is necessary to effectively monitor, manage, and increase production. A seismic inversion methodology using a genetic algorithm (GA) and particle swarm optimization (PSO) technique is proposed in this study to characterize the reservoir both qualitatively and quantitatively. It is usually difficult and expensive to map deeper reservoirs in exploratory operations when using conventional approaches for reservoir characterization hence inversion based on advanced technique (GA and PSO) is proposed in this study. The main goal is to use GA and PSO to significantly lower the fitness (error) function between real seismic data and modeled synthetic data, which will allow us to estimate subsurface properties and accurately characterize the reservoir. Both techniques estimate subsurface properties in a comparable manner. Consequently, a qualitative and quantitative comparison is conducted between these two algorithms. Using two synthetic data and one real data from the Blackfoot field in Canada, the study examined subsurface acoustic impedance and porosity in the inter-well zone. Porosity and acoustic impedance are layer features, but seismic data is an interface property, hence these characteristics provide more useful and applicable reservoir information. The inverted results aid in the understanding of seismic data by providing incredibly high-resolution images of the subsurface. Both the GA and the PSO algorithms deliver outstanding results for both simulated and real data. The inverted section accurately delineated a high porosity zone ( > 20 % ) that supported the high seismic amplitude anomaly by having a low acoustic impedance (6000-8500 m/s ∗ g/cc). This unusual zone is categorized as a reservoir (sand channel) and is located in the 1040-1065 ms time range. In this inversion process, after 400 iterations, the fitness error falls from 1 to 0.88 using GA optimization, compared to 1 to 0.25 using PSO. The convergence time for GA is 670,680 s, but the convergence time for PSO optimization is 356,400 s, showing that the former requires 88 % more time than the latter.

Seismic inversion is used routinely in hydrocarbon exploration and development as an effective geophysical tool. Inversion methodologies vary from relative to absolute poststack, prestack, and stochastic seismic inversion. Most of the inversion algorithms suffer from the “nonuniqueness” problem, which means there is more than one possible geologic model that fits the seismic data. Because of this, other information must be provided in the form of a low-frequency initial model. Construction of an accurate initial model requires honoring the subsurface geology of the field, adequate well control, and accurately picked horizons. The lack of all or some of these items makes it hard to build a reliable initial model, and the subsequent inversion will be subject to cumulative errors and accordingly questionable results. To mitigate these limitations, we implement the multiple linear regression algorithm to blend selected seismic attributes with seismic velocities to build a robust 3D initial model. First, we analyze the elastic properties of a given well against the seismic internal attributes and the available processing velocities. Next, we build a set of multilinear equations to convert the chosen attributes into 3D elastic properties. Finally, we smooth the results and produce the elastic initial models. No horizons are needed in this process. To prove its validity, the proposed approach was applied in the offshore Nile Delta, Egypt. The new initial models for the case studies show better consistency with the geology of the area and contain very fine details compared with the conventional initial models. Hence, the proposed approach suits exploration of deep targets as it provides a reliable link between seismic internal attributes, velocities, and elastic logs for the construction of an accurate 3D initial model even though the data may be insufficient.

This study presents a hybrid seismic inversion framework to estimate 3D acoustic impedance volumes in a geologically complex and data-limited environment. The approach integrates physics-informed pseudo-well generation, based on calibrated rock physics modeling and variogram statistics, with a deep feedforward neural network (DFNN) that maps multi-attribute seismic data to acoustic impedance while substantially reducing dependence on low-frequency background models and dense well calibration. The compact DFNN serves as a high-dimensional nonlinear mapper that learns the relationship between seismic attributes and impedance logs, a task for which it is well suited, leading to accurate predictions. We generated synthetic elastic logs (compressional and shear velocities, and bulk density) using calibrated rock physics modeling and variogram-constrained stochastic simulation to supplement real well logs. We produced a lithologically diverse and statistically coherent training dataset. This process, applied to only 3 real wells, generated a robust ensemble of 36 synthetic pseudo-wells, effectively addressing the severe data scarcity and providing sufficient training data. Seismic attributes representing amplitude, phase, and frequency characteristics are selected to facilitate the model's ability to resolve subtle geological heterogeneity. The trained DFNN is validated through a leave-one-well-out strategy yielding a cross-correlation coefficient of up to 95.4% and a normalized relative error below 1% when tested on the blind wells. Combining physical modeling with data-driven learning reduces reliance on low-frequency background models and dense calibration. Rather than replacing conventional inversion, it provides a complementary, geologically consistent, and computationally efficient approach for reliable reservoir characterization in offshore environments. Future work may focus on incorporate uncertainty quantification and volumetric convolutional networks to further improve spatial resolution and model reliability in complex subsurface settings.

The elastic properties of subsurface media, such as density and P- and S-wave velocity, are essential for understanding the mechanical behavior of the subsurface as they reflect its physical and structural characteristics. Accurate elastic properties are increasingly critical for safety assessment, geological risk analysis, and reservoir characterization in subsurface applications. However, current data acquisition methods present complementary limitations: well-log data provide high-resolution and accurate elastic properties but offers only localized information due to spatial constraints, while seismic data covers regional areas but has lower resolution and accuracy. To address these limitations, this study develops a machine learning (ML) framework that integrates seismic and well-log data to improve the prediction accuracy of subsurface elastic properties. This study introduces a novel approach that integrates ML with Gaussian processes (GP), which allows not only elastic property prediction but also quantification of predictive uncertainty. The proposed approach was validated using well-log and post-stack seismic data from the Volve dataset. Five ML models were trained to predict density and P- and S-wave velocities, with the long short-term memory+GP model achieving superior performance among all tested models. The 3-D prediction capabilities were further validated by applying the trained model to a 3-D seismic cube from the Volve field, successfully estimating 3-D elastic properties and quantifying prediction uncertainty. While the approach improves uncertainty assessment, it should be noted that uncertainty quantification is often limited in highly heterogeneous zones, as prediction reliability is inherently dataset-dependent.

The need to understand field-scale reservoir heterogeneity using seismic data requires implementing advanced solutions such as stochastic seismic inversion to go beyond the resolution of seismic data. Conventional seismic inversion techniques provide relatively low-resolution reservoir properties but do not provide quantitative estimates of the subsurface uncertainties. The objective of this study was to carry out a facies dependent geostatistical seismic inversion to generate multi-realization reservoir properties to improve the geological understanding of the two adjacent offshore fields in Abu Dhabi. An integrated approach of rock physics modelling and geostatistical inversion followed by porosity co-simulation was undertaken to characterize the spatially varying lithofacies and porosity of the complex carbonate reservoirs. Necessary checks to ensure highest quality data input included: 1) Rock physics modelling and shear sonic prediction 2) Invasion correction and production effect correction of elastic logs 3) Seismic feasibility analysis to define seismic facies and 4) Six angle stacks optimally defined to preserve AVO/AVA signature followed by AVO/AVA compliant post-stack processing. Subsequently, the joint facies driven geostatistical inversion was conducted to invert for multiple realizations high-resolution lithofacies and elastic rock properties. Finally, porosity was co-simulated and later ranked to map important geological variations. Based on the rock physics analysis, a 4 facies classification scheme (Porous Calcite, Porous Dolomite, Tight Calcite-Dolomite and Anhydrite) was adopted and used as input in the joint facies-elastic inversion. Before the geostatistical inversion, a deterministic inversion was performed that helped in refining the horizon interpretation of the surfaces used as a framework for the inversion. In geostatistical inversion, results are guided by variograms, facies, prior probability density functions, wells, inversion grid and seismic data quality. At start of the joint inversion, the parameters for inversion are defined in an unconstrained fashion aiming to obtain unbiased parameters which are blind to well control. Finally, using elastic properties constrained at the well locations, the joint geostatistical inversion was run to obtain multiple realizations of P-impedance, S-impedance, density and lithofacies. The cross-correlation between seismic and inverted synthetics was high across the whole area for all the partial angle stacks, with the lowest cross-correlation observed in the far angle stack. Lithofacies and elastic properties were used to co-simulate for porosity. The porosity results were then ranked to provide the P10, P50 and P90 models to be used for reservoir property model building. This study is an example of stochastically generating geologically consistent reservoir properties through high-resolution seismically constrained inversion results at 1ms vertical sampling. Lithofacies and elastic properties were jointly inverted, and co-simulated porosity results provided insights into high-resolution reservoir heterogeneity analysis through the ranking of equiprobable multiple realizations.

Subsurface characterization for lithological and fluid properties is important for all aspects of geophysical exploration where estimating a high-resolution elastic property through seismic inversion is vital. Starting with an initial subsurface model, computing synthetic or predicted seismic data, and matching these data with observed seismic data, seismic inversion uses an optimization process to iteratively modify the initial model until the prediction reasonably matches the observation. Routine applications of seismic inversion for subsurface reservoir characterization are currently restricted to amplitude-variation-with-angle inversion, which uses convolution as the basis for forward modeling to compute synthetic seismic data. Although computationally efficient, the inherent convolutional assumption ignores complex wave propagation effects and often fails to estimate subsurface models with sufficient accuracy. Here, we review the current state of the art for seismic inversion, and we discuss a method that uses an analytical wave equation solver for forward modeling and a global method for optimization that can overcome the current limitations of amplitude-variation-with-angle inversion. Using real seismic data, we demonstrate the accuracy of this method. Because this waveform-based method is computationally demanding, we also discuss the current advances of computational technology, including artificial intelligence that can improve its computational efficiency.

This paper introduces a novel matrix-free approach for full waveform inversion in anisotropic elastic media, incorporating density variation through the utilization of the distorted Born iterative method. This study aims to overcome the computational and storage challenges associated with the conventional matrix-based distorted Born iterative inversion method while accurately capturing the subsurface's anisotropic properties and density variations. An elastic integral equation is utilized to account for the anisotropic nature of elastic wave propagation, enabling more precise modeling of subsurface complexities. This integral equation is efficiently solved by a fast Fourier transform accelerated Krylov subspace method. Leveraging the integral equation with the distorted Born approximation, a linear relationship between the scattered wavefield and the model parameter perturbation is formulated for an integrated inversion scheme. To address the inherent ill-posedness of each linear inversion step, we formulate the normal equation with a regularization term. This is achieved by minimizing an objective function using the generalized Tikhonov method. Therefore, we can find an adequate solution for the inverse scattering problem by solving the normal equation. Following the physical interpretation of Green's function, the Fr{\'e}chet and adjoint operators within the normal equation can be employed in a matrix-free manner, allowing for significant improvement of the computational efficiency and memory demand without compromising accuracy. The proposed matrix-free full waveform inversion framework is thoroughly validated through extensive numerical experiments on synthetic datasets, showcasing its ability to reconstruct complex anisotropic structures and accurately recover stiffness parameters and density.

Inverse analysis has been utilized to understand unknown underground geological properties by matching the observational data with simulators. To overcome the underconstrained nature of inverse problems and achieve good performance, an approach is presented with embedded physics and a technique known as algorithmic differentiation. We use a physics-embedded generative model, which takes statistically simple parameters as input and outputs subsurface properties (e.g., permeability or P-wave velocity), that embeds physical knowledge of the subsurface properties into inverse analysis and improves its performance. We tested the application of this approach on four geologic problems: two heterogeneous hydraulic conductivity fields, a hydraulic fracture network, and a seismic inversion for P-wave velocity. This physics-embedded inverse analysis approach consistently characterizes these geological problems accurately. Furthermore, the excellent performance in matching the observational data demonstrates the reliability of the proposed method. Moreover, the application of algorithmic differentiation makes this an easy and fast approach to inverse analysis when dealing with complicated geological structures.

Abstract Gas hydrates and the free-gas beneath them are believed to represent a huge untapped source of energy. However, volume estimates are uncertain since little is known about the elastic properties of hydrated sediments. We apply elastic inversion to angle-dependent P-wave reflections to estimate elastic properties of hydrated sediments. A multi-channel USGS seismic line from the Blake Ridge off the east coast of North America is reprocessed to obtain migrated common-angle aperture data sets, which are then inverted for elastic properties using elastic inversion. The Pimpedance and S-impedance, and Poisson's ratio were obtained from the elastic inversion. The hydrated sediments have high elastic impedance, P-impedance, and S-impedance, but slightly lower Poisson's ratio values than those of the surrounding unhydrated sediments. The sediments containing free-gas have low elastic impedance, P-impedance and Poisson's ratio, but non-anomalous background S-impedance. Poisson's ratio can help identify the free-gas charged layers, but cannot differentiate between the hydrated sediments and non-hydrated sediments when gas hydrate concentration is low, or between the hydrated sediments and free-gas charged sediments when the gas hydrate concentration is high. Introduction Gas hydrates are ice-like solids composed of natural gases (mainly methane) and water molecules that occur under appropriate conditions of high pressure and low temperature. In many deep sea environments, these conditions are satisfied and gas hydrate can form and remain stable. According to Kvenvolden (1998), gas hydrate has also been recognized as a potential future source of energy, and as a factor in global climate change. The estimated amounts of in situ gas hydrates are 1.0?4.0x1016 m3 (Kvenvolden, 1998; Makogon, 1997), which is about a factor of 2 larger than the methane equivalent of all known recoverable and nonrecoverable fossil fuel deposits (coal, oil and natural gas) (Kvenvolden, 1998; Makogon, 1997). Methane has a global warming potential 20 times larger than the equivalent volume of carbon dioxide, so gas hydrate is a potential greenhouse agent. The estimated amounts of in-situ gas hydrates worldwide are highly speculative. Accurate estimates are difficult because knowledge of the distribution and saturation of gas hydrates in sediments is very incomplete, and the elastic properties of hydrated sediments are not well-known. A better understanding of the distribution and saturation of gas hydrate and the elastic properties of hydrated sediments from seismic data is both technically feasible, and potentially economically important. Oceanic natural gas hydrates are generally studied through seismic data that show distinctive bottom simulating reflectors (BSRs). The BSRs are seismic reflections that locally parallel the seafloor reflection, and have reverse polarity compared to the seafloor reflection. They are believed to be caused by the impedance contract associated with free-gas beneath the base of the GHSZ and mark the base of the GHSZ (Tinivella and Lodolo, 2000). Elastic impedance (EI) inversion (Connolly, 1999) has been used widely in estimating EI from multi-channel seismic data. Lu and McMechan (2004) developed algorithms to estimate P-impedance and S-impedance from elastic impedance. Vs/Vp ratio, Poisson's ratio, and Lame parameter terms are estimated from the P- and S-impedance sections. Gas hydrate and free-gas may be identified through their varying effects on different seismic attributes.

Abstract We present a case-study that compares seismic inversion methods for reservoir characterisation on the Mishrif carbonates in the Rumaila field, Iraq. The two methods interrogated - Deterministic Absolute Acoustic Impedance Inversion (DAAII) and BP's One-Dimensional Stochastic Inversion (ODiSI) - were used to predict porosity. This case-study highlights benefits, limitations and uncertainties associated with these methods. Seismic inversions produce non-unique solutions mainly due to the bandlimited input seismic, as several reservoir property profiles can result in the same seismic trace. DAAII is a widely applied Model- Based Inversion technique, which uses a Low Frequency Model (LFM, derived from well-log impedance) as input to the inversion algorithm. The algorithm uses the LFM as a starting point to seek an impedance profile that best minimises residual error between the resulting modelled synthetic, after convolution with an average of the extracted wavelets, and the input seismic trace. It accepts the inverted impedance profile when the error is minimised. DAAII typically results in a relatively smooth deterministic estimate of absolute Acoustic Impedance (AI), and the output is heavily dependent on robustness of the input LFM. This can be a severe shortcoming in reservoirs with poor well control - simply put, if the input LFM is of questionable quality (e.g. poor well-log data, sampling bias, suboptimal interpolation between wells) then the output impedance is

Caractérisation de la subsurface par inversion de forme d'onde complète 3D élastique. Etude d'un jeu de données sismiques de proche surface multi-composantes L'inversion de forme d'onde complète (FWI) est une procédure d'ajustement itératif des données entre les données observées et les données synthétiques. Les données synthétiques sont calculées en résolvant une équation d'onde. La FWI vise à reconstruire les informations détaillées des propriétés physiques du sous-sol. La méthode FWI a été développée au cours des dernières décennies, grâce à l'augmentation de la capacité de calcul et au développement de la technologie d'acquisition. La FWI a également été appliquée à à des échelles variées, allant de l'échelle globale, lithosphérique, crustale, jusqu'à la proche surface, c'est à dire quelques mètres de profondeur.Dans ce manuscrit, nous étudions l'inversion d'un jeu de données de source et de récepteur multicomposantes en utilisant un algorithme d'inversion de forme d'onde complète viscoélastique pour une cible sismique peu profonde. La cible est une ligne de tranchée enterrée à environ 1 m de profondeur. Nous présentons le pré-traitement des données, y compris une correction par déconvolution pour compenser les différentes conditions de couplage de la source et du récepteur pendant l'acquisition, ainsi qu'un procédé d'inversion en plusieurs étape pour la reconstruction des vitesses des ondes P et S. Notre mise en œuvre est basée sur une modélisation viscoélastique utilisant une discrétisation par éléments spectraux pour rendre compte avec précision de la complexité de la propagation des ondes dans cette région peu profonde. Nous illustrons la stabilité de l'inversion en partant de différents modèles initiaux, soit basés sur l'analyse des courbes de dispersion, soit des modèles homogènes cohérents avec les premières arrivées. Nous obtenons des résultats similaires dans les deux cas. Nous illustrons également l'importance de la prise en compte de l'atténuation en comparant les résultats élastiques et viscoélastiques. Les résultats 3D permettent de localiser précisément la ligne de tranchée en termes d'interprétation. Ils montrent également une autre structure de ligne de tranchée, dans une direction formant un angle de 45 degrés avec la direction de la ligne de tranchée ciblée. Cette nouvelle structure avait été précédemment interprétée comme un artefact dans les anciens résultats d'inversion 2D. L'interprétation archéologique de cette nouvelle structure est actuellement en discussion.Nous réalisons également trois expériences différentes pour comprendre l'effet des données à composantes multiples sur la FWI. La première expérience est une analyse de sensibilité de plusieurs paquets d'ondes (onde P, onde S et onde de surface) sur un modèle 3D simple basé sur une direction cartésienne de la source et du récepteur. La seconde expérience est une inversion élastique 3D basée sur des données synthétiques (utilisant la source de direction cartésienne) et de champ (utilisant la source Galperin) avec diverses combinaisons de composants. Seize combinaisons de composantes sont analysées pour chaque cas. Dans la troisième expérience, nous effectuons la décimation de l'acquisition sur la base de la deuxième expérience. Nous démontrons un avantage significatif des données multicomposantes FWI grâce à ces expériences. Dans une échelle sismique peu profonde, les inversions avec les composantes horizontales donnent une meilleure reconstruction en profondeur. En se basant sur la décimation de l'acquisition, l'inversion utilisant des données sismiques 9C fortement décimées produit des résultats similaires à l'inversion utilisant des données sismiques 1C sur l'acquisition complète.

Since the late 1970s, seismic tomography has been emerging as the pre-eminent tool for imaging the Earth’s interior from the meter to the global scale. Significant recent advances in seismic data acquisition, high-performance computing, and modern numerical methods have drastically progressed tomographic methods. Today it is technically feasible to accurately simulate seismic wave propagation through realistically heterogeneous Earth models across a range of scales. When seismic waves propagate inside the Earth and encounter structural heterogeneities with a certain scale, wave propagation speed changes, reflection, and scattering phenomena occur, and interconversions between compressional and shear waves happen. The combined effect of multiple heterogeneities produces a highly complicated wavefield recorded in the form of three-component (vertical, radial, transverse) seismograms. The ultimate objective of seismic imaging is to utilize the full information from waveforms recorded at seismic stations distributed around the globe in a broad frequency range to characterize detailed tomographic images of Earth’s interior by fitting synthetic seismograms to recorded seismograms. The full-waveform inversion technique based on adjoint and spectral-element methods can be employed to maximumly exploit the information contained in these seismic wavefield complexities to determine the fine-scale structural heterogeneities from which they originated across various orders of magnitude in

Abstract Objective/Scope The loss of some important logging data especially lithology data in old oilfields and huge cost in manpower and material on drilling and coring have brought great difficulties to the development of oilfields. A new method using deep learning for rock typing is achieved to classify and count the limited logging data, establish appropriate pore-permeability (PP) relationship and reduce the risk of reservoir prediction, which provides a concise and effective way for carbonate rock prediction. Methods, Procedures, Process In allusion to the existed problems, the paper collects, ranks the correlation between rock types and conventional logging data, which establishes a neural network model based on deep learning, divides the carbonate reservoirs into 4 types, and estimates the pore-permeability relationship for each type. Finally, a pore-permeability cloud simulation was performed based on the geo-statistical inversion to set up a high-precision reservoir static model with perfect well-seismic tie. The reliable permeability property can be obtained which helps to accurately depict the spatial distribution of the reservoirs. Results, Observations, Conclusions The carbonate reservoir of M formation for H oilfield in the Middle East is of complex pore structure with strong heterogeneity and poor relationship of the pore-permeability (PP). The logs DT, GR, Density, Porosity as the input features of deep learning is optimized to train neural network models, wh

simulation, seismic inversion techniques combine well and seismic data to produce multiple equally plausible 3D models of the elastic properties of the reservoir In the oil and gas industry, reservoir modeling involves the construction of a computer model of a petroleum reservoir, for the purposes of improving estimation of reserves and making decisions regarding the development of the field, predicting future production, placing additional wells and evaluating alternative reservoir management scenarios. A reservoir model represents the physical space of t Geological models are created by geologists and geophysicists and aim to provide a static description of the reservoir, prior to production. Reservoir simulation models are created by reservoir engineers and use finite difference methods to simulate the flow of fluids within the reservoir, over its production lifetime. Sometimes a single "shared earth model" is used for both purposes. More commonly, a geological model is constructed at a relatively high (fine) resolution. A coarser grid for the reservoir simulation model is constructed, with perhaps two orders of magnitude fewer cells. Effective values of attributes for the simulation model are then derived from the geological model by an upscaling process. Alternatively, if no geological model exists, the attribute values for a simulation model may be determined by a process of sampling geological maps. Uncertainty in the true values of the reservoir properties is sometimes investigated by constructing several different realizations of the sets of attribute values. The behaviour of the resulting simulation models can then indicate the associated level of economic uncertainty. The phrase "reservoir characterization" is sometimes used to refer to reservoir modeling activities up to the point when a simulation model is ready to simulate the flow of fluids. Commercially available software…

seismic inversion is the process of transforming seismic reflection data into a quantitative rock-property description of a reservoir. Seismic inversion In geophysics (primarily in oil-and-gas exploration/development), seismic inversion is the process of transforming seismic reflection data into a quantitative rock-property description of a reservoir. Seismic inversion may be pre- or post-stack, deterministic, random or geostatistical; it typically includes other reservoir measurements such as well logs and cores. In geophysics (primarily in oil-and-gas exploration/development), seismic inversion is the process of transforming seismic reflection data into a quantitative rock-property description of a reservoir. Seismic inversion may be pre- or post-stack, deterministic, random or geostatistical; it typically includes other reservoir measurements such as well logs and cores. pre-stack or post-stack seismic resolution or well-log resolution The combination of these categories yields four technical approaches to the inversion problem, and the selection of a specific technique depends on the desired objective and the characteristics of the subsurface rocks. Although the order presented reflects advances in inversion techniques over the past 20 years, each grouping still has valid uses in particular projects or as part of a larger workflow.