Dissertations / Theses on the topic 'Petrophysical interpretation'

39 thin sections and petrophysical log data from the Skalle well in the Hammerfest Basin, in the southwestern part of the Barents Sea, have been studied to interpret lithology, and diagenesis and their effect on the reservoir quality, and to compare reservoir properties of the different reservoir units. Petrophysical log data have been calibrated for reservoir description in cases where core material is not available. The studied formations are comprised by the Stø, Fuglen, Hekkingen, Knurr, Kolje and the Lower Kolmule Formations. The Knurr and Kolje Formations have been identified and interpreted only by wire line logs, as core material was not available for those intervals.The Lower Kolmule Formation of sandstones of lithic greywacke composition, and the Stø Formation with sandstones of subarkosic arenite composition are considered as possible reservoir rocks. All the formations are water filled which is reflected by the low resistivity logs responses. The mature sandstones of the Stø Formation show high reservoir quality (high porosity and permeability) compared to the Lower Kolmule Formation. The Hekkingen Formation is a potential source rock for the Lower Kolmule Formation, as well as a seal (cap rock) for the Stø Formation. Cementation, dissolution, compaction, clay mineral authigenesis and stylolitization are the most significant diagenetic processes affecting the reservoir quality. Some other type of processes such as glauconitization and bioturbation are also common in the studied well.

Magister Scientiae - MSc (Earth Science)
The Orange basin is relatively underexplored with 1 well per every 4000km2 with only the Ububhesi gas field discovery. Block 1 is largely underexplored with only 3 wells drilled in the entire block and only well A?F1 inside the 1500km2 3?D seismic data cube, acquired in 2009. This study is a reservoir characterization of well A?F1, utilising the acquired 3?D seismic data and re?analysing and up scaling the well logs to create a static model to display petrophysical properties essential for reservoir characterization. For horizon 14Ht1, four reservoir zones were identified, petro?physically characterized and modelled using the up scaled logs. The overall reservoir displayed average volume of shale at 24%, good porosity values between 9.8% to 15.3% and permeability between 2.3mD to 9.5mD. However, high water saturation overall which exceeds 50% as per the water saturation model, results in water saturated sandstones with minor hydrocarbon shows and an uneconomical reservoir.

A seqüência mesozóica da bacia de Ucayali é a maior produtora de gás e condensado do Peru. A área do trabalho, denominada Grande Camisea, fica na parte sul da bacia e, na atualidade, pertence à companhia Plupetrol Peru Corporation. Neste trabalho, foi aplicado um método de interpretação de perfis de indução em um poço petrolífero no sector San Martin do campo Camisea (QuickLook Interpretation method). O objetivo consiste na caracterização do reservatório de San Martín utilizando um método de interpretação rápida de perfis elétricos e, assim, fornecer uma visão geral no entendimento de parâmetros de poços e reservatórios, de zonas produtivas e suas características petrofísicas de porosidade e de saturação do óleo. Para validar a interpretação, foram utilizadas a descrição geológica de testemunhos e amostras de calha, descrição e informação do sistema petrolífero do campo e a geologia regional da zona de interesse da bacia. Desta forma, foi possível apresentar uma comparação entre os valores obtidos através dos métodos detalhados executados pela Pluspetrol e o método rápido de interpretação aplicado aqui, assim como o desvio entre ambos os resultados.
The Mesozoic sequence of the Ucayali basin is the main producer of gas and condensate of Peru. The work area is called Gran Camisea, located in the south part of the basin, and, in the present time, belongs to the company Plupetrol Peru Corporation. In this work, a well log interpretation method was used in a gas well in San Martin area, part of the Camisea field. The goal is the characterization of the reservoir of San Martín using a Quick Look log interpretation method, and thus to supply a general view in the understanding of well and reservoirs parameters, productive zones and its petrophysics characteristics of porosity and saturation. To validate the interpretation, besides using the geologic description of well cores and cutting sampling, it was used the description and information of the petroleum system of Camisea gas field and its regional geology. It was possible to present a comparison between Pluspetrol values, obtained through detailed methods, and those from the Quick Look log interpretation method used here, as well as an analysis of convergence between both results.

