Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “3D Geologic Model”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 199 records · Page 11

Hydrogeologic Model Based on Reduced Complexity

This data set provides the basic hydrogeologic model used for development of the TOUGH2 reservoir modeling results for the Kimberlina site. The hydrogeologic model is based on the static geologic model provided as a separate data set but has some simplifications regarding the geologic layers. Some geologic layers have been merged for the hydrogeologic model because the hydrogeologic properties were very similar. The data set contains several files. The word file explains the structure of the data files. The pdf contains a table with the hydrogeologic properties. The data set also includes 15 thickness files for 15 geological layers (including 3 Vedder sand layers and 3 Vedder shale layers), three elevations files (with reference elevation at ground surface/model top at 9000 m), and one fault trace file for seven faults. Most users may only need 2-3 faults close to the injection well (also see the figure). The arrangement of these layer thickness and elevations files can be see from the three layer files (for Tecplot). To develop a grid structure, start with the Vedder top elevation (REFVedderTop.dat) and add/subtract the thicknesses of layers above/below.

3D Groundwater Model↗

2D-to-3D image translation of complex nanoporous volumes using generative networks

Image-based characterization offers a powerful approach to studying geological porous media at the nanoscale and images are critical to understanding reactive transport mechanisms in reservoirs relevant to energy and sustainability technologies such as carbon sequestration, subsurface hydrogen storage, and natural gas recovery. Nanoimaging presents a trade off, however, between higher-contrast sample-destructive and lower-contrast sample-preserving imaging modalities. Furthermore, high-contrast imaging modalities often acquire only 2D images, while 3D volumes are needed to characterize fully a source rock sample. In this work, we present deep learning image translation models to predict high-contrast focused ion beam-scanning electron microscopy (FIB-SEM) image volumes from transmission X-ray microscopy (TXM) images when only 2D paired training data is available. We introduce a regularization method for improving 3D volume generation from 2D-to-2D deep learning image models and apply this approach to translate 3D TXM volumes to FIB-SEM fidelity. We then segment a predicted FIB-SEM volume into a flow simulation domain and calculate the sample apparent permeability using a lattice Boltzmann method (LBM) technique. Results show that our image translation approach produces simulation domains suitable for flow visualization and allows for accurate characterization of petrophysical properties from non-destructive imaging data.

58 GEOSCIENCES↗

Starting Conditions for Hydrothermal Systems Underneath Martian Craters: Hydrocode Modeling

Mars is the most Earth-like of the Solar System s planets, and the first place to look for any sign of present or past extraterrestrial life. Its surface shows many features indicative of the presence of surface and sub-surface water, while impact cratering and volcanism have provided temporary and local surface heat sources throughout Mars geologic history. Impact craters are widely used ubiquitous indicators for the presence of sub-surface water or ice on Mars. In particular, the presence of significant amounts of ground ice or water would cause impact-induced hydrothermal alteration at Martian impact sites. The realization that hydrothermal systems are possible sites for the origin and early evolution of life on Earth has given rise to the hypothesis that hydrothermal systems may have had the same role on Mars. Rough estimates of the heat generated in impact events have been based on scaling relations, or thermal data based on terrestrial impacts on crystalline basements. Preliminary studies also suggest that melt sheets and target uplift are equally important heat sources for the development of a hydrothermal system, while its lifetime depends on the volume and cooling rate of the heat source, as well as the permeability of the host rocks. We present initial results of two-dimensional (2D) and three-dimensional (3D) simulations of impacts on Mars aimed at constraining the initial conditions for modeling the onset and evolution of a hydrothermal system on the red planet. Simulations of the early stages of impact cratering provide an estimate of the amount of shock melting and the pressure-temperature distribution in the target caused by various impacts on the Martian surface. Modeling of the late stage of crater collapse is necessary to characterize the final thermal state of the target, including crater uplift, and distribution of the heated target material (including the melt pool) and hot ejecta around the crater.

Pierazzo, E.↗

A robust framework for frictional fault contact in geological formations using a stabilized augmented Lagrangian approach

