Post‐stack stratigraphic inversion workflow applied to carbon dioxide storage: application to the saline aquifer of Sleipner field‡

In the Norwegian North Sea, the Sleipner field produces gas with a high CO₂ content. For environmental reasons, since 1996, more than 11 Mt of this carbon dioxide (CO₂) have been injected in the Utsira Sand saline aquifer located above the hydrocarbon reservoir. A series of seven 3D seismic surveys were recorded to monitor the CO₂ plume evolution. With this case study, time‐lapse seismics have been shown to be successful in mapping the spread of CO₂ over the past decade and to ensure the integrity of the overburden. Stratigraphic inversion of seismic data is currently used in the petroleum industry for quantitative reservoir characterization and enhanced oil recovery. Now it may also be used to evaluate the expansion of a CO₂ plume in an underground reservoir. The aim of this study is to estimate the P‐wave impedances via a Bayesian model‐based stratigraphic inversion. We have focused our study on the 1994 vintage before CO₂ injection and the 2006 vintage carried out after a CO₂ injection of 8.4 Mt. In spite of some difficulties due to the lack of time‐lapse well log data on the interest area, the full application of our inversion workflow allowed us to obtain, for the first time to our knowledge, 3D impedance cubes including the Utsira Sand. These results can be used to better characterize the spreading of CO₂ in a reservoir. With the post‐stack inversion workflow applied to CO₂ storage, we point out the importance of the a priori model and the issue to obtain coherent results between sequential inversions of different seismic vintages. The stacking velocity workflow that yields the migration model and the a priori model, specific to each vintage, can induce a slight inconsistency in the results.     

Adaptation to an uncertain climate change: cost benefit analysis and robust decision making for dam dimensioning

Climate models project large changes in rainfall, but disagree on their magnitude and sign. The consequences of this uncertainty on optimal dam dimensioning is assessed for a small mountainous catchment in Greece. Optimal dam design is estimated using a Cost-Benefit Analysis (CBA) based on trends in seasonal temperature and precipitations from 19 IPCC-AR4 climate models driven by the the SRES A2 emission scenario. Optimal reservoir volumes are modified by climate change, leading to up to 34% differences between optimal volumes. Contrary to widely-used target-based approaches, the CBA suggests that reduced rainfall should lead to smaller water reservoirs. The resulting change in the Net Present Value (NPV) of water supply is also substantial, ranging from no change to a large 25% loss, depending on the climate model, even assuming optimal adaptation and perfect foresight. In addition, climate change uncertainty can lead to design errors, with a cost ranging from 0.3 to 2.8% of the NPV, depending on site characteristics. This paper proposes to complement the CBA with a robust decision-making approach that focuses on reducing design-error costs. It also suggests that climate change impacts in the water sector may reveal large, that water reservoirs do not always provide a cost-efficient adaptation strategy, and that alternative adaptation strategies based on water conservation and non-conventional water production need to be considered.     

UVolution, a photochemistry experiment in low earth orbit: Investigation of the photostability of carbonates exposed to martian-like UV radiation conditions

The detection and identification of carbonates on Mars are of prime importance to establish the evolution of its atmosphere, correlated to the history of the liquid water, or even to determine the existence of a possible ancient biological activity. Till date, no large deposits of carbonates have been found. In fact, their detection is specific to local areas and in very low amounts. The absence of such deposits is commonly attributed to the harsh environmental conditions at the surface of Mars. Additionally, the presence of UV radiation has been proposed to explain their photodecomposition and hence their absence. However, contradictory results from laboratory experiments mimicking Mars’ surface UV radiation did not resolve the behaviour of carbonates in such an environment, which is why we exposed, in low Earth orbit and in laboratory experiments, both abiotic and biotic calcium carbonates to UV radiation of wavelength above 200 nm, the same spectral distribution as the one reaching the surface of Mars. For low Earth orbit (LEO) exposure, this was done for the UVolution experiment on board the BIOPAN ESA module, which was set outside a Russian Foton automated capsule, and exposed to space conditions for 12 days in September 2007. The targeted carbonates are biominerals and abiotic samples. Our laboratory results mainly show that the exposed carbonates appear to be stable to UV radiation if directly exposed to it. The LEO experiment results tend to the same conclusion, but the integrated exposition time to Solar UV during the experiment is not sufficient to be conclusive. However, the stability of the biominerals derived from the laboratory experiment could strengthen the interest to explore deeper their potential as life records at Mars. Hence, they should be considered as primary targets for in situ analyses during future missions.     

