Deposits of the sandy braided South Saskatchewan River: Implications for the use of modern analogs in reconstructing channel dimensions in reservoir characterization

Estimation of the dimensions of fluvial geobodies from core data is a notoriously difficult problem in reservoir modeling. To try and improve such estimates and, hence, reduce uncertainty in geomodels, data on dunes, unit bars, cross-bar channels, and compound bars and their associated deposits are presented herein from the sand-bed braided South Saskatchewan River, Canada. These data are used to test models that relate the scale of the formative bed forms to the dimensions of the preserved deposits and, therefore, provide an insight as to how such deposits may be preserved over geologic time. The preservation of bed-form geometry is quantified by comparing the Alluvial architecture above and below the maximum erosion depth of the modem channel deposits. This comparison shows that there is no significant difference in the mean set thickness of dune cross-strata above and below the basal erosion surface of the contemporary channel, thus suggesting that dimensional relationships between dune deposits and the formative bed-form dimensions are likely to be valid from both recent and older deposits. The data show that estimates of mean bankfull flow depth derived from dune, unit bar, and cross-bar channel deposits are all very similar. Thus, the use of all these metrics together can provide a useful check that all components and scales of the alluvial architecture have been identified correctly when building reservoir models. The data also highlight several practical issues with identifying and applying data relating to cross-strata. For example, the deposits of unit bars were found to be severely truncated in length and width, with only approximately 10% of the mean bar-form length remaining, and thus making identification in section difficult. For similar reasons, the deposits of compound bars were found to be especially difficult to recognize, and hence, estimates of channel depth based on this method may be problematic. Where only core data are available (i.e., no outcrop data exist), formative flow depths are suggested to be best reconstructed using cross-strata formed by dunes. However, theoretical relationships between the distribution of set thicknesses and formative dune height are found to result in slight overestimates of the latter and, hence, mean bankfull flow depths derived from these measurements. This article illustrates that the preservation of fluvial cross-strata and, thus, the paleohydraulic inferences that can be drawn from them, are a function of the ratio of the size and migration rate of bed forms and the time scale of aggradation and channel migration. These factors must thus be considered when deciding on appropriate length:thickness ratios for the purposes of object-based modeling in reservoir characterization.

Reservoir characterization of the Cardium Formation in the Ferrier Oilfield, west-central Alberta, Canada.

This thesis has the main aim of defining the lithostratigraphy, depositional architecture, post-depositional modifications and reservoir characteristics of the Cardium Formation in the Ferrier Oilfield, and how these characteristics can have great impact over production rates, GOR and produced fluid discrimination. In the Ferrier area, the Cardium Formation is composed by a NE prograding clastic sequence made up of offshore to shoreface deposits sealed by marine shales. The main reservoir is composed by sandstones and conglomerates interpreted to have deposited in a shoreface depositional environment. Lithofacies and net reservoir thickness mapping led to more detailed understanding of the 3D reservoir architecture, and cross-sections shed light on the Cardium depositional architecture and post-deposition sediment erosion in the Ferrier area. Detailed core logging, thin section, SEM and CL analyses were used to study the mineralogy, texture and pore characterization of the Cardium reservoir, and three main compartments have been identified based on production data and reservoir characteristics. Finally, two situations showing odd production behaviour of the Cardium were resolved. This shed light on the effect of structural features and reservoir quality and thickness over hydrocarbon migration pathways. The Ferrier example offers a unique case of fluid discrimination in clastic reservoirs due both to depositional and post-depositional factors, and could be used as analogue for similar situations in the Western Canadian Sedimentary Basin.

Simulated-annealing-based conditional simulation for the local-scale characterization of heterogeneous aquifers

