Investigating CO<Subscript>2</Subscript>-enhanced oil recovery potential for a mature oil field: a case study based on Ankleshwar oil field,Cambay Basin,India

CO₂-enhanced oil recovery (EOR) is an upcoming technology in India. At present, no Indian field is under CO₂-EOR and implementation of this technique to a mature oil field needs a rigorous study. In the present work, we made an attempt to investigate the CO₂-EOR potential of a mature oil field, situated in Cambay Basin, India. The field was put on production in 1961, and it has produced approximately 65.36 MMt oil during massive water flooding, leading to residual oil reserves of 6.49 MMt. The operator of the field is interested in incremental oil recovery from this field by injecting CO₂. This requires estimation of incremental oil recovery potential of the field by carrying out systematic study. We, therefore, developed a conceptual model inspired by Ankleshwar oil field of Cambay Basin using available information provided by the field operator and carried out systematic studies to establish an optimized strategy for CO₂ injection. To achieve this goal, we investigated the effect of various operational parameters on oil recovery efficiency of our conceptual model and selected optimum parameters for reservoir simulations. Simulation results clearly indicate that the field can be a good candidate for CO₂-EOR, and an additional oil recovery of 10.4% of hydrocarbon pore volume is feasible. Major outcome of the study is an optimized black-oil simulation model, which is in good agreement with the fine grid compositional model of high accuracy. The proposed black-oil model can easily be implemented and updated compared with compute intensive finer compositional simulation model.     

Temporal sampling strategies and uncertainty in calibrating a conceptual hydrological model for a small boreal catchment

How much data is needed for calibration of a hydrological catchment model? In this paper we address this question by evaluating the information contained in different subsets of discharge and groundwater time series for multi‐objective calibration of a conceptual hydrological model within the framework of an uncertainty analysis. The study site was a 5·6‐km² catchment within the Forsmark research site in central Sweden along the Baltic coast. Daily time series data were available for discharge and several groundwater wells within the catchment for a continuous 1065‐day period. The hydrological model was a site‐specific modification of the conceptual HBV model. The uncertainty analyses were based on a selective Monte Carlo procedure. Thirteen subsets of the complete time series data were investigated with the idea that these represent realistic intermittent sampling strategies. Data subsets included split‐samples and various combinations of weekly, monthly, and quarterly fixed interval subsets, as well as a 53‐day ‘informed observer’ subset that utilized once per month samples except during March and April—the months containing large and often dominant snow melt events—when sampling was once per week. Several of these subsets, including that of the informed observer, provided very similar constraints on model calibration and parameter identification as the full data record, in terms of credibility bands on simulated time series, posterior parameter distributions, and performance indices calculated to the full dataset. This result suggests that hydrological sampling designs can, at least in some cases, be optimized. Copyright © 2009 John Wiley & Sons, Ltd.     

Radio–optical alignments in a low radio luminosity sample

We present an optically based study of the alignment between the radio axes and the optical major axes of eight z ∼0.7 radio galaxies in a 7C sample. The radio galaxies in this sample are ≈20 times less radio‐luminous than 3C galaxies at the same redshift, and are significantly less radio-luminous than any other well-defined samples studied to date. Using Nordic Optical Telescope images taken in good seeing conditions at rest frame wavelengths just longward of the 4000-Å break, we find a statistically significant alignment effect in the 7C sample. Furthermore, in two cases where the aligned components are well separated from the host we have been able to confirm spectroscopically that they are indeed at the same redshift as the radio galaxy. However, a quantitative analysis of the alignment in this sample and in a corresponding 3C sample from HST archival data indicates that the percentage of aligned flux may be lower and of smaller spatial scale in the 7C sample. Our study suggests that alignments on the 50-kpc scale are probably closely related to the radio luminosity, whereas those on the 15‐kpc scale are not. We discuss these results in the context of popular models for the alignment effect.     

The DNSC08GRA global marine gravity field from double retracked satellite altimetry

Satellite radar altimetry has been monitoring the earth’s oceans from space for several decades. However, only the GEOSAT and ERS-1 geodetic mission data recorded more than a decade ago provide altimetry with adequate spatial coverage to derive a high-resolution marine gravity field. The original geodetic mission data suffer from degradation in quality and coverage close to the coast and in Polar Regions as well as the occasionally wrongly retracking of these, even in the open ocean. In order to improve the quality of these geodetic mission data and to derive a new improved global marine gravity field called DNSC08GRA, a new double retracking technique for analyzing the waveform data has been developed. Multiple retracking allows the system to retrack more data to increase the spatial coverage of the data. Subsequently, a second retracking run is used to enhance the SSH determination by using information from the first fitting to inform the second set of retrackers about smoothly varying sea state parameters. The development of the new global marine gravity field DNSC08GRA is described in this paper. Besides application of new retracking techniques the radar altimetry has been processed using EGM2008 as reference and augmented with ArcGP gravity data and laser altimetry from ICESat to close the Polar gap. DNSC08GRA is seen to perform significantly better than previous global marine gravity field like KMS02. The improvement in accuracy is better than 20% in general, but in coastal regions, the improvement is in many places of the order of 40–50% compared to older global marine gravity field KMS02.     

