Estimate CO2 storage capacity of the Johansen formation: numerical investigations beyond the benchmarking exercise

In this paper, Shell’s in-house reservoir simulator MoReS is applied to a recently introduced CO₂ sequestration benchmark problem entitled “Estimation of the CO₂ Storage Capacity of a Geological Formation” (Class et al. 2008). The principal objective of this benchmark is the simulation of CO₂ distribution within a modeling region, and leakage of CO₂ outside of it, for a period of 50 years. This study goes beyond the benchmarking exercise to investigate additional factors with direct relevance to CO₂ storage capacity estimations: water and gas relative permeabilities, permeability anisotropy, presence of sub-seismic features (conductive fractures, thin shale layers), regional hydrodynamic gradient, CO₂-enriched brine convection (due to brine density differences), and injection rates. The effects of hydrodynamic gradients and gravitationally induced convection only become significant over 100 s of years. This study has thus extended simulation time to 1,000 years. It is shown that grid resolution significantly impacts results. Vertical-grid refinement results in larger and thinner CO₂ plumes. Lateral-grid refinement delays leakage out of the model domain and reduces injection pressure for a given injection rate. Sub-seismic geological features such as fractures/faults and shale layers are demonstrated to have impact on CO₂ sequestration. Fractures located up-dip from the injector may lead to more leakage while the opposite may happen in the presence of fractures perpendicular to the dip. Thin shale layers produce stacked CO₂ blankets. They should be explicitly represented instead of being upscaled using a reduced vertical to horizontal permeability ratio. Results are seen to be far more sensitive to gas relative permeability and hysteresis than to variations in the water relative permeability models used. For a multi-injectors project, there is scope to optimize the phasing of injections to avoid potential fracturing near injectors.     

The surface energy balance at the Huygens landing site and the moist surface conditions on Titan

The Huygens Probe provided a wealth of data concerning the atmosphere of Titan. It also provided tantalizing evidence of a small amount of surface liquid. We have developed a detailed surface energy balance for the Probe landing site. We find that the daily averaged non-radiative fluxes at the surface are 0.7 W m^?2, much larger than the global average value predicted by McKay et al. (1991) of 0.037 W m^?2. Considering the moist surface, the methane and ethane detected by the Probe from the surface is consistent with a ternary liquid of ethane, methane, and nitrogen present on the surface with mole fractions of methane, ethane, and nitrogen of 0.44, 0.34, and 0.22, respectively, and a total mass load of ～0.05 kg m^?2. If this liquid is included in the surface energy balance, only a small fraction of the non-radiative energy is due to latent heat release (～10^?3 W m^?2). If the amount of atmospheric ethane is less than 0.6×10^?5, the surface liquid is most likely evaporating over timescales of 5 Titan days, and the moist surface is probably a remnant of a recent precipitation event. If the surface liquid mass loading is increased to 0.5 kg m^?2, then the liquid lifetime increases to ～56 Titan days. Our modeling results indicate a dew cycle is unlikely, given that even when the diurnal variation of liquid is in equilibrium, the diurnal mass variation is only 3% of the total liquid. If we assume a high atmospheric mixing ratio of ethane (>0.6×10^?5), the precipitation of liquid is large (38 cm/Titan year for an ethane mixing ratio of 2×10^?5). Such a flux is many orders of magnitude in excess of the photochemical production rate of ethane.     

Evaluation of image analysis methods used for quantification of particle angularity

Angularity is an important parameter in the characterization of particle morphology that is used to interpret the transport history of particles in sedimentary deposits. In the past, visual classification using silhouette charts was widely used to determine particle angularity, but this approach is subjective and time‐consuming. With advances in modern image analysis techniques and low‐cost software packages, it is possible to rapidly quantify particle angularity more objectively than using visual classification methods. This study re‐examines the performance of three existing image analysis methods and one new image analysis procedure, applied to six rock and sediment samples that were visually different in angularity. To facilitate comparison between the angularity results, measurements were reduced to rankings for each aggregate sample. These results show that the four image analysis methods rank the angularity of the samples differently, and that none rank the mean angularity index in the same order as the angularity ranking using visual classification. Therefore, further research is needed to develop an image analysis method that can quantify the angularity of sedimentary particles more precisely.     

Late glacial palaeoceanography of Hinlopen Strait, northern Svalbard