Simulated-annealing-based conditional simulations provide a flexible means of quantitatively integrating diverse types of subsurface data. Although such techniques are being increasingly used in hydrocarbon reservoir characterization studies, their potential in environmental, engineering and hydrological investigations is still largely unexploited. Here, we introduce a novel simulated annealing (SA) algorithm geared towards the integration of high-resolution geophysical and hydrological data which, compared to more conventional approaches, provides significant advancements in the way that large-scale structural information in the geophysical data is accounted for. Model perturbations in the annealing procedure are made by drawing from a probability distribution for the target parameter conditioned to the geophysical data. This is the only place where geophysical information is utilized in our algorithm, which is in marked contrast to other approaches where model perturbations are made through the swapping of values in the simulation grid and agreement with soft data is enforced through a correlation coefficient constraint. Another major feature of our algorithm is the way in which available geostatistical information is utilized. Instead of constraining realizations to match a parametric target covariance model over a wide range of spatial lags, we constrain the realizations only at smaller lags where the available geophysical data cannot provide enough information. Thus we allow the larger-scale subsurface features resolved by the geophysical data to have much more due control on the output realizations. Further, since the only component of the SA objective function required in our approach is a covariance constraint at small lags, our method has improved convergence and computational efficiency over more traditional methods. Here, we present the results of applying our algorithm to the integration of porosity log and tomographic crosshole georadar data to generate stochastic realizations of the local-scale porosity structure. Our procedure is first tested on a synthetic data set, and then applied to data collected at the Boise Hydrogeophysical Research Site.

Characterization of geothermal reservoirs’ parameters by inverse problem resolution and geostatistical simulations

BTES (borehole thermal energy storage)systems exchange thermal energy by conduction with the surrounding ground through borehole materials. The spatial variability of the geological properties and the space-time variability of hydrogeological conditions affect the real power rate of heat exchangers and, consequently, the amount of energy extracted from / injected into the ground. For this reason, it is not an easy task to identify the underground thermal properties to use when designing. At the current state of technology, Thermal Response Test (TRT) is the in situ test for the characterization of ground thermal properties with the higher degree of accuracy, but it doesn’t fully solve the problem of characterizing the thermal properties of a shallow geothermal reservoir, simply because it characterizes only the neighborhood of the heat exchanger at hand and only for the test duration. Different analytical and numerical models exist for the characterization of shallow geothermal reservoir, but they are still inadequate and not exhaustive: more sophisticated models must be taken into account and a geostatistical approach is needed to tackle natural variability and estimates uncertainty. The approach adopted for reservoir characterization is the “inverse problem”, typical of oil&gas field analysis. Similarly, we create different realizations of thermal properties by direct sequential simulation and we find the best one fitting real production data (fluid temperature along time). The software used to develop heat production simulation is FEFLOW 5.4 (Finite Element subsurface FLOW system). A geostatistical reservoir model has been set up based on literature thermal properties data and spatial variability hypotheses, and a real TRT has been tested. Then we analyzed and used as well two other codes (SA-Geotherm and FV-Geotherm) which are two implementation of the same numerical model of FEFLOW (Al-Khoury model).

AMPLITUDE-VERSUS-ANGLE ANALYSIS AND WIDE-ANGLE-INVERSION OF CROSSWELL SEISMIC DATA IN A CARBONATE RESERVOIR

