Markov chain Monte Carlo methods for conditioning a permeability field to pressure data

Generating one realization of a random permeability field that is consistent with observed pressure data and a known variogram model is not a difficult problem. If, however, one wants to investigate the uncertainty of reservior behavior, one must generate a large number of realizations and ensure that the distribution of realizations properly reflects the uncertainty in reservoir properties. The most widely used method for conditioning permeability fields to production data has been the method of simulated annealing, in which practitioners attempt to minimize the difference between the ’ ’true and simulated production data, and “true” and simulated variograms. Unfortunately, the meaning of the resulting realization is not clear and the method can be extremely slow. In this paper, we present an alternative approach to generating realizations that are conditional to pressure data, focusing on the distribution of realizations and on the efficiency of the method. Under certain conditions that can be verified easily, the Markov chain Monte Carlo method is known to produce states whose frequencies of appearance correspond to a given probability distribution, so we use this method to generate the realizations. To make the method more efficient, we perturb the states in such a way that the variogram is satisfied automatically and the pressure data are approximately matched at every step. These perturbations make use of sensitivity coefficients calculated from the reservoir simulator.     

The use of multi‐frequency acquisition to significantly improve the quality of fibre‐optic‐distributed vibration sensing

Unlike conventional sensors that measure the passage of seismic waves at a single position, distributed vibration sensing systems, also known as distributed acoustic sensing systems, detect the passage of seismic waves by averaging a measurement of strain over a section of fibre‐optic cable. Distributed vibration sensing systems work by transmitting pulses of light down the fibre and measuring the phase of the Rayleigh backscatter. At random positions along the fibre, however, fading occurs; this is where the amplitude of the backscattered signal is very small due to cancellation of the scattered electric fields, resulting in anomalously noisy traces in a common source gather. This paper addresses the problem of fading in a particular form of distributed vibration sensors: a new optical arrangement of the instrumentation is described that allows the measurement to be carried out quasi‐simultaneously at multiple optical interrogation frequencies. The interrogation frequencies are chosen to be sufficiently different that their fading properties are distinct and the diversity thus obtained is used to aggregate the data obtained to substantially reduce the noise caused by fading. As well as reducing the effects of fading, the aggregation of the independent results can also help to reduce the overall noise of the measurement and improve the linearity of the distributed vibration sensing system.     

Solar Cycle Predictions (Invited Review)

Solar Physics - Solar cycle predictions are needed to plan long-term space missions, just as weather predictions are needed to plan the launch. Fleets of satellites circle the Earth collecting many...     

The role of oxyradicals in intracellular proteolysis and toxicity in mussels

Studies were performed on the common mussel, M. edulis L., to determine whether copper (Cu) exposure can affect the extent to which digestive cell proteins are oxidised and whether such oxidative damage is mediated by free radicals. Three age groups of mussels were exposed for 6 ^-days to environmentally realistic concentrations of Cu and then digestive gland homogenates were examined for evidence of protein carbonyl formation. Significant increases in carbonyls relative to untreated control mussels were seen for the youngest (2–4 year-old) and oldest (≥ 10 year-old) mussels only after exposure for 6 days, followed by recovery from exposure for a further 6 days. Untreated mussels also showed an age-related difference in protein oxidation, with a significantly lower concentration in the youngest animals (2–4 year olds). Copper did not affect the levels of modified tryptophan or tyrosine residues or the extent of total lipid peroxidation in digestive gland homogenate. Significant depletion of total vitamin E (a-tocopherol) was seen only in young and medium-aged mussels following exposure for 6 days. The levels of protein carbonyl groups were increased in digestive cell cytosol and lighter lysosomes but not in heavier lysosomes or digestive gland microsomes following 5 days exposure to Cu. Dihydrohodamine-123 was converted to fluorescent rhodamine-123 following sequestration into digestive cell lysosomes. The results suggest a link between the lysosomal sequestration of copper, a concomitant increase in the production of oxyradicals and the potential for intracellular oxidative damage, as well as an increased capacity for oxidative damage in older animals.     

The Sweet Home rhodochrosite specimen mine,Alma District,Central Colorado: the porphyry molybdenum–fluorine connection

