Gravimetric monitoring of the first field‐wide steam injection in a fractured carbonate field in Oman – a feasibility study

Gas‐Oil Gravity Drainage is to be enhanced by steam injection in a highly fractured, low permeability carbonate field in Oman. Following a successful pilot, field‐wide steam injection is being implemented, first of this type in the world. A dedicated monitoring program has been designed to track changes in the reservoir. Various observations are to be acquired, including, surface deformation, temperature measurements, microseismic, well logs, pressure and saturation measurements to monitor the reservoir. Because significant changes in the reservoir density are expected, time‐lapse gravimetry is also being considered. In this paper we investigate the feasibility of gravimetric monitoring of the thermally enhanced gravity drainage process at the carbonate field in Oman. For this purpose, forward gravity modelling is performed. Based on field groundwater measurements, the estimates of the hydrological signal are considered and it is investigated under what conditions the groundwater influences can be minimized. Using regularized inversion of synthetic gravity data, we analyse the achievable accuracy of heat‐front position estimates. In case of large groundwater variations at the field, the gravity observations can be significantly affected and, consequently, the accuracy of heat‐front monitoring can be deteriorated. We show that, by applying gravity corrections based on local observations of groundwater, the hydrological influences can to a large extent be reduced and the accuracy of estimates can be improved. We conclude that gravimetric monitoring of the heat‐front evolution has a great potential.     

Strain parameters of discontinuities in rock for finite element calculation

The bulk strain induced by excavation of rooms for storage of highly radioactive waste is usually unimportant, but the change in aperture and shear displacement of permeable fractures may be of practical importance. Hence, the normal and shear of such discontinuities have to be considered and are commonly predicted by using special deformation moduli termed “normal” and “shear” stiffnesses in numerical calculations. For the use of some finite element methods, these have to be converted to compression and shear moduli ,which, in turn, requires that the thickness of the discontinuity is known. The normal and shear strain can be expressed in terms of soil mechanical parameters, yielding the compression modulus M and the shear modulus G, with the required form for finite element calculations. By definition they are functions of the normal and shear stiffness K_n and K_s. It is concluded that calculation of the normal and shear strain of discontinuities with clastic fillings by the use of finite element methods and deformation moduli derived from stiffness numbers is very uncertain except when the geometry of the weaknesses can be accurately defined.     

Post-seismic relaxation following the great 2004 Sumatra-Andaman earthquake on a compressible self-gravitating Earth

The   M _w γ 9.0  2004 December 26 Sumatra-Andaman and   M _w = 8.7  2005 March 28 Nias earthquakes, which collectively ruptured approximately 1800 km of the Andaman and Sunda subduction zones, are expected to be followed by vigorous viscoelastic relaxation involving both the upper and lower mantle. Because of these large spatial dimensions it is desirable to fully account for gravitational coupling effects in the relaxation process. We present a stable method of computing relaxation of a spherically-stratified, compressible and self-gravitating viscoelastic Earth following an impulsive moment release event. The solution is cast in terms of a spherical harmonic expansion of viscoelastic normal modes. For simple layered viscoelastic models, which include a low-viscosity oceanic asthenosphere, we predict substantial post-seismic effects over a region several 100s of km wide surrounding the eastern Indian Ocean. We compare observed GPS time-series from ten regional sites (mostly in Thailand and Indonesia), beginning in 2004 December, with synthetic time-series that include the coseismic and post-seismic effects of the 2004 December 26 and 2005 March 28 earthquakes. A viscosity structure involving a biviscous (Burgers body) rheology in the asthenosphere explains the pattern and amplitude of post-seismic offsets remarkably well.     

Geometry and chronology of growth and drowning of Middle Triassic carbonate platforms (Cernera and Bivera/Clapsavon) in the Southern Alps (northern Italy)

The depositional architecture and the geometric relationships between platform-slope deposits and basinal sediments along with paleontological evidence indicate the time interval of the younger Anisian Reitziites reitzi ammonoid zone to largely represent the main stage of platform aggradation at the Cernera and Bivera/Clapsavon carbonate platforms. Published and new U-Pb age data of zircons from volcaniclastic layers bracketing the stratigraphic interval of platform growth constrain the duration of platform evolution to a time span shorter than 1.8±0.7m.y., probably in the order of 0.5-1m.y., reflecting fast rates of vertical platform aggradation exceeding 500 m/m.y. In the range of growth potentials for shallow-water carbonate systems estimated in relation to the time span of observation, this high rate is in agreement with values for short intervals of 10⁵-10⁶yrs (e.g., Schlager 1999). After drowning, the platforms at Cernera and Bivera/Clapsavon were blanketed by thin pelagic carbonates. On the former platform flanks the draping sediments in places comprise red nodular pelagic limestones (Clapsavon Limestone) similar in facies to the Han Bulog Limestones occurring elsewhere in Middle Triassic successions of the Mediterranean Tethys. The drowning of vast areas of former carbonate platforms possibly triggered the onset of bottom-water circulation in adjacent basins as suggested by the abrupt transition from laminated to bioturbated pelagic nodular limestones in the Buchenstein Formation which occurred close to the time of initial platform submergence. During the Late Ladinian the topographic features of the drowned platforms were onlapped by rapidly deposited, predominantly clastic successions including coarse breccias and volcanic rocks sealing and preserving the peculiar stratigraphic setting.     