Crosswell data set contains a range of angles limited only by the geometry of the source and receiver configuration, the separation of the boreholes and the depth to the target. However, the wide angles reflections present in crosswell imaging result in amplitude-versus-angle (AVA) features not usually observed in surface data. These features include reflections from angles that are near critical and beyond critical for many of the interfaces; some of these reflections are visible only for a small range of angles, presumably near their critical angle. High-resolution crosswell seismic surveys were conducted over a Silurian (Niagaran) reef at two fields in northern Michigan, Springdale and Coldspring. The Springdale wells extended to much greater depths than the reef, and imaging was conducted from above and from beneath the reef. Combining the results from images obtained from above with those from beneath provides additional information, by exhibiting ranges of angles that are different for the two images, especially for reflectors at shallow depths, and second, by providing additional constraints on the solutions for Zoeppritz equations. Inversion of seismic data for impedance has become a standard part of the workflow for quantitative reservoir characterization. Inversion of crosswell data using either deterministic or geostatistical methods can lead to poor results with phase change beyond the critical angle, however, the simultaneous pre-stack inversion of partial angle stacks may be best conducted with restrictions to angles less than critical. Deterministic inversion is designed to yield only a single model of elastic properties (best-fit), while the geostatistical inversion produces multiple models (realizations) of elastic properties, lithology and reservoir properties. Geostatistical inversion produces results with far more detail than deterministic inversion. The magnitude of difference in details between both types of inversion becomes increasingly pronounced for thinner reservoirs, particularly those beyond the vertical resolution of the seismic. For any interface imaged from above and from beneath, the results AVA characters must result from identical contrasts in elastic properties in the two sets of images, albeit in reverse order. An inversion approach to handle both datasets simultaneously, at pre-critical angles, is demonstrated in this work. The main exploration problem for carbonate reefs is determining the porosity distribution. Images of elastic properties, obtained from deterministic and geostatistical simultaneous inversion of a high-resolution crosswell seismic survey were used to obtain the internal structure and reservoir properties (porosity) of Niagaran Michigan reef. The images obtained are the best of any Niagaran pinnacle reef to date.

Evolution and sedimentology of a channel fill in the sandy braided South Saskatchewan River and its comparison to the deposits of an adjacent compound bar

The depositional stratigraphy of within-channel deposits in sandy braided rivers is dominated by a variety of barforms (both singular `unit' bars and complex `compound' bars), as well as the infill of individual channels (herein termed `channel fills'). The deposits of bars and channel fills define the key components of facies models for braided rivers and their within-channel heterogeneity, knowledge of which is important for reservoir characterization. However, few studies have sought to address the question of whether the deposits of bars and channel fills can be readily differentiated from each other. This paper presents the first quantitative study to achieve this aim, using aerial images of an evolving modern sandy braided river and geophysical imaging of its subsurface deposits. Aerial photographs taken between 2000 and 2004 document the abandonment and fill of a 1 3 km long, 80 m wide anabranch channel in the sandy braided South Saskatchewan River, Canada. Upstream river regulation traps the majority of very fine sediment and there is little clay (<1%) in the bed sediments. Channel abandonment was initiated by a series of unit bars that stalled and progressively blocked the anabranch entrance, together with dune deposition and stacking at the anabranch entrance and exit. Complete channel abandonment and subsequent fill of up to 3 m of sediment took approximately two years. Thirteen kilometres of ground-penetrating radar surveys, coupled with 18 cores, were obtained over the channel fill and an adjacent 750 m long, 400 m wide, compound bar, enabling a quantitative analysis of the channel and bar deposits. Results show that, in terms of grain-size trends, facies proportions and scale of deposits, there are only subtle differences between the channel fill and bar deposits which, therefore, renders them indistinguishable. Thus, it may be inappropriate to assign different geometric and sedimentological attributes to channel fill and bar facies in object-based models of sandy braided river alluvial architecture.

Proposta metodológica para o imageamento digital e modelagem virtual 3d de um bloco de rochas travertinas