Intermediate sulfidation veins containing quartz–sphalerite–tetrahedrite–rhodochrosite–fluorite in the Sweet Home Mine, Alma District, Colorado were originally mined for silver starting in 1873. For the last 13 years up until 2004, however, the mine has produced world-class rhodochrosite specimens. Some of these specimens are considered to be among the finest mineral specimens ever produced and the finest of their species with values well over $1 million US dollars. The extraction, preparation, and marketing techniques pioneered at the Sweet Home operation have revolutionized the minerals specimen industry. The Sweet Home deposit is interpreted as a single pulse variant of a Climax-type hydrothermal system. Evidence for this includes (1) an age of mineralization (25.8 ± 0.3 Ma) that coincides with the age of the end stages of mineralization of the Climax molybdenum deposit approximately 7.5 km to the northeast; (2) a geochemical (Mn, W, F) and mineralogical (topaz, fluorite, hubnerite, greisen muscovite) signature typically associated with Climax-type systems; (3) the presence of porphyry rhyolite dikes, a breccia dike, and local quartz–molybdenite veins in the nearby area; (4) a small pegmatite within the mine with an age (25.9 ± 0.3 Ma) coincident with mineralization, which also contains minor amounts of disseminated molybdenite; and (5) the presence of similar-appearing gemmy, red rhodochrosites at Climax and other high-silica rhyolite systems. A significant difference is that unlike Climax-type systems, the Sweet Home hydrothermal system appears to have consisted of a single, relatively small pulse of magmatic fluid that slowly cooled and diluted with groundwater. This is inferred to have occurred at moderate depths in the order of 1.5–2.5 km below the surface. The fluids that formed the Sweet Home veins were dilute (salinity in the order of 2–4 wt% NaCl equivalent), high-temperature (temperatures of homogenization up to 370°C), and initially of magmatic origin. Gem quality ruby-red rhodochrosite at the Sweet Home Mine is nearly pure manganese carbonate with minimal solid solution with Fe⁺², Ca, or Mg. It formed at higher temperatures and salinities in comparison to lower value, pink rhodochrosite. Gemmy, ruby-red rhodochrosite is distinctly associated with highly evolved silica-rich igneous/hydrothermal systems. The high fluorine content typical of such systems suggests that Mn was transported in solution as fluoride complexes, which, in turn, favored rhodochrosite deposition at above-average temperatures and with minimal cation contamination. An erratum to this article can be found at     

Subject Index

Subject Index

Subject Index     

Experimental phase and melting relations of metapelite in the upper mantle: implications for the petrogenesis of intraplate magmas

We performed a series of piston-cylinder experiments on a synthetic pelite starting material over a pressure and temperature range of 3.0–5.0 GPa and 1,100–1,600°C, respectively, to examine the melting behaviour and phase relations of sedimentary rocks at upper mantle conditions. The anhydrous pelite solidus is between 1,150 and 1,200°C at 3.0 GPa and close to 1,250°C at 5.0 GPa, whereas the liquidus is likely to be at 1,600°C or higher at all investigated pressures, giving a large melting interval of over 400°C. The subsolidus paragenesis consists of quartz/coesite, feldspar, garnet, kyanite, rutile, ±clinopyroxene ±apatite. Feldspar, rutile and apatite are rapidly melted out above the solidus, whereas garnet and kyanite are stable to high melt fractions (>70%). Clinopyroxene stability increases with increasing pressure, and quartz/coesite is the sole liquidus phase at all pressures. Feldspars are relatively Na-rich [K/(K + Na) = 0.4–0.5] at 3.0 GPa, but are nearly pure K-feldspar at 5.0 GPa. Clinopyroxenes are jadeite and Ca-eskolaite rich, with jadeite contents increasing with pressure. All supersolidus experiments produced alkaline dacitic melts with relatively constant SiO₂ and Al₂O₃ contents. At 3.0 GPa, initial melting is controlled almost exclusively by feldspar and quartz, giving melts with K₂O/Na₂O ~1. At 4.0 and 5.0 GPa, low-fraction melting is controlled by jadeite-rich clinopyroxene and K-rich feldspar, which leads to compatible behaviour of Na and melts with K₂O/Na₂O ≫ 1. Our results indicate that sedimentary protoliths entrained in upwelling heterogeneous mantle domains may undergo melting at greater depths than mafic lithologies to produce ultrapotassic dacitic melts. Such melts are expected to react with and metasomatise the surrounding peridotite, which may subsequently undergo melting at shallower levels to produce compositionally distinct magma types. This scenario may account for many of the distinctive geochemical characteristics of EM-type ocean island magma suites. Moreover, unmelted or partially melted sedimentary rocks in the mantle may contribute to some seismic discontinuities that have been observed beneath intraplate and island-arc volcanic regions.     

Asymmetric spreading across the Juan de Fuca and Gorda rises as obtained from a detailed magnetic survey

Marine magnetic data collected along 37 east-west lines off the Washington and Oregon coast between 42°N and 48°N westward to 144°W have been interpreted. Our interpretation shows that decrease in the offsets of the anomalies across the Surveyor fracture zone and the presence of undisturbed north-south lineations east of it are the result of continuous asymmetric spreading along the Surveyor fracture zone. The survey has delineated the extension of the Blanco fracture zone north-westward to about 133°W. Movement along the Blanco fracture zone was initiated about 15 million years ago. The undisturbed north-south trend of the magnetic anomalies between latitudes 42°N and 48°N and longitudes 133°W and 136°W is interpreted as the interval (22 to 15 my) during which the Juan de Fuca and Gorda rises were one continuous structure. West of 137° the Surveyor, Sedna and three minor fracture zones are mapped.