Hybridization of the probability perturbation method with gradient information

Geostatistically based history-matching methods make it possible to devise history-matching strategies that will honor geologic knowledge about the reservoir. However, the performance of these methods is known to be impeded by slow convergence rates resulting from the stochastic nature of the algorithm. It is the purpose of this paper to introduce a method that integrates qualitative gradient information into the probability perturbation method to improve convergence. The potential of the proposed method is demonstrated on a synthetic history-matching example. The results indicate that inclusion of qualitative gradient information improves the performance of the probability perturbation method.     

A coupling of mixed and continuous Galerkin finite element methods for poroelasticity II: the discrete-in-time case

In this paper, we formulate a finite element procedure for approximating the coupled fluid and mechanics in Biot’s consolidation model of poroelasticity. Here, we approximate the pressure by a mixed finite element method and the displacements by a Galerkin method. Theoretical convergence error estimates are derived in a discrete-in-time setting. Of particular interest is the case when the lowest-order Raviart–Thomas approximating space or cell-centered finite differences are used in the mixed formulation and continuous piecewise linear approximations are used for displacements. This approach appears to be the one most frequently applied to existing reservoir engineering simulators.     

Estimating groundwater evapotranspiration rates using diurnal water-table fluctuations in a semi-arid riparian zone

In semi-arid climates, phreatophytes draw on shallow aquifers, and groundwater evapotranspiration (ET_G) is a principal component of groundwater budgets. Diurnal water table fluctuations, which often are a product of ET_G, were monitored in the riparian zone of Red Canyon Creek, Wyoming, USA. These fluctuations were higher in a riparian wetland (2–36 mm) than a grass-covered meadow (1–6 mm). The onset and cessation of water-table fluctuations correspond to daily temperatures relative to freezing. Spatial differences were due to vegetation type and specific yield, while temporal changes were due to vegetation dormancy. Ratios of ET_G to potential evapotranspiration (PET), K _c,GW, were similar to ratios of actual evapotranspiration (ET) to PET, K _c, in semi-arid rangelands. Before vegetation senescence, K _c,GW increased between precipitation events, suggesting phreatophytes pull more water from the saturated zone as soil moisture decreases. In contrast, K _c decreases with soil moisture following precipitation events as ET becomes increasingly water-limited. Error in ET_G is primarily from estimates of specific yield (S _y), which is difficult to quantify in heterogeneous sediments. ET_G values may be more reliable because the range of acceptable S _y is smaller than K _c and S _y does not change with vegetation type or soil moisture.     

Recent plumbing system of the Krakatau volcano revealed by teleseismic earthquake distribution

Spatial and temporal analysis of global seismological data 1964–2005 reveals a distinct teleseismic earthquake activity producing a columnar-like formation in the continental wedge between the Krakatau volcano at the surface and the subducting slab of the Indo-Australian plate. These earthquakes occur continuously in time, are in the body-wave (m _b) magnitude range 4.5–5.3 and in the depth range 1–100 km. The Krakatau earthquake cluster is vertical and elongated in the azimuth N30°E, suggesting existence of a deep-rooted fault zone cutting the Sunda Strait in the SSW-NNE direction. Possible continuation of the fault zone in the SW direction was activated by an intensive 2002/2003 aftershock sequence, elongated in the azimuth of N55°E. Beneath the Krakatau earthquake cluster, an aseismic gap exists in the Wadati-Benioff zone of the subducting plate at the depths 100–120 km. We interpret this aseismic gap as a consequence of partial melting inhibiting stress concentration necessary to generate stronger earthquakes, whereas the numerous earthquakes observed in the overlying lithospheric wedge beneath the volcano probably reflect magma ascent in the recent plumbing system of the Krakatau volcano. Focal depth of the deepest events (~100 km) of the Krakatau cluster constrains the location of the primary magma generation to greater depths. The ascending magmatic fluids stress fault segments within the Sunda Strait fault zone and change their friction parameters inducing the observed tectonic earthquakes beneath Krakatau.     

Determining the optimum cell size of digital elevation model for hydrologic application

Scale is one of the most important but unsolved issues in various scientific disciplines that deal with spatial data. The arbitrary choice of grid cell size for contour interpolated digital elevation models (DEM) is one of the major sources of uncertainty in the hydrologic modelling process. In this paper, an attempt was made to identify methods for determining an optimum cell size for a contour interpolated DEM in prior to hydrologic modelling. Twenty-meter interval contour lines were used to generate DEMs of five different resolutions, viz., 30, 45, 60, 75, and 90 m using TOPOGRID algorithm. The obtained DEMs were explored for their intrinsic quality using four different methods, i.e., sink analysis, fractal dimension of derived stream network, entropy measurement and semivariogram modelling. These methods were applied to determine the level artifacts (interpolation error) in DEM surface as well as derived stream network, spatial information content and spatial variability respectively. The results indicated that a 90 m cell size is sufficient to capture the terrain variability for subsequent hydrologic modelling in the study area. The significance of this research work is that it provides methods which DEM users can apply to select an appropriate DEM cell size in prior to detailed hydrologic modelling.     