In this paper we present the methodological procedures involved in the digital imaging in mesoscale of a block of travertines rock of quaternary age, originating from the city of Acquasanta, located in the Apennines, Italy. This rocky block, called T-Block, was stored in the courtyard of the Laboratório Experimental Petróleo "Kelsen Valente" (LabPetro), of Universidade Estadual de Campinas (UNICAMP), so that from it were performed Scientific studies, mainly for research groups universities and research centers working in brazilian areas of reservoir characterization and 3D digital imaging. The purpose of this work is the development of a Model Solid Digital, from the use of non-invasive techniques of digital 3D imaging of internal and external surfaces of the T-Block. For the imaging of the external surfaces technology has been used LIDAR (Light Detection and Range) and the imaging surface Interior was done using Ground Penetrating Radar (GPR), moreover, profiles were obtained with a Gamma Ray Gamae-spectômetro laptop. The goal of 3D digital imaging involved the identification and parameterization of surface geological and sedimentary facies that could represent heterogeneities depositional mesoscale, based on study of a block rocky with dimensions of approximately 1.60 m x 1.60 m x 2.70 m. The data acquired by means of terrestrial laser scanner made available georeferenced spatial information of the surface of the block (X, Y, Z), and varying the intensity values of the return laser beam and high resolution RGB data (3 mm x 3 mm), total points acquired 28,505,106. This information was used as an aid in the interpretation of radargrams and are ready to be displayed in rooms virtual reality. With the GPR was obtained 15 profiles of 2.3 m and 2 3D grids, each with 24 sections horizontal of 1.3 and 14 m vertical sections of 2.3 m, both the Antenna 900 MHz to about 2600 MHz antenna. Finally, the use of GPR associated with Laser Scanner enabled the identification and 3D mapping of 3 different radarfácies which were correlated with three sedimentary facies as had been defined at the outset. The 6 profiles showed gamma a low amplitude variation in the values of radioactivity. This is likely due to the fact of the sedimentary layers profiled have the same mineralogical composition, being composed by carbonate sediments, with no clay in siliciclastic pellitic layers or other mineral carrier elements radioactive

Caracterização de reservatórios siliciclásticos neoaptianos: um estudo do membro carmópolis no campo de Camorim, sub-bacia de Sergipe, Brasil

The Camorim Oilfield, discovered in 1970 in the shallow water domain of the Sergipe Sub-basin, produces hydrocarbons from the Carmópolis Member of the Muribeca Formation, the main reservoir interval, interpreted as siliciclastics deposited in an alluvial-fluvial-deltaic context during a late rifting phase of Neoaptian age, in the Sergipe-Alagoas Basin. The structural setting of the field defines different production blocks, being associated to the evolution of the Atalaia High during the rift stage and subsequent reactivations, encompassing NE-SW trending major normal faults and NWEW trending secondary faults. The complexity of this field is related to the strong facies variation due to the interaction between continental and coastal depositional environments, coupled with strata juxtaposition along fault blocks. This study aims to geologically characterize its reservoirs, to provide new insights to well drilling locations in order to increase the recovery factor of the field. Facies analysis based on drill cores and geophysical logs and the 3D interpretation of a seismic volume, provide a high resolution stratigraphic analysis approach to be applied in this geodynamic transitional context between the rift and drift evolutionary stages of the basin. The objective was to define spatial and time relations between production zones and the preferential directions of fluid flow, using isochore maps that represent the external geometry of the deposits and facies distribution maps to characterize the internal heterogeneities of these intervals, identified in a 4th order stratigraphic zoning. This work methodology, integrated in a 3D geological modelling process, will help to optimize well drilling and hydrocarbons production. This methodology may be applied in other reservoirs in tectonic and depositional contexts similar to the one observed at Camorim, for example, the oil fields in the Aracaju High, Sergipe Sub-basin, which together represent the largest volume of oil in place in onshore Brazilian basins

High-resolution aquifer analog of fluvial-aeolian sediments of the Guarani aquifer system

The Guarani aquifer system (GAS) represents one of the biggest aquifers in the world and is the most relevant groundwater resource in South America. For the first time, by combining field and laboratory measurements, a high-resolution aquifer analog model of fluvial-aeolian sediments of the GAS in São Paulo State (Brazil) is constructed. Three parallel sections of frontal outcrops, 28 m × 5.8 m, and two parallel sections of lateral outcrops, 7 m × 5.8 m, are recorded during open-pit mining of sandy sediments and describe in detail the three-dimensional distribution of the local lithofacies and hydrofacies. Variations of hydraulic conductivity, K, and porosity, n, are resolved on the centimeter scale, and the most permeable units of the fluvial-aeolian facies association are identified. The constructed aquifer analog model shows moderate hydraulic heterogeneity and a mean K value of 1.36 × 10-4 m/s, which is greater than the reported range of K values for the entire GAS in São Paulo State. The results suggest that the examined sedimentary unit constitutes a relevant portion of the GAS in São Paulo State in the context of groundwater extraction and pollution. Moreover, the constructed aquifer analog is considered an ideal basis for future numerical model experiments, aiming at in-depth understanding of the groundwater flow and contaminant transport patterns at this GAS portion or at comparable fluvial-aeolian facies associations. © 2013 Springer-Verlag Berlin Heidelberg.