>Magister Scientiae - MSc
The evaluation of petroleum reservoirs has shifted from single approach to an integrated approach. The integration, analysis and understanding of all available data from the well bore and creating property models is an exceptional way to characterize a reservoir. Formulating, implementing, and demonstrating the applicability of the joint inversion of seismic and well-bore related observations, and the use of information about the relationship between porosity and permeability as the key parameters for identifying the rock types and reservoir characterization is a vital approach in this study. Correlating well and seismic data, potential reservoirs can be delineated and important horizons (markers) can be pointed out to better characterize the reservoir. This thesis aims to evaluate the potential petroleum reservoirs of the Wells K-A1, K-A2, K-A3 and K-H1 of the Shungu Shungu field in the Orange Basin through the integration and comparison of results from core analysis, wireline logs and seismic and attempt to produce a good reservoir model and by additionally utilizing Petrophysics and seismic and trying to better understand why the area has dry wells. Different rock types that comprise reservoir and non reservoirs are identified in the study and five Facie types are distinguished. Tight, fine grained sandstones with low porosity values ranging from 3% - 6% where dominant in the targeted sandstone layers. Porosity values ranging from 11% - 18% where identified in the massive sandstone lithologies which where hosted by Well’s K-A2 and K-A3. Low permeability values reaching 0.1mD exist throughout the study area. Areas with high porosity also host high water saturation values ranging from 70 – 84%. An improvement in the porosity values at deeper zones (3700m -3725m) and is apparent. Poroperm plots exhibit quartz cemented sandstones and density with neutron plot suggest that the sandstones in the area contain quarts and dolomite mineralization.Well K-A3, consist of a cluster by quartzitic sandstone, meaning there is a large amount of sandstone present. There are apparent high porosity values around the sandstone. What is apparent from this plot is that there are many clusters that are scattered outside the chart. This could suggest some gas expulsions within this Well. Sandstones within the 14B2t1 to 14At1 interval are less developed in the vicinity covered by well K-A2 than at the K-A1 well location. The main targeted sandstones belong to the lower cretaceous and lie just below 13At1. The four wells drilled in this area are dry wells. The areas/blocks surrounding this area have shown to possess encouraging gas shows and a comparative study could reveal better answers. At deeper zones of the well at an interval of 5350m -5750m, there are more developed sandstones with good porosity values. The volume of shale is low and so is the water saturation. The main target sandstones in the study area are the Lower Cretaceous sandstones which are at an interval 13At1. These sandstones are not well developed but from the property model of the target surface it can be seen that the porosity values are much more improved than the average values applied on all the zones on the 3D grid.

Natural gas hydrates store a substantial portion of the Earth's organic carbon, although their occurrence is restricted by thermobaric boundaries and the availability of methane-rich fluids. The complexity of geological systems and the multiphase flow processes promoting hydrate formation can result in a mismatch between the predicted and the observed hydrate distribution. The purpose of this research is to achieve a better comprehension of the factors that influence the distribution of gas hydrates and the mechanism of fluid movements beneath and across the gas hydrate stability zone (GHSZ). Therefore, this study integrates seismic, petrophysical and geochemical data from different gas hydrate provinces. This work provides evidence that hydrates can occur below bottom-simulating reflectors, in the presence of sourcing thermogenic hydrocarbons. The relationship between fluid-escape pipes and gas hydrates is further explored, and pipe-like features are suggested to host a significant volume of hydrates. The host lithology also represents a critical factor influencing hydrate and free gas distribution and, in evaluating a natural gas hydrate system, needs to be considered in conjunction with the spatial variability in the methane supply. The three-dimensional distribution of gas hydrate deposits in coarse-grained sediments, representing the current target for hydrate exploration, is shown to be correlated with that of the underlying free gas zone, reflecting sourcing mechanisms dominated by a long-range advection. In such systems, the free gas invasion into the GHSZ appears controlled by the competition between overpressure and sealing capacity of the gas hydrate-bearing sediments. Globally, the thickness of the free gas zones is regulated by the methane supply and by different multi-phase flow processes, including fracturing, capillary invasion and possibly diffusion. In conclusion, this research indicates that geological, fluid flow and stability factors interweave at multiple scales in natural gas hydrate systems.