Operator-oriented CRS interpolation

In common‐reflection‐surface imaging the reflection arrival time field is parameterized by operators that are of higher dimension or order than in conventional methods. Using the common‐reflection‐surface approach locally in the unmigrated prestack data domain opens a potential for trace regularization and interpolation. In most data interpolation methods based on local coherency estimation, a single operator is designed for a target sample and the output amplitude is defined as a weighted average along the operator. This approach may fail in presence of interfering events or strong amplitude and phase variations. In this paper we introduce an alternative scheme in which there is no need for an operator to be defined at the target sample itself. Instead, the amplitude at a target sample is constructed from multiple operators estimated at different positions. In this case one operator may contribute to the construction of several target samples. Vice versa, a target sample might receive contributions from different operators. Operators are determined on a grid which can be sparser than the output grid. This allows to dramatically decrease the computational costs. In addition, the use of multiple operators for a single target sample stabilizes the interpolation results and implicitly allows several contributions in case of interfering events. Due to the considerable computational expense, common‐reflection‐surface interpolation is limited to work in subsets of the prestack data. We present the general workflow of a common‐reflection‐surface‐based regularization/interpolation for 3D data volumes. This workflow has been applied to an OBC common‐receiver volume and binned common‐offset subsets of a 3D marine data set. The impact of a common‐reflection‐surface regularization is demonstrated by means of a subsequent time migration. In comparison to the time migrations of the original and DMO‐interpolated data, the results show particular improvements in view of the continuity of reflections events. This gain is confirmed by an automatic picking of a horizon in the stacked time migrations.     

LINET—An international lightning detection network in Europe

During the past years a VLF/LF lightning detection network (LINET) was developed at the University of Munich, which provides continuous data for both research and operational purposes. In particular, the network introduces five new features: a) total lightning capability: both cloud-to-ground strokes (CG) and cloud lightning (IC) are measured; b) low-amplitude reporting: weak lightning events from discharge channel with currents well below 5 kA are detected within the central part of the network, whereby IC events dominate; c) new 3D-discrimination: a time-of-arrival method is utilized to separate CG from IC with good reliability, provided that the sensor baseline does not exceed ~ 250 km; d) IC emission height: for each cloud event a height is determined which is thought to reflect the central region of the involved channel; and e) optimised location accuracy: due to precision and combined action of all influential network components, complemented by site-error corrections, the position accuracy of strokes reaches an average value as small as ~ 150 m, whereby false locations (‘outliers’) rarely occur. During international co-operations LINET has been deployed in four continents: Europe (initially Germany), South America (area of Bauru, Brazil), Australia (around Darwin), and Central Africa (Benin). Since the features quoted above could be verified in the tests, a 65-sensor network was established in Europe and started on May 1, 2006, in co-operation with the service company nowcast. LINET covers a wide area approximately from longitude − 10° to 25° to latitude 35° to 66°; it is available for scientific projects and officially utilized by the German Weather Service for operational purposes. Meanwhile, the network was extended by deployment of additional sites so that it comprises about 90 sensors in 17 countries.     

Deep crustal source of gneiss dome revealed by eclogite in migmatite (Montagne Noire,French Massif Central)