Análise de resolução em tomografia interpoços em 3-D

A aplicação de métodos de inversão aos dados de múltiplos levantamentos sísmicos interpoços permite a reconstrução de modelos de vagarosidade em 3-D de alta resolução adequados para monitoramento de processos de recuperação avançada de petróleo e caracterização de reservatórios. Entretanto, a falta de cobertura volumétrica uniforme de raios de levantamentos interpoços exige informação adicional ao sistema tomográfico para obtenção de soluções estáveis. A discretização do modelo em uma malha 3-D com células prismáticas triangulares e a decomposição em valores singulares são utilizadas para avaliar a reconstrução tomográfica em 3-D. O ângulo da projeção de modelos-alvo no subespaço ortogonal ao espaço nulo efetivo da matriz tomográfica é um critério adequado para se otimizar a malha de discretização do modelo interpretativo e a geometria de aquisição dos dados de modo a melhorar o condicionamento da reconstrução tomográfica. Esta abordagem pode ser utilizada durante as iterações lineares para redefinir a malha ou avaliar a necessidade de informação a priori adicional ao sistema tomográfico.

Architectural characterization of a delta-front reservoir analogue combining ground penetrating radar and electrical resistivity tomography: Roda Sandstone (Lower Eocene, Graus-Tremp basin, Spain)

Three-dimensional reconstruction of reservoir analogues can be improved combining data from different geophysical methods. Ground Penetrating Radar (GPR) and Electrical Resistivity Tomography (ERT) data are valuable tools, since they provide subsurface information from internal architecture and facies distribution of sedimentary rock bodies, enabling the upgrading of depositional models and heterogeneity reconstruction. The Lower Eocene Roda Sandstone is a well-known deltaic complex widely studied as a reservoir analogue that displays a series of sandstone wedges with a general NE to SW progradational trend. To provide a better understanding of internal heterogeneity of a 10m-thick progradational delta-front sandstone unit, 3D GPR data were acquired. In addition, common midpoints (CMP) to measure the sandstone subsoil velocity, test profiles with different frequency antennas (25, 50 and 100MHz) and topographic data for subsequent correction in the geophysical data were also obtained. Three ERT profiles were also acquired to further constrain GPR analysis. These geophysical results illustrate the geometry of reservoir analogue heterogeneities both depositional and diagenetic in nature, improving and complementing previous outcrop-derived data. GPR interpretation using radar stratigraphy principles and attributes analysis provided: 1)tridimensional geometry of major stratigraphic surfaces that define four units in the GPR Prism, 2) image the internal architecture of the units and their statistical study of azimuth and dips, useful for a quick determination of paleocurrent directions. These results were used to define the depositional architecture of the progradational sandbody that shows an arrangement in very-high-frequency sequences characterized by clockwise paleocurrent variations and decrease of the sedimentary flow, similar to those observed at a greater scale in the same system. This high-frequency sequential arrangement has been attributed to the autocyclic dynamics of a supply-dominated delta- front where fluvial and tidal currents are in competition. The resistivity models enhanced the viewing of reservoir quality associated with cement distribution caused by depositional and early diagenetic processes related to the development of transgressive and regressive systems tracts in igh-frequency sequences.

Characterization of a microcystin and detection of microcystin synthetase genes from a Brazilian isolate of Nostoc

