Upscaling: Effective Medium Theory, Numerical Methods and the Fractal Dream

Upscaling is a major issue regarding mechanical and transport properties of rocks. This paper examines three issues relative to upscaling. The first one is a brief overview of Effective Medium Theory (EMT), which is a key tool to predict average rock properties at a macroscopic scale in the case of a statistically homogeneous medium. EMT is of particular interest in the calculation of elastic properties. As discussed in this paper, EMT can thus provide a possible way to perform upscaling, although it is by no means the only one, and in particular it is irrelevant if the medium does not adhere to statistical homogeneity. This last circumstance is examined in part two of the paper. We focus on the example of constructing a hydrocarbon reservoir model. Such a construction is a required step in the process of making reasonable predictions for oil production. Taking into account rock permeability, lithological units and various structural discontinuities at different scales is part of this construction. The result is that stochastic reservoir models are built that rely on various numerical upscaling methods. These methods are reviewed. They provide techniques which make it possible to deal with upscaling on a general basis. Finally, a last case in which upscaling is trivial is considered in the third part of the paper. This is the fractal case. Fractal models have become popular precisely because they are free of the assumption of statistical homogeneity and yet do not involve numerical methods. It is suggested that using a physical criterion as a means to discriminate whether fractality is a dream or reality would be more satisfactory than relying on a limited data set alone.     

Effects of economic recession and local weather on climate change attitudes

How are large companies responding to the challenges of reducing their greenhouse gas (GHG) emissions? An analysis of the published climate change policies and performance of 125 large European companies is presented. The results suggest that most large European companies have now developed the management systems and processes necessary for them to effectively manage their GHG emissions and related business risks. However, there is a significant disconnect between the targets that companies set for themselves and the more ambitious targets being set by the European Union (which has committed to a 20% reduction in its emissions by 2020 against a 1990 baseline). Of the companies surveyed, just over one-third had stabilized or reduced their total GHG emissions over the period 2002–2007, and fewer than one-third expected their emissions to stabilize or reduce in the coming years. The relationship between the quality of corporate policies and performance outcomes (in terms of GHG emissions) suggests that while companies with stronger policies are likely to have relatively better performance, only a minority of those companies with the highest-quality policies are committing to absolute reductions in their GHG emissions.     

Formation of Carbon and Hydrogen Species in Magmas at Low Oxygen Fugacity

Studies of iron-bearing silicate melt (ferrobasalt) + iron metallicphase + graphite + hydrogen equilibria show that carbon andhydrogen solubilities in melts are important for the evolutionof the upper mantle. In a series of experiments conducted at3·7 GPa and 1520–1600°C, we have characterizedthe nature (oxidized vs reduced) and quantified the abundancesof C- and H-compounds dissolved in iron-bearing silicate melts.Experiments were carried out in an anvil-with-hole apparatuspermitting the achievement of equal chemical potentials of H₂in the inner Pt capsule and outer furnace assembly. The fO₂for silicate melt–iron equilibrium was 2·32 ±0·04 log units below iron–wüstite (IW). Theferrobasalt used as starting material experienced a reductionof its iron oxides and silicate network. The counterpart wasa liberation of oxygen reacting with the hydrogen entering thecapsule. The amount of H₂O dissolved in the glasses was measuredby ion microprobe and by step-heating and was found to be between1 and 2 wt %. The dissolved carbon content was found to be 1600ppm C by step-heating. The speciation of C and H componentswas determined by IR and Raman spectroscopy. It was establishedthat the main part of the liberated oxygen was used to formOH^– and to a much lesser extent H₂O, and only traces ofH₂, CO₂ and     

The relationship between lightning activity and surface wet bulb temperature and its variation with latitude in Australia

Summary Convective Available Potential Energy (CAPE) is the driving force for thunderstorm development. CAPE is closely controlled by wet bulb temperature. The lightning activity measured by a network of ten lightning flash counters widely distributed across continental Australia was studied as a function of wet bulb temperature. At each of the stations, the monthly total of lightning ground flashes, N, increased sharply with the increase of the monthly mean daily maximum wet bulb temperature, T_w, max. The dependence was strongest in the tropics and became less pronounced at temperate latitudes. In Darwin (latitude 12° S), the lightning ground flash activity increased by over three orders of magnitude over a 7 °C range of T_w, max. The corresponding increases for Coffs Harbour (latitude 30° S) and for Melbourne (latitude 38° S) were about one and a half orders of magnitude and about half an order of magnitude, respectively, each over a 10 °C range of T_w, max. Power law approximations were derived for each of the ten stations and showed that the logarithm of N was directly proportional to the power, P, of T_w, max. The value of P showed a sharp exponential decrease with increasing latitude away from the equator.     

Energetic Particles Observed in the CUSP Region During a Storm Recovery Phase

In this paper we report energetic ion behavior and its composition variations observed by the Cluster/RAPID instrument when the spacecraft was travelling in the high latitude magnetospheric boundary region on the day of the 31 March, 2001, strongest magnetic storm in the past 50 years. The Dst index reached −360 nT at about 09:00 UT. During its early recovery phase, large amounts of oxygen and helium ions were observed; the ratio of oxygen to hydrogen in the RAPID energy range reached as high as 250%, which suggests that the observed energetic particles might be of magnetospheric origin. The observations further show that enhanced energetic electron fluxes are confined in a very narrow region, while protons have occupied a larger region, and heavy ions have been observed in an even larger region. The flux of energetic electrons show a slight enhancement in a region where the magnetic field magnitude is around zero. These observed energetic ions could be quasi-trapped by the current sheet in the stagnation region of the cusp.     