The GeoForschungsZentrum Potsdam/Groupe de Recherche de Gèodésie Spatiale satellite-only and combined gravity field models: EIGEN-GL04S1 and EIGEN-GL04C

The recent improvements in the Gravity Recovery And Climate Experiment (GRACE) tracking data processing at GeoForschungsZentrum Potsdam (GFZ) and Groupe de Recherche de Géodésie Spatiale (GRGS) Toulouse, the availability of newer surface gravity data sets in the Arctic, Antarctica and North-America, and the availability of a new mean sea surface height model from altimetry processing at GFZ gave rise to the generation of two new global gravity field models. The first, EIGEN-GL04S1, a satellite-only model complete to degree and order 150 in terms of spherical harmonics, was derived by combination of the latest GFZ Potsdam GRACE-only (EIGEN-GRACE04S) and GRGS Toulouse GRACE/LAGEOS (EIGEN-GL04S) mean field solutions. The second, EIGEN-GL04S1 was combined with surface gravity data from altimetry over the oceans and gravimetry over the continents to derive a new high-resolution global gravity field model called EIGEN-GL04C. This model is complete to degree and order 360 and thus resolves geoid and gravity anomalies at half- wavelengths of 55 km at the equator. A degree-dependent combination method has been applied in order to preserve the high accuracy from the GRACE satellite data in the lower frequency band of the geopotential and to form a smooth transition to the high-frequency information coming from the surface data. Compared to pre-CHAMP global high-resolution models, the accuracy was improved at a spatial resolution of 200 km (half-wavelength) by one order of magnitude to 3 cm in terms of geoid heights. The accuracy of this model (i.e. the commission error) at its full spatial resolution is estimated to be 15 cm. The model shows a reduced artificial meridional striping and an increased correlation of EIGEN-GL04C-derived geostrophic meridional currents with World Ocean Atlas 2001 (WOA01) data. These improvements have led to select EIGEN-GL04C for JASON-1 satellite altimeter data reprocessing. Electronic Supplementary Material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Lake Sediment Core Records of Sulphur Accumulation and Sulphur Isotopic Composition in Central Ontario, Canada Lakes

Stable isotopic compositions and concentrations of total sedimentary sulphur (S) were determined in cores from 6 lakes in the acid-sensitive Muskoka-Haliburton region of south-central Ontario. The isotopic composition of S in deep sediment (> ~ 20 cm) was approximately constant in all lakes, and indicated a pre-industrial δ ³⁴S value between +4.0 and +5.3‰, which is similar to current bulk deposition. Similarly, total S concentrations in deep sediment were relatively low (1.9–5 mg S g⁻¹ dwt) and approximately constant with depth within cores. All lakes exhibited up-core increases in total S and decreases in δ ³⁴S at a depth corresponding to the beginning of industrialization in the Great Lakes region ( ~ 1900), resulting in a generally reciprocal depth pattern between total S concentration and δ ³⁴S ratios. While initial shifts in total S and δ ³⁴S were likely due to enhanced SO₄ reduction of newly available anthropogenic SO₄, both the magnitude and pattern of up-core S enrichment and shifts in δ ³⁴S varied greatly among lakes, and did not match changes in S deposition post 1900. Differences between lakes in total S and δ ³⁴S were not related to any single hydrologic (e.g., residence time) or physical (e.g., catchment-area-to-lake area ratio) lake characteristic. This work indicates that sediment cores do not provide consistent records of changes in post-industrial S deposition in this region, likely due to redox-related mobility of S in upper sediment.     

The potential use of OH-defects in enstatite as geobarometer

In this study, single crystals of pure enstatite (Mg₂Si₂O₆) were synthesised under water-saturated conditions at 4 and 8 GPa and 1,150°C with variable silica activity, leading to phase assemblages enstatite + forsterite, enstatite or enstatite + coesite. Run products were investigated using an FTIR spectrometer equipped with a focal plane array detector enabling IR imaging with a lateral pixel resolution of 2.7 μm. IR spectra within the OH-absorption region show two different groups of absorption bands: group 1 (wavenumbers at 3,592 and 3,687 cm^?1) shows strongest absorptions for E||n _β, whereas group 2 (wavenumbers at 3,067 and 3,362 cm^?1) shows strongest absorptions for E||n _γ. The groups are related to different defect types, group 1 to tetrahedral defects (T-site vacancies) and group 2 to octahedral defects (M-site vacancies). The intensity ratio of the bands within one group (i.e. A ₃₀₆₇/A ₃₃₆₂ and A ₃₅₉₂/A ₃₆₈₇) and the intensity ratio of E||n _γ and E||n _α in group 2 bands remain constant within error. In contrast, the intensity ratio of group 2 to group 1 absorption bands [e.g. (A ₃₃₆₂)/(A ₃₆₈₇)] is sensitive to the SiO₂ activity and pressure. On the basis of the results of this and previous studies, a barometer for pure orthoenstatite coexisting with forsterite can be formulated:\( P\,[{\text{GPa}}] = 1.056 + \sqrt {{\frac{{1.025 - A_{{\left( {3362} \right)/\left[ {(3362) + (3687)} \right]}} }}{0.009}}} , \) where A ₍₃₃₆₂₎ and A ₍₃₆₈₇₎ are the integral absorbances of the component E||n _γ of the absorption bands at 3,362 cm^?1 and the component E||n _β of the absorption band at 3,687 cm^?1, respectively.