New models for heater wells in subsurface simulations, with application to the in situ upgrading of oil shale

Downhole electrical heating can be used to achieve the high temperatures required for in situ upgrading of oil shale or oil sands. Heater-well models are needed if this process is to be simulated accurately. The traditional Peaceman approach used for fluid injection and production wells may not be applicable because it does not capture transient effects, which can be important in downhole heater models. Standard models also neglect the effects of heterogeneity and temperature dependence in the rock properties. Here, we develop two new models for representing heater wells in reservoir simulators. The first model is applicable for homogeneous systems with properties that are not temperature dependent. For such cases, we develop a semi-analytical procedure based on Green’s functions to construct time-dependent heater-well indexes and heater-block thermal transmissibilities. For the general case, which can include both fine-scale heterogeneity and nonlinearity due to the temperature dependence of rock properties, we present a numerical procedure for constructing the heater-well model. This technique is essentially a near-well upscaling method and requires a local fine-scale solution in the near-well region. The boundary conditions are determined using a local-global treatment. The accuracy of the new heater-well models is demonstrated through comparison to reference solutions for example problems. The approach is then applied for the coarse-scale modeling of the in situ upgrading of oil shale, which entails a thermal-compositional simulation with chemical reactions. The model is shown to provide an accurate and efficient solution for this challenging problem.     

Diachronic metamorphic and structural evolution of the Connecticut Valley–Gaspé trough,Northern Appalachians

The Connecticut Valley–Gaspé (CVG) trough represents a major, orogen-scale Silurian–Devonian basin of the Northern Appalachians. From Gaspé Peninsula to southern New England, the CVG trough has experienced a contrasting metamorphic and structural evolution during the Acadian orogeny. Along its strike, the CVG trough is characterized by increasing strain and polyphase structures, and by variations in the intensity of regional metamorphism and contrasting abundance of c. 390–370 Ma granitic intrusions. In southern Quebec and northern Vermont, a series of NW–SE transects across the CVG trough have been studied in order to better understand these along-strike variations. Detailed structural analyses, combined with phase equilibria modelling, Raman spectrometry, and muscovite ⁴⁰Ar/³⁹Ar dating highlight a progressive and incremental deformation involving south–north variation in the timing of metamorphism. Deformation evolves from a D₁ crustal thickening event which originates in Vermont and progresses to southern Québec where it peaked at 0.6 GPa/380°C at c. 375 Ma. This was followed by a D₂ event associated with continuous burial in Vermont from 378 to 355 Ma, which produced peak metamorphic conditions of 0.85 GPa/380°C and exhumation in Quebec from 368 to 360 Ma. The D₃ compressional exhumation event also evolved from south to north from 345 to 335 Ma. D₁ to D₃ deformation events form part of a continuum with an along-strike propagation rate of ~50 km/Ma During D₁, the burial depth varied by more than 15 km between southern Quebec and Vermont, and this can be attributed to the occurrence of a major crustal indenter, the Bronson Hill Arc massif, in the New England segment of the Acadian collision zone.     

Revisiting Multi-Gaussian Kriging with the Nataf Transformation or the Bayes’ Rule for the Estimation of Spatial Distributions

A multivariate probability transformation between random variables, known as the Nataf transformation, is shown to be the appropriate transformation for multi-Gaussian kriging. It assumes a diagonal Jacobian matrix for the transformation of the random variables between the original space and the Gaussian space. This allows writing the probability transformation between the local conditional probability density function in the original space and the local conditional Gaussian probability density function in the Gaussian space as a ratio equal to the ratio of their respective marginal distributions. Under stationarity, the marginal distribution in the original space is modeled from the data histogram. The stationary marginal standard Gaussian distribution is obtained from the normal scores of the data and the local conditional Gaussian distribution is modeled from the kriging mean and kriging variance of the normal scores of the data. The equality of ratios of distributions has the same form as the Bayes’ rule and the assumption of stationarity of the data histogram can be re-interpreted as the gathering of the prior distribution. Multi-Gaussian kriging can be re-interpreted as an updating of the data histogram by a Gaussian likelihood. The Bayes’ rule allows for an even more general interpretation of spatial estimation in terms of equality for the ratio of the conditional distribution over the marginal distribution in the original data uncertainty space with the same ratio for a model of uncertainty with a distribution that can be modeled using the mean and variance from direct kriging of the original data values. It is based on the principle of conservation of probability ratio and no transformation is required. The local conditional distribution has a variance that is data dependent. When used in sequential simulation mode, it reproduces histogram and variogram of the data, thus providing a new approach for direct simulation in the original value space.     

Alternative technologies for rural areas–what about the ‘alternative’ dimension of Wi-Fi?

Summary  The question of alternative technologies for high-speed Internet access lies at the heart of rural development problems. In this paper, we focus on one of those technologies, the system combining satellite technologies and Wi-Fi. Based on an empirical study carried out in three rural areas, we analyze the dynamics of the use and appropriation of that technology by the companies and organizations participating in the experiment. Considered both from a technical and social standpoint, the technology ‘in use’ appears here in its structuring dimensions.

Modelling crack propagation in concrete structures with a two scale approach

A simplified computational technique based on a refined global–local method is applied to the failure analysis of concrete structures. The technique distinguishes the scale of the structure, modelled with large size finite elements, from the scale at which material non‐linearity occurs due to progressive cracking and macro‐crack propagation. The finite element solution is split into two parts: a linear elastic analysis on a coarse mesh over the entire structure and a non‐linear analysis over a small part of the structure where a dense finite element grid is employed. In the non‐linear calculation, a non‐local damage model is implemented. These two computations are coupled with the help of an iterative scheme. The size and location of the region where a non‐linear analysis is performed, is adapted to follow the development of the damage zone. Numerical examples of mode I fracture of concrete specimens with straight and curved cracks are presented. Copyright © 2003 John Wiley & Sons, Ltd.     