Numerical simulations are essential to evaluate the performance and safety of engineered subsurface systems such as geological carbon storage sites, enhanced geothermal fields, and oil and gas reservoirs. A key challenge lies in accurately modeling the frictional contact behavior along fault surfaces. This problem involves inequality constraints that arise from the physics of frictional slip, requiring specialized numerical methods to handle the resulting highly nonlinear and path-dependent behavior. Here, in this work, we address this challenge using an Augmented Lagrangian Method (ALM) implemented via the Uzawa algorithm. The formulation employs mixed finite element spaces, combining low-order piecewise linear displacements within the 3D domain cells with piecewise constant tractions defined on the fault surfaces. Furthermore, to ensure stability and satisfy the inf-sup condition, the discrete displacement space is enriched with face bubble functions on both sides of the contact interfaces. This approach offers several advantages over other stabilization techniques that rely on additional terms, and it integrates naturally in the Uzawa framework.

58 GEOSCIENCES↗

TOGA: A TOUGH code for modeling three-phase, multi-component, and non-isothermal processes involved in CO 2 -based Enhanced Oil Recovery

TOGA is a numerical reservoir simulator for modeling non-isothermal flow and transport of water, CO 2 , multicomponent oil, and related gas components for applications including CO 2 -enhanced oil recovery (CO 2 -EOR) and geologic carbon sequestration in depleted oil and gas reservoirs. TOGA uses an approach based on the Peng-Robinson equation of state (PR-EOS) to calculate the thermophysical properties of the gas and oil phases including the gas/oil components dissolved in the aqueous phase, and uses a mixing model to estimate the thermophysical properties of the aqueous phase. The phase behavior (e.g., occurrence and disappearance of the three phases, gas + oil + aqueous) and the partitioning of non-aqueous components (e.g., CO 2 , CH 4 , and n-oil components) between coexisting phases are modeled using K-values derived from assumptions of equal-fugacity that have been demonstrated to be very accurate as shown by comparison to measured data. Models for saturated (water) vapor pressure and water solubility (in the oil phase) are used to calculate the partitioning of the water (H 2 O) component between the gas and oil phases. All components (e.g., CO 2 , H 2 O, and n hydrocarbon components) are allowed to be present in all phases (aqueous, gaseous, and oil). TOGA uses a multiphase version of Darcy’s Law to model flow and transport through porous media of mixtures with up to three phases over a range of pressures and temperatures appropriate to hydrocarbon recovery and geologic carbon sequestration systems. Transport of the gaseous and dissolved components is by advection and Fickian molecular diffusion. New methods for phase partitioning and thermophysical property modeling in TOGA have been validated against experimental data published in the literature for describing phase partitioning and phase behavior. Flow and transport has been verified by testing against related TOUGH2 EOS modules and CMG. The code has also been validated against a CO 2 -EOR experimental core flood involving flow of three phases and 12 components. Results of simulations of a hypothetical 3D CO 2 -EOR problem involving three phases and multiple components are presented to demonstrate the field-scale capabilities of the new code. This user guide provides instructions for use and sample problems for verification and demonstration.

97 MATHEMATICS AND COMPUTING↗

Improved image of intrusive bodies at Newberry Volcano, Oregon, based on 3D gravity modelling

Beneath Newberry Volcano is one of the largest geothermal heat reservoirs in the western United States and it has been extensively studied for the last 40 years. Several magmatic intrusions have been recognized at depths between 2.5 and 8 km and some of them identified as suitable targets for enhanced geothermal energy and tested during two previous EGS campaigns. These subsurface structures have been intersected by three deep wells and imaged by various geophysical methods including seismic tomography and magnetotellurics. Although three high quality gravity surveys were completed between 2006 and 2010 as part of various projects, a complete synthesis and interpretation of the gravity data has not yet been performed. Regional gravity data also exist in the vicinity of the Newberry volcano and have been added to these surveys to constitute a dataset with a total of 1418 gravity measurements. When coupled with existing geologic and geophysical data and models, this new gravity dataset provides important constraints on the depth and contours of the magmatic bodies previously identified by other methods and thus greatly contributing to facilitate any future drilling and stimulation works. Using the initial structures discovered by seismic tomography, inversion of gravity data has been performed. Shape, density values and depths of various bodies were allowed to vary and three main bodies have been identified. Densities of the middle and lower intrusive bodies (~2.6-2.7 g/cm3) are consistent with rhyolite, basalt or granites. Modeled density of the near-surface caldera body match that of a low density tephra material and the density of the shallow ring structures contained in the upper kilometer correspond to that of welded tuff or low-density rhyolites. Modeled bodies are in reality a composite of thin layers; however, average densities of the modeled gravity bodies are in good agreement with the density log obtained in one well located on the western flank (well 55-29). Final gravity data residuals show that most of the observed gravity anomalies at the surface can be explained by the modeled gravity bodies and are consistent with other site characterization information.