The estimation of porosity, water saturation, kerogen concentration, and mineral composition is an integral part of unconventional shale reservoir formation evaluation. Porosity, water saturation, and kerogen content determine the amount of hydrocarbon-in-place while mineral composition affects hydro-fracture generation and propagation. Effective hydraulic fracturing is a basic requirement for economically viable flow of gas in very-low permeability shales. Brittle shales are favorable for initiation and propagation of hydraulic fracture because they require marginal or no plastic deformation. By contrast, ductile shales tend to oppose fracture propagation and can heal hydraulic fractures. Silica and carbonate-rich shales often exhibit brittle behavior while clay-rich shales tend to be ductile. Many operating companies have turned their attention to neutron capture gamma-ray spectroscopy (NCS) logs for assessing in-situ mineral composition. The NCS tool converts the energy spectrum of neutron-induced captured gamma-rays into relative elemental yields and subsequently transforms them to dry-weight elemental fractions. However, NCS logs are not usually included in a well-logging suite due to cost, tool availability, and borehole conditions. Conventional well logs are typically acquired as a minimum logging program because they provide geologists and petrophysicists with the basic elements for tops identification, stratigraphic correlation, and net-pay determination. Most petrophysical interpretation techniques commonly used to quantify mineral composition from conventional well logs are based on the assumption that lithology is dominated by one or two minerals. In organic shale formations, these techniques are ineffective because all well logs are affected by large variations of mineralogy and pore structure. Even though it is difficult to separate the contribution from each mineral and fluid component on well logs using conventional interpretation methods, well logs still bear essential petrophysical properties that can be estimated using an inversion method. This thesis introduces an inversion-based workflow to estimate mineral and fluid concentrations of shale gas formations using conventional well logs. The workflow starts with the construction and calibration of a mineral model based on core analysis of crushed samples and X-Ray Diffraction (XRD). We implement a mineral grouping approach that reduces the number of unknowns to be estimated by the inversion without loss of accuracy in the representation of the main minerals. The second step examines various methods that can provide good initial values for the inversion. For example, a reliable prediction of kerogen concentration can be obtained using the ΔlogR method (Passey et al., 1990) as well as an empirical correlation with gamma-ray or uranium logs. After the mineral model is constructed and a set of initial values are established, nonlinear joint inversion estimates mineral and fluid concentrations from conventional well logs. An iterative refinement of the mineral model can be necessary depending on formation complexity and data quality. The final step of the workflow is to perform rock classification to identify favorable production zones. These zones are selected based on their hydrocarbon potential inferred from inverted petrophysical properties. Two synthetic examples with known mineral compositions and petrophysical properties are described to illustrate the application of inversion. The impact of shoulder-bed effects on inverted properties is examined for the two inversion modes: depth-by-depth and layer-by-layer. This thesis also documents several case studies from Haynesville and Barnett shales where the proposed workflow was successfully implemented and is in good agreement with core measurements and NCS logs. The field examples confirm the accuracy and reliability of nonlinear inversion to estimate porosity, water saturation, kerogen concentration, and mineral composition.
text

Undulating well trajectories are often drilled to improve length exposure to rock formations, target desirable hydrocarbon-saturated zones, and enhance resolution of borehole measurements. Despite these merits, undulating wells can introduce adverse conditions to the interpretation of borehole measurements which are seldom observed in vertical wells penetrating horizontal layers. Common examples are polarization horns observed across formation bed boundaries in borehole resistivity measurements acquired in highly-deviated wells. Consequently, conventional interpretation practices developed for vertical wells can yield inaccurate results in HA/HZ wells. A reliable approach to account for well trajectory and bed-boundary effects in the petrophysical interpretation of well logs is the application of forward and inverse modeling techniques because of their explicit use of measurement response functions. The main objective of this dissertation is to develop inversion-based petrophysical interpretation methods that quantitatively integrate logging-while-drilling (LWD) multi-sector nuclear (i.e., density, neutron porosity, photoelectric factor, natural gamma ray) and multi-array propagation resistivity measurements. Under the assumption of a multi-layer formation model, the inversion approach estimates formation properties specific to a given measurement domain by numerically reproducing the available measurements. Subsequently, compositional multi-mineral analysis of inverted layer-by-layer properties is implemented for volumetric estimation of rock and fluid constituents. The most important prerequisite for efficient petrophysical inversion is fast and accurate forward models that incorporate specific measurement response functions for numerical simulation of LWD measurements. In the nuclear measurement domain, first-order perturbation theory and flux sensitivity functions (FSFs) are reliable and accurate for rapid numerical simulation. Albeit efficient, these first-order approximations can be inaccurate when modeling neutron porosity logs, especially in the presence of borehole environmental effects (tool standoff or/and invasion) and across highly contrasting beds and complex formation geometries. Accordingly, a secondary thrust of this dissertation is the introduction of two new methods for improving the accuracy of rapid numerical simulation of LWD neutron porosity measurements. The two methods include: (1) a neutron-density petrophysical parameterization approach for describing formation macroscopic cross section, and (2) a one-group neutron diffusion flux-difference method for estimating perturbed spatial neutron porosity fluxes. Both methods are validated with full Monte Carlo (MC) calculations of spatial neutron detector FSFs and subsequent simulations of neutron porosity logs in the presence of LWD azimuthal standoff, invasion, and highly dipping beds. Analysis of field and synthetic verification examples with the combined resistivity-nuclear inversion method confirms that inversion-based estimation of hydrocarbon pore volume in HA/HZ wells is more accurate than conventional well-log analysis. Estimated hydrocarbon pore volume from conventional analysis can give rise to errors as high as 15% in undulating HA/HZ intervals.
text