Changes in nutrient uptake of phytoplankton under the interaction between sunlight and phosphate in the Changjiang （Yangtze） River Estuary

We conducted ship-board incubation experiments to investigate changes in nutrient uptake of phytoplankton under different phosphate concentrations and irradiances in the Changjiang River Estuary and its adjacent waters in China. Under 100% natural irradiance the uptake rates of phosphate, silicate, and nitrate were accelerated at high phosphate levels （1.84 μM）, while under low irradiance （about 50% natural irradiance） their uptake rates were restrained at the high but stimulated greatly at the intermediate phosphate concentrations （1.26 μM）, as the growth of phytoplankton, changes in nitrite and ammonium uptake didn＇t follow an obvious pattern. Our results also showed that there were linear relationships between nitrate, silicate and phosphate uptake at different phosphate concentrations under low and high irradiances, and the growth period of phytoplankton was prolonged both at the high phosphate concentrations under high irradiance and at the intermediate concentrations under low irradiance, suggesting that the limitation of phytoplankton growth mainly reflected changes in its growth period, and because no such environment （low irradiance and low phosphate concentrations） actually existed in a high turbidity zone, phytoplankton blooms hardly occurred there. In the absence of irradiance, denitrification occurred readily and phytoplankton was kept decreasing, which resulted in phosphate regeneration.     

Hydrothermal pattern of frozen soil in Nam Co lake basin,the Tibetan Plateau

Hydrothermal processes and the regimes of frozen soil formed in alpine regions with glaciers and lake area are complex and important for ecological environment but have not been studied in Tibet. Based on soil temperature and moisture data from October 2005 to September 2006 collected in the Nam Co lake basin, Tibetan Plateau (TP), those questions were discussed. The mean annual air temperature was −3.4°C with 8 months below 0°C. Air and soil temperature varied between −25.3~13.1°C and −10.3~8.8°C, respectively. Soil moisture variations in the active layer were small with the minimum value of 1.4%, but were influenced greatly by snowmelt, rainfall and evaporation, varying up to 53.8%. The active layer froze later, thawed earlier and was thinner, however, the lower altitude limit of permafrost is higher than that in most areas of TP. The effects of soil moisture (unfrozen water content) on soil temperature, which were estimated through proposed models, were more significant near ground surface than the other layers. The surface soil temperature decreased with snowcover, the effect of cold snow meltwater infiltration on soil thermal conditions was negligible, however, the effect of rainfall infiltration was evident causing thermal disruptions.     

Proximal and distal styles of pegmatite-related metasomatic emerald mineralization at Ianapera,southern Madagascar

The Ianapera emerald deposit is located in the Neoproterozoic Vohibory Block of southern Madagascar. The local geology consists of intercalated migmatitic gneissic units and calcareous metasedimentary rocks, containing boudinaged metamorphosed mafic/ultramafic lenses, all intruded by pegmatite veins. These units occur near the hinge of the tightly folded Ianapera antiform, within a few kilometers of the Ampanihy shear zone. Emerald mineralization is hosted by metasomatic phlogopite veins, and bodies developed within the mafic/ultramafic rocks. Based on field and textural relationships, we distinguish proximal and distal styles of mineralization. Proximal mineralization occurs at the contact of pegmatite veins with mafic/ultramafic units; in the distal style, pegmatites are not observed. Three types of emeralds could be distinguished, mainly on the basis of color and mineral zoning. Some of these emeralds have the most Al-depleted and Cr-rich composition ever recorded. Another characteristic feature to the Ianapera deposit and, to our knowledge, yet unreported, is the association of some emeralds with scapolite in metasomatised mafic rocks. Mineral inclusions are common in most emeralds and include phlogopite, carbonates, barite, K-feldspar, quartz, pyrite, zircon, monazite, bastnaesite, phenakite, plus Fe and Cr oxides. However, feldspar and rare earth element-bearing minerals occur predominantly in proximal emeralds, which also have a more incompatible trace-element signature than distal emeralds. We propose a model related to syn- to post-tectonic magmatic-hydrothermal activity. Pegmatitic bodies intruded units of the Ianapera antiform probably during tectonic relaxation. Exsolution of fluids rich in halogens and incompatible elements from the cooling pegmatites caused hydrothermal metasomatism of Cr-bearing mafic/ultramafic rocks in direct contact with the pegmatites. Local fracturing favored fluid infiltration, permitting the formation of distal mineralization. Emerald composition was controlled by the chemistry of the host rock. The presence of carbonate mineral inclusions in the emeralds and the high F-activity indicated by elevated F-contents in newly formed minerals suggest transport of Be as a fluoride-carbonate complex. It seems likely that beryl formation was triggered by precipitation of F-rich phlogopite, which removed the complexing ligand from the fluid.