Cascade Range↗

Multi-Scale Seismic Measurements for Site Characterization and CO2 Monitoring in an Enhanced Oil Recovery/Carbon Capture, Utilization, and Sequestration Project, Farnsworth Field, Texas

To address the challenges of climate change, significantly more geologic carbon sequestration projects are beginning. The characterization of the subsurface and the migration of the plume of supercritical carbon dioxide are two elements of carbon sequestration that can be addressed through the use of the available seismic methods in the oil and gas industry. In an enhanced oil recovery site in Farnsworth, TX, we employed three separate seismic techniques. The three-dimensional (3D) surface seismic survey required significant planning, design, and processing, but produces both a better understanding of the subsurface structure and a three-dimensional velocity model, which is essential for the second technique, a timelapse vertical seismic profile, and the third technique, cross-well seismic tomography. The timelapse 3D Vertical Seismic Profile (3D VSP) revealed both significant changes in the reservoir between the second and third surveys and geo-bodies that may represent the extent of the underground carbon dioxide. The asymmetry of the primary geo-body may indicate the preferential migration of the carbon dioxide. The third technique, cross-well seismic tomography, suggested a strong correlation between the well logs and the tomographic velocities, but did not observe changes in the injection interval.

Energy & Fuels↗

Determining the 3D Subsurface Density Structure of Taurus Littrow Valley Using Apollo 17 Gravity Data

Surface gravity surveys can detect subsurface density variations that can reveal subsurface geologic features. In 1972, the Apollo 17 (A17) mission conducted the Traverse Gravimeter Experiment (TGE) using a gravimeter that measured the local gravity field near Taurus Littrow Valley (TLV), located on the south-eastern rim of the Serenitatis basin. TLV is hypothesized to be a basaltfilled radial graben resulting from the impact that formed Mare Serenitatis. It is bounded by both the North and South Massifs (NM and SM) as well as other smaller mountains to the East that are thought to be mainly composed of brecciated highland material. The TGE is the first and only successful gravity survey on the surface of the Moon. Other more recent satellite surveys, such as NASA's Gravity Recovery and Interior Laboratory (GRAIL) mission (2011- 2012), have produced the best global gravity field to date (approx. 13km resolution). However, these satellite surveys are not sensitive enough to detect fine-scale (<1km) lunar subsurface structures. This underscores the value of the data collected at the surface by A17. In the original analysis of the data a 2D forward-modelling approach was used to derive a thickness of the subsurface basalt layer of 1.0 km by assuming a simple flat-faced rectangular geometry and using densities derived from Apollo lunar samples. We are investigating whether modern 3D modelling techniques in combination with high-resolution topographical and image datasets can reveal additional fine-scale subsurface structure in TLV.

TGE↗

De-Risking Exploration for Geothermal Plays in Magmatic Environments Through Open-Source Tools: An Open-Source Python Framework for 2D and 3D Play Fairway Analysis

The De-Risking Exploration for Geothermal Plays in Magmatic Environments (DEEPEN) project seeks to accelerate superhot geothermal development by reducing exploration risk through advanced open-source modeling tools. This work presents a novel Python-based framework, geoPFA, for conducting 2D and 3D play fairway analysis (PFA) tailored to superhot geothermal systems. Building on previous methodologies, the framework integrates thermo-hydro-mechanical-chemical simulation outputs from TReactMech, resulting in improved representation of subsurface properties that are critical to superhot resource producibility. The workflow has been applied to the Nesjavellir field in Iceland, a candidate site for the third Iceland Deep Drilling Project's superhot production scenarios. This application demonstrates the value of modular, transparent, and extensible workflows for integrating geological, geophysical, and simulation-derived datasets in high-enthalpy environments. Preliminary results indicate favorable zones consistent with known hydrothermal activity and suggest possible upflow from the Hengill volcanic system. The geoPFA library is publicly available, offering a scalable and reproducible approach to geothermal exploration across varied geological contexts.

15 GEOTHERMAL ENERGY↗

A robust deep learning workflow to predict multiphase flow behavior during geological C O 2 sequestration injection and Post-Injection periods