In orogens worldwide and throughout geologic time, large volumes of deep continental crust have been exhumed in domal structures. Extension-driven ascent of bodies of deep, hot crust is a very efficient mechanism for rapid heat and mass transfer from deep to shallow crustal levels and is therefore an important mechanism in the evolution of continents. The dominant rock type in exhumed domes is quartzofeldspathic gneiss (typically migmatitic) that does not record its former high-pressure (HP) conditions in its equilibrium mineral assemblage; rather, it records the conditions of emplacement and cooling in the mid/shallow crust. Mafic rocks included in gneiss may, however, contain a fragmentary record of a HP history, and are evidence that their host rocks were also deeply sourced. An excellent example of exhumed deep crust that retains a partial HP record is in the Montagne Noire dome, French Massif Central, which contains well-preserved eclogite (garnet+omphacite+rutile+quartz) in migmatite in two locations: one in the dome core and the other at the dome margin. Both eclogites record P ~ 1.5 ± 0.2 GPa at T  ~  700 ± 20°C, but differ from each other in whole-rock and mineral composition, deformation features (shape and crystallographic preferred orientation, CPO), extent of record of prograde metamorphism in garnet and zircon, and degree of preservation of inherited zircon. Rim ages of zircon in both eclogites overlap with the oldest crystallization ages of host gneiss at c. 310 Ma, interpreted based on zircon rare earth element abundance in eclogite zircon as the age of HP metamorphism. Dome-margin eclogite zircon retains a widespread record of protolith age (c. 470–450 Ma, the same as host gneiss protolith age), whereas dome-core eclogite zircon has more scarce preservation of inherited zircon. Possible explanations for differences in the two eclogites relate to differences in the protolith mafic magma composition and history and/or the duration of metamorphic heating and extent of interaction with aqueous fluid, affecting zircon crystallization. Differences in HP deformation fabrics may relate to the position of the eclogite facies rocks relative to zones of transpression and transtension at an early stage of dome development. Regardless of differences, both eclogites experienced HP metamorphism and deformation in the deep crust at c. 310 Ma and were exhumed by lithospheric extension—with their host migmatite—near the end of the Variscan orogeny. The deep crust in this region was rapidly exhumed from ~50 to <10 km, where it equilibrated under low-P/high-T conditions, leaving a sparse but compelling record of the deep origin of most of the crust now exposed in the dome.     

Low-cost inundation modelling at the reach scale with sparse data in the Lower Damodar River basin,India

Abstract

Data unavailability is the main reason for limited applications of hydrodynamic models for predicting inundation in the developing world. This paper aims to generate moderately high-resolution hybrid terrain data by merging height information from low-cost Indian Remote Sensing satellite (IRS) Cartosat-1 stereo satellite images, freely-available Shuttle Radar Topograph Mission (SRTM) digital elevation model (DEM) data, and limited surveyed channel cross-sections. The study reach is characterized by anabranching channels that are associated with channel bifurcation, loops and river islands. We compared the performance of a simple 1D–2D coupled LISFLOOD-FP model and a complex fully 2D finite element TELEMAC-2D model with the hybrid terrain data. The results show that TELEMAC-2D produced significantly improved simulated inundation with the hybrid terrain data, as compared to the SRTM DEM. LISFLOOD-FP was found unsuitable to work with the hybrid DEM in a complicated fluvial environment, as it failed to efficiently divert water in the branches from the main channel.

Editor D. Koutsoyiannis; Associate editor A. Viglione     

New constrains on the thickness of the Semail ophiolite in the Northern Emirates

Near-critical angle and refraction studies were performed at IFP as piggyback studies during a wider programme of crustal imagery operated by WesternGeco on behalf of the Ministry of Energy of the United Arab Emirates. The main objective is to illuminate the base of the Semail Ophiolite along part of a regional transect (D1) crossing the Northern Emirates from the Gulf of Oman in the east up to the Arabian Gulf in the west. Results confirm that the sole thrust of the ophiolite has been folded during the Miocene stacking of the underlying Arabian Platform. The thickness of the ophiolite grades from zero in the core of the Masafi tectonic window, up to a maximum of 1.7 km below the axial part of a successor basin which has been preserved on top of the serpentinite west of the current exposure of the main ultramafic bodies. Apatite grains extracted from plagiogranites of the Semail ophiolite also provide evidences for an early unroofing of the gabbros and plagiogranites during the Late Cretaceous, with cooling ages of 72–76 Ma at the top of the ophiolite in the east (not far from the Fujairah coast line), which are coeval and also consistent with the occurrence of Late Cretaceous paleo-soils, rudists and paleo-reef deposits on top of serpentinized ultramafics in the west. Younger cooling ages of 20 Ma have been also found at the base of the ophiolite near Masafi, in the core of the nappe anticline, thus providing a Neogene age for the refolding of the allochthon and stacking of underlying parautochthonous platform carbonate units. These results, together with the occurrence of a thick sedimentary pile illuminated below the metamorphic sole along the north-trending, strike-profile D2 running parallel to the axis of the Masafi window, should stimulate a renewal of the exploration in the central part of the Emirate foothills, where the ophiolite thickness is currently limited, and was already drastically reduced by the end of its Late Cretaceous obduction.