A nostocalean nitrogen-fixing cyanobacterium isolated from an eutrophic freshwater reservoir located in Piracicaba, Sao Paulo, Brazil, was evaluated for the production of hepatotoxic cyclic heptapeptides, microcystins. Morphologically this new cyanobacterium strain appears closest to Nostoc, however, in the phylogenetic analysis of 165 rRNA gene it falls into a highly stable cluster distantly only related to the typical Nostoc cluster. Extracts of Nostoc sp. CENA88 cultured cells, investigated using ELISA assay, gave positive results and the microcystin profile revealed by ESI-Q-TOF/MS/MS analysis confirmed the production of [Dha(7)]MCYST-YR. Further, Nostoc sp. CENA88 genomic DNA was analyzed by PCR for sequences of mcyD, mcyE and mcyG genes of microcystin synthetase (mcy) cluster. The result revealed the presence of mcyD, mcyE and mcyG genes with similarities to those from mcy of Nostoc sp. strains 152 and IO-102-I and other cyanobacterial genera. The phylogenetic tree based on concatenated McyG, McyD and McyE amino acids clustered the sequences according to cyanobacterial genera, with exception of the Nostoc sp. CENA88 sequence, which was placed in a clade distantly related from other Nostoc strains, as previously observed also in the 165 rRNA phylogenetic analysis. The present study describes for the first time a Brazilian Nostoc microcystin producer and also the occurrence of demethyl MCYST-YR variant in this genus. The sequenced Nostoc genes involved in the microcystin synthesis can contribute to a better understanding of the toxigenicity and evolution of this cyanotoxin. (C) 2009 Elsevier Ltd. All rights reserved.

Serotype and mating type characterization of Cryptococcus neoformans by Multiplex PCR

Cryptococcus neoformans is an encapsulated yeast, etiological agent of cryptococcosis. The species is commonly associated with pigeon droppings and plant materials. The aim of the present work was to verify the presence of the yeast in pigeon droppings, and to identify the isolates obtained in serotypes and mating types (MAT). Ten samples of pigeon droppings were collected in the rural area of the city of Alfenas, Brazil. Samples were inoculated in agar Niger medium for fungal isolation and 22 isolates with characteristics of C. neoformans were obtained. The serotypes and MAT were determined by multiplex PCR using specific primers. Serotypes were also determined by using the Kit Crypto Check. Among the 22 samples evaluated, eight were identified as C. neoformans by classic identification tests. These samples were characterized as serotype A by the Kit Crypto check and as serotype A MAT alpha by the multiplex PCR. The present study reinforces the evidence that pigeon droppings are a reservoir for C. neoformans and confirms the prevalence of C. neoformans var. grubii (Aalpha) among environmental isolates. It also demonstrates that multiplex PCR is an acceptable alternative for serotype analysis because it reduces the costs for each reaction and analyses serotype and MAT simultaneously.

Molecular characterization of environmental Cryptococcus neoformans isolated in Vitoria, ES, Brazil

Cryptococcus neoformans is the major cause of fungal meningitis, a potentially lethal mycosis. Bird excreta can be considered a significant environmental reservoir of this species in urban areas, thirty-three samples of pigeon excreta were collected within the city of Vitoria, Brazil. Cryptococcus neoformans was isolated and identified using standard biochemical assays in ten samples. PCR amplification with primer M13 and orotidine monophosphate pyrophosphorylase (URA5) gene-restriction fragment length polymorphism (RFLP) analysis discerned serotypes and genotypes within this species. All isolates were serotype A (C. neoformans var. grubii) and genotype VNI. The two alternative alleles a and α at the mating type locus were determined by PCR amplification and mating assays performed on V8 medium. All isolates were MAT α mating type but only 50% were able to mate in vitro with the opposite mating type MAT a tester strains (JEC20, KN99a and Bt63). This study adds information on the ecology and molecular characterization of C. neoformans in the Southeast region of Brazil.