Simulation of multiphase flow in porous media is essential to manage the geologic CO 2 sequestration (GCS) process, and physics-based simulation approaches usually take prohibitively high computational cost due to the nonlinearity of the coupled physics. This paper contributes to the development and evaluation of a deep learning workflow that accurately and efficiently predicts the temporal-spatial evolution of pressure and CO 2 plumes during injection and post-injection periods of GCS operations. Based on a Fourier Neural Operator, the deep learning workflow takes input variables or features including rock properties, well operational controls and time steps, and predicts the state variables of pressure and CO 2 saturation. To further improve the predictive fidelity, separate deep learning models are trained for CO 2 injection and post-injection periods due to the difference in primary driving force of fluid flow and transport during these two phases. We also explore different combinations of features to predict the state variables. We use a realistic example of CO 2 injection and storage in a 3D heterogeneous saline aquifer, and apply the deep learning workflow that is trained from physics-based simulation data and emulate the physics process. Through this numerical experiment, we demonstrate that using two separate deep learning models to distinguish post-injection from injection period generates the most accurate prediction of pressure, and a single deep learning model of the whole GCS process including the cumulative injection volume of CO 2 as a deep learning feature, leads to the most accurate prediction of CO 2 saturation. For the post-injection period, it is key to use cumulative CO 2 injection volume to inform the deep learning models about the total carbon storage when predicting either pressure or saturation. The deep learning workflow not only provides high predictive fidelity across temporal and spatial scales, but also offers a speedup of 250 times compared to full physics reservoir simulation, and thus will be a significant predictive tool for engineers to manage the long-term process of GCS.

58 GEOSCIENCES↗

Mortandad Canyon seismic modeling (Progress Report)

The goal of this project is to assess site amplification for LANL through modeling of site and topography effects. Such assessment is difficult due to the lack of data documenting site effects for large normal faulting earthquakes such as the Pajarito Fault System that was found to be governing the seismic hazard at LANL. LANL site is also characterized by unique geologic settings, a series of mesa at the edge of the huge volcanic complex of the Jemez caldera. The latter has produced a series of successive and intertwined deposits of volcanic flow and sedimentary layer. Such layering of material with different elastic properties is expected to create large amplification effects magnified further by the presence at the surface of tuff materials. In 2014, a seismic campaign was performed in the Mortandad Canyon with a network of sensors on and near a local mesa. That campaign was designed to provide experimental data to estimation and prediction of amplification and site effects at LANL. An analysis on the seismic noise records of that campaign shows evidence of amplification of the seismic amplitude for the stations on the top of the mesa. Our goal is to reproduce the effect observed (both retrieving the resonance frequencies and amplification levels). We also plan to test if amplification effects observed for seismic noise can be extended to shaking produced by earthquakes. We are using a 3D modeling tool based on the Spectral Element Method (SEM) that is a particular type of high-order finite element modeling. Two types of seismic sources will be considered: (1) earthquakes for purpose of seismic hazard evaluation, (2) ambient noise. We show 2 series of numerical computations: (1) seismic waves are generated by random sources to simulate ambient noise. (2) seismic waves are planar waves to simulate a distant earthquake. We plan in the future to model rupture scenario on the Pajarito Fault System.

58 GEOSCIENCES↗

Geophysical Imaging of Geothermal Systems Spanning Various Geologic Settings

Despite the wide range of geologic settings in which geothermal systems reside, many systems share characteristics that can be targeted by geophysical techniques. Over the last few years, the authors have used potential field and magnetotelluric (MT) methods to jointly image and characterize numerous geothermal systems in various geologic and tectonic settings. This includes systems in the Basin and Range, Washington Cascade Range, Northeastern Oregon, and the Snake River Plain. The principal geothermal targets are areas of elevated temperatures near fault interactions that promote permeability and fluid flow. Gravity, magnetic, and MT data are all sensitive to different physical parameters that together form a complementary set of tools for characterizing subsurface structure and identifying favorable geothermal targets. In general, potential fields delineate important geologic structures and MT locates fluids and zones of hydrothermal alteration. Joint interpretation of these complementary methods, along with other geologic information, provides a robust framework for geothermal prospecting and understanding a hydrothermal system. Specifically, common geothermal targets lie along potential field gradients and are associated with an electrically conductive anomaly that resides in the less-dense, less-magnetic material.

15 GEOTHERMAL ENERGY↗

Comparative Assessment of U-Net-Based Deep Learning Models for Segmenting Microfractures and Pore Spaces in Digital Rocks