Pulsed-neutron borehole measurements involve a physical process in which a source emits energetic neutrons that lose energy upon collisions with formation nuclei, and are eventually captured by a nucleus to form a heavier, excited state. The excited nucleus decays to its ground state by the emission of gamma rays. Both thermal-neutron and gamma-ray populations decay with time at a rate defined by Sigma, which is a nuclear property that quantifies a material’s ability to capture thermal neutrons. The large contrast in Sigma between hydrocarbon and salty connate water enables calculations of water saturation directly from pulsed-neutron measurements. Sigma logs have proven useful in the assessment of thinly bedded formations because they exhibit a small volume of investigation, and have been deemed superior to resistivity logs in the petrophysical evaluation of carbonate formations. The recognized potential of Sigma logs in formation evaluation initiated the development of multi-detector Logging-While-Drilling (LWD) Sigma measurements. These measurements are acquired using one thermal-neutron and two gamma-ray detectors at different spacings from the source. Such a design is aimed at providing distinct radial depths of investigation to detect filtrate invasion in the near-wellbore zone. Despite their formation-evaluation potential, multi-detector time-decay measurements commonly remain affected by invasion, shoulder-bed, and well-deviation effects. The purpose of this dissertation is to develop a fast-forward simulation method to reproduce multi-detector time decays and combine the method with inversion techniques to improve the petrophysical interpretation of LWD Sigma measurements. First-order perturbation theory and a library of pre-calculated Monte Carlo detector-specific sensitivity functions and time decays are used to numerically simulate borehole Sigma measurements in realistic logging environments. The new simulation method is one hundred thousand times faster than rigorous Monte Carlo calculations and remains within two capture units of disparity. Next, the fast-forward simulation method is embedded within inversion algorithms to estimate layer-by-layer radial length of invasion and formation Sigma corrected for shallow invasion, shoulder-bed, and well-deviation effects. Both fast-forward and inverse modeling algorithms are benchmarked against laboratory and synthetic time decays. The improvement of formation Sigma obtained with inversion-based interpretation leads to an improvement in the estimation of Sigma-derived water saturation. Likewise, the estimated radial length of invasion is combined with neutron and density measurements to correct the latter for invasion effects. Results indicate that the inversion-based interpretation method is well suited for the evaluation of high-porosity formations invaded by salty mud filtrate. Inversion-based interpretation of field LWD time decays enables the estimation of lower values of water saturation when compared to conventional Sigma interpretation or resistivity methods. Estimated values of water saturation are as much as fifty percent lower than predicted by conventional interpretation of Sigma logs in the case of measurements affected by shoulder-bed effects, and as much as one hundred percent lower than predicted by the conventional interpretation method for measurements additionally affected by salty filtrate invasion. The key attributes of the combined petrophysical interpretation of multi-detector Sigma, neutron, and density measurements developed in this dissertation are that it explicitly enforces the physics of all nuclear measurements, honors the pressure and temperature dependency of reservoir fluid nuclear properties, and takes into account a-priori information such as mud-filtrate salinity, connate-water salinity, and bed-boundary locations.
text