Experimental analysis of compaction of concrete and mortar

Compaction of concrete is physically a collapse of the material porous microstructure. It produces plastic strains in the material and, at the same time, an increase of its bulk modulus. This paper presents two experimental techniques aimed at obtaining the hydrostatic response of concrete and mortar. The first one is a uniaxial confined compression test which is quite simple to implement and allows to reach hydrostatic pressures of about 600 MPa. The specimen size is large enough so that concrete with aggregate sizes up to 16 mm can be tested. The second one is a true hydrostatic test performed on smaller (mortar) specimens. Test results show that the hydrostatic response of the material is elasto‐plastic with a stiffening effect on both the tangent and unloading bulk moduli. The magnitude of the irreversible volumetric strains depends on the initial porosity of the material. This porosity can be related in a first approximation to the water/cement ratio. A comparison of the hydrostatic responses obtained from the two testing techniques on the same material show that the hydrostatic response of cementitious materials cannot be uncoupled from the deviatoric response, as opposed to the standard assumption in constitutive relations for metal alloys. This feature should be taken into account in the development of constitutive relations for concrete subjected to high confinement pressures which are needed in the modelling of impact problems. Copyright © 2001 John Wiley & Sons, Ltd.     

Heating and cooling of magnetars with accreted envelopes

We study the thermal structure and evolution of magnetars as cooling neutron stars with a phenomenological heat source in an internal layer. We focus on the effect of magnetized (   B ≳ 10¹⁴  G) non-accreted and accreted outermost envelopes composed of different elements, from iron to hydrogen or helium. We discuss a combined effect of thermal conduction and neutrino emission in the outer neutron star crust and calculate the cooling of magnetars with a dipole magnetic field for various locations of the heat layer, heat rates and magnetic field strengths. Combined effects of strong magnetic fields and light-element composition simplify the interpretation of magnetars in our model: these effects allow one to interpret observations assuming less extreme (therefore, more realistic) heating. Massive magnetars, with fast neutrino cooling in their cores, can have higher thermal surface luminosity.     

A seafloor experiment to monitor vertical deformation at the Lucky Strike volcano, Mid-Atlantic Ridge

Decades of cruise-based exploration have provided excellent snapshots of the structure of mid-ocean ridges and have revealed that accretion is a mixture of steady-state and quantum events. Observatory-type studies are now needed to quantify the temporal evolution of these systems. A multi-disciplinary seafloor observatory site is currently being set up at the Lucky Strike volcano, in the axial valley of the slow spreading Mid-Atlantic ridge as a part of the MoMAR (monitoring of the Mid-Atlantic Ridge) initiative. The aim of this observatory is to better understand the dynamics of the volcano and the hydrothermal vents hosted at its summit as well as their plumbing systems. In August 2006, the GRAVILUCK cruise initiated an experiment to monitor the deformation of Lucky Strike volcano. A geodetic network was installed, and seafloor pressure, gravity and magnetic data were collected. In this paper, we present the method used to monitor volcanic deformation, which involves measuring relative depth difference between points within a seafloor geodesy network. We show that, taking into account oceanographic variability and measurement noise, the network should be able to detect vertical deformations of the order of 1 cm.     

Trace element distribution in peridotite xenoliths from Tok, SE Siberian craton: A record of pervasive, multi-stage metasomatism in shallow refractory mantle

Spinel peridotite xenoliths in alkali basalts at Tok, SE Siberian craton range from fertile lherzolites to harzburgites and wehrlites; olivine-rich (70-84%) rocks are dominant. REE patterns in the lherzolites range from nearly flat for fertile rocks (14-17% cpx) to LREE-enriched; the enrichments are positively correlated with modal olivine, consistent with high-permeability of olivine-rich rocks during melt percolation. Clinopyroxene in olivine-rich Tok peridotites typically has convex-upward trace element patterns (La/Nd_PM < 1 and Nd/Yb_PM ? 1), which we consider as evidence for equilibration with evolved silicate liquids (with higher REE and lower Ti contents than in host basalts). Whole-rock patterns of the olivine-rich xenoliths range from convex-upward to LREE-enriched (La/Nd_PM > 1); the LREE-enrichments are positively correlated with phosphorus abundances and are mainly hosted by accessory phosphates and P-rich cryptocrystalline materials. In addition to apatite, some Tok xenoliths contain whitlockite (an anhydrous, halogen-poor and Na-Mg-rich phosphate), which is common in meteorites and lunar rocks, but has not been reported from any terrestrial mantle samples. Some olivine-rich peridotites have generations of clinopyroxene with distinct abundances of Na, LREE, Sr and Zr. The mineralogical and trace element data indicate that the lithospheric mantle section represented by the xenoliths experienced a large-scale metasomatic event produced by upward migration of mafic silicate melts followed by percolation of low-T, alkali-rich melts and fluids. Chromatographic fractionation and fractional crystallisation of the melts close to the percolation front produced strong LREE-enrichments, which are most common in the uppermost mantle and are related to carbonate- and P₂O₅-rich derivatives of the initial melt. Reversal and gradual retreat of the percolation front during thermal relaxation to ambient geotherm (“retrograde” metasomatism) caused local migration and entrapment of small-volume residual fluids and precipitation of volatile-rich accessory minerals. A distinct metasomatic episode, which mainly produced “anhydrous” late-stage interstitial materials was concomitant with the alkali basaltic magmatism, which brought the xenoliths to the surface.