Segmentation of high-resolution X-ray microcomputed tomography (µCT) images is crucial in digital rock physics (DRP), affecting the characterization and analysis of microscale phenomena in the porous media. The complexity of geological structures and nonideal scanning conditions pose significant challenges to conventional image segmentation approaches. Motivated by the recent increasing popularity of deep learning (DL) techniques in image processing, this work undertakes a comparative study of DL models, specifically U-Net and its variants, for segmenting multiple targets with distinguished features in digital rocks, including discrete fracture networks (DFNs), pore spaces, and solid rock. Particularly, DFNs have a smaller volumetric fraction over others, bringing in a substantial challenge of imbalanced segmentation. The primary focus is to evaluate the architecture and feature enhancement strategies of various DL models, including U-Net, attention U-Net, residual U-Net, U-Net++, and residual U-Net++. The models were designed as 2.5D, utilizing a central 2D image and its two adjacent upper and lower 2D images as input to provide a pseudo-3D context. In addition, because the ground truth of segmentation was unknown for real-world digital rocks, we created a benchmark data set following the inverse operations of segmentation. The data synthesis started from the label images (i.e., solid rock, pore spaces, and DFNs), followed by simulating partial volume blurring, adding random background noise, and introducing ring artifacts to mimic real raw X-ray µCT images. The data set, which included various rock types (i.e., sandstone and artificial data), scanning resolution, and magnitudes of noise and artifacts, was divided into training and testing data sets with a 90% and 10% ratio, respectively. Moreover, in addition to the conventional pixel-wise evaluation metrics, the physics-based metric of the lattice-Boltzmann method (LBM) simulated permeability provided more comprehensive assessments. The results demonstrated that the residual connections, nested architectures, and redesigned skip connections contribute to the model performance and give the residual U-Net++ the highest accuracy. The improvements were mainly on the boundaries and small targets, especially the DFNs, which dominate the interconnectivity and therefore affect the permeability greatly. This study also rigorously evaluated the efficiency and generalization of each model, demonstrating that the sophisticated architectures achieved excellent practicability and maintained robust performance on completely unseen data, ensuring their suitability for diverse and challenging DRP applications.

58 GEOSCIENCES↗

Exploring for Superhot Geothermal Targets in Magmatic Settings: Developing a Methodology

This paper presents preliminary results from a subset of work carried out as part of a multinational research project entitled DErisking Exploration for multiple geothermal Plays in magmatic ENvironments (DEEPEN). One objective of DEEPEN is to develop a customized approach to exploration for superhot geothermal plays in magmatic systems. This paper summarizes key geologic components, risk factors, and exploration methods for geothermal plays in magmatic settings based on a review and comparative analysis of international training sites. As part of a Play Fairway Analysis (PFA) approach to exploring for multiple play types in a single magmatic system, training data were compiled and weights assigned to various evidence layers. Two different approaches for weighting exploration datasets are described in this paper - one based on expert opinions and the other using statistical learning. Weights produced by both approaches will be input into a 3D PFA workflow that combines multiple exploration datasets to generate 3D geothermal favorability models, which will be applied to two international demonstration sites.

GEOTHERMAL ENERGY↗

The character and amplitude of ‘discontinuous’ bottom-simulating reflections in marine seismic data

Bottom-simulating reflections (BSRs) identified in seismic data are well documented; and are commonly interpreted to indicate the presence of gas hydrates along continental margins, as well as to estimate regional volumes of gas hydrate. A BSR is defined as a reflection that sub-parallels the seafloor but is opposite in polarity and cross-cuts dipping sedimentary strata. BSRs form as a result of a strong negative acoustic impedance contrast. BSRs, however, are a diverse seismic phenomena that manifest in strikingly contrasting ways in different geological settings, and in different seismic data types. We investigate the characteristics of BSRs, using conventional and high resolution, 2D and 3D seismic data sets in three locations: the Terrebonne and Orca Basins in the Gulf of Mexico, and Blake Ridge on the US Atlantic Margin. The acquisition geometry and frequency content of the seismic data significantly impact the resultant character of BSRs, as observed with depth and amplitude maps of the BSRs. Furthermore, our amplitude maps reinforce the concept that the BSR represents a zone, over which the transition from hydrate to free gas occurs, as opposed to the conventional model of the BSR occurring at a single interface. Our results show that a BSR can be mapped in three dimensions but it is not spatially continuous, at least not at the basin scale. Rather, a BSR manifests itself as a discontinuous, or patchy, reflection and only at local scales is it continuous. We suggest the discontinuous nature of BSRs is the result of variable saturation and distribution of free gas and hydrate, acquisition geometry and frequency content of the recorded seismic data. Lastly, the commonly accepted definition of a BSR should be broadened with careful consideration of these factors, to represent the uppermost extent of enhanced amplitude at the shallowest occurrence of free gas trapped by overlying hydrate-bearing sediments.