A petrophysical study of the upper grainstone/packstone reservoir of the Oxfordian Smackover Formation in Little Cedar Creek Field was conducted, integrating core description, thin section analysis, log interpretation and cathodoluminescense to characterize controls on oil production in the upper reservoir. Little Cedar Creek Field produces approximately 2.4 million barrels (bbls) of oil annually and is currently in secondary recovery. By analyzing petrophysical characteristics such as porosity and pore type and correlating them to facies changes, better predictions can be made to optimize secondary recovery. The diagenetic history of the ooid-peloid grainstone records six separate events. Early marine phreatic dogtooth sparry rim cement helped create the framework that allowed it to maintain a good portion of its depositional porosity as it underwent subsequent compaction, dissolution and cementation events. The most common porosity types are vuggy, oomoldic and intercrystalline. The Smackover Formation ooid-peloid grainstone/packstone unit consists of multiple alternating ooid-peloid grainstone and peloid packstone/wackestone facies with varying porosity types. The most common types are oomoldic and vuggy with a range of preserved intergranular porosity. Porosity in the grainstone facies averages 17% and 5.6% in the packstone/wackestone facies. The number of facies changes within the upper reservoir does not play a significant role in controlling well production. Facies changes are too thin to be identifiable utilizing well logs alone, although neutron and density well logs do trace a close relationship between log values and core plug analysis values of porosity. Core reports indicate that porosity and permeability correlate strongly with pore size and facies. Areas with thicker accumulations of grainstone facies have higher porosity and permeability values and have higher oil production. Isopach maps of the cumulative grainstone facies indicate thick build-ups parallel to strike for the formation, consistent with a shoal environment. The strongest predictor of well production is the cumulative thickness of grainstone facies within the grainstone/packstone unit of the Smackover Formation. The grainstone is thickest in the southwest part of the field and pinches out updip in the northwest. Secondary recovery gas injection would be most effective if applied in the southwestern portion of the field because it could effectively sweep the oil updip towards the stratigraphic trap.

Rock type is a key concept in modern reservoir characterization that straddles multiple scales and bridges multiple disciplines. Reservoir rock classification (or simply rock typing) has been recognized as one of the most effective description tools to facilitate large-scale reservoir modeling and simulation. This dissertation aims to integrate core data and well logs to enhance reservoir description by classifying reservoir rocks in a geologically and petrophysically consistent manner. The main objective is to develop scientific approaches for utilizing multi-physics rock data at different time and length scales to describe reservoir rock-fluid systems. Emphasis is placed on transferring physical understanding of rock types from limited ground-truthing core data to abundant well logs using fast log simulations in a multi-layered earth model. Bimodal log-normal pore-size distribution functions derived from mercury injection capillary pressure (MICP) data are first introduced to characterize complex pore systems in carbonate and tight-gas sandstone reservoirs. Six pore-system attributes are interpreted and integrated to define petrophysical orthogonality or dissimilarity between two pore systems of bimodal log-normal distributions. A simple three-dimensional (3D) cubic pore network model constrained by nuclear magnetic resonance (NMR) and MICP data is developed to quantify fluid distributions and phase connectivity for predicting saturation-dependent relative permeability during two-phase drainage. There is rich petrophysical information in spatial fluid distributions resulting from vertical fluid flow on a geologic time scale and radial mud-filtrate invasion on a drilling time scale. Log attributes elicited by such fluid distributions are captured to quantify dynamic reservoir petrophysical properties and define reservoir flow capacity. A new rock classification workflow that reconciles reservoir saturation-height behavior and mud-filtrate for more accurate dynamic reservoir modeling is developed and verified in both clastic and carbonate fields. Rock types vary and mix at the sub-foot scale in heterogeneous reservoirs due to depositional control or diagenetic overprints. Conventional well logs are limited in their ability to probe the details of each individual bed or rock type as seen from outcrops or cores. A bottom-up Bayesian rock typing method is developed to efficiently test multiple working hypotheses against well logs to quantify uncertainty of rock types and their associated petrophysical properties in thinly bedded reservoirs. Concomitantly, a top-down reservoir description workflow is implemented to characterize intermixed or hybrid rock classes from flow-unit scale (or seismic scale) down to the pore scale based on a multi-scale orthogonal rock class decomposition approach. Correlations between petrophysical rock types and geological facies in reservoirs originating from deltaic and turbidite depositional systems are investigated in detail. Emphasis is placed on the cause-and-effect relationship between pore geometry and rock geological attributes such as grain size and bed thickness. Well log responses to those geological attributes and associated pore geometries are subjected to numerical log simulations. Sensitivity of various physical logs to petrophysical orthogonality between rock classes is investigated to identify the most diagnostic log attributes for log-based rock typing. Field cases of different reservoir types from various geological settings are used to verify the application of petrophysical rock classification to assist reservoir characterization, including facies interpretation, permeability prediction, saturation-height analysis, dynamic petrophysical modeling, uncertainty quantification, petrophysical upscaling, and production forecasting.
text