Bioavailability of PAH in effluent water from an aluminium works evaluated by transplant caging and biliary fluorescence measurements of Atlantic cod (Gadus morhua L.)

Fixed wavelength fluorescence (FF) of bile was used as a biomarker for PAH exposure in fish caged adjacent to an aluminium works. After 30 days of caging, a 20–50 fold increase of pyrene fluorescence was found in groups caged adjacent to the works as compared to two control groups. The PAH uptake occurred mainly from suspended particles rather than from PAH contaminated sediments at the caging sites. Our results demonstrates the combined use of fish caging and bile-fluorescence measurements as a powerful field model for the monitoring of coastal waters subjected to PAH pollution.     

Ocean urea fertilization for carbon credits poses high ecological risks

The proposed plan for enrichment of the Sulu Sea, Philippines, a region of rich marine biodiversity, with thousands of tonnes of urea in order to stimulate algal blooms and sequester carbon is flawed for multiple reasons. Urea is preferentially used as a nitrogen source by some cyanobacteria and dinoflagellates, many of which are neutrally or positively buoyant. Biological pumps to the deep sea are classically leaky, and the inefficient burial of new biomass makes the estimation of a net loss of carbon from the atmosphere questionable at best. The potential for growth of toxic dinoflagellates is also high, as many grow well on urea and some even increase their toxicity when grown on urea. Many toxic dinoflagellates form cysts which can settle to the sediment and germinate in subsequent years, forming new blooms even without further fertilization. If large-scale blooms do occur, it is likely that they will contribute to hypoxia in the bottom waters upon decomposition. Lastly, urea production requires fossil fuel usage, further limiting the potential for net carbon sequestration. The environmental and economic impacts are potentially great and need to be rigorously assessed.     

Architecture of growth basins in a tidally influenced,prodelta to delta-front setting: The Triassic succession of Kvalpynten,East Svalbard

World-class examples of fault-controlled growth basins with associated syn-kinematic sedimentary fill are developed in Upper Triassic prodelta to delta-front deposits exposed at Kvalpynten, SW Edgeøya in East Svalbard. They are interpreted to have interacted with north-westerly progradation of a regional delta system. The syn-kinematic successions consist of 4 to 5 coarsening-upward units spanning from offshore mudstones to subtidal heterolithic bars and compound tidal dunes, which were blanketed by regional, post-kinematic sandstone sheets deposited as laterally continuous, subaqueous tidal dune fields. The rate of growth faulting is reflected in the distribution of accommodation, which governs sedimentary architecture and stacking patterns within the coarsening-upward units. Fully compartmentalized basins (12, 200–800 m wide and c. 150 m high grabens and half grabens) are characterized by syn-kinematic sedimentary infill. These grabens and half-grabens are separated by 60–150 m high horsts composed of pro-delta to distal delta-front mudstones. Grabens host tabular tidal dunes (sandwaves), whereas half-grabens bound by listric faults (mainly south-dipping) consist of wedge-shaped, rotated strata with erosive boundaries proximal to the uplifted fault block crests. Heterolithic tidal bars (sand ridges) occur in narrow half-grabens, showing migration oblique to the faults, up the dipslope. Structureless sandstone wedges and localized subaqueous slumps that formed in response to collapse of the block crests were only documented in half-grabens. Late-kinematic deposition during the final stages of faulting occurred in partly compartmentalized basins, filled with variably thick sets of continuous sandstone belts (compound tidal dunes).     

U–Pb TIMS age constraints on the evolution of the Neoproterozoic Meatiq Gneiss Dome,Eastern Desert,Egypt