03 NATURAL GAS↗

Coupling fracture facies with in-situ permeability measurements to generate stochastic simulations of tight carbonate aquifer properties: Example from the Lower Cretaceous aquifer, Northern Provence, SE France

The relationships between fracture facies and permeability distribution have been explored using a well-constrained 3D implicit structural model of hemipelagic low porosity/low permeability carbonate from the Northern Provence region, SE France. Fracture and permeability facies were determined using an extensive dataset of 99 hydrogeological wells. Data processing was undertaken using a step-by-step approach, involving: (i) identification of the fracture facies based on well logs detailing structure and stratigraphy; (ii) determination of permeability facies from the a priori correlation between the dimension of the hydraulic radius of influence (deduced from slug test analyses) and the type of reservoir heterogeneity (fissure, fracture, fault zone, etc.); (iii) three dimensional plot of fracture and permeability facies in the geological model using a variographic analysis of data. Thirty-three sequential indicator simulations (SIS) based on geostatistic analyses were realised on both fracture and permeability facies. Finally, a connectivity algorithm was developed to compute the probability of connection between selected infiltration areas and the major aquifer springs via moderate- to high-permeability geological bodies. Key results are: (i) fault zones have a greater role in controlling permeability facies distribution than on fracture facies repartition; (ii) there is little correlation between permeability and fracture facies distributions; (iii) connectivity results highlight the compartmentalisation of aquifers in the Cadarache area, the extensions of permeable geological bodies being limited by the N130 faults.

58 GEOSCIENCES↗

Real Time 3D Observations of Portland Cement Carbonation at CO2 Storage Conditions

Depleted oil reservoirs are considered a viable solution to the global challenge of CO 2 storage. A key concern is whether the wells can be suitably sealed with cement to hinder the escape of CO 2 . Under reservoir conditions, CO 2 is in its supercritical state, and the high pressures and temperatures involved make real-time microscopic observations of cement degradation experimentally challenging. Here, we present an in situ 3D dynamic X-ray micro computed tomography (μ-CT) study of well cement carbonation at realistic reservoir stress, pore-pressure, and temperature conditions. The high-resolution time-lapse 3D images allow monitoring the progress of reaction fronts in Portland cement, including density changes, sample deformation, and mineral precipitation and dissolution. By switching between flow and nonflow conditions of CO 2 -saturated water through cement, we were able to delineate regimes dominated by calcium carbonate precipitation and dissolution. For the first time, we demonstrate experimentally the impact of the flow history on CO 2 leakage risk for cement plugging. In-situ μ-CT experiments combined with geochemical modeling provide unique insight into the interactions between CO 2 and cement, potentially helping in assessing the risks of CO 2 storage in geological reservoirs.

58 GEOSCIENCES↗

Bayesian seismic inversion based on rock-physics prior modeling for the joint estimation of acoustic impedance, porosity and lithofacies

We propose a Bayesian approach for seismic inversion to estimate acoustic impedance, porosity and lithofacies within the reservoir conditioned to post-stack seismic and well data. The link between elastic and petrophysical properties is given by a joint prior distribution for the logarithm of impedance and porosity, based on a rock-physics model. The well conditioning is performed through a background model obtained by well log interpolation. Two different approaches are presented: in the first approach, the prior is defined by a single Gaussian distribution, whereas in the second approach it is defined by a Gaussian mixture to represent the well data multimodal distribution and link the Gaussian components to different geological lithofacies. The forward model is based on a linearized convolutional model. For the single Gaussian case, we obtain an analytical expression for the posterior distribution, resulting in a fast algorithm to compute the solution of the inverse problem, i.e. the posterior distribution of acoustic impedance and porosity as well as the facies probability given the observed data. For the Gaussian mixture prior, it is not possible to obtain the distributions analytically, hence we propose a Gibbs algorithm to perform the posterior sampling and obtain several reservoir model realizations, allowing an uncertainty analysis of the estimated properties and lithofacies. Both methodologies are applied to a real seismic dataset with three wells to obtain 3D models of acoustic impedance, porosity and lithofacies. The methodologies are validated through a blind well test and compared to a standard Bayesian inversion approach. Using the probability of the reservoir lithofacies, we also compute a 3D isosurface probability model of the main oil reservoir in the studied field.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