Nuclear logging-while-drilling (LWD) measurements acquired in high-angle and horizontal (HA/HZ) wells are influenced by tool, geometrical, and petrophysical effects. Reliable interpretation of petrophysical and geometrical properties from LWD measurements acquired in thinly-bedded formations requires that gamma ray, density, photoelectric (PEF), and neutron measurements be quantitatively integrated with explicit consideration of their effective volume of investigation (EVOI). One of the effects of different tool EVOIs is false gas density-neutron crossovers across thinly-bedded formations. Also, in the presence of tool eccentricity, azimuthally-varying standoff gives rise to an azimuthally-varying effective depth of investigation (EDOI), which introduces errors in the inference of formation dip. Conventional Monte Carlo simulations of nuclear measurements are computationally expensive in reproducing multi-sector LWD responses in HA/HZ wells. Using linear iterative refinement of pre-calculated flux sensitivity functions (FSFs), we introduce a fast method for numerical simulation of LWD nuclear images in the presence of tool eccentricity along any well trajectory. Our investigation of measurement responses from FSFs motivates techniques to explicitly consider the EVOI of LWD nuclear measurements. Simple radial DOI and standoff corrections suffice for interpretation of gamma-gamma images but are inadequate for neutron responses due to larger EVOI and azimuthal aperture. We introduce a new azimuthal deconvolution method of neutron images to improve bed-boundary detection. Neutron DOI varies significantly with porosity, whereby we correct neutron images for penetration length due to changes of porosity along the well trajectory. In addition, we implement a new method of separate linear iterative refinement on neutron thermal group responses to improve the resolution of neutron images across heterogeneous and thinly-bedded formations. The method reduces shoulder-bed effects and false neutron-density gas crossovers. We corroborate these techniques with rigorous Monte Carlo simulations in vertical and deviated wells. A field example of application conclusively indicates that numerical simulation of LWD nuclear measurements is necessary for reliable estimation of petrophysical properties.
text

Diagenetic and structural aspects of tight gas sandstones must be addressed concurrently in order to fully understand low-permeability sandstones and to better predict their reservoir quality attributes that arise from a combination of pore-scale and fracture distribution characteristics. This dissertation focuses on aspects of rock evolution that are germane to concurrent structural and diagenetic evolution, such as loading and thermal history, rock mechanical property evolution, and fracture timing. I tested the hypothesis that the cement precipitation step, governed by thermal exposure and grain surface attributes, governs how sandstone attributes evolve using observations from the Late Cretaceous Williams Fork sandstones from the Piceance Basin, Colorado. My research shows that essential information for predicting and understanding fracture patterns in sandstone can be obtained by unraveling cement precipitation (diagenetic) history. Fractures depend on the mechanical properties existing during fracture growth. I show that key rock mechanical properties (subcritical crack index, Young's modulus and Poisson's ratio), petrophysical behavior, and reservoir quality depend in a systematic way on time-temperature history and the intrinsic grain surface attributes of these sandstones. I classified the Williams Fork lithofacies petrographically and correlated those with log responses to create a model that can be used to predict reservoir quality and diagenesis directly from well logs. I determined rock mechanical characteristics by measuring the subcritical crack index (SCI), a mechanical property that influences fracture distribution characteristics, and by examining log-derived bulk mechanical properties. To quantify the influence of quartz cementation on the SCI and to determine the range of SCI values for sandstone of given framework composition at different diagenetic stages, I measured SCI on Williams Fork core samples and their outcrop equivalents. Diagenetic modeling is applied to determine the sandstone characteristics during fracturing.
text