Ages are used to constrain the temporal evolution of the Meatiq Gneiss Dome, Eastern Desert, Egypt, by dating (ID-TIMS) pre-, syn-, and post-tectonic igneous rocks in and around the dome. The Um Ba’anib Orthogneiss, comprising the deepest exposed structural levels of the dome, has a crystallization age of 630.8 ± 2 Ma. The overlying mylonites are interpreted to be a thrust sheet/complex (Abu Fannani Thrust Sheet) of highly mylonitized metasediments (?), migmatitic amphibolites, and orthogneisses with large and small tectonic lenses of less-deformed intrusives. Two syn-tectonic diorite lenses in this complex have crystallization ages of 609.0 ± 1.0 and 605.8 ± 0.9 Ma, respectively. The syn-tectonic Abu Ziran diorite, cutting across the tectonic contact between mylonite gneisses of the Abu Fannani Thrust Sheet and a structurally overlying thrust sheet of eugeoclinal rocks (“Pan-African nappe”), has a magmatic emplacement age of 606.4 ± 1.0 Ma. Zircons from a gabbro (Fawakhir ophiolite) within the eugeoclinal thrust sheet yielded a crystallization age of 736.5 ± 1.2 Ma. The post-tectonic Fawakhir monzodiorite intrudes the ophiolitic rocks and has an emplacement age of 597.8 ± 2.9 Ma. Two other post-tectonic granites, the Arieki granite that intrudes the foliated Um Ba’anib Orthogneiss, and the Um Had granite that cuts the deformed Hammamat sediments, have emplacement ages of 590 ± 3.1 and 596.3 ± 1.7 Ma, respectively. We consider formation of the Meatiq Gneiss Dome to be a young structural feature (<631 Ma), and our preferred tectonic interpretation is that it formed as a result of NE–SW shortening contemporaneous with folding of the nearby Hammamat sediments around 605–600 Ma, during oblique collision of East and West Gondwana.     

Caprock interaction with CO2: A laboratory study of reactivity of shale with supercritical CO2 and brine

Crushed rock from two caprock samples, a carbonate-rich shale and a clay-rich shale, were reacted with a mixture of brine and supercritical CO₂ (CO₂–brine) in a laboratory batch reactor, at different temperature and pressure conditions. The samples were cored from a proposed underground CO₂ storage site near the town of Longyearbyen in Svalbard. The reacting fluid was a mixture of 1 M NaCl solution and CO₂ (110 bar) and the water/rock ratio was 20:1. Carbon dioxide was injected into the reactors after the solution had been bubbled with N₂, in order to mimic O₂-depleted natural storage conditions. A control reaction was also run on the clay-rich shale sample, where the crushed rock was reacted with brine (CO₂-free brine) at the same experimental conditions. A total of 8 batch reaction experiments were run at temperatures ranging from 80 to 250 °C and total pressures of 110 bar (∼40 bar for the control experiment). The experiments lasted 1–5 weeks.Fluid analysis showed that the aqueous concentration of major elements (i.e. Ca, Mg, Fe, K, Al) and SiO₂ increased in all experiments. Release rates of Fe and SiO₂ were more pronounced in solutions reacted with CO₂–brine as compared to those reacted with CO₂-free brine. For samples reacted with the CO₂–brine, lower temperature reactions (80 °C) released much more Fe and SiO₂ than higher temperature reactions (150–250 °C). Analysis by SEM and XRD of reacted solids also revealed changes in mineralogical compositions. The carbonate-rich shale was more reactive at 250 °C, as revealed by the dissolution of plagioclase and clay minerals (illite and chlorite), dissolution and re-precipitation of carbonates, and the formation of smectite. Carbon dioxide was also permanently sequestered as calcite in the same sample. The clay-rich shale reacted with CO₂–brine did not show major mineralogical alteration. However, a significant amount of analcime was formed in the clay-rich shale reacted with CO₂-free brine; while no trace of analcime was observed in either of the samples reacted with CO₂–brine.     

Using the pseudospectral method on curved grids for 2D elastic forward modelling1

When applying the conventional Fourier pseudospectral method (FSM) on a Cartesian grid that has a sufficient size to propagate a pulse, spurious diffractions from the staircase representation of the curved interfaces appear in the wavefield. It is demonstrated that these non-physical diffractions can be eliminated by using curved grids that conform to all the interfaces of the subsurface. Methods for solving the 2D acoustic wave equation using such curved grids have been published previously by the authors. Here the extensions to the full 2D elastic wave equations are presented. The curved grids are generated by using the so-called multiblock strategy which is a well-known concept in computational fluid dynamics. In principle the sub-surface is divided into a number of contiguous subdomains. A separate grid is generated for each subdomain patching the grid lines across domain boundaries to obtain a globally continuous grid. Using this approach, even configurations with pinch outs can be handled. The curved grid is taken to constitute a generalized curvilinear coordinate system. Thus, the elastic equations have to be written in a curvilinear frame before applying the numerical scheme. The method implies that twice the number of spatial derivatives have to be evaluated compared to the conventional FSM on a Cartesian grid. However, it is demonstrated that the extra terms are more than compensated for by the fewer grid points needed in the curved approach.