Timing and structure of the Late and post-glacial development of the northern Svalbard margin, together with the initial influx of Atlantic water into the Arctic Ocean are still very poorly constrained. We investigated a sediment core (NP94-51) from a high accumulation area on the continental shelf north of Hinlopen Strait with the purpose of resolving the timing and structure of the last deglaciation. Detailed analyses of ice-rafted detritus, benthic and planktonic foraminiferal fauna, diatom flora, grain size and radiocarbon dates are used to reconstruct the palaeoceanographic evolution of the area. Our results indicate that the disintegration of Hinlopen Strait ice and possibly the northern margin of the Svalbard Ice Sheet commenced between 13.7 and 13.9 ¹⁴C Ky BP. Influx of subsurface Atlantic waters into the area (12.6 ¹⁴C Ky BP) and the retreat of the sea ice cover, with the accompanying opening of the surface waters (10.8 ¹⁴C Ky BP), happened at different times and both much later than the disintegration of the ice sheets. The transition into the Holocene shows a two-step warming.     

Macroscopic diffusivity in concrete determined by computational homogenization

Effective moisture and chloride ion diffusivity coefficients for concrete are determined by computational homogenization, where concrete is modeled on the mesoscale as a heterogenous three‐phase composite material. By imposing moisture and chloride ion gradients on a representative volume element, effective macroscale properties are obtained through finite element analysis. A parametric study of the effects of the ballast content was carried out. The numerical results correspond well with an estimate of the Hashin–Shtrikman type available in the literature. The computational homogenization strategy proposed here also includes the interfacial transition zone, and its influence on the effective diffusivity coefficients is assessed. Copyright © 2012 John Wiley & Sons, Ltd.     

A Holocene sea‐level curve and revised isobase map based on isolation basins from near the southern tip of Norway

We report a new Holocene relative sea‐level curve based on the stratigraphy in five closely located isolation basins near Lista in southernmost Norway. The results detail the progress and timing of the mid‐Holocene Tapes transgression, the peak of which in this region represents the highest postglacial sea level, as well as the rate of land emergence since then. One additional cored basin is situated above the marine limit. All the basins have bedrock sills that were levelled using a differential GPS. Isolation and ingression boundaries were identified by macrofossil analysis and radiocarbon dated on terrestrial plant remains. In most cases several dates were obtained from each transition. Relative sea level rose with a mean rate of 7 mm a^?1 during the last part of the Tapes transgression 8600?8200 cal. a BP and then gradually slowed to a mean rate of 1 mm a^?1 from 8200?7000 cal. a BP. Mean sea level reached ～5 m higher than the present level when the transgression culminated. Land emergence took place after this, first slowly at a mean rate of 0.4 mm a^?1 until ～3900 cal. a BP before it increased to 2.6 mm between 3900 and 3400 cal. a BP. Since then it has slowly decreased until today and has been ～0.2 mm a^?1 for the last 2000 years. Based on the new curve we present updated Tapes isobases for the region that are displaced by ～20 km in relation to the existing model. From one basin we also report a 5–10 cm thick layer of sorted, sandy gravel, embedded in a more than 5‐m‐thick deposit of homogeneous shallow‐marine mud. The gravel was deposited ～5500 cal. a BP, which is the same age as a tsunami deposit previously mapped in Shetland. As several typical characteristics of tsunami facies deposits are lacking, the origin of the gravel layer remains inconclusive.     

Improved confidence in climate change projections of precipitation further evaluated using daily statistics from ENSEMBLES models

Probability density functions for daily precipitation data are used as a validation tool comparing station measurements to seven transient regional climate model runs, with a horizontal resolution of 25 km and driven by the SRES A1B scenario forcing, within the ENSEMBLES project. The validation is performed for the control period 1961–1990 for eight predefined European subregions, and a ninth region enclosing all eight subregions, with different climate characteristics. Models that best match the observations are then used for making climate change projections of precipitation distributions during the twenty-first century for each subregion separately. We find, compared to the control period, a distinct decrease in the contribution to the total precipitation for days with moderate precipitation and a distinct increase for days with more intense precipitation. This change in contribution to the total precipitation is found to amplify with time during all of the twenty-first century with an average rate of 1.1% K^?1. Furthermore, the crossover point separating the decreasing from the increasing contributions does not show any significant change with time for any specific subregion. These results are a confirmation and a specification of the results from a previous study using the same station measurements but with a regional climate model ensemble within the PRUDENCE project.     

Late Archean basement in the Bangenhuken Complex of the Nordbreen Nappe, western Ny-Friesland, Svalbard

The rocks of western Ny-Friesland, northern Svalbard, are part of a tectonostratigraphy including four thrust sheets, each composed mainly of orthogneisses overlain by younger metasedimentary rocks. Previous geochronological studies have shown that the orthogneisses are dominated by ca. 1750 Mya granitoids. This study of a quartz-monzonite in one of the thrust sheets, the Nordbreen Nappe, yields a single-zircon U-Pb ion-microprobe age of 2709 Ø 28 My. This is the oldest rock unit so far reported in the Svalbard Caledonides. However, age-determinations on detrital zircons in the metasediments of western Ny-Friesland have shown that Late Archean rocks were prominent sources. The new ages presented here provide the first evidence of a local source for these sedimentary rocks.