Lithology-derived structure classification from the joint interpretation of magnetotelluric and seismic models

Magnetotelluric and seismic methods provide complementary information about the resistivity and velocity structure of the subsurface on similar scales and resolutions. No global relation, however, exists between these parameters, and correlations are often valid for only a limited target area. Independently derived inverse models from these methods can be combined using a classification approach to map geologic structure. The method employed is based solely on the statistical correlation of physical properties in a joint parameter space and is independent of theoretical or empirical relations linking electrical and seismic parameters. Regions of high correlation (classes) between resistivity and velocity can in turn be mapped back and re-examined in depth section. The spatial distribution of these classes, and the boundaries between them, provide structural information not evident in the individual models. This method is applied to a 10 km long profile crossing the Dead Sea Transform in Jordan. Several prominent classes are identified with specific lithologies in accordance with local geology. An abrupt change in lithology across the fault, together with vertical uplift of the basement suggest the fault is sub-vertical within the upper crust.     

Site classification of Turkish national strong-motion stations

Since 1973, the General Directorate of Disaster Affairs of Turkey has deployed several strong-motion accelerographs at selected sites. Within the framework of the project entitled Compilation of National Strong Ground Motion Database in Accordance with International Standards, site conditions were investigated within the upper 30-m depth by surface seismic and standard penetration tests. Preliminary characterization of the sites is made by making use of both geophysical and geotechnical criteria of NEHRP Provisions and Eurocode-8 site classification systems. The liquefaction susceptibility of those sites which comprise saturated cohesionless deposits is also determined. Mean shear-wave velocity, mean penetration resistance, site class, and liquefaction susceptibility of each site are tabulated. The Turkish strong-motion database supplemented by detailed information on site conditions is a valuable source of information particularly for those studies that put emphasis on the relationship between site conditions and strong-motion parameters.     

Spectral classification methods in monitoring small local events by the Israel seismic network

We use the dense Israel Seismic Network (ISN) to discriminate between low magnitude earthquakes and explosions in the Middle East region. This issue is important for CTBT monitoring, especially when considering small nuclear tests which may be conducted under evasive conditions. We explore the performance of efficient discriminants based on spectral features of seismograms using waveforms of 50 earthquakes and 114 quarry and underwater blasts with magnitudes 1.0–2.8, recorded by ISN short-period stations at distances up to 200 km. The single-station spectral ratio of the low and high-frequency seismic energy shows an overlap between explosions and earthquakes. After averaging over a subnet of stations, the resolving power is enhanced and the two classes of events are separated. Different frequency bands were tested; the (1–3 Hz)/(6–8 Hz) ratio provided the best discriminant performance. We also estimated normalized r.m.s. spectral amplitudes in several sequential equal frequency windows within the 1–12 Hz band and applied multiparametric automatic classification procedures (Linear Discrimination Function and Artificial Neural Network) to the amplitudes averaged over a subnetwork. A leave-one-out test showed a low rate of error for the multiparametric procedures. An innovative multi-station discriminant is proposed, based on spectral modulation associated with ripple-firing in quarry blasts and with the bubbling effect in underwater explosions. It utilizes a distinct azimuth-invariant coherency of spectral shapes for different stations in the frequency range (1–12 Hz). The coherency is measured by semblance statistics commonly used in seismic prospecting for phase correlation in the time domain. After modification, the statistics applied to the network spectra provided event separation. A new feature of all the above mentioned procedures is that they are based on smoothed (0.5 Hz window), instrument-corrected FFT spectra of the whole signal; they are robust to the accuracy of onset time estimation and, thus well suited to automatic event identification.     

Melting of Mg2GeO4 under pressure

Melting point of germanate forsterite, Mg₂GeO₄, was raised by compression at the rate of 30°C/GPa. The triple point, at which three phases of olivine- and spinel-type solids and liquid coexisted, was fixed at 1950°C and 3.5GPa. Wen these results are combined with the thermodynamical data of forsterite, Mg₂SiO₄, it is estimated that the triple point of forsterite lies in a region ranging from 2700° to 3000°C in temperature and from 20 to 30GPa in pressure.     

The Variscan tectonic inheritance of the Upper Rhine Graben: evidence of reactivations in the Lias,Late Eocene–Oligocene up to the recent

Processing of gravity and magnetic maps shows that the basement of the Upper Rhine Graben area is characterized by a series of NE–SW trending discontinuities and elongated structures, identified in outcrops in the Vosges, Black Forest, and the Odenwald Mountains. They form a 40 km wide, N30–40° striking, sinistral wrench-zone that, in the Visean, shifted the Variscan and pre-Variscan structures by at least 43 km to the NE. Wrenching was associated with emplacement of several generations of plutonic bodies emplaced in the time range 340–325 Ma. The sub-vertical, NE–SW trending discontinuities in the basement acted as zones of weakness, susceptible to reactivation by subsequent tectonism. The first reactivation, marked by mineralizations and palaeomagnetic overprinting along NE–SW faults of the Vosges Mountains, results from the Liassic NW–SE extension contemporaneous with the break-up of Pangea. The major reactivation occurred during the Late Eocene N–S compression and the Early-Middle Oligocene E–W extension. The NE–SW striking basement discontinuities were successively reactivated as sinistral strike-slip faults, and as oblique normal faults. Elongated depocenters appear to form in association with reactivated Variscan wrench faults. Some of the recent earthquakes are located on NE–SW striking Variscan fault zones, and show sinistral strike-slip focal mechanisms with the same direction, suggesting also present reactivation.